ArcObjects Runtime Binding 10.0

I installed the ArcOjbects SDK (if I remember correctly) and in visual studio I have referenced the needed assemblies. Namely ESRI.ArcGIS.Geoprocessor. When in stepping through in debug mode I hit this routine:               void BindingArcGISRuntime(object sender, EventArgs e)             {                      //ProductCode[] supportedRuntimes = new ProductCode[] { ProductCode.Desktop };                 //foreach (ProductCode c in supportedRuntimes)                 //{                 //    if (RuntimeManager.Bind(c))                 //        return;                 //}                           if (!RuntimeManager.Bind(ProductCode.Desktop))                 {                     // Failed to bind, announce and force exit                     Console.WriteLine("Invalid ArcGIS runtime binding. Application will shut down.");                     System.Environment.Exit(0);                 }             }      But I don't hit the "Console.WriteLine("Invalid ArcGIS runtime binding. Application will shut down.");" piece Then straight to this auto-generated routine:               if (RuntimeManager.ActiveRuntime == null)                 {                     EventHandler temp = ResolveBindingEvent;                     if (temp != null)                     {                         temp(this, new EventArgs());                     }                 }                      return Initialize();      Any ideas? Thanks.


How to calculate stream order from vector data in ArcGIS 10

I know that NVision Solutions has a tool for calculating stream order from vector hydrography data for ArcGIS 9.3.1. Is there something similar that works with ArcGIS 10?


Symbolizing Layers from a PostGIS Database to many Client Computers

I was hoping someone could help with a symbolizing issue. I have a PostGIS database on a server with a number of layers. The layers of a series of polygons, lines and points that all share the same symbology. I want to symbolize the layers with the same symbol schema on computers in lab that are connected to the database. I have used the MultiQML tool to symbolize the layers and saved them on the server, which did not load the correct symboled features on the connected lab computers. Do I have to have the QML files on each lab computer, or is the a way to have a QML file on the server that each lab computer will be able to read via the database connection and symbolize the layers correctly without needing to store the QML locally? Thanks, Ryan


How to reverse geocode without Google?

I have a basic tracking system i've written with OpenLayers using google and openstreetmap and I want to use geocoding to get street names etc However, I'm a bit confused. I was hoping there was a easy/free way with OL and OSM however, i'm not sure there is I've tried using google, but it doesn't work with bulk data. It allows me to get about 3 or 4 back and then doesn't work for the rest..... i think they require a wait period between each request anyway, can anyone recommend a free way to geocode without google? please be patient with my lack on knowledge in this areas Thanks


Mosaicking a selection from a raster catalog and returning the mosaicked raster

I have pretty much given up on trying to get help from ESRI forums so I thought I would try my luck with the good people here... Here is what I have (in C#)... I have an IFeatureCursor that is giving me a selection from an IRasterCatalog (so basically it is a cursor full of IRasterCatalogItems). What I would like to do is take these rasters and mosaic them inside a method and return the mosaicked raster. I have tried what I thought would be obvious beginning with an IRasterCollection...               IMosaicRaster m_mosaic = (IMosaicRaster)m_rasColl;     m_mosaic.WhereClause = "";          ISaveAs m_saveAS = (ISaveAs)m_mosaic;     IRasterDataset rasDataset =             (IRasterDataset)m_saveAS.SaveAs("tmp_" + selectedFeatureClass.FeatureClassID.ToString(),                                            tmpWS, "TIFF");           IRaster extractedRaster = rasDataset.CreateDefaultRaster();          // Compute the statistics for the raster.     IRasterCalcStatsHistogram m_calcStats = new RasterCalcStatsHistogramClass();     IStatsHistogram m_stHisto = new StatsHistogramClass();          m_calcStats.ComputeFromRaster(extractedRaster, 1, m_stHisto);          return extractedRaster;      This gave me what seemed to be a properly mosaicked raster, and the values seemed to be correct when I used the identify tool in ArcMap, however, the display looked like a raster that has not had statistics calculated and I can't seem to figure out why. So I started looking at less obvious ways to do it but I haven't had much luck so far. Suggestions? Thanks a bunch!


How to get rpy2 for QGIS on OSX?

I am a new QGIS user, have installed QGIS on a Mac running Snow Leopard OSX 10.6.8 following all the instructions. QGIS seems to work well for everything other than when I boot up I get the message "Unable to load the plugin: Unable to load required package rpy2. Please ensure that both R, and the corresponding version of Rpy2 are correctly installed". I have installed Rpy2 several times and reloaded R. I have R version 2.15 loaded (but also RStudio). Does anyone have any ideas?


How to get FGDB support in GDAL 1.9 without compiling?

I was hoping that I could use the 'master' (nightly-build) version of QGIS to get fgdb support without compiliing (I'm not a developer and frankly have wasted days of productive time trying to get things to compile with no luck). Since it seems that gdal 1.9 support fgdb, I thought I could just download the developer version. However, the latest version of QGIS (191a229), even with gdal 1.9 does not seem to allow me to add my gdb file. Is the gdal version included with the master build just not compiled with fgdb support? Does anyone know if there's a way to do this without compiling a new version?


Georeferencing a whole heap of raster imagery

I have a whole heap of satelite imagery coming in that is not georeferenced. They have a grid overlaid so I do know the bounding box/coordinates of the images (and projection). I just need them georeferenced for use in ArcGIS. I have manually done a few by referencing to a point dataset I traced over the top - and it works great. I just need to automate this process. Each image is exactly the same, so the same points overlaid will always match up (even if it doesnt, the application of this does not need to be accurate). Basically I need to apply the 'link table' information from the example I've already done, and apply it to all the images... Possible in python/model builder?


Naming your Output feature class?

I am working in a work space in python on ArcGIS 10.1 on a remote desktop. Here is my code               import arcpy     from arcpy import env     env.workspace = "E:\Example"          PST = arcpy.ListFeatureClasses ("", "Point")     LST = arcpy.ListFeatureClasses ("", "Polyline")          for point in PST:         for polyline in LST:             arcpy.SplitLineAtPoint_management (polyline, point, point + "500" + ".shp", 1)      I want the output feature class to have the same file name as the point feature class, but with 500 appended to the end. Currently the way the code is now, it will name the output "Example.shp500.shp" which is giving me an error of failing to execute. How can I name my output "Example500.shp"?


Need software to create Info-graphic map

I want to make an info graphic map for history theme hence I am researching on Google for suitable software but not able to find suitable software and very less information regarding available related with info-graphic maps. We are using Open source platform ,while searching I found this map and looking similar kind of maps. Found following useful links : Making info-graphics in ArcMap Blog :Infographics, GIS, and Mapping So does anyone knows about info-graphic maps ? & its software/reference sites ?


How to get the result of a geometry service relation method?

I'm using ArcGIS 10.0 flex API 3.0. I'd like to test if two geometries are equal. For that, the relation method is called like follows:                       var relationParameters:RelationParameters = new RelationParameters();             relationParameters.geometries1 = [geom1];             relationParameters.geometries2 = [geom2];             relationParameters.spatialRelationship = RelationParameters.SPATIAL_REL_RELATION;             relationParameters.comparisonString = "IDENTICAL (G1, G2) = TRUE";             _geometryService.addEventListener(GeometryServiceEvent.RELATION_COMPLETE, onGeomCompared);             _geometryService.relation(relationParameters);      The service is supposed to return a boolean that tells if the two geometries are equal or not. Isn't it? The GeometryServiceEvent returned by the service doesn't contain any boolean and the result object of the event is empty. How can I get the response of the service for my request?


How to calculate distance in meter between QgsPoint objects in PyQgis?

I tried to calculate distance between QgsPoint objects in qgis, i used QgsPoint `StartPoint.sqrDist(EndPoint)`. but i don't find an expected result, i also tried               import math     math.sqrt(StartPoint.sqrDist(EndPoint))      but still getting a very high result. does somebody know a way how to calculate distance in meter between 2 Qgspoint objects.


Wrong object size after projecting spatial data from one coordinate system to another

I am having problems reprojecting items in ArcMap 10.1 from British National Grid To WGS84. (WKID: 27700 Authority: EPSG to WKID: 102100 Authority: ESRI) All Polygon Layers I project [I use Project (Data Management)] from BNG to WGS are basically twice as big in metres as they were originally. Even the Shape Area and Shape Length inside the attribute table changes. Any ideas on this would be greatly appreciated. Thank you.


Extracting random point coordinates using ArcGIS 10 model builder

I have a question about using model builder in ArcGIS 10.1. In my model, I am creating a processing extent, and within that I am creating a random point to be used in further processing. This model is iterative, and for each iteration, I am creating a new random point (effectively overwriting the values for the point).   What I would like to be able to do is to write the x-y coordinates of these random points off to a file before the next iteration is performed (either all together or as separate files for each point).   In a way, what I want is the inverse of the "Get Cell Value" tool and to be able to write that value to a file.   I guess it wouldn't be too difficult in python, but I was hoping to be able to do this in model builder.


How do I truncate a fieldname in model builder using field iterator

I have about 9,000 tables with rather long fieldnames, I need to chop off part of the fieldname, so for instance I want to change 'GP_CodeA82010' to 'A82010'. I am doing this using field iterator and the calculate field, but the calculate field does not seem to like the iteration, so for instance instead on 'A82010' it gives me an output value of '1'. I have tested a few python expression on the calculation field; Right( [FieldName], 6) and slice function- "%Name%"[7:13] Does anyone know a work around this?


"car Icon moving in a polyline from start to end destination"

I am new in google maps api programming and I want to make a map which is the same as this map example "http://econym.org.uk/gmap/example_cartrip2.htm)" does anybody here can help me to start and build this kind of application for android? tnx... :)


Is there a way to find out which pixel values are around the current pixel?

I would like to write a cycle which takes a current pixel, look at the surrounding pixels, choose the pixel with max/min value and skip to that pixel...and so on. Is there in arcpy such a function which can work with pixel's neighborhood?


Do ArcGIS clients that direct-connect to postgresql/postgis databases override configuration settings?

I went in to modify some settings on a Postgresql database via the PGAdmin config editor and some of the settings don't have their default values and they are higher because they've been overridden (e.g. `maintenance_work_mem` ). I know PGAdmin will do this. And that I can do it manually in PSQL. Do the ArcGIS clients ever override postgresql configuration settings? And if so, what triggers it? Is it tool-based? Or when the software starts up? Or is Esri entirely hands-off?


GeoTIFF Mosaic/Crop/Translate: Easiest Method in QGIS?

I routinely have to mosaic multiple GeoTIFFs together, extract a rectangular area, and save that area as a JPG with a world file. Typically I've been using commercial software for that (Geographic Imager) but would like to switch that task over to QGIS. I can eventually get the task completed using multiple tools, but the process takes much longer and involves creation of intermediate GeoTIFFs along the way. In Geographic Imager I simply load the images into Photoshop, crop, and SaveAs and I'm done. So what's the easiest, quickest workflow to:   1. Load multiple GeoTIFFs   2. Mosaic them together   3. Select an area to extract (preferably by drawing a box)   4. Save selected area out to JPG with accompanying world file Thanks!


How to convert a string from GetParameterAsText to uppercase?

This is my first question here and I'm sure it's ridiculously simple. But I have Googled and no luck. I am trying to make a tool in ArcToolbox using a Python script so users can input a street and lot number and it will zoom to the correct location/s in the cadastre. It's working but I have to type the name of the street in upper case in the toolbox and I'd like it work for uppercase or lower. I assume I should use something like .upper() but I don't know where to put it. I've tried and failed in a few places. Here's my snippet:               RoadName = "'" + arcpy.GetParameterAsText(0) + "'"


What are the parameters for my connection QGIS

I am new QGIS, I was wondering what the parameters are for the New MSSQL Connection in the dialog window? Provider/DNS: `? when local ?` tried leaving empty Host: `? when local ?` tried my computer user name/server instance name Database: `full path to file I assume` ![Dialog Window](http://i.stack.imgur.com/djIky.png) Error: ![Error](http://i.stack.imgur.com/j1GdA.png)


Is there a way to automatically select a feature when another feature in a different layer shares attribute data?

I have a mapset displaying foraging areas, each with a unique ID. I spatially joined the mapbooks layer, so each mapbook sheet identifies the foraging area it is depicting (1:1), so sheet 5 is foraging area 18. So when sheet 5 draws, I want it to select foraging area 18, which has a matching attribute. I want to do this so I can change the symbology of the target feature to make the map stand out more. I'm using ArcGIS 9.3.x with the mapbooks extension.


How can I remove dipoles from rasterlayer in QGIS 2.2?

I work with geomagnetic survey data. When convertet into a rasterlayer, every object with a high magnetic potential (for exampel an iron-nail on the surface) is shown as very sharp anomalie called dipole (a black and a white circle tied together). These dipoles disturb the further analysis of the rasterdata. Normally I mark them with dots in a second vectorlayer, but this is a very timeconsuming work. Is there any extension, with wich I can filter these dipoles directly from the rasterlayer?


How to map factors and their relative importance in ArcGIS 10.1 for Desktop?

I want to know how I can calculate a map from a set of factors (in the form of raster layers) and Their relative importance for example: slope 20% lithology 10% . . .


exporting segmentation results from ecognition to other software for anaylsis

I want to create vegetation indices in ERDAS but based on segmentation derived in eCognition. I have tried different ways of exporting my segmentation results to a TIFF file which i can further use for analysis, however I am not getting undesired results.


Bounding Box in CuRL Tool(REST API)

I have set my time interval using CuRL tool for automatic Geowebcache.. How to set bounding box in curl command?


gdalwarp needs bounding box

I try to mask a raster with gdalwarp and I always get the answer: cannot compute bounding box of cutline. what is wrong with my cutline shapefile? or what requirements should it have to work?


User defined Polygon disappears with featureadded and selectfeature

I am trying to achieve the following functionality:   1. User draws a polygon.   2. Once complete, fetch features in that polygon( Using WFS).   3. On Select polygon, show the features that intersect that polygon. My issue is: Once user draws a polygon, I use featureadded to fetch the features within the polygon. Once fetched, I can see the details using alert with no issue. But the circle/polygon disappears once I start zooming in /out. I want the polygon to be still visible. Again on select, the polygon disappears as soon as the polygon is selected. I want the polygon to be visible at all times and only be removed when I delete it. I am using this example as reference : http://openlayers.org/dev/examples/draw-feature.html Any inputs on how I can have the polygon to not disappear? Here is the code I am using:               var map, drawControls, polygonLayer;     var url = "http://localhost:8080/geoserver/wms";     var wfsurl = "http://localhost:8080/geoserver/wfs";     var bounds = new OpenLayers.Bounds(-20037508.34, -20037508.34, 20037508.34,20037508.34);     var EPSG4326 = new OpenLayers.Projection("EPSG:4326");     var EPSG900913 = new OpenLayers.Projection("EPSG:900913");          function init() {      var options = {          controls: [],          maxExtent: bounds,          restrictedExtent: bounds,          maxExtent: bounds,          numZoomLevels: 4,          maxResolution: 156543.0339,          projection: EPSG900913,          units: "m",          displayProjection: EPSG4326           };      map = new OpenLayers.Map('map', options);      // Using OSM as base map      var layer = new OpenLayers.Layer.OSM(          "OpenStreet Maps");                states = new OpenLayers.Layer.WMS(          "NA states", url, {          layers: 'topp:states',          tiled: true,          styles: '',          transparent: true,          format: 'image/png'      }, {          isBaseLayer: false,          displayOutsideMaxExtent: false,          ratio: 1,          buffer: 0,          yx: {              'EPSG:4326': true          }      });           polygonLayer = new OpenLayers.Layer.Vector("Polygon Layer", {          displayInLayerSwitcher: false      });      circleLayer = new OpenLayers.Layer.Vector("Circle Layer", {          displayInLayerSwitcher: false      });      map.addLayers([layer, states, polygonLayer, circleLayer]);      map.addControl(new OpenLayers.Control.PanZoomBar({          position: new OpenLayers.Pixel(2, 15)      }));      map.setCenter((new OpenLayers.LonLat(4, 47).transform(EPSG4326, EPSG900913)), 3);      map.addControl(new OpenLayers.Control.Navigation());      map.addControl(new OpenLayers.Control.MousePosition());      map.addControl(new OpenLayers.Control.Scale($('scale')));      map.addControl(new OpenLayers.Control.ScaleLine({          geodesic: true      }));      map.addControl(new OpenLayers.Control.LayerSwitcher());           selectControl = new OpenLayers.Control.SelectFeature([polygonLayer, circleLayer], {          onSelect: onFeatureSelect,          // onUnselect: onFeatureUnselect      });      drawControls = {          polygon: new OpenLayers.Control.DrawFeature(polygonLayer,          OpenLayers.Handler.Polygon),          circle: new OpenLayers.Control.DrawFeature(circleLayer,          OpenLayers.Handler.RegularPolygon, {              handlerOptions: {                  sides: 40                   }          }),          select: selectControl      };      for (var key in drawControls) {          map.addControl(drawControls[key]);          drawControls[key].events.register('featureadded', drawControls[key], function (f) {              onFeatureAdded(f.feature);          });      }      map.events.register('click', map, function (e) {          var bounds = map.getExtent().clone();          bounds = bounds.transform(EPSG900913, EPSG4326);          mouseLoc = map.getLonLatFromPixel(e.xy);               OpenLayers.Request.GET({              url: url,              params: {                  REQUEST: "GetFeatureInfo",                  EXCEPTIONS: "application/vnd.ogc.se_xml",                  BBOX: bounds.toBBOX(),                  SERVICE: "WMS",                  INFO_FORMAT: 'text/html',                  QUERY_LAYERS: [map.layers[1].params.LAYERS],                  FEATURE_COUNT: 50,                  Layers: [map.layers[1].params.LAYERS],                  WIDTH: map.size.w,                  HEIGHT: map.size.h,                  x: Math.floor(e.xy.x),                  y: Math.floor(e.xy.y),                  buffer: 0,                            //format: format,                  styles: map.layers[0].params.STYLES,                  srs: map.layers[0].params.SRS              },              success: setHTML          });               OpenLayers.Event.stop(e);           });     }          function setHTML(response) {           document.getElementById('nodelist').innerHTML = response.responseText;          }            function toggleControl(element) {        for (key in drawControls) {            var control = drawControls[key];            if (element.value == key && element.checked) {                control.activate();            } else {                control.deactivate();            }        }        if (element.value == "delete") {          polygonLayer.destroyFeatures();          polygonLayer.redraw();          circleLayer.destroyFeatures();          circleLayer.redraw();      }           }          function onFeatureAdded(e) {           var wfsProtocol = new OpenLayers.Protocol.WFS({          url: wfsurl,          featureType: 'states',          featureNS: "http://www.openplans.org/topp",          srsName: EPSG4326,               geometryName: "the_geom"      });                var pfilter = new OpenLayers.Filter.Spatial({          type: OpenLayers.Filter.Spatial.INTERSECTS,          value: e.geometry.transform(EPSG900913, EPSG4326)      });           wfsProtocol.read({          filter: pfilter,          callback: processSpatialBboxQuery,          scope: new OpenLayers.Strategy.BBOX()                });                  }        function onFeatureSelect(e) {           var cql = "INTERSECTS(the_geom," + e.geometry + ")"      alert(cql);      states.mergeNewParams({          'CQL_FILTER': cql      });          }     // Process the features fetched. For the time being we are just showing  POS Name of first feature         function processSpatialBboxQuery(e) {      alert("hi");      };

             <input type="radio" name="type" value="none" id="noneToggle"                 onclick="toggleControl(this);" checked="checked" />                 Navigate
                 <input type="radio" name="type" value="polygon" id="polygonToggle"                 onclick="toggleControl(this);" />                 Draw Polygon
                 <input type="radio" name="type" value="circle" id="circle"                 onclick="toggleControl(this);" />                 Draw Circle
                 <input type="radio" name="type" value="select" id="selectToggle"                 onclick="toggleControl(this);" />                 Select feature
                 <input type="radio" name="type" value="delete"                 id="delete" onclick="toggleControl(this);" />                 Delete feature
                <input type="checkbox" name="allow-pan" value="allow-pan" id="allowPanCheckbox" checked=true onclick="allowPan(this);" />                 allow pan while drawing


Difficulty in running the script

I have developed a script tool for adding a layer in the current map document. It did not work.. The message is "Completed script addlayer...". But there were no layers added in the TOC. Another issue is when i type the same script in the python console window within Arcmap, the layer gets added. After adding the layer through console window, when i execute the script tool, it works and the layer gets added. I don't know what is the problem... pls help.. this is my script..               import arcpy     mxd = arcpy.mapping.MapDocument("CURRENT")     df = arcpy.mapping.ListDataFrames(mxd, "*")[0]     addLayer = arcpy.mapping.Layer(r"D:\Academic\M.E_thesis\data\ascii files\maxent\soil.lyr")     arcpy.mapping.AddLayer(df, addLayer, "BOTTOM")


Can Python toolboxes be used with ArcObjects GeoProcess.AddToolbox? I'm trying to execute a tool in a Python toolbox (.pyt) from .NET. First I call:               GP.AddToolbox(@"C:\path\to\MyToolbox.pyt")       And that call doesn't get an error. Then I build the parameter list and then I execute it like this:               result = (IGeoProcessorResult2)GP.Execute("ToolName", parameters, null);      This immediately throws the "Error HRESULT E_FAIL has been returned from a call to a COM component." error, so nothing is returned to result. My next step is to create a dummy script and build it into a .tbx toolbox and leave the .NET code the same to see if it's something else, but it did make me wonder if this is even possible. I couldn't fine any examples of using AddToolbox with a .pyt file.


Randomly shuffle values from a column into a series of new columns in a feature class in ArcGIS for Desktop?

I am trying to run a Monte Carlo simulation on a point feature class. In order to do so, I am trying to take the values from a column of observed data and randomly "shuffle" them between the points over n additional "simulated" columns. In my case, I'd like to add 199 new columns, each with the values of the the observed data randomly redistributed among the points. My observed data are binary (1 and 0) I have been able to do this in R previously and bring the table back into ArcGIS, but the tables have become so unwieldly that I am having difficulty bringing them back into ArcGIS, and would like to run the analysis in ArcGIS, rather than moving in and out. Does anyone have any idea how to do this? I have tried using cursor, but have failed.


Comparing two maps with nice web interface?

Is it possible in Leaflet or OpenLayers to compare two maps with a bar to slide between the two? Examples include:   * http://www.esri.com/services/disaster-response/hurricanes/hurricane-sandy-the-aftermap   * http://www.abc.net.au/news/specials/qld-floods/ (not a map) I know it's possible, I don't know how to do it though. Anyone got any ideas? Examples?


Problem in splitting and removing Polyline

![Two polylines intersecting](http://i.stack.imgur.com/rD4Mj.png) I am stuck at a problem for weeks, Please Help!!! I am making program in C# Here is what I need to do --- I have two polylines that are intersecting with each other at some points. I want to delete those portion of one polyline which are going above the other polyline. I am attaching an image explaining the problem I am going through. In this image, I have two Polylines Line1 and Line2. So according to that image, I want to remove the **dashed** portions of Line2. Please Help!!! Thanks a ton in advance Here is what I have done till now.... First I am finding out all the points of intersection of the two lines and saving them in a layer.... Then I am Splitting these lines at these points using this code...               while (pPointF != null)     {         IPoint pPoint = (IPoint)pPointF.Shape;         ISpatialFilter pSF = new SpatialFilterClass();              pSF.Geometry = pPoint;         pSF.GeometryField = "Shape";         pSF.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects;              IFeatureCursor pLineCursor = pFeatureclass1.Search(pSF, true);         IFeature pLineF = pLineCursor.NextFeature();              while (pLineF != null)         {             IPolycurve pPolycurve = (IPolycurve)pLineF.Shape;             IPoint pToPoint = pPolycurve.ToPoint;             IPoint pFromPoint = pPolycurve.FromPoint;                  if (pFromPoint.X != pPoint.X && pFromPoint.Y != pPoint.Y && pToPoint.X != pPoint.X && pToPoint.Y != pPoint.Y)             {                 IFeatureEdit pFeatureEdit = (IFeatureEdit)pLineF;                 pFeatureEdit.Split(pPoint);             }             pLineF = pLineCursor.NextFeature();         }     }      The problem with this code that it is giving different results than expected. As I can see from attribute table, it is generating more than expected parts of line i.e. if there are 5 intersection points then there should be 6 parts of a line, nut it is creating 10-12 parts. Also the number of parts created for both lines are different, which logically should be same. Due to these problems, when I am trying to delete unwanted parts of line, the arcGis crashes..... Please help with this problem. . . . Now I am using a Geoprocessor to do splitting, following is the code:                       SplitLineAtPoint pSplitter = new SplitLineAtPoint();             pSplitter.in_features = pPolylayer1;             pSplitter.point_features = pPointsFeatureLayer;             string outfile = txtDirName.Text + "\\" + "splitted";             pSplitter.out_feature_class = @outfile;                  Geoprocessor GeoPro = new Geoprocessor();             GeoPro.Execute(pSplitter, null);      But on the last line, where Geoprocessor is being executed, the ArcGIS crashes... Please tell me where is the problem...


Replicating QGIS style in R with rgdal

I want to use rgdal to plot a shapefile graduating the colour of features according to the values of one field (`density`). I couldn't find online detailed indications on how to do it and the package documentation doesn't explain the options for `plot()`. This is what I composed so far               library(rgdal)            shapefile = "shapefile.shp"     map <- readOGR(dsn = shapefile, layer = "shapefile", verbose=FALSE)          plot(map, xlim = c(6.70, 18.32), ylim = c(35.2, 47.6), border=NA, add=TRUE)      This is the style in QGIS I am trying to replicate: ![enter image description here](http://i.stack.imgur.com/hihky.png)


Field calculator fails when computing First X coordinate in ArcGIS10

I have this line in python to compute First X coordinates into the Field Calculator from each polyline:               !Shape!.firstPoint.X      Then, when I inspected the attributes I see this rare thing: ![Example](http://i.stack.imgur.com/JAiwX.png) You may see that I have two fields ( _StartX_ and _InicioX_ ), both with the same Python code, one is a string field and the other a Float one. But **the computed VALUE is wrong** :( **You can clearly see that "spatially" the features are situated in different positions, but on its computed X Value are the same.** **Could someone explain me that?** Does anybody know any other way to compute this X,Y start and end coordinates from lines? Thanks for your time.


Is it possible to connect to FEMA's WMS and WFS using QGIS?

I have just found out about the wonderful world of WMS and WFS. I have added the FEMA versions of both to QGIS. WMS and WFS. The WMS seems to connect and adds something to QGIS, yet it is not visible. The WFS connects, but I get an error that it is not a valid layer. So... I guess that my question is - Is anyone willing to test these and maybe give some guidance on how to better use these? I have also tried NYS WMS and the Nationwide Inventory WMS `http://107.20.228.18/ArcGIS/services/FWS_Wetlands_WMS/mapserver/wmsserver?` with no luck. Similar results to FEMA's WMS.


Counting road intersections and blocks

I have a shp file that consists of the road system of Mexico City. I would like to count the number of road intersections (crossings) and blocks (which I would define as a segment starting from one intersection and ending in another). This within a buffer of a certain size: a circle centered in a polygon within another shp file. I am mainly using ArcGIS but am quite familiar with R, Python and QGIS if the solution is simpler on any of those. I realize their are a few similar questions around here but havent been able to use any of the information provided there.


Layers with same coordinates not lining up [QGIS]

I have a project where I have a few layers (State of Texas, cities, roads, etc.). They are all resaved and projected to NAD83 / UTM zone 14N EPSG: 26914. The project CRS is also NAD83 / UTM zone 14N EPSG: 26914 and on the fly projection is enabled. I added a text delimited file and cannot get them to line up. I checked the attribute table of the delimited file and the coordinates are correct. I cross-checked this with a point close to where a few of the points from the delimited file are and they are very close so in theory they should be lining up. I tried to save the delimited file and save it with the same projection and add the shp file back in and that did not work. Anyone have any ideas? I can't think of anything else to try.


pyscripter out of memory

I'm running a loop that create feature classes at a geodatabase. After creating 221 datasets I get the error message "pyscripter out of memory". Is there a way to clear the physical memory of computer in order to resume with the script ? I've alredy tried :                 gc.enable       clear = lambda: os.system('cls')       clear()       gc.collect()      ![Windos Task Manager](http://i.stack.imgur.com/jTEJp.png) ![enter image description here](http://i.stack.imgur.com/5KQhI.png) ![enter image description here](http://i.stack.imgur.com/Eh94I.png) ![Out of Memory in pyscripter](http://i.stack.imgur.com/UBUV9.png)


Measuring distance between discharge point, along river, until next intersection with another river?

i hope someone can answer me my simple question. I have two sets of data   1. a river network --> line feature   2. discharge points along those rivers --> Point Data The discharge points represent wastewater treatment plants. I try to find those treatment plants which don't have a inflow of another river, for a minimum of 2 km downstream. This means the river from this point on does not intersect with another river for at least 2km. How can i find those plants? Are there any tools for this? What about the Arc Hydro tools?


CartoDB table seems to be corrupt

Have 218 rows in a table in the free version of CartoDB. Suddenly behaving strangely - not accepting edits to fields in existing rows, only showing 200 rows when sorted Descending on row number (but showing 218 when Ascending), etc. Is there a way to refresh the database or at least find out where the problem is?


how to refresh a layer using openlayer

in geoserver I configured a layer containing positions that change in time (they are smart phone positions indeed). I'am able to insert fresh position in geoserver datastore (a postgis db), but I'm not able to automatically reload the layer in client side using openLayers. In addition I wish to load from the layer only the information contained in the visualized BBOX (ie what the user will see on its screen). Could anyone provide me a suggestion?


Develop 3D buildings with roof structure from raw data

I recently received a task to develop 3D buildings with roof structure. All I have now is satellite images and aerial photographs covering the entire urban area. After reading some previous methods, it seems possible to achieve this with building/roof edges delineated and height information (e.g. a digital surface model). Is it possible to perform these tasks in ArcGIS? If so, with which tools? I am pretty new to GIS so I am not familiar with tools involved in such tasks.


GeoServer setup for intercepting WFS-T requests and handle them via PHP

I'm developing a web app which uses GeoServer side by side with a PHP framework. My goal is to be able to intercept WFS-T requests to update some records whenever a WFS transaction is executed (insert/update/delete). Eg: I want to be able to keep track of modifications to the features saving `user_id` of who made the change and duplicating the affected record before the changes take place to keep an history of changes on features. My actual setup is as follows:   * apache proxies all `/geoserver/*` requests to port 8080 (Tomcat/GeoServer) except for those directed to `/geoserver/wfs/*`   * these last (wfs) are managed by my php framework (they are first checked for session authentication) which reads the request XML body and backs up the interested records accordingly before passing on the request via php cURL I'm quite satisfied of the working for now but...   * I'm not sure it's a good approach (seems a bit complex): I wonder if there's a built in way of managing record backups in GeoServer   * I wonder if there could be security implications in leaving all requests (except for WFS ones) go directly to GeoServer Any suggestion either on php or on GeoServer would be greatly appreciated.


Binding licence for using ArcEngine

I have just started developing with arcengine and some how I am stick in first steps.. When I want to bind licence to start appplication I get an error which means I haven't any licence..(I think so)..Here is program.cs code                   using System;     using System.Collections.Generic;     using System.Linq;     using System.Windows.Forms;     using ESRI.ArcGIS;     namespace WindowsFormsApplication4     {         static class Program         {             /// <summary>             /// The main entry point for the application.             /// </summary>             [STAThread]             static void Main()             {                // ESRI.ArcGIS.RuntimeManager.BindLicense(ProductCode.EngineOrDesktop);                           if (!RuntimeManager.Bind(ProductCode.EngineOrDesktop))                 {                     MessageBox.Show(                         "Unable to bind to ArcGIS runtime. Application will be shut down.");                     return;                 }                 Application.EnableVisualStyles();                 Application.SetCompatibleTextRenderingDefault(false);                 Application.Run(new Form1());             }         }     }      I checked Arcgis Administrator and here is licences that I have..Can you please help me about binding licence?![ARCGIS Licence ](http://i.stack.imgur.com/EbYEt.png)


Parcel ID Mapping (brainstorming)

In my project I have a table full of residential properties by Parcel_ID. Table 2 is a TIGER streets shapefile. I would like to build a query in MS Access displaying all properties within a mile radius of my selected property from table 1. ![project structure](http://i.stack.imgur.com/icQDh.jpg) My questions are :   1. Should I use open-source postgis (spatial databases) and build this project in SQL 2008 to query off map location rather than tabular data?   2. Or Should I use a FileGeodatabase to import the shape file into MS ACCESS?  Thanks for any help or reply, PS - I have the shapefile loaded in ArcGIS explorer - but not sure how to export. 6/17/2012 update* _*_ I posted the state parcel tables below: **County Mapping Parcel Table**   OBJECTID PARCEL_TAG **PARCEL_I** PARCEL_C LOTNUM SUBDIVISIO SUBDIVNUM STNUMBER PRE_DIR STREET_NAM SUFFIX SUF_DIR FULL_STNAM CITY COUNTY STATE TOWNSHIP ZIPCODE ESTSQFT BLOCK MODDATE STATUS **SHAPE_AREA** **SHAPE_LEN** | 1 1251683 8011477 8011477 7809 1723 4015 N KENWOOD AVE KENWOOD AVE INDIANAPOLIS MARION IN WASHINGTON 46208 5280 1004 **** _5279.952157_ *** **_343.9990607_**


How to use a spectral library

I want to classify an Aster image using the ASTER Spectral Library - Version 2.0 Anyone here ever done this? Any suggestions? Note: I'm using opensource GIS and Remote Sensing software Thanks.


Why does the File Geodatabase that I created not appear to any other users on the network?

I created a File Geodatabase to store files that I use in maps for my coworkers. The folder is in a location where all users should have access and permissions. Strangely enough, after I created the gdb, I'm the only one that can see or access it. It appears in my Catalog view, but not in windows explorer. The map document that accesses the fcs in the gdb works fine for me but shows broken links for my colleagues. Any reason why this might be happening?? ![Where's my database?](http://i.stack.imgur.com/7YPlW.png)


How to Export a Shapefile as KML Using R

All, I have exported an ESRI shapefile (polygon) to KML using R, but the resulting layer shows only the outer extent of the polygon. How can I fill in the area with a color? Any help is greatly appreciated. The code is as follows:               library(rgdal)     setwd("G:\\GIS_SJR\\GIS Shapefiles\\Boundaries")     polygon <- readOGR(".", "20km_buffer_albers")     polygonWGS <- spTransform(polygon, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))       writeOGR(polygonWGS, dsn="polygonWGS.kml", layer="polygonWGS", driver="KML")


How to add a graphic table element in Python without Production Mapping extension?

I'd like to add a table in the Map Layout that includes only items visible at the current extent of the dataset. I know that the Production Mapping Extension lets you insert a Graphic Table Element filtered by the dataset current extent (I suppose). But we haven't this extension. Has anyone ever tried to do something like this in Python?


Matching Landsat 8 with previous Landsat images

I am working with Landsat 8 in ENVI. My image's dimensions are 8921 x 13131. I need to apply a forest mask to this image with dimensions of 8805 x 12321. I try to resize an image, but it doesn't work. (After that I need to calculate forest disturbance index for 2013 and subtract it from DI for 2000, so I really need an 2013 image (Landsat 8) to be resized). I resize it using basic tool - spatial subset - subset using - file - and then pick my forest mask. However, as a result I get a file of 8806 x 12321 instead of 8805. If I adjust the number of samples manually, it gives me 8805. But then I try to apply a mask and that's what I get (see a screenshot). They just don't match. And when I try to overlay my mask (which I made based on Landsat 7 image) with my Landsat 8 image, those two images have different frames. Does anybody have any good advice? Thanks a lot. ![enter image description here](http://i.stack.imgur.com/OOFMu.png)


Simplify Multigeometry polygons in QGIS

I am really new to QGIS and KML so this might be obvious but I am too ignorant to know what to do. I have some KML files I would like to use in an application but they are much more complicated than what I need. They contain multigeometry with a number of polygons and I would like to merge these regions to simply get the outline of the total region. In addition, I would like to reduce the number of vertices. ![Displaying layer](http://i.stack.imgur.com/ATYxc.png) I saved this as a shapefile so I could access the Merge Selected Features tool and I merged all of the features: ![Merged features](http://i.stack.imgur.com/fiGlh.png) But when I try to simplify the feature I get the error: > This feature cannot be simplified. Check if it has enough vertices to be > simplified In the documentation for Simplify Feature I saw it mentioned: > If a feature cannot be simplified (e.g. MultiPolygons), a message shows up. So I thought this was likely the problem. Is there any way I can simplify this, or convert the multipolygon to a simple polygon? Edit: During my searches I found mention of a plugin called fTools that seems to have some possibly useful functions, such as "Multipart to singleparts" and other vector tools that would be useful, however I do not seem to have it by default in QGIS 2.2, nor can I find it on the qgis plugins page. Is this plugin gone?


WFS protocol - Update method

I have a question related to the Update method, in case the feature doesn't exists, the update method will generate a new feature? (i.e. will behave as Insert)? Or maybe, in contrast, an insertion using an existing fid (feature id) will perform an update over this feature? Thanks in advance.


evt.layers undefined while trying to display Legend widget in Javascript API

I cannot get an ArcGIS Javascript API Legend widget to show up. An excerpt of my code is as follows. I believe I'm failing somewhere near the console.log that I call...undefined comes up in the console at this point. So for some reason evt.layers is not getting populated. Have I provided enough code? If you need more snippets I can copy/paste more.                       parser.parse();                  //initialize the map object             map = new Map("map", {                 basemap:"gray",                 center: [-96.53, 38.374],                 zoom: 13             });                  //access the hydro layer as a feature layer             var hydro = new FeatureLayer("http://sampleserver3.arcgisonline.com/ArcGIS/rest/services/Hydrography/Watershed173811/MapServer/0", {                 mode: FeatureLayer.MODE_ONDEMAND,                 outFields: ["*"]             });                  //add hydro to map             map.on("layer-add-result", function (evt) {                 console.log(evt.layers);                 var layerInfo = arrayUtils.map(evt.layers, function (layer, index) {                     return {layer:layer.layer, title:layer.layer.name};                 });                 if (layerInfo.length > 0) {                     var legendDijit = new Legend({                         map: map,                         layerInfos: layerInfo                     }, "legendDiv");                     legendDijit.startup();                 }             });             map.addLayers([hydro]);         });


Why aren't the info windows showing up on my cartodb layers?

I made a cartodb map with 2 layers. The layers draw from the same table but I use filters to display some of the data in one layer and some of the data in another. Without the layers my info windows were working and now they're not. I'm not sure how to start troubleshooting this. The SQL expression on layer 1 is `SELECT * FROM table_name WHERE useful_column > 0` and on layer 2 it is `SELECT * FROM table_name WHERE useful_column IS NULL` -- the goal is to use a different info window altogether on the regions where `useful_column` has no data.


How to read out the geoPDF data downloaded from USGS ned database?

I am new to GIS, so this might be a stupid question. Google foo din't get me anywhere, so here goes.. How do I read the geoPDF data that I download from the National elevation database from USGS? All I am looking for is lat/long/height for each point as much as the resulution would allow for a specific area in the map. Are there ways to get this data in a better/more friendly format than geoPDF? (I am good with python, so any solutions that are py-friendly would be awesome, feel free to toss some search words at me that would help me further my knowledge) thank you for your help


Does Collada allow to georeference the 3D coordinates within a DAE file?

Can I assign lat/long/height values to each vertice of a Collada 3D model?


How to use Time as a parameter in Network Analyst Service Area and Closest Facility

I have a road network downloaded from Ordnance Survey (Meridian 2, with A roads, B roads and Minor roads Merged into one shapefile and then built into a Network Dataset for use in ArcGIS' Network Analyst) and would like to use it to determine Service Areas for 4 hospitals (preferably with drive times of 30 minutes and 60 minutes) and Closest Facilities (same 4 hospitals and directions from 539 points). I've used the Network Dataset to find the Service Area and the Closest Facility but the only parameter I can use is distance (metres). Question: **How do I set up the network to use drive time, not distance, and work out Service Area by, say, 30 minutes?** When I try to create the Network Dataset from the downloaded road network, I try to add Drive Time as an Attribute but I can't get the Network Dataset Wizard to accept the evaluators. I've tried to follow the ESRI help pages and I've Googled around but I'm lost. I'm not sure what the time setup should be or if it's something I can just do or something I need to set up from the start. ![enter image description here](http://i.stack.imgur.com/5B2DG.png)


How to measure the accuracy of latitude and longitude?

I am very new to GPS stuff. Suppose I have latitude and longitude as `19.0649070739746` and `73.1308670043945` respectively. In this case both coordinates are `13` decimal places long, but sometimes I also get coordinates which are `6` decimal places long. Does fewer decimal points affect accuracy? What does every digit after the decimal place signify?


How can I tell the accuracy of coordinates conversion by cs2cs?

For example, when I try this:               cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +to +proj=tmerc +lat_     0=22.31213333333334 +lon_0=114.1785555555556 +k=1 +x_0=836694.05 +y_0=819069.8 +     ellps=intl +towgs84=-162.619,-276.959,-161.764,0.067753,-2.24365,-1.15883,-1.094     25 +units=m +no_defs      When I put > 114.1689 22.4518 It outputs > 835447.18 834705.40 3.11 Can I tell from the result, the accuracy and possible numeric error of the result? Some meters or some percent?


Precision of x decimal points in meters?

Is there a handy reference that corelates the number of decimal points in a coordinate to the precision in meters? I had discussions like "4 decimal places will place you within a 500 meter radius" and "6 decimal places will place you in a house". I want a reference that can clear up such facts and requirements.


Can you suggest some well-written introductory texts about coordinate system projections?

Even after some time working in this GIS field, I still don't have a solid understanding of spatial references, how they are projected, and most importantly, **if you still need to worry about the fact that coordinates are geographic once they have been projected.** In which cases should one pay attention when handling coordinates in a projected coordinate system? For example, assume I am working in a projected coordinate system. I have a polygon whose vertices are now specified with X / Y coordinates. Let's say I want to manually calculate the polygon's area -- could I now take any normal algorithm (for Cartesian/Euclidean coordinate spaces) for computing the polygon's area, or do I still have to account for the fact that these are _geographic_ coordinates? What would I have to do to get an accurate result for my area computation? It's questions such as the above example that I still don't understand. Any cues to introductory texts explaining such details are welcome. (I was already considering looking at the PROJ.4 catrographic projection library and how it is used to deduce from that how to deal with projected coordinates, but that would've probably been too early.)


Publish Vector and Raster (GeoTiff) layer in Geoserver SLD

![enter image description here](http://i.stack.imgur.com/iJdMH.png) ![enter image description here](http://i.stack.imgur.com/b2VzS.png) Base layer and overlay layer are same when opening in ArcMap But When publishing the same layers in geoserver , distraction of 27 meter. I have attached an example image to show the problem.


Fixing address on TeleAtlas, how long before Google Maps has this information?

I recently submitted information to TeleAtlas and they told me it was accepted and has been added. When I search on Google Maps the address still does not appear. Does anyone know how long Google Maps take to update with new information or it just a matter of TeleAtlas making the data available for Google?


Calculate area with FME

I have coordinates in decimal degrees and need to calculate area for polygon objects. I tried to use AreaCalculator and set up criteria to find all these bigger than 100 sqm. But result is not as I expected. Should I convert first from decimal degree to metric?


QGis Save Raster as Rendered Image

In QGIS 1.9.0 Master, when you right click on a raster in the Layers Panel and select "Save As", you can select the output mode to be "Raw Data" or "Renedered Image". When selecting the Rendered Image option is saves the raster with the layers current styling. How would I do this via the python console? I am able to script gdal_translate, however I am not sure how to preserve the current styling for the layer. Any help would be much appreciated. Regards Gavin


Get Intersecting roads (ArcGIS Server10 with silverlight API)

I am developing a tool where when I click on a road in the map I have to list the roads intersecting the selected road using ArcGIS Server10 with silverlight API to develop this tool I got the roadID of the selected road using IdentifyQuery Then I got the envelope of that road using QueryTask I passed that road envelope to IdentifyTask again and queried the intersecting roads but its lists me all the roads inside the envelope but I need only intersecting roads also its very very slow :( I am stuck here Please guide Please guide


How to design a database for continents, countries, regions, cities and POIs?

I'm brand new to GIS programming and I am designing a GIS application. Target is to create system with continents, countries, regions (including states, sub-regions, provinces), cities and places in cities. Each of this elements will contain some text information and related stuff. As database we are going to use PostgreSQL with PostGIS. My question is how to design database for this system? I was thinking of 2 tables `polygons` and `points`, but I'm not sure if it's good way of thinking. What we need is some hierarchical base for relationship between countries, cities, regions, etc. Main goal of the application will be collecting tourist data from many sources about specific cities, regions, countries and so on, and store it in database. Let's say we will have city Rome, we collect some info about this city into database from couple sources. And we need to know that Rome is in province Rome, sub-region Lazio in region Lazio, country Italy. So system should be flexible to allow create such relation from real world.


Conversion of coordinates of "Bochum coordinate system"

I need to convert historic coordinate data into more recent Gauß-Krueger data. Specifically, I need to convert "Bochum coordinates" used till 1960 in the local coal mines into DHDN/Gauß-Krüger zone 2 (EPSG 31466, ID 2646). Is there in QGis a tool (similar to that that seems to exist is ESRI ArcGIS) that I can use for this job?


Workspace or data source is read only

I am using ModelBuilder in ArcCatalog 10.1 to convert an ASCII files to TIFF and add spatial reference. I have hundreds of files to convert so I'm using the Iterators. After doing 10-20 files successfully I get the following error.               Executing (ASCII to Raster): ASCIIToRaster "D:\Users\Documents\#Data\# New Data\LIDAR\todo\nz4626_dsm_1m.asc" "C:\Users\Data\# New Data\Lidar2.gdb\nz4626_dsm_1m" FLOAT     Start Time: Thu Aug 15 14:59:06 2013     ERROR 999999: Error executing function.     Workspace or data source is read only.     Workspace or data source is read only.     No spatial reference exists.     ERROR 010067: Error in executing grid expression.      Failed to execute (ASCII to Raster).     Failed at Thu Aug 15 14:59:06 2013 (Elapsed Time: 0.00 seconds)      After using command prompt to remove read only tag from the source and destination folders and restarting ArcCatalog the process can be run again for a few more files before failing again. What is causing the folders to become read only and how do I stop it?


Display 3d building using OpenLayers

Is there any way to display a 3d building model as a layer under OpenLayers ? i have the building models drawn using Google sketchup , is that useful ?


Divide map into features with invisible divisor lines

I'm working with QGIS 1.8.0. Is there a way to divide a map with a grid but make the divisor lines invisible? I have a shapefile as the one in the left and I want it to see as the one in the right but still having the divisions (with the 82 different features). Please let me know. Regards! ![actual images](http://i.stack.imgur.com/3N3l9.png)


Adding a scale using SDK GXP OpenGeo Suite

Followed the tutorial using the OpenGeo suite SDK, and looking at adding a scale to my maps however the tutorial doesn't cover adding widgets only plugins. So I am not sure about how I would add http://gxp.opengeo.org/master/doc/lib/widgets/ScaleOverlay.html wherever I place it fr example the tools section I seem to just get a white screen, which tells me something is wrong. I believe I should call it using `* @require widgets/ScaleOverlay.js` and then as with plugins something like (below) in the tools part of app.js                {     ptype: "gxp_ScaleOverlay",      actionTarget: "map.tbar"     },


ArcGIS in PhoneGap

I am trying to set up a PhoneGap application with an ArcGIS Map using a Windows Phone template in Visual Studio 2012. I have tried using an example from the ArcGIS Javascript API samples to incorporate a Mobile ArcGIS.com map in my application: http://help.arcgis.com/en/webapi/javascript/arcgis/jssamples/mobile_arcgis.html This has not worked. I am getting the following error in Visual Studio when I try to link to the external CSS files: > Unable to edit > 'http://serverapi.arcgisonline.com/jsapi/arcgis/3.4/js/esri/css/esri.css'. > The style sheet must be part of the current project and the project must be > a Web Site or Web Application project. I found a similar example of setting up a PhoneGap project for Android with an ArcGIS Map, but this has also not worked for me. Does anyone know of a working sample PhoneGap project to incorporate a map, using the ArcGIS Javascript API? Alternatively, does anyone know how I can get around the error with linking the stylesheet?


Why doesn't dissolve finish? It just stops without error message

I have a shapfile of the whole UK that contained Postcode area, district, and full postcode. I was going to try and dissolve the regions into seperate business areas by joining a lookup table that identified, for each postcode which business area it's assinged to. But i couldn't get the business area to appear in the dissolve popup window. So instead i modified the DBF file, adding the new business area field. Now when i go into the dissolve function, the dissolve field does appear there but i can't get it to succesfully finish. It just stops. No error messages. Nothing. I think what i'm trying to do is simple but QGIS just doesn't seem to be happy. This is the very first thing i'm trying to do with QGIS so please be patient with me.


How to make line with arrow head in QgsVectorLayer using python?

I need to draw line with arrow head or after drawn the line i need to make it as arrow head line by using python. or can we set any arrow head symbol for line features by using python? Kindly help me with some sample codes.


Changing field type of Excel attribute table for a Join

I have a string field type of shapefile already in ArcMap. I need to import an excel file in and join with the shapefile. The only problem is the field type of the shapefile is a String and the field type for the excel file is a Double. The excel field is all numbers, but the shapefile field is mostly numbers and a few letters mixed in, I think thats why it came out as a string. I know you can't add a field or use the field calculator into a table exported from excel, I'm currently using it as a .csv.


Why do I get a NameError: name 'gdbSuccessCount' is not defined?

I'm using a 412-line script that gathers and merges online tile data. However, I'm receiving a Traceback error, NameError: name 'gdbSuccessCount' is not defined. I've included the portions of the code that gdbSuccessCount is included in (only 3 times it appears as is). I, and my work, will appreciate the help to find the resolve the problem. I apologize for the code mess, the indented code is below it (first time on StackExchange). The problem is the last occurrence, in line 406 as follows: if merge and gdbSuccessCount > 1: Here are lines 333-412:                   r = 0         if GDB:             gdbSuccessCount = 0             for g in gdbList:                 tileGDBID = g                 q250k = "q" + string.split(tileGDBID,"-")[0]                 srcGDB = dataPath +"/" + q250k + "/geodatabase/"+ tileGDBID +".gdb.zip"                 destGDB = dest+"\\geodatabase\\"+ tileGDBID+".gdb"                 destGDBzip = dest+"\\geodatabase\\"+ tileGDBID+".gdb.zip"                 if not os.path.exists(destGDB):                     arcpy.AddMessage("Copying Geodatabase " + tileGDBID)                     if not os.path.exists(destGDBzip):                         try:                             urllib.urlretrieve(srcGDB, destGDBzip)                             gdbSuccessCount += 1                         except:                             arcpy.AddMessage("\t WARNING: Geodatabase for tile " + tileGDBID + " not found - moving on to next one...")                      if os.path.exists(destGDBzip):                         zf = zipfile.ZipFile(destGDBzip)                         zf.extractall(dest+"\\geodatabase")                         del zf                         os.remove(destGDBzip)                 else:                     arcpy.AddMessage("\tGeodatabase " + tileGDBID + " exists, moving on to next tile...")              if LAS:             # this is a batch file that will be used to unzip the LAZ files....             f1 = open(dest + "\\las\\unziplaz.bat", "w")             # now go through the list and download the LAZ files...             lasSuccessCount = 0             for l in lasList:                 q250k = "q" + string.split(l,"-")[0]                         srcLAZ = dataPath + "/" + q250k + "/laz/"+ l +".laz"                 destLAZ = dest + "\\las\\"+ l +".laz"                 destLAS = dest + "\\las\\"+ l + ".las"                 if not os.path.exists(destLAZ):                     arcpy.AddMessage("Copying LAZ " + destLAZ + "\n")                     try:                         urllib.urlretrieve(srcLAZ, destLAZ)                         lasSuccessCount += 1                     except:                         arcpy.AddMessage("\t WARNING: LAZ file for tile " + srcLAZ + " not found - moving on to next one...")                                 if not os.path.exists(destLAS) and os.path.exists(destLAZ):                     f1.write(dest + "\\las\\laszip.exe -v -i " + destLAZ + " -o " + destLAS + "\n")                  #close the batch file                 f1.close()             if lasSuccessCount > 0:                 # add a message telling the user what we are doing...                 arcpy.AddMessage("Unzipping LAZ Files....please be patient")                 # download LASZIP.EXE from the mngeo site...                 lazsrc = "ftp://ftp.lmic.state.mn.us/pub/data/elevation/lidar/tools/lastools/laszip.exe"                 lazdest = dest + "\\las\\laszip.exe"                 urllib.urlretrieve(lazsrc,lazdest)                      #call the batch file to unzip the LAZ files...                 os.system(dest + "\\las\\unziplaz.bat")                      # if the user has ArcGIS 10.1 then build a LAS dataset.                     if arcpy.GetInstallInfo()['Version'] == '10.2.2':                     lasDName = dest + "\\" + "las_data.lasd"                     if not os.path.exists(lasDName):                         arcpy.AddMessage("Creating LAS Dataset - " + lasDName)                         try:                             arcpy.CreateLasDataset_management(dest+"\\las",lasDName,"NO_RECURSION","#","#","COMPUTE_STATS","RELATIVE_PATHS")                         except:                             arcpy.AddWarning("LAS files downloaded but attempt to bring into LAS dataset failed...")                     else:                         arcpy.AddMessage("LAS Dataset - " + lasDName + " found....")             else:                 arcpy.AddWarning("No LAS files found in requested area....")              # if the user wants to merge the downloaded data then do that now....         if merge and gdbSuccessCount > 1:             arcpy.AddMessage("Merging tiles to "+ dest + "\\ELEVATION_DATA.GDB")             merge_tiles(dest)         elif gdbCount == 1:             arcpy.AddMessage("Only one GDB found...no merging necessary")         elif gdbCount == 0:             arcpy.AddWarning("No Geodatabase Tiles found in area requested....")


Joining a table with thousands of transactions to one point (location) in Arc Map

I have a large transaction data set in the millions and around a 100 positional points. I need to join the transaction data set to these points. Arc map won't let me join them without errors or if it does allow it null values are applied to the join. Any help would be appreciated. They have a common field, a meter number for the machines. I have tried using a simple join with the attribute tables. My end goal is to be able to upload this data to an enterprise and integrate it into an operations dashboard. This dashboard would be able to tell the user the number of transactions with the machine within a time period. The transaction data set also has the date attached to it. ArcGIs-10.0 _update_ I can't comment so this will have to do. I found a problem the one web service map I am using has the meters as strings so in the transactions attribute table I added a new field and made it a text. Using the field calculator I added the meter number I was then finally able to perform the join, but only 153 values are in this join meaning it didn't join everything I wanted. How can I get all of the values to join?


only 1 tile visible in GE after gdal_translate -of KMLSUPEROVERLAY

I have a 3-band ASTER scene in any of several formats (PCIDISK, GeoTIFF, Imagine), already reprojected to the geographic coordinate system. I can view and manipulate the scene successfully in PCI Geomatica and ESRI ArcMap. I want to create a KMZ version of it, and following in part this other GIS StackExchange post, I tried using:               gdal_translate -of KMLSUPEROVERLAY myscene.tif myscene.kmz      but the result is invisible in Google Earth (although it zooms to the correct area). It remains invisible when zooming way out and way in, except that the lower-left tile of the 16 tiles it created **does** show up when zooming in. It, however, is completely black as it is in the no-data collar around the image data. Also, if I unzip the KMZ and open any of the individual tile KMLs that were created, **they** open just fine, and the ones that are not in the collar show the correct image data. So I suspect there's something wrong in the way the output KMZ was formed (which uses Regions and NetworkLinks, KML elements I have not worked with before). But I don't know whether this is due to input data or format problems, to incorrect or missing parameters to gdal_translate, or something else. I've tried several input formats with the same results each time. Tips and suggestions (if not outright answers!) much appreciated!


Is there any way to select only the layer of which the features have been selected? That is, if I have selected some roads then I want only the road layer to be selected.


How to fill Field Mappings parameter in script tool on updateParameters?

My goal is to get a list of fields choosen by an user from existing fields in feature class. I think that I can show all fields in Field Mappings dialog, an user deletes unwanted fields and I get a required list.   The problem is how to populate Field Mappings when adding a feature class as a first parameter. I can't select anything in 'Obtained from' property.   The code `self.params[1].parameterDependencies = [0]` in updateParameters function does't work in this case (it worked when 'Field' parameter data type was chosen). What should I use to automatically populate Field Mappings when filling first parameter? * * * EDIT:   Fields list is applied as follows: ![enter image description here](http://i.stack.imgur.com/HGDsE.png)


calculate building density in ArcMap

I have two questions related to calculate building density.   1. I need to calculate the building density in a map which Data is in shape format. The buildings and the area are in two different layers.   2. Same question as above but the map/area is raster format. Thank you.


Land Administration Domain Model example implementations

Does anybody know of a sample implementation? I am currently trying to adapt it to my country but the LADM is rather abstract. Anybody have any idea how to go about this?


How to perform volume calculation in QGIS

EDIT 07/12/12 Ok, meanwhile I came to another problem I want to discuss. I just tried around with the volume calculation with some layers and came to a point where I receive the following error. Screenshot: ![error](http://i.stack.imgur.com/XTQNg.jpg) I did some research and discovered this explanation: http://blog.webforefront.com/archives/2011/02/python_ascii_co.html The problem is, that I have absolutely no Idea about python programming. So my question is: Does anyone know how to solve this described error problem? Note: I left the older questions below in order to see the process I made. * * * * * * * * * I m currently trying to make a volume calculation in QGIS. I realised, that using GRASS or Sextante is the right choice to do this. But since there are barely useful tutorials, that go into depht nor any other guides or hints that may help in this context, I want to seek help here. So I want to ask these questions: 1)I want to calculate the volume of a specific map area. So I discovered the algorithm "r.volume." 1.1.)How does it work exactly?The Help text in that menue is not very useful for me as a Grass beginner. 1.2) Why do I need raster layers for this and no vector layers? 1.3) What certain criteria do these rasters have to fulfill? 1.4) I added a screenshot. What exactly are these "clumps" and what rasters have to be chosen in the drop down list? 1.5) Can I have an accurate example of how to perform such an r.volume algorithm?I d be thankful for any detailed help/ raster files to test this. 2)I learned, that GRASS (in QGIS) was integrated into the Sextante Plugin for QGIS. 2.1) Does this mean I can only use Sextante and still have all the functions of Grassgis?I d welcome that because I always had trouble with importing data into Grassgis. Sextante doesn t (obiously) require this.(?) 2.2) A strengh of **gvSIG** is Sextante, so I read; Is the Sextante Plugin for QGIS equal to gvSIG Sextante concerning performance/ possibilities? thanks a lot, Shepherdjo * * * EDIT:07/11/12 8.08am: I know this is not easy to answer but it would be great too, if you can help me out with the other questions (2.1/2.2). * * * ![r.volume menue](http://i.stack.imgur.com/XUIH7.png)


Error when importing shapefile to Postgis database

I was trying to import shp to postgis database and i get the following error:               Importing with configuration: adresy92, public, geom, C:\Users\tomek\Desktop\adresy92, mode=c, dump=1, simple=0, geography=0, index=1, shape=1, srid=0     Shapefile type: PointZ     PostGIS type: POINT[4]     Failed SQL begins: "SET CLIENT_ENCODING TO UTF8;     SET STANDARD_CONFORMING_STRINGS TO ON;     BEGIN;     CREATE TABLE "public"."adresy92" (gid serial,     "obrot" int2,     "dane_kod" int2,     "xy_dane" varchar(20),     "status_kod" int2,     "xy_status" varchar(50),     "ulica_kod" int2,     "miejsc_kod" int2"     Failed in pgui_exec(): ERROR:  AddGeometryColumns() - invalid SRID     CONTEXT:  SQL statement "SELECT AddGeometryColumn('', $1 , $2 , $3 , $4 , $5 , $6 )"     PL/pgSQL function "addgeometrycolumn" line 4 at SQL statement     Shapefile import failed.      The stange thing is that I have already put a couple shapefiles to this database without any problems. I am importing shapefile through SPIT tool from Quantum GIS 1.8.0. I have also tried to do the same via shp2pgsql plugin in pgAdmin with the same result. PS I found a similar question on gis.stackexchange but it didn't solve my problem. I would be appreciate for your help.


raster algebra in python with rasters of different extents

I am trying to find out how to use GDAL & numpy modules to average a set of rasters....which are Sigma0 values from satellite passes at different times within a year. Each raster has been mosaicked from a number of smaller images using GDAL_merge. Because the orbits differ each time they pass over the area of interest the merged datasets are different extents and not regular squares. So I want to average the values from each pass....baring in mind that when theoretically stacked on top of one another, there is sometimes overlaps of one/two/three images. I imagine the best way of me getting around this is making all rasters the same extent ... and using a 'no data' value for areas of the raster where there is no data...and then ignoring these values during the calc of the average... If this indeed the best way, how do I go about making them all the same extent when they are not regular (rectangle)? ...or is there a better way than this? I am new to using GDAL/numpy ....my impression is that calculations using multiple rasters is best done in numpy arrays? again, is this correct. Thank you in advance Becky


How to avoid colliding parallel lines in QGIS Composer

I am preparing a map which is concerned with showing two line features and there relative location to each other. The line features are 20 m apart and at a scale of 1:50 000 or bigger on A3 I want to show which the lines a little separated from each other. The features follow a common path, but change sides irregularly, i.e. one feature is on the south for while, then the north, then th common path changes north south and the line features then alternate between sitting on the east or west of the path. I would like to do this in QGIS. Is there some sort of automagical dynamic offset control for lines? I can imagine how to rotate symbols based on an adjacent feature, but this is a bit more complicated. I can imagine a method of breaking one line up into segments and then creating an index of those segments compared to the other line acting as an anchor. The index might then be used to specify the offset, however this is a sort one time run solution. Rather just looking for the elegant cartographic solutions, if they exist. ta


Get XY coordinates in UTM from web-mercator where UTM zone is not known

How to get XY(northing/easting) coordinates in UTM where data is WGS_1984_Web_Mercator_Auxiliary_Sphere(WKID: 3857) coordinate system and UTM zone is unknown. For e.g. India falls under multiple UTM zones and map is in Web_Mercator_Auxiliary_Sphere. I am trying to replicate Google earth's UTM coordinate display using ArcGIS Javascript API. Thanks. Nagesh


How to get a collection of georeferencing control points in an ArcMap document?

I have created control points (i.e. links) using ArcMap's Georeferencing toolbar. How can I programmatically retrieve a list of these points? I'm looking for something like an enum of the links you would see in the Link Table (which is visible through the Georeferencing toolbar). Preferably, I would then be able to step through these points and delete any one. I am working with the ArcObjects SDK for C# .NET.


choosing the right parameters to interpolate a DEM form contours

this is my first question on the forums, so be gentle :) I have a 5 meter interval contour lines feature class which i would like to interpolate to a DEM. I am well aware of at least 2 tools ("tin to raster" and "topo to raster") which can help me with that. In order to produce the proper DEM i need to define first its purpose. Long story short - i need to track the progression\advancement of small streams (AKA: gullies). they are typically 1 meter wide and 2 m long. My algorithm which will be using the DEM as a basis, must have as much as accurate DEM as possible. how accurate? - detailed enough to see an advancement of a 30 cm length for example. To do that i must interpolate (from the mentioned above contours) a most pleasing DEM. The questions is: what course of action would you suggest me to implement? in order to choose the correct parameters? I know trial and error would be the most advisable answer, yet there are so many parameters (smoothing, cell size, method type, enforce drainage, toelrance) to consider it might take me weeks until i get the correct DEM. Cheers, David


Cutting image using gdalwarp to extent of data

I am currently generating a raster from point data and applying inverse distance algorithm. This image is then cut with :               gdalwarp -cutline <<src shape file>> -crop_to_cutline -srcnodata 0 results_in.tiff results_out.tiff      Is there a way to cut the output to the extremities of the original data source file rather than a shape file? Many thanks J EDIT : gdal_grid command gdal_grid -a invdist:power=2:smoothing=0.001:radius1=0.0005:radius2=0.0005 -of GTiff -zField z_value -ot Float64 -l data VRT.vrt EM_RESULTS.tiff


Is it possible to add a gvSIG document on Html?And if yes, how?

I have completed a project in gvSIG and now I would like to be able to open it from my webpage, like a hyperlink or path or something.Can I really do that somehow? And if yes, could you please explain how? Thanks a lot. Eva


Looking for average temperature data for Afghanistan

Hey all I am looking for a average temperature data for Afghanistan monthly over a year, precipitation would work too. In a shapefile, if anyone has any or knows anywhere that I could access one please let me know.


What format to use for highly compressed vector data?

Does anyone have information about possible formats for highly compressed vector data? Some basic requirements.   * Must be possible to read using the larger commercially available platforms. (ArcGIS, Luciad, MapInfo etc.)   * Must have tools that allows us to create files/format from other data formats (i.e. Shape)   * Speed when reading the format is important as the data is among other things used as a live map in mobile platforms. Compression speed is less important (although nice)   * Must be possible to use a map product in a wide scale range, a 50k map shall be possible to read at a 1 000 k zoom level   * Ability to do selective updates on area, feature-ID or other search criteria would be nice, but not a must.   * Memory footprint must be small, it needs to work with Windows XP and/or Windows 7 allocating less than 1GB RAM so some kind of smart cache management must be implemented To give an idea of the required compression and as input for memory footprint requirements, I need to store approximately 1TB of shape data to a 100GB disk. Does anyone have suggestions for good products? Best regards, Johan


How to summarize population by neighborhood?

I'm a novice with PostGIS and am currently working on a lab assignment where I have two geometric tables - one containing demographic data by blockgroup (blockgrp2010) while the other identifies various neighborhoods (nbo_hood). I'm trying to write a query that lists the population for each neighborhood. For now, I'm assuming that if the centroid of the blockgroup is inside the neighborhood, then all the population lies within that neighborhood. I drafted the query below:               select nbo_hood.name, blockgrp2010.pop10     from blockgrp2010     inner join nbo_hood on ST_Within ((ST_Centroid (blockgrp2010.geom)), nbo_hood.geom)     group by nbo_hood.name, nbo_hood.geom, blockgrp2010.pop10;      This returns a table listing the neighborhoods and corresponding population for each blockgroup. However, many of the neighborhoods are listed under multiple rows with a different population assigned (obviously because there are multiple blockgroups within the neighborhood). My question: how do I modify the above query (or should I start from scratch) to list each neighborhood only once with the population sum of all blockgroups? I've tried variations of the SUM function with the SELECT statement but nothing as worked. Thank you,


Is co kriging appropriate for interpolation of three spatial co variates?

I want to estimate an interpolation model(s) for daily minimum temperature, daily maximum temperature and daily rainfall. I have samples of approximately the same size for all three variables at different locations although for certain days and locations only 1 or 2 of the 3 covariates are sampled . All three variables have some correlation, especially minimum temperature and maximum temperature. I read that co-kriging is especially adequate in cases where:   * one wishes to interpolate one variable that is relatively sparsely sampled   * one has available another covariate which is more densely sampled This is not my case. Rather I have available samples for three covariates and want to interpolate all three. I have the following questions:   1. Is it better to develop a separate ordinary/universal kriging model for each variable or to develop a co-kriging model for all three.   2. I am using days and locations with no missing values to estimate the parameters of the interpolation model , but want to use it to interpolate missing values in different days and locations. Is it possible to do this in cases where:     * certain locations have all three variables missing     * certain locations have one o two out of the three variables missing I am working in R in case there are any special considerations to heed.


How to reset ObjectID after running Delete_Features on SDE feature table?

I have a python script that deletes all features on a SDE feature table and then inserts new features to it. I noted that the ObjectID keeps incrementing from where the last feature was with before deleted. Just wondering if there is any reset function or tool to make ObjectID start with 1 every time the script runs? Thanks.


OpenLayers.Filter.Logical not working

I have the following code:                  var map = {};        var googleMercator = new OpenLayers.Projection("EPSG:900913");        var wgs84 = new OpenLayers.Projection("EPSG:4326");        var lyr_osm = {};        var texas = {};        var lll = {};        var query = {};        var fin = {};        var fout = {};     function init(){         OpenLayers.ProxyHost = "proxy.cgi?url=";     //--------------------------------------------------map         $(document).ready(function () {         map = new OpenLayers.Map('map',                         {                         numZoomLevels: 18,                         projection: googleMercator                         }            );         });     //--------------------------------------------------layers definition         var lyr_osm =new OpenLayers.Layer.OSM("OSM");     //--------------------------------------------------setting wms layer            texas = new OpenLayers.Layer.WMS("Texas",                "http://54.228.242.237:80/geoserver/us_admin/wms",                {                'VERSION' : '1.1.0',                layers: "us_states_gen",                format: "image/png",                transparent: true,                cql_filter: "state_name LIKE 'Texas'",                },                     {isBaseLayer: false,                     visibility: true,                     sphericalMercator: true,                                             }                );     //---------------------------------------------------filters                                  fin = new OpenLayers.Filter.Comparison({                                          type: OpenLayers.Filter.Comparison.GREATER_THAN,                                          property: "store_id",                                          value: 500                                  });                                  fout = new OpenLayers.Filter.Comparison({                                          type: OpenLayers.Filter.Comparison.LESS_THAN,                                          property: "store_id",                                          value: 100                                  });                                  query = new OpenLayers.Filter.Logical({                                         type: OpenLayers.Filter.Logical.AND,                                         filters: [fin, fout]                                  });              var awfspc = new OpenLayers.Protocol.WFS({             version: "1.1.0",             url: 'http://loclahost/geoserver/stores/ows',             featurePrefix: 'stores',             featureType: 'stores',                 defaultFilter: fin,                 srsName: "EPSG:3857"                 });              lll = new OpenLayers.Layer.Vector("Points", {                  externalGraphic: 'http://localhost/black_marker32.png',                   graphicHeight: 20,                   graphicWidth: 20,               strategies: [new OpenLayers.Strategy.Refresh(), new OpenLayers.Strategy.Fixed()],               protocol : awfspc,         });          //--------------------------------------------------------------------map controls         map.addControl(new OpenLayers.Control.LayerSwitcher());         map.addControl(new OpenLayers.Control.MousePosition());     //--------------------------------------------------------------------add Layers         map.addLayers([lyr_osm, texas, lll]);     //--------------------------------------------------------------------map position         var mapCenter = new OpenLayers.LonLat(-99, 31.3);         mapCenter.transform(wgs84, googleMercator);         map.setCenter(mapCenter,6);     //  map.zoomToExtent();     }      When I use fin filter or fout filter, the filters are applied and they return good. When I use querry filter, nothing is returned, no error and no data. Any clue? Thanks!


Maplex: try different fields

Let's say, a road has the name "William-Bill-Gates-Road". If the road is too short, the label can only be abbreviated by using a dictionary in Maplex. Is it possible to make Maplex try several label fields, if one is not fitting? Structure: name | name_1 | name_2 | ... William-Bill-Gates-Road | William-Bill-Gates-Rd. | W.-Bill-Gates-Rd. | W.-B.-Gates-Rd.


Ecuador Census Data

I am looking for Ecuador census data (digital GIS if possible) at the sub- parroquia (parish) level. The geographies and the tabular data would be great but at this point I am mostly in need of the digital boundaries. I am looking for the full demographic suite of data with attributes if possible. 2010 would be great but I am also interested in 2001 and 1990. Any digital GIS format will suffice for the boundaries and any tabular format for the tabular data. Non-free options are fine as well. The geographies are known as Sector and Zona (Zone) in the Ecuadorian Census.


Calculate end point Height above the MSL of the target point given a starting point, distance, azimuth and elevation

We need to find the end point height (from the earth surface) given the start point lat, lon, azimuth (270 deg), distance (500 km) and elevation above the horizon of start point (10 deg). The question **How To Calculate Height of Point B?** does not seem to consider the earth's curvature:   * Here the angle (elevation) at point a is 10 deg    * so the calculation referred in the above question is only usual cartesian calculation (i.e h).  I would like to know the mathematical formula for the geodetic calculation of the height POINT B above the MSL at POINT B (i.e h +H) in the picture shown below. Also I would like to program the formula using VB.NET on a independent GIS application. Accuracy in the order of meter would be good enough. Please provide the mathematical formula to calculate the height above earth surface at the end point? ![enter image description here](http://i.stack.imgur.com/2X7kO.jpg)


How to debug QGIS installation problems on XP?

I'm installing Qgis (Standalone 1.8) on all the PCs of the town planning office where I work. All good for 4 computers, on the fifth installation does not work. Computers have all XP with SP3. In particular: the installation seems to proceed normally until the end, then unlike the other pc, ends without the window appears as Microsoft Visual C++ and bypasses all the post-install. After installation, the desktop icons are missing Grass, MSYS and OSGeo4W. There are Qgis Qgis desktop and browser, but will not start. The same problem occurs with the previous version. I can not figure out what is missing on this computer, which is mine. Is there anyone who has had the same problem and fixed it?


How To Do OpenLayers.Mobile Popup

It seems that the OpenLayers.mobile.js library (v 2.12) doesn't have a popup class. Does anyone know how to do a popup with the mobile version? I'm using some code that i ported over form normal OL web map. When i call new OpenLayers.Popup.Anchored i get an undefined error on OpenLayers.Popup class. It's not in the OpenLayers.mobile.js. I'm lost


Creating a custom, local Route Task

I'm trying to replicate something similar to the ArcGIS WPF Runtime SDK (C:\ArcGIS\WPF1.0\SDK\Samples\Sample Application\Sample Application.exe) "Routing with Barriers Online" sample with two major modifications ... 1. I have my own network dataset (created from a feature class using ArcObjects classes) 2. I want it to run locally and not on an ArcGIS server instance. Can I exploit LocalServer for this? First, can this be done? If so, what is the best approach? I have both the WPF runtime and ArcObjects SDK to work with.


merge shapefile using ogr2ogr and retrieve the attribute of both

Ciao, I have two shape files that I want merge (overlay/intersect) having in the destination the attributes of the two sources file.               ogr2ogr file_merged.shp tl_2010_us_uac10_conus_albersUA.shp     ogr2ogr -update -append file_merged.shp tabblock2010_10_pophu.shp -nln file_merged      In this way I obtain the polygons of the two shapes, but I lose the attributes of the `tabblock2010_10_pophu.shp`. My final aim is then make an sql statement where I can retrieve polygon in common in the both shapes. Any suggestions?? Thanks in advance


ArcGIS 10: convert labels to points

is there a way to convert labels, created by Maplex labeling engine to point features? In MapInfo and CAD software, it's possible to convert labels to points. Those point features include a text and a rotation for this point. When text is horizontal, it includes all the text. When the text is curved, it included one char for each point with the rotation in degrees. That makes it possible to use the labels in other software or store them in a database. ArcGIS 10 only allows to convert labels to annotations, which result in polygons... Thanks for help :)


Interacting with the Map after selecting a button/tool

I am using the ESRI Javascript API 2.0 and Dojo dijits. Does anyone know how to wire up a dojo button (or toggle button) in a toolbar to a mouse click event on the map? I want the user to have to select a tool/button before clicking on the map and NOT just have the map listen to any old mouse click event as this is too open i.e. I want to constrain the functionality to the selecting a tool from the UI. Many thanks


Detect whether there is a spatial polygon in a spatial extent

Given a spatial extent (for example `cbind(c(0,0),c(10,10)`),I want to see whether there is a part of a SpatialPolygons falls in this extent? How can I do this? Thank you very much!


Python 3 with QGIS on Mac?

I need to start learning Python, and I want to learn with version 3.x. I understand that QGIS--at least some plug-ins--relies on Python and Tcl/Tk, and probably on specific versions. Assuming that QGIS is dependent on Python 2, can anyone tell me how to set up parallel installations of Python 3 and Tcl/Tk such that QGIS will have what it needs and I can use Python 3 for other purposes?


Converting Coverage files into shapefiles?

I am attempting to convert coverage files into shapefiles using arcpy. This is not the e00 files, but the bin files. ArcGIS help suggests that it can be done using **FeatureClassToFeatureClass_conversion** , but I am getting an error saying "parcels does not exist or not supported".               fileName = 'C:\\May_2014_Parcels' + "\\" + filename + "\\" + subfilename + "\\parcels"                        arcpy.FeatureClassToFeatureClass_conversion(fileName,            r'C:\May_2014_Parcels\Output\parcels.gdb', filename + "_" + subfilename)      How do I actually get the featureclass out of the coverage file so I can perform this task?


Show M values on mappoint /mouseover on a polyline - ArcGIS Server 10

I want to show the M values i.e the chainage of the road when I move my mouse over it. I am using mappoint to show some static attributes but I also want to show the dynamic M values along the road when I mouse my mouse pointer over it.. Note: I am using Visual Studio 2010 with Silverlight API 2.2 and ArcGIS Silverlight 10 Please guide...


how to filter point features based on the filterted polygons with OpenLayers

I have one WFS polygon layer and one WFS point layer. The polygon layer has an attribute user_id. The web map has two requirements.   1. Filter the polygons based on the attribute user_id. That is, a user can only see polygons which are assigned to him. This can be done like this: filter: new OpenLayers.Filter.Comparison({ type : OpenLayers.Filter.Comparison.EQUAL_TO, property : "user_id", value : user_id })   2. Filter the points based on the filtered polygons. That is, a user can only see points within the filtered polygons from step1. I checked the  OpenLayers.Filter.Spatial.INTERSECTS. it only takes OpenLayers.Bounds || OpenLayers.Geometry. I have no idea how to achieve this. Please help. Thanks.


How to return the first row of feature class as a list in ArcPy?

I need to return the first record of a feature class as a list. Is this possible? I can't seem to figure it out using `row.getValue()` as it does not allow me to put a wildcard in there.


How to get route plan with certain distance with pgRouting?

I am working for a project. I have a road network which has been imported in PostgreSQL/PostGIS database. I can query shortest path between two points. I want to get route plan with certain distance. Suppose, I want to go from point A to point B. There may have different way to go from point A to B. But, I want to fix the length, suppose, 10 km. That means, I want to know the route plan to go from point A to point B but in 10 km range. Is it possible to get this kind of query using pgRouting? Or I should use any other tools? If it is possible how should I write the query? Note: in my database table 'route', i have 'gid', 'source', 'target', 'cost' etc columns. Thanks, in advance.


Vehicle Routing Problem - Cost Calculation According to the furthest point in route

Certain logistics companies do not calculate their costs according to the total distance travelled but according to the furthest point in the route. In our problem, Vehicle Routing Problem do not really help as almost all of the cases, furthest point is always have less travel distance than the whole route. Do you have any idea on how to tweak Network Analyst in ArcGIS, or know any open-source solution that I can solve this problem with? Thanks a lot in advance, Anil. EDIT:![enter image description here](http://i.stack.imgur.com/mf8rc.png) As seen on the graphic above, it really does not matter how many points you stopped between start and end point. However, we still would have other constraints, such as travel time, vehicle capacity, multiple routes and depots.


Field Collection Apps for Windows 8

I was wondering if any of you out there have experience with any apps that work on windows 8 tablets free or paid?We have a client that wants us to use their tablets for data collection. i was wondering if some of you out there have used apps for the windows store that you can recommend. Thank you


QGIS style save issues

I have a map that I want to do graduated chloropleth maps for each year over 10 years, and I would like to use the same graduations for each year. I know this was asked earlier and therefore I know that I save it as a style and then load it for each new year, but I'm having problems with the save style/load style function in qgis - every time I load the style the data converts back to the original year. Has anyone else had this problem and figured out a solution? I'm using qgis 2.2.0 valmiera (desktop).


Geoserver shutdown after load big imagemosaic

catalina.log: 24.05.2013 15:33:38 org.apache.tomcat.util.net.JIoEndpoint processSocket SEVERE: Error allocating socket processor java.lang.OutOfMemoryError: unable to create new native thread at java.lang.Thread.start0(Native Method) at java.lang.Thread.start(Thread.java:657) at org.apache.tomcat.util.net.JIoEndpoint$Worker.start(JIoEndpoint.java:513) at org.apache.tomcat.util.net.JIoEndpoint.newWorkerThread(JIoEndpoint.java:744) at org.apache.tomcat.util.net.JIoEndpoint.createWorkerThread(JIoEndpoint.java:723) at org.apache.tomcat.util.net.JIoEndpoint.getWorkerThread(JIoEndpoint.java:757) at org.apache.tomcat.util.net.JIoEndpoint.processSocket(JIoEndpoint.java:789) at org.apache.tomcat.util.net.JIoEndpoint$Acceptor.run(JIoEndpoint.java:355) at java.lang.Thread.run(Thread.java:679) ...................... and Caused by: java.lang.NoClassDefFoundError: Could not initialize class sun.font.FontManager at sun.java2d.SunGraphicsEnvironment$2.run(SunGraphicsEnvironment.java:179) at java.security.AccessController.doPrivileged(Native Method) at sun.java2d.SunGraphicsEnvironment.(SunGraphicsEnvironment.java:163) at sun.awt.X11GraphicsEnvironment.(X11GraphicsEnvironment.java:253) at sun.reflect.GeneratedConstructorAccessor109.newInstance(Unknown Source) at sun.reflect.DelegatingConstructorAccessorImpl.newInstance(DelegatingConstructorAccessorImpl.java:45) at java.lang.reflect.Constructor.newInstance(Constructor.java:532) at java.lang.Class.newInstance0(Class.java:374) at java.lang.Class.newInstance(Class.java:327) at java.awt.GraphicsEnvironment.getLocalGraphicsEnvironment(GraphicsEnvironment.java:82) at java.awt.image.BufferedImage.createGraphics(BufferedImage.java:1152) at javax.media.jai.TiledImageGraphics.getBogusGraphics2D(TiledImageGraphics.java:266) at javax.media.jai.TiledImageGraphics.(TiledImageGraphics.java:192) at javax.media.jai.TiledImage.createGraphics(TiledImage.java:755) at javax.media.jai.ROIShape.getAsImage(ROIShape.java:1148) at com.sun.media.jai.opimage.MosaicOpImage.(MosaicOpImage.java:374) at com.sun.media.jai.mlib.MlibMosaicOpImage.(MlibMosaicOpImage.java:49) at com.sun.media.jai.mlib.MlibMosaicRIF.create(MlibMosaicRIF.java:74) ................... and more ERROR [geotools.map] - Call MapContent dispose() to prevent memory leaks .............. After run layer preview(big imagemosaic) geoserver(and other opengeo) send error 502.