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

learn more… | top users | synonyms

2
votes
1answer
123 views

Purpose of “myPlugin_dialog.py” file created by QGIS Plugin Builder

I usually use QGIS Plugin Builder to create working templates for QGIS plugins. While a generic Python plugin needs the files below; __init__.py metadata.txt resources.py (translated from ...
0
votes
2answers
345 views

Add joined data in fiona with new schema

Looked more closely at problem on how to add data to a shapefile using fiona. Docs for fiona very unclear on how to do this. I know how to add a schema after that looks very unwieldy and can not ...
6
votes
1answer
4k views

How to find and delete identical points?

I recently performed Integrate on a merged point feature class in ArcGIS 10.1. I specified "1 meter" XY tolerance for the Integrate command. The resulting point feature appears to be what I want ...
4
votes
1answer
676 views

How to use the OGR/GDAL configure options in QGIS python environment

I’m looking for a solution how to use the OGR/GDAL configure options (--config) in the qgis python environment. http://trac.osgeo.org/gdal/wiki/ConfigOptions Setting special configure options is ...
0
votes
2answers
2k views

Get points from a kml LineString

I have a KML file with a LineString. I'm trying to load the points in Python. I'd like to use the FastKml library, but am open to other solutions. FastKml has one example for reading KML, but I have ...
0
votes
1answer
213 views

How to identify a selected graphic element(s) in python

I have a situation where I am inserting items from Excel into maps a lot. When I insert the items it always makes them a fixed width and then adjusts height to maintain aspect ratio. What I would ...
0
votes
0answers
165 views

Convert raster_stats list output to array / csv file

I usually program in R, but recently I have started to use python for most of my spatial programming. My python skills are still limited, and please excuse the simplicity of the question. I am ...
1
vote
1answer
337 views

ogr2ogr/python overwrite PostGIS table

Using ogr2ogr command line I can overwrite or append to an existing PostGIS table like: ogr2ogr -f "PostgreSQL"-overwrite PG:"host=myserver port=54321 user=user dbname=mydbname password=mypassword" ...
0
votes
1answer
52 views

Is there an equivalent to arcpy.ImportToolbox() in 9.3.1?

Is there an equivalent to arcpy.ImportToolbox() using ArcGIS 9.3.1? I'm assuming it would use arcgisscripting instead of arcpy but I can't find any documentation.
7
votes
2answers
2k views

How do I use QGIS python console to iterate through simple layer queries and png outputs?

I'm really new to QGIS coding and I'm having a hard time finding (especially) the api classes/functions to carry out actions of the Query Builder window. I'm basically trying to do this: Iterate ...
5
votes
2answers
883 views

Focal statistics to calculate saddle points from DEM

I have a conundrum with Focal Statistics and kernel files. I am been trying to identify saddle points in DEMs. I am attempting to rewrite this (thanks whuber) in a python script for 10.1. I have ...
0
votes
0answers
88 views

I could use some help with a bit of code optimization

See edit below. I posted here a couple weeks ago looking for some advice on the best way to iterate through data. You can find that post here. I have since done some tweaking to the code, and ...
0
votes
1answer
136 views

automatic save outputs variables in the folder workspace using pyhon for arc map 10.1

I want to create python script with aim to import him in environment arcmap 10.1 and to make a hydrologic analysis because the scipt meybe can have many outputs raster file or feature class would ...
1
vote
2answers
191 views

Python script to check node / edge coincidence

I have written a script with arcpy to test node and edge coincidence. Each edge (line) should have a node (point) at either end. If this is true, the script does not write anything. If there is node ...
6
votes
3answers
3k views

Convert .mxd with embedded VBA command button to ArcMap 10

I'm trying to convert an .mxd file from an older version to ArcMap 10. The old version had two simple VBA scripts attached to toolbar buttons. I am trying to convert this .mxd to ArcMap 10. For ...
1
vote
0answers
305 views

PyQGIS Dialog Plugin not resetting

I have written a PostGIS search plugin that works pretty well, except two things that I just cannot figure out a solution 1) When the dialog closes and then is reopened the previous query is still in ...
0
votes
1answer
94 views

Model builder to call python script

So I'm looking to write a script that can be called by model builder. What I would like this script to do is to select three records from a feature class that have the highest field values ...
2
votes
1answer
1k views

Can anyone see what I did wrong with this Make Feature layer function? (Arcpy) [closed]

So I am just testing this function by itself before I put it in my loop. So I run the test (Python 2.7) and I get a result saying that it was create, but I do not see it in folder I assigned it to. I ...
0
votes
1answer
314 views

How to schedule a model to run at a certain time?

I'm trying to schedule a model to run at a certain time and have found this ESRI Blog post which is providing a good start. What I can't figure out though is where in the importToolbox function I tell ...
6
votes
3answers
3k views

PostGIS: parse geometry wkb with OGR

I'm trying to pull a LineString geometry out of PostGIS and parse it with OGR (python bindinds). from osgeo import ogr import psycopg2 connection = psycopg2.connect("...") cursor = ...
1
vote
0answers
198 views

How can I draw features on GeoTIFF in Python?

I have a set of data (coords + value) and want to turn it into a heatmap in GeoTIFF. Should I calculate pixels and set them as array, or there are shortcuts? This is what I do now: driver = ...
1
vote
2answers
227 views

What output colour format is expected by the QGIS data defined properties widget?

I'm trying to write a python function that takes a BGSREF attribute from a BGS geology shapefile and uses it to colour up the geology polygons of a shapefile. I posted my non-code working as a part ...
5
votes
1answer
6k views

Arcpy cursors, WHERE clauses, and date/time fields

I'm having trouble making a selection of rows that meet certain date criteria using my arcpy.UpdateCursor(). I need to update records in a field in a geodatabase table with an ascending value based on ...
2
votes
2answers
6k views

GDAL: Convert specific RGB (GeoTIFF) color to transparent

I have an RGB GeoTIFF file and I want to convert a specific RGB colour to transparent. I would like to use GDAL for doing this, but external Python scripts are also an option (though I haven't learned ...
3
votes
1answer
475 views

How to get output value from Feature Compare (schema) and append to existing csv?

I am working on a data integration project and am trying to single out feature classes which have a schema that varies from two base feature classes (one for point, one for polyline features). I've ...
1
vote
0answers
36 views

Able to circumvent QgsCredentials script called from Pre-execution

I've got a stored project that, when opened is connecting to my PostgreSQL DB, and is loading several layers. All good. I'm using a script that is called from Pre-execution (Processing/Options and ...
7
votes
3answers
4k views

How to replace non-English characters in attribute tables using Python?

I have a few shapefiles where some of the attributes contain the non-English characters ÅÄÖ. Since some queries doesn't work with these characters (specifically ChangeDetector), I tried to change them ...
0
votes
0answers
220 views

Nearest road in the vicinity of a point without queries (openstreetmap)

I know that many similar questions about map fitting have been asked. Everywhere there were answers about sql queries like this: SELECT ST_AsText (ST_ClosestPoint (pt, line)) AS cp_pt_line, ...
1
vote
1answer
93 views

Tick boxes for a qgis plugin written in python

in a python plugin for qgis I want the user to be able to select one or more of several options via tick boxes. I imagine it to look like e.g. Qgis 2.2 tool "Vector" -> "Data Management Tools" -> ...
1
vote
0answers
247 views

ArcGIS Geoprocessing Service fails on data loading

I have a Desktop Script tool that works fine when run on desktop, but does not work when the results are turned into a geoprocessing serice. Running ArcGIS 10.2 on Desktop and 10.2 on server. I have ...
1
vote
1answer
175 views

Programmatically modify attributes for vector layer in QGIS

I am trying to programmatically modify attributes for a vector layer, using calculated values currently stocked in a dictionary. Any pointers on where i'm going wrong? my dictionary has the feature ...
2
votes
0answers
158 views

How to detect broken hyperlinks with arcpy?

I have zero experience with python but am trying to search for broken hyperlinks in a large gdb (temporarily converted to a shapefile)- Used code from a previous question (Hyperlink broken path - ...
2
votes
0answers
225 views

Getting “Segmentation fault: 11” with ogr in Python

This is the code I'm trying to run: In [1]: from osgeo import ogr In [2]: datasource = ogr.Open('testdata.shp') In [3]: datasource.GetLayerCount() Out[3]: 1 In [4]: layer = ...
2
votes
2answers
385 views

How do I load mapinfo mid/mif into geopandas?

I suspect this is a trivial question to answer but how do I load mapinfo mid/mif files into geopandas? While it's trivial to load and manipulate .shp files, I can't work out what code to use: Here's ...
7
votes
1answer
2k views

Sharing a Folium Interactive Map w/ Embedded Markers

I have been experimenting with a few Python packages that have great potential for my dissertation research. Since I use pandas for pretty much everything, geopandas was a natural choice for working ...
2
votes
1answer
129 views

How to Export multiple selected features to PDFs with Python?

I have 77 polygons that mark areas of a satellite image that I want to export as separate PDFs. Each of the PDFs should contain just the polygon (and as much as the surroundings to fill out the data ...
-1
votes
1answer
470 views

NameError: name 'QImage' is not defined [closed]

Traceback (most recent call last): File "", line 1, in File "F:/QGIS/bin/test.py", line 6, in img = QImage(QSize(800,600), QImage.Format_ARGB32_Premultiplied) NameError: name 'QImage' is not ...
0
votes
1answer
129 views

Python Built-in Map and Zip function in QGIS 2.0 Processing Script--syntax error

I am attempting to run my Python script using the Processing Toolbox in QGIS 2.0. I've used the built-in Python function map() and zip() several times and am now receiving a syntax error even though ...
0
votes
0answers
140 views

GEOS assertion failed

The following script fails, if I import the 'rasterstats' library. import rasterstats import ogr geom1 = ogr.CreateGeometryFromWkt('POLYGON ((415924.806155724043492 ...
1
vote
1answer
80 views

How to refer to a variable from the conditional statement in the ExtractByAttributes function in ArcGIS 10?

Trying to use the ExtractByAttributes for rasters so that I have multiple attribute values that I want to extract. This works fine: myraster = arcpy.sa.ExtractByAttributes("source_raster", "VALUE IN ...
1
vote
1answer
264 views

Is there any python function to rotate map layers?

Is there any python function allowing a user to rotate map layers in the map canvas. If so, could you please point this out for me? I want to rotate a map layer, using a factor of rotation in ...
0
votes
1answer
527 views

How to use for loop to add and calculate new field?

I am trying to use for loop to add and calculate a newfield for 2 shapefiles and a different new field for 2 other shapefiles. When I tried to do these 2 loops as one, it did not work so I separated ...
4
votes
1answer
2k views

Python Mapping in Matplotlib Cartopy Color One Country

I have plotted a map of the world using matplotlib Cartopy. Now I want to select a specific country in the map in this case the United States and change the color. I think this is possible but not ...
1
vote
1answer
118 views

Gdalinfo and QuantumGIS does not recognize Spatial Reference even though Python does

I apply the following Spatial Reference to a file: spatialRef = osr.SpatialReference() spatialRef.ImportFromEPSG(3411) outraster.SetProjection(spatialRef.ExportToWkt() ) geotransform = ...
0
votes
1answer
80 views

How to restrict read only *.mxd file?

I have an *.mxd file and I don't want to save any edits to it after each session. Is there a way in 'arcpy' or python in general to restrict edit of *.mxd file? Any suggestions would be appreciative.
2
votes
2answers
2k views

What is the most efficient way to iterate through folders and geodatabases?

I am looking at cataloging some data, and I have code that works, it just seems to take a lot of time to run through all of the items. The way this data is structured is a bit of a hindrance, ...
1
vote
0answers
37 views

How to replace accented strings in attribute tables using Python or VBA ArcGIS 10.2 [duplicate]

I need to replace accented words, into the Field Calculator, to the same non-accented ones. For example: María ---> Maria Need some help with this VB Code: Dim data as integer Dim contador Dim ...
1
vote
1answer
723 views

Python package for processing GRIB data in windows

I would like to know what python packages can be used to process GRIB (Gridded Binary) data on a Windows system. I know there is something called pygrib, but to use that I have to use cygwin windows ...
1
vote
1answer
220 views

How to “Apply Edits” to an ArcGIS Server edit session using Python

ArcGIS Server allows the user to apply edits in an editable layer. Eg the sample layer at http://sampleserver6.arcgisonline.com/arcgis/rest/services/Military/FeatureServer/3/applyEdits will accept a ...
1
vote
0answers
15 views

ArcGIS 000210 error copy to feature class [duplicate]

I have a script tool that: Creates a feature dataset within a file geodatabase Performs a Select by Attributes on a cadastral layer Copies the selected polygons into the feature dataset. The copied ...