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

learn more… | top users | synonyms

2
votes
2answers
352 views

How can I close the identify result dialog programmatically?

I am trying to close the identify result dialog which appeared when I select to identify a feature from the active layer. Does someone know how can I do it?
2
votes
2answers
584 views

Error when looping Polygon to Raster with python through different fields ArcGIS 10

This is the first time I have posted a question on here so please have some patience! I've been trying to create a python script to loop through numerous shapefiles converting various fields to raster ...
2
votes
2answers
532 views

Python-OGR: nested loop only loops once

I have a polygon file buffer_stands.shp with 20 features and I have a line file OSMroads.geojson with 9 features. I expected the script to loop through the buffers and print out for each buffer loop ...
2
votes
1answer
573 views

How to extract vegetation polygons from Imagery using Python?

I have been trying to roughly polygonize vegetation areas from Imagery (RGB). I tried using Iso Cluster Unsupervised Classification (Spatial Analyst) but it is still too detailed for my purpose. I ...
2
votes
1answer
432 views

Embedded python script module import issue

I have two ArcGIS python scripts. One is a primary python script that calls a number of functions from a second python script. Everything works well, i.e. primary tool calls and uses functions from a ...
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 ...
2
votes
2answers
377 views

Python/GDAL: Best way to read in raster for use in A* search algortihm

I want to create a least cost path program with Python and GDAL since I couldn't find any open source cost path additions for use in Python. As recommended here the best approach is to use A* ...
2
votes
1answer
1k views

How to install QGIS 2.0 on Mint 15?

Trying to use QGIS 2.0 on Mint 15 (based on Ubuntu 13) I installed, uninstalled and re-installed QGIS from various repositories (ubuntugis-unstable, then qgis.org/debian raring - not clear for me ...
2
votes
1answer
1k views

Subsetting a curvilinear netCDF file (ROMS model output) using a lon/lat bounding box.

I am attempting to get salinity data from one of the ocean models in our region and have run into a roadblock. This is ROMS model output in NetCDF form, but the grid is curvilinear (lon and lat arrays ...
2
votes
1answer
212 views

Is there a way to reference a selection in model builder or python?

Is there a way to reference a selection in model builder or python so that I can edit my python script to only run if there is a selection. And even better if there is only one row selected? I made a ...
2
votes
2answers
508 views

Trying to use Python code and user defined inline parameters to build a SQL expression

I am trying to build a tool where the end user can run a model and make a selection from a feature class based on 3 different fields, calculate the total shape length, & export the results to a ...
2
votes
1answer
917 views

Updating the value of a field with pre-defined values

I am trying change the value of a row for a field whose type is "short" in arcmap using a python script. with arcpy.da.UpdateCursor(FILENAME, FIELD) as cursor: for row in cursor: row = ...
2
votes
2answers
436 views

How to add non spatial sqlite table to map using python in qgis?

Is it possible to load pure sqlite 2.0 (without spatial extension) table into qgis using python? I would like to get the same result as in the case of using Add Vector Layer -> All files. I have ...
2
votes
2answers
456 views

Erase method using ogr

I'm trying to find the erase geoprocessing method as i am familiar with in ARC, but i can't seem to find it in OGR using Python. I have looked through the module ogr documentation but perhaps i missed ...
2
votes
1answer
900 views

Distance of points in polygon to nearest polygon edge

I have a polygon shapefile with about a dozen polygons (Area.shp). I have an equation that calculates a value for points (or small areas) as a function of distance to the edge of these polygons. I ...
2
votes
1answer
197 views

Stumped on python looping and data dictionaries

I have been stumped getting my python code to work for a while now. I've been using python for the last several years on and off for simple things, but I am far from an expert in python. I am trying ...
2
votes
1answer
4k views

How do I select by location based on selected features with Python?

Is there a way with Python to do Select by Location with the Use selected feature option? I'm trying to automate a process and am running into barriers. I need to be able to select the 2 lines ...
2
votes
1answer
1k views

How to perform a cascaded combine from python script in QGIS?

I'm working with a large dataset, and I need to be able to combine up to 400,000 polygons. I have been using the QgsGeometry::combine(QgsGeometry* geometry) function and simply looping over all ...
2
votes
2answers
158 views

How to reproject all featureclasses in a database to another database?

I am very new to arcpy and have ArcGIS 10.0 I would like to make a script where I would only have to change two variables in a script (database_input and database_output) so that I can reproject all ...
2
votes
1answer
2k views

Reverse geocoding to get city, state, country data separately?

I have longitude and latitude data. I need to reverse geocode to get city, state, country data separately. I have checked google service. Google provides complete address to a micro level. ...
2
votes
1answer
168 views

How to change VB syntax to Python in Field Calculator?

I want to calculate a field, this is my working example in VB: "C:\Users\Admin\Desktop\Data\" & !Work! & "\1000_\" & !Data! & ".pdf" I'm slowly changing all my scripts to Python, ...
2
votes
2answers
462 views

How to determine the coordinates of the corners of the raster layer

How can we determine the coordinates of the corners of the raster layer? The image opens function QgsRasterLayer.
2
votes
2answers
392 views

Weird File Path Returning from Print Statement Causing ERROR 000732

I have created a python script in ArcGIS 10.1. My first parameter type is a 'Feature Layer' set to as a multivalue input. I was running my script quite well off of data from one of my drives (my ...
2
votes
1answer
536 views

Am I over-complicating adding and then calculating fields? Suggestions for a better way? (Python)

I'm attempting to re-position a question I asked yesterday here. I am starting to think that the way I wrote my code just isn't the best way. Am I overcomplicating things? I know for a fact that ...
2
votes
2answers
694 views

Save a layer for each unique value using Python in QGIS

I have a map, in which each region is defined as an attribute of the map (I can select each region from the attribute table). I am trying to create new layers each one containing only 1 region using ...
2
votes
1answer
743 views

How to use QgsOverlayAnalyzer class in pyQGIS?

I am very new to Python scripting in QuantumGIS, so bear with me please. I am trying to use the QgsOverlayAnalyzer class in QuantumGIS (http://www.qgis.org/api/classQgsOverlayAnalyzer.html#details) ...
2
votes
1answer
997 views

How to refresh a Combobox Add-in in ArcGIS 10.1

I need to refresh the value in the combobox after an event from other tool in the same add-in (ARCGIS 10.1. Python). I've changed the value of the combobox with the code: combobox.value = val and I ...
2
votes
1answer
347 views

Updating and inserting a field within the same loop

I am working with an empty feature class where I am writing a script to create polylines within the featureclass and attach an associated text value with the polyline. I am using a for loop as noted ...
2
votes
1answer
1k views

how to install networkx python in QGIS

How to install and how to compile networkx python on QGIS? I always get an error that no module named networkx in the console of QGIS I copied the file networkx in C:\Program Files\Quantum GIS ...
2
votes
2answers
310 views

How to automate performing analysis (intersection) between objects within different feature classes, for large number of objects?

I have two tables of roads - the same roads, from different sources, and am trying to compare how similar their positions are. I'm wanting to buffer each road in one of the feature classes and ...
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, ...
2
votes
1answer
644 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
1answer
527 views

Can you script publishing a GP service?

I have a geoprocessing service in ArcGIS for Server 10.1. I would like to share this service with a number of external organisations that also have ArcGIS Server 10.1. i.e I dont want them hitting ...
2
votes
3answers
2k views

An error has occured while executing Python code while doing home range analysis (Qgis 1.8.0)

It keeps appearing: Traceback (most recent call last): File "C:/Users/Worten/.qgis//python/plugins\HomeRange_plugin\HRplugin_dialog.py", line 80, in QObject.connect(self.buttonBox, ...
2
votes
1answer
1k views

Batch intersect from a text file using Python (ArcGIS)

I am attempting to modify my code below which reads from a text file of shapefiles and executes an intersect and outputs them to a user specified location: import arcpy, sys, os from arcpy import env ...
2
votes
2answers
1k views

GDAL Clipping (hack)

I am trying to use the python script from http://www.sparkgeo.com/labs/2010/07/07/gdal-clipping-hack/ to clip a raster based on many shapefiles. The script is using a mysql database an i would like to ...
2
votes
2answers
1k views

With Python and ArcGIS, how do I add features to a featureclass?

In ArcGIS Python code: CopyFeatures_management copies features to the same featureclass - that means overwrite to the feature class. How can I add new features to one featureclass or namely different ...
2
votes
2answers
324 views

Is code completion available in QGIS Python Console?

I wanted to know if the autofill functionality is available in the python console for QGIS as the one available in ArcGIS Desktop's Python Console or am I supposed to hardcode it.
2
votes
1answer
829 views

ArcPy - Dataframe rotation - Getting a polyline length wise in the dataframe in layout view in ArcMap

EDITED: I'd like to rotate my dataframe until a polyline fits dataframe-width wise in the layout view of a .mxd. If the lines runs north/south I would like the dataframe rotated until the line is ...
2
votes
2answers
3k views

Help with ArcGIS10 Python and SearchCursor reading Null values in field

I'm using a Python script to count the number of unique string values in a feature class. My script works if there are values in the field but if there are any NUll values in the field it crashes. Any ...
2
votes
1answer
558 views

How can I, in a general way, rescale a raster in ArcGIS Model Builder

What I want to do is rescale a few rasters from whatever they are now (i.e., values between -99 and +1000 OR -256 - +256 OR 0 - 3.23) to a normalised value (i.e. values between 0 and 1). I would ...
2
votes
2answers
3k views

How to automate point-in-polygon queries in ArcView?

Layer1 = point shapefile (multiple points, geocoded) Layer2 = polygon (1 attribute; (1,2,3...)) I want a non-gis person to be able input a single address and be told the attribute of the polygon ...
2
votes
2answers
300 views

How to use any ArcMap tool from a python script?

I'm trying to develop a tool that uses some of arcmap default tools like the Sketch Tool or Trace Tool. Is it possible? Let me explain a little more, I need to draw a feature, ask for user input, and ...
2
votes
1answer
65 views

PyQgis: During editing, add new feature attributes

Workflow: Start edit mode -> New feature is created -> On the console (pyqgis), new attributes are set -> Exit edit mode. Problem: The attributes are not saved. Only after completion of the editing, ...
2
votes
1answer
40 views

Why is the first element of my WFS layer None instead of a feature?

I have made a WFS connection and got a layer from it. I want to use the first feature from that layer. driver = ogr.GetDriverByName("WFS") wfs_url = "WFS:http://example.com/wfs" wfs = ...
2
votes
1answer
56 views

Suppressing rasterio verbosity

I have the following script using rasterio: if not os.path.isfile(out_dir + os.sep + name_adm1 + '.tif'): with rasterio.drivers(): with rasterio.open(constants.ras_cmask) as src: ...
2
votes
1answer
77 views

QGIS corrupted processing tool

I've been dealing with various problems after re-installing QGIS and other GIS-related software (see QGIS corrupted plugins upon reinstallation). The solution worked for a bout a day and now I lost ...
2
votes
1answer
53 views

How to multiply a user-defined input parameter with a field value within the graphical modeler?

I have the same Question as asked in this post (How to connect input parameter to the field calculator algorithm within the graphical modeler?) and the answer from Joseph already helped a lot: ...
2
votes
2answers
101 views

PackageLayer_management failing within script, needs to run outside of an ArcMap session

I'm automating a workflow that requires outputting a layer package, it throws the general ERROR 999999 when executing and breaks. It does, however, generate a folder in the correct directory, with ...
2
votes
1answer
72 views

Getting script file name as object in Python? [closed]

I'm quite new to Python and programming in general. I have created a few scripts using Python and as part of my framework/best practices for my scripts, I'd like to display when a script begins and ...