title: How to set Eclipse and PyDev for QGIS when ArcGIS is installed?

text: I have got a lot problems with the configuration of Eclipse with PyDev in order to develop QGIS plugins. The situation is very weird and frustrating. I can import qgis module using Pyscripter without any problems but in the case of Eclipse I have already spend a few hours to find out solution. I check a lot of tutorials from the Internet but I have not get answer to my question yet. This is not a duplicate question because I have checked all solutions form other questions. In my case I have got ArcGIS 10.2 and QGIS 2.2 installed on my machine. I use the script below to run PyScripter               @echo off     call "%~dp0\o4w_env.bat"     call "%OSGEO4W_ROOT%"\apps\grass\grass-6.4.3\etc\env.bat     @echo off     path %OSGEO4W_ROOT%\apps\qgis\bin;%OSGEO4W_ROOT%\apps\grass\grass-6.4.3\lib;%PATH%     set QGIS_PREFIX_PATH=%OSGEO4W_ROOT:\=/%/apps/qgis     set GDAL_FILENAME_IS_UTF8=YES     rem Set VSI cache to be used as buffer, see #6448     set VSI_CACHE=TRUE     set VSI_CACHE_SIZE=1000000     set PYTHONPATH=C:\OSGeo4W\apps\qgis\python     Set PATH=C:\OSGeo4W\apps\qgis\bin;%PATH%     set QGISHOME=C:\OSGeo4W\apps\qgis\     Start c:\PROGRA~2\PyScripter\pyscripter.exe --python27 --pythondllpath=C:\OSGeo4W\bin      In that case I can import qgis or qgis.core without any problems. When I check paths using :               import sys     for item in sys.path:         print item      I get this result:               C:\OSGeo4W\apps\qgis\python     C:\OSGeo4W\bin\python27.zip     C:\OSGeo4W\apps\Python27\DLLs     C:\OSGeo4W\apps\Python27\lib     C:\OSGeo4W\apps\Python27\lib\plat-win     C:\OSGeo4W\apps\Python27\lib\lib-tk     C:\OSGeo4W\apps\Python27     C:\OSGeo4W\apps\Python27\lib\site-packages     C:\OSGeo4W\apps\Python27\lib\site-packages\PIL     C:\OSGeo4W\apps\Python27\lib\site-packages\win32     C:\OSGeo4W\apps\Python27\lib\site-packages\win32\lib     C:\OSGeo4W\apps\Python27\lib\site-packages\Pythonwin     C:\OSGeo4W\apps\Python27\lib\site-packages\Shapely-1.2.18-py2.7-win32.egg     C:\OSGeo4W\apps\Python27\lib\site-packages\wx-2.8-msw-unicode      As the results I get ready working with QGIS Pyscripter but using this software I can not debug my plugin. In the case of Eclipse I use this bat file               @echo off     SET OSGEO4W_ROOT=c:\OSGeo4W     call "%OSGEO4W_ROOT%"\bin\o4w_env.bat     call "%OSGEO4W_ROOT%"\apps\grass\grass-6.4.3\etc\env.bat     @echo off     SET GDAL_DRIVER_PATH=%OSGEO4W_ROOT%\bin\gdalplugins     path %PATH%;%OSGEO4W_ROOT%\apps\qgis\bin;%OSGEO4W_ROOT%\apps\grass\grass-6.4.3\lib     set PYTHONPATH=%PYTHONPATH%;%OSGEO4W_ROOT%\apps\qgis\python;     set PYTHONPATH=%PYTHONPATH%;%OSGEO4W_ROOT%\apps\Python27\Lib\site-packages     set QGIS_PREFIX_PATH=%OSGEO4W_ROOT%\apps\qgis     "C:\Progra~2\eclipse\eclipse.exe"      When I use piece of code mentioned above to check paths I get:               C:\Python27\ArcGIS10.2\lib\site-packages\setuptools-2.2-py2.7.egg     C:\Python27\ArcGIS10.2\lib\site-packages\sphinx-1.2.1-py2.7.egg     C:\Python27\ArcGIS10.2\lib\site-packages\jinja2-2.7.2-py2.7.egg     C:\Python27\ArcGIS10.2\lib\site-packages\docutils-0.11-py2.7.egg     C:\Python27\ArcGIS10.2\lib\site-packages\pygments-1.6-py2.7.egg     C:\Python27\ArcGIS10.2\lib\site-packages\markupsafe-0.18-py2.7.egg     c:\OSGeo4W\apps\qgis\python     c:\OSGeo4W\apps\Python27\Lib\site-packages     c:\OSGeo4W\apps\Python27\DLLs     c:\OSGeo4W\apps\Python27\lib     c:\OSGeo4W\apps\Python27\lib\lib-tk     C:\OSGeo4W\bin     c:\OSGeo4W\apps\Python27     c:\OSGeo4W\apps\Python27\lib\site-packages\PIL     c:\OSGeo4W\apps\Python27\lib\site-packages\win32     c:\OSGeo4W\apps\Python27\lib\site-packages\win32\lib     c:\OSGeo4W\apps\Python27\lib\site-packages\Pythonwin     c:\OSGeo4W\apps\Python27\lib\site-packages\Shapely-1.2.18-py2.7-win32.egg     c:\OSGeo4W\apps\Python27\lib\site-packages\wx-2.8-msw-unicode     C:\OSGeo4W\bin\python27.zip     C:\Python27\ArcGIS10.2\Lib     C:\Python27\ArcGIS10.2\DLLs     C:\Python27\ArcGIS10.2\Lib\lib-tk     C:\Python27\ArcGIS10.2     C:\Python27\ArcGIS10.2\lib\site-packages     C:\Program Files (x86)\ArcGIS\Desktop10.2\bin     C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy     C:\Program Files (x86)\ArcGIS\Desktop10.2\ArcToolbox\Scripts      Of course I set Forced Builtins for qgis and PyQt4 but still It does not work. When I try to run import qgis I get this error and I get error               from qgis.core import QgsFeature, QgsGeometry     ImportError: DLL load failed The specified module could not be found.      Any idea how to solve this problem?


title: KML Balloon Not Showing Up In <Folder>

text: I've created a perl script to create multiple line files from points. It iterates through, creates everything as expected. When you hover over a line, I get a balloon with all of the attribute information that I'm pushing. Only issue, I don't want 300 separate kml files. So, I modified the script to write to a kml "folder". All of the lines display as expected. Looks just like the individual files, but I can open the one file and see everything. EXCEPT, now the attribute information does not show up when I hover over a line/click on a line/even when I go to properties. It's the same exact setup that is in the individual file...Is it impossible to show several balloons in one file, even if they are separate "documents"? Thanks for your help! ![Balloon Tag Snip](http://i.stack.imgur.com/bj5xQ.jpg)


title: How to manage the osmosis command line?

text: This refer to this post : Read Here I'm in trouble when using osmosis command line to export data from planet- latest.osm.bz2 directly to osm2pgsql using some criteria to select what to transfer. I think is more easy to export planet-latest.osm.bz2 to an XML with osmosis and then use osm2pgsql to import to Postgre, but I can't go further. So I need help in osmosis parameters to show how can I to export data using filters like               --place='city' --admin_level='2'      and import to Postgre using osm2pgsql. Sorry about this very similar post, but the user was helping me was gived up.


title: Calculating line direction in MapInfo

text: I have a large file of a road network. I need to calculate the value of direction (in degrees or radians) of the line segments, and save this value to a new field. I need to use MapInfo. I have found this excellent example using qgis, but need to undertake the work using MapInfo. Does anyone have a solution. I tried the basic idea of using the function listed in the solution and copy and pasting that in, but it does not work. How to add Direction and Distance to attribute table? Hope this all makes sense. To conclude, I need to calculate direction for each line segment and add them to a new field in MapInfo.


title: How can I calculate Moran's I in 3D?

text: This may be a question to post on a GIS forum, but I was wondering if among the spatial data experts here, someone could help me calculate the extension of Moran's I to 3D? Our data is about 3D habitat complexity, specificlaly, the 3D complexity of the habitat provided by algae to other marine organisms. The idea is to estimate Moran's I using 3D data we have in which each voxel (3d pixel) is either water or algae. So, the neighbours for each point are distributed in 3d space. Norbert Marwan used this to calculate bone complexity (medical journal) but there is something about how poitns are accounted for in the index calculation avoid repetition of neighbour pairs that I don't understand Any help would be most appreciated. All the best, Ana


title: Openlayers GetFeatureInfo not working

text: I installed apache2 and GeoServer, then set up the proxy.cgi. I imported shp based multipolygon file into geoserver. I can access via geoserver openlayers viewer but when I want to do it same in a html file I cant do it, because the getfeatureinfo dont want to work properly. I added to the file:               OpenLayers.ProxyHost = "proxy.cgi?url=";               map.events.register('click', map, function (e) {                         document.getElementById('nodelist').innerHTML = "Loading... please wait...";                         var params = {                             REQUEST: "GetFeatureInfo",                             EXCEPTIONS: "application/vnd.ogc.se_xml",                             BBOX: map.getExtent().toBBOX(),                             SERVICE: "WMS",                             INFO_FORMAT: 'text/html',                             QUERY_LAYERS: map.layers[0].params.LAYERS,                             FEATURE_COUNT: 50,                             Layers: 'csongrad_demo:csongmat',                             WIDTH: map.size.w,                             HEIGHT: map.size.h,                             format: format,                             styles: map.layers[0].params.STYLES,                             srs: map.layers[0].params.SRS};                                   OpenLayers.loadURL("http://xxx.xxx.xxx.xxx:80/geoserver/csongrad_demo/wms", params, this, setHTML, setHTML);                         OpenLayers.Event.stop(e);                     });       after I tried this:               info = new OpenLayers.Control.WMSGetFeatureInfo({                         url: 'http://xx.xx.xx.xx:80/geoserver/wms',                          layerUrls: 'http://xx.xx.xx.xx:80/geoserver/csongrad_demo/wms',                         title: 'Identify features by clicking',                         queryVisible: true,                         eventListeners: {                             getfeatureinfo: function(event) {                                 map.addPopup(new OpenLayers.Popup.FramedCloud(                                     "chicken",                                      map.getLonLatFromPixel(event.xy),                                     null,                                     event.text,                                     null,                                     true                                 ));                             }                         }                     });                     info.activate();                     map.addControl(info);      none of them wanted to work, I always get this error "Loading... please wait..." please help Thanks


title: QGIS 2.2 saving layer as CSV with GEOMETRY=AS_WKT not working

text: Just upgraded to QGIS 2.2 and it looks like saving a layer as CSV doesn't support GEOMETRY=AS_WKT anymore. The new Layer Options portion of the save dialog is cool but it looks like GEOMETRY and LINEFORMAT for CSV have some swapped options (CRLF and AS_WKT are in opposite option lists). I tried custom layer options:               GEOMETRY=AS_WKT     SEPARATOR=SEMICOLON     LINEFORMAT=CRLF      Those worked under QGIS 2.0 but don't appear to work with QGIS 2.2.


title: data defined labeling with postgis layer

text: How to use the expression based (edit: data defined) labeling with postgis? I'm wondering, whats the best way to use expression based (edit: data defined) labels with postgis (I'm new to postgis). Just by adding new columns for x and y and so on, or by using a relational table with that columns? Or is there a other way? Thanks! EDIT: The image below shows the common way to use data defined labels for shape files: Add a column for x and y position for the label (and/or color, size, etc.) to the attribute table and set the correct field in the layer preferences. Now I would like to know if its possible to store the x and y (etc.) columns into a second relational table (with postgis). Because I would like to use the Data-Defined-Function without modifying the attribute table itself. ![enter image description here](http://i.stack.imgur.com/vwnyT.png)


title: How do I stop Kernel Density from cutting off symbology?

text: I need help configuring Kernel Density so that is does not cut off symbology on hot spots. (See pic below) I'm assuming it has something to do with the extent but I could not get it to display the way I want it. Any help would be appreciated. ![enter image description here](http://i.stack.imgur.com/RMmgk.jpg) ![enter image description here](http://i.stack.imgur.com/C8Pq8.jpg) ![enter image description here](http://i.stack.imgur.com/42MDK.jpg)


title: Is there a bug when creating a relationship class with Python using a GlobalID?

text: I'm trying to create a relationship class with Python using a GlobalID. When I use Python it gives me the error: > Failed to execute. Parameters are not valid. ERROR 000800: The value is not > a member of GLOBALID. Failed to execute (CreateRelationshipClass). When I try to use the Create Relationship Class tool through the Toolbox, I get the same error. If I select GlobalID as the primary key, it will ONLY let me select a GlobalID as the foreign key. It will not let me select a GUID of my choice. But when I'm in ArcCatalog and right click the Dataset and go to `New > Relationship Class`... It works! What's wrong here? Thanks.


title: Where can I find GDAL bindings for Python 3.3?

text: I have been using GDAL 1.9.x on Python 2.7 and took the bindings from here. Is there a place to get GDAL bindings for Python 3.3 or is GDAL not officially supported in Python 3.3 yet?


title: Calculate maximum and minimum elevation for a polygon layer from a DEM

text: I have a polygon layer and would like to know the min and max elevation within each polygon. My elevation layer is a digital elevation model. I'm not quite sure how to extract that information. I am using arcMap 10 Any help would be apreciated. thanks!


title: How to determine aerosol model value for i.atcorr in GRASS?

text: The `i.atcorr` module in GRASS requires aerosol model data to perform atmospheric correction using the 6S algorithm. It seems like I need **a value from one of these parameters** to continue:   * Estimate of the **meteorological parameter visibility** v [km], or   * The **aerosol optical depth (AOD) at 550nm** **Is there an easy (and free) way to get such data**? * ***EDIT (new findings):** It seems that a NASA website called Space-Based Measurements of Ozone and Air Quality in the Ultraviolet and Visible got the **free aerosol model data**. It is organized by date and by the instruments that collected the data. Basically:   * **OMI (Ozone Monitoring Instrument)** is for data dated July 2004 onward   * **TOMS Earth Probe** is for data dated 7/25/96 - 12/31/05   * **TOMS Nimbus-7** is for data dated 11/1/78 - 5/6/93 I'm processing a Landsat scene dated in year 2000-2001, so Toms Earth Probe data is the way to go. **But as stated by @Tobi in an answer below, TOMS Earth Probe data is actually an index, and unitless** (not really AOD at 550nm which I want in the 6S algorithm). But I've stumbled upon this called A Long-Term Record of Aerosol Optical Depth from TOMS Observations and Comparison to AERONET Measurements, and they _hinted_ how the data in TOMS observation can still be used. I need more time to learn how to do this. Welcome any further input, always..


title: How to remove/disable a toolbar via python?

text: How do I remove or disable a toolbar e.g., Raster, programmatically using Python?


title: What is the difference between EPSG:900913 and EPSG:3857?

text: I'm using the QGIS google layers plugin to digitize a parcel as a shapefile. Later I want to import it to PostGIS. I know that google maps use a special 'google mercator' projection. What is the correct code of this projection: 900913 or 3857? Or are these two codes describing the identical projection equations?


title: Postgis geometry type = Geometry

text: Simply I've got a number of PostGIS layers that have a geometry type of Geometry, as a posed to point, line, polygon.               wkb_geometry geometry(Geometry,27700)      QGIS 2.0+ doesn't like this very much and wants me to set the geometry type and in general makes my PostGIS connect in QGIS slow. How do I change the Geometry type to point for instance? Si


title: How to rectify polygon shapes in QGIS?

text: I have some polygon-shapes in qgis-1.7.0 and want to ensure that every polygon is a rectangle. See image below. ![enter image description here](http://i.stack.imgur.com/aEYy8.png) Is there any way I can do this in QGIS (or any other f/oss)? Or how can I improve my polygon-drawing to make sure for example that every building consists of rectangles (all angels = 90°)? Any hints appreciated! Thanks!


title: how to specify -t critical dowloading raster with raster2pgsql

text: I'm trying to download a raster into postgres/postgis. So far I wrote:               raster2pgsql -s 3006 -I -C -M -R filename.tif -F -t 463.313x463.313 schema.rastername | psql -h localhost -U postgres -d db      The raster is dowload that way successfully but my problem comes with the fact that the output raster file (once download read in postgres) is smaller than the original one. My raster has 2960rowsx2558columns, each pixel has 463.313x463.313 for dimension, and I have a scale factor of 1. I don't understand how to specify the dimensions. I have read these posts _"Difference in raster size using raster2pgsq"_ and _"Loading a raster into a PostGIS 2.0 database on Windows"_ but I did not get how the conversion is done.


title: I'd like to accelerate the refresh of a Open Layers vector layer throught a jqueryui slider in Firefox

text: I want to develop a kind of timemap combining the jquery-ui slider component and the OpenLayer library. I'm simply using the Vector layer of OpenLayers with Rules style.                map = new OpenLayers.Map('map', {maxResolution:'auto'});     var wms = new OpenLayers.Layer.WMS( "OpenLayers WMS",                                        "http://vmap0.tiles.osgeo.org/wms/vmap0",                                        {layers: 'basic'} );     map.addLayer(wms);     map.setCenter(new OpenLayers.LonLat(0, 0), 0);          // create 20 random features with a random type attribute. The     // type attribute is a value between 0 and 2.     var features = new Array(20);     for (var i=0; i<20; i++) {         features[i] = new OpenLayers.Feature.Vector(             new OpenLayers.Geometry.Point(Math.random()*360-180, Math.random()*180-90),             {type: parseInt(Math.random()*80)}         );     }            var style = new OpenLayers.Style(         {                                        graphicHeight: 20,             graphicYOffset: -19,             rotation: "${angle}"         },         {             rules: [ new OpenLayers.Rule({                          filter: new OpenLayers.Filter.Comparison({                             type: OpenLayers.Filter.Comparison.LESS_THAN,                             property: "value",                             value: 20                         }),                                                      symbolizer: {                             externalGraphic: "http://dev.openlayers.org/releases/OpenLayers-2.13.1/img/marker-blue.png"                         }                     }),                     new OpenLayers.Rule({                             filter: new OpenLayers.Filter.Comparison({                                 type: OpenLayers.Filter.Comparison.BETWEEN,                                 property: "value",                                 lowerBoundary: 20,                                 upperBoundary: 40                                  }),                                                          symbolizer: {                                 externalGraphic: "http://dev.openlayers.org/releases/OpenLayers-2.13.1/img/marker-green.png"                             }                     }),                     new OpenLayers.Rule({                             filter: new OpenLayers.Filter.Comparison({                                 type: OpenLayers.Filter.Comparison.GREATER_THAN,                                 property: "value",                                 value: 40                             }),                                                          symbolizer: {                                 externalGraphic: "http://dev.openlayers.org/releases/OpenLayers-2.13.1/img/marker-gold.png"                             }                    })             ]         }     );          layer = new OpenLayers.Layer.Vector('Points', {         styleMap: new OpenLayers.StyleMap(style)     });          layer.addFeatures(features);     map.addLayer(layer);      I change the displayed value though the slide event of the jquery-ui slider.                $( "#slider" ).slider({             value:0,             min: 0,             max: 100,             step: 1,             slide: function( event, ui ) {                 var end = false;                 //launched.push(ui.value);                 $( "#value" ).val( ui.value );                      function redraw(){                     //launched.pop();                     for(var i = 0 ; i < features.length; i++){                         features[i].attributes["value"] = Math.random()*60;                                  }                     map.getLayersByName("Points")[0].redraw();                      };                 redraw();                       }         });      It work very well with few elements. But if I set the number of displayed element to 200 and open my script with Firefox, the slider is very slow and do not follow the mouse cursor, although it work well with Chrome. You can find my source code on this jsfiddle : http://jsfiddle.net/jeje/5XPcc/6/ Anybody have an idea of what I did wrong ? Do you think there is a way to improve it ? Already asked in stackoverflow https://stackoverflow.com/questions/21626115/how-to-improve-openlayer-vector- layers-refresh-for-200-elements


title: Determining gaps in polylines and identifying continuing sections?

text: Is there a relatively quick way to identify where gaps < a certain size are present between polyline? E.g. would there be an easy way to join line likes the ones in the below image if the gap was <10m? ![enter image description here](http://i.stack.imgur.com/5T6l1.png) I have access to ArcMap, MapInfo and QGIS.


title: Nearest point on a line ( Spherical / Mercator projection )

text: I have a **line (Ax,Ay - Bx,By)** over a mercator projection (google maps) and **a random point (Cx,Cy)** nearest to that line, i would to know the closest point (transparent blue on the image) over that line to point (blue in the image) **EDIT: to clarify that this is in a Mercator projection (spherical projection)** ![enter image description here](http://i.stack.imgur.com/ZSpF0.jpg)


title: Georeferencing CAD Files and Creating Geodatabase

text: I'm hoping this is an easy question to answer. I have georefrenced a CAD data set to an orthophoto that was projected in NAD 83, and turned it into a GDB. Even though the GDB is not projected, is this still considered to be Georefrenced since its based off of the CAD files? When I bring in a file from the GDB I get a warning saying "The following data sources you added are missing spatial reference information. This data can be drawn in ArcMap, but cannot be projected" In order to put it into a GDB, I used the select by attributes to only take the certain features I was working with, and exported them into their own Polygon or Polyline file. I then created a GDB and imported the files into it.


title: Running a Focal Mean or Low Pass Filter or 3 by 3 Mean Kernel in GDAL \ gdal_calc

text: Hi I am looking to run the equivalent of Focal Statistics in gdal_calc. Is it possible? It is a simple focal I require. I need to run a 3 by 3 kernel over the data (mean) and if possible only on NoData (but if on all that is okay as well). I can do this in other software but specifically asking if a GDAL_calc solution exists. This is the exact thing I am trying to do on a singular raster Con(IsNull("C:\temp\SST_Final_C_NoData.tiff"), FocalStatistics("C:\temp\SST_Final_C_NoData.tiff", NbrRectangle(5,5, "CELL"), "MEAN"), "C:\temp\SST_Final_C_NoData.tiff") GDAL 1.11.0 and Python 3.4.1


title: Latitude & Longitude values from Simple Voyage Data Recorder [SVDR]

text: I am receiving values in the form of `0113.7276,N,10339.0121,E` from a Ship Data Recorder. I am not a GIS/ GPS expert and am unable to decipher what do these values mean, because they are out of range for standard latitude and longitude values. What could these values mean?


title: How to get output value from Feature Compare (schema) and append to existing csv?

text: I am working on a data integration project and am trying to single out feature classes which have a schema that varies from two base feature classes (one for point, one for polyline features). I've played around with the Feature Compare tool to compare schemas, but find it tedious to use this tool, even in batch mode, because it outputs a single text file with a few values for _each comparison performed_. I started writing a script in Python to loop through a list of features in a list and compare their schema based on geometry type. However, the problem of either integrating the output into a single text file persists. Is there a way to obtain this value and append it to an existing text file? (Messy) code so far as follows:               # Import arcpy module     import arcpy          # Import environments; set workspace     from arcpy import env          # Set workspace     env.workspace = "Data_Comparison.gdb"          # List features to be compared     fclist = arcpy.ListFeatureClasses()     print fclist          # Loop through list to separate point and line feature classes     for fc in fclist:         desc = arcpy.Describe(fc)         name = desc.file         geo = desc.shapeType              print ("Feature name: " + name)         print ("Feature geometry: " + geo)              if geo == "Point":             arcpy.FeatureCompare_management("swr_net_nodal", fc, "ENTID", "SCHEMA_ONLY", "", "", "", "", "", "", "CONTINUE_COMPARE", "fc + '_compare'")             output = arcpy.GetMessages()             print output


title: occasionally empty search cursor looping through a raster list with arcpy

text: I'm looping through a raster list, to export cell counts (from "Count" field in the table) to a text file, I'm only interested in cases where the field "Value" = 1, so I have included this in my search cursor (inR is from the raster list):               Rows = arcpy.SearchCursor(inR, '"Value" = 1', "", "", "")      However, some rasters in the list do not have value = 1 in the table, in this case, I want to return 0 for this count. I've tried to include this in the following loop through the search cursor               for row in Rows:         countR = row.getValue('Count')         if "countR" not in locals():             countR = 0          #append count to list to export later     value_list.append(countR)     del(countR)      When the search cursor is "empty", the get.Value line does not return any value and the object is not created, however, the following line does not work, the idea is that if countR does not exist, then it will be given the value of 0 for this raster in the list (inR). Edit: this is where I went wrong, the above is not correct, when the search cursor is empty (like an empty list) , there is no loop and none of the statements within that loop take place. Any help would be appreciated, or alternatives. Thanks


title: How to publish GeoTIFF using a WCS service with QGIS Server

text: I need to publish some GeoTIFFs using QGIS Server (2.2.0-1+precise1). I already have installed and setup QGIS Server to serve those geotiffs via WMS. But I need to serve rather the raster data to be able to use it in raster calculations. I don't know how to set up or call the service from an URL. I tried with the parameter SERVICE=WCS instead of SERVICE=WMS, to no avail. If QGIS server can't handle WCS yet, what other open source alternatives are there? Thank you for your time and help.


title: How to determine if an object is a circle

text: I have a polygon shapefile that consists of many different shaped features. I would like to identify only those shapes that are circles. Is there a way to do this? Is it a calculation I could do in the field calculator? Thanks... ![enter image description here](http://i.stack.imgur.com/xiuWr.jpg)


title: Make sidebar div active on Marker click (and vice versa)?

text: What's the best way to change an element to active when you click on a marker. For example, when you click on a marker on a Leaflet map, how to make it so the sidebar element that's related to it will be active? I can write the sudo code but I'm missing pieces.               markers on click     $(#sidebar div)        make (related) li element active      edit: Basically, I'm trying to mimic what a lot of map directory sites like Yelp or Foursquare do. You click on the sidebar list item and the marker becomes active, and vice versa. This makes it so a user and can read info about the marker while seeing where the place is. Mapbox had a cool (now deprecated) Foursquare example. I went through the code for that but barely understood it. That said, the example gets the point across. ![Foursquare Mapbox example](http://i.stack.imgur.com/sT9tz.png) When you mouseover a listitem in yelp, the marker goes black. ![Yelp example](http://i.stack.imgur.com/VY13w.png)


title: Is "puntal" an accepted and defined geometric term?

text: There are definitions for the geometric qualities of polygons and lines: >   * **polygonal** : Having many angles; hence characteristic of a polygon. > Ref: http://en.wiktionary.org/wiki/polygonal >   * **lineal** : of, relating to, or consisting of lines; linear. > Ref: http://oxforddictionaries.com/us/definition/american_english/lineal However, I have been unable to find a scholarly reference for " **puntal** ", that is, having the geometric qualities of a point. The only solid references to "puntal" (in English) I have found are in the discussions and source code of projects like GEOS and JTS Topology Suite, where there might be classes or interfaces corresponding to **Polygonal** , **Lineal** , and **Puntal**. Is there a canonical definition of the geometric sense of this term?


title: Label polygons in MapInfo like QGIS?

text: I was wondering if there is a way to angle Polygon labels in MapInfo like it's done in QGIS (see image link below). ![enter image description here](http://i.stack.imgur.com/UA368.png) Thanks, Alan


title: PostGIS is rejecting an srid code for my projection. I've found a nearly identical projection w/ a legit srid. Will it work?

text: My projection is this: `NAD 1983 StatePlane North Carolina FIPS 3200 Feet`, which has a proj4 string that looks like this:               +proj=lcc +lat_1=34.33333333333334 +lat_2=36.16666666666666 +lat_0=33.75 +lon_0=-79 +x_0=609601.2199999999 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs      Its SRID is 102719. I've tried to create a db column in PostGIS with an srid option for this projection, but its rejecting the srid, calling it invalid. I've found a similar projection with a valid srid code, `NAD83 / North Carolina (ftUS)` (srid=2264), which has a proj4 string nearly identical to the above:               +proj=lcc +lat_1=36.16666666666666 +lat_2=34.33333333333334 +lat_0=33.75 +lon_0=-79 +x_0=609601.2192024384 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs       As you can see, the only difference between these two is that the `lat_1` and `lat_2` degrees are swapped. They even use the same three degrees, only the first two are differently ordered. Can I use this second projection? UPDATE: The answer seems to be leaning towards 'yes', I just ran projection- to-projection transformations from a single projection to the NC projection using the "legitimate" and "illegitimate" proj4's and got similar results down to the 7th decimal place


title: Google route like functionality in Arcmap

text: Drawing routes in Google Earth or Google Maps looks very cool. As point move along the route with specific speed. I am looking same functionality in Arcmap. Can anyone figure out how?


title: Split data for OpenLayers

text: (I gave up. I already spend 10 hours reading internet and I am not any closer to get the job done.) If I start at the beginning:   * I have country map divided into 12 regions.   * I have country map divided into 56 areas.   * I have country map divided into 211 townships.   * I have country map divided into 2700 local areas.   * I have buildings map approx. 500 thousands records. ![enter image description here](http://i.stack.imgur.com/FDjnX.png) I want to build interactive map, but on old computers OpenLayers(geojson) with 211 records is kind of slow. My idea was to show user only regions, then a user select a region and show areas for selected region, user select an area, and then I show user only townships for selected area and so on. I have no problem manually connect townships to areas and areas to regions, just need to split them. My first idea was to split polygons in layer (Areas) into separate files with QGIS and load correct file in OpenLayers. But then when searching the web how to auto split polygons in layer and save them in separate files I asked myself should I use Geoserver or something similar. I hope I described my problem good enough to get some advice, should I use GeoServer or use separate files and how can I save polygons in layer to separate files? How should I approach the problem? Thanks.


title: Weighted Overlay with "nodata" cells

text: I´m currently conducting a suitability analysis with the weighted overlay function. One relevant criteria is the proximity to highways. If the wind turbine is located up to 300 m away from the highway, the cell gets 10 points (out of 10). If it is farther away,the cell shouldn´t get a score, because the criteria is not relevant. So I ran the "euclidean distance" function with a maximum distance of 300 m. All cells farther away than 300 m are declared as "nodata". Then I ran the "weighted overlay" function with 3 more criteria (all have the same weights). The problem is that the study area is bigger than the 300 m buffer around the highways, but after I ran the "weighted overlay" function, only the area within the 300 m buffer was weighted. The rest of the study area was declared as "nodata". The goal is that, when conducting the weighted overlay function, a cell farther away than 300m is not influenced by the criteria "proximity to highway". example: There are 4 criteria with equal importance. C4 is the "proximity to highways" criteria. A cell within the 300 m radius should be calculated like this: c1: 8 points (25%), c2: 6 points (25%), c3 4 points (25%) and C4: 10 points (25%) --> weight: (8+6+4+10)/4 = 7 points A cell outside of the 300m radius should be calculated like this: c1: 8 points (25%), c2: 6 points (25%), c3 4 points (25%) and C4: no input (25%) --> weight: (8+6+4)/3 = 6 points Is this possible? I don´t want to give the "nodata" cells a scale value. Thanks for you help!


title: How do I use variables from one Python script in other Python scripts?

text: I want to set up a single python file containing variables that have the locations of all my data sources. This would then be used by all my other scripts, then as a data source changes, I only have to edit the one file. My Data_sources.py would look something like this:               BC_BEC = "Database Connections\\BC.sde\\FOREST_BC.BC_Eco\\BC.BC_BEC"     BC_TFL = "Database Connections\\Forest_BC.sde\\BC.BC_Admin\\BC.TFL_Boundary"      and so on... How do I get all my other scripts to load or access these variables in the Data_sources.py? Thanks,


title: Can I make tiles with all countries but one transparent?

text: I would like to have one country in a given style and all other countries in another grayed out one. But styles should be changeable in the production webmap, so I need them to be independent of each other. My idea for this was to render the whole world in the background style, and the desired country as an overlay. By overlay, I mean render all areas outside of this country transparent. How can I do this in TileMill?


title: How can I get the field index of the new layer in my script using Python?

text: I am trying to use `layer.changeAttributeValue(fid, fieldIndex, value)` to add values to a field in QGIS Basically, I'm creating a plugin which will, if a checkbox is ticked, create a new field, then iterate through the dataset, filling that field with the x (easting) coordinate of each point (I also have another checkbox for the same for the y coordinate). As such, I cannot set the field index beforehand because I don't know which will be ticked. So - how can I get the field index of the new layer in my script using Python? Also (sorry, I'm new to this), how can I get the `fid` of the current feature as I iterate through the dataset?


title: customize handlerOptions

text: when i customize handlerOptions of drawfeature, savestrategy not work. Thats my code :               var handoptions = {             handlerOptions: {                 freehand: true,                 style: {                     fillColor: "yellow",                     fillOpacity: 0.4,                     strokeColor: "green",                     strokeOpacity: 1,                     strokeWidth: 1,                     pointRadius: 6                 }             }         };              draw = new OpenLayers.Control.DrawFeature(            wfs_layer, OpenLayers.Handler.Point,handoptions,            {                //  title: "Draw Feature",                //  displayClass: "olControlDrawFeaturePolygon",                multi: true,            }            );      i dont know why ?


title: How to migrate a specific workspace, store and layer in Geoserver

text: Geoserver is running on a development and production environment. I would like to migrate just one workspace, store, layer, style, etc from development to production. How can I do this? I know I can copy the entire data-dir to production and restart geoserver; but than all workspaces, stores, etc are migrated. The point is I want just one workspace, store, etc. In the development environment I have several test workspace, store, and layers, which I do not want on production. In the data-dir I see a directory gwc-layers, with a lot of xml-files. But how do I know which xml is which layer?


title: How to add a new column in attribute table with attribute with indexes of another max-min ordered attribute field?

text: I've an attribute column named e.g. `VALUE` with different values I would like to make a column named e.g. `ORDER`such that the greatest value has the attribute `1` the second value `2` etc. etc. so order express the order from the biggest to the smallest value like this:               VALUE   ORDER     12000     1     8000      2     600       3     500       4     200       5      So far I've tried with QGIS 1.8 to order the `VALUE` column (that please note is not in the right order by default) from `max`to `min`simply clicking over the attribute name and than use `$rownum`in field calculator but this seems use the default rows numeration so I don't have the right order. Please edit my question if you think you can explain it better!


title: Customize search feature

text: I'm trying to create a map with CartoDB to show the spanish population between...at a municipality level. I can do everything perfectly, but I have just one problem: I would like to know if is there any possibility to customize the "search" feature, connecting it to the geolocalized names that I have in my table. In the way that if I put "Guadalajara" -that I have in my table- it goes directly to Guadalajara, Spain (the one in my table) and not to Guadalajara, México. And, in the same way, If I put "Madrid" right now it goes to the city, just a bit more north than where my point is. Hope I've been clear. Maybe my solution is to add a filter feature to a map, but I still don't know how to do it. Thanks!


title: How to export a shapefile from ArcGIS 10 into Autodesk Ecotect?

text: I am trying to export a shapefile from ArcGIS 10 into Autodesk Ecotect Analysis 2011 for solar analysis. I have a shapefile with polygons, and the heights of buildings is stored in the attribute table. I tried saving it as a KML file but Ecotect does not import KML files. Instead, Ecotect would accept any DXF or 3D files (such as 3D Metafile and 3D Studio). Is there a way of saving 3D files from GIS other than KML? I tried importing KML into Sketchup and Autocad without success. Any other idea for exporting shapefiles into a software that can extrude them and save as 3D would help too! Cheers!


title: Global solar radiation data

text: I need a solar radiation data in central Asia. Is there any global dataset of solar radiation. If someone know about it, please send to me the link. thanks a lot


title: ArcGIS Find task performance over SDE

text: We are having a problem in our Search feature for our app. We are using arcgis published REST services and SDE database. When we try to do our search manually using the REST url, it takes roughly 9 minutes for it to finish. Although it is searching an LSD layer, i don't think it should take that much time. My question is: Is there anything we can do to improve the performance of the find task? We already did the compress, analyze and rebuild indexes in our SDE database. Thanks.


title: Code works on python console in QGIS but same doesn't work in stand alone script

text: I'm trying map rendering with stand alone script, so I have small script which works fine in python console provided QGIS, however when I run same code in stand-alone script it executes without error, but output is not produced. Here is my standalone script:               from qgis.core import *     from qgis.gui import *     from qgis.utils import *     from PyQt4.QtCore import *     from PyQt4.QtGui import *     from PyQt4.QtXml import *          QgsApplication.setPrefixPathddd("C:\Program Files (x86)\QGIS Valmiera", True)          QgsApplication.initQgis()          mapinstance = QgsMapLayerRegistry.instance()     mapinstance.removeAllMapLayers()          layer = QgsVectorLayer("F:\scripts\map\wb.shp", "test", "ogr")          if layer.isValid():         print "shplayer succesfully loaded - adding to mapinstance"         mapinstance.addMapLayer(layer)     else:         print "shplayer failed to load!"          print mapinstance.mapLayers()          # create image     img = QImage(QSize(800,600), QImage.Format_ARGB32_Premultiplied)     # set image background color     color = QColor(255,255,255)     img.fill(color.rgb())     # create painter     p = QPainter()     p.begin(img)     p.setRenderHint(QPainter.Antialiasing)     render = QgsMapRenderer()     # set layer set     lst = [ layer.id() ] # add ID of every layer     render.setLayerSet(lst)     # set extent     rect = QgsRectangle(render.fullExtent())     #rect.scale(1.39800703)     rect.scale(1.1)     render.setExtent(rect)     # set output size     render.setOutputSize(img.size(), img.logicalDpiX())     # do the rendering     render.render(p)     p.end()          # save image     img.save('F:\scripts\out.png',"png")              QgsApplication.exitQgis()


title: Are there PhD programs that specifically focus on GIS, and if so, what is the typical outcome?

text: I know of the University of Pittsburgh where you can get a PhD in information science that focuses on Geoinformatics, which seems kind of cool. I've been looking around but it seems that most PhD programs in Geography focus on either human or physical and don't really get into studying methods of GIS, which is what interests me. I'm sure there are some that I'm missing though.


title: File representing extents of coordinates systems and their wkid

text: I'm looking for a file that represents the major coordinate systems throughout the world, and the extent that they cover. Ideally they'll also have the wkid as an attribute. Something like the shapefile found here of UTM zones, but with wkid's http://earth-info.nga.mil/GandG/coordsys/grids/universal_grid_system.html Cheers!


title: ArcObjects Workspaces with spaces

text: I am programming with ArcObjects and trying to output a raster to a folder path that contains spaces. Before creating the raster I am trying to validate that the destination path is valid. The following call to IsWorkspace() returns false when the path has spaces. Is this by design or a problem? How should one test that a folder is a valid file-based raster workspace?               Dim aType As Type = Type.GetTypeFromProgID("esriDataSourcesRaster.RasterWorkspaceFactory")     Dim obj As System.Object = Activator.CreateInstance(aType)     Dim pWSF As ESRI.ArcGIS.Geodatabase.IWorkspaceFactory = obj     return pWSF.IsWorkspace("C:\MyGISData\My Output Folder")


title: Hiding Map Elements Outside of Neatline in ArcMap

text: I have a map layout in ArcMap 10.1 that incorporates a neatline that is offset from the edges of the layout by about an inch on the inside. Certain elements of my map extend beyond the neatline and layout view. Does anyone know of a way to clip or hide the elements that extend between the neatline and edge of the layout when exporting the map? Thanks!


title: Help me to understand how i.fusion.brovey works in Grass 6.4.3

text: I'm using Grass 6.4.3 in Linuxmint 16. I don't have a very big problem but I'm trying to understand what is Grass doing, if I'm applying the process correctly and if there is something wrong with my image files. I'm trying to pansharpen a SPOT image. I have two files, one the panchromatic band and the other is the multiband file. I import them to grass as, for example, 'spotimage.pan' (panchromatic) and the other file is automatic named by grass as 'spotimage.1', 'spotimage.2', 'spotimage.3' and 'spotimage.4'. I think that in the original file, the bands order is inverted, I mean that the first band is NIR, the second one is Red and the third one is Green. This is because original image composite was done in order to obtain a false color image when is displayed as R,G,B = B1, B2, B3 So, if I do               d.rgb -o red=spotimage.1 green=spotimage.2 blue=spotimage.3      I obtain a false color view of the image. Then, to pansharpen the image I use               i.fusion.brovey -s ms1=spotimage.3 ms2=spotimage.2 ms3=spotimage.1 pan=spotimage.pan outputprefix=spotimage.pansh      I inverted the order, where I should use the SPOT Band 1, I'm using the band 3 of my image and where I should use the SPOT Band 3, I'm using the band 2 of my image. This because, in the dialog box of `i.fusion.brovey`, where it says _"Name of input raster map (green: tm2 | qbird_green | spot1)"_ I suppose I have to use the name of the SPOT Band 1 and, where it says _Name of input raster map (MIR; tm5 | qbird_red | spot3)"_ I suppose I have to use the name of the SPOT Band 3. Is it correct? The output are three raster files: _pansh.red, pansh.green_ and _pansh.blue_. But I cannot manage to obtain again a false color view with any combination of these pansharpened bands. What am I doing wrong? Are pansh.red the "pansharpened version" of spotimage.2, pansh.green the pansharpened version of spotimage.3 and pansh.blue the pansharpened version of spotimage.1? Remember thar the bands order is inverted. Should I use `i.landsat.rgb`? Can I use it with a SPOT image? Any help? Thanks!


title: QGIS - Calling interpolation extension from python console

text: I would like to call the QGIS interpolation extension function (TIN method) (Raster->Interpolate) from python console. I can't find the corresponding function within QGIS API or within the processing algorithm list. I found the SAGA Triangulation algorithm, wich work fine but is 5-10 x slower and speed is important in my case. Any idea how to execute it? thank you Vincent


title: Geometry type invalid occurs while trying to view data from oracle to arcmap connection

text: I have been trying to connect oracle databse in to arcmap 10.1. It got connected but while trying to add feature type ( route network file) to the map it says "Coodinate system identifier is invalid". The tables which are not associated with any geometry are displayed but the geometry related files are not showing up. I would really appreciate on any comment and thoughts.


title: Insert geometric data to a column from coordinate values in double precision

text: Screenshot of database http://wp.me/aIW9p-i0 For table “households”, I want to create two new columns: 1) “Geom”, data type – geometry; 2) “WKT”, data type ‐ text; I am intending to do the following tasks: 1. For each household, I have the location, indicated by x_coord and y_coord, respectively in the table. Then I have to create a geometric object – point, based on the x_coord and y_coord for each household. Insert these geometric data to the column “Geom” accordingly.   1. Now, after the first step, I have a geometric object for each household. Translate the geometric point data to the well‐known‐text,showing the longitude and latitude of the location of each house. Insert these values to the second column “WKT”, accordingly For the Step 1 I wrote the following function but I am getting the error: **Query returned successfully with no result in 54 ms.** CREATE OR REPLACE FUNCTION my_geom() RETURNS SETOF clean.households AS $BODY$ DECLARE x clean.households%ROWTYPE; pst geometry; BEGIN FOR x IN SELECT * FROM clean.households LOOP pst := ST_GeomFromText('POINT(' || x_coord || ' ' || y_coord || ')'); UPDATE clean.households set geom = pst; RETURN NEXT x; END LOOP; RETURN; END $BODY$ LANGUAGE 'plpgsql' ;


title: Need spatial data for very specific administrative boundaries in Ontario

text: I was wondering if anyone might know where I can download the rectangular administrative boundaries shown in the image below. I have highlighted one of them in red. On the image, you can see that the naming convention consists of some sort of numeric administrative identifier, as well as a name for each area (which seems to correlate roughly to the nearest city or town of any significance). I require the spatial file in order to provide a quote for a large project and I haven't been able to find anything. The only information the client has is in paper format. The title on the overview map shown below is **Provincial Series and Crown Land Use Atlas Index - 1:100,000**. I have searched the Provincial data available from LIO, the Federal data available at GeoGratis, and I've tried numerous Google searches with no luck. Any assistance would be greatly appreciated as I'd like to avoid digitizing if possible. Thanks in advance, and here's the image: ![enter image description here](http://i.stack.imgur.com/FwPuz.jpg)


title: Is there an interface that provides access to a property/method that shows whether a Layer is visible?

text: Repost--I deleted this yesterday when I discovered the answer and remembered this on meta. Had second thoughts this morning, since this might be useful to someone in the future. I'll let this go for a week and provide my own answer if it's unanswered. * * * I came across some hack code that I wrote from a few years back that I need rewrite in c#. At the time, I couldn't find property/method that would tell me if a layer was actually visible on the map. As opposed to the ILayer.Visible property which is the potential for a layer to be visible. IMapLayers.IsLayerVisible looked promising but it didn't return the expected values. My ESRI incident didn't have a related bug and the resolution to that incident was my own hack that I provided. How are others obtaining this information about a layer's true visibility within a map? **Update**   When I said true visibility, I meant that if a layer was visible and it's parent/grandparent/etc are visible too. 'True visibility' was a poor choice of words on my part, because it could be interpreted as a visible layer that has features within the current extent (which is not what I need).


title: How do I set the timezone in GeoServer w/ Jetty on Windows?

text: GeoServer is set to the local server timezone, but I need to have it use UTC. I believe that to make this happen I need to set the Java's user.timezone variable. My understanding is that if user.timezone is not set, Java (and so GeoServer) will use the server's local timezone. I'm running it with Jetty on Windows, talking to SQL Server 2008. We have a field Following these instructions, I tried creating a start.ini and putting               --exec     -Duser.timezone="GMT"      in it, but that file doesn't seem to be used. I next tried adding               wrapper.java.additional.4=-Duser.timezone=GMT      to the GeoServer 2.2.4\wrapper\wrapper.conf file, but that didn't seem to work either. I'm not sure what to try next and my searches are coming up empty.


title: How can I get the set of points (cloud) that are a given distance in SpatiaLite?

text: I have a set of points and need to get the groups of points that are within a distance x. Ie, determining the point clouds that are up to the distance x from each other. Any help would be appreciated. Thanks. Additional information: Example of a table with coordinates (in EPSG:4674) ![Table with points EPSG:4674](http://i.stack.imgur.com/V8LOq.jpg) The extension of the data is variable. Will be implemented in several counties, separately. I need anything like this: ![clustering buffer](http://i.stack.imgur.com/FZJLG.jpg) And I need make this in SQL from SpatiaLite. I can make this using buffer and union but I get just one polygon because I don´t have a "group by" attribute, and I need polygons separately for each group. I´ve been trying this:               CREATE TABLE "T_bufEndR" AS            SELECT "PK_UID" as "ID_end",                  ST_BUFFER(geom_ce, 0.00025) AS "geom_bu"             FROM "T_CoordEnd"            WHERE "situacao" IN (4, 5, 6, 7, 8)      and               CREATE TABLE "T_unbufEndR" AS           SELECT ST_UNION(b.geom_bu, a.geom_bu) AS "geom_bu"             FROM "T_bufEndR" AS "b", "T_bufEndR" AS "a"            WHERE ST_OVERLAPS(b.geom_bu, a.geom_bu) = 1 AND                  "b"."ROWID" IN (SELECT rowid from SpatialIndex                                    WHERE f_table_name = "t_bufendr" and                                         search_frame = "a"."geom_bu")      Is that!


title: OpenLayers with Yandex layer position issue

text: I am kind of new for openlayers and I am working on a project that have many different layers. one of this layer is Yandex, I added this layer few days ago and it seems to work fine, untill QA :) discovered that when you switch from this layer to another it not position the same as the other layer, I added this layer base on the following article (couldn't find anything else that show how to add the Yandex) http://xbb.uz/openlayers/i-Yandex.Maps when I have lines layer on the map and change the base layer from google, bing or any other layer to yandex the lines or the yandex layer are not positioned right. YANDEX layer ![enter image description here](http://i.stack.imgur.com/eEq5I.png) OPEN STREET MAP LAYER ![enter image description here](http://i.stack.imgur.com/7nVSQ.png)


title: Using Quickbird to discriminate vegetation response and spectral signatures - which season?

text: I have a study area in N. Florida loblolly pine plantation with 3 treatments (drought, fertilization, drought+fertilization, control) and 4 replicate plots. All plots will be in the same scene. These treatments are expected to alter the growth and response of pine trees. The only available existing imagery at the initial time of treatment installation is from quickbird (one in March and one in December) - this will be the baseline image. Ideally, it would be nice to acquire imagery every year and season, it's not feasible due to budget constraints. So, I will probably purchase one, maybe two more images over the next 1-2 years. They should be the same season as the "baseline" scene. **My Question:** Which season, if any, is best to discriminate spectral signatures between treatments/plots to really tease out differences in growth due to treatments and estimate NDVI, LAI, and other vegetation indices? I would suspect winter, when everything is not as lush and green and somewhat stressed (relative to summer). What is your expert opinion? Thanks in advance.


title: Labelling a multipoygon scattered in different areas

text: I am using **mapnik** and **tilecache** to generate tiles and **openlayers** to view map in the browser. In one of my **shapefiles(multipolygon)** with administrative boundaries, I found that more than one records have their **polygons displayed in two areas in the map**. In the map, **label is placed in center position** of the two polygons. When I opened the shapefile in **Qgis** , I found that label is displayed for each of the two polygons correctly. How will this type of labelling be done in mapnik. Thanks in advance.


title: Is it necessary to use a third party service for serving this case?

text: I have been assigned a task to help build a mobile app, which only shows terrain view of map with a fixed (min/max) zoom level based on mobile phone position and moveable markers which are provided by our server in json format (vector with speed, directions, etc) . Now.. since i haven't a big experience dealing with GIS data, how could i serve the tiled images to client without having to query a 3rd party service? I only need the static images plus a custom json to client, should i build a server or use a 3rd party service(gmap,etc). Given the amount of expected load and the restrictions set by 3rd party services on limitation of storing or even caching data on local machines is what i try to avoid. Contribution to the 3rd party such as a Google map is something i personally dislike, OpenStreet map though is better from my point of view. **Example** I have the London Territory which it isn't required to maximize Zoom, need to serve its Terrain point of view, based on OSM data for highways/rivers etc, a special json is generated from server-side which both tiles+custom json are supposed to be delivered to client, in order to sync and present a response to clients action. Since it rather unlikely for the terrain view (tile) to change for like the next year, and i can't cache from 3rd parties, how can i serve this? * * * > **Conclusion** _just for the record_ Well after a review and research we decided to go with Leaflet at client side since it looks mature, build our tilesets with TileMill and serve PNGs. The application will use own tiles for rendering terrain view of major cities, with availability for zoom (at fixed levels) , if any user wants to view other cities the application will try to connect to server and start caching tiles that it receives, but in favor of NOT consuming all of phone memory or ram, it will limit its caching on <100MB (guess that is enough).


title: Opening specific area only

text: I have a Shp file loaded of a County boundary, but i would like to open another shp file but only show the lines etc within that area and automatically delete anything outside. Thank-you


title: Are there any good libraries for visually creating simple Google Maps?

text: I'm working on a simple webapp that lets the journalists I work with easily construct styled Google maps. In short, I need an interface that lets users place markers, set zoomlevel and choose a center lat/long — all without writing any JavaScript. I know Google Maps has a "My Maps" functionality that provides this already, but I'm needing that kind of functionality as a way of adding data to my app. ![Something like this...](http://i.stack.imgur.com/LxGrS.jpg) _Something like this..._ Anyone know of such a library? I'm currently using a series of text fields using LocationPicker for the markers, but I still need to be able to set center and zoom — and I can't shake the feeling there's probably already something out there that does this... Thanks! **Edit:** This seems close, if a little clunky... http://www.birdtheme.org/useful/v3tool.html


title: Calculating raster values from center to edges to represent normal distribution

text: Edit for simplification/clarification: What I want to do is take a raster and calculate cell values from center to edges such that the center has an arbitrarily determined high value which decreases down to near zero as you get closer to the edges. The values of said raster would have a Gaussian distribution, peaking in the center and tapering down near the edges. Any suggestions are welcome, conceptual or technical. I am most comfortable with ESRI products, but not an absolute requirement.


title: What type of remote sensing methodology is commonly used to locate and track the movement of woodland caribou herds?

text: I am most curious whether it is possible to use Landsat or other freely available imagery (but info on any other method would great as well) to locate and track the movement of woodland caribou herds over time. I only have some experience doing this type of work with forest depletions and vegetation analysis but am curious whether similar techniques exist for locating and tracking animal herds.


title: Automated attribute update geometric network

text: I have a Network Dataset of a storm water system containing two feature classes of points and lines. (All data is topologically sound with pits and pipes snapped and linework captured in direction of flow.) The point feature class contains an attribute "ASSET TYPE", the line feature class has two attribute fields which hold this value as well (ASSET_US and ASSET_DS) referring to the attached asset of the upstream and downstream pits. My question is, could anyone suggest a way to automate updates (through ESRI's ModelBuilder or Python or FME) to the ASSET_US and ASSET_DS fields in the pipe layer when the ASSET TYPE field has been changed in the pit layer?


title: How do I convert a CSV of latitude/longitude coordinates with timestamps to a GPX tracklog?

text: Converting from a GPX track to CSV is easy, but I'm looking to go the other way. I have a table of point locations logged with a handheld GPS, with timestamps. I need to convert them to a GPX file with a tracklog containing those points so they can be used with some other software to geotag some photos. I can convert them to GPX waypoints easily enough, but not to a GPX tracklog. Here's the header of my CSV:               Lon,Lat,Altitude,NewTrack,TimeStamp,Datum,InternalTime


title: ArcGIS API using maps from ArcGIS Online

text: I am new to the ArcGIS API for JavaScript but I would like to create an application which uses an existing map of ArcGIS Online. Creating an application where just a map is shown works fine. But now I would like to add the posibility to toggle the different layers. I researched about that but I am not sure how to implement that without having layers as map services. I have a free ArcGIS online account where I can upload shapefiles to the map. Can I access those layers from the API?


title: Qgis Server, show Oracle Georaster dynamically

text: I have an Oracle Spatial database and I use rasters (SDO_GEORASTER). I can render some rasters on a web page with Qgis Server : when I render the maps on the webpage, it come from a qgis project (*.qgs) 2-3 new rasters are automatically added in the database each day. I want to create a web page to render the rasters dynamically : the user choose the raster he wants to show in a combobox. Is it possible to connect Qgis server directly to the Oracle spatial database ? Sorry for my english Thanks


title: How to combine Boolean and Fuzzy logic?

text: So I have a set of 5 Fuzzy Memberships and a polygon feature class, to which I already added a yes/no field. For the final map, I need to combine/overlay all Fuzzy Memberships and the rasterized polygon feature class in order to map suitability, and don't know how to combine it. Any advice?


title: How to disperse clusters of proportional markers or polygons?

text: I am mapping a set of points in ArcMap that are symbolized using proportional marker symbols based on a numeric field value. The scale of the map has lead to alot of the point markers overlapping. ![current cluster](http://i.stack.imgur.com/U5HYF.png) I would like to disperse the markers so that they are not overlapping. A random dispersement would be fine, preferably with control over the spacing between the dispersed markers. It seems that proportional symbols can't be converted to representations (else I could have used the Cartographic Refinement - Disperse Markers tool). I could buffer the points by the numeric field to create polygons instead of using point markers to get around the representation issue, but is there a way to disperse polygons away from each other so they are not overlapping (similar to the functionality of the Disperse Markers tool for points)? I'm using ArcMap with an Advanced License. Thanks


title: Why can´t Geowebcache communicate with geoserver?

text: I have installed geoserver 2.5 with integrated GWC on a windows server 2008 R2 in tomcat6. I am not able to create map caches with GWC. GWC gives me "404: cannot communicate with server. I have a reverse proxy set up and a url rewrite in IIS to provide internet access. What do I have to do to establish communication again between Geoserver and GWC?


title: Cause of MODIS Reprojection Tool errors?

text: I have installed the **MODIS Reprojection Tool** (mrt) and it says the installation completed successfully but it does not mosaic my images properly. However, there seems to be some projection errors. I want to open the file in ArcGIS for Desktop. Does anyone know the cause of the MODIS Reprojection Tool errors below? ![enter image description here](http://i.stack.imgur.com/B1Zeq.png)


title: Running x64 bit geoprocessing Python in Python IDLE 2.7.6 x64 error

text: I am trying to set the Python IDLE 2.7.6 x64 to run the arcpy module from ArcGIS x64 bit geoprocessing (10.2.2). I copied the file DTBGGP64.pth from C:\Python27\ArcGISx6410.2\Lib\site-packages to C:\Python27\Lib\site-packages that is written as:               C:\Program Files (x86)\ArcGIS\Desktop10.2\bin64     C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy     C:\Program Files (x86)\ArcGIS\Desktop10.2\ArcToolbox\Scripts      when I then open the python 2.7.6 x64 IDLE and run the following command : import Arcpy, I get this error               Traceback (most recent call last):       File "<pyshell#0>", line 1, in <module>         import arcpy       File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\__init__.py", line 24,     in <module>         from arcpy.toolbox import *       File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\toolbox.py", line 356,     in <module>         from management import Graph, GraphTemplate       File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\management.py", line 22,     in <module>         import _management       File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\_management.py", line 14,     in <module>         import _graph       File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\_graph.py", line 27, in     <module>         import numpy      ImportError: No module named numpy      it seems to me that the 10.2.2 x64 bit python is pointing to the 32 bit python library. I have tried substituting the code from: http://resources.arcgis.com/en/help/...00000008000000               c:\Program Files\ArcGIS\Desktop10.2.2\arcpy     c:\Program Files\ArcGIS\Desktop10.2.2\bin     c:\Program Files\ArcGIS\Desktop10.2.2\ArcToolbox\Scripts      I get the error that it cannot even find the arcpy module. My ultimate goal is to get the arcpy module to run in spss 22 which is built on the Python 2.7.1 x64 framework. However even when I put the DTBGGP64.pth file in the site- packages folder of that module, I get the same exact traceback error that it cannot find Numpy Any thought on how to solve this? Thanks in advance!


title: How to access the datasource from a REST service

text: I have a ArcGIS REST service, is it possible to access the data source of the layer that is within the rest service. For example I have a layer called "Conservation Area" within the service, that references the dataset "PLANNING.CONSERVATION_AREAS". It is "PLANNING.CONSERVATION_AREAS" that i want to retrieve back. Is this possible **EDIT** I want to get access to the actual feature class name which is held within the database (oracle in my case). I have read about the Data Source object REST Data Source Object which does suggest it is possible ? (but not sure how). Thanks


title: How wrap GRIB data in Mapserver 6.4

text: I'm trying to draw a combination of layers with Mapserver. One of the layer takes a shapefile as input and the other one is based on GRIB data. The issue is that the GRIB layer (in green) is drawn from 0 deg to 360 deg. The shapelayers draws from -180 to 180. And indeed, when I examine the GRIB file in Panoply, the meta data of sample shows, among other things, this:               netcdf file:/D:/gribs/sample.grb {       dimensions:         lon = 360;         lat = 181;         isobaric = 4;         time = 3;       variables:         float lat(lat=181);           :units = "degrees_north";              float lon(lon=360);           :units = "degrees_east";            ...       other variables       ...       }      Which leads me to conclude that the coordinate system runs from 0 to 360 degrees in the east ward direction and from 0 to 180 in the north direction. So Mapserver is correct when it draws this image: ![incorrect drawn grib layer](http://i.stack.imgur.com/Omxkw.png) **So my question is how can I wrap the GRIB image so that the right half is placed over the left half of the earth?** The url with parameters i use: http://mapserver-dev.de.meteogroup.net/cgi- bin/mapserv?map=/var/www/mapserver/routeguard/poc_test.map&mode=map&layer=grib_sample&layer=world_outline_red My mapfile:                   MAP         EXTENT      -180 -90 360 90 #max x should be only 180, but then only half of the grib layer is drawn         CONFIG      "CPL_DEBUG" "ON"         CONFIG      "MS_ERRORFILE" "/var/www/mapserver/routeguard/poc.log"         SYMBOLSET   "symbols.sym"         FONTSET     "../fonts/fonts.txt"         IMAGETYPE   PNG8         IMAGECOLOR  255 255 255         SIZE        1600 900          # ------------------------------------------------------------------------------     # the Map Projection     # ------------------------------------------------------------------------------         PROJECTION             #which to use????             "init=epsg:4326" #54004 #4326             "units=m"         END          # ------------------------------------------------------------------------------     # the Web object     # ------------------------------------------------------------------------------              WEB             IMAGEURL    '/mapserver/routeguard/' #do not include the the var/www             TEMPLATE    'poc_template.html'         END          #--------------------------------------------------------------------------------     # the country mask     #--------------------------------------------------------------------------------         #red world, defined in shp file         LAYER             NAME "world_outline_red"             TYPE POLYGON             DATA "/var/www/mapserver/data/routeguard/static/world_outline.shp"             CLASS                 COLOR 255 0 0             END                  PROJECTION                 "init=epsg:4326"             END         END          #----------------------------------------------------------------------------     #grib example     #----------------------------------------------------------------------------         LAYER             NAME "grib_sample"             TYPE LINE             CONNECTIONTYPE CONTOUR             DATA "/var/www/mapserver/data/grib/routeguard/sample.grb"                  PROCESSING "BANDS=2"             CLASS                 STYLE                     WIDTH 1                     COLOR 0 100 0                 END             END                  PROJECTION                 "init=epsg:4326"             END             METADATA                 #"wms_srs"       "epsg:54004"             END              END               END #end map file


title: What does this error mean- 'Underlying Connection was closed'

text: I am receiving this error message on one of our server while publishing a mobile project in mobile project center to ArcGIS server (10.1). ![enter image description here](http://i.stack.imgur.com/qurHj.png) Please advise - what connection is this one referring and how to resolve this? Thanks


title: How to check that a point is outside of a route?

text: I am newbie in PostGIS and also GIS particularly. I need to check that a vehicle is outside of a route. This route has been defined in advance, and of course it stores in a table: RROUTEDETAIL. This table contains: ROUTE_ID, LONGITUDE and LATITUDE. So, when vehicle leaves the route they should be on, the system will give a warning message and then record it into other table. I imagine to check the point I get with every point store in the RROUTEDETAIL, one by one with ST_DWITHIN. If all check get FALSE value, then the system give a warning message. Is this the best approach or is there a better idea ? Thank you very much. Best regards, Dion. * * * thank you. For example, this is the point where vehicle leaves: -7.539079, 109.322891 with value tolerance: 0.5. Then I write:                   select st_dwithin((select st_makeline(ST_Point(-7.539079, 109.322891))),             (select st_makeline(ST_Point(foo."TRACKERDATA_LATITUDE",foo."TRACKERDATA_LONGITUDE"))                 from (select ht."TRACKERDATA_LATITUDE",                              ht."TRACKERDATA_LONGITUDE"                             from "HTRACKERDATA" ht                             left join "RVEHICLE" rv on ht."VEHICLE_ID" = rv."VEHICLE_ID"                             where ht."VEHICLE_ID" = 106 and                              ht."TRACKERDATA_GPSTIME" BETWEEN '2012-01-01 08:00'                             AND '2012-01-01 17:00'                             order by ht."TRACKERDATA_GPSTIME", ht."TRACKERDATA_REPORTTIME") as foo),              0.5);      It always return FALSE. Then I check again with ST_DISTANCE, it gives 99999999999.9 (double precision). How to fix that ? Thank you.


title: How do I use polygon centroids within a GeoQuerySet?

text: I have two models: neighbourhoods and boroughs, each containing a polygon field. I am trying to construct a QuerySet that will give me the neighbourhoods within a given borough. I expected to use _within_ or _contained_ : `neighbourhoods_in_borough = Neighbourhood.objects.filter(poly__within=Borough.poly)` However, it is not always the case that a neighbourhood lies cleanly within a borough. Neighbourhood polygons roughly align with a Borough, but may weave in and out slightly. As such, filters like within or intersects either leave out neighbourhoods or include too many surrounding a Borough which happen to slightly overlap. Because I only want the Neighbourhoods which are predominantly inside of a given Borough, I thought I would test to find Neighbourhoods whose centroids are within a given Borough: `neighbourhoods_in_borough = Neighbourhood.objects.filter(centroid__within=Borough.poly)` However, you can spot the issue here: centroid is not a field on the model, and so the above does not work. I did think adding a preceding query (e.g. Neighbourhood.objects.centroid() ) would make the centroid available to be used in some form, but it fails for the same reason as the example above. Is it possible to filter based on a polygon's centroid rather than the polygon's shape?


title: How to print pie charts in QGIS?

text: I added pie charts via text-file to my map. When printing the map they get lost. Does anyone know how to print them?   QGIS also seems to delet some charts, when there are to many at one place. Does anyone have a solution for that?


title: How can I use polygon grid to export raster

text: I have a polygon grid which consists of 100 square polygons. For each square I need to export the raster overlay. I have ArcGIS 10.0, FME 2014, Saga-GIS. Is there a suitable tool for this or do I need to script something with Python or Model builder?


title: Idiots guide to expressions for QGIS

text: I have no formal training in GIS and have learned on the fly with ArcGIS 9.3 at work and have been dabbling with QGIS on my own time for my business. Is there any resource recommendations ingeneral for getting past the basics of GIS and going into the deep dark secrets? Specifically in relation to experession coding for labeling, queries, etc.


title: Calculating area of MultiPolygons within 250m buffer of points

text: I am using using PostgreSQL 9.2 and PostGIS 2.0. I have two tables: 1) a geocoded table of home addresses (address_list) and 2) a table containing data on tree density in the same geographical area. I have included examples of the relevant columns in each table below. address_list:               house_id (serial) |  geom(Point, 4283)     ------------------+-----------------------             1         |              2         |      tree_density:                 gid (serial)  |  tree_dens |  geom(MultiPolygon, 4283)         ------------+------------+----------------------------               1     |  scattered |               2     |  medium    |               3     |  dense     |      What I would like to do is compute the area of forest (m2) by tree_density categories within 250m of each home address for my whole address table (approx 6,000 entries). (E.g for house_id = 1: area_scattered_within250m = 50m2, area_med_within250m = 0m2, area_dense_within250m = 20m2) Can this be done easily? Even if in 3 or more seperate queries? Does anyone have a moment to assist me with the syntax? I guess I need to combine ST_Area, ST_Buffer and perhaps ST_Contains or St_Within, but I'm new to PostGIS and although I've been ok with other (simpler) spatial joins, I'm lost here and haven't found any similar examples. For the area of my study I transform to SRID 3112 (meters) for measurements. Many thanks Mike


title: QGIS connection to PostGIS showing wrong TYPE

text: I've got a user that it trying to edit his polygon layer but gets a message `PostGIS error while adding features: ERROR: Geometry type (Polygon) does not match column type (MultiPolygon)` When you have a look at his PostGIS connection details, the TYPE reads as polygon but if I log as an administrator into PostGIS via QGIS and edit the same layer i dont get an error message but see the TYPE as multipolygon. I'm not sure where the error lies, in PostGIS or with QGIS not updating itself!


title: Installing pandas with Esri's python?

text: Trying to get pandas to import, but I'm getting what looks like an error caused by a failed dependency of the version of my Esri-provided numpy's compile, based on this stackoverflow posting. Here's what I get.               >>> import pandas     numpy.dtype has the wrong size, try recompiling     Traceback (most recent call last):       File "<interactive input>", line 1, in <module>       File "C:\Python27\ArcGIS10.2\lib\site-packages\pandas\__init__.py", line 6, in <module>         from . import hashtable, tslib, lib       File "numpy.pxd", line 157, in init pandas.hashtable (pandas\hashtable.c:22315)     ValueError: numpy.dtype has the wrong size, try recompiling      I check my numpy version:               >>> import numpy     >>> numpy.version.full_version     '1.6.1'      Based on this, I believe that this is the right version of numpy for the version of pandas I'm using (tried 13.1, 14.0, and 14.1). I'm running the standard ESRI-shipped python for ArcGIS 10.2 (Python 2.7.3 (default, Apr 10 2012, 23:31:26) [MSC v.1500 32 bit (Intel)] on win32).               >>> import arcpy     >>> arcpy.__file__     'C:\\Program Files (x86)\\ArcGIS\\Desktop10.2\\arcpy\\arcpy\\__init__.py'      Has anyone succeeded in getting pandas going this way or have ideas on what's wrong with my environment? Should I update my numpy, which is supposed to be backwards compatible (according to user333700's text at the bottom of the stackoverflow posting), or is that likely to break my ArcGIS functionality? Thanks for any info!


title: How to populate the geometry_columns table

text: I have a SQL Server table of vector data (containing latitude, longitude, UComp and VComp columns), that I am trying to add as a layer to QGIS. Here is a sample of the data in the SQL Server data table: ![enter image description here](http://i.stack.imgur.com/28a8o.jpg) And the table design: ![enter image description here](http://i.stack.imgur.com/rAjEb.jpg) I have successfully connected to my database from QGIS and also created the geometry columns table. What data should I be adding to the [f_geometry_column], [coord_dimension] and [geometry_type] columns in the geometry_columns table?


title: Where to execute sdetable command in Windows?

text: I am trying to fire some Sde Command in Arcgis Desktop 10.1 but I am not aware where to execute the command in windows operating system.Thank you in advance.


title: Why Do Colors Look Darker When Printed on Paper?

text: I made a map with a terrain-colored relief-shading in the background. I then exported them as PDF from Print Composer in QGIS. But **when I print the PDF onto paper:**   * the shaded-relief color in the **background appeared _way_ darker** than it supposed to.   * the colors of roads, symbols and labels in the **foreground looked all right** , though..  All of the colors in the PDF itself appeared correct. **Is there a way to ensure my print colors look right on paper?** Appreciate any tips/documents/links that can help me get my prints right, thanks!


title: Is there a "Select by Attributes" plugin/widget for ArcGIS Server JS API?

text: I'm hoping to find some existing code which emulates the "Select by Attributes" dialog in ArcMap, within the ArcGIS Server JavaScript API. Writing a hard-coded QueryTask is easy - I'm looking for some code which will automatically handle any layer on the map, forming a Query on-the-fly using any attributes, ideally returning the results into a datagrid. Clever people have written Table of Contents and Identify widgets which handle any layers in the map. Has anyone written something similar for "select by attributes"? Thanks


title: GeoWebCache tiles with specific scales

text: Is it possible to generate some tiles by GeoWebCache with specific scales eg. 1 000 000 , 500 000 etc. By default they are eg. 611.4962261962891, 305.74811309814453, 152.87405654907226, 76.43702827453613 what is very confusing for me actually ;/ Thanks for any help Alex


title: What's a term for Sea Level Distance which is usable by novice users?

text: I've been using the term Sea Level Distance to describe the total length of a line string where I ignore the elevation and just calculate the distance from point to point along the ellipsoid and sum the total. Sea Level Distance has always seemed to be an artificial term for this. Slope Distance seems to naturally describe the distance when elevation is included. Is there a term for Sea Level Distance that is more natural and usable for novice users.


title: VBA Convert decimal degree to Meters button

text: I have 2 data frames one which has WGS data in it and the other OSGB1936. I am trying to create a button where the user can input decimal degrees (WGS) and it will zoom them to that location in British National GRid (OSGB1936) I have done the code but when it zooms to OSGB data frame it does not convert the decimal degree to meters. It returns incorrect values for Northing and Eastings in meters. The code does work only if i have done it as a ToolControl_MouseDown but that uses pMxDoc.CurrentLocation.x, pMxDoc.CurrentLocation.y. which i dont want to do because i want the user to enter the value from a text box. Any ideas please?               Private Sub CommandButton2_Click()               Dim pSREnv As ISpatialReferenceFactory2     Set pSREnv = New SpatialReferenceEnvironment     Set m_PCS = pSREnv.CreateProjectedCoordinateSystem(esriSRProjCS_BritishNationalGrid)          Dim pMxDoc As IMxDocument     Set pMxDoc = ThisDocument          Dim x As Double     Dim y As Double           y = CDbl(Texty.Text)      x = CDbl(Textx.Text)          Dim pEnvelope As IEnvelope     Dim pMap As IMap     'Set pMap = pMxDox.FocusMap          Dim strName As String         strName = "Layers"          Set pMap = GetMap(strName)      If Not pMap Is Nothing Then         If pMxDoc.ActiveView Is pMxDoc.PageLayout Then             Set pMxDoc.ActiveView.FocusMap = pMap         Else             Set pMxDoc.ActiveView = pMap         End If              End If              Dim pPoint As IPoint         Set pPoint = New Point                   pPoint.PutCoords x, y     'pMxDoc.CurrentLocation.x, pMxDoc.CurrentLocation.y              Dim pGeom As IGeometry         Set pGeom = pPoint         Set pGeom.SpatialReference = pMxDoc.FocusMap.SpatialReference         pGeom.Project m_PCS                   Set pEnvelope = pMxDoc.ActiveView.Extent         pEnvelope.CenterAt pPoint         pMxDoc.ActiveView.Extent = pEnvelope         pMap.mapScale = 5000         pMxDoc.ActiveView.Refresh                       Application.StatusBar.Message(0) = "meters east: " & pPoint.x & ", meters north:" &     pPoint.y          End Sub      * * * The below code is the updated but still does not do what it should. When i enter the decimal degrees 53.47574, -2.22553 it returns a value of -512464, -2174656 which is in the incorrect position. So somewhere in the conversion between WSG and OSGB it is going very wrong.               `Private Sub CommandButton2_Click()               Dim pMxDoc As IMxDocument     Set pMxDoc = ThisDocument          Dim pEnvelope As IEnvelope     Dim pMap As IMap          Dim pSREnv As ISpatialReferenceFactory2     Set pSREnv = New SpatialReferenceEnvironment          Dim pProjCoordSystem As ISpatialReference      Dim pProjCoordSystem2 As ISpatialReference     Dim pGeoTrans_A As IGeoTransformation     'Dim pGeoTrans_B As IGeoTransformation               Set pProjCoordSystem2 =   pSREnv.CreateProjectedCoordinateSystem(esriSRProjCS_BritishNationalGrid)     Set pProjCoordSystem = pSREnv.CreateProjectedCoordinateSystem(esriSRProjCS_WGS1984UTM_14N)                    Set pGeoTrans_A = _     pSREnv.CreateGeoTransformation(esriSRGeoTransformationType.esriSRGeoTransformation_OSGB1936_To_WGS1984_1)                Dim strName As String         strName = "Layers"           Set pMap = GetMap(strName)      If Not pMap Is Nothing Then         If pMxDoc.ActiveView Is pMxDoc.PageLayout Then             Set pMxDoc.ActiveView.FocusMap = pMap         Else             Set pMxDoc.ActiveView = pMap         End If              End If               Dim x As Double     Dim y As Double                y = CDbl(Texty.Text)      x = CDbl(Textx.Text)                  Dim pPoint As IPoint         Set pPoint = New Point         'Set pPoint = New esriGeometry.Point         Set pPoint.SpatialReference = pProjCoordSystem2                   pPoint.PutCoords x, y     'pMxDoc.CurrentLocation.x, pMxDoc.CurrentLocation.y             Dim pGeom As IGeometry2         Set pGeom = pPoint         Set pGeom.SpatialReference = pProjCoordSystem2         pGeom.ProjectEx pProjCoordSystem2, esriTransformForward, pGeoTrans_A, 0, 0, 0                   Set pEnvelope = pMxDoc.ActiveView.Extent         pEnvelope.CenterAt pPoint         pMxDoc.ActiveView.Extent = pEnvelope         pMap.mapScale = 5000         pMxDoc.ActiveView.Refresh                   Application.StatusBar.Message(0) = "meters east: " & pPoint.x & ", meters north:" & pPoint.y             End Sub'


title: How to import a non-spatial Table into QGIS?

text: Can I import tabular data in Qgis, without any information about X and Y only the table? I need to join them after importing them.


title: How can I apply "natural breaks" classification to a raster symbology in QGIS?

text: Anyone knows how can I apply "natural breaks" to a raster symbology in QGIS? Cause I only see the "Equal interval" option.. ![enter image description here](http://i.stack.imgur.com/v49Jr.png) I'm using QGIS 1.9


title: How to perform SQL query using .STContains in MS SQL Server

text: I have a table that stores polygons and a table that stores points with a value field. Both sets of data are stored using the Geometry type in MS SQL Server. I have a simple requirement - I need to retrieve the records of points that are contained within the polygons and then sum the result based on a value field in the points table. So, the end result of the query is two columns: - polygon names - sum of the value field in the point table               DECLARE @area GEOMETRY  = (SELECT Shape FROM GREENFIELDAREAS WHERE AREANAME = 'Aerodrome')     SELECT NETDWELLINGS, Shape FROM ORA_NETDWELLINGCOMPLETIONS WHERE @area.STIntersects(Shape) = 1      The above query returns point geometries but only for 1 polygon (Aerodrome). What I need is to sum NETDWELLINGS and return that value for all the polygons.


title: how to easily serve vector layer with extra attributes for openlayers

text: On a website I'd like to show an openstreetmap with buildings colored by a certain feature (tag). So far I was able to extract the buildings for a bounding box, save them as an osm file and load it via an openlayers vector layer. Next I'd like to load the buildings dynamically and based on the current bounding box. I could probably achieve that using overpass api and data from openstreetmap.org etc. However in the future I'd like to add a custom tag - e.g. sshpa (sunshine hours per year) - to each building and conditionally style the buildings according to the tag. I could probably write my own server returning vector data to openlayers upon ajax calls based on the bounding box, but I wonder if there already exists a solution, that does this for me. After all this is probably a standard problem for osm users. I am looking for an application that serves vector data to openlayers.js and where I can easily manipulate the data of the layers, e.g. by extending them with tags.


title: REtrieve field info from web map services in ArcMap

text: I am trying to work out what fields are associated with a map service in ArcMap. The service is very slow in drawing and I wondered if it was possible to retrieve attributes associated with a layer another way?


title: Broken PostGIS support GDAL

text: I'm having some problems with Postgres/PostGIS after rebuilding my GDAL library:( I wanted to add MrSID support in GDAL so rebuilt it with the MrSID SDK. Everything compiled ok and I am now able to open Raster .sid files within QGIS. BUT I seem to have screwed my Postgres 9.3 install as now I have no PostGIS support. If I try and open previously functioning postgis databases in QGIS I get a ton of error messages indicating that my database has no working PostGIS support. When I look at my databases in PGADMIN3 I get the following error:` ERROR: could not load library "/usr/lib/postgresql/9.3/lib/postgis-2.1.so": /usr/lib/postgresql/9.3/lib/postgis-2.1.so: undefined symbol: GEOSHausdorffDistance` The postgis-2.1.so library is definitely present on my system and was working fie before I recompiled GDAL. I'm running Ubuntu 13.10, Postgres 9.3, PostGIS 2.1, and GDAL 1.10. This is how I compiled GDAL 1.10:               ./configure --with-geos=yes --with-sqlite3 --with-python --with-spatialite --with-fgdb=/usr/local/FileGDB_API --with-postgres --with-odbc --with-curl=yes --with-mrsid=/usr/local/MrSID_DSDK/Raster_DSDK --with-mrsid_lidar=/usr/local/MrSID_DSDK/Lidar_DSDK --with-poppler=yes --with-jasper=yes --with-opencl=no --with-ogdi=no      Anyone have any idea how to fix this? Do I need to recompile Postgis against Gdal?! Really appreciate any help with this for as usual I urgently need to use my postgis data! Thanks Marty _**EDIT_** Looking at libraries for postgis-2.1.so gives:               $ldd /usr/lib/postgresql/9.3/lib/postgis-2.1.so     linux-vdso.so.1 =>  (0x00007fffcb1fe000)     libgeos_c.so.1 => /usr/local/MrSID_DSDK/Raster_DSDK/lib/libgeos_c.so.1 (0x00007f84d821e000)     libproj.so.0 => /usr/lib/libproj.so.0 (0x00007f84d7fcc000)     libjson-c.so.2 => /lib/x86_64-linux-gnu/libjson-c.so.2 (0x00007f84d7dc1000)     libxml2.so.2 => /usr/lib/x86_64-linux-gnu/libxml2.so.2 (0x00007f84d7a5b000)     libc.so.6 => /lib/x86_64-linux-gnu/libc.so.6 (0x00007f84d7693000)     libpthread.so.0 => /lib/x86_64-linux-gnu/libpthread.so.0 (0x00007f84d7475000)     libgeos.so.2 => /usr/local/MrSID_DSDK/Raster_DSDK/lib/libgeos.so.2 (0x00007f84d717d000)     libstdc++.so.6 => /usr/lib/x86_64-linux-gnu/libstdc++.so.6 (0x00007f84d6e79000)     libm.so.6 => /lib/x86_64-linux-gnu/libm.so.6 (0x00007f84d6b74000)     libgcc_s.so.1 => /lib/x86_64-linux-gnu/libgcc_s.so.1 (0x00007f84d695e000)     libdl.so.2 => /lib/x86_64-linux-gnu/libdl.so.2 (0x00007f84d675a000)     libz.so.1 => /lib/x86_64-linux-gnu/libz.so.1 (0x00007f84d6540000)     liblzma.so.5 => /lib/x86_64-linux-gnu/liblzma.so.5 (0x00007f84d631e000)     /lib64/ld-linux-x86-64.so.2 (0x00007f84d86f0000)


title: Model Builder - import arcinfo licence when required by tool

text: Is it possible to check licences in/out as needed within modelbuilder? or will the following work -write a script to turn arcinfo licence on when the tool requires it and off once it's done? I need to run a complex model on 100+ segments of data on multiple machines concurrently and we have only one arcinfo licence which is required only by the eliminate tool. Each segment takes about an hour which is approx 100hrs of processing time. We don't have this much time so need to split the process into 6 machines running the arcview licence...the blocker is one-step needing arcinfo. As the machines will need the licence for only 10mins or so at differeing times, I want to get the model to run in arcview upto the eliminate, check-out arcinfo (or if unavailable keep checking every 5 mins), finish the tool and then check-in the info licence so another machine can take it and continue. Thanks guys,


title: CORS Printing Issue with ArcGIS Javascript

text: I'm developing an application using the ArcGIS Javascript API and everything has been working fine until I've tried to incorporate the print widget using printtask. When only the base layers are accessed in the code, the map is printed just fine. But when I include any layer running on my local instance of ArcGIS Server 10.2 (dynamic or feature layers) the print operation fails. I assume this is because I need a "proxy" file although I thought I read it's possible to utilize CORS. I'm using the latest FireFox browser, ArcGIS Server 10.2 and I have configured my code to include:   esri.config.defaults.io.corsEnabledServers.push("10.111.123.234:6080"); esri.config.defaults.io.corsEnabledServers.push("sampleserver1a.arcgisonline.com"); esri.config.defaults.io.corsEnabledServers.push("sampleserver6.arcgisonline.com"); esri.config.defaults.io.corsEnabledServers.push("tasks.arcgisonline.com"); esri.config.defaults.io.corsEnabledServers.push("server.arcgisonline.com"); The IP address above is my local PC pointing to the ArcGIS Server port of my services. I have tried just localhost:6080 as well. Either way, the printing functionality still fails so I'm not sure if CORS is valid for this operation and that I have to use a proxy file anyway? Or perhaps I'm doing something wrong? Any thoughts or suggestions would be appreciated.


title: Merging parcels but keeping Geometry seperate

text: I am working on updating a parcel layer we received for a project, and I have duplicate parcelID's for parcels and I have been merging these to create 1 parcel as to not have duplicate ID's. My question is, when i look at the town online viewer some parcels have the same Parcel ID but they still maintain a separate geometry for these adjacent parcels. The information that is pertinent to my project is the same for each parcel, so I would like to merge them but maintain the separate property boundaries. This is an SDE Polygon feature class. P.S. I found a post that is pretty much the same question entitled How to Create Touching Multi-Part Polygons? How to Create Touching Multi-Part Polygons? Have there been any updates on this in the last year? Thank you.


title: Convert 800 .hdf4 files to Tiff using GRASS/QGIS/SAGA/GDAL/R

text: I've have MODIS ndvi data (MODIS_Grid_16 days NDVI) from the iberian peninsula, from February of 2000 to August of 2012 (really big data). How do I batch convert about 800 .hdf MODIS files to .tif ? I can use QGIS/GRASS/SAGA/GDAL/R but I have no experience with bash scripts.


title: Calculate polygon area in ArcMap

text: Please please please help me figure out what the deal is with my area units!!When calculate geometry for census tracts in km^2, I am given back some absurd, meaningless number such as 0.000696186


title: PostGIS - certain MultiPolygons cause “BOOM! Could not generate outside point!”

text: I've just posted this at StackOverflow and someone suggested I try here. Original: http://stackoverflow.com/q/5863859/246265 I'm trying to represent a rectangular area which crosses 180 degrees longitude. For more background see http://stackoverflow.com/q/5737506/246265 Here's my test case:               from django.contrib.gis.geos import Polygon, MultiPolygon     from my_project.my_app.models import Photo          a = Polygon.from_bbox((30, -80, 180, 80))  # the part to the east of 180     b = Polygon.from_bbox((-180, -80, -160, 80))  # a part to the west of 180     c = Polygon.from_bbox((-180, -80, -100, 80))  # a larger part to the west of 180          ok   = MultiPolygon(a,b)     ok2  = MultiPolygon(c)     boom = MultiPolygon(a,c)          # This works     Photo.objects.filter(location__coveredby=ok)[:1]     # This also works so c is ok     Photo.objects.filter(location__coveredby=ok2)[:1]     # This gives "BOOM! Could not generate outside point!"     Photo.objects.filter(location__coveredby=boom)[:1]          # splitting c doesn't help     c1 = Polygon.from_bbox((-180, -80, -140, 80))     c2 = Polygon.from_bbox((-140, -80, -100, 80))     test = MultiPolygon(a,c1,c2)     Photo.objects.filter(location__coveredby=test)[:1]     # BOOM! Could not generate outside point!      By changing the numbers I can make this error come and go. (-180, -80, x, 80) works where x <= -140 for example. For every number there is a threshold like this but I can't find a pattern. For boxes with the same area, some will work and others not. For boxes with the same width some will work and others not. I can look at the SQL being generated but the areas are represented in binary (EWKB) and I'm not sure how to read it. Can anyone explain this?