title: What causes AttributeError: 'module' object has no attribute 'da'?

text: I have a shapefile "selected_polygons.shp" that contains 2000 polygons. I would like to select each polygon contained in the shapefile by using a search cursor in order to apply the function "CostDistance" for each polygon. Here is my code:               import arcpy     from arcpy import env     from arcpy.sa import *          arcpy.CheckOutExtension("Spatial")     arcpy.env.overwriteOutput = True     arcpy.env.workspace = "C:\Users\LCP"     arcpy.env.extent = "costs.tif"          with arcpy.da.SearchCursor("selected_polygons.shp",("*")) as rowsInPatches_start:          for rowStart in rowsInPatches_start:                   ID_patch_start = rowStart.getValue("GRIDCODE")            ## Define SQL expression for the fonction Select Layer By Attribute      expressionForSelectInPatches_start = "GRIDCODE=%s" % (ID_patch_start)            # Process: Select Layer By Attribute in Patches_start      arcpy.MakeFeatureLayer_management("selected_patches.shp", "Selected_patch_start", expressionForSelectInPatches_start)           # Process: Cost Distance      outCostDist=CostDistance("Selected_patch_start", "costs.tif", "", "outCostLink.tif")           # Save the output       outCostDist.save("outCostDist.tif")       The problem is that I obtain a error message when I use:               with arcpy.da.SearchCursor("selected_polygons.shp",("*")) as rowsInPatches_start:      Here is my error message:                AttributeError: 'module' object has no attribute 'da'


title: How to edit line and polygon layers at the same time?

text: I have two layers, lines and polygons. Lines are basically the boundaries of the polygons. Is there any way that when I change or reshape the lines the polygons get updated automatically? I normally change the lines first and since they have to match the associated polygons, I convert them back to polygons again. This way, I need to update the attributes of all the polygons or I have to reshape the polygons according to the change I did to the line layer. When I use topology check I am able to catch the issues and correct them using the correcting topology options. Nevertheless, I was wondering if there is a more efficient way/ approach to do it.


title: Does CartoDB support Retina screens?

text: We are using CartoDB as our production mapping solution. Just wondering does it offers support for **Retina** screens and if so how hard can be to toggle on?


title: How do I convert open street map export to geojson?

text: I am trying to render this output from open street map into geojson to add into a map. I only need to export the data one time. How do I do that? OSM lists it as a use case. I see a js library to do this (but I need a backend solution). It also seems like you can export it as geojson from the api somehow. How do most people usually do this?


title: Warning References to embedded interop assembly in arcobject.net

text: I am including arcobject.net in vb.net project, but i am plagued with following warning. A reference was created to embedded interop assembly 'ESRI.ARCGIS.Geodatabse' because of an indirect reference to that assembly from assembly ESRI.arcgis.carto.Consider changing the 'Embed Interop Types' property on either assembly. Kindly explain the reason behind this warning and how to solve it???


title: How to edit a QGIS plugin GUI?

text: I have inherited a QGIS system from a former employee and my predecessor built a custom plugin for it. I now need to remove a single checkbox option completely from the UI for this plugin. Is there a simple way to do this? Thanks!


title: Radius of curvature for curves in road segments for entire road network

text: I have a road network shapefile for a US state and I am trying to derive the radius of curvature for road segments in this shapefile. I need to be able to perform this task within just a few functions for thousands of road segments. The road segments vary between 0.2 miles and maybe 5 or so miles. Each segment has several curves, but I am only interested in curves small enough to force a vehicle to make a reasonable turn. This should reduce the number of curves of interest per segment to a max of 1 or 2, with many segments not having such curves at all. I have not determined the maximum radius of the curve (the threshold) I would like to use. For segments with only 1 such curve, I would like to obtain the values of those curves for those segments, and for those with more than 1, an average will do.


title: Converting ArcGIS style (symbology) to font for MapInfo/OpenGeo?

text: This is related to Retain ESRI symbology in Mapinfo TAB As a workaround I am wondering if I can save all my ESRI symbology from the lyr files created using the script below to an ESRI style sheet and then convert that to a font set which can be opened in mapinfo/opengeo? ![enter image description here](http://i.stack.imgur.com/N5yOg.jpg) The image above shows some of the [complex] symbology that we have in our ESRI mxd's that need to go into Mapinfo and/or Opengeo to be available via WMS. Tests show that KML created from the lyr files handle the symbology better than tab but not every time. Code so far...               # Author:               import arcpy, string, datetime, shutil, os     import arcpy.mapping as MAP          #Read input parameters from script tool     MXDList = string.split(arcpy.GetParameterAsText(0), ";")     ProjectPath = arcpy.GetParameterAsText(1)          #MXDList=r'P:\2012\183_TownPlanning_Symbology\Working\test.mxd'     #ProjectPath=r'P:\2012\183_TownPlanning_Symbology\Working'          count=0     #Loop through each MXD and print     #for MXDPath in MXDList:     #MXDFile=r'P:\2012\183_TownPlanning_Symbology\Working\Strategic_Plan_Ravenshoe_8.mxd'     count=count+1     print count     ##try:     mxd=arcpy.mapping.MapDocument(MXDFile)     outPath = ProjectPath+'\\'+mxd.filePath[mxd.filePath.rfind('\\')+1:mxd.filePath.rfind('.')]     os.mkdir(outPath)     msd = outPath+'.msd'     dfList = arcpy.mapping.ListDataFrames(mxd)     for df in dfList:     #arcpy.mapping.ConvertToMSD(mxd, msd, df, "NORMAL", "NORMAL")     #arcpy.AddMessage(str(count)+'...' +str(outPath)+'...'+'\n')         lyrList=arcpy.mapping.ListLayers(mxd, "", df)         outPath = ProjectPath+'\\'+mxd.filePath[mxd.filePath.rfind('\\')+1:mxd.filePath.rfind('.')]+'\\'+str(df.name)         os.mkdir(outPath)         for lyrFile in lyrList:             if lyrFile.isFeatureLayer == True:                 if lyrFile.visible == True:                     arcpy.AddMessage(str(lyrFile)+' is visible')                     print str(lyrFile)+' is visible...exporting dataset'                     outFileN=str(arcpy.ValidateTableName(lyrFile.longName[lyrFile.longName.rfind('\\')+1:]))                     try:                         arcpy.FeatureClassToFeatureClass_conversion(lyrFile.dataSource, outPath, outFileN)                         arcpy.SaveToLayerFile_management(lyrFile,outPath+'\\'+outFileN+'.lyr', "ABSOLUTE")                         arcpy.LayerToKML_conversion(outPath+'\\'+outFileN+'.lyr', outPath+'\\'+outFileN+'.kmz')                     except:                         print 'Exists...continuing'                       arcpy.AddMessage(str(lyrFile)+' is visible')             print str(lyrFile)+' is not visible'     ##except:     ##    errorM=arcpy.GetMessages()     ##    arcpy.AddMessage(str(count)+str(errorM)+'\n continuing...\n')     ##    print errorM     ##    #MAP.PrintMap(MXD, printer)     #Remove variable reference to file     del mxd, msd      BTW// I tried the msd route (change to .zip and open the xml) but it is very convoluted.


title: How to cluster kml features in OpenLayers

text: Is it possible to cluster KML features in OpenLayers, and are there any samples? I can find separate samples for clusters and KML but nothing which handles both together. Thanks


title: Why isn't valid geometry not valid in NetTopologySuite when close to poles?

text: I have a Polygon which crosses the south pole: # Polygon in WKT:               POLYGON((5.86826073458095 75.6638564054071,-16.5436836399021 63.2220742891438,-19.3664667294083 59.846221066959,-21.7009816876343 56.4249469806136,-23.6790769742369 52.9700261291777,-25.3903613659649 49.4893115487179,-26.8975641104604 45.9882370271908,-28.2457890135821 42.4706937235238,-29.4682588331783 38.9395620194517,-30.5899861803996 35.3970453090933,-31.630184235186 31.8448854932665,-32.6038903163747 28.2845058942699,-33.5230853693438 24.7171076354231,-34.3974830618665 21.1437357309945,-35.235097920163 17.5653249360032,-36.0426631616959 13.9827314320988,-36.8259448200183 10.3967548983886,-37.5899837597522 6.80815365993192,-38.3392875692471 3.21765455357868,-39.0779879379052 -0.374040045499291,-39.8099755194243 -3.96624846033375,-40.5390214169723 -7.55830680566397,-41.2688934869779 -11.1495673752539,-42.0034749866572 -14.7393968944145,-42.746893755425 -18.3271748440658,-43.5036712696367 -21.9122899273912,-44.278903931454 -25.4941353358308,-45.0784929992426 -29.0721003174899,-45.909447175354 -32.645558114512,-46.780293059812 -36.213846909671,-47.701647966567 -39.7762416028947,-48.6870411529217 -43.3319115133841,-49.7541238657166 -46.8798560026433,-50.9265051564415 -50.4188052017007,-52.2366277028469 -53.9470627335416,-53.7304380696443 -57.4622479660316,-55.4752902472004 -60.960856256507,-57.5739735146974 -64.4374691598322,-60.1910293117637 -67.8832485002322,-63.6053872681201 -71.2828461995608,-68.3234045243048 -74.6074769312702,-75.3289981196023 -77.7943686089884,-86.6900534465662 -80.7085557183446,-106.466508019701 -83.0072017934104,-137.471660384609 -83.9399343949619,-168.142998855139 -82.93831462336,172.499462776297 -80.6056412046207,51.3335176218797 -89.5237677724837,-38.5840495672926 -86.5116366536384,-42.0040683719443 -83.0274485267546,-42.8274013155417 -79.5453266210156,-43.0035444063044 -76.0523116952844,-42.917334253526 -72.5396798562258,-42.697949818896 -69.0250062191839,-42.4008670872053 -65.5077449057289,-42.0533899183215 -61.9874560043951,-41.670338636158 -58.4637908672109,-41.2602920462013 -54.9364821865153,-40.8284114298265 -51.4053395233125,-40.3778451376451 -47.870247975915,-39.910479310127 -44.3311678973098,-39.4273617133736 -40.7881352756541,-38.9289507663875 -37.2412619762738,-38.4152654790081 -33.6907357058719,-37.8859761502671 -30.1368196717061,-37.3404576571625 -26.5798516624587,-36.7778177194076 -23.0202428850507,-36.1969072178556 -19.4584768387702,-35.5963164344505 -15.8951075671059,-34.9743593186678 -12.3307588680448,-34.3290462758966 -8.76612330851596,-33.6580452149769 -5.20196287319082,-32.9586293654512 -1.63911024976531,-32.2276096954905 1.92152745183098,-31.4612484461671 5.47896267724487,-30.6551491136678 9.03211742916007,-29.8041161495169 12.5798092697371,-28.9019756079644 16.120730246266,-27.9413441489322 19.6534181565979,-26.9133294556999 23.1762163144254,-25.8071382075058 26.6872181913455,-24.6095580924228 30.1841914502993,-23.3042664371652 33.6644728848271,-21.8708974452426 37.1248226005016,-20.2837701533619 40.5612194847865,-18.5101356399938 43.9685706157486,-16.50773911296 47.3402932700171,-14.2214065316975 50.6677038885502,-11.5782579077352 53.9391116046634,-8.48105249750914 57.1384547629944,-4.79920779006745 60.2432297260212,-0.357592488403802 63.2213436295558,5.07464609178979 66.0264264522405,5.86826073458095 75.6638564054071))      # Problem: When I view it on google earth, it looks and acts fine. But when I try to do an intersecttion with NetTopologySuite (a C# library) it gives me an error basically saying it's an invalid polygon.               geom = (Geometry)swathGeom.Intersection((IGeometry)_WESTERN_HEMISPHERE);      # Error: > "found non-noded intersection between LINESTRING(-86.6900534465662 > -80.7085557183446, -106.466508019701 -83.0072017934104) and > LINESTRING(-168.142998855139 -82.93831462336, 0 -81.7868941512966) [ > (-101.975528861504, -82.4852087548192, NaN) ]" # Question: How can I get NetTopologySuite to see this as a valid polygon so I can perform an intersection? Do I need to do something with projection? # Polygon in Google Earth: ![enter image description here](http://i.stack.imgur.com/n2OLI.png) ![enter image description here](http://i.stack.imgur.com/MKHeo.png) # Polygon in OpenLayers: ![enter image description here](http://i.stack.imgur.com/sRMkE.png)


title: Why can't I download the georeferencing plugin for QGIS 1.7.4?

text: I used the georeferencing plugin with QGIS 1.6 and now that I've upgraded to 1.7.4 I'd like to do it again. But when I search for the georeferencer plugin, it's not available! Am I doing something wrong?... I also seems that there's no georeferencing option included with version 1.7.4.


title: offsetting symbology position in arcgis 10

text: My data of polylines has lots of attributes to be rendered with a symbology, I made it such that each attribute of interest would be a layer. I would like to display each of these layers at the same time so that the attributes could be seen at the same time. Since the top layers cover the lower ones, I set the symbology so that the layers would have increasing different thickness, with the most bottom layer having the thickest line. Is there a better way to do this, like offsetting a layer's symbology a little bit to the left or right so that each layer can be displayed side by side? Visually, something like this. I'm using ArcGIS 10. Thanks!


title: QGIS toolbox not functional on OSX

text: Has anybody had any luck getting the tools in the QGIS2 Toolbox to run? So far I have tried SAGA and TauDEM but both throw up errors. I eventually got TauDEM to produce something from the D8 tool by feeding it empty rasters to fill but the results were complete nonsense. Outputting to temporary files did not help at all. It's probably yet another path issue, GIS on OSX seems to be one long series of path issues.


title: Why are there no line breaks in my labels even though I added \n?

text: I'm having an issue with creating multi-line labels in QGIS. I believe I have the syntax right in the label expression because the output preview looks correct. However, on the map, everything displays on the same line. If anyone can help, let me know - I'm new to both QGIS and StackExchange. Screenshots are below. ![Label expression and output preview](http://i.stack.imgur.com/enxNS.png) ![Map rendering](http://i.stack.imgur.com/ld2lP.jpg)


title: Reproject a vector layer (EPSG:27700) using WFS onto a googlemaps background (EPSG:900913)

text: I have the following HTML-javascript code and i have problems reprojecting the vector layer (WFS) from EPSG:27700 (British National Grid) to EPSG:900913 (googlemaps). I managed to project the WMS layer but i have difficulties with the WFS. With what i have read so far, to perform this reprojection you need "GeoJSON" and "Proj4js". There is no much documentation and examples using these 2 features. One more thing to add is that using the existing WFS coding i managed to project the "tasmanian" layers that already exist in Geoserver. But these are on a different projecting system (EPSG:4326) Any help would be much appreciated. George Code is the following: html,body { height: 99%; width: 99%; } #map { width: 80%; height: 80%; border: 1px solid black; }               <title>OpenLayers: Google Layer Example</title>     <link rel="stylesheet" href="http://openlayers.org/api/theme/default/style.css" type="text/css" />     <link rel="stylesheet" href="http://openlayers.org/dev/examples/style.css" type="text/css" />          <script src="http://openlayers.org/api/OpenLayers.js"></script>         <script src="http://www.openstreetmap.org/openlayers/OpenStreetMap.js"></script>     <script src="http://maps.google.com/maps/api/js?v=3.2&sensor=false"></script>         <script src="http://ecn.dev.virtualearth.net/mapcontrol/mapcontrol.ashx?v=6.2 mkt=en-us"></script>         <script src="http://developers.cloudmade.com/attachments/download/58/cloudmade.js"></script>          <script src="http://proj4js.org/lib/proj4js-compressed.js"></script>          <script type="text/javascript">              Proj4js.defs["EPSG: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 +no_defs";         Proj4js.defs["EPSG:27700"] = "+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs";         Proj4js.defs["EPSG:4326"] = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs";            function init() {              //set up projections              // World Geodetic System 1984 projection                 var WGS84 = new OpenLayers.Projection("EPSG:4326");                      // WGS84 Google Mercator projection         var WGS84_google_mercator = new OpenLayers.Projection("EPSG:900913");              // OS British National Grid         var os = new OpenLayers.Projection("EPSG:27700");                      //Initialize the map                 var map = new OpenLayers.Map ("map", {                         controls:[                                 //allows the user pan ability                                 new OpenLayers.Control.Navigation(),                                 //displays the pan/zoom tools                                 new OpenLayers.Control.PanZoomBar(),                                 //displays a layer switcher                                 new OpenLayers.Control.LayerSwitcher(),                              ],                                 projection: WGS84_google_mercator,                             displayProjection: WGS84                 } );                      //base layers                 var openstreetmap = new OpenLayers.Layer.OSM();                      var cloudmade = new OpenLayers.Layer.CloudMade("CloudMade", {                     key: 'BC9A493B41014CAABB98F0471D759707'}                 );                      var google_maps = new OpenLayers.Layer.Google("Google Maps", {             numZoomLevels: 20}         );                 var google_satellite = new OpenLayers.Layer.Google("Google Satellite",                          {                         type: google.maps.MapTypeId.SATELLITE,                         'sphericalMercator': true,                         buffer:0,                         numZoomLevels: 20}                 );                 var bing_shaded = new OpenLayers.Layer.VirtualEarth("Bing Maps",                  {                          type: VEMapStyle.Shaded,                         'sphericalMercator': true,                          numZoomLevels: 20}                 );         var bing_hybrid = new OpenLayers.Layer.VirtualEarth("Bing Satellite", {             type: VEMapStyle.Hybrid,                    'sphericalMercator': true,             numZoomLevels: 20}             );         var bing_aerial = new OpenLayers.Layer.VirtualEarth("Bing Aerial", {             type: VEMapStyle.Aerial,                         'sphericalMercator': true,             numZoomLevels: 20}             );                      // WMS Layers for the projection of the shapefiles                 var layer1 = new OpenLayers.Layer.WMS( "National Parks",                             "http://localhost:8080/geoserver/wms",                           {                             srs: 'EPSG:4326',                                 layers: 'National_Parks',                                 transparent: 'true',                                 format: 'image/png',                                 buffer:0                           },                           {                                 'singleTile': true, 'ratio': 1, 'opacity': 0.75, 'isBaseLayer': false, 'wrapDateLine': false,'transitionEffect':'resize', displayInLayerSwitcher : true, 'alwaysInRange': true, 'displayOutsideMaxExtent':true                           }                           );                      // WFS Layers for the projection of the shapefiles                 var wfslayer = new OpenLayers.Layer.Vector("National Parks in GB",                 {                 strategies: [new OpenLayers.Strategy.BBOX()],                 protocol: new OpenLayers.Protocol.WFS({               url:  "http://localhost:8080/geoserver/wfs",                           //featureType: "tasmania_state_boundaries",                           featureType: "National_Parks",                           //featureNS: "http://www.openplans.org/topp",                           featureNS: "http://opengeo.org/protec_are",                           version: "1.1.0",               srsName: "EPSG:900913",               })                         });                                      //transformation between coordinate systems (27700 to 900913)                                 var format = new OpenLayers.Format.GeoJSON({                                 'internalProjection': new OpenLayers.Projection("EPSG:900913"),                                 'externalProjection': new OpenLayers.Projection("EPSG:27700")                                 });                                           var newfeatures = format.write(wfslayer.features);                                 var newlayer = new OpenLayers.Layer.Vector("newlayer");                             newlayer.addFeatures(format.read(newfeatures));                           // Display the layers on the map                 map.addLayers([openstreetmap, cloudmade, google_maps, google_satellite,bing_shaded, bing_hybrid, bing_aerial, layer1, newlayer]);                      var mapextent = new OpenLayers.Bounds(-16.08,49.06,6.76,58.74).transform(map.displayProjection, map.baseLayer.projection);         map.zoomToExtent(mapextent);              }          </script>      OpenLayers: Google Layer Example * * * Thanks a lot for the information guys. Senthil i made the changes that you suggested and i managed to get the layer on the screen. As quite lots of people might have the same problem i would like to summarise some important things. --> If the web page(HTML) containing the map (comment from 'iant') that you are creating exists on the same server as WFS/Geoserver (e.g when you use the same PC) then you dont need to use proxy. --> For conversion between projection systems (e.g. from EPSG:27700 to EPSG:900913) you dont need Proj4Js or GeoJson. WFS Version 1.1.0 will do the trick. --> Do not run the HTML web page as a file by double clicking on it (e.g file:///E:/GISWork/GISApplication/Geoserver-2.0.3/data_dir/www/GIS.html) but run it through the server (e.g `http://localhost:8080/geoserver/www/GIS.html`) otherwise you wont be able to see the layer on the screen. G


title: How to activate Upload Layers and other functions in Geoexplorer using OpenGeoSuite 4.x

text: We installed recently the Community Edition of Opengeosuite 4 on our system. Compared to the OpenGeoSuite 3 (which run perfectly) we encounter many problems. We need to be able to edit and upload shapefiles via the GeoExplorer. Unfortunately in in OpenGeoSuite 4.x there is no Upload Layer button any more. When we try to edit a shapefile that we added to Geoexplorer (via the geoserver admin interface) the edit fails with the notion that the Layer is Read-Only (even though we have set the right permissions in Geoserver and Ubuntu). How can we enable the functionality of uploading layers and edit shapefiles via Geoexplorer in the OpenGeosuite 4.x version.


title: How to post an existing shapefile datastore using cURL?

text: I read geoserver doc, and already Uploaded a shapefile to geoserver using cURL. But I can't change layer name by REST. So I think that there is another way to post a layer. If I can post an existing shapefile datastore, this problem will be solved. and is it possible change layer name(not title) by cURL? This is my command line:               F:\\new\\curl-7.22.0\\src\\DLL-Release\\curl.exe -v -u admin:geoserver -XPUT -H "Content-type:text/xml" -d "file:///data/shapefiles/states.shp" http://********:8081/geoserver/rest/workspaces/myworkspace/datastores/16-43/external.shp


title: Location-based game using GPS

text: I am not sure if this forum is proper place for that kind of question but i hope to get some advice... I looking for some scenario (examples, ideas) for urban game (location-based game) with usage of GPS devices. The thing is that the game should touch on spatial planning subject. It will be for teenager audience.


title: Conversion huge raster (50 000 000 cells) in to .xyz

text: I do have DEM approx. 40x5km with cell size 2x2m, which gives approx. 50 000 000 points. I need to export it to .xyz (1pt = 1 line in format x y z) file in order to use it in MIKE21. It seemed like there is loads of tools that will do it (I do have access to Arcinfo 9.3.1 and 10), though I'm stuck for 2 days already. I tried some solutions and all have failed. I did: 1. SA/EXTRACTION/SAMPLE, receiving table, saving it as .dbf, and trying to import it in to Access 2007. I could not import .dbf produced in Arc into Access. Exporting table to .txt gives me always error.   1. Raster Calculator and Sample(). This worked great. Raster converted in 5 minutes producing 1.2GB file, but in order: ZXY not XYZ. MIKE do not allow to change order of columns, so I need to change order of columns. No clue how to do it. Tried Access 2007 again, but failed on import data.   2. Script Raster2XYZ ver. 1.1 from ESRI website (in VB). Can not install it properly now (gives an error), I used it before and it worked 6 months ago. However I remember that it worked really slow back then. I am not much of a programmer and can not find the bug.   3. Grass GIS and r.out.xyz have bug. Problem with double or float writing(can't remember now).   4. I have R script which does export .asc to .xyz. It take 2.5h to export 40 000 000 cells. Output file is 4.3GB. MIKE will not swallow it. Any ideas how to export such raster to .xyz? Or how to reorder file from pt.2 to receive XYZ instead ZXY? I always used Access to convert data (ie. .dbf to .txt(.xyz)), but now it seems like no option. Is there any other solution? I do have HP workstation Z600 (Xeon 2.4Ghz, 8GB RAM, GF Quadro FX 3800). Any sugestions?


title: Calculate distances between series of points in postgis

text: I have a postgres/postgis table containing a series of sequential points. I want to add a new column to the table which contains the distance between each point and the first one, and so forth until all the points are processed (the first point should have a distance of 0). Can this be done with a query or do I need to use Python or equivalent to create a loop? Many thanks


title: SOAP url of arcgis service

text: I basically use REST URL to display feature layer or ArcGISDynamic Service layer and it is commonly promoted by ESRI by their api example. I was wondering since the ArcGIS server also provide SOAP URL along with the REST URL so is it possible to use SOAP URL to configure feature layer or ArcGISDynamic Service using javascript. So I can put SOAP URL to good use also SOAP protocols obviously has some good advantages.


title: arcGIS ADF - map disappear when I add graphic layers

text: I'm trying to add a graphic layer in a web ADF application with arcGIS for asp.net, but every time I add a graphic layer to the MapResourceManager, the map disappears! To better understand this: here's the simple app with the only element in the MapResourceManager being the map service (OGC (WMS) Service) which provides me the map http://imageshack.us/f/715/pro2w.jpg/ And here's the MapResourceManager with another graphic layer, the map misteriously disappears http://imageshack.us/f/543/pro1y.jpg/ What's wrong with that?


title: Looking for a Provider of UK boundary data

text: I'm working on a web project centered around the input and mapping of GIS data. We have a need to display polygons of town boundaries. The RightMove.co.uk website has and uses this data: http://www.rightmove.co.uk/property-for- sale/map.html?locationIdentifier=REGION^20140&insId=1 ![RightMove overlaying town boundary on Google map](http://i.stack.imgur.com/I7jix.png) **Does anyone know a provider of this information (town/city boundaries)?** I've contacted Google sales, but they've said they don't have this data. So it appears RightMove got it from another source.


title: Saving Google Earth Location: OpenLayers Plugin in Quantum GIS

text: I use the OpenLayers plugin to load Google Earth Layers in Quantum GIS to digitize/identify certain features. As my work is large occassionally I have to quit Quantum GIS to resume later. The problem I face is though I save my project as a qgs file I am not able to get the dynamically updated layer of Google Earth in the next session. The Google Earth layer just appears in the top left corner of the map area in QGIS and it shows only a portion of the area I work with. I can't zoom/pan further to other places. Also I am not in favour of using Google Earth plugin for this work as it just takes a snapshot of the scene I have zoomed onto in Google Earth.


title: Determining accessibility of U.S. street address

text: I'd like to determine whether a given U.S. street address is publicly accessible to vehicle/foot traffic. This would include addresses on both publicly and privately maintained roads that are not within a gated community or another access-restricted area. For instance, 8541 Lake Windham Ave, Orlando, FL 32829 is in a gated community. The openstreetmap tag for this portion of the road is tagged as access=private. However, since OSM data does not include address ranges for each segment (AFAIK), I'm not sure how to associate the address with the private portion of the road. The U.S. Census Tiger data does include address ranges. However, as far as I can tell from the documentation, there's nothing to indicate accessibility of the segment. Is there a freely-available way to determine whether a U.S. street address is on a publicly accessible road?


title: Cannot log on to ArcIMS Web Manager, first time logon

text: My current problem with arcims web manager dot net in arcims 9.3: after installing arcims and configuring it successfully and then installing web manager dot net, I cannot log on to web manager. Does anyone have a solution to this problem? Screenshot More information:   * OS: Windows xp sp3.   * Arcims 9.3, with java jre5 u13, ServletExec5.0ISAPI, IIS6   * arcims web adf dot net 9.3, integrated VB.Net 2008. Thank you.


title: Is it possible to access ArcGIS's table operations "Summarize..." and "Statistics..." in python?

text: ArcGIS has a few nifty table tools, including Summarize and Statistics. Is there a way to access them through arcpy, or will I have to hack my own methods via math and lists?


title: Where can I find Poverty shapefile?

text: Does anyone know where I can download shapefile file that has poverty dataset for California?


title: Switch line direction for PostGIS-lines in QGIS or Pgadmin

text: I am trying to switch the line direction in a PostGIS-layer. The postgis- geometry column is e.g. "LINESTRING(0 0, 1 1, 2 1, 2 2)". The reason is that my symbology is direction-dependent, but I drew some lines in the wrong direction. I have tried using my usual approach that seems to only work for shapefiles. It is the code by Nathan W. that is in the answer that is linked below. Link to similar question:   * How can I switch line direction in QGIS?


title: Chart Production - chart series, etc

text: We do subsea surveys - seabed mapping etc. What software do you use to produce charts? I am especially interested in producing chart series. It would be interesting to hear how you guys produce charts... Best regards, W


title: Can I tell QGIS to load/show only a part of an image?

text: can I tell QGIS to load/show only a part of an image without to extract this region as an extra/ new image?


title: Merging DEMs of different grid size

text: I'm trying to merge two grids with different cell size. One has a cell size of 10m and the other 50m. The 10m one covers a subarea of the 50m one. I know that the most accurate one is the 10m one and hence I would like to keep its values in the overlapping areas, and not do any kind of interpolation. My problem is that the resulting merged grid (after resampling the 50m grid to 10m) has obvious step changes as seen in the picture ![Merged grid](http://i.stack.imgur.com/mneK8.png) Is there a way I can feather out these discrepancies without editing the grids cell by cell? I'm using Qgis and Saga. Thanks


title: return to newline arcgis/arcSDE/postgres labeling

text: i have two tables :   * parcel and owners with many to many relation ship so i want to print owners name in my parcels (type polygon), so i coded using .net a function that create field in parcel table and put a concatenation from all parcel's owner.   * the problem is when i want to add labeling to parcel , all owners placed in one line : ![enter image description here](http://i.stack.imgur.com/8D1Sa.jpg)   * and what i want is to have one line for one owner like this ![enter image description here](http://i.stack.imgur.com/IYlMH.jpg)     * Function                          private string get_names( DataRow row, int id_pers)          { StringBuilder st = new StringBuilder(); st.Append("select nom_pll, nom_fr ,prenom_fr from parcelle pa,personne p,pll_personne pll "); st.Append("where p.objectid = pll.id_pll "); st.Append("AND pa.objectid= pll.id_personne "); st.Append("AND nom_pll= " + row[0].ToString()); st.Append(" GROUP BY nom_fr , prenom_fr, nom_pll");                          DataTable dt = pcnx.Select(st.ToString());         string pers = " ";         if (dt.Rows.Count > 0)         {             pers +=  dt.Rows[0][1] + " " + dt.Rows[0][2] + " , ";             for (int j = 1; j < dt.Rows.Count; j++)             {                 pers +=  dt.Rows[j][1] + " " + dt.Rows[j][2] + " , ";             }             pers +=   dt.Rows[dt.Rows.Count - 1][1] + " " + dt.Rows[dt.Rows.Count - 1][2] ;         }         //resultat ligne = new resultat(row[id_pers].ToString(), pers, row[id_pers + 1].ToString());          return pers; } reference to : Problem with ArcGIS annotation


title: The coordinates measure are out of bounds in Android route task

text: I am a beginner in ArcGIS and I was following this tutorial. I am trying to find a route between two points. I created a sample application, in which first, i get the two points on map single tap as follows               @Override             public boolean onSingleTap(final MotionEvent e)             {                 Point point = map.toMapPoint(new Point(e.getX(), e.getY()));                 points.add(point);                 if (points.size() == 2)                 {                     new FindRouteTask().execute();                 }                 return true;             }      and my FindRouteTask is as follows               private class FindRouteTask extends AsyncTask<Void, Void, RoutingResult>         {                  ProgressDialog  dialog  = new ProgressDialog(FindeRouteActivity.this);             RoutingResult   mResults;                  @Override             protected void onPreExecute()             {                 dialog.setMessage("Finding routes...");                 dialog.show();             }                  @Override             protected RoutingResult doInBackground(Void... params)             {                 RoutingParameters routeParams = new RoutingParameters();                 routeParams.setImpedanceAttributeName("Length");                      String routeTaskURL = "http://tasks.arcgisonline.com/ArcGIS/rest/services/NetworkAnalysis/ESRI_Route_NA/NAServer/Route";                 RoutingTask routeTask = new RoutingTask(routeTaskURL, null);                      NAFeaturesAsFeature naFeatures = new NAFeaturesAsFeature();                 // Create the stop points from point geometry                 StopGraphic startPnt = new StopGraphic(points.get(0));                 StopGraphic endPnt = new StopGraphic(points.get(1));                 // set features on routing feature class                 naFeatures.setFeatures(new Graphic[] { startPnt, endPnt });                 // set stops on routing feature class                 routeParams.setStops(naFeatures);                 try                 {                     mResults = routeTask.solve(routeParams);                 } catch (Exception e)                 {                     e.printStackTrace();                     return null;                 }                 return mResults;             }                  @Override             protected void onPostExecute(RoutingResult result)             {                 dialog.dismiss();                 points.clear();                 if (result == null || !result.getRoutes().isEmpty())                 {                     Toast.makeText(FindeRouteActivity.this, "There was some problem in fetching the routes. Please try again!", Toast.LENGTH_SHORT).show();                 } else                 {                     findRoute(result.getRoutes().get(0));                 }             }              }      But the problem is that `routeTask.solve(routeParams)` always throws a exception               Failed to convert the input propertyset into a recordset to load into the NAClass "Stops".  The coordinates or measures are out of bounds.        Could anyone please tell me what am I doing wrong and what the problem is??   -Thanks in advance


title: WMS Geometry Vector point label Styling

text: How can we add category based styling to geometry point my application show geoserver rendered map, and all the map POI are directly fetched from database and shown to the map.i need know how we can add different style for different points please help me database used is POSTGIS


title: How to handle multiple space separator in Grass?

text: I try to load a Lidar point cloud in .xyz format to Grass (in QGIS 2.0.1). I tried to input the following options to the 'separator' field: `tab`, `'tab'`, `'\t'`. Just for the unlikely case that it worked I tried alse `' '`, `'space'`, `'newline'`, `'comma'` and `';'`. Two things have changed. First: empty space between two apostrophes didn't work yesterday, but it works now (I don't know why). Phrasing it as `'space'` still doesn't work. Second: I found the data are actually separated by two or three spaces, not by a tab. When I set space as the separator, there is not the old wrong separator error, but Y coordinate seems wrong: it looks for it between the two spaces separating the fields. Writing `' '` (double space) to the separator field doesn't help. So I rephrased this question to fit my problem better: **how to handle multiple space separator in Grass?** First few lines of the .xyz file look like this:                 3453939.450  5633900.030   361.670       3453941.750  5633900.020   361.720      Other options are the default ones, just the "statistic" is "median" instead of "n".


title: Find Distance within each Raster contiguous cell block

text: I was asked to find for each Individual Contiguous Forest Block the distances INSIDE of each where it would be the farthest (aka the most isolated) from each forest cell block's respective edges. I know this should involve Euclidean Distance and possibly followed by Reclassify, but perhaps not? I tried performing EucDist on the forest cell block layer with a mask of that same forest cell block layer, but the result didn't seem to be what I wanted. ![Image of individual forest cell blocks that are all > 1000cells](http://i.stack.imgur.com/rH8PD.jpg)


title: Grouping points by location

text: This is ArcGIS specific. I do have surveyed points (sections). There will be 200 sections, each section 20-30 points. I want to create polylines (sections) out of the points. In order to do so, they have to be in right sequences (ordered from right to left or so) so program (ie. ETWizzard) would know how to connect the dots. Points in sections do have random IDs. I want to write script, where logic is: assign points from .shp for each section to specific group (ie. S1, S2...). Then try to give names (sort) to points in section based on x,y. Preferably, I want to do this in Python. I guess the most important question at this stage is: Is there a tool to select and group points based on the location/proximity or any one would hear about any similar solution? There has to be one, since sections are quite commonly used. Any thoughts?


title: OpenLayer - center map based on lon/lat of POI

text: Is it possible to center the map to the coordinates of a POI, imported by using a text file? (OpenLayers.Layer.Text)               var pois = new OpenLayers.Layer.Text( "My Points",                     { location:"./textfile.txt",                       projection: toProjection                     });     map.addLayer(pois);          map.setCenter (????, zoom);


title: In need of WFS-T support in a WPF/Silverlight application

text: I want to build a wpf/silverlight mapping application that enables real-time visualization of features in a map, and that allow different users to add features (points, polygons) to the map. My initial thought is to use WFS-T for this purpose, but unfortunately, I cannot find any free/cheap SDK that provide support for WFS-T within WPF. My question is therefore: Does anyone know about a free/cheap SDK that enable support for WFS-T within WPF applications? p.s. there's no need to mention ArcGIS products as an alternative, I can't afford them. Thank's for your attention.


title: ways to access and download the metadata of OpenStreetMaps

text: Does anyone know if there are any good methods or apps to access and download the metadata of OSM from desktop? Many thanks!


title: How do you programmatically read the coordinate display format in ArcMap?

text: I'm looking for a property (or properties) that will give me the display units and the formatting. ![enter image description here](http://i.stack.imgur.com/zwzoS.png) I assume I need `IMap.DistanceUnits` but I can't find the formatting. Is this exposed? I am also interested in accessing the additional formatting in this window: ![enter image description here](http://i.stack.imgur.com/ESoPm.png) I'm hoping there's a `INumberFormat` property someplace that I can hook into.


title: Comparison of Java-based GIS Programs

text: Hi I want to develop GIS Application in java I found some of useful software like UDIg,Openjump,KOsmo & gvsig but i am unable to decide from which should i start as i want to create my application which i can give to my clients Thanks in Advance


title: GIS data - Geoprocessing Service to receive data

text: I would like to move 8 shapefiles from one intranet to another intranet via arcgisserver (or an other method). Creating a FTP site is not an option. I think I saw somewhere on how to send data to a geoprocessing service. Is this feasible? Sending a shapefile to a geoprocessing service outside of the api with say python? Thanks. My basic need is to move 8 shapefiles from one intranet to another intranet.


title: Is there any way to suppress the calculate field dialog box in arcpy script?

text: I've got a python script that I inherited. The python script takes input from a variety of feature classes and tables and spit it out to a single table. When I run this script via a scheduling utility at night, the script fails when it prompts for user input on whether it's OK to proceed with a truncated field name with a field calculation. For a variety of reasons, I can't change the field name. Is there a way to suppress the input box via code? I looked in arcpy.env but I didn't see anything...any ideas are welcome! thanks!               arcpy.CalculateField_management(Merged__2_, "L1_Parcel_Nbr", "[PropID]", "VB")      where Merged_ _2_ _ is the output (file geodatabase) feature class from a merge geoprocessing function.


title: Problem while adding the store/layer in geoserver

text: I am new to geoserver..Yesterday I worked with that and added stores and layers. Its working fine. But, Today I try to add a new store, i am getting this error message.. _**"Sorry, your session timed out... It looks like you waited too long to make that last change.if this continues to happen, you should get in touch with your system administrator. Go back to the home page and try again."_** How to resolve this error.. Could u please help me out.. Thanks in advance


title: How to best model GPS tracks for storage, visualization and analysis?

text: I am thinking about writing software to deal with GPS Tracks and Waypoints (mostly storing, displaying and calculating metrics such as speed, grade, and some simple statistics). I wonder what should be the most conceptually robust data model regarding trackpoints, and here are some "candidates":   1. Considering Tracks as sequences of Trackpoints: 1.1. Tracks are considered "2D", since map projections are 2D. Trackpoints might or not have elevation, might or not have timestamp. Elevation and timestamp are consideres "extras", "optional". For terrestrial applications, elevation is a direct function of lat/lon (obtainable via DEM); 1.2. Tracks are considered "3D" since geographic space is, indeed 3D, and the trajectory of the receiver is 3D (2D projection is thus a form of data reduction). Timestamp might or not be present (the track could have been drawn by hand). 1.3. Tracks are considered "4D" (3 spatial + time). Thus, a hand-drawn map is a special case where elevation and timestamp are `null` or otherwise not present, but the Trackpoint properties are always "there".   2. Tracks are considered dictionaries of streams, where all the streams have equal lengths. There is a list of latitudes, a list of longitudes, a list of elevations, one of timestamps, etc. This makes easy to calculate statistics of each property, and the concept of Trackpoint becomes "virtual" in a sense, since it is a cross-section of many streams. If I understood right, GPX format adopts 1.1., KML adopts 1.2. (with no support for timestamp), and Strava API adopts 2. (in JSON format), but in the end these are just FILE formats for serialization and storage, not necessarily for modeling, computational representation and numbercrunching. Is there any form that is preferrable, in an object-oriented sense, and why? (I believe that strong typing and sensible modelling at least would avoid operations that doesn't make sense). EDIT: some "intriguing" additional questions:   * Is a hand-drawn track CONCEPTUALLY the same thing that a device-recorded tracklog? Should they be of different data types?   * Should it be considered "correct" that KML stores null elevations as zero? Zero IS an elevation, and if you don't know the elevation you shouldn't assign a numeric zero to it, isn't it?   * Should it matter, in a track with elevation, if the elevation is extracted from DEM data ("offline") or from GPS data or barometric data ("in the field")? Should this be flagged in the Track object? Saved to different Trackpoint properties? Ignored? Should they be different collection datatypes?   * If I edit a device-recorded track in a map editor (adding, moving and removing points), or combine tracks from different dates, how should the timestamps in the trackpoints be handled? Should they be "resetted" to null? Should an object (trackpoint collection) of a different type be created from the former ones?


title: ArcGIS API for Silverlight: Retrieve all Data From a Raster layer?

text: Is it possible to use ArcGIS api for silverlight to get the data in each cell of a raster layer and store them in an array?


title: Survey Data to GIS database

text: This is a General Question on a possible workflow for Survey Data directly to ArcGIS In most GIS/Survey workflows there is a field collection team and backoffice GIS team. When the availability of high-end mapping devices are there, why can we directly edit a Geodatabase in a Survey device or rather create a GIS feature in Survey device itself, eliminating the CAD and GIS feature creation ? What my Question is , is there any Survey Device ( Totalstation etc) currently supporting data editing in GIS format ? Is there any Survey device capable to support a full fledged Geodatabase ? An answer from the community is really a headstart ! Thanks in Advance


title: Point in polygon analysis fails with "TypeError: arguments did not match any overloaded call"

text: I am using QGIS 1.8 on Windows 7 I want to perform a point in polygon analysis using the Vector > Analysis Tools menu. My point layer is the result of a SpatiaLite query which returns 10'000 points. My polygon layer is a simple shapefile with about 150 different polygons. Both layers are in WGS84 (EPSG:4326). When I click the "OK" button in the Point in Polygon tool, just nothing happens. The progress bar stays at 0%, but I cannot click "Cancel" (i.e. nothing happens). I can click on the Close-Button in the upper right though, the window then disappears. **EDIT:** This is the output from the message log:               An error has occured while executing Python code:          Traceback (most recent call last):       File "C:\OSGeo4W\apps\qgis\python\plugins\fTools\tools\doPointsInPolygon.py", line 215, in run         tmpGeom = QgsGeometry(pntFeat.geometry())     TypeError: arguments did not match any overloaded call:       QgsGeometry(): too many arguments       QgsGeometry(QgsGeometry): argument 1 has unexpected type 'NoneType'          Python version:     2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)]               QGIS version:     1.8.0-Lisboa Lisboa, f440b60          Python path: [...]      I have also tried saving the point layer coming from a SpatiaLite DB as a shapefile and reimporting this shapefile. I have also tried to do it with a much smaller subset of points. No success.


title: How to debug third party JavaScript libraries such as openlayers?

text: I am able to work with web development tools like firebug and dev tools and solve bugs in my javascript files. But when it comes to bugs that produce errors such as "typeError: e is null" which prompts me to the third party js library I am using, how am I supposed to debug that issue? I am creating a django app to integrate with a bigger django project. I am using openlayers as simple as possible and even though I am not getting any errors, apart from a css warning, my map division is not loading any maps...How am I supposed to work around this and of course any future similar problems? Any advice or link would be much appreciated. Kind regards, A


title: Mapserver 6.2 concatenate attributes in FILTER

text: Back in version 5.0 I had:               FILTER ("<tcode:expression/>" eq '[STATE][COUNTY]')      final version = FILTER ("11000" eq '[STATE][COUNTY]')      so I would filter out the features with state = 11 and county = 000 Now, in 6.2, it does not work anymore. Does anyone know how to make it work with concatenated attributes? I know that it would work like this:               FILTER ("11" eq '[STATE]' and "000" eq '[COUNTY]')      but it would really help me to find a way for the concatenation to work like:               FILTER ("11000" eq ('[STATE]'.'[COUNTY]')) -no error, but not display anything      Been browsing the documentation for a while but I just can't make this work.


title: ESRI Versioning: A table doesn't match base table schema

text: ArcGIS 10.1 sp1 with Oracle 11g I've got a versioned feature class that I have discovered has a base table with a different schema than the corresponding A table. (A + registration_id from sde.table_registry) The only possible explanation I can think of is someone went in and manually added/removed fields from the base table or the A table. I don't BELIEVE this has happened because there aren't many people with access to this database and I would hope they all know better than to do this on a versioned feature class/table... Are there any other ways that this could happen? The only other item of interest is that this was just upgraded from 9.3.1 to 10.1. I'm in the process of trying to get a 9.3.1 backup used for the upgrade restored to see if that contains this problem... I was initially alerted to this problem while trying to export data. The export to xml workspace fails completely without any info, but the copy feature class via Ctrl+C in ArcCatalog throws the error: [ORA-00904: "A"."DUCT_MATERIAL":invalid identifier][GEOSPATIAL.SPAN] I guess my question should be how did this happen and how do I fix it. I could probably go through and write some sql statements to add these fields, but this may not be the only table with this problem and it would be tedious. Any automated way of doing this?


title: How does ST_Azimuth work?

text: I do not understand the function behaviour. I am giving the function a point B that is North of point A so I would expect something around 90 degrees but I get 263 degrees instead. In another example, I keep the same point A and I am giving the function a point B that is on the South of point A and this time I do get a correct value of 279 degrees. The function definition states that it is using the horizontal of the vector between A and B. But I don't understand the part which says: > Angle is computed clockwise from down-to-up: on the clock: 12=0; 3=PI/2; > 6=PI; 9=3PI/2. Here are the examples, maybe I am doing something wrong with the SQL point b NE of point A:               SELECT degrees(                     ST_Azimuth(                         ST_Transform(ST_GeomFromText('POINT(127.702 38.169)',4326),32652),                         ST_Transform(ST_GeomFromText('POINT(126.728 38.178)',4326),32652)                                )                     )      This returns 271 degrees whilst it should return ~45 degrees Someone care to explain? Thanks in advance!


title: Python: find a method to calculate the “inner centroid” (also known as labelPoint in arcmap) of a polygon

text: i am looking for a method to compute the X,Y coordinate of a inner point (also known as labelPoint in arcmap) of a (irregular) polygon. I am try to find a solution without the use of ArcGIS module ![enter image description here](http://i.stack.imgur.com/uisSD.jpg) the difference between a "centroid" and a labelPoint is that the centroid may be outside the polygon. I am using shapely module from shapely.geometry import Polygon               ref_polygon = Polygon(points)     # get the x and y coordinate of the centroid     ref_polygon.centroid.wkt     'POINT (558768.9293489187300000 6361851.0362532493000000)'      using the following data (the vertex of my polygon) thanks in advance for help and suggestions               points = [(560036.4495758876, 6362071.890493258),               (560036.4495758876, 6362070.890493258),               (560036.9495758876, 6362070.890493258),               (560036.9495758876, 6362070.390493258),               (560037.4495758876, 6362070.390493258),               (560037.4495758876, 6362064.890493258),               (560036.4495758876, 6362064.890493258),               (560036.4495758876, 6362063.390493258),               (560035.4495758876, 6362063.390493258),               (560035.4495758876, 6362062.390493258),               (560034.9495758876, 6362062.390493258),               (560034.9495758876, 6362061.390493258),               (560032.9495758876, 6362061.390493258),               (560032.9495758876, 6362061.890493258),               (560030.4495758876, 6362061.890493258),               (560030.4495758876, 6362061.390493258),               (560029.9495758876, 6362061.390493258),               (560029.9495758876, 6362060.390493258),               (560029.4495758876, 6362060.390493258),               (560029.4495758876, 6362059.890493258),               (560028.9495758876, 6362059.890493258),               (560028.9495758876, 6362059.390493258),               (560028.4495758876, 6362059.390493258),               (560028.4495758876, 6362058.890493258),               (560027.4495758876, 6362058.890493258),               (560027.4495758876, 6362058.390493258),               (560026.9495758876, 6362058.390493258),               (560026.9495758876, 6362057.890493258),               (560025.4495758876, 6362057.890493258),               (560025.4495758876, 6362057.390493258),               (560023.4495758876, 6362057.390493258),               (560023.4495758876, 6362060.390493258),               (560023.9495758876, 6362060.390493258),               (560023.9495758876, 6362061.890493258),               (560024.4495758876, 6362061.890493258),               (560024.4495758876, 6362063.390493258),               (560024.9495758876, 6362063.390493258),               (560024.9495758876, 6362064.390493258),               (560025.4495758876, 6362064.390493258),               (560025.4495758876, 6362065.390493258),               (560025.9495758876, 6362065.390493258),               (560025.9495758876, 6362065.890493258),               (560026.4495758876, 6362065.890493258),               (560026.4495758876, 6362066.890493258),               (560026.9495758876, 6362066.890493258),               (560026.9495758876, 6362068.390493258),               (560027.4495758876, 6362068.390493258),               (560027.4495758876, 6362068.890493258),               (560027.9495758876, 6362068.890493258),               (560027.9495758876, 6362069.390493258),               (560028.4495758876, 6362069.390493258),               (560028.4495758876, 6362069.890493258),               (560033.4495758876, 6362069.890493258),               (560033.4495758876, 6362070.390493258),               (560033.9495758876, 6362070.390493258),               (560033.9495758876, 6362070.890493258),               (560034.4495758876, 6362070.890493258),               (560034.4495758876, 6362071.390493258),               (560034.9495758876, 6362071.390493258),               (560034.9495758876, 6362071.890493258),               (560036.4495758876, 6362071.890493258)]


title: 'Select by attribute' using python in ArcMap 10

text: I need to select a feature using SQL before zooming to it in python. What I'm looking for are replacements for the commented lines in this code:               import arcpy     parcel_num = "123456"     mxd = arcpy.mapping.MapDocument('CURRENT')     df = arcpy.mapping.ListDataFrames(mxd, "Layers") [0]     # df.selectFeatures("SELECT * FROM LandRecords.GISDBA.Parcels "     #                   "WHERE PARCEL_NUM='%s'" % parcel_num)     df.zoomToSelectedFeatures()     arcpy.RefreshActiveView()


title: FieldPyculator Plugin - How to concatenate fields?

text: How can i concatenate two fields using pyculator plugin? There is any document with all pyculator functions that i can use in "field expression"?


title: why is a gdb query WHERE clause invalid according to the ArcGIS REST API?

text: **Update:** This only occurs when the layer being queried has a **join**. It doesn't matter if it has a join to a different database or not. And the query that fails is not trying to access columns in the joined table. In order to more accurately express the problem, I've used fully-qualified (4-part) column names below where required. See the updated WHERE clause below. Original question: We are using ArcGIS Server 9.3.1, Java, and ArcSDE with SQL Server 2005 sp2. It would be very useful for us to be able to monitor our data periodically for defects, with a query such as               select * from sdeDLM.SDE.Lang_Polygons     where ISO_LANGUAGE_CODE_ID <> LEFT(ID , 3)      In ArcMap, I can do that query using Select by Attributes. I enter the above where clause and click 'Verify', and the expression is verified successfully. But in order to monitor the data we need to use the REST API. We can successfully do REST queries with a where clause such as               WHERE sdeDLM.SDE.Lang_Polygons.ISO_LANGUAGE_CODE_ID IS NULL      or just               WHERE ISO_LANGUAGE_CODE_ID IS NULL      Those work just fine. FWIW, the URL for that REST query is:               http://ourServer:8399/arcgis/rest/services/DLM21_query/MapServer/2/query?text=&geometry=&geometryType=esriGeometryEnvelope&inSR=&spatialRel=esriSpatialRelIntersects&where=ISO_LANGUAGE_CODE_ID+is+NULL&returnGeometry=false&outSR=&outFields=&f=html      But when I try a REST query with a where clause like                   WHERE sdeDLM.SDE.Lang_Polygons.ISO_LANGUAGE_CODE_ID <> LEFT(sdeDLM.SDE.Lang_Polygons.ID , 3)      the REST API returns the message:               Unable to perform query. Please check your parameters      As best I can tell, it's the LEFT() function that makes the difference. I have been able to successfully use where clauses with everything else in them, including '<>'. The REST API docs for querying a layer say that "Any legal SQL where clause operating on the fields in the layer is allowed." The URL for the query that gives an error is:               http://ourServer:8399/arcgis/rest/services/DLM21_query/MapServer/2/query?text=&geometry=&geometryType=esriGeometryEnvelope&inSR=&spatialRel=esriSpatialRelIntersects&where=sdeDLM.sde.Lang_Polygons.ISO_Language_Code_ID+%3C%3E+Left%28sdeDLM.sde.Lang_Polygons.ID%2C+3%29&returnGeometry=true&outSR=&outFields=&f=html      Is there a way to get the REST API to accept a where clause using Left()? Or a way to find out why it's not working? I tried browsing the latest logs in C:\Program Files\ArcGIS\server\user\log. The file Server-20100929-143218.dat shows the following output from the time when the query failed:               String request received. Request size is 2629 characters.     QueryFeatureData2 has started.     Preparation for method: QueryFeatureData2 has started.     Preparation for method: QueryFeatureData2 has completed.     Query Operation for method: QueryFeatureData2 on layer 'Language Polygons' has started.     **GeoDatabase Error :An invalid SQL statement was used.**     QueryFeatureData2 has completed.     Method failed.HRESULT = 0x80040207 : This is a FACILITY_ITF error that is specific to the interface that returned this error. See the documentation of the interface that returned this error for information about this HRESULT.     ERROR INFO = An invalid SQL statement was used.      But why is the SQL statement invalid? Clearly the WHERE clause has been tested as valid elsewhere. And why would some of my WHERE clauses be accepted as valid, but not the one with LEFT() in it? Thanks for any help.


title: Moving marker smoothly along two GPS coordinates

text: I am trying to animate a marker from one GPS coordinate to another. I am currently using Google Maps to achieve this and, although my algorithm seems right, Google Maps is doing a poor job at converting the GPS position to the right pixel on the map, which makes the animation not so smooth when the map is zoomed out. Zoomed in, it is very smooth and that's why I suspect the "zig zag" is caused by Google doing a poor LatLng->pixel conversion. The best solution would be for me to deal directly with the markers (x,y) pixel coordinates but there is unfortunately no public method for assigning pixel location to a marker. The code is available here: http://dev.syskall.com/map/ How would you make the animation smoother when zoomed out? Alternatively, I would be interested in using an equivalent to Google Maps that would let me accomplish what I'm trying to do if it isn't possible with Google Maps. PS: I am not sure if this is the right place to ask my question, let me know if it isn't the case.


title: Speed up plug-in datasource in ArcObjects

text: Are there any tricks involved in getting plug-in data sources to run at higher speeds, similar to the built in types? I've noticed that a simple plugin-data source will run many times slower than the equivelent data loaded from a shapefile or personal geodatabase, and the extra time seems to be mainly in the rendering phase. Are there any tricks available to improve this in ArcObjects?


title: Opacity causes tiles to 'flash' in

text: I have a map sitting on top of a dark gray background. I have the canvas basemap set to initially load as semi-transparent:                var canvas = new esri.layers.ArcGISTiledMapServiceLayer("http://services.arcgisonline.com/ArcGIS/rest/services/Canvas/World_Light_Gray_Base/MapServer", {                 "opacity":0.28,                 "id": "canvas"});      When I zoom in/out, the tiles seem to 'flash' in. Easier to see what I mean by watching this short screencast. I am guessing it might be to do with the fact that the layer has some kind of fade in transition for new tiles? If you agree that this is the likely cause, is there any way to turn this effect off? For info, elsewhere in my application I have a transparency slider for this layer, but this has the same value as the initial transparency of the basemap. The reason I mention this, is that if there is a way the tiles can have their transparency altered perhaps by CSS, this will not be a workaround for me, as the tiles can have their opacity changed by the slider.               $('#slider').slider({                 min: 0,                 max: 100,                 value: 28,                 range: "min",                 slide: function(event, ui) {                   map.getLayer("canvas").setOpacity(ui.value/100);               }               });


title: How can I load a single WFS vector feature where the user clicks the map on OpenLayers?

text: I am creating a web app that allows users to easily edit feature attributes as well as geometry. I have editing working by loading a WFS layer, but it slows the app to a crawl when there's more than a couple dozen features on the screen. Instead of loading all of the features as vectors, **I want to display them as a raster layer, then when a feature is clicked, grab the associated vector for that feature** , in order to allow for editing. Has anyone done this before, or know of the best way to go about this?


title: Importing Leica total station .GSI file data into ArcGIS Desktop 9.x?

text: I have a total station data file (test.gsi) that was collected with a Leica 1200. How can I import the data to ArcGIS 9.x and plot points with the point name given by the total station observer. After importing the points i need to draw polygons around the points to define land parcels.


title: ecw on windows 7 qgis 1.8

text: I used 1.7.4 qgis release on Windows 7, by which it was possible to manage ECWfiles; after passing to 1.8 release, I couldn't deal with an ECW file: it tells me that .ecw is not a supported kind of file on QuantumGis. What can I do? Is there anyone that could help me?


title: View an .ecw-file in QGIS 1.8.0?

text: I have got a .ecw-file that I want to somehow import into QGIS 1.8.0. Is there a plugin that I can install that can help me do this or is there another way? The .ecw-file contains rasterdata.


title: GDAL/proj4 equirectangular projection coordinate system units assumed by OGRSpatialReference and GDALDataset::GetGeoTransform()?

text: How do I configure my OGRSpatialReference object to specify an equirectangular projection system http://en.wikipedia.org/wiki/Equirectangular_projection where the affine geo-transform matrix (provided by GDALDataset::GetGeoTransform()) produces angular (WGS84 geodedic latitude & longitude) coordinates? I did not find an answer here: http://www.gdal.org/gdal_datamodel.html or here: How do I convert affine coordinates to lat/lng? or in the GDALDataset and OGRSpatialReference class documentation. My understanding of the equirectangular projection is that it is an angular (latitude & longitude) coordinate system by nature, possibly with scale factors involved. So given the affine transform matrix provided by GDALDataset::GetGeoTransform() for an equirectangular projected dataset, what is the nature of the projection's coordinate system? Is it angle-based (latitude & longitude) or based on virtual linear easting/northing distances that must then be converted to angles? Note I am using a VRT dataset in order to query a UTM projected image as a virtual equirectangular (geographic) projected image. My code:               // open a dataset that uses UTM projection, WGS84 coordinate system     if ((rawDataset = (GDALDataset *)GDALOpen(filename, GA_ReadOnly))==NULL)      (throw exception)      (Output from rawDataset->GetProjectionRef(), etc. all look good. Confirm dataset uses UTM projection.)                // specify a equirectangular projection (also known as geographic or uprojected lat/lon) for the VRT     OGRSpatialReference geographicRef;     geographicRef.SetProjection(SRS_PT_EQUIRECTANGULAR);     // geographicRef.SetEquirectangular(0,0,0,0); // gives same results as above     geographicRef.SetWellKnownGeogCS("WGS84"); // standard WGS84 coordinate system     geographicRef.SetProjCS("Equirectangular / WGS84");     //geographicRef.SetAngularUnits(SRS_UA_DEGREE,0.0174532925199433); // has no effect     //geographicRef.SetAngularUnits(SRS_UA_RADIAN,1.0); // has no effect      geographicRef.exportToWkt(&pszDstWkt);      (Print the projection/coordinate specification string (pszDstWkt). Confirm projection is Equalrectangular and coordinate system is WGS84.               // create the VRT     if ((vrtGeographicDataset = (GDALDataset *)GDALAutoCreateWarpedVRT(rawDataset,(const char *)NULL,pszDstWkt,GRA_Bilinear,10,NULL))==NULL)     (throw exception)      Now examine the original image->UTM and VRT_image->equirectangular geo transform matrices: rawDataset->GetGeoTransform(adfGeoTransform) gives values of:               387080.000000 0.100000 0.000000     3595399.000000 0.000000 -0.100000      The image has 0.1 meter pixel size. So the output of this transform is probably UTM easting/northing units of meters. However, it would be nice to find an API call that tells me the projection coordinate system units are meters. However, vrtGeographicDataset->GetGeoTransform(adfGeoTransform) gives values of:               -11822338.617165 0.108005 0.000000     3616869.349192 0.000000 -0.108005      Clearly the X and Y offsets are too high to imply angular units of degrees or radians. And because 0.108005 is so close to 0.10000 I suspect the units are still proportional to meters. Either I am doing something wrong or perhaps there is an additional scaling factor necessary to convert to angles (degrees or radians). **Update on 2013-12-21** I've learned several things... First GDALAutoCreateWarpedVRT() appears to be ignoring the destination coordinate system defined in char *pszDstWKT. OGRSpatialReference::SetGeogCS() and SetWellKnownGeogCS() have no effect. The absence of SetGeogCS() or SetWellKnownGeogCS() has no effect. Only SetProjection() and SetLinearUnits() appear to have an effect on the VRT's geo transform matrix. Secondly, the VRT's geo transform is closer to what it should be relative to the source file's geo transform. The VRT geo transform has the correct signs while the original file's transform has signs consistent with the file's UTM projection. So GDALAutoCreateWarpedVRT() appears to be changing the projection to Equirectangular or something close to it. Thirdly, via the following code I can obtain a geo transform that provides latitudes and longitudes that are somewhat close to expected values:               OGRSpatialReference geographicRef;     geographicRef.SetProjection(SRS_PT_EQUIRECTANGULAR);     geographicRef.SetLinearUnits("My Degrees Lat,Lon",6378137.0*3.14159265358979/180.0); // equatorial radius     geographicRef.exportToWkt(&pszDstWkt);     vrtGeographicDataset = GDALAutoCreateWarpedVRT(rawDataset,NULL,pszDstWkt,...     vrtGeographicDataset->GetGeoTransform(adfGeoTransform);      Note pszDstWkt contains the definition: 'PROJCS["unnamed",PROJECTION["Equirectangular"],UNIT["My Degrees Lat,Lon",111319.4907932735]]' The source file's definition is: 'PROJCS["WGS 84 / UTM zone 13N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.2572235630016, The VRT's geo transform is now:               -106.20187473836      9.70225388411317e-07                         0     32.4908901704265                         0     -9.70225388411317e-07      At pixel (0,0), the longitude of -106.20187473836 is exactly correct for my source image, verified by both the Lizardtech GeoViewer and Google Maps. However, the (0,0) pixel's latitude should be 32.490109, not 32.490890. And both diagonal coefficients of 9.70225388411317e-07 are somewhat too large. So all other pixels besides the upper left corner are in error. Every other configuration of OGRSpatialReference I've tried gives worse results. And I don't see a more appropriate projection than equirectangular. I did try to force a spherical definition (instead of ellipsoidal) by setting a high inverse flattening in geographicRef.SetGeogCS(). But again GDALAutoCreateWarpedVRT appears to be ignoring the PROJCS section of the definition altogether. So I am not sure how to achieve this. Please help!! My goal is to extract raster sub-images (ROIs) with an affine geo transform that converts pixel (column,row) coordinates to WGS84 geodedic (longitude,latitude) coordinates. And I need to do this with GDAL's C++ API. The source files are too large and too numerous to use gdal_translate or gdalwarp.


title: ArcGIS 10.2 Installation for 64-bit geoprocessing

text: I am trying to install 64-bit processing to ArcGIS 10.2, but the ESRI Patches and Service Packs website only has installations for Service pack 1 that goes with ArcGIS 10.1. One similar question on this exchange shows a screenshot of an installation window that shows an option to include this, but my setup installation program (which I downloaded) does not show this installation window, but instead goes straight to selecting features to install as shown below. I can't find any other documentation for running 64-bit geoprocessing on ArcGIS version 10.2. Can anyone show me how to do this? or is it best to just try to install the service pack 1 version of it? ![enter image description here](http://i.stack.imgur.com/TYmlx.png)


title: Call function inside dojo.ready

text: I have code ArcGIS Javascript like this:               <script type="text/javascript">           require(["dojo/ready",           "dgrid/OnDemandGrid",           "dgrid/Selection",           "dojo/store/Memory",            "dojo/_base/array",           "dojo/dom-style",           "dijit/registry",           "esri/map",            "esri/dijit/Popup",           "esri/dijit/PopupTemplate",           "esri/layers/FeatureLayer",           "dojo/dom-class",           "dojo/dom-construct",           "dojo/on",            "esri/InfoTemplate",           "esri/symbols/SimpleLineSymbol",           "esri/symbols/SimpleFillSymbol",           "esri/tasks/QueryTask",           "esri/tasks/query",           "dojo/_base/declare",            "dojo/number",            "dojo/on",           "dojox/charting/Chart",           "esri/renderers/UniqueValueRenderer",           "dojo/parser",           "dojo/_base/Color",            "dijit/layout/BorderContainer",           "dijit/layout/ContentPane"           ], function (             ready,             Grid,              Selection,              Memory,              array,             domStyle,             registry,             Map,             Popup,             PopupTemplate,              FeatureLayer,             domClass,             domConstruct,             on,              InfoTemplate,             SimpleLineSymbol,             SimpleFillSymbol,             QueryTask,             Query,             declare,              dojoNum,              on,             Chart,              UniqueValueRenderer,             parser           ) {             ready(function() {                    parser.parse();                    var popup = Popup({                 titleInBody: false               },domConstruct.create("div"));                    // create the dgrid               window.grid = new (declare([Grid, Selection]))({                 // use Infinity so that all data is available in the grid                 bufferRows: Infinity,                 columns: {                   "FID": "FID",                   "PROV": "PROV"                 }               }, "grid");               // add a click listener on the ID column               grid.on(".field-id:click", selectState);                    var initialExtent = new esri.geometry.Extent({"xmin":89.971,"ymin":-12.472,"xmax":144.301,"ymax":9.511,"spatialReference":{"wkid":4326}});               map = new esri.Map("map", { extent: initialExtent,infoWindow: popup});                    //var basemap = new esri.layers.ArcGISTiledMapServiceLayer("http://services.arcgisonline.com/ArcGIS/rest/services/Reference/World_Boundaries_and_Places/MapServer");               var basemap2 = new esri.layers.ArcGISDynamicMapServiceLayer("http://localhost:6080/arcgis/rest/services/BaseMap/Peta_Dasar_Indonesia_Colors/MapServer",{"opacity":0.4});               map.addLayer(basemap2);               //map.addLayer(basemap);                         domClass.add(window.map.infoWindow.domNode, "myTheme");                 var template = new PopupTemplate({                 title: "Data Luas Lahan Sawah {PROV}",                 description: "of starters from {PROV} finished",               });                    window.statesUrl = "http://localhost:6080/arcgis/rest/services/BaseMap/Peta_Dasar_Indonesia_Colors/MapServer/8";               window.outFields = ["FID","PROV"];                    var fl = new FeatureLayer(window.statesUrl, {                 id: "stat",                 mode: 1, // ONDEMAND, could also use FeatureLayer.MODE_ONDEMAND                 infoTemplate:template,                 outFields: window.outFields               });                    fl.on("load", function() {                 fl.maxScale = 0; // show the states layer at all scales                 fl.setSelectionSymbol(new SimpleFillSymbol().setOutline(null).setColor("#AEC7E3"));               });                    fl.on("click", selectGrid);                    // change cursor to indicate features are click-able               fl.on("mouse-over", function() {                 map.setMapCursor("pointer");               });               fl.on("mouse-out", function() {                 map.setMapCursor("default");               });                         map.addLayer(fl);                    map.on("load",                  function addFeatureLayer() {                   var defaultSymbol = new SimpleFillSymbol().setStyle(SimpleFillSymbol.STYLE_NULL);                   defaultSymbol.outline.setStyle(SimpleLineSymbol.STYLE_NULL);                        //create renderer                   var renderer = new UniqueValueRenderer(defaultSymbol, "PROV");                        //add symbol for each possible value                   renderer.addValue("Banten", new SimpleFillSymbol().setColor("#00FA1D"));                   renderer.addValue("Dki Jakarta", new SimpleFillSymbol().setColor("#FA0000"));                   renderer.addValue("Jawa Barat", new SimpleFillSymbol().setColor("#FAFA00"));                   renderer.addValue("Jawa Tengah", new SimpleFillSymbol().setColor("#00FA1D"));                   renderer.addValue("Daerah Istimewa Yogyakarta", new SimpleFillSymbol().setColor("#FAFA00"));                   renderer.addValue("Jawa Timur", new SimpleFillSymbol().setColor("#FAFA00"));                   var template = new PopupTemplate({                     title: "Data {PROV}",                     description: "Data Luas Lahan Sawah {PROV}",                   });                   var featureLayer = new FeatureLayer("http://localhost:6080/arcgis/rest/services/BaseMap/Peta_Dasar_Indonesia_Colors/MapServer/8", {                     mode: 1, // ONDEMAND, could also use FeatureLayer.MODE_ONDEMAND                     infoTemplate:template,                     outFields: window.outFields                   });                        featureLayer.setRenderer(renderer);                   map.addLayer(featureLayer);                 } ,                 function( evt ){                   // show the border container now that the dijits                    // are rendered and the map has loaded                   domStyle.set(registry.byId("container").domNode, "visibility", "visible");                   populateGrid(Memory); // pass a reference to the MemoryStore constructor               });                    function populateGrid(Memory) {                 var qt = new QueryTask(window.statesUrl);                 var query = new Query();                 query.where = "1=1";                 query.returnGeometry = false;                 query.outFields = window.outFields;                 qt.execute(query, function(results) {                   var data = array.map(results.features, function(feature) {                     return {                       // property names used here match those used when creating the dgrid                       "FID": feature.attributes[window.outFields[0]],                       "PROV": feature.attributes[window.outFields[1]]                     }                   });                   var memStore = new Memory({ data: data });                   window.grid.set("store", memStore);                 });               }               // fires when a row in the dgrid is clicked               function selectState(e) {                 // select the feature                 var fl = map.getLayer("stat");                 var query = new Query();                 query.objectIds = [parseInt(e.target.innerHTML)];                 fl.selectFeatures(query, FeatureLayer.SELECTION_NEW, function(result) {                   if ( result.length ) {                     // re-center the map to the selected feature                     window.map.centerAt(result[0].geometry.getExtent().getCenter());                   } else {                     console.log("Feature Layer query returned no features... ", result);                   }                 });               }                    // fires when a feature on the map is clicked               function selectGrid(e) {                 var profinsi = e.graphic.attributes.PROV;                 var id = e.graphic.attributes.FID;                 // select the feature that was clicked                 var query = new Query();                 query.objectIds = [id];                 var states = map.getLayer("stat");                 states.selectFeatures(query, FeatureLayer.SELECTION_NEW);                 // select the corresponding row in the grid                 // and make sure it is in view                 grid.clearSelection();                 grid.select(id);               }                    function test(){                 alert("alert");               }             });         });     </script>      I want call function inside dojo.ready function. For example, I want call function test() with onclick with this:               <input type="submit" value="Submit" onclick="test();">      but still not work.   Can anybody help me explain how to call function inside dojo.ready?


title: Using a bounding polygon, instead of a bounding box

text: I've got a large set of LAT/LON coordinates. Up until now, when I needed to select out those that fell within a certain region, I'd just create a bounding box using some MINx, MAXx, MINy, MAXy values and try to get it as close as possible to the area I was interested in. I'd like to get more precise by using more than 4 points as the extent of my search area, in this specific case, 14 points. I'll probably be using Python or SQL for this, so the specific syntax on any solution isn't that important, I'm just interested in the formula for figuring this one out.


title: How to implement pgRouting in a web GIS based on MapServer?

text: I have 3d lines in PostGIS database (representing hiking trails). I'm using phpPgAdmin administration tool. I want to use Dijkstra method for finding shortest way. I inserted all needed columns to execute this function, for example:                   SELECT * FROM shortest_path('        SELECT gid AS id,               start_id::int4 AS source,               end_id::int4 AS target,               cost_length::float8 AS cost        FROM network',     1,     135,     false,     false)      I have .html file and .map file. I alredy have some functions and layers installed in my web application (like layer for shaded relief map). Also, I already have defined layer for 3d hiking trails from PostGIS database and it is shown on the web map by checkbox. I never worked with pgRouting and MapServer before and I need help how to implement this function in my web application. I need detailed answer on **how to implement it in my .map file and .html file**. I need to enable users to choose one starting node (I will name them, but I don't know if I need to name it all?) from the combobox I want to define in my .html file, and one end node, also from the combobox, and then enable showing shortest path, or maybe list of nodes where they need to go for shortest path from start point to end point. Don't give me links like [1]: http://users.qgis.org/planet/user/4/tag/pgrouting/ because it doesn't give me an explanation about working with MapServer and pgRouting. It is very important and emergency, because it is about my graduation thesis. Thank you! * * * This is error when I'm trying to open my application:               writeLog(): Unable to access file. mslog.log msDrawMap(): Image handling error.      Failed to draw layer named 'routing'. msPostGISLayerWhichShapes(): Query error.      Error (ERROR: syntax error at or near "SELECT" LINE 1: ...) as geom,"gid" from      (SELECT * FROM shortest_path(SELECT the... ^ ) executing query:      select encode(ST_AsBinary(ST_Force_2D("the_geom"),'NDR'),'hex') as geom,"gid" from      (SELECT * FROM shortest_path(SELECT the_geom as the_geom, gid AS id, %start_id%::int4      AS source, %end_id%::int4 AS target, cost_length::float8 AS cost FROM network, 1, 135,      false, false)) AS foo where the_geom && GeomFromText('POLYGON((561819.89148581          5072537.66277129,561819.89148581 5093307.66277129,589524.782971619      5093307.66277129,589524.782971619 5072537.66277129,561819.89148581      5072537.66277129))',4326)


title: gdalbuildvrt error when mosaicking SRTM data with different NoData values

text: I'm attempting to mosaic together 9 SRTM tiles downloaded from CGIAR to produce a DEM of Borneo. I'm using:               gdalbuildvrt Borneo.vrt BorneoSRTM/*.tif      where BorneoSRTM contains the 9 tif files to be combined. The command runs without error, however the resulting DEM looks like the below image and all the elevations are less than 13m (some of the input files have elevations close to 4000 m). ![enter image description here](http://i.stack.imgur.com/FnfeU.png) Clearly there's a problem with the tile corresponding to the big white square (5-10N, 110-115E). Looking at the original tif files and I've noticed that NoData is 255 for the problem tile and -32768 for all others. I tried changing NoData to -32768 using:               gdal_translate -a_nodata -32768  srtm_59_11.tif  out.tif      but this didn't solve the problem. Calling gdalinfo on the problem tile gives:               Driver: GTiff/GeoTIFF     Files: srtm_59_11.tif            srtm_59_11.tfw     Size is 6001, 6001     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (109.999583817610983,10.000416763514124)     Pixel Size = (0.000833333333333,-0.000833333333333)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=BAND     Corner Coordinates:     Upper Left  ( 109.9995838,  10.0004168)      Lower Left  ( 109.9995838,   4.9995834)      Upper Right ( 115.0004172,  10.0004168)      Lower Right ( 115.0004172,   4.9995834)      Center      ( 112.5000005,   7.5000001)      Band 1 Block=6001x1 Type=Int16, ColorInterp=Gray       NoData Value=255      And, calling it on another tile gives:               Driver: GTiff/GeoTIFF     Files: srtm_59_12.tif            srtm_59_12.tfw     Size is 6001, 6001     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (109.999583817610983,5.000416642442191)     Pixel Size = (0.000833333333333,-0.000833333333333)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=BAND     Corner Coordinates:     Upper Left  ( 109.9995838,   5.0004166)      Lower Left  ( 109.9995838,  -0.0004167)      Upper Right ( 115.0004172,   5.0004166)      Lower Right ( 115.0004172,  -0.0004167)      Center      ( 112.5000005,   2.5000000)      Band 1 Block=6001x1 Type=Int16, ColorInterp=Gray       NoData Value=-32768      When I call gdalbuildvrt without this problem tile, everything works fine. I realize that this tile only includes a tiny piece of Borneo, but if possible I'd like to figure out a way to include it. Any help would be appreciated.


title: How to use Patch in Geoserver

text: I have got a patch which supports Nvarchar and Nchar types when we need to publish sql server table to Geoserver. But I dont know, How to use this patch OR where i put this patch so that it can work with Geoserver ? Can anybody explain, how i use this patch with Geoserver?


title: How to generate discrete areas from multiple, overlapping features?

text: I have a point dataset that I want to append some information from a polygon to. My polygon shapefile contains over 500 features. Whilst each of these features represents individual areas they are not discrete in their extent. For example, (see highly sophisticated Paint mock up below) a point (point 1) in Figure A may actually overlay all three zones (Point 1, Figure B). A spatial join or query will then result in an incorrect zone being appended. ![Simple overview of the problem](http://i.stack.imgur.com/jBWQt.jpg) What I’m wanting is a dataset like that shown in Figure C i.e. discrete polygons that do not overlap. I have tried splitting my original polygon by an attribute/location and unioning these ~500 shapefiles together but to no avail. Namely, the Union tool would like to do this process two shapefiles at a time (I ideally wanted to point the tool at a folder containing the split shapefiles and union them all together). This should be a really easy problem to solve but I want to avoid doing this manually (i.e. through the Editor 'Clip' tool). **Does anyone have any ideas?**


title: Why does Add Field to Layer File give "ERROR 000840: The value is not a Feature Layer"?

text: I'm trying to add a field to a layer file, but it always gives an error. What I did was making a closest facility layer and save the result as a layer and after that I tried to add a field, but then it gives an error. It says it cannot add a field if it isn't a raster layer, ... I tried to avoid this by making a shapefile (manually) and this can work, but is there a method to do this in Python? So convert a route from the closest facility layer to a shapefile using Python? It always says: ExecuteError: Failed to execute. Parameters are not valid.   ERROR 000840: The value is not a Feature Layer. Thanks! Regards


title: ESRI GDB Support in GDAL / ogr2ogr on Mac OS X 10.8

text: How do you incorporate support for ESRI's fgdb format into GDAL / ogr2ogr on Mac OS X 10.8?


title: SQL Server STpointfromText comes from PHP into javascript as garbage

text: Sorry for the sloppy title. Here's the scenario - I have a table in `SQL Server 2008 R2` It has a geographyColumn, created thusly:               add geographyColumn as geography::STPointFromText('POINT('+convert(varchar(20),x_coord)+'          '+convert(varchar(20),y_coord)+')',4326)      When I viewed it in the table, it looks fine:               POINT (-116.477654 33.764695)      I have some javascript, and a PHP script, to grab points from said table. I can get any text field without a problem. But I want to get the X,Y vales from the geography column. Instead of a nicely formed "X,Y" or even a long string of numbers, I get this mess:               ����zU���f�C@      My php is the following:               <?php $q=$_GET["q"];      $serverName = "ServerIP";      $connectionInfo = array( "Database"=>"geo_spat2", "UID"=>"u", "PWD"=>"p");      $con = sqlsrv_connect($serverName, $connectionInfo);      if (!$con){         die( print_r( sqlsrv_errors(), true));      };      $sql="SELECT * FROM EMPTY_REPS WHERE Division_label = '".$q."'";      $result = sqlsrv_query($con, $sql);      while($row = sqlsrv_fetch_array($result)) {        echo $row['geographyColumn'] . ", ";      }      sqlsrv_close($con);      ?>      Am I missing something simple? Thanks, Evan


title: Export to new projection in QGIS 1.7

text: Suppose I have layers A and B. I don't know their projections but by the looks of the map I know they are not the same. In QGIS 1.6 I would use one of my favorite features:               Vector> Data Management Tools> Export to New Projection      What was nice about that feature is that you could export B to A's projection by selecting "choose projection from existing layer". Hence you need not know A or B's projection. In QGIS 1.7 the `Data management Tools` menu is much reduced and the `Export to New Projection` is missing. Yes, I can do projections on the fly, and I can do `save as`. The former is not useful if I want to convert B's projection to A's. The latter does not give option to set projection from open layer. I have tried to find different ways to access the `Export to New Projection` dialog without success. Any tips?


title: How to dissolve adjacent parcels by same owner but include those parcels located across the street with the same owner?

text: In trying to determine the total contiguous stream length per landowner, I dissolved parcels by parcel owner name attribute and intersected these dissolved parcels with streams. However, dissolving by parcel owner doesn't include parcels located across the street (right-of-way) when 'Create multipart features' is unchecked since they don't share common edges; and when I do check 'Create multipart features', it creates the desired result across the street but also includes piecemealed parcels per owner that are spread out in space which are not considered contiguous. I was even thinking of modeling this somehow in Spatial Analyst. I'm using ArcGIS 10.1 with a basic license. Any ideas of how to achieve the desired result?


title: Using GDAL and Python to replace the Spatial Analyst TabulateArea method

text: Has anyone tried to phase out the use of ESRI's Spatial Analyst extension by replacing it with GDAL in Python? I have a script that only uses the TabulateArea and Statistics analysis but not much else from SA. Would be great to take ESRI out of the picture there and not have to pay for the SA extension. Anybody have any luck with this?


title: ToolBar icon images missing

text: very early build ~ zip file I am getting blank outputs for button images in a custom toolbar. Although the button information populates the dockable toolbar, they lack respective icons [see Fig.1] Thanks for any help or clues! Much of the following code was referenced from Add toolbar via Python               def initGui(self):         # Add toolbar              self.toolbar = self.iface.addToolBar("Graphic Groundwater")         self.toolbar.setObjectName("Graphic Groundwater")              # Create actions         self.getDEM = QAction(QIcon(":/plugins/GraphicGroundwater/icons/getDEM.png"),                                QCoreApplication.translate("Groundwater Modeling", "Get DEM Map"),                                self.iface.mainWindow())         self.getTopo = QAction(QIcon(":/plugins/GraphicGroundwater/icons/getTopo.png"),                                QCoreApplication.translate("Groundwater Modeling", "Get Topo Map"),                                self.iface.mainWindow())         self.addBase = QAction(QIcon(":/plugins/GraphicGroundwater/icons/addBase.png"),                                QCoreApplication.translate("Groundwater Modeling", "Add Base Area"),                                self.iface.mainWindow())         self.addNewBase = QAction(QIcon(":/plugins/GraphicGroundwater/icons/addNewBase.png"),                                QCoreApplication.translate("Groundwater Modeling", "Add New Base Area"),                                self.iface.mainWindow())         self.addRiver = QAction(QIcon(":/plugins/GraphicGroundwater/icons/addRiver.png"),                                QCoreApplication.translate("Groundwater Modeling", "Add River"),                                self.iface.mainWindow())         self.addStress = QAction(QIcon(":/plugins/GraphicGroundwater/icons/addStress.png"),                                QCoreApplication.translate("Groundwater Modeling", "Add Stress"),                                self.iface.mainWindow())         self.addGrid = QAction(QIcon(":/plugins/GraphicGroundwater/icons/addGrid.png"),                                QCoreApplication.translate("Groundwater Modeling", "Add Grid"),                                self.iface.mainWindow())         self.refreshGrid = QAction(QIcon(":/plugins/GraphicGroundwater/icons/refreshGrid.png"),                                QCoreApplication.translate("Groundwater Modeling", "Refresh Attributes"),                                self.iface.mainWindow())         self.runModel = QAction(QIcon(":/plugins/GraphicGroundwater/icons/runModel.png"),                                QCoreApplication.translate("Groundwater Modeling", "Run Model"),                                self.iface.mainWindow())              # Connect action signals to slots         self.getDEM.triggered.connect(self.run)         self.getTopo.triggered.connect(self.run)         self.addBase.triggered.connect(self.run)         self.addNewBase.triggered.connect(self.run)         self.addRiver.triggered.connect(self.run)         self.addStress.triggered.connect(self.run)         self.addGrid.triggered.connect(self.run)         self.refreshGrid.triggered.connect(self.run)         self.runModel.triggered.connect(self.run)              # Add actions to the toolbar         self.toolbar.addAction(self.getDEM)         self.toolbar.addAction(self.getTopo)         self.toolbar.addAction(self.addBase)         self.toolbar.addAction(self.addNewBase)         self.toolbar.addAction(self.addRiver)         self.toolbar.addAction(self.addStress)         self.toolbar.addAction(self.addGrid)         self.toolbar.addAction(self.refreshGrid)         self.toolbar.addAction(self.runModel)          def unload(self):         # remove toolbar on plugin unload         del self.toolbar      [Fig.1] ![blank toolbar](http://i.stack.imgur.com/FhzpS.png) The associated resource_rc file might also be of some diagnosis assistance...               <qresource prefix=":/plugins/graphicgroundwater/" >         <file>GraphicGroundwater.png</file>         <file>getDEM.png</file>         <file>getTopo.png</file>         <file>addBase.png</file>         <file>addNewBase.png</file>         <file>addRiver.png</file>         <file>addStress.png</file>         <file>addGrid.png</file>         <file>refreshGrid.png</file>         <file>runModel.png</file>     </qresource>