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

text: 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


title: How to select distinct values from a column in Mapbasic?

text: 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?


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

text: 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.


title: List of free WFS 2.0 clients

text: 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?


title: QGIS opening shp file: why does it crash?

text: 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


title: How to Properly Share as Image Service?

text: 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?


title: st_linemerge does not work

text: 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?


title: need help for script checking topology using ogr in python

text: 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)


title: Join field from table to a raster layer

text: 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?


title: script for automating insert scale bar, text, north arrow

text: 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.


title: What map providers allow export xml data like openstreetmap?

text: 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


title: Export KML ExtendedData with OGR2OGR

text: 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


title: How to efficiently read a kml file into R

text: 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"


title: Copying schema of Postgis or Spatialite DB to blank one

text: 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.)**


title: Looking for open source software similiar to ESRI City Engine

text: 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 : )


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

text: 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.


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

text: 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?


title: Does ITopologicalOperator.Simplify call IGeometry.SnapToSpatialReference?

text: 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.


title: Comprehensive Intro to GIS Data Structures

text: 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.**


title: Best way to animate radar data in openlayers

text: 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.


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

text: 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?


title: Vertex and Nodes Tool - Migrating from 9.3.1 to 10.2.1

text: 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.


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

text: 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.


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

text: **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).


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

text: 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)


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

text: 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.


title: Error adding coordinates in CSV files

text: 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.


title: How to achieve better Dissolve results?

text: 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)


title: Is it possible to modify the ArcObjects IProgressDialog2?

text: 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;              }


title: OpenLayers: Showing WKT data via JSON from MVC

text:
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;                 }             })         })     });


title: openlayers trigger hover event

text: 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.


title: Efficient computation of nearest distance for each point to closest polygon

text: 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)


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

text: 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.


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

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


title: ArcGIS: How to programmatically refresh selected features

text: **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?


title: How can i generate parcel extracts on demand through a webservice

text: 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.


title: Node Visibility Within Layer

text: 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.


title: PostGIS- GROUP BY for continuous data values?

text: 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)?


title: pair value comparison to generate new raster using arcpy/python

text: 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)


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

text: 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!


title: How to label lines in Tilemill?

text: 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)


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

text: 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?


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

text: 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)


title: About the shortest path returned by pgrouting

text: 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?


title: TileCache error when using import arcpy

text: 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.


title: display WMS layer from geoserver in mobile app with openlayers mobile

text: 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.


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

text: 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.


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

text: 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.


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

text: 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.


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

text: 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);                        }                                              });


title: Map Service Tile Cache: ArcGIS Server 10.2

text: 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?


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

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


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

text: 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.


title: Determine if GPS coordinate in Canada

text: 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?


title: Routing issues with pgrouting

text: 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.


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

text: 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.


title: ArcPy: Initialize zero valued grid Raster and iteratively increment

text: 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))


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

text: 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.


title: Spatial joining syntax

text: 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


title: Raster projection problems

text: 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 -


title: ArcGIS Server Manager - access problem

text: I setup ArcServer 10 with postinstall and restart and I am getting the following error when I try accessing ArcServer Manager through URL: > >     Forbidden >   >     You don't have permission to access /arcgismanager on this server. >   I can admire that i am a member of agsadmin group.


title: What cartographic line symbol could make a coastline represent water like this old map?

text: I live in a coastal area with lots of islands. Sometimes with a black and white map it's hard to distinguish land from sea. So I'm trying to replicate this old map's style. ![Old Map with offset coastline](http://i.stack.imgur.com/9svio.jpg) I'm trying to do this in QGIS but if it's easier in ArcGIS I'd be happy to hear how to do it in there as well. At first I tried to make multiple simple lines and offset the duplicates but it offsets them inside as well as outside. So then I made an SVG with offset lines and I got closer to what I'm aiming for but since the SVG is a rectangle shape there are gaps and overlaps at every angle. ![QGIS attempt](http://i.stack.imgur.com/2D0tN.jpg) Does anyone have any tips? Thanks in advance.


title: INSPIRE, GEOSS and OGC

text: (as simply as you can) what are the differences between INSPIRE, GEOSS and OGC? I will explain what I understood: (1) The OGC provides standards for web-mapping that like WMS, WFS (are these standards or services? is there any difference?) Does OGC provide data? (2) If INSPIRE provide data to EU countries and GEOSS/OGC does that but for all the world so what is the the point of INSPIRE, then? is there and differences in SDI architecture between INSPIRE and GEOSS? Sorry for asking such question but I'm really confused about it :S


title: How to merge output shapefiles in model builder

text: I have a layer with several polygons and a road map. I need to take each polygon and use it to clip the road map, dissolve all of the road records into one, and then merge the polygons into one shapefile. The goal is to have a shapefile with each polygon as one record that contains the total roads within that polygon. I know the easiest solution would be to dissolve the roads layer before hand and then clip using the polygons but I'm working on a continental scale and my computer keeps running out of memory for that. If I clip the road map and then dissolve everything is fine, the issue is the merge. Using row iterator the model only spits out one shapefile at a time after running a clip and dissolve on it. How do I collect all of those outputs and merge them? Current Model: I'm using row iterator to select out each polygon, then use the clip tool to clip the road map, then using the dissolve tool to dissolve all of the roads inside that polygon into one record. Now I'm stuck on how to merge it all back together. Thanks!


title: GeoTIFF raster layer, c#

text: We have a map rendering engine (renders from MapInfo Interchange Format). This engine uses Mercator Projection. Now we have to insert a raster layer (orthophotos of a city), we have GeoTIFF files with the .tfw files. The idea: we create ~200x200px png files for each zoom level and the client will ask a WCF service for the images with parameters: zoom level, top left coordinate, bottom right coordinate. My problem is: how to make the vector and raster layers work together? I don't really want to change the aspect ratio of the raster images. Thanks in advance, Mart


title: Multiple Conditional Statement with in a Con Tool

text: Well i want to apply multiple conditional statements within a Conditional tool. My aim is to loop through a list of rasters and on the basis of a specified range assign each range with a constant 1, 2, 3, 4 or 5 value so that all the final rasters have pixels having a DN of 1,2,3,4, or 5. I copied the template from the following Desktop Help website and adjusted to my requirements. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Conditional_evaluation_with_Con/009z00000008000000/ My code is as follows               import arcpy     from arcpy import env     # Set the current workspace     env.workspace = r"C:\Experiments\Test_Jan03-11_A\Masked"     # Get a list of Rasters     rasters = arcpy.ListRasters("*", "ALL")     for raster in rasters:        output = arcpy.sa.Con(raster <= -10, 1, Con(raster <= 73, 2, Con(raster <= 196, 3,     Con(raster <= 403, 4, Con(raster <= 600, 5, 0)))))     output.save("C:/zE/ResearchWork/UnderWorks/Experiments/Test_Jan03-11_A/Masked/Classify")     Print = "Done"     Runtime error <type 'exceptions.NameError'>: name 'Con' is not defined      The error indicates that i am making a mistake with the conditional statement that i cant identify. Any help will be highly appreciated. I am new at this and enjoy this work but hate when get stuck like this :) * * * I am very thankfull to you guys for your help The following code has worked. But with some issues stated below               import arcpy     from arcpy import env     arcpy.env.extent = "MAXOF"     # check out spatial analyst extension     if arcpy.CheckExtension('Spatial') == 'Available':          arcpy.CheckOutExtension('Spatial')     else:         arcpy.AddMessage("Error: Couldn't get Spatial Analyst extenstion, exiting")     sys.exit(1)     from arcpy.sa import *     # Set the current workspace     env.workspace = r"C:\zE\ResearchWork\UnderWorks\Experiments\Test_Jan03-11_A\Masked\CMasked"     # set cell size     arcpy.env.cellSize = r"C:\zE\ResearchWork\UnderWorks\Experiments\Test_Jan03-11_A\Masked\CMasked\GW1AM2_20130100_01M_EQMA_L3SGSMCLA1110100"     # Get a list of Rasters     rasters = arcpy.ListRasters("*", "ALL")     for raster in rasters:         inRas = Raster(raster)     for raster in rasters:         output = arcpy.sa.Con(inRas <= -10,1,Con(inRas <= 73,2,Con(inRas <= 196,3,Con(inRas <= 403,4,Con(inRas <= 600,5,0)))))     output.save("C:/zE/ResearchWork/UnderWorks/Experiments/Test_Jan03-11_A/Masked/CMasked/MaskedClassify")     Print = "Done"      The problem is that the result is only one raster i was expecting each of the input raster to be classified and saved seperately. Also it is giving me the following error for saving data               Runtime error <type 'exceptions.RuntimeError'>: ERROR 010240: Could not save raster dataset to C:/zE/ResearchWork/UnderWorks/Experiments/Test_Jan03-11_A/Masked/CMasked/MaskedClassify with output format GRID.


title: Open KML Photos In Default Desktop Application

text:
I have created a KML file containing several thousand placemarks. The description contains a simple HTML tag referencing a location on local disk, pointing to a JPG photograph. The photos display in the popup balloons properly. I would like for users of Google Earth 7 to be able to click on the displayed photo and have their default JPG file type handler application open and display the photo file. I have tried enclosing the `<img>` tag with an `<a href>` tag, using the notation `href="file://photos\technician01\photo1.jpg"`, but I can't get the behaviour I am looking for (no external viewer is ever spawned). Is there a way to accomplish this? **Update** Here's a complete block of HTML I have attempted to achieve the described functionality:

     `


title: How to save a joined spatialite table when each table contains PKUID fields?

text: QGIS 1.8.0 -- QSpatialite 6.0.7 -- Windows XP I have spatialite tables joined together that I need to save as a single table, but an error-message says that QGIS can't save the table because there are two PKUID field columns. Spatialite_GUI can delete most spatialite table columns, when I need to do that to manually adjust sqlite column names, but my Qspatialite or QGIS always generates PKUID columns for every table it makes, and that column is always locked (I suppose as a primary key column) The PKUID column can't be deleted by Spatialite_GUI since it's locked, and so I get stuck and have to export everything to CSV files. That's the only workaround I know of, to save the table out as CSV file, manually delete the column from CSV and reload it, then join and "save as". Is there an easier way to save joined tables in QGIS (that have duplicate columns) or locked columns, rather than manually deleting or renaming identically named columns beforehand with CSV files. This is a chore since sqlite db aren't designed to have their field columns renamed or deleted, and QGIS doesn't automatically rename duplicated column names. What's the easiest way to save a joined table with duplicate column names in QGIS? Thanks in advance!


Does ArcGIS 10 have a Geoprocessor Programming Model by using Arcpy? is it possible to use ArcGIS 9.3 Geoprocessor Programming Model within ArcGIS 10?


title: Smallest GIS library for embedded devices?

text: I need to use very simple geometry operations (read SHP or other standard format and simple searches) within an embedded device running a simple linux based os. What is the smallest and simplest API I could use? It could be written in C/C++ or Python in order to allow the integration with my system. Thanks, Samuel


title: "Geometry type does not exist" on a Postgis installation on Fedora 18

text: Just installed PostgreSql and PostGIS over a brand new installation of Fedora 18. I did it this way:               sudo yum install postgresql postgresql-server postgrsql-contrib postgresql-devel pgadmin3      Once installed PostgreSql, I created the postgres user as usual. Then, I installed PostGIS:               sudo yum install proj-devel gdal postgis      As a last step, I've created a new user (mydbuser) and a new db owned by that user:               psql -U postgres -W     CREATE DATABASE mydb OWNER mydbuser;     \q     createlang plpgsql mydb      Now, I've tried to execute the following line to use PostGIS on mydb:               psql -d mydb -U mydbuser -W -f /usr/share/pgsql/contrib/postgis-64.sql      but I get a long list of errors of the same type, this is the first line:               psql:/usr/share/pgsql/contrib/postgis-64.sql:7751: ERRORE:  il tipo "geometry" non esiste      In English, the message is 'ERROR: "geometry" type does not exist'. I've also tried to execute the code as the postgres user, but nothing changed. Any ideas? Diego


title: When should I be calling a 'rectangular grid of pixels' a raster or an image?

text: Reading through GIS literature, it sometimes feels that authors use the two terms with a specific meaning in mind. I always treat the terms interchangeably. Am I wrong? Read the first line in the link. See how it says an image service provides access to raster (and image), hinting at different meanings for the two terms.


title: How to get started on web GIS development with Java?

text: I am very new for GIS, I know java web development. I want to develop a GIS web application in Struts, can you suggest me how to start? My requirement are as -   1. Any API/Tool should be open source.   2. It should work on HTTPS. Any help or tutorial will be appreciate. Thanks in advance!


title: Is there any batch process to get a single band image from stacked image

text: Can you help me to extract a single band image from stacked multiband image. Totally, 260 image has been stacked. When I do to this, it will take a long time. So I'd like to know there is any batch process to extract single band image in ARCGIS, ERDAS ENVi Thanks all


title: Problem with NDVI image value

text: I'm really happy to find this website. It is very helpful for the learners of Arcgis. Now I wanna ask a question. When i open NDVI image in ENVI, its value is still within the range (-1 to 1). But when I open this NDVI image in ERDAS and Arcgis, its value is 0-255. What does it mean? Is the NDVI image some problem? best regards


title: Intersecting two shapefiles from different folders based on their name using scripting in ArcGIS Desktop?

text: How can i script for intersecting the two shape files from two different folders based on their names? i have two folders which contains the shape files with same name i want to intersect the shapefiles based on their name, scripting should pick the files based on their name from two different folders and output file should store in another folder. I'm using ArcGIS Desktop 9.3 and arcgisscripting (that preceded ArcPy).


title: How to display an arrow pointing to a specific bearing in QGIS

text: I have a number of point features which contain an attribute I call "bearing" and another one called "slope". Each point has a bearing between 0 to 359 degrees. And the slope is either "+15" or "-15", etc, such as +15% or -15% depending on the steepest slope of the ground. What I am wanting to do is display an arrow on the map with the arrow pointing in the bearing direction. With a negative slope% pointing outward from the point and a positive slope% pointing inward to the point. I'm using this in a forestry context with cruise plots. I'm fairly certain that ArcGIS can do this (As I've seen it on maps) but I'm hoping that QGIS can as well.


title: Displaying Feature Info for a user defined polygon

text: To retrieve feature info for a polygon selection, I have the OpenLayers code below:                       var featureTypes  = ["nig_admin1", "nig_admin2", "nig_admin3", "nig_landscan"];             var populationCtrl = new OpenLayers.Control.GetFeature({                 protocol: OpenLayers.Protocol.WFS.fromWMSLayer(nig_landscan, {                     url: "http://localhost:8081/geoserver/wfs",                     featureType: featureTypes                 }),                 hover: false             });                  populationCtrl.events.register("featureselected", this, function(e){                 select.addFeatures([e.feature]);             });      I'd like to view the features retrieved in a pop up, or something similar. I've tried adding a pop up code within the `featureselected` event, but I can't get anything to show. Any pointers on how I can proceed will be highly appreciated. `nig_landscan` in the code above refers to the corresponding WMS layer, which is a layer group in GeoServer. My environment setup: GeoServer 2.1.3, PostGIS 1.5, OpenLayers 2.11.


title: How to export highlighted features from attribute table to an excel file?

text: When i selected some portion of the map in QGIS, in the attribute table the selected features gets highlighted... How to export the highlighted features from attribute table to an excel file automatically..?


title: Ogr2ogr driver error when convert shapefile to mif

text: I use the following script to convert shapefile to mif:               ogr2ogr -f "MapInfo File" D:\docs\arcgis_data\shp_mif\BOU.mif D:\docs\arcgis_data\shp\BOU.shp      But I got the error:               ERROR 6: Unable to open D:\docs\arcgis_data\shp_mif\BOU.mif      What's the problem?


title: Geometry conversion problem.Data loss

text: I was trying to convert a roadnetwork table.Its geometry contains both MULTILINESTRING and LINESTRING.I converted all the MULTILINE string to Linestring using the query update all_india_merge set the_geom=ST_GeometryN(the_geom,1) where astext(the_geom) <> 'LINESTRING'; After conversion Some data has been lost.How can I recover those geometries. Please Help me.