Drop down menu not populating for ArcGIS JavaScript API Application

I am using this template http://localgovtemplates2.esri.com/publicinfocenter/default.htm When editing the config.js which is in JSON format, the parameters are `{ Key: "Key", LayerURL: "http:"layerurl", OutFields: "*", RequestTypeFieldName: "FieldName", CommentsLayerURL: "url", CommentsOutFields: "*" },` Does any additional configuration need to take place in the HTML file? The drop down is not js. What might key indicate here? EDIT: Sample below is from the Local Government template of the serviceRequest.js from the template.                var singleClickFlag = true;          //function to add service request layer on map     function AddServiceRequestLayerOnMap() {         var serviceRequestLayer = new esri.layers.FeatureLayer(serviceRequestLayerInfo.LayerURL, {             mode: esri.layers.FeatureLayer.MODE_SNAPSHOT,             outFields: [serviceRequestLayerInfo.OutFields],             id: serviceRequestLayerInfo.Key,             displayOnPan: false,             visible: false         });         map.addLayer(serviceRequestLayer);              if (defaultServiceTab == 1) {             ShowLoadingMessage("Loading...");         }         var serviceRequestLayerHandle = dojo.connect(serviceRequestLayer, "onUpdateEnd", function (features) {             AddServiceLegendItem(this);             PopulateRequestTypes(serviceRequestLayer.fields);             HideLoadingMessage();         });              dojo.connect(serviceRequestLayer, "onClick", function (evt) {             if (singleClickFlag) {                 ShowServiceRequestDetails(evt.mapPoint, evt.graphic.attributes);                 singleClickFlag = false;             }             setTimeout(function () {                 singleClickFlag = true;             }, 1000);             //For cancelling event propagation             evt = (evt) ? evt : event;             evt.cancelBubble = true;             if (evt.stopPropagation) {                 evt.stopPropagation();             }         });              dojo.connect(serviceRequestLayer, "onMouseOver", function (evt) {             map.setMapCursor('pointer');         });              dojo.connect(serviceRequestLayer, "onMouseOut", function (evt) {             map.setMapCursor('crosshair');         });              var serviceRequestCommentLayer = new esri.layers.FeatureLayer(serviceRequestLayerInfo.CommentsLayerURL, {             mode: esri.layers.FeatureLayer.MODE_SELECTION,             outFields: [serviceRequestLayerInfo.CommentsOutFields],             id: serviceRequestLayerInfo.Key + "Comments"         });         map.addLayer(serviceRequestCommentLayer);     }          //function to showinfowindow     function ShowServiceRequestDetails(mapPoint, attributes) {         infoWindowDescriptionFields = [];         selectedServiceStatus = attributes.STATUS;         if (!attributes.REQUESTID || attributes.REQUESTID == showNullValueAs) {             attributes.REQUESTID = showNullValueAs;             map.infoWindow.setTitle(" Service Request ID:" + attributes.REQUESTID + "");         }         else {             map.infoWindow.setTitle(" Service Request ID: #" + attributes.REQUESTID + "");         }         var mainTabContainer = CreateServiceRequestTabContainer(attributes);         map.infoWindow.setContent(mainTabContainer.domNode);         mainTabContainer.resize();         var windowPoint = map.toScreen(mapPoint);         map.infoWindow.resize(310, 200);         setTimeout(function () {             map.infoWindow.show(mapPoint, GetInfoWindowAnchor(windowPoint, 310));             mainTabContainer.resize();             for (var index in infoWindowDescriptionFields) {                 CreateScrollbar(dojo.byId(index), dojo.byId(infoWindowDescriptionFields[index]));             }                  CreateScrollbar(dojo.byId('divDetailsContainer'), dojo.byId('divDetailsContent'));         }, 1000);              FetchRequestComments(attributes.REQUESTID);         FetchAttachmentDetails(attributes[map.getLayer(serviceRequestLayerInfo.Key).objectIdField]);         RemoveChildren(dojo.byId("divAttachmentsData"));         RemoveChildren(dojo.byId("divCommentsContent"));              CreateRatingWidget(dojo.byId('commentRating'));         ToggleCommentsView(false);     }          //function to fetch attachment details     function FetchAttachmentDetails(objectID) {         map.getLayer(serviceRequestLayerInfo.Key).queryAttachmentInfos(objectID, function (files) {             var fileTable = document.createElement("table");             var fileTBody = document.createElement("tbody");             fileTable.appendChild(fileTBody);                  for (var i = 0; i < files.length; i++) {                 if (files[i].contentType.indexOf("image") != 0) {                     fileTBody.appendChild(CreateData(files[i].name, files[i].url, files[i].size, files[i].contentType));                 }             }                  for (var i = 0; i < files.length; i++) {                 if (files[i].contentType.indexOf("image") >= 0) {                     fileTBody.appendChild(CreateData(files[i].name, files[i].url, files[i].size, files[i].contentType));                 }             }             fileTable.appendChild(fileTBody);                  if (files.length == 0) {                 var trNoAttachments = document.createElement("tr");                 var tdNoAttachments = document.createElement("td");                 trNoAttachments.appendChild(tdNoAttachments);                 tdNoAttachments.innerHTML = "No attachments found.";                 fileTBody.appendChild(trNoAttachments);             }             dojo.byId("divAttachmentsData").appendChild(fileTable);                  if (file.length > 0) {                 CreateScrollbar(dojo.byId("divAttachments"), dojo.byId("divAttachmentsData"));             }         });     }          //Function for creating data for the attachments     function CreateData(text, attachmentURL, fileSize, contentType) {         var filetr = document.createElement("tr");         var filetd = document.createElement("td");         if (contentType.indexOf("image") >= 0) {             var filePreview = dojo.create("img");             filePreview.style.width = "275px";             filePreview.style.height = "275px";             filePreview.style.cursor = "pointer";             filePreview.src = attachmentURL;             filePreview.onclick = function () {                 window.open(attachmentURL);             }             filetd.appendChild(filePreview);         }         else {             var filespan = document.createElement("span");             filespan.id = "fileTab";             filespan.innerHTML = text;             filespan.onclick = function () {                 window.open(attachmentURL);             }                  filespan.className = 'spanFileDetails';             filetd.appendChild(filespan);         }         filetr.appendChild(filetd);         return filetr;     }          //function to fetch service request comments     function FetchRequestComments(requestID) {         ShowDojoLoading(dojo.byId("divComments"));         var query = new esri.tasks.Query();         query.where = "REQUESTID = '" + requestID + "'";         query.outFields = ["*"];         //execute query         map.getLayer(serviceRequestLayerInfo.Key + "Comments").selectFeatures(query, esri.layers.FeatureLayer.SELECTION_NEW, function (features) {             var commentsTable = document.createElement("table");             commentsTable.style.width = "95%";             var commentsTBody = document.createElement("tbody");             commentsTable.appendChild(commentsTBody);             dojo.byId("divCommentsContent").appendChild(commentsTable);             if (features.length > 0) {                 features.sort(SortResultFeatures);      //function to sort comments based on submitted date                 for (var i = 0; i < features.length; i++) {                     var trComments = document.createElement("tr");                     var commentsCell = document.createElement("td");                     commentsCell.className = "bottomborder";                     commentsCell.appendChild(CreateCommentRecord(features[i].attributes, i));                     trComments.appendChild(commentsCell);                     commentsTBody.appendChild(trComments);                     CreateRatingWidget(dojo.byId('commentRating' + i));                     SetRating(dojo.byId('commentRating' + i), features[i].attributes.RANK);                 }             }             HideDojoLoading();             CreateCommentsScrollBar();         }, function (err) {             HideDojoLoading();         });     }          //function to create comment record     function CreateCommentRecord(attributes, i) {         var table = document.createElement("table");         table.style.width = "100%";         var tbody = document.createElement("tbody");         var tr = document.createElement("tr");         tbody.appendChild(tr);              var td = document.createElement("td");         var td3 = document.createElement("td");         td.innerHTML = "Importance: ";         td.style.width = "25%";         td3.appendChild(CreateRatingControl(true, "commentRating" + i, 0, 5));         td3.style.width = "100px";         var td1 = document.createElement("td");              var date = new js.date();         var utcMilliseconds = Number(attributes.SUBMITDT);              td1.innerHTML = "Date: " + dojo.date.locale.format(date.utcToLocal(date.utcTimestampFromMs(utcMilliseconds)), { datePattern: "MMM dd, yyyy", selector: "date" });         td1.style.width = "150px";              tr.appendChild(td);         tr.appendChild(td3);         tr.appendChild(td1);              var tr1 = document.createElement("tr");         var td2 = document.createElement("td");         var divComments = dojo.create("div", { "class": "wordBreakComments" }, td2);         divComments.style.width = "270px";         td2.colSpan = 3;         if (attributes.COMMENTS) {             divComments.innerHTML = attributes.COMMENTS;         }         else {             divComments.innerHTML = showNullValueAs;         }         tr1.appendChild(td2);         tbody.appendChild(tr1);              table.appendChild(tbody);         return table;     }          //function to add service request comment     function AddRequestComment() {         var text = dojo.byId('txtComments').value.trim();         if (text == "") {             dojo.byId('txtComments').focus();             ShowSpanErrorMessage('spanCommentError', messages.getElementsByTagName("textComment")[0].childNodes[0].nodeValue);             return;         }         if (text.length > 250) {             dojo.byId('txtComments').focus();             ShowSpanErrorMessage('spanCommentError', messages.getElementsByTagName("textCommentLimit")[0].childNodes[0].nodeValue);             return;         }         ShowDojoLoading(dojo.byId("divComments"));         var commentGraphic = new esri.Graphic();         var date = new js.date();         var attr = {             "REQUESTID": dojo.byId('spanInfoTitle').title,             "COMMENTS": text,             "SUBMITDT": date.utcMsFromTimestamp(date.localToUtc(date.localTimestampNow())),             "RANK": Number(dojo.byId('commentRating').value)         };         commentGraphic.setAttributes(attr);         map.getLayer(serviceRequestLayerInfo.Key + "Comments").applyEdits([commentGraphic], null, null, function (msg) {             if (msg[0].error) {             }             else {                 var table = dojo.query('table', dojo.byId("divCommentsContent"));                 if (table.length > 0) {                     var tr = table[0].insertRow(0);                     var commentsCell = document.createElement("td");                     commentsCell.className = "bottomborder";                     commentsCell.title = attr.COMMENTS;                     var index = dojo.query("tr", table[0]).length;                     if (index) {                         index = 0;                     }                     commentsCell.appendChild(CreateCommentRecord(attr, index));                     tr.appendChild(commentsCell);                     CreateRatingWidget(dojo.byId('commentRating' + index));                     SetRating(dojo.byId('commentRating' + index), attr.RANK);                 }             }             ToggleCommentsView(false);             HideDojoLoading();             CreateCommentsScrollBar();         }, function (err) {             HideDojoLoading();         });         ResetCommentFields();     }          //Function for sorting comments according to date     function SortResultFeatures(a, b) {         var x = a.attributes.SUBMITDT;         var y = b.attributes.SUBMITDT;         return ((x > y) ? -1 : ((x < y) ? 1 : 0));     }          //function to create tab layout for Service request     function CreateServiceRequestTabContainer(attributes) {         var tabContent = document.createElement('div');         tabContent.id = 'tabContent';              for (var i in attributes) {             if (!attributes[i]) {                 attributes[i] = "";             }         }              var dtlTab = new dijit.layout.ContentPane({             title: "Details",             content: CreateServiceRequestDetails(attributes)         }, dojo.byId('tabContent'));         var cmntTab = new dijit.layout.ContentPane({             title: "Comments",             content: CreateCommetsContainer()         }, dojo.byId('tabContent'));         var attTab = new dijit.layout.ContentPane({             title: "Attachments",             content: CreateAttachmentContainer()         }, dojo.byId('tabContent'));              dojo.connect(cmntTab, "onShow", function () {             CreateCommentsScrollBar();         });              dojo.connect(attTab, "onShow", function () {             CreateScrollbar(dojo.byId("divAttachments"), dojo.byId("divAttachmentsData"));         });              var tabContainer = document.createElement('div');         tabContainer.id = 'divTabContainer';         var tabs = new dijit.layout.TabContainer({             style: "width: 300px; height: 162px; vertical-align:middle;",             tabPosition: "bottom"         }, dojo.byId('divTabContainer'));         tabs.addChild(dtlTab);         tabs.addChild(cmntTab);         tabs.addChild(attTab);              tabs.startup();              return tabs;     }


"Poor" quality when zooming out a WMS-Layer based on GeoTiff data

i am having some issues when using the Geoserver Web Mapping Service to display GeoTIFF-Files. When zooming out (in an OpenLayers-Client or in the Geoserver Layer Preview for instance) the quality seems to be very poor. When opening the GeoTIFF in a picture viewer the quality seems to be ok ![screen shot](http://i.stack.imgur.com/6HCQu.jpg) I am working on a project where i have to display very old maps in an webbased application, so the quality is quite important. I have alredy tried to alter the compression levels of the wms in Geoserver, in vain. I'd be grateful if somebody could come up with a hint or an advise. Thank you.!


How can I validate a table name on SDE with Oracle backend?

_Second title: "Is arcpy.ValidateTableName broken on Oracle?"_ I have an Oracle SDE database 9.3.1 (service connection). My arcGIS client tools are version 10.0 SP5). I have some names of feature classes(fc) that I would like create. I want to make sure that the name of the feature class is valid before I create it so I use the `arcpy.ValidaTableName` function. However when I try to validate a long fc I do NOT get a valid table name for Oracle (Oracle has a table identifier character limit on 30).               import arcpy     oracleWorkspace = "F:\shr016\all_is_sde_databases\esri_mogpros_5167_gisowner.sde"          print arcpy.ValidateTableName("very_long_table_name_with_too_many_characters_for_an_oracle_database", oracleWorkspace)          # prints: very_long_table_name_with_too_many_characters_for_an_oracle_database      Is there way I can validate table names on SDE with a Oracle backend? **Update** It seems like this feature is indeed not working correctly. Esri support writes in an email: _I have tested with arcpy.ValidateTableName with an Oracle database and it does behave as you have mentioned in your previous email. Currently, I am not aware of any alternate way of checking if the feature class name is valid from the Oracle side. I will be investigating this issue further and will update you with my findings._ **Update 2** Esri support has followed up: _This tool was designed to check for naming conventions, i.e. presence of any invalid characters or use of any reserverd keywords specific to the RDBMS being used. For checking the character length, I am contemplating logging an enhancement request. In the meantime you could include something like the following in your script to make sure the table name does not exceed 30 characters:_               import string, os          test = "tableName"     oracleWorkspace = "F:\shr016\all_is_sde_databases\esri_mogpros_5167_gisowner.sde"          if len(test) <=30:       print arcpy.ValidateTableName(test, oracleWorkspace)      else:        print "Invalid table name"      **Conclusion:** If you need a function that returns a valid table name for an Oracle SDE table name you have to write it yourself. **Update 3** A final update from esri support which confirms the conclusion: _I have logged an enhancement request for this incident: NIM090928-arcpy.ValidateTableName should check that table names do not exceed the workspace's maximum number of characters._


convert GPS to x,y using reference node

i am using sensor nodes and their antenna to calculate GPS coordinates as notation of (degrees/minutes/seconds), i have some question of how to do these things:   1. how can i convert latitude and longitude to X,Y coordinate, then plot X,y points on a scale, also need to specify one node as reference one to calculate xy related to it.    2. what is the Euclidean distance between GPS position, or it is applicable to calculate distances between the converted XY points ?   3. i need equations to convert GPS to XY, please any helpful links !!


ComException -2147467259 Arcgis 10 VB.net +Geoprocessing

I have a problem when i execute GeoProcessing in ArcEngine 10 with VB.net   1. i had create the toolbox : Models.tbx   2. and Model Process : disolveRoute with 2 parameters  ![enter image description here](http://i.stack.imgur.com/2wniJ.png) My Function is :                Sub diso()         Dim dissolveTable As ITable = CType(featureLayerRoute, ITable)         ' Create geoprocessor. Overwrite true will replace existing output.         Dim gp As IGeoProcessor2 = New ESRI.ArcGIS.Geoprocessing.GeoProcessor()         Dim trac As ITrackCancel = New CancelTracker         Dim Result As ESRI.ArcGIS.Geoprocessing.IGeoProcessorResult2         Dim parameters As IVariantArray = New VarArrayClass()         gp.SetEnvironmentValue("workspace", "Database Connections\IFE connexion.sde")         gp.OverwriteOutput = True         gp.AddToolbox("data\Models.tbx")         parameters.Add(dissolveTable)         parameters.Add("typeroute")         Try             Result = CType(gp.Execute("disolveRoute", parameters, trac), ESRI.ArcGIS.Geoprocessing.IGeoProcessorResult2)         Catch ex As COMException             Console.WriteLine(ex.ErrorCode)         End Try          End Sub      For the inpute feature for disolve --> featureLayerRoute :                 Dim Sp As IWorkspace = ConnectToTransactionalVersion("localhost", "sde:postgresql:localhost", "sde", "....", "....")             featureLayerRoute.FeatureClass = featureWorkspace.OpenFeatureClass("piste_route")             Dim layerrouteAs ILayer = featureLayerDroit             layerroute.Name = featureLayerRoute.FeatureClass.AliasName      but when i execute this function it returned a Code ERROR -2147467259.. Please it's one Week i'm looking to solve it but Nothing :'( Error HRESULT E_FAIL has been returned from a call to a COM component


Geoserver proxy extention

i got geoserver 2.1.2 and when i wanna add the jar file ( for proxy extention for geoserver) in the desired file "geoserver 2.1.2/webapps/geoserver/WEB- INF/lib" and clicked to "Proxy Admin Page" in the GeoServer Web UI i got this error : HTTP ERROR: 500 class org.geoserver.web.proxy.HostnameProvider overrides final method model.(Ljava/lang/Object;)Lorg/apache/wicket/model/IModel; RequestURI=/geoserver/web/;jsessionid=p5z1wj01ackj Caused by: java.lang.VerifyError: class org.geoserver.web.proxy.HostnameProvider overrides final method model.(Ljava/lang/Object;)Lorg/apache/wicket/model/IModel; at java.lang.ClassLoader.defineClass1(Native Method) at java.lang.ClassLoader.defineClassCond(Unknown Source) at java.lang.ClassLoader.defineClass(Unknown Source) at java.security.SecureClassLoader.defineClass(Unknown Source) at java.net.URLClassLoader.defineClass(Unknown Source) at java.net.URLClassLoader.access$000(Unknown Source) at java.net.URLClassLoader$1.run(Unknown Source) at java.security.AccessController.doPrivileged(Native Method) at java.net.URLClassLoader.findClass(Unknown Source) at org.mortbay.jetty.webapp.WebAppClassLoader.loadClass(WebAppClassLoader.java:366) at org.mortbay.jetty.webapp.WebAppClassLoader.loadClass(WebAppClassLoader.java:337) at java.lang.Class.getDeclaredConstructors0(Native Method) . Im looking forward this error already one week and i cannot find any suggestions or solutions to handle this error :s Need HELP please !!


How do I transform a raster image to one with rounded corners (cylindrical to pseudocylindrical)?

I've been working with the altimetry data from topex.ucsd.edu/WWW_html/mar_grav.html Was having difficulties extracting the data, so I wrote a simple tool to do that. Creates a bitmap form the data source (result http://imgur.com/MV3ivxl) The problem is: I want an image with rounded corners (is there a name for that?) ![rounded corners](http://i.stack.imgur.com/jRGM1.jpg) but when I use gdalwarp, I'm getting a rectangle ![rectangle](http://i.stack.imgur.com/cqXOD.jpg). (actual result from gdalwarp command below: http://imgur.com/PZ2TEOe so you can see the projection is changing) Trying to go from Mercator to Robinson:               # data readme says source is spherical Mercator, from 0,80.738 to 360,-80.738     $ gdal_translate -a_srs '+proj=merc +datum=WGS84' -a_ullr 0 80.738 360 -80.738 grav.bmp grav.tif     # save http://spatialreference.org/ref/esri/54030/ogcwkt/ to robinson_ogcwkt.txt     $ gdalwarp --debug on -wo SOURCE_EXTRA=125 -t_srs robinson_ogcwkt.txt grav.tif grav_robinson.tif      output from last command               GDAL_netCDF:      =====     Open(), filename=[grav_robinson.tif]     GDAL_netCDF:      =====     Open(), filename=[grav_robinson.tif]     GDAL: GDALOpen(grav.tif, this=0x1e00740) succeeds as GTiff.     OGRCT: PROJ >= 4.8.0 features enabled     OGRCT: Source: +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs      OGRCT: Target: +proj=robin +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs      OGRCT: Source: +proj=robin +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs      OGRCT: Target: +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs      OGRCT: Source: +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs      OGRCT: Target: +proj=robin +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs      OGRCT: Source: +proj=robin +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs      OGRCT: Target: +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs      Creating output file that is 6157P x 3148L.     GDAL: GDALDriver::Create(GTiff,grav_robinson.tif,6157,3148,3,Byte,(nil))     Processing input file grav.tif.     WARP: Copying metadata from first source to destination dataset     GDAL: GDALDefaultOverviews::OverviewScan()     GDAL: GDALWarpKernel()::GWKNearestNoMasksByte()     Src=0,0,2825x4320 Dst=0,0,3078x3148     0...10...20...30...40...GDAL: Potential thrashing on band 1 of grav.tif.     GDAL: GDALWarpKernel()::GWKNearestNoMasksByte()     Src=2574,0,2826x4320 Dst=3078,0,3079x3148     50...60...70...80...90...100 - done.     GDAL: Potential thrashing on band 1 of grav_robinson.tif.     GDAL: GDALClose(grav.tif, this=0x1e00740)     GDAL: 8640 block reads on 4320 block band 1 of grav.tif.     GDAL: GDALClose(grav_robinson.tif, this=0x1e05610)     GDAL: 12592 block reads on 3148 block band 1 of grav_robinson.tif.      additional info:               $ gdalinfo grav.tif      Driver: GTiff/GeoTIFF     Files: grav.tif     Size is 5400, 4320     Coordinate System is:     PROJCS["unnamed",         GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,             AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]],         PROJECTION["Mercator_1SP"],         PARAMETER["central_meridian",0],         PARAMETER["scale_factor",1],         PARAMETER["false_easting",0],         PARAMETER["false_northing",0],         UNIT["metre",1,         AUTHORITY["EPSG","9001"]]]     Origin = (0.000000000000000,80.738000000000000)     Pixel Size = (0.066666666666667,-0.037378703703704)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  (   0.0000000,  80.7380000) (  0d 0' 0.01"E,  0d 0' 2.63"N)     Lower Left  (   0.0000000, -80.7380000) (  0d 0' 0.01"E,  0d 0' 2.63"S)     Upper Right (     360.000,      80.738) (  0d 0'11.64"E,  0d 0' 2.63"N)     Lower Right (     360.000,     -80.738) (  0d 0'11.64"E,  0d 0' 2.63"S)     Center      ( 180.0000000,   0.0000000) (  0d 0' 5.82"E,  0d 0' 0.01"N)     Band 1 Block=5400x1 Type=Byte, ColorInterp=Red     Band 2 Block=5400x1 Type=Byte, ColorInterp=Green     Band 3 Block=5400x1 Type=Byte, ColorInterp=Blue          $ gdalinfo grav_robinson.tif      Driver: GTiff/GeoTIFF     Files: grav_robinson.tif     Size is 6157, 3148     Coordinate System is:     PROJCS["World_Robinson",         GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,             AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]],         PROJECTION["Robinson"],         PARAMETER["longitude_of_center",0],         PARAMETER["false_easting",0],         PARAMETER["false_northing",0],         UNIT["metre",1,         AUTHORITY["EPSG","9001"]]]     Origin = (0.000000000000000,78.093065584271017)     Pixel Size = (0.049619632572081,-0.049619632572081)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  (   0.0000000,  78.0930656) (  0d 0' 0.01"E,  0d 0' 2.63"N)     Lower Left  (   0.0000000, -78.1095378) (  0d 0' 0.01"E,  0d 0' 2.63"S)     Upper Right (     305.508,      78.093) (  0d 0'11.64"E,  0d 0' 2.63"N)     Lower Right (     305.508,     -78.110) (  0d 0'11.64"E,  0d 0' 2.63"S)     Center      ( 152.7540389,  -0.0082361) (  0d 0' 5.82"E,  0d 0' 0.00"S)     Band 1 Block=6157x1 Type=Byte, ColorInterp=Red     Band 2 Block=6157x1 Type=Byte, ColorInterp=Green     Band 3 Block=6157x1 Type=Byte, ColorInterp=Blue


ArcGIS 10, Converting from date to number

I am trying to transform a date field like `"00:05:40"` into an integer field (i.e. 340 sec). Does anyone know how I could do that conversion please?


Find out tight corners on roads

## Problem: How can I pinpoint tight corners on a road dataset that might cause longer/larger vehicles difficulty in driving around. ## Assumptions:   * Ignore turns into junctions. Just interested in corners along a single road.   * Road network is a polyline layer. No arcs involved.   * Road width is an attribute.   * A large vehicle can be defined as a semi-truck.  ## Ideal Solution: The end goal would be for a map showing potential areas on a road network that would be best avoided by larger vehicles. Ideally this would have a score value as opposed to a yes/no value, to signify how tricky the particular corner is, depending on the size of vehicle. I have likely not provided enough information, so feel free to ask more Qs via comments and ill supplement my Q with more info. For info, this was a question I was posed, and I had some ideas, but would love to hear what your thoughts are on a potential workflow to solve this question.


How do I stop polygons from a .kml layer disappearing as I zoom in on the layer?

I am relatively new to GIS stuff and have a query about .kml files in QGIS. I have downloaded the .kml filefrom https://www.google.com/fusiontables/DataSource?docid=1qorOc15jWIA- OFlNjLCbviAn6F4X8VX8MOvGu7A I import this to QGIS as a new vector layer. When I zoom in on the layer, the polygons that are close to the egde of the screen disappear. How do I stop this from happening? Thanks in advance, Tom


How to calculate the total ascent along a path

I am working on a series of footpath routes and have been asked if it possible to calculate the total ascent that a walker would have climbed if the travelled along each route? I know that it is simply not a case of max-min because the footpath is not just one single climb but has many stages where the path ascends and descends. Presently, I have the data as a series of polylines and an underlying dtm in ArcGIS but do not have a spatial analyst extension. If anybody has any suggestions that would be great because I am quite stuck on this and never like to be beaten! I have GRASS GIS and QGIS too so quite happy to consider external methods if this will help answer the question. Many thanks!


ArcMap Commands .NET code equivalent

What is the .NET code equivalent for this general ArcMap command syntax:               Application.Document.CommandBars.Find(ArcID.Query_ClearSelection).Execute()      I have been using this one below as a work around, however it is much more code than what I would like to use.               Dim pUID As New UID     Dim pCmdItem As ICommandItem     'Use the GUID of the Save command     pUID.Value = "{AB073B49-DE5E-11D1-AA80-00C04FA37860}"     'or you can use the ProgID     'pUID.Value = "esriArcMapUI.MxFileMenuItem"     pUID.SubType = 3     pCmdItem = m_App.Document.CommandBars.Find(pUID)     pCmdItem.Execute()      Thanks


WFS - MapServer - ASP.NET

Nowadays i'm working with Mapcript / Mapserver. I would like to know how to used WFS (Web Feature Service) with ASP.NET. I think i need to use WCF, right? Is there any tutorial or framework? Thanks!


Where to get 2010 Census Block data?

This might be a stupid question but I would like to download 2010 Census Block Data for several towns and I can't seem to find it on American Factfinder. Is it there? And if not, does anyone know where I could find it?


Where could I download Income data at the Block Group level... census data

Where could I download Income data at the Block Group level, actually I'm not sure where or how to download it at the tract level either. I would be interested in the most recent data that is available.


QGIS renders outlines where no should be

QGIS always renders the borders of a feature polygon. How can I disable this? Even if I set Border style to "No Pen" and the colour is exactly the same for different categories there is a small grey border between adjoining features. Unfortunately this border is even rendered when using the map composer:![qgis renders ghost borders](http://i.stack.imgur.com/g0W3S.png)


Issue with Oracle Connection on Windows 7 machine

I have an extension for Arcmap 931 that connects to Oracle 10 client for data access. In windows XP, it works fine. If I install the extension on Windows 7 machine. It give this error. I am not sure what is missing here. I made sure that we have oracle client for 64 bit installed on this machine. _Attempt to load Oracle client libraries threw badImageFormatException. This problem will occur when running in 64 bit mode with the 32 bit Oracle client component installed._ Looking for any feedback. Thanks Jay ![enter image description here](http://i.stack.imgur.com/Irw4v.jpg)


Point input for polyline not acceptable?

Input QgsPoint(s) are returning an error in a Polyline creation routine. Any clues?!                   processing.runandload("qgis:polygoncentroids", Base, "memory:CenterNode")         processing.runandload("qgis:extractnodes", Base, "memory:OuterNodes")         Center = ftools_utils.getMapLayerByName(unicode("memory:CenterNode"))         Outer = ftools_utils.getMapLayerByName(unicode("memory:OuterNodes"))              # create line mesh points         line_center = []         featuresC = Center.getFeatures()         for elem in featuresC:             xy1 = elem.geometry().asPoint()             line_center.append(xy1)         print line_center              line_ends = []         featuresO = Outer.getFeatures()         for elem in featuresO:             xy2 = elem.geometry().asPoint()             line_ends.append(xy2)         print line_ends              Mesh = QgsVectorLayer("LineString", "Mesh", "memory")         prM = Mesh.dataProvider()         Mesh.startEditing()         for i in range(len(line_ends)):             print line_ends[i]             print line_center[0]             seg = QgsGeometry.fromPolyline(line_ends[i], line_center[0])             prM.addFeatures([seg])             Mesh.updateExtents()         Mesh.commitChanges()         QgsMapLayerRegistry.instance().addMapLayers([Mesh])         self.iface.mapCanvas().refresh()      ![enter image description here](http://i.stack.imgur.com/E3A0C.png)


50cm imagery - aerial or satellite?

The area of interest is 1450 km2. My options are: 50cm newly collected aerial natural color imagery for $28,300 OR 50cm recent archival satellite imagery (includes near IR band) for $20,600 The primary use will be land appraisal, geology(rock outcrop identification), and general GIS work such as digitizing of hydrology and transportation features, orientation maps, etc. The company hired to fly the area and process the imagery is having various technical difficulties so our product has not yet been delivered. Recent satellite imagery is available for the entire area. I am aware of the pan-sharpening multispectral imagery process and that aerial image (when available) will be show more detail on the ground but what other pros and cons are there? Which type of imagery is better for the uses I listed above?


Standard Style for Infrastructure

Hi: I'm trying to locate or find out if there is a standard style for municipal infrastructure. Specifically Canadian (Provincial BC, AB) standards for symbolizing water, sewer, valves, hydrants etc. Thanks, -Lori


Is it possible to read the map style and symbols from ArcGIS mxd

We used to use ArcGIS File Geodatabase(data source) and MXD file(for styling the map) to edit,style, and generate the final tiled map. Now we want to try with mapnik, but there are so many types of features,it is almost impossible to edit all the styles for mapnik. So I wonder if I can read the Arcgis MXD file, and extract the styles and symbols, then I can transfer them to the style format of mapnik. Is this possible? BTW, it would be better if this can be implemented without any ESRI components.


How to use Geokettle to load shapefiles into oracle spatial

I'm looking for an example how to load a shapefile from the file system into an oracle spatial database. What I have so far is:   * I load the shapefile with the 'Shapefile File input' step   * I use the calculator step to transform the WKT from the shapefile into a Geometry   * I use the Table output to write the data into the database. I get an error from the Table output step which says 'Error setting value #6 [GEOM Geometry] on prepared statement (Geometry), java.sql.SQLException: invalid datatype The target column is of type SDO_GEOMETRY. Thanks D3


Join Fields from Multiple Tables to One Table in Modelbuilder

I have 40 tables with about 2000 row with a shared ID that I want to merge into 1 table with 2000 rows. I only want to retain a distance field calculated within each table. The distance field is uniquely named Dist_miles_[name] with "[name]" as a unique value for each table. This way in the final table I will have: Key value, Dist_miles_[nameofAdvisor1], Dist_miles_[nameofAdvisor2], Dist_miles_[nameofAdvisor3]... Dist_miles_[nameofAdvisor40]. I would like to use iterate table and join field, and delete field to remove all the extraneous fields. However, the specificity of the variable I want to keep poses problems because I can't select it in the join field tool. I was hoping the join field tool would take a wildcard and be able to pull only the Dist_miles* into the final output table, but it doesn't. Is there a workaround this? UPDATE: Used Polygeo's tip to create modelbuilder. Works fine for single table. Tried to iterate on tables but doing that created an error at Calculate tool: **"A field name was not found or there were unbalanced quotation marks. ERROR 999999: Error executing function. Failed to execute (Calculate Field)."** Thanks!


How to parse ArcGIS Server generated output XML?

I have a use case in which I need to parse the XML generated by an ArcGIS server, then push that data to another API. The ArcGIS server sits behind the firewall at our customer's site and exposes a REST through which we can query the database. The REST URL looks something like this:               http://SERVER_NAME/arcgis/rest/services/APP_NAME/MapServer      I plan on stepping through the DOM and extract the records in Java (probably using Xerces). All the fields except the geometry fields are pretty straightforward. The geometry fields are defined as type `esriFieldTypeGeometry` with the data being stored as base64 binary string. My question is: how do I parse the binary data in Java and convert it to a set of coordinates? Here is an extract of the XML:               <esri:Workspace xmlns:esri='http://www.esri.com/schemas/ArcGIS/10.1'         xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance'         xmlns:xs='http://www.w3.org/2001/XMLSchema'>         <WorkspaceDefinition xsi:type='esri:WorkspaceDefinition'>             ...         </WorkspaceDefinition>         <WorkspaceData xsi:type='esri:WorkspaceData'>             <DatasetData xsi:type='esri:TableData'>                 ...                 <DatasetType>esriDTFeatureClass</DatasetType>                 <Data xsi:type='esri:RecordSet'>                     <Fields xsi:type='esri:Fields'>                         <FieldArray xsi:type='esri:ArrayOfField'>                             ...                             <Field xsi:type='esri:Field'>                                 <Name>SHAPE</Name>                                 <Type>esriFieldTypeGeometry</Type>                                 <IsNullable>true</IsNullable>                                 <Length>4</Length>                                 <Precision>0</Precision>                                 <Scale>0</Scale>                                 <Required>true</Required>                                 <DomainFixed>true</DomainFixed>                                 <GeometryDef xsi:type='esri:GeometryDef'>                                     <AvgNumPoints>0</AvgNumPoints>                                     <GeometryType>esriGeometryPoint</GeometryType>                                     <HasM>false</HasM>                                     <HasZ>false</HasZ>                                     <SpatialReference xsi:type='esri:ProjectedCoordinateSystem'>                                         <WKT>PROJCS[...]</WKT>                                         <XOrigin>-20037700</XOrigin>                                         <YOrigin>-30241100</YOrigin>                                         <XYScale>10000</XYScale>                                         <ZOrigin>0</ZOrigin>                                         <ZScale>1</ZScale>                                         <MOrigin>0</MOrigin>                                         <MScale>1</MScale>                                         <XYTolerance>0.001</XYTolerance>                                         <ZTolerance>0.001</ZTolerance>                                         <MTolerance>0.001</MTolerance>                                         <HighPrecision>true</HighPrecision>                                         <WKID>102656</WKID>                                         <LatestWKID>2234</LatestWKID>                                     </SpatialReference>                                     <GridSize0>3575.0162431554045</GridSize0>                                 </GeometryDef>                                 <AliasName>SHAPE</AliasName>                                 <ModelName>System-maintained Geometry for a feature</ModelName>                             </Field>                         </FieldArray>                     </Fields>                     <Records xsi:type='esri:ArrayOfRecord'>                         <Record xsi:type='esri:Record'>                             <Values xsi:type='esri:ArrayOfValue'>                                 ...                                 <Value xsi:type='esri:PointB'>                                     <Bytes>AQAAAMAnD4uRGzJBAFYObYdpJ0E=                                     </Bytes>                                 </Value>                             </Values>                         </Record>                     </Records>                 </Data>             </DatasetData>         </WorkspaceData>     </esri:Workspace>      Thanks in advance for any suggestions you guys can give.


Can two versions of GRASS be installed in a computer and still run independently?

GRASS 6.4.0 is currently installed in my computer. The script that I am trying to run is present in 6.4.1. It is stated in my previous post. If I installed 6.4.1, will my 6.4.0 still run? I don't want to risk uninstalling 6.4.0 since my project is currently linked to it.


How to unproject WKT to WKT in .NET?

I have WKT in this format:               GEOMETRYCOLLECTION(POLYGON((491990.250718 6212873.59755,491492.160493 6212904.94117,491265.64715 6212790.54445,491319.517251 6212550.36334,491579.624164 6212637.87553,491611.164518 6212718.61909,492055.359904 6212718.66908,491990.250718 6212873.59755)))      It is currently in Euref98 projection and I need it converted to LongLat. I'm a .NET developer so I prefer a .NET solution. Is there a .NET library that can parse WKT text into an in-memory model and write the in-memory model back to WKT?


Import shapefile is not working in pyshp?

I am just trying to use "pyshp" package to open shapefile in Python and use the Attribute data for further development. Firstly I have downloaded pyshp package for python 2.7 and installed it in my computer. But, whenever I want to import shapefile, I got following error:               Runtime error      Traceback (most recent call last):       File "<string>", line 1, in <module>     ImportError: No module named shapefile


Vector Data for Disputed Island Areas (Paracel and Spratley Islands)?

I am looking for a free vector data set for the Paracel and Spratley Island archipelagos. Natural Earth, suggested on StackExchange, does not include these features, due to their small size. Has anyone seen such a layer, ideally including island names? My interest is not in their sovereignty, but rather in showing their location and potentially spatial relationship with coral reefs. As such, data could be in any language (although English is preferred). Other island groups in the region are easy to find, but these two seem to be "missing" from many data sets because they are not assigned to any one country.


How to deploy a QGIS Python Standalone Application with minimum installation size?

I am trying to create a Standalone QGis Application using python. I was able to create the nice looking shapeviewer application by following the steps here. ![enter image description here](http://i.stack.imgur.com/xxFMy.png) Now I would like to deploy this in another system. Considering this is a small application, how can i keep the installation to a smaller size in the target machine? I would like to explore the possibility of not having to install QGis to run such a simple application. Is this possible?


Copying Features from One ArcMap window to another, without losing Related Tables

In our current workflow we are working in a Personal Geodatabase with few features and related tables , many persons will be working in different schemes. When we need to merge these schemes\Pgdb to a single output , we are copying feature class to feature class . Obviously this changes the Object ID in the resulting output Geodatabase and we lose the related table Information. Is there any way to copy feature classes from different inputs Pgdbs and still preserve the related table ?


ArcObjects Java Calculate Line Length

I'm relatively new to programming and ArcObjects. I'm trying to simply calculate the length of a line (test.shp) using ArcObjects and Java. I've tried the code below with no luck. This is simple in ArcMap because one just needs to use 'calculate geometry' on the field within the attribute table. But I need code to do this. Any suggestions would be greatly appreciated!               FeatureClass fc1 = null;             ShapefileWorkspaceFactory shapefileWorkspaceFactory = new ShapefileWorkspaceFactory(); Workspace shapefileWorkspace = new Workspace(shapefileWorkspaceFactory.openFromFile("c:/PFRAS/TEMP", 0));                   fc1 = new FeatureClass(shapefileWorkspace.openFeatureClass("test.shp"));          Field fe = new Field();     fe.setName("Dist");     fe.setType(esriFieldType.esriFieldTypeDouble);     fe.setScale(2);     fe.setPrecision(16);     fc1.addField(fe);          IFeature feature = fc1.getFeature(0);     Line line = new Line();          //Assume that the underlying feature is a Line.     IGeometry geom = feature.getShape();         if (geom instanceof Line)         line = (Line)geom;                       double length = line.getLength();          int fieldIndex = fc1.findField("Dist");          FeatureCursor fCursor = new FeatureCursor(fc1.IFeatureClass_update(null, false));             Feature feature2 = (Feature) fCursor.nextFeature();             while (feature2 != null)              {                               feature2.setValue(fieldIndex, length);              fCursor.updateFeature(feature2);              feature2 = (Feature) fCursor.nextFeature();               }             fCursor.flush();             Cleaner.release(fCursor);


Setting Band Statistics on Raster Image (GeoTIFF)

After doing calculations on a raster and saving it, the image statistics are wrong and need updating. I find examples saying that this should do the job, like in this example:                dsband.WriteArray(data)      dsband.FlushCache()      dsband.GetStatistics(0,1)      However this does not seem to update my band statistics correctly, they seem ok in Python with data.max() and data.min(), but opening the file in QuantumGis (image properties>style>load min/max) shows the wrong stats. I only get new image statistics correctly by separately opening and setting statistics like this:               driver = gdal.GetDriverByName('GTiff')     driver.Register()     ds = gdal.Open(infile, gdal.GA_Update)     dsband = ds.GetRasterBand(1)     (newmin, newmax)= dsband.ComputeRasterMinMax(0)     (newmean, newstdv) = dsband.ComputeBandStats(1)     dsband.SetStatistics(newmin, newmax,newmean, newstdv)     dsband.FlushCache()            dsband = None       ds = None      So what's the difference between both and why does QGis only show the correct stats with the second code? thanks for any advice!


How to compare schemas from two File geodatabases?

We have template geodatabase, which we send out to our regional Data creators. They update these geodatabases and digitize all the data from their region in the predefined featureclasses. Now that we have got the File geodatabases, we find that often many of the featureclasses have been modified. Maybe a field might be added or removed. In some cases, the Feature-class itself is deleted or a new one created. I would like to get a report indicating which feature-classes and tables are changed, as well as the differing Fields in the common tables & featureclasses. I have already looked at several questions such as:   * Tools to compare unversioned copies of geospatial data   * Is there a good way to get a diff between an old and a new version of an SDE database? But the answers given in these Questions have not been helpful.


How to configure various client applications to display AVL data?

Given x number of GPS units transmitting NMEA formatted strings back to a server (xxx.xxx.xxx.xxx:xxxxx) is there a way to listen in on the traffic and parse NMEA formatted strings for display on a map? Could a combination of a COM port data emulator and the standard tools in the ArcMap GPS toolbar on a client machine get me there? Any thoughts generally (or specifically in regard to ArcGIS) on what I might do (if anything) if I know where the data is coming in?


I need help to find out more information about Topology Errors toolbox or Extension?

I need someone out there who can answer my question. I do not remember if I add this Topology Editor Extension or toolbox to my ArcMap. I would like to know more about this toolbox. Thanks. (Changed from Errors to Editor, apologized) Please do not confuse this with the Topology toolbox that comes with the ArcMap Thanks ! ![Topology Errors ](http://i.stack.imgur.com/ymZfo.png)


Steps to polygonize raster KML image

I suspect I'm trying to do the same thing as in Convert image in KMZ to polygon. Given a kml displaying an image within a given lat/long boundary, what are the steps I follow to use Quantum GIS to convert this KML to a vectorized polygon? The image is simply a solid color map -- I.E. the associated PNG simply has each pixel as a set color or as a transparency, and I'm intrerested in being able to query if an arbitrary lat/long is within the bounds of the vectorized polygon. Edit: It doesnt look like I can open the kml directly. Its raster image is a PNG. How do I set the lat/long bounds of the image if I simply open this PNG as a raster layer?


Simple Features and ArcGIS Domains

I am working with a number of detailed ArcGIS Geodatabases with domains used to enforce specific workflows, etc. This is all good and great, but I would like to move my data into either a PostgreSQL/PostGIS or SQL Azure datastore and would like the ESRI client APIs (Mobile SDK, Web APIs, and Desktop templates) to enforce workflows in a similar manner to Geodatabase domains. Obviously a domain is nothing more than a table constraint, but there is lots of built in ArcGIS functionality which makes life easier for the end user (mainly desktop users). This seems very timely since it appears that at 10.1 ESRI is expanding support for native spatial support (i.e. Simple Features). Will ESRI provide domain support within a native spatial DB from ArcGIS Desktop? Thoughts, suggestions, etc.


Is the ArcGIS Server JavaScript API free to use?

The 2008 blog post announcing the ArcGIS Server JS API says: > Everyone can use the ArcGIS JavaScript API. There is no fee for using the > ArcGIS JavaScript API or deploying an application built with it. But Esri's Web Mapping Pricing page says "Contact us for pricing". Which is correct?


Converting GRIB to Geotiff with gdal_translate in Python

I want to extract each band from a GRIB (gridded binary from http://data- portal.ecmwf.int/data/d/interim_daily/) in a single Geotiff image. I´m using Python with gdal_translate in a QGIS backround, but the dgal_translate command doesn´t work. The import of the GRIB-data works well, but in my loop about all bands, where I want to execute the gdal_translate command I get an error. The loop also works well, but it tells me, that the algorithm for the GDAL command is unknown (I imported the tool in my script). I get the information (in the QGIS console) while running: "Algorithm not found" (I get this information for each band, but the program continues). There is no other explicit error warning. My code:               import os      from osgeo import gdal, ogr, osr from          filename='C:\\Users\\mypath\\test1.grib' rasterInput =     gdal.Open(filename, GA_ReadOnly)          if rasterInput is None:         print 'Grib import failed'          print 'Grib imported'          bandCount = rasterInput.RasterCount print bandCount               format = "GTiff" driver = gdal.GetDriverByName(format)          for a in range(bandCount):                   outputFile = 'C:/Users/mypath/Geotiffs/%s.tif' %(a+1)              translate = processing.runalg('translate, -b %s, -co COMPRESS=NONE, -co BIGTIFF=IF_NEEDED, -of GTiff, %s, %s' %(a+1, rasterInput, outputFile))          del bandCount      I copied and pasted this command from the QGIS tool raster --> conversion --> translate, but it just works for one band not for all in range automatically. How can I do this per command-line? I´m a beginner at Python and GDAL and I can´t find any satisfying result. I hope that is explicite enough for some answers. Thanks a lot Mari


How to estimate the average width of a sidewalk?

I have a table with two polygons representing _urban block_ with and without the _sidewalk_ : `geomBlock` and `geomBlockWithSidewalk`. On my city we can approximate the polygons of `geomBlockWithSidewalk` to a `geomBlock` buffer of 2 meters. The standard OGC function that generates this buffer is               geomBlockWithSidewalk_false = ST_Buffer(geomBlock,2.0)      It is "false" because reality is very very different than this ideal scenery. My problem now is estimate the sidewalk width `w`, that is **How to estimate`w` on the best approximation for `st_buffer(geomBlock,w)`?** Here an illustration of typical urban blocks shapes and sidewalk variations. ![Urban blocks](http://gis.cambridgema.gov/DataDictionary/Datapages/..%5Cimages%5Cpavedsurfaces.jpg). * * * NOTICE: **If** the sidewalk representation are stored as "rectangular strips", the problem of this question is similar to " _How can I calculate the average width of a polygon?_", because both solutions need to solve the problem of " **estimation of theaverage width of a _rectangular strip_** ". * * * **Another usual representation for sidewalks** , used here, is aggregating the shape of the sidewalks around a block with the block... The _rectangular strip_ here can be also a ring (a polygon with a hole): the sidewalk is the polygon resulted from the spatial difference between geometries of the blockWithSideWalk (`g2`) and the block (`g1`) "`g2 - g1`", that is `ST_Difference(g2,g1)`. ![g1](http://i.stack.imgur.com/NzJad.png) ![Sidewalk as a rectangular ring strip](http://i.stack.imgur.com/pkhhT.png) Imagining cut the ring: the cut make it a "rectangular strip"... But not exactly, it is more precisely a "trapezoidal strip". ![enter image description here](http://i.stack.imgur.com/VkpHm.png) It is a new question. After @whuber closed this question as "duplicated", I posted the general question: "Is there a st_buffer inverse function?", where we can see that all the questions are very closely related,   * the "width of the beach of a little island" problem,   * the calculation of "the average width of a street" problem, and   * this (here) "average width of a sidewalk represented with the block" problem. * * * The solution can use the "hypothesis of retangular blocks with sides `L` and `H` >>`w`". It is a function,               CREATE FUNCTION sidewalk_width(geometry, geometry) RETURNS float AS $$         -- calculations using geomBlockWithSidewalk and geomBlock         -- ST_Within(geomBlock,geomBlockWithSidewalk) is true.     $$ LANGUAGE SQL immutable;


OpenLayers: Point icons disappearing

I have a very simple setup where I have a feature vector layer on top of an OSM layer in OpenLayers 2.13.1. However, I notice that the icons of my points in the vector layer begin to 'become invisible' despite their SVG DOM bounds being positioned where the icon should be, making them sort of invisible as I zoom in. The layer that holds my features is not overly exciting:               var ooiLayer = new OpenLayers.Layer.Vector('OOI Layer');      And here is how I create my feature vectors:               // Declared elsewhere but added for completeness sake     var ooiStyle = $.extend({}, OpenLayers.Feature.Vector.style['default'], {         graphicWidth: 25,         graphicHeight: 25,         fillOpacity:.7,         externalGraphic: 'img/ooi.png'     });                    function add(ooi) {         var wkt = new OpenLayers.Format.WKT(ooi);              // The following line looks more complicated than it is. Assume the variable holds something like "POINT (31.123 34.321)" after evaluation.         var wktData = ooi.entityInstancesGeometry[0].geometry.geometry.wellKnownText;         var vector = wkt.read(wktData);         if (ooi.entityTypeId == 14 && vector.geometry instanceof OpenLayers.Geometry.Point)             this.createArea.call(this, vector.geometry.x, vector.geometry.y, ooi);         else {             var actualVector = new OpenLayers.Feature.Vector(                 latLon(vector.geometry.x, vector.geometry.y), // just a helper function that applies the correct projection if required; works as intended, so no worries here                 ooi,                 $.extend({}, ooiStyle, {                     title: ooi.entityName,                     externalGraphic: graphicFor(ooi.entityTypeId) // returns an URL to a 25x25 PNG                  })             );                  ooiLayer.addFeatures(actualVector);         }     }      From the looks of it, it almost appears as if the first icon on the map sets the visual bounding box for all the other icons. I've created a small example to illustrate this behavior: http://youtu.be/_axFRI6Pw1U I am not sure what is wrong here. Any ideas are appreciated :( Regards,   Manny


How to implement 1:n joins with string attribute concatenation?

I have a question concerning pyQGIS. I'm a QGIS long-time (self-learned) user, but just a pyQGIS newbie. I want to implement 1:n spatial join with attribute concatenation in the table. I have a vector points layer with multiple points on the corresponding polygons, I already made a spatial join of the points with the polygons, so I have now a point layer with the key field of the polygon (column 10 in the table). Now I just want to search and compare the table attributes in column 10 and if they are the same, concatenate the text in column 5. Should be easy, but for me it isn't. at first I just wanted to compare the column 10 and if it is the same - print how many time it repeats, so my first approach:               #!/usr/bin/env python     # -*- coding: cp1250 -*-          layer = qgis.utils.iface.activeLayer()     features = layer.getFeatures()     count=1     for feature in features:         # fetch attributes         attrs = feature.attributes()         attrs1 = features.next()              # attrs is a list. It contains all the attribute values of this feature         print attrs[10]         print attrs1[10]         print "FID:"+str(feature.id())         if attrs[10] == attrs1[10]:             count+=1             print count         else:             count=1             print count      the program works, but it skips each 1 row. the next approach (with help of Krzysztof Maślanka):               #!/usr/bin/env python     # -*- coding: cp1250 -*-     layer = qgis.utils.iface.activeLayer()     features2 = layer.getFeatures()     features = layer.getFeatures()     powt=1     for feature in features:         attrs = feature.attributes()         for feature2 in features2:             # fetch attributes             attrs1 = feature2.attributes()             # attrs is a list. It contains all the attribute values of this feature             if attrs[10] == attrs1[10]:                 powt+=1                 print "FID:"+str(feature.id())                 print attrs[10]                 print attrs1[10]                 print powt             else:                 print "FID:"+str(feature.id())                 print attrs[10]                 print attrs1[10]                 powt=1                 print powt      The problem with the above program is that it should compare every row to each other, but the primary loop exits after the sub-loop is finished, so it just compares the 1st row to the table, then exits. Can anybody help with this? It would cover the lack of this in QGIS. OR maybe somebody knows how to make a multiple 1:n spatial join with attribute concatenation in GRASS GIS or some other OS4GEO software (I didn't find it, but why should one ram against open door if it exists :-)?


List Joined tables

I've a question about arcpy List* modules. I have a quite large amount of data on SDE Server. Earlier I made a list of the feature classes , and tables in the mxd-s placed there using da.Walk and ListFeatureClasses/ListTableViews. Now I'd like to list the joined tables too, but I've no further idea how to do it. Actually is it possible to get theese relations? Thank you for your help Reka


ArcGIS Server Data Extraction - Job does not exist or is inaccessible

I have followed the documentation and video to set up a data extraction geoprocessing service, which I am accessing with the ArcGIS Viewer for Flex (v. 2.3.1). When I run the Extract Data widget in the viewer I get an error: Job '(job folder name)' does not exist or is inaccessible. When I go to the arcgisjobs folder on the server the folder is there, along with the scratch folder and the output .zip file containing the data that was clipped. According to the messages.xml file in the scratch folder, the job finishes successfully. I have checked the ArcGIS Server and Windows logs and haven't found anything obvious that would point to an issue. Folder permissions should be correct, since the task can create the job folder, clip the data, and create the output .zip folder. The task will run successfully when used through ArcGIS Desktop. So has anyone experienced similar behavior with the data extraction tool or any other GP tasks? Any ideas of where I can look next? Environment:   ArcGIS Server 10, .NET, SP2 - Distributed installation   ArcGIS Desktop 10, SP2   Windows Server 2008 x64 R2      Thanks, Tim **Update:**   I did some testing today (in Firefox) and found a few more things:      1. The GP Service does _occasionally_ work, it seems. Every few tries I do get a prompt to download the extracted data, only to get a security sandbox violation:   http:// **SERVER IP** /Samples/DataExtract/Flex/widgets/DataExtract/DataExtractWidget.swf cannot load data from http:// **SERVER NAME** /arcgisjobs/sumner.dataextract_gpserver/j77643b1b114e458f82e6f416262aaf58/scratch/output.zip.      2. When I am prompted to download the file, Firebug shows that it cannot download my crossdomain.xml file. It shows the address as http:// **SERVER NAME** /crossdomain.xml      In both cases, it is trying to access a http://servername URL, and can't get there. Problem is, I can't figure out where it's getting this URL. I tried modifying the virtual directory URL's and republishing the services using the ArcGIS Server IP instead of machine name, with no change in behavior.


FME - Convert DWG to JPG

I am a beginner in Safe FME and i need to know how to Convert DWG files to JPG Files. Is it possible?


Address or reverse geocode a point shapefile

I have a shape file that contains points. I want to do a reverse geocode and give these points an address. Can someone please help me?


Can I convert a kmz file to kml?

I have a kmz file. Is there an easy way to convert it to KML? Can I just "unzip"?


Save user specified settings on flash drive

With some help in this forum I already managed to successfully setup Quantum GIS 1.8 on a flash drive (see my post on Running QGis from flash drive). QGis starts and runs stable on different computers, but I still have problems to port also user speficied settings. To port the settings I use the `[--configpath]` option as suggested in the QGis-Manual. According to this I changed the qgis.bat batch file this way to start QGis with a qgs-project file:               start "Quantum GIS" /B %OSGEO4W_ROOT%\apps\qgis\bin\qgis.exe --configpath %DRV_LTR%\qgisconfig\ %DRV_LTR%\Daten\TESTPROJEKT.qgs      In general this works fine - QGis starts with the determined project file on different working stations. And installed Plugins, for example, also appear installed when the flash drive is connected to another working station. But it seems, that the `[--configpath]` option does not port user speficied settings themselves: Installed Plugins do not appear checked when set to another workstation, removed toolbars on one machine appear on the other, so do by „Customization“ hidden parts of QGis. It seems, that these UI-settings still are saved locally. When removing all entries in the Windows registry, keys appear from new in `HKEY_CURRENT_USER\Software\QuantumGIS` after running QGis.   I was hoping that the `[--configpath]` option could also save user specified settings on the flash drive. The overall goal of my work is to supply non-GIS experienced project partners with an easy to handle GIS-Data-Viewer, that can be used without installation. Especially with the new „Customization“ dialog QGis 1.8 seems to be perfect fort that, if not used Menus and Toolbars could just be turned off or hidden. I am thankfull for any suggestions!


Does QGIS DB Manager plugin lock PostGIS transactions?

I do most of my geospatial analyses in PostGIS 2.0 using PgAdmin III. I use QGIS 1.8 for visualizing my data but also for querying my database using the DB Manager plugin. Sometimes (quite often actually...) **when I import tables from PostGIS into QGIS using DB Manager and then I run queries on these layers in pgAdmin, they keep on running forever**. The only way to bring the queries to stop is to shut QGIS down. Not sure if this is a bug in the plugin or if my PostGIS settings are wrong but I must say I'm well annoyed by this behaviour. Many thanks in advance.


ArcObjects: Scale text not changing map units from miles to feet

I am adding scale text to a page layout. The code works, but the map units on the text scale still appear as miles and not feet as coded. Here is the code snippet to add scale text....                           IEnvelope stEnvelope = new EnvelopeClass();                 stEnvelope.PutCoords(14.1, 1.85, 15.37, 2.04);                 IUID uid3 = new UIDClass();                 uid3.Value = "esriCarto.ScaleText";                 IScaleText pScaleText = new ScaleTextClass();                 INumericFormat pNumericFormat = new NumericFormatClass();                 pNumericFormat.UseSeparator = false;                 pNumericFormat.RoundingValue = 2;                 ITextSymbol stTextSymbol = new TextSymbolClass();                 stTextSymbol.Color = pRgbColor;                 stTextSymbol.Size = 10;                 pScaleText.NumberFormat = pNumericFormat as INumberFormat;                 pScaleText.PageUnits = esriUnits.esriInches;                 pScaleText.PageUnitLabel = "inches";                           //Assign the Map Units to be feet on the Text Scale                 pScaleText.MapUnits = esriUnits.esriFeet;                 pScaleText.MapUnitLabel = "feet";                                //pScaleText.Style = esriScaleTextStyleEnum.esriScaleTextCustom;                 pScaleText.Symbol = stTextSymbol;                 IFrameElement stFrameElement = pGraphicsContainer.FindFrame(map);                 IMapFrame stMapFrame = stFrameElement as IMapFrame;                 IMapSurroundFrame StMapSurroundFrame = sbMapFrame.CreateSurroundFrame(uid3 as UID, null);                 IElement stElement = StMapSurroundFrame as IElement;                 stElement.Geometry = stEnvelope;                 stElement.Activate(activeView.ScreenDisplay);                 pGraphicsContainer.AddElement(stElement, 0);      Any coding help would be appreciated.


What is the method for finding the percentage of common cells among 2 or more DTMs for the same area?

I have a study area and 6 different DTMs for that area from different sources and I would like to know the percentage of pixels with the same value. I suppose I will have to extract some statistic and make some calculations in Excel, but I don't know what to extract and how to do it. If somebody knows some methods or something similar please give some indications.


How can I plot GPS coordinates on Google Map

What is the best way to plot GPS coordinates on a Google Map? Are there any tools that can do that (written in python specifically)?


Is it possible to use more than two control point links when georeferencing CAD files in ArcMap?

I'm trying to georeference CAD to a building footprint in ArcMap 10 (although I've had this same problem in earlier versions). The problem is, it won't let me use more than two control point links, so I have to update the display before adding two more. Every time I update the display, the CAD is stretched and placed properly where the control points were, but the other ends of the CAD are off. Is this simply something I'll have to live with? Or is there a way I can use more control points, or possibly convert the CAD to a less restricted file format?


How do I load layer, shape or raster map files using the ESRI ArcGIS Silveright API?

Is there a way to use the silverlight API to load common maps formats off of my local machine?


Can I get some recommendations for an Open Source Trouble Ticket system?

When I research trouble ticket systems a lot of them are for programming. I am interested in something more general purpose. Some things that would be a bonus would be: -Integrated with Drupal -PostgreSQL based -web interface for Clients Opinions and advise is needed.


How to access Esri ArcPy modules from Enthought Canopy?

I have ArcGIS 10.1 running on my windows system and I have also installed the new 32-bit Canopy Express Enthought Python Distribution. Is it possible to load and use the Esri ArcPy modules from Canopy?


How to delete an edge in QGIS

I have used the Lines to polygons tool to turn a city boundary linestring feature into a polygon. However, the tool has connected some nodes more than one time, resulting in multiple edges ![city boundary with multiple edges](http://i.stack.imgur.com/Cj2aA.png). I would like to remove some edges but not the nodes, therefore I cannot use the Node tool. Is there any way to do this (or to keep the tool from connecting node more than one time)? Thanks


Reduce processing time in Erase function using Arcpy in IDLE

I am processing 18,000 feature classes (saved in a file GDB), each with ~500 records (polygon features) in ocean areas, and I am running a script to remove area overlapping with the coast, using the Erase function. The script is running successfully, but literally taking days to complete. Some files take seconds to process (I'm assuming this is due to the fact, which I've verified, that they are not overlapping with coastline) while some feature classes take over an hour to process. I am working locally on my machine, running the script in IDLE, with ArcGIS 10.1 (all licenses) on a 32-bit Windows 7, 3.33GHz Intel DuoCore, 4GB RAM computer. I'm learning the ropes with Python and would really appreciate tips on my script if there are steps I can take to speed up processing time. Thank you!               # Import modules     import arcpy     import sys     import os     from arcpy import env          arcpy.env.overwriteOutput = True     print arcpy.env.overwriteOutput          #Erase coastline from selected buffers.     workspace = "C:\\DIR\\points_buff\\"     arcpy.env.workspace = workspace     erase_output = workspace+"erased.gdb\\"     erasedirectorypath = "C:\\DIR\\points_buff\\selectedbuff.gdb"     fileList = []          for dirpath, dirnames, erases in arcpy.da.Walk(erasedirectorypath,                                                   datatype="FeatureClass",                                                   type="Polygon"):         for erase in erases:             fileList.append(os.path.join(dirpath, erase))             filenameparts = erase.split('_')             if int(filenameparts[1]) > int(1967):                      try:             #Set variables for Erase tool                     in_features = erasedirectorypath+"\\"+erase                     coastline_feat = "C:\\DIR\\Basic_NE.gdb\\Coastline_Eastern_US"                     out_erase_feat = erase_output+erase[:-4]+"_er"                     arcpy.Erase_analysis(in_features, coastline_feat, out_erase_feat, '')                     print arcpy.GetMessages()                     print "Coastline has been erased from " + erase + " and saved to " + erase_output                      except:                     print arcpy.GetMessages()


Select / Deselect and Identify Features in QGIS

Sorry if this is all a bit 101, but I'm new to QGIS. I basically have a large Arc project that I can access in QGIS for reference, so will only really be using it to identify features. For example to check the classification of a road line (which may be correct even if a lozenge on a map is incorrect, for example.) I have my project open, but am unable to select any features from layers that are checked. I only have the option to Deselect Features from All Layers. However, when I try to do that, nothing seems to happen. So I am unable to select any features I want to identify. This must be something pretty simple, but some help would be much appreciated for a newbie! Thank you.


Increase max zoom level in OpenLayers

I'm developing a map in an Android WebView, using OpenLayers and a WMS Service. I can only zoom a certain amount of times (which I would like to increase). My code is:               map = new OpenLayers.Map({div:"map"});     baseLayer = new OpenLayers.Layer.WMS("WMS", link, {layers:"basic"});     map.addLayer(baseLayer);      I've tried from the page Setting Zoom Levels               var options = { numZoomLevels: 25 };     map = new OpenLayers.Map( {div:"map"} , options);       but then the map doesn't even show. And this:               map = new OpenLayers.Map({div:"map"});     baseLayer = new OpenLayers.Layer.WMS("WMS", link, {layers:"basic", numZoomLevels: 25});     map.addLayer(baseLayer);      but in this case, nothing changes (I can only zoom in 16 times, as before). Where is my mistake?


How to toggle Fusion Tables layer on/off based on map scale?

In the Google Maps API v3 I'm displaying points derived from a Fusion Table by specifying a FusionTablesLayer Is it possible to set a minimum/maximum scale on this layer, so it turns off when zooming out? EDIT: I'm presuming that I'll need to listen to the map's zoom_changed event, and handle the visibility manually. **But how can I change the visibility of a Fusion Tables layer?** A **marker** has the boolean property Visible and the method SetVisible() but I can't see an equivalent for a layer. An **overlay** can have show/hide methods - do I need to build my Fusion Tables layer as an overlay? (I'm a Google Maps novice so I don't know much about overlays vs layers) Thanks


Retrieve speed and number of Lanes in Google Maps API, OSM data or any other stream?

Is it possible to return the posted speed and number of lanes of a highway segment using Google Maps API, OSM data, or any other stream (paid for or free)?


Publish ArcSDE FeatureService for Versioned Editing

**Scenario** : I need to Publish a data (from ORACLE ArcSDE) for versioned editing. Each user should be able to connect to his version and edit the data throu web application (Silverlight/JAVA/Flex API). **Question** : Do I need to create a FeatureService for each user (with ArcMap document containing data from version)? Or is there a simpler way to do it? I would prefer if I have to make only one ArcSDE connection in ArcMap and publish that resource as a featureservice.


Problem with editing representation

How we can edit the representation of a layer programmatically in DotNet ? ''Representations allow you to customize the appearance of features by storing symbol information with the feature geometry inside feature classes. This additional control can help you meet stringent cartographic specifications or simply improve the display.''


How can I track/graph arcserver usage statistics?

I’m currently using the ArcGIS Server Properties dialog in ArcCatalog to periodically check the usage statistics for our services. I’d like to be able to chart those to get a big picture view of what’s going on. Does anyone know the location on an ArcServer machine where a file with these statistics is being written, ideally something that could be pulled into a spreadsheet? Thanks for any help. ![enter image description here](http://i.stack.imgur.com/8xQZ7.jpg)


Finding TOC Layers Type (layer,grouplayer,..) and layerindex of a layer

How can i find the type of the layers in TOC ? these types are very important to me : 1.layer (that has a valid datasource) 2.layer (that lose its datasource because of changing the address of its datasource ) 3.GroupLayer And i want to know about the index of the layer that we have with ILayer. A function like this : Public Function GetLayerIndex(mLayer as ILayer) as long


Renaming Columns in Geotools FeatureType

I am attempting to build a simple java program that loads shapefiles into a Postgis 2.0 database. Everything is working out great except that I would like to be able to modify the FeatureType in the following ways: 1) change the geometry column name from "the_geom" to "geom" 2) set all columns to be lower case. 3) change the primary key from serial to bigserial. The following is my current code which basically is just transferring the shpSchema to the pgSchema. I am fairly certain that what I am trying to do is fairly simple and common, but I can find no obvious way to do it or any documentation --or even any examples of someone trying to do this.                   Map shpConnectParams = new HashMap();         shpConnectParams.put("url", shapefile.toURI().toURL());              DataStore shpStore = DataStoreFinder.getDataStore(shpConnectParams);              String shpTypeName = shpStore.getTypeNames()[0];         FeatureSource shpFS = shpStore.getFeatureSource(shpTypeName);              SimpleFeatureType shpSchema = (SimpleFeatureType) shpFS.getSchema();         CoordinateReferenceSystem crs = shpSchema.getCoordinateReferenceSystem();         Map pgParams = new HashMap();         out.println(shpSchema.getGeometryDescriptor().getType());         for (AttributeDescriptor attrDesc : shpSchema.getAttributeDescriptors()) {                  out.println(attrDesc.getName());         }              pgParams.put("dbtype", "postgis");                   //must be postgis         pgParams.put("host", "localhost");      //the name or ip address of the machine running PostGIS         pgParams.put("port", new Integer(5432));             //the port that PostGIS is running on (generally 5432)         pgParams.put("database", "test");              //the name of the database to connect to.         pgParams.put("user", "postgres");                    //the user to connect with         pgParams.put("passwd", "test!");                 //the password of the user.         pgParams.put("schema", "scratch");   //the schema of the database         pgParams.put("create spatial index", Boolean.TRUE);                        SimpleFeatureTypeBuilder builder = new SimpleFeatureTypeBuilder();              builder.addAll(shpSchema.getAttributeDescriptors());         builder.setCRS(crs);         builder.setName(shpTypeName);              DataStore pgDatastore = DataStoreFinder.getDataStore(pgParams);         pgDatastore.createSchema(builder.buildFeatureType());              SimpleFeatureSource pgFS = pgDatastore.getFeatureSource(shpTypeName);              SimpleFeatureType pgSchema = pgFS.getSchema();                   Transaction transaction = new DefaultTransaction("create");                   if (pgFS instanceof SimpleFeatureStore) {             SimpleFeatureStore pgFeatureStore = (SimpleFeatureStore) pgFS;                                      SimpleFeatureCollection features = (SimpleFeatureCollection) shpFS.getFeatures();             pgFeatureStore.setTransaction(transaction);             try {                 pgFeatureStore.addFeatures(features);                 transaction.commit();                  } catch (Exception problem) {                 problem.printStackTrace();                 transaction.rollback();                  } finally {                 transaction.close();             }             System.exit(0); // success!         } else {             System.out.println(shpTypeName + " does not support read/write access");             System.exit(1);         }              pgDatastore.dispose();


How to use CropCoverage Render Transformation in Geoserver SLD (WPS)?

I have yet to find an example on the web (or in code) of anyone using the "Crop Coverage" Render Transformation from a Geoserver SLD. Does anyone have an example or at least be able to point me to the syntax of using an existing Geoserver hosted layer as a the "crop geometry?" **gs:CropCoverage** is one of the WPS Processes used to crop raster layers. It is also possible to chain more WPS processes in a SLD but I cannot find an example of how to use gs:CropCoverage in a SLD. At the moment I am using Geoserver 2.4.3 on Tomcat 6.


Is it possible to see polyline as route or track data on GPS?

I would like to see polyline on GPS. First, I try to convert polyline to GPX. I uploaded GPX on GPS. SO far, I could not see the line on GPS. First of all,is it possible to see polyline as route or track data on GPS?


Overview of statistical maps/charts types

I know of a few ways to visualize statistical information on a map, like choropleth maps, prism maps or simply displaying little bar-charts or pie- charts near locations on a map. But there must be a lot more different ways to display statistical information on a map. I can only find lists with specific subsets, but no clear overview of all commonly used methods. Does an overview exist of all commonly used ways to map data to geographic location?


Getting geometry of drawing in ArcGIS Server JavaScript API

To me, ArcGIS Help and Online Resources is always full of unneccessay information and missing some important points. For example there are lots of samples for drawing but I couldn't find any thing for steps post drawing like saving on Database or Map! By the way, I am trying to Save some user inputs (Points or Multipoints) on the Map to MySQL database. I have a Table in database as :               create table Points (      name VARCHAR(20) PRIMARY KEY,      location Point NOT NULL,      description VARCHAR(200),      SPATIAL INDEX(location)      );       Now, I would like to save Geometry of the drawing point(s) in to this table. Can you please point me to an sample, tutorial or give me an idea how to get this information from the drawing classes?


Demographic Information of U.S. Census Blocks

I am struggling to find the blockwise (not blockgroup or tract) data like employment status, household income, commuting distance, list of place of interest, industry information. By using American Factfinder I get easily the updated data of these factors by tract and 2010 version of blockgroup. ACS Employment Status Data by Block Group, 2006-2010 is given here. Any help is appreciated.


Issue with recently cached ArcGIS services

ArcGIS for Server 10.0 SP5 Virtual Windows Server 2008 R2 Enterprise 16.0 GB RAM, 4 Processors (2.00 GHz each) Hello All, I am having issues with two services that I have cached as basemaps for my web applications. Here's a little back story - These services have been in use for a couple years, and every so often (about 6 months) these caches are rebuilt to address changes in basemap items (i.e. street centerline changes). We use a variety of different basemaps in our web map applications, and this weekend I recently updated the cache on two of them since it takes so long to do them in our current setup (one machine environment). These services worked just fine, then after recaching them they appear to be fine, until a user spends a little more than 2 minutes on an application using the service. Then the basemap tiles begin to lag and not load. Network debugging shows the pending processes, and SOMETIMES the tiles do load after about 1-2 minutes. MOST times they do not load for at least 4 minutes. Remaining times they just stop loading, and process stay pending until they timeout with a 500 server error. I have an open ticket with Esri, but they have yet to pinpoint the issue. Some changes we have tried include:   * Editing the impersonation for the rest and services config files from true to false   * Running post installation   * Changing services to low isolation and increasing instances to 8 Results for each change are as follows respectively:   * Caused Bad 404 errors for REST page.   * This was done in to try and fix above 404 error, it did not. Impersonate values changed back to true to fix 404 error.   * No changes seen with services. Same lag issues as before. In addition to the above, the ArcGISServiceAppPool process periodically will cause the CPU on the server to peg out at 100%. I have seen the KB articles (w3wp.exe, although it is for 9.x versions, and lsass.exe) related to this issue. Implementing the suggested fixes are what caused the 404 issue. I read through this post and it relates to some of the issues I am also experiencing. I have checked the log and do not see anything that helps identify the issue at hand. If anyone has experienced a similar situation and can provide any assistance, I would greatly appreciate it. Please let me know if you need additional details. Thanks in advance! EDIT - question is also posted to ArcGIS forum. User asked if service caches were being created on demand. Answer: They are not, however both with and without this option have been tested and the results are similar. UPDATE - I believe I may have found why the issue is occurring only at some levels. It appears that some tiles are missing at these levels. Since some tiles are missing, users still pan/zoom which then become network processes that never complete. If enough of these add up, it bogs the server down (since most other requests are fulfilled, it's hard to monitor performance). I think from here I am going to recache the map at the last three levels, then copy the updated levels to the original cache and see if this fixes the issue.


landsat with C# or Java

I have imaged data about landsat8 satellite. These are raster data, i want to process to take out the temperature of each region (geometry). I can do this in postgresql, but i don't want to store the raster data because they are too enormous. I can use c # or java to process that problem or not? From images of landsat8 through java software, can i collect vector data corresponding to each zone in database and then import postgrerql? I want to use 2 languages, but if they are impossible, are there any languages supporting for that? Thank you very much


Buffer along a route

I'm using Googlemaps service to plot some marker on the map. All markers are stored in a postgresql table, named "locations" and that has a location column (data type to geometry).   When the user clicks on the markers a route is created between them. My question is: how can I generate a buffer (with a limit of 10km along route) and select all my locations inside the buffer. I`m really new in postgis. P.S. This link inspired me but without success   Thanks!


Which function for creating a POINT in PostGIS?

When defining a Point in PostGIS, when do you decide to use which of the following?   * `ST_SetSRID(ST_MakePoint(lon,lat),4326)`   * `ST_SetSRID(ST_Point(long,lat),4326)`   * `ST_SetSRID(ST_GeomFromText('POINT(lon lat)',4326)`   * `ST_GeomFromEWKT('SRID=4326;POINT(lon lat)')` If its essentially a difference in performance, which will be the fastest?


How to manage the osmosis filters?

I am a complete beginner in OSM and osmosis, and I am trying to build a stripped down version of a local map, in order to add nodes of public lights that would be used for maintenance. I have used the following command line in order to keep only streets and roads, but there are still plenty of one way signs (arrows) cluttering the narrow streets in the center of the village, and I would like to get rid of them. I could not find any way to eliminate them (or to make them invisible in the output) I am also trying to get rid of the names of the streets, since most tend to consume too much space in some areas where I will need to add nodes. Any comment would be appreciated as well in order to refine or correct the command.               osmosis --read-xml file="Source-Map.osm.xml" --tag-filter accept-ways highway=* --tag-filter reject-ways highway=track --tag-filter reject-nodes highway=bus_stop --tag-filter reject-nodes highway=one_way --tag-filter reject-nodes barrier=* --tag-filter reject-nodes highway=untagged --tag-filter reject-relations --used-node --write-xml file="FilteredMap.osm"


Restricting base map on specific extent in ArcGIS JSAPI 3.3/3.4

Currently I am using ArcGIS JS API ver 3.3/3.5 in my application. I want to restrict the basemap for US only .For example one can use default navigation features like zoom in, zoom out for the US Country only. Trying following things: 1.I have set the US extent code in my application but when I zoom out then I can see the whole world map I want to restrict it to certain extent only 2.Going through this question Restrict pan and zoom in ArcGIS Javascript API? 3.Restricting the streets base map using this question but unable to restricting the map How do I use my own basemap with Javascript API? 4.I am using this sample map in my application. Any help regarding this will be great!!! Thanks in advance !!!


GeoTIFF without projection?

I have encountered a GeoTIFF raster files with EPSG:4326 (WGS84), i.e. no projection defined! Is this normal? I would expect raster to be in Projected Coordinate System. Raster is a grid of pixels so I expected this would be defined in the plane.


Rf-Classification seems to work but gives an error 15 seconds later

This is my code. It classifies an imagery-stack(xvars) and reads a shapefile with training-points. If I run single rows they work. Even if I run them all, they work and start to predict but after 15s it stops working and gives an error. As far as I understand, the classification doesn't need the raster-names because they are read automatically. My training-shapefile (sw_trainshape.shp) inherits a table with columns               FID    Shape    Class      ObjectID       x              y      1     Point     bush         1      481791,2429   5626286,6397      2      ...       ...        ...            ...      My tif-files are named _band1_ , _band2_ , _band3_ and so on. I have 6 bands. Yet I'm not experienced enough, so I could use some help as to why my code doesn't classifies. ERROR: > "Loading required package: tcltk > > Error in predict.randomForest(model, blockvals, ...) : variables in the > training data missing in newdata" Code:               setwd("D:/BA-Workspace/DOP_10/orthophotos_abcd/R/test_run_R/test_other")               library(sp)     library(rgdal)     library(raster)     library(randomForest)          # create list of rasters     rlist=list.files(getwd(), pattern="tif$", full.names=TRUE)           # CREATE RASTER STACK     xvars <- stack(rlist)                # READ Raster TRAINING DATA     sdata <- readOGR(dsn=getwd(), layer="sw_trainshape")          # ASSIGN RASTER VALUES TO TRAINING DATA     v <- as.data.frame(extract(xvars, sdata))     sdata@data = data.frame(sdata@data, v[match(rownames(sdata@data), rownames(v)),])          # RUN RF MODEL     rf.mdl <- randomForest(x=sdata@data[,3:ncol(sdata@data)],     y=as.factor(sdata@data[,"Class"]),                        ntree=501, importance=TRUE)          # CHECK ERROR CONVERGENCE     #plot(rf.mdl)          # PLOT mean decrease in accuracy VARIABLE IMPORTANCE     #varImpPlot(rf.mdl, type=1)          # PREDICT MODEL     predict(xvars, rf.mdl, filename="RfClassPred.img", type="response",          index=1, na.rm=TRUE, progress="window", overwrite=TRUE)      ![added sdata@data](http://i.stack.imgur.com/3IOgB.jpg) Console:               > setwd("D:/BA-Workspace/DOP_10/orthophotos_abcd/R/test_run_R/test_other")     >      >      > library(sp)     > library(rgdal)     > library(raster)     > library(randomForest)     >      >      > # CREATE LIST OF RASTERS     > rlist=list.files(getwd(), pattern="tif$", full.names=TRUE)      >      > # CREATE RASTER STACK     > xvars <- stack(rlist)           >      > # READ Raster TRAINING DATA     > sdata <- readOGR(dsn=getwd(), layer="sw_trainshape")     OGR data source with driver: ESRI Shapefile      Source: "D:/BA-Workspace/DOP_10/orthophotos_abcd/R/test_run_R/test_other", layer:     "sw_trainshape"     with 256 features and 10 fields     Feature type: wkbPoint with 2 dimensions     >      > # ASSIGN RASTER VALUES TO TRAINING DATA     > v <- as.data.frame(extract(xvars, sdata))     > sdata@data = data.frame(sdata@data, v[match(rownames(sdata@data), rownames(v)),])     >      > # RUN RF MODEL     > rf.mdl <- randomForest(x=sdata@data[,3:ncol(sdata@data)],     y=as.factor(sdata@data[,"Class"]),     +                        ntree=501, importance=TRUE)     >      > # CHECK ERROR CONVERGENCE     > #plot(rf.mdl)     >      > # PLOT mean decrease in accuracy VARIABLE IMPORTANCE     > #varImpPlot(rf.mdl, type=1)     > #setOldClass(SpatialPointsDataFrame)     > # PREDICT MODEL     > predict(xvars, rf.mdl, filename="RfClassPred.img", type="response",      +         index=1, na.rm=TRUE, progress="window", overwrite=TRUE)     Error in predict.randomForest(model, blockvals, ...) :        variables in the training data missing in newdata      ![enter image description here](http://i.stack.imgur.com/ozm81.jpg) Solution, thanks to TimSalabim:               setwd("D:/BA-Workspace/sw_west_aug/reduced_size/")               library(sp)     library(rgdal)     library(raster)     library(randomForest)               # CREATE LIST OF RASTERS     rlist=list.files(getwd(), pattern="tif$", full.names=TRUE)           # CREATE RASTER STACK     xvars <- stack(rlist)            # CREATE RASTER STACK     xvars <- stack(rlist)       x <- coordinates(xvars)[, 1]     y <- coordinates(xvars)[, 2]          x_rst <- y_rst <- xvars[[1]]     x_rst[] <- x     y_rst[] <- y          xvars <- stack(x_rst, y_rst, xvars)     names(xvars) <- c("X", "Y", "focal_1", "focal_2", "focal_3")     # READ Raster TRAINING DATA     sdata <- readOGR(dsn=getwd(), layer="training_west")          # ASSIGN RASTER VALUES TO TRAINING DATA     v <- as.data.frame(extract(xvars, sdata))     sdata@data = data.frame(sdata@data, v[match(rownames(sdata@data), rownames(v)),])          sdata@data  <- sdata@data[-c(5,6)]           # RUN RF MODEL     rf.mdl <- randomForest(x=sdata@data[,3:ncol(sdata@data)],   y=as.factor(sdata@data[,"class"]),                        ntree=501, importance=TRUE)          # CHECK ERROR CONVERGENCE     #plot(rf.mdl)     #sdata@data           # PLOT mean decrease in accuracy VARIABLE IMPORTANCE     #varImpPlot(rf.mdl, type=1)     #setOldClass(SpatialPointsDataFrame)     # PREDICT MODEL     predict(xvars, rf.mdl, filename="RfClassPred.img", type="response",          index=1, na.rm=TRUE, progress="window", overwrite=TRUE)


Using Spatial Analyst Classification in ArcGIS Engine?

I tend to use ML classification of Spatial Analyst ToolBox in ArcGIS Engine, I wrote this sample code:                       ESRI.ArcGIS.RuntimeManager.Bind(ESRI.ArcGIS.ProductCode.Desktop);             IAoInitialize aoInit = new AoInitialize();             esriLicenseStatus licStat = aoInit.Initialize(esriLicenseProductCode.esriLicenseProductCodeEngine);             licStat = aoInit.CheckOutExtension(esriLicenseExtensionCode.esriLicenseExtensionCodeSpatialAnalyst);                       Geoprocessor gp = new Geoprocessor();                  gp.AddToolbox(@"C:\Program Files (x86)\ArcGIS\Desktop10.1\ArcToolbox\Toolboxes\Spatial Analyst Tools.tbx");                  //--------------------------------------ENVI                  IVariantArray parameters = new VarArrayClass();             parameters.Add(inputtxt.Text);             parameters.Add(sigtxt.Text);                  parameters.Add(outputtxt.Text);                  IGeoProcessorResult result;             result = (IGeoProcessorResult)gp.Execute("MLClassify ", parameters, null);      but unfortunately the code stop with this error: Error HRESULT E_FAIL has been returned from a call to a COM component. I do not know how to handle it?


How to handle sensitive geographic data?

You use sensitive geographic data every day. What are your strategies to protect them in your geographical information systems? What kind of architecture do you use? What encryption method do you use? What do you do for users who export sensitive data from your databases? Thank you very much.


Do any geodatabase schema and data migration libraries for .NET exist?

I'm performing code refactoring and database normalization in my project and looking for any .NET library for geodatabase schema and data migrations. Does any solution exist for ArcObjects? During revision of my geodatabase schema, it may be that I discover that two of my tables contain an identical bunch of columns and they should be extracted into a separate table. So, I change my test geodatabase and program logic to conform to the new schema. But there are many production deployments with the old one. When a new version of the program is deployed to the client, it should perform a database upgrade. And I`m looking for .NET libraries for AO that are designed to write geodatabase upgrades (migrations) with less code (like Migrator.NET).


Improve gdaldem hillshade output (lighten up flat parts)

I have a hill-shading that I made using `gdaldem hillshade` which is already pretty good. I used the option `-alt 60` to get nice contrast, but the problem is, that with lowering the light-source, also parts in the valley that are quite flat, are not really bright (white). How could I improve that? I've seen that `gdaldem hillshade` doesn't have many options. What other program could I use to improve that? I'd like to say for example: Everything below 20% incline should be a bit lighter.