Python is an open source interpreted programming language used in many GIS programs.
0
votes
0answers
7 views
Using Python & GDAL to reproject, merge, and finally clip a series of rasters
I'm trying to learn GDAL in conjunction w/ python to write what I believe is a simple script. My goal is to write a script that will look at a list of .tif files, reproject them to a single spatial ...
1
vote
0answers
6 views
Python grass : Region iterations for i.segment for hundreds of images
I'm using Grass 7 on Windows Vista on one computer and Windows 7 in another. I'm relatively new to Python and Grass.
I'm performing object segmentation (i.segment) on hundreds of 3 band tif groups. ...
4
votes
0answers
16 views
GP service (from Python script tool) time-out: how to execute finally: statement
I am trying to build more fail-management logic into some of the GP services I develop. As a part of the business logic within the script tool I build, I write some data into the tables and thus lock ...
0
votes
0answers
12 views
My own base map in Python MatPlotLib
Is it possible to load my own base map for plotting data with matplotlib ?
I have some problems trying to freeze a script with matplotlib Base map included.
Then I think that I can use my own base ...
2
votes
1answer
61 views
ArcGIS 10.0, Field Calculator Nested “if” statement with strings
I need to assign a value depending on other two columns "FNP" and "Teilort". I wrote this script in phyton for field calculator. It actually works (no error messages) but the results in every cell is ...
1
vote
0answers
16 views
Correcting addresses with Python or Model Builder
I am trying to make a script in Python (or Model Builder in ArcGis) for correcting addreses.
I have a point-shapefile with the addresses(it has X and Y, and fields for "number","street name" and an ...
4
votes
0answers
26 views
Why are my calls to SAGA routines in QGIS not working cross platform?
The ability to write and run cross platform scripts for QGIS using python is OUTSTANDING! However, when I try to run the SAGA:MERGESHAPESLAYERS algorithm from a python script it works fine in windows ...
1
vote
1answer
8 views
Converting a large set of GRIB2 files to a spreadsheet format
I have a large set of .pgrbf.grb2 files downloaded from the Climate Prediction Center, and I need to run statistical analyses comparing them to data that was produced locally in .csv format. I'm not ...
2
votes
1answer
36 views
problem transforming CRS coordinates EPSG:4269 to EPSG:3857 in python
source_srs = get_srs(shp_abs_path)
source_epsg = source_srs.GetAttrValue("PROJCS|GEOGCS|AUTHORITY", 1)
print source_epsg #prints 4269
target_epsg = 3857 # because these are headed to leafletjs
...
1
vote
1answer
32 views
Flip multiple vectors
I have thousands of vectors (out of several million) which I need to flip, for directionally purposes (Network Analyst).
The vectors in question are annotated that their directionality has been ...
3
votes
1answer
79 views
How to convert polygon singleparts to multipart?
I am interested in converting a polygon featureclass created via Create Fishnet (Data Management) from singleparts to multipart. QGIS has a very handy tool called singleparts to multipart that does ...
1
vote
1answer
32 views
Open Source solution to batch process shortest distance between points along road networks?
I have the locations geocoded (a fairly small dataset, ~200), but I haven't found a straightforward way to load these data into QGIS and perform network analysis between nodes.
Any suggestions?
...
0
votes
1answer
52 views
Creating Map from database
I'd like some general advice on creating dynamic maps. I have a bunch of geospatial points saved on SQL Server and PostgreSQL. People will enter new data into the databases from time to time. I want ...
2
votes
0answers
31 views
Nested conditional statement in ArcGIS on multiple rasters
I have a list of rasters in ArcGIS and I want to create a conditional (Con) nested statement based on all the rasters of the list. I want this to be reiterate through many lists (the number of rasters ...
6
votes
2answers
94 views
How to delete last part of a string starting from the 1st non numeric character? [on hold]
I'm looking for a Python or VB statement to remove the last part of a string, starting where the 1st non Numeric character appears, using the Field Calculator in ArcMap 10. Here's an example of what ...
2
votes
0answers
38 views
ESRI's arcpy.asShape Bug Workaroud
I've read a few questions\answers about converting a GeoJSON Shape into an ESRI Object. Unfortunately most of the answers pointed to decoding the dictionary's manually. Im a bit more lazy ^_^
The ...
4
votes
0answers
44 views
Focal Statistics for Special Neighborhood
I'm looking to calculate focal statistics for each cell of a raster, within a neighborhood of a specified criteria.
Background - I have three binary rasters, each representing a single vegetation ...
0
votes
1answer
40 views
Arcpy MakeFeatureLayer won't loop
The arcpy snippet below is intended to generate a list of technicians scheduled to work, and then generates a unique feature layer for each one's stops.
fcStops = "TaskRoute1016"
with ...
1
vote
1answer
39 views
Fiona + Shapely: Loading a set of LineStrings and writing their centroids to a shapefile, including original properties
I am using Fiona to load a set of links. Then I put every record in the rows of a pandas.DataFrame. Then I get the centroid of each LineString using Shapely. Afterwards I copy the schema of the ...
1
vote
1answer
22 views
To count the number of grid cells in a polygon
Is there a way ArcGIS could count the number of 15 seconds interval grid cells within a polygon after digitizing that polygon. This should appear in a field in the attribute table.
3
votes
1answer
45 views
Track progress of running GP tools (in arcpy script) programmatically?
I am looking for a way to access a GP tool running in a Python script to be able to see how much of the job has already been done.
I have found SetProgressor function, but this one is available only ...
1
vote
2answers
20 views
Store result from gdallocationinfo as variable in python?
I'm trying to find the cell value at various locations in a set of rasters so I can count the number of real values (aka not 'NoData') that occur at all the locations within each raster. The various ...
1
vote
1answer
27 views
arcpy geometry difference method fails randomly, not reproducible
Has anyone used the arcpy geometry difference operator and seen unexplained errors like what I describe below?
This is a weird one. I have a script that calculates wave fetch distance in the ocean ...
2
votes
1answer
44 views
Changes not persisting in ArcGIS server from ArcPy script
I have a ArcPy script (10.2) which connects to an ArcGIS server (10.0) versioned geodatabase, via an .sde connection file in the Database connections directory through direct connect on a Microsoft ...
3
votes
1answer
36 views
shapefile record has no attribute named bounds
This is the code that I am using to print the bounds of my shapefile.
import fiona
import shapely
import shapely.geometry
fc = fiona.open("filename.shp")
shapefile_record = fc.next()
shape = ...
2
votes
0answers
22 views
How to move feature using python console in qgis
i am building a qgis plugin that can create as much copies of a selected feature as i want and paste it to my required place.I have code to copy and paste feature
select = ...
2
votes
1answer
33 views
Subset multiband tiff file using Python in ArcGIS for Desktop?
I have a multiband (1140) raster file in tiff format that I want to subset. Each band represents a prediction of monthly rainfall (12 bands for each year) over a 95 year period (12*95=1140). I need to ...
3
votes
3answers
72 views
Proj4 String for NAD83(2011) / Louisiana South (ftUS)
I am trying to convert NAD83(2011) / Louisiana South (ftUS) coordinates on survey maps to lat/lon using Proj4 in Python. Could someone please point me in the right direction? I have verified the ...
3
votes
1answer
48 views
How to automate the execution of Python console in QGIS?
I want to automate the execution of Python console in QGIS.
When QGIS start, Python console starts at the same time, and Python predetermined code wants to be executed.
Currently, in order to ...
2
votes
0answers
31 views
Kartograph - AttributeError: _Shape instance has no attribute 'parts'
I am new to GIS and open source tools so bear with me.
I am attempting to convert a shapefile to svg by running kartograph as a python module. My source file (provided by the local GIS ...
1
vote
0answers
16 views
How to freeze basemap using cx_freeze
I am working on script and I´d like to freeze it for use in computers who haven´t got python installed-
The script is running fine but when I freeze it using cx_freeze an error appears.
at the end ...
3
votes
1answer
55 views
How to interpolate point set into value contours?
The task at hand is creating 'isobar' value contours that are interpolated from a point field of either the x/y or lat/lon type. I really have no idea where to begin mainly because this needs to be ...
1
vote
1answer
29 views
Why does QGIS 2.4 prematurely terminates on startup at “Starting Python”?
After reinstalling Python (3.4) as well as OSGeo4W and the QGIS Standalone, my QGIS always crashed when "Starting Python" on my Windows 7 64 Bit environment. It simply did not proceed to the QGIS-GUI ...
2
votes
1answer
43 views
How to highlight a selected feature in qgis plugin using python [duplicate]
I am working on a plugin that would highlight selected features by changing its color. Can anyone help me regarding this? The plugin searches for required features and once found, I want that those ...
2
votes
1answer
44 views
Grass python shell can't find my data
I have a problem, the grass 7 python shell can't find my rasters. When I try to use the python shell to set the region I get "region not found". I can set the region in the grass shell (and other ...
3
votes
0answers
51 views
How to use Python's subprocess module in a for loop ArcGIS?
I am new to scripting with ArcGIS. Although this might be an easy thing for those that know Python well, I haven't figured out how to use the subprocess module in a script.
I actually think I'm ...
3
votes
1answer
29 views
Fast editing in raster calculator — adding lots of images together
I need to add pictures together to make one raster.
In one session there are about 150 - 200 images (named randomly - see syntax example)
There are about 50 sessions (named randomly)
Is there any ...
1
vote
0answers
19 views
Error message retrieval in grass.read_command
I am writing a script in python and would like to get the error messages of the function grass.read_command for seeing why it is not working.
I found a webpage ...
2
votes
0answers
27 views
Validating whether PDF file Uploaded as part of Geoprocessing Service of ArcGIS for Server?
I'm developing a Python Script that will be a Geoprocessing Service in ArcGIS Server.
The Script has a Button for Users to OPTIONALLY Upload a PDF.
If the PDF is Uploaded it gets inserted into the ...
2
votes
2answers
56 views
Delete Files in a folder and disregard an existing folder
I did a lot of geoprocessing steps and now I like to delete all files including result, result, _join, _interand reclass in their name. I can do this with the following
import glob
import os
for fl ...
2
votes
0answers
42 views
Setting up pycharm and pyqgis
I wanted to set up PyCharm the way this blog post explains: Setting up PyCharm for PyQGIS and Qt.
I did everything as instructed but when it comes to adding project interpreters I got this error ...
1
vote
1answer
66 views
GDAL Merge Py on ubuntu many tif tiles into one
Trying to batch mosaic 1000 Gtiff rasters (all in one folder) into 1 tif file via gdal merge on linux ubuntu.
Having trouble with the linux command to bring in all files in a folder into the output ...
1
vote
0answers
34 views
Clip Management won't work outside of ArcMap/Catalog
Okay so I have a short script that works great when copied into the python window in ArcMap or ArcCatalog, but refuses to work in an IDE for whatever reason. The error that is thrown is:
Start Time: ...
1
vote
0answers
18 views
PyQGIS add WMS layer with CQL_FILTER
Dear PyQgis developers!
I'd like to add a wms layer through my QGIS plugin with custom CQL_FILTER parameters. My problem is, that the CQL_FILTER parameters are not passed to the GeoServer (I can see ...
2
votes
1answer
75 views
PostgreSQL PL/Python Function for “IF” condition
I am using OSM data to run pgr_dijkstra() for weighted-path calculation. But after importing OSM data using osm2pgsql I populated some additional columns into my ways table. So now it has gid serial, ...
2
votes
1answer
61 views
Georeferencing a raster using GDAL and Python
I want to georeference a raster using python and GDAL. My current approach is to call gdal_translate and gdalwarp using os.system and an ugly list of ground control points. I'd really like a way to do ...
0
votes
2answers
89 views
How to replace “NULL” value with Zero in geodatabase table
I am having some difficulties to change "NULL" values to zero in geodatabase. I have tried the following python code in ArcGIS 10.1 ("Field Calculator") to change it but did not work.
Pre-Logic ...
2
votes
2answers
61 views
ArcGIS Field Calculator passes European-formatted floats as tuples
I have the problem that the ArcGIS Field Calculator doesn't read a float number when it is included in Python by referring to the field that contains it, in this way:
def giveResult(floatValue):
if ...
0
votes
2answers
61 views
How to calculate cumulative decreasing sum in ArcGIS-10.2
I need to calculate the cumulative decreasing sum as a field in a rasters attribute table. The cumulative decreasing sum is the accumulation of the last values to the first, so it is basically the ...
2
votes
2answers
86 views
Is there function to know raster's size (KB) with python?
I am doing a script with python. There I have many raster in canvas (qgis) which I would like to know size in MB. I am searching function in python that give me this result. Anyone can help me?
Thank
...
