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

learn more… | top users | synonyms

0
votes
0answers
31 views

Does somebody has a Python script for a Multiple Flow Algorithm?

I´m looking for for a python script for a multiple flow algorithm assigning flow to all lower cells from one central cell. This algorithm is part of the r.terraflow and would need it to start at a ...
0
votes
0answers
77 views

I cannot use Gdal Python but Gdal binary, how to fix?

**UPDATE: I can use GDAL with Python now by copy file from Gdal 1.11 compiled folder to Python package folder. But the problem with GDAL binary and GDAL library that they cannot read the HDF5 file. ...
0
votes
0answers
33 views

QGIS Init script: How to get active layer and GID of edited object?

I would like to ask, how to get the GID (from Postgresql+Postgis) and active layer of edited object in QGIS Init script (Python)? Code Layer=qgis.utils.iface.activeLayer() works in QGIS console, but ...
0
votes
0answers
28 views

Error when providing “Python init function” from Python script in QGis

Problem I am using QGis and Python scripting to programmatically set a layer's edit form and "Python init function". It appears that the Python init function is called correctly in one situation ...
0
votes
0answers
17 views

maptool doesn't receive key pressed events after QMenu.popup() function

I am developing a qgis python plugin which uses a maptool receiving key pressed events by keyPressEvent(self, event) function, mouse movements events by canvasMoveEvent(self, event) function and mouse ...
0
votes
0answers
34 views

Suggested Workflow for Generating NITF (w/headers) using Nitro in Python

It seems that there is no real suggested workflow for generating an NITF file from a georeferenced, validated tiff image. Here are my inputs: Pyramid Tiff Image - georeferenced MBR (minimum ...
0
votes
0answers
58 views

What is GDAL equivalent to ArcGIS Lookup (Spatial Analyst) tool?

I use the Lookup (Spatial Analyst) tool to create several raster layers from a combined raster layer. I know ArcPy can do the same thing. However, I would like to know whether GDAL (or another ...
0
votes
0answers
63 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: ...
0
votes
0answers
30 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.. ...
0
votes
0answers
37 views

How to display image on the Graphics View widget in QGIS

The idea is to create a plugin in QGIS in which on mouse clicking of the layer show up a popup with Graphics View widget on it. I found this pure code: self.scene = QtGui.QGraphicsScene() ...
0
votes
0answers
24 views

Correctly drawing wrapping circles using Basemap

I've been using the functions from here to draw circles on maps using Matplotlib Basemap. This works well, but breaks badly at the map edges: I'm wondering whether there's a general approach for ...
0
votes
0answers
28 views

QGIS custom action and button toggle

I've added a custom toolbar with some buttons to my QGIS 2 plugin. I'm wondering how to get the button toggled (pushed in-looking effect) using Python. I've gotten this to work for existing buttons ...
0
votes
0answers
30 views

How can I define coordinate system of feature class when overlaying with coordinate zone shapes in python?

I have generated bounding boxes for CAD drawings that do not have defined coordinate system files. I know that each drawing could be in one of 6 different systems that do not overlap. My idea was to ...
0
votes
0answers
58 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 ...
0
votes
0answers
57 views

How to write a signal for 'QgsMapCanvas::xyCoordinates (const QgsPoint & p)'

I am very new to QGIS and python. I am writing some python codes and I want to capture the x,y coordinates when the cursor moves on the map layer. (same as the coordinates show in the QGIS status bar) ...
0
votes
0answers
20 views

Plugin Ui interation with user and Qgis

A downsized hydrogeology plugin ( GraphicGroundwater = GGW ), instilled the need of a reduced combination of Carson's 'fTools doVectorGrid' and Pryet's 'Qgridder' plugin(s). The GGW Ui is displayed ...
0
votes
0answers
23 views

ArcObjects open report wizard

ArcMap provides a function to create report and after the creation can save as different formats. Now I want to allow user to skip the intermediate steps. I want do the report function using ...
0
votes
0answers
36 views

SST process using just GDAL and python

I am trying to process AVHRR 4km SST using a module that uses only GDAL and Python. OSGEO4W, GDAL, Python 3.4.1. I borrowed the download code from http://grasswiki.osgeo.org/wiki/AVHRR. The rest I ...
0
votes
0answers
28 views

Use python to execute GRASS module v.distance

I have successfully use python script to import a line shapefile and a point shapefile into GRASS, and the following are their attribute table in QGIS, line shapefile point shapefile Now I want ...
0
votes
0answers
47 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
0answers
49 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
0
votes
0answers
37 views

Assigning a changing value to a variable in ArcMap using Python

I have two lists/strings/sets (the exact format can easily be changed) of values (u and v). I would like to do a cost distance analysis between each point in u and all points in v. For this, I need to ...
0
votes
0answers
42 views

QGIS r.cost / r.drain and errors

I'm trying to find the optimal route of several archaeological sites with Qgis 2.2. I've a DEM and r.walk or r.viewshed executed successfully, but r.cost always fails creating the output, however ...
0
votes
0answers
45 views

Error: Algorithm not found

I wrote the following code to start using QGIS APIs. in a very short script seen below, the "Error: Algorithm not found" is thrown and "print processing.alglist()" is showing none. Any idea why this ...
0
votes
0answers
68 views

PyQgis: How can I “append” features to a Vector Layer instead of “stacking” them?

As I read data from a text file and add points to a QgsVectorLayer (as well as related attributes), each new record gets 'stacked' on top of the previous one in the Attribute Table (and within the ...
0
votes
0answers
42 views

QGIS Python API: SAGA difference algorithm does not split parts

I am using the QGIS Python API in an external script (platform: ubuntu). Everything works as expected. But when I am using the SAGA difference algorithm within the script, there is an issue. The ...
0
votes
0answers
20 views

'QgsGraph' object does not support indexing

Traceback (most recent call last): File "/home/kasun/.qgis2/python/plugins/floogin/floogin.py", line 246, in findRoutes path=self.astar.astarSearch(graph, start, end) File ...
0
votes
0answers
48 views

Draw polygon shape via Qbutton

Ideally; a Ui Qbutton is needed that will activate map canvas editing to add a polygon via the classic mouse click rubber band method. Upon a concluding right click with the active tool, editing ...
0
votes
0answers
54 views

gdal commands cannot import gdal — possibly ArcGis conflict with separate Python?

I installed ArcGIS 10.1 yesterday (Windows 7 64-bit), but have a separate Python 2.7 installation. Now my gdal does not work correctly anymore. It looks like all commands without the ending "*.py" ...
0
votes
0answers
34 views

Running algorithm from Processing takes much much longer time than doing it from Qgis GUI

I have a road layer with 32000 features, and I want to find all the intersections. when I use processing->SAGA -> convert polygon/line vertices to points, it is all done in less than 5 seconds. ...
0
votes
0answers
25 views

QgsFeatureRequest and QgsVectorLayer are missing in the QGIS

I am trying to use part of my friend's code in eclipse written in Python using QGIS APIs. I have installed the QGIS 2.2.0. now before running the code, it seems that something is missing here. All ...
0
votes
0answers
27 views

Using GDAL and Python to replace the Spatial Analyst TabulateArea method

Has anyone tried to phase out the use of ESRI's Spatial Analyst extension by replacing it with GDAL in Python? I have a script that only uses the TabulateArea and Statistics analysis but not much else ...
0
votes
0answers
21 views

Is it possible to compile a grass addon into a commandline app?

As the title says, I'm wondering if it's possible to take the source of a grass add-on and compile it into a commandline app. Or, alternatively, use the add-on in a python program. Basically, I ...
0
votes
0answers
36 views

Shapely tolerance overlapping layers

Is it possible to define a tolerance by applying a overlapping layers function in Shapely?. I'm trying to make an intersection between lines but the results are different with other GIS. regards.
0
votes
0answers
43 views

when using arcpy.mapping.ListDataFrames, python.exe crash

Has anybody had this problem? this drive me crazy. what is wrong with my code? import arcpy mxd = arcpy.mapping.MapDocument(r"D:\Project\Project.mxd") df = arcpy.mapping.ListDataFrames(mxd, ...
0
votes
0answers
110 views

GRASS “Argument list too long” with run_command

I have a python script that uses the grass.script.run_command function to read in, process, and export a bunch of raster files. On projects where I have a large number of rasters (>1,000), it fails ...
0
votes
0answers
17 views

Is it possible to override default QGIS error messages?

From my QGIS plugin I would like to override/replace one of the the default QGIS error message ("The layer layer_url is not a valid layer and can not be added to the map") with my own message. I have ...
0
votes
0answers
52 views

leaflet.js ImageOverlay projection bug

We are trying to display a image over the leaflet map using "ImageOverlay" layer, but it does not work according to expected. In Google Maps works perfectly. Apparently Leaflet just stretches the ...
0
votes
0answers
22 views

Using file directories from a table

I have a toolset that will be using a set of file directories as input (shapefiles, raster etc.).The table can be stored as CSV or DBF. Is there a way to read a specific cell (row,column) of the ...
0
votes
0answers
53 views

Accessing CRS info with pyqgis in a standalone application

I'm currently having trouble accessing and defining CRS information in a standalone application. I've read the API documentation and have tried it in the python console in QGIS 2.2.0. The python ...
0
votes
0answers
52 views

How to pass a layer to a thread in PyQGIS?

In my plugin, a user can select one of the loaded layers from a combobox, and then proceed to do further calculations. These calculations are done in a separate thread, "WorkerThread". My question is ...
0
votes
0answers
25 views

Call post-execution command manually in QGIS Processing

I wrote myself a Processing Script and it works fine. Now i want to style the outcome (a new raster layer) with a QgsColorRampShader(), but the problem is that Processing adds the output after the ...
0
votes
0answers
61 views

MODIS products reprojection problem

I am working with Modis Data (MOD02, MOD11, MOD03, MOD35). In particular I am looking for a solution to reproject the subdatasets below, into albers equal. Because I am using NLCD dataset which uses ...
0
votes
0answers
49 views

error when running a script for clipping a group of rasters on QGIS 2.2 OS X 10.9

I'm trying to clip a group of raster files on QGIS 2.2 (KyngChaos) running on OSX 10.9. When I run the script I get an error related with the location of the files. I'm wondering if this is something ...
0
votes
0answers
27 views

Cross platform gis api wrapper architecture

I have a python class library for arcpy which is just a couple of class definitions with my common code attached as methods and properties so that I can do code centralization for stuff that I have ...
0
votes
0answers
47 views

need help creating layer with geom type wkbGeometryCollection with gdal python

I am having trouble creating a layer with geom type as geometry collection. Thus I am unable to output to a shapefile. I have attached code below. dstLayer is null. This does not happen if geom type ...
0
votes
0answers
43 views

Import GRD and count days with temperature below zero, using python

Minimum temperature of a year is given in GRD file. I have to read that file and then calculate the frost day (less than zero degrees celcius) from it in python.
0
votes
0answers
79 views

QGIS custom form and autocompleted data entering with “validation check”

In Qt-Designer I created a custom form where I can choose serveral species names in two different combo boxes. This combo boxes are linked via "value relation" to the accordant non spatial tables ...
0
votes
0answers
48 views

argument 1 has unexpected type error when parse xml with ElementTree

I try to make plugin in QGIS Dufour to open and read xml metadata and then to show some information in metadata i try to parse it with ElementTree and then show it in QLineEdit, this code i use to ...
0
votes
0answers
77 views

aggregate counts on a grid and determine grid adjacency

So here is the problem I am working on. I have a bunch of lat, long coordinates for event occurrences. I want to create a grid and then count how many points fall into each grid pixel. Once I have ...