title: DescribeProcess response in Geoserver

text: The Encoding and the Schema tag for ComplexData in a DescribeProcess-Response document is an optional parameter (defined in the WPS Specification). In GeoServer, the Response document doesn't contain these parameters. Also, for a literal input parameter, the parameter "DataType" isn't shown. What do I have to change in GeoServer, that these parameters are shown in the response document? (Respectively, how can I edit the parameters?) Thanks for helping.


title: How to best alter surface from LIDAR for designing low water crossing and its approaches?

text: I was wondering if you might be willing to help me work out the best way to use LiDAR data for my thesis project. It involves altering the terrain to add low water crossings at various angles and seeing what effect they have on the floodplains and sediment transport (-the analysis portion should be able to be done using Hec-ras, using TINs or DEM input files). I seem to keep getting stuck on which route is best to alter the data. I've been altering it via 1m contours and creating new TINs, however it seems to make a very choppy TIN and I am afraid it is altering more terrain than necessary to add this 5m wide crossing. I would appreciate any and all insight you can provide me, I'm at a wits end here.


title: Mosaic to New Raster() tool seems to interpolate data gaps on one computer but not another, what's going on?

text: i have a group of raster files in .tif format that I want to mosaic together into one .tif. I've noticed that I get data gaps along the seams when I do Mosaic to New Raster() on my home computer (ArcMap 10.1). Also, these data gaps are not necessarily where there is no data because when I display the original .tif files, there is coverage where the final mosaic-ed raster shows "no data" values. There should be some gaps because the original rasters do not completely overlap, but there is more coverage than what is shown here as gaps. Here is a screen shot of the mosaic-ed raster with the data gaps:![enter image description here](http://i.stack.imgur.com/rTINt.png) Then, when I do the same operation on the same rasters in the same order, but on a different computer (ArcMap 10), there are no gaps at all. Here is a screen shot of the mosiac-ed raster in the same region as above: ![enter image description here](http://i.stack.imgur.com/NTWYM.png) Ideally, I would like a final single raster with no data gaps so I would like to use the no gap raster, but is it reliable? I've read that Mosaic to New Raster does not interpolate the no data areas, so what could be going on? How do I know which final layer is reliable?


title: How to use ArcGisServer from Engine?

text: MapServerQuerySublayer (Version 10) does not seem to support setting a definition query, although according to documentation it should using ImapServerSublayer2.LayerDescription as IMapserverSublayer. Nor does it support selections, i.e. rendering objects matching a query or selection differently. What are our options when making an Engine application which contains a MapControl consuming only arcgisserver layers? In the application we would like to change an individual layer's definition query to only show object that matches the user's interest. We would also like to render some objects as "selected" to reflect user selections. It would also be nice to be able to get the selected objects to the client as features, which also does not seem to be supported out of the box in ArcGisServer. Any ideas on how to accomplish this? **Edit:** I don't think that we can wait for ESRI to fix this bug, so I'm looking for an alternate solution, which may include fetching the selection as features and use a local layer. Engine seems to lack support for fetching features from a map service without having to connect to the context itself (not using the Layer's connection) or fetching them one by one using identify. Using the Feature service seems like a good idea, but Engine has no built in functionality to from a Layer fetch IFeatures from the FeatureService. I think that there would be a more graceful way than having to involve extra connections to the server, which implies more configuration in the application.


title: How to access map service from GeoServer and can we use that service in ArcGIS API?

text: I have installed GeoServer and publish one shapefile on it. Now my question is, How can i access that shapefile from GeoServer? and can i work with that resource (shapefile) with ArcGIS API. As we can access any resource from ArcGIS sample server like the code below: queryTask = new esri.arcgis.gmaps.QueryTask("sampleserver1.arcgisonline.com/ArcGIS/rest/services/…;); Similarly, which piece of code i need to access my published file from GeoServer ?


title: find the centroid of a cluster of points

text: While searching the web, solutions for finding centroids of polygons come up rather often. What I'm interested in is finding a centroid of a cluster of points. A weighted mean of sorts. I would appreciate it if someone could provide some pointers, pseudo code (or even better, an R package that has already solved this) or links of how this issue can be tackled. _EDIT_ Convergence has been afoot (again). iant has suggested a method to average coordinates and use that for the centroid. This is exactly what crossed my mind when I saw the right picture on this web page. Here is some simple R code to draw the following figure that demonstrates this (× is the centroid):               xcor <- rchisq(10, 3, 2)     ycor <- runif(10, min = 1, max = 100)     mx <- mean(xcor)     my <- mean(ycor)          plot(xcor, ycor, pch = 1)     points(mx, my, pch = 3)      ![enter image description here](http://i.stack.imgur.com/1ErHM.png) _EDIT 2_ `cluster::pam()$medoids` returns a medoid of a set of cluster. This is an example shamelessly stolen from @Joris Meys:               library(cluster)     df <- data.frame(X = rnorm(100, 0), Y = rpois(100, 2))     plot(df$X, df$Y)     points(pam(df, 1)$medoids, pch = 16, col = "red")


title: Open Source Map Composing Tool

text: Are there any open source tools which can be used for map composing, preferably GUI based. I use Quantum GIS sparingly but I am interested to know if there are others. The tool must be able to read raster/vector layers in standard formats and allow exporting in various formats. Something like GeoPDF export would be a great plus.


title: OpenLayers on mobile device: possible to use local map tiles?

text: Is that possible for an OpenLayers app (running on mobile device) to load the map tiles that are deployed in the device, if the app is able to access these files (e.g., using PhoneGap)? from this question: Using local mbtiles in openlayers? with code:               map.addLayer(new OpenLayers.Layer.TMS('overlayWorld',     'http://localhost:8888/2.0.0/overlay_62ee43/{z}/{x}/{y}.png', {     maxResolution: 156543.0339,     type: 'png',     layername: 'overlayWorld',     isBaseLayer: false     }));      it seems work. but it employs local web server. In the case of mobile device, does it still need a web-server? any advice will be appreciated!


title: getExtent() in OpenLayers returns strange values

text: I'm trying to get the bounding box of the map area displayed by OpenLayers. I tried to use map.getExtent() because it should be exactly what I need, but the returned values are "strange": For example, if I center the map around Grenoble (France, 45.19 5.71) with zoom level 14, the getExtent returns : 634409.225396,5647894.800465,643008.391077,5652672.114732 I guess I'm doing something wrong, but I can't figure out what. The getExtent() is supposed to return a bound object with Lat/Lon values... The map is created with :               var map = new OpenLayers.Map('map', {       projection: new OpenLayers.Projection("EPSG:900913"),       displayProjection: new OpenLayers.Projection("EPSG:4326")     });      I'm using OSM as my base layer.


title: Solar radiation calculation from GTOP030 DEM image

text: Can you help me to calculate solar radiation from GTOPO30 DEM image at the 30 m spatial resolution. I'm trying to do this many times using spatial analysis- solar radiation calculation tool in arcgis. But there always appeared error message saying unknown problem. what does it mean? Is there anyone who has a experince with this. your suggestion will help me a lot.


title: About FAPAR calculation

text: I want to calculate FAPAR in different land cover types using max and min value of NDVI. I have a raster image with 9 Land cover types and max, min value of NDVI in each land cover type. But I don't know the way how to calculate it in different LC types from NDVI image. what kind of function can be used to do this. what software is more useful for this kind of calculation?


title: In QGIS, has FGDB support improved this year, 2013?

text: I'm aware that there are users who successfully use Esri File Geodatabases with QGIS 1.8 Lisboa. FGDB support in QGIS requires the Esri API for FGDB, which is not open source (as I understand), and QGIS does not support it by default. Has this year (2013) brought any new significant improvements in support for the FGDB format, for improved read/write functionality for instance, or is something forthcoming for next year (2014) to lessen the setup time for using this format?


title: How stable is WPF in Arcmap?

text: I'm tempted to propose porting some windows forms based arcmap tools to WPF. Before doing that though, I'd like to hear from others who have done something similar, and what sorts of issues were encountered. **Compared to Windows Forms, how stable is WPF in Arcmap?**


title: OpenLayers: cache WFS Request for Vector Layer

text: In our mobile application, we'd like to use the map offline. Actually, we don't care about tiles, but we want the user to be able to see some Vector layers. Here is a piece of code:                           layer = new OpenLayers.Layer.Vector("Block", {                     strategies: [new OpenLayers.Strategy.BBOX()],                     filter: VMMS.GIS.GetBlockFilter(i),                     protocol: new OpenLayers.Protocol.WFS({                         url: self.urlProxyGeoserver,                         version: "1.1.0",                         featureType: "block",                         featureNS: self.urlGeoserver + self.dataStore,                         featurePrefix: self.dataStore,                         srsName: "EPSG:900913",                         geometryName: "block_geom",                         schema: self.urlGeoserver + "wfs/DescribeFeatureType?version=1.1.0&typename=" + self.dataStore + ":block"                          }),                     projection: "EPSG:900913",                     styleMap : self.StyleBlockWithOpacity                 });      I would like to be able to load this layer before to be offline, then use it. I tried to do it, and cache it with localstorage (HTML5), but some problems occurs because of javascript cyclic object. I could make the request (with url parameter) with my own API, and then cache the response. But I have the feeling that we have to use OpenLayers API. Otherwise, do you have any idea? On google, I have only seen solution to cache tiles, but not about my problem. I'm kind of newbie to Openlayers, so I might have miss something. In that case, I would be happy to know what :) Thanks Edit: By the way, I tried the recommandation there: WFS Caching / Back button But it's not working (or doesn't match with my needs for some)


title: Using Python to execute Spatial Ecology GME commmands

text: I think I just need a second set of eyes on this. I'm just using Python to capture some user inputs to pass down to GME command. I'm only using the export to CSV to hammer out the behavior, but the final outcome will be an ArcPY tool for creating sample points. Fun stuff. (Please ignore my imported modules).                   import arcpy, os, sys, subprocess          inFeatures = r"D:\TEMP\junk.shp"     OutputLocation = r"D:\TEMP"     OutputBasename = os.path.splitext(inFeatures)[0]     OutputCSV = OutputBasename + ".csv"     GME = r'C:\Program Files (x86)\SpatialEcology\GME\SEGME.exe'     ARG = ' export.csv(in=\"' + inFeatures + '", out=\"' + OutputCSV + '");'     CMD = GME + ARG     print CMD     subprocess.call(CMD)      The printed command is as follows: > C:\Program Files (x86)\SpatialEcology\GME\SEGME.exe > export.csv(in="D:\TEMP\junk.shp", out="D:\TEMP\junk.csv"); And that looks to be correct I have excluded the '-c' switch for the sake of debugging. When GME is launched, the following output is provided. > > Found R here: C:\Program Files\R\R-3.0.1\bin To switch to a different > version of R use the 'Set R executable folder' item in the File menu, or use > the r.setpath() command. export.csv(in=D:\TEMP\junk.shp, > out=D:\TEMP\junk.csv); Error: Variable name could not be resolved > (D:\TEMP\junk.shp) or quotes were omitted when providing a text parameter; > command cannot be interpreted The actual command within GME should look like this:               export.csv(in="D:\TEMP\junk.shp", out="D:\TEMP\junk.csv");      So, for some reason my quotes are not getting passed to GME? The only clue I have related to this in the documentation is as follows: > Please note that you need to add a backslash before all quote marks in the > commands text, otherwise it is the same syntax as GME. For instance, here is > a Python call to GME that tells it to run the GME script called > "E:ndatangmescript.txt": subp.call(r'C:nProgram > FilesnSpatialEcologynGMEnSEGME.exe -c run(in=n"E:ndatangmescript.txtn");'); Also, please note I'm still learning the ins & outs of subprocess. Aside from 'is my syntax correct?', I'm not sure how to roll this into a more direct question. Thanks! Z EDIT: I'm adding the only Python reference in the help document: > You must rst load the subprocess Python library (using: import subprocess as > subp) at the beginning of your Python session. To call GME from the command > line you use the following generic syntax: os.system(r'pathnSEGME.exe > commands') where path is the full directory path to the SEGME command, > usually: subp.call(r'C:nProgram FilesnSpatialEcologynGMEnSEGME.exe > commands') and commands is the list of commands you wish to execute. Note > that the r at the beginning of that path name is not a mistake: it tells > Python to interpret the following text literally. Please note that you need > to add a backslash before all quote marks in the commands text, otherwise it > is the same syntax as GME. For instance, here is a Python call to GME that > 17 tells it to run the GME script called "E:\data\gmescript.txt": > subp.call(r'C:\Program Files\SpatialEcology\GME\SEGME.exe -c > run(in=\"E:\datangmescript.txt\");'); (This command has been displayed on > two lines here because it is too long for a single line, though in Python it > would be a single line). Note that by including -c directly after SEGME.exe, > this will force the GME interface to close when it has nished running the > specied commands.


title: A software architecture for viewing/editing roads and routes

text: I have a routing software module which uses PostGIS/pgRouting to calculate routes. The road network is generated from an .osm file with osm2pgrouting tool. I need to enhance the functionality of the software to meet the following requirements:   1. The user should be able to edit the roads network in a graphical editor (add/edit/delete roundabouts, lanes, road signs, highway tags, etc.);   2. The edits should be immediately updated in the geometry columns - basically, if there is a new road in the OSM map, we have a new line in the respective PostGIS column. pgRouting will use the updated roads network);   3. The user should be able to view the underlying OSM area and the edited roads network as two layers in the graphical editor. Note: I tend to think, that JOSM is a very convenient tool to edit roads, but it can not save_to/query a PostGIS-enabled database. On the other side, there is QGIS, OpenJUMP and uDig, but the support for OSM editing is far less mature than in JOSM. I thought about having a bash script which runs osm2pgrouting each time the .osm file is updated. This way I can have the .osm file edited via JOSM and the roads which are inside this .osm stored in a PostGIS-enabled database and viewed with QGIS (the most recent roads network would be accessible by QGIS). But:   * I would have to force the user to use two graphical editors;   * and QGIS does not render OSM data in a very attractive way. How can I meet the requirements 1, 2 and 3 in a more elegant way?


title: Running ARCgis 9.3.1 and 10.0 on the same machine

text: I have the full installation Media for 10 (9.3/SQL05 already loaded). And two editor licenses (one loaded and a new license file). What manner of wizardry is involved in getting around the conflict detection where I can run 9.3/SQL05 and 10.0/SQL08 on the same machine? I have a gdb that I sync with at a far remote local that can't be helped at this time. But I need some of the new functionality in 10 for other projects. their must be a way to get both working on the same machine huh? I find no real help at all on the net. Most people just upgrade to 10. I found this horrible post: http://blogs.esri.com/Dev/blogs/arcgis/archive/2009/10/23/Installing- ArcGIS-9.4-_2800_and-can-you-leave-9.3.1-on-the-machine_2900_.aspx But their must be a way.


title: How do I add a layer and a feature with one click

text: In Openlayers, I want on a single click on the map for a layer to be added and a feature to be dropped where I click.               vectorLayer = new OpenLayers.Layer.Vector("layer", {styleMap: currentStyle,'displayInLayerSwitcher': true});     map.addLayer(vectorLayer);     drawFeature = new OpenLayers.Control.DrawFeature(vectorLayer, OpenLayers.Handler.Point);     map.addControl(drawFeature);      I think this should probably be done with events somehow. But every option I try invovles a double click - one to add the layer, and one to add the feature. Anybody know how to solve this crisis? Thanks!


title: Dividing a polygon into two OpenLayers

text: I'm using OpenLayers to draw a polygon to my users. I want the ability to implement through a line, split the polygon into two, with the cutoff line that the user drew. Does anyone know if this is possible to do with OpenLayers? Maybe indicate some link or sample stuff?


title: Unencoding JSON output from OSRM

text: I set up Open Source Routing Machine (OSRM). The documentation does a pretty good job of laying everything out and I was able to get a server set up and running on my machine. My question is about handling the JSON output from my server. As explained in the documentation the routing information is encoded. My application was built accepting JSON from Mapquest's Open API which is unencoded. My question is do I need to perform a special function to unencode the data? Or better yet, is their some feature baked into OSRM I am not using properly to get the unencoded results? _**edit_** **Solved** : Adding `?compression=false` to the request query will output unencoded data.


title: Automatically discard nodata area of raster

text: Given a raster with an overly large nodata surround, how can I automatically discard the surround? Said another way: how does one shrink the raster extent to only occupy the area containing values? I know one can interactively draw a box and clip the raster, or zoom and export only visible extent, or set an analysis mask. I'm after a geoprocessing tool or simple python script. Tools at hand are ArcGIS with Spatial Analyst, QGIS, and GDAl/OGR, but really anything is fair game. The example raster is 16bit signed integer, with `-32768` defined as `NoData`, but the null value or bit level really shouldn't matter. ![nodata extent is many times larger than the actual data extent](http://i.stack.imgur.com/jtcO6.png)


title: Which GIS software or plugins can perform a Cellular Automata Simulation?

text: I am interested in running cellular Automata to run a simulation on city growth. Which GIS software can run this out of the box, or via plugins?


title: Dynamic SLD styling - unique colours for linestrings

text: I have a postGIS table that looks like this: (plus a 'track' geometry column, a linestring for the vessel's track)               sampling_interval | voyage_name |   start_time_utc    | id |      ship_name       | start_port | end_port  |           description           |   last_update_utc     ------------------+-------------+---------------------+----+----------------------+------------+-----------+---------------------------------+---------------------                     5 | ss2011_e01  | 2011-02-25 23:15:35 | 73 | RV Southern Surveyor | Hobart     | Sydney    | Transit to drydock              | 2011-03-02 02:26:55                     5 | ss2011_v05  | 2011-09-21 06:36:20 | 74 | RV Southern Surveyor | Geraldton  | Geraldton | GA Hydrocarbons                 | 2011-10-17 23:02:25                     5 | ss2012_v01  | 2012-04-20 00:02:45 | 82 | RV Southern Surveyor | Brisbane   | Brisbane  | Sustained monitoring of the EAC | 2012-04-29 05:48:45                     5 | ss2012_t02  | 2012-05-02 03:01:35 | 86 | RV Southern Surveyor | Brisbane   | Lautoka   | Transit Brisbane Lautoka        | 2012-05-09 22:14:00      I have a WMS layer in geoserver pointing to this. What I would like to do is ensure that each vessel track is a distinguishable colour. New voyages are added to this table regularly; if at all possible I would rather not have to edit the SLD next time (and if I do have to edit it, it'd need to be scriptable). Using a SLD with some logic in it to do this styling seems like the obvious approach, but from looking at the guides/docs, I'm not sure if this is possible (this is really my first exposure to SLDs). Is this possible/sensible using SLDs? I'm open to other approaches, of course, if there is something that makes more sense. Any pointers appreciated - like I said, I'm very much a newbie with SLDs (but comfortable with postGIS/Geoserver/Java etc)


title: Start session edit on a specific layer on the map?

text: I created a GIS using ArcObjects, my application needs to be connected with a geodatabase stored in SQL Server. I used ARCSDE to connect my application to the database and when the database is connected, my application loads two tables - the first one is called `Commune` and the second one is called `Parcelle`. Their `esriDatasetType` is: `esriDatasetType.esriDTFeatureClass`. A Parcelle belongs to only one Commune; and a Commmune can contain many Parcelles. What I want to do is start the session and edit the Parcelle layer without the possibility to edit the Commune layer. When I use my code, I can edit both layers in the map:               public static void Start(IEngineEditor _engineEditor, IMapControl4 map, string targetLayerName)     {         //If an edit session has already been started, exit.         if (_engineEditor.EditState == esriEngineEditState.esriEngineStateEditing) return;              //Start editing the workspace of the first feature layer found.         for (int layerCounter = 0; layerCounter <= map.LayerCount - 1; layerCounter++)         {             ILayer currentLayer = map.get_Layer(layerCounter);             if (currentLayer.Name != targetLayerName) continue;             try             {                 if (currentLayer is IFeatureLayer)                 {                     IFeatureLayer featureLayer = currentLayer as IFeatureLayer;                     IDataset dataset = featureLayer.FeatureClass as IDataset;                     IWorkspace workspace = dataset.Workspace;                     _engineEditor.StartEditing(workspace, map.ActiveView.FocusMap);                          ((IEngineEditLayers)_engineEditor).SetTargetLayer(featureLayer, 0);                     break;                 }             }             catch (Exception ex)             {                 Dialog.info("Info", ex.Message);             }         }     }      My function has 3 parameters: `IEngineEditor`, `IMapControl4` and the name of the layer that I want to edit, which is the layer Parcelle in my case.


title: What's the best way to display a huge Aerial image using Geoserver

text: I have an Aerial image - Almost 1.5 GB , GEOTIFF - and i have been trying to add it as a layer in Geoserver .. Adding the layer went successful , but the performance were very bad . I tried to cache the layer using the Geowebcache , and still the performance is bad . Any tips ?


title: How to import a postgis database in ubuntu (osgeo live 7)

text: How does one import a postgis database in the ubuntu/osgeo live environment? In windows I would copy the backed up database file into the postgres bin folder and then use cmd to restore the databse. As I am new to the ubuntu/osge live environment I am not sure how to go about restoring a postgis database.


title: Is it possible to disable an Extension Addin on Startup

text: I am currently creating an Addin, that can be toggled on/off from the extension manager. What I want to do is when ArcMap is being loaded, to check the license level being used (which I have done), what I then would like to happen is that the extension in the Manager is disabled and the text changed to something like "xxx ArcEditor Required". Is this possible to do? I know I can change the name of the Extension via the config file, but cant seem to find the property to do it outside of it. thanks!


title: How to config geoserver to read style of vector layer from postgis/SQL store?

text: I've imported vector layers from map info pro to both SQL server and post gis using easy loader. In both cases I've included style column in the exported table. When I create new store in geoserver and see the list of layers , each layer has the mi_style column in the list of features at the bottom of the screen . But when I publish the layer it still have the default style of geoserver . How can I set it to read the style from the column ?


title: How to integrate XML markers with MarkerCluster?

text:
I'm having a problem integrating XML markers with MarkerClusterer. I need to integrate these two and show only the markers within bounds. I'm able to load the results in the XML to the map, but can not integrate the MarkerClusterer.               var customIcons = {           Commercial: {             icon: 'http://cdn1.iconfinder.com/data/icons/free-business-desktop-icons/32/Company.png',           },           Condo: {             icon: 'http://cdn1.iconfinder.com/data/icons/REALVISTA/project_managment/png/32/agency.png',           },           Home: {             icon: 'http://cdn1.iconfinder.com/data/icons/free-business-desktop-icons/32/Home.png',           },           Lot: {             icon: 'http://cdn1.iconfinder.com/data/icons/Real-Estate-png/32/Lands.png',           }         };               function load() {       var map = new google.maps.Map(document.getElementById("map"), {         center: new google.maps.LatLng(<?php echo $lat ?>, <?php echo $long ?>),         zoom: <?php echo $zoom ?>,         mapTypeId: 'roadmap'       });       var infoWindow = new google.maps.InfoWindow;            downloadUrl("mapxml.php", function(data) {         var xml = data.responseXML;         var markers = xml.documentElement.getElementsByTagName("marker");         for (var i = 0; i < markers.length; i++) {           var name = markers[i].getAttribute("name");           var address = markers[i].getAttribute("address");           var type = markers[i].getAttribute("type");           var point = new google.maps.LatLng(               parseFloat(markers[i].getAttribute("lat")),               parseFloat(markers[i].getAttribute("lng")));           var html = "" + name + "
" + address;           var icon = customIcons[type] || {};           var marker = new google.maps.Marker({             map: map,             position: point,             icon: icon.icon,           });           bindInfoWindow(marker, map, infoWindow, html);         }       });     }


title: How to fix PJ_healpix.c in PROJ.4?

text: Michael Speth and i are in the process of fixing the PROJ.4 implementation of the HEALPix map projection (see http://trac.osgeo.org/proj/browser/trunk/proj/src/PJ_healpix.c for the current buggy implementation) and are running into some difficulties in the ellipsoidal case. Here's the issue. We've implemented the spherical version of the HEALPix forward projection via               FORWARD(s_healpix_forward); /* sphere */         return healpix_sphere(lp, P);     }      and that works fine. We now want to reuse healpix_sphere() in the ellipsoidal version of the HEALPix forward projection by simply replacing latitude `lp.phi` with authalic latitude and replacing the ellipsoid `P` with its authalic sphere. What's the best way to do this in PROJ.4? We tried               FORWARD(e_healpix_forward); /* ellipsoid */         lp.phi = auth_lat(lp.phi, P->e, 0);         P->a = auth_radius(P->a, P->e);         P->e = 0;         return healpix_sphere(lp, P);     }      where auth_lat() and auth_radius() are our homemade authalic latitude and radius functions. But the numbers coming out aren't correct. We think that our attempt at changing P into its authalic sphere via the P->a and P->e assignments is failing. Any tips? Michael also asked this question on the PROJ.4 forum (http://lists.maptools.org/pipermail/proj/2012-November/006465.html) last week but hasn't received a response. Thanks for your attention.


title: Export Web Map Task for 3076 features

text: I am trying to use Export Web Map Task but it has a limit of 1000 records, how can I use this service or replicate this service to execute my print command. Any suggestions? I am working on ArcServer 10 not 10.1 and so cannot use the print task and print widget.


title: pgRouting - Make fails on Mac

text: I am having trouble getting started with pgRouting. 1) from binary The Kyngchaos binary just won't install and says it requires PostgreSQL 9.1. I do have it, but probably not where the binary expects it:               herve:pgrouting-master hsenot$ psql --version     psql (PostgreSQL) 9.1.4     herve:pgrouting-master hsenot$ which psql     /opt/opengeo/pgsql/9.1/bin/psql      Any advice on making it available to the installer, by tweaking the path maybe?               herve:~ hsenot$ echo $PATH     /opt/boost_1_52_0:/Library/Frameworks/GDAL.framework/Programs:/Library/Frameworks/UnixImageIO.framework/Programs:/Library/Frameworks/PROJ.framework/Programs:/usr/bin:/bin:/usr/sbin:/sbin:/usr/local/bin:/usr/local/git/bin:/opt/opengeo/pgsql/9.1/bin:/opt/opengeo/opengeosuite-2.5-sdk/bin:/usr/X11/bin:/opt/local/include:/opt/local/bin:/opt/local/sbin      2) from source I've tried with the pgRouting 1.05 (tar) and the lastest GitHub version. In both case, it throws the same error. After "cmake .", I had to manually set the PostgreSQL directories in CMakeCache.txt:               //Path to a file.     POSTGRESQL_INCLUDE_DIR:PATH=/opt/opengeo/pgsql/9.1/include/postgresql/server     //Path to a library.     POSTGRESQL_LIBRARIES:FILEPATH=/opt/opengeo/pgsql/9.1/lib      Without this, cmake (at the next step) says it can't find postgresql.h The 'cmake' trace is:               herve:pgrouting-master hsenot$ cmake .     -- Boost version: 1.52.0     Boost headers were found here: /opt/boost_1_52_0     Output directory for libraries is set to sh:     /opt/opengeo/pgsql/9.1/bin/postgres/pg_config: Not a directory     Installation directory for libraries is set to sh:     /opt/opengeo/pgsql/9.1/bin/postgres/pg_config: Not a directory and for SQL files is set to /usr/share/pgrouting     -- Configuring done     -- Generating done     -- Build files have been written to: /tmp/pgrouting-master      So there is already something fishy here with the "Not a directory" messages. The 'make' error trace is:               herve:pgrouting-master hsenot$ make     Scanning dependencies of target routing     [ 16%] Building C object core/src/CMakeFiles/routing.dir/dijkstra.o     [ 33%] Building C object core/src/CMakeFiles/routing.dir/astar.o     [ 50%] Building C object core/src/CMakeFiles/routing.dir/shooting_star.o     [ 66%] Building CXX object core/src/CMakeFiles/routing.dir/boost_wrapper.o     /opt/boost_1_52_0/boost/graph/detail/adjacency_list.hpp: In instantiation of ‘boost::detail::adj_list_any_edge_pmap::bind_<boost::adjacency_list<boost::listS, boost::vecS, boost::directedS, boost::no_property, Vertex, boost::no_property, boost::listS>, Vertex, boost::edge_weight_t>’:     /opt/boost_1_52_0/boost/graph/detail/adjacency_list.hpp:2687:   instantiated from ‘boost::detail::adj_list_choose_edge_pmap<boost::edge_weight_t, boost::adjacency_list<boost::listS, boost::vecS, boost::directedS, boost::no_property, Vertex, boost::no_property, boost::listS>, Vertex>’     /opt/boost_1_52_0/boost/graph/detail/adjacency_list.hpp:2690:   instantiated from ‘boost::detail::adj_list_edge_property_selector::bind_<boost::adjacency_list<boost::listS, boost::vecS, boost::directedS, boost::no_property, Vertex, boost::no_property, boost::listS>, Vertex, boost::edge_weight_t>’     /opt/boost_1_52_0/boost/graph/properties.hpp:215:   instantiated from ‘boost::detail::edge_property_map<boost::adjacency_list<boost::listS, boost::vecS, boost::directedS, boost::no_property, Vertex, boost::no_property, boost::listS>, boost::edge_weight_t>’     /opt/boost_1_52_0/boost/graph/properties.hpp:233:   instantiated from ‘boost::property_map<boost::adjacency_list<boost::listS, boost::vecS, boost::directedS, boost::no_property, Vertex, boost::no_property, boost::listS>, boost::edge_weight_t>’     /tmp/pgrouting-master/core/src/boost_wrapper.cpp:77:   instantiated from here     /opt/boost_1_52_0/boost/graph/detail/adjacency_list.hpp:2653: error: forming reference to void     /opt/boost_1_52_0/boost/graph/detail/adjacency_list.hpp:2654: error: forming reference to void     /opt/boost_1_52_0/boost/graph/detail/adjacency_list.hpp:2658: error: forming reference to void     /opt/boost_1_52_0/boost/graph/detail/adjacency_list.hpp:2661: error: forming reference to void     /tmp/pgrouting-master/core/src/boost_wrapper.cpp: In function ‘int boost_dijkstra(edge_t*, unsigned int, int, int, bool, bool, path_element_t**, int*, char**)’:     /tmp/pgrouting-master/core/src/boost_wrapper.cpp:77: error: no matching function for call to ‘get(boost::edge_weight_t, boost_dijkstra(edge_t*, unsigned int, int, int, bool, bool, path_element_t**, int*, char**)::graph_t&)’     make[2]: *** [core/src/CMakeFiles/routing.dir/boost_wrapper.o] Error 1     make[1]: *** [core/src/CMakeFiles/routing.dir/all] Error 2     make: *** [all] Error 2      Has anyone encountered this one before? More generally, has anyone compiled pgRouting on Mac recently? Any suggestion welcome. Cheers


title: Check if lat/long point is within a set of polygons using google maps

text: I know how to do the above using Esri (Query task on shapefile), but can this also be done using Google Maps? Do I have to query each polygon or is there one method to query a set of polygons?


title: Looking for vector data for the country Turkey

text: specifically roads, railroads, and rivers shape files.... no cost would be nice, but it needs to be more detailed and accurate than mapcruzin's free data. Something under $500 would also be considered. Anyone here have experience / recommendations with locating low cost shape file data for Turkey?


title: What Algorithm is used by ArcGIS Watershed tool?

text: Does anyone know what type of algorithm is used in the ArcGIS Watershed tool (in the Spatial Analyst package)? Very little information given on Esri's website ... but I suspect it may be some kind of depth/breadth search. I've looked at these ArcGIS Online Help pages:   * **Watershed (Spatial Analyst)**   * **How Watershed works** So yes, it uses the flow direction raster, but what algorithm is it using to traverse the raster? Please note, I'm not looking for answers along the lines of 'it uses D8..'...D8 is not really an algorithm, but a model to help define the algorithm you would use. I.E you could implement the D8 scheme within a depth- first search algorithm and/or a breadth-first search algorithm


title: Populate feature grid with database query results Using JSON, JSP

text: I convert database query result into json format.But this json data not populate on feature grid, Not showing any errors.not able to understand the error.please tell me what is wrong with this code I am using Geoserver,Openlayers,PostGIS,GeoExt,Gxp This is Js page                   addOutput: function() {                    var vecLayer = new OpenLayers.Layer.Vector("vector");             map.addLayers([vecLayer]);               return gxp.plugins.fsi_query.superclass.addOutput.call(this,{              xtype: 'form',              border: false,              bodyStyle: "padding: 10px",              width:300,              items:[{xtype: 'textfield',                       fieldLabel: 'Enter From Date',                       id:'Dftext'                  },                  {xtype: 'textfield',                      fieldLabel: 'Enter To Date',                       id:'Dttext'                  },             {xtype: 'textfield',               fieldLabel: 'Enter  From depth',                id:'De1text'               },          {xtype: 'textfield',                fieldLabel: 'Enter To depth',               id:'De2text'              }],              bbar: [{                  text: "Query FSI",                  iconCls: "gxp-icon-fsiquery",                  handler: function() {                           var fdate = Ext.getCmp('Dftext').getRawValue();                      var tdate = Ext.getCmp('Dttext').getRawValue();                      var fdepth = Ext.getCmp('De1text').getRawValue();                      var tdepth = Ext.getCmp('De2text').getRawValue();                           // create feature store, binding it to the vector layer                        var store = new GeoExt.data.FeatureStore({                             layer: vecLayer,                             fields: [ {name: 'depth', type: 'string'},                                     {name: 'species_name', type: 'string'},                                     {name: 'total_catch', type: 'integer'},                                     {name: 'number', type: 'integer'}                                 ],                             proxy: new GeoExt.data.ProtocolProxy({                                 protocol: new OpenLayers.Protocol.HTTP({                                     url: "http://localhost:8070/FSI/fsi_query.jsp?fdate="+fdate+'&tdate='+tdate+'&fdepth='+fdepth+'&tdepth='+tdepth,                                          format: new OpenLayers.Format.GeoJSON()                                 })                             }),                             autoLoad: true                         });                        alert("Hai111");                         // create grid panel configured with feature store                       var gridPanel = new Ext.grid.GridPanel({                             title: "Feature Grid",                             region: "east",                             store: store,                             width: 320,                             columns: [{                                 header: "depth",                                 width: 100,                                      dataIndex: "depth"                             }, {                                 header: "species_name",                                 width: 40,                                      dataIndex: "species_name"                             }, {                                 header: "total_catch",                                 width: 40,                                      dataIndex: "total_catch"                             }, {                                 header: "number",                                 width: 40,                                      dataIndex: "number"                             }],                             sm: new GeoExt.grid.FeatureSelectionModel()                          });            }                   }]          });          }      This is JSON out Put                    ["{"total_catch":120,"geometria":"POINT(92.111.2333333333)","species_name":"Thunnus albacares","number":2,"depth":1100}","{"total_catch":120,"geometria":"POINT(92.1 11.2333333333)","species_name":"Carcarhinus limbatus","number":2,"depth":1100}"]


title: Copy a spatial table from one PostGIS database to another PostGIS database

text: I have two PostGIS databases on two different servers. The first server is PostGIS 2.0 and the second server is PostGIS 1.5. On the PostGIS 2.0 I have a table "lakes" with MULTIPOLYGON geometry column. My task is to copy my table from the PostGIS 2.0 database on the first server to a PostGIS 1.5 on the second server. I can't access the first database through remote connection, and therefore I can only export the table by exporting it to a .sql file. What command do I use to export my "lakes" table to an .sql file to ensure that the geometry column information and the spatial index is correctly recreated in the second database? I used the tool shp2pgsql before. Is there a similar tool for "PostGIS to PostGIS" export?


title: gdal_merge Sum/Agregation files into single output band

text: I need to create a single tiff files from several others (each have the same extent and resolution) and the output file need to have a single band containing the aggregated/sum of the values in the separate files. I was playing with gdal_merg, - in the the gdal_merge documentation the definition include _"In areas of overlap, the last image will be copied over earlier ones."_ wich is not what I need. I got somewhat over that by using the _-separate_ parameter that give me a single file with multiple bands, each representing one of the source files. Is there a python/command line way how I can sum these multiple bands programmatically? (similar to the Raster Calculator GUI in QGIS) Or even better do it in a single step? Thanks in advance. g


title: ArcGIS and Citrix

text: We are currently running ArcGIS 10 in our company through Citrix as well as installs on user workstation. We always seem to have niggling problems with Citrix hence local installs. For example:   * Poor display of smoothed grids in Citrix showing very blurry but then perfect display when exported to a PNG   * Mapped drives not being able to connect but once the session is closed down and re-opened and it's fine again   * ArcScene has terrible performance when it has a couple of stacked surfaces to render   * ArcGIS Extension licenses not always being released when switched off in extension dialog (user has to shut down their session) There are lots more issues that can be found on the forums of course and I have read the strategy documents too but I'm looking for more detail on how to actually deploy it without these niggling problems! But obviously as a concept it is great as it saves so much time for IT and managing the application centrally is the way forward, etc. But the users are suffering at the moment. So my question is what are peoples opinions with deploying ArcGIS through Citrix? Have people had "strategists"/Citrix experts to come in and show us how to set it up specifically for ArcGIS to get the best performance and display? Basically the users should not see a difference between ArcGIS on their machines and ArcGIS on Citrix. Otherwise it won't get used... Thanks Rob


title: How to create a drop down list?

text: I want to show a drop down list in the feature attributes dialog. **Example:** I have 3 different building zones. The hot Wohnzon A, B residential zone and residential zone C. Now I have an attribute called the Zonenart. Now I track a new zone. And instead of the attribute "Zonenart" write something purely by hand, I want to be able to select it directly. Thank you very much


title: Implementing Vehicle Routing Problem in ModelBuilder?

text: I am trying to solve the "Vehicle Routing Problem", using ModelBuilder. For example, using Exercise 7: Servicing a set of orders with a fleet of vehicles from the Network Analyst Tutorial:   1. add a "Make Vehicle Routing Problem Layer" tool    2. add a "Add Locations" tool and set "stores" as "input location"    3. add a"Add Locations" tool and set "Distribution Center" as "input location"    4. add a "solve"tool to calculate.  That's all ... and I have a error "Routes class is empty." I think I have to add that something connected with the route class, but I did not understand how to do it in ModelBuilder ![enter image description here](http://i.stack.imgur.com/eIEWd.png) Click here for higher resolution image How to make the scheme correctly? PS Sorry for my english


title: 1944 Cyprus Grid coordinate data projection problem in MapInfo

text: I have coordinate data from the 1944 British Army General Staff, Geographical Section (GSGS) Cyprus map which was cast on the Cassini projection, but was referenced to the Clarke 1880 ellipsoid where `a = 6,378,249.145 m and 1/f = 293.465`. I gather from Clifford Mugnier's 2006 article that the latitude of origin should be `(φo) = 35º 00’ 00”N` and the central meridian `(λo) = 33º 19’ 00”E`, with a `False Easting = 200 km`, and a `False Northing = 150 km`. I've added these data to the prj file for the Cassini-Soldner projection `("Cassini-Soldner Cyprus", 30, 1000, 7, 33.316667, 35.00000, 200000, 150000)`, which puts the points in roughly the right location, but their locations are distorted slightly, resulting in some sites appearing in the sea, etc. Is there any way of getting the points in the exact right location? Have I forgotten to include some information? Thanks


title: How to avoid overlapping SVG icons and labels?

text: I've dropped a point and then selected a svg icon. But now I'm right next to the other points I want to label. Now you have the svg icon, but be on the label, that is, the label may no longer be seen. Is this a default setting that you can change? Thank you very much


title: Tracking a vehicle without GPS device

text: We are planning to write an application which have to track the user driving in his truck in the cities like New York city accurate upto which lane he was driving on. In the cities with those tall buildings, its not possible to get good GPS signal. So we are planning to use hardwares with accelerometers and magnetometers sensors to make up for the lost GPS signals. We tried with some Android tablets, but the sensors are not accurate enough. Right now we are going to try some hardwares from Trimble with geo position correction and motions sensors. So any kind of suggestions on tracking a vehicle without a GPS device is welcome.


title: ogr2ogr when importing from fgdb and use in python

text: I want import data from my GDB to PostGIS. It works fine on some datasets (all records -data- is imported or nothing is imported) from a gdb but just doesn't import others. I can't find any documentation on why from sites like http://www.bostongis.com/PrinterFriendly.aspx?content_name=ogr_cheatsheet So I have ogr2ogr -f "PostgreSQL" -a_srs "EPSG:4283" PG:"dbname=tp_test user=postgres port=54321 password=XXXXXXXXXXX host=localhost" d:/scratch/arcgis/scratch.gdb It only seems to report errors but if it works then you just get the command line again. When I check the postgis database, some of the datasets from the gdb are imported and others aren't. They are all in the same EPSG. Any ideas? This is part of a larger script/tool that I have already written to read my mxd's, export and georeference the data to a single format, save the style layer, convert it to SLD and ideally the next step is to save it to a postGIS database. This can then be used in Geoserver with the correct ESRI symbology. See Is it possible to convert an ESRI style sheet to a font set?


title: How do I prevent Qgis from being detected as "not responding" when running a heavy plugin?

text: I use the folowing line to inform the user about the status:               iface.mainWindow().statusBar().showMessage("Status:" + str(i))      The plugin takes about 2 min to run on my data set but windows detects it as "not responding" and stop showing the status updates. For a new user this is not so good since it looks like the program have crashed. Is there any work around so the user is not left in the dark regarding the status of the plugin?


title: Mosaic Aerial Photography Tiles to a Seamless Raster

text: I've got about 6 different Aerial Photography datasets I am managing at the moment. All of them are made up of tiles covering a large area. Some of these datasets contain up to 2000 tiles. All the tiles are slightly overlapping so there re no gaps in the data. Tiles work fine for me, but other users are requesting Mosaic's of the datasets. I've been trying over the last few days to Mosaic these tiles into a seamless raster. I've been using the ArcGIS Mosaic to New Raster Tool. It worked great on the first (smallest) dataset, but since then ArcGIS crashes after running for about 8 hours. I've tried the tool on a more power machine and with less tiles, but ArcGIS still keeps crashing. I'm not sure what else to do here. I need to mosaic these tiles into one raster, but despite any changes to the options I keep getting the same result. I'm wondering if the Moasic to New Raster is the correct tool for this job, or should I be using something else? I'm looking for a either a ArcGIS or python programming solution to this issue. The options I am using for the Mosaic to New Raster Tool:   * Input: 500 to 2000 input tiles that are JP2   * Output: Save the new raster as a JP2   * Coordinate System: NAD83 17N to NAD83 17N   * Pixel Type: 8 Bit Unsigned to 8 Bit Unsigned   * Cellsize: 0.2 to 0.2 (varies, but the the cellsize stays the same from the input to the output)   * Band: 3 (4 in two dataets)   * Mosaic Method: I've tried pretty much every option here, the preferred is Blend.   * Colourmap: First (I haven't changed this at all yet)


title: How to divide a polygon into specific sizes using ArcGIS?

text: I have several thousand irregularly shaped polygons in a shapefile. I want to be able to split each polygon into three areas, and to specify what the size of those areas (they sum to the previous total area). It does not matter what the sub-polygon shape is, as this is for visualization purposes. How can I go about doing this? Is there a standard algorithm that exists that I can use? * * * One approach I considered was to get all the points that make up the polygon. Then, I would join two randomly together using a straight line, split the polygon, then check if the area was within a satisfactory tolerance. If it was too small, I would change the point in one direction; if it was too large, I would change to a point in the opposite direction.


title: What are the implications of invalid geometries

text: I've imported some data in a Postgis database and some of the geometries are reported invalid (ST_IsValidReason reports self-intersection or ring self- intersection). The queries I am performing don't seem affected by the invalid aspect of these geometries (i'm only using ST_Distance queries). What are the things that break when geometries are invalid? Is fixing these geometries "automatically" (buffer(geom, 0) or ST_SimplifyPreserveTopology(geom, 0.0001)) an option?


title: How to get the envelope of all the features in a feature class in Arcmap?

text: I would like to get a single envelope of all the features in a feature class. Is there any way to do that in ArcMap?


title: Create centreline from a polygon

text: I am working on a project to create a centerline from a polygon automatically. The user select a polygon and click a button, the centerline will create for him. I have tried the build-in function “Collapse Dual Lines to Centerline”, but it cannot perform well so I decide not use this function in this project. In further study, I find out some terms such as “Medial Axis”, “Skeleton”, “Voronoi polygon”, those are some concepts related to centerline formation. However, I find difficult to understand those algorithm. Does anyone have a suggestion to me for this project? Does anyone have a suggestion for how to implement those concepts in ArcMap? Is there any sample script or algorithm I can refer to as reference? Development language: C#, VB.Net ArcGIS : ArcMap 9.3.1 Best regards, Gary Sham


title: postgres command to save the line coordinates in openlayers

text: I am working on project that we create one map ant editing toolbar. Then by using that editing tool bar we can draw point, line, polygon on the map. Next we are converting the line, point, polygon objects into coordinates using geometry.toString() function. After that i want to save that coordinates in postgres. I am not able to come out from this dig. If you know the route please let me know. I will follow the steps. Regards, Kiran


title: QGIS combine Vectorlayer and Rasterlayer to one shapefile

text: I have to create some Thematic Maps, especially Choropleth Maps on Matlab. Therefore i've found a shapefile which I can use, but need to draw some lines there. I've imported my shapefile as a Vector Layer and my .jpeg file as a Raster layer and want to combine both to one shapefile. I did "Polygonize" my Raster layer and combined both Vector layers to one shapefile, but it doesnt look the way as i want it. Here are some screenshots: My Vector Layer: ![Vector Layer](http://i.imgur.com/MsC6tF9.png) My Raster layer: ![Raster Layer](http://i.imgur.com/C055QLy.png) I want this result: ![Result](http://i.imgur.com/6tJVbdC.png) Does anyone know, how I can achieve this result with this transparency? Best regards.


title: Enabling CORS in OpenLayers

text: **Is it possible to enableCross-Origin Resource Sharing (CORS) in OpenLayers? If so, how?** I have a web map in OpenLayers, and I'm using an HTML canvas to render it as a PNG image. My script works great, but due to the browser's same-origin policy, I can only convert my canvas to an image if all the map tiles I've drawn to it come from the same host (e.g. `sub.domain.com`). To boost my map's loading speed, I pull tiles from an array of subdomains: `a.domain.com`, `b.domain.com`, etc. This greatly improves performance, working around the browser's limit of four simultaneous connections per host, but prevents me from ever converting my canvas to an image, as the tiles come from multiple hosts. I'm trying to use CORS to work around this problem. I've set the appropriate `Access-Control-Allow-Origin` response header for the tileset I serve with PHP, but this doesn't have an effect. I suspect this is due to not having the `Origin` header in the request (as shown in the example on the Wiki page). It would seem that some OpenLayers configuration is required. But what? Has anyone else done this successfully? Thanks in advance!


title: How to change coordinates of entire geometries in QGis

text: There is similar question but for ArcGIS. So how to do it in QGis? I have dataset from WWF which use different shoreline than GHHSG, in short all geometries need shifted down and left slightly.


title: How can I define hierarchical levels in eCognition

text: I'd like to define classes with 2 levels. Firstly I'm gonna classify some broad classes such as forest, grassland in a level 1. Furthermore, I'd like to classify these classes into some sub-classes such as Evergreen forest, deciduous forest etc. Grassland will also be classified into 3 sub-classes as well. How can I define this hierarchical levels in eCognition.


title: How to create Map Document for DotSpatial?

text: How can I create a Map Document (including layers, symbology info etc.) for DotSpatial just like we have MXD File in ESRI ARCGIS.


title: How to split a feature using GDAL

text: I want a feature to be cut into 2 or more based on deviation of a segment i.e., the angle between 2 successive segments. At places, I may also require to merge 2 or more features.


title: Pyshp + Mapserver: How to add label at specified coordinate?

text: **EDIT** : solved, went for the solution as described below I have a set of polygons defined in a xml file. I create a shapefile with Pyshp. This works fine. But now i want to add labels to the polygons. The labels are defined in my xml like this:                 <shape id="1" description="WINTER SEASONAL ZONE II" color="#FF99D9EA">         <polygon>           <location lat="50" lon="-68.5"/>           <location lat="40" lon="-68.5"/>           ....           <location lat="54.666" lon="-132.666"/>           <location lat="54.666" lon="-105"/>         </polygon>         <polygon>           <location lat="-47.1" lon="-65.9"/>           <location lat="-34" lon="-50"/>           ....           <location lat="-43.333" lon="-74.333"/>           <location lat="-45.75" lon="-74.333"/>         </polygon>         <label lat="41" lon="-40" text="WINTER SEASONAL ZONE II\n1 Nov - 31 Mar - Winter\n1 Apr - 31 Oct - Summer" />         <label lat="41" lon="172" text="NORTH PACIFIC WINTER SEASONAL ZONE\n16 Oct - 15 Apr - Winter\n16 Apr - 15 Oct - Summer" />         <label lat="-40" lon="77" text="SOUTHERN WINTER SEASONAL ZONE\n16 Apr - 15 Oct - Winter\n16 Oct - 15 Apr - Summer" />       </shape>      The labels have their specific positions. Also, the number of polygons and labels per shape are not related in anyway. What's the best way to go? I thought of creating a second POINT layer with the positions from the label tags and draw this layer on top of the polygon layer. is there a simpeler way?


title: how can i correct atmospherically landsat 8 images?

text: I am trying to correct atmospherically an image from Landsat 8 with 6s model from grass 7 and 6.4. Because there are no parameters for Landsat 8 yet, i am trying to give parameters from Landsat 7 images. when i use i.atcorr i get this message : wavelength less than 0.25 micron: let's take s(l)=s(0.25) and the command ends. How can i correct my Landsat 8 images? Is there someone who knows when Landsat 8 images will officially be in 6s model?


title: Any clue how to read a .style (ArcGIS) file in Python, for converting into another platform-specific format?

text: The problem I am face with is that I have all my data as shapefiles but my organization works with only 1 ArcGIS license. Migrating to QGIS to prominent amount of work, I need to convert my ArcGIS .style file to the .sld format of QGIS, to ensure uniformity of data capture and maps. I am working with Python for this, Could anyone give me any suggestions?


title: How to determine on what side of a line is a point is?

text: I want this in order to generate house numbers, odd or even depending on wether a feature is on one side of a street line or not . Thanks EDIT: In the end i think the most portable way is to use this( BTW line is actually 'linestring' not segment): ST_LineCrossingDirection( linefromAtoB,streetline) <0 where A is the point , B is ST_ShortestPoint(thepoint,streetline); if <0 then linefromAtoB enters through the left of The streetline; if >0 the it enters from right side of streetline). You have to determine what street line is closest to each point first. EDIT2: I made a mistake: C is the symmetric of A in relation to B= ST_ShortestPoint(thepoint,streetline). C's coordinates are `XC=2*(XB-XA) -XA` and for y `YC=2*(YB-YA) -YA`.   ST_LineCrossingDirection( linefromAtoC,streetline) <0 where A is the point , ; if <0 then linefromAtoC enters through the left of The streetline; if >0 the it enters from right side of streetline). You have to determine what street line is closest to each point first.


title: Problem with custom projection in OpenLayers

text: I'm using OpenLayers 2 with Proj4JS 1.10. When trying to change the map's projection using               Proj4js.defs['STEREOGRAPHIC:180,90'] = "+proj=stere +lat_0=90 +lon_0=180 +a=3396190.0 +b=3396190.0 +over=F +geoc=T";     map.baseLayer.projection.proj = new Proj4js.Proj('STEREOGRAPHIC:180,90');      and then zooming in to refresh the map, I get a JavaScript error tracing back to the Graticule control: `Uncaught TypeError: Cannot read property 'x' of null`. Is there an issue with my proj4 string or could this be a bug with OpenLayers or Proj4JS? Thanks


title: Multiprocessing Errors - ArcGIS implementation

text: I was wondering if anyone else in the community here has attempted to use multi-processing for spatial analyses. Namely I am trying to iterate through a series of rasters, create a multiprocessing job for each and run them through a number of geoprocessing steps within one def function. Something along the lines of               def net(RasterImage, OutFolderDir):         arcpy.env.overwriteOutput = True           arcpy.env.workspace = OutFolderDir          DEM_Prj = DEM_Prj.tif              try:             arcpy.ProjectRaster_management(RasterImage, DEM_Prj....             FocalStatistics(DEM_prj....)             ...          if __name__ == '__main__':           InputFolder = r'C:\test\somepath'           Output = r'C:\test\somepath2'           arcpy.env.workspace = InputFolder           arcpy.env.scratchWorkspace = r'C:\test.gdb'                  fcs = arcpy.ListRasters('*')         pool = multiprocessing.Pool(4)            jobs = []                          for fc in fcs:             rIn = os.path.join(InputFolder,fc)             rOut = os.path.join(Output,fc[:-4])             jobs.append(pool.apply_async(net,(rIn, rOut)))          Now the multiprocessing does run, usually for the first batch! However, I keep running into several different errors when attempting several datasets(more than 4 files - i.e. 4 core multiprocessing) including:               ERROR 010302: Unable to create the output raster: C:\somepath\sr6f8~1\FocalSt_srtm1     ERROR 010067: Error in executing grid expression.     Failed to execute (FocalStatistics).      and               ERROR 999999: Error executing function.     Failed to copy raster dataset     Failed to execute (ProjectRaster)      Notice in the first error the strange folder that is created (in the OutFolderDir location) associated with the focal statistics that nearly creates an exact replica of the final output. My question is based off your experience is it impossible to create several step geoprocessing within one multiprocessing function? Or do I need to tile these steps into their individual geoprocessing steps? Kind Regards, Bjorn _**UPDATE_** Still encoutering similar errors - moving the import functions to the def function has shown that               import arcpy      from arcpy.sa import *      cannot create an output with an added syntax warning that of import * is not allowed. _**UPDATE #2_** I know this is a late reply but I thought it might benefit someone else for future reference to my workaround that allows multiprocessing to work with arcpy. The main problem I found after returning to this problem is not the competition of the arcpy modules but rather competition over the scratchWorkspace that the ArcObjects utilize to save the temporary files. Therefore consider running a counter into the multiprocessing parsing argument to make a unique scratchWorkspace for each process i.e.               Counter = 0      for fc in fcs:                       rIn = os.path.join(InputFolder,fc)                       rOut = os.path.join(Output,fc[:-4])                             jobs.append(pool.apply_async(net,(rIn, rOut,Counter)))                     Counter += 1      Then in the main function make a specific temporary directory and assign an unique scratchWorkspace to each multiprocessing task.               def main(RasterImage,OutFolderDir,Counter)               TempFolder = os.path.join(os.path.dirname(OutFolderDir),'Temp_%s'%  (Counter))               os.mkdir(TempFolder)               arcpy.scratchWorkspace = TempFolder               ...       Hope that helps and thanks to Ragi for the inital suggestion to use separate temp workspaces - still baffled by why it originally did not work. _**Additional Resources_** ESRI Multiprocessing Blog Python,Gis and Stuff Blog


title: How to programmatically Summarize an ITable?

text: Is there a way in ArcMap plugin to make Summarize-function like functionality programmatically? There is a feature layer, lets say _Building_ , And it has fields like _OID, FloorNumber, Color, Height_ etc... So I want summary like:               Floors TotalHeight Heightest      3      120         20 ( 4 floor buildings in total have 120 meters height )      4      ...         ...


title: Python Postgres Multiprocessing speed increase

text: I'm trying to speed up my PostGIS queries using multiprocessing. My current setup is using python and psycopg2 such as below. Although this has given a speed increase, it still seems like there are bottlenecks preventing rapid speed increases and I'm not sure where to go next. I've increased a lot of the postgres parameters as suggested in 'Performance Tuning Postgres', but when I run this in AWS I never seem to be anywhere near maxing out RAM or I/O which should be the limitations for DB activity apparently. Can anyone suggest other methods for speeding this up?               > import os, sys, psycopg2, multiprocessing, time      > start = time.time()      > conn = psycopg2.connect("dbname=template_postgis_20 user=postgres")      > cur = conn.cursor()      > cur.execute("""SELECT count(*) FROM sites""")      > count = cur.fetchall()     > def getOidRanges(rownums, count):     >      >     oidranges = []     >     for row in rownums:     >         minoid = int(row[0])     >     return oidranges     >      > def mp(rownums, whereclause):     >      >     for row in rownums:     >         if row[0] == whereclause:     >             gid1 = int(row[0])     >             cur.execute("""     >             UPDATE sites SET postcode=(SELECT field62 FROM (SELECT field62, COUNT(field62) FROM addressbaseplusbh1_2     >             WHERE ST_Within(addressbaseplusbh1_2.geom, (select geom from sites where gid={0})) GROUP BY field62 ORDER BY count DESC)     >             as postcode LIMIT 1) WHERE gid = {0};""".format(gid1))     >             conn.commit()     >      >     return     >      > if __name__ == "__main__":     >      >     nbrquery = ("""SELECT gid FROM sites ORDER BY gid;""")     >     cur.execute(nbrquery)     >     rownums=cur.fetchall()     >      >     cores = (multiprocessing.cpu_count()-1)     >     procfeaturelimit = 1     >      >     oidranges = getOidRanges(rownums, count)     >      >     if len(oidranges) > 0:     >         pool = multiprocessing.Pool(cores)     >      >         for oidrange in oidranges:     >      >             whereclause = oidrange[0]     >             jobs = pool.apply_async(mp, (rownums, whereclause))     >      >         pool.close()     >         pool.join()     >         jobs.get()     >      >         try:     >             conn.commit()     >             cur.close()     >             conn.close()     >             end = time.time()     >             print end - start     >         except:     >             pass      EDIT: @Craig, would it work then just to have this as the executed block?               curs.execute("""UPDATE sites SET postcode=(SELECT field62 FROM (SELECT field62, COUNT(field62) FROM addressbaseplusbh1_2 WHERE ST_Within(addressbaseplusbh1_2.geom,      (select geom from sites where gid={0})) GROUP BY field62 ORDER BY count DESC)     as postcode LIMIT 1) WHERE gid = {0};""".format(whereclause))     return


title: How do I get the vector output segmentation algorithms to work properly in QGIS 2.0?

text: I'm trying to run the Watershed Segmentation algorithm, specifically, but I have been receiving the same error message using any of the available segmentation algorithms. I'm new to QGIS and have been running it in OSX 10.9. The error message I receive says: > Oooops! The following output layers could not be open Output vector file: > /var/folders/mq/01l7wc357319d3bxlsfp8n840000gn/T/processing/606c5f7baef84d2c9d91b5fc023a5947/-mode > vector -mode.vector.out.shp The above files could not be opened, which > probably indicates that they were not correctly produced by the executed > algorithm Checking the log information might help you see why those layers > were not created as expected After checking the log, the OTB execution console output says: > OTB execution console output Segmentation -filter watershed -in /Volumes/New > Volume/Thesis/soterra_LCmap/soterra_mosaic -filter.watershed.threshold 0.01 > -filter.watershed.level 0.1 -mode.vector.inmask /Volumes/New > Volume/hancockseg/hancockrst/hancockmask.shp -mode.vector.stitch true > -mode.vector.minsize 1 -mode.vector.simplify 0 -mode.vector.layername layer > -mode.vector.tilesize 1024 -mode.vector.startlabel 1 -mode.vector.outmode > ulco -mode vector -mode.vector.out > /var/folders/mq/01l7wc357319d3bxlsfp8n840000gn/T/processing/606c5f7baef84d2c9d91b5fc023a5947/-mode > vector -mode.vector.out.shp > > ERROR: option -mode.vector.out.shp does not exist in the application. I'm really at a loss as to how I can overcome this issue. As I said, I'm new to QGIS and don't have any experience with open source GIS/RS software in general. My experience has mostly been with ESRI, ERDAS, and IDRISI. I'm very interested in performing some basic segmentation and classification but with a very specific goal in mind which I doubt I could accomplish using the Semi- Automatic Classification Plugin. I'm tempted to give it a try. Any help, including links to resources addressing this error message would be much appreciated.


title: arcpy.SearchCursor strangely throwing IO Error

text: I've been trying to use arcpy.SearchCursor on an xls file. Something I've done countless times before. I had a script running and then yesterday it started throwing: > IOError: "H:\Test_Folder\Report.xls\'Titles Dataset$'" does not exist This only happens when I run the script from the Python window. If I were to open the Python Shell and execute it line for line:               import arcpy     for row in arcpy.SearchCursor(r"H:\Test_Folder\Report.xls\'Titles Dataset$'"):         print row      The xls is recognized and proceeds to print each row. I'm boggled by this. I have no idea why the python window cannot recognize the xls. Is there something I'm missing?


title: How to merge multiple shapefile via ArcEngine 9.3 Java?

text: I try to merge the some shapefiles into one by using ArcEngine java platform. But it seem like that parameters of my merge class is wrong ,i google it a long time and got nothing. Pls figure my problem out .The Error is :Invalid value type for parameter inputs.               public static  IGeoProcessorResult mergeTool(String shpPath,String[] shpNames,  String outFile)      {           //Initialize the geoprocessorResult object         //IGeoProcessorResult vGeoPR = null;         // Initialize the geoprocessor.         GeoProcessor vGeoP = null;         IGeoProcessorResult vGeoPR = null;         try {             vGeoP = new GeoProcessor();             IWorkspaceFactory vWSpaceFactory =new ShapefileWorkspaceFactory();             //IFeatureClass vFeaCls = vFeatureWSpace.openFeatureClass(ShapeName);            IWorkspace vWorkSpace;             IFeatureWorkspace vFeatureWSpace = null;            vWorkSpace = vWSpaceFactory.openFromFile(shpPath, 0);         vFeatureWSpace = (IFeatureWorkspace)vWorkSpace;             IFeatureClass vFeaCls = null;                   ArrayList<IFeatureClass> vInArray = new ArrayList<IFeatureClass>();             for(int i=0;i<shpNames.length;i++)             {                      vFeaCls= vFeatureWSpace.openFeatureClass(shpNames[i]);                 vInArray.add(vFeaCls);             }             //vGeoP.addMessage("I will overwrite your file if name exits!");           // vGeoP.setEnvironmentValue("workSpace", vFeatureWSpace);             vGeoP.setOverwriteOutput(true);               Merge  vMerge = new Merge(vInArray,outFile);                 vGeoPR=vGeoP.execute(vMerge,null);                 return vGeoPR;


title: Arcpy ChangePrivileges_management not working

text: I was trying to revoke access for a user using Python script, but it does not seem to work.   I mean the script runs successfully without any error, but does not what it's intended to... Here're the steps I follow:   1. I am connecting with the schema owner, and revoking access for another user, say 'B' for a particular dataset [Python]   2. The script runs OK...without any error (exits with code 0)   3. Connect with user B to chk if it's worked...and NO, I can see the data still Any idea why its hapeening?   I've tried with exclusive lock; and ArcCatalog...even then the same result :( Ujjwal


title: Bentley Microstation in comparison to ESRI ArcMap : Filetypes and Capabilities

text: I'm curious if there is a side-by-side comparison of Bentley's Microstation **v8i (Select series 3)** with ESRI's ArcMap **10.1** showing capabilities and filetype acceptance. I have heard some conversations lately where a Bentley user stated unequivocally that Microstation could do more and open more filetypes natively than ArcMap could. I'm not a fan of either over the other, but I viewed Microstation as a CAD system with some GIS capabilities whereas ArcMap was a GIS system that could do CAD work. And that compared side-by-side Microstation would be the lower ranking application as ArcMap just had more capabilities. But now I'm wondering if there is something I am unaware of or whether someone else is just quite clueless.


title: Looking for good looking maps with nice colors for presentations

text: I'm currently working on a project where we need some good looking maps with country boundaries displaying the world. We do not need the details of the maps, but the colors are important. It will we used in a Silverlight application that also should be available in offline mode, meaning no server connection. Can anybody help me finding some map resource we could use?


title: FME - How to send list of files to reader

text: I have an Excel spreadsheet which contains the file paths of MapInfo files that I want to convert to a File Geodatabase using FME. Is it possible to send this spreadsheet to the MapInfo reader so that it only converts those files? So the spreadsheet might look something like this:               File Path     C:\Temp\File1.TAB     C:\Temp\File2.TAB     D:\Folder\File3.TAB     G:\Other Folder\File4.TAB      All file paths are in column A of the spreadsheet.


title: Label in the middle of linestring in OpenLayers

text: Is there anyone can put label in the middle of the linestring dynamically? I've used labelAlign option. I've changed this value from bc, cc, cm, etc. But it doesn't work at all. There's also option for x offset and y offset. But I don't know how to calculate the right position so the label can show at the middle of linestring. Not at the first line string. Do you have any idea how to do it? Thanks,


title: How can I preserve image quality when combining aerials of different resolutions?

text: I have two different rasters that were flown at different times during the same month. One has a resolution of 4-inches, the other has a resolution of 6-inches. Using ArcGIS tools, how can I combine these rasters while preserving as much image quality as possible? I intend to use the combined output as a basemap.


title: processing MODIS hdf files with gdal

text: I'm trying to open an MODIS file (hdf format) with gdal in a python script under windows OS. I use the HDF4 driver, as read at http://www.gdal.org/frmt_hdf4.html.   But when I run it I got this error message:               ERROR 4: `MOD10A1.A2000064.h12v04.005.2007164045648.hdf' not recognised as a supported file format.       I have tried with the gdal.ALLRegister() statement with the same result. My goal is to obtain the properties of the image using GeoTransform, RasterXSize, etc. and reproject it afterward. If you could provide me with any help it would be much appreciated. ps. I know I can use the Modis Reprojection Tool (MRT) but I need it to be included in my script as it will iterate over a folder containing a lot of images. thanks


title: How can I draw features on GeoTIFF in Python?

text: I have a set of data (coords + value) and want to turn it into a heatmap in GeoTIFF. Should I calculate pixels and set them as array, or there are shortcuts? This is what I do now:                   driver = gdal.GetDriverByName('GTiff')         dst_ds = driver.Create('%s-%s.tif' % (name, zoom), height, width, 1, gdal.GDT_Float32)         raster = numpy.zeros((width, height))              cur.execute('select * from image')         for i, row in enumerate(cur):             raster[row[0], row[1]] = row[2]      I'd like to be able to do something like this:                band = dst_ds.GetRasterBand(1)      for row in cursor:          band.drawCircle((row[0], row[1]), 20, row[2])      The important part is saving the depth of data.


title: How to extract all Way elements tagged as highways from OSM data using QGIS?

text: I want to extract/filter in all Way elements within a specific region tagged as 'highway' from an OSM file. I am currently thinking of using a Python library called osmread (by A Dezhin) that I stumbled upon online. It's pretty straightforward. Here's an example usage that I got from his repository.               from osmread import parse_file, Way          highway_count = 0     for entity in parse_file('foo.osm.bz2'):     if isinstance(entity, Way) and 'highway' in entity.tags:         highway_count += 1          print "%d highways found" % highway_count      This is good enough, but I was wondering if qGIS also has a similar function that allows users to filter elements (Node, Way, Relation) given a specific set of criteria. I would prefer to use qGIS for the 'instant' visualization. In the end, I hope to build a road network and conduct some analysis. Thanks, all, in advance.


title: How to get only part of a WMS GetCapabilities request

text: Is it possible to modify the WMS GetCabilites request URL in order to get back a certain part? For instance I am only interested in this part of the XML:               <LatLonBoundingBox minx="16.6176088621753" miny="44.782391955068526" maxx="53.325886544307956" maxy="62.785875486994094"/>     <BoundingBox SRS="EPSG:32632" minx="1000000.0" miny="6000000.0" maxx="4000000.0" maxy="7000000.0"/>


title: Why does a raster with resolution 5x5m not create a grid with these dimensions?

text: I created a 5m x 5m raster using ArcGIS 10.0 Spatial Analyst from a point file. As can be seen in the screenshots of the image and the table, the resulting grid has a higher resoluation than is actually stated. Can anyone explain why this is? Here is the image of the points and raster: I expected the raster grid to be the same size as the distance between the points...ie 5m! ![Raster](http://i.stack.imgur.com/dQwjO.jpg) and here is the table ![enter image description here](http://i.stack.imgur.com/XTNbz.jpg)


title: How to create point features from 1-Many relationship?

text: I have an address point shapefile for buildings in a city. The point shapefile shows only a single point for each building; apartments are listed in a table called additional addresses. I need to create a many-to-one relate from the additional addresses table to the address point shapefile using the unique street address. I then need to make a copy of the address point for each instance in additional addresses table. I need to end up with a stack of points, one for each apartment at the address. Has anyone seen a Python script that does something similar?


title: Has anyone experienced problems editing between ArcGIS 10.1 and an older ArcSDE geodatabase 9.3.1?

text: We are unable to edit feature classes in an ArcSDE 9.3.1 Versioned GDB using an ArcGIS 10.1 .mxd map. The feature classes are part of a geometric network. Previously we were able to edit feature classes in an ArcSDE 9.3.1 Versioned GDB using an ArcGIS 10.0 .mxd map. It seems that ArcGIS Advanced (ArcInfo level) for 10.1 is incompatible with our 9.3.1 SDE.Has anyone else experienced this? Comment 7/24/2013: Thank you for the responses so far. We are unable to update the 9.3.1 GDB because 1 editor, who is deaf, utilizes special customized tools for performing her work tasks. We are getting a contractor to update these tools to work in 10.1. Then we will convert the GDB to 10.1 and this problem should dissappear. Yes, we did grant privileges to all editors. The editing problem only occurs if the feature class is part of a geometric network. Otherwise, if the feature class is in a dataset that does not have a geometric network, then editing can be performed. It appears then that ESRI's claim that 10.1 is backwards compatible is invalid for 9.3.1 with geometric networks. We will thus use only ArcGIS 9.3.1 and 10.0 until this specialized conversion is complete. Thank you for your assistance.


title: How to calculate raster class metrics by polygon zone?

text: I have a raster of morphology data where each bin has a code. I want to count the number of each code that fall under a polygon from another file. Zonal stats as table only gives me the min/max occuring codes and not count of each code. Can anyone suggest how to get count of each code for each polygon?


title: Why do our .shx files disappear?

text: We are using QGIS in our workgroup, and have recently updated to 1.8.0 Lisboa. To share our file structure we are using egnyte similar to dropbox. We all use Mac OS, and have used QGIS for years, and egnyte for several months now without any problems. After updating to 1.8.0 we have experienced a frightening problem: The shx- files from a vast amount of our shape-folders are disappearing mysteriously, leaving the data useless. Fortunately we have been able to recover most of the shx-files from back ups, but obviously we can not rely on that for the future. The shx-files that disappear seem to be from folders we have recently worked with, but I'm not able to see any other tendencies in what is gone and what is still intact. Very grateful for any help on this issue. All the best


title: Calculating distance between points using ArcGIS Desktop?

text: I want to write Python syntax to calculate the distance from point to other points in the water network and so on for others, I have 21 points. I need distances have to follow the network That means I will have 441 records. I don’t want to use route method it’s not useful.


title: How to filter the MouseWheel event in a QGIS C++ plugin?

text: I’m developing a C++ plugin for Qgis and I want to filter/override the mouse- wheel event.      I created a “maptool” to select a polygon feature on a vector layer, the feature has two attributes : an id & an integer H.   I want to override the mouse wheel event so that once I select a polygon, the wheel will allow me to change (increment/decrement) the value of attribute H and not zoom/unzoom the mapcanvas!      I created and “eventFilter” on my plugin’s base class which I install on the mapcanvas (QgsMapCanvas::installEventFilter()) once my maptool is activated, works fine to filter a key press for example but when it comes to the mouse- wheel it seems that it bypasses the eventFilter and does the usual zoom/unzoom action !! Any ideas on how to solve this problem? Thank you !


title: Extensions not installed in ArcGIS 10

text: I have the trial version of ArcGIS 10 (ArcEditor licence). The trial automatically installs and authorises all extensions for the trial period. However, in ArcMap I can only see Maplex, 3D Analyst and Spatial Analyst in my Extensions dialog box. When I check ArcGIS Administrator, it lists all the extensions as "authorised", but only the extensions I mentioned above as "installed". The others are "not installed". I am a bit confused about this, as the extensions are supposed to be automatically included as indicated here.


title: GeoExt 2 - Extjs 4

text: As a complement to my previous post : http://stackoverflow.com/questions/11309998/extjs-4-with-geoext-libraries- issue i can't find a css file for GeoExt 2 here : https://github.com/geoext/geoext2 i just could've find GeoExt javascript in src folder i really need to use GeoExt 2 so anyone can find me a good source ?


title: Wrong data with FEATURE_COUNT=1 for GetFeatureInfo

text: When I am querying the WMS service with FEATURE_COUNT=1 then I get the wrong meta data. Once I make FEATURE_COUNT=2 then first row contains the correct data. I am executing the URL directly in Browser (Internet Explorer) What is the reason for this. https://www2.landgate.wa.gov.au/ows/wmspublic?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetFeatureInfo&LAYERS=LGATE-071&QUERY_LAYERS=LGATE-071&SRS=EPSG:4283&BBOX=117.747,-24.152,118.687,-23.329&HEIGHT=400&WIDTH=400&X=118.232&Y=-23.735&FEATURE_COUNT=2 Thanks gsalunkhe


title: How to Load RDF Data to PostGIS?

text: I want to load some sample RDF (Relational Data Format) to POSTGIS. Are there any tools or methods to do this.


title: FME - How to add an attribute when schema definition is dynamic?

text: I have an FME workspace which reads a folder of MapInfo TAB files and writes them to a File Geodatabase. The reader and writers are both using dynamic schema as I need to run this workbench on multiple folders each containing several files. What I would like to do is add an attribute to the feature class which contains the file path of the input TAB file. I have exposed the **`fme_dataset`** attribute on the reader and I have added it to the writer as well. That didn't work. I've also tried copying the **`fme_dataset`** attribute to another attribute called **`File_Path`**. That didn't work either. I'm assuming that this isn't working because the schema definition is coming from the reader. But, why won't it work even when I expose the **`fme_dataset`** attribute? I've looked at **this similar post** on FMEPedia, but I don't see how I can use a schema mapper for this purpose. **EDIT 07/01/2014** I thought this might be relevant as well, so I'll add the following: I have created a separate attribute called **`FeatureClassName`** which is just the source file name which has been validated (removes non-word characters such as brackets). In the dynamic properties for the destination geodatabase, I have used a custom schema definition: ![enter image description here](http://i.stack.imgur.com/7NVxi.png) I am setting the Feature Type Name from that attribute: ![enter image description here](http://i.stack.imgur.com/bqmD4.png) The Schema definition is set to default (from Feature Type Name - above). Is it possible that this could be the reason I'm not able to add the attribute?


title: How to create labels using 2 columns?

text: I want to label a layer with 2 columns/attributes. It´s not a problem in Mapinfo, but can it be done in QGIS as well?


title: Find area with highest relative density of one type of object

text: I'm using PostgreSQL/PostGIS and I have an advanced type of query. My database contains two types of object - let's say pubs and churches :) Is there a way that I can find the area with the densest concentration of pubs compared to churches, and vice versa? I can obviously display this visually with a heatmap, but I just wondered if there was a programmatic way to do this with PostgreSQL? Perhaps the point that is closest to the most churches relative to pubs? Thanks.


title: Geoprocessing service - empty memory after calling

text: I'm following this sample for ADF libraries http://edndoc.esri.com/arcobjects/9.2/NET_Server_Doc/developer/samples/Web_Applications/ArcGIS_Buffer_Geoprocessing/ee05ac98-384b-4dbe-b8e3-8d741180dc9a.htm I'm practically drawing a circle around a point with a geoprocessing service. The radius is about 20 km. The code is working fine, but the "Clear" routine isn't. Every time I clear the map from the previous data, the memory of the buffer service isn't cleared too and I get this ![http://img814.imageshack.us/img814/101/errorrc.jpg](http://i.stack.imgur.com/IhTzp.jpg) The code of the "Clear" routine is the same of the sample, but isn't working:               // Clears features from all graphics layers in the resource specified by _graphicsResourceName     protected void ClearGraphics()     {         // Retrieve the resource and clear its graphics dataset         ESRI.ArcGIS.ADF.Web.DataSources.Graphics.MapResource graphicsMapResource =             Map1.GetFunctionality(_graphicsResourceName).Resource as             ESRI.ArcGIS.ADF.Web.DataSources.Graphics.MapResource;         graphicsMapResource.Graphics.Clear();          // This won't work too     GeoprocessingResourceManager1.GetResource(0).ClearState();          // Refresh the resource and copy the map's callback results to the callback results collection so     // the graphics are removed from the map     Map1.RefreshResource(graphicsMapResource.Name);     _callbackResultCollection.CopyFrom(Map1.CallbackResults);     }      This should be simple but I can't figure out how to solve this.. the only solution I found is to restart the IIS server which is obviously very bad Can someone please help me out?


title: arcgis server 10, vb.net, adding dynamic data

text: i am using arcgis server 10, with visual web developer 2010 express edition. i downloaded "ArcGIS_AddDynamicData_VBNet\Default_ADF.aspx" project form ArcGIS Add dynamic data for ArcGIS Server 10. everything works fine, now instead of the simple display they have in there i wanted to add the "web mapping application template". so what i really want is to have the same functionality shown in the project but with the default "web mapping application template". plz can you tell me how i could possibly integrate the default web mapping application template into that project? thanks a lot in advance!!!


title: What technique is used to acquire client location?

text: So far I know how to grab the client location to a application or web application by using the users geolocation or GPS location. I remember the possibility of finding the users approximate location by using their IP address. How is this done if you know how it is done can you link me to a tutorial or spec in addition to a short description? **Edit (Added Picture):**   ![enter image description here](http://i.stack.imgur.com/nMnNQ.png)    **Edit:** Possibly found through ISP?


title: Refreshing "layer render order" widget

text: i made changes in the render order of my map canvas using a :               qgis.utils.iface.mapCanvas().setLayerSet(a)      "a" being a list of **QgsMapCanvasLayer** My code is working as expected (i see result on map canvas) **but** the "layer render order" widget (dunno the right english name...) **is not refreshing** ... How can i programmatically do that ? Thanks !


title: How do I find the resolution of a scale range from the Arcgis Server REST services?

text: What's the best way to determine a LOD's resolution from the esri REST API? Currently, I make a request for a specific extent to the server from the flex api, and from that I divide the width of the image in pixels by the width of the map it represents in map units, but aside from feeling like a bit of a hack, it has trouble when dealing with services that are represented in degrees rather than meters.


title: WPS "Hello World" example in GeoServer?

text: I've been searching for a good example of how to publish a Web Process Service (WPS) on the GeoServer platform. The page where you'd think they'd document this process is actually pretty thin. There's some documentation on implementing an OWS service, but I am unsure what the conceptual and implementation overlap is between OWS and WPS in GeoServer. Is anyone aware of a good example of implementing a WPS on GeoServer?


title: How to group distinct layers together in ESRI Flex app?

text: I'm wondering if anybody has idea for how to best group distinct layers -- even layers of differing types like DynamicLayer, WFSLayer, etc. The goal in mind is that they could be visually represented together in a TOC component and they could be toggled on and off as a group in the TOC. The TOC component which comes with the out of the box Flex Viewer is setup to do this for ArcGISDynamicMapServiceLayers -- which have the layerInfo property that informs of the parent/child relationships. The question is -- has anybody done this before and if so, what is the best way to accomplish this? I was considering creating a LayerGroup object which extends the base Layer object and has a layerInfo property that stores a custom LayerInfo object -- similar to how the WMSLayerInfo is used to store parent/child relationships for WMSLayers.


title: How to drag a feature in Open Layers 3

text: I'm trying to replicate something that was fairly easy in Open Layers 2 in Open Layers 3, which is: dragging an existing feature to a new location on the map. In Open Layers 2 I just added a DragFeature control to my map like this:               this.mapRasterDragFeature = new OpenLayers.Control.DragFeature(_self.rasterVectorSource,          {'onComplete': function(e) { console.log(e);}});     this.mapRaster.addControl(_self.mapRasterDragFeature);     this.mapRasterDragFeature.activate();      But this control seems to have disappeared in Open Layers 3. I've searched far and wide in the documentation for the new way to this, but I can't seem to find it. I hope I'm overlooking something. I'm adding the feature through a click event on the map like this:               this.mapRaster.on('singleclick',function(e){         var iconFeature = new ol.Feature({             geometry: new ol.geom.Point(e.coordinate),         });         iconFeature.setStyle(_self.iconStyle);         _self.rasterVectorSource.addFeatures([iconFeature]);     });      All I've found looking like a drag interaction is ol.interaction.DragAndDrop(), but I think it does something different than just dragging and dropping a feature... Can someone help me out on this one?


title: ArcSDE register with gdb vs register as versioned

text: When wanting to update tables within a sde gdb, is there a difference between registering as versioned vs registering with the gdb vs registering with sde? I've read the question below and it hints that versioning is a subset of gdb registration. I am hearing that Unregistering a table as versioned will allow the use of sql to update the table, I assume directly from say SSMS. We are at ArcGIS for Server 10.1 and Microsoft SQL Server 2008 R2. Difference between registering with ArcSDE and registering with Geodatabase?


title: How to check ArcGIS REST Services Directory for date and time of publication?

text: Is there a way to check quickly within the service description when (date, time) the service had been published?


title: How to Krige with expectation of 0

text: I am trying to Krige in ArcGIS from point samples of a chemical concentration. Towards the center, most of the points are larger, and on the periphery, all of the points are 0. The generated raster looks good in the center and towards the edges, but beyond the outside points of zero, the raster values trend back up. I set the major range of the semivariogram to be 50, and the trending upwards is beyond thhis range of influence for any of the points that are greater than 0. So it seems that Arc is setting an expectation somewhere around the average of all the points. Is there are way to set the expectation to 0 in ArcGIS. If not, what is a simple way to Krige in R or QGIS where you can set the expectation and the parameters of the semivariogram? Here's a screenshot: ![Krigging](http://i.stack.imgur.com/SrKQd.png) See how near the periphery, especially upper left and lower right corners, the value trends back up, when it should be zero.


title: How to add a dynamic timeline to the QGIS canvas?

text: I am using the plugin Timemanager which works really great but there is one thing I'd like to improve. Instead of displaying the date and time on each frame, I would like to display a timeline stretching from *date_start* to *date_end* which shows the current frame time within this timespan. I understand this possibility is not (yet?) offered by the plugin and therefore I am looking for an amazing idea on how to implement this. A core Qgis functionality that could be a workaround? An already existing plugin that I have missed? A small python script? Obviously, in order for this to work nicely with the timemanager plugin, this timeline should be displayed in the map canvas, not in the composer. A solution for Qgis 2.0 would be nice but I could also live with something for 1.8. If someone has an idea or (even better) has already implemented something similar, please raise your hand!


title: Use of Copy / Paste when editing attribute tables in QGIS

text: I have a number of polygon attribute tables with columns that I want to update with new information. Is there a way of (relatively) quickly editing info in attribute tables in QGIS eg such as using Control C and V such as in MapInfo. Currently seem to need to double click/right click/copy, move to cell, double click/right click/paste each individual cell I want to change - would like something simple rather than going down the field calculator, filter, etc route. With points no problem; I usually throw them out to Excel, edit and then bring back in again. Any pointers appreciated.


title: unable to modify records from one shapefile to another

text: I am trying to update the atrribute data of a feature class with respect to another one by matching a common field. But it is returning an error for the setValue command. My code looks like:               import arcpy     fc="D:/Geonames_CA.gdb/Export_Output"     fs="D:/abc.gdb/Export"     rows = arcpy.SearchCursor(fc)      fldList=arcpy.ListFields(fc)     rws = arcpy.SearchCursor(fs)     flList=arcpy.ListFields(fs)     for row in rows:         for fld in fldList:             a = row.getValue("PLACE_NAME")             a_xmin=row.getValue("XMIN")             a_xmax=row.getValue("XMAX")             a_ymin=row.getValue("YMIN")             a_ymax=row.getValue("YMAX")     for rw in rws:         for fl in flList:             b = row.getValue("PLACE_NAME")             b_xmin=row.getValue("XMIN")             b_xmax=row.getValue("XMAX")             b_ymin=row.getValue("YMIN")             b_ymax=row.getValue("YMAX")             if b==a:                 cursor = arcpy.UpdateCursor(fs)                 for r in cursor:                     r.setValue(b_xmin, a_xmin)                     r.setValue(b_xmax, a_xmax)                     r.setValue(b_ymin, a_ymin)                     r.setValue(b_ymax, a_ymax)                     cursor.updateRow(r)                 del r, cursor      I am only getting this error: > Traceback (most recent call last): File "C:\Python27\ArcGIS10.1\Lib\site- > packages\pythonwin\pywin\framework\scriptutils.py", line 326, in RunScript > exec codeObject in **main**._ _dict_ _ File "D:\Script2.py", line 25, in > r.setValue(b_xmin, a_xmin) File "C:\Program > Files\ArcGIS\Desktop10.1\arcpy\arcpy\arcobjects\arcobjects.py", line 1007, > in setValue return > convertArcObjectToPythonObject(self._arc_object.SetValue(*gp_fixargs(args))) > RuntimeError: ERROR 999999: Error executing function. Please help..