title: Is it safe to store coordinates on NFC tags always in one CRS and rely on conversion for others

text: I'm involved in NFC (Near Field Communication) standardization of a record structure for geographic point locations. That will allow to put coordinates on NFC memory tags which can be read with a smartphone or to exchange between phones and the like. A first proposal made was to always include the coordinate tuple (2D or 3D) and the CRS, so a coordinate tuple for a specific CRS could be stored. To this, an alternative opinion claimed that the CRS could set implicitly (for example WGS-84) and coordinates for any other CRS be computed, of course only for the same dimensions. It may be worth to note that a compact encoding format is important because NFC tags, especially the cheaper ones, do not provide a lot of memory. So any byte less is good if functionality isn't reduced (one envisioned use case is that NFC tags could be affixed to landmarks and allow the location be read electronically - this was thought to sometimes require other CRS than WGS-84). So the question is, if a coordinate tuple is given in one CRS is it safe to assume that it can be converted to any other CRS of practical relevance? To make this more open-ended, would there be any drawbacks if a data format for geographical coordinates is fixed to only one CRS and requires conversion to and from other CRSs if needed. And as further explanation, the intention is not to derive a format for a specific application but for general use by anyone who wants to store or read a location from an NFC tag/card or share with another smartphone when touched (like to show the location on a map application). Similar to what RFC 5870 (Geo URI) provides but more compact.


title: Custom Tools not Saving in Toolbar on Startup/ Normal.mxt?

text: I am trying to save some tools that I have created into the default toolbar layout that is present when ArcMap 10.1 is launched. It is possible to change the arrangement of toolbars, add new ones and even add in existing tools, but any custom tools I set in the commands category "geoprocessing tools" are gone and do not appear in the toolbar. How can I save these custom tools to the default layout that appears on startup? **Update:** Saving a document with the custom tools added and reopening it preserves the tools, but opening a template (10.1 mxd) containing the tools, then opening and older document will cause the tools to disappear. Oddly enough, any other changes made after the custom tools were added are preseved. **Update (03/18/2014):** I have upgraded to 10.2 and this has still not been fixed. I did not see this issue in the 10.21 changelog so I doubt it will address it. For those interested in this issue, I have started a forum post on the ESRI forums here.


title: QGIS and ArcGIS on the same PC

text: I work for a small non profit conservation organization and we have a reasonably high turnover of people... so it stands to reason that some people don't want to change from a preferred program. this in mind... Can i install ArcGIS 10.1 and QGIS on the same PC without problems?


title: Convert MapInfo raster file (.tab) to use it in QGIS and/or ArcGIS

text: I am trying to convert raster files from MapInfo (.ecw, .jpg, .TAB) in order to use them with QGIS or ArcGIS. I have access to all 3 softwares. I have tried Data interoperability tools in ArcGIS but it is not working ("Tool execution failed."). I have also used a MapInfo tool (.mbx) designed for this purpose (http://georezo.net/wiki/main/logiciels/mapinfo/outils/ut4rt) with more success. I can then open the file with QGIS, but the output file is a bit off (even with the right coordinate system). Any ideas ? Sincerely


title: JS API WebApp Zoom Problem in IE8

text:
I implemented a zoom to feature extent function in my code, it works fine in new FireFox, Chrome, and IE10, but the tester's are having troubles with it in IE8. I noticed that all the data shows up on the map and in the table, yet it doesn't zoom. What am I doing wrong? Here is where I call my zoom function.               function displayResults(results, messages) {         var featureset = results[0].value;              var simplePolySymbol = new esri.symbol.SimpleFillSymbol();         simplePolySymbol.setOutline(new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID, new dojo.Color([0, 96, 170, .75]), 1));              dojo.forEach(featureset.features, function (feature) {             feature.setSymbol(simplePolySymbol);                  resultsLayer.add(feature);         });              myMap.addLayer(resultsLayer);         myMap.Search.handleStats(results, messages);              var extent = esri.graphicsExtent(featureset.features);         myMap.setExtent(extent, true);     }      Also my table doesn't populate in IE8, althought im seeing the values I need in the console. Here is my table code:               myMap.Search = {         gpGetStats: new esri.tasks.Geoprocessor(gpURL),         getRouteStats: function () {             console.log(dojo.byId('resultsLayer').value);             var params = {                 "ZIP_CRID": dojo.byId('resultsLayer').value             };             resultsLayer.Search.gpGetStats.execute(params, resultsLayer.Search.handleStats);         },              handleStats: function (results, messages) {             console.log(results);                  var content = '';                  content += '
';             content += '
 ZIP CRID 	 Bus 	 Res 	 Tot 	 <200 	 DS Key
';                  dojo.forEach(results[0].value.features, function (crid, i) {                 console.log(crid.attributes.ZIP_CRID + ', ');                 content += '
' + crid.attributes.ZIP_CRID;                 content += '	 ' + crid.attributes.BUS_CNT;                 content += '	 ' + crid.attributes.RES_CNT;                 content += '	 ' + crid.attributes.TOT_CNT;                 content += '	 ' + crid.attributes.LT_200_IND;                 content += '	 ' + crid.attributes.DS_KEY + '
';             });             content += '
';             dojo.byId('details').innerHTML = content;         },


title: Default Opengeo suite print module- scales given in selector differ when using the zoom slider

text: Using a default installation of the 3.0.2 suite, opening the Composer view(with the MapQuest OSM tiles as background) and clicking print (to pdf) to activate the print preview there appears to be a couple of issues with the scale given in the scale selector box. _One_ \- if one slides the zoom slider so that the scale selector box gives 1:512 000, and then opens the scales list in the box and chooses 1:512 000 the scale visibly changes in the print preview map (it goes to the next smaller scale). _Two_ \- the scales in the scales list in the scale selector box do not entirely match the list of scales given in config.yaml: the smallest (1:524288000) is not listed. _Three_ \- there are 17 zoom levels in the default Composer map, 18 in the print preview's zoom slider, and 19 in the scale selector box. The scales given by the zoom slider in the default Composer map roughly match the scales given in Config.yamls, i.e 1:68k something becomes 1:64k etc, but this should not be reason for this default behaviour? I have tried changing the number of scales config.yaml as well as changed the maxZoomLevels for the map, but besides zoom slider in the print preview following the number of zoom levels (scales) of maxZoomLevels, even if the main map's zoom slider does not, I cannot get things to precisely line up. A search here and abroad did not reveal this issue described, but if parts of this (or all) is due to a lack of understanding of certain fundamentals on my side, then pray tell where I can go read :) Many thanks for reading.


title: Extract BoundingBox of Cities

text: For simple geocoding purposes I want to extract all cities from a single country of openstreetmap. They only queries I want to do are:   * ZipCode -> (List of) BoundingBox of City/Village   * City/Villagename -> (List of) Boundingbox of City/Village I'm not interested in streetnames, rivers etc. Is this possible using `osmfilter`. Can anyone point me towards the tags of the data i'm looking for?


title: Programmatically manipulate GPX data

text: I need to do the following in Python   * open OSM & GPX files (I have packages for this)   * transform points from GPX (like stick a track to roads)   * calculate the results (distances, cumulative distances, etc.) I need to do this repeatedly, so a Python script is much more preferred. Desktop software with plugins is not suitable. I'll prefer some Python & C modules than a quest of installing plugins. PostGIS may be an option too. Shapely (Python package) seems to not be able to do this (it works only in 2D on a plane, and mentions it has no projections). I don't need 3D, but I have the input as lat&lon coordinates, and need to do geometric transformations (project a point on a polyline) and calculate distances in metres. What modules should I use?


title: Using Custom Map Projection On Leaflet

text: I have an application displaying a EPSG:3414 map running on Openlayers - jsfiddle working sample that works. I am trying to implement the same on Leaflet - jsfiddle non-working, it is not working as I may have left out something or may have gotten some concept wrong. Unfortunately, I am quite unfamiliar and am in the process of figuring out how to get this to work. Any help would be welcome. Thanks! **Updates:**   1. map is ESRI type using ZYX and SVY21 projection):   2. updated code here (JS fiddle not updated) **Leaflet Sample:**               var proj = new L.Proj.CRS.TMS('EPSG:3414', "+proj=tmerc +lat_0=1.366666666666667 +lon_0=103.8333333333333 +k=1 +x_0=28001.642 +y_0=38744.572 +ellps=WGS84 +units=m +no_defs",         [4257.9355, 11470.3186, 64359.3062, 49067.5413],          {             resolutions:[76.4372195411057,38.2186097705529,19.1093048852764,9.55465244263822,4.77732622131911,2.38866311065955,1.19433155532978,0.597165777664889,0.298450596901194]         }     );          var map = L.map('my_map', {         crs: proj,         continuousWorld: true,         worldCopyJump: false     });          map.addLayer(new L.Proj.TileLayer.TMS('http://{s}.onemap.sg/ArcGIS/rest/services/basemap/MapServer/tile/{z}/{y}/{x}', proj, {         subdomains: ['t1', 't2'],         tms: true,         continuousWorld: true     }));      Update, more information to help. I used openlayers to implement, have not had time to try out leaflet.               Spatial Reference: PROJCS["SVY21",GEOGCS["SVY21[WGS84]",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",28001.642],PARAMETER["False_Northing",38744.572],PARAMETER["Central_Meridian",103.8333333333333],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",1.366666666666667],UNIT["Meter",1.0]]          Single Fused Map Cache: true          Tile Info:              Height: 128         Width: 128         DPI: 96         Levels of Detail: (# Levels: 9)             Level ID: 0 [Start Tile, End Tile]                 Resolution: 76.4372195411057                 Scale: 288896              Level ID: 1 [Start Tile, End Tile]                 Resolution: 38.2186097705529                 Scale: 144448              Level ID: 2 [Start Tile, End Tile]                 Resolution: 19.1093048852764                 Scale: 72224              Level ID: 3 [Start Tile, End Tile]                 Resolution: 9.55465244263822                 Scale: 36112              Level ID: 4 [Start Tile, End Tile]                 Resolution: 4.77732622131911                 Scale: 18056              Level ID: 5 [Start Tile, End Tile]                 Resolution: 2.38866311065955                 Scale: 9028              Level ID: 6 [Start Tile, End Tile]                 Resolution: 1.19433155532978                 Scale: 4514              Level ID: 7 [Start Tile, End Tile]                 Resolution: 0.597165777664889                 Scale: 2257              Level ID: 8 [Start Tile, End Tile]                 Resolution: 0.298450596901194                 Scale: 1128          Format: JPEG         Compression Quality: 80         Origin:             X: -5878011.89743302             Y: 10172511.897433         Spatial Reference: PROJCS["SVY21",GEOGCS["SVY21[WGS84]",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",28001.642],PARAMETER["False_Northing",38744.572],PARAMETER["Central_Meridian",103.8333333333333],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",1.366666666666667],UNIT["Meter",1.0]]          Initial Extent:              XMin: -3250         YMin: 13872.0314312493         XMax: 68250         YMax: 56958.47469375         Spatial Reference: PROJCS["SVY21",GEOGCS["SVY21[WGS84]",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",28001.642],PARAMETER["False_Northing",38744.572],PARAMETER["Central_Meridian",103.8333333333333],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",1.366666666666667],UNIT["Meter",1.0]]          Full Extent:              XMin: -3250         YMin: 13872.0314312493         XMax: 68250         YMax: 56958.47469375         Spatial Reference: PROJCS["SVY21",GEOGCS["SVY21[WGS84]",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",28001.642],PARAMETER["False_Northing",38744.572],PARAMETER["Central_Meridian",103.8333333333333],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",1.366666666666667],UNIT["Meter",1.0]]          Units: esriMeters          Supported Image Format Types: PNG24,PNG,JPG,DIB,TIFF,EMF,PS,PDF,GIF,SVG,SVGZ,AI,BMP


title: Satellite/Aerial Imagery for Nepal for mapping in OpenStreetMap

text: I was trying to map places in Nepal using satellite imagery but there is no high resolution imagery available at zoom level 15-16 . I was looking to map places like Illam, Dhankuta. Is there a free source of aerial imagery for these places so that i can use it with JOSM to map these places in openstreetmap?


title: Android app to save offline map + POIs?

text: I don't know much about GIS and would like to find an application for travelling with the following features:   * support Points of Interests (sights, restaurants, etc.); Clicking on any POI should display a pop-up with a bit of information (eg. address, history of a building, etc.)   * download maps for offline use, so that no data connection is required   * should display the user's current location using GPS. Apparently, a data connection is not required to use GPS, although it might speed things up when enabling GPS Are there Android applications that can do this? That would make it easier to travel rather than carrying infos on paper. Thank you. \-- Edit: Android Google Maps doesn't support layers for offline maps, so adding one's locations through the web Google Maps won't work. Users recommended MapsWithMe which is supposed to support "bookmarks" in the Pro version. If someone uses the Pro version, can you confirm that locations are downloaded and saved in offline maps? Also, is it possible to build a list on a PC through eg. Google Maps/MapsEngine, and tell MapsWithMe Pro to download it? E-mails to the author (info@, support@, and android@) returned no reply.


title: How do I open the GeoDB files on NaturalEarthData's site without ArcGIS?

text: I'm trying to use the admin1 dataset from the NaturalEarthData website. The site offers this data as both shapefiles and GeoDB. The shapefile attribute data is apparently encoded in Windows-1252 and a lot of the characters in some of the attribute data I need are just '?' characters. The NaturalEarthData site mentions that:               The GeoDB format preserves local scripts better than SHP format.      I built GDAL against the FileGDB SDK from ArcGIS and verified the GDAL FGDB driver is working correctly using this link: http://trac.osgeo.org/gdal/wiki/FileGDB But I can't open files from NaturalEarthData. I always get the following error:               Unable to open datasource `ne_10m_admin_1_states_provinces.gdb' with the following drivers     ...     FileGDB     ...      I don't have access to ArcGIS. Am I doing something wrong? Is there any way I can access the GeoDB data?


title: ArcGis File Geodatabase with Null values (Field Calculator)

text: I want to use the field calculator to get the sum of a couple of columns. The problem is that some of them has Null Values. I am not able to find a solution. I tried this Python functions:               stack( !basic.Blue_Thickness!, !basic.Green_Thickness!, !basic.Red_Thickness!, !basic.F1_Thickness! )      Version 1:               def stack(*args):       myList = list(args)       myList = [item for item in myList if (item is not None)]       return sum(myList)      Version 2:               def stack(*args):         return sum(filter(None, args))      Version 3:               def stack(*args):       myList = list(args)           for item in myList:         if item is None:            myList.remove(item)       return sum(myList)      The result of all this functions is that I get always NULL back, only if there is a row with no NULL then I get a result. Background: I have a File Geodatabase with one main table and some tables which joins this table. The calulation takes place only in the main table (source columns and write column). I am using ArcGIS 10.0. The Source of these functions is this discussion: calculating a field in which null values may be present I would be really happy if somebody could help me. cheers Immanuel


title: Tool/Algorithm for generalizing OSM-Files for printing maps

text: Is there a tool or algorithm to (topological) generalizes OSM-shapes? I would like to build a printable map in different scales. In the image below you can see the same divided road (yellow) at different scales. ![enter image description here](http://i.stack.imgur.com/PPrTd.png)


title: How to Change the schema of SDE dataset in Connected environment

text: We continuously use to receive changes in Schema ( field +/-, domain changes etc ) for our Utility dataset.. Since we have swicthed to a multiuser Geodatabase ( ArcSDE on Oracle ) , we would like to make apply this schema changes without affecting our current workflow. How can we change\update the new schema on our existing SDE dataset , with minimal stoppage of the workflow ? One easy way would be to remove the existing dataset and configuring the new dataset ( which has new schema ) but wouldnt that make the existing completed jobs incompatible with the new schema ? A quick response from experts highly appreciated.


title: ARCGIS10: updatelayer results in error - LayerObject: Get attribute renderer does not exist

text: I'm quite new in using python in Arcgis 10. I wanted to use the symbology of a geostatistical layer in the TOC in another geostatistical layer in the TOC. Therefore, one should use arcpy.mapping.updatelayer. I saved the geostatistical layer as Layer file (TensorStyle.lyr). The code to copy the symbology looks like               import arcpy          mxd = arcpy.mapping.MapDocument("CURRENT")          df = arcpy.mapping.ListDataFrames(mxd)[0]          updateLayer = arcpy.mapping.ListLayers(mxd, "*", df)          sourceLayer = arcpy.mapping.Layer("TensorStyle.lyr")          for lyr in updateLayer:              arcpy.mapping.UpdateLayer(df, lyr, sourceLayer)      If I run the script I get the message: > Runtime error : LayerObject: Get attribute renderer does not exist Could you please help me. Thank you in advance, Ronny


title: Compile a plugin in 1.9 Qgis-dev

text: I am actually learning qgis and how to devellop a plugin. I have to write this plugin in C++ and not python because people coding in C++ want to modify it after. I have found some plugins in the git repository in C++ and I manage to compile it with 1.9 librairies on Qgis. But it's seems they change the constructor of the QgisPlugin. QgisPlugin( name_, description_, version_, type_) all the parameters are QString in 1.9 Qgis but "type_" was a QgisPlugin parameter before. I put a string on the constructor and I could compile the plugin in ".dll" but qgis 1.9 does'nt recognize my plugin. Do you know where I can found the source of a plugin made in 1.9 ? Thank you if you can answer my question Regards


title: How to display partially a WFS layer? (just the content within a box)

text: I need to know how to display a WFS layers, but not the whole content, only the part within a rectangle. I have this code:               var SMap = new OpenLayers.Layer.Vector("Status map", {         strategies: [new OpenLayers.Strategy.BBOX()],         protocol: new OpenLayers.Protocol.WFS({             url: "http://demo.opengeo.org/geoserver/wfs",             featureType: "states"         }),     });     SMap.setVisibility(true);     map.addLayers([SMap]);      And I also have two coordinates from this code:               controlS = new OpenLayers.Control();     OpenLayers.Util.extend(controlS, {         draw: function () {             this.box = new OpenLayers.Handler.Box(controlS,                 {"done": this.notice},                 {keyMask: OpenLayers.Handler.MOD_SHIFT});                 this.box.activate();         },         notice: function (bounds) {             var ll = map.getLonLatFromPixel(new OpenLayers.Pixel(bounds.left, bounds.bottom));              var ur = map.getLonLatFromPixel(new OpenLayers.Pixel(bounds.right, bounds.top));              alert(ll.lon.toFixed(4) + ", " +                    ll.lat.toFixed(4) + ", " +                    ur.lon.toFixed(4) + ", " +                    ur.lat.toFixed(4));         }     });     map.addControl(controlS);      What can I do to display only the WFS layer's content within the imaginary box formed by 2 pairs of coordinates? Thanks and regards, Rafael


title: Graphic drag and drop using ESRI Javascript API 3.7 in IE10

text: I am using ESRI JS API 3.7 and Dojo 1.9. I have a problem with graphic DnD in IE9 and 10, but not IE8.               this._listeners.connect("dropListener",this.esriMap, "onMouseUp", d_lang.hitch(this, function(evt) {         // some code     }));      I registered an event handler for onMouseUp to retrieve the coordinate of the location where the graphic is dropped. In IE, onMouseUp is fired only when I clicked on the map. When I dragged the graphic, the _pointerUp method in PointerEvents.js was called, but then it stopped there. No further events are fired until I clicked on the map, where the _pointerUp method will be called and in the end, onMouseUp will be fired. Any ideas how I could retrieve the location information without first clicking on the map? Thanks in advance.


title: How do I find the intersection of two polylines in PostGIS?

text: I have a set of road conditions reported for small segments of road, and I am hoping to determine the specific subset of each segment which intersect with a route along the road so that I can draw the intersection of the route and road conditions on a map. Is there a way to determine the specific points of intersection of two polylines in PostGIS in a way suitable for plotting on a Google Map? I am new to PostGIS so please accept my apologize if I am getting any terminology wrong.


title: Career advice in Gis field

text: I am a B.Tech graduate(2013) in Geo Informatics and Geo Technology.Presently, i am working i an company as Junoir Executive in Gis . I also know some programming languages like "V.B.net" and C#.My doubt is how can i improve my Career strength in this field so that i can earn handful salary ?


title: OpenLayers - How do I change the WMS URL in Code

text: I want to be able to change the url for a layer in code. The reasons is that I am implementing local caching. When local caching is Not enabled, we use 4 servers to deliver tiles to run parallel http requests. However, when the caching is enabled, it caches tiles for the four servers. There is no guarantee that the request for the same tile will be for the same server, so, potentially, the local cache may end up four times the size it needs to be. (and four copies of the same tile will be in the cache) To counter this, we have set up another server to be used only for users who are using the local cache to guarantee that the cache will always return a tile if it is there (simply because the server is always the same) I have experimented with layer.getFullRequestString( { }, url ); I am able to actually retrieve the request string if I use this syntax: layer.getFullRequestString( { } ) (note, it errors if you do not put the {} ) but have not been able to actually cause a change in the layers fetch url. The function seems to ignore the supplied url completely. I have tried a single url and a url array to no avail. Openlayers Documenmtation is as follows:               getFullRequestString:function(  newParams, altUrl   ) Combine the layer’s url with its params and these newParams.      > Add the SRS parameter from projection -- this is probably more eloquently > done via a setProjection() method, but this works for now and always. > > Parameters newParams {Object} altUrl {String} Use this as the url instead of > the layer’s url Returns {String} Any thoughts or better ways to address this very much appreciated Thanks


title: How to work with tables in MapBasic?

text: I try to get data from table use MapBasic:               Dim nCol = "ColumnName"     Fetch First From Untitled     Do While Not EOT(Untitled)       colValue = Selection.nCol     Note colValue     Fetch Next From Untitled     Loop      Compilation ends success but when i start application in MapInfo get error:               Not defined Selection.nCol      How to work with columns in this case?


title: How to select a polygon feature on mapservice?

text: I manage to publish a mapservice with lots of polygon. I wish to select a polygon and see its attributes. May I know how to I perform that with javascript API? There doesn't seems to be a straight forward solution on their sample page.


title: Best way to group datasets by country/region

text: I 'm using ArcGIS 10.0 and have a folder with misc GIS datasets from all over the world. I wrote a script that finds all of these datasets and writes them to a .csv file, but I also want to identify what country/region that dataset falls into and write it to the CSV. I thought of pulling the extent and comparing that to country specific extents would be one way to group the datasets but there would be overlap causing some data to grouped incorrectly. Any suggestions besides doing desc.extent need to be very timely since I'll be sorting through lots of data. Thanks for your help.


title: Why does ArcGIS for Desktop legend not always show symbol as desired?

text: I have a problem with the legends in the ArcMap because the patch does not show the marker symbol inside it well. Whenever, I move the legend it changes but never to what I want for every patch. I have tried legend properties, and so on but nothing helped me. You can see my legend **here**.


title: How to return values of its order?

text: I have a matrix like this:                  10 0.646519601   0   0.085929349 0.255477717 0.012073334        40 0.138403594   0   0.065979943 0.776693742 0.018922721      I want to keep the first column, and get rest numbers sorted by row. The result should be like this:                  10 5 1 2 4 3  #0.64 is the largest; 0 is the smallest in the row#        40 4 1 3 5 2


title: How to create filled contour lines in QGIS?

text: How can I create contour lines that are filled with a graduated symbol usinq QGIS? In ArcGIS for Desktop I can do **this**. What about in QGIS?


title: How to make grid lines opaque so just lat/long coords are on map border?

text: Can any one help me with adding graticules to my map? I have managed to add a lat/long grid over my map, however I want the grid lines to be opaque with just the lat/long coords on the map border - is this possible? I am using QGIS 1.8


title: How to add a new data frame to an MXD using C#

text: Does anyone have code in C# to add a new data frame (map) to an MXD? I have an MXD created in ArcMap 10.1 that has an existing data frame. I need to programmatically add a new data frame and populate it with several layers from a personal geodatabase. I'm using VS2010 with .Net 3.5. I have created a MultiItem that reads the names of folders that contain the PGDB files. After selecting a name from the list (which is populating correctly) it triggers OnSelChange. In that routine I need to create a data frame named based on the folder name selected from the MultiItem list. The name comes in just fine, but I'm stuck creating the new dataframe. I tried using IMaps.Create() but it doesn't add anything to the count of maps. I have tried then using IMaps.Add() but VS reports that Add is not defined for IMaps contrary to being described in the ArcObjects reference. The only procedure I have found is the one written in VBA using ArcObjects 10. It seems rather convoluted given there are supposed to be functional methods for IMaps to do exactly what I need. Any help would be greatly appreciated.


title: Splitting a bounding box into n amount of bounding boxes

text: Firstly, I've never worked in the GIS area before nor do I have much knowledge of the mathematics behind some of the stuff I am trying to implement. Perhaps i'm getting in over my head. **My Situation** I am creating a geospatial map viewer in Qt. So far I have a MapTile class which stores its center location (lat and long) and the width of the box in kilometers, which may or may not be a good idea. From this information I can construct a bounding box from the center point and the width by using the algorithm found here. The MapTiles are all square, no irregular shapes or rectangles. Is this a bad implementation? **My Problem** Now I am stuck with implementing my MapTileMatrix class which stores a multidimensional array of MapTiles[rows][cols]. I figured I could tell the MapTileMatrix what its bounding box is, and then have it generate the MapTiles position and width from this data. I haven't been able to find any algorithm online or mention of how to split a bounding box down into n number of sub- boundingboxes. Is there a method of dividing a bounding box into say 4 equal parts? What about 9, 16 or 25? If this is not possible perhaps I should reconsider my design. **On a side note:** Are there any guides / tutorials to implementing a map viewer application? I have read over the Web Map Tile Service and am using this as a guide. This has been my introduction to GIS programming and I have only been working on it for a few days so have little to no knowledge of the domain.


title: SpatiaLite: dumpkml with name and description

text: How does one get the optional arguments in .dumpkml working? According to the .help .dumpkml should work like this: > .dumpkml table_name geom_column kml_path [precision] [name_column] > [desc_column] I can get it working with > .dumpkml table geom table.kml But every point is just called 'name' and when I try: > .dumpkml table geom table.kml name_column=column All I do is lose the geometry (x,y=0,0) and everything still gets called 'name' rather than the value within the specified column. Any assistance appreciated, Phil.


title: Asynchronous Geoserver WPS Example

text: Anyone have a Java based example of how to implement an asynchronous Geoserver web processing service (WPS)? Based on some presentations that I found, async GeoServer WPS support was added in 2.1.3. I've searched the internet and also the GeoServer and GeoTools websites with no luck. Thanks in advance. My environment = PostgreSQL 9.2.7 + PostGIS + GeoServer 2.4.5 + JDK 1.7


title: ArcGIS basemap Url and querystring parameters

text: This is very basic question but might quite important.While creating Arcobjects addin we give Url for basemap as Following: `http://services.arcgisonline.com/ArcGIS/rest/services/World_Topo_Map/MapServer?f=lyr&v=9.2` At the end of the url there are few querystring parameters such as f=lyr and v=9.2. For this example they are **f** stand for Format as Lyr and **v** is I think is gdb version that is 9.2 . My question there will be more such querystring parameter and values that might determine datum and projection etc.So if anybody has any help regarding this then please shed some light.


title: Export "Zonal Statistics As Table" output in Modelbuilder

text: I have an extremely simple model that takes as input a shapefile for zones along with a value raster. All the model does is produce a table of my desired zonal statistic for each zone (the mean). **Next, I hope to get Modelbuilder to export those results** into a txt or other file I can easily manipulate later in Excel. I have about 48 rasters for which I would like to do this, and I can copy and paste the desired columns side-by-side into the Excel file myself, so although it will be a bit tedious I'm not worried about that. But it _would_ be nice if there were a way to automatically export the table. Is this possible in Modelbuilder? _Other background info, if helpful: I found the following exchange intriguing (Iterate through multiple Raster files and compute statistics using Zonal Statistics), and had hoped to be able to iterate my rasters and have them automatically append all in one table which I could later export myself in one easy click. I was unsuccessful in getting this to work, perhaps due to my misunderstanding of properly calling the target files, column names, etc. (I got the models described here to work perfectly well on my data... except that I could not successfully automate adding and naming new fields... my result had one field with the final raster/statistic having overwritten the others.)_ ETA: I'm using ArcGIS 10.2.1


title: clip linestring by length

text: Is there a way to clip a linestring based on distance? For example, linestring from a to b is 150 meters. I would like to get the linestring from a to a point that is is 75 meters away from a following the a-b linestring. What I have thought until now is to get the line NpointCount and start adding together small line segments until I get close to 75 meters, but still, how am I supposed to add the final line segment? I hope I am making sense here.. thanks in advance!


title: GDAL tiling DTED (dt2) dataset as gdal2tiles for 3D Maps

text: i am searching some webgl technologies to create 3d map with digital elevation data. for web tech. i have found ReadyMap SDK to create 3d map from TMS... in ReadyMap for elevation data it uses sth below to create a json data from tile which consists of tifs.               map.addElevationLayer(new ReadyMap.TMSElevationLayer({                 name: "Elevation",                 url: "http://readymap.org/readymap/tiles/1.0.0/9",                 args: "json=true",                 tmsType: "google"             }));      in TMS mechanism one tile scheme is sth like this: tms scheme : 10/294/752.jpg in google scheme : 10/294/271.jpg and for elevation data ReadyMap uses tiles as below: http://readymap.org/readymap/tiles/1.0.0/9/7/157/94.json?json=true i think in above link in zoom 7, folder name 157 and tile name 94 is both jpg and json format . jpg formats for tiles is in other storage or server.... and my question is that for creating digital elevation how can i tile my dted files like TMS scheme for creating digital elevation data? can i tile dted 2 with gdal2tiles or sth? thank you...


title: How to calculate tortuosity of a coastline?

text: I would like to calculate the **tortuosity** (or fractal dimension) of a coastline within a 1km buffer surrounding a point on that coastline. I have hundreds of points and therefore need a tool or scripts that will do the calculation on a mass scale. Can someone help me please? I am using ArcMap10


title: How to access .shp file from local computer or some other server

text: Below is the code in which we can access a feature layer from ArcGIS sample server. queryTask = new esri.arcgis.gmaps.QueryTask("http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/3"); But I want to access my own .shp file residing on my local computer or on another server, e.g this.. queryTask = new esri.arcgis.gmaps.QueryTask("http://172.20.100.12/Demo/Resources/ABC.shp"); In this case, queryTask does not return featureSet and so .shp file could not display on Map. Can anybody help me .. where is problem and how i can fix it.


title: ArcGIS Server 10.1 publish geoprocessing tools

text: I have been trying to publish a geoprocessing tool to our new ArcGIS for Server 10.1 with SP1. I configured the serves to not to copy data to server when publishing services. The tool ran well and was shared as a service definition file successfully from ArcCatalog. It however always failed in publish server from ArcGIS Server Web Manager due to the following error:               ERROR 001488: The service you are trying to publish includes data but the site has been configured such that data cannot be copied to the server. ERROR 001364: Validation of service definition failed. Failed to execute (Publish Service Definition). Failed.      I have been kind of fed up with the cumbersome and doggy way of publishing map/geoproc services introduced in the new ArcGIS Server 10.1. Version 10 has a much easier and clearer way of doing it. Can anyone help with any hints? Thanks Alex


title: Why is my polyline on the wrong bearing?

text: I am trying to overlay duplicate segments of polylines on a map via javascript. I am using latlon.js and geo.js from here The idea is that a user can input the distance they have travelled and when the map is refreshed, a second multipolyline will be added. To calculate the the point which represents where the user currently is, I am using a known previous location and known 'next location' to calculate a bearing between the two and then using that bearing, combined with the distance the user has travelled to calculate where they are now, before finally creating a new polyline based on that data. The problem is that this new polyline is drawn incorrectly on my map and I don't understand why. You can see the problem at [removed]. Log in as User2/pass2 and then select 'challenges' from the nav bar. The most relevant code from map.js follows. I am using leaflet maps.                 var userDistance = GymloopChallenge.user_distance;       var lastPoint;       var comparisonArr = []; //used to calculate bearing       if (userDistance <= dist1) {           comparisonArr = [rioLL,laLL];           distToAdd = Math.abs(userDistance);           userPoly = [];           lastPoint = rioLL;           console.log('if 1');            } else if(userDistance >= dist1 && userDistance <= dist2 ) {           comparisonArr = [laLL,nyLL];           distToAdd = Math.abs(dist2 - userDistance);           userPoly = [line5];           lastPoint = laLL;           console.log('if 2');       } else if (userDistance > dist2 && userDistance <= dist3) {           comparisonArr = [nyLL,lonLL];           distToAdd = Math.abs(userDistance - dist3);           userPoly = [line5, line1];           lastPoint = nyLL;           console.log('if 3');       } else if (userDistance > dist3 && userDistance <= dist4) {           comparisonArr = [lonLL,sydLL];           distToAdd = Math.abs(userDistance - dist4);           userPoly = [line5,line1,line2];           lastPoint = lonLL;           console.log('if 4');       } else if (userDistance > dist4 && userDistance <= dist5) {           comparisonArr = [sydLL,tokLL];           distToAdd = Math.abs(userDistance - dist5);           userPoly = [line5,line1,line2,line3];           lastPoint = sydLL;           console.log('if 5');       }          console.log("Distance to add: "+ distToAdd)     //console.log("COmparison array: "+comparisonArr);       var bearing = comparisonArr[0].finalBearingTo(comparisonArr[1]);       var userLocation = lastPoint.destinationPoint(bearing, distToAdd);       //console.log(userLocation);            userLineStartPoint = new L.LatLng(lastPoint._lat,lastPoint._lon );       userLocationMapPoint = new L.LatLng(userLocation._lat,userLocation._lon );       userPoly.push(new Array(userLineStartPoint ,userLocationMapPoint));       console.log(userPoly);         var userPolyline = new L.MultiPolyline(userPoly, {           color: 'red',           id:'poly2',           weight: 3,           opacity: 0.8            });       map.addLayer(userPolyline);


title: How can i calculate average slope and length of slope of a multiple field farm in Arcmap 10?

text: I always calculate slope on farms i review with 2' contour lines. Some fields are difficult since i need the average slope, but the field is rolling hills. I usally take several slopes and take the average. I want to have a way to do this in ArcMap 10. Most farms have 6 - 10 fields and every field needs its individual slope (in degree) and slope length.


title: How to start creating a map and exporting it as a raster MBTiles file

text: I am new to this whole GIS thing (my background is software-development, but not GIS) but I am super interested in making maps. So I started playing around with TileMill and was making my first steps there. Now I know the basics of where to get data from, CartoCSS, Elevation Data, GeoTiffs and so on... My goal: A outdoor map as an MBTiles file with raster data in it. The data should cover parts of Europe. Alright... I guess I should start by asking the following questions: **Tools:** Which tool would you suggest for that job? Would TileMill work? Will I have problems with TileMill2 if I want to render raster data? (Would this tool help?) I need to add GeoTiff overlays as well. TileMill2 wont be able to do that right? I also don't want to add the raster data on top of the rest, but as a layer in the middle. I am on a Mac, so Tools that would run on my machine would be best :) But if there is like _THE_ tool for it that only runs on Windows or Linux I should be able to get my hands on it as well... Oh, and inexpensive/free tools are preferred **Data:** Where do I best get my data from? When playing in TileMill1 I downloaded data from geofabrik. Now it seems as if the data isn't complete in there. Especially information about mountains... where there is forest and so on is pretty limited it seems... So I went over to TileMill2 and started playing around there. The "Mapbox Streets V4" layer has quite some information in there... But then I get very little information about f.e. where buildings are. Are there any other remote sources I can add to TileMill2? Some that are more suitable for an outdoor map? **Style my map:** In MapQuest when you zoom out you get this nice world map with sand in Africa, green in Europe, Ice at the poles.... how do they draw that? Is that the Humanitarian Data Model? I already learned how to get elevation data and how to make it visible with gdal. Is this the way to go? But then I would need a tool (like TileMill1) where I can add this raster layer. TileMill2 wont support that, right? So many questions! Thanks in advance for any tips!


title: Importing multiple GeoTiffs in TileMill

text: I am a software engineer that is very new in the geo world. So excuse my ignorance if this question seems simple or if terms I use are totally wrong. **Problem:** I have been tasked with taking multiple geotiffs and serving them from a map tile system on our servers. From my research so far, I believe that I need to convert these geotiffs to mbtiles. I've got the serving of the mbtiles part figured out, what I really need help on is coming up with a quick way to convert the geotiffs into mbtiles. **Things I've tried:** I currently have a script that can do this conversion of 56GB of Geotiffs in about 7.25 hours using gdalbuildvrt, gdalwarp, and gdal2tiles. Here is what it does:               ls -1 *.tiff > files.txt     gdalbuildvrt -input_file_list files.txt project.vrt     gdalwarp -t_srs EPSG:3857 -of VRT project.vrt project_mercator.vrt     gdal2tiles.py project_mercator.vrt project.mbtiles      It works but I'd like to improve on this. I've heard that it might be faster to use TileMill to export, but so far my attempts to import the geotiffs into TileMill have failed. I have 158 geotiffs and I don't want to have to "Add Layer" for each one using the interface. I thought I might be able to import the .vrt file I generated in my own script but when I try to add it in TileMill I get an error that says, "GDALWarpOptions.Validate() hSrsDS is not set". **Question to the community:** What is the best way to import multiple geotiffs into a TileMill project?


title: How to get a set of all buildings of an OpenStreetMap section in form of single points / Multiple polygons as one point

text: The objective is to get a set of all buildings of an Open Street Map section in form of single points in a shapfile.   1. I downloaded the section from osm.    2. Imported the section into QGIS (and later imported it into a PostGis db using osm2pgsql)   3. To transform the polygons (representing buildings) into points I used the Polygon Centroids function Unfortunately some of the buildings are represented by multiple polygons subsequenty I did get multiple points for one building which I is not what I wanted. Searching the web I found out that these polygons should be multipolygons but neither qgis nor osm2pgsql does import these relations so that I could merge the poygons of a building Any ideas how i could get single point list of the buildings?


title: Howto create complex polygon with pyshp

text: using the pyshp library (https://code.google.com/p/pyshp/) I am trying to create a "complex" polgon with one hole inside the main-poly plus one island- poly. It should look like this: ![enter image description here](http://i.stack.imgur.com/fgahs.jpg) I used the following code from the docu, only changing the number of points and using anticlockwise rotation for the hole:               import shapefile     w = shapefile.Writer(shapefile.POLYGON)     w.poly(parts=[[[0,50],[50,50],[50,0],[0,0],[0,50],[10,40],[10,10],[30,10],[30,40],[10,40],[70,20],[100,20],[100,0],[70,0],[70,20]]])     w.field('FIRST_FLD','C','40')     w.field('SECOND_FLD','C','40')     w.record('Poly','PolyTest')     w.save('TESTPOLY')      I get this: ![enter image description here](http://i.stack.imgur.com/JRQbD.jpg) Obviously I am making a stupid python-rookie-error. Any hints are very appreciated. Thanks.


title: Why can't I save a line shapefile?

text: I'm working with qgis v. 2.0.1. I have some data points and I just wanted to create a new shapefile with line geometry. In the new vectorlayer-window, I added some attributes (approx. 10) and klicked ok but it didn't create a shapefile (nothing happened). I tried again without adding any attributes and it worked. Why is there a problem when I want to add some attributes? Of course I could add attributes with the Table Manager or manually in the attribute table but I'm just wondering why it's not working this way. Thanks.


title: Failed to create signature file

text: I need some help, because I've no clue how solve this issue the exercise that I'm trying to perform. The point where I'm stuck is 4. "Create signature file as LandsatCompS.sgs". I've done what it says in the exercise. I got a error message "Failed to create signature file. The aim of the exercise is to carry out a supervised image classification of the Northern half of Jutland. 1\. Start ArcMap and add the Landsat images from the B series (summer images).   1. Create a composite image composite image for all bands except band 6 – the thermal band. Save the image as LandsatCompS (S for summer).   2. Open the orthophotos for land (not urban). Add eventually more orthophotos from the Geodata library. Try to identify the necessary number of different land cover classes to cover all areas on the orthophotos. Not more than 15!   3. Open the Image Classification window, and start digitising training areas. Two for each class. **Create the signature file as LandsatCompS.gsg.**   4. Carry out maximum likelihood classification using LandsatCompS as input raster and LandsatCompS.gsg as signature file. Output is ML_class_NJ1. Use default values for the other inputs.   5. Put the right colours on the new map and compare with the orthophoto and other maps available for Northern Jutland.   6. Where are the problems? Try to make an improved classification by increase the number of classes and digitise new training areas. Carry out the classification and assess the result. Which one is best? ![Classes](http://i.stack.imgur.com/uj7N0.jpg) ![Error message](http://i.stack.imgur.com/OWrS6.jpg) ![Landsat Image](http://i.stack.imgur.com/OQfRO.jpg) ![Training samples](http://i.stack.imgur.com/mNU3E.jpg) ![enter image description here](http://i.stack.imgur.com/wNgdZ.jpg)


title: How to properly assign spatial reference to converted AutoCAD data?

text: I know this is a frequent problem, I myself confronted with this problem many times, yet somehow each Autocad dataset I received is different (although they're from the same area). To some data, the projection I gave them worked perfectly, and in this case it didn't. So, I've recently received some AutoCad data. Converted them in ArcMap. They are missing spatial reference. But I know their location. They're suppose to be in a specific city in Romania. So first thing I did, I gave them the Stereo 70 projection (Romania's national projection), I assumed this is the projection they were created in, the data gets somewhere in the middle of the Mediterranean Sea. Second thing I tried is is giving them WGS84, and they got somewhere in the Atlantic Ocean ;) PS: I use ArcGIS 10.1 PS2: I've tried all 3 conversion methods (I know) in ArcMap. PS3: The autocad data came with no spatial reference. I only know they should be in a specific city.


title: 404 Error with Openlayers WMS layer

text: I'm trying to create a map with the combination of openlayers and extjs. The map works fine and everything seems to be fine except for a a few 404 errors which I get in the browser console. Although this doesn't seem to affect any functionality, I had to get rid of this. The code and the error is as below: Layer creation:               var wmsLayer = new OpenLayers.Layer.WMS("vmap0",                               "label",                               {'displayInLayerSwitcher':false                               });      Error:               "NetworkError: 404 Not Found - http://localhost:8080/myproject/label?DISPLAYINLAYERSWITCHER=false&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&FORMAT=image%2Fjpeg&SRS=EPSG%3A4326&BBOX=-180,-16.857142857143,-106.85714285714,56.285714285714&WIDTH=256&HEIGHT=256"      I get 15 error of same type except that the value in the `BBOX` parameter alone changes. Please help to get rid of this errors. Code to draw layers:               var vecLayer = new OpenLayers.Layer.Vector("vector", {                           rendererOptions: { zIndexing: true },                           styleMap: new OpenLayers.StyleMap({                               'default': allCountryStyle,                               "select": new OpenLayers.Style({                                   strokeWidth: 3,                                   strokeColor: "#000000",                                   graphicZIndex: 2                               })                           }                           ),                           protocol: new OpenLayers.Protocol.HTTP({                               url: url,                               format: new OpenLayers.Format.GeoJSON()                           }),                           strategies: [new OpenLayers.Strategy.Fixed()]                       });


title: Accessing GeoServer on my own domain

text: I am running GeoServer 2.4.3 on Ubuntu Server 12.04 and tomcat7 and I want to access it on mydomain.com. I have searched a lot and found the all the solutions suggest to use a ProxyPass. One of the posts on GIS Stack Exchange that has a similar problem is this one but the website to the solution provided in the post is down! This is what I have been doing so far: Added the follwing to `/etc/apache2/sites-available/default`                ServerName mydomain.com      ProxyRequests Off      ProxyPreserveHost On           <Proxy *>          Order deny,allow           Allow from all       </Proxy>      ProxyPass /geoserver http://servername:8080/geoserver      ProxyPassReverse /geoserver http://servername:8080/geoserver      Then I have edited `/var/lib/tomcat7/webapps/geoserver/WEB-INF/web.xml`               <context-param>           <param-name>PROXY_BASE_URL</param-name>           <param-value>http://mydomain.com:8080/geoserver/</param-value>         </context-param>      Finally, I have added In GeoServer > Global Settings > Proxy Base URL > http://mydomain.com:8080/geoserver/ Now, when I am trying to enter mydomain.com:8080/geoserver I get Unable to connect error! What I have been doing wrong? Any ideas?


title: How to use Advanced Symbology in ArcGIS for Android?

text: I try to use the advanced symbology in order to draw nato symbols in JAVA I need little help about it,   1. What the parameters _WKID mean like in this command : lmessage.setProperty("_WKID", "3857");   2. What the parameter UniqueDesignation mean?   3. THE MOST IMPORTENT QUESTION: How can send parameters to symboles llike the symbole "Target" (TACGRP.FSUPP.PNT.TGT.PTGT) (2.X.4.1.1.1) with sic G*FPPTS--- _*_ *X parameters like "Target Alt" , "Description" and "Name"?


title: How to access the button on the form?

text: I have a form that appears when you start the plugin.How do I properly apply to the button on the form, so that when you click on the button performed "function".               import sys,os     import win32api     from PyQt4.QtCore import *     from PyQt4.QtGui import *     from qgis.core import *     import ctypes     from PyQt4 import uic          class HelloWorld:        def __init__(self, iface):             self.iface = iface           self.canvas = iface.mapCanvas()             def initGui(self):             pushButton=QPushButton("pushButton")           self.action = QAction(QIcon(":/plugins/"), "&HelloWorld", self.iface.mainWindow())           self.action.setWhatsThis("Configuration for test plugin")           self.action.setStatusTip("This is status tip")           QObject.connect(self.action, SIGNAL("activated()"), self.hello_world)                 QObject.connect(pushButton, SIGNAL("clicked()"),self.function)                self.iface.addPluginToMenu("HelloWorld", self.action)                def unload(self):               self.iface.removePluginMenu("HelloWorld",self.action)           self.iface.removeToolBarIcon(self.action)             def function(self):           QMessageBox.information(self.iface.mainWindow(), QCoreApplication.translate('HelloWorld', "HelloWorld"), QCoreApplication.translate('HelloWorld', "HelloWorld"))           return             def hello_world(self):                 path = os.path.dirname( os.path.abspath( __file__ ) )           self.optionsDialog = uic.loadUi(os.path.join(path,"untitled.ui"))           self.optionsDialog.show()          if __name__ == "__main__":        pass


title: Where does Geo-Platform fit in a GeoStack?

text: I am totally new to GIS development and I am trying to grasp the concept. From my understanding the stack looks as follows (Please correct me): Back End GIS Server: GeoServer/GeoNode/ArcGIS Server (Some of the back end servers do have their own spatial database - however, PostGIS can still be used in conjunction with them) Back End Spatial Database (optional): PostgreSQL/PostGIS database can be used in conjunction with back end GIS server (ArcGIS for example) Front End Library (to display the shapes/points on top of a map(Google Maps/Open Street/etc.)): Openlayers, lefletjs and GeoExt Now I am reading about Geo-Platform and I am not sure where exactly does it fit in the above?


title: Field-based hyperlinks in an ArcMap map package (.mpk)?

text: I have a point shapefile where each point has multiple hyperlink fields in its attribute table as per this useful question and answer. These hyperlinks are to photos stored locally. Ideally I want to bundle up the map (shapefiles and all) and the photos and export them as a map package (*.mpk) so that the receiver can open it up and click on the hyperlinks in the shapefiles. Is there any way of doing this? I know the map package export function has the option to include additional files but I am not sure how to get this to work for me. Attachments are not an option, as my licence is Basic, and geodatabases or dynamic hyperlinks are not options as I need each point feature to link to multiple photos. Thanks!


title: Examples of Animated maps

text: I am currently looking examples of animated maps. Creating animated maps with the help of this post and similar post checking Animaps site also I am presenting this map in local school so those students should take interest to learn and create various simple interactive maps. I am in process to create animated map contains Gardens and amusement parks across the city. So any help will be great !!! Thanks in advance : )


title: Loading data onto PostGIS using KML files with attibute data

text: So, I have a kml file which contains 317 data points (317 polygons and the respective 317 names). The kml file also contains attribute data which pops up in a bubble on Google earth with about 12 custom data points (Code, name, locality, no. of houses etc.) which could have been loaded onto the kml file using ArcGIS (I'm assuming this, cause I received just the kml file from my source) or could have been custom coded using a text editor. Now, I'm trying to build a PostGIS db using this data along with the custom data points (attribute data) which is stored in the kml file. I tried converting the kml file to a shapefile to load it onto PostGIS using shp2pgsql, but this ignores all the attribute data. Only the geometry and name gets converted. I also tried converting the KML file to xml/csv/xls format to atleast try and merge the attribute data with the main db on postgresql, but I failed cause there is no converter out there which pulls out all the data from a KML file. My question is, is there a way to load this kml file with attribute data onto PostGIS? Is there a simpler converter for KML to CSV, or similar spreadsheet formats? What's the best way to do this?


title: Vertical Datum Conversion Utilities

text: I have some data (points and rasters) that I would like to be able to convert to and from several vertical datums, including MLLW in coastal Florida. I would appreciate it if people could provide some insight into accuracy, pitfalls, requirements, etc. of using available utilities such as VDatum or corpscon. Is one more user-friendly than the other? Are there any others that I should consider? Thanks!


title: Check Text Mode when Typing text MapInfo

text: I'm writing a code MapBasic How do I know text mode when I add text to Map. Can help me!, check cursor texting mode


title: Error 14: Could not open CRS database

text: I am using Python Console and loading a Vector Layer using python commands. When the layer is added, QGIS first asks for the CRS to use, but I don't have all the CRS listed out to select the Reference system from and the "Error Log" shows the error as: **_> "Could not open CRS database C:/Program Files/QGISDufour/./resources/srs.db. Error(14): unable to open database file"_** Now what shall I do get this cleared. I cannot load my vector layer because of CRS problem.


title: Openlayers create, edit and save/export to kml

text: I'm trying to create, edit and save a point vector as kml file in an Openlayers app. I have managed the create and edit by adding and empty vector layer and then adding the control OpenLayers.Control.EditingToolbar("name of the layer") But how could I save/export this afterwards in kml format? I just need this on the browser side a popup to let me save the file locally not on the server side. Thanks Aris


title: ArcGIS Silverlight Web Mapping Application has two logins

text: I am trying to implement a Silverlight Web Mapping application with ArcGIS server and I'm trying to use the Token authentication with active directory as the user/roles store. I have the service running so that I can successfully obtain tokens but when I try to initialize a map service the ArcGIS server challenges me again to authenticate with widows even though I have successfully created a token. The IIS is set for windows authentication but in any other configuration (Anonymous, Impersonate) the map service never loads.                       //Import  Base Data Map Service             ArcGISTiledMapServiceLayer BaseMap = new ArcGISTiledMapServiceLayer();             BaseMap.Url = DataSourceStrings.ProductionServer + DataSourceStrings.I_GeneralWesternCanadaURL;             BaseMap.Token = Token; // This is a valid token             BaseMap.ID = "BaseMap";             BaseMap.Visible = true;             //BaseMap.ImageFormat = ArcGISDynamicMapServiceLayer.RestImageFormat.PNG32;             BaseMap.Refresh();             BaseMap.InitializationFailed += ArcGISDynamicMapServiceLayer_InitializationFailed;             BaseMap.Initialized += ArcGISDynamicMapServiceLayer_Initialized;             Map.Layers.Insert(0, BaseMap);      Solution...? So I managed to fix this by setting the Server.dat file to disable security but the services still require a token in their web.config files. Seems Arcgis server like to handle the windows authentication but the design was to handle the authentication in silverlight using a wcf service. I've tested this and everything seems to work the way i want it too. Kinda of a hack way of doing it but everything works right.


title: raster data retrieval from postgis using c++

text: I have raster data in a postgis server (I am able to see it in qgis). When I retrieve it and form gtiff image out of it using gdal, the generated image is not correct (there is a coloring problem). I am able to retrieve the metadata correctly. For raster I am doing `st_AsBinary(rast)`. And then passing it to gdal raster io buffer parameter Metadata is read first successfully. Image was stored as tiles **Code:**               for(int row=0;row<99;row++){     int xoff,yoff;             yoff=row*metadata.tilesizeY ;             for(int col=0;col<99;col++){                   xoff=col*metadata.tilesizeX;                   query<<"select st_asbinary(rast) from public.tkptest where rid ="<<(row*99)+col+1;                   retval = PQexec(options.connection,query.str().c_str());                   res= PQresultStatus(retval);                   if(res!=PGRES_TUPLES_OK){                         std::cout<<"We did not get any data\n";                         return res;                   }                   dstDS->RasterIO(GF_Write,xoff,yoff,metadata.tilesizeX ,metadata.tilesizeY ,PQgetvalue(retval,0,0),metadata.tilesizeX *metadata.datasize,metadata.tilesizeY *metadata.datasize  ,metadata.dtype,metadata.numBands,panBandMap,0,0,0);                   PQclear(retval);                   query.clear();                   query.str("");                  }     }     GDALClose(dstDS);


title: Features to 3D from 2D in QGIS

text: Is there a way to convert 2D features to "3D" features in QGIS? I know there is a way to do this in ArcMap by using an attribute value in the attribute table, but I am looking to do this with shapefiles in QGIS. Thanks...


title: How to compute the bounding box of multiple layers in lat/long?

text: I am writing an application to performance test all kinds of map services, primarily AGS 9.x, AGS 10, and WMS 1.x. Part of the application involves generating random bounding boxes for individual requests, within the service's full extent. This part is working fine for geographic and projected coordinate systems when the service's full extent is known (e.g. through an AGS service's fullExtent property). My problem is with WMS: each layer in a GetCapabilities response can define its bounding area in >= 1 CRSs. Some parts of the application need to know if a service's CRS is geographic or projected, so in order to remove ambiguity in WMS I am always using the layer's LatLonBoundingBox which is always defined and in EPSG:4326. I then have to calculate a full service bounding box based on all layers that go into an individual request (which are randomised). This is where it gets tricky. I am getting lost because for each lat / lon bounding box the LLx (lower left longitude) could be a larger or smaller number than the URx (upper right longitude), depending on which meridians it spans. Every time I start drawing square or circular diagrams I think I have an approach figured out, and then find a case that ruins it and my brain turns to mush. I am going to keep thrashing at it until it works, and if I get a solution post it here, but I'm sure there has to be an accepted and fully tested approach which would make my life easier. I just can't find it right now.


title: GeoDjango - How to secure tile map services?

text: I use geodjango to create and serve tiles that I usually display into OpenLayers as openLayers.Layer.TMS I am worried that anybody could grab the tms URL and plug it into their own map without asking permission, and then consume a lot of the server's CPU and violate private data ownership. On the other hand, I want the tile service to be publicly available without login, but from my website only. Am I right to think that such violation is possible? If yes, what would be the way to be protected from it? Is it possible to hide the url?


title: Black & White / Greyscale Maps and Illustrations

text: I think greyscale maps & illustrations (with geographic data) seem to be abandoned and the colour maps are mostly used nowadays. This is understandable and personally I also need less time to produce a colour map than a grey scale map, although in the end I fell much more satisfied with the grey scale results.   * What are your thoughts?   * Is it still feasable to elaborate grey scale maps and illustrations?   * Where can I find some good tutorials on this(e.g. advises on grey color selection, how to deal with such a limited color range,etc)? I am not looking for "converters" from color to grey scale, like the one from ESRI (that seems to work ok) that can easily and with soem quality export a colour map to grey scale: http://mappingcenter.esri.com/index.cfm?fa=ask.answers&q=2163


title: Calculate min distance between points in PostGIS

text: I have point layer in PostGIS database. Help please with query which add for each points in layer distance in attribute table between this point and another nearest point. Thanks! **Update**   Calculate distances between series of points in postgis there is similar question, but there is in answer:               UPDATE my_table SET dist=ST_Distance(my_table.geom, b.geom)      FROM (SELECT geom FROM my_table WHERE **gid = 1**) b;      distance added for points - from **gid = 1** point to another points. For me each point object in attribute table should have distance from point to another nearest point. **Update** ![enter image description here](http://i.stack.imgur.com/nACZY.jpg) For example, I need add distance between point and another **nearest** point in attribute table for point: for point 1 - distance between 1 and 2; for point 2 - distance between 2 and 3; for point 3 - distance between 3 and 2; for point 4 - distance between 4 and 3; for point 5 - distance between 5 and 4; for point 6 - distance between 6 and 3.


Does a geoprocessing tool exist that can find all DTED tiles that intersect an AOI? I have a area of interest represented by a polygon and DTED tiles (some of which intersect with the AOI). Is there a tool that I can use to find which tiles intersect? Thanks!


title: Using OpenLayers - seeking example for saving feature to database using Ajax

text:
I'm trying to save a point to the database using a balloon that has an input field and a save button. The button triggers the code and sends the feature info to the PHP page. It works, but the save action refreshes the page. I would like it to save in the background without refreshing the page, like Ajax. Can someone advise how to do this ?                  var mySaveStrategy = new OpenLayers.Strategy.Save({ auto: true });              var myProtocolHTTP = new OpenLayers.Protocol.HTTP({                  url: "proc/newPoint.php",                  params: {                          m: 'r',                          srs: 'EPSG:4326'                  },                  callback: function() { alert("User callback called!"); },                  format: new OpenLayers.Format.GeoJSON({                          ignoreExtraDims: true,                          internalProjection: map.baseLayer.projection                 })          });       I edited the code to reflect my latest code (16.12.2012):                       function handler(request) {                 alert(request.responseXML);                             }                                 var myProtocolHTTP = OpenLayers.Request.POST({                 url: "api.php?action=insertPOI",                 data: "",                 params: {                     m: 'r',                     srs: 'EPSG:4326'                 },                 headers: {                     "Content-Type": "text/plain"                 },                 callback: handler             });                   function formAction(action) {             var feature = mytoolbar.selCtrl.openpopup.feature;             switch (action) {                 case 'delete':                     feature.state = OpenLayers.State.DELETE;                     vectorLayer.strategies[0].save();                     break;                 case 'cancel':                     break;                 case 'save':                     var attrs = feature.attributes;                     attrs.name = $('fname').value;                     attrs.name_lat = $('fname_alt').value;                     attrs.type = $('ftype').value;                     attrs.notes = $('fnotes').value;                          if (feature.mystate) {                         feature.mystate = null;                         feature.state = OpenLayers.State.INSERT;                         var ctrl = mytoolbar.controls[1];                         ctrl.layer.addFeatures([feature]);                         ctrl.featureAdded(feature);                         ctrl.events.triggerEvent("featureadded", {'feature': feature});                     } else {                         feature.state = OpenLayers.State.UPDATE;                         vectorLayer.strategies[0].save();                     }                    break;                  }              function getPopupForm(feature, context) {                  var gid      = feature.attributes.gid      || '';             var name     = feature.attributes.name     || '';             var name_alt = feature.attributes.name_alt || '';             var notes    = feature.attributes.notes    || '';             var html = "
" +                 "<form action='' name='baloon'>" +                    "<input id='fgid' type='hidden' value='" +                       gid + "' />" +                    "<button type='submit'  id='baloon-save' value='Save' onClick=\"formAction('save')\" />"+                    "	" +                    "
"+                    "<button type='submit' id='baloon-cancel' value='Cancel' onClick=\"formAction('cancel')\" /> " +                    "

"+                    "<button type='submit' id='baloon-delete' value='Delete' onClick=\"formAction('delete')\" />
"+                    "

</form>" +                "
";                  return html;         }      '


title: Cannot intercept Alt button when user drags a box (OpenLayers)

text: I am using OpenLayers. When a user holds Shift, clicks on map and moves the mouse, the box is shown. When user releases the mouse button, I intercept the event and use the coordinates of the box.                       var boxControl = new OpenLayers.Control();             OpenLayers.Util.extend(boxControl, {                 draw: function () {                     // this Handler.Box will intercept the shift-mousedown                     // before Control.MouseDefault gets to see it                     this.box = new OpenLayers.Handler.Box(boxControl,                             { "done": this.notice },                         { keyMask: OpenLayers.Handler.MOD_SHIFT });                     this.box.activate();                 },                      notice: function (bounds) {                          var ll = map.getLonLatFromPixel(new OpenLayers.Pixel(bounds.left, bounds.bottom));                     var ur = map.getLonLatFromPixel(new OpenLayers.Pixel(bounds.right, bounds.top));                     doSomething(ll, ur);                 }             });             map.addControl(boxControl);      It works when with Shift button. But when I change MOD_SHIFT to MOD_ALT, my code is ignored, the bounding box is never shown on screen and "notice" function is never called. Could you help to fix it (for either Alt or Control buttons), please? PS. I cannot use OpenLayers.Control.Select


title: problems with r.sun from GRASS

text: I´ve problems with r.sun from GRASS. I need to calculate beam radiation, reflect radiation and diffuse radiation wwith this tool. my problem is I do not understand that you have to put mandatory files as "imput" because I put MDT, aspect layer and slope layer, number of cells fromm MDT, ando MDT coordenates. when I put these file the result is: > Oooops! SEXTANTE could not open the following output layers Output > irradiation layer [Wh.m-2.day-1]: > C:\Users\Leticia\AppData\Local\Temp\sextante\irrad_layer.tif Output > insolation time layer [h]: > C:\Users\Leticia\AppData\Local\Temp\sextante\insol_time.tif Outpu diffuse > irradiation layer [Wh.m-2.day-1]: > C:\Users\Leticia\AppData\Local\Temp\sextante\diff_rad.tif Output ground > reflected irradiation layer [Wh.m-2.day-1]: > C:\Users\Leticia\AppData\Local\Temp\sextante\reflec_rad.tif Output global > (total) irradiance/irradiation layer [Wh.m-2.day-1]: > C:\Users\Leticia\AppData\Local\Temp\sextante\global_irrad.tif The above > files could not be opened, which probably indicates that they were not > correctly produced by the executed algorithm Checking the log information > might help you see why those layers were not created as expected This > algorithm requires GRASS to be run. A test to check if GRASS is correctly > installed and configured in your system has been performed, with the > following result: GRASS seems to be correctly installed and configured" Can anyone help me? Please!!!!!!


title: How to display labels from many-to-many relationship in ArcGIS 10 using VB.NET?

text: I have two tables in a geodatabase:   * parcel; and    * person  1 or + parcel, can have 1 or + person (relationship is many-to-many) Using VB.NET I can display all the parcel labelled with the "id" field (see picture on left). What I want is to display the parcel with labels for the Names of their persons (see picture on right): ![enter image description here](http://i.stack.imgur.com/yRYjj.png) I think that I have to use join between two tables I can get the name with this code                Sub addNom()             ' Open the origin and destination classes.             Dim featureWorkspace As IFeatureWorkspace = CType(Sp, IFeatureWorkspace)             Dim originClass As IFeatureClass = featureWorkspace.OpenFeatureClass("parcelle")             Dim destinationClass As IObjectClass = CType(featureWorkspace.OpenTable("personne"), IObjectClass)                  ' Create a memory relationship class factory and open a relationship class with it.             Dim relClassFactoryType As Type = Type.GetTypeFromProgID("esriGeodatabase.MemoryRelationshipClassFactory")             Dim memRelClassFactory As IMemoryRelationshipClassFactory = CType(Activator.CreateInstance(relClassFactoryType), IMemoryRelationshipClassFactory)             Dim memRelClass As IRelationshipClass = memRelClassFactory.Open("pll_personne", originClass, "id", _                                                     destinationClass, "id", "Contains", "Is in", esriRelCardinality.esriRelCardinalityOneToMany)                  ' Find the counties related to the state with an ObjectID of 1.             Dim originFeature As IFeature = originClass.GetFeature(6)             Dim relatedSet As ISet = memRelClass.GetObjectsRelatedToObject(originFeature)             relatedSet.Reset()                  ' Iterate through the set while finding the names of the counties.             Dim nameIndex As Integer = destinationClass.FindField("nom_fr")             Dim setObject As Object = relatedSet.Next()             While Not setObject Is Nothing                 Dim row As IRow = CType(setObject, IRow)                 Console.WriteLine("OID: {0}, NAME: {1}", row.OID, row.Value(nameIndex))                 setObject = relatedSet.Next()             End While         End Sub      The Question is: **How can I display the labels using this result?**


title: Submit empty string via POST wfs in GeoExt.form.FormPanel

text: I have a two textfields in FormPanel. When one of them is empty, after pressing the search button the form sends the empty string in WFS query. How can I disable the empty field which sends the WFS request after I pressed the search button?                   formPanel = new GeoExt.form.FormPanel({         ref: "formPanel",         width: 300,         height: 200,         region: "west",         protocol: protocol,         buttons:[{          text: "search",          handler: function() {              formPanel.getForm().search();          },          //scope: formPanel      }],         items: [{             xtype: "textfield",             name: "trans_id__eq",             value: "",             fieldLabel: "ID",             disabled: true,             ref: 'myTextField'         },              {             xtype: "textfield",             name: "nabywca__eq",             value: "",             fieldLabel: "nabywca",             submitEmptyText: false          },         ],         listeners: {             actioncomplete: function(form, action) {             // this listener triggers when the search request             // is complete, the OpenLayers.Protocol.Response             // resulting from the request is available             // through "action.response"                 features = action.response.features;                 app.featureGrid.store.loadData(features);   //załadowanie zwróconego wyniku wyszukiwania do "app.featureGrid.store"                 vm=app.mapPanel.map.getLayersByName("Results");                 if(vm.length==0){                     vecLayer = new OpenLayers.Layer.Vector("Results");                    app.mapPanel.map.addLayer(vecLayer);                     app.featureGrid.store.bind(vecLayer);                     app.featureGrid.getSelectionModel().bind(vecLayer);                 //app.featureGrid.getSelectionModel().bind(vectorLayer);             }         }         }                         })                   );      When I leave blank first textfield and put the value to the second my wfs request send empty string and I get the error. My protocol:               var protocol= OpenLayers.Protocol.WFS({                      version: "1.1.0",                 url: "http://localhost:8080/geoserver/wfs",                 featureType: "transakcje_punkty",                 featureNS: "http://geoserver.org/transakcje",                  srsName: "EPSG:900913"                       });      Wfs is working properly when I have one field, but I would like to add 7 "textfields". The user fills the fields with the attributes he want to find and rest leaves empty. But empty field send some empty value to WFS and I get an error.... instead results


title: How to solve routes to offset points in Network Analyst?

text: The scenario is that there are points on a grid with intervals of 1000m covering an area of hundreds of square kilometres. I want to find routes on a road network to any one of these points. A road network built with Network Analyst works fine and the points are also built into the network. However the points do not intersect with the roads and sometimes offset 10s to 1000s of metres so no solution/route can be found to these points. How can routes be solved to the closest part of a road nearest to these points if they are loaded as stops?


title: Osgeo Welcome window

text: Is it possible to disable OsGeo Live 7.0 Welcome window? It annoys me a bit appearing each time I Run the system. Solved by: menu->settings->settings manager->session and startup->"Application Autostart"->uncheck the message appearance


title: How to get features in Openlayers vector layer

text: What is the exact API to get all features in a vector layer in OpenLayers? Suppose I have the "map" object, within which is a "Vector Layer": how can I obtain features in the layer?


title: Why does ArcMap 10 export circles instead of diamonds?

text: My default symbol for building is a small square. A jpg export or any other export renders it as a small circle. ## It's a killer for hydrology maps which use circles for many other items. Solutions? From the posts i went back and looked at the default symbols that the open- source layer provided. It LOOKED like a square rotated. But it was something else...yes a circle! Not sure why this happens in Arc...However I simply changed my symbol and voila! it exported just fine. While i got what i needed the probelm persists. When i went back and looked at the square it said it was a polygon! Odd. Anyway thanks to these posts i now have a work around; many thanks ![compostite of parameters for default symbol](http://i.stack.imgur.com/Te9s5.gif)


title: What does "Input geometry has unknown (-1) SRID" mean?

text: I am trying to run the following query to determine the area and number of points of a polygon:               Select osm_id, st_npoints(the_geom) as numPoints, st_area(transform(the_geom, 3395))/1000000 as areakm2 from water_table where st_isvalid(the_geom);      I receive the following error:               Input geometry has unknown (-1) SRID      Anyone know what the problem is?


title: arcgis server 10 tif mxd problem

text: i am using arcgis server 10. first i opened arcmap and added a.tif file. then i saved it using                File ->  Save As       as a ArcMap Document (*.mxd) file, let's say b.mxd. Now i open arccatalog and added arcgis server let's say "test". then by right click on this server i added new service (my b.mxd). when i press finish, it says the following:                Server Object instance on machine arcgis-server is created with errors.         The Layer: 'a.tif' in Map: 'Layers' is invalid.      Failed to open raster dataset      Can you plz explain what am i doing wrong? thanks in advance!!!


title: How to add an identify button with Arcgis JS API?

text: GIS stackexhange is the first website that i've come across that has come near to answering my question, but some how i can't work the answer out in my head/or the answers here have confused things more. I'd like to put a simple identify button on a web map. In simple terms, i want a map where the identify task is disconnected onload with an equally simple toggle on/off identify button. i understand that you need to dojo.disconnect the onclick identify function and that you then need to dojo.connect it back in the onclick of a button but i cant work out how it all comes together. I've read and re-read the arcgis js api 'working with events' page but this too has lost me. any help would be great.


title: Append and Merge tools issues in Arc 10

text: Append and Merge tools issues in Arc 10. keep getting error 000278 - failed on input OID ...... (could not write value to one field or the other). Does anyone know how to fix this? Never had this problem in Arc 9.3, where if the value did not fit, it would get chopped off.


title: GPS - Garmin 60CSX connectivity issue

text: I have a Garmin 60csx, about 2 years old and I have attempted to connect it to my laptop/ArcView (GPS toolbar) using the Detect GPS Port to no avail. There is no port number being displayed and that is what it states in the popup box. the baud rate, parity etc... are accurate. Should the Comm Port have a number displayed upon connection if it is working? Is there a patch or a specific usb cable to use or specific Garmins will only function properly in ArcVIEW? * * * now it is saying the USB is not recognized. I'm thinking a new USB cable as one solution.


title: Where deploy Custom Provider DLL for ArcGIS Server if GAC is not an option

text: I have a client that has developed a custom Membership Provider for ArcGIS server and cannot deploy the provider to the GAC. Does anyone know where we could put the DLL such that the ArcGIS Security and Token applications can "see" it? Thanks, Carl


title: Google Maps 3d with geoserver wms

text: I am using opengeo suite with openlayers and google maps as my base layer. I have points in a postgis database that I am using to generate wms tiles. I'm pretty sure I have the projections right because the points are spot-on for the most part. I run into problems when I get in to the last two zoom levels when the google maps turn into that 3d view in some areas. The markers then show up as in the wrong places and if I pan the map they pan at a different pace. If I don't let the map zoom in that far, its fine, but we have to be able to zoom in closer. Can I turn off the 3d view on google maps or adjust the projection?


title: Detect city from coordinates

text: I need to know how to detect which city belong certain coordinates. Wich is the best way to do this? Thanks in advance. Pd: excuse my english


title: Struggling with ArcGIS REST API and tokens

text: I'm doing my best to secure my AGS (9.3.1) services, but have struggled with tokens. As you can see, my _nw_ account belongs to the _agsadmin_ and _agsusers_ roles: ![enter image description here](http://i.stack.imgur.com/t9DSq.png) I have assigned these same roles to my service and its containing folder: ![enter image description here](http://i.stack.imgur.com/GIfmp.png) Generating a token is no problem: ![enter image description here](http://i.stack.imgur.com/alJ6B.png) But alas, upon logging in to the REST interface, I am presented with this _Unauthorized access_ message: ![enter image description here](http://i.stack.imgur.com/WRX60.png) I found a post on the ESRI forums describing the same problem, but responses are inconclusive. ("Please make sure that the Security settings for this map service is correct.", blah blah.) Has anyone encountered this? What could I be doing wrong? **Update: we have abandoned AGS-based security because of its many flaws and capricious inconsistencies. Beware!**


title: labeling in qgis 2.0: if one letter in string is "X"

text: I have a shapefile and I want to label all shapes with a "B" or "A" in its name (string). All others should be ignored. Like: _If [field] contains "B" draw label_ Im using qgis 2.0.


title: The simplest way to uninstall any and all ArcGIS products?

text: What it is the simplest and most straightforward way to uninstall any and all ArcGIS products more or less in one go that might be present on a Windows machine? perhaps scriptable? Add/Remove Programs from the Control Panel is straightforward, but when there a several products installed it's a process of click-click-click, wait for product 1 to be uninstalled, which could take some time, then click-click- click for #2, and so on. This is okay for one machine but when there are several it's real pain.


title: Why am I losing part of my raster when I reproject?

text: QGIS 2.0 Working with tiff raster file. Reprojections are in WGS84 Mercator. Here is my problem: I am trying to use buffers or urban areas to analyze light intensity on a raster layer. I want 30 miles, 60 miles, and 100 miles for the buffer sizes. I have reprojected the vector layers to WGS84 / World Mercator. When I ran a zonal statistics, everything went smoothly until I was double- checking the output. I saw "null" on a few of the rows of information in the data table. When I looked at the raster image, which I also reprojected to WGS84 / World Mercator, a large chunk of the raster layer was missing. Images below. Non-reprojected raster layer: ![enter image description here](http://i.stack.imgur.com/YrygV.jpg) Reprojected raster layer: ![enter image description here](http://i.stack.imgur.com/PRuph.jpg) I am assuming there is a problem with my reprojections. I am still new at GIS, so please bear with me. Thanks in advance.


title: Lat/lon of pixel inside Landsat scene

text: How do I find the lat/lon of each pixel inside a landsat scene. Either the equations for converting i,j to lat/lon or just a bunch of files, one for each scene. Background: I have downloaded a bunch of scenes from Greenland and now I want to compare them to some MODIS data. MODIS provides files with lat/lon for every pixel.


title: How to write to Sql server from MapBasic and insert rows

text: I have created a table in Sql server 2014 with the columns: Id (int), Geom (geometry), RouteId (int), Desc (varchar(MAX)) I have succeeded in making the table mappable through MapInfo I have checked .TAB file for the table, that it has IsReadOnly = false Stil when i try to execute an insert query on the table I get the error messsage that the table is read only. Anyone knows how to fix this? Also I cannot figure out how to write a correct insert query I have tried doing that on a copied version of the table code is as follows               Dim rId As Integer     Dim oPline as Object     rId = 20     Dim t As String     t = "test"     Create pline into variable oPline 2 (1,2)(1,3)          Insert into TestGeom1 (Geom, RouteId, Desc) VALUES (oPline, rId, t)      Writing like above I get the error messsage field Geom does not exist in table TestGeom1 If i replace Geom with obj or just remove it completely i get the error Expression does not evaluate to a column or a table name. So please who can tell me what is the correct syntax for this?


title: How to cluster Google Maps kml overlay

text: I have several thousand markers stored in a kmz file that I load into a Google Maps (v3 API). I've looked around to find a way to cluster the overlay data, but to no avail.   MarkerClusterer looks nice, but it only seems to work with a marker array?               function initialize() {         geocoder = new google.maps.Geocoder();         var latlng = new google.maps.LatLng(25,-20);         var myOptions = {             zoom: 2,             center: latlng,             mapTypeId: google.maps.MapTypeId.ROADMAP         };         map = new google.maps.Map(document.getElementById("map_canvas"), myOptions);         var ctaLayer = new google.maps.KmlLayer('thedata.kmz');         ctaLayer.setMap(map);              //cluster         var mcOptions = {gridSize: 50, maxZoom: 15};         var mc = new MarkerClusterer(map, ctaLayer, mcOptions);     }      How can I cluster the markers in the overlay, just like markerclusterer does it?


title: solar radiation calculation

text: I have tried to calculate solar radiation (whole year with monthly interval) for Portugal from an ESRI grid (27m DEM) in Arcgis10. It seems to work, however it takes much more than a week and in the meantime there was an energy failure. I use ArcToolbox/Spatial Analyst Tools/Area Solar Radiation. I have cut the dem so as to have only the south of Portugal in the hope that the processing would be faster. However, ArcGIs 10 crashes after a few minutes and I get the error: "Background processing has encountered a serious error and is unable to continue. Contact ESRI". I have also tried in Arcgis 9.2 and I get the same message. It simply seems not to like the cut dem. Could anyone help? Teresa calvão


title: Opticks 4.11.0 (Ubuntu 64bit) - unable to load TIFF file

text: I am trying Opticks 4.11.0-1 - have official release opticks_4.11.0-1_amd64.deb but unlike in Windows I am totally unable to load Landsat 7 TIFF files. In Windows I can do this using auto importer. In linux it reports The dataset is empty" although I am able to load the same file on windows machine without problem. The package looks to be OK in Muon Package Manager (so not corrupted), gdal is installed (I also have QGIS) and Opticks do not show any error message during startup. Do you have any idea what could be problem? thanks


title: Which data providers support loading attribute tables without geometry?

text: I have a table in Memory from a cx_Oracle query (a list of rows, and a list of field names and their types) in my python plugin and I want to load it as a geometry-less layer (Attribute table) in qgis. Which data Provider will Support that? QgsOracleProvider has a constructor that Returns false, disconnects and sets valid to False if the table queried has no geometry column. QgsMemoryProvider only allows Point, Linestring, Multipoint, MultilineString, and Multipolygon geometries. Any other geometry is interpreted as WKBUnknown, which causes the `QgsMemoryProvider.isValid()` to return False. If I cheat and create a `QgsVectorLayer("Point", "myLayerName", "memory")`, and add features with no geometries, then I get a seg fault when I add the layer to the map layer registry. Also important to note, is that when a layer is imported using QgsMemoryProvider, it is not possible to Join against it (A Join from the layer does work.) Since mostly people want to load non-geometric tables into Qgis in order to join them to an existing vector, this reduces the Memory provider's suitability for non geometric tables. (Unless I am wrong. I would be delighted to be wrong.) QgsDelimitedTextProvider requires an Input file, and I would much rather not write out a csv to the file and then Import it again when I already have the data. But maybe I have to? Thanks


title: Compare two tables. Show values only in one table

text: I have two tables. Each table has a column with a reference number. Most of the records in these tables are the same but there are some which are in one table but not in the other. How can I compare these columns and pick out the records which are only in one table? I could write a MapBasic program that loops through each record and looks it up in the other table but I'm wondering if there is a simpler SQL-like way of doing this.


title: Open software tool or QGIS plug-in to make land cover classification

text: Is there any tool or plugin for QGIS/GRASS to make an automatic land cover classification from satellite imagery --- based on CLC nomenclature for example? Of course, it will not be 100% sure but it would be a first step of the classification process. I have seen a similar tool expressed into a pdf file on the QGIS blog, but it seems the plugin has not been implemented (at least, not listed into the plugin repository). I have also the existence of the Molusce plugin. I haven't tested it yet but the approach seems to be based on land use changes.


title: QGIS layers no longer overlay when Enable 'on the fly' CRS transformation is on

text: My QGIS project has two layers - one retrieved from GADM (administrative boundaries of Guinea Bissau) and another one of the villages collected in the field using GPS Decimal Degrees - D.ddddd. If I have in Settings - Project Properties - CRS - Enable 'on the fly' CRS transformation clicked on, these layers no longer overlay. If it is off it is not possible to know distances in Km. It would be great if someone has a clue about what is going on... Thanks! Ila


title: ArcGIS Full-Text search

text: I've set up Full-Text searching on a table in my SDE database (ms sql backend). I can perform the full-text commands from SQL queries, but it doesn't seem to work in ArcMap or any rest services. I full-text supported at all in arcgis server? ie: if i run this SQL query on my server, it returns all trail names that have trail in it (trail, Trail, Trails, etc) SELECT * FROM dbo.TRAILS WHERE CONTAINS(Webname, '"trail"'); if i add CONTAINS(Webname, '"trail"'); to the "WHERE" section in arcmap, querying the same table, it gives me an error, saying full-text indexing is not set on that table. EDIT: i've got it so it doesn't give me the error anymore, but the CONTAINS(Webname, '"trail"') doesn't search for any iteration of trail in the field, it only returns the exact text i was looking for (webnames with "trail" in it, not "trail", "Trail, "Trails", etc) and ideas?


title: What's the proper way of deleting versioned data from an SDE geodatabase?

text: I haven't really read this anywhere. What is the proper way of deleting versioned data from an SDE geodatabase to make sure nothing is left behind in the A & D tables? Do I need to unversion the data before I delete it? Do I need to remove the user versions? Or does it matter? Thanks.


title: Call MRT Swath from IDL

text: Do you know how to call MRT Swath from IDL? I would like to work with Modis L1B (MOD02HKM), read it, reproject, some processing and writing the results. Any help highly appreciated.


title: Specifying a file for VirtualShape in SpatiaLite

text: I'm ~~fairly~~ very new to spatialite, and am stuck on something very basic. I'm creating a virtual table from a shapefile using virtualshape. This works fine when the shapefiles are in the same folder as spatialite.exe. `CREATE VIRTUAL TABLE test USING VirtualShape(test, CP1252, 28356);` I can't figure out how to specify an absolute path ie `C:\1. Project Data\1814\test` `CREATE VIRTUAL TABLE test USING VirtualShape("C:\1. Project Data\1814\test", CP1252, 28356);` Fails. What's the syntax that's required here?


title: Identify all layers on the map without pre-designing the infoWindow

text:
This applies to the ArcGIS Server JavaScript API 2.4. The Identify Task sample shows how IdentifyTask can be used to find features in multiple layers, as specified in the IdentifyParameters:               identifyParams.layerIds = [0,2];      In this sample, the infoWindow already contains elements for the expected layers, which are used to display multiple features in separate tabs of the infoWindow:

          //In the JavaScript:     dijit.byId("bldgTab").setContent(layerTabContent(bldgResults,"bldgResults"));     dijit.byId("parcelTab").setContent(layerTabContent(parcelResults,"parcelResults"));      **Is there a way to handle this dynamically, and create the required tabs at run-time?** For example, I wish to run IdentifyTask on all visible layers, and return the results in a separate tab of the infoWindow. How can I create the tabs dynamically?