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

learn more… | top users | synonyms

3
votes
1answer
147 views

A python tool like GNU's diff for shapefiles?

I often use the GNU diff tool to compare two files or even two folders for differences. I was wondering if there is a Python based tool out there that does basically the same for shapefiles? In the ...
4
votes
1answer
403 views

PostGIS table to network analysis

I use PostGIS 2.0 and Python API (GeoAlchemy) to work with it. I want to convert topological informations (points, lines, faces) from PostGIS tables to topological network (points, lines) to study it ...
0
votes
0answers
16 views

Problem with full extent after using python code [duplicate]

I delete my previous post on the same topic because it was not clear. I made a python code for ArcGIS 9.3 which : copy a shapefile every feature in the shapefile are shift for an x and an y given. ...
1
vote
1answer
493 views

Add georeferenced png layer to map automatically

Using ArcGIS Desktop v10.2 with an Advanced level license we have georeferenced our engineering drawings in PNG format. Next we have created a polygon layer for every single drawing to see the Area of ...
2
votes
1answer
403 views

PyScripter crashes at “import processing”

I set up PyScripter latest version, 2.7) to work with QGIS (2.2, OSGeo4W install) following http://chelahmy.blogspot.ca/2011/02/setting-up-pyscripter-for-quantum-gis.html. It generally works well. I'm ...
7
votes
4answers
1k views

Is there any way to get an email when ArcGIS Server services are stopped/down?

Occasionally, but all too often unfortunately, my AGS services just stop for no reason. I don't usually notice this until the users call me and ask why their maps have no data. Is there any way (I ...
1
vote
2answers
883 views

GDAL and Python - Don't print GDAL Error messages

According to the GDAL Python Documentation, GDAL functions do not raise exceptions or errors but do print errors to the console via stdout and stderr. I have some code in my script that calls ...
0
votes
1answer
88 views

How to list..only the Field names for a group of files in a folder?

I am having...for eg, 27 files in a folder (may be a shp or tab). Is it possible to list the Field_Names (ID, Type, Date, Allotment, etc..) alone of all those 27 files in a csv format or xls format?
4
votes
1answer
413 views

How to automatically check for lidar flight line overlaps?

I have a series of lidar-based raster products created via FUSION. Since I am working with low density lidar data, many of the raster images show significant striping where lidar returns overlap on ...
0
votes
1answer
109 views

Error on rendering linestring

I am using Mapnik 2.0.1rc0 on a Windows 7 with python 2.7 32bits inside a GeoDjango App using xml. The datasource come from spatialite 4.0. My tiling script can render polygon and point features ...
1
vote
1answer
949 views

How do you set up XY Event Layer Parameters on ArcGIS Script Tool? (asking for user input)

I've made the XY Event Layer work through python script, but can't get it to work correctly as a python script tool on arc, asking for user input. Does anyone have any suggestions? This is the script ...
0
votes
0answers
27 views

Creating separate Polylines from Points using ModelBuilder and Python Script Tool? [duplicate]

I have a list of coordinate points (long and lat) of 4 locations (points 0, 1, 2, 3) saved as a feature class shapefile. I want to create a polyline that will connect points 0 and 3 and create a ...
1
vote
2answers
399 views

I need an arcPy script for choosing objects from a table at random

I am trying to write an ArcPy script that will select X number of rows from a table with N number of objects in it, then output that selection to a new layer or shapefile. I need to be able to ...
0
votes
1answer
1k views

How to add several layers to the table of contents in ArcMap using Python in one shot?

I am already able to add individual layers to ArcMap table of contents using the following codes: # Adding the layer to the table of contents mxd = arcpy.mapping.MapDocument("current") df = arcpy....
5
votes
1answer
169 views

Automate loading a table into SDE

How can I automatically load a table stored on a network share into SDE when the original table is edited? Environment is ESRI 10.1 ArcSDE on SQL Server. Preference for python or existing GP tools. ...
0
votes
1answer
92 views

arcpy: how to access element in list? [closed]

I would like to know how to get the elements from the list, maybe using for loop?... here are the code that can save all the angles into myList[] myList = [] for num in range(1,totalStop+1): ...
0
votes
2answers
105 views

Naming your Output feature class?

I am working in a work space in python on ArcGIS 10.1 on a remote desktop. Here is my code import arcpy from arcpy import env env.workspace = "E:\Example" PST = arcpy.ListFeatureClasses ("", "Point")...
1
vote
1answer
539 views

QgsGeometryAnalyzer.dissolve problems

I am completely new to the QGIS python api and I am running into an error which I cannot solve when trying to use the dissolve class. I used the basics of the PyQGIS Cookbook. In my understanding ...
23
votes
4answers
1k views

Recommend topics to be included in a Computer Science for Geospatial Technologies course

I will be instructing a course at the local university titled Computer Science for Geospatial Technologies. This is an introductory course meant to introduce computer science concepts to geospatial ...
1
vote
1answer
937 views

Export map to geotiff?

I try to export my map (layout or data view doesn't really matter) to a georeferenced tiff-file (geotiff). However, this presumably easy task is proving rather difficult to accomplish. I first tried ...
0
votes
1answer
236 views

Error when calculate field using VBscript or Python (slice a field)

I have an expression in VBscript that works well for labeling: Function FindLabel (field) if (INT(RIGHT(field, 2)) <> 1) then FindLabel = INT(MID(field, 8, 3)) &"-"& INT(RIGHT(field,...
3
votes
1answer
1k views

Script to “extract by mask” a list of raster files

I'm trying to mask a group of raster files and save the new rasters in a different folder using the same name of the original files. It doesn't sound as a difficult task, but I'm having trouble ...
0
votes
0answers
141 views

need help creating layer with geom type wkbGeometryCollection with gdal python

I am having trouble creating a layer with geom type as geometry collection. Thus I am unable to output to a shapefile. I have attached code below. dstLayer is null. This does not happen if geom type ...
0
votes
1answer
61 views

Changing band for greyscale in plugin

I am trying to make a qgis plugin to change which band in a multiband data set (netCDF where the band represents time steps) is used for a grey scale rendering. The file works fine with manual changes,...
1
vote
1answer
363 views

reproject vector layer with pyqgis with new CRS does not worl

In QGIS I am running a script in th python toolbar to export all visible vector layers to a GEOJSON with this snippet: canvas = qgis.utils.iface.mapCanvas() allLayers = canvas.layers() for i in ...
0
votes
1answer
135 views

arcpy to perform “select by attribute” and display the result on the map in a table [duplicate]

I am thinking to use arcpy to select by attribute, after that I want to display the selected features attribute in a table on the map, how can I make it? Thanks in advance.
0
votes
0answers
147 views

Import GRD and count days with temperature below zero, using python

Minimum temperature of a year is given in GRD file. I have to read that file and then calculate the frost day (less than zero degrees celcius) from it in python.
1
vote
2answers
270 views

Keep leading zeros when converting .txt into dbf using CopyRows management

I am writing a python code in order to bring my.txt file into a dbf file in arcgis. I am using CopyRows_management which seemed to be working just fine until I realized that one of the 40 fields from ...
2
votes
0answers
21 views

How to access annotation feature classes in ArcPy? [duplicate]

I want to batch replace paths in a series of .mxd files using Python. Many of the layers have moved and in some cases the name has changed or its dataset is different. I've written this code which ...
0
votes
0answers
192 views

QGIS custom form and autocompleted data entering with “validation check”

In Qt-Designer I created a custom form where I can choose serveral species names in two different combo boxes. This combo boxes are linked via "value relation" to the accordant non spatial tables (...
3
votes
1answer
1k views

After ArcGIS 10.2.1 update, Arcpy 64 bit import Error

I've updated my ArcGIS desktop version from 10.2 to 10.2.1 and have noticed that when I run the 64 bit version of IDLE I now get an error when trying to import arcpy: Python 2.7.3 (default, Apr 10 ...
1
vote
1answer
564 views

Secure credentials for scheduled python script targeting ArcSDE 10.1 on AWS

I have a python script that will be used as a scheduled task. Our dba team is wanting the credentials to be as secure as possible. I would say that using the os security system as far as locking ...
3
votes
2answers
425 views

Tabulate Area crashes Python and ArcGIS

I have a script that runs Tabulate Area on a single feature class multiple times with different rasters. In total Tabulate Area runs 5 times. 2 of the input rasters are provided as arguments to the ...
0
votes
2answers
319 views

Have I created a valid Spatialite database?

I didn't manage to find the better way to populate the exisiting Spatialite database with new computed data than just creating new columns and filling them in the Python script. You can see that the ...
2
votes
1answer
213 views

How to use the List field function for a toolbox script?

i´m trying to write a simple script for a toolbox, similar to the Addfield_datamanagement toolbox. The script adds a new field to a table of a feature class or a shapefile with attribute tables. So ...
0
votes
0answers
82 views

Error 14: Could not open CRS database

I am using Python Console and loading a Vector Layer using python commands. When the layer is added, QGIS first asks for the CRS to use, but I don't have all the CRS listed out to select the Reference ...
-1
votes
2answers
647 views

How to use raw_input() with while-loop? [closed]

my task is to execute multiple buffers with the list-loop! import arcpy arcpy.env.overwriteOutput = 1 arcpy.env.workspace="C:\\salzburg.gdb" sbgRivers="sbg_rivers" buff_name=raw_input("Please ...
1
vote
0answers
273 views

Problem with polygon overlay in raster stats using gdal

I am using the raster_stats module in python to calculate the mean carbon value for different sets of polygons. Here is the code that I am using: africa_sat = 'E:/Research_Work/data/processing_data/...
1
vote
1answer
384 views

My class using ogr to calculate distance is returning incorrect results away from the equator

In an answer to this question ustroetz kindly provided some sample code for use of ogr to do lat/lon distance and point in polygon calculations. I've created a Python class based on his example to do ...
1
vote
1answer
144 views

Excel/Google Docs > Plss Map [closed]

Is there a way to map out PLSS (Public Land Survey System) TWP (township),RNG (range) and Section based off a field in google docs? I was thinking if you had your first column with the PLSS area in ...
0
votes
1answer
77 views

Calculate text field based on codeblock query of numeric field

I am trying to populate a field with a L/R bank designation based on a route measure (as part of a larger script), but I am struggling. I have a shapefile created (tir_pts_lyr), I have a blank column (...
3
votes
1answer
661 views

Are the functions in OGRPolygon accessible through the python bindings?

I'm currently trying to use the OGR python bindings to process some shapefiles. I want to obtain the number of interior rings in a polygon feature. Unfortunately I'm unable to use the documented ...
2
votes
1answer
138 views

What is causing sr.lock on tables after looping?

The following code puts an sr.lock on each table in my List. What’s causing my tables to lock up? I understand my script is reading or opening my tables temporarily to perform the “...
3
votes
1answer
804 views

Can ArcGIS Desktop 10 Add-in Use IronPython?

Coming from the same background as the original poster in this thread, I wonder if anyone has had any success creating an Arc10 Add-In using IronPython. I have been able to access ArcObjects via ...
0
votes
0answers
259 views

MODIS images set CRS and reproject

I am developing a plugin which reads MODIS images, and I know that MODIS images are in sinusoidal projection and I need to set that CRS and project the image to WGS84 (epsg 4326) in python. I search ...
2
votes
4answers
2k views

Offline Python or C++ library for plotting locations on map?

I'd like to plot locations using longitude and latitude coordinates using Python or C++, is there a good library for such a task? I know that matplotlib is an option, I'm wondering if there are any ...
4
votes
1answer
406 views

How to read BMP file location of image used for picture symbol in LYR file?

I have a number of *.lyr files which are pointing to various .bmp files which are used as marker symbols in ArcGIS for Desktop 10.1 I was wondering if there is a way to extract the file paths of the *...
9
votes
2answers
2k views

How to access Esri ArcPy modules from Enthought Canopy?

I have ArcGIS 10.1 running on my windows system and I have also installed the new 32-bit Canopy Express Enthought Python Distribution. Is it possible to load and use the Esri ArcPy modules from ...
1
vote
0answers
229 views

How can a script tool share (dynamic) variables among different modules?

I'm running a python script from the toolbox window (as a script tool). I have a second python (tkinter) module that I use as an input interface. Now, I know that the script tool can call the ...
0
votes
1answer
122 views

USGS Earthquake GPS Metadata: Buffer Zones & Python

I'm making a python web crawler that gets USGS earthquake GPS data, and I'm wondering how I should process the latitude longitude pairs to see if it falls within a specified radius of a specified GPS ...