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

learn more… | top users | synonyms

11
votes
1answer
2k views

Is it possible to take automatically a “screenshot” of my selected feature in QGIS using Python?

I am implementing a tool which allows the user to create a profile of a feature (in my case urban parcels). Is it possible to automatically take a screenshot of this feature (displayed in a QGIS ...
1
vote
2answers
846 views

Python module “shapely.wkb” missing

I am trying to follow a tutorial from here http://qgistips.spatialthoughts.com/2012/12/find-neighbor-polygons-in-layer-in-qgis.html#more I have downloaded the relevant python script and edited it as ...
0
votes
0answers
193 views

How to access the rules of QgsRuleBasedRendererV2 using python?

I can't get access to the rules of rule-based renderer, python says there's no attribute "rootRule", though documentation says it should be (http://www.qgis.org/api/classQgsRuleBasedRendererV2.html) ...
0
votes
1answer
689 views

Parsing list object file paths for file name

I am trying to automate a naming process that takes the file name from a list of file paths to name other folders, geodatabases, zipfiles, etc. The problem I'm running into is that, I assume, I am ...
3
votes
0answers
66 views

Programatically reorder layers in qgis? [duplicate]

Possible Duplicate: Sort layers in QGIS table of contents? I write a qgis python plugin ,and i try to connect keybord keys to reoder layers in qgis . Are there a method in qgis allows me to ...
4
votes
1answer
224 views

How to select points within polygon constrained by field?

I've created a polygon from a raster dataset: arcpy.RasterToPolygon_conversion(RasterInt, "C:/arc/poly_int.shp", "NO_SIMPLIFY", "Value") The 'value' field has numbers 1-37, which map to the ID ...
2
votes
1answer
1k views

How to make a 'Make Table Query' script in ArcGIS?

I am trying to get the make table query running in a script. IF I simply export the script from a model, it works without a problem. However, I do have to change the variable names, as this will be ...
1
vote
1answer
74 views

Styling map elements based on their properties

I am currently generating tiles based on OSM data using Mapnik and Python and am quite happy using the XML styling but I am at a loss as to how to vary the styling on an element by element basis. I ...
3
votes
1answer
297 views

Join GPS points obtained from .GPX file in Python Console

Having loaded the .gpx file into a new layer from python console, I'm unable to access the gps points in the attributes table. Essentially, i'm unable to join these points to form a polygon ...
11
votes
2answers
3k views

How do you delete an in-memory Table View?

In a Python script, I'm cycling through thousands of Excel files and using arcpy.MakeTableView to create in memory table views from the first spreadhseet in each file. How do you delete the reference ...
1
vote
1answer
279 views

Why won't QGIS let me write to a file?

I am writing a QGIS Python Plugin. I am trying to write data to a file using with open("filename.txt", "w") as f: f.write("file text.") QMessageBox.information(None, "Test", "Written to file.") ...
10
votes
1answer
454 views

Is there a way to use ArcPy and PyGRASS tools in the same Python script?

I´m trying to write a Python script, that includes different ArcPy / ArcGIS commands and want to combine it with GRASS geoprocessing tools. Unfortunately importing the GRASS libraries doesn´t work. ...
3
votes
2answers
567 views

Create Table (Data Management), Unable to Save to Feature Dataset

I'm trying to create a table in a feature dataset and get the error: Parameters are not valid. The code works when I write only to the gdb. Am I not allowed to write to a feature dataset or am I doing ...
2
votes
0answers
76 views

Arcpy: Counting number of segments where altitude rises or decreases based on DEM?

I have a DEM of Sudety mountains and linear .shp with several hiking paths. I'm trying solve one problem: How to make a script based on input DEM and input line shapes that would count number of ...
2
votes
0answers
365 views

problem in read_shp() to generate a networkx.DiGraph from point and line shapefiles

I tried to use the Python module networkx and specifically the read-shp () to ensure the transformation of a graph to a file sharpe but I always find this error noting that this is the first time I ...
3
votes
2answers
1k views

Can't Select and Export from Python Script

I am a novice in python programming but I am trying to learn by doing. I got stuck with this function a while back and can't seem to find a solution or a workaround. I am writing a script that is ...
2
votes
1answer
388 views

Saving rasters in a Python for loop fails only on last iteration

I'm a newcomer to Python scripting and I'm having what's seemingly a very simple problem but one that's stumped me enough to send me posting to these forums for the first time. I'm using a for loop ...
3
votes
2answers
1k views

Checking if a point is in a circle

I'm using shapely to check if a point is in a polygon, but I also need to check if a point is in a circle. Checking the docs, I don't see how do that... Maybe I could create a circle? circle = ...
6
votes
2answers
366 views

How can I turn a vector layer into a graph G = (V, E)?

I have a number of rail networks, in the form of shapefiles of the "line" type. I want to turn them into a more abstract graph format, i.e. the kind with vertices and edges, not chart-type graphs. ...
7
votes
1answer
196 views

Reading shapefile tolerance in python

Is there a way to get the cluster or xytolerance from a shapefile into python using free/open source software? I'm currently using the pyshp package and can't see methods for this anywhere. Also it ...
0
votes
1answer
756 views

Invalid data source issue when trying to use lyr.replaceDataSource

When I run the following code in an open mapdocument, it executes and does what I want it to, replace the datasource of the Indianacounty_shp layer with the County_Boundary feature class. import ...
5
votes
2answers
1k views

Create subsetted multiband image from multiple files using gdal or python?

Given a time series of an area (single band), how you subset them and combine the subsets into a single multiband file (eg. vrt) using gdal or gdals' bindings for python?
0
votes
1answer
2k views

Runtime error <type 'exceptions.TypeError'>: 'unicode' object is not callable

I'm working on a script to replace a layer datasource, and I'm getting an error I cannot figure out. When I run the following code in a python window inside ArcMap, I get a "Runtime error : 'unicode' ...
3
votes
1answer
815 views

Relative Hyperlinks in QGIS without end user's software dependency

I am trying to work out how can I add a relative path hyperlink to a shapefile in QGIS. I have found several great tips on various forums but only one has so far potential to be implemented. The ...
2
votes
1answer
926 views

Raster operation with grass.mapcalc

I need to iteratively add each raster map of viewshed analysis obtained from a series of single observer locations. For this i used grass.mapcalc with Python (GRASS 7 for Windows). However, my ...
4
votes
2answers
298 views

“Hello World” introduction for ubuntuGIS?

Is the GIS equivalent of "Hello world!" = "You are here?" I installed ubuntuGIS recently, but I'm a little overwhelmed navigating the network of many packages to begin accomplishing simple things. ...
5
votes
1answer
209 views

GDAL/OGR Determine if dataset is raster or vector

I am writing a script using the GDAL/OGR Python bindings and it is possible that the user may provide a raster dataset or a vector dataset as input. If the user provides a vector dataset, the script ...
7
votes
1answer
212 views

Unable to load .gpx file from Python Console

I'm trying to create a plugin which will extract the gpx files and load them onto the active layer window in QGIS. While I'm able to load it using the QGIS GUI from Layer>New Vector Layer, I'm not ...
5
votes
1answer
226 views

Add “x” in attribute table if raster is georeferenced (ArcGIS and Python)

I have a new question and I hope somebody can help me with this. I have got a shapefile which shows the extent of every raster I have. Some of these rasters are georeferenced but some not. Now I ...
1
vote
0answers
113 views

Is it possible to script ISIS Mapper?

Does anyone have any experience of scripting or batch processing with ISIS Mapper? I am working on a large modelling and mapping project that will require the processing of significant numbers of 1d ...
5
votes
1answer
97 views

Global Basin Delineation Across Dateline

I'm trying to delineate all river basins on a global 1/16 degree grid. The algorithm that the Basin tool in the Spatial Analyst Toolbox works well (fast) except it doesn't know how to handle the ...
2
votes
1answer
606 views

Why does Merge then Dissolve return ERROR 000732: Input Features: Dataset … does not exist or is not supported?

I'm not sure why I'm running into this problem. I thought it might have something to do with the merged file being locked or open, but I don't see a "LOCK" file being created in Windows Explorer. ...
6
votes
3answers
350 views

Are there any open source Python modules compatible with arcpy that will allow complex SQL queries?

I am trying to handle string (non-geometry) duplicates in a large attribute table of a shapefile without having to convert to gdb and back or connect to a postgres db. I need to be able to use a count ...
2
votes
1answer
714 views

Integrated viewshed analysis

I need to do integrated viewshed analysis of a target with several observer locations. I am going to use r.viewshed analysis function in GRASS. With this function, it is possible to perform analysis ...
3
votes
4answers
715 views

Create a Menu Containing a List of Layers That Can be Added to an MXD Using arcpy.mapping

Is it possible to create a list of layers (residing in something like a menu) that can be selected and added to an mxd using arcpy.mapping. I've been using the bit of code below (there's someone I ...
4
votes
2answers
775 views

OSM Python Utils: Modifying generate_tiles.py to create 512 x 512 tiles?

[Edited to include example image outputs] SE folks, I'm wanting to modify generate_tiles.py to create 512 x 512 tiles. At om_henners' recommendation, I tried replacing 12 of the 13 instances of 256 ...
13
votes
3answers
5k views

How to programatically check for a mouse click in QGIS

I want to know how to check for a mouse click in QGIS. I am trying to write a python plugin and want to provide functionality similar to the "Select Single Feature" tool that already exists in QGIS. ...
3
votes
1answer
896 views

How to pragmatically upload a shapefile to geoserver using cURL? I'm having problems

I'm building the curl command in a python script, but when I run the script I get this error: * Connected to localhost (127.0.0.1) port 8081 * Server auth using Basic with user 'admin' > PUT ...
2
votes
1answer
1k views

rtree python polygon index

Sorry for a very basic questions about rtree. I would like to check whether the point is inside the polygon and the distance (point to the nearest point in polygon) as well. I was thinking about ...
4
votes
1answer
469 views

Is it possible to run an arcpy execute SQL method without connecting to any database or server?

I want to simply process a layer by writing a tool script in a toolbox without having to connect any database, so I think I may have to set up a workspace environment in ArcMap, but I'm not sure if I ...
3
votes
1answer
890 views

GeoDjango Changing My Coordinates

I'm running into a problem when saving a coordinate into one of my models in (Geo)Django. My point field is setup as: location = PointField(blank=True, null = True) So while running some tests in ...
2
votes
1answer
319 views

How to programatically replace ArcGIS Online image with new imagery in mxd's?

I have many mxds that previously used ArcGIS Online Bing imagery. I want to replace those with newly purchased imagery. I know I need to use the arcpy.mapping.MapDocument module to do this though not ...
3
votes
1answer
174 views

Show feature classes as a checklist in a workspace in ModelBuilder

In my model, I would like to show all the feature classes inside a file gdb as a "checklist", similar to how the tool dialog window for Make Query Table looks: So instead of the input being a ...
4
votes
1answer
1k views

How to write values from arcpy statistical output to table?

The attached script is my stalled attempt to perform the following workflow: Create a table and create fields Loop through a folder of point shapefiles and perform average nearest neighbor analysis ...
1
vote
1answer
550 views

ImportError: cannot import name QgsVectorFieldSymbolLayer

I am using QGIS 1.7.4., downloaded April 2012, but receive an import error loading the class QgsVectorFieldSymbolLayer. As i understand this class has been added in the Python bindings on 2011-11-15. ...
0
votes
1answer
185 views

How to code Data Driven Process using ArcPy?

I need to come up with a Python code to automate a map creation process but I am a little bit of lack of python knowledge. Any help would be much appreciated. My problem is that I have a 20 mxd file ...
2
votes
2answers
320 views

Is code completion available in QGIS Python Console?

I wanted to know if the autofill functionality is available in the python console for QGIS as the one available in ArcGIS Desktop's Python Console or am I supposed to hardcode it.
2
votes
1answer
82 views

remember QGIS instance

it would be absolutely great if someone could help me with the following problem: I have a simple Python-Plugin that starts another executable in a subprocess. From within this application (which ...
1
vote
1answer
156 views

arcpy.sa.CostDistance from Python gives error about number of args but executes fine through Toolbox

I am attemping to call arcpy.sa.CostDistance to calculate the cost distance to a source point as well as generate the optional Cost Back link raster: arcpy.sa.CostDistance("input", cost, "costDist", ...
11
votes
2answers
4k views

How to split shapefile per feature in Python using GDAL?

is it possible to split a shapefile per feature in python? (best would be a solution where i can temporarily save the resulting vector objects to the memory instead to disk). The reason: I want to ...