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

learn more… | top users | synonyms

1
vote
3answers
483 views

Merging many to one data rows in arcpy [duplicate]

The dataset 'left_right_merge' at the bottom of my screen shot shows the dataset that I'm currently working with. The NODE values indicate a single point feature, but is duplicated with many ...
1
vote
2answers
235 views

Conditional Calculate Field Many to One

Anyone have experience calculating a field based on a many to one scenario? I have a feature class in a geodatabase with a field of unique names. Feature Class Table (let's call this "MAINTABLE")...
9
votes
1answer
2k views

Importing GDAL from Python window of ArcGIS for Desktop?

I am trying to run some Python code using the Python window in ArcGIS 10.1 that uses both the arcpy and gdal modules. However, when I try and import the gdal module I get an error: ImportError: No ...
0
votes
3answers
465 views

Finding coordinate pair at any given distance along calibrated route?

Goal: given a route as a set of ordered coordinates, measure the overall distance of the route and be able to access the coordinate pair at any given distance along the route. Application: I'm ...
1
vote
0answers
178 views

How to keep all the unique columns of the shapefiles after merge?

Writing a script where the user is asked to select a directory which contains Folder_1 and Folder_2, both containing shapefiles and run various algorithms. The basic logical process of the script is: ...
3
votes
1answer
320 views

How to implement 1:n joins with string attribute concatenation?

I have a question concerning pyQGIS. I'm a QGIS long-time (self-learned) user, but just a pyQGIS newbie. I want to implement 1:n spatial join with attribute concatenation in the table. I have a vector ...
3
votes
0answers
84 views

How to find the minimum number of inspection points?

I am working on a pipeline inspection project and there is a high liklihood that some inspection tool deployment/retrival station will need to be installed on the pipe network. the problem is ...
0
votes
2answers
906 views

Any clue how to read a .style (ArcGIS) file in Python, for converting into another platform-specific format? [duplicate]

The problem I am face with is that I have all my data as shapefiles but my organization works with only 1 ArcGIS license. Migrating to QGIS to prominent amount of work, I need to convert my ArcGIS ....
1
vote
2answers
767 views

Calculating distances between points two feature datasets

I will do my best to explain my problem. I have two feature datasets. One set has about 100 points that are "city centers." The other dataset has points for each census tract in each of these 100 ...
7
votes
1answer
2k views

Is there a python option to “join attributes by location”?

I am trying to perform the function join attributes by location as found on the QGIS menu Vector>Data Management Tools. I am seeking an open source python option for this. I know arcpy has a spatial ...
1
vote
3answers
228 views

How can I save a column from the attribute table of a vector layer to a variable?

I am trying to write a python plugin for QGIS that can convert three vector layers into an object File for an urban climate model. Unfortunately I am stuck again with a possibly easy question: How ...
-2
votes
1answer
649 views

Creating a python script to import an excel spreadsheet, create shapefile and do some analysis [closed]

I am trying to write a script that does the following. Imports an xls file containing x, y, and other columns with data. I want to make this into a point shape file. I want to take that point ...
2
votes
1answer
216 views

Can not import qgis.core in python

I installed OSGeo4W. I have no other Python version on my computer. I can import osgeo. But when I try to import qgis.core I am getting the error: >>> import qgis.core Traceback (most ...
7
votes
1answer
2k views

Add layer to open mxd using ArcObjects from outside the map document crashes ArcMap

I have an Access database (MS Access 2010) that allows users to enter a variety of land treatment information including associated points, lines, or polygons which users can select and the database ...
0
votes
0answers
191 views

Save map canvas as image with drawings

Is there any way to save the map canvas as an image, with the drawings on map? I draw a line on map, and saved it as image, but only the map saved(without line). I used the following code.. iface....
1
vote
1answer
927 views

Creating a Simple Join Script in Python

I'm trying to create a simple join script in python from an arcgis model, where the only steps are joining a map to a table to produce a map with the table's properties. The attributes joined together ...
7
votes
2answers
3k views

Are there any Python scripts for creating .mxd files?

I am a newbie to both ArcGIS and Python. My requirement is to automate the below MANUAL process: Creating a layer in ArcGIS for Desktop. To put it in another words, creating an ArcMap document (.mxd)...
1
vote
1answer
162 views

assigning census block demographics to postal data (zip/zip+4)

Using Python and R, what would be the best strategy to tie Census Block-level data to zip/zip + 4? I understand we're not matching polygon boundaries (zip codes are not geographies), but somehow ...
2
votes
0answers
421 views

Split a Polygon that Crosses International Date Line, ArcPy [closed]

I have a lots of existing data in WGS-1984 (too much to manually edit). I need to split these lines and polygons by the international date line. More projection information: GEOGCS["GCS_WGS_1984",...
0
votes
1answer
155 views

How to turn WKBtype integer into text for memory layer creation in QGIS

I'm creating a memory layer in QGis python and I need to match the geometry of the memory layer that I am building to the geometry of another layer. After some investigation it seems to me that, as ...
2
votes
2answers
1k views

Is there a way to use gdal functions from OSGeo4W out of their shell?

I would like to use some gdal functions from OSGeo4W (gdal2xyz, and gdal_merge) for my work. I can use them in the OsGeo4W shell provided after installing the package. But I want to use them in ...
1
vote
1answer
341 views

Can i use CreateMapSDDraft in the arcpy.mapping library to set the scale range on tiled services at 10.1?

I know that esri has exposed many properties to the acrpy lib for publishing services. I need to published a tiled/cached map service but i need to set the min/max scale when i do that. Can this be ...
0
votes
1answer
146 views

IOError trying to change PICTURE_ELEMENT sourceImage

I recently ran into this issue after I published a Geoprocessing Task to a production server. It appears to be the same problem discussed in this question, but there is no answer posted there yet. ...
1
vote
0answers
199 views

Generating Custom Map image using Python with small installation footprint? [duplicate]

As a part of reporting requirements there is a need to generate region maps. These highlight various regions (countries or states/provinces) with different colors (or shades of a color) depending on ...
3
votes
1answer
972 views

How to find the id of the last added feature using pyQGIS?

Is it possible to programmatically (in python) access the attribute info of the last added feature in QGIS 2.0 before edits are saved? I'm failing to find a way to get the ID and subsequently other ...
0
votes
1answer
780 views

Creating records in output shapefile using Python

I have a list (output from a previous loop) containing attributes from an input shapefile. What I am stuck with is inputting those attributes to a new shapefile. I have plotted the point data using w....
1
vote
2answers
295 views

How to show attribute table in the form

Is there any way to show the attribute table of a layer in a QtextBrowser or a QtableView on form? I did this but,(shows attribute table in a separate dialog) I have no idea, how to show it in my form....
1
vote
0answers
1k views

How can I get all the latitude, longitude coordinates/points inside a rectangular polygon drawn in google maps?

Google maps lets you obtain a valid address by clicking on houses or buildings one by one, but it would be much easier if you could draw a polygon on the map and get a list of all valid addresses ...
0
votes
0answers
37 views

Copy queried attributes of a shapefile to a new shapefile [duplicate]

I have a polyline-shapefile, which I have successfully read into python using pyshp module. My problem is to save the 2nd node of every line feature into a separate, single point file. I have ...
1
vote
1answer
199 views

Fix Black Border in Python Batch Raster Conversion

I have a Raster image like this as you can see the road are not connected between 2 raster (up and down) its because the scanned raster use different projection (dgn utm 50s) than my data frame (utm ...
3
votes
2answers
1k views

Python Search cursor: Error handling

I have a search cursor that using a SQL query to match addresses and then returns the associated roll number with that address. However, some addresses do not match the roll number and then I believe ...
0
votes
1answer
614 views

Updating records based on value in another field using arcpy.da.UpdateCursor?

I am trying to write file names and paths to two separate fields in a feature class, and the value for each record is dependent on finding a match in a list. When I run my code it writes the same ...
3
votes
1answer
555 views

How to set drawing style for QgsRasterLayer in QGIS 2?

I am using PyQGIS. In QGIS1, I can use: rlayer.setDrawingStyle ( QgsRasterLayer.PalettedColor ) but In QGIS2, I got some errors. I follow pyqgis developer cookbook, put some demo scripts in ...
0
votes
1answer
632 views

How to use python to set a variable to a value from a vector attribute table for GRASS

I'm adding to a python script for GRASS v.transects.py. I want to set a variable to equal a value from a field in the vector's attribute table using a for loop to cycle through each row. This is one ...
0
votes
1answer
226 views

PtDistWithin() with Geometry from query

I'm trying to run two Spatialite queries in Python — the first one picks a single point from a table, and the second tries to find all points that are within a given distance of it. So my queries look ...
1
vote
1answer
211 views

Plugin of Multiple buffers in QGIS

I hope to create a plugin to do multiple buffers in QGIS. This is my script: """/*************************************************************************** izisse_test ...
2
votes
1answer
216 views

Dissolving stream orders using Python rather than VB script?

I've been trying to dissolve streams, using orders programatically and I came across this VB script. The reason is that if I use ArcGIS dissolve, it merges the confluence of the same orders. I would ...
1
vote
2answers
464 views

How to make a routable graph from a Point and a Line shapefile, using Networkx?

I have one shapefile which consists of all the streets in the study area, and another shapefile which has points representing origins and destinations. I would like to know if it's possible to find ...
3
votes
2answers
1k views

GDAL_CALC works but I get a python error at the end of each process that prevents automation

Python 3.4.1, GDAL 1.11.0, OSGEO4W, Windows 7, OSGEO4W Windows Shell All 64-bit (ArcGIS is installed as well 10.2 with Python 2.7) So gdal_translate works fine. gdal_calc provides perfect ...
0
votes
4answers
438 views

Why QgsMapCanvas only appear after lost focus?

I created the simple frame .UI file, compiled with pyuic4, and tried the ShapeViewer.py to load the shapefile "world_borders.shp". However, the map doesn't display at the first time. I have to click ...
2
votes
3answers
130 views

How to analyze and combine overlaping elevation maps into one?

I have 30 overlapping maps that contain tree density heights (measures between 1-100%) of a large area, and I need to get the highest point by pixel and create a new map with all these highest points ...
2
votes
1answer
246 views

Merge from arcpy.ListTables() producing duplicate rows?

I was trying to manually merge file geodatabase tables and was getting Null values for the output file rows and decided to write a python script to merge the files. Here's my script: import arcpy ...
1
vote
0answers
51 views

ExtractByMask with one mask and multiple layers ArcMap

I have a multiple cost distance rasters and multiple feature layers. I need to find and join the feature layers (or associated shapefiles) that lie within a threshold of costs on the cost rasters. For ...
0
votes
1answer
2k views

“EOFError: EOF when reading a line” Using raw_input() in ArcMaps Python Console [duplicate]

I am writing a python script to automate some ArcMaps data processing. I am getting hung up at raw_input(). I am using raw_input to get input from a user. A folder will be created with the folder ...
1
vote
1answer
897 views

Error getting spatial reference as string in arcpy

I am trying to get the string associated with a spatial reference to print into an xml using arcpy and xml.cElementTree. Using the snippet: ... spatial_ref = arcpy.GetParameterAsText(3) ... Prj = et....
1
vote
0answers
3k views

How do I install PyProj module for Python 2.7 on Windows 7?

I am having trouble installing the pyproj module for Python 2.7 on Windows 7. I have tried both the 32 and 64 bit installers from this list but the install crashes each time I run it. I have the ...
1
vote
0answers
123 views

Disable customized context menu for some cases using python

I create a python add-in tool to add new item to the context menu, but I want the new item enabled only when user right click on map when within a buffer of polyline feature. I haven't figure out how ...
1
vote
1answer
200 views

Using os.walk to update broken links and save mxd

I've got code to walk through an entire directory and sub-directories repair broken links in each mxd and save the mxd. It works for the mxd's in the top directory but fails when trying to enter any ...
1
vote
2answers
661 views

Edit Field Value of Feature Class Using Python in ArcGIS 10.2

I have a feature class whose attribute table looks like this: Latitude | Longitude | Area For the area column, all rows are null. I would like to put in a specific value for all rows using Python. ...
1
vote
0answers
179 views

How to export the legend (Layers Panel) as an image in QGIS using the Python console?

I want to export the Legend (the Layers Panel) as an image using the python console.