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

learn more… | top users | synonyms

2
votes
1answer
109 views

Subset multiband tiff file using Python in ArcGIS for Desktop?

I have a multiband (1140) raster file in tiff format that I want to subset. Each band represents a prediction of monthly rainfall (12 bands for each year) over a 95 year period (12*95=1140). I need to ...
2
votes
1answer
66 views

Cut polygon shapefile with 2 parallel lines to derive new inside polygon?

Using Python scripting: How can I clip a polygon shapefile with parallel lines (polyline shapefile) to derive new inside polygon?
2
votes
1answer
75 views

Compiling custom plugins for QGIS on Windows

I've made my first plugin using plugin builder and I want to compile it. I read that I need two commands to do it under windows: pyuic4 and pyrcc4. But when I use any of these commands I get this: ...
2
votes
1answer
83 views

ERROR 001243 preventing automatic publishing of geoprocessing service

I am trying to automate the process of publishing a geoprocessing service to ArcGIS for Server 10.2.2, using a python script: import arcpy result = r"C:\data\ContourLineFunction.rlt" ...
2
votes
2answers
73 views

csv to shp missing file

I transform a csv file into a shape file using the following script import shapefile as shp import csv import tkinter.filedialog import urllib.request import os a = ...
2
votes
1answer
81 views

Exporting layout elements and constants to Python script?

Is it possible to export the elements and constants used in a layout to a Python script? Something similar to creating a model and exporting it to Python but instead of being a model, it's a layout? ...
2
votes
1answer
61 views

Adding new Field and writing values for all features

I need to create new Field in Attribute table of vector layer, and afterwards to calculate values and store them in that Field for all features I have in that layer. I have solution which works, but ...
2
votes
1answer
196 views

how to convert kml to wkt using geodjango or python?

in geodjango i found that i can use WKT to insert geometry such as: pnt = GEOSGeometry('POINT(5 23)') but it doesn't show how to do it if the data is KML, as my data is in kml format, so my plan is ...
2
votes
1answer
62 views

How to check shape length before running code in Python/ArcPy?

I am using a code to convert a polygon into evenly spaced parallel lines for point in range(0,numPoints+1): shortPoint1 = currentLines[firstLong].positionAlongLine(lineSpaceNum*point) ...
2
votes
2answers
91 views

Problem with accessing the layer legend using pyqgis in QGIS 2.4

I'm trying to clear the layer selection using pyqgis. It appears to me that: lgd = iface.mainWindow().findChild(QtGui.QTreeWidget, "theMapLegend") Doesn't work in QGIS 2.4 to access the layer ...
2
votes
1answer
95 views

How to detetmine county of GPS coordinate?

I have a SHP file that contains all counties in the USA. I also have roughly 100 million GPS locations. I need to tally the number of GPS locations per county. A linear search through the counties ...
2
votes
1answer
70 views

Python and gdalwarp: change resolution of raster

I am new to GDAL and have been having some hard time figuring out how to change a raster's resolution and how to align it to another existing raster. I have come to the point where I have a raster ...
2
votes
1answer
95 views

Basin numbering according to Pfafstetter Method?

I would like to divide my basin into subbasins and numbers according to the Pfafstetter method. Is there an easy way to implement this in Python using ArcPy? Or are there some freely available ...
2
votes
3answers
279 views

How to increase raster extent by filling extent with no data using ArcGIS Desktop without Spatial Analyst?

I am using feature to raster to rasterize water polygons. It only makes a raster the extent of the features in the polygons, whereas I would like it to be the extent of another existing .tif image. ...
2
votes
1answer
132 views

How to batch process near feature distance in ArcGIS

I have 2000+ point shape files. These are all GPS tracks. I want to calculate the distance from each point to the nearest feature (road, train line, restaurant etc) within a 30 meter radius. I am ...
2
votes
1answer
105 views

Python console not working at all

My python console won't work at all. I can't open it. My QGIS version is 2.4. OS is win7 64. I've reinstalled QGIS a few times but I always get the same result. Strange. Something went wrong with the ...
2
votes
1answer
67 views

Use point shape file to spatial query multi polygon shape files with python?

I have an address point layer and would like to spatially query this layer using various polygon shape files to intersect this information so that each address contains the attributes from these ...
2
votes
1answer
69 views

Filtering MSSQL layers in QGIS doesnt work

In my QGIS Plugin I am filtering layers with this Python code: FeatureIterator = layer.getFeatures( QgsFeatureRequest().setFilterExpression ("Fieldname='SomeValue'") ) for feature in FeatureIterator ...
2
votes
1answer
56 views

Convert polyline to presence absence raster in arcpy

I have a polyline that I am converting to a raster in arcpy. This is occuring within a loop and at each iteration I change the processing environment to a different area of interest. In this case it ...
2
votes
1answer
625 views

Conditional Statements in QGIS raster calculator?

What is the equivalent statement of the following expression in the QGIS raster calculator? Con(((Abs("DEM" - "FSMin") > 120) | (Abs("DEM" - "FSMax") > 120)), 1, 0) I have no problem ...
2
votes
1answer
58 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
63 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
104 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
152 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
116 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
2answers
253 views

Writing shapefile with projection defined crashes fiona

I've just started trying to use Fiona to write out shapefiles. I found the example code below, which seems to work fine from shapely.geometry import mapping, Polygon import fiona # Here's an example ...
2
votes
1answer
233 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
228 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
204 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
192 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
1answer
633 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
401 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
164 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
74 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
872 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
203 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
242 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
177 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
297 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
236 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
101 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
136 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
200 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
416 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
177 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
406 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
211 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
489 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
548 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
291 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 ...