ArcGIS 10.1 - c# custom geoprocessing function tool - how to populate a polyline output parameter?

I'm at 10.1SP1 with VS2010. In my tool's ParameterInfo property I create an output GPParameter object of IGPParameterEdit.DataType GPLineTypeClass. Then in the Execute() method I do work that creates an IPolyline that I want to stuff into that parameter's IGPParameter.Value property as an IGPValue. **PROBLEM 1** : Using the typical               IGPValue outValue = gputils.MakeGPValueFromObject(myIPolyline)       using the GPUtilities class's IGPUtilities interface causes an E_FAIL COM error. huh. How does one debug that? Well, try another way - So I tried this instead:               IGPLine2 outValue = new GPLineClass();     outValue.Polyline = myIPolyline;     gpUtils.PackGPValue((IGPValue)outValue, outputParam);      **PROBLEM 2** : But I get an invalid cast exception on that first line: `System.InvalidCastException: Unable to cast COM object of type 'ESRI.ArcGIS.Geoprocessing.GPLineClass' to interface type 'ESRI.ArcGIS.Geoprocessing.IGPLine2'. This operation failed because the QueryInterface call on the COM component for the interface with IID '{1248339D-7B76-4A43-8FE0-752B1BF9D39A}' failed due to the following error: No such interface supported (Exception from HRESULT: 0x80004002 (E_NOINTERFACE))` GPLineClass cannot be cast to IGPLine2. What! Documentation says this class implements this interface. I CAN cast it to IGPLine, but that interface is ridiculous - I'd have to disassemble the IPolyline and assign each vertex, and it won't work for multipart geometries. SO - what is the secret sauce to stuff an IPolyline into an output parameter's value? And what the heck is up with that IGPLine2 thing? Thanks!


What is the unit of measurement for buffer calculation?

What is the unit of measurement used by Geometry class for it's buffer calculation ? Api doc currently defines it as " distance - the width of the buffer (may be positive, negative or 0)" I need to draw polygons around a given Geometry with X miles away. Wondering if there is an easy way to do it with current set of GeoTools API. Api : http://www.vividsolutions.com/jts/javadoc/com/vividsolutions/jts/geom/Geometry.html#buffer(double)


How to decode Magellan URL coordinates to lat/lon

Magellan uses this URL to designate a coordinate in New York:               mgn:///address/?coord=merwFvlqbM,name=Viva%20Herbals%20Pizzeria,number=179,street=2nd%20Avenue,city=New%20York,country=US      How can I decode the "coord" part (`merwFvlqbM`) of the URL? It looks like an encoded latitude longitude value. Other NY coords have similar values: glp **wF** |fs **bM** ubr **wF** r~p **bM** (common characters bold)


OGR GetX() returns always zero

I am trying to get the first x-coordinates of a polygon. However I only get zero in return. What could cause that?               import os, sys, gdal, ogr     from gdalconst import *          driver = ogr.GetDriverByName('ESRI Shapefile')     testarea = driver.Open('testarea.shp', 0)     testarealyr = testarea.GetLayer()     testareafeature = testarealyr.GetNextFeature()          # get the x,y coordinates for the point     testareageom = testareafeature.GetGeometryRef()     print testareageom          x = testareageom.GetX()     print x               >>>      POLYGON ((-124.1963006330602 43.006410659375554,-124.1861022086067 43.006647759060762,-124.1858958821004 43.00274627515271,-124.19612378176909 43.002422936639086,-124.19612378176909 43.002422936639086,-124.1963006330602 43.006410659375554))     0.0     >>>


Looping based on a spreasheet using python (in ArcGIS 9.3)

Here's what I am trying to do: I want to move feature classes from one `.gdb` (OLD) to another (NEW), but only the feature classes that are present in an excel spreadsheet.   1. read from excel spreasheet (old name column)   2. compare to old .gdb feature classes to find match   3. move matched feature class to new .gdb   4. rename feature class in new .gdb based on new name in excel spreadsheet (new name column)   5. continue doing the above until the end of the spread sheet


How to populate polygon extents to an attribute table in QGIS

Is there a way to calculate the XMIN, YMIN, XMAX and YMAX of a polygon and write it to a field in an attribute table? The Geometry functions available in the Field Calculator are for points. Thanks...


How to insert google maps lat/lng into postgis in geometry datatype

I came across: https://developers.google.com/maps/articles/phpsqlinfo_v3 I've changed the code a bit and all process is well done for postgis database. I get the data from infobox when user clicked. ..                var map = new google.maps.Map(document.getElementById("map_canvas"), options);           var html = "
" +                      "
Name:	 <input type='text' id='name'/>
" +                      "
Address:	 <input type='text' id='address'/>
" +                      "
Type:	 <select id='type'>" +                      "<option value='bar' SELECTED>bar</option>" +                      "<option value='restaurant'>restaurant</option>" +                      "</select>
" +                      "
	<input type='button' value='Save & Close' onclick='saveData()'/>
";         infowindow = new google.maps.InfoWindow({          content: html         });      I get the data from another php via url;               // Insert new row with user data     $query = sprintf("INSERT INTO markers " .              " (id, name, address, lat, lng, type ) " .              " VALUES (DEFAULT, '%s', '%s', '%s', '%s', '%s');",              pg_escape_string($name),              pg_escape_string($address),              pg_escape_string($lat),              pg_escape_string($lng),              pg_escape_string($type));          $result = pg_query($query);      But the problem is i want to store google maps lat/lng not as float, but as geometry(propably point) data type in a column. How should i change the code or philosophy ?


Build a model sequence or merge models in ArcGIS

This is surprisingly difficult to search for, so I'll ask though it may be a silly question. Is it possible to chain several model-builder models in sequence, such that one will automatically run with the precondition that a prior model has executed? There are plenty, so it seems more sensible than a merge, but I would be interested in that as well.


How to read elevation data from WMS response?

I am trying to get elevation data from the tile received from USGS server. I query the data using Image Format image/Png. When I read the data I get color values. I need to know how to read the elevation values from it. Is there any documentation or sample available.


Can I modify a layer which is uploaded by myself?

I first uploaded a geojson type polygon to the OpenLayer. Can I still edit the feature of this layer? Here is the code I am using: (It can upload the geojson layer but can not edit it)                       var map, layer;                  function init(){             map = new OpenLayers.Map( 'map');             layer = new OpenLayers.Layer.OSM( "Simple OSM Map");             map.addLayer(layer);             map.setCenter(                  new OpenLayers.LonLat(0.8,52.472).transform(                     new OpenLayers.Projection("EPSG:4326"),                     map.getProjectionObject()                 ), 7             );                          var featurecollection = {"type":"Polygon","coordinates":[[[-248052.78111938,6928922.6476283],[-193018.12076172,6933661.7433813],[-167335.27926147,6885812.1636815],[-184915.79576462,6855084.4783151],[-228484.9018811,6854931.6042586],[-260588.45375641,6872512.1207617],[-266703.41601837,6898194.962262],[-248052.78111938,6928922.6476283]]]};            var geojson_format = new OpenLayers.Format.GeoJSON();            var vector_layer = new OpenLayers.Layer.Vector({             projection: new OpenLayers.Projection("EPSG:4326")            });                  map.addLayer(vector_layer);            vector_layer.addFeatures(geojson_format.read(featurecollection));          //Until Here, No Problem                      var select_control = new OpenLayers.Control.SelectFeature(vector_layer);             select_control.activate();                  var modify_control = new OpenLayers.Control.ModifyFeature(vector_layer,{                 mode:OpenLayers.Control.ModifyFeature.RESHAPE             });             map.addControl(modify_control);             modify_control.activate();              }      Is there something wrong with the editing code? Thanks in advance.


Access multiple feature servers through ArcGIS REST API

we are adapting a web app to support an ArcGIS server back end through the use of Javascript ajax calls which in turn display the data using Leaflet. The current version of the application runs queries against a PostGIS database which allows us to join tables and perform geospatial operations in a single server call. For example, given a lat/lng point we will need to find an intersecting polygon in one layer, then find polygons in other layers that are within a set distance of that polygon. In PostGIS we are able to accomplish this with one query, but it appears we will have to send out multiple requests when using ArcGIS server. Is there any way to achieve a similar level of expressiveness with the ArcGIS REST API, or will we have to send out a series of requests to get the same results?


ArcGIS geoprocessor memory leak: what's the underlying cause and why is it so hard to fix?

blah238's very nice answer to this question reminded me of all the time I've spent over the past two years battling against the dreaded geoprocessor memory leak in ArcGIS Desktop. Ah, happy days! I only have personal experience of problems with version 9.3.1, but there are numerous threads on the ESRI support pages (e.g. here, here) which suggest that the problem goes back to at least v9.1. blah238's answer and this thread indicate that the problem still isn't fixed at v10. I know that workarounds are possible, usually involving the use of subprocesses to force the gp to release the memory it's hogging, but they're clunky and often very slow. This question isn't meant as an attack on ESRI, but I'd be really interested to know: **Why has this problem proved so difficult to fix?** For my work at least, it's a major issue and one that's pushed me away from using ArcGIS for spatial modelling. This has been a great learning experience, but I'm surprised that fixing the problem isn't higher on ESRI's priority list. Any thoughts or comments welcome. Cheers!


Masking WMS (Esri Oceans Basemap) layer in ArcGIS for Desktop?

I would like to display ESRI's Oceans basemap web map service (WMS) in one of my maps, but not the whole layer, just the landmass depicting Canada as I find it a nice fit for the look I'm trying to achieve. Is there any way to mask the features around Canada's landmass in that layer? I tried using the Advanced Drawing feature, but the Oceans basemap layer is not listed, probably because it's a WMS. I am working with ArcGIS 10.2. Are there any suggestions as to how I could do this? Otherwise, if anyone can suggest a symbolized layer of Canada's landmass (and only containing Canada's landmass) ressembling that depicted on the ESRI Oceans Basemap that I could download or connect to somewhere, I would be grateful. I could make my own. I tried with a few DEM's and hillshades but didn't get good results.


Tools to make buffer around the polygon on a web map?

We're starting web map app and still deciding what tools is the best to use. We would need option that user add some points and polygons, and also, to automatically get buffer around those polygons. We considered Geoserver as the best solution, also Postgresql as the database. Can you advice us regards to 'buffer' functions, which tools can we use for buffering? Should we use OpenLayers or Leaflet is 'good enough'? Thanks in advance.


How to add layer/raster with specific position in the QGIS layer list?

I use addRasterLayer(rast,'Name') but this adds the raster to the top. It is possible to **set position** when I add raster at once? Is it possible??


Looking for Python module to do lat lon calculations

I'm looking for a Python module that can provide some / all of the following services using latitude/longitude to define points:               Distance between two points     Point in polygon      Distance along a multi-point path      Identifier of closest point from a dict of points with identifiers      Distance from a point to each point in a dict containing points with identifiers     etc      You get the idea... Any suggestions?


Error reported when trying to import prj for data points in New Brunswick Canada into QGIS

I have data points recorded with a Trimble GIS. The projection is New Brunswick Sterographic Double Projection, which as near as I can find via a Google search is EPSG:2953. However, when I try to import the proj (see below) I get the error "not a valid or recognised data source". The same files import just fine into ARCMAP 10. How do I get these files to read correctly into QGIS 1.7.3 (Windows 7) or does this mean there is a bug in the projection data base of QGIS? I guess I could drop the projection and specify it in the program but I am concerned that the projection may be wrong if the files don't import correctly in the first place. Thanks Alyre Projection file:               PROJCS["NAD_1983_CSRS_New_Brunswick_Stereographic",GEOGCS["GCS_North_American_1983_CSRS",DATUM["D_North_American_1983_CSRS",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Double_Stereographic"],PARAMETER["False_Easting",2500000.0],PARAMETER["False_Northing",7500000.0],PARAMETER["Central_Meridian",-66.5],PARAMETER["Scale_Factor",0.999912],PARAMETER["Latitude_Of_Origin",46.5],UNIT["Meter",1.0]]


Problem Installing pgRouting 2.0; libcgal8 dependency

I've been trying to install PostGIS and pgRouting on an AWS Ubuntu server (first tried 13.10, then reverted to 12.04, both 64bit "small" instance) for several days now but no matter what I do I keep getting hung up by the dependency libcgal8. Among many other things, I've tried using exactly the steps @dkastl recommended in this post:               sudo add-apt-repository ppa:ubuntugis/ppa     sudo add-apt-repository ppa:georepublic/pgrouting     sudo apt-get update     sudo apt-get install postgresql-9.1-pgrouting postgresql-9.1-postgis-2.0      After executing the last command the system complained that               Package 'postgresql-9.1-postgis-2.0 is unavailable, but  is referred to by another package      Checking the ubuntugis ppa: I determined that the correct name is `postgresql-9.1-postgis`. Substituting that package name got me to the roadblock I've encountered all along:               The following packages have unmet dependencies:     postgresql-9.1-pgrouting : Depends : libcgal8 but is not installable     E: Unable to correct problems, you have held broken packages.      I've tried everything to get around this issue including the following:   * Separating the postgis and pgrouting install commands and switching the order they were executed   * Tried using `$ sudo apt-get install postgresql-9.1-pgrouting --no-install-recommends` but this didn't work either (still complained about unmet libcgal8 dependency)   * Tried downloading libcgal8 and its dependencies and installing them directly, which was a lost cause   * Using `$ sudo apt-get -f install` which was suggested by the system at some point, and seemed to do something but it didn't fix the problem   * So many other things that I either can't remember or were too desperate to mention Hoping one of you can point to something really silly that this newb is doing or not doing. Thanks in advance for your help.


QGIS: field calculator is not working with shapefile created manually

I have an issue the field calculator of GIS 2.2 field calculator. I've not encounter problem to update values with the field calculator (I'm using the $row_num function) if the shapefile was created with another GIS software or if for example I used the "Regular Points" tool to create a new shapefile. However, if I've created the shapefile "manually" ("New shapefile layer") with QGIS, I can't update field values with $row_num with the field calculator (but it's working fine with the exact same shapefile created manually with QGIS in ArcGIS 10). Am I missing something or it's a bug? My colleague has the same problem on his computer and with completly different shapefiles.


ArcGIS Python - Convert feature class names to lower case

Can anybody anyone point me to directions on how to convert file geodatabase feature class names to lower case by using python? Thanks


Oracle Spatial support QGIS 1.9-0-Master (Ubuntu)

I just installed QGIS 1.9.0-Master (Build 480) on a Ubuntu virtual machine. I am running/compiled against GDAL/OGR 1.9.2. My question is why I cannot add Oracle Spatial Layers as I am able to do using QGIS-Master for Windows. Is this under development for Linux systems or may I missed some package during installation? Best Regards Daniel


QGIS 1.8.0 on Linux treats NetCDF 4 as HDF5

I've been generating some large (28800x13920) grids of BYTE data representing the NALCMS2005 North American Land Cover dataset (original 250m projected data resampled to a WGS84 grid at 15 arc sec). I store the grid as a NetCDF 4 file with GZIP compression, using gzip level 6 and CF conventions; i.e. I use a "crs" variable to define the coordinate reference system, "lat" and "lon" dimensions and coordinate variables to define the geodetic latitude and longitude axes, and "lat_bounds" and "lon_bounds" variables to define the bounds of each row/column of the grid. When I select one of these files with qgis-1.8.0-9_0.el6.x86_64 on a CentOS system or with qgis-1.8.0-6.fc16.i686 on a Fedora Core 16 system (these are my two Linux options right now), the "Select raster layers to add..." dialog that lists the layers found in the file displays them as HDF5:"NA_LandCover_2005_WGS84_500m.nc"://Data QGIS is able to load the data and display it, but the coordinate reference system and grid geometry are completely ignored. I realize that NetCDF4 files are technically also HDF5 files, but read on. When I select one of the same files with QGIS on a Mac (Mountain Lion 10.8.4), the "Select raster layers to add..." dialog show the layers in the file as NETCDF:"NA_LandCover_2005_WGS84_500m.nc"://Data and successfully recognizes the CRS in the file. Unfortunately, this solution is not doing a lot for me because reading and displaying the data on this system takes about 30 minutes, and good luck changing anything or doing any navigation. So... if I do "gdalinfo --formats | grep -i df" on the different systems I get exactly the same result: GMT (rw): GMT NetCDF Grid Format netCDF (rw+): Network Common Data Format HDF4 (ro): Hierarchical Data Format Release 4 HDF4Image (rw+): HDF4 Dataset NDF (rov): NLAPS Data Format HDF5 (ro): Hierarchical Data Format Release 5 HDF5Image (ro): HDF5 Dataset "gdalinfo --version" on the three systems:   * Mac: GDAL 1.9.1, released 2012/05/15   * FC16: GDAL 1.7.3, released 2010/11/10   * CentOS 6.3: GDAL 1.8.1, released 2011/07/09 Maybe those older GDAL libraries on the Linux systems are the problem and I just need to build QGIS from source using GDAL 1.9 on one of the Linux machines.


ArcGIS 401 error during initialization of javascript map

I have a dumbed-down application to only reference a geometry service, and a feature service. When the map is being initialized, it sends out a GET request of: `http://my.server.com/ArcGISTEST/rest/info?f=json` but when using Chrome debug tools / fiddler I see it is returning HTTP 401 error when I open this url in its own window/tab it produces the following output: `{"currentVersion":10.02,"soapUrl":"a url is in here","secureSoapUrl":"a url is in here","authInfo":{"isTokenBasedSecurity":false}}` It appears this is an object to determine the status/capability of the arcgis server to assist the javascript api. Why am I getting this error? How do I resolve it? Could this be a symptom of my bigger problem (unable to edit features via javascript api).


Need help using FieldMappings object with Append_management in Python

I have a fairly large script that I'm reworking. The script is manipulating tables with 50 or so fields. (I will be cleaning some of these up to reduce the number.) Toward the end of my script I call the Append_management tool and and trying to simplify the field mapping component of that process by utilizing the FieldMappings object. **Question:** Is it possible to use the FieldMappings object only for the fields where the input and target field names differ? In other words; I would like to tell the script that if it finds matching fields, to use those, and for all other that where the field names differ, refer to the FieldMappings object. Does this make sense? Below is a sample of my Append_management section of code.               arcpy.Append_management(PZSPJoin, out_fc, "NO_TEST", "Join_Count \"Join_Count\" true true false 4 Long 0 10 ,First,#,Jun12May13PZoneSPJoin,Join_Count,-1,-1; \                                                                             TARGET_FID \"TARGET_FID\" true true false 4 Long 0 10 ,First,#,Jun12May13PZoneSPJoin,TARGET_FID,-1,-1; \                                                                             Status \"Status\" true true false 1 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,Status,-1,-1; \                                                                             Score \"Score\" true true false 2 Short 0 5 ,First,#,Jun12May13PZoneSPJoin,Score,-1,-1; \                                                                             Match_type \"Match_type\" true true false 2 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,Match_type,-1,-1; \                                                                             Stan_addr \"Stan_addr\" true true false 79 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,Stan_addr,-1,-1; \                                                                             Match_addr \"Match_addr\" true true false 103 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,Match_addr,-1,-1; \                                                                             ARC_Street \"Street\" true true false 60 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,ARC_Street,-1,-1; \                                                                             ARC_Zone \"Zone\" true true false 40 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,ARC_Zone,-1,-1; \                                                                             ADDRESS \"ADDRESS\" true true false 37 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,ADDRESS,-1,-1; \                                                                             CITY \"CITY\" true true false 20 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,CITY,-1,-1; \                                                                             ZIP \"ZIP\" true true false 5 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,ZIP,-1,-1; \                                                                             STSUB \"STSUB\" true true false 16 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,STSUB,-1,-1; \                                                                             POSTDIR \"POSTDIR\" true true false 3 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,POSTDIR,-1,-1; \                                                                             UNITID \"UNITID\" true true false 16 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,UNITID,-1,-1; \                                                                             XCOORD \"XCOORD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,XCOORD,-1,-1; \                                                                             YCOORD \"YCOORD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,YCOORD,-1,-1; \                                                                             JUN_12CCF \"JUN_12CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,JUN_12,-1,-1; \                                                                             JUL_12CCF \"JUL_12CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,JUL_12,-1,-1; \                                                                             AUG_12CCF \"AUG_12CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,AUG_12,-1,-1; \                                                                             SEP_12CCF \"SEP_12CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,SEP_12,-1,-1; \                                                                             OCT_12CCF \"OCT_12CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,OCT_12,-1,-1; \                                                                             NOV_12CCF \"NOV_12CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,NOV_12,-1,-1; \                                                                             DEC_12CCF \"DEC_12CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,DEC_12,-1,-1; \                                                                             JAN_13CCF \"JAN_13CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,JAN_13,-1,-1; \                                                                             FEB_13CCF \"FEB_13CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,FEB_13,-1,-1; \                                                                             MAR_13CCF \"MAR_13CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,MAR_13,-1,-1; \                                                                             APR_13CCF \"APR_13CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,APR_13,-1,-1; \                                                                             MAY_13CCF \"MAY_13CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,MAY_13,-1,-1; \                                                                             JUN_12MGD \"JUN_12MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,JUN_12MGD,-1,-1; \                                                                             JUL_12MGD \"JUL_12MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,JUL_12MGD,-1,-1; \                                                                             AUG_12MGD \"AUG_12MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,AUG_12MGD,-1,-1; \                                                                             SEP_12MGD \"SEP_12MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,SEP_12MGD,-1,-1; \                                                                             OCT_12MGD \"OCT_12MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,OCT_12MGD,-1,-1; \                                                                             NOV_12MGD \"NOV_12MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,NOV_12MGD,-1,-1; \                                                                             DEC_12MGD \"DEC_12MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,DEC_12MGD,-1,-1; \                                                                             JAN_13MGD \"JAN_13MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,JAN_13MGD,-1,-1; \                                                                             FEB_13MGD \"FEB_13MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,FEB_13MGD,-1,-1; \                                                                             MAR_13MGD \"MAR_13MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,MAR_13MGD,-1,-1; \                                                                             APR_13MGD \"APR_13MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,APR_13MGD,-1,-1; \                                                                             MAY_13MGD \"MAY_13MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,MAY_13MGD,-1,-1; \                                                                             ACCTCLASS \"ACCTCLASS\" true true false 30 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,ACCTCLASS,-1,-1; \                                                                             LOC_DESC \"LOC_DESC\" true true false 254 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,LOC_DESC,-1,-1; \                                                                             CUSTOMER \"CUSTOMER\" true true false 30 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,CUSTOMER,-1,-1; \                                                                             SUBGROUP \"SUBGROUP\" true true false 10 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,SUBGROUP,-1,-1; \                                                                             PRESSZONE \"PRESSZONE\" true true false 4 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,PRESSZONE,-1,-1; \                                                                             PZONE_A_ID \"PZONE_A_ID\" true true false 4 Long 0 10 ,First,#,Jun12May13PZoneSPJoin,PZONE_A_ID,-1,-1; \                                                                             AREA_ \"AREA_\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,AREA_,-1,-1; \                                                                             PRES_ZONE \"PRES_ZONE\" true true false 4 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,PRES_ZONE,-1,-1; \                                                                             PZ_DESC \"PZ_DESC\" true true false 50 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,PZ_DESC,-1,-1; \                                                                             HGL_FT \"HGL_FT\" true true false 4 Long 0 10 ,First,#,Jun12May13PZoneSPJoin,HGL_FT,-1,-1; \                                                                             WELL_CAP_GPM \"WELL_CAP_GPM\" true true false 4 Long 0 10 ,First,#,Jun12May13PZoneSPJoin,WELL_CAP_GPM,-1,-1; \                                                                             B_PUMP_GPM \"B_PUMP_GPM\" true true false 4 Long 0 10 ,First,#,Jun12May13PZoneSPJoin,B_PUMP_GPM,-1,-1; \                                                                             GND_STOR_MG \"GND_STOR_MG\" true true false 4 Long 0 10 ,First,#,Jun12May13PZoneSPJoin,GND_STOR_MG,-1,-1; \                                                                             ELEV_STOR_MG \"ELEV_STOR_MG\" true true false 4 Long 0 10 ,First,#,Jun12May13PZoneSPJoin,ELEV_STOR_MG,-1,-1; \                                                                             CREATEUSER \"CREATEUSER\" true true false 30 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,CREATEUSER,-1,-1; \                                                                             CREATEDATE \"CREATEDATE\" true true false 36 Date 0 0 ,First,#,Jun12May13PZoneSPJoin,CREATEDATE,-1,-1; \                                                                             MODUSER \"MODUSER\" true true false 30 Text 0 0 ,First,#,Jun12May13PZoneSPJoin,MODUSER,-1,-1; \                                                                             MODDATE \"MODDATE\" true true false 36 Date 0 0 ,First,#,Jun12May13PZoneSPJoin,MODDATE,-1,-1; \                                                                             TotalCCF \"TotalCCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,TotalCCF,-1,-1; \                                                                             TotalMGD \"TotalMGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,TotalMGD,-1,-1; \                                                                             AnnualAVG_CCF \"AnnualAVG_CCF\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,AnnualAVG_CCF,-1,-1; \                                                                             AnnualAVG_MGD \"AnnualAVG_MGD\" true true false 8 Double 8 38 ,First,#,Jun12May13PZoneSPJoin,AnnualAVG_MGD,-1,-1; \                                                                             Download \"Download\" true true false 36 Date 0 0 ,First,#,Jun12May13PZoneSPJoin,Download,-1,-1", "")


Wrong size when use Mapproxy to serve Mapnik tiles as ESPG: 900913

I have been stuck at this problem for a whole day. Really would like to ask for some help. Problem: In MapProxy Demo TMS page, EPSG:4326 looks fine, except the map is tilted and I don't want that. Tried to make EPSG:900913 work, However, no matter what I try. the map is too small, at zoom level 18 I can not see street, the whole city remains a little dot on screen for all zoom levels. Or sometime I got nothing. Background:   * Street data stored in Spatialite in Lat Lng format, and all geometry has srid 4269   * Uses MapProxy 1.5 and Mapnik 2.2.0   * tried to append 900913 to mapnik\share\proj\epsg file (seems mapnik don't support 900913 natively). randomly give value such as      * <900913> +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs <>     * <900913> +proj=longlat +datum=WGS84 +no_defs <>   * randomly change all srs string in mapnik xml config to srs strings above   * Map looks fine in TileMap with SRS "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", map looks small and tilted with SRS "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over"   * I thought Lat Lng data is absolute format, any projection should works fine with my data.   * Clear Cache and restart server every time.


How to serve OSM Tiles on Windows?

Was hoping I could just extract tiles into a local file system, point IIS at it and be done but I guess that is not the reality of Open Street Map. I'm getting buried in all the possible ways of self-hosting OSM tiles but unable to find simple step-by-step instructions. I've identified PostgreSQL, PostGIS, Mapnik as key players but can't get all the pieces together. Any step-by-step instructions on creating a custom OSM tile server?


Why the GIS industry is always ruled by two options only, namely ESRI products and OSGeo projects?

I always wonder, why in GIS industry there are merely two major options available. Although there are many other proprietary software products for GIS, their market share is well below ESRI products (here, here and here although quite back dated). My opinion became stronger after joining GIS SE. Most of the **questions** , as well as the leading **users** (with highest reputation), are either for/interested in **ESRI line of products** or **FOSS** (of course apart from some questions related to data and theory). However this not the case in other streams of information technology, databases and programming languages in broad sense. I personally feel one reason could be the lack of interest by the software giants like Microsoft, Google or Oracle. While the first two have WebMaps, these are not WebGIS, and the spatial extension of Oracle is not fully GIS too. What are the major reasons for such a trend? There is nice article somewhat related to this, " The Five Major Roadblocks to GIS Gaining Corporate Market Share". But this is not the answer I am looking for. **EDIT:** I am afraid. Most of the answers seem to prove open source tools are better than ESRI products! But that was not the question. Because it is difficult to compare a matured but expensive tool with a bunch of comparatively newer open source tools with great usability. My question is why no industry standard product from any established software giants? Hope this helps...


Mobile location determination

I was wondering how Google Maps on my Android phone can locate me with a good accuracy when my GPS is disabled. At first I though it was with the cell tower but I'm not sure it can achieve this precision and it is way better when the wireless is enabled.   1. If Wifi is ON and GPS is ON = Perfect location   2. If Wifi is ON GPS OFF = Still Perfect location   3. If Wifi is OFF and GPS is ON = location is not presisce. Tested this with 5 -6 Phones Any Idea why ?


How are housenumbers and streets linked in a OSM based PostGIS DB?

I'm trying to understand how OSM is linking housenumbers and streets with each other. I downloaded an osm file for germany and imported it into a PostGIS DB. To understand how the linking is working I just took a street with a housenumber and exported it as xml from osm.org. I opened the xml file and saw the entry for the housenumber and the street and it looked like this:               <way id="141827749" user="MichaH" uid="8464" visible="true" version="1" changeset="10149457" timestamp="2011-12-18T18:41:33Z">      <nd ref="1005720791"/>      <nd ref="1552517039"/>      <nd ref="1552517025"/>      <nd ref="1005720152"/>      <nd ref="1005720791"/>      <tag k="addr:city" v="Frankfurt"/>      <tag k="addr:housenumber" v="16"/>      <tag k="addr:postcode" v="60326"/>      <tag k="addr:street" v="Wörsdorfer Straße"/>      <tag k="building" v="yes"/>     </way>      Obviously I took the id=141827749 and found one entry in planet_osm_ways. Result:               id        | nodes                                                    | tags     141827749 | {1005720791,1552517039,1552517025,1005720152,1005720791} | {building,yes,addr:housenumber,16}      So I got the housenumber. But where is the street name? When I search the tags for the street name i get 3 entries pointing to other nodes. The correct entry for it is as I know (because its a residential area):               id      | nodes                                                                                | tags     6235469 | {52082395,1552517028,1552516972,1552516886,1552516856,1552516836,1552516769,1552516705,1552516652,1552516571,52082405} | {oneway,yes,name,"Wörsdorfer Straße",highway,residential}      But where is the link between all this? Does anybody know? I didnt find any info on the OSM wiki pages. Thanks in advance.


How do you import a .SHP into Sketch-Up?

Does anyone know how to use a .SHP in Trimble Sketch-Up? I have tried exporting to KML and CAD from QGIS, but the KML file doesn't read in Sketch-Up. Thanks! -m


Can I have a layer in my map that gets its imagery from say google? (ArcGIS Mobile 10)

I've been messing with ArcGIS Mobile and I get that it loads and caches layers from an ArcGIS Server for disconnected use. I'm wondering if it supports layers where data is loaded from a layer in realtime? I can't find any examples where I create a layer, give it a URL and have it load imagery data as I pan the map.


How can I import layers into a database?

According to the QGIS User Guide (release 1.8.0, page 203/204) I only have to drag a layer from the QGIS Browser to the DB Manager. But for me this does not work.


Intersect multilinestring (A) with multilinestring (B) retaining linear sequencing in resultset?

I have a PostGIS table **testtin2** which is a sort of pseudo TIN consisting of a single multilinestring geometry, which I am intersecting with another table **align** , also a multilinestring, to generate a set of coordinates along an alignment (or alignments) via the following query...               SELECT ST_X(q1.int_pnt), ST_Y(q1.int_pnt), ST_Z(q1.int_pnt)     FROM     (SELECT ST_Intersection(align.geom, testtin2.geom) as int_pnt      FROM align INNER JOIN testtin2 ON ST_Intersects(align.geom, testtin2.geom)     WHERE ST_isvalid(align.geom)='t' AND ST_isvalid(testtin2.geom)='t') As q1;      Which works fine except for the fact that the resulting output is seemingly in a 'random' order. Is there a way this query might be modified or re-written such that the coordinates come out in the correct linear sequence along the alignment linestring(s) from the **align** table? I need this to happen so that the drafting software I am passing the coordinates to will join the dots up in the correct order. Thanks in advance.


How to transform Lambert coordinates in C#?

I want to transform Lambert (EPSG:3035) coordinates to EPSG:4326 in C#. I have tried this:               var lambert = cf.CreateFromWkt("PROJCS[\"ETRS89 / ETRS-LAEA\",GEOGCS[\"ETRS89     \",DATUM[\"European_Terrestrial_Reference_System_1989\",SPHEROID[\"GRS      1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],AUTHORITY[\"EPSG\",\"6258     \"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree     \",0.01745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4258     \"]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"latitude_of_center     \",52],PARAMETER[\"longitude_of_center\",10],PARAMETER[\"false_easting     \",4321000],PARAMETER[\"false_northing\",3210000],UNIT[\"metre\",1,AUTHORITY[\"EPSG     \",\"9001\"]],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH],AUTHORITY[\"EPSG\",\"3035\"]]");      Error: Projection Lambert_Azimuthal_Equal_Area is not supported. As far as I understand, ProjNet should support Lambert. Can someone help me? Or show me another way how to achieve my goal?


What event is fired when a feature is edited?

When the user draws a new feature (line, point, etc.) in ArcMap, there are null attributes. If the user enters values in these null fields what is the event that is fired? _I want to control the entered values._


How to get WMO SYNOP data?

I would like to get WMO weather stations observations daily reports from the web in file or web-service form. According to WMO Resolution 40 -info\- this data should be availiable in text form with SYNOP codes. Does anyone know a place to look for this data? I've looked through NCEP and NASA sites, and didn't find the dataset itself. I've found only http://www.ogimet.com site, that provides data in html and txt form, but you have to specify a WMO station index or a region - so no service for global coverage (all WMO stations). I suppose this site gets information from the resource, that I'm looking for.


How to include mapnik lib in python 2.7

I am using python 2.7 now in want mapnik to be included in python. I have use exe file to include other library file like GDAL, and other. But for Mapnik can any one suggest me any exe like that so it can be easy or library file but how to install it, as i found many help not usefulness. PLEASE any one help me out in this issue


invalid KML representation on ST_GeomFromKML

I exported a table with a `geography(MultiPolygon,4326)` column using **ogr2ogr** to a KML file. When trying to re-import any of the shapes using the `ST_GeomFromKML` function I get the following error:               ERROR:  invalid KML representation          ********** Error **********          ERROR: invalid KML representation     SQL state: XX000      A sample of the KML shape I'm passing as argument is:               <Polygon><outerBoundaryIs><LinearRing><coordinates>-3.106121,53.524247,0.0 -3.102785,53.522356,0.0 -3.099633,53.521059,0.0 -3.097637,53.519834,0.0 -3.096472,53.518869,0.0 -3.095176,53.518035,0.0 -3.092952,53.517169,0.0 -3.090586,53.516356,0.0 -3.089088,53.515628,0.0 -3.08863,53.515296,0.0 -3.088306,53.514773,0.0 -3.088044,53.513979,0.0 -3.088025,53.513251,0.0 -3.087312,53.51078,0.0 -3.085799,53.508947,0.0 -3.084285,53.506689,0.0 -3.083638,53.506559,0.0 -3.083456,53.506456,0.0 -3.082374,53.504838,0.0 -3.081294,53.50436,0.0 -3.080699,53.503849,0.0 -3.079781,53.50176,0.0 -3.080216,53.501094,0.0 -3.080577,53.500779,0.0 -3.080403,53.500497,0.0 -3.080611,53.500357,0.0 -3.080569,53.500262,0.0 -3.079844,53.499901,0.0 -3.078279,53.500012,0.0 -3.077074,53.499674,0.0 -3.077043,53.499568,0.0 -3.077502,53.499601,0.0 -3.077802,53.499724,0.0 -3.078244,53.499719,0.0 -3.078791,53.499417,0.0 -3.07865,53.499375,0.0 -3.078458,53.499047,0.0 -3.078545,53.498628,0.0 -3.078439,53.498512,0.0 -3.078565,53.4983,0.0 -3.07818,53.49812,0.0 -3.078418,53.497947,0.0 -3.078275,53.497887,0.0 -3.078502,53.497779,0.0 -3.078095,53.49778,0.0 -3.077909,53.497649,0.0 -3.076642,53.495559,0.0 -3.074094,53.49256,0.0 -3.073583,53.492224,0.0 -3.071519,53.491591,0.0 -3.069933,53.490344,0.0 -3.067873,53.487934,0.0 -3.066784,53.486822,0.0 -3.066341,53.485681,0.0 -3.064965,53.48452,0.0 -3.064071,53.483344,0.0 -3.063361,53.48261,0.0 -3.062358,53.481779,0.0 -3.060913,53.480844,0.0 -3.059644,53.479436,0.0 -3.058885,53.479084,0.0 -3.058551,53.478809,0.0 -3.057679,53.477685,0.0 -3.056766,53.477329,0.0 -3.056209,53.477321,0.0 -3.054533,53.476541,0.0 -3.054158,53.476071,0.0 -3.053637,53.475686,0.0 -3.050002,53.473287,0.0 -3.049804,53.473026,0.0 -3.049445,53.471779,0.0 -3.049255,53.471421,0.0 -3.048791,53.470924,0.0 -3.04766,53.470316,0.0 -3.046093,53.469765,0.0 -3.042068,53.466788,0.0 -3.041516,53.466488,0.0 -3.041086,53.466535,0.0 -3.041495,53.466257,0.0 -3.041753,53.465722,0.0 -3.041636,53.465445,0.0 -3.041106,53.464933,0.0 -3.028521,53.455075,0.0 -3.028182,53.454884,0.0 -3.027857,53.45489,0.0 -3.026988,53.455444,0.0 -3.025183,53.455136,0.0 -3.023849,53.454794,0.0 -3.034479,53.451557,0.0 -3.028883,53.44732,0.0 -3.028206,53.446716,0.0 -3.024574,53.443565,0.0 -3.02148,53.440243,0.0 -3.019852,53.437712,0.0 -3.019207,53.436152,0.0 -3.008783,53.438358,0.0 -3.008371,53.43831,0.0 -3.007707,53.438365,0.0 -3.00601,53.438769,0.0 -3.005977,53.438903,0.0 -3.00612,53.43908,0.0 -3.005502,53.439219,0.0 -3.005371,53.439045,0.0 -3.005054,53.439002,0.0 -3.000359,53.439912,0.0 -3.000583,53.440287,0.0 -2.997346,53.440855,0.0 -2.9974,53.440947,0.0 -2.993993,53.441549,0.0 -2.992722,53.440295,0.0 -2.991689,53.440711,0.0 -2.991084,53.440188,0.0 -2.990493,53.440356,0.0 -2.990529,53.440588,0.0 -2.979872,53.442502,0.0 -2.979102,53.443306,0.0 -2.977634,53.442823,0.0 -2.974937,53.443332,0.0 -2.976007,53.444715,0.0 -2.97629,53.445262,0.0 -2.974999,53.446484,0.0 -2.974857,53.447385,0.0 -2.973437,53.447316,0.0 -2.973428,53.447785,0.0 -2.973775,53.449643,0.0 -2.972951,53.449692,0.0 -2.97298,53.449867,0.0 -2.972666,53.450457,0.0 -2.971614,53.452187,0.0 -2.971732,53.456193,0.0 -2.973435,53.46015,0.0 -2.973839,53.462714,0.0 -2.969575,53.465769,0.0 -2.969663,53.465859,0.0 -2.969578,53.46592,0.0 -2.969853,53.465944,0.0 -2.967359,53.467739,0.0 -2.966448,53.467817,0.0 -2.966339,53.467743,0.0 -2.963667,53.467923,0.0 -2.962821,53.468017,0.0 -2.961919,53.468241,0.0 -2.961735,53.468119,0.0 -2.960256,53.467949,0.0 -2.960341,53.467635,0.0 -2.959328,53.467643,0.0 -2.958559,53.468798,0.0 -2.95817,53.468718,0.0 -2.9572,53.47075,0.0 -2.956333,53.472971,0.0 -2.955776,53.472797,0.0 -2.955724,53.472403,0.0 -2.955166,53.472404,0.0 -2.954904,53.472659,0.0 -2.953994,53.472542,0.0 -2.953036,53.472757,0.0 -2.947897,53.473048,0.0 -2.947174,53.472178,0.0 -2.946174,53.472959,0.0 -2.939856,53.473354,0.0 -2.939865,53.47356,0.0 -2.939051,53.473625,0.0 -2.938246,53.47359,0.0 -2.933399,53.473813,0.0 -2.933291,53.473095,0.0 -2.929877,53.473624,0.0 -2.927723,53.473805,0.0 -2.928387,53.473444,0.0 -2.92765,53.47332,0.0 -2.922593,53.474976,0.0 -2.921212,53.475508,0.0 -2.919153,53.476761,0.0 -2.918951,53.477114,0.0 -2.919146,53.477806,0.0 -2.919106,53.478046,0.0 -2.918442,53.478125,0.0 -2.917992,53.478305,0.0 -2.917247,53.479099,0.0 -2.916975,53.479146,0.0 -2.916838,53.479276,0.0 -2.917187,53.480115,0.0 -2.919585,53.48205,0.0 -2.917227,53.482278,0.0 -2.914885,53.482398,0.0 -2.913795,53.482649,0.0 -2.912705,53.483164,0.0 -2.912205,53.483657,0.0 -2.91225,53.484216,0.0 -2.911695,53.485177,0.0 -2.911211,53.485751,0.0 -2.910712,53.48663,0.0 -2.90895,53.488823,0.0 -2.907143,53.490033,0.0 -2.905172,53.492011,0.0 -2.903678,53.49169,0.0 -2.899819,53.491659,0.0 -2.899466,53.491739,0.0 -2.899473,53.492192,0.0 -2.899809,53.492365,0.0 -2.899347,53.493669,0.0 -2.898825,53.494125,0.0 -2.89797,53.494441,0.0 -2.897963,53.494837,0.0 -2.897303,53.494941,0.0 -2.896775,53.4954,0.0 -2.896588,53.495934,0.0 -2.896753,53.496265,0.0 -2.896411,53.496679,0.0 -2.897077,53.497553,0.0 -2.897848,53.498079,0.0 -2.897991,53.49851,0.0 -2.898604,53.499143,0.0 -2.898842,53.499649,0.0 -2.898815,53.499843,0.0 -2.898551,53.500151,0.0 -2.898488,53.500744,0.0 -2.897817,53.501163,0.0 -2.897766,53.50138,0.0 -2.89545,53.502327,0.0 -2.894304,53.502613,0.0 -2.893806,53.502916,0.0 -2.893012,53.502967,0.0 -2.892321,53.503149,0.0 -2.891336,53.502909,0.0 -2.889697,53.502816,0.0 -2.888014,53.50381,0.0 -2.888206,53.504217,0.0 -2.888924,53.505141,0.0 -2.889016,53.505735,0.0 -2.889729,53.507136,0.0 -2.892103,53.510805,0.0 -2.892387,53.510989,0.0 -2.892414,53.511124,0.0 -2.887744,53.514539,0.0 -2.885885,53.517333,0.0 -2.885558,53.518151,0.0 -2.884605,53.519129,0.0 -2.884081,53.519539,0.0 -2.883361,53.519554,0.0 -2.881233,53.519981,0.0 -2.881741,53.520522,0.0 -2.88203,53.520483,0.0 -2.882823,53.5206,0.0 -2.883971,53.52063,0.0 -2.886327,53.521058,0.0 -2.887464,53.521476,0.0 -2.887935,53.521739,0.0 -2.888207,53.521738,0.0 -2.888384,53.521876,0.0 -2.88898,53.522035,0.0 -2.889111,53.522189,0.0 -2.889475,53.522174,0.0 -2.889647,53.522267,0.0 -2.889868,53.522126,0.0 -2.89013,53.522158,0.0 -2.890382,53.521935,0.0 -2.890946,53.52195,0.0 -2.891189,53.522144,0.0 -2.891389,53.522182,0.0 -2.891688,53.522129,0.0 -2.892062,53.521928,0.0 -2.892502,53.522046,0.0 -2.893365,53.521789,0.0 -2.893989,53.52174,0.0 -2.894328,53.521587,0.0 -2.894696,53.521607,0.0 -2.894768,53.52144,0.0 -2.895137,53.521333,0.0 -2.895712,53.521565,0.0 -2.895746,53.521739,0.0 -2.896003,53.521832,0.0 -2.896163,53.52178,0.0 -2.896188,53.521531,0.0 -2.896511,53.521538,0.0 -2.896713,53.521737,0.0 -2.896699,53.52195,0.0 -2.896828,53.52197,0.0 -2.897153,53.521843,0.0 -2.897205,53.521667,0.0 -2.897472,53.521655,0.0 -2.897878,53.522095,0.0 -2.898161,53.522042,0.0 -2.898259,53.521892,0.0 -2.898534,53.521775,0.0 -2.898493,53.521585,0.0 -2.898596,53.521567,0.0 -2.899122,53.521828,0.0 -2.89916,53.521972,0.0 -2.899779,53.522045,0.0 -2.901564,53.521616,0.0 -2.906293,53.521708,0.0 -2.907245,53.521982,0.0 -2.907787,53.522417,0.0 -2.908567,53.522469,0.0 -2.908775,53.522351,0.0 -2.909224,53.522514,0.0 -2.909375,53.5227,0.0 -2.910289,53.522666,0.0 -2.910842,53.522775,0.0 -2.910843,53.522904,0.0 -2.910974,53.522957,0.0 -2.911315,53.522746,0.0 -2.911555,53.522776,0.0 -2.911736,53.522903,0.0 -2.911448,53.523031,0.0 -2.911344,53.523185,0.0 -2.91158,53.523321,0.0 -2.912084,53.523161,0.0 -2.912207,53.523256,0.0 -2.91221,53.523447,0.0 -2.912701,53.523369,0.0 -2.912653,53.523506,0.0 -2.91274,53.523537,0.0 -2.912991,53.523422,0.0 -2.913317,53.523603,0.0 -2.913476,53.523365,0.0 -2.913845,53.523476,0.0 -2.913907,53.523412,0.0 -2.914092,53.523556,0.0 -2.914376,53.523563,0.0 -2.914359,53.523368,0.0 -2.914468,53.523281,0.0 -2.914931,53.52353,0.0 -2.915257,53.523468,0.0 -2.915454,53.523638,0.0 -2.915397,53.523704,0.0 -2.915685,53.52368,0.0 -2.915904,53.523976,0.0 -2.916348,53.523853,0.0 -2.916361,53.523681,0.0 -2.916782,53.5235,0.0 -2.916742,53.523182,0.0 -2.916918,53.523145,0.0 -2.917163,53.523478,0.0 -2.916958,53.523817,0.0 -2.917012,53.523952,0.0 -2.919085,53.524357,0.0 -2.920792,53.52493,0.0 -2.921877,53.525451,0.0 -2.922459,53.525514,0.0 -2.922845,53.525376,0.0 -2.922996,53.525235,0.0 -2.92386,53.525569,0.0 -2.924108,53.525764,0.0 -2.924118,53.525919,0.0 -2.923966,53.526027,0.0 -2.924865,53.526221,0.0 -2.925641,53.526888,0.0 -2.925808,53.527293,0.0 -2.925899,53.528769,0.0 -2.926315,53.529478,0.0 -2.926962,53.529989,0.0 -2.928123,53.530485,0.0 -2.92905,53.530515,0.0 -2.929344,53.530596,0.0 -2.929931,53.531072,0.0 -2.930518,53.531263,0.0 -2.930855,53.531168,0.0 -2.931259,53.531225,0.0 -2.931357,53.531422,0.0 -2.932688,53.53245,0.0 -2.9329,53.532377,0.0 -2.933093,53.532427,0.0 -2.933222,53.532764,0.0 -2.933432,53.532891,0.0 -2.933813,53.533394,0.0 -2.933836,53.533586,0.0 -2.933692,53.533789,0.0 -2.933754,53.534075,0.0 -2.933611,53.534415,0.0 -2.934054,53.534883,0.0 -2.934576,53.535045,0.0 -2.934452,53.5355,0.0 -2.934659,53.53573,0.0 -2.93463,53.535837,0.0 -2.934841,53.535943,0.0 -2.935127,53.535921,0.0 -2.935247,53.536026,0.0 -2.935282,53.536166,0.0 -2.93509,53.536224,0.0 -2.935315,53.536563,0.0 -2.935081,53.536589,0.0 -2.934963,53.536489,0.0 -2.934775,53.536664,0.0 -2.935156,53.537025,0.0 -2.935404,53.537138,0.0 -2.935901,53.536999,0.0 -2.936029,53.537059,0.0 -2.936384,53.537029,0.0 -2.936162,53.537491,0.0 -2.936402,53.537538,0.0 -2.936412,53.537746,0.0 -2.93598,53.537887,0.0 -2.935957,53.537973,0.0 -2.93612,53.538238,0.0 -2.936081,53.538396,0.0 -2.936773,53.538697,0.0 -2.937655,53.538535,0.0 -2.937918,53.538841,0.0 -2.938232,53.538508,0.0 -2.938416,53.538714,0.0 -2.938607,53.538759,0.0 -2.938744,53.538502,0.0 -2.939063,53.538737,0.0 -2.939084,53.539092,0.0 -2.939651,53.538995,0.0 -2.939926,53.539304,0.0 -2.94041,53.5394,0.0 -2.940709,53.539608,0.0 -2.940873,53.539495,0.0 -2.941262,53.539877,0.0 -2.941567,53.539685,0.0 -2.94145,53.540326,0.0 -2.941606,53.540666,0.0 -2.941888,53.540679,0.0 -2.942013,53.540994,0.0 -2.943086,53.541068,0.0 -2.943685,53.541191,0.0 -2.943925,53.540939,0.0 -2.944725,53.54122,0.0 -2.9449,53.541339,0.0 -2.944613,53.541604,0.0 -2.944753,53.541723,0.0 -2.944743,53.541876,0.0 -2.945134,53.542046,0.0 -2.945223,53.542508,0.0 -2.945478,53.542808,0.0 -2.945925,53.542987,0.0 -2.94718,53.54376,0.0 -2.947876,53.544319,0.0 -2.948319,53.544356,0.0 -2.950904,53.544269,0.0 -2.951182,53.544502,0.0 -2.952022,53.544896,0.0 -2.953665,53.545277,0.0 -2.95408,53.545569,0.0 -2.954628,53.545665,0.0 -2.95504,53.546054,0.0 -2.956041,53.545892,0.0 -2.956182,53.545901,0.0 -2.95637,53.546072,0.0 -2.956607,53.546083,0.0 -2.9569,53.545839,0.0 -2.95747,53.54575,0.0 -2.957769,53.545604,0.0 -2.958192,53.545693,0.0 -2.958665,53.545571,0.0 -2.959714,53.545547,0.0 -2.96054,53.545948,0.0 -2.961249,53.546025,0.0 -2.963006,53.546393,0.0 -2.963226,53.546349,0.0 -2.963271,53.546705,0.0 -2.963343,53.546759,0.0 -2.963539,53.54663,0.0 -2.963623,53.546685,0.0 -2.963579,53.547005,0.0 -2.96381,53.547071,0.0 -2.964109,53.546999,0.0 -2.964372,53.547029,0.0 -2.965169,53.547648,0.0 -2.965653,53.547562,0.0 -2.965746,53.547607,0.0 -2.965873,53.547868,0.0 -2.966788,53.547542,0.0 -2.967211,53.547467,0.0 -2.967509,53.547607,0.0 -2.967918,53.5475,0.0 -2.967978,53.547339,0.0 -2.968453,53.547421,0.0 -2.968628,53.547314,0.0 -2.96912,53.547356,0.0 -2.969138,53.547252,0.0 -2.969943,53.54722,0.0 -2.970178,53.547074,0.0 -2.970416,53.547107,0.0 -2.970549,53.547242,0.0 -2.970206,53.547297,0.0 -2.970322,53.547503,0.0 -2.970612,53.54748,0.0 -2.970795,53.547718,0.0 -2.971294,53.54762,0.0 -2.971415,53.547706,0.0 -2.971427,53.547906,0.0 -2.971586,53.547986,0.0 -2.972061,53.547882,0.0 -2.972575,53.54794,0.0 -2.972958,53.547791,0.0 -2.973302,53.547869,0.0 -2.973263,53.547673,0.0 -2.974094,53.547551,0.0 -2.974178,53.547393,0.0 -2.974468,53.547363,0.0 -2.974569,53.547202,0.0 -2.974943,53.54727,0.0 -2.975188,53.547231,0.0 -2.975542,53.546973,0.0 -2.97558,53.54682,0.0 -2.976509,53.546612,0.0 -2.977229,53.546262,0.0 -2.977396,53.546018,0.0 -2.977188,53.545831,0.0 -2.977641,53.545902,0.0 -2.977794,53.545299,0.0 -2.978119,53.545004,0.0 -2.978487,53.543975,0.0 -2.978348,53.542732,0.0 -2.977547,53.542199,0.0 -2.977282,53.541717,0.0 -2.976838,53.541361,0.0 -2.976634,53.54079,0.0 -2.97612,53.540353,0.0 -2.975621,53.540257,0.0 -2.975523,53.540122,0.0 -2.974719,53.539817,0.0 -2.974391,53.539611,0.0 -2.974068,53.538664,0.0 -2.973702,53.538015,0.0 -2.97281,53.537102,0.0 -2.972846,53.536123,0.0 -2.971895,53.535162,0.0 -2.971278,53.534807,0.0 -2.970501,53.534131,0.0 -2.97033,53.533334,0.0 -2.970048,53.532851,0.0 -2.969679,53.532476,0.0 -2.969211,53.531515,0.0 -2.968591,53.530664,0.0 -2.968481,53.53002,0.0 -2.967938,53.529576,0.0 -2.967764,53.529299,0.0 -2.967134,53.528191,0.0 -2.966849,53.527372,0.0 -2.965091,53.526349,0.0 -2.96483,53.526007,0.0 -2.964184,53.525603,0.0 -2.96405,53.525381,0.0 -2.963275,53.524815,0.0 -2.963233,53.524316,0.0 -2.962991,53.523684,0.0 -2.963025,53.522881,0.0 -2.963464,53.521956,0.0 -2.96364,53.520497,0.0 -2.964197,53.51964,0.0 -2.964048,53.519186,0.0 -2.96417,53.518584,0.0 -2.964699,53.518508,0.0 -2.965234,53.518571,0.0 -2.96724,53.517954,0.0 -2.968171,53.517885,0.0 -2.968382,53.517688,0.0 -2.970205,53.51684,0.0 -2.97121,53.51628,0.0 -2.97169,53.516082,0.0 -2.972473,53.516032,0.0 -2.972971,53.515424,0.0 -2.973772,53.515505,0.0 -2.976101,53.51528,0.0 -2.976457,53.51547,0.0 -2.977009,53.516319,0.0 -2.980611,53.518771,0.0 -2.982029,53.51922,0.0 -2.982559,53.519523,0.0 -2.983437,53.520671,0.0 -2.984,53.521052,0.0 -2.984761,53.521371,0.0 -2.985724,53.52155,0.0 -2.986923,53.521547,0.0 -2.988192,53.521364,0.0 -2.989623,53.521625,0.0 -2.989813,53.521729,0.0 -2.990118,53.522532,0.0 -2.990351,53.522759,0.0 -2.99155,53.52297,0.0 -2.991738,53.52267,0.0 -2.995469,53.519657,0.0 -2.994792,53.519202,0.0 -2.994826,53.519062,0.0 -2.996748,53.518037,0.0 -2.99727,53.517587,0.0 -2.997622,53.517567,0.0 -2.99785,53.517475,0.0 -2.998144,53.516609,0.0 -2.999059,53.515647,0.0 -2.999584,53.514869,0.0 -3.001604,53.513293,0.0 -3.00351,53.512459,0.0 -3.003661,53.511795,0.0 -3.004208,53.511303,0.0 -3.004672,53.511159,0.0 -3.00497,53.510895,0.0 -3.005114,53.510539,0.0 -3.005069,53.510264,0.0 -3.005763,53.510426,0.0 -3.006054,53.510415,0.0 -3.006916,53.510222,0.0 -3.009225,53.510039,0.0 -3.009696,53.509752,0.0 -3.012328,53.508992,0.0 -3.013712,53.508723,0.0 -3.014604,53.509456,0.0 -3.019221,53.512277,0.0 -3.043811,53.525895,0.0 -3.045598,53.526426,0.0 -3.04884,53.528352,0.0 -3.052484,53.529852,0.0 -3.053508,53.530078,0.0 -3.054894,53.530191,0.0 -3.058178,53.530912,0.0 -3.058396,53.531036,0.0 -3.06001,53.530376,0.0 -3.0595,53.530035,0.0 -3.058535,53.528273,0.0 -3.058617,53.528126,0.0 -3.059851,53.527832,0.0 -3.062978,53.527274,0.0 -3.063148,53.526595,0.0 -3.063203,53.526485,0.0 -3.063385,53.526459,0.0 -3.108843,53.525486,0.0 -3.106121,53.524247,0.0</coordinates></LinearRing></outerBoundaryIs></Polygon>      Am I missing anything? The documentation for the ST_GeomFromKML is not very useful for troubleshooting this. **Note:** The same happens when exporting to KML using ST_asKML, i.e.: the following script also throws an "invalid KML representation" error `SELECT ST_GeomFromKML((SELECT ST_asKML((SELECT shape FROM geographies LIMIT 1))));` **Edit:** Added full KML shape


Are polygons created in Arcpy built in any order?

My script uses a search cursor to evaluate polyline coordinates and build polygons based on those coordinates. The script works well, but I am curious as to whether anyone knows the order (if there is any) that polygons are built. I know that the output polygon feature class reads the rows (and their coordinates) in order, but I do not know if they _build_ them in order. I could see the polygons being built according to which order they were read, or based on size, etc. Any input would be appreciated. Thanks in advance.


Is there any tool to reproject x,y coordinates, apart from Proj4?

I have a JEE project (Java, JBoss) and I'm looking for a library that allows x,y coordinates reprojection. I know about Proj4. Do you have any other option? Is it possible to find on the net a reprojection algorithm that I could use? The reprojection I'm interest in is the French NTF Lambert II Etendu to RGF Lambert 93.


Downloading multiple files from TIGER FTP site

![TIGER FTP Site](http://i.stack.imgur.com/zpIOq.png) Hello everyone! I've been downloading TIGER data per county from the FTP site one at a time, with some states having hundreds of files. I was curious if there was a faster way to download maybe all the files for each county for each state at once? Thanks for your time and help! Attached is a picture of what i'm dealing with, would like to eventually download all the files for each state at once versus one at a time.


Joining features based on location and attribute at the same time

I have a task of assigning data from database that I represent as point layer, and containing information about characteristics of soils through the world, and also soil type. I also have polygon layers for each lattitude, that contain area of each type of soil (same ones as point layer, therefore can be ussed as join feature). Inside each lattitude layer, each soil type is only one, multipart polygon. I need to assign values from point layer in a way that the soil type always match, but also that it uses point closest geograpicly. If I use just join by attributes, it doesnt take into account spatial location, and vice versa, and I really have trouble figuring how I could approach that. Each layaer has around 40 soil types, and so doing it manually would take ages. I can also add that majority of points are way off the polygons that I intend them to join to, but as long as the soil types match and points is closest, its what I need. I have basicly no knowledge about python scripting. I tried making model to handle that, but its also beyond me. Thanks in advance


TileMill Exports Failing

I'm having trouble exporting my map in TileMill (v0.10.1). I'm using OSM Bright with a PostGIS database connection to import OpenStreetMap data. When exporting an SVG or PDF through the user interface in TileMill, the export fails more often than not with a SIGABRT error. The console shows the following error:               [tilemill] [tilemill-ui] Rendering file     [tilemill] [tilemill-ui] libc++abi.dylib: terminating     [tilemill] Export process failed with signal 'SIGABRT'      I then tried exporting through the command line by running `./index.js export OSMBright ~/Documents/MapBox/export/cityhall2.svg --format=svg --width=2000 --height=2300 --bbox="-79.8786965, 43.2604005, -79.8702965, 43.2516005"`, with a nearly identical error. It usually exports a blank file, although once in a while it does export nearly what I want (with zoom levels and symbology slightly off, but still usable). Has anyone had this problem before? Thanks for the help!


Ecoregions shapefiles: subsetting, extracting and summarizing data from raster

I am attempting to subset and project an ecoregions shapefile from : EPA ECOREGIONS my data set is found here: SHAPEFILES and RASTER When I run this code, two major problems occur: 1) I cannot get the subset of ecoregions to plot on top of my raster file and 2) I cannot extract the values (summary information) from the subset of ecoregions from the raster. Can anyone help with this?               require (raster) raster<-raster('FG_RCP26_2050.tif')      plot (raster)      ![Raster with Data](http://i.stack.imgur.com/5UQMj.jpg) ecoregions.all<\- readOGR("NA_CEC_Eco_Level2.shp",layer="NA_CEC_Eco_Level2") head (ecoregions.all) plot (ecoregions.all) ![ecoregions](http://i.stack.imgur.com/DthVm.jpg)                   #this is the subset of ecoregions that I need        NA_L2<-c('5.1','5.2','5.3','8.1','8.3','8.4','8.5')                     #subset the ecoregion.all to just those in the northeast      #this must be where the problem is      ne.ecoregions = coregions.all[match(toupper(NA_L2),toupper(ecoregions.all$NA_L2CODE)),]                    #plot ne.ecoregions on top of raster and extract summary statistics for each ecoregion      #summary statistics =distribution of values within each ecoregion     plot (raster) plot (ne.ecoregions, add=T)          ecoregions.raster<-mask (raster, ne.ecoregions)      summary (ecoregions.raster)


JSTS : Union on polygons with holes

![JSTS union](http://i.stack.imgur.com/CcUPL.png) Hi all, I'm trying to union two polygons (see above) with JSTS. My first polygon contains an hole, and my second polygon fill partially this hole. The expected result is a unique polygon with two hole (or one hole with two points at the same location? I've tried with PostGIS and the result is two holes). The actual result is what you can see on C: my first polygon has no hole anymore. Is it a defect of the JSTS tool? Do you have an idea on how can I fix it? Or is it just some configuration to add somewhere? I'm using the JSTS union function like that: firstGeometry.union(secondGeometry); Please ask if you need more information. Thank you for time reading me! Sébastien EDIT: Thanks for your answers! @WhiteboxDev : I don't think I have topological errors. I'm easily able to reproduce this case by drawing polygons on OpenLayers, and my polygons seems to be correct. Please note that I'm using JSTS and not JTS. I've tried with the buffer(0,0) function, but I have unfortunately the same result. Here is an example of a union that fails: POLYGON((335179.0532258 328356.69669355,335178.1782258 328437.63419355,335285.3657258 328437.19669355,335284.9282258 328355.82169355,335179.0532258 328356.69669355),(335213.6157258 328373.32169355,335198.3032258 328400.88419355,335225.4282258 328426.69669355,335267.4282258 328415.75919355,335262.1782258 328371.13419355,335213.6157258 328373.32169355)) POLYGON((335213.6157258 328373.32169355,335225.4282258 328426.69669355,335267.4282258 328415.75919355,335213.6157258 328373.32169355)) And the result of a JSTS union() between these two geometries or a buffer(0,0) on a new geometry collection containing these geometries is the following polygon : POLYGON((335179.0532258 328356.69669355,335178.1782258 328437.63419355,335285.3657258 328437.19669355,335284.9282258 328355.82169355,335179.0532258 328356.69669355)) EDIT (2): It happens in each case where I have   * 1 to many hole(s)   * 1 to many polygon(s) within this/these hole(s), creating "sub-holes"   * Sub-holes are sharing a vertex Doesn't seems to happen in other cases. Geometries are topologically correct. EDIT (3): I've found a defect about that on GitHub: https://github.com/bjornharrtell/jsts/issues/135 I'm still looking for a way to fix it!


Editor and projections ArcGIS

When opening an edit session Arc is telling me that some of my shape files' coordinate systems do not match the data frame. Here is the data frame coordinate system: > Stereo_70 > > Projection: Double_Stereographic > > False_Easting: 500000.000000 > > False_Northing: 500000.000000 > > Central_Meridian: 25.000000 > > Scale_Factor: 0.999750 > > Latitude_Of_Origin: 46.000000 > > Linear Unit: Meter > > GCS_Dealul_Piscului_1970 > > Datum: D_Dealul_Piscului_1970 And here is the coordinate system that all the shape files in question have: > Projected Coordinate System: Stereo_70 > > Projection: Double_Stereographic > > False_Easting: 500000.00000000 > > False_Northing: 500000.00000000 > > Longitude_Of_Origin: 25.00000000 > > Scale_Factor: 0.99975000 > > Latitude_Of_Origin: 46.00000000 > > Linear Unit: Meter > > Geographic Coordinate System: GCS_Dealul_Piscului_1970 > > Datum: D_Dealul_Piscului_1970 > > Prime Meridian: Greenwich > > Angular Unit: Degree Now is it just me or is anyone else having difficulty seeing the difference? Doe this mean I'm going to have to reproject the files in question FROM their original projection, TO their original projection? I'm not even sure that's possible but will try now. And short of that I guess I just wipe the projections and define them again.


ArcGIS Server 9.3 - Javascript API - Selecting Feature Problem

I have ArcGIS Server 9.3, and I'm trying to use the JS API. I successfully displayed the map, and I can locate and zoom to a feature based on ParcelID value. The problem is I can't highlight this parcel. I mean I tried to configure the symbol, but its not being displayed. Here's the code:               <script>              dojo.require("esri.map");         dojo.require("esri.layers.FeatureLayer");         var map,featureLayer,query,sfs;            function init() {         map = new esri.Map("mapDiv", {            sliderOrientation: "horizontal"         });                    var imageParameters = new esri.layers.ImageParameters();         imageParameters.format = "PNG24";  //set the image type to PNG24, note default is PNG8.              //Takes a URL to a non cached map service.         var dynamicMapServiceLayer = new esri.layers.ArcGISDynamicMapServiceLayer("http://*******/*****/rest/services/****/MapServer", {           "opacity":1,            "imageParameters":imageParameters         });         dynamicMapServiceLayer.setVisibleLayers([4,10,11,16]);         map.addLayer(dynamicMapServiceLayer);                   //Locate Parcel--------         sfs = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_SOLID,new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_DASHDOT,new dojo.Color([255,0,0]), 2),new dojo.Color([255,255,0,0.25]));              featureLayer  = new esri.layers.FeatureLayer("http://*****/*****/rest/services/*****/MapServer/4", {                 outFields: ["*"],                 mode: esri.layers.FeatureLayer.MODE_SELECTION             });         featureLayer.setSelectionSymbol(sfs);         map.addLayers([featureLayer]);              query = new esri.tasks.Query();                 query.where = "ParcelID_NS = '101010101'";              var deferred = featureLayer.selectFeatures(query, esri.layers.FeatureLayer.SELECTION_NEW, function (selection) {             var stateExtent = selection[0].geometry.getExtent().expand(3.0);             map.setExtent(stateExtent);         });         //----------------------           }                 dojo.ready(init);          </script>      Any help? Thanks.


Why does the attribute Table not maintain input data?

for a longer time I have a problem with the attribute table. I create the polygons on the map, I fill up the table with attributes. Then, if I click "save", all the data from the table are being deleted and in all fields of the table a word "error" does appear. After the restart of the program, only that columns containing the number values are being kept. That ones which contain the words are deleted and the word "null" does appear in all the fields. Is somebody able to help me?


Which SRS is computational more effective to store continental scale data? (given some conditions)

For a Precision Agriculture application I need to make calculations in UTM projection. The data will be stored in Postgis and is made of soil samples results (point data) and field boundaries of many farms. But the farms are scattered through many utm zones (ex. a farm is in Z23 and the other is in Z22, so I can't store the data in UTM, it needs to be stored in a SRS without zones. Since the calculations will be made in UTM the data will always need to be converted back and forth to the storage SRS. So which of those inputs is cheaper in terms of computations of SRS conversion when the output SRS is UTM Wgs84:   -Geometry vs. Geography Postgis types   -Geometry projected (ex World Mercator) vs. Unprojected (latlon)  I think that it involves the internal Proj4 math but I couldn't understand it myself...


Is 45005000 a recognisable location?

Hello and sorry if this is not a relevant question to this SE. Please if it is not, mitigate it to where you see fit. I have recently purchased an item online and I have been provided by the seller with a tracking number applicable to this site: http://track-chinapost.com/startairmail.php . So when I check where my item is, I get Location: 45005000 Any ideas what this could mean?


Geometry issues when creating table from existing spatial table

In short, I want to split an existing spatial table into two parts. The split itself is straight forward. I start with this (fully functioning) table:               CREATE TABLE table     (       my_identifier smallint,       geom geometry,       CONSTRAINT enforce_dims_geom CHECK (st_ndims(geom) = 2),       CONSTRAINT enforce_srid_geom CHECK (st_srid(geom) = 4326)     )     WITH (       OIDS=TRUE     );      and create the new tables by:               CREATE TABLE public.table_part1 AS     SELECT * FROM public.table WHERE idmy_identifier = 1;          CREATE TABLE public.table_part2 AS     SELECT * FROM public.table WHERE my_identifier = 2;      However, the problem is that although the geom column gets transferred to my new tables, the geometry (?) itself doesn't get transferred correctly:               CREATE TABLE table_part1     (       my_identifier smallint,       geom geometry     )     WITH (       OIDS=FALSE     );      I tried running `SELECT Probe_Geometry_Columns();` but it didn't help. Do I need to insert data into the geometry_columns table? If so, exactly what?


Delete polygon underlap?

I am editing a layer in QGIS with severe polygon overlaps. Is there a way to delete the underlapping portions (i.e. the invisible parts) using QGIS? If I import the layer into GRASS using the GRASS-plugin, I can get a very useful layer which contains the overlaps (polygons_2), but if I use the difference geoprocessing tool to slice the overlaps from the original file, it naturally cuts away the overlap as well as the underlap, leaving holes instead of overlaps. I have also tried the v.clean functions in grass, but the overlaps are really too big to use snapping without distorting the data.


How to transform Invalid Date format in Data Interoperability

My workbench is running into some date data it is rejecting. Here is the error I'm receiving... The invalid value '20100419000000' was supplied to the field 'install_date_orig' I believe the feature is rejected because of the '000000' for the time portion of the date. Any suggestions of which transformers may be useful to get around this issue. Here's an list of of what I've tried so far.   * DateFormatter with this format, FME Date (%Y%m%d)   * SubStringExtractor to split the field into the year, month and day parts, StringConcatenator to recombine those parts, the DateFormatter as per above, and finally AttributeRemover to exclude the bad field.   * TestFilter and Tester to filter out the bad records.  I've noticed that when it reaches the bad record it fails before it reaches the first transformer. Is there anyway to catch a record with a bad field value before it reaches any of the transformers? I appreciate any suggestions.


Is the "sdetable" command available with ArcGIS Workgroup, or only the Enterprise edition? We have ArcGIS Server Workgroup 9.3.1 and want to update some of the feature class tables with data from another system using SQL UPDATE commands. From taking a look at the documentation, we must first create a multi-versioned view, and then run the SQL UPDATE commands against this view. The problem is that we do not have the "sdetable" command required to create the multi- versioned view. We think that these command line tools may be available only in the Enterprise edition. Is there another way that we can perform basic data updates with SQL commands in the Workgroup Edition? I understand that the correct way to do this is not by editing the default version directly. EDIT:   Our data is versioned for different uses, including some field applications. We were planning on creating a "sql-updates" version just for making updates to attributes in a feature class. If I use SQL UPDATE commands to update this version, my concern is that the changes will not get replicated to the other versions, or the data could become corrupt. This seems like it would be bypassing delta table entries. From what I have read, we should use "sdetable" to create a multiversioned view, and then use stored procedures create_version, set_current_version, edit_version (1 - start session), run my updates, edit_version (2 - close session), and finally delete_version (after reconciled). I see these stored procedures in the database, but can't do the first step of creating the multiversioned view. Maybe it would be possible to perform these steps on the "sql-updates" version without creating the view. **EDIT 2:**   Since the original question has been answered, I opened a new question:   How can we use SQL to update feature class attributes in a versioned geodatabase in ArcGIS Workgroup


Displaying road outlines but not segments i.e. no lines at intersections?

![enter image description here](http://i.stack.imgur.com/uuE0F.png) I have the roads that I want to just display at outlines but not segments (I.e not disconnected from each road). I am trying dissolve but its not working. Attached is a link to a package file of the data. Please help! https://www.dropbox.com/s/55hqr9471v9qxwf/Road%20Test%20package.mpk


QGIS composer - How can I eliminate striping across transparent raster layers?

I have downloaded QGIS dev version 1.7.0-114 and I have come across some weird behaviour in the composer. I create a map with a hillshade and a transparent DEM ontop of it, pretty standard stuff. When I export this map to an image format I get two strips across my output, this happens for a number of formats (png, jpg, bmp, tiff) and the stripes seem to be in the same place everytime. I have attached an example below (sorry for the size of the attachment). ![enter image description here](http://i.stack.imgur.com/RKFhO.jpg) The stripes appear in the same position on the page every time and this happens with different raster datasets, so I have ruled out the data. Also depending on the page size you may get 1 stripe (A5 as shown above), 2 stripes (A4) or no strips (custom dimensions h100 x w100). Is this something that can be overcome via some settings which I may have overlooked? Or might this be something funny in the software?


how to create Multiple layers?

I'd like to create multiple layers for asset management using openlayers, could anybody help?


Can ArcGIS layer be shown on Google Earth?

If we have a table in SQL Server which was used as layer in ArcGIS. Can the same thing be done in Google Earth? I mean does Google Earth support adding custom layers? If so how?


Open Layers 2.13.1 Image Layer Zoom Animation Issue

I updated my web app to use 2.13.1 from 2.12 and see that Image layers now animate zoom level changes? I don't believe it previously did this and investigating the code only briefly, I'm not sure how this even happens. I also do not see any obvious way to specify that zoom level changes should animate or not with Image Layers. The reason I care is that I'm using a Google base layer with an Image layer on top, and zoom changes animate for the Image Layer but not really the same way with the Google layer. The result is that the Image overlay does not appear to be properly attached to the Google layer through the zoom animation and it looks wrong. I experimented with the "animationEnabled" option on the Google layer but it doesn't improve the situation I'm not sure that I've explained this very well so I created a quick screen cap move of it in action. Unfortunately screen recorded slowed my computer down so the movie isn't as smooth as what I see without screen recording but it still shows the issue. http://youtu.be/i53UoxOmVQc I simply want all layers to have synchronized zoom animations. Any ideas how to fix?


Geoserver WPS - Tiled Barnes Surface with geospatial index

I'm using Barnes Surface to compute a set of points from my PostgreSQL db. The problem is that my whole table is processed for each tile. I added a geospatial index on my points, how can I tell geoserver to compute only the points present in the bounds of the tile?


PGRouting Layer plugin error in QGIS 2.2

I am trying to do some basic computations using the PGRouting plugin for QGIS, but I am getting this error all time:               An error has occured while executing Python code:          Traceback (most recent call last):       File "C:/Users/Gonzalo/.qgis2/python/plugins\pgRoutingLayer\pgRoutingLayer.py", line 407, in run         if db and db.con:     UnboundLocalError: local variable 'db' referenced before assignment          Python version:     2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)]               QGIS version:     2.2.0-Valmiera Valmiera, c3a2817          Python path: ['C:/Users/Gonzalo/.qgis2/python/plugins\\processing', 'C:/PROGRA~1/QGISVA~1/apps/qgis/./python', u'C:/Users/Gonzalo/.qgis2/python', u'C:/Users/Gonzalo/.qgis2/python/plugins', 'C:/PROGRA~1/QGISVA~1/apps/qgis/./python/plugins', 'C:\\PROGRA~1\\QGISVA~1\\bin\\python27.zip', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\DLLs', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\plat-win', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\lib-tk', 'C:\\PROGRA~1\\QGISVA~1\\bin', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\PIL', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\win32', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\win32\\lib', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\Pythonwin', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\Shapely-1.2.18-py2.7-win-amd64.egg', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\six-1.3.0-py2.7.egg', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\xlrd-0.9.2-py2.7.egg', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\xlwt-0.7.5-py2.7.egg', '.', 'C:\\PROGRA~1\\QGISVA~1\\apps\\qgis\\python\\plugins\\fTools\\tools']      I made a local installation of postgres 9.3.4, later I added postgis using the stack builder, and finally I followed the instructions posted here: PGRouting 2.0 for windows. I also installed python-psycopg2 for my version of python. The SQL query runs well in pgAdminIII, but I can not make that QGIS plugin work, so any help would be really appreciated.


Shoehorning hierarchial data into MS SQL ArcGIS SDE

So, I've been tasked with creating a set of tables to store hierarchical data within an SDE. This is in reference to utility data, i.e, a set of valves, control panels, actuators, etc (each potentially having their own attributes) belong to a pump, the pump belongs to a process (also having its own attributes), which belongs to a pumpstation (which can have its own attributes) etc etc. The only way I can think of to do this is to basically have a table for each asset (a table each for valves, control panels, actuators, pumps, process, pumpstation) and a relationship class for each one to link to the table higher up in the hierarchy for it, all the way to the top, and then a final relationship class to link the pump station table to an already existing pump station point feature class. I've already googled some to find examples of using hierarchical data in ArcGIS and i'm not having much luck. Anyone have any suggestions if there's a better way of doing this?


small difference in projection between tilemill and postgis

I would like to create a world map with interaction on continents created directly from a postgis database. For the moment, I have a POSTgis database with the data. I created a background with TileMill. I get data from Postgis in SVG and put them on the map using Raphael. All seem well except that I have a small difference between both map :s I use default projection of TileMill (and I think we can't change it :)). And here is the POSTGIS command I'm executing to get the SVG:               SELECT st_asSVG(            st_scale(                st_translate(                    St_Intersection(                        ST_Transform(ST_SetSRID(ST_GeomFromText('POLYGON((-151.875 80.6470,175.7812 80.6470,175.7812 -63.2336,-151.875 -63.2336,-151.875 80.6470))'), 4326), 3857),                    (SELECT st_union(geom) FROM geo.continentboundaries))                , 20037508.72, -20037508.72)            , 0.0000199043047125,0.0000199043047125)         , 0, 0)      Here is the result: http://i.stack.imgur.com/7lS0F.png Do you know where I have a mistake? Edit: some precisions: I'm using the same data for Tilemill and the svg. I added in the tilemill project the example of tilemill there is no problem of projection. In the layer in tilemill I'm using 900913. st_srid give me 3857 in Postgis which seems ok? I'm a beginner in SIG


Exploding overlapping to new non-overlapping polygons?

Given multiple polygons that overlap in multiple ways, I would like to export from these features all polygons that don't overlap with others, iteratively. The product would be a number of features with no overlap that, when summed together, make up the original. The products could then be used as input to Zonal Statistics, and this would be much faster than iterating Zonal Statistics over each polygon. I have been trying to code this in ArcPy without success. Does code to do this already exist?


What is geom_way column in osm2po generated data?

Could anybody tell me the meaning of geom_way column generated by osm2po?


Looking for a functionality in QGIS to find all neigbours in a polygon layer

I am looking for a functionality in QGIS to find all neigbours for each polygon in a polygon layer. The result should be a table with 2 ID-Columns: one for the polygon and the second for the neighbours. I tried "spatial query" in QGIS with the layer as source and its copy as reference, but to no avail. In Mapbasic of Mapinfo it would be a simple expression:               select t1.id1, t2.id2 from t1, t2 where t1.obj intersects t2.obj       (where t1 and t2 are two copies of the polygon layer). Thanks for any help!


Is it possible to automatically set CRS of Rasters from PRJ file?

I am trying to learn my colleagues to display measured data (ASCII) in maps relatively easily. I can prepare project, color and symbology style for the point layer and background map. I also found good way to get background map (PNG + world file) for offline use via MOBAC and OSM. MOBAC uses EPSG:3395 WGS 84 / World Mercator and I have the appropriate PRJ file. However, I found out, that QGIS does not use the PRJ for rasters and load them always using default CRS (in my case EPSG:4326). Is it possible to make QGIS load the CRS settings from the PRJ (or other file) without need to set the CRS manually for every raster? thanks


How to produce "SpatialPixelsDataFrame" from stacked RasterLayer?

How to produce "SpatialPixelsDataFrame" from staccked RasterLayer ? I have a stacked rasterlayer and I want to convert them to SpatialPixelsDataFrame. RasterLayer information is               > RAD.all[[1]]     class       : RasterLayer      dimensions  : 180, 240, 43200  (nrow, ncol, ncell)     resolution  : 0.01, 0.01  (x, y)     extent      : -123.305, -120.905, 36.705, 38.505  (xmin, xmax, ymin, ymax)     coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0      data source : C:\Users\Michael\Documents\1998bil\1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil      names       : X1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc


Near analysis by matching attribute in two different layers

I need to perform a near analysis to assign cases in one feature class to their nearest match on an attribute in a different feature class (ex: Tutors and Tutees need to be paired by same language spoken). I have not been able to find a way to do this. The closest I have reached to simplifying this process is to perform a split by attribute on the values for the matching variable. So instead of having all tutors in one feature class, I now have n language number of feature classes containing tutors and n language number of feature classes containing tutees. I then plan on performing a near analysis between each feature class. However this is extremely tedious and may result in 100+ near analyses needing to be performed because there are other criteria I will need to match independently on in a similar way. Is there a better approach? Much thanks in advance. Edit: The near must be performed between two different layers to preserve the different roles of individuals being matched. I am new to learning modelbuilder so any ideas on conditionals, or iterators that could determine the "nearest" match on a shared attribute for points housed into two different layers would be much appreciated!


Mapping elements in a warehouse

I am new to GIS. For a small project I want to map elements in a warehouse and want to find the location of the elements later (e.g. nails are located at x,y). Now I am wondering what kind of position I should choose. Should I choose an absolute location coordinate system like `1,3`, or should I use real geocoordinates like `-7.023, 32.2423` for my use case? Where are the benefits/disadvantages with each approach?


Onclick event with polygons

When I click on any polygon in the openlayer I need to show the respective information in the popup. Right now I am using the below code.               eventListeners: {         getfeatureinfo: function(event) {             window.open("http://mydomain.com/polygoninfo.php?polygonid=" + event.feature.attributes.id, "window_name","location=1,status=1,scrollbars=1,resizable=no,width=650,height=650");         }     }      The issues with the above code are:   1. The popup is opening on complete layer, meaning that the popup is opening only on polygons if I click some other parts of the layer where the polygon is not present.   2. I need to capture the unique ID of the polygon when the polygon is clicked. I am using sld for labels, etc. Is there a way to create a uniqueId of polygon inside sld like:                       <sld:Label>                <ogc:PropertyName>CAGE_NO</ogc:PropertyName>              </sld:Label>      Sorry if any thing in my previous question is unclear or confusing. So, please let me elaborate on the issue as below:   1. As per my understanding, the actual requirement/task is to provide a layer which represents various states in the India map. Once the user clicks on any of the polygons/states, then the request should be redirected to the respective state's home page.   2. My concern is when a user clicks on any of the polygons/states, how to capture or send the respective state's ID/name so that I can redirect to the respective state's home page based on ID/name.


Is it possible to create map packages for ArcMap version 10 using 10.1?

Is it possible to create a map package in an older version? Currently I am using 10.1 and the client uses 10.0. Can you create a map package for different versions? I see you can create layer packages using Python for different versions but you cannot specify the version when creating a map package. As I don't have 10.0 I can't test whatever I create these. When you unpack the file using ArcCatalog does it create it v.10. and 10.1? Can you create the original file in 10.0?


Draw lines between two points in QGIS

I have an excel (Fig) which has both X & Y in projected coordinates (WGS84/UTM). Is it possible to draw lines between the X & Y in QGIS ? I've tried to create points based on **X,From_y** and **X,To_y**. I've created two points by just loading them using add delimiter layer (pasted To_Y below From_Y and loaded it). Then tried with Points 2 paths. But it didn't work out. Anyone know how to fix this ? All i need is traverse line between the points ![enter image description here](http://i.stack.imgur.com/hWmhd.png) Edit: This is what i am looking for. I've coordinates(X1, a1 and X1, a2) for n number of points which has to be connected by a vertical line ![enter image description here](http://i.stack.imgur.com/o7gRz.png)


Parse Field to Select Records

I wish to select records based on the NRN attribute in my dataset : `"NRN" IN ('1','2','23','31')` Some records have multiple values in the NRN field separated by the underscore character.   For example : `20_21, 14_23` etc How do I select any records that contain the values that I want ignoring the records that have part of the string but not the whole number? Example : `select 14_23 but not 20_21`.


How to select two points on a map using QGIS?

I need to select two points on a map in QGIS. I selected the element selection tool in menu view and I want to select two points. How can I accomplish this? Are there other tools or methods to select points in QGIS?


How to create one shapefile from multiple rasters?

Ok, I have about ten maps that connect along their edges to each other linking rivers and properties. When I link them all in Inkscape I end up with a large file that won't save (and a smaller one with 5 of the properties saves but is too big and slow in QGIS). So I am thinking that I will load them, about three linked together at a time, into QGIS and create my shapefiles. Is it possible in QGIS to then merge/link the shapefiles along their edges so that the rivers and properties line up and create one large shapefile? additionally, is it possible to then take this large shapefile and georeference it to a base map or can that only be done with rasters? thanks


On the usage of Proj4Js->tmerc.js

I am using "proj4js & open layers" for my mapping. I got projection & datum (which is based on "Transverse Mercator") of a local plane and this plane is part of a grid which returns coordinates of local points in the form of northings and eastings. Now, as we know that Open Layers accepts lattitude and longitude hence I am doing transformation of northings,eastings of a point to lat,long using transverse mercator (tmerc.js) script of proj4js...So my question is will this be enough for the right transformation or there is any thought or suggestion you would like to give regarding the transformation script I am using ? Thanks


Merge nearby items into a single representation

I'm investigating OpenLayers and I'm trying to see if it's possible to filter down the number of displayed feature items. Is it possible to display many nearby items as a single point item on a map when zoomed out. So instead of hundreds of nearly overlapping points I can filter them down a few dozen consolidated points? Or is this something that has to be in the Spatial server?


How do I set custom spatial reference/coordinate sys to my shapefile

I have a python code which converts JSON text file to a shapefile. I have to set coordinate system to this shapefile created. I do not want to access the .prj file from my computer as this code will be embedded in a web toolbox.


Google Maps - Use fitBounds on KML polygon fusion table layer

I have seen example here how to use map.fitBounds on fusiontablelayer here and here , but I think it works with markers table only. Its impossible to parse the lat and long when the vertex coordinates are in the geometry column as in a kml polygon layer. Is there a solution to use fitBounds or other function to zoom to the extent of a kml polygon fusion table layer?


ArcMap: Calculate Quarter using Date field

I am using ArcMap 10.0 and making a map where I symbolize on annual quarter (i.e. Q1 2013, Q2 2013, and so on by quarter through 2015). If a date field contains this: 01/15/2013, and I wanted to create a new field called QUARTER, that would contain the quarter in which the date falls, how could I do that? Right now I'm just definition querying a range of dates, but surely there is a way to just calculate Quarter, right? Thanks in advance!


How do I use st_dwithin with meters?

I am trying to follow PostGIS ST_Buffer Radius Help to use a 50 meter distance but postgis is returning an error. I am trying the following SQL,               SELECT id FROM table      WHERE ST_DWithin(ST_GeogFromText('SRID=4326;POINT(-3.165356 55.926665)'),            geography(geometry), 50);      The hints indicated at ST_GeogFromText are,   * function st_geogfromtext(unknown) does not exist at character 92   * No function matches the given name and argument types. You might need to add explicit type casts. Any help/alternatives would be appreciated.


Natural neighborhoods terminology

Sometimes the hardest part of an analysis is knowing what something is called. What are some R packages, but more importantly, terminology I should be searching for to define contiguous city neighborhoods based on home price changes - basically a topographic map treating these price changes as elevations (and ignoring how real estate people have defined neighborhoods). I'm not sure if this is called "natural neighborhood delination" or "spatially restricted clustering" or just "grouping analysis".


How to change data type in Excel file before creating feature class?

I am trying to import data from Excel spreadsheets where each point coordinate has a 'datetime' timestamp (e.g. "2012-06-01 14:45:32") associated to it. When I check the table properties ArcGIS 10 always recognizes the timestamp field as 'date' (i.e. "2012-06-01") and drops the time-of-the-day info. Is there a way to import the timestamp field as a full datetime type? And if yes, how do I go about it?


Oceans Baseline Data

I am working on a volunteer project and we are in desperate need of a baseline layer. I have been unable to find it anywhere (there are plenty of shoreline and coastline layers out there) does anyone have a source for this data? Thanks, Dave   **Added:**   The baseline from which the territorial sea is measured is mainly the low- water line along the coast as marked on large-scale charts officially recognized by the coastal state.


VB.Net Field Population Methods

Disclaimer: I am basically brand new to VB.NET and ArcObjects. I am computing some statistics as a component of a larger arcobject script. I write these components to an ITable. This ITable is added to the map display as a standalone table. I find that I am duplicating a fair bit of code to assign a value to each field in the table. This is not something I would do in Python and wonder if a better methods applies in the VB.Net world. VB.NET CODE:                  'Populate the stats table with the iteration.         Dim row As IRow = tTable.CreateRow         Dim fields As IFields = row.Fields              Dim fThres As Double = fields.FindField("ThreshDist")         row.Value(fThres) = lStats(0)         row.Store()              Dim fCountR As Double = fields.FindField("CountR")         row.Value(fCountR) = lStats(1)         row.Store()      where lStats is a list of doubles, and tTable is an ITable object. In python I would simply use a dictionary to map the index (as key) to the fields name, as value. A simple for loop would take care of population using an insert cursor. Does a similar paradigm exist using VB.Net? Am I missing something obvious in the documentation?