title: make a selection by attributes of the fourth part of a range of numbers with arcpy

text: I am using Arcmap 10.1 and I need to create a script (for use on shapefiles) where I have a input field named "sequence" which is numerical data. I need to select the fourth part of a range of numbers. For example, 1-100 depending on the amount of data in the shapefile and I need to make a selection by attributes of "sequence" where I divide 100/4 = 25 and the selection is the first 1-25, the second selection of 26-50, the third 50-75, etc. The number of records can change and be 150/4 = 37 and the selection would be 1-37 the second third etc 38-75![enter image description here](http://i.stack.imgur.com/KrLZL.png) I have the ide of using arcpy.GetCount_management for me to total the records but do not know how to use the resuldato to divide by 4 and then put the result and make the selection with arcpy.SelectLayerByAttribute_management but not how to do this and if this is the right way to do


title: How to copy or export list filter from one model to another?

text: I have a model tool with a variable parameter that is filtered by a string list (county names) I manually created this list filter when I created the model tool. I am now creating a new model tool and would like to use the same list of county names is there anyway to import the list filter from my other model so I do not have to manually create the list again. I am an amateur at best with python so if it is possible to keep the solution within the realm of ModelBuilder that would be best but if python is neccessary bear with me. This is the value list I would like to reuse ![enter image description here](http://i.stack.imgur.com/5qiq2.png)


title: How to get the feature set for geoprocessing service?

text:
I am calling a geoprocessing service in a web app. The geoprocessing service has one parameter (a feature set). But I am not sure how to find the parameter in the javascript code to send it to the geoprocessing since I am using an editor widget to first create the polygon and then run the geoprocessing service. Can you guys take a look at my code and tell me how I reference a feature set from the editor widget results? Thank you!               ], function(         Map, BasemapGallery, arcgisUtils, Edit,          ArcGISTiledMapServiceLayer, FeatureLayer,         Query, QueryTask,         SimpleMarkerSymbol, SimpleLineSymbol,          Editor, TemplatePicker,          esriConfig, jsapiBundle,         arrayUtils, parser, keys,                borderContainer, contentPane, titlePane,         dom, on       ) {         parser.parse();                     // snapping is enabled for this sample - change the tooltip to reflect this         jsapiBundle.toolbars.draw.start = jsapiBundle.toolbars.draw.start +  "
Press ALT to enable snapping";              // refer to "Using the Proxy Page" for more information:  https://developers.arcgis.com/en/javascript/jshelp/ags_proxy.html         esriConfig.defaults.io.proxyUrl = "/proxy";                       esriConfig.defaults.geometryService = new esri.tasks.GeometryService("http://tfsgis-iisd01:6080/arcgis/rest/services/Utilities/Geometry/GeometryServer");              map = new Map("map", {            basemap: "hybrid",                center: [-98.57, 30.98],           zoom: 6,           slider: true          });              map.on("layers-add-result", initEditor);                   var ActivityArea = new FeatureLayer("http://tfsgis-iisd01:6080/arcgis/rest/services/MyMapService2/FeatureServer/1",{           mode: FeatureLayer.MODE_ONDEMAND,            outFields: ['*']         });              var ActivityPoint = new FeatureLayer("http://tfsgis-iisd01:6080/arcgis/rest/services/MyMapService2/FeatureServer/0",{           mode: FeatureLayer.MODE_ONDEMAND,            outFields: ['*']         });              var stewardship = new FeatureLayer("http://tfsgis-iisd01:6080/arcgis/rest/services/MyMapService2/FeatureServer/2",{           mode: FeatureLayer.MODE_ONDEMAND,            outFields: ['*']         });              map.addLayers([ActivityArea,ActivityPoint,stewardship]);                             function initEditor(evt) {                var templateLayers = arrayUtils.map(evt.layers, function(result){             return result.layer;                });                     var templatePicker = new TemplatePicker({             featureLayers: templateLayers,             grouping: true,             rows: "auto",             columns: 3           }, "templateDiv");           templatePicker.startup();                var layers = arrayUtils.map(evt.layers, function(result) {             return { featureLayer: result.layer };           });                var settings = {             map: map,             templatePicker: templatePicker,             layerInfos: layers,             toolbarVisible: true,             createOptions: {               polylineDrawTools:[ Editor.CREATE_TOOL_FREEHAND_POLYLINE ],               polygonDrawTools: [ Editor.CREATE_TOOL_FREEHAND_POLYGON,                 Editor.CREATE_TOOL_CIRCLE,                 Editor.CREATE_TOOL_TRIANGLE,                 Editor.CREATE_TOOL_RECTANGLE               ]             },             toolbarOptions: {               reshapeVisible: true             }                     };                                    var params = {settings: settings};               var myEditor = new Editor(params,'editorDiv');           //define snapping options           var symbol = new SimpleMarkerSymbol(             SimpleMarkerSymbol.STYLE_CROSS,              15,              new SimpleLineSymbol(               SimpleLineSymbol.STYLE_SOLID,                new Color([255, 0, 0, 0.5]),                5             ),              null           );           map.enableSnapping({             snapPointSymbol: symbol,             tolerance: 20,             snapKey: keys.ALT           });                myEditor.startup();                        }                         on(dom.byId("button"),"click", function() {         gp = new esri.tasks.Geoprocessor("http://tfsgis-iisd01:6080/arcgis/rest/services/CalcFeaturesStewardshipOneParameter3/GPServer/CalcFeaturesStewardshipOneParameter")         //gp.outSpatialReference = map.spatialReference                   var features= [];              var selectGraphics = featureLayer.getSelectedFeatures();              features.push(selectGraphics);         var featureSet = new esri.tasks.FeatureSet();              featureSet.features = features;                   var params = {Stewardship:featureSet}              gp.submitJob(params, completeCallback , statusCallback)                        });              function completeCallback(jobInfo){            }            function statusCallback(jobInfo) {         var status = jobInfo.jobStatus;         if(status === "esriJobFailed"){           alert(status);              }         else if (status === "esriJobSucceeded"){           alert(status);         }       }


title: New to ArcGIS and want to compute distance between two points "as the fish swims"

text: We have a layer named Hydrography which is the water area for a given assignment. I have a station point and a target point. I need to compute the distance between those two points but inside of the hydrography. (boats cannot go over land) The DefenseService distance computation looks passable but it seems to be in a straight line since it does not take any Raster as input. The Euclidean seems to almost be there but I cannot find anything that is straight forward. I am using ArcGIS Explorer 2500, ArcGIS for Server 10.1 SP1 Enterprise, and we are writing our services and tools in C#. I have the loading of the hydrography, the creation of the IPoints, but at a loss as to what to do next! Any help appreciated, please.


title: Can't get MapServer WFS to display in OpenLayers

text: I'm new to OpenLayers and MapServer and I'm trying to serve a shapefile through a WFS but can't seem to get the data to display. When I use Firebug and HttpFox in Firefox I seem to be getting some data coming back, but it just won't get rendered on my web page. It is entirely possible that I have missed something simple but I can't see what. I've got the following in terms of my setup. MapFile:               MAP          NAME "MapServerWFS"     STATUS on     SIZE 600 600     EXTENT 0 0 700000 1300000     UNITS METERS     SHAPEPATH "shp"     IMAGECOLOR 240 240 240      PROJECTION          "init=EPSG:27700"     END          WEB         TEMPLATE "template.html"         IMAGEPATH "C:/DevProjects/fgsv2/output/"         IMAGEURL "../output/"         METADATA           "wfs_title"          "WFSDemo" ## REQUIRED           "wfs_onlineresource" "http://localhost/cgi-bin/mapserv.exe?map=C:/DevProjects/fgsv2/mapfile.map&" ## Recommended           "wfs_srs"            "EPSG:27700" ## Recommended           "ows_srs"            "EPSG:27700" ## Recommended           "wfs_abstract"       "This text describes my WFS service. This should be some counties." ## Recommended           "wfs_enable_request" "*"  # necessary           "wms_enable_request" "*"  # necessary           "wms_title"           "WMSDemo"         END     END          #this connects to a shapefile     LAYER         NAME "County"         DATA "County"              #CONNECTIONTYPE OGR         #CONNECTION "C:/devprojects/fgsv2/shp/County.shp"         STATUS ON         PROJECTION             "init=epsg:27700"         END         METADATA             "title" "County"             "wfs_title"         "County" ##REQUIRED             "wfs_srs"           "EPSG:27700" ## REQUIRED             "gml_include_items" "all" ## Optional (serves all attributes for layer)             "gml_featureid"     "db_id" ## REQUIRED             "wfs_enable_request" "*"         END              EXTENT 0 0 700000 1300000         TYPE POLYGON         CLASS             NAME "County"             STYLE                 color 100 250 250                 outlinecolor 50 50 250             END         END          PROCESSING "CLOSE_CONNECTION=DEFER"         dump true     END #end of layer                    END # end of map file      The data I'm using is English and Welsh Counties, so the projection is set to EPSG:27700. The OpenLayers map and layers bit is as follows:                               var farm_style = OpenLayers.Util.extend({}, OpenLayers.Feature.Vector.style['default']);                 farm_style.strokeWidth = 3;                 farm_style.strokeColor = "#ff0000";                 farm_style.fillOpacity = 0.5;                                 var mymap = new OpenLayers.Map({ div: "map",                     layers: [                     new OpenLayers.Layer.Vector("County", {                         isBaseLayer: true,                         strategies: [new OpenLayers.Strategy.BBOX()],                         styleMap: new OpenLayers.StyleMap({ "default": farm_style}),                          protocol: new OpenLayers.Protocol.WFS({                             url: "http://localhost/cgi-bin/mapserv.exe?map=C:/DevProjects/fgsv2/mapfile.map",                             version: "1.0.0",                             featureType: "County",                             featurePrefix:"ms",                             featureNS: "http://mapserver.gis.umn.edu/mapserver"                         })                     })                 ,                 ],                     center: [-1, 53],                     zoom: 8                 });      When I view the POST and Content through my Firefox addIns I can see POST:               <wfs:GetFeature service="WFS" version="1.0.0" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.0.0/WFS-transaction.xsd">     <wfs:Query typeName="ms:County">     <ogc:Filter>     <ogc:BBOX>     <ogc:PropertyName>the_geom</ogc:PropertyName>     <gml:Box srsName="EPSG:4326">     <gml:coordinates decimal="." cs="," ts=" ">     -2.64794921875,51.35205078125 0.64794921875,54.64794921875     </gml:coordinates>     </gml:Box>     </ogc:BBOX>     </ogc:Filter>     </wfs:Query>     </wfs:GetFeature>      The reference to EPSG:4326 is abit confusing as I've made no mention of it elsewhere. Is it a default setting somewhere? The Content tab contains the following:                   <wfs:FeatureCollection xsi:schemaLocation="http://www.opengis.net/wfs     http://schemas.opengis.net/wfs/1.0.0/WFS-basic.xsd     http://mapserver.gis.umn.edu/mapserver http://localhost/cgi-bin/mapserv.exe?map=C:/DevProjects/fgsv2/mapfile.map&SERVICE=WFS&VERSION=1.0.0&REQUEST=DescribeFeatureType&TYPENAME=ms:County&OUTPUTFORMAT=XMLSCHEMA">     <gml:boundedBy>     <gml:Box srsName="EPSG:27700">     <gml:coordinates>     294170.798523,90299.411200 655604.702242,588517.298279     </gml:coordinates>     </gml:Box>     </gml:boundedBy>     <gml:featureMember>     <ms:County fid="County.1">     <gml:boundedBy>     <gml:Box srsName="EPSG:27700">     <gml:coordinates>     320262.500400,378383.398000 361791.001700,422761.904601     </gml:coordinates>     </gml:Box>     </gml:boundedBy>     <ms:msGeometry>     <gml:MultiPolygon srsName="EPSG:27700">     <gml:polygonMember>     <gml:Polygon>     <gml:outerBoundaryIs>     <gml:LinearRing>     <gml:coordinates>     330996.596900,394516.299901 338224.102200,379434.497000 328536.798900,380475.199200     326182.199000,378383.398000 320262.500400,388522.196000 330996.596900,394516.299901     </gml:coordinates>     </gml:LinearRing>     </gml:outerBoundaryIs>     </gml:Polygon>     </gml:polygonMember>     <gml:polygonMember>     <gml:Polygon>     <gml:outerBoundaryIs>     <gml:LinearRing>     <gml:coordinates>     330871.099200,405417.000200 335362.999800,402517.699900 337285.798200,405717.803700     344062.801400,398931.696000 348231.000800,404146.903200 351662.604000,402905.498700     361791.001700,394518.799500 355128.200600,395322.399401 354161.500000,387836.602400     350551.103500,389713.904500 344678.600001,382585.295100 334583.800700,388194.599500     326895.800900,406184.799700 336963.201600,422761.904601 330871.099200,405417.000200     </gml:coordinates>     </gml:LinearRing>     </gml:outerBoundaryIs>     </gml:Polygon>     </gml:polygonMember>     </gml:MultiPolygon>     </ms:msGeometry>     <ms:OBJECTID>1    </ms:OBJECTID>     <ms:COUNTY_STR>Merseyside    </ms:COUNTY_STR>     <ms:PMS_REGION>7    </ms:PMS_REGION>     <ms:COUNTRY>ENGLAND    </ms:COUNTRY>     <ms:SNAC_GOR>North West    </ms:SNAC_GOR>     <ms:SHAPE_LEN>2.19563981295e+005    </ms:SHAPE_LEN>     <ms:EARegion>3    </ms:EARegion>     <ms:db_id>1    </ms:db_id>     <ms:updateID>0    </ms:updateID>     <ms:ET_ID>0    </ms:ET_ID>     <ms:ET_X>344814.69641    </ms:ET_X>     <ms:ET_Y>389730.08240    </ms:ET_Y>     <ms:Distance>0.00000000    </ms:Distance>     </ms:County>     </gml:featureMember>     <gml:featureMember>     <ms:County fid="County.2">     <gml:boundedBy>     <gml:Box srsName="EPSG:27700">     <gml:coordinates>     411883.903568,379684.895552 475119.498243,418312.098600     </gml:coordinates>     </gml:Box>     </gml:boundedBy>          These co-ordinates are correct for the data I'm trying to serve so I'm happy that things are appearing to be read ok. Panning and zooming in the map div on my test page will result in the content being altered based on a different bounding box, I just can't see anything. Am I missing some styling information somewhere? I've searched for solutions to this but nothing has worked. There's no proxy set up but my understanding is that for a completely local setup, a proxy is not required. Please can someone take a look at the setup I have and let me know if I've missed something obvious.


title: Azimuthal Equidistant - 40,000 Km range

text: I've been looking at making an azimuthal equidistant map. There are several web sites that offer to make one, based on my location, but they all have the same problem. They assume my antenna rotates 360 degrees - so they map the entire world, staarting from my location, and finishing at 20,000km (1/2 earth circumference). However, my antenna setup only covers 180 degrees. As such, I'm looking for a "half circle" shaped map, that ranges to 40,000 km. Does such a map exist? Is "azimuthal equidistant" simply the wrong search term, leading me astray?


title: Custom Symbols are too small in KMZ File from ArcMap

text: We are working with the local Police department on an iOS and Android app which will have an embedded Google map showing crime data from the past month. This map is generated via a .KMZ file exported from ArcGIS. We created a set of custom .EMF icons for this map and the Crime Analyst at the Police Department imported them fine, but they are extremely small in the KMZ layer that gets added to the Google map. It looks like the icons get saved to the .KMZ file as 16x16 .PNGs. When the Crime Analyst increases the size of the icons in ArcGIS, it makes them pixelated and blurry. I've also tried unarchiving the .KMZ file and inserting larger versions of the custom symbols into the folder and then zipping it up and converting it to a .KMZ file, but then the Google Map wouldn't display any of the custom symbols. Is there a simple way to retain the original size of the imported custom icons or to increase their size without pixelating them? We have the custom icons as .SVG, .PNG, .EMF, and as a custom icon font. Is it possible to use either the .SVG or the icon font as a custom marker in ArcGIS so that scaling will not be an issue? Let me know if you need more information or to see some of the files. I don't know much about the Arc suite of software and am just the app developer, but I can get information from the Crime Analyst. **EDIT** : Here are the answers to your questions.   * What tool are you using to export the kmz? - In arctoolbox, conversion tools, "layer to kml" tool   * What is the projection of your map document? - NAD_1983_StatePlane_South_Carolina_FIPS_3900_Feet Lambert_Conformal_Conic   * What is the unit of the projection? - Feet   * Imported them into what? - Select a generic symbol. Edit symbol in **Symbol Property Editor**. Symbol Type is **Picture Marker Symbol**. Then select the .EMF file.


title: How do I get ArcGIS tiles to display a single layer when multiple layers are available?

text: I am using Bing Maps v7 to display layers from an ArcGIS server. The server displays multiple layers as one, but I would like to only display a single layer. Suppose a map service has the layers: Topographic Info (0), Places Info (5), and Citations(10) I write code to get a tile like so:               var source = new Microsoft.Maps.TileSource({             uriConstructor: function( tile ) {                 return baseUrl + "/tile/" + tile.levelOfDetail                              + "/" + tile.y + "/" + tile.x;             }         });      Which would procdure a URL like this: http://server.arcgisonline.com/ArcGIS/rest/services/World_Topo_Map/MapServer/tile/12/1528/1218 Is it possible to have the generated tile only display information regarding a single layer, such as Topographic Info (0)? I tried variations of things like ?layers=0 in the query string, but nothing seems to work.


title: How to create QGIS theme map?

text: I need to make theme map using diagrams like in MapInfo, can I do it in QGIS? Or can I create a surface using attributes?


title: Azimuth and Distance, NumericalDigitize, Numerical Vertex Edit, XY tools and other plugins are not on QGIS Plugin Manager version 2.0.1

text: Will these plugins no longer be available in QGIS 2.0.1?


title: QGIS Cloud in QGIS 2.0 missing and 2.0 Crashing

text: Is QGIS CLoud still missing off QGIS 2.0. I am having a lot of issues with 2.0. on another matter I cannot load print templates without the system crashing. Are templates non compatible. When i am also loading projects into 2.0 from 1.8 some raster layer properties have reset to default values....eg greyscale back to colour.


title: QGIS 2.0 Plugins missing

text: At the risk of a series of downvotes but I have just installed QGIS 2.0 after using 1.7.4 and 1.8 previous. It looks great but there are a lot of plugins missing. I understand from the install that I may need to install them again for fear of system conflicts. However some really good pluggins are missing. The following are the main ones:   * Numerical Vertex Edit    * Easy custom Labeling    * Open Layers Plugin    * Open Streetmap    * QGIS Cloud    * Remove Empty Layers    * Spatial Query This list is not exhaustive but I am sure they are built in somewhere and can't find them. Thanks.


title: Colorify polygons in a pythonic way in Qgis

text: I want to addd a Layer and give the layer a color for the border and for the filling. With the code below i am possible to achieve a color for the border, but a did not found a solution to give a color for the filling. In addition, is it possible to change the width of the border line?               uri.setDataSource("shema_name", "layer_name", "geom","","layer_name_id")     lyr = QgsVectorLayer(uri.uri(), "layer_name_you_like", "postgres")          props = { 'color_border' : '255,165,0,125', 'style' : 'no', 'style_border' : 'solid' }     s = QgsFillSymbolV2.createSimple(props)     lyr.setRendererV2( QgsSingleSymbolRendererV2( s ) )          QgsMapLayerRegistry.instance().addMapLayer(lyr)


title: Csv To Shape Error

text: I have found a usefull tool to convert csv files into shape files with geotools: http : // docs.geotools.org/latest/userguide/tutorial/feature/csv2shp.html#create-a- shapefile-from-a-featurecollection But i have a problem running it, i get following error: Failed to execute goal org.codehaus.mojo:exec-maven-plugin:1.2.1:exec (default-cli) on project tutorial: Command execution failed. Process exited with an error: 1 (Exit value: 1) The error comes from this line:               final SimpleFeatureType TYPE = DataUtilities.createType(                     "Location",                   // <- the name for our feature type                     "location:Point:srid=4326," + // <- the geometry attribute: Point type                     "name:String"         // <- a String attribute             );      My POM file is this:               <project         xmlns="http://maven.apache.org/POM/4.0.0"         xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"         xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">         <modelVersion>4.0.0</modelVersion>         <groupId>org.geotools</groupId>         <artifactId>tutorial</artifactId>         <version>0.0.1-SNAPSHOT</version>         <packaging>jar</packaging>         <name>tutorial</name>         <url>http://maven.apache.org</url>         <properties>             <project.build.sourceEncoding>UTF-8</project.build.sourceEncoding>             <geotools.version>11-SNAPSHOT</geotools.version>         </properties>         <dependencies>             <dependency>                 <groupId>junit</groupId>                 <artifactId>junit</artifactId>                 <version>3.8.1</version>                 <scope>test</scope>             </dependency>             <dependency>                 <groupId>org.geotools</groupId>                 <artifactId>gt-shapefile</artifactId>                 <version>${geotools.version}</version>             </dependency>             <dependency>                 <groupId>org.geotools</groupId>                 <artifactId>gt-swing</artifactId>                 <version>${geotools.version}</version>             </dependency>             <dependency>                 <groupId>org.codehaus.mojo</groupId>                 <artifactId>exec-maven-plugin</artifactId>                 <version>1.2.1</version>             </dependency>         </dependencies>         <repositories>             <repository>                 <id>maven2-repository.dev.java.net</id>                 <name>Java.net repository</name>                 <url>http://download.java.net/maven/2</url>             </repository>             <repository>                 <id>osgeo</id>                 <name>Open Source Geospatial Foundation Repository</name>                 <url>http://download.osgeo.org/webdav/geotools/</url>             </repository>         <repository> <!--Add the snapshot repository here-->             <snapshots>                 <enabled>true</enabled>             </snapshots>             <id>opengeo</id>             <name>OpenGeo Maven Repository</name>             <url>http://repo.opengeo.org</url>         </repository>       </repositories>         </project>      I hope someone can help me!


title: How to register as versioned using SDELayer?

text: This is the normal process using ArcCatalog : ![enter image description here](http://i.stack.imgur.com/9Tiu4.png) I want to automate this process using sdeLayer or any tools - does anybody have an idea as to how to do it? **there is an application sdelayer (sdelayer.exe) in "C:\Program Files\ArcGIS\ArcSDE\pgexe\bin" i want to create a command line to register featurelayer**               sdelayer -o register -l <table,column> -e <entity_mask> -t <storage_type>     {[-C <row_id_column>[,{SDE|USER}[,<min_ID>]]]} [Spatial_Index]     [-E <{empty | xmin,ymin,xmax,ymax}>] [{-R <SRID> | [Spatial_Ref_Opts]}]     [-S <layer_description_str>] [-k <config_keyword>]     [-i {<service> | <port#> | <direct connection>}] [-s <server_name>]     [-D <database_name>] -u <DB_user_name> [-p <DB_user_password>] [-q]


title: Google map showing Overlapping Boundaries

text: I am working in process of building an online Geocoding and Reverse Geocoding Service .For That I have been preparing a database by creating a centroid point in a shapefile and sending the point to Google API and getting address polygons.After that i will merge similar address polygons. my problem is this: By reverse geocoding process,I took an example of Two locality polygons say : **_Amjikarai,Bala Vinayaga Nagar_** :the obtained output polygons of two area is: ![enter image description here](http://i.stack.imgur.com/vnbe9.jpg)                       where there is disjoints and Intersections .      So I thought to check the their respective boundaries in Google map.I was shocked that Google itself showing overlapping boundaries of the two areas: ![enter image description here](http://i.stack.imgur.com/H9XcV.jpg) So what is the problem or how can I rectify this problem?


title: Sum many netCDF rasters

text: For the first time, I am using NetCDF data files containing daily data about soil moisture at a global level. I want to be able to obtain monthly mean value of soil moisture for each grid cell. **There is one NetCDF file for each day, and each file contains four bands. I only need the first band which contains the data I want.** Using QGIS 1.8, I encounter the following problems: _1. Loading the rasters_ When I want to load the 31 rasters related to January (example), I can select the 31 files, but then, I have to select the band I am interested in for each layer, one by one. **Is there a way to load all of them, but automatically select the band I want?** ![enter image description here](http://i.stack.imgur.com/u9EAc.png) _2. Converting the layers_ The data comes in a floating number format. I detected some issues when it comes to sum values. However, when I convert the layer to a .tiff raster, it seems to transform the value into an integer. This is what I want, but GDAL doesn't want to batch convert all the files into .tiff, either because of the file names of because of the multiband issue. **Is there a way to batch convert NetCDF files into .tiff?** ![enter image description here](http://i.stack.imgur.com/CquXN.png) _3. Batch summing the layers_ Finally, even if I have the layer in .tiff format, I wonder if there is an easier way to sum all the values than by using RasterCalculator. **Is there a tool than batch sums up rasters within a folder?** Also, I noticed that I cannot use RasterCalculator directly on the NetCDF layers because of the file names. If the solution is on this trail, **is there a way to batch change the name of a layer within a NetCDF file?** That's a lot of questions, so many thanks in advance for your help! Damien


title: GeoServer preventing overwriting html files

text: I have GeoServer installed on Windows XP as part of the OpenGeo bundle. It's working well and I have just started creating my own app using it. I have created an html file in \WebApps\apps and can view it via localhost:8080. My problem is that once I have viewed the file in a browser it becomes locked and I can't save over it. Has anyone else had this problem? Is there a fix?


title: Convert surveyor's calls in excel table to polygon shapefile

text: Do you know of a way to take surveyor's calls from an excel file (for example: N35-45-15W 60.0 ft, etc) and have them generate polygon shapefiles point by point in GIS? I know you can easily do this with point features that have coordinates, but I am trying to figure out how to build full polygons (without XY coordinate information, just surveyor's calls).


title: Working out parcel direction

text: I am looking at applying a direction indicator to land parcels and houses (in MasterMap). effectively, I am looking at applying a direction indicator to the land and the buildings associated to the land (NW-SE, N-S, E-W etc). Is there an easy way to do this? I guess you'd need to also find the direction of the road it was situated on for completeness?


title: Separate tilestache cache for each mapnik

text: I'm using tilestache with mapnik source and redis cache. I'm wondering if there is a possibility for defining seperate cache for each layer? I want to have unique "key prefix" for each layer, therefore redis administration would be simpler (for me). What is more I want to have all tilestache definitions in one config/cfg file.


title: Watershed Segmentation SAGA-GIS

text: I'm trying to segment a topographic map to vectorize it. I'm using the watershed algorithm of SAGA-GIS. It's working fine with polygons, but with roads or contour lines is not working properly. Here is an example about I refer... the countour lines have differents seed values...   * Anyone knows how could I improve it to obtain clear lines?   * Anyone knows how to remove the square appearance after vectorization? Thank you very much! ![Contour Lines after Watershed segmentation](http://i.stack.imgur.com/hoOBm.png)


title: How to get ESRI Metadata into GeoNetwork?

text: Has anyone had success in transferring metadata created in ArcGIS 10.2 into GeoNetwork? Currently I am creating metadata with ArcGIS using the Item Description menu item in ArcCatalog. I would like to store and publish that metadata in GeoNetwork but I have not been able to get this to work. I have tried to import that xml into GeoNetwork using a variety of formats and style sheets but I either get a validation error or incomplete information.


title: What causes an Out of Memory error in ArcServer?

text: I'm creating a web app in the ESRI JS API and have recently added a 2008 FIRM service to the app. I cached the service. When I would try and click and identify certain flood zones I get an error in the console that says `ERROR: Out of Memory` I thought maybe this was due to the cache since I believe its around 7gb. So I created a new service with the data without a cache. Same Error. Is there a good way to troubleshoot this? Is there a typical reason why this error occurs? The data is only for 1 county so its not an overly large dataset by any means. I'm using ArcServer 10.1 **SERVER SPECS** 4x 2.0GHz 4th Generation Intel Xeon E5-2650 CPUs w/ 1600MHz FSB, Hyper- Threading, and Virtualization Technology enabled 8GB DDR3 1600MHz quad-channel ECC RAM VMXNET3 Ethernet Adapter connected at 10Gbs Windows Server 2008R2 – Standard Edition Virtual Hardware version 9 (vmx-09) w/ the latest VMware tools running on ESXi 5.1 Enterprise Plus **UPDATE** The memory error only happens in IE. The app works fine in Chrome and Firefox.


title: Frequency analysis and overlay analysis

text: Am new to gis and I hope to find help here. I am working on a research that requires that I overlay conditions of drainage systems and quality of healthcare in a coastal state to get areas which have potential for increase in the incidence of malaria. Each of these two attributes of areas in the state have more than one value. For instance, under quality of health care, I have data on:presence of healthcare service, drug provision and attitude of healthcare personnel. Under conditions of drainage system, data on: presence of open drainage, effectiveness of waste management. The data is derived from answers from respondents in a household survey. The data is in numeric form Am using arcgis 9.2 To my question, can I use frequency analysis in gis to summarize the data and go on to overlay the two attributes


title: How to disassociate roughness from slope

text: How to measure topographic roughness irrespective of slope in ArcGIS? I have used many methods to generate roughness layers, including TRI, VRM (which is supposed to be less correlated with slope), slope aspect index, slope variance, slope standard deviation etc. So far, all of the output layers are extremely similar to each other and to my slope layer. They do not seem to be clearly differentiating between areas of steep and rugged terrain and those which are steep but smooth. How can I generate a roughness layer which is not so dependent on slope?


title: Is it possible to create a standalone model in ArcMap Modelbuilder?

text: Is it possible to create an algorithm in ArcMap Modelbuilder that could then become a fully standalone applet? I would like to be able to provide this processing algorithm even to users who do not own an ArcMap license.


title: Esri basemap render speed abnormally slow in Chrome

text: **_Summary of Issue:_** In my project, using the oceans basemap from Esri is very slow to render when zoomed in very close (for Chrome only). Topo is also slow at high zoom levels, but it is bearable in comparison to oceans. Using Chromes developer tools to look at the Network timing for the tiles reveals some are waiting 400+ ms, but this is also the case for topo so that's not too helpful. Profiling using the Timeline tool shows that painting takes considerably longer with the oceans basemap. You can try it yourself quickly in the sandbox by replacing "topo" with "oceans": http://developers.arcgis.com/javascript/sandbox/sandbox.html?sample=map_simple **_Question:_** Is there anything on my end that I can do to optimize this in Chrome?


title: Open Raster stored in ESRI V10.0 Geodatabase in QGIS V2.4

text: Is there a way to open a raster that is stored in ESRI V10.0 Geodatabase in QGIS V2.4. I can open the vector layer, but there does not seems to be a way to open a raster.


title: Set a split policy/rule on a field without domain in ArcMap

text: Can you define a split policy or split rule on a field without a domain? When a domain is involved, you can set the default value of the field and chose a rule for when a feature is split or merged. The setSplitRule exists, but the only two options are 'None' or 'Ratio'. 'None' sets the attributes to the copy of the original value. http://resources.arcgis.com/en/help/main/10.2/index.html#//018z0000004v000000 I am looking to set a field to the Default Value (null) instead of copying the value from the original feature.


title: How to query one shape file using another shape file?

text: How to query a shapefile using another shapefile in ArcGIS? For example: I had one shapefile with states, another shapefile with towns. Now I want to query towns using the states file like towns of particular states. Note: I can extract the data using _clip_ , _union_ or _intersect_ tools in ArcGIS. Here, I don't want to extract data, I need to view the data.


title: C++ parser for OGC Filter Encoding 1.0

text: I'm developing a C++ application and I have to parse a OGC Filter in order to implement querying functionality on a metadata repository. Does anyone know any C++ library to parse OGC Filters? Thanks.


title: Can I simplify polygons in QGIS by Size?

text: Just wondering if it's possible in QGIS to simplify polygons based on size? I've done it many times using the tolerance level to reduce the overall size however in this case each polygon cannot be bigger than a certain size? Or can I simplify only those polygons which have more than a certain amount of points (which would determine size). So the very small polygons wont be affected. I'm attempting to create a google map of every polygon so I can't exceed a certain size. Thanks


title: Using gdalwarp with multiple elevation models

text: I'm trying to use `gdal_warp` to slice and convert some SRTM data. The data I'm interested in are in S30E152.hgt and S30E153.hgt. What I'd like to do is:   1. Convert these DEM files from EPSG:4326 to EPSG:3857   2. Select a subsection of the data at coordinates 17008000,-3590694 17051000,-3503623   3. Write this subsection to a GeoTIFF file for further processing Firstly, I add the two DEM files to a VRT: `gdalbuildvrt combined.vrt S30E152.hgt S30E153.hgt` Then I check the VRT to ensure the height values are sane: `gdalinfo -mm combined.vrt [snip] Band 1 Block=128x128 Type=Int16, ColorInterp=Undefined Computed Min/Max=-14.000,1501.000 [snip]` Then I use `gdal_warp` to convert and slice the VRT: `gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -te 17008000 -3590694 17051000 -3503623 -r bilinear combined.vrt sliced.tiff` Then I check the GeoTIFF to ensure the height values are still sane: `gdalinfo -mm sliced.tiff [snip] Band 1 Block=449x9 Type=Int16, ColorInterp=Gray Computed Min/Max=0.000,0.000 [snip]` It appears that whenever I use `gdalwarp` with a VRT containing multiple DEMs, it always clobbers the height values. For comparison, I created a VRT with only a single DEM and tried the same operations (using a smaller subsection that fits within range). `gdalbuildvrt test.vrt S30E152.hgt gdalinfo -mm test.vrt (produces Min/Max=-14.000,1501.000) gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -te 16931600 -3388300 16976200 -3439400 -r bilinear test.vrt test.tiff gdalinfo -mm test.tiff [snip] Band 1 Block=446x9 Type=Int16, ColorInterp=Gray Computed Min/Max=53.000,1471.000 [snip]` This seems to work OK. Maybe the problem was with my coordinates. So I thought I'd warp the entire VRT, and then select the subsection: `gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 combined.vrt converted.tiff gdalinfo -mm converted.tiff [snip] Band 1 Block=2328x1 Type=Int16, ColorInterp=Gray Computed Min/Max=-14.000,1501.000 [snip] gdaltranslate -projwin 17008000 -3503623 17051000 -3590694 converted.tiff sliced.tiff` Although I've double-checked the coordinates, I always get this warning: `Warning: Computed -srcwin 914 1337 449 910 falls partially outside raster extent. Going on however.` Now this gives me a file with reasonable height values: `gdalinfo -mm sliced.tiff [snip] Band 1 Block=449x9 Type=Int16, ColorInterp=Gray Computed Min/Max=0.000,369.000 [snip]` However when I check to see if the data is sane, I get an empty image: `gdal_translate -scale 0 369 0 255 -outsize 100 200 -of PNG sliced.tiff sliced.png` At this point, I'm totally lost. I think I have data, but am not sure.


title: "Definition Query on a versioned Geodatabase"

text: I would like to implement Specific views to the SDE users, based on Definition Queries of common attributes , like Date created > = October 1. But on first look, the definition queries are for "feature Display setting" of each featureclass - Our aim is to implement it for the Feature dataset in a particular version The below article did not help, as I could not locate , choose Production Cartography In the Categories list. http://resources.arcgis.com/en/help/main/10.1/index.html#//0103000000s5000000 What is the procedure to display features on a common attribute, across the feature dataset in a Version.


title: How to capture contour lines from OSM cycle map as vector?

text: In QGIS, there is a plugin to insert a OSM cycle map layer, that brings as an image from http. But I can't convert this OSM Plugin layer to a raster image. I can't even imagine how to get the contour lines that are showing by the cycle map as a vector. Is this possible? If is not possible through OSM, there is another open resource to get this?


title: ArcMap 9.3 Hangs on trying to save

text: I am having trouble with a machine running Arc GIS 9.3. When one tries to save an mxd, ArcMap stops responding and just hangs. I have even tried with an blank new mxd but the moment one clicks save it hangs without even opening the dialogue box that allows you to navigate to the location to save. Please help


title: Projection ETRS89 (Ellipsoid: GRS80) 6-Degree Zone 33 to UTM33N - WGS84

text: Hi I'm having trouble to reproject a shape file in Qgis. The projection is supposed to be "ETRS89 (Ellipsoid: GRS80) 6-Degree Zone 33" and I would like to project it to **UTM33N - WGS84 (EPSG 32633)**. When I set the Layer CRS to EPSG 25833 or EPSG 4258 and save it as UTM33N-WGS84 I don't get the result as intended... This is the file I would like to reproject: http://www.projektidee.org/images/extern/landschaftsschutz.zip It should fit to this area in the end: http://www.projektidee.org/images/extern/gebiet_UTM33.dbf.zip Maybe one of you could give me a hint what I'm doing wrong... Best regards desputin


title: How to transfer attributes to thousands of local points from ArcSDE feature class with millions of polygons?

text: I have a set of points and a huge feature class with topographic information (OS MasterMap for who knows it). I'm talking millions of polygons huge. I want to transfer attributes of the polygons to the points where they intersect. I've tried spatial join, but it ran out of memory. I've tried identity and it's giving me no clue at all about how long it might take, though I'm going to let it run over the weekend to see where it ends. I then thought of making a script to take one point at a time and use that to select the underlying polygons. Doing a select layer by location with the little dialog window in ArcMap works great, it takes less than a second. However doing just that from the python window takes many minutes. I've started writing a full ArcPy script to do a select by location one by one point and then update the point's attributes with the selected polygons' attributes. Updating one point took about 15 minutes. And I've got about 8000 points to do. Unfortunately my boss won't let me sit idle and wait for a script to finish if it takes 80 days to run :( Anyone has an idea about how this might be done more quickly, and why (why ESRI, why???) SelectByLocation works so much faster in ArcMap itself? edit: maybe should mention the points feature class is local and the polygons feature class is SDE ============================================================================= SOLUTION: Among many suggestions, the only feasible one for me, being unable to install anything or use multiple workstations at the same time, was Polygeo's solution, please see the answer below. The key thing to do is: set a processing extent to a very small area around each point (do this in a cursor), make a one-polygon copy of the large dataset using this processing extent and get the attributes from the small copy. I'm bound to try to do a query within the oracle database itself as well, but I'll need some assistance from a colleague with that. As well I'm trying to get PostGIS installed and use that. I'll report back on these as well as soon as I have results. Thanks for your help!


title: Count points inside polygons using PostGIS

text: I have 2 vector data in my Postgis DB: one is a point layer (punti) with about 200000 elements and the other is a polygon (province) with 4 elements. They cover the same area and I need to find the number of points for each polygon. This is the query:               province.gid, count(*) AS totale FROM punti, province WHERE      st_intersects(province.the_geom,punti.the_geom) GROUP BY province.gid;      This is the EXPLAIN ANALYZE > HashAggregate (cost=35.35..35.36 rows=1 width=4) (actual > time=23660.787..23660.789 rows=4 loops=1) -> Nested Loop (cost=0.00..35.34 > rows=1 width=4) (actual time=0.307..23532.479 rows=199977 loops=1) Join > Filter: _st_intersects(province.the_geom, punti.the_geom) -> Seq Scan on > province (cost=0.00..1.04 rows=4 width=36) (actual time=0.003..0.009 rows=4 > loops=1) -> Index Scan using punti_the_geom_idx on punti (cost=0.00..8.31 > rows=1 width=100) (actual time=0.038..76.254 rows=69541 loops=4) Index Cond: > (province.the_geom && punti.the_geom) Total runtime: 23660.850 ms What do you think about this total runtime? Thx!


title: Looking to combine temporal data from multiple shapefiles to produce a days per year file for each feature

text: I have many shapefiles in which there is a start date and an end-date per feature. I am hoping to produce a single shapefile that has each enclosed area with a "days per year" kind of calculation. Trying to visualize this... if feature a is [ __ ] and feature b is ( __ ) and they overlap... **[ feature a 10/1/2011-10/3/2011 _( overlap ]** feature b 10/2/2011-10/8/2011 )_ into **[ feature a = 3 days ]** [ feature c= 8 days ] _[ feature b = 6 days ]_ So, overlaps calculate properly and become their own feature. Not sure if this is making sense? Any help would be great! I've tried using unions, looked at rasters, checked out the toolboxes... Not sure what route to take. My final map I'm hoping to produce should be a map that shows the duration differences. Thanks!


title: Extract Data widget flex viewer 2.5

text: I am working on the data extract widget for the flex 2.5 viewer. It has installed properly and works fine but the output is not what I expected. The tool will not divide or split polygons. For example when I draw the AOI the return will omit any polygons that have been split or divided. I was expecting to get all data in the AOI. Is this the way the tool works? Or is there something I can do to fix this? My .xml file for the widget is below: -  --> - http://s2-cor- pangaea:8399/arcgis/rest/services/demo1/GPServer/Extract%20Data%20Task draw


title: ArcGIS 10 Raster Mosaicing versus Lizardtech Express Server

text: Trying to determine if I still need to use Lizardtech Express Server to mosaic rasters together or just use the Mosaic Raster Dataset function of ArcGIS 10? I know for earlier releases of Express Server and ArcGIS you could find info on tying the 2 together but I can't seem to find anything about it now. I've learned that when the silence is deafening on a topic it generally means that "they" want it to just die and go away. Help me StackExchange GIS, you're my only hope.


title: Best way to build my own simple web-based map server

text: I've got a requirement for a project where I need some relatively simple mapping functionality on a web page. I'm familiar with the Google Maps APIs but unfortunately this is a project that will be hosted on a closed network with no internet access. Is there a relatively straight-forward way to host my own map server with the following functionality (preferably on Windows but I would consider Linux):   1. I only need to cover a map area no bigger than the size of Texas   2. I don't need streets, terrain, or any other complex layers. At most I would needs outlines of the counties in the region with a few text markers for large cities.   3. I need to be able to display custom data driven markers on the map when it is rendered (I don't care if that is on the server side or if I need to do it in javascript on the client)   4. 2 or 3 levels of zooming would be great, but no necessary I've tried SUAS server v3 with XAMPP and ran into numerous issues, including low-level PHP errors.


title: Finding a route with desired cost not minimum cost

text: The goal of routing is usually finding the route with minimum cost. That cost can be distance, time, calories or some other value but still our desire is to find the route that costs least. But what if I want to find a route with desired cost? Not the fastest route but the route that will take around an hour (let's assume I have an hour to walk from point A to B and I want to see as much on my way there as I can). Or that I want a route that will be around 2 kilometres long. Using a standard routing algorithm and rejecting the proposed routes until one with desired cost is found will probably lead to something like "Head straight from A to B and then circle around B a few times". Is there any previous work on this? Does anyone have any ideas?


title: How to smooth the polygons in a contour map

text: Here is a contour map for which all the polygons of levels are available. Let ask **how to smooth the polygons keeping all vertices preserved in their exact locations?** Indeed the contour is made on top of a grid data, you may suggest then to smooth the grid data and hence the resulting contour will be smoother. Note that this is not working as my desire since the smoothing function such as Gaussian filter will remove small packs of data and will change the range of the third variable e.g., the height which are not allowed in my application. Actually I'm looking for a piece of code (preferably in **Python** ) which can do the smoothing of 2D polygons (any type: convex, concave, self-intersecting etc) reasonably painless (forget pages of codes) and accurate. FYI, there is a function in ArcGIS that does this perfectly, but using third party commercial applications is not my choice for this question. ![enter image description here](http://i.stack.imgur.com/nAeyI.png) **Updates**   **1)** Scipy.interpolate: ![enter image description here](http://i.stack.imgur.com/ZLMhT.jpg) As you see the resulting splines (red) are not satisfactory! **2)**   Here is the result using the code given in here. It is not working well!   ![enter image description here](http://i.stack.imgur.com/JL9aP.png) **3)**   To me the best solution should be something like the following figure in which a square is being smoothed gradually by changing only one value. I hope for similar concept for smoothing any form of polygons. ![enter image description here](http://i.stack.imgur.com/xHyxR.png) Satisfying the condition that spline passes the points:   ![enter image description here](http://i.stack.imgur.com/5cmVu.png) **4)**   Here is my implementation of "whuber's idea" line by line in **Python** on his data. There are possibly some bugs since the results are not good. ![enter image description here](http://i.stack.imgur.com/b0E2b.png)   K=2 is a disaster and so for k>=4. **5)**   I removed one point in the problematic location and the resulting spline is now identical to whuber's. But it is still a question that why the method does not work for all cases?   ![enter image description here](http://i.stack.imgur.com/3isPI.png) **6)**   A good smoothing for whuber's data can be as follows (drawn by vector graphics software) in which an extra point has been smoothly added (compare with update 4):   ![enter image description here](http://i.stack.imgur.com/4cVxi.png) **7)**   See the result from Python version of whuber's code for some iconic shapes:   ![enter image description here](http://i.stack.imgur.com/803BS.png)   **Note** that the method seems doesn't work for polylines. For the corner polyline (contour) green is what I want but got red. This needs to be addressed since contour maps are always polylines although closed polylines can be treated as polygons as in my examples. Also not that the problem emerged in update 4 has not been addressed yet. **8) [my last]**   Here is the final solution (not perfect!):   ![enter image description here](http://i.stack.imgur.com/YO8TF.png)   Remember that you will have to do something about the area pointed by stars.There is perhaps a bug in my code or the proposed method needs further development to consider all situations and to provide desired outputs. Thanks and Good luck!


title: Looking for market share of GIS desktop tools

text: What are the relative market shares of the various desktop GISes? I'm interested in this within the open-source and non-opensource realm. One article that I read, is that the market share of ESRI is estimated to be approximately 30% of the global GIS market in 2009 (Daratech, 2011 - dead link unfortunately) but that they have a much higher perceived 'mind-share' (Batty). I started thinking about this as I wondered which was the most popular open- source GIS program, and assumed it was QGIS but I'm finding it hard to locate any numbers. To clarify (in response to a comment), I am interested in GIS Desktop tools that allow quantitative analysis (which would exclude Google Earth).


title: OpenStreetMap or Google Map Maker: where should I donate and volunteer my geographical knowledge?

text: I definitely know a lot about the places I visit and want to share it with everyone. I want to know if donating and volunteering at OpenStreetMap—which is where I _have_ been helping—is a waste of time: it is possible that it is all for naught to help on such a relatively low-trafficked site. On the other hand, I am concerned about the closed, commercial nature of Google Maps, even though the traffic is high and my knowledge would likely get more reuse. Where should I help with collaborative mapping?


title: NDVI from QGIS is giving me NaN values

text: NDVI,TTVI maps are nan values. when I do the NDVI calculation using the raster calculator it looks fine. what's going on?


title: qgis-web-client stuck on Loading Map

text: I'm trying to publish a map created in qgis desktop via qgis-web-client. The left hand bar loads with the layers I would expect visible but I never actually see any mapping data just a grey screen with "Loading Map....". If I make the call to the qgis server manually via a web-browser e.g.               http://adddress/qgis/qgis_mapserv.fcgi?map=/var/www/qgis-web-client/projects/test/filename.qgs&&_dc=1394809120649&SERVICE=WMS&VERSION=1.3&REQUEST=GetCapabilities      I get what looks like valid XML returned. In the chrome dev/javascript console I'm seeing error: ` Uncaught TypeError: Cannot read property 'abstract' of undefined ` From webgisinit.js line 1680 ` var thisAbstract = wmsLoader.projectSettings.service.abstract; ` Anyone seen this before or got any clues how to fix?


title: Working with GDB in QGIS 2.0 - plain language instructions please!

text: I'm a QGIS user running v2.0 (Dufour) on Mac OSx10.8.5. I have some data in .gdb format and need to access this and work with it in QGIS. I'm a relatively savvy computer user, but not at all a programmer or developer. I've been reading the forums and have seen some steps outlined for GDAL installs etc., and compiling etc., but am a bit daunted and don't want to mess something up on my system. Can some patient soul please spell out for me the steps I need to take to get .gdb files functioning in QGIS? Alternately I could try to use GDB Flee or something like it, but not sure where to start. Thank you so much. I'm learning so much from reading these forums but still have a long ways to go!


title: One-To-Many Table Relate Viewed Through ArcGIS Server

text: I want to create a Table Relate with a one-to-many relationship between a layer and a non-spatial table based on a common field. My question is, if I publish this data to ArcGIS Server, will an identify operation on that layer return _all_ of the related values in the non-spatial table? I want to try this, but I don't have sample data that meets these requirements (if anyone knows where I can find some, it would be very helpful). Here is an example from ESRI of an Identify operation in ArcMap on two tables with a one-to-many relationship. I need to know if the results in ArcGIS Server would be consistant. I am using ArcGIS 10 and SQL Server 2008 R2. Thanks, Seth


title: Monitoring programmes

text: I am not sure if it is a right place to ask this one but let me try. Monitoring programmes are known for various targets and for various spatial levels. They may include regions, countries or even continents. It is widely accepted that monitoring systems are of increasing importance as economic activities and their impacts on nature become more and more global. I am looking for a various monitoring programmes. The best will be that related to forestry (but in fact i am seeking for anyone) like this one I am interested how do you evaluate the selected level I and level II indicators used in the ICP forest ecosystem monitoring programme? Do you consider them exhaustive to get a comprehensive view of the state of the forest ecosystems under the influence of air pollution? Also if you have some examples in mind it will be great to share it!


title: Find AB_CSRS.gsb for NAD_1983_TO_WGS_1984_8 Transformation

text: I am trying to change my WGS84 data to NAD83 for Alberta data. My research tells me that I need to use the NAD_1983_TO_WGS_1984_8 Datum Transformation which is in the AB_CSRS.gsb. I can not find this gsb file for download on the internet. I did manage to find and obtain the NTv2_0.GSB but that seems to only be for NAD27 To NAD83. Does anyone know where I can get the AB_CSRS.GSB file? Thanks


title: Split Lines at Intersection of Other Lines

text: I have two line layers (roads and contours in this example) and I would like to split the roads layer every time it crosses/intersects a feature on the contour layer. I can't seem to find a tool that does this, as the intersect tool is expecting a polygon for at least one of the layers. It also appears that GRASS' v.clean only works on a single layer. It seems like a simple task, but I'm at a loss!


title: How to make a different font and color of text in labels?

text: Is it possible to make a bold or italic, and change the color of individual words Label? Example: the first line is bold in red text, font 12, second row 10 green italic font ![enter image description here](http://i.stack.imgur.com/0y3Gr.png)


title: Misalignment when generating centroids with Qgis

text: I was wondering if anyone could shed a light on why this behavior is happening in QGIS ? ![alt text](http://i.stack.imgur.com/mziOb.jpg) I'm pretty confident that the issue here is not due to "centroids within polygon" **Edit** Added a detailed image ![alt text](http://i.stack.imgur.com/JndXW.jpg)


title: Issues using Arc 10.0 extension recompiled for 10.1

text: I have a written an ArcGIS extension in C# .Net for version 10.0. I have attempted to recompile this extension for 10.1 and that works. But running the extension in 10.1 is causing ArcGIS to crash. While debugging the code I'm getting an error saying "ExecutionEngineException was unhandled". This is not always happening on the same part of the code. One time its a cursor, another time its setting symbology, etc. its different every time. I have read through the documentation here on how to upgrade and none of these solutions are helping. Has anyone ever experienced this before? I have no idea what else to try.


title: Network Analyst hierarchy routing issues. (ArcMap)

text: My friend and I have created a walking network. We created a network without a hierarchy evaluation and the routing worked fine, but when we added a hierarchy to the equation many routes won't solve. We created a hierarchy class and attributed each of our walk types with a hierarchy of 1, 2, or 3. We thought that the cause of the routes not working was when a route would jump from a hierarchy of 3 to 1, or 1 to 3, skipping 2. But routes do work in some of those situations. Does anyone know why it breaks, or if there is a way to fix it? (We know that all of our paths are snapped together)


title: Query to select MIN and MAX values by groups

text: In ArcMap 10 i'm working in a Personal Geodatabase where the points (see screenshot below) represent the number of each house in the street. Each point has an ObjectID. This ObjectID is an unique value for each "link" of the street to which all the points are related.![enter image description here](http://i.stack.imgur.com/swdTz.jpg) I have a total of 35,000 points (house numbers), which are all related to 14,000 different links. What I want is to select both the MAX and MIN value in each side of the street. Each side of the street also has a unique ID, which appears in my attribute table as CVEFT. On every link, the CVEFT (for both sides) is different. This is what my Table looks like: ![enter image description here](http://i.stack.imgur.com/pBK0H.jpg) As you can see in the screenshot above, there are several points (NUMEXT) related to the same link (OBJECTID). You can also see that all the points related to a single link may have 2 different street-side codes (CVEFT). I want to build a SQL query that allows me to select the MIN and MAX values for every Link, for each side of the street. Does anyone have any idea of how to do this? Any ideas would be appreciated! Thanks!


title: Longitudinal and Latitudinal coordinates of interstate exits

text: I need to find the coordinates of all the interstate exits for some research. I'm not looking to code up a program that will do it for me, but I was hoping that someone had already done it and would be willing to share that information with me. Any help would be much appreciated.


title: Should Python Add-In installed in map which is packaged be part of *.mpk?

text: Using ArcGIS for Desktop 10.1 SP1, I just created a map package (.mpk) from a map (.mxd) which had a Python Add-In (.esriaddin) installed. When I unpacked the .mpk the .esriaddin was neither installed nor present in the folder that was created in C:\Users\Graeme\Documents\ArcGIS\Packages. Intuitively, I would expect it to be at least present in the package, so am wondering if it that capability has been omitted by accident or design?


title: R, GeoJSON and Leaflet

text: I recently learned about leafletjs.com from an R-Bloggers.com post. One such tutorial that I would like to achieve is to create interactive choropleth maps with leaflet. I have been using the rjson package for R to create the data.js file to be read by leaflet. Although I have had success with using the provided shape file as a readable JSON file in leaflet, I am unable to repeat the process when trying to merge additional properties from the data frame ("data.csv") to the JSON file. In this case, I have done rGIS to attach data on the number of cans in each school listed in the data frame. What I would please like to achieve is to create a choropleth map in leaflet that displayes high school district (as identified by the NAME variable), and the sum of "cans". The issue, I believe, is that writeOGR exports the information as points, rather than polygon?               {         "type": "Feature",         "properties": {             "name": "Alabama",             "density": 94.65         },         "geometry": ...         ...     }              ###load R scripts from dropbox         dropbox.eval <- function(x, noeval=F) {         require(RCurl)         intext <- getURL(paste0("https://dl.dropboxusercontent.com/",x), ssl.verifypeer = FALSE)         intext <- gsub("\r","", intext)         if (!noeval) eval(parse(text = intext), envir= .GlobalEnv)         return(intext)         }              ##pull scripts from dropbox          dropbox.eval("s/wgb3vtd9qfc9br9/pkg.load.r")             dropbox.eval("s/tf4ni48hf6oh2ou/dropbox.r")              ##load packages         pkg.load(c(ggplot2,plyr,gdata,sp,maptools,rgdal,reshape2,rjson))              ###setup data frames         dl_from_dropbox("data.csv","dx3qrcexmi9kagx")         data<-read.csv(file='data.csv',header=TRUE)              ###prepare GIS shape and data for plotting         dropbox.eval("s/y2jsx3dditjucxu/dlshape.r")              temp <- tempfile()         dlshape(shploc="http://files.hawaii.gov/dbedt/op/gis/data/highdist_n83.shp.zip", temp)         shape<- readOGR(".","highdist_n83") #HDOE high school districts           shape@proj4string               shape2<- spTransform(shape, CRS("+proj=longlat +datum=NAD83"))                   data.2<-ddply(data, .(year, schoolcode, longitude, latitude,NAME,HDist,SDist), summarise,                 total = sum(total),                 cans= sum(cans))                   ###merging back shape properties and data frame         coordinates(data.2) <-~longitude + latitude         shape2@data$id <- rownames(shape2@data)              sh.df <- as.data.frame(shape2)              sh.fort <- fortify(shape2 , region = "id" )              sh.line<- join(sh.fort, sh.df , by = "id" )                        mapdf <- merge( sh.line , data.2 , by.x= "NAME", by.y="NAME" , all=TRUE)              mapdf <- mapdf[ order( mapdf$order ) , ]                   ###exporting merged data frame as JSON     mapdf.sp <- mapdf     coordinates(mapdf.sp) <- c("long", "lat")     writeOGR(mapdf.sp, "hssra.geojson","mapdf", driver = "GeoJSON")      However, it appears that my features are repeating itself constantly. How can I aggregate the features information so that it looks more like the following:                 var statesData = {"type":"FeatureCollection","features":[     {"type":"Feature","id":"01","properties":{"name":"Alabama","density":94.65},"geometry":    {"type":"Polygon","coordinates":[[[-87.359296,35.00118],[-85.606675,34.984749],  [-85.431413,34.124869],[-85.184951,32.859696],[-85.069935,32.580372],[-84.960397,32.421541],[-85.004212,32.322956],[-84.889196,32.262709],[-85.058981,32.13674],[-85.053504,32.01077],[-85.141136,31.840985],[-85.042551,31.539753],[-85.113751,31.27686],[-85.004212,31.003013],[-85.497137,30.997536],[-87.600282,30.997536],[-87.633143,30.86609],[-87.408589,30.674397],[-87.446927,30.510088],[-87.37025,30.427934],[-87.518128,30.280057],[-87.655051,30.247195],[-87.90699,30.411504],[-87.934375,30.657966],[-88.011052,30.685351],[-88.10416,30.499135],[-88.137022,30.318396],[-88.394438,30.367688],[-88.471115,31.895754],[-88.241084,33.796253],[-88.098683,34.891641],[-88.202745,34.995703],[-87.359296,35.00118]]]}},          {"type":"Feature","id":"02","properties":{"name":"Alaska","density":1.264},"geometry":   {"type":"MultiPolygon","coordinates":[[[[-131.602021,55.117982],[-131.569159,55.28229],[-131.355558,55.183705],[-131.38842,55.01392],[-131.645836,55.035827],[-131.602021,55.117982]]],[[[-131.832052,55.42469],[-131.645836,55.304197],[-131.749898,55.128935],[-131.832052,55.189182],[-131.832052,55.42469]]],[[[-132.976733,56.437924],[-132.735747,56.459832],[-132.631685,56.421493],[-132.664547,56.273616],[-132.878148,56.240754],[-133.069841,56.333862],[-132.976733,56.437924]]],[[[-133.595627,56.350293],


title: How to select multiple polygons or points from WFS or WMS and getFeatureInfo

text: Is there a way to select multiple polygons and/or points from a WFS or WMS service and then get the feature information for each selected item? There are numerous examples on the OpenLayers website but none of them show how to multi-select AND get the attributes for the selected items. For instance I can get the feature info for a single selected item by using the getFeatureInfo option for a WMS service but this does not allow me to select multiple items. I absolutely have to have the ability to select multiple items! Also, I can select multiple items by drawing a box, a freehand polygon, or multi-selecting individual features using a number of different methods but none of those methods have the ability to get the feature info. What good is being able to select multiple items but not be able to get the attributes for the features? Am I missing something really obvious? My goal is to be able to select multiple features and then extract my unique identifier (point number or polygon number) from each of the features info in order to query data from a SQLServer database. I don't need to do anything fancy to the spatial data all I want to do is select multiple items and take the unique id from each of them so I can query my database for info that relates to the selected items. I am using OPENLAYERS to build my map interface and GEOSERVER to serve my WMS/WFS services. I can connect to my point and polygon shapefiles using either WMS or WFS but I can not find a way to select multiple items and then get the attribute info for the items. Is there a way in OpenLayers to select multiple items and then loop through the list of items? I have been stuck on this problem for four very long and frustrating days and if anyone has any ideas on how to do this I would be eternally grateful!! 10/11/2012 - HERE IS THE CODE I ENDED UP USING //THIS FUNCTION UPDATES THE TEXT BOX WITH THE LIST OF SELECTED ATTRIBUES WHEN THE USER SELECTS A NEW POLYGON OR MULTIPLE POLYGONS var update = function(feature) { names.push(feature.attributes['yourLayersAttributeNameHere']); showStatus(names.join(", ")); };               //THIS FUNCTION UPDATES THE TEXT BOX WITH THE LIST OF SELECTED ATTRIBUTES WHEN THE USER UNSELECTS A POLYGON     var clear = function(feature) {       name=feature.attributes['yourLayersAttributeNameHere'];          for(var y in names){            if (name === names[y]){              names.splice(y,1);             }         }         showStatus(names.join(", "));     };                         //INSTANTIATE A NEW CONTROL THAT ALLOWS THE USER TO SELECT FROM THE AVAILABLE WFS FEATURE      // var control = new OpenLayers.Control.SelectFeature( nameOfyourLayerHere, /{ clickout: true, toggle: true, multiple: true, hover: false, //if hover is set to true it takes away the selection functionality.   toggleKey: "ctrlKey", // ctrl key removes from selection /multipleKey: "shiftKey", // shift key adds to selection /onSelect: update, //ON SELECT CALL THE UPDATE FUNCTION onUnselect: clear, //ON UNSELECT CALL THE CLEAR FUNCTION // // } // );               //ADD THE CONTROL TO THE MAP CONTROLS COLLECTION      // map.addControl(control); //ACTIVATE THE NEW CONTROL //control.activate();


title: Select only features within the layout view

text: I want to create a tool layer in ArcGIS10 that will select only features within the current layout view so that I can add the attribute table to the layout and only show attributes for the selected features. Has anyone got any ideas on how to do this? I want the selection to be done automatically each time the mxd is opened, and not have to do it manually.


title: What is the Difference between Geometric and Geographic columns?

text: I am a beginner in GIS and Learning PostGIS. I came across this two terms Geom and Geog, Could anyone explain in brief. What is the best to use with map data and Do all functions support both Geom and geog.


title: Qgis; How can I add only mean values instead of coun, sum and mean from QgsZonalStatistics (zonal statistic)?

text: I am using Qgis 2.2.0 Valmiers. I have ASCII grid, total 864 and I have a polygon layer as a zone. after zonal static batch run, attribute fields stop at 2046, fields number shown from "Table Manage" plug-in. Question 1. Is there any way I can create only mean column field instead of default three fields? Question2: Does QGis have a limit on number of column fields which I don't think? Here is a my code and it runs in a Python console. from qgis import core from qgis.utils import iface from qgis.analysis import QgsZonalStatistics import os,sys,time print "Process started: ", time.ctime() s_tm = time.time() f_list = [] LOCTN = "C:\1SFO\Tye\Data\2011\ASC_QGIS\" polygonLayer = QgsVectorLayer(u'C:/1SFO/SHP/WS_SFO_epsg4326_2011.shp','WS EPSG4326_2011','ogr') QgsMapLayerRegistry.instance().addMapLayer(polygonLayer) iter = polygonLayer.getFeatures() for files in os.listdir(LOCTN): extension = os.path.splitext(files)[1][1:] if extension == "asc": f_list.append(files) for y in range(len(f_list)):                 FLNAME = f_list[y]            INFILE = LOCTN+f_list[y]            RSTRNAME = FLNAME.replace(".asc","")[29:36]            rasterFile=LOCTN+ f_list[y]            zoneStat = QgsZonalStatistics (polygonLayer, rasterFile, RSTRNAME, 1)            zoneStat.calculateStatistics(None)      print "Process Ended: ", time.ctime() Thanks a lot


title: error using the QGIS "terrain profile" plug in

text: I have searched the home page for the "terrain profile" tool but have not found why I am getting error when trying to utilize it. Any ideas? I have a polyline and 4 dems ready to run the tool.![enter image description here](http://i.stack.imgur.com/u40u2.png) Would like to export the profile as an image and excel file.


I am using the WGS 84/UTM coordinate sysytem. When I go to > Project Properties / General / Canvas Units to Degrees, Minutes, Seconds my coordinates are wrong. I am using QGIS 2.0.1 ![enter image description here](http://i.stack.imgur.com/fQetf.png)


title: Feature class to Feature class not working

text: I am trying to use the Featureclass to Featureclass conversion tool to make a copy of an existing layer (in a geodatabase). My primary motive is to copy over only a few (2 out of 18) fields to the new feature class. I have put a simple script to do this, but this copies all the fields even tough I specify the field names that should be copied. Can someone please let me know what am i doing wrong. I want the output feature class (Test) to have two fields A and B.               # Import arcpy module     import arcpy               # Local variables:     Input = "C:\\Work\\Projects\\0a_IL_Naperville\\Segmenter_v8\\1 00001 IL Naperville 2013.gdb\\Streets"     Out_location = "C:\\Work\\Projects\\0a_IL_Naperville\\Segmenter_v8\\1 00001 IL Naperville 2013.gdb"          keep = ['ST_NAME','FUNC_CLASS','Shape_Length','OBJECTID','Shape','A']          arcpy.FeatureClassToFeatureClass_conversion(Input,Out_location, "Test",field_mapping=keep)


title: How to reclassify raster data by their sign?

text: I want to do the following: I have two data raster sets with different units and values. One results from measured metric data, the other from modelled grid values. What I want to compare is simply whether the data sets differ in where they show positive and negative values. So - all negative data on both sets should be "-1", all data between -0.001 and 0.001 should be "0", and all positive values simply "1". Do you have an idea how to do this best on ArcGIS 10? I suspect raster calculator...? Any hints are appreciated.


title: Is it possible to grow QGIS symbol size over time?

text: I've got a data set of people movements within my municipality. (between points of places) I've managed it with the plugin for QGIS 'time manager'. Over time people come and go, and based on a count of movements of a specific point at a specific time (time manager), i want to let my point symbol grow. Is this possible with QGIS rule based symbol?


title: Changes made using wfs-t and GeoExt don't get saved to PostGIS

text: I'm using GeoExt to display and edit a vector layer from GeoServer (running with Apache Tomcat) but if I modify the shape of an existing polygon (or create a new one but I'm trying not to leap ahead of myself) the change doesn't get written back to PostGIS through GeoServer. I'm using code obtained from a tutorial at http://workshops.opengeo.org/geoext/wfs/wfst.html. I see a few people have had problems with the same code but can't find any solutions that solve my problem. The code for the save button action:                   handler: function() {       app.featureGrid.store.protocol.commit(        vectorLayer.features, {         callback: function() {             var layers = app.mapPanel.map.layers;             for (var i=layers.length-1; i>=0; --i) {               layers[i].redraw(true);             }             app.featureGrid.store.reload();         }       });     }      My map layers get redrawn and the polygon changes persist but if I reload the page, the original shape is back. There is no error in Chrome developer tools but if i look in the GeoServer log file, I get:               2012-05-10 14:51:31,311 INFO [geoserver.wfs] -      Request: getServiceInfo     2012-05-10 14:51:31,321 ERROR [geoserver.ows] -      org.geoserver.wfs.WFSTransactionException: Feature type 'dev_apps' is not available:          at org.geoserver.wfs.Transaction.execute(Transaction.java:222)         at org.geoserver.wfs.Transaction.transaction(Transaction.java:108)         at org.geoserver.wfs.DefaultWebFeatureService.transaction        (DefaultWebFeatureService.java:162)      I'd really appreciate some help as I don't know where to look anymore!


title: Import and Plot MS Excel file in ArcGIS (v9.3) ModelBuilder

text: I am trying to create a tool using ArcGIS ModelBuilder (v9.3) that will allow me to import an Excel file, and output it as a shapefile in a coordinate system of my choice - or, if choosing the CRS is not an option, in WGS84. I have set the input dataset (as part of the model), but can't seem to find a tool to facilitate the import. I have an ArcView licence.


title: How to migrate SVG-based styles for QGIS to other computers?

text: I'm working on a project that has involved creating many recolourable SVG icons and bringing them into QGIS. These are then used in layer symbologies, with these subsequent layer styles exported as .qmls and .slds. When I open these .qmls on another computer the svg points are not attached - rather there is just a blank spot. I've tried to get just the svgs to load by exporting the style library from the computer that all of the svgs were created on. However, when I imported this into another computers style library everything was once again blank. Putting the folder with all of the svgs in the C:\Program Files (x86)\Quantum GIS Lisboa\apps\qgis\svg directory is not working either...they show up in the option to select svg icon, however when a layers qml file is loaded it is once again blank. I'm trying to get my head around this, so if there's an obvious answer please let me know. How can I move these svg/style to other computers so that when I load a file with a set qml layer the proper svgs load as well? Thanks anyone for your help.


title: Querying intersecting areas areas - Oracle Locator or Spatial?

text: My spatial database is Oracle. We are wondering if we can use basic Oracle Locator or need to get the more expensive, full, Oracle Spatial. **My application's requirements include having to be able to query "which geometries in the spatial db intersect this geometry X".** The whole locator / spatial issue is a bit confusing since it looks like Oracle initially installs all Oracle Spatial features even if you are just using basic Oracle db and do not have license for full Oracle Spatial. And you may end up using the Spatial features by mistake: http://www.spatialdbadvisor.com/oracle_spatial_tips_tricks/55/oracle-locator- vs-enterprise-licensing-the-sdogeom-package What I understand of Oracle documentation http://docs.oracle.com/cd/E11882_01/appdev.112/e11830/sdo_locator.htm#SPATL340 is that geometry function SDO_GEOM.SDO_INTERSECTION is **not** supported by Oracle Locator. So does this mean that we really need full Oracle Spatial license to fulfill the requirements? Or is there some other way then "SDO_GEOM.SDO_INTERSECTION" of finding about whether areas intersect? We are using GeoServer on top of Oracle (and accessing the geometries with WFS-T). Does this change anything, will the intersection be somehow calculated in GeoServer (doubt it?) or will the intesection-calculation have to be done in Oracle (probably?)?


title: identify the fraction of raster cells and listing all the cell values in every polygon

text: In Arcgis 10 I have a netcdf file that maps the fractions of cropland in the world from 1700 to 2007. I extracted the years that i needed as raster layer format (continous data from 0 to 1) and now, using a shapefile of the world's countries, I would like to identify every raster's cell inside every country (possibly in a table format where I find id of the cell, country name, fraction cultivated). Is it even possible to do it? Moreover I would like to have also the fraction of cell belonging to a country and the fraction belonging to another one (i.e. in the boundary). Please I need your help, I am new in Arcgis and Gis in general!


title: How to view values from stacked raster maps of a particular pixel at the same time (ArcGIS)

text: Re: ArcGIS v10.1 Hi. I have a stack of raster maps. I would like to click on a pixel, and view ALL of the values for that particular pixel at the same time. The only tool that I could find is the pixel inspector, which only displays the value for one map at a time. I want to compare all values from my layered maps for this particular location at the same time. Has anyone done this with ESRI software? Cheers


title: How to add graphs - pie charts, bar diagrams in geoexplorer?

text: Does any one know how to add interactive graphs in Geoexplorer as shown in this video ? Do we have to incorporate Google fusion tables ? Any kind of help will be appreciated. https://www.youtube.com/watch?v=d9ADUHGUXL0


title: How to make elements size fixed (Arcobjects)?

text: I'm using arcobjects .Net 10 to add graphic elements to my map. My problem is that when the map scale changes, the element size change too. How to fix element size at all scales?


title: Save GPX from drawn feature in OpenLayers

text:
I am trying to modify the Vector Formats example to only display the GPX in the textarea rather than all the options in the original example but I can't get it working. My code below currently displays [object Object] is the text area. The next stage would be to save it as a GPX file but I can't even get this part to work. Help appreciated.

Draw route to update the output below.

<textarea id="output"></textarea>


title: How to set custom fullextent?

text: I have a WPF Map control whose extent I would like to restrict to a certain area (e.g. to the full extent of some layer). That is, a user should not be able to zoom out too far, or shift the map to an area outside that full extent. How can I restrict Map display to a certain extent? This is my code:               <esri:Map x:Name="Map" SizeChanged="Map_SizeChanged">         <esri:ArcGISTiledMapServiceLayer x:ID="LB" ... />         ...     </esri:Map>      And here is the code-behind (this code does not work):               void Map_SizeChanged(object sender, SizeChangedEventArgs e)     {         if (currentExtent > Map.Zoom(Map.Layers["LB"].FullExtent)) // check the current extent of the map. If currentExtent is more than FullExtent         {             // jump back to full extent             Map.Zoom(Map.Layers["LB"].FullExtent);         }     }      I'd like to get something like that:               private void Map_SizeChanged(object sender, SizeChangedEventArgs e)     {         … // get current extent of the map              … // check the current extent of the map.           // if currentExtent is more than FullExtent jump back to full extent     }      I'm using the ESRI API for WPF and C# in Visual Studio 2010.


title: openlayers cluster strategy - exhibiting custom icons for clusters of 1 feature

text: i have a page that displays a **vector layer** imported from a **kml file** (basically like the the sundials example). it works perfectly well, both using fixed strategy ONLY, as well as using both fixed and CLUSTER strategies simultaneously. i would like to create a **hybrid display** , whereas **"clusters of 1 feature"** would be displayed with the original **custom icon** contained in the KML file under `<Style><IconStyle><Icon><href>img/arrowRed.png`, for example. right now, **if i use cluster strategy, clusters of 1 (feature) are displayed using the default icon** (yellow disk). i would rather not use any plugins or libraries that aren't standard openlayers. **any suggestions?** below a part of the original javascript code using clusters (remove the cluster strategy declaration `new OpenLayers.Strategy.Cluster()` and custom icons are displayed just fine):               var urlKMLClient = 'KMLClientsAll.kml';      var layerKMLClient = new OpenLayers.Layer.Vector("Clients", {             strategies: [new OpenLayers.Strategy.Fixed(), new OpenLayers.Strategy.Cluster(), refresh],             protocol: new OpenLayers.Protocol.HTTP({                 url: urlKMLClient,                 format: new OpenLayers.Format.KML({                     extractStyles: true,                      extractAttributes: true,                     maxDepth: 2                 })             })         });       ![WITHOUT CLUSTER](http://i.stack.imgur.com/3muMR.jpg) see Client custom icon (green marker) gives way to cluster default icon (yellow disk) below ![WITH CLUSTER](http://i.stack.imgur.com/OTnYY.jpg)


title: Geoserver + MongoDB: connection OK but layers list empty

text: I setted up Geoserver (ver. 2.5.2) with the MongoDB plugin (ver. 11.2 - http://repo.opengeo.org/org/geotools/gt-mongodb/11.2/). The connection between Geoserver and MongoDB is OK and when I create a new datastore in geoserver->stores->Add new Store, it creates the database on mongo as expected. The problem is that, if I go in layers->Add a new resource and select the correct resource, the layer's list remain empty. On the top i see only "You can create a new feature type by manually configuring the attribute names and types. Create new feature type..." but i cannot find documentation on how to use that. Any help? Thanks in advance **EDIT:** in the goserver's log i've just found this error:               ERROR [data.mongodb] - getLayers error; command failed [mapreduce]:          { "serverUsed" : "[ip:port]" ,          "errmsg" : "exception: ReferenceError: mapfields_recursive is not defined near 'ction()         { mapfields_recursive (\"\", this);' " , "code" : 16722 , "ok" : 0.0}       but I can't still figure out what is the problem and how to solve it. Any ideas?


title: How can I extract variable information from an OpenLayers 3 object?

text: I have tried to make a simple map with OL 3's Geolocation module. I succeeded with a resource wasting method. There must be a simpler way to extract variable information from an OpenLayers 3 object. My code is the following:               var view = new ol.View({         center: [0, 0],         zoom: 16         })          var map = new ol.Map({         target: 'map',         layers: [             new ol.layer.Tile({                 source: new ol.source.OSM()             })         ],         view: view         });          var geolocation = new ol.Geolocation({         projection: view.getProjection()     });          geolocation.setTracking(true);          var geolocation2 = new ol.Geolocation({         projection: view.getProjection()     });          geolocation2.setTracking(true);          var featureSource = new ol.source.Vector({     });          var feature = new ol.Feature();          feature.bindTo('geometry', geolocation, 'position')         .transform(function() {}, function(coordinates) {           return coordinates ? new ol.geom.Point(coordinates) : null;         });               featureSource.addFeature(feature);          var featureLayer = new ol.layer.Vector({         source: featureSource,         name: 'Your position'     });          view.bindTo('center', geolocation2, 'position');          map.addLayer(featureLayer);      The code works, but not the way it should operate. It uses two `Geolocation` instances and because the `bindTo` method the view periodically resets. I couldn't manage to extract the `position` variable of the `Geolocation` object. I've tried `Geolocation.getPosition()`, `Geolocation.get('position')` and `Geolocation.p.position`. Every one of them acted like a void method and returned an `undefined` value. I use Google Chrome version 35.0.1916.153 on a Debian Wheezy. To complete the story, if I use only one `Geolocation` instance and bind it to the `feature` and the `view` objects, the point follows the center of the map with every change in the view.


title: Non-spatial Data Driven Pages

text: Is it possible to utilise the Data Driven Pages (DDP) tool in ArcMap or ArcPy to create a set of maps with the same extent but showing different layers? To clarify I regularly make large sets of maps for projects all in the same location often with the same baselayers but modified to focus on different aspects for the relevant section of a report such as monitoring locations, land-use, key features etc. which all require separate .mxds to be created and maintained. A DDP type system would allow me to utilise the dynamic text and ensure a standard template / output throughout, it would also make life a lot easier when updating a general figure number or comment as a report is written in one .mxd rather than multiple for a client. I'm using ArcMap 10.1, but having read through the ArcGIS DDP online help I assume this would have to be done with ArcPy as I've been unable to figure out how to do it within ArcMap. I'm aware that QGIS has the Atlas Plugin, although I'm not as familiar with its capabilities but would happily use it as an alternative.


title: How can you create a depth raster from a DTM and a polygon in QGIS?

text: I have a DTM raster and a flood outline in polygon format - how do I create a depth raster for the flood outline (i.e. flood level minus DTM level)?


title: Looking for system specifications for a modern GIS workstation

text: My current employer has given me the opportunity to design my own GIS workstation. The problem is I have no idea what is important when it comes to hardware! I have been given a datasheet from Fujitsu which includes the latest CELSIUS M720 series workstations but the configuration includes about 30 different processors, graphics cards, and hard drives. I will be working a lot with the ESRI desktop suite (spatial analyst and 3D analyst) and doing a fair amount of raster processing (DEMs, visibility analysis, etc.), but I also use the GDAL library, SAGA-GIS, Quantum GIS and others. Does ArcGIS 10 actually use multiple cores? I have the choice from 4 cores with 8 threads up to 8 cores with 8 threads - 16 threads? Is SSD SATA III better than HDD SATA III? The Celsius R920 supports up to 512 GB RAM. But what do I need? It will probably be the last time I get a new workstation for the next 5 years, so I am curious what others have ordered or feel will be necessary for the coming years. Obviously the cloud isn't taking over as quick as they say!


title: System Requirements for Quantum GIS

text: > **Possible Duplicate:**   >  Looking for system specifications for a modern GIS workstation I have been using the program to make maps to include in my reports over the past year or so. It is excellent. However, my computer seems to not be able to handle all the data sometimes, and the program crashes. It happens when I do too much editing to data points that I have manipulated. My question is what system requirements are recommended to run this program with a lot of data points? Or does it simply depend on what I want to do with the program? I am probably going to get a faster computer primarily based on wanting to use this program. Thanks Marcus


title: Which processor for GIS Desktop?

text: > **Possible Duplicate:**   >  Looking for system specifications for a modern GIS workstation I am picking a new pc configuration for my GIS workstation. I use ESRI ArcMap and ArcCat geoprocessing tools, together with PG Admin III, SAGA-GIS. Quantum GIS, Adobe Photoshop/Illustrator as well as command line tools like FW-Tools, Putty and text editors....Not necessarily all at the same time, but more often than not there are between 5 and 10 programs running at the same time. I am configuring a new workstation and have a list of processors which are recommended by Fujitzu...namely the Intel Xeon E5 range. What I am not sure about is how many cores/threads I need. Does each application use a thread?...or even a core? Should I get as many cores as I can afford? any advice most welcome, Yours,


title: Covering a Polygon with Minimum number of circles with fixed locations

text: In QGIS, I have a vector polygon of a fixed size and shape. I also have a set of circles of a fixed radius and fixed location, which cover part of the polygon. What I would like to do is remove as many of the circles as possible while leaving as few gaps as possible, a gap being an area not covered by the circles. No circle should have overlap with another circle of more than 90% (or 80% or 70%- I am trying to run a few scenarios) of its area. I would also like to treat circles which go off the edge of the polygon in the same way as overlap, i.e. a circle with more than 90% of its area not covering the polygon would be eliminated. Inevitably there will be times when two or more circles could be chosen, I would like to prioritize them based on specific attributes. Here is a picture of the complete set of circles(pink) and the polygon I need them to cover (blue). ![enter image description here](http://i.stack.imgur.com/dkDlj.jpg) I am new to GIS and any extra elaboration on how to implement is greatly appreciated.


title: cast from IFeatureClass to ISchemaLock

text: I'm programming with the ARCGIS Java SDK, in Eclipse. I need to cast an **IFeatureClass** to an **ISchemalock**. This cast is supposed to be allowed. There is an example of this cast in: http://resources.esri.com/help/9.3/arcgisengine/java/doc/b22267cb-642c-11dc-9ca3-0b35f906bb2e.htm However, when I try to perform this cast in my code, an exception appears: _Exception in thread "main" java.lang.ClassCastException: com.esri.arcgis.geodatabase.IFeatureClassProxy cannot be cast to com.esri.arcgis.geodatabase.ISchemaLock_ Are you able to make this cast on an IFeatureClass ? Here is the code to perform this:               //Attempt to acquire an exclusive schema lock on the feature class.     IFeatureClass fc = fw.openFeatureClass("featureclass name");     ISchemaLock schemaLock = (ISchemaLock)fc ;      I would appreciate if you could try this simple line of code in one of your programs. Keep me posted... **EDIT:** I wanted to see if an IFeatureClass was, at least, recognized as an instance of ISchemaLock. Believe it or not, it was NOT recognized as an ISchemaLock! The execution of the code below did not enter the **if** statement.               if (fc instanceof ISchemaLock)      {       System.out.println("all ok");       ISchemaLock schemaLock = (ISchemaLock)fc;     }      Maybe is there something not correctly installed with Eclipse ?


title: Styling a WFS layer, how to control title in legend?

text: I have a problem with styling a WFS in OpenLayers. I tried different ways (from the examples) to style the WFS, but I could only get the default style in orange. Then I found one way, which worked for me. The problem is that in the legend is written "Untitled 1" next to the now correct symbol as with the default style. How can I change this title or remove it? Has anybody have a solution on how to style the WFS with an external SLD? I also tried different examples found on the web but I couldn't get it to work. Here's the code I am using to style the WFS but that displays the wrong title:               var StyleObjekt = new OpenLayers.StyleMap({                                 "default": new OpenLayers.Style({                                                               fillColor: "#fff000",                                                               strokeColor: "#666666",                                 }),                                 "select": new OpenLayers.Style({                                                                 fillColor: "#ffffff",                                                                 strokeColor: "#666666",                                                                 }),                                 "temporary": new OpenLayers.Style({                                                                 fillColor: "#ffffff",                                                                 strokeColor: "#666666",                                                                 })                 });            var gwg_f_wfs = new OpenLayers.Layer.Vector("Flurstuecke_WFS", {                 styleMap: StyleObjekt,                 strategies: [new OpenLayers.Strategy.BBOX()],                 protocol: new OpenLayers.Protocol.WFS({                     version: "1.1.0",                     url: "http://localhost:8080/geoserver/wfs",                     featurePrefix: "NewGeo",                     featureType: "Flurstuecke_GWG",                     featureNS: "http://www.newgeo.de/GeoServer",                     srsName: "EPSG:2398"                 })                           });


title: How to install PostGIS correctly with OpenGeo Suite 4.0 community edition on Windows 7?

text: I have installed the OpenGeo Suite community edition, version 3, onto Windows 7 multiple times in the past, with virtually no problems. With the latest version of OpenGeo Suite (4.0.1) (which includes PostGIS 2.1), I am able to install almost all of the components successfully, except PostGIS. With earlier editions, PostGIS was installed as part of the full Opengeo installation process, but now, although PostGreSQL 9.3 is successfully installed, there are no PostGIS or opengeo database templates installed, nor are the spatial functions installed. Nevertheless PostGIS 2.1 is installed. It seems that OpenGeo Suite 4 intalls all the component but no longer creates a default database or user account in Postgresql/PostGIS. What steps are necessary to configure OpenGeosuite (windows installation) in order to be able to fix this issue? thanks


title: Extract points from regular data.frame with coordinates based on SpatialPolygonsDataFrame in R

text: I have a regular data.frame called datos.merged > class(datosmerged)               [1] "data.frame"      The last two columns are the corresponding geographical coordinates > datosmerged[1,c(23,24)]                  XCOORD   YCOORD     1 2473520   668648.7      I have a file called shape > class(shape)               [1] "SpatialPolygonsDataFrame"     attr(,"package")     [1] "sp"      I would like to use this very simple shape to subset the previous data.frame. How can this be done?


title: How to create links between points and nearest vertex of polyline feature class?

text: I have a point feature class and for each point I want to find the nearest vertex on a polyline feature class. I then want to create a line between each point and it's nearest vertex on the other layer. How is this done? If it can't be done with a polyline layer, I can also use a point layer in lieu.


title: map() Function and IE8

text: I am using Javisantana example of Autocomplete an everything works inside my map, until I tested my map with IE8 and experienced this error:               Object doesn't support property or method 'map'       After some research it seems as if IE does not support that function some of the hacks/patches that I have tried fix that error but seem to cause errors within cartodb js. Has anyone ran into this problem? If so any tips or advice on a work around?


title: does cartodb works offline

text: i had installed cartodb on a local machine. It works fine. But when i tried to go off-line, and then i was able to login, nothing can be done. It seems page is looking for some on-line data's. Is it possible to completely work **cartodb in an off-line mode**