title: When is a projected length calculated via a pythagorean distance more accurate than an ellipsoidal length?

text: I'm using spatialite 3.0.1 with the following test data (in EPSG: 4326):               Name        astext(geometry)         ----------  --------------------     Belitsa     POINT(23.563 41.947)     Malko Tarn  POINT(27.533 41.983)     Kotel       POINT(26.4752 42.938      My client demands I provide the projected st_length distance instead of geodesiclength. (compare the following queries)               select a.name, b.name, st_length(transform(makeline(a.geometry, b.geometry), 32635)) from ProjCist a, ProjCist b;          Belitsa     Malko Tarn  329056.08890108      ...      versus               select a.name, b.name, geodesiclength(makeline(a.geometry, b.geometry)) from ProjCist a, ProjCist b;          Belitsa     Malko Tarn  329092.535472995      ...      The 50 meter difference in length, for our purposes, is _hugely_ significant. Which is more correct? Why?


title: Problems editing geometry in JS script using a find task in a web application

text: I've generated a javascript code to my app based on the sample "Show find task results in a DataGrid" of ArcGIS API for JavaScript and I am using editable layers (WFS-T) and it works, but when after use it I try to select the data to change the geometry and edit it only allows to me edit the attributes, but I can't select the vertexs/nodes of the geometry by double-clicking. This trouble only appears when I make a search, but not before use it. How can I make my data editable again? I've eliminated the graphics with map.graphics.hide() that I've generated before for the query and to store the data on the grid, and it works, but I've the problem with editing geometries. If anyone can help me it would be very helpful! Thanks


title: OpenLayers : Multiple layers popup issue

text:
I'm actually having a big issue that I cannot solve so far. Here is my problem. I actually have a map on which I'm displaying markers (WMSGetFeatureInfo) depending on some categories that I want to show or hide. But here is the thing, I actually get those markers as a ***.map** file which is actually a PNG file containing all the markers. When I'm clicking on a marker, I'm getting the position of the click and then add a popup. But the more layers I have display on the map, the more the popup is pointing far away from the related marker. Here some code :                   var layer = new OpenLayers.Layer.WMS(             category.label,             "/mapserv?Map=" + category.url,             {                 layers  : category.layer,                 format  : 'image/png',                 version : '1.3.0',                 srs     : 'ESPG:3163'             },             {                 isBaseLayer : false,                 singleTile  : true,                 visibility  : visibility,             }         );              var info = new OpenLayers.Control.WMSGetFeatureInfo({             title          : 'get details by clicking',             layers         : [ layer ],             infoFormat     : "text/plain",             queryVisible   : true,             eventListeners : {                 getfeatureinfo : function( event ){                     document.body.style.cursor = 'auto';                          var getId = function( text ) {                         result = context.settings.pattern.exec( text );                         if( result == null ) return;                         return result[1];                     }                          var id = getId( event.text );                          request = jQuery.ajax({                         url     : "/cartoweb/FicheTheme",                         type    : "get",                         data    : "idGeoEad="+id,                         success : function(response, textStatus, jqXHR){                             if (response == null || response['data'] == null) {                                 return;                             }                             context.generateAndShowPopup( event.xy, response['data']['metaData'] );                         }                     });                 },                 beforegetfeatureinfo : function( event ){                     document.body.style.cursor = 'wait';                 },                 nogetfeatureinfo : function( event ){                     document.body.style.cursor = 'auto';                 }             }         });              this.map.addControl( info );         info.activate();               generateAndShowPopup : function( latlong, text ) {         var lonlatfrompx = this.map.getLonLatFromViewPortPx( latlong );         var anchor       = {             'size'      : new OpenLayers.Size(0,0),              'offset'    : new OpenLayers.Pixel(-36, 6),              'keepInMap' : true         };              // Hide by default popup actually open         if(this.popup !== undefined) {             this.popup.hide();         }              // Create a new popup         this.popup = new OpenLayers.Popup.Anchored(             "chicken",             lonlatfrompx,             new OpenLayers.Size( 2000, 2000 ),             '

' + text + '
',             anchor,             false,              function(){}         );              // Popup settings         this.popup.setBackgroundColor( 'transparent' );         this.popup.panMapIfOutOfView         = true;         this.popup.calculateRelativePosition = function () { return 'tr'; }              // Add it on the map         this.map.addPopup( this.popup );              var that = this;         setTimeout(function(){             that.popup.updateSize();         }, 50);              this.map.setCenter( lonlatfrompx );     }      Not sure if I'm doing something wrong, but if someone already faced the same issue, that will be really helpful to know where I did something wrong. Thanks a lot **EDIT**   Ok so I've tried to reproject the latlong on the layer but it still not working. However, I'm making progress :               // Layer is the layer associated to each WMSGetFeatureInfo     // this.layer is the base layer used to display the map     lonlatfrompx.transform( layer.projection, this.layer.projection );     this.map.setCenter( lonlatfrompx );      **EDIT 1**   Ok so I'm making progress. I think I know where the problem comes from but I still don't see/know how I'm gonna solve it (yet). When intanciating the `info` variable, I'm making an ajax call inside which will show a popup when click. The thing is that the more I add layer, the more this ajax request is trying to get information. Let's say, I show on the map 5 differents layers (IT, design, startup, computer and desktop). Then when I'm clicking on a markers (WMSGetFeatureInfo), it actually try to make an ajax `request * [number of layers visible]` and then change at some point (sometimes), the position value. In order to solve this, I will need to avoid this ajax request to be executed more than one time. Any ideas?


title: Why do I lose data after converting from .xlsx to .shp?

text: I have a fairly large excel file that I need to convert to a shapefile, but everytime I export the data, I lose a bunch of it. The excel file has about 277,000 records, but the shapefile will only contain about 260,000 records. I've tried to convert the .xlsx to .csv, which only lost about 2,000 records the first time I exported it, but lost about 20,000 the second time. Why am I losing data during this conversion?


title: QGIS 2.0 Python Error on startup in Ubuntu 10.04

text: I upgraded QGIS to 2.0.1 through the official repository in Ubuntu 10.04 64 bit version. It would start but would not load ftools and many other plugins. It would also not display the "Get more" option under the "Manage plugins" menu. Assuming its to do with Python I upgraded it to 2.6.5. Ever since, QGIS throws up an error message on startup and exits. Somebody suggested that it might be to do with permission of the .qgis2 folder. However, I have checked and verified that the owner of this folder is 'my user' and NOT root. Any tips/advice on how to resolve this would be highly appreciated. This is the message I get during startup: Couldn't load plugin 'processing' from ['/usr/share/qgis/python', '/home/sachin/.qgis2/python', '/home/sachin/.qgis2/python/plugins', '/usr/share/qgis/python/plugins', '/usr/lib/python2.6', '/usr/lib/python2.6/plat-linux2', '/usr/lib/python2.6/lib-tk', '/usr/lib/python2.6/lib-old', '/usr/lib/python2.6/lib-dynload', '/usr/lib/python2.6/dist-packages', '/usr/lib/python2.6/dist-packages/PIL', '/usr/lib/python2.6/dist-packages/gst-0.10', '/usr/lib/pymodules/python2.6', '/usr/lib/python2.6/dist-packages/gtk-2.0', '/usr/lib/pymodules/python2.6/gtk-2.0', '/usr/lib/python2.6/dist- packages/wx-2.8-gtk2-unicode', '/usr/local/lib/python2.6/dist-packages', '/usr/share/qgis/python/plugins/fTools/tools'] Traceback (most recent call last): File "/usr/lib/python2.6/dist- packages/qgis/utils.py", line 182, in loadPlugin **import** (packageName) File "/usr/lib/python2.6/dist-packages/qgis/utils.py", line 453, in _import mod = _builtin_import(name, globals, locals, fromlist, level) File "/usr/share/qgis/python/plugins/processing/ **init**.py", line 20, in from processing.tools.general import runalg, runandload, alghelp, alglist, algoptions, load, extent, getobject File "/usr/lib/python2.6/dist- packages/qgis/utils.py", line 453, in _import mod = _builtin_import(name, globals, locals, fromlist, level) File "/usr/share/qgis/python/plugins/processing/tools/general.py", line 29, in from processing.core.Processing import Processing File "/usr/lib/python2.6/dist- packages/qgis/utils.py", line 453, in _import mod = _builtin_import(name, globals, locals, fromlist, level) File "/usr/share/qgis/python/plugins/processing/core/Processing.py", line 44, in from processing.grass.GrassAlgorithmProvider import GrassAlgorithmProvider File "/usr/lib/python2.6/dist-packages/qgis/utils.py", line 453, in _import mod = _builtin_import(name, globals, locals, fromlist, level) File "/usr/share/qgis/python/plugins/processing/grass/GrassAlgorithmProvider.py", line 34, in from processing.grass.GrassAlgorithm import GrassAlgorithm File "/usr/lib/python2.6/dist-packages/qgis/utils.py", line 453, in _import mod = _builtin_import(name, globals, locals, fromlist, level) File "/usr/share/qgis/python/plugins/processing/grass/GrassAlgorithm.py", line 31, in import importlib File "/usr/lib/python2.6/dist-packages/qgis/utils.py", line 453, in _import mod = _builtin_import(name, globals, locals, fromlist, level) ImportError: No module named importlib Python version: 2.6.5 (r265:79063, Sep 26 2013, 18:58:43) [GCC 4.4.3] QGIS version: 2.0.1-Dufour Dufour, exported Python path: ['/usr/share/qgis/python', '/home/sachin/.qgis2/python', '/home/sachin/.qgis2/python/plugins', '/usr/share/qgis/python/plugins', '/usr/lib/python2.6', '/usr/lib/python2.6/plat-linux2', '/usr/lib/python2.6/lib-tk', '/usr/lib/python2.6/lib-old', '/usr/lib/python2.6/lib-dynload', '/usr/lib/python2.6/dist-packages', '/usr/lib/python2.6/dist-packages/PIL', '/usr/lib/python2.6/dist- packages/gst-0.10', '/usr/lib/pymodules/python2.6', '/usr/lib/python2.6/dist- packages/gtk-2.0', '/usr/lib/pymodules/python2.6/gtk-2.0', '/usr/lib/python2.6/dist-packages/wx-2.8-gtk2-unicode', '/usr/local/lib/python2.6/dist-packages', '/usr/share/qgis/python/plugins/fTools/tools'] An error occured during execution of following code: pyplugin_installer.initPluginInstaller() Traceback (most recent call last): File "", line 1, in File "/usr/share/qgis/python/pyplugin_installer/installer.py", line 43, in initPluginInstaller pluginInstaller = QgsPluginInstaller() File "/usr/share/qgis/python/pyplugin_installer/installer.py", line 58, in **init** repositories.load() File "/usr/share/qgis/python/pyplugin_installer/installer_data.py", line 329, in load url = settings.value(key+"/url", "", type=unicode) TypeError: QSettings.value(QString, QVariant defaultValue=QVariant()): 'type' is not a valid keyword argument Python version: 2.6.5 (r265:79063, Sep 26 2013, 18:58:43) [GCC 4.4.3] QGIS version: 2.0.1-Dufour 'Dufour', exported Python path: ['/usr/share/qgis/python', '/home/sachin/.qgis2/python', '/home/sachin/.qgis2/python/plugins', '/usr/share/qgis/python/plugins', '/usr/lib/python2.6', '/usr/lib/python2.6/plat-linux2', '/usr/lib/python2.6/lib-tk', '/usr/lib/python2.6/lib-old', '/usr/lib/python2.6/lib-dynload', '/usr/lib/python2.6/dist-packages', '/usr/lib/python2.6/dist-packages/PIL', '/usr/lib/python2.6/dist- packages/gst-0.10', '/usr/lib/pymodules/python2.6', '/usr/lib/python2.6/dist- packages/gtk-2.0', '/usr/lib/pymodules/python2.6/gtk-2.0', '/usr/lib/python2.6/dist-packages/wx-2.8-gtk2-unicode', '/usr/local/lib/python2.6/dist-packages', '/usr/share/qgis/python/plugins/fTools/tools'] Regards Sachin


title: Feature Extraction from LiDAR in Erdas 2012 or ArcGIS10.1

text: I know it's possible but I just can't find or figure out a workflow for it. Erdas Webinars has the workflow for ver 2011 but this is very different in 2013 as it uses a raster img file instead. I also looked at the workflow files from Imagine Help. In ESRI it's a similar issue. Any advise on workflows for this?   We own both these software so I would like to manually extract defects from a terrestrial LiDAR scan using either of these.


title: How to gain the facility of snapping when multiple points are selected

text: I have created a fishnet label(square 961 points) centering a point in other layer in arc map but the fishnet is deviated from the intended center. Now i want to place the midpoint of labels on the intended center by selecting all the labels. But when i select all label points, i do not get the facility of snapping whereas i put the midpoint of the labels on top of the intended center accurately. Whenever i zoom the map i found midpoint of fishnet labels and intended center did no match.                                               Present situation      ![Present situation](http://i.stack.imgur.com/DnSuD.png)                                               Wanted situation      ![Wanted situation](http://i.stack.imgur.com/6Mv42.png) Both are in the same coordinate system. How can i get rid of this situation . Thanks in advance


title: Select and zoom python script

text: I am attempting to make a simple select and zoom python script for my address points. Currently the script will complete successfully but will not select or zoom to the entered point. I am using two parameters that have been set in the Script Tool. Address_Number (Long) and Address_Name (String) which both have identical spellings/data types in Properties>>Parameters Here's the code:               import arcpy     mxd = arcpy.mapping.MapDocument("CURRENT")     df = arcpy.mapping.ListDataFrames(mxd, "*")[0]     addLayer = arcpy.mapping.Layer("G:\Layers\GIS_SDE.AddressPoints.lyr")     arcpy.mapping.AddLayer(df, addLayer, "TOP")     Address_Number = arcpy.GetParameterAsText(0)     Address_Name = arcpy.GetParameterAsText(1     arcpy.SelectLayerByAttribute_management("GIS_SDE.AddressPoints", "NEW_SELECTION", "A_NBR = %Address_Number% AND ADDR = %Address_Name%")     df.zoomToSelectedFeatures()      Any help is appreciated. Thanks, Jordan


title: Trouble with accesing API code: "ReferenceError: esri is not defined"

text: I am developing a web map application in the ArcGIS API for Javascript. I two have ArcGISDynamicMapServiceLayer layers that I am trying to run a simple query task on. The code I am using works in other applications doing the same task with different data. My problem is, when I load the app (tried both Firefox and Chrome) in a browser, the map itself won't load (indefinite load screen). Console in both browsers shows this error: > ReferenceError: esri is not defined. Comment #13 on this Esri forum post alludes to the DOM issue, but I need help understanding how to handle the issue. Also, it's possible that the DOM thing is not the issue at all. So, does anyone have experience with this error? How have you worked around it? Thanks so much for your input. UPDATE: This is the code that is actually throwing the error: > var scalebar = new esri.dijit.Scalebar({ map: map, attachTo:"bottom-left" > }); Note that I have included the dojo require statement for that class: > dojo.require("esri.dijit.Scalebar"); And if I simply comment out that function, the same undefined error comes with the next esri object class referenced (which also has its appropriate dojo require statement in the doc).


title: Problem with PostGIS in Opengeo Suite in Ubuntu 13?

text: Everything is working fine after installing the Opengeo Suite in Ubuntu except PostGIS issues. In the dashboard the PostGIS manage is disabled(grey) but PostgreSQL is working via pgadmin3. And also when we create a new database the tables "Geometry Columns" & "Spatial Ref_System" isn't getting created. Also we are not able to import any shape files.


title: Need an Algorithm to Determining Polygon Borders

text: I'm currently writing a program in Java and I'm dealing with polygon objects. I want to create a different color border between touching objects depending on their class (i.e. land, ocean, etc). I need an algorithm for detecting which segments of my ocean border another ocean or land. Right now I'm trying to do this by detecting which object share coordinates. And creating a total count of each class, and whichever class has the highest count, that is what that segment is bordering. It's too simple and doesn't work in about 40% of the cases, so I need something better. Can any one recommend a different idea, or a source to look at? Thanks in advance.


title: QGIS 1.8.0-Lisboa DBmanager (OSx) can't connect to remote PostGIS 2.0.2 database (Debian)

text: The PostGIS database is remote but I have made changes to the Postgresql config files and can access it remotely via pgadmin. When I open DBmanager in QGIS I see a note on the bottom that says "postgis:No module named psycopg2." I can't find anyway to connect to the database through the GUI. I have been able to connect to the database in QGIS using the SPIT plugin, though when I try to import layers it returns the following error: "Problem inserting features from file /path/to/shape/file.shp Invalid table name." I can get the data from QGIS to the server by ftping the files to the server and then importing them with the shp2pgsql command. I can then download the data back into PostGIS without a problem, so I know the database is working properly. What I want though is to be able to import QGIS layers to the PostGIS database from within QGIS. Please forgive my ignorance. I'm learning this stuff as I go. Any clues you can provide will be much appreciated. I've looked at the answer here: psycopg2 error with QGIS 1.5. I might be missing something, but it didn't help me. That answer explains how to config Windows registry keys, but I'm running QGIS on OSX and the database server is running on Debian (stable), so I don't have registry keys.


title: Is it possible to create a data frame graphic using arcpy?

text: I'm interacting with an MXD using arcpy, and wish to create a graphic within a data frame - as with the Draw toolbar when working manually. The GraphicElement page says... > The GraphicElement object provides access to properties that enables its > repositioning on the page layout as well as methods that allow for > duplicating and deleting existing graphic elements. ...which doesn't mention the ability to create a new graphic element. Is there another area I should be looking at in arcpy, or is this not currently supported?


title: Creating Pie Chart with radius based on another field?

text: I have to symbolize my ares polygon feature class and create pie Chart based on the following values in ArcGIS 10.2 for Desktop:               Area    Tds_Calc    Free_Chlorine   Ph    Fluoride  Turbidity   Failures     ========================================================================          Area A  39.87           23.47       5.14     0          2.57     5          Area B  28.54           54.88       10.24   3.41        1.95    19          Area C  72.05            5.12        0.39   0.79        1.57    12          Area D  1.61             61.29       8.06   14.52       1.61    25          Area E  15               66.67       3.33      5        6.67    30          Area F  55.56            27.78       16.67     0         0       9      Pie Chart will be created based on the 5 fields and the last field, i.e. "Failures" should be used as radius of the pie chart. Thus, we will have size of the Pie based on the "Failures" and rest of the fields will be depicting the share of parameters. How can I achieve this task?


title: How to use QGIS from command line on Windows

text: I have QGIS running on a windows operating system. I am willing to use QGIS from command line to follow this procedure. Unfortunately, i cannot manage to have qgis declared in the command windows. I have added "C:\Program Files (x86)\Quantum GIS Lisboa\" to the path variable. but it does not seem to be sufficient. Would you have any tip on that ? **EDIT** As proposed by Andre Joost, declaring the \bin subdirectory is a first step. Following procedure, the command I am trying to run is:                C:\...>qgis --project myproject.qgs --snapshot image.png       --width 1500 --height 1000 --extent 831421.30, 6496828.19, 860046.84, 6536520.86      QGIS loads the project and then the cursor freeze/standby. and Qgis answers no more for several minutes and I kill it. the snapshot image.png and companion world file are generated but are empty (7Ko png). the cmd output is:               _path=C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\Scripts     OSGEO4W home is C:\PROGRA~2\Quantum GIS Lisboa      **END of EDIT** best regard, Pascal


title: Can ArcGIS Online overlay raster data?

text: I'm starting to explore ArcGIS Online as my first introduction to web mapping (perhaps not the best choice of cloud services but thought I'd give it a shot). I was intrigued by their discussion of "story maps" and wanted to create one that overlaid raster imagery that I had on my system. However, I can't figure out how to add raster data. It seems to only accept .shp, .csv etc. I can add rasters if they're converted to layers, but can't seem to actually overlay them on the map. In looking through the help it seems the option is using a map service, the more I dig into this the more complicated it seems. Is it not possible to simply overlay raster data as one would do in ArcMap Desktop?


title: QGIS WMS didn't show georeferenced tiff

text: # I try to reference image and use them via WMS.   1. Chouse project CRS to EPSG:900913    2. Add OpenStreetMap source using gdal wms              <GDAL_WMS>       <Service name="TMS">           <ServerUrl>http://tile.openstreetmap.org/${z}/${x}/${y}.png</ServerUrl>       </Service>       <DataWindow>          <UpperLeftX>-20037508.34</UpperLeftX>          <UpperLeftY>20037508.34</UpperLeftY>          <LowerRightX>20037508.34</LowerRightX>          <LowerRightY>-20037508.34</LowerRightY>          <TileLevel>18</TileLevel>          <TileCountX>1</TileCountX>          <TileCountY>1</TileCountY>          <YOrigin>top</YOrigin>       </DataWindow>       <Projection>EPSG:900913</Projection>       <BlockSizeX>256</BlockSizeX>       <BlockSizeY>256</BlockSizeY>       <BandsCount>3</BandsCount>       <Cache />     </GDAL_WMS>        3. Reference some image using osm background, using EPSG:900913 as target crs   4. Import generated tiffs into project. By default they was imported not in EPSG:900913 but in EPSG:3395   5. Select EPSG:900913 for layer with imported image - after that image appear there it must be.   6. Save project and try to get layer with referenced image via wms (qgis-mapserver was installed according http://hub.qgis.org/projects/quantum-gis/wiki/QGIS_Server_Tutorial tutorial) I didn't see my referenced image. I check layer codes and so on. I try to get osm background (write both layers for wms layers parameter in different order) and get's only osm background. Next thing I try - save gcp's and use gdal VRT               <VRTDataset rasterXSize="1800" rasterYSize="668">         <GCPList Projection="EPSG:900913">         <GCP Id="1" Info="a" Pixel="64.0462" Line="589.346" X="4.17039e+06" Y="7.5036e+06" Z="0.0" />         <GCP Id="2" Info="b" Pixel="64.8845" Line="72.4501" X="4.17039e+06" Y="7.50362e+06" Z="0.0" />         <GCP Id="3" Info="c" Pixel="1746.68" Line="72.6177" X="4.17057e+06" Y="7.50365e+06" Z="0.0" />         <GCP Id="4" Info="d" Pixel="1746.64" Line="589.724" X="17057e+06" Y="7.50362e+06" Z="0.0" />       </GCPList>       <VRTRasterBand dataType="Byte" band="1">         <ColorInterp>Gray</ColorInterp>         <SimpleSource>           <SourceFilename relativeToVRT="0">/home/dkiselev/desktop/123.jpg</SourceFilename>           <SourceBand>1</SourceBand>               </SimpleSource>       </VRTRasterBand>     </VRTDataset>      In this case I have nothing, maybe GCP's didn't applyed to source image. Any way It was a try to find workaround. Here is gdlainfo for resulting tif file:               gdalinfo /home/dkiselev/desktop/123_modified.tif      Driver: GTiff/GeoTIFF     Files: /home/dkiselev/desktop/123_modified.tif     Size is 1920, 556     Coordinate System is:     PROJCS["Google Maps Global Mercator",         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"]],         PROJECTION["Mercator_1SP"],         PARAMETER["central_meridian",0],         PARAMETER["scale_factor",1],         PARAMETER["false_easting",0],         PARAMETER["false_northing",0],         UNIT["metre",1,             AUTHORITY["EPSG","9001"]]]     Origin = (4170378.687966782134026,7503652.570609360933304)     Pixel Size = (0.103459182316438,-0.103459182316438)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       COMPRESSION=DEFLATE       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  ( 4170378.688, 7503652.571) ( 37d27'47.34"E, 55d54' 9.61"N)     Lower Left  ( 4170378.688, 7503595.047) ( 37d27'47.34"E, 55d54' 8.56"N)     Upper Right ( 4170577.330, 7503652.571) ( 37d27'53.76"E, 55d54' 9.61"N)     Lower Right ( 4170577.330, 7503595.047) ( 37d27'53.76"E, 55d54' 8.56"N)     Center      ( 4170478.009, 7503623.809) ( 37d27'50.55"E, 55d54' 9.09"N)     Band 1 Block=1920x1 Type=Byte, ColorInterp=Red       NoData Value=0     Band 2 Block=1920x1 Type=Byte, ColorInterp=Green       NoData Value=0     Band 3 Block=1920x1 Type=Byte, ColorInterp=Blue       NoData Value=0      And qgis generated gdal script for it               gdal_translate -of GTiff -gcp 61.699 588.844 4.17039e+06 7.5036e+06 -gcp 67.0641 71.1088 4.17039e+06 7.50362e+06 -gcp 1750.37 71.1088 4.17057e+06 7.50365e+06 -gcp 1755.74 594.209 4.17057e+06 7.50363e+06 "/home/dkiselev/desktop/123.jpg" "/tmp/123.jpg"     gdalwarp -r near -order 1 -co COMPRESS=DEFLATE -dstalpha "/tmp/123.jpg" "/home/dkiselev/desktop/123_modified.tif"      So I try to use vrt as out format for gdal_translate and gdalwarp with rather same result: I can see correctly referenced image in qgis work window, but get's nothing via wms.


title: Model Builder - iterate through columns in a table

text: This seems like it should be a no-brainer, but I am unable to find any info on how to use model builder to run a certain spatial analysis function on each column in a table. Problem: I have a table to sample locations and at each location there are a series of measures (e.g. temperature, flow rate, number of certain species). For each of those attributes (columns in the table), I was to run the IDW function from the surface tools. This would be easy if I had each attribute in a separate shp/feature class, but not very efficient. I tried Iterate Feature Selection, but was unable to pass the field name as the Z value in the IDW function. EDIT: this is for ArcGIS 10.1 (Basic and Spatial Analyst license) Any help would be greatly appreciated! Thanks!


title: ArcObjects running in Addin slower?

text: I have created a class library which does some geoprocessing. The addin calls a class that is an async process. I have ensured the thread is STA and the arcobjects are thread safe (ie. not passed from the addin). All arcobjects are created within the thread. Because this is a class library I wrapped it in a winforms UI and also as an addin. Both sets of code are exactly the same and test were performed using the exact same data. Both the winforms and addin complete the code with the desired results and no memory leaks are evident. For the case of the addin there is no interaction with the map period at this point and there are no mapping or display elements in the winforms code either. the only UI updates are the update of a progress dialog in both the addin and UI. The addin is using a dockable window (user control ui). The issue I am seeing is when the library is called from the addin, code execution is 5x slower than the same code called through the winforms application. Any ideas on where I might look to see why this is occurring?


title: Postgres Performance with Radius Searches

text: Im porting an application from mysql to postgres, and have several tables with 500k+ rows in each of geo points (lat, lon). In mysql i used a bounding box query with haversine to select rows within a given distance of a center point. i understand i can use either postgis's st_dwithin or the earthdistance module (or i could probably port my mysql query over too). Just wondering if anyone has a recommendation on which method would be best performance wise with these large data sets.


title: What are definitions of watershed segmentation, local maxima, variable window search, etc (about Lidar for vegetation)?

text: I am not familiar about Lidar. I want to measure trees height and number of trees. I guess, watershed = to locate the trees, local maxima = to measure trees height. I used these methods in SAGA GIS. I want to compare with Lidar Analyst. What are methods in Lidar analyst so that it can automatically count trees height, trees position?


title: Gdal tiling error for EPSG: 3857

text: i have lots of data which was coordinated with gdal some programatic way at below:               gdal_translate -of GTiff  -a_ullr %s %s %s %s -a_srs EPSG:3857 %s/%s  %s/%s.tif      then for tiling them all i am building virtual dataset with this:               gdalbuilvrt some.vrt *.tif      for getting its geographic info with gdalinfo i am getting this result:               Coordinate System is:     PROJCS["unnamed",         GEOGCS[,             DATUM["unknown",                 SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],             PRIMEM["Greenwich",0],             UNIT["degree",0.0174532925199433]],         UNIT["metre",1,             AUTHORITY["EPSG","9001"]],         AUTHORITY["EPSG","3857"]]     Origin = (.......,.......)     Pixel Size = (4.777314268465985,-4.777314267784062)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  ( .....          Band 1 Block=2816x1 Type=Byte, ColorInterp=Red     Band 2 Block=2816x1 Type=Byte, ColorInterp=Green     Band 3 Block=2816x1 Type=Byte, ColorInterp=Blue      then with this virtual dataset for gettin geodeitc tile i am calling gdal2tiles.py with profil geodetic but i am getting this result:               ERROR 1: No PROJ.4 translation for source SRS, coordinate     transformation initialization has failed.     Generating Base Tiles:     Generating Overview Tiles:      where i am wrong i havent found yet. but when running gdal2tiles.py for only one tif file it is working and in global mapper it returns true result... how can i tile my image pls i need your helps? thank you


title: How to select a Subset of PolySet data in R?

text: I am using the `PBSmapping` package in R to model data from diva-gis. I can read in the shapefiles and plot them with no issue.               library(PBSmapping)     df = importShapefile("/path/df.shp")     plotPolys(df)      But the problem is I can't figure out how to subset the PolySet data. For example, for Bhutan's administrative areas file, runnning `str(df)` shows:               Classes ‘PolySet’ and 'data.frame': 59294 obs. of  5 variables:      $ PID: int  1 1 1 1 1 1 1 1 1 1 ...      $ SID: int  1 1 1 1 1 1 1 1 1 1 ...      $ POS: int  1 2 3 4 5 6 7 8 9 10 ...      $ X  : num  91 91 91 91 91 ...      $ Y  : num  27.2 27.2 27.2 27.2 27.2 ...      - attr(*, "PolyData")=Classes ‘PolyData’ and 'data.frame': 201 obs. of  19 variables:       ..$ PID       : int  1 2 3 4 5 6 7 8 9 10 ...       ..$ ID_0      : int  35 35 35 35 35 35 35 35 35 35 ...       ..$ ISO       : Factor w/ 1 level "BTN": 1 1 1 1 1 1 1 1 1 1 ...       ..$ NAME_0    : Factor w/ 1 level "Bhutan": 1 1 1 1 1 1 1 1 1 1 ...       ..$ ID_1      : int  486 486 486 486 486 486 486 486 472 472 ...       ..$ NAME_1    : Factor w/ 20 levels "Bumthang","Chhukha",..: 15 15 15 15 15 15 15 15 1 1 ...       ..$ ID_2      : int  9280 9281 9282 9283 9284 9285 9286 9287 9141 9142 ...       ..$ NAME_2    : Factor w/ 201 levels "Aemeri","Athang",..: 8 19 73 126 129 138 163 186 4 29 ...       ..$ VARNAME_2 : Factor w/ 0 levels: NA NA NA NA NA NA NA NA NA NA ...       ..$ NL_NAME_2 : Factor w/ 0 levels: NA NA NA NA NA NA NA NA NA NA ...       ..$ HASC_2    : Factor w/ 0 levels: NA NA NA NA NA NA NA NA NA NA ...       ..$ CC_2      : Factor w/ 0 levels: NA NA NA NA NA NA NA NA NA NA ...       ..$ TYPE_2    : Factor w/ 1 level "Unknown": 1 1 1 1 1 1 1 1 1 1 ...       ..$ ENGTYPE_2 : Factor w/ 1 level "Unknown": 1 1 1 1 1 1 1 1 1 1 ...       ..$ VALIDFR_2 : Factor w/ 1 level "Unknown": 1 1 1 1 1 1 1 1 1 1 ...       ..$ VALIDTO_2 : Factor w/ 1 level "Unknown": 1 1 1 1 1 1 1 1 1 1 ...       ..$ REMARKS_2 : Factor w/ 1 level "Second level is probably 3rd level": 1 1 1 1 1 1 1 1 1 1 ...       ..$ Shape_Leng: num  0.82 0.706 0.433 1.053 0.751 ...       ..$ Shape_Area: num  0.01852 0.01736 0.00877 0.04298 0.01871 ...      - attr(*, "parent.child")= num  1 1 1 1 1 1 1 1 1 1 ...      - attr(*, "shpType")= int 5      - attr(*, "prj")= chr "GEOGCS[\"GCS_WGS_1984\",DATUM[\"D_WGS_1984\",SPHEROID[\"WGS_1984\",6378137.0,298.257223563]],PRIMEM[\"Greenwich\",0.0],UNIT[\"D"| __truncated__      - attr(*, "projection")= chr "LL"      How can I subset out PolySet data where "Name_1" is "Bumthang"? Incidentally, in `maptools`, this can be done by:               library(maptools)     df = readShapeSpatial("/path/df.shp")     df.data = as.data.frame(df) #For understanding the data structure     df.Bumthang = admin.mpt[as.character(admin.mpt$NAME_1)=='Bumthang',]


title: How to use/interpret/visualise a distance matrix with unknown coordinates?

text: How would you go about using a distance matrix where points have no coordinates, i.e. what useful statistics can be made and how to visualise it in some useful way?


title: Hydrologic analysis with sextante-saga in qgis

text: (rough translation): Hi, I want to do hydrologic analysis of a DEM with sextante-saga in qgis. I found the plugin and the needed algorithms in SAGA, enabled them in the sextante toolbox, but when I run them, they error out saying SAGA is not installed (not where sextante is looking for them). My question is therefore where to find these SAGA binaries? Thanks > Salve, dal nome penso che sia italiano come me. Spero possa darmi una > mano... Voglio fare l'analisi idrologica da dem su qgis (1.8) con sextante. > Ho scaricato il plugin e messo in C:\OSGeo4W\apps\qgis\python\plugins .. > l'algoritmo di flow direction (catchment area) e gli altri algoritmi sono in > SAGA, che ho attivato dalla toolbox di sextante tramite interfaccia > qgisdesktop, quando lancio l'alforitmo mi viene detto che devo prima > installare saga, ho capito che devo scaricare a parte _SAGA binaries_ e > metterlo sotto la folder di sextante... Può indicarmi dove posso scaricare > quest'applicazione? In un suo post passato ho visto che lo ha chiamato _msw > SAGA package_! Grazie INFINITE per la dritta! Cordialmente laura guarnieri


title: Calculate 'Rectangle' Coordinates Given 2 Points and width

text: I want to create a rectangular polygon using two points as guides. So let's say a journey starts in Egypt and ends in London, my polygon should have 4 points:   * 10 miles further from London than Egypt is, following the line between them (roughly south in this example).   * Halfway between the two cities but 50 miles at right angles from the line that joins them.   * Like the above point but 50 miles in the other direction.   * 20 miles further from Egypt than London is, following the line between them (roughly north in this example). I'll end up with a rough diamond shaped polygon that would completely contain the straight line journey from Egypt to London. I hope this makes sense; any help for how I can calculate the 4 points is appreciated. ### Rough Map Diagram ![http://i.stack.imgur.com/ESp88.jpg](http://i.stack.imgur.com/hn96h.jpg)


title: What is most accurate coordinate system for calculating areas of polygons?

text: I have multiple polygons, all located offshore around the UK. I am trying to calculate the areas of these polygons in square kilometres, without splitting them by UTM zone. Is it possible to do this accurately? I have been using World Mercator, but there are big distortions at this scale. If the polygons were limited to the east coast, say, I would use UTM31N. However, as they are all over the UK waters, I am unsure what to use. Is there a coordinate system that covers this area? British National Grid does not extend far enough out to sea to be of use.


title: Calculating area of an EPSG:4326 Polygon

text: > **Possible Duplicate:**   >  What is the most accurate coordinate system for calculating areas of > polygons? In our Postgres/Postgis database we store polygons in a column defined as lonlat/EPSG:4326, as it may contain data from anywhere in the world. We have a requirement to calculate the surface area (i.e. in square meters or hectares). The method ST_Area can perform the calculation but returns the answer in the used coordinate systems units. That is, it returns the area as what I believe is "square degrees". As the size of a degree is variable depending on the location in the world, this does not make much sense. One solution is to transform (ST_Transform) the polygon to another crs before calculating the area, perhaps the appropriate UTM32 zone, but that will require us to dynamically determine the zone first. Alternatively we could always transform to the Google Mercator (EPSG:900913), but we fear for the precision. Any inputs?


title: Projecting polygon when calculating its area?

text: I have a broad range of polygons (in vector format and in SpatialPolygonsDataFrame class), some going east-west, others going south- north, while even others span polar regions or very large ranges as mid-to- high/low-to-mid latitudes/longitudes. In addition, as a single polygon can consist of several merged sub-polygons, they do not necessarily have to be aggregated or continuous, but can instead be chopped-up across large ranges. I have until now calculated the area without projecting the data using `gArea {rgeos}` as follows `gArea(SpatialPolygons(spdf@polygons))`. Since I am not interested in visualization of any maps but only comparing polygon area and range sizes, I reckoned that the unit (square degrees vs. square meters) did not matter... but now I am not so sure. So my questions are; `(i)` Calculating polygon areas, would I need to project my data using an e.g. `equal area projection`? If yes, `(ii)` Since I do not care about any map visualizations and my polygons can span very broad ranges making it hard to find a single best projection, can I use the same projection for all polygons e.g. `Lambert azimuthal equal-area projection` when calculating the area? Sorry for asking maybe obvious questions, but having searched the web I only got more unsure and being novice in the area does not help. Any explanations or pointers would be highly appreciated, thanks!


title: Running python 2.6 and 2.7 in ArcGIS 10.0

text: I figure that there must be a way for me to run python 2.6 and 2.7 on the same computer. I've got a program that requires 2.7, but ArcGIS requires 2.6. I'd like to be able to run both programs at the same time. I've installed both and in the geoprocessing settings set the default editor to C:\Python26\ArcGIS10.0\Lib\idlelib.idlepyw. Now when I right click in ArcToolbox, select edit, the IDLE 2.6 editor comes up but the script won't appear like it formerly did. I can view the script in IDLE 2.7 but that doesn't do me any good. -Tom Dilts


title: Convert polyline to polygon Using ArcObjects

text: I have a polyline and I want to convert it to a single polygon feature,I used this code but I get error while I can convert the same shape file using feature to polygon on ArcGIS Where is my mistake?                ESRI.ArcGIS.DataManagementTools.FeatureToPolygon ftp = new ESRI.ArcGIS.DataManagementTools.FeatureToPolygon(path + "lines.shp", path + "polygons.shp");             ftp.attributes = "NO_ATTRIBUTES";             GP.Execute(ftp, null);      The path is set to               string path= System.IO.Path.GetDirectoryName(Application.ExecutablePath) + @"\temp\Conversations\";      and the error that raises is:               Error HRESULT E_FAIL has been returned from a call to a COM component.      I tried to change output path to               "E:/MASTER/Project/Semirom/file gis semirom/New Folder/gisiran.gdb/tester"      But I get the same error..Do you thing I must set Environment workspace or sth like that?I also use `ESRI.ArcGIS.ConversionTools.GPXtoFeatures` and `ESRI.ArcGIS.DataManagementTools.PointsToLine(path + "points.shp", path + "lines.shp");` with the same settings and they works fine..:( Is there any better approache to creat a polygon? Thanks in Advance


title: PGrouting encapsulated queries error [quotes problem]

text: I'm facing a quote problem on a pgrouting query. What I'm thinking is wrong is that my query only works if I use one, then four, then eight quotes on encapsulated queries. Is there a way to simplify this ?               SELECT *     FROM pgr_pointsaspolygon(     'SELECT id::int4, st_x(ways_vertices_pgr.the_geom) as x, st_y(ways_vertices_pgr.the_geom) as y        FROM pgr_drivingDistance(         ''''SELECT gid as id, source, target, ST_Length_Spheroid(the_geom, ''''''''SPHEROID["WGS 84",6378137,298.257223563]'''''''') as cost FROM ways'''',          7,          550,          false,          false)     INNER JOIN ways_vertices_pgr on (ways_vertices_pgr.id = id2)');      Thanks in advance


title: creating isochrones in postgis/osm2po/pgrouting and then saving the isochrones as polygons

text: So I am creating some isochrones in pgrouting based upon a table from osm2po. Now I have a query that will get all of the edges that are `c` driving distance away from a point. My question is, how can I save this query or these individual isochrones into a table of polygons? Perhaps a little background might help, but I want to show how the isochrones change over time as different events occur or as event intensities change. So I wanted to be able to create a video that shows how the contours of the isochrones changes over time. Hence I need to create a table that has each isochrone polygon as well as a timestamp for the date. Any help you can provide in converting the isochrones into polygons--hopefully in postgis--would really help.


title: Why do I get multiple GetFeatureInfo responses from ArcGIS Server WMS when the layers are not visible?

text: I am serving 3 identical layers in ArcGIS Server 10.0 as one wms service. Layer 1 is visible from scales 1:2000,000 to 1:250,010, layer 2 from 1:250,000 to 1:50,010 and the 3rd layer from 1:50,000 to 1:1000. Each Layer has a slightly different symbol to give the effect of a reference scale. (The reason for doing it this way is because our wms services are being integrated into a geoportal which uses other 3rd party GIS Servers which use SLDs to symbolize the data). My problem is that when I click on one of the point in the client application, ArcGIS Server returns responses from all 3 layers even though 2 layers are defined as invisible. Is this a known bug in ArcGIS Server 10.0? Is there a way to combat this?


title: Can I prevent QGIS raster clipper from changing output cell/pixel size?

text: I have performed a DEM extraction using the clipper tool in QGIS 1.9.0 Master, using a mask layer, the original resolution of the DEM is 20x20 but in the resultant DEM the new resolution is 19.98x19.98. Both in meters. Is it normal? How can I prevent QGIS clipper tool to change the pixel size?


title: How to convert a column of numeric IDs to letters?

text: I'm working with qgis 1.8.0. I divided a map into 500 features (small rectangles), and now I have the doubt if there is any way to change the ID of the attributes from numbers to letters? When I did the division of the original map, each of the 500 divisions got one numeric ID (please refer to column 3 of image) and now I'm trying to load that map to STATPLANET PLUS and when working with that program, the ID may not be numerical. Any ideas? Regards!! actual attribute table ![image of attribute table with column 3 showing integer values.](http://i.stack.imgur.com/Hy6m4.png) ![configuration settings used on field calculator](http://i.stack.imgur.com/NxEsE.png)


title: SQL Data: How to draw circles with a radius based on lat/long data in DB

text: I have a . sqlite file with many columns, two of those columns being `lat` and `long`. I've been searching around for a bit now (don't think I want to be using `BuildCircleMbr` from spatialite?) but I'm not entirely too sure how to go about it. I'm thinking it has something to do with `GeomFromText` but can't get the query set up correctly. Final product would be having these circles and display them on a map using MapServer/OpenLayers. Any input is appreciated.


title: How to select duplicate polygons?

text: Although it may seem this has been answered before, there doesn't seem to be a solution. I am trying to simply select duplicate polygons. I don't want to delete them or anything like that. I need them selected for analysis. I'm open to solutions in QGIS or ArcGIS, doesn't matter.


title: Replicas: Synchronize changes for individual feature class (not entire replica)

text: I have two ArcSDE databases, running with a two-way replica.One serves our interactive mapping website which gives certain users the ability to update, add points, edit attributes, etc. Then we have our master which is updated frequently by power users. There are only a few feature classes that are currently being edited, and I would like to synchronize these updates ONLY if these feature classes have been edited. I've developed a python script using the compare features tool, which compares the feature datasets with the master. I'd like to only synchronize changes from these three layers, update the master (so child to parent) and then have the master sync any new changes to the edit database (parent to child). Is there any way to synchronize changes and schema for individual datasets or is there a better way to go about this?               import arcpy         aburg_master_DBO_TownOwnedJan2013 = "Database Connections\\Connection to       gisserver.sde\\aburg_master.DBO.Geocortex\\aburg_master.DBO.TownOwnedJan2013"     Aburg_Edit_DBO_TownOwnedJan2013 = "Database Connections\\Connection to gisserver(2).sde\\Aburg_Edit.DBO.Geocortex\\Aburg_Edit.DBO.TownOwnedJan2013"     compare_txt = "C:\\Test\\compare6.txt"          # Process: Feature Compare     result = arcpy.FeatureCompare_management(aburg_master_DBO_TownOwnedJan2013,   Aburg_Edit_DBO_TownOwnedJan2013, "GlobalID", "ATTRIBUTES_ONLY", "IGNORE_M;IGNORE_Z", "0.001 Unknown", "0", "0", "", "OBJECTID;created_user;created_date;last_edited_user;last_edited_date", "NO_CONTINUE_COMPARE", compare_txt)     changes = result.getOutput(1)     print "Comparision Result is: " + changes     change = str(changes)     if change == 'false':         print "LOL"     else:         print "Nope"


title: ArcGIS & historical versions performance expectations

text: Short background: ArcGIS Server 10.1.1 Enterprise edition using MSSQL for our RDBMS. Currently we utilize versioning only to support replication between 2 databases (edit to DEFAULT). Of late, I'm being asked more and more about the potential for producing 'snapshots' in time. Of course, ESRI supports historical based versions, but I get the sense that this is more for project- level management...not long term management, mostly because for performance reasons, you'll need to be reconciling versions, etc. We also conduct many many edits even on a given day; this DB is constantly updated. Questions therefore are:   1. Am I interpreting the ESRI documentation correctly in terms of historical versions not being long-term data management strategies?   2. Are there other options I should be looking at to support this request for historical data modeling/tracking/archiving?   3. For those using ESRI's archiving: what has your performance been like? How long have you sustained historical versions for? Thanks in advance!


title: Is it possible to do a side by side install of ArcMap 10 & 10.1 using xp mode in windows 7?

text: My organization is running Arcgis 10.1 throughout the office. Staff currently is starting a project in which we need to use FEMA's HAZUS software. The problem is that HAZUS is only compatible with ArcGIS 10 SP 1&2\. I was wondering if anyone here had attempted to do an install of 10.0 in Windows 7's xp mode virtualization. I have never used xp mode before and was wondering if it uses a separate drive partition from windows 7 or if this type of installation scenario is even possible. If anyone has any experience with this type of scenario let me know how it went.


title: building a progress dialog with pythonaddins in ArcMap

text: As a component of my add-in, I compute a bunch of summary statistics when the user clicks a button. This can take a while to process and it can appear as though the program has stopped responding. I would like to add a dialog box which can indicate to the user that the program is still processing. I was interested in using wxPython, but I don't necessarily want to rely on users installing dependencies that aren't included with the standard bundle. Any good ideas? I was considering something like this (except for my concern above):


title: how are range domains validated when editing a feature service in ArcGIS Online?

text: I have a feature class and a range domain is applied to a field. I know we can use the "Validate Features" command for validation. My question is that when the feature class is published to ArcGIS Online and exposed as a feature service, how is validation applied? Is it automatic?


title: PostGIS Raster functions not working

text: Hopefully I have a simple question that I'm not phrasing correctly in google. I'm trying to learn about storing and producing rasters in PostGIS. There are tons of examples and documentation online to learn from, but whenever I attempt to use _ST_AsRaster_ or _ST_AsPNG_ on a vector geometry table, I get the error: > ERROR: function st_asraster(geometry) does not exist > > LINE 1: SELECT ST_AsRaster(public.cspd_patrol_grid.geom) > > HINT: No function matches the given name and argument types. You might need > to add explicit type casts. I wasn't sure if I had a recent enough version of PostGIS with 2.0.3, so I upgraded both PostgreSQL and PostGIS to the latest table versions (9.3.5 and 2.1.3 respectively), but it did not resolve the problem. I suspected it had something to do with this: http://www.postgis.net/docs/manual-dev/postgis_gdal_enabled_drivers.html where the GDAL drivers are disabled by default, but going through the PostGIS installation, it has prompts asking if you want to enable them. I do not know how to verify the settings of postgis environmental variables. So does anyone have any idea what I'm missing from a vanilla install of postgres&postgis to be able to use raster functions?


title: Java how to create geotiff from points

text: I have an array of objects where each object has three fields: lat, lon, and category. It is essentially pixels that have been given a land use category. I cannot figure out how to turn this data into a geotiff. I assume I should use geotools. The coordinates are in dec degrees lat long. I know I can use q or arcgis to do this but I need to do it as part of an automated process. Just need a push in the right direction.


title: How to display MousePosition outside map container?

text: I am a new to OpenLayers. By Default the MousePosition displays the coordinates within the map. I need to display the coordinates outside the map is a specific place. How do I do this?


title: MapXtreme Java Edition

text: I'm learning GIS in Java with MapXtreme Java Edition but I can't find MapXtreme Java Edition file setup for Window (trial). Can you share link download or software for me?


title: How does Jenks optimisation method assist in ESRI mapping?

text: Why would someone use the Jenks optimisation method in ESRI mapping? What would be easily located or discovered with this analysis?


title: Changing the CRS of a shapefile

text: Im trying to convert shapefile from one CRS to another in QGIS. On the layer menu, I right click on the shapefile, save it as with other name and there I change the CRS to the desired one (WGS 84). Then I check the properties of the file and under the tab of METADATA and it appers under Layer Spatial Reference System the desired CRS but when I project it, its not on the desired CRS. I also can see on the lower central coordinate display on the main screen that its not in WGS84 (because of the values). Any idea whats happening? Thanks. ![enter image description here](http://i.stack.imgur.com/qhLK0.png) ![enter image description here](http://i.stack.imgur.com/iN3gk.png)


title: Move points along a line, 100m left or right

text: In ArcMap 10, I have several hundred points taken along a zig-zag transect (line) but I need to move these points 100m to the left or right, perpendicular to the line. The points are stored as a single shape file with a 'SIDE' field, which is either 'left' or 'right' indicating which side of the line to move the point. I have tried linear referencing but got totally lost. Is that the tool to use or is there another simple method of achieving this?


title: Converting .fdshape to .shp or other QGIS readable format

text: I have been searching for a solution to being able to read .fdshape extension files in QGIS with no luck. It seems that the data had to be converted from .shp to .fdshape in the first place so I was wondering if there is an easy way of getting this information converted back?


title: What does "Computed -srcwin falls outside raster size" mean when clipping?

text: What does this mean, when I try to clip and save the result as a sqlite or other format please? "Computed -srcwin falls outside raster size of 7807x4771" cheers, Bruce


title: Determine when a route enters and leaves a polygon

text: I have all of the segments of a route from White Horse, YK to Miami, FL. Now I want to determine how many miles I travel in each State / Province. I have Polygons for each State and have been able to use Ray Casting to test whether or not each Segment is in a State but is there a faster way to determine when a route enters and leaves each State without testing ever single segment? Doing each segment would require over one hundred thousand lookups.


title: ArcObjects: How to get the Cell Size from the raster.MeanCellSize call?

text: When I setup the call to MeanCellSize -- for which ESRI documentation states will return the mean average of the cell size -- all I get is an IPnt object. I've played with a ton of the various casts but cannot figure out what to cast it to so that I can actually get a value for the mean cell size!! Has anyone figured this out?


title: How to create an editable custom popup or infoTemplate

text: Does anyone have a resource with an example of how to create a custom popup/infoTemplate in JSAPI? I found this link: close example in .net but I want the user be able to just add text to a popup created after clicking on the map. I know how to create an infoTemplate, etc. but not one that's editable. That would be the perfect scenario as it shows an anchor location. Big thanks in advance.


title: Points from python directly to ArcMap

text: I have a python code that is generating coordinates one by one. Is there a way to display them directly in ArcMap as they are generated? I am aware of the insert cursor but the lock on the table prevents python from writing into it. **Ok, now that I used the Python window in ArcMap new error has occured. Could it be because the "czk" is a list?** Here is the code:               czk=[413960.0, 47130.0]     cursor = arcpy.da.InsertCursor("C:/GIS/AlJeba/Tocke/Tocke.shp", ("SHAPE@XY"))     for row in czk:         cursor.insertRow(row)     del cursor      And here is the error: Runtime error Traceback (most recent call last): File "", line 333, in TypeError: argument must be sequence of values


title: How to map out commuting patterns by county?

text: I am looking to create a map that shows the commuting patterns in between counties. I am using data from Table 1 at this **link**. Obviously it will be formatted in a different way. Any ideas? Also, here is a picture of a previous map done that is achieving the same effect![enter image description here](http://i.stack.imgur.com/EIB5e.jpg)


title: extend the size (area) of raster

text: My problem is that I would like to multiply two rasters using raster calculation tool. Both rasters have values 0 and 1, but the covered area of first raster is bigger (100 x 100 miles) then the area of second (80 x 80 miles). How can I extend the second raster to the same size? I would just like to add columns and rows with 0. At first I thought that this would be simple, but after searching on internet I didn't found any solution. :( Thank you in advance


title: Merge String fields in ModelBuilder by undefined ArcGIS user?

text: I'm using ArcGIS 10. I have a doubt, because I already have a merged shapefile, with two string Fields, some Null and some filled (comes from the merge tool). So, now I want to get the values from the two string fields into a NEW string Field but with ModelBuilder. So, I do know how to do it manually and get these values with the FIELD CALCULATOR, because I could choose the field names, but **how could I get this prepared for an user with ModelBuilder?** VB Code:               Dim X     if IsNull( [ETIQUETA] ) then        X = [Data_Hamlet_Galicia_POI_NAME]     else        X = [ETIQUETA]     end if     __esri_field_calculator_splitter__     X      The FIELDs above may not have the same name, so it will not be right to let this code into the "Field Calculator" Tool within ModelBuilder. Any suggestions? Here is my model. ![Model](http://i.stack.imgur.com/IjfbD.png) \-------- E D I T ------ Well, @Hornbydd , I've tried but fails, maybe I'm doing something wrong. Here the new model and the error log: ![Model2](http://i.stack.imgur.com/EdTy6.png)


title: Is it possible to decrease detail level (decrease file size) in Spatialite OSM export?

text: I currently use Spatialite exports from OSM data (together with Google maps like styles by Underdark) which are very useful for me as the files are easily portable with the project and I do not need Internet connection. However, sometimes the result is a huge file - for example, my current area:               ogr2ogr -f "SQLite" -dsco SPATIALITE=YES -spat 14.20622 50.57070 17.32913 48.81776 Jeseniky.db czech-republic-latest.osm.pbf      results in about 1.7 GB file which loads slowly even on powerful machines. In fact I do not need so much detailed map data for common use. Is it possible to create similar Spatialite file from OSM data, but with data reduced to equivalent of 1:100 000 or 1:50 000 scale? It would be perfect and I suppose I could have whole country in such detail level for common use and sometimes only would export "full detail" area as I currently do. is it possible? thanks


title: How can I replace all null values from attribute table to zeros?

text: I have a new column in my attribute table with a lot of null values and I want to replace all them to zeros. How can I do that in calculator? I'm using QGis 1.8 Thank you!


title: Sorting points to flight path route

text: Original question: I would like to sort my points in QGIS. See image Restated question: The grid tool currently numbers by starting at R1C1 and going to R1Cx, then R2C1 to R2Cx. How do I number the points from R1C1 to R1Cx, then R2Cx to R2C1? Can the solution be applied after a regular grid has been destroyed (points removed) or does it need to be applied to a whole grid, assuming no option for this ordering scheme exists in the grid creation tool itself? The point is to avoid the long 'row reset' or 'back to start of row' line as this wastes time and fuel, rather than moving directly down to the next row. ![enter image description here](http://i.stack.imgur.com/poXQO.jpg) ![enter image description here](http://i.stack.imgur.com/MrnTO.jpg)


title: How does geoserver determine what SRID to pass to PostGIS [OPERATION ON MIXED SRID GEOMETRIES]

text: We are using Geoserver to render some data onto maps in a web application. The map shapes are read by Geoserver from a postgis database. The shapes are imported into postgis using the ogr2ogr utility. An event on the map makes geoserver query the postgis database for data specific the area clicked. Our application does not specify what query geoserver should run against the database. What is happening is that, for some reason we do not know yet, geoserver fires the correct query when postgis is populated with one set of shapefiles but fires a wrong query when postgis is populated with another set of shapefiles. We know that while geoserver is able to recognise the projection of the first set of Shapefiles as EPSG:900913 when forming the postgis query, it is not able to recognise the projection of the second set. As a result, the query throws an "Operation on mixed SRIDs" postgis error when run against the second set. We compared the two sets of shapefiles and found that the .prj file of the first set of shapefiles specifies many more parameters than the same file in the second set of shapefiles, important of which is the PROJECTION["Transverse_Mercator"] . Here is the query geoserver fires.               SELECT "id","district","poll_id","deployment_id","top_category","description","ogc_fid",encode(ST_AsEWKB("wkb_geometry"),'base64') as "wkb_geometry","dname_2010" FROM (select geoserver_pollcategorydata.*, ssd.ogc_fid as ogc_fid, ssd.wkb_geometry as wkb_geometry, ssd.name_3 as dname_2010 from southsudan_districts2013 as ssd join geoserver_pollcategorydata on geoserver_pollcategorydata.district=ssd.name_3 where geoserver_pollcategorydata.poll_id = 231 and geoserver_pollcategorydata.deployment_id = 1) as "vtable" WHERE  "wkb_geometry" && ST_GeomFromWKB($1, -1) AND ST_Intersects("wkb_geometry", ST_GeomFromWKB($2, -1)) LIMIT 1      Any help figuring out what is causing Geoserver to query with an SRID -1 (ST_GeomFromWKB($2, -1)) on the second set of files will be highly appreciated. Thanks


title: QGIS - move a point object to specific coordinate

text: Is there a way of select an object (point, polygon) and move it to a specific coordinate (manually inserted)? (with Numerical Digitize plugin can not make a move of the object but rather it creates me a new object) thank you,


title: Trouble joining table in arcgis. Joining field is 'greyed out'

text: I'm trying to join a shape file with a excel file I converted into a CSV. I had done this before with the same shape file but with a different CSV but this time when I try to select the column it's removed from the join list. I have tried exporting the CSV as DBF then turning it back into a CSV as well which did not work. I have included two screenshots below. I want to use the column PCODE3 because it shares the same data as MUN_P_6. Screenshots: Full list of fields: http://i.imgur.com/lhEHTrk.png Limited lis of fields: http://i.imgur.com/G2WI2w0.png


title: How to set the output extents of rasterization in QGIS

text: I am trying to rasterize a bunch of vector based shape files in QGIS. All these files are covering certain types of ground use in a city (houses, street, open field,...). I need these files to be rasterized because I will be running an algorithm on them to determine the least cost path. And for this I need them to be rasterized with a certain fixed pixel size (1m x 1m in the real world). My vector files are 6 km x 10 km but when I rasterize one of the files (for ex: railroads) which isn't reaching over the total 6 km in x direction and over the 10 km in y direction but only over the first 3 km x and 10 km y, the rasterized file is also only 3 km x 10 km. Somebody any idea how to solve this? Because in this manner, when I choose a pixel size of 6000 + 10000 I'll get pixels with a size of 0,5 m x 1 m in stead of 1 m x 1 m.


title: How to plot coordinates from *.csv file on Google Maps?

text:
I want to highlight a mobile trace on Google Maps and GIS data from my network must also be visible. I am able to plot hard coded coordinates on Google Maps but not from a *.csv file. Note: *.csv file contains Lat & Long values to highlight traced route Hi....I am able to highlight the route now but if i enter too many coordinates then the page responses very slowly....kindly can any ony let me know the limits ....I am post my code...plz mention if i am going in wrong direction.....Thanks               <>                  <!DOCTYPE html>                         <meta name="viewport" content="initial-scale=1.0, user-scalable=no">          <meta charset="utf-8">     <title>MySiteMap</title>
 This is My web Page

Data in Excel is : **3652 https://sites.google.com/site/turbhegiskmlfiles/gis/shirur_to_ahmednagar.kml **

 Data is as follows:
                                                 <>


title: ArcGIS Android Reverse geocoding, Cannot perform query. Invalid query parameters

text: I'm trying to write reverse geocoding. And i copied code given here... Here is my code sample..                     public class MagProActivity extends Activity {                      private MapView mMapView;                 private GraphicsLayer locationLayer;                            /** Called when the activity is first created. */                 @Override                 public void onCreate(Bundle savedInstanceState) {                     super.onCreate(savedInstanceState);                     requestWindowFeature(Window.FEATURE_NO_TITLE);                     setContentView(R.layout.main);                          mMapView = (MapView) findViewById(R.id.map);                     mMapView.addLayer(new ArcGISDynamicMapServiceLayer("http://services.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer"));          mMapView.setOnLongPressListener(new OnLongPressListener() {                      private static final long serialVersionUID = 1L;                      @Override                 public void onLongPress(float x, float y) {                     Point mapPoint = mMapView.toMapPoint(x, y);                          new ReverseGeocoding(MagProActivity.this, mMapView).execute(mapPoint);                      }             });                }         }      ReverseGeocoding class               import java.util.Map;          import android.app.Activity;     import android.os.AsyncTask;     import android.widget.Toast;          import com.esri.android.map.MapView;     import com.esri.core.geometry.Point;     import com.esri.core.geometry.SpatialReference;     import com.esri.core.tasks.ags.geocode.Locator;     import com.esri.core.tasks.ags.geocode.LocatorReverseGeocodeResult;          public class ReverseGeocoding extends AsyncTask<Point, Void, LocatorReverseGeocodeResult>{              private Activity curActivity;         private MapView mapView;              public ReverseGeocoding(Activity _curActivity, MapView _mapView){             this.curActivity = _curActivity;             this.mapView = _mapView;         }              @Override         protected LocatorReverseGeocodeResult doInBackground(Point... params) {             // create results object and set to null             LocatorReverseGeocodeResult result = null;             // set the geocode service             Locator locator = new Locator();             try {                      // Attempt to reverse geocode the point.                 // Our input and output spatial reference will be the same as the map.                 SpatialReference mapRef = mapView.getSpatialReference();                 result = locator.reverseGeocode(params[0], 50.0, mapRef, mapRef);             } catch (Exception e) {                 e.printStackTrace();             }             // return the resulting point(s)             return result;         }              protected void onPostExecute(LocatorReverseGeocodeResult result) {                  String resultAddress = "";                  // Construct a nicely formatted address from the results             StringBuilder address = new StringBuilder();             if (result != null && result.getAddressFields() != null) {                 Map addressFields = result.getAddressFields();                 address.append(String.format("%s\n%s, %s %s",                         addressFields.get("Street"), addressFields.get("City"),                         addressFields.get("State"), addressFields.get("ZIP")));                      // Show the results of the reverse geocoding in a toast.                 resultAddress = address.toString();                 Toast.makeText(curActivity, resultAddress, Toast.LENGTH_LONG).show();             }         }          }      When i long click on mapview this error will come com.esri.core.io.EsriServiceException: Cannot perform query. Invalid query parameters and toast will show all parameters null except city.. I think i'm doing some mess with Tiled service on Dynamic service.. Please help what i did wrong.. Thanks in Advance


title: Rasterization in different bands

text: I have a set of vector based shape files, representing different types of ground coverage in a city (for ex: roads, buildings, open field,...). For the algorithm that I want to run on it, I need just one raster file with these different ground coverages represented by different bands. The final file needs to be from IMG file type. Can somebody help me with this?


title: Very poor performance of QGIS and PostGIS2. User error or feature?

text: I have a large table in Postgis2. It has >500 million polygon, line and point features all happily living in there and the table performs brilliantly in other GIS. I am now trying to just add that table to a project in QGIS but here is the thing, why does it take so long? When i try to 'Add PostGIS overlay' it says detecting for about 4 hours to tell there are 3 geometry types. OK, i can take that if i only have to do it once... Once it finishes that i select all three parts click ok.. so far im 2 days into trying to get it to load into the canvas. What is going on? I am zoomed in so only a few polygons should be rendered on screen. Is this an error or a 'feature' of QGIS? Is there something i need to specifically do to enhance the performance for QGIS? Many thanks for your help people. Update: Thanks for the help.   I am using MapInfo and Cadcorp SIS both of which love the "everything in one place" model rather than the splitting it up, join it back later with queries model. And thanks @AndrewJoost, the estimated metadata saved about an hour but still it takes a few hours too long and misses one of the geometry types off the list... QGIS not supporting mixed geometries in one overlay is no problem, as long as it supports them in one table.   And the way it does this with a filter is great for small tables but just fails when you have a large one it seems. I am unable to seperate the table into three parts as this does not make sense for the other systems using the data, and would also break them. @NathanW How can i get 3 records for a single table into the geometry_columns view now that it is no longer a physical table? any ideas? Sorry about all the newbie questions Any other ideas to stop QGIS trying to be helpful and scanning a table before allowing me to add it? I dont mind choosing geometry types by hand for every overlay if thats what it takes? Cheers


title: Can a joined-layer have its attribute changed depending on what layers are selected/deselected?

text: Using QGIS, I have a polygon-grid layer which covers everything, **Layer_All** joined with 2 other grid layers: **Layer_1** and **Layer_2** , both of which intersect at various points. I would like **Layer_All** to show a different colour if **Layer_1** & **Layer_2** are either selected or both selected with the intersection being another new colour etc. Is this level of interactivity possible via the Field Calculator? EDIT: Should probably add a bit more info. Each layer has a common attribute column with data defining what colour that polygon should be. I'm looking at something which can add up the data if another layer is maunally selected. Hope this helps.


title: How do I link an existing install of GDAL to FileGDB in Linux?

text: The distribution I use (CentOS) only has GDAL 1.8 available through elgis. I've downloaded the FileGDB API from ESRI and I'd like to be able to link it after an install. How can I do this?


title: Using QGIS with EPSG:404000

text: I'm using QGIS with a flat-world approximation dataset and for my specific use, geometric coordinates on a 1 unit = 1 metre scale are perfect. Unfortunately, QGIS doesn't have EPSG:404000 built-in, and complains if I try to create a custom CRS with +proj="" as the instructions on http://spatialreference.org/ref/sr-org/7589/ say. Is there any way to teach QGIS about EPSG:404000 and get it to work basically without projection ?


title: Is there any open-source object-oriented GIS?

text: Most GIS are layer based. But I am looking for open-source object-oriented GIS. Is there any out there? By object-oriented I mean logical data model for vector GIS. Here's explanation of my definitions.


title: What is the difference between the astar_sp and shortest_path_astar  functions?

text: I want to ask about what different of function : "astar_sp" and "shortest_path_astar" ? Do you think the query's result from both of them is same or different?


title: How to prevent user from deleting a feature?

text: How can I prevent the user from deleting a feature?               void m_editEvents_OnDeleteFeature(IObject obj)             {                 // Code             }


title: Quantum GIS Garmin Maps62 device and GPS utility - making a polygon for random point

text: I have a Garmin GPSMaps62 unit and GPS utility. I have logged tracks and need to get these to load as polygon layers in Quantum so that I can then ask Q to give me random sample points and upload these back to the unit and use the find option. Any help in simple steps (I am learning) would be fantastic please. I am particularly struggling with getting the tracks downloaded to quantum via GPS utility to show as polygon layers.


title: In Arcpy, Erase function error: "XML being loaded could not be parsed"

text: Here we go! A error in my previously-smooth-running arcpy script has me confounded. The error occurs where the script selects individual features in two different feature classes, and then performs an Erase function. It seems to be working fine for the first one or two files that it erases, but then stops with "Error 000210: Cannot create output." Using arcpy.GetMessages gives greater detail... > ERROR 000210: Cannot create output The XML being loaded could not be parsed. > [GDB_Items] The XML being loaded could not be parsed. [XML document must > have a top level element. XML document must have a top level element. The only useful info online I've found suggested restarting the script each time it fails, which works - in that I restart, and the script moves on to the next item in line to be processed . This is (most of) the offending script, and I've highlight where the Error is coming from:               import arcpy     from arcpy import env     import time     from time import strftime          #Set overwrite output     arcpy.env.overwriteOutput = True     print arcpy.env.overwriteOutput          donutlog = "C:\\BOEM\\Donuts\\donut_log.txt"     workspace = "C:\\BOEM\\pt_buff2\\"     fileList = []     buffersd4 = []     buffersd3 = []          buffs = workspace+"buffer.gdb"     directorypath = buffs               #Insert test to see if ID was done around here:     finishedIDs = {}          FClist = []     with open(donutlog) as f:         for line in f:             if line.startswith("Executing"):                 FCname=line.split(' ')                 FClist.append(FCname[2][:-4])          for dirpath, dirnames, fgdb in arcpy.da.Walk(directorypath,                                                   datatype="FeatureClass",                                                   type="Polygon"):         for fc in fgdb:             fc_split=fc.split('_')             StartTime=time.clock()             if fc.endswith("d4"):                 if fc in FClist:                     print fc + " has already been looked at!"                 else:                          #list features by APR16_ID in that fc                     APR16_ID4 = [row[0] for row in arcpy.da.SearchCursor(buffs+"\\"+fc,['APR16_ID'])]                          # make feat layer of fc_endswith_d4                     in_feat4 = buffs+"\\"+fc                     out_layer4 = fc+"_lyr"                     arcpy.MakeFeatureLayer_management(in_feat4,out_layer4)                     print "Make Feature Layer for " + fc                          #make feat layer of fc_endswith_d3                     in_feat3 = buffs +"\\"+ fc[:-2]+"d3"                     out_layer3 = fc[:-2]+"d3"                     arcpy.MakeFeatureLayer_management(in_feat3, out_layer3)                     print "Make Feature Layer for " + in_feat3     ### ERROR IN THIS SECTION ####                     for ID4 in APR16_ID4:                         where_clause4=""" "APR16_ID" = """ + "'"+ID4+"'"                         arcpy.SelectLayerByAttribute_management(out_layer4,"NEW_SELECTION",where_clause4)                         arcpy.SelectLayerByAttribute_management(out_layer3,"NEW_SELECTION",where_clause4)                         erase_out= "C:\\BOEM\\Donuts\\Donuts_d4_output.gdb\\"+"don_d4_"+ID4                         ### Error in the Erase function here:                          arcpy.Erase_analysis(out_layer4, out_layer3,erase_out,"")                         print "Erase "+ out_layer3 +" from "+out_layer4+" to produce "+erase_out                         print arcpy.GetMessages()                         print "___"      Thanks for any and all ideas.


title: How to use na.Solve() on layers already created?

text: I am fairly new to arcpy and I am having troubles recreating the following set of instructions with a python script.   1. Open the KeyChainTest.mxd   2. Enable the Network Analyst extension.   3. Click solve on the Network Analyst toolbar. The green route line should go through the large point, which represents a locked door. If it avoids the door by going the long way around the block, the test has failed.   4. Double-click the Route layer in the TOC to open the Layer Properties. Click the Attribute Paramters tab. For the KeyIDs parameter, remove to essentially indicate you don't have the key for the door anymore.   5. Click OK to close the Layer Properties dialog box. Solve again.   6. The route should should take the long way around to avoid the door. If it goes through the door, the test has failed. I have all the data setup, I just don't understand how I would go about solving with arcpy and then check whether the solve worked correctly. Thanks for your time! (If I need to narrow down something or add additional information don't hesitate to ask)


title: How to access ArcGISDynamicMapServiceLayer with credentials (token)

text: I know I simply need to append a token=XXX to the export query string but how, in code, do we initialize "esri/layers/ArcGISDynamicMapServiceLayer" to use credentials? In my case I only want to provide a token value. Here's my code to show how I've hacked around this...probably good enough for my needs but seems wrong. Should I be creating a Credentials instance and associating it with the layer?               function createDynamicLayer(options, map) {         var dataset = options.dataset;         var query = ioQuery.queryToObject(acme.right(dataset.CacheURL, "?"));         var imageParameters: esri.layers.ImageParameters = new esriImageParameters();         var credentials: esri.Credential;              imageParameters.format = query.format || "png";              var url = acme.left(dataset.CacheURL, "/export");         var layerName = acme.left(acme.right(dataset.CacheURL, "/export/", true), "?");         if (layerName) {             layerName = layerName.split(",");             layerName = array.map(layerName, "return parseInt(item, 10)");             lang.mixin(imageParameters, {                 layerIds: layerName,                 layerOption: esriImageParameters.LAYER_OPTION_SHOW             });         }              if (query.token) {             credentials = new IdentityManagerBase.Credential({ token: query.token });             url += "?" + ioQuery.objectToQuery({ token: query.token });         }              var layerOptions = lang.mixin({}, lang.getObject("layers.esri.dynamic", esriConfig), {             id: options.id,             imageParameters: imageParameters         });              //Takes a URL to a non cached map service.         var dynamicMapServiceLayer = new ArcGISDynamicMapServiceLayer(url, layerOptions);         return dynamicMapServiceLayer;     }


title: Why can't QGIS connect to MS SQL Server?

text: All, Up until yesterday, QGIS connected to my SQL Server 2012 database no problem. However, yesterday I installed SQL Server Management Studio (SSMS), and the SQL Server Client tools, and now I get this error:![enter image description here](http://i.stack.imgur.com/pQcEl.png)   * 'test connection' succeeds   * I am connecting to this account with a valid User   * I am connecting successfully through SSMS   * I am connecting successfully through ArcMap _**Edit: This seems to be a 32bit / 64bit issue, as I can connect to our 32bit SQL Server 2012 Any thoughts???_** Thanks! -m


title: Count the number of characters in a field (Python)

text: I have a field that contains street names and a second field called COUNT. I need to perform a Python field calculation on the COUNT field that counts the number of characters in the STREET_NAME field. For example if the street name is Main St then the COUNT field would be 7.


title: Raster Clipping and Manipulation

text: I have two raster tiff files one is a raster that I created from 2 Foot contour lines, the other raster file is a small area that was surveyed to half a foot accuracy. My question: Is there away to essentially clip the area out of the 2ft raster where my survey raster resides and then mosaic or merge the survey raster with the 2ft raster?


title: SQL Server geography.Intersects Yields Unexpected Results (arch shape)

text: I've coupled ESRI's ArcGIS with SQL Server 2008 R2 in an attempt to allow users to search for geographical points using a polygon (with ArcGIS' drawing toolbar). I've run up against an issue that's not making any sense to me. It seems easiest to illustrate it visually:   1. I draw a rectangle along the top of the northern hemisphere and get results I wouldn't expect (I know there are more points up there): ![result 1 with a rectangle](http://i.stack.imgur.com/ljROb.png)   2. I draw a smaller square within that rectangle and get the results I was expecting: ![result 2 with a square](http://i.stack.imgur.com/FAPbr.png) Can anyone explain why this might be happening and if there's some way around it? It seems like a SRID mismatch, but I'm using 4326 for ArcGIS and for SQL Server. If I switch to using `geometry` instead of `geography` the rectangle works fine. I assume this is because geometry isn't taking into account the curvature of the earth, which isn't present on the flattened-out map. `geometry` has problems of its own though (my map wraps infinitely left/right, and so some of my polygons refer to points outside 180 degrees, which `geometry` doesn't seem to support). For reference, here are the results in SQL Server: For the rectangle: ![enter image description here](http://i.stack.imgur.com/dBHEE.png) For the square: ![enter image description here](http://i.stack.imgur.com/t6EKI.png) And `UnitsSearchByLocation` looks like this:               ALTER PROCEDURE [dbo].[proc_UnitsSearchByLocation]         @polyText varchar(4096),         @wkid int     AS         declare @g geometry         set @g = @polyText              declare @poly geography              begin try             -- try the polygon as provided, although it has the potential             -- to be oriented incorrectly             set @poly = geography::STPolyFromText(@polyText, 4326)             select * from Unit where Location.STIntersects(@poly) = 1         end try         begin catch             -- assume the error was a result of the orientation being wrong             -- and attempt to correct it             declare @validPolyText varchar(4096)             set @validPolyText = @g.MakeValid().STUnion(@g.STStartPoint()).ToString()             set @poly = geography::STPolyFromText(@validPolyText, @wkid)             select * from Unit where Location.STIntersects(@poly) = 1         end catch     GO      Although it seems irrelevant to me, I'm using this Imagery from ESRI (: \- https://server.arcgisonline.com/ArcGIS/rest/services/ESRI_Imagery_World_2D/MapServer Thanks for any help you can provide..


title: gdaltools in qgis v2.4?

text: Where's the GDALTools plugin for QGIS 2.4? It's not listed in list of available plugins, and it's seemingly not installed out of the box. Is it just that it's not yet compatible?


title: Ortho-rectification of Aerial Image

text: I took an aerial image of our farm using RC quadcopter. How am I going to orthorectify this image? Please can you help me with detailed instructions/methods involving the process.


title: Convert series of WKT to Shapefile

text: I have seen am ability in `SQL Server 2008 Spatial Tools` here that saves result of query as shapefile..I am intrested in how to do this task using C#..Is there any way to do this?I found ways to convert WKT to WKB and so on but I could not figur out how to convert them to .shp.I found a tool called `ogr2ogr` but I dont know how to use it in my C# code


title: Is there a JavaScript SDK for Vehicle Tracking?

text: Does exist a JavaScript SDK for Vehicle Tracking? Like this one:Map Suite


title: Adding attachments to features

text: Does anyone have an elegant solution to adding/viewing **attachments** , like a table/photos for example from within ArcMap using only the ArcGIS basic license level? Or any suggestions? I'm looking for a way to deliver an .mxd or file geodatabase with features that are clickable and link to tables of info and photos. The above feature of 10.0 is great; however I am limited by the license we have. Thanks!


title: How to use gdal utilities to subset from a raster?

text: How to you use gdal to subset a smaller raster from a bigger one given the bounding coordinates? EDIT: Can you subset a region from a vrt ?


title: Not able to create a geoprocessing package in ArcMap using my own Network Dataset and associated data

text: I went through the following two examples to get an idea of how to create a custom, local Routing Task . . . This is a walkthrough that takes you step by step through the creation of your first routing geoprocessing package: http://resources.arcgis.com/en/help/...00000m5000000/ And here is the second part, that completes the picture and shows you how to run it from code: http://resources.arcgis.com/en/help/...00000m6000000/ I am able to create the gpk using the above mentioned link instructions. However, when I attempted to create a gpk using my own data, I got an error message each time (General function failure ERROR 999999: Error executing function). When I researched this error, I didn't find anything helpful. I am using a map document that contains a Network Dataset created from a Shapefile Feature class. The settings that I use to create the package are . . . Support ArcGIS Runtime enabled (other two options are disabled) Summary and Tags detail filled in Is there any way to find out more info on the error? I don't know why it's failing, but I keep changing little details and trying it again. Thanks, Renee


title: Render specific vector data from database with openlayer

text: I want to render some data on a map. They are stored as points in my database (postGiS), each one is associated with a name. I want to render points with the same place name in the same color. However, I don't want to put the name of the places directly in my .xml style file. (Too many places) Do you know how can I do ? Maybe an idea is to try to return json data in OpenLayers js, but I don't know how I can get it. Does GeoServer enable to return json to OpenLayers ? I am totaly lost. Please help !!! Thks


title: How to create a point along a line given a distance from another point

text: I'm trying to figure out a point(D) on B-C that is a specific distance from A.   I really don't have any idea on how to proceed with the ellipsoidal math so I figured I'd ask here. How would I figure out the longitude, latitude of D which is 5mi(26400ft) away from A, on the line B-C? The bearing of A-D would be just as useful as i could figure out the point 5mi away from A with the bearing. I wouldn't mind if the formula lined up with the Vincenty model of the earth. Using the Vincenty formula's I've already calculated the following. Points (longitude, latitude)   **A** : -86.32899412316736,41.12507719625815   **B** : -86.29237110757369,41.16845596588609   **C** : -86.19689279712584,41.16923568742915   A-B Distance: 18750ft   B-C Distance: 26289ft   A-B Initial Bearing: 32.5331889647   B-C Initial Bearing: 89.3493884487   ![enter image description here](http://i.stack.imgur.com/ghjab.png) Thanks for the help!


title: Pyqgis vector layer rollback and commitChanges

text: I have problems with the vector layer rollback and commitChanges function, it doesn't matter wich one i call the vector layer is is Changed.Here is my code i hope someone can help me. This is my firt post and im new with qgis. from qgis.core import QgsApplication,QgsVectorLayer,QgsVectorDataProvider from PyQt4.QtCore import QVariant if **name** == " **main** ":               QgsApplication.setPrefixPath('/usr', True)     QgsApplication.initQgis()           vLayer = QgsVectorLayer(                 "/home/diademuertos/scripts/qgis/data/point.shp",                 "vLayerTest",                 "ogr")          vLayer.startEditing()          caps = vLayer.dataProvider().capabilities()          if caps & QgsVectorDataProvider.ChangeAttributeValues:         attrs = { 0: QVariant(2),1: QVariant("B"),2:QVariant("Dont!!!")}         vLayer.dataProvider().changeAttributeValues({ 1 : attrs })     else:         print "OOPS"          r = False     if r:         print "commitChanges"         vLayer.commitChanges()     else:         print "rollback"         vLayer.rollBack()          QgsApplication.exitQgis()


title: Updating file geo database table

text: I have an input csv file that I process using a python script and generate a new output csv file. I then, in the same python script use Table to Table (Conversion) in order to store the data from the output csv file in a file geo database. The information in the input csv file will change but I'm currently only able to create a new and not update the table. So the question is, how can I update the table or if there's a smarter solution to this problem than generating and converting csv files?


title: Mapbuilder can't see nvarchar2 to use it in info columns

text: When I open mapbuilder and create a theme, I want to create ifo columns alias, like block name, block code, but the fields which datatypes are nvarchar2 does not apears in the list to add it to the list, any help? many thanks,


title: Is there a simple way to clip multiple layers and orthophotos?

text: (*disclaimer, I don't know any Python/scripting) Ok, so I have clipped several layers in the past using Arcmap 10.+, but I am having a hard time wrapping my head around how hard it can be at times. I used to own a photography/design business, so I'm used to manipulating crops via Photoshop. I think because of this, I have issues with simple clipping when composing a map. For this current clip/crop, I would like to freehand a polygon around several connected streams and "crop" out everything around that polygon. Within the polygon, there are several layers having many polylines and points, and a Ortho_tile_tiff catalog layer aerial view. I believe it's this Ortho_tile geodatabase that's throwing me here. If someone could break it down for a semi_green GISer I would be much appreciative. Thanks!


title: D3 adding topojson - no data visible

text: Yesterday I asked what would be the best method of creating an overlay in a web map that would allow for some interactivity. Create interactive map with no server Thanks to a variety of people who recommended I go down the D3 and Topojson route. The example that looks like the perfect solution is http://bl.ocks.org/mbertrand/5218300 So I created a topojson from an ESRI Shapefile using the following command               topojson --id-property full_name -o southeast.json southeast.shp      I have since tried every possible permutation and example to try and get the topojson loaded and I have hit a brick wall. I just do not understand what I am doing wrong. I have used the following examples http://www.schneidy.com/Tutorials/MapsTutorial.html http://stackoverflow.com/questions/14492284/center-a-map-in-d3-given-a- geojson-object Conversion from Shapefile to TopoJSON format results in incorrect map when rendered using d3js http://milkator.wordpress.com/2013/02/25/making-a-map-of-germany-with- topojson/ and Mike's original example http://bost.ocks.org/mike/map/ Maybe it is because i am not using EPSG:4326, but this example says that you do not have to use a projection http://bl.ocks.org/mbostock/5557726 The data can be found here - https://www.dropbox.com/s/dfxxi4xwrsdfa8t/southeast.json And this is my code               <script src="http://d3js.org/d3.v3.min.js"></script>     <script src="http://d3js.org/topojson.v1.min.js"></script>     <script>          var width = 500,         height = 500;          //null projection as my ESPG is not supported     var path = d3.geo.path()         .projection(null);          var svg = d3.select("body").append("svg")         .attr("width", width)         .attr("height", height);          svg.append("rect").attr('width', width).attr('height', height)         .style('stroke', 'black').style('fill', 'none');          //d3.json("southeast.json", function(error, south) {       //console.log(south);     //});          d3.json("southeast.json", function(error, uk) {       svg.append("path")           .datum(topojson.feature(uk, uk.objects.southeast))           .attr("d", d3.geo.path().projection(null));     });      I first tested with sending the data to the console and that worked fine and told me that the object is "southeast". If I try this I get an n is not defined error in Firebug. I then tried               // d3.json("southeast.json", function(error, south) {     //   svg.append("path")     //       .datum(topojson.feature(south, south.objects.southeast))     //       .attr("d", d3.geo.path().projection(null))     //       .style("fill", "red")     //       .style("stroke-width", "1")     //       .style("stroke", "black");     // });      Which seems to work as there are no Firebug errors, but no data is visible Then I tried               // d3.json("southeast.json", function(error, uk) {     //   svg.selectAll(".southeast")     //       .data(topojson.feature(uk, uk.objects.southeast).features)     //     .enter().append("path")     //       .attr("d", path)     // });      And the same thing happens With no errors I am struggling to understand what I have done wrong If I could just make this example work in EPSG:27700 with my topojson it would be amazing http://bl.ocks.org/mbertrand/5218300 Thanks # EDIT One of the other answers was to use a simply getJSON function and add the data as geojson in OpenLayers. I have just tested it using the same topojson:               $.getJSON('southeast.json', function (data) {                 fusiegem = topojson.feature(data, data.objects.southeast);                 var geojson_format = new OpenLayers.Format.GeoJSON();                 var vector_layer = new OpenLayers.Layer.Vector();                 osMap.addLayer(vector_layer);                 vector_layer.addFeatures(geojson_format.read(fusiegem));             });      And that successfully loaded, so it seems there is nothing wrong with my topojson file, it is the D3 bit that does not seem to work correctly


title: Update/refresh edited PostGIS table in QGIS

text: I imported a PostGIS table (with geometry) from my PostgreSQL database as a layer in QGIS. Later, I edited the table using PostgreSQL `UPDATE` commands. Is there a way to quickly update (like a **refresh ↻** command?) the layer in QGIS with the new attributes or geometry without having to remove the layer and reconnect with the database?


title: Merge 30000 polygon shapefiiles in a logical manner in ArcGIS model Builder

text: I have 30000 shapefiles (a1.shp to a30000.shp). They are all polygons. Each polygons has one and only one polygon feature in it. I need to Merge the polygons into 1000 shapefiles. My output would be 1000 polygons (each with 30 sequential polygons in it) So output b1.shp = Merged polygons a1.shp to a30.shp So output b2.shp = Merged polygons a31.shp to a60.shp So output b3.shp = Merged polygons a61.shp to a90.shp and so on up to... b1000.shp = Merged polygons a29970.shp to a30000.shp Looking for an ArcGIS model builder or GDAL solution.


title: Software to detect intersections (within a tolerance) between 3d polylines

text: Is there any way to detect intersections (within a tolerance) for a 3d set of polylines? ArcGIS Topology Check seems to project everything into 2d (x-y plane) before checking, so shows a lot of false positives. Autocad Map 3d Drawing Cleanup suffers from the same problem.