Layer set visible

I add this layer like this:               var POIData = new esri.layers.ArcGISDynamicMapServiceLayer("https://mysite.com/ArcGIS/rest/services/POI_Data/MapServer");                         map.addLayer(POIData);                         POIData.hide();      I want to show it from a button: This is the button:
POI
      And this is the function behind it:               function showPOI(){     POIData.show();      } What am i doing wrong? Any ideas. I just now cant see it, and cant set it to visible. Thanks in advance


Spatial join to kernel density / kernel density with polygon/point output in Arc

I'm looking for a way to join polygons to a density function for some other feature of interest. Current workflow   * make density raster from feature set of interest, using point/line/kernel density   * sum raster inside output geometry (for polygons) or just sample it (for points) This seems inefficient as I have to calculate a density raster for many points I won't even use. Is there a tool that does the same thing but without the raster? This would be like a kernel-density-aware spatial join, or a kernel density function that outputs to a preexisting geometry rather than to a raster.


GeoExt (ExtJs + OpenLayers): bug in IE9 when try to change map_panel from container A to container B

Techonology: GeoExt, which combines ExtJS with OpenLayers I am trying to dynamically change the app layout by changing a map_panel from container A to container B, and vice verse; basically, for each operation (double click panel header in following codes), the original container A will hide and remove this map_panel from it, then the new container B will add this map_panel object as an item and show up. The following code snippet works fine in Chrome15 and FF8, but in IE9, it has following errors:               SCRIPT5007: Unable to get value of the property 'removeChild': object is null or undefined       OpenLayers.js, line 825 character 412 I am wondering what else I need to add or config in order to make this work in IE as well... any advice are much appreciated...... codes (copy&run):               <!DOCTYPE HTML>     <html>     <head>         <meta http-equiv="Content-Type" content="text/html; charset=utf-8">          <script type="text/javascript" src="http://extjs.cachefly.net/ext-3.2.1/adapter/ext/ext-base.js"></script>         <script type="text/javascript" src="http://extjs.cachefly.net/ext-3.2.1/ext-all.js"></script>         <link rel="stylesheet" type="text/css" href="http://extjs.cachefly.net/ext-3.2.1/resources/css/ext-all.css" />         <script src="http://www.openlayers.org/api/2.10/OpenLayers.js"></script>         <script type="text/javascript" src="http://api.geoext.org/1.0/script/GeoExt.js"></script>     </head>     <script type="text/javascript">     var map, layer;          map = new OpenLayers.Map( ); //define map object     layer = new OpenLayers.Layer.OSM( "Simple OSM Map");     map.addLayer(layer);     map.setCenter(         new OpenLayers.LonLat(-71.147, 42.472).transform(             new OpenLayers.Projection("EPSG:4326"),             map.getProjectionObject()         ), 4     );                   Ext.BLANK_IMAGE_URL = 'img/s.gif';          Ext.onReady(function() {         Ext.QuickTips.init();              var mapPanel = { //put map object in a Panel, using GeoExt class 'gx_mappanel'             xtype : 'gx_mappanel',             id : 'MAP_PANEL',             map : map,             zoom: 6         };              //define viewport         var viewPort = new Ext.Viewport({            layout: 'border',             items: [             { //east container for mapPanel                 title: 'Map',                 id: 'mapEast',                 region: 'east',                 layout: 'fit',                 width: 200,                 border: true,                 frame: true,                 hidden: true,                 items: []             }, { //center container for mapPanel                 title: 'Map',                 id: 'mapCenter',                 region: 'center',                 frame: true,                 layout: 'fit',                 items: [mapPanel]             }]         });              var mapCenter = Ext.getCmp("mapCenter");         var mapEast = Ext.getCmp("mapEast");              mapCenter.header.on("dblclick",function(){ //double click header of mapCenter, map goes to east Container, and centerContainer hides             mapCenter.hide();             mapCenter.removeAll();                  mapEast.add(mapPanel);             mapEast.show();             viewPort.doLayout();         });              mapEast.header.on("dblclick",function(){ //double click header of mapEast, map shows in center Container, east container hides             mapEast.hide();             mapEast.removeAll();                  mapCenter.add(mapPanel);             mapCenter.show();             viewPort.doLayout();         });     });      </script>     <body>     </body>     </html>


Deploying to ArcGIS 10 Server .NET/Java

We're in the process of transforming a lot of 9.1 ArcServer/ArcObjects etc applications to ArcGIS Server 10. However, we're unsure as to licensing on 10. As I understand it, I cannot, according to license restrictions, deploy a toolbox that has been written using ArcObjects 10 .NET SDK on the ArcGIS 10 Server for Java platform, for licensing issues. We _can_ deploy a .NET application to an ArcGIS Java Server 10 instance, if we bundle up all of the application object references in the installer, but this goes against the license agreements (you're not allowed to bundle them up in the installer apparently), so these have to be deployed on an ArcGIS 10 .NET Server instance. Is that pretty much the long and the short of it? I am sure we can't have two different flavour instances of ArcGIS 10 Server installed on the same machine (one Java, one .NET) so they need to be on separate boxes (or VM's)? Thanks in advance for any words of calm...


Test Attribute with 1 or more possible values in FME?

This question continues story on link What's the syntax for an AttributeFilter in FME?. Now, how can I test attribute value if there are more than one value and proceed it to the Clipper (uses as a multiple clippers)? Use some looping method? In xml file is for example   So there can be maybe only one or can be many of them. I tried and Tester output is only one value (last of possible values, in my example above Tester output is only areaID=1253). Tnx,


How do I simplify the straight sections of (multi)Linestring in PostGres/PostGis leaving the curved sections untouched

I have (multi)linestrings containing upwards of 5000 points and I want to perform two simplifications possibly using Ramer Douglas Peucker or alternative methods in PostGis The (multi)linestrings consist of straight lines, 90 degree curves, 180 degree curves and loops (circles) I understand that ST_simplify may remove the loops and possibly some of the bends depending on the tolerance set. I need to have a Zoomed Out version which can probably be done satisfactorily done with ST_simplify However, on the zoomed in version, it is important to keep the loops, and curve points. I would like to remove all the points from the straight segments (within a tolerance) which constitute 90% of the (multi)linestring and keep the points on the curve section. How can I achieve this in PostGis? Thanks Mark


Why are map scales not used instead of Zoom-levels in web mapping clients such as OSM or Google Maps?

Google Maps and OpenstreetMap and many other web mapping products utilize the concept of zoom-levels to change the size of the map in a web client. Most people are however more familiar with the concept of map scales(i.e 1:1000, 1:50000 etc) when zooming in or out of a map or considering measurements. When using zoom-levels the scales are not very intuitive (http://api.geoext.org/1.1/examples/zoom-chooser.html). Even desktop GIS software maintains the concept of scales when zooming in or out of digital maps. So why is the concept of map scale not used instead of Zoom-levels in web mapping clients such as OSM or Google Maps? This would enable the maps to be much more functional for administrative purposes for example in land use planing where the use of set scales like 1:10,000 are required when printing to paper.


Fading symbology at edges

A general problem for me with digitized material is that it can be too precise. Sometimes I just need to show a line or surface without sharp borders. I don't regularly need to visualize the exact extent of, say, a wetland to meter or dm level, only show that it exists in approximately that area. Now I have lines (that I also have buffered 50 m) which I need to "smooth" at the edges. I want the center of the line/buffer to be strong in color. Further away from the center it should be more transparent, and eventually fade away. How can I accomplish this? (I seem to remember seeing a question related to this here a while back, but now I can't find it.)


How to merge street layers in QGIS for exporting to PostGIS?

Is there any way of merging lots of street data in QGIS and getting one layer for exporting it to Postgresql? I know one method using mmqgis "Merge Layer" tool. It works in Ubuntu but not in Windows. When I want to open my MapInfo .TAB file in Windows it opens but this time mmgqis tools do not work. In Ubuntu it gives " _*_.TAB is not a valid or recognized data source "... Can I solve this with some Python code or any other way? Thanks


arcpy: Compare extents of rasters

I am writing a script that takes an input of an Area of Interest (AOI) and then walks through the base directories to find any file that has an extent that intersects with the AOI, clip it and generate a text output file. I've got most of it working except comparing the raster extent to the AOI extent. It is not giving errors but just doesn't work... Any advise? I am also having issues with trying to clip the extent to a polygon, it clips but maintains the MAX/MIN and doesn't clip to the mask. I have tried three methods as in the comments in the script. Extract by Mask says I don't have Spatial Analyst but it is on. Thanks guys, ![enter image description here](http://i.stack.imgur.com/9nJyO.png) ERROR MESSAGE when ExtractByMask is used. > > > > Spatial Analyst license is AVAILABLE > P:\2011\Job_154_PythonScript_for_AOI\Working\Orthophotomosaic\310000_8105000.tif > 310000 8104000 311000 8105000 NaN NaN NaN NaN outside extent Traceback (most > recent call last): File > "P:\2011\Job_154_PythonScript_for_AOI\Working\SubsetGenerator.py", line 80, > in outExtractByMask = ExtractByMask(File, AOI) File "C:\Program > Files\ArcGIS\Desktop10.0\arcpy\arcpy\sa\Functions.py", line 6726, in > ExtractByMask in_mask_data) File "C:\Program > Files\ArcGIS\Desktop10.0\arcpy\arcpy\sa\Utils.py", line 47, in swapper > result = wrapper(*args, **kwargs) File "C:\Program > Files\ArcGIS\Desktop10.0\arcpy\arcpy\sa\Functions.py", line 6722, in wrapper > out_raster) File "C:\Program > Files\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing_base.py", line 474, in > return lambda *args: val(*gp_fixargs(args)) ExecuteError: Failed to execute. > Parameters are not valid. ERROR 000824: The tool is not licensed. Failed to > execute (ExtractByMask). UPDATED/WORKING CODE               # Subset Generator     # Author: George Corea, Atherton Tablelands GIS     # georgec@atgis.com.au; info@atgis.com.au     # Licence:     #AddMessage, AddWarning, AddError          import arcinfo     import arcpy, glob, os, sys, string     from arcpy import env     from arcpy.sa import *          path = os.getcwd()     os.chdir(path)          class LicenseError(Exception):         pass          try:         if arcpy.CheckExtension("Spatial") == "Available":             arcpy.CheckOutExtension("Spatial")             print "Spatial Analyst license is AVAILABLE"         else:             # raise a custom exception             #             raise LicenseError              env.workspace = path+'\\junk'               except LicenseError:         print "Spatial Analyst license is unavailable"     except:         print arcpy.GetMessages(2)               RootDirectories = [r'P:\2011\Job_154_PythonScript_for_AOI\Working\Orthophotomosaic']     RasterTypes = ['tif','jpg', 'ecw'] #boolean, valuelist               #AOI = arcpy.mapping.Layer('AOI.shp')     #AOIextent=AOI.getExtent()     AOI = 'AOI.tif'     AOIobj=arcpy.Raster(AOI)     AOIextent=AOIobj.extent          #print AOIextent     SR = 0     x=0     Subset=0     FileArea=0     f = open(path+'\\AOI_Clip\\AOI_SubsetGenerator.txt', 'a')     f.write("UniqueID, baseName, extension, dataType, path, \     OrigName, NewName, size(uncompressed bytes), area(m2)"+"\n")     f.close()          for RootDirectory in RootDirectories:         for root, dirs, files in os.walk(RootDirectory):             for RasterType in RasterTypes:                 FileList = [os.path.join(root, fi) for fi in files \                             if fi.endswith(RasterType)]                 for File in FileList:                     FileDesc=arcpy.Describe(File)                     SR = FileDesc.spatialReference                     if SR.name != "Unknown":                         FileObj = arcpy.Raster(File)                         FileExtent = FileObj.extent                         print "Processing: " + File                         Extent_Overlap = str(FileExtent.overlaps(AOIextent))                         Extent_Within = str(FileExtent.within(AOIextent))                         Extent_Touches = str(FileExtent.touches(AOIextent))                         Extent_Crosses = str(FileExtent.crosses(AOIextent))                         if Extent_Overlap == 'True':                             Subset =1                             print "Extent_Overlaps"                         elif Extent_Within == 'True':                             Subset =1                             print "Extent_Within"                         elif Extent_Touches == 'True':                             Subset =1                             print "Extent_Touches"                         elif Extent_Crosses == 'True':                             Subset =1                             print "Extent_Crosses"                         else:                             #print "OUTSIDE extent"+str(Subset)                             Subset =0          #Process files that are subset                              if Subset == 1:                             #print Subset                             x=x+1                             outRaster=path+'\\AOI_Clip\\'+str(FileDesc.baseName)+"^AOI_Clip"+str(x)+"."+str(FileDesc.extension)                                  try:                                 outExtractByMask = ExtractByMask(File, 'AOI.shp')                                 outExtractByMask.save(outRaster)                                 print "Created: " + str(outRaster)                             except:                                 print "Raster is outside of extent"                                 #FileArea=float(arcpy.GetRasterProperties_management(File, "ROWCOUNT")*arcpy.GetRasterProperties_management(File, "COLUMNCOUNT")*arcpy.GetRasterProperties_management(File, "CELLSIZEY")*2)                             FileArea=\                                   (string.atol(str(arcpy.GetRasterProperties_management(outRaster, "ROWCOUNT")))\                                   *string.atol(str(arcpy.GetRasterProperties_management(outRaster, "COLUMNCOUNT"))))\                                   *(string.atol((str(arcpy.GetRasterProperties_management(outRaster, "CELLSIZEY")))\                                   *string.atol((str(arcpy.GetRasterProperties_management(outRaster, "CELLSIZEX"))))))                             print FileArea                             f = open(path+'\\AOI_Clip\\AOI_SubsetGenerator.txt', 'a')                             f.write(str(x)+","+str(FileDesc.baseName)+","\                                     +str(FileDesc.extension)+","\                                     +str(FileDesc.dataType)+","\                                     +str(FileDesc.path)+","\                                     +str(SR.name)+","\                                     +str(outRaster)+","\                                     +str(FileObj.uncompressedSize)+","\                                     +str(FileArea)                                     +"\n")                             f.close()                         else:                                 print "No Subsets"                  print "Changing Directory to: " + str(dirs)          arcpy.CheckInExtension("Spatial")     print "Checked Out extensions and Exiting"


How to check whether a layer is vector or not in OpenLayers?

am not sure whether this is a stupid question or not.. but I am wondering how to dynamically check whether a layer object belongs to OpenLayers.Layer.Vector class. Thanks for any replies!


How to solve an installation problem involving QGIS, Python, and a plug-in?

In an attempt to get the plugin "Profile Tool" working, I installed geoslib and python Shapely 1.2.10 module as requested using following:               sudo port install geos      Then install Shapely once downloaded:               python setup.py install      Now whenever I load a Python plugin in QGIS, the whole app crashes. Loading with "--noplugins" or deleting plugins from ~/.qgis is the only way to get QGIS to load. I suspect I have messed up the Python installation. Which version of Python should QGIS being using? The version in /usr/bin/python or does it come with its own version? From within QGIS, the python console displays path as:               >>> import sys      > > > sys.path > ['/Applications/GIS/QGIS.app/Contents/MacOS/../Resources/python', > '/Users/jon/.qgis/python', '/Users/jon/.qgis/python/plugins', > '/Applications/GIS/QGIS.app/Contents/MacOS/../Resources/python/plugins', > '/Library/Frameworks/SQLite3.framework/Versions/3/Python', > '/Library/Frameworks/GDAL.framework/Versions/1.8/Python/site-packages', > '/Library/Python/2.6/site-packages/Shapely-1.2.10-py2.6.egg', > '/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python26.zip', > '/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6', > '/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/plat- > darwin', > '/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/plat- > mac', > '/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/plat- > mac/lib-scriptpackages', > '/System/Library/Frameworks/Python.framework/Versions/2.6/Extras/lib/python', > '/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib- > tk', > '/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib- > old', > '/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib- > dynload', '/Library/Python/2.6/site-packages', > '/System/Library/Frameworks/Python.framework/Versions/2.6/Extras/lib/python/PyObjC', > '/System/Library/Frameworks/Python.framework/Versions/2.6/Extras/lib/python/wx-2.8-mac- > unicode', > '/Applications/GIS/QGIS.app/Contents/Resources/python/plugins/fTools/tools'] Crash log below. Any help greatly appreciated. I am stuck with a non-functioning QGIS app and can't do much at the minute without some plugins.                   Process:         QGIS [1320]     Path:            /Applications/GIS/QGIS.app/Contents/MacOS/QGIS     Identifier:      org.qgis.qgis     Version:         1.7.0 (1.7.0-Wroclaw (exported))     Code Type:       X86-64 (Native)     Parent Process:  launchd [227]          Date/Time:       2011-08-15 17:12:28.401 +1000     OS Version:      Mac OS X 10.6.8 (10K540)     Report Version:  6     Sleep/Wake UUID: 6C49B351-8735-4919-8292-F2C1CEE58FFB          Exception Type:  EXC_BAD_ACCESS (SIGSEGV)     Exception Codes: KERN_INVALID_ADDRESS at 0x0000000000000000     Crashed Thread:  0  Dispatch queue: com.apple.main-thread          Thread 0 Crashed:  Dispatch queue: com.apple.main-thread     0   ???                              000000000000000000 0 + 0     1   org.python.python                0x000000011bb9fd3b _PyImport_LoadDynamicModule + 189     2   org.python.python                0x000000011bb9eaf8 PyImport_ReloadModule + 1347     3   org.python.python                0x000000011bb9ef78 PyImport_ReloadModule + 2499     4   org.python.python                0x000000011bb9f59a PyImport_ImportModuleLevel + 1265     5   org.python.python                0x000000011bb879cb _PyBuiltin_Init + 14264     6   org.python.python                0x000000011bb8d187 PyEval_EvalFrameEx + 15317     7   org.python.python                0x000000011bb8ecce PyEval_EvalCodeEx + 1803     8   org.python.python                0x000000011bb308e1 PyClassMethod_New + 1748     9   org.python.python                0x000000011bb0eff3 PyObject_Call + 112     10  org.python.python                0x000000011bb889db PyEval_CallObjectWithKeywords + 175     11  org.python.python                0x000000011bb8c5ea PyEval_EvalFrameEx + 12344     12  org.python.python                0x000000011bb8ecce PyEval_EvalCodeEx + 1803     13  org.python.python                0x000000011bb8ed61 PyEval_EvalCode + 54     14  org.python.python                0x000000011bb9c8f2 PyImport_ExecCodeModuleEx + 258     15  org.python.python                0x000000011bb9d248 PyImport_AppendInittab + 1026     16  org.python.python                0x000000011bb9e1d0 _PyImport_FindModule + 379     17  org.python.python                0x000000011bb9eaf8 PyImport_ReloadModule + 1347     18  org.python.python                0x000000011bb9ef78 PyImport_ReloadModule + 2499     19  org.python.python                0x000000011bb9f59a PyImport_ImportModuleLevel + 1265     20  org.python.python                0x000000011bb879cb _PyBuiltin_Init + 14264     21  org.python.python                0x000000011bb8d187 PyEval_EvalFrameEx + 15317     22  org.python.python                0x000000011bb8ecce PyEval_EvalCodeEx + 1803     23  org.python.python                0x000000011bb308e1 PyClassMethod_New + 1748     24  org.python.python                0x000000011bb0eff3 PyObject_Call + 112     25  org.python.python                0x000000011bb889db PyEval_CallObjectWithKeywords + 175     26  org.python.python                0x000000011bb8c5ea PyEval_EvalFrameEx + 12344     27  org.python.python                0x000000011bb8ecce PyEval_EvalCodeEx + 1803     28  org.python.python                0x000000011bb8ed61 PyEval_EvalCode + 54     29  org.python.python                0x000000011bb9c8f2 PyImport_ExecCodeModuleEx + 258     30  org.python.python                0x000000011bb9d248 PyImport_AppendInittab + 1026     31  org.python.python                0x000000011bb9eaf8 PyImport_ReloadModule + 1347     32  org.python.python                0x000000011bb9ef78 PyImport_ReloadModule + 2499     33  org.python.python                0x000000011bb9f56a PyImport_ImportModuleLevel + 1217     34  org.python.python                0x000000011bb879cb _PyBuiltin_Init + 14264     35  org.python.python                0x000000011bb8d187 PyEval_EvalFrameEx + 15317     36  org.python.python                0x000000011bb8ecce PyEval_EvalCodeEx + 1803     37  org.python.python                0x000000011bb308e1 PyClassMethod_New + 1748     38  org.python.python                0x000000011bb0eff3 PyObject_Call + 112     39  org.python.python                0x000000011bb889db PyEval_CallObjectWithKeywords + 175     40  org.python.python                0x000000011bb8c5ea PyEval_EvalFrameEx + 12344     41  org.python.python                0x000000011bb8ecce PyEval_EvalCodeEx + 1803     42  org.python.python                0x000000011bb8ed61 PyEval_EvalCode + 54     43  org.python.python                0x000000011bb9c8f2 PyImport_ExecCodeModuleEx + 258     44  org.python.python                0x000000011bb9d248 PyImport_AppendInittab + 1026     45  org.python.python                0x000000011bb9eaf8 PyImport_ReloadModule + 1347     46  org.python.python                0x000000011bb9ef78 PyImport_ReloadModule + 2499     47  org.python.python                0x000000011bb9f56a PyImport_ImportModuleLevel + 1217     48  org.python.python                0x000000011bb879cb _PyBuiltin_Init + 14264     49  org.python.python                0x000000011bb8d187 PyEval_EvalFrameEx + 15317     50  org.python.python                0x000000011bb8ecce PyEval_EvalCodeEx + 1803     51  org.python.python                0x000000011bb308e1 PyClassMethod_New + 1748     52  org.python.python                0x000000011bb0eff3 PyObject_Call + 112     53  org.python.python                0x000000011bb889db PyEval_CallObjectWithKeywords + 175     54  org.python.python                0x000000011bb8c5ea PyEval_EvalFrameEx + 12344     55  org.python.python                0x000000011bb8d2e1 PyEval_EvalFrameEx + 15663     56  org.python.python                0x000000011bb8ecce PyEval_EvalCodeEx + 1803     57  org.python.python                0x000000011bb8ed61 PyEval_EvalCode + 54     58  org.python.python                0x000000011bba665a Py_CompileString + 78     59  org.python.python                0x000000011bba6805 PyRun_StringFlags + 119     60  org.python.python                0x000000011bb8b4ca PyEval_EvalFrameEx + 7960     61  org.python.python                0x000000011bb8ecce PyEval_EvalCodeEx + 1803     62  org.python.python                0x000000011bb8d35e PyEval_EvalFrameEx + 15788     63  org.python.python                0x000000011bb8d2e1 PyEval_EvalFrameEx + 15663     64  org.python.python                0x000000011bb8d2e1 PyEval_EvalFrameEx + 15663     65  org.python.python                0x000000011bb8d2e1 PyEval_EvalFrameEx + 15663     66  org.python.python                0x000000011bb8ecce PyEval_EvalCodeEx + 1803     67  org.python.python                0x000000011bb8ed61 PyEval_EvalCode + 54     68  org.python.python                0x000000011bba665a Py_CompileString + 78     69  org.python.python                0x000000011bba6805 PyRun_StringFlags + 119     70  libqgispython.1.7.0.dylib        0x000000011bafe496 QgsPythonUtilsImpl::evalString(QString const&, QString&) + 122     71  libqgispython.1.7.0.dylib        0x000000011bafdd21 QgsPythonUtilsImpl::startPlugin(QString) + 111     72  org.qgis.qgis                    0x0000000100110dc6 QgsPluginRegistry::loadPythonPlugin(QString) + 496     73  org.qgis.qgis                    0x00000001001129f7 QgsPluginRegistry::restoreSessionPlugins(QString) + 2193     74  org.qgis.qgis                    0x000000010002ae6e QgisApp::QgisApp(QSplashScreen*, bool, QWidget*, QFlags<Qt::WindowType>) + 3050     75  org.qgis.qgis                    0x0000000100012057 main + 7683     76  org.qgis.qgis                    0x000000010000feac start + 52          Thread 1:  Dispatch queue: com.apple.libdispatch-manager     0   libSystem.B.dylib                0x00007fff8585ac0a kevent + 10     1   libSystem.B.dylib                0x00007fff8585cadd _dispatch_mgr_invoke + 154     2   libSystem.B.dylib                0x00007fff8585c7b4 _dispatch_queue_invoke + 185     3   libSystem.B.dylib                0x00007fff8585c2de _dispatch_worker_thread2 + 252     4   libSystem.B.dylib                0x00007fff8585bc08 _pthread_wqthread + 353     5   libSystem.B.dylib                0x00007fff8585baa5 start_wqthread + 13          Thread 2:     0   libSystem.B.dylib                0x00007fff8585ba2a __workq_kernreturn + 10     1   libSystem.B.dylib                0x00007fff8585be3c _pthread_wqthread + 917     2   libSystem.B.dylib                0x00007fff8585baa5 start_wqthread + 13          Thread 3:     0   libSystem.B.dylib                0x00007fff8585ba2a __workq_kernreturn + 10     1   libSystem.B.dylib                0x00007fff8585be3c _pthread_wqthread + 917     2   libSystem.B.dylib                0x00007fff8585baa5 start_wqthread + 13          Thread 4:     0   libSystem.B.dylib                0x00007fff8585ba2a __workq_kernreturn + 10     1   libSystem.B.dylib                0x00007fff8585be3c _pthread_wqthread + 917     2   libSystem.B.dylib                0x00007fff8585baa5 start_wqthread + 13          Thread 0 crashed with X86 Thread State (64-bit):       rax: 0x000000011fd5b840  rbx: 0x000000011f740408  rcx: 0x0000000000000b31  rdx: 0x000000011fd5b840       rdi: 0x000000011bc9ed00  rsi: 0x000000011f740730  rbp: 0x00007fff5fbf80b0  rsp: 0x00007fff5fbf8098        r8: 0x0000000000000000   r9: 0x000000011cb84ad0  r10: 0x000000011cb7afe0  r11: 0x000000011bb3fe7d       r12: 0x000000011bc9ed00  r13: 0x0000000000000000  r14: 0x00007fff5fbf8630  r15: 0x00007fff5fbf8150       rip: 0x0000000000000000  rfl: 0x0000000000010246  cr2: 0x0000000000000000          Binary Images:            0x100000000 -        0x100491fff +org.qgis.qgis 1.7.0 (1.7.0-Wroclaw (exported)) <1D91FB61-31DA-06DE-DA64-5AB29A3C5F07> /Applications/GIS/QGIS.app/Contents/MacOS/QGIS            0x1005ee000 -        0x10067bff7 +libqwt.dylib 5.2.1 (compatibility 5.2.0) <A4BDC10C-06CC-D7B8-BEB8-0037CEDEA165> /Applications/GIS/QGIS.app/Contents/MacOS/lib/libqwt.dylib            0x1006cc000 -        0x100703fff +QtSql 4.7.3 (compatibility 4.7.0) <F5B9EBFD-65ED-35D6-26E3-769F1E0DEF2B> /Applications/GIS/QGIS.app/Contents/Frameworks/QtSql.framework/Versions/4/QtSql            0x10071a000 -        0x10096cfe7 +libqgis_core.1.7.0.dylib 1.7.0 (compatibility 1.7.0) <9FC26325-F285-0791-0D3E-2F4AE77DC454> /Applications/GIS/QGIS.app/Contents/MacOS/lib/libqgis_core.1.7.0.dylib            0x100a5f000 -        0x100bc9ff7 +libqgis_gui.1.7.0.dylib 1.7.0 (compatibility 1.7.0) <15EC96DB-4381-4023-CA0F-EEAD93420BA6> /Applications/GIS/QGIS.app/Contents/MacOS/lib/libqgis_gui.1.7.0.dylib            0x100c7e000 -        0x100cbffef +libqgis_analysis.1.7.0.dylib 1.7.0 (compatibility 1.7.0) <85BFF9C9-A84F-D989-75E3-C803CA5B5C37> /Applications/GIS/QGIS.app/Contents/MacOS/lib/libqgis_analysis.1.7.0.dylib            0x100cd9000 -        0x100cf5fff +libpq.dylib 5.3.0 (compatibility 5.0.0) <1FB622DB-BAB7-F4CD-D9A9-7E024478EBB7> /Applications/GIS/QGIS.app/Contents/MacOS/lib/libpq.dylib            0x100cff000 -        0x10105cfff +org.sqlite.sqlite3 3.7.4 (SQLite3 3.7.4-2) <C90C1409-CAD8-8662-6047-667AA6108EB8> /Library/Frameworks/SQLite3.framework/Versions/3/SQLite3            0x10107b000 -        0x1010b8fff +QtXml 4.7.3 (compatibility 4.7.0) <811907C3-6FE0-EA37-67F0-8BE955F0BA1B> /Applications/GIS/QGIS.app/Contents/Frameworks/QtXml.framework/Versions/4/QtXml            0x1010d2000 -        0x10135aff7 +QtCore 4.7.3 (compatibility 4.7.0) <1A2CF13B-BB74-85EA-5A4D-1C477E02DF6E> /Applications/GIS/QGIS.app/Contents/Frameworks/QtCore.framework/Versions/4/QtCore            0x1013ee000 -        0x101dcbff7 +QtGui 4.7.3 (compatibility 4.7.0) <1F568542-C5A4-8A3C-5006-EA65FEDE23F6> /Applications/GIS/QGIS.app/Contents/Frameworks/QtGui.framework/Versions/4/QtGui            0x102075000 -        0x10219cff7 +QtNetwork 4.7.3 (compatibility 4.7.0) <691D3311-87BE-C23E-49F3-FC300D428C61> /Applications/GIS/QGIS.app/Contents/Frameworks/QtNetwork.framework/Versions/4/QtNetwork            0x1021f5000 -        0x102244fef +QtSvg 4.7.3 (compatibility 4.7.0) <A7789926-A9A6-9BE3-4C11-E8D58E816DFC> /Applications/GIS/QGIS.app/Contents/Frameworks/QtSvg.framework/Versions/4/QtSvg            0x102265000 -        0x102294fff +org.maptools.proj 4.7.0 (PROJ 4.7.0-2) <7FCE23B7-503C-8C93-6BD2-5D08C5B262DB> /Library/Frameworks/PROJ.framework/Versions/4/PROJ            0x1022a0000 -        0x1023abfff +org.osgeo.geos 3.3.0 (GEOS 3.3.0-1) <BF9394BA-B8CD-4150-250A-98A38060BEB8> /Library/Frameworks/GEOS.framework/Versions/3/GEOS            0x102506000 -        0x103265fe7 +org.gdal.gdal 1.8.0 (GDAL/OGR 1.8.0-2) <92C8A957-8799-F90E-FC59-05DE1B59FCA6> /Library/Frameworks/GDAL.framework/Versions/1.8/GDAL            0x10388a000 -        0x1038cdff7  libiodbc.2.dylib 4.18.0 (compatibility 4.0.0) <9E10F33D-B573-3B98-E756-392037B9229E> /usr/lib/libiodbc.2.dylib            0x1038d8000 -        0x1038e6fff  libiodbcinst.2.dylib 4.18.0 (compatibility 4.0.0) <32500035-5777-D45C-E5AD-54A5D0C4F461> /usr/lib/libiodbcinst.2.dylib            0x1038ec000 -        0x1039f9fef +com.kyngchaos.UnixImageIO 1.3.0 (UnixImageIO 1.3.0) <9F2021D8-AA8A-DD27-0BFB-96335704C28E> /Library/Frameworks/UnixImageIO.framework/Versions/D/UnixImageIO            0x103a31000 -        0x103a69fe7 +libogdi.dylib ??? (???) <F2BF367E-21A2-CEC2-1B5B-64E0508DBB1F> /Library/Frameworks/GDAL.framework/Versions/1.8/Libraries/libogdi.dylib            0x103a75000 -        0x103a7dfef +libintl.dylib 10.1.0 (compatibility 10.0.0) <937F7B46-EC9C-4441-43FD-4CD508BE99E2> /Applications/GIS/QGIS.app/Contents/MacOS/lib/libintl.dylib            0x103cc3000 -        0x103cd0fff +com.rogueamoeba.audio_hijack_server.hermes 2.2.8 (2.2.8) <DBE64033-014D-0508-024D-EE1BDA4C8528> /usr/local/hermes/modules/Instant Hijack Server.hermesmodule/Contents/MacOS/Instant Hijack Server            0x103ffa000 -        0x103ffbff7  _functools.so ??? (???) <97A294C7-290C-9746-B294-275516490F40> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/_functools.so            0x11675b000 -        0x116796ff7 +libqjpeg.dylib ??? (???) <2D99F9D4-B313-6024-82E7-42B9EB233A47> /Applications/GIS/QGIS.app/Contents/PlugIns/imageformats/libqjpeg.dylib            0x11679f000 -        0x1167acfff +libcoordinatecaptureplugin.so ??? (???) <38B5DF22-D8D7-16FF-44C1-0DA3BA3E0BF7> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libcoordinatecaptureplugin.so            0x1167b5000 -        0x1167c1fff +libcopyrightlabelplugin.so ??? (???) <659E091F-9BCD-47DC-7A04-1FA7B197D9DC> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libcopyrightlabelplugin.so            0x1167fa000 -        0x1167fbfff  cStringIO.so ??? (???) <0C7D1D15-4241-99A9-7671-A3CF105C2EBB> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/cStringIO.so            0x116ed6000 -        0x116ee9ff7 +libdelimitedtextplugin.so ??? (???) <4C8C4C33-BDB6-663A-381B-4AE3B9F766B0> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libdelimitedtextplugin.so            0x116ef3000 -        0x116ef5ff7  math.so ??? (???) <54D066FD-A9F9-EA69-8C3C-CE83DF68C58A> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/math.so            0x116efb000 -        0x116efcfff  _random.so ??? (???) <FC54824D-3F1F-D538-DBF8-05E393A696B6> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/_random.so            0x117f80000 -        0x117f8cfef +libdelimitedtextprovider.so ??? (???) <7344FDA8-1F61-8093-6087-B970C1B3601A> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libdelimitedtextprovider.so            0x117f94000 -        0x117fb7ff7 +libdiagramoverlay.so ??? (???) <0C636064-256B-F93C-9D28-86EC28330DFC> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libdiagramoverlay.so            0x117fcf000 -        0x117fe3ff7 +libdisplacementplugin.so ??? (???) <003F036F-251B-DB07-E119-4AEA19C2075B> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libdisplacementplugin.so            0x117ff2000 -        0x117ff5fff +_gdalconst.so ??? (???) <20E606D9-4BB8-8821-CC78-C31B17CC804C> /Library/Frameworks/GDAL.framework/Versions/1.8/Python/site-packages/osgeo/_gdalconst.so            0x117ff9000 -        0x117ffbfe7  binascii.so ??? (???) <D7D60C49-DEFA-5DB1-83B0-F91347966BC8> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/binascii.so            0x119461000 -        0x119482fef +libdxf2shpconverterplugin.so ??? (???) <E210C964-F20A-0C79-8E18-BD4A2EB79A99> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libdxf2shpconverterplugin.so            0x119494000 -        0x1194c9fef +libevis.so ??? (???) <F818705D-59FA-A057-6B4A-4E2C8101D81D> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libevis.so            0x1194e0000 -        0x1194f1ff7 +libgdalprovider.so ??? (???) <E4501494-A2FD-D14A-C63E-B5515D10EE6F> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libgdalprovider.so            0x1194fb000 -        0x11956bfff +libgeorefplugin.so ??? (???) <5183AEB7-CBA0-AC52-104F-DF4127CD45F4> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libgeorefplugin.so            0x11959d000 -        0x11971cfef +org.gnu.gsl 1.14 (GSL 1.14-1) <7A5D6DBB-2673-BE9B-EA2A-5FA502828B83> /Library/Frameworks/GSL.framework/Versions/1/GSL            0x119775000 -        0x119778fff  _csv.so ??? (???) <ACE39A2B-E3BC-EC0C-D825-1812240F2B06> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/_csv.so            0x11b701000 -        0x11b730fe7 +libgpsimporterplugin.so ??? (???) <3B10A0BE-C0F2-1B82-40A5-A945F309E882> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libgpsimporterplugin.so            0x11b742000 -        0x11b758fe7 +libgpxprovider.so ??? (???) <7D83BB8A-090F-92BA-AD8F-5727FA698EAA> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libgpxprovider.so            0x11b767000 -        0x11b77bff7 +libinterpolationplugin.so ??? (???) <519BF936-025C-4D4B-BA9F-1AFC827A3C5B> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libinterpolationplugin.so            0x11b789000 -        0x11b795fff +libmemoryprovider.so ??? (???) <1ADC4FEA-C6C9-4F36-EDFA-72123FC375E7> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libmemoryprovider.so            0x11b79e000 -        0x11b7abfe7 +libnortharrowplugin.so ??? (???) <83F8A703-C231-D010-FEA4-A538337BA922> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libnortharrowplugin.so            0x11b7b5000 -        0x11b7d3fe7 +libofflineeditingplugin.so ??? (???) <48F67327-1A7C-3268-98E9-C8B3BAE3EDE9> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libofflineeditingplugin.so            0x11b7e5000 -        0x11b800ff7 +libogrprovider.so ??? (???) <444694DF-5A9C-42D5-E0EC-8E05046F16C7> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libogrprovider.so            0x11b80d000 -        0x11b81efef +liboracleplugin.so ??? (???) <5C7297D3-EBD5-1573-DFDD-3C81DB9A2782> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/liboracleplugin.so            0x11b828000 -        0x11b847fef +libosmprovider.so ??? (???) <87F44AF8-B191-C2F5-BA78-0F11584787EB> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libosmprovider.so            0x11b857000 -        0x11b894fe7 +libpostgresprovider.so ??? (???) <D6F01231-8B79-963F-9E8D-0E58BCF9178A> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libpostgresprovider.so            0x11b8a5000 -        0x11b8b1fef +librasterterrainplugin.so ??? (???) <5F8F4162-E5B0-DA3B-7747-DFBEBABFD1B3> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/librasterterrainplugin.so            0x11b8ba000 -        0x11b8dcff7 +libroadgraphplugin.so ??? (???) <8852024C-6520-BBDC-DC11-19DA4EDB56C0> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libroadgraphplugin.so            0x11b8f3000 -        0x11b901fff +libscalebarplugin.so ??? (???) <C2070C09-CD50-B394-F7BF-E940C76A29A3> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libscalebarplugin.so            0x11b90b000 -        0x11b92bfef +libspatialiteprovider.so ??? (???) <7DCD0C0C-EFC4-89E4-C027-60AD5D24CBC6> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libspatialiteprovider.so            0x11b936000 -        0x11b953fff +libspatialqueryplugin.so ??? (???) <10E21ED9-FEEC-9E88-28D0-E3072EEE4E2E> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libspatialqueryplugin.so            0x11b964000 -        0x11b9b3fef +libspitplugin.so ??? (???) <D729BCBF-BE8A-7A18-C3DB-B5A7D8C1689D> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libspitplugin.so            0x11b9c6000 -        0x11b9f4fe7 +libsqlanywhereplugin.so ??? (???) <12FC3AA4-E9EE-7E27-6825-1BBAB532ED6C> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libsqlanywhereplugin.so            0x11ba0c000 -        0x11ba1aff7 +libqgissqlanyconnection.1.7.0.dylib 1.7.0 (compatibility 1.7.0) <49430729-C774-59EE-AC9F-AE6A2EC92C1A> /Applications/GIS/QGIS.app/Contents/MacOS/lib/libqgissqlanyconnection.1.7.0.dylib            0x11ba20000 -        0x11ba3efef +libsqlanywhereprovider.so ??? (???) <F74D1F73-A392-74BF-046E-D68CB47C482C> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libsqlanywhereprovider.so            0x11ba49000 -        0x11ba6dfef +libwfsplugin.so ??? (???) <2E749ABF-A44F-D28F-92F5-A74EE81C74B2> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libwfsplugin.so            0x11ba7e000 -        0x11baa7fff +libwfsprovider.so ??? (???) <1FD25BE5-6AA6-6FEB-144A-3AB15D99ACBD> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libwfsprovider.so            0x11bab9000 -        0x11bae8fe7 +libwmsprovider.so ??? (???) <E98A3441-B816-DC87-0CCB-39F74406627A> /Applications/GIS/QGIS.app/Contents/PlugIns/qgis/libwmsprovider.so            0x11baf8000 -        0x11bb03ff7 +libqgispython.1.7.0.dylib 1.7.0 (compatibility 1.7.0) <F590EC1D-0215-D020-0BBA-00EE88F142C7> /Applications/GIS/QGIS.app/Contents/MacOS/lib/libqgispython.1.7.0.dylib            0x11bb08000 -        0x11bc18ff7  org.python.python 2.6.1 (2.6.1) <126DA8FF-5BC2-8788-51E3-D7A29A3F9F0F> /System/Library/Frameworks/Python.framework/Versions/2.6/Python            0x11bd68000 -        0x11bd7bfef +sip.so ??? (???) <A207F0A4-B4B2-A7C3-54FB-CC2A0EE2CAEC> /Applications/GIS/QGIS.app/Contents/Resources/python/sip.so            0x11bd84000 -        0x11bf25fe7 +QtCore.so ??? (???) <BFB7B36F-C825-98EE-3ECE-DD72A17DF992> /Applications/GIS/QGIS.app/Contents/Resources/python/PyQt4/QtCore.so            0x11c0f9000 -        0x11c0faff7  time.so ??? (???) <80513398-F49E-79D1-5014-514361869D40> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/time.so            0x11c200000 -        0x11c719ff7 +QtGui.so ??? (???) <9DD97B5D-46EA-89A8-F4B2-C7E0C0E24FF1> /Applications/GIS/QGIS.app/Contents/Resources/python/PyQt4/QtGui.so            0x11cae1000 -        0x11cae6fff  itertools.so ??? (???) <9287854F-7F2B-D4AF-FCA3-EB69DA821DA9> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/itertools.so            0x11caee000 -        0x11caf1ff7  strop.so ??? (???) <4A91CDB0-6E91-DA0E-8E6B-38BE29105EA0> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/strop.so            0x11caf6000 -        0x11cafaff7  _struct.so ??? (???) <B4DD710D-1511-99F0-66AE-329185E2B88C> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/_struct.so            0x11cc00000 -        0x11cd81ff7 +core.so ??? (???) <50F0DE0C-1095-D401-4D07-A4D41D35F0E8> /Applications/GIS/QGIS.app/Contents/Resources/python/qgis/core.so            0x11ce87000 -        0x11cef9fef +QtNetwork.so ??? (???) <B7A8946F-0C04-7A26-986F-771060CFCF25> /Applications/GIS/QGIS.app/Contents/Resources/python/PyQt4/QtNetwork.so            0x11cf51000 -        0x11cf87fff +QtXml.so ??? (???) <0E038330-445A-CFCD-9999-ECF0792C36BD> /Applications/GIS/QGIS.app/Contents/Resources/python/PyQt4/QtXml.so            0x11cff1000 -        0x11cff4ff7  zlib.so ??? (???) <647721E3-67B5-8CD0-3A78-060FF9C80924> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/zlib.so            0x11cff9000 -        0x11cffafff +Qt.so ??? (???) <7FA143FB-1433-4AB0-7E5B-05A5A912637F> /Applications/GIS/QGIS.app/Contents/Resources/python/PyQt4/Qt.so            0x11d100000 -        0x11d1e0fff +gui.so ??? (???) <E6263E74-A069-A45D-502D-D7AB4D2EE335> /Applications/GIS/QGIS.app/Contents/Resources/python/qgis/gui.so            0x11d27d000 -        0x11d2c5fe7 +QtWebKit.so ??? (???) <3068F58D-058F-4AA8-28B1-A1ABF5323755> /Applications/GIS/QGIS.app/Contents/Resources/python/PyQt4/QtWebKit.so            0x11d2ff000 -        0x11e5b4fe7 +QtWebKit 4.7.3 (compatibility 4.7.0) <0722A3F1-7C5D-2DE4-B431-8DF4587AF8B5> /Applications/GIS/QGIS.app/Contents/Frameworks/QtWebKit.framework/Versions/4/QtWebKit            0x11ed4f000 -        0x11ed97fef +phonon 4.4.0 (compatibility 4.4.0) <BB166737-664C-0BA3-684B-55C8E936958A> /Applications/GIS/QGIS.app/Contents/Frameworks/phonon.framework/Versions/4/phonon            0x11edbe000 -        0x11ee3dff7 +QtDBus 4.7.3 (compatibility 4.7.0) <2082C4C0-0173-10A0-512B-1F4ACDD65F2E> /Applications/GIS/QGIS.app/Contents/Frameworks/QtDBus.framework/Versions/4/QtDBus            0x11ee61000 -        0x11eea8ff7 +_gdal.so ??? (???) <326CAECE-6D5A-E600-EA6B-FB05F31C771F> /Library/Frameworks/GDAL.framework/Versions/1.8/Python/site-packages/osgeo/_gdal.so            0x11eefa000 -        0x11efe0ff7 +gdal_ECW_JP2ECW.dylib ??? (???) <3D045D3F-3501-BA5F-2156-7E903351CF15> /Library/Application Support/GDAL/1.8/PlugIns/gdal_ECW_JP2ECW.dylib            0x11f08f000 -        0x11f0daff7 +_ogr.so ??? (???) <B78AB528-01A1-7056-FA76-7D933CDFB486> /Library/Frameworks/GDAL.framework/Versions/1.8/Python/site-packages/osgeo/_ogr.so            0x11f0e9000 -        0x11f0f4fe7  datetime.so ??? (???) <C8DD63B4-9685-C134-E60A-0ED11A0CCAF7> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/datetime.so            0x11f240000 -        0x11f274fff +_osr.so ??? (???) <D4B2287A-3AD5-90FA-8422-D515ECAF0A5B> /Library/Frameworks/GDAL.framework/Versions/1.8/Python/site-packages/osgeo/_osr.so            0x11f5bf000 -        0x11f5c8fff  _sqlite3.so ??? (???) <9FBFA469-EEA4-67CA-B92D-5C30D91887E9> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/_sqlite3.so            0x11f5d2000 -        0x11f64ffef  unicodedata.so ??? (???) <27EF63BF-90E9-3D71-FE24-1EA634E99498> /System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/lib-dynload/unicodedata.so            0x11f7e4000 -        0x11f7f1fff  _sort.so ??? (???) <AC4F6F41-7C43-A370-7044-4DF77A34C674> /System/Library/Frameworks/Python.framework/Versions/2.6/Extras/lib/python/numpy/core/_sort.so      etc etc


How to select distinct values from a column in Mapbasic?

I have a table in a mapinfo .tab file. Is there any way to select distinct values from the column? I tried :-               SELECT distinct (some_column) from table_name into temp_table      But this query would select all values. I wanted to select only distinct value from that column. Is there any way to achieve it?


(How) do the modifiable areal unit problem and modifiable temporal unit problem interact?

I am reconstructing historical time series of land use and other related data from a variety of sources with different spatial and temporal resolutions and would be interested in what possible (hidden) problems could arise when aggregating and combining the data. I am especially curious about the ecological fallacy and whether choices in the aggregation of temporal data affect change of support in the spatial dimension and vice versa? References to literature are fine, but I am more interested in examples of how those interactions would arise if at all.


List of free WFS 2.0 clients

I have the task to provide my point observations data in the "WFS 2.0" (Web Feature Service 2.0) standard format. Are there any free, easy-to-use, desktop GIS clients or online viewers for showing and querying my WFS 2.0 layers? I have only found the "QGIS WFS 2.0 Client plugin" Do any other popular free GIS viewers (Kosmo, MapWindow, ArcGIS explorer..) have support for showing and querying WFS 2.0 web feature service?


QGIS opening shp file: why does it crash?

I am trying to open a shp file from the EEA (Urban Atlas), and with some of the shp files qgis keeps crashing as soon as I try to load it. QGIS 1.8.0 is running on Windows 7; I reinstalled QGIS, reloaded the shp file from the EEA website. Metadata of the shp file is: Geometry type of the features in this layer: Polygon The number of features in this layer: 15271 Editing capabilities of this layer: Add Features, Delete Features, Change Attribute Values, Add Attributes, Delete Attributes, Create Spatial Index, Fast Access to Features at ID, Change Geometries Extents: In layer spatial reference system units : xMin,yMin 4208112.15,2360649.92 : xMax,yMax 4271746.52,2392865.69 Layer Spatial Reference System: +proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs


How to Properly Share as Image Service?

I have 5 historical aerial sets for my county each consisting of several hundred SID files. Previous employees have created Image Catalogs for each. Now I want to turn these into services I can share with ArcGIS for Server 10.1. Running the mosaic tool is taking forever, I started at 830am and as of 4pm it is only through 30 of 508 files. There has to be an easier way? Is there a workflow or a way that will be quicker?


st_linemerge does not work

I have two LINESTRING's which intersect by checking with st_intersects. Then i try to collect/union both linestrings into one Multistring, then try to merge with st_linemerge, but it fails with:               "ERROR:  Geometry type (MultiLineString) does not match column type (LineString)"      I check both geometries by using st_astext and get this:               LINESTRING(5.8836331 51.0196776,5.8841849 51.0190699,5.8843962 51.0188219,5.8849422 51.0182963)     LINESTRING(5.8841849 51.0190699,5.8840283 51.0190173,5.8839951 51.0189857,5.8839898 51.0189518,5.8840082 51.0189217,5.8840578 51.0188784,5.8841361 51.0188532,5.8842     492 51.0188324,5.8843962 51.0188219)      I see both have at least one point which intersects (be the same), so it should work. But when i try to merge them it fails or simply returns the multilinestring again. Failing query:               gis=# select st_astext(     st_linemerge(       st_collect(         st_geomfromtext('LINESTRING(5.8836331 51.0196776,5.8841849 51.0190699,5.8843962 51.0188219,5.8849422 51.0182963)'),         st_geomfromtext('LINESTRING(5.8841849 51.0190699,5.8840283 51.0190173,5.8839951 51.0189857,5.8839898 51.0189518,5.8840082 51.0189217,5.8840578 51.0188784,5.8841361 51.0188532,5.8842492 51.0188324,5.8843962 51.0188219)')       )     ));           RESULT:      MULTILINESTRING((5.8836331 51.0196776,5.8841849 51.0190699,5.8843962 51.0188219,5.8849422 51.0182963),(5.8841849 51.0190699,5.8840283 51.0190173,5.8839951 51.018985     7,5.8839898 51.0189518,5.8840082 51.0189217,5.8840578 51.0188784,5.8841361 51.0188532,5.8842492 51.0188324,5.8843962 51.0188219))     (1 row)      What is wrong here?


need help for script checking topology using ogr in python

My main goal is to create a new feature indicating if the line features intersect each other as well as if they are not connected to a feature. My problem is how to identify the intersecting and not connected points and converting them to another geometry layer. The only thing I came up with is:               import ogr          driver = ogr.GetDriverByName('ESRI Shapefile')          inDatasource = driver.Open('test.shp', 0)     inLayer = inDatasource.GetLayer()          #loop through features in the layer     feature = inLayer.GetNextFeature()     while feature:         #condition         feature.Destroy()         feature = inLayer.GetNextFeature()          #check topology          outData = driver.CreateDataSource('TopologyResult.shp')     outLayer = outData.CreateLayer ('TopologyResult', geom_type = ogr.wkbPoint)          fieldDefn = ogr.FieldDefn ('id', ogr.OFTInteger)     outLayer.CreateField(fieldDefn)          featureDefn = outLayer.GetLayerDefn()     outFeature = ogr.Feature(featureDefn)     outFeature.SetGeometry(point)


Join field from table to a raster layer

I'm trying to join a field in a table to a raster layer. I'm using the tool 'Join Field'. The table is in a excel sheet, but have been imported to a file geodatabase. The table has an ObjectID field which i'm using as join field. When i run the tool there is no error message, but the field from the table hasn't been joined to my raster layer. Anyone know what the problem could be?


script for automating insert scale bar, text, north arrow

I've always wanted an automatic process for these tools instead of having to click insert again and again. It seems easy but I can't do it in modelbuilder and I can't seem to find anyone making a script for it on the web. I want to use a certain north arrow and scale bar each time.


What map providers allow export xml data like openstreetmap?

I'm looking for map providers that return xml data about an specific area like openstreetmap export functionality. Do you know any service opensource or comercial that works in that way? Thanks in advance, Pedro


Export KML ExtendedData with OGR2OGR

I would like to export with ogr2ogr the data I have in a kml as ExtendedData                       <Placemark>             <ExtendedData>                 <SchemaData schemaUrl="#15540_punti">                     <SimpleData name="myid">9</SimpleData>                     <SimpleData name="IDENTITA">15540_IX</SimpleData>                     <SimpleData name="NOME">Fontana</SimpleData>                     <SimpleData name="AMP_CONO_V">80</SimpleData>                     <SimpleData name="PROF_1SSP">5,2</SimpleData>                     <SimpleData name="PROF_1P">26</SimpleData>                     <SimpleData name="PROF_2P">78</SimpleData>                     <SimpleData name="PROF_3P">156</SimpleData>                     <SimpleData name="PROF_SFOND">245.26</SimpleData>                     <SimpleData name="ASSE_DIR">117g39p8s</SimpleData>                     <SimpleData name="UNITA PAES">Piazza Santa Croce</SimpleData>                     <SimpleData name="QUOTA_SLM">276 + 1.50</SimpleData>                 </SchemaData>             </ExtendedData>             <Point>                 <coordinates>11.26078815736627,43.76921971994336,0</coordinates>             </Point>         </Placemark>         <Placemark>      I have tried with gdal 1.9.2 and libkml to export in CSV, but I obtain a file that does not contain the fields of the original kml (myid, IDENTITA, NOME, etc.). This is the output               Name,description,timestamp,begin,end,altitudeMode,tessellate,extrude,visibility     ,,,,,,-1,-1,-1     ,,,,,,-1,-1,-1     ,,,,,,-1,-1,-1     ,,,,,,-1,-1,-1     ,,,,,,-1,-1,-1     ,,,,,,-1,-1,-1      Is it possible with ogr2ogr to export KML ExtendedData in tabular form? Thank you, Andrea


How to efficiently read a kml file into R

Hi I'm working with huge .kml files (up to 10 Gb) and need an efficient way to read them into R. Until now I've been converting them to shapefiles via QGIS and then back into R with readShapePoly and readOGR (the latter, by the way, is ~1000 faster than the former). I'd ideally like to cut-out the QGIS intermediary stage as it is cumbersome and slow. So, the question is actually very simply: how to read .kml files in directly. I see this can be also be done with readOGR. Unfortunately, I cannot see how to implement the worked example (after lengthy preparation of .kml file: `xx <- readOGR(paste(td, "cities.kml", sep="/"), "cities")`). It seems that "cities" here is the name of the spatial objects. Roger Bivand admits that "How one discovers this name is not obvious, since the KML driver in OGR needs it to access the file. One possibility is:               system(paste("ogrinfo", paste(td, "cities.kml", sep="/")), intern=TRUE)      " But this does not work for me either. Here's a test .kml file to try it on. With it in my working directory, `readOGR("x.kml", "id")` generates this error message:               Error in ogrInfo(dsn = dsn, layer = layer, encoding = encoding, use_iconv = use_iconv) :        Cannot open layer .       And `system(paste("ogrinfo", "x.kml"), intern=TRUE)` generates:               [1] "Had to open data source read-only."   "INFO: Open of `x.kml'"                    [3] "      using driver `KML' successful." "1: x (3D Polygon)"        , which I simply don't understand. P.s. would `getKMLcoordinates` {maptools} be a valid alternative? I've also tried this:               tkml <- getKMLcoordinates(kmlfile="x.kml", ignoreAltitude=T)     head(tkml[[1]])     tkml <- SpatialPolygons(tkml,                              proj4string=CRS("+init=epsg:3857"))      The coordinates are generated correctly, but my attempt to convert them back into a polygon object failed with the following message:               Error in SpatialPolygons(tkml, proj4string = CRS("+init=epsg:3857")) :        cannot get a slot ("area") from an object of type "double"


Copying schema of Postgis or Spatialite DB to blank one

I want to copy a **Postgis or Spatialite DB sctructure** (features, tables) to a **blank DB**. **In ArcGIS** i can easily handle it by using **XML Workspace Document**. Fisrtly extact it from a **Geodatabase** and than import it to **blank/empty/new** whatever you say. By the way the **stucture of first DB copied to second DB (without data)**. Features (point, line, polygon etc), tables comes with this XML document. So how can i do this for **Postgis or Spatialite DBs** (opensource dbs). **(I dont want data just structure.)**


Looking for open source software similiar to ESRI City Engine

I am looking for open source software similar to ESRI City Engine and currently going through open source software directory. **Update :** About ESRI City Engine (Quick introduction) We want to make few historical places in 3 D format & I heard that in some cases its difficult to integrate with other software (i.e. ESRI City engine to Google earth and vice verse etc) Most of the spatial data is in PostgreSQL/PostGIS . Hence looking for alternative so any open source software included 3 D editing facility will be great. thanks in advance : )


Projection WGS 84 is slightly displaced referring to a hillshade in QGIS

I have a hillshade and a vector layer. The vector layer is slightly displaced referring the hillshade. About lat=100m and long=175m The CRS of the project is EPSG:3397 - PD/83 / Gauss-Kruger zone 4. The CRS from the vector layer: `+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs`. How can I place the vector layer on the right place in QGIS I tried it with the `+towgs84=...` attribute. But this didn't work.


Problem opening GRASS GUI: libgeos-3.3.6.so: cannot open shared object file:

I have been using GRASS to process some Landsat imagery but very recently, and as far as I am aware - after having done nothing different or changing the settings, I am trying to use GRASS and on start-up I receive the message:               GRASS 6.4.2 (newLocation):~ > 3D view mode: libgeos-3.3.6.so: cannot open shared  object file: No such file or directory     3D view mode: libgeos-3.3.6.so: cannot open shared object file: No such file or directory     3D view mode: libgeos-3.3.6.so: cannot open shared object file: No such file or directory     g.list: error while loading shared libraries: libgeos-3.3.6.so: cannot open shared object file: No such file or directory     g.list: error while loading shared libraries: libgeos-3.3.6.so: cannot open shared object file: No such file or directory     GRASS module 'g.region' not found. Unable to start map display window.      Consequently the GRASS GUI does not open once I have chosen my location (I have tested a range of locations that previously opened fine). I am using NX client to run GRASS. Would anyone be able to help with this error message?


Does ITopologicalOperator.Simplify call IGeometry.SnapToSpatialReference?

Does the ITopologicalOperator.Simplify method call IGeometry.SnapToSpatialReference as part of the simplify process? I thought I saw some ESRI documentation that indicated that it does, but I cannot find it now. I would prefer not to make a redundant call if it is not needed. If you have a link to ESRI documentation indicating that it does, please include it.


Comprehensive Intro to GIS Data Structures

I am looking for a comprehensive introduction to design and analysis of spatial data structures. **If possible, please only link to online resources unless the offline resources is especially of note.**


Best way to animate radar data in openlayers

I have a project using OpenLayers which displays radar data on a map. The radar data itself is just a bunch of polygons each of which is filled with one of 6 possible colors. The goal is to animate the radar on the map. Each file contains radar data for a given time and the files are separated by about 5 minutes so my current approach is to loop through all the files and load them one by one into new separate layers. Once each layer is created it has its visibility set to false and it's added to the map. I then animate the layers using a timer which turns on the visibility of one layer and turns off the visibility of the preceding layer. Currently the layers are all vector layers and the data is loaded from KML files although the data files can be pretty much any format that will work best for this project. The problem with this approach is that once I hit a fairly large time span (around 3 hours or so) of data (equating to approximately 36 layers) the memory consumption gets pretty high (around 250mb). The final product is supposed to be able to handle up to 18 hours of data in one loop which based on the above number would require over 1GB of memory just for the animation and would likely crash the browser or at least make it very sluggish. I attempted the same thing using WMS layers for each layer but the redrawing was too slow (the animation changes layers every 100ms) and the memory consumption wasn't much better than the vector version. I have scoured the net trying to find some example of how to do this but so far am coming up empty handed. Does anyone out there have any suggestions on how to do this? I'll take any ideas you have at this point as I've been stuck on this for weeks now with no answer in sight.


Leaflet, defining a specific region for a canvas to draw upon

Using leaflet, you can draw on a given map by adding a Canvas layer. Some instructions are given here, http://leaflet.cloudmade.com/reference.html#tilelayer-canvas               var canvasTiles = new L.TileLayer.Canvas();          canvasTiles.drawTile = function(canvas, tilePoint, zoom) {         var ctx = canvas.getContext('2d');         // draw something on the tile canvas     }      However, I'd like to to draw at a specific lat-long location on the map, can that be done?


Vertex and Nodes Tool - Migrating from 9.3.1 to 10.2.1

A previous employee built a custom tool for users in our organization that displays nodes and vertices when editing. Essentially the purpose is:               The Display Topology Nodes button will symbolize the Topology Nodes of the active     (selected) Feature Class in a Map Topology depending on the number of line segments that share that node.          The Display Vertices button will display the vertices of the active (selected) layer in the Table of Contents.     Nodes will be portrayed as such:          RED = End Node     PURPLE = Where two (2) features intersect          CYAN = Where three (3) features intersect          GREEN = Where four (4) features intersect          MAGENTA = Where five (5) features intersect          YELLOW = Where six (6) features intersect      My question is: Is there a built in tool using ArcGIS Desktop 10.2.1 that will provide this functionality? I am getting a little pushback on upgrading these users and I need to provide this service, however, it looks like the tool as it is now, does not work on 10.2.1.


How to enable user to specify field mappings from a table joined to a feature class?

I have a feature class and a new table every few days. I need to have the table's attributes added to the feature class based on a join. What would be a good interface in ArcMap to let the user map the table's fields to the feature class so that after the join, the tables' attributes can be moved into the feature class. If a field is not mapped then a new field in the feature class would be created and the attributes moved over. I am looking at building a script if another solution doesn't exist. I am looking at the fieldmappings objects but not sure if it will work.


How to import r.in.gdal from WCS (only region extent)

**Can I import WCS-Data only for the current region by`r.in.gdal`?** I tried to import raster data from my WCS into GRASS GIS, using r.in.gdal as proposed in http://grasswiki.osgeo.org/wiki/Global_datasets#OGC_WCS_- _Albedo_example . Then I added an BBOX definition to the XML file as follows:               <WCS_GDAL>     <ServiceURL>http://localhost/mywcs?</ServiceURL>     <CoverageName>blw.eros_z</CoverageName>     <Timeout>90</Timeout>     <GetCoverageExtra>&BBOX=2596000,1181000,2597500,1182500</GetCoverageExtra>     </WCS_GDAL>      However, r.in.gdal just starts to load the whole dataset in 2MB chunks, ignoring my BBOX argument. r.in.gdal also lacks a `-r` option, what would be even better. Altough, when I try to manually download my region extent with a browser, it works well and I get the wanted GeoTiff file for my extent:               http://localhost/mywcs?SERVICE=WCS&VERSION=1.0.0&REQUEST=GetCoverage&coverage=blw.eros_z&FORMAT=image/tiff&CRS=EPSG:2056&BBOX=2596000,1181000,2597500,1182500&RESX=2&RESY=2      Background: I have large raster files (12 x 1 GB, erosion risk data of CH http://s.geo.admin.ch/bb8ef7c7f) with overlapping areas. Instead of patching everything together with r.patch (did work well for 2 rasters, less for >10...), I want to get only the area of interest per WCS (Mapserver 6.4.0).


How to determine why feature is invalid according to Validate features?

I have a feature class with a lot of domains and a lot of different properties e.g. null values, subtypes, and so on. I am using Validate features (Editor tool) and get the selection of the invalid feature(s), but I cannot figure out why some features are invalid. Is there a way to get a report or something? ![enter image description here](http://i.stack.imgur.com/QDSuZ.png)


Why is QChainage not assigning linear referencing sequentially from start to finish?

I have a simple linear feature representing a pipeline that needs linear referencing every kilometer. When I add the specifics to the tool it runs fine but my result does not count sequentially from the start of the line to the end. My results begin at zero and go up to 32km but then jump up to 79km. the end of the line starts at 33km and then counts sequentially backwards along the line back to the 79km marker. I am dealing with a simple single feature shapefile, not a multipart or anything and I have tried dissolving again but that doesn't seem to do anything ( as it shouldn't) and have tried copy to a new file but still end up with the same results. Any help or thoughts on this would be appreciated.


Error adding coordinates in CSV files

I'm wondering if anyone can help me. I added a CSV file on ArcGIS 10.2.1, I displayed the coordinates and everything worked fine. However, when I duplicated and edited the new CSV file to change the numbers in one column it's not working. Instead of my XY columns I keep getting the choice of lots of fields instead. Does anyone know how to fix this? It has been an intermittent problem through my project. ![enter image description here](http://i.stack.imgur.com/4Eo51.png) Thanks for your reply. I'm editing the file in Excel and saving as to create a new file. I then restart ArcMap and add the file but the same error happens. Any ideas? The original file still works as normal. The names of the columns are ID, X, Y, OFFSETA, OFFSETB, RADIUS2. The headings don't change with the new files I create, just the figures in one of the columns.


How to achieve better Dissolve results?

I applied the Vector>Geoverarbeitungswerkzeuge>Auflösen (Dissolve) on a zoning plan (see picture zoning plan). I used the attribute village name as the dissolving criteria. The resulting geometry for the village (yellow colored) is sprinkled with small polygons or also larger polygons are not dissolved (see picture sprinkled geometry). In the attribute table the information is in one row. Are there any approaches to obviate this? Changing any settings? I am using QGIS 1.7.2-Wroclaw. ![zoning plan](http://i.stack.imgur.com/bYAce.png) ![sprinkled geometry](http://i.stack.imgur.com/ucopG.png)


Is it possible to modify the ArcObjects IProgressDialog2?

It turns out that the best way to run long tasks in ArcMap is to use the IProgressDialog2. The BackgroundWorker in combination with ArcObjects (STA Threads) is a mess. However, it seems that the only implementation from ESRI is the IStepProgressor. In my case I do not need and do not want to show a progress bar. Moreover I do not need the cancel button. Is there a way to remove both things? I want only the text and the circle as you can see in my screenshot ![my progress dialog](http://i.stack.imgur.com/VRIFS.png)               public static void ShowProgressDialog<T>(string message, Action<T> action, T arg)         {                  ESRI.ArcGIS.esriSystem.ITrackCancel trackCancel = new ESRI.ArcGIS.Display.CancelTrackerClass();                  ESRI.ArcGIS.Framework.IProgressDialogFactory progressDialogFactory = new ESRI.ArcGIS.Framework.ProgressDialogFactoryClass();                  // Set the properties of the Step Progressor             System.Int32 int32_hWnd = ArcMap.Application.hWnd;             ESRI.ArcGIS.esriSystem.IStepProgressor stepProgressor = progressDialogFactory.Create(trackCancel, int32_hWnd);                       ESRI.ArcGIS.Framework.IProgressDialog2 progressDialog2 = (ESRI.ArcGIS.Framework.IProgressDialog2)stepProgressor; // Explict Cast                  progressDialog2.CancelEnabled = false;             progressDialog2.Description = message;             progressDialog2.Animation = ESRI.ArcGIS.Framework.esriProgressAnimationTypes.esriProgressSpiral;                  action.Invoke(arg);                  trackCancel = null;             progressDialog2.HideDialog();             progressDialog2 = null;              }


OpenLayers: Showing WKT data via JSON from MVC

I am trying to read geographic data from MS SQL Server 2012 and display it using the JavaScript library OpenLayers in an MVC4 view. I can confirm that the data leaves MVC OK but I have problem is getting the data to be processed by OpenLayers. I have tried two approaches: a) using OpenLayers.Protocol.Script and b) using $.getJSON to read the data and adding using addfeatures. Both fail for different reasons. I will describe the method using OpenLayers.Protocol.Script as it is the preferred method. My MVC View looks like this:
Open Layers - Two layer from SQL, Massive

          @section scripts {                         }      And my MVC action that responds to the OpenLayers.Protocol.Script looks like this (note that the field 'Location' is of type DbGeography and .ToString() returns the WKT version of the geography). The service simply returns a IQueryable of type IGeoItem.               public JsonResult GetGeoLayerData( int layerId, IServiceGetSpecificList<IGeoItem, int> service)         {             var formattedData = service.GetSpecificList(layerId).Select(x =>                  new { wktgeography = x.Location.AsText(), value = x.Value, info = x.Info });             return Json(formattedData, JsonRequestBehavior.AllowGet);         }      What happens is:   1. A request to the MVC action happens and a valid JSON result is handed back (as seen by FireBug NET console).   2. The read function never gets called (checked with Firebug breakpoint and Console.Log message).   3. ... Therefore no vector features are added to layer1. I would appreciate help on solving this as I have tried a number of approaches and none of them work. Maybe I have made a silly mistake but I can't see it. Thanks in advance. **UPDATE with solution:** Following the advice from cavila (see below) I have fixed this problem. The change was in the definition of the Vector layer protocol from Script to HTTP. I have included just that section of the view below in case it helps anyone else               var layer1 = new OpenLayers.Layer.Vector("Layer 1", {         projection: proj_4326,         strategies: vector_strategies,         protocol: new OpenLayers.Protocol.HTTP({             url: "GetGeoLayerData?layerId=1",             format: new OpenLayers.Format.JSON({                 read: function(jsonp) {                     var features = [];                     $.each($.parseJSON(jsonp), function(i, item) {                         var feature = wktReader.read(item.wktgeography);                         feature.attributes = { value: item.value, info: item.info };                         features.push(feature);                     });                     return features;                 }             })         })     });


openlayers trigger hover event

I have a bunch of tracks on screen in a vector layer. I used some hover and select styles in the layer styleMap and these styles are properly triggered by:                   var highlightCtrl = new OpenLayers.Control.SelectFeature(tbTrack.vector, {             hover: true,             highlightOnly: true,             renderIntent: "hover",         });              var selectCtrl = new OpenLayers.Control.SelectFeature(tbTrack.vector,             {clickout: true, renderIntent: "select"}         );              tbMap.map.addControl(highlightCtrl);         tbMap.map.addControl(selectCtrl);              highlightCtrl.activate();         selectCtrl.activate()      All these tracks are listed in a table under the map. I want to trigger the corresponding hover and select styles for a track if I click or hover on them in the table. How to do this? I allready have the coupling between the table and the lines, I just don't know the right way to trigger the hover or select events.


Efficient computation of nearest distance for each point to closest polygon

I'm looking for a more efficient way to compute the nearest distance between features from 2 layers, in this example for each point the nearest distance to the closest polygon. Is there a more efficient way compared to below where I compute all polygon distances for each point? (Python 2.7.2, GDAL/OGR 1.9.1)               nPts = pointsLayer.GetFeatureCount()     nPolys = polysLayer.GetFeatureCount()               for pt in range(0,nPts):           minDist = 1000000           pointFeature = pointsLayer.GetFeature(pt)           for poly in range(0,nPolys):             polyFeature = polysLayer.GetFeature(poly)             Dist = pointFeature.GetGeometryRef().Distance( polyFeature.GetGeometryRef() )             if (Dist < minDist):                   minDist = Dist      print (pt,minDist)


WMS layers don't draw in QGIS or ArcMap but do in OpenLayers

WMS Map was published by myself in AutoDesk Infrastructure Map Server 2013. I can connect and load detials (properties / metadata) in QGIS Browser, ArcCatelog. I can load the wms into QGIS desktop and ArcGIS desktop but the layers do not draw (in the mapframe) in either program. I can isolate and draw the layers in OpenLayers html file just fine. The only thing that concens me is a SQL Spatial layer in the Map created using ArcSquirrel. Any help would be great.


How to retrieve the feature layer metadata from ESRI Geoportal Server?

Is there any option to retrieve the metadata information of a particular layer (within the ArcGIS Server Service) registered with ESRI Geoportal Server?


ArcGIS: How to programmatically refresh selected features

**Background:** I have C#, ArcGIS command that 1. lists selected features, and 2. allows you to delete selected features from this list out of a geodatabase. **My problem:** After I delete features, ArcGIS doesn't update the selected features list. When my program asks ArcGIS for the selected features, it gets back references to the features it already deleted. It shouldn't. **My Question** : What should I do to tell ArcGIS to check the selected features list for non-existent features? Is there a "refresh" command that will do this, or do I have to reselect the features (excluding the deleted ones) programmatically?


How can i generate parcel extracts on demand through a webservice

Ok, i need to generate parcel extract documents on demand. The backend is Arcgis 9.3.1 on SQL server 2008. I would like to build a webservice (C# .Net) on the server that generates these extracts. My experience is writing with arcobjects to extend Arcmap. But now i just want to access the data and generate without this UI. What are the options i have? Any suggestions would be great. (p.s. it is not necessary to stick to .NET) I only need some advice or direction. I'll try to describe my problem in general. The extracts will be requested over the internet. So the only thing the system should do is returning an jpg or pdf. That's all. The problem i have now is that i don't now how to do it. Any concept or advise would be great. Sorry if this is still to generic. I have to format the extract from a template so a grid and a north arrow is added. I know how to do this with arcobjects in arcmap but i don't know how to use arcobjects from a webservice and i even don't know if this is possible.


Node Visibility Within Layer

Is there a way within QGIS to hide/display (freeze/thaw) an individual node feature within a layer or does this control only exist at the layer level.


PostGIS- GROUP BY for continuous data values?

In the sample code for the PostGIS ST_ConvexHull function, multiple polygons would be created based upon different, discrete values found in the field "disease_type".               --Get estimate of infected area based on point observations     SELECT d.disease_type,         ST_ConvexHull(ST_Collect(d.the_geom)) As the_geom         FROM disease_obs As d         GROUP BY d.disease_type;      Lets say there was a numeric field (lets call it "myfield") with continuous values from 0-5. How would you produce similar results as GROUP BY but defining your own breakpoints (for example: 0-1.25, 1.25-3.5, 3.5-5)?


pair value comparison to generate new raster using arcpy/python

I want to write a function which takes in two arguments (from two rasters) and compare those 2 values against a set of hard coded pair values inside function. Then return a value which i want to use to generate a new grid. I am not quite there to complete the function. Any suggestions?               def grid(val1, val2):         table = {'10','1':'75','10,2':'75','10,3':'83','10,4':'87'}         for i, j, k in table.iteritems():             if val1 == i:                 if val2 == j:                     val = k         return val          t = grid(10,1)


How can I draw squares in PostGIS layers at different angles?

I am trying to draw two squares on a 45-degree angle and then merge them using PostGIS to edit. QGIS can draw squares but only on a perpendicular plane (eg, 0, 90, 180), and I need to be able to draw at different angles (eg, 45-degrees). Thank you!


How to label lines in Tilemill?

I am designing my map with tilemill. Now I am struggling with data visualization. I would like to visualize the road name like the 1st image instead of the second image. I have loaded the data source of tilemill from shapefile which contain Road name that field I want to visualize. Can anybody help me on this issue, whether it is possible or not to make the annotation from data? ![I like to visualize](http://i.stack.imgur.com/Y29qR.jpg) ![But it is showing like this](http://i.stack.imgur.com/paLD7.jpg)


TileCache SVN trunk on IIS 7 - An error occurred: No section: 'cache'

I was having tile generation problems with TileCache 2.11 and decided to update to the trunk. I repeated the instructions necessary to get TileCache working in IIS but now I'm seeing this error:               An error occurred: No section: 'cache'       File "C:\inetpub\wwwroot\tilecache\TileCache\Configs\File.py", line 147, in read         self.cache = self._loadFromSection(config, "cache", TileCache.Cache)       File "C:\inetpub\wwwroot\tilecache\TileCache\Config.py", line 74, in _loadFromSection         type  = config.get(section, "type")       File "C:\Python25\lib\ConfigParser.py", line 511, in get         raise NoSectionError(section)      I've made the necessary change to tilecache.py:               from TileCache import Service, cgiHandler, cfgfiles          if __name__ == '__main__':         svc = Service.load("c:\\inetpub\\wwwroot\\tilecache\\tilecache.cfg")         cgiHandler(svc)      And I've confirmed my configuration file exists and does contain a [cache] section:               [cache]     type=Disk     #base=c:\inetpub\wwwroot\tilecache\Cache     base=g:\tilecache\Cache      Using SysInternals ProcessMonitor I've determined that python.exe never attempts to access tilecache.cfg. Renaming tilecache.cfg results in the same error. Why isn't it loading the file? Is the path incorrectly formed?


Link/unlink buttions in the Attribute Table's Form View

Toggling the edit button, in the Form View/Table View, buttons for adding/deleting a record is available. Two additional buttons are also available (Fig. 1). From the look of the icons they perform some sort of link and unlink function between records. What I can't figure out is what I can "gain" by linking records. Is it for the purpose of editing several records in a similar manner simultaneously? Is it a purely temporary UI link or is it persisted in the DB if I save afterwards? And I can't seem to get unliking to work after linking several records. ![Link buttons](http://i.stack.imgur.com/T8Mnn.png)


About the shortest path returned by pgrouting

I have a couple of questions related to the shortest path returned by pgrouting. I have this result for example.               edge_id   cost            reverse_cost  x1           y1         x2           y2     1         0.0016861729    1000000       -93.2518843  44.9718147 -93.2527683 44.9707752          160       0.0034403359    1000000       -93.2544283 44.9686054  -93.2527683  44.9707752      Now as you can see the first two results of the shortest path query are showing edges which are one way. It is because the reverse cost is high 1000000. Now if I want to trace out the path of the shortest path then I will join the coordinates like this               (-93.2527683  44.9707752) -> (-93.2527683 44.9707752) -> (-93.2527683 44.9707752)      and soon. Now my question is the coordinate pair (-93.2527683 44.9707752) is common in the first and second edge which are both one ways. Now when joining the coordinates from the second pair, I am joining it like               (-93.2527683 44.9707752) -> (-93.2527683 44.9707752)      which is in opposite direction, so is it fine? I mean the edge is one way but I am joining the coordinates from like (x2,y2) to (x1,y1). Is this correct or I am missing something?


TileCache error when using import arcpy

I recently ran the 10.1 sp1 and my scripts stopped working. I can't even right click on the scripts and check properties without the following popping up:               InitializeParameters Syntax Error:  Runtime error     Traceback (most recent call last):       File “#ANNOBO2.InitializeParameters.py”, line 1, in module     File “d:\arcgis\desktop10.1\arcpy\arcpy\__init__.py”, line24, in module     From arcpy.toolbox import*     File “d:\arcgis\desktop10.1\arcpy\arcpy\toolbox.py”, line 485, in <module>     From management import ExportTileCache as ExporttileCache_management     ImportError:  cannot import name ExportTileCache      I did a little digging and found all the tools in python script in the management.py. What was missing was the scripts for the Tile Cache toolset. Not sure what happened when I installed sp1 but now I just want it fixed. Any help would be greatly appreciated. Thanks in advance. Currently using Python 2.7 Ah, ok. Ill look and see if there are any reminants of older python on my system. If not ill try to reinstall everything. Thanks for the quick replies.


display WMS layer from geoserver in mobile app with openlayers mobile

i try to display WMS layer from geoserver with openLayers Mobile,and its working great in the navigator (chrome) but its not appeering in the simulator device (android)...   & with another WMS layer like open street map , its working perfectly in both, web navigator and android device , so i can't figure out where the problem is , its been more then 10 days and i can't handel this problem.. this is my code:               var map;         OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";              function init() {             map = new OpenLayers.Map('map', {                     controls: [                         new OpenLayers.Control.Navigation(),                         new OpenLayers.Control.PanZoomBar(),                         new OpenLayers.Control.LayerSwitcher({'ascending':false}),                         new OpenLayers.Control.ScaleLine(),                                             new OpenLayers.Control.MousePosition()                                            ]                        });                      var wms = new OpenLayers.Layer.WMS(                     "OpenLayers WMS", "http://vmap0.tiles.osgeo.org/wms/vmap0",                     {layers: "basic"}                 );                     var osgeo = new OpenLayers.Layer.WMS("OpenLayers WMS",                     "http://vmap0.tiles.osgeo.org/wms/vmap0",                     {layers: 'basic'},                     {isBaseLayer: true, transitionEffect: 'resize'}                 );                    com = new OpenLayers.Layer.WMS(                 "communes",                 "http://localhost:8080/geoserver/wms",                 {layers: 'cite:matcommune_view', format: 'image/gif',transparent: true}             );                 map.addLayers([osgeo, com]);          map.setCenter(new OpenLayers.LonLat(-104, 42), 3);      im looking forward your help, tks.


How to Add Shapefile with Symbology from Layer File to ArcMap MXD using ArcPy/Python?

Say I have a Shapefile: C:\DATA\points.shp A layer file: C:\SYMBOLOGY\symbology.lyr And a map document: C:\MAPS\map.mxd How can I add points.shp to map.mxd with the properties of symbology.lyr all using python? Again I want to thank everyone in advance for their help.


Is it possible to recover ArcGIS Online tiles of a zone, then reproject and republish them in a Map Service?

I'd like to use ArcGIS Online basemaps but not is the Web Mercator projection. Because of that I want to know if it's possible to reproject and republish them in a Map Service.


Which methods are most suited to classifying urban areas from Landsat imagery?

I am trying to use Landsat data to map urban extents over time. I need to perform a classification into built up / non built up areas. I am new to image classification. I have access to the following software:   * ArcGis   * ENVI   * Ecognition   * GRASS / any other open source software you can recommend. Which methods and applications are most suited to this analysis? I would appreciate if you could help me understand which methods or applications yield the highest classification accuracy.


ArcGIS Server tiled layer doesn't update on pan in OpenLayers

I'm trying to add an ArcGIS Server tiled/cached layer to OpenLayers. This is working using code modified from these samples. The only problem is that when panning the map, the layer does not update to show data which was previously off the screen. In this example I've panned the map to the west, and the overlay data hasn't updated: ![enter image description here](http://i.stack.imgur.com/8N8QQ.png) However, if I zoom (rather than pan) the layer does update. Do I need to explicitly call an update method on the layer when panning? Note that the first example script correctly updates after panning (but this is a basemap, not an overlay - does that make a difference?) My sample code is available here. I need to hard-code the layerInfo parameters, rather than retrieving them at run-time (I obtained the values using ..../MapServer?f=json&pretty=true) Thanks for any advice. [EDIT: the layer _does_ update when panning in Chrome or IE9, but not in FireFox or Opera.] [EDIT2: if I switch the layer off then on in the moveend function, the layer correctly updates after panning in FireFox and Opera:               map.events.register("moveend", map, function() {         var lyrs = map.getLayersByName("geology");        if(lyrs.length > 0) {           lyr = lyrs[0];           lyr.setVisibility(false);           lyr.setVisibility(true);                        }                                              });


Map Service Tile Cache: ArcGIS Server 10.2

If I delete my old map services, will the associated Map Service Tile Cache be deleted as well? For some reason the Map Server Delete Tile Cache tool is not working for me, and our tiles are cached to our server at this moment. We need to make space on the server for new services, so I am affraid to delete the old services and have their tile cache remain on the server. Am I safe to just delete the old services? Or is there a way to access the cache in the server manually so I can delete it before proceeding?


can QGIS for Mac OSX 10.6 handle (both read and generate) GeoPDF files?

I would like to know if QGIS for Mac OSX 10.6 can handle (both read and generate) GeoPDF files. Can anybody help ?


Interpolating and infilling a raster from 5km pixels to 100m pixels in Arcmap 10

I am trying to refine a raster dataset of annual rainfall from a 5km grid down to a 100m grid. Before I continue, I am aware that this could result in some uncertainty regarding the data. However, the data is only available in the 5km grids, and this is of no use to me, as I am looking to compare rainfall to landuse, and the landuse map is to a 100m grid. To get around this, I thought I would be able to infill the data through the use of contours (as rainfall data is generally gradient based, rather than discrete isolated data points) I contoured the 5km raster file using Contour(Spatial Analyst), and then tried to use the contour file to create a new raster file (Feature to Raster) at a 100m resolution. However, because the contour feature is a polyline feature, there isn't any data between the contours. This means that anywhere the contour lines fall more than 100m apart, the raster file leaves a blank pixel. As I see it, there are three options that I may be able to look into (no idea how, or if, the first two would work) Option 1: Directly refining the resolution of the original raster file, by increasing the resolution from 5km to 100m pixels and infilling the new pixels with values interpolated from the original file. Option 2: Some sort of interpolation tool, (I saw several in the Geostats toolbox, but don't understand how they work) that could transform the contour polyline feature into a continuous dataset, which could then be converted to raster. Option 3: Create a contour file with ridiculously fine contours (I have already tried 0.5 mm, over a range from ~ 500mm to 3000mm, and this didn't do the trick in places) Any help would be appreciated. Cheers.


Determine if GPS coordinate in Canada

I have GPS tracking solution and most my customers in US. My reverse-geocoding solution pretty basic:   1. Table with US ZIP database   2. Procedure in SQL Server where I search for city/state closest to reported lat/lon Now I need to add Canada to my service coverage. I will buy Canadian postal codes and plan to store it in separate table. So, straight-forward solution will be to:   1. Check if there is close US city/state and report   2. If not found - check in Canada table  This solution will put extra load on server for those canadian users since US table will be queried first. For various reasons I don't want to merge tables. However, if you look at map - you can tell that Canada is above USA :) So, I was wondering if I can use some hard-coded numbers to know immedeately where this point is. Right? Wrong? Other ideas on how to tackle this?


Routing issues with pgrouting

I was using the pgrouting plugin to display the shortestpath. I used the plugin as suggested by underdark to display the shortest path in qgis. My source node is 283 and my target node is 172. When I used the pgrouting plugin to give me the shortest path I got one as attached in the screenshot. ![enter image description here](http://i.stack.imgur.com/Hnrv7.png) However, when I queried the database of the graph which was generated by the tool osm2po. I found only one edge which was connected to node 283. The result of the query is given in the figure below ![enter image description here](http://i.stack.imgur.com/Y5TMf.png) My question is if you look at the results of the query, the edge is from node 279-283 and is one way as well with 283 being the target. Now since I used qgis to show the shortest path from 283-172. There is no starting from 283 as source. So my assumption was it shouldn't have shown any path. So I was wondering if there are some issues with the pgrouting plugin. How come it gave me the shortest path when there is no edge with 283 as source? As I have explained there is only one edge connected to 283, but 283 is target and also the edge is one way from source to target.


How to rename feature classes based on field values within the attribute table?

My goal is to iterate through several feature classes and rename each feature class based on a field value within the feature class's attribute table. These feature classes are a single polygon feature, thus they only have a single row in the attribute table. I was thinking I could use the Rename GP tool to do this so I've got that part figured out. Next I can use the SearchCursor to find the field value that I want to use to rename the feature class with. So I'm struggling here in two places. First, I'm unsure how to iterate through each of the feature classes to get their current names and then apply that to the first parameter of the arcpy.Rename_management tool. I know I can get the name of the feature class using arcpy.Describe but I'm stuck after that. Second, I'm struggling with how to assign each of the field values from the attribute table to a variable to be used for the second parameter of the Rename tool. I can easily print the name using `rows = arcpy.SearchCursor(fc, "", "", "name") for row in rows: print row.name` but I'm struggling with how to assign its value to a variable. Let me know if I need to explain this any better.


ArcPy: Initialize zero valued grid Raster and iteratively increment

I'm jumping right into `ArcPy` to help out a friend who needs to automate some ArcGIS work. The problem at hand can be reduced to the following pseudocode:               accumulationRaster = SomeZeroedRaster     for each Raster pair, R, S, pertaining to rain and snow:           accumulationRaster += Con(S == 0, Con(R >= 1, 1, 0), 0)      In other words, iterate through all data of rain and snow of each day. Perform a logical operation on each cell pair, so that the corresponding cell in the `accumulationRaster` gets incremented if no snow (`S == 0`) and some rain `(R >= 1)`. Having no luck so far with various attempts, I thought I'd ask here, as the problem should have an easy solution. Full working code with the missing gaps pointed out below. Most of it is just trivial file path set up, and the above pseudocode should suffice: The missing pieces are numbered `1`, `2`, and `3`.   1. How to initializase a zero raster w/proper dimensions.   2. How to make a copy of this zero raster. If `1` is simple, just do this instead.   3. Increment accumulation raster. * * *               import os     import arcpy     from arcpy import env     from arcpy.sa import *          in_workspace = "C:/Data/ProjectFolder"     os.chdir(in_workspace)     arcpy.env.workspace = in_workspace          maaneder = ["01_januar", "02_februar", "03_mars", "04_april",                 "05_mai", "06_juni", "07_juli", "08_august",                 "09_september", "10_oktober", "11_november", "12_desember"]          ################ 1 ###################     baseAccum = #1: ?? (zero Raster w/correct type & dimensions)     ######################################          outputPath = "rainDaysInGrowthSeason/"     for year in range(1997,2015):         for month in range(1, 13):             ################ 2 ###################             monthAccum = #2: ?? (Copy of baseSum)             ######################################                  for day in range(1,32):                 # e.g 1997/09_september/                 fpath = str(year) + "/" + maaneder[month-1] + "/"                 # e.g. 1997_09_08.asc                 filesuffix = str(year) + "_" + str(month).zfill(2) + "_" + str(day).zfill(2) + ".asc"                 # e.g. Snoedybde/1997/09_september/sd_1997_09_08.asc                 snoefile = "Snoedybde/" + fpath + "sd_" + filesuffix                 # e.g. rrl/1997/09_september/rrl_1997_09_08.asc                 rrlfile = "rrl/" + fpath + "rrl_" + filesuffix                      hasSnoe = os.path.isfile(snoefile)                 hasRrl = os.path.isfile(rrlfile)                 if (not hasSnoe or not hasRrl):                     continue                      snoeRaster = Raster(snoefile)                 rrlRaster = Raster(rrlfile)                      ############## 3 #####################                 monthAccum = #3: ??                              # something a'la:                              # += Con(snoeRaster == 0, Con(rrlRaster >= 1, 1, 0), 0)                 ######################################                  # e.g. rainDaysInGrowthSeason/1997             outputDir = outputPath + str(year)             if (not os.path.isdir(outputDir)):                 os.mkdir(outputDir)             monthAccum.save(outputDir + "/r_" + str(year) + "_" + str(month).zfill(2))


ArcMap trace the creation/modification/deletion events by windows user for all changes (not only the latest change)

I want to check whether there is a way to trace all the creation, modification and deletion events in the editing session. What I aware now is how to get the latest updated information of an object in the attribute table of a layer. I am using ArcMap 10.2.2.


Spatial joining syntax

A standard PostGIS point-in-polygon query might look like this. (example tutorial)               SELECT p.*,b.polyname      FROM points p, shapes s      WHERE ST_Within(p.the_geom,s.the_geom)      Imagine this was _not_ a spatial join. Then it would look like this:               SELECT i.*,c.categoryname      FROM items i, categories c      WHERE i.zipcode = c.zipcode      Most SQL folks would advise rewriting that query using a variety of the `JOIN` keyword syntax:               SELECT i.*,c.categoryname      FROM items i      JOIN categories c ON i.zipcode = c.zipcode      **Does PostGIS offer this sort of syntax for spatial joins?** I can't find it. All I see are a lot of functions that return `true` or `false`, leading to the workaround syntax. Ideally, I would think a spatial join would look like:               SELECT p.*,b.polyname      FROM points p      JOIN shapes s      ON       p.the_geom IS WITHIN s.the_geom


Raster projection problems

How would I go about specifying the following projection in Arc and QGIS? After defining this, how would I then go about changing the spatial reference to ESPG 27700? The following are the details of a European raster - I need to define the projection and then transform it - every time I've tried it the UK ends up very far away from where it's supposed to be. Format: ESRI ascii grid Map projection: Lambert Azimuthal Equal Area, ellipsoid WGS84 Center point: 48° N, 18° E. Grid cell size: 1000m North: 2800000m South: -1700000m West: -2480000m East: 2520000m Rows: 4500 cells Columns: 5000 cells Missing value: -9999 Thanks -