title: Calculation based on features in shapefile

text: I hope I get the terminilogy right! Correct me if I mess something up. I need to perform some calculations for a scoring system for different habitat types. The Problem is that I have two different tables in two different shapefiles. They need to be divided through each others somehow (yet again, all needs to be done automatically in the modelbuilder). For the calculation I need to Sum up the `Size` of all the `Types` with a common/matching `Name` in `Table A` (I tried to clarify this through the field `Sum`).               Shapefile A with Table A     Name    Type    Size   Sum     A       BO      10     ?     A       BE      20     ?     C       BA      30     ?     D       BO      10     ?     D       FW      50     ?      The `Sum` for `Name A` would be `30 (10+20)`, for `D` it would be `60 (50+10). In my second table, I do need not to do any calculations. It looks likes this:               Shapefile B with Table B     Name    Type    Size     A       BO      20     A       BE      50     C       BA      60     D       BO      90     D       FW      10      I know the `Size` of each `Type` in in each `Name`. The `Size` of `Table B` needs to be divided through the `Sum` of all the `Types` for each `Name` in `Table A`. I hope you can see the problem I am running in. How do I tell the Modelbuilder to summarize all the `Size` of all `Types` that have a matching `Name`? And how would I populate `Fields` with this outcome for further calculations?


title: Pgrouting can't find postgresql (on Mac)

text: I've installed both postgreSQL 9.0 and postgreSQL 9.1 on my mac, using binaries from http://www.kyngchaos.com. I also have postgis 1.5 installed. I tried to install pgrouting, again using binaries from aforementioned website, but when I ran the core sql it gave the Error: could not access file “$libdir/librouting (How to fix Error: could not access file "$libdir/librouting"?) error. So, instead I tried to build it from the source. but now when I run: cmake . , I get" -- PostgreSQL not found. FATALERROR Please check your PostgreSQL installation" I tried to append it to path using "export PATH="$PATH:/usr/local/pgsql-9.0/bin"" but it still didn't work. Any help is appreciated! Thanks P.s. I should mention that I want to tell the system to install pgrouting for postgresql 9.0 and NOT 9.1 version.


title: Is MapMaker interoperable?

text: Has anyone used / come across MapMaker GIS Software http://www.mapmaker.com? A consultant has stated that they will use this for a piece of work involving GIS analysis and that it is compatible with ArcGIS. The work will involve analysing many datasets, predominantly in SHP format to produce a land management strategy - the outputs of which will need to feedback into our corporate system (PostGIS server / ArcMap desktop). From the website it seems that it will import SHP files, but no mention of export. It looks like a cartographic tool rather than capable of GIS analysis. I would be grateful if anyone has used it and can comment on suitability for task and interoperability. Thanks


title: How do I use QGIS python console to iterate through simple layer queries and png outputs?

text: I'm really new to QGIS coding and I'm having a hard time finding (especially) the api classes/functions to carry out actions of the Query Builder window. I'm basically trying to do this:   1. Iterate from n=1 to 200...   2. Query layerA where field_id = n   3. Query layerB where field_id = n   4. Query layerC where field_id = n   5. Save as img_n.png Pretty simple, I think, but it has become a real pain to do by hand in the GUI.


title: How to label a WMS layer from an Online Server, over which I do not have control?

text: I have discovered that the Western Australia Dept of Mines and Petroleum offers WMSs. On adding these to ArcMap, however, I notice that they do not have labels. For example, this means one must use the identify tool to get the name of a petroleum permit under application. This is a bit dismaying, as I was hoping to use these WMSs in place of downloading the data manually once a month. Without labels in the WMS tiles, I find them less useful for making maps or putting them in ArcReader. Is there a way to force the labels to appear?


title: problem in using schema.ini to force ArcMap recognize field types

text: I have a folder of csv files (500+), and created a very long schema.ini file to force ArcMap recognize field types as I want. However, for some reason, ArcMap can read some CSV files properly; but cannot read the rest correctly. The problem with the rest CSV files is the fields are not recognized, although schema.ini has the needed information. Is there a length limit issue in schema.ini? I guess it may not work well when the schema list is getting long.


title: Adding a boolean yes/no field in Arcgis

text: How does one add a field to an Arcgis feature class with a boolean data type? That is an attribute where the allowed value is only one of a pair such as 1/0, Yes/No, On/Off, Present/Not-present, etc. A boolean data type is not available in the Data Type pick-list.


title: How to render points from CartoDB with Leaflet and Leaflet Vector Layers?

text: I have a table in CartoDB containing some points. I want to render these points dynamically with Leaflet. There is cartodb.js for such tasks, but it does not support full customization of Tooltips, as I have seen. Instead I found a Leaflet-Plugin for that. In my HTML I have the following JS-Code, but it does not work except showing the Base Map. What am I doing wrong? Are there other solutions for Leaflet?               var map = new L.Map('map').setView([51.0, 10.0], 6);     var defaultLayer = L.TileLayer.provider('MapQuestOpen.OSM');     map.addLayer(defaultLayer);          var baseLayers = {     "Mapquest OSM": defaultLayer,     "OpenStreetMap German Style": L.TileLayer.provider('OpenStreetMap.DE'),         };          // Use Leaflet Vector Layers plugin to show data from CartoDB          var cartodblayer = new lvector.CartoDB({         user: "user",         table: "example_data",         scaleRange: [6, 20],         symbology: {             vectorOptions: {                 icon: new customIcon({                     iconUrl: "/icon/icon.png"                 })             }         },         popupTemplate: '<p>{homepage}</p>',         singlePopup: false     });          cartodblayer.setMap(map);          map.addControl(new L.Control.Layers(baseLayers, '', {collapsed: true}));


title: Converting .kml file to .tab (Map info) format

text: If I have saved polygons as a (.kml) Google Earth File and I want to draw it in MapInfo, is there a tool to do this ? If not, how can be it be made using MapBasic?


title: PROJ4 / PostGIS transformations between WGS84 and NAD83 transformations in Alaska

text: I'm trying to implement the recommended WGS84 - NAD83 transformation for Alaska using PostGIS. This post: Geographic Transformations in Alaska Between NAD83 and WGS84 indicates the proper transformation in ArcGIS to be WGS_1984_(ITRF00)_To_NAD_1983. I haven't been successful at finding this in `spatial_ref_sys`. I tried adding my own transformation into `spatial_ref_sys` using the following proj4text:               +proj=longlat +ellps=GRS80 +datum=NAD83     +towgs84=-0.9956,1.9013,0.5215,-0.025915,-0.009246,-0.011599,-0.00062     +no_defs      The `towgs84` values were copied from the listed settings that ArcMap shows when you choose this transformation (reversing the sign because it's to WGS84 in PROJ4 and from WGS84 in ArcMap). In PostGIS, I use it like this, transforming WGS84 to NAD83 using my custom projection, then transforming it again into 4269 (or whatever NAD83 projection I really want, like 3338).               ST_Transform(ST_Transform(geom_wgs84, 900999), 4269)      The problem is that the transformed values don't seem to match the transformations that ArcGIS are using, so if I have a mix of WGS84 data and custom-transformed into NAD83 data, ArcGIS's WGS84 -> NAD83 transformation is different from PostGIS and the data doesn't line up. The errors are small (1-2 meters), but large enough to see. Can anyone suggest an appropriate PROJ4 string to replicate the results of the ArcMap transformation, or something better for the state of Alaska such as using a NAD grid? Or tell me what I'm doing wrong with my custom projection?


title: Displaying raster attributes on data driven pages

text: I have a feature class of 50 properties that are ranked. I have a photo for each property stored in a raster attribute column (in the feature class) called "Picture". I am trying to create a Dynamic text attribute for the picture column to use with my Data Driven Pages. I have struck out so far on getting the images to appear in the layout. Any ideas? Help would be greatly appreciated. -Adam


title: When is Geographic Weighted Regression used and why?

text: I am currently trying to better understand what geographic weighted regression is, and what it is useful for. I have been reading up on it, but would like to get others people's views on it and see what they have to say. I would like to know when you would want to use, and why is it useful? I graduated with my degree in GIS last December, and have been in the GIS workforce for almost a year now. I started out in utilities, but am now in the 911 field, so I am trying to learn more about hot spots, and using weighted regression.


title: How can I transform polylines into points every n metres in PostGIS?

text: Question in pretty much the same as this, but I need to use PostGIS. I have a number of polylines, and want to trasform them into points every n metres along these lines. **edit:** Thanks a lot to Stefan. My final query was a bit different, probably due to older version of PostGIS (1.5), but in the end I got what was needed. Points were combined into a heatmap of roads density. ![enter image description here](http://i.stack.imgur.com/3V3E4.jpg)


title: Geoserver SLD using PropertyIsEqualTo And PropertyIsLike

text: I have two SLD rules one using PropertyIsEqualTo and another using PropertyIsLike **Using PropertyIsEqualTo**               <Rule>               <Name>Vendor FirstVendor</Name>               <Title>FristVendor</Title>               <ogc:Filter>                 <ogc:PropertyIsEqualTo>                 <ogc:PropertyName>Vendor</ogc:PropertyName>                 <ogc:Literal>firstVendor</ogc:Literal>                 </ogc:PropertyIsEqualTo>               </ogc:Filter>               <PointSymbolizer>                 <Graphic>                   <Mark>                     <WellKnownName>circle</WellKnownName>                     <Fill>                       <CssParameter name="fill">#0033CC</CssParameter>                     </Fill>                   </Mark>                   <Size>8</Size>                 </Graphic>               </PointSymbolizer>                        </Rule>      **Using PropertyIsLike**               <ogc:Filter>                 <ogc:PropertyIsLike  wildCard="*" singleChar="." escape="!">                 <ogc:PropertyName>Vendor</ogc:PropertyName>                 <ogc:Literal>firstVendor*</ogc:Literal>                 </ogc:PropertyIsLike>     </ogc:Filter>      The SLD using Like works fine where as using equal does not work at all.I have seen the string in database well there is not leading or trailing spaces.


title: How to make routable OSM map

text: I have some routable data for testing my pgRouting functions, and everything is fine. Now I have to download some more countries from OSM. I downloaded Germany shape files, but OSM maps are not routable. OK, next step, I found osm2pgrouting Import Tool but it works only on Linux, I need Windows version. Is there any possibility to make OSM map routable on Windows?


title: GIS used to analyse micro-environment

text: Can anybody give me some articles or show me some link where GIS is used to analyse and manage micro-environment, such as a college campus or a small locality. I need to do a project on micro-environment analysis using GIS, so i need some references from where i can get some ideas. Please help.


title: Legend showing info for world map

text: I'm using ArcGIS 10.1 API for JavaScript. I'm showing dynamic map service layer containing 8 layers and I want to show legend for the currently visible layer (only 1 layer is visible at a time). The legend looks ok, but there's one little problem. The legend contains info even for the world map, which is under my visible layer, but I don't want to show this info about world map. Code snippet:               var legendDijit = new esri.dijit.Legend({             map:map,         id: "fluLegend"     }, tsLegDiv);      How can I achieve to ommit info for world map in the legend ? I haven't found any answer in the other topics containing questions about legend. Thanks g The code after Darksanta answer looks like this: The problem is, that the legend loads correctly the first time when layers are added to map, but when I change the layer visibility by changing the 'Year' filter (each layer is visible only when certain year is selected), the legend is as it was and does not reflect the new visible layer. I know, that the Darksanta example (by ArcGIS) is ok, but there must be something fishy in my code, that causes the legend to not work correctly. Maybe the code is completely incorrect. I'm not familiar with js or arcgis, it's my first js + arcgis experience.               <script language="Javascript">         dojo.require("dijit.layout.BorderContainer");         dojo.require("dijit.layout.ContentPane");         dojo.require("dijit.layout.AccordionContainer");         dojo.require("esri.map");         dojo.require("esri.dijit.TimeSlider");         dojo.require("esri.tasks.query");         dojo.require("esri.layers.FeatureLayer");         dojo.require("esri.dijit.Legend");              var map;         var fluLayers, featureLayer;         var timeExtent;         var timeSlider;         var queryTask, query;//, featureSet;         var infoTemplate;              //konstanty         var defaultLayerIndex = 7;         var defaultYear = 2006;         var defaultAgeStart = 6;         var defaultRangeAge = 8;              function init() {             map = new esri.Map("map", {             basemap: "streets",             center: [19.63, 48.60],             zoom: 8             });                  fluLayers = new esri.layers.ArcGISDynamicMapServiceLayer("<map service address not accessible from outside>");             fluLayers.setMinScale(0);             fluLayers.maxRecordCount = 100;             fluLayers.setOpacity(0.5);                  fluLayers.setVisibleLayers([defaultLayerIndex]);             setLayerDefinitions(defaultAgeStart, defaultRangeAge);             map.addLayers([fluLayers]);                  dojo.connect(map, "onLayersAddResult", updateSlider);             dojo.connect(map,'onLayersAddResult',function(results){                         var layerInfo = dojo.map(results, function(layer,index){                                     return {layer:layer.layer,title:layer.layer.name};                         });                         if(layerInfo.length > 0){                                     var legendDijit = new esri.dijit.Legend({                                                 map:map,                                                 layerInfos:layerInfo                                     },"legendDiv");                                     legendDijit.startup();                         }             });         }              //Funkcia, ktora nastavi definiciu aktualne viditelnej vrstvy.         function setLayerDefinitions(ageStart, ageRange){             var layerDefinitions = [];             var layerDefString;             var visibleLayerIndex = getVisibleLayerIndexByYear();                  var startDate = getStartDateByYear();             var endDate = getEndDateByYear();                  if (ageRange === 0)                 layerDefString = "Chripka.STARTAGE = " + ageStart + " AND Chripka.RANGEAGE IS NULL AND Chripka.PERIODSTART >= date \'" + startDate + "\' AND Chripka.PERIODSTART < date \'" + endDate + "\'";             else                 layerDefString = "Chripka.STARTAGE = " + ageStart + " AND Chripka.RANGEAGE = " + ageRange + " AND Chripka.PERIODSTART >= date \'" + startDate + "\' AND Chripka.PERIODSTART < date \'" + endDate + "\'";                  layerDefinitions[visibleLayerIndex] = layerDefString;             fluLayers.setLayerDefinitions(layerDefinitions);         }              //Funkcia, ktora aktualizuje timeSlider v pripade nejakej zmeny.         function updateSlider() {             if (dijit.byId('timeSlider')) {                 dijit.byId('timeSlider').destroy();             }                  var tsDiv = dojo.create("div", null, dojo.byId('timeSliderDiv'));                  timeSlider = new esri.dijit.TimeSlider({                 style: "width: 1000px;",                 id: 'timeSlider'             }, tsDiv);                  map.setTimeSlider(timeSlider);             timeSlider.setThumbCount(2);                  if (!timeExtent)                 setTimeExtentByYear();             timeSlider.createTimeStopsByTimeInterval(timeExtent, 1, 'esriTimeUnitsWeeks');             timeSlider.setThumbIndexes([0, 1]);             timeSlider.setThumbMovingRate(4000);             timeSlider.startup();             dojo.byId("daterange").innerHTML = "<i>" + timeSlider.timeStops[0].toLocaleDateString() + " - " + timeSlider.timeStops[1].toLocaleDateString() + "<\/i>";                  //add labels for every other time stop             var labels = dojo.map(timeSlider.timeStops, function(timeStop, i) {                  if ( i % 2 === 0 ) {                     return i + 1;                 } else {                     return "";                 }             });                  timeSlider.setLabels(labels);                  dojo.connect(timeSlider, "onTimeExtentChange", function(timeExtent) {                 var startValString = timeExtent.startTime.toLocaleDateString();                 var endValString = timeExtent.endTime.toLocaleDateString();                      dojo.byId("daterange").innerHTML = "<i>" + startValString + " - " + endValString + "<\/i>";             });         }              function getStartDateBySlider(){             var startDate;                  if (timeSlider){                 if (timeSlider.getCurrentTimeExtent()){                     var tex = timeSlider.getCurrentTimeExtent().startTime;                     startDate = (tex.getUTCMonth() + 1).toString() + "/" + tex.getUTCDate().toString() + "/" + tex.getUTCFullYear().toString();                 }             }                  if (!startDate)                 startDate = getStartDateByYearExt("2006");                  return startDate;         }              function getEndDateBySlider(){             var endDate;                  if (timeSlider){                 if (timeSlider.getCurrentTimeExtent()){                     var tex = timeSlider.getCurrentTimeExtent().endTime;                     endDate = (tex.getUTCMonth() + 1).toString() + "/" + tex.getUTCDate().toString() + "/" + tex.getUTCFullYear().toString();                 }             }                  if (!endDate)                 endDate = getEndDateByYearExt("2006");                  return endDate;         }              //Funkcia, ktora je vykonana potom, ako dojde ku zmene roku.         function yearChange() {             //if the slider is currently playing stop and reset it             if (timeSlider.playing) {                 timeSlider.pause();             }                  setTimeExtentByYear();                  var ageStart = getSelectedStartAge();             var rangeAge = getSelectedRangeAge();             onUpdateLayerDefinitions(ageStart, rangeAge);         }              //Funkcia, ktora nastavi timeExtent podla zvoleneho roku.         function setTimeExtentByYear(){             var year = dojo.byId('year').value;                  var startDate, endDate;             startDate = getStartDateByYearExt(year);             endDate = getEndDateByYearExt(year);                  timeExtent = new esri.TimeExtent();             timeExtent.startTime = new Date(startDate + " UTC");             timeExtent.endTime = new Date(endDate + " UTC");                  dojo.byId("timeSliderYear").innerHTML = year;             dojo.byId("daterangeYear").innerHTML = "<strong>" + timeExtent.startTime.toLocaleDateString() + " - " + timeExtent.endTime.toLocaleDateString() + "<\/strong>";             dojo.byId("ageRange").innerHTML = getSelectedAgeRangeString();         }              //Funkcia, ktora vrati zaciatocny datum pre zvoleny rok.         function getStartDateByYear(){             var year = dojo.byId('year').value;             return getStartDateByYearExt(year);         }              //Funkcia, ktora vrati koncovy datum pre zvoleny rok.         function getEndDateByYear(){             var year = dojo.byId('year').value;             return getEndDateByYearExt(year);         }              //Funkcia, ktora vrati zaciatocny datum pre zvoleny rok, ktory pride ako stringovy parameter.         function getStartDateByYearExt(year){             return "09/29/2006"; //SIMPLE         }              //Funkcia, ktora vrati koncovy datum pre zvoleny rok, ktory pride ako stringovy parameter.         function getEndDateByYearExt(year){             return "12/29/2006"; //SIMPLE         }              //Funkcia, ktora vrati zvoleny zaciatocny vek.         function getSelectedStartAge(){             return defaultAgeStart; //SIMPLE         }              //Funkcia, ktora vrati ku zvolenemu zaciatocnemu veku jeho rozsah.         function getSelectedRangeAge(){             return defaultRangeAge; //SIMPLE         }              //Funkcia, ktora vrati zvoleny vekovy rozsah ako retazec.         function getSelectedAgeRangeString(){             var radios = dojo.query(".list_item");                  for (var i=0, il = radios.length; i < il; i++) {                 if (radios[i].checked) {                     return radios[i].value;                 }             }                  return "";         }              //Funkcia, ktora aktualizuje definiciu aktualne viditelnej vrstvy a spravi vsetko ostatne spojene s touto aktualizaciou vrstvy.         function onUpdateLayerDefinitions(ageStart, ageRange){             //default je vekove rozpatie 6-8             if (typeof ageStart === 'undefined')                 ageStart = defaultAgeStart;                  if (typeof ageRange === 'undefined')                 ageRange = defaultRangeAge;                  var visibleLayerIndex = getVisibleLayerIndexByYear();             fluLayers.setVisibleLayers([visibleLayerIndex]);                  setLayerDefinitions(ageStart, ageRange);                  updateSlider();         }              //Funkcia, ktora vrati index viditelnej vrstvy, tzn. zavisi to od zvoleneho roku.         function getVisibleLayerIndexByYear(){                 return 6; //SIMPLE             }         }              //Funkcia, ktora je vykonana potom, ako dojde ku zmene filtracneho kriteria.         function filterChange() {             //if the slider is currently playing stop and reset it             if (timeSlider.playing) {                 timeSlider.pause();             }                  setTimeExtentByYear();                  var ageStart = getSelectedStartAge();             var rangeAge = getSelectedRangeAge();             onUpdateLayerDefinitions(ageStart, rangeAge);         }              dojo.ready(init);          </script>


title: Aiming to retrieve source and destination from Android, use pgrouting dijkstra to find shortest path and display it back on Android

text: I have created an Android application that will accept source and destination from the user. It then calls a php page which takes in the lattitude and longitude of the destination and the source. I am using postgis in the backend, and have created the topology required for pgrouting. Following is my php code: First I find the all the source integers thats near to my current location using the following function:-                   function         getSourceId($currentPositionLongitude,$currentPositionLattitude,$approximation,$dbConnectionString){         // the query will return the source id of the user.        $db_connection = pg_connect($dbConnectionString);      $queryForSource = "select source  from mumbai_roads     where ST_Dwithin(geom,ST_GeomFromText('POINT( $currentPositionLongitude $currentPositionLattitude)'),$approximation)";           $source=array();             $myresult = pg_query($db_connection,$queryForSource);             while($myrow = pg_fetch_assoc($myresult)) {                   $resultSource = $myrow['source'];                array_push($source, $resultSource);                  }         //print_r($source);         return $source;     }      Now i find all the target integers that are near my destination location:-                   function getTarget($destinationLongitude,$destinationLattitude,$approximation,$dbConnectionString){        $db_connection = pg_connect($dbConnectionString);       $queryForDestination = "select target  from mumbai_roads     where ST_Dwithin(geom,ST_GeomFromText('POINT( $destinationLongitude $destinationLattitude)'),$approximation)";       $target = array();         $myresult = pg_query($db_connection,$queryForDestination);         while($myrow = pg_fetch_assoc($myresult)) {             $resultTarget = $myrow['target'];                array_push($target, $resultTarget);         }         //print_r($target);         return $target;     }      Finally as i have both the source and target integers, I try to get the shortest path for all the sources and the targets:-               if(!(empty($source) || empty($target))){           /*looping through all the source and the targets */           foreach ($source as $sourceId) {             foreach ($target as $targetId) {               $json = getShortestPath($sourceId,$targetId,$dbConnectionString);               if($json!=""){                 $result=$result.$json;                   echo "$result";               }             }                }         }          function getShortestPath($source,$target,$dbConnectionString){        $db_connection = pg_connect($dbConnectionString);       $queryForShortestDistance = "select seq, id1 AS node, id2 AS edge, di.cost AS cost1, ST_AsGeoJSON(geom) AS geom       FROM pgr_dijkstra(         'SELECT id, source, target, cost FROM mumbai_roads',         $source, $target, false, false       ) as di       JOIN mumbai_roads pt       ON di.id2 = pt.id";       $geom="";          $myresult = pg_query($db_connection,$queryForShortestDistance);         while($myrow = pg_fetch_assoc($myresult)) {             $geom = $myrow['geom'];              }         return $geom;     }      I am getting the path as a multiline which i am mapping on android, but using the above code, i am not getting the whole path, but just a small fraction of the whole path. Please help!


title: Displaying a dynamic KML file on Google Maps as it is on Google Earth

text: how can i display my kml files through network link on Google Maps instead of showing it on Google Earth. I have tried the step mentioned on this link https://developers.google.com/kml/articles/pagesforkml but stucked at step 3 where it is saying for uploading kml file as no option is available there. Actually ,I'm uploading KML maps because , my KML file size is more than 2 MB and as i refresh it on google earth , it got hanged.. So what i need is a network link on Google maps which will show my my.kml file placemarks and this network link is having a refresh interval time ..my.kml file is generating at every 10 seconds by a java program... my network link code is : // pasted only the network link part               <NetworkLink>     <Link>     <href>/home/placeMarkers.kml</href>     <refreshmode>onChange</refreshmode>     <refreshInterval>60</refreshInterval>     <viewRefreshMode>onStop</viewRefreshMode>     <viewRefreshTime>20</viewRefreshTime>      </Link>     </NetworkLink>     </kml>


title: Leaflet + GeoServer + British National Grid

text: I have followed along some examples and was able to display BNG 27700 data nicely on the map, however now I want to perform getFeatureInfo function on click. Again the code I used is as below: `var BBOX=map.getBounds().toBBoxString(); var WIDTH=map.getSize().x; var HEIGHT=map.getSize().y; var X=map.layerPointToContainerPoint(e.layerPoint).x; var Y=map.layerPointToContainerPoint(e.layerPoint).y; var SRS=map.options.crs.code; var URL="http://localhost:8080/geoserver/OSM/wms?INFO_FORMAT=application/json&REQUEST=GetFeatureInfo&EXCEPTIONS=application/vnd.ogc.se_xml&SERVICE=WMS&VERSION=1.1.1&WIDTH="+ WIDTH+"&HEIGHT="+HEIGHT+"&X="+X+"&Y="+Y+"&BBOX="+BBOX+"&SRS"+SRS+"&LAYERS=OSM:BuildingBoundaries&QUERY_LAYERS=OSM:BuildingBoundaries";` However no feature is returned. If I make use of sample data from Geoserver in leaflet, the click event works perfectly. Am I missing something?


title: How to create feature type from xml attribute value?

text: Here is a sample from an xml child element:               <element>     <aspects>                              </aspects>     <aspects>                              </aspects>     </element>      I want to extract the value "Grundlose Kuhle" from the aspects.a.id="1".v element. How can I get FME to put all "v"-values from ID="1" elements into a list so that I can write them to a featuretype?


title: Focal statistics to calculate saddle points from DEM

text: I have a conundrum with Focal Statistics and kernel files. I am been trying to identify saddle points in DEMs. I am attempting to rewrite this (thanks whuber) in a python script for 10.1. I have used focal statistics and kernel files for irregular neighborhoods in the past with no issues. Now, for some reason, I continuously get the following error: ERROR 010308: Invalid Irregular neighborhood mask My script is below:               # Import system modules     import arcpy, os, codecs     from arcpy import env     from arcpy.sa import *          arcpy.CheckOutExtension("spatial")          # Obtain parameters     dem = arcpy.GetParameterAsText(0)     outFolder = arcpy.GetParameterAsText(1)          arcpy.env.workspace = outFolder     DEM = Raster(dem)     kernel = os.path.join(outFolder, "kernel.txt")          # Compute the neighborhood sums involved in estimating the second partial deriviatives.          arcpy.AddMessage('Attempting FocalStat 1...')     n1 = NbrRectangle(3, 3, "CELL")     g1 = FocalStatistics(DEM, n1, "SUM", "NODATA")          arcpy.AddMessage('Attempting FocalStat 2...')     n2 = NbrRectangle(1 , 3, "CELL")     gX = FocalStatistics(DEM, n2, "SUM", "NODATA")          arcpy.AddMessage('Attempting FocalStat 3...')     n3 = NbrRectangle(3, 1, "CELL")     gY = FocalStatistics(DEM, n3, "SUM", "NODATA")          arcpy.AddMessage('Attempting FocalStat 4...')     text = codecs.open(kernel, "w", encoding='ascii')     text.write("3 3\n1 0 0\n0 0 0\n0 0 1")     text.close()     n4 = NbrIrregular(kernel)     gM = FocalStatistics(DEM, n4, "SUM", "NODATA")          arcpy.AddMessage('Attempting FocalStat 5...')     text = codecs.open(kernel, "w", encoding='ascii')     text.write("3 3\n0 0 1\n0 0 0\n1 0 0")     text.close()     n5 = NbrIrregular(kernel)     gP = FocalStatistics(DEM, n5, "SUM", "NODATA")          # Delete kernel txt     os.remove(kernel)          # Compute the determinant of the Jacobian.     # It is multiplied by 10^6 to bring it into a range that SA can display.     arcpy.AddMessage('Calculating determinate of the Jacobian...')     xCell = arcpy.GetRasterProperties_management(DEM, "CELLSIZEX").getOutput(0)     g = (g1 - (gX * 3)) * (g1 - (gY * 3)) / 9.0 - (((gP - gM) / 4.0)**2) / (int(xCell)**4 / 1000000)     g.save("JacobDet.img")      I have tried using the focal stat tool in arc with a stand alone txt file with the same results. Has anyone encountered this problem before? Is there something wrong with my code/kernel file? Thanks in advance for any comments.


title: Extract tiff from WPS execute response

text: I am very new with OWS, specially WPS, and I am testing 52North WPS for a project. I have run the Sextante visibility process, which result format is a tiff image. The execute response starts as follows:               <?xml version="1.0" encoding="UTF-8"?>     <ns:ExecuteResponse xmlns:ns="http://www.opengis.net/wps/1.0.0"         xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"         xsi:schemaLocation="http://www.opengis.net/wps/1.0.0 http://schemas.opengis.net/wps/1.0.0/wpsExecute_response.xsd"         serviceInstance="http://localhost:8081/wps/WebProcessingService?REQUEST=GetCapabilities&SERVICE=WPS"         xml:lang="en-US" service="WPS" version="1.0.0">         <ns:Process ns:processVersion="1.0.0">             <ns1:Identifier xmlns:ns1="http://www.opengis.net/ows/1.1">visibility</ns1:Identifier>             <ns1:Title xmlns:ns1="http://www.opengis.net/ows/1.1">Visibilidad</ns1:Title>         </ns:Process>         <ns:Status creationTime="2013-04-17T12:13:16.298+02:00">             <ns:ProcessSucceeded>Process successful</ns:ProcessSucceeded>         </ns:Status>         <ns:ProcessOutputs>             <ns:Output>                 <ns1:Identifier xmlns:ns1="http://www.opengis.net/ows/1.1">RESULT</ns1:Identifier>                 <ns1:Title xmlns:ns1="http://www.opengis.net/ows/1.1">Visibilidad</ns1:Title>                 <ns:Data>                     <ns:ComplexData mimeType="image/tiff"><![CDATA[MM?*?????????????????????      To extract the image to a file (to store it in Geoserver) I get the CDATA, decode with Base64Decoder to a byte array and the write it to a file, like visibility.tiff. But this image seems to be "corrupted", I can't preview it and Geoserver fails when I try to create the datastore with the error: "Could not list layers for this store, an error occurred retrieving them: El argumento "value" no debe ser nulo." The last sentence means the argument "value" must not be null. Am I doing something wrong? Thanks, Santi


title: Problem loading PostGIS (Postgres) View layer in QGIS

text: Why is Postgres user XX via command „Add PostGIS Layers“ in QGIS enabled upload only View tables which owns? Or how to achieve it to the user XX via PostGIS can load View table which owns the user YY? Main table from the database with different users are assigned by rolls, however the View tables we could not find that there is a possibility of assigning rolls. Any idea how solve this problem?


title: Convert cursor output to list of geometries

text: I want to convert all the geometries returned by search cursor to a list.               import math     import arcpy     from arcpy import env          env.workspace = r"C:\GIS Data\GIS data_for Maeenul vai"          desc = arcpy.Describe("River.shp")     shapefieldname = desc.ShapeFieldName          rows = arcpy.SearchCursor("River.shp")          featureList = []          for row in rows:              feat = row.getValue(shapefieldname)         featureList.append(feat)              print "%i %i" % (feat.firstPoint.X, feat.firstPoint.Y)         print         print "%i %i" % (feat.lastPoint.X, feat.lastPoint.Y)         print         print          del row, rows     print "---------------------------------------------------------------"          for feat in featureList:         print "%i %i" % (feat.firstPoint.X, feat.firstPoint.Y)         print         print "%i %i" % (feat.lastPoint.X, feat.lastPoint.Y)         print         print         print      Output is > 3610930 2135882 3611593 2134453 > > 3611806 2134981 3611593 2134453 > > 3614160 2136164 3617432 2131734 > > 3611593 2134453 3617432 2131734 > > 3617432 2131734 3620568 2127591 > > 3620568 2127591 3620785 2127423 > > 3617980 2126657 3620568 2127591 > > 3616768 2129454 3617948 2126649 > > 3617948 2126649 3617980 2126657 > > 3615102 2128889 3617587 2126510 > > 3617587 2126510 3617948 2126649 > > 3617624 2126416 3617980 2126657 > > 3613129 2128176 3615155 2125617 > > 3615155 2125617 3617587 2126510 > > 3615086 2125515 3615155 2125617 > > * * * > > **List Data** > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 > > 3615086 2125515 3615155 2125617 So basically the append methods of list only copies the last element returned by cursor to the whole list. How can I resolve this?


title: Merge multiple geometry type TAB files in QGIS/OGR

text: I have a number of vector layers with all types of geometry (point, line/multiline, poly) and I want to combine them into a single layer in QGIS. I know MapInfo can do it but with the other 99% of my workflow in QGIS it's the last step to having a complete pipeline. It seems there are a number of tools to merge shapefiles together (assuming equivalent geometry types) in QGIS but this doesn't help me. There's ogr2ogr but running the following statement gives me a failure:               ogr2ogr -f "MapInfo File" merge.tab points.tab     ogr2ogr -update -append merge.tab lines.tab -f "MapInfo File"     FAILURE:     Unable to open existing output datasource 'merge.tab'      It'd be fine if I wanted to merge a bunch of layers with the same geometry (use SHP, merge, save result as TAB), but multiple geometry types are eluding me. I could write a little program using Tatuk Developer Kernel but I'm not that happy with its handling of projections (using EPSG:28355 it seems to create a projection that is 'custom' when viewed in QGIS or MapInfo).


title: postgis problem with shortest distance calculation

text: while working with POSTGIS pgrouting, for calculateing the distance between two roads(lines) i got the shortest_path function. But the logic is based on Start_point(Start_id) and end_point(end_id) but in my data the linestring contains so many internal points like ('linestring(1 1,2 2,3 3,4 4,5 5)' just for example..) it is taking start point (1 1) endpoint(5 5) if other line starting with (5 5) it is showing as route...like ('linestring(5 5,6 6)') But line which crossing the point inside the linestring like(2 2,3 3,4 4) which is not telling as connected.. example table roads: id name way 1 A linestring(1 1,2 2,3 3,4 4,5 5) 2 B linestring(5 5,6 6) 3 c linestring(2 1,2 2,2 3) if i am applying shortest_path function from point(1 1) to (6 6) its showing the way but for (1 1) to (2 3) it is not showing anything...but there is a route for this (1 1,2 2,2 3) can anyone please help me out for finding the solution.. Regards Deepak M


title: How to meassure the length of common edges between polygons?

text: I have a shapefile that contains several polygons which in fact represent the limits of several counties. My task is to find each length of each shared border. So in other words I want to know how many km of border share the A county with the B county and so on. Of course I am looking for a tool to do this automatically. Measuring them manually would take ages. I envision some sort of matrix, but not sure how to do that in ArcGIS. I use ArcGIS Desktop 10.0 (ArcEditor) . 10x


title: QGIS how to compute a trendline of two point lines?

text: Does anyone know how to compute a trendline for two lines of points? I have two lines of points of a 10 km long urban road. 1st line is one driving direction and the 2nd line of points ist the other driving direction. Is there a way to get a trendline of these both point lines also as a point line? Thanks!!


title: Create a style file using ArcEngine?

text: Is it possible to create a style file to complement a shapefile using ArcEngine? I'm using C# VS2010 ArcEngine 10. Thanks.


title: Data interpreted incorrectly when opening NetCDF file in ArcMap

text: I'm trying to open a NetCDF file with some current data on a polar stereographic projection in ArcMap, but the data are not displayed correctly. The only obvious problem I can see is that ArcMap seems to be interepreting some of the global attributes incorrectly. For example, in ncdump it says               :false_easting = 2384800.f ;      Whereas in ArcMap it says               False_Easting 1.890518526342814e-273      Is there anything I can do to ensure the attributes are read correctly? I thought the point of a format like NetCDF was that it would specify both the data and the data type in a standard way.


title: How to access user credentials in database connection?

text: I am trying to automate a couple of things when a QGIS plugin is activated:   1. If Layer A is not in the TOC, add it. This requests username and password to access a PostGIS table using the inbuilt credentials dialog.   2. If Layer A now exists in the TOC, access the same PostGIS table to run a query using the same user credentials. Apparently as I'm not storing the password, I can't access it to use in connecting to the database (using psycopg). Is there any way around this? I want to avoid having the user type credentials twice, and also I want to avoid storing the credentials permanently. I have tried to reimplement the QgsCredentialDialog class but that adds so much complexity (e.g. when a connection error is thrown). What other ways can I add a layer using a non-stored username and password, and access the credentials later to query the same database? The easiest option I can think of is to access the inbuilt credentials dialog request, but I'm not sure how to do that. **EDIT** My implementation of this was to create a function which I call at the point my layer has been added. It then sets my plugin class attributes `username` and `password` correctly which I can use in my db query.               def get_credentials(self, lyr):         connInfo = QgsDataSourceURI(lyr.dataProvider().dataSourceUri()).connectionInfo()         success, username, password = QgsCredentials.instance().get(connInfo, None, None)         if success:             # QgsCredentials.instance().put(connInfo, username, password)             self.username = username             self.password = password


title: Convert Surface to MultiSurface in SQL Server

text: I need to create view where geometry column is MultiSurface datatype. But, column in my table is Surface datatype. Is there any chance how to convert single Surface column to MultiSurface column (MSurface is collection, so I guess it can be collection of single feature). Thanks in advance!


title: Why does gdalbuildvrt not allow these things?

text: I don't know why does gdalbuildvrt not allow:   * Unreferenced images   * Rotated images   * South up images Could someone throw some light? Many thanks in advance


title: How can I associate one point layer to one line layer?

text: **How can I associate in QGIS one point layer to one line layer** A sample scenario could be this one:   * A point layer "Valves" and line layer "Pipelines"   * Each layer should have its own set of attributes   * A "pipeline" has to begin and end in one "valve"   * A change in the "Valves" geometry implies a change in the "Pipelines" geometry, i.e., if I move a "valve", the "pipelines" connected to that node should follow the node wich are connected to


title: Polar grids in ArcGIS

text: I am looking automate the creation of a polar grid shapefile (polygon). I am currently creating this in a many-step editing process, starting from a concentric buffer. Here is how I think I should be approaching the issue. Given a polar grid of this shape: ![enter image description here](http://i.stack.imgur.com/5Mjww.png) Where V1 represents the origin, and the length of the line segment (v1,v4) represents the buffer distance. V2 and V3 represent points along the concentric ring that intersects the polar lines. The amount of polar lines can be changed. Any thoughts on how to do this?


title: How to get US Highway exit geocodes using Google maps API or through some other means?

text: For our application, we need to get all the US interstate highway exits along with their Geocodes. Can you please explain how to get these using Google Maps API or Google Local search or otherwise? So far I have no clue as to how to proceed. I greatly appreciate if you can point me in right direction. I am looking for a free or a paid solution that is around $1000. I know there are people who provide geocodes. But they are out of our budget. Several people are suggesting TIGER data with ArcGIS. I know nothing of them. I appreciate if you can help in this regard.


title: installing module in embeded python in qgis

text: I'm newbie to Quantum GIS and Python. I'm using QGIS in windows platform. How I can install module for python embedded in QGIS? Also, how can I run embedded python from windows directly? Thanks


title: Moving raster to different mapset GRASS GIS

text: I have a raster DEM that I created in one mapset but I would like to copy it to another mapset for another purpose. ![alt text](http://i.stack.imgur.com/thAb0.png) Is this possible within GRASS or can I just copy some files over in the filesystem to the FloodModels mapset folder on my hard drive?


title: Load One Layer from ArcGIS Server MapServer

text: I am pretty new to ArcGIS Server and am building a map applications using the Javascript API. I can easily load the map layer below by adding it as a ArcGISDynamicMapServiceLayer http://www.parkcounty.org/ArcGIS/rest/services/Roads_Map/MapServer/ As you can see there are many layers as part of the above url. What if I just want to load the roads layer as a ArcGISDynamicMapServiceLayer? I've tried replacing the above url with the road specific URL below and that does not seem to work http://www.parkcounty.org/ArcGIS/rest/services/Roads_Map/MapServer//2 Is there something else I can do as the front end developer or is it up to those that are hosting and managing the server to create a unique url for the roads?


title: Is it common to save geographic coordinates in minutes instead of degrees?

text: Are there any communities or fields where it is common to save geographic coordinates in minutes instead of degrees? (I just received a batch of addresses referenced that way and it made me wonder.) Are there any advantages in saving lat/lon in _'decimal minutes'_? For example:               2894.511,982.8775      instead of               48.2418,16.3812


title: Draw great circle rays from a single point along specific directions

text: How can I draw great circle rays starting from a single given point P (for example, from coordinate [34 N, 116 W] in this map) and along directions D1, D2, ... (for example, 30, 45 and 135 degrees) in QGIS? I do not have any shape files or other vectors for the lines. All I have is a point which can either be imported or selected interactively. Note that I'd like to plot _rays_ , not paths, since I don't have a destination point (only the direction). However, the map will finally be clipped to its extents, so if it makes things simpler, I can also generate points outside the map with the desired direction and import them.


title: GeoDjango: order PointField by distance?

text: I've got a GeoDjango model as follows:               class City(models.Model):         id = models.IntegerField(primary_key=True)          county = models.ForeignKey(County)         location = models.PointField(null=True, blank=True)      How do I get all the cities within a particular bounding box, ordered by distance? I've got the relevant latlngs as strings, but I'm not really sure how to construct the bounding box, or the central point, or the query :) EDITED: have now worked out how to construct the bounding box and get places - but can anyone help me sort by distance?               swLat = request.GET.get('swLat')     swLng = request.GET.get('swLng')     neLat = request.GET.get('neLat')     neLng = request.GET.get('neLng')     centreLat = request.GET.get('centreLat')     centreLng = request.GET.get('centreLng')     centre = fromstr('POINT(%s %s)' % (centreLat, centreLng))     boundaries = fromstr('POLYGON((%s %s,%s %s,%s %s,%s %s,%s %s))' % \            (swLng, swLat, swLng, neLat, neLng, neLat, neLng, swLat, swLng, swLat))     cities = City.objects.filter(location__within=boundaries)     # how to order cities by distance from centre?      Thanks for helping out a beginner!


title: Can I use a join in Mapviewer

text: I've got two tables, parcels and permits. I'd like to display permits on a map without having to copy the geometry data from parcels. Is there a way to create a join in Mapviewer? Do I have to create a view?


title: How to reproject a raster in QGIS?

text: I am a history teacher and I like to make my own maps (I can never find exactly what I want). I am using QGIS 1.8 and Natural Earth data. I want to reproject their cross-blended hypso rasters into projections that are more familiar looking to my students, like Lambert Conic Conformal. I have tried changing the crs for the whole project, and the warp and assign functions from the raster menu, and changing the crs for the raster layer itself. None of these seem to work properly - the resulting maps are unrecognizable and extremely large (file size). I feel like I am making some kind of basic mistake but I don't know where to begin with this. Any help would be appreciated.


title: Recommended Server Specs for GeoServer

text: we serve our Geographic data via our flagship application Dorset Explorer. This uses OpenLayers and GeoServer to serve our Geographic data, including OSM and OS Basemaps and a mass of other data such as Aerial Photography, Schools etc. My question is, what sort of Server Specs should we be using for GeoServer. We will only have available to us one server for delivery of tiles via WMS (and perhaps in the future we may be using some WFS on it as well). We get between 100 and 200 users per day, but this will hopefully start rising as soon as we are confident that the servers will hold up! Is it just a matter of stick in the highest specs we can afford or are there some limits or theoretical maximums. Any pointers would be much appreciated!


title: Select polygons contained inside a polygon and assign ID

text: I have a polygon feature class `Base_FC`, where each poly has a unique alphanumeric ID. My second polygon feature class `Type_FC` has the details for the various polygons which would fall inside `Base_FC` when intersected. These polys have unique alphanumeric IDs as well. For display purposes, I would like to assign each `Type_FC` poly that falls within a `Base_FC` poly a number, along with its neighbouring polys. For example:   * `Base_FC` polygon B1 _contains_ three `Type_FC` polys T5, T8 and T11.   * `Base_FC` polygon B2 _contains_ four `Type_FC` polys T20, T21, T25 and T30. The numbers should ideally be assigned spatially from left to right, so if B1 has the polys going from left to right as T8, T5 and T11, they would be labelled 1, 2 and 3 respectively. I realise that there are a number of ways to do this - the method I was thinking of was to create a model which iterates through each `Base_FC` feature, selects by location all the `Type_FC` features it contains, then use the Python field calculator to number them - but I'd like to know of other methods, or if there is a tool out there which can do this. I'm also not sure how to handle the issue of numbering them so that visually it looks good. Maybe using the coordinates of the centroids to ensure numbering goes from left to right?


title: Resolving ArcSDE for Oracle 11g Layer not found (-20) error?

text: I am using ArcSDE for Oracle 11g, with Oracle spatial database. I use the sdelayer command ot register my feature class. Surprisingly I get this error of **Layer not found(-20)** error. The table does exist in Oracle and I had configured feature classes using the same data source. An **Esri Technical Article** suggests to increase tablespace, but I do have free spaces in my tablespace. Please advise me on this?


I am exporting OSM data (in .osm) files to use in qGIS. I am interested in having the street names displayed. I know they are in the OSM data set because they are displayed when I browse the data online, but when I view my downloaded files and open the attribute tables for streets, most of them simply have "null" or "label" as values. What am I doing wrong? How can I extract data from OSM which includes the street names?


title: finding the long axis of an irregular polygon

text: Hi need some help in how to find the long axis of an irregular shaped ellipse using ArcGIS. I have found the centroid but don't know how to determine the longest distance.


title: Why does QGIS add decimals when a csv is saved as a shapefile?

text: When I save a csv as shapefile the fields with decimals are saved adding decimals, then the new table in some case has new values. Old Table (CSV)               ALTURA ELI  ALTURA ORT     1.94        39.6     -11.19      28.22     -36.88      2.59     -33.64      6.08      New Table (DBF/SHP)               ALTURA ELI  ALTURA ORT     1.940000000000000,39.600000000000001,     -11.190000000000000,28.219999999999999,     -36.880000000000003,2.590000000000000,     -33.640000000000001,6.080000000000000,


title: QGIS html annotation relative paths

text: I find the QGIS html annotation tool is very useful for quickly accessing a combination of key data attributes and photos. I am having trouble using relative paths for images. This doesn't work:                      a.k.a.                     I'm wondering if it's possible to employ a QgsProject.instance() in the QGIS Macro option under project properties. I don't know html very well, so I'm wondering how I could link the relative pathname to the project path using something like this:               proj = QgsProject.instance()     UriFile = str(proj.fileName())     pathname = str(os.path.dirname(UriFile))      How can I concatenate this in the html display option?


title: Imagine files between ERDAS and ENVI

text: I have ENVI 5.0 , ERDAS 2013 SP2, and ArcGIS 10.1 SP1. Every time I save an Imagine file from ENVI 5.0 and open in ERDAS 2013 and it seems like ERDAS does not like it and it would not open and would close the program. Does anyone out here have the same experience ? or am I the odd man out here ?


title: Cropping or removing borders (not NoData) from multiple scanned and georeferenced maps (ArcMap or QGIS)

text: I have several hundred scanned maps in TIFF format, each correctly georeferenced and forming a large grid of images. The problem is that the border for each map is present - a white band, with a black edge (not NoData), containing references etc, and overlaps the adjacent tiles. I want to remove these borders from all maps (see below - the white cross through the middle and the edges around the outside). ![4 scanned maps showing white, overlapping borders through the middle and around the edge](http://i.stack.imgur.com/DPcCB.png) Is there a way to remove X number of pixels from the top, left, bottom, and right of rasters? And is there a way to do it for multiple rasters in the same process? I.e., rather than using the extent (as with the Clipper tool). It would be very easy in an image processing program like IrfranView, but I want to retain the geo-data (World files etc) which would be thrown off as the extent and position of the corners would change. I can use QGIS or ArcMap.


title: Looking for free tools to automatically generate spatial data catalogue

text: I've recently inherited a whole bunch of spatial data (mostly shapefiles but some .gdbs) which relate to a bunch of different cities. I'm looking for a way to take an inventory of what data I've got for each city, preferably automated! Does anyone know of any (free) tools that would allow me for example, to create a catalogue of my data (and metadata) from the files contained within my directory structure? I'm using Linux (Ubuntu 14.04) and tend to use QGIS/PostGIS.


title: How do I install the OpenGeo Suite on a web server?

text: I've been doing some research over the past few days to learn just how I can bring my parcel, right of way, easement etc... GIS data out to the field with me on my Android mobile devices with real-time GPS capabilities. I'm not in a position to spend hundreds of dollars a month for an ArcGIS subscription, so I'm trying to learn OpenGeo Suite instead. I've very limited experience with web servers, linux, etc... so I know I've got a lot of reading ahead of me, but I'm hoping someone here can put me on the right path or let me know if I'm just too far out of my league to do this without a lot of prior experience. I do have a website and host running cpanel with a bunch of things already installed including Apache 2.2.26. If I'm understanding the OpenGeo instructions correctly, I need to install Apache Tomcat (or something similar) onto my web server in order to install OpenGeo's production WAR bundle? I see they also have installers for Ubuntu and Red Hat, however I have no idea if this is for a web server or desktop computer or either. Am I on the right track with this line of thinking or way off?


title: Append ArcGIS feature layer data using GDAL/OGR API

text: ArcGIS servers are configured with a maximum records count that they're allowed to return following a REST query. This maximum is usually 1000 records, which means that GDAL/OGR clients should implement an iterator to go through all the data in a layer. My attempt to solve this problem was to use returnIdsOnly=true option to ArcGIS, which results in the list of object ids in the current layer. Next, I read these layers using a query on object id range, and this works fine as far as ArcGIS is concerned. However, GDAL/OGR does not seem to have an easy way of concatenating results from different ArcGIS queries belonging to the same layer. Every time a data source is opened with OGRSFDriverRegistrar::Open function, a separate OGRDataSource is created with its own set of information. I wonder if anyone could guide me on how to merge results of OGRSFDriverRegistrar::Open calls into the same OGRLayer object? I've tried using layer union operation, but it seems to be failing. Any hints?


title: Help me with Trace tool in ARCGIS10

text: Iam digitizing this map when I use trace is there any way to trace black lines? (not digitized) bmp file it was Easy way to make polygons backed is bmp map image


title: How to plot approximate geospatial data on a historical map?

text: Our library is undertaking a project to map historical data regarding 18th and 19th century British and American printers and publishers. One of the problems that we have encountered is the different granularity of location data. For example, one publisher may list a street number, street name, and city whereas another publisher may only mention a street name or even only the city. In such cases, we are wondering the best way to georeference the location and present this approximate data. Does the spatial humanities community have rules or conventions for dealing with these cases? Also, what are the most effective ways for visualizing “approximate vs precise” geospatial data on the same map? Any suggestions or resources that you can offer would be greatly appreciated!


title: For multiple rows of values on singular route feature, how does ArcMap symbolize by quantities?

text: So here is my situation: I have a road segment that is one feature. Along this road I have data that is captured via GPS pings (which generally occur every 5 seconds), and are delineated by 'START_DATE' and 'END_DATE'; so there are multiple pings that make up the information per 'capture' or row. ![Road Segment \(one feature\)](http://i.stack.imgur.com/tkMRh.jpg) ![START//END DATE](http://i.stack.imgur.com/4FyJU.jpg) The question I have is with the fact that I have multiple rows of values that encapsulate the entire segment and therefore occupy the exact spatial area of the line. One of the columns in the table is 'QUANTITY_SOLID', which contains the total amount of solid materials that were disseminated during that 'START_DATE / 'END_DATE' instance. ![Quantity_Solid](http://i.stack.imgur.com/4dcCu.jpg) I then attempt to symbolize by 'Quantities > Graduated Colors' ![Graduated Colors](http://i.stack.imgur.com/UOFqw.jpg) I understand that because it is only one feature class, and all of the information contained in the rows along that line segment all cover the same singular feature class, that there will not be any differentiation between the colors when symbolized. However, this leads to my actual question: **_How does ArcMap choose what color to make the line if there are multiple values over the same spatial area? Does it take the sum of all the 'QUANTITY_SOLID' and average them? Or something else?_**


title: How to find and replace a singe word in a string?

text: How do i find and replace a single word of a string? e.g. Calle Iparraguirre to Iparraguirre I've tried:               replace( 'Calle\\s', 'Calle\\s', '')               replace( "name",'calle','')      and all combinations of code, but nothing is working :/ Thanks in advance! *i've read other q&a's but no result


title: Copy ArcGIS symbol files to QGIS

text: Is there a way to copy the ESRI symbol files for use in QGIS. if not id there a resource to get various svg files to use in QGIS? DIfferent linetypes as well?


title: How can I assess the accuracy of a reprojection?

text: I'm reprojecting a set of measured GPS co-ordinates from WGS84 to OSGB36, using two methods: the pyproj.transform method, which uses PROJ.4, and a "homemade" method from here. The latter method is accurate enough, but I'd like to know how best to measure its accuracy compared to the PROJ.4-based transform (both with and without the OSTN02 datum shift).


title: How join attributes for point feature from nearest points in road network

text: I have two layers: 'linear' - road network and 'point' - points in each linear features intersections. Point layer have in attribute table name for each point features: A,B,C,D,E,F, etc. How add in point layer names of nearest points of road network? For example:               for point A in attribute table should be attributes B,E;     for point B in attribute table should be attributes A,C;     for point C in attribute table should be attributes B,D;     for point D in attribute table should be attributes C,H,N,I;      In QGIS using Distance Matrix possible similar task, but there is "nearest distance", not road network. I think it's done in PostgreSQL/PostGIS, but don't find solution. Thanks! ![enter image description here](http://i.stack.imgur.com/p0v4X.jpg)


title: Popups on Layer (OpenLayers)

text: I have the following situation: i display points and point informations in popups. the points are added to a layer. The popups are added to the map, with map.addPopup function. It works great. The problem is that i want to control the visibility of the popups. For example, if i should be able to add the popups to a layer, then i should be able to hide/show that layer with LayerSwitcher (OpenLayers.Control.LayerSwitcher). But i haven't found any example when somebody adds a popup to a layer. Every example uses the map.addPopup function. Any idea?


title: Create new regions in layer by filling in between polylines on existing layer

text: I'm not sure if this is possible but would save me hours of work. I have a layer which consists of polylines. Most of those polylines bound areas; houses, roads, etc. Currently I am required to colour each of those areas to define what they represent in the real world. My current technique is drawing a region in a new layer; as you can guess this takes a lot of time, especially for complicated areas. Is there a feature in QGIS whereby I can create a new region on a layer by selecting a point within the bounded area and it creates the region shape? EDIT I have tried converting the whole layer of polylines to regions but the results were not as desired, there are some lines that are not bound to an area (partial fences, walls etc) so the results look odd. It would be something similar to 'hatching' in an CAD program, to fill the area bound by lines. Only I want to do this to create the region.


title: How to generate random locations nearby my location?

text: I am trying to create random locations nearby my location. What i want is to create random latitude/longitude pairs inside a 200 meters circle surrounding my location. This is the formula i came up with (with the help of people at StackOverFlow): (Random number between -1 and 1)*radius + (old longitude) = new longitude within radius of old longitude (Random number between -1 and 1)*radius + (old latitude) = new latitude within radius of old latitude The thing is that something weird is happening with my implementation because all the random locations are too near of my location center, it seems that the formula does not cover the whole radius. Any idea of what could be wrong with my formula? Edited to show the current java implementation:               public static Location getLocation(Location location, int radius) {         Random random = new Random();              // Convert radius from meters to degrees         double radiusInDegrees = radius / METERS_IN_DEGREES;              double x0 = location.getLongitude() * 1E6;         double y0 = location.getLatitude() * 1E6;         double u = random.nextInt(1001) / 1000;         double v = random.nextInt(1001) / 1000;         double w = radiusInDegrees * Math.sqrt(u);         double t = 2 * Math.PI * v;         double x = w * Math.cos(t);         double y = w * Math.sin(t);              // Adjust the x-coordinate for the shrinking of the east-west distances         double new_x = x / Math.cos(y0);              // Set the adjusted location         Location newLocation = new Location("Loc in radius");         newLocation.setLongitude(new_x + x0);         newLocation.setLatitude(y + y0);              return newLocation;     }      I am not sure what i am doing wrong, because the new locations are created in the middle of the sea :) Any idea? thank you


title: Is it possible to convert a jpg file into lyr

text: I am a beginner with GIS 10.0. IS there any way to convert a jpg file into a layer? I am working on a problem that needs to have the entire map of my city as a layer but I just have a jpg file of the map.


title: WMS application in Apache

text: To host a map application in real world, developed using GeoServer, we need a server with Tomcat/JSP rather than Apache, right? Correct me if I am wrong since my background is as GIS not development. Is there any options to host GeoServer in Apache Server or any similar applications like GeoServer that can connect with an Apache server. Because to maintain a good Apache Server is economical compared to Tomcat. Update: In GeoServer, as localhost, if we login and check Server Status, Memory Usage shows 376 MB. Therefore, if we plan to hire a host plan in JSP/Tomcat, the heap size offered should be greater than 376 MB, isn't it?, which will escalate the server plan cost. Need help on this aspect.


title: Add buffer only on one specific side of polygons

text: I have a shapefile with multipart geometry (multipolygons), namely communities (blue on the picture below). There are also lakes (white), but for these I don't have a geometry, they're just "empty". I want to conduct a point-in-polygon analysis with another point dataset. Each point should be assigned to a community or be deleted when it falls outside a community (e.g. in a lake). The problem is that there are many points for lake-bordering communities that are very close to the shore. Because of the strong generalization, these points get accidentally deleted even though they would actually belong to a community. ![enter image description here](http://i.stack.imgur.com/GsSY0.png) The only thing I can think about is somehow to extend the shore line of lake- bordering communities by about 100m, but I'm wondering whether there is some way of doing this automatically in QGIS. I have no clue though I would approach this problem.


title: Comparing choropleth maps statistically

text: I have several sets of choropleth maps showing ranges of poverty by census tract, each using a different poverty measure. All of the maps use the same categorical ranges to define poverty status. The maps clearly differ visually in the distribution of poverty, but I am having trouble finding a way to provide a statistical comparison between maps. I have tried using contingency tables and the comparison of chi squared values - however, this really only compares the counts of areas that fall into each category, but provides no indication of spatial coincidence. Are there any statistical measures to compare these types of maps both numerically and spatially? Edit: Examples of maps are below ![Measure 1 - Adjusted per capita income \(by standard score\)!\[\]\[1\]](http://i.stack.imgur.com/kw3kB.jpg) ![Measure 2 - Income Ratio \(by standard score\)](http://i.stack.imgur.com/2rcXC.jpg)


title: Where can i download a good World image in GeoTIFF?

text: I can either be from MERIS (ESA) or from MODIS (Nasa). I already reached these sites: http://www.unearthedoutdoors.net/global_data/true_marble/download http://www.eosnap.com/public/mosri/2009/world/world-mosaic-final-100.jpg But from the first I was not able to download, and from the second it is not available in tiff. Thanks


title: Sum up area left "underneath" polygons in QGIS?

text: I have an area shape layer and multiple other area shape layers. Now I would like to calculate how much space (square kilometers or percentage of area) of the polygons of this area layer is NOT covered by the other layers. Some of the other layers are overlapping so it isn't as easy as to substract these areas. Is this possible and if yes - does anyone here know how? best regards desputin


title: Does ArcGIS.com support dynamic map services (or only feature services and tiled services)?

text: Using ArcGIS Server it is possible to publish a map service, which has a URL like http://services.arcgisonline.com/ArcGIS/rest/services/Demographics/USA_1990-2000_Population_Change/MapServer/1 In this case, the "1" refers to the particular polygon layer. It's possible to create a featureLayer from this REST endpoint. When I publish a map to ArcGIS.com the options are "Feature Access" or "Tiled Mapping". This means that the REST endpoint has **/FeatureServer/** rather than **/MapServer/** in the URL. Am I correct in the assumption that ArcGIS.com doesn't support dynamic map services, and only supports feature services and tiled services?


title: "Invalid feature class type" for input feature while using "Locate Features Along Routes"

text: I am getting ERROR 000810 (Invalid feature class type) when I add my stream shapefile as an in_routes. That is really weird to me. ![enter image description here](http://i.stack.imgur.com/7L2nl.jpg)


title: Calculating point proximity in ArcMap 10.1

text: I have a geocoded point file which (theoretically) contains one point for every zip code in the U.S. I want to create a smaller list of zip codes from this file. The list should include only those zip codes that are 100 miles from one another. In other words, if I were to look at this on a map, I would only see one point for every 100 miles, and any points in between would not be included. I'm having trouble getting this part of my project started in ArcMap 10.1. How would I go about creating this output? I thought that I could use the Near or Point Distance tools, but I always run into trouble when trying to use the same file as both my inputs. Also, I thought that I could drape a 100x100 mile grid on top of the point file, and then grab the closest one to the center of each panel, but I'm not sure how to set that up either. Any ideas?


title: QGis doesn't start anymore on Debian Jessie

text: I'm running Debian Jessie (Testing) on my laptop. I installed QGis 2.0.1 from QGis stable repos almost a month ago and had no problems till yesterday; now QGis doesn't start, only splash screen starts and then terminate. The output is (I'm italian):               andrea@debian:~$ qgis     Warning: loading of qt translation failed [/usr/share/qt4/translations/qt_it_IT]     QGIS died on signal 11Annullato     andrea@debian:~$       and when I post the command as root the output is:               root@debian:/home/andrea# qgis          (process:12038): GConf-WARNING **: Client failed to connect to the D-BUS daemon:     Did not receive a reply. Possible causes include: the remote application did not send a reply, the message bus security policy blocked the reply, the reply timeout expired, or the network connection was broken.          (qgis.bin:12038): GConf-WARNING **: Client failed to connect to the D-BUS daemon:     Did not receive a reply. Possible causes include: the remote application did not send a reply, the message bus security policy blocked the reply, the reply timeout expired, or the network connection was broken.     Warning: Qt: Session management error: None of the authentication protocols specified are supported     Warning: loading of qt translation failed [/usr/share/qt4/translations/qt_it_IT]          (qgis.bin:12038): GConf-WARNING **: Client failed to connect to the D-BUS daemon:     Did not receive a reply. Possible causes include: the remote application did not send a reply, the message bus security policy blocked the reply, the reply timeout expired, or the network connection was broken.     QGIS died on signal 11Annullato      2 days ago some libqt4-* were updated and QGis depends on them, but I don't how to fix it


title: QGIS 1.8 in Linuxmint 14 (Ubuntu 12.10)

text: I have Linux Mint 14 KDE which is based on Ubuntu 12.10 Quantal. In the QGIS download page, there's no "Quantal" repositiries. I'm trying to use the "Precise" repositories, but I'm not allowed to install the latest QGIS version (1.8) and I received many messages of broken packages. Any suggestion? Thanks! Lucas


title: Why doesn't QGIS accept my manually altered geospatial column?

text: I have a table with a geospatial column with srid 28992 but noticed that if i try to show it in QGIS (v1.9) on any openlayer map (through the OpenLayers plugin) it is off by a few hundred meters. Then i converted it to 900913 (remembering something that with an old plugin this projection was required) and then it is shown on the correct spot. However when i look in the metadata of this layer in QGIS (right-click -> Properties -> Metadata) everything is blank and when i do Zoom to Layer Extent the map jumps to water near Africa (suggesting it doesn't understand the projection). Somehow the metadata gets lost to QGIS during my manual conversion. The script i use for the manual conversion is:               alter table <table_name> drop constraint "enforce_srid_<column_name>";     update <table_name> set <column_name> = ST_Transform(<column_name>,900913);     alter table <table_name> add constraint "enforce_srid_<column_name>" CHECK (st_srid(<column_name>) = 900913);     update geometry_columns set srid=900913 where f_geometry_column = '<column_name>';      The result of this seems no different (except that now i have an extra column) than if i use:               SELECT AddGeometryColumn('<table_name>','<new_column_name>',900913, 'LINESTRING', 2);     UPDATE <table_name> SET <new_column_name> = ST_Transform(<column_name>,900913);      Does any of you know what god particle i am looking for?


title: Problem with GDAL_WMS QGIS connecting to an ArcGIS Server?

text: I've been working, trying to get QGIS ( _v2.0 32bit_ ) to connect to an ArcGIS Server service using the GDAL_WMS xml file. This is the service I'm trying to connect to: http://geonb.snb.ca/ArcGIS/rest/services/GeoNB_Basemap_Imagery/MapServer/?SERVICE=WMS&REQUEST=GetCapabilities This is the xml file I currently have:               <GDAL_WMS>        <Service name="TMS">             <Version>1.1.1</Version>           <ServerUrl>http://geonb.snb.ca/ArcGIS/rest/services/GeoNB_Basemap_Imagery/MapServer/tile/${z}/${y}/${x}.jpg</ServerUrl>             <Layers>Bathurst_ortho_final_graphic,SDEDO.FTON_ORTHOS_Clip2,MiramichiOrthos,SDEDO.snbortho,Year of Photography</Layers>             <ImageFormat>image/jpeg</ImageFormat>        </Service>        <DataWindow>             <UpperLeftX>2282341.95</UpperLeftX>             <UpperLeftY>7698181.3</UpperLeftY>             <LowerRightX>2737347.05</LowerRightX>             <LowerRightY>7252058.7</LowerRightY>             <TileLevel>13</TileLevel>             <TileCountX>1</TileCountX>             <TileCountY>1</TileCountY>             <YOrigin>top</YOrigin>         </DataWindow>         <Projection>EPSG:2953</Projection>         <BlockSizeX>512</BlockSizeX>         <BlockSizeY>512</BlockSizeY>         <BandsCount>3</BandsCount>         <MaxConnections>10</MaxConnections>         <Cache />     </GDAL_WMS>      And this is what I get: ![http://i.imgur.com/yCJTaFO.png](http://i.imgur.com/yCJTaFO.png) This is not a WMS issue. The service doesn't support WMS connections and that's the reason for using the XML with GDAL. GDAL then creates a "WMS" like service on top of the images using the XML. It basically fetches images and puts them in the right spot. What I think is happening is it misinterpreting what tiles to fetch and how to stretch them. (Maybe) Here's an example of a tile fetched directly from the service: http://geonb.snb.ca/ArcGIS/rest/services/GeoNB_Basemap_Imagery/MapServer/tile/0/64/64.jpg ![http://geonb.snb.ca/ArcGIS/rest/services/GeoNB_Basemap_Imagery/MapServer/tile/0/64/64.jpg](http://geonb.snb.ca/ArcGIS/rest/services/GeoNB_Basemap_Imagery/MapServer/tile/0/64/64.jpg) This all works so there's something wrong with my XML. The documentation for the XML is here(www.gdal.org/frmt_wms.html), but it's not very descriptive. The specific section that relates to what I want to accomplish is the TMS service. Any help would be greatly appreciated!


title: How to order points of a polygon?

text: I want to create a polygon, which represents the administrative boundary of a city. It is specified by a relation, which consists of ways               <member type="way" ref="25657490" role="outer"/>     <member type="way" ref="209124836" role="outer"/>     <member type="way" ref="119880388" role="outer"/>     <member type="way" ref="119880389" role="outer"/>     <member type="way" ref="119983752" role="outer"/>     <member type="way" ref="136380499" role="outer"/>     <member type="way" ref="120024176" role="outer"/>     <member type="way" ref="209124834" role="outer"/>     <member type="way" ref="77909008" role="outer"/>      In order to get all buildings inside that boundary using the method described here, I need to convert these ways to an ordered list of points such that every point is connected to its successor. **Example 1** I have way 119880388, whose last point is `132.100403,43.3959161`, hence, it must be followed by a way, which starts there (way 209124836 in this case).               119880388;132.1469949,43.504065|132.1439962,43.4906939|132.1400282,43.4747441|132.1402881,43.4739524|132.1397559,43.4733796|132.139471,43.4710427|132.1324332,43.4650797|132.1321001,43.4625553|132.1273123,43.4602785|132.1268996,43.4611024|132.1262186,43.4615517|132.1250473,43.4606504|132.1253362,43.4590926|132.1199706,43.4567557|132.1181958,43.4534|132.1165449,43.4536997|132.1151003,43.4558569|132.1122111,43.4550479|132.1120048,43.4568456|132.1106349,43.456653|132.1074084,43.4565518|132.1045214,43.4568074|132.1018878,43.457799|132.1045918,43.45418|132.1045214,43.4533315|132.1049016,43.4526109|132.0974329,43.4515698|132.0970943,43.4510349|132.0962379,43.4508758|132.0956005,43.4499071|132.0965565,43.4495456|132.100208,43.4491561|132.096659,43.4434507|132.098349,43.4423872|132.0988278,43.440833|132.0981518,43.4397695|132.099194,43.4388083|132.0939234,43.4347791|132.0934023,43.4339098|132.0915574,43.4311281|132.0898252,43.4291849|132.0888394,43.4281416|132.0866256,43.4271254|132.0844248,43.4271615|132.0837273,43.4275697|132.0837273,43.4262298|132.0852483,43.4242149|132.0850934,43.4232227|132.0822767,43.4233557|132.0820373,43.4226294|132.0811501,43.4219032|132.083854,43.4190391|132.0911773,43.4176683|132.0920223,43.4173205|132.0863045,43.4153769|132.0873467,43.4135355|132.0895155,43.4113666|132.0937405,43.4059441|132.0943038,43.4036522|132.0952051,43.4011965|132.0967132,43.3988428|132.0991074,43.3975125|132.100403,43.3959161|          209124836;132.100403,43.3959161|132.0973611,43.3915157|132.0945444,43.3876676|132.0947698,43.3862962|132.0962344,43.3844129|132.0891928,43.3825706|132.0861226,43.3815675|132.0843482,43.3835532|132.0826018,43.3824273|132.0823202,43.381158|132.0812181,43.3793895|132.0799157,43.3777247|132.0783487,43.3754872|132.0775037,43.3725185|132.076884,43.3713105|132.0734759,43.3689559|132.0715324,43.3696725|132.0686313,43.3693654|132.0671666,43.3700206|132.0670821,43.3692425|132.0649696,43.3681778|132.0646316,43.3658435|132.0619981,43.3665295|132.0613643,43.3651985|132.0623502,43.3648913|132.0608851,43.3614618|132.0596304,43.360689|132.0584397,43.359838|132.0568436,43.3605808|132.0524573,43.3620336|132.0516256,43.3627178|132.0467659,43.3583303|132.0505899,43.3557236|132.0502911,43.3542031|132.0457676,43.3518866|132.0457395,43.350668|132.0445706,43.3502788|132.0457676,43.3470836|132.0464436,43.3469607|132.0466971,43.3455064|132.046345,43.3446358|132.0472089,43.342042|132.0438833,43.3408081|132.043113,43.3408202|132.0429138,43.3413996|132.0405483,43.3413749|132.0393135,43.3412952|132.0381384,43.3409041|132.0377749,43.3393052|132.0369876,43.3385224|132.0374383,43.3363304|132.0392269,43.3357363|132.0392019,43.3313722|      **Example 2** There are ways 25657490 and 77909008, which both end in point `132.0638277,43.3156504`. Hence the correct order is `25657490, 77909008'` where `77909008'` is the inverted list of points of way 77909008 (`132.0638277,43.3156504; 132.0641092,43.3159; ...`).               25657490;132.0392019,43.3313722|132.0440197,43.3311274|132.0473289,43.3308512|132.0508932,43.3302733|132.0528244,43.3296177|132.0532452,43.3290625|132.0534252,43.3291183|132.0536528,43.3295151|132.0540176,43.3294683|132.0547127,43.3293056|132.0556996,43.3292972|132.0561718,43.3297738|132.0578026,43.3296177|132.0592617,43.3291807|132.0595621,43.3283378|132.0597767,43.327651|132.0608925,43.3269642|132.061665,43.3273076|132.0623516,43.3275573|132.0628666,43.3273388|132.0628666,43.3264023|132.063442,43.3254975|132.063507,43.3251849|132.0634183,43.3248517|132.0633476,43.3245859|132.0636823,43.3243424|132.0639065,43.3243493|132.0640514,43.3245921|132.0643604,43.3245921|132.0643422,43.3244722|132.0645664,43.3243798|132.065056,43.3230585|132.0654247,43.322644|132.065562,43.3219696|132.0655105,43.3216699|132.0657594,43.3212203|132.0653989,43.3213951|132.0650127,43.3217448|132.0648496,43.3223005|132.0648239,43.3225503|132.064223,43.3218947|132.0638111,43.3208893|132.0635333,43.3200936|132.0635103,43.3191278|132.0638536,43.3185658|132.0647549,43.3177852|132.0635103,43.315943|132.0638277,43.3156504|          77909008;132.1112088,43.3053453|132.1109191,43.305439|132.1103183,43.3058684|132.1100501,43.3059386|132.1093419,43.3059308|132.1090415,43.305962|132.1088377,43.3061182|132.10843,43.3068755|132.1073678,43.3069848|132.1070567,43.3069302|132.1067563,43.3067662|132.106134,43.3064617|132.1052113,43.3061963|132.1050397,43.306165|132.1042136,43.3064695|132.1039132,43.3070238|132.1035055,43.3071722|132.1032372,43.3073439|132.1030763,43.3075704|132.1030012,43.3078358|132.1030012,43.308312|132.1027223,43.3084682|132.1024004,43.3085853|132.1019927,43.3086399|132.1014992,43.308679|132.1011237,43.3088351|132.1002761,43.3091708|132.0999757,43.3092333|132.0997611,43.3094753|132.0996967,43.3096704|132.0989886,43.3104043|132.0984736,43.3104589|132.0975939,43.31074|132.0971647,43.310982|132.0966712,43.3111225|132.0953623,43.311185|132.09474,43.3111225|132.0941177,43.311021|132.0936886,43.311021|132.0932594,43.3111147|132.0930556,43.3112474|132.092959,43.3114192|132.0929054,43.3116143|132.0927337,43.3117549|132.0923475,43.311872|132.0918754,43.3122076|132.0917467,43.3122935|132.091339,43.312395|132.0909313,43.312395|132.0907274,43.3124262|132.0904485,43.3124184|132.0902983,43.3124262|132.0899764,43.3123638|132.089912,43.3124418|132.089397,43.3125902|132.0891181,43.312598|132.088528,43.3124965|132.08841,43.3125121|132.0878521,43.3127229|132.0876053,43.3127385|132.0872942,43.3126916|132.0869831,43.3127229|132.0861677,43.3130273|132.0860496,43.3132147|132.0857278,43.3133474|132.0855561,43.3133318|132.0851055,43.3135113|132.0848158,43.3137377|132.0844403,43.3138626|132.083818,43.3142061|132.0836786,43.3143388|132.0836035,43.3145418|132.0836464,43.3146589|132.0834211,43.3148228|132.0833996,43.3149945|132.0830027,43.3151585|132.0828095,43.315338|132.0824555,43.3153458|132.0820907,43.3152209|132.0817903,43.3152756|132.0810286,43.3153146|132.0808354,43.315377|132.0805994,43.3155566|132.0804921,43.31558|132.0803741,43.3155644|132.0802883,43.3155019|132.0801917,43.3153146|132.0801059,43.3152443|132.0799664,43.3152131|132.0797733,43.3152365|132.0796231,43.3152834|132.0794192,43.3154239|132.0792583,43.3154629|132.0791295,43.3154629|132.0790115,43.3154083|132.078797,43.3149243|132.0786575,43.3148618|132.0785073,43.314893|132.0784,43.3149555|132.0781103,43.3150179|132.0777455,43.3151428|132.0774344,43.3153146|132.0771018,43.3155253|132.076619,43.3156893|132.0762328,43.3159078|132.0757607,43.3160093|132.0751277,43.3160562|132.0745161,43.316142|132.0739261,43.3163996|132.0717588,43.31711|132.0699349,43.3178203|132.0687011,43.3184448|132.0679286,43.3187726|132.0674244,43.3189053|132.0670596,43.3189522|132.0666412,43.3189522|132.0663408,43.3188585|132.0661369,43.3186868|132.0660511,43.3184838|132.0660082,43.3182106|132.0658043,43.3178359|132.0650748,43.3172427|132.0641092,43.3159|132.0638277,43.3156504|      Are there any ready-made routines (ideally in Java), which do that (append points of ways to each other, inverting them where necessary).


title: Errors when building OpenGeo Suite with Maven 3.1.1

text: I am trying to build the OpenGeo(Boundless) Suite with Maven 3.1.1 I am following the instructions on GitHub. When building I am getting the next error : GitHub : https://github.com/boundlessgeo/suite               % mvn clean install          [ERROR] Failed to execute goal org.apache.maven.plugins:maven-antrun-plugin:1.7:     run (clean) on project data_dir: An Ant BuildException has occured: The following error occurred while executing this line:     [ERROR] java.io.FileNotFoundException: D:\boundless\suite-master\suite-master\geoserver\data_dir\externals\suite-data\default\build.xml (The system cannot find the file specified)     [ERROR] around Ant part ...<ant dir="externals/suite-data/default" target="clean"/>... @ 4:59 in D:\boundless\suite-master\suite-master\geoserver\data_dir\target\antrun\build-main.xml      When building the GeoExplorer application on itself I get the next error : GitHub : https://github.com/boundlessgeo/suite/tree/master/geoexplorer               % mvn clean install          [ERROR] Failed to execute goal org.apache.maven.plugins:maven-antrun-plugin:1.7:     run (default) on project geoexplorer: An Ant BuildException has occured: Java returned: -1     [ERROR] around Ant part ...<java fork="true" classname="org.ringojs.tools.launcher.Main" failonerror="true">... @ 40:84 in D:\boundless\suite-master\suite-master\geoexplorer\target\antrun\build-main.xml      Does anyone have similar experiences ? Is this problem caused by the Maven version ?


title: arcpy array: trouble with reassigning inf to 0

text: I've converted rasters to arrays via arcpy.RasterToNumPyArray() then divided them by 1. Some resulting values are zero and I need to reassign them after the division to zero. The array types are float64. I've tried quite a bit of versions of the following code (from other SE searches assigning zero to infinity though) without luck. Any advice?               array = array[array == float("inf")] = 0          TypeError: 'float' object does not support item assignment


title: Publishing Multiple Shapefiles

text: I'm using OpenGeo Suite 3.1.I have stored large number of shapefiles(100) in a store()Directory of spatial files. How to Publish all the files at a time in geoserver?


title: Why does vector grid creation crash?

text: Im trying to create a vector grid for a municipality. Do you know why does it crash? The scale I'm using is 1:1000000. Any tips? Thanks! ![enter image description here](http://i.stack.imgur.com/5l4Y1.png)


title: WMS orthophotos not visible on screen

text: Using QGIS on Windows XP based desktop. I have been loading orthophotos from a WMS server up to yesterday morning. Today QGIS correctly connects to the server and apparently loads some bytes, however I only get overlapping yellow frames and the codes of the orthophotos on a white screen. No error messages appear. What I have done up to now: I rebooted the pc, I checked the WMS URL (even if QGIS connects to the server), I have installed the latest version of QGIS (now 2.0.1). Any clue? I can correctly view images and videos on my PC.


title: Can Esri ASCII Input files be Aligned with Track Heading in QGIS?

text: Using QGIS to input Esri ASCII data files, I uploaded many small files, each file representing a second of data: 16 sensors X 25 lines. Each file uploaded creates a small rectangle on my larger map layer. The default MAP Orientation of Esri ascii files is with the Top edge of the rectangle facing North and Right/Left sides pointing Directly North. I want to ROTATE this MAP rectangle, during INPUT, so the Right/Left edges would be at an Angle to North. Is there a command line option that allows this?


title: Problem in removing dangles in road network

text: Im using qgis2.0 , In a given road network i have to remove dangles, for that what should be the tolerance to be added generaly. please help me resolve this issue


title: Last moment a layer in Geoserver was viewed

text: My company publishes >1,000 layers using several geoservers. We are publishing them on several different websites. Because there are so many layers, it is impossible for us to know which layer is being used and which one is not. Therefore, I was wondering whether it is possible to view the last time a layer in Geoserver was requested. Is there any option to know this?


title: Interpolating vertically in 3D analyst

text: Here is my problem: I am working on mapping pollutant concentrations in a pond from point-sampled data. There are 50 sampling locations in the pond (randomly dispersed), and at each location there are three sampling events (=150 total); one at the surface, one in the middle, one at the bottom. I have created a TIN extent based on depth of pond, and now I would like to "fill" the pond by interpolating volumes of pollutant concentrations from the 150 sample points. I am trying to do this by applying a spline interpolation through the three sampling events at each of the 50 locations, in order to produce 50 interpolated lines. From there i could apply kriging for each pixel down through the pond profile (although this way does ignore important spatial correlation). I am stuck at the spline interpolation of the lines. Is there a way to interpolate between points when the X,Y coordinates are all the same, but the Z coordinate changes? In other words, how does one do vertical interpolation in 3D analyst? Thank you, Michael


title: Grid question for Fire Station MAP

text: How can I generate a grid that would have a quadrant like naming? Or is there a way just to rename quadrants![enter image description here](http://i.stack.imgur.com/0Wpvx.jpg) through data driven pages. Sort of Newbish but basically creating a new firestation map book and the typical row names (as see attached), does no suit us too well. I would like to see (2nd attachement with 10) more of a quadrant. How would I generate a data driven template/grid with that? Thanks! ![enter image description here](http://i.stack.imgur.com/CU4Tm.jpg)


title: How to normalize multiple PNG rasters with differing colormaps

text: I have a series of adjacent rasters which comprise a digital elevation model. Unfortunately they have been given to me as PNG files, rather than Esri Grids. The only key that I have to the actual elevation values is a separate legend PNG. Each PNG file has an inbuilt colormap file, and they display correctly when drawn in ArcMap using the automated Colormap symbology option. The red line in this screenshot shows the boundary between two adjacent images: ![enter image description here](http://i.stack.imgur.com/bP8KR.png) However, Identifying within the dark green sections of each image shows that while the colors match (a,r,g,b = 255,0,225,0 on both the upper and lower images) the Color Index values do not match. This means that if I try to Mosaic these images, there are discrepancies between the values. I'm using the Mosaic To New Raster tool with the Match Colormap option, but this seems to match the Color Index values rather than the colors(?). Exporting the colormaps and comparing them in a text editor confirms that the Color Index vs a,r,g,b values do not match between the images, and that there isn't a simple shift between the values in each image: ![enter image description here](http://i.stack.imgur.com/bWAV0.png) Does anyone have any suggestions on how I can normalize (if that's the correct term) these adjacent images so they are using the same values? I can then mosaic and reclassify this single image to match with the legend. (A cumbersome workaround would be to reclassify each individual image based on the legend, but I'd like to avoid this as it'll be hard to scale). Thanks.


title: Can QGIS alone can be used without installing grass?

text: can QGIS alone can be used without installing grass? i want to use functions of grass and i am not interested to create mapset and all. just with raster and vector files, i want to study temporal data. will it be possible with qgis alone? what is the procedure.


title: Trimming Specific Text Values in QGIS Field Calculator

text: Sorry if this is a dumb question, but I am wondering if there is a way to copy the text from one field to a new field, but eliminate specific words. For example:               **Old Field**                         **New Field**          Killbear Provincial Park               Killbear     Quetico Provincial Park                Quetico     La Cloche Provincial Park              La Cloche     Wildcat Hill Provincial Park           Wildcat Hill      I thought of trying a trim function, but not all the values have the same start location or trim length. Thanks...


title: Using GeoWebCache to cache google maps tiles

text: I'm new to GIS and I would like to know if anybody could give me a quick example of using geowebcache with google maps tiles. Thanks.


title: Help in GRASS/QGIS/Python

text: I have Idea about ArcGIS, and a bit of QGIS and zero programming knowledge, which I am trying for thesedays. I saw GRASS GIS and have started working with shapefiles, which are not being loaded... :( I saw the toolbox and could not find v.in.shape tool. Using file management -> import vector into grass-> v.in.ogr.qgis I cod load the file but could not go farther. I now wanted to access this via python console. Any guide on how should I carry on with my task through which I may be able to use the tools. I can get good help in programming. I tried to get my shapefile read in python console in QGIS but could not. I am very basic in this so please give easy solutions. thanks a lot.


title: How to convert accuracy of DEM from 5 decimal places to 2

text: I had a csv dataset comprising of regular points containing elevation "z-values". I interpolated this using ESRI´s spatial analyst using the IDW method to produce a continuous raster surface. I set the resolution at 12,5 meters, because the points were 12,5 meters apart. I now have a raster image which gives me values with 5 decimal places. The original data only had 2 decimal places. My question is...should I round this up to 2 decimal places, because having this level of accuracy doesn´t reflect the absolute accuracy of the image. If I should round up the values, how/with which tool would I do it? thanks, Rob


title: Reordering Category Layers in QGIS

text: ![enter image description here](http://i.stack.imgur.com/GlaOa.jpg)I was wondering if there is a way to reorganize the layers when applying a style based on a categorizing value? I have many values, but I would like to rearrange them alphabetically. In ArcMap there are arrows to promote, or demote the values in the list, but I do not see anything like this in QGIS 1.8.0. Any help would be appreciated. Thanks...


title: Grass 7 .las data import problem

text: I compiled Grass7 on Ubuntu 12.10, with liblas in order to run r.in.lidar and v.in.lidar. I now want to import .las data; the option in file -> import raster data -> LAS LIDAR points import is given and enabled, but when I select it nothing happens. Any idea where the bug could be?


title: What GIS Data should Amazon add to its Public Data Sets?

text: Amazon says: > If you have a public domain or non-proprietary data set that you think is > useful and interesting to the AWS community, please submit a request below > and the AWS team will review your submission and get back to you. Typically > the data sets in the repository are between 1 GB to 1 TB in size (based on > the Amazon EBS volume limit), but we can work with you to host larger data > sets as well. You must have the right to make the data freely available. Given that Amazon is a gold sponsor at the Esri Dev Summit, maybe this would be a good time for the community to recommend data sets for Amazon to put in their cloud. Having the data on Amazon allows fast access by things like Server Object Extensions and geoprocessing services running on EC2. Currently I just see 3 datasets for geography. **Update** I dropped by Amazon's booth at the Dev Summit and spoke a bit with an AWS rep. He said they are still sorting out how they will keep public dataset current. He said there has been a lot of interest from the Human Genome community. He said to expect something in last half of 2011 regarding changes to community data sets that will make it easier for GIS. I also visited with Microsoft. While Azure looks promising, it doesn't look like you'd be able to do server-side geometric network traversal, as you would with AWS.


title: Google Fusion Tables vs Amazon Web Services: which is better for GIS data?

text: What are the pros/cons for storing storing spatial data in Google Fusion Tables versus Amazon Web Services? Factors would include: Cost, performance, ease of update.