What's long binary data?

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?


How to add ECW support to GDAL bundled with PostGIS Windows

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


Looking for a template project for OSM data

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


Standards for creating a RESTful Geocoding webservice

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.


How to filter wrong GPS points?

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)


Enlarge/roughen the Vector Size of Tiles

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.


Last edited information missing when used QGIS - Dxf to Shp

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.


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

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.


QGIS identify features button issue

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


Map printing using OpenLayers

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.


How do you treat partial 3d features in PostGIS?

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?


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

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


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

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)


ArcGIS Linear Referencing successive overlay operations

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.


Adding georeferencing data into tiff image

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


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

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.


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

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?


Time not configuring correctly?

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.


Conditional Statement for attribute data

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?


Joining Multiple Features to One Feature

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?


Clipping between features within a single vector?

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!


How to automatically vectorize data from WMS?

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.


QGIS polygon areas inside a polygon

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.


How to align directions of connected lines?

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?


How to Interpolate very large datasets in QGIS

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?


Vector layer and google map do not line up in PDF

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.


FME incorrectly reading SHP file schema in reader

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?


How to install ArcGIS Code Assist in Aptana 3?

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?


FME assign specific node id with OSM writer

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.


wireless sensor network integration with GIS

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


Multiple labels in QGIS QML

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?


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

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


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

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.


Leaflet Custom menu links

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?


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

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


How to extract small portion of GIS database from entire database

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.


File GDB to SDE Database using FME

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.


Some data Loss when Import PGDB Table from Excel

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)


How to add a coordinate after add layer?

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


OpenLayers change feature's style

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'                 })             });


Labelling Kriging result (contour)

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)


Visualizing Layers with lyr-files in ArcGIS API for JavaScript

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?


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

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?


CartoDB st_makeline with geography datatype

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!


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

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 ?


Problem creating point shapefiles programmatically using python in QGIS 2.2

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!


Caching WMS data

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?


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

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


Registering local folders on ArcGIS 10.2.2

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.


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

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


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

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)


Using GeoServer REST Configuration to add a datastore

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.


Problem with deleting file locks in ArcGIS 10

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.


Trouble creating a query table

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


How to Buffer with intersections in QGIS

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)?


Where is the open Java client for ArcGIS Server?

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.


Hover/Popup feature problem

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";     }                    }


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

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?


WPF Runtime Arcgis - How to open MXD?

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


WPF Runtime Arcgis - How to create a MPK?

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


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

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!


Find gap in / endpoints of linestring

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)._


How to add a label of a point to a polygon

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.


Measurement units not returned in project CRS units

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?


FME autocad to csv

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?


All UK census geographies in one file irrespective of scale

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?


Myanmar Datum 2000 information

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?


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

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.


How to add topology to PostGIS 2.0?

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.


processing with xslt-

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" />


PostGIS - incorrect result distance between two points

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.


Measurements in QGIS 1.8.0

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


Modifying File GDB Feature Class Features Using QGIS

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


Proper way to rbind SpatialPolygonsDataFrames with identical polygon IDs?

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"


Database Connections folder is empty when viewed with ArcGIS tool

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)


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

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


Edge column in pgr_dijkstra

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)


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

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?


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

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?


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

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".


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

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?


QGIS on MAC cannot be opened because of an error

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


CartoDb - Aggregate and Infowindow

I'd like to use CartoDb to map historic addresses of a city (link to demo). Individual addresses are mapped as points on the map. Detailed information for each address point should be displayed in the infowindow The geocoded addresses are stored like this:               NAME   PROFESSION ADDRESS ... (+ CartoDb fields like cartodb_id, the_geom_webmercator)     Smith  shopkeeker grand-rue 10     Martin baker      grand-rue 11     Wagner grocery    grand-rue 11      As long as there is only one person living at an address, I get the desired infowindow by simply doing a _select_ * _from_ _table_name_. When two or more persons are associated with the same address I obviously would need some aggregate function in PostgreSQL. I am experimenting with:               SELECT  address, string_agg(name || ', ' || profession, '; ' ORDER BY name, profession) FROM table_name group by address      I am running into two problems:   1. In order to map, CartoDb expects the columns 'cartodb_id' and 'the_geo_webmercator' to be included in the SELECT statement. When doing so, the result of the query is not grouped by 'address' in the desired way; instead the unique 'cartodb_id' is used for grouping. So, the following is not grouping according to 'address':              SELECT  address, string_agg(name || ', ' || profession, '; ' ORDER BY name, profession), the_geo_webmercator, cartodb_id FROM table_name group by address, the_geo_webmercator, cartodb_id        2. The output of 'string_agg' is treated as one field in the design for the infowindow.              Grand-rue 11     Martin, baker; Wagner, grocery      What I would like to have, though, is a separate treatment of 'name' and 'profession' in the infowindow, like:               Grand-rue 11     1. Martin        baker     2. Wagner        grocery      Any hint is appreciated. Thanks! Peter


.dbf assocaited with aerial photos

I have aerial photos that are .jp2 format and are compiled into a .dbf. In ArcGIS I can simple add the .dbf file and all the aerials automatically appear. Can this be done in QGIS, rather than insert each individual .jp2 file?


Autodesk Infrastructure Map Server - WFS 1.1.0 GetFeature Request error

This is a duplicate request posted on the AutoDesk MapGuide Enterprise Discussion Group. If I take a working request...               <wfs:GetFeature service="WFS" version="1.1.0" maxFeatures="10"     xmlns:wfs="http://www.opengis.net/wfs">       <wfs:Query typeName="Water:WA_LARGE_VALVE" SRSNAME="EPSG:26910">         <Filter xmlns="http://www.opengis.net/ogc">           <Intersects xmlns="http://www.opengis.net/gml">             <PropertyName>Geometry</PropertyName>             <Envelope srsName="EPSG:26910">               <lowerCorner>491784.0 5456759.0</lowerCorner>               <upperCorner>491809.0 5456805.0</upperCorner>             </Envelope>           </Intersects>         </Filter>       </wfs:Query>     </wfs:GetFeature>      And replace Envelope with Polygon it fails...               <wfs:GetFeature service="WFS" version="1.1.0" maxFeatures="1"     xmlns:wfs="http://www.opengis.net/wfs">       <wfs:Query typeName="Water:WA_LARGE_VALVE" SRSNAME="EPSG:26910">         <Filter xmlns="http://www.opengis.net/ogc">             <Intersects xmlns="http://www.opengis.net/gml">             <PropertyName>Geometry</PropertyName>             <Polygon srsName="EPSG:26910">               <exterior>                 <LinearRing>                   <posList>489093.629 5456057.158 489092.66826402                   5456047.4034839 489089.822976626 5456038.02382838                   489085.202480615 5456029.37948835 489078.984339059                   5456021.80266094 489071.407511651 5456015.58451938                   489062.763171618 5456010.96402337 489053.383516101                   5456008.11873598 489043.629 5456007.158 489033.874483899                   5456008.11873598 489024.494828382 5456010.96402337                   489015.850488349 5456015.58451938 489008.273660941                   5456021.80266094 489002.055519385 5456029.37948835                   488997.435023374 5456038.02382838 488994.58973598 5456047.4034839                   488993.629 5456057.158 488994.58973598 5456066.9125161                   488997.435023374 5456076.29217162 489002.055519385                   5456084.93651165 489008.273660941 5456092.51333906                   489015.850488349 5456098.73148062 489024.494828382                   5456103.35197663 489033.874483899 5456106.19726402 489043.629                   5456107.158 489053.383516101 5456106.19726402 489062.763171618                   5456103.35197663 489071.407511651 5456098.73148062                   489078.984339059 5456092.51333906 489085.202480615                   5456084.93651165 489089.822976626 5456076.29217162                   489092.66826402 5456066.9125161 489093.629                   5456057.158</posList>                 </LinearRing>               </exterior>             </Polygon>             </Intersects>         </Filter>       </wfs:Query>     </wfs:GetFeature>      I receive is this error:               An exception occurred in FDO component.         String does not represent a valid filter      Is there something wrong with this request or maybe both requests?


Selecting lines occupying identical positions in MapInfo or QGIS

I have two shapefiles from a client, one is their entire network (complete with unique IDs) and another a selection I need. However the selection file does not have unique IDs but the entire one does. How would I be able to select the lines that appear in both files from the entire network shapefile? Thanks, Alan


Openlayers -- Is there a way to delay server requests on zoom for a particular layer?

I have a rather complex issue where I have 10 TMS Layers (I'm using OpenLayers) that animate. When I zoom in and out quickly I can quickly rack up 1000's of hits to the same server. Is there a way I can delay the request of that layer until zooming has stopped for x amount of time? Thanks!


Spatial Data with Unity

I was wondering if anyone has used spatial data as an input source to build 3D environments and interactive applications with Unity? The geospatial 3D world is very young, and some what immature, while the gaming industry is much more robust and developed. Unity is very developed and provides a lot of great options for building 3D environments. However, the typical data sources used in the geospatial world (vectors, imagery, DEMs) are foreign to Unity. I have seen some round about ways to use this data in a Unity project, but I was wondering if anyone has experience or success with this process, and would be willing to share. Thanks...


check if mapquest geocoding return no values

I use mapquest's geocoding api. I want to check if the request is sucessfull, but returned no data. I type `vdfsbvdf54vdfd` (just a silly string) in a form field and nothing happens. I expect an alert like so `Sorry, wrong input`. Alert never happens This is my code snippet               $.ajax({         type: "POST",         url: "http://open.mapquestapi.com/geocoding/v1/address?key=012",         data: { location: ad, maxResults: 1}         })               .done(function( response ) {             var geoclng = response.results[0].locations[0].latLng.lng;             var geoclat = response.results[0].locations[0].latLng.lat;                  if (geoclng=="") {alert("Sorry, wrong input");}                   //create magic using lon/lat                               )}      If it helps, when I alert the `response` I get `object Object`. I dont know if its the API or I handle wrong the request's response Thanks in advance


Refering to 0° of longitude/latitude in a range

This is more of a question in diction... In scientific writing, is there a preferable way to write a range, bounded by 0° longitude, in a way similar to "5° W to 10° E" if there is no cardinal direction for `0`. The two ways that seem most reasonable to me are "0 to 5 degrees longitude" or "0° longitude to 10° E". Thanks!


Rebuilding Indexes in SDE Database throws "Could not get a lock." error

While trying to do some maintenance on our SDE database but i am always getting a "Could not get a lock." error on the details section of the process. Do i need to setup something for me to successfully run the index rebuilding tool? Btw, i am using ArcMap 10.1 to rebuild the indexes. Thanks.


Make opengeo suit online,Need some free web server

I have installed opengeo suit and it works fine on localhost,How can I make it online?Which kind of web hosting do I need?Do I need any server or a linux webhost is ok? Is there any documention that I can use it to learn how to do this? Do you know some free web servers to test opengeo on it? Thanks a lot for your help


Software needed to cross-reference geocodes to a Shapefile

I do not have any GIS software except Microsoft MapPoint. I have State district KML and Shapefiles of the districts. I am assuming these are for ArcGIS. I do not want to draw maps. Instead, I am looking for a program that, if I give it one of 500,000 geocodes, it will tell me from a database perspective which District it is in. I am looking for analytical information not visual.


Flow Mapping Lines by Magnitude

![enter image description here](http://i.stack.imgur.com/dZ7V3.png)I am attempting to map the flow of remittances into Japan using the FlowMapper plug in in QGIS. When I use the Style-->Advanced--->Size Scale field--->magnitude command the lines become so thick and varied that they are impossible to distinguish from one another (in fact my map becomes one giant block of "lines"). Does anyone have any advice on how to fix this other than manually changing the size of lines classified under equal intervals? Is there a way to label the flow lines by country of origin so that the text is in and oriented with the flow line, not located on top (that is, once I'm able to get the flow lines to their desired thickness? How do I make the lines curved so they are more aesthetically pleasing? Thank you!


Cannot merge my map on CartoDB

I am trying to create a map on CartoDB using this kml file   http://cdn.journalism.cuny.edu/mapFiles/NYC_Community_Districts.kml and a CSV file. I have uploaded both to CartoDB but there the name column comes up empty so I cannot merge the two. The CSV file contains a community district column that I want to merge the map with but can't. Any ideas?


Batch merge large numbers of.asc grid files

I have >2.5k .asc grid files which I need to merge into a single overarching dataset. I would like to avoid having to convert them into esri grids first so was looking to use QGIS (v.2.0). I am aware of the command line length limitation with the QGIS merge tool so was looking to use the OSGeo4W shell to automate the process. Can someone help me with some code to cycle through the files and merge them? I was looking at this thread but the suggested solution doesn't work - it results in an "index out of range error" The files are all of the same format, in a single folder and have the name convention SDXXX.asc - where X is a single digit.


arcpy - ascii encoding error when updating field content in SQL Server table

I know there are plenty of posts and blog explanations about unicode errors, but I still can't figure out to handle it in my particular case. So here is my problem: I am writing a Python script to update records in an SQL Server table, with field content from a shapefile that has been edited in ArcPad. I use an `arcpy.da.UpdateCursor` to update existing records. I've put `# -*- coding: cp1252 -*-` on top of my script as I work with French characters. Just doing so, I get the ascii encoding error: > 'ascii' codec can't encode characters in position 0-1: ordinal not in > range(128) I've tried putting `u` before the field content:               ...     with arcpy.da.UpdateCursor(DBtable, DBFields, where_clause) as DBCur:         for DBrow in DBCur:             ...             DBrow[8] = u"{} - {}: {}".format(DBrow[8], date , AProw[8]) # AProw comes from a SearchCursor reading the shapefile.      I can then print a message with the text content, and it doesn't return an error message, but all characters are replaced (I get something like `???4???5????????>???5???` in my database field). Any accentuated or punctuation character can be present in the fields (there's a free text comment field) so I don't want to check for every possible non ascii character and replace it. I work with ArcGIS 10.2.2, the shapefile is edited in ArcPad 10.2 and the database is SQL Server 2008 R2. What am I missing? Thanks for any help. **EDIT: This only occurs if the destination database is SQL Server. No problem with a file gdb.** I have to add that the SQL Server table already contains non ascii characters.


How to replicate ArcGIS Intersect in PostGIS

I'm trying to replicate this ArcGIS process in PostGIS: http://blogs.esri.com/esri/arcgis/2012/11/13/spaghetti_and_meatballs/. This describes how to break buffered points down into polygons based on their intersections, counting the number of layers, and attributing that to the polygons in order to classify them. I'm using it to create a rough point density map with vectors, and the results were surprisingly nice for my data set in ArcGIS. However, I am struggling to come up with something workable in PostGIS where I need it for producing dynamic point density layers for a web map. In ArcGIS, I simply ran the Intersect tool on my buffered points layer to create the shapes I needed. In PostGIS, I ran this query: `CREATE TABLE buffer_table AS SELECT a.gid AS gid, ST_Buffer(a.geo,.003) AS geo FROM public.pointTable a;` `CREATE TABLE intersections AS SELECT a.gid AS gid_a, b.gid AS gid_b, ST_Intersection(a.geo,b.geo) AS geo FROM public.pointTable a, public.pointTable b WHERE ST_Intersects(a.geo, b.geo) AND a.gid < b.gid;` `DELETE FROM intersections WHERE id_a = id_b;` The output looks pretty much identical to the ArcGIS output, except that it is not breaking the polygons down to the same extent that is required for a meaningful density map. Here are screenshots of what I mean: ![ArcGIS](http://i.imgur.com/JDY0KO5.jpg) ![PostGIS](http://i.imgur.com/MBRj6dP.jpg) ArcGIS is on the left, and PostGIS is on the right. It is slightly difficult to tell, but the ArcGIS image shows the 'interior' polygon created where all 3 buffers intersect. The PostGIS output, on the other hand, does not create that interior polygon and instead it keeps its components intact. This makes it impossible to provide a classification for just that interior area with 3 layers on top of each other compared to just 1 for the outer parts. Does anyone know of any PostGIS function to break the polygon down to the extent I need? Alternatively, does anyone know of a better way to produce a point density map with vectors in PostGIS?


Transfering a lot of features from geodatabases

I'm having a problem where I need to export a large amount of content from 2 geodatabases of different nature (SDE -> File). I have three options: Spatial Select, a regular SQL select and the full feature class/table. What are the best practices to handle these large transfers? I'm having issues with the following code:                   public void ExportAll(IGisEntity entity)         {             ITable table = _Context.GetEntityTable(entity.EntityName);             IFeatureWorkspace workspace = _OutWorkspace as IFeatureWorkspace;             ITable outTable = workspace.OpenTable(entity.EntityFeatureClass);                  if (table == null)                 throw new ArgumentException("A tabela de origem não existe.");                  if (outTable == null)                 throw new ArgumentException("A tabela de destino não existe.");                  ICursor selectCursor = table.Search(null, false);             IRow selectRow = null;                  ICursor insertCursor = outTable.Insert(true);             IRowBuffer buffer = outTable.CreateRowBuffer();                  int recordCount = table.RowCount(null);             int count = 0;             int step = recordCount;             if (recordCount > 10000)                 step = Convert.ToInt32(Math.Round(recordCount * .1, 0));                  int remainder = 0;                  while ((selectRow = selectCursor.NextRow()) != null)             {                 buffer = selectRow;                 insertCursor.InsertRow(buffer);                 count++;                 Math.DivRem(count, step, out remainder);                 if (remainder == 0)                     insertCursor.Flush();             }             insertCursor.Flush();                  CustomComReleaser.ReleaseComObject(insertCursor);             CustomComReleaser.ReleaseComObject(selectCursor);         }              public void ExportSelect(IGisEntity entity, int sectorNo)         {             ITable table = _Context.GetEntityTable(entity.EntityName);             IFeatureWorkspace outWork = _OutWorkspace as IFeatureWorkspace;             ITable outTable = outWork.OpenTable(entity.EntityFeatureClass);                  if (table == null)                 throw new ArgumentException("A tabela de origem não existe.");                  if (outTable == null)                 throw new ArgumentException("A tabela de destino não existe.");                  IQueryFilter filter = new QueryFilterClass();             filter.WhereClause = String.Format(@"{0}_GR_IDENTIFICADOR LIKE '1{1}%'", entity.EntityCode.Substring(1, 3), sectorNo.ToString());                  ICursor selectCursor = table.Search(filter, false);             IRow selectRow = null;                  IRowBuffer buffer = outTable.CreateRowBuffer();             ICursor insertCursor = outTable.Insert(true);                  int recordCount = table.RowCount(filter);             int count = 0;             int step = recordCount;             if (recordCount > 10000)                 step = Convert.ToInt32(Math.Round(recordCount * .1,0));                  int remainder = 0;                  while ((selectRow = selectCursor.NextRow()) != null)             {                 buffer = selectRow;                 insertCursor.InsertRow(buffer);                 count++;                 Math.DivRem(count, step, out remainder);                 if (remainder == 0)                     insertCursor.Flush();             }             insertCursor.Flush();                  CustomComReleaser.ReleaseComObject(insertCursor);             CustomComReleaser.ReleaseComObject(selectCursor);         }      When one operating on ArcSDE systems these methods tend to hang and gives me I/O errors. I look into IGeodatabaseBridge interface and IGeoDbDataTransfer, but none of them gives me the option to export a cursor (product of a select, for example). **This is now a major issue: everytime I'm getting a RPC_ERROR_SERVERFAULT with this snippet.**                       while ((row = select.NextRow()) != null)             {                 buffer = row as IRowBuffer;                 insert.InsertRow(buffer);                 count++;                 if (count == 4999)                 {                     insert.Flush(); // this is where the error explodes at my face!                     count = 0;                 }             }      Any choices?


Get features in order they are selected using QGIS

Is there a way of establishing the order that a user selected features using QGIS and Python? By example: we have a network of line features representing a road network, and the user defines a route from point A to point B. Currently, using layer.selectedFeatures() returns the correct features, but in the order they were drawn rather than the order they were selected. I've tried a few things to do this manually:   * QgsMapToolEmitPoint(canvas) obviously doesn't work - I want the feature not the point clicked. But I guess in worst case I could hack a spatial query to do a look up of the closest feature. Not ideal.   * QgsMapToolSelect(canvas) doesn't seem to be available as part of the Python accessible GUI or CORE, although I can find mention of it existing elsewhere. Can anyone suggest a best case solution for this? All suggestions/tutorials/related questions appreciated. Thanks in advance for your time, Dan


With just a web browser we want our field staff to collaboratively make complex maps with secured private data--looking for options?

**Simple Question** With just a web browser we want our field staff to collaboratively make complex maps with secured private data. What are our options? Note: Please provide specific product examples to help me and others start the research process and ensure we dont overlook a great solution. I will give a concrete example in response to the comments ( @brad and @iant ) that my earlier version of the question above is too general. **Concrete Real-World Scenario Example** For example, consider there are 100 local people spread around a province in Afghanistan and they only have a web browser. The objective is to have them collaboratively making a map of the safety and danger zones in their province. **Requirements** They do NOT require advanced spatial analysis such a routing etc. Nevertheless, the final map-data product needs to be structured (geoJSON, shapefile, etc) so some one can later use the data for more complex spatial analysis. The key functional requirement is that these people need to upload to the collaboratively edited base map the following new features with associated labels:   * point locations of violent incidents with a labels   * polygon zones of dangerous and safe zones with labels   * line roads with labels  **More elaborate question** What are the specific products I should look into where the data will be secure, away from the public eyes --- _ideally_ something that connects with Google's Datastore, which is approved for SBU data (sensitive but unclassified data) --- _most likely_ some product where spatial data is stored and served internally ?


What is the difference between a Custom Layer & a plugin data source?

You can add custom data to ArcGIS desktop in two different ways:   * **Plug-in Data Sources**   * **Custom Layer** What is the difference between these two? I'm looking at the possibility of creating an editable layer from some custom format. I already know that it is not possible to create an editable Plug-in Data Source. Is an Editable custom Layer possible? Or is there some other way of achieving what I want?


How to draw ploygon in leaflet that use event?

Is there only one way to draw a polygon through the Leaflet.draw plugin? If you want draw anything, you must add the control ("map.addControl(drawControl);"). I want to draw a polygon on leaflet map by Event that I can trigger through any DOM event. Does anybody have a good idea how I can do this? Thank you for your time.


Simplify shapefile parcel layer

I have a very large shapefile with over a million parcels within it. The problem I am having is that the polygons representing each parcel are very jagged. What I am looking for is a simple and fast way of essentially "simplifying" these polygons and thus getting cleaner lines out of things. Naturally, in doing so I need to ensure all parcels are still linked together so there are no spaces between the parcel boundaries. What's the best way of achieving these objectives?


How to highlight and zoom to features after applying a definition query programmatically in Arcgis silverlight api

A definition query is applied to a bunch of layers after a user queries for a feature( with an ID). The def query is applied but I'm looking for a way to highlight those features that are a result of the query and zoom to them. Any ideas?


WGS72 to NAD83- how to find correct transformation

I am trying to get this raster file to convert from WGS72 to NAD83, but ArcMap10 doesn't appear to have the correct transformation listed. Where can I find this transformation, and/or how do I import ones that are not already listed as available?


How can I get my GPS track data into a graphics program?

I'm using a Garmin Edge 705 bike GPS/computer unit, which produces .gpx files. I can use the standard software it comes with to plot on Google Maps etc, but I'd like to be able to output just a vector polyline for use in Adobe Illustrator. Any ideas? I'm guessing .svg is going to be the best format, and I could use Inkscape or similar if that makes it any easier. * * * Sorry, other projects are taking up my time but I will try the suggestions out as soon as I can.


Incremental spatial clustering algorithm

I'm in search of a incremental spatial clustering algorithm. Here is my use case:   * users create entries with an initial position   * users may change the positions of existing entries I now want to implement a decoupled service that provides clustering information of these data. The service will be notified each time a new entry has been added or an existing entry has been moved. What is a good clustering algorithm therefore? Ideally, it should scale well up to high amounts of data and if there's a trade-off between cluster quality and runtime complexity, I'm okay with degrading results, and _eventually consistent_ results (stale results are okay for some time). To summarize my requirements:   * spatial clustering based on positions   * incremental modifications on changes   * add new positions    * change existing positions   * good runtime performance  Thanks in advance!


Need symbology that allows for symbols based on multiple fields, sometimes simultaneously (for certain features)

I've got a shapefile containing a large number of O&G parcel features. Each feature parcel has an associated status (things like "Held by Production," "Shut-in," "Leased, Non-producing") and an associated lease expiration date. I need to create a map that depicts this information by color-coding each parcel by its expiration date AND its status. For example, I need to color all parcels with leases expiring in 2014 red, 2015 orange, 2016 yellow, etc. I also need to identify those parcels with Status = "Paid Up" with diagonal lines overlaying the color fill designating expiration year. I hope my question is clear, and thanks in advance for any feedback you can offer! EDIT Throughout the process of working with this problem, I've come up with another tripping point: How does one assign symbology based only on the YEAR of a date field? For instance, my shapefile has a DATE field with the following format: 01/01/2014. I would like to symbolize each feature based on the YEAR only. The symbology tool allows me to assign symbology for each individual value (1/1/2014, 1/2/2014, etc.), but not based on only the YEAR portion of the DATE field. In fact, when I populate the symbology tool with all unique values of the DATE field and sort, the list sorts alphabetically rather than chronologically. This leads to weirdness like 1/1/2015 appearing before 1/2/2014.