title: Create a named map in CartoDB

text: I'm trying to create a template map following this tutorial from the cartoDB's maps API. My template.json file               {       version: "0.0.1",       name: "template_1",       // embedded authorization certificate        auth: {        method: "open"        },       placeholders: {           color: {             type:"css_color",             default:"blue"           },           del_s: {             type:"sql_literal",             default: "19420210"           },           del_n: {             type:"number",             default: 19420210           }       },       layergroup: {         "version": "1.0.1",         "layers": [{           "type": "cartodb",           "options": {             "cartocss_version": "2.1.1",              "cartocss": "#layer { polygon-fill: <%= color %>; }",             "sql": "select * from table1 WHERE del = <%= del_s %>",             "interactivity": ["cartodb_id","del","gk"]           }         },         {           "type": "cartodb",           "options": {             "cartocss_version": "2.1.1",              "cartocss": "#layer { polygon-fill: red; }",             "sql": "SELECT * from table2 WHERE del = <%= del_n %>",             "interactivity": [ "cartodb_id","del"]            }         }]       }      }      When I try to make the request to create the template..               curl -X POST \        -H 'Content-Type: application/json' \        -d @template.json \        'https://docs.cartodb.com/api/v1/map/named?api_key=568kaka8bc89'      I can not connect, I get an error:               curl: (7) Failed connect to docs.cartodb.com:443; Connection timed out      anyone knows what I'm doing wrong? UPDATE Try with my username, and get an "internal server error"               curl -X POST \        -H 'Content-Type: application/json' \        -d @template.json \        'https://myusername.cartodb.com/api/v1/map/named?api_key=myapikey'


title: Is it possible to share temporary (InMemoryWorkspace) FeatureClass between Threads?

text: I'm trying to run some operations (functions) in seperate threads so that the UI is not blocked. One of this operation is the generation of a feature class copy from an existing feature class into a temporary in memory feature class (used for processing later). > In cases where you must pass ArcObjects components from the main thread into > a worker thread, serialize the object into a string, pass the string to the > target thread, and deserialize the object back. I'm using XMLPersistedObject to pass an serialized IName object to the copy function. The FeatureClass is then opened via IName.Open() and copied (structure + rows) into a temporary Workspace (created via the InMemoryWorkspaceFactory). The copied FeatureClass is then serialized back into an IName object again an returned to the main thread. There IName.Open() fails with a generic COM Error (-2147467259). The whole thing works when I use a different Workspace (e.g. ScratchWorkspaceFactory), but the copy process (see code below - maybe this can be improved?) and subsequent operations take longer than. So, is there a possibility to share temporary FeatureClasses between Threads?               public class DataManagement     {         #region Public Methods and Operators              public static string CreateAndCopyToTemporaryFeatureWorkspace(string name)         {             IFeatureClass featureClass = ArcObjectHelper.CreateNameFromXML(name).Open();                  return ((IDataset)CreateAndCopyToTemporaryFeatureWorkspace(featureClass)).FullName.ToXMLString();         }              #endregion              #region Methods              private static IFeatureClass CreateAndCopyToTemporaryFeatureWorkspace(IFeatureClass featureClass)         {             IFeatureClass localFeatureClass = featureClass;                  // Copy FeatureClass structure             IFeatureClass featureClassNew = CreateTemporaryFeatureClassCopy(localFeatureClass);                  // Copy rows             using (ComReleaser comReleaser = new ComReleaser())             {                 IFeatureCursor featureCursor = localFeatureClass.Search(null, false);                 comReleaser.ManageLifetime(featureCursor);                      IFeature feature;                      IFeatureBuffer featureBuffer = featureClassNew.CreateFeatureBuffer();                 comReleaser.ManageLifetime(featureBuffer);                      IFeatureCursor insertCursor = featureClassNew.Insert(true);                 comReleaser.ManageLifetime(insertCursor);                      while ((feature = featureCursor.NextFeature()) != null)                 {                     featureBuffer.Shape = feature.Shape;                          for (int i = 0; i < localFeatureClass.Fields.FieldCount; i++)                     {                         if (localFeatureClass.Fields.Field[i].Type == esriFieldType.esriFieldTypeOID)                         {                             continue;                         }                              featureBuffer.Value[i] = feature.Value[i];                     }                          insertCursor.InsertFeature(featureBuffer);                 }                      featureCursor.Flush();             }                  return featureClassNew;         }              private static IFeatureClass CreateTemporaryFeatureClassCopy(IFeatureClass featureClass)         {             IWorkspaceFactory workspaceFactory = new InMemoryWorkspaceFactory();             IFeatureWorkspace featureWorkspace = ((IName)workspaceFactory.Create("", "InMemoryWorkspace", null, 0)).Open();                  IFields fields = new Fields();             IFieldsEdit fieldsEdit = (IFieldsEdit)fields;             fieldsEdit.FieldCount_2 = featureClass.Fields.FieldCount;                  for (int i = 0; i < featureClass.Fields.FieldCount; i++)             {                 if (featureClass.Fields.Field[i].Type == esriFieldType.esriFieldTypeOID)                 {                     IField field = new Field();                     IFieldEdit fieldEdit = (IFieldEdit)field;                     fieldEdit.Name_2 = "ObjectID";                     fieldEdit.Type_2 = esriFieldType.esriFieldTypeOID;                     fieldsEdit.Field_2[i] = field;                 }                 else                 {                     IClone clone = (IClone)featureClass.Fields.Field[i];                     IField clonedField = (IField)clone.Clone();                     fieldsEdit.Field_2[i] = clonedField;                 }             }                  IUID uid = featureClass.CLSID;                  Random random = new Random((int)DateTime.Now.Ticks);                  return featureWorkspace.CreateFeatureClass(random.Next(1000000, 10000000).ToString(CultureInfo.InvariantCulture),                 fields, (UID)uid, null, esriFeatureType.esriFTSimple, "shape", "");         }              #endregion     }


title: Get the distance between 2 WGS84 GEOSGeometries

text: I currently want to check if 2 given `GEOSGeometry` objects of type `Point` are within 100m of each other. (length flexible) Looking at previous asked questions on this topic, most of the answers were totally not specific and even worse sent the author to other frameworks. My goal is to stay with GEOS / GDAL (geo django) to solve this problem. The current code I use:               from django.contrib.gis.geos import Point, GEOSGeometry     pnt = GEOSGeometry('SRID=4326;POINT(40.396764 -3.68042)')     pnt2 = GEOSGeometry('SRID=4326;POINT( 48.835797 2.329102  )')     pnt.distance(pnt2) * 100      Found on gis.SO. It is very difficult for me to confirm that this code works, because I would have to check the distances between to points and check the results against each other. **The main goal is to take 2 points and get the distance in meters based on the WGS84 projection and the GEOS or GDAL framework**


title: ArcGIS Export Selected Raster Graphics

text: I'm having a problem exporting raster graphics using ArcGIS 10. It seems simple enough, but for some reason I cannot get ArcMap to work with me. Here are my steps:   1. Add a raster to ArcMap   2. Open the raster's Attribute table and select some values   3. Right click on the raster in the ArcMap table of contents --> Data --> Export Data   4. I am given 3 options for the data extent (a) Data Frame, (b) Raster Dataset, (c) Selected Graphics. The last option (which is the one that I want) is grayed out.  I have not been able to activate the "Selected Graphics" option in the Export Raster Data dialog. The strange thing is that my colleague has the exact same version of ArcMap and when repeating the steps above, the option to export Selected Graphics is not grayed out. I'm assuming that there is some setting that I have to adjust. Does anyone know what can be causing this?


title: ArcMap, ArcCatalog very slow to open on new laptop with ample resources?

text: I have a recent installation or ArcGIS Desktop on a 64 bit Windows 7 Machine, with 8GB ram, intel core i7 processor, and a solid state drive. It takes ~2 min to open either ArcMap or ArcCatalog from the start menu (e.g. without loading any data). As a new installation, I have only the default folder connections. Is this normal, or is it possible that there is something wrong with my installation?


title: Need help adding OpenLayers to Geonode

text: I just installed GeoNode and really love the manner in which it helps me arrange the contents related to the map but I am not sure about adding OpenLayers to it. Well here is where I am faltering,   1. I am not sure where the OpenLayers code goes? I felt that I could use the code to customize GeoNode's look and feel and make use of the customization features on the Web Page, but to publish the Map you receive a embed link which only gives you the Map and not the other features.   2. If we have to use OpenLayers separately, like in the HTML, I am not sure, How to run that static html file, since the localhost is being used by the GeoNode.


title: PostGIS Uploader - Need Many .shp files in One table

text: How do I put many shape files into 1 table? The PostGIS plug-in shape file uploader makes one seperate table each time it uploads a shape file.


title: Order of KML Layers in Google Maps API V3

text: I'm working on a website that when finished will allow the user to turn on and off a number of different KML layers using checkboxes, the layers are parsed through geoXML3 and I want to know how to control the order which the layers appear. Currently if I have 2 polygon layers (layer A and B) and I click the checkbox to load layer A then I click the checkbox to load layer B layer B will be on top as it should be. However if I then uncheck both boxes, and select B first and then A, B is still drawn on top. I assume this is because I parse the data only the first time the boxes are checked and afterwards I just show or hide them, but is there a way to manage which layer is on top? Below is an example of the function I use to parse the layer and then turn it on or off.               function boundaryToggle_FOS() {     if (boundaryIsLoaded_FOS == 0)     {     boundaryIsLoaded_FOS = 1;         if (firstLoadFOS == 0)             {             kml_FOS.parse('FOS_Boundary_v03.kml');             firstLoadFOS = 1;             }         else             {             kml_FOS.showDocument();             }     }     else     {         kml_FOS.hideDocument();         boundaryIsLoaded_FOS = 0;     }}


title: OGR: Meaning of "2.5D extension as per 99-402"

text: I'm looking at the enum defs in `ogr_core.h`. I see a bunch with comments that read: > 2.5D extension as per 99-402. What in the world does that mean? Presumably it's referencing a spec or RFC?               typedef enum      {       ...       wkbPoint25D = 0x80000001, /**< 2.5D extension as per 99-402 */       wkbLineString25D = 0x80000002, /**< 2.5D extension as per 99-402 */       wkbPolygon25D = 0x80000003, /**< 2.5D extension as per 99-402 */       ...     }


title: How to fix the "Snap tolerance" error message in QGis?

text: I am trying to use the "Delete part" tool on a shapefile. While I am starting to use it, an error message shows up : "Could not snap segment. Have you set the tolerance in Settings > Project Properties > General?". I tried to set different tolerances in the snapping options, the error still shows up and I can't delete any part of my polygon. How to fix that? I use QGis 1.8.0-Lisboa and Windows 7. Thanks


title: trying to make a point to polygon attribute join with a many to many match

text: I am trying to match US charter schools to noncharter schools in QGIS. I have constructed a buffer of 10 miles around each charter schools and then tried a join attributes by location match with non charter schools. This has been partly successful. I get matches but only one charter school buffer zone is associated with each traditional school. What I really want to have is a match between each traditional school and all the charter school buffer zones the traditional school is located in. The result will be a table with multiple records for some/many traditional schools because they are located in more than one charter school buffer zone. In short I want to create a many point to many polygon join? thanks in advance, Nat


title: How do I convert Eastings and Northings projection to WSG84 in geopandas?

text: I'm very new to projections and struggling to convert a map dataset from OSGB36 (Eastings and Northings) to WGS84 (lat/longs). I've used geopandas to import a shapefile. Now I want to export it as GeoJSON so I can then use Folium to map it on leaflet.js. If I check the crs variable with .crs I get:               {u'datum': u'OSGB36',     u'k': 0.999601272,     u'lat_0': 49,     u'lon_0': -2,     u'no_defs': True,     u'proj': u'tmerc',     u'units': u'm',     u'x_0': 400000,     u'y_0': -100000}       Trouble is if I'm a bit lost on the attributes to set to convert it to WGS84:               staffs.to_crs({'datum':"wgs84"})          RuntimeError: projection not named      Can you help?


title: Distance raster from one raster value (A) to another value (B) while ignoring similar values

text: I have two different patch types: A and B I want to get a distance raster from all A types to the closest Type B patch but I don't want the raster to calculate from A to A or B to B. Just A to B. Therefore in this simple line below you have 3 A patches and 2 B patches. A patch is a group of like types. A '-' represents a null pixel. In the example below there is a line of 17 pixels.               A A A - - - A A - - - B B B - - A - - - B B      The only way I can think of doing this is to create a distance raster for each patch of type A to patches B and run a min algorithm on them. Therefore in the above example I would create 3 different rasters and find the minimum value at each pixel. There has to be a better way that I'm not thinking about.               Distances from AAA patch     A A A 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19     Distances from AA patch     6 5 4 3 2 1 A A 1 2 3 4 5 6 7 8 9 10 11 12 13 14     Distances from A patch     16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 A 1 2 3 4 5      Then use something like if original raster = B or original raster = A than 0, else min(distance rasters)               Produces output     A A A 1 2 1 A A 1 2 3 B B B 2 1 A 1 2 3 B B      Thanks!


title: How to use arcpy.mapping to batch export layer to KML?

text: I am trying to export map layers from an MXD to KML using arcpy with the following code:               import arcpy, os, sys, traceback          arcpy.env.overwriteOutput = True          mxd = arcpy.mapping.MapDocument(r"C:\Layer_KML\MXD\Test.mxd")          try:         for lyr in arcpy.mapping.ListLayers(mxd):             i = 1 # forgive the wonky naming convention, I just wanted something quick for test purposes             outKML = "C:\\Layer_KML\\KML\\" + "Test" + str(i) + ".kmz"                  # These lines are where the error/s is occuring             newLayer = arcpy.SaveToLayerFile_management(lyr, "C:\\Layer_KML\\Layers\\" + str(lyr) + ".lyr")                     arcpy.LayerToKML_conversion(newLayer, outKML, "20000", "false", "DEFAULT", "1024",     "96")                  i = i + 1          except:             ...      Though my traceback keeps telling me : {ERROR 000623: Invalid value type for parameter in_layer.} for both (I commented out each indivudally to test) the .SaveToLayerFile, and LayerToKML routines. I am assuming these routines want to except a full file path, while the arcpy.ListLayer() only yields a basename. Maybe? I am really not sure. Any thoughts as to why this isn't working? I have already seen this question : How to batch export layers to kml in Arcmap and this is not what I am looking for. I want to do this in Python rather than through a toolbox tool. Any help is appreciated. Thanks!


title: Adding raster image to Raster field with Python?

text: I have a feature class with with field that is of Raster data type. I want to add an image to each record in the feature class. Is it possible to use Python to add the image to the Raster field for each record? I've been searching but can't find anything. Thanks.


title: Use of OpenStreetMap France in another country

text: I've just stumbled across OpenStreetMap France and was amazed to find that it has 2 extra zoom levels, 19 and 20, whereas OpenStreetMap only goes up to 18. In our application we display potentially masses of geographic data over the top of basemaps. We use Ordnance Survey and OpenStreetMap. The Ordnance Survey mapping we use we can render off at higher zoom levels, which will (once we've upgraded our GeoServer) allow our application to zoom closer and see greater details in the overlaid data. But by default we use OpenStreetMap as it is quicker, looks better (especially far zoomed out) and we like trying to get people to use OSM. Sadly this means we can only go to zoom level 18 (unless I do the OpenLayers thing to 'force' it to zoom in which drops the resolution). However, if we were to use OSM France we could get those extra zoom levels, and as far as I can see it is (almost) as quick and doesn't appear to have any language diffrences (apart from London being displayed as Londres, which I am not concerned about as we are based in Dorset). I guess my question is twofold. 1 - Can anyone think of a good reason for me not to use OSM France 2 - Why doesn't standard OSM have these two extra zooms (capacity?) and will they ever?


title: Using NetCDF for point time-series observations?

text: I have a large growing dataset of water and climate observations for several years in regular hourly time-step. The number of locations is cca. 1000 locations with latitude, longitude coordinates. At each location up to 5 parameters (water level, discharge, rainfall, snow depth, temperature) are measured at regular hourly time-step (gaps in time-series are labeled with a special "nodatavalue" code). Currently I am using a Microsoft SQL Server database with a following main tables:               stations (id, name, latitude, longitude, elevation)     parameters(id, name, units, scalefactor, nodatavalue)     stations_parameters (station_id, parameter_id, start_time, end_time)     waterlevel (station_id, time_utc, value)     discharge (station_id, time_utc, value)     temperature (station_id, time_utc, value)     rainfall (station_id, time_utc, value)     snow (station_id, time_utc, value)      The size of my database is 2 GB. For various reasons mainly due to database size limitation on my webhost and interoperability issues I need to move away from MS SQL Server. I would like to use the NetCDF format because I read that it is suitable for multi-dimensional data, and that it has built-in support for space and time queries. Specifically I need to run the following types of queries very fast: Query 1: For a time-range give me the average, max, min or sum of a parameter for all stations. For example create a map of maximum temperature between 2013-07-25 and 2013-07-26. Query 2: For a time-range, point location and parameter, give me the time- series. For example hourly time-series of temperature at station 777 from 2013-01-01 to 2013-07-25. My questions:   1. Is it possible to use NetCDF for this type of multi-dimensional data? (note that the space is not a regular grid but rather sparse point-locations)   2. Is there any free software tool that has built-in support for visualizing a point time-series NetCDF data file? I have used Panoply and IDV but only with regular NetCDF grids.   3. If yes - What would be the structure of my NetCDF data file (dimensions, variables and their order)?


title: Qgis crash when I attempt to open ECW raster format.

text: I added Ubuntugis unstable repo and installed Qgis, libgdal-ecw-src, libecwj2-3.3 but still do not work. When I attempt to open an ecw file with Qgis it crash and tells me that there are problems with this libraries: > /usr/local/lib/libNCSUtil.so.0.0.0   >  /usr/local/lib/libNCSEcw.so.0.0.0   >  /usr/local/lib/libNCSCnet.so.0.0.0 I am using ubuntu 12.10 and when I do:               gdalinfo --formats|grep -i ecw      it seems to work.               ECW (rw): ERDAS Compressed Wavelets (SDK 3.x)     JP2ECW (rw+v): ERDAS JPEG2000 (SDK 3.x)      What is the problem?


title: Load geographical data administrative boundaries for Europe by levels

text: How can I load all administrative boundaries by levels for Europe in coordinates? I need .osm dates,but not all data inside Europe,just administrative boundaries. I import .osm for Europe in postgis base and by filter 'administrative' read data, but I am not sure is it right query. Thanks in advance! PS. Can I read those data from natural-earth and how?


title: A shapefile / boundary coordinates for Bangladesh divisions

text: Are there any shapefiles on the boundary coordinates for each of the Bangladesh administration divisions? Essentially what I want to achieve is to have the minimum and maximum of the longitudes & latitudes of three divisions: Khulna, Barisal & Chittagong so I can create a rough mask for these regions. I have used this link to obtain the Bangladesh shapefile but does not give me the individual division borders, or that I can see when using ArcGis Explorer. I was thinking, if I had the data for the border coordinates I could extract the min & max values. Thanks


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

text: I'm having trouble with the standard proxy.cgi written in Python for OpenLayers. Here's my code:               #!C:/Python27/python.exe -u          """This is a blind proxy that we use to get around browser     restrictions that prevent the Javascript from loading pages not on the     same server as the Javascript.  This has several problems: it's less     efficient, it might break some sites, and it's a security risk because     people can use this proxy to browse the web and possibly do bad stuff     with it.  It only loads pages via http and https, but it can load any     content type. It supports GET and POST requests."""          import urllib2     import cgi     import sys, os          # Designed to prevent Open Proxy type stuff.          allowedHosts = ['www.openlayers.org', 'openlayers.org',                      'labs.metacarta.com', 'world.freemap.in',                      'prototype.openmnnd.org', 'geo.openplans.org',                     'sigma.openplans.org', 'demo.opengeo.org',                     'www.openstreetmap.org', 'sample.azavea.com',                     'v2.suite.opengeo.org', 'v-swe.uni-muenster.de:8080',                      'vmap0.tiles.osgeo.org', 'www.openrouteservice.org',                     'dcdu.municipiodeslp.gob.mx:8082', 'www.ssn.unam.mx',                     'localhost','localhost:8082', 'localhost:80','127.0.0.1:8082','127.0.0.1']          method = os.environ["REQUEST_METHOD"]          if method == "POST":         qs = os.environ["QUERY_STRING"]         d = cgi.parse_qs(qs)         if d.has_key("url"):             url = d["url"][0]         else:             url = "http://www.openlayers.org"     else:         fs = cgi.FieldStorage()         url = fs.getvalue('url', "http://www.openlayers.org")          try:         host = url.split("/")[2]         if allowedHosts and not host in allowedHosts:             print "Status: 502 Bad Gateway"             print "Content-Type: text/plain"             print             print "This proxy does not allow you to access that location (%s)." % (host,)             print             print os.environ              elif url.startswith("http://") or url.startswith("https://"):                  if method == "POST":                 length = int(os.environ["CONTENT_LENGTH"])                 headers = {"Content-Type": os.environ["CONTENT_TYPE"]}                 body = sys.stdin.read(length)                 r = urllib2.Request(url, body, headers)                 y = urllib2.urlopen(r)             else:                 y = urllib2.urlopen(url)                  # print content type header             i = y.info()             if i.has_key("Content-Type"):                 print "Content-Type: %s" % (i["Content-Type"])             else:                 print "Content-Type: text/plain"             print                  print y.read()                  y.close()         else:             print "Content-Type: text/plain"             print             print "Illegal request."          except Exception, E:         print "Status: 500 Unexpected Error"         print "Content-Type: text/plain"         print          print "Some unexpected error occurred. Error text was:", E      The odd thing is that I am able to receive a GeoRSS feed through the same proxy but not a WFS from GeoServer as shown below: ![enter image description here](http://i.stack.imgur.com/p5sDO.png) ![enter image description here](http://i.stack.imgur.com/WLb7V.png) ![enter image description here](http://i.stack.imgur.com/u6SrH.png) Even though 'localhost:8082' is added to my 'allowedHosts', I receive a '502 Bad Gateway' error. I've looked at many questions on proxy.cgi and a tutorial but I still have no idea what's wrong. I'm running Apache on XAMPP on port 80 while GeoServer is on port 8082. **UPDATE** Here's what my "Workspace" in GeoServer is: ![enter image description here](http://i.stack.imgur.com/uXEpI.png) and the layer is registering as WFS: ![enter image description here](http://i.stack.imgur.com/LKHRv.png) And if I browse to the WFS at: `http://localhost:8082/geoserver/wfs/DescribeFeatureType?version=1.1.0&typename=zac:gcolonia` I receive xml:               <?xml version="1.0"?>     <xsd:schema xmlns:cite="http://www.opengeospatial.net/cite" xmlns:gml="http://www.opengis.net/gml" xmlns:it.geosolutions="http://www.geo-solutions.it" xmlns:mx="mx" xmlns:nurc="http://www.nurc.nato.int" xmlns:sanluis="sanluis" xmlns:sde="http://geoserver.sf.net" xmlns:sf="http://www.openplans.org/spearfish" xmlns:slp="slp" xmlns:test="http://test" xmlns:test_php="http://test_php" xmlns:tiger="http://www.census.gov" xmlns:tlx="tlx" xmlns:topp="http://www.openplans.org/topp" xmlns:xsd="http://www.w3.org/2001/XMLSchema" xmlns:zac="http://localhost:80/Atlas_TLX/zac" elementFormDefault="qualified" targetNamespace="http://localhost:80/Atlas_TLX/zac">     <xsd:import namespace="http://www.opengis.net/gml" schemaLocation="http://localhost:8082/geoserver/schemas/gml/3.1.1/base/gml.xsd"/>     <xsd:complexType name="gcoloniaType">         <xsd:complexContent>             <xsd:extension base="gml:AbstractFeatureType">                 <xsd:sequence>                     <xsd:element maxOccurs="1" minOccurs="0" name="gid" nillable="true" type="xsd:long"/>                     <xsd:element maxOccurs="1" minOccurs="0" name="edo_cve" nillable="true" type="xsd:int"/>                     <xsd:element maxOccurs="1" minOccurs="0" name="regcat_cve" nillable="true" type="xsd:int"/>                     <xsd:element maxOccurs="1" minOccurs="0" name="mun_cve" nillable="true" type="xsd:int"/>                     <xsd:element maxOccurs="1" minOccurs="0" name="zoncat_cve" nillable="true" type="xsd:int"/>                     <xsd:element maxOccurs="1" minOccurs="0" name="loc_cve" nillable="true" type="xsd:int"/>                     <xsd:element maxOccurs="1" minOccurs="0" name="col_id" nillable="true" type="xsd:long"/>                     <xsd:element maxOccurs="1" minOccurs="0" name="col_area" nillable="true" type="xsd:decimal"/>                     <xsd:element maxOccurs="1" minOccurs="0" name="catcol_nombre" nillable="true" type="xsd:string"/>                     <xsd:element maxOccurs="1" minOccurs="0" name="col_geom" nillable="true" type="gml:GeometryPropertyType"/>                 </xsd:sequence>             </xsd:extension>         </xsd:complexContent>     </xsd:complexType>     <xsd:element name="gcolonia" substitutionGroup="gml:_Feature" type="zac:gcoloniaType"/>     </xsd:schema>      I'm beginning to think it's a problem with OpenLayers. Here's how I'm adding it to the map:               gcolonia = new OpenLayers.Layer.Vector("gcolonia", {                                     strategies: [new OpenLayers.Strategy.Fixed()],                                     protocol: new OpenLayers.Protocol.WFS({                                         srsName: new OpenLayers.Projection("EPSG:4326"),                                         url: "http://localhost:8082/geoserver/wfs",                                         featurePrefix:"zac",                                         featureType: "gcolonia",                                         featureNS: "http://localhost:80/Atlas_TLX/zac",                                         geometryName: "col_geom"                                          }),                                     visibility: false                             });      **Newest Update** Browsing to `http://localhost/cgi-bin/proxy.cgi?url=http://www.ssn.unam.mx` -> fine Browsing to `http://localhost/cgi- bin/proxy.cgi?url=http://localhost:8082/geoserver/web` ->               This proxy does not allow you to access (localhost:8082).      Yet, my `proxy.cgi` contains that location as shown above. My `cgi-bin` is located in `xampp/cgi-bin`. I've also cleared my browser cache. Any thoughts?


title: Access traffic data using Google Maps API

text: The following screenshot shows a traffic filter, adopted in the new Google Maps. ![enter image description here](http://i.stack.imgur.com/OsNCV.png) Can we implement the similar function using Google Maps API? It looks like that `google.maps.TrafficLayer` can only display the live traffic data.


title: The scale in Qgis is not setting to feet

text: I am working in QGIS 2.4. I have changed the CRS of all the layers and in the project properties to `EPSG:102739 - NAD_1983_StatePlane_Texas_Central_FIPS_4203_Feet`. I am in print composer trying to create a scale bar and I can not figure out why the scale is so off. I am trying to go from feet to miles and it just creates a really long scale bar. Basically `1 mi = 5280` is what I am going for I have already checked each layer and they are all in state plane 4203. I am not sure what I am missing. I checked other threads about scale in QGIS and I am still at a loss. I have tried flipping the conversion around and its still not working. Edit 2: Sounds crazy is there a difference between Qgis for Windows and for Linux. I loaded the exact same data(just copied and pasted the file with the map and data) to a windows computer. Opened it and everything works fine. Edit: The answer given How to get correct scale values? does not address my issues. I have already done everything in that answer. All of my layers have been resaved with projection and the project projection has been changed. Problems with scale bar units in qgis composer What do I need to change to get scale bar units other than degrees? Why does QGIS scalebar show a wrong scale? ![enter image description here](http://i.stack.imgur.com/YCtd9.png) ![enter image description here](http://i.stack.imgur.com/AVdsF.png) ![enter image description here](http://i.stack.imgur.com/4LqsG.png)


title: Using QGIS to create a table of correspondence between two area layers

text: I have two area layers over a map of Australia, and I want to know if it's possible to use QGIS to create a table of correspondences between them. ![map areas](http://i.stack.imgur.com/dv2GV.png) The Australian Bureau of Statistics provides tables of correspondences between areas, and lists the ratio of one area type that falls into another area type. The example below show the correspondence between Statistical Area 2 and federal electorates: ![ABS correspondence](http://i.stack.imgur.com/QIC4g.jpg) Is there a way to create a table like this using, say, the Field Calculator function in the Attribute Table once I've joined the two layers?


title: Reclassify a rasterList and sum the reclassified value from each Pixel

text: I'm trying to write a script in arcpy which reclassified a List of Rasters. This rasters are made by another tool before so there is already a loop. Actually I get the error 000582 and before this the script could reclassify but overwrote the rasters. In the end I like to count the pixel with the value 1 so that I create one raster with the sum of the value 1 as pixelvalue. How can I do this and what is wrong with the following script? Do I need the tool “plus” at the end?               import arcpy     import os     from arcpy import env     from arcpy.sa import *          arcpy.CheckOutExtension("Spatial")     env.overwriteOutput = True     env.workspace = r"D:\Julia_T\projekt\aufgabe3\testdatenA"     Dir = env.workspace     #folder to list     ascFileList = arcpy.ListFiles("*.asc")     print ascFileList               for ascFile in ascFileList:     #geoprocessing steps     ascFileName = os.path.splitext (ascFile)[0]     print ascFileName     #works     rastFile = ascFileName + "_output.img"     print rastFile     #works     #run the tool     arcpy.ASCIIToRaster_conversion (env.workspace + "\\" + ascFile, env.workspace + "\\" +     rastFile, "INTEGER")     # Process: Projektion definieren     arcpy.DefineProjection_management (env.workspace + "\\" + rastFile,     "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRI    MEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]")     #list raster     rasters = arcpy.ListRasters("*", "IMG")     for ras in rasters:         print "working..." + ras #works         # new name         name = "class_" + ras         # Recalssify         reclassField = "Value"         remap = RemapValue ([["8",1],["2",0],["5",0]])         outReclassify = Reclassify(ras, reclassField, remap, "NODATA")         #save         outReclassify.save (Dir + "\\" + name)          # Execute Plus     #outPlus = Plus(inRaster1, inRaster2)      I'm working with ArcMap 10.2.1. Thanks!


title: How can a column in a feature class be automatically updated as soon as another column is changed?

text: I have a feature class in a file geodatabase and it contains a field called [curing] in short integer type and another field called [IsPopulated] in short integer as well. What I want to achieve is to make [IsPopulated] automatically populated with "0" or "1" depending on whether the [curing] column is filled with a numeric value. "Field Calculator" does not help as it's just a one-off operation and NOT automatic.


title: How to convert unix time to qt time type?

text: How can I convert unix timestamp to qt time type using qgis field calculator? I would like to avoid importing my table into excel and convert the timestamp to qt time type.


title: How to export polygon coordinates to CSV (not WKT)?

text: I am new to QGIS, and I am now trying to export polygon coordinates from a shape file to a csv file, but so far all my attempts have failed. I have found numerous questions and answers about this, but none have solved my problem. I have downloaded the file www.naturalearthdata.com/http//www.naturalearthdata.com/download/10m/cultural/ne_10m_admin_0_countries.zip which is a shape file containing a basic world map. I have imported the shape file into QGIS by adding a vector layer, and then I try to export the polygons with coordinates for each country. The tips that I have found suggest either saving the layer by right-clicking it and then selecting Save as or opening the attribute table and copying it into Excel. If I try these I only get master information with one row for each country in the map. Does anyone know how I can retrieve detail information for each country containing the polygons with coordinates?


title: Possible to use Google Earth to generate "skyline" planes?

text: I was really interested to notice that my city appears to have gotten the 3D Buildings treatment on Google Earth. I am trying to figure out how I might these 3D buildings to generate 'skylines' for a given street. Speaking more abstractly, I think what I want is to generate/export a projection of an arbitrary subset (let's say: 150 > n > 1) of these 3D buildings onto an arbitrary 2D plane of my choosing. More specifically to the goal at hand, one of that 2D plane's dimensions would run normal to the earth's surface, from the surface to some arbitrary height, and the other dimension would run for some arbitrary length of street. (I guess I'll assume for the moment that the street is effectively straight.) More practically speaking, I envision being able to (in the Google Earth app) place a near clipping plane (i.e. take out all the buildings/terrain/etc between the camera and this plane) and a far clipping plane (i.e. take out all the buildings/terrain/etc behind that plane to infinity), then position the camera and export what the camera sees (again, ideally in a variety of ways, like a standard raster of what I see in Google Earth, or an orthographic wireframe projection, or an orthographic hidden-line projection, or... etc.) I see that it's possible to download single models, but it seems like it'd be a long haul to get from that to the desired goal, assuming for the moment I even knew what tool I'd use to do the clipping planes/camera/export once I got the models all downloaded. Anyone have any ideas on how I might go about this?


title: How to compare two string field values within a given string value in ArcGIS10

text: I have a polygon shapefile which displays villages [ **NAMING** ]. On its attribute table I have, in one Field, their names and for each feature and I have its municipality too [ **POLYGON_NM** ]. Now what I want to do is to find if I have features repeated into the same municipality. I wanted to work with ModelBuilder, and I've tried with some scripting code into Field Calculator, but it's a disaster more than success. Here is a screenshot: ![example](http://i.stack.imgur.com/HGowR.png)


title: Halt QGIS on startup

text: I have a custom made login form in Python, which launches on QGIS startup. But before the user has a chance to enter the credentials, the program starts normally, bypassing the form. From then on the user can interact with QGIS like the form isn't even there (which defeats the purpose of the form!! ). Is there a way to halt QGIS startup until the user is authenticated?


title: Problem with multiple feature selection

text: In an attempt to retrieve all features (point geometry) that fall in a selected polygon, I get the error `this.layer.addFeatures is not a function`. My code is as shown below:               function createWFSVector(map, WFS_HOST) {     var population = new OpenLayers.Layer.Vector("Population Vector", {         strategies: [new OpenLayers.Strategy.BBOX()],         protocol: new OpenLayers.Protocol.WFS({             version: "1.1.0",             url: WFS_HOST,             featureType: "population_layer",             featureNS: "pe_ns",             srsName: "EPSG:4326"         })     });          addMapControls(map, population);     map.addLayers([population]);     }          function addMapControls(map, wfsLayer) {     var select = new OpenLayers.Control.SelectFeature([wfsLayer], {         box: true,         multiple: true,         onSelect: addSelected,         onUnselect: clearSelected     });     map.addControl(select);     map.addControl(new OpenLayers.Control.EditingToolbar(select));     select.activate();     }          function addSelected(feature) {     selectedFeatures.push(feature);     }          function clearSelected(feature) {     selectedFeatures = [];     }      I can't seem to spot the problem with my code. Any assistance will be appreciated.


title: Why does QGIS Polygonize ignore the no data value?

text: I'm running into an issue georeferencing a number of PNG image files. The files simply contain an image where each pixel is either transparent or a single color, forming a coverage mask. In the past, when I've georeferenced these images, the 'use 0 for transparency' checkbox masked off the transparent parts of the images and I was able to polygonize the colored parts. However, with a number of images in a recent patch QGIS will simply not mask off the transparent sections. I've opened the files with an image editor and resaved them to no effect, so it's not a formatting problem with the input png. I can somewhat work around this by setting the value of the transparent section to the no data value, thus masking off the image in the raster layer. However, when I polygonize the raster, the no data value still gets polygonized into the .shp file. How do I polygonize the raster and specify only the single band I want to polygonize?


title: Trying to make QT + QGIS working

text: I'm mostly a Java programmer, I worked on C++ many years ago, and knowing this my project manager has asked me to create a little app in Qt that can open a shape file, using QGIS. But I'm encountering some difficulties. I have downloaded the QGIS tutorials of Tim Sutton from GitHub (https://github.com/qgis/QGIS-Code-Examples) and I began working on this: 4_adding_rasters_to_canvas. I have downloaded QGIS with all the sources, version 1.8, and Qt4 (first I tried QGIS 2.0 and QT5, but I have read that there is no compatibility, so I have tried to reproduce the original environment of the tutorials). I have open the project, edited the .pro file to include folders and libraries. I have downloaded OSGEO4W64 (I am on Win7) with all the sources. I have configured the kit under Qt4, using as a compiler Microsoft Visual C++ 10.0 (the only one that doesn't make appear the red error icon in the kits list). I have configured CMake with the path to the CMake 2.8 executable. When I build the project, I see the following errors:               main.obj:-1: error: LNK2019: unresolved external symbol "public: virtual __thiscall QgsApplication::~QgsApplication(void)" (??1QgsApplication@@UAE@XZ) referenced in function _main     main.obj:-1: error: LNK2019: unresolved external symbol "public: __thiscall QgsApplication::QgsApplication(int &,char * *,bool,class QString)" (??0QgsApplication@@QAE@AAHPAPAD_NVQString@@@Z) referenced in function _main      I have done many things for solve this issues, but I want to start again from this point to be sure. My (first for now) question is: I need to compile the QGIS sources under Visual C++ before building the project, because Qt can't see the QGIS dll? Thanks for your time. Bye, Ivan A. **UPDATE** Ok, it seems that I've managed to solve the issue: I've recompiled the 1.8 sources with cmake-gui, and linked the libraries in the build folder. In this way QT sees the libraries. But now I have another issue. There is this code in the tutorial               QgsMapLayerRegistry::instance()->addMapLayer(mypLayer, TRUE);      (mypLayer is a QgsRasterLayer*) and when building it causes the following error:               error: LNK2019: unresolved external symbol "public: class QgsMapLayer * __thiscall QgsMapLayerRegistry::addMapLayer(class QgsMapLayer *,bool,bool)"      I have read that the addMapLayer method is deprecated, so I have replaced it with the method addMapLayers (according to Tim Sutton: http://osgeo- org.1560.x6.nabble.com/Display-raster-layer-on-own-Map-Canvas-td5002935.html ). So I used:               QgsMapLayerRegistry::instance()->addMapLayers(QList<QgsMapLayer*>() << mypLayer);      and the error slightly change but doesn't go away...               error: LNK2019: unresolved external symbol "public: class QList<class QgsMapLayer *> __thiscall QgsMapLayerRegistry::addMapLayers(class QList<class QgsMapLayer *>,bool,bool)"      It's the proper syntax (I found examples on the internet)? There is another way to do this?


title: Super large scale maps with Leaflet, how?

text: I have a map/architectural drawing of two floors of a huge building with a lot of rooms. Now I'd like to make an interactive slippy map using preferably Leaflet showing all the rooms and use markers/polygons showing what happens in each room and who's sitting where, etc. These markers should quite easily be editable. Does anyone have any ideas of how to make this happen? First I was thinking about just pretending house covers the whole world and then use standard web mercator with `generate_tiles.py`. Another idea is to use zoomify together with Leaflet. The problem with both methods is that it's difficult to get right coordinates. Best would be to use some kind of a local projection.


title: PGRouting 2.0 for RHEL5

text: Is there a PGRouting 2.0 rpm package available for Redhat Enterprise Linux 5.8 and has anyone managed to install it successfully?


title: LINECAP and PATTERN on PHP mapscript

text: I am trying to create a railroad line pattern by using LINECAP and PATTERN. The manual says it can be done by using both paramaters on STYLE inside mapfile. Unfortunately, i couldn't find any method or option to do that on php mapscript. Closest i can find is by using setpattern(). But that method is inside symbolObj. Any idea how to do LINECAP and PATTERN with php mapscript ? Thanks. :)


title: osm2pgsql: Import data in custom projection

text: I have osm-file and I want to import it to the postgresql via osm2pgsql with custom projection EPSG:20022. I look to the help and found out any keys like -E EPSG: but when I try to do with it I get the data in SRID=4326. Is there a way how to fix it?


title: Geodjango insert statement: psycopg2.IntegrityError: violates check constraint

text: Why is the following INSERT statement throwing this error?               INSERT INTO domes_place (placeid, structidx, oscode) VALUES ('10', '1', 'POINT(51.121282945772 1.3138398501222)');     psycopg2.IntegrityError: new row for relation "domes_place" violates check constraint "enforce_srid_oscode"      This is the table definition in postgres:                          Table "public.domes_place"        Column   |          Type          | Modifiers      ------------+------------------------+-----------      placeid    | integer                | not null      structidx  | integer                |        oscode     | geometry               |      Check constraints:         "enforce_dims_oscode" CHECK (st_ndims(oscode) = 2)         "enforce_geotype_oscode" CHECK (geometrytype(oscode) = 'POINT'::text OR oscode IS NULL)         "enforce_srid_oscode" CHECK (st_srid(oscode) = 4326)      Isn't this POINT in WGS84 - i.e. in SRID 4326?


title: How do I create a static popup that displays a local web page?

text: I am using OpenLayers 2.11 and would like to place a marker on my map which, when clicked on, would open a popup which then loads a local web page. I would like the popup to be located in the lower, left hand corner of the screen regardless of zoom level or map scrolling by the user. I have been unable to locate an example of how to perform this. Any assistance is appreciated. Thank you!


title: GeoServer raster tiles are individually contrast-balanced: how to turn this off?

text: I'm running GeoServer locally as a servlet (in JBoss), and I'm serving a GeoTIFF file through the WMS protocol. In a webpage running an OpenLayers map, I have the GeoTIFF-WMS layer overlaid above an OpenStreetMaps base layer. The problem is that it looks like this: ![GeoTIFF over OpenStreetMap](http://i.stack.imgur.com/mnfR5.png) The GeoTIFF is clearly being tiled (good) and each tile is normalized to a given brightness level (bad). When I zoom in, I can see a similar effect taking place on the smaller tiles, and I know that the original GeoTIFF doesn't have this blockiness: it's a diagonal strip of 16-bit greyscale surrounded by zeros that ought to be transparent. How do I turn off this brightness-normalization feature (or what is it called, so that I can Google it)? Even better, how do I tell GeoServer that pixels with zero values ought to be transparent?


title: Mapping Trails with GPS - Hardware Suggestions

text: I have been hired to map about 100km of XC ski trails. I am super proficient in GIS but have not done any detailed GPS work since college when I worked with some sort of Trimble explorers. This is a complex network of trails and each segment needs to be a unique feature. I own a Garmin etrex 30. My guess is that this will not be a great unit for this sort of work. I am thinking of buying a used Trimble GeoExplorer and an antenna that I can attach to my hat so that the GPS can be in a hip pack. I do want my polylines to be accurate and if I remember correctly the Trimble points can be nicely corrected after collection. Any other suggestions for hardware? Do I even need a device this sophisticated?


title: How can I split polygons in the Natural Earth 10m land shapefile?

text: This problem, that I have also been struggling with, in QGIS 1.7.3 as well. Here are the steps to reproduce the problem: **1)**Download Natural Earth 10m land shapefile from original website. * _2)_ *Open the file as vector shapefile in QGIS. * _3)_ *On layer browser, right click and enable "Toggle Editing" * _4)_ *Press button "Select single feature by click", then click on the piece of land (becomes yellow) * _5)_ *Click "Split Feature" scissors button. At this point the land layer disappears, (the map display goes white, but land layer is still on layer browser though) and the "Toggle Editing" is disabled without me doing anything. As soon as I re-enable editing, the map comes back. Other features like drawing polygons function normally. Thanks in advance.


title: Using iface.showAttributeTable() on QgsVectorLayer retrieved by getObject()

text: I'm unsure if i'm properly using the processing.getObject() function. I'm creating a python processing script for QGIS 2.0 and 2.2, but can't seem to get an attribute table to open programmatically. The code below works fine in the Python Console, but causes QGIS to crash when I tried as a Processing Script. The only difference was that the layer is called by QgsVectorLayer() in the Python Console and processing.getObject() in the Processing script, thus my reasoning behind my question.               from qgis.core import *     from PyQt4.QtCore import *     from PyQt4.QtGui import *     from processing.core.VectorWriter import VectorWriter     import numpy     import processing     import csv     import qgis.utils          #Here we define the input and outputs     #====================================     ##[Example scripts]=group     ##Couche_compteurs_individuels=vector     ##Couche_compteurs_generaux=vector     ##Couche_canalisations_sur_compteur=vector          #And here is the body of the algorithm     #=======================================          #input layers values are always a string with its location.     #That string can be converted into a QGIS object (a QgsVectorLayer in this case))     #using the processing.getobject() method     indivLayer = processing.getobject(Couche_compteurs_individuels)     grosLayer = processing.getobject(Couche_compteurs_generaux)     resLayer = processing.getobject(Couche_canalisations_sur_compteur)          CG_IDlist = []          def AEPstats():         Gros_Feat = grosLayer.getFeatures()         grosIDselection = []         grosidNull = False         for feat in Gros_Feat:             ID_idx = grosLayer.fieldNameIndex('ID_CG')             #selection list of null values             if feat.attributes()[ID_idx] == NULL:                 grosIDselection.append(feat.id())                 grosLayer.setSelectedFeatures(grosIDselection)                 grosidNull = True             #Putting existing attributes in list             else:                 CG_IDlist.append(feat.attributes()[ID_idx])         #Open msg box, attribute table with selected items, and end function         if grosidNull == True:             QMessageBox.warning(                 None,                 'Probleme de saisie',                 "Completez identifiant(s) de reference ('CG_ID') de Compteurs Individuel")             qgis.utils.iface.showAttributeTable(grosLayer)             return      Thanks in advance


title: Convert vector data with irregular points (no grid) to raster with QGIS?

text: I am aware that how to convert a vector layer to a raster layer is a quite frequent problem and that there are many questions and answers dealing with it. However, none was helpful concerning my specific problem, so I hope it's okay to ask again. I have a huge dataset consisting of fish catch data - different species and the location where they were caught. The original data is a csv file, but I'm working with it in QGIS and have created a vector file from the spreadsheet. Now I want to use another program (Zonation) that only accepts raster files to calculate species richness in certain areas - is there a possibility to create a raster from a vector layer with irregular points? I have tried the Rasterize tool, but it only gives me a solid rectangle and all individual data points are lost. I'm working with CRS WGS 84 and I have tried to set the cell size to 0.4, but still it's only one solid square. I've never worked with raster files before and I'm very grateful for your help - thanks!


title: Arcgis CopyFeatures tool extremely slow when exporting to SDE

text: I have about 10 points in a file gdb that I am attempting to copy over to an sde feature set. In arcpy and in model builder it takes about a minute to execute. Do you know of any tips to speed this up? In SSMS I can query the table almost instantaneously, so there must be a bottleneck somewhere that I am not aware of, since the network is close to zero latency. In the progress bar, it reads "Writing Features 1" for the majority of the time, then seems to complete the rest of the task almost immediately. If it is really supposed to be this slow, I would at least like to have some idea of what it is actually doing during this time. The file is only a few kilobytes total, so I assume this cannot be normal. **Update** : I was trying this operation over a site-to-site VPN. It is a fast connection, so there is still no real reason it should take a full minute to upload 10 points. I tried this on a local connection (local to the SQL Server), and the same operation took 2 seconds. It seems that CopyFeatures has some sort of bug when working through a VPN connection. Any ideas on how I could find a workaround (I need to be able to copyfeatures over the network with a VPN connection)?


title: Windows 7 ArcGIS 9.3 with TOPO!Pro Runtime Error

text: A new, clean install of ArcGIS on Windows 7 SP1 and National Geographic TOPO! Pro Suite -runtime error "Runtime Errror (UpdateGroupLayer)#-2147217384 Failed to open raster dataset" TOPO! Pro was designed with XP not Win 7. unable to get Tech Support (closed last year - bummer) Help the fireguy please!!


title: Setting ProductCode for python script (Citrix environment)

text: I have set up a python script in ArcToolbox for our users. However when running it with ArcView 10.0 licensing the script fails with: > ERROR 000824: The tool is not licensed I think it's an Identity_analysis command triggering this error. I've googled the topic and it seems there was a SetProduct(code) command, but I can't seem to get it to work. Can someone provide the syntax to get this running, or is it a limitation running ArcGIS in a Citrix environment?


title: Is a proxy required to use OpenLayer's GetFeatureInfo?

text: My domain is hosted on HostGator and apparently they don't allow the use of proxies. I have a map that has a WMS point feature hosted from my GeoServer using OpenLayers. I want to be able to select a point and get information from it. Can I do this without using a proxy?


title: What is the relationship between Map Service Instances and SDE connections?

text: We have a GPS project based on ArcGIS Mobile 10.0. There are two web map services and each were kept as Nonpooled with Minimum instances 15 and maximum instances 60 by original developer. The database is ORACLE and connection is made through SDE. Recently whole enterprise moved to 10.1 except GPS project which has some crashing issues with 10.1. But that led hosting of this GPS project on an old Test server. Two incidents happened in last two weeks when the connections for GPS projects reached alarming high in few minutes (122 and 58). First time it affected the SDE server and caused total chaos- we had to kill the processes and bounce the SDE server back. Second time, we restarted the services before things out of control. I investigated the first time - only 20 GPS devices were connected that time. I am new to ArcGIS mobile and don't know SDE in detail(It is still under DBA control and I can only investigate the appserver). How will the number of SDE connections being used one time transforms to map services? I want to go to management with the reasons and recommendations.


title: Interpolation outside the contour of a set of points

text: I have a field contour in a shape layer. I have another layer compound by a set of points which are the locations where soils samples have been taken within the field. I have in the attribute table of this second shape layer the Phosphorus content of the soil in those points. As I have said, the contour of the set of points is smaller than the field contour, being the contour of the set of points enclosed within the field contour. I have interpolated this set of points by creating a TIN in order to have a value of Phosphorus content for any point in the field. But through this procedure, the area between the contour of the set of points and the contour of the field is excluded from the TIN. I have also tried the Kriging method (this method allows to enter the field contour shape as the Processing Extent in the Environments section, being the referred problem solved) but through this interpolation method a new problem arises: the sample points are lying over classes which value ranges don’t include the real values of the sample points. Then, the original data is in some way adulterated. I guess the IDW method isn’t good for this purpose: there isn’t any point in the field contour, and thus the value will go down towards the field contour in an unreal way due to this fact. Any idea or solution for this issue (interpolating the area excluded from the TIN)?


title: Rasterize difference between QGIS and FWTools

text: I am wanting to use FWTools to rasterize a large number of Shapefile polygons. I have used the Rasterize tool in QIS and it works great. But when I use the generated code in FWTools the process does not work. Here is the generated code:               gdal_rasterize -a COVTYPE -ts 5000 3705 -l LCC001k_dis     C:/GlobalData/Vectors/LCC/LCCRasters/LCC001k_dis.shp     C:/GlobalData/Vectors/LCC/LCCRasters/LCC001k.tif      Unfortunately I am not a programmer, so I am not able to create my own script to perform this task. Can anyone assist me with why this code doesn't work in FWTools, or what changes I need to make to utilize the Rasterize functionality.


title: Transparency in ArcMap Graphic Elements

text: Is it possible to apply transparency to graphic elements in an ArcMap layout? I'm specifically refering to a neatline -- a box containing my legend, source info, scale bar, etc.


title: Workflow process to style vector data using GeoServer

text: I'm starting to get a better understanding of map-servers in general, but I am still unclear on how best to use GeoServer. At the moment, I use web-hosting to store the web-browser related parts of my site, and `php` code to connect to a postgresql/postgis database. All of my data is stored in a desktop server, which I chose to use as it is a very low cost way of storing a large amount of data.               1] web-browser <--> database         This works fine also, as I do not need to style vector data; I am only returning numerical results with geojson vector overlays. However, I would like to incorporate GeoServer into this workflow so that I can apply better cartographic styles to data in the future. I understand the work flow should be similar to the following:               2] web-browser <-->  GeoServer <--> database      My difficulty in understanding this piece, is how to connect GeoServer into the process. Can I run GeoServer on my desktop and connect to it to the browser via `php`, similar to the database approach? Or does GeoServer need to be installed on the web-host and managed there?


title: Fast Offline OSM Routing Library to Replace MoNav (.NET)

text: I'm looking for a .NET replacement to MoNav. I'm currently calling a C++ MoNav DLL via COM and it works great--it's fast and memory efficient, and can both load itself and route across the country almost instantaneously. A pre- processed graph for a 5GB OSM file is only ~400MB, which is another plus as I do not need to distribute the source OSM files. However, the COM dependency is a slight issue with certain customers' setups. Its turn-by-turn directions leave a lot to be desired as well. So I was wondering if anyone knows of a .NET replacement? One that preferably used Contraction Hierarchies, has a low memory footprint, and can reliably route across the United States in under a second. That may be asking a lot, I know. My searches have only turned up OsmSharp thus far, and while it looked promising initially, in practice the load times were abysmal for even pre- processed CH graphs (though I am still playing with it and tweaking parameters to try and resolve this). Just wanted to know if anyone knew of a suitable replacement before I attempt to write a C++/CLI wrapper for the original MoNav library and spend a few days (weeks?) on a better turn-by-turn system. So far I haven't found one. Thanks!


title: maintain different style for each feature in same layer

text: on feature selection, I am trying to change feature style of selected feature.so that assign different style to each selected feature that exist in same layer. I have done this successfully but when I zoomed the map, feature redraw and assign the last symbology to all selected features. Please resolve my issue Here is my code: var style = OpenLayers.Util.extend({}, OpenLayers.Feature.Vector.style['default']);                   style.strokeColor = color;        style.strokeWidth = width;        selectItem.style = style;        polygonLayer.drawFeature(selectItem);


title: Edit a classified grid file (.grc) to improve accuracy

text: I have a .grc file (classified grid file) which contains some clutter information, forest, sea, woodland, etc... I need to alter the data in some area to improve accuracy. Basically a larger forest or possibly create a new value like "tall tree forest", etc.... Is there a way to do this? Thanks in advance


title: Cannot install node-mapnik on Windows 2008 x64

text: I have been spending the last few days trying to get node-mapnik running but with no success. Initially inspired by Wilhelm Berg's demo (https://www.youtube.com/watch?v=GlTaRpg3CbQ) which renders pure vector maps using node-mapnik I was interested in getting this running on a Windows Server. I have worked through all of the StackExchange issues on these topics including Wilhelm's checklist for getting this working on Windows 7 x64 but I am at a dead end and my knowledge is really limited in being able to further diagnose the problem. My environment is:   * Windows 2008 Server R2 x64   * Visual Studio Express 2012 Desktop Edition (32 & 64 bit)   * Node.js v0.10.26 (64 bit)   * Python27 (32 bit)    I have taken the following steps:   1. Installed Mapnik 2.2.0 SDK downloaded from git and placed in C:\apps\mapnik-v2.2.0      2. Added to PATH: C:\apps\mapnik-v2.2.0\lib C:\apps\mapnik-v2.2.0\bin      3. Installed Python v2.7 (32bit) downloaded and placed in C:\apps\python27      4. Added to PATH: C:\apps\python27 C:\apps\python27\scripts      5. Added to PYTHONPATH: C:\apps\Python27\Lib; C:\apps\Python27\DLLs; C:\apps\Python27\Lib\lib-tk; C:\apps\mapnik-v2.2.0\python\2.7\site-packages;      6. Rebuild node-gyp using VS 2012 node-gyp rebuild --msvs_version=2012      7. Downloaded and installed protobuf (32 bit) and placed in c:\apps\protobuf      8. Added to PATH: c:\apps\protobuf      9. Installed protobuf using NPM npm install protobuf --msvs_version=2012      10. Downloaded Cairo using git git clone git://anongit.freedesktop.org/git/cairo      11. Copy Cairo files to mapnik: Copy cairo\src\cairo-svg-surface.c to C:\apps\mapnik-v2.2.0\include Copy cairo\src\cairo-svg.h to C:\apps\mapnik-v2.2.0\include Copy cairo\src\cairo-svg-surface-private.h to C:\apps\mapnik-v2.2.0\include      12. Make proj4 epsg files available:  SET PROJ_LIB=C:\apps\mapnik-v2.2.0\share\proj      13. Attempt to install node-mapnik npm install node-mapnik --msvs_version=2012    I receive the following errors: > node-pre-gyp ERR! Source compile required: 404 status code downloading > tarball ... > > mapnik_grid.cpp c:\users\administrator\node_modules\node- > mapnik\build\release\obj\global_interm ediate\vector_tile.pb.h(9): fatal > error C1083: Cannot open include file: 'google/protobuf/stubs/common.h': No > such file or directory (C:\Users\Administrator\n ode_modules\node- > mapnik\build\Release\obj\global_intermediate\vector_tile.pb.cc) > [C:\Users\Administrator\node_modules\node-mapnik\build\mapnik.vcxproj] > c:\users\administrator\node_modules\node-mapnik\src\utils.hpp(11): fatal > error C1083: Cannot open include file: 'mapnik/unicode.hpp': No such file or > director y (..\src\mapnik_geometry.cpp) > [C:\Users\Administrator\node_modules\node-mapnik \build\mapnik.vcxproj] > mapnik_layer.cpp mapnik_cairo_surface.cpp > c:\users\administrator\node_modules\node- > mapnik\src\mapnik_featureset.hpp(6): f atal error C1083: Cannot open include > file: 'mapnik/datasource_cache.hpp': No s uch file or directory > (..\src\mapnik_featureset.cpp) [C:\Users\Administrator\no de_modules\node- > mapnik\build\mapnik.vcxproj] mapnik_image.cpp > C:\Users\Administrator\node_modules\node-mapnik\node_modules\mapnik-vector- > tile \src\mapnik3x_compatibility.hpp(4): fatal error C1083: Cannot open > include file : 'mapnik/version.hpp': No such file or directory > (..\src\mapnik_map.cpp) [C:\U sers\Administrator\node_modules\node- > mapnik\build\mapnik.vcxproj]   >  mapnik_datasource.cpp C:\Users\Administrator\node_modules\node- > mapnik\node_modules\mapnik-vector-tile \src\mapnik3x_compatibility.hpp(4): > fatal error C1083: Cannot open include file : 'mapnik/version.hpp': No such > file or directory (..\src\mapnik_color.cpp) [C: > \Users\Administrator\node_modules\node-mapnik\build\mapnik.vcxproj]   >  node_mapnik.cpp ..\src\mapnik_grid_view.cpp(9): fatal error C1083: Cannot > open include file: 'm apnik/grid/grid.hpp': No such file or directory > [C:\Users\Administrator\node_mo dules\node-mapnik\build\mapnik.vcxproj] > mapnik_projection.cpp C:\Users\Administrator\node_modules\node- > mapnik\node_modules\mapnik-vector-tile \src\mapnik3x_compatibility.hpp(4): > fatal error C1083: Cannot open include file : 'mapnik/version.hpp': No such > file or directory (..\src\mapnik_layer.cpp) [C: > \Users\Administrator\node_modules\node-mapnik\build\mapnik.vcxproj]   >  mapnik_feature.cpp ..\src\mapnik_grid.cpp(9): fatal error C1083: Cannot > open include file: 'mapnik /version.hpp': No such file or directory > [C:\Users\Administrator\node_modules\n ode-mapnik\build\mapnik.vcxproj] > mapnik_palette.cpp ..\src\mapnik_image_view.cpp(11): fatal error C1083: > Cannot open include file: 'mapnik/color.hpp': No such file or directory > [C:\Users\Administrator\node_modu les\node-mapnik\build\mapnik.vcxproj] > mapnik_vector_tile.cpp c:\users\administrator\node_modules\node- > mapnik\src\utils.hpp(11): fatal error C1083: Cannot open include file: > 'mapnik/unicode.hpp': No such file or director y > (..\src\mapnik_cairo_surface.cpp) > [C:\Users\Administrator\node_modules\node-m apnik\build\mapnik.vcxproj] > mapnik_expression.cpp   >  mapnik_memory_datasource.cpp ..\src\mapnik_memory_datasource.cpp(3): fatal > error C1083: Cannot open include file: 'mapnik/version.hpp': No such file or > directory [C:\Users\Administrator\n ode_modules\node- > mapnik\build\mapnik.vcxproj] c:\users\administrator\node_modules\node- > mapnik\src\utils.hpp(11): fatal error C1083: Cannot open include file: > 'mapnik/unicode.hpp': No such file or director y (..\src\mapnik_feature.cpp) > [C:\Users\Administrator\node_modules\node-mapnik\ build\mapnik.vcxproj] > ..\src\mapnik_image.cpp(7): fatal error C1083: Cannot open include file: > 'mapni k/color.hpp': No such file or directory > [C:\Users\Administrator\node_modules\no de-mapnik\build\mapnik.vcxproj] > C:\Users\Administrator\node_modules\node-mapnik\node_modules\mapnik-vector- > tile \src\mapnik3x_compatibility.hpp(4): fatal error C1083: Cannot open > include file : 'mapnik/version.hpp': No such file or directory > (..\src\mapnik_datasource.cpp ) [C:\Users\Administrator\node_modules\node- > mapnik\build\mapnik.vcxproj] C:\Users\Administrator\node_modules\node- > mapnik\node_modules\mapnik-vector-tile \src\mapnik3x_compatibility.hpp(4): > fatal error C1083: Cannot open include file : 'mapnik/version.hpp': No such > file or directory (..\src\mapnik_projection.cpp ) > [C:\Users\Administrator\node_modules\node-mapnik\build\mapnik.vcxproj] > C:\Users\Administrator\node_modules\node-mapnik\node_modules\mapnik-vector- > tile \src\mapnik3x_compatibility.hpp(4): fatal error C1083: Cannot open > include file : 'mapnik/version.hpp': No such file or directory > (..\src\mapnik_palette.cpp) [ C:\Users\Administrator\node_modules\node- > mapnik\build\mapnik.vcxproj] C:\Users\Administrator\node_modules\node- > mapnik\node_modules\mapnik-vector-tile \src\mapnik3x_compatibility.hpp(4): > fatal error C1083: Cannot open include file : 'mapnik/version.hpp': No such > file or directory (..\src\node_mapnik.cpp) [C:\ > Users\Administrator\node_modules\node-mapnik\build\mapnik.vcxproj] > c:\users\administrator\node_modules\node-mapnik\src\utils.hpp(11): fatal > error C1083: Cannot open include file: 'mapnik/unicode.hpp': No such file or > director y (..\src\mapnik_vector_tile.cpp) > [C:\Users\Administrator\node_modules\node-map nik\build\mapnik.vcxproj]


title: how to calculate the average temperature of each band of a .grd file?

text: I am having a GRD file in which the temperature at different longitude and latitude are given in a 32x35 array for one day(band),i am having 365 such bands and want to calculate the average of all individual bands.


title: Calculation to convert db data to Lat/Lng degrees

text: In a database I have the following Lat/Lng data. I wish to convert it to a readable format.               177514245:34726053  =>  S 32 7.3947, E 115 51.4245     177513828:34726948  =>  S 32 7.3052, E 115 51.3828     177508930:34829424  =>  S 31 57.0576, E 115 50.8930     177508084:34829842  =>  S 31 57.0158, E 115 50.8084     177504234:34831118  =>  S 31 56.8882, E 115 50.4234     177504474:34831039  =>  S 31 56.8961, E 115 50.4474     177507746:34829955  =>  S 31 57.0045, E 115 50.7746     177504623:34830989  =>  S 31 56.9011, E 115 50.4623     177507014:34830191  =>  S 31 56.9809, E 115 50.7014     177506702:34830295  =>  S 31 56.9705, E 115 50.6702      I was given the formula for converting the db data o the readable format but lost the information after a hhd failure. The number before the ':' becomes the Lng and the second the Lat. I don't want to through any red herrings into confuse anyone, but this is what I can remember. Apparently this is a common format. The raw Lat and Lng number was subtracted by a constant. I can't remember if they were the same, and then divided to give degrees and minutes. But take that all with a grain of salt. Thanks in advance, Simon   1. I have no idea. The data in the db is the 18 varchar figure. The OEM software represents the coordinates in the format as shown above.   2. Have been Google ing it for 3 days and have found nothing. Can't remember the site I posted the question on about 3 or 4 years ago so I could go back and read it again. At the time I was told it was a pretty common method of storing coordinates.   3. The db is sqlite, I use RAD Studio. The OEM GUI that displays the data doesn't allow me to access the data through it. But I can read/write the db itself and I wish to map the coordinates. Hence the need to be able to make sense of how:               177514245:34726053      becomes               S 32 7.3947, E 115 51.4245


title: How to get for a given pixel the values from different overlapped rasters?

text: I have a few rasters representing different soil properties, like soil pH and sand content. The rasters overlap perfectly, have the same dimensions and same coordinate reference system. I would like to be able to click on any pixel of my map and get the pH and sand values for this location. The "Identify Features" tool seems to return only the value of the selected raster, pH or Sand. It doesn't give the values of each rasters displayed. Are they any equivalent of the "Identify Features" tool that work on multiple rasters? I am QGIS 2.2.0 user on a Windows computer.


title: how to include part of a map extent in layout view

text: I have a map of Scotland and within the map layout I'd like to include Shetland within a box at the top so that I don't have to minimise the full extent of the map - I was wondering what the best way to do this would be. Previously i created a duplicate layer in the data view and zoomed it to only fit in Shetland in the layout view but this seemed a clunky way of doing this? Anyone have any advice?


title: Google earth incorrectly rendering coordinates?

text: My circles are rendering funky on the poles, and along the equator. ### Here are my coordinates:               158.05553767833,-78.151355482613     152.68654993135,-76.924955252783     147.72303124605,-75.713526831152     143.06274497048,-74.524895329406     138.63332012974,-73.365493766502     134.38233416968,-72.240817945152     130.27100948566,-71.155701792104     126.27011840531,-70.11448759429     122.35726277966,-69.121132703464     118.51502708319,-68.179276670913     114.72969830909,-67.292283062004     110.9903610723,-66.463264689784     107.28824575653,-65.695097793865     103.61625023098,-64.990428769143     99.968582435602,-64.351675871499     96.340488248312,-63.781027589248     92.728040190375,-63.280438894499     89.127969904661,-62.851626275065     85.537532303713,-62.496062233228     81.954392666042,-62.214969784592     78.376530286024,-62.009317375003     74.802153890611,-61.879814541074     71.229625140461,-61.826908561053     67.657387272327,-61.85078227176     64.083896401248,-61.951353160458     60.507553237337,-62.128273775806     56.926633008616,-62.380933438     53.339211222379,-62.708461163863     49.743082525748,-63.109729656798     46.135669299876,-63.583360142432     42.513915669083,-64.127727755403     38.87416120926,-64.740967096403     35.211986616615,-65.420977473287     31.522020664857,-66.165427203563     27.797693492691,-66.971756167889     24.030914923935,-67.837175533906     20.211647006764,-68.758663165469     16.327325438369,-69.732952609442     12.362061986503,-70.756512568067     8.2955243178175,-71.825512172637     4.1013320161292,-72.935764736584     -0.25528729388528,-74.082638186513     -4.8209969442254,-75.260912557308     -9.6596917885399,-76.464550923909     -14.860832875137,-77.686324272452     -20.55272215938,-78.917181740079     -26.923940644354,-80.145162642771     -34.257840694438,-81.353462558233     -42.985222847202,-82.516927422837     -53.749206371382,-83.595765562688     -67.414565058046,-84.525376429557     -84.75320323954,-85.205991865779     -105.36439351601,-85.51452525361     -126.6349523656,-85.372691897859     -145.37174704566,-84.818624018911     -160.36597119513,-83.969463910769     -172.11759253603,-82.937731923535     178.47360903516,-81.800519081463     170.67339362215,-80.605603613521     163.97799351745,-79.38273857854     158.05553767833,-78.151355482613      ### And here is the result: ![Result](http://i.stack.imgur.com/QskM9.jpg) **What could be causing this and how can can I fix it?**


title: How to update a CartoDB table at a certain time of day

text: I have a table in CartoDB 'inland' that I would like to have automatically update each day at 8pm PST. I know that PostgreSQL has no internal function to run an UPDATE COLUMN function at a certain time of day. Is there a way to use Cron or PgAgent to do this? Would I need to setup CartoDB on my local machine instead of only using the webservice?


title: arcgis python window VS pywin32

text: I have problems while using the ZonalHistogram function in Python. It seems to give strange results, while, with the same code run in arcgis' python window everything works fine. I am trying to study the snow cover per elevation ranges, thus including in the function a MODIS image (showing snow cover) as the input raster, and a reclissified DEM as the input value raster. I tried to modify the reclassification many times, but still not solving the problem... Moreover, inverting the input raster with the input value raster (meaning, inverting the function parameters), everything works fine, but it's a mess since I'm working on thousands of images and having the MODIS not always the same number of columns (while being the DEM always the same). Hope someone could help me. Here is my script, so that you can have a look at it in order to better understand my request:               import arcpy, os, sys, shutil          newpath= r'G:\\Chile_2012\\MODIS10A2.V005\\provvisoria\\Zonal_Histogram\\temporary'     if not os.path.exists (newpath):         os.makedirs (newpath)          arcpy.env.workspace= 'G:\\Chile_2012\\MODIS10A2.V005\\Cobertura_nival\\Cuencas_Copia\\campione'     arcpy.env.overwriteOutput= True     from arcpy.sa import *          listRaster= arcpy.ListRasters ()     elev_ranges= 'G:\\Chile_2012\\MODIS10A2.V005\\DEM\\elev_ranges'     out_dir= 'G:\\Chile_2012\\MODIS10A2.V005\\provvisoria\\Zonal_Histogram\\'          try:         for raster in listRaster:             arcpy.CheckOutExtension ('Spatial')             ZonalHistogram (raster, 'Value', elev_ranges, newpath+ '\\'+ raster)              arcpy.ClearEnvironment ('workspace')         arcpy.env.workspace= newpath         listTable = arcpy.ListTables ()         print listTable              arcpy.Merge_management (listTable, out_dir+ 'Rio_Aconcagua_Merge.dbf')              shutil.rmtree (newpath)          except:         print arcpy.GetMessages ()      I would add that my idea would be to work with layers (by means of functions that could transform rasters in layers) since the table returned in arcgis has a "LABEL" field (coming from the classification in symbology) which is the thing I would like to include in the final table and which is not included while working out of arcgis (in pywin), in order to know on what I am working on.


title: Is there an exhaustive, searchable catalog of all GIS web services (e.g. OGC WMS/WFS, REST, SOAP)?

text: Many of us have struggled to get a list of GIS services that are available on the internet. Specifically map services that expose OGC based interfaces like WMS, WFS or Web Service interfaces like REST, SOAP etc. For instance, just for oceanic and marine data there are several online catalogs such as   1) FDGC   2) NOAA NowCoast   3) IMO Mest   4) NASA If you start with weather, transportation, demographics, agriculture, this list becomes huge! These sites are not easy to use, don't have the most updated info and don't have a good taxonomy/ classification of services towards easy discovery of such services. Is there a better attempt out there to create an exhaustive, searchable catalog of all such services? I would like to know their usage rights, free or pay etc. Uptime, performance etc. would also be nice to know.


title: Reprojecting Natural Earth raster into natearth projection with gdalwarp causing black squares?

text: I'm trying to reproject the Natural Earth 2 raster layer (http://www.naturalearthdata.com/downloads/10m-raster-data/10m-natural- earth-2/) into the natural earth projection (http://en.wikipedia.org/wiki/Natural_Earth_projection) using gdalwarp. The whole thing looks pretty good except I get 4 black squares by the edges of the map as seen below. gdalwarp options used [gdalwarp included in QGIS: GDAL 1.11.0, released 2014/04/16]: `gdalwarp -overwrite -s_srs EPSG:4326 -t_srs "+proj=natearth +wktext" -r lanczos -of GTiff NE2_HR_LC_SR_W_DR.tif reprojected.tif ` Any ideas? ![natearth reprojection with black squares](http://i.stack.imgur.com/p5wTs.png)


title: How to programmatically work with featureclasses from multiple geodatabases?

text: I have produced a Python script which creates multiple buffers of specific (listed) featureclasses within a geodatabase. This is achieved by running each feature through a `for` loop. I would like to know if/how it would be possible to list featureclasses from more than one geodatabase? Currently I am using `env.workspace` to specify the single geodatabase. Listing the features and then running them through the 'for loop'. A possible problem might be that many of the featureclasses in the various geodatabases have the same name.               import arcpy     import os.path     import tempfile     from arcpy import env          # sets the workspace to my gdb     arcpy.env.workspace = "C:/GIS Home/Mining Features (MATTHEW).gdb"          # Output workspace and temp location is defined     output_workspace = "C:/GIS Home/project_1"     temp_workspace = tempfile.mkdtemp()          # list of features to buffer in gdb     fcList = ["Adit_Metalliferous", "Shaft_Metalliferous", "Quarry_Metalliferous", "Dump_Metalliferous", "Engine_House", "Goffan_Working", "Lode_Surface_Outcrop", "Lode_Undefined_Elevation", "Mine_Level", "Mine_Raise", "Portal", "Reported_Subsidence", "Secured_Feature", "Tunnel"]          #create empty lists for output     l_250 = []     l_500 = []     l_750 = []          for featureClass in fcList:         b250 = os.path.join(temp_workspace, featureClass + "_b250.shp")  # make temp output file         b500 = os.path.join(temp_workspace, featureClass + "_b500.shp")  # make temp output file         b750 = os.path.join(temp_workspace, featureClass + "_b750.shp")  # make temp output file         arcpy.Buffer_analysis(featureClass, b250, 250, "", "", "ALL")   # carry out buffer         arcpy.Buffer_analysis(featureClass, b500, 500, "", "", "ALL")   # carry out buffer         arcpy.Buffer_analysis(featureClass, b750, 750, "", "", "ALL")   # carry out buffer         l_250.append(b250)  # append the buffer to the l_250 list         l_500.append(b500)  # append the buffer to the l_500 list         l_750.append(b750)  # append the buffer to the l_750 list          merge_250 = os.path.join(temp_workspace, "250_merge.shp")     arcpy.Merge_management(l_250, merge_250)     merge_500 = os.path.join(temp_workspace, "500_merge.shp")     arcpy.Merge_management(l_500, merge_500)     merge_750 = os.path.join(temp_workspace, "750_merge.shp")     arcpy.Merge_management(l_750, merge_750)          dissolve_250 = os.path.join(output_workspace, "250_dissolve.shp") # create 250_dissolve.shp for final output     arcpy.Dissolve_management(merge_250, dissolve_250)     dissolve_500 = os.path.join(output_workspace, "500_dissolve.shp") # create 500_dissolve.shp for final output     arcpy.Dissolve_management(merge_500, dissolve_500)     dissolve_750 = os.path.join(output_workspace, "750_dissolve.shp") # create 750_dissolve.shp for final output     arcpy.Dissolve_management(merge_750, dissolve_750)


title: Creating Legends for QGIS 2.0 Feature Blend Mode Outputs

text: I am working on visualizing ~4100 buffered point features. The shading of the buffers is keyed to a field in the data. Because the buffers overlap, I turned on multiply in the new feature blend mode options. I like the visual result, but what seems to be missing is an ability to generate a meaningful legend for the resulting image. Are there any workarounds for this? Is the only alternative to venture down the path of polygon intersects and symmetrical difference operations? ![buffered points with multiply feature blend mode](http://i.stack.imgur.com/ajhrQ.pngb)


title: How to use the raster calculator in the processing modeler of QGIS

text: I am trying to create a easy way to do raster analysis by building a model for two way raster calculations. I have two raster layers that I want to combine into a single raster. For now I've done this manually by doing the raster calculation myself but for the future I want to be able to input these layers into a model without having to go through the trouble. What is the proper annotation of this formula? How I've tried it is not working. ![enter image description here](http://i.stack.imgur.com/eZDZE.png)


title: Arcgis Server failed to create site?

text: Failed to create the site. Admin URL 'http://my.site.url:6080/arcgis/admin' is not reachable from 'LOCALHOST': 'Connection timed out'. Anyone have any idea why i am facing this error? i have installed arcgis 10 server of Red Hat and now i am trying to create site it is giving me this ? i have tried restarting serve but problem won't go away?


title: Cleaning up data in Excel for use in QGIS

text: I am trying to clean up a large .txt dataset in Excel. I think that a macro or another application may be able to clean up the data they way I want but I have limited programming skills. I need to create fields from different sections in the .txt file so I can then geocode the points. I am conducting a crime analysis for my local police department and up until now I have been using wildcards and the Find & Replace tool to do this, but I have nearly 2 years worth of police calls. Fields needed:"Address", "City", "State", "Type", "Location", "Date", "Day of Week", "Time Received", "Time Dispatched", "Time Arrived", "Time Cleared". I put up the data from January 2012 on Github in repo called SLO-Crime under my name, chadbunn. I will be using CartoDB to geocode the data and then QGIS to perform the analysis. I will probably make some maps in TileMill as well if I have enough time to. Any input would be appreciated because it is getting tedious! Thanks! https://github.com/chadbunn/SLO-Crime/tree/master Text sample below showing a single incident (plain text file, each section composed of a dispatch header summary followed by type / address / officer and other information):               ===============================================================================                                  120101001 01/01/12 Received:00:05 Dispatched:00:06 Arrived:00:12 Cleared:00:53                                   ===============================================================================                                  Type: MENTAL SUBJ                                             Location:LZ1                                   As Observed:                                           Mental Health                                       Addr: 1257 RICH, San Luis Obispo, CA             Clearance Code:Report Filed in                                       Responsible Officer: Hyman, J                                    Units: 4235  ,4241                                    Des: incid#=120101001 AP/KAVANAGH,KOLLEEN 060455 166 PC clr:RTF oc:MENH call=1l


title: What is difference between GeometryService areasAndLengths method and Measurement widget in ArcGIS JavaScript API?

text: What's the difference between the GeometryService areasAndLengths method and the Measurement widget to make a caculation in the ArcGIS JS API? Accuracy?


title: Calculate population coverage

text: I am relatively new to QGIS and hope I explain this right. I am looking at the percentage of an population that has geographical access to healthcare. Therefore I have one shapefile with village boundaries including population number of each village. Another shapefile is all Healthcare facilities in the area buffered with 5km zones. I now want to calculate the percentage of the village area that is covered by the Healthcare facility buffers and multiply this by the my population numbers. I did manage to create an intersect shapefile of the two formal shapefiles that created correct new polygons of intersection per village. I now however do not know how to sum the areas (Area_1) of all these buffer polygons who are coded equally in the attributable table (under P_06_ID) and divide this by the total village area (AREA). And than multiply this number by the population (Population). All help welcome! ![enter image description here](http://i.stack.imgur.com/W1cbF.png) ![enter image description here](http://i.stack.imgur.com/EQVqh.png) * * * Thank you so much for the help. The command does indeed calculate the percentage of the population that overlaps with that specific polygon. However it still does not sum the different polygons. ![output after command](http://i.stack.imgur.com/W1cbF.png) So what I now want to do is to sum those numbers with the same parish ID automatically (I do have around 500 parishes). So in this example sum PopPerc from cell 699- 709. Any suggestions?


title: How to use pgRouting's dijkstra_sp_delta_directed with WHERE clause?

text: Right now I am using PgRouting for vehicle routing and I have problem with Dijkstra shortest path. My query is:               SELECT * FROM shortest_path('SELECT gid  AS id,           source::int4 AS source,           target::int4 AS target,           cost::float8 AS cost,           reverse_cost::float8 AS reverse_cost     from network ',                                257027,     276521,     true,     true)      But, my table network is very large, so I should use "dijkstra_sp_delta_directed" wrapper function with bounding box to make query faster.               Select gid from dijkstra_sp_delta_directed('network', 5700, 6733, 0.1,true,true);      But I don't know how to write some where clause like "from network where road_class=1" Also, I don't know where and how to write (road_class * cost)::float8 AS cost Thank you in advance.


title: Updating the value of a field with pre-defined values

text: I am trying change the value of a row for a field whose type is "short" in arcmap using a python script.               with arcpy.da.UpdateCursor(FILENAME, FIELD) as cursor:              for row in cursor:             row = 'Center'             cursor.updateRow(row)      with this i receive the error "RuntimeError: The value type is incompatible with the field type." I know I am giving a string type to a "short" type field. Is there a way around this? There are pre-set values in this field. Do I need to run some other kind of code to create a list of possible values for this field, then select the value based on some kind of list interaction eg `row = fieldlist[1]` Please ask me to clarify anything, but this is everything that seems relevant to me at this point.


title: Switching from ArcDesktop to ArcGIS Online

text: I have many GIS files that I was using to create some maps for my work. Our desktop version of ArcGIS expired and we decided to move to the basic online version. I am unable to upload any of my databases and all the shape files that I have are said to have over 1,000 features. Is there any way I can get this to upload? Maybe splitting the data in two parts and uploading them separately? Any help would be appreciated as all my data that I've been working on is completely useless right now.


title: How do I generate a random raster (or probabilistic cellular automaton) in QGIS?

text: I'm trying to create something like a cellular automaton land cover change model using QGIS, in which cells randomly change land cover based on probabilities determined from various input rasters. I've already have a map of probabilities of conversion, and I need to take the next step. I found this question, How do I randomly alter a raster map of habitat types?, which looks very similar to what I want to do. The suggested approach is to use a Markov chain or probabilistic cellular automaton, and the commenter writes: "The transitions can be implemented using Map Algebra statements, which makes this method practicable in any raster-based GIS, even those without direct or quick access to cell-level data." Can anybody tell me how to do probabilistic map algebra in QGIS or one of its extensions? I don't see any probability tools in the Raster Calculator. Failing that, is there a tool that will simply generate random rasters to be later combined using the Raster Calculator?


title: gdal_retile - Out of range band requested

text: I'm using gdal_retile to create ImagePyramid from a set of GeoTIFF images. I'm using following command:               gdal_retile.py -v -r bilinear -levels 8 -ps 2048 2048 -co "TILED=YES" -co "COMPRESS=JPEG" -targetDir F:/tiles --optfile F:/tileslist.txt      After couple of hours(source images are 320GB) I get following error message:               Traceback (most recent call last):       File "C:\PROGRA~2\FWTOOL~1.7\bin\gdal_retile.py", line 940, in ?         sys.exit(main(sys.argv))       File "C:\PROGRA~2\FWTOOL~1.7\bin\gdal_retile.py", line 857, in main         dsCreatedTileIndex = tileImage(minfo,ti)       File "C:\PROGRA~2\FWTOOL~1.7\bin\gdal_retile.py", line 335, in tileImage         createTile(minfo, offsetX, offsetY, width, height,tilename,OGRDS)       File "C:\PROGRA~2\FWTOOL~1.7\bin\gdal_retile.py", line 467, in createTile         dec.uly+offsetY*dec.scaleY)       File "C:\PROGRA~2\FWTOOL~1.7\bin\gdal_retile.py", line 246, in getDataSet         s_band = sourceDS.GetRasterBand( bandNr )       File "C:\PROGRA~2\FWTOOL~1.7\pymod\gdal.py", line 654, in GetRasterBand         raise ValueError, 'Out of range band requested: %d' % i     ValueError: Out of range band requested: 2      I tried to retile separately images that execution stops on, but then I get no error message. Everything works well, so my guess that it's something wrong with the source image doesn't seem to be valid. Does anyone have any other ideas of what could be wrong or how I should troubleshoot the problem?


title: Fastest way to remove matched points

text: Hoping someone can see a faster way of doing this, because at the moment I'm averaging 9 hours per run for what I'm trying to do and well basically getting my behind chewed off by the boss :-) **The Scenario** We have a database running under PostgreSQL/postgis that's used to manage some mapping data for management of a transport fleet. Within this database there are several tables but for my question here, we are interested only in "built up areas" (aka limitedzones) and "Movement Speeds" (aka speeds) - I can't go too much into detail because of NDA's etc, but "limitedzones" is a table that contains a number of MULTIPOLYGONS each representing a geographic area, and has an associated speed limit that a vehicle is allowed to travel whilst inside the bounds of that polygon. "speeds" is a table of POINTS each point following the path of an allowed road and sampled with an allowed speed limit at 1 second intervals, representing a road that a vehicle is allowed to travel on. At present some of these "allowed roads" cross through "Limited areas" where both the "allowed road" and "Limited Area" speed limit are the same. I've been tasked with removing these extra points from the "Allowed Roads" table where they match the conditions set out above. **What I've tried so far** My first attempt was to simply select the points, so lets say for 40 Mph and put them into a memory table, then loop over that table comparing them to the 'limitedareas' table using "ST_Intersects" this didn't work as expected, and failed to remove all the points. It did however give me an idea of how many points I was looking at removing, and the count was easily in excess of 2 million. After a bit more fiddling around, and a few more attempts I finally arrived at the following:               DELETE FROM speeds s     USING limitedzones lz     WHERE ST_Within(s.geometry,lz.geometry)     AND s.speed = lz.speed      and in a similar fashion used the following to interrogate & count:               SELECT s.* FROM speeds s     JOIN limitedzones lz ON ST_Within(s.geometry,lz.geometry)     WHERE s.speed = lz.speed     --LIMIT 10      Both of these work exactly the way I expect, but the select (When not limited) takes about 3 hours to run, and the delete over 9 hours!! If at all possible I need to find a faster way of stripping these points out, as there is likely to be further work to add this into a regular update system, used ad-hoc on a frequent basis to keep the database in order. I have good indexes both B-Tree (gist) and regular on all the columns, and when individual tables are used on there own, everything is nice and fast, so I know it's the join on 2 tables that's causing me the slowdown. Any suggestions are welcome, and I'm even open to do this using a custom function / procedure if need be.


title: A quick-fix solution for geocoding?

text: I am working on a prototype for a web-based trip planner and need to have a geocoding function (for the start and destination adresses), up and running the "easy" way. My area that the prototype is covering is apprx.2*4 sqkm. Which option is the fastest and easiest (to learn) as I have very limited time? I have been looking into gisgraphy (downloading and geocoding in the database), nominatim/Mapquest etc. I prefer not be constrained to use a specific map provider, like Google or yahoo. However, if necessary it is possible to change geocoder later on... for now I need the "quick fix" Thanx!


title: How to symbolize more than one attribute simultaneously in QGIS?

text: Each symbol can show different properties such as temperature and pressure, etc. I want to change the symbol color according to one property and size with another one. How can I do it in QGIS? Thanks!


title: How do you scan/quantify vegetation using satellite images?

text: I am relatively new to GIS (about 4-5 months in) I did not receive any formal education about this so I am completely clueless about how to go about measuring/quantifying vegetation using satellite images. Can somebody please help? Is there any specific kind of image file I need to do this, if there is, how/where can I get them? Thanks! Enrique Roa Undergrad biology student