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

learn more… | top users | synonyms

8
votes
5answers
1k views

Will there be a Python API for ArcGIS Server?

It seems like an obvious next step given the big push towards Python in ArcGIS Desktop at 10.0. Are there any limitations of Python compared to javascript, flex, etc that would prevent ESRI from ...
0
votes
2answers
40 views

remove values in list using logical operators

Just would like to ask how can I masked or remove the values in my list based on logical operators. Here is the case: I converted an NDVI.img file into a 1 dimensional array. and as part of the ...
0
votes
1answer
39 views

Problems with Geometry collections. (Fiona: write output of Multipoints and Points)

Updated: I am getting the following error on this function: ValueError: Record's geometry type does not match collection schema's geometry type: 'Point' != 'MultiPoint' The problem is that I have ...
2
votes
2answers
187 views

Running Python script in Task Scheduler— Script will not run

I have a script which I wish to run on a regular monthly basis without me needing to open and run the it. I have referred to every web page I can find on Windows Task Scheduler and followed their ...
1
vote
2answers
75 views

Using Python to export multiple layers to .pdfs

I have found a solution to the issue below. The problem had to do with coding GetParameterAsText without actually using it as an input in a tool. That was the reason the raster layers wouldn't turn ...
2
votes
2answers
301 views

Python problem visualization .tif image

I am trying to read a hyperspectral image in Uint16 .tif format and visualize it. I use GDAL in Python to read the data and everything works perfectly when I use matplotlib to show one band of the ...
3
votes
1answer
243 views

How to clip a Oracle query layer directly in a Python script?

I have a query layer directly running from ORACLE database in a MXD. I want to clip this layer with a layer being created in another MXD in a Python based tool. My question is to do this clip ...
0
votes
0answers
19 views

Bulk Converting Ordnance Survey to Google Map lat lng

Apologies if this is a simplistic question, I am learning! I am attempting to use pyproj with Python. The coordinates I wish to convert came from the Ordnance Survey Integrated Transport Layer. I ...
1
vote
1answer
58 views

Python encoding error in ArcPy

I am writing a JSON file, as string input, using python, into a .json file in the file system. However I am getting the following error sometimes and I don't know why. Any ideas? ...
0
votes
2answers
45 views

How to measure distance using shapely

Reading the documentation...it looks like I should have this right, but I don't. The example is: Point(0,0).distance(Point(1,1)) 1.4142135623730951 Based on this example I wrote this block ...
2
votes
2answers
59 views

What is the “best” way to automate the GIS slope and aspect algorithm?

I am currently obtaining slope and aspect of roofs using ArcGIS as follows: Mosaic several Lidar .asc files to new raster Extract by mask using Mastermap building outlines to get hold of just the ...
7
votes
1answer
6k views

issues importing arcpy in Python IDLE

I would want to use arcpy in my python code. I am able to import it through the ARCGIS desktop python console. But I am not able to import it into Python IDLE. I get the following error Traceback ...
1
vote
0answers
51 views

Loop exits after 100 count while creating domains using python [on hold]

I am trying to create domains on feature class or table class by reading a table which holds all look up information. Following is example code for what I am doing: tableIndex = ...
1
vote
1answer
124 views

using postgis function in pyqgis

I'd use Postgis function in Pyqgis script. This code: qsql="(select gid, id, St_Buffer(the_geom,200) as geom from archi)" uri = QgsDataSourceURI() uri.setConnection("ip", "port", "x", "x", "x") ...
2
votes
1answer
133 views

Accessing CRS info with pyqgis in a standalone application

I'm currently having trouble accessing and defining CRS information in a standalone application. I've read the API documentation and have tried it in the python console in QGIS 2.2.0. The python ...
2
votes
1answer
36 views

Why does projection WKT differ between osr.SpatialReference and gdalinfo

I am writing a simple geotiff using the GDAL python bindings, but the projection WKT reported by gdalinfo differs from the source. The script from osgeo import gdal, osr import numpy # Create ...
0
votes
0answers
16 views

Programmatically Add Layer Definition Files (QGIS Python)

I would like to add the same Layer Definition File to QGIS desktop multiple times with a Python script rather than having to make multiple clicks each time I want to add that layer. Can someone ...
0
votes
1answer
49 views

ArcPy Script to update LYR data source not identifying all layers in folder

I've put together a script that for the most part works as intended - to walk through a series of folders, find SDE layer files and update their data source to direct connects if the layers support ...
2
votes
0answers
45 views

How to use the multiprocessing module on python to execute a spatial join?

I have several shapefiles of the entire Amazon rainforest each with millions of polygons (they are fairly detailed and take some time to draw on ArcMap). I need each of these shapefiles to be ...
2
votes
0answers
40 views

Converting csv that it can be displayed

I have this Excel sheet here (far to right side there is a box where the file can be downloaded, the very first one). http://www.bafg.de/GRDC/EN/02_srvcs/21_tmsrs/211_ctlgs/catalogues_node.html Im ...
2
votes
1answer
132 views

How do I link OGR with spatialite?

Problem & Evidence: I'm trying to create a spatialite file with the code snippet below. However the outDataSource object gets assigned NoneType and all subsequent (unshown) code fails. I'm ...
1
vote
1answer
44 views

How to call 'shapely' .coords on a fiona filtered list?

I have to filter a shapefile by attribute then compare it to another shapefile using .coords...however I am a bit perplexed how to do this as my filter is returning a list type...(and .coords requires ...
0
votes
1answer
22 views

ImportError: no module named site

I am developing a standalone app using qgis librarys and run the python inside qgis. I import the librarys with no errors but I have one last error: ImportError: no module named site. I define the ...
1
vote
1answer
392 views

Calculating Stream Length-gradient Index (Hack Index) with the field calculator

I'm trying to calculate the Stream Length-Gradient Index of a stream (example data below). Just started to get into python, but am struggling with a solution to this problem. I have 42 streams to run ...
0
votes
3answers
64 views

Using arcpy outside the arcmap

I want to use some arcpy function outside arc map i used following code but it ends up with error no module named arcpy can anyone help import sys sys.path.append("D:\Installed ...
3
votes
1answer
281 views

How come “Zoom to Polygon” always puts me in the middle of Nebraska?

This is an issue related to this question I'd asked a few days ago. I've got that code working, to an extent. Every time, though, that I double click on one of the Polygons in the attribute table that ...
0
votes
0answers
38 views

Arcpy save command giving errors

I've been trying to get this to work but for some reason it just keeps failing. Here is a very simplified version of my script: import arcpy arcpy.CheckOutExtension("spatial") #Spatial Analyst ...
0
votes
1answer
44 views

Print values of a field in arcgis using python [closed]

I am pretty new at python scripting and I am trying to print values of a field to screen. My code is that import arcpy layer = arcpy.GetParameterAsText(0) arcpy.AddMessage('Layer Name : '+layer) ...
12
votes
3answers
755 views

Industry-based examples of using Python in geoprocessing and arcpy

I was asked recently by some students what types of operations GIS Analyst/Developer are usually trying to automate by using Python with geoprocessing in ArcGIS and ArcPy site package. I guess it ...
0
votes
0answers
14 views

script to Determine Parcel Neighbors [duplicate]

I have a layer of Plots Urban planner I want to fill the neighborhood north and south, east and west of each piece fields i need script to Determine Parcel Neighbors using PYTHON
1
vote
1answer
45 views

Setting PythonPath

I have a question which might be a little bit odd, but I want to make sure I get that right.... I want to use this script here Ortho Projection produces artifacts (unfortunatley I cannot download the ...
1
vote
0answers
26 views

Kartograph - converted svg not the same as before

I'm using the following script to convert a map from NaturalEarthData. But the result is very sobering. It seems like something went wrong with the scaling, but I got no errors when executing the ...
2
votes
1answer
120 views

Python Error while installing Cartogram

I am quite new to GIS and qGIS especially. I would like to install cfarmers cartogram plugin. From what I understand, this has to be done manually - so I 've followed the instructions from ...
14
votes
5answers
375 views

Working as a Python Development Project Team with ArcGIS?

We have a development project in Python (ArcGIS 10). This project involves a mixture of toolboxes, map templates, layer files, file geodatabase templates (acting as templates being imported into a ...
0
votes
0answers
20 views

Average cell value within polygon for each raster in a mosaic dataset - ArcGIS

I am working with a mosaic dataset that has the following properties and attributes... It consists of 365 rasters, each representing a day of the year. I am needing to get the average of the ...
2
votes
1answer
2k views

How to fix python ImportError: DLL load failed: The specified module could not be found?

During import fiona, I am getting following error. Error is different from other user so I asked this question may be it is repeated. import fiona Traceback (most recent call last): File ...
1
vote
1answer
33 views

The truth value of a raster is ambiguous. Invalid use of raster with Boolean operator or function.

I'm having trouble with an if/else statement in my for loop. I have a raster Vmax with a constant value, another raster V which has a value of 0 in the beginning and will change during execution of ...
2
votes
1answer
232 views

Setting up pycharm and pyqgis

I wanted to set up PyCharm the way this blog post explains: Setting up PyCharm for PyQGIS and Qt. I did everything as instructed but when it comes to adding project interpreters I got this error ...
0
votes
1answer
28 views

Converting JSON to GeoJSON MultiPoint using `ogr.CreateGeometryFromJSON`

I am trying to convert a JSON array to a geom in Python, by first converting it to a GeoJSON. I have point coordinates in a JSON array, [{'lat':100.0, 'lng':0.0},{'lat':101.0, 'lng':1.0'}] that ...
0
votes
2answers
51 views

Model Builder Error: Python cannot be found

Building some models in Model Builder for a training session, and it's not working out so well. Model Builder can't seem to find Python. Model Builder Error 000574: Python 2.7 is not installed ...
1
vote
1answer
35 views

Parameters not valid in batch Alter field Arcpy

I am trying to perform a batch of Altering field names through a stand alone python script using Arcpy but it seems I am making a basic mistake in my python code. Using a list I refer to the original ...
4
votes
2answers
925 views

How to change the color of a vector layer in pyQGIS?

Is it possible to change the color while adding a vector layer? I have a vector layer with water areas and want to load it with a blue color style. The layer is added with: ...
-1
votes
0answers
58 views

Creating multiple rows from entries in a single field?

This is probably an easy fix but I have a problem where I would like to use a script or model to export the contents of a feature class's attribute table into a separate table. The separate table ...
2
votes
1answer
126 views

Loading features attributes into numpy array

I am currently working on an add-on that relies on the user choosing fields of tables (that can be from shapefiles, csv or dbf) that are open in QGIS and then making computations with such data. ...
1
vote
2answers
62 views

How to make a REST get request to Geoserver in Python?

I was wondering how to make a GET request to geoserver 2.6 in Python. import urllib2, json url = 'http://localhost:8080/geoserver/rest/workspaces.json' authKey = 'admin:geoserver' headers ...
0
votes
3answers
70 views

How to get the python code for a ArcToolbox call

I was wondering if it is possible to get the exact python code corresponding to a specific run of a tool in ArcToolbox? The reason I ask is that I want to replicate the function in a Tool in Python, ...
1
vote
0answers
17 views

QGIS plugin : connect to polygonselect Maptool signals

Working on a new plugin I am trying to connect a custom function to the mouse right-clic signal while using the polygon selecting maptool. I have activated the maptool through : def run(self): ...
2
votes
1answer
51 views

Python tool won't stop running— Not an infinite loop issue

I'm having trouble with a python script. It is essentially a near duplicate of another script which is running fine-- it is using different files, but apart from that the tools used and the order in ...
0
votes
0answers
27 views

How to read Sentinel-1 tiff images from ESA scihub with Python (and GDAL)?

I can read georeferenced tif images of Sentinel-1 downloaded from PolarView with Python / GDAL but I cannot read the SAFE format from ESA. I need it because PolarView (very nice web interface) keeps ...
0
votes
2answers
44 views

how to make clusters from small polygons with constraint of boundary touching?

I want to make some clusters from a bunch of small polygons with the constraint that each basic unit (polygon) in the resultant cluster must share a portion of its boundary with at least one other ...