ArcGIS Android - Get LatLon position of the graphic

I'm adding marker to my map using the following code:               mvHelper.addMarkerGraphic(Double.parseDouble(friendData.getLatitude()), Double.parseDouble(friendData.getLongitude()), friendData.getFriendName(), friendData.getStatus(), null, icon, false,                     0);      This is my listener:               mvHelper.setOnGraphicClickListener(new OnGraphicClickListener() {             @Override             public void onGraphicClick(Graphic graphic) {                 Log.d("MVHELPER GRAPHIC LISTENER TAG",graphic.toString());                 calculateDirection((float)((Point)graphic.getGeometry()).getX(),(float)((Point)graphic.getGeometry()).getY());             }         });        I'm trying to get the lat long position of the clicked marker. But the listener doesn't provide it. So I tried the above way and got this error: `com.esri.core.io.EsriServiceException:Unable to complete operation. Location "Location 2" in "Stops" is unlocated. Need atleast 2 valid stops. "Stops" does not contain valid input for any route.` and this is the method where I get the route:               calculateDirection(float x, float y){     Log.d("Marker X and Y",String.format("%f // %f",x,y));     final Point loc = mapView.toMapPoint(x, y);     Log.d("Marker X and Y in MapPoint",String.format("%f // %f",loc.getX(),loc.getY()));     Thread t = new Thread() {             @Override             public void run() {                 try {                          RouteParameters rp = mRouteTask                             .retrieveDefaultRouteTaskParameters();                     NAFeaturesAsFeature rfaf = new NAFeaturesAsFeature();                          Point p = (Point) GeometryEngine.project(loc, wm,                             egs);                     Log.d("Marker X and Y after projection",String.format("%f // %f",p.getX(),p.getY()));                     StopGraphic point1 = new StopGraphic(mLocation);                     StopGraphic point2 = new StopGraphic(p);                     rfaf.setFeatures(new Graphic[] { point1, point2 });                     rfaf.setCompressedRequest(true);                     rp.setStops(rfaf);                          rp.setOutSpatialReference(wm);                          mResults = mRouteTask.solve(rp);                     mHandler.post(mUpdateResults);                 } catch (Exception e) {                     mException = e;                     mHandler.post(mUpdateResults);                 }             }         };              t.start();      } **Update 1:** After trying @jdONeill answer. I get the Lat and Long values.This is what I did.               Point temp = (Point)graphic.getGeometry();          Point p = (Point) GeometryEngine.project(temp, wm,egs);     calculateDirection((float) p.getX(), (float) p.getY());      This is the result:               04-30 16:17:54.408  17908-17908/com.test.esrisample D/Marker X and Y﹕ -93.203285 // 44.881104     04-30 16:17:54.408  17908-17908/com.test.esrisample D/Marker X and Y in MapPoint﹕ -10379311.931913 // 5607101.545301     04-30 16:17:54.458  17908-17908/com.test.esrisample D/dalvikvm﹕ GC_FOR_ALLOC freed 484K, 22% free 10018K/12760K, paused 16ms, total 16ms     04-30 16:17:54.719  17908-18235/com.test.esrisample D/Marker X and Y after projection﹕ -93.238945 // 44.908331      Why is there a slight change in the latitude and longitude values?


Can I create a polyline with one long segment and three shorter underlying segments?

We currently have a street centerline polyline feature class where the street segments correspond to our pavement management database. Each segment has a unique ID that matches a record tracking pavement type/year. Segment breaks occur at intersections, but can also occur mid-block if a subdivision was built in multiple additions. We would like to start adding Left From/To and Right From/To fields to the street centerline feature class for geocoding purposes. I would like to be able to have a street segment start at an intersection and end at an intersection for address coding, but still keep the underlying multiple segments that correspond to the breaks in the pavement management database. Is there an easy way to do this? Do I need to have two seperate line files or is there way to have this in one feature class?


Calculate distance between X/Y coordinates by individual in R

In R, I'd like to calculate the straight-line distance between each capture of an animal in my dataset - no need to take into account the curvature of the earth. I have a dataset that looks like this:               record.nbr<- c(1:8)     individual<-  c(1,2,1,2,1,3,1,2)     x<-c(167685,167945,167685,153985,167685,158675,167645,167667)     y<-c(9876548,9879248,9876838,9596548,9926548,9878578,9876548,9166548)     julian.date<-c(125,126,127,127,128,129,130,130)     captures<-data.frame(cbind(record.nbr, individual, x, y, julian.date))      Here each "individual" represents a unique animal, and x/y points represent locations where the animals were captured. x/y points are UTM coordinates. Some individuals are caught many times, while other individuals are only captured once. I'd like to calculate the distance (and the mean date) between recaptures by individual. In the output table, each row would summarize the distance between consecutive recaptures of an individual, plus the mean date between recaptures, like so:               individual first.rec  second.rec  distance  mean.date     1          1             3        V         126     1          3             5        W         127.5     1          5             7        X         129     2          2             4        Y         126.5     3          4             6        Z         128


How to use Scale tool when Georeferencing Raster in ArcMap?

I'm georeferencing a bunch of maps in ArcMap 10.0. The Shift and Rotate commands are pretty easy to use, no problems there. But whenever I use the Scale command, the image completely changes position. Is there a way to put an anchor down (similar to the rotate tool)? I feel like I must be missing something but haven't figured out what it is yet. I don't want to rubbersheet them into place, so I'm trying not to use control points until I have it in position, but first I need to make sure the scale bar on the map matches the distance it's supposed to. Edit: Asking a question does wonders in finding an answer, even when nobody has answered yet. I found the Rescale_management tool and will use that to set the scale properly first, then position. If the length of the scale bar is 1256ft and it should be 1200ft, I'll use a factor of .9554 in the X and the Y directions. Is there a better way?


How to assign node IDs to links in a network?

I have two shapefile layers in ArcGIS for Desktop 10. I have one line layer and node layer. My node layer(shp) has a "code_1" field and line layer has a "code_node" field. The intersection point of the line layer has only one node. I want get the code from the "code_1" field of the node layer for the line's "code_node" field For example, for one line, the first node has code_1=1 and the node at the end of line has code_1=2. In this case the `code_node` attribute in that line must be "12",. How can I automatically get the code from the nodes for each line? ![enter image description here](http://i.stack.imgur.com/vVegb.gif)


Can i access UTF grids without a server?

I'm creating a slippery map that will run off of a usb stick or cd (client request). I'm able to load up the map, display pre-cooked tiles and some some simple vectors. The client would like some more interactivity and the level of detail they want (cell values from rasters) is far to large for json files. I've made UTF grids using tilestache and have created a directory for them in the same folder as the map html. However, the UTF grids wont display and i suspect it is because it violates the same-origin policy, so my questions are as follows:   1. Can i access and dispaly UTF grids in local storage without a server?   2. Is there another route i should be taking other than UTF grids? Thank you


What is the difference between Sub and a Function in Mapbasic?

In MapBasic it seems like Functions can have side effects so I don't really see the difference between using the following: `Declare Sub MyProcedure (myParameter As String)` and then use `myParameter` after `Call MyProcedure(myParameter)` or `Declare Function MyProcedure () As string` and do something like `myParameter = MyProcedure()` Am I missing something?


I am trying to set up Geoportal to harvest services hosted by Geoserver. I tried follow this instructions:   1. Open your Geoportal Server.    2. Go to the Administration tab Choose the 'Add' option    3. Choose 'Register resource on the network' and click Proceed    4. In the URL resource form, paste the OGC GetCapabilities URL (including all parameters)    5. Keep the options for allowing the resource to be found and for it to be synchronized at an interval    6. Choose the interval.    7. Select if you want to be notified of synchronization results.    8. Click Save.    9. In the Manage page, approve this newly registered resource. Now I did the following in Geoserver:   1. Registered a ArcSDE shapefile in Geoserver that I can see in "Layer Preview" and I can get the common formats like OpenStreet map working for the layer. But what is an OGC GetCapabilities URL??? Where do i declared and link it to the layer I registered? Many thanks for any help Carlos.


Moving from MapInfo to web mapping

I use MapInfo and have a lot of data in .tab format that I would like to put into a database in order to access it through a website. I would like to continue using MapInfo as my GIS programme and export to the web database with the minimum of fuss. What would be a good type of database for this? I was considering MySQL as I have used it a lot before and it has native spatial data types. However I have read that MapInfo can't connect directly to MySQL and that its spatial data support is a bit patchy. Alot of discussion on the subject is several years old though, so can anyone confirm or deny this? If not MySQL, what is a good choice? PostgreSQL/PostGIS?


Is ANZLIC Metadata supported using the new approach of ArcGIS 10.x?

I like the new approach to metadata at ArcGIS 10 because I fall into the camp that prefers minimal metadata for many purposes, but like to have the option to work with more comprehensive standards when they are appropriate. However, I'm keen to hear whether people think ArcGIS 10 works well for those trying to adhere to ANZLIC Metadata Profile 1.1? My understanding is that the ANZLIC Metadata Profile 1.1 (based on ISO 19115) is very similar to the ISO 19139 Metadata Implementation Specification so can the ArcGIS 10 metadata style called ISO 19139 Metadata Implementation Specification simply be used? Or is the ArcGIS/ANZLIC community awaiting the release of an ANZLIC specific metadata style for ArcGIS 10 by Esri or a third party? A reply to an ArcGIS Discussion Forum thread suggests that it might be.


Maps, black and white/ gray scale printing

I'm working on Mapserver to do beautiful maps with AGG rendering Until now, I've only thought about the best way to get color maps. Now, I have to think how to make color ones compatible with a grey scale printing http://en.wikipedia.org/wiki/Grayscale How do you deal with this? Any advices, tips or resources will be welcome. * * * Edit : My question is really more relative to visual perception than Mapserver. It seems I will have to read Semiology of Graphics by Jacques Bertin (reedited this year http://www.esri.com/news/releases/10_4qtr/bertin.html ) or look on Edward Tufte various book/articles http://www.edwardtufte.com/tufte/


Cutting up a shape file into smaller parts

I've downloaded the 350MB coastline shape file from Openstreetmap, but I really would like to work with a small part of them, is there a way to generate a shapefile of a small area. I'm using MacOSX or Unix for this. Edit: Apparently the coastline file from Openstreetmap wasn't one big polygon so everything worked fine. Without cutting polygons.


Updating a layer's bounding box in GeoServer after WFS-T insert

I'm creating a data store in GeoServer from an OpenLayers DrawFeature control, using WFS-T to save the geometry into GeoServer. This works great, except the data store's bounding box doesn't get set in GeoServer. However, if you create a data store from a shapefile, the bounding box gets automatically set. I know you can go into the GeoServer admin page and tell it to calculate the bounding box, but I need a way to do it programmatically.


Enable FIND in ArcGIS Server 10.1

Trying to get a FindTask working via this sample and the only difference between my code and the samples that I can see is that the sample data has a find operation enabled. I can't find in my server manager where to enable this. And is it even necessary? Right now the search isn't returning any results. ![ESRI Service](http://i.stack.imgur.com/9amAl.png)               <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN" "http://www.w3.org/TR/html4/strict.dtd">      Find Task               <style type="text/css">       html, body { height: 100%; width: 100%; margin: 0; padding: 0; overflow:hidden; }       body { font-size: 0.9em; font-family: Geneva, Arial, Helvetica, sans-serif;       } .details { font-weight: bold; } #grid { border:       1px solid #333;}     </style>     <script type="text/javascript">       dojoConfig = {         parseOnLoad: true       };     </script>     <script type="text/javascript" src="http://serverapi.arcgisonline.com/jsapi/arcgis/?v=3.2"></script>     <script type="text/javascript" language="Javascript">       dojo.require("dojox.grid.DataGrid");       dojo.require("dojo.data.ItemFileReadStore");       dojo.require("esri.map");       dojo.require("esri.tasks.find");       dojo.require("dijit.layout.BorderContainer");       dojo.require("dijit.layout.ContentPane");            var findTask, findParams, map;            function init() {         var initExtent = new esri.geometry.Extent({"xmin":-9199621.530456403,                 "ymin":3365235.5214724857,                 "xmax":-9077475.159256855,                 "ymax":3443965.6606061114,                 "spatialReference":{"wkid":102100}});         map = new esri.Map("map",{extent:initExtent});         var censusMapLayer = new esri.layers.ArcGISDynamicMapServiceLayer("http://gisprod2/arcgis/rest/services/Dynamic/Zoning/MapServer", {         });         map.addLayer(censusMapLayer);              //create find task with url to map service         findTask = new esri.tasks.FindTask("http://gisprod2/arcgis/rest/services/Dynamic/Zoning/MapServer");              //create find parameters and define known values         findParams = new esri.tasks.FindParameters();         findParams.returnGeometry = true;         findParams.layerIds = [4];         findParams.searchFields = ["ZONECLASS"];              dojo.connect(map, 'onLoad', function(theMap) {           //resize the map when the browser resizes           dojo.connect(dijit.byId('map'), 'resize', map,map.resize);         });       }                 function execute(searchText) {         //set the search text to find parameters         findParams.searchText = searchText;         findTask.execute(findParams, showResults);       }            function showResults(results) {         //symbology for graphics         var markerSymbol = new esri.symbol.SimpleMarkerSymbol(esri.symbol.SimpleMarkerSymbol.STYLE_SQUARE, 10, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID, new dojo.Color([255, 0, 0]), 1), new dojo.Color([0, 255, 0, 0.25]));         var lineSymbol = new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_DASH, new dojo.Color([255, 0, 0]), 1);         var polygonSymbol = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_NONE, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_DASHDOT, new dojo.Color([255, 0, 0]), 2), new dojo.Color([255, 255, 0, 0.25]));              //find results return an array of findResult.         map.graphics.clear();         var dataForGrid = [];         //Build an array of attribute information and add each found graphic to the map         dojo.forEach(results, function(result) {           var graphic = result.feature;           dataForGrid.push([result.layerName, result.foundFieldName, result.value]);           switch (graphic.geometry.type) {           case "point":             graphic.setSymbol(markerSymbol);             break;           case "polyline":             graphic.setSymbol(lineSymbol);             break;           case "polygon":             graphic.setSymbol(polygonSymbol);             break;           }           map.graphics.add(graphic);         });         var data = {           items: dataForGrid         };         var store = new dojo.data.ItemFileReadStore({           data: data         });         grid.setStore(store);       }       dojo.addOnLoad(init);     </script>      Find State/City/River: Layer Name Field Name Value


How to publish my map from geoserver

> **Possible Duplicate:**   >  How to develop a map with layers in Geoserver? I'm trying to create a simple map for my city. I have a shape file my city, but I can't create on the web. I tried to take from geoserver to use geojson format but I didn't. Can you help me please.


How to Display Live Incoming Smartphone Images with Ushahidi on a Map

I have not yet used Ushahidi and I could not find the answer on their website to my question. Can I easily do the following?   1. setup an Ushahidi server/map,    2. go outside to make photos with my Android/iPhone smartphone   3. upload these photos to Ushahidi which should be saved in a database   4. display the photos on the map according to the coordinates. Also few lines of comments would be nice. When I now visit my site and open the map, I would like to see some dots (actual position of an entry/uploaded image). By double clicking on the dot, the original uploaded image+comments should be displayed(pop up). Is this possible with Ushahidi? edit: nobody working with Ushahidi?


British National Grids projection does not add Ordinance Survey grids to shapefile

I have created a shapefile in ArcCatalog by right-clicking the folder I want it in, then selecting new, then shapefile, then changing the type to "polygon", and using the British National Grids under the projected co- ordinate system. I then tried to add the shapefile in ArcMap by clicking the Add Data button and selecting the shapefile I created. My problem is that when I open the shapefile in ArcMap, the map is still blank and I cannot see the projected shape of the British National Grids I need. Please help.


How to add more events like Control.Click and Handler.Point for a single variable in OpenLayers

I have two points which are marked dynamically on the map and when I mark the start point it is returning lon,lat, but I want to assign this coordinates to the variable that marked/handled the start symbol and a point. Ex:               measureControls={                   start: new OpenLayers.Control.DrawFeature(start, OpenLayers.Handler.Point),                   stop: new OpenLayers.Control.DrawFeature(stop, OpenLayers.Handler.Point),      for the same start and stop, I would like to add event 'click' and assign the point when clicked. Thank you


Is it possible to use symbols instead of colours to visualise raster cell values?

I have a raster image (GeoTIFF) with _integer_ values. I know how to visualise raster values in QGIS using colormaps or similar, but I would like to assign symbols to each colour value instead - just like a colormap, but with symbols. **Example:** ![enter image description here](http://i.stack.imgur.com/IRWBO.png) **Question:**   Is this possible and in that case how?


Only return the NumberOfFeatures in a WFS query

I'm running a GetFeature WFS query in Geoserver (looking for features within a polygon) but all I actually want is the total number of features. Is there a way of just returning this info (in the hope that the query will run a little faster) rather than all the records? The query is returning over 10,000 features so I hope it would make a difference if possible. Many thanks, Steve


QGIS project on QGIS server -- cannot load wms layer

I've created a qgis project on QGIS 1.8, copied it onto the QGIS Server, and am trying to load the project as a WMS Layer back into QGIS. There is one vector layer in the project. Paths have been set to absolute and the CRS for both the layer and project is EPSG 27700. In the "Add Layer(s) from a Server" dialog box, the project is showing, but not the vector layer. Clicking on the project, at the bottom of the dialog box it says "No Common CRS for selected layers". Any help is very welcome -- I have searched the forums and tried the solutions suggested QGIS Server Capabilities Response doesn't contain any Layers and QGIS Server Capabilities Response doesn't contain any Layers (Fedora 15 64-bit) but still cannot get it to work. Many thanks!


Greyscale Raster (tif) in ArcMap

I've searched all over google and have found nothing. I want some tiff files I have in ArcMap to be in greyscale with the rest of the layers in colour, similar to the style override > greyscale option in MapInfo. I will be printing A1 posters and do not want the background imagery in colour. I've tried changing the symbology to stretched (as suggested on a forum), but it looks weird, like this: ![enter image description here](http://i.stack.imgur.com/PuzKv.png)


Cropping elements from one layer to an element from another

I'm trying to make a hex map of a country (Russia, in this case). To do this, I've opened a map of Russia in QGIS (from Natural Earth), set up my projection, and created a hex grid with MMQGIS. So far, so good. Now I want to crop the hexes with the shape of Russia. I've tried using the intersect option, but I get this error:               "Input CRS error: One or more input layers missing coordinate      reference information, results may not be as expected."      and an empty shapefile for my troubles. What am I doing wrong? UPDATE: Ok, so now I have them both in the same coordinate system (based on an Albers projection using meters), but I'm still not able to get an intersection to work: > Warnings: Input CRS error: One or more input layers missing coordinate > reference information, results may not be as expected. GEOS geoprocessing > error: One or more input features have invalid geometry. Some output > geometries may be missing or invalid. The output shapefile is empty. I _was_ able to generate a version that almost works using WGS84 coordinates, but that had some problems with the 180° meridian... ![](http://f.cl.ly/items/3E0p2V3Y3d1F2X0C3T24/QGIS.png)


automatic templates KML

I have imported several points from QGIS into Google Earth. When I click over an icon point, I get a ballon with all the attributes, OK. One of the attributes is called Text. Ok, so what if i want to take the Text value to be used in ot modify the balloon? I don't want to have to manually write, point by point, the value. i have seen in this URL that you can use extendedData to create a template. The KML generated by QGIS, instead of this example snippet (used in Google tutorials)               <Placemark>     <styleUrl>#golf-balloon-style</styleUrl>     <ExtendedData>       <Data name="holeNumber">         <value>1</value>       </Data>     </ExtendedData>      QGIS creates this...               <Placemark>     <styleUrl>#golf-balloon-style</styleUrl>     <ExtendedData>     <SchemaData schemaUrl="#punts_R">       <SimpleData name="Text">my text</SimpleData>         </SchemaData>     </ExtendedData>      I wonder if this QGIS created schema there is some way the parameters replacement may work (I didn't get it work, but yes with the Google 'standard'). It would save me a lot of work...


Why isn't StackBuilder installation of PostGIS working?

I've installed PostGIS (ver. 2.0.1) on my Mac using EnterpriseDB's StackBuilder, which I also used for my PostgreSQL (9.2.2) installation. Although StackBuilder says that PostGIS is now installed, and although I have a PostGIS directory in /usr/local/PostgreSQL/9.2/, it doesn't appear to be actually installed (see queries below).               gregory=# select postgis_full_version();     ERROR:  function postgis_full_version() does not exist     LINE 1: select postgis_full_version();                    ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.          gregory=# create extension postgis;     ERROR:  could not open extension control file "/usr/local/PostgreSQL/9.2/share/postgresql/extension/postgis.control": No such file or directory      and `select * from pg_available_extensions;` does not return a record for PostGIS. Does anyone know what's up with this and how to get it to work?


Unable to display FeatureLayer using ArcGIS server

I can't figure out why this FeatureLayer doesn't show up on top of my base map. I'm using ArcGIS Server 10.0 and the ArcGIS JS API 2.6. Here's the information from the service's REST parameters page:               Feature Layer: Inspection Locations (ID: 5)     Display Field: NAME     Type: Feature Layer     Geometry Type: esriGeometryPoint          Min. Scale: 250000     Max. Scale: 0          Extent:      XMin: -8879044.1973     YMin: 4938432.3627     XMax: -7999181.0204     YMax: 5622923.0498     Spatial Reference: 102100 (3857)     .     .     .     Supported Interfaces:   REST          Supported Operations:   Query   Add Features   Update Features   Delete Features   Apply      Edits ...and here's the code:                 <script type="text/javascript"> var djConfig = { parseOnLoad: true }; </script>       <script type="text/javascript" src="http://serverapi.arcgisonline.com/jsapi/arcgis/?v=2.6"></script>       <script type="text/javascript">djConfig = { parseOnLoad:true };</script>            <script type="text/javascript">            dojo.require("esri.map");       dojo.require("dijit.layout.ContentPane");       dojo.require("esri.layers.FeatureLayer");       var map;            function init() {              var initExtent = new esri.geometry.Extent({           "xmin": -8879044.1973,           "ymin": 4938432.3627,           "xmax": -7999181.0204,           "ymax": 5622923.0498,           "spatialReference": {             "wkid": 102100           }              });         map = new esri.Map("map", {extent: initExtent });              var basemap = new esri.layers.ArcGISTiledMapServiceLayer("http://server.arcgisonline.com/ArcGIS/rest/services/World_Topo_Map/MapServer");              var featureLayer = new esri.layers.FeatureLayer("http://gis4/ArcGIS/rest/services/NY/NYeFind/FeatureServer/5",{          mode: esri.layers.FeatureLayer.MODE_SNAPSHOT,          outFields: ["*"]        });              map.addLayers([basemap, featureLayer]);              var resizeTimer;         dojo.connect(map, 'onLoad', function(theMap) {           dojo.connect(dijit.byId('map'), 'resize', function() {             clearTimeout(resizeTimer);             resizeTimer = setTimeout(function() {               map.resize();               map.reposition();             }, 500);           });         });              } // end init              dojo.addOnLoad(init);      Oddly enough, if I hit the "View in ArcGIS.com Map" link, it works perfectly.


How to increase RAM allocation to QGIS?

I am using QGIS 1.8.0 and it runs slowly (with huge layers the size of NSW, Australia) and it appears that QGIS is limited in its RAM allocation to 128.9MB. How can I increase the speed with which it deals with such big files by allocating a larger amount of my computers power to QGIS? I am running Windows 8 Pro (64-bit) on: \- Intel Pentium P6200 @ 2.13GHz (dual core); \- 2.00GB RAM \- AMD Radeon HD 6300M Series   * Engine clock speed: 500-750 MHz   * Processing power (single precision): 80-120 GigaFLOPS   * Polygon throughput: 125-187.5M polygons/sec   * Data fetch rate (32-bit): 16-24 billion fetches/sec Texel fill rate (bilinear filtered): 4-6 Gigatexels/sec Pixel fill rate: 2-3   * Gigapixels/sec Anti-aliased pixel fill rate: 8-12 Gigasamples/sec   * Memory clock speed: 800-900 MHz Memory data rate: 1.6-1.8 Gbps Memory bandwidth: 12.8-14.4 GB/sec


count lines in polygons

I found the following link to a post saying that a count lines in polygons had been added to ftools for QGIS: http://hub.qgis.org/issues/4712 However I can't find how to access this feature nor any online discussion of how to do this - can anyone help? Alternatively, are there any other tools that will allow this? I have a hexagonal grid and want to count the number of animal track paths that go through each cell for a migratory distribution map. Thanks


Switching career from GIS Technician/Analyst to GIS Web Map Developer/Programmer?

I am a GIS Analyst at work that dabbles in python scripting in ArcGIS and the Google Maps JavaScript API V3. My professional goal is to either be in web mapping development or desktop/server GIS application development. I have seem to hit a plateau and I am not sure how to progress with my learnings or professionally to reach my goal. I ask those of you out there that made the jump from technician/analyst how did you do it? How can I land a job in web mapping/application development with little to no experience.


Binding IActiveViewEvents::ItemReordered to method in my ATL class

May you help me please to know how I can wire this event to my method? I have an ATL class and I need to know when a layer added to map object so that I can update my list of layers in this class.My class implements IExtension and I've put this code in Startup method of this interface: //Some code to obtain a reference to Map object __hook(&IActiveViewEvents::ItemReordered,pMap/ _pointer to map_ /,&MyClass::MyEventHandler); and in Shutdown method of IExtension I unhook it,but I receive compile error says the source and handler method not compatible while both signatures,arguments and names are exactly same. Is it true a all this method?any simpler method for this wiring?or any reference on net? Regards


How to interpolate daily climate data using ArcGIS for Desktop?

I want to interpolate daily climate data available from 3 locations. I am not sure how to interpolate the data for daily values for 30 years from the available 3 locations.


How to display coordinates of a point in a projected cooridante system?

I want to check the coordinates of a point. On selecting this point in ArcGIS it displays a latitude and Longitude . In ArcGIS I found an option in 'Identity Tool' that the Latitude and Longitude can be changed into desired units say for eg: Decimal degrees, decimeters etc. But when the same file is opened in QGIS and the same point is selected it shows me a Latitude and Longitude point but its unit can not be adjusted. In QGIS, is there a method to change the display units of the point, similar to how ArcGIS does this?


No difference between distances measured in projected and non-projected data

I'm performing measurements of distances in Km between a lot of localities by using the "Point Distance" in ArcGis 9.3. I realized that the results obtained after using projected coordinate system are the same of those obtained without projecting data. I change projected coordinate system in "Data frame properties". Where am I wrong? I am a newbe in Gis and in ArcGis too. Any advices will be invaluable. Thanks in advance.


Route optimisation for multiple vehicles

I have 100 destinations and 5 vehicles and I need to code a solution that efficiently routes the vehicles to each destination, so each destination is visited by 1 of the vehicles. Some of these destinations could have time windows to be visited during as well. I'm using PostgreSQL and PostGIS, and was looking to use pgrouting but I'm not sure it's suitable for multiple vehicles - from my limited knowledge of Dijkstra etc I think they're designed for a single vehicle. Any thoughts on whether pgrouting can solve this problem, and if so, any code examples? If not, are there any open-source alternatives that can do it?


How to sensibly digitise this raster, scriptably, in ArcGIS

I have the raster below which I want to convert to a set of vector lines: ![Raster Lines](http://i.stack.imgur.com/m1b8k.png) As you can see, there are small breaks in the lines and some bits are double- width. I want to be able to reliably import these sort of images to a vector in ArcGIS. I've spent a lot of time trying to get tools like Integrate, Extend Lines, Trim Lines and Unsplit Lines to do what I want, but I just can't seem to achieve it. I guess this problem must have been experienced before, does anyone have any ideas for a workflow that would do what I want? I thought of maybe trying to clean up the raster before-hand, but I'm not sure how to go about that. Again, any ideas would be much appreciated.


Maximum number of Tiled Layers in Open Layers

I'm developing an OpenLayers site where the client wants a dozen layers that can be individually turned on and off, which is a bit of a bummer in terms of performance. The data is PostGIS and static enough to cache tiles. The area is quite small geographically but requires 9 zoom levels to allow detailed views. I'm using a cached tiled basemap which is fast and efficient. I've been using Single Tile WMS layers for the other dozen layers, but am wondering if I should used Tiled, cached versions of these. There are other efficiency gains to be made by limiting the buffer, and adding sub-domains such that images can come from tiles.example.com, tiles2.example.com, etc, however all other efficiency strategies aside, I'm wondering at what point the OpenLayers performance difference is for 12 untiled layers, verses 12 tiled layers? Said another way, at what number of tiled layers in OpenLayers will the performance degrade such that Single Tile layers are preferable? If I have 12 tiles per view, and 12 layers, can OpenLayers handle 144 tiles without bogging down?


Creating fixed width polygons in ArcMap

I've been adding easements, which typically have a fixed width. Are there any existing tools or scripts that would allow me to define one edge and automatically create the other? The buffer tool is poor, the line offset tool is for lines. I want to be able to use the same tools (sketch, intersect, endpoint arc) and have actual curves,not point approximations (what buffer does).


Where does "Update measurements" put its updated measurements? (MeemsTools)

I have several .shp files loaded into MapWindows GIS. I go to plug-ins > Meems Tools and click on "UpdateMeasurement Plugin", and it loads. I click on the Update measurements, and in the the Update measurements window that appears, I click the "Update measurements" button, without changing any of the options in the window. The "select attribute" section grays out, and the bottom of the window lists "Finished updating measurement values", next to the "MapWindow Wiki help" link. > Where does it put these updated measurement values?    possibly relevant information: I'm using Windows 7 version 6.1 .   The files are loaded as multiple layers in the same window of MapWindow GIS.   The files are simplifications using MapShaper.org of the "district" files from "Boundary" on http://hawaii.wr.usgs.gov/maui/data.html .


Multi-Colored Polygon in Thematic Map

I have to make a thematic map with an attribute table with duplicated features but different attributes. See below:               FEATURE   ATTRIBUTE     -------------------     1         bl     1         gr     2         gr     2         re     3         re     4         bl     5         ye     5         ye      Does anyone knows how to paint a feature (polygon) with the corresponding colors? I mean like for feature 1 paint half blue and half green, for feature 2 half green half red... Thanks.


Need gps trackers for social science research

I'm interested in studying people's movement over time (24 hours), and would like to do so with devices that also measure elevation very accurately. I understand there's stuff out the where I can buy these for about $20 a piece, as I need about 30 of them. The data will be used in arcgis. Can someone recommend a specific brand to buy?


Bing Map Aerial Question

I am aware of Microsoft Bing Maps Aerial licensing changes with ArcGIS for Desktop at the end of this year. My question here is why I have two different ArcMap I saved had them in it. One was able to show it and the other would not allowed me to show me? What is the reason for that? ArcGIS 10.1 SP 1 Here is screenshot: ![Here is the Sample of the screenshot](http://i.stack.imgur.com/oZGlQ.jpg)


How to perform Shapefile to Raster with ArcMap 10.1 while preserving resolution and values?

I am new to ArcMap and I would really appreciate some help. I want to create a new raster using a value in a column of a point shapefile. I have a point shapefile with points equally spaced and each placed at the center of each pixel of a raster layer. (I created this shapefile using 'Raster to point' tool) Now I want to create a raster with the same resolution of the original raster but using as value the values of column in the shapefile table. How can I do that? I tried running 'Point to raster'(black squares)but it's not doing what I am looking for. ![enter image description here](http://i.stack.imgur.com/qaPra.png) This is what happens running 'point to raster' setting the cell size as the original raster: ![cell size set equal to original raster](http://i.stack.imgur.com/Agr2q.png)


ArcSDE installation; "sdesetup: command not found"

I am trying to install ArcSDE component for postgresql database from A DVD file. I have copied files of following folders on the server :   1. `command_64`   2. `pg_64` in both of theses files there are `install` files. for both I used :               ./install -load      yes     mount point: DVD location or wherever you are installing from     install directory: enter to accept /home/sde     Package numbers to load: enter to accept all     yes     no      now I have sdeexe101 file which contains bin, etc, liband locale folders in it. I also As the root user, copied the st_geometry.so, libsg.so, and libpe.so files from $SDEHOME to the lib directory in the postgres location.               [root@machine pgsql]# cp /home/sde/sdeexe101/lib/st_geometry.so /usr/lib/pgsql     [root@machine pgsql]# cp /home/sde/sdeexe101/lib/libsg.so /usr/lib/pgsql     [root@machine pgsql]# cp /home/sde/sdeexe101/lib/libpe.so /usr/lib/pgsql      when I tried to execute sdesetup to create geodatabase and authorize it using:               sdesetup -o install -d POSTGRESQL -u sde -p sde...      I got this: **`sdesetup: command not found`** what am I missing here? Are pg_64 and commands_64 files that I needed? I appreciate any help!


GeoServer and GeoExplorer issues setting up per-layer security

I can't get per-layer security working... I have the following set up:   1. A user USER, with no group membership and one role EXAMPLE_ROLE and a password   2. the role has no parent role and no properties   3. under data security I have:      * _._.r set to *      * workspace.mylayer.w set to EXAMPLE_ROLE   4. All the services security is set to _._ * But, when USER logs in to geoexplorer's saved map with mylayer in there, USER can edit, but saving returns the error message "{server}mylayer is read-only. If I set the second secirity to workspace.mylayer.w * then that user can edit. What's going on? What am I missing? EDIT: using version 3.1 of OpenGeoSuite, GeoServer 2.2-SNAPSHOT EDIT2: on looking into it further, I noticed that I don't get the login button on geoexplorer. I also get this error on Chrome's console:               Failed to load resource: the server responded with a status of 404 (Not Found) http://tobler-v:8080/geoserver/rest/imports      and I wonder if this is something related to this: GeoExplorer provokes "401 Unauthorized" request and lacks login


How do I get an accurate time attribute when importing a GPX layer as a vector

When I import a GPX file in QGIS (as a vector as track points) and I check the time attribute it appears to be 4 hours later than it should be. Interestingly when I take this same GPX file and use GPSBabel to convert it to a CSV file I find that the time listed in the CSV is the correct time. So am I doing something wrong or is this a bug in QGIS. Thanks for the help. David


How to generate an image from OSM on command line?

Is there a GNU/Linux command line tool that allows to generate an image file from OSM, preferably as an SVG?


Which framework is best to implement an application robustly, using ArcGIS Server JS API?

Over on the Esri forum, user Matt Moyles suggested that the approach used in the Esri JS samples isn't suitable for robust development of a web mapping application using HTML, JavaScript and CSS: > ESRI's suggested approach to javascript application development is ancient > and outdated. I would not recommend following the samples. Dojo 1.7 supports > AMD with Asynchronous dependency loading. I would start with the dojo > boilerplate template and "work" the arcgis api into that. The samples are > not suitable for serious application developers. They are mostly just proof > of concept snippets. > > A serious application should be developed with the aid of some sort of > framework to help structure things. I have been using dojox.mvc with great > success! But other options include things like backbone.js, spine.js, or > even javascript MVC. > > Dojo Boilerplate - https://github.com/csnover/dojo-boilerplate   * does anyone agree/disagree with this statement?   * are there any online examples of better approaches to web application development using the ArcGIS Server JS API?


Extract data from QueryLayer to Shapefile

I have a QueryLayer that I need to export as a shape file. This was previously implemented as an ArcSDE feature class which we converted to a QueryLayer (running against SQLServer) in order to have "live" access to data coming in from a variety of systems. We have an export tool built in model builder that worked with the ArcSDE feature class using the ExtractData toolbox, however this doesn't seem to work against the QueryLayer. I would appreciate any insight or ideas! Thanks! Below is the code for the data extract:               import arcgisscripting, os, sys, traceback, zipfile     gp = arcgisscripting.create(9.3)          def getPRJFile(inputCoordSysString, prjDir):         inputCoordSysString += ".prj"         # walk through the dirs to find the prj file         if os.path.exists(prjDir):             for x in os.walk(prjDir):                 if inputCoordSysString in x[2]:                     return True, os.path.join(x[0], inputCoordSysString)         else:             return False, ""              # if we got to here then it didn't find the prj file         return False, ""          def zipUpFolder(folder, outZipFile):         # zip the data         try:             zip = zipfile.ZipFile(outZipFile, 'w', zipfile.ZIP_DEFLATED)             zipws(str(folder), zip, "CONTENTS_ONLY")             zip.close()         except RuntimeError:             # Delete zip file if exists             if os.path.exists(outZipFile):                 os.unlink(outZipFile)             zip = zipfile.ZipFile(outZipFile, 'w', zipfile.ZIP_STORED)             zipws(str(folder), zip, "CONTENTS_ONLY")             zip.close()             gp.AddWarning("  Unable to compress zip file contents.")          def zipws(path, zip, keep):         path = os.path.normpath(path)         # os.walk visits every subdirectory, returning a 3-tuple         #  of directory name, subdirectories in it, and filenames         #  in it.         for (dirpath, dirnames, filenames) in os.walk(path):             # Iterate over every filename             for file in filenames:                 # Ignore .lock files                 if not file.endswith('.lock'):                     #gp.AddMessage("Adding %s..." % os.path.join(path, dirpath, file))                     try:                         if keep:                             zip.write(os.path.join(dirpath, file),                             os.path.join(os.path.basename(path), os.path.join(dirpath, file)[len(path)+len(os.sep):]))                         else:                             zip.write(os.path.join(dirpath, file),                             os.path.join(dirpath[len(path):], file))                           except Exception, e:                         gp.AddWarning("    Error adding %s: %s" % (file, e))              return None          def createFolderInScratch(folderName):         # create the folders necessary for the job         folderPath = gp.CreateUniqueName(folderName, gp.scratchworkspace)         gp.CreateFolder_management(gp.scratchworkspace, os.path.basename(folderPath))         return folderPath          def getTempLocationPath(folderPath, format):         # make sure there is a location to write to for gdb and mdb         if format == "mdb":             MDBPath = os.path.join(folderPath, "data.mdb")             if not gp.exists(MDBPath):                 gp.CreatePersonalGDB_management(folderPath, "data")             return MDBPath         elif format == "gdb":             GDBPath = os.path.join(folderPath, "data.gdb")             if not gp.exists(GDBPath):                 gp.CreateFileGDB_management(folderPath, "data")             return GDBPath         else:             return folderPath          def makeOutputPath(raster, inLayerName, convert, formatList, zipFolderPath, scratchFolderPath):         outFormat = formatList[1].lower()              # if we are going to convert to an esri format on the clip, put the output in the zipfolder         # else put it in the scratch folder in a gdb         if convert:             outwkspc = getTempLocationPath(zipFolderPath, outFormat)         else:             outwkspc = getTempLocationPath(scratchFolderPath, "gdb")              if inLayerName.find("\\"):             inLayerName = inLayerName.split("\\")[-1]              # make sure there are no spaces in the out raster name and make sure its less than 13 chars         if outFormat == "grid":             if len(inLayerName) > 12:                 inLayerName = inLayerName[:12]             if inLayerName.find(" ") > -1:                 inLayerName = inLayerName.replace(" ", "_")              # make the output path         tmpName = os.path.basename(gp.createuniquename(inLayerName, outwkspc))         tmpName = gp.validatetablename(tmpName, outwkspc)              # do some extension housekeeping.         # Raster formats and shp always need to put the extension at the end         if raster or outFormat == "shp":             if outFormat != "gdb" and outFormat != "mdb" and outFormat != "grid":                 tmpName = tmpName + formatList[2].lower()              outputpath = os.path.join(outwkspc, tmpName)              return tmpName, outputpath          def clipFeatures(lyr, where, featureFormat, zipFolderPath, scratchFolderPath, convertFeaturesDuringClip):         global haveDataInterop         try:     ##        if not convertFeaturesDuringClip and not haveDataInterop:     ##            raise "LicenseError"                  # get the path and a validated name for the output             layerName, outputpath = makeOutputPath(False, lyr, convertFeaturesDuringClip, featureFormat, zipFolderPath, scratchFolderPath)                  # do the clip             gp.select_analysis(lyr, outputpath, where)             gp.AddMessage("  clipped %s..." % lyr)                  # if format needs data interop, convert with data interop             if not convertFeaturesDuringClip:                 # get path to zip                 outputinzip = os.path.join(zipFolderPath, layerName + featureFormat[2])                 if featureFormat[2].lower() in [".dxf", ".dwg", ".dgn"]:                     gp.AddWarning("..using export to cad..")                     gp.ExportCAD_conversion(outputpath, featureFormat[1], outputinzip)                 else:                     if not haveDataInterop:                         raise "LicenseError"                     diFormatString = "%s,%s" % (featureFormat[1], outputinzip)                     # run quick export                     gp.quickexport_interop(outputpath, diFormatString)              except "LicenseError":             gp.AddWarning("  failed to export to %s.  The requested formats require the Data Interoperability extension.  This extension is currently unavailable." % featureFormat[1])             pass              except:             errorstring = gp.GetMessages(2)             if errorstring.lower().find("failed to execute (quickexport)") > -1:                 gp.AddWarning("  failed to export layer %s with Quick Export.  Please verify that the format you have specified is valid." % lyr)             elif errorstring.lower().find("failed to execute (clip)") > -1:                 gp.AddWarning("  failed to clip layer %s..." % lyr)             else:                 gp.AddWarning("layer %s failed..." % lyr)                 gp.AddWarning(gp.GetMessages(2))             pass          def clipAndConvert(lyrs, where, featureFormat, coordinateSystem):         try:             # for certain output formats we don't need to use Data Interop to do the conversion             convertFeaturesDuringClip = False             if featureFormat[1].lower() in ["gdb", "mdb", "shp"]:                 convertFeaturesDuringClip = True                  # get a scratch folder for temp data and a zip folder to hold             # the final data we want to zip and send             zipFolderPath = createFolderInScratch("zipfolder")             scratchFolderPath = createFolderInScratch("scratchfolder")                  # loop through the list of layers recieved             for lyr in lyrs:                 # temporary stop gap measure to counteract bug                 if lyr.find(" ") > -1:                     lyr = lyr.replace("'", "")                 describe = gp.describe(lyr)                 dataType = describe.DataType.lower()                      # make sure we are dealing with features or raster and not some other layer type (group, tin, etc)                 if dataType in ["featurelayer", "rasterlayer"]:                     # if the coordinate system is the same as the input                     # set the environment to the coord sys of the layer being clipped                     # may not be necessary, but is a failsafe.                     if coordinateSystem.lower() == "same as input":                         sr = describe.spatialreference                         if sr != None:                             gp.outputcoordinatesystem = sr                          clipFeatures(lyr, where, featureFormat, zipFolderPath, scratchFolderPath, convertFeaturesDuringClip)                 else:                     gp.AddWarning("  Cannot clip layer: %s.  This tool does not clip layers of type: %s..." % (lyr, dataType))                  return zipFolderPath              except:             errstring = "Failure in clipAndConvert..\n"             tb = sys.exc_info()[2]             tbinfo = traceback.format_tb(tb)[0]             pymsg = "ERRORS:\nTraceback Info:\n" + tbinfo + "\nError Info:\n    " + \                     str(sys.exc_type)+ ": " + str(sys.exc_value) + "\n"             errstring += pymsg             raise "clipAndConvertError", errstring          if __name__ == '__main__':         try:             # Get the Parameters             layers = gp.getparameterastext(0).split(";")             where = gp.getparameterastext(1)             inputFeatureFormat = gp.getparameterastext(2)             coordinateSystem = gp.getparameterastext(3)             outputZipFile = gp.getparameterastext(4)                  if gp.CheckExtension("DataInteroperability") == "Available":                 gp.CheckOutExtension("DataInteroperability")                 haveDataInterop = True             else:                 haveDataInterop = False             # Do a little internal validation.             # Expecting "long name - short name - extension             # If no format is specified, send features to GDB.             if inputFeatureFormat == "":                 featureFormat = ["File Geodatabase", "GDB", ".gdb"]             else:                 #featureFormat = inputFeatureFormat.split(" - ")                 featureFormat = map(lambda x: x.strip(), inputFeatureFormat.split("-"))                 if len(featureFormat) < 3:                     featureFormat.append("")                                          coordinateSystem = 'same as input'                  # Do this so the tool works even when the scratch isn't set or if it is set to gdb/mdb/sde             if gp.scratchworkspace is None or os.path.exists(str(gp.scratchworkspace)) is False:                 gp.scratchworkspace = gp.getsystemenvironment("TEMP")             else:                 swd = gp.describe(gp.scratchworkspace)                 wsid = swd.workspacefactoryprogid                 if wsid == 'esriDataSourcesGDB.FileGDBWorkspaceFactory.1' or\                    wsid == 'esriDataSourcesGDB.AccessWorkspaceFactory.1' or\                    wsid == 'esriDataSourcesGDB.SdeWorkspaceFactory.1':                     gp.scratchworkspace = gp.getsystemenvironment("TEMP")                  # clip and convert the layers and get the path to the folder we want to zip             zipFolder = clipAndConvert(layers, where, featureFormat, coordinateSystem)                  # zip the folder             zipUpFolder(zipFolder, outputZipFile)              except:             tb = sys.exc_info()[2]             tbinfo = traceback.format_tb(tb)[0]             pymsg = "ERRORS:\nTraceback Info:\n" + tbinfo + "\nError Info:\n    " + \                     str(sys.exc_type)+ ": " + str(sys.exc_value) + "\n"             gp.AddError(pymsg)


Determine what value used in Editor.OnCreateFeature

How can I determine what Value-Layer (not the layer itself, but the different colors for different attributes) was used in ArcGIS CreateFeature Event, when I insert a new feature using the Editor?


How do I generate a "granule" index shapefile for ImageMosaic in Geoserver?

I am experiencing problem with generating granule index shapefile for the ImageMosaic plugin in Geoserver; there are no clear instructions to follow. Is there anyway that I can generate it automatically using raster dataset? I have a huge raster dataset to serve as WMS. Thanks in advance.


Is there a way to export metadata to HTML with FME

I'm trying to export metadata from ArcSDE to HTML. I've been trying to use the XSLTProcessor transformer with the ArcGIS.xsl stylesheet. I can't seem to make it work. Any ideas. Is there another way to achieve this? Thanks


How do I overlay a KML on leaflet 0.4.4

I am using Leaflet 0.4.4, I have a KML I would like to overlay on the map, how do I go about doing this? **Edit:** I am already aware of https://github.com/CloudMade/Leaflet/issues/14, which directs you to download a plugin from http://psha.org.ru/b/leaflet-plugins.html. I was hoping to have a demonstration of the implementation of how to add a KML overlay. **Edit [Solved]:** I found the demo, here it is:               <script type='text/javascript'>         var map = new L.Map('map', {center: new L.LatLng(58.4, 43.0), zoom: 11});         var osm = new L.TileLayer('http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png');         var track = new L.KML("KML_Samples.kml", {async: true});              track.on("loaded", function(e) { map.fitBounds(e.target.getBounds()); });              map.addLayer(track);         map.addLayer(osm);              map.addControl(new L.Control.Layers({}, {'Track':track}));     </script>       You also need to make sure that you need to enable cross-origin resource sharing, I did this in PHP with:                 <?php           header("Access-Control-Allow-Origin: *");       ?>      The problem I was having was the Access-Control-Allow-Origin. Thanks all


ArcGIS 10.1 Crashes when updating address locator

I keep getting a "ArcGIS has encountered a serious error and needs to close" crash dialogue box whenever I try to update an address locator imported from one server to another (updating the location of the reference database). I've tried closing ArcMap and just using ArcCatalogue, still doesn't work. I'm on a 3Ghz/4G RAM relatively new (<1 year) PC at work. Can anyone help me with a work-around for this?


Has anyone had any success using the Stack Profile tool in ArcGIS with a terrain dataset as a profile target?

Has anyone had any success using the Stack Profile tool in ArcGIS for Desktop [Advanced] with a terrain dataset as a profile target? I have successfully used this tool with a raster dataset as a profile target but when I use a terrain dataset the tool fails. It writes out an empty table and returns > **Error 999999** with the message: > **table not found**


Creating Contours from a Classified Raster - ArcGIS

I have a group of classified raster files, showing point densities at 3 intervals. I want to create contour lines based on this data at these same intervals. I'm guessing this should be fairly easy, but how? I'm on ArcGIS 10.1. EDIT: I have about 100 raster files, each of them is a group of points which I've performed a Kernel Density Estimation on, so is now a raster. I've classified these into 3 classes 15-50, 50-100, and 100-250 (there is also a 0-15, but I don't need this included). How would I go about creating contours at each of the limits of the 3 classes I want to keep? Martin - Deleting extra contours would be fairly time-consuming. The operation would preferably be done as a batch - they all have the same class symbology and they would all have the same contour limits too.


Creating Layouts for selected windows in MapInfo

This is probably a very basic question, so I apologise in advance. I have a workspace from which I need to create a lot of different layouts, which need to be kept. Each time I create a layout which has five windows (2 of which are the same, copyright and key), I then minimize them all and open a new map window for the next site. When I then open a new layout window for all currently open windows, every window appears in the layout and I have to delete all but the ones that I need. Is there a way of specifying which ones I want to be included each time?


Plotting 1M+ intensity reading as custom shapes with Google Maps

I currently have a database (Mongo, atm) holding roughly 1M points of data corresponding lat/lon/intensity values that I'd like to plot over a region using the Google Maps API. Ideally, I'd like to have full control over the representation of the points, as the intent is to allow the end-user to control color/opacity/shape through a series of controls. I also need the ability to specify a color range which maps back to the intensity of a given point. That said, I've investigated using Google Fusion Tables, but for a heavily customized app I found that display styles limiting. I've also looked a bit into marker clustering, as well as a number of external plugins, but none seem to provide the option to generate the marker type dynamically. Any help/advice/direction is greatly appreciated. Thanks!


FME Server SOAP or REST

I would like to use SOAP or REST with FME Server 2011. But for SOAP, Safe propose to send an email for more information. Does anyone have experience with this Api ? Any suggestions much appreciated. Thanks _Tools: Visual Studio 2010 / CSharp / Wpf_


How to append geographic features into a .dgn file?

How can I export all the geographic features of a featureclass into .dgn (Microstation) file using C#/ArcObjects? I was able to create a new .dgn file and I want to export all the geographic features of a selected feature class into that newly created .dgn file.


How to fix pgRouting error "starting vertex not found"?

I've posted on here before about my battles with drivetime polygons and got some sensible suggestions back, so having Googled myself stupid over this I'm turning to you once again in the hope of a nugget of wisdom. I'm on a Windows 7 PC. I've followed the excellent guide at underdark.wordpress.com and got postgreSQL installed, with PostGIS and pgRouting. I've downloaded OSM data for the UK from Cloudmade (and Geofabrik for that matter) and I've used the wonderful osm2po to get that data into the postgreSQL database. So far, so good. No errors, no problems, looks like a database to me. I've got a table called network with obvious column names like id, source, target, kmh, cost etc. Just what I would have expected if I'd defined them myself. So I've fired up the SQL query panel, banged in the SELECT * FROM Driving_Distance statement and fed it the id, source, target and cost columns. I was trying to find 30 minutes around Northampton and the output did indeed some roads within Northampton(shire) - but the rest were in North Wales, which can't possibly be right. At 3:30am this morning it occured to me that the table also had columns named osm_source_id, osm_target_id and osm_id so perhaps it might be a bit more sensible if I used those. Made sense to me as most things do at that time of the morning, so I ran it again and now I just get **"starting vertex not found"** every single time. So to summarise. This works, but doesn't find the right roads:               SELECT network.id,network.x2,network.y2,network.geom_way      FROM network      JOIN      (SELECT * FROM driving_distance('       SELECT id,     source,     target,          cost/60::float8 AS "cost"       FROM network',       1516905,       0.025,       false,       false)) AS route      ON      network.id = route.vertex_id;      This doesn't work at all:               SELECT network.id,network.x2,network.y2,network.geom_way      FROM network      JOIN      (SELECT * FROM driving_distance('       SELECT osm_id::int4 AS "id",     osm_source_id::int4 AS "source",     osm_target_id::int4 AS "target",         cost/60::float8 AS "cost"       FROM network',       1516905,       0.025,       false,       false)) AS route      ON      network.id = route.vertex_id;      The same message happens whether or not I'm using Cloudmade or Geofabrik data, so I'm at a bit of a loss. Any suggestions at this point would be extremely welcome as my deadline is long passed and I don't think I can come up with any new reasons to the boss why it's not done! Thanks Michael


Split Grid data into 4

I am using ArcGIS software I have a Feature Class layer that represent a Grid network, which consists of rectangles that covers all the areas. I want to use it to generate another feature class layer with smaller grid, by dividing each rectangle into 4 quarters. What is the best way to do it with ESRI software? is there a geoprocessing tool that can do that?


Best tutorials on raster data formats and raster manipulation

I am looking to improve my understanding of raster formats and how to manipulate them. I don´t really want to go into the mathematics of it, but really want to know the differences between the different formats available, best practices, performance optimization etc. If anyone has any good sources which they would like to share, I´d be most grateful. In particular I would like to know the benefits of using 3 band TIFFs, JPEG2000, ECW against ESRI formats such as image mosaic, gdb raster, or raster catalog. I am having huge problems holding large amounts of aerial imagery in ArcGIS 10x. I cannot create mosaics or raster catalog at the moment and would also like to cut down the amount of file space used by these rasters. I get the images in TIF format @350mb per tile. I first wanted to hold them in ecw format, but then found that ecw are unclean around the edges making it impossible to change areas outside of the image into nodata. There are obviously many questions here. Sorry


OGC filter intersects polygon

I am trying to get a selection of points in my dataset, by selecting a polygon and returning all points inside it. I am using an OGC fiter, and currently I am testing it directly in the OpenLayers view of the data on geoserver. When I use the construct:               <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">         <ogc:Intersects>             <ogc:PropertyName>Geom</ogc:PropertyName>             <gml:Envelope xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:900913">                 <gml:lowerCorner>-161672.38316912 6970834.4247659</gml:lowerCorner>                 <gml:upperCorner>-161290.19802775 6971216.6099073</gml:upperCorner>             </gml:Envelope>         </ogc:Intersects>     </ogc:Filter>      only the points inside the box are returned, as expected. However, I want to use a polygon, so as a test I constructed the same box with the gml:               <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">         <ogc:Intersects>             <ogc:PropertyName>Geom</ogc:PropertyName>             <gml:Polygon xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:900913">                 <gml:outerBoundaryIs>                     <gml:LinearRing>                         <gml:coordinates decimal="." cs="," ts=" ">                             -161672.38316912,6970834.4247659 -161290.19802775,6970834.4247659 -161290.19802775,6971216.6099073 -161672.38316912,,6971216.6099073 -161672.38316912,6970834.4247659                         </gml:coordinates>                     </gml:LinearRing>                 </gml:outerBoundaryIs>             </gml:Polygon>         </ogc:Intersects>     </ogc:Filter>      When I use this filter, however, no points are returned (neither points inside the polygon nor outside). The eventual polygon I want to use is more complex than this but I thought I'd start with a simple one. Does anyone know why this is not returning any points?


Shapefile of Israeli sub-districts

is anyone aware of freely available shapefile of Israeli sub-districts? The file I found only contained data at the district level.


How do I query XY Event Layer in Model Builder?

I am trying to build a model to input table data, make xy event layer and then select features by an SQL expression and make that a shapefile. What tool should I use for the SQL expression to select the features?


estimate travelling time for the route using pgRouting and OSM data

Could someone share with me some ideas or your experience how to estimate travelling time based on the provided route within pgRouting? The problem is that roads inside and outside cities have the same type (for example "highway=primary") but different max speed allowed. But this max speed attribute is not always set in real OSM data. The only possible solution I see at the moment is to preprocess all roads before calculations and set for each some expected speed. For example, if "highway=primary" and it's inside the city then expected speed is 40km/h, if it's outside then expected speed = 90 km/h. But how to detect that the road is inside some city bounds or not? Any ideas are welcomed!


ArcSDE - use native geometry types for later migration away from SDE?

**Background:** We have a non-spatial reporting web application with a SQL Server backend. We would like to add a spatial component to this application: serve up some basemaps, serve up points and polygons for view and editing, build reports for selected polygons. To that end, the plan is to use ArcSDE with SQL Server, and use ArcServer to serve up base layers and vector geometries. (something like WMS/WFS-T) **Problem:** We may, at some point, decide we want to move away from expensive proprietary software, and use a FOSS alternative (e.g. GeoServer, or TinyOWS if/when it supports SQL Server). If possible, I'd like to avoid building the application irrevocably around ESRI technology. My understanding of SDE is that it can either use ESRI specific spatial formats, or spatial formats native to the database system it's installed on. (I presume there is some performance penalty for using the native formats.) **My Question(s):** If we were to set up SDE on this SQL server database using native SQL Server spatial formats, would we be able to set up say, GeoServer, to run against the tables that SDE was managing? What would it take to rip SDE out and use something else? (Also, a meta-question: is this a reasonable question to be asking? Is there something important I'm missing?)


Having difficult time to implement a Rule using OpenLayers.Filter.Function

I'm trying to set up a Rule using the Function Filter. But it does not work: the lines on that layer are just black. And more importantly, the filtering function is never executed. Can someone please point out what is the mistake that i'm making? Here is the code. Thank you for your time and kind concern.               var my_filter = new OpenLayers.Filter.Function(         function(attributes) {            console.log(attributes);            var x=0;             return true;              });     var ruleLow = new OpenLayers.Rule({         filter:my_filter ,       symbolizer: {pointRadius: 10, fillColor: "green",                    fillOpacity: 0.5, strokeColor: "green"}       });               var my_style=new OpenLayers.Style( null,ruleLow);          var my_style_map=new OpenLayers.StyleMap({         "temporary":my_style,         "default":my_style,         "select":my_style          });     this.vectors= new OpenLayers.Layer.Vector(         "Vector Layer",         {             styleMap:my_style_map,             renderers:this.renderer,         }     );


Merging Layers and assigning different values

> **Possible Duplicate:**   >  Merging Two or More Layers I am working with three clipped layers that are shown in the TOC. I created these layers from a large watershed shape-file. However, each clipped layer includes the layer before it and only expands on the one before it. I was wondering is there a way that I could combine the three layers without using the same features more than once and also at the same time assigning them different values. Pretty much as it is shown in the image that I have provided? How can I assign each clipped layer different values without repeating a feature? I have tried using the merge tool but then I don't know how to give each features different values. This may be a little confusing please let me know if you have any questions ![enter image description here](http://i.stack.imgur.com/0m35D.jpg) Thanks


How do I set a CRS when it isn't listed in QGIS? (Mac, not PC)

I want to work with a shapefile that is in a Coordinate System that I can't locate in QGIS for Mac. When I load this shapefile data, it shows up fine, on- the-fly has generated the CRS data most likely from the .prj file that has the info in it. However, if I add a coastline from WGS-84 coordinates, or any other data from another coordinate, my data ends up in Africa or even farther than that! So I try to set a coordinate system for it so QGIS will have a reference point. I type in the filter 102631 because that is the number for the Projection data, but that projection isn't listed. What do I need to do? Can I manually alter something so QGIS has a reference point for this projection data? The shapefile was created in a NAD-83 Alaska Coordinate System (NAD 1983 StatePlane Alaska 1 FIPS 5001 Feet ESRI 102631 http://spatialreference.org/ref/esri/102631/). On QGIS for the Mac this coordinate reference system is not listed. How can I correct this on the Mac? Is there a way to download more Projections?


PostGIS Raster: Outputting raster into GDAL supported file?

I have import a raster into my postGIS database via raster2pgsql. My raster tiled to 100x100. Now I'm trying to export (output) the raster to GDAL supported file (for example: PNG). Following the instruction in postGIS manual dev 2.1 (page 68). However, I don't get any output file in destination folder. Here's what I type in SQL tool in pgAdmin III               <!-- language: lang-sql -->     SELECT lo_export(demo.oid,'D:\demo_rast.png')     FROM       (SELECT oid,               lowrite(lo_open(oid, 131072), png) AS num_bytes        FROM (              VALUES (lo_create(0),                      ST_AsPNG(                                 (SELECT rast                                  FROM landsat                                  WHERE rid=1)))) AS v(oid,png)) AS demo     WHERE demo.oid = 1         WHERE demo.oid = 1


Given a known batch of mobile-obtained gps data with timestamps, what algorithms can be used to clean it?

I'm working on the backend for an iOS/Android app that collects locations every hour and syncs with the server every three hours. The mobile apps themselves are not under my control and so I can't tweak or even read the code (although I could ask the developers about specifics if need be). The server, given those "batches" of location data enhanced with timestamp and vertical/horizontal accuracy information, tries to determine visits to states/countries: if, for example, it gets some locations that can be placed in new jersey for a given day, it is said that the user visited new jersey once that day. It's based on **reverse geocoded** addresses to determine the places visited. The problem is that some of those locations can be noise: bad accuracy, "jumps" or readings from the very boundaries of states can affect the algorithm (which is very naïve: every location it finds, if it's on a different place, is considered as the beginning of another visit, unless it's too unique or too inaccurate). I was reading earlier about kalman filters, but, as I understand, they not only warrant a very good understanding of the mathematical model involved, they are more suited for real-time measurements. In my case, I already have enough inputs to not need so much prediction as I need smoothing. I'm starting to read about the **Douglas–Peucker algorithm** to maybe find a solution, but I don't know if I'm on the right track. What do you guys suggest?


ArcMap10: Symbology of Overlapping Polygons (Oil & Gas)

I am mapping an area where there are a lot of overlapping features because of certain Oil & Gas wells being drilled within the same geographic area, but to different depths and maintained by different companies. There might be one large Oil & Gas unit (polygon) that is drilled to 5,000' and another smaller unit (polygon) within the same geographic boundary that is drilled to 10,000'. Each polygon is symbolized based on company, so I have been messing with symbol levels to draw certain ones on top, but I don't feel like this is the best way to get ArcMap to automatically draw smaller ones on top and keep the larger ones on bottom. Is this possible? Thanks.![enter image description here](http://i.stack.imgur.com/fxdED.png)


Error after installing QGIS 1.8

I have a problem with the installation of QGIS 1.8.0. After finishing the installation two icons show on the desktop: Qgis desktop and Qgis browser. When I click on these, it shows the same message, that `there is no file msvcr70.dll` in the computer and that it can help to reinstall the program (it didn't). Thanks for any ideas.


Select KML feature with specific attribute

Im using Openlayers, and importing polygon data to the map with KML files. How could i select the feature which has the attribute name="TEST" ? KML def:               <Placemark>     <Name>TEST</Name>     <Polygon><outerBoundaryIs><LinearRing><coordinates>//Some coordinates\\</coordinates></LinearRing></outerBoundaryIs></Polygon>     </Placemark>      Iknow that you can loop trough them to check it. But i was hoping for something more efficient.


Is it possible to get the EPSG value from an OSR SpatialReference class using the OGR Python API?

When reading a layer from a OGR PostGIS connection I can get the SpatialReference of the layer, but is it possible to get the EPSG value? Is there any documentation on this? For example:               lyr = conn.GetLayerByName(tbl) # Where conn is OGR PG connection     srs = ly.GetSpatialRef()     print srs      Returns:               PROJCS["OSGB 1936 / British National Grid",     GEOGCS["OSGB 1936",         DATUM["OSGB_1936",             SPHEROID["Airy 1830",6377563.396,299.3249646,                 AUTHORITY["EPSG","7001"]],             AUTHORITY["EPSG","6277"]],         PRIMEM["Greenwich",0,             AUTHORITY["EPSG","8901"]],         UNIT["degree",0.01745329251994328,             AUTHORITY["EPSG","9122"]],         AUTHORITY["EPSG","4277"]],     UNIT["metre",1,         AUTHORITY["EPSG","9001"]],     PROJECTION["Transverse_Mercator"],     PARAMETER["latitude_of_origin",49],     PARAMETER["central_meridian",-2],     PARAMETER["scale_factor",0.9996012717],     PARAMETER["false_easting",400000],     PARAMETER["false_northing",-100000],     AUTHORITY["EPSG","27700"],     AXIS["Easting",EAST],     AXIS["Northing",NORTH]]      So how do I get the EPSG value for the projection? E.g.:               srs.GetEPSG()     print srs     27700      I've tried srs.GetAttrValue('AUTHORITY'), but this just returns 'EPSG'. Any help or pointers to documentation would be most appreciated. Thanks, Tom


How to re-project ease( Equal Area Scalable Earth) grid with a ~25 km cylindrical projection to WGS84 0.25 degree?

I have nc (NetCDF) files for global soil moisture. From the metadata, the projection is cylindrical and the resolution is 25 km (it is based on authalic sphere based on International 1924 ellipsoid). As I want to compare with other data, I have to make them identical. My other data are in WGS84 with 0.25*0.25 degree resolution and extent -180, 180, -90, 90. So I want to re-project the data I have from: > EASE grid, cylindrical, 25 km, 1383 pixel 586 lines to: > WGS84, 0.25*0.25 degree, 1440 pixel 720 lines


How to add layer in Arcscene using script tool

import arcpy     mxd = arcpy.mapping.MapDocument("CURRENT")     df = arcpy.mapping.ListDataFrames(mxd, "*")[0]     addLayer = arcpy.mapping.Layer(r"F:\project\Final\Ground floor\steps.lyr")     arcpy.mapping.AddLayer(df, addLayer, "BOTTOM")      this code is for adding a layer in arcmap. when i use the same code to add a 3D layer in arcscene i get the following error               cannot open map document          Failed to execute (addlayer).      please help. its important for my thesis work!


Arcpy: Counting number of segments where altitude rises or decreases based on DEM?

I have a DEM of Sudety mountains and linear .shp with several hiking paths. I'm trying solve one problem: How to make a script based on input DEM and input line shapes that would count number of segments of hiking path where the terrain rise is higher than 50 meters or lower than -50 meters and sum these that are >50m to one field and that are <50m to another field. I was trying to solve this problem with 3D line stack profiles and `arcpy.UpdateRowCursor`, but don't know how to order an iteration that would compare row.HEIGHT with row.HEIGHT of the previous row. Another question is how to name stack profile's output tables by a field name of particular row?


QGIS server installation and configuration for Ubuntu 14.0.4

Can anyone provide a step by step for installing and configuring QGIS server on Ubuntu? I have spent hours trying everything I could and installing everything I could but still can not get it to work.


ArcGIS Server -> Internet Service works well, Local Service fails

I have a web ADF application which connects to a map service (either via WMS or Internet Service using http). I'd like to use DCOM connections (faster as the documentation states), so I would like to use a 'Local Service' for the map service.. but only WMS or Internet Service types work, when I choose the Local Service type the map disappears and my app is no longer able to retrieve the map Has anyone ever experienced such an issue?


Is there a way to center an IPointCollection?

I built an IPointCollection and can display it, but I would like to display multiple copies in different locations. Is there a way to center a point collection on another IPoint?


KMZ Quality: Desktop vs Plugin

Is there a reason why the quality of KMZ images, when added to a website using the Google Earth Plugin, drop in quality (it appears blurry) compared to when the same KMZ is seen through Google Earth Desktop?


PGR_CreateTopology maintain Source and Target ID's

Is there a way to set the 'source' and 'target' manually and not have it overridden in `PGR_CreateTopology()`? I have a predefined node ID's (INODE, JNODE) in my table which is the same to a other tables in my database that I want to maintain for simplicity and consistency. This is happening:               GID  | INODE  | JNODE  | SOURCE  | TARGET     ---- | ------ | ------ | ------- | ------     1    | 25     | 6      | 1       | 2           and I want this:               GID | INODE | JNODE | SOURCE | TARGET     --- | ----- | ----- | -----  | ------     1   | 25    | 6     | 25     | 6


Force ArcCatalog to refresh TreeView and View after programmatically creating a table

I'm programmatically creating a table in a geodatabase and am trying to force ArcCatalog to refresh the display to show the new table. I've tried the following which does not work:               Dim gxApp As IGxApplication = TryCast(m_application, IGxApplication)     ... 'Create table here          'Then trying to force tree view to show new table...     gxApp.Refresh(0)          'And trying to force view to display...     gxApp.View.Refresh()      Missing something simple here.


Help using the US American Community Survey Data

I want to use the USA American Community Survey Data for spatial analysis, but I'm having difficulty understanding how the survey tables are constructed so that I can relate them to block-group shapefiles. I've successfully made a shapefile that contains all the US blockgroups, but I can't figure out how the census data is related to it using the STFID code. I have a few thousand text files (named `e20095ak0001000.txt` for example). When I read them into `R` they appeared with the following headers:               [1] "ACSSF"    "X2009e5"  "ak"       "X000"     "X0011"    "X0000001"     [7] "X4269"    "X2211"    "X224"     "X157"     "X189"     "X142"         [13] "X77"      "X269"     "X247"     "X169"     "X389"     "X146"  ... up to 200 or so      (I'm assuming the blockgroup code does not match the filename, as I have 12,465 files and 212,083 blockgroups) I realize I need to figure out what the the ID code is, to link these to the block groups, but I can't find documentation that explains this information (even after reading the technical documentation). Where are these column headers explained? I am using the 5 year ACS survey.


Whats the status for OSM in print composer QGIS 2.0?

I have a Stamen Toner OSM background on a report am doing. It seems impossible to use "Print" or "Save as image" functions since the OSM is being projected wrong with 10 km or so.. Right now I'm taking screenshots from there as a work around but it feels kind of wrong?      Screenshot with correct OSM projection:   ![Screenshot with correct OSM projection](http://i.stack.imgur.com/JzF1g.png)      Save from print composer with wrong OSM projection: (Buffers should be located at the lower right corner   ![Save from print composer with wrong OSM projection](http://i.stack.imgur.com/S04y2.jpg)


How to read arc/info e00 vector files (no binaries) with Python?

I can easily read e00 files with QGIS and it seems as if in older versions of ogr there was driver support for the arc/info e00 file format. Nowadays that seems to be deprecated, as running `ogrinfo --formats` doesn't list arc/info. I also found this useful link, a general description of the file format. Based on that information I would be able to write my own program to put all the pieces together. But hey, this must have been done before, right? Can somebody point me into the right direction?


Python crashes during large cursor process

In a current project, I am attempting to extract several fields for each record of a table in a GDB that has ~8.5 million records. My current strategy is to create a search cursor, fetch the data I want from each row, and pass it to a `csv.writerow` to write a line with the desired information. Unfortunately, after about 340k records, I get a message from Windows telling me Pythonw.exe has stopped working, and the process stops. Why would this happen? Is there a better method for extracting large amounts of data to text? DBF isn't an option because these files could become larger than the 2GB limit. Interestingly enough, after a few more tests, it appears that it **always** crashes on the 341051st record, and always at the same point on the line. **Edit** Now attempting to split up the cursor to read only a subset of the data each time. It is apparently a memory issue. I created a list for each subset query eg:               squery = ['OBJECTID >=0 AND OBJECTID <=300000',     'OBJECTID >=300001 AND OBJECTID <=600000']          And run the search cursor as an iteration:               for query in squery:         InRows = gp.SearchCursor(InFC,query)         for row in InRows:             val = row.val             val2 = row.val2             items = [val,val2]             OutCsv.writerow(items)         del row         del InRows


Can I have extra row in datatable not used for mapping

Is there a way that I can disregard 1 row in the table when generating a map and its Jenks values? We have the need to create separate tables and charts from the data we are mapping. There needs to be an extra row of data in the data-table for the tables/charts. But this won’t be included in the map or the Jenks values we’re using for a custom legend.


Most efficient method to join multiple fields in ArcGIS

I have as set of shapefiles with approximately 300000 polygons each, as well as corresponding point shapefiles and tables with approximately the same number of row (equal or less rows than the polygons FC). I need to join about 30 fields from the table of the point FC and from the stand alone table. It works, but it is quite slow. I've tried different method : AddJoin + claculateField, JoinField and arcpy.da.UpdateCursor. Join Field is the slowest and the Cursor is the fastest but I found it tricky to manage when the two tables are not matching exactly. Nevertheless, I was surprised that a spatial join is sometimes faster than a table join: this makes me think that I a missing something. I saw on a Website that I should work with layer instead of the FC, but I used "makeFeatureLayer" in my script and it didn't help much. Can I improve this with indexing, with other FC types or with another approach ? Note that this is part of a broader Python script with arcpy, so I need a solution that can be used within this script.


OGCServer as WMS application in Apache

I'm a really new in Mapnik OGC Server. Could you please help me with configuring OGCServer as WMS Server on Apache2? I use Windows platform and have already configured Mapnik,Python27. I've also successfully started OGCServer as python service. What should I do else so that WMS of OGCServer works as cgi or fastcgi or wsgi application? Let me way how to fix my issue or any step by step manuals how to configure Apache Server. Any ideas would be really useful for me!


when using arcpy.mapping.ListDataFrames, python.exe crash

Has anybody had this problem? this drive me crazy. what is wrong with my code?               import arcpy     mxd = arcpy.mapping.MapDocument(r"D:\Project\Project.mxd")     df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]      python crash at the last line


How to convince an organization to move towards open source GIS?

After spending the better part of the year on this site I am seeing how big the open source GIS community is. I used to just think it was ESRI and nothing else. That said, how do you begin to nudge or convince others that open source programs might be more beneficial than ESRI products? I'm not specifically asking about features in one vs the other, but actually trying to convince higher ups that a move to a similar product that would provide comparable results in terms of productivity, efficiency and end product. I guess what it boils down to in the end, is, at the end of the day, does the cost savings outweigh the time that it might take to get GIS staff acclimated to open source software?


Oracle Linestrings and MultiLinestrings

We have a very complicated setup for a project we are doing, and it involves 40+ developers and two oracle databases, one with sde and one only with oracle spatial and it's quite big. We'll here are my issues: In the SDE database, ArcSDE is creating my linear tables as Linestrings, while in the other Oracle database, we have Multilinestring tables. I'm not sure how that got by, and there is a ton of code built on top of the oracle only schema. How can I transform a simple linestring in a multilinestring? Here is what we've tried:   1. Using oracle append function to see if it returned a multi* geometric field. it does return a multi* geometry, but some parameters in the metadata are different   2. using sdelayer -o delete and sdelayer -o register to force a multilinestring. The procedure was fine, but even after new inserts, the oracle db is still seeing linestrings.   3. Merging two geometries inside arcgis. Still linesting geometry. Any thoughts? Hacks? Solutions? Changing the datamodel or code, specially code inside our oracle only db is an issue. Not possible. Thanks for your attention and insights! george


Importing DWF files into QGIS

I regularly receive files from a service authority in DWF format, which shows the location of their assets in there area. I note that there is no option to have the data in an alternate format as it is automatically generated. I cant seem to convert this file into an alternate format for import into QGIS (for example DXF). I note that it can be opened in Autodesk Design Review but this cannot convert the data (other than printing to a PDF). I would appreciate any guidance on how others incorporate this data format into their work. Thanks


Alternative cartographic software programs?

I am looking for an alternative software product for cartography other than ArcGIS and Adobe Illustrator. The graphic quality in ArcGIS outputs is substantially poor when compared to outputs in illustrator, yet illustrator is not a true cartographic software package. Any open source? Any extensions to ArcGIS or Illustrator that would enhance performance? Read a couple referrals to MAPublisher in other forums, yet this package is too expensive ($2000+). Thanks for any help you are able to provide.


How to style geojson polygon using L.mapbox.featureLayer()

I have this code:                   var SupIncend = L.mapbox.featureLayer().addTo(map);         var ArbolesAntes = L.mapbox.featureLayer().addTo(map);         var ArbolesX = L.mapbox.featureLayer().addTo(map);         var ArbolesDespues = L.mapbox.featureLayer().addTo(map);              L.control.layers({             'Superficie incendiada': SupIncend.loadURL('http://www.mapatraslasierra.com.ar/ProyectoForestal/geojson/zona_incendio_unida.geojson')         },{             'Superficie arbolada Agosto 2013': ArbolesAntes.loadURL('http://www.mapatraslasierra.com.ar/ProyectoForestal/geojson/ArbolesPerennes_AntesIncendio.geojson'),             'Superficie arbolada Octubre 2013': ArbolesX.loadURL('http://www.mapatraslasierra.com.ar/ProyectoForestal/geojson/ArbolesPerennes_Sobrevivantes.geojson'),             'Superficie arbolada Julio 2014': ArbolesDespues.loadURL('http://www.mapatraslasierra.com.ar/ProyectoForestal/geojson/ArbolesPerennes_Rebrote.geojson')             }).addTo(map);      It is working but I need to give differents colors and opacity to each layers (geojson representing polygons).


arcgis geodatabase replica between different versions

We have 2 geodatabase corresponding to 2 systems. There is a replica between these 2 geodatabase all use version 10. Recently there is a version upgrade in the parent geodatabase. It will be upgraded to version 10.2.2 (the child veision will remain version 10). I want to confirm can I create a replica between version 10 and version 10.2.2. Is there any potential danger? We have a scheduler to synchronize changes of the 2 geodatabase use the replica now. with the upgrade of the parent geodatabase, will the replica still work?


Max cell value from a stacked raster

How can I find a max cell value from a stacked raster.               Rmax <- maxValue(RAD1998.all[[1]])       works fine but               Rmax <- maxValue(RAD1998.all[[2]])       gives NA. Of course not in a stacked raster. Here is my code:               RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = ""))     list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = T, pattern = ".asc.bil"))      RAD1998.all <- stack(list.ras)


Speeding up arcpy's TextElement.delete()

I'm making tables to go on data driven pages. For each page there are a series of tables that are produced. The elements of the tables are made using the clone method. The problem is that it takes ages to delete the cloned elements before I go onto the next page. The code I'm using to delete the cloned elements is:               for elm in arcpy.mapping.ListLayoutElements(mxd, wildcard="_clone"):       elm.delete()      Does anyone know how I can speed up the deletion of these elements?


What is the interface/property in ArcObjects that corresponds to a layer's dataSource property in arcpy?

I am trying to get the "Catalog Path" of a Feature Layer's data source in the table of contents of ArcMap using ArcObjects. This is very straightforward in arcpy but seemingly not so much in ArcObjects. In arcpy, the Layer.dataSource property returns the "Catalog Path" of the layer. For example, for an SDE feature class in a feature dataset, it returns something like: `Database Connections\MyConnection.sde\SDE.MyFeatureDataset\SDE.MyFeatureClass` The help says that the `dataSource` property is just the `workspacePath` and `datasetName` properties combined, but this is not exactly true. For one, the `workspacePath` property expands the "Database Connections" shortcut, and secondly it does not include the feature dataset name. Thus `workspacePath` returns something like: `c:\Documents and Settings\user\Application Data\ESRI\Desktop10.0\ArcCatalog\MyConnection.sde` Incidentally this is also the value returned by IWorkspaceName.PathName. So is there an ArcObjects property that corresponds directly to the `Layer.dataSource` property, or is additional logic needed to build the same "Catalog Path" from the available ArcObjects interfaces?


Median value for a case within ArcGIS 10.1

Similar to this post (ArcGIS median value from a layer) I am looking for a toolbox to get the median value for a data field by a case. For example, say I have dataset with percent tree cover for each county for each state; I would like to summarize the data by calculating the median percent forest cover by state. The only reason the answer to the previous post doesn't work for my needs is that I am using ArcGIS 10.1 and the toolbox/script (link below) doesn't work in 10.1 I have attempted to update the script associated with the previous tool for use in 10.1, but haven't been successful. Is there an easy way to do this? The tool created by Ryan Dalton is available here: http://resources.arcgis.com/gallery/file/geoprocessing/details?entryID=A0CD1751-1422-2418-882E-001EE0DC0D35 ....


Scrollbars not visible in InfoWindow

I am having problems getting the scrollbars to show in my InfoWindow. Code is below. Any help would be appreciated. XAML:                <DataTemplate x:Key="MyFeatureLayerInfoWindowTemplate">                         <Grid>                             <Grid.ColumnDefinitions>                                 <ColumnDefinition Width="Auto"/>                                 <ColumnDefinition Width="*"/>                             </Grid.ColumnDefinitions>                             <ItemsControl ItemsSource="{Binding Keys}" Grid.Column="0" />                             <ItemsControl ItemsSource="{Binding Values}" Grid.Column="1" />                         </Grid>                     </DataTemplate>               <esri:InfoWindow x:Name="MyInfoWindow"                                  CornerRadius="20"                                  ContentTemplate="{StaticResource MyFeatureLayerInfoWindowTemplate}"                                  Map="{Binding ElementName=MyMap}"                                  MouseLeftButtonUp="MyInfoWindow_MouseLeftButtonUp"                                   ScrollViewer.VerticalScrollBarVisibility="Visible"                                      ScrollViewer.HorizontalScrollBarVisibility="Visible"                                      Background="{DynamicResource {x:Static SystemColors.WindowBrushKey}}"     MaxHeight="200" MaxWidth="200"/>      C#:               if (clickpoint != lastPoint)                 {                     MyInfoWindow.Anchor = clickpoint;                     MyInfoWindow.Content = e.Graphic.Attributes;                     MyInfoWindow.IsOpen = true;                     lastPoint = clickpoint;                 }      Picture of what I am seeing: ![enter image description here](http://i.stack.imgur.com/t52Ec.png) The content appears in the InfoWindow, but the scrollbars are nowhere to be found. Hopefully it is something simple, I feel like I have tried everything. Thanks in advance.


Text fields vs. Numeric fields

When designing geodatabases I have followed industry standards and used text and numeric fields. Some people have different ways of designing geodatabases and only use text fields. What are the schools of thought on text vs. numeric fields? Why would you use a text field when it is always a numeric value?


Spatial join between point files using FME

I have two files: the first has x,y,z information and is a (big) point cloud file, and second has only x,y information and is a (small) point shp file. The two files occupy the same area. For each of the x,y points in the shp file, I would like to find the nearest x,y point in the point cloud file and pass the corresponding point cloud z value to the attribute table of the shp file. I would also like to avoid converting the point cloud to any other format (due to file sizes) I am fairly new to FME and I'm not sure what is the right tool for the job. After some research, I tried the SpatialRelator but this seems to be the wrong tool (at least by itself) because it does not seem to support a "near" (within a certain distance) calculation. Maybe I'm too used to ArcGIS, but FME does not seem to have a "spatial join" type tool, or at least they are calling it something else...