title: Have Atlas sorting and filtering controls been removed in 2.4?

text: In this diagram ![enter image description here](http://i.stack.imgur.com/RzOe3.png) the atlas generation tab of the print composer has several controls that I don't see in my version (2.4) of QGIS. In particular the feature sorting and feature filtering, both of which sound like they might solve a problem I am having with Atlas generation see :How do I control rendering order of items in composer?. Is this a screen shot from an older version and these controls have been removed?


title: Extremely slow processing time to map lat/lon data onto census shapefile with ArcPy?

text: I am a newcomer to ArcGIS and Python so I would greatly appreciate any and all input on this issue. What I am trying to do is very simple: map 400K business locations using lat/lon coordinates onto a shapefile of California census tracts, and then export the data from the spatial join into a text file so that I have a spreadsheet containing the census tract that each business falls into. The input file of lat/lon data is a csv with 400K records, 5 variables and is ~15MB. I have successfully completed this task interactively in ArcMap 10.2.1, but I would like to have a record of my process for future replication, hence my foray into ArcPy. In ArcMap, this process takes a few minutes. However, my ArcPy script takes multiple hours to complete. It appears to be the CopyFeatures step that requires significant computer resources and so I am not including ArcMap details or ArcPy code beyond this step. Note also that I have run my full ArcPy script on a subset of my data (the first 20 records of my 15MB file) and it works fine and finishes in just a few seconds. Given that, I know the issue is likely with the size of my file, but I have not been able to determine or successfully implement a work around to decrease the memory usage of my code. Are there any ways I can modify my script to decrease processing time for this simple task? I'm hopeful there's a way to make this as efficient with ArcPy as it is with ArcMap. Please let me know if you have any ideas of how to do this. **ArcMap steps:**   1. Import shapefile containing census tracts     * Add data -> choose shapefile   2. Import and prepare lat/lon data for spatial join     * File -> Add Data -> Add XY Data -> select my 15GB csv -> X Field (longitude) -> Y Field (latitude) -> Coordinate System (WGS 1984)     * Right click on XY data -> Data -> Export data -> All Features -> Output feature class is a .shp saved in my GIS workspace -> upon completion of this step, a pop up displays asking whether want to add to map as layer and I choose Yes (don't think this is significant however) **ArcPy script:**               # Import modules     import arcpy, csv     import os     from arcpy import env          # Set environment settings     env.workspace = "C:/Users/ccwright/Desktop/GIS_workspace"     arcpy.env.overwriteOutput = True          # Set local variables:     in_table = "infousa_latlon_small.csv"     x_coords = "longitude"     y_coords = "latitude"     out_layer = "lat_lon_layer"     feature_class = "lat_lon_feature_class.shp"     census_tracts_shp = "census_tracts_CA_only/US_tract_2000_CA_only.shp"     spatial_join = "spatial_join.shp"     spatial_join_output_txt = "spatial_join_data.txt"          # Set coordinate system to WGS-1984:     WKID = 4326      sr = arcpy.SpatialReference()     sr.factoryCode = WKID     sr.create()     env.outputCoordinateSystem = sr          # Process: Make XY Event Layer     arcpy.MakeXYEventLayer_management(in_table, x_coords, y_coords, out_layer, sr, "")          # Process: Copy Features     arcpy.CopyFeatures_management(out_layer, feature_class, "", "0", "0", "0")


title: Exporting USGS contour line shapefile to DWG or DXF?

text: How do I assign the Z value based on the Elevation attribute using ArcGIS 10.1 for Desktop? I've been messing with this for an hour or two now and every time the exported DWG has ZERO for z value for all lines.


title: geocoding addresses specifically into NAD1983 or converting X,Y into 1983 after geocoding

text: I'm going bananas trying to geocode addresses into a NAD1983. I have ArcMap 10.1 and I'm using the default address locator (North American 10.0 from arcgis online). No matter how low I set the sensitivity, it gives me "There was an error trying to process this table" when I specifically choose NAD1983. If I allow it to create the X,Y coordinates in WGS1984, then I get coordinates...BUT my other coordinates (I'm going to measure distance/proximity) are in NAD1983. Any suggestions? I can't figure out how to change them from WGS1984 to NAD1983 when they're in a layer, although I have changed the data frame successfully it seems... Any guidance is very appreciated!


title: qgis printcomposer: can I wrap text in attr table?

text: I have some longish text attributes and when I try and put these on the map using the attribute table they get truncated. Is there a way to get these long lines wrapped.


title: bulk upload shp to postgis

text: > **Possible Duplicate:**   >  Bulk load multiple shapefiles into PostGIS Is it possible to bulk upload shp files into postgis. Im currently playing around with postgis, qgis and geoserver (as a potential arcgis server, arcmap and arcsde replacements) and want i want to be able to do is run a batch file (I guess) that will upload a number of shp files (exported from arcsde) that overwrite existing shp/tables files in postgis. In effect postgis would be a slave database (for the time being). Med


title: Calculating mean slope: Harmonic or Arithmetic Mean?

text: I have to calculate the average upslope percent slope for a large dataset, the basic method is detailed here. However, I've begun to wonder if harmonic mean might be more appropriate than the standard arithmetic mean, since it is technically a rate of change. I haven't seen this turn up in any of the other discussions on averaging slope over points, areas, lines, etc. It should be fairly straightforward to accomplish. edit: The purpose of calculating average slope in this case is to generate one parameter (of many) to be used in modeling channel initiation thresholds. I have a set of field-collected channel head locations that I will collect the flow accumulation, various average upslope parameters, etc, and will use multiple linear regression to try to describe accumulation thresholds in terms of the other parameters.


title: Reference Error: selectZIP is not defined

text: In Firebug I keep getting this error when I try to click a button I made.               <button type="button" onclick="selectZIP()">Execute</button>      I declared a function called selectZIP in a .js file and put it in my HTML file. Yet am still having this error. SelectZIP is a geoprocessor which tries to...select ZIP codes.               function selectZIP(){         gp = new esri.tasks.Geoprocessor("******/arcgis/rest/services/EDDM/selectZIP/GPServer/selectZIP");         myMap.selectZIP.setOutputSpatialReference({             "wkid": 102100         });      In my .js file I have this, before my button statement               <script src="js/EDDM_API - Copy.js"></script>      I don't know where to go from here.


title: OGR VRT MSSQL example

text: I'm trying to import vector data into SAGA GIS from a SQLExpress 2008 server. Firstly is this possible? I would have thought using an OGR VRT layer would have been the way to do it, I've attached a sample of my vrt file, any help would be great.               <OGRVRTDataSource>     <OGRVRTLayer name="tempGeom">         <SrcDataSource>MSSQL:server=(local)\SQLEXPRESS;database=test;trusted_connection=yes</SrcDataSource>          <SrcSQL>select *, geometry::STGeomFromText('POINT (' + CAST(x AS VARCHAR) + ' ' + CAST(z AS VARCHAR) +' '+ CAST(y AS VARCHAR) +')',0) as geom  from dbo.testTable</SrcSQL>     </OGRVRTLayer>


title: How to convert lines to polygons using NetTopologySuite?

text: I am trying to convert lines to polygons using NetTopologySuite(a C# port of Java Topology Suite JTS). I ran the sample console code succesfully and it generated two polylines. I then tried to run the code using three line strings below which are noded at endpoints.               (0 0,10 0)     (10 0,10 10)     (10 10,0 0)      but this move is bringing the demons out of the code. A bunch of errors are popping up as listed below.               C:\Users\Home\Desktop\polgon\polcreate\polcreate\bin\Release>polcreate.exe        at NetTopologySuite.Geometries.LinearRing`1.ValidateConstruction() in E:\Programming     \JTS\JTS2.11\v2.11_cleanup\NetTopologySuite\Geometries\LinearRing.cs:line 119        at NetTopologySuite.Geometries.LinearRing`1..ctor(ICoordinateSequence`1 coordinates,      IGeometryFactory`1 factory) in E:\Programming\JTS\JTS2.11\v2.11_cleanup\NetTopologySui     te\Geometries\LinearRing.cs:line 43        at NetTopologySuite.Geometries.GeometryFactory`1.CreateLinearRing(ICoordinateSequenc     e`1 coordinates) in E:\Programming\JTS\JTS2.11\v2.11_cleanup\NetTopologySuite\Geometrie     s\GeometryFactory.cs:line 552        at NetTopologySuite.Geometries.GeometryFactory`1.CreateLinearRing(IEnumerable`1 coor     dinates) in E:\Programming\JTS\JTS2.11\v2.11_cleanup\NetTopologySuite\Geometries\Geomet     ryFactory.cs:line 535          and so on....      The full code is shown below. Can anyone guide me the right way?               #define simple     using System.Text;     using System;     using System.Collections.Generic;     using GeoAPI.Geometries;     using GeoAPI.IO.WellKnownText;     using NetTopologySuite.Geometries;     using NetTopologySuite.Operation.Polygonize;          #if simple     using NUnit.Framework;     using Coord = NetTopologySuite.Coordinates.Coordinate;     using CoordSeqFac = NetTopologySuite.Coordinates.CoordinateSequenceFactory;     #else     using Coord = NetTopologySuite.Coordinates.BufferedCoordinate;     using CoordSeqFac = NetTopologySuite.Coordinates.BufferedCoordinateSequenceFactory;     #endif          namespace polcreate     {         /// <summary>           /// Example of using Polygonizer class to polygonize a set of fully noded linestrings.         /// </summary>           [TestFixture]         public class PolygonizeExample         {             [STAThread]             public static void Main(String[] args)             {                 PolygonizeExample test = new PolygonizeExample();                 try                 {                     test.Run();                 }                 catch (Exception ex)                 {                     Console.WriteLine(ex.StackTrace);                 }             }                  [Test]             public void Run()             {                 IGeometryFactory<Coord> geoFactory =                     new GeometryFactory<Coord>(new CoordSeqFac());                 WktReader<Coord> rdr                     = new WktReader<Coord>(geoFactory, null);                 List<IGeometry<Coord>> lines                     = new List<IGeometry<Coord>>();                      // old values which ran successfully                 //lines.Add(rdr.Read("LINESTRING (0 0 , 10 10)"));                 //lines.Add(rdr.Read("LINESTRING (185 221, 100 100)")); //dangling edge                 //lines.Add(rdr.Read("LINESTRING (185 221, 88 275, 180 316)"));                 //lines.Add(rdr.Read("LINESTRING (185 221, 292 281, 180 316)"));                 //lines.Add(rdr.Read("LINESTRING (189 98, 83 187, 185 221)"));                 //lines.Add(rdr.Read("LINESTRING (189 98, 325 168, 185 221)"));                           //i changed the values here..                 lines.Add(rdr.Read("LINESTRING (0 0 , 10 0)"));                 lines.Add(rdr.Read("LINESTRING (10 0, 10 10)"));                 lines.Add(rdr.Read("LINESTRING (10 10, 0 0)"));                      Polygonizer<Coord> polygonizer                     = new Polygonizer<Coord>();                 polygonizer.Add(lines);                      IList<IPolygon<Coord>> polys = polygonizer.Polygons;                      Console.WriteLine("Polygons formed (" + polys.Count + "):");                 foreach (object obj in polys)                 {                     Console.WriteLine(obj);                 }                 Console.ReadLine();             }         }     }


title: Parse XML files in Python (ElementTree)

text:
I am trying to parse some xml and rss feeds to extract some of their data in order to store it in a PostGIS database. The file I want to parse is here: http://earthquake.usgs.gov/earthquakes/catalogs/1hour-M1.xml and looks like this:               <?xml version="1.0"?>     <feed xmlns="http://www.w3.org/2005/Atom" xmlns:georss="http://www.georss.org/georss">       <updated>2012-03-05T19:57:55Z</updated>       <title>USGS M 1+ Earthquakes</title>       <subtitle>Real-time, worldwide earthquake list for the past hour</subtitle>       <link rel="self" href="http://earthquake.usgs.gov/earthquakes/catalogs/1hour-M1.xml"/>       <link href="http://earthquake.usgs.gov/earthquakes/"/>       <author><name>U.S. Geological Survey</name></author>       <id>http://earthquake.usgs.gov/</id>       <icon>/favicon.ico</icon>       <entry><id>urn:earthquake-usgs-gov:ak:10425830</id><title>M 2.5, Alaska Peninsula</title><updated>2012-03-05T19:37:00Z</updated><link rel="alternate" type="text/html" href="http://earthquake.usgs.gov/earthquakes/recenteqsww/Quakes/ak10425830.php"/><![CDATA[
Monday, March  5, 2012 19:37:00 UTC
Monday, March  5, 2012 10:37:00 AM at epicenter

Depth: 185.70 km (115.39 mi)
]]><georss:point>58.5578 -155.6727</georss:point><georss:elev>-185700</georss:elev><category label="Age" term="Past hour"/></entry>       <entry><id>urn:earthquake-usgs-gov:nc:71742560</id><title>M 1.5, San Francisco Bay area, California</title><updated>2012-03-05T19:31:16Z</updated><link rel="alternate" type="text/html" href="http://earthquake.usgs.gov/earthquakes/recenteqsus/Quakes/nc71742560.php"/><![CDATA[
Monday, March  5, 2012 19:31:16 UTC
Monday, March  5, 2012 11:31:16 AM at epicenter

Depth: 0.20 km (0.12 mi)
]]><georss:point>37.9882 -122.4550</georss:point><georss:elev>-200</georss:elev><category label="Age" term="Past hour"/></entry>       <entry><id>urn:earthquake-usgs-gov:ak:10425819</id><title>M 2.2, Central Alaska</title><updated>2012-03-05T19:28:00Z</updated><link rel="alternate" type="text/html" href="http://earthquake.usgs.gov/earthquakes/recenteqsus/Quakes/ak10425819.php"/><![CDATA[
Monday, March  5, 2012 19:28:00 UTC
Monday, March  5, 2012 10:28:00 AM at epicenter

Depth: 112.30 km (69.78 mi)
]]><georss:point>63.2167 -150.5241</georss:point><georss:elev>-112300</georss:elev><category label="Age" term="Past hour"/></entry>       <entry><id>urn:earthquake-usgs-gov:nc:71742550</id><title>M 1.8, Northern California</title><updated>2012-03-05T19:18:22Z</updated><link rel="alternate" type="text/html" href="http://earthquake.usgs.gov/earthquakes/recenteqsus/Quakes/nc71742550.php"/><![CDATA[
Monday, March  5, 2012 19:18:22 UTC
Monday, March  5, 2012 11:18:22 AM at epicenter

Depth: 2.40 km (1.49 mi)
]]><georss:point>38.8177 -122.8205</georss:point><georss:elev>-2400</georss:elev><category label="Age" term="Past hour"/></entry>       <entry><id>urn:earthquake-usgs-gov:ak:10425806</id><title>M 2.1, Southern Alaska</title><updated>2012-03-05T19:14:44Z</updated><link rel="alternate" type="text/html" href="http://earthquake.usgs.gov/earthquakes/recenteqsus/Quakes/ak10425806.php"/><![CDATA[
Monday, March  5, 2012 19:14:44 UTC
Monday, March  5, 2012 10:14:44 AM at epicenter

Depth: 17.10 km (10.63 mi)
]]><georss:point>60.5011 -145.1175</georss:point><georss:elev>-17100</georss:elev><category label="Age" term="Past hour"/></entry>     </feed>      I am quite new to Python and to XML but I believe that the right direction to choose is ElementTree. Thus I have started with the following code:               #-*- coding: utf-8 -*-          import os     import urllib     import xml.etree.ElementTree as ET          def main():       feed = urllib.urlopen("http://earthquake.usgs.gov/earthquakes/catalogs/1hour-M1.xml")            try:         tree = ET.parse(feed)         print "Download ok"         root = tree.getroot()         print root         event = root.find("entry")         for e in event:           print e.attrib       except Exception, inst:         print "Unexpected error opening %s: %s" % (tree, inst)          if __name__ == "__main__":       main()      but an error is thrown... Can anyone point me the best direction to follow for this? I do not really understand how to extract the data from the various tags. And then, what's the ideal strategy to "store" this data: an array of dictionaries? Thanks in advance!


title: Testing ArcMAP 10.0

text: I want to test performance of Arc Desktop based applications. I tried GDBT but it not give me any information about request time etc. I tried also a Workflow Profiler, but this tool usually crash my ArcMap instalation. I think about use of arcPy but I am not sure how to do that properly. Any Ideas? Maybe someone of you made this kind of tests. For Manual tests I can use AutoIT but I am not sure what use for loading tests of Arc desktop. And how to get any queries which are send by arc map to Database?


title: How to access an wms layer that requires authorization via OpenLayers?

text: Is there a way to send username and password parameter to a wms layer that needs authorization? If not, I guess it's perhaps a non-spatial thing to solve so that you before the javascript sends a wms request you initiate a httprequest to the server where the wms service is hosted. If anyone have an good example of how to do that from java, html or javascript I would really appreciate it. If you use firefox you could add http:\username:password@wmsurl and its solved and the information is the stored in the session/cache ... Thanks for the help Mike


title: Is there any tool/plugin for QGIS to easily select feature and edit attributes

text: I need to select polygons and change some attributes in a shapefile. I'm using Quantum GIS - I select the polygon, open attribute table and edit the attributes I want. I know there must be an easier and faster way to do this. I was thinking in some tool which when the feature is selectd pop-up a form with the respective attributes and lets you edit them.


title: How to change color of a feature in openlayers?

text: I'm loading a geojson file into openlayers to show some polygons. From this file I also create a list of all names of these polygons. Now I want to change the color of a polygon when selected from the list (clicked on the name). What I've tried is to create a style (http://docs.openlayers.org/library/feature_styling.html) but I couldn't find out how to add this style to the polygon. How can I do that?


title: ArcGIS 10.0 - How to have GP Service make rest calls instead of Silverlight client?

text: I have a Silverlight client that densifies a user drawn polyline. The returned vertices from the polyline are then sent to esri elevation sample server for each vertex to find an elevation value. The user defines how dense the line should be. My question is how could I have my python GP Service make the rest calls to the ESRI elevation sample server once it receives, from the silverlight client, all the densified vertices? EDIT: Would there have to be a python client or perhaps have a slave silverlight client make the calls and return them to python and then back to the master silverlight client. The master silverlight client would message to the user that their variables requires this to be a long running process. Makes no sense? Thanks.


title: Making Multipart Features to Unique features

text: I have a layer with two lines crossing over each other. I create intersection point over them and then create a buffer for a specified distance and then using difference tool, I delete the buffer area. The resultant of the above process is four individual lines. Although the lines are four, the features are only two since they have been just separated by the above process. What is the way to make the feature of each one unique so that I get a total distinct feature of four?


title: What are key things to test in ArcMap 10.1, service pack 1?

text: I have been asked to test service pack 1 for ArcMap 10.1 before it is provided to other users at my university. Are there any key things I should do for any SP test, or should I try to look at things that have specifically been addressed in this SP? My plan was to create a new .mxd file and run a few basic geoprocessing commands. Realistically I'm not going to able to check everything, but are there other key things I should do?


title: How to reduce all value above certain value using Raster Calculator?

text: I have a raster with elevations above 2000m that I want to "cut" to 1800m. I cannot quite figure out the correct string and most simple way to do this in the ArcGIS 9.3 raster calculator. I tried to ceiling function without much luck; would an expression similar to IF([my_raster]>2000,1800,[my_raster]) do the trick?


title: Create multiple tables from distinct attribute category values using Model Builder

text: ![parcel map with book values](http://i.stack.imgur.com/8i5Q3.jpg)I've got a shapefile I would like to extract distinct values from an attribute field and have the output as tables based on those values. The shapefile is made of polygons representing parcels and each parcel is categorized by a "book" value. I can extract a table for each book using the Table Select tool in Model Builder, however, I need to produce a table for each of the dozens of "book" values that make up this shapefile. This is basically a "select by attribute" problem but I will need to sort and produce output in the form of tables. I've attached a pdf of the Layer view and contents to help illustrate the problem. My question is, what sql expression can I use in the Model Builder Table Select tool that will:   1. Select and sort values based on an attribute field (in this case, by "book")   2. Auto create tables for each category of values (i.e., 200, 209, 235) I'm open to other model builder tools or methods but the "Table Select" method is the most direct approach, however, I can only create one table at a time - I need to create many in one shot. Any help is greatly appreciated. Galindo37


title: Programmatically Load Hyperlink File in ArcMap 10.2

text: I am using an arcpy script tool to create an XY Event Layer from a text file and then save the layer to a .lyr. This works, but I also want to load a hyperlink file and apply it to the layer automatically. I tried using UpdateLayer, but this just replaces the entire layer, it doesn't keep the same source data. So, does anyone know if there is a way to load a .hlk file programmatically? The solution does not have to be in Python.


title: GDAL - How to get gemetry type of a MapInfo TAB file?

text: I am using the OGR MapInfo File driver to read a MapInfo TAB file. However, the function GetGeomType() returns 0 which means 'wkbUnknown` even if the features in the TAB are only multi-polygons. Is it because MapInfo TAB can store many different GeomTypes in the same file? If I take for granted that only one geometry type would always be in the TAB, how could I retrieve it?                       driver = ogr.GetDriverByName("MapInfo File")             datasource = driver.Open(os.path.join(dirname,shapefileName))             layer = datasource.GetLayer(0)             geometryType = layer.GetGeomType()                  >>>geometryType             >>>0


title: QGIS python location problem

text: It seems that MacPorts installs Python27 in /opt/local, but that MacPorts QGIS updates and uses the Apple Python27 in /Library. The result is that MacPorts additions to Python27 such as OpenCV don't work because they are added to the /opt/local Python27 whereas with QGIS installed all Python links point to the /Library Python27 so things like OpenCV can't be used. Has anyone found a way of keeping QGIS working while either 1. making it use the MacPorts Python27 at /opt/local, or 2. enabling MacPorts to add things to the /Library Python27? The MacPorts people agree it's a problem, but are focussed on finding a 100% MacPorts solution (none yet after 5 months); they aren't interested in cooperating with the native QGIS or Apple Python. I've no response from the OpenCV people as to whether it can be linked to the /Library Python27. (Not being a Unix guru, I'd prefer to stay with MacPorts if that's possible.)


title: ArcGIS Server 10 not publishing map services

text: In the past week, I managed to set up ArcGIS Server on the server but subsequently encountered errors when publishing a web map in ArcMap. I have set up the ArcGISSOC and ArcGISSOM accounts and given rights to the specific folders. However I encountered the following problems that I have attached screen shots and highlighted in RED. The error reads, `"The SOC account cannot resolve the path to the data, or the SOC account has not been granted read permission to the data. In the case of ARCSDE, direct connect may not be configured properly or the SOC account may not have permission to access the database."` I also noticed under the (ArcGISSOM Inheritance) screen shot, under the inherited from tab, the status reads as `<not inherited>`, what could be the cause of this because from my understanding, the webmapping folder should have rights from the server drive. I also tried publishing from the ArcGIS server manager but it was unsuccessful. Could you kindly assist in solving the problems that I am encountering so that I can effectively do the web mapping.


title: How to change feature class and field aliases in bulk?

text: I have over a hundred FCs, each with 10 or 20 attributes to add or change the aliases for, two or more times a year. Needless to say, this is not something I'm going to grunt my way through. How can I automate this process? Python solution preferred but will use anything that works. I have access to Arcgis 9.3.1 and 10 (ArcInfo license level).


title: Street block and Census Block

text: Study Area: Manhattan Software: ArcGIS 10.0 I'm measuring the street density using Street Network Grain size measure = Number of Street blocks per square mile. I have the NY street network but I do not have the map of Blocks. Is there any place I can get the Street Block maps? Alternatively - Is the census block equivalent to a street block in Manhattan Borough. I have the census block polygon shape file from TigerLine. Thanks for the inputs!


title: model builder iteration

text: I am trying to develop a model in modelbuilder that will allow for a number of polygon shapefiles to be paired with each other. Say a group of 3 shapefiles (1,2,3) located in the same folder would be processed producing individual shapefiles 1-2, 1-3, 2-3. Is there a way to do this in modelbuilder? Thanks in advance


title: How to check if two qgsPoints are equal in PyQgis

text: I have two points:               a=QgsGeometry.fromPoint(QgsPoint(227739,908343))     b=QgsGeometry.fromPoint(QgsPoint(227739,908343))      but when I run:               a==b      it returns False. What am I doing wrong? **Update:** it returns True when I only use:               a=(QgsPoint(227739,908343))     b=(QgsPoint(227739,908343))


title: Am I over-complicating adding and then calculating fields? Suggestions for a better way? (Python)

text: I'm attempting to re-position a question I asked yesterday here. I am starting to think that the way I wrote my code just isn't the best way. Am I overcomplicating things? I know for a fact that adding a bunch of new fields and then calculating those fields based on the new fields is a common occurrence in GIS, so I just can't understand why it's so difficult to find out how to do it in a Python script. The following code is a snipet of what I came up with, but despite hours and hours of troubleshooting, I keep getting an error (although it did work on a small sample size). I'll paste the code first, and then the error below. I'm receiving the error just after all of the fields are successfully added. I've done a number of things in an effort to troubleshoot, such as import all of my files into a geodatabase and run it from there; I've tested small samples of data (which does actually work, so I know the code "works"), and I've switched out fields I thought were problematic with other fields that were not. So far, I'm unable to run this code on all 108 fields that I'm trying to add and then calculate. Again, I receive an error just before it attempts to calculate the fields, and it seems to be referring to the line: `value = row.getValue(field)* row.anotherField` Any suggestions for a better way to do this?? Thanks! Code snipet:               ....          fc = arcpy.CopyFeatures_management(layerName, outFeatures)          #---------------------------------------------     # Add fields to a feature class:     #---------------------------------------------     fieldList = ["field1", "field2", "field3", ............., "field99"]  #<-- fieldList     fieldListx = []     for f in fieldList:         fieldx = f+"x"         arcpy.AddField_management(fc, fieldx, "DOUBLE")         fieldListx.append(fieldx)          #---------------------------------------------     # Calculate fields:     #---------------------------------------------     rows = arcpy.UpdateCursor(fc)     for row in rows:         for fieldx in fieldListx:             field = fieldx.rstrip('x')             value = row.getValue(field)* row.anotherField  #<--This is the line the error is referring to             row.setValue(fieldx,value)         rows.updateRow(row)     del row, rows      Error:                Traceback (most recent call last):        File "C:\GIS\1GIS_DATA\CUR_GIS\Election_Districts\Python and Tools\Join_Voter_Results_to_Allocation_Layer_and_Add_New_Fields.py", line 125, in <module>          value = row.getValue(field)* row.UNpctOfEd #<--GIVING ME SHIT WHEN I ADD THE LARGER LIST!!!        File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\arcobjects\arcobjects.py", line 944, in getValue          return convertArcObjectToPythonObject(self._arc_object.GetValue(*gp_fixargs(args)))      RuntimeError: ERROR 999999: Error executing function.


title: How do I add an offset to a lat/lon to get a new point using PostGIS

text: Using PostGIS I have a lat/lon point and want to add an east/west and north/south offset in metres to get a new point. Is that possible?


title: Are there any .Net opensource library for geometry operations?

text: Are there any opensource librarys for doing geometry operations (union, intersect, relation...) in .Net?


title: How to create a time series with a NetCDF raster from Ensembles E-OBS using GDAL, R, GRASS GIS or SAGA GIS

text: I downloaded a NetCDF file that contains daily precipitation values from 1-1-1950 to 30-6-2013. The raster has 23191 bands. http://eca.knmi.nl/download/ensembles/download.php My goal is to analyze how precipitation values varied and relate the changes with the conservation of a forest ecosystem. I'm unsure on how to proceed, but I think the first step should be to create a time series from the NetCDF raster. How can I do this using GDAL, GRASS GIS or R?


title: Any way to use a local background map in ArcGIS Mobile

text: My customer requirement is to support a MrSID image directly in ArcGIS Mobile. The details are sketchy but I think the intent is that they want a local background raster Map while the ESRI.ArcGIS.Mobile.Map gets used to handle other tasks (?vector display, updating GDB?). The Map control seems designed to consume Map Services via the MapCache object. Is there any way to get the MrSID image to render, pan and zoom in ArcGIS Mobile? A "How to" reference is perhaps too much to hope for. But I'll appreciate any related suggestions. For example, how might we use a local TIFF or JPEG image as the background map?


title: How to use tiles in this strange tiling scheme?

text: I'm looking at cell coverage maps such as this one, and trying to figure out how to get a GIS program to use it as a source of information. Ultimately, I need to be able to combine the output with other cell coverage maps from different sources. But in order to do that, I need to understand how each one works and how to align them with each other. Strangely, each of the five zoom levels offered has a different tile size:   1. 160x160 pixels   2. 128x128   3. 145x145   4. 198x198   5. 182x182 I've never seen a map tile service like that before. Usually the tile sizes are all the same. Do map-service-consuming tools (I'm thinking of open-source tools like Tilecloud) have the ability to consume maps whose layers have different - practically random - tile sizes? Additionally, I'd like to know how the tiles are georeferenced. Maybe that's not exactly the right way to say it ... but I need to be able to combine tiles from different maps so that they're aligned properly. Maybe I could get the information needed for that if I knew what kind of map tile API is being implemented by http://gsma.streetmap.co.uk/custom/barts/gsma/. Maybe they have a REST URL that gives data about spatial reference and so on. But I haven't been able to figure out what API they're implementing. Using the browser debugging tools, I see URLs like http://gsma.streetmap.co.uk/custom/barts/gsma/gmap.srf?o=1&i=502/951/Zoom3/Zoom3_502_Tile007_9511638.png for the tiles. The mobileworldlive.com page loads several javascript files, including two from gsma.streetmap.co.uk: http://gsma.streetmap.co.uk/custom/barts/gsma/SMJ_GSMA.srf?id=502&layer=951&nd=http://mobileworldlive.com/maps/images/map_na.jpg and http://gsma.streetmap.co.uk/SMJ_ObsEmo.js (The first one has an "include" that loads the second.) They contain a bunch of custom map display code. I suppose I could try to learn about the API through reading that code, but that seems like a time- consuming and error-prone way to do things if this is implementing an API that's documented somewhere. Thanks for any suggestions on how to attack this.


title: Compare slopes between adjoining line segments

text: Using a DEM, I assigned slope values to individual segments in a street network feature class. In some areas the slope values seem suspect, such as where a highway overpasses a local road. Is there a method to find those segments whose slope is radically different from adjoining segments (e.g. 5 to 25%) to locate possible inaccuaracies?


title: How to clip a raster by a shapefile in pyqgis?

text: Currently I am working with some satellite images (grib files) in pyqgis. As I want to work only with the cells in my study area, I need to clip the raster by a shapefile. I know I can do this by Raster tools in QGIS, but since I have a lot of rasters to work with, I wanted to use python and save time. Perhaps somebody knows the pyqgis commands to do this? I am new to QGIS and python, so I would really appreciate some help.


title: ESRI JSON Polygon Ring Orientation

text: I am having a hard time finding any detailed documentation defining ESRI's JSON specs. I am hoping someone can shed some light on the following two questions.   1. Similar to the GeoJSON spec, the first array in the ' ** _paths_** ' collection is always an exterior ring, and all subsequent arrays in this collection are interior rings (holes)?   2. Is above is true, how does the ESRI JSON spec deal with multi-polygons?


title: Add primary key to all tables in PostgreSQL automatically

text: I have a database with 10 schema and every schema has 10 tables which could be points, lines, etc.). I have created a folder called data and sub-folder with the same name of schema.. I use FME to detect new tables in the folders and transfer data to postgis (data from administrative folder will be written it in administrative schema). Problems; With generic/dynamic workspace, FME write geometry in all allowed geometry in postgis (geometry must be changed manually in postgis from geom type geometry to geom type linestring,., etc. ). The second problem is to have a primary key in postgis for all tables (did not work with FME). thanks Is it possible to create a function/Trigger that can look for all tables (also new tables if I want to add) and add a primary key for each table?


title: How to use Canvas element with Openlayers?

text: I would like to use element canvas of HTML5 in Openlayers to render maps. I am searching some code or reference in which I can check it, can you help me in this? Thanks in advance


title: "EOFError: EOF when reading a line" Using raw_input() in ArcMaps Python Console

text: I am writing a python script to automate some ArcMaps data processing. I am getting hung up at raw_input(). I am using raw_input to get input from a user. A folder will be created with the folder name based on the user input. I don't understand why this isn't working. I've viewed tutorials, read the manual, etc... and I still can't figure out why I am getting this error for raw_input(). I've managed to be able to get user to select files and directories suing the Tkinter toolkit should I look into using this toolkit to get text from a user prompt as well?


title: List of central coordinates for all countries?

text: I need a list of the central points for all countries: China: lat/long (coordinates of the most central point in China)   France: lat/long (coordinates of the most central point in France)   etc...


title: Bounding Boxes around nearby points

text: I've got a table with 200.000 points. My aim is to draw bounding boxes around points where this point have the same name and that are in a maximum distance of 500 meters to each other. I guess I can do with `st_envelope` but I do not know how to write the SQL query. Thanks for your help! Best Regards, Carsten


title: How to perform Noise Pollution Mapping using ArcGIS for Desktop?

text: I've been trying to do a noise pollution map, and I really don't know where to start with the information I have. The end result needs to be like in the link below. I've seen the question, How to create noise maps? (open source or scripted in ArcGIS) but it doesn't give an accurate metod, and I don't know how to go about it. I've searched for weeks and weeks now about a simple method, and I'm none the wiser. Data I have:   * **parameter** - _speed_ (roads)   * _environments_ (hard ground=lakes and soft ground= grassland)   * _sound data_ i.e. dB accounting to number of vehicles What are the tools I can use in ArcGIS for Desktop to achieve a map like the link above? I only have access to ArcGIS for Desktop.


title: Converting DEM data to height information in meters using QGIS?

text: I'm relatively new with QGIS and currently I'm working with ASTER and SRTM DEMs. The goal is to create vector heightcontours with labeling in meter. Can you tell me how I can convert the grayband data in height information in m?


title: Calculate coordinates of centroid in arcgis using python

text: I have polygon features and want to calculate coordinates of it's centroids in arcgis using python - all goes fine, but it gives me coordinates in meters and not in degrees as i want. Here is the code:               arcpy.CalculateField_management(folder+"\\"+f+"_fakela_shp.shp", "X", "!SHAPE.CENTROID@DECIMALDEGREES!.split()[0]", "PYTHON")      arcpy.CalculateField_management(folder+"\\"+f+"_fakela_shp.shp", "Y", "!SHAPE.CENTROID@DECIMALDEGREES!.split()[1]", "PYTHON")       **@DECIMALDEGREES** should done the job, but it is broken somehow.


title: get IMap in arcgis 10.2

text: i'm updating a SOE from pre 10.1 arcgis server to arcgis server 10.2.2. The orginial code included the following part, which is not working anymore:               IMapServer3 mapServer = serverObjectHelper.ServerObject as IMapServer3;     IMapServerObjects3 mapServerObjects = mapServer as IMapServerObjects3;     IMap map = mapServerObjects.get_Map(mapServer.DefaultMapName);      I found that pre 10.1, this code would only work for MXD and not for MSD documents. ESRI says this: > All map services now use the fast drawing engine that was associated with > MSDs in versions 10.0 and previous. So it seems logic this piece of code would not work anymore. Anyone has an idea on how I get my IMap variable in 10.2? Thanks


title: syntax error when using gdal in python console QGIS

text: When I run the following code in the python console in QGIS I get a syntax error:               import osgeo.gdal     gdal_rasterize [-b 1] [-burn 255] [-a Pixelkost] [-te 170000.2 168000.27 177999.8 177999.8] <Adp322.shp> <work.tif>      b command for burning it in band number 1 burn command for the color red a command for using the attribute pixelkost te command for x,y min and x,y max the already loaded Adp322.shp file (I also tried it with the whole system path towards the file but this didn't work either) work.tif is the name the new file should have. (Should this maybe be an already existing file? but how should I produce this?)


title: changing attributes or names of shapefiles

text: I am trying to join a subset of two similar shapefiles, one with states from the US and the second with provinces from Canada. The datasets are available here: shapefiles from Geocommons I noticed that in the two files the attributes are slightly different. In the US map, states are labeled as "STATE_NAME" but in the Canadian map the attributes are simply "NAME". This is a problem because I cannot join the two shapefiles into a single one. Does anyone have a work around for this? here is my code so far:               require (raster)          #load in boundaries for plotting      state <- readOGR(dsn = '/usa_state_shapefile.shp', layer = "usa_state_shapefile")     projection(state) <- CRS("+proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs")     # Subset US shapefile by desired states     nestates <- c("Maine", "Vermont", "Massachusetts", "New Hampshire" ,"Connecticut",                   "Rhode Island","New York","Pennsylvania", "New Jersey",                   "Maryland", "Delaware", "Virginia", "West Virginia", "North Carolina")     state.sub <- state[as.character(state@data$STATE_NAME) %in% nestates, ]     summary(state.sub)          provinces<-readOGR (dsn = '/canadian_provinces.shp', layer = "canadian_provinces")     projection(provinces) <- CRS("+proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs")     canprov<- c ("Quebec", "Ontario", "Newfoundland  & Labrador", "New Brunswick", "Prince Edward Island", "Nova Scotia")     provinces.sub <- provinces[as.character(provinces@data$NAME) %in% canprov,]     summary (provinces.sub)      Im guessing that if I rename the attributes to the same thing like "NAME" then I should be able to merge the two shapefiles using some sort of `rbind` or `cbind` function.


title: maximum circle inside an irregular polygon from a random point

text: I am trying to work out an algorithm for creating the maximum radius circle within an irregular polygon (a census tract) based on a given center of the circle. The motivation for this is to obscure the location of an individual who has answered a survey. Their actual location is known, however it needs to be obscured in the analysis, in order to release the data to the public, for further analysis. We want to have a donut shaped polygon for each survey respondent that has some inside radius (easy), bounded by an outside radius that is constrained by the census tract that the individual is in. Their final location will be randomly placed within the donut polygon. I have seen many answers to similar questions here, but not this specific one, which in this case starts with a SPECIFIC location. Once we have the donut established, we can randomize the location of the individual response within the polygon. That's relatively easy... Thanks for your ideas, mine so far have seemed pretty brute force, and computationally "expensive" or inefficient...


title: Type Error with ArcPy Update Cursor

text: filter = 'LinkID IS NOT NULL'      cur1 = arcpy.UpdateCursor("ODDS_Links", filter)           high_id = 0     for row1 in cur1:      if high_id < row1.LinkID:       high_id = row1.LinkID                    filter2 = 'LinkID IS NULL'      cur2 = arcpy.UpdateCursor("ODDS_Links", filter2)          i = high_id     for row2 in cur2:      i+=1      row2.LinkID = str(i)       cur2.updateRow(row2)       The error I am getting is below. I know this works without the counter but something is wrong with counter and it won't add by one. If I comment out the counter it displays the error. I know I should use da but really don't feel like changing things around when I am so close. > TypeError: coercing to Unicode: need string or buffer, int found


title: how to manage xml fragments in fme

text: I am trying to flatten an xml file with fme 2014. I´m having problems understanding how to manage the following xml schema. One of the XML child elements is called "aspects" and contains the following ...               <aspects>     <a id="1" l="de" v="Alter Schafstall und Eichen beim Südsee-Camp"/>      How do I go about flattening the attributes? FME seems to flatten most of the table, but the above gets defined as an "XML Fragment". In the table the field "XML_Fragment" then contains the above xml structure with tags etc...ie is not being flattened.


title: How do I plot data stored in an attribute table (really in a database) in qgis

text: I have a csv of UK postcodes from Ordnance Survey that looks like               SE1 0AA,10,531869,179427,E92000001,E19000003,E18000007,E09000028,E05000536      When I import to QGIS using 'Add Delimited Text Layer', I can declare columns 3 and 4 as x and y and everything plots nicely with a new layer created. I want to store the same data in PostGIS so I can join other data on the postcode and then plot resulting attributes in QGIS. The postcode table, with the same fields, has no geometry and I can't see a way to declare the x and y columns as coordinates, so it is shown in QGIS as just an attribute table. Can anyone explain what I have missed? I apologise for any mangled terminology - this may be why I haven't found any solution on the web yet.


title: Adding a Vertex To Exisiting Feature(linear) Using OGR Python

text: I am working on a project that involves finding the intersection of linears in a shapefile, if that vertex(the found intersection) doesn't already exist for each(or one) of the features, then that vertex(the found intersection point) should be added to both features. I'm using the shapely and OGR modules for python. Can somebody just give me a general example of adding a vertex to an existing feature in shapefile using either or and saving it. Thanks.


title: Using SpatiaLiteLayers in QGIS

text: I've two questions regarding the handling of SpatialLite databases in QGIS.   * I've a SpatiaLite database with a table "Aufschlüsse" and the user defined columns "Name", "Typ" and "Art". If I do load this table into QGIS I do get the column "Name" only (byside the predefined "PKUID"). How can I change the SQL-Statement to get all my columns or e.g. "Name" and "Typ"?   * I've tried to delete the column "Art" using the Database Manager plugin. But I'm unable to delete any of the columns. The "delete column" buttom is inactive. Why? Sincerely, Thomas


title: Trouble using saveACopy with group layers

text: This code works fine on layers in the TOC that are not group layers, but errrors out if there are any group layers in the TOC.               import arcpy, re, os     arcpy.env.overwriteOutput = True          #loop through layers and save a copy          def layer_save(mxd_path,path):              mxd = arcpy.mapping.MapDocument(mxd_path)         df = arcpy.mapping.ListDataFrames(mxd)[0]         fcs = arcpy.mapping.ListLayers(df)         for item in fcs:             formatted = str(item).replace("/","-")             item.saveACopy(str(path) + "\\" + formatted + ".lyr")                 if __name__ == '__main__':              #variables         mxd_p = arcpy.GetParameterAsText(0)         path_lyrs = arcpy.GetParameterAsText(1)         #kickoff         layer_save(mxd_p,path_lyrs)


title: merge 405 Aerial Photos into one big file

text: I have tried to merge my 405 Aerial Photos into one. I have FME and ArcInfo 10.0 at my disposal on Windows 7 64bit. I have tried the following: create a mosaic; export all images from a rastercatalog; create new empty raster and then loaded the rasters, and in FME I have tried the RasterMosaiker. Everything starts ok, but then due to the file size creates an error after many hours. Does anyone have any experience creating such large single raster data files and can offer some tips on how I could achieve this, My only restriction is that I don´t wan´t to decrease resolution. Thanks for any comments, Rob


title: Adding EPSG:3857 / 900913 to Oracle's mdsys.cs_srs

text: I am working with an Oracle Standard database at version 11.2.0.1.0. I also have access to XE version 11.2.0.2.0. The later version has both SRIDs 900913 (Google Mercator) and 3857 (official Web Mercator), but the earlier version does not. I read that SRSs are easily inserted into Oracle's mdsys.cs_srs table, so I selected the values for 900913 from the later DB and inserted into the earlier DB. The insert reported success, but no matter how many times I update it the WKTEXT3D field in 11.2.0.1.0 remains null. If I try a coordinate transformation from a single point (both 2D and 3D) in 4326 (WGS84 Lat/Lon) to 900913 I am told `"Parsing error on source or destination coordinate system WKT, or incompatible coordinate systems"`, then `"Check the validity of the WKT for table entries, and check if the requested transformation is valid"`. I can use the same transformation query (below) in 11.2.0.2.0 without issue. Does anyone know   1. Why the WKT extracted from 11.2.0.2.0 doesn't seem valid in 11.2.0.1.0?   2. How I can adjust it to make it valid?   3. If there is some fundamental issue in 11.2.0.1.0 that makes this transformation invalid? Any input on these issues verrrrrry much appreciated. **Edit** : Following MelitaKennedy's suggestion I started trying to track down the 'Popular Visualisation Pseudo-Mercator' projection in my 11.2.0.2.0 DB and I think I found it in MDSYS.SDO_COORD_OPS. However, it references a COORD_OP_METHOD_ID (MDSYS.SDO_COORD_OP_METHODS table) which defines several 'methods'. I assume that each of these methods relates to a function within Oracle, so maybe it simply isn't possible to create the entire CRS / Projection package in a pre-11.2.0.2.0 DB? Transform Query:                 select         sdo_cs.transform(           sdo_geometry(             2001,             4326,             null,             sdo_elem_info_array(1,1,1),             sdo_ordinate_array(-132.525665, 45.432552)           ),           sdo_dim_array(               sdo_dim_element('LONG', -180.0, 180.0, 0.5),               sdo_dim_element('LAT', -90, 90, 0.5)           ),         900913)       from         dual;      900913 Insert:               insert into mdsys.cs_srs (cs_name, srid, auth_srid, wktext, wktext3d)       values (         'Google Mercator',         900913,         900913,         'PROJCS["Google Mercator", GEOGCS [ "WGS 84", DATUM ["World Geodetic System 1984 (EPSG ID 6326)", SPHEROID ["WGS 84 (EPSG ID 7030)", 6378137.0, 298.257223563]], PRIMEM [ "Greenwich", 0.000000 ], UNIT ["Decimal Degree", 0.0174532925199433]], PROJECTION ["Popular Visualisation Pseudo-Mercator"], PARAMETER ["Latitude_of_origin", 0.0], PARAMETER ["Central_Meridian", 0.0], PARAMETER ["False_Easting", 0.0], PARAMETER ["False_Northing", 0.0], UNIT ["Meter", 1.0]]',         'PROJCS["Google Mercator", GEOGCS["WGS 84", DATUM["World Geodetic System 1984", SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG", "7030"]], TOWGS84[0, 0, 0, 0, 0, 0, 0], AUTHORITY["EPSG", "6326"]], PRIMEM["Greenwich", 0.000000, AUTHORITY["EPSG","8901"]], UNIT["degree (supplier to define representation)", 0.0174532925199433, AUTHORITY["EPSG", "9122"]], AXIS["Lat", NORTH], AXIS["Long", EAST], AUTHORITY["EPSG", "4326"]], PROJECTION ["Popular Visualisation Pseudo-Mercator"], PARAMETER ["Latitude_of_origin", 0.0], PARAMETER ["Central_Meridian", 0.0], PARAMETER ["False_Easting", 0.0], PARAMETER ["False_Northing", 0.0], UNIT["metre", 1.0, AUTHORITY["EPSG", "9001"]], AXIS["X", EAST], AXIS["Y", NORTH], AUTHORITY["", "900913"]]'       );


title: Indexing GIS data associated with user with timestamp

text: I have data for each user in this format user1:(timestamp,lattitude,longitude) user1:(timestamp1,lattitude1,longitude1)..... Here the 1 at the end of timestamp/lat/long is to specify a different point in 3d space. I want to have a indexing structure where I can associate a point in 3d space (spatio-temporal) with a user and be able to do nearest neighbour queries based on either spatial data or temporal data or both. Is there a good indexing scheme for something like this or an algorithm in literature to look at?


title: ArcMap 10.1 Layout toolbar disabled/greyed out

text: For some reason, when I removed and replaced toolbars in an .mxd document the Layout toolbar is greyed out. I've tried opening/closing the document and ArcMap, in addition to readding the toolbar several times, and also tried going into Data Frame Properties and ensured that automatic radio button was checked. Any suggestions for what I'm sure is a simple, yet annoying, problem would be much appreciated!


title: Can you do ST_Transform on geographies?

text: I am brand new to postgis. I just learned that geometries model simple cartesian objects and that geographies model shapes on the surface of the earth. I have a bunch of data in columns X and Y of a table -- where X and Y are state plane coordinates in a medium-sized US city. It seems like you use `ST_Point()` to get the coordinates into a point and then `ST_SetSRID` to get them into a certain coordinate system and then ST_Tranform to get them into a lat/long format. But these are points on the surface of the earth -- so I really want a want a geography not a geometry, right? Am I describing the correct method to do this? Do you do ST_Transform on geometries only?


title: Method to generate points in any projection

text: I am writing an application in ArcMap to generate a grid of points. These points have a very defined layout and I have worked out an algorithm that works great under a projected coordinate system. However, as soon as I try it under a geographic coordinate system, my math doesn't stand up as well and I am not an expert in this type of math (but will to learn if pointed in the right direction) Are there any suggested methods, libraries, or anything else that might make this easier? A little bit more information about the dots themselves: In a circle of 144 foot radius, the dots will be generated in a hex type pattern. The hex is made up of an equilateral triangle with sides of 9 meters (in meters because ArcMap IPoint X,Y properties are in meters).


We are going to let our users add points, polylines and polygons to a map. Are there any international standards for map symbology? For example for emergency symbols?


title: How can I auto-increment output file-names in Toolbox ESRI ArcMap 10

text: I have a simple tool which outputs a file. I would like the filename to be a number which increments every time the tool is run. How would I go about implementing this **using the toolbox in ESRI ArcMap 10?** Thank you for your time


title: How to prevent text labels from overlapping other markers?

text: Is their a way within QGIS to prevent text labels from overlapping icons used to label a point layer? In the example map shown below, there are two layers:   1. Roads layer which is labeled using the new labeling engine.   2. A point layer which uses svg icons to show highway shields. I have other use cases for this such as preventing overlapping of rail station icons. In this example, I could simply move the point, however, in other cases, such as a subway station, that is not a workable solution. A possible solution which I am considering would involve:   1. Adding a column to the point layer's dbf file.   2. Filling the added column with either _Null_ or a space.   3. Enabling QGIS's [new style] labeling for the added column. Is this the way to go or is there a better method? ![Sample Map](http://i.stack.imgur.com/VbV59.png) QGIS info: I am using versions 1.7.4 and a recently compiled version of 1.8 on Mac OS 10.7.4. Linux and Windows versions are also available to me. * * * **Progress Update** : The method I described above works though two modifications to the process are required.   1. Instead of using a space character or a _Null_ , a _non-breaking space_ has to be used. While testing I found, fields with _Null_ are not rendered and fields with only a space(s) are converted to _Null_. (See this Wikipedia page for information on how to type a _non-breaking space_.)   2. While setting the Layer Labeling Settings, on the Advanced tab, set the **Placement** to _Over Point_. Also set the **Priority** such that it is higher than any other labeled layers. Depending on the size of the icon used, some experimentation may be required to get things to work as expected. Depending on the circumstance, it may be easier to use multiple characters rather than attempting to increase the font weight and/or size. ![Sample Map 2](http://i.stack.imgur.com/Qbx0u.png) Tested with QGIS 1.7 and a recently compiled of 1.8 on Mac OS 10.7.4.


title: ArcPy Rename Loop

text: Does anybody have an ArcPy script which points to a directory and renames files within that directory? For example, I have approximately 1500 raster files with names along the lines of m_3609922_nw_14_1_20100425 which need to be subsetted to: "3609922nw". Thanks much! -Aaron


title: PostGIS query with Lon Lat

text: Need SQL syntax for ... I have a table with individual polygons like a jigsaw puzzle. Given a Lon Lat, I need to find the polygon where the lon lat in located. I am trying: select * from "parcels_06003" where "geom" @ ('-120.0077874119654', '38.71168052297362') geom is my polygon column in the db.


title: iPhone app to record 3D location information

text: Does anyone know of an app that would run on an iPhone or iTouch that would record the X,Y and Z information from an external GPS? I am using a Dual XGPS150 Bluetooth receiver. The GPS will record X,Y,Z info, but I am having issues with an app that will record the info, and provide it back to me in a useful format. I have tried GISRoam, but it will not record the Z information. I have looked at other apps, but I am not sure if it will record the Z info and allow me to export it into a format that can be used in a GIS solution like QGIS. Does anyone have any experience with the WolfGIS, Map My Run, or any of the other IOS apps? Thanks...


title: Getting Attachments with token in a map service

text: I have successfully built the functionality to load an image stored as an attachment in silverlight by retrieving the url of the attachment and binding to it: (http://Server/ArcGISServer/rest/services/PublishedMap/MapServer/0/3/attachments/4) However i have implemented token authentication and i was wondering if its still possible to retrieve the attachment. I cant seem to get it now. do i need to embed the token to the url or something? (https://Server/ArcGISServer/rest/services/PublishedMap/MapServer/0/3/attachments/4?token=WCQfGMQZ7SmTProN6fX4SZRBr_0pRFDCvYGzkO2I9BlY2vdsKtQRXoHj8pVoFL_izGuLKiBzd3zZlqFv5U767Q..)


title: How to split polygons based on other polygons?

text: In ArcGIS, I have a layer of polygons (created from the Service Area function) and want to split each one of those based on another polygon layer (the boroughs of London). The purpose is that I want to measure 30-minute accessibility however the density of accessible opportunities various at the borough level. (Dealing with data as far back as 1860s, better data unavailable). I want to measure the area for each polygon (generated by the Service Area function) that lies within each borough. The Intersect function seems to not suffice as it not only splits each "service area" polygon at the borough level, but also by other "service area" polygons. Perhaps another step is required, or a better way to look at the problem altogether. Any help appreciated.


title: Help writing script using a shapefile - select by attribute

text: I am trying to take an existing shapefile and narrow down the attributes using the 'select by attribute' The problem I am trying to solve is: Find parcels that have an acreage >= 100. The code I have is not narrowing down the search or showing up in arcmap. Below is the code I have so far. Any help would be greatly appreciated!               import arcpy          #Get Parcel layer     countyParcels = "C:\\Users\\Krystle\\Desktop\\WCGIS\\Geog485         \\final\\TarrantCoParcels\Parcels2013.shp"     nameField = "TAXPIN"          try:         #Make a feature layer of parcels2013         arcpy.MakeFeatureLayer_management(countyParcels, "AllParcelsLayer")          except:         print "Could not create feature layer"          try:         #Create a seach cursor         rows = arcpy.SearchCursor("AllParcelsLayer")         #Call SearchCursor.next() to read the first row         rows = rows.next()              #Start a loop that will exit when there are no more rows available         while row:             arcpy.SelectLayerByAttribute_management("AllParcelsLayer", '"CALCULATED" >=100.00')             #Print vaules in the current row             print row.getValue(nameField)             #Call SearchCursor.next() to move to the next row             row = rows.next()              #Clean up cursor and row objects         del row         del rows         #Write the selected features to a new featureclass         arcpy.CopyFeatures_management("AllParcelsLayer", "parcels_100plus")          except:         print arcpy.GetMessages()


title: How to draw bar diagrams on the map?

text: I need to draw simple histograms in QGIS version 1.8 from an imported DB. I need a double bar diagrams (eg: male/female - in/out) but I only find the possibility to set pie charts. Is it possible to draw histograms with high proportional to the value?


title: How are service areas created in ArcGIS Network Analyst?

text: I want to create a service area map of Burma in order to spot the areas that are within xx kilometers from a set of cities. The problem is that the generated polygons aren't exactly what I expected. In order to understand why, I tried to generate service areas on a very simple network where the greens dots are my three facilities : ![enter image description here](http://i.stack.imgur.com/a6wVF.png) The first network is a simple triangle. We can see that computing service areas on this shape completely failed. However, if I complexify a little bit the network by adding edges inside the triangle, it runs smoothly : ![enter image description here](http://i.stack.imgur.com/jEoiH.png) In front of these results, I wonder if someone can explain to me how a service area is generated. More precisely, how does ArcMap interpolate the values of the polygons in the areas where there are no streets? Thanks a lot for your help!


title: QGIS python - reference a layer in the table of contents

text: I'm trying to create a plugin for QGIS using python, But need some help. I want to 1.Zoom to a layer and 2. make the layer visible(switch off or on). But I have managed to zoom to an active layer by setting the active layer in the map canvas by using.               import qgis     canvas = qgis.utils.iface.mapCanvas()     acl = canvas.layer(1)     qgis.utils.iface.setActiveLayer(acl)     > True     qgis.utils.iface.zoomToActiveLayer()      But this only work as long as the layer is viable(obviously since I use canvas). So I need another way to "point" to a layer and then zoom to and switch off. Thanks.


title: Data driven page error when adding world_file

text: my python script for data driven pages causes a runtime error when I set `world_file=True`, or `geoTIFF_tags=True`. Else it runs and shots the ddp - in this case a grid of my city. I'm searching the web for a while but don't find a propper hint, so I ask for your help. Do I have to edit the dataframe? The Problem looks like this:               import arcpy     export_path = r"L:\Bilder\Tile"        field_name = "PageName"      mxd = arcpy.mapping.MapDocument("current")      for i in range(1, mxd.dataDrivenPages.pageCount + 1):        mxd.dataDrivenPages.currentPageID = i        row = mxd.dataDrivenPages.pageRow        print row.getValue(field_name)        arcpy.mapping.ExportToTIFF(mxd, export_path + "_" + row.getValue(field_name) + ".tif", world_file=True)      del mxd      > A1 Runtime error Traceback (most recent call last): File "", line 9, in File > "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\utils.py", line 181, > in fn_ return fn(*args, **kw) File "c:\program files > (x86)\arcgis\desktop10.2\arcpy\arcpy\mapping.py", line 1374, in ExportToTIFF > layout.exportToTIFF(*args) File "c:\program files > (x86)\arcgis\desktop10.2\arcpy\arcpy\utils.py", line 181, in fn_ return > fn(*args, **kw) File "c:\program files > (x86)\arcgis\desktop10.2\arcpy\arcpy\arcobjects\mixins.py", line 596, in > exportToTIFF return self._arc_object.exportToTIFF(*args) TypeError: > PageLayoutObject: Fehler beim Ausf\xfchren von ExportToTIFF


title: How to convert Python (ArcPy) script into tool or toolbox script?

text: I've written a simple python script that allows me to quickly create a copy of a template Mxd, pointing to different shapefiles. The user specifies the new location, a name for the new mxd, then the folder that the shapefiles are stored in and the name of a point, polyline, and polygon shapefile. This script works fine.               newFilePath = r"X:\folder\outputfolder"     mxdName = "New_MXD"     shpFolder = r"X:\folderpath\templates"     shape_points = "input_P"     shape_polylines = "input_L"     shape_polygons = "input_G"                         template_points = r"X:\folderpath\template_P.shp"     template_polylines = r"X:\folderpath\template_L.shp"     template_polygons = r"X:\folderpath\template_G.shp"          import arcpy      MyMXD = r"X:\folder\template.mxd"     mxd = arcpy.mapping.MapDocument(MyMXD)       outMxd = newFilePath +"\\"+ mxdName+".mxd"     for lyr in arcpy.mapping.ListLayers(mxd):         if lyr.supports("DATASOURCE"):         if lyr.dataSource == template_points:             print lyr.dataSource             lyr.replaceDataSource(shpFolder, 'SHAPEFILE_WORKSPACE', shape_points)             print "data replaced"         if lyr.dataSource == template_polylines:             lyr.replaceDataSource(shpFolder, 'SHAPEFILE_WORKSPACE', shape_polylines)                          if lyr.dataSource == template_polygons:             print lyr.dataSource             lyr.replaceDataSource(shpFolder, 'SHAPEFILE_WORKSPACE', shape_polygons)          mxd.saveACopy(outMxd)      I would like to turn this script into a tool that my colleagues can run without having to open the script. I've seen Esri guides for adding scripts to a toolbox, and for creating python toolboxes, but its not clear which I need to do. Any suggestings would be great. I've tried the first option, and replaces the first block with:               NewFilePath = arcpy.GetParameterAsText(0)     mxdName = arcpy.GetParameterAsText(1)     shpFolder = arcpy.GetParameterAsText(2)     shape_points = arcpy.GetParameterAsText(3)     shape_polylines = arcpy.GetParameterAsText(4)     shape_polygons = arcpy.GetParameterAsText(5)      When I tick "Run Python script in process" the script runs but nothing happens. When I untick this and tick "show command window when excecuting script" the script fails and I get ERROR 000714: Error in script ViewSismageExport. Is setting up parameters more complicated than this? I first tried setting the folders to be folders and the shapefiles to be shapefiles, but I've not got everything as string, with it still not working. Does anyone know what I'm doing wrong?


title: Why do the raster band values have different values than the raster statistics?

text: I am working on a raster file (1 Band). According to the raster metadata, my minimum value is 0 and maximum is 7458. This seems to be confirmed by the Style options that return some category values in the value range mentioned above. In addition, these numbers make total sense with the nature of my data. However, when I use the "Identify Features" tool, I get some values completely out of the range (these values look like this one generally: 4.5761380462188e-315). I don't understand why and it's a problem because I would like sometimes to get the value of a single pixel. Do you have an idea why I have this discrepancy between the Identify Feature tool and the raster statistics? I use QGIS 2.2.0-Valmiera on Windows 7. Thank you


title: Configuration with geoserver openlayers

text: I did the following: 1 - I installed geoserver by the installer for Windows as a service and standard folders . 2 - Open `http://localhost:8080/geoserver/www/ol-demo.html` to test my geoserver WMS layers and not appear the Google , Yahoo , Earth and WMS layers from my geoserver . 3 - I copied into the www folder with examples of open layers begginers , this folder contains the img folder , theme , metadata , file openlayers.js OpenLayers and my code examples from the book . So I opened and ran `http://localhost:8080/geoserver/www/code/chapter_5_ex3_google_layer_mashups.html` Google Layer , the overlay of WMS ( example , not my local layer ) . questions : 1 - I used the Windows installer , it installed the GeoServer is all right ? Thus no Tomcat or other container . 2 - What can be wrong in the settings of my GeoServer to not be serving their layers ? Installed according to the instructions from the official website , to run manually as a service , with personalized folder in the default folder and no attempt worked. 3 - In this example `http://localhost:8080/geoserver/www/ol-demo.html` he not only calls the WMS layer (which is on my localhost ) , but also calls Google , Yahoo and Earth these should appear on the map . Why do not I ? My initial goal is to make the overlay of my WMS layers provided by PostGIS on Google Maps, where I did not succeed . So I decided to use the demo to see if my geoserver is working. Any help is welcome , as I am with this problem for days and I do not know that attempts to do.


title: How can I measure area in square kilometres from lat / lon gps coordinates?

text: I have code from openlayers that I cannot quite handle because I am so new for these kind of codes. My main goal is to measure area from n amount of lat/lon gps coordinates and convert the result into square kilometres. I have adopted some code from the post here and it seems almost ready but I was not able to figure out how to use it with my test coordinates that are               var lats = [25.767368, 34.088808, 40.727093]     var lons = [-80.18930, -118.40612, -73.97864];      Which I need to use in code like this to get the area and finally convert it into square kilometres. Any help with this? Please make it simple as this is pretty much my first attempt on geagraphic coding.               var area = 0.0;     var len = ring.components && ring.components.length;          if (len > 2) {          var p1, p2;          for (var i=0; i<len-1; i++) {              p1 = ring.components[i];         p2 = ring.components[i+1];         area += OpenLayers.Util.rad(p2.x - p1.x) *             (2 + Math.sin(OpenLayers.Util.rad(p1.y)) +             Math.sin(OpenLayers.Util.rad(p2.y)));          }          area = area * 6378137.0 * 6378137.0 / 2.0;          }


title: Not able to bulk download MODIS images

text: I am trying to download Modis NDVI and EVI images. I have used the following R-code:                   setwd("/home/karthik/Valmiki/Modis_bands/")         library(raster)         library(RCurl)         source("ModisDownload.R") ##ModisDownload.R is in the setwd directory#         X="MOD13Q1"         getMODIS(x=X, h=25, v=6, dates = c('2000.02.18', '2014.07.31'),  version='005')      I am getting the following error:                   Error in readChar(con, 5L, useBytes = TRUE) : cannot open the connection         In addition: Warning message:         In readChar(con, 5L, useBytes = TRUE) :         cannot open compressed file 'ModisLP.RData', probable reason 'No such file or directory'      I am using Ubuntu 12.04. Please let me know me where I am going wrong.


title: How to add tree height values to raster / vector data?

text: I am new to working with raster data and I am trying to teach myself a few things. I have what I think is a simple problem but I can't seem to figure this one out after a few hours of searching Arc Help etc. I have raster with ground elevations to which I want to add tree heights. I have a Forest Cover polygon vector layer which has tree heights for my study area. I would like to add the heights of trees to the elevation data. My first guess was to convert the Forest Cover polygon to a raster based on the tree height field and then Merge/Combine the two rasters? I have esri's 3d Analyst so I think I would have all the available tools? I am using ArcMap 10.0


title: Why does the QgsVectorFileWriter have probelms creating vector files?

text: I am using QGis 2.0.1 and have tried to create a new GIS feature file (SQLite, Spatialite, and Shapefile) of a subset of the data. I am then saving a copy of the style of this layer and applying it to the new layer. Also, I am starting with a Spatialite database. All three have different issues: SQLite - The data is copied, but the fields loose their edit widget type (If you go in to the properties of the layer, then the Fields menu on the left; it is the Edit widget drop down.) Shapefile - The file has no features in the table, but everything else is fine. Spatialite - QGis locks up. Shapefile code:                   writer = QgsVectorFileWriter(layer_path,                                      'CP1250',                                      fields,                                      QGis.WKBPoint,                                      crs,                                      'ESRI Shapefile')              if writer.hasError() != QgsVectorFileWriter.NoError:             raise IOError('Can\'t create the file: {0}\n'.format(layer_path))              for feature in features:             writer.addFeature(feature)              del writer      SQLite code:                   writer = QgsVectorFileWriter(layer_path,                                      'CP1250',                                      fields,                                      QGis.WKBPoint,                                      crs,                                      'SQLite')              if writer.hasError() != QgsVectorFileWriter.NoError:             raise IOError('Can\'t create the file: {0}\n'.format(layer_path))              for feature in features:             writer.addFeature(feature)              del writer      Spatialite code:                   spatialite = ["SPATIALITE=YES"]              writer = QgsVectorFileWriter(layer_path,                                      'CP1250',                                      fields,                                      QGis.WKBPoint,                                      crs,                                      'SQLite', spatialite)              if writer.hasError() != QgsVectorFileWriter.NoError:             raise IOError('Can\'t create the file: {0}\n'.format(layer_path))              for feature in features:             writer.addFeature(feature)              del writer      I have also tried `spatialite = ["SPATIALITE=YES", ]` I just want to know if there is anything else that I can do to get this to work the way I want.


title: How to get label text for each feature in vb.net in Arcview

text: I am trying to get text for the label that Arcview is going to place for each feature when the labeling is on, so that I can write it to a file. Since the labels may have complex operations involved, I don't want to go to the attributes. Any help would be appreciated. I am working in vb.net Thanks


title: Consume netcdf Data with ESRI Javascript API?

text: We've got loads of raster data that can be organized by time, so we thought packaging them into netCDF would be a good a idea. But we also want to serve out the data to a web application that leverages the ArcGIS javascript api. I've seen the ncWMS page, but the info there doesn't exactly make it clear how to consume the ncWMS. Anybody have experience with this or/could point me to an example?


title: Extract all the intersections from a polyline shapefile

text: I would like to extract all the intersections from a polyline shapefile and print their addresses as text using either ArcDesktop or Quantum GIS, is there a straightforward way to do this? My polyline has address information. Its a "roads" file, I don't need to compare it to other layers.


title: C# ArcObjects RubberBand Tracking Cancel

text: I've written a function to create a polygon based on geometry collected from the IRubberBand interface as follows:               ITrackCancel esc = new CancelTrackerClass();     esc.CancelOnKeyPress = true;     esc.Reset();     ArcMap.Document.ActiveView.ScreenDisplay.CancelTracker = esc;     //MessageBox.Show("This will allow a key press to cancel the tracking if uncommented");     IRubberBand rubberBand = new RubberPolygonClass();     IScreenDisplay screenDisplay = ArcMap.Document.ActiveView.ScreenDisplay;     IGeometry geometry = rubberBand.TrackNew(screenDisplay, null);      I am attempting to cancel the tracking of the polygon (i.e. cancel the drawing operation) on key press. The above code will not cancel the tracking unless I insert a MessageBox display somewhere before the tracking begins. My guess is that this problem is somehow related to this issue: How to use ITrackCancel to escape IRubberBand track? . I'm also a bit skeptical of how I've hooked into the active view. If anyone could provide some insight on how to best cancel this operation without an awkward messagebox being involved I would certainly appreciate it.


title: table of polygon centroids

text: Given a polygon layer, census tracts, how can I generate a table of centroid coordinates in QGIS, or add columns with the centroid coordinates or generate a new layer of polygon centroid points?


title: Trouble with making shapefile from xls file

text: Using DIVA GIS, each time I attempt to import some species data from an xls file to be plotted onto a map, the program becomes very slow and eventually becomes unresponsive. My dataset originally included 15,000 separate points so I tried importing only 2,000 of them and the same error occurred. Am I simply trying to import too much data at once? If so, can this be overcome? I then intended to analyse my species distribution against the imported climatic data to discover a statistical relationship.


title: PHP MapScript API Development for MapServer

text: PHP MapScript Tutorial I'd like to use the PHP MapScript API and am curious how others have successfully used it to produce a site that serves up remote sensing data. I've got land surface temperatures, emissivities, and reflectances I'd like to share via a mapserver. I've already installed mapserver on an apache webserver and have also installed the MapScript shared library. What is a good development environment to use for PHP MapScript development? Any example sites that use PHP MapScript people know of? How useful is the MapScript API and how easy is it to use? Thanks in advance!


title: Streamlining Python Code for Big Data

text: I have Python code which is designed to take point shapefiles through the following workflow:   1. Merge points   2. Integrate points, such that any points within 1 m of each other become one point   3. Create Feature layer, where points with z < 10 are selected   4. Buffer Points   5. Polygon to raster 1m resolution   6. Reclassify, where 1 - 9 = 1; NoData = 0 Each shapefile has approximately 250,000 to 350,000 points covering ~ 5x7 km. Point data used as inputs represent tree locations. Each point (i.e. tree) has an associated "z" value which represents crown radius and is used in the buffer process. My intent is to use the final binary output in a separate process to produce a raster describing canopy cover. I ran a test with four shapefiles and it produced a 700MB raster and took 35 minutes (i5 processor & 8GB RAM). Seeing that I will need to run this process on 3500 shapefiles, I would appreciate any advice to streamline the process (see attached code). **Generally speaking, what is the best way to deal with geoprocessing big data? More specifically, are there any tweaks to the code or workflow that could help increase efficiency?** **Edit** : Time (% of total) for geoprocessing tasks:   * Merge = 7.6%   * Integrate = 7.1%   * Feature to Lyr = 0   * Buffer = 8.8%   * Poly to Raster = **_74.8%_**   * Reclassify = 1.6% ## ![enter image description here](http://i.stack.imgur.com/rWD9M.png)               # Import arcpy module     import arcpy          # Check out any necessary licenses     arcpy.CheckOutExtension("spatial")          # Script arguments     temp4 = arcpy.GetParameterAsText(0)     if temp4 == '#' or not temp4:         temp4 = "C:\\gdrive\\temp\\temp4" # provide a default value if unspecified          Reclassification = arcpy.GetParameterAsText(1)     if Reclassification == '#' or not Reclassification:         Reclassification = "1 9 1;NODATA 0" # provide a default value if unspecified          Multiple_Value = arcpy.GetParameterAsText(2)     if Multiple_Value == '#' or not Multiple_Value:         Multiple_Value = "C:\\t1.shp;C:\\t2.shp;C:\\t3.shp;C:\\t4.shp" # provide a default value if unspecified          # Local variables:     temp_shp = Multiple_Value     Output_Features = temp_shp     temp2_Layer = Output_Features     temp_Buffer = temp2_Layer     temp3 = temp_Buffer          # Process: Merge     arcpy.Merge_management(Multiple_Value, temp_shp, "x \"x\" true true false 19 Double 0 0 ,First,#,C:\\#########omitted to save space          # Process: Integrate     arcpy.Integrate_management("C:\\gdrive\\temp\\temp.shp #", "1 Meters")          # Process: Make Feature Layer     arcpy.MakeFeatureLayer_management(temp_shp, temp2_Layer, "z <10", "", "x x VISIBLE NONE;y y VISIBLE NONE;z z VISIBLE NONE;Buffer Buffer VISIBLE NONE")          # Process: Buffer     arcpy.Buffer_analysis(temp2_Layer, temp_Buffer, "z", "FULL", "ROUND", "NONE", "")          # Process: Polygon to Raster     arcpy.PolygonToRaster_conversion(temp_Buffer, "BUFF_DIST", temp3, "CELL_CENTER", "NONE", "1")          # Process: Reclassify     arcpy.gp.Reclassify_sa(temp3, "Value", Reclassification, temp4, "DATA")


title: How to check which version of GDAL do I have and how many?

text: I have installed in my computer and running the following software: PosgresSQL 9.3.2.3 PostGIS 2.1.1.1 QGIS 2.0.1 Dufour (this one installed through the OSGeo4w package) QGIS 1.7.0 Wroclaw I am concerned as far as I know that all of them somehow operate, or use, or access or read the GDAL library. does it mean I must have different versions of GDAL running depending on which one I am using? If I open the OSGeo4W Shell it says GDAL 1.10.1, how can I check the other versions if any? Perhaps, all of them use only the latest one. Ideally I would like to know how to check this and by which means.


title: How to find out the co-ordinate system using .prj file

text: I'm very new to GIS concepts. Is there a way or procedure of finding out whether the co-ordinate system is in UTM, latitude/longitude or in stateplane using .prj file? My .prj file contains the following text: > GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS > 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]] Thanks in advance for the help.


title: Implementaion of GSCoverageEncoder for WorldImage in Java

text: I want to upload a raster image to GeoServer as WorldImage from my Java application. I'm using geoserver-manager v1.4.2 - a Java-based implementation of GeoServer's REST API. It provides convenient methods to work with the GeoServer REST API. The problem is that the support for the raster image types is not implemented (only Image Mosaic is supported , which doesn't suite me). Does anyone has an example implementation of GSCoverageEncoder for an image with world file?


title: How to select only Spatial tables from the postgres database?

text: My database contains both Spatial and Non-Spatial tables, but I want to retrieve only Spatial table from the query. Any suggestions to select only Spatial tables. 'the_geom' is the geometry column in the spatial table. Otherwise, is it possible to select tables from its column name. I tried with this code `select relname from pg_stat_user_tables WHERE schemaname='public'` ; but from this we get all table names. Thanks in advance..


title: Is it possible to have a multi-modal network in QGIS?

text: I would like to design a multimodal network in QGIS and connect it with real- time traffic data (for road, rail and inland waterway). Is it possible to have a multimodal network in QGIS? And how can I get the real-time traffic data in Austria? I already have a multimodal network in ArcGIS and I know that there is DTF files from NAVTEQ for roads. Is there any similar file formats for rail and inland waterway?


title: Find directly joined neighbour points using a different line layer in QGIS

text: I have two layers of data:   1. A set of points showing every significant point on the UK rail network (Network Rail Stanoxes). These are a mixture of stations, junctions, timing points and signals.   2. An OS OpenData line layer showing the UK rail network. ![Layers](http://i.stack.imgur.com/0GEla.png) My aim is to use QGIS to find how the points in layer 1 are joined. I need to take each point in layer 1 and travel down each line in layer 2 until I reach another point. The link between points then needs to be saved. A point can be linked to one or more other points, creating the whole UK network. A fairly simple example is shown below for the south east of the UK. The output I'd like to generate is the direct links between points: A-B, B-C, C-D, C-E ![QGIS Screenshot](http://i.stack.imgur.com/dbFR0.png) I'm a complete QGIS novice, so apologies if this is either completely obvious or impossible!


title: can we add the feature layer from code

text: I need to add a feature layer to arcgisonline account from the code. i have created a developer account where if i want to create a feature layer i have to upload a CSV file and it makes the feature layer. > Problem is: i need to create the sane feature layer using my code. like is > their any REST API or some API which ArcGIS provide this option?


title: Are there any tools for editing Sql Server spatial types in Silverlight?

text: I see from this post that I can tweak RIA Services to transmit Sql Server spatial types. What I'd like to do is edit these spatial types. Does anyone know of any tools that would run in Silverlight that would support this?


title: GIS file format for space and time data?

text: I need some advice with GIS file format choice. Here is my problem: I have about 100 stations. At each station water temperature is measured daily for 1 year I want to load this data to QGIS and be able to quickly display a water temperature map for a user-selected day. I came up with following file structure: A shapefile with lat, lon, time, value columns. For each combination of lat, lon, time there is a new point shape. So for each station I have 365 "overlapping shapes". Then I can use the "query" in QGIS to filter out shapes that have the time attribute equal to the selected day. This is OK for 365 days but becomes little inefficient for several years (shapefile becomes huge and QGIS displays it very slowly) Is there any other GIS format (like Spatialite or Postgis) that I could use to store the time-space data more efficiently and still be able to visualize it in QGIS without extra programming?


title: Using arcpy.mapping.ExportToTIFF() to create dataframe geotiff

text: Has anyone created a dataframe geotiff using arcpy.mapping.ExportToTIFF() with satisfactory results? My output has poor quality and the georeferencing is significantly off. This example is the geotiff at 50% transparency on top of my data. The shift between "South El Monte" is the best shift example. ![enter image description here](http://i.stack.imgur.com/Q2qiZ.jpg) Here is an example of the code:               currentMxd = arcpy.mapping.MapDocument(r"F:\Test\amp_Edit.mxd")     dataFrame = arcpy.mapping.ListDataFrames(currentMxd, "Layers")[0]      arcpy.mapping.ExportToTIFF(currentMxd, r"F:\Test\outtif.tif", dataFrame, df_export_width = 46000, df_export_height = 29000, resolution = 1000, geoTIFF_tags = True)


title: WFS Mapguide Maestro. Can't query features?

text: I'm using MapGuide Maestro 5.0 and Geoserver 2.4.0. I have added a new WFS feature source to the viewer, and used the following url:                http://localhost:8080/geoserver/scotwat/wfs?      In the vector layer def it shows the fields and these have been checked to be visible. The viewer shows the features, but when I click on them to select, nothing happens. I have looked through the Properties of both the feature and layer, and cannot find an option to turn it on. I have also ensured that the box 'Features on the layer are selectable (if visible)' is checked. I also added the wfs to qgis, and can select features through there. EDIT: I have progressed slightly, but am still stuck. The problem lies with the primary key of the WFS service. Within geoserver the option 'Expose primary key' is checked. MapGuide produces the following warning the feature source is validated:                No primary key defined for class: scotwat:paris_hetero, features will not be selectable and feature joins on this class may produce unexpected results      The data source the WFS is generated has a primary key, but this has not been seen by MapGuide Maestro. Any help? If I know what the primary key how can I define it in MapGuide Maestro?


title: WMS behind proxy

text: I have problem with some kind of WMS. It serv images like WMS but only images - no xml file on GetCapabilities - and require additional params like KEY and PERM because is behind PHP proxy.               https://geo.powiatwalecki.pl/is/map.php?     SERVICE=WMS     &VERSION=1.1.1     &REQUEST=GetMap     &LAYERS=dz,dz_nr     &SLD_BODY=dz/rgb%28255,255,255%29/rgb%28255,255,255%29/x/0/1px/x     &KEY=915a087d219ca2c19160a8d2bd7817c5     &PERM=f;f;f;f;f;u1,u2,u3,u4,u5,u6,u7,u8,u17,u9,u10,u11,u12,u13,u14,u15,u16     &SRS=EPSG:2176     &BBOX=5597964.471,5905671.349,5598075.999,5905704.819     &WIDTH=1916     &HEIGHT=575      Here is link from code above Image KEY parameter change every day so here is link to OpenLayers client using this service Question is how to add this to QGis to use it like wms layer. I'm trying through gdal xml file, but without effect.


title: ArcGIS 10.1 edit PostgreSQL layer

text: I've been working with QuantumGIS/PostgreSQL/postGIS for a while now, without problems. But now I have to change QuantumGIS for ArcGIS 10.1. I'm having trouble adding a layer from PostgreSQL and being able to edit that layer. Basically I want to work directly with the DB, if I change something in the map it should be reflected on the database. Can I do this without a Geodatabase? With a Geodatabase? I'm having trouble finding precise information. Any help would be appreciated, Thanks


title: get value from attribute table and change symbology with QGIS and Python

text: I'm newbie developing plugins for QGIS using Python, this is the first one and it is very frustrating hehe. Ok, this is the situation, I have a 1 point layer in QGIS. Its attribute table has 4 fields and 1 row:               ID|Estation|Data|Estate          1     1     20  "Green"      What I'm trying to do with my plugin is to get the string data from "Estate" field and check the string data: if string is "Green", the one point layer symbology will turn to green. if string is not "Green", the one point layer symbology will turn to red. The last part of this plugin is to refresh every X seconds to check "Estate" field, and change the symbology of this point layer if "Estate" field has changed. This is what I have done so far in my python script, I put from the run method:               # run method that performs all the real work     def run(self):              # create and show the dialog         dlg = PlayMonitoringDialog()          # show the dialog         dlg.show()         result = dlg.exec_()          # See if OK was pressed         if result == 1:            # do something useful (delete the line containing pass and           # substitute with your code           def h():                  self.iface.mapCanvas().refresh             layer = self.iface.activeLayer()             Estado = layer.attributeDisplayName(3)             if Estado == "green":               props = { 'color' : '0,255,0' }               sl = QgsSymbolLayerV2Registry.instance().symbolLayerMetadata("SimpleMarker").createSymbolLayer(props)               s = QgsMarkerSymbolV2([sl])               layer.setRendererV2(QgsSingleSymbolRendererV2(s))             else:               props = { 'color' : '255,0,0' }               sl = QgsSymbolLayerV2Registry.instance().symbolLayerMetadata("SimpleMarker").createSymbolLayer(props)               s = QgsMarkerSymbolV2([sl])               layer.setRendererV2(QgsSingleSymbolRendererV2(s))             if hasattr(layer, "setCacheImage"): layer.setCacheImage(None)             layer.triggerRepaint()           lanza()                 def lanza():         t=threading.Timer(3.0,h)         t.start()       lanza()      It doesn't work :(, so I hope someone can help me. I have seen the small examples in http://www.qgis.org/pyqgis-cookbook, I have checked http://qgis.org/api/index.html and I have looked for examples in the Internet, but I still don't know how to solve the problems in my script. Thank you very much in advance and sorry for my english, I haven't practice for a while. Lucas


title: "Hardclip" a TIN

text: How do you go about hardclipping a TIN in Arc10? Previously there was a modify TIN tool, but I cannot seem to find it. So my question remains, how do I go about clipping a TIN with a polygon layer?


title: Algorithm to determine if there are any landlocked polygons

text: Imagine that I have the need to determine landlocked countries. Is there a to do it automatically? Any known implementations? The target platform on this will be ESRI.


title: How do I create a geodata service from a certain Version?

text: This is with ArcGIS for Server and Desktop 10.1 I want to host a geodata service, but I don't want it to be from the DEFAULT version. I want it to be from an alternate version, that I can check the changes on before I post it to the DEFAULT. I don't see any options to do this. Is it possible?


title: Which Geographic Transformation to use for WGS_84 to D_Sphere?

text: I have a shapefile for world countries and I want to reproject it from GCS_WGS_1984 to Sphere_Eckert IV. However, it asks for Geographic Transformation; it says this is optional but it is required in order to proceed. I saw that the datum is D_Sphere but I cannot find a combination between this and WGS in the provided list. Can please someone tell me what the right choice is? I use ArcGIS 10.1. Thanks!