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

learn more… | top users | synonyms

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" ...
1
vote
1answer
72 views

How to use in_field_map parameter of CreateCompositeAddressLocator_geocoding?

I would like to know if anyone who has worked with the CreateCompositeAddressLocator_geocoding method knows how the "in_field_map" parameter is formed?
1
vote
0answers
195 views

Table record counts change when using python tool

I have a python tool which investigates mxd files and clips/projects and exports all visible datasets. Recently I've found that the record counts (even on non-clipped data) is different. It's ok if I ...
1
vote
2answers
217 views

What is the best quick way to populate a field (in ArcGIS 10)? (for rare language/font e.g. Guzrati numerals)

I have polygons generated from shapefile (esri) format I have a tiff image of plot/parcel map(this tiff has plot id scribed in somewhat consecutive in manner) of that polygon feature. I would ...
2
votes
2answers
223 views

How to sort a field according to value range?

I have to check an attribute table field (MY ID NO) of polygon feature against a tiff image of parcel map where fields are written as 12-58,65,87,2565-89999,58648. Example: Parcel map tiff(image) ...
-3
votes
1answer
405 views

How to calculate altitude with given lat,long position?

I'm new to GIS. I want to calculate the elevation of any point on earth giving lat, Long (for example like google elevation Api, or the earth tools: http://www.earthtools.org/webservices.htm # ...
2
votes
4answers
794 views

Error When Trying to Run Slope Analysis Script

I have written a script which is supposed to loop through all the rasters in a folder and perform a slope analysis on each raster and save the results in another folder. I keep getting the following ...
6
votes
1answer
846 views

pyproj, proj4 doesn't return same coordinates after inverse conversion

I'm trying to interpolate between coordinates, as a step in that process I convert latlong to UTM so I can work in metric. I don't get the expected latlong when converting back, reproduced below: ...
2
votes
1answer
209 views

Dissolving stream orders using Python rather than VB script?

I've been trying to dissolve streams, using orders programatically and I came across this VB script. The reason is that if I use ArcGIS dissolve, it merges the confluence of the same orders. I would ...
6
votes
2answers
303 views

Python scripting- compute statistics

I'm new to Python but I've heard is a good way to automate tasks. I have thousands of point shapefiles and need to extract the maximum value of each one and put them in a new table. I'll appreciate ...
0
votes
1answer
1k views

Invalid topology error when running script

I am trying to run a python script and I get an invalid topology rule when I get to the feature to line tool. I am using arc 10.1 but I think the script was written in 10.0. I know it was built in ...
2
votes
1answer
297 views

ArcGIS Desktop registry value under software key? [duplicate]

Possible Duplicate: How can I programmatically get the path of “Python.exe” used by ArcMap I am trying to get some information from registry about ArcGIS. I am using ArcGIS Desktop 10.0 and ...
4
votes
1answer
353 views

Why is the Frequency Tool in ArcGIS so Fast and Efficient?

Why is the Frequency tool in ArcGIS is so fast and efficient? Is it possible to make my own Frequency Tool with similar performance (in python)? I tried to do some based on python dictionary and ...
1
vote
2answers
847 views

Python module “shapely.wkb” missing

I am trying to follow a tutorial from here http://qgistips.spatialthoughts.com/2012/12/find-neighbor-polygons-in-layer-in-qgis.html#more I have downloaded the relevant python script and edited it as ...
19
votes
6answers
3k views

How to make a GIS inventory?

My office will be seeing a big change in its GIS section. This section has been operational since the 1980's and has a huge collection of GIS data (i.e., shapefiles, raster files, data, etc) but never ...
8
votes
2answers
711 views

How to create a solar intensity map for a garden plot?

I'm interested in creating a raster surface that depicts the average solar intensity per unit area for a small garden plot. I imagine the inputs into the model would be the following: Plot ...
0
votes
1answer
692 views

Parsing list object file paths for file name

I am trying to automate a naming process that takes the file name from a list of file paths to name other folders, geodatabases, zipfiles, etc. The problem I'm running into is that, I assume, I am ...
2
votes
1answer
670 views

Update Layer with python, error attribute renderer

I am attempting to automate symbology using python, and ArcMap 10.0. I am importing this script into a toolbox, which I run from arcmap Currently I am getting two errors: 1) LayerObject: Set ...
0
votes
0answers
193 views

How to access the rules of QgsRuleBasedRendererV2 using python?

I can't get access to the rules of rule-based renderer, python says there's no attribute "rootRule", though documentation says it should be (http://www.qgis.org/api/classQgsRuleBasedRendererV2.html) ...
4
votes
1answer
224 views

How to select points within polygon constrained by field?

I've created a polygon from a raster dataset: arcpy.RasterToPolygon_conversion(RasterInt, "C:/arc/poly_int.shp", "NO_SIMPLIFY", "Value") The 'value' field has numbers 1-37, which map to the ID ...
3
votes
0answers
66 views

Programatically reorder layers in qgis? [duplicate]

Possible Duplicate: Sort layers in QGIS table of contents? I write a qgis python plugin ,and i try to connect keybord keys to reoder layers in qgis . Are there a method in qgis allows me to ...
2
votes
3answers
1k views

Couldn't load plugin fTools due an error when calling its classFactory() method

When I start qgis I get a python exception: Couldn't load plugin fTools due an error when calling its classFactory() method ... 309, in _import mod = _builtin_import(name, globals, locals, ...
4
votes
1answer
65 views

How to deliver an additional Python Library with own Plugin?

For my own plugin I need a special library (pyexiv2) to write a comment into a jpg file (PIL does not do that). When I want to share my plugin in one of the repositories in the future, what would be a ...
2
votes
1answer
1k views

How to make a 'Make Table Query' script in ArcGIS?

I am trying to get the make table query running in a script. IF I simply export the script from a model, it works without a problem. However, I do have to change the variable names, as this will be ...
1
vote
3answers
172 views

Applying rotation before FeatureToRaster_conversion?

I am converting polygon shapefiles to rasters. With arcpy I am using the function FeatureToRaster_conversion. gp.FeatureToRaster_conversion(InFeatures, InField, OutRaster, InCellSize) This function ...
1
vote
1answer
279 views

Why won't QGIS let me write to a file?

I am writing a QGIS Python Plugin. I am trying to write data to a file using with open("filename.txt", "w") as f: f.write("file text.") QMessageBox.information(None, "Test", "Written to file.") ...
1
vote
1answer
74 views

Styling map elements based on their properties

I am currently generating tiles based on OSM data using Mapnik and Python and am quite happy using the XML styling but I am at a loss as to how to vary the styling on an element by element basis. I ...
1
vote
1answer
258 views

calculating terrain elevation line

I found too many hits on the web to even see the right direction to go. I have to provide a small proxy server serving requests for a line of terrain elevation data between two coordinates. It is a ...
3
votes
2answers
568 views

Create Table (Data Management), Unable to Save to Feature Dataset

I'm trying to create a table in a feature dataset and get the error: Parameters are not valid. The code works when I write only to the gdb. Am I not allowed to write to a feature dataset or am I doing ...
2
votes
0answers
76 views

Arcpy: Counting number of segments where altitude rises or decreases based on DEM?

I have a DEM of Sudety mountains and linear .shp with several hiking paths. I'm trying solve one problem: How to make a script based on input DEM and input line shapes that would count number of ...
2
votes
2answers
1k views

How to re-project the EASE (Equal Area Scalable Earth) grid with a ~25 km cylindrical projection to WGS84 0.25 degree?

I have nc files. from the metadata ,the projection is cylindrical and the resolution is 25 km: indrical NC_GLOBAL#ease_resolution=25 ...
10
votes
1answer
454 views

Is there a way to use ArcPy and PyGRASS tools in the same Python script?

I´m trying to write a Python script, that includes different ArcPy / ArcGIS commands and want to combine it with GRASS geoprocessing tools. Unfortunately importing the GRASS libraries doesn´t work. ...
2
votes
0answers
366 views

problem in read_shp() to generate a networkx.DiGraph from point and line shapefiles

I tried to use the Python module networkx and specifically the read-shp () to ensure the transformation of a graph to a file sharpe but I always find this error noting that this is the first time I ...
1
vote
1answer
951 views

Slow tiles generation for leaflet js

I am currently trying to use leaflet to display images (jpeg or png) instead of maps, I have followed the instructions on this blog post to create the tiles using gdal2tiles.py here is the code we ...
1
vote
0answers
344 views

Morphometric and hydrological analysis with Python and ArcGIS

How can I write a program using python, that calculates some morphometric and drainage indexes using numeric attributes in the attribute table with ArcGIS? The data I have in my table are area and ...
2
votes
1answer
388 views

Saving rasters in a Python for loop fails only on last iteration

I'm a newcomer to Python scripting and I'm having what's seemingly a very simple problem but one that's stumped me enough to send me posting to these forums for the first time. I'm using a for loop ...
19
votes
1answer
7k views

How do I get the pixel value of a GDAL raster under an OGR point without NumPy?

I am working on a computational model of the abundance of wild pollinators across a landscape. The model itself is complete, and I am now struggling with a post-processing step. I have my GDAL ...
3
votes
1answer
965 views

How Can I Clip a Basemap with an ascii file?

I'm new here but with many doubts.I work with ArcGIS 10 and I add a basemap layer from the ArcGIS online but I just want the continent coverage to my map,and the area that I wanted I have an ascii ...
7
votes
1answer
196 views

Reading shapefile tolerance in python

Is there a way to get the cluster or xytolerance from a shapefile into python using free/open source software? I'm currently using the pyshp package and can't see methods for this anywhere. Also it ...
9
votes
1answer
1k views

How to prompt for user input during ArcGIS Python script tool execution?

I want to prompt for user input after all parameters have been specified for the script tool. I've seen many forum posts and stackexchange questions about this, and they all conclude that it's ...
4
votes
1answer
398 views

Calculating Service Area for each feature in feature class using ArcGIS Network Analyst?

I am trying to implement a python script that creates in a first step a Network Dataset and that in a second step gets a point feature class for which the tool should create for each feature in the ...
0
votes
1answer
757 views

Invalid data source issue when trying to use lyr.replaceDataSource

When I run the following code in an open mapdocument, it executes and does what I want it to, replace the datasource of the Indianacounty_shp layer with the County_Boundary feature class. import ...
5
votes
2answers
1k views

Create subsetted multiband image from multiple files using gdal or python?

Given a time series of an area (single band), how you subset them and combine the subsets into a single multiband file (eg. vrt) using gdal or gdals' bindings for python?
0
votes
1answer
2k views

Runtime error <type 'exceptions.TypeError'>: 'unicode' object is not callable

I'm working on a script to replace a layer datasource, and I'm getting an error I cannot figure out. When I run the following code in a python window inside ArcMap, I get a "Runtime error : 'unicode' ...
3
votes
2answers
1k views

Checking if a point is in a circle

I'm using shapely to check if a point is in a polygon, but I also need to check if a point is in a circle. Checking the docs, I don't see how do that... Maybe I could create a circle? circle = ...
5
votes
1answer
658 views

How to interpolate polyline or line from points?

We have a lot of points, but with no timestamps to create ways. For now, I've created a sample of the way the points traverse. It's in here, this way can be compared with a Google one. I haven't ...
2
votes
2answers
537 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 ...
3
votes
1answer
815 views

Relative Hyperlinks in QGIS without end user's software dependency

I am trying to work out how can I add a relative path hyperlink to a shapefile in QGIS. I have found several great tips on various forums but only one has so far potential to be implemented. The ...
4
votes
2answers
298 views

“Hello World” introduction for ubuntuGIS?

Is the GIS equivalent of "Hello world!" = "You are here?" I installed ubuntuGIS recently, but I'm a little overwhelmed navigating the network of many packages to begin accomplishing simple things. ...
5
votes
1answer
209 views

GDAL/OGR Determine if dataset is raster or vector

I am writing a script using the GDAL/OGR Python bindings and it is possible that the user may provide a raster dataset or a vector dataset as input. If the user provides a vector dataset, the script ...