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

learn more… | top users | synonyms

1
vote
0answers
132 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
492 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
323 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
355 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
421 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 ...
12
votes
2answers
2k views

SP1 for ArcGIS 10.1 has extension for 64-bit Python for background geoprocessing — implications for other Python use

If the 64-bit Background Geoprocessing for ArcGIS Desktop 10.1sp1 is installed how does that affect Python configuration for general use? Esri has released SP1 for the ArcGIS 10.1 product, and for ...
2
votes
1answer
887 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
220 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
720 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
239 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
905 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 ...
4
votes
1answer
1k views

book/tutorial on using rpy2 in python for kriging?

I am proficient in Python, but have no knowledge of R. Can you suggest a good tutorial in using rpy2 or any other R library to perform Block/Cokriging "Kriging"
2
votes
1answer
207 views

How to construct a line from points in pyshp?

I am reading a shapefile for a pipeline and what I've noticed is that if I just draw a line from one point to the next (in the order I read them), there will be points where the line terminates, and ...
6
votes
2answers
2k views

ArcObjects + comtypes at 10.1 and newer

Has anyone played around with using ArcObjects in Python at version 10.1? I am having all sorts of problems, such as TypeErrors when wrapping the ESRI OLB modules and AttributeErrors when attempting ...
1
vote
0answers
191 views

zonal statistics in Python crashes

I have a python script that works (it's worked for several other projects in the past and it worked for my current project initially). But whenever I run it, Python crashes sporadically. It would go ...
2
votes
1answer
679 views

Convert a List to a file list to use in an arcpy tool

I just can't figure out how to feed a list of files into an arcpy tool... I have List_tif.append(os.getcwd()+"\\Stage1\\"+name+"\\GeoTif\\"+name+"_"+str(ddp.pageRow.PageName) + ".tif;") That ...
2
votes
2answers
169 views

How to use model-builder iterator with a custom bit of arbitrary code?

How do I integrate a in-line variable from modelbuilder into a python script? Specifically, I'm using an iterator, and I'm trying to include the 'name' variable in a python script. I believe the ...
6
votes
1answer
536 views

How to draw shapes with highest possible decimal accuracy

I am trying to read input of lat and long to draw the shape based on condition.The input is standard array of lat and long with 15 decimal accuracy. I am realizing - I am loosing some decimal accuracy ...
5
votes
2answers
301 views

FindLabel not Defined

I'm trying to make an advanced label expression using the python parser and I keep getting the following error: def FindLabel ( [Direction] ): return "{}° {}' {}\"".format( [Direction].split("-")...
1
vote
0answers
103 views

How to get the geometry of a polyline to a specific precision?

In python, I use point = geom.asPolyline() but it looks like I am not getting any precision just whole numbers? How can I get the geometry and a certain precision?
3
votes
1answer
258 views

How to use the QgsDiagramRendererV2?

I'm writing a plugin for QGIS in python in the QGIS API doc, I can see the QgsDiagramRendererV2 class, under the qgs.core section http://www.qgis.org/api/classQgsDiagramRendererV2.html What I'm ...
4
votes
1answer
1k views

raster2pgsql: Strange error when creating a table

Hi GIS Stackexchange people, When I try to load a asc file into my psql setup, I get a strange error message. I have created the sql input file via raster2pgsql -s 4236 -I -C -M gis_data.asc -F ...
2
votes
0answers
41 views

How can I set python path and version for QGIS? [duplicate]

Possible Duplicate: Couldn’t load SIP Module When I'm starting QGIS I'm getting the problem: "Couldn't load SIP module. Python support will be disabled." And there is also: Python version: 2....
9
votes
1answer
649 views

How to split two interscting polygons by the average linestring of the shared area (QGIS / python shapely)

I have two polygons, POLYGON((-5 -5,-5 5,5 5,5 -5,-5 -5)) POLYGON((1 1, 1 7, 7 7, 7 1, 1 1)) and wish to split them evenly. The result should be a split along the red marked line. I'm using ...
5
votes
2answers
1k views

How to improve Python script processing speed (performance)?

I've been working on a script which essentially reads in a large amount of point and polygon data and generates frequency values for a featureclass of 2.5 acre hexes covering California. I have found ...
3
votes
3answers
2k views

How to measure minimum distance between objects in two layers using arcpy?

I am working on an alien species analysis looking at the spread of a species from the mainland to islands in an archipelago. I am doing the analysis using arcpy and I have succeeded in isolating risk ...
5
votes
2answers
3k views

How to show a message box in QGIS?

I am starting to develop a plugin for Qgis using python, which i am very new at, and would like to view some results of variables by displaying a message in QGIS. How can I do this?
7
votes
3answers
539 views

How to improve positioning of images taken with GPS and Gyro info?

I'm sort of 'recreating' google street view with MIT right now. We have a car with a camera on top and we're driving at 10mph down a road. As we drive we're shooting pictures at 30fps and each photo ...
10
votes
1answer
744 views

How to see what lines have changed in an updated Shapefile?

I'm using my county's GIS data to update the roads in OpenStreetMap. I have a copy they published last year and one that was just published, and I'd like to find all the LineStrings that have either ...
2
votes
0answers
147 views

Convert geojson labels MBR to MapInfo text object

Im using Dymo (https://github.com/migurski/Dymo/) to generate automatic point labels placement. The output of Dymo is a geojson containing MBRs of the positioned labels with the names of the labels in ...
2
votes
1answer
447 views

Edit Shapefiles Attribute data in GeoServer Programatically

I'm rendering a shapefile using GeoServer and using OpenLayers library to render it in the web browser. Now I have created the option of editing the attribute data in web browser. I want the edited ...
9
votes
2answers
5k views

Python scripts that run inside ArcMap vs ones that run outside?

I'm just starting to get into Python scripting for work.. I'm currently creating a script to automate a process. Basically, it asks the user the client name, gets a projection if available, creates ...
3
votes
1answer
206 views

How to get raster layer pixels contained in a shape defined in a vector layer

I am new in QGis python plugins dev and I have a "maybe stupid" question (in the title :-) ). I try to retrieve pixels (from a raster layer) which are "contained" in a polygon defined in a vector ...
3
votes
1answer
741 views

How do I programmatically export a single Data Driven Page to JPEG?

I am trying to write a python script to export a individual pages to a jpeg. I am getting the error "PageLayoutObject: Error in SetDataDrivenPageID". Can anyone see where I am going wrong? import ...
4
votes
2answers
1k views

Programmatical Setting of 'Store Relative Pathways' for multiple MXDs ArcGIS 10.1

I am looking to move large number of mxds to a different location on the same drive. Many mxds haven't got 'Store relative pathways' option checked therefore this operation will result in broken links ...
15
votes
3answers
10k views

Writing numpy array to raster file

I'm new to GIS, and I'm lost with something I feel ought to be relatively simple. I have some code that converts infrared images of Mars into thermal inertia maps, which are then stored as 2D numpy ...
1
vote
1answer
506 views

Does arcpy have an issue with directory names starting with a number?

I have a script to export DDP extents to tif and when I change arcpy.mapping.ExportToTIFF(mxd,"tests\BA\BA10k_" + str(ddp.pageRow.PageName) + ".tif",df,9306,9306*ar,1800,True) to arcpy.mapping....
0
votes
1answer
163 views

Extending ArcGIS Desktop using Python vs Visual Basic vs Java? [closed]

I have to extend the features of ArcGIS. I've been asked to create a model written in Python that can run out of ArcGIS 10. This model will have a simple user interface where the user can drag/drop ...
4
votes
2answers
479 views

How to pass an output stored in a single variable into another script using arcpy

I am running set of processes each with different outputs in a python script. Is there a way to set a variable as global so that just the output of that single step can be accessed from other scripts?
6
votes
3answers
886 views

How to extract area programmatically?

is there a possibility to extract the area covered by one geometry of a shapefile (polygon)? I tried with thepoly.area() but get "nan" as result. Ok. I think you are right and there is ...
4
votes
1answer
1k views

Python script/Model to populate cross street names at intersections

Has anyone created a Python version/model of this VBA script from ArcScripts which populates the to and from street names at street intersections? A logical workflow is outlined quite nicely in this ...
4
votes
1answer
6k views

How to convert WMS to SHP?

I'm looking for tool that will be convert WMS to vector layer. Are they ready tool?(for example in python)
4
votes
3answers
2k views

Copy dbf data into a shapefile in Arcpy

Basically it's to perform a one-to-many join in a python script, from a dbf to a shapefile. I saw (Greg Krakow's method) using QGIS, but I can't use it (I'm working on it, but QGIS is not yet ...