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

learn more… | top users | synonyms

2
votes
1answer
484 views

UnicodeDecodeError: While installing gdal in a new Python Virtual Environment

I was trying to create a Python Virtual Environment for some initial GIS development. I installed pip and virtualenv and then tried the following commands from the Linux-Terminal to accomplish the ...
2
votes
1answer
351 views

Pyproj: converting WGS84 to Robinson

I have a WGS84 coordinate pair (37.7749295, -122.4194155). I have an image that is a Robinson Map Projection of the world. I'm trying to write a python script to locate my WGS84 coordinate on my ...
2
votes
1answer
110 views

Problem with Python codeblock recognizing a field

'<'First Post Alert> I've been through multiple iterations / gyrations, and while there may be multiple problems here (and I'll appreciate feedback on all of them), the recurring issue is that the ...
2
votes
1answer
715 views

Is there any way to fast read pixel value using gdal in C#

I create with C# a series of lat/lon value pair in the code(using i, j loop) and I want to extract the pixel value based on either lat/lon pair or pixel location (i, j). I know there is one gdal ...
2
votes
1answer
73 views

raster overlap count

Would it be possible to simply detect and count overlaps between raster files, i.e. not polygons against rasters, using arcpy? By default I don't have the spatial analyst extension for ArcMap 10.1 but ...
2
votes
1answer
793 views

Python support with QGIS 2.4 on Ubuntu 14.04

I recently updated my Ubuntu system to version 14.04. I added the sources and installed QGIS as per the directions on the QGIS website. Everything looks good, except for I do not have all the ...
2
votes
1answer
131 views

How to find missing dependencies in a distributed ArcGIS Python toolbox?

I am developing a Python Toolbox that uses a few external libraries (ogr/gdal, matplotlib, shapely, reportlab among others) which files I have directly copied in the same folder as my toolbox .py, as ...
2
votes
1answer
210 views

How to use the MMQGIS plugin through the Python console?

I know how to make Hub Lines using the GUI, I want to know how to do the same using the Python Console, a tutorial would be great. Thank you
2
votes
2answers
241 views

How to switch to edit mode vector layers, in pyqt application based on qgis?

I do QGIS based application, I need to implement the ability to modify the geometry of vector layers in manual mode (editing nodes), how to do it?
2
votes
1answer
60 views

Zoom option error on gdal_tiler.py

I'm using the tiler-tools library from: https://code.google.com/p/tilers-tools/wiki/QuickStart And specifically this command: ./gdal_tiler.py --zoom=9-15 my_map.vrt To generate tiles. When I use ...
2
votes
2answers
741 views

GDAL_CALC works but I get a python error at the end of each process that prevents automation

Python 3.4.1, GDAL 1.11.0, OSGEO4W, Windows 7, OSGEO4W Windows Shell All 64-bit (ArcGIS is installed as well 10.2 with Python 2.7) So gdal_translate works fine. gdal_calc provides perfect ...
2
votes
2answers
219 views

How to calculate monthly maximum values from daily data in a 365 band grd?

I have a MINT2005.grd file which has 365 bands representing days of year from 2005. Each band is a 32x35 matrix representing the temperature at different latitude and longitude pair. How do I ...
2
votes
3answers
630 views

convert a shapefile into graph with nodes and edges

Is there a way to convert a shapefile into graph with nodes and edges(except networkx). I'm using qgis and I want to convert a road vector layer into a graph, so that I can use it to apply shortest ...
2
votes
1answer
307 views

QGIS Python Scripting - Syntax definition

After playing around with the Model Builder and looking at the python scripts of what I make, I noticed some syntax which I can't find the definition anywhere. Here is a snippet of the python script ...
2
votes
2answers
101 views

Issues Reading a table from Geodatabase

I have scripted a process to read a CSV file into a file geodatabase. Then, I want to read this imported file and work on it by adding new fields and calculating values. However, while using the ...
2
votes
2answers
789 views

How to create a Shapely LineString from two Points

If have two points, from which I want to create a straight LineString object: from shapely.geometry import Point, LineString A = Point(0,0) B = Point(1,1) The Shapely manual for LineString states: ...
2
votes
1answer
78 views

Finding Code in QGIS Select Features Toolset

Can anyone tell me in what directory the pre-installed Python scripts for the QGIS selection tools (i.e., Select by Point, Select by Polygon) are found? I am using on a Windows 7 installation of QGIS ...
2
votes
1answer
518 views

Determine Equal-Area Slope of Longest Water Course in Catchment

I'm still very new to python and need some direction in developing a Arcpy\Python script to determine the Equal-Area Slope of the Longest Water Course within a catchment. The Longest Water Course is a ...
2
votes
1answer
273 views

How can I extract only band 4 from a geotiff multi bands to create geotiff or AAIgrid with this band using Python?

I have a geotiff file with many bands... I will like to create a .asc file or tiff file by converting only one band from the geotiff file. I don't find the solution by using python... With ...
2
votes
2answers
118 views

ArcGIS - Extract by Attibutes in Model Builder

I have a model which at one point uses the Extract by Attributes tool to produce a new raster with just the areas containing the maximum values of the first. I would think that the SQL expression ...
2
votes
2answers
238 views

Update raster value to reflect total number of points within cell- python

I have what should be a fairly simple task, but I can't seem to figure out how to do it. I have to do it in python/arcpy/gdal/whatever because I'm trying to automate this process for a hundred or so ...
2
votes
1answer
178 views

How to Format Description when writing Help for Python script tool?

I try to format the description in script tool but I don't know how it works so i don't manage to add a line break for example. Does it works with HTML tags ? When editing my script the description ...
2
votes
2answers
236 views

Why is arcpy.StackProfile_3d -> arcpy.SaveGraph_management creating empty graph?

In a python script I'm trying to save a graph to a .png file with graph output from arcpy.StackProfile_3d but the .png file shows an empty graph. The table created from arcpy.StackProfile_3d shows ...
2
votes
1answer
389 views

Using 'construct points' on all lines in a shapefile (arcpy)

I need to use the editing tool 'Construct Points' on several shapefiles that each have ~1300 polylines. I need to do this to create points every 1km along each of the polylines, which represent the ...
2
votes
2answers
274 views

How to save python ogr2ogr shapefile to memory, and stream to web client as zip

I'm trying to accomplish something similar to the site http://ogre.adc4gis.com/ I want to be able to manipulate a geojson file, and send it to the client as a shapefile, in a single zip file. I'm ...
2
votes
1answer
1k views

import qgis.core problem with QgsFeature and QgsGeometry

I've tried several previous answers to similar questions - no luck. I run Win7. I installed QGIS 2.2.0 using the OSGeo4W 32-bit installer. Problem 1) from qgis.core import * failed. Solution 1) ...
2
votes
2answers
274 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 ...
2
votes
2answers
917 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, ...
2
votes
1answer
132 views

If Then Modelbuilder, Project error

I have a model where I want to check the projection. If it is GCS_WGS_1984, I want it to go one path, and if it isn't, I want it to go another path. I created a "Check Coordinates of Input" script ...
2
votes
1answer
326 views

Transferring data between spatial databases

I would like transfer spatial data from PostGIS to Oracle Spatial. I want to use Python's modules psycopg2 and cx_oracle. Transfer Oracle > PostGIS is OK, but in opposite version I have problem with a ...
2
votes
2answers
243 views

Can I use the LEFT (string) operator in ExecuteSQL()?

In QGIS I can use the LEFT operator on a dataset like this: LEFT( "Name2",3) ='031' What is the syntax for using this expression with ogr in python? Isn't this the right syntax? ...
2
votes
2answers
604 views

How to create equidistant points in QGIS along a polygon?

I am looking for a solution to distribute with a python script equidistant points in QGIS along a polygon. I tried to use the following code example: ...
2
votes
1answer
450 views

Automatic calculation in one attribute table based on values in another attribute table using python

I am trying to fill the fields of one gdb table based on some statistics over the attributes of another shape file. As it was not possible in ARCgis directly, I am trying to use python coe for this ...
2
votes
1answer
698 views

How to make autoIncrement number start with (001…009, 011…099)?

i have a challenge with my time in my work, i have feature class called "Parcel" storing many data, i generated an autoIncrement number by using field calculator, in field "parcel code" Expression: ...
2
votes
1answer
110 views

Irrational feature ID behavior in a PyQgis script

I've written a script to determine the coordinates of a point shapefile and insert them into its attribute table. It works perfectly when the ID sequence of the indiviual features begins with 0, yet ...
2
votes
1answer
319 views

Which data providers support loading attribute tables without geometry?

I have a table in Memory from a cx_Oracle query (a list of rows, and a list of field names and their types) in my python plugin and I want to load it as a geometry-less layer (Attribute table) in ...
2
votes
3answers
120 views

How to analyze and combine overlaping elevation maps into one?

I have 30 overlapping maps that contain tree density heights (measures between 1-100%) of a large area, and I need to get the highest point by pixel and create a new map with all these highest points ...
2
votes
2answers
45 views

Assigning ID to features with similar attributes

I have a polygon feature with three attribute columns. The values in each of the attribute columns are number coded. I want to generate an ID (fourth column) for the unique combination of the values ...
2
votes
2answers
623 views

Find index location of a lat/lon point on a raster grid in ArcPy

This question has been confusing me for a while, and I'd love y'all's help. I'm trying to find the index location on a raster (in terms of column and row#) of a specified lat/lon point. My first ...
2
votes
1answer
373 views

How to assign values to field using if/then logic involving multiple fields?

I want to run a field calculation on a string field named [territory] that first looks at a string field [majr_src] to see if there is a value other than "NONE", if so [territory] = [majr_src]. If ...
2
votes
1answer
603 views

How to refresh Catalog window from ArcMap using Python?

I am running my script from ArcMap and I have the Catalog window open on the side. I delete some extra files in the end of my script but I have to refresh the Catalog window to see them go. So I ...
2
votes
1answer
668 views

ArcGIS 10.0 get sde information out of .mxd

The overall goal here is to iterate through several directories of .mxd files and document which SQL Server instance of SDE they are pointing to, with the database name and feature class names. All ...
2
votes
2answers
74 views

Accessing features via arcpy

i try to use arcpy to accessing the features of a shapefile. I tried to use a searcheCursor, like that: cursor = arcpy.da.SearchCursor(Points, ["SHAPE@XY"]) for row in cursor: x, y = row[0] ...
2
votes
2answers
651 views

Interpolation iteration using ArcGIS- arcpy

My current project requires interpolation of state level weather data. Approximately 5 years of daily data needs to be interpolated. I still need to figure out the best interpolation method. However, ...
2
votes
1answer
231 views

Model or Code to query all feature classes in geodatabase by selecting by attribute

I have a geodatabase that has 60+ feature classes that each contain records for 60+ locations. I need to export all records from each feature class for a specific location. In each FC there is a ...
2
votes
2answers
274 views

Group and export to pdf based on an attribute field using DDP - Python Script

Using Python, I want to export a map series. Each series needs to be grouped based on an attribute field which is not being used in DDP settings (page num, page index etc...). For example: my ...
2
votes
1answer
134 views

GDAL local files VS ftp files

I can successfully get raster information when the raster file is located on my local file system: import gdal from gdalconst import * filename= "MOD11A1.A2012193.h11v11.005.2012196013534.hdf" ...
2
votes
3answers
2k views

Python script for identifying duplicate records (follow up)

I am trying to calculate a field that marks duplicates and sequentially numbers them. I am using the code below that I found here: Generating sequential numbering for duplicate values in a field ...
2
votes
1answer
878 views

How to import sextante(processing) for QGIS 1.9 (2.0)

In the current QGIS development snapshot the SEXTANTE toolbox has been merged into QGIS and is now known as "processing". This is a good step, but the previous import command doesn't work anymore. ...
2
votes
1answer
70 views

arcgisscripting in CGI script

I'm trying to write a Python script that imports arcgisscripting (stuck on 9.3.1 here) which is then run as CGI by a web app. It runs fine from the Windows command line but as soon as I drop it into ...