Python is an open source interpreted programming language used in many GIS programs.

learn more… | top users | synonyms

2
votes
1answer
59 views

How can i add script with math functions to model builder?

I am very new to model builder + python. (Sorry for my english too) I am trying to make a model. Don't know if it matters to upload it. If you want i can. I want to calculate variance from a ...
2
votes
1answer
72 views

Error reading large virtual raster with rasterio: “Cannot invert degenerate transform”

I've created a virtual raster using gdalbuildvrt, which covers a large extent (England) at a high resolution. The data is fairly sparse, focused around the coastlines with no files covering inland ...
2
votes
1answer
129 views

SetFeature creates infinite loop when updating sqlite feature using OGR

when using the SQLite driver in python - setting a feature after changing an attribute field value creates an infinite loop on the first record. the same code works fine if i substitute the sqlite ...
2
votes
1answer
185 views

How would I group values in a field based on range of number using python, and then update the attribute table?

An example of what I'm trying to accomplish is described below, however they don't describe how to use the data update cursor to update the attribute column. data - Group values based on range of ...
2
votes
1answer
147 views

How to add an interpolation raster from the python console in QGIS?

I'm trying to add an interpolation layer based on an existing vector layer by using the Python Console (initially in the menu Raster/Interpolation), but can't find anything related in the ...
2
votes
1answer
264 views

Does not show shortest path in networkx (in QGIS)

I used the following commands (in qgis console) to get the shortest path with astar algorithm. But the result shows 'NetworkXNoPath' every time. Why was that? Please help me to find a solution to this ...
2
votes
1answer
267 views

FME - Use Python Caller to Set Search Envelope in Another Workspace

I have an FME workspace which has a number of different workspaceRunner transformers. The first workspaceRunner launches a workspace which buffers a pipeline and outputs the bounding coordinates to a ...
2
votes
1answer
445 views

Union multiple shapefiles with overlapping shapes

I have several shapefiles that overlap and I would like to union their geometries to get rid of those overlaps. I know how to do this with OGR and Python for single geometries. from osgeo import ogr ...
2
votes
1answer
234 views

How do I plot a line on a map, or obtain the latitude and longitude coordinates of such a line, that is 120km north of the US-Canada border?

How do I plot a line on a map, or obtain the latitude and longitude coordinates of such a line, that is 120km north of the US-Canada border, including parts of the border that cross water? (I managed ...
2
votes
1answer
235 views

Can we access the ArcMap Image Analysis tools using Python?

Esri recently added a new Image Analysis window to ArcMap - it's found under > Windows in the main toolbar: The Image Analysis window's Processing section simplifies the experience of applying ...
2
votes
2answers
2k views

Qgis - Couldn't load plugin 'processing' - name 'QsciLexerCustom' is not defined - Manjaro

UPDATED I receive the following error when I start Qgis: Couldn't load plugin 'processing' from ['/usr/share/qgis/python', '/home/v/.qgis2/python', '/home/v/.qgis2/python/plugins', ...
2
votes
1answer
696 views

Problem creating point shapefiles programmatically using python in QGIS 2.2

I am trying to create a point shapefile in the memory buffer (for a temporary label) and then zoom to the extents on this buffer. The working (in version 2.1) test code for this is: def ...
2
votes
1answer
454 views

Hydrology tool dataset in ArcGIS 10.1 for calculating slope?

I am quite new to the hydrology toolset in ArcGIS and thought it may be useful for extracting slope data from a raster. I was wondering if anyone had experience trying to extract slope data upslope ...
2
votes
2answers
214 views

Is there a way to do a python batch script for Analyze and Compress on SDE Databases?

I created and modified python scripts based off of ERSI's example (http://resources.arcgis.com/en/help/main/10.2/index.html#//003n000000v7000000), however, what would be great is to be able to batch ...
2
votes
1answer
75 views

Labeling only the latest feature

I need to label only the latest feature in a particular feature class. I actually would like to automate it, so that when I will add any new feature to this feature class, the label will automatically ...
2
votes
1answer
954 views

How to use PyQGIS ( Python QGIS API ) to open/close Update Drawing Order checkbox in QGIS2's Layers (Legend Interface)

Thanks qgis 2, the Legend Interface (layers dockwidget) now has Update Drawing Order check box on right-click mouse context menu. I think this feature can let user temporally closes rendering, then ...
2
votes
1answer
225 views

How to use Python console under QGIS 2 with gdal module?

I am trying to move my plugin from QGIS 1.8 to QGIS2.0. However, the gdal module cannot use under QGIS2.0. I try to use Python console, typing: import gdal #it's OK gdal.GetDriverByName( ...
2
votes
1answer
255 views

Clustering polygons by sum of values

I would like to group polygons that (1) are contiguous and (2) when the combined values of field X >= Y. Specifically I am grouping zipcodes that are connected. a group is when 1 or more zipcodes ...
2
votes
2answers
182 views

How to create simple select script in ArcPy?

I am new to Python and want to create a script inside the ArcMap Python window. My current workspace is a series of rivers portrayed in the state of Ohio. I have a list of fish and the rivers (from my ...
2
votes
1answer
317 views

Problem with running Python scripts in ArcMap [closed]

I'm having a problem to get my python scripts running in ArcMap. I keep getting this message no matter what script i use. Executing: Script1 Start Time: Thu Dec 05 15:37:50 2013 Running script ...
2
votes
1answer
261 views

GeoDjango - return SRID of geometry in GeoQuerySet

I'm wondering if there is a function in GeoDjango which would allow me to determine the SRID of a geometry column in a GeoQuerySet. I know I could do it with a raw SQL query (I'm using PostGIS) but I ...
2
votes
1answer
108 views

How to determine line numbers within an LAS file?

I am working with mulipulse lidar data that collects points along a number of lines within the flight path. From what it appears there are four scan lines within the data, at different intervals but ...
2
votes
1answer
143 views

Arcpy.GetRasterProperties_management values do not match values found via ArcCatalog or ArcView

I am using Arcpy to extract the mean and standard deviation values of tiff Landsat satellite images and I get very high values. Values in the thousands rather than in the low hundreds! If I look at ...
2
votes
1answer
226 views

How to join USDA SSURGO soil table data to features in QGIS?

I am trying to make a soil map using QGIS. The only problem is that I cannot join table data to features in QGIS. In ArcGIS for Desktop you just open the folder and then join the tabular table ...
2
votes
2answers
499 views

Calling variables within ArcPy functions

I have a defined function (ex. def my_function(x,y)). def my_function(x,y): """algorithm""" my_variable = my_function(filex, filey) How can I call this variable within the code block section of ...
2
votes
1answer
196 views

Does spherical equilateral triangle T contain point P?

I am working on a Python script to: compute the vertex coordinates of a geodesic sphere/icosahedron, project the triangles onto a sphere, then find which spherical triangle contains an arbitrary ...
2
votes
1answer
448 views

How to create a symbol with graduated marker size and line width in pyQGIS?

I am using QGIS + pyQGIS + Python. Layer geometry is line. Each line has a magnitude and a direction (from 1st node to 2nd node). What exactly I want to perform is "graduated width marker symbology ...
2
votes
1answer
226 views

How to smoothly join two rectangles

I am looking for the Math and steps involved not using a tool. We have PostGIS but its an old version so we are working with the data in Python. Given three coordinates: -107.37, 52.46 -103.24, ...
2
votes
2answers
546 views

How to open OSM file in Python with GDAL

I am trying to open a *.osm file in Python with GDAL. But I doesn't work. What am I doing wrong? import ogr driver = ogr.GetDriverByName('OSM') datasource = driver.Open('map.osm') layer = ...
2
votes
1answer
597 views

How to create a QGIS action which copies features between layers?

i tried (i never tried before...) to write a simple python code like action in QGIS on my "XXX" vector layer that: open a layer "YYY" copy the selected feature on my "XXX"layer paste the selected ...
2
votes
1answer
306 views

Python ArcGIS 10.1 Search Cursor

I have a file geodatabase with a polygon feature class named "AOI". The "AOI" feature class currently holds two records. I want to loop through each record and run a geoprocessing task, for example ...
2
votes
1answer
1k views

GDAL loads in Python, but not IPython (Linux)

I am trying my hardest to utilize the work of the good folks at the GeoDa center, but I cannot seem to get all of the components to load. GDAL, in particular, has been giving me some difficulty. The ...
2
votes
2answers
237 views

Change python + ogr2ogr.exe to python + ogr

I have written several complicated python scripts that end up using ogr2ogr to load or translate data. I have been converting these scripts to standalone exe files meaning the team do not require ...
2
votes
1answer
65 views

Multiple Selections from one shapefile

I have two shapfiles, Country and Blocks. I need to run a very basic selection to find out what blocks intersect the country shpe. But, the country shape has 250+ countries and I need to do this for ...
2
votes
1answer
374 views

reclassify vector data from multiple source data fields (ArcGIS) (ModelBuilder)

I am using ArcGIS 10.1 and ModelBuilder. I am not familiar with Python. My question is very similar to ReClassify Vector Data Between Values . However that question is about reclassification using ...
2
votes
1answer
437 views

Only perform Near analysis on selected rows

I am writing a Python script to perform a near analysis on features from two different shapefiles with common attributes. This is happening in a Search Cursor. However, whenever I tried to use the ...
2
votes
1answer
106 views

Is there a quicker way of checking if a feature exists with OGR

I'm writing lot's of features to a FileGDB using OGR. The slow point in my script is checking that the feature isn't already in the Feature Class, based on the ID. Currently I'm using a function ...
2
votes
2answers
569 views

Python Code: “Compare two fields and select missing records”

The code I wrote is not right!! It's selecting more than what is actually missing. I think I'm wrong on the matched code line. I need the script to make a list of the "FLAG" field and a list of the ...
2
votes
1answer
1k views

How to use CURRENT map document in Python tool?

I am attempting to use this code in a tool that I am working on: mxd = arcpy.mapping.MapDocument("CURRENT") df = arcpy.mapping.ListDataFrames(mxd)[0] XMIN = df.extent.XMin YMIN = df.extent.YMin XMAX ...
2
votes
2answers
334 views

Tabulate Area crashes Python and ArcGIS

I have a script that runs Tabulate Area on a single feature class multiple times with different rasters. In total Tabulate Area runs 5 times. 2 of the input rasters are provided as arguments to the ...
2
votes
1answer
311 views

Arcpy Network Analyst solver error

Hello GIS Stack Exchange, I am finding it difficult to load locations to sublayers using python in ArcGIS 10.0 network analyst Vehicle routing problem. I'm also having some trouble debugging my ...
2
votes
2answers
787 views

Field calculator or update cursor?

I have a feature class in which one field (let's call it "Field_A"), as yet empty, is to be populated by the contents of another field ("Field_B") multiplied by a certain number (I'll call it the ...
2
votes
1answer
394 views

How to interpolate GPS data correctly?

I'm working on a project that involves using linear interpolation to fit lines to lat/lng points. I'm pretty new at scientific programming related stuff, and I was wondering if people could steer me ...
2
votes
1answer
321 views

Search cursor stops selecting after two iterations

I have a script which selects an extent from a grid and exports the data frame. import os, arcpy mxd = arcpy.mapping.MapDocument('CURRENT') df = arcpy.mapping.ListDataFrames(mxd,"Layers")[0] ...
2
votes
1answer
263 views

Writing python plugin: How terminate plugin from script without closing QGIS?

I'm writing a QGIS conversion plugin in python to be able to create Obstacle-files for a model from shapefiles. What I want to do now is just to close the plugin cleanly. This means, to close the ...
2
votes
1answer
152 views

Script to either rotate existing text element or insert graphic element

I have about 300 maps that contain a text element that needs to be rotated. I am pretty new to python but I have heard how great it is. This seemed like the perfect task for python but I have had no ...
2
votes
1answer
485 views

How to add one or more points to an existing polyline or polygon?

Using the following code, I make a polyline with two points. vl = QgsVectorLayer("LineString", layerName, "memory") pr = vl.dataProvider() vl.startEditing() fet = QgsFeature() ...
2
votes
1answer
395 views

View HTML field data with QT4 browers as a QGIS Run Feature Action Tool

I've successfully populated a sqlite table with an HTML-encoded-table as a field attribute, and so there is literally a Table within a single field value in sqlite, for each record. This is going to ...
2
votes
1answer
400 views

need help for script checking topology using ogr in python

My main goal is to create a new feature indicating if the line features intersect each other as well as if they are not connected to a feature. My problem is how to identify the intersecting and not ...
2
votes
1answer
313 views

OGR CreateLayer returns TypeError

I'm trying to create a shapefile layer based on the spatial projection of an already existing raster. (Ultimately, this will be used in gdal.Polygonize to get a shape from my raster data.) However, I ...