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

learn more… | top users | synonyms

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
480 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
229 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
397 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
620 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
557 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
800 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
481 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
248 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
389 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
460 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
396 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
242 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 ...
14
votes
5answers
429 views

Working as a Python Development Project Team with ArcGIS?

We have a development project in Python (ArcGIS 10). This project involves a mixture of toolboxes, map templates, layer files, file geodatabase templates (acting as templates being imported into a ...
10
votes
2answers
6k views

Generating GeoJSON with Python

Hi I want to programatically create a GeoJSON file using polygons from a shapefile but adding attributes from my own application. This is easily done for a shapefile: def create_data_dayer(self,...
1
vote
3answers
450 views

Can one edit xls, and convert xls to csv - with a portable python or Perl module? [closed]

I am on a network and computer where nothing can be installed. On this computer i run portable apps through a USB, however I need also to to edit xls, and save as csv. Since i cannot install anything ...
13
votes
1answer
442 views

How to use GDAL with PyPy?

I am trying to optimize my python code using pypy. While I am able to import most libraries and run succesfully, I am not able to import GDAL libraries. The same libraries when imported in Python 2.7 ...
4
votes
0answers
6k views

How to inserting script into ArcGIS ModelBuilder? [closed]

I am looking for some tips about inserting a script in the ModelBuilder... I am actually building a model in the ArcGIS ModelBuilder. At one point I would like to insert a script that crosses ...
16
votes
4answers
4k views

How to buffer in pyQGIS?

I'm looking for some examples to do simple python scripts in QGIS. How would I do a buffer analysis on a dataset? I can't seem to find much in the manual and QGIS python that quite matches the Esri ...
5
votes
1answer
163 views

Where do I append point values with pyshp?

I would like to generate a .shp with five points using shapefile available through pyshp. However, when I loop through my list I get only the last value. This makes sense as I did not append any of ...
7
votes
3answers
4k views

Extract raster by mask output does not completely overlap the mask

I have a flow accumulation grid and a shapefile containing catchment areas. Now I want to do some calculations on them so I need to extract data from cells in the raster adjacent to certain catchment ...
3
votes
1answer
303 views

pair value comparison to generate new raster using arcpy/python

I want to write a function which takes in two arguments (from two rasters) and compare those 2 values against a set of hard coded pair values inside function. Then return a value which i want to use ...
4
votes
4answers
648 views

How to compute geodesic area in GeoDjango?

I have some polygons which can be anywhere in the world. I would like to compute their areas. The GeoDjango docs for GEOSGeometry.area don't specify whether it is geometric or geodesic area, but ...
9
votes
2answers
1k views

Python, GDAL and building raster attribute tables

I have an integer raster for which I would like to build a raster attribute table using Python and GDAL. I can create a GDAL raster attribute table in Python as follows: >>> rat = gdal....
6
votes
1answer
105 views

How to get info from selected object?

I am new to the Quantum Plugin development and wonder if there is more information about the following. I have added a shape file and am (off course) able to select a point from this shape file. In ...
4
votes
3answers
910 views

Programmatically creating exports from PostGIS database into different file formats (e.g. shapefile, dxf, dgn,…)?

Is there a kind of wrapper/library written in python which would enable me to programatically create exports from postgis database into different file formats (most common being, shapefile, dxf, dgn,.....
6
votes
2answers
2k views

Describing a coordinate system in arcpy

I have a python script in ArcGIS 10.1 which takes in a few inputs from the user. One of the inputs is a coordinate system. I would like to test whether the coordinate system the user chose was ...
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 ...
2
votes
0answers
80 views

Removing FID as an option from ArcGIS Tool Validation Parameters

I am currently building a script tool which requires the user to select a sorting field from a list (built using a dependency on the featureclass). This Element is working fine, but I would like to ...
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
333 views

Cannot add layer to map in arcpy

This is related to my question ArcPy replaceDataSource, but can anyone explain to me why the following code fails to add a valid layer from a .lyr file, and I get the same problems as in that question ...
5
votes
1answer
499 views

Can you synchronize tables between ArcSDE(MSSQL) and PostgreSQL(PostGIS) regularly?

I have spatial tables in ArcSDE 10.1 (MSSQL 2008 R2) geodatabases and need to update the same spatial tables in PostgreSQL 9.1(PostGIS 2.0.1) geodatabases. Basically I use my ArcSDE environment for ...
1
vote
1answer
176 views

keeping the coordinate system of raster files in the resulting raster file after operation with numpy

I'm trying to read as array two raster files using ReadAsArray in GDAL/python, perform some map algebra operation on it using numpy and then write the resulting raster file in GTiff format using ...
2
votes
1answer
360 views

Index for data driven page created atlas?

We have a parcel atlas created from a data driven page document where each page represents a particular 16 mile square township. I have created my own index of polygons instead of ArcGIS index ...
3
votes
1answer
1k views

Calculate attributes to table based on other attributes

I am going through thousands and thousands of road lines and giving them a value based on QC with imagery. So for instance for the add_data field I give it one of these five numbers: 0 = not ...
5
votes
1answer
3k views

ArcPy replaceDataSource

I'm attempting to replace a datasource in an outline map frame using ArcPy. The rest of the script seems to run OK, but after the script has run the datasource the datasource of that layer shows an ...
5
votes
1answer
1k views

Where can I find GDAL bindings for Python 3.3?

I have been using GDAL 1.9.x on Python 2.7 and took the bindings from here. Is there a place to get GDAL bindings for Python 3.3 or is GDAL not officially supported in Python 3.3 yet?
3
votes
1answer
437 views

Problem editing raster and saving with projection and transformation intact

I am trying to make a custom tool to Arcmap using Python. It is supposed to read in a raster DEM, a point shape file and a distance, and output another raster. The output raster should be ...
3
votes
3answers
6k views

How can I fix “ImportError: No module named qgis.core”?

First of all, I'm using Python 2.7.3, QGIS 1.8 (Lisboa) on Windows 7 pro 64 bits. Now let's go to the actual question: I have written some code (Python) that works well within the QGIS built in ...
5
votes
1answer
708 views

looping the mask in python for arcgis 10

I am working with raster data on a worldwide level with a resolution of 30 arc seconds. I need to export the data into a table format but given the high resolution you can not run the "Sample" (...
6
votes
2answers
1k views

How to add GDAL to existing Python installation?

I would like to download and install GDAL for Python version 2.6 [MCS v.1500 32 bit(Intel)] on win32 that supports reading FileGDB's with Python. But I am not sure how to go about this. I found one ...