title: How can I add multiple rules in a GetLegendGraphic request in Geoserver?

text: I am using following GetLegendGraphcs request https://develop.blaa.nl/data/test/wms?LAYER=nccw:variabele_woz&request=GetLegendGraphic&format=image/png&styles=nccw_variabele_woz&rule=one With this I am able to get Legend for Phase 1. I have rules for Phase 2 (rule2) , Phase 3 (rule3) and phase 4 (rule4) too. If I want rules for Phase 1 and Phase 2, how shall I use GetLegendGraphics request where I can include more rules. Thanks & Regards


title: How to Create Touching Multi-Part Polygons?

text: I have a set of parcels as illustrated in the image below. There are currently no multi-part features in dataset. What I would like to do is to merge (not the right word) all of the parcels with same address, say 421 for example. I want to have these three polygons as one multi-part polygon, that still has the same geometry, as in it doesn't dissolve the lines. ![Image 1](http://i.stack.imgur.com/lFJvw.png) When I have tried this thus far, the lines dissolve resulting in only the boundary being preserved. ![Image 2](http://i.stack.imgur.com/wOSWi.png) I feel like there is a super easy was to do this, but I can't find it. I have full access to ArcGIS Advanced and use it to manage the topology of this dataset.


title: WFS filter query does not return points (selected from polygon)

text: I'm a newbie to OpenLayers (and JavaScript for that matter) and hoping you can offer some advice. I have been referring to examples http://openlayers.org/dev/examples/wfs-filter.html and http://osgeo- org.1560.x6.nabble.com/Select-WFS-Features-using-intersect-of-polygon- td4467419.html but have not had any luck getting any features returned from a spatial filter (a user-drawn polygon). The server response is as follows, and I'm not sure what it means:               <ows:Exception exceptionCode="NoApplicableCode">         <ows:ExceptionText>             org.xml.sax.SAXParseException; lineNumber: 1; columnNumber: 294;             Element or attribute do not match QName production: QName::=(NCName':')?NCName.              Element or attribute do not match QName production: QName::=(NCName':')?NCName.         </ows:ExceptionText>     </ows:Exception>      Anyway, my code is below. Please let me know if there is something I am missing.               var wfsPIF = new OpenLayers.Layer.Vector("PIFs", {         strategies: [new OpenLayers.Strategy.BBOX()],         protocol: new OpenLayers.Protocol.WFS({             version: "1.1.0",             url: "http://new-raptor.weatherfusion.com/geoserver/wf_data/wfs",              featureType: "erie_pif_wfs",             featurePrefix:"",             geometryName: "geom"         }),         projection: "EPSG:4326",         styleMap: new OpenLayers.StyleMap({            pointRadius: 2,            fillColor: "red"         })     });      map.addLayer(wfsPIF);               var drawControl, polygonLayer;          polygonLayer = new OpenLayers.Layer.Vector("Query Polygon");     map.addLayer(polygonLayer);     drawControl = new OpenLayers.Control.DrawFeature(polygonLayer,OpenLayers.Handler.Polygon);     map.addControl(drawControl);     drawControl.activate();          polygonLayer.events.on({         featureadded: function(event) {             var geometry = event.feature.geometry;             wfsPIF.filter = new OpenLayers.Filter.Spatial({                 type: OpenLayers.Filter.Spatial.WITHIN,                 value: event.feature.geometry             });             wfsPIF.refresh({force: true});             return false;         }     });


title: Values in attribute table not showing up

text: I have values in an attribute table I want to label my layer with. But when I label it, a totally different set of values show up. In the attribute table, my values are between 0 and 200. But the labels have values as extreme as 300,000. What's going on? Thanks!


title: Problem with performance when loading PostGIS layer using pyQGIS

text: I'm having a question on loading a PostGIS layer into the project from my QGIS plugin. The code is like this:               uri1 = qgis.core.QgsDataSourceURI()      uri1.setConnection(str(host), str(port), str(database), str(user), str(password) )     uri1.setDataSource(str(schema), "layer", "geom")     qgis.utils.iface.addVectorLayer(uri1.uri(),"layer", "postgres")      It actually works fine but the performance is causing me lots of troubles. It takes up to 10 minutes, since there are about a 100,000,000 entries in the table. When I load the layer manually using the PostGIS button on QGIS, there is a checkbox in the settings which says: "Geschätzte Tabellenmetadaten nutzen". In english this would be: "Use estimated table metadata". When checked, the layer loads in realtime. My question is: Is there a way to load a PostGIS layer into QGIS using the estimated table metadata, and thus increase performance significantly?


title: Why does XYToLine_management give ExecuteError ERROR 999999 The table was not found?

text: I've written this script to import a csv file using ArcGIS' flavor of python.. ArcPy.. It seems to run fine up until the XYtoLine tool runs and then it fails.. This is the error message...               Traceback (most recent call last):     File "C:\temp\GIS_projects\Sample_tripData\test.py", line 26, in <module>     arcpy.XYToLine_management(input_table,fc,"TSLONGDD","TSLATDD","TELONGDD","TELATDD","GEODESIC","TRIPSEQNUM")     File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\management.py", line 2913, in XYToLine     raise e     ExecuteError: ERROR 999999: Error executing function.     The table was not found. [tripLines]     Failed to execute (XYToLine).      this is the code I am running...                   # Import system modules     import arcinfo     import arcpy     import sys     if arcpy.CheckProduct("ArcInfo") == "Available":         print "ArcInfo available, processing script..."         print "importing modules and set environment...."              #overwrite pre-existing files         arcpy.env.overwriteOutput = True              #delete old feature class         if  arcpy.Exists("c:\\temp\\gis_projects\\sample_tripdata\\test.gdb\\tripLines"):             arcpy.Delete_management ("c:\\temp\\gis_projects\\sample_tripdata\\test.gdb\\tripLines")              #set workspace         arcpy.env.workspace = "c:\\temp\\gis_projects\\sample_tripdata\\test.gdb"         fc = "tripLines"              # Set local variables         print "setting data input and output variables..."         input_table = "C:\\temp\\GIS_projects\\Sample_tripData\\A6_TEMP_ATL_TRIP_HAULx.csv"              #XY To Line         print "running XY to Line Geoprocessing tool..."         arcpy.XYToLine_management(input_table,fc,"TSLONGDD","TSLATDD","TELONGDD","TELATDD","GEODESIC","TRIPSEQNUM")              # Process to delete identical line features...         # Set the field upon which the identicals are found         fields = "Shape;TRIPSEQNUM;Shape_Length"              # Set the XY tolerance within which to identical records to be deleted         xy_tol = "0.0002 Miles"              # Set the Z tolerance to default         z_tol = "0"              print "Deleteing identical features..."         # Execute Delete Identical          arcpy.DeleteIdentical_management(fc, fields, xy_tol, z_tol)              print "Dropping leading character..."              cursor = arcpy.UpdateCursor(fc)         field = "TRIPSEQNUM"         row = cursor.next()         while row:           # set value of field to only those values after the zero value location           # I want to strip off the first character within the value of the field           row.setValue(field, row.getValue(field)[0:])           cursor.updateRow(row)           row = cursor.next()              del fc, input_table, row, cursor, out_lines, input_table, fields, xy_tol, z_tol, field              print "Processing Complete!"          else:         print "ArcInfo license not available"         sys.exit("ArcInfo license not available")      I've been over this a bunch of times and I still keep getting this error.. anyone have any ideas as to what is going on? thanks in advance


title: vb.net arcobjects remove a buffer

text: I have a function which can select and zoom to a feature, at the same time I will create a buffer around the selected feature so that further more function provided only for the selected feature. The problem is when I clear the selection, the buffer outline is still there until another feature is selected. The buffer outline makes the feature inside looks different as other features in the same layer as shown in the following picture (in buffer outline is in red): ![enter image description here](http://i.stack.imgur.com/0WrcZ.jpg) I wonder is there a solution to change the buffer outline colour? Or is there a way to remove the buffer? I am using arcmap 10.2 and vb.net arcobjects.


title: How to show a specific geometry on the map and hide the others?

text: I have a FeatureLayer, i used the method `IFeatureLayer.Search()` and the interface `IQueryFilter` to get a specific Feature, and i want to show only this Feature on the Map, how can i do that ?


title: QGIS software Layer Attribute Name Length Limited to 10 characters

text: In QGIS software, the "vector layer" has "attribute" but the name length is limited to only 10 Characters. Is there any method to name an "attribute" to a longer name? (e.g. attribute name " **building name** " become " **building n** " in QGIS layer attribute table) ![enter image description here](http://i.stack.imgur.com/4adPm.png) My layer define dialog window is as below ![enter image description here](http://i.stack.imgur.com/ZgCHP.png)


title: Put color for raster data(aspect.tif) on QGIS1.8.0

text: I create aspect.tif from DEM file. I could put the values(0 to 360)from contrast enhancement. I could set color from color map(freakout or Pseudocolor). In here, the data is shown by only 4 colors. I could not change colors and add color type. First of all, I could not select colormap tab. My ideal color setting is 8 colors with 8 direction(0 to 45, 45 to 90, ....315 to 360) Please tell me how to change color.


title: Delete points on a polygon using QGIS

text: I'm new at this, a novice at best. I'm using QGIS 1.7.3. I simply want to select a large number of points on a single polygon and delete them. Is this possible?


title: ArcGIS 10 Python split string

text: I am trying to use Python to split a string field into three new fields at the separator " - ". The field is in the format "a - b - c". So I initially had this in the field calculator for newFieldA:               !myField!.partition(" - ")[0]      This gave me the value a in newFieldA. Similarly, I calculated newFieldB and newFieldC. After a more experienced Python user took a look at it, this was the final result:               ' ' if !myField! == 'NULL' else !myField!.split(' - ')[0]      The test for NULL was put in because there were some NULL values throwing errors. Each of the three fields were calculated using this code, just the index was changed (0, 1 or 2). My problem now is, after changing spelling errors in myField, I tried to split the field again, using the exact same code which worked previously. It now fails to run, giving the error of "wrong field name or unbalanced quotation marks". I am definitely using the correct field, and the quotation marks are not unbalanced. Also, I literally copied and pasted this from the code I used previously, so I don't see why it would be failing now. I feel like I'm missing something really simple. Does anyone have any ideas? **edit** I've placed the following in the codeblock               define splitField(value):     ' ' if value Is None     else:       value.split(' - ')[0]      and in the expression window               splitField(!myField!)      I am getting a syntax error on line 1, I have not been using Python long enough to see what it could possibly be, as that looks right to me. Does anyone have a suggestion to make the code better? **edit** I'm working with crime data, so in the field SUBACTION I have "Policing - Theft - Fraud". I then created three fields MANDATE, CATEGORY, and OFFENCE. The end result should be "Policing" in MANDATE, "Theft" in CATEGORY and "Fraud" in OFFENCE. I changed the codeblock to               def splitField(value):       ' ' if !SUBACTION! == 'NULL'           else:             value.split(' - ')[0]      It is now giving a syntax error in line 2.


title: Khmer Unicode script in ArcGIS 10

text: I am volunteering for the wonderful Cambodian urban NGO Sahmakum Teang Tnaut as a GIS and data analysis advisor for 10 months. My Khmer counterparts have found that ArcGIS does not support Khmer Unicode for producing Khmer script in maps, such as for labels and titles. Does anyone have any experience with getting Khmer Unicode to work? I, and my colleagues have it installed on my computer so can use it in other programs, just not ArcGIS it seems. At the moment we have to use any of the Limon typefaces, which are much more cumbersome to write in as it takes up to three characters per Khmer character.


title: osm2po "show option"

text: The "show expansion" option in the web interface, how do i invoke that from the command line. I am trying something like this...with "debugExpansion=true"... http://:8888/Osm2poService?cmd=fr&source=486477&target=0&findShortestPath=false&ignoreRestrictions=false&ignoreOneWays=false&routerClassId=0&heuristicFactor=0.0&maxCost=1.0&debugExpansion=true&key=Key&value=Value&format=geojson When I run this, I get a "ERR.1 No route found" return value. I want to get the convex hull that shows up in the web interface.


title: reusing cached tiles with leaflet, mbtiles and mbtiles-php

text: I am using Leaflet and infostream's fork of mbtiles-php to serve MBTiles stored on my server. I'm using a Leaflet layer group (and a layer control) to serve two MBTiles files. When I toggle between the two layers using the layer control at the same zoom level, my local cache is rebuilt rather than reusing the existing tiles. In the Leaflet layer group example, toggling between the two layers reuses the locally cached tiles, although I assume that the example tileset is not in MBTiles format. In the examples that I've seen through MapBox (using MBTiles but with the ModestMaps library), the locally cached tiles are also being reused. My information on reuse of caching come from clearing the cache (using Chrome) and then examining the developer tools network tab to see the tiles being built (or not built). Am I missing some customization of mbtiles-php that would allow reuse of tiles? I've seen reference to mbutils, but I'd prefer to keep the tiles in MBTiles format.


title: Join GPS points obtained from .GPX file in Python Console

text: Having loaded the .gpx file into a new layer from python console, I'm unable to access the gps points in the attributes table. Essentially, i'm unable to join these points to form a polygon (buildings) or a line segement (roads) depending on the number of gps points in a gpx file. It's possible to achieve the same in using the QGIS GUI, but I need to perform the same from the python console (to be implemented in a plugin). Any help or suggestions would be much welcomed. Thanks


title: Split/Erase/Differentiate polygon by polygon within a same layer

text: Is it possible to split a polygon with a polygon within the same layer in QGIS ? As in the image, i want the Polygon 2 to be split & the pointed area has to be retained as a separate polygon. Yes..i can use Split tool. But the problem is....if the polygon is so large & if it has so many nodes, we have to snap the split line with the nodes properly, rite ? Or is there any options available to trace the whole split part line...instead of snapping with each & every node. (i don't want to duplicate the layer & use geoprocessing tools like Clip/Difference/etc.. & then remove the duplicate geometry. Coz it brings lot of confusion due to the spatial nature of geological features we are dealing with) Grateful for any help.. ![enter image description here](http://i.stack.imgur.com/bvGek.jpg)


title: Why is this formatted string "n. def."?

text: I'd like to show the user the current map scale. I update the current map scale in a textblock everytime the map size/scale has changed. This is my code, to show the map scale (in C#):               private void Map_SizeChanged(object sender, SizeChangedEventArgs e){         MapScaleTextBlock.Text = String.Format("Scale {0:N}", Map.Scale); //Note: It has to be formatted     }      and here is the xaml code:               <TextBlock Name="MapScaleTextBlock"></TextBlock>      The output is always this: > Scale n. def. Whats wrong? I searched already but I found nothing helpful. I'm using ESRI API for WPF and C# in Visual Studio 2010.


title: D3 gdal Unable to open datasource problem

text: I am following the tutorial for d3.js here http://bost.ocks.org/mike/map/ I've download gdal using               brew install gdal      after which I checked using               which ogr2ogr      it returned               /usr/local/bin/ogr2ogr      however when I run this command               ogr2ogr \       -f GeoJSON \       -where "ADM0_A3 IN ('GBR', 'IRL')" \       subunits.json \       ne_10m_admin_0_map_subunits.shp      I get               FAILURE:     Unable to open datasource `ne_10m_admin_0_map_subunits.shp' with the following drivers.      so I tried                  /usr/local/bin/ogr2ogr \ -f GeoJSON \ -where "adm0_a3 IN ('GBR', 'IRL')" \ subunits.json \ ne_10m_admin_0_map_subunits.shp      it returned               Unable to open datasource `GeoJSON' with the following drivers.                     -> ESRI Shapefile       -> MapInfo File       -> UK .NTF       -> SDTS       -> TIGER       -> S57       -> DGN       -> VRT       -> REC       -> Memory       -> BNA       -> CSV       -> GML       -> GPX       -> KML       -> GeoJSON       -> GMT       -> SQLite       -> PCIDSK       -> XPlane       -> AVCBin       -> AVCE00       -> DXF       -> Geoconcept       -> GeoRSS       -> GPSTrackMaker       -> VFK       -> PGDump       -> OSM       -> GPSBabel       -> SUA       -> OpenAir       -> PDS       -> WFS       -> HTF       -> AeronavFAA       -> EDIGEO       -> GFT       -> SVG       -> CouchDB       -> Idrisi       -> ARCGEN       -> SEGUKOOA       -> SEGY       -> XLS       -> ODS       -> XLSX       -> ElasticSearch       -> PDF      which I don't understand as I downloaded both gdal and topojson, and geojson is obviously in the list of driver? Any help would be appreicated Edit: I agree it is the exact duplicate of the question listed, however if you closely examined the answer to that question, you will find that the answer did not solve the question but rather the issue goes away mysteriously                If you're in the same position I was at the time of asking this, I can only advise you to keep re-installing all the stuff you've already installed several times, reboot occasionally and hope.      I reboot and reinstall gdal and topojson couple of times, so I believe this question is in essence a repeat of an unsolved question I also tried the comment               ogrinfo ne_10m_admin_0_map_subunits.shp      and it returned               Unable to open datasource `ne_10m_admin_0_map_subunits.shp' with the following drivers.


title: markers disappear on map

text: I do not understand why the markers disappear, zoom 6, only in this area map. Do you have an idea? Map


title: Mask image based on no-data pixels in another image with different projection - GDAL or Python

text: I've got a land cover image which covers the whole globe, and a satellite image which covers a very small part of the globe. These two images have different projections, co-ordinate systems and pixel sizes. Also, the satellite image is not fully 'square' within the image, it is on a slant, and there are significant no-data areas around the edge of the actual image. I want to extract the pixels in the land cover image which intersect the satellite image. All of this needs to be completely automatic - either via a Python script or a GDAL command-line. In the past when I've had a polygon outline of the image I've used a `gdalwarp` command like               gdalwarp -overwrite -of GTiff -cutline polygon.shp -crop_to_cutline landcover.tif output.tif      But I can't find a good automatic way to produce a polygon of the footprint of my satellite image, given than I don't want the footprint of the _whole_ image, but just the part of the image that isn't No Data. So it doesn't look like that method will work. Theoretically I should be able to do a simple mask, but I can't work out how to do this when my images are of different resolutions. Does anyone have any ideas how I should go about doing this?


title: I had an error in Gdalwarp and i don't know the reason

text: I was trying to warp a raster using a mask and shift it before clipping but when i exectute the command bellow in Qgis/PyQt using QProcess:               gdalwarp -s_srs "+proj=merc +a=6378137 +b=6356752 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs" -t_srs "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs" -dstnodata -999 -q -cutline C:\Users\Mohamed.fkili\.qgis\statsRectangle2.shp -crop_to_cutline -of GTiff C:/Users/Mohamed.fkili/.qgis//0.0.3//data/rasters/UsableWaterBClipped.tif C:/Users/Mohamed.fkili/.qgis//0.0.3//data/rasters/UsableWaterBClipped.tif.tmp      i got the error bellow:               ERROR 1: Translating source or target SRS failed:     "+proj=merc +a=6378137 +b=6356752 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"      but in the dos commande ligne i have no problem and my generated raster is perfect. anyone can help identifying the problem?


title: No Popup for bottom layers, only top layer works

text: I have an OpenLayers Map. On this map, I have three overlay layers. The first 2 layers come from one GeoServer Layer/data source, with two alternate styles. The 3rd layer is a separate GeoServer Layer / data source. For example, I have a GeoSever Layer called 'Roads'. For my 'Roads' Layer, I have a 'Highway' style to show features with a highway attribute, and a 'Street' style to show features with a street attribute, but both are from my 'Roads' layer. The third layer is a called 'Lights', and is a separate set of features. The issue is that in OpenLayers, If I have only the 'Roads/Highway' Layer visible, I can successfully do a WMSGetFeatureInfo for my popup. If I only have the 'Roads/Street' Layer visible, I can do a WMSGetFeatureInfo for my popup. But, when I have both the 'Roads/Highway' and the 'Roads/Street' layers simultaneously visible, I can only get a response from the top most drawn layer (painters model). The other layer beneath the resonsive layer, doesn't give me featureInfo. Note: when I have either the 'Roads/Highway' or the 'Roads/Street' layer visible with 'Lights' layer, then I can get a featureInfo response from both the 'Lights' layer and whichever 'Roads' layer I have on. One further note which I am not sure is relevant, is that I have my GeoServer styles connected to a Workspace.


title: How to query annotations from a GDB file in .Net

text: I have annotations in a GDB file that I need to extract to load into a non- ESRI spatial store. It's my understanding that there isn't any tool to export annotations to a shape file. I need to copy points, text, and style info. I have not worked with annotations before. My prior experience with ESRI products involves writing a shape file parser in C#.Net and writing .Net code that interacts with the C API for ArcSDE to query ArcSDE. It's my understanding that ArcSDE SE_annotation methods don't really return very much detail about annotations, so I would be open to working with ArcObjects as well. I've been reviewing that API, however, and I'm not clear on how I would query annotations. If anyone has experience with this sort of thing and can suggest ideas or point me to relevant documentation I'd be very grateful.


title: Issue on Setting Coordinates for SVG Map

text: Can you please help me to figure out how to set coordinates for BC Canada , based on following snippet. This code is setting coordinates for French as the `43.15710, 8.17199` is the lower right point of French Box boundary               getCoords : function (lat, lon) {         if (lat < 43.15710 && lon > 8.17199) {                  var xfactor = 43.64246;             var xoffset = 181.34520;             var x = (lon * xfactor) + xoffset;                  var yfactor = -65.77758;             var yoffset = 3346.37839;             var y = (lat * yfactor) + yoffset;             }              else {             var xfactor = 45.48385;             var xoffset = 220.22005;             var x = (lon * xfactor) + xoffset;                  var yfactor = -65.97284;             var yoffset = 3371.10748;             var y = (lat * yfactor) + yoffset;             }                  return {x : x, y : y};             },      now in case of BC Canada which the lower right point is `49.044281, -114.031232` how can I calculate the xfactor, yfactor, xoffset and yoffset? Thanks


title: Spatial Analysis - Identifying Regions with Common Boundaries

text: I have the shape file for DC PSA's as defined in 2012 (see below). I would like to be able to identify in semi-automatic fashion which PSA's are neighbors - that is, which PSA's are nearest neighbors to one another (for example, 201 and 401 are nearest neighbors). With this information, one could define a sort of spatial distance metric between PSA $i$ and PSA $j$, call it $d(i,j)$, populating a distance matrix, call it $D$. Then $d(i,j)=0$ is pure autocorrelation (it is the diagonal of $D$). If $d(i,j)=1$, you have nearest neighbor spatial autocorrelation. One can then handle $d(i,j)>1$ separately. With certain physical systems this seems like a simplistic but natural way of thinking about distance. Does anyone know of a simple way to be able to calculate $d(i,j)$ as define here, without doing it manually? ![DC PSA Map 2012](http://www.neighborhoodinfodc.org/psa12/psa_2012.jpg)


title: Loading a layer into the same position as an existing layer in QGIS

text: I understand that it is not possible to directly reload a shapefile layer in QGIS using the `reload` method in Python, since this method is not implemented for shapefile data providers. Thus, the workaround is to unload the layer and then reload the layer. However, if one does this naively, the new layer will show up at the top of the map legend, rather than at the position that the original layer occupied. Is there a straightforward way to save the position of the old layer and then reload a new layer at that same position?


title: Is it possible to place a point at a coordinate in QGIS?

text: Is it possible to place a point at a coordinate in qgis? I have a map, using irish grid reference and wish to place a point at an exact 6 digit coordinate. Is it possible to do it. Previous to this i have placed the point roughly and then moved it whilst zoomed to the exact location, surely there is an easier way.


title: Mercator projection that shrinks toward pole

text: I've noticed that lots of map viewer (especially those nautical chart viewer) will shrink a mercator projected chart while panning towards poles. Is this some compensation for mercator distortion? Is this a dedicated projection? How to implement that?


title: Can geoserver limit layer access to clients by credentials?

text: I have 2 clients which don't want to host geoserver on their systems and use my server. Since I can encode their js , is there a way to limit the workspace for each client do other can't call it using wms ?


title: Testing software (like QTP ,Selenium) for ESRI platform (Desktop and Web based)

text: Does anybody know any testing software/tools for ArcGIS platform (version 10.1). Specifically : Customized ArcGIS Desktop application & Web application. I researched on Google but am not able to find any specific software. Few experts told me to try manual testing so I'm a little bit confused. Any help/suggestion will be great !! Thanks in advance : )


title: Why are layers not overlapping to each other inspite of having same CRS?

text: I am new to QGIS and working with QGIS 1.7.3 version, I have my all layers in WGS 84 CRS but even though they are not overlapping to each other, What should I do ?


title: Drop down menu not populating for ArcGIS JavaScript API Application

text: 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;     }


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

text: 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.!


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

text: _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._


title: convert GPS to x,y using reference node

text: 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 !!


title: ComException -2147467259 Arcgis 10 VB.net +Geoprocessing

text: 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


title: Geoserver proxy extention

text: 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 !!


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

text: 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


title: ArcGIS 10, Converting from date to number

text: 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?


title: Find out tight corners on roads

text: ## 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.


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

text: 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


title: How to calculate the total ascent along a path

text: 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!


title: ArcMap Commands .NET code equivalent

text: 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


title: WFS - MapServer - ASP.NET

text: 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!


title: Where to get 2010 Census Block data?

text: 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?


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

text: 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.


title: QGIS renders outlines where no should be

text: 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)


title: Issue with Oracle Connection on Windows 7 machine

text: 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)


title: Point input for polyline not acceptable?

text: 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)


title: 50cm imagery - aerial or satellite?

text: 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?


title: Standard Style for Infrastructure

text: 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


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

text: 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.


title: How to use Geokettle to load shapefiles into oracle spatial

text: 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


title: Join Fields from Multiple Tables to One Table in Modelbuilder

text: 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!


title: How to parse ArcGIS Server generated output XML?

text: 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.


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

text: 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.


title: How to unproject WKT to WKT in .NET?

text: 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?


title: Import shapefile is not working in pyshp?

text: 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


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

text: 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.


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

text: 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?


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

text: 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 ?


title: ArcObjects Java Calculate Line Length

text: 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);


title: Setting Band Statistics on Raster Image (GeoTIFF)

text: 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!


title: How to compare schemas from two File geodatabases?

text: 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.


title: How to configure various client applications to display AVL data?

text: 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?


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

text: 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)


title: Steps to polygonize raster KML image

text: 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?


title: Simple Features and ArcGIS Domains

text: 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.


title: Is the ArcGIS Server JavaScript API free to use?

text: 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?


title: Converting GRIB to Geotiff with gdal_translate in Python

text: 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


title: How to estimate the average width of a sidewalk?

text: 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;


title: OpenLayers: Point icons disappearing

text: 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


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

text: 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 :-)?


title: List Joined tables

text: 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


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

text: 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.


title: FME - Convert DWG to JPG

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


title: Address or reverse geocode a point shapefile

text: 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?


title: Can I convert a kmz file to kml?

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


title: Save user specified settings on flash drive

text: 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!


title: Does QGIS DB Manager plugin lock PostGIS transactions?

text: 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.


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

text: 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.


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

text: 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.


title: How can I plot GPS coordinates on Google Map

text: 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)?


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

text: 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?


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

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


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

text: 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.