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

learn more… | top users | synonyms

6
votes
1answer
2k views

Python script to convert Lat Long to ITM (Irish Transverse Mercator)

I am currently working on a python script that converts Lat Long (WGS84) to ITM. The formulas used are based on a paper published by OSGB (Appendix B): http://www.leica-geosystems.co.uk/downloads123/...
7
votes
2answers
266 views

How to use Python modules as data providers in QGIS?

I don't know if it has been discussed in the past. Qgis data providers can only load native compiled libraries (.so, .dll, etc) that offer the required interface. Is there anything blocking to the ...
1
vote
1answer
550 views

Dealing with overlapping grid (and graticule) labels using ArcGIS Desktop?

I do have 3 grids made by Grids menu in data frame properties. I do need them to display in certain layout and also to behave accordingly to couple of rules (overlapping issues), so I can automate my ...
3
votes
2answers
3k views

Is there a way to automatically populate attribute values in Qgis?

We are using a Qgis/PostGis combo, where Qgis is used for adding features with attributes into a PostGis database. Is there a way to automatically fill the attribute data values with some project-...
1
vote
0answers
630 views

Looping multiple feature classes in python for ArcGIS 10

I currently have a large list of feature datasets and dbf tables in a geodatabase. I would like to join each shapefile with its matching DBF file using the "Make Query Table" tool and then "Copy ...
7
votes
3answers
4k views

Calculating Field based on name of shapefile

I have multiple shapefiles to which I have added a field called "gname1". For each shapefile, I want "gname1" to be set to the name of that particular shapefile, for every feature within the shapefile....
5
votes
1answer
7k views

How to reference to my in QT designer created buttons in python?

I have created an own QDialog (using QT Designer) with three push buttons. The names are pButton1, pButton2 and pButton3. After saving the file I created the ui file by using: pyuic4 ui_xy.ui > ...
3
votes
3answers
3k views

How do I access a multivalue output parameter from a python script in modelbuilder in ArcGIS 10.0?

How do I access a multivalue output parameter from a python script in modelbuilder, in ArcGIS 10.0? That is, I have a python script that loops through several input layers (working through a script ...
4
votes
1answer
279 views

Using “Time” in PostgreSQL/PostGIS

I'm loading gtfs into a PostGIS database but I have a problem when loading my stop times. The data is in time format "00:00:00" and for 1:30 AM, my data is "25:30:00". I'd like to keep it that way ...
7
votes
2answers
2k views

Is it possible to set the project CRS via a Python command?

Is it possible to set the project CRS of an project in QGIS via a Python command? In my case the user should open QGIS, start my written plugin and this plugin will load spatial data which has the ...
2
votes
2answers
352 views

Import 3D Analyst Functions With ArcPy

This really seems like I should be able to do this, but I am running into some problems. I want to import 3D Analyst with python, in the same way I import Spatial Analyst, but it returns an error. ...
0
votes
2answers
1k views

How to auto populate multiple fields in attribute table?

i have Feature Class Buildings, it has many fields like this(Region_id,Region_Name,Governorate_Name,Village_id,Village_Name,Building_id,Building_Hight,Building_Area,Floor_No,...
1
vote
0answers
224 views

Reclassify ESRI raster using variables from a SQL Server table (connecting with pyodbc)

I am trying to reclass a raster file (representing vegetation and land cover types) using values stored in a SQL Server Express table. I have a set of ~400 species records in the table, with each row ...
2
votes
1answer
861 views

Passing Python List into Layer Definition Query

I am trying to cursor search through items in a field in a feature class. I then append them to a list called 'distList'. I want to pass the list into layers definition query as follows: "Field1" in ...
6
votes
1answer
269 views

Is it possible to install ArcPy on a server for use in a geoprocessing web application?

I have a website on my intranet and would like to add some geoprocessing tools to it with Python. Many Users will not have Arc installed on their machines, so I think I need to have the scripts set ...
11
votes
1answer
2k views

Is it possible to take automatically a “screenshot” of my selected feature in QGIS using Python?

I am implementing a tool which allows the user to create a profile of a feature (in my case urban parcels). Is it possible to automatically take a screenshot of this feature (displayed in a QGIS ...
1
vote
1answer
474 views

python qgis.core import error caused by spatialite libraries

I have an error when trying to import qgis.core to python: >>> import qgis.core Traceback (most recent call last): File "<stdin>", line 1, in <module> ImportError: /usr/lib64/...
3
votes
1answer
263 views

Make feature layer based on date, problem in python

I have a script that runs, that takes the current date and puts into a variable, which is then used to calculate the date field for the Feature class. Further in the script, the updated data is loaded ...
10
votes
2answers
5k views

Finding mid-point of line using ArcPy?

How can the mid-point of a polyline be calculated using the geoprocessing framework, e.g. in a Python script? The polyline.centroid property returns The true centroid if it is within or on the ...
2
votes
2answers
1k views

Conditional combining of multiple rasters to produce one output using Python?

I have two rasters for same spatial extent and need to compare both rasters cell by cell and produce one output raster. Each has 10 land cover classes. Based on those 10 classes the comparison should ...
8
votes
1answer
4k views

GDAL RasterizeLayer Doesn't Burn All Polygons to Raster

I am trying to burn a shapefile to a raster using GDAL's RasterizeLayer. I pre-create an area of interest raster from a different shapefile, given a specific pixel size. This AOI then serves as a base ...
11
votes
1answer
2k views

Copying file geodatabase on ftp site to local disk using Python?

There is a file geodatabase sitting on an ftp site that I would like to download with a Python script. Right now I'm thinking one way to do this is to copy the ftp geodatabase to a geodatabase on my ...
5
votes
3answers
2k views

Add 3rd party repositories button not in version 1.8.0 Lisboa

The available manuals and other help sources all refer to the 'Add 3rd Party Repositories button' in Plugins>Fetch Python Plugins and then the Repositories Tab. My problem is that the installation of ...
5
votes
1answer
2k views

Export each layer in map to a separate image

I have a mxd file with multiple layers I have to export to PNG/JPG. I would like to create a script which turns on layers one at the time before they are exported. So far I have come up with the ...
9
votes
1answer
2k views

Importing GDAL from Python window of ArcGIS for Desktop?

I am trying to run some Python code using the Python window in ArcGIS 10.1 that uses both the arcpy and gdal modules. However, when I try and import the gdal module I get an error: ImportError: No ...
5
votes
2answers
4k views

How to create reports in ArcGIS using Python?

Does anyone have experience in creating reports from your ArcGIS data (geodatabase, map exports) using a python reporting engine? In my current project, we are using Jasper Reports in the following ...
0
votes
2answers
703 views

Why does TableToTable give ERROR 000732 when published to ArcGIS Server?

I have a python code that uses ArcPy command "arcpy.TableToTable_conversion" to convert a database table to a dbf file. arcpy.TableToTable_conversion The input table is from a SQL Server database ...
3
votes
1answer
1k views

How to implement signals on add or modify feature?

I've created a test plug-in in QGIS successfully, but have had a hard time implementing a signal on a vector layer event (add or modify feature). I've add this to my initGUI method: QObject.connect(...
10
votes
1answer
948 views

How to visualize azimuthal data with uncertainties?

I am trying to make a figure showing azimuthal data with a different range of uncertainties at each point. This oldschool figure from a 1991 paper captures the "bowtie plot" idea that I'm aiming for: ...
1
vote
1answer
1k views

installing module in embeded python in qgis

I'm newbie to Quantum GIS and Python. I'm using QGIS in windows platform. How I can install module for python embedded in QGIS? Also, how can I run embedded python from windows directly? Thanks
5
votes
1answer
783 views

raster algebra in python with rasters of different extents

I am trying to find out how to use GDAL & numpy modules to average a set of rasters....which are Sigma0 values from satellite passes at different times within a year. Each raster has been ...
7
votes
1answer
621 views

Arcpy crashes in loop geoprocessing (general behaviour)

I have several standalone scripts that create a list of shapefiles to process (or coverages in one of the scripts). Sometimes my script crashes after it processes three files, sometimes ten, ...
2
votes
2answers
840 views

Script to extract UTM Zone using Range and Meridian information

I am looking for a python or php script that will accept the Range and Meridian values of a location and use it to identify the UTM zone. I have a table that comprises of Easting, Northing, Range and ...
2
votes
3answers
2k views

An error has occured while executing Python code while doing home range analysis (Qgis 1.8.0)

It keeps appearing: Traceback (most recent call last): File "C:/Users/Worten/.qgis//python/plugins\HomeRange_plugin\HRplugin_dialog.py", line 80, in QObject.connect(self.buttonBox, SIGNAL("...
2
votes
2answers
520 views

Using osgeo modules in ArcGIS10.1 and QGIS 1.8 python shells

This is related to an old question of mine - Issues with osgeo - qgis.core and arcpy modules Now in new versions of all the software QGIS Lisboa (1.8), ArcGIS 10.1 and pyScripter 2.5.3. They all use ...
5
votes
2answers
2k views

arcpy.UpdateCursor for multiple fields

I've hit a wall trying to loop through multiple fields with arcpy.UpdateCursor. I have a list of files (2) with multiple fields (10-15'ish) that have its own error check in SQL, and an error value to ...
3
votes
1answer
425 views

ArcMap 10.1 Python split inconsistent string

I converted a kml file to a feature class. All of the fields from the kml were combined into two fields in the feature class. One of them contains information delimited by "/" and the other contains ...
2
votes
1answer
941 views

Write spatial operations to PostGIS from QGIS or ArcGIS

I need to write my spatial operations to a (new) postgis layer (as the data is to large for a shp file). Can this be done in the QGIS browser or from ArcGIS without SDE? When you go (for example) --> ...
6
votes
2answers
4k views

How can i obtain the data of one selected feature in QGIS using python?

I have a problem / question regarding python and QGIS. Is it, and yes how, possible to get the data of one selected feature in QGIS using python? In my example I have a layer which shows some parcels ...
10
votes
1answer
4k views

When iterating over a vector layer with pyqgis, how do I check whether a feature is selected?

While iterating over a vector layer using the following code (summerized from the example in the pyqgis cookbook), is there a way to check whether a feature is selected? provider = vlayer....
1
vote
1answer
325 views

How to store AutoCAD data in ArcSDE database?

Does anyone have an idea how i can create a database instance on my ArcSDE geodatabase to store my .dwg files? I have over 1 million of these files and would be happy to automatically extract and ...
2
votes
3answers
359 views

Simple script in ArcPy to analyze building proximity to restaurants

I am trying to make a script that will go through buildings in a shapefile row-by-row and do a count of points within .25, .5,...1 mile of it. The other points are restaurants. import arcpy arcpy....
8
votes
2answers
2k views

How to add a toolbar via python?

Through tutorials I learned how to add a toolbutton to the plugins-toolbar via python. Now I wonder how to add a complete toolbar with toolbarbuttons via python. Can anybody give some example-code? ...
3
votes
1answer
839 views

Calculate bounding box + buffer for rendering a map tile with mapnik library

Mapnik library has a Projection object (I'm using Python bindings): spherical_mercator = '+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +...
3
votes
1answer
984 views

GeoDjango Changing My Coordinates

I'm running into a problem when saving a coordinate into one of my models in (Geo)Django. My point field is setup as: location = PointField(blank=True, null = True) So while running some tests in ...
2
votes
0answers
65 views

How do I make a path cross the most possible classes in a raster?

I want to create a function and possibly use Arcmap tools to create something similar to Cost Path. Instead of just choosing the LCP, it must cross through the most different integers. Let's say I ...
3
votes
4answers
8k views

Trying to write a python script to crawl a directory and list all datasources

I'm editing a scrip from ESRI to try and create a list of all datasources within MXD's. I've got the script to crawl the file tree and find all the info, but I would like it to output the data for ...
11
votes
3answers
4k views

How to smooth / “reinterpolate” a raster with GDAL?

Short: I want to change the resolution of a raster and smooth the gray ramp like shown in the images bellow. The preference is to use GDAL, PIL or Numpy. Description: I'm kriging points into rasters ...
8
votes
2answers
816 views

How to catch PostgreSQL/PostGIS errors in GDAL/OGR Python bindings?

Is it possible to catch PostgreSQL output error messages in OGR after a query has failed? For example conn.ExecuteSQL('SELECT * FROM non_existing_table;') In Postgres this returns an error message:...
11
votes
4answers
5k views

Programatically splitting polygon with line?

In looking at this question, it occurred to me that I am not sure how one would go about using geoprocessing tools in ArcGIS to split a polygon with a line programatically. Manually, you would use ...