title: LSA-SAF data reprojection with GDAL and R

text: The Land Surface Analysis Satellite Applications Facility (LSA SAF) distributes several products derived from the EUMETSAT satellites data. I am working with `R` to read and process the data (the file of this example is available from the Static Auxiliary Data page of the LSA-SAF service.)               library(rgdal)     library(raster)          r <- raster("HDF5_LSASAF_USGS_DEM_Euro.hdf")      Although the file is read correctly, the projection info is missing.               > r     class       : RasterLayer      dimensions  : 651, 1701, 1107351  (nrow, ncol, ncell)     resolution  : 1, 1  (x, y)     extent      : 0, 1701, 0, 651  (xmin, xmax, ymin, ymax)     coord. ref. : NA      data source : /home/oscar/temp/HDF5_LSASAF_USGS_DEM_Euro.hdf      names       : HDF5_LSASAF_USGS_DEM_Euro      values      : -32768, 32767  (min, max)      If I use `gdalsrsinfo` to parse the file and obtain the projection, I get `ERROR 1: ERROR - failed to load SRS definition from ...` which is more or less the same if I use `GDALspatialRef` in `R`:               > GDALSpatialRef('HDF5_LSASAF_USGS_DEM_Euro.hdf')     [1] ""      As far as I know these files use the GEOS projection. However, if I set this PROJ.4 string and the file is transformed to a latitude-longitude projection, the results are erroneous (look at the extent of the `projectExtent` output):               projection(r) <- CRS("+proj=geos +lon_0=0 +h=35785831 +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs ")          > projectExtent(r, CRS(projLL))     class       : RasterLayer      dimensions  : 651, 1701, 1107351  (nrow, ncol, ncell)     resolution  : 8.983153e-06, 9.043695e-06  (x, y)     extent      : 0, 0.01528034, 0, 0.005887445  (xmin, xmax, ymin, ymax)     coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0       It seems that the LSA-SAF data are not yet geo-corrected and the images are expressed in the "raw" Column-Row system, so this could be the root of the problem. There are Latitude and Longitude static files which could be used for remapping the rasters. However, I would like to use a more direct approach with a correct PROJ.4 string and a projection transform with `gdalwarp` or `rgdal::spTransform` in `R`. Thanks in advance for any help.


title: LayerUploadPanel not works

text: how can i add a layer upload panel tool in my suite sdk? The main problem is when i click "upload layers" button. Appear errors in the firebug. I don't know the reason why those errors appear TypeError:b[e] is not a constructor... TypeError:this.dataStore is undefined:this.dataStore.emptyText = this.dataStoreEmptyText; (in LayerUploadPanel.js line 327). i'm not sure if i haven't added the necessary code to app.js. i want to upload mi own shapes from layeruploadpanel, but i don't how to do it. Anyone could help me please?


title: How can I access a raster dataset's attribute table?

text: I need to loop through the rows in the raster attribute table and add a date for each row. Thanks,


title: contour lines from point shape attribute

text: I have a point shape file and one of the attribute refers to the altitude. Is is possible to obtain a contour lines shape file in GRASS?


title: Is a WMS map always a constant scale?

text: If I get a map back from a WMS server is the scale always constant over them map regardless of the reference system? So if I have a map which is 100 pixels wide, will the point at 50 pixels always be half way between the minimum and maximum values of the sides of the bounding box that defines the map? So can I always linearly interpolate the map coordinates based on the pixel position? Or can there be cases where 50% of the way across the map in pixels is not 50% of the way between the minimum and maximum of the width of the bounding box for some reason, like the projection is not constant?


title: how to intersect a raster subset with a buffered vector layer (points) in QGIS?

text: I am a bit new in the GIS and qGIs sphere so maybe my question will look a bit naive, but at the same time it means it should probably be easy for most of you guys ;-) I work in QGIS.I have 2 layers: (1) raster: grid representing global land cover / CRS "WGS 84" (2) vector: a series of point representing factories / CRS "WGS 84" I would like to select a subset from the grid only pixels with value 11 & 14 and intersect these selected pixels with a buffer of 50km around each factory. The result I would like to save as a new layer (preferably raster). How can I do that in a simple but correct way? Thanks Xavier


title: Traversing a PostGIS table tree

text: I'm developing a geospatial db for a company that operates at a number of large facilities. Its primary purpose is to help staff locate themselves and equipment within each facilities' multi-storey buildings. I'm using ArcGIS's Building Interior Space Data Model (BISDM) as a schema template and implementing it in PostGIS. The table structure looks something like this:               facility       building         floor           boundary           interior_space             equipment            transport         road         pedestrian_path            perimeter         fence         vehicle_gate         pedestrian_gate      Each parent table has a one-to-many relationship with it's children, so the schema defines a simple hierarchy. Each table has `short_name`, `long_name` and `geom` columns (as well as others, depending on the table). My problem is this: I want to build a reverse-geocoder API that takes lat/long as input and returns a string representing the address of the matching locations. I've gotten stuck on how to construct this string. My idea was to create a query that takes as its starting point the row containing the geom closest to the lat/long point and then traverse up the table hierarchy, collecting the value of short_name for each parent record until it reaches `facility`. So the assembled address in this case might be _'Room 323A, Floor 4, Building 12, Facility Q'_ (`interior_space->floor->building->facility`). However, because I don't know which table my query will start from (since it could be a point anywhere within the facility) I don't know what the parent table names will be in advance, and so I can't construct any `JOIN` queries. I'm familiar with approaches like _adjacency lists_ and _closure tables_ to model hierarchal systems, but these all require the data to be in a single table. The PostGIS requirement (well, strong recommendation) for different geometry types to be kept in separate tables means I can't use these approaches. Can anyone suggest how I can construct this string without having to create a bunch of hard-coded SQL queries?


title: How to put data in remote geoserver?

text: I created some layers and upload it to PostGIS and then GeoServer and then styled it using GeoExplorer. I have done all this in localhost. I used the iframe generated in GeoExplorer and used it in my webpage. So I can open it in my computer. But I don't understand how to keep this in remote server so that everyone can access it. I used OpenGeoSuite for doing all this. Is the data I uploaded in PostGIS and GeoServer just located in my computer or in remote server? Please help!! ## Question Extension This is the code that I am trying to run. But it's not working even in the same machine:               var map;          function init() {        map = new OpenLayers.Map('map_element', {});                  var road_name = new OpenLayers.Layer.WMS(                 "IIRS:river_canal_project - Tiled",     "http://localhost:8080/geoserver/IIRS/wms",                 {                     LAYERS: 'IIRS:river_canal_project',                     STYLES: '',                     format: format,                     tiled: true,                     tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                 },                 {                     buffer: 0,                     displayOutsideMaxExtent: true,                     isBaseLayer: false,                     yx : {'EPSG:4326' : true}                 }              );             Var industrial_park = new OpenLayers.Layer.WMS(                 "IIRS:future_industrial_project - Tiled",      "http://localhost:8080/geoserver/IIRS/wms",                 {                     LAYERS: 'IIRS:future_industrial_project',                     STYLES: '',                     format: format,                     tiled: true,                     tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                 },                 {                     buffer: 0,                     displayOutsideMaxExtent: true,                     isBaseLayer: false,                     yx : {'EPSG:4326' : true}                 }              );             var land_use = new OpenLayers.Layer.WMS(                 "http://localhost:8080/geoserver/IIRS/wms",                 {                     LAYERS: 'IIRS:current_landuse_project',                     STYLES: '',                     format: format                 },                 {                    singleTile: true,                     ratio: 1,                     isBaseLayer: true,                    yx : {'EPSG:4326' : true}                 }              );        map.addLayers([land_use, industrial_park, road_name]);        if(!map.getCenter()){            map.zoomToMaxExtent();        }     }      Along with the code, I inserted openlayers.js in script and only that. I found myself some mistakes and now I that can load the map but..the images are not loading. This is my edited code :               map = new OpenLayers.Map('map_element', {});                             // setup tiled layer                 var tiled = new OpenLayers.Layer.WMS(                     "earth:ne_10m_railroads - Tiled",           "http://localhost:8080/geoserver/earth/wms",                     {                         LAYERS: 'earth:ne_10m_railroads',                         STYLES: '',                         format: OpenLayers.Format.WMSGetFeatureInfo,                         tiled: true,                         tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                     },                     {                         buffer: 0,                         displayOutsideMaxExtent: true,                         isBaseLayer: true,                         yx : {'EPSG:900913' : false}                     }                  );                           map.addLayer(tiled);     var bounds = new OpenLayers.Bounds(                     -135.326858520508, -30.5615100860596,                     179.357788085938, 69.604377746582     );                   map.zoomToExtent(bounds);     }


title: Qgis - Raster clip Problem

text: I tried to use the Raster (Extration) tool "clipper" and I always get this error log: What is my problem? Fehler bei der Ausführung von Python-Code:               Traceback (most recent call last):       File "C:/Dokumente und Einstellungen/user/.qgis//python/plugins\GdalTools\GdalTools.py", line 342, in doClipper         self.runToolDialog( d )       File "C:/Dokumente und Einstellungen/user/.qgis//python/plugins\GdalTools\GdalTools.py", line 370, in runToolDialog         dlg.show_()       File "C:/Dokumente und Einstellungen/user/.qgis//python/plugins\GdalTools\tools\doClipper.py", line 47, in show_         BasePluginWidget.show_(self)       File "C:/Dokumente und Einstellungen/user/.qgis//python/plugins\GdalTools\tools\widgetPluginBase.py", line 44, in show_         self.onLayersChanged()       File "C:/Dokumente und Einstellungen/user/.qgis//python/plugins\GdalTools\tools\doClipper.py", line 71, in onLayersChanged         self.inSelector.setLayers( Utils.LayerRegistry.instance().getRasterLayers() )       File "C:/Dokumente und Einstellungen/user/.qgis//python/plugins\GdalTools\tools\GdalTools_utils.py", line 161, in getRasterLayers         return filter( self.isRaster, LayerRegistry.layers )       File "C:/Dokumente und Einstellungen/user/.qgis//python/plugins\GdalTools\tools\GdalTools_utils.py", line 156, in isRaster         if layer.usesProvider() and layer.providerKey() != 'gdal':     AttributeError: 'QgsRasterLayer' object has no attribute 'usesProvider'      Python-Version: 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)] QGIS-Version: 1.9.0-Master Master, 298f0c0


title: How to obtain to Max zoom levels for Google layer using Open layers

text: I was trying to obtain maximum zoom level ( level 22) for Google street map using Openlayers.layer.google script. The following code does not seem to work.It works to the maximum zoom level when Var zoom is set to 15. below or above 15, it set centre to world extent. Any idea what went wrong ? we tried using Google maps API v2 and v3 , same result for Maxm zoom, except for value = 15. var position = new OpenLayers.LonLat(lon, lat).transform(fromProjection, toProjection); var zoom = 22; map.setCenter(position, zoom);


title: How to Install the ArcGIS JavaScript 3.2 API library on an IIS Web Server?

text: I would like to deploy Server API for JavaScript 3.2 in my server to use it in my asp.net project instead of connect to server What are the steps to be followed in deploying it to a local WebServer?


title: How could I compare polygons with appropriate circle?

text: My problem is that I have a polygon shape file with 1034 records which are the Greek municipalities. My goal is to compare the shape of every record with the circle. In other words, I want to observe the deviation of the polygons from the perfect shape-circle. Thank you!


title: Using ArcGIS API for Javascript with Chrome on info window popup calls back two instances?

text: Using the ArcGIS API and harvesting external web map services. When using the application in Chrome when I click on a feature the info window calls back 2 instances of the same layer feature info in the popup. I just wondered if this was a known issue and if there's a way to remedy it?


title: FME API Module null?

text: I try to convert my geodatabase to Excel. And i have a strange message: > FME API version of module 'XLS_ADO' matches current internal version (3.6 > 20101014) Trying to find a DYNAMIC plugin for reader named `NULL' Loaded > module 'NULL' from file 'C:\apps\FME\plugins/NULL.dll' FME API version of > module 'NULL' matches current internal version (3.6 20101014) FME API > version of module 'XLS_ADO' matches current internal version (3.6 20101014) > FME API version of module 'XLS_ADO' matches current internal version (3.6 > 20101014) MS Excel Writer: Opening`null' for write operation MS Excel > Writer: Connection failed. Connection string > `Provider=Microsoft.Jet.OLEDB.4.0;Data Source=null;Extended > Properties="Excel 8.0;HDR=YES"'. Provider error`La table externe n'est pas > dans le format attendu.' Failed to open writer The table is not in the expected format. An idea ? Thanks


title: Can WMS be used to get a local copy of raster data?

text: I have tried reading up on WMS on wiki and other sources but haven't managed to get my head around a couple of points. 1) Using WMS is it possible to get a local copy of the raster data (say a tiff image tile). The reason, I ask is, I need to use his raster data in a secondary process that needs a local copy of the file. 2) I am primarily interested in allowing users the option of downloading orthoimagery and elevation from the USGS repository. I've tried using the Seamless viewer but I am thinking more along the lines of user specifying a coordinate (or four corner coords) and getting the orthoimage and elevation data corresponding to that. I thought WMS (Geoserver being my first option) is the way to go for this, but not certain. So wanted to get your advice. Thanks in advance. Sincerely, ST.


title: Error: Algorithm not found

text: I wrote the following code to start using QGIS APIs. in a very short script seen below, the " **Error: Algorithm not found** " is thrown and " **print processing.alglist()** " is showing none. Any idea why this is happening?               from PyQt4.QtCore import *     from PyQt4.QtGui import *     from qgis.core import *     from qgis.gui import *     import nltk     import sys, os     import shapely      import math     from collections import defaultdict     from shapely.geometry import LineString     import json     import datetime     import time     import itertools     import glob     import sys               app = QApplication(sys.argv)     qgis_prefix = "C:/Program Files/QGIS Valmiera/apps/qgis"     QgsApplication.setPrefixPath(qgis_prefix, True)     QgsApplication.initQgis()     layer1 = QgsVectorLayer("D:/Ehsan/Test/InputData/SHP/OSM_roads_no_psv_latin1_id.shp",     "test" , "ogr")          if not layer1.isValid():         print "Layer failed to load! :("     print layer1.isValid()          import processing     processing.alglist()     processing.runalg("qgis:clip","D:/Ehsan/Test/InputData/SHP/OSM_roads_no_psv_latin1_id.shp","D:/Ehsan/Test/InputData/temporary/tmpotile.shp","D:/Ehsan/Test/LM_Feat_test.shp")      Notice: the imports list is to show that I have almost everything imported and added to the path and pythonpath... Please help me to figure out how to resolve it. So far, the QGIS has seemed to be so buggy and problematic outside its python console. @Lukas Graf


title: time processing

text: I am computing NDVI for Landsat TM in grass 6.4.3 in raster map calculator, and time processing is many hours and finally tke message is "failed". trying to find a solution I discovered that it happens in all processes in raster map calculator. any advice? I have 3 images of landsat TM, TM3, TM4 and TM5. I started with GRASS 6.4.3 to correct these images as: Imagery→Satellite image tools→Atmospheric correction using also a DEM of srtm. The outcome seems excellent. Then, I tried to compute               NDVI=1.0 * float ((TM4-TM3))/float ((TM4+TM3))      and the process takes so many hours without results. I also tried to do the simplest TM4+TM3 and also the same thing.


title: SQL expression error for Select Layer By Attribute

text: arcpy.MakeFeatureLayer_management("N05-12_Station2_JGD.shp","station_12")          year=2012          arcpy.SelectLayerByAttribute_management("station_12","#",'"closing">=%d AND "opening"<=%d' % year)      > Traceback (most recent call last):File "", line 1, in > arcpy.SelectLayerByAttribute_management("station_12","#",'"closing">=%d AND > "opening"<=%d' % year) TypeError: not enough arguments for format string I want to select records that closing>=2012 and opening<=2012. Above script, it occurs an error. How can I write right script?


title: Arcgis Javascript code InfoWindow template

text:
I am creating an InfoWindow in javascript code like following.               var htmlForms = "City ID: ${CITYID}"+"
";     htmlForms += "City Name: ${CITYNAME}

";     htmlForms += "
";     htmlForms += "<button id='imgOut'>Hide</button><button id='imgIn'>Show</button>";     htmlForms += "
";          var formContents = {         title: "City Form",         content: htmlForms     };     var formTemplate = new InfoTemplate(formContents);      I wanna hide and Show image in this html script. So I added two button **imgOut** and **imgIn**. When I add dojo/fx wipeOut and wipeIn properties. But not working.                       coreFx.wipeIn({                 node: "cityImage"             }).play();


title: How send Geoserver WMS GetFeatureInfo query to an area

text: I am using GetFeatureInfo query to WMS service. When I click a WMS object on the map, the response is coming.               http://localhost:8080/geoserver/wms?     request=GetFeatureInfo     &service=WMS     &version=1.1.1     &layers=topp%3Astates     &styles=     &srs=EPSG%3A4326     &format=image%2Fpng     &bbox=-145.151041%2C21.73192%2C-57.154894%2C58.961059     &width=780     &height=330     &query_layers=topp%3Astates     &info_format=text%2Fhtml     &feature_count=50     &x=353     &y=145     &exceptions=application%2Fvnd.ogc.se_xml      But I want to send a request a polygonal area. Not a single x,y location. How can I do this? Is possible?


title: How can I create a "setback" buffer in QGIS?

text: I've been going a bit crazy in trying to create a "setback" buffer on QGIS for coastal lines in a map projecting this buffer outwards towards the sea. I've managed to create a normal Bidirectional buffer, but that's not the outcome that I'm looking for. Does anyone know how to go about this?? Kind Regards, Enrico


title: KML to Dat Converter?

text: Are there any free KML to DAT(destinator compatible) converter available ? I tried to use global mapper but was in vain. Also used ITN converter but the line was converted to points.


title: Digital Elevation Model Fusion

text: I want to fuse Digital elevation Models extracted from Carto-1 satellite. I already extracted DEMs & tried few algorithms in Matlab. Unfortunately as they are not Georeferenced i am not getting accurate results. So is their any tool or model which can use for Georeferencing & fusing DEMs?


title: How to use %value% (from For iterator) in SQL expression of Calculate Field in ModelBuilder?

text: I have a shapefile with data of fish specimens and the coordinates where they were tagged for research purposes. One of the fields in the attribute table of these shapefiles is the fish length in centimeters (ReLenCM). I am building a model that will create three new fields in the shapefile with weight values calculated with three different estimations of the length-weight relationship for that species. To make my enquire simpler, instead of using those functions yet I´m trying to just have the length value multiplied by an iteration. My For iterator is From 2, to 4, by 1. Meaning my three iterations are 2, 3 and 4. Therefore, for instance for the first row in the table (see pic) I am aiming to have something like this: ReLenCM = 79; W_2 = 158; W_3 = 237; W_4 = 316 This is the configuration of the **ArcGIS 10.1 for Desktop** model: ADD FIELD:   * Field Name_Add = W_%value% (resulting in W_2, W_3 and W_4)   * Field type = Double   * Input table = Model parameter (shp of tagged fish releases) CALCULATE FIELD   * FieldName = %Field Name_Add% (not entirely sure about this one)   * Expression = Here is where I want to multiply length by each of the iterations. You can see my attempt in the picture ([ReLenCM] * [W_%value%]), but with that expression all I get are "0" values... Basically, I just need to know how to insert the %value% variable from the iterator in the SQL expression. I have found some posts describing a similar issue, but always with String variables, not numeric. ![Model diagram](http://i.stack.imgur.com/a38sT.png)


title: Matching GPS tracks

text: I'm looking for a proper way to **process multiple GPS tracks** (so time ordered list of lat/lon points) and match/classify that segments, that are covered by multiple tracks: ![enter image description here](http://i.stack.imgur.com/WCygS.png) So it's not about recording, but using known algorithms/open source to do the analysis.   Sadly I didn't found a ready-to-use solution at the web, just scientific conceptual papers. My ideas would be:   * use spatial index to scan just in point subsets (but what with points at borders?)   * scan point sets with a buffer (to skip accuracy jitter)   * classify points that sit in the buffer But I like to have a real processing that is known to work. I guess this is a very common requirement to fleet monitoring solutions, so there must be some existing solutions that I just don't know :/   Can anyone suggest an existing tool (in QGIS or via a Python lib) that can do this matching? _BTW: Happy new year 2014 to all GIS and SE folks :)_


title: How to import label classes to make symbology classes?

text: In Answer to an earlier Question I was able to write SQL Expressions to create label classes. Now, I wish to use those label classes to create symbology classes. It is possible to make label classes from symbology classes, but I want to do it the other way around.


title: ArcPy Result Object Associated with Currently Running Script

text: We are at the end of a Python script using ArcPy. We can crawl around ArcPy to get information, but is there any way to get a result object for our currently running script? My motivation is to use existing code to process the same result object we might get from executing a separate script tool earlier in the script. To clarify, I want this object: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Result/000v000000n7000000/ I was wishing I could get a Result object for the currently running process. Here are my notes for accessing ArcPy instead:               #Result.inputCount and Result.outputCount     outParameters = []     inParameters = []     for parameter in arcpy.GetParameterInfo:         if parameter.direction == "Input":  #unsure of exact string             inParameters.append(parameter)         else:             outParameters.append(parameter          inputCount = len(inParameters)       outputCount = len(outParameters)           #Result.MaxSeverity     maxSeverity = arcpy.GetMaxSeverity()           #Result.MessageCount     messageCount = arcpy.GetMessageCount()


title: ArcObjects ISQLSyntax with Oracle DB esriSQL_DelimitedIdentifierPrefix issue

text: New to stack exchange, heading over from the esri forums, it is great! I am trying to produce a query that I can use across all workspaces, naturally using ISQLSyntax. The problem I'm having is that when using esriSQL_DelimitedIdentifierPrefix on a Oracle SDE instance it is returning a speech mark, whereas I'd expect it to return a blank string. When I use this code on a FileGDB or an MDB it works as expected, but on SDE I get the error: Underlying DBMS error[ORA-00904: "table.field_name": invalid identifier][OWNER.TABLE] The error goes away if I remove the prefix and suffix, so there are no typos in the query. I don't want to do a select case SDE because it seems against the point of having ISQLSyntax. Any Ideas?               Dim Prefix As string, Suffix As string     Dim upper As String, sQualColName As String     Prefix = pSQLSyntax.GetSpecialCharacter(esriSQLSpecialCharacters.esriSQL_DelimitedIdentifierPrefix)     Suffix = pSQLSyntax.GetSpecialCharacter(esriSQLSpecialCharacters.esriSQL_DelimitedIdentifierSuffix)          sQualColName = pSQLSyntax.QualifyColumnName("table", "field_name")          sVal = Prefix & sQualColName & Suffix & " = " & Chr(39) & Trim(strFindString).ToUpper & Chr(39)


title: Merge long/lat points within a bounding box depending on a radius in MySQL

text: Here is a pic of what I basically want to achieve: http://i.stack.imgur.com/gjbfa.png So as the title says, I want to merge long/lat points which they're radius (of 25Kms for example) touch inside a bounding box of long/lat points. Here is my very simple DB structure:               +-------+-------------+------------+     | id    |        long |        lat |     +-------+-------------+------------+     |     1 |   -90.27137 |   50.00702 |     |     2 |   -92.27137 |   52.00702 |     |     3 |   -87.27137 |   48.00702 |     |     4 |   -91.27137 |   51.00702 |     +-------+-------------+------------+      Here is my query so far:               set @bottom_lat = -100.27137;     set @bottom_lon = 40.00702;          set @top_lat = -80.27137 ;     set @top_lon = 60.00702 ;     ;          SELECT AVG(latitude), AVG(longitude)     FROM destination     WHERE latitude > @bottom_lat AND longitude > @bottom_lon AND latitude < @top_lat AND longitude < @top_lon      So my query just merging all points inside an imaginary bounding box without considering radius. I know that I would prorably have to work with the Haversine formula but I'm crap at maths and MySQL which make things a little bit difficult. Indeed, I could eventually merge points if I had just one radius but each points have its own radius and I'm struggling. This is for a student project and any help will be much much apreciated. References: -My query on SQL Fiddle: http://sqlfiddle.com/#!2/3a42b/2 ( contains a SQL Fiddle exemple for the Haversine Formula in comment ) -The Haversine formula in MySQL query: (work for checking all points inside a given radius)               SELECT*, ( 6371* acos( cos( radians(     @my_lat) ) * cos( radians(      destination.latitude ) ) * cos( radians(      destination.longitude ) - radians(     @my_lon) ) + sin( radians(     @my_lat) ) * sin( radians(      destination.latitude ) ) ) ) AS distance      FROM destination     ORDER BY distance limit 1     ;


title: Create ERDAS IMG Statistics using GDAL C++ API?

text: I want to create very many ERDAS .img files using GDAL under C++. I don't want to have to process them all in ArcMap to create statistics for them (mainly because ArcMap fumbles when calculating statistics for very many input images in batch mode). I want to create the .img files complete with statistics calculated for each. I see python utilities for this, but nothing under the C++ API. Does anyone have a solution?


title: What is difference between GeoTIFF and Esri ASCII grid file formats for SRTM data?

text: what is the difference between the ArcInfo ASCII format and the geotiff file format when downloading the SRTM data. I downloaded the tiles (21) for the Tibetan Plateau and imported thme into ArcMap 10. It looks quite different. The ArcInfo ASCII looks like a puzzle and the geotiff ones are very grey... but no puzzle. Has anyone an idea?


title: Using field calculator to delete specific contours

text: I have 3 contour shapefiles with 25m, 50m and 100m equidistance. Now I want to delete in the 25m-shp-file every **50m and **00m contour and in the 50m-shp- file every **00 contour. How can I show them using the field calculator of qgis? Thanks!


title: GeoServer and OpenLayers in production - Securing WMS calls?

text: I've a mapping site where users have to login in with username and password. Once logged in the user gets a map of their data using WMS and WFS to display and edit data. This data is served using GeoServer and displayed using OpenLayers. The WMS on the layer is secured using a username\password. I've seen similar questions to this but basically I want to access the secured layers via WMS without exposing my username\password. Any advice\clarity would be appreciated - my own thoughts are as follows:   1. User logs in to website.    2. WMS request using HTTPS for their specific information - uses some filter based upon who they are.   3. How do I send this filter without exposing it in the WMS call?   4. How do I send username\password without exposing it in WMS call? Theres lots of talk about reverse proxies - do I send the WMS request with no username\password then when the request reaches the proxy I validate the user then modify the request by adding the WMS username\password using the reverse proxy. If so dont I still need to send some identifier to determine who sent the request or will the proxy knows this due to the users session? Basically any help would be appreciated


title: How to store attributes of a large shapefile in a list?

text: I am trying to store the attributes of a shapefile into a list. With following code it works for smaller shapes:               addA = []     While provider.nextFeature(feat):         attrs = feat.attributeMap()         for (k,attr) in attrs.iteritems():             addA.append(attr.toString())      If I try to store a large shape my pc runs out of memory or freezes. I tried to put the attributes in a blist and in a numpy array, but nothing really worked (or probably I did not understand how to set up a proper numpy array...). Is there a solution for this problem?


title: Set a Required Field in QGIS Edit Form

text: Is it possible to make a field required in a QGIS Edit Form?


title: Which is better for tiling: raster2pgsql or GeoWebCache?

text: We have some raster image data stored as GeoTIFF which we would like to make available online as an overlay using the OpenGeo framework. As the data covers only a spatially restricted area and we don't need too much zoom levels, we consider precomputing ("seeding") the tiles as applicable with respect to memory consumption. Using OpenGeo, there seem to be two modules which can handle tiling:   1. raster2pgsql: a command line tool which converts a raster file into a series of SQL commands that can be loaded into postGIS. It offers the option to specify the tile size and the zoom levels (in their syntax "overviews")   2. GeoWebCache: assists GeoWebServer with cache functionality and tile generation. The latter can be done prior to client requests or on the fly. I'm not sure about how the two differ with respect to the precomputed tiles, can anybody clarify this? At first glance I tend to GeoWebCache because of the convenient access through its web interface.


title: Is it possible to put Geoserver behind a proxy and access the WMS/WFS using OpenLayers

text: I just start to work with Geoserver and Openlayers and one of the main concerns of the project now is do not expose the geoserver web interface outside of the internal network. No access to www.mydomain.com/geoserver/ However the application that is available in the Internet needs to get data (wms/wfs) from the Geoserver. This all through OpenLayers. Is this possible? Thank you all


title: Handle add new feature event and/or access feature before commit?

text: Is it possible with PyQGIS to catch "add new feature" event before commiting changes on the layer? or at least to access new features before commit.


title: ASTER mineral indices with VINR and SWIR data

text: I need to run a supervised classification on ASTER data in order to map a rock formation, which is composed of chert. There is no shortage of indices which can be used to highlight silica rich formations, see http://www.ga.gov.au/image_cache/GA7833.pdf. However, all of the indices I have found thus far utilize the TIR bands, which have a spatial resolution of 90m, and I need to map at a much finer scale than this permits. The either of the the 30m and 15m resolutions of SWIR and VNIR respectively should suffice. The problem is that I can't find any silica indices which use these bands. I understand that any such indices might be sub optimal in comparison to those which use the TIR bands, but in this case resolution is of uptmost importance. So, the question is; does anyone know of any established chert/silica indices which use the VNIR or SWIR bands?


title: What is best way to reorganize osm into PostGIS database entity?

text: I'm just wondering how to reorganize the OSM given entity to more flexible and schematic purpose. The OSM gave 4 entities the planet_osm_line, planet_osm_point, planet_osm_polygon and planet_osm_roads. Is there a good practice to reorganize them like into city_building, city_road and etc. I'll be reorganizing a city's spatial objects. Thank you in advance


title: Generating polygons with shared boundaries using QGIS?

text: Essentially the problem is that if you digitise each polygon separately the boundaries will not match exactly. My instinct is to have a table of line segments and generate the polygons from that using postgis but then all the line segment end points have to join exactly which is painful as there are a lot of line segments. I have thought about writing a script that uses heuristics to match up lines and then add extra segments to bridge the gaps as suggested in a previous question I asked: How to automatically convert line features into polygons using QGIS/PostGIS? I feel this is such a common task that there must be some "best practice". If I were starting from scratch then how would I do it?


title: Batch fix overlap errors

text: I have a feature class with about 250,000 contiguous polygons. ArcEditor turns up around 9,000 overlap errors when I validate the topology. I want to merge each of these errors into one of the polygons they source from - they're all very tiny overlaps, so it really doesn't matter which one. All 9,000 errors can be fixed in exactly the same way. In the Error Inspector, the Merge option is only available when I select one error at a time - selecting a bunch of them greys it out. This means I need to do the same repetitive task for every single error. I would love to be able to avoid doing these one at a time, might save me from carpal tunnel and wasting the better part of a day. Any suggestions?


title: Clipping points using buffer in Arcmap takes FOREVER

text: I'm clipping points (crimes) by using a buffer around a shapefile of city parks. There are quite a few records involved (2 million+), but it is still taking forever - as of now, the clip has been running for an hour and ten minutes. I'm also using a subset of the parks data, so I'm not clipping crimes with every polygon in the shapefile. My questions: How does Arcmap go about it's clipping? Does it have to check every point, even those that are nowhere near the polygons? Is there a better method? I have QGIS at home, but I'm not sure if it's as fast as Arcmap. Thanks, Ward


title: How to hide infoWindow while using measure tool and showing it back after using it ArcGis JavaScript API

text: I have an infoWinfow set on the click event of the map. I want to disable this event when I'm using the measurement tool and then enable it back after the measurement tool is used. So far I've done this code, where the click event for the infoWindow is disable when the measurement tool is used. I want to enable the click event back after the measurement tool is used. Any ideas on how this could be done? Thanks.               function initToolbar(mymap) {     var measurement = new esri.dijit.Measurement({         map: mymap     }, dojo.byId('measurementDiv'));     measurement.startup();     measurement.setTool("area", false);          measurement.areaToggleButton = function () {         this.clearResult();         this.setTool("area");         dojo.disconnect(clickHandler);     };     measurement.distanceToggleButton= function () {         this.clearResult();         this.setTool("distance");         dojo.disconnect(clickHandler);     };     }


title: How to remove polygons in shapefile based on attributes in R

text: I have a shapefile where some attribute data is blank. I can delete those rows attribute table using na.omit. However, the polygons remain when I map the shapefile. I'm using R to make these changes. My shapefile contains over 900 records/polygons and 450 of these records have NO attribute information attached. These 450 are the ones that I want to remove. In the example, below I will need to remove both the shape and the record of ID:4.               For example:     ID Attribute1 Attribute2     1     data       data     2     data       data     3     data       data     4      NA         NA      How do I go about removing these shapes in R? Thanks for your help.


title: ArcGIS 10.1 SP 1 install on Windows Server 2012 questions

text: I will be setting up ArcGIS server 10.1 SP 1soon and I have a few questions before I start this. I have looked at the Web Adapter and the ArcGIS install guides but wanted to get some input. 1.Has anyone been installing on Windows Server 2012? I will be installing ArcGIS 10.1 SP1 and the Web Adapter on a separate 2012 servers. Have there been any problems with this that people have ran into? Should I just stick with 2008 R2 or just go with 2012? 2.I see the ArcGIS server needs IE 7 min, and looks like it can go up to 9. Will this software work with IE 10 though? 3.What ports am I going to need open between the Web Adapter server and the ArcGIS server, 6443? Is 6080 just for the server to World if I am doing the web and GIS all on 1 server? 4.If I am using SSL and HTTPS over port 6443, should I be using certs?


title: 3D elements buffer extrusion

text: I'm using ArcScene10.1 and I have a dataset containing of 3D points (x,y,z coordinates) with extrusion values (Evalue). Now I would like to create a square or ring buffer around the points and extrude that buffer by the extrusion value. The goal is to "layer 3D to feature class" those geometries and colour them by an other value (Cvalue). So far I was able to create a buffer around those points but I can't extrude them. eg.:               x y z Evalue Cvalue     1 2 3 5       3     2 4 1 3       15      ... ![draft](http://img594.imageshack.us/img594/741/extrudebufferdraft.jpg)


title: Get raster expression

text: In Arcmap I have made a new raster with ArcToolBox -> Spatial Analyst Tools -> Map Algebra -> raster calculator. Now, months later, I have forgotten the source formula expression. How can I display?


title: PostGIS Trigger doesn't work

text: The trigger below don't works and I need that when a new pipeline (adutoras) is inserted or updated, the procedure check what are the municipalities (municipios) that intersects the pipes and put the names of the municipalities inside the field "municipios_abastecidos".               CREATE OR REPLACE FUNCTION topology_fields()      RETURNS "trigger" AS     $$     BEGIN     NEW.municipios_abastecidos:=(                           SELECT string_agg(municipios.nome,', ')                           FROM municipios, adutoras                           WHERE ST_Intersects(adutoras.geom, municipios.geom)                           AND adutoras.gid = NEW.gid                         );     RETURN NEW;     END;      $$     LANGUAGE 'plpgsql';               CREATE TRIGGER fill_topology_fields     BEFORE INSERT OR UPDATE ON adutoras     FOR EACH ROW EXECUTE PROCEDURE topology_fields();      Thanks in Advance


title: leaflet fitBounds with multiple marker

text:
I have created two checkboxes. If I check one it will zoom to a point with marker. If I uncheck the checkbox the marker removed. Similarly second one is working. But I want to do if the two checkboxes are checked at a time it will zoom at a position where both are visible. Again when one of it checkbox is unchecked it will zoom to the points best fit position. I know that fitbounds can solve this but couldn't imliment it. Here is my code I used. Any suggestion will be a great help. Thanks               <input type="checkbox" name="One" id="one" >One</input>

<input type="checkbox" name="Two" id="two">Two</input>
           $('#one').click(function(){          if($("#one").prop('checked') === true) {              var marker = L.marker([32.304255,-64.752098]).addTo(map).bindPopup('One');              map.setView(new L.LatLng(32.304255,-64.752098), 16);              marker.on('dblclick', function(){              marker.openPopup();              });          }      });      $('#two').click(function(){          if (this.checked) {              var marker= L.marker([32.319343,-64.728066]).addTo(map).bindPopup('Two');              map.setView(new L.LatLng(32.319343,-64.728066), 16);              marker.on('dblclick', function(){              marker.openPopup();              });           }      });          var xLayer =[];     function clearLayer(){     for(var i =0; i < xLayer.length; i++)         map.removeLayer(xLayer[i]);         xLayer = [];     };          $('#one').click(function(){         if($("#one").prop('checked') === false) {         clearLayer();         }     });     $('#two').click(function(){         if($("#two").prop('checked') === false) {         clearLayer();         }     });


title: Creating seperate layers out of one big polygon shape by using attributes

text: I try to explain may problem. I have a big polygon shape (> 200 fields). I have the owners of these fields in a column. I need the seperated shapes for each owner. Is there a tool to do that? Is it also possible to get an information about the whole area? Like a summary: 201 polygones area: 1450 ha In Saga I have a tool "Split shapes by attribute". This works fine with the expected result (10 owners: 10 new layers). This tool doesnt work with the sextante plugin. Thanx!


title: field calculation in arcgis

text: * i have a data which contain a code of area an value of area in shapefile, ex.: area|value   1|33   1|43   1|31   2|12   2|43   3|66   3|76   the question is how to calculate total value per area with arcgis, i heve been tried with arcgis calculator and it doesnt work? what i expect for result is      area|value|total   1|33|107   1|43|107   1|31|107   2|12|55   2|43|55   3|66|142   3|76|142


title: Working with ambiguous AddGeometryColumn functions

text: I'm using PostGIS to add geometry columns to tables, but keep receiving the error:               ERROR: function addgeometrycolumn(character varying, character varying, integer, character varying, integer) is not unique     SQL state: 42725     Hint: Could not choose a best candidate function. You might need to add explicit type casts.      The statement I am trying to call is: SELECT addGeometryColumn('schema', 'table', 'geom', 3857, 'POINT', 2); I have tried type-casting each argument but I received the same message. I have also removed and re-installed the postgis extension, which made no difference. I have used this function countless times before and am somewhat puzzled as to why it isn't working all of a sudden. I have also tried a few other postgis functions which gave me similar errors. I am considering reinstalling postgresql/postgis but I will have to spend time reconfiguring the set-up so I am hoping there is an easier way to fix/get around this. Does anybody have any advice as to what may be causing this? Thanks in advance


title: Looking for PyQGIS Syntax Examples/Guide

text: I am a long-time user of arcgis with intermediate experience python scripting for automation of geoprocessing workflows. I am moving towards using QGIS on the regular and I would like to dive into python scripting for QGIS. However, I have previously relied on esri's help screens which contain python syntax as well as examples (ESRI Example). Is there an equivalent of this python documentation for QGIS? I have found the PyQgis Cookbook limited and the API Documentation beyond what I'm able to undertake. Any suggestions would be much appreciated. thanks in advance, mike


title: Census Demographic Data from DMA ID

text: How can I go about finding demographic data by DMA ID? Is there a data file I can download from somewhere that I can use to import into my web application?


title: Making Tiles - Varying the Zoom Level by Population Density

text: I want to create map tiles for all of the United States at high zoom level (12 or 13). I am using TileMill and mapping census tracts (income). However, I only need the high zoom level for populated areas. How could I exclude the areas with low population? I want to do this to minimize tile generation time and webhost disk space. For instance, I'm going to treat Alaska separately as most of it has low population density. I can also exclude oceans and large lakes. I could probably exclude a large portion of WY, ND, SD, and several other western states.


title: joining polygons

text: I am wondering how to join spatial polygons using R code ? I'm working with census data where certain areas change over time and I wish to join the polygons and the corresponding data and simply report on the joined areas. I am maintaining a list of polygons that have changes census to census and that I plan to merge. I'd like to use this list of area names as a look up list to apply to census data from different years. I'm wondering what r function to use to merge selected polygons and respective data. I ahve googled it but simply become confused by results. Hoping you all can help :)


title: Insert into X coordinate separately?

text: I have a postgis table with a point geometry for lat/long but i want to run a query that inserts a new row with lat, long values separately instead of concatenating them into a geometry. so is there a way to insert the X and Y values separately ?


title: Arcobjects: DisplayTransformation.FromMapPoint() transforms wrong

text: I want to transform map coords into page coords in ArcObjects in Python. So I do, say:               >>> pApp = NewObj(esriFrame.AppROT, esriFrame.IAppROT).Item(0)     >>> pDoc = pApp.Document     >>> pMxDoc = CType(pDoc, esriMapUI.IMxDocument)          >>> pDisplay = pMxDoc.ActiveView.ScreenDisplay # something wrong here?     >>> pDisTrans = pDisplay.DisplayTransformation          >>> pageX, pageY = pDisTrans.FromMapPoint(ptInRealCoords) # coords are ~389700,~316671m     >>> pageX     6523059 # which should be ~4.73 in cm      Any clues whats wrong?


title: Openlayers Linestring over edge of map

text: im trying to create a route displayer for ships. But im struggling with the edges of the map. When going from ie. 179* to 2* it displays a line across the entire map, instead of displaying it on the continuation of the map. (Map repeats) I do see why it would do that. But i was hoping there was some support in openlayers for making this more or less automatic instead of breaking up the route. Any ideas on how to fix this? **Expecation:** ![Expectation](http://i.stack.imgur.com/1ttwm.png) **Current display:** ![Current display](http://i.stack.imgur.com/PH5J9.png)


title: Find end of arc given start/end angles, start point and radius

text: Using PostGIS, how do I work out the end point of an arc given the start point, start and end bearings (relative to true north) and radius of the arc? Eg. say I need draw an arc of 4000m radius from bearing 180 degrees to 90 degrees with the start (the 180 degree point) at POINT(10 20). How do I calculate the end point (other than approximately drawing the arc, point by point)? This is only for short distances, so working on a plane should be OK. I can work out the distance between the start and end point using the length of chord formula and if I could just get the bearing between them as well then I could use ST_Project - but I don't know how to work out the bearing. (In the simple example above it should be 45 degrees.) **Solution:** It ended up being quite simple in the end. Thanks to **whuber** for explaining it very well. I wanted to put the PostGIS code here in case anyone else wants it:               centre_point := ST_Project(start_point, radius_m, normalise_angle(start_angle + pi() / 2 * CASE WHEN is_clockwise THEN 1 ELSE -1 END));     end_point := ST_Project(centre_point, radius_m, normalise_angle(end_angle + pi() / 2 * CASE WHEN is_clockwise THEN -1 ELSE 1 END));      `normalise_angle` is a simple function that ensures the angle is between 0 and 2PI (as required for ST_Project), `start_point` is a `geography(Point)` and `radius_m` is the radius in metres. Angles are in radians.


title: Create a color scale/legend for choropleth map using geotools (or other open source Java library)?

text: I'm experimenting with choropleth maps using geotools. Is there an easy way to add a color scale to the map (using geotools itself or some other os library)? By a color scale I mean a bar that shows how the feature attribute values map to color, with labeled ticks either automatically generated, or at break points of my choosing. Ideally, I'd like to generate a separate color scale that could be rendered independently of the choropleth map, so that I can use a single color with layouts containing multiple maps. I'm not absolutely wedded to geotools, and would consider other open source libraries, though I'd much prefer to stick with JVM languages (Clojure, Groovy, Java, ...).


title: Map source for iLife '11?

text: The new iLife has map functionality, what's the source of the map data used? ![ilife map image](http://i.stack.imgur.com/KoMRU.jpg) Further pictures on engadget, which show a city-level view.


title: Is there a geodetic benchmark that is fixed with respect to plate tectonics?

text: I see that NGS benchmarks are used to monitor the movement of geologic plates. Is one of the benchmarks assumed to be in a fixed location, with movements of all other benchmarks measured with respect to it, or are all benchmarks moving as a result of plate tectonics?


title: geopandas error when writing to_file() "ValueError: 'long' is not in list"

text: I'm having a problem getting the "to_file()" function to work with geopandas               import geopandas as gpd     f  = gpd.read_file('fin.shp')     f.to_file('fout.shp')      Error Output:                   ---------------------------------------------------------------------------      ValueError                                Traceback (most recent call last)     <ipython-input-14-6fdee04dae37> in <module>()           1 f  = gpd.read_file('fin.shp')     ----> 2 f.to_file('fout.shp')          C:\Python27\lib\site-packages\geopandas-0.1.0.dev-py2.7.egg\geopandas\geodataframe.pyc in to_file(self, filename, driver, **kwargs)         308         filename = os.path.abspath(os.path.expanduser(filename))         309         with fiona.open(filename, 'w', driver=driver, crs=self.crs,      --> 310                         schema=schema, **kwargs) as c:         311             for i, row in self.iterrows():         312                 c.write(feature(i, row))          C:\Python27\lib\site-packages\fiona\__init__.pyc in open(path, mode, driver, schema, crs, encoding, layer, vfs)         154         c = Collection(path, mode,          155                 crs=crs, driver=driver, schema=this_schema,     --> 156                 encoding=encoding, layer=layer, vsi=vsi, archive=archive)         157     else:         158         raise ValueError(          C:\Python27\lib\site-packages\fiona\collection.pyc in __init__(self, path, mode, driver, schema, crs, encoding, layer, vsi, archive, **kwargs)         129             self.encoding = encoding         130             self.session = WritingSession()     --> 131             self.session.start(self, **kwargs)         132             self.encoding = encoding or self.session.get_fileencoding().lower()         133           C:\Python27\lib\site-packages\fiona\ogrext.pyd in fiona.ogrext.WritingSession.start (src/fiona/ogrext.c:15304)()          ValueError: 'long' is not in list      I'm using fiona 1.1.3 and geopandas 0.1 with Python 2.7 on Windows. It may be a problem with my ogrext.pyd dll in fiona - any help is much appreciated


title: How to use ISymbol QueryBoundary

text: I have a balloon callout ISymbol and I want to get a polygon object so I can see if the user clicks it. I am trying to use the function QueryBoundary http://help.arcgis.com/en/sdk/10.0/arcobjects_net/componenthelp/index.html#//001w00000372000000 The function has four arguments:               public void QueryBoundary (         int hDC,         ITransformation displayTransform,         IGeometry Geometry,  //what can I put in here?         IPolygon boundary     );      One of my major hangups with the third argument is simply the fact that I am using this function to try to get the geometry... so what sense does it make to pass in the thing that I want? Also since the function is called like this :               thissymbol.QueryBoundary(screenDisplay.hDC, displayTransformation, thisGeometry, polygon);      Presumably it doesn't need the geometry of "thissymbol" since "thissymbol" is the caller.


title: which is the suitable Geographical Coordinate System and Projected Coordinate System for India

text: I want to create a spatial data for entire India, usually I use WGS84 as GCS and UTM zones for PCS in the places in India, but when I have to create the spatial data for entire India, I get confuse which geographical coordinate systems and projected coordinate system will be suitable for India, please help me to identify the suitable one and explain me why it is suitable.


title: How to return end of string up to a certain character using Field Calculator?

text: I am trying to accomplish this task: http://stackoverflow.com/questions/15851568/how-to-get-the-last-part-of-a- string-before-a-certain-character within the Field Calculator. I have the below but cannot get it to run. What am I missing? Any suggestions or assistance would be much appreciated. ![enter image description here](http://i.stack.imgur.com/sMxbT.jpg)


title: Obtaining a handle for newly created map in GeoExt

text: Hopefully a straightforward question... I am following the example based on this page here http://www.geoext.org/lib/GeoExt/widgets/MapPanel.html to create a map within a window, here's my code:               new Ext.Window({         title: "GeoExt MapPanel Window",         height: 400,         width: 600,         layout: "fit",         items: [{             xtype: "gx_mappanel",             id: "mappanel",             layers: [basemap],             extent: "-5,35,15,55"         }]     }).show();      So the above example create a new map in a floating window which is exactly what I'm after... My question is how exactly do I get a handle to the map within this new window so I can add controls and layers after the map has been created? Something along the lines of:               map = grab_the_map_somehow_from_Ext.Window_above();     map.addControl(my_new_control);     map.addLayers([my_new_layer]);      Thanks for your help.


title: Spatially join polygons (target) to polylines (join) on longest intersecting polyline

text: I am curious if I can, without having to write a python or an ArcObjects script, create a spatial join that will result in a polygon feature class containing the attributes of the LONGEST intersecting join polylines. The input polylines I have roughly correspond to the polygons to which I want to transfer the attributes. I've determined that it is almost always the longest intersecting polyline whose attribute I need to transfer into the polygon. Currently, using the figure below as an example, the spatial join hits 3 polylines (join_count = 3 for this particular polygon) but instead of joining on the longest intersecting polyline it seems to choose the first one it encounters. (which is also contrary to what the tool's description says about the JOIN_ONE_TO_ONE operation which is suppose to aggregate the join features into the output feature. The description speaks of SUM merge rule in its example which I do not see in this particular tool as an option) I am using INTERSECT because there is no suitable spatial relationship I can think of in this "polyline to polygon" situation. I understand what is going on but don't know whether a workaround exists. I just wish there would be another set sub-parameters allow me to control what happens when more than one join features exists; in this case selecting the longest polyline by default. I am using ArcGIS 10.2.1 Standard. ![enter image description here](http://i.stack.imgur.com/ESjIZ.jpg)


title: Display linestring on Cartodb with SQL

text: I have a Cartodb api for a separate map, I need to display a linestring for each time the sql is run, but not store the linestring in the database.


title: Which kind of grid format is this?

text: in a project I got some raster files with the following content (header and data are in one txt file):               HOUR 17     MINUTES 30     YEAR 2012     MONTH 5     DAY 4     CITY Emilia;Bologna     NROWS 99     NCOLS 156     ULXMAP 11.2087224106     ULYMAP 44.5539582648     XDIM 0.0016666668     YDIM 0.0012500001     CELLE_TRUE 2999     CELLE_FALSE 1     0.15478748  0.15478748  0.15478748  0.15478748  0.15478748       0.15478748  0.15478748  0.15478748  0.15478748  0.15478748       0.15478748  0.15478748  0.077375025 0.018038042 0.01576243     [...]      I tried to identify the format but did not get far. It seems to be a "kind of" ESRI ASCII grid, but it features non-square pixels (as you can see), which is not supported by ESRI ASCII grid (AFAIK). It might be an Arc Grid format (I found this link: http://landcover.usgs.gov/faq_importing_arc.php) but then I am lost. GDAL does not recognize the format, QGIS and ArcGIS 10 don't like the format either. Does anybody know this particular format (coming from an Italian source)? Or give me a hint how to process it to an more "readable" format? Regards and thanks in advance Jan


title: Accented characters in csv file appear as '?' when importing to QGIS

text: The tabular data (.csv file) and the data in the attributes table of the vector (.shp file) I am trying to join includes accented characters - they are place names in Pamana. For example 'Río Abajo' and 'Curundú'. The problem is that when I import the csv file as a vector file QGIS can't seem to read the accented characters and puts a little '?' instead for those characters. So I just get NULL values appearing in the joined attributes table - all the data for the unaccented names turns up perfectly. I have tried changing the names to unaccented latin characters, from the earlier example 'Río Abajo' becomes 'Rio Abajo'. But then the names don't match up with the place names already in the attributes table of the shapefile I am using (administrative areas for Panama) so the NULL values appear again. Hope you can help.


title: Fw tools using error

text: I am using Fw tools. And I am writting cmd screen:               ogr2ogr -f \"GeoJSON\" -t_srs WGS84  D:\\electric_Data1\\direk.json  D:\\electric_Data1\\ag_direk.TAB      tried. And I have a error. error code:               unable to find driver '"GeoJSON"'      The following drivers are avaible:           esri shape file     mapinfo file     bla bla bla     GeoJSON     bla bla bla


title: "Easy Custom Labeling" arrows: How do I get them?

text: I have installed the Easy Custom Labeling plugin, but haven't figured out how to make the "Generate lines" function work. When I click on the tool in the toolbar, it opens the dialogue to specify my CRS. Ok, that's done. So having moved a label that was too crowded, how do I make the connector line? THANKS TVZ (I'm running QGIS 1.8)


title: Save Multiple Selection from Attribute Table

text: I am working with a shapefile that has all the address points for a municipality - it has over 500,000 features. I would like to search the attribute table for specific addresses and save those somehow so I can sort and only show the ones I want - likely around 1000. After searching for one property I can move it to the top of the list, but when I search for the next and move it to the top, it replaces my first result. How can I flag or select a series of search results and then only save those in my attribute table? P.S. I am using QGIS v. 2.4 Thanks!


title: Class Probability Tool in ArcGIS Desktop 10.2

text: I am trying to create a multi-band raster using the class probability tool in ArGIS Desktop 10.2. I used the classified image I created based on my signature file as the input raster, and my training classes (.gsg) extention as my input signature file. All I get as an output when the tool runs is a mostly grey layer with no actual information attached to it. Anyone have any ideas about what I am doing wrong.


title: Any way to view/create/edit Google My Maps KML files on Android?

text: I have two large (several hundred features) maps created by Google My Maps that I use for my field research. They contain important information required by the field crew for their work on site. Currently these maps are accessed via iPhone by using a free app called My Maps Editor. This allows each individual to view all the features as well as create and edit new features as needed, that the rest of the crew will view. However, I have been wanting to switch from iPhone to Android for a very long time. Currently the ONLY thing holding me back is that I can't get the same functionality on Android that I can on Iphone. Yes, Google's native maps and Google Earth apps allow you to VIEW the My Maps, but there are several limitations which make this essentially useless. What I need is to be able to:   * View multiple maps overlaid (Google only shows one at a time)   * View ALL features on each map (Google caps off at around 200 and arbitrarily decides which ones to show)   * View and edit these features (Google does not match the names with the features, and does not allow existing features to be edited). Is anyone aware of an app that allows this on Android? Alternatively, is anyone aware of another strategy that might allow KML files to be accessed via mobile phones over 3G?


title: How to make 3D building models from lidar data?

text: I have Lidar data, ground points and first returns. I have processed it and have ready tiles from which I am trying to build 3D building models. So far I was able to just extrude the buildings in the ArcScene and that is it. I am a beginner in GIS/Remote Sensing and any help would be highly appreciated. Also, I am trying to work this out in Arc10.


title: How to register table with geodatabase with a command line?

text: Is there a way to register with geodatabase a table stored in an Oracle 11g database, without using ArcMap or Python script? I'd like to do it with a command line like when registering with ArcSDE (sdelayer -o register ...)


title: How can I create 1 seamless tab file for master map?

text: This is my first post on here so please excuse if it's in the wrong place. I've processed some mastermap files which display absolutely fine. However, during the processing it would not let me create just 1 seamless tab file as the index was over 2Gb. I therefore had to split it in to attributes, i.e. Topo_Point_Seamless, Carto_Area_Seamless etc. Now although these work if you load them all separately, it makes it more awkward to work with. Ideally I'd like just 1 seamless tab file i can turn on and off in MapInfo and geoset's etc. However i've tried the seamless creator in MapInfo (11.5) but it says I cannot create a seamless file from other Seamless files. Does anyone know how i could go about this? Many thanks, Nick


title: arcgis server performance comparison

text: Are there any recent comparisons/studies of arcgis server vs. open source gis servers in common tasks like providing a wms/wfs, filtering data/search performance, and so on? We received a lot of critic about featuring this thesis on our blog as the thesis is outdated and a little bit picky on serving WFS only. Citation from reddit: "This is like comparing two cars based on cup holder design." There is an older presentation given here but like written: It's from 2010


title: How to use wildcard when appending PDFs using ArcPy?

text: I'm working on a script that (as written) produces a large number of individual pdf pages (>100), instead of hardcoding the names of each of the outputs, I was thinking of using a wildcard. However, neither of the following seem to work:                pdfDoc.appendPages("C:\\temp\\*.pdf")      or                pdfDoc.appendPages("C:\\temp\\"*".pdf")      Below is the full script:               mxd = arcpy.mapping.MapDocument("C:\\temp\\Graffiti_DM3.mxd")     df = arcpy.mapping.ListDataFrames(mxd, "*")[0]     graffitilyr=arcpy.mapping.ListLayers(mxd, "*", df)[0]          newquery = '"Incident_Z" =' + "'10003'"     graffitilyr.definitionQuery=newquery     features = arcpy.SearchCursor(graffitilyr)     count = 1               for feature in features:             df.rotation = "0"             arcpy.SelectLayerByAttribute_management(graffitilyr, "NEW_SELECTION", "FID=" + str(count))             #df.extent = graffitilyr.getSelectedExtent(True)             df.zoomToSelectedFeatures()             df.scale = df.scale * .001             arcpy.RefreshActiveView()             df.panToExtent(graffitilyr.getSelectedExtent())             #df.scale = "1000"             outpdf="C:\\temp\\" + str(count) + ".pdf"             arcpy.mapping.ExportToPDF(mxd,outpdf)             count = count+1          pdfpath="C:\\temp\\pdfbook.pdf"     pdfDoc = arcpy.mapping.PDFDocumentCreate(pdfpath)     for pdf in pdfs:              pdfDoc.appendPages("C:\\temp\\*.pdf")     pdfDoc.saveAndClose()     del features     del pdfDoc      Any help would be much appreciated


title: Classic United Kingdom Projection?

text: When you visit a classroom in the United States and see a map of the contiguous states on the wall, it's usually in NAD83 Albers Equal Area with the central meridian right at 96 West. Is there an equivalent for the UK? Or to phrase it another way, what's the most common projection when viewing a map of the UK?


title: Is there a way to identify features on OpenStreetMap.org?

text: Is there an identify features option for OpenStreetMap.org (OSM.org)? Stemming from this question: How to get shapefile of river from OpenStreetMap? I'm looking for the ability to use OSM.org to figure out the tags and values for features I want to download through the Overpass API, or even QGIS for that matter. If there was a way to use OSM.org to browse the data of OSM through the map interface, then click on a feature to get the information about it. Currently, if you know what a feature is called, you can search for it, and you get the information I'm referring to as the output: ![enter image description here](http://i.stack.imgur.com/mBYYe.png) This would be the exact information I would expect if I were to use an Identify button and click on a feature on OSM.org. This would greatly help, as the first referenced post suggests, the ability to pull down specific data from Overpass API or QGIS. Thanks!


Annotation Published to my ArcGIS Server is being displayed as a point when I build my application with ArcGIS Viewer for Silverlight. How do I get the annotation to display as text and not a point feature?


title: Applying distance/area dependent procedures in geographic coordinates over large areas?

text: I'm currently doing some niche modelling with data in geographic coordinates spanning all of Latin America. I am using R for my analysis. Some of the tasks I have to perform are area/distance dependent and yield incorrect results when applied in geographic coordinates. For example, spatial sampling (sp::spsample) leads to a bias towards more points in southern South America and buffering points and polygons (rgeos::gBuffer) leads to stretching/skewing of the buffer away from the equator. What's the best way to address this issue? Should I project my data, apply the necessary procedure, then project back? If so, what projection would work for such a large area? Ideally, I'm looking for a general approach that could be applied to any region, but I fear this doesn't exists.


title: Where can I find block (sub-district) level boundaries for India?

text: I'm looking for block (sub-district) level shapefiles for India. The GADM database only has files for the 2300 Tehsils, whereas I'm looking for the about 5500 Census blocks (sub-districts) for any year possible. Thanks.


title: Info Template not getting displayed

text: I have a layer like,               realTimeLayer = new esri.layers.GraphicsLayer();      I have an infoTemplate like,               infoTemplate=new esri.InfoTemplate("/*Content goes here*/");      I am creating a graphic, adding it to the layer and plotting the point like,                var graphic= new esri.Graphic(candidatePoint,candidateSymbol,attr,infoTemplate);      realTimeLayer.add(graphic);      map.addLayer(realTimeLayer);      points.addPoint(candidatePoint);      In this case the info template works fine. But during some other event, I am adding the graphic like               var graphic= new esri.Graphic(candidatePoint,candidateSymbol,attr,infoTemplate);     map.graphics.add(graphic);     points.addPoint(candidatePoint);      But this time, I don't get the info template in IE9 even after clicking on the push pin. What could be the reason? And how should I correct this?


title: Mapping Kansas City based on X and Y Coordinates in R

text: I have a csv file in R that has coordinates of Kansas City, however, the values at these points aren't in any traditional longitude/latitude formatting that I am aware of. My data looks like this:               Tract            X               Y     218.01       2797574          1143595     218.01       2799346          1144802     133.02       2786207           998301     ...      Is there anyway that I can generate a map from this type of data that will also eventually allow for the visualization of other variables such as population or crime rate?


title: MapServer and OpenLayers: Properly using MAXSCALEDENOM in mapfile (google baselayer)

text: My base layer: google map My points on the map: from an airports.db file, set up as a WFS layer in the OpenLayers code. As of now: map loads, points in the correct places on the map. My goal: when the page loads in the browser, the points should not display on the initial zoom. When you zoom in to the second level, they should. I've successfully done this before using just an OpenLayers map, but never a google map as a base layer. For `MAXSCALEDENOM` I've tried all sorts of values, but the points are always displaying. Here is the Airports `LAYER` in the mapfile.               LAYER                 NAME Airports                 DEBUG 5                 METADATA                         "wfs_srs" "EPSG:3857"                         "wfs_title" "Airports"                         "gml_featureid" "Name"                         "gml_include_items" "all"                         "wfs_enable_request" "*"                         "wfs_encoding" "UTF-8"                         "wfs_extent" "-180 -90 180 90"                 END                 TYPE point                 STATUS ON                 MAXSCALEDENOM 20000                 LABELITEM 'Name'                 CONNECTIONTYPE OGR                 CONNECTION '/tmp/navhome/navdata/airports.db'                 DATA 'select "POINT(" || LONG || " " || LAT || ")" as WKT_GEOMETRY, * FROM airports WHERE name LIKE "S%"'                 #DATA 'select "POINT(" || LONG || " " || LAT || ")" as WKT_GEOMETRY, * FROM airports WHERE ICAO != "ZZZZ"'                 PROJECTION                         "proj=latlong"                         "ellps=WGS84"                         "datum=WGS84"                 END # PROJECTION         END # LAYER      Here is the OpenLayers code for the WFS layer:               var wfs = new OpenLayers.Layer.Vector("Airports", {                 styleMap: new OpenLayers.StyleMap({                         externalGraphic: "red_pin.png",                                 graphicWidth: 18,                                 graphicHeight: 18                 }),                 wrapDateLine: true,                 strategies: [new OpenLayers.Strategy.Fixed()],                 protocol: new OpenLayers.Protocol.WFS({                         srsName: "EPSG:3857",                         version: "1.0.0",                         url: "http://my_url/cgi-bin/mapserv.exe?map=/tmp/navhome/navcore/airports2/airports_viewer.map&SERVICE=WFS&VERSION=1.1.0&srsName=EPSG:3857",                         featureType: "Airports"                 })         });         map.addLayer(wfs);      Any input on this matter is appreciated. I assume it has something to do with the extent, but I'm not totally sure.


title: Geoserver WFS Point Coordinates

text: I'm using SharpMap's desktop client within a C# WinForms application and GeoServer/PostGis as datasource. When i add a WFS layer that contains Multi-Line it works like charm. But if i add a point layer, SharpMap.Utilities.Wfs.ParseCoordinates function parses coordinates with the wrong order. I meant that sharpmap assigns X coordinate to Y and Y to X. ParseCoordinate function is like that:               protected Coordinate[] ParseCoordinates(XmlReader reader)     {         if (!reader.Read()) return null;              string name = reader.LocalName;         string coordinateString = reader.ReadElementString();         var vertices = new List<Coordinate>();         string[] coordinateValues;         int i = 0, length = 0;              if (name.Equals("coordinates"))             coordinateValues = coordinateString.Split(_Cs[0], _Ts[0]);         else             coordinateValues = coordinateString.Split(' ');              length = coordinateValues.Length;              while (i < length - 1)         {             double c1 = Convert.ToDouble(coordinateValues[i++], _FormatInfo);             double c2 = Convert.ToDouble(coordinateValues[i++], _FormatInfo);                  if (name.Equals("coordinates") || name.Equals("pos"))                 vertices.Add(new Coordinate(c1, c2));             else                 vertices.Add(new Coordinate(c2, c1));         }              return vertices.ToArray();     }      I've changed order of c1 and c2 values to fix that problem but not sure why it's acting like that. Here goes the GML output of my geoserver:               <gml:coordinates xmlns:gml="http://www.opengis.net/gml" decimal="." cs="," ts=" ">     51.52672002,25.33049886     </gml:coordinates>      X = 51 and Y = 25 but reader.ReadElementString(); statement at ParseCoordinates function reads them like that: 25 and 51 I'm not sure if it's due to geoserver's configuration or sharpmap's bug.


title: CartoDB preferred file type

text: I was wondering if anyone knew the preferred upload file type is for CartoDB.I upload zipped .shp files for polygons, and sometimes .csv file for point data. I noticed that when I uploaded a zipped .shp file it hit the 50mb account limit even though the file only ended up being 9mb in cartodb when I upgraded accounts. Is this becasue it is duplicating the table during processing? Is there a preferred file type for cartoDB? GeoJSON? thanks,


title: Questions regarding the processing of mosaic Landsat 8

text: **Approximate Translation of Original Question:** A question of theoretical and practical about the mosaic of satellite images. I have 3 different months of Landsat images, which need to accomplish the end MaxVer classification. At what point should I make the mosaic of 3 scenes? before or after sorting? Follow the steps ... Acquisition >> Georeferencing >> Atmospheric Correction >> Classification >> Mosaic or Acquisition >> Georeferencing >> Atmospheric Correction >> Mosaic >> Classification **Original Question:** Uma dúvida de ordem teórica e prática a respeito do mosaico de imagens orbitais. Possuo 3 imagens Landsat de meses distintos, a qual necessito realizar ao final uma classificação MaxVer. Em que momento devo realizar o mosaico das 3 cenas? antes ou depois da classificação? Segue as etapas... Aquisição>>Georeferenciamento>>Correção atmosférica>>Classificação MaxVer>>Mosaico ou Aquisição>>Georeferenciamento>>Correção atmosférica>>Mosaico>>Classificação MaxVer


title: Interpretation of ArcGIS Kernel Density legend parameters

text: In ArcMap 9.3 I've used Kernel Density to map various incidents, but the resulting shapefile doesn't display any units of measurement. Is there a good, not-to-technical source that would explain in lay terminology the interpretation of the output values in terms of the input cell size and search radius?


title: Get number of rules in geoserver style using webservice

text: how can I know the default style for a layer in geoserver, and the number of rules it contains, using a REST call or service? Thank you bye stefano


title: How to pre-populate JavaScript forms in ArcGIS Web Application?

text: I have an ArcGIS web application and would like to pre-populate order numbers. Say if order 1 has been entered for a point on a map, when the next point is added, order 2 would populate automatically. Of course these numbers are a series of letters and numbers, say abcd1234, next order would automatically populate abcd12345. I am using SQL server to submit new features with feature access capabilities. The web application is using the JavaScript API. I have the code hosted on an IIS web server in order to access the actual JavaScript.


title: OpenLayers Clustering Adds Features to Selected Features Array on Zoom

text: This is my first posted question despite using StackExchange for a few years now. I've been stuck for quite some time trying to figure out this issue with OpenLayers. The problem is as follows: I have a bunch of features plotted on a map and am able to select features in various ways (e.g. drawing a box around them, clicking on individual features, entering a zip code that selects all features within that zip code, etc.). This all works fine. The issue is that I have a clustering strategy on the features layer and when I zoom out, currently selected features get clustered together with unselected features and those unselected features will become selected when I zoom back in. What I'm trying to accomplish is to NOT have the previously unselected features become selected after clustering. Here is a link to a visual hosted on my website (http://korynorthrop.com/uploads/Visual.jpg) I could post some code as well, but thought I'd simplify this question first.


title: Continuous color band for a raster legend in QGIS composer legend

text: I'm looking for a way to insert continuous colour legend for a raster layer. I managed to create discrete symbols with values along the bar, but I would like to have a continuous bar with some values on the bar, similar to the old plugin "One-Band Raster", that it doesn't work with the Qgis version higher than 1.7 I tried to set in properties->style, where I selected the Render Type "Singleband pseudocolor" and see the "generate new color ramp" box. I tried both "Continuous" or "Equal interval". However, in the Composer the legend always is displayed with discrete symbols. Is there any way that can make QGIS Composer to display a nice continuous band for a raster layer?


title: Help, I have 10,000+ .shp files covering 3 UTM zones with unknown projections

text: I was handed a ton of property files with no projections attached but was told they are all in UTM zone 21S. Well once I merged them all together and projected them I realized that the data covers 3 UTM zones 20S, 21S, and 22S. So basically each property file could be one of 3 UTM zones. Is there a way to figure out what UTM zone each property file belongs too without checking each one at a time? I can python code if need be. I'm using ArcMap 10.2 I tried that. Here is a more detailed explanation of my situation. These files were original downloaded by UTM zone. So all properties in UTM z 20S were downloaded together with a projection of UTM z 20S, all properties in UTM z 21S were downloaded together with a projection of UTM z 21S ect. The shape files I received didn't come with the .prj files and I have no way to obtain them. So when I project them all with the same coordinate system say UTM z 21S, they overlap, all of the property files from UTM z 20S and 22S are overlapping the correctly referenced properties of z 21S. The misplaced properties are located in the correct area of the UTM zone, just the wrong one. Does that mak sense?


title: How do I get GeoJSON data from GeoServer into show up on my Leaflet map?

text: I'm trying to create a simple "Green" map for my city. The idea being that it is a website with a map and layers you can turn on and off. I am serving the data through GeoServer and am using GeoJOSN to add the data to the map. I checked the URL for the data and GeoServer is serving up the file in GeoJOSN format. Eventually I would like to be able to let the user to click on the points on each layer to get more information but I need to get the data showing up on the map first. The leaflet code below is pretty simple but I can't seem to get the data to show up on the map. I refereed to the Leaflet documentation but can't figure it out. I'm guessing I am just leaving something simple out due to my lack of knowledge about Leaflet. Leaflet web mapping ninja's, a little help please.


title: How <-> operator work?

text: I'm doing a search to find links close to a point.   The problem is if i limit to 20 result i found the closer link at 6mts   but if i limit to 10 result that link doesnt appear and the closest link is at 17mts Shouldn't the top 10 closest link to my point appear in both querys? I'm using PostgreSQL 9.3.2, POSTGIS=2.1.1 ![query result image](http://i62.tinypic.com/2a65agm.jpg)                   SELECT              S.Seg_ID,              S.Link_ID,              S.Course,              ST_Distance(ST_GeomFromText('POINT(-66.879670 10.497950)',2309), S.geom  ) as distance         FROM vzla_seg S         ORDER BY S.geom <-> ST_GeomFromText('POINT(-66.879670 10.497950)',2309)         LIMIT 20


title: Extract raster values under under each polygon buffer in qgis

text: In qGIS: I have raster file and >100 points. I constructed buffers polygons (vector-->geoprocessing-->buffer) around each of the points. Than I clipped (raster-->extraction-->clipper) the raster file with the buffer polygons, resulting in a further raster layer. How do I get the raster values under each of the polygon buffers? My tiff as pixels, which can have five different values (222, 333, ...) and I want to know how often each value appear under each polygon (e.g. polygon I: 10 x values 222, 100 x value 333...; polygon II: 35 x values 222, 50 x value 333...; )


title: Finding the Default version in a large linage tree

text: I have a very large versioning linage tree. I am using the arcGIS 10.0 GDBT tool. How can I find the node with the default version? ![enter image description here](http://i.stack.imgur.com/Ptz2a.png)


title: ArcGIS online live data

text: Does ArcGIS online have the capability of providing real time (live) tracking of assets? If I have a staff member with a cell phone can ArcGIS online show me where that person is at any given moment? I have ArcGIS desktop. I do not have ArcGIS Server.


title: Merge lines with the same 'name' attribute on QGIS?

text: I have a shapefile with a few thousand lines, each is a street segment. I'd like to merge all features with the same 'streetName' attribute into a single attribute. I only really care about the name of the street and it's geometry. What would be the query necessary to do such a thing? I'm new to GIS in general and have run into a problem I can't figure out how to fix. Excuse my noobness. Thanks in advance.


title: Taking attributes from state road data and adding to over lapping route data

text: Ok so... I am working on a project that is looking at the regional public transportation within my area and comparing it to the overlapping road systems. I want to display the different attributes that these routes run on to give recommendations as to what kind of development should go in certain areas. So the attributes that would be in play would be speed, sidewalk access, bike lane access, as well as a couple of other things. I have created the shapefile of the routes and I am going to compare them to the ISRN shapefile that I got from DOT. How would I go about taking these ISRN attributes and moving them over the shapefiles that I have created as well as the bike paths and sidewalks.


title: ArcPy Create Point Along Line A Specified Distance From A Specified Point

text: In ArcGIS Basic I am trying to take a road polyline(one segment for the entire length of road) and get a point placed on the map a specified distance from an intersection. I have a layer of intersection points and the layer of roads. I need to tell it to place a point on the polyline where the starting intersection point is. Then travel the specified distance down the line and place a point for a sign. Any suggestions? I am currently trying to write a script to select a road and points and create a temporary layer of each. Then from that temporary road I CreateRoutes_lr. Then take the intersection point and LocateFeaturesAlongRoutes_lr. From locating the point I would like to take the "fmp" from distance as a variable. Is there any way to extract it?


title: Layer tick boxes embedded in PDF map

text: Maybe this is more of a PDF question than a GIS question. Today, I received a layered PDF map shown as a landscape A3, with a square map on the left, and with the map title, the legend, the scale, a tiny overview map etc. on the right-hand side. However, each legend item had two tick boxes beside it: one for switching on/off the symbols and one for switching on/off the labels. It also contained radio buttons to switch between various background maps. I thought this was really elegant, and would like to create something similar. I know how to create a layered PDF in ArcGIS, and then go to the Layers toolbar in Adobe Reader to switch on/off layers in the map. However, I'd like to know how I can embed the layer control inside the map. (Note: I am not looking on how to draw individual tick and radio boxes, as I feel it's done by some application (semi-)automatically) EDIT: Here is a picture of it, cut down to the bare minimum for copyright reasons ![Tick boxes and radio button in PDF map](http://i.stack.imgur.com/Ie2BS.jpg)


title: Create a raster from file geodatabase featureclass based on a field from a related table

text: I have a USDA Statsgo soil polygon featureclass stored in a file geodatabase along with an attribute table; I have a one-to-one relationship class defined based on the soil type field. When I bring the featureclass into ArcMap I can view the attributes in the related table when I use the identity tool. What I would like to do is create a raster from this featureclass using values from the related table. This seems simple enough, but when I used the Polygon to Raster conversion tool the only fields available to use are those which are stored in the featureclass itself. If I define a join on the same featureclass and table I can then select the appropriate field in the joined table, which makes me wonder why I bothered creating the relationship class in the first place (resulting raster values are numeric - I lose my classification codes of A through D). Can anyone explain 1) why I can't use the relationship class to access fields in the related table, and 2) how to keep the alphabetic soil classification codes?


title: Fastest way to transform WKT in python

text: Is it possible to use pyproj to transform WKT without creating geometry objects such shapely geometry, GEOS geometry or OGR geometry. If not what is the fastest way to do it?


title: esri ArcGIS javascript api compatibility with Internet Explorer 7 and 8

text: I've developed a web map using esri arcgis javascript api version 3.8 and its not compatible with IE 7 or 8. I've advised the people who have reported this to update their version of IE but, if it doesn't involve rewriting the website, I'd like to offer something that can work with these versions. Here's the website: http://il.water.usgs.gov/data/mchenry/map.php As you can see, in IE 7 or 8 the map doesn't even load which throws javascript errors. _EDIT_ I have added some script that changes the markup and css to handle HTML in IE versions lower than 9 found here: http://stackoverflow.com/questions/5449502/html5-and-css3-for-ie7-and-ie8 Since doing this, I've noticed that the map loads on zoom in/out but none of my points load.


title: Fixed Zoom Extend For Overview Map on ArcGIS Server JavaScript API 3.8

text: Can some one please let me know if it is possible to have the Map Overview with Fixed extend and none zoom in or out functions? I know this was not possible on v.3 but how about v3.8? Thanks


title: Edit web.xml in Tomcat 8

text: Need to edit the web.xml tomcat 8 to increase the size allowed for me to do deploy war file geoserver. I was in C: \ Program Files \ Apache Software Foundation \ Tomcat 8.0 \ webapps \ manager \ WEB-INF and open the notepad in the sublime but I can not edit it says is being used by some program. Have stopped tomcat and I closed the program and web administration page. Where this file may be in use?


title: Change from polygon to multipolygon in postgres

text: can somebody describe please how to do the transformation from polygon to multipolygon in the postgres DB, is not possible for me to save my edited data. All data loaded just need to do some editing upon my polygons Thanks in advance regards


title: Using 'construct points' on all lines in a shapefile (arcpy)

text: I need to use the editing tool 'Construct Points' on several shapefiles that each have ~1300 polylines. I need to do this to create points every 1km along each of the polylines, which represent the orbits of a satellite. I've tried several other methods such as densify but none of them solve my problem as well as 'Construct Points'. The only problem is that you can only edit one line at a time, which would obviously take too long to do manually. So I'm wondering, is there a way to implement 'Construct Points' in python so that it loops through all polylines and applies the tool with the same parameters? I feel like it should be possible, but I've never used editing tools in arcpy so I'm not sure. The output of 'construct points' is a point shapefile for the points that were constructed along the line used as an input. I'm wondering if I could do it in python such that it constructs the points on the line, then appends the output points to a larger point shapefile. The following is pseudo code of what I have in mind, although the syntax is of course not correct. If anyone could shed any light, that would be great. I really would like to know if this is even possible with arcpy before I spend time trying to pursue it. Thanks!               # using just one shapefile as an example (pseudo code, incorrect syntax):          orbitshapefile = (shapefile with ~1300 polylines)     pointshapefile = (empty shapefile to store points)          StartEdit(orbitshapefile) # open editing on the shapefile     for orbitline in orbitshapefile:         editor.ConstructPoints(orbitline, 1km, ..(other parameters).., outputpoints) # call editing function         pointshapefile.append(outputpoints) # append output points to 'big' shapefile      EDIT: I just did construct points manually on just 2 lines in editor and found that using the same point shapefile on the second line that was used with the first adds those points to that shapefile, doesn't overwrite them. However, I need an indication of which line these points came from. There is an 'Id' attribute in the resulting table, but everything is 0. Would there be a way to carry over the line's FID to all the points that came from it (alt. I could do a separate script after that uses intersect b/w line and points to come up with a lineID attribute). An updated version of the above code:               for orbitline in orbitshapefile:         lineID = orbitline.FID         editor.ConstructPoints(orbitline, 1km, ..(other parameters).., outputpoints)         outputpoints[pts that were just created in this iteration].SetAttribute(lineIDField = lineID)               Some additional thoughts: I'm not sure how to subset the points that were just created-- maybe I could subset by saying something like: points in outputpoints s.t. lineID == 0 (because all points that were created in previous iterations will have their lineID attribute set to something other than 0) Any thoughts on this rationale? Thanks so much PS: Can anyone AT LEAST tell me if using an editor tool iteratively in python is even possible?


title: GeoDjango Changing My Coordinates

text: I'm running into a problem when saving a coordinate into one of my models in (Geo)Django. My point field is setup as:               location = PointField(blank=True, null = True)      So while running some tests in the django shell I noticed my coordinates were changing...               obj.location = 'POINT(-120.18444970926208 50.65664762026928)'     obj.save()     print(obj.location.x)      **-120.1840742**               print(obj.location.y)      **50.6569673** Another thing I've noticed If I set the point using an UPDATE in postgres and then...               print(obj.location.y)      **50.65664762026928** but if I run a save and check the coordinate again I get...               obj.save()     print(obj.location.y)      **50.6569673** Anybody have any ideas where I can start looking to sort this out or any suggestions what may be causing this. Thanks, Dustin Software:   * Django 1.3.1    * GDAL 1.9.1    * Python 2.7.2    * PostGIS 1.5.3    * GEOS 3.3.3   * PostgreSQL 9.1.4


title: Openlayers WMS layer always sends "EPSG:900913" as srs parameter

text: Why the wms layer ignoring srs parameter srs:"EPSG:4326" and always sending "EPSG:900913" as srs value..? this the code,                   map.projection= new OpenLayers.Projection("EPSG:4326");              var wmslayer = new OpenLayers.Layer.WMS("Cities",      "http://{208.184.113.91}/ArcGIS/services/BurkeMcDowell2013_06_26FORGOPI/MapServer/WMSServer?service=WMS",                                                          {layers: '0',             srs:"EPSG:4326",             transparent: true            },{            projection: new OpenLayers.Projection('4326')           });


title: How to display a table with more than one geometry column in QGIS?

text: In PostGIS, it is possible to use addgeometrycolumn more than once on the same table. Example: I have one table "scenes" with many attributes, and 2 geometry columns: the_geom1 and the_geom2.               select addgeometrycolumn('scenes', 'the_geom1', 0, 'POLYGON', 2);     select addgeometrycolumn('scenes', 'the_geom2', 0, 'POLYGON', 2);      In QGIS, when I display "scenes", it only displays the_geom1. To display the features using the second geometry, I found a solution with :               create view v2 as select id,the_geom2 from scenes;      But I'd like to access all other fields; something like:               create view v2 as select * (except the_geom1),the_geom2 from scenes;      Any ideas ? Thanks, Eric.


title: Creating a Raster Filter

text: I am looking to create a raster filter which reclassifies single or contiguous pixels less than, for example, 30 pixels (shown in white below) into another class (shown in black). The filter, majority filter, aggregate and other spatial analyst generalization tools seem to be too over reaching in this case. Any Python or Arc solutions and ideas would be appreciated. ![enter image description here](http://i.stack.imgur.com/JmyJ4.jpg)


title: Extract portion of WMS raster layer

text: Using QGIS on Windows XP based desktop. I am trying to use USDA orthophoto data from this WMS server: "http://raster.nationalmap.gov/arcgis/services/Orthoimagery/USGS_EDC_Ortho_NAIP/ImageServer/WMSServer" I can get the layer into QGIS with no problem but the only layer on the WMS server is of the entire United States. I want to extract from that the area I am interested in (a city). I have a vector layer of that city to use as a clip coverage but GDAL can only clip from a local raster layer. I have probably left out some pertinent information so please just ask if you need more information.