title: Should 2D Positions and 3D Positions be considered conceptually different entities?

text: Let's suppose we want to model GeographicPositions, described by coordinates or coordinate tuples inside some Coordinate Reference System. Also, these would be put inside some sequence to represent Lines (Polylines, Tracklogs, etc.) If we consider a Mapping Software Application, where we can draw over a map image with the mouse (such as Google Maps, OziExplorer, etc.), we are defining horizontal positions based on Latitude and Longitude, so these positions are inherently 2D. Now if we consider a GPS device, it tracks spatial positions based on Latitude, Longitude and Elevation, so these positions are inherently 3D. We could, of course, use DEM to get terrain elevation based on horizontal position, but that should not (as I understand) be considered an INTRINSIC PROPERTY of that position, since it is defined by an external function (an operation over the coordinates AND the DEM dataset). In the other hand, a GPS device can be moved not only along the terrain surface, but it can be inside an airplane, and thus its elevation is intrinsically defined and independent of horizontal position, not being "extractable" by any operation over DEM dataset. So, my questions are:   * Is it common practice in GIS Systems to treat 2D and 3D positions (and collections of positions) as conceptually different entities?   * Should a 2D position be considered a limited type of a more general 3D position?      * With Elevation set to None/Null?     * With Elevation set to 'NaN'?     * Without Elevation property at all?   * Should a 3D position be considered an extended type, derived from a more basic 2D position by adding an Elevation property?   * How to represent a 3D tracklog where some of its positions' elevations are missing? How should it be mathematically processed regarding its missing values? Should a tracklog instance like this exist in the first place? I would be glad to update the question if anyone thinks it must be clarified!


title: XML file created when .mxd is saved

text: I have ArcGIS 10.1 and Xtools 10.0, each time .mxd is saved there is a .xml file created. How do you prevent the creation of the .xml file with the same name as the .mxd? Where in the Xtools menu is that an option?


title: Determining the Coordinate System of shapefiles

text: I have 8 shapefiles which have a common coordinate system but none of them have it defined i.e. "unknown" is written when I look in the properties dialog box. However, I can see the extent which I am sharing here: Top: 42.349145 Bottom: 41.710962 Right: -82.430463 Left: -83.126193 I need to define the projection of these shapefiles so that I could do further analysis. I know that the features these shapefiles represent are located in Windsor, Ontario i.e.they should have a projection of WGS 84 UTM Zone 17N. Can I directly define this projection without knowing about the geographic coordinate system?


title: error when running a script for clipping a group of rasters on QGIS 2.2 OS X 10.9

text: I'm trying to clip a group of raster files on QGIS 2.2 (KyngChaos) running on OSX 10.9. When I run the script I get an error related with the location of the files. I'm wondering if this is something related with the way OSX locate files or if it's something related with Python. Here is the script (taken from _fromgistors.blogspot.com_ )               # input variables     # where to save clipped rasters     outputDir = "/Volumes/Win2TB/SIG/Azuero/LC80120552014029LGN00_converted_crop"     # upper left point     UX, UY =  504259 , 915459     # lower right point     LX,  LY = 628983 , 791491          # actual script     import os     lddLrs = qgis.utils.iface.legendInterface().layers()     for lyr in lddLrs:      if (lyr.type()==QgsMapLayer.RasterLayer):       os.system("gdal_translate -projwin " + str(UX) + " " + str(UY) + " " + str(LX) + " " + str(LY) + " -of GTiff " + str(lyr.source()) + " " + outputDir + "/" + str(lyr.name()) + ".tif")       qgis.utils.iface.addRasterLayer(str(outputDir + "/" + lyr.name() + ".tif"), str(lyr.name() + "_clip.tif"))      This is the error I'm getting: GDAL provider: Cannot open GDAL dataset /Volumes/Win2TB/SIG/Azuero/LC80120552014029LGN00_converted_crop/RT_LC80120552014029LGN00_B7.tif: `/Volumes/Win2TB/SIG/Azuero/LC80120552014029LGN00_converted_crop/RT_LC80120552014029LGN00_B7.tif' does not exist in the file system, and is not recognised as a supported dataset name. Raster layer: Provider is not valid (provider: gdal, URI: /Volumes/Win2TB/SIG/Azuero/LC80120552014029LGN00_converted_crop/RT_LC80120552014029LGN00_B7.tif It seems that QGIS is looking for the files in the folder where I want them to be saved. Any help is appreciated. Thanks


title: Using Multiple ShapeFiles as a single address locator

text: Let's say I have 5 shapefiles, and I want to use them as a single address locator via the ILocator interface. What would you recommend? 1) Merge the shapefiles into a single shapefile then convert to a feature class?   2) Merge multiple feature classes into a single feature class?   3) Get multiple locators and merge them into a single locator? C# VS2010 ArcEngine 10.


title: Overlap between Google Earth and Topographic OSM

text: First, I have Google Earth and Garmin BaseCamp with OSM maps. I would be satisfied if the solution is only available on one of these software. A third software that combine both view can be interesting if I can export to one of these. Today I can draw a route on pedestrian road by using the routing of BaseCamp and I can as well export it to Google earth. But I did not find any option to merge both views. Is it possible to have all pedestrian road with satellite view on any of these software and how?


title: Displaying point data to reflect time/age

text: I have a layer of 3,000 points for a state, with each indicating the age of a specific type of building. What is a recommended visual display that effectively and equally highlights the age of these structures. I thought of using a heat map, but that could lead to some bias, as in a small vicinity you can have 3 new buildings and 1 old one or vice versa taking away from the truthfulness. However, I'd like a visual aid that when looked at briefly or a snapshot can tell a story of the age of these structures.


title: Why are the coordinates transformed with proj.4 off compared to http://valdemar.kms.dk/trf/?

text: I am trying to transform from epsg:4326 (wgs84) to epsg:25832 (utm ETRS89) using Proj.4. I am comparing the results with results from http://valdemar.kms.dk/trf/ . it seems like the results calculated by proj.4 is off by 2m northing and 4m in easting in Copenhagen. The results farther west are more correct. Am i missing some sort of correction attributes in my calculation?


title: Using Terraformer to convert to WKT

text: I'm using Terraformer to parse some ArcGIS JSON into a Terraformer primitive. I'd like to store the geometry in PostGIS for other purposes. So I'm trying to use the WKT.convert tool to convert this primitive into Well Known Text. Here is the JSON primitive given to me by ArcGIS.parse():               {         "type": "Polygon",         "coordinates": [             [                 [                     -111.53992413334949,                     40.54679752920996                 ],                 [                     -111.27075909428711,                     40.54471043695149                 ],                 [                     -111.32019757084959,                     40.43400149811995                 ],                 [                     -111.21033428959964,                     40.34195348225386                 ],                 [                     -111.35864971928709,                     40.222522475606446                 ],                 [                     -111.56464337163074,                     40.31263909907781                 ],                 [                     -111.47400616459953,                     40.40890997045407                 ],                 [                     -111.534430969287,                     40.54053605732489                 ],                 [                     -111.53992413334949,                     40.54679752920996                 ]             ]         ],         "bbox": [             -111.56464337163074,             40.222522475606446,             -111.21033428959964,             40.54679752920996         ]     }      When I calle WKT.parse(), I get:               Error: Unable to parse         at Error (<anonymous>)         at Object.parse (/Users/jim/Sites/utilisync-server/node_modules/terraformer-wkt-parser/terraformer-wkt-parser.js:571:13)      That simply refers to this line in terraformer-wkt-parser:               function parse (element) {         var res, primitive;              try {           res = parser.parse(element);         } catch (err) {           throw Error("Unable to parse", err);         }              return Terraformer.Primitive(res);       }      There's nothing left in the stack trace to go on. Is there a better way to get PostGIS geometry from ESRI JSON? (that can be implemented in code, NOT a program that I'd have to run) Does anyone see anything in my JSON that should cause the WKT parser to choke? It looks fine to me.


title: How to return FeatureLayer from gp service(using custom function tool) (.NET)?

text: My basic problem is that I have a custom function tool HLZCalulation that I created. It takes in a string to a file gdb , opens it, does some calculations, and then returns a featurelayer with the results of said calculations. The tool works as both a simple tool and a model, returning the correct features and not reporting any errors. However, when I publish the model as a service and try to run it from ArcMap, it reports success message from the server but the output feature class is empty. When I check the server job folder, it has created the shapefile with the features inside it but fails to return the featurelayer without throwing any specific errors (the REST portal returns success message, then a fail message – see error message below). Packing the parameter is done at the end of execute (once the shapefile is created) rather than validate as with the calculate area sample on the esri site, which _could_ be a problem. The output parameter definition is included below, as well as the lines of code for packing the parameter, and the output from running the rest portal(in browser) any ideas? //////////////////////////////////////////////////////////PARAM DEFINITION                           // Output parameter (Derived) and data type is GPFeatureLayerTypeClass                 IGPParameterEdit2 outputParameter = new GPParameterClass();                 outputParameter.DataType = new GPFeatureLayerTypeClass();                 // Value object is GPFeatureLayerClass                 outputParameter.Value = new GPFeatureLayerClass();                 // Set output parameter properties                 outputParameter.Direction = ESRI.ArcGIS.Geoprocessing.esriGPParameterDirection.esriGPParameterDirectionOutput;                     outputParameter.DisplayName = "HLZ Output FeatureClass";                     outputParameter.Name = "hlz_featureclass";                     outputParameter.ParameterType = ESRI.ArcGIS.Geoprocessing.esriGPParameterType.esriGPParameterTypeDerived;                     parameters.Add(outputParameter);      //////////////////////////PACKAGE PARAM                   IGPValue resultFeatureLayer = m_gpUtilities.CreateFeatureLayer(resultFeaturePath + ".shp");           IGPParameter parameter = (IGPParameter)userParamValues.get_Element(OUTPUT_PARAM);                 m_gpUtilities.PackGPValue(resultFeatureLayer, parameter);      //////////////Rest Portal output esriJobMessageTypeInformative: Submitted. esriJobMessageTypeInformative: Executing... esriJobMessageTypeInformative: Executing (Model): Model # D:\data\test\TGD_CLIP.gdb esriJobMessageTypeInformative: Start Time: Mon Oct 11 17:29:18 2010 esriJobMessageTypeInformative: Executing (HLZ Calculation): HLZCalculation D:\data\test\TGD_CLIP.gdb # esriJobMessageTypeInformative: Start Time: Mon Oct 11 17:29:18 2010 esriJobMessageTypeInformative: Intersect completed, should return esriJobMessageTypeInformative: Executed (HLZ Calculation) successfully. esriJobMessageTypeInformative: End Time: Mon Oct 11 17:32:08 2010 (Elapsed Time: 2 minutes 50 seconds) esriJobMessageTypeInformative: Executed (Model) successfully. esriJobMessageTypeInformative: End Time: Mon Oct 11 17:32:08 2010 (Elapsed Time: 2 minutes 50 seconds) esriJobMessageTypeError: Failed.


title: Greek labels in tilemill

text: I have a shapefile with greek characters. When i import it in tilemill, they appear unreadable. Is there any way to show them properly?


title: How to include external libraries with my QGIS plugin?

text: I'm developing a plugin which makes use of " **networkx** " library extensively. However, networkx is not part of Python libraries pre-installed in QGIS. Is there a way to include this library with my plugin installation, so that my plugin can be portable? Thanks!


title: Converting xyz batch to ESRI grid format using XYZ2GRID 2.1 - batch problems

text: I am trying to use Min-Lang Huang's XYZ2GRID2.1-Batch program downloaded from ArcScripts and I keep getting the error C:\users\public...........xyz Error-Object required Does anyone know what the problem is? I have created an empty ASCII folder and an empty GRID folder and all the rest of it but not sure what the problem is. Does this script need to sit anywhere special? How does it know that my spatial extension is activated? I am using ArcGIS 9.3.1 and Windows vista. Any help is appreciated thanks DG


title: Raster To Ascii cmd-line tool

text: Where can I find a tool like ArcGis: Raster_to_ASCII Provided that the tool **MUST** be accessible through command line interface.


title: How to add pins (graphics) to ArcGIS map in wpf?

text: I have a wpf app with an ArcGIS map, I used the extent property to zoom and pan to Canada (the location I want) and added the pins (images) by adding a graphics layer, see code below. currently I'm using this map: http://server.arcgisonline.com/ArcGIS/rest/services/ESRI_StreetMap_World_2D/MapServer I'd like to use this instead: http://server.arcgisonline.com/ArcGIS/rest/services/World_Topo_Map/MapServer the problem is that the pins show up in the middle of the map, the 0,0 location, when I switch to the new map. Are graphics specific to a type of maps or am I adding graphics the wrong way? What's the best way to add images to map? Thanks, T               <Grid x:Name="LayoutRoot">         <Grid.Resources>             <esri:PictureMarkerSymbol x:Key="PushPin"                                       OffsetX="11"                                       OffsetY="39"                                       Source="pushpin.png" />         </Grid.Resources>              <esri:Map x:Name="MyMap"                   Extent="-125, 40, -60, 60">             <esri:Map.Layers>                 <esri:ArcGISTiledMapServiceLayer ID="StreetMapLayer"                                                  Url="http://server.arcgisonline.com/ArcGIS/rest/services/ESRI_StreetMap_World_2D/MapServer" />                 <esri:GraphicsLayer>                     <esri:GraphicsLayer.Graphics>                         <!-- Halifax -->                         <esri:Graphic Symbol="{StaticResource PushPin}">                             <esri:MapPoint X="-63"                                            Y="45" />                         </esri:Graphic>                              <!-- Toronto -->                         <esri:Graphic Symbol="{StaticResource PushPin}">                             <esri:MapPoint X="-79.5"                                            Y="44" />                         </esri:Graphic>                              <!-- Vancouver -->                         <esri:Graphic Symbol="{StaticResource PushPin}">                             <esri:MapPoint X="-123"                                            Y="49.5" />                         </esri:Graphic>                     </esri:GraphicsLayer.Graphics>                 </esri:GraphicsLayer>             </esri:Map.Layers>         </esri:Map>     </Grid>


title: Is there any way to tell whether a park from GeoNames is small or big?

text: I have this query for parks in Alabama for example: http://api.geonames.org/search?featureCode=PRK&username=genadinik&country=US&style=full&adminCode1=AL&maxRows=100 But it returns all sorts of parks like ballparks, playground, national forests. Is there a way to distinguish city-block parks from stadiums and from big parks for outdoors? Also, there is a value that appears named score like this: 34.15653610229492 What is that value? Thanks!


title: Using Python to Switch Layout Templates

text: Using ArcMap 10.1: Since there is not a way to control portrait or landscape of the layout view through python, is there a way to switch layout templates without opening a new mxd through python? Since there is a change layout template button it would be nice to have a python function to do the same..


title: What is the most efficient way of displaying GPS in real time using OpenLayers?

text: i'm would like to display GPS traces (not only a position) in realtime on OpenLayers. What is the most efficient way to do it? Sending KML or GPX files to browser or sending data via GeoJSON? br d


title: How to Query Parks in GeoNames and Combine it With a Polygon Layer?

text: I run a hiking site and would like to get parks where you can at least have a nice walk. That means they should not be playgrounds, stadiums, and probably be bigger than 3 square city blocks at least. Is there a way to get this information from GeoNames? I can get parks, but how do I combine it with the polygon size? Thanks!


title: QGIS: GPS Tracking with improved Altitude

text: I am preparing for archeological fieldwork in Sudan. I recently came up with RTK equipment and software using RTKlib to export the real time corrected data (cm level) to QGIS via a simulated NMEA device. However I checked the NMEA data flow and it includes elevations down to mm; in QGIS only one decimal is given. Is there any way to improve this and have the altitude information shown in cm as well? Also QGIS applies geoid corrections but I wish to use ellipsoidal elevation. And finally right now I have to include the elevation data manually in the datamask for every point while saving. Is there a way QGIS can save the elevation in points and polylines?


title: GeoServer Layer Preview Stopped Working

text: I got GeoServer up and running today and created a store that is linked to my postgis database, with a layer corresponding to one of the tables in the database. I then previewed the layer using the GeoServer Layer Previewer and it all worked nicely. I began messing around with GeoServer, trying to get this layer plotted on my own OpenLayers map, and now none of the layers work in the Layer Previewer. When I click "go", a new tab opens with a blank box where the map should be, and the tab indicates that it is loading. However, it does not ever appear to load correctly. Does anybody have any idea what I might have done to mess this up? Thanks, -Rob


title: Map disappear after a click button event

text: I have a problem in my project that uses both **openlayers** and **asp.net**. The problem is that my map disappears when I click in the button on my web interface. Anyone already encountered like this problem?


title: GetFeature request giving variable results

text: I have a map with a tilecache baselayer and a wms property layer. By clicking on a property a GetFeature request is done using OpenLayers.Protocol.WFS.fromWMSLayer for the protocol. A test version of the map is at http://geonet.allerdale.gov.uk/testmapapp/allerdalelookup/?uprn=100110685115 When clicking on Property 5 Market Place (the one highlighted on load) a list of possible properties for that feature should appear. Sometimes this list shows as: ![enter image description here](http://i.stack.imgur.com/lpisZ.jpg) However on refreshing the page and clicking on the same feature it gives me: ![enter image description here](http://i.stack.imgur.com/ndoVP.jpg) The GetFeature request I am using is                 var control = new OpenLayers.Control.GetFeature({                        protocol: OpenLayers.Protocol.WFS.fromWMSLayer(llpg_layer,                         {  outputFormat: "JSON",                            readFormat: new OpenLayers.Format.GeoJSON()                        }),                        single: false,                     clickTolerance: 10,                     filterType: OpenLayers.Filter.Spatial.CONTAINS                   });                control.events.register("getfeatureinfo", this, function(e) {                feature_layer.removeAllFeatures();                feature_layer.addFeatures([e.feature]);                      showInfo(e.feature);                        });             control.events.register("featuresselected", this, function(e){                var features = e.features;                  if (features.length > 1){                     multipleFeatures(features);                  }             });             control.events.register("featureunselected", this, function(e) {                 feature_layer.removeFeatures([e.feature]);             });              map.addControl(control);           control.activate();      From what I can see on firebug the first request posts this:               <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.1.0" outputFormat="JSON" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">     <wfs:Query typeName="AllerdaleOnly:abc_llpg">     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">     <ogc:Contains>     <ogc:PropertyName></ogc:PropertyName>     <gml:Envelope xmlns:gml="http://www.opengis.net/gml">     <gml:lowerCorner>300587.38521012 528552.83828508</gml:lowerCorner>     <gml:upperCorner>300588.26715409 528553.72022905</gml:upperCorner>     </gml:Envelope>     </ogc:Contains>     </ogc:Filter>     </wfs:Query>     </wfs:GetFeature>      And the second one posted this:               <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.1.0" outputFormat="JSON" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">     <wfs:Query typeName="AllerdaleOnly:abc_llpg">     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">     <ogc:Contains>     <ogc:PropertyName></ogc:PropertyName>     <gml:Envelope xmlns:gml="http://www.opengis.net/gml">     <gml:lowerCorner>300587.38521012 528544.23051142</gml:lowerCorner>     <gml:upperCorner>300588.26715409 528545.11245538</gml:upperCorner>     </gml:Envelope>     </ogc:Contains>     </ogc:Filter>     </wfs:Query>     </wfs:GetFeature>      **UPDATE** I have noticed the following I have added mousePosition to a local version of the page, and when the feature (5 Market Place) selects as expected the top left and bottom right as follows: Top Left: 300578.68/528556.80 Bottom Right: 300591.56/528550.80 When on reload the feature selects is not working as expected, the 2 corners' position has changed. Top Left: 300578.86/528548.57 Bottom Right: 300591.38/528542.75 The following function is called in the initial map creation function, adding a highlighted feature and centring and zooming the map on that feature:                 var searchBLPUs = function (uprn){                   var filter = '<Filter xmlns:ogc="http://www.opengis.net/ogc">';                 filter += '<PropertyIsEqualTo>';                 filter += '<PropertyName>a_uprn</PropertyName><Literal>';                 filter += uprn.trim();                 filter += '</Literal></PropertyIsEqualTo></Filter>';                      OpenLayers.Request.GET({                         url: that.wfsurl,                         params: {                                 typeName: "AllerdaleOnly:abc_llpg",                                 service: "WFS",                                 version: "1.1.0",                                 outputFormat: "JSON",                                 readFormat: new OpenLayers.Format.GeoJSON(),                                 srsName: "EPSG:27700",                                 request: "GetFeature",                                 filter: filter                         },                         success: function(reply) {                             var format = new OpenLayers.Format.GeoJSON();                             var feature = format.read(reply.responseText)[0];                             that.centerOnFeature(feature.geometry);                             var address = feature.attributes.b_address;                             showAddress(address, feature.attributes.a_uprn);                         },                         failure: function(reply) {                                 alert("failed");                         }                 });        }


title: How can I create a button to add a POI layer?

text: I am making a map of a city with openlayers. I have the shapefile data of the POI (hospital,bank,school etc.) I want to add the layers with buttons not in layerswitcher. Like whereis.com's quick search buttons ... Can someone please help?


title: Custom Form with QTDesigner: Combine signals and slots?

text: I dare to ask that question here and not in the programmers section, cause i'm actually not looking for a coding answer, but a solution through the QtDesigner GUI (if possible). Custom QT-forms combined with QGIS can be a really powerful solution for data- entry tasks, so i try to tweak things further, but already reached my limits. ![my form](http://i.stack.imgur.com/ZKNug.jpg) Little explanation on the form: I have lists of bird observations with start and end times. I already connected some time fields through signals and slots, so when changing the start field, the end fields take over that value as well (only have to change the minutes afterwards -> less typing -> less errors) When an entry is finished and the next item is digitized, the previous values will be loaded as default. As most of the time intervals are contigous, i would like to be able to copy the former end-value to the new start-value when pressing the button E -> A. In the example that means, that the reloaded end-value 15:03 would be inserted instead of 15:00. It's quite easy to send values form here to there via signals and slots, but i didn't understand so far if it's possible to trigger such an action also with a third partner (the button)


title: Allocate places according to distance

text: I've asked this over on Stack Overflow but thought that someone may have a GI point of view that may help: As an over-simplified example I have a list of events that have a maximum attendance:                event    | places     ===================      event_A  |   1      event_B  |   2      event_C  |   1      And a list of attendees with the distance to the events:                attendee    | event_A dist | event_B dist | event_C dist     ==========================================================      attendee_1  |      12      |      15      |      12            attendee_2  |      11      |      15      |      11      attendee_3  |      10      |      11      |      12      Can anyone suggest a (hopefully) simple method to produce a set of options providing the best case allocations based on shortest total distance and also on shortest mean distance? I currently have the data held in Oracle Spatial database, but I'm open to suggestions. In my actual data there will be more places than attendees, so the events do not have to be full, but they cannot be "over-full". For the possible solutions I'd like 2; one showing the lowest total distance (i.e. sum of all attendees' distances), the other showing the mean distance (i.e. sum of all attendees' distances / No. of attendees) _Update:_   * Each atendee should be assigned to exactly one event    * Each event has a limit as to how many atendees are assigned to it    * Underfull or even empty events are no problem    * Each assignment between an event and an atendee corresponds to a given distance    * I want to minimize the overall distance for all assignments


title: How to read the app.config settings in ESRI ArcObjects ArcMap Add-In

text: I am developing custom toolbar and buttons for ArcMap 10.1 utilizing VS2010 and ArcObjects SDE 10.1 extensions. I put some application settings in the app.config file and I want to make them available in the runtime. When I execute               private static readonly string appSetting1 = System.Configuration.ConfigurationManager.AppSettings["appSetting1 "];      the value is null. When I execute               AppDomain.CurrentDomain.SetupInformation.ConfigurationFile      the value is > "C:\Program Files (x86)\ArcGIS\Desktop10.1\bin\ArcMap.exe.Config" Should I put my application settings is the ArcMap.exe.Config in order make them available to the ArcMap 10.1 custom toolbar or there is another way to do this? Thanks


title: Get coordinates from current view using arcobjects

text: I'm wondering if there's a way to get the longitude/latitude coordinates from the center of a data frame? The problem I'm faced with is this: I want to set a bookmark and have ArcMap tell me the coordinates of the center of that bookmark. (Without using a mouse. This all needs to be automated.) This is akin to hovering one's mouse over the center of a map and reciting the long./lat. numbers that are shown at the bottom of the ArcMap window. I'm pretty "green" when it comes to VBA programming in ArcMap 10.0. At first I thought I could utilize the "Info" tool, but now I'm guessing that I need to get the use the IEnvelope interface to calculate the center of the map, then grab the coordinates from that. Any ideas?


title: Help with formatting ArcGIS text elements with Python

text: I'm struggling with a formatting issue that comes from inserting text into a Text Element object in ArcGIS using Python. I want to print a table of values of which all but the last column are numbers. My python script loops through a list and applies formatting to each element like this:               #Text element is referenced:     tableTextQuad1 = arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT", "tableTextQuad1")[0]     #All of the work occurs within a function that takes the table, x postion, y position, and text element reference (txtBox1) as its arguments)     #A clone of the element is created and used by the newCellText     newCellText = txtBox1.clone("_clone")       #A for loop is used to loop through lines in a table, and each line is formatted like so:     reformat1 = '{:>3}      {:>6.2f}      {:>3}      {:>7.2f}      {:>7.2f}     {:<10}'.format(t[0], t[1], t[2], t[3], t[4], t[5])     myList.append(reformat1)     myList.append('\n')     #To get rid of commas I do the following:        alpha = "', '".join(myList).replace("', '", "")     #Then the alpha string is sent to the cloned text object         newCellText.text = '{:}'.format(alpha)     newCellText.elementPositionX = x1     newCellText.elementPositionY = y1      When I don't include the last element (t[5]) the table looks great...like this: ![enter image description here](http://i.stack.imgur.com/iprk4.jpg) But as soon as I try adding the string t[5] (which are directions like North, South...) to the line everything in the line gets shifted and seems to ignore the width and character spacing allocation I set in .format(). The result looks like this: ![enter image description here](http://i.stack.imgur.com/5iw5O.jpg) I've tried many different formatting parameters for t[5] but nothing seems to work the way I want it. Can anyone tell me why all the numbers would adhere to formatting rules like right/left justify and width but the inclusion of letters on the same line messes it all up? What is the correct formatting I should use? Many thanks in advance for help! *Update - I have also tried adding a carriage return to the \n, so that '\r\n' is appended to each line, but this does not fix my problem. Could it be a bug?


title: How to remove shapefile polygon vertex marker points in QGIS 1.8.0

text: I have created a shapefile layer and drawn a polygon to show the boundary of a woodland. How can I edit the polygon to remove the crosses that denote the intersection of two points on the boundary?


title: IBufferConstruction create buffer with different geometry types (C# ArcObjects)

text: I am using the following method in C#               public void ConstructBuffers (         IEnumGeometry pInputGeometries,         double distance,         IGeometry CollectionpOutputBuffers);      to create buffers. In the method above if I have a `IEnumGeometry` with same types of multiple geometries I can dissolve the overlaps to get a single geometry. But when I have a `IEnumGeometry` with different geometry types (points, lines and polygons) and try to dissolve the overlaps I don't get a single geometry but my output (`CollectionpOutputBuffers`) is the same number of geometries as input. Anyone has any idea what may be wrong.


title: Specify a field in attribute table with raster calculator

text: How can I specify the attribute field to use in the raster calculator of ArcMap 10.2? In 9.x, I used the dot notation. i.e. To use the "New_Value" field, I could do the following: [MyRaster].New_Value + 1 = new raster This doesn't appear to work in 10.2?


title: What programs would allow for the mapping of a geoid in 3d?

text: I would like to find a method for plotting an exaggerated geoid (e.g. topography): ![enter image description here](http://i.stack.imgur.com/0xarn.jpg) I primarily use R for my work but would be interested in other open source software options (e.g. QGIS w/ GRASS). I found this wonderful tool for MATLAB (link) but haven't discovered other options.


title: Editor tracking on geometric network feature classes

text: Editor tracking is not working on my geometric network feature classes. Although I enabled the editor tracking - dates and times are not recorded in layers while editing. Is it impossible? or a local problem with my geodatabase? ![EditorTrackingEnabeld](http://i.stack.imgur.com/vWfiB.png)


title: How to edit attribute tables using Python/ArcPy?

text: I am new to Python, so sorry for such a trivial question. There is a shapefile (let's call it Structures.shp) and the script I'm writing needs to update the value of only one record in only one field (let's call the field structuretype). The field is the subtype field and the default value is (for example) 3, but the script needs to change the value for this particular record to 4. What is the easiest way for the script to update the value of this one cell?


title: Display list of OpenLayers Features which are currently visible on screen

text: I have an OpenLayers map a few thousand Points of Interest. These POIs are Vector Features rendered in a Vector Layer. We are loading the POIs once in order to minimize load on the server. Next to the map will be a box which will display the names of POIs currently in view. I'm not clear on the best way to update this list of names. I believe that I can store an array of the POIs (Features). Then, when the map is panned/searched/zoomed, loop through them, calling the "onScreen" function.               map.events.register("moveend", map, function(evt) {           var output = [];           for (var i=0;i<myFeatures.length;i++) {               var poi = myFeatures[i];               if (poi.onScreen()){                  output.push(poi.attributes['name']);               }           }          // then, use output to display the list in html...          myDOMbox.innerHTML = output.join();          });      Is there a more efficient way to do this? Thanks in advance.


title: Running python script outside of Arcmap's python console

text: My codes works fine in Arcmap's python windows. But when I try to run it in a command promt windows, It says that;               ERROR 000837: The workspace is not the correct workspace type.      But it works fine in Arcmap's python console. My code is:               path_to_sde = "Database Connections\\my_conn_file.sde"          arcpy.ReconcileVersions_management(path_to_sde, "ALL_VERSIONS", "SDE.DEFAULT", "USER.ADMIN", "LOCK_ACQUIRED", "NO_ABORT", "BY_OBJECT", "FAVOR_TARGET_VERSION", "POST", "KEEP_VERSION", Reconcile_Versions_Log__2_)      PS: I tried full path of sde file instead of "Database Connection/filename.sde" Do you have any suggestion?


title: registering or unregistering as versioned a future dataset with python

text: As you know, you can store feature classes in arcsde in two ways. **Firts** , in the root path: sde_conn_file.sde/FeatureClass **Second** , in a dataset: sde_conn_file.sde/ **DataSet** /FeatureClass How can I register as a versioned the feature class that given **second** example? I tried the code given below but, it doesn't work.                   FCList = arcpy.ListFeatureClasses("*","All","DATASETNAME")         for fc in FCList:             arcpy.RegisterAsVersioned_management(fc,"NO_EDITS_TO_BASE")


title: Data-defined Styles in QGIS

text: I have downloaded a shapefile of bedrock geology. I need it to be displayed as it should be with the style applied. The supplier will no give me a qml. but i have noticed in the .dbf file the Cyan Magenta yellow and black numbers are stated. There are hundreds of entries and i dont have the time or patience to go through each one individually. Is there a way i can apply the colours onto the shapefile (which is polygons and the extracted csv file creates points which have no xy). QGIS1.8.0 on windows 7.


title: QGIS join data not transferring

text: I'm reanalyzing data from a VA Late Woodland Palisaded Village. The site is mapped in CAD and exported as a shapefile. QGIS brings up the 1059 features as polygons. I am using a CSV file for the attribute data for the decorated pottery. I've joined the vector layer with the attribute file. All of the headings are there, but not all of the data is present. It appears that where any feature has more than one entry, the attribute table stops at the first one. Is there a way to get it to list all of them? The DB is organized into: Feature # in which not all have decorated pottery, ID# for each piece in each feature that has decorated pottery and then lists the various attributes. Thanks in advance.


title: Transform GeoRSS from WGS84 to UTM 33N with ogr2ogr

text: Is this really not supported or do I have missed something in my code? If it is not possible with ogr2ogr, do you know of any other open source methods to reproject from GeoRSS to GeoRSS or from shp to GeoRSS?               srcfile=D:/GIS_DataBase/QGIS_Projects/BIK-Stm/BIK_STAINACH_new.xml     destfile=D:/GIS_DataBase/QGIS_Projects/BIK-Stm/BIK_STAINACH_UTM.xml          ogr2ogr -f "GeoRSS" $destfile -t_srs "EPSG:25833" $srcfile          echo check result:     ogrinfo -al $destfile          $ sh reproject.sh     ERROR 6: For a non GML dialect, only WGS84 SRS is supported     ERROR 1: Terminating translation prematurely after failed     translation of layer georss (use -skipfailures to skip errors)     ...


title: ArcGIS10: Publishing GPtools etc. in AddIns

text: I'm creating an ArcGIS 10 add-in in .NET C# which contains some tool commands and an extension class. This extension is being launched on application start- up. If I create any custom GP Tools, custom symbology renderers, custom toolboxes etc. in code do these still need to be registered during an install or can I launch them (as in have a custom symbology class in a layer properties dialog) dynamically from the extension class during run time? Assume that the user of the tools does not have admin rights to their machine. Hopefully I've made myself clear if not just ask! Thanks Rob


title: Georeference images together

text: I have taken a number of images from a very large rock overhang that contains ancient paintings. I want to piece the images together to create a single image of the entire rockshelter. Through ARCGIS I tried to do this by georeferencing the images, which have no spatial data. It worked quite well. However, I don't know how to save it so that when I move onto the next image, the previous georeferencing remains. How can I save it?


title: CartoDB: Question about map loading indicator

text: We are using the CartoDB.js API to create our visualizations at runtime using a pre-exsiting Google map. Everything seems to be working fine so far, except that we would like to show some sort of visual loading indicator while the cartodb layer tiles are loading to let the end user know that the map is still rendering. I see that a visual loading indicator is displayed on the (google) map when I use the cartodb.createVis() function with a pre-defined visualization, however I do not see how I can access/use this with the API for my map visualization that I create at runtime. Does anyone know how to display a loading indicator on the google map or have an example of doing so? Thanks, Greg


title: Whats going on with Mapinfo these days?

text: I have not touched Mapinfo for a good 2-3 years. (Think it was 7.5 back then, around the time it got aquired by Pitney Bowes) I found it a great package, particularly with Vertical Mapper and utilising MapBasic to automate tedious tasks. I would like to get an update on peoples recent experiences with the product. I saw some nice screens of v10 where the interface looked a lot shinier. Has it changed much over the past few years? What about the web side of MapInfo, was it Map Extreme or something? I am pretty sure there are no similar web platforms like what ESRI offer (Silverlight/JS/etc) - Are there any plans to catch up on the web side of things?


title: Programmatically access HTML Popup in Layer Properties ArcGIS 10.2

text: Is it possible to access HTML Popups (found in layer properties) programmatically in ArcGIS 10.2? I can't find any reference but thought I'd ask to double check.


title: Using URL parameters in ArcGIS Server 10 Viewer for Flex

text: Anyone know if you can use URL parameters to turn on layers. I have had success using URL parameters to go to specific coordinates.


title: 2010 Census Tiger County shp file with defined coastal land area

text: I am looking for a 2010 Census Tiger Shape file that hugs the coastline and doesn't include the ocean in county boundaries. I have a pre 2010 file that did this and I can't seem to find a 2010 version that follows the same pattern.   For Reference, I have been looking here: http://www2.census.gov/geo/tiger/ Below are visuals of the difference **Pre 2010 - Hugs land** ![Pre 2010 - Hugs land](http://i.stack.imgur.com/zkk0z.png) **2010 Includes Ocean** ![enter image description here](http://i.stack.imgur.com/3OJRU.png)


title: Is there a way to clip areas that are intersecting in Cartodb?

text: I am new to cartodb and wish to know how can one clip only the areas that are intersecting?


title: Great GIS, maps & cartography quotes

text: Programmers' list has been a pleasure to read. Stats' list is growing fast as well. So perhaps there is a place for separate one on this site? Please share your nuggets of wisdom / humour / inspiration from the world of GIS, cartography and mapping.


title: USA County level, detailed OpenStreetMap basemap in tiff?

text: Is there any way to get detailed (close zoom) OSM tiffs at the county level in the USA? I would like to use them as basemaps for a class I am teaching in GIS, using QGIS. I know about (and love) the OpenLayers plugin, but it can be flaky (http://hub.qgis.org/issues/5827) and it depends on a decent network connection. I don't care if the maps are a year or so old, and I also don't care if they are really big. I would like full coverage of Thurston County WA with full detail, though. I would most like to have an FTP site to download from, not some app or on- the-fly generated thing. EDIT: Here is the workflow I am using, which I think will satisfy my needs:   1. Download pbf file for my state from Geofabrik   2. Use Maperitive to generate an MBtile for the region I want (a few square miles only currently, for demonstration, will try the county later today). I can run this overnight if necessary (generat-mbtile minzoom=1 maxzoom=18)   3. Load the mbtile into QGIS as a raster format (yay GDAL!) Thanks to everyone for your help, even though I didn't take any of your ideas in total they helped me figure out the pieces. I keep thinking I want to produce an "atlas" (metaphorical) of US Census tracts, and one of the things I would include for each tract would be an mbtile and tiff, I think. It might be a good service to provide these, maybe at the Census Tract level too, for us desktop mappers... If I am missing something critical, please clue me in....


title: Using IRelationalOperator ArcObjects interface with feature classes stored in Oracle?

text: Could you give an example how you’d use the IRelationalOperator ArcObjects interface with features from different feature classes stored in oracle database? In my visual studio project I have an xml file which reference all my actual data and I have created variables in C# to hold the values of these data referenced in the xml file. ESRI's example only works if you know the actual features and their index # e.g.               // Get the marsh and highway features from their respective classes.     IFeature marshFeature = vegFeatureClass.GetFeature(518);     IFeature highwayFeature = roadsFeatureClass.GetFeature(39);          // Get the geometries of the two features.     IGeometry marshGeometry = marshFeature.Shape;     IGeometry highwayGeometry = highwayFeature.Shape;          // Cast the highway's geometry to IRelationalOperator and determine if     // it crosses the marsh's geometry.     IRelationalOperator relationalOperator = (IRelationalOperator)highwayGeometry;     Boolean crosses = relationalOperator.Crosses(marshGeometry);     Console.WriteLine("Highway crosses marsh: {0}", crosses);      I want to be able to use this interface without knowing what the exact index value of GetFeature() function is.


title: How can I shift an entire .las file in Z?

text: I have a set of .las tiles. I needed to edit one of the tiles in Trimble's Inpho DTMaster to remove some noise. The coordinates are correct in the program, showing a my desired elevation. On export, however, the program shifts the entire file up in Z by 90 feet. It is not a spatial reference/coordinate problem, I've ruled that out already (unless the software is doing something unbelievably dumb). This is all done in a local spatial rectangular system. ![enter image description here](http://i.stack.imgur.com/6hl2c.jpg) Is there an easy way to subtract 90 from the elevation field of every point in the .las? I'd like to try to offset the Z. I've tried to convert it to an ascii file to do it semi-manually, but it is too large and fails after it grows to about 2.5 GB. I have access to various software, including ArcGIS, QGIS, Inpho Photogrammetry, Global Mapper.


title: Word cloud labels in QGIS from database

text: I would like to do some word clouding as labels in QGIS from a series of tables in a database. Examples of word clouds can be found here: Word Cloud - Google Search I have two tables in a database: locations and codes. The locations are cities that have geometry. The codes are the words I would like to use as labels. I have a count with the codes, which I hope to use as text font size. Does anyone know if there is a way to label in QGIS from a related table. I have drawn a schema and example of what I would like to achieve. Any thoughts are welcomed. ![enter image description here](http://i.stack.imgur.com/4jS2z.png)


title: Is there a big (if any) performance difference between Geoserver SQL Layer and OpenLayers geojson Layer?

text: I have a table with multiline geography items in SQL Server 2008 (Routes). I want to show them on OL Google Map but I'm having trouble creating SQL Data source in geoserver (something related to TCP/IP although its opened and enabled in SQL Config). As a workaround , I've create a web service which brings GeoJSON from the SQL and add it as a layer. My question is , How "bad" performance wise is my solution if I want to show couple of routes that each is made from around 7k points ?


title: How to fix OpenLayers plugin for QGIS Master?

text: I can't get the OpenLayers plugin to work in QGIS Master (1.9), I loaded the repository successfully, but the plugin doesn't come up as an option anywhere. I'm new to this, but have tried to make sense of several fixes (here: http://hub.qgis.org/issues/6979). I don't know where/how to save the patch/file. Any help is appreciated, I have to have some kind of streets basemap by the morning. Thanks!! Brian


title: How to create sequential field based on sorted data? (or extract pop-up data?)

text:
I downloaded a kmz field and converted it to a GDB fc on arcmap. However, the field I need to label the layer by is stored in the pop-up html field... It looks like this:                         <META http-equiv="Content-Type" content="text/html">     <meta http-equiv="content-type" content="text/html; charset=UTF-8">

     Benchmark: 100

     Point ID	     100

     Type	     Vertical

     Northing	      2119938

     Easting	     533366

     Elevation (NAVD88)	     6.1 ft

     Elevation (NGVD29)	     7.18 ft
      All I need is the benchmark number (Benchmark: 100) from each field Is there a way to either extract fields from the popup field OR create a new field based on the sorted popup field? Tried a bit of the code here: How to populate Unique ID field after sorting? but did not get anywhere. Any advice would be greatly appreciated


title: What do different realizations of the same datum have in common?

text: Firstly, am I correct in assuming that the process of determining the shape and the size of the spheroid is independent from the process of determining the origin and the orientation? For example, for NAD27, it seems that the following equation -- called the Clarke 1866 ellipsoid (reference 1) -- was fitted to the surveyed points: ![enter image description here](http://i.stack.imgur.com/kEEtX.gif) where the parameters to be determined were x_0, y_0, and z_0 instead of a general equation of an oblate spheroid: ![enter image description here](http://i.stack.imgur.com/draeR.gif) Suppose we have the complete equation (with parameters already determined), and the origin and the orientation relative to some arbitrary standard coordinate system, of each different realization of all the datums. Further suppose we do not know the name of each. For example, NAD83(CORS94) may be called "1" and WGS84 (G730) perhaps "2". So they are all jumbled up. Is it possible to sort them into different groups so that each group consists of the different realizations of the same datum? * * * References:   1. ARSITECH "Constants for Reference Ellipsoids used for Datum Transformations" http://www.arsitech.com/mapping/geodetic_datum/


title: Converting SBET (Smoothed best trajectory) files to an ASCII format

text: I have LiDAR trajectory data in an SBET (Smoothed best trajectory) format, but to the analysis I wish to carry out I need it in an ASCII format. Does anyone know if there's a way to convert it? Thanks for your help :)


title: Spatial Join of 2 line layers offset by minuscule amount

text: Let's say I have 2 line layers (green & red, in the image below) and I need to transfer attributes from one layer (red) to another (green) by doing a Spatial Join. The layers are identical, but offset by a minuscule distance, b/c of some reprojection & geographic transformation -- something like 0.005 meters (5 mm). Of course this is not even a meaningful offset for any real-world analysis issue, but for a Spatial Join, it seems it would prevent the lines from actually overlapping perfectly & thus the Intersect operation would transfer the wrong attributes. Using "Within a Distance" w/ some similarly tiny distance would invariably capture not only the parallel green line that the red line was meant to be on top of, but the perpendicular green line as well. It's hard to believe I'm really out of luck on such a simple operation & have to resort to some much more complex procedure, simply due to a less than 100% absolutely perfect overlap between 2 otherwise completely identical line layers. And I can't really imagine any workaround. Any ideas?? Thanks very much! Image below!: ![enter image description here](http://i.stack.imgur.com/nSZLN.jpg)


title: Unsupported operand error writing to file (Python, arcgis)

text: In the code below, I get a "TypeError: unsupported operand type(s) for +: 'NoneType and 'str'" on the f.write line. If I don't concatenate, it runs fine, but all the layer names run together. I understand what the error means, but not why I'm getting a NoneType. I've tried casting lyr.name to str, same result. All I'm trying to do is write out a text file of layer names. Any idea what's causing this and how to get around it? I've tried testing for type is None, isinstance(None), no luck.Thanks.                   f = open(r"C:\ListOfLayers.txt", "w")         count = 0         lyrs = arcpy.mapping.ListLayers(mxd, "", df)              for lyr in lyrs:             f.write(str(s)) + "\n"             count += 1              f.write("\n \n There are " + str(count) + " layers.")              f.close()


title: Can you geocode to a polygon?

text: I have a postal code file, I want to geocode to that polygon (not point) so that I can perform some analysis on each of the polygons. Sometimes more than one point will share the same polygon so it can't be a 1 to 1 join.


title: Resources for mapping and displaying spatial data in R

text: I'd like this to be a repository for all resources available specifically for mapping and displaying spatial data in R. I've seen (and contributed) to this post, but would like this post to focus specifically on data visualization and mapping in R (as opposed to spatial analysis in general).


title: Is 64-bit background geoprocessing native to ArcGIS 10.2?

text: I cannot seem to find any documentation as to whether or not 64-bit background geoprocessing is native to ArcGIS 10.2. Or, is a separate download and install required as in 10.1? There were rumors prior to the 10.2 release that 64-bit background geoprocessing comes pre-installed with the upgrade/download. Is there an authoritative reference on this subject? No relevant system information is available on my Arc Administrator configuration information.


title: Including polygons as part of network for ArcGIS Network Analyst?

text: I am creating a pedestrian walking network of a downtown area. I am digitizing line segments for paths and sidewalks, but I need to incorporate parking lots into my analysis. So far, I have been creating multiple lines through parking lots to represent possible routes that a pedestrian might take through the lot. This is tedious, and somewhat random. I have my parking lots as a polygon feature class. Is there any way for me to include these polygons in the network by snapping my walking network lines to the parking lot polygons?


title: is there any way to run a GIS (Esri, Qgis or mapinfo) on a Windows8 RT tablet?

text: I am looking for a good idea to run a GIS program on a Dell XPS 10 with the Wind.8 RT.Thanks for any


title: Use a Python script to stack layers from multiple directories and produce multiple outputs in ERDAS

text: I am using ERDAS Imagine 2013 and am trying to create a Python script that takes lists created in another script containing the path to each scene's raw data folder and base image ID's (for example LC80230282013183LGN00). I then use a for-loop to go to each data folder, selecting the layers to stack to produce my stacked images. Also, I'm working with ETM+ and OLI data so I have a conditional statement to determine which bands to stack (1-7 for OLI and 1-5 & 7 for ETM+). What I have so far will stack the layers for the first data folder in the list, but then I get an error when it attempts to stack layers in the next folder in the list. I'm fairly inexperienced with Python and coding in general and I feel like either my logic is wrong or I'm just way off with my code. What I have:               from imagine import modeler          from getpath import imageName, imagePath, outputPathStack          m = modeler.Model()          numImages = len(imageName)          for i in range (0,numImages):              #check for landsat 8 image         if "LC8" in str(imageName[i]):             ri1 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B1.tif")             ri2 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B2.tif")             ri3 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B3.tif")             ri4 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B4.tif")             ri5 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B5.tif")             ri6 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B6.tif")             ri7 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B7.tif")                  StackBands = m.StackLayers(ri1, ri2, ri3, ri4, ri5, ri6, ri7)             ro = m.RasterOutput(StackBands, outputPathStack + \                                 imageName[i] + "_stack.img")         #check for landsat 7 image            elif "LE7" in str(imageName[i]):             ri1 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B1.tif")             ri2 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B2.tif")             ri3 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B3.tif")             ri4 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B4.tif")             ri5 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B5.tif")             ri6 = m.RasterInput(imagePath[i] + "/" + imageName[i] + "_B7.tif")                  StackBands = m.StackLayers(ri1, ri2, ri3, ri4, ri5, ri6)             ro = m.RasterOutput(StackBands, outputPathStack + \                                 imageName[i] + "_stack.img")         m.Execute()      The error I get is: File "C:\Workspace\stacklayers.py", line 35, in m.Execute() RuntimeError: erdas::sblib::Operator::SetErrorMessage failed Spatial Model failed in Raster Output 2. Error: erdas::sbsupportlib::LayerOperator::OnExecute failed erdas::sbsupportlib::LayerOpOnExecute::WriteBlock failed Bad response from exchange.. Thank you in advance for any help or advice


title: Selecting features to project to different coordinate system?

text: I've got data from all over Europe in one feature class projected in Web Mercator Auxiliary Sphere. Due to the fact, that length and area values are not usable I'm now working on a python script (ArcGIS 10.0) that projects the features to the corresponding UTM zones and writes the calculated values back to the orginal fc. First the script assigns the corresponding UTM zone to each feature via spatial join with an UTM zone fc. Within the next step I'd like to select all features within one UTM zone and project them (e.g. all features with value 32 will be projected to UTM ETRS89 zone 32 and so on). My problem is, that project_management does not have a "where_clause" so what will be the best way to make this selection? I've got lot's of data in one feature class and many feature classes, so speed does matter... Does anyone know a better way to get precise length and area values for data in Web Mercator projection?


title: SDK GXP Legend OpenGeo Suite Scroll?

text: Is it possible to have the legend scroll using the OpenGeo Suite SDK and GXP. Following this tutorial Tutorial LinkIt is possible to have it any size if external as the height and width can be changed but when embedded into the west container, if too many layers are selected the legend doesn't fully display each layer? Suggestions or ideas?


title: Counting and sum by two layers

text: I´ve got to layers, both polygons. Layer 1 consists of propertiers, where ID is individually names. Layer 2 consists of lakes. The query I want is: Total count of squaremeters of lakes for each property in layer 1. And an output with the results for each output. Maybe as add to the already existing table in layer 1? Tablestructures with area-info are ready go. I just need s tool and guide how to do it... Thank you


title: How to assign IDs to villages which may have duplicate names?

text: I have more than 150 village points, each with a village name attribute `village`. The same village name may appear on multiple points. I want to assign individual IDs for each village name, say `vill_A = 1`, `vill_B = 2`, and so on. I tried to calculate `[village]`, but it didn't work. Any suggestions?


title: How do I know if I am running 64-bit background geoprocessing on ArcMap 10.2?

text: I'm running ArcMap 10.2 and just set up 'ArcGIS for Desktop Background Geoprocessing (64 bit)' and now I have both the ArcGIS10.2 folder and ArcGISx6410.2 folder inside my python folder. I have background processing enabled as well. When I open the python window inside of ArcMap and do import sys, print sys.version it shows that it is running the 32 bit version. How do I get the 64 bit version to run inside of ArcMap?


title: How can I fix incorrect links in my QGIS install?

text: Today I installed QGIS 1.8.0 and I have the same error as this one: How can I fix QGIS 1.8 startup issues due to missing qwt5.dll library?. The sollutions how ever don't seem to work. What I have found out about it is that certain files and subfolders are in folders where they are not supposed to be. For example: the QGIS-executable is in c:\program files\quantum gis lisboa\apps\qgis\bin\ while it should be in c:\program files\quantum gis lisboa\bin\ where also the qwt5.dll-file can be found. By copying the files from the first folder to the latter I am able to start QGIS. I still get errors though, for example that plugins can't be found. This again, I can solve be copying folders. Now I can use QGIS more or less, but it doesn't run like it is supposed to. Does anyone know how I can fix this problem?


title: Find polygon centroid with Python

text: I'm using ArcGIS 9.3 to find centroid of some polygon features. I'd like to do this in Python, but so far I can only find solutions with VBA. Can anyone provide some leads here?


title: Interactive maps for a survey - Are there any tools yet?

text: i am working on a draft for an online-based survey. I would like to enable the participants of the survey to draw polygons, lines and points (several) on a basemap (e.g. google maps) and (extra feature) additionally to comment on these elements. Each element drawn by one participant should be stored. I am looking for a ready-to-use method/tool to embed it in the survey. Right now, a fully costumized solution is not an option. Unfortunately, i couldn't find a suitable tool yet. The whole idea is more or less a mental map, an interactive tool to highlight topographic elements on a basemap. If you know some possible solutions or tools to meet my requirements, i would appreciate any suggestion. Thanks, J


title: RARE II Roadless Data

text: I am looking for data showing RARE II Roadless Areas on USFS Lands. RARE II was an early inventory of Roadless areas that was completed in the late 70s. This dataset is different than areas that are currently recognized as


title: Query Layer to access data on Sharepoint

text: I want to connect to a Sharepoint list from ArcMap. Without knowing very much about these things I thought that it would be possible to read this non- spatial data to create lookup tables in my file geodatabase. In searching for an answer, I can't seem to get away from ArcGIS for Sharepoint, which isn't what I'm looking for at all. I'm using ArcGIS 10.0


title: Feature Class To Geodatabase => ERROR 000210: Cannot create output

text: error messages:               Executing: FeatureClassToGeodatabase E:\jags20120321.mdb\ch_cs "Database Connections\Connection to localhost.sde" "Database Connections\Connection to localhost.sde"     Start Time: Sat Jun 16 13:32:28 2012     Running script FeatureClassToGeodatabase...     Failed to convert E:\jags20120321.mdb\ch_cs. ERROR 000210: Cannot create output Database Connections\Connection to localhost.sde\jags.DBO.ch_cs     Failed to execute (CopyFeatures).          Completed script FeatureClassToGeodatabase...     Executed (FeatureClassToGeodatabase) successfully.     End Time: Sat Jun 16 13:32:37 2012 (Elapsed Time: 9.00 seconds)      Those are the error messages. I don't know what's the reason. I use arcgis 9.3, arcgis sde 9.3, and sql server 2008.


title: Is ArcGIS 10 compatible with Windows Server 2012 OS?

text: Is ArcGIS 10 for Desktop, ArcSDE 10 for Windows, and/or ArcGIS Server 10 for Enterprise compatible with Windows Server 2012 OS? Or do we have to install 10.1 on Windows Server 2012?


title: Running Random Forest Classification on WorldViev-2 image in R - problem

text: I have a problem with performing RF in R for a WorldView-2 Scene. I've never ever before used R for remote sensing classification, so I simply followed what is written here. I've loaded scene (8-band TIFF) with raster::brick and called it abcd, loaded simple 4-point, 2-class shapefile with OGR, run randomForest model and called it wv2rf and tried to use predict:               predict(abcd, wv2rf, "rftes.img", index=1, na.rm=TRUE)      R returned an error, that I cannot handle with:               Error in `[.data.frame`(blockvals, , f[j]) : undefined columns selected      I guessed that it has to do something with index value, but I tried to change it and nothing worked whole code:               #loading data     abcd <- brick("13AUG20094646-M2AS-13EUSI-1283-01.tif")     sdata <- readOGR(dsn="mypath\\poligony testowe", layer="training")     v <- extract(x=abcd, y=sdata, df=TRUE)     sdata@data = data.frame(sdata@data, v[match(rownames(sdata@data), rownames(v)),]     sdata@data[3] <- NULL     #classification     wv2rf <- randomForest(x=sdata@data, y=(sdata@data[,"klasa"]), ntree=10, importance=TRUE)     predict(abcd, wv2rf, "rftes.img", index=1, na.rm=TRUE)      Here are files with my sdata, wv2rf and abcd variables, should work with load("filename") in R. Hope they help you help me. http://speedy.sh/9f7cf/rf.zip


title: How to execute gp.execute ArcSDE tool in C#

text: So I created an ArcGIS 10 toolbox in ArcSDE and added a tool to it. There are enough examples out there on how to execute a tool.tbx that resides outside SDE or GeoDB. But in C#, how do you use gp.execute on a tool that resides **within ArcSDE**? thanks!


title: Open Layers, Geoserver WMS and OSM Cloudmade?

text:
While I have some familiarity with GIS, this is my first attempt using OpenLayers. I am trying to display the following layers in Openlayers: 1) Cloudmade Base Map 2)WMS Data overlay from my Geoserver (`http://54.244.249.54/geoserver/PalmBeach/wms?service=WMS&version=1.1.0&request=GetMap&layers=Florida2&styles=&bbox=-82.97059340790229,24.320616400434016,-80.1393763936513,26.907996410185536&width=512&height=467&srs=EPSG:4326&format=application/openlayers`) I have used the Cloudmade sample as a starting point (http://developers.cloudmade.com/wiki/21/CloudMade_Tiles), and have successfully changed the parameters to center and zoom the map, but I cannot get the WMS to display as a layer, despite it being listed as a layer in the layer selector. Im not certain if its a projection issue (I changed the projection in Geoserver to match the Open Layers projection EPSG:4326) or just my infamiliarity with Openlayers, but would appreciate any guidance. My code is below:                   <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN"  "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">                   <title>Using CloudMade Tiles with OpenLayers</title>         <link rel="stylesheet" href="openlayers/theme/default/style.css" type="text/css" />

                                         Thank you for any insight and suggestions!


title: ArcMap crashes after using "geocode addresses" tool

text: I'm working with ArcGIS 10.2 in a project where I need to publish my own geoprocessing service in order to Geocode a table of addresses. I tried creating models in model builder then publishing them and everything was working fine. But after trying to publish a model wher the input is an Excel file that I convert to a tableau and then geocode it, ArcMap crashed and since then everytime I try to run a model containing "Geocode addresses tool" ArcMap crashes. I tried to run a simple "Geocode addresses" from the toolbox and the toolbox crashes. I tried deleting temp files and changing the "normal.mxt" and rebooting and the problem persists. I find this behavior strange since the same models were working before even in ArcGIS Server and now everytime a model is using "Geocode addresses" ArcMap crashes in the "geocoding addresses" step.


title: How to find Latitudes and Longitudes of projections of the vertices of a rectangular plane below earth's surface?

text: I want to find out the latitudes and longitudes of projections of the vertices of a rectangular plane inside the earth's surface. I know dimensions of rectangle, angles of orientation and latitude and longitude of the projection of the center of the rectangle. Earth's curvature is also to be considered I added an Image which is a 2D representation of what i want. The plane actually is having 2 angles (its trace makes an angle 'A' with North and the plane makes an angle 'B' with the horizontal). I want to find lat longs of all 4 projections Note: The dimensions shown are just for depiction. they are not to scale ![enter image description here](http://i.stack.imgur.com/uiiMH.png)


title: How to create multi-image single band single file geotiff?

text: I have a geotiff file that is single band but contains multiple subdatasets. I use gdal with numpy to open a specific subdataset and mask some of its values. Example:               $ gdalinfo sample.tif     ...     Subdatasets:     SUBDATASET_1_NAME=GTIFF_DIR:1:sample.tif     SUBDATASET_1_DESC=...     SUBDATASET_2_NAME=GTIFF_DIR:2:sample.tif     SUBDATASET_2_DESC=...          # pseudocode     sds = gdal.Open('GTIFF_DIR:2:sample.tif')     sds_mask = apply mask to sds      I want to put sds_mask back into the original geotiff or a copy of it. Is this possible with gdal? I tried writing out all the subdatasets, modifying one of the subdatasets, and then merging them back together with gdal_translate and gdal_merge.py with the -separate option but it does not create a geotiff with subdatasets.


title: Missing hstore column in TileMill!

text: I made a PostgreSQL with PosGIS Extension DB with the pgsnapshot_schema and loaded an .osm.pbf-file with OSMOSIS into it. Now i connected it to TileMill, but if i look up into the table i just see columns with the datatype of integer(id, version, User_id,changset_id). But i would like to use Tilemill with the hstore column to select the data and use it for the labels! In TileMill i wrote into the **Table or subquery:`ways`**; **Geometry field:`linestring`** Does anyone know, how i can fix it easily? ![This is how it looks like](http://i.stack.imgur.com/JFH1D.png)


title: How to use "this" Levenshtein Python code into ArcGIS 10

text: I need some help with this: First of all, I wanted to "Compare strings from two fields", and asking before at gis.stackexchange, in this post someone help me with this topic, referring me to a Python code, but, unfortunately, from the deepest ignorance of this kind of issues, I wanted to know if it is possible to run this code within ArcGIS 10 (field calculator). The code: PyLevenshtein , Thanks in advance, and sorry if I wasted your time.


title: Finding nearest raster cell value based on vector point?

text: I have two layers in ArcGIS: One is a vector point and the other is a raster layer. I want to append the value of the raster data grid to the vector point. The problem is the vector point is away from any existing raster grids with value, so the command "Sample" return no value for those vector point. So my question is: How can I find the nearest raster cell given a vector point and extract the value from the cell? Actually there are two things I need to do:   1. Identify the outliner within a ranger (Usually 1-2 cells away from my existing raster)   2. Assign a value from the nearest raster to them, if they are in a range * * * Edited: I have about 3000 vector points and my problem is this: ![enter image description here](http://i.stack.imgur.com/jsYWZ.png) I need to append the raster value based on locations to the vector points. I used the tool "Sample" and it works well for most of the points. However, there is a situation like this ![enter image description here](http://i.stack.imgur.com/tE3OP.png) The right point can get the value from the raster ("Sample" works) but the left point cannot due to the alignment issues. Raster to polygon is not working because I get error message that I cannot turn them into vector because they are out of the domain. Also I have points like this![enter image description here](http://i.stack.imgur.com/7uwMQ.png) Those upper points should not get any values because they are too far away from the raster. I thought of someways to work around: Step 1. Do the "sample" first Step 2. Select the null value after the "sample" operations Step 3. construct buffer zones (within 1 decimal degree) based on the null points Step 4. use **zonal statistics??**? or other tools to loop through all the raster cells in the buffer, find the nearest one, extract its value and put it in the vector point. I am stucked at step 4. I don't know which tool in ArcSDK I can use to achieve this function. Or else...Do you have better idea to deal with that? I have 20 raster layers like this and I would like to make it automatically (using model builder and arcobject).


title: Opening many shape at the same time

text: It may seem like a simple question but I haven't found a solution. I have a structure of directories with many shapefiles inside. Is there a method to open all these shapefiles at the same time with ArcMap? The function "add data" forces me to do it for every directory separately.


title: Any restrictions on the use of ArcGIS REST API?

text: I'm building a PHP script which ends up with an array of latitude, longitude pairs (around 2,500), and I want to be able to reverse geocode them into full adresses, then recording them on my database. I ran a search on Reverse Geocoding and have narrowed my options to Google and ESRI. I'm leaning towards ESRI because they have a REST API, whereas Google works via Javascript, which would give me much more trouble to implement. But before I start coding I was curious whether the ArcGIS REST API has any restrictions at all: is it free or do I need to subscribe? Do I need an API key? If so, where do I get it? Are there restrictions regarding requests per minute / hour / month? I've looked around their docs, but can't seem to find anything non-technical.


title: By converting into a shp file points are reduced

text: I am new here, so sorry if I missed something, but when I try to converta data from an xls file into a shp file by export data, all the points I had before get reduced to only one?! I also want to convert Eastings and Northings into Lat/Lon and have unfortunately two coordinate systems that I need to align thus GDA1994 and WGS1984. I am also having the problem that the points and a coastline dont show up, but separately they do. Thanks a lot already. A little more info :) Sorry for that ... I am using ArcGIS 10.1 and what I did was this so far: in ArCatalog export to a gdb (single)+pulled it into ArcMap-> Display X,Y and then I saw nicely all the points that were originally in the table (a few thousand bird positions), but then if I export the data as a shapefile, only one point is left instead of the few thousand positions?! This would be my first question, What am I doing wrong here? The coordinate system should be WGS84, but the that aim would to convert the Eastings and Northings of all these points (as I only have this data) afterwards into Longitude and Latitude, which would be my second question (tried to do it via "add fields"+calculate geometry)...The underlying map for all these points has got the coordinate system GDA94. I hope this helps a little more :) Thanks again, highly appreciated!!! PS. Could it be due to different different sings and formatting in the excel table? When do I use a csv file?


title: Ortho Projection produces artifacts

text: I`m trying to create a sphere-like view using qgis and the "world from space"-projection http://spatialreference.org/ref/sr-org/6980/ (essentialy an ortho-projection). ArcGIS wraps the shapes correctly but QGIS (2.01) produces nasty artifacts. ![enter image description here](http://i.stack.imgur.com/1VNJA.png) I have to produce the globes on a regular basis with different angles so does anyone out there have an idea how to fix this problem?


title: Trouble with opening/using i.atcorr in GRASS

text: I'm relatively new to GRASS. I'm running version 6.4.3 on a Windows 7 machine. I'm trying to process Landsat 5 imagery from multiple states, so I was hoping to approximately follow the workflow here: Is there a way to automate the preprocessing of LANDSAT data? I have no problem running i.landsat.toar, but when I try to use i.atcorr I get errors: When I start to type into the command console, usually I type the command, then a space, the enter and it opens the GUI window. When I type i.atcorr, as soon as I press the space, this displays: > Traceback (most recent call last): File "C:\Program Files\GRASS GIS > 6.4.3\etc\wxpython\gui_core\prompt.py", line 1133, in OnKeyPressed > > self.cmdDesc = gtask.parse_interface(GetRealCmd(cmd)) File "C:\Program > Files\GRASS GIS 6.4.3\etc\python\grass\script\task.py", line 470, in > parse_interface > > tree = etree.fromstring(get_interface_description(name)) File "C:\Program > Files\GRASS GIS 6.4.3\etc\python\grass\script\task.py", line 448, in > get_interface_description > > "\n\nDetails: %(det)s") % { 'cmd' : cmd, 'det' : decode(cmderr) } > grass.script.core . ScriptError : Unable to fetch interface description for > command 'i.atcorr'. Details: [Note that no details appear after it says "Details:"] Then if I press "Enter", an error pops up with the same text: > "Unable to fetch interface description for command 'i.atcorr'. Details: " I thought perhaps the problem was that I didn't actually have the function installed on my computer, so I went through g.extension and get this: > g.extension extension=i.atcorr svnurl=http://svn.osgeo.org/grass/grass- > addons/grass6 which: wget: unknown command WARNING: The GRASS program path: > [C:/Program Files/GRASS GIS 6.4.3] contains a space. The >build which > follows will likely fail. Fetching from GRASS-Addons SVN (be patient)... > svn: E170000: URL 'http://svn.osgeo.org/grass/grass- > addons/grass6/imagery/i.atcorr' doesn't exist ERROR: GRASS Add-on i.atcorr > not found in repository or no network connection or another >problem" Google tells me that I'm not alone in this "contains a space" error, but I couldn't find any clear information on how to solve this or work around it. Any insight or advice would be much appreciated. I am not sure if I can circumvent the problem by just using the command console, but as I looked through the wiki on the topic (http://grasswiki.osgeo.org/wiki/Atmospheric_correction) it wasn't clear to me if I was supposed to build my own icdn text file? That's not something that came when I downloaded the data from USGS. I was hoping the GUI version of the function would make this more clear to me but since I can't get it to open, I'm still confused. Thanks in advance.


title: PyQGis - Working with Spatialite

text: I wrote a Vector layer to a spatialite file/db using the dedicated API tool, using valid driver names supported formats by OGR. When connecting to that db, I can't retrieve properly the geometry objects when using SQL command over the created db (other regular SELECT/INSERT/... SQL commands are perfectly working) Creating the Spatialite File using API Tool               QgsVectorFileWriter.writeAsVectorFormat(layer, layer.name() + ".sqlite",                                                             "utf-8", None, "SQLite")      I noticed it gives a default "GEOMETRY" name to the geometry column within my table Then, when trying to execute some SQL comman over that file:               from pyspatialite import dbapi2 as db          db_filename = 'db.sqlite'     conn = db.connect(db_filename)     c = conn.cursor()     ex = c.execute          rs = c.execute('SELECT sqlite_version(), spatialite_version()')     for row in rs:         msg = "> SQLite v%s Spatialite v%s" % (row[0], row[1])         print msg          rs = c.execute('SELECT IsValid(GEOMETRY) FROM table LIMIT 10' )     for row in rs:         msg = row[0]         print msg      returns: SQLite v3.7.17 Spatialite v3.1.0-RC2 -1 for every selected row, which mean that the geometry is invalid. Other function such as AStext(GEOMETRY), etc return None or various error code Of course, if I reload that file wthin Qgis using the API command, geometry is correct. I suspect I misunderstand what QgsVectorFileWriter is actually doing, please advice. (note: I'm using GDAL 1.9)


title: PostGIS Column has Z dimension but geometry does not

text: I am able to create a geomtry column for 3d multiline using this:                   CREATE TABLE trajectories     (       gid serial NOT NULL,       vehid character varying(10),       state integer,       geom geometry(MultiLineStringZ),       CONSTRAINT trajectories_pkey PRIMARY KEY (gid)     )     WITH (       OIDS=FALSE     );     ALTER TABLE trajectories       OWNER TO postgres;      But when I tried to alter an existing table by modifying an existing geom column or adding a new one, I am getting he following error: > > ERROR: Column has Z dimension but geometry does not > > ********** Error ********** > > ERROR: Column has Z dimension but geometry does not SQL state: 22023 ERROR: > Column has Z dimension but geometry does not               ALTER TABLE cl ALTER COLUMN geom TYPE geometry(MultiLineStringZ)        1. How do I modify and existing geometry column to geometry(MultiLineStringZ)?   2. How do I add new geometry(MultiLineStringZ) column to an existing table?


title: Using newly created features in the next part of a python script

text: Good Afternoon I am trying to write a python script which buffers a large number of features from a geodatabase (but not all) and then dissolves the newly created buffer shapefiles together. Currently I have only been able to achieve this by listing the specific locations of the newly created buffer shapefiles (which is OK for a couple of features but a nightmare with more). I am new to python and unsure how to references too (or ideally make a list) of the newly produced buffer shapefiles without individually listing there locations. This is also further complicated by producing 3 different size buffers for each feature class. Here is the first bit of my script before I carry out the dissolve. Many thanks in advance for any pointers.               # imports arcpy and os     import arcpy, os     from arcpy import env          # sets the workspace to my gdb     arcpy.env.workspace = "C:/GIS Home/Mining Features (MATTHEW).gdb"          # Output workspace location is defined     output_workspace = "C:/GIS Home/project_1"          # list of features to buffer in gdb     fcList = ["Shaft_Metalliferous", "Quarry_Metalliferous"]          for featureClass in fcList:          arcpy.Buffer_analysis(featureClass, output_workspace + "/" + featureClass[:-4] + "_b250.shp", 250, "", "", "ALL")     arcpy.Buffer_analysis(featureClass, output_workspace + "/" + featureClass[:-4] + "_b500.shp", 500, "", "", "ALL")     arcpy.Buffer_analysis(featureClass, output_workspace + "/" + featureClass[:-4] + "_b750.shp", 750, "", "", "ALL")


title: How can I accurately measure area of polygons using field calculator in qgis?

text: I have several polygon shapefiles, both set to WGS84 CRS. I have been using the google basemap with the Openlayers plugin, using the mercator CRS. On the fly CRS transformation is enabled. Should my shapefiles be set to the same thing as the google basemap? When I use the measure area tool (View>Measure>measure area), the areas of the polygons in my two shapefiles are accurate. However, when I try to calculate the area using the Area function of the field calculator, the values are incorrect. I have read posts saying that the openlayers plugin distorts the areas. I tried saving the shapefiles as wgs84, opening them in a new project, and calculating the areas, but the values were the same. I am in Pennsylvania, USA. If not the mercator CRS, what should I use? Should find an alternative basemap?


title: spatial data filtering without GIS?

text: I have access to a database here (description of the columns are also available) containing spatial data in text format and with WGS84-EPSG:4326 coordinates. I would like to only select data in my country (Czech Republic) and save them to our database. Is there any easy way to do this without need to import it to spatial DB and then filter it out? I am thinking about very simple calculation using simple conditions etc. Is there a better solution than to draw a rectangle in QGIS covering whole country and then use it to define range for longitude and latitude and use it as condition for filtering? Thanks


title: Merging multiple .tab files?

text: I was wondering if there was a quick way to merge multiple .tab files together? I have about 100. Would appending one tab to another with the MapBasic window open, then changing the name of the table appended to TableInfo(1,1) work? Thanks, Alan


title: Many problems with QGIS processing

text: I'm having a lot of problems with QGIS processor. Most of the algorithms show me this QGIS log error: > La fuente de datos no es válida. (Data source not valid) This happens with SAGA and GRASS algorithms mostly. This is pretty annoying because I can't work properly. The log processor only tells me this: > > processing.runalg("grass:v.overlay","/home/luismi/Documentos/Cartografia/arqueologia/Mapa_IMperio_Romano/vectoriales/2/roman_empire_ad_200_provinces.shp",1,"/home/luismi/Documentos/Cartografia/arqueologia/Mapa_IMperio_Romano/pruebas/vectoriales_prov/Imperio_Extension.shp",1,False,"-9.48732,43.12373,22.74626,55.10117",-1,0.0001,0,None)** What should I do? Where I should look? I'm using QGIS 2.4 on Lubuntu 12.04, processor plugin 2.2.0-2. Ask for more info if it's needed. EDIT: I've installed OSGEO Live-7.0, it happens the same thing. When I run "Spatial point Pattern analysis" this is what I get:               Algorithm Spatial point pattern analysis starting...     geostatistics_points "Spatial Point Pattern Analysis" -POINTS "/tmp/processing/ed1a459a603a4f20927861d30a5c06b2/ne10mpopulatedplacessimple.shp" -STEP 5 -CENTRE "/tmp/processing/7c03f63d9e4a48418f4cb1dbadb2afac/CENTRE.shp" -STDDIST "/tmp/processing/f1ee992e4f0040449e2523621fb31c13/STDDIST.shp" -BBOX "/tmp/processing/8e4efd5dd12546a98417de775a2d4b56/BBOX.shp"     Error: library     _____________________________________________          ##### ## ##### ##     ### ### ## ###     ### # ## ## #### # ##     ### ##### ## # #####     ##### # ## ##### # ##     _____________________________________________               68 loaded tool libraries (678 tools):     - garden_fractals     - pj_proj4     - ta_slope_stability     - tin_viewer     - garden_games     - io_table     - io_gps     - grid_calculus     - sim_erosion     - db_odbc     - garden_webservices     - grid_spline     - shapes_tools     - table_tools     - grid_gridding     - shapes_polygons     - imagery_photogrammetry     - ta_profiles     - io_esri_e00     - ta_preprocessor     - io_grid_grib2     - sim_ecosystems_hugget     - sim_hydrology     - docs_pdf     - ta_morphometry     - imagery_classification     - db_pgsql     - statistics_kriging     - io_virtual     - pointcloud_tools     - statistics_points     - sim_cellular_automata     - contrib_perego     - imagery_svm     - io_shapes_las     - grid_tools     - io_grid     - shapes_lines     - imagery_rga     - sim_fire_spreading     - statistics_regression     - docs_html     - grid_calculus_bsl     - shapes_points     - table_calculus     - io_grid_image     - io_shapes_dxf     - shapes_grid     - ta_compound     - io_gdal     - grid_visualisation     - ta_lighting     - imagery_segmentation     - tin_tools     - ta_channels     - sim_ihacres     - grid_filter     - garden_learn_to_program     - garden_3d_viewer     - ta_hydrology     - statistics_grid     - io_shapes     - shapes_transect     - grid_analysis     - climate_tools     - pj_georeference     - pointcloud_viewer     - imagery_tools          type -h or --help for further information          Converting outputs     Loading resulting layers     Algorithm Spatial point pattern analysis finished          processing.runalg("saga:polygoncentroids","/usr/share/qgis/python/plugins/processing/tests/data/polygons.shp",False,None)      And QGIS log says this again:               La fuente de datos no es válida ()              Error 4 de OGR [3]: Unable to open /usr/share/qgis/python/plugins/processing/tests/data/polygons.shp or /usr/share/qgis/python/plugins/processing/tests/data/polygons.SHP.      Second EDIT: I have installed a new Xubuntu 14.04 in another partition in my PC, I've installed SAGA, GRASS and QGIS again since ubuntugis repositories, this problem persist!!! Anyone has any idea why is happening this? Some algorithms are working fine other don't!! Third Update: ![Look at he error 4 in the log. This happen when I run only the dialog.](http://i.stack.imgur.com/0t27y.png) ![And here is the error. Look at the log part.](http://i.stack.imgur.com/QMSMA.png) File of random points created with QGIS algorithm.


title: Truncate operation not working in ModelBuilder?

text: ![enter image description here](http://i.stack.imgur.com/FWApm.png) I'm having trouble with a model that is supposed to truncate data from the attribute table of a feature class (which is in a file geodatabase not ArcSDE) and then append it with new data. The truncate operation returns a success result but the geometry remains intact after the truncate. I also tried the delete features and delete rows tools but they each return a generic function error. I'm not sure if this problem is due to the fact that the feature class being truncated is associated with a published map service. I should note that this custom geoprocessing tool is being published to a server from where it will be invoked programatically in a flex arcgis web application. What worries me is whether truncating data directly from a feature class which is associated with a published map service is a bad idea. Also, this model worked fine initially until the first time that it was invoked from the flex application through the published geoprocessing tool. I am using ArcGis 10.1, ArcGIS Viewer for Flex 3.6.


title: How to configure tilecache on nginx

text: I am running a **geodjango** application with **mapnik** and **tilecache**. Currently tilecache is run on another server to decrease the load on the application server. The application server has **nginx** as the webserver while tilecache is run on the **apache** server. Now we have a requirement that necessitates us to put both the map and application parts in the same server. Could you please tell me how to configure **tilecache url** in nginx server. Thanks in advance.


title: What are the differences between the photometric options for GDAL_merge?

text: I am using GDAL_merge to merge several geotiff tiles into one large tile. One of the creation options I can add is photometric. There are several options for this.   * MINISBLACK   * MINISWHITE   * RGB   * CMYK   * YCBCR   * CIELAB   * ICCLAB   * ITULAB What are the differences between these different outputs? I understand the MINISBLACK and MINISWHITE options, but the others are less clear. Thanks in advance.


title: QGIS - Continuous buffer image [Heat Map]

text: Is it possible to have the representation like the one the heat-map like the figure below using QGIS? The way seems to use buffer options, but I do not know how. Does exist any plugin that helps? I am using QGIS1.8 on windows 7 ![power](http://i.stack.imgur.com/RE59C.jpg)


title: Changing basemaps with different spatial references (ArcGIS Javascript)

text: I've got a new requirement to switch basemaps that have different spatial references. The app is using ArcGIS Javascript, and switching to Silverlight is not an option at this point in development. I saw this question from 2 years ago and was hoping there was an alternative. Is there an easier way to change basemaps that have different spatial references in 3.6?


title: SVG markers with PyQGIS

text: How do you SVG-marker layers in PyQGIS? I'm importing a vector layer and want to display it with SVG-markers. It seems that QgsSvgMarkerSymbolLayerV2 does not exist in PyQGIS...


title: OpenLayers take user input of XY and create a point

text: I'm looking to take input from a form (X and Y WGS84 coordinate) and then draw a point on my map when the user clicks execute. I could not locate something like this in the OpenLayers examples.


title: What information is needed to determine the section of the earth's surface visible in a photograph?

text: What information (either provided by a user or taken from sensor readings) would be required to determine the geographic area encompassed by a photograph? It seems like the area visible from a photograph would be a rectangle projected onto the surface of the Earth, but I'm unsure what information you'd need to figure out exactly what area a picture encompasses. It seems clear that you would at least need the location, height, and orientation of the camera (but maybe I'm wrong about that?). There's probably also some information needed about the focal length, but I'm rapidly getting out of my depth here (pun regretfully intended). As an alternative, suppose that you could get the user to geotag `n` points in the picture -- how many points would need to be tagged to compute the geographic bounds of the image? The use case is that I would like to be able to present various statistics to a user about the area depicted in a photograph; ideally the user wouldn't have to do too much (e.g., if this were done from within a smartphone maybe all necessary information would be available from the phone's sensors), but if that's not possible it's okay to ask the user for necessary information. Apologies in advance if this is too vague or posted in the wrong area; I'm happy to move it or add clarifying details as necessary.


title: Jsonp call to geoserver

text: If I try, like this way:                       $.ajax({            jsonp: false,            jsonpCallback: 'getJson',            type: 'GET',            url: "http://79.141.150.167/geoserver/xxx/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=xxx:business_service&maxFeatures=50&outputFormat=text/javascript&format_options=callback:getJson",            async: false,            dataType: 'jsonp',            success: function(data){               console.log(data);            }         });      My browser tells me following:               Resource interpreted as Script but transferred with MIME type text/xml: "http://79.141.150.167/geoserver/xxx/ows?service=WFS&version=1.0.0&requ…tputFormat=text/javascript&format_options=callback:getJson&_=1393412236410".      Uncaught SyntaxError: Unexpected token < (ows:1) ...what I'm I doing wrong? Thanks in advance!


title: bathymetry for a UK location

text: I'm wondering whether I can get water depth for the mudflats at Lindisfarne, NE England, UK: http://maps.google.com/?ie=UTF8&hq=&ll=55.65,-1.821&z=13 The mudflats are intertidal, so I'm hoping to get water depth according to state of tide at any particular location on the mudflats.


title: Why can ArcGISDynamicMapServiceLayer read a cache service?

text: I received a cache folder created from a MXD that I don't possess. Then I was asked to publish that cache as a map service. To do that, I created a map service using a MXD (I created) that contains one layer. That layer has no relation with the cache. Then I configured the map service to use the cache I was given. In a !flex application, when I use the map service using esri:ArcGISDynamicMapServiceLayer, it works : I get on the map what the cache contains, not what the MXD contains. I'm a little confused. ArcGISDynamicMapServiceLayer is not supposed to read from the cache of the service! Only ArcGISTiledMapServiceLayer is supposed to be able to read from the cache. Am I wrong? Why do I get the cached data on the map instead of the MXD data (which is different from the cache)? Here is a simple graphic explaining the situation. ![enter image description here](http://i.stack.imgur.com/AkGxv.jpg) Regards.


title: Where can I get Beijing GIS data?

text: I was meandering so many places to find GIS data for Beijing. but it seems like I cannot get to any of useful resources. At least I got road, water shape files and outline of beijing city. specifically **what i need is buildings and parcel(lots) maps** to investigate site analysis Beijing CBD area. Thank you in advance for your kind consideration


title: Map Algebra in VB.NET Add-In Tool

text: I'm trying to run a map algebra expression in a vb.net add-in, and can't seem to get it to work. It's a roving window that interpolates noData areas in a DEM. The expression works in the Python window, and in a Python stand alone script (after replacing the bound raster [inputDEM], with "inputDEM"). In vb it throws a syntax error. Any thoughts? Thanks so much, DR               Try                     Dim demPath As String = demLayer.FilePath.Remove((demLayer.FilePath.Length - demLayer.Name.Length), demLayer.Name.Length)                     Dim inputDem As IRasterDataset = OpenRasterDataset(demPath, demLayer.Name)                     Dim mapAlgebra As IMapAlgebraOp = New RasterMapAlgebraOpClass()                     Dim environment As IRasterAnalysisEnvironment = mapAlgebra                     Dim workspaceFactory As IWorkspaceFactory = New RasterWorkspaceFactoryClass()                     environment.OutWorkspace = workspace                     mapAlgebra.BindRaster(inputDem, "inputDEM")                     Dim outputDem As IRaster = mapAlgebra.Execute("Con(IsNull([inputDEM]),FocalStatistics([inputDEM],NbrRectangle(16,16),'MEAN'),[inputDEM])")                     Dim saveDem As ISaveAs2                     saveDem = outputDem                     saveDem.SaveAs("fillGaps", workspace, "GRID")                 Catch ex As Exception                     MsgBox(ex.Message, MsgBoxStyle.Exclamation, "Error")                 End Try