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

learn more… | top users | synonyms

1
vote
1answer
480 views

Removing Shadows from Aerial Imagery in ArcGIS

I am creating an automatic change detection tool in Python for ArcGIS 10.2 using aerial imagery. The tool works but it also detects shadow change. I have been able to create a classified raster of ...
1
vote
0answers
39 views

Automating Raster Calculator [duplicate]

I am seeing mixed answers on other posts so I thought I would ask here. Is there a way to automate the raster calculator? I am trying to create a series of inundation maps at fine increments. It ...
1
vote
0answers
45 views

QGIS / Python differences in result

I am working on a script that automatically calculates something for me. To be sure, I checked the result of the calculation with QGIS. Here are the statistics: STATISTICS PYTHON-CODE: MAXIMUM = 0....
1
vote
2answers
335 views

Accessing Rest Api of Geoserver using python

is there anyway to access Rest Api of Geoserver using python other than gsconfig. since it doesn't support for me.i'm using python3.4 Like creating workspace,adding layers programmatically to ...
1
vote
1answer
206 views

How to use subprocess module to run QGIS toolbox modules

I've been trying to use subprocess for a while now, but can't get it to work when calling QGIS modules. Here's my script: import subprocess call("processing.runalg('gdalogr:merge', [rasterin1,...
0
votes
1answer
74 views

How to use wildcards for characters inside the name [closed]

I have shapefiles in that form RCs4s3s2c10_S2C10 RCs4s3s2c20_S2C20 RCs4s3s2c30_S2C30 RCs4s3s2c40_S2C40 RCs4s3s2c50_S2C50 RCs4s3s2c60_S2C60 RCs4s3s2cb_S2CB RCs4s3s2c10_1_S2C10 RCs4s3s2c20_1_S2C20 ...
2
votes
1answer
85 views

Replace Expression in Field Calculator to include a Numeric Range

There has to be an easier way to clean up address data, I hope? GeoCoding, I need to remove APT# and any number or letter that follows.
2
votes
1answer
65 views

Determining the cosine of the angles between two lines using dot product

SCENARIO:I have a certain buffered region (lets say a one way road segment), which will have only one allowed direction represented by vector A. And, there are several other lines which may or may not ...
0
votes
1answer
65 views

Error 000177: Failed to create relationship class in Python

I am trying to use python to create a relationship class between a feature class and dbf. I keep getting, "ERROR 000177: Relationship classes can only be created in the same database that ...
1
vote
0answers
42 views

Error in installing gsconfig python library for manipulating a GeoServer instance via the GeoServer RESTConfig API

I like to install gsconfig python libraries I got error like this Collecting gsconfig Downloading gsconfig-0.6.13.tar.gz Traceback (most recent call last): File "<string>", line 20, in <...
1
vote
0answers
142 views

Use Gdal from Python IDE

I am just starting to experiment with writing python scripts. I have been using QGIS and the Osgeo4w shell for sometime now. I would like to be able to call Gdal from the python IDE i have installed (...
2
votes
1answer
183 views

problems overwriting files in ArcSDE with python script

I'm having a great deal of difficulty overwriting files when using ArcSDE. The usual 'arcpy.env.overwriteOutput' method doesn't seem to work. When I run this script: import arcpy from arcpy import ...
0
votes
1answer
68 views

Can you force cartoDB to read a file as a synced table

I'm writing a script that downloads from a data source every hour. I then do some processing on it and convert it to a kml saved in my google drive. As the data is updating every hour I want to be ...
1
vote
1answer
156 views

Field Calculation Python acceleration of gravity [closed]

In my problem i have a column DISTANCE and a column TIME and i want to calculate with python in field calculator how much time it needs to make that distance (for example 1000m) if you are moving with ...
1
vote
0answers
336 views

Update Data Driven Pages Index layer via Python

I am trying to find a way to update the index layer for data driven pages through arcpy. I understand that to access DDP from arcpy.mapping you need to be referencing a MXD that has DDP enabled. My ...
2
votes
2answers
165 views

Creating a really large shapefile without eating all the virtual memory

I am doing something similar to below to create a really large (large than memory on my machine) shapefile: shapes = shapefile.Writer(shapefile.POINT) shapes.autoBalance = 1 shapes.field('MyId','N') ...
1
vote
1answer
401 views

How to update raster color ramp with new min/max values using PyQGIS?

I'm trying to update a single band psuedo color renderer with new min/max values from the python console. All the code samples i have been able to find deal with constructing a renderer from scratch, ...
1
vote
0answers
84 views

Checking if Terrain Model has pyramids or data sources using ArcPy/Python?

I asked this question in the ESRI forums and got nothing but a lesson on LAS point cloud classifications — not even remotely what I had asked. I’m hoping someone here can point me in the right ...
-2
votes
1answer
127 views

NumPy Array to Feature Writes only first array object to feature class [closed]

I am consuming a json web-service with NumPy Array to Feature, I have many nested try/except blocks to consume the data/find errors/etc. The script is functional, however I am having trouble seeing ...
1
vote
1answer
43 views

How do I add packages to Standalone Qgis's Python?

I have installed the QGIS 2.8 on my Windows 8.1 system, using the Standalone Installer. I now wish to install the VectorBender plugin. The docs for this plugin say, that I need to install matplotlib ...
3
votes
1answer
247 views

Is there a quick way to see the fields in a table via arcpy?

I frequently find myself in Python debug mode, trying to examine the fields of a layer or table. arcpy.Describe(lyr).fields returns an array of <geoprocessing describe field object> which isn't ...
0
votes
0answers
106 views

How to reshape a polygon with a hole in qgis automatically?

I am new to qgis and appreciate any help/advice. I have a polygon shp file and some of my polygons are within another polygon (again:they are on the same layer). So I want to split the big polygon ...
1
vote
0answers
87 views

I am building QGIS 2.6.1 from source code on windows 7, 32 bit and after successfully building all projects I get this error

Couldn't load PyQGIS. Python support will be disabled. Traceback (most recent call last): File "", line 1, in ImportError: No module named qgis.core Python version: 2.7.4 (default, Apr ...
3
votes
0answers
286 views

Python vs R for analzing spatial data [closed]

I am relatively new to both the python and R programming languages, having a small amount of exposure to both. Which language has the edge when working with spatial data and when would you use R over ...
0
votes
1answer
355 views

Vector to Raster using GDAL (gdal.RasterizeLayer) : Error in Output Raster(generating output raster with NAN values)

For rasterizing a vector layer, I have tried gdal.RasterizeLayer method for converting vector(shapefile) to raster(tiff). But its giving raster output with NAN values(complete black image). I need ...
0
votes
0answers
200 views

Can I use Redfearn's formula to convert between WGS84 and UTM?

I have GPS data at 1 point per second for several weeks in WGS84 latitude/longitude form. It was collected using a UBLOX GPS in 2015 in Australia (Zone 55s) Would it be correct to use Redfearn's ...
2
votes
0answers
82 views

How to clip or erase a csv file with a non-rectangular shapefile?

I can clip a csv file with a ogr2ogr having 4 extents of a rectangular polygon. Is it possible to clip a csv file with a non rectangular shapefile?
2
votes
2answers
326 views

How to make a movie from a ascii grid files?

Lets say I have several ascii grid files for each time step for a particular location. I would like to spit out the image for each of these ascii grid files and legend on the side bar. I also would ...
-1
votes
1answer
222 views

How to use a wildcard for a field to sort many shp

I want to use a wildcard for a field so as to sort a group of shapefiles. All of them have the same columns but there is one column in each of them with only the first 4 characters common. Is ...
1
vote
1answer
138 views

How do you use Windows folders as acrpy Parameters

I am having problems trying to perform simple folder management operations in an ArcGIS python toolbox. First here's some code snippets (I have left out all of the stuff that does not pertain to this ...
0
votes
0answers
129 views

google geocoding api usage restrictions and other free api options

I would like to use the google geocoding api to geocode a few thousand addresses. I understand the daily limit for a free account is 2,500. I want to retrieve the coordinates and then map them in ...
2
votes
1answer
295 views

Identifying Left and Right Mouse Button Clicks in PyQgis application

I am relatively new to PyQt. I am trying create a custom plugin for Qgis which enables the user to select some features by drawing polygon on the canvas using mouse clicks and then performs ...
1
vote
1answer
61 views

Can I mimic the “open” action in QGIS with python and manipulate the field content used for the action?

I am using QGIS actions (type: "open") to link ECM searches to objects on the map. The basic syntax is very simple: protocol://searchterm1=[%"field1"%]&searchterm2=[%"field2"%] However, some ...
1
vote
0answers
127 views

Quick way to select by attribute and zoom to selection in QGIS

I need to visually control hundreds of points in a layer on top of aerial photographs. Using the attribute table, selecting (with column filter expression) and zooming manually is cumbersome - so I ...
0
votes
1answer
305 views

Batch Mosaic Rasters: arcpy.MosaicToNewRaster_management (ERROR 000628: Cannot set input into parameter coordinate_system_for_the_raster)

I am trying to batch mosaic rasters using arcpy.MosaicToNewRaster_management in arcpy: #This is for mosaicking a lot of rasters import arcpy from arcpy import env arcpy....
1
vote
0answers
75 views

Importing tif image from PostgreSQL Database to webpage

I am trying to import tif images from PostgreSQL Database to my website. In the website, I need to make a button control which will open the images. I think I have to use Python as a server side ...
5
votes
2answers
707 views

Advanced python labeling in ArcMap - working with null values

I'm trying to label certain features using a specific field (osm_name_58_en). However, when this field is null, I want to label the feature with (som_english_32_name). If both are null then I want ...
1
vote
1answer
124 views

Auto-update a feature class with Python + .csv?

I am hoping to find out how to automatically update a feature class. We have semi-regular updates to an ArcGIS Viewer for Flex website. Looks like a Python script may be the way to go about it (and I ...
3
votes
1answer
669 views

Anything Python related causes crash to ArcGIS for Desktop?

Problem: ArcMap 10.3 crashes to desktop anytime I use a Python related function (Toolboxes, Field Calulator, Etc) How I got to this point: I was trying to gain some familiarity with Python ...
0
votes
0answers
113 views

How to create a Shapefile of a subset of input features in PyQGIS?

I'm having some trouble applying my Python skills to develop some scripts in PyQGIS. For me it's far simpler to write it in ArcPy, but then again I wouldn't get that much satisfaction after finishing ...
1
vote
1answer
252 views

QGIS Python Scripting permission denied

I've just started a Python Scripting tutorial on QGIS. The problem arrives right at the end, when I have to use this script: output_file = open('c:/Users/Ujaval/Desktop/airports.txt', 'w') for f in ...
3
votes
4answers
219 views

Removing first and third zeros of a fieldname in label expression

The fieldname only has three characters. I dont want to display zeros in the first and third values. example below. 010 = 1 100 = 10 01Y = 1Y 10Y = 10Y I have figured out how to strip out zeros ...
2
votes
0answers
127 views

GDAL equivalent to ArcGIS Union?

I have a single shapefile full of overlapping polygons. The ArcGIS Union tool will generate individual, non-overlapping polygons with that geometry. Is there a way to do this with GDAL/OGR or another ...
2
votes
2answers
200 views

Cutting a polygon out of another polygon with python [duplicate]

Licence: ArcEditor Version: ArcMap 10.0 i have some issues to separat my data with python. To illustrate the problem, i have simplifyed the data with figures. I want to cut the blue polygon out ...
1
vote
0answers
50 views

Cloud and cloud shadow identification with MERIS

I am looking at MERIS images to use for marine applications (Chl-a. SPM) and would like to confidently remove cloud and cloud shadow pixels from the data. I have a few ideas: Using the blue band ...
0
votes
1answer
237 views

Can't get the correct version of GDAL installed to use ReadAsArray without getting error in Python

Using the instructions here . . . https://pythongisandstuff.wordpress.com/2011/07/07/installing-gdal-and-ogr-for-python-on-windows/ I downloaded the following files: GDAL-1.11.1.win-amd64-py2.7.msi ...
3
votes
2answers
256 views

Python labeling expression using numeric field not working

I'm trying to label parcels from a shapefile in the Python expression below in ArcGIS 10.2. All parcels should be labeled with the last 8 characters of their [PARCEL_NUM], and parcels 10 acres or ...
0
votes
1answer
21 views

vector clip based on chronological value

I would like to know if there is any algorithm in python that can clip vector features according to the "date" value? To be more specific the newest date is the dominant value and clipping will start ...
1
vote
1answer
66 views

Multiply two lists of raster from two subfolder

I have two subfolders, which contains rasters for each country but which express different variables. I am trying to figure out a way to Loop a multiplication of each raster into a folder with its ...
0
votes
1answer
71 views

Convert DOY to mm-dd format [closed]

I have a problem with my date format. how can I convert day of year values to mm-dd format? For instance if my doy have a value of 161 this will be converted to 06-10 (June 10 for instance). And this ...