title: GDAL Installation help in windows

text: I'm trying to build a GIS application in Django Framework.   1. I completed the GeoDjango installtion (postresql, spatial libraries).    2. I installed OSgeo4W in my C drive.   3. I also setup environment variables.    4. I'm using PyCharm as an IDE and have a virtual env set up for my project.  But I cannot get the osgeo import to work in my django project. It's not recognized but the import works fine in python shell. Do I need to install all these gis components into my virtualenv site- packages? How do I get around this?


title: Converting a JOSM indoor map to QGIS

text: I have created an OpenStreeMap layer of the indoor of a bulding by following this tutorial: https://wiki.openstreetmap.org/wiki/IndoorOSM, since that is only a proposed tagging schema for indoor mapping when I have tried to load it into QGIS the only thing I got was the building and not the indoor nodes. How can I turn those rooms into something that can be interpreted by more advanced tools like QGIS or a postgis database? PS: My project consists on providing routes between two different rooms of a building, any help or suggestion on how to do it or how to render osm's indoor layers is welcome.


title: How to add .shp file to ArcMap?

text: Ive spent the last 12 hours trying to figure out how to use ArcMap. Ive made about 2% progress. I attached the link to the .shp files at the bottom. Ive been exporting .shp files via Pathfinder Office. Every time I export this file it either exports in the middle of Atlantic or somewhere in the middle of Ohio. I am working in southern PA. This is the play by play. In Arc Map, I add new base map. Save the projection- WGS_1984Web_Mercator_Auxiliary_Sphere. From there I export my data from terrasync to PFO run the dif cor. and export it to a .esri shapefile. when I export the .shp file I make sure to select properties and browse to the projection file I just saved. If I can get past this step. I need to use the kriging tool to add colors to all 5000 points I logged. I would like to have a key something like Green 0-12000 Yellow 12001-18000 Red 18001-100000 Ive done quite a bit of research on kriging but there is several different tools to use and im not sure which one to use. https://sites.google.com/site/strubegps/kmlfiles


title: Is it possible to calculate a value based on values of another row?

text: I have a point shape file in qgis. For one of the attributes, I would like to offset the values, i.e. the value for point 1 goes to point 10, the value for point 2 goes to point 11 and so on. I see the $rownum attribute in the field- calculator, but I cannot see any examples on how (if possible at all) I can use that to access values in another row... Is this possible using the field calculator or do I have to make a python script? (Yes, I know for the example, rows 1 - 9 will get null- values and the last 9 values will "fall off the end") (Using qgis 2.0.2)


title: Requesting a JSON format from Geoserver using Openlayers and GeoExt

text: I am requesting a vector layer from geoserver using WFS request in JSON format. I have this url http://myip.com:8001/geoserver/workspace/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=workspace:mylayer%20layers&maxFeatures=50&outputFormat=application/json. I want to populate my grid/table in my client side bind to the vector layer. What happens is that I won't be able to see my layers. My question is that, am I still going to use a server proxy? I am using Tomcat 6 in Windows 7. Try to see my code:               var vector = new OpenLayers.Layer.Vector("My Layer");     store = new GeoExt.data.FeatureStore({         layer: vector,         fields: [             {name: 'name', type: 'string'},             {name: 'description', type: 'string'}         ],         proxy: new GeoExt.data.ProtocolProxy({             protocol: new OpenLayers.Protocol.HTTP({                 url: "http://myIP:8001/geoserver/workspace/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=workspace:mylayer%20layers&maxFeatures=50&outputFormat=application/json",                 format: new OpenLayers.Format.GeoJSON()             })         }),         autoLoad: true     });          var gridPanel = new Ext.grid.GridPanel({         //title: "Feature Grid",         region: "east",         store: store,         width: 'auto',         height: 100,         columns: [{             header: "Name",             width: 200,             dataIndex: "name"         }, {             header: "Description",             width: 100,             dataIndex: "description"         }],         sm: new GeoExt.grid.FeatureSelectionModel()      });


title: How to do geometric correction of remote sensing data in R?

text: My challenge is to replace most of the remote sensing commercial tools with R code (I'm teaching assistant at the University and I would like to move to open source R code to show the added values to the students). I'm currently stuck as I cannot find a way to perform **geometric correction** of remote sensing data. The exercise should be:   1. Select **Ground Control Points** on both the georeferenced image and the image to correct    2. Use polynomial mapping for image correction (students should try to test different **polynomial order** and calculate the **RMSE** for each step)    3. **Warp** the image in a defined projection using the interpolation method (either nearest neighbour, bi-linear interpolation, cubic convolution) The 2 last steps could be processed using GDAL but it is not really interactive with R. Does someone has any idea or experience with such processing?


title: No module named pyspatialite- plugin:Qspatialite, OpenSUSE 13.1

text: I cannot run the plugin because of module missing. I can't find that module in the OpenSuse repository, and python community repository. I'm using QGIS 2,0 and OpenSuse 13,1 Where can I find module pyspatialite with dependencies, because I've found that module at: packages.altlinux.org/en/Sisyphus/srpms/python-module- pyspatialite/get but if I try install it I get: nothing provides python-modules needed by python-module- pyspatialite-3.0.1-alt1.i586. Module psycopg is installed in system - with no effect. Thanks for answer, Tom


title: How can I add height data to a polyline in QGIS?

text: I want to place a polyline at the correct altitude in QGIS. Is that possible and how do I do this? What I have is a polyline (railroad) and I want to add information on the electrics above the railroad. Say, the electrics are z meters above the rails, how do I get this data into the attribute table? I want to use this information to calculate the number of trees that can reach the electrics if they fall. And yes, I get that they would reach the rails even if they aren't tall enough to reach the electrics. EDIT: I forgot to add one more thing - I'm pretty new to the world of GIS. Any ideas on how to do this? Thanks. EDIT: I've tried to use the tip I got in the comments below, thank you for trying to help, but it was a bit to complicated for me, I simply don't understand how to do it. I just want to add values to the z-coordinat (of a polyline). Please, help...


title: Reshape tool stops working in QGIS 2

text: I am using QGIS 2 on Windows 7. I am reshaping some polygons which I am supposed to correct with the reshape tool in the "advancing digitization" bar (![reshape_tool](http://i.stack.imgur.com/Ast2c.png)). Somehow, after some changes are made successfully, the tool sometimes stops working, and I have to delete the polygon and create a new one in order to finish my work. This is time consuming as I have plenty of shapes to correct. Any clue why this occurs? Is it a bug which someone else is experiencing as well? Any help would be very appreciated.


title: Clipping without cutting features

text: I want to create a multimodal network dataset for the city of Montreal. I have roads and transit lines features for Quebec province (where Montreal is). Because I need the features limited to Montreal only, I clipped them using Clip tool in geoprocessing menu with Montreal census tracts layer as the clipping feature. Now, the problem is that Montreal has water bodies which are not a part of census tract layer. The areas of water bodies are therefore blank spaces between census tracts. The clipping procedure cuts the roads and transit lines where they pass over the water bodies. This, obviously, can't be used for creating an accurate network datset. Which tool can I use to clip the road and transit lines feature classes to city of montreal census tract layer only, without cutting them? Alternatively, is it possible that I create network datset for the whole province and then clip it?


title: how to add orthophoto layer in qgis

text: I am not an advanced user, but I try to be oneday :) here is my issue: I have tracks (just waypoints) from my GPS in gpx format. I also have orthophoto pictures in tif format, that I would like to make a layer with I use Quantum GIS (1.8.0) actually i know how to put it in, but it shows wrong (coordinate system = WGS 84/EPSG 4326 - europe) my goal is to put together all waypoints and export as a shp file... thank you so much!


title: Can I convert an ISegmentCollection to an IEnumFeature?

text: I would like to try to call the function, `IFeatureConstruction.SplitPolygonsWithLines( )`; The 4th arg is a type `IEnumFeature`, I have an `ISegmentCollection`. Can I convert ISegmentCollection to an IEnumFeature?


title: How to prevent a geoprocessing task from being run twice?

text: If I run twice a geoprocessing task (using the REST API), what does happen? Does the geoprocessing service execute the two tasks simultaneously, ...? How can I configure the geoprocessing service from executing two tasks at the same time?


title: Group Layers Do Not Display in GeoExt LegendPanel

text: I have a GeoExt.LegendPanel which is working just fine for layers which are not grouped. However, the grouped layers do not show any icon in the LegendPanel. I have tried to create the group layer various ways, using the exact name of the layer group in Geoserver and aldo by specifying a layer string in code. I've been looking around for a while now, but I don't see what I am missing.. Here is relevant code:               //Layer Group     var geography2 = new OpenLayers.Layer.WMS(                         "Geography2", "http://localhost:8080/geoserver/wms", {layers: 'cite:geog_multipol_vw,cite:geog_multiline_vw', transparent: true}, {isBaseLayer: false, displayInLayerSwitcher: true, visibility: true, minScale: 27734017.045258757}                    );          //LegendPanel config     var legend = new GeoExt.LegendPanel({                     region: "east",                     title: "Legend",                     width: 200,                     autoScroll: true,                     layerStore: mapPanel.layers                 });          //Explicitly point LegendPanel to GetLegendGraphic on WMS Geoserver for GeoWebCached Layers                 var LayerRec;                 for(var l = 4; l < map.layers.length; l++)                 {                     LayerRec = mapPanel.layers.getAt(l);                     LayerRec.set("legendURL","http://localhost:8080/geoserver/wms?TRANSPARENT=true&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetLegendGraphic&EXCEPTIONS=application%2Fvnd.ogc.se_xml&FORMAT=image%2Fgif&LAYER=" + map.layers[l].params["LAYERS"]);             }      Thanks for your time in trying to help me :), elshae


title: What's the easiest way to move tile-based services?

text: We're still running ArcGIS 9.3.1 for a short while longer, and have been tiling more of our data, which is great until we have to deploy it. Deploying currently consists of transferring the entirety of the tileset to the production machines, which takes about half of forever, especially on a deadline. We're currently using tar.gz, but would welcome a better approach. Tiling in place is not an option as there is not a "spare" server at the production facility. Edit: I should mention that the destination server is a _long_ way away, and behind a firewall, meaning that our access to it is limited, and mapping a drive or a share is Right Out. I'd be tempted to use external drives and FedEx, but that would require a physical person at the destination, which tends not to be the case for months at a time.


title: Assessing OpenStreetMap coverage for international routing

text: I posted the following question on the main StackOverflow website, but was informed that I might get a better response here. Hence I'm copying it below. Note however, that I have an active bounty on the original posting. * * * I have been using a commercial solution for route distances and travel times for North America and Western/Mid Europe. I am considering expanding the project to cover other countries - and perhaps the entire world. A very limited budget and patchy regional coverage from individual commercial providers, probably make locally-hosted OpenStreetMap the only viable option. Before someone suggests an online solution, my application requires a lot of intensive route calculation - something which would cost a lot or be very impolite (and probably banned) if performed using a web service. The results of the calculations are put back in the public domain, so rediting OpenStreetMaps is not a problem. My problem is how do I assess the routing data coverage for individual countries in the OpenStreetMap database? Such an assessment could determine if the project is viable, and a suitable order for processing (ie. do the countries with the best coverage first). High-end commercial data providers can typically supply statistical descriptions, as well as regional descriptions of surveyed coverage. OpenStreetMap is much more patchy - an area typically includes some roads, but not all roads. Individual location errors of a few metres of even 10-20m will not be a problem for my application (I'm looking at city-city distances), but route graph connectivity is. Ie. the road vectors must logically meet correctly at a junction. Has anyone attempted to create statistics describing data coverage of the OpenStreetMap database? If not, how would you go about it? The best I can think of is to take a random sampling of places (eg. cities), and then attempt to calculate routes. There would have to be an assumption that the major roads will tend to be added before the minor roads. Therefore a route between two distant cities would use the logical major road, and not a minor road (which is typically longer/slower) because the major road is missing. Another problem would be that it is physically not possible to drive between many towns. Often this is due to the presence of islands (where ferries could be used) but often there is no surface route (eg. settlements in Nunavut). So how would such statistics be used when comparing between (say) Tonga and Afghanistan. Afganistan probably has very low data coverage. Tonga is probably better but the settlements are spread out across an archipelago. Some details about my application: All start and end points are towns and cities with locations taken from the Geonames database. Typically I am looking at the 1000 largest cities in a country that also have a population of at least 1000. Routes are currently calculated in duplicate as both fastest routes and shortest routes. Reasonable road speeds vary according to broad road categories. Estimated travel times are computed alongside road distances. These details are preferences for consistency- they are not set in stone.


title: Examples of beautiful statewide maps that have a lot of layers (20+ layers)

text: I'm in the process of building a statewide map of all the transportation modes, waterways, and other information for a whooping total of 40 layers. I was wondering if you have seen print maps that have many displayed layers in a professional and beautiful manner, and if you can share it, or at least provide some of the design aspects that made the map so "pretty". My question is different from beautiful maps as I'm only interested in printing them and not putting them on the web


title: send viewparams geoserver using leaflet

text: I need to pass my longitudes and latitudes of origin and destination for a sql view in geoserver to get a WMS layer with the route. I've performed in qgis to view and is working perfectly. This is the view in geoserver: ![enter image description here](http://i.stack.imgur.com/OqW7N.jpg) This my code to add the layer on the map, I'm using leaflet.               function findRoute(){        route = L.tileLayer.wms("http://localhost:8888/geoserver/pgrouting/wms", {          layers: 'pgrouting',          viewparams: [x1:paramOrigemLng;y1:paramOrigemLat;x2:paramDestino1Lng;y2:paramDestino1Lat]         });     }      Can anyone give me an example of layer viewparams the key format: value with leaflet because I'm not hitting the syntax. I'm doing as the link: http://docs.geoserver.org/stable/en/user/data/database/sqlview.html


title: How to select polygon close to each other?

text: I have grass patches and I put a buffer of 250m around them. I did a "dissolve to get larger polygons that overlap, but the attribute table emitted a lot of information, which are important such as the hectare och area. How do I dissolve, and still keep these vital information. Now the next step I want to accomplish is to "merge" grass patches that are less than 500m apart into a larger polygon. I'm having a very hard time doing this. I'm using ArcGIS for Desktop 10.1. What tools would be recommended?


title: maximum limit for plotting lat-long on google maps

text: I want to highlight a route on google maps along with our network sites. I have successfully done it but when I load data in large proportion then my plugin becomes unresponsive. Can someone suggest me a way to avoid this error.


title: How to select random adjacent census tract

text: I have set of census tracts spaced out across a city and I would like to randomly choose an adjacent census tract to explore nearby demographic data. I'm not sure if there is an easy process for getting adjacent census tracts for selected census tracts.


title: Whats an appropriate map projection for georeferencing aerial photos?

text: I have some aerial photos that were shot from directly overhead a 5x5 mile property. They were then mosaiced to create a continuous overhead tiff of the entire property. I am now going to georeference in ArcMap using satellite imagery as my base reference. Is there one map projection that would be more appropriate than others for using as reference? Edit: By appropriate I mean - which projection will results in the most visually similar look to an overhead, aerial photo? Using WGS84 is too wide and stretched.


title: How to project pixel coordinates in leaflet?

text: I am trying to create a map based on an image which is 16384x16384 px, but I would also like to add markers at specific locations using pixel coordinates of this given image. I created a tile layer, a map element and set maximum boundaries, so that I can't scroll out of the image, using this code:               var map = L.map('map', {         maxZoom: 6,         minZoom: 2,         crs: L.CRS.Simple     }).setView([0, 0], 2);     var southWest = map.unproject([0,16384], map.getMaxZoom());     var northEast = map.unproject([16384,0], map.getMaxZoom());     map.setMaxBounds(new L.LatLngBounds(southWest, northEast));          L.tileLayer('tiles/{z}/{x}/{y}.png', {         tms: true     }).addTo(map);      The points i'd like to add to the map are stored in an external file as GeoJSON, and they look like this http://www.de- egge.de/maps/terranigma/terraPoints.js I load them using this snippet of code:               var terraPoints = new L.geoJson(points, {         onEachFeature: function (feature, layer) {             layer.bindPopup(                 "<b>Ort: </b>" + feature.properties["points[, 1]"]              );         }     });          map.addLayer(terraPoints);      But of course, they don't show up, because the reference systems don't match. The points use pixel coordinates, but my map uses geographic coordinates. Is there a way to "unproject" the points while loading them? Any help is much appreciated and thanks in advance!


title: How to integrate info Z into a point shape, no attribute table

text: Firstly, I use ArcGIS-9.3 and I have a point shape, but when I try to get info about its value Z (Sketch Properties) the result is "0" for ever points, like this: X:274988,260 Y:4126481,508 Z:0 M:NaN My question is how can I pass Z info to my points from a DEM for example??. But I would like to note that I don´t want the Z info on my attributes table, but I want it in the Sketch properties. Thanks.


title: Select by SQL in QGIS Python

text: In the GUI there is the Query Builder, but that does not seem to be easily available in Python. The processing (SEXTANTE) module has the `qgis:selectbyattribute` algorithm, but that only supports selecting by a single attribute. There is also `QgsFeatureRequest.setFilterExpression`, but apparently that is not available in current stable version (2.0.1 Dufour). Do I have to write my own?


title: Is Read Only nature of SQLite support in ArcGIS for Desktop 10.2 documented?

text: I've recently started creating **SQLite** geodatabases using **ArcGIS for Desktop 10.2** with both **ST_GEOMETRY** and **SPATIALITE** storage data types and was impressed by how easy this was to do. I then tried to **Start Editing** them only to be blocked by an error " **No editable layers** " and a warning " **The workspace containing this data cannot be edited** ". ![enter image description here](http://i.stack.imgur.com/K5cfl.png) I looked and could not find this limitation mentioned anywhere. Is the **Read Only** nature of ArcGIS support for SQLite documented?


title: Can Arcmap 10 do change detection?

text: I have an assignment to find out changed areas of the two different maps (2000 and 2010) by using Arcmap 10. Is it possible to use it?


title: Point count within borders, in two different layers, using QGIS

text: I have a layer of georeferenced points, and a layer of vectorized borders. Is there a way I can use QGIS to calculate how many points fall within each area of the borders layer?


title: How to add report to layout using Arcpy

text: I have a script that creates MXDs and layouts dynamically, i need to to add to each data frame in the layout report using arcmap report wizard rather then export to pdf, is it possible ?


title: Plot a EPSG:4326 KML File using Google Maps V3

text: I have a KML file that uses EPSG:4326. Is there any way to plot this onto a Google Map (v3) which uses Mercator? Here is the page I am working on The district data comes from one source and is Google friendly. The outage data is generated by the power company using their GIS system. It shows a list of power outages. However nothing is shown on the map. Have seen this file work in openlayers but I really have no clue why it would not be showing up.


title: Is there a way to set rotation on inset maps to match data driven maps

text: I have been playing with the new map book functionality in ArcGIS 10 to re- create maps I have made in the past. The new Data Driven maps with page queries is a pretty powerful tool. Is there a way to match the rotation of inset map, secondary data frame, with the main layer being used for the data driven maps? I found it is possible to match the extent in the 'Data Frame' tab of the 'Data Frame Properties', is this also the place to find the rotation option? If this is something that you need to code, is it possible to build it into the mxd and to have it launch on refresh? Not sure I understand how Python is supposed to work within a mxd. Thanks for your time, Guido


title: selecting by attribute and specific value replace

text: I want to select all of the features by attribute where the attribue in myattrib is like... Selection "criteria H". That is not a problem. However then I need to replace (myattrib," H", "") except I only want to replace the H at the end of the attribute, not within, like "toast is Hot H". In this record I want the value to end up as "toast is HotH" There are 15k records but only 4800 have the value I am replacing. And it is always the same value. edit response: No. All values I want to fix end in H with a space in front of it. I am depending on search replace but not sure how to leave the space in front of Hot. In my example I am able to select with the "% H" which specifies there is nothing after the H so the select won't get Hot even though there is a space before it. However when it comes to replace it does a literal replace on " H" and doesn't care about the ot. I think if I knew some regex and "IF" arcmap supports that, It might do the trick (from what I have seen on the FMEtalk forum). I'm sure someone can talk me through this.


title: Flow direction generated for river is incorrect

text: I am trying to generate a hydrologic model and I can not really proceed on the next processes. I have this problem on generating flow direction. I have already executed fill sinks and the next is flow direction of the grid. However, the generated flow direction grid is not the actual flow of the river. it should be flowing from the highest elevation to the lowest, right? The river should flow out towards the sea, but the direction of the flow is the opposite. So I was thinking what could be instances that could affect why did this happen. I figure out people from this community could be of help to me. I am using ArcMap 10.1. Here are my screenshots: fill sinks ![fill sinks](http://i.stack.imgur.com/PnrV9.png) flow direction ![flow direction](http://i.stack.imgur.com/HxoHr.png) flow direction values ![flow direction values](http://i.stack.imgur.com/2tmUv.png) flow accumulation (flow break) ![flow accumulation \(flow break\)](http://i.stack.imgur.com/2y6bK.png) The yellow arrow tells the direction of the actual flow of the river. The red dashed arrow tells the direction of the generated flow direction. ![enter image description here](http://i.stack.imgur.com/s5egW.png)


title: Positioning and tracing polygons over an image

text: On QGIS 2.0 - I want to add an image (TIFF or whatever) then resize and correctly position it over a satellite image (openlayersplugin). Then, the task is to create a new shape file layer and draw over the lines in the image which represent the TPUs (Tertiary Planning Units) of Hong Kong. This data is actually for sale from the Hong Kong Planning Department, unfortunately, specific licensing codes require any work done with the data to be completed within Hong Kong. Hence why I plan to simply use the PDFs they posted online to create my own. Which they don't mind me doing, in fact, they even described a way to extract the DWG imbedded in the PDFs on their website for guidence. If you're interested, this will be used to complete some public transport accessibility on the city to be compared with a number of other cities (www.snamuts.com). We have a draft submission deadline coming up soon so I need to complete this task asap. I have tried adding a raster layer and selecting the TIFF but the image just doesn't appear... Even if that did work, is there a function to shift and transform the image so it aligns with the satellite base map below?


title: GDAL: Update spreadsheet from polygons in a shapefile

text: I am a beginner with GIS data, and have been looking for an answer but have not yet found it. I have a spreadsheet with lat/long points. I also have a shapefile with several numbered polygons over a geographical area. I would like to do the following:   1. For each row in the spreadsheet (lat/long point) determine which polygon it is within   2. Update the spreadsheet (or generate a new spreadsheet/shapefile/etc.) with each row having its polygon number in a new column. I suspect I could possibly use ogr2ogr to do this, or possibly write some Python code to do it, but I really don't even know where to start, so any help at all would be appreciated.


title: How to apply Mapnik style for OSM data in ArcMAP?

text: I am loading OpenStreetMap vector data into ArcMAP (using the Download and Symbolize OSM Data tool of ArcGIS Editor for OpenStreetMap 1.1). I'm sure there must be a simple way to apply the style/symbols of the Mapnik baselayer as rendered on the OSM website (example), rather than using the ArcMAP defaults or manually adjusting symbology for each feature type, but I just can't figure it out. Can someone point me to a potential solution? Cheers


title: Why isn't my ArcGIS JavaScript QueryTask working?

text: var map;         require([                 "esri/map",                  "esri/layers/ArcGISDynamicMapServiceLayer",                  "esri/layers/ImageParameters",                  "esri/tasks/QueryTask",                  "esri/tasks/query",                  "esri/symbols/SimpleMarkerSymbol",                 "esri/InfoTemplate",                  "dojo/_base/Color",                  "dojo/domReady!"                 ], function(                 Map,                  ArcGISDynamicMapServiceLayer,                  ImageParameters,                  QueryTask,                  Query,                  SimpleMarkerSymbol,                  InfoTemplate,                  Color                 ) {              map = new Map("mapDiv", {             center: [-123.8425, 45.4552],             zoom: 16,             basemap: "streets"             });              var imageParameters = new ImageParameters();              imageParameters.layerOption = ImageParameters.LAYER_OPTION_SHOW;         imageParameters.layerIds = [4];         imageParameters.transparent = true;              var dynamicMapServiceLayer = new ArcGISDynamicMapServiceLayer("http://host:6080/arcgis/rest/services/Data/MapServer", {               "opacity" : 1.0,               "imageParameters": imageParameters             });              map.addLayer(dynamicMapServiceLayer);              queryTask = new QueryTask("http://host:6080/arcgis/rest/services/Data/MapServer/4");              query = new Query();         query.returnGeometry = true;         query.outFields = ["OBJECTID"];              infoTemplate = new InfoTemplate("${OBJECTID}", "OBJECTID : ${OBJECTID}");              symbol = new SimpleMarkerSymbol();         symbol.setStyle(SimpleMarkerSymbol.STYLE_SQUARE);         symbol.setSize(1000);         symbol.setColor(new Color([255,0,0,1]));              query.where = "OBJECTID = 1";              queryTask.execute(query, showResults);              function showResults(featureSet) {             //remove all graphics on the maps graphics layer             map.graphics.clear();                  //Performance enhancer - assign featureSet array to a single variable.             var resultFeatures = featureSet.features;                  //Loop through each feature returned             for (var i=0, il=resultFeatures.length; i<il; i++) {                 //Get the current feature from the featureSet.                 //Feature is a graphic                 var graphic = resultFeatures[i];                 graphic.setSymbol(symbol);                      //Set the infoTemplate.                 graphic.setInfoTemplate(infoTemplate);                      //Add graphic to the map graphics layer.                 map.graphics.add(graphic);             }            }          });      I followed this ESRI example: https://developers.arcgis.com/javascript/jshelp/intro_querytask.html But when I test it I get an error:               Resource interpreted as Script but transferred with MIME type text/plain: "http://host:6080/arcgis/rest/services/Data/MapServer/4/q…Fields=OBJECTID&callback=dojo.io.script.jsonp_dojoIoScript3._jsonpCallback". init.js:495     Resource interpreted as Script but transferred with MIME type text/plain: "http://services.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer?f=json&callback=dojo.io.script.jsonp_dojoIoScript1._jsonpCallback". init.js:495     TypeError {stack: (...), message: "Cannot read property 'clear' of null"}      "TypeError: Cannot read property 'clear' of null         at showResults (http://test/test/script.js:61:15)         at e._successHandler (http://js.arcgis.com/3.10/init.js:490:193)         at e._handler (http://js.arcgis.com/3.10/init.js:1424:385)         at http://js.arcgis.com/3.10/init.js:174:23         at Object.d.load (http://js.arcgis.com/3.10/init.js:1420:411)         at http://js.arcgis.com/3.10/init.js:630:478         at c (http://js.arcgis.com/3.10/init.js:74:221)         at d (http://js.arcgis.com/3.10/init.js:74:10)         at resolve.callback (http://js.arcgis.com/3.10/init.js:75:350)         at c (http://js.arcgis.com/3.10/init.js:74:436)"       I know there should be a query result because I tested it through REST. I entered the query OBJECTID = 1 through REST and it returned a single record. Any idea what could be wrong? I'm stumped.


title: Cannot add layer Openlayers 3

text:
I have been trying for a couple of days to perform three very simple tasks combining Geoserver, Opencycle maps base map and open layers 3. I am following the pg_routing workshop as I am using a function that let the user click on the map and retrieve the shortest path with a ol.overlay command. I can plot an opencycle map with my data in a different web page but when I try to add the layer in the model of the pg_routing workshop I get the error 'undefined is not a function' immediately after the declaration of the new WMS layer `source: new ol.source.ImageWMS()` I paste the entire web page code here, I have searched extensively everywhere, adding a layer is a very simple operation and I can make in another web pages without the click function. I think it is related to this if I can't easily add a new layer but I can't find where the problem is, can someone help me? thanks.               <!DOCTYPE html>                     <title>ol3 pgRouting client</title>         <meta charset="utf-8">         <link href="ol3/ol.css" rel="stylesheet">

start

final
           <button id="clear">clear</button>


title: Is it possible in uDig to use Features from map to create a profile in profile View?

text: Its great to have the Profile View in uDig. The user can create a profile withe the profile Tool. Its required to create a line in the Map ![enter image description here](http://i.stack.imgur.com/Qj6Zf.png) I was wondering if there is a possibility to select a feature in the map and use the geometry to calculate the profile ![enter image description here](http://i.stack.imgur.com/j0f1Q.png) Is there a chance to to so, How can I solve this requirement? is it possible to switch measurement unit and use the same is used for the distance tool?


title: What is the Scan Direction of a satelite image?

text: There are two possible values that I am aware of for describing the scan Direction of a satellite capture:   1. Forward   2. Reverse     * What is the meaning of Forward and Reverse scan direction?      * What are the pros and cons of each direction?      * Provided I have a have choice should I favor one from the other? Which one?


title: OpenLayers Vector Features with common properties in multiple Layers

text: I am trying to manage features across different layers. Features are displaying weather parameters at some points on the map. Hence there is a wind layer, a temperature layer, etc. which all are vector layers styled differently. One point connects one of the features per layer, shaping something like a "column" of different features: a point connects to a wind feature on the wind layer, a temperature feature on the temperature layer, and so on. Data for the features come from a common GeoJSON file, containing a geographic part and a properties part. It looks like this:                   { "type": "FeatureCollection",        "features": [            { "type":     "Feature",              "geometry": {                   "type": "Point",                   "coordinates" : [-26.142, 48.341]                   },             "properties": {                   "idx" : 0,                   "name" : "Fetone",                   "old" : "626",                   "date" : "2014-06-27 05:50:00",                   "dbid" : "10370210",                   "obsid" : "1",                   "zoom" : "80",                   "type" : "1",                   "temp" : "27",                   "speed" : "6",                   "dir" : "180",                   "hum" : "75",                  "pres" : "1011.00"                  }            },           // ... ... ...           { "type":     "Feature",              "geometry": {                   "type": "Point",                   "coordinates" : [20.125, 62.011]                   },             "properties": {                   "idx" : 200,                   "name" : "Tohature",                   "old" : "1826",                   "date" : "2014-06-27 05:30:00",                   "dbid" : "10367123",                   "obsid" : "512",                   "zoom" : "1",                   "type" : "4",                   "temp" : "19",                   "speed" : "0",                   "dir" : "",                   "hum" : "99",                  "pres" : "1015.00"                  }            }        ]     }      The way that I distribute these data is by essentially copying the properties to all the layers:               var layer1 = new OpenLayers.Layer.Vector("Wind Layer"); // uses "speed" and "dir" properties to style its features     var layer2 = new OpenLayers.Layer.Vector("Temperature Layer"); // uses "temp" property      var layer3 = new OpenLayers.Layer.Vector("Humidity Layer"); // uses "hum"     var layer4 = new OpenLayers.Layer.Vector("Pressure Layer"); // uses "pres"     var layer5 = new OpenLayers.Layer.Vector("Position Layer", { // uses "type", "name" and other properties                         strategies: [                              new OpenLayers.Strategy.Refresh({interval:48000}),                             new OpenLayers.Strategy.BBOX(),                             new OpenLayers.Strategy.Fixed()                              ],                         styleMap: myStyles,                         protocol: new OpenLayers.Protocol.HTTP({                             url: someURL,                             async: true,                             handleResponse: function(r) {                                 var x=r.priv.responseText;                                 if(/^\s*$/.test(x)) { alert("Empty response!"); }                                     var i;                                 var arr=[layer1, layer2, layer3, layer4];                                 var j=new OpenLayers.Format.GeoJSON({'internalProjection': mercator,'externalProjection': wgs84 });                                 for(i in arr) {                                     arr[i].destroyFeatures();                                     arr[i].addFeatures(j.read(x));                                     }                                 // this is needed for selecting the corresponding feature when a select event happens on layer1 or layer5                                 for(i=0;i<layer5.features.length;i++) { layer1.features[i].idx=layer5.features[i].idx; }                                 }                              }),                         renderers: ["SVG","Canvas"], rendererOptions: {zIndexing: true}                         });      However, this seems like taking up too much space, because each layer takes the same copy of the data, and only uses them differently for the style. How can I reduce this amount of data? The next problem is the selection. I try to keep a property "idx" (from "index") with the same number between connected features, which is used to highlight as selected the associated feature whenever I select a feature. Is there something more efficient in this way, handling a collection of connected features between layers? I don't want an arbitrary selection of features on multiple layers, I want the selection of specific pre-arranged "columns" of features. Any help is welcome.


title: How can I correctly delete a join from featurelayer using arcobjects?

text: If I create a join using ArcMap dialog definition query can be changed if joined layer contains a definition query. For example, layer A has definition query "A", layer B has definition query "B", etc if I create join A-B it's definition query will be "(A) AND (B)". if I create join A-B-C-D it's definition query will be "(A) AND (B) AND (C) AND (D)". if I remove join C using ArcMap dialog, the definition query will be restored to (A) AND (B). How can i remove a join using arcobjects and restore definition query?


title: Failed to use InsertCursor to add point features

text: I have following python script to create a feature class in an SDE database and try to add point features to it.               sde_conn = "d:/sde.sde"     st_fc = arcpy.CreateFeatureclass_management(sde_conn, "LIVE_ST", "POINT","","DISABLED","DISABLED",arcpy.SpatialReference(4326))     rowInserter = arcpy.da.InsertCursor(st_fc, ("SHAPE@X", "SHAPE@Y"))          # FOR loop starts here         lon = {a number}         lat = {another number}         rowInserter.insertRow(([lon, lat]))     # FOR loop ends here          del rowInserter      When I ran the script, I received the following error:               rowInserter.insertRow(([lon, lat]))     TypeError: expect single [X,Y[,Z[,M]]]      Are you able to point me to the right direction to get it fixed? Thanks


title: Wrong line segments are being matched by Spatial Join in ArcPy?

text: I'm working on an analysis in which I want to compare 2 polyline-datasets. In one of the first steps, I use a spatial join to create a sub-dataset with only those line segments that have an exact geometric match between the 2 basic datasets.               arcpy.SpatialJoin_analysis(oldlines, newlines, output_commonlines, join_type="KEEP_COMMON", match_option="ARE_IDENTICAL_TO")      If I do this, the titles in the resulting attribute table are :               OID - ID - Attr_a - Attr_b - Attr_c - ID_1 - Attr_a_1 - Attr_b_1 - Attr_c_1      So far so good, the values for almost every coupled segment correspond with the original values, but there is one single exception (that I know of, it is possible that there are more similar cases). The following image indicates two line segments (one straigth and one curved line) , which remain unchanged between the two datasets: ![lines](http://oi61.tinypic.com/x1wx3l.jpg) If I perform the spatial join (which I discussed above), the resulting combined line dataset links the two segments like this: segment1_joined = segment1_old + segment1_new (which is fine). segment2_joined = segment2_old + segment1_new (which is not fine, because 1 and 2 are not identical). Has anyone encountered a similar problem? And if so, did you find a solution? It might also be good to know that segment 1 has no intermediate vertices (so it's perfectly straight).


title: How to control opacity for kmz files in google maps api

text: I need to load several kmz files at one time via Google Maps API, and want to control the opacity for some of the kmz files. Is there a way to do it? I checked http://stackoverflow.com/questions/4651930/using-a-kmllayer-object- within-the-google-maps-api-is-it-possible-to-allow-a-us, but the methods do not work for me. By editing opacity via the css, the whole map's opacity will be changed...


title: Stream Trace to Select Attributes

text: I'm trying to work on a school project, but I only need streams up stream from a certain point. As you can see in the picture, I've run the stream trace tool upstream from a point i've created. I'm trying to find the total length of the red streams. However, you cannot select attributes from a stream trace. How can I accomplish this task? ![Upstream trace results](http://i.stack.imgur.com/of55G.jpg)


title: Polyline frequency to raster

text: I have some animal tracking data in the form of a polyline shapefile and want to convert this to a raster so that each pixel is given a value of 1 for every instance in which a line crosses it. I've tried using Feature to Raster in ArcMap, with the Field specified to a field in which every line has a value of 1, but this only produces a raster with all pixels of "1" value, there is no consideration of the density of lines. I've also used Line Density, but this is taking into account lines within the _radius_ of each pixel rather than only _within_ the pixel. I hope this is clear and appreciate any help.


title: Visualize Data as 3D taken from Excel file on a map

text: My problem is that I have to do 3D visualization of data on map. I have Excel file with three fields: Latitude, Longitude and Value,as shown here Excel file. I need to represent Value field on map as elevated 3D bar or chart with color varying according to value in third field of excel file. It's More like this  3D Visualization of Value. I have checked the technique shown in "Best Practise for visualization of speed data" using Sporttracks Software.It's a third application to be used. I want to do it with some other technique getting same result. Or how can I use the excel file to get the result using Sporttracks Software.


title: Splitting discontiguous vector feature in QGIS

text: I'm a QGIS novice. I'm trying to correct a map of Europe in QGIS, and it's incorrectly assigned Corsica to Germany. I can't seem to split Corsica off of Germany because when I use "split features" to draw a line around Germany or Corsica, it says "No feature split done" because (I think) the line didn't actually cross a feature. How do I split them?


title: How to Edit FileGDB attribute using QGIS?

text: I am able to attach a file GDB in QGIS (version 1.9). But I cann't edit it's content available in attribute table and 'Toggle editing mode' button was also disabled in Attribute table. Does anyone knows how to edit FileGDB attributes using attribute Table ?


title: What Information Does the Blue Band Provide?

text: I am playing around with ERDAS and was implementing image difference on two rasters separated by interval of time. I was using a blue band raster and I created a default threshold for percentage change to be 20% and also did a standard deviation of 2.0. The image gives two completely different outputs for the two types of measuring the difference, using the percentage change and standard deviation. This is the image. The green describes the lower extreme, describing the reduction in change and the red describes the higher extreme, indicating the increase in change. Though I am not able to understand the difference between the standard deviation output and the percentage output and what is the significance of using one over the other? Hope I was able to explain the question well. In the image the upper one is the percentage change and the lower one is the standard deviation.


title: Why can't QGIS open my KML file?

text: I have a kml file with multiple absolutely identical polygons but with different names. I need to load it to qgis and find points in these polygons. But I guess qgis thinks that multiple duplicated polygons is bad and not letting me even load it. Saying that there is an error in the structure. I even found a polygon that he doesn't like and there is no error in it. QGIS can load let's say 3 duplicated polygons but not loading others. Looks like a critical mass for it. So I have two questions. How can I load this file? 2) If I load this file would Point in polygon work for duplicated polygons with different names? I need a system to find points in all duplicated polygons. Thank you!


title: Distance Between Two Points vs Length of Linestring

text: So I have two points that I want to calculate the distance between. If I use:               SELECT       ST_Distance_Spheroid(point_a, point_b, 'SPHEROID["GRS 1980",6378137,298.257222101]') / 1000     FROM (            VALUES (              ST_GeomFromText('POINT(116.330533333333 39.976266666667)'),              ST_GeomFromText('POINT(116.34415 39.97875)')            )          ) AS query (point_a, point_b);      I get the value 1.1953967900687876km. I think this is the correct value since when I use any online distance calculator I get this result; however, if I put the two points into a LINESTRING and then compute the length I get a different number:               SELECT       ST_Length2D_Spheroid(trajectory, 'SPHEROID["GRS 1980",6378137,298.257222101]') / 1000     FROM (            VALUES (              ST_GeomFromEWKT('SRID=4326;LINESTRING(39.976266666667 116.330533333333, 39.97875 116.34415)')            )) AS query(trajectory);      I get the value 1.5228519672706637km. I don't understand why this is happening. In the case of just two points shouldn't these distance bet the same? Thanks in advance.


title: Web Adaptor : Internal vs. External Web Apps

text: We are getting ready to install the **ESRI Web Adaptor** in order to make a select group of Flex maps available externally. The Web Adaptor configuration page says to enter the **GIS Server URL** and the web adaptor will "For GIS Server URL, type the URL to one of the GIS servers in your ArcGIS Server site. This URL will be used to discover all of the machines participating in your site and register them with the Web Adaptor." What if we want one set of Flex maps to be Internal ONLY and one set to be EXTERNAL? What is the configuration required to do that? It appears that if I point the Web Adaptor to my current INTERNAL ONLY GIS Server that everything I put there will be available through the adaptor? Do I have to put my EXTERNAL Flex Maps on the Web Adaptor Server? Or (as I understood it) can I just point the Web Adaptor to a folder and tell it, those apps are EXTERNAL and everything else is INTERNAL ONLY...?


title: Hydrology tool dataset in ArcGIS 10.1 for calculating slope?

text: I am quite new to the hydrology toolset in ArcGIS and thought it may be useful for extracting slope data from a raster; I was wondering if anyone had experience trying to extract slope data upslope from a point? I basically need to find a location 10 meters (on the ground, not as the crow flies) up the steepest slope from a known X,Y location. As a 10 m buffer will not take into account surface geomorphology, I was wondering if anyone had any ideas regarding tools that could be useful? Seems like drainage tools could be handy in this case, but I'm not sure which ones to use or how to create a polyline from drainage raster datasets. This would eventually be put into a python script to automate slope extraction for landslides. Any help/advice is appreciated - thanks!


title: What is the file structure of an ArcGIS Place File?

text: Does anybody know format of .dat file to help with bookmark automation?


title: Customizing popups in ArcGIS Online for GeoRSS services?

text: I'm curious if anybody has successfully customized ArcGIS Online popups from a Geo-RSS source. Specifically, we'd like to remove the stuff in the header, the FID field and also the 'Points' text shown in the layer name (see below). I've read in the documentation that this isn't possible (yet) with AGOL but I'm curious if anybody has found a work-around (perhaps by modifying the CSS) ? ![enter image description here](http://i.stack.imgur.com/3nHtS.jpg) At the moment we're toying with the idea of an ETL process using either ESRI scripted methodology (link), the data interoperability extension for ArcGIS Server, or FME Server. Has anybody else dealt with a similar scenario or can you suggest another option? Thanks in advance for any help you offer. ~j~


title: how to label specific name using maplex

text: I have a water tank data(.shp) and I need to label it using field NAME. However I just want to label water tank that have names like "Tangki Air Sembulan". If in the field NAME only states "Tangki Air", I do not want the label to appear. Im using Arcgis 10.2 and label using Maplex. I try to use SQL Query like below:               NAM NOT Like '%Tangki Air%'       AND NAM NOT Like '%TANGKI AIR%'       but it's not working. I cannot edit the field name so is there a way that can help me solve this problem?


title: Automating polygon mergers with python

text: This is the first real world python program am trying to write. I am trying to create a scrip that takes a shapefile that has n polygons with an integer formatted date (e.g. yyyymmdd or 6/24/2014 is 20140624) attribute and returns a shapefile which has a single polygon for each date, and that polygon is the summation of all polygons with that date or earlier. Additionally I need the polygon to retain the the attributes of the most recent polygon. My initial thought was to manually sort the attribute table form most recent to oldest and have something like this to do the rest:               import arcpy     form arcpy import env     import shapefile          w = shapefile.Writer()     # Create a list of all off the features in the shapefile.     def generate_working_list(a):         env.workspace = "a"         wl = arcpy.ListFeatureClasses()          # copy the oldest polygon, merge it with the next oldest, copy the newly created polygon, repeat      def aggregate(wl):         twl = []         fl = []         wl[-1].append(twl)         wl[-1].append(fl)         wl.pop()         while len(wl) > 1:              wl[-1].append(twl)              '''Dissolve/ merge twl'''              twl[0].append(fl)              wl.pop()         swl[0].append(twl)         arcpy.Append_management(twl, fl) '''I think that this is the right tool to use'''         twl[0].append(fl)         return fl         w.save("Merged")


title: Permanently attaching Image file to Feature Class in 10.1

text: I need to forward on a File Geodatabase, and I am using the attachment manager to attach the file images associated with the different features in my Feature class. If I forward on this FGD, will the images be permanently there? I have not been able to find anything about whether it is permanent, or I would need to send the source files along, and the attachments would have to be redirected on the end of the recipient? If not is there another way to permanently attach an image file to a Feature class in a File Geodatabse? These photos are not geotagged, nor do they have a time stamp on them so neither of the photo tools would work on them. There is no possibility that this can be in SDE. Thanks.


title: QGIS Widget for data entry. Many fields based off 1 combobox

text: I have used QGIS for a while now but have just started getting into creating widgets for data entry. I want to know if there is a way to create a widget that will auto populate 2 or 3 other columns based on 1 combobox. I am doing this for street signs and want the Engineers to enter a point, select the type of sign and that will populate the column an autocad block name as well. For those familiar with ArcMap, you can symbolize points based on unique values, many fields and then start an edit session and create points based on these categories. That's what I want to do with QGIS since not everyone in our office has ArcMap. Thank you.![enter image description here](http://i.stack.imgur.com/iDsQh.png) ![enter image description here](http://i.stack.imgur.com/hggSB.png)


title: Distributing processed SRTM data legally?

text: I'm working on different maps, based on SRTM data. The goal is to lay out the maps at a park entrance against a little money for the printing and my personal costs. The SRTM 4 data from CGIAR I origanally used prohibits the commercial use, so I'm not sure, if I'm allowed to do so. Is there any SRTM data, free to use for that kind of purpose?


title: Map printing using ArcGIS silverlight API

text: I am using a GP task to print my map and its legend. Well it works fine. But it is really slow. Also I am not able to print the route event layer that I create. It does not show on the map. Please suggest me an efficient way of printing the map from client side. I am using ArcGIS 10 Server with Silverlight API on Visual Studio 2010. Thanks in advance


title: How to smooth/generalize a polygon in qgis

text: Is it possible to smooth or generalize a polygon in qgis? Ìs there a plugin available like the generalizer-plugin for lines? I could convert my polygons to lines, smooth them and convert them back to polygones. But this is not working for my polygons because of sveral reasons. Any idea? Thanks!


title: Issue converting DEM in ASCII with LO17 (apparent) projection

text: I'm a newbie where GIS is concerned, so please don't hesitate to point me to a manual. We've got DEM files without any projection information. One is called `000426_LO17_200M_BD16_25.ORT` and contains lines like the following:               -75925.00 2710900.00   834.88     -75925.00 2710800.00   832.98     -75925.00 2710700.00   830.77     -75925.00 2710600.00   832.38     -75925.00 2710500.00   830.83      So we think this is a LO17 projection, at a 200M DEM spacing, and we're guessing the 3rd column is elevation in m. Given that we know the data is from South Africa (or close to it), this seems about right. We're trying to get this data into Google Fusion, but need to convert it into a format that Google Fusion will understand, or equivalently, into a gdal format that we can then convert to Google Fusion. We've tried a gdal_grid conversion, but when attempting to import into Fusion, we get Fusion is complaining about no mapping for TMSO to proj.4 If anybody has any ideas or suggestions, it would be greatly appreciated. All the best, Craig


title: How to create report using arcpy in arcgis 10.2 (result the same as using view -> Reports -> Create Report)

text: ![enter image description here](http://i.stack.imgur.com/RVUGt.png) I need to create reports in arcgis 10.2. because I am not intent to request user to go through all the steps for generating reports, so I want to create an add-in tool to do this. The result is the same as the picture, Can it be realized? Thanks in advance.


title: Arcpy Max Value per Group in One Field

text: I have a point shapefile that have two attributes that I want to extract from them. These two fields are "Group" and "Distance" I want to grab the maximum "distance" from each unique number in the field "group". ![Max Value From Each Group](http://i.stack.imgur.com/ppkUy.jpg) I was wondering if any knew of a method of using arcpy searchcursor to find this max value? Writing something like               rows = arcpy.SearchCursor(input)       Counter = 0       for row in rows:         if row.NEAR_DIST > Counter:           Counter = row.NEAR_DIST      will only give me the max value of the entire dataset as a whole? Any ideas?


title: ArcGIS JS API: Onhover with multiple layers

text: I have a map with two layers: one contains the regions of a country, say, France, as polygon features. And another one contains french cities as point features. I just want to show the **name of the feature on a map tip onHover**. However, can this be done using ArcGIS JS API? When hovering on a city, what the map tip will show? Region name or city name? I know that you can do this onClick as seen here: ArcGIS JS Api: Getting information by location from multiple layers I think if there are overlapping features, onHover only will show the values of the layer at the top. Is this right? Thanks!


title: How to compare string fields in the same table with ArcGIS 10

text: I don't know exactly if there is a tool for this or a maybe a python script could be more suitable. I want to compare string elements from a field (POI_NAME) with another string element in another field (ROTULO). Find attached an screenshot. My goal, AT LEAST, is to have a YES/NO matching, but It could be fantastic having a % of matching. Do anyone know anything about? ![enter image description here](http://i.stack.imgur.com/zYw65.png) I tried something like this: ![Python code](http://i.stack.imgur.com/Xs7ZG.png) But it doesn't work. After editing the code, to: ![New code](http://i.stack.imgur.com/UGiLH.png) We can see that it does not working well (BECAUSE IT HAS SPECIAL CHARS AND SPACES): ![enter image description here](http://i.stack.imgur.com/zsdGT.png)


title: Spatial query in PostgreSQL to find a route to a specific point?

text: I drew 2 shapefiles - point.shp and line.shp. Line.shp has some lines intersecting points in point.shp. I converted the shapefile to PostgreSQL, but how can I make a spatial query to select a line from one point to another point ? Thanks all.


title: How to Turn Standard Web Server into GIS Capable?

text: I have been pointed to help our IT dept to embeded GIS ability in our custom ERP system. The server that host the system has the following standard settings : Apache, PHP, MySQL and JavaScript. The problem is our IT guys not allowed me in "messing" with their system, so to suggesting installing mapserver or geoserver or other gis capable server is highly restricted. Their current solution is converting our native mapinfo tab files into svg files. IMHO this solution will not answered the management needs. I am thinking to use pure openlayers, another problem is the only way to read our tab files is by convert them into geojson. The concept on their mind is all of GIS data is inserted into MySQL, I've heard about MySQL-GIS but I don't know how the plugin is supported today. Basicly I don't do the job, I just give them the right direction, but with this limitation I can't figure out how to achieve this. So to cleared out my question is, I have a standard Apache, PHP, MySQL and JavaScript system, I'm hoping my GIS data is read and write from or to MySQL. How can I turn that standard web system into GIS capable system without changing any setting on my server? OpenLayers are welcome. Any suggestion are most welcome. Regards, Zery


title: How to perform transformations using ProjNet?

text: While writing some unit tests I thought it would be convenient/possible to use ProjNet to convert from EPSG:4326 into several other coordinate systems (3875,3395). I figured it would be as easy as grabbing the WKT values from spatialreference.org:   * http://spatialreference.org/ref/epsg/4326/ogcwkt/   * http://spatialreference.org/ref/epsg/3395/ogcwkt/   * http://spatialreference.org/ref/sr-org/6864/ogcwkt/ Add a reference to ProjNet and running the following code:               var cf = new ProjNet.CoordinateSystems.CoordinateSystemFactory ( );     var f = new ProjNet.CoordinateSystems.Transformations.CoordinateTransformationFactory ( );     var sys4326 = cf.CreateFromWkt ( wkt4326 );     var sys3857 = cf.CreateFromWkt ( wkt3857 );     var sys3395 = cf.CreateFromWkt ( wkt3395 );     var transformTo3875 = f.CreateFromCoordinateSystems ( sys4326, sys3857 );     var transformTo3395 = f.CreateFromCoordinateSystems ( sys4326, sys3395 );      But 3857 and 3395 could not be parsed. I googled around and found a workable 3857 though I haven't confirmed it produces a correct result. I continue to get an error parsing 3395 ( _Expecting ('PROJECTION') but got a 'UNIT' at line 1 column 296._ ). Perhaps the "OGC WKT" can be modified in a predicable way to work with ProjNet? Or perhaps ProjNet is not the best .NET solution for re-projecting coordinates?               const string wkt4326 = "GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.01745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]]";     const string wkt3857 = "PROJCS[\"Popular Visualisation CRS / Mercator\", GEOGCS[\"Popular Visualisation CRS\", DATUM[\"Popular Visualisation Datum\", SPHEROID[\"Popular Visualisation Sphere\", 6378137, 0, AUTHORITY[\"EPSG\",\"7059\"]], TOWGS84[0, 0, 0, 0, 0, 0, 0], AUTHORITY[\"EPSG\",\"6055\"] ], PRIMEM[\"Greenwich\", 0, AUTHORITY[\"EPSG\", \"8901\"]], UNIT[\"degree\", 0.0174532925199433, AUTHORITY[\"EPSG\", \"9102\"]], AXIS[\"E\", EAST], AXIS[\"N\", NORTH], AUTHORITY[\"EPSG\",\"4055\"] ], PROJECTION[\"Mercator\"], PARAMETER[\"False_Easting\", 0], PARAMETER[\"False_Northing\", 0], PARAMETER[\"Central_Meridian\", 0], PARAMETER[\"Latitude_of_origin\", 0], UNIT[\"metre\", 1, AUTHORITY[\"EPSG\", \"9001\"]], AXIS[\"East\", EAST], AXIS[\"North\", NORTH], AUTHORITY[\"EPSG\",\"3785\"]]";     const string wkt3395 = "PROJCS[\"WGS 84 / World Mercator\",GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.01745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],PROJECTION[\"Mercator_1SP\"],PARAMETER[\"central_meridian\",0],PARAMETER[\"scale_factor\",1],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],AUTHORITY[\"EPSG\",\"3395\"],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH]]";


title: How to download only selected attributes using Overpass turbo?

text: does anyone know of a query to download select attributes using Overpass turbo? I am trying to download very granular road dataset for an entire country, but I only want a few columns. Thanks!


title: Getting zone from UTM in order to convert to latlong

text: I am trying to convert UTM to a latlong point. To get that done I need to know the zone where the UTM coordinates are. how can I find the zone from UTM?. the place coordinates that I work on lay on zone 39 and 40.


title: Web site for plan upload and measure

text: I've got plans of my house and a warehouse on paper. I would like to scan them in JPEG, import them on my PC and measure them. I'm looking for a website where I can upload those tiff, jpeg or pdf files, define a scale (measuring something of known measure to give a correct scale to the entire image), measure the rooms by drawing polygons over them, and print out the images with my measured polygons or send them to my email with the total area reported. I can't find anything like this on the web. Is this kind of utility uncommon? Can any body suggest a simple solution for measuring my planimetry without installing any software?


title: openlayers map inside a dojo dialog

text: My application is currently displaying some search results in a dojo grid inside a dialog. I've been asked to embed a map inside the dialog but am struggling to do so. Has anyone done this? I've got a simple map inside a jsp and have tried               dialog.set("href","url_for_my.jsp")      Also tried putting a button on the page to with map.render("mapdiv") after the dialog is displayed. Dialog.refresh() has no effect. There are no errors and no requests made for map tiles. Can anyone suggest a better way of displaying a map in a popup? Thanks!


title: How to convert WMS layers to raster layers QGIS 1.8.0

text: I need to export a series of maps for a report. The layers imported into QGIS are all WMS files provided by the GeoBC website. As I have now discovered, I cannot export these maps using the print composer because of the limitations of WMS. I also cannot add a scale bar to these maps as the WMS CRS is WGS 84 not UTM. If I understand the problem correctly, the solution would be to convert the WMS layers in to raster layers initially, and then to vector layers? Does this sound correct? If so, is there a tool to convert WMS to raster I can use as the first step?


title: How to derive point height from contour data?

text: I got vector layer that represents height over sea level. Then, I also have different vector layers. Is there any good way of using QGiS or things related to calculate sea level for those points? Edit:   Those points that are in different layers only have location. They are names of hills, rivers and so forth. What I need is to merge the info of those two layers. So I know how high over sea level the items are. Edit two; (Sorry not subject that I am strong in) I got line vector layer. With line every 20m of elevation. Over all of Iceland. With about 60.000 points that I need to find hight over sea level on. I dont mind leting the computer work on this but not by hand :)


title: Merge overlapping DEMs with irregular boundaries

text: I have a 10m DEM covering part of a state, a 20m DEM covering the whole state, and a 90m DEM covering the whole country. How can I merge them together, getting the best quality available in each area? The 10m and 20m are provided in ADF format, but I've converted them to .tif using ogr2ogr. Using gdal_merge doesn't work, because the boundaries of each area are irregular. This is the result: ![enter image description here](http://i.stack.imgur.com/v5U5h.png) (The northeastern edge of the state is irregular, but there obviously shouldn't be the gap in the middle.) **UPDATE** If it's useful, I also have a shapefile for each 10m DEM, marking its extent. **UPDATE 2** Ok, here's what the whole state looks like with the 20m DEM. ![enter image description here](http://i.stack.imgur.com/HKZiL.png) The process I'm going through so far is this:               for f in dtm20m dtm10m_nw dtm10m_e; do     export GDAL_CACHEMAX=1000     echo -n "Re-projecting: "     gdalwarp  -s_srs EPSG:3111 -t_srs EPSG:3785 -r bilinear ./vmelev_${f}/${f} ${f}-3785.tif           echo -n "Generating hill shading: "     gdaldem hillshade -z 5 $f-3785.tif $f-3785-hs.tif     echo and overviews:     gdaladdo -r average $f-3785-hs.tif 2 4 8 16          done


title: Cost Attributes- Network analyst for Arc GIS 10.2

text: Hi can anyone tell me whether Network analyst can consider multiple cost attributes? I am looking for most ecologic route(not shortest), therefore I have assigned a value of pollution to each road segment. I have 2 layers in my analysis road and cycle layer, both have these pollution values for each segment (cycle routes have mostly lower values) and I selected this pollution attribute as a default cost attribute but would also like to consider another attribute in the best route analysis. Can this be done or does network analyst always consider just one cost attribute?


title: Hex grid value interpolation in POSTGIS

text: I have a table of tessellated hexagons that I created as a hex binning grid to cover an area. I've then overlaid a point set to give some of the hexagons a value, see diagram:- ![enter image description here](http://i.stack.imgur.com/rvzTI.png) The shapes with value "0" have no value and I'd like to interpolate these from the values of the neighbours as per Kriging. How can I accomplish this in postgres given my table looks like:-               hex_grid_data     (      gid SERIAL,      wkb_geometry (geometry),      value INTEGER     )


title: How do I export only additions when running export-diff to a shapefile in geogit (geogig)?

text: I currently have two commits in my GeoGit (GeoGig) repository. I'm running the following command on the two commits to get the differences exported to a shapefile. `geogit shp export-diff [commit1] [commit2] [path] differences.shp` I know doing the `diff` command that there are X number of Additions, X modifications, and X removals. The shapefile creating is providing me both additions, removals, and modifications, but what I would like is a shapefile showing only the additions, a shapefile showing only the removals, and a shapefile showing the modifications, or possibly add a column to the shapefile showing what the status of the feature is (A, R, M). Is this even possible with GeoGit (GeoGig) currently, and if so what arguments should I add to this statement?


title: Lat Long to Useable values in QGIS

text: I am using both QGIS 1.8 and 2.0. So either is helpful. A client has given me locations of trial pits but i cannot load them into QGIS no matter how hard i try. For example, 50°24.711 North 7°25.529 West. The site is near Enniskillen Co. Fermanagh Northern Ireland. Any system i can find to convert degrees minutes and seconds will not allow 24.711, and seconds only go to 60. Any help is greatly appreciated.


title: Change Map layers when zooming in (JavaScript API)

text: I have created a web map application (Javascript API) :               require([       "esri/map",       "esri/dijit/InfoWindowLite",       "esri/InfoTemplate",       "esri/layers/FeatureLayer",       "dojo/dom-construct",       "dojo/domReady!"],     function(       Map,       InfoWindowLite,       InfoTemplate,       FeatureLayer,       domConstruct     ){       map = new Map("map", {         basemap: "topo",         center: [-85.416, 49.781],         zoom: 5,         logo: false});            TestLayer = new esri.layers.ArcGISDynamicMapServiceLayer(         "http://Myserver/arcgis/rest/services/TestMAP/MapServer");            query = new esri.tasks.Query();       query.returnGeometry = true;       query.outFields = ["CCSNAME", "DAUID20111" ,"DAUID201_5" ,"DAUID201_9" , "DAUID20_14",  "DAUID20_19", "NAME" , "CITY",  "POST_CODE", "HRUID2007.csv.GPF", "HRUID2007.csv.GPE", "HRUID2007.csv.SPECF", "HRUID2007.csv.SPECE" ];            symbol = new esri.symbol.SimpleFillSymbol(         esri.symbol.SimpleFillSymbol.STYLE_SOLID, new esri.symbol.SimpleLineSymbol(           esri.symbol.SimpleLineSymbol.STYLE_DASHDOT, new dojo.Color([0,255,0]), 2),new dojo.Color([0,255,0,0]));            map.addLayer(TestLayer);      In the default zoom level set, 3 map layers from `**http://Myserver/arcgis/rest/services/TestMAP/MapServer**` are showing up. How can I change my code so that in the default zoom level I get another 3 layers from this url:               "http://Myserver/arcgis/rest/services/TestMAP2/MapServer"      and then when I zoom in I get the current ones.


title: Set the parameter values retreved from database in Arcgis script tool in toolbox using python

text: I am now thinking a function to allow user to modify his information stored in the database. To do this, I developed 2 script tools in ArcGIS(10.2), One is the ask user to input his ID, the other one is to retrieve the data from database and allow user update the information before update it in database. The 2 tools are shown below in sequence. ![enter image description here](http://i.stack.imgur.com/RBHcs.png) ![enter image description here](http://i.stack.imgur.com/0H1Ed.jpg) I have 3 questions:   1. How to set the parameter values in the script tool in toolbox, I mean when open the tool, all the existing value retrieved from database has already displayed in the tool. I tried to search the solution but haven’t got a useful solution.   2. When I tried to use the first tool to call the second tool, is it possible to close the first one? It may look messy to the user.   3. Since the user ID is unique and cannot be modified, how to prevent this parameter from editing.  Thanks for your attention and advice in advance. I am so sorry the picture not that clear.


title: GPS waypoint mapping with android phone

text: I am a bat enthusiast and volunteer regularly with the local bat groups on bat surveys. In them we record bats with a bat detector - making a note of the different species (there are 18 in the UK), the activity and its location on a map (which can be inaccurate). Once home I will then take this information and on word put different coloured dots (depending on the species) on a digital map to be sent away to be filed. I was wondering if there a way to do the recording part on my android phone. I've downloaded various waypoint recording app which give me a gpx or kml file. With this i need a mapping system that can take this information and lay it onto a chosen map. also, the program needs to be able to distinguish between the various colours (species) and these colours need to show up on the map. Also, if there are more than one species recorded on the same location the dot needs to be twice its size. so to reiterate the program needs to give me flexibility to modify what the dots look like. hope someone can help. thanks.


title: QGIS 2.0 portable or silent installation

text: Is it possible to install QGIS 2.0 unattended (silent) on Windows-Client (XP,7) via commandline ? Or does there exists a portable version for 2.0 ?


title: Shortest path to nearest... something

text: I have a list of addresses and a map including train stations, schools etc., and I would like to calculate the **_shortest_** distance from each address to the **_nearest_** train station, school etc. I'm aware that the `google maps API` can calculate the shortest distance to any specific point, but my problem is complicated by the fact that I also need to find the nearest train station etc. Is there an easy way to do this in QGIS?   * I use QGIS    * Windows 7 64bit Thanks!


title: Postgis | Problems Reprojecting a Point from UTM ED50 to GEOGRAPHIC ETRS89

text: I've got a POINTS PostGIS table in the 23030 (ED50 UTM 30N) SRID. I'm having a lot of problems trying to reproject those points from ED50 UTM coordinate system to ETRS89 Geographic (lat/long) coordinate system. The aim is to convert the points for Google Maps. I take an example point for my tests:               x: 499462.1     y: 4799340.56       1.- If i just transform the point from UTM (23030) to Geographic lat/long (4230) using the `ST_TRANSFORM("GEOMETRY",4230)` PostGIS function. I get the following coordinates:               long: -3.00663639334493     lat: 43.3461345057409      You can see in the picture where it is located in Google Maps. It has sense to be located in a different place to the original point, because although being in Geographic projection system it is not in ETRS89/WGS84. 2.- So I reproject the original point first to ETRS89 (nearly the same as Google's WGS84) using `ST_TRANSFORM("GEOMETRY",25830)` and I get the following point:               x=499462.124763252     y=4799253.63249274      Ok, if I load that point in a GIS, I can see there is a difference between the previous point which is ok as I have changed the SRID. 3.-I thought that if I reprojected the ETRS89 UTM point to ETRS89 lat/long it would be enough for getting the coords valid for Google Maps. The surprise came when I got the same point of the step1!!! (just doing the 4230 transformation) I used a `ST_TRANSFORM(ST_TRANSFORM("GEOMETRY",25830),4258)` UPDATE query applied to the original point. ![https://www.dropbox.com/s/6r82qtimboxnf9z/ejemplo_foros.png](http://i.stack.imgur.com/SRGEf.jpg) I don't really know how to solve this mess, I'm really confused.


title: How to convert all values in all fields of table to upper case using ArcPy & Python?

text: Does anyone know how to use the .upper() in a stand alone python script so it can be used to change all characters in a table into UPPER case? I don't want to use the field calculator since I have multiple tables, etc... I currently am using a cursor to correct mistakes in my tables, but can't figure out the correct syntax to use it to convert all the strings to UPPER case!


title: arcpy.mapping.Export creates some labels but not others

text: I've created a simple python toolbox (pyt) to export map layouts to multiple formats. My problem is that some labels which are defined in the layout are exported, but others are not. For example point labels export but contour labels dont. When I export from within ArcMap all labels are exported. Question: Is there any argument missing from my "arcpy.mapping.Export" statements that would affect labeling? I had problems posting the whole toolbox but it can be downloaded from here: Export MXD Toolbox Here's the code that does the exporting.               # mxd name     mxd = arcpy.mapping.MapDocument(parameters[1].valueAsText)     # workspace     env.workspace = parameters[0].valueAsText     # workspace text to use in paths     workSpace = parameters[0].valueAsText     AddMessage("The workspace is " + env.workspace)          # resolution     dpi = parameters[2].valueAsText     # filename     name = parameters[3].valueAsText     # option for formats to export     createPDF = parameters[4].valueAsText     createTIFF = parameters[5].valueAsText     createEPS = parameters[6].valueAsText     createAi = parameters[7].valueAsText          # Export PDF     if createPDF == "true":       AddMessage("Exporting PDF")       arcpy.mapping.ExportToPDF(map_document = mxd, out_pdf = workSpace + "\\" + name + ".pdf", resolution=dpi,                convert_markers=True, layers_attributes="LAYERS_AND_ATTRIBUTES")     else:       AddMessage("You chose not to create a PDF")          # Export TIFF     if createTIFF == "true":       AddMessage("Exporting TIFF")       arcpy.mapping.ExportToTIFF(map_document = mxd, out_tiff = workSpace + "\\" + name + ".tif", resolution=dpi, tiff_compression="LZW")     else:       AddMessage("You chose not to create a TIFF")          # Export EPS     if createEPS == "true":       AddMessage("Exporting EPS")       arcpy.mapping.ExportToEPS(map_document = mxd, out_eps = workSpace + "\\" + name + ".eps", resolution=dpi, image_quality="BEST",                picture_symbol="VECTORIZE_BITMAP", convert_markers=True)     else:       AddMessage("You chose not to create an EPS")          # Export to Ai     if createAi == "true":       AddMessage("Exporting Ai")       arcpy.mapping.ExportToAI(map_document = mxd, out_ai = workSpace + "\\" + name + ".ai", resolution=dpi, picture_symbol="VECTORIZE_BITMAP",                convert_markers=True)     else:       AddMessage("You chose not to create an Ai")


title: How to add an interpolation raster from the python console in QGIS?

text: I'm trying to add an interpolation layer based on an existing vector layer by using the Python Console (initially in the menu Raster/Interpolation), but can't find anything related in the documentation. I saw the QgsIDWInterpolator in the API, but can't get it working. It seems it expects a list of LayerData but it doesn't work neither. Here my current code :               ld=qgis.analysis.QgsIDWInterpolator.LayerData     ld.vectorLayer=vlayer     ld.zCoordInterpolation=0     ld.interpolationAttribute=5     ld.mInputType=0     ldlist=[ld]     itp=qgis.analysis.QgsIDWInterpolator(ldlist)      I get the following error for the last command :               Traceback (most recent call last):       File "<input>", line 1, in <module>     TypeError: arguments did not match any overloaded call:       QgsIDWInterpolator(list-of-QgsInterpolator.LayerData): argument 1 has unexpected type     'list'       QgsIDWInterpolator(QgsIDWInterpolator): argument 1 has unexpected type 'list'      Thanks !


title: How to create mask file from shapefile using ArcGIS for Desktop?

text: I have a shapefile of Lakes over my study area. Each polygon has attribute information of each lake, name, shape etc. Actually I want to eliminate all lakes from analysis. In other word, I want lakes to have 0 value and the area that out of boundary of lake to have their orginal pixel values. when I do subset using this shapefile, the output is not I want because my shapefile deoesn't have any value outer boundary of lake. So how can I built a proper mask file (inner part of lake has 0 value, outer is 1) using this shapefile in Arcgis?


title: Postgis 2.0 geographic polygon has max longitudinal range

text: My problem is that I have a load of geographic point data (srid 4326) in a table. A user may do a bounding box search - essentially create a polygon as follows:               l_poly geography(Polygon,4326):= 'POLYGON((-120 -89, -120 81, 59.999 81, 59.999 -20, -120 -20))';      ... then query using the ST_Intersect predicate               WHERE ST_Intersects(my_tab.point_data, l_poly)      In the above example, I can return a point at lon=0, lat=0. However, once my eastern-most longitude goes above 60degrees (currently 59.999), then the longitudinal range exceeds 180 degrees, and my polygon wraps the opposite way around the earth. In this case I can not return data at 0 longitude, but I would at e.g. 180 longitude. It seems as though a geographic polygon cannot be defined that exceeds a longitudinal extent of 180 degrees. Is this true? Can anyone point me in the direction of a solution? I have tried casting a geography from a geometry defined with enforced Right Handedness, which I thought might define the interior of the polygon to the right hand side of a clockwise drawn polygon. But this doesn't work either. Thanks in anticipation. :)


title: GeoWebCache: Layer is not served by a WMS backend

text: I try to get legends for my layers in a little GIS application while using GeoWebCache. Since I turned GeoWebCache in GeoServer on I do not get the legend image for my layers anymore. I looked around and found this thread here: How to get WMS LegendGraphics using Geoserver and GeoWebCache?. However, my attempt to set the 'legendURL' by using '...LAYER" + map.layers[6].params["LAYERS"])' fails. The only thing I am getting back is that 'undefined' is not an object (evaluating 'mapPanel.layers[6].params'). Can somebody help me here? The application is temporarily here available. Thanks a lot, Thomas


title: How to resolve "premature end of script headers" error with QGIS Server?

text: I attempted installing QGIS Server following these instructions. Installation was successful with no errors reported. However, when I try to use `GetCapabilities` I am getting this error: > **Internal Server Error** > > The server encountered an internal error or misconfiguration and was unable > to complete your request. > > Please contact the server administrator, `admin@local.domain` and inform > them of the time the error occurred, and anything you might have done that > may have caused the error. > > More information about this error may be available in the server error log I looked for the error in Apache's logs and found this entry: `[Mon Oct 28 16:02:40 2013] [error] [client 127.0.0.1] Premature end of script headers: qgis_mapserv.fcgi.exe` Can anyone suggest how to troubleshoot this issue?


title: Looking for a MAP to SLD converter

text: Does anyone know where I can find the converter from .map to .sld? I have a big mapfile of styles for OSM and I have the styles mapnik.xml for OSM but I use Geoserver for my work. Or does anyone have some sld styles of OSM for PostGIS layers planet_osm_point and planet_osm_polygon? (I already have the style for planet_osm_line.)


title: How to set default values for reference grid in ArcMap?

text: I currently use ArcMap 10.2 and find it quite handy to be able to set default values for settings such as Labeling Properties, ensuring new maps use relative paths and to always build pyramids for a raster in the **Customize** --> **ArcMap Option** window. Is it possible to set the default values for a new reference grid in ArcMap? Specifically, is it possible to set the initial default value for the number of decimal places for a new reference grid, instead of clicking through several windows to set number from 6 decimal places to 0 every time a new reference grid is created? I ask because the majority of the maps I make have no use for any decimal places in the reference grid.


title: GRASS 6.5 or 7.0 - Need Help with Loop Script

text: I am new to Python and scripting in general. This is also my first time a posting question. My overall goal is to create a script that I can cycle through all the raster maps in a mapset and then save the univariate stats to a .csv file for each raster map. I am working on the first part: cycling through all the maps and printing the stats to the terminal window. Later, I will try and save the data to a .csv file. I have seven maps in my practice mapset and I get seven errors, but my problem is that I do not know what variable to assign to the parameter MAP (hopefully, I stated that correctly). Here is what I have thus far:               import sys     import os     import grass.script as grass     import grass.script.setup as gsetup     gisbase = os.environ['GISBASE']     gisdb="/Users/kc/GrassData"     location="Data_2013"     mapset="methods_check"     gsetup.init(gisbase, gisdb, location, mapset)     grass.run_command("g.list", _type="rast")      mymaps = grass.parse_command("g.list", _type="rast")      print len(mymaps)      for items in mymaps:         grass.run_command('r.univar', map='items', separator=",")      What should I put for map = 'items' so that it cycles through all the raster maps of a mapset? I am using GRASS 7.0, but the same thing happens in GRASS 6.5. In Grass 6.5, the parameter separator must be changed to fs. Both versions label the input file "map". Any help would be greatly appreciated. Thank you.


title: Maximum number of record returned by the server. Why 1000?

text: Why in AGS "maximum number of record returned by the server" is set to default 1000? What happens if I increase limit to 100000 or 1000000 in 100 services?


title: Eliminate 0 Values From Raster

text: Arcgis 10.2 I currently have a raster file that I want to turn into a set of polygons. This will help me to create specific buffers around certain points. The extent of my values are 0, and 1. When I change from raster into polygon, I end up with polygons for both 0 AND 1. I only want polygons for pixels with value 1. If I have polygons for values of 0, I end up capturing a bunch of empty space that it totally not necessary to the project. Here is an image of what I am working on (white = 1, black = 0): ![enter image description here](http://i.stack.imgur.com/eJXm5.jpg)


title: Assign Custom Band Combination To Image Catalogs In ArcGIS

text: I have an image catalog composed of 4 band rasters. I would like to use the same image catalog, and set custom renderers for a true color and false color combination and save each of those catalogs off as a layer file. In addition I would like to turn off the default stretch for each catalog. The code I am using is below (where lyr is a reference to a raster catalog):                               IRasterCatalogLayer rcl = lyr as IRasterCatalogLayer;                          IRasterRGBRenderer pRGB = new RasterRGBRendererClass();                     IRasterStretch pRasStretch = pRGB as IRasterStretch;                     rcl.Renderer = pRGB as IRasterRenderer;                          pRasStretch.StretchType = esriRasterStretchTypesEnum.esriRasterStretch_NONE;                     int r, g, b;                     pRGB.QueryBandIndices(out r, out g, out b);                     pRGB.UseRedBand = true;                     pRGB.RedBandIndex = 3;                     pRGB.UseGreenBand = true;                     pRGB.GreenBandIndex = 0;                     pRGB.UseBlueBand = true;                     pRGB.BlueBandIndex = 1;                          int r1, g1, b1;                     pRGB.QueryBandIndices(out r1, out g1, out b1);                          rcl.Renderer.Update();      The stretch is being applied without problem, but I cannot change the band combinations. Thoughts? Is there a better way to do this without resorting to ArcObjects to create layer files.


title: GeoExt Vs Openlayers

text: OK: [ Openlayers + ExtJS = GeoExt] I am new to webmapping and still shopping around. I bumped into a Case Study with a gentle learning curve, matches my needs and it uses OpenGeoSuite. I got to learn of GeoExplorer and I found it to be a neat interface as I intend to serve geodata to gis novices-- some eye candy would be necessary. I however have issues with the feel of the 'Identify' pop-up box. My questions: Is it wise to adopt this 'interface'? (i) does it have a future - future releases/development (need 2nd opinion after reading http://opengeo.org/technology/geoext/ ) (ii) easy of customisation -- little programming/ scripting experience [ remember 'Identify'] (iii) shall I instead focus on mastering OpenLayers as it is?


title: Organizing PostGIS data of a large database for fast dynamic user searching

text: I have a large amount of GIS data, roughly 30 layers, each with 50,000 to 100,000 (Geo)records. Some of this data is GPS collected points, others are POI, some is government provided data (roads, rivers, water sheds, etc), So far I'll either be using mapfish, geoserver, or feautureserver with jquery/ajax backend. How should I organize my GIS Data to allow for fast and dynamic webside client searching? as well as feel free to provide your experience from projects especially working on large data sets with heavy user quering


title: Open layers plugin problem

text: I am working in open layers plugin.i need to project the satellite data with respect to vector data shape file i am having .for that Andre Joost suggested to create vector grid with vector data reference system .The new problem is when i change the coordinate system to degree minutes seconds the corrdinate system does not changes and it shows only the cordinates in meters itself .what is the mistake i have done?


title: Model Builder not completely honoring lyr file

text: I have created a series of map templates for one of our departments that each have an associated model. As part of the model I also created a lyr file for symbology and labeling (the subject parcel layer portion). ![enter image description here](http://i.stack.imgur.com/1k1EG.jpg) This is referenced in the output feature parameters of my model. ![enter image description here](http://i.stack.imgur.com/GzE5O.jpg) The tool runs correctly and the Subject Parcel_layer gets added to the map template. The overall symbology from the layer file is correct, however, the labeling preferences I've set up are not (font size, type, etc.) As a side note, the Subject Parcel Layer does not exist in the template, it gets added through the model. Any ideas on why the labeling portion is not being honored?


title: QGIS SimpleSymbolRenderer label for Symbol in Legend

text: Is there an ability to set a label for the symbol (in QgsSingleSymbolRendererV2) which is showing in the legend?


title: change symbology of an graduated renderer in qgis

text: In my plugin(c++) I add a vectorlayer to map, then I change the symbology of the QgsGraduatedSymbolRendererV2 of this layer. After this I refresh the map with `mQGisIface->mapCanvas()->refresh();` In legend the symbols were all correct shown, but not in map. Just after press F5 or click the refresh-button the symbology will be right rendered in map. What is the different between the map-refresh button and the command in code? what could I do to refresh in code? Thx for any hint.


title: Extracting area with a range of lat long?

text: Is there any online tool to extract a certain area with a range of lat long? I'd like to overlay some locations on a map. To do so, I need a map whose range is specified so that I can determine the locations. I found openStreetMap does it but it's got some extra information on a map, including streets, shops and so on. For now, I need only country boundaries and names of country in Europe. Is this that hard to do this simple thing? Could you help?


title: Create lines from points arcpy

text: I have a points shapefile, and a file that tells me which of the points have a line between them. I am using arcpy to create a new shapefile to hold the lines and then create the lines based on the specifications in the file. For creating the lines, I am reading the coords of the points from the point shapefile itself.               out_path = arcpy.env.workspace     out_name = arcpy.ValidateTableName("the_lines")     geometry_type = "POLYLINE"     has_m = "DISABLED"     has_z = "DISABLED"     spatial_reference = arcpy.Describe(points_shp).spatialReference # Use Describe to get a SpatialReference object     # Execute CreateFeatureclass     arcpy.CreateFeatureclass_management(out_path, out_name, geometry_type,spatial_reference = spatial_reference)          array = arcpy.Array()     point = arcpy.Point()     tie_shp_file = arcpy.env.workspace + "/" + out_name     insert_cursor = arcpy.InsertCursor(tie_shp_file)     feat = insert_cursor.newRow()          for line in all_lines:   #all_lines is of the form array([[0, 3],[0, 7],[1, 3]) indicating line [from, to]         #from         point.x=coords_pts[str(pt[line[0]])][0]         point.y=coords_pts[str(pt[line[0]])][1]         print "FROM: " + str(pt[line[0]]) + " COORDS_X: " + str(coords_pts[str(pt[line[0]])][0]) + " COORDS_Y: " + str(coords_pts[str(pt[line[0]])][1])         array.add(point)         #to         point.x=coords_pts[str(pt[line[1]])][0]         point.y=coords_pts[str(pt[line[1]])][1]         print "TO: " + str(pt[line[1]]) + " COORDS_X: " + str(coords_pts[str(pt[line[1]])][0]) + " COORDS_Y: " + str(coords_pts[str(pt[line[1]])][1])         array.add(point)         print "ARRAY: " + str(array)         polyline=arcpy.Polyline(array)# Create a Polyline object based on the array of points         array.removeAll() #Clear the array         feat.shape = polyline         insert_cursor.insertRow(feat)      coords_points is a dictionary that has the form:                {'10': (8594242.514466941, 3330642.8703512326)...}     i.e. {'Oid_str': (x,y)}      The code is creating the required number of entries in the shapefile but no lines are being shown and the lengths of the line is shown as 0. Where am I going wrong?


title: Intersect and MBR in POSTGIS and Mysql

text: I would like to know if POSTGIS and MySql with Spatial Extension use true geometry to find intersections/overlaps between geospatial objects. But i don't find it precisely in the documentation. I know that doing an ST_Intersects is more precise than only using the && command. I want to know if in this example the round geospatial objects would intersect/overlap: https://dl.dropboxusercontent.com/u/13064345/example.jpg The squares are the minimum bounding boxes, they overlap. But i am not sure of the ouput in POSTGIS and MySQl. Hope you can help me.


title: Exporting Maps as PDFs in R

text: In R, trying to export choropleth maps with density and frequency histograms as legends; have had success exporting as an image but not as a pdf. Exporting as an image (.png) works well; the map, including the placement of the histogram within the frame, looks the same in the exported image as it does within the Plot window in R. However, when exporting to a PDF, the histogram changes size and placement and often overlaps with the map (which renders the histogram unreadable). Any ideas as to why the two exported maps look different? And how to export to a PDF so that the exported map looks like the map in R? Thanks!


title: Model Builder Error populating attributes with Calculate Field

text: I would like to run a model that iterates over all feature classes in a feature dataset and uses the feature class name to populate an attribute field (see model below). However, when running the model I get the below error: ![enter image description here](http://i.stack.imgur.com/Q5vW2.jpg) Might the problem be in the calculate field expression? It gives me the option to use "Name" from the iterator (hence the arrow connecting name to the calculate field tool). But I can't get the model to run. The resource cited here has not provided much insight. ![enter image description here](http://i.stack.imgur.com/tFilT.jpg) thanks in advance I also tried to split the models in two (see below). The first one executes no problem. The second gives me the same previous error above. ![Adds text field successfully](http://i.stack.imgur.com/frbDR.png) ![enter image description here](http://i.stack.imgur.com/zc35K.jpg)