title: Is ArcMap 10.0 Query Layer passing invalid SQL to Oracle 11gR2?

text: In two different environments I'm having trouble using ArcMap 10.0 Query Layers to view SDO_GEOMETRY (point) data in Oracle 11.2. This suggests I'm doing something wrong... Points are in Oracle stored under SDO_GEOMETRY, 2001 type (2D point), 4326 SRID (WGS84 lat/lon). The table / column has a spatial index and is entered in user_sdo_geom_metadata. I can successfully query the table in a way that uses its spatial index (SDO_WITHIN_DISTANCE), so I know this is functional. When I add the table to ArcMap as a query layer the query validates fine, however nothing is displayed when I zoom to the layer. I have verified that the zoom extent is correct. If I try to show labels I see a drawing error:               Underlying DBMS error[ORA-29902: error in executing ODCIIndexStart() routine     ORA-22060: argument [2] is an invalid or uninitialized number     ORA-06512: at "MDSYS.SDO_INDEX_METHOD_10I", line 333     ]      Not understanding how this involves me I decided to trace the database. ArcMap makes the following query to Oracle:               select LOCATION_GEOMETRY     from (select * from TEST.BC_POINT_TEST_1) a     where sdo_filter(LOCATION_GEOMETRY,sdo_geometry( to_blob( :i1 ), :i2),'querytype=window') = 'TRUE'      I filled in the blanks to create the following query:               select LOCATION_GEOMETRY     from (select * from TEST.BC_POINT_TEST_1) a     where sdo_filter(LOCATION_GEOMETRY, SDO_GEOMETRY(2003, 4326, NULL, SDO_ELEM_INFO_ARRAY(1, 1003, 3), SDO_ORDINATE_ARRAY(53.276, 23.232, 53.587, 23.451)),'querytype=window') = 'TRUE'      Which worked fine. However, if I replace `4326` with `NULL` it throws a similar error (but not exactly the same) to what I'm seeing in ArcMap. This makes me think possibly ArcMap isn't passing the correct SRID for Oracle to know what to do. I set the data frame's coordinate system to WGS84 but this didn't make a difference. The error reported through ArcMap references `argument [2]`. If this is actually the second bind parameter in the traced query that would be the `SDO_ORDINATE_ARRAY`, but as this is an array it would make no sense to complain that this is `an invalid or uninitialized number`. If anyone has prior experience with this issue or suggestions on where to take it I would be most appreciative! **Edit** : so far nothing has helped solve this problem, but I recently had an opportunity to test with 10.1 and everything worked as expected. This suggests there's nothing wrong with the Oracle instance.


title: Create multiple vector layers by reading a GeoJSON file once in OpenLayers

text: I am a begginner to OpenLayers, and my question is about making things more efficient. My work is about displaying time-dependent weather elements on several point locations on the map. All these are initially stored in a MySQL database, and I use a PHP file to convert a (parameterized) query result into GeoJSON format. All the features are points with a pair of coordinates, but many properties (like name, temperature, dewpoint, wind direction, wind strength etc.) With a refreshing strategy, I manage to keep the elements up-to-date. Each property is supposed to be displayed in a different layer on the map, so that the user can show or hide whatever finds or not necessary. There are a names' layer, a temperatures' layer, a dewpoints' layer and so on. The way to create two of these layers can be presented at the last lines of the following code:               var wgs84 = new OpenLayers.Projection("EPSG:4326");     var mercator = new OpenLayers.Projection("EPSG:900913");          var myURL="http://localhost:8080/wxproj/wx-json.php";          var clustering = new OpenLayers.Strategy.Cluster({distance:40});     var refreshing = new OpenLayers.Strategy.Refresh({interval:1000});          var clusters = new OpenLayers.Layer.Vector("LIGHTNINGS", {          strategies: [clustering, refreshing, new OpenLayers.Strategy.Fixed()],         styleMap: myLightningStyles,         protocol: new OpenLayers.Protocol.HTTP({             url: myURL,             format: new OpenLayers.Format.GeoJSON({                 'internalProjection': mercator,                 'externalProjection' : wgs84 }),             async: false             })         });          var stations = new OpenLayers.Layer.Vector("STATIONS", {          strategies: [new OpenLayers.Strategy.Fixed()],         styleMap: myStationStyles,         protocol: new OpenLayers.Protocol.HTTP({             url: myURL,             format: new OpenLayers.Format.GeoJSON({                 'internalProjection': mercator,                 'externalProjection' : wgs84 }),             async: false             })         });           map.addLayers([clusters, stations]);      The obvious drawback is that I have to call the same file "wx-json.php" several times, one for each layer, getting in almost all cases the same GeoJSON-formatted result exactly. While it is sufficient to call it once, every time needed, and then distribute the properties (and its values) to the layers (and their features, accordingly). One thought about bypassing this result, is appending parameters to the URL, so that I get the GeoJSON-formatted features with one property at each call:               var myURL="http://localhost:8080/wxproj/wx-json.php?property=station";     var stations = new OpenLayers.Layer.Vector("STATIONS", { ... });          var myURL="http://localhost:8080/wxproj/wx-json.php?property=lightnings";     var clusters = new OpenLayers.Layer.Vector("LIGHTNINGS", { ... });      But I find it also consuming. Is there any more efficient way to get the GeoJSON data only once and then create the layers with features according to each property? Edit: I found a similar question that restates my own demand, that is a single GeoJSON file backing multiple Vector Layers. But, there is no code to understand how this is being done. * * * Well, it has been a month without some answer, so I suppose my question needs some refrase: What I really need with OL is filling a number of overlay layers with features by **reading once** (in a period) a GeoJSON file. I get each half-an-hour a different GeoJSON file, which corresponds to a lot of point features, each with a pair of coordinates and some properties. For each one of the properties, I create one overlay layer in OpenLayers. Essentially, there are the same points in all layers, but a different property is hightlighted for every layer. The way I went so far, was reading the same GeoJSON file as many times as many are the layers. This way has been expressed in the first block of code above, at my first question one month ago: Each layer manages to get its own features reading for itself from the PHP file, which produces a GeoJSON format output. However, this way is not desirable, because of the too many independent calls to the webserver, while the layers display essentially the same thing; what's more, if the PHP file changes its output in the meantime of layers' reading, the displayed result is prone to confusion. My question is therefore about another way of loading my features. I want to read once, in a period say half-an-hour, the GeoJSON (output from the PHP) and get all my layers filled from this.


title: How do I join an MS Access (accdb) table to a Feature Layer in ArcGIS 10.1?

text: I would like to join a table from a microsoft access database (.accdb) to a featurelayer in ArcGIS 10.1. How do I associate tables from that database with existing layers that have a matching attribute?


title: Select a single entry using SearchCursor and create variable

text: I am trying to select a single entry of a floating point entry in a dbf table that contains 9 entries (sampled from a raster image), some of which is NULL using SearchCursor. I want to select one of the entries that is not NULL and assign it as a new variable for use later in my script. I can do a query that excludes NULL values (includes all values greater than 0) but I don't know how to select a single entry that is not NULL and assign as a new variable. The field name is:               Rt_b1     <Null>     0.138817     0.136027     0.136027     <Null>     <Null>     0.136027     0.138817     0.136027      Any help would be much appreciated. I'm using Windows XP and ArcGIS 10.1               with arcpy.da.SearchCursor(Sample_TIF2, ("Rt_b1"), '"Rt_b1" > 0') as cursor:          for row in cursor:                print str(row[0])


title: Cannot get changeAttributeValues() to work in a loop

text: I'm trying to use QGIS to loop through a set of points, find their elevation from a DEM and use that value to update the point layer's attribute table. The loop works fine until I try to use changeAttributeValues(), and I can change the attribute values outside of the loop, but when I try to put the two together I cause a dump.               while provider.nextFeature(feat):         geom = feat.geometry()         x = geom.asPoint()         res, ident = rlayer.identify(QgsPoint(x))              for (k,v) in ident.iteritems():             elevation = float(v)         attrs = { 2 : elevation }         caps = vlayer.dataProvider().capabilities()         if caps & QgsVectorDataProvider.ChangeAttributeValues:             vlayer.dataProvider().changeAttributeValues({ fid : attrs })      This causes the following error in Linux:               *** glibc detected *** python2: free(): invalid pointer: 0x00000000024655b8 ***      Edit: The problem isn't the loop, it is now I'm calling changeAttributeValues(). I've simplified the script down to just update a single point and it causes the same core dump. Here's the updated script:               #!/usr/bin/env python2     # -*- coding: utf-8 -*-          from qgis.core import *     import qgis.utils          #initializes QGIS and points python to where it can find QGIS's stuff     QgsApplication.setPrefixPath("/usr/", True)     QgsApplication.initQgis()          vlayer = QgsVectorLayer("/gis/vector/elev_pts.shp", "elev_pts", "ogr")     if not vlayer.isValid():         print "Layer did not load!"          vlayer.dataProvider().changeAttributeValues({ 0: 2: 432.1 })          QgsApplication.exitQgis()


title: Geopaparazzi Export

text: I've been using geopaparazzi to collect some data in the field. During the collection I use the notes option to store some valuable information. But after the export of the data in gpx, there is no tags showing my notes.               <?xml version="1.0" encoding="UTF-8"?>     <gpx       version="1.0"       creator="Geopaparazzi - http://www.geopaparazzi.eu"       xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"       xmlns="http://www.topografix.com/GPX/1/0"       xsi:schemaLocation="http://www.topografix.com/GPX/1/0 http://www.topografix.com/GPX/1/0/gpx.xsd">     <time>2014-04-03T12:18:48Z</time>     <bounds minlat="0.0" minlon="0.0" maxlat="48.786244" maxlon="11.383844"/>     <wpt lat="48.786228" lon="11.38384">       <ele>419.7365417480469</ele>       <name>building</name>       <cmt></cmt>       <desc></desc>     </wpt>     <wpt lat="48.786244" lon="11.383844">       <ele>419.5691833496094</ele>       <name>Sol</name>       <cmt></cmt>       <desc></desc>     </wpt>     </gpx>      The different waypoints were supposed to get more information. i.e wpt Building was supposed to have also Heigt, name, number of floors and so on. Exporting the data in spatiallite seems not to be supported. Please, has someone any idea, how can export the data with all the relevant information? THX


title: Regarding MCD43 products

text: I have a special MCD43 daily product with 250 m resolution in .hdf format. Whenever I try to open it in the MODIS Re projection Tool, the error i get is "error in module:read hdr file". What should I do? I want to convert these files to .tiff format. Please help!


title: Should Delete Field work in ModelBuilder while there is a Join present that appears unrelated?

text: I am having an **issue** with a very simple Model, both in ArcGIS 9.3 and ArcGIS 10.0. The part of the Model that is causing trouble is shown below, along with ModelBuilder's error message (in ArcGIS 9.3): ![Screenshot of the ArcGIS model in question](http://i.stack.imgur.com/litvv.png) Let me summarise (from left to right):   1. I have a feature layer and a table view, which I joined using the _Add Join_ tool. This happens via a field pair named `ID` / `FKID` (which both have type `LONG`).    2. I am then exporting ("materialising") the result of the join as a new shapefile ( _Joined Feature Layer.shp_ ) using the _Copy Features_ tool.   3. Right afterwards, I am removing the join via _Remove Join_ , because I know from experience that joins in ModelBuilder easily cause all sorts of trouble.   4. At the same time, I am trying to delete a field from the exported shapefile via _Delete Field_ … but ModelBuilder won't let me connect the tool to the input shapefile, instead giving me above error message. This doesn't make sense. I've been able to track the problem down. This issue seems to be caused by the join between the originating table view and feature layer. But I cannot see how this join is still relevant to the exported shapefile, which by all means should no longer be connected to its originating dataset at all, or be influenced by it! My thinking is that connecting the _Delete Field_ tool to a feature class or shapefile should not be a problem. * * * A workmate has found the following **workaround** which involves an additional _Copy_ : ![Screenshot of a possible workaround](http://i.stack.imgur.com/dvPeN.png) Is there something wrong in my first Model, or is this a bug in ModelBuilder? Perhaps even a documented one?


title: GXP Hidden layer switcher - composer.html

text: I am working with GEOSERVER, GeoExplorer, GXP. I am creating a WPS client and need to generate a vector layer that serves to make operations on vector data, such as buffer, split polygon. I know how to create the vector layer but I want that is not on the menu layers GeoExplorer, I want this hidden. I try modify this script wpsdemo.js from (http://suite.opengeo.org/docs/latest/processing/wpsclient/index.html ) where: In wpsdemo.js replace:               ....     target.on('ready', function() {                 // Get a reference to the vector layer from app.js                 this.layer = target.getLayerRecordFromMap({                     name: 'sketch',                     source: 'ol'                 }).getLayer();                 // Some defaults     ....      for code that does the following:   1. create dynamically OpenlayerLayer.Vector (with `name: ''sketch', source: 'ol', displayInLayerSwitcher: false`) <---- I know how to do it   2. add this layer to the map <---- I dont know (my map is defined in composer.html)   3. GIS operations do but with the hidden layer in the layers menu switch I try the following settings (in the composer.html) using the property **displayInLayerSwitcher: false** but does not work :               {                     // A vector layer to display our geometries and processing results                     source: "ol",                     name: "Res",                     type: "OpenLayers.Layer.Vector",                     fixed: true,                     visibility: false,                     selected: false,                     displayInLayerSwitcher: false                        }      In what way I can hide the layer in the layer switcher of geoexplorer? Thanks in advanced


title: Clip tool not working: 000735 error?

text: I am trying to use the clip tool. When i select it from the list, a window opens however it is greyed-out. i cannot enter any imput or clip features. 'Environments...' does not allow me to select any items. ![enter image description here](http://i.stack.imgur.com/ZDw0I.png)


title: Multiple map layers : TPK

text: Working on ARCGIS for Offline mapping ,Can a tpk file can have multiple map layers inside it?


title: arcGIS observer point analysis: Exceeding maximum allowed points amount not exceeded

text: I need to make an observer points analysis (3D analyst tool). I have a shapefile of 22 points, so I divided it in 2 layers of 11 each as the maximum points allowed is 16. However, as I launch the process, the program says: "ERROR 010155: Too many points. Maximum number allowed is 16. Failed to execute (ObserverPoints)" ... But I have only 11! Could anyone help me please? Thank you very much!


title: How does Spatialite handle concurrent access?

text: I understand that Spatialite is meant to be used by a single user/machine at a time. However if the data is stored on server, there is no guarantee that someone on another machine will not also open the same database and use it. Has anyone got any experience with how QGIS and spatialite react to being in a multiuser environment. In most cases, the layer would be a read only layer. QGIS allows multiple people to read the same shapefiles, which is a similar situation. What happens if multiple people have concurrent editing sessions open?


title: How To Export Selected Data Driven Pages to PDF?

text: I have a map book with data driven pages that shows the fire boundaries, roads, etc for our fire rescue department. I want to automate the process of creating individual mini map books for each fire station that only include the pages within their boundary. I started off by creating a model that will iterate through the district boundary names and select all pages with those boundaries: ![enter image description here](http://i.stack.imgur.com/tNwtV.png) Now, I'd like to take the selected pages and export them to a pdf in individual folders. My goal is to have a folder for each district name and within the folder have individual pdfs of the pages that make up that district. I'm assuming the best way to do this portion is with `arcpy.mapping` so I have exported the model above to a python script:               # Import arcpy module     import arcpy          # Load required toolboxes     arcpy.ImportToolbox("Model Functions")               # Local variables:     Fire_Boundaries = "Fire Boundaries_exportPDF"     Grid = "HydrantGrid_exportPDF"     Selected_District = "I_Fire Boundaries_exportPDF_NAME"          # Process: Iterate Feature Selection     arcpy.IterateFeatureSelection_mb(Fire_Boundaries, "NAME #", "false")          # Process: Select Layer By Location     arcpy.SelectLayerByLocation_management(Grid, "INTERSECT", Selected_District, "", "NEW_SELECTION")      I'm not sure where to go from here with arcpy/python? Any help would be appreciated.


title: How to analyze and combine overlaping elevation maps into one?

text: I have 30 overlapping maps that contain tree density heights (measures between 1-100%) of a large area, and I need to get the highest point by pixel and create a new map with all these highest points if it makes sense. Basically, say in point (pixel) "A" in map 1, 2...30 the maximum elevation is 25 m, and in the next pixel "B" (in all the other maps) is 20. I would need to create a new map (a new matrix) that condenses all these elevations/densities (and will look like 25, 20, x, y, etc). Someone has suggested making a program to analyze all of the maps by pixel, but I'm new to this, so any help, or suggestion will be appreciated. Ok, based on a comment I got, I'll add more detail here about the data: I have the data in three different formats (the same data) these being tif, .dtm. and ASCII. I have a section of a river, and the 30 maps overlap on top of each other exactly. They have measurements of three density (0, where there is no "vegetation" in the pixel, to 100, all is vegetated), now each of these overlapping maps contain this same density (0-100) for different heights, namely 1-30m. For example: in map 1, I would have calculated already from another program the points between the "bare ground" and 1m, so my map "1" contains everything ABOVE 1m, my map 2, will contain everything after 2m, etc until reaching 30m. From these maps, I'll need to produce a map with the highest three height by pixel. Thanks


title: How to keep initial zeros when importing csv to dbase in ArcMap

text: I have a CSV table with millions of records. When I imported in into dbase via ArcMap, those inital zeros in a certain field just went away. E.g. one field show be "01001", but it displayed as "1001". How can I keep the inital zeros?


title: Generating sequential numbering for duplicate values in a field sorted by ObjectID

text: I am trying to figure out a solution to an field calculation problem. What I'm trying to do is automatically (using either ArcPy or the Field Calculator) generate sequential numbers, beginning at 1, in a new blank integer field (called 'Point_ID') for every recurrence of a value in a second field (called 'Line_ID'). The sorting of the sequential values in the second field ('Line_ID') will be based on the order in a third field ('FID'). Can anyone help me do this? I am very green when it comes to advanced field calculations and ArcPy (and Python in general, for that matter). So your explicitness is much appreciated. Another way to put it: there are duplicate values in the Line_ID field. I would like to create a new field that counts up for every duplicate occurrence in the Line_ID field, with the sort order based on the FID field. So, if there are nine values in the Line_ID field that have the value "A," the new Point_ID field will go from 1 to 9, with the order based on the sort on FID.


title: How to shift vertices of shapes in shapefile using FOSS python library?

text: I need to shift each vertice of shape in a shapefile (of any type:point, line, polygon, multiline etc.) with Python. I need to loop in all shapes, get the vertices, get it's coordinates and for each X,Y update with new values. Which Free and Open Source Software (FOSS) library should I use (ogr, fiona shapely etc.)? Any example?


title: Texture overlay in d3 map

text: D3 map can be drawn and colored as shown in this example : http://techslides.com/demos/d3/d3-world-map-zoom-points.html. I need to overlay a texture or image within the boundary of a given country. Is there a way to overlay textures in d3 map ?


title: Geoserver optional symbology with variable substitution in SLD

text: I am trying to write an SLD that optionally takes a symbology for a PolygonSymbolizer and I can't work out how to do default to no symbology if the symbology env parameter is not set. I have used the Variable Substitution example at http://docs.geoserver.org/stable/en/user/styling/sld- extensions/substitution.html as my basis.               <?xml version="1.0" encoding="UTF-8"?>     <sld:StyledLayerDescriptor xmlns="http://www.opengis.net/sld"     xmlns:sld="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc"     xmlns:gml="http://www.opengis.net/gml" version="1.0.0">       <sld:NamedLayer>         <sld:Name>Configurable Polygon Style with Defaults</sld:Name>         <sld:UserStyle>           <sld:Name>Configurable Polygon Style with Defaults</sld:Name>           <sld:Title>Configurable Polygon Style with Defaults</sld:Title>           <sld:FeatureTypeStyle>             <sld:Name>configurable_polygon</sld:Name>             <sld:Rule>                    <sld:PolygonSymbolizer>                           <sld:Fill>                   <sld:GraphicFill>                     <sld:Graphic>                       <sld:Mark>                         <sld:WellKnownName>                           <ogc:Function name="env">                             <ogc:Literal>symbology</ogc:Literal>                             <ogc:Literal>solid</ogc:Literal>                           </ogc:Function>                              </sld:WellKnownName>                         <sld:Stroke>                           <sld:CssParameter name="stroke">                             <ogc:Function name="env">                               <ogc:Literal>fill-color</ogc:Literal>                               <ogc:Literal>#1A4F57</ogc:Literal>                             </ogc:Function>                           </sld:CssParameter>                           <sld:CssParameter name="stroke-opacity">                             <ogc:Function name="env">                               <ogc:Literal>fill-opacity</ogc:Literal>                               <ogc:Literal>0.8</ogc:Literal>                             </ogc:Function>                           </sld:CssParameter>                           <sld:CssParameter name="stroke-width">                             <ogc:Function name="env">                               <ogc:Literal>fill-width</ogc:Literal>                               <ogc:Literal>0.5</ogc:Literal>                             </ogc:Function>                           </sld:CssParameter>                         </sld:Stroke>                       </sld:Mark>                       <sld:Size>                         <ogc:Function name="env">                           <ogc:Literal>fill-size</ogc:Literal>                           <ogc:Literal>10</ogc:Literal>                         </ogc:Function>                       </sld:Size>                     </sld:Graphic>                   </sld:GraphicFill>                   <sld:CssParameter name="fill">                     <ogc:Function name="env">                       <ogc:Literal>fill-color</ogc:Literal>                       <ogc:Literal>#37747d</ogc:Literal>                     </ogc:Function>                   </sld:CssParameter>                   <sld:CssParameter name="fill-opacity">                     <ogc:Function name="env">                       <ogc:Literal>fill-opacity</ogc:Literal>                       <ogc:Literal>0.8</ogc:Literal>                     </ogc:Function>                   </sld:CssParameter>                 </sld:Fill>                 <sld:Stroke>                   <sld:CssParameter name="stroke">                     <ogc:Function name="env">                       <ogc:Literal>stroke-color</ogc:Literal>                       <ogc:Literal>#1A4F57</ogc:Literal>                     </ogc:Function>                   </sld:CssParameter>                   <sld:CssParameter name="stroke-opacity">                     <ogc:Function name="env">                       <ogc:Literal>stroke-opacity</ogc:Literal>                       <ogc:Literal>0.8</ogc:Literal>                     </ogc:Function>                   </sld:CssParameter>                 </sld:Stroke>               </sld:PolygonSymbolizer>             </sld:Rule>           </sld:FeatureTypeStyle>         </sld:UserStyle>       </sld:NamedLayer>     </sld:StyledLayerDescriptor>      I want to be able to essentially test whether the 'symbology' env parameter in the WMS request has been set and apply the enclosed styling otherwise default to having no `<sld:GraphicFill>`. Does anybody have any pointers?


title: Plugin error in Qgis

text: I just now downloaded a plugin named "XY tools"from fetch plugins option . I am a user of Qgis v 1.9.0.It shows me the following error: An error has occured while executing Python code:               Traceback (most recent call last):       File "C:\Documents and Settings\15491/.qgis2/python/plugins\xytools\xytools.py", line 123, in xyToolClick         if self.layerInfo.has_key(self.layer) or self.getXyColumns(self.layer):       File "C:\Documents and Settings\15491/.qgis2/python/plugins\xytools\xytools.py", line 296, in getXyColumns         dlg = XyToolsDialog(self.XY_COMBO_MSG, layer.dataProvider().fields(), layer.name())       File "C:\Documents and Settings\15491/.qgis2/python/plugins\xytools\xytoolsdialog.py", line 35, in __init__         for field in fieldMap.values():     AttributeError: 'QgsFields' object has no attribute 'values'          Python version:     2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)]               QGIS version:     1.9.0-Master Master, 1a723b4          Python path: ['C:\\Documents and Settings\\15491\\.qgis2\\python\\plugins\\GeoCoding', 'C:\\OSGeo4W\\apps\\qgis-dev\\python\\plugins\\sextante', 'C:\\Documents and Settings\\15491\\.qgis2\\python\\plugins\\GeoCoding', 'C:\\Documents and Settings\\15491/.qgis2/python/plugins\\sextante_taudem', 'C:\\Documents and Settings\\15491/.qgis2/python/plugins\\sextante_animove', 'C:\\Documents and Settings\\15491/.qgis2/python/plugins\\sextantelwgeomprovider', 'C:\\Documents and Settings\\15491/.qgis2/python/plugins\\sextanteexampleprovider', 'C:\\Documents and Settings\\15491/.qgis2/python/plugins\\permaclim', 'C:\\Documents and Settings\\15491\\.qgis2\\python\\plugins\\GeoCoding', 'C:/OSGeo4W/apps/qgis-dev/./python', 'C:\\Documents and Settings\\15491/.qgis2/python', 'C:\\Documents and Settings\\15491/.qgis2/python/plugins', 'C:/OSGeo4W/apps/qgis-dev/./python/plugins', 'C:\\OSGeo4W\\bin\\python27.zip', 'C:\\OSGeo4W\\apps\\Python27\\DLLs', 'C:\\OSGeo4W\\apps\\Python27\\lib', 'C:\\OSGeo4W\\apps\\Python27\\lib\\plat-win', 'C:\\OSGeo4W\\apps\\Python27\\lib\\lib-tk', 'C:\\OSGeo4W\\bin', 'C:\\OSGeo4W\\apps\\Python27', 'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages', 'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\PIL', 'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\win32', 'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\win32\\lib', 'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\Pythonwin', 'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode', 'C:\\Documents and Settings\\15491\\.qgis2\\python\\plugins\\mmqgis/forms', 'C:\\Documents and Settings\\15491\\.qgis2\\python\\plugins\\pdokbaggeocoder/forms', 'C:\\OSGeo4W\\apps\\qgis-dev\\python\\plugins\\fTools\\tools']      what is the problem?


title: Changing the scale (generalisation) in ArcGIS 10

text: I have 2 polygon layers, accordingly 1:10 000 and 1:500 000. I want to generalize the first layer from 1:10 000 to 1:500 000. I tried using Simplify polygon tool, and since these scales relate to each other as 1:50, I used these options in the tool:               Simplification Algorithm : POINT_REMOVE     Maximum allowable Offset : 50 Meters      Is this the right approach, because the result does not look like what I was looking for. There is only a subtle change in detail level.


title: Comparison of Open Source Desktop GIS Packages

text: I'm curious as to what comparisons are out there between different Open Source Desktop GIS software packages. I'm aware that there has been a QGIS and gvSIG comparison on this site but I am looking for comparisons between any OpenSource Desktop GIS. A substantial list of Open Source GIS software can be found in response to the question: What are some Free and Open Source GIS Desktop packages I don't expect a single comparison with everything in it, but comparative subsets. Also any comprehensive benchmark/testing that may have been performed between any grouping of these products.


title: Missing tools (icons) QGIS

text: I lost selection and navigation tools from the GGIS 1.8 toolbar when I tried drag then to different place. I can't get them back anymore. Where they disappear? All these tools are checked but are still missing.


title: Raster processing with Python or Arc raster calculator

text: I have 4-band NAIP imagery that require the following processing steps:   1. Convert to 8 bit unsigned B/W   2. Reverse ( _note_ not inverse) the image, so that dark values appear light and light values appear dark. Moderate values should remain relatively unchanged. Could someone please point me in the right direction for a Python or Arc raster calculator approach?


title: canvec: why so quiet? where are all the users?

text: It has now been a little over 3 years since Natural Resources Canada released CanVec, a national coast to coast to coast spatial dataset at 1:50,000 or better nominal scale. In the interverning time there have been several updates, we now sit at Edition 6, and there are scheduled updates every 6 months, indefinitely. Canvec replaces the National Topographic Database (NTDB) which was both expensive, ~$240 per NTS tile if memory serves, and definitely not libre. With the NTDB we had to pay royalties for every map that went to another party, even if we gave it away. The royalty was reasonable, on the order of 25 cents each, but we had to pay in advance, in $5,000 chunks, with a minimum of $8,000. Canvec is a well modelled internally consistent dataset of good quality, better than the NTDB, and it's FREE, both libre and gratis. What puzzles me, greatly, is the virtual silence online* about this fantastic resource. And this where we come to you, dear GIS professionals. Why? Why is there so little chatter about Canvec? Is this scale of data just not that relevant south of 60°N and you're using something different? (what?) Are the naming conventions so opaque people just can't drum up the energy to figuring out what's what? Or is it in use all over and I'm just not looking in the right corners? And finally we come to the question behind the questions: Where are the resources and tools to make using this stuff easier (or even possible)? Thank you for your thoughts.


title: How to calculate Image boundary / footprint of satellite images using open source tools?

text: I'm referring to this question&discussion: How do I create a shapefile showing footprints of Rasters? I need to create the polygon outlines of several single layer raster images, not the extent/bounding box, but the area without the nodata values, as shown here: How do I create a shapefile showing footprints of Rasters?. In the answer to to post above, the _Image Boundary plugin_ was mentioned, however, I don't find it in the QGIS 1.8.0 Lisboa, running on Ubuntu? Is the tool still available? If not, is there a way to do this using opensource tools, either gdal, R, QGIS, GRASS, or similar, not with ArcMap? Thanks!


title: WMSGetFeatureInfo returns blank

text: I've been working on a WMS with a QuantumGIS server. So I've added the QGIS map as a layer to my OpenLayers map on my webpage. Now I want to use the WMSGetFeatureInfo function of OpenLayers to get information about the feature I've been clicking on. So I've followed some tutorials / examples on the web. Unfortunately it returns blank whenever i click on the layer. Here is the code:               $(document).ready(function() {              //wenn geolocation erfolgreich ermittelt wurde, setze die map auf die koordinaten              function success(position) {             lat = position.coords.latitude;             lon = position.coords.longitude;             jumptolonlat(lat, lon);              }              function error(msg) {             console.log( typeof msg == 'string' ? msg : "error");         }              if (navigator.geolocation) {             // getCurrentPosition ruft die Funktion success auf und übermittelt die Position Werte             // error wird ausgeführt wenn es einen Fehler beim ermitteln der Position gibt             navigator.geolocation.getCurrentPosition(success, error);         } else {             alert("GeoLocation API ist NICHT verfügbar!");         }         //document.getElementById("lat").value = "";         //document.getElementById("lon").value = "";         var options = {             projection : new OpenLayers.Projection("EPSG:900913"),             displayProjection : new OpenLayers.Projection("EPSG:4326"),             units : "m",             //resolutions : [156543.03390000001, 78271.516950000005, 39135.758475000002, 19567.879237500001, 9783.9396187500006, 4891.9698093750003, 2445.98490468750],             //maxResolution : 156543.0339,             maxResolution : 0.17578125,             minResolution : 156543.03390000001,             maxExtent : new OpenLayers.Bounds(-20037508.34, -20037508.34, 20037508.34, 20037508.34),             controls : [new OpenLayers.Control.Navigation(), new OpenLayers.Control.LayerSwitcher(), new OpenLayers.Control.PanZoomBar(), new OpenLayers.Control.Attribution()]         };         // Karte initiieren         map = new OpenLayers.Map('mapdiv', options);              // Layer hinzufügen         layerMapnik = new OpenLayers.Layer.OSM("Mapnik");         var wms = new OpenLayers.Layer.WMS("Karte", "http://localhost/qgis/qgis_mapserv.fcgi.exe", {             layers : 'gk-1000',             format : 'image/png',             map : 'C:/mymaps/karte.qgs'         });         map.addLayers([layerMapnik, wms]);         // Zentrierung der Karte beim initiieren         var centerLonLat = new OpenLayers.LonLat(9.869824084639358, 51.25325741806755).transform(new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject());         if (!map.getCenter())             map.setCenter(centerLonLat, 6);              // Click Control zur Karte hinzufügen         click = new OpenLayers.Control.Click();         map.addControl(click);              // Aktivieren (kann auch mit click.deactivate() wieder deaktiviert werden)         click.activate();              var popup = null;         info = new OpenLayers.Control.WMSGetFeatureInfo({             url : 'http://localhost/qgis/qgis_mapserv.fcgi.exe',             title : 'Identify features by clicking',             queryVisible : true,             layers : [wms],             eventListeners : {                 getfeatureinfo : function(event) {                     if (!event.text.match(/<body>\s*<\/body>/)) {                         if (popup) {                             popup.destroy();                             popup = null;                         }                         popup = new OpenLayers.Popup.FramedCloud("popup", map.getLonLatFromPixel(event.xy), null, event.text, null, true);                         document.getElementById('ergebnis').innerHTML = "Ergebnis: " + event.text;                         map.addPopup(popup);                     }                 }             }         });         map.addControl(info);         info.activate();          });      Any help would be appreciated.


title: Free application or web service for reverse geocoding

text: I'm looking for a fast and free web service or local application (this it would be better, a local service with database of street, numbers, city, etc.) for reverse geocoding. I could use google but is limited to 5000 query per day... Thank you in advance ;) EDIT: What I m trying to do is convert positions acquired from a gps (1 per second) in addresses. These positions could be more than 7000 per track, and if I have 10 track for example, the mole of data to "translate" becames "huge"..


title: How to resize a polygon in km?

text: I have a polygon on my map and I want to resize it in that way I add X km to all of the vertexes. Any idea?


title: How can I combine addr:* tags in JOSM to create the addr:full tag?

text: I have about 400 addresses that I need to update in OSM. They have already been conflated using Josh Doe's Conflation plugin, http://wiki.openstreetmap.org/wiki/JOSM/Plugins/Conflation, which my version of JOSM no longer supports for some reason. I have updated the **addr:street:name** and **addr:street:type** for the dataset in JOSM, but now I need to update the **addr:full** tag. Is there a way to perform a calculation in JOSM that would concatenate something like this: **addr:full** = **"addr:housenumber"** + **"addr:street:name"** + **"addr:street:type"** If there is a way in iD, Potlatch2 or Merkaartor that would be great too.


title: Very high resolution stereo-imagery

text: I'm looking for very high resolution stereo imagery to perform a DEM extraction. The only place I can find sample data is http://www.digitalglobe.com/product-samples but the form does not work for me. Any idea where I can find the data I'm looking for? QuickBird 2, IKONOS, WorldView, GeoEye, RapidEye, or some other very high resolution source and of mountainous terrain would be ideal.


title: How can I make my labels scale dependent?

text: I present labels from a textlayer and make the visibility dependant from the zoom. However when I save the project and reopen it, the zoom is grayed out. The project does not keep the inserted zoom level.


title: Convert shape file to kml without losing symbols

text: I am using QGIS 2.2 and want to convert a shape file containing point coordinates with userdefined svg-symbols to a kml file. Everything seems to work. The orginal file is in WGS 84 / UTM zone 32N and I convert it to kml EPSG 4326. When I open the new kml file in GoogleEarth the position of the points is correct but as a symbol I get GoogleEarth-Pins. Is it possible to include the svg-symbol and the label as well in the kml file? Thanks!


title: How to serve both WFS and WFS-T from one server?

text: In GeoServer, how can I serve two types of WFS from one server? Setting the type of the service is done once in global settings. Is it possible to have WFS from workspace A and WFS-T from workspace B?


title: Attribute Inspector for points - JavaScript API

text: I am attempting to get an application working where editors can click on a point feature and edit a form. I need the form to have a save button that when clicked will apply the edits to the feature. I am selecting the point/s using an envelope. Is there a way I can tell what feature I need to apply the edits to on my save click event? I would like to store the edits to a feature as fields are changed and then only save the edits to the feature class with the save button. I am trying to use this example to understand the edit saving and then the feature save. Attribute Inspector Here is what I am working with right now. I am grabbing the point features within an envelope:               function initEditor() {               var featureLayer = layersFS[1].layer;                            var selectQuery = new Query();           map.on("click", function (evt) {                    var centerPoint = new esri.geometry.Point                (evt.mapPoint.x, evt.mapPoint.y, evt.mapPoint.spatialReference);               var mapWidth = map.extent.getWidth();                    //Divide width in map units by width in pixels               var pixelWidth = mapWidth / map.width;                    //Calculate a 10 pixel envelope width (5 pixel tolerance on each side)               var tolerance = 10 * pixelWidth;                    //build query filter               query = new esri.tasks.Query();               query.returnGeometry = true;               query.outFields = ["*"];                    //Build tolerance envelope and set it as the query geometry               var queryExtent = new esri.geometry.Extent                       (1, 1, tolerance, tolerance, evt.mapPoint.spatialReference);               query.geometry = queryExtent.centerAt(centerPoint);               featureLayer.selectFeatures(query, FeatureLayer.SELECTION_NEW, function (features) {                   if (features.length > 0) {                       //store the current feature                       for (var i = 0; i < features.length; i++) {                           updateFeatures.push(features[i]);                                                 }                       //map.infoWindow.setTitle(features[0].getLayer().name);                       map.infoWindow.show(evt.screenPoint, map.getInfoWindowAnchor(evt.screenPoint));                   } else {                       map.infoWindow.hide();                   }               });           });      Then I am planning to update the specific feature with the save button:                saveButton.on("click", function () {               //need to figure out how to get the current feature with the save click,                //the envelope may return 1 to many records and I need the current feature               updateFeatures.forEach(function(entry) {               updateFeature.getLayer().applyEdits(null, [updateFeature], null);                           }            )});                ///////Need to apply similar logic as save button click but this will update the current feature in the updateFeatures[]           attInspector.on("attribute-change", function (evt) {               //store the updates to apply when the save button is clicked                updateFeature.attributes[evt.fieldName] = evt.fieldValue;           });      What is the best way to determine what my current feature is in the info window that is being edited? Thanks for any advice!


title: How to query points that fall outside their designated regions?

text: I have a shapefile with 250 points. Each point has a `region` attribute. In another shapefile I have all the possible regions with their names (it is polyline geometry). Unfortunately the coordinates of the points are sometimes wrong. I would like to query which points don't lie in the region specified by their `region` attribute.


title: how trace polyline between two points?(Esri)

text: How i can draw a line between two points...e.g.: I have initial point and endpoint and i need trace a line between this points...how i can do that?


title: Where can I find annual temperature & relative humidity data of every town in India?

text: Need the historical data of temperature & relative humidity of a place called Bhiwadi in Rajasthan state of India


title: How to add the Delimited Text Plugin for QGIS 1.8?

text: I can't find the Delimited text Plugin- everything keeps telling me it's in the core build, but which core build? the Mac OS variant comes in multiple installers so without programatically unpacking all of those files to discover it where the hell is it? Can anyone please point me in the right direction, or post a link to the file; so that I can get it installed on my macs. Or alternatively state that the plugin is NOT available for the mac version [unhappy face]


title: Question about subs in MapBasic

text: For some reason my sub isn't 'linking' in MapBasic, I'm pretty new to this, but why is the sub assuming it begins on line 110 rather than where it actually says sub? ![enter image description here](http://i.stack.imgur.com/hTiKW.png) ![enter image description here](http://i.stack.imgur.com/xg2sK.png) Thanks, Alan


title: How to build Android spatialite?

text: I have been trying to build spatialite for Android for a few days now without success. I have even tried the recent documentation from the Google project page here: https://code.google.com/p/spatialite-android/. Build documentation is is quite vague and there are not many resources on the interweb to help me out. Just wondering if this project has seen adoption or is simply dead. If its dead its too bad... it was a good initiative.


title: What is the easiest way to get points Info from OSM

text: I need to parse some OSM XML data to get some nodes with their tags. Does a program or script already exist to do this? The result can be any format that I can parse to DB, e.g. XML, CSV, txt file with tabs. The solution must work with ~400GB file (full planet), and it would be extremely good, if it can split the output files by size. A solution which works on Debian is preferred, although I will take anything which already exists.


title: How to calculate smallcircle on sphere

text: I use this python code to calculate a greatcircle: https://code.google.com/p/pyroms/source/browse/trunk/pyroms/greatcircle.py?r=39 After initializing the greatcircle I can calculate points on the line. I would like to do the same with an arbitrary smallcircle. I would like to define the smallcircle using its center location (in latlon) and one location on the line (also latlon), so not using a diameter. Say as an input I would use a certain point spacing in meters, how would I then proceed to calculate the other locations on the smallcircle line? Thanks!


title: How to randomize colors for large datasets in QGIS new symbology?

text: Someone has already asked how to randomize colors by class in QGIS new symbology, but the answer there does not work with datasets containing more than 100 unique values since the Random option only allows you to create only 100 colors. Using that technique, the features get categorized using the equal breaks method making the color the same for many _consecutive_ features instead of assigning one per feature and then starting over with the 100 colors again after they have all been used once. My goal is to differentiate features in a polygon layer with many similarly located and named features, e.g., Woodlands 1 and Woodlands 2 are next to each other and so need to be different colors. Is there a way to make sure that each consecutive feature is a different color from the immediately preceding and following feature? In the old symbology, you can just use "Randomize" to achieve this. Is there a way to do it in the new symbology? Thanks for any ideas or help!


title: How to rescale data ranges in ERDAS?

text: I have an Normalized Difference Vegetation Index(NDVI) and Enhanced vegetation index (EVI) layer and I want to re-scale their data ranges to value ranges from 0–to–1, in order to optimally stretch the contrast (information contained in the image). I have analysed the histogram of the distribution of grey values and need to omit the lowest 0.5% and the highest 0.5% of the data, however I am not too sure how to go about doing this. I tried doing this in excel (manually deleting 0.5% min and 0.5% max values) and then completing the step in ERDAS model builder. The formula used in model builder (NDVI- new min)/(new max-new min), however my results were unsuccessful.


title: Error while importing CSV File using PyQGIS

text: Currently I am trying to import a CSV file using the python extension in QGIS and I have run into a problem that I can not seem to solve. I am trying to do a basic import of a CSV file from a python code I created (See below) and every time the file loads the values in my attribute table become re-arranged (see image). When I load the file from the "Create a Layer from a Delimited Text File" option it loads without a problem. Here is a sample of my code:               uri = "file:///D:/MLB Stadiums/Test_R2/Player_Files/Espinosa_2014.csv?            delimiter=%ss&xField=%s&yField=%s" % (",", "field_3", "field_4")        vlayer = QgsVectorLayer(uri, "Espinosa", "delimitedtext")        caps =vlayer.dataProvider().capabilities()      Here is a comparison of what it should look like after import and what it does look like. ![enter image description here](http://i.stack.imgur.com/1b9HP.png) Does anyone know why it is rearranging the values?


title: Leaflet.js and CartoDB.js event conflict

text: Little problem here. I’ve got a layer from a CartoDB viz that I’m listening to the featureclick event. When a feature is clicked I highlight it by applying cartoCss. However, I want to unhighlight it when I click on an area of map with no feature. I did this using the Leaflet map.on(‘preclick’) event. I chose preclick so it fires before the featureclick and applies the default cartoCss first. Problem is the two events don’t play nicely. The map click event always fires, but the layer featureclick event only fires some of the time. So it sometimes takes two or three clicks before the feature is highlighted – even when the cursor has changed to a pointer so the map knows there is a feature there. What is the best way to unhighlight a feature after a mouse click on the basemap? Can I ignore my map click event when a layer feature is clicked? Can I change which event fires first i.e. layer featureclick before map click? Code:                  map.on('preclick', function(e) {               //remove project highlighting and close popup               J.Dashboard.highlightProject(null);        });               layer.getSubLayer(0).on('featureClick', function(e, latLong, pos, data, layerIndex) {               J.Dashboard.highlightProject(data.name);        });             highlightProject : function(projectName) {               //Update the style               J.Map.projectsLayer.getSubLayer(0).setCartoCSS(J.Style.getSelectStyle(projectName));      } Cheers, James


title: R- alphahull - duplicate points error help

text: I came across R function -alphahull- and believe it is the best option to turn gps points into a polygon (3 or 4 days worth of gps points across residence/neighborhood for each of 300 subjects). However, I am just starting to use R and I would appreciate any help with removing duplicate points, as I am getting the following error: "Error in tri.mesh(X) : duplicate data points" after running:               library(sp)     library(alphahull)     library(maptools)          #define projection     project2<-"+proj=utm +zone=17 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"     # read point shapefile     data.shape <- readShapePoints("D:/clayera/UWO/msthesistake2/tesisGIS/D1001_GPS_WalkExtSel.shp",proj4string=CRS(project2))     plot(data.shape)     summary(data.shape)     # extract coordinates     x.coords <- coordinates(data.shape) #this way we get duplicate coordinates     # alpha-shape: 100 meter threshold     x.as <- ashape(x.coords[,1], x.coords[,2], alpha=100)      Apologies for such a basic question and thanks for any comment/help.


title: Using android phone to plot GPS points?

text: I was wondering if I could use my android phone (HTC Wildfire) to plot gps points and export these into GIS software (I have QGIS and MapInfo). Thanks, Alan


title: How to create a network with all the track points

text: I have created a sample network through the use of gps visualizer. Than I used the ogr2ogr function to import this data in my database. Now I have a table which holds all the track points of these tracks. I would like to build a network between all these points so that at the end I can calculate the shortest path from a point to any point. In all I have 28 track points where 26 of these points are distinct. I would like some suggestions how can I determine the source and the target of each point in order to create a network. Any suggestions please?


title: Using PHP and JavaScript to dynamically create shapefile or KML

text: While I have used GIS software (ArcGIS and just started playing with QGIS at home), I have never looked into the code or file structure of shapefiles. I am working on a project for my work, a web-based application running on a LAMP server and written in PHP and JavaScript. The idea is a search box that will be used to search the "geographical hierarchy" of a specific location. This would allow my users to type in a city, and it would return the state, country, and continent in which the city is located. Google can provide this information and more, but this will be used on a closed network so that's not an option. I am using a MySQL database containing all continents down to the city level. I still have a little way to go, but I know what needs to be done to finish this part. I recently thought it would be nice to add the ability for my users to enter a geo-coordinate and get the same result as if they had entered the city the coordinate is located in. Right now I do not have any coordinates in my database, but I know I have to reference them somehow. I believe there is an SDE database on this network, but i have not tried to access it or see what it would take to pull data from it. I would like to just reference it to get the geo data rather than duplicating the data in my database, if possibly (to keep my database as small as possible). Another feature I would like to provide (now, if it is not too difficult, or as a future update) is the ability for my users to dynamically build a KML (as a lot of users would be using GoogleEarth) or a shapefile based on the users' specific request (single country/multiple countries, states or cities within countries, etc). To build the KML, I know that I would need a list of decimal lat/long for the shape of the polygon. I am not sure what would be needed for the shapefile. My reasoning for this feature is to allow GoogleEarth users the ability to quickly/easily generate a KML and, for GIS users, my tool could just be another front end to the SDE to generate specific shapefiles. So, specifically: 1. What do I need to do to be able to look up a location based on a coordinate located within that location? 2. What is the best way to dynamically create a KML? 3. How can I dynamically create a shapefile keeping the data from the original database? Since this is a closed network, I am limited to the libraries and scripts that other people have already added. I am hoping that i can code most of this by hand or transcribe it, and I would be most comfortable using PHP and/or JavaScript. However, if there is a specific script or library that I would need, there is a process to get stuff moved to this network, but it is not a simple, and then it has to be approved, which is its own process. As this is a work in progress, I am open to suggestions to improve what I have done so far or what my thinking for the rest of my program. I apologize for this being so long, but I wanted to make sure I included my entire thought process to provide the background for my questions. Thank you in advance for your input.


title: meaning of createObject(Point.getClsid())

text: What actually does the `createObject(Point.getClsid())` do?Does it create a new Point Object,and return the class id of the `Point`.Does it mean that each point created will have the same class id. I went through ESRI documentation but all it said was Identify class ID for an object.


title: How To Add Values to an Attribute Field Based on a Selection - Python

text: Alright another problem that I have. I am trying to give a priority ranking to lots based on various criteria. Basically I have 7 areas that are criteria, if lots overlap each other the criteria is added up and given a sum value (so the highest number could be 7). Now I am trying to create a python script that basically goes through my criteria and adds a value to a field added to an attribute table For example, I would do a select by location for lots that are within 30m of wetlands - if YES then the field "CRITERIA" would be given a 1, if no it would be a 0 this same selection idea would be done for the rest of my criteria - then what I need is to have it total up the 1's that are in each row into a separate "RANKING" field.


title: How do I find least cost path over a DEM based on slope magnitude and direction?

text: What sort of algorithms can I use to find a least cost path through a DEM where it costs more to traverse a cell going downhill or uphill than it does to follow the contour of the land? All the examples I've seen involve creating a least cost surface where the _magnitude_ of the slope is part of the cost - but the _direction_ isn't. Additionally, the examples assume the cost to travel through a cell does not depend on the direction of travel.


title: Including two different maps in one QGIS Print Composer

text: How do I use print composer in QGIS 2.0.1 to produce a final product with two different maps on the same sheet? I can't figure out how to use the 'add a map' button to add anything other than a duplicate copy of what I already have. Your help is much appreciated!


title: Create a raster with georeferenced information in R

text: Referring to my previous question, how do I "geocode"/convert - using R - a generated raster image to a format that can be shown as a layer on a map? In short, I'm aiming to recreate the results of the commercial SpatialKey (screenshot) software. I'm using a GeoServer instance to host the data for the maps, but how can I incorporate the required projection and coordinate info into the image so it will align correctly with its true geographical position. The density heatmap was produced using spatial data.


title: Excel Nad83 table to gpx-Total noob

text: I have some waypoint locations that I have entered into a spread sheet in Nad83 the format of table looks like this NAME, -, -, UTM ZONE, Easting, Northing, -, Comment 1184236, 1, CP, 15, 451596, 5664617, 20, No post ["-"'s Are information only needed in the spreadsheet] I have to get these coordinates into some software that will let me export the points as .GPX. I use to do this using a wizard in mapinfo but I no longer have the software. Basically I'm looking for something that will let me define the coordinates as NAD83 and let me select which columns are to be used for Name/Zone/Easting/Nothing.


title: Fleet tracking software options

text: Before I sit down and make the mistake of start writing my own, I wanted to know my options for good, **modern** , fleet tracking options that are not too expensive. What I am looking for is something that doesn't look like it was built in 1995 and that will not cost me an arm and a leg as I scale up. Brownie points if the project is Open Source, but using a hosted solution will do just fine provided that it is not too expensive.   * What are people using to track on the field assets in realtime?    * If price was not an issue, what would you use? **Update:** Hardware is not much of an issue for this project. There are several spare servers with plenty of bandwidth and redundant drives connected to the Internet 24/7. As far as client hardware, there are several of-the- shelf devices with GPS and connection to the internet, not to mention that some of these devices will be modified smart phones. My concentration is now more related to the server and web front end software I will use for this. **Update 2:** I was asked for a use case - so I will give a simplified version. A company that does deliveries wants to know where their fleet is. There are several managers and each of them manage a chunk of delivery trucks. Each of the managers (there are about 100) is only interested in their delivery trucks (each has around 20) and wants to see them in realtime at - the same time -on a web map. Each driver wants to see where their partners are, too - but this is secondary. Privacy issues are just that each manager will be able to see their own drivers only. I have any other data required (streets, geofences, etc - but that is a separate discussion - I am looking only for this basic functionality). Yes, some time-based reports about where the trucks are/were will be necessary.


title: How can I translate Modis HDF To GeoTiff

text: I've a MODIS04 L2, and I want to extract band 8:Image_Optical_Depth_Land_And_Ocean. I use Quantum GIS to open this band, because Quantum GIS see this is not coordinates so this allow me to select projection as below: +proj=utm +zone=50 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +towgs84=0,0,0. All I know is that I want to translate to WGS84. This open this band seems good but when I open this band with ENVI this is better (differents output images). But I see the Modis Conversion Toolkit use same warp parameters as Quantum Gis expect it has a parameters ("Numbers of Warp Points"). I think that I miss something with Quantum Gis and GDAL too. This is Envi Conversion Toolkit ![enter image description here](http://i.stack.imgur.com/zwRQL.png) This is open file HDF with QuantumGis ![enter image description here](http://i.stack.imgur.com/AqDbW.png) ![enter image description here][3]


title: Show geoJSON properties in jQuery Mobile Panel

text:
First off I'm a novice at JS so please bear with me if this is a simple task to complete. I'm currently using the stock LeafletJS popup bubble to display the properties from my geoJSON. However, I would like to utilize the jQuery Mobile Panel to display the properties from my geoJSON file when it is clicked on, it is similar to the OpenLayers example here. Thanks in advance and below is my code for loading the geoJSON!               $.getJSON(dict[geoJSON Variable Goes here],           function(data){             L.geoJson(data, {               style: style,                onEachFeature: function (feature, layer){                  if (feature.properties.SOMETHING == 0){                      return layer.bindPopup("" + feature.properties.Name + " LABEL
LABEL");}                     else if (feature.properties.Attacks == 1){                              return layer.bindPopup("" + feature.properties.Name + " LABEL
" + feature.properties.SOMETHING + " label");}                     else {                           return layer.bindPopup("" + feature.properties.Name + " LABEL
" + feature.properties.SOMETHING + " label");                                     }                                 },                             }).addTo(GROUPLAYER FOR LAYER ON/OFF);                         });


title: ArcGIS map export is exporting some features but not others

text: I feel this is pretty simple but I haven't had much success searching for an answer and I'm also running out of time so I'm resorting to asking here. I am trying to export a map with the exact same settings as I always do, but it is choosing not to export the more complicated features such as contours, water and our vegetation inventory while successfully exporting roads, annotation and all the layout components. The problem does not seem to be confined to a certain geometry type. I get the same result no matter the format or settings. Some of the data is on a local server while most of it is from a shared corporate server. The main culprit is probably the fact that I'm using data driven pages, but I get the same result when I turn them off. Any ideas? It's not even attempting to export them as it completes in a few seconds while normally I would expect it to take longer. I unfortunately do not have any graphics programs to mess around with. Thanks in advance.


title: GeoExt (ExtJS + OpenLayers) not displaying Google Map

text: I've just started exploring GeoExt. It displays all the map controls, Google Logo and "Terms of Use". But the map is not rendered. I'm using firebug, and no JavaScript errors are generated. I'm using ext 3.4, OpenLayers 2.11, and GeoExt 1.1.                   <link rel="stylesheet" type="text/css" href="libs/ext/resources/css/ext-all.css">         <script type="text/javascript" src="libs/ext/adapter/ext/ext-base.js"></script>         <script type="text/javascript" src="libs/ext/ext-all.js"></script>         <script src="libs/openlayers/lib/OpenLayers.js"></script>         <script type="text/javascript" src="libs/geoext/lib/GeoExt.js"></script>         <script src="http://maps.google.com/maps/api/js?sensor=false"                 type="text/javascript"></script>              <script type="text/javascript">              Ext.BLANK_IMAGE_URL = "ext/resources/images/default/s.gif";         var app, items = [], controls = [];              Ext.onReady(function() {             app = new Ext.Viewport({                 layout: "border",                 items: items             });         });              items.push({             xtype: "gx_mappanel",             ref: "mapPanel",             region: "center",             map: {                 numZoomLevels: 20,                 controls: controls             },             extent: OpenLayers.Bounds.fromArray([                 -122.911, 42.291,                 -122.787,42.398             ]),             layers: [new OpenLayers.Layer.Google(                     'Google Streets',                     {isBaseLayer: true}             )]         });         controls.push(             new OpenLayers.Control.Navigation(),             new OpenLayers.Control.Attribution(),             new OpenLayers.Control.PanPanel(),             new OpenLayers.Control.ZoomPanel()         );              </script>     </head>     <body>     </body>      Any assistance will be appreciated.


title: How to use GeoServer Web Processing Service (WPS) with OpenLayers and WFS?

text: I wanted to try creating WPS service in geoserver on web application. I wanted to know reference links or some help regarding using WPS on web.


title: pgrouting with TRSP in LINUX

text: > **Possible Duplicate:**   >  Why are start and end segment not partially skipped? Im using turn_restrict_shortest_path algo for routing. For start and end segment it is not skipping the partial part of segment as im giving the point value on edge (i.e. 0.5) for start nad end segment. Is that working properly or bug? Need help. Regards, Amol.


title: Problem using qgis when I upgraded to osx 10.8.1

text: I have been running qgis in Mountain Lion for a few weeks now, but when I upgraded to OS X 10.8.1, qgis would no longer work. I get an error saying: **QGIS cannot be opened because of a problem.** Check with the developer to makes sure QGIS works with this version of OS X. You may need to reinstall the application. Be sure to install any available updates for the application and OS X. I have tried re-installing qgis (I have even tried installing the nightly build) and I have also updated all of the packages that qgis requires. Has anyone had a similar problem or more importantly, does anyone have a solution? Thank you!


title: Can I run tools from QGIS using shell commands?

text: Is it possible to script rutines in QGIS (python) and then runing them remotely? To illustrate, can I through a shell (ssh or telnet) ask to QGIS to transform a shapefile into a gjson format? Mostly, I'm thinking in the interaction of QGIS with another software like R for example, where I can generate command lines. (I'm not talking about QGIS internal command shell)


title: How do I add a polyline segment created with Points to Lines tool to a feature class and apply that feature classes attributes to the line segment?

text: I am a new ArcMap user transitioning from using AutoCAD. I have created lines from a point shapefile and have broken them into segments. These segments need to have a specific attribute table attached to them with unique ID's. I have created all of the line segment types as new polyline feature classes and have the appropriate attribute table correctly formatted for each line type. They are empty at this point. What I need to do is to take the segments that I have and move them into the various feature classes so each segment has the 'classes attribute info. The ultimate goal is to have these segments correctly attributed for submission to a government agency. I hope this description makes some sense. I know that I can 'create features' and trace these segments but I was hoping to find a more efficient way.


title: Unable to create WFS successfully

text: I know there has been lots of WFS assistant and I did follow this one word by word, which I found was the best in the entire gis stack but still no clue why I am unbale to render a WFS on base map. Here is the code,               var map = new OpenLayers.Map("map-id");                 var bio = new OpenLayers.Layer.WMS(                   "bio11",                   "http://localhost:8080/geoserver/iPlant/wms",                   {layers: 'iPlant:bio11'}                 );                 var imagery = new OpenLayers.Layer.WMS(                     "Global Imagery",                     "http://maps.opengeo.org/geowebcache/service/wms",                     {layers: "bluemarble"}                 );                      /*var QuebecSpecies = new OpenLayers.Layer.WMS(                   "AllQuebecSpecies",                   "http://localhost:8080/geoserver/iPlant/wms",                   {layers: 'iPlant:AllQuebecSpecies'},                   {'opacity': 0.4, 'isBaseLayer': false, 'visibility': true}                 );*/                      var QuebecSpecies = new OpenLayers.Layer.Vector(                 "AllQuebecSpecies",                 {                     strategies: [new OpenLayers.Strategy.Fixed()]                     , projection: new OpenLayers.Projection("EPSG:4326")                     , protocol: new OpenLayers.Protocol.WFS({                         version: "1.1.0",                         url: "http://localhost:8080/geoserver/wfs",                         featurePrefix: 'iPlant', //geoserver worspace name                         featureType: "iPlant:AllQuebecSpecies", //geoserver Layer Name                         featureNS: "http://geoserver.org/iPlant", // Edit Workspace Namespace URI                         geometryName: "the_geom" // field in Feature Type details with type "Geometry"                     })                 });                      map.addLayers([imagery, bio, QuebecSpecies]);                 map.addControl(new OpenLayers.Control.LayerSwitcher());                 map.addControl(new OpenLayers.Control.MousePosition({element: $('location')}));      As you can see in the code i have commented out the part where I am using the same layer as WMS, where it works perfectly fine, but then it does not overlay on the base map (meaning Quebec layer is a point shapefile had so the rest of the space is white background which covers up the base raster, that is why I have set the transparency to be 0.4). Help me out guys, really tired of finding and trying and having no luck. **UPDATE** I moved the shapefile folder of mine to www folder and also changed the store location folder to www. Workspace : iPlant Store Name : SpeciesShp FeatureLayer : AllQuebecSpecies Is there some other way to add shapefile to Geoserver for using it as a WFS? **UPDATE 2.0** Since there were suggesstions of maybe my html file is not in the right folder, I now have my HTML file in /var/www/test.map/html as well as /var/lib/tomcat6/webapps/test.map.html. Still no luck, WFS does not show up :( This is the error I am getting in develoepr tools,               XMLHttpRequest cannot load http://localhost:8080/geoserver/wfs. Origin http://localhost is not allowed by Access-Control-Allow-Origin.      **UPDATE 3.0** Ok thanks to more suggesstions now I have my map.html in /var/lib/tomcat6/webapps/geoserver/data/www/map.html. Now I can still work with WMS but now if i try to use WFS, the error I can see in the developer tool is               Uncaught TypeError: Cannot read property 'Fixed' of undefined      Also the map container itself disappears. HELP!!!


title: Finding intersection density for weighted lines connecting points across the continental USA

text: I'm working on an one-time problem that I'm hoping to solve with a computer program using free GIS libraries (preferably for VB.Net, C#, or Python... Java or C/C++ could work too). Alternatively, I wouldn't mind using a tool like the free MapWindow GIS. _The problem:_   * Given knowledge about the ultimate origin/destination for passengers flying, draw lines weighted by passenger volume connecting origin and destination airports. (i.e. passengers flying from San Francisco to New York can travel via many routes, but we only care about the total volume of passengers flying between the two cities).   * Find intersections for all the lines. Each intersection will have a 'weight' - the sum of the passenger volume on the two intersecting lines.   * Given 3-digit ZIP codes, calculate a density for each ZIP code by summing the weights of each intersection within X miles of the centroid. The idea here, being to find good natural locations for hub airport given current passenger volume. _Tools/Accuracy:_ I found Shapely for Python that looks useful. It includes a lot of spatial analysis tools (like centroid, buffer, intersection, etc.) that would make solving this problem fairly easy. However, it works in the Cartesian coordinate system. It seems like that would be rather inaccurate on a continental scale. It seems like, from an accuracy standpoint, it would be best to do spatial analysis in the geographic coordinate system (lat/lon on a datum) itself rather than re-projecting the data. So my questions:   1. Since we only want a general picture, is the Cartesian coordinate systems accurate enough?   2. Can I do spatial analysis directly in a geographic coordinate system/datum?   3. What tools exist that would allow me to convert the data (if necessary), and do the required spatial analysis on it (with minimal setup effort since this is an one-time thing). (Feel free to re-tag as appropriate... thanks!)


title: How to configure custom scale in Leaflet?

text: I have a custom WMS tileset generated using ESRI ArcGIS Server. I have successfully configured an OpenLayers page to display this basemap. The code is as follows:               var map = new OpenLayers.Map('map', {         projection : "EPSG:3857",         resolutions : [ 132.29193125052919, 92.604351875370426, 67.733468800270941, 33.866734400135471, 16.933367200067735, 8.4666836000338677, 4.2333418000169338, 2.1166709000084669, 1.0583354500042335, 0.52916772500211673, 0.26458386250105836 ],         maxExtent : new OpenLayers.Bounds(-8352920.48698, 4887220.912586, -8082558.309632, 5067712.988343)     });          var basemap = new OpenLayers.Layer.WMS("basic_wm", "http://myhostname/geowebcache/service/wms/", {         layers : 'basic_wm',         format : 'image/jpeg'         }, {         tileOrigin : new OpenLayers.LonLat(-20037700, 30241100)     });     map.addLayer(basemap);      I am trying to reproduce the equivalent in Leaflet, but have so far been unsuccessful. I see that I can override the scale function and set a max zoom, but I am not certain how to override the scale function. I'm also wondering if I shouldn't be able to provide my DPI, tile origin, max extent, etc? Here is what I have so far in Leaflet:               var map = new L.Map('map', {          center: new L.LatLng(40.693726386087945, -73.9856945509795),          zoom: 1     });          var basemap = new L.TileLayer.WMS("http://myhostname/geowebcache/service/wms/",          {             layers: 'basic_wm',             attribution: "Copyright 2012 The City of New York",             maxZoom: 10     });     map.addLayer(basemap);


title: Linux Tool to visualize CityGML file

text: I am looking for a linux program that allows me to visualize CityGML models (3D). I have tried a whole bunch so far (most of them via wine but also native ones) but I couldn't get a single one of them running. I am using Ubuntu 13.04 and would appreciate if somebody could answer from personal experience instead of just posting a bunch of links.


title: how to keep the ids of waypoints of a shape file when i convert it to a gpx IN QGIS?

text: I tried to convert a shape file to a gpx in QGIS, and I do not know how to keep the ids (attributes! Any solution? Thanks a lot!!


title: Semi-circle representation in legend - ArcGIS!

text: I have a little problem that needs solving. How to represent the legend in ArcMap as in the picture above? ![legend](http://i.stack.imgur.com/fuPjG.png) Or as semi-circle instead of complete circles?


title: ArcGis 10 - create polygons as % of original

text: Could someone help me with the following? I have a polygon shapefile and need to create new polygons, within the originals, as % of original area. I used buffer by field (Field= new area as %) but the buffers are created in circles around the polygons. What am I doing wrong? Could you suggest smething else?


title: Where can I find shape file (free/open source) for all country . at admin 2/3 should cover major districts/towns

text: I need to download shapefiles for all country from which i need to genereate topojson for around all countries. currently I am using data from naturalearthdata.org which only provides Admin 0 ,1 for all countries. I need for Admin 2 ,3 which cover towns/cities/districts level details. I have tried diva-gis.org i had to download country by country from dropdown which is taking lot of time. 1 . Is there any site i can download shape file for admin2,3 for all countries as single shape file or country specific shape files .(free/open source) PS:I will integrate this map/ topojson inside d3.js /java /j2ee project.


title: Spatial references do not match when editing geometric network

text: I have been editing a geometric network. Since yesterday, one thing changed. I am now not able to edit the pipes in the geometric network. I get a "spatial references do not match" error. Nothing else has changed since yesterday. The mxd is the same, the database is the same, all the spatial references of the mxd, network dataset and network feature classes look identical. Oddly, I am able to move the manholes without getting an error. Does anyone know what could cause this sudden behavior and how to fix it? I am on ArcGIS 10.2 and don't have attribute assistant installed. Thanks for your answers.


title: Move shapefile from one folder to another using python

text: I'm working with shapefiles and using the ArcPy library. My script is creating a new shapefile based on a join and I want to know how to move the newly joined-shapefile from the environment location to a new folder. Any ideas on how to do this?


title: Can't solve "OperationNotSupported" error on WFS GetFeature request

text: I'm working on a web map viewer using OpenLayers, ExtJS, GeoExt and GXP. I've added a LayerTree (gxp_layertree) that shows layers from various WMS servers from different GeoServer versions. I've also configured a FeatureManager (gxp_featuremanager) with 'autoLoadFeatures = true' and a FeatureGrid (gxp_featuregrid). My problem appears when I click on a layer tree node and a WFS GetFeature request fires: the response is always "OperationNotSupported" "No such operation wfs null WfsDispatcher". I've checked the xml body content with the embedded 'demo requests' in GeoServer and it's Ok. That is:               <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.1.0" maxFeatures="100" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">       <wfs:Query typeName="feature:mce_prov_escuelas" srsName="EPSG:4326" xmlns:feature="mce">         <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">           <ogc:BBOX>              <ogc:PropertyName>the_geom</ogc:PropertyName>              <gml:Envelope xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:4326">                 <gml:lowerCorner>-62.406926496825 -26.770412516161</gml:lowerCorner>                 <gml:upperCorner>-57.557744198657 -22.49170214764</gml:upperCorner>              </gml:Envelope>           </ogc:BBOX>          </ogc:Filter>       </wfs:Query>     </wfs:GetFeature>      It's a POST request and here's a copy of the headers (from Firebug):               REQUEST HEADERS     Accept      text/html,application/xhtml+xml,application/xml;q=0.9,*/*;q=0.8     Accept-Encoding gzip, deflate     Accept-Language es-ar,es;q=0.8,en-us;q=0.5,en;q=0.3     Content-Length  723     Content-Type    application/xml; charset=UTF-8     Host        localhost     Referer     http://localhost/visor/index.html     User-Agent      Mozilla/5.0 (X11; Ubuntu; Linux i686; rv:24.0) Gecko/20100101 Firefox/24.0     X-Requested-With    XMLHttpRequest      Any clue?? I've tried everything I know and googled it without many results. Thanks in advance, and sorry for my grammatical mistakes. Greets from Argentina. Cristian


title: How can I use wms service layers to edit in openlayers?

text: I'm new to openlayers. I'd like to edit a feature layer of wms service in openlayers. I created a wms service using ArcGIS server 10. I couldn't use the wms layers to edit using the openlayers editor. Is there any solution to set layers of the wms service (ex: layer 1 of wms service) to edit in openlayers? In other words, how can I set the wms layer to `OpenLayers.Control.EditingToolbar()`?               <!DOCTYPE html>                   <title>Creating a simple map</title>         <meta http-equiv="Content-Type"               content="text/html; charset=UTF-8">


title: How to display direction of track with arrows using OpenLayers

text: I would like to display track on the map using OpenLayers in the way like this: http://jsfiddle.net/ _DR_ /hgCUq/1/light/ How to show direction of the track with arrows here? Any other approach would be great as well. I googled https://github.com/Sirtrack/PointTrackWithArrows but this project looks buggy and conflicts with OpenLayers.Layer.PointTrack


title: What statistic / process should I be using?

text: I am trying to determine the best process to analyse some predicted water levels against observed values at specific points. I can do this simply as 1 to 1 relationship to get a r square. I would however like to analyse this relationship with a third variable such as the distance from the closest surface water (ie lake). Is there a ARCGIS 10.1 statistics tool or process that anyone can recommend? Thank you for any suggestions.


title: Draw an item from coordinates in QGIS?

text: Is it possible to draw a line or to draw a triangle or another figure by digit the coordinates of the vertex? I use QGIS 1.8 on windows 7 system.


title: How identify Layer name when removing it from legendInterface qgis?

text: I am trying to connect the removing layer action with a slot for doing some processing after that. I already try this:               QObject.connect(self.iface.legendInterface(), SIGNAL("itemRemoved()"),self.RemoveScenario)      but i need to identify the removed layer name. anyone have an idea how to do that?


title: Enabling driving distance in Pgrouting

text: I would like to enable pgr_drivingDistance function fo Pgrouting. I followed this manual: http://pgrouting.org/docs/1.x/install.html, but I can't execute the last step:               psql -U postgres -f /usr/share/postlbs/routing_dd.sql routing     psql -U postgres -f /usr/share/postlbs/routing_dd_wrappers.sql routing      These files don't exist in that adresses. And when I try to change the adress from "/usr/share/postlbs/routing_dd.sql" to the path to routing_dd.sql in pgrouting, I get these errors:               psql:/home/user/pgrouting/src/driving_distance/sql/routing_dd.sql:27: ERROR:  permission denied for language c     psql:/home/user/pgrouting/src/driving_distance/sql/routing_dd.sql:37: ERROR:  permission denied for language c     psql:/home/user/pgrouting/src/driving_distance/sql/routing_dd.sql:80: ERROR:  must be owner of function pgr_pointsaspolygon      Until this step, everything goes well (and I really don't forget to compile Pgrouting with DrivingDistance functionality). Does anyone have an idea what I'm doing wrong? EDIT: Pgrouting with driving distance should be installed and enabled now, but when using drivingDistance function, I'm getting this error:                File "/usr/local/lib/python2.7/dist-packages/django/core/handlers/base.py", line 114, in get_response         response = wrapped_callback(request, *callback_args, **callback_kwargs)         cursor2.execute("SELECT id1::int4,cost::float8 from pgr_drivingDistance(%s, %s::int4, 10.0,false,false)", [query, source])       File "/usr/local/lib/python2.7/dist-packages/django/db/backends/util.py", line 69, in execute         return super(CursorDebugWrapper, self).execute(sql, params)       File "/usr/local/lib/python2.7/dist-packages/django/db/backends/util.py", line 53, in execute         return self.cursor.execute(sql, params)       File "/usr/local/lib/python2.7/dist-packages/django/db/utils.py", line 99, in __exit__         six.reraise(dj_exc_type, dj_exc_value, traceback)       File "/usr/local/lib/python2.7/dist-packages/django/db/backends/util.py", line 53, in execute         return self.cursor.execute(sql, params)     DatabaseError: error with no message from the libpq     [06/Apr/2014 09:05:46] "POST /map/findroute/ HTTP/1.1" 500 16972     Traceback (most recent call last):       File "/usr/lib/python2.7/wsgiref/handlers.py", line 86, in run         self.finish_response()       File "/usr/lib/python2.7/wsgiref/handlers.py", line 129, in finish_response         self.close()       File "/usr/lib/python2.7/wsgiref/simple_server.py", line 36, in close         SimpleHandler.close(self)       File "/usr/lib/python2.7/wsgiref/handlers.py", line 257, in close         self.result.close()       File "/usr/local/lib/python2.7/dist-packages/django/http/response.py", line 308, in close         signals.request_finished.send(sender=self._handler_class)       File "/usr/local/lib/python2.7/dist-packages/django/dispatch/dispatcher.py", line 185, in send         response = receiver(signal=self, sender=sender, **named)       File "/usr/local/lib/python2.7/dist-packages/django/db/__init__.py", line 94, in close_old_connections         conn.close_if_unusable_or_obsolete()       File "/usr/local/lib/python2.7/dist-packages/django/db/backends/__init__.py", line 462, in close_if_unusable_or_obsolete         if self.is_usable():       File "/usr/local/lib/python2.7/dist-packages/django/db/backends/postgresql_psycopg2/base.py", line 192, in is_usable         self.connection.cursor().execute("SELECT 1")     InterfaceError: connection already closed     [06/Apr/2014 09:05:46] "POST /map/findroute/ HTTP/1.1" 500 59     ----------------------------------------     Exception happened during processing of request from ('127.0.0.1', 53783)     Traceback (most recent call last):       File "/usr/lib/python2.7/SocketServer.py", line 582, in process_request_thread         self.finish_request(request, client_address)       File "/usr/lib/python2.7/SocketServer.py", line 323, in finish_request         self.RequestHandlerClass(request, client_address, self)       File "/usr/local/lib/python2.7/dist-packages/django/core/servers/basehttp.py", line 126, in __init__         super(WSGIRequestHandler, self).__init__(*args, **kwargs)       File "/usr/lib/python2.7/SocketServer.py", line 638, in __init__         self.handle()       File "/usr/lib/python2.7/wsgiref/simple_server.py", line 124, in handle         handler.run(self.server.get_app())       File "/usr/lib/python2.7/wsgiref/handlers.py", line 89, in run         self.handle_error()       File "/usr/lib/python2.7/wsgiref/handlers.py", line 304, in handle_error         self.finish_response()       File "/usr/lib/python2.7/wsgiref/handlers.py", line 127, in finish_response         self.write(data)       File "/usr/lib/python2.7/wsgiref/handlers.py", line 210, in write         self.send_headers()       File "/usr/lib/python2.7/wsgiref/handlers.py", line 267, in send_headers         if not self.origin_server or self.client_is_modern():       File "/usr/lib/python2.7/wsgiref/handlers.py", line 280, in client_is_modern         return self.environ['SERVER_PROTOCOL'].upper() != 'HTTP/0.9'     TypeError: 'NoneType' object has no attribute '__getitem__'


title: Zonal Statistics Freezing QGIS 2.2

text: Running on Windows; 8gb memory, 500gb HD, i5 cpu; Raster file: 250mb I am having a problem with zonal statistics advancing past 17%. I am trying to run statistics on a crop grown in certain counties. In the image below, you can see the counties of Washington state in the vector buffer. The raster data is dispersed throughout the state. I have run zonal statistics on much larger raster files. The raster image comes from CDL NASS portal, which contains many more pixel values other than the one I have isolated. But it allows for a researcher to isolate the variable (crop) they are interested in and download that image without the rest of the data. Any ideas? ![enter image description here](http://i.stack.imgur.com/UA5Q2.png)


title: Visualizing large GML files on the web browser

text: There is a GMl file on the "C" location. When I want to visualize the gml data on the browser with using XSLT in c# .Net envoirenment, I am not able to visualize files which is not too small. My visualizing code is; Xml1.DocumentSource = filepath2; Xml1.TransformSource = Server.MapPath("~/Xslt/transform.xsl");


title: PyQGIS: QgsMapTool to add point feature

text: I am having issues programmatically adding point features using QGIS 2.0. I am trying to add a point feature, select another layer's nearby feature, and open the feature form to enter attribute data for the point. I don't get any errors, but the features are not adding even after I validate the edit. I definitely have a bug somewhere around here, because sometimes adding the feature causes QGIS to crash.                   self.fuitesLayer.startEditing()         fuites_caps = self.fuitesLayer.dataProvider().capabilities()         self.feat = QgsFeature()         pntGeom = QgsGeometry.fromPoint(point)         self.feat.setGeometry(pntGeom)         self.fuitesLayer.addFeatures([self.feat], False)              pntBuff = pntGeom.buffer((self.canvas.mapUnitsPerPixel() * 2), 0)         rect = pntBuff.boundingBox()         Nature_idx = self.resLayer.fieldNameIndex('Nature')         Diam_idx = self.resLayer.fieldNameIndex('Diametre')         self.selectList = []         self.natureList = []         self.diametreList = []         iter = self.resLayer.getFeatures()         for feature in iter:             if feature.geometry().intersects(rect):                 self.selectList.append(feature.id())                 self.natureList.append(feature.attributes()[Nature_idx])                 self.diametreList.append(feature.attributes()[Diam_idx])                                self.resLayer.setSelectedFeatures(self.selectList)              self.iface.openFeatureForm(self.fuitesLayer, self.feat, False)


title: How to overlay polygons that are in the same layer

text: I need to INTERSECT many polygons that are in the same layer. Actually each polygons from my layer represent an area of influence from one specific bird species. Since I used to use ArcMap, I found the following entrance in this blog that has helped me very much. But now I need to do the same analyse but using sqlite and QGIS (or spatialite).


title: How do I avoid this osm2pgsql duplicate key error in Windows?

text: I am attempting to load an OSM .bz2 file into a PostGIS database using osm2pgsql but am stuck at a duplicate key value in the nodes which violates a unique constraint - any suggestions? output is below: > > C:\Temp>osm2pgsql -s -C 3000 -U postgres -W -d POSTGIS -S > C:\temp\default.style us-west-latest.osm.bz2 osm2pgsql SVN version > 0.69-21289M > > Password: Using projection SRS 900913 (Spherical Mercator) Setting up table: > planet_osm_point NOTICE: table "planet_osm_point" does not exist, skipping > NOTICE: table "planet_osm_point_tmp" does not exist, skipping Setting up > table: planet_osm_line NOTICE: table "planet_osm_line" does not exist, > skipping NOTICE: table "planet_osm_line_tmp" does not exist, skipping > Setting up table: planet_osm_polygon NOTICE: table "planet_osm_polygon" does > not exist, skipping NOTICE: table "planet_osm_polygon_tmp" does not exist, > skipping Setting up table: planet_osm_roads NOTICE: table "planet_osm_roads" > does not exist, skipping NOTICE: table "planet_osm_roads_tmp" does not > exist, skipping Mid: pgsql, scale=100, cache=3000MB, maxblocks=384001*8192 > Setting up table: planet_osm_nodes *** WARNING: intarray contrib module not > installed *** The resulting database will not be usable for applying diffs. > NOTICE: table "planet_osm_nodes" does not exist, skipping Setting up table: > planet_osm_ways NOTICE: table "planet_osm_ways" does not exist, skipping > Setting up table: planet_osm_rels NOTICE: table "planet_osm_rels" does not > exist, skipping > > !! You are running this on 32bit system, so at most !! 3GB of RAM can be > used. If you encounter unexpected !! exceptions during import, you should > try running in slim !! mode using parameter -s. > > Reading in file: us-west-latest.osm.bz2 Processing: Node(113650k) Way(0k) > Relation(0k)COPY_END for COPY planet_osm_nodes FROM STDIN; failed: ERROR: > duplicate key value violates unique constraint "planet_osm_node s_pkey" > DETAIL: Key (id)=(2147483647) already exists. CONTEXT: COPY > planet_osm_nodes, line 98676608 > > Error occurred, cleaning up Thanks!


title: Collapse dual carriageway lines in PostGIS

text: I have OpenStreetMap data, which represents larger roads (e.g. dual carriageways, aka divided highways) as parallel lines. The question is answered for ArcGIS or Skeletron here, for example using ESRI's Merge Divided Roads or Collapse Dual Lines to Centerline tools. Does anyone have suggestions for how to do this within PostGIS? I am working on a server without ArcGIS, and would like to avoid the overhead of transferring the data back and forth across machines. **_Update_** The linked question isn't really answered, and is looking for a solution that is not PostGIS. This should be kept open.


title: how to aggregate the rates, represented by points, within X radius

text: the point represents rates like 10%, etc. (this was also pretty clear from 2 years ago) -- been trying to do this for 2 yrs... the points dont represent _numbers of events_ since they are rates (Finding largest number of observations within a given radius) but the points have geo ids so they are attached to location, not distance (whatever _distance_ is to imply) -- update: it represents z-value but we're just dealing with rates


title: Python Error while Publishing mxd ArcGIS 10.1

text: I was in the process of publishing my mxd , I had reached till the analysis step where I encontered the error at staging the service. Before this below mentioned code I had created a connection ags file.               analysis = arcpy.mapping.CreateMapSDDraft(MapDoc, SDDraft, serviceName, "ARCGIS_SERVER")     if analysis['errors']=={}:          arcpy.StageService_server(SDDraft, sd)          arcpy.UploadServiceDefinition_server(sd, connectionFilePath)      This is the error I encountered at the arcpy.StageService_server(SDDraft, sd) step: ERROR 001272: Analyzer errors were encountered (codes = 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3). Failed to execute (StageService). What should I do now?


title: identify sidewalks only at one side of the roadway

text: I got a task which I don't know how to save some labor work. So I have two shapefiles in hand. One is roadway shapefile and the other is sidewalk shapefile. My task is to identify where there is only one-side sidewalk. The problem is when I created the sidewalk file, I didn't record the direction. Even if I did, I still can't come up with a way to automatically identify the roadway parts with sidewalks on one side.Is there anyway to make the identification easier? The only thoughts I have right now is to go through the sidewalk file and pick out the one-side sidewalks one by one. Thanks for any of your thoughts. Here is a screenshot.The red line represents sidewalks and the blue line represents the center line of roadway. ![enter image description here](http://i.stack.imgur.com/aXLi5.jpg)


title: Adding dynamic picture elements to mapbook

text: I have a feature class full of traffic signal data including up to 3 pictures of each traffic signal. I want to create a series of data driven pages that include 3 picture/data frames that dynamically update the pictures based on HTML link in the attribute table.