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

learn more… | top users | synonyms

1
vote
1answer
629 views

ReportLab integration into QGIS returns error message “no module named reportlab”

I want to use ReportLab in my own plugin which I have written for QGIS. I have installed ReportLab for the python versions 2.4, 2.5 and 2.6. The package is respectively installed into the ...
2
votes
2answers
478 views

Using osgeo modules in ArcGIS10.1 and QGIS 1.8 python shells

This is related to an old question of mine - Issues with osgeo - qgis.core and arcpy modules Now in new versions of all the software QGIS Lisboa (1.8), ArcGIS 10.1 and pyScripter 2.5.3. They all use ...
5
votes
2answers
2k views

arcpy.UpdateCursor for multiple fields

I've hit a wall trying to loop through multiple fields with arcpy.UpdateCursor. I have a list of files (2) with multiple fields (10-15'ish) that have its own error check in SQL, and an error value to ...
3
votes
1answer
375 views

ArcMap 10.1 Python split inconsistent string

I converted a kml file to a feature class. All of the fields from the kml were combined into two fields in the feature class. One of them contains information delimited by "/" and the other contains ...
1
vote
1answer
809 views

Write spatial operations to PostGIS from QGIS or ArcGIS

I need to write my spatial operations to a (new) postgis layer (as the data is to large for a shp file). Can this be done in the QGIS browser or from ArcGIS without SDE? When you go (for example) --> ...
5
votes
1answer
2k views

How can i obtain the data of one selected feature in QGIS using python?

I have a problem / question regarding python and QGIS. Is it, and yes how, possible to get the data of one selected feature in QGIS using python? In my example I have a layer which shows some parcels ...
9
votes
1answer
2k views

When iterating over a vector layer with pyqgis, how do I check whether a feature is selected?

While iterating over a vector layer using the following code (summerized from the example in the pyqgis cookbook), is there a way to check whether a feature is selected? provider = ...
29
votes
2answers
2k views

Streamlining Python Code for Big Data

I have Python code which is designed to take point shapefiles through the following workflow: Merge points Integrate points, such that any points within 1 m of each other become one point Create ...
1
vote
1answer
268 views

How to store AutoCAD data in ArcSDE database?

Does anyone have an idea how i can create a database instance on my ArcSDE geodatabase to store my .dwg files? I have over 1 million of these files and would be happy to automatically extract and ...
2
votes
3answers
298 views

Simple script in ArcPy to analyze building proximity to restaurants

I am trying to make a script that will go through buildings in a shapefile row-by-row and do a count of points within .25, .5,...1 mile of it. The other points are restaurants. import arcpy ...
8
votes
1answer
1k views

How to add a toolbar via python?

Through tutorials I learned how to add a toolbutton to the plugins-toolbar via python. Now I wonder how to add a complete toolbar with toolbarbuttons via python. Can anybody give some example-code? ...
0
votes
1answer
586 views

How to create a report of all feature classes in a geodatabase?

I am not getting any response from this script and no error displayed at all I am writing this script to help me display the feature classes in my geodatabase Here is my script import arcpy ...
0
votes
2answers
989 views

trouble with UpdateCursor where clause in arcgis 10.1

I'm trying to use UpdateCursor in arcpy, but the where clause I'm using seems to be causing trouble. However, when I open up the table in ArcMap and use Select By Attributes with the exact same clause ...
3
votes
1answer
687 views

Calculate bounding box + buffer for rendering a map tile with mapnik library

Mapnik library has a Projection object (I'm using Python bindings): spherical_mercator = '+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null ...
2
votes
1answer
586 views

Field Calculate Degree Minute Second in Different Format

I have a field that is in degree minutes seconds, populated from ET Geowizards COGO Inverse tool. It populates the values as 90-12-28.15 where 90 is the degree, 12 is the minutes and 28.15 is the ...
3
votes
1answer
715 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
0answers
63 views

How do I make a path cross the most possible classes in a raster?

I want to create a function and possibly use Arcmap tools to create something similar to Cost Path. Instead of just choosing the LCP, it must cross through the most different integers. Let's say I ...
1
vote
1answer
2k views

Setting a query loop

I am looking to apply a Select By Location protocol to my data, but need to first query individual IDs in my target and source layer (the same ID in each table, one at a time). How do I loop through a ...
3
votes
4answers
6k views

Trying to write a python script to crawl a directory and list all datasources

I'm editing a scrip from ESRI to try and create a list of all datasources within MXD's. I've got the script to crawl the file tree and find all the info, but I would like it to output the data for ...
6
votes
3answers
2k views

How to smooth / “reinterpolate” a raster with GDAL?

Short: I want to change the resolution of a raster and smooth the gray ramp like shown in the images bellow. The preference is to use GDAL, PIL or Numpy. Description: I'm kriging points into rasters ...
7
votes
2answers
621 views

How to catch PostgreSQL/PostGIS errors in GDAL/OGR Python bindings?

Is it possible to catch PostgreSQL output error messages in OGR after a query has failed? For example conn.ExecuteSQL('SELECT * FROM non_existing_table;') In Postgres this returns an error ...
11
votes
4answers
4k views

Programatically split a polygon with a line

In looking at this question, it occurred to me that I am not sure how one would go about using geoprocessing tools in ArcGIS to split a polygon with a line programatically. Manually, you would use ...
3
votes
2answers
527 views

Python find string with two decimal characters

I have two columns in a personal geodatabase table that I'm trying to add up (to a third field) using a python function in field calculator, however some of the string values are not real numbers. ...
0
votes
1answer
265 views

how to integrate external image processing programs within ModelBuilder

I want to use python libraries to do image processing for very large files (over 2GB), but the current python libraries can not read 16 or 32 bit image correctly. At least, I tried Python Image ...
5
votes
1answer
154 views

How to automatically place parcels on suitability map?

I need to make a map that randomly places parcels (groups of pixels or polygons) on a suitability map by first filling the most suitable areas, then the next most suitable and so on, until a ...
9
votes
1answer
553 views

Is there a way to properly have GDAL raise exceptions in Python?

I am working with GDAL in Python and have wanted to get GDAL to raise exceptions so that I can do something useful with them. For instance if a file cannot be opened with gdal.Open() I would like to ...
3
votes
4answers
2k views

Arcpy: Calculate new field from user input and existing data

Can't seem to figure out what is wrong with this create fields from a calculation. The code creates the new field in the shapefile as selected by the user but I keep getting a 9999 error on every ...
13
votes
3answers
5k views

How to find and merge duplicate points?

I have two point shapefiles: Red and Blue. Each point shapefile has a z value which represents a radius (i.e. represented as buffers in the example). A duplicate point is defined as a point that ...
0
votes
0answers
170 views

How to build QGIS with Python on Redhat?

I'm attempting to run QGIS on Redhat 5 64 bit architecture. Version of Python is 2.4.3. Version of QGIS is 1.7.4 (wroclaw). When I run qgis, the splash screen appears, followed by a QMessageBox that ...
7
votes
1answer
1k views

How to perform a table look-up using Python in ArcGIS 10.0?

Using ArcGIS 10.1 the arcpy module has a TableToNumPyArray command, which converts a table into an array. Unfortunately I don't have access to 10.1 on my work PC. Is there a workaround which will ...
1
vote
1answer
830 views

Is getting wrong attribute values from a feature in qgis

I make two polylines features(in my example both consists of two points) with the add feature tool. Then I (with code) add some attribute values on the first feature, then changes the geometry of the ...
1
vote
1answer
129 views

Selection on Feature layer (geoprocessing)

I have a feature layers with states and counties fields. I want to make a selection tool. Tool dialog should have input for feature layer and two drop downboxes for states and counties. When user ...
0
votes
1answer
358 views

MapProxy installation

I'm tryign to install MapProxy on a windows manchine. My hope is that I can install it and test it on my own computer before moving it to a server. At the moment I can seem to find mapproxy-util and ...
3
votes
3answers
3k views

“Dataset xxx does not exist or is not supported” using Table to Domain in Python

I have a directory of hundreds of *.csv files, containing code/description pairs for geodatabase domains. I'm trying to automate the creation of these domains using Python. If I manually run the ...
2
votes
1answer
936 views

Importing QGIS core into Python on OSX

I am interested in learning how to write standalone scripts using QGIS. An answer to my previous question directed me to the PyGIS CookBook which has the following instructions for importing the ...
6
votes
1answer
4k views

How to write standalone Python scripts using QGIS?

I've been learning how to use the Python Console in QGIS using references from qgisworkshop.org. I'm familiar with writing standalone scripts in ArcGIS and want to learn how to do the same with QGIS. ...
5
votes
2answers
972 views

Python - Check for number of shapefiles in a directory

I want to test: to see whether or not ONE AND ONLY ONE shapefile (which will be based on a filename wildcard) exists in a directory. Add these shapefiles to a list Use this list variable as an ...
1
vote
0answers
163 views

produce images with shapefiles and OpenStreetMap

I have a python file http://pastebin.com/jzwE6c6P which iterates over the features of a shapefile and renders a simple image by zooming to the bounding box of each feature. What I need is ...
1
vote
1answer
150 views

ArcServer Geoprocessing Service

I originally posted this in the ArcGIS forums, but I haven't been getting any bites. I've developed a set of ArcToolbox Python scripts, which I've been using on certain projects for clients. The ...
0
votes
1answer
258 views

ArcGIS Python Error Raster Domain

I´m always getting an Error with my Python-Code. Even though it is almost the same as on the ArcGIS Help Website. The Error occurs in line 19: Traceback (most recent call last): File ...
3
votes
2answers
2k views

Python in ArcGIS: How to use a variable as part of the output path

I´m trying to use a variable as part of the output path of a shutil command. Unfortunately I couldn´t work it out how to do it. Can anybody help? import arcpy import shutil Name = ...
0
votes
1answer
1k views

how to execute a python tool in ArcGIS 9.3

I have a problem executing the code below. The file is not found in split_row_0717.py Please kindly help any modification and thank you. The error said that Traceback (most recent call last): ...
2
votes
0answers
150 views

Why would an ArcGIS 10 Python custom tool fail without throwing a “failed to execute” or “completed script” message at the end?

I have written a custom tool in Python for ArcGIS that converts in batch SHP files to GML files. It processes one file at a time in a folder, file, or ZIP file and then moves on to the next. If the ...
1
vote
2answers
255 views

Python syntax error…cannot see the error!

I'm very new to python and have created a script for enabling a user to list file types in a workspace making use of the parameters function in ArcGIS so it's more user friendly. However I cannot get ...
4
votes
1answer
1k views

Pass a variable from Modelbuilder into a Python script

I'm new to modelbuilder and python. I have used modelbuilder to iterate through all rows in a shapefile (Points), process them, and produce a seperate output shapefile for each row. The filename is ...
4
votes
1answer
501 views

ogr2ogr when importing from file geodatabase and using in Python?

I want import data from my GDB to PostGIS. It works fine on some datasets (all records -data- is imported or nothing is imported) from a gdb but just doesn't import others. I can't find any ...
1
vote
0answers
527 views

Python batch text replace crashes some ArcGIS10 MXD on reopen

In an attempt to teach myself some Python, I am writing a script that will loop through a netwrok directory for all MXDs and when found, replace layout text elements with dynamic text element ...
0
votes
1answer
181 views

Qgis - Running SDA4PP [duplicate]

Possible Duplicate: How to install SDA4PP plugin for Kriging in QGIS? I'm trying to use the SDA4PP plug-in for QGIS but it doesn't work. The software says to me that I need the rpy2 plugin ...
3
votes
1answer
395 views

Python: shp file to GML using OGR

I am trying to convert a shapefile to GML 3.1.1 . I am using gdal 1.9.1 however the GML conversion it generates for a geometry feature is GML 2 and it also swaps the long/latitude coordinates. Are ...
6
votes
3answers
4k views

Select features by attribute if in Python list

I believe this will amount to largely a sql syntax question. Essentially, I am trying to complete a select by attribute in Python but based on the query of whether an attribute is present in a list. ...