title: POSTGIS: function st_collect(geography) does not exist

text: I am trying to use some Geography functions, but I am getting an error that (all) functions are not existing. I have checked other Stackexchange posts (e.g. this one), but the solutions did not work for me. Since no function is available, I assume that something is wrong with the geography tables. I am using POSTGIS 2.0.1, therefore the `geography` functions should be supported.               database# select postgis_full_version();      POSTGIS="2.0.1 r9979" GEOS="3.3.5-CAPI-1.7.5" PROJ="Rel. 4.8.0, 6 March 2012" GDAL="GDAL     1.9.1, released 2012/05/15" LIBXML="2.7.3" LIBJSON="UNKNOWN" TOPOLOGY RASTER      I can see the geography tables in the database and the database owner also owns the geography tables.                public   | geography_columns | view     | developer      public   | geometry_columns  | view     | developer      public   | raster_columns    | view     | developer      public   | raster_overviews  | view     | developer      **What could cause the error below?**               function st_collect(geography) does not exist     LINE 1: SELECT ST_xMin(ST_Collect(location)) AS id, ST_xMax(ST_Colle...                            ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.


title: Issue On Having HomeButton and Overview Together on ArcGIS Sever JavaScript API

text: Using ArcGIS 10.1 and ArcGIS Server Javascript Api I am having issue on having Home Button and Map Overview on the map. for what ever reason the home button is not showing up the map Here is the Script which I have is:               <script>       require([         "esri/map",          "esri/dijit/OverviewMap", "dojo/parser","dijit/layout/BorderContainer", "dijit/layout/ContentPane",         "esri/dijit/HomeButton", "esri/layers/FeatureLayer", "dojo/dom-construct", "dojo/domReady!"       ], function(           Map, OverviewMap,       parser, HomeButton, FeatureLayer                 ) {         parser.parse();          var map = new Map("mapDiv", {           basemap: "topo",           center: [-126.416, 55.781],           zoom: 6         });       //===================================================== Overview         var overviewMapDijit = new OverviewMap({           map: map,           visible: true         });         overviewMapDijit.startup();           //===================================================== Shapefile         //add a layer to the map         var featureLayer = new FeatureLayer("http://somewhere/1", {           mode: FeatureLayer.MODE_ONDEMAND,         });      map.addLayer(featureLayer);           //===================================================== Home Button     var home = new HomeButton({     map: map     }, "HomeButton");      home.startup();            });      </script>      and CSS for Home Button is:               #HomeButton { position: absolute;  top: 120px; left: 50px;  z-index: 2; }      Can you please let me know why this is happening?


title: How to use Button with Behavior in Silverlight Legend?

text: I'm trying to use a button with behavior in the MapLayerTemplate of an Esri Silverlight SDK 2.1 legend. When I click the button and then drag on the map, the map pans instead of showing a selection rectangle. The behavior code is pulled from this sample, while the legend is from this sample. If I put the same xaml code within the layoutroot, it works as expected - I see a selection rectangle and the selection completes successfully.               <esri:Legend x:Name="MyLegend" Map="{Binding ElementName=MyMap}"                   LayerIDs="ProjectPolys,AssetPoints,AssetLines,BaseLayers,WaterInfrastructure,SewerInfrastructure"                 LayerItemsMode="Tree"                  ShowOnlyVisibleLayers="False"                 Refreshed="Legend_Refreshed">         <esri:Legend.MapLayerTemplate>             <DataTemplate>                 <StackPanel Orientation="Horizontal">                 <CheckBox Content="{Binding Label}"                         IsChecked="{Binding IsEnabled, Mode=TwoWay}"                         IsEnabled="{Binding IsInScaleRange}">                                                     </CheckBox>                     <Button Content="Select">                         <i:Interaction.Triggers>                             <i:EventTrigger EventName="Click">                                 <esriBehaviors:SpatialQueryAction                                              DrawMode="Rectangle"                                              LayerID="MyGraphicsLayer"                                              Url="http://hq-gis-01/ArcGIS/rest/services/WaterInfrastructure2/MapServer/15"                                     Symbol="{StaticResource GraphicsLayerFillSymbol}"                                              TargetName="MyMap" />                             </i:EventTrigger>                         </i:Interaction.Triggers>                     </Button>                     <Slider Maximum="1" Value="{Binding Layer.Opacity, Mode=TwoWay}" Width="50" />                 </StackPanel>             </DataTemplate>         </esri:Legend.MapLayerTemplate>         <esri:Legend.LayerTemplate>             <DataTemplate>                 <CheckBox Content="{Binding Label}"                 IsChecked="{Binding IsEnabled, Mode=TwoWay}"                 IsEnabled="{Binding IsInScaleRange}" >                 </CheckBox>             </DataTemplate>         </esri:Legend.LayerTemplate>     </esri:Legend>      **Update:** Following Seth's suggestion I tried wiring the trigger in the code-behind. I get the same results. The code behind approach fails for the button in the datatemplate, but works for the one on a stackpanel sibling to the map.               private void Assign_Behavior_Click(object sender, RoutedEventArgs e)     {         var button = sender as Button;         ESRI.ArcGIS.Client.Actions.SpatialQueryAction sqa = new ESRI.ArcGIS.Client.Actions.SpatialQueryAction();         sqa.DrawMode = DrawMode.Rectangle;         sqa.LayerID = "MyGraphicsLayer";         sqa.Url = @"http://hq-gis-01/ArcGIS/rest/services/WaterInfrastructure2/MapServer/15";         sqa.Symbol = LayoutRoot.Resources["GraphicsLayerFillSymbol"] as Symbol;         sqa.TargetName = "MyMap";              var et = new System.Windows.Interactivity.EventTrigger("Click");         et.Actions.Add(sqa);         var triggerColl = System.Windows.Interactivity.Interaction.GetTriggers(button);         triggerColl.Add(et);     }      I'm beginning to think this is not an Esri issue - I can't find an example anywhere showing how to attach trigger actions from buttons in a datatemplate.


title: How to compare geometry in PostGiS?

text: I have a two tables.   1 table have srid=4326   2 table have srid=70066 I want to compare geometry columns of this tables. I find function `ST_Equals(geometry, geometry)` but its not what i want. I need something what can compare geometries which slightly different from each other like on picture. ![enter image description here](http://i.stack.imgur.com/qEPZh.jpg) **UPDATE**   i have test polygon in tables with area=9098.79m^2 in first table and with area=9721.82m^2 in sacond table. I try this query               SELECT   st_area(ST_Intersection(ST_TRANSFORM(filedata.the_geom,70066),filedata_temp.the_geom))     FROM filedata, filedata_temp     Where filedata.num=1      srid=70066 is rectangular coordinate system and i get query result `8910.5625` why so much?


title: Does a projection's WKT reflect the computer's locale?

text: This should be quick - we're trying to determine if there are problems parsing some WKT text because the system locale considers the decimal separator to be a comma instead of a period. Anyone know? References would be handy so I can show the devs. Our specific implementation is using Esri's Web ADF but I am curious to know how other systems handle it, and whether or not the OGC has established an official opinion.


title: esri pgdb in qgis

text: Using osgeo4w 1.5, when I go to add vector layer and select file type, navigate to my esri pgdb. I get a not a valid or recognized data source. If I use the database type I get a page that defaults back to type odbc even though I just chose esri personal geodatabase in the previous dialog, it also requires host, database, port. There seems to be no way to navigate to a file. The help file is also blank and says new ogr database connection dialog. Do I not have all the ogr support installed? Has anyone looked into adding support for FME in qgis? I use it to connect to many formats in all my other software. Thanks


title: Vector boundary generation -similar to Image Boundary

text: I want to automate the creation of boundary polygons from vector files. This is not footprint or extent creation (I already have scripts to do this). I want to create a tool similar to the QGIS, Image Boundary plugin which creates a single polygon which describes only the areas which have data. It should work on points/lines and polygons. For the points/lines, I guess it would need to create a buffer and then use that in generating a multipart boundary polygon (so it won't be exact but based on the buffer you set). The purpose to use the output in a metadata catalogue where when you click on a location you only see a list of all files which have data in that spot and not just that the extent/footprint is within. I found a tool from a post here that claims to do this but it is circa 2008 and isn't python based- http://blog.shoutis.org/2008/10/announcement-boundary- generator-v-01.html An opensource QGIS/GDAL solution would be ideal but arcpy is fine as well.


title: How to deal with empty corners in raster when calculating fragmentation index

text: I have clipped buffer area around coordinate from one big raster (forest cover). Such created raster was loaded in `R` where `PatchStat` function (library `SDMTools`) gives **fragmentation index**. After detailed look of such created raster I noticed that I have **equal classification for non-forest area (0) and corners (also 0)**. > My question is whether such surplus of empty area make **an influence on the > value of fragmentation index**. _I think solution bellow may be helpful, but > I do not know how to do it._ ![enter image description here](http://i.stack.imgur.com/uPHp8.jpg)


title: GRASS module g.mremove returns irrational outputs

text: I want to use **g.mremove** to remove multiple files in my mapset,               import os     import sys     sys.path.append(os.path.join(os.environ['GISBASE'], 'etc', 'python'))     import grass.script as g     import grass.script.setup as gsetup     gisbase = os.environ['GISBASE']     gisdb = 'C:\Users\Heinz\Documents\grassdata'     location = 'newLocation'     mapset = 'TC'     gsetup.init(gisbase, gisdb, location, mapset)          g.run_command('g.mremove', rast = 'b_*')     print g.read_command('g.list', _type='rast')      but the output still has b-beginning files,               *** Remote Interpreter Reinitialized  ***     >>>      ----------------------------------------------     raster files available in mapset <TC>:     acc        b16        b23        b30        b7         bas3       basin_test     b10        b17        b24        b31        b8         bas4       dra     b11        b18        b25        b32        b9         bas5       hbas     b12        b19        b26        b33        bas        bas6       hc     b13        b20        b27        b4         bas1       bas7       str     b14        b21        b28        b5         bas10      bas8       vis     b15        b22        b29        b6         bas2       bas9          raster files available in mapset <PERMANENT>:     hc          ----------------------------------------------      How to fix this problem？ I am using GRASS 6.4.3, python 2.7.6 under Windows 8.1 x64. * * * **Update#1** I have edited the g.run_command line above to this,               g.run_command('g.mremove', flags = 'f', rast = 'b_*')      the output is the same, then I try to do it in the python shell of the GRASS layer manager, and here are the code and output,               >>>g.run_command('g.mremove', flags = 'f', rast = 'b_*')     0     >>>g.read_command('g.list', _type = 'rast')     ----------------------------------------------     raster files available in mapset <TC>:     acc        b16        b23        b30        b8         bas5       hc     b1         b17        b24        b31        b9         bas6       str     b10        b18        b25        b32        bas        bas7       vis     b11        b19        b26        b33        bas1       bas8     b12        b2         b27        b4         bas10      bas9     b13        b20        b28        b5         bas2       basin_test     b14        b21        b29        b6         bas3       dra     b15        b22        b3         b7         bas4       hbas          raster files available in mapset <PERMANENT>:     hc      Still give me irrational outputs.


title: ogr2ogr MIF to KML fails on writing kml in linux

text: I have a bunch of MIF files which I would like to translate to KML format. ogr2ogr -f "KML" geo-mif geo-kml FAILURE: Unable to open datasource `geo-kml' with the following drivers. seems that ogr2ogr has issues writing KML ...? do i need to install some additional stuff? other formats like Mapinfo file works like a charm e.g.: ogr2ogr -f "Mapinfo file" geo-kml geo-mif works perfectly thx for any direections... hugo


title: Problem with finding point within the box in MongoDB

text: I'm doing a search within box. It works perfect except only one situation. Let me describe it. Suppose the point has coordinates: (lng:-155, lat:19) If we have box: sw - (-170, 10), ne (29,69) that everything is fine and we can find this point because it's inside of this box. But if we move the box a little bit and now it's coordinates are sw - (144, 10), ne (-75,69) that results are zero although the point is still inside of this box. Please see the picture describing. Maybe I misunderstand how the search "within" box but in this case could you please tell me how to do it properly to find the point in the second situation. Thanks!! ![enter image description here](http://i.stack.imgur.com/wLQeW.png)


title: Guidance when to style your features in Tilemill vs. Leaflet?

text: I've found this generally to be dependent on each project's goals (or what the client wants), but when you have the options, what any rules of thumb or guidelines do use you in a project to determine whether to:   * use/style a base map in Tilemill/tilemill2; create a separate layer within TM of your unique data that you're intending to overlay the basemap; style the data within TM2.   * use a base map in Tilemill/tilemill2, load it in leaflet, add the features from a geojson file in leaflet with the features' styles in the geojson file itself.    * use a base map in Tilemill/tilemill2, load it in leaflet, add the features from a geojson file in leaflet and then style the features in the javascript on the page.  In cases where you're styling a lot of features (more than a few hundred?), I would add all of it within Tilemill since you're not dealing with the client side rendering.. In cases where your features are than 10 or so, I'd style them within the html file, I'd assume. What other factors should I keep in mind when I'm weighing these approaches? (I'm using TM and tm2 interchangeably here; I'm aware they're separate programs)


title: Grayscale Tiff Images

text: I have a gray scale tiff image. I can extract the values from this tiff using numpy and make a new image with blue scale. But how can I convert it to shades of blue instead of shades of gray directly from the Tiff? Is it even possible? Edit: Well I have two sets of tiff images(data), one with color-palette and another with gray scale ,I could generate tiles from color-palette one using generate_tiles.py, and show them on a map, but with gray scale one I can't do that, so I am wondering if I could change this gray scale tiff to color- palette one (with custom colours) so I can generate tiles and visualize the data.


title: Changing raster cells values inside of a polygon to the closest cells just outside the polygon!

text: I have been wrestling with this for a while and couldn't find out how I can do this. I want each cell of a raster inside of a polygon boundary to have the value of the closest cell just outside of the polygon. I think the closest answer was the one at this question, but I couldn't apply it! My plan was:   1. Use "Polygon to Raster" to change the polygon to a raster. The cells inside the polygon become 1 and the points outside of the polygon are nodata.   2. Use Con(IsNull(Raster),0,NULL) and assign it to Raster2   3. Add Raster1 and Raster2 to get nodata inside the polygon   4. It should be easy from here with nodata inside the polygon. We just use "Euclidean Allocation" and it does the job. I'm stuck at step 2 and get an error: > [SELECT * FROM VAT_Extent_RAS WHERE IsNull( "Value" )] Failed to execute > (Con). Do you have any thoughts? Is there an easier way to do this? **Update:** I'll try to summarize my progress here. I want to change all raster cells inside the polygon to "nodata" while the raster cells outside of the polygon remain the same. I used polygon to raster function first. This gives me a raster with 0 cell inside the polygon and "nodata" outside of the polygon. What do I do next?


title: Find maximum and minimum value in a list of raster

text: I trying to write a standardisation script that standardises each raster dataset in a list of rasters by the highest and lowest values found in the list of rasters. The standardisation code is as follows:               for raster in rasters:       output = (Raster(raster) - min_value) *100 / (max_value - min_value) + 0      Where the max_value and min_value are the maximum and minimum values in the raster list. I initially though I would be able to use the GetRasterProperties_management function on the raster list but this doesn't work because the raster list is just a list! The other idea I had was to loop through each raster in the list getting the max and min values and then outputting them as a table of values. Then finding the max and min value in each table to use in the standardisation equation. But I'm struggling to figure out how to output max and min values as a table. So I would like help if anyone has got a better idea for the whole procedure or knows how to output max and min values from a looped GetRasterProperties_management into a table. I'm using ArcGIS 10. Thank you, Jack


title: Keep synchronization when simplify features from different data source

text: We have two shapfile data source: 1 `water-riverbank.shp` With geometry type of `Polygon` which stands for the area of water 2`water-riverline.shp`. With geometry type of `LineString` which stands for the direction of the water. We use these features to lable the water. ![enter image description here](http://i.stack.imgur.com/Jb1FK.png) Now we use mapshaper to simplify these data,however we may get this kind of result: ![enter image description here](http://i.stack.imgur.com/wqOAW.png) The simplified riverbank and riverline are not synchronized. Any way to avoid this?


title: How to get column names of layer (WMS) using OpenLayers?

text: Is it somehow possible to get layer's column name using OpenLayers? any example and documentation would be really useful? Edit: I have already tried WMS_describeLayer `topp:states` layer:               <?xml version="1.0" encoding="UTF-8"?>     <!DOCTYPE WMS_DescribeLayerResponse       SYSTEM 'http://localhost:8080/geoserver/schemas/wms/1.1.1/WMS_DescribeLayerResponse.dtd'>     <WMS_DescribeLayerResponse version="1.1.1">         <LayerDescription           name="topp:states"           owsType="WFS"           owsURL="http://localhost:8080/geoserver/wfs/WfsDispatcher?"            wfs="http://localhost:8080/geoserver/wfs/WfsDispatcher?">             <Query typeName="topp:states"/>         </LayerDescription>     </WMS_DescribeLayerResponse>


title: Displaying Interpolated surface

text: > **Possible Duplicate:**   >  Displaying/Accessing interpolated data I successfully carried out IDW interpolation procedure on rainfall data, using the spatial extent, same pixel size, and same pixel registration (in the Environment Settings window) as the NDVI image I want to process the rainfall data with. But it didn't dispaly in ArcMap even after I zoomed to layer and layer resolution. So I couldn't export the interpolated data for further image processing. What is the problem and how can I solve it? Could it be that the points in my input points feature are too scanty? Thank you very much as I eagerly await a solution.


title: Show feature classes as a checklist in a workspace in ModelBuilder

text: In my model, I would like to show all the feature classes inside a file gdb as a "checklist", similar to how the tool dialog window for `Make Query Table` looks: ![Make Query Table dialog window](http://i.stack.imgur.com/0ylg5.png) So instead of the input being a table, and the field names are populated below so the user can select which fields are needed, I would like the input to be a file gdb and the feature classes inside it are populated in the window, so only the required feature classes can be selected. Is this possible inside modelbuilder via a tool or some Python code (maybe `ListFeatureClasses` but it must display the output)? It would also be helpful to know what that type of window is called. I know the top section is the `Multiple Value` window, so what would the "checklist" window be called?


title: Problems installing QGIS 2.2 on Ubuntu 14.04

text: After major upgrade to 14.04 I tried to purge and reinstall QGIS. Unfortunately I get and error (full listing below) and cannot finish clean installation. Information in the log points errors to appear with pyspatialite and/or python-pyspatialite but I cannot find any solution to this issue. Also - QGIS actually manages to start after installation but points to another error with PyQGIS this time: ![enter image description here](http://i.stack.imgur.com/jOlUn.png) Any ideas how to solve this issue?               rdk@hal:~$ sudo apt-get install qgis python-qgis qgis-plugin-grass     Reading package lists... Done     Building dependency tree            Reading state information... Done     The following extra packages will be installed:       grass grass-core grass-doc grass-gui libmotif-common libmrm4       libopenscenegraph99 libosgearth3 libosgearthannotation3 libosgearthfeatures3       libosgearthqt3 libosgearthsymbology3 libosgearthutil3 libqgis-analysis2.2.0       libqgis-core2.2.0 libqgis-gui2.2.0 libqgis-networkanalysis2.2.0       libqgisgrass2.2.0 libqgispython2.2.0 libspatialindex3 libtcl8.5 libtk8.5       libuil4 libxm4 osgearth-data python-opengl python-pyspatialite       python-qgis-common qgis-common qgis-plugin-globe qgis-plugin-globe-common       qgis-plugin-grass-common qgis-providers qgis-providers-common tcl8.5 tk8.5     Suggested packages:       grass-dev e00compr avce00 gnuplot gpstrans python-rpy2 xml2 grass-dev-doc       openscenegraph libgle3 tcl-tclreadline     The following NEW packages will be installed       grass grass-core grass-doc grass-gui libmotif-common libmrm4       libopenscenegraph99 libosgearth3 libosgearthannotation3 libosgearthfeatures3       libosgearthqt3 libosgearthsymbology3 libosgearthutil3 libqgis-analysis2.2.0       libqgis-core2.2.0 libqgis-gui2.2.0 libqgis-networkanalysis2.2.0       libqgisgrass2.2.0 libqgispython2.2.0 libspatialindex3 libtcl8.5 libtk8.5       libuil4 libxm4 osgearth-data python-opengl python-pyspatialite python-qgis       python-qgis-common qgis qgis-common qgis-plugin-globe       qgis-plugin-globe-common qgis-plugin-grass qgis-plugin-grass-common       qgis-providers qgis-providers-common tcl8.5 tk8.5     0 to upgrade, 39 to newly install, 0 to remove and 0 not to upgrade.     Need to get 62.2 MB of archives.     After this operation, 235 MB of additional disk space will be used.     Do you want to continue? [Y/n] Y     Get:1 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main qgis-providers-common all 2.2.0-2~trusty1 [1,539 kB]     Get:2 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe libmotif-common all 2.3.4-5 [11.8 kB]     Get:3 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libqgis-core2.2.0 amd64 2.2.0-2~trusty1 [2,607 kB]     Get:4 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libqgis-gui2.2.0 amd64 2.2.0-2~trusty1 [1,913 kB]     Get:5 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main qgis-providers amd64 2.2.0-2~trusty1 [1,770 kB]     Get:6 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libosgearth3 amd64 2.5.0+dfsg-3~trusty1 [775 kB]     Get:7 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe libxm4 amd64 2.3.4-5 [967 kB]     Get:8 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libosgearthsymbology3 amd64 2.5.0+dfsg-3~trusty1 [244 kB]     Get:9 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libosgearthfeatures3 amd64 2.5.0+dfsg-3~trusty1 [246 kB]     Get:10 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libosgearthannotation3 amd64 2.5.0+dfsg-3~trusty1 [172 kB]     Get:11 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libosgearthutil3 amd64 2.5.0+dfsg-3~trusty1 [633 kB]     Get:12 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libosgearthqt3 amd64 2.5.0+dfsg-3~trusty1 [176 kB]     Get:13 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libqgis-analysis2.2.0 amd64 2.2.0-2~trusty1 [930 kB]     Get:14 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libqgis-networkanalysis2.2.0 amd64 2.2.0-2~trusty1 [794 kB]     Get:15 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libqgisgrass2.2.0 amd64 2.2.0-2~trusty1 [844 kB]     Get:16 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main libqgispython2.2.0 amd64 2.2.0-2~trusty1 [789 kB]     Get:17 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main osgearth-data all 2.5.0+dfsg-3~trusty1 [3,946 kB]     Get:18 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe libmrm4 amd64 2.3.4-5 [56.2 kB]     Get:19 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main python-qgis-common all 2.2.0-2~trusty1 [2,882 kB]     Get:20 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe libspatialindex3 amd64 1.8.1-2 [218 kB]     Get:21 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/main libtcl8.5 amd64 8.5.15-2ubuntu1 [684 kB]     Get:22 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main python-qgis amd64 2.2.0-2~trusty1 [2,671 kB]     Get:23 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/main libtk8.5 amd64 8.5.15-2ubuntu3 [641 kB]     Get:24 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main qgis-common all 2.2.0-2~trusty1 [6,335 kB]     Get:25 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe libuil4 amd64 2.3.4-5 [126 kB]     Get:26 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe grass-core amd64 6.4.3-3 [7,199 kB]     Get:27 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main qgis amd64 2.2.0-2~trusty1 [4,890 kB]     Get:28 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main qgis-plugin-globe-common all 2.2.0-2~trusty1 [780 kB]     Get:29 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main qgis-plugin-globe amd64 2.2.0-2~trusty1 [856 kB]     Get:30 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main qgis-plugin-grass-common all 2.2.0-2~trusty1 [1,198 kB]     Get:31 http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/ trusty/main qgis-plugin-grass amd64 2.2.0-2~trusty1 [1,508 kB]     Get:32 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe grass-doc all 6.4.3-3 [6,303 kB]     Get:33 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/main tcl8.5 amd64 8.5.15-2ubuntu1 [13.9 kB]     Get:34 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/main tk8.5 amd64 8.5.15-2ubuntu3 [12.0 kB]     Get:35 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe python-opengl all 3.0.2-1 [364 kB]     Get:36 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe grass-gui amd64 6.4.3-3 [1,818 kB]     Get:37 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe libopenscenegraph99 amd64 3.2.0~rc1-4 [5,283 kB]     Get:38 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe python-pyspatialite amd64 3.0.1-4 [29.4 kB]     Get:39 http://mirror.switch.ch/ftp/mirror/ubuntu/ trusty/universe grass all 6.4.3-3 [6,148 B]     Fetched 62.2 MB in 5s (10.8 MB/s)     Extract templates from packages: 100%     Selecting previously unselected package libmotif-common.     dpkg: warning: files list file for package 'gettext' missing; assuming package has no files currently installed     dpkg: warning: files list file for package 'texinfo' missing; assuming package has no files currently installed     (Reading database ... 500015 files and directories currently installed.)     Preparing to unpack .../libmotif-common_2.3.4-5_all.deb ...     Unpacking libmotif-common (2.3.4-5) ...     Selecting previously unselected package libxm4:amd64.     Preparing to unpack .../libxm4_2.3.4-5_amd64.deb ...     Unpacking libxm4:amd64 (2.3.4-5) ...     Selecting previously unselected package libmrm4:amd64.     Preparing to unpack .../libmrm4_2.3.4-5_amd64.deb ...     Unpacking libmrm4:amd64 (2.3.4-5) ...     Selecting previously unselected package libspatialindex3:amd64.     Preparing to unpack .../libspatialindex3_1.8.1-2_amd64.deb ...     Unpacking libspatialindex3:amd64 (1.8.1-2) ...     Selecting previously unselected package libtcl8.5:amd64.     Preparing to unpack .../libtcl8.5_8.5.15-2ubuntu1_amd64.deb ...     Unpacking libtcl8.5:amd64 (8.5.15-2ubuntu1) ...     Selecting previously unselected package libtk8.5:amd64.     Preparing to unpack .../libtk8.5_8.5.15-2ubuntu3_amd64.deb ...     Unpacking libtk8.5:amd64 (8.5.15-2ubuntu3) ...     Selecting previously unselected package libuil4:amd64.     Preparing to unpack .../libuil4_2.3.4-5_amd64.deb ...     Unpacking libuil4:amd64 (2.3.4-5) ...     Selecting previously unselected package qgis-providers-common.     Preparing to unpack .../qgis-providers-common_2.2.0-2~trusty1_all.deb ...     Unpacking qgis-providers-common (2.2.0-2~trusty1) ...     Selecting previously unselected package libqgis-core2.2.0.     Preparing to unpack .../libqgis-core2.2.0_2.2.0-2~trusty1_amd64.deb ...     Unpacking libqgis-core2.2.0 (2.2.0-2~trusty1) ...     Selecting previously unselected package libqgis-gui2.2.0.     Preparing to unpack .../libqgis-gui2.2.0_2.2.0-2~trusty1_amd64.deb ...     Unpacking libqgis-gui2.2.0 (2.2.0-2~trusty1) ...     Selecting previously unselected package qgis-providers.     Preparing to unpack .../qgis-providers_2.2.0-2~trusty1_amd64.deb ...     Unpacking qgis-providers (2.2.0-2~trusty1) ...     Selecting previously unselected package grass-core.     Preparing to unpack .../grass-core_6.4.3-3_amd64.deb ...     Unpacking grass-core (6.4.3-3) ...     Selecting previously unselected package grass-doc.     Preparing to unpack .../grass-doc_6.4.3-3_all.deb ...     Unpacking grass-doc (6.4.3-3) ...     Selecting previously unselected package tcl8.5.     Preparing to unpack .../tcl8.5_8.5.15-2ubuntu1_amd64.deb ...     Unpacking tcl8.5 (8.5.15-2ubuntu1) ...     Selecting previously unselected package tk8.5.     Preparing to unpack .../tk8.5_8.5.15-2ubuntu3_amd64.deb ...     Unpacking tk8.5 (8.5.15-2ubuntu3) ...     Selecting previously unselected package python-opengl.     Preparing to unpack .../python-opengl_3.0.2-1_all.deb ...     Unpacking python-opengl (3.0.2-1) ...     Selecting previously unselected package grass-gui.     Preparing to unpack .../grass-gui_6.4.3-3_amd64.deb ...     Unpacking grass-gui (6.4.3-3) ...     Selecting previously unselected package libopenscenegraph99.     Preparing to unpack .../libopenscenegraph99_3.2.0~rc1-4_amd64.deb ...     Unpacking libopenscenegraph99 (3.2.0~rc1-4) ...     Selecting previously unselected package libosgearth3.     Preparing to unpack .../libosgearth3_2.5.0+dfsg-3~trusty1_amd64.deb ...     Unpacking libosgearth3 (2.5.0+dfsg-3~trusty1) ...     Selecting previously unselected package libosgearthsymbology3.     Preparing to unpack .../libosgearthsymbology3_2.5.0+dfsg-3~trusty1_amd64.deb ...     Unpacking libosgearthsymbology3 (2.5.0+dfsg-3~trusty1) ...     Selecting previously unselected package libosgearthfeatures3.     Preparing to unpack .../libosgearthfeatures3_2.5.0+dfsg-3~trusty1_amd64.deb ...     Unpacking libosgearthfeatures3 (2.5.0+dfsg-3~trusty1) ...     Selecting previously unselected package libosgearthannotation3.     Preparing to unpack .../libosgearthannotation3_2.5.0+dfsg-3~trusty1_amd64.deb ...     Unpacking libosgearthannotation3 (2.5.0+dfsg-3~trusty1) ...     Selecting previously unselected package libosgearthutil3.     Preparing to unpack .../libosgearthutil3_2.5.0+dfsg-3~trusty1_amd64.deb ...     Unpacking libosgearthutil3 (2.5.0+dfsg-3~trusty1) ...     Selecting previously unselected package libosgearthqt3.     Preparing to unpack .../libosgearthqt3_2.5.0+dfsg-3~trusty1_amd64.deb ...     Unpacking libosgearthqt3 (2.5.0+dfsg-3~trusty1) ...     Selecting previously unselected package libqgis-analysis2.2.0.     Preparing to unpack .../libqgis-analysis2.2.0_2.2.0-2~trusty1_amd64.deb ...     Unpacking libqgis-analysis2.2.0 (2.2.0-2~trusty1) ...     Selecting previously unselected package libqgis-networkanalysis2.2.0.     Preparing to unpack .../libqgis-networkanalysis2.2.0_2.2.0-2~trusty1_amd64.deb ...     Unpacking libqgis-networkanalysis2.2.0 (2.2.0-2~trusty1) ...     Selecting previously unselected package libqgisgrass2.2.0.     Preparing to unpack .../libqgisgrass2.2.0_2.2.0-2~trusty1_amd64.deb ...     Unpacking libqgisgrass2.2.0 (2.2.0-2~trusty1) ...     Selecting previously unselected package libqgispython2.2.0.     Preparing to unpack .../libqgispython2.2.0_2.2.0-2~trusty1_amd64.deb ...     Unpacking libqgispython2.2.0 (2.2.0-2~trusty1) ...     Selecting previously unselected package osgearth-data.     Preparing to unpack .../osgearth-data_2.5.0+dfsg-3~trusty1_all.deb ...     Unpacking osgearth-data (2.5.0+dfsg-3~trusty1) ...     Selecting previously unselected package python-qgis-common.     Preparing to unpack .../python-qgis-common_2.2.0-2~trusty1_all.deb ...     Unpacking python-qgis-common (2.2.0-2~trusty1) ...     Selecting previously unselected package python-pyspatialite.     Preparing to unpack .../python-pyspatialite_3.0.1-4_amd64.deb ...     Unpacking python-pyspatialite (3.0.1-4) ...     Preparing to unpack .../python-qgis_2.2.0-2~trusty1_amd64.deb ...     Unpacking python-qgis (2.2.0-2~trusty1) ...     dpkg: error processing archive /var/cache/apt/archives/python-qgis_2.2.0-2~trusty1_amd64.deb (--unpack):      trying to overwrite '/usr/lib/python2.7/dist-packages/pyspatialite/dump.py', which is also in package python-pyspatialite 3.0.1-4     dpkg-deb: error: subprocess paste was killed by signal (Broken pipe)     Selecting previously unselected package qgis-common.     Preparing to unpack .../qgis-common_2.2.0-2~trusty1_all.deb ...     Unpacking qgis-common (2.2.0-2~trusty1) ...     Selecting previously unselected package qgis.     Preparing to unpack .../qgis_2.2.0-2~trusty1_amd64.deb ...     Unpacking qgis (2.2.0-2~trusty1) ...     Selecting previously unselected package qgis-plugin-globe-common.     Preparing to unpack .../qgis-plugin-globe-common_2.2.0-2~trusty1_all.deb ...     Unpacking qgis-plugin-globe-common (2.2.0-2~trusty1) ...     Selecting previously unselected package qgis-plugin-globe.     Preparing to unpack .../qgis-plugin-globe_2.2.0-2~trusty1_amd64.deb ...     Unpacking qgis-plugin-globe (2.2.0-2~trusty1) ...     Selecting previously unselected package qgis-plugin-grass-common.     Preparing to unpack .../qgis-plugin-grass-common_2.2.0-2~trusty1_all.deb ...     Unpacking qgis-plugin-grass-common (2.2.0-2~trusty1) ...     Selecting previously unselected package grass.     Preparing to unpack .../archives/grass_6.4.3-3_all.deb ...     Unpacking grass (6.4.3-3) ...     Selecting previously unselected package qgis-plugin-grass.     Preparing to unpack .../qgis-plugin-grass_2.2.0-2~trusty1_amd64.deb ...     Adding 'diversion of /usr/bin/qgis to /usr/bin/qgis.bin by qgis-plugin-grass'     Adding 'diversion of /usr/bin/qbrowser to /usr/bin/qbrowser.bin by qgis-plugin-grass'     Unpacking qgis-plugin-grass (2.2.0-2~trusty1) ...     Processing triggers for man-db (2.6.7.1-1) ...     Processing triggers for bamfdaemon (0.5.1+14.04.20140409-0ubuntu1) ...     Rebuilding /usr/share/applications/bamf-2.index...     Processing triggers for desktop-file-utils (0.22-1ubuntu1) ...     Processing triggers for gnome-menus (3.10.1-0ubuntu2) ...     Processing triggers for mime-support (3.54ubuntu1) ...     Processing triggers for menu (2.1.46ubuntu1) ...     Processing triggers for hicolor-icon-theme (0.13-1) ...     Processing triggers for doc-base (0.10.5) ...     Processing 1 added doc-base file...     Registering documents with scrollkeeper...     Processing triggers for shared-mime-info (1.2-0ubuntu3) ...     Errors were encountered while processing:      /var/cache/apt/archives/python-qgis_2.2.0-2~trusty1_amd64.deb     E: Sub-process /usr/bin/dpkg returned an error code (1)


title: Is there a method in ArcObjects that is equivalent to ArcPy's: SelectLayerByLocation_management?

text: I have 2 `IFeatureClass` and I want to do an INTERSECT between the two. What is the best way to do this from ArcObjects?


title: Batch Geocode Postcodes/Addresses in UTM

text: My work has a database of previous jobs (within housing plots) which are meant to have been georefernced with UTM coordinates, but this has not been completed for a number of years. New jobs can use data from old ones, and there is a facility in the database to check for properties within a certain radius of the new job using UTM coordinates. I have batch georeferenced addresses using Lat Long before using sites like GPS Visualizer, but cannot find a means to georeference in UTM, and converting lat long to UTM seems difficult. Does anyone know of a site which will help me do this, or know of a site which will convert directly to UTM coordinates? The UTM Zone is 17N (Winnipeg).


title: How to save a vector layer from QGIS to Sql Server 2008 or 2012?

text: using QGIS, i've loaded up some administrative boundary data of New Zealand. I wish to now import this vector data from QGIS into my Sql Server 2012. Is there any way to do this with QGIS 1.8.0-Lisboa?


title: Unable to deploy add-in updates (ArcGIS Explorer build 2500)

text: I'm having trouble deploying updates for an add-in developed for AGX Desktop (build 2500). The initial deployment is fine. If I copy my addin.eaz file to the ‘C:\Users[user]\AppData\Roaming\esri\arcgis explorer\Addins' folder, then start AGX, it works as expected. Replacing the original add-in with an updated version doesn't seem to work as the described in the **documentation** : > When a new add-in file is found in this location by ArcGIS Explorer, it will > unpack the contents of the file into a cache. ArcGIS Explorer will check if > any add-ins have been updated in this location on startup, and if necessary > unpack the updated files so that the latest version is used. If I just overwrite the original .eaz file with an updated version, then start AGX, the new version is ignored. Confusingly, if I rebuild the project in Visual Studio then start AGX, it works fine. The build information in Visual Studio shows that the new .eaz file was copied to the AddIns folder mentioned above. Why this works after a build and not when I manually copy the .eaz file is a mystery (on a different machine). The only way I can get this to work is to manually remove the add-in within AGX, then start and stop AGX, then copy the updated .eaz into the AddIns folder, then start AGX again. Does anyone know of any restrictions when deploying updates, or have I missed something? Thanks, Glenn


title: How to choose suitable pyramid levels and tile size?

text: I have a big tiff file and i want to create tiles with gdal_retile.   However I am still concerned about the pyramids levels and the tiles' size.   how can I choose suitable pyramids levels (for example 5, 6 or 18 zoom level) and the pixel size for each of my tiles?


title: python multiprocessing error in QGis (Windows)

text: I have this error when I try to use Manager() from python multiprocessing library. (Just start QGIS open Python console and execute the two following lines)               >>> import multiprocessing     >>> mgr = multiprocessing.Manager()     Traceback (most recent call last):       File "<input>", line 1, in <module>       File "C:\DEV\osgeo4w_dev\apps\Python27\lib\multiprocessing\__init__.py", line 99, in Manager         m.start()       File "C:\DEV\osgeo4w_dev\apps\Python27\lib\multiprocessing\managers.py", line 524, in start         self._process.start()       File "C:\DEV\osgeo4w_dev\apps\Python27\lib\multiprocessing\process.py", line 130, in start         self._popen = Popen(self)       File "C:\DEV\osgeo4w_dev\apps\Python27\lib\multiprocessing\forking.py", line 252, in __init__         cmd = get_command_line() + [rhandle]       File "C:\DEV\osgeo4w_dev\apps\Python27\lib\multiprocessing\forking.py", line 339, in get_command_line         if process.current_process()._identity==() and is_forking(sys.argv):     AttributeError: 'module' object has no attribute 'argv'     >>>       It happens in both 1.8 and 1.9 dev version on Windows. Of course everything is ok using QGIS on Linux!... :P It seems something related to the main Python/QGis process startup and the way the library tries to start a new process.... but I can't figure out any workaround. Any idea? (already posted this on qgis dev list but had no response...)


title: How to convert shapefile to SVG keeping internal data?

text: I am very new to GIS. I found a tutorial on how to make a choropleth using Python tools. I want to follow it but use my own data and map. However the tutorial uses a US map SVG file and I need to use a CA zipcode map. The only map that I have found in a shapefile from the US Census. I have tried using kartograph.py to convert it to SVG, but it loses the zipcode within the file when I do so. How can I convert the file to SVG and keep each path id'ed as the appropriate zipcode. I am trying to follow **this tutorial**. That should give you an idea of what I am trying to do, only I don't have a SVG formatted map. I have seen references to javascript utilities. I have never used Javascript but suppose I could use something simple.


title: Batch Fuzzifying with ModelBuilder

text: I am trying to run the Fuzzy Membership on simultaneous rasters in ModelBuilder and I am running into the problem (or bug) where it allows me to change the "Membership Type" but ends up defaulting to "MSSmall" irrespective of whether I change it in the batch table.   To put this in perspective, the "Remap" tool can be used to batch reclassify rasters, but the Fuzzy Membership tool defaults back to MSSmall even though it deceptively shows it processing with the chosen membership type. Attached are two pics showing the ModelBuilder diagram I built and the resulting interface whence double-clicked in ArcCatalog. ![Diagram of Fuzzifying Multiple Rasters](http://i.stack.imgur.com/vUllB.jpg) ![Table that allows me to change Membership Type but ends up using "MSSmall" no matter what is set.](http://i.stack.imgur.com/DuocV.jpg)


title: Import GPX to QGIS with WGS84

text: When loading a GPX-File in QGIS it is not referenced on WSG84. It´s shown with wrong coordinates. Settings are WSG84.


title: ArcGIS 10 (SP 3) - Geometric Network - which features are moved?

text: In ArcGIS 9.3 Geometric Network, you can select which features can be moved when snapping. In ArcGIS 10 (SP 3) this functionality has disappeared and you either choose to have snap tolerance or not. Does anyone know why this functionality is no longer available? When snapping 3 feature classes to create a geometric network, how do we know which feature class has moved? We only want 1 of the 3 feature classes to move. Does anyone know how we can define which feature class is to move (and snap) to the other feature classes?


title: How to remove only donut holes from polygons using ArcGIS Desktop 10.0?

text: I have a trick for removing the donut holes, which is converting the polygon to line and then creating polygon out of the lines and in the end merging all of them. The drawback of this strategy is that in removes other holes and empty areas which in reality are empty and should be kept empty. See below image: ![enter image description here](http://i.stack.imgur.com/6HDRz.jpg) Any thoughts/ideas would be appreciated.


title: geoserver wms issues

text: My WMS on a local test server was working fine with GeoExplorer but suddenly stopped working and now I don't get the LocalServer in the datasources list. I have tried reboot/restart. I tried adding it (`http://localhost:8081/geoserver/ccc_beta/wms`) manually but that doesn't work. I can however view the data in Openlayers without any issues. However Geoxplorer and other software access to WMS doesn't work. In QGIS -add WMS I get               Could not understand the response.  The wms provider said:     Could not get WMS capabilities: error occurred while parsing element at line 5086 column 11     This is probably due to an incorrect WMS Server URL.     Response was:          <?xml version="1.0" encoding="UTF-8"?>     <WMS_Capabilities version="1.3.0" updateSequence="1019" xmlns="http://www.opengis.net/wms" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.opengis.net/wms http://localhost:8081/geoserver/schemas/wms/1.3.0/capabilities_1_3_0.xsd">       <Service>         <Name>WMS</Name>         <Title>GeoServer Web Map Service</Title>      I also cleared all my stores/cache etc. Version2.2-SNAPSHOT Git Revisionadbea2ddf643dce951945d291174a4714fdd801 ==== UPDATE I put the data on to a linux install and it was all working fine till now when I restarted the machine and get the same error only in Geoexplorer (can't find local data) but everything works fine in Openlayer view. The log is               2013-07-24 09:30:58,239 INFO [geoserver.web] - Unable to find property: 'format.wms.text/html' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,239 INFO [geoserver.web] - Unable to find property: 'format.wms.application/xml' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,241 WARN [geoserver.web] - Unable to reproject layer worldto spherical mercator     2013-07-24 09:30:58,241 INFO [geoserver.web] - Unable to find property: 'format.wms.image/png; mode=8bit' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,242 INFO [geoserver.web] - Unable to find property: 'format.wms.application/vnd.google-earth.kmz' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,242 INFO [geoserver.web] - Unable to find property: 'format.wms.text/html' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,242 INFO [geoserver.web] - Unable to find property: 'format.wms.application/xml' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,244 WARN [geoserver.web] - Unable to reproject layer pop_heatmapto spherical mercator     2013-07-24 09:30:58,244 INFO [geoserver.web] - Unable to find property: 'format.wms.image/png; mode=8bit' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,244 INFO [geoserver.web] - Unable to find property: 'format.wms.application/vnd.google-earth.kmz' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,244 INFO [geoserver.web] - Unable to find property: 'format.wms.text/html' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,244 INFO [geoserver.web] - Unable to find property: 'format.wms.application/xml' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,246 WARN [geoserver.web] - Unable to reproject layer volcano_clusterto spherical mercator     2013-07-24 09:30:58,246 INFO [geoserver.web] - Unable to find property: 'format.wms.image/png; mode=8bit' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,246 INFO [geoserver.web] - Unable to find property: 'format.wms.application/vnd.google-earth.kmz' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,247 INFO [geoserver.web] - Unable to find property: 'format.wms.text/html' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:30:58,247 INFO [geoserver.web] - Unable to find property: 'format.wms.application/xml' for component:  [class=org.geoserver.web.demo.OpenGeoMapPreviewPage]     2013-07-24 09:31:26,874 INFO [geoserver.wms] -      Request: getServiceInfo     2013-07-24 09:31:26,922 INFO [geoserver.wms] -


title: Showing custom markers, shadows and popups from json data

text: The examples given in the documentation of Open Layers show very basic ways of utilising the different features. This example shows how to add shadows to the marker. This example shows how to create basic popups. Looking at the sample code, I have no idea how to combine the two together, nor display different icons for the markers. Can anyone show me a tutorial or example that combines the various aspects of custom markers, shadows and popups using json data?


title: How to install Sourcepole zonal statistics

text: I am trying to install the zonal statistics tool created/supplied by sourcepole from http://www.sourcepole.com/2010/12/16/zonal-statistics-plugin- for-qgis. The instructions area as follows: * ***To install: unpack zonal_statistics.tgz under /src/plugins, insert zonal_statistics as new subdirectory in /src/plugins/CMakeLists.txt and recompile.**** I have downloaded and unzipped the file, but I am now not sure how to add it to the QGIS software. I am using a Mac OS, and QGIS version 1.7.4. I can't find the folder /src/plugins which the instructions speak of.... and in general am very confused by the instructions (more due to my semi-computer illiteracy than the instructions themselves). Any help I would greatly appreciate. Regards Becky


title: How to draw a square from 2 coordinates?

text: I have a CSV file with the following output               id, x min, y min, x max, y max, value     1 , 0    , 0    , 100  , 100  , 5     2 , 0,   , 100  , 100  , 200  , 10     3 , 100  , 100  , 200  , 200  , 15     4 , 100  , 0    , 200  , 100  , 20      I can plot the co-ordinate points onto my basemap fine. What I'm trying to do is tell QGIS that row 1 of the CSV has the bottom left and top right corners of a square, to draw it, and to assign it a value of 5. I'm not even sure this is the right way to achieve my goal - which is to display a grid with different values attached to each square, and then produce a heat-map based on those values. Thanks.


title: GIS ANALYST + cadastral support

text: I have many points that were picked by handheld GPS and I want to fit them to their relevant cadastral plots automatically. The attribute information (plot number) for these points are exactly similar to the attribute information (plot number) for the cadastral plots. However, I am unable to automate the process of moving the many points which are outside their relevant plots. Instead I am forced to do them manually. I have tried spatial join but the points can't be moved to their plots. Any idea on how to do this in arcgis9.3?


title: gdallocationinfo from C++?

text: I need a C++ API to convert x,y Image Coordinates in their native projection to output image pixel,line coordnates. I see the gdallocationinfo documentation, but no API for it. Where in online GDAL doc I look ...or do I have to write my own?


title: Are long lines in shapefiles to be considered geodesics or straight lines in the 2D latlong space?

text: Is there a definition on how lines connect vertices in the shapefile format? In the simplest case, imagine a line with only 2 points, from 40,-118 to 40,-112 - that's somewhere random in the US - with standard WGS84 geographic coordinate system. Here is the content of the .prj file:               GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]      Is the point say 40.1,-116 north or south of the line?   1. If we consider lines are linearly interpolated in the latlong space, it follows the 40 deg parallel (small circle) and the point is north of the line.   2. If we consider lines are shortest paths on the Earth surface, it's a geodesic (great circle) with a maximum latitude at the middle of the line, higher than 40.1 deg. Then the point is south of the line.   3. Or is it simply undefined? Since shapefile format has no notion of curves, but only straight segments connecting lines. Line needs to be densified (points added along the line) to clarify this answer. If I create such scenario in QGIS, the line follows the 40 deg parallel, and would tell me answer is 1. But I wouldn't take this as a definite answer and would like to hear a more solid one.


title: GDAL - least cost path analysis height data

text: I would like to use Pythons GDAL bindings to find the least cost path through an array and then show the route on top of the original array. I am using DEM data and would like the slope to be taken into consideration for the least cost path. I have used this code GDAL - Perform Simple Least Cost Path Analysis but it returns a seperate path array and the path x,y, but height data can't be extracted. Any help? ![The DEM file and the route attached](http://i.stack.imgur.com/J6lHf.png) Secondly, any ideas of how to obtain a 3D path array layer? I've tried to create one but the resulting array has no z values


title: ArcGIS: Export GRID to TIFF keeping its (classified) symbology

text: I want to export an esri grid-raster into a *.tiff without losing the symbology I created for the grid (the symbology is classified). ![enter image description here](http://i.stack.imgur.com/VTedF.jpg) No matter how I try to solve it I get a stretched symbology in a greyscale. What I've tried so far:   * Copy/export raster via: toolbox; TOC, rightclick --> data --> export data. When checking the 'Use Renderer' and 'Force RGB' boxes on the data export menu a lot of values turn into 'noData'. Furthermore, the colors are slightly different.   * Created a layer-file from the original raster. Copied the grid and tried to set the source of the layer file to the newly copied raster. However, I could not set the new source.   * Imported the symbology using a *.lyr file (Chris W's comment) with the folder button next to the save button upper right of the image: nothing happened.   * Imported the symbology directly from my symbolized grid (Chris W's comment): nothing happened. There must be a way to solve this but I cannot find a solution... Any hints? Edit: Added clipped grid and *.lyr-file here.


title: How to show shapefiles with Bing map in asp.net

text: I want to display my shapefile data on bing map using ASP.NET. Its an web application. I don't want to use any third party tool for shapefile conversion.


title: Adding marker on a leaflet map

text: I'm very new to GIS and for learning I started with building a simple web app with GeoDjango. I am using dango-leaflet. Since I have a very limited knowledge I am facing many challenges. I'm trying to put a specifc marker on my map as a test # models.py               from django.db import models     from django.contrib.gis.geos import point     from django.contrib.gis.db import models as gismodels          class Points(gismodels.Model):           data_lat = 44.915223          data_long = -93.209741                    @property               def geom(self):              return point(self.data_lat, self.data_long, srid=4326)          def __unicode__(self):             return self.title      about.html               <body>           <h1> Map  </h1>          {% leaflet_map "main" callback="main_map_init" %}                <script type="text/javascript">                    function main_map_init (map, options){                          var dataurl = '{% url "data" %}';                         $.getJSON(dataurl, function(data) {                                L.geoJson(data).addTo(map);                    })               }           </script>          </body>      I'm having issue with data not being passed on to map. It would be nice if someone could point me to the right direction.


title: How to calculate the area between two line shapes?

text: I'm a newbie using QGIS (2.0). I'm trying to calculate the sidewalk area of a city. For that I have some shapes, extracted from Mricostation, with the topology (facades, fords, fences, etc.), but they are in line shapes. I have merged them in only two shapes: exterior line and interior line limits. I've tried to convert them to polygon, with the idea to use the "difference" tool after that, but I wasn't able to do it. I tried to use: "Geometry tools --> Lines to polygon" but it doesn't work, causing lots of holes and mistakes. "Geoprocessing tools --> Convex hull(s)" is not the way because I will need to do so many polygons and it will be an incredible waste of time. Any ideas? Thank you so much. ![enter image description here](http://i.stack.imgur.com/FxHxG.jpg)


title: OpenStreetMaps routing - starting far from a road

text: I am writing an web application, that needs routes between some markers (A, B, C, D). So far I used the this api with PHP, where I pass each pair of points (A and B, B and C, C and D .... and so on) and loop through each pair to get the coordinates and build polyline. But some of the points are not on the street, like the starting point in this example. So, technically is there a way to find if any of the points is not on the road?


title: how to change marker angle in openlayers?

text: I want show one marker with normal angle 90 degrees(it is default of openlayers) and other marker with 120 degrees in openlayers. To differentiate with the normal ones. how can i do this.


title: Cannot insert simple geometry data

text: I installed PostgreSQL 9.2.4 and PostGIS 2.0.3 on my Windows 8 machine. Then I created a table using phpPgAdmin and added a geometry column using the following statement.               SELECT AddGeometryColumn('tablename', 'area', -1, 'POLYGON', 2);      That resulted in these table columns.               +------+---------------------+--------+----------------------------------+-----------+     |Column|Type                 |Not Null|Default                           |Constraints|     +------+---------------------+--------+----------------------------------+-----------+     |id    |integer              |NOT NULL|nextval('simple_id_seq'::regclass)|PRIMARY KEY|     |name  |character varying(50)|        |                                  |           |     |area  |geometry(Polygon)    |        |                                  |           |     +------+---------------------+--------+----------------------------------+-----------+      After some simple text queries I tried to play with geometry data. This is the query I tried to execute.               INSERT INTO tablename (name, area)     VALUES (         'Lake',         GeometryFromText('POLYGON((84 104, 212 48, 321 122, 243 179, 275 238, 131 240))', -1)     );      Unfortunately, it throws an error message.               ERROR:  function geometryfromtext(unknown, integer) does not exist     LINE 4:     GeometryFromText('POLYGON((84 104, 212 48, 321 122, 243 ...                 ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.      Since `AddGeometryColumn` was successful, I assume that PostGIS is set up properly. How can I insert geometry data in my table? What am I doing wrong?


title: How to merge many shape with the same name in different folders?

text: I've many (hundreds) shape (3D points and 3D polygons) that have the same name, saved in different folders; and I've tried to merge them using batch with **Merge tool** in ArcGIS... but inexplicably ArcGis say that _cannot process file with the same name_! I'Now it's clear that I've too many files to rename them. I've then used QGis and the process works, but lose the elevation information (Z) that I need. I've seen other analogue questions but nothing directly usable. Some suggestions for a simply solution?


title: GeoServer - OSM .shp file and aerial shifted

text: I have a problem with two layers that I have imported into geoserver. One is retrieved via WMS from a Geoserver (aerial tiles) I don't have access to (I'm just allowed to use it via WMS) the other one is a shapefile retrieved from geofabrik and uploaded to my geoserver. The Problem is that when I create a grouplayer with both layers they are somehow shifted (pls see picture attached). Both layers use the same Projection and the shift seems to be the same everywhere. Is there a possibility to "move" the shape layer a little bit - or how can I get the roads from the tiles to fit inthe pictures. Btw: I do know that there were similar questions here unfortunately I could not find an answer to my problem. Thx a lot ![enter image description here](http://i.stack.imgur.com/PiM1Z.jpg)


title: Processing HDF4 data in R

text: I am trying to process MODIS files in R, which are in HDF4 format (link to data). I have experience reading in HDF5 and netCDF files in R, but HDF4 seems to be a common problem in R. I am using a MAC (Mountain Lion OS X 10.8.5), and have found the following webpages (link1, link2) helpful in setting up GDAL and the r package, `rgdal`. My problem is that is doesn't look like I have any HDF4 drivers, and I'm wondering what step I have missed?               library(rgdal)     gdalDrivers()          ...     ...     44              GTX                                 NOAA Vertical Datum .GTX   TRUE FALSE     45              GXF                             GeoSoft Grid Exchange Format  FALSE FALSE     46             HDF5                       Hierarchical Data Format Release 5  FALSE FALSE     47        HDF5Image                                             HDF5 Dataset  FALSE FALSE     48              HF2                               HF2/HFZ heightfield raster  FALSE  TRUE     49              HFA                              Erdas Imagine Images (.img)   TRUE  TRUE     ...     ...      So, nothing listed for HDF4 as shown in point 6.2 of the example. Does anyone know what my problem might be? The issue of reading HDF4 data into R seems to be a common one. I would have no problem doing some sort of batch conversion to HDF5 or netCDF if someone has a suggestion for this.


title: Negative Values on Vegetation Index over water, exclusion of pixels?

text: I am having problem with som buffer zones where I try to extract NDVI's from MOD13files on coordinates that are closeby water. My buffer zones around these coordinates include water annd therefore gives me negative mean values in these areas. Is there a way of getting around this problem by excluding pixels, eliminating negative values from the NDVI? I am using ArcGis9.2 and I am kind of new to this although with a tight schedule so any help would be greatly appreciated! Thanks!


title: why distanceTo gives different results?

text: i am trying to get distance between two polygons and somehow i get different distance, depending on how I call distanceTo a.distanceTo(b) != b.distanceTo(a) should it be like that?               layers[0].features[0].geometry.distanceTo(         layers[0].features[1].geometry,         {details:true,edge:true}     );     Object {distance: 10.48962511638198, x0: 2684839.287184146, y0: 7750111.194373621, x1: 2684849.7767009, y1: 7750111.146694}      and               layers[0].features[1].geometry.distanceTo(         layers[0].features[0].geometry,         {details:true, edge:true}     );     Object {distance: 19.05365206557588, x0: 2684857.7419808176, y0: 7749989.676176046, x1: 2684838.7291617, y1: 7749988.4294338}      **Some additional information...** Main map object:               var map = new OpenLayers.Map('map', {         units: 'm',             projection: 'EPSG:3857',             layers: [                 new OpenLayers.Layer.Google(                     "Google Roadmap",                     {type: google.maps.MapTypeId.ROADMAP, numZoomLevels: 22}                 ),                 conflictLayer,markerLayer,             ],             center: new OpenLayers.LonLat(1787495.895083,  5755860.4398742),             zoom: 18         });      later I just add more layers and keep these layers in array, so 'layers[0]' is one of them. code to generate new layer:                 var newLayer = new OpenLayers.Layer.Vector(o.name, {           styleMap: styleMap,           renderers: renderer,           zone_id:o.id         });         map.addLayer(newLayer);         layers.push(newLayer);


title: MOD11A2, falsely negative nightly LSTs

text: I am having troubles with extracting data in ArcGIS buffer extractions from the MOD11A2files after processing them into GEOtiff, most of them seem fine but when looking at nightly temperature on some of my coordinates of interest my values go way below zero. Some down to -273. I have multiplied the data with 0,02 and then subtracted 273 to get them as Celsius. Is this a common problem? I can't manage to find anything about it in forums etc and unfortunately I have a tight schedule and am stressed to get the results out. Anyone familiar with this problem?


title: Problem importing ESRI shapefile information to MySQL

text: I recently updated to Qgis 2.0 and therefore the "exporttomysql" plugin does not work anymore. What I want to do is exporting the geometry information of a shapefile(Polygon) into a MySQL database. I found some information on how to that at http://www.bostongis.com/PrinterFriendly.aspx?content_name=ogr_cheatsheet and adopted their code to fit my needs               ogr2ogr -f "MySQL" MYSQL:"CSV_DB,host=myhost,user=user,password=pwd,port=3306" -nln "world" -a_srs "EPSG:4326" ~/Documents/Fields.shp       However, I get the error "Unable to open the datasource ~/Documents/Fields.shp with the following drivers.... When I open the file with "ogrinfo" it tells me the the opening with the driver "ESRI Shapefile" was successful. 1; FieldsFinal(Polygon) Any ideas how I get the data into my database? Thanks


title: Best format for maps on mobile device with routing functionality

text: I am trying to choose the best file format for maps that meets these requirements:   1. There are command line tools to convert from .osm data   2. Are compact enough for mobile devices and internet download   3. Have the ability to be fed into a routing engine. I am aware of .shp, geojson, spatialite and the proprietary ctm1. Any other formats I should check out or opinions on the ones I mentioned?


title: How to "spatially normalize" rasters of different extents to sum values of overlapping cells

text: How do I add rasters together that are of different spatial extents? I have many, many features that overlap over a large area (entire east coast) and I am looking to sum them "stacked" where they overlap. My strategy is to: 1) Convert each feature to a raster (using Polygon to Raster tool) 2) Set Null values for feature to 0 using Con(IsNull("raster"),0, "raster") in Raster Calculator 3) "Normalize" each raster (spatially) to have matching spatial extent (of the entire study area) - do this by adding each raster to the "Constant Raster" I created, which has the study area's extent, the same projection, and the same cell size of the individual rasters. The result _should_ be a new raster to the extent of the study area, with the same correct values at the location of the converted-feature in question. 4) Use Raster Calculator to simply add the normalized rasters together. Something's not working at step #3. When I add the projected, constant raster (all values = 0) to the individual feature-converted-to-raster (with null values set to 0) in the Raster Calculator Tool, with Processing Extent set to that of the Constant Raster, the resulting raster output is of the correct extent but the values of the feature-converted-to-raster are incorrect. Instead of an expected range of, for example, [0, 77] (where the feature that was converted had the value 77, and any "background" value = 0) the value is something extreme, such as [-1.79e+308, 1.79e+308] _[Edited to explain unexpected output]_ Can someone explain what it is I am missing here? My searches on StackExchange only revealed this question but I must have missed something. Thank you.


title: Using Labels from Maplex engine inside ArcGIS Explorer?

text: I'm trying to display data in the upper left and upper right corners of this parcel data. I can get it to work inside ArcMap with the Maplex engine, but when I save a layer file and open it in arc explorer it puts the two labels back in the middle. Any ideas how I can get those to stay in the corner of each parcel? ![enter image description here](http://i.stack.imgur.com/p2kby.png) ![enter image description here](http://i.stack.imgur.com/96bvD.png)


title: How to compare two Geodatabases

text: I would like to compare the data of two geodatabases, in SDE or a file geodatabase, to detect the changes occured between them , i've tried to compare the Xml workspace document of both databases using notepad++ "compare plugin" , but it's not easy to check and interpret the differences. i'm looking for a better solution.


title: Changing user's password and adding users in Geoserver using REST service

text: I want to change user's password (e.g. for `admin` user) and create a new one using `Geoserver`'s `REST` `API`. Is this possible? I tried to look for some solution in their DOCUMENTATION, but nothing found about users and passwords...


title: Finding adjacent tiles in Google Maps

text: I built a NodeJS API that returns locations and returns paths and stores them locally in MongoDB for fast responses. Also another part of the API is that it returns compositions of tiles (2-dimensional grids, between two points). I got it to work quite nicely, only thing I'm struggling with is the vertical adjacent tile on zoom levels lower than 10. I made quite a filthy way of finding the adjacent tile I just found out, and I'm also not really sure why it works so good on zoom levels that are 10 or higher than 10. Weirdly enough all the horizontal tiles seem to line up perfectly around the world. This is what one of the inefficiencies looks like: ![inefficiency](http://i.imgur.com/u6Cf9wt.jpg) I kind of know what I'm doing wrong (I think, not the best at calculating with longitudes and latitudes), but I have to take into account that the earth's radius is becoming more when you go further down to the equator. And I'm not taking that into account with my calculations. This is how I calculate my adjacent tile:               adjacent: function (α, width, height, zoom) {         this.bounds(width, height, zoom);              var sideX = Math.abs(this.geometry.bounds.northWest.lon - this.geometry.bounds.northEast.lon),             sideY = Math.abs(this.geometry.bounds.northWest.lat - this.geometry.bounds.southWest.lat);              α = Math.abs(α) % 360;         var αMax = Math.atan(sideY / sideX),             opposite, adjacent, lat, lon;              αMax = αMax.toDegrees();              // Right side of tile         if (α < αMax || α >= (360 - αMax)) {             opposite = Math.tan(α.toRadians()) * sideX;             lat = this.lat + opposite;             lon = this.lon + sideX;         }              // Bottom side of tile         if ((180 - αMax) > α && α >= αMax) {             adjacent = sideY / Math.tan(α.toRadians());             lat = this.lat + sideY;             lon = this.lon + adjacent;         }              // Left side of tile         if ((180 + αMax) > α && α >= (180 - αMax)) {             opposite = 2 * (Math.tan(α.toRadians()) * (sideX / 2));             lat = this.lat + opposite;             lon = this.lon - sideX;         }              // Top side of tile         if ((360 - αMax) > α && α >= (180 + αMax)) {             adjacent = sideY / Math.tan(α.toRadians());             lat = this.lat - sideY;             lon = this.lon + adjacent;         }              var adjacentPoint = new Point([lat, lon]);              return adjacentPoint;     }      And this is how I calculate my bounds:               bounds: function (width, height, zoom) {         // Calculate distances to adjust by         var Δx = Math.round(width / 2), Δy = Math.round(height / 2);              // Seed wind quarters         var quarters = {             north: this.adjustLatByPixels(Δy * -1, zoom),             south: this.adjustLatByPixels(Δy, zoom),             west: this.adjustLonByPixels(Δx * -1, zoom),             east: this.adjustLonByPixels(Δx, zoom)         };              // Calculate geometry         this.geometry = {             quarters: quarters,             bounds: {                 northWest: new Point([quarters.north, quarters.west]),                 northEast: new Point([quarters.north, quarters.east]),                 southWest: new Point([quarters.south, quarters.west]),                 southEast: new Point([quarters.south, quarters.east])             }         };              return this.geometry;     }      Hopefully there is a geography expert around that can help me with this problem.


title: How to delete individual labels?

text: How do I delete an individual label without converting the entire data frame into annotation?


title: OpenLayers GeoJson style

text: I have a GeoJson that I'm displaying in OpenLayers. I created a style map and added it to a layer. The problem is that in GeoJSON I have two types of features - Points and Linestring. And they both are rendered the same. I would like to render points in one style, and linestring as another. How can I do this?


title: System variable PATH overwritten in Qgis

text: The system env variable 'PATH' seems to be overwritten in Qgis, replacing it with their own values. Not sure if it was intended or not. Anyway in version 2.2, I could solve it by go into settings and appending the PATH variable with the values copied from system properties => environment variables. Apparently for version 2.4, the above solution does not work anymore and Qgis will always overwrite the PATH variable with their own values. I needed the PATH variable as my project requires the execution of external programs from Qgis by making use of python subprocess module. Now in v2.4, external programs cannot be executed. I do not wish to define absolute path for all the external programs as the project will frequently be ported to different systems. Any solutions or feedback to this will really be appreciated!


title: How to make an editable City Map from OSM data?

text: I want to make a map of Cardiff, Wales, showing features useful to tourists, e.g. roads with road names, parks, museums, other POIs, etc. I want the final map to be clear to read - simple roads and names, basic outlines/boundaries - so editing (moving/manipulating) the data from its true georeferenced position would be ideal. Using data from OpenStreetMap, I can't work out how to make a simple map in QGIS, and using Maperitive I can't find a way to edit the data from its original, simple view, despite being vector data. A lot of the information online seems to be quite old, for example the OSM plugin is no longer available in the latest version of QGIS. Is anyone able to offer advise on making a simple, editable map of a city?


title: What is the relationship between Autodesk and the MapGuide Open Source project?

text: In the year 2006 Autodesk started MapGuide as an open source project (MGOS). At the same time they began to develop the commercial version of MapGuide (MapGuide Enterprise MGE /Autodesk Infrastructure Map Server AIMS). With respect to the software development: What is the relationship between Autodesk and the MapGuide Open Source project today (and in the future?)? How does Autodesk support the MapGuide Open Source community today? Are developers from Autodesk writing code for MGOS? How benefits MGOS from new features of AIMS? Regards Jens


title: How to produce a singlepart featureclass with unique IDs for each spatially non-connected polygon?

text: What is the best method to convert a singlepart shapefile (i.e. one record that contains several non-connected polygons) into a mutlipart shapefile? --EDIT For clarification the original question asked how to produce a singlepart featureclass with unique IDs for each spatially non-connected polygons. The explode method in advanced editing or Multipart to Singlepart would give this result...


title: How to find the Shortest Path between two coordinates points using pgrouting A-Star?

text: In PostGIS I have two tables: an osm road network and a table that contains the start and end coordinate points of each path (i.e. userid, start_x1, start_y1, end_x2, end_y2). I would like to know: 1) how to create an osm road network that can be used for implementing the Astart algorithm with pgrouting (I actually have an osm road network which I obtained using osm2po and that works fine Dijkstra shortest path algorithm using pgrouting layer plugin) 2) how to find the shortest path of each start and end coordinate points using A-star. Thank you, Selini


title: How to show a popup on mouse-over, not on click?

text: Is it possible in Leaflet that popup opens on mouse-over, not on click? This is working for just one marker at a time, but I need it for a bigger number of markers:               marker.on('mouseover', function(e){         marker.openPopup();     });


title: Leaflet geoJson in front of markers in a FeatureGroup

text: I am loading some markers and json data to a leaflet map. I need to get my json in front of the markers. I have tried placing the markers and json in seperate FeatureGroups and then tried bringToFront() and bringToBack() to no avail               var markerHolder = new L.FeatureGroup().addTo(map).bringToBack();     var jsonHolder = new L.FeatureGroup().addTo(map).bringToFront();      I also tried using this panes work around but when using a Feature group in place of the tileLayer as in the example below, an error was thrown at getContainer()               var topPane = map._createPane('leaflet-top-pane', map.getPanes().mapPane);         var topLayer = L.mapbox.tileLayer('lxbarth.map-vtt23b1i').addTo(map);         topPane.appendChild(topLayer.getContainer());         topLayer.setZIndex(9);      Here is a fiddle illustrating my problem - could anyone demonstrate how to get my json features in front of the markers?


title: GRASS GIS using a loop with r.neighbors command

text: I'd like to do a **trend analysis with ndvi data** from 1982 til 2005. I am using the `r.neighbors` command to get first the spatial variance (then I would use `r.series` method=slope to get the actual trend) and in order to speed things up I'd like to do it using a loop... however it is not running :(   Here is my code... can someone help me?               for a in        1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993        1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 ; do            for i in 311 ; do         r.neighbors           input=AVHRR_NDVI_$aDOY$i           output=AVHRR_NDVI_$aDOY$i.sp.var           method=variance           size=5          echo "loop ends" ;        done ;      done


title: How to publish a layer with standalone geowebchache in geoserver?

text: I have a multiple geoservers on two different systems with same port(one is backup for another) and geowebcache in on of the system with different port . Now i want to publish my layers in geoserver with standalone geowebcache. is there any material to refer? Please help me?


title: ECW on Mac Osx 10.8 Mountain Lion

text: I'm having some problems with GdalTools on QGIS 1.8, with MacOsx 10.8: I've the GDAL complete library 1.9 but ECW is still a not supported format. What should I do? Thanks


title: Can QGIS give GeoJSON layers a meaningful default name?

text: When I add a GeoJSON layer to my QGIS project, it always gets the name `OGRGeoJSON`. This is in contrast with shapefile layers, which seem to be named after the file name by default. Is there a way to either configure QGIS or change my GeoJSON files so that they get a meaningful default name when I first add them to a project?


title: Has anyone compiled QGis for 64-bit Windows 7?

text: It seems to me that a 64-bit build would be a much desired version. Has it been done or are there Qt 64-bit compiler issues blocking the way?


title: OpenLayers Vector Layer callback does not execute with GeoJSON

text: Why does this callback get ignored:               new OpenLayers.Layer.Vector("My layer",           {             projection : "EPSG:4326",             strategies : [ new OpenLayers.Strategy.Fixed(), refreshStrategy ],                  protocol : new OpenLayers.Protocol.HTTP(             {                           url    : "myphpurl.php",               params : { rand : Math.random() * 1000 },               format : new OpenLayers.Format.GeoJSON(),               callback: function(){alert('test');}             }),      I'm trying to get access to how many features were returned...


title: Fast data representation on the vector layers

text: I have found an interesting article about how `HTML5` supports representation on large datasets. Article can be found here. In the article it is said that `HTML5` and feature simplification can increasingly improve visualization speed. Does anybody knows about the source code (or just explain what is needed to be done to reach such results) or suggest how `OpenLayers` supports feature simplification?


title: Why is MapServer adding horizontal artifacts when projecting EPSG:4326 data to EPSG:3857?

text: We're running MapServer with a map that has just one layer, political boundaries from a SHP file. This layer is using a basic EPSG:4326 projection, and the map itself is using EPSG:3857. The resulting images generated by MapServer CGI queries have horizontal lines at certain latitudes: ![enter image description here](http://i.stack.imgur.com/a8p69.png) Beyond this problem, the map is okay and lines up perfectly with the Open Street Map in OpenLayers. The MAPFILE we're using:               MAP       WEB               METADATA             "wms_title" "Map Server"             "wms_enable_request"  "*"             "wms_srs"             "EPSG:3857"                END       END            PROJECTION         "init=epsg:3857"       END            IMAGETYPE      PNG       EXTENT -20037508.34 -20037508.34 20037508.34 20037508.34       SIZE           400 300       SHAPEPATH      "../data"       IMAGECOLOR     71 245 242            LAYER # States polygon layer begins here         NAME         Borders         DATA         ne_10m_admin_1_states_provinces_shp         STATUS       OFF         TYPE         POLYGON            PROJECTION         "init=epsg:4326"        END              CLASS           NAME       "Borders"                STYLE             COLOR        43 156 51             OUTLINECOLOR 32 32 32           END         END       END # States polygon layer ends here          END       I've verified that the proj.4 EPSG file has the proper definition of EPSG:3857, and don't understand what's happening. Any insight would be greatly appreciated.


title: How to use Cross Origin Resource Sharing (CORS) from ArcGIS API Javascript Application?

text: I am getting the CORS error                 XMLHttpRequest cannot load http://services.arcgisonline.com/ArcGIS/rest/info?f=json.        No 'Access-Control-Allow-Origin' header is present on the requested resource.        Origin 'http://myserver.com' is therefore not allowed access.      when I use the ESRI topo base map in my ArcGIS JSAPI application. The application is hosted on IIS8 (Windows Server 2012 R2) I am reading the ESRI ArcGIS API for JavaScript Samples CORS buffer There is a suggestion to add the root URL for CORS enabled servers to               esri.config.defaults.io.corsEnabledServers      Can you provide example how to do that? I am not sure if the problem is on my server/application or the ESRI map server


title: Can you run GIS Squirrel next to SDE?

text: My company has SDE to run our webmapping system. We want our arcbasic license users to be able to edit this data without dropping 30k for 6 users on upgrading to standard. Can we run squirrel next to SDE to edit the same database? I'm assuming we can't, but wanted to be sure. I'm using SQL server 2008.


title: is MySQL spatial ready for production?

text: currently I'm using PostgreSQL/PostGIS for a small spatial project. Last week the project manager asked me if we could migrate to MySQL with Spatial Extension. I know MySQL but have no experience with the Spatial Extension. I heard that the extension is not very fast and might contain some bugs. Now my question - is MySQL Spatial extension stable and fast enough to be considered for an online database with about 20000 spatial operations a day? The operations are just the standard stuff like: - export to KML - point in polygon - overlapping polygon The database has one GB of spatial data. Best Regards, D3


title: Will QGIS 2.0 work on Mac 10.8.5

text: I've downloaded QGIS 2.0 and have opened a project. All the layers are there, but it will not show the selections on the map. The attribute table shows the selected items. Am I missing something or is there some other cause?


title: How to interpolate a point vector layer inside a boundary (polygon layer) QGIS/GRASS/SAGA?

text: I have a to create a heatmap of water temperatures inside a lake for which I have sample points. The map should look much like the climate maps generated by climate change models. I have the temperature points and they were attached as attributes to a lat/long point shape file layer. I would like to interpolate these points using temperature as an interpolation attribute, but I'd like to contain the values within the lake. I would prefer not to do an interpolation with the lake area polygon as a masking/clipping layer, but rather the interpolation should take in account the boundaries much in the way ANUDEM seems to do it for hydrological maps, if my understanding is correct. That is borders should influence the value of the interpolated value. Please let me know if such interpolation is possible in one of: QGIS, GRASS, SAGA. If ArcGIS can provide such functionality while the others don't I can get a temporary license to do this. Thanks Edit: After doing more research I think that what I need is called interpolation with fault lines or barriers. I found an Open Source program _surfit_ , which apparently does that. I am curios if in any of the software mentioned above there is something that does that?


title: Recommended programming language for remote sensing?

text: I am beginning studies that hopefully will lead to a long career as a remote sensing specialist. I am currently working with ArcGIS for some applications and learning ENVI for others. I have realized that it is imperative that I learn a programming language, and am stuck facing a choice between IDL and Python. I would love to hear which programming language the community recommends for the remote sensing of land surface processes. Thanks! Emily


title: Get coordinate of second to last vertex for python ArcGIS

text: I need to obtain the x and y attributes of the second to last vertex for each segment and add to the attribute table using Arcpy (or some other python module).   To get the Y co-orodinate of the first vertex I can use !Shape.FirstPoint.Y!               arcpy.CalculateField_management(Flow_Paths, "Y_Cord_A", "!Shape.FirstPoint.Y!", "PYTHON_9.3", "")      In QGiS I can simply call Xat(-2) and Yat(-2) to get the second to last co- ordinate, is there an equivalent command for ArcPy?


title: In Python, reading a GDAL raster from memory instead of a file

text: I'd like to analyze some raster data in Python using NumPy. This data doesn't necessarily exist as a file on the filesystem, it could be the result of a query against PostGIS, or something transferred over the network. Is it possible to get GDAL to read data from a buffer in memory, instead of requiring a filename? I realize one solution could be to open up a named temp file, write the data to that, and pass it as a filename to GDAL, but that seems inefficient, and more likely to incur I/O bottlenecks. Additionally, in the case of PostGIS, it would be nice to read the native format directly, as opposed to converting it to some other file format and reading that. Edit: to clear up a bit of confusion, just using the PostGIS driver directly won't work either, since the query could be for a raster analysis, rather than just data stored in the DB.


title: Debugging errors of SAGA interpolation algorithms in QGIS 2.0 Processing

text: I encounter recurrent problems with the geoalgorithms that perform interpolations. No matter technique I use I end with the following errors. As state in OGR tab QGIS is unable to open a file in C drive. I have very little IT background so the message error is meaningless. Is there anybody able to decrypt the message and solve the problem? thank you ![enter image description here](http://i.stack.imgur.com/BlLnX.jpg) ![enter image description here](http://i.stack.imgur.com/gjkqe.jpg) here this is the algorithm i used ![enter image description here](http://i.stack.imgur.com/kfYAk.jpg)


title: problem on loading kriging outputs

text: My problem is the following: everytime i run ordinary kriging or other interpolation techniques i obtain this error message from QGIS ![enter image description here](http://i.stack.imgur.com/37v9X.jpg) Notice that these are the details for SAGA ![enter image description here](http://i.stack.imgur.com/sLxaP.jpg) I would like a clarification about: why i am getting this error message and what i should change (path were SAGA giu is located?) thank you


title: How to make multiple qgis mapbook pages in 1 pdf?

text: Is it possible to export a mapbook in QGIS 2 with multiple pages all in one PDF? Thanks


title: ArcGIS 10.1 calculating area out of Aspects

text: I have used the aspect tool to visualize the 8(+ Flat) slopes of the different directions, now I would like to get the numbers of the areas for each different direction (+Flat). Anybody knows how that works?


title: Creating an OVF file for connecting Oracle and QGIS

text: I'm trying to create an OGR file that can be opened in QGIS that pulls in vector data from an Oracle Database. I've tried a few combinations of things but can't seem to get anything to work. I've few questions if anyone could help:   1. Does the file need to be .ovf or a .vrt? Does this matter?   2. Can the OGR file use oracle spatial geometry or do you have to encode it from data in the columns?   3. The connection string I've used is in the format username/password@databasename. I have Oracle instant client installed (which works correctly with other applications). Is this all you need in the connection string or do you need the full URL + Port number for the database and also the schema?  Thanks very much for any help!               <OGRVRTDataSource>         <OGRVRTLayer name="TEST_LAYER_NAME">             <SrcDataSource>OCI:testuser/testpass@testDB</SrcDataSource>              <SrcLayer>TEST_POINTS</SrcLayer>             <GeometryField encoding="PointFromColumns" x="EASTINGS" y="NORTHINGS"/>         </OGRVRTLayer>     </OGRVRTDataSource>


title: ArcGIS 10.2 Extract Polyline from Start/End Coordinates

text: I have two files. A polyline shape file of highways and a spreadsheet that contains beginning and ending coordinates w/ descriptions of segments of the highways. For example:               BegX, BegY, EndX, EndY, DescriptionField1, DescriptionField2     -80.0, 30.0, -80.1, 30.1, Hwy1Partial, 1     -80.1, 30.1, -80.2, 30.2, Hyw1Partial, 2      I need to create a new polyline shape file that breaks up the hwy segments from the original polyline file into a new one that has the beginning and ending coordinates from the spreadsheet. I've attempted to join the data but can't figure out how to get it to extract the polyline that starts at the begin coordinate and ends at the end coordinate. I've used the XY to Line tool, but that just gives me a line connecting the begin coordinate to the end coordiante. I need to match it up with the polyline from the hwy shape file. Thank you, any help is appreciated. **Edit/Update:** The shapefile I am using does not have M values. The spreadsheet I have contain latitude/longitude coordinates that look to have originally came from the shape file (they fall directly on top of it when plotted). I need a polyline file that matches the spreadsheet exactly. Meaning the spreadsheet would become the new shape files attribute table, and each record would represent a polyline that was extracted between the beginning coordinate and end coordinate. ![Example](http://i.stack.imgur.com/ubizn.jpg) The black line represents the shape file that I have. The two coordinates are examples of beginning and ending lat/lon coordinates that are in the spreadsheet. What I want is to extract the blue line from the shape file into a new polyline file with the attributes of the original spreadsheet. I have played around with making XY point layers out of my beginning and ending coordinates then splitting the file based on those new points. The problem is that my shapefile has more lines than I actually need. The resulting file doesn't split it completely the way I need, and I cannot find a way to relate my attributes from the original spreadsheet back to it. What I really need is a custom query that takes in a begin lat/lon, then an end lat/lon and queries the shapefile to get the resulting line. After it finds the new line it would take it to a new shapefile; copying the attributes over from the spreadsheet to be associated with the newly created line.


title: How do I change the colors of my layer file?

text: I have a ILayer which I load from a .lyr file which I created from a shapefile. How can I change the colors of the elements? I want to make the roads red instead of green for example. I need to do this programatically with ArcEngine 10 and C# VS2010.


title: How to get count of non-NA raster cells within polygon using R

text: I've been running into all sorts of issues using ArcGIS ZonalStats and thought R could be a great way. Saying that I'm fairly new to R, but got a coding background. The situation is that I have several rasters and a polygon shape file with many features of different sizes (though all features are bigger than a raster cell and the polygon features are aligned to the raster). I've figured out how to get the mean value for each polygon feature using the raster library with extract:               #load packages required     require(rgdal)     require(sp)     require(raster)     # ---Set the working directory-------     datdir <- "/test_data/"          #Read in grid of water depth     ras <- raster("test_data/raster/pl_sm_rp1000/w001001.adf")          #read in polygon shape file     proxNA <- shapefile("test_data/proxy/PL_proxy_WD_NA_test")           #calc mean depth per polygon feature     #unweighted - only assigns grid to district if centroid is in that district     proxNA$RP1000 <- extract(ras, proxNA, fun = mean, na.rm = TRUE, weights = FALSE)          #plot depth values      spplot(proxNA[,'RP1000'])      The issue I have is that I also need an area based ratio between the area of the polygon and all non NA cells in the same polygon. I know what the cell size of the raster is and I can get the area for each polygon, but the missing link is the count of all non-NA cells in each feature. I managed to get the cell number of all the cells in the polygon `proxNA@data$Cnumb1000 <\- cellFromPolygon(ras, proxNA)`and I'm sure there is a way to get the actual value of the raster cell, which then requires a loop to get the number of all non-NA cells combined with a count, etc. BUT, I'm sure there is a much better and quicker way to do that! If any of you has an idea or can point me in the right direction, I would be very grateful! Many Thanks!


title: How do you publish map caches from ArcGIS Data and Maps for Server?

text: I haven't published ESRI-supplied data before, but I followed the instructions that came with the Data and Maps for ArcGIS Server 10. When I view the supplied MXD file for a map cache, it's just a black background. The help docs say it's supposed to be different from what you expect, so I followed the rest of the steps:   1. Copy map cache folders to the arcgisserver\arcgiscache folder.   2. Copy the map documents to a shared location accessible to your GIS server. (Also copy the data folder from the documents folder).   3. Create a new map service using the MXD as the source document. I followed the steps that come with the data but the service never displays any data. I suspect that the imprecise step of "Copy the map documents to a shared location accessible to your GIS server" requires something I haven't done. Has anyone published map services from the ArcGIS data discs before?


title: Spatial database containing coordinates for stars in the sky?

text: My goal is to map the Skydome ( stars in the sky projected in earth surface ) using GIS. Has anyone done this before? I can't find any project close to this. What I just need to know is if there is a dataset -not the paper or software, just the points (geometry) - of the relative position of the stars in the sky. I mean if I look at the sky I can see a kind of surface where the stars are. This is the Skydome. If this surface is over my table, I have a kind of "MAP" of the sky of my hemisphere. If I have both hemispheres I have the Sky Map of the stars. ![enter image description here](http://i.stack.imgur.com/vMJdA.gif)


title: Formula for Enhanced Vegetation Index (EVI) in ERDAS modelbuilder?

text: I am trying to derive Enhanced Vegetation Index (EVI)from Landsat 5 TM data using ERDAS model builder.d I can't seem to get the output I want as it comes out blank. This is the formula I used 2.5 * (("%Band 4%" - "%Band 3%") / ("%Band 4%" + 6 * "%Band 3%" - 7.5 *"%Band 1%" + 1))can someone assist me with the exact formula i need to type into ERDAS model builder. Does anyone see something in here that doesn't make sense-- or even better--a formula I could copy/paste?