title: When creating a new raster, how to determine a raster's upperleftx, upperlefty?

text: I'd like to create a raster from point data, averaging the point data over a given area and placing the results into a given pixel/cell. With each pixel/cell as 50m x 50m. I'm new to this and have a beginner question. I see that in PostGIS 2.0 there is a ST_MakeEmptyRaster() function, but what are the units for **width** , **height** , and **pixelsize**. And how do you determine the **upperleftx** , **upperlefty** points? Is this just random, or is there a convention used to determine this point?               raster ST_MakeEmptyRaster(integer width, integer height, float8 upperleftx, float8 upperlefty, float8 pixelsize);


title: Android GIS - Reverse Geocode using Shapefile (.shp)

text: I want to develop an Android application. I'm still a newbie here. What I want is to develop my app just like a GIS that can implement reverse geocoding. So I can put in lat / lon coordinates and then text is returned, for example the city name. But I don't want to implement Reverse Geocode using Geocoder from Google, all I want is for the coordinates that I entered to be checked from a shapefile (.shp) in the database, then return a city name. And I want to use MySQL as a database for my shapefile (shp). Is it possible, and what should I do to make it happen? Could you give a tutorial? Thanks in advance.


title: Qgis Python code library documenation.

text: I'm looking for python documentation for specific operations in the Qgis python console. I know there is no suggestions simular to Arcpy but I would like a pure list of operation code lines as "layer.featureCount()" Or is this taken from python standard librarys?


title: Creating an INFO table using ITable in Arc (Java)

text: I'm trying to set up an info table in an addin for ArcGIS 10. I'm using Java and Eclipse IDE. The code I have is here:               public class Tables {          String database = "C:\\GIS";     String nameOfTable = "getOn";     public ITable createInfoTable(String database, String nameOfTable) throws IOException  {                            IPropertySet propertySet = new PropertySet();             propertySet.setProperty("DATABASE", database);                            IWorkspaceFactory workspaceFactory = new ArcInfoWorkspaceFactory();                            IArcInfoWorkspace workspace = (IArcInfoWorkspace)workspaceFactory.open( propertySet, 0);                  //Fails here (or above actually)  Check out Factory.open                       IArcInfoItems arcInfoItems = new ArcInfoItems();             ITable  myTable = workspace.createInfoTable("myTable", arcInfoItems);                       return myTable;      I've adapted this code from C# as there are no examples on the Arc pages for Java. You can see with the comment where it is failing. If you can see what's going wrong please let me know. And I will repost if I crack it. (see here for original C# example: http://resources.esri.com/help/9.3/ArcGISEngine/ArcObjects/esriDataSourcesFile/IArcInfoWorkspace.CreateInfoTable_Example.htm)


title: Serving GroundOverlay image data from a KML file via WMS

text: I've got a KML file with a collection of sparsely distributed (and small) GroundOverlays, with each overlay containing a distinct LatLonQuad and Icon. I would like to serve the this image data via a WMS (ideally GeoServer). The Icons (images) in the KML are aerial photos of real-world objects, so it is imperative that they scale correctly according to the zoom level of the view. I was thinking of writing a script to convert the KML file to one or more GeoTIFF files or something similar (as described here: Export georeferenced image overlay from KML?), but it feels like that is not the most efficient approach, particularly considering that the GroundOverlays are so sparsely populated. Another option I was looking at was trying to just create a Shapefile from the KML file and use SLD with a dynamic symbolizer (as described here: http://stackoverflow.com/questions/2589181/display-image-over-map-layer-in- geoserver) to avoid creating a raster for such sparse data, but I suspect that the icon images will be tiled instead of stretched at different zoom levels, which is not the intention. Is there a quick and easy way to serve this kind of data via a WMS?


title: python code to create an intersection loop with a geodatabase

text: I am trying create a loop that will intersect each individual township with a statewide section layer. I am having trouble with the syntax of the intersect tool I think, but maybe it is something else. Any help is greatly appreciated.               # Import system modules     import sys, string, os, arcgisscripting     import arcpy, os, sys, string     from arcpy import env     from arcpy.sa import *     gp = arcgisscripting.create()          arcpy.OverWriteOutput = 1     arcpy.env.overwriteOutput = True               input_folder = "C:\\state\\towns.gdb"     arcpy.env.workspace = input_folder     out_wksp = "C:\\state\\township_sections.gdb\\"          env.workspace = input_folder          fcl = arcpy.ListFeatureClasses()     print fcl          infeatures = ["C:\\state\\towns.gdb\\sections", fcl]     print infeatures               #ExtractByMask_sa (in_raster, in_mask_data, out_raster)      for fc in input_folder:         print fc         arcpy.Intersect_analysis(infeatures, out_wksp + fcl, "ALL", "", "INPUT")         #arcpy.Intersect_analysis(["C:\\state\\towns.gdb\\sections" fcl], out_wksp + fcl, "ALL", "", "INPUT")     #    fc = fclist.Next()         print 'clipped', fcl          print 'done'


title: Find Selected Element

text: Using the esriControls.ControlsSelectTool I select a text element that I have previously placed on the map. In code, how can I find this element on the map? I know I can iterate using the IGraphicsContainer but how can I tell that the element is the selected element?


title: SQL Server Geometry to geography conversion

text: I have imported a number of Shape files from the ordnance survey using shape2SQL these have been imported as geometry however i would like to convert them to Geography and specifically to srid 4326 the shape when imported contains a number of polygons ultimately i'm trying to ascertain if a number of lat/long coordinates i hold fall within one of these polygons (so i suppose i also need to alter the coordinates to geog (which i can see has been answered before)) finally the plan is to display this on an SSRS report or onto a Bing Map display (hence why at this stage i believe i require SRID 4326) so i guess my question is, how do i go about importing the shape as SRID 4326


title: How to create a 50x50 m hexagonal grid with mmqgis?

text: I'm trying to use the mmqgis plugin to create a hexagonal grid, and I am really perplexed on how to fill in the variables required. I would like to have hexagons around 50m x 50m in size (real life), and my project is in long/lat. I hope that someone can help. Thanks Sal


title: Can Torque animate markers along polylines?

text: I found that torque is currently designed to animate markers along a table of points. I'd like to extend torque to animate markers along the polylines of a table's row. This would allow torque to animate more than one marker per table. Can anyone point me in the right direction?


title: How can i obtain the data of one selected feature in QGIS using python?

text: I have a problem / question regarding python and QGIS. Is it, and yes how, possible to get the data of one selected feature in QGIS using python? In my example I have a layer which shows some parcels for urban planning. If I now select one of the parcels I want to get the data for exactly this parcel. This data should be read out by an python script and then added to an automatically created PDF file. The last step is easy, but I didn´t get the data from this selected feature... It is possible to get data from the active layer, but is it also possible to get the data from one selected feature of this layer? Hope you can help me. Regards


title: Google Maps Ground Overlay: How to calculate bounds?

text: I'm working with the Google Maps API for the first time and I'm trying to append a simple image of a Building on top of my current map view. This seems to be very straight forward when I look at the GMapsAPI here … https://developers.google.com/maps/documentation/javascript/overlays#GroundOverlays However, I have no idea how to calculate the `bounds` and the exact position of my overlay. As you can see in the documentation the position of the ground overlay is done via `LatLng`               var imageBounds = new google.maps.LatLngBounds(         new google.maps.LatLng(40.716216,-74.213393),         new google.maps.LatLng(40.765641,-74.139235));      In my case I have no idea how to calculate those bounds. What are those values exactly? Why are there two seperate `LatLng` positions needed? I have a simple .png of a building that I want to layer on top of an exact position. I get this position by entering the destination in Google Maps and reading the `Lat` and `Lng` values from the embed-code. Is there a better way of doing this? Any ideas, tipps or tricks on how I could do this? Or is it possible to calculate those `bounds` when I know the center of the image-position and the dimesions of my image file? (My image is 50x50px and a .png) Thank you for your help. Kind Regards, Matt


title: Is PostGis suitable for Global data routing?

text: I'm searching quite some time now... Is it possible to load the global routing data into PostGis? What can I expect from a performance point of view? It's for routing purposes only. I've downloaded the global data, extracted it, saw the 300Gig size and started wondering... the .osm is 300Gig, but still has to be converted to pgRouting. Searching on the internet ends up mostly in examples whereas pgRouting is used for boundaries like cities... hence my question over here. Thanks for answering, any type of help is useful! Johan


title: Steps to creating a GeoTIFF

text: First off, I'd like to mention that I'm pretty new to working with GIS, so please forgive my inexperience. My situation is as follows: I'm tasked with georectifying a set of imagery and importing those images into a Google Earth Enterprise server as an overlay. From what I've researched, it sounds like a good route to do this would be to convert the images into a .tif, install GeoTIFF metadata, and import the image directly (Google Earth will figure out the rest). However, I'm not exactly sure where to go from there to actually create the GeoTIFF metadata (and install it). Here's the data I have: lat/lon of the center of the image, orientation (i.e. degrees to rotate to orient north), and resolution (e.g. 90m). Does anyone know of a good method/set of tools to derive any other data and install it into the .tif? Is something other than a GeoTIFF a good idea? Thanks for the help!


title: Transform SPOT5 images to natural color images

text: We're currently using PixelSense to convert raw SPOT5 images to natural color ones for display and printing. Do you know a good algorithm to perform the same task as we want to ingest raw images in an image server and provide an image service with on the fly natural color conversion.


title: How to add custom events to OpenLayers-based API?

text: I am woking on a GIS API built on top of OpenLayers. I have a class called xgis.group, which represent a group of layers. I want to add an event that can be triggered when a group is added. For example, in the map, I have tried this without success:               xgis.map.prototype.addGroup = function(name) {              /* CODE FOR ADDING A GROUP */              /* Trigger when group is added */         this.events.triggerEvent('groupadded', group);     };      I don't know how to deal with custom-made events. Could you lend me some help, please? Thanks


title: Geocode from mile marker/milepost to latitude/longitude?

text: I have a dataset that describes a few million incidents along Texas roadways. The vast majority of incidents only list a highway name and a milepost. For example, milepost 120 on I-20 or milepost 263 on US 60. How would I convert these milepost/highways to latitude and longitude? This is for a low-budget academic project, so trying to stick with free tools as much as possible. My Google searching hasn't turned up much fruit, although some indications are that TxDOT may have data I could request through the Texas Public Information Act, but not sure yet.


title: Editing layers in geoserver with web application

text: I wanted to edit layers in geoserver where they are store in PostGIS database. I have refereed http://www.gistutor.com/openlayers/22-advanced-openlayers- tutorials/47-openlayers-wfs-t-using-a-geoserver-hosted-postgis-layer.html link for that, which works. But now i also wanted to edit previous feature and also add attributes of feature that has been added (note: from this link edited feature cannot be edited again once web application is refresh and attributes cannot be added) . I also wanted to correct edited mistake without entirely removing the feature. I also refereed links like Ref Link 1 and Ref Link 2. Please help me out for same. From this i wanted to create an web base application for editing layers from geoserver. If any other opensource software if you can suggest for this web application please suggest. Please help me out for same.


title: Calling LineString.closestPoint or LineString.project methods from Python?

text: I'm trying to use the GEOS library from Python, specifically via Django (`django.contrib.gis.geos`). Looking at the GEOS documentation, there are two methods on `LineString` I want to call that aren't available from Python: `project` and `closestPoint`. Is there any reason why they're not available? Can they be added to the Python bindings? As far as I can tell, they're in `/usr/lib/libgeos_c.so.1.6.0`, is there a way to call them directly if they can't be added to the library?


title: How to Apply a Stretch to Multiple Rasters?

text: In ArcGIS 10.0, how can I apply a stretch, which I would normally do by selecting 2 standard deviations in the "Symbology" tab of Layer Properties, to 200 rasters simultaneously?


title: How to install add-ons in GRASS?

text: How do I install the r.inund.fluv script in GRASS? https://svn.osgeo.org/grass/grass-addons/grass6/raster/r.inund.fluv/


title: download ESRI World Imagery basemap for offline use?

text: So our search and rescue team would like to be able to view ESRI's World Imagery basemap through the free tool ( _ArcGIS Explorer_ ). Thing is, we will be offline when we are in the field. So - is there a way to download our region of the ESRI World Imagery basemap, for use offline? Note, I'm not a GIS professional, so am probably mixing some terminology. I found the URL where the World Imagery basemap data lives, and the query form to get back a flat image (jpg I think). But, we would really like to have it integrated as a basemap so we can use it seamlessly in `ArcGIS Explorer`. (Is a basemap a different beast than just a raster layer?) Is there any way to get this as a 'Data Appliance' from ESRI? link Our county staff uses the full paid ArcGIS tools so they might be able to get them for our SAR team's use - the issue is that we would want to view it offline in the free tool on our laptops in the field. We really do want to use ESRI World Imagery, or an equivalent (`it is 1 foot = 1 pixel in our area`). Thanks **UPDATE - 2013-10-01 - SOLVED:** gmapcatcher does the trick perfectly. It doesn't read from ESRI World Imagery, but it does download tiles from Bing and others (0.3m in our area) for offline viewing. It's not very feature-rich but that is a big part of the appeal! Learning curve is just about zero. There's probably a fairly easy way to turn those downloaded tiles into a correct directory structure for serving a tileset on MapServer, in order to use it with more feature-packed tools. If anyone finds a quick way to do that tile-file-rearranging please post here, thanks!


title: Two rasters created from the same polygon have different extents?

text: My problem is that i am trying to export some rasters to Rockworks so that i can create a sand layer from them. Rockworks requires the grid files i use to have the exact same XY co-ordinates. So i have a polygon of my sand layer - i create one raster of the surface (by using extract by mask from some AHN data i have) and the second one i create myself with a depth of -10. When i check their layer properties however they are different! e.g. the top extent of the lower raster is 503002.235094748 and the top extent of the upper raster is 503000.235094748. The left and bottom extents are exactly the same.. Theyre both in the same co-ordinate system..?! Thank you!


title: How to use the multiple attributes feature in ArcGIS 10

text: I am trying to make a map using US census data showing vacant housing to total housing on the same map. I came across the "multiple attributes" feature under the symbology tab and wondered if this is the way to go and if so how to use it? Ive only been able to get one type of symbol, the gradient feature to appear. Thanks for the help!


title: Building FGDB support for GDAL 1.10 on Ubuntu

text: Because data.gc.ca is now distributing some of its data in FGDB format (example: Plant Hardiness Zones), I would like to be able to use it with gdal/QGIS. The compile/install is not going as planned, however. These are the steps I've taken:   1. Enabled the source repository for ubuntugis-unstable.   2. Downloaded the Linux 64-bit File Geodatabase API 1.3, installed it in `/usr/local/FileGDB_API`, added the necessary ldconfig setting, and successfully built and run the test programs. (As suggested procedure here: [Qgis-user] How to compile gdal with ESRI FileGDB)   3. Installed the necessary source dependencies: `sudo apt-get build-dep gdal`   4. Pulled down the source: `apt-get source gdal`   5. Added the line `\--with-fgdb=/usr/local/FileGDB_API \` to the configure command in the `debian/rules` file.   6. From the source directory, issued `dpkg-buildpackage -rfakeroot -uc -b` (as advised here: How do I get the source code of packages installed through apt-get? - Ask Ubuntu) The build goes well, but fails late in the package compilation stage, complaining about a missing Ruby 1.9 tree in the fakeroot filesystem:               …     mv /home/scruss/gdal-1.10.0/debian/tmp/usr/local/lib/site_ruby/1.9/* /home/scruss/gdal-1.10.0/debian/tmp/usr/lib/ruby/1.9/.     mv: cannot stat ‘/home/scruss/gdal-1.10.0/debian/tmp/usr/local/lib/site_ruby/1.9/*’: No such file or directory     make: *** [install] Error 1     dpkg-buildpackage: error: fakeroot debian/rules binary gave error exit status 2      As this is unrelated to FGDB, it looks like I've got everything else in place. But it also likely means that gdal wouldn't build anyway. I'd rather stick with the ubuntugis packaging, as Ubuntu/Debian gets difficult with source packages built outside its control.


title: Openlayers selection

text: Is it possible to select a polygon from a predefined Lat & Lon programmatically? eg location point: x='1234' y='5678' intersects polygon '1234', therefore SELECT POLYGON '1234' So far I have managed to get the lat lon of the drop point:               $( "#cow" ).draggable({ revert: true});             $( "#map" ).droppable({             drop: function( event, ui ) {                   var point = map.getLonLatFromPixel( map.events.getMousePosition(event, ui   ))                  OpenLayers.Util.getElement("pos").innerHTML = point                       }             });      Now I just neeed to select the polygon that x and y intersect. Ideas?


title: Watershed other than stream ( between two points)

text: Is there any method available to demarcate watershed for a location other than stream? I mean, if I want to find out the watershed for a place, where no stream available but, sheet of water can flow (like a causeway / Irish bridge). Thanks....


title: Why does same image in two ArcGIS ArcMap configurations show different extent numbers?

text: I would like to get hear from you about the situation I am seeing here. What It got my eye was that when I download an ASTER from Reverb/ Echo from USGS website. At work, I load my ASTER file to my ArcMap and the ASTER shows the extent is in decimal degrees and the spatial reference shows the UTM Zone 13 Northern Hemisphere.So, I was not able to overlay with my known Projected Coordinate System in NAD 1983 UTM 13 Zone to overlay each other. I have been having some projection issues because I want it to reprojection to NAD 1983 UTM Zone 13N. Now I backup the file and took it home and load it to my home ArcMap and it shows the extent is in meters and the spatial reference shows UTM Zone 13 Northern Hemisphere. I am wondering if this has to do with ArcGIS 10.1 at work was not installed properly or did I did do something that cause the different coordinate system on my ArcGIS at work? I would like to find the solution because I know someone else had the same situation I am in. I revised the error from home to work


title: Select Single Feature

text: I'm making a plugin in quantum gis 1.8 Lisbon. And I intend to select the value of fields in a table, when I click the mouse over a point. (Ie something similar to Single Feature Select tool). Does anyone know how to do? I also need to do queries in the Postgres database table, but I can not. I can only get the name of the tables and fields, but not the values. Can someone help me? thank you


title: Using Python to create new point location based on various azimuth and distances

text: I have a feature class with thousands of points. In the same table with these point locations, I have two additional fields. One for azimuth, and one for a length (or distance) to move along that azimuth in order to create a new point location. I am looking for guidance on a python script that would accomplish this task? table example x y length_fornew_point azimuth 400460.99 135836.76 5150.04 55.41 .... .... ultimately the script would run through my table, look at the existing coordinates and create a new end point at corresponding specified distance/azimuth in the table. These two coordinates would then be used to create a new two point polyline. Thanks much for the help


title: Iterate features on SpatialPolygonsDataFrame in R

text: I have a `SpatialPolygonsDataFrame` object. Is there way to use `apply` on this object rather than using a for loop to iterate features in this object?


title: How to abort a layer event in Openlayers

text: I have a problem ( How to stop features to redraw unwillingly in OpenLayers ) that might be solved by aborting the beforefeatureremoved or beforefeaturesremoved event. What I would like to to is to stop the layer from refreshing (removing and the adding again) features from a vector layer. If I for example make a javscript error in those methods the features isn't removed or refreshed. But causing a javascripts error causes other issues so I hope that there is a controlled way of doing this in Openlayers I have tried the               layer.events.on({        'beforefeatureremoved':function(evt){            OpenLayers.Event.stop(evt);         }     });      But it doesnt seem to do anything, Im I Br Mike


title: Where can I find a shapefile of all the Montreal Neighborhoods?

text: I am trying to map out Montreal and was wondering if there is already an existing polygon shapefile of all the neighborhoods. Here is a link to a perfect example of what I am looking for. http://en.wikipedia.org/w/index.php?title=File:ArrondissementsQuartiersMontreal.svg&page=1 If I can get that in shapefile format, it would save me a lot of time! Thank you!


title: How to use getfeatureinfo with MapServer and OpenLayers?

text: I'm trying to use the `getfeatureinfo` function of Openlayers with `MapServer`. When I click on a point on my map, I get the following error in Firebug:               <?xml version='1.0' encoding="ISO-8859-1" standalone="no" ?>     <!DOCTYPE ServiceExceptionReport SYSTEM "http://schemas.opengis.net/wms/1.1.1/exception_1_1_1.dtd">     <ServiceExceptionReport version="1.1.1">     <ServiceException code="LayerNotQueryable">     msWMSFeatureInfo(): WMS server error. Requested layer(s) are not queryable.     </ServiceException>     </ServiceExceptionReport>      I've included the javascript and mapfile that I am using. As far as I can tell, I have the correct arguments in my Mapfile but it does not appear to enable querying using `getfeatureinfo`. I'm also a little unclear on how my html templates should be structured or if they are essential for this. This is the javascript (I am using OpenLayers 2.12):                   info = new OpenLayers.Control.WMSGetFeatureInfo({         url: 'http://mysite.com/mapserv.cgi?MAP=points.map',          title: 'Identify features by clicking',         queryVisible: true,         eventListeners: {             getfeatureinfo: function(event) {                console.log(event.text)                               }         }     });     map.addControl(info);     info.activate();      And this is my mapfile, `points.map`:               MAP       IMAGETYPE      PNG       EXTENT         -9780000 2770000 -8850000 3650000         SIZE           500 600       SHAPEPATH      "/home/data"       IMAGECOLOR     255 255 255       PROJECTION         "init=epsg:900913"       END       CONFIG          "MS_ERRORFILE" "/home/tmp/ms_error.txt"       DEBUG 5                WEB         METADATA             "wms_title"           "WMS Demo Server"             "wms_onlineresource"  "http://mysite.com/mapserv.cgi?MAP=points.map&"             "wms_srs"             "EPSG:900913"             "wms_enable_request"  "*"             "ows_enable_request"  "*"         END       END            ##### symbol definition       SYMBOL         NAME "circle"         TYPE ellipse         FILLED true         POINTS             1 1         END       END            ####### layer       LAYER          NAME            points         CONNECTIONTYPE  POSTGIS         CONNECTION      "host=localhost dbname=data user=postgres port=5432"         DATA            "geom from points using unique gid using SRID=3086"         HEADER          'head.html'           FOOTER          'foot.html'         METADATA                     "wms_enable_request" "*"                 "gml_include_items" "all"             "wms_include_items" "all"                    END           STATUS          ON         TYPE            POINT         PROJECTION             "init=epsg:3086"         END         CLASS           NAME 'point symbols'           STYLE             SYMBOL "circle"             COLOR         100 100 100             OUTLINECOLOR  10  10  10             SIZE 50           END         END       END         END       My template files are very basic - just `<html><body>` for the `HEADER` and `</body></html>` for the `FOOTER`. Where am I going wrong in this overall process?


title: Routing with multiple points

text: I am using postgreSQL, postGIS and pgrouting. Now I am really confused on how am I going to find the shortest path between multiple points. I know how to find it between two points through dijkstra but I am confused how to accomplish the same result with more than two points. I have done some research and found this: https://github.com/pgRouting/pgrouting/wiki/One_to_many-Dijkstra---To-review But still I did not understand if its the right way and how to install it.


title: QGIS - Add a combobox to a toolbar

text: Is there a way to add a combobox to a toolbar for QGIS? So far I can only find examples that add buttons. I found this tutorial which seems to explain how to do this with general PyQT but haven't had any luck at all doing this with QGIS. So any help or guidance would be greatly appreciated.


title: How to insert a "basemap" into ArcScene?

text: I want to model a 3D scene, showing layers floating above a flat map. Currently, I am able to show the floating layers, but I cannot add the bottommost map layer. I just want to add a "basemap", like ArcMap allows, but none of the files seem to be compatible. Just saving a basemap as a lyr file in ArcMap and importing it in ArcScene does not work, and raster images obtained from the Internet do not line up properly when added to the scene. I know that the online basemaps that are available in ArcMap do not work in ArcScene. I was wondering if there was any way to add a simple, flat map to a scene in ArcScene? Otherwise, just a simple US map DEM or something dimensional like that would work. All I really need is a reference for the floating layers. Is there somewhere I can download ArcScene maps? The illustration below illustrates what I want as my final product: ![The intended final product](http://i.stack.imgur.com/A3Jm5.png)


title: Obtaining bearings to polygon extents in ArcGIS?

text: In the context of a wider question, (Generating multiple visibility bearings in ArcGIS as automated procedure?) is it possible to extract bearing information from a particular point or set of points, to extreme tangents of a polygon? For example, I have a point shapefile with a set of viewpoints (small red dots in image below). I also have a **polygon** shapefile containing the _buffer_ of another set of points, which outline the exact extents, in m, of my development points (shown in blue below.) I want to know the extents, in bearing format, of the furthest left tangent point to the furthest right tangent point of the buffer shapefile of my development points. (the red arrows below, i.e. 275 - 282 degrees.) ![Example: the red dashed lines show extents of development visible from each VP.](http://i.stack.imgur.com/ggzuF.jpg). Can I do this in GIS, and if so how / what tools should I be looking at?


title: How to Migrate web maps from ArcGIS Flex to JavaScript API?

text: A client has several web maps that were made using ArcGIS Flex. But the client needs all the current web maps to be changed to ArcGIS JavaScript API. What steps need to be taken to make this happen efficiently?


title: Creation of tiled grid in JAVA with geotols

text: I am trying to create a grid using GeoTools in JAVA, covering specyfic area with resolution 50 m (distance between two points). To do so I have two points: the most nortwestern (let's call it NW) and the most southeastern one. I created a table and put the NW point in tab[0][0]. Starting from tab[0][0] using the Geodetic Calculator (link to the doc) I am calculating the next points and putting them in right place into my table. Here is the code I developed:               public void createGrid(double granularity) {         logger.info(GridManager.class + " creation of grid called");              int defaultCellID = 0;         float defaultSigLevel = -1000;              /*  __ __ __ __ __           * |NW|__|__|__|NE|           * |__|__|__|__|__|          * |__|__|__|__|__|          * |SW|__|__|__|SE|          */              // mind, this is lon/lat         Point2D.Double nePoint = new Point2D.Double();         nePoint.setLocation(sePoint.getX(), nwPoint.getY());              Point2D.Double swPoint = new Point2D.Double();         swPoint.setLocation(nwPoint.getX(), sePoint.getY());              // calculations of a width of the grid         this.gc.setStartingGeographicPoint(nwPoint);         this.gc.setDestinationGeographicPoint(nePoint);         double d = this.gc.getOrthodromicDistance();         int width = (int) Math.ceil(d / granularity);              // calculations of a height of the grid         this.gc.setDestinationGeographicPoint(swPoint);         d = this.gc.getOrthodromicDistance();         int height = (int) Math.ceil(d / granularity);              // grid creation         grid = new TileCartesian[height][width];         grid[0][0] = new TileCartesian((float) nwPoint.getY(), (float) nwPoint.getX(), defaultSigLevel, defaultCellID, 0, 0);         double dist;              for (int i = 0; i < height; i++) {             this.gc.setStartingGeographicPoint(grid[i][0].getLon(),grid[i][0].getLat());             dist = granularity;                  for (int j = 1; j < width; j++) {                      this.gc.setDirection(90 /* azimuth */, dist);                      Point2D.Double p = (Point2D.Double) this.gc                         .getDestinationGeographicPoint();                      grid[i][j] = new TileCartesian((float) p.getY(),                         (float) p.getX(), defaultSigLevel, defaultCellID, i, j);                 dist += granularity;             }                  if (i != height-1) {                 this.gc.setDirection(180 /* azimuth */, granularity);                 Point2D.Double p = (Point2D.Double) this.gc                         .getDestinationGeographicPoint();                 grid[i+1][0] = new TileCartesian((float) p.getY(), (float) p.getX(), defaultSigLevel, defaultCellID, i+1, 0);             }         }                   logger.info(GridManager.class + " Grid created");     }      My problem is that the most SE point of the calculated grid is far from the assumed one. And when I am exporting the result to GeoTIFF this cause deformation. I really do not know how to solve this problem. Maybe there exists any other way to create such a grid, but I couldn't find any solution which would fit my needs.


title: tools and ETLs for migrating data benchmarking - from shapefile to SQL SERVER

text: i need to migrate my shapefiles to my database SQL SERVER, and i don't know what's the best tool or ETLs to choose, i already tried FME Desktop and i succeded to load the shapefiles on my db, but while searching i found other tools (like : artend, data migration tool kit, ETL integrator, Scriptella ...). Now i am confused, i don't know what to choose amoung many open source and licenced solution. Plz help.


title: Can I create a TIN from contours lines (shp) in QGIS?

text: I need to do a slope analysis, and would like to make a TIN from 2m contours lines. Or do a slope analysis directly from contours? The only DEM available is not precise enough. I do not have any elevation raster files. And how do I then make slope classes? The TIN properties does not seem to be able to classify anything other than elevation.


title: Why can't I import arcpy?

text: >>> import arcpy     Traceback (most recent call last):       File "<interactive input>", line 1, in <module>       File "C:\Python27\lib\arcpy\__init__.py", line 21, in <module>         from arcpy.geoprocessing import gp       File "C:\Python27\lib\arcpy\geoprocessing\__init__.py", line 14, in <module>         from _base import *       File "C:\Python27\lib\arcpy\geoprocessing\_base.py", line 14, in <module>         import arcgisscripting     ImportError: DLL load failed: The specified module could not be found.      This is the error I get when i try to import arcpy. Any suggestions?


title: How is IValidate.GetInvalidRules supposed to work?

text: we are currently trying to programmatically validate features against geodatabase-defined rules. This includes attribute rules as well as connectivity and relationship rules. We have consulted the documentation for IValidate, IValidation and the Validating features article on EDN to learn how it's done. After setting up the appropriate connectivity rules, invocation of the Validate Features in ArcMap correctly yields the number of invalid features in a selection. If I select a single feature, invocation of this command correctly displays what is wrong about that particular feature in a message box. Now, calling IValidate.Validate on that feature works as well. Its return value indicates that the feature is indeed invalid, and the error message is identical to that displayed by the Validate Features command in ArcMap. However, calling this method has two disadvantages which are problematic for our purposes:   1. It performs the validation in a short-circuit manner. As soon as one of the assigned rules is invalid, it stops evaluating the remaining rules. We need to evaluate all the rules so that we can display the violations at once.   2. It provides the error information as a string. For our purposes, we would ideally like to examine the information in a more structured way, which means we need to access the IRules that turned out to be invalid. So, this is where IValidate.GetInvalidRules comes into play. The documentation for this method states: > The GetInvalidRules method returns an enumerator of rules that are currently > violated by the row. which is exactly what we are after. However, the enumeration returned by this method does not seem to ever contain anything. What is particularly weird is that:   1. Calling `IValidate.Validate` on a feature yields the appropriate validation error message as expected.   2. Calling `IValidate.GetInvalidRules` on the very same feature always returns an empty enumeration as if no connectivity or relationship rules were invalid. (It does not matter whether the call is made within an edit session or not.) How can I validate a row against all the assigned geodatabase rules at once? Has anyone used `IValidate.GetInvalidRules` successfully? **UPDATE:** `IValidate.GetInvalidRules` seems to return invalid rules which pertain to attributes. However, invalid network connectivity rules are still never returned from this method.


title: Synchronize data edits between ArcMap and ArcGIS Online

text: What, if any, workflows are available to round trip data between ArcGIS Desktop and an _ArcGIS Online for Organizations_ (AGOL) web map with editable feature layers? For example, we have a web map with an editable feature service layer. Some staff might add features to the map using ArcGIS.com Map Viewer, others with ArcGIS Online Explorer, and still others with ArcGIS Desktop. Nothing needs to be done for the first two, as soon as the map is saved the other application will see it. However for ArcGIS Desktop, the so-called powerhouse in this equation, there doesn't seem to be a way to synchronize or upload changes back to ArcGIS Online. Is this really the case? There aren't any SDE or databases in the mix, just ArcGIS Desktop with file geodatabase and AGOL. **Update:** ArcGIS Server is not involved at all (disregarding that's what Esri uses power AGOL behind the scenes, as there is no meaningful interaction with it from a user perspective).


title: Regenerate Geowebcache Tiles on PostGIS change

text: I am using Geoserver with Geocache to get WMS tiles and a PostGIS backend. There are 2 systems using the PostGIS database: one that is with Geoserver and one that is an old system that just uses the database. All changes made to the geometry using either system need to cause the tiles to be recreated, and we can't change the code in the old system. I think I need to use a trigger in PostGIS. I understand how to create the trigger, but how do I make it regenerate the tiles that include the point that was changed?


Following the pgRouting workshop here, I have the routing query below:               $sql = "SELECT gid, start_id AS source, end_id AS target, the_geom,                        distance(the_geom, GeometryFromText(                            'POINT(".$lonlat[0]." ".$lonlat[1].")', 4326)) AS dist                       FROM ".TABLE."                        WHERE the_geom && setsrid(                            'BOX3D(".($lonlat[0]-0.1)."                                    ".($lonlat[1]-0.1).",                                    ".($lonlat[0]+0.1)."                                    ".($lonlat[1]+0.1).")'::box3d, 4326)                       ORDER BY dist LIMIT 1";      This query fails with the error: `Query failed: ERROR: Operation on two geometries with different SRIDs`. Now, I've been looking around and came across this. I've tried the suggestions on this post, but the error persists. Any idea on how to fix this will be appreciated. **UPDATE:** Running `SELECT st_srid(the_geom) FROM table_name` against my database table returns `-1`. In my query and map, the assumed external projection is `EPSG:4326`. Is changing this value likely to help? Is there any harm in changing the srid value directly in the database? How can I go about changing the srid value in the database if it doesn't present any risk? **UPDATE:** I had to change the srid of the column from the `UNKNOWN` indicator of `-1` to `4326`. The answer at this StackOverflow post was very helpful in this. Also, the question in this SE GIS post was instrumental in helping me in the right direction.


title: How to use WHERE CLAUSE to filter the features by its attribute in PyQgis?

text: How to Get features by their attributes (similar to Iqueryfilter in arcobjects) in Qgis using python? Instead of getting all features and filter it manually, is there any option to use whereclause to filter it out? For Example: I have a field name called 'Counties'. It has more than fifty thousand of features.i.e not possible to fetch all features and filter it because of time consuming. So i can query it by using iqueryfilter.whereclause = 'Counties = Norwich' in arcobjects. Similar thing i need in PyQgis. kindly let me know whether its doable. If so, please give me an example


title: Two different ways of labeling in QGIS Layouter

text: I've made a map with points and labels related to those points. In one area, the points are very dense, therefore, I'd prefer to create one small scale map and in the corner a large scale map, to get more detail. Now, my problem is that - for the dense part - labels should only be put on the large scale map, for the other points, the label should be present on the small scale map. How would you solve this? I can choose labels for one map (like in: How to only show labels for an arbitrary selection of items?), make a layouter for the map, and subsequently do the same for the other map. But than, if you want to change something to the first, labels are updated again to the actual state of your qgis project.


title: display multiple vector layers from single PostGis table in GeoServer using OpenLayers -- REWORDED

text: Prompting suggestions I’m going to attempt to reword my question I have a PostGis table that contains overlapping vector polygons. I need to display a selection in Geoserver and provide functionality for the user to toggle on & off unwanted polygons. I need to add multiple layers from a PostGis Table. I’ve made a start with the code below but I’m stuck with the following:   * Why are all the layers Base Layers? Land should be the only base layer the rest overlays?   * Can all overlay layers be on at once?   * How do I move the layerswitcher control outside the map area? When I click inside the map how do I show the attribution?   * How do I create a new layer by script either with a polygon given a   co-ordinate string, or a point and radius?   * How do I create a list of all the selected layers? ` body { margin: 1em; } #map { width: 800px; height: 475px; border: 1px solid black; } OpenLayers.IMAGE_RELOAD_ATTEMPTS = 3;                  var map;         function init(){             map = new OpenLayers.Map('map');             var Base = new OpenLayers.Layer.WMS(                 "Land",                 "http://localhost:8080/geoserver/cite/wms",                 {LAYERS: 'cite:land',                  format: 'image/png',                  isBaseLayer: true                      }             );                  var L10299 = new OpenLayers.Layer.WMS(                 "10299",                 "http://localhost:8080/geoserver/cite/wms",                 {                 VIEWPARAMS: 'kathy:10299',                 LAYERS: 'cite:Q_Results',                 isBaseLayer: false                 }             );                       var L43352 = new OpenLayers.Layer.WMS(                 "43352",                 "http://localhost:8080/geoserver/cite/wms",                 {                 VIEWPARAMS: 'kathy:43352',                 LAYERS: 'cite:Q_Results',                 isBaseLayer: false                 }             );                       var L25177 = new OpenLayers.Layer.WMS(                 "25177",                 "http://localhost:8080/geoserver/cite/wms",                 {                 VIEWPARAMS: 'kathy:25177',                 LAYERS: 'cite:Q_Results',                 isBaseLayer: false                 }             );                       map.addLayers([Base, L10299, L43352, L25177]);                            map.addControl(new OpenLayers.Control.LayerSwitcher());                  map.zoomToExtent(                 new OpenLayers.Bounds(12.996, 52.172,13.931, 52.958)             );         }     </script>      ### OpenLayers PostGis Layers `


title: How to update feature class by iterating through values with branching?

text: I'm really struggling with getting this to work. Basically, I want my script to get the "newID" value and use that as a look up value. If the "newID" is found in the OldID list then do one thing, else do something else. The problem I'm having is it only grabs the first ID value and then does loops the rest of the code based on whatever category the first ID fell under. I need it to get the value of each field in the new data table, and write it to the corresponding field in the old data table. It does this just fine, but obviously it isn't checking the ID's after the first one, or updating/appending the geometry. Here's what I've got so far. Just showing applicable code....               rows = gp.SearchCursor(new_data)     row = rows.next()          while row:         newID = row.getValue('ID')         newShape = row.SHAPE         y = 0              while newID in OldID_List:             gp.AddMessage('')             print 'Starting the updating for ratios'             where = 'ID' + ' = ' + str(newID)             new_rows = gp.UpdateCursor(old_data, where)             new_row = new_rows.next()                  while y <= len(Field_List):  # checks to see if ID exists -                                           # this means just update                 for each in Field_List:                     newValue = row.getValue(each)                          if newValue is None:                         new_row.setValue(str(each), 0)                     else:                         new_row.setValue(str(each), str(newValue))                          y += 1                          if y > len(Field_List):                         new_row.Shape = newShape                         new_rows.UpdateRow(new_row)                         new_row = new_rows.next()                         row = rows.next()                         newID = row.getValue('ID')                         newShape = row.SHAPE                         y = 0              while newID not in OldID_List:  # these are new objects to be appended             flag = 1             nsert_Rows = gp.InsertCursor(old_data)             nsert_Row = nsert_Rows.NewRow()                  while y <= len(Field_List):                 for each in Field_List:                     newValue = row.getValue(each)                     print 'Starting appends'                          if newValue is None:                         nsert_Row.setValue(str(each), 0)                     else:                         nsert_Row.setValue(str(each), str(newValue))                          y += 1                     if y > len(Field_List):                         nsert_Row.ID = newID                         nsert_Row.Shape = newShape                         nsert_Rows.insertrow(nsert_Row)                              row = rows.next()                         newID = row.getValue('ID')                         newShape = row.SHAPE                              appends += 1                         y = 0     del rows     del new_rows     if flag == 1:         del nsert_Rows      Any help or suggestions would be appreciated. Thanks.


title: Optimizing osm2pgsql imports for OSM data

text: I'm currently building an instance on EC2 on which to import the entire Planet.osm snapshot of the whole Earth's worth of data for some projects we're working on. I've spun up a large Ubuntu x64 instance and attached plenty of separate storage on an EBS volume for the Postgres database and modified it to house the PGSQL data there. Now the server is having trouble using `osm2pgsql` to import the snapshot... After a couple of attempts with different memory configs and whatnot, the process keeps outputting "Killed" after getting most of the way through; once it was killed while "going over pending ways" and the next time, after slightly adjusting the slim cache, it reached "processing ways" before crashing out. From what I've read, this is generally due to memory issues. Here's my latest attempt to run the import:               osm2pgsql -v -U osm -s -C 4096 -S default.style -d osm /data/osm/planet-latest.osm.bz2      And here are the specs for a Large instance on EC2: > Large Instance 7.5 GB of memory, 4 EC2 Compute Units (2 virtual cores with 2 > EC2 Compute Units each), 850 GB of local instance storage, 64-bit platform My question is -- are there some good benchmark resources to determine the tuning requirements for osm2pgsql and Postgres? Speed of import isn't even that important to me, I'd just like to be able to make sure the process completes safely, even if it takes 4 or 5 days... I've read through Frederick Ramm's "Optimising the rendering chain" (PDF) document from last year's SOTM, but are there other good opinions / resources?


title: ogr2ogr selecting features by attributes

text: I am writing a shapefile after some processing; however, in the final step I need to do some selection by attributes from the the shapefile. I am using the command in a shell with the aim of using it in a Python script when it works.               ogr2ogr -f "ESRI Shapefile" -select * where ID="1" outfile.shp infile.shp      I am getting the error message:               FAILURE:      Unable to open datasource `Downloads' with the following drivers.      What could I be doing wrong? Your assistance is much appreciated.


title: Shapefile to KML automatic updates (ModelBuilder)

text: I am trying to set a model to where it can automatically update a KML only when there are new features added or modified to my shapefile. Trying to avoid having to convert periodically. Can anyone help?


title: Following path in OpenStreetMap

text: I want to visualize a road present in the OpenStreetMap. Using OpenLayers, I try to do in the following way:   1. Take a set of co-ordinates which should represent the two sides of the road.   2. Join both the points on the map. (leads to a straight line) Now, the problem with this approach is that it is just a straight line (obviously), what is the way of following the path along the road ? Just a note that I'm not tied up with OpenLayers or anything. I can use any mapping library as long as it achieves this functionality.


title: How to change network analysis settings in a service from a web app

text: I'm a total n00b with ArcGIS Server, so I'd like to know if I publish a network analysis service using ArcGIS for Server, will the user be able to change analysis settings as one would in ArcMap e.g. changing impedance from distance to travel time, adding time options etc? I've looked at the server help for NA as well as the network analysis service help, but it doesn't really say. After some looking around, it seems that creating a model with the Network Analyst tools I need (such as `Make Route Layer`), making the attributes and other relevant components model parameters, and publishing it as a geoprocessing service would be the way to go. Again, I am just getting started with this topic, so these are just my initial observations.


title: How to highlight kml features in openlayers

text: I'm trying to make the features from a kml layer highlight on hover, but it doesn't seem to work for me. On the other hand, if i create a vector layer with a polygon feature it works just fine. I made an example with kml and polygon feature just to try to figure out what the difference is. Does anyone know what I am missing? This is the .js code:               var map;          function init(){             map = new OpenLayers.Map('map');                  var kml = new OpenLayers.Layer.Vector ('Australia Boundaries',{isBaseLayer:true, protocol: new OpenLayers.Protocol.HTTP ({url: 'kml/australia.kml', format: new OpenLayers.Format.KML({ extractStyles: true, extractAttributes: true})}),projection: new OpenLayers.Projection("EPSG:4326"), strategies: [new OpenLayers.Strategy.Fixed()]});             map.addLayers([kml]);              //create vector layer with a polygon feature                  var vector = new OpenLayers.Layer.Vector('Vector', {projection: new OpenLayers.Projection("EPSG:4326") });                  var point = new OpenLayers.Geometry.Point(90.752725,10.475584);             var pointList = [];             for(var p=0; p<6; ++p) {                 var a = p * (2 * Math.PI) / 7;                 var r = Math.random(1) + 20;                 var newPoint = new OpenLayers.Geometry.Point(point.x + (r * Math.cos(a)),                                                              point.y + (r * Math.sin(a)));                 pointList.push(newPoint);             }             pointList.push(pointList[0]);                  var linearRing = new OpenLayers.Geometry.LinearRing(pointList);             var polygonFeature = new OpenLayers.Feature.Vector(                 new OpenLayers.Geometry.Polygon([linearRing]));             vector.addFeatures([polygonFeature]);             map.addLayer(vector);             //create highlight and select controls for the vector and kml layers                  var highlightCtrlVector = new OpenLayers.Control.SelectFeature(vector, {                 hover: true,                 highlightOnly: true,                  renderIntent: 'temporary'});                  var selectCtrlVector = new OpenLayers.Control.SelectFeature(vector,                 {clickout: true}             );             var highlightCtrlKml = new OpenLayers.Control.SelectFeature(kml, {                 hover: true,                 highlightOnly: true,                 renderIntent:'temporary'});                  var selectCtrlKml = new OpenLayers.Control.SelectFeature(kml,                 {clickout: true}             );             map.addControl(highlightCtrlVector);             map.addControl(selectCtrlVector);             map.addControl(highlightCtrlKml);             map.addControl(selectCtrlKml);             highlightCtrlVector.activate();             highlightCtrlKml.activate();             selectCtrlVector.activate();             selectCtrlKml.activate();                  map.addControl(new OpenLayers.Control.LayerSwitcher());             map.zoomToMaxExtent();         }      And here is the example: http://geoteste.comli.com/app/indexTest.html


title: Arc GIS 10.1 - How to sort Layers ascending at TOC in "List by drawing order"

text: I have some 100 features in my ArcMap Project. I choose at TOC the View mode "List by drawing order". How can I sort this list ascending or by name alphbetical? Under the options field, is nothing like this. The drawing order should be ascending by name. Or how can I sort them on other way? Thanks for any help. Luis


title: how to determine drive times like those available in google maps

text: I'm very new to GIS but have played around a bit with postGIS and qGIS. I'm wondering if there's a way to generate drive time estimates like those given by google maps using open source tools. For a small number of routes I can query google maps, and scrape results but I'm sure there's a better way. Thanks


title: How to join 2 shape files and create a new shape file using ArcObjects?

text: Newbie in arcobject programming. Using C#. Developing a desktop application. I have 2 shape files. One is a point shape file Called **Addresses** and another Polygon shape file which is called **Area**. The Address shape file has an attribute called **Address_Number** while the Area shape file has a attribute called **Area_Number**. I need to select only the features that match "Address_Number" (from Addresses) to "Area_Number" (from Area) and then create a new shape file. How can I do this? Any help with some code example is greatly appreciated.


title: Google maps marker is not getting displayed

text: I am trying to learn google map API. I see a map but I am unable to get a marker. Here is my code **Activity**               public class HomeActivity extends Activity {              @Override         protected void onCreate(Bundle savedInstanceState) {             super.onCreate(savedInstanceState);             setContentView(R.layout.home_page);             GoogleMap map = ((MapFragment) getFragmentManager().findFragmentById(R.id.map)).getMap();             map.addMarker(new MarkerOptions()                     .position(new LatLng(37.7750, 122.4183))                     .title("San Francisco"));         }          }      **home_page.xml**               <?xml version="1.0" encoding="utf-8"?>     <fragment xmlns:android="http://schemas.android.com/apk/res/android"               android:id="@+id/map"               android:layout_width="match_parent"               android:layout_height="match_parent"               android:name="com.google.android.gms.maps.MapFragment"/>      What am I missing out?


title: Python, comtypes and ArcObjects: Error creating AppROT object

text: I am experimenting with `comtypes` and ArcObjects under Python 2.6.5 and ArcGIS 10 SP1. I'm using the pure Python method to wrap the ArcObjects OLBs described in this answer, but getting an error in the `comtypes.CoCreateInstance` method. Here is the code I am running:               def WrapModules():         #force wrapping of all ArcObjects libraries (OLBs)         import os         import comtypes.client         # change com_dir to whatever it is for you         com_dir = r'C:\Program Files\ArcGIS\Desktop10.0\com'         coms = [os.path.join(com_dir, x) for x in os.listdir(com_dir) if os.path.splitext(x)[1].upper() == '.OLB']         map(comtypes.client.GetModule, coms)          def GetApp():         """Get a hook into the current session of ArcMap"""         from comtypes.gen import esriFramework         pAppROT = NewObj(esriFramework.AppROT, esriFramework.IAppROT)         if pAppROT is not None:             iCount = pAppROT.Count             if iCount == 0:                 print 'No ArcGIS application currently running.  Terminating ...'                 return None             for i in range(iCount):                 pApp = pAppROT.Item(i)  #returns IApplication on AppRef                 if pApp.Name == 'ArcMap':                     print "ArcMap found"                     return pApp             print 'No ArcMap session is running at this time.'         print "No AppROT found"         return None          def NewObj(MyClass, MyInterface):         """Creates a new comtypes POINTER object where\n\         MyClass is the class to be instantiated,\n\         MyInterface is the interface to be assigned"""         from comtypes.client import CreateObject         import traceback         try:             ptr = CreateObject(MyClass, interface=MyInterface)             return ptr         except:             print traceback.format_exc()             return None          if __name__ == "__main__":         WrapModules()         pApp = GetApp()         if pApp is not None:             print "HWND: %d" % pApp.hWnd         else:             print "No ArcGIS application found!"      And here is the output from the script:               Traceback (most recent call last):       File "C:\temp\ComHelpers.py", line 35, in NewObj         ptr = CreateObject(MyClass, interface=MyInterface)       File "C:\Python26\ArcGIS10.0\lib\site-packages\comtypes\client\__init__.py", line 235, in CreateObject         obj = comtypes.CoCreateInstance(clsid, clsctx=clsctx, interface=interface)       File "C:\Python26\ArcGIS10.0\lib\site-packages\comtypes\__init__.py", line 1145, in CoCreateInstance         _ole32.CoCreateInstance(byref(clsid), punkouter, clsctx, byref(iid), byref(p))       File "_ctypes/callproc.c", line 925, in GetResult     WindowsError: [Error -2147221231] ClassFactory cannot supply requested class          No AppROT found     No ArcGIS application found! Thanks for any insights you might have!


title: Interpolation Technique

text: We have following solar metro logy parameter through satellite data with resolution of 100 Km^2 (10Kmx10Km)               Insolation on horizontal surface                 Diffuse radiation on horizontal surface                  Direct normal radiation             Like this we can have about 4-5 more parameters which are associated with Solar Power generation. These same parameters are also measured at 23 ground stations. We also have same parameter data which are taken at ground station at respective long-latitude locations. **Our objective** will be to estimate all those solar parameters at every 1 Km x 1Km area By interpolating satellite data and ground station data in given vicinity area with best possible accuracy My question is Which is best interpolation technique we can used please provide some references and the software which can be used ..it will be great help to me.


title: How to properly convert EPSG:3857 to EPSG:4326 tif?

text: I tried using the answer from "Transforming EPSG:3857 to EPSG:4326" to do this but gdalinfo tells me that this didn't really work properly.               Driver: GTiff/GeoTIFF     Files: geo-map.tif     Size is 1024, 512     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (-179.999999999999545,85.051100000000005)     Pixel Size = (0.351538008801935,-0.332291793794902)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  (-180.0000000,  85.0511000) (180d 0' 0.00"W, 85d 3' 3.96"N)     Lower Left  (-180.0000000, -85.0822984) (180d 0' 0.00"W, 85d 4'56.27"S)     Upper Right ( 179.9749210,  85.0511000) (179d58'29.72"E, 85d 3' 3.96"N)     Lower Right ( 179.9749210, -85.0822984) (179d58'29.72"E, 85d 4'56.27"S)     Center      (  -0.0125395,  -0.0155992) (  0d 0'45.14"W,  0d 0'56.16"S)     Band 1 Block=1024x2 Type=Byte, ColorInterp=Red       Mask Flags: PER_DATASET ALPHA      Band 2 Block=1024x2 Type=Byte, ColorInterp=Green       Mask Flags: PER_DATASET ALPHA      Band 3 Block=1024x2 Type=Byte, ColorInterp=Blue       Mask Flags: PER_DATASET ALPHA      Band 4 Block=1024x2 Type=Byte, ColorInterp=Alpha      The corner coordinates still projects to 180, 85.0511 which is for EPSG:3857. For EPSG:4326, isn't it supposed to be 180, 90? so I tried translating it again using:               gdal_translate -a_srs EPSG:4326 -a_ullr -180 90 180 -90 geo-map.tif new-geo-map.tif      and got this:               Driver: GTiff/GeoTIFF     Files: new-geo-map.tif     Size is 1024, 512     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (-180.000000000000000,90.000000000000000)     Pixel Size = (0.351562500000000,-0.351562500000000)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  (-180.0000000,  90.0000000) (180d 0' 0.00"W, 90d 0' 0.00"N)     Lower Left  (-180.0000000, -90.0000000) (180d 0' 0.00"W, 90d 0' 0.00"S)     Upper Right ( 180.0000000,  90.0000000) (180d 0' 0.00"E, 90d 0' 0.00"N)     Lower Right ( 180.0000000, -90.0000000) (180d 0' 0.00"E, 90d 0' 0.00"S)     Center      (   0.0000000,   0.0000000) (  0d 0' 0.01"E,  0d 0' 0.01"N)     Band 1 Block=1024x2 Type=Byte, ColorInterp=Red       Mask Flags: PER_DATASET ALPHA      Band 2 Block=1024x2 Type=Byte, ColorInterp=Green       Mask Flags: PER_DATASET ALPHA      Band 3 Block=1024x2 Type=Byte, ColorInterp=Blue       Mask Flags: PER_DATASET ALPHA      Band 4 Block=1024x2 Type=Byte, ColorInterp=Alpha      I realized that it's still not the right thing. The thing is, I need the .tiff to use in WhirlyGlobe but I always end up with a discrepancy. I tried Natural Earth's .tiff on the globe and it worked fine:               Driver: GTiff/GeoTIFF     Files: NE2_LR_LC_SR_W.tif            NE2_LR_LC_SR_W.tfw     Size is 16200, 8100     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (-180.000000000000000,90.000000000000000)     Pixel Size = (0.022222222222220,-0.022222222222220)     Metadata:       AREA_OR_POINT=Area       TIFFTAG_DATETIME=2012:07:16 09:57:53       TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)       TIFFTAG_SOFTWARE=Adobe Photoshop CS5 Macintosh       TIFFTAG_XRESOLUTION=72       TIFFTAG_YRESOLUTION=72     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  (-180.0000000,  90.0000000) (180d 0' 0.00"W, 90d 0' 0.00"N)     Lower Left  (-180.0000000, -90.0000000) (180d 0' 0.00"W, 90d 0' 0.00"S)     Upper Right ( 180.0000000,  90.0000000) (180d 0' 0.00"E, 90d 0' 0.00"N)     Lower Right ( 180.0000000, -90.0000000) (180d 0' 0.00"E, 90d 0' 0.00"S)     Center      (  -0.0000000,   0.0000000) (  0d 0' 0.00"W,  0d 0' 0.00"N)     Band 1 Block=16200x1 Type=Byte, ColorInterp=Red     Band 2 Block=16200x1 Type=Byte, ColorInterp=Green     Band 3 Block=16200x1 Type=Byte, ColorInterp=Blue      please tell me what I'm doing wrong? :(


title: Postgis csv and error loading into mysql

text: I want to load spatial data into MySQL using a CSV I have generated from a Postgres table and when loading it into MySQL database using phpmyadmin it complains "Cannot get geometry object from data you send to the GEOMETRY field." How best can I export PostGIS geometry into MySQL using a csv without generating errors? I have created a MySQL table as               CREATE TABLE cadastre (       gid int(11) auto_increment NOT NULL,         geom MultiPolygon NOT NULL,         lpkey varchar(50),          SPATIAL INDEX(geom),          PRIMARY KEY (`gid`)        ) ENGINE=MyISAM  DEFAULT CHARSET=latin1;


title: Geoexplorer not displaying geoserver layers

text: I have a problem with geoexplorer not displaying my geoserver layers anymore. It gives me an error of "Non-existing source 'local' referenced in layer config." I really dont know hwta went wrong or what to do.


title: How to use gdal module in python

text: I am trying to import gdal modules from osgeo installed in my win64x machine but I get the following error: > Traceback (most recent call last):   >  File "", line 1, in   >  from osgeo import ogr   >  ImportError: No module named osgeo I am using the python shell from arcgis 10.1 that comes with python 2.7.2


title: Arcpy - deleted my entire workspace?

text: I've got a script which runs through some GPS data and creates several outputs (tables, shapefiles etc.). However, when I've just run it, the folder where all the input data lies (along with most of the outputs) has just vanished without a trace. It's not even the recycling bin. The only references to 'Delete' are deleting fields, deleting a table if it exists (*arcpy.Delete_management(gpsStreetSummary)*) and at the start deleting a GDB if it exists. There's nothing I can find in the script to explain why it would just delete everything?? Does anyone know if there is a bug in arcpy that sometimes causes this to happen?


title: DynamicCacheLayerManager.Connect deletes my tiles

text: You are supposed to be able to connect to your existing tiles with this command.               dynamicCacheLayerManager.Connect(tileCacheDir, tileCachePath);      When I do that my tiles are destroyed and rebuilt. Clearly this isn't the desired behavior, does anyone have any tips? I'm using ArcEngine 10 in C# with VS2010.


title: Qgis 1.8.0 throws error "Can not find the entry point sqlite3_rtree_geometry_callback"

text: Loading QGIS on Window 8 causes this problem: `Can not find the entry point sqlite3_rtree_geometry_callback in the procedure of dynamic loading C:\\....\Quantum GIS Lisboa \ bin \ spatialite.dll` Everything worked properly for many months. Thanks.


title: Data organization strategies for .csv and .shp

text: I'm organizing a large amount of data for a project and the data is stored in .csv format, with a code to relate each row to a polygon in a .shp. I was planning on converting the .csv data into .dbf and importing it into a geodatabase, all using ESRI tools. Is this the best way of organizing such data? (The alternative I thought of was to join the .shp and the .csv and to export the data so I have a permanent link, but it seems like I would have a lot of duplicated data) I'm interested in learning how best to improve my work flow and organize this data. I am planning on doing it all using ArcGIS 10 and using python where possible.


title: How do I set up the ArcGIS Viewer for Flex?

text: I am still a newbie at the Flex Viewer; I have been using the Out of the Box Web Application for Server 10.0. I tried to put a simple service as an additional basemap (I also tried it as an additional operational layer) and now the viewer will not fully open up. After I refresh the page multiple times it finally comes up but doesn't work properly. The map service I have added is not cached but I don't cache it for the Out of the box web application. The map service is relatively small geographically, and currently has no large feature classes in it or imagery. Would caching the map service change it that much?? I'd like to add imagery into but I want to get it working without it first. I would like to replace the web mapping application with the Flex viewer as there are more options. The data server and the web server are on the same machine. The map service is also using UNC paths as I am creating the original map on a different machine but if that is huge problem that could be adjusted. Trouble is the machine that has ArcGIS Server on it does not have Desktop and vice versa. Kind of limited with a single use license. Any help on any of this would be greatly appreciated.


title: Zoom function for add in in arcmap

text: I am completely new at gis but I have to make an add in to work with arcmap to make it zoom to a place I want it to go. (I use a textfield in a php website to tell where to go in arcmap) Is there a zoom function I can use in my case?


title: GDAL Create( ) for a netCDF file segfaults

text: I have a problem with the creation of NetCDF dataset using GDAL. When I use "gdalinfo" and "gdal_translate" utilities on a NetCDF format, they work fine, but when I create a new dataset using the function Create( ) member of the class GDALDriver (http://www.gdal.org/classGDALDriver.html) and passing the string format "netCDF", it says "Segmentation Fault". This instruction returns NULL: poDriver = GetGDALDriverManager()->GetDriverByName(pszFormat); with : const char *pszFormat = "netCDF"; It isn't a code error, because with other formats it works fine. I suppose it is a problem of libraries versions, I have: - gdal-1.10.1 - libhdf5-0-1.8.9-21.1 - libhdf5_hl0-1.8.9-21.1 - HDF4 4.2.9 - NETCDF 4.3.0


title: Calculate distances between polygon centres in QGIS

text: Is there a practical way to calculate distances between polygon centroids in QGIS? I have a UK city region map divided into administrative units (1 ward = 1 polygon). I want to represent their population growth in function of the distance to the core city and the best way to do it seems to be looking for the geometrical 'centroid' of each polygon and measure its distance to the 'centroid' of the core city polygon. I can do it manually with the ruler tool but it takes a long time and is not very accurate. Thanks


title: ArcGIS difference reclassify tool properties symbology

text: Is there a difference between the reclassify tool and just changing the ranges of a layer within the properties>symbology?


title: Run FME Script in batch file using .NET

text: Encountered another roadblock. I have a working FME script that I want to automatically run using a console app. I created a simple batch file which accepts a MDB file path as parameter:                 fme.exe E:\sample.fmw --SourceDataset_GEODATABASE_MDB "%1"      when I run this batch file, everything works well. But when I try to do the same thing only in my console application using the following code:               processInfo = New ProcessStartInfo("cmd.exe", "/c C:\executeFME.bat C:\sampleFile.mdb")     processInfo.CreateNoWindow = True     processInfo.UseShellExecute = False     processInfo.RedirectStandardError = True     processInfo.RedirectStandardOutput = True          process__1 = Process.Start(processInfo)     process__1.WaitForExit()          Dim output As String = process__1.StandardOutput.ReadToEnd()     Dim [error] As String = process__1.StandardError.ReadToEnd()          Dim exitCode = process__1.ExitCode                        Console.WriteLine("ExitCode: " & exitCode.ToString(), "ExecuteCommand"     process__1.Close()      I get the error " **Python Exception : No module named arcgisscripting** " I have a Startup python script that will check the geometries of the supplied mdb file and this is where it errors out. Do I still need to configure something to make it to work? Thanks. **EDIT** I noticed in the log files that the manually executed batch file using the correct Python Interpreter: **C:\Windows\SysWOW64\python27.dll** while the batch file executed by the console app uses a different one: **C:\apps\FME\fmepython27\python27.dll**. I tried modifying the Interpreter settings by using:               fme.exe APPLY_SETTINGS SYSTEM "Python/Python Interpreter" c:/windows/SysWOW64/python27.dll      Ran it again and it was still using the old interpreter. I think this is the cause of this but I'm out of ideas on how to fix it. Please help!


title: How to write in only a file .csv (or .txt) the values contained in several rasters?

text: I am an beginner in Python programming but I would like to write in only a file .csv (or .txt) the values contained in several rasters. For example, I built several rasters by using a loop. The rasters are like this:                outCostPath.save('P_' +  str(int(ID_patch_1)) + '_' + str(int(ID_patch_2)) + ".tif")      For each raster, I would to save the value that is situated at the second line and the column "PATHCOST" in a file .csv. In the end, this file .csv should contain all values contained in all created rasters. Here is my algorithm:               i = 1          for ID_patch_1 = 1 to 100             for ID_patch_1 = ID_patch_1 + 1 to 100              ## Building a raster         fileRaster = outCostPath.save('P_' +  str(int(ID_patch_1)) + '_' + str(int(ID_patch_2)) + ".tif")              ## Writing in a file .csv         fileCsv(i,1) = ID_patch_1 ## put value of ID_patch_1 in fileCsv at line i and first column         fileCsv(i,2) = ID_patch_2 ## put value of ID_patch_2 in fileCsv at line i and second column         fileCsv(i,3) = fileRaster(2,3) ## put value situated at line 2 and third column in the created raster         i = i + 1              end      end       The problem is that I don't know to code this algorithm in Arcpy. Thanks very much for your help.


title: How to add field to postgres layer with PyQgis?

text: I'm fetching a posgtres layer with PyQgis and then I'm trying to add field with this:               lyr.startEditing()     lyr.addAttribute(QgsField("class", QVariant.Int))      But It is not adding the field and also not giving any error. I manually added the field and inserted values programatically and it worked. So what am I doing wrong?


title: Poster or Map Presentation Guidelines?

text: Is there a website that will guide me through for various kind of paper size to create a/ or several presentation for digital or hardcopy maps/ or poster with different fonts, display area, title and so on. I am more interesting in reading on what is the recommend for the font size and the paper size.


title: How to find annotation overlaping file in arcmap 9.3?

text: I have all my textstring(Annotation) at distance 1.5 feet? I want to find all the annotation which at distnce less than 1.5? Also want to find the overlapping file? Can I make any file like shape file of so that I can zoom to location?


title: How often does Google change the geo-coordinates for a given address?

text: I geocoded the address "50 Hurt Plaza Southeast, Atlanta, GA 30303" twice. Once on July 28th, 2014 when I got the coordinates: (33.754208, -84.387272) And again today, August 18th, 2014, and got the coordinates: (33.754129, -84.38721). Should I expect google to be constantly updating its geo-coordinates for given addresses through the year? Or is this a rare occurence? It's a concern for me because I confirm the identity of buildings in the database I'm working with by comparing the geo-coordinates, and that's not going to work so well if they're in flux. Note the results above have been rounded to the sixth decimal. Any thoughts would be appreciated. Thanks!


title: Accessing NOAA CO-OPS water levels using SOS with bounding box and time extent

text: I'm tring to access all water level data within a bounding box and temporal extent using the OGC SOS from the NOAA CO-OPS. On the page here: http://opendap.co-ops.nos.noaa.gov/ioos-dif-sos/ there is a very useful SOS query builder tool, and I can formulate a query to get all the station data in a specific region: http://opendap.co-ops.nos.noaa.gov/ioos-dif- sos/SOS?service=SOS&request=GetObservation&version=1.0.0&observedProperty=water_surface_height_above_reference_datum&offering=urn:ioos:network:NOAA.NOS.CO- OPS:WaterLevelActive&featureOfInterest=BBOX:-74.4751,40.389,-73.7432,40.9397&responseFormat=text%2Fcsv&eventTime=2012-10-26T00:00:00Z/2012-11-01T00:00:00Z&result=VerticalDatum%3D%3Durn:ogc:def:datum:epsg::5103&dataType=VerifiedSixMinute But when I execute this query it tells me max of 1 day can be requested in a collection. So I guess there are two approaches. **Approach 1** : Do a query on collection with bounding box, but only the start time requested: http://opendap.co-ops.nos.noaa.gov/ioos-dif- sos/SOS?service=SOS&request=GetObservation&version=1.0.0&observedProperty=water_surface_height_above_reference_datum&offering=urn:ioos:network:NOAA.NOS.CO- OPS:WaterLevelActive&featureOfInterest=BBOX:-74.4751,40.389,-73.7432,40.9397&responseFormat=text%2Fcsv&eventTime=2012-10-26T00:00:00Z&result=VerticalDatum%3D%3Durn:ogc:def:datum:epsg::5103&dataType=VerifiedSixMinute which returns:               station_id,sensor_id,"latitude (degree)","longitude (degree)",date_time,"water_surface_height_above_reference_datum (m)",datum_id,"vertical_position (m)"     urn:ioos:station:NOAA.NOS.CO-OPS:8518750,urn:ioos:sensor:NOAA.NOS.CO-OPS:8518750:W1,40.7006,-74.0142,2012-10-26T00:00:00Z,0.446,urn:ogc:def:datum:epsg::5103,1.848     urn:ioos:station:NOAA.NOS.CO-OPS:8531680,urn:ioos:sensor:NOAA.NOS.CO-OPS:8531680:W1,40.4669,-74.0094,2012-10-26T00:00:00Z,0.319,urn:ogc:def:datum:epsg::5103,1.624      and then I could strip out the station_id and use them in single station queries. **Approach 2** : Issue multiple collection request for periods of one day and then concatenate the responses together. I guess I would use approach 1 if I had more days than stations, and approach 2 if I had more stations than days. Is there a better way?


title: How to set up an ArcMap Relate with a table from MS SQL Server 2012

text: This is my first time working with data directly from MS SQL Server 2012, and I'm unsure whether I can use it in a Relate in ArcMap. I have a table in MS SQL Server, which I brought into ArcMap via the Database Connections option. When I add this table to ArcMap using the standard Add Data button, a New Query Layer dialog appears: ![enter image description here](http://i.stack.imgur.com/taCqP.jpg) ArcMap requires a unique identifier field, however this table doesn't have one. I chose one of the existing fields, even though it's not unique (the first problem?). I can now view the table in ArcMap, and I'm able to set up a Relate between this table and a featureclass, based on the matching BoreholeID field. However, when I select a borehole from the featureclass, or a row in the table, the Relate is displayed but greyed out: ![enter image description here](http://i.stack.imgur.com/Dja2h.jpg)   * Can I use a table from MS SQL Server Management Studio in an ArcMap relate?   * If so, is the cause of the problem the lack of a unique ID field?   * If not, what else could be causing the problem?


title: How to get point file in 5km Buffer of route generated in pgrouting

text: I am using pgrouting to get shortest route. Now i want to get point that comes in buffer of 5km along router.Point is other shape file. route is using road shape so any one give me idea for same.


title: R with ggplot2/ggmap - How to load a png image as a map

text: I am trying to find a workaround for `ggmap`'s missing support of world maps (i.e. you can't create any maps that show latitudes > 80°, due to idiosyncrasies in the `mapproj` package). To a limited extend however, it seems possible to create empty world maps and save them as an image (png etc.), even if you can't use the ggmap object directly as one normally would in `ggmap(get_map(...))`. That's why I'd like to load a png (ideally, one I created with `ggmap`) into `ggplot2` and use that as a map instead. How exactly can I do that? I am aware that you can load background images im `ggplot2` (see this stackoverflow question). But I'd also like to plot points on my map, so it's important that the latitude/longitude values are mapped correctly.


title: How to change the Transparency of an ILayer?

text: I've created a function which sets the `IColor` to a specific `ILayer`, the value of transparency is between 0 and 255. When I set the value to 0, it gives me a transparent layer which is OK, but when I set it to 1,2,3,4,5,....255 it gives me the color as if it has 255 transparency. I don't understand, what's wrong with my code:               public static void setLayerColor(ILayer layer, IColor fillColor, IColor outlineColor, short width = 1)         {if (layer == null || fillColor == null || outlineColor == null) return;          IGeoFeatureLayer geoFeatureLayer = (IGeoFeatureLayer)layer;          ISimpleRenderer simpleRenderer = (ISimpleRenderer)geoFeatureLayer.Renderer;               ISimpleFillSymbol symbol = new SimpleFillSymbolClass();          ILineSymbol lineSymbol = new SimpleLineSymbolClass();               lineSymbol.Width = width;          lineSymbol.Color = outlineColor;               symbol.Color = fillColor;          symbol.Outline = lineSymbol;          simpleRenderer.Symbol = (ISymbol)symbol;}      My function has 4 parameters:   * `ILayer`, the layer of which I want to change its color.   * `IColor`, fillColor.   * `IColor`, outlineColor.   * `short`, the width of the outline.


title: How to add a sequence to the attribute table?

text: Is there a way to use the field calculator in QGIS to create a range of values in a count sequence with a specified input range? For example, if I wanted to have the values 1 to 5 repeated within my table according to the first record in the table. Thanks...


title: How are geographic coordinates stored for multivariate statistical analysis?

text: In terms of Principal Components Analysis or Multiple Discriminant Analysis, I have used dummy variables to represent a group of like-features that are located in the same geographic region. For example, in the analysis of two metropolitan cities, the census tracts that make up city one are classified as 0 and city two as 1. This binary classification groups the census tracts into two regions. I am wondering how I could store geographic coordinates (centroid?) of these census tracts so that I could introduce a spatial component to the multivariate analysis. This is dissimilar to the approach mentioned above since I am not grouping any features; but rather, I am interested in determining whether or not spatiality is a statistically significant independent variable (in terms of regression), or if spatiality is a significant part of a component/factor (in terms of PCA/FA). This would be simple enough if the variable was a single measured value; but alas, we have x and y coordinates to think about.


title: Contradictories results with Moran's I and Geary's C

text: I have contradictories results for both index. One side Moran's I says there is a positive correlation , but Geary's C says there isn't any correlation. ![enter image description here](http://i.stack.imgur.com/FF1Ci.png) ![enter image description here](http://i.stack.imgur.com/FIUye.png) Geary's C is more sensitive local variations, and Moran index is more global... wich of those should I consider? Supposedly both index should indicate negative correlation. Any idea? Thank you per advance! :D


title: Join based on maximum overlap in PostGIS/PostGresQL?

text: I have two sets of polygons in two tables. The sets overlap each other. For each polygon in set A, I would like to get the ID of the polygon in set B that it overlaps the most. I'm using PostgreSQL with the PostGIS extension. I know just enough about SQL to know that you can only join based on true/false conditions. So this won't work:               SELECT       a.id as a_id,       b.id as b_id,     FROM       a     JOIN       b     ON       max(ST_Area(ST_Intersection(a.geom, b.geom)))      because max() can't be in the ON clause. `ST_Intersects()` is a true/false test, so I **could** join on that, but polygons in set A will often overlap with more than one polygon in set B, and I need to know which one overlaps the **most**. ST_Intersects would presumably just return the first overlapping ID it came across, regardless of the extent of the overlap. This seems like it should be do-able, but it's beyond me. Any thoughts?


title: Is it possible to create filters based on attributes for markers in Leaflet?

text: I need to create filters based on attributes of markers. I have markers defined this way:               var locations = [         ['AB11 5HW','17','A',57.147701,-2.085442 ] ,     ['AB11 8DG','3','B',57.129372,-2.090916 ] ,     ];      If the user selects for 'A' or 'B' all markers containing attribute A or B should become visible on the map. Something like               var markersA =  a[2] contains 'A'       I thought it is maybe possible with Leaflet Marker groups or something similar. I'm new to Javascript so I'm interested if it's even possible.


title: QGIS howto create raster from map canvas easily

text: i have question ho to create raster from QGIS map canvas (what i seen in my qgis) easily. Now i have used print composer and after some luck i get whole dataset into map item, after that i throw lucky guess for dpi and i get something. What i would like to do is : make map look pretty in qgis, select cell resolution and get what i see in georeferenced raster file (geotiff, png etc..) from whole dataset. is this possible or do i just use print composer wrong ?


title: GeoServer Date ftl showing time stamp rather than just the date

text: I had an issue with Geoserver showing my postgresql data in the US (MDY) format rather than the UK / European format (DMY). I fixed this by removing the comment #, and setting the language, to en_GB. In the tomcat6.conf file.               # You can change your tomcat locale here     #LANG="en_US"      The data now displays on my geoserver outputs as DMY, however HH:MM is also shown when i do a get feature request. The postgresql field is Date, and Geoserver lists it as a Date in the layer information panel. However i have 00:00 after the date when i click for a GETfeature information request. **How would I stop 00:00 after the date?**


title: Load locations faster in ArcGis network analyst

text: I am currently using ArcGis Network Analyst to calculate an OD cost matrix. I was wondering whether there is a faster way to load my locations (origins, destinations and barriers) because my current rate is around 10 locations loaded per second (which is quite slow for I have about 7.000 points to load). Moreover, is there a way to save these locations in order to use them later? I think of origins/destinations but also of the barriers. I didn't find anything about this in the ArcGis documentation. Thanks for your help!


title: How to Configure sextante in Qgis 1.9 ubuntu

text: I am using Ubuntu 13.04 and the most up to date version of Qgis 2.1.0-Master. I am trying to get sextante to work correctly and have a hard time. Specifically SAGA 2.1.0 and GRASS 6.4.2 will not execute correctly. When I try to run any of the grass tools qgis hangs then crashes. When I try to run SAGA tools it tells me saga is not installed or needs to be configured. The latest version of Sextante seems to have removed the SAGA folder input field which I am assuming needs to be updated. I am not sure why grass is crashing at all. I have another similar question open HERE. Any help would be greatly appreciated.


title: MMQgis Plugin error in Qgis

text: I am using Qgis 1.9 .Just now I downloaded a plugin "mm qgis"But When I run this plugin I am getting following error: An error has occured while executing Python code: Traceback (most recent call last): File "C:\Documents and Settings\15491/.qgis2/python/plugins\mmqgis\mmqgis_dialogs.py", line 525, in run message = mmqgis_geometry_convert(self.iface, layername, newtype, splitnodes, savename, 1) File "C:\Documents and Settings\15491/.qgis2/python/plugins\mmqgis\mmqgis_library.py", line 1154, in mmqgis_geometry_convert layer.dataProvider().select(layer.dataProvider().attributeIndexes()) AttributeError: 'QgsVectorDataProvider' object has no attribute 'select' Python version: 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)] QGIS version: 1.9.0-Master Master, 1a723b4 Python path: ['C:\Documents and Settings\15491/.qgis2/python/plugins\sextante_taudem', 'C:\Documents and Settings\15491/.qgis2/python/plugins\sextante_animove', 'C:\Documents and Settings\15491/.qgis2/python/plugins\sextantelwgeomprovider', 'C:\Documents and Settings\15491/.qgis2/python/plugins\sextanteexampleprovider', 'C:\Documents and Settings\15491/.qgis2/python/plugins\sextante', 'C:\Documents and Settings\15491/.qgis2/python/plugins\permaclim', 'C:/OSGeo4W/apps/qgis-dev/./python', 'C:\Documents and Settings\15491/.qgis2/python', 'C:\Documents and Settings\15491/.qgis2/python/plugins', 'C:/OSGeo4W/apps/qgis- dev/./python/plugins', '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\bin', '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\wx-2.8-msw-unicode', 'C:\Documents and Settings\15491\.qgis2\python\plugins\mmqgis/forms', 'C:\OSGeo4W\apps\qgis- dev\python\plugins\fTools\tools', 'C:\Documents and Settings\15491\.qgis2\python\plugins\pdokbaggeocoder/forms'] what is the mistake?


title: SpatiaLite: Creating to, from nodes from a normal road layer (LINE)

text: Following this tutorial but using my own data. How do I create the from_node and to_node columns? I am using Spatialite and Sqlite3. I'd prefer to not use a gui as I'm working remotely via ssh on a server. The network is currently held in a spatialite database with a geometry column e.g:               id  |  road_name   |  geometry     1   |  Phil Street |        I wish to, using sql, add column to and from node e.g.:               id  |  road_name   |  from_node  |  to_node  |  geometry     1   |  Phil Street |  1          |  2        |


title: How to calculate distance between point feature to multiple polyline features in ArcGIS

text: I have a GPS track which is basically a point layer. I have different route networks such as - train, tram and bus routes. I want to calculate the distance between each point to nearest features within a search radius. I used Generate Near table to calculate distance. This works fine. However after getting the table I have to sort the distance for each and every point. I am looking for a method which will directly generate the distance result as follows- pointID | near_feature_1_distance | near_feature_1_distance | near_feature_1_distance where pointID is the GPS points, near_feature_1_distance is the distance from that GPS point to a given near_feature_1 (say tram line) and so on. I did this in ArcGIS 9.3 long back. But now I am using ArcGIS 10.1 (however I have 9.3 as well) but I cannot remember how I did it. Any help will be highly appreciated.


title: SLD encoding in GeoSever on Ubuntu

text: I am facing a challenge with the encoding in my SLD file in GeoServer on a Ubuntu server. When I create a SLD file containing the danish letters æ, ø, å in the GeoServer web interface like this filter:               <ogc:Filter>      <ogc:PropertyIsEqualTo>        <ogc:PropertyName>priority</ogc:PropertyName>        <ogc:Literal>Høj</ogc:Literal>      </ogc:PropertyIsEqualTo>     </ogc:Filter>      In the "Style Editor" the "ø" looks fine and when submit I can see (in Firebug) that the posted request contains the letter "ø". But when I try to edit the style again, my filter looks like this:               <ogc:Filter>      <ogc:PropertyIsEqualTo>        <ogc:PropertyName>priority</ogc:PropertyName>        <ogc:Literal>H?j</ogc:Literal>      </ogc:PropertyIsEqualTo>     </ogc:Filter>      The "ø" has been replaced with a "?". My SLD header looks like this:               <?xml version="1.0" encoding="UTF-8"?>      I have tried to use Latin-1 encoding               <?xml version="1.0" encoding="ISO-8859-1"?>      But I the "ø" still shows as a "?". When I run "env | grep LANG" in the terminal on the Ubuntu server I get this output:               LANG=en_GB.UTF-8     GDM_LANG=en_GB     LANGUAGE=en_GB:en      This should indicate the default encoding is utf-8. When I run a "file --mime- encoding myld.sld" I get:               us-ascii      My question is: **How do I get GeoServer to write the sld file in utf-8 encoding or in Latin-1 so I can use the danish letters in my SLDs?**


title: OpenLayers map breaks at certain zoom level when many popups exist in the map

text:
I'm still learning OpenLayers. I found this weird issue in my OpenLayers map. I create FramedCloud popups and add them to each features.               popup = new OpenLayers.Popup.FramedCloud("popup",         OpenLayers.LonLat.fromString(parsedFeature.geometry.toShortString()),         null,         "
date: " + someDate + "
longitude: " + data.lon + "
latitude: " + data.lat + "
",         null,         false     );      But I don't add them until it gets to certain zoom level. Then I add the popups to the map with the following code               olMap.events.register('zoomend', olMap, function () {         var zoom = olMap.getZoom();         var popup;         if (zoom > 11) {             for (i = 0; i < dm.currentLocationsVector.features.length; i++) {                 popup = dm.currentLocationsVector.features[i].popup;                 if (!popup.map) {                     olMap.addPopup(popup);                 }             }         } else if (zoom < 12) {             for (i = 0; i < vector.features.length; i++) {                 popup = vector.features[i].popup;                 if (popup.map) {                     olMap.removePopup(popup);                 }             }         }     });      Actually it works good until zoom 17 then when the map zoom level gets 18 the map breaks. ![enter image description here](http://i.stack.imgur.com/HLmiZ.png) It looks like the popups and the map get messed in some way. I don't think I have problem in creating popups.Thought it might some bug in OpenLayers. Cheers


title: Has Export to AI (Adobe Illustrator) changed at ArcGIS for Desktop 10.2?

text: In the past when I would export vector layers from ArcGIS each layer would be nested in a group with a clipping path over the top. Now when I Export to AI from 10.2 that layer's features are 5-6 layers deep. Did this software behaviour change at 10.2? ![enter image description here](http://i.stack.imgur.com/eDB0Q.jpg)


title: What's the correct transformation for Lambert Conformal Conic for Mexico?

text: I have a shp archive without the .prj file. I don´t know the proyection of this shp, but I think it is in lambert conformal conic. Can anybody helpme to proyect this shp? i´ve tried to proyect it in wgs84 utm 16N and geografic wgs84, but i have no success. this shp is from Quintana Roo, México. I´m working in qgis


title: How to Import Metadata of feature classes with them into ArcSDE?

text: I have a database that has numerous feature classes of different types. None of these feature classes have metadata. I understand that when you import feature classes to a database, the metadata does not come with. I'm trying to import metadata to features that have never had metadata and I'm not seeing a way to do this.


title: OpenLayers Mapserver spatial filter INTERSECTS format problems

text: Trying to get hand-drawn polygon delimiter for spatial query via OpenLayers to mapserver (then postgis). I can get the WKT of the drawn polygon and can substitute the dynamic values, but am having trouble getting the spatial query properly formatted. Filter defined in OpenLayers:               myfilter = new OpenLayers.Filter.Logical({                 type: OpenLayers.Filter.Logical.AND,                 filters: [ new OpenLayers.Filter.Spatial({                    type: OpenLayers.Filter.Spatial.INTERSECTS,                    property: "footprint",                    value: curr_wkt                    }),                    new OpenLayers.Filter.Comparison({                       type: OpenLayers.Filter.Comparison.EQUAL_TO,                       property:  "roll",                       value: roll                     })                   ]              });      What myfilter looks like after created:               (INTERSECTS(footprint, '(POLYGON((-161.6625 71.073751831055,-167.625 69.498751831055,-168.3 65.673751831055,     -166.275 61.848751831055,-157.6125 64.773751831055,-161.6625 71.073751831055))')) AND (roll = '4')      But in the POST output from Firebug when run, the filter gets converted to an undefined Box (or Envelope if using 1.1.0): POST result - footprint undefined undefined undefined undefined roll 4 Is the WKT the wrong input format to use? Why does it change it to Box?