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

learn more… | top users | synonyms

1
vote
1answer
113 views

sa.fill creates NoData raster in standalone script ArcMap 10.0

I am trying to put together a geoprocessing workflow in Python, and I've hit a stumbling block with sa.fill. When run as a standalone script, the tool produces a raster composed entirely of NoData ...
3
votes
0answers
634 views

Batch Converting ASCII into Feature Classes using shapefile lists [closed]

Problem. I have 413 ASCII files totalling 4 billion points, and my system cannot handle converting them to MULTIPOINT so that I can generate Terrains. Therefore I would like to create 12+ subsets of ...
4
votes
1answer
624 views

How to best organize files for a large arcpy script tool?

I am writing a lengthy arcpy tool and am working with multiple python scripts, simply for organization's sake. Everything runs well when I run it from the command line, but now that I am starting to ...
5
votes
3answers
3k views

Why does arcpy.ImportToolbox() not work in ArcGIS 10 (SP5) when calling custom model for topology?

I have an arcpy/python script using ArcGIS10 and using a custom toolbox and model. The toolbox has an alias of defaultTopo and my model within the toolbox is called AddRuleToTopologyModel. This model ...
5
votes
2answers
2k views

Why does this simple Python OGR code create an empty polygon?

Please excuse my ignorance...I am new to Python OGR. I am trying to create a simple square polygon and fail...this produces an empty polygon. Where did I go wrong? from osgeo import ogr driver = ...
7
votes
1answer
637 views

How to detect a GDAL/OGR 'Warning'?

When I run a batch script using a GDAL/OGR program, I can detect if an error occured by checking ERRORLEVEL for the value 1. When using GDAL/OGR in python or .NET, I get an exception when an error ...
6
votes
1answer
3k views

Python: find a method to calculate the “inner centroid” (also known as labelPoint in arcmap) of a polygon

i am looking for a method to compute the X,Y coordinate of a inner point (also known as labelPoint in arcmap) of a (irregular) polygon. I am try to find a solution without the use of ArcGIS module ...
2
votes
1answer
646 views

arcpy.AddMessage and codepage problems

Some times in my tools I need to output Cyrillic text via arcpy.AddMessage function. Previously in 10.0 it was very easy. I had .py file of the tool saved in UTF8 codepage. This code wrote messages ...
7
votes
2answers
2k views

Does QGIS have something like the MS Office macro recorder?

I would like to use python console in QGIS to automate tasks. Is there something like a macro recorder in MS Office to get the python commands?
5
votes
2answers
2k views

How to create a shapefile consisting of 2 lines using Python OGR?

How would you create a shapefile consisting of 2 lines using Python OGR? First line is no problem, but how do you reuse MyLine and Feature objects to create the second line? Ultimately this will ...
6
votes
3answers
397 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 ...
3
votes
4answers
711 views

Is there anything faster than cursor in ArcGIS Desktop to calculate records?

What is the most efficient (fastest) way to calculate records in a layer in ArcGIS? Is there anything faster than the Cursor in ArcGIS?
1
vote
0answers
669 views

How do I find the pgAdmin connection parameters?

I need to create a connection from arcmap to pgadmin using a python module because of an arcpy script, but I am unsure of how to get the pgadmin connection parameters. I think it may be some kind of ...
2
votes
1answer
631 views

How do I select * rows of multiple columns in one query using Select by Attributes within ArcMap?

I am trying to convert a psql query string into arcpy formatting for an UpdateCursor where expression parameter: '''SELECT \"NAME_0\",\"NAME_1\",\"NAME_2\",\"NAME_3\",\"NAME_4\",\"NAME_5\",\"...
3
votes
2answers
528 views

How to call gdal_retile from Python code?

I am trying to call gdal_retile from my Python code like this: import gdal_retile and in main: gdal_retile.main("-v -r bilinear -levels 4 -ps 2048 2048 -co \"tiled=YES\" -targetDir pyramid --...
3
votes
1answer
980 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
2answers
1k views

Python: “Import Script” greyed out

I have a tool that I created from a Python script which I update from time to time. Usually I will right-click the tool and choose Import Script, but lately this has been unavailable, as shown below: ...
4
votes
1answer
1k views

How to hide/show QGisVectorLayer from Python code?

Once created a layer, how can I hide/show it? I can enable/disable rendering of a specific layer by selecting the checkbox through QGIS, but I need to do it programmatically from python code. How can ...
3
votes
2answers
323 views

Where can I get public transportation route information for USA?

I'm trying to find a list of public transit route information, mainly location of bus and train stops for the US. I've looked around at some of the Google products and openstreetmap but can't find ...
1
vote
2answers
132 views

Distribute gridded data to corresponding subbasin boundaries?

I have one ArcGIS shapefile of my river basin consisting of many subbasin. Another shapefile showing the location of gridded points over the whole basin (Its latitude and longitude are used as file ...
3
votes
3answers
855 views

How to make Field Value constant size (length) by appending zeros towards end?

How do I append my field values in my .dbf to correspond with the values of my feature's attribute table. I need to append 0's to the field value if the field value is less than 11 characters in ...
2
votes
1answer
2k 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 ...
8
votes
3answers
3k views

GDAL and Python: How to get coordinates for all cells having a specific value?

I've got an Arc/Info Binary Grid---specifically, an ArcGIS flow accumulation raster---and I'd like to identify all cells having a specific value (or in a range of values). Ultimately, I'd like a ...
4
votes
1answer
521 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 ...
1
vote
0answers
377 views

How does QGIS execute Python plugins?

I am trying to run a python script from within a python plugin. The script runs a CGIHTTPServer in a separate thread. It works as expected when executed from the command line and i can issue requests ...
2
votes
2answers
993 views

How to import qgis.core to python IDLE?

I try to follow the step in How to import qgis.core to python app? to import qgis.core with Python Shell. But when I type set PYTHONPATH=C:\Program Files (x86)\Quantum GIS Wroclaw\apps\qgis\bin it ...
2
votes
2answers
2k views

How do you create a button in ArcGIS and have it run code?

I'm brand new to ArcGIS and have only a year or two of programming experience. I know Java and a little C#. I have a Java app that processes geographical data. I need to have a button in ArcMap that, ...
24
votes
3answers
11k views

How to call gdal_translate from Python code?

Is it possible to use some gdal API to call gdal_translate from Python code? I do not mean simply executing the gdal_translate.exe from the file system, but rather call it somehow in code so I do not ...
2
votes
1answer
1k views

How to change cursor shape in pyQGIS?

I noticed that PyQGis has QgsMapToolPan() and QgsMapToolZoom() functions, and the cursor shape on the canvas will change to a hand (for pan), magnifying-glass(+/-) (for zoomIn / zoomOut). I have a ...
3
votes
1answer
489 views

Why shapefile loaded in canvas always appear in random color?

I created the simple frame .UI file, compiled with pyuic4, and tried the ShapeViewer.py to load the shapefile "world_borders.shp". However, the map is always loaded in different random colors every ...
2
votes
0answers
230 views

projection discrepancy with gdal (1.9) and proj4

I boiled the problem down to that python example: import sys try: from osgeo import gdal from osgeo import osr except: import gdal import osr if __name__ == '__main__': version_num = int(gdal....
12
votes
2answers
9k views

How to export only certain columns to a CSV file in ArcGIS?

I've written a python script using arcpy which outputs a polygon feature class into a File Geodatabase. I have added a function to export the attributes to a separate CSV file. I am using the code ...
7
votes
1answer
3k views

How to remove a vector QgsVectorLayer from QGis python plugin

I managed to add QgsVectorLayer in QGIS through a python plugin, for example: vl = self.iface.addVectorLayer(uri.uri(), layerName, self.dbConn.getProviderName()) But how can I remove it?
2
votes
1answer
405 views

Uploading shapefile to GeoServer using cURL and replace file?

In the sequence of this: How to upload a shapefile to geoserver using cURL? which was properly solved, I now need to be able to replace the existing file, cause my process is to be repeated every ...
6
votes
2answers
623 views

How to redistribute population data using Python in ArcGIS?

I am trying to redistribute population data from census centroids among to housing density.Would you please help me to do that? I do have census area centroids (points) and the boundaries belong to ...
4
votes
2answers
570 views

How to upload a shapefile to geoserver using cURL?

I need to automatically (using python) upload a shapefile into my geoserver. I'm trying to do as described in here but i'm getting this error: java.io.IOException: Wrong magic number, expected 9994, ...
6
votes
2answers
2k views

Get X, Y coordinates in the center of a raster cell

I am trying to get the X, Y coordinates of all the cells of a raster file, using Python 2.6 and ArcGIS 10.0. I would like to get the X, Y coordinates in the center of the cell. Is there any way to ...
8
votes
2answers
2k views

How to buffer a vector shapefile using ogr python?

I am trying to learn how to use ogr in python using the country and populated places datasets from http://www.naturalearthdata.com/downloads/50m-cultural-vectors/. I am attempting to use filters and ...
5
votes
6answers
823 views

How to divide a feature class into ranges?

I have a shapefile consisting of around 4000 polygons. The polygons are created from buffering around a similar number of islands somewhere in the Norwegian archipelago. Due to an error in the ...
2
votes
1answer
1k views

Create new polygon layers from postcodes latitude & longitude

First at all, I want to say that I am geocoding enthusiastic but have no formal background on GIS. I'll do my best to use the right wording and make myself clear what I am striving for. Also the only ...
3
votes
2answers
1k views

How to debug 'Layer is not valid error' from Python plugin?

I'm trying to understand what goes wrong with a QGIS python plugin. Basically I'm trying to render a series of POINT. The query works fine, at least giving it by hand from inside PostGIS DB: "SELECT ...
2
votes
1answer
488 views

NEAR_3D function without the need for 3D Analyst, Arcpy

I am looking for a python script which would do similar as the ArcGIS 3D Analyst NEAR_3D function. Basically I would like to programatically in python find nearest points of one point feature class ...
3
votes
0answers
249 views

How to merge SeaWiFS and MODIS data or use other merged products?

I hope to produce a time-series for sea surface temperature (SST) and chlorophyll-a (Chl-a) for use in MaxEnt without the use of arcgis, preferring Python tools to others. I am new to satellite data ...
4
votes
1answer
401 views

How to set page size in ArcPy?

It is possible to set page size in ArcPy? I'm looking for any tricky way to set the page size, using pure Python or ArcPy. I'm using ArcGIS 10.1 for Desktop.
8
votes
1answer
469 views

How to debug a QGIS python plugin with popups?

I'm trying to write a python plugin for QGis. I'd like to debug some DB queries and I need a way to display the result. I thought about using a popup to display it, but I don't know how to show a ...
2
votes
1answer
397 views

gdal.py as part of bat script in windows 7

I need to run the following bat file from a MS4W shell which has gdal_retile.py as the final option. I have added C:\ms4w\python\gdal to my pythonpath but within a pyscripter shell I can not import ...
4
votes
5answers
1k views

Python loop Polygon to raster not working

I want to convert a large number of different shp polygons to raster in ArcGIS 10.1, using an existing field column "class" for the new raster classes. I am not very fluent in python but I thought ...
2
votes
2answers
1k views

Problems converting CSV to Shapefile using pyshp

I have a python script which converts csv data to shapefiles. It converts but won't get angle from csv. Can anyone help me? The CSV data looks like this: fid_nr;x;y;angle 1;55.474040000;13....
2
votes
0answers
243 views

Localized Linear Histogram Match using python

I would like to use the local linear histogram match for filling the missing pixels of a Landsat 7 image with another older image (SLC Gap-Filled Products). Do you have any idea how can I do it by ...
7
votes
4answers
2k views

Arcpy's replaceDataSource Error

I've checked the other answers here and none of the answers worked for me... It's a "simple" matter of finding a layer by name in a map and changing it's datasource. Here's my script: #-- New ...