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

learn more… | top users | synonyms

3
votes
2answers
64 views

Efficient way to use a CSV to Define a Feature Class

I was provided an Excel Spreadsheet which had the definitions of fields for a feature class. There were several worksheets in the file, each containing field definitions for different feature ...
2
votes
2answers
67 views

Adding layer into ArcMap via Python - .pyt versus .py

In the final step of a script I've built, I'm trying to add a layer into an empty blank map document. This is my block of python code below that adds the layer. Everything works fine if its run from ...
2
votes
2answers
117 views

Is there a function to select a region of a map for further processing with Python?

I would like to select a map's region in QGIS for further processing with python console. Ideally I would like to select the region with a mouse drag (like Gimp crop function) and then e.g. analyze ...
2
votes
2answers
274 views

How to use gdal module in python

I am trying to import gdal modules from osgeo installed in my win64x machine but I get the following error: Traceback (most recent call last): File "", line 1, in from osgeo import ogr ...
2
votes
2answers
155 views

Issue while running grass from subprocess

I have a map based web application using geodjango. I needed to display the crime hotspots within a selected polygon area. I found that v.kernel is very good for generating heatmaps/hotspots. So in my ...
2
votes
2answers
623 views

OSR Python testSR.ImportFromEPSG(4326) returns empty spatial reference

This is empty using Python .ImportFromEPSG(4326): import osr #this fails: testSR = osr.SpatialReference() testSR.ImportFromEPSG(4326) print testSR.ExportToPrettyWkt() #why did the import from EPSG ...
1
vote
2answers
43 views

python crop aerial image

How can i crop an area from an aerial image knowing the latitude and longitude of the image's corners and also the cropped section? I have a LiDAR scan and i and want to crop an image overlay to fit ...
1
vote
2answers
31 views

How to switch to edit mode vector layers, in pyqt application based on qgis?

I do QGIS based application, I need to implement the ability to modify the geometry of vector layers in manual mode (editing nodes), how to do it?
1
vote
2answers
48 views

Getting error after trying to harvest from geoserver into pycsw

I'm using osgeo live to test the pycsw functionality. I need to harvest metadata from an existing Geoserver service (previously configure for this task). I'm quite new with CSW and I'm not sure if my ...
1
vote
2answers
361 views

Using reader Feature Type Name as Writer Feature Type Name in FME

I've search and searched for a way to do this, but no method seems to work. Basically I want to use the input feature class name (feature type in FME) as the output name. This normally works, ...
0
votes
2answers
275 views

problems with py2exe, shapely and fiona

I have a short python script, called test: import fiona with fiona.open("C:\\Essai_compil", 'r') as l_Obs_shp: for shp in l_Obs_shp: print shp['geometry'] I want to compile with ...
0
votes
1answer
71 views

Automating Intersect between two polygon layers

I have a spreadsheet that has county-level occurrences of around 1000 species within the US in the column headings (Species State County). Thus, for some counties there may be multiple species. I ...
0
votes
1answer
39 views

Blank Shapefile output?

When executing the following command string in the active browser, Qgis projects the triangle polygon and associated attribute table correctly. Upon reopening the shapefile in a new browser, the ...
0
votes
1answer
55 views

How to set up pycsw on xampp (win7)?

I'm trying to setting up pycsw on xampp. In the pycsw's default.cfg I've the following configuration: [server] home=C:\Users\masa\xampp-201405\htdocs\pycsw url=http://localhost:8000/pycsw/csw.py ...
0
votes
1answer
244 views

ArcGis Python Script, Calculate Field

Sorry for my English. I am trying to write a Python Script for ArcGis 10.2. #import system modules import arcpy import math from arcpy import env #Set environment options ...
0
votes
1answer
35 views

USGS Earthquake GPS Metadata: Buffer Zones & Python

I'm making a python web crawler that gets USGS earthquake GPS data, and I'm wondering how I should process the latitude longitude pairs to see if it falls within a specified radius of a specified GPS ...
0
votes
1answer
83 views

How to copy attribute values of fields in selected features of one layer to another Qgis programmatically?

I am trying to get the values of attributes of the selected features of a layer field by field and copy to another layer that I created in qgis with pyqgis. Is here a example for this. Every example ...
0
votes
1answer
71 views

Intersecting Polygons

I have 10 shape files containing polygon features. I would like to produce a single polygon layer containing portions of the input polygons where no fewer than two of the aforementioned feature ...
0
votes
1answer
17 views

whereClause syntax for Cost Distance parameter

I am looking to create a sql expression regarding point feature class and selection based on 'Name' field to pass to arcpy.sa.CostDistance() param[0]. similat to: savepath = env.workspace = ...
0
votes
1answer
67 views

Weighted Overlay - passing parameter as text error?

For some reason when hardcoding the save path for the wieghted raster just before runtime it works perfectly. Like such... arcpy.env.cellSize = "90" arcpy.env.workspace = r'G:\Xcel\Route ...
0
votes
1answer
33 views

Changing band for greyscale in plugin

I am trying to make a qgis plugin to change which band in a multiband data set (netCDF where the band represents time steps) is used for a grey scale rendering. The file works fine with manual ...
0
votes
1answer
101 views

visualize bounding box to a polygon

I try to make plugin in QGIS that compare 2 bounding box data, ex: (-32.61442167739459, -32.017601026186576, 13.015663782506067, 13.428145228091406) and (-38.6102467739459, -38.017601026186576, ...
0
votes
1answer
129 views

Polygon Neighbors Analysis script (ArcPy 10.1) producing erroneous results?

I'm struggling with a python script that is supposed to calculate a field for a polygon feature class. The calculation is of the average of a field value (P_TOT_DEV) for each of a given polygon's ...
0
votes
1answer
47 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 ...
0
votes
1answer
190 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
1answer
41 views

How to find zoom levels based on a given bounding box in srs 4326?

I need a function that outputs Zoom Level of a bounding box. for example imagine this bbox: 38.056640499272,22.412109500728,67.060546749272,36.91406262572 How can I found the level?
0
votes
1answer
129 views

Python Geoserver gsconfig create resource from Postgis store table

I've working Workspace and PostGIS store in GeoServer. I would like to use python and gsconfig to create new resource/new layer from existing PostGIS table, with: computing bbox, adding abstract, ...
0
votes
1answer
300 views

Shapely LineString and Polygon intersect?

here is my python code, it does not intersect and when I plot this in google map , I see the intersection. shapely.geometry is in Cartesian plane, how you are solving this? l = "40.65:-114.07 ...
0
votes
1answer
35 views

Access atributte processing algorithm

I have a doubt. In my plugin I use algorithms from processing sextante. For example: Processing.runAlgorithm("grass:v.to.rast.attribute", None, inputLayer, 0, Indexes, extent, cellSize, -1.0, 0.0001, ...
0
votes
1answer
65 views

Plugin QGIS Error “QTextEdit.setPlainText(QString): not enough arguments”

I try to make save as button but i got this error message when try to save file in my directory File "C:\Users\Mr.Pakde/.qgis2/python/plugins\condition\conditiondialog.py", line 80, in saveAs ...
0
votes
1answer
50 views

missing values w/ updatecursor and addfield - gp

I added a feature to my 'overloading' field update script which adds a new properly named field if there is no space left to append the value into the record. This appears to run well, however upon ...
0
votes
1answer
131 views

Manipulating distances in python. Problem with ShapeLength

I have a question more than a problem. I'm getting the distance of a line using python with arcpy.da.SearchCursor or FeatureClassToNumPyArray. The result I get is a double. Example: ...
0
votes
1answer
46 views

Setting up Raster Calculater for QGIS Python Console: ImportError?

I'm trying to follow the instructions here How to evaluate raster calculator expressions from the console? to set up the python console to to use the raster calculator. I tried: from qgis.analysis ...
0
votes
1answer
58 views

Given a BBOX and Zoom Range, Return All TMS Addresses

Given a bounding box using WGS84 lat long pairs (represented any way, not important to me) and a range of zoom levels, I'd like to be able to generate a CSV representing all standard z/x/y TMS ...
0
votes
1answer
400 views

ArcPy SelectLayerByAttribute Error and how to control GP output layer name

I have a rather long script that first references an SDE feature class (non-versioned) to perform an attribute selection based on a SQL selection. I've been unable to figure out why the script will ...
0
votes
1answer
378 views

OSRM API - nearest point fed into routing

I am using the OSRM API to (1) determine the nearest point on a road segment from a give random point with http://server:5000/nearest?loc=_lat_,_lon_ The returned coordinates are used then as an ...
0
votes
1answer
173 views

Orthographic reprojection with mapnik

I am trying to create a google earth like wep application. My problem is that I can't get to render a full view of the earth. The zoom_all() function generates the following error : RuntimeError: ...
0
votes
1answer
485 views

How to clip raster by multiple polygons in multi rasters?

Using QGIS, I want to clip a raster by multiple polygons contained in a single shapefile. I read that gdalwarp can clip my raster by multiple polygons BUT the output is a single merged clipped raster. ...
0
votes
1answer
259 views

How to properly add a Python Module to a Plugin?

I've been trying to use a module that is installed in my local version of python, but whenever I open the plugin itself, it does not manage to locate said module. The module in this case is Noise ...
-1
votes
1answer
55 views

plypython set-up

I am trying to set up plypython on the following system: Windows 8.1 PostgreSQL 9.2 64-bit PostGIS 2.0 Python 2.7.05 I have added plpython2.dll to \PostgreSQL\9.2\lib When I run: CREATE ...
-1
votes
1answer
177 views

Calculating intersection between points and polygon in shapefile

I want accomplish the following using command line python tools only (shapely + fiona) I was wondering if there was a straightforward way for me to do the following: Read in a shapefile (all ...
-1
votes
1answer
188 views

How to Geocode use GeoPy API combined with GDAL Python bindings?

My question is related to the 'How to Geocode 300,000 addresses on the fly?' topic that has been asked by Mehul Nahar. I'm interested in running the processes explained by capooti in response to mehul ...
-1
votes
1answer
80 views

How to find the most east or west longitude around the date line

I've got a python script that is comparing extents of rasters in various projections. Generally it works, but sometimes, especially when the rasters are global in extent or cross the dateline it ...