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

learn more… | top users | synonyms

3
votes
2answers
956 views

finding the index of a field with its name using arcpy

I wan to extract the index of a field with its name using arcpy. it seems what I have written is so dumb. is there any simpler method with out using for-loop. import arcpy fieldList = arcpy....
0
votes
2answers
391 views

creating many shapefiles from one multi-polygonal shapefile

Does anyone know of a way to create a standalone shapefile from every entry (row) in the attribute table of a large set of polygons (4909 in total)? preferably using python in a way that the results ...
3
votes
2answers
477 views

Can I convert the sharp edges of a polygon easily to round edges?

I use the python module "shapely" to make Polygons (like below). Can I convert the sharp edges of a polygon easily to round edges? Such that when I have unified some polygons I could issue a command ...
1
vote
1answer
1k views

Problem loading spatialite extension with sqlite3 Python connector

I am trying to load the Spaltialite extension using SQLite3 python module. I am getting the following error when I run my script. EOFError: [Errno 10054] An existing connection was forcibly closed ...
0
votes
3answers
428 views

Is it possible to publish a map service from an SD File purely over REST (no arcpy)?

I want to automate the publishing of map services, but without using the arcpy library. Is it possible to publish a map service from an .SD File purely over REST (no arcpy)?
2
votes
1answer
4k views

Display a georeferenced DEM surface in 3D matplotlib

I want to use a DEM file to generate a simulated terrain surface using matplotlib. But I do not know how to georeference the raster coordinates to a given CRS. Nor do I know how to express the ...
3
votes
1answer
342 views

Help with “projected XY to screen xy” conversion algorithm?

I'm trying to convert coordinates from any given coordinate system (be they the coordinates of a graph figure or a map, meaning they can span across both negative and positive values across any of the ...
0
votes
0answers
146 views

pyspatialite, cannot update geometry

I'm trying to update geometry column with pyspatialite but it returns error: spatialite.dbapi2.IntegrityError: mypoints.the_geom violates Geometry constraint [geom-type or SRID not allowed] My ...
0
votes
3answers
765 views

Using ArcPy and Python to export Table View as PDF format?

I am using ArcPy to select from a layer and display the result on map as a table view. Now I want to export the table view as a PDF format. Anybody know about it? If it is beyond ArcPy to export ...
1
vote
1answer
1k views

Subquery a list in select layer by attribute management using python

I am new to python and having trouble with my script. I want only the records from the COMTRS field (of the pur_dbf_view) that match values from the ComtrsList to be selected; however, I am getting ...
0
votes
1answer
174 views

gdal_calc.py raster reclass equivalent of ArcGIS reclass or r.reclass in grass

Would like to achieve this in gdal_calc.Py if possible. I have a raster with a value of 0-7 in GTiff I would like one output raster where NoData = 0,1,2,3 1 = 4,5,6,7 So basically an ArcGIS ...
4
votes
1answer
283 views

Can I use a QGIS geometry in a Python PostGIS query?

I have a "Python Init function" for the attribute form in QGIS, which is called when I add a new geometry. I can access the geometry of the newly created feature like this geom = feature.geometry() ...
5
votes
2answers
347 views

How to check if Waste-Water Manhole connect to lines with the right flow direction?

I have a waste-water geometric network with flow direction throw digitize direction. I want to check if each manhole have a line flowing in and a line flowing out. Is there an attribute of line vertex ...
-1
votes
1answer
466 views

Convert Road layer shapefile into a graph

I'm going to use Astar algorithm to find a path between two selected points.And I used the following code for that.Is there any way I could convert my Road Layer shapefile into a graph, which can be ...
0
votes
0answers
327 views

How to delineate Watershed from DEM in python?

My question is that without help of arcpy how can I delineate watersheds from DEM? If anyone help me upto catchment delineation it will be very helpful to me
1
vote
2answers
75 views

Problems searching for attributes

I've come across this problem a few times and have not been able to find a satisfactory answer. I am new to python though so perhaps I have not been looking in the right places! Anyway, In ArcGIS I ...
0
votes
1answer
329 views

How to open QGIS using VBA?

Using VBA, I want to open QGIS (the actual GUI) and execute a simple script through the Python console or the editor.
1
vote
0answers
169 views

QGIS 2.2.0 GRASS not worked

I use QGIS 2.2.0 on Windows 8 64 bit. I try to install use Standalone installer and OSGeo4W and still get the error And then i start GRASS C:\Program Files (x86)\GRASS 6.4.2>grass64 ╬╪╚┴╩└: ...
2
votes
1answer
1k views

How to create a symbol with graduated marker size and line width in pyQGIS?

I am using QGIS + pyQGIS + Python. Layer geometry is line. Each line has a magnitude and a direction (from 1st node to 2nd node). What exactly I want to perform is "graduated width marker symbology ...
0
votes
1answer
125 views

Add attributes to an output shapefile schema (via meta)

I am trying to copy the schema of an existing shapefile and add to it in an output shapefile. Somehow I am not getting anywhere with this...though the docs seem very clear (and the download page : ...
0
votes
2answers
133 views

Generating multiple visibility bearings in ArcGIS as automated procedure?

As a relative newbie to GIS, I'm quite sure there will be a way to do this in GIS (probably most efficiently by writing a python script!) but (to me anyway) this is a fairly complex thing to do and I ...
0
votes
1answer
117 views

Change Attirbute Name (Fiona)

I have a shapefile that has 2 attributes named "X" and "Y". I want to change this to "Lon" and "Lat". My (limited) experience tells me that I should read in the original shapefile, preserve the ...
2
votes
2answers
793 views

How to create a circle vector layer with 12 sectors with Python/pyqgis?

From a given coordinate and radius I need to create a circle with 12 sectors. The sectors should start from north (0 degrees) in 30 degree steps clockwise. I could also imagine to create the sectors ...
0
votes
1answer
72 views

Set the parameter values retreved from database in Arcgis script tool in toolbox using python [closed]

I am now thinking a function to allow user to modify his information stored in the database. To do this, I developed 2 script tools in ArcGIS(10.2), One is the ask user to input his ID, the other one ...
1
vote
2answers
499 views

Open source, Python option for Arcpy's Con tool?

I'm wondering if there is a method to perform a function similar to Arc's Con tool, using open source options in python. I have two rasters, one that shows locations of clouds and the other is just a ...
3
votes
1answer
487 views

How to make a Python script that does a task in QGIS?

When I run this code in QGIS Python console it works fine, but I want want to run it from an external Python script. What I want the script to do is: 1- Open QGIS 2- Run this code 3- Close QGIS from ...
3
votes
1answer
2k views

Converting Matplotlib contour objects to Shapely objects

I would like to convert a Matplotlib contour collection to a Shapely geometry collection, in order to calculate the area contained within closed contour lines. it seems that all the information is ...
3
votes
2answers
8k views

Using conditional statements in Field Calculator

I have two fields, "P_Area" and "PC_Area", whose sizes I'd like to compare in a third field, "Areas_Equal". The "Areas_Equal" field should contain a "Yes" if the other two are equal and a "No" if they ...
1
vote
1answer
1k views

How to add the attributes of selected features to my new layer in pyqgis

I am trying to create a layer from the selected points coordinates. I am drawing a circle to the coordinate of selected points. I also need to add initial attributes of my selected features to my new ...
1
vote
1answer
90 views

Problems running grass.mapcalc

I am using the GRASS Python Scripting Library on Ubuntu 14.04 with the packages available from the Universe repository. I have the following block of code: try: grass.mapcalc("$...
0
votes
1answer
66 views

QGIS Python - Executing functions for multiple directories

I would like the user to select 2 directories, each directory having their own defined algorithms. For example: First_Directory will add a column to all shapefiles in the chosen directory; ...
4
votes
1answer
1k views

How to find the distance to the nearest coastline from a land-based latitude/longitude point?

I'm completely new to GIS and unsure what I need to solve this problem: In the image above turquoise is actually land. Sea is white. I tried using the NASA Distance to the Nearest Coast data dump ...
1
vote
1answer
301 views

how to assign keyboard key to python script in arcmap

I have a python script and I use it through my toolbox but I use it many times in a day. Is it possible to assign a keyboard shortcut to my tool ?
3
votes
1answer
1k views

QGIS says “osgeo [python-gdal]” isn't installed on Mac

I know similar questions have been asked (like this one) -- but I'm encountering a unique problem. I installed GDAL complete 1.11 from KyngChaos, but when I open QGIS, I get the following error ...
2
votes
1answer
997 views

How to open OSM file in Python with GDAL

I am trying to open a *.osm file in Python with GDAL. But I doesn't work. What am I doing wrong? import ogr driver = ogr.GetDriverByName('OSM') datasource = driver.Open('map.osm') layer = datasource....
1
vote
1answer
2k views

How to Properly Loop Through a List With a Cursor?

My team is working on a python script that does various updating tasks to our streets file. A part we are getting stuck on is trying to update a field. Currently we have a field FEDROUTE that is ...
1
vote
1answer
357 views

Is there a way to access the GeoProcessing tools in QGIS, in code

Is there a way to access the GeoProcessing tools in QGIS, in code? I want to find the Difference between two layers and make a new shapefile with the result.It can be manually done by selecting vector-...
3
votes
2answers
204 views

FME - Conditionally make parameters visible

I have an FME workspace in which the user sets several parameters before running. The first parameter is a multiple choice where the user selects which layers they would like to process. Depending on ...
0
votes
2answers
122 views

FME - Filter Parameter choices to only projected coordinate systems

I have an FME workspace in which the user sets several parameters before running. One of the parameters is a Coordinate System. I would like the user to choose only Projected Coordinate Systems as ...
2
votes
2answers
80 views

GetCount raises exception

I opted to use gp.GetCount in my script. This raises an exception, though the process completes "successfully" while returning zero features (which is untrue). I'm wondering how this could be. Using ...
1
vote
1answer
115 views

ArcGIS Geoprocessing “None Type” Second time

I am trying to publish an ArcGIS geoprocessing service in 10.2 I have a script setup that will run multiple times through Pycharm, but when I run it in ArcGIS it will run once and return this error: ...
3
votes
1answer
2k views

Sorting shapefile records on attribute to update new rank field using Python?

Basically, I am trying to sort the records within the shapefile based on an attribute(eg: pop) and then add a new field that indicates ranking based on the sort. eg: the record with the highest pop ...
1
vote
0answers
62 views

Is is possible to run the QGIS python module server side?

I'm pretty new to programming but I want to use the QGIS python module on a server. Basically I am trying to create a JPEG map programatically on a server. I already have a python script to create ...
1
vote
1answer
233 views

Locating 0-points multipart polygons using arcpy

By Running this piece of code: import arcpy arcpy.env.workspace = "C:\working\GEODATABASES\\tests.gdb" fcList = arcpy.ListFeatureClasses() for fc in fcList: with arcpy.da.SearchCursor(fc, ("OID@",...
0
votes
1answer
150 views

PyQGIS: How to send signal to call unsetMapTool()

I've created a map tool to add some functionality to adding point features in QGIS 2.0. I'd like to unset this map tool if editing has manually been terminated by user. I've attempted as shown below, ...
2
votes
1answer
449 views

QGIS Python Scripting - Syntax definition

After playing around with the Model Builder and looking at the python scripts of what I make, I noticed some syntax which I can't find the definition anywhere. Here is a snippet of the python script ...
5
votes
2answers
466 views

Accessing arcpy functions from IronPython

I made a few python scripts to automatically optimize and publish my data.. Now i'm thinking about writing a .NET program to manage and run these scripts: the obvious way is invoking them in a ...
4
votes
2answers
6k views

If/then VBScript/Python code equivalent to SQL IN ('x','y',z') expression in ArcGIS Field Calculator?

we need just like the following, only with multiple records selected... have tried the "xyz" in (1,2,3,4,5) SQL code but of course VB Script wants something else. dim f2 if [field1] = "foo" then f2 ...
2
votes
1answer
2k views

import qgis.core problem with QgsFeature and QgsGeometry

I've tried several previous answers to similar questions - no luck. I run Win7. I installed QGIS 2.2.0 using the OSGeo4W 32-bit installer. Problem 1) from qgis.core import * failed. Solution 1) ...
2
votes
2answers
952 views

Calculate Slope Between Pixels Through Time: A Time Series Trend Analysis

I have a 9 year weekly time series (~500 raster grids of equal cell size and extent). I'm interested in obtaining the regression line slope between pixels (Imagine stacking all 500 grids on top of ...