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

learn more… | top users | synonyms

3
votes
3answers
119 views

Flow Accumulation Python or MATLAB Script

does anyone have any resources for scripting the D8 flow accumulation in Python or MATLAB after obtaining the flow direction values and filling in 'sinks' and 'flats'? I am trying to understand the ...
1
vote
0answers
34 views

Geopandas and Zonal Statistic error

I'm trying to use Geopandas and zonal statistic from rasterstats together. I run this code import geopandas as gpd from rasterstats import zonal_stats geodf = gpd.read_file("SHP/shape.shp") ...
1
vote
0answers
25 views

How to access a feature by specific field value using python ogr in an efficient way?

I have one point layer and one polygon layer with a common identifier attribute field called "N", same objects amounts but Feature IDs (FID) are different. I have to iterate on the first one and what ...
1
vote
0answers
44 views

Write stacks of GeoTIFF image to NetCDF by chunk size

I converted GeoTIFF image stacks to NetCDF file with writing data by Index order (layer by layer) (this was a main solution). The problem is, it writes slowly (for 25 images - 1Gb takes 90 seconds). ...
7
votes
2answers
435 views

Setting in memory workspace in ArcPy?

Is it possible to set env.workspace = "in_memory" in ArcPy using both ArcGIS Pro and the ArcGIS 10.2.2 (or 10.3) architectures? What I am trying to do is get the output of a snap pour point operation ...
1
vote
0answers
24 views

python handling case in address (lowercase,uppercase) [closed]

I'm working on a python QGIS plugin that takes address info from specific parcels and lists it. The data is in various cases and I'm having a hard time figuring out how to handle it. My data looks ...
1
vote
0answers
20 views

Use QgsRasterCalculator inside function

I'm relatively new to PyQGIS and I'm building a plugin to perform some calculations between raster layers. I was seeking to use the QgsRasterCalculator for some operations, however, the function that ...
3
votes
1answer
83 views

PyQGIS couldn't import processing with standalone script

I have read post on that problem (How can I access `processing` with Python?). The problem is even with a sys.path.append it couldn't import processing. It's very frustrated. Here the code I used: ...
1
vote
1answer
185 views

Is it possible to embed python in code for HTML pop-ups in ArcGIS

When clicking a point feature I want the html pop-up to not only contain the data of the feature but also an image of the area around the feature (eg. 2 km^2). The image would have the base map 'World ...
1
vote
1answer
44 views

How can a gui be visible till all the work is done? (QGIS plugin, python)

I have made a gui in QTDesigner for QGIS2.10 (as a dialog). When I press the ok button (if result == QtGui.QFileDialog.Accepted:) everything works just fine but the gui is disappeared. How can the gui ...
0
votes
0answers
25 views

Using shapely and geopandas in pyspark

I'm new to geospatial analysis and trying to start the analysis using python. Since the data is too huge we are trying to make use of hadoop and spark to get the analysis done. I would like to use ...
0
votes
0answers
15 views

Specifying column to use in gdal_translate

I am using gdal_translate to convert a tiff file to ascii. The tiff file is a raster with several columns in the attribute table. How do I specify to gdal_translate, which column to use for creating ...
0
votes
0answers
36 views

How to “INSERT” geotiff raster into postgresql/postgis table using psycopg2

Can someone walk me through how to insert a raster into a prebuilt postgresql table? The table has the postgis extension and I've been able to manually insert a geotiff using raster2pgsql, but I would ...
2
votes
2answers
46 views

How to use a raster dataset to store and retrieve “vectors” using Python?

For a rectangular area of about 10 by 10 km, I have to generate "displacement vectors" with a 1 meter resolution (i.e. a delta-x and delta-y value for every cell of 1 by 1 meter). I will end up with ...
0
votes
0answers
103 views

Dissolve using GDAL/OGR?

I'm trying to dissolve a shapefile according to a field using Python 2.7.10 and I don't want to use ArcPy. Is there any way to do this using GDAL/OGR?
2
votes
1answer
123 views

How to create Multipart Line in ArcPy?

I am developing some script in where I want to create custom multi part line geometry. For example: I have one array of points a,b,c and second array of points e,d,f. I want to create one ...
6
votes
2answers
1k views

Dissolve polygons based on attributes with Python (shapely - fiona)

I have been trying to create a function that does basically the same thing that the QGIS "dissolve" function. I thought it would be super easy but well apparently not. So from what I have gathered ...
0
votes
0answers
15 views

Processing toolbox on Google Cloud compute engine

I used to utilize the QGIS python console to do GIS batch operations on my local machine. Computation times are rapidly increasing and I'm switching to a virtual machine solution in the cloud. How ...
2
votes
1answer
53 views

scripted new items don't show up in the “items” panel on Print Composer

I'm working on a plugin to put items (pictures, arrows, labels) on the Print Composer layout. Everything is fine with the exception that once they are drawn, none of these show up in the "items" panel ...
1
vote
1answer
73 views

Problem intersecting shapefiles with OGR in Python

I'm trying to move from ArcPy for geoprocessing. After searching some posts, I found this solution using the GDAL/OGR Python bindings: Intersect Shapefiles using Shapely My code is the following: ...
3
votes
1answer
43 views

Calculate field interactively in ArcGIS ModelBuilder

I'm working on a closest facility analysis using ArcGIS ModelBuilder and want to enable the user to add locations interactively (i.e. typing the coordinates, not using any table or source layer). I ...
0
votes
2answers
455 views

How do you project a single pointGeometry object with arcpy?

In an attempt to create a python add-in for ArcMap that will pass x and y coordinates to another application via a mouse click on the map, I have written the following code. I need the output ...
1
vote
0answers
14 views

Python error using Vterrain in qgis

Traceback (most recent call last): File "C:/Users/Haimwant/.qgis2/python/plugins\VTerrain\VTerrain.py", line 409, in run os.spawnl(os.P_NOWAIT, directory + "Enviro", ...
1
vote
1answer
223 views

Use subprocess module to run GDAL processes

I am attempting to write a python script that will reproject & merge a series of .tif files then clip them to the boundaries of a .shp file. I have successfully used the subprocess module to do ...
1
vote
0answers
28 views

How to transfer SQL variable back into ArcPy script?

I am trying to teach myself to access sql tables from an ArcPy script as part of a larger QC process. I am fairly certain that I have the sql statement correct (it works in SQL Server, albeit ...
3
votes
1answer
49 views

What is the real distance between positions?

I'm trying to figure out what gives me the most accurate distance between positions: Measuring distance directly (using Geod.inv)? Measuring distance in 2D (after transforming to a known projected ...
2
votes
2answers
106 views

Convert SHP to MIF/MID

I am new to GIS and QGIS software. is there a batch convert function to convert shape files to mif/mid files? I can do this individually but is time consuming? I have seen similar functions on here ...
2
votes
1answer
492 views

Extract data from netcdf using a shapefile polygon in python

I need to subset data from a NetCDF using a specific shapefile. The data are sea surface temperature and ocean color at 1/4 degree resolution. I have 4 polygons describing the US. Northeast ...
5
votes
3answers
490 views

How to iterate over selected layers?

In QGis by pressing shift or ctrl keys a user can select multiple layers in toc. How can a python script iterate over selected layers? I'd expect a iface.selectedLayers() but I can find only ...
1
vote
1answer
67 views

Unable to delete folder with shutil.rmtree() with ArcPy?

Unable to delete folder with shutil.rmtree() with arcpy: del mxd delete_path_edit = temp_dir_extended delete_path_forward = '/'.join(delete_path_edit.split('\\')) shutil.rmtree(delete_path_forward) ...
3
votes
2answers
163 views

QGIS Real-Time Vector Layer Creation

QGIS is running on my Windows PC. I have a CSV file on my PC that is remotely updated with values of a raspberry pi's GPS position and some wifi-status readings. Is it possible to create a vector ...
2
votes
2answers
100 views

QGIS - Edit UI - Python Init function localization with QTranslator

I have interaction with layer defined by custom designed UIs. In layer properties the following settings: Edit UI - customline.ui and Python Init function: customline_set.py. #coding: utf8 from ...
8
votes
2answers
6k views

Using the python shape library: pyshp - how to convert .csv file to .shp

I am trying to understand how I can use the csv module in python to open a csv file in the same folder as the python script, and then create a shapefile using the shapefile module pyshp. The csv file ...
0
votes
1answer
45 views

Clip a GeoTIFF with a shapefile

I have a geotiff file that I made in Python. After I created it, I realized it was for all of North America, and I only need the United States. I then got what I assumed was a shapefile of the US ...
0
votes
0answers
30 views

Using ReportLab with Arcpy

I am having some issues creating a geoprocessing service within arcmap using reportlabs pdf toolkit. The script runs fine outside of arcmap, however when i try and run it in arcmap to publish as a ...
1
vote
1answer
26 views

QGis's Geoferencer plugin isn't able to paste text that has been added to clipboard via mac's “pbcopy” terminal command [closed]

I recently wrote a little script that used mac's pbcopy shell command to dump some text into my clipboard for use in the QGis program. For some reason I'm not able to paste stuff that has been ...
0
votes
1answer
854 views

ArcPy SelectLayerByAttribute Error and how to control GP output layer name

I have a rather long script that first references an SDE feature class (non-versioned) to perform an attribute selection based on a SQL selection. I've been unable to figure out why the script will ...
1
vote
0answers
44 views

My QGIS keeps me returning “ascii” errors

First, whenever I open QGIS, the extension "Quickmap Services" returns an error with this tracker: Impossible to load the extension quick_map_services because of an error on calling InitGui() method ...
17
votes
7answers
16k views

What's the easiest way to install GDAL and OGR for Python?

I am developing some code in Python, and I want to use the GDAL/OGR python bindings. What is the minimum I need to install in order to use the GDAL and OGR python bindings? What is the easiest way to ...
4
votes
1answer
59 views

Troubles Updating Values of Selected Features with PyQGIS

I'm trying to modify the values of some selected features. I have written this small script, but nothing happens when I run it. Is there something that I've missed? Here is my code. lyr = ...
1
vote
1answer
183 views

how to find that two polygon boundaries are connected?

I am working on a Geo database where there is a layer with many polygons. i want to find that which polygons boundaries are connected so that i can create new layer from that polygons and which ...
-1
votes
2answers
30 views

applying equation on numpy array [closed]

I have a 3d array: >>> print arr_u_num [[[ nan -3.90149352 -1.21931281 -1.08366785] [ nan -2.07270503 -1.27621887 -0.89320389] [ nan -2.69618795 -1.98927089 ...
1
vote
0answers
14 views

Check if raster (part of) lies within a pixel - python/gdal

I have a list of the corner-coordinates of different grid (c1-c4). First i sorted these coordinates to ensure that c1 is the upper left, c2 is the upper right, c3 is the lower left and c4 is the lower ...
2
votes
1answer
322 views

How to calculate the bounding box of a geojson object using python or javascript?

I have django app with leaflet on the frontend and I need to calculate the bbox of my geojson so that I can pass the coordinates to map.fitBounds() I have tried some libs as well as trying the ...
1
vote
0answers
39 views

Using row.setValue() function with arcgisscripting?

I'm trying to set column value NULL. I tried different ways to do that, I used CalculateField_management and UpdateCursor but i did not get success. when i used CalculateField_management. It says ...
0
votes
1answer
1k views

Unable to delete .lock files within geodatabase using Python

I am having problems with lock files while deleting and copying geodatabases for backup purposes. print "Now backing up your Geodatabase" distutils.dir_util.remove_tree ...
3
votes
0answers
50 views

Can qgis_core work with external python distribution (such as Anaconda)

I've built a python application that has its own GUI. I am trying to use qgis.core/pyQgis as map display only. Specifically, i have anaconda 2.3 and its python distribution (2.7.11). After research on ...
3
votes
1answer
80 views

Standalone python app that uses qgis.core without installing qgis

I am building a standalone python application that has its own GUI, i would like to use qgis.core/pyqgis as a map display only. But from my searches on this topic, it seems that the full QGis software ...
1
vote
2answers
277 views

QGIS (2.6.1) Processing Toolbox InitGui error after changing /Users folder on Mac OSX

Using: Mac OSX 10.9.4 QGIS 2.6.1 I had previously installed QGIS and had everything running fine. My home directory was /Users/rowanbado, but for another application I had to add an ...
1
vote
3answers
880 views

How to split a line every x meters? [duplicate]

I have a line that I ultimately want to split up into points. The points should be every 100 meter along the line. So I don't want to extract the nodes. Are there any OpenSource (QGIS, Python) tools ...