title: Is there an open source alternative to ArcGIS for Desktop's Integrate tool?

text: I need to snap about 1000 lines in one layer to the lines in another layer (there are about 5000 lines in the second layer). The lines to be snapped are about 15 map units away from the base layer lines. It looks like ArcGIS could do this using the **Integrate** tool in the Data Manager toolbox. Is there a QGIS/GRASS/SAGA/other open source alternative to Integrate or will I need to crack open my Python skills? Thanks for any help/advice on this


title: How can I merge postcode sector polygons into larger territories in QGIS?

text: I'm new to GIS so patience might be required :-) I'm building a basic web-based mapping tool for business use, to be used as a demo to show that we need one but managers don't know it yet! So far I've built a HTML page displaying an OpenLayers OSM base layer. I've been supplied with some postcode sector polygons in shapefile format that I'm wanting to merge into our business territories (wards, districts etc.). Aims were to plot these as layers on top of OSM that the user can switch between etc. Using QGIS, i don't seem to be able to import or link to a table of our geography hierarchy and dissolve the polygons by these fields. Is this best achieved by importing CSV or linking to db etc. I'm a bit lost :-) our geo hierarchy:               Disrict, Ward, Neighbourhood, PC Sector     11, 11A, 11A15, AA1 1 etc.      Any advice would be gratefully received :-)


title: What are some of the OSX GIS tools

text: Most Opensource tools are cross platform; however, some GIS tools such as ArcGIS are platform specific (windows only). What are some of the GIS tools that are usable on the mac? **EDIT** Although thank you for putting the basic apps in your answers, my question specifies that Most if not all of opensource tools are cross platform, I was really hoping to get more obscure libraries, packages, and programs that I didn't know about on OSX similar to dotspatial and sharpmap for windows


title: Elevation graph from GPX file

text: I've got many GPX files in the following format:               <?xml version="1.0" encoding="utf-8"?>     <gpx xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xsd="http://www.w3.org/2001/XMLSchema" version="1.1" creator="c174a4ce-d611-43ef-b4d9-e4122bc8f2a4" xmlns="http://www.topografix.com/GPX/1/1">       <rte>         <name>Withens Clough, Cragg Vale and Stoodley Pike</name>         <rtept lat="53.705918666306154" lon="-2.02137437495311">           <ele>268</ele>         </rtept>         <rtept lat="53.705457119907116" lon="-2.0253609418514311">           <ele>293.8</ele>         </rtept>         ....       </rte>     </gpx>      I use this is a website to generate a route using OS Openspace maps. However, this doesn't create an elevation profile. **Is there a web service or something I can run on my web server that can generate GPX elevation graphs for embedding into webpages?** My web server runs PHP. I've found plenty of Google map ones but don't need the map, just the graph, but in the implementations I've looked at it seems like one depends on the other and don't want to compromise loading times by including two maps or remove the OS maps which are a lot better, in my view.


title: What is the best web editor and server for postgis?

text: I'm looking for a web editor (similar to openstreetmap) that works with an map Server (postGIS). This should permit edit and create new structures on map. What is the best option? Can you help me? Thank you for your time.


title: Geocortex using MS WF

text: I am going to be looking at the new Geocortex software soon, and wondered if anyone else was using it, what they think of it and if they had any online resources they could share. Thanks in advance


title: how to convert UTM northing and longitude to latitude

text: I am doing some work with the Military Grid Reference System (MGRS), which is based on Universal Transverse Mercator (UTM). I find that I need to convert from UTM northing and longitude (actually, delta angle between desired longitude and appropriate UTM central meridian) to latitude. I've searched for days and cannot find this conversion. I need an algorithm. Can anyone point me to a reference to such an algorithm, or to open source software that does this conversion?


title: How to use CQL_FILTER for the highest value in a layer field?

text: I need to query/filter my wms in openlayers to only see the polygons classified with the highest values in a given field. Would be something like:               wms_layer.mergeNewParams({'CQL_FILTER': "ND = max(ND)"});      But this is not working. Any ideas?


title: How to convert OSM polyline data to polygons?

text: I just now downloaded OSM data .I am having the data as polylines. I want that polylines as polygons is it possible in Qgis?below I also posted the screen shot of the polyline also.![enter image description here](http://i.stack.imgur.com/MPkJK.jpg)


title: Splitting polygons in OpenLayers is not working

text: I have some WFS layers in OpenLayers, containing Polygons. I try to make a tool to split these polygnos and save them to db.               //The layer which I want spliting work on     var contract_layer = map.getLayersByName("contracts");     //The spliting control     var split_control = new OpenLayers.Control.Split({            layer: contract_layer,            listeners: {                'beforesplit' : function(){                    alert("before");                }            }     });       and then I add this control to my map and make it active:               toolbarItems.push(new GeoExt.Action({            control: split_control,            map: map,            toggleGroup: 'map'        }));             map.addControl(split_control);        split_control.activate();      And now when I chose the tool, I can draw a line which ends with double-click, but that's it. nothing else.Nothing splits and nothing is sent to geoserver. What do you think I'm missing?


title: Popup for ArcGISDynamicMapServiceLayer point data

text: I was going through this example from ESRI Resource popup. I was able to integrate it to my polygon layer but it is failing for point and line layer.I am not sure but               query.spatialRelationship = esri.tasks.Query.SPATIAL_REL_INTERSECTS;      could be the reason.I tried the other combination of the task but was not successful if anyone has another example or if anyone could suggest something in existing example then it will be very helpful.I am using ArcGIS server 10.2,Javascript API 3.1.


title: GDAL - projection from pixels to geo coordinates

text: I am using GDAL library in C++. On my dataset, I have obtained geo transform with `GetGeoTransform` to GT array (GCP is NULL). Than I use this values for this equation (Source):               Xgeo = GT(0) + Xpixel*GT(1) + Yline*GT(2)     Ygeo = GT(3) + Xpixel*GT(4) + Yline*GT(5)      I also add 0.5 offset to `Xpixel` and `Yline`, to calculate coordinate in the "center" of pixel. If I put [0,0], I got correct result in lat/lon. However, if I change `Yline` to something other than 0, my longitude is messed up - instead of 44.0, I got -349.2... (latitude is still correct). What am I doing wrong ? **My data informations** \- origin, pixel size and top left are taken from GT array used in equation above; projection is `GetProjectionRef`               Projection is `GEOGCS["Coordinate System imported from GRIB file",DATUM["unknown",SPHEROID["Sphere",6371200,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]]'     Origin = (0,0) - GT(2), GT(4)     Pixel Size = (0.5,-33.268) - GT(1), GT(5)     Top Left = (5.75,66.634) - GT(0), GT(3)


title: Using a value from "iterate feature selection' to calculate a field

text: I'm trying to streamline a model that identifies all zones that intersect each parcel in our Local Government Area (eg. Parcel x has Zone B1, B2 and E4). I am trying to make use of iterators in ArcGIS10. The problem is in bold italics below; I have a Zones layer and a Parcels layer. I select all polygons of Zone A (using iterate feature selection) then do a select by location to get parcels that intersect those selected polygons. I now want to check if Field1 is empty, if it is, **_use calculate field to populate Field1 with the "Zone A" value used to select my zone in the first place._** I then continue with some simple if-then-else logic and the model then iterates with Zone B. ![Screen Capture of my model](http://i.stack.imgur.com/uMeNC.jpg) I built a working version of this in ArcGIS 9.3.1 that works but it is extremely cumbersome and making small changes to the process is very time-consuming. Any suggestions are much appreciated. P.S. I have no python experience and only a little VB knowledge but if writing an arcpy script is the only solution to my problem then I'm more than happy to start learning!


title: TINs: What advantage does unequal sampling of data points have over regular sampling?

text: Triangulated Irregular Networks. Under what conditions is one sampling method more appropriate over the other?


title: Split edges at intersections in PostGIS

text: I'm attempting to do a walking distance analysis using OpenStreetMap data (osm2pgrouting) and PostGIS and following underdark's blog posts, however I need to exclude residential and service roads that cross primary and secondary roads, as they are unsafe to cross. To do this I think my first step is to split edges at intersections. I'm guessing I'll want to use the ST_Split function, but I'm not sure how to put it all together in one SELECT statement. (A later question will be how to remove such crossing edges based on the highway=* tag, but I need to figure this one out first.)


title: QGIS 1.7 + MacOS 10.6 : georeferencing from map canvas bugs

text: I am experiencing issues with the gdal plugin to georeference a raster. It works fine when I enter the coordinates manually but it does not work when I try to link the plugin with the map canvas. Does anyone had this issue and was able to solve it? Thanks Guillaume QGIS 1.7 MacOS 10.6.7


title: Newbie PostGIS Geometry and Multilinestring Clarification

text: I'm having some trouble figuring out if I've got the right data loaded for PgRouting. Using the bash script included in PostGIS 2.0, I loaded Tiger2010 data for U.S. State California.   1. The edges table contains the_geom column, data type geometry. Using Underdark's example, it seems I need the roads in a Multilinestring format in order to begin generating routes. Can the data in the edges table be converted into type multilinestring? The table query is below.   2. I altered my geocoded location data to data type geography. PgRouting's shortest_path function needs the data in integer format. How do I convert the geography type point into an integer that shortest_path can use? Thanks for your patience.               CREATE TABLE tiger.edges     (       gid integer NOT NULL DEFAULT nextval('edges_gid_seq'::regclass),       statefp character varying(2),       countyfp character varying(3),       tlid numeric(10,0),       tfidl numeric(10,0),       tfidr numeric(10,0),       mtfcc character varying(5),       fullname character varying(100),       smid character varying(22),       lfromadd character varying(12),       ltoadd character varying(12),       rfromadd character varying(12),       rtoadd character varying(12),       zipl character varying(5),       zipr character varying(5),       featcat character varying(1),       hydroflg character varying(1),       railflg character varying(1),       roadflg character varying(1),       olfflg character varying(1),       passflg character varying(1),       divroad character varying(1),       exttyp character varying(1),       ttyp character varying(1),       deckedroad character varying(1),       artpath character varying(1),       persist character varying(1),       gcseflg character varying(1),       offsetl character varying(1),       offsetr character varying(1),       tnidf numeric(10,0),       tnidt numeric(10,0),       the_geom geometry,       CONSTRAINT edges_pkey PRIMARY KEY (gid),       CONSTRAINT enforce_dims_the_geom CHECK (st_ndims(the_geom) = 2),       CONSTRAINT enforce_geotype_the_geom CHECK (geometrytype(the_geom) = 'MULTILINESTRING'::text OR the_geom IS NULL),       CONSTRAINT enforce_srid_the_geom CHECK (st_srid(the_geom) = 4269)     )     WITH (       OIDS=FALSE     );     ALTER TABLE tiger.edges OWNER TO postgres;     GRANT ALL ON TABLE tiger.edges TO postgres;     GRANT ALL ON TABLE tiger.edges TO gis_group;          -- Index: tiger.idx_edges_tlid          -- DROP INDEX tiger.idx_edges_tlid;          CREATE INDEX idx_edges_tlid       ON tiger.edges       USING btree       (tlid);          -- Index: tiger.idx_tiger_edges_countyfp          -- DROP INDEX tiger.idx_tiger_edges_countyfp;          CREATE INDEX idx_tiger_edges_countyfp       ON tiger.edges       USING btree       (countyfp);          -- Index: tiger.idx_tiger_edges_tfidl          -- DROP INDEX tiger.idx_tiger_edges_tfidl;          CREATE INDEX idx_tiger_edges_tfidl       ON tiger.edges       USING btree       (tfidl);          -- Index: tiger.idx_tiger_edges_tfidr          -- DROP INDEX tiger.idx_tiger_edges_tfidr;          CREATE INDEX idx_tiger_edges_tfidr       ON tiger.edges       USING btree       (tfidr);          -- Index: tiger.tiger_edges_the_geom_gist          -- DROP INDEX tiger.tiger_edges_the_geom_gist;          CREATE INDEX tiger_edges_the_geom_gist       ON tiger.edges       USING gist       (the_geom);


title: OpenLayers: how to save added point in server side

text: * I'm new in OpenLayers.   * I have simple application, which use a picture as map,   * on this picture users can add, move, or delte point, and then save.   * I have implemented adding, moving, deleting functionalities, only except save.   * codes are in below : codes here   * anyone who knows how to save change on server side? save to database, file?


title: How to validate Personal Geodatabase against XML schema?

text: Is their any tool in arcgis model builder which can validate a personal geodatabase and a XML schema to give out an report saying if the gdb is validated or not ...


title: How to use SLD Select control?

text: I have a WMS layer and i want to get fid of polygon in it by click it. How can i do it?   I know that i can create WFS layer from this WMS and get feature attributes from layer or store. But its possible do this use just WMS? I dont want use WFS Layer couse it may contain ~30k features. But if i will use WFS Layer with visible=false selection control will be works? I find SLDSelect control its can help me?   Maybe its possible to do like in gxp library?


title: Invalid Geometry Error after ST_ConvexHull

text: I've a Postgis table of around 80,000 points and want to be able to populate a second table summarizing the data with convex hulls based on points with a common attribute. For example, the data contains the telephone exchange which each point is associated with - I want to create a convex hull for all points in each telephone exchange. It ought to create around 1000 overlapping polygons. I tried using the normally excellent QGis ftools convex hull plugin with the option to create hulls based on an input field but my computer spent 6-hours maximising 3 of the four cores before crashing QGis. And I've tried doing it directly in sql: > INSERT INTO new_table(exchange_name, the_geom)   >  SELECT exchange_name, ST_ConvexHull(ST_Collect(the_geom)) As the_geom   >  FROM first_table   >  GROUP BY exchange_name But this complains about invalid geometry - I've tried geometry constraints of "MULTIPOLYGON" and generic "GEOMETRY" so not sure why this should happen - the source table loads perfectly in QGis. The Grass "v.hull" tool only creates a convex hull for all points and not based on an attribute. Any tips or suggestions would be very welcome!


title: Why is shapefile layer not visible on WMS layer in QGIS?

text: In my QGIS project I have added a WMS layer for which the CRS is WGS84 EPSG:4326. Now I want to show my shapefile layer which has the same CRS. QGIS shows the layer separately but does not show with the above WMS layer. Can any one help me to understand why? I am using Quantum GIS 1.7.0 I made up shape file for my project for a city area in which point locations are there.with some attribute data. I have directly drag and drop shape file in q GIS. is it ok or i made some wrong thing?please guide me..


title: Batch Conversion software for latitude and longitude into shapefile format

text: I have a huge list of latitude and longitude and I want them to be converted into a shapefile format using some batch processing.


title: Editing QGIS point layer: retrieve attributes of intersecting line layer

text: I'm creating a custom editing form using QGIS 2.0 and would like to modify the mouse event when creating a point. I've gotten a start by reading this thread: How to programatically check for a mouse click in QGIS In fact, the code below is copied and modified directly from it, so thank you if you helped me with this code snippet. My goal: When creating a new point, I would like to retreive feature attributes of another (line) layer that intersects with the mouse click, using a rectangle filter. The mouse click opens the form when creating the point, while automatically retrieving data from another line layer intersecting with the rectangle buffer. This is in the interest of comparing data in these two layers that are supposed to correspond. Following the lead of the mentioned thread, i'm attempting to make use of the following code without creating a separate Map Canvas. To put it to use, what would I provide, if anything, as the Python Init function under Layer Properties > Fields? Currently I have mypyfile.formOpen as I've added python logic as detailed here: http://nathanw.net/2011/09/05/qgis-tips-custom- feature-forms-with-python-logic/ I may be going about this in the wrong way, so i'm open to suggestions to direct me towards creating this data collection tool.               class PointTool(QgsMapTool):         def __init__(self, canvas):             QgsMapTool.__init__(self, canvas)             self.canvas = self.iface.mapCanvas()             self.iface = iface             self.pointTool = QgsMapToolEmitPoint(self.canvas)              def run(self):             self.canvas.setMapTool(self.pointTool)             QObject.connect(self.pointTool, SIGNAL("canvasClicked(const QgsPoint, Qt::MouseButton)", self.selectFeature))              def selectFeature(self, point):             #resLayer previously defined             if not resLayer or resLayer.type() != QgsMapLayer.VectorLayer():                 QMessageBox.warning(None, "No!", "Select a vectorlayer")                 return                  width = self.iface.mapCanvas().mapUnitsPerPixel()*2             rect = QgsRectangle(point.x() - width, point.y() - width, point.x() + width, point.y() + width)             rect = self.iface.mapCanvas().mapRenderer().mapToLayerCoordinates(layer, rect)                  resLayer.select([], rect)             feat = QgsFeature()             ids = []             while resLayer.nextFeature(feat):                 ids.append(feat.id())             #selecting the feature or features that intersect with mouse click buffer             resLayer.setSelectedFeatures(ids)          #code continues to guide user data entry using python logic


title: what is the relevance of cosine in determining map resolution?

text: "At any latitude and zoom level, you can determine the scale by using the following equation"               Map resolution = 156543.04 meters/pixel * cos(latitude) / (2 ^ zoomlevel)      is a quote from (http://msdn.microsoft.com/en-us/library/aa940990.aspx) Can anyone explain to me the relevance of the cos(latitude) part? I am having difficulties visualising this.


title: IToolbar uncheck items - ArcGIS

text: I currently have an app with two menu bars. One being the AxToolbarControl in a UserControl and one in a separate menu bar containing in-house stuff. I currently have the issue where if I select something in the AxToolBarControl and then select something in the other menu-bar the item in the AxControl is still checked. This obviously results in two commands being executed on the mousedown event handler. I've looked at the IControl and IToolbar but it looks like .Checked is a ReadOnly. Any ideas?


title: How to collect data via WFS and Popup

text: I have a map with a WFS-Layer and I am able to add Polygons to my data stored in a PostGIS-database, served with Geoserver (Openlayers, Geoserver 2.1.4, PostGIS 2.0 for PostgreSQL 9.1).) How can I add a description to the polygon? E.g., after drawing a polygon I want to specify the crop that are grown there, i.e. insert information into the Postgis-table-column "crop". I found this example: http://gis.ibbeck.de/ginfo/apps/OLExamples/OL26/examples/styles_unique_with_group_wfs.html. The Popup type "edit name and description" would be perfect for my means, unfortunately I can't work out how it works as i have a really hard time understanding its code, greenhorn that I am...   * How is it done?   * Or is there an easier example/ way of how to use a popup to edit your attributes? Thanks so much, I appreciate any help I can get! * * * Here's my WFS-Layer and the tools used so far. The desired column to edit/add Information in the postgis-table "test" is "crop".               function init() {     Save-strategy     var saveStrategy = new OpenLayers.Strategy.Save();          //empty map, bounds are test-layer bounds (EPSG:32647)       map = new OpenLayers.Map({         div: "map",         allOverlays: true,         maxExtent: new OpenLayers.Bounds(         653237.69439077,1519879.063165,655229.57939001,1520825.6733868         )         });          //WFS-Layer Test= editable data     var test = new OpenLayers.Layer.Vector("Editable Features", {         strategies: [new OpenLayers.Strategy.Fixed(), saveStrategy],             protocol: new OpenLayers.Protocol.WFS({         url: "http://..../wfs",         featurePrefix: 'testkf',         featureNS: "http://.../testkf",         extractAttributes: true,         featureType: "test",         geometryName: "geom",         })     });          map.addLayer(test);          //Toolbar:         var panel = new OpenLayers.Control.Panel(         {'displayClass': 'customEditingToolbar'}         );              var draw = new OpenLayers.Control.DrawFeature(         test, OpenLayers.Handler.Polygon,         {         title: "Draw Feature",         displayClass: "olControlDrawFeaturePolygon",         multi: true         });              var save = new OpenLayers.Control.Button({         title: "Save Changes",         trigger: function() {         if(edit.feature) {         edit.selectControl.unselectAll();         }         saveStrategy.save();         },         displayClass: "olControlSaveFeatures"         });               panel.addControls([save, draw]);     map.addControl(panel);       map.addControl(new OpenLayers.Control.LayerSwitcher());     map.zoomToMaxExtent();     }      EDIT: Hi, so nobody knows anything about this topic? Or is there another problem with my question? I found another example, however, this is way too complex for me to understand yet: http://dev4.mapgears.com/bdga/bdgaWFS-T.html This is exactly what I am looking for. Can anyone break it down to me what the steps are to achieve this? Or is there any way to improve my question? Thanks for your help!


title: Doubt in Field Calculator

text: I am working in an project.I have a point Shapefile its projection is **NAD83(HARN)Pennslyvaniasouth(ftUS)*** _EPSG3365_ * when I measure its distance(between two points) it gives me the answer in feet .yes its acceptable for me.I also calculated its x-coordinate and y-coordinate for the points in field Calculator it displays me the answer in feet.This is the problem.In my company they are insisting me to show the LAT-LONG value in Decimal Degrees.Is it possible to show the Lat-long in Decimal degrees in QGIS?


title: Looking for free program, similar to 3D Analyst, to visualize tree stands in 3D

text: Heyho, I want to make a 3D stand visualization of trees in a specific area. I need something similar in functionality to ArcGIS 3D Analyst, but I don't have a license for that plugin anymore. The only alternative I have found is the Stand Visualization System from the USDA Forest Service, which works quite well, but has no GIS functionality. I want to visualize the height and clustering of geographically referenced trees. There have to be different models for the tree species (spruce, beech, ... ). EDIT: For all people searching for more tools -> here is another posting on gis.stack


title: View information raster data in postgresql

text: I download and use raster2pgsql.exe import tiff image into postgresql. Can i use the query to retrieve the information of the image does not? Information such as: date of capture, capture time.


title: Add time restriction in ArcGIS Network Analyst

text: I'm trying to use train tracks as a network in ArcGIS Network Analyst. We are trying to create a model for security agents on board our trains and we'd like to optimize their daily "route". For now, it's all made by hand. My approach is to take our gtfs data and create a network using stop times to determine when the train segment is "traversable". When there's a train on the segment, it's open, if not, it's restricted. The problem is I can't set a restriction based on time in Network Analyst (or wasn't able to). I would like to be able to say by a formula or script than you can take that segment only from 7:03 AM to 7:10 AM for example. Is it possible? Thanks!


title: ArcGIS Javascript API zoom out?

text: I am using the ArcGIS API for Javascript and am wondering if it comes with functionality for zooming out using a rectangle, similar to the behavior in ArcGIS Desktop? Specifically, the user draws a rectangle and it zooms out based on that rectangle. I can create the functionality myself, but was wondering if it already existed in the API...


title: Data naming convention - starting with a number

text: My company uses ArcGIS and has a project and data file naming standards in place and (for the most part) followed. Something that has always bothered me about he naming standards is that it mandates starting all project and data file names with the project number - an eight digit _number_. I've always held the belief that naming GIS files starting with numbers is a bad thing, and have had (especially with GRIDS) processes fail because of the file name. I'm looking to amend the corporate standards to drop the project number requirement, however I cannot find much in the way of documentation on why "numbers as a first character" in filename is a bad thing. Can anyone point me in the right direction as far as resources to support this argument? Thanks,


title: How to divide a polyline by intersecting isochrones

text: I have two shapefiles, the first representing a path (a polyline with a single category), the second representing isochrons from the same origin (point 0) which intersect the above mentioned path. I wonder if it is possible to divide the polyline (path) by the intersecting isochrones (at point 1, point 2, etc.), thus obtaining a segmented polyline with different categories. QGIS has an useful tool called "line intersections" (vector - analysis tools - line intersections), but it only creates a further shapefile with points at the intersections of the two layers. Can someone suggest a solution preferably in QGIS or GRASS GIS.


title: Importing AUTOCAD .dwg to Arcgis Symbols

text: My client has some symbol drawn in autocad.. they are saved in .dwg files. How can I import them into ArcGIS symbols? I have FME 2013, I am not sure if that helps.


title: Trigger a python script to run after a row is added to the feature class

text: I have a Python script which runs on a SDE feature class daily (as a scheduled task). I would actually like this script to run whenever a new row is added to the feature class (data is captured in the field using a mobile device and posted to the server). I was considering using a trigger as mentioned in this question as my data is not versioned, but I do not want to overcomplicate things if there is an inbuilt SDE tool which can achieve this or something similar. The data is currently stored in ArcSDE 10/SQL Server 2008, but I do have ArcSDE 10.1/SQL Server 2012 available. Currently when the script runs, it checks for `NULL` values in a certain field. If it finds it, it selects those rows and performs several geoprocessing operations on them. It then populates that field with a value to show that everything has been processed. While I have a lot of experience editing/reconciling/posting/versioning SDE data, I only have limited experience with the database/admin side of it (assigning privileges, basic dabbling in SQL management studio), with a basic understanding of the underlying architecture.


title: How do I create a custom map with subtle details for a localized place?

text: I've worked with google maps for android previously. But, now I need to build **custom maps for a place such as a home** , showing bedrooms, kitchen and living room, **or a factory** , displaying all the departments. I would like to use these maps in an android application. How do I achieve this objective?


title: Performance of Postgres VIEW published in GeoServer

text: Does anyone know what will be the performance difference between publishing a spatial table directly vs. creating a view of the same table and publishing it in GeoServer? Does the spatial index on a table still work on a view derived from it? Thanks a lot!


title: How to correctly render a raster Layer to GTiFF with QGIS python console?

text: I am trying to render a map from a raster layer. As recommended in pyqgis cookbook and here(osgeo.org) I use `QgsMapRenderer()`. I am able to produce a picture but it has not appropriate extent and geo-spatial information. Using QGIS GUI, the rendering works well, (Save As, Rendered image,..): ![enter image description here](http://i.stack.imgur.com/0xWga.png) This is a global map (whole earth) 350x300 displayed in Gimp. The tiff can be loaded to the Geoserver without problem: ![enter image description here](http://i.stack.imgur.com/0h1o1.png) * * * However, using following code in the QGIS python console the picture looks different and the .tif appears not to have any geo-spacial data:               from qgis.core import *     from qgis.gui import *     from qgis.utils import *     from PyQt4.QtCore import *     from PyQt4.QtGui import QImage, QColor, QPainter          modeled = "d:/GISdata/THRASH/layer.tif"           fileInfo = QFileInfo(modeled)     baseName = fileInfo.baseName()     modlayer = QgsRasterLayer(modeled, baseName)     QgsMapLayerRegistry.instance().addMapLayer(modlayer)          # create image     img = QImage(QSize(modlayer.width(),modlayer.height()),        QImage.Format_ARGB32_Premultiplied)          # set image's background color     color = QColor(255,155,200)     img.fill(color.rgb())          # create painter and start to pain to the img     p = QPainter(img)          # create Map Renderer     render = QgsMapRenderer()          # set layer set     lst = [ modlayer.id() ]  # add ID of every layer     render.setLayerSet(lst)           # set extent     rect = QgsRectangle(modlayer.extent())      render.setExtent(rect)          # set output size     render.setOutputSize(img.size(), img.logicalDpiX())          # do the rendering     render.render(p)     p.end()          # save image     img.save("C:/GISdatafolder/render.tif","tif")      ![enter image description here](http://i.stack.imgur.com/LxDAK.png) ![enter image description here](http://i.stack.imgur.com/RVyvz.png) the text is: > Could not list layers for this store, an error occurred retrieving them: > Failed to create reader from file://C:\GIS\render.tif and hints null Also, when I load the .tif back to QGIS it says that CRS was undefined and default WGS84 was applied. How can I produce the picture with desirable extent and with information about CRS?


title: Creating a fully custom map using OSM infrastructure

text: I would like to make a fully custom map (fictional, have to start from scratch basically) and use the OSM tools to make following things possible:   * **crowd editing** \- the source data which is used to generate the tiles (and in the end a large continuous map) will be publicly available, anyone can make a change to it (I don't require a nice web interface like OSM has, simple github repository will suffice)   * **custom rendering** \- I know that this is possible using custom tiles and such so this shouldn't be an issue. Before I dive into this I would like to know if this is possible. I think tools like JOSM should be good enough for drawing the maps, then I can define my own tiles and use one of the rendering tools to display the map. Are there any pitfalls I'm missing? Or maybe a better way to achieve this?


title: AutoDesk Infrastructure Studio cannot access Map Server (Windows + IIS)

text: After restarting the entire server I continue to get this error when accessing any mapagent service (e.g. `http://myserver/mapserver2013/mapagent/mapagent.fcgi?SERVICE=WFS&request=getcapabilities`). Although the error does a good job at suggesting the problem there is no help for how to actually fix the problem. I suspect Map Server is not running...but how do I confirm and fix that? > Problem: Cannot establish connection to the Infrastructure Map Server. > > Possible Causes: -Your session has been idle for too long; -Infrastructure > Map Server is no longer responding; -Internet connection problems. > > Solution: Restart your application / session, restart the Infrastructure Map > Server service, or contact with the server administrator. Also, when I perform Site Authentication in Audodesk Infrastructure Studio I get an error saying: > Unable to connect to Map Server. Please ensure that the server process is > running and firewall rules have been set correctly. I don't know what the server process is called. I see nothing related to autodesk/mapserver/mapguide under services.


title: How to listen for ESRI drawing events using ESRI drawing tools in ArcObjects?

text: I am trying to make use of the built-in ArcObjects drawing tools (ie ControlsNewCircleTool) - but I need to listen for the events so I can add the drawing graphics to the database in our application. Is there an event or way to listen for when the user has drawn a new graphic and also modified it?


title: OpenLayers plugin and QGIS 1.9: Error after install

text: After I install the OpenLayer Plugin (1.0.0) on my QGIS 1.9 this error appears:               Konnte Erweiterung openlayers_plugin nicht laden wegen eines Fehler beim Aufruf seiner initGui()-Methode          Traceback (most recent call last):       File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 198, in startPlugin         plugins[packageName].initGui()       File "/home/map/.qgis//python/plugins/openlayers_plugin/openlayers_plugin.py", line 160, in initGui         if not self.__setCoordRSGoogle():       File "/home/map/.qgis//python/plugins/openlayers_plugin/openlayers_plugin.py", line 211, in __setCoordRSGoogle         if not self.__coordRSGoogle.createFromEpsg(idEpsgRSGoogle):     AttributeError: 'QgsCoordinateReferenceSystem' object has no attribute 'createFromEpsg'          Python-Version:     2.7.3 (default, Sep 26 2012, 21:53:55)      [GCC 4.7.2]               QGIS-Version:     1.9.0-Master Master, exported          Python-Pfad: ['/home/map/.qgis//python/plugins/sextante', '/usr/share/qgis/python', '/home/map/.qgis//python', '/home/map/.qgis//python/plugins', '/usr/share/qgis/python/plugins', '/usr/lib/python2.7', '/usr/lib/python2.7/plat-linux2', '/usr/lib/python2.7/lib-tk', '/usr/lib/python2.7/lib-old', '/usr/lib/python2.7/lib-dynload', '/usr/local/lib/python2.7/dist-packages', '/usr/lib/python2.7/dist-packages', '/usr/lib/python2.7/dist-packages/PIL', '/usr/lib/python2.7/dist-packages/gst-0.10', '/usr/lib/python2.7/dist-packages/gtk-2.0', '/usr/lib/pymodules/python2.7', '/usr/lib/python2.7/dist-packages/ubuntu-sso-client', '/usr/lib/python2.7/dist-packages/wx-2.8-gtk2-unicode', '/home/map/.qgis/python/plugins/qgis2kml', '/home/map/.qgis/python/plugins/qgis2kml', '/usr/share/qgis/python', '/usr/share/qgis/python/plugins/fTools/tools']      Any ideas? Thanks! Martin


title: Data View/Layout View Changes - ArcMap 10

text: I have code that will hide a certain set of points from the view using an ExclusionSet, and then switch to layout view for printing. I noticed in my 9.1 version, when after I am done printing, and I switch back to data view, and the points I had before running the code are shown. This was my intention. In 10 on the other hand, the code runs the same, but when I return to data view, the feature exclusions are still in-tact. I have to manually remove them via the display properties of the feature layer. Is there anyway to have separate properties for both views? Whether it be code or manual? Thanks!


title: How to understand output of ArcGIS System Monitor?

text: We had some performance issues with our enterprise ArcGIS server application, in this regard I have installed **ArcGIS System Monitor v1.0.9**. After installing and configuring the site, they are several variables (below). I would like to know how these variables effect the performance and their definitions. Please throw some light on this. ARCGIS Percent Processor Time   ARCGIS Working Set GB   ARCGIS Active Working Set GB   ARCGIS Instance Count   ARCGIS Active Instance Count   ARCGIS Idle Instance Count   ARCGIS Virtual GB   ARCGIS Private GB   ARCGIS Page File GB   ARCGIS Page Faults/sec   ArcSOC Percent Processor Time   ArcSOC Working Set GB   ArcSOC Active Working Set GB   ArcSOC Active Instance Count


title: FME - How to extract the rightmost folder in a path

text: I have an FME workbench which reads all of the TAB files in a folder and all of the subfolders and writes the paths of those TAB files as well as their type (NATIVE, RASTER, etc) to an Excel file. The folder is a published parameter. What I would like to do is name the Excel file the same as the rightmost folder (I'm not sure what you'd call this) in the parameter. For example, when the workbench is run, the user selects this folder:               C:\Temp\ABC\DEF\GHI      The excel file would be GHI.xlsx or if the user selects this folder:               C:\Folder\ABC\DEF123      The excel file would be DEF123.xlsx I have tried many different ways to get this folder name. I tried creating a new parameter that builds the name but no luck. I've also tried setting the parameter as an attribute using the parameterFetcher then using the SubstringExtractor. However, I can't seem to get the start and end index set properly. This is what I've got in my substringExtractor: ![enter image description here](http://i.stack.imgur.com/VjnNZ.png)


title: how to group interconnected edges inside a table?

text: Given a table and a random set of edges with source and target columns, how can I find out which edges connect with each other and then group them together using PostgreSQL?


title: Python module "shapely.wkb" missing

text: I am trying to follow a tutorial from here http://qgistips.spatialthoughts.com/2012/12/find-neighbor-polygons-in-layer- in-qgis.html#more I have downloaded the relevant python script and edited it as required. When I try to run the script using the Python console in QGIS >Plugins>Python Console I get the following error message. ImportError: No module named shapely.wkb Start of relevant script               from PyQt4.QtCore import *     from shapely.wkb import loads     # Replace the below value with the field containing name or id of the feature     # For example, if your field is called name then change the line below to     # name_filed = 'name'      The module shapely.wkb is called on the second line of the script so the failure appears to be occurring as the second line of code is executed. I have searched for a Python module with that name in Plugins>Fetch Python Modules Is this something that needs to be installed separately? From what I have read it should be installed with an installation of PostGis which I though got loaded when I used StackBuilder to install PostGreSQL. Would installing installing Python 2.7.3 on my MAC Book Pro help?


title: Why is the Frequency Tool in ArcGIS so Fast and Efficient?

text: Why is the Frequency tool in ArcGIS is so fast and efficient? Is it possible to make my own Frequency Tool with similar performance (in python)? I tried to do some based on python dictionary and searchCursor but when I have more than 100,000 records to analyse it takes a long time.


title: Converting Between Similar Classes in ArcObjects and WPF Runtime Libraries (C#)?

text: There seems to be an "impedance mismatch" (if I can borrow the term) between similar classes in ArcObjects and the WPF runtime libraries. For example, there is a ESRI.ArcGIS.Client.Geometry.Polyline in the WPF runtime and a ESRI.ArcGIS.Geometry.IPolyline (and corresponding PolylineClass) in ArcObjects. I want to be able to do something like this (or at least similar):               var poly = new ESRI.ArcGIS.Client.Geometry.Polyline();     var ipoly = new ESRI.ArcGIS.Geometry.PolylineClass();          // stuff happens          poly = ipoly;      Is there an easy way to convert between the two? Or am I going to have to write some deep copy methods?


title: Queries on virtual network in Spatialite only make errors

text: I've installed Spatialite 2.4 in Ubuntu. I try to convert an osm file into a roads network using this tutorial and make searches on it, but the queries only produce "no column" errors.               $ spatialite_osm -o data/mosobl.osm -d test.sqlite -T roads -m     SQLite version: 3.6.22     SpatiaLite version: 2.4.0     using IN-MEMORY database     Loading OSM nodes ... wait please ...     Loading OSM nodes ... wait please ...     Loaded 3261086 OSM nodes     Verifying OSM ways ... wait please ...     Verified 108981 OSM ways     Disambiguating OSM nodes ... wait please ...     Found 112 duplicate OSM nodes - fixed !!!     Loading network ARCs ... wait please ...     Loaded 227809 network ARCs     Dropping temporary table 'osm_tmp_nodes' ... wait please ...     Dropped table 'osm_tmp_nodes'     Dropping index 'from_to' ... wait please ...     Dropped index 'from_to'     exporting IN_MEMORY database ... wait please ...     IN_MEMORY database succesfully exported     VACUUMing the DB ... wait please ...     All done: OSM graph was succesfully loaded      Checking the table structure:               $ spatialite test.sqlite      SpatiaLite version ..: 2.4.0    Supported Extensions:     - 'VirtualShape'    [direct Shapefile access]     - 'VirtualDbf'      [direct DBF access]     - 'VirtualText'     [direct CSV/TXT access]     - 'VirtualNetwork'  [Dijkstra shortest path]     - 'RTree'       [Spatial Index - R*Tree]     - 'MbrCache'        [Spatial Index - MBR cache]     - 'VirtualFDO'      [FDO-OGR interoperability]     - 'SpatiaLite'      [Spatial SQL - OGC]     ...          spatialite> .schema roads     CREATE TABLE "roads" (     id INTEGER NOT NULL PRIMARY KEY AUTOINCREMENT,     osm_id INTEGER NOT NULL,     class TEXT NOT NULL,     node_from INTEGER NOT NULL,     node_to INTEGER NOT NULL,     name TEXT NOT NULL,     oneway_fromto INTEGER NOT NULL,     oneway_tofrom INTEGER NOT NULL,     length DOUBLE NOT NULL,     cost DOUBLE NOT NULL, "geometry" LINESTRING);     CREATE TRIGGER "ggi_roads_geometry" BEFORE INSERT ON "roads"     FOR EACH ROW BEGIN     SELECT RAISE(ROLLBACK, 'roads.geometry violates Geometry constraint [geom-type or SRID not allowed]')     WHERE (SELECT type FROM geometry_columns     WHERE f_table_name = 'roads' AND f_geometry_column = 'geometry'     AND GeometryConstraints(NEW."geometry", type, srid, 'XY') = 1) IS NULL;     END;     CREATE TRIGGER "ggu_roads_geometry" BEFORE UPDATE ON "roads"     FOR EACH ROW BEGIN     SELECT RAISE(ROLLBACK, 'roads.geometry violates Geometry constraint [geom-type or SRID not allowed]')     WHERE (SELECT type FROM geometry_columns     WHERE f_table_name = 'roads' AND f_geometry_column = 'geometry'     AND GeometryConstraints(NEW."geometry", type, srid, 'XY') = 1) IS NULL;     END;      Verifying the structure:                   $ spatialite_network -d test.sqlite -T roads -f node_from -t node_to -g geometry -c cost -o roads_net     SQLite version: 3.6.22     SpatiaLite version: 2.4.0     Step   I - checking for table and columns existence          spatialite-network          ==================================================================        SpatiaLite db: test.sqlite     validating table: roads          columns layout     ==================================================================     FromNode: node_from       ToNode: node_to         Cost: cost         Name: *unused*     Geometry: geometry          assuming arcs to be BIDIRECTIONAL          NETWORK-DATA table creation required: 'roads_net'     Overwrite not allowed if table already exists     ==================================================================          Step  II - checking value types consistency     Step III - checking topologic consistency     Step  IV - final evaluation          Statistics     ==================================================================         # Arcs : 455042         # Nodes: 178470         Node max  incoming arcs: 7         Node max outcoming arcs: 7         # Nodes   cardinality=1: 40228 [terminal nodes]         # Nodes   cardinality=2: 16388 [meaningless, pass-through]     ==================================================================               OK: network passed validation         you can apply this configuration to build a valid VirtualNetwork     OK: validation passed               OK: NETWORK-DATA table 'roads_net' successfully created     OK: table 'roads_net' successfully created      Now I log into spatialite and try to get some routing. Should it work from the original roads table?               spatialite> select * from roads where node_from=4870889 and node_to=4870985;     (nothing)      No. Or from roads_net?               spatialite> select * from roads_net where node_from=4870889 and node_to=4870985;     SQL error: no such column: node_from      Apparently, I can only get a single way with this search: take `node_from` and `node_to` and put them in the query to `roads`. What should I do to run a route query? PS: Here's what tables I have in the db. The tutorial mentions roads_net_data which is absent here.               spatialite> .tables     geom_cols_ref_sys       roads                   views_geometry_columns     geometry_columns        roads_net               virts_geometry_columns     geometry_columns_auth   spatial_ref_sys


title: How do you sort a Standalone table that's already visible within ArcMap?

text: I'm using ITableWindow2.FindViaStandaloneTable--that works. I can get the ITableWindow for the stand alone table and the table from that. Then I try to sort using an ITableSort. No exceptions, but the table doesn't sort. Is ITableSort the right approach or is there a different approach to sorting within a Table Window?               ITableSort sort= new TableSortClass();     sort.Fields = "MDATE";     sort.Ascending["MDATE"] = true;     sort.Table = standAloneTable.Table;     sort.Sort(null);      The goal is to programmatically sort a field. The way an end-user would manually right click a column on an attribute table and sort ascending or descending.


title: How to handle multiple layers when importing a polygon shapefile into GRASS?

text: When I import a shapefile into GRASS, it seems to arrive as 3 layers. I.e. they appear in the qGIS toolbox browser as:               myMap 0_polygon, myMap 1_polygon, myMap 2_polygon      Is there a way to merge the layers into one? I have run:               v.clean tool=snap,break,rmdupl,bpol thresh=10 input=myMap output=myMapCleaned      And searching on this topic suggests that I should now run v.patch to merge the layers. But how do I specify the individual layers in the 'input=...' string? Thanks.


title: Print Maps - OpenLayers, GeoServer, GeoExt, MapFish

text: Can anyone point me towards a clear and unambiguous description of how to install and code a very simple map print function for OpenLayers/GeoServer? I am using GeoServer 2.5.2 and OpenLayers 2.12 on a Windows machine. Javascript handles the OpenLayers map and I send a lot of stuff to VB.NET for heavier backend work. As of now I have added the GeoServer 2.5 printing plugin, instructions here: http://docs.geoserver.org/stable/en/user/community/printing/ So far so good, everything works. After this point there are two options, MapFish Print Module or GeoExt. It appears that GeoExt requires Extjs which I would have to license (my project is not opensource) and I don't want to deal with licensing or license fees. Can GeoExt be used alone with OpenLayers/GeoServer to print? The other option is MapFish - MapFish Print Module - which seems like a much better alternative but I am completely befuddled on how to integrate it into my existing GeoServer/OpenLayers/JavScript/Vb.NET application. For instance it says to build the file. Build it where? In windows command line? In a java editor? Any help would be appreciated! Thanks


title: Google Maps API v3 vs v2 (JavaScript): overlay problems

text: I have been using the Gmap API (v2) for dynamically displaying some data:               oMap.map = new GMap2(document.getElementById("map"));        oMap.geoXml = new GGeoXml('http://darudar.org/var/files/workaroud/gmap/data.php?num=30');      I changed to v3 of the API, using the following code:               var georssLayer = new google.maps.KmlLayer('http://darudar.org/var/files/workaroud/gmap/data.php?num=30');     georssLayer.setMap(oMap.map);      Unfortunately, after the change the data markers do not show correctly or do not show at all. With v2, zooming and scaling allows a kml layer to get a BBOX array and to redraw the markers correctly. Does anyone know how to recreate this functionality with v3 of the API ? PS: Sorry for my eng...


title: Render a map showing lots of individual traces/routes simultaneously

text: Google's "My places" is great for showing lots of traces (in my case, paths representing bike rides I've done). But it only shows 20 or so per page. I'd really like a map (a statically rendered .png would be ok - but prefer dynamic online) with 1-200 traces highlighted over a base map. I'm happy to use online or offline tools (preferably Windows, Mac ok, last resort Linux) to achieve this, but the fewer steps and unfamiliar technologies the better. (I'm not really familiar with any offline GIS tools.) Lots of other sites allow you to upload traces, which they'll display in a list - but I haven't found any that will actually render them all simultaneously.


title: Looking for a tutorial on syncing QGIS with Trimble ProXT

text: Does anyone have any experience/advice to sync a TrimbleProXT gps antenna with the current "Lisboa" 1.8.0v? Additionally did you manage to collect features (points,line,polygons) in the field ok? let me know if there was any gripes or issues so I can be forwarned ? A fools guide step by step guide would be very usefull for this newbie! Best Ian


title: Experience with Geoserver's GML Extension

text: I'd like to know if Geoserver's GML extension works (well) with GeoServer 2.1-RC2 and whether it will enable uploading of GML files via REST API.


title: How to calculate most efficient route passing through all the houses in the world?

text: I'm new to GIS. I need some help in determining the best or most efficient route, using a flying sleigh, through all the houses in the world. One of my coworkers told me that this site would be the best place to ask, because I would find lots of helpful GIS experts. I will need some guidance about what software to use, where to get the data, and how to process it. Since I had some extra expenses this month, I would prefer some Open Source solutions. Thank you all! PS: I'm in a bit of a hurry, as I need this for tomorrow!


title: using geocoder of opentripplanner

text: I am developing a webapp using opentripplanner. I have followed 2 and five minute tutorials at their websites, there is also a opentripplaner- geocoder.war webapp, but i am unable to figure out that how can i use this geocoder.


title: How to link an application to grass?

text: I am creating an application (C#) and I want to connect it to GRASS, mainly run grass scripts from my application. This might seem to be a repeated question however I am thinking of the possibility that I will no longer open GRASS explicitly (double click the application and the GUI will pop-out). Also I know that GRASS can create PNG images but can grass do this too without being ran? I really need GRASS to be hidden from the view so that it would be neat. BTW I know that ofcourse GRASS would still be needed to be installed in the PC and not completely remove it. Are there any available tutorials other than the one on the wiki page? Or a sample code?


title: Keeping SDE for data maintenance, but use postgis for display

text: We build and maintain GIS data in a production environment with SDE on SQL 2008. We also have a need to open up WMS and WFS that will make a complete ESRI solution really expensive. With an ArcGIS Server Enterprise Advanced license, is it possible to replicate the data into a non-SDE postgis database and push deltas from the production environment?


title: How to mix polygon extrusion and different altitude modes in KML?

text: I want to create boxes by simply extruding polygons but with special "floors". I've been able to find several existing solutions for the extrusion problem. The remaining problem is that the "floor" of my boxes do not fit with the altitude modes provided by the KML reference. I want to be able to create floors with specifications like "max (100m absolute, 50m above ground)". I don't see any simple solutions. I could do a grid sampling to get ground altitude (I'm already doing this for something else) and generate only absolute altitude. But I'm afraid that I will have problems in areas with too many mountains if the sampling is too sparse, or I will have a huge KML file if the sampling is too dense.


title: How to findout direction PostGIS

text: I want to find out the direction from the 2 points from postGIS database. Please see the attached image. I want to find out the direction from the mainPt to any point. But want to get the result as NW, NE, SE or SW. Simply if I get the degree of the angle it can easily be solved. Thanks![enter image description here](http://i.stack.imgur.com/ndrEs.png)


title: BIOMOD Error message: argument of length 0

text: I’m using the package BIOMOD in R to modelling habitat in rodent species. The data formatting works but I'm getting this error message when I run Modelling:               ## Conputing the models ##     myBiomodModelOut <- BIOMOD_Modeling(biomodDataL,                                 models = c('GLM','MAXENT','GAM'),                                 models.options = myBiomodOption,                                 NbRunEval = 2,                                 DataSplit = 80,                                 Yweights = NULL,                                 VarImport = 2,                                 models.eval.meth = c('TSS','ROC'),                                 SaveObj = TRUE,                                 rescal.all.models = TRUE)      > > ################################################################################ > Creating suitable Workdir... >   >   >                         ! Weights where defined to rise a 0.5 prevalence !   >   > > -=-=-=- Lemming Modeling Summary -=-=-=-=-=-=-=-=-=- > > 1 environmental variables ( X2 ) > > Number of evaluation repetitions : 3 > > Models selected : GLM MAXENT GAM Total number of model runs : 12 > > > -=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- > > -=-=-=- Run : Lemming_PA1 > > -=-=-=--=-=-=- Lemming_PA1_RUN1 Model=GLM quadratic with no interaction >   >   >         Stepwise procedure using AIC criteriaError in 1:ncol(explVar) : > argument of length 0 >   > > # The explanatory variable is an index of productivity (i.e. NDVI). This raster, with explanatory variable (2'x 2' resolution), and vector with the response variable are in the same projection. Do you know why I have this kind of message? Thanks in advance for your help, Mary


title: Disabling/overriding internal arcmap commands with vba

text: I am creating a template for ArcGis JTX. This is a mxd file, with a vba onload event handler which removes some functionality from vanilla ArcMap. I'm stumped on one problem: double clicking on a featurelayer opens a layer properties screen. I'd like to disable that functionality. An ideal solution would be to replace the command itself with an ICommand that shows a messagebox with the text "Disabled by administrator" Another (less than ideal) solution would be to remove the ICommand from the ICommandsCollection, but I haven't found such an object in the documentation. As an alternative, I'd replace the doubleclick eventhandler from the featurelayer, disable shortcuts and remove buttons from toolbars and menus. It's a lot of programming, and I have yet to learn how to replace the eventhandler. Do you have any idea if one of these solutions can be done? tl;dr: I need a function DisableCommand(UID) in an mxd with vba.


title: Troubleshooting ArcObjects performance problems

text: We have developed a relatively complex geoprocessing routine using a C# / ArcObjects (COM) / ArcGIS Server 9.3.1 / ArcSDE 9.3.1 / Oracle 11g stack. Running the routine on a complete dataset takes many, many days and often does not complete. Troubleshooting is hindered by the opaque nature of COM and the many levels of abstraction. Specifics aside, _could anyone suggest a heuristic for troubleshooting performance problems with ArcObjects?_ Hopefully the answer could be useful to others. **Edit:** Using a profiler (JetBrains dotTrace), I can see that my problem calls are: `ITable.GetRow`, `ITable.Select`, `ICursor.NextRow`, `ITable.Search`, `ISelectionSet.Search`, `ISelectionSet.get_IDs`, `IFeatureClass.Select`, and `IFeatureClass.GetFeature`.


title: What is displayed when I am using a geographic CRS for my project

text: I understand that each geographical dataset has a CRS, and that in order to overlay them, you need to either make sure that all datasets use the same CRS, or re-project some of them (or use on-the-fly projection). However, I do not really understand one thing; when I set the project CRS in QGIS, I am setting the CRS in which I wish to view the overlay. Now, if that particular CRS is a _projected_ one, then I am comfortable viewing it on my (flat) screen. However, what if the CRS is _geographic_ , such as WGS84? Do the axes then simply correspond to latitude and longitude?


title: Limiting the maximum and minimum values of an interpolation surface tool (ArcGIS, Spline)

text: I've been using the spatial analyst spline tool in ArcGIS 10.1 and was wandering whether it is possible to limit the maximum and minimum values to those of the z-attribute. Currently the spline will exceed both these values.


title: How to create a map with OSM and webatlas.de's WMS layer?

text: I spent hours trying to figure out the following problem, I hope someone can help me. You can see here that this WMS has different Projections available: http://www.webatlasde.de/arcgis/services/Maps4BW/MapServer/WMSServer?SERVICE=WMS&request=GetCapabilities (here you can find even more information and an example map http://www.webatlasde.de/ArcGIS/rest/services/Maps4BW/MapServer ) How do I have to set up my map-object or layers so that the webatlas.de and OSM (or Google) are displayed in one Map with the correct projection? Can't I use EPSG:3857 as an substitut for EPSG:900913? In fact I can only get webaltas.de working without spatial distortion with the following configuration:               Proj4js.defs["EPSG:31467"] = "+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs";          var extent1 = new OpenLayers.Bounds(3364573, 5257369, 3635262, 5536211);          function Map_main(){              map = new OpenLayers.Map("mapdiv",{                  projection: new OpenLayers.Projection('EPSG:31467'),             displayProjection: new OpenLayers.Projection('EPSG:31467'),                    controls: [new OpenLayers.Control.Navigation(),                   new OpenLayers.Control.PanZoomBar(),                   new OpenLayers.Control.ScaleLine({ bottomOutUnits:'',bottomInUnits:'',maxWidth:150}),                   new OpenLayers.Control.MousePosition(),                   new OpenLayers.Control.LayerSwitcher()],                   units: "m",             maxExtent: extent1         });              webatlasde = new OpenLayers.Layer.WMS("Webatlas", "http://www.webatlasde.de/arcgis/services/Maps4BW/MapServer/WMSServer", {             layers: 'Halbtonvariante',             format: 'image/png'         });              map.addLayer(webatlasde);          }      If I try this, the map is distorted               Proj4js.defs["EPSG:31467"] = "+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs";          var leftBottom = new OpenLayers.LonLat(3364573, 5257369).transform(new OpenLayers.Projection("EPSG:31467"), new OpenLayers.Projection("EPSG:4326") );     var rightTop = new OpenLayers.LonLat(3635262, 5536211).transform(new OpenLayers.Projection("EPSG:31467"), new OpenLayers.Projection("EPSG:4326") );     var extent2 = new OpenLayers.Bounds(leftBottom.lon, leftBottom.lat, rightTop.lon, rightTop.lat);          function Map_main2(){              map = new OpenLayers.Map("mapdiv",{                  projection: new OpenLayers.Projection('EPSG:4326'),             displayProjection: new OpenLayers.Projection('EPSG:4326'),                    controls: [new OpenLayers.Control.Navigation(),                   new OpenLayers.Control.PanZoomBar(),                   new OpenLayers.Control.ScaleLine({ bottomOutUnits:'',bottomInUnits:'',maxWidth:150}),                   new OpenLayers.Control.MousePosition(),                   new OpenLayers.Control.LayerSwitcher()],                   units: "m",             maxExtent: extent2         });              webatlasde = new OpenLayers.Layer.WMS("Webatlas", "http://www.webatlasde.de/arcgis/services/Maps4BW/MapServer/WMSServer", {             layers: 'Halbtonvariante',             projection: new OpenLayers.Projection('EPSG:31467'),             displayProjection: new OpenLayers.Projection('EPSG:4326'),              format: 'image/png'         });              map.addLayer(webatlasde);         }      Can someone tell me how to correctly change the projection for the webatlas.de WMS? Or do I even have to use its REST-Interface with OpenLayers (I tried but miserably failed)? I hope I could make explain myself. Have a nice day! All the best Gregor


title: Network analysis service areas vs. network analysis

text: i have calculated the distance of certain houses from public parks using both Network analysis service areas and network analysis. However each method provides a different distance. Would this be expected? or should they be the same as they are both calculated along the road/path networks. The Network analysis has given shortest route, could it be that the service area wouldnt have included that short route?!! Thanks


title: How to implement Spatial Autocorrelation using QGIS or PostgreSQL (or any free application)?

text: I am a total newbie on the GIS scene so I hope this is not a stupid question. I have two layers in a PostgreSQL database. The first layer contains points where events happened in Europe. The second contains a shape file of Europe. I have been able to map the points on the map of Europe using QGIS. I would now like to implement Spatial Autocorrealation (Moran i). I have seen a few demos of this being done using **GeoDA** but it appears that it only handles a single shapefile. Can anyone point me in the right direction on how to implement this using QGIS or PostgreSQL or any other free application Thank you for your time


title: ArcObjects identify tool

text: I'm using ArcGIS 10.1 and developing an application using visual studio 2010 and ARCObject SDK 10.1. I would like to use identity tool in my application and if not possible I want to execute a query on the attribute table of the map layer when a user clicks on it. Any ideas on how I can accomplish this task? From my research,I think IQuery interface may be used; however, I would like to do the query based on a mouse click event.


title: Is there an arcpy equivalent to IPixelBlock in ArcObjects?

text: Trying to figure out a way of doing this via arcpy, but don't see anything remotely related to this. My raster dataset is a GRID.                   'pPixBlock is an IPixelBlock              'sets raster to be one column and one row in size         pSizePoint.SetCoords(1, 1)              pPixBlock = praster.CreatePixelBlock(pSizePoint)      Any insights will be greatly appreciated. Thanks!


title: Update database row with points that fall within polygons

text: All, I have a PostGIS/Postgresql DB that has two tables in it. One with point geometries and the other that represents country boundaries as polygons. I would like to add the country name that each point intersects with to each row in my point table. Maybe as one big update query. I think this is possible to do using straight SQL but I don't know where to begin. Any advice on this would be greatly appreciated... Adam


title: adding multiple points to a linestring in PostGIS

text: Though a newbie to PostGIS, I'm having trouble solving something that I think should be farily straightforward. I have two linestring tables, both consisting of simple, 2-node line objects: one has the segments of exploded road centre-lines ('r_lines'), while the other has lines ('a_lines') from an addresspoint to the nearest point on the nearest centre-line segment, which in most cases is neither the start nor end node of that roadline segment. Ultimately I need to have a table consisting entirely of (connected) 2-pt linestrings. I haven't been able to figure out how to introduce the one endpoint of the a_lines into the r_lines to create new segments, and the problem seems to be that one r_line segment might be intersected by several a_lines. So a previously described method for adding vertices to linestrings by breaking and re-merging can't be directly applied, unless I could do it iteratively (which I don't know how to do). I have also tried dropping all the points into one table and then creating a new linestring with the grouping method described here. This did create lines, but not the ones I was expecting, as they frequently seemed to double-back on each oter. Grateful for any suggestions. (And apologies for the purely textual description: as a new user of little reputation I'm prevented from uploading a helpful illustration.)


title: Is it possible to do a query to multiple tables?

text: I have a point a line and a polygon layer. They have the same attribute table. Can I do one query that refers to all the attribute tables? So I wouldn't have to do every query three times. Thanks. I'm working ith qgis, v. 2.0.1. Edit: I joined the layers but how am I doing a query that covers all the attributes? For example I'd like to ask: Which of the attributes have the name "tree"? As I have some trees along a river, I digitized them as a line, but as I have also single trees, I also have "tree"-attributes as points. How can I do just one query for this? Any suggestions?


title: How can I insert images in infowindows in #CartoDB?

text: I followed all the instrux to insert images into a CartoDB infowindow, but it hasn't worked? Can someone help me with how I can make it work? Thanks.


title: ArcMap: Append is really slow

text: I'm writing a Python script for ArcGIS 9.3.1 which is supposed to do the following:   * Create a new feature class in local FGDB   * Loop over a set of IDs (~40,000)   * Select features from an Oracle GDB with that ID   * Append features into the local feature class This was taking about 10 seconds per row which was longer than ideal, so I copied the source data to my hard drive (~250 MB) which somehow made it _even slower_. I added an index to the ID field but that didn't help either. I've narrowed the slowness down to the Append step but I cant think of how to speed it up. Does anyone have any ideas? These are the two key lines of code:               gp.SelectLayerByAttribute_management(landbase_lyr, 'NEW_SELECTION', landbase_where)     gp.Append_management(landbase_lyr, cases_lyr, 'NO_TEST')      Please let me know if I can provide any more information.


title: ArcGIS Server 10.1 - Inconsistent Querying Errors

text: I have recently upgraded to ArcGIS Server 10.1 SP1. Since making the switch, I am getting a lot of inconsistent errors, specifically performing queries against a map service's REST API with the Javascript API. I say inconsistent because an identical query will sometimes work, and other times I get an 'unable to complete operation' error. The queries are extremely simple, for example the where clause is "COUNTRY_ID = 2", and I am positive its not the where clause because it always worked coming from ArcGIS Server 10.0. Testing the query in the REST page directly from the services directory yields the same result. My result alternates between success and error - query 1: success, query 2: error, query 3: success, query 4: error. Has anyone experienced something similar? UPDATE: From info I found in this post, I found that making POST requests resolves my issues.


title: How to generate landmark-based route description using ArcPy?

text: I am now using ArcGIS 10.1 and currently a beginner in ArcPy. There is a project requiring to generate landmark-based route description. I am wondering how to use a buffer to find the landmarks in the neighborhood of the pre-defined route, in additional to provide a relationship between landmark and the route are there any alternatives way to find out this "landmark-based route description" instead of using ArcGIS existing option? The final product of the project requires a description of route in every turn/edge, such as: "Turn right in McDonald to the Manhattan Road" "Walk along the Manhattan Road, and there is a Park on your right" It is not necessary to provide a distance in the description, as the target of the description is pedestrian and he/she doesn't really need this piece of information.


title: ArcMap: Joining Tables for the Labeling Manager

text: I have two tables I would like to join, in an mdb, and have the label manager display columns from the joined table. I am able to do this through the ArcMap Gui (right clicking a layer and making a join through "Joins and Relates"), but I would like to do this programmatically. Is this possible? I joined the two tables by following this example: http://resources.arcgis.com/en/help/arcobjects- net/conceptualhelp/index.html#//0001000002zw000000 but the label manager does not seem to recognize this join and as a result I am missing columns. Thanks Edit: I am working in C#.


title: Is ArcFM or ArcMap leaking memory when switching between page templates?

text: Can anyone else reproduce the memory leak described below? Any insight into what may be causing the leak is appreciated. I have custom code that swaps between landscape and portrait templates and exports a large number of maps to PDF. Unfortunately, there is a memory leak occurring when I swap from template to template that eventually causes the ArcMap to crash. I can reproduce the memory leak using the Minerville data and workflow described below. **Configuration:** ArcFM Solution Desktop 10.1 ArcGIS 10.1 for Desktop 10.1 Windows 7 64-bit or Windows XP 32-bit or Windows 2008 64-bit **Workflow:** 1) Open Minerville database in ArcMap 2) Open “Stored Items” dropdown 3) Select “Page Templates” tab 4) Double-click the “Adjacent Grids” template to load it 5) Repeat steps 2-3 and double-click “Designer Detail Map Sheet” 6) Repeat steps 2 – 5 and watch memory usage climb in task manager Using UMDH, the largest ArcFM-related offenders seem to be these: Stack #1:               ntdll!RtlAllocateHeap+00000274     KERNELBASE!LocalAlloc+0000005F     GDI32!pmfAllocMF+0000018C     GDI32!SetEnhMetaFileBits+00000033     OLEAUT32!_PictLoadUnknownMetaFile+00000088     OLEAUT32!_PictLoadNewImage+000001A0     OLEAUT32!_PictLoadPicture+00000025     OLEAUT32!CPicture::Load+00000027     SystemCore!ObjectStream::LoadObject+000003E9     AfCore!LoadObject+0000005E     RasterGraphicElements!BasePictureElement::Load+000000ED     SystemCore!ObjectStream::LoadObject+000003E9     AfCore!LoadObject+0000005E     AfCore!ExtraSafeReadObject+00000072     PageLayout!PageLayout::Load+00000174     SystemCore!ObjectStream::LoadObject+000003E9     MMStoredStuff!DllGetClassesForCategoryReg+00011B32     MMStoredStuff!???+00000000 : 2DB04186     MMStoredStuff!DllGetClassesForCategoryReg+0000D9E5     MMStoredStuff!DllGetClassesForCategoryReg+0000C297     mmPageTemplatesUI!DllCanUnloadNow+00020402     mmPageTemplatesUI!DllCanUnloadNow+00038EAF     mmFrameworkListObjects!DllGetClassesAndCategories+000C38DC     mscorwks!CLRToCOMWorker+0000019A      Stack #2:               ntdll!RtlAllocateHeap+00000274     MSVCR90!malloc+00000079     MSVCR90!operator new+0000001F     DisplayCore!ClassFactory<CartographicLineSymbol>::CreateInstance+000000D8     AfCore!FactoryCache::CreateInstance+00000095     AfCore!FactoryCacheDispatcher::CreateInstance+00000050     AfCore!AfCreateInstance+00000040     SystemCore!ObjectStream::LoadObject+00000297     AfCore!LoadObject+0000005E     DisplayCore!MultiLayerLineSymbol::Load+000000F8     SystemCore!ObjectStream::LoadObject+000003E9     AfCore!LoadObject+0000005E     GraphicElements!LineElement::Load+00000076     SystemCore!ObjectStream::LoadObject+000003E9     AfCore!LoadObject+0000005E     GraphicElements!GroupElement::Load+00000082     SystemCore!ObjectStream::LoadObject+000003E9     AfCore!LoadObject+0000005E     AfCore!ExtraSafeReadObject+00000072     PageLayout!PageLayout::Load+00000174     SystemCore!ObjectStream::LoadObject+000003E9     MMStoredStuff!DllGetClassesForCategoryReg+00011B32     MMStoredStuff!???+00000000 : 2DB04186     MMStoredStuff!DllGetClassesForCategoryReg+0000D9E5     MMStoredStuff!DllGetClassesForCategoryReg+0000C297     mmPageTemplatesUI!DllCanUnloadNow+00020402     mmPageTemplatesUI!DllCanUnloadNow+00038EAF     mmFrameworkListObjects!DllGetClassesAndCategories+000C38DC     mscorwks!CLRToCOMWorker+0000019A      UPDATE #1: Schneider Electric's support was able to reproduce the problem and has created a bug. I was told that if a fix was made at all, it would not be available until the 10.2 release. I'll post more as I get more information.


title: Cluster GPS Points by set Distance

text: I have X,Y data that I am trying to buffer to find clusters. My definition for a cluster is " 2 or more points with a maximum distance of 200 meters". My problem is: Point A, Point B, and Point C included into a a single cluster. Point A to Point B is 50m, Point A to C is 190m, but point B to C is 225 meters. Therefore I do not believe that all points in my cluster are meeting my definition of a cluster?


title: How to highlight track in KML on hover with OpenLayers

text: I have KML layer with multiple tracks and I would like to highlight them individually on mouse hover. I tried to implement this like in the example found in OpenLayers examples, but couldn't get it work. In the example there is following code that adds polygon as feature. In my case, how do I add KML-tracks as feature so that I can highlight them with SelectFeature?               var feature = new OpenLayers.Feature.Vector(         OpenLayers.Geometry.fromWKT(             "POLYGON((28.828125 0.3515625, 132.1875 -13.0078125, -1.40625 -59.4140625, 28.828125 0.3515625))"         )     );     vectors.addFeatures([feature]);      You can find simplified version of my index.html and KML file here: https://gist.github.com/1379118


title: Reverting QGIS to version 2.2 from 2.4

text: I have QGIS 2.4 installed on a machine running Windows 7 64-bit "by accident". I was using 2.2 and when I added some package(s) it also installed the 2.4 version. Can anyone provide some advice with the following?   1. Remove version 2.4 and ensure there are no hooks, etc remaining on my computer (e.g. best method such as using Revo uninstaller or other)   2. Re-install version 2.2 (I'm assuming it would be download it from here: http://qgis.org/downloads/ , specifically the file QGIS-OSGeo4W-2.2.0-1-Setup-x86_64.exe, dated February 26, 2014)   3. Ensure that adding packages in the future doesn't automatically upgrade my version 2.2 to 2.X.


title: How do you clip a gnomonic or stereographic azimuthal projection?

text: I have polygons comprised of lat/lon coordinates and I want to project them onto a plane that's appropriate for polygon operations (such as checking for intersections). In some cases, the polygons that I want to project might fall outside the valid area that can be projected (such as gnomonic being confined to less than a hemisphere). Given an arbitrary center lon and lat, how can I clip polygons that are partially outside the valid range so I can project them? In the attached picture, note that the polygons and graticules are clipped. The way I do a gnomonic projection currently is directly project coordinates in cartesian (ECEF) out from the center of the Earth onto a tangent plane. But I obviously can't project any points or lines where ray to a coordinate is greater than or equal to 90 degrees from the center lon and lat. I'm looking for a method description as opposed to a pre-existing software package. ![enter image description here](http://i.stack.imgur.com/t05hs.png)


title: How to fix a raster with wrong Eastings?

text: I get an error when trying to warp a tif file in QGIS: ERROR 1: Too many points (35721 out of 35721) failed to transform, unable to compute output bounds. I am using QGIS 1.7.0 Wroclaw with gdal 1.8. I am trying to reproject from ETRS89 to HD72 EOV. There was an error in the data, with the UTM zone added to the eastings. That is, the x coordinates read 34327213.050 when they should read 327213.050 in UTM zone 34N. I have tried changing the .tfw file but this did not help. Then I translated the tif to a .vrt and edited the vrt, removing the rogue 34, and fixing the extent. This also did not help - I still get the same error when I try to warp. In layer properties, this raster now has the same origin, extent and SRS as other similar rasters, but only this one gives the error. How do I work this out? Any suggestions will be greatly appreciated. Cheers, Rudi ps. A similar question was asked 4 months ago, but it was not the exact same problem, and the OP never explained how the problem was resolved. Edit: as suggested, here is a link to the related question


title: OpenGeo Explorer error publishing layers to geoserver

text: I want to use Opengeo Explorer to publish layers to geoserver but i am getting an error report. I create workspaces with no problems and i publish styles with no problems. Can you pls help? > Error while storing uploaded file::Read timed out Traceback (most recent > call last): File > > "C:/Users/tomis/.qgis2/python/plugins\opengeo\gui\explorer.py", line 118, in > run command(*params) File > "C:/Users/tomis/.qgis2/python/plugins\opengeo\gui\confirm.py", line 18, in > publishLayer catalog.publishLayer(layer, workspace, overwrite, name) File > "C:/Users/tomis/.qgis2/python/plugins\opengeo\qgis\catalog.py", line 379, in > publishLayer self.upload(layer, workspace, overwrite, title) File > "C:/Users/tomis/.qgis2/python/plugins\opengeo\qgis\catalog.py", line 247, in > upload self._uploadRest(layer, workspace, overwrite, name) File > "C:/Users/tomis/.qgis2/python/plugins\opengeo\qgis\catalog.py", line 197, in > _uploadRest overwrite=overwrite) File > "C:\Users\tomis.qgis2\python\plugins\opengeo\ext- > libs\gsconfig-0.6.6-py2.7.egg\geoserver\catalog.py", line 383, in > create_featurestore raise UploadError(response) UploadError: Could not save > the layer Naxos_Points, there was an upload error: Error while storing > uploaded file::Read timed out


title: Best practices for sharing feature classes between different geodatabases

text: We're starting to work on different projects. All of them have their own data and classes but some of them are common (ex. parcels, traffic, road definition, etc.). What are the best practices for sharing across projects? I found only these approaches: I could take from the source geodatabase but I cannot create relationships or make a copy in each geodatabase but I have the risk that some of them will not be updated (in any case I have to do manually). My environment is based on ArcGIS and AutoCAD and it could be that we do an upgrade buying ArcGIS server.


title: What file type is used in the FindTask [Local] class?

text: the FindTask Class: http://resources.arcgis.com/en/help/runtime- java/concepts/index.html#/Finding_text_in_the_map/01qv00000079000000/ I tried lpk/lyr/mpk. maybe I am doing something wrong? But I am trying to use that class Locally, and it just doesn't work and there are no errors.


title: How to make variables output to user defined workspace

text: I am running a script in ArcMap and I want to reference the user defined workspace as an output instead of hard coding it in. I have the scratch output set as a user defined parameter. However, I cannot figure out how to reference that when declaring the variable. Below is a section of my code:               import arcpy     from arcpy.sa import*     import sys     from arcpy.mapping import*          #Extension     arcpy.CheckOutExtension("3D")     arcpy.CheckOutExtension("spatial")     arcpy.OverwriteOutput = True          #Argument     arcpy.env.scratchWorkspace = arcpy.GetParameterAsText(0)          #Variable     Survey_Points = "E:\\URSP\\Final\\Output.gdb\\Survey_Points"          #Process     arcpy.MakeXYEventLayer_management(Survey_Data, Survey_X, Survey_Y, "Survey_Layer", spRef, Survey_Z)          arcpy.CopyFeatures_management("Survey_Layer", Survey_Points, "", "0", "0", "0")


title: Location names autocomplete service with custom Nominatim installation

text: I'm looking to host an autocomplete service. I have installed Nominatim and imported the planet file with Nominatim. It took about 3 days to dump data into database and an additional 7 days to run indexes on an Amazon EC2 extra large instance without any database tuning. The problem with Nominatim, is that it is designed to parse a complete address and not isolated letters. I would like to modify the code to accept what we could call partial queries, i. e. if the map is showing France, and the user types the letter 'p', the first result is "Paris, France". https://github.com/twain47/Nominatim/blob/master/lib/Geocode.php#L343 I'm looking for advice from someone who has attempted this before in order to avoid diving into 1500 lines of code and hitting a snag. Geocoding on my installation is not performing well even without type ahead, and I don't know if it will be possible to host a decent service. I found this library https://github.com/komoot/photon which seems to work fast enough, but it doesn't seem to include a bounding box as a parameter to the search. I would like this service to work as good as possible with existing Open Source tools and without a third party hosting service. Is typeahead achievable using available libraries and data, with a small effort, let's say 80 hours of coding/investigation?


title: How to get EPSG with OGRinfo?

text: Is it possible to get the EPSG code of a shapefile with `ogrinfo`? I know I can get layers SRS WKT with `ogrinfo -so -al stands1.shp`, but I am looking for the EPSG code.


title: Is there a way to run a Def Query on an AOI?

text: I want to spatially only display features within an AOI boundary. Is there a way to do this like with the def query tab in a feature classes attributes?


title: ArcGIS JS - Feature Layer selection symbol before layers loaded

text: I have a feature layer that is dynamic and could pull a `geometryType: esriGeometryPoint` or `geometryType: esriGeomertyPolygon` I was wondering how would I know when to set the selection symbol to `new esri.symbol.SimpleMarkerSymbol(...)` or `new esri.symbol.SimpleFillSymbol(...)` Before I was able to just set the selection symbol to `new esri.symbol.SimpleMarkerSymbol(...)` and it would fill the polygon symbol. But now something seems to have changed where it now just places a circle or marker symbol in its place. I've tried to set the selection symbol on `featureLayer.on("load", function (evt) { ... })`, but it just seems to use the default color of the layer. Thanks in advance.


title: Calculating bounding box coordinates based on center and radius

text: I need to calculate a bounding box for osm query. I get coordinates of center of the area (WGS84) and radius (meters). How can I do it? I though about using some kind of projection, calculating projected box coordinates and then convert it back. Is it good idea? If so, which projection will suit the best to this task? I've read that miller's projection is a poor idea and it's better to use state plane projection. But on the other hand in my case bounding box doesn't need to be accurate. The most important for me is to calculate those coordinates quickly. I would be grateful if you could show me some code example that do conversions essential for calculating bounding box. I wanted to make this works on whole planet, but if there are better approximations I can assume that it should work in eastern Europe Latitude: 48.342 - 55.279 Longitude: 13.645 - 25.071


title: Best Practice for maintaining a "Cleanup Floor"

text: I do not have a GIS background and was hoping that someone could point me in the right direction as to the best procedure / software to do the following. Starting with a current LiDAR survey, I need to be able to compare an older survey (contours in a DWG) and create a composite surface that is the lower of the two at any point. These are both too large, and too close in elevation to deal with the triangulated surfaces in AutoCAD (which is my default right now). Also, moving forward I will have to continue to integrate new survey into this surface to maintain a "cleanup floor", which AutoCad does poorly, as it does not easily allow for just integrating the cut portion of any new survey. My gut tells me there must be a simpler method to do this by converting them both into a uniform grid (DEM is the convenient output from AutoCAD), and then comparing them in a GIS software but I have no idea which software to go to, or what the simplest methodology would be. Any help would be great! Thanks.


title: Lagrange Multiplier in lm.LMtests (spdep package for R)

text: Are the Lagrange values comparable in lm.LMtests? If both RLMerr and RLMlag are significant, would the test with the largest Lagrange value be considered 'better' at explaining the spatial relationship for the given model? Or 'better' at explaining which type of spatial autocorrelation is most prevalent in the model? This is my first question, so apologies for probably not following the necessary protocol/formatting/salutation OR if this belongs in some other silo of StackExchange.


title: How do I create a set of equally spaced points along a line using QChainage and export these as latitude and longitude coordinates in degrees?

text: In QGIS 2.2.0, how do I create a set of equally spaced points along a line using QChainage and export these as latitude and longitude coordinates in degrees? I am using the QChainage function in QGIS, but it is generating the following error:               An error has occured while executing Python code:          Traceback (most recent call last):     File "/Users/shafiquejamal/.qgis2/python/plugins/qchainage/qchainagedialog.py", line 126, in accept     selectedOnly)     File "/Users/shafiquejamal/.qgis2/python/plugins/qchainage/chainagetool.py", line 92, in points_along_line     features = create_points_at(startpoint, endpoint, distance, geom)     File "/Users/shafiquejamal/.qgis2/python/plugins/qchainage/chainagetool.py", line 43, in create_points_at     feature.setGeometry(point)     TypeError: QgsFeature.setGeometry(QgsGeometry): argument 1 has unexpected type 'NoneType'      I have consulted the following posts, but was unable to solve this issue: Creating points along lines? How can I represent distance using the dashes or markers of a vector line?   * This post directs me to QChainage, which I am trying to use. http://nathanw.net/2012/08/05/generating-chainage-distance-nodes-in-qgis/   * I don't know how to implement this solution; I tried to run the python code but it just produced a blank Canvas. It may be because I failed to select the lines, which I thought is done by making the layer with the line active.  How to create equidistant points in QGIS?   * Again, I'm not sure how to implement the code in the above post.  How to create points in a specified distance along the line in QGIS？   * This recommends Densify, which doesn't seem relevant to what I want to do. SEXTANTE is deprecated, and not available as a plugin in QGIS 2.2.0 (I think its functions are built in to QGIS 2.2.0). http://wiki.tuflow.com/index.php?title=QGIS_Convert_Lines_to_Points   * The MMQGIS option can do something similar, but it only gets the nodes; it does not produce equally spaced points.  I know this question has been asked and answered in other posts (most or all of which I think are listed above). The difference with my question is that I am getting this strange error when I try to use QChainage. I was able to use QChainage before and it worked, but now I am using it on layers that I generated by adding a buffer to a line, converting the buffer polygon to a line, and then removing parts of the line (see here: How do I plot a line on a map, or obtain the latitude and longitude coordinates of such a line, that is 120km north of the US-Canada border?). If you look at the line in the image below, you can see that the line is broken (it is several lines). Maybe this is the problem? ![enter image description here](http://i.stack.imgur.com/BROqX.png) I would appreciate any help. Thanks in advance! Edit Tues 6 May 2014: I am pasting below a picture of what the line looked like after I added a buffer and converted the polygons to lines, but before deleting parts of the line: ![enter image description here](http://i.stack.imgur.com/8akeP.png) Here is a zoomed in view of the line, showing that it comprises line segments: ![enter image description here](http://i.stack.imgur.com/m8tC3.png)


title: Direction on OS Meridian data

text: I've been looking to use the Ordanance Survey's Meridian data (part of the their OpenData suite) for doing vehicle routing. Whilst I can infer road class from the NUMBER property on each polyline, I cannot see anything that indicates if a polyline is unidirectional, or bidrectional. Am I just being unable to see this, or is this encoded somewhere else?


title: Search Intersections with ArcGIS Server 10 without using Address Locator

text: Is there a way to search for intersections on a centerline layer in ArcGIS Server 10 without using an Address Locator? I can send direct queries to the centerline layer for individual road segments. For example, if I add query?text=Main to the layer URL in a browser window, it will return all roads where the display field matches "Main". Is there a similar query I can do to find intersections? I know this is possible with an Address Locator, but I am trying to find an alternative. Thanks, Seth


title: How to resample a batch of rasters using OGR/GDAL?

text: I'm looking for a function similar to GRASS r.resample. I have a folder full of TIF files and want to reduce their resolution ("scale them down") to reduce file size.


title: How to convert OSM Data to PostgreSQL for Routing Purposes?

text: I'm creating a routing application using OpenStreetMap. Currently, my task is to convert the data from an OSM file to PostgreSQL. I read about osm2pgsql, but the problem is it says there that osm2pgsql is not suited for routing. I read about Osmosis too, but there is something about rails port, which I think means I need to learn Ruby in Rails. I tried studying RoR in the past and I found it too difficult for me to learn. I'm using PHP for my application. How can I able to extract the data from the OSM file and export it to my PostgreSQL. I know I might be missing something here, but I'm a newbie when it comes to mapping application. I hope someone can help me. Thanks! :D


title: ArcGIS Server directory use to store intermediate files

text: I have a tool which does the following work:   1. Connects and fetches data from hdf using impala   2. Create intermediate csv file.   3. from the csv file it creates custom feature classes. Currently when running it as a standalone tool I have used absolute paths in the local disks and it works fine. I have used random number to generate different work directories to store the data. Now I need to publish this tool as a geoprocessing service. So I need to access the server directories to save my intermediate files. As I read, each time a tool is run it creates a folder with a unique id within arcgisjob folder. If I understood right, I need access to that directory. Please suggest how to use the folder. I am using ArcGIS 10.1 Server for Windows.


title: Solve Mysql Spatial error

text: Converting shp2mysql will leads an null value in geom filed. how to rectify the the problem.               shp2mysql -s 32643 landuse landuse LIS >landuse.sql


title: Not able to display fill symbol with other styles except esriSFSSolid

text: I am using a Unique Value Renderer and the only time I see the Simple Fill Symbol display the fill is when I use esriSFSSolid style. Other styles don't display any fill when used with rest api, just the border (outline). Following is the documentation: http://resources.arcgis.com/en/help/arcgis-rest- api/index.html#//02r3000000n5000000 Simple Fill Symbol Simple fill symbols can be used to symbolize polygon geometries. The type property for simple fill symbols is esriSFS. JSON Syntax { "type" : "esriSFS", "style" : "< esriSFSBackwardDiagonal | esriSFSCross | esriSFSDiagonalCross | esriSFSForwardDiagonal | esriSFSHorizontal | esriSFSNull | esriSFSSolid | esriSFSVertical >", "color" : , "outline" : //if outline has been specified }


title: Using ArcGIS Map 10.1 how does one calculate the amount of landcover lost if sea level was to rise to a 1m and 2m Contour Line

text: I'm relatively new to ArcGIS and stuck trying to understand how to calculate the amount of each land cover type lost if sea level was to rise to a 1m and 2m contour threshold. Essentially I have created a basic map of land cover for an area adjacent to a harbor. I have extracted 0m 1m and 2m Contour lines from a Lidar-derived EDM (Using a =<2m query) and plotted them on the land cover map as seen below.![Land cover and Contour Photo](http://i.stack.imgur.com/4TKzP.jpg) I then tried using the intersect tool with both the land cover shp and contours shp but a very similar map was produced outlining essentially the same image as before, as seen below, the attributes table has included some more information but i'm not sure what to use in calcualting the amount of area lost if sea level was to rise. ![enter image description here](http://i.stack.imgur.com/TlJS9.jpg) Any information about what to do would be much appreciated.


title: Google earth flight path

text: I have a polyline defined on google earth. I want to create a flight path along the polyline so that I can create a video of how the topography changes along the polyline. Does anybody know how to do this? Any help is appreciated.


title: threading attribute table in ArcMAP

text: Is there a way to put ArcMAP's attribute table on a different thread or processor?


title: Using rpy2 in ArcGIS

text: I have a 3rd party python geoprocessing script that requires `rpy2`, and thus also r. What is the best way to get these working together? Using R with ArcGIS Desktop says the most useful option for using R with ArcGIS is via the Geospatial Modelling Environment but as near as I can tell GME doesn't use rpy2. The script I'd like to use is not that large, but I didn't write it and don't want to port it to a different module. I don't understand the workings well enough to be confident I could do so without making a mistake.