title: ArcGIS JS API Edit Toolbar permits MOVE but not EDIT_VERTICES or ROTATE - JavaScript error thrown

text: I'm putting together an ArcGIS (v10.1) JS API (v3.2) application and I've hit a snag in drawing / editing functionality. Users can draw a new geometry or edit an existing one. I'm using the draw and edit toolbars to achieve this, but attempting to activate the edit toolbar is producing some strange errors and I don't know where to start debugging. This function works fine for moving a Graphic around on the map:               editGraphic: function(graphic) {              this.editToolbar.activate(             esri.toolbars.Edit.MOVE,             graphic,             {allowAddVertices: true, allowDeleteVertices: true} // used for EDIT_VERTICES action         );     }      However, if I change the action to `esri.toolbars.Edit.EDIT_VERTICES`, `esri.toolbars.Edit.ROTATE`, or any combination I get               Uncaught TypeError: Cannot call method 'add' of null _VertexMover.js:19      I've managed to do some digging (it's not easy as ESRI has everything uglified and Chrome's developer tools crash when I try to pretty print the API to support tracing the error. It seems the `null` object in question is a property called `_scratchGL`, which I assume should be the Graphics Layer, but is null (note it's not `undefined`). The edit toolbar is defined as `this.editToolbar = new esri.toolbars.Edit(this.map)` and obviously the map is what it should be. I discovered while debugging other errors that the Graphic should be passed to my edit function via a Graphics Layer event, so the `graphic` variable passed to my function is sourced in this way (I think this is necessary to ensure key properties are set on Graphic, such as the Graphics Layer it belongs to). Unfortunately due to the size and many files in my application it's just not possible to show a useful sub-section of the code and I can't provide the whole lot. The fact that I can use the MOVE action tells me that I'm doing things correctly, so why can I not edit vertices or rotate? Fundamentally I'm not doing anything different to this example, which uses the same version of the API. However it is perhaps significant that if I add the following line of code it does not respond to clicks, even when all other map click handlers are disabled               dojo.connect(this.map.graphics, "onClick", function(evt){         debugger;     });      Has anyone had a similar experience? I've Googled this quite a bit but to no useful outcome. I don't even know where to start on the problem. Any thoughts very much appreciated.


title: How to hide some parts of an image layer with OpenLayers?

text: # **My Problem** I'm currently displaying a layer image above a _Google Maps_ base layer. I would like, when the user check/uncheck some checkboxes, that **some parts** of this image layer _disappear_ , and the **base layer** _appears instead_ these parts.   **But I want the other parts still visible with the same opacity.** * * * # **Explanations** The image parts to hide are **polygons**. _Those polygons are already generated as layers_ and are above my image layer. I could put a black background color to hide the wanted image parts, but it wouldn't display the base layer instead...   I'm stuck. * * * # **Solutions I thought about** I would like you to say me if that's possible and, if you can, how to do it. _I'm not asking for some code, I can do it myself, I'm just kind of a newbie at OpenLayers and I'd like someone to help me to understand how it works. Plus, I have to add functionnalities to an already written code so it's not very easy._ Here are data I have, I removed the ones which are not useful (this is only an example, empty data are usually filled with correct values). The images :               images[id] = {         name: "Image",         image_date: '2012/08/25',         image_path : 'static/path/to/image',         width : '250',         height : '250',         georef_north : '',         georef_south : '',         georef_east : '',         georef_west : ''     };      And the polygons :               polygons[id] = {         name:'',         longitude:'',         latitude:''     };      I'm getting them like that via WFS protocol like this :               layers['blockLayer'][i] = new OpenLayers.Layer.Vector("Block", {             strategies: [new OpenLayers.Strategy.BBOX()],             filter: VMMS.GIS.GetBlockFilter(i),             protocol: new OpenLayers.Protocol.WFS({                 url:  "../vine/jsGeoserver",                 version: "1.1.0",                 featureType: "block",                 featureNS : urlGeoserver + dataStore ,                 featurePrefix : dataStore ,                 srsName: "EPSG:900913",                 //srsName: "EPSG:4326",                 geometryName: "block_geom",                 schema : urlGeoserver +"wfs/DescribeFeatureType?version=1.1.0&typename=" + dataStore +":block"                  }),             projection : proj_spherMerc,             //projection : proj_WGS84,             styleMap : StyleBlockWithOpacity         });      ## First Solution > Why wouldn't we put the image layer in each polygon layer, and display only > the layers we want to display ? Would that work ? I'm not actually sure about that and I don't want to get into it if it's not realizable. Because I don't know if I can display an image IN a layer and not entire, and it has to be well placed in the map and in the layer... Though it seems to me it's doable. ## Second Solution > Why wouldn't we display the base layer again above the image to hide some > parts ? The image would still be here, but as the base layer would be on top > of it, we wouldn't see it. This solution seems to me even less doable, or, at least, **very** heavy. An example image, it's what I'm currently displaying : There's the google maps base layer, then the image blue and green, and all forms are polygons layers above the image.   Screenshot * * * Thank you all for your answers !


I have a set of coordinates that creates a polygon with holes in it. I'm using qgis with mmqgis to extract. How does MYSQL accept the WKT representation of a polygon with holes? By holes I mean areas within the perimeter of the polygon that belong to that of another polygon, they're voids in the shape. Is it with something like GeomFromText('MULTIPOLYGON((0,0 3,3 0,3 0,0), (1,1 1,2 2,2 2,1))')? The first set of coordinates is the main polygon, the second set is the inner?


title: memory issue when trying to buffer/union large dataset using postgis

text: I am trying to generated a dissolved buffer of a large point dataset (Ideally 29 million points - Address data in Great-Britain, but I receive the data by chunks of 1 million points, so these could be used instead). After trying (and obviously failing!) different methods using ArcGIS, I came across this article: http://blog.cleverelephant.ca/2009/01/must-faster-unions- in-postgis-14.html Since I have been willing to try postgis for a while a gave this a go. Obviously being an absolute beginner I struggled a bit to get the whole thing up and running, but I finally managed to import data and buffer it. I am struggling however with the ST_Union call, because of memory issues (I am using postgis 2.0.1 running on pretty old 32bits windows pc with little scope to change that in the near future). I tried changing the memory settings of postgres, as described in the question below: Queries returning very big datasets in PostGIS I believe I also created a spatial index properly on my geometry column. As far as I can tell, it didn't help. I managed to run ST_Union on a 500 000 points subset, which is already pretty big. I think at that point the best for me might be to try using the Declare/Fetch method mentioned in the above question, however being a complete beginner I am having difficulties. I believe it is in my interest to keep the slices as big as possible (because the ST_Union algorithm is fast and optimised). I think I would therefore just need to find a way to "slice" my ST_Union() queries to output one multipolygon for each 500000 rows. I would then be just one ST_Union on the resulting multipolygons away from the holy grail... It is also possible that my query is sub-optimal. Or, there might by another great simple solution that I just haven't thought about or found online yet... Here is a query I am currently using:                   drop table if exists ab_buffer_unioned;         create table ab_buffer_unioned (         gid serial primary key         , the_geom geometry         );              insert into             ab_buffer_unioned (the_geom)         select             st_Union(the_geom)         from             (select the_geom from ab_buffer limit 500000) as absubset;      Does anybody have any idea on how to fix my issues, or modify that query to add a loop or a cursor? Many thanks **UPDATE** Many thanks to Paul for his suggestion below. The GeoHash trick worked great as can be seen below (see how the cluster size varies with the population density!) ![ST_Union applied to slices of data sorted by GeoHash](http://i.stack.imgur.com/BFQ5w.jpg)


title: OpenLayers: Click event to display position shows Uncaught TypeError: Cannot call method 'getLonLatFromPixel'

text: I have created a map using openlayrs, and put a click event to display its position in screen, when i clicked got an error message like **Uncaught TypeError: Cannot call method 'getLonLatFromPixel' of undefined** , another error like **Uncaught TypeError: Cannot call method 'xy_to_ll' of undefined.** , I've tried OpenLayers.Control.Click but it shows some other errors, like e is undefined. And I've another problem with openlayers google map, on setting the mapCenter wms,yahoo maps shows the exact location, but google won't. How to fix this.I have figured out one solution for this to set sphericalMercator, but this doesn't fixes that problem My code is given below:               //canvasId is the id of html element.      Mapper =             {                 MapCanvas: null,                 position: null,                 MapperInit: function (canvasId) {                     /// <param name="canvasId" type="String">Id of the element on which the map is to be generated</param>                          this.MapCanvas = new OpenLayers.Map(canvasId, {                         layers: [                             new OpenLayers.Layer.Google("Google", { 'sphericalMercator': true,                                 'maxExtent': new OpenLayers.Bounds(-20037508.34, -20037508.34, 20037508.34, 20037508.34)                             }),                             new OpenLayers.Layer.WMS("WMS","http://vmap0.tiles.osgeo.org/wms/vmap0", { layers: "basic" }),                             new OpenLayers.Layer.Yahoo("Yahoo")                             ]                     });                     //Adding controls to map                          this.MapCanvas.addControl(new OpenLayers.Control.LayerSwitcher());                     this.MapCanvas.addControl(new OpenLayers.Control.MousePosition());                     this.MapCanvas.addControl(new OpenLayers.Control.OverviewMap());                     this.MapCanvas.zoomToMaxExtent();                          //sets the map center                     this.MapCanvas.setCenter(new OpenLayers.LonLat( 77.5, 8.5), 5);                          this.MapCanvas.events.register("click", this.MapCanvas, function (e) {                         console.log("Value of E => "+e.xy);                         this.position = this.MapCanvas.getLonLatFromPixel(e.xy);                         alert(this.position);                     });                 }             };      Please help me to find what is wrong in my code.


title: Is there a RESTful way to cancel a GP Job?

text: I have an asynchronous gp service that I can successfully run and cancel from within Arcmap. Using fiddler I see that arcmap communicates with the GP service via SOAP. I'd like to cancel the job using a REST request. However I don't see anything in the documentation for the GP Job describing how to do this. Does anyone know if there's a RESTful way to cancel a job? If not, would there be any bad side effects from writing another web service that scans the log file(s) to find the PID of the gp job then kills the process? (The gp service is configured non-pooled). **Update** ESRI support confirmed no RESTful cancel method (incident #879881). They suggested using the SOAP cancelJob method instead.


title: Batch Assign Coded Domain Values with Multiple SubTypes

text: I have a File GeoDatabase (v10.0) with multiple layers. Each layer has a range of Subtypes (1 to 20+). I have created coded domain values that I would like to assign to their appropriate field. I know the manual approach, but I am hoping to find a more automated system. Does anyone know of a way, tool, plugin, etc. to batch assign the coded domain values to fields that have subtypes? UPDATE: Attached is a graphic of the properties of one of my layers. In the Subtypes section you can see that there are many Codes (205, 206, 207, 224, etc.). In the Default Values and Domains section there are four (4) fields of interest (TypeComp, Type2Comp, Layer, and UserName). In order to have the correct Domains associated to the correct field, I need to associate the Domain to the field for each Subtype code. I want to know if there is an automated way to do this, as I have 13 layers, and each layer has multiple Subtype codes. ![enter image description here](http://i.stack.imgur.com/5e02j.jpg) Thanks...


title: WMTS or TMS Map Labels

text: Is there another WMTS or TMS offering labels for cities, peaks, rivers, lakes, etc. like the stamen-toner-labels? I did some research but could not find anything. It has not necessarily to be for free and I am searching for something already existent (=> atm I don't want to implement it myself with OSM data). Thanks for any hints.


title: Getting Default symbology in a feature class with subtypes

text: I was wondering if there was a way to get default symbology in a GDB feature class with subtypes? I don't want to create a layer file for each feature class. Is there a way to have the symbology built in?


title: OpenLayers: how to style WFS layer with SLD file?

text: I'm experimenting with a simple web mapping service build on PostGIS/GeoServer/OpenLayers stack on Win 2008 Server. I overlay vector layer from PostGIS using:               var Index1v = new OpenLayers.Layer.Vector("Index1 v", {         strategies: [new OpenLayers.Strategy.BBOX()],             protocol: new OpenLayers.Protocol.WFS({                 version: "1.1.0",                 url: "http://localhost:1979/geoserver/wfs",                 featureType: "bern_poly_s100",                 featureNS: "http://www.ispm.ch",                 srsName: "EPSG:900913"             })     });      Layer displays correctly, but in plain, default[?] style from GeoServer. I've got a SLD style file build with uDig and it seems to be correct since WMS service is using it without any problems. Can I use it to style WFS layer as well?


title: Mapping from one plane on the other plane despite of masking regions

text: have a set of data given here, where in the first and second columns there are the sky coordinates (ra,dec) e.g. right ascension and declination, respectively and in the third and forth, the coordinates in a Cartesian system (x,y), e.g. positions on the CCD image. ![enter image description here](http://i.stack.imgur.com/LaGnM.png) I need to make a two-dimensional interpolation surface using coordinates x and y and another using Ra and Dec. The problem is the existence of masked regions, as shown in the figure above. That is what I so far tried and didn't give the right answer:               from scipy.interpolate import griddata     import numpy as np     import matplotlib.pyplot as plt          data = np.loadtxt('test.asc')          ra = data[:,0]     dec = data[:,1]     Xpos = data[:,2]     Ypos = data[:,3]          xi = np.linspace(Xpos.min(), Xpos.max(), 1000)     yi = np.linspace(Ypos.min(), Ypos.max(), 1000)     xi, yi = np.meshgrid(xi, yi, copy=False)     ra_int = griddata(data[:,2:4], ra, (xi.flatten(), yi.flatten()),                       method='cubic')     dec_int = griddata(data[:,2:4], dec, (xi.flatten(), yi.flatten()),                        method='cubic')      Using griddata fails and return just `NaN` values. Is there any way to do this interpolation in order to estimate the values of Ra and Dec from a given x and y coordinates even in the masked regions (map from `x` and `y` to `ra` and `dec`)?


title: Curve from a Douglas-Peucker line

text: When generating a line using the Douglas-Peucker Algorithm the resulting object is a line with steep angles, I would like to "smooth" that line and making it look less "squarish" using curves, like bezier curves: there is an algorithm for that ?


title: 3D analyst Enclose Multipatch

text: I have a point file with elevation data. Some of the points are elevations relating to soil horizons. I've created a tin surface for each soil horizon then used the "extrude between" tool to create a 3D polygon of the area between each surface, basically creating a set of blocks stacked on each other representing the different soil horizons. I now want to select points that fall within each soil profile. However, my multipatch is not closed so using the tool "inside 3D" does not give me all the points that fall inside each soil profile. I have about 5 different profiles, 1 of which is closed and the others are not. I've tried the "enclose 3D" tool but it adds volume in the z direction, misrepresenting the soil profile. Why is the extrude between tool not enclosing all multipatches and how can I fix this?


title: Missing Raster/Vector Layers after converting KMZ using KMLToLayer tool of ArcGIS Desktop?

text: I am trying to write a tool that will let users upload a KML/KMZ file and then be able to view it in a client application (something like an esri js api or Flex app). I am ok with storing the data in a feature class (file, personal, sde) once it is uploaded. I have tried writing a Python script using the KMLToLayer geoprocessing tool (I'm in 10.04) but I'm getting spotty results. My test kmz has two vector layers and two raster layers but my output has only one vector layer and nothing else. Are there other tools I could use?


title: Python resource for remote sensing?

text: I am looking for a good resource (MOOC, book, etc..) to teach Python for remote sensing applications. I am already familiar with http://www.rsgislib.org/ \- are there any other resource out there for a beginning Python programmer?


title: Change the background color of a AxTOCControl

text: I'm having trouble finding a way to change the background color of the AxTOCControl. Is there an easy way to do this? I'm using ArcEngine 10 C# and VS2010. Thanks


title: Download/create address locator styles for Europe?

text: Where can I find some address locator styles for European countries? or how can I create my own address locator that contains street, house number and city, and use their aliases? PS. I have the book: Customizing Locators in ArcGIS 10


title: Using matplotlib to create and export graphs automatically from ArcMap

text: I am going to start a project soon where, after the data has been processed in various ways, an output jpg showing a graph of the data should be produced. The workflow is as follows:               1. A feature class with several fields containing various values is created.     2. A radar/bar chart based on several fields within the feature class should be produced     3. Colour the chart according to whether the field value is within a certain threshold e.g. if the threshold for the field CHEM1 is 4, a feature with CHEM1 = 2 will be coloured red, while a feature with CHEM1 = 6 will be yellow.     4. Each unique feature within the feature class should have its own graph i.e. 10 features = 10 graphs     5. The graphs should be saved as jpgs     6. The jpgs should be added as attachments to the geodatabase the feature class is in - not needed at this point but should be kept in mind      I will be automating this process after the initial run, as the new data will be processed monthly and the jpgs generated automatically, so a python script will definitely used. Before I begin though, I am unsure as to which graphing method I should use. I have researched the following options for my approach:               1. Python + ArcPy graphing     2. Python + matplotlib     3. Python + Excel      My pros and cons for each option:   1. I am used to ArcMap's graphing options, and have been using arcpy and ModelBuilder daily for the past year, so it will be relatively simple for me to use the inbuilt arcpy functionality.   2. I have never used matplotlib before, but my brief research indicates that it is quite simple to use for basic charts, and it doesn't seem to be that much more difficult to create more advanced charts using it. Also, since it's included with ArcGIS, it seems like a logical choice.   3. I have worked heavily with Excel for many years, and am very familiar with creating charts in it. I also did VBA, and it is my understanding that after looking at the VBA code for an action, one can translate that into Python quite quickly.  I am leaning towards the third option, since I feel most comfortable with that combination. However, since this process will have to be automated completely and will eventually run on the server, I do not like the idea of having to install Excel on the server so that this process can be run. I also fear that all the arcpy code will slow the process down (from a combination of not having enough experience with the graphing methods, not being able to optimise my code, slow computers in general, network issues), and that even though matplotlib appears to be the simplest/fastest approach, I will end up taking longer writing the code because I have not used it before. My question is not which combination works the best overall (because that would be subjective/too localised), but which option would be the best combination to achieve the outcome for this particular task?


title: PHP Mapscript get feature / shape by index in layer

text: I upgraded my Mapserver version from 5.6.x to 6.2.0 (yes, big version jump here!) and I am stuck with getting a shapeObj from a layer by simply using its Index (or FID for shapefiles or GID in case of Postgis). In the old version I used the following:   `$layer->open(); $shapeOjbect = $layer->getFeature($fid); $bounds = $shapeObject->bounds;`      which returned to me the rectangular bounds.      Everything has changed in the new version of Mapserver/Mapscript. `GetFeature()` is not enabled anymore, and `$layerObj->getShape($resultObj)` can be used instead BUT getShape requires the `resultObj(nth result)` from a $layerObj->queryByXXX , but there is no XXX = Index (e.g. `$layerObj->queryByIndex($fid)` ).      There is a $mapObj->queryByIndex(layerIndex, tileIndex, $fid, [$addToquery]), but I have not seen any examples using that function.      I hope I explained clearly. Does anybody have any experiences on using this function?      Thank you   Francesco


title: What's the relation between QGIS and GRASS?

text: I have been working with GIS and remote sensing in an university and professional environment some ten years ago (ArcInfo etc), so I have some experience even though somewhat outdated Right now I am looking for an GIS application for a small farming project and I found GRASS and QGIS projects which sound both very interesting. However I do not really get the difference between the two programmes. Which one do I use best for simple thematic maps with different layers? Thank you very much in advance.


title: How do I reproject a arbitrary coordinate system grid to OSGB

text: I have a polygon grid of a certain shape that has some data associated with it. Unfortunately its not properly projected, being in what seems to be an arbitrary scale and coordinate system (it was originally a DXF). I need to reproject this. The grid cells are about 8 * 11 "units" (I'm going to guess metres - it looks like a "projected" system). I need them to be 3000 * 4000 metres (OSGB). How do I reproject this grid into the desired new grid? Once its been "resized", I can then easily offset it to georeference it properly. I have access to FME and ArcGIS. Thanks


title: How do I add polygons to an existing attribute table where the data is stored as points? ArcGis 10

text: I have a shape file of golf courses saved as points, I would like to digitise the boundaries and add the polygons to the existing attribute table. I have opened editor and clicked on the point data shape file but have no option to create a polygon. Thank you


title: How to export vectors from QGIS to MapInfo?

text: I am trying to save some vector lines that I made in QGIS (1.8 Lisboa) to a format that can be added to MapInfo Professional 6.0 (ancient version I know..). When I save the vector to a Mapinfo format it produces two files: .MIF & .MIDI. However in the Mapinfo if appears as though a vector file has 5 separate files: .ID .IND .MAP .MIF .DBF......and sometimes a .MIDI file. Can I not add the vector line without these other files? If so is there any way that I can produce them? Thanks!


title: What's long binary data?

text: I'm trying to turn data in a Microsoft Access file into an ArcGIS-10 layer. The table contains a column labeled "shape." Each line in this column reads "long binary data." What's that?


title: How to add ECW support to GDAL bundled with PostGIS Windows

text: I'm currently trying the new raster functionalities of PostGIS 2.0 and I successfully imported some GeoTIFF file into PostGIS (although it only find one band so it's all black and white for some reason). I used raster2pgsql.exe bundled with PostGIS 2 to store it. Most of the raster data i have is in ECW format and I would like to import them directly with raster2pgsql. I know i can convert my ECW to GeoTIFF but i would like to do it without conversion. When i list the supported input format of raster2pgsql (raster2pgsql -G), it does not list ECW as it can not be distributed with the install. So i downloaded the ERDAS SDK following these instructions: http://trac.osgeo.org/osgeo4w/wiki/pkg-gdal-ecw And copied the DLL in my SysWOW64 system folder (i'm using Windows 7 64-bits). But it still won't list ECW as an input format. When looking up online, i found that PostGIS uses a bundled GDAL found in C:\Program Files (x86)\PostgreSQL\8.4\bin\libgdal-1.dll I don't know where to copy the DLLs so that the GDAL bundled with PostGIS will grab them and accept ECWs rasters. Or if anyone know how to tell PostGIS to use the GDAL install that comes with FWTools (C:\Program Files (x86)\FWTools2.4.7\bin) which is directly compatible with ECWs (path variables?) Thanks a lot for your help! Fabien


title: Looking for a template project for OSM data

text: I have a PostGIS 2 database filled with OSM data for my country via Osmosis. Is there perhaps a template/project (!styling) for QGIS i could download, where i would simply change the provider of data to my own database (different database name and credentials) so i wouldnt have to deal with the tedious styling. Thanks


title: Standards for creating a RESTful Geocoding webservice

text: I am about to start defining a new API which will allow our clients to perform Geo Coding lookups. I would like to follow any standards there are out there. However looking at a few examples:   * The Google Geocoding API   * Yahoo! Maps Web Services - Geocoding API They seem to be setup using their own standards... I was wondering if there is a recognised standard out there.. and if you could direct me to a site that uses the standard that would be great. Off the top of my head I can think of   * GeoRSS   * Atom with open search and geo extensions However I'm not sure if these are recognised as the official standard.


title: How to filter wrong GPS points?

text: I have a database of GPS points. There aren't any tracks, only points. I need to calculate some value for every 100 meters, but sometimes GPS gave a wrong coordinates that lies far from real GPS points, and instead of calculating values for a small square, I have to calculate it for a really big rectangular area. What is the best algorithm to filter wrong GPS points? I made a screenshot to help understand: ![!\[http://content.screencast.com/users/smirnoffs/folders/Jing/media/94624331-db6a-4171-bed9-e2183f953a1d/gps_error.png\]](http://i.stack.imgur.com/AoWjW.png)


title: Enlarge/roughen the Vector Size of Tiles

text: Im currently running a GeoServer and my task is to make the vectors rougher thant they are right now. I really dont know where and how the Geoserver processes the data from the DB and how it will be displayed. I want to do this to speed up the render time for loading the map. If the server uses 20 points right now, i want it to use 10 for example. I dont know where to start. If my information provided is not enough, just tell me and ill tell everything you want to know. Thanks alot.


title: Last edited information missing when used QGIS - Dxf to Shp

text: I am doing as constructed drawings which require the information to be in a GIS system and more specifically drainage as constructed drawings. I must have a separate file for Pipes and Pits. First I use Autocad, save as a DXF then use Quantum GIS to convert the DXF to SHP. So the problem is that the last edited piece of information from the autocad file is missing when i convert it. I can go back and change the last edited line, polygon or polyline and the last thing touched will still be missing. Help? I have looked on the QGIS website but haven't been able to find anything helpful.


title: How to improve classification accuracy for 3-band (RGB) imagery?

text: I am hoping to extract some simple vegetation features (sagebrush) from 1m 3-band (RGB) NAIP imagery. Unfortunately, there is no near-infrared band available for this dataset and I need to use this particular imagery for a time-series analysis, so I am stuck with the 3-bands. If this were 4-band imagery, I would consider adding NDVI and EVI vegetation indices as ancillary data for the classification. I do plan on incorporating texture into the classification. What additional band indices or useful information from widely available data (e.g. NED, landform) can I incorporate into the classification to increase the accuracy? I am flexible in what classification approach that I take.


title: QGIS identify features button issue

text: Is there anyway to get the Identify Features button in QGIS to automatically select a layer and provide the information without having to highlight it first in the layers TOC? MapInfo and ArcGIS automatically highlight, but cant see that QGIS does


title: Map printing using OpenLayers

text: Following OpenLayers map printing demo here, I'm able to get my map in image format. But I'm facing two issues with the export:   1. Google Street base map is not included in the image output;   2. Whenever I activate any WMS layer, the export fails with the error `TypeMismatchError: The type of an object is incompatible with the expected type of the parameter associated to the object` To get an image output of the map, I'm using `Canvas.toDataURL()`, which opens another window with the image of my WFS layers. According to the demo, when using `Canvas.toDataURL()`, images should be requested through a proxy, a point I don't understand as far as Google base layer is concerned. Any clarification will be appreciated. I'm not sure what I may be doing wrong. Below shows how I retrieve Google base map and WMS layers: Google base map:               var gStreet = new OpenLayers.Layer.Google(                 'Google Street',                 {numZoomLevels: 19},                 {visibility: true},                 {isBaseLayer: true},                 {useCanvas: OpenLayers.Layer.Grid.ONECANVASPERLAYER}         );      WMS Layer:               var wmsLayer = new OpenLayers.Layer.WMS("My WMS Layer",             GEOSERVER_HOST + "/wms",              {layers: "wms_layer", format: 'image/png'},             {                 useCanvas: OpenLayers.Layer.Grid.ONECANVASPERLAYER,                 isBaseLayer: false             }         );      `GEOSERVER_HOST` is proxied via Apache web server. I've configured Apache to channel all GeoServer requests as necessary. As such, I'm not making direct request to GeoServer in my client. Any idea on what I may be doing wrong will be appreciated.


title: How do you treat partial 3d features in PostGIS?

text: We have features from survey data that contain partial 3d information. The most common example would be a 2D LineString representing a road, that contains elevation information in certain points where it was surveyed. Other examples include roof shapes - A MultiLineString where some key points have an assigned elevation from the building plan, but not all. Using PostGIS, which data model would you recommend to store this kind of information, to keep it as accessible as possible, without losing or generating interpolated information?


title: Creating a custom Lambert Conformal Conic projection in ArcGIS 10.1 and/or 10.2.2

text: referring to the table below, I am trying to use the values on the left from another program that uses different terms for the parameters of a custom projection. I just wanted to make sure that I've matched up the terms used in ArcGIS correctly. The left side under projection: Lambert Conformal is the information I have and under ArcGIS are the terms I am need to match. The way they are lined up is my best guess, and this is located in California.               Projection: Lambert Conformal           ArcGIS     Reference latitude:       36.00       Latitude of Origin     Reference longitude:    -120.50       Central Meridian     Standard longitude:     -120.50       Central Meridian     First true latitude:      35.00       Standard_Parallel_1     Second true latitude:     65.00       Standard_Parallel_2


title: Datum transformations, do we use wgs84 as an intermediate common reference?

text: I was transforming some data the other day using the 7 param helmert transformation and it got me thinking. I have a set of transformation parameters, but these only describe the transformation from WGS84 datum, to Bessel 1841 datum. In practice... using my software, i can technically transform this data from/to any other datum i want, but does this mean: that all datums that exist actually use WGS84 (or some other) as a common denominator to which they all have known transformation parameters, and that the transformation actually happens (in the software) in the form of Datum A -> WGS84 -> Datum B, because there are no known direct parameters from A to B. In my case, since A is already WGS84 it isnt noticable of course (WGS84 -> System B)


title: ArcGIS Linear Referencing successive overlay operations

text: I have 8 tables with RouteID, BeginMeasure and EndMeasure. I want to successively overlay them starting with the first two and then overlaying the result and the next table then the result of that overlay and the next table until all tables have been run. I am doing UNION as overlay type. For some reason when I get to the 4th overlay the tools stops overlaying and just seems to start appending. If two records have the same measure range and all other values are ok the operation doesn't seem to combine the records. Are Null vs blank strings considered conflicts for unioning records? EDIT: The steps are are to use the Result of the first overlay and as the input table and the next table to be overlayed as the overlay table (second parameter). Is there an ArcMap memory issue? What is the difference between the input table and the overlay table other than just attribute ordering? EDIT 2: Looks like there may be memory issues in Arcmap. After testing/executing 6 or 7 overlays with different datasets it seemed to break down and not want to work. Closed ArcMap and re opened and started working again.


title: Adding georeferencing data into tiff image

text: I want to create a georeferenced `.tiff` file from a `.png` file. So I convert the .png into a normal `.tiff` file using the following command line :               convert image.png image.tiff      Here, all works fine, I have a valid `.tiff` file, which displays fine. Then, I want to add some georeferencing data in it. To do so, I found a tool called `tiff2geotiff`, I'm using it like this :               tiff2geotiff -4 "+proj=latlong +ellps=WGS84 +datum=WGS84 +no_defs" -c none -n "-122.5575664 38.5818201 -122.5449425 38.5896175" notgeoreferenced.tiff georeferenced.tiff      This gives me a **black image** which, imported in QGis, displays fine in the good place, except that it's black everywhere (proof that the _georeferenced data have been well added_ ). Therefore, when I do :               gdalinfo georeferenced.tiff      I obtain :               Warning 1: TIFFReadDirectory:Bogus "StripByteCounts" field, ignoring and calculating from imagelength     Driver: GTiff/GeoTIFF     Files: araujo_tmp.tiff     Size is 1585, 979     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (-122.557563781738281,38.589618682861328)     Pixel Size = (0.000007966552118,-0.000007972639275)     Metadata:       AREA_OR_POINT=Area       TIFFTAG_DOCUMENTNAME=georeferenced.tiff     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  (-122.5575638,  38.5896187) (122d33'27.23"W, 38d35'22.63"N)     Lower Left  (-122.5575638,  38.5818135) (122d33'27.23"W, 38d34'54.53"N)     Upper Right (-122.5449368,  38.5896187) (122d32'41.77"W, 38d35'22.63"N)     Lower Right (-122.5449368,  38.5818135) (122d32'41.77"W, 38d34'54.53"N)     Center      (-122.5512503,  38.5857161) (122d33' 4.50"W, 38d35' 8.58"N)     Band 1 Block=1585x1 Type=Byte, ColorInterp=Red     Band 2 Block=1585x1 Type=Byte, ColorInterp=Green     Band 3 Block=1585x1 Type=Byte, ColorInterp=Blue      **I think the first warning may be the problem...** Is there someone who knows why or another way to do this ? I'm looking for solutions in command line or in `PHP`. Thank you ! * * * # Solution **I found out`gdal_translate` can add georeferencing data into `.tiff` images**, so here is the command line which convert from png and add georeferencing data (you have to replace the extents (`llx ury urx lly`), of course):               gdal_translate -a_nodata 0 -of GTiff -a_srs EPSG:4326 -a_ullr llx ury urx lly pngfile.png tifffile.tiff


title: GFS SIG/SFC file format of siganl.gdas/sfcanl.gdas analyzing/decoding problem

text: I am creating a database from datafiles. I encountered many fails and errors, which I have already fixed (also with the help of this Q&A web). Now I am stucked a bit in decoding of 6 files because I totally don't know in what format they are coded. Those files are:               sigf06.gdas.2010123100     sigf00.gdas.2010123100     siganl.gdas.2010123100          sfcf06.gdas.2010123100     sfcf00.gdas.2010123100     sfcanl.gdas.2010123100      They are binary and when I open them in hexa (in console with od -c "filename"), I see the first 3 files have > "\0 \0 \0 G F S S I G \0 003 017 =" on the first line and the next 3 files have > "\0 \0 \0 G F S S F C \0 003 017 =" on the first line and after that all those 6 files have binary data. All I know is that I downloaded them from directory http://nomads.ncdc.noaa.gov/modeldata/cmd_HIC/2010/201012/20101231/ (any of the 4 .tar.gz files) and the link to this directory I have from http://nomads.ncdc.noaa.gov/data.php#CFSR-data -> High Resolution Initial Conditions HTTP link. All the datafiles contain geo data in a grid, I think that those 6 files store land values. First 3 files have 227MB, the other 3 files have 117MB. I am already quite insane how non-documented the files are, but it's still the best source for the needed data to my database. There is probably some documentation of the files, but I haven't found any. There was about 20 dataset files in those .tar.gzs, the majority I already decoded. If I knew just the format of these files, I would convert them easily, but link to the documentation of these files with structure and filetypes would be the best.


title: Is the GEO tag for QR codes lat/long or long/lat

text: I just found out that you can add a geographic component to a QR code. I also just found you can create one on the Google Chart API using this: https://chart.googleapis.com/chart?cht=qr&chs=400x400&chl=geo:30.3,-97.6 That creates me a QR code that looks like this: ![enter image description here](http://i.stack.imgur.com/bT7F2.png) So far, so good. And when I read that back on my iPhone using the QR Code Reader app, it sends me to a Google map in the correct location (Austin, TX). The thing is, the string used to create this uses lat/long. If I wrote the same data to a KML file it would expect long/lat. So the question is, what spec do I use for QR codes:   1. Is the QR spec actually lat/long?   2. Is the Google Chart API getting it mixed up?   3. Is the QR reader app getting it wrong?  My money is on the QR reader messing up - but I'm not sure. Does anyone have a definitive answer?


title: Time not configuring correctly?

text: I'm trying to match data in a textfile that contains time, value, and fips code to the cntyidfp part of a map of US counties. However, even though I am able to match the two sets of data together, ArcGIS does not recognize I have multiple dates. It either says the date values are from 01/01/0001 (which isn't even close to correct), or that the only date that exists is 01/01/1999. The date range is one week, and I want to create an animation by the hour. Is there an error in how my data is saved or something? This is what the textfile generally looks like :               string_datetime,load,fips     "01/01/1999 0:00:00",14769,"01001"     "01/01/1999 1:00:00",14634,"01001"     "01/01/1999 2:00:00",14481,"01001"     "01/01/1999 3:00:00",14635,"01001"     "01/01/1999 4:00:00",14617,"01001"     "01/01/1999 5:00:00",15225,"01001"     "01/01/1999 6:00:00",15567,"01001"     "01/01/1999 7:00:00",15720,"01001"     "01/01/1999 8:00:00",15864,"01001"     "01/01/1999 9:00:00",15856,"01001"     "01/01/1999 10:00:00",15666,"01001"     "01/01/1999 11:00:00",15267,"01001"     "01/01/1999 12:00:00",14759,"01001"     "01/01/1999 13:00:00",14302,"01001"     "01/01/1999 14:00:00",13974,"01001"     "01/01/1999 15:00:00",13925,"01001"     "01/01/1999 16:00:00",14362,"01001"     "01/01/1999 17:00:00",15504,"01001"     "01/01/1999 18:00:00",15707,"01001"     "01/01/1999 19:00:00",15498,"01001"     "01/01/1999 20:00:00",15214,"01001"     "01/01/1999 21:00:00",14944,"01001"     "01/01/1999 22:00:00",14354,"01001"     "01/01/1999 23:00:00",13535,"01001"     "01/02/1999 0:00:00",13033,"01001"     "01/02/1999 1:00:00",12776,"01001"     "01/02/1999 2:00:00",12648,"01001"     "01/02/1999 3:00:00",12687,"01001"     .     .     .     .     "01/07/1999 18:00:00",7229,"56039"     "01/07/1999 19:00:00",7118,"56039"     "01/07/1999 20:00:00",6752,"56039"     "01/07/1999 21:00:00",6467,"56039"     "01/07/1999 22:00:00",6114,"56039"     "01/07/1999 23:00:00",5594,"56039"      I have also tried resaving the dates in various formats, including without quotes, or in a 1999-01-01 fashion. Thus far, however, nothing seems to work. Where am I potentially going wrong? The county map is the standard shapefile from the tiger website.


title: Conditional Statement for attribute data

text: I'm trying to write a conditional statement that references another field within the same shapefile layer. I want it to only do so if the one field is " Null " and leave is alone if there is already a number there. I would think there is a simple way to do this, but I'm having trouble. Did anyone have a similar problem they figured out?


title: Joining Multiple Features to One Feature

text: I have two layers - an index layer (polygon) and a town boundary layer (polygon). The desired output is a FC with each town represented as only one feature and that FC's table having a field with values including all index layer features that intersect each town. Is there one tool in ArcGIS that can accomplish this or am I looking at a multi-step process?


title: Clipping between features within a single vector?

text: Is it possible (via the Qgis API and friends), to clip polygonal features **within** a single vector layer? The hope is to accomplish this without having to duplicate the feature into a new layer. Brief experimentation with the _processing.runalg("qgis:clip"...)_ command only showed me capabilities to accomplish such between multiple layers rather than multiple features WITHIN one layer. Thanks!


title: How to automatically vectorize data from WMS?

text: I am facing a problem, when I have a historic map in WMS layer and I need to get buildings from this map (approx. 2000) to vector layer fo further analysis. Is there any workflow how to automatically vectorize these objects from WMS (for example with ArcScan)? I don´t want to do it all manually, because it is time consuming very much.


title: QGIS polygon areas inside a polygon

text: I'm having hard time doing the following analysis in QGIS: I have watershed polygons with individual IDs. Then I have four different shape files, each having different land use polygons. I would like to know the areas of each land use inside a individual watershed polygon.


title: How to align directions of connected lines?

text: I have a vector line made of different pieces of lines with different directions. I merged the lines to one. But this merged line still have different direction in its "invisible" parts. I want to plot a profile, so I need ONE direction. Any idea?


title: How to Interpolate very large datasets in QGIS

text: I am working with a very large set of point data, (XYAttribute) which I wish to display as a raster. I've been using Interpolate, this is successful but very slow. I have other programmes that simply display the data according to a colourbar, is there a faster/less data intensive means of converting my dataset to a raster?


title: Vector layer and google map do not line up in PDF

text: I'm using qgis 1.8 with windows 7. When I use print composer to export as PDF, the google map does not line up with the vector layers. Please help.


title: FME incorrectly reading SHP file schema in reader

text: I am trying to read in a SHP file using the FME SHP file reader and some of the fields are not being read in correctly to FME. This is causing random values to be introduced into the dataset. The file I am importing you are able to see the correct schema when viewing the data in QGIS. Where the highlighted field (No. 3) is Listed as a `Real` data type. ![QGIS Schema](http://i.stack.imgur.com/u1dif.png) Though when this data is read into the reader the data type is changed from a `Real` value to a `Number`. This is what i believe is causing th e data issues. See below the FME reader details,the highlighted row corresponds with the highlighted row above. In QGIS this data is correctly represented in the data table with Numeric values only, which is as expected. Then when the data is read in and pushed straight to a Viewer transformer the data is incorrect is some cases. Zero values (0) are represented as NN and in the image below are highlighted incorrect values. ![enter image description here](http://i.stack.imgur.com/UNuE0.png) How do i prevent this issue from occurring that is corrupting these values? Am I able to statically set the schema of a reader transformer?


title: How to install ArcGIS Code Assist in Aptana 3?

text: I have downloaded Aptana and **ArcGIS Code Assist** but I don't know how to "DRAG" the vsDOC into Aptana. > Aptana 3 Installation Instructions Download the zip file that contains the > VSDoc file for the version of the ArcGIS API for JavaScript you are using. > Drag the VSDoc file into your Aptana project. I drag it into my project workspace and it only shows the .js file content. How do I actually use it?


title: FME assign specific node id with OSM writer

text: Is there a way to reassign the node's id value to something I generate? I did try exposing the `uid` attribute in the write transformer however only produced a new attribute and did not affect the `id` value. //The question has been changed. To clarify the question at hand, How can I/Can I assign a node's id with a specific value when fme starts to write an osm file? Thanks.


title: wireless sensor network integration with GIS

text: I am working on a project in which I need to represent data on a map after collecting it from wireless sensors.The wireless sensors are placed in the field and it collects information regarding soil ph,moisture,temperature.It is fitted with a gps which takes the co-ordinates.These attributes are then transferred to a computer in the form of an excel sheet. I want to know if it is possible to take this data directly into Qgis via a pen-drive,bluetooth etc and plot it on the corresponding map and show the co- ordinates from where the data has been collected


title: Multiple labels in QGIS QML

text: I've a GML file from the Ordinance Survey which includes a Cartographic Text layer which I'd like to load and render in QGIS. The layer has attributes which are intended to be use for the font size, font type, its orientation, etc but it doesn't appear to be in an immediately useful format. For example, orientation is in 1/10's of a degree and font is simply an number between 1 and 3. Setting aside loading the GML into PostGIS and replacing the fields with ones that may be more readily interpreted, I was wondering if there's a way within the QML structure to support more than one label type - for example, each symbol links to its own label section?


title: How to programmatically detect a canvas refresh or pan/zoom in QGIS

text: For my python plugin I would like to remove a VertexMarker when the user pans/zooms there map. Which is the best way to detect a canvas refresh or pan/zoom carried out by the user? Can you suggest any plugins which have this functionality I can take a look at? Thanks


title: How to get values of each cell in raster attribute table?

text: How do I get the values of each cell in the raster attribute table. Usually the attribute table gives you the value and number of counts. But, in my case I want each cell listed and corresponding cell values.


title: Leaflet Custom menu links

text: I have an issue where leaflet when embedded in Drupal will not allow the user to click on the map as the Y coordinates have a decimal on for some reason. As a quick fix, I am thinking to put the map into full screen mode outside of Drupal. However I would like to place a link/button on the map that points back to my Drupal site. Any suggestions on placing custom buttons on the map?


title: how to limit points(of particular parameter) in shapefile in WMS?

text: I wanted to filter out the points in shape by passing parameter and display points on screen in WMS that pass the parameter. I am taking shapefile from geoserver in WMS. Can any one help me ot with the same. I was trying this is this write but not working out:               var filter = new OpenLayers.Filter.Comparison({                         type: OpenLayers.Filter.Comparison.EQUAL_TO,                             property: "type",                             value: voltn //voltn is variable                             });      var parser = new OpenLayers.Format.Filter.v1_1_0();      var filterAsXml = parser.write(filter);      var xml = new OpenLayers.Format.XML();      var filterAsString = xml.write(filterAsXml);      mah2.params["FILTER"] = filterAsString;      mah2.redraw();      Can any one help me out with this


title: How to extract small portion of GIS database from entire database

text: Is it possible to extract a small portion of data from an entire database? If possible then how can I do it. I want data to be extracted which is within my map view only, so that other data objects are excluded: only objects contained in the map view should be in the extract.


title: File GDB to SDE Database using FME

text: I'm currently stuck in my FME Workspace script. Ideally, I'm reading from a file geodatabase and will update the SDE database in SQL server. I have CRC calculators to check for the changes in the Attributes and Coordinates and if the 2 CRC values are different then i will use the data from the File GDB to update the SDE database. This is a screenshot of my workspace: ![enter image description here](http://i.stack.imgur.com/BQVOK.png) The issue that i am encoutering is it always updates the data because the geometries being stored in the SDE database is not the same with the geometries in the File GDB. I modified the SDE writer feature type to match the X,Y Origin and scale from the File GDB feature class and disabled the "Simplify geometry" but the coordinates stored in the SDE Database are still changing. ![enter image description here](http://i.stack.imgur.com/IehFn.png) These are the coordinates from the File GDB and the coordinates stored in the SDE Database: ![enter image description here](http://i.stack.imgur.com/hhLBo.png) Also both the SDE and File GDB have the same spatial reference. I am quite stumped about this. Tried alot of things but it still automatically changes the coordinates when saved to the SDE database. Any ideas will absolutely be helpful.


title: Some data Loss when Import PGDB Table from Excel

text: When import table from excel some data were loss. i was try to change field type but i cant solve this error. I have attached an example image to show the problem. ![enter image description here](http://i.stack.imgur.com/SC0Sm.png)


title: How to add a coordinate after add layer?

text: how can i edit a layer that already out of the coordinate? (example: already setting blank layer with appropriate coordinate, but when added a new layer, it look like there is negative coordinate...) can some explain to me how to solve this matter? TQ


title: OpenLayers change feature's style

text: Working with OpenLayers 2.11, I'm desperately trying to change my feature style as the user selects it (fairly simple I agree!). Indeed, It seems easy but I can't see where It fails, perhaps someone experienced with OpenLayers could enlighten me. Below is my code :                  function onFeatureSelect(event)         {              var feature = event.feature;              var styleWaypointSelect = new OpenLayers.StyleMap(         {             externalGraphic: 'img/map_marker.png',             graphicWidth: 16, graphicHeight: 26, graphicYOffset: -24,             title: 'selectedFeature'         })                   var clonedFeature = waypointLayer.getFeatureById(feature.id);         clonedFeature.style = styleWaypointSelect;              waypointLayer.drawFeature(clonedFeature,styleWaypointSelect);         waypointLayer.redraw(); /* refresh() I tried several things to refresh the layer but it doesn't seem to change anything at all*/      Thanks, Gallien LABEYRIE EDIT : I finally figured out that OpenLayers StyleMap class allows different properties, "default", "delete", "temporary", "select". (you can see this when you debug the style object in firebug) I added the "select" style of my vector layer like this :                       var styleWaypoint = new OpenLayers.StyleMap({                   "default": new OpenLayers.Style({                     externalGraphic: 'img/mobile_loc.png',                     graphicWidth: 16, graphicHeight: 26, graphicYOffset: -24,                     title: 'marker'                 }),                    "select": new OpenLayers.Style({                     externalGraphic: 'img/mobile_loc_select.png',                     graphicWidth: 16, graphicHeight: 26, graphicYOffset: -24,                     title: 'marker_select'                 })             });


title: Labelling Kriging result (contour)

text: I've finished Kriging, creating a contour-map, representing groundwater-height in my area. By a reason I don't know, it isn't possible to label my results. Neither can i export my data to vectors (error: "the input layer is not valid") how is that possible and what am i doing wrong? Thanks in advance![enter image description here](http://i.stack.imgur.com/qNNya.png)


title: Visualizing Layers with lyr-files in ArcGIS API for JavaScript

text: I do have a probably very simple question. I'm currently running a Python-Script on my Desktop-machine. In the end of this script the different layers that I've created will be visualized according to a specific Attribute in an Attribute Class. That works just fine. Now I want to use this script on my ArcGIS Server and build a JavaScript- Webapplication. Is there a way that this webapplication does the same Visualization process as the Python Script? Here is a short example of the python-script i'm using for visualization:               for dlr_check in dlr_checks:      if dlr_check.getValue("Attribute_Class") == "Attribute_A":         arcpy.ApplySymbologyFromLayer_management("Feature_Class", "path/to/lyr-file/Attribute_A.lyr")      elif dlr_check.getValue("Attribute_Class") == "Attribute_B":         arcpy.ApplySymbologyFromLayer_management("Feature_Class", "path/to/lyr-file/Attribute_B.lyr")      Or does that Python-Code work in a JavaScript web application as well?


title: Is there no EPSG code to use for Behrmann equal-area projection in ogr2ogr?

text: I'm trying to use ogr2ogr to reproject some shapefiles. The projection I wish to use is Behrmann equal-area. I cannot find an EPSG code to use for the '-t_srs' part of the syntax. Is there any way around this? The best I can find is this - http://spatialreference.org/ref/esri/54017/ Any help much appreciated! Update:               for file in $INPUT/*.shp ; do     ogr2ogr -t_srs PROJCS[\"World_Behrmann\",GEOGCS[\"GCS_WGS_1984\",DATUM[\"WGS_1984  \",SPHEROID[\"WGS_1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.017453292519943295]],PROJECTION[\"Behrmann\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"EPSG\",\"54017\"]] $OUTPUT/$file $file     done      This is now what I have. Apart from inverted commas, i'm not sure what else needs to be escaped in this statement. Any ideas?


title: CartoDB st_makeline with geography datatype

text: Hi I am using an online cartodb account to try and create a shapefile based on flight routes between one airport and another. So far, I have uploaded my lat/long points for each connected airport, and used st_makeline to create a line string between them in geometry datatype. However, this outputs a direct straight line, whereas it should really follow a geographic projection and output as a curved line. Does anyone have any advice on how to do this? Bearing in mind I am using cartodb online so am restricted to the custom SQL query it provides. Current query is as follows: SELECT routeid, boardpoint_offpoint, ST_Makeline(the_geom) as the_geom FROM routemapping group by routeid, boardpoint_offpoint Many thanks for any help that can be given!


title: Do you know a service route calculation to get a route between several points?

text: I need a routing service with an API callable from a web or desktop application. I have several points: a start, an end, and intermediate points in the disorder. I expect the routing service to send me the shortest path between the start and the end, via every intermediate point. Do you know such a service ?


title: Problem creating point shapefiles programmatically using python in QGIS 2.2

text: I am trying to create a point shapefile in the memory buffer (for a temporary label) and then zoom to the extents on this buffer. The working (in version 2.1) test code for this is:               def originalFunction():         vl = QgsVectorLayer("point?crs=epsg:27700", "temporary_points", "memory")         pr = vl.dataProvider()         # Enter editing mode         vl.startEditing()         count = 0         while count < 10:             # add a feature             fet = QgsFeature()             fet.setGeometry( QgsGeometry.fromPoint(QgsPoint((376000+count),(408000+count))) )             vl.addAttribute(QgsField("mytext", QVariant.String, "", 255))             vl.addAttribute(QgsField("mytext2", QVariant.String, "", 255))             fet.initAttributes (2)             fet.setAttribute(0,"hello")             fet.setAttribute(1,"hello")             pr.addFeatures( [ fet ] )             count+=1         # Commit changes         vl.commitChanges()         QgsMapLayerRegistry.instance().addMapLayer(vl)         canvas = qgis.utils.iface.mapCanvas()         qgis.utils.iface.setActiveLayer(vl)#This is done so the zoom to extent works         qgis.utils.iface.zoomToActiveLayer()#Unless this happens the extent seems to #include the origin         canvas.zoomByFactor(1)#this time the zoom works      This code should create a line of points diagonally and zoom to their extent in version 2.1. When I try and run the same code in the new version (2.2 primarily but also 2.3) it does not work (no points are generated, and the fix for zooming to extent no longer seems to be working either). I have looked at the cookbook here: http://www.qgis.org/en/docs/pyqgis_developer_cookbook/vector.html#modifying- vector-layers and tried to use that as a guide to rebuild the function. However my result (below) does not seem to be operational, despite deviating from the example code as little as I can (I found one of the functions in the example did not exist for the object it was being called on so had to replace it as shown in my comments). New Function:               def testFunction():         # create layer         vl = QgsVectorLayer("point?crs=epsg:27700", "temporary_points", "memory")         pr = vl.dataProvider()         caps = vl.dataProvider().capabilities()              if caps & QgsVectorDataProvider.AddAttributes:             res = vl.dataProvider().addAttributes( [ QgsField("mytext", \         QVariant.String), QgsField("mytext2", \         QVariant.String) ] )         print vl.updateFields()         count = 0         while count < 10:             if caps & QgsVectorDataProvider.AddFeatures:                 feat = QgsFeature()                 feat.initAttributes (2)     #            feat.addAttribute(0,"hello")#Did not work, gave error: object has no attribute 'addAttribute'                 feat.setAttribute(0,"hello")                 feat.setAttribute(1,"hello")                 feat.setGeometry( QgsGeometry.fromPoint(QgsPoint((376000+count),(408000+count))) )                 (res, outFeats) = vl.dataProvider().addFeatures( [ feat ] )                 count+=1         QgsMapLayerRegistry.instance().addMapLayer(vl)         canvas = qgis.utils.iface.mapCanvas()         qgis.utils.iface.setActiveLayer(vl)#no longer fixes the zooming to extent bug         qgis.utils.iface.zoomToActiveLayer()         canvas.zoomByFactor(1)         canvas.setExtent(vl.extent())      Thanks for reading, I apologize for the length but I wanted to include code so that others can replicate the problem and see where I am going wrong!


title: Caching WMS data

text: I have an mxd created from a wms server and I am using it as a basemap in a flex viewer app. In order to speed up renedering time I need to cache this. However I have been unsuccessful. After the scales have been created using arc catalog map service properties the manage map server cache tiles fails to manage tiles at all extents and levels. Is there some other process I can do to help speed up rendering?


title: Is it possible to create a custom Editor Object for Arcmap and tie it to a custom tool?

text: I have some spatial enabled tables tables in SQL Server 2008 and add them to ArcMap using custom layers(or a a spatial enabled OLE DB provider) successfully,but as you know,the editor toolbar of ArcMap doen't get enabled for these layers. Is it possible (and if yes how) to create a custom "Editor Object" for ArcMap which implements all necessary interfaces(like IEditor,IEditor2,ISnapEnvironment,...) and use this custom editor for these custom layers?and is it possible to create a "Sketch Tool" which can work with this custom editor and have same properties like snapping? Best Regards Reza


title: Registering local folders on ArcGIS 10.2.2

text: We have all our geodatabases stored on a dedicated ArcGIS server in a folder D:\GeoData. This folder is shared as \\\GISSERVER\GeoData. Any users working from their work stations have to reference this UNC path, however the server can reference this as local path D:\GeoData My question is, should the UNC path always be used, and local path NEVER used, even if publishing from the server? I am proposing that the data be moved to a different server so that there is never any ambiguity, i.e. everyone uses the same path, period. Presently, we have 2x registered folders for the same data source (folder). One is Local and one is UNC. It doesn't make sense to me. We never copy data to the server when publishing, always keeping a single source.


title: WKID (Well Known ID) for national projection for ArcGIS REST API

text: I have to use my national projection system which does not have corresponding WKID listed for ArcGIS REST API. Without it I am unable to use this projection system for my web application. Any pointers which can help me out to use my data in its native projection format? Thanks


title: GRASS - Will r.to.vect work on this scanned map?

text: Total GRASS noob here. I have several very large hand contoured maps to digitize, and this thread intrigues me. In the past I have spent hours upon hours digitizing these maps and revisions in ArcMap. If I can save myself and my client time and money by using `r.to.vect` in GRASS, that would be great. The question is, to those with a trained GRASS eye, would `r.to.vect` even work on a scan like the sample below? I need to get the curved lines out to vectors. I'm wanting some educated eyes on this before I invest the time in learning enough GRASS to get this done. Thanks. ![I need to get the curved lines out of this raster](http://i.stack.imgur.com/FiK0H.png)


title: Using GeoServer REST Configuration to add a datastore

text: I have used the example from the Geoserver REST Docs to add a datastore from at PostGIS database. I have used the following xml doc to add the new datastore:               <dataStore>       <name>my_datastore</name>       <connectionParameters>         <host>localhost</host>         <port>5432</port>         <database>postgis</database>         <user>username</user>         <passwd>secret</passwd>         <schema>public</schema>         <dbtype>postgis</dbtype>       </connectionParameters>     </dataStore>      The datastore was created with the connection parameters as I expected. But the **datastore.xml** created looks different from the datastore that I have created with the **GeoServer Web Admin**. The connection parameters I am missing is:   * Connection timeout   * validate connections   * max connections   * Loose bbox   * Expose primary keys   * fetch size   * Max open prepared statements   * preparedStatements   * Estimated extends   * min connections So my question is: What are the above mentioned connection parameters corresponding xml elements? EDIT: I am using GeoServer 2.1-RC1.


title: Problem with deleting file locks in ArcGIS 10

text: I did not realize that I should have been deleting cursors after accessing attribute tables until today. Up until now I never did this, but I've noticed that there are lots of file locks after my script runs (I presume it was a bad practice on my part to begin with). ArcGIS 10 help tells me that I should use the following code:               # Delete cursor and row objects to remove locks on the data      #      del row      del rows      I'm doing that in the following function:               # calculate areas     def area_calc( shapefile ):         "calculates geometric properties of each building"                    # loop through attribute table             Rows = ARCPY.UpdateCursor( shapefile )               for row in Rows:                  geom = row.Shape                  polyArea = geom.area             row.setValue( perimName, 0 )  # set to zero             polyPerim = geom.length                  row.setValue( areaName, polyArea )   # update area value             row.setValue( perimName, polyPerim ) # update perimeter value                         Rows.updateRow(row)                 # Delete the row and cursor             del row          del Rows              ARCPY.AddMessage( 'Area and Perimeter calculated for each polygon' )      However, when I do that I get the following error message:               Traceback (most recent call last):       File "LoopModuleBuildingAnalysis.py", line 9, in <module>         import ConstructionMaterialPhase1       File "D:\Code\AnalysisPhase1.py", line 101         del row           ^     SyntaxError: invalid syntax     Press any key to continue . . .       Any suggestions what I should be doing? Thanks in advance.


title: Trouble creating a query table

text: I need to join one feature class and one table. This is a one to multiple join, therefore, my only option seems to be by creating a query table. My problem is that the fields I want to join do not seem to be matched properly when I run the query table tool. I've tried using the following expression in the query table expression field to indicate that these are the fields to be joined: 'FishingAreaGroundFishNafoUnit.Code' = 'FishArea_Code_homeports_TotalValue.Code' However instead of returning a table with just a few hundred records, I get one with 42,000 records which makes no sense at all. I think I may be doing something wrong (perhaps the problem is with my query), but I'm not sure what. When I verify my query in the query builder window, I keep getting a message saying 'The expression was verified successfully but no records were returned' which does not make any sense as there should definitely be some matches. I've also been using his thread (bottom post) to guide myself in creating the query table http://forums.esri.com/Thread.asp?c=93&f=1728&t=172344in Thanks in advance for your help! Léa


title: How to Buffer with intersections in QGIS

text: I am buffering a layer of facilities (points). I will then use the buffer to clip and intersect another layer, containing resources. I want to answer the question, how much of the resource lies within X distance of the facility? Buffer facility -> clip resource so only areas within buffer remain -> intersect facility buffer and resource to show amount of resource within distance of facility. Now, I am using ftools to create my buffer. The problem is, when the facilities are already within X distance of each other it does not create a full circular buffer, the polygons created do not overlap. How can I create a buffer in QGIS where the polygons of the buffer are allowed to overlap each other (possible in ArcGIS)?


title: Where is the open Java client for ArcGIS Server?

text: Short of hacking the Android SDK into your Java project, are there any Java clients that can consume and work with ArcGIS Server services? e.g. create Feature Sets, Geometries, perform Queries/Geoprocess tasks, etc.. Similar to what you can use the Javascript or ObjC SDK's. Though I don't need to visualize, only interact with the services. The closest thing I've seen is the Geometry API: https://github.com/Esri/geometry-api-java No, I don't want ArcObjects or any SDK that requires some $$license.


title: Hover/Popup feature problem

text:
im really having an issue here. I am trying to realise a map where i can hover over the features and they get highlighted (works), then i want to be able to click them (works) , i also want a div to change accordingly to the feature (partially works) and i want the popup to close. Im aware that this question sounds kind of easy, though ive found out that multiple events on a feature are not an easy task to implement. The map works as far as : I can hover a field, click on it, get a popup but heres where the fun begins. Once the popup is open i cannot close it anymore. Though the popup can be closed when clicking anywhere outside of the feature , so i think this has got to do with my clickout implementation. so far ive tried some stuff and im not seeing any progress. If the feature im hovering hasnt been clicked on and i hover out of the feature without clicking, the feature stays selected as if the mouse was still hovering. This effect goes away once the feature has been clicked once, but comes back immediately if you start the procedure again (so hover feature 1-hoverStyle,hoverout-still hoverStyle, click feature 1- clickStyle, clickout- standard style, hover feature2- hoverStyle, hoverout-still hoverStyle, hover feature 1- hoverStyle , hoverOut- no Style ). I hope you can get something out of this information. so now to the div change. i have it working on another layer, but in the layer with all the select and hover bla, it just doesnt want to work. Heres some code for you to understand or help me with my Problem:                   {                   function tooltipSelect(event){         var feature = event.feature;         var selectedFeature = feature;         //if there is already an opened details window, don\'t draw the          tooltip         if(feature.popup != null){             return;         }         //if there are other tooltips active, destroy them         if(tooltipPopup != null){             map.removePopup(tooltipPopup);             tooltipPopup.destroy();             if(lastFeature != null){                 delete lastFeature.popup;                 tooltipPopup = null;             }         }         lastFeature = feature;         var tooltipPopup = new OpenLayers.Popup("activetooltip",                 feature.geometry.getBounds().getCenterLonLat(),                 new OpenLayers.Size(80,12),                 " "+feature.attributes.name,                 true );         //this is messy, but I'm not a CSS guru         tooltipPopup.contentDiv.style.backgroundColor='ffffcb';         tooltipPopup.closeDiv.style.backgroundColor='ffffcb';         tooltipPopup.contentDiv.style.overflow='hidden';         tooltipPopup.contentDiv.style.padding='3px';         tooltipPopup.contentDiv.style.margin='0';         tooltipPopup.closeOnMove = true;         tooltipPopup.autoSize = true;         feature.popup = tooltipPopup;         map.addPopup(tooltipPopup);     }     function tooltipUnselect(event){         var feature = event.feature;         if(feature != null && feature.popup != null){             map.removePopup(feature.popup);             feature.popup.destroy();             delete feature.popup;             tooltipPopup = null;             lastFeature = null;         }     }          function onMarkerSelect(event) {         //unselect any previous selections         tooltipUnselect(event);         var feature = event.feature;         var selectedFeature = feature;         var popup = new OpenLayers.Popup.FramedCloud("activeAlarm",                 feature.geometry.getBounds().getCenterLonLat(),                 new OpenLayers.Size(100,100),                 "Loading...",                 null, true, onMarkerPopupClose );         feature.popup = popup;         popup.setOpacity(0.7);         map.addPopup(popup);         //call ajax to get the data         loadDetails(feature.attributes.description);     }     function onMarkerUnselect(event) {         var feature = event.feature;         if(feature.popup) {             map.removePopup(feature.popup);             feature.popup.destroy();             delete feature.popup;         }     }     function onMarkerPopupClose(evt) {         select.unselectAll();     }               //lType: Region-Type     //0=Alle, 1=Land, 2=Reg.Bezirk, 3= Kreis, 4=Gemeinde                  function loadLayersByColour (lType) {              // alert("in loadLaysersByColour");         // 1.) einmaliges Erstellen der Filterliste fuer WebService Request         if (listRun == false){buildFilterLists(6,0);};         // 2.) entfernen eventuell bereits vorhandener alter Layer         if (listRun != false) {removeLayer(0);};              // 3.) aufbauen der Layer              // zunächst  Bundesländer         if (lType==1 || lType==0) {                  // create a wfs layer with the style map             states = new OpenLayers.Layer.WFS(                     "Bundesländer",                     "http://myServer/geoserver/wfs",                     {typename: "myStore:layer"},                     {isBaseLayer: false, extractAttributes: true, styleMap: myStyleMapL}             );                   map.addLayers([states]);                  /*                 ctrlSelectFeatures = new OpenLayers.Control.SelectFeature(                         states,                         {                             clickout: true, toggle: false,                             multiple: false, hover: false,                             toggleKey: "ctrlKey", // ctrl key removes from selection                             multipleKey: "shiftKey" // shift key adds to selection                         }                 );                      map.addControl(ctrlSelectFeatures);                 ctrlSelectFeatures.activate();                      var aktLayer = 1; //index of a vectorlayer                 map.layers[aktLayer].events.register('featureselected',   map.layers[aktLayer], regFeatureSelected);                 map.layers[aktLayer].events.register('featureunselected', map.layers[aktLayer], regFeatureUnselected);                   */                  /*             if (typeof(clickCtrl) !== 'undefined' && clickCtrl != null ) {               alert("clickCtrl defined and not null ");             } else {alert("clickCtrl undefined or null");};               */                     // map.addControl(selectL);             // map.addControl(clickCtrl);                     // map.addControl(selectCtrlClick);             // selectL.activate();             // clickCtrl.activate();             // selectCtrlClick.activate();                  // fuer popup             //selectedControl = new OpenLayers.Control.SelectFeature(states,             //  {onSelect: onFeatureSelect, onUnselect: onFeatureUnselect});                  //AENDERUNG             //add the main select (for regular select)             selectedControl = new OpenLayers.Control.SelectFeature(states, {                 toggle: true,                 clickout: true,                 multiple: false,                 onSelect: alert("hier"),             });             states.events.on({ "featureselected": onMarkerSelect,                  "featureunselected": tooltipUnselect()});             //"featureunselected": onMarkerUnselect});                  //add the second select for tooltips:             var highlightCtrl = new OpenLayers.Control.SelectFeature(states, {                  hover: true,                 highlightOnly: true,                 renderIntent: "temporary",                 eventListeners: { featurehighlighted: tooltipSelect,                  featureunhighlighted: tooltipUnselect } });                  //the order in which you add these lines seems to matter!             map.addControl(highlightCtrl);             map.addControl(selectedControl);             highlightCtrl.activate();             selectedControl.activate();                        };              // Globale Variable für Merken der geladenen Layertypen  setzen         lStat = lType;           };               //  "
Feature: " + feature.id +"
Area: " + feature.geometry.getArea()+"
",         function onFeatureSelect(feature) {         selectedFeature = feature;         popup = new OpenLayers.Popup.FramedCloud("chicken",                  feature.geometry.getBounds().getCenterLonLat(),                 null,                 "
Feature: "  +"
Area: " + "
",                  null, true, onPopupClose);         feature.popup = popup;         map.addPopup(popup);     };     function onFeatureUnselect(feature) {         map.removePopup(feature.popup);         feature.popup.destroy();         feature.popup = null;     };          function onPopupClose(evt) {         selectedControl.unselect(selectedFeature);     };          function regFeatureSelected()     {         aktPopup = create_Popup();         var objF = this.selectedFeatures[0];         var desc   = objF.attributes['gen'] ? objF.attributes['gen'] : "about:blank";         $("nodelist").innerHTML = "Selektiertes Land " + ", " + desc;     }          function regFeatureUnselected()     {         if(aktPopup)             map.removePopup(aktPopup.popup);              document.getElementById("nodelist").innerHTML = "selected feature description";     }                    }


title: How to reclassify raster inside polygon mask, leaving values outside as they are?

text: I am in the final stages of a landcover classification and need to reclassify thematic raster values that are inside different polygon masks, while leaving the values outside as they are. Can I do this with raster calculator? And, if so, how?


title: WPF Runtime Arcgis - How to open MXD?

text: I would like to open MXD with ArcGIS Runtime WPF. There is no example about this. Anybody has a solution ? Thank you


title: WPF Runtime Arcgis - How to create a MPK?

text: I would like to create a Map PacKage with ArcGIS Runtime WPF. Have you an example ? Thank you


title: Why can't I multiply two string fields in the Field Calculator by converting them to int() and float()

text: I'm using ArcGIS 10 and new at python scripting. I have a string field that I need to calculate by multiplying two other string fields. I can do it in python using two arbitrary values...               CINPT1 = "1.00"     CNETVAL = "8500"          print int(float(CINPT1) * int(CNETVAL))      but...I need this as an arcpy fieldcalculation. What I have now is below. I'm basically trying to what is above: the CNETVAL field should equal CINPT1 * CUNITV               import arcpy, sys, string, os          MAPCHAR = "MAPCHAR"     CINPT1 = "CINPT1"     CNETVAL = "CNETVAL"     CUNITV = "CUNITV"          expression = "multiStringFunc(!CNETVAL!)"     codeblock = """def multiStringFunc(netvalue):         value1 = float(CINPT1)         unitvalue = int(CUNITV)         netvalue = int(value1 * unitvalue)              return netvalue"""          # Process: CUNITV = CINPT1 * CUNITV     arcpy.CalculateField_management(MAPCHAR, CINPT1, expression, "Python", codeblock)      What i get in return is "global name 'CINPT1' is not defined". I am pretty new at using python, so i don't know if I'm just using the value wrong within the codeblock or I'm just totally off. Anything you can contribute would be appreciated. Thanks!


title: Find gap in / endpoints of linestring

text: I have a set of features of type `linestring`. I'm trying to polygonise the linestrings and assumed that they were closed (naive me :) ). Of course some were not closed and I'm trying to locate the gaps or the endpoints of the individual linestrings, so I can close them manually. **Question** :   How can I locate the gap in or endpoints of an open linestring? _Available software: QGIS, FME, OpenJump, (ArcMap)._


title: How to add a label of a point to a polygon

text: i have several polygons with in each polygon a point with a label (number). I want to add the label in the polygon to the attribute-table of the polygon it belongs to.


title: Measurement units not returned in project CRS units

text: I am attempting to calculate area of country polygons using the field calculator. The layer in question is unprojected WGS84, but I have set the project CRS to World Equidistant Cylindrical (Sphere), which is listed with the projected CRSs. I've also enabled on-the-fly transformation. From my limited understanding of the program (I am running QGIS 1.8.0), the measurements should be given in the units of the project CRS, and the proj4 string for the World Equidistant Cylindrical seems to indicate that its units are in meters. Nonetheless, the measurements are off by orders of magnitude. For example, Afghanistan is significantly larger than 92 square meters, or even 92 sq km. I can only assume I am getting area in degrees. What must I do to change that?


title: FME autocad to csv

text: I have a few building blocks and land parcels in a dxf file that look like this ![enter image description here](http://i.stack.imgur.com/9hjPF.jpg) I am writing to a csv file using FME the following attributes: **land parcel number | xmin| xmax | ymin | ymax** The last four attributes correspond to the bounding box of each land parcel. I want to add the attribute 'building block' to each land parcel. I.e. in the above example Γ251 (=autocad_string) should correspond to each land parcel shown. **land parcel number | xmin| xmax | ymin | ymax | building block name** Any ideas on how make this work?


title: All UK census geographies in one file irrespective of scale

text: Is it possible to get all UK census geographies in one file? For example, I am currently looking at a table of data titled               Life satisfaction by country, region, unitary authority and county, April 2011 to March 2012      Which turns out to mean unitary authorities or counties, or local authorities in Scotland. Nevertheless each entry has its own code (e.g. `E06000006`). Is there a single shapefile that will match all of these without having to download all the geographies separately?


title: Myanmar Datum 2000 information

text: Does anybody know where I can find the technical details for Burma's latest datum? It's called Myanmar Datum 2000 and I can't find any details online. I need the details to create a custom datum on my Garmin GPS. What I need is the dx, dy, dz, da, and df. If this isn't available, what is the closest datum that could be used instead?


title: How to get full 16 bits range using RasterSymbolizer and ColorMap ramp?

text: I am attempting to get the full 16 bits values of my DEM (stored as TIFF files) rendered in WMS tiles (png format). For this I use GeoServer 2.5 and a RasterSymbolizer with a ColorMap of type "ramp". Then I used the color "#000000" for the value 0 and the color "#003E80" for 15000. My understanding was that GeoServer will interpolate as such: 1 -> #000001 2 -> #000002 256 -> #000100 15000 -> #003E80 etc... Instead it looks like it is taking the 16 bits values and divide it by 255 and apply the coefficient to each color. This in effect reduces my height resolution by a big amount! Is there a way to have the behavior I expected using a ColorMap (or any other sld operator)? Thanks.


title: How to add topology to PostGIS 2.0?

text: I have installed PostGIS 2.0 and need new topology but it doesn't seem installed by default. I found a lot of command line to install postgis 2.0 with topology but is there a way just to add topology? Thanks in advance.


title: processing with xslt-

text:
on opensuse 13.1 i try to do some gis-works with a large file: france- latest.osm.bz2 which i gathered from here: http://download.geofabrik.de/europe.html what do i do with that file france-latest.osm.bz2               bzcat france-latest.osm.bz2      what is aimed? i want to extract all things that belong to the POI restaurant which is long lat name adress etc - etx. i have the following things up and running: package perl-XML-Twig and run xml_split with a command available on openSUSE to split xml files named xml_split (it is part of the package perl-XML-Twig) Now we try to run the following command (I hope we have enough hard disk space since the output is roughly 20GB).               bzcat france.osm.bz2 | xml_split -s 100M -b france -n 3 -      this will result in a bunch of 100 Mb large xml files france-001.xml,france-002.xml and so on. Weu then have the xslt (the name of the root element) and of course we will need a loop in the bash to process the several files and collect all the results together. http://www.w3.org/1999/XSL/Transform'>               <xsl:output method="text" encoding="UTF-8"/>     <xsl:template match="/">                  <xsl:for-each select="xml_split:root/node/tag[@k='amenity' and @v='restaurant']">             <xsl:value-of select="../@id"/>             <xsl:text>	</xsl:text>             <xsl:value-of select="../@lat"/>             <xsl:text>	</xsl:text>             <xsl:value-of select="../@lon"/>             <xsl:text>	</xsl:text>             <xsl:for-each select="../tag[@k='name']">                 <xsl:value-of select="@v"/>             </xsl:for-each>             <xsl:text>
</xsl:text>         </xsl:for-each>     </xsl:template>      question: what do i need to get all the aimed data out of the dataset - i.e. long lat name adress etc - etx. here below we have a data-chunk out of the xml-file that we have parsed: see it               <node id="52768810" lat="48.2044749" lon="11.3249434" version="7" changeset="9490517" user="wheelmap_visitor" uid="290680" timestamp="2011-10-07T20:24:46Z">         <tag k="addr:city" v="Olching" />         <tag k="addr:country" v="DE" />         <tag k="addr:housenumber" v="72" />         <tag k="addr:postcode" v="82140" />         <tag k="addr:street" v="Hauptstraße" />         <tag k="amenity" v="restaurant" />         <tag k="cuisine" v="mexican" />         <tag k="email" v="info@cantina-olching.de" />         <tag k="name" v="La Cantina" />         <tag k="opening_hours" v="Mo-Su 17:00-01:00" />         <tag k="phone" v="+49 (8142) 444393" />         <tag k="website" v="http://www.cantina-olching.com/" />         <tag k="wheelchair" v="no" />


title: PostGIS - incorrect result distance between two points

text: I am quite confused the way PostGIS is returning the result for this basic query. So I am trying to find the distance between two points 45.2714,71.2087 and 42.3739,66.39 and the distance should be ~505 km but when I ran the query (considering 1 degree = 111km)               select st_distance (st_setsrid(st_makepoint(71.2087,45.2714),4326),            st_setsrid(st_makepoint(66.39,42.3739),4326))*111 as d       I get the result as 624 km. Can anyone please tell what I am doing wrong here ? Thanks.


title: Measurements in QGIS 1.8.0

text: I am using qgis version 1.8.0. My project contains tiff and several "shape" layers. The CRS is WGS 84 (EPSG:4326). I am attempting to measure the areas of the layer polygons. However if I click on the symbol "Measure area", I am getting a message. It keeps telling me the measurements are problably wrong. (Testing this statement by measuring leads to enormous but wrong values) It suggests to define a coordinate system in project properties. I have been trying this serveral times - with no results. Furthermore I read "Working with Projections" and "Measuring" in the User guide (http://docs.qgis.org/user_guide/html/en/user_manual/introduction/general_tools.html and http://docs.qgis.org/user_guide/html/en/user_manual/working_with_projections/working_with_projections.html). Also after adjusting all other CRS settings the problem stays. What should I do ? Thank you for your advice


title: Modifying File GDB Feature Class Features Using QGIS

text: I've installed the File GDB API using the OSGEO4W advanced installation when installing QGIS. I've had a bit of trouble with it initially and have posted a few questions with some helpful Stack members. I can now open File geodatabase feature classes in QGIS and start an edit session, but the only thing I can do is create new features using the 'Add Feature' tool. I cannot edit existing features because the 'Move Features', 'Node Tools', 'Delete Selected Features' and the 'Cut, Copy and Paste Features' Tools are all greyed out. I've tested this using both ArcGIS 10.0 and 10.1 file geodatabases with the same result. The tools are greyed out except the 'Add Feature' tool even though the layer is selected and active. I tested to see if I can get the tools to be active by using an ESRI Shapefile and this works which leads me to believe that the API isn't working as I think it should.....or maybe I'm just being a bonehead. Has anyone else experienced this behavior? Should I be able to modify existing features or should I only be able to create new ones? Thanks, Mike


title: Proper way to rbind SpatialPolygonsDataFrames with identical polygon IDs?

text: What's the proper R idiom for rbinding SPDFs together when the IDs overlap? Note that here (as is often the case) the IDs are basically meaningless so it's pretty annoying that I can't just make rbind ignore them....               library(sp)     library(UScensus2000)     library(UScensus2000tract)          data(state) # for state names     states <- gsub( " ", "_", tolower(state.name) )     datanames <- paste(states,"tract", sep=".")     data( list=datanames )     lst <- lapply(datanames,get)          nation <- do.call( rbind, lst )     Error in validObject(res) :        invalid class “SpatialPolygons” object: non-unique Polygons ID slot values          # This non-exported function designed to solve this doesn't seem to work any more.     d <- sp:::makeUniqueIDs( list(arizona.tract,delaware.tract) )     Error in slot(i, "ID") :        no slot of name "ID" for this object of class "SpatialPolygonsDataFrame"


title: Database Connections folder is empty when viewed with ArcGIS tool

text: I'm trying to create an ArcGIS tool which accepts an SDE connection file as an input. I made a new interface with an input parameter and set the type to Database Connections. When I open the tool and click the "browse" icon, the dialog opens to the Database Connections folder but there is nothing in it. If I navigate to Database Connections in the ArcCatalog window however I see all my usual connection files. Does anyone know why the tool would not be displaying those? ![enter image description here](http://i.stack.imgur.com/gIndA.png)


title: How to apply the same style to different columns of the attributes table?

text: How do I apply the same Categorized style, including same custom intervals and respective labels, to more than one column, if I have a shapefile with an extensive attributes table? If I save a style for a column, I can’t apply it to another one. Thank you in advance for you help


title: Edge column in pgr_dijkstra

text: I am using pgr_dijkstra for finding the shortest path between two given nodes. The network map is actually a grid based one and I have created all the necessary columns and working fine. But I do not understand how does the id2 (which is edge) given? ![Query in QGIS](http://i.stack.imgur.com/VGEyz.png) The values given for id2 (e.g., 59721224, 59721496 ...) can not be found any of my tables. So how does this id2 (edge) values are given? Every single time I query pgr_dijkstra, id2 values are changing but none of the other columns are changing. This is how my wp_norway_network1 table looks: ![wp_norway_network1](http://i.stack.imgur.com/6uB7U.png)


title: Why does MAPublisher need to install Office 2010 redistributable x64 version to connect to Personal Geodatabase (MS Access)?

text: I am trying to connect to a personal geodatabase. a .mdb file, with MAPublisher but I get the message immediately to install a Office 2010 redistributable x64 version. I have installed 64 bit office (Access) and the needed driver, but still get the message. Any ideas on this communication issue?


title: Is there a limit to the number of polygon points you can export to CSV in QGIS 2.4?

text: Using QGIS 2.4 I am currently importing a csv file with ~30,000 polygons (polygons in WKT format). I also have a United States shapefile. I have used the Clip function to get only the polygons that are in the US. When I go to save the new clipped shapefile as a CSV file with GEOMETRY set to AS_WKT, sometimes the number of lat/lon points is very long and so it cuts off, not finishing the polygon or multipolygon. Is there a way to set the export to include all points?


title: How do I retrieve start and end point coordinates with Python/arcpy?

text: How do I retrieve start and end point coordinates with Python/arcpy for a polyline feature class? I expect to pass a segment identifier to a subroutine and have it pass back start and end coordinates. The Field Calculator method doesn't work for me, because I need the value for other calculations that can't be performed within it. (I also prefer not to change the data to store these coordinates as attributes.) I am attempting to calculate address ranging for a "center out" addressing scheme. The address value depends on the distance to the "county center".


title: Configure Orfeo Toolbox to work with QGIS v2 Windows 7 from standalone installer

text: I have asked this question on the QGIS mailing list, but have not yet had a satisfactory result. The QGIS manual at http://www.qgis.org/en/docs/user_manual/processing/3rdParty.html is confusing about how OTB and QGIS work together and how to install them. On the one hand, it indicates that the OTB files are installed as part of the standalone installation package. But then goes on to tell us what the file paths (in the Options and Configuration area) for the OSGEO4 installer. My installation, on a Win 7 64 bit machine, was done via the standalone installer. I can see the various algorithms under Processing, but none of the OTB ones work. They all render an error that says that OTB is not configured correctly. Because I haven't installed via the OSGEO4 installer, I can't use the file paths that are outlined on the above page. All I can tell is that 'it should work' out of the box, but that's not my experience. Can anyone help?


title: QGIS on MAC cannot be opened because of an error

text: I know this question has been answered before, 3 months ago, but cannot work out how to join that thread. I have been using Mt Lion with QGIS successfully for some weeks but now get the 'Cannot be opened..' message on attempting to start. I have reinstalled the GDAL framework Complete but without success. Any suggestions most welcome! TW The core of the error message is: Application Specific Information: dyld: launch, loading dependent libraries Dyld Error Message: Library not loaded: /Library/Frameworks/PROJ.framework/Versions/4/PROJ Referenced from: /Applications/QGIS.app/Contents/MacOS/QGIS Reason: Incompatible library version: QGIS requires version 8.0.0 or later, but PROJ provides version 7.0.0