AIS Data Transformed or CVS Format

I am looking for AIS Shipping Route Data. I have not had much luck finding it in either SHP format or even in CVS format does anyone have a real good source of finding this type of data? To add to my problems I need only a specific area of coverage? It is starting to look like I may need to purchase the Data. Thank You


Performance of google map tile creation processes

I know that question is rather vague but please bear w/me. I'm trying to get an idea of what sort of product performance - specifically timing - people have seen for various methodologies they have used to create google/bing map tiles. There are a slew of methods for how to do this (e.g. gdal2tiles, FME, maptiler, etc). An initial attempt at simply taking a large PNG and creating tiles using imagemagick, on a pretty decent linux server, yielded some pretty long processing times and so I wanted to see what other people are using in production. New tiles would need to be generated at least daily and so turnaround time on this is pretty critical. The only real requirement is that it can run on a linux server. Obviously, free is better but I don't want to restrict myself to that. The input can be raw gridded/raster data or a large image. The output needs to be image tiles capable of being used as-is in google or bing maps. Just for the sake of comparison, I'll say that the timings should be for google map's zoom level 7. I appreciate everyone's help and again I want to apologize for how vague this question probably seems. UPDATE: As far as the inputs, I currently have multiple (raw) data sources in various formats: netCDF, GRIB, GRIB2. In addition to the raw data itself, I also have the ability to generate really large images of that data which could then be sliced/tiled. Ideally, I would just be chopping the image up but I am willing to try whatever will get me the fastest results.


How to customize widgets in ArcGIS Viewer for Flex 3.4 to support French?

How can I customize my widgets in flex viewer such that the labels will be in French? I work with flash builder 4.6.


WMS - get point data in computer readable form (other than map tile)

is it possible to get a data from some map services on questions like   * **is there a water on [X,Y]** (i.e. blue colour on map)   * is there a dessert on [X,Y] (i.e. yellow colour on map)   * similarly for any available layer In language of computers, to know, **where on some location (point,radius) is some layer present** and more ideally, what are its properties (river=1,lake=0,waterfall=0,...) I get get these data visually just by looking on some map tile presented on screen, but this is hard for computer calculations. I have read WMS reference http://docs.geoserver.org/stable/en/user/services/wms/reference.html#getfeatureinfo but I am not any more clever from that. Thanks


python equivalent of spatstat's spatial point process functions

Is there a free python package (i.e. PySAL?) that implements the standard summary statistics for spatial point processes, e.g. the F, G, J, and K functions? I'm happy to use R (i.e. spatstat) if that's the best option, but I want to try out some other things I already have implemented in python. The implementations don't have to be nearly as exhaustive as R, just looking for something to get started.


How to get names of all layers in a store in GeoServer using OpenLayers

I have strict requirement to retrieve names of all layers in a GeoServer store using OpenLayers. I've gone through the post here, but it doesn't answer my question. I've also gone through GeoServer's REST API, which would be ideal if I'm using another language, say C#, Java, or PHP. But I want to do this entirely in JavaScript. **UPDATE:** From what I can see in the `WFS GetCapabilities` document, I can retrieve names of all layers by sending a WFS GetCapabilities request using filters. But I'm not sure if its possible to incorporate filters in a `GetCapabilities` request. Neither do I have any idea how such a filter will look like. Any advice on how I can go about achieving this will be highly appreciated. Any alternative approach will also be appreciated. **UPDATE:** I believe I've found a solution here; "WFS GetCapabilities requests can be filtered to only return layers corresponding to a particular namespace. To do this, add the following code to your request: namespace=`<namespace>`". I hope this assists someone else.


To get attribute data of polygon layer by giving Lat/Lon

Is anyone having an idea about how to get attribute data of a polygon layer by giving Coordinates with in that polygon. Here i have tried with these codes.. `SELECT * FROM kar_bbmp_198 WHERE ST_Contains(geom, ST_GeomFromEWKT('SRID=4326;POINT(779322 1440590)'))` `SELECT * FROM public.kar_bbmp_198 WHERE ST_Within(geom,ST_PointFromText('POINT(779322 1440590)',32643));` But I am getting an error message like `ERROR: Operation on mixed SRID geometries ********** Error ********** ERROR: Operation on mixed SRID geometries SQL state: XX000` Any idea... I am using PostGIS 2.0 and PostgreSQL 9.1 Thanks in advance...


Drawing polyline that grows as the car moves by taking data from MySql and PHP in Google API

I wanted to show the polyline on google Map as car moves with animation (like in this site: http://econym.org.uk/gmap/example_cartrip2.htm) by taking data from MySQL and PHP. For that i refered the code from Google API tutorial Polyline. I have also taken data for this from my SQL which is as follow: .html               function load() {         var point;                         var flightPlanCoordinates=new Array();                 var map = new google.maps.Map(document.getElementById("map"), {                     center: new google.maps.LatLng(18.33, 73.55),                     zoom: 7,                     mapTypeId: google.maps.MapTypeId.ROADMAP,                     mapTypeControl: false,                     mapTypeControlOptions: {                     style: google.maps.MapTypeControlStyle.HORIZONTAL_BAR                      },                     navigationControl: true,                     navigationControlOptions: {                         style: google.maps.NavigationControlStyle.SMALL                     }                 });                      downloadUrl("xmltry.php", function(data) {             var xml = data.responseXML;                  var markers = xml.documentElement.getElementsByTagName("marker");             for (var i = 0; i < markers.length; i++) {                 point = new google.maps.LatLng(                   parseFloat(markers[i].getAttribute("lat")),                   parseFloat(markers[i].getAttribute("lon")));                 flightPlanCoordinates[i]=point;               }               var flightPath = new google.maps.Polyline({         path: flightPlanCoordinates,         strokeColor: "#FF0000",         strokeOpacity: 1.0,         strokeWeight: 2       });            flightPath.setMap(map);                });                  }                  function downloadUrl(url, callback) {           var request = window.ActiveXObject ?               new ActiveXObject('Microsoft.XMLHTTP') :               new XMLHttpRequest;                request.onreadystatechange = function() {             if (request.readyState == 4) {               request.onreadystatechange = doNothing;               callback(request, request.status);             }           };                request.open('GET', url, true);           request.send(null);         }         function doNothing() {}      .php file:               <?php     //require("phpsqlajax_dbinfo.php");     include 'common_db.inc';     // Start XML file, create parent node          $dom = new DOMDocument("1.0");     $node = $dom->createElement("markers");     $parnode = $dom->appendChild($node);           // Opens a connection to a MySQL server          $connection=mysql_connect ("$dbserver", "$dbuser", "$dbpass");     if (!$connection) {  die('Not connected : ' . mysql_error());}           // Set the active MySQL database          $db_selected = mysql_select_db($dbname, $connection);     if (!$db_selected) {       die ('Can\'t use db : ' . mysql_error());     }           // Select all the rows in the markers table          $query = "SELECT  `lat`, `lon` FROM current_location Where Device_ID='HGS1005'";     $result = mysql_query($query);     if (!$result) {         die('Invalid query: ' . mysql_error());     }           header("Content-type: text/xml");           // Iterate through the rows, adding XML nodes for each          while ($row = @mysql_fetch_assoc($result)){         // ADD TO XML DOCUMENT NODE         $node = $dom->createElement("marker");         $newnode = $parnode->appendChild($node);         $newnode->setAttribute("lat", $row['lat']);         $newnode->setAttribute("lon", $row['lon']);         }           echo $dom->saveXML();          ?>      But now i wanted to show polyline with animation as show in 1st references site.Can any one suggest me some references or help me out in same. I am getting this result which is correct but want animation in it as in http://econym.org.uk/gmap/example_cartrip2.htm ![enter image description here](http://i.stack.imgur.com/gO1sS.png) Hope this time my title and description part is understood.


Modify WFS with OpenLayers

I want to edit object on WFS layer but have one problem. I use `OpenLayers.Control.ModifyFeature(myVecLayer)` to editing.By the logic i do this.   1. Press button and create WFS layer with chosen object   2. Add `OpenLayers.Control.ModifyFeature` to map controls and activate this   3. Edit object   4. Press button again and remove WFS layer with my object, deactivate `OpenLayers.Control.ModifyFeature` control.   5. Select next object (or same object) and do "1", "2" And object not editing.   But if I do not remove layer after editing, but just add new WFS layer to map `OpenLayers.Control.ModifyFeature` control works My code:               function ff(){     if(counter==0){         //WFS             var myStyle = OpenLayers.Util.extend({},                                     OpenLayers.Feature.Vector.style['default']);                                 myStyle.strokeWidth = 1.5;                                 myStyle.strokeColor = "#ff0000";                                 myStyle.fillOpacity = 0.1;                                 OpenLayers.ProxyHost = "proxy.cgi?url=";                                 var renderer = OpenLayers.Util.getParameters(window.location.href).renderer;                                 renderer = (renderer) ? [renderer] : OpenLayers.Layer.Vector.prototype.renderers;                                            myVecLayer = new OpenLayers.Layer.Vector("Редактируемый участок", {                                     strategies: [new OpenLayers.Strategy.BBOX(), new OpenLayers.Strategy.Save()],                                     style: myStyle,                                     protocol: new OpenLayers.Protocol.WFS({                                         version: "1.0.0",                                         srsName:"EPSG:900913",                                         url: "http://localhost:8080/geoserver/wfs",                                         featureType: "filedata",                                         featureNS: "http://www.opengeospatial.net/cite",                                         defaultFilter : new OpenLayers.Filter.Comparison({                                             type: OpenLayers.Filter.Comparison.EQUAL_TO,                                             property: 'id',                                             value: "162",                                         })                                     }),                                     renderers: renderer                                 })                                   map.addLayers([myVecLayer]);                                           if(ppp==null){                                     if (edit==null){                                         edit = new OpenLayers.Control.ModifyFeature(myVecLayer);                                         map.addControl(edit);                                         ppp=map.getControlsByClass("OpenLayers.Control.ModifyFeature");                                     }                                          alert(ppp);                                 }                             //  ppp.autoActivate=true;                                 for(var i=0;i<ppp.length;i++){                                     ppp[i].activate();                                     alert("2 "+ppp[i].active);                                 }                                      counter=1;                          }      Its will nice if u help.


Find intersection points of two polyline features

I'm currently working with street intersections trying to find point(or points) where 2 given streets intersect(!disjoin). Streets are all Polyline features, and I have to work with 2 features at given time. I think this way it can be determined if they intersect(haven't yet tested it),                   public List<IPoint> RoadIntersections(IPolyline poly1, IPolyline poly2)         {             IRelationalOperator ro = null;             ro = poly1 as IRelationalOperator;                  ISegmentCollection sc = new PolylineClass() as ISegmentCollection;             sc.AddSegment((ISegment)poly2);                  if (ro.Disjoint(sc as IGeometry) == false) //if they intersect             {                 //TODO: find intersection points             }          ...      but what I would like to do is to find point (or points) where those two polylines intersect. Any help would be appreciated.


style loses to default(orange) with zoom in and out in openlayers

i was doing a project using openlayer and did some code to load wfs vector layer with STYLE on loading and during mouse hover which worked well .but when it zoomed in or out for the first time it loses the STYLE defined during onload TO default orange style. please give me some insights on how to maintain the STYLE when zooming i think this is my problem (http://trac.osgeo.org/openlayers/ticket/2281) regards , Zola


How do I catch PostGIS/PostgreSQL RAISE messages in QGIS?

I have a trigger set up in PostGIS/PostgreSQL which will show a `RAISE WARNING` when an inserted or updated polygon is within a certain distance of some predefined planning constraints. I want this warning to be sent to QGIS to alert a user to the proximity warning, but with a `WARNING` this does not occur. On the other hand, using `RAISE EXCEPTION` shows an error in the QGIS message bar but doesn't allow the insert/update to be committed (which I want - it should be a warning, not an error). Example of the messages I want to be sent to QGIS:               CREATE TABLE input_table (         id serial,         geom geometry                  );          CREATE TABLE listed_buildings (         id serial,         geom geometry                  );          CREATE OR REPLACE FUNCTION check_constraint_proximity()       RETURNS trigger AS     $BODY$     declare     dist_lb real;     begin         dist_lb := (select st_distance(new.geom, st_collect(listed_buildings.geom)) from listed_buildings);              if dist_lb <= 3 then             RAISE WARNING 'Check planning boundary - it is within approx %m of a listed building', round(dist_lb);         end if;         if dist_lb <= 1 then             RAISE EXCEPTION 'Check planning boundary - it is within approx %m of a listed building', round(dist_lb);         end if;         RETURN new;     end;     $BODY$       LANGUAGE plpgsql;          CREATE TRIGGER check_constraint_proximity_trigger       AFTER INSERT OR UPDATE       ON input_table       FOR EACH ROW       EXECUTE PROCEDURE check_constraint_proximity();          insert into listed_buildings (geom) values (st_geomfromewkt('polygon((0 0, 0 1, 1 1, 1 0, 0 0))'));          insert into input_table (geom) values (st_geomfromewkt('polygon((2 2, 2 3, 3 3, 3 2, 2 2))'));     insert into input_table (geom) values (st_geomfromewkt('polygon((0 0, 0 1, 1 1, 1 0, 0 0))'));


Problems using geoserver and postgis to perform wfs transaction?

I have been following the tutorial for transaction. I am using apache webserver and jetty powered geoserver for hosting and using proxy.cgi file for bypassing XHR requests. as described in tutorial I have made necessary changes to my file. I can display the layer, add, edit and delete new features but not able to save them. I have checked the state of new features and it says 'INSERT'(using alert). the line saveStrategy.save() is also called as the alert line immidiately after is working (an XHR transaction returns successful response). but whenever I refresh the page the new features do not show up, and they are not saved in postgis also.               OpenLayers.ProxyHost= "/cgi-bin/proxy.cgi?url=";     //variables     var host = "http://127.0.0.1";     var port = ':8888';     var centerX = 85.33333;//9497800;     var centerY = 27.68888;//3212000;     var range = .5;//10000;     var map_bound = [centerX-range,centerY-range,centerX+range,centerY+range];     var map_wd = 390;     var map_ht = 400;     var legend_wd = 100;     var legend_ht = map_ht;     var tbar_wd = map_wd + legend_wd;     var tbar_ht = 25;     var zoom = 15;          var DeleteFeature = OpenLayers.Class(OpenLayers.Control, {         initialize: function(layer, options) {             OpenLayers.Control.prototype.initialize.apply(this, [options]);             this.layer = layer;             this.handler = new OpenLayers.Handler.Feature(                 this, layer, {click: this.clickFeature}             );         },         clickFeature: function(feature) {             // if feature doesn't have a fid, destroy it             if(feature.fid == undefined) {                 this.layer.destroyFeatures([feature]);             } else {                 feature.state = OpenLayers.State.DELETE;                 this.layer.events.triggerEvent("afterfeaturemodified",                      {feature: feature});                 feature.renderIntent = "select";                 this.layer.drawFeature(feature);             }         },         setMap: function(map) {             this.handler.setMap(map);             OpenLayers.Control.prototype.setMap.apply(this, arguments);         },         CLASS_NAME: "OpenLayers.Control.DeleteFeature"     });     Ext.onReady(function() {                  var saveStrategy = new OpenLayers.Strategy.Save();         /*saveStrategy.events.register("start", '', alert('saving'));         saveStrategy.events.register("success", '', alert('save success'));         saveStrategy.events.register("failure", '', alert('save failure'));         */         //map configuration         var map_options = {             allOverlays:true,             maxExtent:new OpenLayers.Bounds(map_bound[0],map_bound[1],map_bound[2],map_bound[3]),             numZoomLevels: zoom+10,             controls: [new OpenLayers.Control.PanZoomBar(),                 new OpenLayers.Control.MousePosition(),                 new OpenLayers.Control.ScaleLine(),                 //new OpenLayers.Control.Scale()             ],             projection: "EPSG:4326"         }         var map = new OpenLayers.Map(map_options);         var navi = new OpenLayers.Control.NavigationHistory();         map.addControl(navi);              //layers configuration              var roadswgs = new OpenLayers.Layer.WMS(             "roadswgs"             ,host+port+"/geoserver/wms"             ,{  layers: "postgis:trans_ln_8"                 ,transparent: true                 //,styles: "line"             },             {isBaseLayer: false}         );              var portwgs = new OpenLayers.Layer.WMS(             "airportwgs"             ,host+port+"/geoserver/wms"             ,{  layers: "RIS:airport_runway_wgs"                 ,transparent : true                 ,styles: "simple_roads"/*,visibility: false*/},             {isBaseLayer: false}         );              //var vector = new OpenLayers.Layer.Vector('editable');              var vector = new OpenLayers.Layer.Vector('editable'             ,{  strategies: [new OpenLayers.Strategy.BBOX(),saveStrategy]                 ,protocol: new OpenLayers.Protocol.WFS({                     url:  "http://127.0.0.1:8888/geoserver/wfs",                     featureType: "airport_runway",                     featureNS: "http://localhost:8888",                     //version: "1.1.0",                     //geometryName: "geom"                 })                 ,styleMap: new OpenLayers.StyleMap({                     pointRadius: 20,                     strokeColor: "#ff0000",                     strokeWidth: 1,                     fillOpacity: 1                 })             }         );              var wfsroads = new OpenLayers.Layer.Vector("PostgisWFS",              {   strategies: [new OpenLayers.Strategy.BBOX()],                 protocol: new OpenLayers.Protocol.WFS({                     url:  "http://127.0.0.1:8888/geoserver/wfs",                     featurePrefix: "postgis",                     featureType: "trans_ln_8",                     version: "1.1.0"                                  }),                 styleMap: new OpenLayers.StyleMap({                     pointRadius: 3,                     strokeColor: "#ff00ff",                     strokeWidth: 3,                     fillOpacity: 0                 })             }         );              var osm = new OpenLayers.Layer.OSM("OSM");              //var gsat = new OpenLayers.Layer.Google("Google Satellite",{type:G_SATELLITE_MAP, numZoomLevels: 22});              map.addLayers([wfsroads,roadswgs,vector,portwgs/*gsat*/]);              //mappanel configuration         var mapPanel = new GeoExt.MapPanel({             renderTo: 'gxmap',             height: map_ht,             width: map_wd,             map: map,             title: 'A Simple GeoExt Map',             centre: [centerX,centerY],             zoom: zoom-2         });         //layer tree configuration         var layerList = new GeoExt.tree.LayerContainer({             text: 'All Layers',             layerStore: mapPanel.layers,             leaf: true,             expanded: true         });              var layerTree = new Ext.tree.TreePanel({             title: 'Layers',             renderTo: 'layerTree',             root: layerList,             width: tbar_wd,             height: '*'         });              //legend configuration         var legendPanel = new GeoExt.LegendPanel({             title:'Legend',             layerStore: mapPanel.layers,             renderTo: "legend",             border: true,             width: legend_wd,             height: legend_ht         });         //toolbars         toolbars(map,navi,vector,mapPanel); //this part puts simple tools like zoomin, zoomout, maxextent and it is working fine         var nav = new GeoExt.Action({             text: "nav",             control: new OpenLayers.Control.Navigation(),             map: map,             // button options             toggleGroup: editor,             allowDepress: false,             pressed: false,             tooltip: "navigate",             // check item options             //group: "draw",             checked: false         });         //editor toolbar         var editor = "editor tools";              var snaproad = new OpenLayers.Control.Snapping({             layer: vector,             targets: [{                 layer: wfsroads,                 tolerance: 15,                 greedy: true             }]         });         snaproad.activate();         var snapvector = new OpenLayers.Control.Snapping({             layer: vector,             targets: [{                 layer: vector,                 tolerance: 50,                 greedy: true             }]         });         snapvector.activate();         var draw = new OpenLayers.Control.DrawFeature(             vector,     //layer             OpenLayers.Handler.Path, //Handler             {   multi:true,                 type: OpenLayers.Control.TYPE_TOOL,             }//options         );         // DrawFeature control, a "toggle" control         draw_ac = new GeoExt.Action({             text: "draw",             control: draw,             map: map,             // button options             enableToggle: true,             tooltip: "select feature",             toggleGroup: editor         });         /*draw.featureAdded = function(feature) {             feature.layer.eraseFeatures([feature]);             // cast to multilinestring             feature.geometry = new OpenLayers.Geometry.MultiLineString(                 feature.geometry             );             feature.style.strokeColor = "#ff0000";             feature.state = OpenLayers.State.INSERT;             feature.layer.drawFeature(feature);         }         */                   var modify = new OpenLayers.Control.ModifyFeature(vector,{type: OpenLayers.Control.TYPE_TOGGLE});         var select = new OpenLayers.Control.SelectFeature(vector,{hover:true});         map.addControl(select);         select.activate();              var modify_ac = new GeoExt.Action({             text: "modify",             control: modify,             map: map,             // button options             enableToggle: true,             tooltip: "modify feature",             toggleGroup: editor         });                   var del = new DeleteFeature(vector, {title: "Delete Feature"});         var del_ac = new GeoExt.Action({             text: "delete",             control: del,             map: map,             // button options             enableToggle: true,             tooltip: "delete feature",             toggleGroup: editor         });                   //add controls to map         mapPanel.map.addControls([draw,modify,del,select]);         /*         // SelectFeature control, a "toggle" control         select_ac = new GeoExt.Action({             text: "select",             control: new OpenLayers.Control.SelectFeature(vector, {                 type: OpenLayers.Control.TYPE_TOGGLE,                 hover: true             }),             map: map,             // button options             enableToggle: true,             tooltip: "select feature"         });              var draw_but = new Ext.Button({                 //associate tool with button             text: 'Draw Line',             enableToggle: true,             toggleGroup: toggleEditor,             handler: function(toggled){                 if (toggled) {                     snaproad.activate();                     draw.activate();                 } else {                     draw.deactivate();                 }             }         });              var modify_but = new Ext.Button({                   //associate tool with button             text: 'Edit Line',             enableToggle: true,             toggleGroup: editor,             handler: function(toggled){                 if (toggled) {                     draw.deactivate();                     snapvector.activate();                     select.activate();                     modify.activate();                 } else {                     modify.deactivate();                 }             }         });                    var delete_but = new Ext.Button({                   //associate tool with button             text: 'Delete Line',             enableToggle: true,             toggleGroup: editor,             handler: function(toggled){                 if (toggled) {                     draw.deactivate();                     select.activate();                     modify.deactivate();                     del.activate();                 } else {                     del.deactivate();                 }             }         });         */         var save = new GeoExt.Action({             text: "save",             //icon: 'icons/icon_zoomnext.png',             //control: saveStrategy.save(),             onClick: function(){                 saveStrategy.save();                 alert('saved');             },             disabled: false,             //deactivateOnDisable: true,             //activateOnEnable: true,             tooltip: "save road"         });              /*var save_but = new Ext.Button({             text: "Save Changes",             enableToggle: false,             toggleGroup: toggleEditor,             click: function() {                 alert('sve atart');                 map.layer[2].commit();                 alert('not saved');             }         });*/              var toolbar = new Ext.Toolbar({                 //embed tool in toolbar             renderTo: 'tools_editor',             width: tbar_wd,             height: tbar_ht,             autoScroll : false,             autoWidth : true,             autoHeight : true,             buttonAlign: 'left',             fieldLabel : 'Editor tools',             items: [draw_ac,'-',modify_ac,'-',del_ac,'-',save]         });              // the part below is just to make the bounds show up on the page         function updateBounds(){             var bbox = map.getExtent().toArray();             code = "Bounds      lower left: "+bbox[0]+','+bbox[1]+"     upper right"+bbox[2]+','+bbox[3];             document.getElementById('bounds').innerHTML = code;         }         // update the bounds with each map move         map.events.register('moveend', map, updateBounds);         // and update the bounds on first load         updateBounds();          });      I am using openlayers and geoext, and this script is called from a php file. there has been some development when this code                   draw.featureAdded = function(feature) {             feature.layer.eraseFeatures([feature]);             // cast to multilinestring             feature.geometry = new OpenLayers.Geometry.MultiLineString(                 feature.geometry             );             alert(feature.state);             //feature.style.strokeColor = "#ff0000";             feature.state = OpenLayers.State.INSERT;             feature.layer.drawFeature(feature);         }      is enabled, after saving the postgis table shows increase in no.of row but the geometry is empty. the request generated is               <wfs:Transaction xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.1.0" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">     <wfs:Insert><feature:airport_runway xmlns:feature="http://geoserver/postgis">     <feature:geom>     <gml:MultiCurve xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:4326"/>       </feature:geom>     </feature:airport_runway>     </wfs:Insert>     </wfs:Transaction>      and the response is               <wfs:transactionresponse version="1.1.0" xsi:schemalocation="http://www.opengis.net/wfs http://localhost:8888/geoserver/schemas/wfs/1.1.0/wfs.xsd" xmlns:ogc="http://www.opengis.net/ogc" xmlns:tiger="http://www.census.gov" xmlns:cite="http://www.opengeospatial.net/cite" xmlns:nurc="http://www.nurc.nato.int" xmlns:sde="http://geoserver.sf.net" xmlns:postgis="http://geoserver/postgis" xmlns:wfs="http://www.opengis.net/wfs" xmlns:topp="http://www.openplans.org/topp" xmlns:geoserver_tomcat_amrit="http://localhost:8080/geoserver" xmlns:it.geosolutions="http://www.geo-solutions.it" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:sf="http://www.openplans.org/spearfish" xmlns:ows="http://www.opengis.net/ows" xmlns:gml="http://www.opengis.net/gml" xmlns:ris="http://localhost:8888/geoserver" xmlns:xlink="http://www.w3.org/1999/xlink">     <wfs:transactionsummary>     <wfs:totalinserted>     1     </wfs:totalinserted>     <wfs:totalupdated>     0     </wfs:totalupdated>     <wfs:totaldeleted>     0     </wfs:totaldeleted>     </wfs:transactionsummary>     <wfs:transactionresults>     <wfs:insertresults>     <wfs:feature>     <ogc:featureid fid="airport_runway.65"/>     </wfs:feature>     </wfs:insertresults>     </wfs:transactionresults>     </wfs:transactionresponse>      if the above code is commented then transaction returns this error               <ows:exceptionreport version="1.0.0" xsi:schemalocation="http://www.opengis.net/ows http://127.0.0.1:8888/geoserver/schemas/ows/1.0.0/owsExceptionReport.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"   xmlns:ows="http://www.opengis.net/ows">     <ows:exception exceptioncode="InvalidParameterValue">     <ows:exceptiontext>     Error performing insert: Error inserting features     </ows:exceptiontext>     </ows:exception>     </ows:exceptionreport>


Extracting all field names of a featureclass using query filter programmatically

I want to get all field names present in a feature class using queryfilter programmatically.


How can I solve "Shortest and Fastest With time" routing issue?

Please see the image below. (If you can't see the image please search "Bangalore to Palakkad" in Google Maps.) I am searching for a route from Bangalore to Palakkad. Google shows two routes. I think the first one is the shortest path and other is fastest. But look at the distance and **time**. **385** for the shortest (time to travel : 6hr **51 min** ) and **410** for fastest (time to travel : 6hr **7 min** ). How are they calculating this? ![enter image description here](http://i.stack.imgur.com/jLvVs.jpg)


Lidar - Automated building extraction (Improve Process)

I have the following process for automatically creating building outlines/footprints from Las files and was wondering if anyone can help improve it or has suggestions for other processes to use? **EXPERIMENTS:** **A. ArcGIS - Lidardataset tools** : while you can see just the building class points you can not create polygon automatically. There is the minimum bounding geometry but as it looks at the dataset as a whole it requires you to manually choose the points for each building. **B. External Lidar Analysis tool by Overwatch** \- building extraction is a core functionality. http://www.featureanalyst.com/lidar_analyst.htm (No longer active - July 28/2014) It took approx 8 minutes to create a bare earth (essential) and then extract the buildings. It gives the best results by far, even capturing building under the trees (1&4). This is without even looking at a tutorial as it is a very intuitive interface. ![enter image description here](http://i.stack.imgur.com/0N1bJ.jpg) Some issues are it has missed the covered corridors (3) but has correctly left out the open roof (2). **C. LASTOOLS -** used lasmerge to merge 4 relevant tiles and then lasclip to extract zone of AOI only. arcpy.lasboundary("AOI.las","1","buildings","true","true","shp","AOI_buildings.shp","","#","true") ![enter image description here](http://i.stack.imgur.com/JSx9j.jpg) The output was very jagged so I used # arcpy.SimplifyBuilding_cartography("AOI","AOI_buildings_simplified.shp","4 Meters","0 SquareMeters","NO_CHECK") ![enter image description here](http://i.stack.imgur.com/c2mhL.jpg) and then ran arcpy.SimplifyPolygon_cartography("AOI_buildings_simplified","AOI_s.shp","POINT_REMOVE","1 Meters","0 SquareMeters","NO_CHECK","KEEP_COLLAPSED_POINTS") ![enter image description here](http://i.stack.imgur.com/vPWpC.jpg) 1 and 2 show the issue areas. Part of it is due to bad lidar classification. This gives good building lines for long straight edged buildings but if there is overlap by branches or it is a complex building with multiple edges, roofs on multiple levels close to each other etc, manual creation is necersary for any detailed site level work. For general estimation roof shape the process is accurate enough. Process takes less than 10mins total once the process has been finalized. **My Research:** http://stackoverflow.com/questions/5037588/simplified-or- smooth-polygons-that-contain-the-original-detailed-polygon and the process suggested by @gromix and What are Definition, Algorithms and Practical Solutions for Concave Hull?


Good practice for generating contours from a DEM

We have just received a large set of DEMs at work and I would like to generate contours from them. The DEMs have a resolution of 1m and a size of 1kmx1km. Output from gdalinfo:               Driver: AAIGrid/Arc/Info ASCII Grid     Files: 380000_6888000_1k_1m_DEM_ESRI.asc     Size is 1000, 1000     Coordinate System is `'     Origin = (380000.000000000000000,6888000.000000000000000)     Pixel Size = (1.000000000000000,-1.000000000000000)     Corner Coordinates:     Upper Left  (  380000.000, 6888000.000)     Lower Left  (  380000.000, 6887000.000)     Upper Right (  381000.000, 6888000.000)     Lower Right (  381000.000, 6887000.000)     Center      (  380500.000, 6887500.000)     Band 1 Block=1000x1 Type=Float32, ColorInterp=Undefined       NoData Value=-9999      I know I can use gdal_contour to generate the contours (my blog post on the topic) but I'm wondering what some **best practices** for generating contours are. Are there any rules you should follow to get the most out of the contours but not make stuff up or loose too much information? Say I want to generate three sets of contours:   * 250mm   * 1m   * 5m **Is there anything I should do to the DEM before each set?** **Is post smoothing of the lines a good way to go or is smoothing the raster a better option?**


How to determine if feature A and B are located in the same room?

I have created a map using MapWindows like the following: ![dummy_house](http://i.stack.imgur.com/zI6nD.png) A little explanation here, the is a map of an elderly home. The yellow boxes are room and the white rectangles inside the yellow boxes are the beds. I have a few things that need to be done using QGIS (preferably integrated with R so that I can automate the process):   1. I need to check whether the beds are "appropriately" placed in the room. For example, in the "lower-left quadrant" room, the "upper-right quadrant" bed is place outside the boundary of the room, this is not appropriate and I need to fix it.   2. I need to count how many beds are there in one room. Can I do it in GIS/R?   3. I need to know if a particular bed is in the same room as another particular bed. I have about 40 maps and a total of >2000 beds, I can't do the 3 things above manually. I wonder if some mechanical means can help. I am using MapWindows for creating the maps, QGIS for visualization, and R for data processing (checking the dbf files for the naming of the beds, creation of csv files for giving "information" for each bed). I wonder if the 3 tools will be sufficient for the 3 tasks above. Thanks in advance.


albers projection to Lat, Lon

Thanks in advance for helping me with my issue. This GIS stuff is using parts of my brain that have remained idle since high school, and some parts are working faster than others. I have land data stored as Albers projection and I want it converted to Lat, Lon so I can display it on Bing maps. The prj.adf file in the folder containing my data is this:               Projection    ALBERS         Datum         NAD83         Spheroid      GRS80         Units         METERS         Zunits        NO         Xshift        0.0            Yshift        0.0          Parameters           29 30  0.0 /* 1st standard parallel           45 30  0.0 /* 2nd standard parallel          -96  0  0.0 /* central meridian           23  0  0.0 /* latitude of projection's origin           0.0 /* false easting (meters)           0.0 /* false northing (meters)      I have been searching the web for said transformations, and have come close with the article by James Murtha (http://www.delphi3000.com/articles/article_4479.asp), but this apparently maps Lat Lon to his 1000 scale map, which is not the same scale as my data. The units described in my prj file show "Meters" - so do I need to calculate the number of meters between the 1st and 2nd standard parallel and assume that this is my unit size in the y direction? Does anyone have any formulas (or links to information) that I can apply to map between (Lat, Lon) and the Albers (X,Y) projections? Sorry, but I do not have ARCGIS or any other tools to help me so I have to do this in native code without the conversion libraries in these tools. Thanks - Steve Tufty


ArcMap Crashes on my PC

my ArcMap app crashes after opening ArcTool box. Any ideas what might cause this problem and the solution please. Thanks


Interpolation using GRASS plugin in QGIS

I have been using the IDW interpolation tool through GRASS plugin in QGIS and the interpolation output I am getting is for the entire mapset region. My question is: Is there a possibiltity that we can interpolate a certain small region within a given MAPSET region and not the entire MAPSET region ?


arcpy.sa.Viewshed with OFFSETA

This is using ArcGIS10 arcpy with spatial analyst. I need to call Viewshed, but with an observer height set. I have tried using a Z-aware point feature set, but that only throws an error. I am feeding an arcpy.PointGeometry object with a Raster object of a dem for my inputs. e.g. viewshed = arcpy.sa.Viewshed(dem, source_point) I am using dynamically generated points as the source for my observer points for the viewshed. Feature classes have to be built in memory. I have tried several ways to build an in-memory feature class with an OFFSETA field, but so far it has not worked. Anyone have any idea how to use the Viewshed OFFSETA, or another way to set observer height when calling Viewshed inside python using an in memory feature class? * * * Here is the solution I came up with, using Matthew Snape's advice. dem, source_point, offseta, and offsetb are all inputs to function this is contained in:                  # Construct viewshed input        arcpy.CreateFeatureclass_management("in_memory", "viewshedinput", "POINT")        arcpy.AddField_management("in_memory\\viewshedinput", "OFFSETA", "FLOAT")        arcpy.AddField_management("in_memory\\viewshedinput", "OFFSETB", "FLOAT")        rows = arcpy.InsertCursor("in_memory\\viewshedinput")        row = rows.newRow()        row.shape = source_point        row.OFFSETA = offseta        row.OFFSETB = offsetb        rows.insertRow(row)        del rows, row        # Calculate viewshed for source point        viewshed = arcpy.sa.Viewshed(dem, "in_memory\\viewshedinput")


How to get GeoNetwork to search through PostGIS fields via GeoServer WMS

I can't seem to get GeoNetwork to search through my PostGIS fields via GeoServer WMS. According to the GeoNetwork Documentation: > it should loop over datasets served by the service and produce also metadata > for each datasets I have a simple test table with simple features in PostGIS for testing:               | id | date_creation | date_publication | Author | location                |     +----+---------------+------------------+--------+-------------------------+     | 27 |  12/24/2006   |    12/29/2006    |  MDA   | POINT(133.8855 -5.6701) |     | 24 |  02/11/1952   |    05/10/1952    |  NRO   | POINT(22.8855 -44.6701) |     | 12 |  05/12/1996   |    12/24/1956    |  CCP   | POINT(33.8855 -34.6701) |     | 99 |  12/24/2008   |    12/24/1956    |  MDA   | POINT(44.8855 -44.6701) |        * I've added the PostGIS store in GeoServer and published the Layer.    * I've added the "http://localhost:8080/geoserver/wms" service in a harvesing entry. ![enter image description here](http://i.stack.imgur.com/rLUWo.png) When i run the harvester i get 2 entries in GeoNetwork. But the only way i can find the entries is if i search for the GeoServer published "Layer Name". I will not get any results if i search for an author like "MDA" in my table. Am i missing something? Is there a setting or config for exposing the fields i want to share via metadata on GeoServer or something of the sort? Is this even possible?


How do you export features in OpenLayers 2.13?

I'm new to OpenLayers (2.13). I saw in the example section that it is possible to draw features like lines, points and polygons. http://openlayers.org/dev/examples/draw-feature.html Quite interesting. But there is no export feature with it. I mean drawing features on the map is useless if I can't export it in KML or GML or GeoJSON. Is there any example about exporting features in KML? Thanks


Unable to resolve constructor for: 'dijit.layout.TabContainer'

I have started using ArcGIS for the first time, actually my first time using GIS. So I have taken the demo codes and changed them to what I needed. So first I took some code for a tab control on the left of the map, then I wanted to fill in the second tab's content with drawing tools. But there is parse used two times so I get errors. This is the code: _Layout.cshtml - the partial in MVC:               <!DOCTYPE html>                   <meta charset="utf-8" />         <meta name="viewport" content="width=device-width" />         <title>@ViewBag.Title</title>              <meta http-equiv="Content-Type" content="text/html; charset=utf-8">         <meta name="viewport" content="initial-scale=1, maximum-scale=1, user-scalable=no">         <title></title>         <link rel="stylesheet" href="http://js.arcgis.com/3.8/js/dojo/dijit/themes/nihilo/nihilo.css">         <link rel="stylesheet" href="http://js.arcgis.com/3.8/js/dojo/dijit/themes/claro/claro.css">         <link rel="stylesheet" href="http://js.arcgis.com/3.8/js/esri/css/esri.css">                                                                    @RenderSection("Head", false)             @Styles.Render("~/Content/css")             @Scripts.Render("~/bundles/modernizr")                               @RenderBody()                  @Scripts.Render("~/bundles/jquery")             @RenderSection("scripts", required: false)                        here I have `<script>var dojoConfig = { parseOnLoad: true };</script>` \- this was in the example with the tab control. and the Index.cshtml that uses the _Layout partial:                   @{         ViewBag.Title = "Map";         Layout = "~/Views/Shared/_Layout.cshtml";     }

             This is the header section
subheader

                     Draw:
                                          <button data-dojo-type="dijit/form/Button">Point</button>                     <button data-dojo-type="dijit/form/Button">Multi Point</button>                     <button data-dojo-type="dijit/form/Button">Line</button>                     <button data-dojo-type="dijit/form/Button">Polyline</button>                     <button data-dojo-type="dijit/form/Button">Polygon</button>                     <button data-dojo-type="dijit/form/Button">Freehand Polyline</button>                     <button data-dojo-type="dijit/form/Button">Freehand Polygon</button>                                          <button data-dojo-type="dijit/form/Button">Arrow</button>                     <button data-dojo-type="dijit/form/Button">Triangle</button>                     <button data-dojo-type="dijit/form/Button">Circle</button>                     <button data-dojo-type="dijit/form/Button">Ellipse</button>

             this is the footer section

               @section Head     {         ` and the tabs don't appear, and the drawing doesn't work. But sometimes, after a few refreshes, the tabs appear and the drawing works. But can't find the cause.


What version of SVG should I be saving my images with for QGIS? Any particular options I should use/avoid?

I have a illustration in Adobe Illustrator that I'm trying to export as an SVG file and add to my map layout using the Print Composer of QGIS 1.7.3. When I add the image, I get strange artifacts that aren't in the Illustrator version. I tried a few different version of SVG in the Save As dialog, but they all end up generating the strange artifact. What version of SVG should I be saving my images with? Any particular options I should use/avoid? The orange block/white dots pattern below is the weird artifact: ![svg artifact in qgis](http://i.stack.imgur.com/Qu79S.png) Here are the Save As options for SVG in Illustrator CS 3: ![enter image description here](http://i.stack.imgur.com/pOUAD.png)


'invalid url' on ArcGIS server REST API

From time to time, the services on our ArcGIS server will not be available, and when requesting tiles, I will get a 'invalid url'. This is okay, sometimes it is scheduled downtime, maintenance or just some of our server scripts not bringing the services up again. We also have a cache server in front of the ArcGIS server, that can serve often requested tiles really fast, w/o sending the request to the server. The cache server is set up to store the successful requests for a period of time. The thing that bothers me, is that the page 'invalid url' is served with the http code `OK 200`, and thus making our cache server store it. Of course, I could make exceptions, probably based on mime-type or something, but that's not the optimal solution. ![Example of the invalid url page](http://i.imgur.com/HLCmT.png) I would actually like to know if it is possible to alter the 'invalid url' error page, with another http code, such as `Bad request 400` or `Not implemented 501`, both which would make Varnish behave. EDIT: I should mention that the ArcGIS server runs in a .NET environment.


Vector Machine Classification in ArcGIS

Is it possible to perform a vector machine classification of raster image using the ArcGIS software package. This supervised classification method is available in the ENVI software but not for ArcGIS. I have integrated the ENVI toolbox into my ArcCatalog Toolbox set. While the ENVI toolbox offers me the tools to perform supervised classification using four different methods, it does not include the Vector Machine Classification Method. Since, I've access to an ArcGIS license but not an ENVI license, I am looking to access this functionality in ArcGIS. Any suggestions on how this can be achieved?


Point markers with icon/picture symbols showing white area around them on ArcMap?

I'm trying to use some icons (pictures) but i haven't been successful yet. This is where i got the icons: http://mapicons.nicolasmollet.com/category/markers/offices/ I want to build a 1 meter by 2 meter tourism map of a specific area of Lisbon using these icons to ilustrate where are the restaurants/bars/banks/supermarkets/etc... ![icons](http://i.stack.imgur.com/3X6P6.jpg) What happens when i add an icon in the arcmap is this:icones2.jpg They have this white area around them but i would really like to take it away. I tried to get some help from some photoshop geeks but whatever the solution, when i upload them to arcmap they just don't change and still have the white edge.


Do SAS Bridge to Esri and SAS/GIS extensions add value to Spatial Analyst and Spatial Statistics?

My organization uses SAS extensively, but has been increasing GIS use and capacity recently (that's why they hired me). I was asked to evaluate the two SAS GIS extensions _Bridge to Esri_ and _SAS/GIS_. I've found some documentation and a lot of mostly sparse and vague marketing materials... see... http://support.sas.com/rnd/datavisualization/BridgeForESRI/V2/ http://support.sas.com/rnd/datavisualization/papers/SASMapping.pdf http://support.sas.com/documentation/cdl/en/apdatgis/62945/PDF/default/apdatgis.pdf I was wondering if anyone here had any experience, insight and input with these extensions? I gather that they provide the ability to run some spatial stats processes in SAS environments, but if we have a spatial analyst license and access to the ESRI spatial statistics tools, is there any real value added here? I've found exporting and importing DBF and CSV files between SAS and ArcGIS fairly fluid.


Change Case in .dbf file QGIS

Asked another question below and this will develop on from this. The Bedrock 625K BGS mapping is all in UPPERCASE. It looks clumsey and wrong in the legend. I wish to change it to Sentence Case. I know the long winded way is to do each entry manually...there is over 7000 of them...so no. Working in either QGIS 1.8.4 or 2.0


Excel Distance Calculation

Is there a way (preferably in Excel) to calculate the distance (based on lat, lon) between two lists of points? My end goal would be, for example, to have a list of Starbucks and a list of McDonald's, and then show the nearest neighbor, if you will. Thanks in advance.


Is it possible to use JTS for calculations on data in geographic coordinates?

I have some spatial data that has all of its coordinates as lat/lon pairs (with about 10 digits decimal precision), it's stored in a database as WGS84 data. Some of the data is represented as polygons which are the resulting union of some smaller polygons whose boundaries are stored. I also have a number of points from which I build line segments (just 2 points in each segment) which I use later for intersection tests with the polygons. I'm using a SpatialIndex to improve my queries so I insert the envelopes of all polygons in a tree (tested with both QuadTree and STRtree). Then, I connect two points into a line segment and I'm using its envelope to query the tree for possible intersections. The problem is that I get pretty much all the polygons as a result which is clearly wrong. To give you some idea about the real scale of my data, I have about 100 polygons that cover the whole North America, each line covers a very very small part of a single polygon. Ideally, I would expect no more than 2 polygons as a result. I'm using JTS to do this calculation and I'm aware that it's not really suited for spherical data so can you suggest me another library/tool to achieve the desired behaviour or possible a workaround (for example, projecting before using JTS)? EDIT: I'm using postgresql + postgis extension, storing the all spatial data as geometry in 4326.


How to connect places (points) without crossing lines?

I have two sets of data. One with placemarks represented by blue dots in image below, and another set of geographic data with lines represented by grey lines. 1) For the first point in drawing, I'm looking for a way for automatically connecting placemarks, by adding new lines geometry between them, but not crossing grey lines (representing streets in fact). In drawing, green connections are connections wanted, and red connections are unwanted ones. 2) This second point is maybe more difficult to achieve. Connections between two dots should not exist if they're not "following" (maybe by having roughly the same angle of) the grey lines. Is it possible to play with data to solve this problem by playing with QGIS, or another tool, and save new lines in database? ![Schema](http://i.stack.imgur.com/89qDw.png) P.S.: I'm new to GIS world.


How to bulk change annotation locations?

I'm working with a large number of labels that I've converted to annotation (using ArcMap 10.1). All of them need to fit within the parcel borders that they correspond to. Is there a way to do it all at once or am I stuck going one by one to fit them in the boundaries?


Difference between ArcObjects SDK for Java and ArcGIS Runtime for Java?

I would like to develop spatial tools and automate my spatial analysis process (including converting XY into shp, proximity analysis, network analysis, geocoding etc on thousands of files) in Java. I was wondering if ArcGIS runtime sdk for Java can do the task or do I need ArcObjects sdk for Java. I can see that they are two different things but I am really confused what is the difference between them and which one does what? I am also looking for some tutorial.


How to add a progress bar to a python plugin?

I would like to add some type of progress bar to my Qgis python plugin. Can someone point me to an example?


Is there a way to scale to 1:1 in QGIS Cloud?

Someone knows if with QGIS cloud Pro Plan, it's possible have a scale till 1:1. I use QGIS to draw small stuff on my gis (archaeological finds, as really small pottery's sherds) and I need a scale smaller than 1:624 as now I can view with QGIS cloud free plan. Thank a lot and best regards


Latitude and Longitudes are replaced with zeros when importing excel 2007 worksheet into personal geodatabase

I am trying to import an Excel 2007 data set into a personal geodatabase so I can convert coordinate systems. The data set includes point locations in Decimal Degrees. When I import the Excel worksheet, the latitude and longitudes in the resulting geodatabase table have been replaced with zeros. Why is this? This didn't occur with a different data set I have. Thanks Jason


How can I georeference a second image without losing georeferencing for first image in ArcGIS 10?

What happens: 1) Take a map with an image that was georeferenced earlier, which works fine. 2) Add a second image to the TOC and fire up the georeferencing toolbar. 3) Select the new, ungeoreferenced image from the drop-down list... 4) ... georeferencing on first image vanishes. How do I georef the second image without losing the first georeferencing?


Near tool producing odd results for batch-snapping; how to best perform batch-snapping?

I am trying to use the _Near_ tool in Arc10 as described in Snapping points to a line using ArcGIS geoprocessing tools for the purposes of batch-snapping a feature class of points to a feature class of lines. The XY coordinates produced by the _Near_ tool do not exactly map onto the lines - there can often be about a foot distance between the line and the new point. The new points also do not appear to be on the shortest-distance path from the original points. The projections are the same for the two feature classes and the data frame, but points are in meters while lines and data frame are in feet. In attempting to fix this problem, I have projected my original points into feet before using _Near_ , with no change in the outcome. Any suggestions on why this method doesn't work for me, and how to get my points to batch-snap to lines, are greatly appreciated!


Reproject TM65 to OSGB in ArcGIS

I have a TIFF file in the TM65 (Transverse Mercator) projection, and want to convert it to the OSGB (British National Grid) projection. I'm trying to do this in ArcGIS using the Project Raster tool. I can select the input and output co-ordinate system fine, but it then asks me to choose a Geographic Transformation. It seems that I have to choose two - one to go from TM65 to WGS84, and then another to go from WGS84 to OSGB. My question is - which one should I choose? What are the differences between them? I have tried choosing some almost at random (and also tried matching them up - so I do TM65->WGS84_2 and then WGS84_2 to OSGB) and each time once the points are reprojected they are displaced from their true locations - the amount varies between fairly small (but large enough to be an issue) and quite large. Is there any way I can get an accurate conversion between TM65 and OSGB? Am I using the wrong tool?


cannot get EPSG of bounding box

I use openlayers 2.12. I created a form where the user can insert a bounding box, to insert to the database a new borough. To check the `EPSG` of the borough I tried the following               //ba bb bc bd came from the form     var thisIsNew=new OpenLayers.Bounds(ba, bb,bc, bd).getProjection();     alert (thisIsNew);      and also               var thisIsNew=new OpenLayers.Bounds(ba, bb,bc, bd);     var newest = thisIsNew.getProjection();     alert(newest);      I never get an alert and on the console I see that "the object does not support that method", for both of the above ways. What am I missing? Thanks in advance


Is the postgis layer option in qgis 2.0.1 Dufour two way? Was the spit plugin built in?

I am trying to install postgresql with postgis on a remotely accessed terminal only installation of ubuntu 12.04 that I want to use to connect to qgis 2.0.1 and hopefully have a two way transfer between editing data and committing it to my database and grabbing it to edit. If so, is there any known compatibility issues if I use postgresql 9.2 and postgis 2.0 given that qgis 2.0.1 is only a month old? Essentially I want to be able to edit data in qgis and have a postgis database connected to geoserver for display.


How to add a geometry to an existing feature that was created with NULL geometry?

I have a simple SpatiaLite geodb loaded in QGIS with one table. There is a POINT geometry field and some alphanumeric data. In some cases I can record data about a specific site, but the coordinates are not available at this time and will be retrieved and recorded later. In some other cases I have the coordinates and I want to insert them right into the geometry field, preferably by clicking on the map canvas in the same way available for adding new features (having a separate text field for the coordinates and generating the geometry from that is a no-go). How do I add a geometry to an existing feature that was created with a NULL geometry field?


Error while placing labels within polygons using ArcGIS JS API and a external Geometry Service

I have published in my ArcGIS server a map service with a table containing the polygons of all countries on Earth. So this table contains more than 200 rows and the geometry field is of type `type: esriFieldTypeGeometry , alias: Shape`. What I would like to do is creating a label inside this polygon, so what I did was modifying this ESRI example (Create points for labeling) that, instead of allowing the user drawing a polygon to be sent to the **external Geometry Service** , expects a click from user's side, retrieves clicked-country geometry with a Query operations and draws a predefined label on the country. Of course, last step is not done yet, but this would be the idea. The problem is that after building a new `Graphic` object with the new geometry retrieved from the layer and call ESRI Geometry Service as in the example, I got an error on Firebug like this:               TypeError: _c94 is undefined     ...eateSegments:function(_764){_764.shape.path=_764.vmlPath;_764.segmented=false;_7...      And well, I don't know why is this happening as **error is not verbose** at all. I noticed that the error is caused by:               http://serverapi.arcgisonline.com/jsapi/arcgis/?v=3.2 (On Line 34)      So these are my doubts:   * The ESRI example comes, out of the box, with the 1.6 version of ArcGIS JS API. And I changed it to 3.2, because if not, I got this error just when starting the application **because I imported the Feature Layer with dojo.require** :  "NetworkError: 404 Not Found - http://serverapi.arcgisonline.com/jsapi/arcgis/1.6/js/esri/layers/FeatureLayer.xd.js" uncaught exception: Could not load cross-domain resources: esri.layers.FeatureLayer   * Is it possible that I need a proxy to pass the country geometry to Geometry Server? I mean, if you check the following image, you can see that there are about 2800 rings. **I am not sure if these data are sent to the geometry service** , but if so, I guess that the **GET request will be huge** and therefore, a proxy should be involved, isn't it? ![enter image description here](http://i.stack.imgur.com/eq1bb.png) Please, can anyone provide some help? Thank you very much!


What changes will be made to arcpy in ArcGIS 10.1?

I heard that two big changes to arcpy will be add-ins and new, faster cursors. Also, new features include "the ability to administer your enterprsie geodatabase (ArcSDE) using Python scripting" and "Dynamic legend creation." I was wondering if anyone knew of other specific changes, maybe something overheard at the 2011 user conference. A changelog of all the changes would be ideal, but I haven't had luck finding anything.


Map and onKeyDown event

I have a web page with a map and search field. Something like Derek's "beer's finder" When I type a search request I somethimes use arrows to fix something. But if mouse is over the map, then I can't use arrows for text editing because map catches them. I tried to fix it the following way (without result):                   dojo.connect(this._currentmap, "onKeyDown", function(ev) {             if (ev.keyCode >= 37 && ev.keyCode <= 40) { //if arrrow                 dojo.stopEvent(ev);             }         });      Do you know how I can fix this problem?


Only display the latest feature added

I'm trying to modify the following sample to only show the latest feature added to the map. The sample is located here: https://developers.arcgis.com/javascript/jssamples/ed_attributeInspectorValidate.html I've tried to change               //Add the editable feature layer to the map             var pointsOfInterest = new esri.layers.FeatureLayer("http://sampleserver5.arcgisonline.com/ArcGIS/rest/services/Notes/FeatureServer/0",{               mode: esri.layers.FeatureLayer.MODE_ONDEMAND,                outFields: ['*']             });      from MODE_ONDEMAND to MODE_SELECTION               //Add the editable feature layer to the map             var pointsOfInterest = new esri.layers.FeatureLayer("http://sampleserver5.arcgisonline.com/ArcGIS/rest/services/Notes/FeatureServer/0",{               mode: esri.layers.FeatureLayer.MODE_SELECTION,                outFields: ['*']             });      Which works when adding a new feature but as soon as the user clicks anywhere else on the map the newly added feature disappears. Is there anyway to modify the code so that the newly added feature remains displayed until another feature is added? Effectively making it so only the latest feature is shown on the map. Thanks!


How to produce density map in square kilometers?

I have a shp file with the location (lat, lon) of lightning spatial distribution. I want to produce a density map of lightning occurrence per square KM. I used density point function and I don't know how to transform the output in square KM. My map's unit is kilometer I should use 1 Km radius? I'm using ArcGis 9,3 I appreciate any answer.


ArcGIS 10.1 Network Analyst: Snapping stops from point-shapefile to route analysis layer doesn't work properly

I'm struggling with a task, which should be not that difficult....but anyway I do not find a solution. **Background** I'm map-matching tracks, generated by a GPS-device, to a road network. Therefore I make use of the ArcGIS 10.1 Network Analyst. I connected the GPS-Point measurements to polylines (1 for each track), created a buffer around them and made an intersection of the buffer-boundary with the street- network to determine point-barriers. According to the paper of _Ron Dalumpines and Darren M. Scott form McMaster University, Ontario, Canada_ I'd like to find the shortest path within the buffered area (confined by the generated point barriers). **Problem** The problem I'm struggling with is the following: I'd like to define the stops (start- and end-point of the track) through a shapefile which contains the start- and endpoint for a certain track. The position of these points already coincidences with network-junctions points. Nevertheless ArcGIS seems not to be able, to link these points with the network and these "stop- points" get marked as "Unlocated". If I use the default search tolerance of 5000 meters, ArcGIS is able to locate the stop-points....but the found ones are not at all the nearest ones to the position of the loaded shapefile. Has anybody an idea what could cause this strange way of snapping? Has anybody already had a similar problem? As a newbie I'm not allowed to upload pictures...otherwise I would have uploaded two pictures explaining my situation. Thanks already in advance for your hints and advices ;) Best Regards, Dom ETH Zurich, Switzerland


Installing ArcSDE engine plugin for Eclipse?

I'm looking for the following `jar`s for the ArcSDE API:   * `jsde_sdk.jar`   * `jpe_sdk.jar`   * `concurrent.jar`   * `icu4_3_2.jar` thanks to this answer I've found that they can be gotten installing the ArcSDE `engine` plugin for Eclipse. But on ESRI's help page it calls for the following directory:               %ARCGISHOME%/java/tools/eclipse_plugin/arcgis_update_site/engine      My installation, instead, only has the following one:               %ARCGISHOME%/java/tools/eclipse_plugin/arcgis_update_site/arcobjects      That is, the `arc Objects plugin`, instead of the `engine plugin`. Where (or how) can I get the ArcSDE `engine` plugin for Eclipse?


Why do I get a gdal_translate error when creating ecw from tif

I have python 2.7 and the latest gdal 32bit core with also GDAL-1.10.1.win32-py2.7.msi and the ecw plugin on windows xp. I have checked that the ecw plugin works by checking the installed formats with gdalinfo --formats. When I try a simple tif to ecw translation I get this error:               Input file size is 157, 157     0ERROR 6: GDALDriver::Create() ... no create method implemented for this      format      Can anyone tell me why?


How do I group multiple polygons into the same record?

I am trying to figure out how to group multiple polygons into the same record (using something like a spatial join tool) using CartoDB. I'd prefer to keep the polygons as separate rather than joining them to create one larger polygon, but I need them to be part of the same record


Using r.in.gdal, how to specify band to extract?

I want to use r.in.gdal to extract a single channel from a multi-channel image. The man page indicates that I can specify which band to extract using the argument "band=" followed by a number in the range (1,4). What is the mapping between (1,4) and (red,green,blue,alpha)?


I want to add a circle which stays the same size

As stated in my title, I want to draw a circle with a fixed diameter. It should not change when I am zooming in or out. How do I do that? (I am using the CartoDB Service). I am trying to visualize how many points on the map are within a 10 mile range for example.


Calculating the closest distance between two points in ArcMap 10.2

I am looking for a script to calculate the closest meter to a set of addresses. I have a feature class of meters and a feature class of addresses. I have populated the address with fields for closest meter and another field for the distance of the closest meter to an address. Do anyone have a python script that would calculate and populate the field in the address feature class with the closest meter to an address and then the distance? I am still learning Python scripting.


Update Point ID as Start Node and End Node to Line through FME

I need your help to perform the below process in FME. I have two layers: Pole Data (Point) and Cable (Line). Cable has no ID. I want to keep the Cable which is going through the Pole. Then make the cable as Single segment between two poles. Then I want to update the Start ID and End ID (Pole ID) to each segment (Cable layer). Please share the design of workbench and what parameters I will use in each transformer.


patch analysis and spatial statistics for Grid data

I have TerraSAR data (raster data) for settlement, I must show the large density building, medium and low density in the map. Is this possible in the spatial analysis or patch analysis? I look forward to hear of you? Regards Neda


WFS request with Openlayers; blank request

I'm having trouble to add a WFS layer over an OSM map. At first I thought that my data had a problem, and I thus decided to adapt one of the example on the Openlayers website. I have no error returned, but the layer is not appearing... Instead, I see in Firebug that there is a WFS request made to the server, but it is empty               http://demo.opengeo.org/geoserver/wfs      No GET request or arguments whatsoever. I can correctly make the request and I get a xml file as a result (in eg Firefox). I thought it was maybe because of the a proxy and I set one up (the one in Openlayers FAQ), but I get the same thing. Here is the code I use:               <!DOCTYPE html>                   <meta charset="utf-8">


How to install shapely on Ubuntu 12.10

I have tried to install shapely on Ubuntu 12.10. Therefore I installed libgeos-c1, pip and shapely with `$ pip install Shapely`. The shapely files are now in /home/user/build/Shapely. When I try to import it in my python script with `from shapely.geometry import Point` gives the following error: `from shapely.geometry import * File "/usr/lib/python2.7/dist- packages/qgis/utils.py", line 309, in _import mod = _builtin_import(name, globals, locals, fromlist, level) ImportError: No module named shapely.geometry` Any idea why it is not working?


How do I choose parameters for the Huff model that best fit my data?

I'm using the Huff model for some analysis I'm doing. For those not aware, the Huff model is a probabilistic model for retail trade area analysis. ![Huff model algorithm](http://i.stack.imgur.com/acXmh.png) `H_ij` is the probability that a person at location `i` will visit a place at location `j`   `A_j` is the attractiveness of location `j`   `D_ij` is the distance the person at location `i` is from location `j`   `\alpha` is the attractiveness enhancement parameter   `\beta` is the distance decay parameter I have a dataset comprised of over two million elements. For each element, I know the person's ZIP code as well as the ZIP code they visited. So, "locations" in this context of the Huff model are ZIP codes. I'm using the number of people that visited a ZIP code as the attractiveness of the location (i.e., `A_j`). What I need is a method for estimating the two parameters, `\alpha` and `\beta`, from data. How can I determine what these parameters should be? [edit]   After thinking about this, I realized that this is really an optimization problem. I think the approach I'm going to take is using a metaheuristic global optimization algorithm:   1. Search will select the `\alpha` and `\beta` parameters randomly, according to its rules.   2. For each entry in the dataset, I select a random number, `r`. I compute the Huff probability, `H`, that the person in the entry visited the ZIP code where I know they visited. If `r < H`, then we got it right.   3. So, the "fitness" for a particular set of parameters, `\alpha` and `\beta`, is the number of cases we get right. This is therefore a maximization problem.   4. We run this lots of times (because each search is stochastic) with reasonable ranges for `\alpha` and `\beta` (something like `0 < \alpha, \beta <= 20`), and the winning `\alpha` and `\beta` are the ones that produce the maximal fitness. Does this seem like a reasonable solution?


GIS Education Guide

I am a Computer Science graduate. I have one year of experience in JAVA programming & I recently started working with GIS applications. I found it so interesting that I would like to do my Post Graduate & PHD in a field related to GIS. Can anybody suggest some online courses, Professional Schools or Universities with a good reputation/curriculum?


FlexViewer 2.1 - Add scroll for widgetgroup

I have a site that is using the ArcGIS Viewer for Flex version 2.1. In the header I have a widget group that contains 25 items. The list is so long that it goes beyond the bottom of the browser window and I can't see the last 5-6 items. Is it possible to add scroll capability to the widgetgroup so I can scroll through the items? Thanks, Tim


Is possible to use ArcView select and query data from Geoserver WFS?

I am looking at implementing a GIS solution using GeoServer. If I implement a GeoServer solution is possible using ArcView via WFS to select and query data? From the research I have done it looks like there is a client plugin cost of around $1500 (Data Interoperability extension). As most of our clients use ArcView it is not practical to ask each of them to purchase this extension. Thanks folks!


Create a geography/geometry column from x and y fields - SQL Server 2008

Hy guys, I'm new in SQL Server 2008 and I hope you will understand my question/need. Thus, I have a table which contains 3 fields (Name, Lat and Long) in my data base (spatial). I want to create a geometry/geography column based on those fields (Lat and Long) but unfortunately without any success. My question is how can I managed to do that? Thank you in advance. Tudor


What tool(s) in the Sextante toolbox in qgis 1.8 can be used for interpolation? Kriging and spline

I have looked at every single one and can not figure out which on can be used.


How to change the color of a tile layer?

If I generate a very large number of map tiles (300,000), is there an easy way to let a Google Maps user change their color? For instance, is it possible to change the color scheme of a map tile layer by overlaying it with other images (red, blue, green) with varying degrees of opacity? My map tiles are showing household income at the census tract level. So that is a single numeric value represented by a color. So I don't risk losing information if I add (or remove) red, blue or green from it.


How do I use the Terrain with Labels basemap in the Map constructor

Seems like all the other basemaps can be used in the esri.Map() constructor, but this one doesn't have a corresponding entry in the list of acceptable values. I've tried just "terrain", and it didn't work.


Setting the number of bars in a histogramm

I am using qgis' 2 "Statist" to build a histogramm for my data. I see that I can set the band width of the bars, but it seems that "Statist" does not provide the option of setting the number of bars. Any ideas on how to set the number of bars in a histogramm would be much appreciated!


Exporting ArcMap produced layer file histogram?

I could simply take a screenshot, or reproduce it in other software, but what I see when I go to layer properties>symbology>classify is perfect for what I need. The only thing is it would be nice to have picture with better resolution, so I'm wondering if I could somehow export this graph. I have ArcMap Basic.


How to write GML and GeoJSON features while retaining attributes in .NET?

I am developing a .NET project [C#] that needs to export/write features as GML and GeoJSON, including, of course, the attributes; data could queried by a database or manually entered. I checked NetTopologySuite and GeoApi, I have also tried DotSpatial and SharpMap, but none seems to be able to include attribute fields and values in the output. Am I wrong? Is there anybody who can help me or suggest the right place to start from?


Mapserver mapscript map inset

What is the recommended way to produce a map similar to the one in http://www.nexusenergy.com.au/6468050/nexus-energy-oil-gas-exploration-and- productio.htm I would like to dynamically create the insets, based on user supplied input. Will be using geodjango to drive mapserver, or is there an alternative ie using qgis (and using geodjango to drive qgis somehow?) Thanks


Unable to load GRASS plugin

I upgraded from QGIS 1.7.0 to 1.8.0, including updating GRASS to version 6.4.3. Since then, GRASS isn't any longer listed in the Fetch Python plugins. I uninstalled and re-installed the whole package without any results. Any tips how to reactivate GRASS in my system (OSX 10.6.8)? Thanks EDIT: Apparently Sextante has to be activated to be able to load the GRASS plugin. However, the "importlib" Python module is lacking to activate Sextante. Any suggestions on how to install this module?


How to rotate u and v components of wind in rotated gribs?

I have a grib message of a rotated hirlam product. The gribmessage states, that the u and v component contained are relative to the rotated grid. How can I unrotate them?


Where do I put pause and balloon visibility commands in KMZ tour?

I am teaching a new course called Global Studies, which looks at many aspects of various regions of the world. As in introduction to each region, I plan on making a Google Earth tour with information, pictures, video, etc... That part I know how to do, but I want to "polish" it. I want to open the balloon once the student arrives at a location and I want the tour to pause. I want the balloon to close when the tour is contninued. From research, I learned that I need to edit the kml file (located inside the kmz file - accessed via zip). I need to add the appropriate code:               (Turning on Balloon)            <gx:AnimatedUpdate>               <Update>                 <targetHref/>                 <Change>                   <Placemark targetId="BLAHBLAHBLAH">                     <gx:balloonVisibility>1</gx:balloonVisibility>                   </Placemark>                 </Change>               </Update>             </gx:AnimatedUpdate>          (Pausing)          <gx:TourControl>         <gx:playMode>pause</gx:playMode>     </gx:TourControl>          (Turning off Balloon)            <gx:AnimatedUpdate>               <Update>                 <targetHref/>                 <Change>                   <Placemark targetId="BLAHBLAHBLAH">                     <gx:balloonVisibility>0</gx:balloonVisibility>                   </Placemark>                 </Change>               </Update>             </gx:AnimatedUpdate>      The issue I have is that I can't seem to figure out WHERE to place it within the doc.kml. The examples show simple "flyto" commands for each location. The doc.kml in the kmz file has way more flyto's than my 16 locations. If I search for the lat/long, there are huge numbers as well. If I search for the location/placemark names, they are only found in the beginning of the file, where there are no flyto statements. Is there a way I can figure out where to place the code to achieve the polished look I am trying to accomplish? My kmz file can be found here. I want to learn how to do this, so in the future I can do it with other tours and other regions.


Is there a database/service that allows me to get the terrain type to a coordinate?

I hope i'm not totally wrong with this kind of question at this stackexchange page. Is there any sort of service avaiable that gives me (or allows me somehow to generate) the areatype to a coordinate? Something like (lat,long) -> position is at a beach. I need this kind of data on a more or less global scale, but "only" at those places, where you have at least a gsm (not satallite) connection - so "only" those areas that are more or less close to populated places. The closest thing i've found would be to analyse map pictures and store the results, which doesn't seem very elegant to me and probably would take alot of time. Is their any better option avaiable?


How to dissolve inner polygons in a multipolygon

I'm trying to get the boundary of a collection of zip codes. Here is my current sql:               SELECT ST_AsText(ST_SimplifyPreserveTopology((ST_Union(geom)), .003)) as poly     FROM zcta5_data     WHERE zcta5 IN ('zipcodes')      Here, is what my result looks like. How do I get rid of the inner polygons? ![enter image description here](http://i.stack.imgur.com/AFgtq.png)


Determine density or cover from a raster layer that has been clipped by a polygon

I am a very new GIS user using ArcMap 10. I have a raster layer of tree canopy that I have clipped so that the layer is in a polygon that is a circle with a 2km radius. I have done this for four circles and I am trying to compare the tree canopy across the four circles. Something like a density or % area covered would be perfect. The problem is... the tree canopy raster layer attribute table only has counts for where the canopy is... and no counts for where it is not. I could get a % cover (for comparisons across the four circles) if I knew what the count of the no tree canopy cells was, but this information is absent. It is almost as if the raster layer is like a point layer, but I cannot count each cell as a point. I hope this makes sense. I am completely lost.. any help would be appreciated! Thanks in advance!


Why does KMZ to GPX reverse the objectid?

I am working on KMZ to GPX conversion and every time I use the data interoperability extension and quick export it seems to reverse the objectid, now my 1st waypoint is the last placemark point of the KMZ file that was used to convert it. I am using ArcGIS desktop advanced license 10.2.2. Does anyone know why this would be happening? Very frustrating!


How to manually reclassify a contiguous region of misclassified pixels?

The attached image is the result of an image segmentation algorithm on high resolution imagery. The smaller clusters of pixels represent trees and brush and the larger cluster of pixels (marked with a red arrow) is a misclassified farm field. In ArcGIS 10.1 is there a way to select the misclassified pixels and reclassify them to 0 or NoData? ![enter image description here](http://i.stack.imgur.com/lamUZ.png)


Contour generation using Interpolation - Correcting edge Effects

I utilised the Interpolation feature in QGIS to prepare contours from a DEM I had generated. Sometimes when contours intersect the edge of the DEM area, they run along the edge for some distance. This doesn't seem to occur all the time, but only at select locations. Is there a method to avoid this? It seems if I cut 1m of the edge of the contour they would be perfect.


Batch Converting ASCII into Feature Classes using shapefile lists

Problem. I have 413 ASCII files totalling 4 billion points, and my system cannot handle converting them to MULTIPOINT so that I can generate Terrains. Therefore I would like to create 12+ subsets of the area, using the lists generated from a shapefile(s).   1. I've made my seperate Feature Classes listing the variable number of files that will make up each subset.   2. I can generate (and print) a list of files from the Feature Class using the SearchCursor function.   3. I'm running into problems using the generated list(s) to actually create the MULTIPOINT Feature Classes. Here is my code:               # Import modules     import arcpy, os, sys, traceback     from arcpy import env               # generate a list of ascii files from the area-subset feature classes     env.workspace=r"F:\RTAMS\Petawawa\11-Library\LEGEO_Aug2012\LEGEO_Aug2012.gdb\LiDAR_ASCII"     fc=arcpy.ListFeatureClasses()          for fc in fc:         cur=arcpy.SearchCursor (fc)              for row in cur:             fileName=row.TILENAME                  #this is the location of the ASCII files we want to import, so we change the workspace             env.workspace=r"F:\RTAMS\Petawawa\11-Library\LEGEO_Aug2012\LiDAR\UTM18_WGS84\ASCII\FEATURES"                  # we should now have a list of files, which we will use to import into area-subset MULTIPOINT files             # we set up some parameters for the ASCII3D to Feature Class tool                  try:                 arcpy.CheckOutExtension("3D")                 outWorkspace=r"F:\RTAMS\Petawawa\11-Library\LEGEO_Aug2012\LEGEO_Aug2012.gdb\LiDAR_ASCII_to_MULTIPOINT"                 inFormat="GENERATE"                 outFC=os.path.join(outWorkspace,fc+"_a")                 outType="MULTIPOINT"                 # find the ascii files listed in the search cursor                 txtList=arcpy.ListFiles(fileName)                 #Execute the ASCII3D_to_Feature_Class                 arcpy.ASCII3DToFeatureClass_3d(txtList,inFormat,outFC,outType)                 print outFC + " is complete."                  except arcpy.ExecuteError:                 print arcpy.GetMessages()                  except:                 # Get the traceback object                 tb = sys.exc_info()[2]                 tbinfo = traceback.format_tb(tb)[0]                 # Concatenate error information into message string                 pymsg = 'PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}'\                     .format(tbinfo, str(sys.exc_info()[1]))                 msgs = 'ArcPy ERRORS:\n {0}\n'.format(arcpy.GetMessages(2))                 # Return python error messages for script tool or Python Window                 arcpy.AddError(pymsg)                 arcpy.AddError(msgs)      Hoping you can help me out.


Running ArcGIS Data Reviewer using ArcPy?

I want to run the ArcGIS Data reviewer in Python code, to check domains. My scenario is check the attribute values if its meet the domain value. Is it possible to run ArcGIS Data Reviewer in Python code ?


pgrouting: Directed Shortest Path calculation using Dijkstra ignores one-way streets

I think I found an interesting bug in pgrouting. When I use the Dijkstra algorithm (no matter if I use the shortest_path() or dijkstra_sp_delta_directed() method) to calculate the directed shortest paths between two particular vertexes (id 72846 and id 85556), it only delivers the correct result for one direction. I experimented a little bit and then I found out that there are two arcs in the roads_po table with opposite direction which connect two vertexes that are part of both calculated shortest paths (gid 700 and gid 169244). When I calculate the shortest path from vertex 72846 to vertex 85556, everything works fine as can be seen here: correct_result: ![enter image description here](http://i.stack.imgur.com/1blzV.png) But when I calculate the shortest path from vertex 85556 to vertex 72846 (the opposite direction), the wrong arc with gid 169244 instead of the correct arc with gid 700 is chosen resulting in a violation of the one-way direction. wrong_result: ![enter image description here](http://i.stack.imgur.com/maQlQ.png) It seems that pgrouting calculates both paths correct, but it has a problem that two vertexes along the paths are connected by two arcs, and chooses the wrong arc for output for one of the paths. Did someone else had such problems with pgrouting? Cheers, Michael


St_buffer width issues

I am using postgres with postgis for spatial data and ArcGIS for visualisation. I am really confused regarding what units the st buffer function uses because when I insert just the value 1 as measurement to be buffered around a line, the line becomes very wide. Then when I enter a value like 0.00001 the line width would be small. Is there a way where I change or set the units of the value? How can I have more control on the value I want the line to be wide? Since I am building roads with various width values.


GeoServing thousands of rasters, is there a limit on the number of images?

I have a web application on which the users do some processing that outputs rasters. Those rasters are transferred to geoserver for later visualization through WMS. The GeoServer is running on TomCat6, this geoserver is feeded daily with hundreds of ~2mb rasters. A few days ago it reached ~3000 datastores (and layers) and it stopped responding every time a new raster is added, then it needs to be rebooted. I followed the optimization guideline on the geoserver documentation and changed to Oracle Java and set the container options (more java memory and throughput garbage collector). That solved the problem and it's running very fast. But the number of rasters continues to increase and my question is if there is a limit on the number of raster stores that geoserver can handle so that I should expect to have the same problem in the future.


Problems with the scalebar not CRS related in QGIS

When I add one layer to the map, open it in the print composer and add a scale bar the scale bar is fine. I can add layers (with the same CRS) and everything is fine.   Once there are two layers in my map and I open a new print composer and add a scale bar it comes out wrong (just a 0 with a line under it). This is still the case when using the same two layers that I had in the original print composer. Am I doing something wrong or is this a bug? Rachel


How to add a floating "crosshairs" icon above leaflet map?

While reading the article about MapBox's map utility to searching for flight MH370, I noticed that they have a stationary "floating" target (or crosshairs) icon that helps users identify the center of the map, presumably to help them get the exact coordinates they are intending to report. I've searched our site, but cannot find anything that relates to how to create this object that (I assume) just sits above the map frame. How would one add a target (or crosshairs) icon centered on their own Leaflet map?


How to correlate Image DN with Spectrometer reflectance?

I am trying to develop a system to measure NDVI from the commercial grade digital camera. I am trying to get two bands NIR and RED from the image and apply it for NDVI measurement using diffent filters. I also want to apply the Spectrometer reading on field to validate the NDVI measured from the Images. Can anyone help me, how to correlate spectrometer reading(Normally Reflectance of the surface) and the image DN(Three Layer Image) ??


Coordinate mix up

I am working on a basic map, pulling coordinates from google earth, converting them to decimal degrees in excel, and then importing them into qgis. I am using qgis 1.8. When converting from degrees, minutes, seconds I get a coordinate in decimal degrees that looks like this (37.01879,-77.15026). This seems to be the correct gps coordinate for this given address. However, I imported my excel file as a delimiyed csv file, and found that none of the coordinates matched up. My project in qgis vis projected in wgs 84. I manually input the location in qgis, and the coordinates that were given back to me were this (11776056.2,3715032.7). I thought this might be because map units were set to meters, wrong projection, etc. But I cannot seem to figure out the issue. I have the same problem any time I try to take coordinates from Google Earth.


Activate the Geography Page programatically in ArcPad 10

In ArcPad Studio, when you are working with the axf file of a feature layer, there is an attribute for labeled "geographypagevisible". This determines whether or not the geography page is visible on an edit form. There is a way to programatically activate a page using Object.Activate, but I cannot figure out how to activate the geography page due to some event. Line 5 demonstrates what I want to achieve, but it does not work. How do I activate the Geography page programatically?   1. Sub GetId()   2. Dim objLayer, objForm   3. Set objLayer = Map.SelectionLayer   4. Set objForm = objLayer.Forms("EDITFORM")   5. objForm.Pages("geographypage").Activate   6. End Sub Thanks!


Converting data from gml.gz to ESRI format

I am using the Productivity Suite v3 Data Converter tool to convert Urban Paths (from OS Master Map ITN layer (UK)) from gml.gz format to ESRI format. However, there appears to be a bug in this application as it is impossible to select the Urban Paths layer within OS Master Map ITN. Also, none of the windows displays properly. I have uninstalled and reinstalled Productivity Suite v3 but the tool still doesn’t work. I have used it in Productivity Suite v2.1 with no problems but now that I have ArcGIS v10.1 I am stuck with Productivity Suite v3 as v2.1 isn’t compatible with ArcGIS v10.1. InterPOSE doesn't work for my data and Interpose is unable to provide technical support as I am a PhD student and get the software free. Are there other ways to decompress and convert gml.gz data to ESRI format? I hope you can help. Many thanks.


Measurement not correct in QGIS 1.8

I have a SHP file of sub-areas of U.S. counties (roughly about the size of a Census tract) that I've reprojected into a 'true distance' projection (World Equidistant Cylindrical) in QGIS 1.8. However, when I measure the length of the sub-area border segments I get results that don't make sense (3 meters for instance). I've done this before in earlier versions of QGIS and my results seemed fine. Any idea what the problem could be? Thanks for your help.


How upload KML from your local disk using Google Maps API?

I am trying to get load my KML from my local disk with given below code :               <script type="text/javascript">          function initialize() {         var mapOptions = {           center: new google.maps.LatLng(43.321675,-75.172631),           zoom: 8,           mapTypeId: google.maps.MapTypeId.ROADMAP         };         var map = new google.maps.Map(document.getElementById("map_canvas"),mapOptions);         var ctaLayer = new google.maps.KmlLayer('C:/Users/Ashutosh_k/Downloads/newyork.kml');         ctaLayer.setMap(map);       }     </script>      Problem : When my google map is displaying its not showing up the elements of KML .


How to align/georeference distorted vector layers in QGIS?

I've got a shapefile which has been digitized from a strongly distorted historical map, and would like to align it as nicely as possible to another precise shapefile layer. The transformation must not be affine : since the map is historical, it has both cartographic and measurements errors. The best I've achieved for now is to use :               ogr2ogr -gcp 01_a_x 01_a_y 01_b_x 01_b_y -gcp 02_a_x 02_a_y 02_b_x 02_b_y ... -tps      The results is not bad at all (see images below), but the deformation is not "precise" enough. I mean that the points I defined do not align perfectly in the resulting geometry. I would like to have those point align perfectly. Of course, points that are not set using -gcp would not have the same precision since they are only interpolated. Is there a way to achieve that ? Thanks a lot !! Olivier Step 1 : misaligned layer in red, reference layer in gray ![Step 1](http://i.stack.imgur.com/AtAO6.png) Step 2 : points that I defined using -gcp ![Step 2](http://i.stack.imgur.com/FX0qU.png) Step 3 : resulting layer. Not bad, but as you see, the points I defined are not perfectly aligned... ![Step 3](http://i.stack.imgur.com/Ig8yJ.png)


ArcGis Javascript api detect one second hover

For reference I am using arcgis_js_v33 api. I am trying to figure out a good way to detect if a user has hovered over the map for at least one second and if they have to pop up the info box identifying what is under their cursor. My issues seems to be that I cannot find a good way to start a timer if they stop moving their mouse and to kill the timer if they move their mouse before the timer reaches its time. Any ideas?


How to Convert Multiple Raster in to NettCDF in Arcmap 10.1?

My name is Paresh Chaudhary and I've worked in the field of Geometics since 2008. I have 30 day by day surface temperature images of MODIS, aqua satellite. Now i want to convert it into a NetCDF raster in ArcMap 10.1. I know one tool, **_Raster to NetCDF_** in **Multidimension Tools** tool set which does this but you can't input multiple images. So is it possible to convert 30 temperature variables into NetCDF raster? I would like to visualise all images with time and also perform temporal analysis. Thanks a lot in advance


Configuring pgRouting on Win7

I'm installing pgrouting (ver 1.03_pg-8.4.2) on my windows machine. I installed PostgreSQL and postgis without any problem, then I downloaded the binary package of pgrouting and I installed it (copy and Paste in postgresql dir). Now I'm trying to configure my first database following the instructions of http://genuchten.blogspot.com/2008/04/untitled.html but I get some errors: 1)               psql -U postgres -f ../share/contrib/lwpostgis.sql routing     ERROR: ../share/contrib/lwpostgis.sql No such file or directory      2)               psql -U postgres -f ../share/contrib/spatial_ref_sys.sql routing     ERROR: ../share/contrib/spatial_ref_sys.sql No such file or directory      3)               psql -U postgres -f ../share/contrib/routing_core.sql routing     ERROR: incompatible library: librouting.dll: verion mismatch      4)               psql -U postgres -f ../share/contrib/routing_core_wrappers.sql routing     ERROR: Type "geoms" does not exist      Can anyone help me? Thanks Laura