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

learn more… | top users | synonyms

8
votes
2answers
2k views

What are pros and cons of Python add-ins vs .NET add-ins to ArcGIS for Desktop?

I found in ESRI .NET help: Python All ArcGIS Desktop applications include an embedded scripting language called Python. Many areas of ArcGIS—particularly geoprocessing—are accessible through ...
0
votes
1answer
874 views

Conditional Calculation of a text Field Value Based on antoher text field value ArcGIS10.1 [duplicate]

I Created a new field to explain some zoning codes in an attribute Table. I basically want the new field "Zone" to have a string value based on the field "Zoning" For example if the "Zoning" field ...
1
vote
1answer
675 views

issues with updating data source, script runs but nothing happens, arcpy

I have seen many questions posted that are similar to mine but still havent found the answer so figured I would post my script. I want to change all the data sources found in all mxds in a location to ...
1
vote
0answers
246 views

GDAL and Python: How to “pad the extent” or otherwise handle band.ReadRaster() for area falling slightly outside the image bounds?

In a nutshell, I'm trying to create a bare-metal tiling utility using GDAL and Python. As you can see in the graphic below, which shows the perimeter of a raster and the OSM-structured tile index, ...
2
votes
1answer
411 views

How to assign values to field using if/then logic involving multiple fields?

I want to run a field calculation on a string field named [territory] that first looks at a string field [majr_src] to see if there is a value other than "NONE", if so [territory] = [majr_src]. If ...
0
votes
1answer
113 views

Looping on the Select_analysis tool

Im trying to perform a loop on the Select_Analysis tool in Arcgis. Code looks like this: inputfile = open("locality.txt", 'r') # locality.txt is a txt file of locality names loc_list = [] for line in ...
4
votes
1answer
5k views

How to edit attribute tables using Python/ArcPy?

I am new to Python, so sorry for such a trivial question. There is a shapefile (let's call it Structures.shp) and the script I'm writing needs to update the value of only one record in only one field ...
1
vote
1answer
587 views

How to avoid setting env.workspace and other paths/settings all the time

I have some scripts that are constantly setting env.workspace of arcpy. Usually, I have a caller script(script1.py) that may call some other script (script2.py). I've created a file globals.py with ...
2
votes
1answer
3k views

How to use ArcPy to copy selected features to an existing feature class?

Is there a way to copy a selected feature into an existing feature class? Sort of like append would do, but only the selected features. I want to iterate through a table and select features with a ...
4
votes
1answer
497 views

GRASS 6.5 or 7.0 - Need Help with Loop Script

I am new to Python and scripting in general. This is also my first time a posting question. My overall goal is to create a script that I can cycle through all the raster maps in a mapset and then ...
1
vote
2answers
304 views

Optimising Min/Max temporal raster search in python/gdal

I need to search across ~100 1200x1200 raster tiles and produce a min and a max raster. Along the way i also want to exclude values that are over 1 or less than 0. What would be the quickest method to ...
3
votes
3answers
908 views

ArcPy nested loop problem

I have a question connected with ArcPy. I was trying to make a nested loop with 2 cursors - 1st with 2 elements and 2nd with 5 elements. The statement should return 10 connections, e.g. 1-1, 1-2, ... ...
3
votes
1answer
2k views

What expressions for QGIS field-calculator?

I do not understand what language I have to use in field-calculator and labels on QGIS: SQL or Python? And where I can find the possible script/expressions to use? in the image you can see the ...
2
votes
0answers
633 views

Large point to Polygon by buffer-join-buffer-dissolve issues

I have a few large datasets (about 10 million points) which I am buffering, then joining, dowgrading and finally buffering and dissolving. This is all in a model which works fine with sample smaller ...
7
votes
3answers
958 views

Antipodal distance (or polygon fetch or polygon diameter) for concave polygons

I am working on some class examples in Python implemented within ArcMap to calculate the antipodal distance within a polygon. This is fairly routine for convex polygons, however, for concave ...
1
vote
1answer
200 views

Raster Smoothing python [duplicate]

I want to perform raster smoothing on a DEM file. I want to to assign the elevation of each pixel using the average elevations of its 8 neighbors. The ESRI's "Python scripting for ArcGIS" has a ...
0
votes
1answer
471 views

How to set up parameters in Python toolbox tools?

I build a model builder and export it to python. How can I define parameters for python script in python properties? May you explain about Data type in script proprties? thanks alot
3
votes
1answer
375 views

QGIS Arch Linux compilation error

Just pulled QGIS from git and am trying to compile it, here's the relevant error: [ 89%] Building CXX object python/CMakeFiles/python_module_qgis_core.dir/core/sipcorepart3.cpp.o In file included ...
1
vote
0answers
329 views

Matplotlib Basemap Resolution with PyQt makes Python Quit Unexpectedly

I'm having difficulties using Matplotlibs Basemap in combination with a PyQt4 interface. Basically I create a Pyqt4 Canvas in which I plot a Matplotlib Basemap and some markers. I call the function in ...
3
votes
3answers
1k views

Problem with arcpy.CalculateField_management

I have a table with a column named "Kategorie", filled with values. I added a new column named "RULE". I want to fill this column with Strings depending on "KATEGORIE". I have a PY-Script and it ...
10
votes
2answers
319 views

Determine minimum number of polygons from shapefile to cover an area of interest

I have a large number shapefiles representing areas of interest for an analysis that will be conducted using various sources of satellite imagery (IKONOS, RapidEye, etc.). Unfortunately, the imagery ...
5
votes
1answer
3k views

If Then Else in ArcGIS model builder

I have a model which iterates through feature classes in a geodatabase. For each feature class, it runs the check Geometry tool. I then use the "Get Count" tool to count the number of rows in the ...
0
votes
2answers
94 views

What queries can I use to select data?

Which queries can I use to select the following information: How many towns are within a distance of 3km of both a road and a railway and are within 6km of the coast? And any towns that have both a ...
1
vote
0answers
244 views

I am getting a blank layer preview from geoserver after Quick Exporting my shapefile to postgis and then getting postgis to point to geoserver [closed]

What could be the potential problems? Looking at the layer's data fields/form the one thing that isn't filled is the Native SRS and its source, but I have already passed in the projection info as it ...
2
votes
0answers
490 views

Finding coordinates from a nested grid as a numpy array

So this is a pretty basic question, but I'm trying to figure it out. I've written a script that allows a user to read in a raster and provide a cell resolution, and outputs the amount of cells ...
3
votes
0answers
597 views

Move shapefile from one folder to another using python [closed]

I'm working with shapefiles and using the ArcPy library. My script is creating a new shapefile based on a join and I want to know how to move the newly joined-shapefile from the environment location ...
1
vote
0answers
161 views

Efficient computation of path profile with Numpy and GDAL

I have numpy array with SRTM elevation data and I need calculate path profiles - vectors with elevation data from one point (transmitter) to every points of my map. Here's my current solution: def ...
8
votes
2answers
230 views

Feeder bus service

First of all, a little background. I work for a regional transit agency. We are making a "diagnostic" about our feeder bus service. We'd like to know what proportion of our users could take the bus ...
1
vote
0answers
168 views

AttributeError: Error in getting tool information when using TableToTable tool

I am trying to use the arcpy.TableToTable_conversion tool in my python script, but seem to be getting the following error on the line arcpy.TableToTable_conversion("FINAL_NLPG", ...
1
vote
1answer
87 views

CalculateField_management

How to put the field from arcpy.GetParameterAsText() into the expression here: CalculateField_management(in_table, field, expression, {expression_type}, {code_block}) Now it is like ...
1
vote
1answer
310 views

Writing array to vrt

I would like to write an array from a single band image to a 3 band vrt where each band has the ColorInterpretation() tag of red, green, blue respectively. My code below; from osgeo import gdal ...
5
votes
1answer
510 views

Create MultiPolygon from QgsGeometry.fromMultiPolygon

I have a list of polygons' points like [[x1,y1],[x2,y2],[x3,y3],[x4,y4]] [[x5,y5],[x6,y6],[x7,y7].[x8,y8]] now I want to create multipolygon from it by using feature = QgsFeature() ...
2
votes
1answer
360 views

Python ArcGIS 10.1 Search Cursor

I have a file geodatabase with a polygon feature class named "AOI". The "AOI" feature class currently holds two records. I want to loop through each record and run a geoprocessing task, for example ...
0
votes
0answers
165 views

Identify segments that share a last vertex with the first vertex of selected segment

This has been marked as a duplicate to Connecting line features and determining the length of the longest line but I believe there are some subtle differences. And to be honest I'm not an advanced ...
1
vote
1answer
258 views

Grass GIS and Python Error

I am trying to use Python with Grass GIS. To get started I try to read in a raster TIFF. The import of grass.script works fine: import grass.script as g But if I am trying to run a command I get ...
1
vote
2answers
1k views

Python - GDAL : Polygon to point shape file

I want to read all the points in a polygon shape-file and use their coordinates to make a new point shape file. then find the min and max of the X coordinates. So here is the code I wrote for the ...
3
votes
2answers
2k views

Python: how to create a point shape-file from a text file

I'm writing a python code to read the points in a polygon shape-file and save them in a point shape file. I'm not sure about my approach so if you have a better idea (even if it is totally different ...
0
votes
1answer
507 views

ArcPy Mark ALL duplicated values

I'm sorry if this is too similar to other posts but I'm an arcpy newbie and I can't quite get what I need. I have accomplished this task in the field calculator (thanks to another post here, but I ...
8
votes
2answers
3k views

Split a feature when intersecting with a feature of another layer using PyQGIS/Python?

I have a buffer layer (green polygon) which I want to split to two polygons whenever it crosses a barrier (blue line). I have been trying to use "splitGeometry" method, but I just can't get it to ...
1
vote
0answers
320 views

Couldn't load PyQGIS, QGIS 2.0.1 on Ubuntu 12.04 Precise

After upgrading to QGIS 2.0.1 I still get the error: Couldn't load PyQGIS. Python support will be disabled. Traceback (most recent call last): File "", line 1, in File ...
1
vote
1answer
88 views

How to select pixel with smallest number in stack raster using Python script in ArcGIS for Desktop?

I have 6 tiff file in ArcGIS which are identical in term of geometry (same frame, same number of row and column). The pixel value of each tiff file is integer and range from 1 to 4. So if I overlay ...
1
vote
2answers
1k views

How do you get all the results from a Selection made by Attribute or Location using ArcPy?

Using arcpy I'm trying to get the results of a selection, or what is currently selected as a list. Basically I'm attempting to build an addin which on click will either perform an operation where I ...
1
vote
1answer
888 views

ArcGIS 10.1 Converting a model tool to a script tool

I have recently created a model tool for a class project and my instructor would like it to be converted to a script tool for practice purposes. I have exported the script from my model and loaded it ...
1
vote
0answers
993 views

arcpy spatial join input parameter problems [closed]

Using ArcGIS 10 SP4 with all extensions... I have a parent directory (J:\...\spatialjoin) containing two sub-directories of feature classes I want to spatially join. Directory 1 (target features) ...
0
votes
1answer
367 views

How to export multiple MXDs to PNG without using VBA?

Okay, so I can feel that I am so close to solving this issue with the help of Jason (below) but I have modified to use as a script. The GetParameterAsText returns the following error when trying to ...
7
votes
2answers
4k views

Getting XY map coordinates on mouse clicks as input variables into a python script

This is my first time coding in python, so bear with me. I wrote a script to calculate an analytical solution to find the orientation of a plane based on three points. I'm using ArcGIS 10.1. # Input ...
3
votes
0answers
238 views

python.multiprocessing and “FATAL ERROR (INFADI) MISSING DIRECTORY”

I have a question which is closely related to this old question here: python.multiprocessing and "FATAL ERROR (INFADI) MISSING DIRECTORY" I have also looked at this, which seems to relate ...
1
vote
0answers
38 views

How to find out which layer had a Selection made on it? [duplicate]

If the user makes a select using arcpy.SelectByAttribute() or arcpy.SelectByLocation() how do you find out which layer they made the select so you can use arcpy.SearchCursor(layer)?
1
vote
0answers
93 views

Obtain virtual output directory path in geoprocessing service

I am publishing a geoprocessing service and I want to be able to access the virtual output directory path configured for the service. Now I know that gp services will magically return virtual paths ...
2
votes
2answers
646 views

How to add values in X column to dictionary or list using ArcPy?

My apologies if this question has been asked before. I've spent some time searching for a similar question with no results. Is it possible to add each value under X column to a list using ArcPy? I ...