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

learn more… | top users | synonyms

5
votes
1answer
490 views

Select and copy attributes features in ArcMap using Python add-in tool?

I'm trying to add a tool to an add-in toolbar to select features from an existing feature class and copy attributes to another feature class,from point to point. TopoHeight_Field_Name='Comment' def ...
3
votes
1answer
487 views

How do I smooth edges of polygons in QGIS?

Is there a built in function or plugin I can use to achieve smooth lines for my polygons instead of the blocky ones I have now?
2
votes
1answer
22 views

AttributeError: 'SpatialReference' object has no attribute 'ImportFromESPG'

I am trying to create a spatial reference using ImportFromESPG and this is what is returned: AttributeError: 'SpatialReference' object has not attribute 'ImportFromESPG' Code and error are copied ...
4
votes
2answers
130 views

How to update raster pixel values based on another raster with conditions in QGIS?

I have 2 land cover rasters. I want to update parts of the first raster with values from the second using QGIS. The workflow should be like this: 1. Select values equal to 30 and 34 in the first ...
1
vote
0answers
32 views

Find and filter points within a variable radius

I have a set of projected spatial points (x,y,id,dist,value) that needs to be filtered. For each point I have to find the other ones falling within a variable distance. For each selection, I need only ...
3
votes
1answer
96 views

How to convert multiple csv files to shp using python and no arcpy?

At my internship, I have been asked to learn Python on the job. Its been a rough road never having coded before so please forgive me if my problem is elementary. The Lightning CSV files I am ...
1
vote
1answer
48 views

Why are there my_form_open() Python functions in my QGIS project file?

I have a project open with 6 layers (MapInfo .tab format), have set the style for each layer and added labels for 1 of them. So for a simple project the .qgs file is about 2000 lines long and among ...
1
vote
2answers
7k views

KeyboardInterrupt doesn't work when arcgisscripting module is loaded

The KeyboardInterrupt exception -- thrown when Ctrl-C, or Ctrl-Break on some machines, are pressed in a Python console window -- is commonly used to gracefully break out of long-process loops, ...
1
vote
1answer
24 views

Updating a raster using a shapefile [closed]

I have a huge raster file and want to increment all the values within a relatively small shapefile by 1 in Python. Speed is of critical importance, as I have to do this many times. I have been ...
1
vote
1answer
78 views

How do I fill a Field Combo Box with vector layer fields in QGIS?

I am new to python and I am trying to write my first plugin for QGIS. Earlier I asked this question to find out how I could populate a combobox with point vector layers (and nothing else). I've been ...
1
vote
2answers
49 views

Rpy2: how do I use the function 'as' (r library methods) in python?

I want to coerce an object of class Raster Layer to class SpatialPixelsDataFrame since that is required as input for the function autoKrige from the r library automap. However, when I try to use as ...
1
vote
2answers
56 views

Renaming shapefiles and raster files with Python

I'm trying to add "_Kiwi" to the name of multiple shapefiles and raster files in a folder C:\data. For the files with multiple extensions (.shp.xml for example), I can't get the files to rename ...
4
votes
1answer
511 views

How do I link OGR with spatialite?

Problem & Evidence: I'm trying to create a spatialite file with the code snippet below. However the outDataSource object gets assigned NoneType and all subsequent (unshown) code fails. I'm ...
2
votes
3answers
636 views

How to add an image to an existing data frame using Python?

This is the first time I am using Python on my own, and I have searched for this answer for a few hours with no luck. I am wanting to take a .JPG image and put it in a data frame. The image has to ...
3
votes
1answer
72 views

Add field and calculate expression with PyQGIS

I want to use PyQGIS to add a new field and calculate the values for each feature. Similar to Field Calculator option. My 'Field Calculator' expression is for example: y(start_point($geometry)) from ...
1
vote
1answer
55 views

How can I add output vector to a list to loop a geoprocessing on it?

I have about 18000 circular buffers that I want to transform in some kind of a wind rose. To do this, I transform the buffer polygons to lines, create points along the lines on the wind directions ...
12
votes
4answers
3k views

What software to use for making digital elevation models from UAV aerial imagery?

I have a an unmanned aerial vehicle (UAV, also called a drone) and would like to use it to create a digital elevation model and orthomosaic over a smaller area (10 km2). My platform will be able to ...
1
vote
0answers
28 views

Qgis standalone script using eclipse IDE -No data provider available issue

I know, the question is asked few times also in this forum but i apply already the presented solutions without any chance to fix the problem. In General the source code looks like the following: from ...
7
votes
1answer
228 views

gdalwarp producing empty rasters when called from python script, but not via command line

I'm trying to aggregate some raster data (all .tif) up to a coarser resolution (from .05 degree to .25 degree) using gdalwarp in python, but the command is not working. Instead of getting an output ...
1
vote
1answer
39 views

Read only geometry from shapefile with python/geopandas

I would like to open a shapefile(shp) without loading the attributes(dbf), as I do not care about attributes and it takes about 25 minutes to load them. I am using Python Geopandas, so any Fiona or ...
2
votes
0answers
42 views

QGIS Override MapTool Python

My goal is to reuse the default edits tool of QGIS in a custom App. But I would like to modify some behavior with those tools. For instance, the MapTool beside the Edit Node Action. iface....
2
votes
1answer
55 views

Cannot update GDAL on Ubuntu

I am trying to use GDAL_Rasterize and found a link via Google that suggested that there were bugs in older versions. I am running Ubuntu 14.04 LTS and 'gdalinfo --version' gives me 'GDAL 1.11.2, ...
0
votes
0answers
34 views

Map Service properties (arcpy): Error parsing recycleInterval, recycleStartTime, keepAliveInterval

Need to publish web service using arcpy tuning up its properties through the parameter names passing them into function parameter_toggle (see code bellow). We are getting error (see picture bellow) in ...
1
vote
0answers
24 views

Using Python to alter the drawing layers visibility within the layer properties?

Does anyone know of a way that I could programmatically access and alter the visibility of each individual drawing layer found within the layers property dialog box? The CAD file that I have brought ...
1
vote
1answer
40 views

Statistical properties of raster image

I want to calculate mean and standard deviation of a raster image. In QGIS zonal statistics, we can calculate it within a polygon or a vector. But I want to calculate for whole image. In Python I'm ...
1
vote
1answer
303 views

How to add mdb geodatabase layer (feature class) to QGIS project using python?

I want to add a feature class of a personal geodatabase mdb to my qgs project by code (a plugin I am working on). I am using a command similar to that used for adding shapefiles, but this is not ...
0
votes
1answer
30 views

How to Convert 3D coordinates to 2D coordinates for R-Tree?

I am trying to implement a bounding box via a spatial index and get an error stating that I must either provide coordinates as either (minx, miny, maxx, maxy) or (x,y) for 2D indexes. I've tried ...
8
votes
3answers
5k views

How can I convert a ascii file to geotiff using Python?

I have a ascii raster format file. For example: ncols 480 nrows 450 xllcorner 378923 yllcorner 4072345 cellsize 30 nodata_value -32768 43 2 45 7 3 56 2 5 23 65 34 6 32 54 57 34 2 2 54 6 35 45 65 34 ...
3
votes
1answer
44 views

Concerning “Error when using r.series in a python script”

Based on the "Error when using r.series in a python script" discussion (Error when using r.series in a python script) and based upon the solution given by Joseph I tried to make the same but by using ...
1
vote
1answer
1k views

Import shapefiles into PostGIS using Python with psycopg2

I want to import shapefiles to PostGIS using Python with arcpy & psycopg2 to connect to PostGIS. I don't want to use the existing shapefile loader nor ogr, just write my own import tool. My steps ...
1
vote
0answers
85 views

Converting Tiles to GeoPackage with Python

I am trying to convert tiles (as PNG) into GeoPackage. In order to do so, I would like to use the geopackage-python tool. However, when executing the script I run into the following error: Layer 14 ...
1
vote
1answer
55 views

PyQGIS Processing - wrong number of parameters

I'm facing a very strange issue from the module Processing using it as a standalone (not in the Qgis python console). I'm trying to create voronoi polygons from a points csv file. I can load my file ...
0
votes
1answer
191 views

GDAL2Tiles crash

GDAL2Tiles seems to crash after generating tiles and no google viewers are generated, but only openlayers. This is the code I run: os.system("gdal_translate -of VRT -a_srs EPSG:4326" + ...
1
vote
2answers
228 views

QGIS Deployment in Citrix XenApp - How to Remove Pyramids from the image “properties”

I'm working on deploying QGIS 2.4 into a XenApp 6.5 desktop environment. We've tested QGIS to see if it will run in our environment on a test server and have found that building Pyramids with ...
2
votes
1answer
67 views

Creating lines from points in QGIS, python error ( 'NoneType' object has no attribute 'asPoint' )

I'm new to Qgis, and i'm stuck with a simple task, I have a csv file with points coordinates and i would like to connect those points in a path, I'm working with wildlife movement data. i've tried ...
3
votes
0answers
50 views

Cannot Open VRT of Sentinel2 image created by AutoCreateWarpedVRT

I am trying to create a warped VirtualRaster of a Sentinel2 image. Here the Python code I tried: from osgeo import gdal, osr srs = osr.SpatialReference() srs.SetFromUserInput('EPSG:32632') gdal....
1
vote
0answers
62 views

Find Properties of Nearest Line Segments to Point in Python

I have been trying to find a workable Python approach to finding nearest line segments to one (or more) points, and came accross @gene's solution to @RustyMagoo's problem here: from shapely.geometry ...
17
votes
14answers
48k views

Full list of ISO ALPHA-2 and ISO ALPHA-3 country codes

I'm searching for a full list of countries with their country codes. Like on this page (need full and valid): http://www.nationsonline.org/oneworld/country_code_list.htm
4
votes
1answer
63 views

How to use joinFieldNamesSubset from QgsVectorJoinInfo class

I have written a small python plugin to refresh a join on a point layer in QGIS. point_layer=None for lyr1 in QgsMapLayerRegistry.instance().mapLayers().values(): if lyr1.name() == "point_layer": ...
1
vote
2answers
53 views

Preserving special chars using osgeo.ogr.Driver to shape in Python?

I am converting kml files to shape and run in to the problem that special characters are not preserved. #!/usr/bin/env python # -*- coding: utf-8 -*- import os, ogr, osgeo, sys ogr.UseExceptions() # ...
2
votes
0answers
25 views

Exiftool.exe in QGIS [closed]

I am developing a GIS application and in my code I need to call the exiftool executable which is installed in a different path. I call in command prompt and in OSGEO4W shell and both are fine. But ...
0
votes
0answers
37 views

How to store large time series of raster images

I've got a large time-series of rasters, derived from satellite images. Each raster is about 2000 x 2000 pixels, and I've got about 10,000 individual images. I'm trying to work out what the best way ...
1
vote
0answers
35 views

Per-pixel Descriptive Statistics through multiple Raster layers

I'm looking for computing different descriptive statistics through a set of raster layers. I need to make it with a python script capable of running into a standalone tool within the QGIS. To do so ...
4
votes
1answer
52 views

Register custom qgsfunction in pyqgis

I'm using pyqgis within an standalone application. Now I wanted to create a custom expression to to use as filterExpression in rule of ruleBasedRendererV2. (I tried the renderer with existing ...
3
votes
1answer
756 views

GDAL/OGR to extract pixel values

I want to use GDAL/OGR to extract pixel values from a raster file. The extraction should be based on polygons, stored in a shapefile. My thought was to iterate over the features in the shapefile, ...
3
votes
1answer
63 views

How do I fill a combobox in QGIS with only points vector layers?

I'm new to writing plugins in python for QGIS. I've been searching on internet for a while and wrote the code below. Even though it fills my combobox, all layers (rasters, vectors of all geometry ...
3
votes
1answer
469 views

WindowsError: [Error 32] The process cannot access the file because it is being used by another process

I am reading images (through Python) and do something. These images are placed in a temporary file, through: tf = tempfile.mkdtemp(filename) I use gdal to read the images (short part of code): ... ...
1
vote
1answer
797 views

Python package for processing GRIB data in windows

I would like to know what python packages can be used to process GRIB (Gridded Binary) data on a Windows system. I know there is something called pygrib, but to use that I have to use cygwin windows ...
4
votes
2answers
902 views

Configuring spatialite database access for Python 2.6.6

I'm working on a GIS project, and I would like to implement and test some geo-spatial algorithms in Python. For this purpose, I will not only need sqlite, but also spatialite, in order to store and ...
4
votes
1answer
67 views

point created on a line not precisely on line

In QGIS I wanted to move points to the nearest lines. I checked the solutions here, but there are using tools which I am not familiar to yet(POSTGIS or GRASSGIS). I figured out my own way with using ...