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

learn more… | top users | synonyms

3
votes
0answers
67 views

Select shape file by their coordinates behavior [closed]

I want to filter shapefiles by their latitude behaviour: If the latitudes of edges A and B are constantly decreasing or increasing then this class marked as correct, otherwise are incorrect. For ...
5
votes
2answers
76 views

intersecting two shapefiles from Python or command line

I'm new to working with geo data and trying to address a similar question to the one posted here. Using command line and/or shapely/fiona, in the most efficient and elegant way. I created two ...
4
votes
2answers
92 views

Feature selection of one vector layer based on another vector layer: Any advice for using PostGIS or Python with PostGIS together?

I have a data set (more than 50,000+) of points (addresses) and Polyline (street) shapefiles of the area of interest. The street shapefile contains Annual Daily Traffic (ADT) data of each street for a ...
1
vote
1answer
42 views

Find Feature class into another feature data

I am trying to find the feature classes from one feature data into another feature classes. But I'm getting nothing at the end as a result of the matc statement. The data in the arc environment is ...
0
votes
1answer
37 views

PythEveryone in ArcGIS virtual campus module: Runtime error + namerror

I'm working on the python for everyone lab in Esri's Virtual campus. I'm running into an issue that I can't seem to solve, I've already done some outside research hence posting here! The issue: 4th ...
7
votes
1answer
129 views

ArcGIS Pro and ArcGIS Desktop on same machine

I'm just trying out ArcGIS Pro for the first time and I had to change my PYTHONPATH to python34 ArcGIS Pro is now working without crashing every time I open the Analysis tab... However now when I try ...
2
votes
1answer
98 views

Huge file size after averaging two rasters

After I average two raster files (singleband gray, initial size ~700MB each) through QgsRasterCalculator() with formula "Raster1@1/2 + Raster2@1/2" the size of the resulting file is 2.9GB. Is there ...
2
votes
0answers
38 views

Setting extent with python in QGIS in new project

Situation: I am trying to write a simple python script for QGIS. The script aims to load one layer, zoom to given coordinates, and 'place a pin' there. The coordinates will change every time I run the ...
3
votes
1answer
86 views

Automating QGIS project composer printing?

I have a large project which consists of around 10x QGIS project files. Each project contains a print composer with an atlas set up to generate up to 10 maps. What is the easiest way to automate the ...
3
votes
1answer
76 views

How can I use GDAL Merge with multiple files in one folder in python 2.7?

I have multiple files in one folder which are called, for example: 049.h11v11.hdf, 049.h11v12.hadf, 049.h11v13.hdf, 065.h11v11.hdf, 065.h12v12.hdf, and so on. The first 3 numbers are the julian day of ...
3
votes
1answer
29 views

Rasterize vector file using gdal.RasterizeLayer() without losing attribute information in python

Is it possible to rasterize a vector file with gdal.RasterizeLayer(), saving attribute as Raster value ? Like gdal_rasterize can do with : gdal_rasterize -a attribute -ts rasterXsize rasterYsize -l ...
1
vote
1answer
26 views

How to determine the appropriate size of a 3 dimesional raster block/chunk

I am processing large MODIS time series (12 years) dataset with an array dimension of 4800 x 4800 with a data type of Int16-Sixteen bit signed integer. I want to build a 3 dimensional array ...
3
votes
2answers
71 views

Trying to use LIKE in an expression in pyQgis

I want to select elements from a layer. The field I use is a string/character field. My code is as follows: expr = QgsExpression(ur'"myField" = \'valueILookFor\'') it = ...
8
votes
1answer
117 views

Mapnik failing through a webstack

I've inherited a pretty large / complex (internal) web portal (run under python 2.6.6), which has some mapping capabilities, using Mapnik 2.2.0. The web service has recently started throwing 500 ...
2
votes
1answer
43 views

Why is QgsFeatureIterator returning bad data?

This is a spin-off of another question I asked previously, but it has grown in scope and deserves its own thread. While iterating over a vector layer (shapefile) I included the following ...
1
vote
0answers
34 views

Arcpy 64 bit, or multiple Python installs? [duplicate]

I am new to Python, and am trying to write a script that would upload a geodatabase to AGOL at a scheduled interval. However, the 32 bit distribution of Python cannot handle the file size (over 1 MB), ...
3
votes
2answers
58 views

merge dbf files - without using arcpy?

I have ~7,000 2 line dbf files (zonal stats output) that I would like merged into one table, all the same format. I have used the arcpy merge tool, as seen below. #merge tables together ...
1
vote
1answer
49 views

QGIS Processing Script will not Validate

I am writing a script for QGIS that will identify the perimeter of only external facing walls. The script is probably a little long to post directly into this post, but it is accessible here. The ...
1
vote
2answers
32 views

Error Installing GDAL python package on a MAC

I'm trying to install the [GDAL python package][1] on a Mac (OSX Yosemite) with python 2.7.9 and GDAL 1.11.3. When I try $ sudo easy_install GDAL I get a number of errors ending with: error: Setup ...
3
votes
1answer
82 views

How can I install pyproj into arcpy?

I am trying to install pyproj. I downloaded the zip folder on this page pyproj. Then I copied it to C:\Python27\ArcGIS10.3\Lib. However, when I try to import the module, I get the message ...
1
vote
2answers
94 views

Calculate the distance between two coordinates (WGS84) in ETRS89

I have an Excel-file with coordinates (WGS84). From these coordinates I create a line (shp) and transform it (ETRS89). After the transformation I calculate the distance from point to point. I do this ...
4
votes
0answers
34 views

How to unproject a UTM ESRI Shapefile over WGS84 extent boundaries

Background: I am working on a collaborative project in Python that uses GDAL/OGR to work with rasters and vectors. There is one function that attempts to project one ESRI Shapefile's coordinate ...
1
vote
1answer
49 views

ArcGIS REST API: using /addItem with secure service

Using the code below I am adding 100's of Map Services to a portal instance. The script successfully runs and adds them however the services are secure so I would like portal to store the username ...
1
vote
2answers
121 views

GDAL Translate error

I am trying to translate a .tif raster to an xyz one in Python The code is this: import sys from osgeo import gdal gdal.AllRegister() driver = gdal.GetDriverByName('Local Disk') dataset = ...
2
votes
1answer
45 views

IPython plugin not loading in QGIS 2.12

I am trying to run the IPython console plugin in QGIS 2.12 After installing the plugin there is this error: You need to install IPython 3.1.0 (and then restart QGIS) before running this ...
1
vote
1answer
31 views

rendering a single band raster - new maxvalue

I have a downloaded raster file, singlebandgray, and want to redefine it's maximum value to a higher integer value. Via the QGIS desktop it works smoothly by hand via the Properties. Via Python I was ...
2
votes
1answer
39 views

Actions undertaken in QGIS somewhere stored as Python code in a history file?

For practioners that are not willing to study all the methods for each class, it would be helpful if every action that was made within the QGIS interface can be looked up afterwards in a history file ...
3
votes
0answers
56 views

Python script to automate saving satellite images for different combinations of lat longs

PS. : This is my first try at a GIS application. Please forgive if my question is pedantic. I would like to get a series of satellite images all at the same scale centered around different points on ...
2
votes
1answer
41 views

Getting script file name as object in Python? [closed]

I'm quite new to Python and programming in general. I have created a few scripts using Python and as part of my framework/best practices for my scripts, I'd like to display when a script begins and ...
4
votes
0answers
35 views

Custom forms feature in QGIS: QMessageBox.warning closes QDialog as well

I am attempting to work with custom forms in QGIS. I have included a validate function in my code that will pop a QMessageBox() if the validation fails. My problem is that when I close the QMessageBox ...
2
votes
1answer
48 views

Fastest way to process OSM

I've been working with imposm and Python, it's fairly easy to parallelise. I've hit a performance problem when trying to associate the coordinates to the ways. Creating a dictionary from ...
3
votes
0answers
29 views

How to fill composer map with canvas data and export composition as png using Python script in QGIS

I'd like to create a Python script in QGIS that will take some data (shp + tif, in my example), use tamplate Map Composer (from a file) and export created composition to a png image. With almost no ...
3
votes
2answers
80 views

How to know if your script is executed in a QGis python console or a standard python console?

Is there any way to detect if a script is computed on QGIS python console or not? Any attribute or test to assess it?
5
votes
2answers
106 views

Field Editor with Python: select by attribute & open table (QGIS)

I want to search on number values (>, = etc.) within a column and open table for the selected ones. It's an easy task within field calculator, but I want to see how it's build up with the python ...
0
votes
1answer
54 views

GeoJSON properties not displaying in popup-leaflet GeoDjango

I want to show the markers in different colours based on column(rastervalues) values in leaflet. My geojson data is {'geojson': "[{'type': 'FeatureCollection', 'features': [{'type': 'Feature', 'ge ...
1
vote
0answers
29 views

Python's logging module not working with ArcPy [duplicate]

I have a Python script imported to an Arc toolbox. It uses Python's logging module for message recording. import logging import time import os, sys import arcpy ...
1
vote
0answers
26 views

Geocoding and SRID - can't get desired SRID

I want to geocode a set of addresses with GoogleMaps, store the results in a different SRID in a PostGIS table and then project it in QGIS over a layer with the identical SRID. So far, it has been a ...
1
vote
0answers
15 views

Subset of spatial data causes slow SQL Spatial STIntersects

I've been having an issue were a specific subset of data causes an SQL spatial intersect to run extremely slowly. I am running sets of 1000 points against a subset of polygon layer to be intersected. ...
0
votes
0answers
5 views

Converting branched multiline into series of single lines in shapely

If I have a multiline string representing a branched network with a common starting point (red network in pic below): Does anyone know how to convert this to a set of individual single linestrings ...
1
vote
0answers
25 views

Interpolating every X distance along line in shapely

If I have a shapely multiline object that contains many lines whose total length each is 50km (when traced from the origin), and I want to interpolate along the multiline every X meters (let's say ...
1
vote
1answer
55 views

How to set conditional imports for a PyQGis script running inner or outer of QGIS?

Related to Cannot run standalone QGIS script I try to detect if my script is run from a QGIS instance (through the plugin) or from a standard python console. This detection will enable conditional ...
3
votes
1answer
85 views

Convert KML to GeoJSON with ogr2ogr and remove KML special attributes

I'm trying to convert KML to GeoJSON (or some other formats) using ogr2ogr (look at this question). It's all working fine, except of one problem: I'm always receiving these metadata attributes with ...
2
votes
1answer
40 views

Working Windows Standalone .EXE From QGIS Python Script Using PyInstaller?

Has anyone been able to create a working standalone Windows .EXE file from a QGIS Python Script Using PyInstaller?
2
votes
3answers
89 views

If else if in ArcMap Field Calculator (Python)

Help! I think I am missing something that must be so obvious but I cannot figure out what I am doing wrong. Here's some background: I have one short integer field (MajAspNum) that has values ...
7
votes
0answers
140 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 ...
-2
votes
0answers
57 views

Running geometry check on all feature classes in personal/file geodatabase?

Anybody can help me out and build a small Python script to do the following: browsing though all feature classes (also include those in the dataset) in an personal or file geodatabase Run the ...
5
votes
1answer
87 views

Sequentially distilling several raster data sets into final map based on first time pixel’s value was greater than zero?

I have 120 daily data sets that contain a value only when a threshold temperature is exceeded. If the threshold has been reached that pixel was assigned a value indicating the day it occurred (the ...
4
votes
1answer
77 views

Python and ArcObjects- Stand Alone open MXD

I wrote a Python function using comtypes to access ArcObjects to programmatically change the size of layer symbols in an MXD. The function works great in the Python window of ArcMap. However, I need ...
2
votes
1answer
39 views

Raster interpolation with gdalwarp, and replicating with python API

My problem is similar to this: Replicating result of gdalwarp using gdal Python bindings. I want to resample a raster using bilinear interpolation, which I've done successfully with gdalwarp by adding ...
3
votes
1answer
93 views

Using Update Cursor, populate field

I am trying to write a script where you preemptively select features and then go to the tool which calls for the feature layer and it populates fields based on certain criteria which I will design of ...