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

learn more… | top users | synonyms

0
votes
0answers
165 views

K-means cluster analysis with initial centers with ArcGIS

I am trying to figure out a way to make a K means cluster analysis using initial centers for the clusters, and only classifying, no iteration, for the features in a feature class. Having a feature ...
0
votes
0answers
58 views

mapserver python mapscript WFS RAW POST data

I installed mapserver-6.2.2 and python mapscript in a Openshift gear (Red Hat cloud). But I can't work with WFS RAW POST data. This is my configure line ./configure --prefix=$OPENSHIFT_DATA_DIR ...
0
votes
0answers
20 views

Adding features to rtree using multiple generators

I am using the Python Rtree library http://toblerity.org/rtree/ I know you can supply a generator function to the Index/Rtree constructor to insert features. How can I add features to an existing ...
0
votes
0answers
41 views

Python scripting error (merge)

I've been trying to use this code to add a field to, merge, and sort multiple feature classes: import arcpy, os arcpy.env.workspace = r'W:\S&P\s&p ...
0
votes
0answers
119 views

qgis python error mac yosemite

This problem started when I installed GRASS. I changed one environment variable and then I tried fixing, then I uninstalled egis, grass gdal, then I upgraded to Yosemite, and I still have the ...
0
votes
0answers
45 views

.tif image file cropping

following image has black portion on all side which value is [-inf], that what was showing in ENVI. How can I crop this black portion using that [-inf] value?
0
votes
0answers
37 views

Parameters for CreateFishnet

I´m about to create a fishnet. The CreateFishnet tool has several parameters, between them the origin_coord - The starting pivot point of the fishnetl. Data type of this parameter is Point Now there ...
0
votes
0answers
40 views

MakeXYEventLayer_management input problems

I have a huge excel data base and want to make with the tool MakeXYEventLayer_management a new Layer for my coordinates in my table. I found this code on ArcGIS help: try: # Set the local variables ...
0
votes
0answers
49 views

Training course for Fixing Data Sources

Is there any good training course or literature that pays some attention to automated processes for fixing map documents' data sources? Esri's little documentation on it is not helpful enough for me.
0
votes
0answers
78 views

Copy features to dataset in a loop

I'm trying to write script, that looks at the feature classes in a geodatabase and depending on what the feature class ends with, it creates a new dataset in the geodatabase and copies the file into ...
0
votes
0answers
30 views

increase gdal_polygonize's capacity

I have a very big PNG (82500x40900) with more than 5000 polygons. I want to use gdal_polygonize but the geojson output shows that gdal_polygonizer stopped after 4770 polygons and the geojson isn't ...
0
votes
0answers
29 views

Make qgsmapcanvas.extent().size equal to qwidget.size()

I am using Qt and QGIS to write a touch device application. Right now, I have a QWidget which is the child of qgsmapcanvas, and both of them are full screen. However, the qwidget size is (1920,1080), ...
0
votes
0answers
75 views

Rename open MXD using ArcPy

Is there a method in ArcPy that will rename the MXD you currently have open? I don't what to use mxd.saveACopy("newFileName") as it will just save a COPY with a new name, not rename the current MXD. ...
0
votes
0answers
94 views

name 'row' is not defined error appears overnight

I'm using arcpy with arcgis 10.0. I strived days to write my first python script, and I finally get what I want until...a new error occured during the night. The script you can see below was tested ...
0
votes
0answers
181 views

convert lat/lon into meters

I have a database of latitude, longitude and radius in the following form: 48180759,11518950,19.0 Both latitude and longitude are expressed in degrees but then multiplied by 1e6 (or 1000000). ...
0
votes
0answers
33 views

Python, KML and Geoserver/MapServer

I have this question: Is it possible to create an script that syncs a Google Drive folder with a MapServer/Geoserver and creates the visualization of an KML/HTML embedded file in Mapserver/Geoserver? ...
0
votes
0answers
51 views

workflow for reading data from a geotif and plot it on a map using python

I have a geotif with 9 different values (0-9). The projection of this is +proj=longlat +datum=WGS84 +no_defs. I want to plot this colors over a map (mercator projection). I tried to use ...
0
votes
0answers
86 views

Python Automated Data Driven Pages - Labeling and Symbolizing Errors

I am trying to use a python script to automate the export of the data driven map document. However when I use the python script to run the DDP export some of the labels are not exported and the ...
0
votes
0answers
190 views

Alternative way to find intersection between two polygon feature classes

I have been working on making a script to find intersection between two polygon ArcGIS feature classes. I can successfully do this in arcpy instances using spatial join or using search cursors on the ...
0
votes
0answers
62 views

gdalwarp hdf problem

I am working with some hdf files. The following command works on my Mac. command = gdalwarp -of GTIFF -tr 1000 1000 -tps -s_srs "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs " -t_srs ...
0
votes
0answers
52 views

Duplicate FID for points

I ran the script below to export feature vertices to points. The code works but after 10,000 points the same FID appears for 5 to 10 points. This continues for the rest of the dataset. There are ...
0
votes
0answers
40 views

Create polygon with predefined size

I need to create a polygon (square) with a predefined size to place it anywhere in a map. The size is given in centimeter (size in cm * scale = size in real). Is that possible with postgis/qgis?
0
votes
0answers
34 views

Help with Parameters in Qgis Python lidartools:lasclip

I am trying to use the python console to batch quite a few processes from lastools. I'm having trouble with the arguments to lasclip. clipped_las = processing.runalg('lidartools:lasclip', True, ...
0
votes
0answers
97 views

How to create a 8Bit unsigned geotiff with gdal?

I'm trying to create a 8 bit geotiff with this code. Is there a way to use unsigned 8 bit? I only found the GDT_UInt16? import numpy as np from osgeo import ogr, osr, gdal from osgeo.gdalconst ...
0
votes
0answers
117 views

Python script to get feature count in a feature class that fall in a specific range

I need to write a Python script to be used in the ArcGIS Field Calculator that will return a feature count for values between a range defined by the values in two other columns. Here is a graphic ...
0
votes
0answers
100 views

How to create shape file and show in map using Python?

I have calculated distance between longitude and latitude by import .csv using python. Now I want to show this points or distance between points in map. So how can i do it by using only python? and ...
0
votes
0answers
70 views

GeoTiff using GDAL

How can I project a tiff using the position of the four corners where the picture have to be placed ? The tiff file has no georeference or projection information I have got a grid with values of ...
0
votes
0answers
110 views

Is there a faster alternative to arcpy spatial join?

I am wondering if there is a more efficient method than using arcpy spatial join in python. I am currently using intersect and closest methods (including a near distance field) so would need ...
0
votes
0answers
25 views

Script Efficiency

I have been working on a script to read pixel values and print the values out to a text file. This is currently working but takes around two hours with a 6.6gb file. This is to replace an older script ...
0
votes
0answers
55 views

Custom ArcGIS Server Backup Python Script

I am running into issues using the ArcGIS REST API import site operation because the server is federated with Portal and only has a domain certificate (i.e. not CA-signed). So I am looking for a ...
0
votes
0answers
142 views

Fitting Time-series Values with a Gaussian Function

I'm trying to fit a stack of NDVI values to a Gaussian model to allow for determining dates of certain NDVI values using Python and NumPy/SciPy. I've attempted to do this with ...
0
votes
0answers
76 views

PyQGIS: Coordinates to TextBrowser

I'm doing tutorial: http://www.qgisworkshop.org/html/workshop/plugins_tutorial.html#tie-qgspoint-output-to-the-gui Right now I at the "Tie QgsPoint Output to the GUI". I did everything similar to ...
0
votes
0answers
169 views

Error selecting geometry from Oracle DB

I'm trying to select a geometry column as well-known text from a feature class in an Oracle 11g/ArcSDE 10.1 geodatabase. (This is in a Python script using pypyodbc as a driver). This is my SQL: ...
0
votes
0answers
75 views

export qgis point symbols per feature

I am currently trying to export the drawn symbol in a QGIS map for a webmap. The current approach is to extract the properties: if i.rendererV2().dump()[0:6] == 'SINGLE' and i.geometryType() == 0: ...
0
votes
0answers
28 views

Using setCrsTransformEnabled in qgis 2.4

In the python console of qgis I try to load layers from shapefiles which have different CRS using the setCrsTransformEnabled() method. To test the procedure I changed the CRS of the shape file using ...
0
votes
0answers
68 views

Export qgis vector layer to PostGIS table using python

I've created script, which does bunch of stuff, and in the end, I would like to export it to postGIS table. What is easiest way to do such export using qgis api or sth. I dont want to write SQL ...
0
votes
0answers
243 views

Looping through each specified layer and exporting PDF map showing shapefile and information taken from each attribute table?

I have been working on an automation script for the purpose of mass pdf exportmaps. I am required to export many pdf maps showing individual boundary files and client information. I believe the ...
0
votes
0answers
199 views

Exporting table of contents (TOC) with group layers to file geodatabase (GDB) with feature datasets?

I have created a python code for exporting shapefiles from TOC from different source folder to a output folder. However, in this case, I have some group layers in the Table of Contents (TOC). ...
0
votes
0answers
119 views

Why will my field calculator python function always return Null if any inputs are Null?

I've written a simple python function that outputs a list based on yes/no fields. It appends a word to a list if it has a "yes" in another column then writes the list, formatted according to the ...
0
votes
0answers
49 views

ArcPy PointToRaster memory locking error solution?

RasterToPolygon ArcPy tool fails on generated raster (PointToRaster) with large(?) dataset. Not sure the size threshold, but fails on 78559 x 33544 raster, completes on smaller datasets so I know ...
0
votes
0answers
86 views

imposm2 to imposm3

I am using imposm2 to import osm data to postgresql database but recently moved to imposm3 which is still in development phase I want to know is there an easy way to convert imposm2 mapping file ...
0
votes
0answers
74 views

ArcPy loops/cursors faster without background processing?

I am running loops via cursors and raw Python in Arcpy. Within the loop several processes take place, selecting shapefiles and converting to raster, clipping, running euclidean distances, extracting ...
0
votes
0answers
203 views

Proxy.cgi for OpenLayers reading from GeoRSS but not from WFS

I'm having trouble with the standard proxy.cgi written in Python for OpenLayers. Here's my code: #!C:/Python27/python.exe -u """This is a blind proxy that we use to get around browser restrictions ...
0
votes
0answers
79 views

Importing OGR from OSGEO Issues

I am trying to run: from osgeo import ogr and I get: Traceback (most recent call last): File "", line 1, in from osgeo import ogr File ...
0
votes
0answers
263 views

arcpy.Clip_analysis() Automatically Affixing .shp.shp to Target Feature Class Name?

I created a Python script that was previously working fine, but is now giving me an error when I try to use the arcpy.Clip_analysis() function in arcmaps 10.2.1. I am not adding a filetype extnsion ...
0
votes
0answers
63 views

geoprocessing service keeps running with default value for choice list, instead of actual value selected

I have published a script tool to a geoprocessing service. The tool is supposed to take user input from a choice list (which contains values from a field), and then that is loaded into select by ...
0
votes
0answers
71 views

Is there a way to take control of the mouse pointer using python in QGIS?

i would like to take control of the mouse pointer in QGIS (for example constraining it to a given extent...) Is there any way to do that using PyQGIS ? Thanks !
0
votes
0answers
96 views

Read large raster into lat/lon value pair to csv file using python and/or gdal?

I have rasters for all of the states in US and they are in 10m resolution. So they can be very large, containing a large number of pixels if you think of TX or AK. I want to create a csv file with ...
0
votes
0answers
86 views

Python + GME (Spatial Ecology) + Subprocess.call

I'm trying to run the following code to use GME commands in Python. import arcpy,os,subprocess from arcpy import env env.workspace=r"D:\Users\myname\Desktop\aviso" count=0 for raster in (f for f in ...
0
votes
0answers
79 views

Is it possible to rotate the geographic projection in QGIS?

I'm trying to create a Qgis plugin that can rotate a raster layer in the canvas so that the north is pointing at another direction. Unlike the thread below, I'm open to defining a custom projection. ...