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

learn more… | top users | synonyms

2
votes
0answers
179 views

How to display images as cross sections in 3d

I want to display a base map and several vertical cross sections in 3D. The base map image is easy to rectify as a geotiff. How can I also display the cross sections? The cross sections are at ...
-2
votes
2answers
415 views

Python script not working after exporting it from ModelBuilder

trying to figure out why my script wont work (I should mention I'm still at an intro level). This is a script I exported from ArcGIS model builder which worked fine using that: # Import arcpy module ...
1
vote
2answers
166 views

Old Windows Installers for python library fiona

Due the the recent bug as presented here for fiona 1.1, i am obliged to use a version of fiona in 1.0.x on windows 8.1. however, i can't find older versions of the installer at here. Are there any ...
2
votes
0answers
70 views

Using the geocode function with cursors

How do I use a cursor with a geocodeaddress function? I want to be able to: 1) read from a table 2) only geocode addresses(from table1) whose street number is not 0 using a cursor 3) then update a ...
1
vote
1answer
3k views

Delete first two letters from a string field

I'm looking for a field calculator python statement to delete only the first two letters of string entries. For Ex: !Name!.strip([0:2]) I know this syntax is incorrect. Thanks a lot
0
votes
1answer
166 views

Installing Lecos plugin for QGis 2.0.1 on linux fedora 20

I try to install Lecos plugin, which requires python imaging library (pil), but the installation fails with the following: "you need to have the image library pil installed" I installed the ...
5
votes
2answers
1k views

how would one get the end points of a polyline?

I am attempting to write a script to get the end points of a series polylines, and then use those end points to create a polygon. I'm not sure which tools to begin to use on this project. I would ...
2
votes
2answers
541 views

How to add an image to an existing data frame using Python?

This is the first time I am using Python on my own, and I have searched for this answer for a few hours with no luck. I am wanting to take a .JPG image and put it in a data frame. The image has to ...
-2
votes
1answer
578 views

QGIS converting to ASCII question

I am attempting to do some analysis on a raster file and I need to convert it to ASCII. In QGIS, I used [Raster -> Conversion -> Translate] with the asc file format, and I suspect it is simply ...
0
votes
0answers
181 views

ArcMap 10.2 SplitLineAtPoint setSplitRule not calculating RATIO

I am trying to set a length field to use the RATION split policy. I set the split policy for the feature layer but when I use the SplitLineAtPoint_management it does not calculate the RATIO length. ...
1
vote
0answers
109 views

header information is missing from image, after writing using opencv

I am doing complement of an image using open-cv python, while writing the image, open-cv is not writing the header information in the complemented image. my image data are in TIFF format. i am doing ...
1
vote
1answer
139 views

Is it a good idea to use GRASS Algorithms for Geoprocessing in a Server Environment?

I´m designing a web-service to run in a server environment that will do some geoprocessing actions to datasets uploaded by users and send the results back. I´m programming in Python, and so far I´ve ...
1
vote
2answers
464 views

Using python to get extent and split coordinates into list

I have the following code: for fc in fcs: desc = arcpy.Describe(fc) ext = desc.extent ext.split(" ") print ext I'm getting error: AttributeError: Point: Get attribute split does not ...
3
votes
2answers
489 views

Python, GDAL, VRT: help wth code efficiency and processing time

I'm going to attempt to make this question as succinct as possible while still including relevant information, as I have trouble with being concise and am trying to avoid writing unnecessary stuff. ...
0
votes
1answer
109 views

python tool parameter: what datatype to list unique values of field?

I'm attempting to convert a 9.3 arcobjects tool to a 10.2 python tool. The form needs to allow user input of a shapefile, a field, and then the user needs to enter numbers corresponding to unique ...
2
votes
2answers
2k views

Can you change the processing extent in environment settings using python to a shapefile?

I have managed to do this manually at Geoprocessing>Environments>Processing Extent and setting the shapefile in the drop down bar. I was wonder if this could be automated with a python script? My ...
0
votes
0answers
194 views

Invalid MXD filename in a python script

I have a script that prints and exports a batch of MXDs. It runs fine as long as the lenght of the path and mxd are shorter. I can't tell what the 'sweet spot' is, but I know it runs fine on my ...
0
votes
1answer
72 views

ArcGIS Geoprocessing - how to set a condition for data

I would like to write a script in Python which takes data and does some things with it. So far, I have written the following code : import arcpy arcpy.env.workspace = "C:\wody.gdb" ...
2
votes
1answer
466 views

Why does OGR driver.Open() for File Geodatabase give AttributeError: NoneType object has no attribute Open?

There are a number of existing Questions that hint at a possible Answer to this Question, but none seem to be providing me with the solution to my problem: How to access feature classes in file ...
1
vote
2answers
464 views

How to form SQL Expressions with mix of integer/numeric and text/string fields/values in ArcPy?

I am trying to create a query for the select by attribute function within a Python script.This is my current code. lot = '96' dp = 756421 lotFile = ...
1
vote
1answer
120 views

Python script to select a point within a cell?

I have two features: points and grid cells (vector). I want to select a point within a grid cell (some cell do not have point) using python script and create a new feature layer with below conditions: ...
1
vote
2answers
153 views

Fill Sinks of multiple rasters in a Python expression

I have 31 rasters that represent different watersheds and I would like to use the Fill (Spatial Analysis) tool to fill the sinks of all rasters at once in Python. All the rasters have the same name ...
2
votes
1answer
389 views

Transferring data between spatial databases

I would like transfer spatial data from PostGIS to Oracle Spatial. I want to use Python's modules psycopg2 and cx_oracle. Transfer Oracle > PostGIS is OK, but in opposite version I have problem with a ...
2
votes
0answers
2k views

smoothing results from scipy griddata interpolatioin

I am using the griddata interpolation package in scipy, and an extrapolation function pulled from fatiando: import numpy as np import scipy from scipy.interpolate import griddata import ...
1
vote
0answers
276 views

How to utilize GPSQLExpression data type in Python Toolbox

I am trying to write my first Python Toolbox for ArcGIS (neat stuff by the way). The tool is simple enough - it connects to SQL Server and downloads data for the user. The trick here is, I want the ...
1
vote
0answers
628 views

Watershed delineation; how to program

I've been trying to roll my own watershed delineation program, taking a flow direction grid and pour points as input. I coded up a quick search algorithm in Python which inspects the flow direction ...
-2
votes
1answer
1k views

How to create a link in the message body of email with Python? [closed]

I want to send an email with python that contains as a link in the body of an email as well. This link is actually a path to a file. If a receiver clicks on link It should open a folder where the file ...
0
votes
1answer
435 views

How can i clip vector layer in qgis and python?

I want to write a python function to create a shape file layer from clipping an input layer using a mask. Does someone know the way how to do it?
0
votes
0answers
110 views

How to bulk change the schema name of layer properties in ArcMap

I would like to change the data source location (schema) in layer property in ArcMap. Our data are stocked in an Oracle Spatial database. Recently, the DBA has changed the schema name and then I ...
1
vote
1answer
165 views

Delete Turned Off (Unchecked) Layers from Arcmap with Python

I have about 100 figures that have been updated countless times and a lot of the old data was left in each .mxd but was only turned off (unchecked) and not removed. I would like to remove these ...
2
votes
1answer
525 views

ListRasters not recognizing raster files in python

I have 2 rasters in my workspace, one in TIFF Format and the other in GRID format: "C:\ExampleFolder" When I use the python script: import arcpy >>> ...
1
vote
0answers
107 views

Automate data entry into Access from ArcMap

I have an Access 2007-2010 database that I'm currently working on by manually searching for roads in a defined set of areas in ArcMap 10.1 and then filling in the area's ID in an empty column in the ...
0
votes
1answer
292 views

How to use the output file of a script as an input file for the rest of the script in arcpy?

If my question does not make sense or needs clarification please let me know. I am struggling to solve it so I was struggling to phrase it into a question. I have this code here: import arcpy,os ...
3
votes
0answers
152 views

Calculate size of area where data values exceed certain limit

I am using Python and Basemap to plot sea ice concentration within the Arctic region. I need to calculate the total area (square meters) where the sea ice concentration exceeds a certain limit (e.g. ...
1
vote
1answer
72 views

Two boundary sets intersect - some are genuine diffferences, others slivers resulting from poor digitising. How to differentiate?

Problem: I have two boundary sets. Intersecting reveals that there are genuine differences between them - as expected. However, there are also a host of slivers due to their being digitised by ...
2
votes
1answer
437 views

“ImportError: No module named gdalsrsinfo” What's the issue?

in the python console I can call import gdal, and gdal is recognized as a module; however, most gdal programs, such as gdalinfo and gdalsrsinfo are not recognized. I checked the GDAL source files to ...
2
votes
0answers
377 views

Zonal Statistics using Gdal and Python

I recently came across a loop_zonal_stats function by @ustroetz found here: Issue Trying to create Zonal Statistics using Gdal and Python.. I am using a (single) polygon shapefile and a raster file to ...
1
vote
1answer
533 views

How to zoom to selection in arcpy.mapping and only see selected features?

I selected data by using select layer by attribute, then I need to only select these selection and zoom in to the selection. my code is this, it can zoom in base on the scale I set, but I see other ...
0
votes
1answer
392 views

How to create a loop for multiple points with cost distance / cost path in ArcGIS 10.2 in python

I would like to know how to simplify the following code in order to execute it more efficiently than copy paste 82 processes (with a loop I guess?) but I don't know how to do it. There are 4 main ...
2
votes
0answers
49 views

How to Complete ModelBuilder Process Quicker?

I've built a script tool in ArcGIS's ModelBuilder to collect parameters from the user and use them as arguments for another python script. The output of this script is a text file. The only GIS ...
4
votes
1answer
908 views

How to calculate total number of data cells in raster?

I am using ArcGIS 9.3 and python 2.5. My code is shown at the bottom of this post. I have a raster containing fire intensity information at 200m cell resolution. I run the script and the first count ...
0
votes
2answers
652 views

Why does TableToTable give ERROR 000732 when published to ArcGIS Server?

I have a python code that uses ArcPy command "arcpy.TableToTable_conversion" to convert a database table to a dbf file. arcpy.TableToTable_conversion The input table is from a SQL Server database ...
1
vote
0answers
182 views

Where are dynamic hyperlinks stored?

A while back (ArcGIS 10.0 times, or maybe even 9.3) a colleague of mine created a mxd with a few hundred manholes and attached dynamic hyperlinks to them (pointing to photos of them). (I know this ...
0
votes
2answers
102 views

Is it possible to alter raster values without Arcpy? [closed]

I'd like to be able to alter raster values, through fairly simple equations, without using the module Arcpy. Is that possible? It seems like gdal, numpy, and rasterlang might make this possible but ...
8
votes
2answers
315 views

Unusual results for geoprocessing speed tests

I have been observing unusual performance with a Python geoprocessing script. The (attached) script performs the following actions: Use a search cursor to look up UTM zone corresponding to polygon ...
3
votes
1answer
596 views

How to save geodatabase attachments to disk with updated exif data

I have a folder of photos without xy information. Using a spreadsheet and lots of patience, I attached the photos to their corresponding GPS point as geodatabase attachments. I now need to produce a ...
0
votes
0answers
247 views

recording pixel value with Python script in ArcGIS 10.x

A scripting newbie question... I'd like to write a Python script for ArcGIS 10.x that would record the pixel value for all pixels in a raster. Furthermore, I'd like to only record a given value once ...
0
votes
1answer
420 views

python script to convert shapefiles in a folder to .kmz's with the SAME shapefile name

Pretty new at working with python and am Stuck...I have a "for" Loop that looks in a folder and, based on a criteria (Homes valued over 500k), creates .kmz's with symbology for each of the .shp files ...
1
vote
1answer
193 views

Arcpy MosaicToNewRaster problem

I can't seem to get the mosaic created for this. import arcpy arcpy.env.workspace = r"C:\projectTemp\FloodAnalysis\output\2yr\raster" arcpy.env.overwriteOutput = True srName = ...
0
votes
1answer
157 views

Split an attribute Python 9.3 in modelbuilder expression? [closed]

Updated question Goal: I need to upgrade an existing Model-builder were the expressions uses Python 9.3. Workproces & Goal: My *.mxd Have 1 attribute field containing the following build: ...