title: Importing GPS coordinates to extract multiple DN values in ERDAS

text: I have a set of GPS points (stored as a csv for now) and would like to import these points to ERDAS to find the DN values associated with each location. I am able to import a SHP to ERDAS but the only way I have found of extracting the DN values is via the use of the Inquire tool. This enables me to input the lat/lon values manually to find out the DN values but I was wondering whether I am missing a more efficient batch method for doing this using a txt / csv upload?! **Update: managed to solve this in QGIS usinghttp://www.digital- geography.com/qgis-plugins-point-sampling-tool/**


title: How to create static image file with geometries from shapefile using ArcGIS Online?

text: I'm trying to create static map images, including overlaid geometry features imported from a shapefile, server-side (in PHP if that matters), in order to use them in generated PDFs. We're using ArcGIS Online and I had been assuming this would be something easy to do with the REST API but I can't find any examples or even discussion of how to do this. It's trivially easy, by contrast, to overlay geometries from shapefiles on maps using their Javascript API but I can't figure out how to do it server- side. Thanks for any advice.


title: How to compute drive access to a station based on vehicle origin points?

text: I have a shapefile that contains station locations and another shapefile that is the geocoded origins of those who drove to the station. I would like to make a buffer around the station that shows the area where 90% of the vehicles came from so I can show the general drive access area for each station without the extreme outliers. Again, from my current files I would like to create a raster or shapefile that shows the general catchment area of each station.


title: What is correct Where Clause when using variable and Select By Attribute tool within loop?

text: I have seen several similar posts related to this technique but I am still having issues. The following script is intended to select all points from one shapefile that are within the boundaries of a second polygon shapefile, then edit the number of points selected into one of the polygon fields. I have not been able to get the where clause statement to function properly although from the error message it appears to be correct. Both shapefiles are within a file geodatabase and feature layers where created for both in order to use the selection functions. Here is the code:               arcpy.MakeFeatureLayer_management("treatment_maps",'maps')     arcpy.MakeFeatureLayer_management("points_2014",'points')     cursor=arcpy.UpdateCursor("treatment_maps")          for row in cursor:         map=row.getValue('map_name')         where='"map_name"' + '=' + '\'' + str(map) + '\''         arcpy.SelectLayerByAttribute_management("maps","NEW_SELECTION",where)         arcpy.SelectLayerByLocation_management("points","COMPLETELY_WITHIN","maps")         count=arcpy.GetCount_management("points")         row.setValue('DI_Count', count)         cursor.updateRow(row)      Here is the error statement:               Executing: SelectLayerByAttribute maps NEW_SELECTION "map_name"='Airport'     Start Time: Wed Jan 15 13:00:56 2014     ERROR 000358: Invalid expression     Cannot acquire a lock.     Cannot acquire a lock.     Failed to execute (SelectLayerByAttribute).     Failed at Wed Jan 15 13:00:56 2014 (Elapsed Time: 0.00 seconds)      The treatment_maps features class contains polygons and counts (DI_Count) for points located within each polygon, which need to be updated. I'm hoping to get this script to work properly as I have around 100 polygons to update. Thank you for any help or advice!


title: How can I set Google V3 layer opacity in OpenLayers 2.13.1?

text: I've noticed that layer.setOpacity doesn't affect the google v3 layer anymore. Is there an alternative way to set the v3 opacity as base layer?


title: ArcGIS API JavaScript, Flex, and Silverlight updated comparision

text: Does anyone know of a site, blog, etc. that constantly updates a comparison between the various ArcGIS APIs for each new release of the API?


title: Calculate the sum of the geometric length of an arcgis shape file

text: I have a shape file having lots of line routes. I projected it and calculated the length using calculate geometry. now i want to do the sum of all these routes in arcgis. What are the possible ways.


title: capabilities URL for open street map

text: I want to use open street map as cascading WMS to my Geoserver. Which capabilities URL should i use? Can any body tell me the capabilities URL for OSM ? Regards


title: Open a GeoJSON file with QGis

text: If I have a GEOJson file, for example a file named test.js with inside the following:               { "type": "FeatureCollection",     "features": [     { "type": "Feature",       "geometry": {"type": "Point", "coordinates": [102.0, 0.5]},       "properties": {"prop0": "value0"}       },     { "type": "Feature",       "geometry": {         "type": "LineString",         "coordinates": [           [102.0, 0.0], [103.0, 1.0], [104.0, 0.0], [105.0, 1.0]           ]         },       "properties": {         "prop0": "value0",         "prop1": 0.0         }       },     { "type": "Feature",        "geometry": {          "type": "Polygon",          "coordinates": [            [ [100.0, 0.0], [101.0, 0.0], [101.0, 1.0],              [100.0, 1.0], [100.0, 0.0] ]            ]        },        "properties": {          "prop0": "value0",          "prop1": {"this": "that"}          }        }      ]     }      how to open it with QGis? I open QGis, press "add vector layer", choose "File" but I am not able to open it. If I choose GEOJson format, my file is not showed in browser window. If I do not specify input type, I get the error: "your file is not valid or recognized data source". How have I to format this file? It seems there a little lack of documentation about this topic...


title: Optimizing OSM tiles download

text: Disclaimer: I don't know anything about maps, please excuse me. In our mobile app, we use Cloudmade as our map tiles provider. For certain places, we download multiple tiles for each zoom level (level 13 through 19) to be used offline. The thing is that sometimes you need to download tiles for thousands of places at once, making it a very lengthy download. We are trying to optimize this process and there are a few things that I've come up with, which I don't know if they even make sense:   1. Cloudmade serves the tiles as png files. We could try to cache those tiles in our server, and reduce their size.   2. Skip some zoom levels, let's say we only download one every two (13,15,17,19).   3. Download a higher quality version on a higher zoom level, let's say 10 and split it up in the client. Do any of these sound reasonable? What else could we do? Thanks a lot!


title: Why a shapefile of a geological map is not the fault lines are displayed in QGIS?

text: I am a new user of QGIS 2.4. I do not display correctrly geological SHP layer. This shape contains polygons with shapes, fault lines, dips, volcanic structures, shoreline types, geological boundaries between other layers. Characters set: utf8. Reference systems: UTM WGS 84, Zone 17 North. This layer was digitized from the geological map of Panama, 1:250,000. When using Arcgis are fault lines in QGIS but can not be observed. This shape contains the following extensions: .shp, .shx, .dbf, .prj, .qix., shp.xml. Help me


title: alignment of shapefiles in GRASS

text: I have two shapefiles below. One is in the state plane projection, and the other is not projected. In ArcMap the shapefiles line up perfectly. When I import them to my mapset (using the override option), they do not line up. Anyone know why this is happening? I've looked through the posts and can't find an answer. Thanks Data Type: File Geodatabase Feature Class Feature Type: Simple Geometry Type: Line Geographic Coordinate System: GCS_North_American_1983 Datum: D_North_American_1983 Prime Meridian: Greenwich Angular Unit: Degree Projected Coordinate System: NAD_1983_StatePlane_Texas_South_Central_FIPS_4204_Feet Projection: Lambert_Conformal_Conic False_Easting: 1968500.00000000 False_Northing: 13123333.33333333 Central_Meridian: -99.00000000 Standard_Parallel_1: 28.38333333 Standard_Parallel_2: 30.28333333 Latitude_Of_Origin: 27.83333333 Linear Unit: Foot_US Geographic Coordinate System: GCS_North_American_1983 Datum: D_North_American_1983 Prime Meridian: Greenwich Angular Unit: Degree


title: geoserver fails to start with jdbcconfig community module

text: I have installed the geoserver JDBC config community module, but now geoserver fails to start. I have set the properties to connect to an oracle database like this:               #Default GeoServer JDBC config driver and connection pool options. Edit as appropriate.     #Mon Aug 18 13:46:46 CEST 2014     enabled=true     initdb=true     initScript=${GEOSERVER_DATA_DIR}/jdbcconfig/scripts/initdb.oracle.sql     import=true     jdbcUrl=jdbc:oracle:thin:@xxx.xxx.local:1521/dbname     driverClassName=oracle.jdbc.driver.OracleDriver     username=username     password=password     pool.minIdle=4     pool.maxActive=20     pool.poolPreparedStatements=true     pool.maxOpenPreparedStatements=100     pool.testOnBorrow=true     pool.validationQuery=SELECT SYSDATE FROM DUAL;      However when I start geoerserver it returns the following stacktrace: > 2014-08-28 17:40:19,377 ERROR > [org.springframework.web.context.ContextLoader] - Context initialization > failed org.springframework.beans.factory.BeanCreationException: Error > creating bean with name 'geoServerLoader' defined in URL > [jar:file:/C:/Program%20Files%20(x86)/GeoServer/webapps/geoserver/WEB- > INF/lib/gs-main-2.5.1.jar!/applicationContext.xml]: Initialization of bean > failed; nested exception is > org.springframework.beans.factory.BeanCreationException: Error creating bean > with name 'JDBCGeoServerLoader' defined in URL > [jar:file:/C:/Program%20Files%20(x86)/GeoServer/webapps/geoserver/WEB- > INF/lib/gs-jdbcconfig-2.5-SNAPSHOT.jar!/applicationContext.xml]: Cannot > resolve reference to bean 'JDBCCatalogFacade' while setting bean property > 'catalogFacade'; nested exception is > org.springframework.beans.factory.BeanCreationException: Error creating bean > with name 'JDBCCatalogFacade' defined in URL > [jar:file:/C:/Program%20Files%20(x86)/GeoServer/webapps/geoserver/WEB- > INF/lib/gs-jdbcconfig-2.5-SNAPSHOT.jar!/applicationContext.xml]: Cannot > resolve reference to bean 'JDBCConfigDB' while setting constructor argument; > nested exception is org.springframework.beans.factory.BeanCreationException: > Error creating bean with name 'JDBCConfigDB' defined in URL > [jar:file:/C:/Program%20Files%20(x86)/GeoServer/webapps/geoserver/WEB- > INF/lib/gs-jdbcconfig-2.5-SNAPSHOT.jar!/applicationContext.xml]: Cannot > resolve reference to bean 'jdbcConfigDataSource' while setting constructor > argument; nested exception is > org.springframework.beans.factory.BeanCreationException: Error creating bean > with name 'jdbcConfigDataSource': FactoryBean threw exception on object > creation; nested exception is org.apache.commons.dbcp.SQLNestedException: > Cannot create PoolableConnectionFactory (ORA-00911: invalid character )   >  etc lots more stack   >  etc   >  etc   >  etc   >  .   >  .   >  .   >  Caused by: java.sql.SQLException: ORA-00911: invalid character at > oracle.jdbc.driver.DatabaseError.throwSqlException(DatabaseError.java:112) > at oracle.jdbc.driver.T4CTTIoer.processError(T4CTTIoer.java:331) at > oracle.jdbc.driver.T4CTTIoer.processError(T4CTTIoer.java:288) at > oracle.jdbc.driver.T4C8Oall.receive(T4C8Oall.java:743) at > oracle.jdbc.driver.T4CStatement.doOall8(T4CStatement.java:207) at > oracle.jdbc.driver.T4CStatement.executeForDescribe(T4CStatement.java:790) at > oracle.jdbc.driver.OracleStatement.executeMaybeDescribe(OracleStatement.java:1037) > at > oracle.jdbc.driver.T4CStatement.executeMaybeDescribe(T4CStatement.java:830) > at > oracle.jdbc.driver.OracleStatement.doExecuteWithTimeout(OracleStatement.java:1132) > at > oracle.jdbc.driver.OracleStatement.executeQuery(OracleStatement.java:1272) > at > org.apache.commons.dbcp.DelegatingStatement.executeQuery(DelegatingStatement.java:208) > at > org.apache.commons.dbcp.DelegatingStatement.executeQuery(DelegatingStatement.java:208) > at > org.apache.commons.dbcp.PoolableConnectionFactory.validateConnection(PoolableConnectionFactory.java:658) > at > org.apache.commons.dbcp.BasicDataSource.validateConnectionFactory(BasicDataSource.java:1558) > at > org.apache.commons.dbcp.BasicDataSource.createPoolableConnectionFactory(BasicDataSource.java:1545) > ... 94 more Bottomline is that it returns an invalid character error.   I have have looked at the provided initdb.oracle.sql, but cannot find anything wrong. Even when I have only 1 line which I typed myself it fails on the ORA-0911 error. For instance the first line:               CREATE SEQUENCE seq_OBJECT;      This also fails on the invalid character error. I also tried this without the semicolon.   Any suggestions?   I am using geoserver 2.5.1 and have downloaded the appropriated JDBCconfig files. I am trying to connect to an oracle 11 database. I am working on a windows 7 desktop. I have only added what I think is the interesting part of the stacktrace, but i can add the rest if needed.


title: Creating centreline from polygon

text: I´´m looking for a tool in Q-GIS that can create a centreline from polygon. Similary to this post: Finding centrelines from polygons? I´m not sure that the Geoprocessing Tool can deal with it. Thx


title: How to support multiple Coordinate systems with ArcGIS Server 10.0 WMS?

text: I have a serve a number of wms Layers with ArcGIS Server (AGS) with the following coordinate system support - EPSG:25832, EPSG:31467, EPSG:4326 and CRS:84. As far as I can tell, AGS only supports 2 coordinate systems in the GetCapabilities document. Is this correct? Do I have another way to server wms layers in multiple Coordinate systems.


title: Limit Zoom on WMTS Layer using Openlayers 2.12

text: I'm implementing a baselayer to be used by a Plone addon (collective.geo.openlayers). My implementation specifically creates a layer and returns it to the addon, which enables many different layers to be configured through a webinterface. I got it working well, but there is one problem I'm facing that I cannot wrap my head around. The map displayed is a small local map which only shows tiles when zoomed in to zoom level 11 and higher. For this reason I would like to limit the user from zooming out farther than that. But when I try to add maxZoomLevel, minZoomLevel or numZoomLevels to the layer definition I get type errors in the getTileInfo function because Openlayers is trying to access the longitude of the tileOrigin property which for some reason is null. I tried different suggestions about limiting the zoom level but to no avail. **Does anyone know how to limit the zoom levels in the following layer definition?** I cannot simply work around the issue by doing something to the map instance itself, as I only control the layer creation.               function() {          var zugmap_matrixids = [{              identifier : 0,             topLeftCorner : new OpenLayers.LonLat(919690.324327, 5987771.047746)         }, {             identifier : 1,             topLeftCorner : new OpenLayers.LonLat(929474.2639475, 5987771.047746)         }, {             identifier : 2,             topLeftCorner : new OpenLayers.LonLat(929474.2639456, 5987771.04773376)         }, {             identifier : 3,             topLeftCorner : new OpenLayers.LonLat(931920.24885072, 5985325.06282864)         }, {             identifier : 4,             topLeftCorner : new OpenLayers.LonLat(931920.24885072, 5984102.07037608)         }, {             identifier : 5,             topLeftCorner : new OpenLayers.LonLat(932531.745078525, 5983490.57415958)         }, {             identifier : 6,             topLeftCorner : new OpenLayers.LonLat(932837.49319014, 5983184.82603666)         }, {             identifier : 7,             topLeftCorner : new OpenLayers.LonLat(932990.36724671, 5983031.95198009)         }     ];           var layer = new OpenLayers.Layer.WMTS({         name : "LuftbildPlus+",         url : [             "http://webdienste.zugmap.ch/OrthofotoPlus_tiled/service.svc/get"         ],         layer : "LuftbildPlusGoogle",         matrixSet : "LuftbildPlusGoogle",         matrixIds : zugmap_matrixids,          tileFullExtent : new OpenLayers.Bounds(919690.324327, 5948635.289264, 978393.96205, 5987771.047746),         maxExtent : new OpenLayers.Bounds(919690.324327, 5948635.289264, 978393, 5987771),         style : "default",         format : "image/png",         isBaseLayer : true,         zoomOffset : -11,         //numZoomLevels: 8     });          return layer     }      You'll notice that numZoomLevels is commented out. I feel like this should be the correct property to set as the matrices only contain 8 levels, but if I uncomment it, I get the type error I was talking about. If I don't use numZoomLevels or any other property that controls the zoom the map works fine.


title: Listing datasets doesn't bring back the correct result

text: If I do a ListDatasets() in python, then I don't get the correct result in the list except for datasets belong to other schema owners, so why can I still see the datasets in ArcCatalog. The code is staightforward:               workspace = r'Database Connections\\enidb.sde'     arcpy.env.workspace = workspace     dslist = arcpy.ListDatasets()      The dslist is empty bar the CULTURAL.ESRI_DATASETS - I should be able to see all the ENIDB ones as I am logged in as that schema owner: ![enter image description here](http://i.stack.imgur.com/Ws9WG.png) ![enter image description here](http://i.stack.imgur.com/L1mg0.png) Any clues folks?


title: Send XML WFS Request via PHP

text:
I'm trying to build a web mapping app that retrieves data from our WFS server. I'd like to avoid JavaScript for the main search scripts. From what's I've read the best way to send the XML POST request to the server is to use CURL. I have tested the following very simple XML in GeoServer and it's fine. Hwever, when I send the request as below I get an error. I have trawled Google but can't seem to find any examples of doing this. I can't believe I'm the only one trying to do this! The PHP Code:               <?php      //Set variables.     $url = "http://squirrel.bis.local:8080/geoserver/wfs";     $post = '<GetCapabilities          service="WFS"          xmlns="http://www.opengis.net/wfs"          xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"          xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0    /wfs.xsd"/>';          $post = str_replace("\r", "", $post); //Clean string.     $post = str_replace("\n", "", $post);     $post = str_replace("\t", "", $post);     $post = str_replace("> <", "><", $post); //Clean any leading spaces between elements.          //What are we sending?     echo "Post Request:
";     echo htmlentities($post);     echo "

";          //Send post request.     $ch = curl_init();     curl_setopt($ch, CURLOPT_HEADER, 0);     curl_setopt($ch, CURLOPT_RETURNTRANSFER,1);     curl_setopt($ch, CURLOPT_URL, $url);     curl_setopt($ch, CURLOPT_POST, 1);     curl_setopt($ch, CURLOPT_POSTFIELDS, $post);     $content=curl_exec($ch);          //Collect result.     $result = curl_exec($ch);     curl_close ($ch);          //What did we get?     echo "Post Response:
";     echo "
" . htmlentities($result) . "
";     echo "

";          ?>      The error is:               <ows:ExceptionReport version="1.0.0"     xsi:schemaLocation="http://www.opengis.net/ows      http://squirrel.bis.local:8080/geoserver/schemas/ows/1.0.0                   /owsExceptionReport.xsd"     xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"     xmlns:ows="http://www.opengis.net/ows">     <ows:Exception exceptionCode="MissingParameterValue" locator="request">     <ows:ExceptionText>Could not determine geoserver request from http request         org.geoserver.platform.AdvancedDispatchFilter$AdvancedDispatchHttpRequest@19f0b0a6     </ows:ExceptionText>     </ows:Exception>     </ows:ExceptionReport>      Thanks Steve


title: Where are the Spatialite settings stored?

text: after moving a spatialite-database to another directory, qspatialite says "Can't connect to Database" / Error unable to open database file. Even no project is open or I try to create a new project. In Qgis browser it lists the spatialite-database but says "database does not exist". I scanned all Qgis project files for the according database, but it does not show up. The questions are:   * where are these settings stored?   * how to change these settings to the correct path?


title: Why is GeoNetwork not displaying any maps?

text: I've just installed GeoNetwork 2.6.4 on CentOS 6.2 using MySQL as the database backend. Now, when I go to the web interface, I see the main web site, but no maps are displayed. On the geonetwork machine (actually it's a VM) I see them, but over the network I don't. Can anyone think of anything I should check or change? I set the jetty.xml file to listen on all interfaces and turned iptables off (it was on by default) so that I could access the web site at all.


title: Resample binary raster to give *proportion* within new cell window

text: I would like to go from a binary raster of forest/non-forest at 30m resolution to a raster at 240m resolution with the value of each cell being the proportion of that cell forested - i.e. the data values for the 240m cells would go from 0 (where all cells in the 30m raster were zero/non-forest) through 0.5 (where half the 30m cells were forest, half were non-forest) to 1 (where all cells in the 30m raster were forested). Although it is has been suggested not to use bilinear interpolation on discrete data, as far as I can tell the result with binary 0/1 data would be to give an average value (i.e. a proportion between 1 and 0). Is this a sensible way to do it, or is there a better method? I can use Arc, QGIS and Idrisi. Thanks for your help.


title: Flattening JP2 files to single band (i.e. averaging bands)?

text: I have created thousands of png tiles using ArcGIS and have merged groups of them into large jp2 files for serving via GeoServer. I need to flatten these to a single band so that I can apply a colour map to achieve transparency. Merge single band and multi band color images is similar but the suggested tool "Composite Bands tool in arctoolbox" only merges single band to 3 band. I need the reverse! This is based on the following question Getting a constant background color for transparency I have arcgis 10.1, qgis/gdal, erdas imagine and global mapper.


title: How to produce a mean raster file based on other raster files?

text: > **Possible Duplicate:**   >  Tools to flatten the jp2 to a single band (average the bands)? I have a number of raster files from a same area, that coordinates of corresponding pixels perfectly match one another. I have to perform simple operations like producing a raster that its pixel values are the average of a 3 raster files corresponding pixel values. Simply speaking, I want a new raster file that is average of 3 other raster files.


title: gdal_translate syntax windows command line

text: I am trying to extract a band from a MODIS hdf file on windows with gdal_translate. Here's my prompt:               gdal_translate -of GTiff "HDF4_EOS:EOS_GRID:\"MOD13A2.A2009001.h17v08.005.2009020092120.hdf\":MODIS_Grid_16DAY_1km_VI:1km 16 days NDVI" NDVI.tif      However I get:               ERROR 4: `'HDF4_EOS:EOS_GRID:MOD11A2.A2009001.h17v07.005.2009014044746.hdf:MODIS     _Grid_8Day_1km_LST:LST_Day_1km'' does not exist in the file system,     and is not recognised as a supported dataset name.          GDALOpen failed - 4     `'HDF4_EOS:EOS_GRID:MOD11A2.A2009001.h17v07.005.2009014044746.hdf:MODIS_Grid_8Da     y_1km_LST:LST_Day_1km'' does not exist in the file system,     and is not recognised as a supported dataset name.      I imagine it's something to do with my syntax being a bit 'iffy' though it seems perfectly fine to me Help much appreciated


title: How to use projection of a shapefile in order to have a better view in OpenGL

text: As you see here, this is the code that I have used to extract vertices of ESRI shapefiles in order to show them with OpenGL.   I have also uploaded my project here, if you want to give it a try and see my problem yourself. ( Compare what is shown in OpenGL with Arcmap, Global Mapper or QGIS ) IT seems that the layer is stretched abit in OpenGL. So two questions rose up for me:   * Can using projection of the layer help me? How? I mean how can I apply the projection of layer to OpenGL?    * As you see the OpenGL is using an ortographic projection (glOrtho). Do you think using (glFustrum) instead will solve my problem and I can ignore projection of layer?


title: Differences between points and lines as poles in an equal area projection

text: I am planning on using an equal area projection in a current project. I used this website to determine which projection I should be using. I hadn't previously considered it, and was wondering if there are any hard set guidelines/rules about the advantages or disadvantages of poles being represented as points or lines (I will eventually be calculating total values of rasters throughout the globe).


title: Is there a world projection that looks like a 3D globe?

text: I would like to try and show data such that it appears as if it was on google earth. For example, in this case, the view of the earth is circular-like and the data shown curves around. Are there any map projections that enable me to do this? (I do understand that geographic projections convert information from a 3D space into a 2D space. What I want to do, just for illustration purposes is to make the 2D space look like 3D!)


title: Where do I find the code for gd_calc.py?

text: I need to create a raster that is the average of a large set of other rasters (many thousands). I've learned from another question that the script gd_calc.py can do a limited version of this, limited to 26 input files. I'd like to modify the script to take an arbitrarily large number of input files, but I don't know where to find it on my system, being new to GRASS, Python, and just about everything else. So, can anyone tell me where to find it?


title: How to export KML from points plotted in javascript Google Maps API

text: I want to export a KML file from points plotted in Javascript using the Google Maps API. Is there a way to do this? I would like the user to be able to log in and see Google Maps, then plot points on the map and finally export those points in a KML format. Thanks.


title: Where are decoration (scale bar, etc) plug-ins in QGIS 1.8?

text: For QGIS 1.8 what happened to the "Decorations" (North Arrow, Scale Bar, etc) plug-in (for the canvas, not composer)? Thanks.


title: Qgis or ArcMap ,find later added data from two layers?

text: i have two layer which consist of 2010 and 2011 years poi. and in qgis i want to find which poi have been added new in 2011 data .. or can we get result from arcgis or any tool ... thnks


title: how to introduce LDAP after installing esri geoportal(windows)

text: I have already installed **ESRI Geoportal** version 1.2.1 and I am upgrading it to 1.2.2. The 1.2.1 version is also customized according to the needs. This version was installed using single administrator but I would like to convert it to be able to use active directory(LDAP). How can this be done without having me to reinstall it from scratch?


title: How to connect QGIS to SDE?

text: I read your amazing solution in this topic Is it possible to connect to ArcSDE using QuantumGIS?. Unfortunately I'm not so expert in open-source products (we use only Esri suite) and I need to ask you some things. - How can I create a .ovr file with my SDE connection? - How can I add it to QGis? Thanks a lot for your help and best regards. Matteo


title: How to build QGIS with Python on Redhat?

text: I'm attempting to run QGIS on Redhat 5 64 bit architecture. Version of Python is 2.4.3. Version of QGIS is 1.7.4 (wroclaw). When I run qgis, the splash screen appears, followed by a QMessageBox that says               Couldn't load PyQGIS     Python support will be disabled          can't import cStringIO      However, here's the kicker. The Python version is blank and the Python path is blank. (QGIS Version is set appropriated). Because of my company's configuration, everything is installed from source. python is on my PATH.               PYTHONPATH  = /usr/lib64/python_2.4:/usr/lib64/python_2.4/site-packages     LD_LIBRARY_PATH = /usr/lib64/python_2.4/site-packages/PyQt4:/apps/qgis_1.7.4/lib/qgis:/apps/qgis_1.7.4/lib      I find it really strange that QGIS doesn't seem to be finding python or the PYTHONPATH environment variable. Any suggestions?? Why is no traceback printed like I see in many other forum questions?


title: Can you permanently set the polygon/polyline segment length?

text: I'm trying to create a several hundred polygons in ArcGIS 10, each of which may have up to several hundred vertices. For consistency sake I would like each line segment in these polygons to be of a set length, say 10m. I know you can right click and set the next segment length (or use the keyboard shortcut "Crtl L") as you are drawing, but this is going to insanely repetitive if I have to do this thousands of times. Is there a way to permanently set the segment length as you are drawing a polygon or polyline? Thanks!


title: How to recreate past landscapes on ArcGIS 10.1?

text: I am trying to recreate the landscape as it may have been around 4000 BC. I have an idea about what tree species existed in particular regions of Scotland (by reading pollen studies). Now I need to somehow add trees to the Land-Form PANORAMA raster (taken from Ordnance survey UK) to visualize the landscape. I am aware of adding layers (with different datasets - for instance geological information) but I am not sure how I can add tree cover in ArcGIS as I don't have this information in the form of a layer which can be added to the raster.I was wondering whether anyone could point me in the right direction. Any help or advice would be greatly appreciated.


I'm using Postgresql 8.4, Postgis 1.5 for My GIS works. Planning to change my Database to Microsoft SQL Server 2012. Can we use all the Postgis functionalities in Microsoft SQL Server 2012. Please guide me...


title: How to position callout bubble leader using maplex?

text: I'm using data driven pages to produce a number of maps. I'm using Maplex to have a callout bubble labeling my features. The problem is that for about 5% of the features, the callout bubble is not touching the feature. I want the tip of the bubble to either touch the boundary or be within the polygon. In the image below, the labels on the right and left (59 and 63) are correct, _while the one in the middle is incorrect because it isn't touching the feature_. Note, these labels are coming from the blue polygons, not the yellow line in the middle. ![enter image description here](http://i.stack.imgur.com/Y86Za.jpg) The relevant settings are shown in the images below: ![enter image description here](http://i.stack.imgur.com/4iDbd.png) ![enter image description here](http://i.stack.imgur.com/aWGNq.png) ![enter image description here](http://i.stack.imgur.com/wXIXn.png) ![enter image description here](http://i.stack.imgur.com/thE00.png) I would think that having the anchor point set as the eroded centre would do the trick. However, it looks like it may be placing the label in the centroid instead. Any thoughts?


title: sum count column in raster attribute

text: I am a beginner to python coding in GIS. I want to calculate the total number of counts in raster attribute table. In order to do that:   1. I would like to build a attribute table of raster data   2. use a loop through attribute table using ARCPY.SearchCursor() Therefore, I wrote a code to build an attribute table. However, I am getting an error while running the code and not able to figure out the error in the code. I would appreciate any help/suggestion! Code:               import arcpy     # Create the geoprocessing object      import arcgisscripting     gp = arcgisscripting.create()          #To overwrite output     arcpy.env.overwriteOutput = True          #Set environment settings     gp.workspace = "C:/Subhasis/Test/raster-read"          #Build Raster Attribute Table     gp.BuildRasterAttributeTable_management("sti-9", "Overwrite")          Error Message:     Traceback (most recent call last):       File "C:/Subhasis/Test/raster-read/test.py", line 15, in <module>         gp.BuildRasterAttributeTable_management("sti-9", "Overwrite")     ExecuteError: ERROR 000049: Failed to build attribute table     Failed to execute (BuildRasterAttributeTable).


title: How to create attribute table for raster data?

text: I have raster dataset for rainfall but its not having attribute table. How can I create tabular data using the same raster file? How to join it back again to it? Software used is ArcMap 10.1.


title: How do I overcome this problem of converting a raster layer to polygon form?

text: I am trying to convert a raster layer of U.S. climatic regions into polygon format using ArcMap 9.3 under Windows XP. The original data are to be found at the following address: Koppen Climate Classification from Idaho's Geospatial Data Clearinghouse Initial attempts to use the Raster to Polygon tool generated "999998 Unexpected Error" and even "Serious Application Error" on a couple of occasions. To investigate further, I then replicated the procedure on another machine (running version 10 under a more recent operating system), which threw up the more meaningful error: "Raster must be integer type". Having tried the "Math - Int" tool, it elicited a "000367 Invalid GP data type” error, so I resorted to the "Math - Times" tool, successively trying multiples from 100 to 1000000, every time receiving the error: "000864 The input is not within the defined domain”. Can anyone suggest a workable procedure that I might use to convert the raster layer of the source data into a polygon format?


title: Missing SRS and NoData when writing TIFF rasters in geotools

text: As I am writing .tif rasters in geotools, the CRS of the written file is not the same of the source (the source raster which I read and processed). My question is : How can I maintain and write .tif rasters under the same CRS of the source raster and how can I register a value as "NoData" for that produced .tif raster? When right-clicking a geotools-produced raster in ArcGatalog and examining its properties, the SRS shows some SRS which is different from the source and "NoData" filed shows -3.40282346639e+038. Can any one advise please?


title: OpenGeo Suite v3.1 Port settings no start.ini file

text: I am trying to change the port settings for the OpenGeo Suite. There is no jetty\start.ini file. I have looked at the two config files that come with the software and tried changing the jetty.port = 8080 in the jetty.xml file. Does anyone know how I can change the port for opengeo suite v3.1?


title: How do I set a single output location for multiple files in ModelBuilder?

text: I have created a model in ModelBuilder. The first step is for the user to specify the location of the File Geodatabase that all the outputs will be saved in. What I want is for all the outputs (from the different tools within the model) to be saved in that FGDB; however, it could be called anything and located anywhere. So how do I cause the output to be saved in the location that is specified in the first step? Thanks!


title: GP Progress Dialog in ArcGIS10

text: I have an ArcGIS custom IGPFunction tool. I want to use the GP progress dialog not the progress bar at the bottom of the ArcGIS window as I don't want people to interact with the map whilst it's processing. In the Execute method you have an instance of ITrackCancel but I can't work out how to show the default GP progress dialog from there. I have this as code at the moment but it still doesn't show the dialog. Any ideas?                       IStepProgressor pStepPro = (IStepProgressor)TrackCancel;             pStepPro.MinRange = 0;             pStepPro.StepValue = 1;             pStepPro.Message = "Starting Import...";             pStepPro.Position = 0;             pStepPro.Show();                  ...                  pStepPro.Hide();      Thanks Rob


title: How to determine line numbers within an LAS file?

text: I am working with mulipulse lidar data that collects points along a number of lines within the flight path. From what it appears there are four scan lines within the data, at different intervals but with identical gps timestamps. I am trying to determine the name and number of individual lines within the las file. I am using liblas module in python. I found this documentation that explains the different fields stored in an las file. It mentions a data field (get_data and set_data) at the very bottom of the page. The 'point data format' and 'point data record length' in the header set aside space for this 'data' field. My header says I have 28 bytes set aside for the data field, and there are 28 values stored in the data field. The 19th value (at least in two datasets from two different sensors) refers to the line number. I have a single value in single pulse data and 4 in multi-pulse data. I was wondering if there is a standard for what is stored in this field or if it is proprietary. Also, as a way to get the name of each scan line, I wrote the following code:               import liblas     from liblas import file as lasfile          # Get parameters     las_file = r"E:\Testing\00101.las"          f = lasfile.File(las_file, mode='r')          line_list = []     counter = 0     for p in f:         line_num = p.data[18]         if line_num not in line_list:             line_list.append(line_num)         counter += 1     print line_list      It results with the following error:               Traceback (most recent call last):       File "D:\Tools\Python_Scripts\point_info.py", line 46, in <module>         line_num = p.data[18]       File "C:\Python27\ArcGIS10.1\lib\site-packages\liblas\point.py", line 560, in get_data         length = self.header.data_record_length       File "C:\Python27\ArcGIS10.1\lib\site-packages\liblas\point.py", line 546, in get_header         return header.Header(handle=core.las.LASPoint_GetHeader(self.handle))     WindowsError: [Error -529697949] Windows Error 0xE06D7363      Does anyone know more about the line numbers stored in the las point/header? Can anyone explain the error?


title: How to digitize polygons from georeferenced raster file

text: I used QGIS to georeference a raster file of a custom US map. The map has features which are not strictly defined by state borders. I would like to digitize the map by creating contiguous borders for each feature. I tried that a few different way, but none of them really worked:   * traced the borders with polylines, and digitizing with the Polygonizer plugin. The plugin kept crashing.   * used Conversion > Raster > Vectorize (Raster to Vector). That one worked, but the result was many overlapping polygons for each feature, each with slightly different borders. I have somewhere between 50-100 different raster images to process, which is why any automated solution is highly preferred. I'm new to GIS software, but could easily pre-process the images if needed, either manually or with a custom script (using Python/PIL, ImageMagick, Photoshop batch processing or similar tools). Thanks for any pointers! **Original** ![Original image](http://i.stack.imgur.com/YELFr.png) **Georeferenced image** ![Screenshot of georeferenced image](http://i.stack.imgur.com/N9Vdw.png)


title: How to find the features, which have an override on the offset field

text: I’m trying to write a code, which will find the features, which have an override on the offset field. It is very easy to find if a feature has any override, but to find if a feature has an override on specific Geometric effect probably needs absolutely different approach. Here is my code: Public Sub RepClass() On Error GoTo ErrorHandler                   Dim pMap As IMap         Dim pActiveView As IActiveView         Dim pEnumLayer As IEnumLayer         Dim pLayer As ILayer         Dim pFeatureSelection As IFeatureSelection = Nothing         Dim pFSelLayerObalkaZR As IFeatureSelection = Nothing         Dim pFeatureLayerML As IFeatureLayer = Nothing         Dim pFeatureLayerZKomOBjektB As IFeatureLayer = Nothing         Dim pFeatureClass As IFeatureClass = Nothing              pMxDocument = m_application.Document              pMap = pMxDocument.FocusMap              pActiveView = pMxDocument.FocusMap                   'Go through all layers in map document         pEnumLayer = pMap.Layers              pLayer = pEnumLayer.Next         Do Until pLayer Is Nothing             If TypeOf pLayer Is IFeatureLayer Then                 'Find the right layer                 If pLayer.Name = "Z_KomObjekt_L" Then                     'MsgBox("Tady 2 " & pLayer.Name, MsgBoxStyle.OkOnly)                     pFeatureLayerZKomOBjektB = pLayer                     'MsgBox("Tady 2 " & pLayer.Name, MsgBoxStyle.OkOnly)                     pFeatureClass = pFeatureLayerZKomOBjektB.FeatureClass                     MsgBox(pFeatureClass.AliasName & " nalezen", MsgBoxStyle.OkOnly)                     Exit Do                 End If             End If             pLayer = pEnumLayer.Next         Loop         pEnumLayer.Reset()         'Exit if not a line  feature class         If Not pFeatureClass.ShapeType = ESRI.ArcGIS.Geometry.esriGeometryType.esriGeometryPolyline Then             MsgBox("Shape type must be a Point")             Exit Sub         End If              Dim pRWE As IRepresentationWorkspaceExtension         Dim pRepClass As IRepresentationClass         Dim pRepresentationRule As IRepresentationRule = Nothing         Dim pRepresentationRules As IRepresentationRules         Dim pRepRuleInit As IRepresentationRuleInit        Dim lCircleID As Long         Dim pDiamond As ISimpleMarkerSymbol         Dim lDiamondID As Long         Dim pQueryFilter As IQueryFilter         Dim pFCursor As IFeatureCursor         Dim lSymFieldIndex As Long         Dim i As Long         Dim pFeature As IFeature              'MsgBox("OK")         pRWE = GetRepWSExtFromFClass(pFeatureClass)  'Function that returns a reference to RepresentationWorkspaceExtension              pRepClass = pRWE.OpenRepresentationClass("ZM50_Z_KomObjekt_L") 'Representation Class belonging to feature  layer         pRepresentationRules = pRepClass.RepresentationRules                   Dim pQFilter As IQueryFilter         pQFilter = New QueryFilter         pQFilter.WhereClause = ""              Dim pFeatureCursor As IFeatureCursor         Dim pMapContext As IMapContext         Dim pDisplay As IDisplayTransformation         Dim pRepresentation As IRepresentation         Dim pOverride As IOverride         Dim pGraphicAttributes As IGraphicAttributes = Nothing         Dim iIndex2 As Integer         Dim iNumberOfRepLayers As Integer         Dim pGeometricEffects As IGeometricEffects = Nothing         Dim pGeometricEffect As IGeometricEffect         Dim pBasicLineSymbol As IBasicLineSymbol         'Dim pGeometry As IGeometry = Nothing         'Dim pRepRule As IRepresentationRule = Nothing              pMapContext = New MapContext         pDisplay = pMxDocument.ActiveView.ScreenDisplay.DisplayTransformation         pMapContext.InitFromDisplay(pDisplay)         pFeatureCursor = pFeatureLayerZKomOBjektB.Search(pQFilter, False)         pFeature = pFeatureCursor.NextFeature         Do Until pFeature Is Nothing             pRepresentation = pRepClass.GetRepresentation(pFeature, pMapContext)             pOverride = pRepClass.GetRepresentation(pFeature, pMapContext)             If pOverride.HasAttributeOverride = True Then                 MsgBox("Feature " & pFeature.OID & " has an override")             Else                 MsgBox("Feature " & pFeature.OID & " hasn not an override")             End If                  AND HERE I’M LOST, I DON’T KNOW WHAT TO DO NEXT, NEXT LINES ARE ONLY SHOOTS IN THE DARK             'Find what rule a feature has             If pRepresentation.RuleID = 62 Then                 MsgBox("RuleId is 62")                 pRepresentationRule = pRepresentationRules.Rule(62)                  ElseIf pRepresentation.RuleID = 64 Then                 MsgBox("RuleId is 64")                 pRepresentationRule = pRepresentationRules.Rule(64)                 iNumberOfRepLayers = pRepresentationRule.LayerCount                 MsgBox("Number of layers in RepRule " & iNumberOfRepLayers)                      If TypeOf pRepresentationRule.Layer(0) Is IBasicLineSymbol Then                               Dim pStroke As ILineStroke                     'pBasicLineSymbol = New BasicLineSymbol                     pBasicLineSymbol = pRepresentationRule.Layer(0)                     pStroke = pBasicLineSymbol.Stroke                     pGraphicAttributes = pStroke                 End If                  End If                  pFeature = pFeatureCursor.NextFeature         Loop      ErrorHandler: MsgBox(Err.Description, vbInformation, "CreateRepresentationClass") End Sub               'GetRepWSExtFromFClass is a function that takes a reference to a FeatureClass and returns the RepresentationWorkspaceExtension      Public Function GetRepWSExtFromFClass(ByVal pFClass As IFeatureClass) As IRepresentationWorkspaceExtension 'Get RepresentationWorkspaceExtension from Feature Class On Error GoTo ErrorHandler Dim pDataset As IDataset Dim pWorkspace As IWorkspace Dim pExtManager As IWorkspaceExtensionManager Dim pUID As IUID                   pDataset = pFClass         pWorkspace = pDataset.Workspace         pExtManager = pWorkspace         pUID = New UID         pUID.Value = "{FD05270A-8E0B-4823-9DEE-F149347C32B6}"         GetRepWSExtFromFClass = pExtManager.FindExtension(pUID)              Exit Function      ErrorHandler: MsgBox(Err.Description, vbInformation, "GetRepWSExtFromFClass") End Function


How to set attributes for writers generic objects, for instance, lines?! Writer generic objects lines are created, but without attribute because I do not know how to set attributes and values.


title: QGIS 2.0 Plug Repository Connection Failed

text: I'm trying to connect to the QGIS Official Plugin Repository but the status within the Plugin Manager window is always"unavailable" for the URL - http://plugins.qgis.org/plugins/plugins.xml. I'm running QGIS 2.0 within Windows 7. Could the problem be that QGIS is unable to detect the proxy settings used on the corporate network?


title: How to view Esri GRID format file in wxpython file browse window?

text: I'm trying to browse for all kinds of vector and raster data using a filebrowsebutton. The problem is I don't know how to let the browse window see the Esri **GRID** files (which contains a folder and a bunch of other files in a Windows standard file browser window) as a **single** file just as the ArcGIS tools see it, so that I can load it in by clicking on it. I searched through the internet but didn't find any discussion on it. I know how to write code to access a specific GRID file, but not through a file browser. Can anyone help? Thanks a lot. Below are my codes about the filebrowsebutton:               wildcard = "Shapefiles (.shp)|*.shp|" \                    "Raster files (GRID, .img, .tif, .bmp, .gif, .jp2, .jpg, .png)|*; *.img; *.tif; *.bmp; *.gif; *.jp2; *.jpg; *.png|" \                    "All files (*.*)|*.*"              self.fileBrowseButton1 = wx.lib.filebrowsebutton.FileBrowseButton(buttonText='Browse',               dialogTitle='Choose a file', fileMask= wildcard,               id=wxID_FRAME1FILEBROWSEBUTTON1, initialValue='',               labelText=u'Locate the GIS dataset:', parent=self,               pos=wx.Point(40, 80), size=wx.Size(496, 40), startDirectory=u'.',               style=wx.TAB_TRAVERSAL,               toolTip='Type filename or click browse to choose file')


title: Does ModelBuilder need to use variable to define projection and then re-project?

text: I am trying to build a Model in ModelBuilder to create two different process raster datasets:   1. One set is to define a raster dataset to set a WGS 84 projection; and then   2. the second set is to take the define projection raster datasets to reprojection to NAD UTM 83 to a new folder. Do I need to create a variable to copy from a define projection raster set and then Reprojection raster dataset to NAD UTM 83? What is the correct variable tool should I use to perform to connect them? BTW, I have to do this with a lot of raster datasets and do the same thing.


title: Syntax for wfs-t request

text: I am in need of the correct syntax for performing transaction on WFS layer using OpenLayers. I have gone through many examples on the OpenLayers site and other examples, but can't figure out how to use the necessary and optional options like geometryname, schema etc


title: Requires Datatype for inserting features

text: I have imported postgis data in geoserver. In geoserver the spatial data column is displayed as MultiLineString, Does it affect when inserting features into geoserver or does it have to be explicitly defined as geometry?


title: Why are points collected with a Trimble FM 750 off by more than a mile when displayed in QGIS?

text: I have a grower that has a FM 750 and is trying to make points in a field. When finished he takes the file out of the display and goes to his computer and he can bring up the points in the right field using Google Maps. When he tries to lay the points in QGIS the points are more then a mile off from his field or show up in another field. Guess what I'm looking for is a support number to get a hold of to visit about this or if the display and the QGIS systems are not compatible.


title: Counting slippy map tiles

text: I would like to calculate approximately, how many tiles of given size would need to be rendered for a given bounding box and zoom level range in a given projection (always 3857). For example, for:   * Tilesize: 256x256 pixels   * Longitude: 15-25°   * Latitude 10-20°   * Zoom levels 5-15 I know that number of tiles at zoom level nx = 4^x. I think the procedure goes like this: add 1 tile for every level in which one tile is still enough to cover the entire bounding box. For every level further down, assume tiles scale based on the equation above. Therefore my question is: How do I know the last zoom level on which I can cover the bounding box with 1 tile?


title: How to remove columns in a spatial polygon data frame in R?

text: My spatial polygon data frame (spdf) contains too many columns (variables) and I want to remove most of the columns entirely. I know how to do this with a regular data frame in R, but I am unsure how to do this when dealing with spdf.


title: How you can repair invalid topology errors within postgis?

text: I've imported the shapefile of corine land cover project to my postgis database and the               select gid, ST_IsValidReason(geom) from gis.clc2000 where not ST_IsValid(geom);      returns in every occasion               gid    |                                st_isvalidreason                                     -------+--------------------------------------------------------------------------------          1 | IllegalArgumentException: Points of LinearRing do not form a closed linestring     .....      I know that my polygons are closed by visually instecting some of the errors. How can I fix the topologies errors from withing the server? edit: I'm using postgresql/postgis2.0.1, Ubuntu12.04


title: Is there a publicly available instance of ArcGIS Server

text: Are there any publicly available instances of ArcGIS Server for learning? I have searched Amazon EC2 but I don't see any public instances available (unless you already have an ArcGIS Server license).


title: How do I combine many ECW files in MapInfo?

text: I have MapInfo 11.5.0.16 (I don't have ArcGIS) and I am trying to combine (in some sense) a large collection of ECW files. I can load one ECW at a time in MapInfo but I have very many. MapInfo creates one TAB file per image, and I need to manually choose a projection for each image. This is time consuming and prone to human error and I would need to open each tab file, which creates one layer per image. I would like to open one TAB file in MapInfo which creates one layer in a map, where that layer will display all the ECW files. I've reviewed the similar questions and Google but I can't seem to find a solution that works for MapInfo. Thanks for your assistance!


title: Using QGIS calculate distance of point along a line

text: I have lots of points, and a road (lines) and I would like to calculate from point A to point B (the distance must route through the roads (lines). I found this QGIS calculate distance of point along a line , but when I'm trying to do a v.distance, I don't see it saving the distance into the DB. I need the DB to show:               Source      Destination       Distance      Point A     Point B          10.0011      Point A     Point C           101.12      point A     Point D            1212       Point B     point c           121.0      point B     Point D         31.0      ...         ...                ...      Any idea? Thanks!


title: Is there a way to load a kml raster super-overlay in QGIS?

text: I have created a set of tiled rasters as a super-overlay in Google Earth, and it performs remarkably well. Is there some way to load a super-overlay as a layer in QGIS? I am wondering if this could be an option to manage large numbers of rasters without setting up a WMS.


title: WPSClient with raster input parameter

text: I'm trying to run a WPS process on a GeoServer instance using OpenLayers WPSClient. One of the parameters of the process is a raster layer. How do I pass this to using the WPSClient? The documentation indicates how to call this with vector data, but I can't seem to figure out how to do it with raster data.


title: How to create an extent indicator in QGIS 1.8 print composer?

text: I'm looking to put an inset map for my main map frame onto my map layout in QGIS print composer. I would like it to indicate the extent and location of the main map on a larger regional map. This something I am familiar with accomplishing in ARCMap but I have yet to discover if QGIS has an equivalent function. The indicator should show the extent of the data in the main frame in another data frame.


title: How to install PostGis on EC2 micro-instance

text: I tried to install PostGIS on my Amazon EC2 microinstance using `sudo yum install postgis2` and got following error message:               $ sudo yum install postgis2                             Loaded plugins: priorities, security, update-motd, upgrade-helper     amzn-main                                                | 2.1 kB     00:00          amzn-updates                                             | 2.3 kB     00:00          Setting up Install Process     No package postgis2 available.     Error: Nothing to do      How can I install PostGIS on an EC2 microinstance?


title: Any QGIS 2.0 updated RPM for Fedora?

text: I'm very excited about the new release of QGIS (2.0), and would like to install it on my Fedora 18 machine. I couldn't find any relevant RPMs, and the documentation on the QGIS main site for Fedora is pretty old (16/17). Did anyone did it and can throw some tips? Thanks!


title: GeoNames query for Canadian parks is not working

text: > **Possible Duplicate:**   >  Geonames querying for parks in Puerto Rico or US Virgin Islands Hi, I am trying this query http://api.geonames.org/search?featureCode=PPL&username=genadinik&country=CA&style=full&adminCode1=AL&maxRows=10 I was trying to get parks in Alberta, Canada. But I got no results. I am not sure where to get the list of abbreviations for Canada. How would I go about finding the state abbreviations that GeoNames requires for any country? Thanks!


title: How to work with Multivalues in arcgis script?

text: I am really beginner in arcgis scripting. I was unable to find a solution of my problem, so I'm writing here. My problem: On my input I want some shapefiles. When I load them from ArcGis Tool, I obtain something like this:               'C:\Users\Steve\Disk Google\BP-Kaiser\In\Cesta_Clip.shp';     'C:\Users\Steve\Disk Google\BP-Kaiser\In\Pesina_Clip.shp';     'C:\Users\Steve\Disk Google\BP-Kaiser\In\SilniceDalnice_Clip.shp'     'C:/Users/Steve/Disk Google/BP-Kaiser/In/Cesta_Clip.shp'      I want to have these shapefiles in a list, so I convert it using split() function to look like this:               list = ["'C:/Users/Steve/Disk Google/BP-Kaiser/In/Cesta_Clip.shp'", "'C:/Users/Steve/Disk Google/BP-Kaiser/In/Pesina_Clip.shp'", "'C:/Users/Steve/Disk Google/BP-Kaiser/In/SilniceDalnice_Clip.shp'"]      And then I would like to work with it. But if I use just a simple function:               arcpy.CopyFeatures_management(list[0], outWorkspace+"/copy.shp")      it doesn't work. But if I enter the path to the file directly, it works. I think that it has to be something simple, but like I said I am really beginner. If you wold like to help me, but you will need more information, write me please.


title: Handling points in openlayers

text: I have a several hundred points that I am showing on a map. At the moment, my method of adding the points seems inefficient. I create an array using the following approach:               var features = [new OpenLayers.Feature.Vector(OpenLayers.Geometry.fromWKT("POINT(-4.2522971 45.3447353)").transform(new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject()),{"Name":"aaa"}),     new OpenLayers.Feature.Vector(OpenLayers.Geometry.fromWKT("POINT(-4.259215 45.344827)").transform(new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject()),{"Name":"bbb"})]      and then add this array to a map-layer:               var vectors = new OpenLayers.Layer.Vector("vector", {isBaseLayer: false, projection:"EPSG:4326", styleMap:myStyle});              vectors.addFeatures(features);     map.addLayer(vectors);      Is there are a more efficient way of displaying the points? I also want to be able to dynamically change the radius of each point separately. I'm not sure if that is something that I should consider when constructing the point. For example, I can choose the radius when constructing, but I do not know how to update it:               new OpenLayers.Feature.Vector(OpenLayers.Geometry.fromWKT("POINT(-4.2522971 45.3447353)").transform(new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject()),{"Name":"aaa"}, {radius: 10})


title: How to create an interactive map with jQuery and HTML?

text: I want to develop a map in JQuery and HTML. It needs to show world country borders. It must be enabled mouse hover effect and also popup (when mouse over) in which I can add my own text. What JS plugin do you suggest? I tried JVectorMap but it has only hover effect and I couldn't make popup window with own text. Thank you.


title: How to add labels outside of features?

text: I'm trying to achieve this kind of labeling using QGIS. Apparently, it doesn't seem to be easily supported. Would you have an idea for a workaround ? ![enter image description here](http://i.stack.imgur.com/e2tcm.png) I managed to have circled text using "Diagrams", which it seems I can even place using an X and Y coordinate field ! But I have no idea on how to draw a line between the circle and the building... I have about 400 circles to place, and they could quite often change, so I can't rely on a manual method. Any idea ? Thanks a lot ! * * * So here's what I've got using the plugin you suggested : ![what i've got](http://i.stack.imgur.com/HET1v.png) Not bad at all !! Thank you VERY much ! I'm only struggling to put the number in the center of the circle. I see there are two fields : LblAlighH and a LblAlignV. Entering "center" in LblAlignH centers the number horizontally, but I tried "middle" or "center" in LblAlignV without anything happening... (in the image, I've entered center for LblAlignH in the one at the bottom) Would that be possible ? :) (Are you by any chance the developer of that plugin ?) * * * Okay, "half" in the VAlign field seems to work fine ;) To create the circles, I used line decorations in the Arrow_Label_XXX layer : ![how to](http://i.stack.imgur.com/aCX4W.png) And here's the final result : ![final result](http://i.stack.imgur.com/dftwz.png)


title: Creating JP2 rasters with FWTools

text: I have been trying to make JP2 raster files with FWTools with limited success. FWTools will not allow for the creation of a JP2 file at a specific file size. However, I am able to make a JP2 with the same data within ArcMap. Is there a file size limit with FWTools? Does anyone know of any work arounds? I am looking to make larger sized GeoTiffs (100+ MB) into JP2 files. Thanks...


title: Is there a Log or Ln function in QGIS Field Calculator?

text: I can't seem to find any way to take the logarithm (natural or base 10) of a field in QGIS. I have gotten to the field calculator but there is not a button for log and I can't seem find any formula that works. I'm using 1.4.0 on linux currently. I couldn't find anything in release notes about this being added since, but does anyone know if it's in a newer release? Many thanks!


title: Printing using arcgis runtime java

text: I have a map viewer that was developed using the arcruntime java sdk. I am trying to make a print button and it works with one issue. I cant seem to figure out how to make it fit on one page.                   private void printMap() {              Paper paper = new Paper();              // ...set up the paper options...         PrinterJob printJob1 = PrinterJob.getPrinterJob();          PageFormat pageFormat = new PageFormat();         HashPrintRequestAttributeSet set = new HashPrintRequestAttributeSet();                       if (printJob1.printDialog(set)) {                 try {printJob1.print(set);}                 catch (PrinterException exc) {                    // System.out.println(exc);                  }              }                 //set up page format [DEFAULT] options and pass the paper settings through         pageFormat.setPaper(paper);              // set page format for the map to print         map.setPageFormat(pageFormat);         printJob1.setPageable(map);              try {            // send job to printer!            printJob1.print();                 // new PrintMe().print(YourJTextComponent);              } catch (PrinterException e) {            // handle the exception         }          }      Im using Eclipse Kepler SR1, resolution is 1920x1080. any ideas or suggestions? This is how it prints ![enter image description here](http://i.stack.imgur.com/OXp2W.png)


title: Zonal statistics vs. random points?

text: Using ArcGIS 10.1 for Desktop, I wish to model the relation between soil productivity (SP) and a number of independent variables such as, e.g., terrain elevation and slope. For this purpose, I am going to use logistic regression with elevation and slope (from DEM) as predictors and SP as dependent categorical outcome. As for SP, I got this information from historical maps (early 1900s) depicting parcels of land that have been given different colors according to SP. I have digitized those parcels, coming up with different polygons (with the SP class entered as attribute). Now, my question is which of the two approaches below would be more sound, and/or if anyone of you has tackled a similar issue: is it better   1. to use 'zonal statistics as table' to get the mean elevation and slope value for each parcel,   2. to get slope and elevation value from random points within each SP class parcel. May be I am over-thinking or I am missing something. My concern with approach (1) is that since some of the parcels are quite huge in size (up to 12 Hectares or 120,000 sq m) I fear that the average slope and elevation for the whole polygon could not adequately capture the within-parcel variability of those variables.


title: Adding customised data into geotiff

text: Not too familiar with the geotiff format and hence would like to raise a query on whether I could embed additional information data into geotiff metadata section? Is the fields type for geotiff fixed? or is there any allowance for customised propriteary data to be included? Thanks


title: Qgis ver.2.2 can not load osm file

text: I have download shapefile osm with vector >> osm >> downdload data in qgis ver. 2.2. When i I open osm data,qgis not responding. I hope someone could help me. Thanks ![enter image description here](http://i.stack.imgur.com/L03rc.png)


title: How can I create Albers projection of Russia map?

text: I need a conic tranformation of map, so this code               var map = L.map('map', {     crs: new L.CRS.proj4js('EPSG:4326',         '+proj=aea +lat_1=52 +lat_2=64 +lat_0=0 +lon_0=105 +x_0=18500000 +y_0=0 +ellps=krass +units=m +towgs84=28,-130,-95,0,0,0,0 +no_defs',         {           resolutions: [8192, 4096, 2048],         }       ),       continuousWorld: true,       worldCopyJump: false     });     L.tileLayer('http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png').addTo(map);      doesn't work. I got the error:               TypeError: this._scales is undefined     return this._scales[zoom];      What proper `EPSG` and `proj4def` do I need to use?


title: Problem with Adding Points ArcGIS Android

text: it's my code debug show noting and i can't see my point anywhere, is it a problem with this code?               public void addthepoint()       {               map = ((MapView)findViewById(R.id.map));          this.pointLayer.removeAll();          this.pointLayer = new GraphicsLayer();             this.map.addLayer(this.pointLayer);               //   HashMap localHashMap = new HashMap();            // localHashMap.put("name", "ENja");          double lat = Double.parseDouble("**.8854587");             double lon = Double.parseDouble("**.2812024");             mapPnt = new Point(lon, lat);             Graphic graphic = new Graphic(mapPnt, PointMarker());             pointLayer.addGraphic(graphic);       }      and this is the point symbol function. please check it and let me know what's wrong.               public SimpleMarkerSymbol PointMarker()          {            return new SimpleMarkerSymbol(Color.GREEN, 35, SimpleMarkerSymbol.STYLE.DIAMOND);          }


title: How can I assign Raster2 pixel value to Raster1 where they overlay?

text: I want to assign Raster2 pixel value to Raster1(Bigger) where they overlay. Both of them have the same cell size. How can I accomplish this process in ArcGIS Desktop? (Like Raster Calculator ) Don't know how. Thanks !!


title: QGIS 2.2 Raster Merge losing color data

text: I'm trying to merge a couple of GeoTiff maps using Raster > Miscellaneous > Merge but the merge is only coming out in grayscale. The geotifs are using palettes color. How do I get the merge to keep the color information?


title: What is the easternmost point in the US?

text: In his blog for The Atlantic magazine, James Fallows claimed that Eastport, ME was the easternmost point in the US. A writer commented that no, the easternmost point in the US was somewhere in Alaska, as Alaska crosses the 180th meridian. Another correspondent (me) wrote that this was silly - that the easternmost point of a landmass is the extreme point of the landmass that can be reached by traveling in an easterly direction, or more precisely, the last point on the country found by sweeping a line of longitude in an easterly direction. (This is more complicated for countries that don't have contiguous territory, but for now let's assume the territory is contiguous (and yes, I know that Alaska borders Canada and not another US state)). Is there any convention in the geography community that establishes the meaning of "easternmost", preferably one that can be cited?


title: Create a polygon shapefile from .dbf table

text: I am using arcgis 10.1 and have a table "dbf" containing information of polygons. This table has fields "OID" "SHAPE" "SHAPE.area" "SHAPE.LEN" there is a way or tool to convert this table as shp polinos


title: ORA-04043: object ST_GEOMETRY_SHAPELIB_PKG does not exist

text: I uninstalled ArcSDE using the esri.exe, and dropped user sde through running the sql statement "drop user sde cascade". Now I want to reinstall ArcSDE. The installation runs successfully, but when I run post installation configuration to "create SDE repository", it failed. Here is the error in sde_setup.log: > [Thu Jul 24 22:19:55 2014] SQL Stmt: [Thu Jul 24 22:19:55 2014] ERROR > setting up st_geometry type. Error: -51 > > [Thu Jul 24 22:19:55 2014] DBMS error code: 4043 ORA-04043: object > ST_GEOMETRY_SHAPELIB_PKG does not exist > > [Thu Jul 24 22:19:55 2014] ERROR setting up St_Geometry, Error = -1 How can I fix this error? I'm using Oracle 11g release 2 as the underlying database, and my ArcSDE version is 10.0


title: Copying feature layers

text: I am trying to copy feature layers using the following code, but get no results. Can anyone help me understand why the code does not work?               import os     import arcpy      arcpy.env.workspace = ws = 'C:\\WorkSpace\\sahie_saipe\\sahie_saipe.gdb'     dictionary_Layer = {'co06_Layer':'co06', 'co07_Layer':'co07', 'cnty08_Layer':'co08', 'co09_Layer':'co09', 'co10_Layer':'co10', 'co11_Layer':'co11', 'co12_Layer':'co12', 'st_Select':'state'}     for layer in arcpy.mapping.ListLayers(mxd):         if layer in dictionary_Layer:             old = os.path.join(ws, layer)             new = os.path.join(ws, dictionary_Layer[layer])             try:                 arcpy.CopyFeatures_management(old,new,"")                 print 'copy %s to %s' %(layer, dictionary_Layer[layer])             except:                 print 'Could not copy %s'


title: Is it possible to create a join between a feature class and a table created with an InMemoryWorkspaceFactory

text: I'm trying to join a featureclass contained inside a pgdb with a table created with an InMemory Workspace. The in-memory table contains results that I want to use to define symbology on the featureclass. The creation of the table goes well, I'm able to see the calculated data, but the IRelationshipClass created returns no records. I tried to use multiple fields to join the data without success. I also added the created table to ArcMap to create the join directcly in ArcMap, without success. It passes all the validations, but returns no records. I tried to join my inmemory table with a shapefile featureclass instead and it worked. Here's a sample of code                           // Create an in-memory workspace factory.                 IWorkspaceFactoryPtr ipWorkspaceFactory(CLSID_InMemoryWorkspaceFactory);                      // Create a new in-memory workspace. This returns a name object.                 IWorkspaceNamePtr ipWorkspaceName;                 ipWorkspaceFactory->Create(NULL, L"ValeursCalculees", NULL, 0, &ipWorkspaceName);                      INamePtr ipName(ipWorkspaceName);                      // Open the workspace through the name object.                 IUnknownPtr ipWorkspaceUnk;                 ipName->Open(&ipWorkspaceUnk);                 IWorkspacePtr ipWorkspace(ipWorkspaceUnk);                      IFeatureWorkspacePtr ipMemoryFeatureWorkspace(ipWorkspace);                      IDatasetPtr ipFeatureClassDataset(ipFeatureClass);                 CComBSTR strNomFC;                 ipFeatureClassDataset->get_Name(&strNomFC);                      IFieldsEditPtr ipFieldsEdit(CLSID_Fields);                                ... Field creation code ...                      ITablePtr ipTable;                 ipMemoryFeatureWorkspace->CreateTable("InMemoryTable", ipFieldsEdit, NULL, NULL, NULL, &ipTable);                               ... Table Filling code ...                      //IRelQueryTable                 IFeatureWorkspacePtr ipFeatureWorkspace(CTools::GetWorkspaceFromFeatureClass(ipFeatureClass));                 IObjectClassPtr ipObjectClass(ipTable);                      IMemoryRelationshipClassFactoryPtr ipMemoryRelationshipClassFactory(CLSID_MemoryRelationshipClassFactory);                 IRelationshipClassPtr ipRelClass;                 HRESULT hr = ipMemoryRelationshipClassFactory->Open(L"Join", ipFeatureClass, L"Numero_AD", ipObjectClass, L"NumeroAD", L"forward", L"backward", esriRelCardinalityOneToMany, &ipRelClass);                      IDisplayRelationshipClassPtr ipDisplayRelationshipClass(ipFeatureLayer);                 ipDisplayRelationshipClass->DisplayRelationshipClass(ipRelClass, esriLeftOuterJoin);                      //Add the table to ArcMap                 IStandaloneTableCollectionPtr ipStandaloneTableCollection(CTools::GetMap(ipMxDocument));                 IStandaloneTablePtr ipStandaloneTable(CLSID_StandaloneTable);                 ipStandaloneTable->putref_Table(ipTable);                 hr = ipStandaloneTableCollection->AddStandaloneTable(ipStandaloneTable);                 ipMxDocument->UpdateContents();      Am-I trying to do something impossible?


title: Create symbology from multiple layers

text: I have four different polygons of soil data that encompass overlapping geographical areas. There is also some overlap of soil series, but each has a few unique ones:   * Polygon A has Series A, B, C   * Polygon B has Series A, C, D   * Polygon C has Series A, B, D, E   * Polygon D has Series A, E, F, G Right now, each of them is symbolized based only on that polygon, meaning I've got different colors for the same soils in different maps. ![Mismatched colors in legend](http://i.stack.imgur.com/JQtxN.png) My goal is consistent symbology for all the layers so each series always has the same color. The end product is about eight maps, and I need the same representation for all. Since no layer encompasses the entire soil range, however, I can't copy symbology from any layer to the others without some soils being absent. To make it even more annoying, the different polygons have different naming structures for the soil types (some are abbreviated, some are broken down into additional detail [e.g. sandy clay loam vs. sandy loam], etc.). Is there any way in ArcMap (10.1) to create a symbology from multiple layers/features at once, to cover all possible categories A though G?


title: PostGIS non functional on Ubuntu 14.04

text: I am unable to create the PostGIS extension on new databases in Ubuntu 14.04:               $  psql pg_testing -c "create extension postgis"     ERROR:  could not load library "/usr/lib/postgresql/9.3/lib/postgis-2.1.so": /usr/lib/liblwgeom-2.1.3.so: undefined symbol: GEOSDelaunayTriangulation      All the software was installed from the Postgres repositories:               $ dpkg -l | grep postgres     ii  libpostgresql-jdbc-java                               9.2-1002-1                                          all          Java database (JDBC) driver for PostgreSQL     ii  pgdg-keyring                                          2013.2                                              all          keyring for apt.postgresql.org     ii  postgresql                                            9.3+154.pgdg14.04+1                                 all          object-relational SQL database (supported version)     ii  postgresql-9.3                                        9.3.4-1.pgdg14.04+1                                 amd64        object-relational SQL database, version 9.3 server     ii  postgresql-9.3-dbg                                    9.3.4-1.pgdg14.04+1                                 amd64        debug symbols for postgresql-9.3     ii  postgresql-9.3-postgis-2.1                            2.1.3+dfsg-3.pgdg14.04+2                            amd64        Geographic objects support for PostgreSQL 9.3     ii  postgresql-9.3-postgis-scripts                        2.1.3+dfsg-3.pgdg14.04+2                            all          Geographic objects support for PostgreSQL 9.3 -- scripts     ii  postgresql-client                                     9.3+154.pgdg14.04+1                                 all          front-end programs for PostgreSQL (supported version)     ii  postgresql-client-9.1                                 9.1.13-1.pgdg14.04+1                                amd64        front-end programs for PostgreSQL 9.1     ii  postgresql-client-9.3                                 9.3.4-1.pgdg14.04+1                                 amd64        front-end programs for PostgreSQL 9.3     ii  postgresql-client-common                              154.pgdg14.04+1                                     all          manager for multiple PostgreSQL client versions     ii  postgresql-common                                     154.pgdg14.04+1                                     all          PostgreSQL database-cluster manager     ii  postgresql-contrib                                    9.3+154.pgdg14.04+1                                 all          additional facilities for PostgreSQL (supported version)     ii  postgresql-contrib-9.3                                9.3.4-1.pgdg14.04+1                                 amd64        additional facilities for PostgreSQL     ii  postgresql-server-dev-9.3                             9.3.4-1.pgdg14.04+1                                 amd64        development files for PostgreSQL 9.3 server-side programming          $ apt-cache policy postgresql-9.3     postgresql-9.3:       Installed: 9.3.4-1.pgdg14.04+1       Candidate: 9.3.4-1.pgdg14.04+1       Version table:      *** 9.3.4-1.pgdg14.04+1 0             500 http://apt.postgresql.org/pub/repos/apt/ trusty-pgdg/main amd64 Packages             100 /var/lib/dpkg/status          9.3.4-1 0             500 http://de.archive.ubuntu.com/ubuntu/ trusty/main amd64 Packages          $ apt-cache policy postgresql-9.3-postgis-2.1     postgresql-9.3-postgis-2.1:       Installed: 2.1.3+dfsg-3.pgdg14.04+2       Candidate: 2.1.3+dfsg-3.pgdg14.04+2       Version table:      *** 2.1.3+dfsg-3.pgdg14.04+2 0             500 http://apt.postgresql.org/pub/repos/apt/ trusty-pgdg/main amd64 Packages             100 /var/lib/dpkg/status          2.1.2+dfsg-2 0             500 http://de.archive.ubuntu.com/ubuntu/ trusty/universe amd64 Packages      And regarding GEOS:               $ dpkg -l | grep libgeos     ii  libgeos-3.4.2                                         3.4.2-4ubuntu1                                      amd64        Geometry engine for Geographic Information Systems - C++ Library     ii  libgeos-c1                                            3.4.2-4ubuntu1                                      amd64        Geometry engine for Geographic Information Systems - C Library     ii  libgeos-dev                                           3.4.2-4ubuntu1                                      amd64        Geometry engine for GIS - Development files          $ apt-cache policy libgeos-3.4.2     libgeos-3.4.2:       Installed: 3.4.2-4ubuntu1       Candidate: 3.4.2-4ubuntu1       Version table:      *** 3.4.2-4ubuntu1 0             500 http://de.archive.ubuntu.com/ubuntu/ trusty/universe amd64 Packages             100 /var/lib/dpkg/status      My guess is that PostGIS is expecting a different version of GEOS. The question is which? And where from to install it? **Update** : There are some hints of this issue being caused by an outdated version of `liblwgeom`. I thus removed all related packages (`apt-get purge liblwgeom*`) and re-installed PostGIS. This is what I have now:               $ dpkg -l | grep liblwgeom     ii  liblwgeom-2.1.3                                       2.1.3+dfsg-3.pgdg14.04+2                            amd64        PostGIS "Lightweight Geometry" library      But the problem remains.


title: Moving Esri grid format and info folder inside working directory?

text: Using the below code, I want to extract a new raster from an original raster using value field and I want to automate the process for 70 more rasters. However, for testing purpose, I'm using only 3 rasters. Out of 3 rasters inside the working directory, this code works only for one raster which has both " ESRI grid" and "info" folder inside the working directory. I was not able to keep the "info" folders of other two raster inside the working directory as the "info" folder of 1st raster is already in the working directory. Could you please let me know how can I keep the "info" folders of other two rasters inside the working directory, so that code will work and I can automate the process? Code:               import arcpy, os     from arcpy import env     from arcpy.sa import *          #To overwrite output     arcpy.env.overwriteOutput = True          #Set environment settings     env.workspace = "C:/Subhasis/Test/Neshanic_Python"          outws="C:/Subhasis/Test/Neshanic_Python/extract"          #checkout ArcGIS spatial analyst extension license     arcpy.CheckOutExtension("Spatial")          inraster = arcpy.ListRasters("*", "GRID")          for i in inraster:         flds = ("VALUE", "COUNT")           dct = {row[0]:row[1] for row in arcpy.da.SearchCursor(i, flds)}           sumcnt = sum(dct.values())         dct1 = {k:v for (k,v) in dct.items() if k >= 15}         sumcnt1 = sum(dct1.values())         percentage=(float(sumcnt1)/float(sumcnt))         print percentage         newraster = ExtractByAttributes(str(i), "VALUE>=15")         outname=os.path.join(outws,str(i))         newraster.save(outname)


title: Arcpy (10.0) and centroids/points shapefile

text: I'm trying to create a point shapefile from a grid index polygon shapefile I created using arcpy. I know this is simple, but could anyone give a few pointers out there. Here is what I've done now: but my error is: Traceback (most recent call last): File "C:\Users\RLong\Documents\NewPoint.py", line 19, in arcpy.FeaturetoPoint_management("GridIndex", "SamplesLoc", "Centroid") AttributeError: 'module' object has no attribute 'FeaturetoPoint_management' Any ideas?               ###Simple point shape file object for testing     ###     ###Import standard modules     import os, sys, math, arcpy, traceback, string, fileinput     ###Setup your workspace or file you're going to use     from arcpy import env     arcpy.env.workspace = "C:/Users/RLong/GIS_Data/Conversions/GridIndex.shp"     ###setup local variables     in_Features = "GridIndex.shp"     point_location = "Centroid"     out_feature_class = "SamplesLoc"          ###execute point     arcpy.FeaturetoPoint_management("GridIndex", "SamplesLoc", "Centroid")


title: Current time from GPRMC or GPZDA?

text: In `GPRMC` there is a time field, called `fix time`. Can I consider this always the same as value from `GPZDA` time? Is there any possibility that when GPS device is configured to stream both `GPRMC` and `GPZDA` sentences, but their value will be different? For my task it would be great if I could get time from `GPRMC`, however, I'm not sure if I can get consistent time value from there.


title: Techniques for offline reverse geocoding on a mobile device?

text: I am working on a mobile mapping application (currently iOS, eventually Android) - and I am struggling with how to best support reverse geocoding from lat/long to Country/State without using an online service. Apple's reverse GeoCoding API depends on Google as the backend, and works great while connected. I could achieve similar functionality using the Open Street Maps project too. What I really want however is to create a C library that I can call even when offline from within my application, passing in the GPS coordinates, and having it return the country and/or state at those coordinates. I've seen examples of how to do this on a server, but never anything appropriate for a mobile device. What are some techniques to accomplish this?


title: Algorithm to move polygon nodes when simplifying

text: I want to simplify polygons that come from a raster analysis. In order to avoid overlap and gaps, I first convert my polygons to lines. The problem is that all algorithm are then keeping the two end points, which are not necessarily at the best position. The image below illustrates my problem (it is difficult to explain by words). The "manual" case is the kind of result that I expect (of course I cannot do this manually for all features because I have far to many of them):   * With Wang's algorithm, I have some bends around the node and it is away from the straight line that I would like   * With DP algorithm, I have the straight lines but they are shifted compared with the line that I would like.  My priority is to keep a good positional accuracy (minimizing the distance from the "middle" of the line). I would like a "T"(though not necessarily right angles) like I have with DP algorithm, but at a location that would make the lines close to the "middle" of the original lines. I prefer solutions in QGIS 2.0 or ArcGIS 10.1 (advanced) but the implementation is a bonus for me. My question is about an algorithm that would help me find the best location of the polygon nodes for "non shifted straight lines". ![enter image description here](http://i.stack.imgur.com/BDlb5.jpg)


title: gdalwarp outline artifact

text: I'm using gdal to warp images into a different map-projections. However the warping processs seem to leave grey marks where the boundary of the original image was (see warped.png below). I figured out this was probably caused by my image having a white background and the 'filler' pixels created by the warping process being transparent but black. And the interpolation is smudging the black and white together (even though it shouldnt as the black is transparent). Is there any way to avoid the grey outline from appearing? Maybe force gdalwarp to fill with opaque white rather than transparent black? I cant do a square crop on the warped image as I'll loose bits of the image, and I cant just make my background black instead of white as this messes up some filtering I have to do beforehand. Theres seems to be lots of options for gdalwarp that sound useful like -srcnodata but I cant figure out how to use them. Thanks input.png warped.png warped.tiff


title: Restricting access to geoserver module

text: I built an extension module for Geoserver 2.5 using Maven and Wicket. It places a Label to the left of the navigation list that calls a HTML-Page in Geoserver when clicked. It looks and works as it should; the problem is that it is visible and usable without login. How can I configure my extension module/GeoServer so that it gets to the plugin only after login?


title: How to obtain overlapping semi-transparent circles effect?

text: I have been trying to obtain the following markers effect. How can I do the same in ArcGIS (or any free software)? This image is from a book.![alt text](http://i.stack.imgur.com/CZSop.jpg)


title: How to add an online basemap in ArcGIS 10

text: How can I add an external online map in ArcMap by inserting the http address, to use as a basemap in order to georeference multiple aerial photos? Or is there a way to georeference these photos using ArcGlobe? P.S.: I don't really care about the precision of the georeferencing. I just want to use the template as a quick-to-use index for the aerial photos. P.S.2: There is an option in ERDAS that you can add any external map to use as a reference in LPS. Thank you in advance!


title: How to set different map layers for different zoom levels using ArcGIS JavaScript API?

text: I have created a layer using "clip" tool in ArcMap. ![enter image description here](http://i.stack.imgur.com/SO3cO.jpg) I have another layer in TOC too. whichis the same area in "clip" Layer with more polygons. ![enter image description here](http://i.stack.imgur.com/zhhnw.jpg) In my web map application (javaScript API) Is there any way to set these layers in a way that when I zoom in I get more details like image2 and by zooming out I get polygons in image1.


title: Postgres/PostGIS Spatial Query Performance Problem

text: I have a query in Postgres/PostGIS that I am a bit confused about. The spatial intersection happens very fast, ~ 2 seconds for 10,000, but the update takes > 15 seconds for the same number of points. a_table contains ~ 66,000,000 points b_table contains ~ 3,000 polygons Both are spatially indexed (gist and clustered) and the b.wanted_value column is indexed as well. Any suggestions for improving the speed of this query would be sincerely appreciated.               UPDATE a_table a     SET a.needed_value = (       SELECT b.wanted_value        FROM b_table b        WHERE ST_Intersects(a.the_geom, b.the_geom)     );


title: Join feature classes from two ArcSDE geodatabases with same schema?

text: I am working with two ArcSDE geodatabases -- one for test, one for production. They both have the same schema. I recalculated an attribute in the test environment and wanted to push that to production. My instinct was to join the two feature classes and use the Field Calculator, but there was no way to distinguish which column belonged to which database. So I appended `_TEST` to the aliases of the test class. When I ran the join again, it automatically appended `_TEST` to the _production_ aliases as well. [grunt] Somehow I figured out which was which and was ready to use Field Calculator but... the option was grayed out in the right-click menu for the column. All columns, in fact. So my main question is: How is it possible to join feature classes from two ArcSDE geodatabases and copy attribute values from one to the other? Secondly, why on earth were my production field aliases being suffixed with `_TEST`? It seems like ArcMap has no way of managing a join between feature classes with the same field names. Is this true?


title: How to make background of EMF symbol transparent in ArcGIS?

text: I have been sent a collection of PNG and SVG images with which I should symbolize my points in ArcGIS 10.0. Unfortunately ArcGIS 10.0 doesn´t support PNG when symbolizing points, so I have converted them to EMF files. The background transparency of the images are now lost and are displayed black. How can I make the background of these emf images transparent in ArcGIS?


title: How to create a custom data viewer on top of QGIS?

text: I was exploring the possibility of creating a simple Data Viewer based on QGIS. I would like to create an simple application that will open QGIS up with a limited set of tools & menus. QGIS would open up with a given Project file and data (in shapefile format). Preferably, I'll like to make a stand alone installation program, that will configure QGIS and copy the data like I have mentioned above. Is something like this possible?


I tried to find some online platform-independent GIS services that could be called by any programming language and doing various GIS operations - like intersection, data format conversion, transformation and so on. WPS standard is a hot candidate for that (like in case of wps-grass-bridge) but unfortunately I cannot find such (or similar) services online - published by some provider. Most of WPS that is possible to find are just for testing purposes or are not working anymore (like from this post: Looking for a Public WPS Service). I am aware that some operations may be very time and processor consuming but still there are many simple operations that could be managed and offered by this way - but they are not. Does anybody know why? There are many libraries that allows to publish WPS... Can be the problem in difficulty of using WPS for end users? Or people just don't want to use such operations online? Anyway it is a pity because it would offer a creation of various new services over existing applications (like Google Maps) without complicated programming and with possibility of services orchestration... So my question is: **Why it is still not common to publish and use WPS?**


title: Web Map Server (geoserver) on CentOS 6

text: I am trying to setup a web map server on a CentOS 6.2 virtual machine for development. I have never used Linux or web mapping frameworks, but I am trying to do this work as low cost as possible. So far, I am having some difficulty determining how to install and configure all the packages I need. Here are the packages I want to install: Apache Web Server Apache Tomcat? - is this necessary? Geoserver PostgreSQL/PostGIS GeoExt? Mapfish? Has anyone had any experience with these packages in a Linux server? Could you provide any advice or links for setting this up? Thanks in advance


title: What to use to divide longitude on smaller parts than time zones?

text: I have a database of earth's grid, I have consolidated elements from 1152x576 to 360°x180° and then those degrees into geographical zones and time zones so I get 5(6) geographical zones (N/S frigid, N/S temperate, (N/S) torrid) and ~24 time zones ((+/-)0, +1 to +12 (or +-12) and -12 to -1) and that's the problem, because it's 5:24 zones (latitude:longitude) and it's bad ratio. I am displaying those zones in a table and it doesn't look well, when the x axis is 5 times longer than y axis. Is there widely used consolidation of longitude into less, but bigger areas? I was thinking for example about 8 areas, 4 on western and 4 on eastern hemisphere, but I could't find if this division of longitude exists. Would it be somthing like WWW, WWE, WEW, WEE, EWW, EWE, EEW, EEE? Is it used normally so the people looking into my database will understand it? I can make it look like "WWW (0°-45°)" etc, I think it will be the best, am I right? EDIT: I encountered a new problem, that this division of longitude has different borders than division on time zones, time zone start on 7.5° + x*15° and this division would start on 0° + x*45°, so I couln't consolidate 3 time zones into one portion of 45° and I consolidated those 2 approaches separately.


title: Looking for a Public WPS Service

text: Can anyone point me to a public WPS (Web Processing Service)? I need to do some testing on it and results in Google and Bing show nothing. Thanks.


title: How to identify image location?

text: ![enter image description here](http://i.stack.imgur.com/BFt28.jpg) I have a satellite image of only one part of some city with a lot of vegetation. I need to find out where is that image taken. I'm interested is there any way to find out where it is taken on Google Maps, Google Earth or something similar? It is .jpg image so I don't know how to start search. Thank you.


title: Selection of lines that lays mostly in defined zone (majority overlap)

text: I have lines of surface runoff and I need to select the ones that are lying mostly in internal buffer of polygon. I'm trying to make some script in python, but it's very important to make the script simpliest as it's possible, because there are 30 milions of this lines (and some about 650 000 polygons) and this is only a part of whole problem. Has somebody any suggestion? I'm using ArcGIS 10.2. ![Sample of my problem - red circles shows lines that are need to be selected](http://i.stack.imgur.com/wXxlE.png,) ![There are the raw data](http://i.stack.imgur.com/cyTUZ.png) There are the raw data.


title: Ordnance Survey 1 km Grid

text: I am wanting to get a copy of the OS 1km Grid layer. I have found these previous questions:     * Ordnance Survey grids - download or generate? and   * Ordnance Survey Grids But neither provide the information I required. The first question does provide some links to the OS web site but these now return the 404 error. A subsequent comment mentions its half way down this page:   * https://www.ordnancesurvey.co.uk/opendatadownload/products.html But which specific item should I be downloading? It is not obvious!


title: Batch Processing 400 MODIS HDF4 Files Using Several Tools In ArcGIS 10

text: This is my first post, I hope I don't screw it up. We have about 400 MODIS HDF4 files we are looking to batch process in the ArcGIS 10 or Python environment. The tools we are looking to use are Extract Subdataset (HDF4 files have subdatasets and we are only interested in one of them, which its index number is 1), ProjectRaster(to UTM 11 NAD 83), convert HDF4 files to TIFF (although I believe this can be done when executing the Extract Subdataset tool as it lets you choose an output file format), Int (we need the raster file to have an attribute table with pixel values). I've tried performing these tasks in modelbuilder in batch mode but i've had no success in properly setting it up. We are interested in both the ModelBuilder and Python approach. It would be great if someone can post some batch processing code that would apply to these tools.


title: ArcGIS Server cache from outside wms service

text: I have an arcgis server web application (version 10.0) that utilizes an external wms service for displaying some information layers on top of a local basemap. This external service has a tendency to be rather slow and buggy so I am wondering if it would be posible to use a remote wms service as a source for creating a local cache from within arcgis server, following its tile scheme. The method would be something like: a) connect to thewms service in ArcGIS Desktop b) save the mapfile and publish it to our local server c) set the appropriate scales via ArcCatalog d) creating the cache from ArcToolbox Has anybody tried this approach and are willing to share their experience?


title: QGIS and gvSIG comparison

text: Many debutant GIS users are asking me which one is better QGIS or gvSIG. Since I use Qgis and I don't have any experience in gvSIG, I aways say that both are very nice softwares and the chose is a matter of liking vanilla or chocolate.   But I want to have a better answer for that question, so I would like to hear from people with experience on both of them, what drove you to choose one or the other as your favorite opensource desktop GIS.


title: Is "Hawth's Tools" available for ArcGIS 10?

text: I am looking for Hawth's Tools for ArcGIS 10. Is it avaialble from here? I found that similar functions are in the new program GME, but Hawth's Tools for ArcGIS 10 is not visible yet. Will anyone help me? thanking you in anticipation