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

learn more… | top users | synonyms

11
votes
1answer
1k views

Is there a way to access QGIS plugins in Python?

Is there a way to access QGIS plugins in Python? I'd like to write a script to do some geoprocessing, and I'd like to use the clip, regular points and convex hull tools in fTools, as well as the Road ...
2
votes
0answers
199 views

qgis python app how to add edit feature button?

I'd like to develop a QGIS Python app for plotting points and editing the attributes. I've started with the tutorials here: http://blog.qgis.org/?q=node/59 and read the cookbook ...
4
votes
1answer
961 views

Using Gui tools and arcpy.mapping

For a project I want to create a little program that does simple things like export a map from the python command line, but also users are able to edit said map (without having to open arcmap) e.g ...
7
votes
1answer
2k views

Python Script for getting elevation difference between two points

I have several stream segments 1000 Km long. I need to find the elevation difference between two consecutive points of distance 1 Km starting from the upstream to downstream. How can I get the ...
0
votes
0answers
858 views

Problem reading attribute with Python using osgeo.ogr

I am working thru "Python Geospatial Programming" and I have a problem in a program in Chapter 6 that I cannot seem to solve. The small program: import osgeo.ogr shapefile = ...
1
vote
2answers
330 views

access attributes of one class for use in another class in python add-in

I want to get combobox attribute 'dem' and use it for further processing in class ExtentData. But I am not sure how to get attributes from one class to another for processing in pythonaddin. In ...
3
votes
2answers
956 views

Csv to Shp conversion? [duplicate]

I am working on a really big script right now where I have a csv file that I have removed rows and columns from, and edited the headers. I need to create one big shapefile for the entire csv file then ...
2
votes
1answer
2k views

Reverse geocoding to get city, state, country data separately?

I have longitude and latitude data. I need to reverse geocode to get city, state, country data separately. I have checked google service. Google provides complete address to a micro level. ...
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
337 views

How to add MS SQL table using python console?

I have MS SQL Server 2008 R2 database. How to add MS SQL Table using python console in QGIS?
0
votes
1answer
111 views

Error: Exception Raised on the “for row in rows” command line?

I am comparing fields in layers. I have the AADT layer and AADTAnno layer and I want it to search the fields FLAG and TFLAG. For Example: If there is a 29C9 in AADT layer in FLAG fields then there ...
2
votes
1answer
508 views

How to “Merge Selected Features” with python?

I have been looking for the answer to this question for some time now, but can't seem to find it! My problem is straight forward: I have one layer It contains lots of features (only polygons), say ...
3
votes
1answer
212 views

Get the total distance from a point to the end of a linestring

I have a list of GPS coordinates that make up a bus route: busRoute = [(x1, y1), (x2, y2), ... (xn, yn)] I also have a point that may not be exactly on that route (due to GPS error): ...
4
votes
2answers
2k views

Is there a way to get at C's OGRLayer.Intersection() or OGRLayer.Clip() in Python?

I am working with OGR in Python and noticed that the C libraries have useful Intersection() and Clip() functions as part of the Layer class. Is there any way to get at these functions in Python? I ...
3
votes
0answers
415 views

Create easy Annotations from point features with python in ArcGIS 10.0

Because in ArcGIS 10.0 there is still no simple toolbox function to convert labels to Annotation, I try to go a substitude way with the callable Toolbox function [TiledLabelsToAnnotation][1] on a ...
4
votes
1answer
474 views

Switching between background 64 bit processing and 32 (so I can use 3rd party ie XTools)

I often run scripts now with 64 bit background processing enabled but when I want to use certain tools like Xtools table to text I have to turn it off as these tools are 32 bit. Is there a way to ...
4
votes
4answers
527 views

How to display a “Save As” dialog for raster files?

I'm writing a QGIS plugin in python. I need to show a "Save as" dialog window in order to allow the user to select the path of the output raster file. Do I need to create a QFileDialog and manually ...
4
votes
1answer
1k views

Populating a Value List in one Parameter, from Values Derived from a Second Parameter, Python

Using the Validator in the ArcGIS 10.1 tool GUI, can I loop through and List feature classes from a workspace type parameter (ie a file geodatabase) and then go on and list the values of a field in a ...
1
vote
0answers
231 views

Averaging satellite data with missing values (using Python)

I am using sea surface temperature, chlorophyll, and other data parameters as covariants in a model, and I would like to average them over a particular time period. I am currently doing this using a ...
6
votes
3answers
4k views

Raster: How to get elevation at lat/long using python?

I was wondering if anyone has some experience in getting elevation data from a raster without using ArcGIS, but rather get the information as a python list or dict? I get my XY data as a list of ...
6
votes
1answer
498 views

ModelBuilder: How to extract data by date, export to separate tables and interpolate tables to rasters

I've been using ArcGIS for a couple years, but when it comes to model builder I've only worked with very simple models. I'm hoping someone that is much smarter than myself can help me in building a ...
2
votes
2answers
121 views

How do I create a network from connected links

I have a bundle of links in 'some' networks. I need to find out which ones are connected to each other i.e. I need to find out which links are not connected to the 'main' network. I have been working ...
0
votes
1answer
204 views

Qgis - taudem plugin - python script not found

I installed the sextante toolbox on qgis. my os is fedora. grass works fine but i've got problems with taudem. it seems like some python files are missing. for instance when i try to run the ...
0
votes
1answer
414 views

openlayers python proxy.cgi script is not taking parameters

I have setup the python proxy.cgi in cgi-bin folder in my server. [myserver]/cgi-bin/proxy.cgi returns the openlayers website [myserver]/cgi-bin/proxy.cgi?url=www.openlayers.org returns list index ...
3
votes
3answers
9k views

Use Python to parse string [closed]

OK, I know that this is Python 101, but I can't seem to figure it out. I have a function setup that will return a string, specifically an azimuth in DMS. I need to take the return values and pull ...
2
votes
1answer
564 views

Python code: Using Search Cursor to find missing values in attribute table between two layers in one MXD

I need help! I am writing a program that will look through the table of a station layer and then go to the annotation layer and list all the missing values. For Example: The station layer has a field ...
3
votes
1answer
122 views

python sqlite giving errors

I have python code that I am executing but it gives me following error: <sqlite3.Cursor object at 0x7f3c7666dc00> when I execute the cur.execute(statement). the code that is below the ...
8
votes
1answer
2k views

Pros and cons of using Ruby/Rails/RGeo or Python/GeoDjango for a web based GIS

We're at the investigation phase of creating a web based GIS. We're assessing options for an ORM to use in all of our backend import/export and processing routines. From my research, GeoDjango appears ...
3
votes
1answer
4k views

Python GDAL: Save array as raster with projection from other file

I have an array of data, and for each datapoint I know the latitude and longitude. I'd like to save it as a GTiff with the same projection as other rasters I have. This is what I've tried so far, but ...
1
vote
0answers
185 views

How to programmatically change an SVG symbol in the Vector Layer?

I'm trying to change the symbol used in a layer via Python. I've been able to tunnel through a vector layer to the symbol layer, but I don't know how to change the symbol it's using. I'd like to ...
7
votes
2answers
2k views

Converting KML files for use with Python library Shapely

I am trying to get a KML-file (wijken.kml) into Shapely. The KML-file is validated against the proper XML schema, so I guess the input is correct. Routes I have tried: 1) Converting to WKT or WKB ...
3
votes
0answers
126 views

traversing KML files with ogr python

i'm new to OGR processing and KML format so be gentle: i've got a series of KML files all containing multiple regions partitioning a larger area in each. i also have a series of geocoded points for ...
0
votes
2answers
1k views

Using Python to parse an XML containing GML tags

I am trying to write a python script which pulls GML tags out of an XML file and formats them into WKT for inserting into a PostGIS database. I have been successful in doing so for an XML containing a ...
1
vote
0answers
90 views

Clipping Buffers

I have buffered a line shapefile where many lines are intersecting but not broken at intersections. I am trying to select each record in the buffer file and clip all intersecting polygons in the same ...
3
votes
2answers
1k views

Delete feature selection using Python (ArcGIS 10)

I'm having some issues with a Python script for ArcGIS 10. The idea is to delete points from a single shapefile according to the values of two columns ("Tracks" and "Tracks_1"). Therefore, I want to: ...
1
vote
1answer
543 views

How do i use gdal2tiles for creating tilecache?

I want to create tiled cache from raster file, so i downloaded gdal2tiles.py. Now i don't get what (and how) i must do with that python script file.
5
votes
3answers
2k views

Interacting between Openlayers and Python

I have a openlayers website like this one that gives the user the ability to draw polygons and then creates a GeoJSON file. I want to process this GeoJSON with a python script. The python script ...
6
votes
4answers
3k views

Efficient arcpy method of points in polygon

Coming from a background of Mapinfo and Manifold, I frequently find myself frustrated with the ArcGIS SQL toolset. In both these former programs it is trivial to write a query asking, for instance, ...
2
votes
4answers
3k views

How can you auto-populate field using spatial location of another feature and its attributes?

I want to auto-populate a field based on the spatial location of another feature's attribute. Add a point on top of a building. The point automatically adds the name of the building. This seems like ...
3
votes
1answer
359 views

OGR - Weird GeomType -2147483645 on polygon shapefile

I am trying to import a shapefile (multipolygon) in an app with ogr but I got an error about invalid geom type. If I inspect the data with ogr like this: datasource = ogr.Open("c:\\temp\\data.shp") ...
1
vote
0answers
71 views

How to identify grid shift file used by osgeo in Python?

I wrote a Python script reprojecting points from a source coordinate reference system (SCRS) to a target one (TCRS) both identified by an EPSG code. I would like my python script to retrieve details ...
2
votes
2answers
322 views

Python Toolbox Raster Parameter Error

I borrowed a simple script I found online to clip a shapefile from the extent of a raster. It works perfectly from the ArcMap built in python command prompt. I modified is slightly to use in a tool. ...
4
votes
2answers
386 views

Moran's I z-score in ArcPy

I have a ArcPy script to calculate Moran's I several times (50 to be precise). However, I am not interested in the Moran's I value, what I need is the z-score returned by the tool. Even though the ...
2
votes
3answers
911 views

Accessing GEOS functions in Python

The issue: I am using Windows 7, 64-bit operating system with Python 3.2. I must continue to use Python 3.2 because of a dependency with another API. I need some functions in GEOS that are not ...
0
votes
0answers
285 views

mapnik error - fatal flex scanner internal error--end of buffer missed

I am using Mapnik 2.0.1rc0 on a Windows 7 with python 2.7 32bits inside a GeoDjango App. I often got this error: fatal flex scanner internal error--end of buffer missed that makes python to crash ...
0
votes
0answers
16 views

Sextante Modeler Error “global name os is not defined” [duplicate]

I want to generate a model whic reclassify radter values. So i added Saga reclassify grid values tool to Modeler and saved the model. When i run the model i got this error both on Ubuntu OSgeolive 6.0 ...
0
votes
1answer
375 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 ...
1
vote
0answers
518 views

How to transfer the attributes of features based on their spatial relationship?

I'd like to check if a point feature class is placed on a line feature using python in ArcGIS 10.0. I have created a set of point features and placed them directly on a group of line features at ...
2
votes
1answer
444 views

Compare two Digital Elevation Model (DEM)

I would like to know what's the best method(s) to compare two DEMs without using any GUI soft like QGIS or ArcGIS. I'd like to do it in Python or C++, with GDAL or other similar lib. Does calculating ...
4
votes
1answer
122 views

How to return text from GP service?

So I've got a website that I want to extend by letting people look up postal/zip codes for an address they enter. I've already got all the addresses (with their postal codes) for the city in an ...