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

learn more… | top users | synonyms

6
votes
1answer
435 views

Histogram Matching Python

I am trying to do histogram matching using Python to improve the mosaicking process of multiple overlapping rasters. I am basing my code on that found at: ...
2
votes
2answers
81 views

Multiply Rasters from two subfolders

I have two subfolders with Rasters named with the same parameters. I am trying to figure out a way to Loop a multiplication of each raster into a folder with its pair of the other subfolder. Looks ...
0
votes
1answer
183 views

Cause of ArcGIS Raster Calculator Parsing error <type 'exceptions.SyntaxError'>: invalid syntax?

This is my first attempt at much code in the raster calculator of ArcMap 10.0, but I think I'm close. I have nine rasters (A-I below) and want a new raster with value = 1 where any of the rasters ...
1
vote
0answers
20 views

Using ST_AsText() without session object

I am storing gps locations as feature in a restfull api. The locations are linked by foreign key to "place" objects in a flask/sqlalchemy/geoalchemy framework. I would like to print a user readable ...
1
vote
1answer
74 views

“ImportError: no module named qgis.core” in standalone app

I am developing a standalone app (in Windows) using libraries from QGIS. When I run the script in cmd, I get the ImportError: no module named qgis.core. I read in several posts that I need to ...
0
votes
1answer
36 views

Missing labels when using arcpy.mapping.ExportToPDF

I have a script that is run on a daily basis to export a folder of MXDs to PDFs and I have noticed some peculiar behavior. Many, but not all of the PDFs will be missing some labels or their label ...
0
votes
1answer
31 views

Using arcpy to separate surface objects mixed into one file into different shapefiles

I'm trying to do some elementary tasks with arcpy and ArcGIS, namely separate some information in a file into separate shapefiles based on values for a field. My initial file has a field called ...
2
votes
2answers
145 views

UpdateCursor on two tables with value from SearchCursor

I try to keep myself busy and learn some python. I have two tables test2.shp (for searchcursor) test.shp (needs to be updated) I want to update test.shp with a UpdateCursor. The column which ...
2
votes
2answers
97 views

Batch add column to layers and add layer prefix in PyQGIS

I have a number of layers each with a prefix like A_ B_ and so forth. I'd like to: Add a column (text) called DENSITY to each shp Add the prefix letter from the layer name to each row in the ...
3
votes
1answer
139 views

ArcGis Python Security Concerns

I’m working for a large organization which are using ArcGis Desktop on some (windows) clients. The security requirements for the system is very strict, and a part of that is prohibiting the users the ...
1
vote
1answer
30 views

Trouble calculating attribute values in pyqgis [closed]

I recently posted another question about calculating unique IDs for features using python and qgis here Calculating unique ID value for a field in QGIS using python. I am trying to use the same ...
5
votes
1answer
145 views

How to open QGIS attribute table in form view using Python?

I've got a plugin where I load some layers. After the last layer loaded, I show the attribute table. However, I would like to select the form view and not the table view (see image). How can I do ...
1
vote
1answer
65 views

Calculating NDVI with Rasterio

I'm using Rasterio to read GeoTIFF files from Landsat 8 and calculate NDVI into a new GeoTIFF file. My code looks like this: import numpy import rasterio import subprocess with ...
7
votes
2answers
2k views

How to Refresh open attribute table in ArcMap using ArcPy rather than ArcObjects?

I am running the Calculate Field GP tool from the arcpy in ArcMap. When a row value is updated, it does not show up untill I use the Reload Cache tool (Table Options > Reload Cache). How do I refresh ...
1
vote
1answer
67 views

Running Python toolbox tool both standalone from other tools within same Python toolbox?

I have created a python toolbox for work that has two tools inside. These tools use the toolbox template and so have a init, getParameterInfo, etc. I want to be able to run one tool both stand alone ...
0
votes
1answer
168 views

Copying XY to clipboard using Python Add-In (ArcMap)

I am trying to create a simple tool to copy XY from where i click the mouse button. What i have done already is: create add-in files using Python Add-In Wizard create toolbar and tool inside it in ...
2
votes
2answers
165 views

How to calculate monthly maximum values from daily data in a 365 band grd?

I have a MINT2005.grd file which has 365 bands representing days of year from 2005. Each band is a 32x35 matrix representing the temperature at different latitude and longitude pair. How do I ...
1
vote
1answer
618 views

Linux Client for OpenGTS server

I've installed an OpenGTS server and successfully linked it with the GPS2OpenGTS_Pro application for android. The buses I'm tracking offer 3G with a wireless AP, however whenever the bus stops, the ...
2
votes
1answer
348 views

Can you access ArcScan from a python script?

I'm using a python script with ArcPy in ArcGIS 10, and want to access some ArcScan functions. Is this possible? I can't see anything in the documentation about it, but I can't really find any very ...
0
votes
1answer
55 views

Is it possible from the list or np.array, updating the table using arcpy.da.UpdateCursor?

I have created array. Is it possible to update an existing table if the number and order of columns and rows identical? In the array do not have a column headers. for example: array = ...
2
votes
1answer
127 views

How to print 2 maps in the same PDF page using PyQGIS?

I am trying to generate a one-page PDF document containing the following items: 1) On the left half, a full state map (say, New York State). 2) On the top right quadrant, a map that zooms in to one ...
2
votes
2answers
59 views

Is it possible to sort the features by an attribute programmatically?

I have written a PyQGIS script that loops through an attribute table and finds duplicated values in a target attribute field. My technique requires that the attribute table is first sorted by the ...
3
votes
2answers
71 views

How to access checked layers in PyQGIS?

qgis.utils.iface.activeLayer() accesses the current active layer, but I would like to get all checked layers in the QGIS ToC. How can I do it?
3
votes
1answer
64 views

How to order points along a line?

I have a bus line with bus stops along the route. The sample data looks like this: LINESTRING (-70.57243474999999933 -33.37452829000000065, -70.57308152999999606 -33.37470645999999874, ...
3
votes
2answers
115 views

How to calculate a value based on another field under arcpy.da.UpdateCursor and cursor.updateRow?

I am trying to do a calculations under arcpy.da.UpdateCursor. I want to calculate the values for the records that have FIPS=06037. I want to take the values from "kg" column and populate the new value ...
6
votes
4answers
982 views

Code to open web browser crashes ArcMap when run from a Python add-in

I created an ArcMap Python add-in using ArcGIS 10.1 to open a webpage after the map is clicked. It opens the page but then ArcGIS crashes. import arcpy import pythonaddins import webbrowser class ...
2
votes
1answer
43 views

PyQGIS jump to row N

Is it possible to go to a specified row number (or perhaps the next/previous) row using PyQGIS? This code gets the value of the first attribute in the first row of the table features = ...
3
votes
2answers
177 views

VB Script or Python for Field Calculator when using greater/less than?

What I have: I have 4 fields for my address ranges. L_ADD_FROM, R_ADD_FROM, L_ADD_TO & R_ADD_TO (all in DOUBLE field type). What I need: The lesser value from L_ADD_FROM and R_ADD_FROM to ...
3
votes
3answers
1k 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 ...
0
votes
1answer
127 views

40k+ points from Shapefile, clustering them and populating on Google Maps/Webapp

I just started a new Python/Django webapp job related to GIS (and i'm totally new to GIS). Here is my task and what i have found so far, and i am stuck as to finding out what's the best way or the ...
0
votes
0answers
28 views

How to query a WMTS layer with a point layer in QGIS?

With the "identify" tool it is possible to query WMTS layers. Then, you can copy the acquired attributes and afterwards the clicked coordinates, then create a CSV file with attributes/coordinates. Is ...
1
vote
0answers
44 views

How can I extract the measurelength, highmeasure and lowmeasure of routes in ArcGIS

I'm trying to create road Routes in ArcGIS, which traditionally have been created using workstation. I've managed to write a Python script to do this, however, unlike ArcInfo, the attributes ...
3
votes
1answer
310 views

Finding optimal network with python and network analyst in arcgis 10.0

Hello GIS StackExchange, I'd like to know if there's any way to find the shortest distance across the diagonals of a polygon. I am creating a network analyst solution to determine the most cost ...
1
vote
0answers
23 views

Odd Error in Metadata Script

So I developed a script that traverses through a root directory, and any subdirectories and catalogs layers files. Then exports these layer files' FGDC metadata and outputs them to a folder. The code ...
3
votes
2answers
61 views

How to trigger a Python script from an ArcGIS Desktop “save” operation?

Similar to the question Run a Python Script when opening an MXD, I'm looking for a way to run a Python script that is triggered when the "Save" command/button is pressed. I've read the Creating a ...
5
votes
1answer
156 views

How to create a python plugin in QGIS 2.4.0?

While I am learning to create a python plugin in QGIS 2.4.0, I have followed the proceedure expalined in ...
1
vote
3answers
57 views

Calculate field management not working within customized Python tool in ArcGIS

I am trying to make a custom made Python Script tool which I would like to use in Model builder. I am using the arcpy.Calulate_field_management function (simple reclass function as an example), but it ...
1
vote
1answer
38 views

ArcMap: Unknown Spatial Reference fix with GDAL/Python

I have a .tif file that I created with Python and when I open it in ArcMap it says it doesn't have spatial reference information. Is there a way to set the spatial reference with GDAL if I have a .prj ...
3
votes
2answers
165 views

Spatial Classification

I am not sure what to call my problem. From this image, I have a set of classified polygons (colored). I would like a way to classify the white polygons in the image. I don’t have any attributes to ...
1
vote
0answers
83 views

How to create an expression using the arcpy.SearchCursor? [closed]

Example: On the basis of the table which is not defined by the number of columns and rows, you need to evaluate expressions as shown, using the arcpy.SearchCursor .... First, create a list of ...
0
votes
0answers
23 views

How to merge grd files programmatically, ideally in Python?

We currently get LIDAR data from the field in grd format. Sometimes this data is overlapping. This data gets fed into a python program that does some analysis. What I'm after is some way of ...
5
votes
1answer
809 views

Python equivalent of gdalbuildvrt

Is there a way to perform the same task as the gdalbuildvrt utility using GDAL Python bindings? So far I have not found any way to do this other than creating a vrt of a single dataset and manually ...
1
vote
0answers
33 views

QGIS actions - get currently selected attribute column (with python?)

I'm trying to open files dynamically through actions. A simple example would be: Action name: OpenLink from PyQt4.QtCore import QUrl; from PyQt4.QtWebKit import QWebView; myWV = QWebView (None); ...
0
votes
1answer
48 views

Running Grass7 in a conda environement

I have Anaconda Python Distribution (Python 3.4.3 :: Anaconda 2.1.0 (64-bit)) installed on Ubuntu 14.04. I would like to create a conda environment for running Grass 7. To do that I used: conda ...
0
votes
1answer
60 views

Inserting row into table and then updating it in sum of each column using ArcPy?

My problem is the following: I have a table with an unknown number of columns (n) Use arcpy.ListFields create a list of columns and generates Field name list Then, using arcpy.da.InsertCursor create ...
2
votes
1answer
172 views

Is it possible to run a search cursor on a multipart feature class?

I have a multipart feature class of cities. I want to run a search cursor on a cities feature class where: one city is selected and a select layer by location is run to select all parcels that are ...
0
votes
2answers
77 views

ArcGIS field calculator: Calculating a new field using multiple fields with python

Basically, I have a large polygon data set that has soil data for an area and I want to assign numeric values to the data based on soil texture text values (e.g. SIL=4, CLY=7, ...). However, the soil ...
0
votes
0answers
23 views

Password Caching in ArcMap using Oracle accounts

The issue I'm looking for help with is that we have several Oracle accounts through which we host data. For example we have ACCOUNT@ORAC and ACCOUNT1@ORAC (these are made up names as I don't to ...
5
votes
2answers
237 views

Cancel event in Python Gp tools with arcpy

In GP-tools written in C# (ArcObjects) we can listen for the Cancel-event and react on that. I wonder if there is anything similar in Python GP-tools applying arcpy. I googled that issue and the only ...
1
vote
1answer
18 views

Trouble updating attributes of a layer based on another layer's attributes if features intersect in python, QGIS

I am trying to assign attributes from features in one layer that intersect features in another layer. When I run my code, I do not get an errors but it does not update the attribute field as I am ...