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

learn more… | top users | synonyms

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 ...
14
votes
5answers
431 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
454 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
450 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
168 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
304 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
658 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
928 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
345 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
517 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
181 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
367 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
438 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 ...
4
votes
4answers
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
714 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 ...
1
vote
2answers
427 views

Where to find PySAL toolbox for ArcGIS?

I installed PySAL but i can't find the toolbox for ArcGIS.
1
vote
0answers
136 views

Is it possible to script ISIS Mapper?

Does anyone have any experience of scripting or batch processing with ISIS Mapper? I am working on a large modelling and mapping project that will require the processing of significant numbers of 1d ...
3
votes
1answer
509 views

Putting coded domains into list to be created in Python script?

I have a large list of coded domains. They are numbers like: 01.6530.01, 05.7991.08, 05.7991.01, 05.7991.06 The code will be the exact same as the description. I've tried to make a list using: ...
3
votes
3answers
324 views

Setup a map server that runs from USB to whatever computer it is connected?

I would like to have a USB that I can bring with me to any computer. That USB should hold a map server, so that the client computer wont need internet and avoid using OSM or Google maps or any other ...
3
votes
2answers
359 views

How to generate Buildings ID numbers inside parcels

In order to decrease manual data entry, I would like to automate a task to generate numbers for all buildings inside a parcel starting with 1. For example, I have 4000 parcels and each parcel has ...
3
votes
2answers
436 views

Process every 100 records in list

An arcpy merge process is failing due to the size of the datasets and I need to just merge every 100 tiles before it goes on to the next 100 tiles. I wrote the following snippet but it is just ...
8
votes
1answer
3k views

How to make a new raster using con() script in Python?

I'm new to Python and have been trying to code a raster calculator style script. Basically, I want to convert a raster that has values 0 to -10,000 into a raster that is populated by 1s between ...
2
votes
1answer
911 views

How to add a progress bar to a python plugin?

I would like to add some type of progress bar to my Qgis python plugin. Can someone point me to an example?
3
votes
1answer
224 views

BV LISA using PySal?

I want to run a LISA with a time series spatially lagged variable using PYSAL. I see that LISA will allow the definition of y & w but not Wy(t-k). Any thoughts how I can get the lagged variable ...
9
votes
4answers
8k views

Using ogr2ogr to convert GML to shapefile in Python?

I am trying to convert a GML to an ESRI shapefile using ogr2ogr utility in a python script. I have successfully installed the GDAL/OGR package via osgeo but I am now struggling to find/understand ...
4
votes
1answer
744 views

How to use the OGR/GDAL configure options in QGIS python environment

I’m looking for a solution how to use the OGR/GDAL configure options (--config) in the qgis python environment. http://trac.osgeo.org/gdal/wiki/ConfigOptions Setting special configure options is ...
2
votes
0answers
61 views

Is there a way to distinguish between digitized and pasted features?

Is there a way to tell if an added feature in QGIS is in fact a newly digitized one or a pasted feature? I think QGIS is treating both cases in the same way, emitting the same signal: ...
5
votes
2answers
1k views

Problem with the where clause in ArcGIS when dealing with variables

I am writing a python script for ArcGIS. I have an input feature class which contains multiple pipelines which is saved as a variable InputPipeline. I would like to iterate through the pipelines and ...
3
votes
2answers
242 views

How can I create a portable script for coordinate conversions to run from USB?

I would like to be able to convert Lat long values in wgs84 to the swedish RT90, the values are in CSV files under headings X and Y. I would like to have this portable, an eg runned in a script in ...
5
votes
1answer
119 views

Is it possible show the feature attributes dialog when I add a feature in code?

I want to know if it is possible to make the feature attributes dialog in Qgis show up when I add a feature by code using: QgsVectorLayer::addFeature( tFeature ) This is available when I add ...
5
votes
2answers
5k views

GEOS and Shapely installation on Windows

When trying to install Shapely on my Windows 64bit computer, I cannot get the GEOS library to work. So far, I have run the OSGeo4W installer from which I installed GDAL (I believe the geos library ...
3
votes
2answers
917 views

Using 'gdalwarp' in Python script for Grass

I'm writing a Python script that aims to re-project various LANDSAT tiles into a common projection and then import the tiles into Grass. The script at present is shown below: #!/usr/bin/python ...
5
votes
1answer
2k views

How to update shapefile data source in multiple dataframes and mxds using lyr.replaceDataSource and Python dictionary?

After reorganizing the file structure for some of our GIS base data files I need a script to update the data sources for ~ 100 shapefiles in about 50 mxd files. Some of the mxds have multiple data ...