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

learn more… | top users | synonyms

4
votes
2answers
1k views

Method for finding circle center points that cover an extent

If I have a geographic extent, is there any way to find all of the center points of all of the 5km circles that would be needed to cover the entire extent? I'm working with an API that will return ...
7
votes
2answers
239 views

How to identify a FRONT LAND?

I work at a town hall. I need to identify the front of the land. Here in my country we call it of "TESTADA" or "FRENTE DE LOTE". I did not find the exact term in English. THE FRONT OF A LAND is a ...
6
votes
1answer
230 views

Why do ECWs go missing when exporting MXD to PDF with ArcPy?

I'm trying to export some MXD maps to PDF with arcpy (ArcGIS 10.1). Everything works fine and dandy except for ECW images, which don't appear in the PDF. However, when I execute the same code from the ...
1
vote
0answers
186 views

Create custom “Export To” context menu tool using selected feature in ArcCatalog

I'm trying to create a custom "Export To" script similar to the Right-Click Export Content menu item "To Geodatabase... (multiple)". I have a 3 feature classes in a file geodatabase. I selected two of ...
2
votes
2answers
788 views

How to create equidistant points in QGIS along a polygon?

I am looking for a solution to distribute with a python script equidistant points in QGIS along a polygon. I tried to use the following code example: ...
2
votes
1answer
565 views

Automatic calculation in one attribute table based on values in another attribute table using python

I am trying to fill the fields of one gdb table based on some statistics over the attributes of another shape file. As it was not possible in ARCgis directly, I am trying to use python coe for this ...
-1
votes
2answers
250 views

How to get layer extension info via Arcpy

I want to get file format of a layer with using Arcpy. I am writing a tool that needs to take raster or feature extension like SHP, TIFF etc. Which functions should i use ? Note: I use ArcGIS 10
1
vote
1answer
195 views

trouble deleting data in an fgdb using python

I am trying to use delete_management to delete files in a gdb based on whether or not they are in a list. I keep getting an error with the delete_management tool that says the file either does not ...
1
vote
2answers
148 views

How to transfer multi attributes in to features?

i have 2 feature classes, i want to transfer all attributes by one click to another attribute feature class, i have used Spatial Adjustment tools, but the challenge with me is to, how to transfer all ...
2
votes
1answer
1k views

How to make autoIncrement number start with (001…009, 011…099)?

i have a challenge with my time in my work, i have feature class called "Parcel" storing many data, i generated an autoIncrement number by using field calculator, in field "parcel code" Expression: ...
0
votes
1answer
908 views

Where Clause problems when all parts are user input variables

I have written some generic user ipnut code for a toolbox script (ArcGIS 10.0) but because everything is a variable I don't know how to write the SQL whereclause. Please can someone help and suggest ...
1
vote
1answer
99 views

python scripts for accessing info from a postgresql database in case of a qgis application?

I am beginner and I need to know if I could access a postgresql database from python script made for a qgis application
0
votes
2answers
320 views

extract edge of a line of a vector layer in a new point layer [closed]

I'm using python to build a plug-in in Qgis. I need to extract the edges of each line feature of a vector layer to create a new point layer of all edges. Someone can help me? My data input is a vector ...
0
votes
1answer
79 views

How to Add Field Between Two Existing Field in Shapefile

Can you please let me know if it is possible to add a new Field between two existing field using ArcPy? For Example if I have "Name", "Area", fiels already in my shapefile can I add a new field called ...
3
votes
2answers
799 views

Create point geometry in a feature class with all the attributes prepopulated

So here is the problem. I have a text file with space delimited values. A standard text file with the headings separated by spaces and then on each next line values for these headings (columns) also ...
3
votes
1answer
329 views

How to apply a color gradient to a shapefile from QGIS console?

Anyone knows how to do a color gradient from the console to a shape file in Qgis (orin Python)? Like the image below.
3
votes
3answers
1k views

Formula for coordinates (lat, lon), azimuth and distance.

I would like to know how to find a set of coordinates on a small circle some distance from another point. For example, lets say A = (39.73, -104.98) #Denver B = (39.83, -106.06) #point ...
3
votes
3answers
286 views

How to convert Midpoint of Polyline into Point Feature using ArcPy?

I am trying to write a python script to determine flow direction on a water pipe network. I've used this script to find the mid point along each polyline, now i need to convert this data into a point ...
0
votes
1answer
487 views

Better way to handle HDF and GeoTIFF files

I am using GDAL and MODIS reprojection tool in my python project and want to know better way for math new indices of HDF data (MODIS data). Also I want to visualize them in web. Which format is more ...
0
votes
1answer
359 views

Python: Using re and update cursor over three columns [closed]

I'm working on this code and got a little to excited...only to find out it wasn't perfect! I'm trying to eliminate the suffix (row1) of roadnames(row0) and parse just the roadname to a new ...
8
votes
4answers
5k views

How to convert geometry to WKT using ArcPy?

I would like to be able to convert a single feature (or potentially many features) in an ArcGIS layer into Well Known Text (WKT) format during a geoprocessing operation using Python and ArcPy (no ...
1
vote
1answer
266 views

GDAL reports no bands for a 12-band raster and no projection either

When opened in ArcGIS, the raster correctly unfolds 12 bands, but when I try to open it with GDAL through Python, it appears as a raster that has zero bands. Here is the code: import gdal from ...
3
votes
3answers
3k views

How to batch clip to single-part shapefile using ArcGIS 9.3 and Python 2.5?

How do I clip a shapefile to a single-part shapefile using ArcGIS 9.3.1 and Python 2.5? Example: I have a set of lines (lines.shp) I am trying to clip to each record in a shapefile called Grids.shp ...
7
votes
1answer
477 views

Draw middle ordinate perpendicular to chord, out to curve, with arcpy

I've been assigned to find the radius of road curves. My boss found this article describing a pretty straightforward method of doing so: find the start/end points of the curve; create a chord; ...
7
votes
3answers
2k views

How to Receive an Email after a script has run?

I'm working my way through my first attempt at scheduled tasks with modelbuilder and python. So far I have: Created a model that reconciles our SDE versions and then compresses the database. ...
11
votes
3answers
1k views

How to convert LAS file to numpy array?

I have begun learning how to manipulate LAS data in python and wanted to see how others handle LAS files. I would like to read the points (I am using a numpy array), and filter out classes 1 and 2 ...
6
votes
2answers
647 views

How do I get the projection of an LAS file using liblas?

I am using the python api of liblas to read/write las files. I would like to get the spatial reference information of the las file using liblas, but have been unable to figure it out. Can anyone ...
0
votes
0answers
80 views

Loading QGIS and files from Python's use of os.system?

I am trying to open QGIS and load two files using the following code in Python: os.system('qgis D:\data\vector.shp D:\data\raster.tif') However it returns this error message: 'qgis' is not ...
1
vote
1answer
209 views

Polygons not appearing when creating shapefile based on criteria using GDAL and Python (area calculated properly)

I am looping over the features of a shapefile containing multiple polygons. If certain criteria are met, I want to save the current feature in a new shapefile. I am doing this using GDAL 1.10.0 with ...
1
vote
1answer
339 views

Flow Direction Arrows in python

Any help would be appreciated. I am trying to write a python script to determine flow direction on a water pipe network. I have calculated the mid point of each section of pipe and created a point ...
1
vote
1answer
754 views

using GDAL/python to stack georeferenced images of different sizes

I am in the process of porting a code I wrote in IDL (interactive data language) to python but am running into a bit of a problem that I am hoping someone can help me with. The code goes like this: ...
5
votes
1answer
741 views

Create point between two points with given distance or length

What I have: Coordinates of Point A and Point C and therefore vector and azimuth from A to C What I want: Coordinates of Point B (between A and C) with a distance from A to B of 35.2 meters. ...
0
votes
1answer
1k views

Finding differences between two shapefiles in Python

I'm trying to find out whether there were any border changes between two yearly shapefiles. I've heard a lot of praise for Python's open source GIS tools like Shapely/Fiona, and thought I would try ...
0
votes
1answer
3k views

How to set PYTHONPATH in OS X to use PyQGIS?

I am trying to use PyQGIS in custom applications. I am on Mac OS X 10.9 with QGIS 2.0.1-Dufour. I know I have to set the the PYTHONPATH environment variable. This is the path to QGIS on OS X: ...
1
vote
1answer
520 views

Convert measurement geometry from ground reference frame to satellite reference frame

I need to convert solar_zenith_angle, solar_azimuth_angle, platform_zenith_angle, and platform_azimuth_angle from a reference frame at the WGS84 surface to a reference frame at the satellite observer, ...
0
votes
0answers
64 views

Why legendLayerActions is not available in QgsLegendInterface

There's a reason why QgsLegendInterface (the python interface to QgsLegend) doesn't expose: QList< LegendLayerAction > QgsLegend::legendLayerActions( QgsMapLayer::LayerType type ) const method? ...
1
vote
0answers
205 views

How to run PostgreSQL stored procedures from QGIS? [closed]

I'm trying to figure out if it's possible to run a stored procedure in a Postgres database from QGIS, more specifically I would like to do this from within a python plugin. In this case I'm trying to ...
1
vote
1answer
1k views

Insert Cursor fails to write results to table

Insert Cursor fails to write results to table. The Search Cursor is returning the results, but errors when trying to insert into field. # Insert Cursor for GDB table to GDB table # Import modules, ...
0
votes
0answers
89 views

Populate a new field with list of the top 5 values of all fields for that record

I have 15 fields. I would like to add a new field that is populated with the field name and value of the top 5 highest values of the 15 fields for each record in the table. The new field would be a ...
1
vote
1answer
259 views

Unsupported operand error writing to file (Python, arcgis)

In the code below, I get a "TypeError: unsupported operand type(s) for +: 'NoneType and 'str'" on the f.write line. If I don't concatenate, it runs fine, but all the layer names run together. I ...
6
votes
3answers
699 views

How to Run Viewshed on one feature at time in shapefile?

I have over 100 features in a shapefile and I want to create individual viewsheds for each and then run processes on each of those resulting viewsheds in ArcGIS using Python. Is it possible to write ...
6
votes
1answer
2k views

Why close a dataset in GDAL Python?

I see quite often in Python GDAL code that people close datasets at the end of their script. Why does it makes sense to close a dataset in Python GDAL? Are there any consequences if I don't do it? ...
5
votes
2answers
2k views

Are there Python packages for Landsat 8 processing?

While there are a number of python projects for MODIS acquisition and processing (pyModis and get_modis), I have yet to encounter a toolset for automating Landsat 8 processing in Python. There's ...
1
vote
1answer
473 views

Map Layout Elements In New Layer Automated Using Python

I was wondering if there was a way to put a title and integer data from an attribute table in a new layer created in Python. I was able to create a new layer using copy features, zoom to the selected ...
2
votes
1answer
475 views

Getting SAGA-GIS Python Bindings without compiling from source

I am very new to SAGA-GIS, but I recently developed a workflow using SAGA for processing LiDAR data. Now I need to batch the process using Python. I installed SAGA from the Windows installer and I ...
0
votes
0answers
675 views

RasterXSize from gdal Nonetype error

I have this code to create a raster file in python: fileInf = QFileInfo(userIntervals) baseNam = fileInf.baseName() rlay = QgsRasterLayer(userIntervals, baseNam) gdalRast = ...
4
votes
1answer
1k views

Select by attribute expression invalid in python

I'm not sure how my expression is failing. Any ideas? inPoints = "points.shp" i = 0 arcpy.MakeFeatureLayer_management(inPoints, "pts") arcpy.SelectLayerByAttribute_management("pts", "NEW_SELECTION", ...
1
vote
1answer
805 views

Why does ListLayoutElements give IndexError: list index out of range?

When trying to create a legend on ArcMap use the following Python/ArcPy code I received an error message: import arcpy mxd = arcpy.mapping.MapDocument(r"E:\cities\MyProject.mxd") df = ...
9
votes
6answers
8k views

ArcGIS/ArcPy Python, how to get latitude and longitude of a projected point?

I have a point feature in a feature class that is being accessed by ArcPy. The point is projected but I need to find an efficient means to get the unprojected latitude and longitude for this point. ...
1
vote
2answers
242 views

Difference between OGRGeometry and Geometry in C++ and python

According to the documentation, there is a python version of OGRGeometry's C++ class (http://www.gdal.org/ogr/classOGRGeometry.html) which is http://gdal.org/python/osgeo.ogr.Geometry-class.html ...