title: QGIS and Ubuntu 'browse' dialogue causing crashing

text: Below is an image of 'the suspect dialogue' in Ubuntu that when you encounter it in QGIS you are most definitely going to crash the application. http://ow.ly/i/ZGNc/original You find this particular dialogue - which doesn't appear with every menu - when you use: \- polygon centroids \- raster extraction / projection \- create new shapefile (from desktop when you browse to location of new file) ubuntu 12.04 qgis 1.8


title: How can I split a segment into two segments?

text: The problem is how to split segments or in other words how to insert new points into a polygon so that the segment where I created the new point is now two segments connected by said point. Every time I create a new point within a segment, the segment stays the same length as before and is not being split by the new point. (I don't know why the following snapshots have double columns on the left side but just ignore that. was not like this in the program). First I have the segment I want to split: ![segment is dark red between two points](http://i.stack.imgur.com/sDfi6.jpg). Then I go to edit mode for the segment layer, click on the knot tool and select the segment: ![location where the segment should be split](http://i.stack.imgur.com/7ictB.jpg). I double-click on the segment where i want a new point. The "new place for the point" appears: ![enter image description here](http://i.stack.imgur.com/u1idA.jpg). I then go to pointlayer, activate the edit mode, click on "add new element" and create a new point near the segment: ![new point is shown near the "x"](http://i.stack.imgur.com/kJYZl.jpg). After I activate the snap-function and define the radius, I place the point on the "x" with the knot-tool: ![point is now in the right place](http://i.stack.imgur.com/N74vp.jpg). Now if I go on the info-pointer and click on the segment, the segment is still as one piece and is not split by the new point: ![still one segment](http://i.stack.imgur.com/5qQbW.jpg). Why?


title: Technique using GDAL to overlay / fit one PNG map to directly overlay on another?

text: I have one .PNG basemap with known North edge of Latitude = 52.48, a South edge of Latitude = 21.94, an east edge of Longitude = -61.87 and a west edge of Longitude = -129.37. And I have another PNG image overlay with a known North edge of Latitude = 50.03, a South edge of Latitude = 22.02, an east edge of Longitude = -67.24 and a west edge of Longitude = -124.99. Is there a way using GDAL or other simple scripting tools to process the second PNG image to be in an output that directly overlays (same x/y dimensions) the basemap? I have no trouble generating map tiles that overlay each other beautifully - but I seem to be having no luck figuring out the syntax of the seemingly simpler case. Help in the form of a simple example is much appreciated! To be clear: the goal is an output .png file that has the exact same x/y dimensions as the basemap, and which can directly overlay on top of it.


title: Selection and Popups on Overlapping/Coincident Vector Layer Polygons in OpenLayers

text:
I have a GeoJSON file of polygons that I am bringing into my OpenLayers session. Some of the polygons occupy the same space on the map. I am using a select control and popup control like:-               selectControl = new OpenLayers.Control.SelectFeature(geojson_layer,     {         onSelect: onFeatureSelect, onUnselect: onFeatureUnselect,         clickout: true, toggle: true,         multiple: false, hover: false,         toggleKey: "ctrlKey", // ctrl key removes from selection         multipleKey: "shiftKey" // shift key adds to selection     });     map.addControl(selectControl);     selectControl.activate();     geojson_layer.events.on(     {         "featureselected": function(e)          {             showStatus("selected feature "+e.feature.attributes.APPREF+" on  Applications");         },         "featureunselected": function(e)              {             showStatus("");             }         });                function onFeatureSelect(feature)              {                 selectedFeature = feature;                 popup = new OpenLayers.Popup.AnchoredBubble                 (                        "chicken",                      feature.geometry.getBounds().getCenterLonLat(),                     null,                     "
App Reference: " + feature.attributes.APPREF +"
Proposal:
" + feature.attributes.PROPOSAL+ '
Further Details'+"
",                     null,                      true,                      onPopupClose                 );                 feature.popup = popup;                 popup.autoSize = true;                 popup.minSize = new OpenLayers.Size(150,50);                 popup.maxSize = new OpenLayers.Size(250,150);                 //popup.panMapIfOutOfView = true;                 popup.closeOnMove = true;                 map.addPopup(popup);             }             function onFeatureUnselect(feature)              {                 map.removePopup(feature.popup);                 feature.popup.destroy();                 feature.popup = null;             }      Where more than 1 polygon occupies the same space - can anyone advise how I allow for drilling down and/or making the other features at the location available for a popup (and therefore the link out for further information). I suspect the answer lies somewhere in Z-indexing, but can't wrap my head around how I would integrate that from the example thats provided. Presumably there is a way of handling this type of scenario - if someone could enlighten me that would be very much appreciated Thanks


title: Extent Intersection when WrapAround180 enabled

text: Allowing WrapAround180 allows map extents to go beyond the true Web Mercator extent. Once this happens I am not able to use map.extent.intersects(point). My work around was to   1. Normalize extents in lat/long   2. Convert point to lat/long   3. Split extents to two parts if crossing date line.   4. Test intersection with extent(s) Is there an easier way? Here's my code: var extent = map.extent; var pointToTest = someGraphic;                    var normMin = wmUtils.xyToLngLat(extent.xmin, extent.ymin);                var normMax = wmUtils.xyToLngLat(extent.xmax, extent.ymax);                    //test if minx is greated than maxx          if (normMin[0]> normMax[0]) {                var pointMin1 = new Point([normMin[0],normMin[1]],new SR({ wkid:4326 }));           var pointMin2 = new Point([180,normMax[1]],new SR({ wkid:4326 }));           var wmMin1 = wmUtils.geographicToWebMercator(pointMin1);           var wmMin2 = wmUtils.geographicToWebMercator(pointMin2);                var pointMax1 = new Point([-180,normMin[1]],new SR({ wkid:4326 }));           var pointMax2 = new Point([normMax[0],normMax[1]],new SR({ wkid:4326 }));           var wmMax1= wmUtils.geographicToWebMercator(pointMax1);           var wmMax2= wmUtils.geographicToWebMercator(pointMax2);                     var normExtent1 = new Extent(wmMin1.x,wmMin1.y,wmMin2.x,wmMin2.y, new SR({ wkid:102100 }));           var normExtent2 = new Extent(wmMax1.x,wmMax1.y,wmMax2.x,wmMax2.y, new SR({ wkid:102100 }));              } else {               var pointMin = new Point(normMin,new SR({ wkid:4326 }));          var wmMin = wmUtils.geographicToWebMercator(pointMin);          var pointMax = new Point(normMax,new SR({ wkid:4326 }));          var wmMax = wmUtils.geographicToWebMercator(pointMax);          var normalizedExtent = new Extent(wmMin.x,wmMin.y,wmMax.x,wmMax.y, new SR({ wkid:102100 }));             }             if (normMin[0]> normMax[0]) {             if ((!normExtent1.intersects(pointToTest.geometry) && !normExtent2.intersects(pointToTest.geometry))) return;           } else {                         if (!normalizedExtent.intersects(pointToTest.geometry)) return;        }


title: QGis: Difference between Mac OSx and Windows

text: i try to import a KML-Placemark file with QGIS (--> Add Vectorlayer, UTF-8). When i import the File with my Mac (OS 10.6) i don´t have any problems viewing all the attributes stored in the KML-File (they are stored in the Tag "ExtendedData") ![enter image description here](http://i.stack.imgur.com/CIDao.png) but when i import the same KML-File in the same QGis Version on a Windows 8 OS i only see the attributes "Name" and "Description" (The points are shown correctly) ![enter image description here](http://i.stack.imgur.com/btjyo.png) this is referring to my question here: KML in QGis with additional Data i thought the problem is solved but it is only to 50 % :) here is my updated kml-file:               <?xml version="1.0" encoding="UTF-8"?><kml xmlns="http://www.opengis.net/kml/2.2">     <Document>     <Placemark id="PM1">     <TimeStamp>             <when>2007-01-14T21:05:02Z</when>           </TimeStamp>     <name>Data1</name>     <description>AA</description>     <ExtendedData>     <Data name="Orientierung">     <displayName>Orientierung</displayName>     <value>1</value></Data>     <Data name="Brennweite">     <displayName>Brennweite</displayName>     <value>279/100</value>     </Data>     <Data name="Datum">     <displayName>Datum</displayName>     <value>2013:04:04 13:20:49</value>     </Data>     <Data name="Pfad">     <displayName>Pfad</displayName>     <value>file:///C:\Temp\20130404_132050.JPG</value>     </Data>     <Data name="Kamera">     <displayName>Kamera</displayName>     <value>GT-P3110</value></Data>     </ExtendedData>     <Point>     <coordinates>15.448561,47.11721,0</coordinates>     </Point>     </Placemark>          <Placemark id="PM2">     <TimeStamp>             <when>2007-01-14T21:05:02Z</when>           </TimeStamp>     <name>Data2</name>     <description>AA</description>     <ExtendedData>     <Data name="Orientierung">     <value>5</value></Data>     <Data name="Brennweite">     <value>400/400</value>     </Data>     <Data name="Datum">     <value>2013:04:04 13:20:50</value>     </Data>     <Data name="Pfad">     <value>file:///C:\Temp\HUGO.JPG</value>     </Data>     <Data name="Kamera">     <value>GT-P3110</value></Data>     </ExtendedData>     <Point>     <coordinates>16.458561,48.15721,0</coordinates>     </Point>     </Placemark>     </Document>          </kml>      does anyone know why this problem is appearing?? hope you can help me!   thanks michik


title: Geoserver won't cache wms server with style in workspace

text: I am using Geoserver 2.5 on ubuntu and am trying to cache a layer that has it styles in workspaces. When I try to start the cache with seed/truncate nothing happens and I get things like this in the logs:               2014-04-07 10:56:03,939 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - Converted URL to lowercase, from: '/gwc/rest/seed/wis_gent:maaiveld'; to: '/gwc/rest/seed/wis_gent:maaiveld'  and httpMethod= GET     2014-04-07 10:56:03,939 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - ~~~~~~~~~~ antPath= /** methodList= [GET]     2014-04-07 10:56:03,939 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - Candidate is: '/gwc/rest/seed/wis_gent:maaiveld'; antPath is /**; matchedPath=true; matchedMethods=true     2014-04-07 10:56:03,939 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - returning ADMIN     2014-04-07 10:56:03,940 TRACE [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/gwc/rest/seed/wis_gent:maaiveld]     2014-04-07 10:56:03,941 DEBUG [org.geoserver.filters] - Compressing output for mimetype: text/html;charset=ISO-8859-1     2014-04-07 10:56:03,942 DEBUG [org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1] - SecurityContextHolder now cleared, as request processing completed     2014-04-07 10:56:04,028 TRACE [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/gwc/rest/web/geowebcache_logo.png]     2014-04-07 10:56:04,029 DEBUG [org.geoserver.filters] - Not compressing output for mimetype: image/png     2014-04-07 10:56:04,034 DEBUG [org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1] - SecurityContextHolder now cleared, as request processing completed     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - Converted URL to lowercase, from: '/gwc/rest/seed/wis_gent:maaiveld'; to: '/gwc/rest/seed/wis_gent:maaiveld'  and httpMethod= POST     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - ~~~~~~~~~~ antPath= /** methodList= [GET]     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - Candidate is: '/gwc/rest/seed/wis_gent:maaiveld'; antPath is /**; matchedPath=true; matchedMethods=false     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - ~~~~~~~~~~ antPath= /** methodList= [POST, DELETE, PUT]     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - Candidate is: '/gwc/rest/seed/wis_gent:maaiveld'; antPath is /**; matchedPath=true; matchedMethods=true     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - returning ADMIN     2014-04-07 10:56:05,611 TRACE [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/gwc/rest/seed/wis_gent:maaiveld]     2014-04-07 10:56:05,613 DEBUG [org.geoserver.gwc.layer] - --> GWC Seeder Thread-6 submitting getMap request for meta grid location [0, 0, 0] on org.geoserver.gwc.layer.GeoServerMetaTile@27897fc     2014-04-07 10:56:05,614 INFO [org.geoserver.wms] -      Request: getServiceInfo     2014-04-07 10:56:05,616 DEBUG [org.geoserver.ows] - Getting layers and styles from LAYERS and STYLES     2014-04-07 10:56:05,616 ERROR [org.geoserver.ows] -      org.geoserver.platform.ServiceException: No such style: wis_bestemming         at org.geoserver.wms.map.GetMapKvpRequestReader.parseStyles(GetMapKvpRequestReader.java:1386)         at org.geoserver.wms.map.GetMapKvpRequestReader.read(GetMapKvpRequestReader.java:318)         at org.geoserver.wms.map.GetMapKvpRequestReader.read(GetMapKvpRequestReader.java:80)         at org.geoserver.ows.Dispatcher.parseRequestKVP(Dispatcher.java:1412)         at org.geoserver.ows.Dispatcher.dispatch(Dispatcher.java:622)         at org.geoserver.ows.Dispatcher.handleRequestInternal(Dispatcher.java:263)         at org.springframework.web.servlet.mvc.AbstractController.handleRequest(AbstractController.java:153)         at org.geoserver.gwc.GWC.dispatchOwsRequest(GWC.java:1162)         at org.geoserver.gwc.layer.GeoServerTileLayer.dispatchGetMap(GeoServerTileLayer.java:572)         at org.geoserver.gwc.layer.GeoServerTileLayer.getMetatilingReponse(GeoServerTileLayer.java:517)         at org.geoserver.gwc.layer.GeoServerTileLayer.seedTile(GeoServerTileLayer.java:711)         at org.geowebcache.seed.SeedTask.doActionInternal(SeedTask.java:136)         at org.geowebcache.seed.GWCTask.doAction(GWCTask.java:76)         at org.geowebcache.seed.MTSeeder.call(MTSeeder.java:36)         at org.geowebcache.seed.MTSeeder.call(MTSeeder.java:25)         at java.util.concurrent.FutureTask.run(FutureTask.java:262)         at java.util.concurrent.ThreadPoolExecutor.runWorker(ThreadPoolExecutor.java:1145)         at java.util.concurrent.ThreadPoolExecutor$Worker.run(ThreadPoolExecutor.java:615)         at java.lang.Thread.run(Thread.java:744)     2014-04-07 10:56:05,617 ERROR [org.geowebcache.seed.MTSeeder] - Problem communicating with GeoServer     org.geowebcache.GeoWebCacheException: Problem communicating with GeoServer         at org.geoserver.gwc.layer.GeoServerTileLayer.getMetatilingReponse(GeoServerTileLayer.java:522)         at org.geoserver.gwc.layer.GeoServerTileLayer.seedTile(GeoServerTileLayer.java:711)         at org.geowebcache.seed.SeedTask.doActionInternal(SeedTask.java:136)         at org.geowebcache.seed.GWCTask.doAction(GWCTask.java:76)         at org.geowebcache.seed.MTSeeder.call(MTSeeder.java:36)         at org.geowebcache.seed.MTSeeder.call(MTSeeder.java:25)         at java.util.concurrent.FutureTask.run(FutureTask.java:262)         at java.util.concurrent.ThreadPoolExecutor.runWorker(ThreadPoolExecutor.java:1145)         at java.util.concurrent.ThreadPoolExecutor$Worker.run(ThreadPoolExecutor.java:615)         at java.lang.Thread.run(Thread.java:744)     Caused by: java.lang.IllegalStateException: Expected: RenderedImageMap, got null         at org.geoserver.gwc.layer.GeoServerTileLayer.dispatchGetMap(GeoServerTileLayer.java:575)         at org.geoserver.gwc.layer.GeoServerTileLayer.getMetatilingReponse(GeoServerTileLayer.java:517)         ... 9 more      Using a style without a workspace works great. Does anyone have any idea what could be wrong and if I can use a style in the workspace?


title: python.multiprocessing and “FATAL ERROR (INFADI) MISSING DIRECTORY”

text: I have a question which is closely related to this old question here: python.multiprocessing and "FATAL ERROR (INFADI) MISSING DIRECTORY" I have also looked at this, which seems to relate to something similar: Multiprocessing Errors - ArcGIS implementation just like in the linked question, I run into the error: "FATAL ERROR (INFADI) MISSING DIRECTORY" at random times, but never the BND-error I have multiple folders with thousands of shapefiles in each, where I want to convert two of the shapefile fields to two different rasterlayers. the shapefiles were converted from csvs via multiprocessing, with no issues whatsoever. I have tried most of the workarounds from the previous question, such as creating temporary and unique scratchworkspaces and workspaces, check that the infile exists and create a string in memory and flush it, but to no avail. I have also tried splitting the workload, so that it will first write all rasters for one field, then do all the rasters for the second field. so my question is, would anybody know of another potential workaround? I was thinking that it might be possible to somehow make sure that when a new raster is written to the folder, that no other processes try to access the space, but I'm unsure how to code it. relevant code excerpts below:               import arcpy     import os     import sys     import arcgisscripting     import multiprocessing     from functools import partial          def multiproc(infile,path,savepath1,savepath2):           try:                 arcpy.Exists(infile)                  # Set local variables             inLayer = path + "sample.xml"             fieldList = arcpy.ListFields(infile[:-4]+".dbf", "field1")             inData = infile + " field1"             outLayer = infile[:-4]                  # Check out the ArcGIS Geostatistical Analyst extension license             arcpy.CheckOutExtension("GeoStats")             # Execute CreateGeostatisticalLayer             arcpy.GACreateGeostatisticalLayer_ga(inLayer, inData, outLayer)                arcpy.GALayerToGrid_ga(outLayer, savepath1 + "/" + infile[:-4])             arcpy.Delete_management(outLayer)                  #starting second raster creation process             arcpy.Exists(infile)             fieldList = arcpy.ListFields(infile[:-4]+".dbf", "field2")             inData = infile + " field2"             outLayer = infile[:-4]                  # Check out the ArcGIS Geostatistical Analyst extension license             arcpy.CheckOutExtension("GeoStats")             # Execute CreateGeostatisticalLayer             arcpy.GACreateGeostatisticalLayer_ga(inLayer, inData, outLayer)             arcpy.GALayerToGrid_ga(outLayer, savepath2 + "/" + infile[:-4])             arcpy.Delete_management(outLayer)                   except:             import logging             LOG_FILENAME = 'c:/error.txt'             logging.basicConfig(filename=LOG_FILENAME,level=logging.DEBUG,)                  logging.debug(infile)                  logging.exception('Got exception on main handler')             pass               if __name__=='__main__':     # Set environment settings     #setting paths for the directories           ##Finding paths and list of files                #multiprocessing part         Pool=multiprocessing.Pool         workers=multiprocessing.cpu_count()         pool=Pool(processes=workers)              partial_multiproc=partial(multiproc,path=path,savepath1=savepath1,savepath2=savepath2)         pool.map(partial_multiproc,filelist)         pool.close()         pool.join()      EDIT1: after having moved the data and tried on an internal harddrive (the previous attempts have been on an external HD), I also got the BND-error on rare occassions, whereas the INFADI is much more prevalent. EDIT2: I have been unable to find a full solution so far, but a potential workaround lies in creating temp-folders that are deleted with arcpy at the end of the multiproc-function. even though the error persists, the faulty folders are - as far as I can tell - not deleted, meaning that at least it is possible to know which rasters were not completed and redo these later.


title: Coarse heat map placement using gdal

text: I have generated a coarse heatmap of where data is available by 1 degree x 1 degree cells that I want to overlay on top of a world image. Using a python script I merged two images (heat map and a world map -- generated with gdal_rasterize). Reviewing the picture I have a vertical registration issue (I am guessing a flattening issue). The picture can be sent if needed but I cannot post it at this time (I do not have a reputation >= 10). How can I use gdal to merge these images correctly? Thanks


title: Off by about 85 meters when converting NAD27 UTM to WGS84 lat/lon

text: Looking through other answers on converting to or from NAD27 CONUS, I don't see an answer that seems to fit my problem. I am creating paper topographic maps by layering information derived from USGS DEM data (NAD83) for contours and slope shading. I am using OpenStreetMap (OSM) data (WGS84) for trails, roads, water features, etc. Finally I am overlaying a UTM grid. All this is being formatted using Mapnik. Everything is aligning well (trails and road appear correct with respect to the contours. The waterways from OSM also nicely match up with the drainages as shown by contours derived from USGS DEM data, etc.) And if I use WGS84 for my UTM grid it seems to line up well too, at least the positions of map features are the same within my measuring capability as a USGS topo map. But if I overlay a NAD27 UTM grid my easting grid lines appear to be off by about 85 meters. The northing grid lines appear to be correct. I could live with 5 or 10 meters as that is the accuracy of my GPS but not 85 meters. I create the UTM grid by determining the UTM measurements for the edges of my map and then going to points the nearest even 1000 meters off paper and putting ways/lines between the points so located. The position of the points is converted to from UTM to lat,lon and the points and ways are written to a OSM XML file which is used as the source for the UTM grid layer in Mapnik. Conversion of UTM to lat,lon is using pyproj in a Python script. My instantiation of the projection objects is as follows:               wgs84utm = pyproj.Proj(proj='utm', zone=v['zone'], datum='WGS84', ellps='WGS84')     nad27utm = pyproj.Proj(proj='utm', zone=v['zone'], datum='NAD27', ellps='clrk66')      When I convert a point from UTM to lat,lon I use the following:               if (proj == 'wgs84'):         point = wgs84utm(easting, northing, inverse=True)     elif (proj == 'nad27'):         point = nad27utm(easting, northing, inverse=True)      One thing that confuses me is I only see one NAD27 UTM setup in the pyproj espg list. Based on my GPS and other mapping programs I was expecting to see several (CONUS, Alaska, Canada, etc.). There is a NAD27 CONUS but it is an Albers projection, not a UTM projection. Based on values for that I tried               nad27utm = pyproj.Proj(proj='utm', zone=v['zone'], datum='NAD27', ellps='clrk66', lat_1=29.5, lat_2=45.5, lat_0=23, lon_0=-96, x_0=0, y_0=0)      But the results were the identical. I am a newbie with python and map projection conversions, so I'm pretty sure that I've made some mistake in my usage of the projection library but don't know what. I'd appreciate any suggestions on where I have gone wrong. Thanks!


title: Building real-time point/polyline connection in GIS

text: In the first instance I need to create the network by constructing the basic elements, the nodes and links. This is relatively easy work as SharpMap has provided geometry tools that allow users to directly create point/polyline via simple mouse clicks e.g. capturing coordinate to generate point and two coordinate to construct line. These functions have been encapsulated and can be directly used. For this task, however, every time I construct a new polyline, I need to keep it connected with two nodes that I manually created, so then the network is connected otherwise the node and link will be separated. I really want to know the programming procedure on how to do this in real- time. That is, every time I create a line, I hope it will search the existing points and display if there are existing points that can be connected. If it can be connected, e.g. by displaying different colour, and then I can snap one of nodes (starting node or ending node) of the line on that point. If I do this way, I think I can construct the connected network as required. I am wondering if anyone has such experience to build such function, or you can provide me some guidance in this aspect? Thanks for the helps.


title: GeoServer map layers/tiles Loaded in the screen or Not

text: I'm using OpenLayers - 2.11,GeoServer 2.13, ExtJS 3.4, and Java for developing GIS based web applications. Please guide me to check if all my GeoServer map layers/tiles Loaded in the screen or not. While dragging map is taking time to load next tiles. So that time i want to mask all the map body and unmask after loading all the map tiles. ![enter image description here](http://i.stack.imgur.com/lOVyq.png)


title: how to remove and change qgsvertexmarker in QGIS using c++ code?

text: In QGIS following code is used to set marker at point A.                 QgsVertexMarker * marker = new QgsVertexMarker(mPlugin->iface()->mapCanvas());       marker->setColor(Qt::black);       marker->setIconType(2);       marker->setCenter(A);      How to delete this marker using C++ code? In python it is                canvas.scene().removeItem(marker)      I am using this marker in a loop where point A will change continuously so i would like to delete marker at previous position and insert in new. when i use below code marker is deleted permanently.                 mPlugin->iface()->mapCanvas()->scene()->removeItem(marker);      so how to use this marker to plot moving point?


title: How to work with the STRtree Index when object in the index change it's position?

text: I'm work with the NTS(NetTopologySuite), I get large amounts of Point geometries, for fast querying the point object, I make the STRtree Index like this, _STRtree tree = new STRtree(); foreach (var p in allPoint) { tree.Insert(p.EnvelopeInternal, p); }_ While some point objects in the "allPoint" change it's position. And I find I can't query the point which change position anymore.. How to work with it? Thanks.


title: Access Stamen tiles flavors with OpenLayers

text: i'm using OpenLayers (v2.13.1) and would like to use Stamen Toner tiles as background. The example is simple to understand : http://maps.stamen.com/test/openlayers.html But i'd like to use the "lite" variant of the toner tiles, how is it possible with OpenLayer ? Thanks


title: CalculateField_management

text: How to put the field from arcpy.GetParameterAsText() into the expression here: CalculateField_management(in_table, field, **expression** , {expression_type}, {code_block}) Now it is like **"!Name!.replace(\"о\",\"r\")"** but here the field _!Name!_ is already firmly defined and isn't depend on arcpy.GetParameterAsText() thank you


title: Can we have a local BingMap server or Google Earth server?

text: Is it possible to have all the map data come from local server instead of Google Earth or Bing Map but still use their tools as front end. For example Bing Map control (from their SDK) would be talking to some local server in the network. I want this kind of setup in a place where there is no internet and the client wants to use Bing Map toolset. Same question applies to Google Earth. OR is it that its only possible with something like ArcGIS like server software but in that case Google Earth or Bing Map control won't be its front end. Please excuse me for my lack of understanding of GIS in general


title: Field Calculator If/Then in ArcGis 10.1

text: I've been looking through these threads, i found something similar and thought it could be the answer, but it doen't seem to work for me... I could do this the slow way with Select By Attributes, but It's time consuming and I will have to repeat it several times for other works. I have a field called "Kultur" (String) with different codes (DG, A, FF, ...). I want to change this information into a numeric information, so i added a new field called "Kultur_Cl" (Short Integer). In my case if "Kultur" is DG, "Kultur_Cl" should be 1. If "Kultur" is A, "Kultur_Cl" should be 2, ... Here is what i tried : ![enter image description here](http://i.stack.imgur.com/JFKae.png) And here is what it looks like in ArcMap: ![enter image description here](http://i.stack.imgur.com/d1nDV.png) Once it's done i would like to do : If "Kultur_Cl" = 1 then "Kultur_Pkt" = 15, if "Kultur_Cl" = 2 then "Kultur_Pkt" = 18, ... Would this go the same way, or is it different because there would only be numeric data? Thanks for your help !!


title: Geometry gets modified when updating only attribute data

text: Using ArcGIS Runtime for Windows Mobile 10.1. We are trying to update attribute data of features and somehow the geometry gets also updated when we save the changes to the mobile cache and synchronize them back to the database. If I update multiple features, all previously updated features get the geometry of the feature that was last updated. Here is an example code which reproduces the issue. The firstRow gets the geometry of secondRow when the data is synchronized back to the database. In mobile cache everything seems to be correct, atleast I can see the features correctly on the map and select them. And also the attribute data is correctly synchronized.                       int firstId = 2105650708;             int secondId = 2105650709;                  MobileCache mobileCache = new DatabaseService().mobileCache;                  FeatureSource source = mobileCache.FeatureSources["TABLE"];                  FeatureDataTable firstDt = source.GetDataTable(new QueryFilter(                 String.Format("{0} = {1}", Columns.OBJECTID, firstId)));             FeatureDataRow firstRow = firstDt[0];             firstRow[Columns.DESCRIPTION] = "FIRST";             firstDt.SaveInFeatureSource();                  FeatureDataTable secondDt = source.GetDataTable(new QueryFilter(                 String.Format("{0} = {1}", Columns.OBJECTID, secondId)));             FeatureDataRow secondRow = secondDt[0];             secondRow[Columns.DESCRIPTION] = "SECOND";             secondDt.SaveInFeatureSource();      Has anyone experienced anything similar?


title: How to get turn angle, line orientation and line length in ArcMap 9.3.1?

text: I am mapping and analysing animal movements using arcgis 9.3.1 along with Hawthorne Beyer's GME, and am trying to utilise the correlated random walk commands, but in order to do this I need to find out the turn angles, step length and line orientation. Is there a way to easily do this in arcmap for an entire vector layer? Such as is done with calculating total or portion line length by adding a new field and calculating geometry? I am not adept at using scripts, so if the solution requires this, can you also include instruction on how to run it? Very many thanks in advance,


title: Creating weighted list using Select tool in ModelBuilder?

text: I've done a raster to point conversion and have a column of the grid codes. Let's say it is composed of 10 numbers [9,5,7,1,7,3,7,2,8,1]. Right now I am doing a selection on those numbers based on the top 50% of the Max value(Select values greater than 4.5, which gives {9,8,7,7,7,5}) What I would like to do is take the original array and select all of the values that make up greater than 50% of the total(Total is 68/2 = 34, so SUM(9,8,7,7,7) > 34) and have those as my selection points. I think I can sort the data and use another 1 or 2 calculate fields to get the sum of the column and the 50% value, but how would I go about adding the values to get above 50% and save them as the selection?


title: Locating shapefile of European railways for goods transportion?

text: We are trying to make a project about European goods trade so we need railway data in shapefile format. More specifically, we are trying to find the potential of mine products being exported from Turkey to Europe by using Geomarketing tools. Basically, we need a European railways layer for ArcGIS (excluding passenger- only railways).


title: How to release Events

text: I want to use the OnCreateFeature-Event to combine it with other actions. So I use               Events.OnCreateFeature += delegate(IObject obj) { LoopTools_OnCreateFeature(obj); };      to listen to this event. That works very well. But how can I stop listening? I tried               Events.OnCreateFeature -= delegate(IObject obj) { LoopTools_OnCreateFeature(obj); };      but this didn't release the event. So when I looped through some sweeps, it happened that I have wired the event two or more times and the corresponding actions are done to often. I also tried               Events.OnCreateFeature -= LoopTools_OnCreateFeature;      but that has the same effect. What am I doing wrong? I'm thankful for any suggestions.


title: Counting self intersections of linestring using PostGIS?

text: Is there a way of using PostgreSQL + PostGIS for finding the number of self intersections in a linestring? I have a linestring that goes forward but in circles, and I need to find the number of times it crosses itself.


title: How do you memorize the meaning of longitude and latitude?

text: Longitude and latitude are concepts I usually don't need an then I sometimes suddenly need to learn their meaning for a period of time b/c I am using them. Is there an aide-memoire to finally remember their meaning until the end of my days? What is the etymology of these two terms?


title: Can OGR handle Esri's Geodatabase extensible markup language (XML) files?

text: The USDA Forest Service distributes GIS data as Esri Geodatabase extensible markup language (XML) files. I downloaded one in the hopes that I would be able to convert it to something less esoteric using OGR. However:               $ ogrinfo ~/Downloads/S_USA.National_Trails.xml     FAILURE:     Unable to open datasource `/Users/martijnv/Downloads/S_USA.National_Trails.xml' with the following drivers.      Unfortunately, neither the FileGDB driver nor the OpenFileGDB driver's pages explicitly mention this XML format. It seems to be defined here on ESRI's website. Any way to read / convert this with FOSS4G tools?


title: Cannot see the shapes in arccatalog or arcmap after creating spatial indexes

text: I created a spatial index on a feature class, but when i try to load it in arcmap, i cannot see the shapes. When i remove the spatial index, i can see the shapes in arcmap/arccatalog.


title: Can you use wildcards for column headers in Field Calculator?

text: I know you can use wildcards for fields like so (in the Field Calculator):               case when "column" ILIKE '%example%' then 1     else 0     end      Can you do something similar for the column headers themselves? I ask this because I'm joining several layers (using the **Joins** properties) and it would make it easier for me to edit one set of filters instead of editing several filters since the column names will have to include the joined layer name. I'm using QGIS 2.2.


title: BBOX and Projections of WMS services in Geoserver

text: we have tried to overlay a raster layer on Google map and have some difficulties to align the raster images onto the base map. When it is on the overview (Zoom level 1), it looks like this: ![Overview of the raster](http://i.stack.imgur.com/d7qjr.png) On a lower resolution, you could still spot the mis-alignment in high altitude area such as North Europe and Canada: ![Zoom in](http://i.stack.imgur.com/O4m6G.png) BUT when we zoom into the highest level, the problem has gone: ![enter image description here](http://i.stack.imgur.com/UOUrq.jpg) We have tried these approaches to address this problem:   1. Change the projection of the raster layer on Geoserver from WGS 84 to Web Mercator, does not even shown the image on Google Map.   2. We try to change the bbox option in the request but it returns strange results that we cannot even zoom in. Do you have any idea why this happen? What can we do to address this problem? Thank you very much!


title: How to calculate wait time on a one-lane road with two-way travel?

text: I want to know how to calculate the wait time on a single-track/one-lane road (only room for one car at a time) for one car waiting for another car to pass by. For example, two cars drive on the one-lane road in opposition to each other. When the two cars meet in the road, one car should wait for the other car to pass. I wonder if Network Analyst (especially, Vehicle Routing Problem) can solve this problem.


title: How to zoom to an ILayer?

text: I used `IFeatureLayerDefinition.CreateSelectionLayer()` to create a new `IFeatureLayer` then i cast this one to `ILayer` and set it to the Map, my problem is that i can't zoom programmatically to this `ILayer`, i tried to call `ICommand.ControlsMapFullExtentCommandClass` but it didn't work. This is my code :               public static void myFunction(ITable table_commune, string code_commune)     {                  if(table_commune == null || code_commune == null) return;                  ILayer layer = (ILayer)table_commune;             IFeatureLayer featureLayer = (IFeatureLayer)layer;                  IFeatureLayerDefinition featureLayerDef = (IFeatureLayerDefinition)featureLayer;                  IFeatureLayer new_featureLayer = featureLayerDef.CreateSelectionLayer("Commune", false, null, "NUM_COMMUNE = " + (code_commune));             ILayer new_layer = (ILayer)new_featureLayer;                  Map.AddLayer(new_layer);             Map.ActiveView.Refresh();                  ICommand command = new ControlsMapFullExtentCommandClass();             command.OnCreate(Map.Object);             if(command.Enabled) command.OnClick();     }


title: Facebook Place Api

text: Do you know if there are a requirements or restriction on Facebook Place API? Is it possible to download the data? Is it open-source?


title: Can I change the font type in all Text_Elements in a page?

text: I want to change all the font styles in a number of mxd's from "Lucinda Sans" or whatever the current one is to "Arial" without changing the font size. I wrote a simple script to find/replace text but can't figure out how to change the font except by adding html code as in the example below (see image as well). But this doesn't change the root font and also I want it to change the font of all text strings in the Text and Legend elements. No need to search/replace any text. Thanks in advance. Screen capture and Code are below ![enter image description here](http://i.stack.imgur.com/IUBtV.gif)               import os, arcpy, arcpy.mapping     oldText = arcpy.GetParameterAsText(0)     newText = arcpy.GetParameterAsText(1)     pdf = arcpy.GetParameterAsText(2)     mxd=arcpy.mapping.MapDocument("Current")          # Loop through each text element in the map document     for textElement in arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT"):              # Check if the text element contains the out of date text         if textElement.text == oldText:              # If out of date text is found, replace it with the new text             textElement.text = '<FNT name="Arial">' + newText + '</FNT>'          #while PDF = 1:     # Export the updated map to a PDF     #path = r'P:\2011\Job_031_TownPlanning_SeriesProduction\Working\draft_output\esri\2\wetland\'     #arcpy.mapping.ExportToPDF(mxd, newText + "-TownPlanning_Wetlands_Overlay_Ver2.pdf")          #arcpy.mapping.ExportToJPEG(mxd, outFile)     #outFile = newText + df.name + ".pdf"          # Clean up the MapDocument object by deleting it     del mxd          arcpy.RefreshActiveView()     arcpy.RefreshTOC()


title: Out of Range error when creating thumbnails of MXDs

text: I am tryting to create thumbnails of all the mxd's in a directory but no mater how high I set the value of the "While" loop it gives me an "Index out of range" error... Also can I get it to look at subdirectories as well? best,               #Make thumbnails from all mxd's in current directory          import arcpy, glob, os, sys     from arcpy import env     from arcpy import mapping     env.overwriteOutput = True          #path = sys.argv[1] # for using in toolset, use directory parameter     path = os.getcwd() # Script in same directory as files being processed     mxdList = glob.glob(path + "\*.mxd")     env.workspace = path     x=0          print 'Reading mxd files from ' + path          mxd_List = glob.glob('*.mxd')     count_Mapdocs = len(mxd_List)     print mxd_List,count_Mapdocs          while x < 50:      mapDocs=mxd_List[x]      mapDoc=arcpy.mapping.MapDocument(mapDocs)      framesList=arcpy.mapping.ListDataFrames(mapDoc)      mapDoc.activeView=framesList[x]      mapDoc.makeThumbnail()      mapDoc.save()      print "Created thumbnail for " + mxd_List[x]      print "x= "+ str(x)      x=x+1      print "x updated =" + str(x)


title: Looking for a region growing algorithm

text: I work with a raster set representing land values and I would like to automatically create polygons with random shape but with certain size and/or land value given a centroid point of each polygon. The centroid point will act as the initial raster cell (a pixel called as a seed cell) that will belong in that certain polygon and then a region growing algorithm will begin adding pixels around a centroid pixel until the constraint of size and/or land value of a polygon will be satisfied. Is anybody who has any idea about an existing region growing algorithm? or at least which is the class of ArcObjects or Interfaces that I may work with, so as to read the value of each pixel, create a new polygon by adding new pixel to it? I hope that the above make sense and somebody may help me. ### Edit Basically I intend to automatically create polygons with regular shapes and a series of constraints such as size, land value etc. using a genetic algorithm (GA). But first I need to create a set of random solutions in order to feed the GA.So, for instance on land block which is an area enclosed by roads I have already saying 4 centroid points. Each centroid point represents the approximate location of each new land parcel. In addition, each centroid is associated with attributes of each parcel i.e. size and land value. So, I want to begin by each centroid as a seed point of the region growing algorithm to starting creating a random shape for each parcel based on each centroid. I hope that the above make sense. I look forward how can I create this region growing algorithm or if there is already one in VBA and ArcObjects.


title: Display Feedback "Stutters" When Moving Features

text: I have a custom tool for moving features that implements IDisplayFeedback and geometry-specific CoClasses to show the user where they are moving the feature.                   IMovePolygonFeedback polygonFeedback;         bool feedbackStarted;              public override void OnMouseDown(int Button, int Shift, int X, int Y)         {             feedbackStarted = false;                  //A BUNCH OF STUFF HERE TO SELECT A FEATURE ON MOUSEDOWN                  geometry = feature.Shape;                  polygonFeedback = new MovePolygonFeedbackClass();             polygonFeedback.Display = doc.ActiveView.ScreenDisplay;             polygonFeedback.Start(geometry as IPolygon, doc.CurrentLocation);             feedbackStarted = true;         }              public override void OnMouseMove(int Button, int Shift, int X, int Y)         {             IPoint point = doc.ActiveView.ScreenDisplay.DisplayTransformation.ToMapPoint(X, Y);             if (polygonFeedback != null && feedbackStarted)                 polygonFeedback.MoveTo(point);         }      What sometimes happens when the mousedown and mousemove happen quickly is that the feedback display stutters for a second, leaving a (what I assume to be) graphic, then continues on with a new graphic moving around with the mouse until mouseup. Looks kinda like this: http://i.imgur.com/99fj8.png At mouseup I change the feature shape to the new position and store the feature. I union the original feature extent with the new feature extent and pass that envelope to a partial refresh. Problem is sometimes parts of the "stutter" aren't in that envelope so I get dangly bits: http://i.imgur.com/v1oxx.png I'm trying to avoid straight refreshes but these dangly bits are not pretty and don't go away until a refresh happens (forced by user or move/pan, etc.) Questions:   1. Does anybody know why the "stutter" is happening and how can I prevent it?   2. Is there a way to partially refresh the dangly bits?   3. What is the feedback? Is it a graphic? Can I trap it during the onmove to add it to the partial refresh envelope?   4. Is there an infintely better way to do this that I'm mising?  Help!


title: Table Borders won't print when added to display

text: I've encountered this problem every so often over the last few years. When adding tables to the display in ArcMap, the following occurs when printing the map:   1. All table borders visible on the display but do not print   2. All Table borders visible on the display, but only top and left border prints   3. Only top and left table borders visible, yet all borders print   4. All table borders visible and all print (as it should be) What frustrates me about this is that it is extremely inconsistent - I can have 4 tables on my display and 2 will print perfectly while the other 2 won't. I can move the tables slightly, then the original 2 won't print right while the other 2 will. Today I can print it perfectly, tomorrow I can open that same mxd and print without making any changes, and the tables will print wrong. As a workaround I've always had to draw a rectangle around the table, align the borders and group each table with its border (and if its more than 3 tables, this quickly becomes annoying). This problem has "haunted" me since version 9.2.


title: Why does ESRI make me store my DBF doubles as strings?

text: As far as I can tell I have to store all numeric values in a shapefile's associated DBF file as strings (yes, the N type is essentially a string). Can anyone confirm that ESRI does not support the use of Double (O) nor Long (I) data types in the .DBF file? I'm using this DBF reference **Based on the response** thus far I apparently did not provide enough information. My question is not one of precision but rather one of optimization. When manually setting a dbf field to type O (for double) in the dbf file, I receive this message when trying to load the standalone table into ArcMap. ![Field type fail](http://i.stack.imgur.com/yXcEp.png) Storing the number 1.23456789123434 as type N (defined in linked document as: "Number stored as a string, right justified, and padded with blanks to the width of the field.") costs me 15 bytes while storing the same number as type O (double) would cost me only 8 bytes. Please correct me if I am wrong on this.


title: How to overlay map onto a corner of an image, using that image's GPS exif data?

text: I have several hundred JPG images taken during a time lapse of a drive, with each image containing embedded GPS information. Is there an OS X friendly workflow that can allow me to:   1. create a highway map image (say a 300x300 pixel Google map representing approx a 3mi x 3mi area) centered around the geo coordinates taken from the embedded GPS exif data of that image   2. overlay that map image onto the corner of the original image, and    3. repeat for all images in the time lapse sequence. I think with a combination of ImageMagick and shell scripting, I can probably do (2) and (3), but I don't know how to do (1).


title: Problem with geotiff in qgis only in some PCs

text: Helllo I'm having some weird problem opening a 16bits geotiff with qgis. In my pc it opens and displays perfectly, but in another pc (with the same version of qgis) it displays all wrong the image getting pixels of several colors without any sense. ![enter image description here](http://i.stack.imgur.com/9RqaR.jpg) Any ideas? Thanks


title: Georeferencing an image in C# with GDAL

text: I want to upload a jpeg or some raster image into a C# form and then on entering the latitude and longitudes I need to be able to create an GEOTIFF file of the same image. How can this be done in C#.NET? you had asked me to use GDAL libraries i did add references, i'm using Visual studio 2010, i have included library files, i have computed llx lly lux luy also, but how to tag it to an image, is there any inbuilt function or anything? this is the code i have used               using System;     using System.Collections.Generic;     using System.ComponentModel;     using System.Data;     using System.Drawing;     using System.Linq;     using System.Text;     using System.Windows.Forms;     using OSGeo.GDAL;     using OSGeo.OGR;     using OSGeo.OSR;          namespace geo_tagging     {         public partial class Form1 : Form         {             public Form1()             {                 InitializeComponent();             }                  private void button2_Click(object sender, EventArgs e)             {                 int orx = pictureBox1.Image.Height;                 int ory = pictureBox1.Image.Height;                 int ctrx, ctry, lux, luy, llx, lly, bufx, bufy;                 bufx = orx / 108000;                 bufy = ory / 108000;                 ctrx = int.Parse(textBox2.Text);                 ctry = int.Parse(textBox3.Text);                 lux = ctrx - bufx;                 luy = ctry - bufy;                 llx = ctrx + bufx;                 lly = ctry + bufy;                      Gdal.GCPsToGeoTransform(             }                  private void button1_Click(object sender, EventArgs e)                  {                 OpenFileDialog dlg = new OpenFileDialog();                           if (dlg.ShowDialog() == DialogResult.OK)                 {                                      PictureBox PictureBox1 = new PictureBox();                 string Chosen_File = "";                 Chosen_File = dlg.FileName;                 textBox1.Text = dlg.FileName;                 pictureBox1.Image = Image.FromFile(Chosen_File);                 int X = pictureBox1.Image.Height;                 int Y = pictureBox1.Image.Width;                 label2.Text = X.ToString();                 label4.Text = Y.ToString();                      }                      dlg.Dispose();             }                  private void button3_Click(object sender, EventArgs e)             {                 textBox1.Text="";                 textBox2.Text="";                 textBox3.Text="";                 pictureBox1.Image.Dispose();                  }                 private void pictureBox1_Click(object sender, EventArgs e)             {                 OpenFileDialog dlg = new OpenFileDialog();                           if (dlg.ShowDialog() == DialogResult.OK)                 {                     PictureBox PictureBox1 = new PictureBox();                     string Chosen_File = "";                     Chosen_File = dlg.FileName;                     textBox1.Text = dlg.FileName;                     pictureBox1.Image = Image.FromFile(Chosen_File);                     int X = pictureBox1.Image.Height;                     int Y = pictureBox1.Image.Width;                     label2.Text = X.ToString();                     label4.Text = Y.ToString();                      }                      dlg.Dispose();             }                           private void textBox1_Click(object sender, EventArgs e)            {                textBox1.Text = "";            }                  }      }      * * * **never mind got it guys, we need to use a exe file to call it as a process and pass arguments into it!** thanks for the help


title: Does anyone know how to get the QGIS Web Client measure tool to show units in feet and miles?

text: Does anyone know how to change the QGIS Web Client measure tool (distance and area) so that it shows feet or miles rather than meters?


title: Why are GeoJSON coordinates different exported from QGIS vs. ogr2ogr?

text: I have a shapefile with a Google Mercator projection (900913) that I need converted to GeoJSON. I need to do this on a regular basis amongst other things so I wrote a shell script to do it using ogr2ogr. I am also converting this to a 4326 projection so the coordinates are more usable. This is how I am using ogr2ogr to do that: `ogr2ogr -f "GeoJSON" -s_srs EPSG:900913 -t_srs EPSG:4326 output.geojson my-shapefile.shp -lco WRITE_BBOX=YES -lco COORDINATE_PRECISION=6 --debug on` For some reason all of my coordinate arrays have 3 values, the last one always being 0.0. Like so: `[ -122.633291, 48.628477, 0.0 ]` Where is the extra 0.0 coming from?? However, if I export from qgis the coordinates come out as expected.. `[ -122.633291, 48.819157 ]` I also notice that the longitude is slightly different between the two. Any ideas on how I could improve my ogr2ogr command? Or perhaps this is a bug? Pete **Here are some more details about the shapefile being converted.** Storage type of this layer:ESRI Shapefile   Geometry type of the features in this layer: Polygon   The number of features in this layer: 431   Editing capabilities of this layer: Add Features, Delete Features, Change Attribute Values, Add Attributes, Create Spatial Index, Fast Access to Features at ID, Change Geometries      In layer spatial reference system units : xMin,yMin -19951913.23,2139102.40 : xMax,yMax 20021888.10,11554793.57      +proj=merc +lon_0=0 +lat_ts=0 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs


title: How to write an address search plugin?

text: I am migrating some users from MapInfo to QGIS, and before they sign off on it, they want a very simple address search in the same way that we currently have one working in MapBasic. I;ve been having a look around, but can;t seem to get my head around modifiying the quick find and find by attribute plugins which already exist. Essentially they currently have the following: A mapbasic routine which presents the user with portions of the address which they can fill in Such as               field 1: Building name/number     field 2: Street Name     field 3: Locality     field 4: Town     field 5: Postcode      The process then queries a specific layer using like statements, returns a list of matching records which when clicked zoom to the relevant address. They are adamant that using the QGIS query builder is going to be 'too complicated' for these users and want me to simplify the process in QGIS to a similar one to the above. So we have a vector layer containing addresses available in the project - is it possible to mirror this MapInfo functionality - and if so, can anyone provide me with some guidance or point to a tutorial for doing this type of task. Thanks Edit: 15:22 11/06/2013 To Clarify: The addresses are held in a MapInfo tab file, which is loaded in to the project. All of the fields that would be queried are of type string - there is an additional useful field called searchAddr which has field 2, 1 and 3 concatenated (split by commas). The data also contains the coordinates as integers for each location. So in terms of code (hopefully, if I use how it works in SQL to explain that will make it clear), we would be doing soemthing like:               Select * from <layer>      where searchAddr like '%<field1>%'      and searchAddr like'%<field2>%'      and searchAddr like '%<field3>%'      and Town like '%<field4>%'      and Postcode like '%<field5>%'      Therefore if the user puts in 42 as the House Number in field 1 and postcode AB12 3CD in field 5 the query would be               Select * from <layer>      where searchAddr like '%42%'      and searchAddr like'%%'      and searchAddr like '%%'      and Town like '%%'      and Postcode like '%AB12 3CD%'


title: OGC WMS Validation

text: We're trying to validate our WMS server with OGC standards,to do so OGC has a test page available. Here,only stacked on a single error                wms:profiles-basic-1 (View Details): Failed      When we look into detail http://cite.opengeospatial.org/test_engine/wms/1.1.1/exec_summary page requires following conditions in order to validate WMS 1.1.1 standards: >   1. It must support image/png or image/gif for GetMap requests.   > >   2. It must contain the CITE dataset in layers that are subsetable and > resizable and support EPSG:4326.   > >   3. The EPSG:4326 support must be precise up to the ten-thousandth of a > degree.   > >   4. The default style for the cite:Lakes layer must fill in the polygon > with non-white pixels. > Condition 1 is clear and is supported on the system.The problem is I believe on the others (2,3,4)   What causes validation test to fail ? Best Regards


title: Modify default options in QGIS

text: I would like to change some default QGIS options. For example, the default font used by QGIS composer on my system is ubuntu. I would like to change that. I would also like to get no border when creating an object. Is there a config file I can create or modify in my .qgis user folder? A user option file (like in scite editor) for example, would be valuable if it does not exists yet.


title: How to split a polygon into multiple small parts by grid

text: we delineated a river basin in shape file and generated a line-featured shapefile using Hawths tools. now our task is that we want to divide this basin into small polygons by gird in order to compute the area within the grid at the boundary. Does anyone know how to implement this using arcgis 9.3? Your help will be greatly appreciated.


title: Can I say WGS84 is more accurate than NAD83 in terms of their reference ellipsoids?

text: Other than its application, can we say the reason of common adoption of WGS84 is due to its accuracy?


title: How to load and view .gri and .grd files in QGIS

text: An online GIS Repository called DIVA-GIS (http://www.diva-gis.org/gdata) offers data in three formats - ERSI files, "DIVA gridfiles", and DBF files The DIVA gridfiles have extensions .gri and .grd How can these file types be opened in qGIS? The "Add Raster" menu option does not seem to work


title: What does feature mean in QGIS?

text: I cant understand the meaning of feature in QGIS `QgsFeature` class. How it is related to layer and geometry? I have read qgis documentation but i cant understand what is feature.


title: Interpolation for timeseries in ArcGIS

text: I'm trying to calculate average rainfall value from a number of points using Inverse Weighted Distance tool in ArcGIS 9.3. My problem is that: each point has it own time series, therefore the interpolation process should be able to carry out for all years (kind of iteration so to speak). Following is a sample attribute table: ID X Y Name Rain1990 Rain1991 Rain1992 Rain1993 .... Rain2010 1 xx1 yy1 AA 1210 1189 1863 1269 ......   2 xx2 yy2 BB 1492 1502 2187 1923 ...... ...... Could anybody show me how to do that. Any suggestion would be much appreciated Nice weekend every1 thecatalyst Edit 1: I finally did this by using C++ code which required ArcGIS mask grid, data files & locations of all the points. If someone needs the code, contact me and I'd love to share! Edit 2: I recently used R to do this interpolation task. You can use either hydroTSM, gstat or spacetime packages. Few example links below: http://spatial- analyst.net/wiki/index.php?title=Spatial_interpolation_exercises_%28NL%29 http://www.geostat-course.org/Topic_Bivand_2012


title: Geowebcache bounding box to tile conversion is not working properly (EPSG:4326)

text: I am trying to make Geoserver and built in Geowebcache work together. I have a page with google maps and I am rendering images from a published WMS layer using getTileUrl function. Unfortunately I am getting wrong tiles when geowebcache is activated. Here is a request which uses caching (caching is enabled for png8):                   http://myserver/geoserver/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.1.1&TILED=TRUE&TRANSPARENT=TRUE&WIDTH=256&HEIGHT=256&FORMAT=image/png8&LAYERS=layername&SRS=EPSG:4326&bbox=28.7841796875,40.94671366508002,28.828125,40.979898069620155      response headers:                   Cache-Control:no-cache         Content-Disposition:inline; filename=layername.png8         Content-Type:image/png         Date:Fri, 30 May 2014 08:19:17 GMT         ETag:91591a0d685881053a9d4761f42086f8         geowebcache-cache-result:HIT         geowebcache-crs:EPSG:4326         geowebcache-gridset:EPSG:4326         geowebcache-layer:layername         geowebcache-tile-bounds:28.7841796875,40.95703125,28.828125,41.0009765625         geowebcache-tile-index:[4751, 2980, 12]         Last-Modified:Thu, 29 May 2014 16:31:31 GMT         Server:Apache-Coyote/1.1         Transfer-Encoding:chunked      Here is a second request which hits directly to geoserver (caching is enabled for jpeg):                   http://myserver/geoserver/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.1.1&TILED=TRUE&TRANSPARENT=TRUE&WIDTH=256&HEIGHT=256&FORMAT=image/jpeg&LAYERS=layername&SRS=EPSG:4326&bbox=28.7841796875,40.94671366508002,28.828125,40.979898069620155      response headers:                   Content-Disposition:inline; filename=layername.jpg         Content-Type:image/jpeg         Date:Fri, 30 May 2014 08:19:01 GMT         geowebcache-cache-result:MISS         geowebcache-miss-reason:no tile cache for requested format         Server:Apache-Coyote/1.1         Transfer-Encoding:chunked      Although I am requesting for the same bounding box, the first request which uses geowebcache returns wrong portion of the map while the second works fine. Also, the first request finds a tile in the 12th zoom level while my map is actually in zoom level 13. I guess, bbox -> tms tile mapping conversion is not working properly. Anyone can propose suggestions on how I am doing wrong? Note: my wms layer has native & declared srs:4326 Note 2: I am not using openlayers. Edit: I can check I am using the correct bounding box using the following url:                  http://www.maptiler.org/google-maps-coordinates-tile-bounds-projection/        Google: (4751,3072)        TMS: (4751,5119)        Zoom:13             == Tile Boundaries ==        WGS84 datum (longitude/latitude):        28.7841796875 40.94671366508002        28.828125 40.979898069620155      But webgeocache finds [4751, 2980, 12]. In fact when I look at webgeocache directory tile x=4751 exists only in level 12 but not in level 13. Edit 2: Example for topp:states layer Google: (301,382) TMS: (301,641) Zoom:10                   http://localhost:9090/geoserver/gwc/service/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.1.1&TILED=TRUE&TRANSPARENT=TRUE&WIDTH=256&HEIGHT=256&FORMAT=image/png8&LAYERS=topp:states&SRS=EPSG:4326&bbox=-73.828125,41.50857729743935,-73.125,42.032974332441405              response headers:         Cache-Control:max-age=3600, must-revalidate         Content-Disposition:inline; filename=geoserver-dispatch.image         Content-Type:image/png         Expires:Fri, 30 May 2014 10:47:17 GMT         geowebcache-cache-result:HIT         geowebcache-crs:EPSG:4326         geowebcache-gridset:EPSG:4326         geowebcache-tile-bounds:-73.828125,41.484375,-73.125,42.1875         geowebcache-tile-index:[151, 187, 8]         geowebcache-tile-index:[151, 187, 8]         Last-Modified:Fri, 30 May 2014 09:43:38 GMT         Server:Jetty(6.1.8)         Transfer-Encoding:chunked              Correct tile should be:         http://localhost:9090/geoserver/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.1.1&TILED=TRUE&TRANSPARENT=TRUE&WIDTH=256&HEIGHT=256&FORMAT=image/jpeg&LAYERS=topp:states&SRS=EPSG:4326&bbox=-73.828125,41.50857729743935,-73.125,42.032974332441405


title: Extract boundary from detailed layer in ArcGIS or QGIS

text: Is it possible to return the boundary (outside) permiter of a polygon layer? Let's say you have the map of the USA, for simplicity, broken down to the county level, is it possible to return the boundary map of the whole country from that map?


title: Transactions in Geoserver / WFS - Is it possible to know which "action" failed?

text: I'm trying to implement transactions for some queries to Geoserver. I'm using plain HTTP Post to make the requests, from Java. One of those queries looks something like :               <wfs:Transaction service='WFS' version='1.0.0'                                 xmlns:wfs='http://www.opengis.net/wfs'                                                        xmlns:gml='http://www.opengis.net/gml'                                     xmlns:ogc='http://www.opengis.net/ogc'                                     xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance'                      xsi:schemaLocation='http://www.opengis.net/wfs                                                 http://schemas.opengis.net/wfs/1.0.0/wfs.xsd'>                                                                       <wfs:Insert>                                                                  [[FEATURE_TO_INSERT_1]]                                              </wfs:Insert>                                                              <wfs:Insert>                                                                  [[FEATURE_TO_INSERT_2]]                                              </wfs:Insert>                                                              <wfs:Insert>                                                                  [[FEATURE_TO_INSERT_3]]                                              </wfs:Insert>                                                               </wfs:Transaction>          If all three inserts succeed, they are commited. If one fails, none are commited and something like this is returned :               <?xml version="1.0" ?>     <ServiceExceptionReport        version="1.2.0"        xmlns="http://www.opengis.net/ogc"        xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"        xsi:schemaLocation="http://www.opengis.net/ogc http://schemas.opengis.net/wfs/1.0.0/OGC-exception.xsd">        <ServiceException>           java.lang.RuntimeException: Parsing failed for DIAMETRE: java.lang.NumberFormatException     Parsing failed for DIAMETRE: java.lang.NumberFormatException     null     </ServiceException></ServiceExceptionReport>      This is fine, but I'd like to know _which one of those three inserts actually failed_! Is there a way to give each action of the transaction an `id`, and make Geoserver output it somewhere in the returned `ServiceExceptionReport` so I'm able to identify which action the error is related to? Update : I'm currently using WFS 1.0.0. Maybe there are more options for transactions with a newer version?


title: recover toolbar in visualizer qgis 1.8.0 lisboa

text: I accidentally removed the toolbar in the printer visualizer. Anyone knows how to recover it? I would expect a button of show toolbar or something, but I can not find it anywhere... Thanks!


title: Leaflet rendering TileMill tiles in the wrong order

text:
I'm stumped, and quite a lot of searching has turned up nothing helpful. I'm on the first step of what will ultimately be a map with three raster layers and three layers of markers (with controls to toggle between the rasters and their associated data, with popups/tooltips for the data), but I can't get the Tilemill tiles to render properly. (I created the rasters in QGIS, ran them through GDAL to colorize and re-project them, and then put them into TileMill.) To wit: ![enter image description here](http://i.stack.imgur.com/xgFxM.jpg) That's in Chrome. This is Safari: ![enter image description here](http://i.stack.imgur.com/WU4Wn.jpg) And here's the code, which looks quite innocuous to me (I've commented out everything but the tiles at this point, so that's all I've included):                          <title> - Leaflet demo</title>                                Map

                Has anyone ever encountered this problem? If so, how did you overcome it? I'm eager to start into the rest of the process of making this map (by which of course I mean "uncovering the next problem to solve"). _Edited to add:_ I tried re-exporting a new MBTiles file, uploading it as new data, and creating a new Mapbox project from that data. No dice. Next step in isolating the issue would be to set up the MBTiles file on a server, skipping Mapbox, but that seems more likely to create problems than to fix them.


title: How can the scale be projected correctly in Quantum GIS when it is imported from google earth (KML) files?

text: While importing KML files from google Earth to Quantum GIS, the scale is not being projected properly and the area of a particular selected feature is not being calculated in the attribute table by the field calculator. How can the area of a selected feature be calculated by the field calculator, as right now it is showing up zero as the value?


title: ArcGIS JSAPI draw toolbar tooltips not displayed in IE10

text: I am experiencing an issue where draw toolbar tooltips are not displayed in IE10 when using the ArcGIS JSAPI. Try the following example in IE10 and chrome (draw a polygon to see what I mean) http://developers.arcgis.com/en/javascript/samples/graphics_add/ Has anyone else experienced this issue and is there any way to get the tooltips to display in IE10?


title: What is Spatial Reference System in Laymans Language?

text: I am trying to understand the concept behind the Spatial Reference System, Why is it used and what is so specific about the number. I went through many definations and almost all went over my head. So Could any one help me understand what is SRID, and its importance in GIS. I thought, it was a reference to draw maps from a real world place, expressed in map. But then I inserted a point in PostGIS, with SRID 300 and then inserted the same point with SRID 350, I found the two points showing the same place, in QGIS viewer. I am totally confused.


title: Sending GPS Location to server

text: Is it possible to send my current location on iPhone or Android to the server so that others can view my current location?


title: Does ArcGIS (arcObjects SDK for java) support multiline labels using TextSymbol?

text: Is it **possible** to add multline labels using textsymbols in ArcGIS 10.1 (SDK for Java) in a 3d model? In my system, even the spaces are getting omitted, for a String `s="this is \n war"` Label produces output:`thisiswar` any help or alternative would be much appreciated. thank you


GetFeature with WFS layer is not working as the click co-ordinates are returned in spherical co-ordinates. using geoserver 2.1.3 and openlayer 2.11 layer information is in EPSG:4326 WMS works fine WFS trasaction works fine trying to GetFeature information on click from the WFS layer as the total number of features are small and it is a custom vector layer.               control = new OpenLayers.Control.GetFeature({         protocol: wfs_layer_protocol,         box: false,         hover: false,         click: true     });          control.events.register("featureselected", this, pickNodeId);     map.addControl(control);     control.activate();      the WFS post request has the following bounds -14710846.457064 4317003.9227542-13710750.910779 6517099.4690396


title: ArcIMS vs ArcGIS Server

text: having never used ArcIMS I am not able to answer the following question. "Why don´t we use ArcIMS anymore?" I have recently taken over a position which used to administer ArcIMS clients but now use ArcGIS Server architecture. Colleagues who have been here a while still ask why ArcIMS needs to be replaced and the only answer I can give is that the software is not supported or developed by ESRI anymore. But regarding stability, functionality,supported formats etc I can´t say much. Could anyone give me a concise little list of reasons why ArcIMS is "old" and should be replaced by ArcGIS Server? It would also help me to understand how Server GIS in ESRI has developed. Thanks for any comments, Robert


title: How can i detect the number of polygon is drawed in a edit session in QGIS

text: I am working with qgis 1.8 and i want to know how many polygon drawed in one edit sessions is there any signal emitted with qgis ???


title: Validate a gml file

text: I receive a gml-like response from a RESTful Web service. It looks good for me but I can not imported in QGIS. Is there an online validation tool for a valid gml file? Or am I missing something for the schema? The response I received is as               <Polygons>     <DriveTime duration="30">     <Polygon>     <exterior>     <LinearRing>     <posList>41.0296010971069 -73.9398550987244 41.0188937187195 -73.9155864715576 41.0581398010254 -73.8332319259644 41.017370223999 -73.8251423835754 40.9897971153259 -73.7459421157837 40.9789824485779 -73.6955809593201 40.9320545196533 -73.7493538856506 40.862123966217 -73.8021397590637 40.8192944526672 -73.8084053993225 40.8174920082092 -73.8013029098511 40.8047890663147 -73.6684370040894 40.7836532592773 -73.5890650749207 40.7572603225708 -73.6199641227722 40.7433986663818 -73.6648964881897 40.7084226608276 -73.6869120597839 40.6833171844482 -73.6838436126709 40.6573748588562 -73.6869549751282 40.6232357025146 -73.7346339225769 40.6431698799133 -73.7853598594666 40.6136870384216 -73.8205289840698 40.6405091285706 -73.8775634765625 40.629415512085 -73.8853740692139 40.6116271018982 -73.897111415863 40.5840969085693 -73.9238047599792 40.5747413635254 -73.9734792709351 40.58349609375 -73.9834785461426 40.6026148796082 -74.015429019928 40.5769944190979 -74.0807676315308 40.5689263343811 -74.0921187400818 40.5981302261353 -74.1254639625549 40.6105327606201 -74.1812968254089 40.5980658531189 -74.2316579818726 40.6317758560181 -74.238395690918 40.680570602417 -74.2200064659119 40.6869649887085 -74.2736077308655 40.714066028595 -74.295859336853 40.7720446586609 -74.269380569458 40.8061408996582 -74.2862892150879 40.8088231086731 -74.290087223053 40.8478331565857 -74.2046213150024 40.8865427970886 -74.2311215400696 40.8948683738709 -74.2362928390503 40.9156179428101 -74.1825199127197 40.9354877471924 -74.133939743042 40.950014591217 -74.1370940208435 40.9720087051392 -74.080445766449 40.9833598136902 -74.0716052055359 40.9472250938416 -73.9906239509583 40.9620308876038 -73.9610767364502 41.0296010971069 -73.9398550987244</posList>     </LinearRing>     </exterior>     </Polygon>     </DriveTime>     </Polygons>


title: Can't remove the little islands when creating a svg map with kartograph.py

text: I'm new to python and to kartograph. I'm using it to create a map of France, Italy, Portugal and Spain. But I don't want the little islands (mainly the French ones ), causing the generation of a big useless map, for me. I used a python script to create my map ( thanks to rmbianchi ):               # file: makeMap.py     from kartograph import Kartograph     K = Kartograph()          def myfilter(record):         return record['iso_a3'] in ["FRA","ITA","ESP","PRT"]               config = {         "layers": {             "FTC": {                 "src": "ne_50m_admin_0_countries.shp",                 "filter": myfilter,                 "attributes": {"iso_a3":"iso_a3", "name":"name", "id":"iso_a3"},             "bounds": {         "mode": "polygons",         "data": {             "layer": "FTC",             "filter": {                 "filter": ["ISO_A3", "in", ["FRA", "ITA", "ESP", "PRT"]]             },             "min-area": 0.50         }     }             }         },      }          K.generate(config, outfile='ftc.svg')      But, little islands are still displayed, how can I avoid that ? Thanks in advance, if I forgot to give some informations, don't hesitate to ask ! Regards.


title: How to calculate/record attachments' file names using Field Calculator?

text: I have a point feature class containing attachments corresponding to jpg files (sometime more than one jpg per feature). I would like to use the "field calculator" to calculate/record a new field that would contain file paths or names of attachment files. Does someone know how to do that? I am new to Python and VB script...


title: British National Grid to KML results in offset polygons

text: I have a shape file of English wards from ONS. It is similar to the files available here: http://www.ons.gov.uk/ons/guide- method/geography/products/census/spatial/2011/index.html except I had to explicitly ask them for ward shapes, for which they have just released census data. They sent me a shapefile, and I would like a KML. The shapfile lines up perfectly (I use MapWindowGIS to view it) with The other shape files in that set, suggesting there is nothing wrong with the coordinates of this one file. When I convert it to KML, all the coordinates are offset. I use Google Earth to view it, and it does not line up with the other KML files (for example MSOA). In both cases I am looking at the clipped to the coastline versions. I am specifically looking at Lambeth, and I can see that the border does not line up with the Thames and is in land too far east and and possibly a bit too far south. To convert from shapefile to kml I have tried postgis, a Shape2Earth plugin for map window GIS and an appplication called shp2kml. In all cases I get the same result. I have tried converting the Shapefile from BNG to WGS84 then converting that to KML. I also tried converting to ETRS89, and then to WGS84 as suggested in one of the answers here: Converting from British National Grid (27700) to WGS84 (4326) to Web Mercator (102100) My prj file looks like this:               PROJCS["British_National_Grid",     GEOGCS["GCS_OSGB_1936",         DATUM["D_OSGB_1936",             SPHEROID["Airy_1830",6377563.396,299.3249646]         ],         PRIMEM["Greenwich",0.0],         UNIT["Degree",0.0174532925199433]     ],     PROJECTION["Transverse_Mercator"],     PARAMETER["False_Easting",400000.0],     PARAMETER["False_Northing",-100000.0],     PARAMETER["Central_Meridian",-2.0],     PARAMETER["Scale_Factor",0.999601272],     PARAMETER["Latitude_Of_Origin",49.0],     UNIT["Meter",1.0]     ]      Any idea what I am doing wrong?


title: Clip for each feature

text: I want to clip a polyline feature class by a polygon feature class (buffers created from points). Several polygons overlap, so I want to get each polygon feature clipped with each polyline feature. Using the Clip_analysis tool, I get the orange lines: ![enter image description here](http://i.stack.imgur.com/Nh9g5.png) whereas I want the following red lines: ![enter image description here](http://i.stack.imgur.com/7jxe9.png) In a python script, I tried the following code to clip every polyline with each polygon:                       temp = arcpy.ValidateTableName("clip_temp")         bufferCursor = arcpy.SearchCursor(buffer)         n = 0         for pRow in bufferCursor:             arcpy.AddMessage(n)             feat = pRow.Shape             arcpy.Clip_analysis(inputPolylineLayer, feat, temp, "#")             if n == 0:                 arcpy.Copy_management(temp, clip)             else:                 arcpy.Append_management(temp, clip)             n=n+1      But the execution is very very slow and it crashes before finishing. I also tried the Intersect_analysis tool but I didn't get the right geometries: ![enter image description here](http://i.stack.imgur.com/ohJnD.png) Do you have any idea of a better methods to do it?


title: Reprojecting raster data from a rotated pole grid

text: I have some raster data (ASCII grid files) that use a custom rotated-pole projection, that I'd like to reproject into something more useful, such as lon/lat WGS84. A description of the grid projection is here, but I can't see how I can use this to create a projection string for use in GIS. http://ukclimateprojections- ui.defra.gov.uk/ui/docs/grids/prob_land_25km_rotated/index.php This SE question seems to discuss using rotated pole data, but I can't work out how to apply it to my problem: Lon/lat transformation I found a thread here asking exactly same question (about the same dataset), but with no replies: http://lists.osgeo.org/pipermail/gdal- dev/2011-March/028174.html The header of the ASCII grid data looks like this:               ncols 39     nrows 52     xllcorner -0.50     yllcorner -0.50     cellsize 1.00


title: Make two layers editable in leaflet draw

text: I use Leaflet draw to draw markers and polygons on a map. The created features are saved to two separate feature groups. One feature group for the markers and one for the polygons. Now I want to edit the features. How can I add two feature groups as editable layers:               'edit': {       'featureGroup': featureGroupMarker, featureGroupPolygon     }


title: Problem when executing Dissolve method in ArcEngine

text: I'm always having a problem in GP.execute() when trying to execute this method. Here is my code:                Friend Function DissolveFeatures(ByVal pInFLayer As ESRI.ArcGIS.Carto.IFeatureLayer, ByVal DissolveField As String, ByVal OutName As String) As ESRI.ArcGIS.Carto.IFeatureLayer                      Dim GP As New ESRI.ArcGIS.Geoprocessor.Geoprocessor                 Dim Dissolve As ESRI.ArcGIS.DataManagementTools.Dissolve = New ESRI.ArcGIS.DataManagementTools.Dissolve                 Dim pFLayer As ESRI.ArcGIS.Carto.IFeatureLayer = New ESRI.ArcGIS.Carto.FeatureLayer                      Try                     GP.AddOutputsToMap = True                     GP.OverwriteOutput = True                     GP.SetEnvironmentValue("workspace", "Database Connections\IFE connexion.sde")                          Dissolve.in_features = pInFLayer                     Dissolve.dissolve_field = DissolveField                     Dissolve.out_feature_class = "Database Connections\IFE connexion.sde\" & OutName                          RunTool(GP, Dissolve, Nothing)                               pFLayer.FeatureClass = CType(GP.Open(OutName), IFeatureClass)                     pFLayer.Name = pFLayer.FeatureClass.AliasName                      Catch ex As Exception                     MsgBox(ex.Message & "Dissolve Features")                 End Try                      Return pFLayer                  End Function               Private Shared Sub RunTool(ByVal geoprocessor As ESRI.ArcGIS.Geoprocessor.Geoprocessor, ByVal process As IGPProcess, ByVal TC As ITrackCancel)                  ' Set the overwrite output option to true             geoprocessor.OverwriteOutput = True                  Try                 geoprocessor.Execute(process, Nothing)                 ReturnMessages(geoprocessor)                  Catch err As Exception                 Console.WriteLine(err.Message)                 ReturnMessages(geoprocessor)             End Try         End Sub                        ' Function for returning the tool messages.         Private Shared Sub ReturnMessages(ByVal gp As ESRI.ArcGIS.Geoprocessor.Geoprocessor)             ' Print out the messages from tool executions             Dim Count As Integer             If gp.MessageCount > 0 Then                 For Count = 0 To gp.MessageCount - 1                     Console.WriteLine(gp.GetMessage(Count))                 Next             End If         End Sub      This returns an error: `Error HRESULT E_FAIL has been returned from a call to a COM component` Please help.


title: Problems closing dbf-table edited with ArcObjects

text: To change values and add new values to a dbf-file, I connect to the dbf and use the ITable interface to make the edits. At the end I release the COM- objects of the table and the workspace. But when I try to open the changed file, it says the file is read-only. Did I forget to release something else? Or is there another reason for the read-only status? Thanks and regards Fab


title: How to update/delete data from a datagridview?

text: I wanna be able to update/delete data from a datagridview, and when I do that of course the data in my PostgreSQL database should be edited too. So I've thought about giving the user the ability to select a layer, choose an attribute, choose a logical operator (=, < or >) and then type the criteria he wants for his search. So for example he could make a search like this one : > Layer : Activity > > Attributes : Id > > Logical Operator : = > > Criteria : 15 And when he clicks on Search, the result should display on the datagridview The thing here is, if I want to update a row or a number of rows, how can I be able to do that? Same thing for deleting? Is it possible to edit data in the datagridview and then I don't know do like an update button and save the changes accordingly but how am I gonna do that ? # **EDIT** So now I have this idea : I have two forms : > **Form 1 : Activity** contains : > >>   * the datagridview linked to PostgreSQL database (and the DB is linked to ArcCatalog) >>   * a button named Update on which I can click and show the Form2 where I can edit my data (I pass the values from the datagrid view to my textboxes fields in Form2) >>   * a button named Delete on which I can click and delete the IFeature. >> > > **Form 2 : UpdateActivity** contains : > >>   * the fields of the selected Ifeature, with the values already passed >>   * A Validate button on which I can click and update the data after I finish editing, also the form should close and then I see the updated data in the datagridview in form 1 >> >>> **P.S** : How can I call the datagridview in the second form? >>> Here's some code : On Activity Form:               // This is the Update button that pop up the Form2, and it's where I pass the values     private void btnUpdate_Click(object sender, EventArgs e)         {             updf = new UpdateForm();                  //Loading data from datagridview to UpdateForm             updf.txtIdActivity.Text = this.TableAttributaire.CurrentRow.Cells[0].Value.ToString();             updf.txtNameActivity.Text = this.TableAttributaire.CurrentRow.Cells[1].Value.ToString();             updf.txtAddressActivity.Text = this.TableAttributaire.CurrentRow.Cells[2].Value.ToString();             //Update Data             updf.ShowDialog();               }      On Activity Form:               //Clicking on this button should I create another form like I did with the update ?     //Or just delete straight from the datagridview with a messagebox?     private void btnDelete_Click(object sender, EventArgs e)         {             //Haven't done it yet?         }      On UpdateActivity Form:               //This is the validate click button on which I can click and save and update the data     private void Validate_Click(object sender, EventArgs e)         {             IFeatureClass pFeatClass;             IQueryFilter pQF = new QueryFilter();             // assume pFeatClass is an IFeatureClass object             // and pQF is a queryfilter that allows to select the features in question.             IFeatureCursor pCur = pFeatClass.Update(pQF, true);              IFeature pFt = pCur.NextFeature();             pFt.set_Value(indx, pValue); // I don't know what to do here ??             pCur.UpdateFeature(pFt);       //But where is the datagridview in all this? How can I call it in this second form and update tha data ?     //Also like I said before I should be able to close this form after I validate and refresh the datagridview in Activity Form         }      **Updating/Deleting from the map** : > How can I select an object on my map in Activity Form, and display the > Update Form accordingly and be able to edit the data? Same thing for > deleting?


title: Calculating Map Resolution Without Pre-Defined Scale Levels

text: I am trying to determine the formula necessary to calculate map resolution (change in lat,long per screen unit). I am interested in only WGS84 coordinate space. Every formula I have found thus far seems dependent on a map that has predefined zoom levels, i.e. a tiled map. For example, the common formula I find is: Map resolution = 156543.04 meters/pixel * cos(latitude) / (2 ^ zoomlevel) How can I adjust this formula for a map without zoom levels?


title: ArcGIS 10 Add-In compile problem - Why are new edits to code igonred by ArcMap?

text: I've been compiling and debugging my Add-In (a toolbar with a few menus, custom commands, custom selection tool and a context shortcut menu) with no issues until **suddenly** the edits and changes I do in Visual Studio are reflected in the Add-In. There are no errors and compile success. I first noticed this when I was wanting to debug but application was ignoring new edits and additions. The Add-In is still there but some commands don't seem to not respond to current changes. What am I doing wrong?


title: Geoserver 2.4.0 via the web admin interface, not adding layers in Internet Explorer, but working in Firefox

text: For some reason I cannot seem to add layers via Internet Explorer, but I can via Firefox 23. In addition the styles that I upload also appear broken in Internet Explorer 8, but do not in Firefox. Is there any sort of modification to enable it to work in any browser?


title: Statistics to shades in a world map

text: Is there any API or web service to create statistical world maps where the values are different shades of a colour? Something like http://www.visualizingeconomics.com/wp-content/uploads/sachs- gdpdensity.png


title: Compiling a list of units that border each other

text: I have a shapefile of U.S. counties. I am interested in compiling a list of those that border each other. In particular, I am interested in those that border each other but are in different states. For example, I am interested in knowing that Rabun County borders counties in North Carolina and which ones and not the counties in Georgia that it may also border. I work with ArcGIS 10.1 and use ArcMap. I see that there is a feature called Polygon Neighbors. However, I don't fully understand how to use it. I'm not sure what to put into the "Input Features" or how to proceed from here. I would very much appreciate any help!


title: How can I change the CSS style of the geoserver rendered map produced when I click on "Open Layer "Tab under "Common format"?

text: I have tried many times with left no stone unturned, but all efforts came in vain. I might be of less knowledge in this area. Then through surfing I found this site as a newbie. I will be very obliged if I get a solution to the below mentioned scenarios. In GeoServer when I click on "layer preview" button then "Open Layer" for preview, rendered map aligns to the left. I need to align center and give some CSS style and others to the rendered page as below.   1. Left side panel for navigational purpose through layers made from PostGIS data.   2. Middle side is intended map.   3. Right side is for showing some attribute table for the left panel chosen layer ![enter image description here](http://i.stack.imgur.com/LuTA3.jpg) How can I do these?


title: Wavelet vs Discrete compression for imagery

text: What is the main reasons to choose wavelet (jp2, sid, ecw) instead of discrete (JPEG) compression methods for raster data? Lets limit this question to aerial imagery and lossy compression. Why should we give preference to wavelet compression instead of "classic" JPEG? It will give smaller size but may be there any drawbacks, what do you think? Will it lead to more CPU processing for decomression (during viewing)?


title: Difference between Albedo and Surface Reflectance

text: Can anyone tell me what is the difference between Albedo and surface reflectance? I know for sure that they are not the same. Any links / comments are most welcome. Thanks GBh


title: qgis: module crashed or killed

text: using Grass 6.4.3. via Qgis 2.2.0 Valmiera (standard installation package): Why would Grass modules suddenly produce a "module cashed or killed" error. specifically in relation to r.surf.idw2 and r.surf.contour (i havent managed to go further into the workflow) (it was fine a few hours ago)


title: Calculated distance doesn't match Google Earth

text: I have a google earth route that measures 190.1km in Google Earth. I then pulled down the .kml file, extracted the points and measured them using Vincenty's formula (ellipsoid earth) and got 190.2. BUT, if I measure the SAME points along the SAME route to a different distance along that route, I'm way off. Sometimes as much as 2km. For example, If I want to know the lat and lon of the point closest to 100km, I measure the distance between all the points along the route until the measurement is as close to 100km as it's going to get. When I compare that point to the route in Google Earth, it's around 98km. Any ideas of why this is the case?


title: How to create atlas along line/route in QGIS

text: I'd like to create a mapbook along a given route that is designated by a single linear feature. Is it possible to have QGIS (2.2) automaptically compose a series of pdfs that over lap by a given percent and rotate the map to optimize space? Thanks.


title: Is there a way to force ArcCatalog to use a specific tessellation scheme for a SQL Server spatial index?

text: If I Recalculate a spatial index from ArcCatalog, the tessellation scheme is always set to "Geometry Grid". Even when the previous index was Geometry auto grid. ArcCatalog always changes it back. Is there a way to specify the tessellation from ArcCatalog? Alternatively, can SQL Server default to a specific tessellation scheme? _Microsoft SQL Server 2012 - 11.0.2100.60 (X64)   ArcGIS Desktop 10.1 SP1_


title: Using GDAL to read dted files

text: How do I get the elevation at a certain point using gdal to read in a dted file? I don't see much useful information out on the web. I'm sorry I should have been more specific. How would I accomplish this in python?


title: What causes ModelBuilder error occurred in script on this page MdDlgContent.htm?

text: I've been creating a model in Arc 10.0 for a few days now without problems. Then I added the Calculate Field Tool. I selected Make Variable --> From Parameter --> Expression. When I check off "Model Parameter" on the expression variable once it's been added to the model, save and close the model, then try to open it in Toolbox, I receive two popup errors, both saying "An error has occurred in the script on this page." The first error says Error: Expected ')' in Line 79 Character 70 (Code 0): ![enter image description here](http://i.stack.imgur.com/BSFW5.jpg) The second says The value of the property 'IntForm' is null or undefined, not a Function object in Line 465 Character 1 (Code 0): ![enter image description here](http://i.stack.imgur.com/Aey0I.jpg) Both messages are referencing the following URL: file:///C:/Users/MYUSERNAME/AppData/Roaming/ESRI/Desktop10.0/ArcToolbox/Dlg/MdDlgContent.htm. I located this article in ESRI support and followed the instructions, to no avail. (It must not be the exact same issue, as my message box doesn't reference Internet Explorer in the title). I also tried downloading SP5. I opened the MdDlgContent.htm to take a look but that wasn't much help either because, despite some familiarity with Python and JavaScript, I am hardly a programmer. Any help troubleshooting or resolving this issue would be greatly appreciated!! I posted this same question in the ESRI Forums as well: http://forums.arcgis.com/threads/79448-Model-Builder-An-error-has-occurred-in- the-script-on-this-page.-MdDlgContent.htm


title: Having problems creating VRT in QGIS from png files

text: This is similar to a question asked by someone else, re OS tiles, however, there is a difference! I am a bit of a amateur in GIS and python, but if I take it slow I can usually just about work it out! This one though: I have a large amount(over 10k) of geo-refd png files which I have mosaicked in Arc, but also want to be able to access them in other programs such as QGIS over the web and within a private network using Geoserver - so I have tried creating VRT's with them, and sometimes it works, but sometimes it doesn't: it seems that on selecting over 200-odd png's, the vrt stops being created: is there a limit on the amount of rasters for a VRT, or can some-one give me some help please!! This is so frustrating ;o) The error message is: " **The process failed to start. Either the invoked program is missing, or you may have insufficient permissions to invoke the program.** " and the code is:                   gdalbuildvrt -allow_projection_difference          "F:/Historic Buildings/GIS/Data/Mapping/Landranger/Raster/15k/Data/SX/SX.vrt"          "F:/Historic Buildings/GIS/Data/Mapping/Landranger/Raster/15k/Data/SX/1.png"          "F:/Historic Buildings/GIS/Data/Mapping/Landranger/Raster/15k/Data/SX/2.png"          "F:/Historic Buildings/GIS/Data/Mapping/Landranger/Raster/15k/Data/SX/3.png"          "F:/Historic Buildings/GIS/Data/Mapping/Landranger/Raster/15k/Data/SX/4.png"          "F:/Historic Buildings/GIS/Data/Mapping/Landranger/Raster/15k/Data/SX/5.png"          "F:/Historic Buildings/GIS/Data/Mapping/Landranger/Raster/15k/Data/SX/6.png"          "F:/Historic Buildings/GIS/Data/Mapping/Landranger/Raster/15k/Data/SX/7.png"      with an additional +200 similarly addressed png's I am using QGIS 1.8.0 Lisboa, on W7 64bit Any help would be gratefully received


title: OpenLayers - vector boxes and paths

text: I'm using OpenLayers to make warehouse map. How to prevent drawing lines (paths) over some vector objects already on the map?


title: Arctoolbox and python

text: I am trying to use arctoolbox and python so that the user can enter in parameters. In my script I have the following lines               _StationInterval = arcpy.GetParameterAsText(0)      I just want the user to be able to enter in a value for the station interval. However I get the error:               <type 'exceptions.NameError'>: name '_StationInterval' is not defined     Failed to execute (Final3Test).      Before I used arctoolbox, I had the station interval set to a value in the script, i.e. _StationInterval = 10000 Am I missing something in my script??