title: PostGIS - Geometry or Geography find points within distance in meters

text: I have a requirement where I need to store lat/long coordinates for locations all around the world. I then need to query for those locations that are within a specified distance (in meters) from another location. These distances are usually within the 500-1000 meter range so I don't need crazy precision. Following various examples I've found, I've been using SRID 4326 with my geometry column in a PostGIS database. The column simply holds lat/long coordinates. The ST_DWithin method takes a distance in meters for geography types, but I don't know what the unit is for 4326. It doesn't seem to be meters with some queries that I tried. I'm assuming it's degrees. I've read that for performing Spatial queries on small areas, it's best to use geometries since operations are quicker at the price of precision. How can I use these methods to provide meters for either `ST_DWithin` or `ST_Distance`?


title: Efficient way to store multiple values of an x,y location in Python?

text: I am trying to read values from a stack of images and store them in a structure that will allow for easy readability. My current idea is to create a loop               for x in width:         for y in length:             for z in height:      and try appending the z value to a list at [x, y]. The issue with this technique is that I do not know how to prevent losing the x and y position through the iteration and it seems rather inefficient. The input data will consist of z values extracted from tiff images. I am reading them by first creating a stack through a loop statement               for i, fname in enumerate(image_files):     img = cv2.imread(fname, cv2.CV_LOAD_IMAGE_UNCHANGED)     img_stack[:,:,i] = img      which I found here: http://stackoverflow.com/a/20056062/1342580 (CV2 is only returning nan for me but that's a different question) After the stack is created I used to loop at the beginning of my question to go through all of the pixels (which at the moment is 3858 x 3742) and append the z value of an x,y location from each image to the x,y value (for example if I had 4 images in the stack the result at position 1, 1 would be [1, 1, Z1, Z2, Z3, Z4]. Then I append that to another list that will contain a list for every x,y position Is there a better technique for doing this?


title: Cannot access Earth layer in Google Maps V3 API?

text: The "new" google maps has the awesome feature of accessing the earth globe in the browser (on Mac and Windows), without the need for installing the google earth plugin. I know that you can access the new look and feel of the maps in Version 3 of the API using google.maps.visualRefresh global property set to true. However I have not been able to access the earth globe. I had a look at the available mapTypes and there is no "earth", or "satellite_3d", and nothing really in the documentation which leads me to think it is available. Does anyone know if this is at all possible? And if so how do you go about it? Thanks Ando


title: ArcGIS Javascript example for geometry services?

text: I have a client that wishes to show him all the geometry services example in Javascript API http://help.arcgis.com/EN/ARCGISSERVER/10.0/APIS/REST/geometryserver.html ... I'm new to Javascript and ArcGIS, is there a sample with all the examples for geometry services? https://developers.arcgis.com/en/javascript/jssamples/util_buffergraphic.html Yes I have been here https://developers.arcgis.com/en/javascript/ but they don't have all example.


title: 'Example for displaying attributes' Unable to see attribute data

text: 'Example for displaying attributes' available with geoExt package are not working in my system. What could be the possible reason. I can see attribute table in web page windows, but no attribute entry is there. Is there some proxy setup issue? I have not done any proxy setup.


title: What's the best way to add a projection to GDAL?

text: I have a number of files in Maryland State-Plane (US foot) that have the following wkt: "PROJCS["NAD_1983_StatePlane_Maryland_FIPS_1900_Feet",   GEOGCS["NAD83",   DATUM["North_American_Datum_1983",   SPHEROID["GRS 1980",6378137,298.2572221010002,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],   PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4269"]],   PROJECTION["Lambert_Conformal_Conic_2SP"],   PARAMETER["standard_parallel_1",38.3],   PARAMETER["standard_parallel_2",39.45],   PARAMETER["latitude_of_origin",37.66666666666666],   PARAMETER["central_meridian",-77],   PARAMETER["false_easting",1312333.333333333],   PARAMETER["false_northing",0],   UNIT["US survey foot",0.3048006096012192,AUTHORITY["EPSG","9003"]]]" OGRSpatialReference::autoIdentifyEPSG fails to recognize this. As far as I can tell, the only difference in this wkt and the one for 2248 on spatialreference.org is some trailing significant digits, some omitted Authority codes, and a reversal in the order of the standard parallels. [Edit] I've been looking at the pcs.csv file and while I don't know exactly how to read the columns, it appears that the parameters in the file for 2248 don't match the above wkt or the wkt on spatialreference.org. The standard parallels differ, as does the latitude of origin. This ups the ante on what the right action is. What is the preferred method for adding this to the data files so that it will be recognized? Do I edit my pcs.csv file? The pcs.override.csv? Or something else? Changing the content of the file is impractical in this case.


title: Ubuntu Newbie I am frustrated trying to install (all things GEOS) WEB-GIS

text: I am new to ubuntu SAUCY and I must confess the transition from windows has not been really easy for me. on windows everything seems pretty much straight forward but I am getting frustrated because It seems when I try to install either pyproj, gdal, mapnik, proj4, Postgre/postgis and other Web gis packages they probably don't work together and their dependencies aren't the same. my question is isn't there an easy way to get all these packages to work together? when I install QGIS it removes something else same goes to mapnik. I install mapnik the other day it even crash and messed up my system, had to do a clean OS install all over again. here's the error im getting now.               libmapnik-dev: Depends: libmapnik (= 2.2.0-ubuntu1~saucy1) but 2.2.0-ubuntu1~saucy1 is to be installed                    Depends: mapnik-utils (= 2.2.0-ubuntu1~saucy1) but 2.2.0-ubuntu1~saucy1 is to be installed                    Depends: libpng-dev but it is a virtual package                    Depends: libtiff-dev but it is a virtual package                    Depends: libltdl7-dev but it is a virtual package                    Depends: libicu-dev (>= 4.2) but 4.8.1.1-12ubuntu2 is to be installed      here is the last link(help) have have used. https://github.com/mapnik/mapnik/wiki/UbuntuInstallation all the GIS package I have tried installing is throwing this same error message. I really need help on how to go about getting these installed because I cant work without all of them. Thanks Anthony.


title: Where can I get detailed Bathymetry data for Ontario

text: I need detailed Bathymetric data covering water bodies (including smaller ones) in the Province of Ontario. Our Provincial Government has it, but the use of the data are restricted due to liability issues. Does anyone know place to download this data, or can anyone suggest a way to interpolate lake depth contours (ie - would it be feasible to get a Google Earth Pro licence for this)? I have access to 3D Analyst and Spatial Analyst. Thanks in advance for your suggestions


title: Programmable buffers using GDAL and python

text: Is there a way to programmatically generate buffer layers for vector layers using GDAL and python?


title: Selection problem in WMS layer using OpenLayers

text: I added 4 WMS layers (three polygons and one have point geometry) from GeoServer on my webpage My requirement is:   1. I have to select a single polygon of one wms layer using a mouse click. After the selection, the polygon color should be changed. If you select another polygon, the previous selected polygon color is to be changed to original color.   2. If I select a polygon of one wms layer, I want to get the `NAME` of the selected polygon as an alert. Any help is greatly appreciated.. Thank you..               layer = new OpenLayers.Layer.WMS(           "Layer", "http://********:8080/geoserver/POC/wms",                      {                     LAYERS: 'POC:CIRCLE_VIEW',                     STYLES: '',                      wrapDateLine: true,                      format: format,                     transparent:true                      },                 {                    singleTile: true,                     ratio: 1,                     isBaseLayer: false,                    yx : {'EPSG:42106' : false}                      }                  );       ` I added four layers in this way..               var infoControls = {              click: new OpenLayers.Control.WMSGetFeatureInfo({                  url: 'http://********:8080/geoserver/POC/wms',                   title: 'Details',                  layers: [layer,layer2,layer3,layer4],                  queryVisible: true              }),              hover: new OpenLayers.Control.WMSGetFeatureInfo({             url: 'http://********:8080/geoserver/POC/wms',                  title: 'Details',                layers: [layer,layer2,layer3,layer4],                  hover: true,                  queryVisible: true              })          };         for (var i in infoControls)          {                  infoControls[i].events.register("getfeatureinfo", this, showInfo);             map.addControl(infoControls[i]);          }             infoControls.click.activate();     function showInfo(evt) {                 if (evt.features && evt.features.length) {                           highlightLayer.destroyFeatures();                      highlightLayer.addFeatures(evt.features);                      highlightLayer.redraw();                 } else {                     highlightLayer.destroyFeatures();                      highlightLayer.addFeatures(evt.features);                      highlightLayer.redraw();                     document.getElementById('responseText').innerHTML = evt.text;                 }             }      In this I am showing attribute data of selected polygon.


title: How to find SRID for random polygon?

text: I'm testing spatial database performance on a dataset. Spatial objects in this dataset (stored as WKT in a flat text file) are random valid polygons which have vertex coordinates (x,y) where 100,000> x >0 and same for y. Currently I'm having problem with inserting these data into PostGIS which has a geometry constraint check and I got an error like "violates check constraint "enforce_srid_polygon ". So, how can I find/create a srid for my dataset ? BTW: the polygons are valid polygons which are extracted from images by segmenting objects.


title: How to preserve the original direction of NHD segments when creating a route using the gp.CreateRoutes_lr tool?

text: I am trying to figure out how to preserve the original line direction of NHD High Resolution data that I am selecting out and creating a route with. I'm using ArcGIS 9.3.1 and this is a python script. I am selecting out some of the NHD streams by the ReachCode field into a new feature class (I have also tried this with a feature layer in the context of the issue I'm posting about but creating a separate feature class allowed me to debug further), and then am taking that feature class and dissolving it. I then take that dissolved NHD section and am creating a route with the gp.CreateRoutes_lr tool. I've narrowed down the issue to the step when the script creates the route and it seems that the direction of the line feature gets flipped when I create the route (I have definitively verified this is happening). Is this what is supposed to happen? I've never worked with the linear referencing tools before this so am unsure if this is what should be happening. Each of the line segments have a defined FROM and TO node direction which is important for streams in terms of correct flow direction. The reason the direction matters to me for my script is that I am creating the vertices via the gp.FeatureVerticesToPoints_management tool get the elevations of the FROM and TO node to determine the stream slope for points located along the NHD section. This impacts the slope calculation obviously. I just compared a site that had a negative slope (with assumption that a negative slope implies an instance where the route direction got flipped as I describe) with a site that had a positive slope (with assumption that a positive slope indicates an instance where the route direction did not get flipped). Both assumptions were true, and it seems for the former the direction got flipped, and the latter it did not. So it seems that this tool does not flip the directions consistently... Is there any way to preserve the original direction of the NHD segments when I create a route? Or, what would influence the direction of the output route line feature?? thanks, Tom


title: How do I add precise border with longitude and latitude for any country using mapinfo?

text: I have a plain image that I need to add international border with precise longitude and latitude but I have no experience with Mapinfo. How is it possible? Is there any way to import maps from Bing or Google Earth with that information?


title: Adjust Picture Marker Symbols with Python? (ArcGIS 10 SP1)

text: I am working toward a map cache, and therefore have about 12 versions of the same datasets to be shown at different scales. Invariably, the symbol sizes I have set on highway shields (for example) need to be fussed with in order to make them look right. In order to change the symbol size (scale to text size doesn't fit my application), you have to click four times from the label manager, and you have to do it over and over and over until all eight label classes are done. And then you want it just a tad larger.... Can python (or anything else) be used to adjust picture marker symbol size across multiple label classes? Thank you!


title: Find distance from a point to linestring

text: I have table (shapetable) with my coastline shape file loaded. I have another table (pointtable) with lat long points. Wondering how would i go about writing a sql query to find the distance between the lat long points and the shapefile. The shapetable has a geography type attribute, from my reading it appears that geometry type linestrings and geography are one and the same. Appreciate the help thanks vic **EDIT:** Being relatively a newbie at this, I'm not using PostGIS. I am using sql server 2008. I imported the shape file using a shape2sql utility. The table has a bunch of shapefile attributes and also the geography data type. If i understand this correctly, the geography binary data is nothing but linestring of geometry data type(I'd verify this by running a stgeoemtrytype() on the column) I have another table with a bunch of lat long. I want to implement a solution where by i want to get distance from the shapefile, using stdistance. One confusing part i found in going through a bunch of documents was that, given a geography data type, how could i find out the linestring coordinates. Thanks again


title: How to stop values from rounding to whole numbers in ArcGIS attribute tables?

text: I am trying to edit an attribute table in ArcGIS 10.1 and so far have been unable to keep it from rounding to whole numbers. Even when I manage to "pad with zeroes" using the numeric settings, it will round to 6.0. I have tried the Editor > Options > General method as well. My decimal places are currently set to round to 6, but any data I enter rounds to a whole number. My field is set to "Long Integer". What obvious solution am I missing?


title: QGIS: On the fly projection from WGS84, 19 N, to Web Mercator. Do I need to define a new CRS with transformation parameters?

text: Workin in QGIS, I have a shapefile called "protected areas" of Dominican Republic, that was digitized from topomaps 50k. Shapefile is projected to WGS84, 19N zone. My problem is that when I load openlayers, the limits of almost all protected areas, don't get aligned with openlayers. In the following image there is an example:![Image 1](http://www.geografiafisica.org/wp-content/uploads/2012/12/screenshot1.jpg) The limit of the protected area (in red) is the coastline; but the real coastline plotted in Openlayers (google satellite) doesn't align with my shapefile (this area has rocky coastline, that haven't changed since topo was made). There is a ca. 45m shift in the X axis, and ca. 50m in the Y axis. In the following image: ![Image showing Shift](http://www.geografiafisica.org/wp- content/uploads/2012/12/screenshot2.jpg) the limit of the protected area is prefectly aligned with topomap (coastline, road, etc.). As I know, OpenLayers is projected using "Web Mercator". So, do I need to define a new CRS with transformation parameters between WGS84 19N and Web Mercator? I believed they use the same ellipsoidal reference, so I shouldn't need to define a transformation. Any help would be appreciated. Thanks. Jose


title: How to perform a spatio-temporal query in attribute table?

text: I am not good at programming, so I would ask for your understanding concerning my question. I have an attribute table that contains five columns which are: (Date, ID, Latitude, Longitude, Rainfall). Rainfall for different periods of time (which are defined in column "Date") have been recorded in certain stations. In the column "ID" is the ID of the station on which the recording has been made for a specific period of time. Every station has a unique ID, but at the same location, rainfall has been recorded for different time periods. On the column "Rainfall" appears the amount of rain that was recorded on the specific time and the specific station. My question is, how can I get the total rainfall in a certain station that has been recorded, using all the records in column "Date" that refer to that specific station. Could I somehow use the coordinates of its station to get the total rainfall that has occurred on that specific station? The software I use is QGIS 2.4 on Windows machine.


title: Why does QGIS crash when merging shapefiles?

text: This is a repeatable problem. I have created a shapefile using a csv. I am trying to merge it with an existing shapefile, both of which contain only points. (Using the Vector | Data Management | Merge Shapefiles option.) I have double-checked that the shapefile type is correct, and the coordinate reference system matches. As soon as I click OK to accept the result, the message appears "qgis.exe has stopped working" etc., the usual Windows crash message. Things I have done to reslove this: 1. I have a friend who has the same software and the same project files (we are collaborating on this project). I sent him my original csv file with the new data in it, and he was able to successfully merge the files at his end. So I know that the csv file is OK. 2. Uninstalled and re-installed the program. Identical result, program crashes every time as soon as I go to accept the merging of the shapefiles. Using Windows 7. Does anyone know how to resolve this issue? If not, is there another way to bring new data into the project? I don't want to use the "Add Feature" option as it is too difficult to get the co-ordinates correct with just hovering over the map.


title: Returning string from OpenLayers.Format.WKT write?

text:
How to return a string from OpenLayers.Format.WKT write               this.innerHTML = new OpenLayers.Format.WKT().write(         layer.selectedFeatures[0]     );      my full code is like this please suggest me...               //OpenLayers.ProxyHost = "proxy.cgi?url=";     OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";          var wps = "http:",         capabilities, // the capabilities, read by Format.WPSCapabilities::read         process; // the process description from Format.WPSDescribeProcess::read          // get some capabilities     getCapabilities();          // create the UI     var layer = new OpenLayers.Layer.Vector("Scratchpad");     var toolbar = new OpenLayers.Control.EditingToolbar(layer);     toolbar.addControls([new OpenLayers.Control.ModifyFeature(layer, {         title: "Select feature"     })]);     var map = new OpenLayers.Map('map', {         controls: [             toolbar,             new OpenLayers.Control.ZoomPanel(),             new OpenLayers.Control.PanPanel()         ],         layers: [             new OpenLayers.Layer.WMS(                 "OSM", "",                 {layers: "openstreetmap", format: "image/png"}             ), layer         ]     });     map.zoomToMaxExtent();          // add behavior to html elements     document.getElementById("processes").onchange = describeProcess;          // using OpenLayers.Format.WPSCapabilities to read the capabilities     function getCapabilities() {         OpenLayers.Request.GET({             url: wps,             params: {                 "SERVICE": "WPS",                 "REQUEST": "GetCapabilities"             },             success: function(response){                 capabilities = new OpenLayers.Format.WPSCapabilities().read(                     response.responseText                 );                 var dropdown = document.getElementById("processes");                 var offerings = capabilities.processOfferings, option;                 // populate the dropdown                 for (var p in offerings) {                     option = document.createElement("option");                     option.innerHTML = offerings[p].identifier;                     option.value = p;                     dropdown.appendChild(option);                 }             }         });     }          // using OpenLayers.Format.WPSDescribeProcess to get information about a     // process     function describeProcess() {         var selection = this.options[this.selectedIndex].value;         OpenLayers.Request.GET({             url: wps,             params: {                 "SERVICE": "WPS",                 "REQUEST": "DescribeProcess",                 "VERSION": capabilities.version,                 "IDENTIFIER": selection             },             success: function(response) {                 process = new OpenLayers.Format.WPSDescribeProcess().read(                     response.responseText                 ).processDescriptions[selection];                 buildForm();             }         });     }          // dynamically create a form from the process description     function buildForm() {         document.getElementById("abstract").innerHTML = process["abstract"];         document.getElementById("input").innerHTML = "
Input:
";         document.getElementById("output").innerHTML = "";              var inputs = process.dataInputs, supported = true,             sld = "text/xml; subtype=sld/1.0.0",             input;             //var change=inputs.innerHTML.toFixed(2);         for (var i=0,ii=inputs.length; i<ii; ++i) {         alert("Hello! I am an alert box!!");             input = inputs[i];             if (input.complexData) {             //alert("Hello! I am an alert box in complexdata !!");                 var formats = input.complexData.supported.formats;                 if (formats["application/wkt"]) {                 //alert("Hello! I am an alert box in app/wkt!!");                     addWKTInput(input);                 } else if (formats["text/xml; subtype=     -collection/1.0"]) {                 //alert("Hello! I am an alert box in text/xml!!");                     addCollectionInput(input);                 } else if (formats["image/tiff"]) {                         // addRasterInput(input);                     //alert(truncator(5.1231231, 1));                     //alert("Hello! I am an alert box in image/tiff!!");                     addRasterInput(truncator(input,2));                 } else if (formats[sld]) {                 //alert("Hello! I am an alert box in formatssld!!");                     addXMLInput(input, sld);                 } else {                     supported = false;                 }             } else if (input.boundingBoxData) {             alert("Hello! I am an alert box in bounding box!!");                  // var change=input.innerHTML.toFixed(2);                 addBoundingBoxInput(input);             } else if (input.literalData) {             //alert("Hello! I am an alert box in literal box!!");                 addLiteralInput(input);             } else {                 supported = false;             }             if (input.minOccurs > 0) {             //alert("Hello! I am an alert box in minoccurs box!!");                 document.getElementById("input").appendChild(document.createTextNode("* "));             }         }              if (supported) {             var executeButton = document.createElement("button");             executeButton.innerHTML = "Execute";             document.getElementById("input").appendChild(executeButton);             executeButton.onclick = execute;         } else {             document.getElementById("input").innerHTML = '' +                 "Sorry, the WPS builder does not support the selected process." +                 "";         }     }          // helper function to dynamically create a textarea for geometry (WKT) data     // input     function addWKTInput(input, previousSibling) {     alert("Hello! I am an alert box in dynamically create a textarea for geometry data input!!");         var name = input.identifier;         var container = document.getElementById("input");         var label = document.createElement("label");         label["for"] = name;         label.title = input["abstract"];         label.innerHTML = name + " (select feature, then click field):";         previousSibling && previousSibling.nextSibling ?             container.insertBefore(label, previousSibling.nextSibling) :             container.appendChild(label);         var field = document.createElement("textarea");         field.onclick = function () {         alert("Hello! I am an alert box in onclick!!");             if (layer.selectedFeatures.length) {             alert("Hello! I am an alert box in selected features function for geometry data input writing!!");                 this.innerHTML = new OpenLayers.Format.WKT().write(                          //layer.selectedFeatures[0]                     //layer.selectedFeatures[0]                     //name.length();                     //alert(trun.toFixed(2));                     //layer.selectedFeatures[0]                     layer.selectedFeatures[0].geometry.toString()                      );             }             createCopy(input, this, addWKTInput);         };         field.onblur = function() {             input.data = field.value ? {                 complexData: {                     mimeType: "application/wkt",                     value: this.value                 }             } : undefined;         };         field.title = input["abstract"];         field.id = name;         previousSibling && previousSibling.nextSibling ?             container.insertBefore(field, previousSibling.nextSibling.nextSibling) :             container.appendChild(field);     }          // helper function for xml input     function addXMLInput(input, type) {         var name = input.identifier;         var field = document.createElement("input");         field.title = input["abstract"];         field.value = name + " (" + type + ")";         field.onblur = function() {             input.data = field.value ? {                 complexData: {                     mimeType: type,                     value: this.value                 }             } : undefined;         };         document.getElementById("input").appendChild(field);     }          // helper function to dynamically create a  collection reference input     function addCollectionInput(input) {         var name = input.identifier;         var field = document.createElement("input");         field.title = input["abstract"];         field.value = name + " (layer on demo server)";         addValueHandlers(field, function() {             input.reference = field.value ? {                 mimeType: "text/xml; subtype=-collection/1.0",                 href: "http://",                 method: "POST",                 body: {                     : {                         version: "1.0.0",                         outputFormat: "GML2",                         featureType: field.value                     }                 }             } : undefined;         });         document.getElementById("input").appendChild(field);     }          // helper function to dynamically create a raster (GeoTIFF) url input     function addRasterInput(input) {         var name = input.identifier;         var field = document.createElement("input");         field.title = input["abstract"];         var url = window.location.href.split("?")[0];         field.value = url.substr(0, url.lastIndexOf("/")+1) + "data/tazdem.tiff";         document.getElementById("input").appendChild(field);         (field.onblur = function() {             input.reference = {                 mimeType: "image/tiff",                 href: field.value,                 method: "GET"             };         })();     }          // helper function to dynamically create a bounding box input     function addBoundingBoxInput(input) {         var name = input.identifier;         var field = document.createElement("input");         field.title = input["abstract"];         field.value = "left,bottom,right,top (EPSG:4326)";         document.getElementById("input").appendChild(field);         addValueHandlers(field, function() {             input.boundingBoxData = {                 projection: "EPSG:4326",                 bounds: OpenLayers.Bounds.fromString(field.value)             };         });     }          // helper function to create a literal input textfield or dropdown     function addLiteralInput(input, previousSibling) {         var name = input.identifier;         var container = document.getElementById("input");         var anyValue = input.literalData.anyValue;         // anyValue means textfield, otherwise we create a dropdown         var field = document.createElement(anyValue ? "input" : "select");         field.id = name;         field.title = input["abstract"];         previousSibling && previousSibling.nextSibling ?             container.insertBefore(field, previousSibling.nextSibling) :             container.appendChild(field);         if (anyValue) {             var dataType = input.literalData.dataType;             field.value = name + (dataType ? " (" + dataType + ")" : "");             addValueHandlers(field, function() {                 input.data = field.value ? {                     literalData: {                         value: field.value                     }                 } : undefined;                 createCopy(input, field, addLiteralInput);             });         } else {             var option;             option = document.createElement("option");             option.innerHTML = name;             field.appendChild(option);             for (var v in input.literalData.allowedValues) {                 option = document.createElement("option");                 option.value = v;                 option.innerHTML = v;                 field.appendChild(option);             }             field.onchange = function() {                 createCopy(input, field, addLiteralInput);                 input.data = this.selectedIndex ? {                     literalData: {                         value: this.options[this.selectedIndex].value                     }                 } : undefined;             };         }     }          // if maxOccurs is > 1, this will add a copy of the field     function createCopy(input, field, fn) {         if (input.maxOccurs && input.maxOccurs > 1 && !field.userSelected) {             // add another copy of the field - we don't check maxOccurs             field.userSelected = true;             var newInput = OpenLayers.Util.extend({}, input);             // we recognize copies by the occurrence property             newInput.occurrence = (input.occurrence || 0) + 1;             process.dataInputs.push(newInput);             fn(newInput, field);         }     }          // helper function for adding events to form fields     function addValueHandlers(field, onblur) {         field.onclick = function() {             if (!this.initialValue) {                 this.initialValue = this.value;                 this.value = "";             }         };         field.onblur = function() {             if (!this.value) {                 this.value = this.initialValue;                 delete this.initialValue;             }             onblur.apply(this, arguments);         };     }          // execute the process     function execute() {         var output = process.processOutputs[0];         var input;         // remove occurrences that the user has not filled out         for (var i=process.dataInputs.length-1; i>=0; --i) {             input = process.dataInputs[i];             if ((input.minOccurs === 0 || input.occurrence) && !input.data && !input.reference) {                 OpenLayers.Util.removeItem(process.dataInputs, input);             }         }         process.responseForm = {             rawDataOutput: {                 identifier: output.identifier             }         };         if (output.complexOutput && output.complexOutput.supported.formats["application/wkt"]) {             process.responseForm.rawDataOutput.mimeType = "application/wkt";         }         OpenLayers.Request.POST({             url: wps,             data: new OpenLayers.Format.WPSExecute().write(process),             success: showOutput         });     }          // add the process's output to the page     function showOutput(response) {         var result = document.getElementById("output");         result.innerHTML = "
Output:
";         var features;         var contentType = response.getResponseHeader("Content-Type");         if (contentType == "application/wkt") {             features = new OpenLayers.Format.WKT().read(response.responseText);         } else if (contentType == "text/xml; subtype=-collection/1.0") {             features = new OpenLayers.Format.T.v1_0_0().read(response.responseText);         }         if (features && (features instanceof OpenLayers.Feature.Vector || features.length)) {             layer.addFeatures(features);             result.innerHTML += "The result should also be visible on the map.";         }         result.innerHTML += "<textarea>" + response.responseText + "</textarea>";     }


title: High Resolution Mineral Maps

text: Australia has recently published high resolution surface mineral maps using data from ASTER see here. Does anyone know of similar free products available for the entire globe?


title: Area of an irregular polygon using side lengths

text: I have only lengths for the sides of an irregular polygon, can anyone tell me how I can measure the area of the polygon? Remember only lengths of all the sides , no angles or coordinates. Few forums mention about trangulation of the polygon etc But I only have side lengths. Does anybody has any feedback? Thanks.


title: Save folder connection in ArcGIS 10

text: In ArcCatalog one has the option to save locations of folders via the "Connect To Folder" button. I have the problem, that ArcGIS don`t save each folder connection which I would like to have in Catalog. Everytime I fail when I want to save the path **U:\GROUPS\GIS\Geo_Data** or **U:\GROUPS\GIS**. Whereas **U:\GROUPS** or **C:\home\data** is no problem. Anyone encountered the same problem?


title: How to cluster points based on both proximity to each other and proximity to another layer (roads)?

text: Using QGIS I have managed to created cluster for points based on their proximity of each other (using buffer and spatial joins etc), but I would also like to add another variable to it: proximity to a road (another specific layer). Is there any idea if this can be done with plugins and/or an algorithmic procedure of some sort? Three attached images describe the problem.   1. First image is the points (green dots) that I want to cluster, depending on their proximity to each other and their proximity to a road (black/white lines).   2. Second image is the result I get, which is only one big cluster, entirely depending on the distance between all points.   3. The third image is the result I want: two clusters, separated because of their proximity to the "closest coherent/connected road" I first thought of some solution with adding a buffer on the road layer, but that will probably not work, partly because the road lines are very often segmented. ![image 1](http://i.stack.imgur.com/FmdCS.png) ![image 2](http://i.stack.imgur.com/6weBs.png) ![image 3](http://i.stack.imgur.com/h4rA1.png)


title: Which software to use for designing city plans?

text: Which software would you suggest I use for designing city plans? The result needs to be an image (jpeg) so georeferencing is not necessary. I have options like ArcGIS, AdobeInDesign, Adobe Illustrator. I don't know which software is the best for maps like these in example: ![enter image description here](http://i.stack.imgur.com/tZVcC.jpg) Also, I would need to add some 2d effects on the buildings. Thank you for suggestions.


title: Why can't I open QGIS projects via windows explorer?

text: We have upgraded from QGis 1.8 to QGis 2.2. Qgis 2.2 is first installed, then 1.8 removed. Since then it is impossible to open .qgs-projects directly from windows explorer. We have checked explorer, file extensions and relinked to C:\Program Files\QGIS Valmiera\bin\qgis.bat . But it doesn't help. We get the error 'qgis_core.dll cannot be found, please reinstall'. But that is rather drastic. I use windows xp.


title: SLD for OSM data in GeoServer

text: > **Possible Duplicate:**   >  Standard SLDs for OpenStreetMaps? I can't find a good SLD for OSM. I downloaded the OSM data and imported it to PostGIS. With GeoServer I use it as a data source to create my layers. Is there no open source SLD for OSM out there?


title: reproject vector layer with pyqgis with new CRS does not worl

text: In QGIS I am running a script in th python toolbar to export all visible vector layers to a GEOJSON with this snippet:               canvas = qgis.utils.iface.mapCanvas()     allLayers = canvas.layers()     for i in allLayers:               if i.type() != 0 :             print(i.name() + " skipped as it is not a vector layer")           if i.type() == 0 :             exp_crs = QgsCoordinateReferenceSystem(4326,QgsCoordinateReferenceSystem.PostgisCrsId)             print i.crs()             print exp_crs             qgis.core.QgsVectorFileWriter.writeAsVectorFormat(i, '/home/' + str(i.name()) + '.js', 'utf-8', exp_crs, 'GeoJson')      If my vector file has a different CRS it keeps using the vector file CRS and doesn't use the defined CRS EPSG:4326... What am I doing wrong? The whole project context can be found on GITHUB


title: How to use matplotlib with ArcGIS 10.1 without Microsoft Visual C++ Runtime Error?

text: I have problems using the matplotlib for python with ArcGIS 10.1 If I run the code below as a script tool in ArcMap. No pdf is generated and when I close ArcMap this error is shown. ![enter image description here](http://i.stack.imgur.com/9NEDb.png) If I run the code from IDLE a similar error is shown for C:\Python\ArcGIS10.1\pythonw.exe but the pdf is generated correct. If I run the code in a cmd prompt no error code is shown and the pdf is generated correct. Does anybody know how to make the script work as a script tool in ArcMap?               from matplotlib import pyplot          fig = pyplot.figure(figsize=(8.27,11.69))          axis1 = fig.add_subplot(311)     axis1.plot(range(10))          fig.savefig('multipleplots.pdf')


title: Error exporting data from Oracle to shp using ogr2ogr

text: I have an oracle table which I would like to query/make a selection. I would then like to export the selection to a shp file. I'm using gdal and ogr2ogr. Following the guidelines here and here I have come up with the following:               ogr2ogr -f "ESRI Shapefile" c:/Testing/Scratch/GDAL/symbol.shp OCI:dbusername/db_password@tns_name -sql "SELECT* FROM tablename where VERSION_DATE >= CURRENT_DATE -8" -progress      I then get the following error: "Unable to open datasource 'OCI:dbusername/db_password@tns_name' with the following drivers" There is then a long list of drivers including "ESRI Shapefile". I have now added " after OCI as suggested.               ogr2ogr -f "ESRI Shapefile" c:/Testing/Scratch/GDAL/symbol.shp OCI:"dbusername/db_password@tns_name -sql "SELECT* FROM tablename where VERSION_DATE >= CURRENT_DATE -8" -progress      I no longer get the errors about it being unable to open the datasource with the drivers. I now just get an "Access in denied". I have a working connection to the database in ArcGIS using the same username & password. Does gdal generate its own user?


title: Ocean flow models - predict the path of an object dropped in the ocean

text: This question is from a user who wants to calculate the predicted path of a bottle dropped from a ship. He has a spreadsheet of hundreds of these bottles and wishes to display the predicted vs actual landing points of the bottles. Is there a publicly-accessible model which will show the predicted path of an object, anywhere in the earth's oceans? For example, NOAA published a series of maps during the Gulf of Mexico oil crisis, showing the predicted path of the oil spill. I'm presuming that this is a fairly specialised branch of GI-science and the answer will be that such a model depends highly on local variations, the date, wind and other weather conditions, etc. As a next-best option, is information on the general patterns of the ocean currents publicly accessible in a GIS-ready format? It may be possible to create a geoprocessing model from this?


title: Looking for XSL templates using Google Charts API

text: Hi all im searching for examples on usage of Arcgis Desktop html popups in general, being more specific on xsl templates using google charts api on it. Any references and comments are welcome. Thanks Jorge Vidinha


title: Delay on Marker Title in Leaflet

text: I am trying to put title on leaflet markers. This explains the way to give title to marker which would appear on hover. Is there any way we can set time on that i.e. make the title appear after 50ms of the mouse hover?


title: Automate import of multiple rasters as layers into tilemill

text: Could anyone please suggest a way of automating the import of raster tiffs into tile mill? It's driving me mad as I have well over 100 files. I could do it manually but I'm on a tight deadline... I have many warped tiff files ready to be imported as individual layers into tile mill so they can be saved out as a single mbtiles file. Many thanks in advance!


title: Find out whether geography point lies within geography bounds

text: I have three points:   * Point A   * Point SW   * Point NE I want to know if Point A resides within the bounds determined by points SW (for south-west coordinates) and NE (for north-east coordinates). I've so far tried both `ST_Intersects` and `&&` to determine this, which works fine for small regions, but fails whenever the bounding area determined by SW and NE points covers more than half the globe. In this case it seems to instead use a different, smaller, bounding area using the same four points. As an example, if I'd want to cover everything from the SW point of South Africa to the NE point of Ireland it would instead catch points where the the NW point is in Ireland and the the SE point is in Madagascar. That is, China would not considered to be within the bounds. Here are some SQL queries to illustrate this:               select      ST_GeographyFromText('SRID=4326;POINT(0 0)') && ST_GeographyFromText('SRID=4326;POLYGON((180 -80,  -80 -80, -80 80, 180 80, 180 -80))') as expectTrue_1,      ST_Intersects(ST_GeographyFromText('SRID=4326;POINT(0 0)'), ST_GeographyFromText('SRID=4326;POLYGON((180 -80,  -80 -80, -80 80, 180 80, 180 -80))')) as expectTrue_2,      ST_GeographyFromText('SRID=4326;POINT(0 0)') && ST_GeographyFromText('SRID=4326;POLYGON((180 -80,  180 80, -80 80, -80 -80, 180 -80))') as expectTrue_3,      ST_Intersects(ST_GeographyFromText('SRID=4326;POINT(0 0)'), ST_GeographyFromText('SRID=4326;POLYGON((180 -80,  180 80, -80 80, -80 -80, 180 -80))')) as expectTrue_5,      ST_GeographyFromText('SRID=4326;POINT(-90 0)') && ST_GeographyFromText('SRID=4326;POLYGON((180 -80,  -80 -80, -80 80, 180 80, 180 -80))') as expectFalse_1,      ST_GeographyFromText('SRID=4326;POINT(0 0)') && ST_GeographyFromText('SRID=4326;POLYGON((1 -1,  -1 -1, -1 1, 1 1, 1 -1))') as expectTrue_5;      I'd want for this to return `t t t t f t` but it returns `f f f f t t`. I've also tried `ST_Intersects(ST_MakeBox2D())` as per the following SO post: http://stackoverflow.com/questions/24212355/query-by-coordinates-takes-too- long-options-to-optimize/24218861?iemail=1&noredirect=1#24218861 I'm using PostGIS 2.1.3 with the `Geography` data type. What functions or approaches should I use to have PostGIS treat my more-than- half-globe-spanning polygons as more-than-half-globe-spanning ones instead of 'optimizing' them to be less? I understand I could just split the bounds in, say, two parts, and check for each, but there must be a cleaner solution.


title: web map service - does it update by default when data is amended in Oracle source tables?

text: I have a WMS which source data is held in Oracle. Every year an amendment to the view in Oracle is run to update records. Does the web map service need to be re-published once these amendments are made or are these updates to the data automatically reflected in the WMS? Any help would be appreciated.


title: How to Geocode a Google SpreadSheet?

text: I have a spreadsheet in Google Docs with about 20-30 rows, with the following 6 columns : `ID, Name, Address, Cell_No, Latitude, Longitude` where the last two columns are empty. What is an easy way, for a novice user, who is a non-programmer, to geocode the address and populate the Latitude and Longitude fields with result of Geocoding? I just need to populate the Latitude & Longitude fields. I do not need to see any data on a map.


title: How to color lat/lon points in given RGB data?

text: I am using QGIS on a Mac. I am working with hydroacoustic data for mapping and interpreting the seafloor. Data are collected, processed and interpreted. The result is a txt file with Lat and Lon in the first columns followed by R, G and B data, e.g. [54.20,7.78,0.94,0.62,0.37]. The files can have 100,000 lines or more. I am not really good in GIS but I need to use it for some reasons. I seem to be unable to plot these data in GIS. There is obviously no symbology/style way to tell the system to take the color values from a file (i.e. the next 3 values in the same line of the file). It is so easy in MatLab: scatter(lat, lon, symbolsize, [R G B]) but I can't work easily with multiple layers in Matlab.


title: Is it possible to do server-side (WMS) rendering / caching of vector data

text: I have to develop a webapp that will have to display a large number of polygons on a webpage, with complex filtering options. The source data never changes (I will simply reload a new version of the DB every few months). I will probably use OpenLayers on the client-side. I developed a prototype that drawed the feature as vectors in javascript, but I'd like to switch to server-side rendering for better performance. I could probably do my app connecting directly OpenLayers to a few different WMS layers in GeoServer (postgis backend) with CQL filters, but I have the feeling GeoServer is not very flexible when using a combination of complex filters and dynamic styling of my WMS and the solution could become quite complex. From the developer point of view, a more confortable/flexible approach seems to use GeoDjango and its spatially-enabled data models to generate a dynamic (content and styling) KML file of my data. A (probably naive) elegant solution would be to integrate WMS rendering (and possibly caching) layer between the OpenLayers client and the GeoDjango app that generates KML output. The workflow would be: 1) Openlayers sends user query to "intermediate layer" 2) This layers forward query to the GeoDjango app, get a KML file and render it as WMS and return it to OpenLayers. (alternatively, OL could also talk directly to GeoDjango it sould use a middleware to postprocess its KML result before returning it). As my data almost never changes, I could then benefit from a caching tool such as GeoWebCache between OL and the "rendering layer". I know Google Maps API do server-side rendering of KML files when they are large, but it's a black box tied to their JS library, and I'd really like to do something similar with OL. I've done a few searches, but was not able to find any tool for the "intermediate layer". WMS servers such as GeoServer/Mapserver can use a KML data source, but it seems limited to a static file somewhere in the filesystem, and not from a dynamic URL. Does a tool like this exists ? Or is my approach simply wrong ? Thanks in advance, Nicolas


title: Performance of Picture Marker Vs Character Marker (font)

text: We have a point data type to symbolise. The feature classes can often contain more than 6000 points up to 20000. There are about 30 different symbol catergories to these points. My question is what's more performant - creating EMF pictures and symbolising using picture markers or create a font and use character marker symbolisation? We would be using this dataset in both ArcGIS Desktop and ArcGIS Server (silverlight API). Thanks! Rob


title: ArcPy or FME for adding history to an object?

text: I have a very specific problem and I'm looking for a solution for it. I need to add some sort of history to an object. The history will be stored like this:   * begin_date   * end_date   * previous_id Suppose the object that I'm wanting to add the history to is X. We have a central database (at the moment all mdb files, we will switch soon to gdb though) that contains the object X with the history. Workers on the field have a local version of the database and they also fill an object X_i that has an attribute changed. X and X_i have an id to track them. Now when we want to merge the local version of the database a script should run and fill the history of the central db like this:               if X_i.changed = true:         search central_X.id, fill end_date with today         search local_X.id, copy this attribute to central_X and:             fill begin_date             fill previous_id with central_X.id             generate new unique id for the central_X      We also add new fields to the central db and update the unique ids:               for every i in local_X:         if i not in central_X:             copy i to central_X             generate new unique id for i             enter begin_date      For deleted items we do something similar but then we set the end_date. At the end of the merge. We create a new X_i where all the ids are correctly filled and the changed attribute is empty everywhere. My question about this is: _I was thinking about doing this in ArcPy or FME because I have some experience with both. What do you think would be best, ArcPy or FME (or something else)._ _Can you point me in the right direction; what are the toolboxes, functions you would use for this?_ I have looked at alternatives like an ArcSDE with archiving, but that would mean we would need to invest in an ArcSDE, change a significant part of our workflow and not all items that we change will have X_i.changed, so the archiving would contain way too much information to be really useful. Edit: added workflow to clarify. ![workflow](http://i.stack.imgur.com/zPA29.png)   1. I can use changedetector for this.   2. If X_i.changed is true   3. If X_i.id is not empty   4. set end_date for all processed data   5. **no idea how to do this**   6. **no idea how to do this**   7. **no idea how to do this**   8. find the biggest id so far and proceed from this id.


title: How to clear in_memory workspace within ModelBuilder?

text: I have been using `in_memory` workspace within ArcGIS 10.1 ModelBuilder as part of a workflow to mosaic and clip raster datasets to study area bounds (Figure 1). The model runs fine the first round but crashes the program on the second or third run. I suspect two issues may be occuring: 1) Memory limitation are being exceeded or 2) There is conflict when trying to overwrite temp files `in_memory`. **Is there a way to clear`in_memory` workspace from within ModelBuilder?** Figure 1. ![enter image description here](http://i.stack.imgur.com/gJaRX.jpg) **Edit:** The final model that is now stable ![enter image description here](http://i.stack.imgur.com/zJntN.jpg)


title: Preventing symbols from different layers from interfering with each other in ArcGIS online

text: I am tackling a project that contains several new elements for me since I am a high school social studies teacher and not a full-time GIS professional. For this particular project I am trying to show the various US interventions during the Cold War around the world. I have created individual points for each intervention. These points are in a shapefile that contains a list of 49 different interventions. The fields include the country of the intervention, the year, a short description and whether or not a certain type of action occurred. Example: Guatemala, 1954, (description) and the following values for the six types of action possible: Command=1, Troops=0, Bombing=1, Naval=0, Nuke=1, Jets=0. What this means is that in 1954 we used Command Operations, Bombing and Nuclear Threats to intervene in Guatemalan affairs. On this particular point then I would like to see the three distinct symbols for those types of action. Currently I can achieve this by creating 6 layers but the symbols are all stacked on top of each other. This has all been done in Desktop 10.1 Advanced. I plan to create a feature service on our AGOL account once I have solved this problem. Thanks for any advice.


title: How to Find points in a Kilometre radius

text: I am new to PostGIS and GeoData. I need to find where the lightning has stroke by a period of time and kilometre radius based on a given location. I have given lon/lat and my table is:               latitude numeric not null,     longitude numeric not null,     time timestamp(0),     geom(point,4326)      I do not know how to use within or radius query. I basically need to find how many points (lightning i have in a radius of a given location(lon/lat). I am using PostGreSQL 9.2 and PostGIS 2.0. Thank you.


title: Can GeoEXT 1.0 be used along with ExtJS 4.0?

text: On the GeoExt website their tutorial page has some instructions on how to setup up a application to use GeoEXT. This is the page I am referring to. It mentions that you should download GeoEXT 1.0 and ExtJS 3.2. I have since noticed that Sencha has released ExtJS 4. Can ExtJS 4 be used along side GeoEXT 1.0 without any issues? Or should I go by the book and install ExtJS 3.2? Thanks Ando


title: How do you Convert a .sp1 file into a Shapefile

text: I was previously able to convert .sp1 files into shapefiles, but that was last February. I can only remember how to do some of it. How do you convert them into shapefiles ? These are preplot receiver / source points as I am using ArcView 9.3


title: Finding optimal network with python and network analyst in arcgis 10.0

text: Hello GIS StackExchange, I'd like to know if there's any way to find the shortest distance across the diagonals of a polygon. I am creating a network analyst solution to determine the most cost effective route between 2 points. The solver will travel along the roads edge for the most part to determine optimal routes but there may be situations where cutting across certain polygon features will be more efficient than continuing along the edge of the road. The question here is what is the best way to cut across the traversable polygon features? ![Sample route optimizer](http://i.stack.imgur.com/gkOK6.png) The image is a sketch i did to illustrate how i'd like the route optimizer to function. The blue point feature is the route origin point while the green feature is the destination point. The green polygons are parcels that can be traversed, while the yellow polygons cannot be cut across. The blue line is the road network which is the main edge feature to determine the optimal route, the purple diagonal(ish) line is the result of the router determining that cutting across the polygon would indeed be more optimal than following the road. The idea i have for cutting across the polygon features is to first determine the shortest distance within all traversable polygons and then create a line feature class containing the shortest distances for all the traversable polygons in the network dataset. This line feature class will then be connected to the same connectivity group as the roads feature class to ensure the router can connect both features and possibly determine that cutting across the polygon is indeed more optimal than following the road edge in certain situations. One problem i have encountered in my research is that polygons may be made up of countless number of vertices and as such it would be difficult to determine the shortest path across the polygon. This question is related to my thesis project, it is to be written in python within arcgis 10.0. I have done quite a lot of research already so don't feel obliged to remind me about the need to learn this for myself because i have and i'm stuck. My biggest fear is designing the model half way only to find out that the approach i took was wrong from the start. That's why i need some help. I'd appreciate any documentation that can help or any complete resource that will prevent me re- inventing the wheel, or even ideas about how to go about certain aspects of building the router. Thanks as always gis stackexchange. Apologies for the lengthy post.


title: How do I implement a custom IFeatureDatasetExtension?

text: **How can I implement IFeatureDatasetExtension and register it so that ArcGIS discovers and initializes it?** In arcobjects, I see that Esri has a number of coclasses that implement IFeatureDatasetExtension2. Using Categories.exe, I see there is a COM category called "Esri Workspace Dataset Extensions". ![enter image description here](http://i.stack.imgur.com/QizOF.png) I don't see anything on the FeatureDataset coclass that would let me register extensions, the closest thing is IFeatureDatasetExtensionContainer, but it only lets me access extensions that are already defined. There is no `Add` method.


title: Trying to use the GeoWebCache "gmaps" service

text: I am using GeoServer 2.4.5 with built-in GeoWebCache 1.5.0. I'm integrating GeoServer with a Google Maps V3 api frontend, and I got it to work reasonably well using the WMS service. However, I noticed that there's a builtin 'gmaps' service, as described here: http://geowebcache.org/docs/1.5.0/services/gmaps.html I've tried to use it, with a url like:               http://myserver:8080/geoserver/ok/gwc/service/gmaps?layers=ok:OK_Twp&zoom=9&x=114&y=200&format=image/png      however, I get the following response:               This XML file does not appear to have any style information associated with it. The document tree is shown below.     <ows:ExceptionReport xmlns:ows="http://www.opengis.net/ows" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" version="1.0.0"     xsi:schemaLocation="http://www.opengis.net/ows http://myserver:8080/geoserver/schemas/ows/1.0.0/owsExceptionReport.xsd">     <ows:Exception exceptionCode="InvalidParameterValue" locator="service">     <ows:ExceptionText>No service: ( gmaps )</ows:ExceptionText>     </ows:Exception>     </ows:ExceptionReport>      For some reason, it can't find the 'gmaps' service. Any ideas? I've made sure to enable "Global Services", but can't find any configuration regarding the gmap service in the GeoServer web interfact.


title: How to offset/displace overlapping lines in a web map?

text: I have a dataset of points representing a To/From status - something like distribution centres and delivery addresses. When the user clicks on a point I'd like to show the relationships between this point and its related points using lines. Eg when clicking on an address, show the path to the delivery centre. This works fine, except when there are multiple overlapping lines. In this situation I'm looking for a method to separate the lines automatically, so that it's clear there are multiple lines, and each can be interrogated. I've put a simple example on JS Fiddle - this sample creates 5 overlapping lines. Is there any way to displace them automatically, while retaining the start and end nodes? ![enter image description here](http://i.stack.imgur.com/emU4b.png) This is the effect I'm aiming for, ideally handling multiple overlapping lines automatically: ![enter image description here](http://i.stack.imgur.com/Pro8R.png) This example uses Google Maps, but an answer in ArcGIS Server's JS API would be acceptable too.


title: How to get polyline out of multipart feature using ArcObjects C#

text: I've got a polyline feature class and I've identified the first feature that is a multipart feature. I now want to find out which of the paths of this multipart feature intersects with a point. Up till now I managed to fill a geometry collection with the parts (paths) of the polyline. My method to check for the intersection requires a polyline input. And this is where I'm stuck. How can I convert/cast/? one of the two paths in the geometry collection to a polyline? Here is some code, but `pPolylineCSPart` is null:               //right polyline into geometry collection     IGeometryCollection pGeometryCollection = new PolylineClass();     pGeometryCollection = pRightPolyline as IGeometryCollection;     //loop through geometry collection     for (int intPartCount = 0; intPartCount < pGeometryCollection.GeometryCount; intPartCount++)     {         IPolyline pPolylineCSPart = new PolylineClass();         //cast polyline part to polyline         pPolylineCSPart = pGeometryCollection.get_Geometry(intPartCount) as IPolyline;         //check if part intersects river         bMultipartIntersectRiver = PolylineIntersectioCheck.Are2PolylinesIntersecting(pPoint, pPolylineCSPart);         if (bMultipartIntersectRiver)         {             pSplitFeat.Shape = pPolylineCSPart;             SplitCSFClass.CopyAttributes(pTransectSplit, pSplitFeat);             pSplitFeat.Store();         }     }      Many Thanks!


title: What is gin_trgm_ops of postGIS

text: I am using a query where the following line written               CREATE INDEX trgm_idx ON public.roads_r USING gin (name gin_trgm_ops);      Its returning the following error > operator class "gin_trgm_ops" does not exist for access method "gin" How can I overcome this?


title: Tracking field personnal using android app and showing in ArcGISonline

text: Planning to equip field personnel with GPS enabled phones (android) and collect the location and photograph of assets. Can we use ArcGIS for android in the phone, use this to upload GPS and photos , and the use ArcGIS online to show these locations ? Do we need ArcGIS for server or simple we can do it using ArcGISOnline and android SDK ? Adarsh


title: PostGIS ERROR: function x(geometry) does not exist

text: I am using the following postgis query for updating a pgrouting table from my custom road data. In the `x(ST_StartPoint(part_geom)) as x1,` line its returning an error mentioned > ERROR: function x(geometry) does not exist               INSERT INTO public.roads_r (     SELECT DISTINCT ON (the_geom)     nextval('public.roads_r_gid_seq') AS gid,     gid AS parent_gid,     name,     no,     type,     NULL as source,     NULL as target,     ST_Length(part_geom, true) as cost,     ST_Length(part_geom, true) as reverse_cost,     x(ST_StartPoint(part_geom)) as x1,     y(ST_StartPoint(part_geom)) as y1,     x(ST_EndPoint(part_geom)) as x2,     y(ST_EndPoint(part_geom)) as y2,     NULL as to_cost,     NULL as rule,     ST_Multi(part_geom) AS the_geom     FROM (         SELECT              a.gid,              a.geom,              (ST_Dump(ST_Union(b.geom))).geom AS part_geom,              a.road_name as name,             a.road_no as no,             a.road_type as type         FROM public.roads as a, public.roads as b         WHERE ST_Intersects(a.geom, b.geom)         GROUP BY a.gid, a.geom, a.road_name, a.road_no, a.road_type     ) AS res     WHERE ST_Contains(geom, part_geom) = true OR ST_Contains(ST_Buffer(geom, 0.00001), part_geom) = true)      Any help will be highly appreciated.


title: Install shapely for QGIS on Mac

text: I know there are a number of similar questions, but none of them have answers that allow me to use shapely in QGIS. In particular I am trying to install the Contour plugin, but this fails with the message "The plugin depends on some components missing on your system. You need to install the following Python module in order to enable it: shapely.geometry." The reason for this error is explained well here: How to install the Python module shapely for QGIS use? These are the steps I have taken to solve the issue: I first installed shapely via pip. This allowed me to `import shapely` in python 2.7, which is the default on my system, but installing the plug in still failed with the same error. This question How to install shapely for python 2.6 (Mac) NOT 2.7? explained that qgis uses python 2.6, not 2.7, so I downloaded the tar file and did a ` sudo python2.6 setup.py install`. Now I can `import shapely` in python 2.6 too. But the installation of the plugin still fails with the same warning! I've also tried `import shapely` from the python console within QGIS, and there is "no module named shapely." Is there something I've overlooked? EDIT. Thanks to the comments from gene below, it turns out that QGIS on Lion does use python 2.7, not 2.6. But still it doesn't work. It turns out though that I have a few different versions of 2.7: From the terminal:               $ python     Python 2.7.2 (default, Mar 8 2012, 15:29:59)     [GCC 4.2.1 (Based on Apple Inc. build 5658) (LLVM build 2336.1.00)] on darwin     Type "help", "copyright", "credits" or "license" for more information.     >>>        From the QGIS python console:               >>> import sys     >>> print sys.version     2.7.1 (r271:86832, Jun 16 2011, 16:59:05)      [GCC 4.2.1 (Based on Apple Inc. build 5658) (LLVM build 2335.15.00)]     >>>       From the terminal again:               $ python2.7     Python 2.7.3 (default, Oct 22 2012, 06:12:32)      [GCC 4.2.1 Compatible Apple Clang 3.1 (tags/Apple/clang-318.0.58)] on darwin     Type "help", "copyright", "credits" or "license" for more information.     >>>       So it seems that QGIS uses 2.7.1, Apple defaults to 2.7.2, and I also have 2.7.3. So I need to install shapely for 2.7.1?


title: XML workspace document import missing data

text: Ok, building on my last post on this topic, export/import from/to SDE. I export my SDE using export "XML Workspace Document..." and I'm using the Z File format because of the large size of the database (exported Z files still approaching 10gb in size). Export process runs successful now, and the import process into a new SDE runs successfully but, when I start checking through all the feature classes, some are missing their data. Schema comes across fine though. Can seem to find any pattern with the feature classes missing data. Truly puzzling. Any help appreciated Thanks


title: How to calculate volume of bund around pit from TIN and external feature edge

text: I need to do some volume calculations of bunds. I suppose I should use the Surface Volume tool in 3D Analyst of ArcGIS. Here is the tool description. I am uncertain about the plane height parameter, however. From the tool webpage there is a diagram showing a plane height that intersects the surface. This seems to match our circumstance, where the bund rises above the polygon describing the external edge of the feature, but then dips below the plane of this polygon on the internal face of the bund going down to the base of the pit. Here is the complicated part. The elevation of the external edge of the feature differs. Since, therefore, it is not a perfectly horizontal plane, it doesn’t seem sensible to me to put in an elevation value for the plane height. If a constant value were used, then the bund may not be cut at its external base where the elevation on the boundary varied above or below that value. I see there is also the cut/fill tool. To me this would require a raster of the polygon describing the outer, bottom edge of the bund to be supplied. I'm not sure how to produce a valid raster representing this plane, however. Thoughts?


title: What are the essential elements of a cartouche, if any...?

text: I'm trying to design artistic cartouches for some maps that show traditional knowledge locations from specific First Nations in Canada so the cartouche's will reflect their Nation's identity, culture, etc... I'm not having trouble with ornamental features or colours. I struggle trying to fit more modern elements like titles, descriptions, data sources, and maybe even projection/datum information. **Are there some guidelines I should follow when trying to decide what other elements "must" be in the cartouche?** In this case, the audience will be the knowledge holder (person who contributed to the information on the map) so that they can have a hardcopy of their work to share with whom they please. The second set of maps will be for the community showing the comprehensiveness (everyones information) of the study. Possibly displayed on a wall in the community building for example. _Below is a small image of a **simple** cartouche from a historical map I have for those who don't know what I mean:_ ![Cartouche on Klondike Gold Rush Map showing the All Canadian Route](http://i.stack.imgur.com/YUe5U.jpg)


title: How do you add file name to an attribute field?

text: Does anybody know how to add the shape file name to a field in the file's attribute table? I've found a description: > 'You have to use inline variable in CalculateField like %Name% in > expresssion so that whatever the name of the file is, the string will be > inserted based on Name of the dbase file !' I've tried to CalculateField tool with %Name%, but I got message:               Executing: CalculateField "2011-11-6 15_41_point" DBFName %Name% PYTHON #     Start Time: Thu Dec 08 23:34:23 2011     WARNING 000405: No records within table     Succeeded at Thu Dec 08 23:34:24 2011 (Elapsed Time: 1.00 seconds)      Where is my mistake?


title: Datatypes to return in GeoServer WPS

text: Is it possible, that my created WPS (in Geoserver) returns an image or a pdf (general: returns a datatype, that isn't a supported format by Geoserver). Is it sufficient to expand the ProcessParameterIO class? Can anybody give me a general overview, how to do that?


title: Advice on choosing ideal software for consultancy

text: I am taking the step to start a small business consultancy using digital map technology to assist retail business to better understand their environment and make informed location based decision. My focus is retail business and service planning, however, my challenge is the choice of software. ESRI business analyst is not available to the part of the world where i live due to inadequacy of data but i don't want this to deter me from starting the business. I know there is always a way around obstacles if you work at it long enough. While i work on solving the data adequacy issue i want to know if map-info professional will be a good software to purchase (I've never used it before). Any advice?


title: Intersection of convex polygons on a spherical surface

text: There are a number of algorithms available for the intersection of two convex polygons, but I am interested in an algorithm to find the intersection of two convex polygons on the surface of a sphere, where the edges of the polygons are great circles. The polygons might be arbitrarily large (i.e. no approximation for small areas). Is there a description of such an algorithm available? Better, is there any Python or C code out there that would do this, and has a liberal open source license or is public domain?


title: How to treat missing values as null

text: I have an excel file with missing values encoded as blanks and I want this excel file to be converted into DBF but with the conversion process, the blanks (missing values) become 0. The reason I am doing this is because I want to perform interpolation but because of the missing values becoming 0 the interpolation becomes wrong. My question is how to treat missing values as without deleting the point shapefile and still proceed with interpolation.


title: Is there any random data generator supporting geometry in PostGIS?

text: I found few random data generators for PostgreSQL (such as this one), but none supporting PostGIS geometry type, other than point. Making points would be quite easy, but for lines with multiple nodes and for polygons (I need more complex polygons than simple rectangles) it would be hell to write the function myself. I would prefer using others' tools (free and opensource, if possible) instead of making them myself. EDIT: there are algorithms generating random points within a polygon. For polygons, I would prefer a different approach: to generate quasi-centroids first and then make polygons around them. I have some idea how to do it (random starting direction or random angle from previous node + random distance from quasi-centroid). Although if I found a tool generating polygons from scratch, making a centroid or a random point within would be easy. If no one points to any such a tool in a few days, I'll make it myself and post the code as an answer. I don't need complex line generator so much now, so I'll either make it later, or wait for some other answers. Working complex line generators will get at least +1 from me even if I found I won't really need them, and if combined by some good complex polygon generation they should merit acceptation.


title: Expression based labeling: Look for points in defined radius/buffer

text: How can I do that with qgis expression based labeling: CASE WHEN _within a radius/buffer of 1km (outgoing from every point with attribute 'place') is a point with attribute 'museum',_ THEN 'M' ELSE WHEN _within a radius/buffer of 1km is a point with attribute 'historic'_ THEN 'H' ELSE NULL END


title: Is there an opensource way of making 'streetview' panoramas?

text: I'm interested in making _streetview_ panoramas and am curious if there are any opensource approaches that I should consider for both making and publishing these panoramas. I've been playing with the Android functionality on my phone to make Android Photospheres (using a camera app SDK - here for the curious). These can be published in a streetview-like panorama via Google. This works, and is nifty but it requires Google to approval the panorama for it to be published online. Another approach is to generate my own tiled grid of photos and use the google stretview api: https://developers.google.com/maps/documentation/javascript/streetview#StreetViewOverlays I'd much rather produce these using some opensource tools, so I'm not completely locked into Google. In this interview with Steve Coast from OpenStreetMap, he refers to a device you can put on your phone which will result in _streetview_ like imagery. Is there a recommended way for this imagery be collected, stored, stitched together and shared?


title: How to synchronize the marker's rotation field with the angular distortion in the map projection

text: I'm Using QGIS 1.8.0 to plot a global dataset showing the positions of aircraft in flight. The data I have are .csv files with three columns: Latitude, Longitude, and angle of travel. I assign the "angle of travel" column of data to the rotation field in the "Style" tab of the "Layer Properties" dialog window. This works fine except that that the rotation angle does not take into account the angular distortion of the map projection. This means that in places where the angular distortion is high, the flight path looks like a row of dominoes: ![Flight paths in an area with high angular distortion](http://i.stack.imgur.com/grgkb.jpg) Whereas where the angular distortion is low the flight paths looks like a straight line, which is what I want. So my question is: **In QGIS is there a way to synchronise the marker's rotation field** (now called "angle" in QGIS 2.0.1) **with the angular distortion in the map projection?** _EDIT:_ Since posting this question I have started using the new QGIS 2.0.1.   In QGIS 1.8.0 I dont think it was possible to do what I am asking, but in the new version it looks like it is. However, it requires a thorough understanding of the maths behind map projections, which I do not have. In the "Style" tab of the "Layer Properties" dialogue window, when assigning "data defined properties" to the angle of the marker, it is now possible to open an "Expression string builder" dialogue window. ![enter image description here](http://i.stack.imgur.com/ShzE0.jpg) Here, I presume, it will be possible to enter a mathematical equation which modifies the angle of the marker so that it accounts for the angular distortion in any particular map projection. The map projections I am interested in are **mollweide** (or other similar elliptical projection) for a global overview, and also the **UTM zones**. The reasons for my map projection choices are: Mollweide:   I like the look of mollweide for showing global trends. UTM:   To make a complete visualization using all the data and revealing some of the lesser used routes clearly I will have to import and rasterize the data in geographic portions as I don't have enough RAM to deal with the entire dataset in one go.   I was thinking of using UTM zone strips, a bit like a set of globe gore, which I can then stitch back together once rasterized using Geocart. The mercator projection idea below is not ideal as it creates a pronounced narrowing of features near the poles, and there is a fair amount of air traffic in the north.   The conformal projection idea represents a lot more work for me, as many more tiles will be needed than if I use pole to pole UTM strips. So, Im still holding out for someone to tell me at least whether or not it's theoretically possible to type in a formula that will affect all the points in the correct way, if not the formulae themselves. Many thanks


title: Which Landsat version is appropriate to imply a supervised classification with GRASS GIS?

text: I need to imply a supervised classification with GRASS GIS to find build-up areas. I want to avoid Landsat 7 images due to SLC problem.   1. Can I use images of Landsat 8 or it's better to use images of Landsat 5?    2. Are images of Landsat 8 full operational or are they in calibration/test mode?   3. Which version is best supported from Grass GIS (6.4 or 7)?


title: How to get a list of attributes in FME when there are no records in table?

text: How to get a list of attributes in FME when there are no records in the table? I try use a `PythonCaller` but when there are no records, FME passes over this table. In `PythonCaller` I write:               import fmeobjects     def FeatureProcessor(feature):         print feature.getAllAttributeNames()      Is it possible to get attributes when there is no record?


title: Dialogue window for ground to grid tool in arcmap 10

text: I am using ground to grid correction tool with some success but in the window I cannot see all the print. Is there a way to size the window? Does anyone know how to apply ground to grid correction ? For this is self taught and I see the results but not sure if the workflow I am using is the best approach. I am using Arcmap 10 on a network.


title: Does ArcMap 10.2.1 no longer listen for DELETE key?

text: We've noticed that in ArcMap 9.3 we could select a feature on the map (while editing of course) and just press the delete key on the keyboard to delete the feature. ArcMap 10.2.1 does nothing when you press delete. If we click the delete button (the X button on the toolbar) it will delete just fine. Did ESRI take away the keyboard event listener? Is there an option we aren't seeing?


title: How can I access the AddShapefile method in HookHelper?

text: I'm new to ArcEngine and ArcObjects. I'd like to develop an extension for ArcMap where I have to use `IHookHelper` instead of `AxMapControl`. I am looking for the `AddShapefile` method, which is available for `AxMapControl`, but apparently missing from the `IHookHelper` interface.               public partial class conversionRtoP : Form     {         public IHookHelper phookhelperRtoP;         public IRaster praster;              public conversionRtoP()         {             InitializeComponent();         }              private void conversionRtoP_Load(object sender, EventArgs e)         {             for (int i = 0; i < phookhelperRtoP.FocusMap.LayerCount; i++)             {                 if (phookhelperRtoP.FocusMap.get_Layer(i) is IRasterLayer)                 {                     comboBox1.Items.Add(phookhelperRtoP.FocusMap.get_Layer(i).Name.ToString());                 }             }         }              private void button1_Click(object sender, EventArgs e)         {             Geoprocessor gp = new Geoprocessor();             IConversionOp pconversionop;             pconversionop = new RasterConversionOpClass();                  for (int i = 0; i < phookhelperRtoP.FocusMap.LayerCount; i++)             {                 if (phookhelperRtoP.FocusMap.get_Layer(i).Name == comboBox1.Text)                 {                     praster = (IRaster)phookhelperRtoP.FocusMap.get_Layer(i);                 }             }                  ESRI.ArcGIS.ConversionTools.RasterToPolygon prastertopoly = new RasterToPolygon(praster, outputpathtext);             ESRI.ArcGIS.Geoprocessing.IGeoProcessorResult gpresult;             gpresult = (IGeoProcessorResult)gp.Execute(prastertopoly, null);             int j;             string gppath;             gppath = outputpathtext.Text;             j = gppath.LastIndexOf("\b");             // this is my problem:             ***phookhelperRtoP.(where is addshapefile method)***         }              private void button3_Click(object sender, EventArgs e)         {             SaveFileDialog savedialog = new SaveFileDialog();             savedialog.CheckFileExists = true;             savedialog.Filter = "Shapefile (*.shp)|*.shp";             savedialog.Title = " Save";             savedialog.RestoreDirectory = true;             savedialog.FileName = comboBox1.Text;             DialogResult dr = savedialog.ShowDialog();             if (dr == System.Windows.Forms.DialogResult.OK)             {                 outputpathtext.Text = savedialog.FileName;             }


title: Geometry Conversion Problem

text: How do I convert from an ESRI.ArcGIS.Geometry.IPoint to an ESRI.ArcGIS.Client.Graphic point object? In other words, I'm wanting to convert from an IGeometry to a runtime geometry, I guess.


title: Simplify polygon in OGR

text: We use ArcObjects pretty heavily and simplify polygons using something similar to the `ITopologicalOperation::Simplify(void)`. No tolerance is passed in as the parameter (it uses the XY tolerance of the associated spatial reference). Now, we're using OGR for some functionality and would like to simplify using this library. What's the equivalent? I'm aware of `OGRGeometry::SimplifyPreserveTopology(double dTolerance)` but I'm not sure where to get the tolerance from the in OGR library.


title: Add TileCache Schema Layer to Leaflet Map

text: I am trying to get Leaflet to display a TileCache Schema tileset. One example: http://example.com/tiles/Admin/06/000/000/357/000/000/158.png As far as I can tell, Leaflet only supports TMS style tiling schemas such as: http://a.tile.openstreetmap.org/15/16938/11166.png I tried to find an extension of a L.tileLayer that supports the TileCache tiling schema, but no such luck.


title: Problem with SPIT (error while importing layers to POSTGIS)

text: I was attempting to import a polyline roads layer into postgis, using both SPIT in QGIS 1.7 and postGIS GUI for this purpose. I have received errors through both processes as described here:   * SPIT > The database gave an error while executing this SQL: INSERT INTO > "public"."highways"("STREET","FROMLEFT","TOLEFT","FROMRIGHT","TORIGHT","STREETCODE","OBJID","USERID","FJUNCTION","TJUNCTION","LENGTH","ROA... > (rest of SQL trimmed) The error was: ERROR: value too long for type > character varying(31)   * postGIS GUI ![Error msg. from postGIS](http://i.stack.imgur.com/TZAN6.png) What can be the problem and how can I fix it? Thanks for your response. Dor


title: I need QGIS routing plug-in

text: I am new to pgrouting. How does it work? What benefits can I get from it? I tried the examples, read some articles, but could not find the QGIS routing plugin for Windows. It is not showing up in the plug-in manager. Finding a route in a given dataset by mentioning a start and end node works, but it is not performing on my local road network shapefile. The query executes, but doesn't show any results. Also it says > Postgres with PostGIS and pgAdmin QGIS with PostGIS Manager and RT Sql Layer > plugins I don't find PostGIS Manager and RT Sql Layer plugins. Have seen osm2po articles, how do they work? What do I do next?


title: Error while trying to view the attribute table

text: I have a shapefile which displays fine in ArcMap. When I view the attribute table I get the following error. However the same shapefile loads fine and can view its attribute table to in QGIS. Any ideas ?? ![enter image description here](http://i.stack.imgur.com/ziCSM.jpg)


title: Downloading Geodatabase to ArcCatalog

text: I think this is a fairly easy question so apologies. Someone sent me a .gdb in Google Drive. When I download it to my computer, it is a file folder with all the .gdb contents in it. I can't figure out how to add it as a .gdb in ArcCatalog - I know it must be straightforward but I'd be very grateful if someone could explain it to me. Thanks in advance!


title: Examples of GIS Visualisation

text: **Examples of GIS Visualisation** [This is a community wiki post, **you are inviting other members of the community to edit it** - not just to correct spelling or improve formatting] Example: Dr. John Snow (1813-1858), a legendary figure in the history of public health, epidemiology and anesthesiology created one of the first uses of geographical analysis on a static map to solve a problem of stopping the spread of cholera through Soho, London, England. ![enter image description here](http://i.stack.imgur.com/L8vfX.jpg) http://upload.wikimedia.org/wikipedia/commons/thumb/2/27/Snow-cholera- map-1.jpg/643px-Snow-cholera-map-1.jpg The first two cholera cases in 1848 occurred in the same room but a week apart at No. 8, New Lane (in red) near Gainsford Street in the community of Horsleydown... water pumps were causing the spread and over a period of time it was found after collating deaths and visualising them. **Good Examples can be ancient, historical or cutting edge.** ![enter image description here](http://i.stack.imgur.com/F2uLN.jpg) source: http://www.ph.ucla.edu/epi/snow/1859map/snowindex.html


title: Field Calculator decimal places - how can I limit them?

text: I've added a new column to an attribute table and performed an operation with the field calculator to populate the new column. The result is a column with a percentage. My problem isn't major but I can't seem to limit the number of decimal places so that I just have solid numbers. Right now I have 13 places after the decimal, I would like it to be 0.


title: Why are the formats of these coordinates different?

text: I'm really sorry if this is a silly question (and I'm sure it is) but I have tried searching but I am a complete newbie to coordinates and maps and I don't know what I'm looking for so am not finding my answer! I have a gpx file which has coordinates like as follows: > 51 47.322,-2 47.218 > > 51 47.340,-2 47.071 I need to add some new coordinates but I am taking these from google maps,but the coordinates are in a different format. ie: > 51.804979,-2.738793 I know there is a logical reason for this as I've read about different coordinate types, but I don't understand why and what types they are and what I need to do to convert from one to the other? Thanks for any help! Bex


title: 'Select' and 'hover over' features are not working although no error is displayed (openlayers)

text:
On the following html code i get no errors when i check it on firebug. When i select or hover over features i want to see them "highlighted". When i move my cursor onto a feature i can see that something is loading but it is not highlighted. At the same time firebug shows that "post wms" has worked fine and i can see that it retrieves the total size of the feature (Kbs). I can not figure out what the problem is... I am using Openlayers and Geoserver for displaying the shapefiles Code:                                           <title>Welcome to Choros</title>         <link rel="stylesheet" href="http://openlayers.org/api/theme/default/style.css" type="text/css" />         <link rel="stylesheet" href="http://openlayers.org/dev/examples/style.css" type="text/css" />              </style>
Chώros

                 coordinates

             Click on the map to get feature info
           * * * Hello Geographica and sorry for the late reply. When i step through the code line by line i can see that firebug skips the following 2 lines: hover.addFeatures([e.feature]); hover.removeFeatures([e.feature]); As you mentioned before i think that e.feature is not valid. Attached is the image that shows the response that i get from firebug. Can you actually figure out why the e.feature is not valid. I can not see any reason why it shouldnt be. ![firebug Response](http://i.stack.imgur.com/cR83t.jpg)


title: Registering a SDE Service in 10.1 with Command Line

text: I am installing ArcSDE on a new development box. I have executed the setup.exe, created an Enterprise database called "GDB123" and created a new sde service called "sde_123". I need some clarification on how to register the service using the administrative commands. The install guide says if your are using dbo-schema register the service with this command `sdeservice -o register -d SQLSERVER,MYSERVER\SSINSTANCE -r SDE_DBA_USER -v sa -i sdesvc` I do not understand what value I am suppose to type in for the `-r` and `-v` options. I think `-r` needs to be the database name "GDB123" but I have no idea what to put for the `-v` option. Heres what I have so far `sdeservice -o register -d SQLSERVER -r ?? -v ?? -i sde_123` How do I make this work? ArcSDE 10.1 install for SQL Server on Windows 2008 RS box Install Guid: http://resources.arcgis.com/en/help/install-guides/arcsde- application-server-microsoft-sql- server/10.1/index.html#/Creating_an_ArcSDE_service_on_Windows_using_the_sdeservice_command/00950000003v000000/ SDE Service Commands Options: http://resources.arcgis.com/en/help/install- guides/arcsde-application-server-microsoft-sql- server/10.1/index.html#/The_sdeservice_command/00950000004m000000/


title: How to locate events at intersections using linear referencing?

text: I am having problem in locating intersection crashes. Actually I have intersecting legs with their specified begin and end mileposts. I am trying to locate crashes based on the milepost values. Everywhere it seems fine using arc map 10.1 linear referencing methods except the places where crashes lie exactly in between the intersection legs i.e. exactly on the intersecting point of segments. The crashes are giving milepost based on different leg segment rather than the actual leg. For example in a intersection of street A and B, a crash lies exactly at the intersecting point. The attributes of that crash says it occurred in street A but the linear referencing gives the milepost value based on Street B. Please share your knowledge or ideas if anybody has solution to this problem. Many Thanks!


title: Stereo Compilation Software with PostGIS

text: Does anyone know of stereo compilation software that will allow for the writing into a PostGIS database? Both Datum and Stereo Analyst will write to an ESRI Geodatabase. Standard stereo compilation software will write to a DGN, but I am unaware of anything that will write directly into a PostGIS database. Thanks...


title: Skimage graph route to multiple destinations

text: I want to route through an array with `skimage.graph.route_through_array`. I have one start point and multiple destination points. I am wondering if it is possible to insert the multiple destination into the function. I know I could also loop over the function, but I am looking for something faster.


title: How to create/export multiple maps in one layout using ArcMap?

text: I have an air pollution data that consists monthly mean values of 12 monitoring stations. I made interpolation and want to export the 12 different map in one layout with one legend. Since I am a newbie with ArcMap, I wonder if it is possible. I am using ArcGIS 10.1 for Desktop.


title: Multispectral image segmentation for natural-resources applications using R

text: R has the capacity for image segmentation, although all of the examples I have come across utilize a single band for the segmentation (example). I am interested in combining R's capacity for random forest image classification with an object-oriented segmentation approach. What functionality does R have for multispectral image segmentation that is suitable for natural resources-based analysis?


title: Update Existing Feature Class by Replacing Features with Updated Versions from a Shapefile

text: I have a feature class in a versioned ArcSDE geodatabase with archiving enabled. I receive regular updates for this data in shapefile format. I want to be able to replace any features that have changed with the new version from the shapefile. How can I do that without recreating the feature class and thus losing my archived history? I have gotten as far as identifying which features have changed, been added, or been deleted. But how to I select all of the features in the existing feature class that have the same ID as the changed/deleted features in the updated data and replace them with the updated features?


title: PyQGIS esri ascii grid to .bil conversion

text: Esri ASCII GRID can be converted to esri .bil format by following way in a desktop From QGIS desktop menu-->Raster-->conversion-->Translate(convert format) How can I write a python code because I have 864 esri ascii grid data?


title: QGIS Sextante Python UnicodeEncodeError

text: I have some issues with Quantum Gis 1.8. When I try to use shapes to grid (sextante-saga tool) I get the following error (comes up when I try to define the new grid file name): Fehler bei der Ausführung von Python-Code: Traceback (most recent call last): File "C:/Users/Max/.qgis//python/plugins\sextante\gui\OutputSelectionPanel.py", line 94, in saveToFile settings.setValue("/SextanteQGIS/LastOutputPath", os.path.dirname(str(filename))) UnicodeEncodeError: 'ascii' codec can't encode character u'\xfc' in position 19: ordinal not in range(128) Python-Version: 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)] QGIS-Version: 1.8.0-Lisboa Lisboa, f440b60 Python-Pfad: ['C:/Users/Max/.qgis//python/plugins\sextante_taudem', 'C:/Users/Max/.qgis//python/plugins\sextantelwgeomprovider', 'C:/Users/Max/.qgis//python/plugins\sextanteexampleprovider', 'C:/Users/Max/.qgis//python/plugins\sextante', 'C:/Users/Max/.qgis//python/plugins\permaclim', 'C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python', 'C:/Users/Max/.qgis//python', 'C:/Users/Max/.qgis//python/plugins', 'C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python/plugins', 'C:\PROGRA~2\Quantum GIS Lisboa\bin\python27.zip', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\DLLs', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\plat-win', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\lib-tk', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\qgis\bin', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\PIL', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32\lib', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\Pythonwin', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\wx-2.8-msw- unicode', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\qgis\python\plugins\fTools\tools', 'C:/Tempordner_GIS/J\xc3\xbchnde']


title: Spatial Join problems in ArcGIS 10.1 - intersect polygon and lines

text: 1. I have hurricane tracks (Lines) and Florida county boundary (polygon) and I want to extract only those hurricane tracks that intersect Florida.    2. In ArcMap 10.0 I used Spatial Join (from Analysis/Overlay) to join the hurricane tracks and Florida polygon map. This created a new shp file of only those tracks that intersected Florida. Importantly it kept the entire length of each of the hurricane track lines that intersected Florida.   3. I have tried repeating this in ArcMap 10.1 using Spatial Join and while it seems to intersect the hurricane tracks with Florida okay - the new shp file does not have the entire length for each of the hurricane tracks - only short bits of lines outside the Florida polygon boundary. Can anyone provide guidance on why this may be the case - it is critical for this exercise to extract the those hurricane tracks intersecting Florida polygon but we need the entire length of each of the tracks for each of these hurricanes.


title: Four band raster - ColorInterp of the IR band

text: I have a four band aerial photography raster. The bands contain Red, Green, Blue, and IR data. Using GDALInfo we get something like this currently:               Band 1 Block=2933x1 Type=UInt16, ColorInterp=Undefined     Band 2 Block=2933x1 Type=UInt16, ColorInterp=Undefined     Band 3 Block=2933x1 Type=UInt16, ColorInterp=Undefined     Band 4 Block=2933x1 Type=UInt16, ColorInterp=Undefined      The first three bands can be set to have a colour Interpretation of RGB respectively (GIS' get flaky if they're not set it seems). But what should the **ColorInterp** be for the fourth, InfraRed band? Left `Undefined`? Or is there a setting for IR? Thanks.