title: Hover on Airport Data using OpenLayers + Mapfile

text:
I need to create a simple hover using openlayers for worldwide airport data. I'm thinking maybe the associated address, maybe a picture in regards to the airport, ...anything really. It's just an html file with some js, nothing too crazy but I've been out of the GIS/coding game for a while. This is the current code for the html file. As you can see, there is an associated map file. Since there are so many airports, I'm guessing you would want to put the code in the mapfile (in regards to the hover(s))

                  Here is the map file. Apologies in advance for a little screwed up indentation. So I guess my main question is: do I need to be messing with the map file or just the html? I found another example where they were referencing the hover in the map file but that example was not also a straight html file...if that makes sense. I looked through the openlayers documentation but it was pretty tough to grasp. Any input is greatly appreciated.


title: Create spatial view that references multiple databases

text: Based on this link, it sounds like I can create a view from separate spatial databases. > Cross-database queries are permitted and ArcSDE views can reference data > from multiple databases provided the CROSS_DB_QUERY_FILTER parameter in the > SDE_dbtune table is set to 0. But they don't show how. Does anyone know how I can go about doing this? I'm using MSSQL and Oracle databases. ![enter image description here](http://i.stack.imgur.com/rVsES.png)


title: Add raster paths to shapefile with python (ArcGIS)

text: I have a problem and maybe someone is able to help me or give me some hints how to start. Situation: I have a lot of raster and of every raster I have the center coordinate. Out of this coordinate I build a shapefile to have the information about the approximate location of each raster. Now I included the path for each raster in the attribute table of the shapefile with the field calculator. Problem: Some of the raster are exactly georeferenced with a word-file. I want to add a new attribute field to show if there is a georeference or not. For this I can use the following code (at the moment a input cursor is missing):               from arcpy import env          arcpy.env.overwriteOutput = 1          env.workspace = r'Z:\path'     path = env.workspace          # Create a table and add fields     table_name = "TempTable.dbf"     arcpy.CreateTable_management(path, table_name)     arcpy.AddField_management(path + "\\" + table_name, "RasterName", "TEXT")     arcpy.AddField_management(path + "\\" + table_name, "cs", "TEXT")          # Get the list of raster files     rasters = arcpy.ListRasters()          # Loop through the files and perform the processing     for ras in rasters:              # Determine if the input has a defined coordinate system         dsc = arcpy.GetRasterProperties_management(ras, "TOP")                       if float(str(dsc.getOutput(0)).replace(",", ".")) > 100:                  row.cs = "y"     # Adds "x" to the row where the image has a georeference              else:                  row.cs = "n"     # Adds "n" to the row where the image has no georeference      Now I do not know how to add the information in the right way to the attribute table (maybe a join?). I hope there is someone who is able to help me with this. Thanks!


title: Do transfer edges add to the distance cost of a route?

text: Transfer edges are commonly used in a network dataset to carry a time cost in one direction only (either to-from or from-to). This is fine if the impedance used is time, i.e. find the route that will take the shortest time. However, if the impedance is now set to distance (find the shortest route), will the length of the transfer edges unfairly impact the length of the route?


title: prevent popup from appearing when clicking outside wms layer

text: Any help with this would be appreciated. I am projecting a WMS-layer on a map using OpenLayers and Geoserver. The next thing I want to do is to display info on any given feature in my WMS-layer in a popup whenever a user clicks it. This is all working fine, but my problem is that whenever I click outside the WMS-layer, I get an empty popup. However, I only want a popup to appear when a user clicks an actual feature. Does anyone know how to fix this? Thnx in advance!


title: Why is pgrouting driving distance skipping edges?

text: I have a road network where I do a driving distance operation through pgrouting. This is my SQL:               SELECT * FROM driving_distance('SELECT gid AS id,source,target, cost FROM elveg',676663,1000,false,false);      From this query I get a strange result: http://bildr.no/view/1373984 Black is my road network and green is the result from the above query. From the image you see that some edges are not included(not green) while edges which can only be connected through this not included edge in fact is included? How can this be? Am I doing something totally wrong here?


title: What are the requirements for building a Web Mapping and analysis Site?

text: I want to make a site like bhuvan in which it can provide different layers of information and also tools for performing GIS analysis I know three main requirements :   * client side scripting (openLayers)   * Webservices (Mapserver) + Tiling/Caching Services   * Database (PostGIS) I have ArcGIS desktop installed on my computer...and we can digitise many of the layers. My questions are :   1. How can I put these layers I created in ArcGIS in my website i.e., how to publish these maps to any open webservices.   2. For what purpose we need Database like postGIS. How are webservices differ from databases. With out database, can we build a webGIS? Please correct me if I am wrong anywhere...Thank you!!!


title: Semi-Automated Classification Plugin ROI Creation Failure

text: I can't get the Region-growing to work when running the ROI creation module. What am I missing? Please see the "AttributeError" line halfway down: # An error has occured while executing Python code: Traceback (most recent call last): File "C:/Users/reidi/.qgis2/python/plugins\SemiAutomaticClassificationPlugin\semiautomaticclassificationplugin.py", line 1029, in pointerClickROI self.createROI(self.pntROI) File "C:/Users/reidi/.qgis2/python/plugins\SemiAutomaticClassificationPlugin\semiautomaticclassificationplugin.py", line 1245, in createROI self.regionGrowing(dictBands, point.x(), point.y(), self.rngRad, self.minROISz, tmpShp) File "C:/Users/reidi/.qgis2/python/plugins\SemiAutomaticClassificationPlugin\semiautomaticclassificationplugin.py", line 1323, in regionGrowing dtSrc = drvr.CreateDataSource(outputVector) AttributeError: 'NoneType' object has no attribute 'CreateDataSource' Python version: 2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)] QGIS version: 2.2.0-Valmiera Valmiera, c3a2817 Python path: ['C:/Users/reidi/.qgis2/python/plugins\processing', 'C:/OSGEO4~1/apps/qgis/./python', u'C:/Users/reidi/.qgis2/python', u'C:/Users/reidi/.qgis2/python/plugins', 'C:/OSGEO4~1/apps/qgis/./python/plugins', 'C:\OSGEO4~1\bin\python27.zip', 'C:\OSGEO4~1\apps\Python27\DLLs', 'C:\OSGEO4~1\apps\Python27\lib', 'C:\OSGEO4~1\apps\Python27\lib\plat-win', 'C:\OSGEO4~1\apps\Python27\lib\lib- tk', 'C:\OSGEO4~1\bin', 'C:\OSGEO4~1\apps\Python27', 'C:\OSGEO4~1\apps\Python27\lib\site-packages', 'C:\OSGEO4~1\apps\Python27\lib\site-packages\PIL', 'C:\OSGEO4~1\apps\Python27\lib\site-packages\win32', 'C:\OSGEO4~1\apps\Python27\lib\site-packages\win32\lib', 'C:\OSGEO4~1\apps\Python27\lib\site-packages\Pythonwin', 'C:\OSGEO4~1\apps\Python27\lib\site-packages\Shapely-1.2.18-py2.7-win- amd64.egg', 'C:\OSGEO4~1\apps\Python27\lib\site-packages\six-1.3.0-py2.7.egg', 'C:\OSGEO4~1\apps\Python27\lib\site-packages\wx-2.8-msw-unicode', 'C:\OSGEO4~1\apps\Python27\lib\site-packages\xlrd-0.9.2-py2.7.egg', 'C:\OSGEO4~1\apps\Python27\lib\site-packages\xlwt-0.7.5-py2.7.egg', 'C:\OSGEO4~1\apps\qgis\python\plugins\fTools\tools', 'C:/Users/reidi/Downloads/SACP'] #


title: geoserver - publishing on a public site

text: I have set up geoserver and have it running well on my LAN. I finally set up a domain and want to host my own geoserver instance publicly. I have a jetty version of geoserver running on an OS X, Mavericks Server. What do I have to do to point my geoserver instance publicly? Any tips or nudges in the general direction will be greatly appreciated. Thanks


title: shape file into DataSet in postGIS using openJUMP, but getting error

text: ![enter image description here](http://i.stack.imgur.com/wXVvz.jpg) I am new to this field and I want to save my shape file into DataSet in postGIS using openJUMP, but I'm getting an error...can anybody tell me solution for this?


title: which map projection should be used for global map with geography and human demographics

text: I am working on a project to try to identify the most threatened populations (on a global scale) of having a limited access to freshwater. I want to compare those most at risk currently and then add in climate change projections to see how it is expected to change. I am trying to justify what projection system to use. With such a global scale what is the best option, if focusing on freshwater geography as well as human demographics? I am a little overwhelmed, but looking at mostly WGS 1984 coordinate systems. How would you justify a certain projection for such a broad analysis area and concentrations.


title: Empty FeatureCollection using TinyOWS and OpenLayers

text: I'm using TinyOWS and OpenLayers following this sample but when I open the page I can't see anything. I suppose that is becouse I get an empty FeatureCollection from the response server (`localhost/cgi- bin/proxy.cgi?url=http%3A%2F%2F127.0.0.1%2Fcgi-bin%2Ftinyows`) as shows below:               <?xml version='1.0' encoding='UTF-8'?>     <wfs:FeatureCollection      xmlns:tows='http://127.0.0.1'      xmlns:wfs='http://www.opengis.net/wfs'      xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance'      xmlns:gml='http://www.opengis.net/gml'      xmlns:xsd='http://www.w3.org/2001/XMLSchema'      xmlns:ogc='http://www.opengis.net/ogc'      xmlns:xlink='http://www.w3.org/1999/xlink'      xmlns:ows='http://www.opengis.net/ows'      xsi:schemaLocation='http://127.0.0.1         http://127.0.0.1/cgi-bin/tinyows?service=WFS&version=1.1.0&request=DescribeFeatureType&Typename=tows:tows:fran   http://www.opengis.net/wfs        http://schemas.opengis.net/wfs/1.1.0/wfs.xsd        http://www.opengis.net/gml        http://schemas.opengis.net/gml/3.1.1/base/gml.xsd'     >     </wfs:FeatureCollection>      The openlayers script is:               var map, wfs;     OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";     var DeleteFeature = OpenLayers.Class(OpenLayers.Control, {       initialize: function(layer, options) {           OpenLayers.Control.prototype.initialize.apply(this, [options]);           this.layer = layer;           this.handler = new OpenLayers.Handler.Feature(               this, layer, {click: this.clickFeature}           );       },       clickFeature: function(feature) {           // if feature doesn't have a fid, destroy it           if(feature.fid == undefined) {               this.layer.destroyFeatures([feature]);           } else {               feature.state = OpenLayers.State.DELETE;               this.layer.events.triggerEvent("afterfeaturemodified",                                              {feature: feature});               feature.renderIntent = "select";               this.layer.drawFeature(feature);           }       },       setMap: function(map) {           this.handler.setMap(map);           OpenLayers.Control.prototype.setMap.apply(this, arguments);       },       CLASS_NAME: "OpenLayers.Control.DeleteFeature"     });     function showMsg(szMessage) {       document.getElementById("message").innerHTML = szMessage;       setTimeout(           "document.getElementById('message').innerHTML = ''",2000);     }     function showSuccessMsg(){       showMsg("Transaction successfully completed");     };     function showFailureMsg(){       showMsg("An error occured while operating the transaction");     };     function init() {       map = new OpenLayers.Map('map', {           projection: new OpenLayers.Projection("EPSG:31467"),           units: "m",           maxResolution: "auto",           maxExtent: new OpenLayers.Bounds(3427000,5788000,3444000,5800000),           controls: [               new OpenLayers.Control.PanZoom()           ]       });       var osm = new OpenLayers.Layer.WMS(           "OSM by Omniscale WMS",           "http://osm.omniscale.net/proxy/service",           {layers: 'osm', format: 'image/jpeg'},           {projection:"EPSG:31467",  units: "m", maxResolution: "auto", maxExtent: new OpenLayers.Bounds(3427000,5788000,3444000,5800000)}       );       var saveStrategy = new OpenLayers.Strategy.Save();       saveStrategy.events.register("success", '', showSuccessMsg);       saveStrategy.events.register("fail", '', showFailureMsg);       wfs = new OpenLayers.Layer.Vector("Editable Features", {           strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],           projection: new OpenLayers.Projection("EPSG:31467"),           protocol: new OpenLayers.Protocol.WFS({               version: "1.1.0",               srsName: "EPSG:31467",               url: "http://127.0.0.1/cgi-bin/tinyows",               featureNS :  "http://127.0.0.1",               featureType: "fran",               geometryName: "the_geom",               schema: "http://127.0.0.1/cgi-bin/tinyows?service=wfs&request=DescribeFeatureType&version=1.1.0&typename=tows:fran"           })       });       map.addLayers([osm, wfs]);       var panel = new OpenLayers.Control.Panel(           {'displayClass': 'customEditingToolbar'}       );       var navigate = new OpenLayers.Control.Navigation({           title: "Pan Map"       });       var draw = new OpenLayers.Control.DrawFeature(           wfs, OpenLayers.Handler.Polygon,           {               title: "Draw Feature",               displayClass: "olControlDrawFeaturePolygon",               multi: true           }       );       var edit = new OpenLayers.Control.ModifyFeature(wfs, {           title: "Modify Feature",           displayClass: "olControlModifyFeature"       });       var del = new DeleteFeature(wfs, {title: "Delete Feature"});       var save = new OpenLayers.Control.Button({           title: "Save Changes",           trigger: function() {               if(edit.feature) {                   edit.selectControl.unselectAll();               }               saveStrategy.save();           },           displayClass: "olControlSaveFeatures"       });       panel.addControls([navigate, save, del, edit, draw]);       panel.defaultControl = navigate;       map.addControl(panel);       map.zoomToMaxExtent();     }      This is the tinyOWS config file:               <tinyows online_resource="http://127.0.0.1/cgi-bin/tinyows"          schema_dir="/usr/local/share/tinyows/schema/">            <pg host="127.0.0.1" user="postgres" password="pass1234" dbname="test_gis" port="5432"/>            <metadata name="TinyOWS Server"             title="TinyOWS Server - Demo Service" />            <layer retrievable="1"          writable="1"          ns_prefix="tows"          ns_uri="http://127.0.0.1"              name="fran"              title="French Administrative Sub Boundaries (IGN - GeoFLA Departements)" />     </tinyows>      What am I doing wrong?


title: mile distance around each feature ArcGIS

text: I have a school layer information and a route map layer. In ArcMap, how can I show the x mile distance around each school base on route information. I did some research and I found different approaches such as Network Analysis tool, Proximity Tool, Buffer Geo-processing tool,... . I am not sure what approach I need to use.


title: SImple model works, resulting script doesn't

text: The goal is to convert an XY layer (from Excel) into a feature class. The "MakeXY" portion works flawlessly, the arcpy.FeatureClassToFeatureClass_conversion portion inevitably returns a 999999 "Table was not found" message. I can build a model and run it and it all works. The exported script, recycled through the Python window, also fails on a 999999. I figure it's got to be an acrpy.env.workspace issue (but may be wrong) and so have added that line to the exported script. Still get the 999999. So what's this newbie doing wrong?               import arcpy     arcpy.env.workspace = r"C:\\Temp\\MuleDeer\\AllKLM5.gdb"          # Local variables:     X12XY = "X12XY"      AllKLM5_gdb = "C:\\Temp\\MuleDeer\\AllKLM5.gdb"          # Process: Feature Class to Feature Class     arcpy.FeatureClassToFeatureClass_conversion(X12XY, AllKLM5_gdb, "CDM_9c_FC", "", "Freq \"Freq\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Freq,-1,-1;Date \"Date\" true true false 8 Date 0 0 ,First,#,SOB$_Features,Date,-1,-1;Time \"Time\" true true false 255 Text 0 0 ,First,#,SOB$_Features,Time,-1,-1;TTF \"TTF\" true true false 8 Double 6 15 ,First,#,SOB$_Features,TTF,-1,-1;Northing \"Northing\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Northing,-1,-1;Easting \"Easting\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Easting,-1,-1;Zone \"Zone\" true true false 255 Text 0 0 ,First,#,SOB$_Features,Zone,-1,-1;Altitude \"Altitude\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Altitude,-1,-1;Maxsnr \"Maxsnr\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Maxsnr,-1,-1;Hdop \"Hdop\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Hdop,-1,-1;Vdop \"Vdop\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Vdop,-1,-1;Satt \"Satt\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Satt,-1,-1;Fix \"Fix\" true true false 255 Text 0 0 ,First,#,SOB$_Features,Fix,-1,-1;V1 \"V1\" true true false 255 Text 0 0 ,First,#,SOB$_Features,V1,-1,-1;V2 \"V2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,V2,-1,-1;T1 \"T1\" true true false 8 Double 6 15 ,First,#,SOB$_Features,T1,-1,-1;T2 \"T2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,T2,-1,-1;H1 \"H1\" true true false 255 Text 0 0 ,First,#,SOB$_Features,H1,-1,-1;H2 \"H2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,H2,-1,-1;S1 \"S1\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S1,-1,-1;S2 \"S2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S2,-1,-1;S3 \"S3\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S3,-1,-1;S4 \"S4\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S4,-1,-1", "")      (The fuller script is a bit more elegant in its use of variables, etc. but fails in the same place with the same error.) = T = P.S. Get Spatial, here's the current script which fails as is. The larger script which creates the XY Layer is also available if you want me to post it but it's a pretty plain vanilla "arcpy.MakeXYEventLayer_management" call. I continue to wonder if the "X12XY" reference incomplete. I've also tried a string with contents like "C:\Temp\CDM\X12.xlsx" and "C:\Temp\CDM\x12.xlsx\X12XY" with no success.               arcpy.env.workspace = r"C:\Temp\MuleDeer"          # Local variables:     X12XY = r"X12XY"          AllKLM5_gdb = r"C:\Temp\MuleDeer\AllKLM5.gdb"          # Process: Feature Class to Feature Class     arcpy.FeatureClassToFeatureClass_conversion(X12XY, AllKLM5_gdb, "CDM_9c_FC", "", "Freq \"Freq\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Freq,-1,-1;Date \"Date\" true true false 8 Date 0 0 ,First,#,SOB$_Features,Date,-1,-1;Time \"Time\" true true false 255 Text 0 0 ,First,#,SOB$_Features,Time,-1,-1;TTF \"TTF\" true true false 8 Double 6 15 ,First,#,SOB$_Features,TTF,-1,-1;Northing \"Northing\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Northing,-1,-1;Easting \"Easting\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Easting,-1,-1;Zone \"Zone\" true true false 255 Text 0 0 ,First,#,SOB$_Features,Zone,-1,-1;Altitude \"Altitude\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Altitude,-1,-1;Maxsnr \"Maxsnr\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Maxsnr,-1,-1;Hdop \"Hdop\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Hdop,-1,-1;Vdop \"Vdop\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Vdop,-1,-1;Satt \"Satt\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Satt,-1,-1;Fix \"Fix\" true true false 255 Text 0 0 ,First,#,SOB$_Features,Fix,-1,-1;V1 \"V1\" true true false 255 Text 0 0 ,First,#,SOB$_Features,V1,-1,-1;V2 \"V2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,V2,-1,-1;T1 \"T1\" true true false 8 Double 6 15 ,First,#,SOB$_Features,T1,-1,-1;T2 \"T2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,T2,-1,-1;H1 \"H1\" true true false 255 Text 0 0 ,First,#,SOB$_Features,H1,-1,-1;H2 \"H2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,H2,-1,-1;S1 \"S1\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S1,-1,-1;S2 \"S2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S2,-1,-1;S3 \"S3\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S3,-1,-1;S4 \"S4\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S4,-1,-1", "")      P.P.S Here's the script creating the X12XY object. Works first time, every time. The eventual goal is to batch-process a set of (identical) Excel files into feature classes. The only workbook in the target directory is named XL12.xlsx; the worksheet of interest is labeled "SOB" (Note: The indentation is off as I had to tweak each line manually after pasting it here).               import arcpy, os, sys, traceback          arcpy.env.workspace = r"C:\Temp\CDM"     outPath = r"C:\TEMP\CDM"     arcpy.env.overwriteOutput = True          # Local variables     sWorksheetName = "SOB$"     sProjection = "PROJCS['NAD_1983_UTM_Zone_10N',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['false_easting',500000.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',-123.0],PARAMETER['scale_factor',0.9996],PARAMETER['latitude_of_origin',0.0],UNIT['Meter',1.0]];-5120900 -9998100 10000;-100000 10000;-100000 10000;0.001;0.001;0.001;IsHighPrecision"          # List all files in folder then list all tables in file     fileList = arcpy.ListFiles("*.xlsx")     print 'Listing .xlsx files in workspace'          try:          # Loop through each file     for fileT in fileList:         fileName = os.path.splitext(fileT)[0]         sDest = fileName  + "XY"         sSrc = fileT + "\\" + sWorksheetName         print fileT + " : " + fileName + " : " + sSrc           # Toolbox / Data mgt / Layers and Table vuews"          arcpy.MakeXYEventLayer_management(sSrc, "Easting", "Northing", sDest, sProjection, "")               except:          tb = sys.exc_info()[2]     tbinfo = traceback.format_tb(tb)[0]     pymsg = "PYTHON ERRORS:\nTraceback Info:\n" + tbinfo + "\nError Info:\n     " +        str(sys.exc_type) + ": " + str(sys.exc_value) + "\n"     msgs = "ARCPY ERRORS:\n" + arcpy.GetMessages(2) + "\n"          arcpy.AddError(msgs)     arcpy.AddError(pymsg)          print msgs     print pymsg          arcpy.AddMessage(arcpy.GetMessages(1))     print arcpy.GetMessages(1)      P.P.P.S Spatial guy, I get the concept of the MakeXY table existing only in memory. Less clear on the datatype issue (love non-typed languages - NOT!). The following code puts it all within the "for" loop yet still fails -- but perhaps for a new reason. The Python output is: Listing .xlsx files in workspace                   X12.xlsx : X12 : X12.xlsx\SOB$              Reality check: C:\Temp\MuleDeer\AllKLM5.gdb is there              X12XY is there               About to call FeatureClassToFeatureClass_conversion()              ARCPY ERRORS:         ERROR 999999: Error executing function.         The table was not found. [X12FC]         Failed to execute (FeatureClassToFeatureClass).                   PYTHON ERRORS:         Traceback Info:           File "<string>", line 42, in <module>              Error Info:              <class 'arcgisscripting.ExecuteError'>: ERROR 999999: Error executing function.         The table was not found. [X12FC]         Failed to execute (FeatureClassToFeatureClass).      I'm assuming (dangerous move) that if output feature class does not already exist the FeatureClassToFeatureClass() will create it. That's sort of the point, isn't it? So am I now battling a syntax error of some sort or just fundamentally unclear about FeatureClassToFeatureClass? Finally for the night, this system seems happy with the "str(sys.exc_value)" line. Don't ask, don't tell, just keeping my fingers crossed.               import arcpy, os, sys, traceback          arcpy.env.workspace = r"C:\Temp\CDM"     outPath = r"C:\TEMP\CDM"     arcpy.env.overwriteOutput = True          # Local variables          sXColName = "Easting"     sYColName = "Northing"     sWorksheetName = "SOB$"     sProjection = "PROJCS['NAD_1983_UTM_Zone_10N',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['false_easting',500000.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',-123.0],PARAMETER['scale_factor',0.9996],PARAMETER['latitude_of_origin',0.0],UNIT['Meter',1.0]];-5120900 -9998100 10000;-100000 10000;-100000 10000;0.001;0.001;0.001;IsHighPrecision"     sXLSMapping = "Freq \"Freq\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Freq,-1,-1;Date \"Date\" true true false 8 Date 0 0 ,First,#,SOB$_Features,Date,-1,-1;Time \"Time\" true true false 255 Text 0 0 ,First,#,SOB$_Features,Time,-1,-1;TTF \"TTF\" true true false 8 Double 6 15 ,First,#,SOB$_Features,TTF,-1,-1;Northing \"Northing\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Northing,-1,-1;Easting \"Easting\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Easting,-1,-1;Zone \"Zone\" true true false 255 Text 0 0 ,First,#,SOB$_Features,Zone,-1,-1;Altitude \"Altitude\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Altitude,-1,-1;Maxsnr \"Maxsnr\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Maxsnr,-1,-1;Hdop \"Hdop\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Hdop,-1,-1;Vdop \"Vdop\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Vdop,-1,-1;Satt \"Satt\" true true false 8 Double 6 15 ,First,#,SOB$_Features,Satt,-1,-1;Fix \"Fix\" true true false 255 Text 0 0 ,First,#,SOB$_Features,Fix,-1,-1;V1 \"V1\" true true false 255 Text 0 0 ,First,#,SOB$_Features,V1,-1,-1;V2 \"V2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,V2,-1,-1;T1 \"T1\" true true false 8 Double 6 15 ,First,#,SOB$_Features,T1,-1,-1;T2 \"T2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,T2,-1,-1;H1 \"H1\" true true false 255 Text 0 0 ,First,#,SOB$_Features,H1,-1,-1;H2 \"H2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,H2,-1,-1;S1 \"S1\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S1,-1,-1;S2 \"S2\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S2,-1,-1;S3 \"S3\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S3,-1,-1;S4 \"S4\" true true false 255 Text 0 0 ,First,#,SOB$_Features,S4,-1,-1"     sMasterGDB = "C:\\Temp\\MuleDeer\\AllKLM5.gdb"     # List all files in folder then list all tables in file     fileList = arcpy.ListFiles("*.xlsx")     print 'Listing .xlsx files in workspace'          try:              # Loop through each file         for fileT in fileList:             sFileName = os.path.splitext(fileT)[0]             sDestXY = sFileName + "XY"             sDestFC = sFileName + "FC"             sSrc = fileT + "\\" + sWorksheetName             print fileT + " : " + sFileName + " : " + sSrc                  # out with the old (if any)             if arcpy.Exists(sDestXY):                 arcpy.Delete_management(sDestXY)                  # out with the old (if any)             if arcpy.Exists(sDestFC):                 arcpy.Delete_management(sDestFC)                  if arcpy.Exists(sMasterGDB):                 print 'Reality check: ' + sMasterGDB + " is there"                  # Toolbox / Data mgt / Layers and Table views"              arcpy.MakeXYEventLayer_management(sSrc, sXColName, sYColName, sDestXY, sProjection, "")                  # Reality check             if arcpy.Exists(sDestXY):                 print sDestXY + " is there "                       print "About to call FeatureClassToFeatureClass_conversion()"             # Toolbox / Conversion / To Geodatabase             arcpy.FeatureClassToFeatureClass_conversion(sDestXY, sMasterGDB, sDestFC, "", sXLSMapping, "")                  del sXColName, sYColName, sWorksheetName, sProjection, sXLSMapping, sMasterGDB             del fileT, sFileName,  sDestXY, sDestFC, sSrc          except:              tb = sys.exc_info()[2]         tbinfo = traceback.format_tb(tb)[0]         pymsg = "PYTHON ERRORS:\nTraceback Info:\n" + tbinfo + "\nError Info:\n     " +        str(sys.exc_type) + ": " + str(sys.exc_value) + "\n"         msgs = "ARCPY ERRORS:\n" + arcpy.GetMessages(2) + "\n"              arcpy.AddError(msgs)         arcpy.AddError(pymsg)              print msgs         print pymsg              arcpy.AddMessage(arcpy.GetMessages(1))         print arcpy.GetMessages(1)


title: Crop collar from chart

text: I need advice on how to complete the following task: I have aeronautical charts covering parts of the US (example: http://aeronav.faa.gov/content/aeronav/sectional_files/Albuquerque_89.zip - Each zip contains ~2 files covering both the northern and the southern part of the area). Unfortunately for me, these maps have collars that I want to remove so that I could create a new geotiff file with only the real content of the original geotiff. I have tried Global Mapper's "Crop collar" command, but that didnt gave me the results I wanted to have. In my understanding I should use a tool like GDAL warp to complete this task in bulk for all of the images (~106 images in total), but the situation is that I am completely inexperienced using that tool. I have already tried running `gdalwarp -srcnodata 0 -dstalpha input.tif output.tif` but that didn't gave me good results as well, therefore I would be glad to receive some help! Thanks in advance


title: GDal problem "no suitable image found" when I run “ogr2ogr” command

text: I am trying to run "ogr2ogr" command               macbook$ sudo ogr2ogr -f "PostgreSQL" PG:"dbname=test user=postgres" /Users/macbook/planning.gdb      But I am getting this error:               dyld: Library not loaded: /Library/Frameworks/GDAL.framework/Versions/1.10/GDAL       Referenced from: /Library/Frameworks/GDAL.framework/Programs/ogr2ogr       Reason: no suitable image found.  Did find:         /usr/local/lib/GDAL: not a file      I following this answer : Error Updating Gdal with Homebrew               macbook$ export PATH=/Library/Frameworks/GDAL.framework/Programs:$PATH     macbook$ which ogr2ogr      /Library/Frameworks/GDAL.framework/Programs/ogr2ogr          But I am getting this error whenever I run the "ogr2ogr" command           macbook$ sudo ogr2ogr -f "PostgreSQL" PG:"dbname=test user=postgres" /Users/macbook/planning.gdb     dyld: Library not loaded: /Library/Frameworks/GDAL.framework/Versions/1.10/GDAL       Referenced from: /Library/Frameworks/GDAL.framework/Programs/ogr2ogr       Reason: no suitable image found.  Did find:         /usr/local/lib/GDAL: not a file      When I run : macbook$ ls /usr/local/lib/ I can find the folder :               dtrace                  libopencv_highgui.dylib     gcc                 libopencv_imgproc.2.4.6.dylib     gdal                    libopencv_imgproc.2.4.dylib


title: Is there a significantly faster way to buffer than arcmap?

text: I hope this question is appropriate here and I apologize if it is more of a discussion question than a specific technical problem. I have to buffer 150,000 polygons for a real estate project. While it's not impossibly slow, it got me wondering about 'the fastest way to geoprocess'. What packages, tips, or tricks could be used to make this faster? Does anyone have an idea how arcmap, arcpy, shapely, postgres compare in performance time? EDIT: After asking around, I heard 'arcpy is slow!' a lot.


title: How to create a point layer from a one-to-many join?

text: I have two layers:   * One with geographical coordinates containing every address in the county and a unique identifier.    * The second contains a list of incidents with the unique identifier from the address layer.  I could join the incidents layer to the address layer and then query only those addresses that have an incident. However there could be multiple incidents at on address, so I need to join the address layer to the incidents table. This gives me a record for each incident and the X Y coordinates from the address layer. How do I specify that this layer now contains points without writing an intermediate file and re importing? The two original layers will change with time and I don't want to save intermediate files.


title: Examples of Web Maps

text: As part of an internal "learn at work month" program, myself and a colleague have volunteered to do a presentation: **"The Wonderful World of Web Maps"**. To this end, we're now looking for examples of wonderful, weird, functional, beautiful, and just plain awesome web maps that are out there. Any online map that is either interactive or dynamically updates itself would qualify. Part of this question is already answered with the Examples of Beautiful Maps question which has some-web maps as does Examples of badly designed maps but we're more interested in a diverse range of maps rather than just the beautiful and bad. I've done a little searching and come up with a couple (that just happen to have a transport theme) like: The Live Ship Map Flight Tracker (Heathrow Airport) What other web-maps are out there that might warrant demonstrating or could use wider exposure?


title: Create Geotiff from Tiff using java

text: Respected sir/madam, I want to add geo tags to existing TIFF file and make it to geotiff file. Otherwise, for bufferedimage is it possible to add geotags to make Geotiff file. Thanks in advance...


title: Is is possible to have a sld with a rotated point symbolizer along a line in geoserver?

text: Geoserver 2.1.3 I would like to symbolize a line with a Square but use the `<Rotation>45</Rotation>` tags to make the symbol look like a diamond. My sld is below and unfortunately the symbols are not rotated...Is this not possible or is my SLD wrong? thanks, Rob               <?xml version="1.0" encoding="ISO-8859-1"?>     <StyledLayerDescriptor version="1.0.0"     xsi:schemaLocation="http://www.opengis.net/sld StyledLayerDescriptor.xsd"     xmlns="http://www.opengis.net/sld"     xmlns:ogc="http://www.opengis.net/ogc"     xmlns:xlink="http://www.w3.org/1999/xlink"     xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">       <NamedLayer>         <Name>line-symbol rotation</Name>         <UserStyle>           <Title>line-symbol rotation</Title>                 <FeatureTypeStyle>                     <Rule>                     <Name>line-symbol rotation</Name>                                         <LineSymbolizer >                     <Stroke>                   <GraphicStroke>                     <Graphic>                <Mark>                  <WellKnownName>square</WellKnownName>                  <Fill>                    <CssParameter name="fill">#009900</CssParameter>                  </Fill>                </Mark>                <Size>12</Size>                <Rotation>45</Rotation>              </Graphic>                              </GraphicStroke>                               </Stroke>               </LineSymbolizer>                              </Rule>                   </FeatureTypeStyle>                  </UserStyle>       </NamedLayer>     </StyledLayerDescriptor>


title: Is there an easy common gis abstraction available in .Net

text: Is there an easy common gis abstraction available in the .Net programming world that cuts the dependency to a concrete gis implementation, e.g. to the esri silverlight/wpf sdk. ... for example OpenLayers is such a library for several map providers in JavaScript. By using an IoC container (like Spring.Net) the dependency could be hidden. I noticed the GeoAPI buth this seems too complex for an easy abstraction level, I think.


title: add wkt representation of the geometry as field in qgis

text: I have the shapefile of boundaries loaded in qgis. I want to add a field in the file which is WKT representation of the feature. **Update** As Nathan W answered I used `geomToWKT( $geometry )` and got result. The thing is the coordinates are in lon lat format. I found out that that i need them in lat lon format. Is it possible to change the order of coordinates?


title: How to find locations for hydro power using raster analysis?

text: I have a stream network and DEM file both in raster format. I need to extract values from DEM to stream network. For the analysis in stream network after that, I need to find out the places for hydro power based on various criteria of head, stream order. So how can i start my work? I am planning to convert both raster file into text file using numpy array and extracting head for all the cells of stream and after that i will start processing. Is that a good idea to start or is there any other way? I am quite new to python.


title: How to extend a custom NASolver in ArcGIS Server 10.1?

text: The past year we developed a custom NASolver with C# , it works fine and that job has finished. Now we are planning to use it from ArcGIS Server. The problem is that we have no information for that. The creation of a new Server Object Extensions could be a solution, but it will take time... We want to know if exists another way to do that by extending the iNAServer object. Thanks in advance for your help.


title: Comparing census tracts

text: I am working with Arcmap 10.1 and census tract level data for 1990 and 2011. I am trying to compare different options for tracking changes between the tracts. The tracts change between the census polls, and there are some tracts from 1990 that merged to form a single tract in 2011, and those that split up. I would like to know what are some options I have to keep track of the changes in the tracts and how I can compare demographic changes from 1990 to 2011?


title: assign a number for each group in field calculator

text: I have several thousand features each of which has a name. The names are not unique. what I want to do with the field calculator is assign a number to each name so I can refer to a number and always get the same name. For example if I had:               chris, dave, chris, tom, dave, chris , dave, tom      I'd like my number column to show               1, 2, 1, 3, 2, 1, 2, 3      Is this something I can do using python in field calucator?


title: delete duplicate

text: Delete duplicate. Delete duplicate.


title: Problem with georeferencer in Qgis 2 64bit installation

text: I tried to use the georeferencer in the new Qgis 2 64 bits for windows. I tried with the standallone installer and with the OSgeo4w installer. The problem is when I load the raster that I want to georeference, all I get is a gray square with a bunch of blurred colors at the top. I tried to use the georeferencer in Vista 64 (8gb of RAM) and in win 7 64 (16 gb of RAM). Also running Qgis as an administrator did not solve the problem. I used the same image in a 32 bit installation of Qgis 2 and I have not a problem loading it or working with it. I would appreciate any hint in what can I do to solve the issue


title: How to update to QGIS 2.0.1 on Mac?

text: I'm trying to upgrade to the newest version of qGIS but get the following error message when I open the installer:               QGIS requires the GEOS framework version 3.4 or newer.      I've had no luck maintaining my version of geos via the geos download page, and have instead installed it via homebrew. Is it possible to tell qGIS where the homebrew maintained version of GEOS lives? I would maintain osgeo via http://trac.osgeo.org/geos/ but keep getting error messages when I run the newest installer (3.4.2) and all advice I've found elsewhere on debugging the installer suggest that I shouldn't bother messing around w/ all the dependencies and should instead use homebrew. So this brings me back to where I started: is it possible to have homebrew maintain a version of osgeo that qgis can continue to locate? Or, am I forced to keep trying to debug the geos installer?


title: Need help in Symbology for ArcGIS JS API 3.2/3.3

text: I am displaying various symbols in my application and initially all are displaying at once (looks like bulky) so how do I set symbol levels in ArcGIS so that according to zoom level symbol size will increase or it will auto-fit according to resolution. Please check the screenshot for various levels. ![Level 1](http://i.stack.imgur.com/uzrTr.jpg) ![Level 2](http://i.stack.imgur.com/8NO5n.png) ![Level 3](http://i.stack.imgur.com/0miRJ.png) I am researching on my side and ArcGIS help. (I think visible scale range will work but I am not sure and total layers are five so should I set same levels for all layers? **Update :** I can not change the default extent (because that layers data is coming dynamically (We have data across the world so setting to specific extent is not possible)) Software : ArcGIS Server/Desktop 10.1 , ArcGIS JSAPI 32./3.3 Any help will be great!! ! Thanks in advance : )


title: Calculate geographic size in meters given lower left coordinate and upper right coordinate

text: How do you find the geographic size in meters given a lower left coordinate and upper right coordinate. The projection would be WGS84 projection. Ex. 33,-119,35,-117. What would be the area in meters that covers?


title: QGIS 2.2 - possible to do a JOINS from the Action list?

text: I have a set of layers that need to be joined ( _Layer Properties_ > _Joins_ ) to one specific layer. The reason for this is that I will be updating and adding new layers into my project and would like a single **Action** command for the specific layer to automatically have the selected layers join with it. I want to avoid having to run joining/merging algorithms altogether. Is this easily achievable? After spending time searching on how to create Actions, most of what I found describe the same command which is how to open up a web browser to search for whatever is in the specified attribute. Unfortunately, the default actions don't seem to contain anything similar to what I am looking for.


title: Global timer java script in KML

text: i have a timer script which shows the time of system in Description tag               <script type='text/javascript'>       setInterval('tick()',50);   //tick is a plugin and 50 is in ms delay time      function tick()      {        document.getElementById('clock').innerHTML=new Date();      }     </script>      i want that if one placemark is placed in India , then same script should show the current time of India and if one placemark is in America , it should show America Time only regardless of machine in which KML has been loaded. How can i do this ?? how can i sync all time


title: fTools "Join Attributes by Location" missing in SEXTANTE/Processing Toolbox

text: I am working with the new Processing Framework in QGIS 2.0.1. I really like the graphical modeler, because it allows my to make my own tools without having to be able to write pyhton or any other programming language. But it seems like some of the fTools (the Vector-Tools in QGIS) are missing? For example "Join Attributes by Location". Is this an error, or is there a reason for that tool being missing in SEXTANTE? To make things a bit more clear, I attached 2 Screenshots to show which tool I mean. Here you see all the tools avaliable in the Processing Framework ![Here is the list of vector tools in the framework](http://i.stack.imgur.com/QHToE.png) Here you see the tool I am looking for.. ![Here is the tool in the vector toolbar in QGIS](http://i.stack.imgur.com/C4uPK.png) So is there any way to get the Tool in the Framework? If not, has anyone an idea how to get a similar result with another tool? (Futher information on the task i want to accomplish: I've got a .shp with landlot, and a .shp with 3 polygons, and I want to have a field in the attribute table of the landlots which shows which of the polygons is covering which landlot.) Thanks in adavance! Gunnar


title: How to select pixel with smallest number in stack raster using Python script in ArcGIS for Desktop?

text: I have 6 tiff file in ArcGIS which are identical in term of geometry (same frame, same number of row and column). The pixel value of each tiff file is integer and range from 1 to 4. So if I overlay them on top of each other, for each pixel I will have 6 DN number (think like an RGB color composite). Now I need an output raster with pixel value that take the smallest number in the stack. So for 2 pixel that have 123456 and 334662 it will take 1 and 2 as output respectively. I would like to do it in Python but if it can be done via Raster Calculator or ModelBuilder it would also be good.


title: Identify segments that share a last vertex with the first vertex of selected segment

text: This has been marked as a duplicate to Connecting line features and determining the length of the longest line but I believe there are some subtle differences. And to be honest I'm not an advanced enough GIS user to be able to fully understand the solution provided. Short version, using python for ArcGIS, what is the quickest method to identify all segments that share a last vertex with the first vertex of a selected segment. I also need to extract associated data of the identified segments from the dbf. Currently I'm using cursors and x/y co-ordinate data! It takes a while. Long version, I'm looking to calculate the maximum upstream length of all line segments in a shape file and add that information to the dbf. I have a number of shapefiles to do this for with between 5000 and 10,000 segments per file. See image below for an example of what I am trying to do. I have a python script that works but it takes a long time to run. The script is below. In summary, it starts of by identifying all the unique Order numbers. The Order numbers relate to how many segments are upstream of a segment. The script then sets all segments with an Order number of 1 to the length of the segment. These are the most upstream segments. The final step is to iterate over the remaining order numbers. For each segment (Si), the length is obtained and the coordinates of the first vertex are identified. The cumulative length of any segments (S1, S2 etc) that share a common last vertex are also identified and the maximum cumulative length obtained. The cumulative length for Si can then be assigned to the segment The part of the script that is taking up most the time is lines 47 to 49, here I look for any connected upstream segments using a search cursor. The problem is I loop through the entire dbf and recreate the cursor for almost every row!                           import arcpy                 from sys import argv                      #Add path to location where flow paths are stored.  Notes for the attribute table                 #The length of each segment must be located in a column called "Length"                 #The order of each segment must be located in a column called "Order",                 # Saga calculates this when it generates channel networks                 #The x and y co-ordinates of the first and last vertex must be calculated and                 # saved in columns "X_Cord_F", "Y_Cord_F", "X_Cord_L" and "Y_Cord_L" respectively                 # I used the field editor in QGiS to do this.                 #The order of each segment must be located in a column called "Order", Saga calculates this when it generates channel networks                 Flow_Paths = "C:\Users\chcgis\Desktop\WELLINGTON\WellingtonOLF\Results\Channel12500.shp"                      #Finds all unique Order numbers in the shape file                 rows = arcpy.UpdateCursor(Flow_Paths)                  output = set()                 for row in rows:                     if row.Order not in output:                         output.add(row.Order)                 output.remove(1)#remove '1' from list as it is a special case addressed below                 output = sorted(output)                 outputList = list(output)                 del row                 del rows                      #Initialises the first Order number with the length                 rows = arcpy.UpdateCursor(Flow_Paths)                  Length = 0                 for row in rows:                     if row.Order == 1:                         row.Cum_Length = row.Length                         rows.updateRow(row)                 del row                 del rows                      #Loops over remaining order numbers and finds the longest upstream length                 for a in outputList:                      print a #A check to see how far along we are                     rows = arcpy.UpdateCursor(Flow_Paths)                      for row in rows:                         Cum_Length = 0                         if row.Order == a: #Loop over current order number                             X_F = row.X_Cord_F # Firs X Cord                             Y_F = row.Y_Cord_F # Firs Y Cord                             Length = row.Length                             rows_i = arcpy.SearchCursor(Flow_Paths,fields="X_Cord_L; Y_Cord_L; Cum_Length")                              for row_i in rows_i: #Find upstream path lengths and extract maximum lengths                                 if X_F == row_i.X_Cord_L:                                     if Y_F == row_i.Y_Cord_L:                                                    Cum_Length = max(Cum_Length,row_i.Cum_Length+Length)                             row.Cum_Length = Cum_Length #Update row with max length                             rows.updateRow(row)                     del row                     del rows              ![enter image description here](http://i.stack.imgur.com/yXLFn.jpg) I also have access to Saga and QGiS. All help is much appreciated! Cheers,


title: Exporting data from Collect Values output in ArcGIS ModelBuilder?

text: I'm currently playing with ModelBuilder. I've added a picture of the model I have at the moment. This model is currently iterating 6 shapefiles, so each of the 'output values' contain 6 figures in a list. I can't find a way to extract the values from these lists into a table/text file or similar. Is there a way to do this? ![model](http://i.stack.imgur.com/G5Knv.jpg) This is what is shown when I open 'output values' after running the model: ![enter image description here](http://i.stack.imgur.com/wwGle.png) I just want to grab those 6 numbers somehow...


title: PgRouting routing algorithm from start to end via nodes with restrictions

text: I have the following situation: I want to find the shortest route from `A to B through nodes (C,D,E, etc.)` with turn directions and road type filtering. Is there a way to accomplish this using pgRouting? Has anyone done this? I haven't begun developing anything yet, I'm just trying to find out if pgRouting is the way to go. Thanks. Update: The algorithm should determine the order of nodes (C,D,E, etc.)


title: Cut many lines with many polygons: Too much for my CPU/RAM?

text: I have a line-shp (80mb) and a polygon-shp (5mb) and want to intersect them. In result, the lines (contours) should be cut by the polygons. After using one of the tools `"Vector -> Geoprocessing -> Dissolve/Clip/Intersect"` to do so, the loading bar keeps at 0%. Same by marking just 10 of the 6000 contours to intersect them. The CPU usage is at 50% (one core 100%) all the time, even for more then 30minutes... What is the problem? How to solve it? `But dont tell me, I need a new computer! :)` > My config: AMD II X2 240 ; 4GB RAM; xubuntu 12.04; Qgis 1.8


title: QGIS on Thin Client and Windows 7/8

text: We are migrating from windows XP in the coming months, potentially to a windows 7 VDI environment. Does anyone know answers to the following: • Is QGIS supported running in a thin client environment? • If the answer above is yes, which thin client solutions is your product supported within? (i.e. citrix, Microsoft, vmware) • Is QGIS capable of running in a 64 bit windows 7 / 8 environment? • A recommended client specification, physical and VDI. Cheers Thomas


title: GeoServer + ArcMap

text: I need to setup GeoServer and access GeoTIFF images in the ArcMap. I have started with installing Geoserver and created a layer with GeoTIFF. Now, what are the next steps to publish the map(using WMS) and access them in the ArcMap. Can anyone please help me in this?


title: load base layer from database

text: I have svg maps of floor plans which I am using OpenLayers to project and display interactive features on. I would like to load the files from my database, (where they are stored as blobs) instead of specifying a path in my file system. Is this possible?


title: (How) can I write string data to a raster with Python/GDAL?

text: In Python, I have a NumPy array of string values (categories of some variable) that I'd like to write to a raster file (.tif or .img) using GDAL. The impression I get is that it's not possible to write a string value into the grid cell of a raster file, meaning I'd need to create an integer lookup table and write the integers into the grid cells instead. Then I'd need to write the character data into the column of the raster's attribute table so that you could know what string values the integers represent. Is my reasoning correct? If so, I'd like to know the basic process for building a raster attribute table using Python, GDAL. I haven't had much success in the past, and, as I understand it, not all GDAL raster formats are able to have attribute tables at all anyway. Thankyou in advance!


title: to extract 3d points from a shp file

text: I'm trying to work with the shapefile on python, the file " donnees.shp" contains 3d points I need to extract attributes from this file but my code returns an error:               >>> import shapefile        >>> r = shapefile.Reader("donnees.shp")      >>> r.shapes()[0].shapeType      ... 1       >>> r.shapes()[0].points         ...[[841001.0, 6399999.0]]     >>> r.shapes()[0].z          ... Traceback (most recent call last):     File "", line 1, in      r.shapes()[0].z     AttributeError: _Shape instance has no attribute 'z'      what is the nature of this error? what am I missing? thank you for you answer


title: Joining Census Data

text: I have collected census data for a county in Pennsylvania. Data varies from 1999-2010. For my data before 2010, a significant number of blocks do not match up (image below) when I join the data (based on GEOID). I have joined to 1990, 2000 and 2010 census block groups. There are roughly the same number of block groups, but the GEOID's are always off for maybe 30% of the block groups.![enter image description here](http://i.stack.imgur.com/e9Oj3.jpg) I feel like i'm missing something here... Anyone had trouble with this before? So I'm offering a bounty now. Basically I just need to be able to join this data with the shapefile block groups from any year... But my guess is that they would best align with 2000.


title: Difference between lon and lng

text: When browsing the web and documentation I see the abbreviation lon and lng used for the longitude part of geometry. I was wondering if there is any difference between lng and lon.


title: KML with Google Maps : Google violates same origin policy?

text: Doesn’t Google Violate the same origin policy when it tries to access the KML files from a globally accessible server? The condition for Google maps to display the KML is that it should be globally accessible. Now their are 2 problems with this:   1. Security of the KML files. (This can be avoided by using geoXML3)   2. Violating the Same-Origin Policy. Is their a way around to the same-origin policy? How does google parse the KML file? Thankyou.


title: How to save table from PostgreSQL to MapInfo File?

text: I have two tables: 1 with geometry in WGS   2 with geometry in Local System I tried this:               ogr2ogr -f "MapInfo File" output2.tab PG:"host=localhost user=postgres password=123456 dbname=postgis" -sql "SELECT * from filedata WHERE num=1"      and get files. However, after I open the files in QGIS, I see that file from the table in WGS coordinates is correct but the file from the table in the local system has a shift. WGS ![enter image description here](http://i.stack.imgur.com/pPd0w.jpg) Local system ![enter image description here](http://i.stack.imgur.com/uiUqG.jpg) The upper part of the polygons are different but its ok. I add proj line to file `GDAL\projlib\epsg` but it did not help. What am I doing wrong?


title: While developing stage, Is Proxy.cgi required to get Openlayers WMSGetFeatureInfo?

text: I have published the data in Geoserver (Tomcat and OS Ubuntu) and I am developing a map application in my local machine (OS windows xp) using asp.net. In my application, I can get WMS feature information in IE8 with the confirmation msg "This page is accessing information that is not under its control. This poses a security risk. Do you want to continue?" given Yes. But in Chrome, i cannt get feature info. It says               XMLHttpRequest cannot load http://xxx.xxx.x.xxx:8080/geoserver/cite/wms?LAYERS=cite%3Aplym_p1_examplec…e%2Fpng&INFO_FORMAT=application%2Fvnd.ogc.gml&SRS=EPSG%3A27700&X=596&Y=160. Origin http://localhost:3890 is not allowed by Access-Control-Allow-Origin.      I have found similar questions Probleme with WFS protocol in OpenLayers but its for after hosting the web application in geoserver.My case is still its in under developing in my local machine. My code is:               roadidPicker = new OpenLayers.Control.WMSGetFeatureInfo({     url: 'http://xxx.xxx.x.xxx:8080/geoserver/cite/wms',         title: 'identify features on click',                 queryVisible: true     });     roadidPicker.infoFormat = 'application/vnd.ogc.gml';     roadidPicker.events.register("getfeatureinfo", this, pickRoadid);     osMap.addControl(roadidPicker);     roadidPicker.activate();          function pickRoadid(e) {         if (e.features.length > 0) {         var val = e.features[0].attributes['theme'];         alert(val);     }     }      I have configured the proxy.cgi file in Geoserver (usr/lip/cgi-bin/proxy.cgi) and added my local ip address in               allowedHosts = ['xxx.xxx.x.xxx', 'www.openlayers.org', 'openlayers.org', 'labs.metacarta.com', 'world.freemap.in', 'prototype.openmnnd.org', 'geo.openplans.org', 'sigma.openplans.org', 'demo.opengeo.org', 'www.openstreetmap.org', 'sample.azavea.com', 'v2.suite.opengeo.org', 'v-swe.uni-muenster.de:8080', 'vmap0.tiles.osgeo.org','www.openrouteservice.org']      and added the below line in my script.               OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";      but it says another error in chrome console:               GET localhost:3890/usr/lib/cgi-bin/….…T%3Dapplication%252Fvnd.ogc.gml%26SRS%3DEPSG%253A27700%26X%3D682%26Y%3D284 404 (Not Found)       I dont know whether its right or wrong. Proxy.cgi should be copied in Geoserver or my machine? in which path? Please anyone suggest me to fix this issue


title: How to detect if a point lies inside a building?

text: My main motive is to detect if there's in any infrastructure nearby to a given point in a certain radius. For this I am using Places api and using nearbySearch method of it. The problem in this is that it only detects any infrastructure that is already marked by someone on Google Map's. Eg: http://web.iiit.ac.in/~jayesh.lahori/map.png , blue point is itself inside the building, this will not be detected because nothing has been marked at this spot. Any suggestions ??


title: Error with Gdal 1.10 install on Ubuntu 13.10

text: I've installed GDAL 1.9 then upgrade to 1.10 on Ubuntu 13.10 and then now it's error (gdalwarp can work but gdal_translate is fail) gdal_translate ERROR 1: gdal_translate was compiled against GDAL 1.9 but current library version is 1.10 How can I fix this problem? Thanks


title: What is the best way to display point data (with standard deviation) using QGIS/SAGA/Manifold?

text: 1. I would like to display point data values (indicating biomass at a point) but also to show the standard deviation value for each point. What would be the best method to do this in QGIS (I'm a basic user!). I do not have ±100 points. I am imagining a point with size = biomass and a buffer showing SDev, but I don't know how to do this.   2. The next step would be to show the breakdown of values in each point by species (about 5 species biomass values making up the total biomass at each point). Would the best way be with a pie graph per point? And how to show each StDev? I am using QGIS, Manifold and (less often) SAGA. Any suggestions welcome, many thanks


title: How to set visibility of columns of the attributetable?

text: Is there a possibility in QGIS setting a field(column) of the attribute-table to not visible? In the properties-dialog of the layer i couldn't find something. Thanx for any idea.


title: mapserver not displaying anything

text: I've been reading through Tyler Mitchel's Web Mapping Illustrated and also downloaded a copy of the itasca sample from maptools.org. Then created some shapefiles in OpenEV and modified the itasca.map file to include those layers and removed all the other layers, but it seems that there's no effect - a blank display I've also tried to play with the outline color, but nothing. Find the map file below: find the map file here: http://kompisen.se/xqSjGadQsN Update: forgot to specify some details, sorry. OS: Ubuntu Mapserver: 5.6.1.1 OpenEV: 1.8 installed via FWTools 2.0.6


title: Import Error:No module named shapely.wkb

text: OpenERP GeoEngine I have installed **OpenERP 6.1** version, works fine .,When i'm trying to add **CamptoCamp GeoEngine Modules with OpenERP** getting below find eror: **Installation Prerequisites:**               OpenERP 6.1,CamptoCamp GeoEngine modules(sales,partner) ,Python 2.6,Shapley 1.2.18,Geojson 1.0.1,   POSTGIS 1.5      When Starting server in Service it throws below Error: 12:43:00,283 1904 INFO ? openerp: OpenERP version 6.1-20130918-233418 12:43:00,283 1904 INFO ? openerp: addons paths: E:\OpenERP 6.1-20130918-233418\Server\server\openerp\addons 12:43:00,283 1904 INFO ? openerp: database hostname: localhost 12:43:00,283 1904 INFO ? openerp: database port: 5432 12:43:00,283 1904 INFO ? openerp: database user: openpg 12:43:00,315 1904 INFO ? openerp.service.netrpc_server: starting NET-RPC service on 0.0.0.0:8070 12:43:00,315 1904 INFO ? openerp.netsvc: Starting 1 services 12:43:00,361 1904 INFO ? openerp.wsgi.core: HTTP service (werkzeug) running on 0.0.0.0:8069 12:43:00,440 1904 INFO ? openerp.addons.web: embedded mode 12:43:00,642 1904 CRITICAL ? **openerp.modules.module: Couldn't load module web** 12:43:00,642 1904 CRITICAL ? **openerp.modules.module: No module named shapely.wkb** 12:43:00,642 1904 ERROR ? openerp: Failed to load server-wide module `web`. The `web` module is provided by the addons found in the `openerp-web` project. Maybe you forgot to add those addons in your addons_path configuration. Traceback (most recent call last): File "openerp-server", line 252, in File "E:\OpenERP 6.1\Server\server.\openerp\modules\module.py", line 419, in load_openerp_module File "E:\OpenERP 6.1\Server\server\openerp\addons\web__init__.py", line 29, in wsgi_postload File "E:\OpenERP 6.1\Server\server\openerp\addons\web\common\http.py", line 453, in **init** File "E:\OpenERP 6.1\Server\server\openerp\addons\web\common\http.py", line 523, in _load_addons File "E:\OpenERP 6.1\Server\server.\openerp\modules\module.py", line 139, in load_module File "E:\OpenERP 6.1\Server\server\openerp\addons\base_geoengine__init__.py", line 8, in File "E:\OpenERP 6.1Server\server\openerp\addons\base_geoengine\geo_model.py", line 11, in File "E:\OpenERP 6.1\Server\server\openerp\addons\base_geoengine\geo_field.py", line 9, in **ImportError: No module named shapely.wkb** 2013-09-24 12:43:00,642 1904 INFO ? openerp: OpenERP server is running, waiting for connections...


title: How to export Polygons to CSV with coordinates?

text: I am very new to QGIS, so please bear with me. I have a ".shp" file, which is contains a set of boundaries for suburbs of Australia. I opened it to QGIS and I would like to export it into a CSV file so that the final file would be like following. POSTCODE, LATITUDE, LONGITUDE When I select "Layer > Open attribute table" I can see the "postcode" column. I had a look at a similar question here Get list of coordinates for points in a layer? Which told me to simply go to "Layer > Save as" and then type "GEOMETRY=AS_XY" or "GEOMETRY=AS_XYZ" in the "Layer" section of the "OGR creation options". I followed exactly, but the latitude and longitude values are not being returned. I feel like I'm quite close! Could someone please help me. Thank you.


title: ArcGIS 9.3.1 cannot delete a version

text: We are having trouble deleting a version made in our 9.3.1 sde. When we attempt to delete it, ArcCatalog returns a message saying Unable to Delete - Target State not found in States table. Upon looking at the states table there is no reference to that particular version. Is there any other way to delete this version? We will be upgrading to Version 10.1 soon, will this cause a problem or can it be deleted during the migration process? Thanks!


title: How to remove decimal places from currency formatting in ArcGIS?

text: For years I have put up with the 2 decimal places that ArcGIS adds to the "currency" formatting, but most of the time I do not want them and end up removing them in my exported map images. I cannot see a way to edit the number of decimals - does anyone have any ideas how to get rid of these decimals?


title: Openlayers, valid zoomlevels are offset by 1 level. Why?

text: I have an openlayers based web app, map object is initalized like this:               p4326 = new OpenLayers.Projection('EPSG:4326');     p3857 = new OpenLayers.Projection('EPSG:3857');          //initialize map object, bind it to map_canvas div     map = new OpenLayers.Map('map_canvas',          {         //allOverlays: true,         projection: p3857,         displayProjection: p4326,         maxExtent: new OpenLayers.Bounds(-20037508.3427892,-20037508.3427892,20037508.3427892,20037508.3427892),          numZoomLevels:20,          maxResolution:156543.0339,         units:'m',          projection: p3857});      I have a prerendered map from levels 8-16, however, in order to view zoom level 16 of the map, i have to initialize its layer like this:               var osm = new OpenLayers.Layer.XYZ("OSM Mapnik", "http://localhost/osm_tiles/osm/${z}/${x}/${y}.png", { isBaseLayer: true, sphericalMercator: true, numZoomLevels: 17});     map.addLayer(osm);      Notice, that i set the **zoom level 17** (or i cant zoom and view level 16 of the map). Can someone tell me why?


title: Using IGeoProcessor2 to update Oracle Spatial

text: In the below example the input and out are shape files read from folders.Is it possible to read and write Oracle Spatial Geometry tables? I have ArcGIS,ArcSDE and Oracle Spatial 11g.   How to set parameters to create a new Geometry table if it does not exist?   How to append the data if a geometry table already exists ?               IGeoProcessor2 pGeoProcessor = new GeoProcessorClass();     IVariantArray parameters = new VarArrayClass();     parameters.Add(pFlayer);     parameters.Add(OutLayerPath + "\\" + pFlayer.Name + "MBM.shp");     parameters.Add("CIRCLE");     parameters.Add("NONE");     parameters.Add("");     parameters.Add("MBG_FIELDS");     pGeoProcessor.Execute("MinimumBoundingGeometry_management", parameters, null);


title: gxp_wmsgetfeatureinfo GetFeatureInfo on a non-advertised layer possible?

text: So, as the title states, I will like to know if it is possible to use gxp_wmscsource to send an GetFeatureInfo query to a non-advertised layer that was added to my opengeo SDK app using lazy loading. I made my app and added a gxp_wmsgetfeatureinfo to it, but if a click over a non-advertised layer, the app doesn’t even sends a GET, I suppose it is because there is not a list of available layers... I know it is possible to send a WMSGetFeatureInfo to a layer even if it is not advertised but I don’t know how to tell gxp_wmsgetfeatureinfo to do it.


title: Is it possible to incrementally update added/edited datasets to ArcGIS Online?

text: My organization would like to run change detection and update our ArcGIS Online data incrementally. Rather than re-uploading the entire data set every time we want to update the map. We are a city government that has constant edits being made to all of our data. We'd like to use FME to detect changes to the data and push these out to ArcGIS Online on a schedule. We have been using Google Maps Engine to do this and are looking in to purchasing a license. I'd hate to spend the money if there is a way to do this in AGOL that we are somehow missing.


title: Apend child to a layertree with nodeType: "gx_layer"

text: Im building an application with geoext and opengeo suite SDK. I´m trying to create a layer tree with some "grouped" layers like the last node here I saw somewhere that I should set the node Type to "gx_layer" to get that effect, but if I do that then the entire tree won’t load. I only get the root node and when I click on it , the icon changes to the “loading” image, but then nothing happens


title: OSM exports at different scale to what is on screen

text: I am using ArcGIS 9.3.1 I have added OSM as a WMS layer. I have a map set up how I want it to look in the layout view. When I go to export the map, whether as an image file, or in AI, which is the format I want to use, the map comes out looking different to how I have it on screen, basically the text on the OSM layer is much smaller. Does anyone have any suggestions for how I can stop this happening?


title: Project point to map with Arcpy using GPS Coordinates

text: I am trying to create a new layer with just a single point using arcpy. All that is provided is a gps lat and long. It seems like it should be a simple problem but I just can't get it to work. I've determined the error is that I'm not using the correct spatial reference (the point I get right now is next to the equator/meridian) but I'm not sure how to fix that.               import arcpy          arcpy.env.workspace = 'c:\\Scrap'     arcpy.env.overwriteOutput = True          latitude = 41     longitude = -71          p = arcpy.Point(longitude, latitude)     pg = arcpy.PointGeometry(p)     arcpy.CopyFeatures_management(pg, 'test.shp')


title: How can I do more flexible matching in conditional expressions than provided by LIKE?

text: With LIKE '%' works but I have had no luck with character ranges, i.e. [0-9] to match a digit. Oh for a regular expression ;)


title: creating shortest route using line segment layer as network locations

text: Using Arc10.2.1, trying to create a truck route that begins at one location on a county street grid and ends up in the same place, obeying a feature class of guard rail line segments throughout an entire county - we need the shortest possible. I need to use these line segments as network locations. However, the only line-based NA objects that can be defined as network locations are line and polygon barriers. But I'm not trying to avoid each guardrail, just follow it. I tried creating centroids and vertices from each segment to use as stops (on a single road as a test), but found out it sends me on this 1000 mile rat race to connect 3 points in a neighborhood road. Is there a problem with how the tool is interacting with the LRS itself? This seems like it would be a common problem, yet I haven't been able to find anything on it. Any help would be appreciated. Thanks!


title: Batch geocoding 1.5 million data points in SQL Server

text: I am currently attempting to batch geocode 1.5 million points that are stored in a sql database. The points are not very specific (ie only city, state points). I have used ArcGIS in the past but due to the new ArcGIS online credit system, it is restrictively expensive to geocode with their database. For that reason, I was able to access a composite geocoder that I attempted to connect to ArcGIS via the Create a geocoder feature, however, this system is not working as well (either it freezes with my database, or if I use a test dataset, it cannot match any addresses). I am seeking suggestions on:   1. How to best setup the geocoder I have to run in ArcGIS or    2. Other methods that may work better (I am really open to anything!) Let me know if there is more helpful information that I can provide.


title: How to calculate polygon areas and perimeters?

text: If I load a polygon layer and make a new column, how do I fill that column with the 1) area of each polygon and another column with 2) the length of the perimeter of each polygon?


title: How to calculate area in km2 Qgis 2.0

text: I have a simple question. I have calculated the area of the polygons in my shapefile, but the number seems confusing for me. I need the area to be shown in km2 but instead I get numbers like:0.004008 What are the units of this number? When I used the "measure area" icon, that is, when I measured the same polygon manually the size of the polygon shows me: 50.342km2 See, the two numbers are completely different. How is that?? Thanks for any help!!


title: Stamp out and transfer values from one overlayed vector layer to another

text: I have a vector layer with administrative boundaries and another vector layer with postal code boundaries, which are more precise/smaller entities/objects. I now want to transmit the values of one column/attribute of the first layer to the second. Is it possible to overlay them and kind of stamp out this attribute and thus transfer it to the postal code layer? Thanks for any information Marc


title: Examples of 3-D Maps

text: Similar to beautiful maps example is there any similar collection for 3-D maps? Thanks in advance : )


title: Is there a way to perform zonal statistics-like operations with feature zone data without doing feature-to-raster transformation?

text: I want to use zonal statistics on Chinese counties (which can be very small) from a fairly coarse (0.5 degree) raster grid. The problem with zonal statistics is that it returns a much smaller number of counties after the feature-to-raster transformation. This means a lot of my counties won't be assigned any data. Can this be overcome somehow? For example, I would like for example to have a list of raster cells that lie within each county, along with the cell's value and the share of the cell that lies within the county. This way, I would be able to construct a nice measure.


title: How to migrate PostGIS spatial data to MsSql spatial data?

text: am trying to migrate PostGIS geometry to SQL server Spatial data. which ETL i have to use. it's the first time i'll work with SQL server Spatial. is there any Data ETL model that do this migration. thanks for help


title: Is it possible to use USGS .shd files in QGIS?

text: Has anyone successfully used QGIS with US Geological Survey shd files to map geology - specifically to display colors and geological symbols (dip/strike etc) stored in a USGS supplied shd file.   Looking for a non arc solution.


title: Google maps indoors submission

text: Does anyone here have experience with submitting Google indoor maps? I submitted one two weeks ago where the submission page stated "submitted for processing". (I submitted the indoor map on behalf of those I am delveloping an app for.) I have not received any confirmation email which has worried me a little since after searching Google, I can see people in 2012 appeared to get such emails. Is there anything I need to do beyond uploading the map? (matching the building etc. of course) - i.e. are there hidden requirements that an email or mail is sent to the address mapped or something like that? If so, I need to contact the one I am developing an app for :)


title: Geodjango - User's uploaded shapefiles and model creation

text: I have a big interrogation about how it could be possible to have a GeoDjango app that can store in database the uploaded shapefiles (or other geospatial files) from users. When you load shapefiles in GeoDjango, you have to create a new model that matches the fields of the shapefile. How could it be possible to create a new model automatically each time a user uploads a files? You can't append the new model definition to your models.py each time? For me, that does not make any sense. Does it means the best way is to use a generic model definition and only allow files upload that match this definition,and then to create a new instance of this generic model? Can someone explain me what is the best practice to create an app that can store users uploaded spatial files?


title: Convert multiple lines into single line

text: I've no idea if this is basic or what but I have this data that contains multiple lines (forming a polygon shape). I want to convert those lines (multiple - each with a start and end point) to a single line. I tried dissolve, multipart to single part etc. I cannot use topology since it influences precision. I have zero knowledge of programming. any solutions will be appreciatable. thanks.


title: pgsql2shp error when query is introduced

text: I'm trying to export a table from PostGIS using pgsql2shp. The basic query to export the entire table works fine using the following format:               "C:\Program Files\PostgreSQL\9.1\bin\pgsql2shp.exe" -f "D:\path\shapefile.shp" -h     localhost -p 5432 -u postgres -P xxxx mapdata planning.dm_planning-apps      However when I add a query (which works fine in PGAdmin III) the command fails, due to the hyphen in the table name. Working query in PGAdmin:               SELECT * FROM planning."dm_planning-apps" WHERE "dm_planning-apps"."Modified_Reference"     ILIKE '%np%'      Included in command:               "C:\Program Files\PostgreSQL\9.1\bin\pgsql2shp.exe" -f "D:\path\shapefile.shp"  -h     localhost -p 5432 -u postgres -P xxxx mapdata "SELECT * FROM planning."dm_planning-apps"     WHERE "dm_planning-apps"."Modified_Reference" ILIKE '%np%'"      Produces a "syntax error at or near "-" and highlights the hyphen in "...FROM planning.dm_planning-apps WHERE...". The double quotes obviously aren't doing the same job they were within PDAdmin where the query worked fine. Can anyone suggest a way to make this query work without renaming tables?


title: SSHFS won't read / write ArcMap documents (.mxd)

text: I am trying to use SSHFS to access our file server and work on geographic data. I can mount the drive, get arc catalog to connect to the folder, and I can access and edit geodatabases and shapefiles. However, when I try to open a Map document, it says "Failed to open document", click ok, and "The specified file does not contain a valid ArcMap Document". I also cannot open a blank map, add all my data and then save it to the mounted drive, it gives a similar error. I know the documents work, I can download them and they work just fine. Is this a problem with ArcGIS or with my sshfs connection?


title: Leaflet Polyline with custom fill pattern

text: I'm drawing Lines with leaflet.js an have to create custom styled polylines. So i got various image files which should be like a repeated background image in HTML. So is it possible to give a drawn L.polyline a Custom background image (texture)? I know SVG has Support for that but is It also possible via leaflet API?               <defs>         <pattern id="bad" patternUnits="userSpaceOnUse" x="0" y="0"          width="10" height="10" viewBox="0 0 10 10" >             <image x="0" y="0" width="10px" height="10px"              xlink:href="img/non-sprite/grade_0.png" />         </pattern>          <defs>      And tried to use the pattern by assign its ID as color/fill without any luck.               L.circle(lat, lng, 200, {fill : 'url(#best)'}).addTo(map);      The circle is shown on map but only in color 0033ff (Leaflet default line color)


title: Can ArcMap do HTML pop up with a chart?

text: am using Arcmap 10.0 to create a map. Need to display a chart for my data when users click on the marker. Must I use plugin like google chart? Can I do it directly in Arcmap?


title: PgRouting-How to calculate only walking routes

text: I am using osm2pgrouting to import data to my PostGis database from OSM and I was wondering if there is some way to specify to PgRouting to calculate paths with only ways that can be walked by people. I've been able to make paths with any kind of ways and the paths use the roads, something that I don't want to do. Thanks in advanced!


title: How to setup PostGIS for use with multiple machines?

text: I'm a beginner user of postgresql/postgis. I am starting to use Linux more and more but still rely on windows (on different machines). How would I setup a postgis database that I could work with under different systems? Thanks, Eelke


title: How to convert bearing, distance, xy coordinate measurements to point shapefile?

text: I have a natural resources dataset where individual seedlings in canopy gaps were recorded from a center "legacy tree" within each canopy gap. Each sapling has the following measurements associated with it:   1. Legacy tree GPS xy coordinates (i.e. a center location in the canopy gap)   2. Bearing to sapling from the legacy tree (e.g. 218 degrees)   3. Distance in meters from the legacy tree There are thousands of these individual tree locations, so an automated approach would be best. How can I convert this dataset into a point shapefile/featureclass? ArcGIS, Arcpy/Python and R solutions would be especially useful. Thanks.


title: Unable to edit data in an ArcSDE geodatabase

text: I'm dealing for the first time with an ArcSDE geodatabase so I'm facing simple problems to you. ESRI help wasn't of much help (to a begginer at least). I can't edit my data in ArcMap, most probably because I've already ""registered as versioned" my featuredataset. I did use "Unregister as versioned" from ArcCatalog, but still can't edit my data. And then I got some "schema lock error"" :) So any clues? Ps: I'm the only user of the geodatabase. Thank you very much, Ligia


title: How to control visibility of cached layer in ArcGIS Server map service?

text: I am consuming mapservice form ArcGIS Server 10 using ArcGIS Silverlight API. I have created cache (fusion layers) for my services by right clicking on the mapservice in ArcCatalog and clicking create cache. It works fine in my application but the problem is after viewing the cache when do map full extent the cache is still visible when it had to hide since I have set zoom visibility. Please guide me as to how I can keep the map cache visibility off when it is out of the zoom visibility.


title: How to view/edit attributed relationship data in Identify/Attributes dialog in ArcMap?

text: Is there any extension or standard option to view/edit attributed relationship records attributes using standard ArcMap attribute editor dialog? Same to question about identify dialog :) Thanks


title: How to protect my application ESRI FGDB data from authorized access using other applications?

text: I have ArcObjects app which uses file geodatabase for storing data and performing CRUD operations. This file geodatabase needs to be protected from reading data in other applications such as ArcMap, Civil 3d and other 3rd party's. Is there any way to do this? Thanks


title: How to turn off Google map data copyrighted window in Google Street layer (Qgis - Openlayers Plugin)?

text: > **Possible Duplicate:**   >  How to close watermarks in openlayers plugin After a few weeks not using Qgis, this morning I opened it up, accepted a request to update Openlayers Plugin. Right now, I found a window saying "Map data ©2012 Google, Mapabc" at the center of the screen, following me every where I drag the map. Any one experienced this, please tell me how to remove that screen? I am using Qgis 1.7.4 on Windows XP, Openlayers Plugin version 0.90. ![enter image description here](http://i.stack.imgur.com/38gLN.png)


title: Trouble with Python plugins on QGIS

text: This question has been bugging me for a while now, and I still haven't managed to solve it, so I'm giving it another go, hopefully someone may help me out. I've been having problems loading several QGIS plugins. I've just done a reinstall of QGIS 1.7.4., followed all the instructions, so just did a default installation (I've tried the osgeo4w installer already and it only made matters worse...). Whenever I try to use, let's say the "1 band raster color table" plugin, I get error messages. First of all, when I load the plugin, I get the message "Cannot continue. Problems!". If I close this message, the plugin stays open anyway and let's me select the options I want. When I click "Create", I get yet another error message, that starts off with this:               "Traceback (most recent call last):       File "C:/Users/Tiago/.qgis/python/plugins\bcccoltbl1\bcctbl.py", line 826, in applyTBL         myStats = self.myGRD.getStats()     AttributeError: 'str' object has no attribute 'getStats'"      Now, I've been told that this might have to do with me having non-Latin characters in the path... I've been trying to understand this, but I'm getting really confused... When I used QGIS 1.6, I had no problems of this sort. But I didn't change a thing in my disk's organization when I upgraded to QGIS 1.7, so how can this be, that the plugins ran smoothly before but not anymore? If I had non-latin characters, they would be there before, anyway, so the plugins should have had these same errors... This is getting me really frustrated, because it will make it very hard for me to keep using QGIS if I can't even change the looks of rasters! Please, I really need help with this. If necessary, I can send some more information. Thanx!!!


title: GeoServer automate "reload feature type"

text: Using GeoServer, if you add a feature type bound to a SQL table that is empty, when you then populate that table no data shows up in layer preview. Clicking the "reload feature type" button corrects the problem. Is there any way to do this through the REST API, or some other way to automate it? It looks like the resource reset command could achieve this, but it seems to do so for all stores. I would rather only do it for the store in question.


title: Why does QueryFeatureIDs return an empty array when querying a Route Event Source?

text: I'm struggling with using the ESRI SOAP api to perform a spatial filter against a layer with a datatype of 'Route Event Source' Because the number of features returned by a single spatial query could easily exceed the configured ArcGIS server max, I first grab all the features Ids of the features within the selection envelope, and then retrieve the feature attributes in chunks. No matter what I do, the array of Feature Ids returned by `QueryFeatureIds` is _always_ empty. With the same `SpatialFilter`, though, I receive a correct count from `QueryFeatureCount` and I receive the correct records from `QueryFeatureData`. Here's the test case that illustrates the problem I'm having. Am I missing something, or is that the intended behavior? (I've changed the selection envelope in the sample)               var proxy = new MapServerProxy("http://mymapserver/ArcGIS/services/myservice/MapServer");          var envelope = new EnvelopeN();     envelope.XMax = 100010.0;     envelope.XMin = 100000.0;     envelope.YMax = 200010.0;     envelope.YMin = 200000.0;          var spatialFilter = new SpatialFilter();     spatialFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects;     spatialFilter.GeometryFieldName = "Shape";     spatialFilter.FilterGeometry = envelope;          int routeEventLayerIndex = 0;     var count = proxy.QueryFeatureCount("Layers", routeEventLayerIndex, spatialFilter);     Assert.AreEqual(2, count); // PASSES, 2 features within envelope          var records = proxy.QueryFeatureData("Layers", routeEventLayerIndex, spatialFilter);     Assert.AreEqual(2, records.Records.Length); // PASSES, data for 2 features is returned          var fidSet = proxy.QueryFeatureIDs("Layers", routeEventLayerIndex, spatialFilter);     Assert.AreEqual(2, fidSet.FIDArray.Length); // FAILS!!! Length of FIDArray is 0      **Edit**   I've tested this against ArcGIS Server 10.0 and 10.0 sp 4; both return the same results.


title: Find centroids of polyline and create point using Arcobject

text: I have one road feature class and one area feature class .What i have to do is i have to create a point at the centroids of the polylines of (road feature class) and find intersection of point with the area boundary. How can i get the centroids of polylines in arcobject through vb.net code and create point with that (x,y) and then find intersection with area boundary????


title: How can save and restore my map's extent with GeoExt?

text: I'm using GeoExt with OpenLayers and trying to save the extent of the map. while exploring the page, users can be able to save the extent and restore via their browsers. Here is the point that I'm stuck: **how can I catch mousewheel's actions to save the extent of the map and keep it in cookies?** any ideas?


title: Openlayers TMS default image

text: In OpenLayers, when a tile is missing in a TMS layer, a semi-transparent pink rectangle shows in its place. Is there a way to load some default image, like one saying "No data for this region under this zoom" instead? One possibility is to request the tiles through something like a PHP script or a Java servlet. If a nonexistent tile is requested, this servlet would then send the default tile. However, I'm currently taking the tiles directly from the filesystem, so I'm looking for a solution to be done only with JS / OpenLayers.


title: Zoom to Polygon in Arcmap from Access

text: I am attempting to Zoom from an Access record to a specific polygon in ArcMAP but I can't seem to get the function to work. The field in Access is identical to the field in arcmap. I.E. they both have a unique polygonID (custom field) identifier. I have been able to open the specific record in Access from ArcMAP but would like the functionality to work both directions. Does anyone know if there is any examples of this online and if it is even possible/worth doing with VBA being phased out. Thanks


title: Is the loop in this python code breaking because of duplicate values in the column?

text: This code is meant take a shapefile in the current mxd, loop through all of the values in a collumn of said shapefile. For each value it will create a definition query, zoom to features then save a copy of the mxd. When I run the code it completes, but only generates one mxd when I'm expecting about 25. Is the loop breaking due to duplicate values in the Field_1 collumn? If so how can I work around, if not what is breaking the loop after one single iteration?               import arcpy, os     from arcpy import env     #Overwrite features     arcpy.env.overwriteOutput = True     shape = r"C:\NewProjects\shapes\Meters.shp"     output = r"C:\NewProjects\shapes"     layer = r"C:\NewProjects\shapes\newlyr"     mxd = arcpy.mapping.MapDocument(r"C:\NewProjects\shapes\Try.mxd")     df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]     inset = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]     lyr1 = arcpy.mapping.ListLayers(mxd, "Meters", df)[0]  # pick the layer you want to "extract" data from     #For each work order:     rows = arcpy.SearchCursor(lyr1, "", "", "FIELD_1", "FIELD_1 A")     for row in rows:         name = str(row.getValue('FIELD_1'))         sqlExp = """"FIELD_1" ='"""  + str(row.getValue('FIELD_1'))+"""'"""    # Create a clause to select only the current record         lyr1.definitionQuery = sqlExp         df.extent = lyr1.getSelectedExtent(True)     newMxd = "Figure1_" + name     mxd.saveACopy(r"C:\NewProjects\shapes\\" + newMxd + ".mxd")     print "Saving ", newMxd + ".mxd"     print "All done."     del mxd


title: Filtering a Topojson on server-side. Which part of topology I need?

text: I have a big Topojson which I want to avoid loading it at page load. Instead I would like to parse on server-side at user request, sending back chunks of this original Topojson to be ploted. So I want to filter by a list of attributes in topojson. Initially I thought that just sending its corresponding geometries with its arcs would be OK (sorry, i am used to GeoJSON) but then I see below the Topology section, which I am afraid is quite dependent on the previously selected arcs (in "Geometry"). Topology section is quite a very long part of this TopoJSON. The question is, as I cannot find a link between the arcs in "Geometry" and "Topology", do I have to load all Topology part, even if I only select a few of "Geometries"? it would not solve my problem at all… any alternative? I want to filter by COD_POSTAL                 "geometries": [         {           "type": "Polygon",            "properties": {             "COD_POSTAL": "08630"           },            "arcs": [             [               0,                1,                2        ]     ]     },         {           "type": "Polygon",            "properties": {             "COD_POSTAL": "08256"           },            "arcs": [             [               8,                9,                10,                11,                12,                13             ]           ]         }


title: Additional Basemap Menu in ArcGIS Viewer for Flex?

text: I am using ArcGIS Viewer for Flex 3.5 to create an application. As it covers more than one areas of interest, I need to provide an option to the user to only select base maps from his own desired domain. In order to do that, I need more than one Basemap dropdown menus in the application whereas, the default map switcher widget only provides options for one menu. I have not been able to find any worthwhile support or guidance for the same. The only relevant source of information which I could find is a thread related to similar task in sample flex viewer in old ESRI forums at this link: http://forums.esri.com/Thread.asp?c=158&f=2421&t=282643. However, the flex code when compared the 3.5 version is no more relevant and I cannot figure out how to make use of that.


title: GeoGratis API, joining categories in a single request

text: I am a developer trying to work with the GeoGratis web API, provided here: http://geogratis.gc.ca/api/en I am having difficulty with the use of catagories in the API, and am wondering who I could talk to about this. Specifically, I want to make a request like this: http://geogratis.gc.ca/api/en/nrcan-rncan/ess- sst/-/OR/(urn:iso:type)map/(urn:gc:subject)climate but no response is given. This request does work: http://geogratis.gc.ca/api/en/nrcan-rncan/ess- sst/-/AND/(urn:iso:type)map/(urn:gc:subject)climate The difference between the two requests is the use of 'OR' in the first and 'AND' in the second. I think I should be allowed to use 'OR' as well as 'AND', based on the documentation here: http://geogratis.gc.ca/api/en (you have to open the 'documentation' subheading, under 'Developers', and then scroll all the way down to the bottom of the page.) Is there anyone that can help me with this?


title: Atmospheric correction for RapideEye without atmospheric parameters

text: I'd like to perform an atmospheric correction of my RapidEye images but metadata do no include atmospheric parameters required to use the 6s algo implemented in the Orfeo ToolBox. Is there a solution to still apply the corrections? Browsing on the web, I found that in grass it was possible not to specify them (http://grass.osgeo.org/grass64/manuals/i.atcorr.html) but rather define if the site is midlatitude summer, midlatitude winter, subarctic summer, subarctic winter, etc. Do those areas have predefined parameters and if yes where can I find them?


title: OpenStreetMap tiles not loaded despite being received

text: While using OpenLayers to load OpenStreetMap, I observe that part of the map would sometimes be covered with blank pink tiles. I searched and found similar complaints with solutions suggesting the modification of `OpenLayers.IMAGE_RELOAD_ATTEMPTS` parameter to an integer number greater than 0. I did exactly that, but it doesn't solve the problem. Using Firebug, I am able to verify that the tiles were indeed downloaded successfully. However, for some unknown reason, they just do not load onto the map. I have no issue with Google Maps as a base layer. The problem seems to be restricted to OpenStreetMap only. Any idea how to solve this pink tile of death? ![OpenStreetMap with pink tiles](http://i.stack.imgur.com/WGCeX.png)


title: Advantages of R-trees in comparison to geohashes

text: Geohashes are being widely used in products like: Lucene, mongodb, etc and have become one of the most important technology of present day. Have Geohashes replaced the good old R-trees or do R-trees have any advantages as compared to Geohashes. I'll be grateful if someone can help me understand the power of R-trees over Geohashes.