Using WMTS Stream in OpenLayers

I am trying to write some OpenLayers javascript that can display a WMTS stream from a publicly available REST service located at: http://www.mchenrycountygis.org/arcgis/rest/services/2012aerials/MapServer The problem I'm having is that while I am able to see OpenLayers generating Ajax requests to the endpoints, the tiles it is asking for do not exist. By playing around with the settings, I'm able to change these URLs somewhat, but I can't get the right combination of settings to actually generate valid URLs. My current progress is documented in this jsfiddle: http://jsfiddle.net/v4WuW/2/ Code is also posted below Any help would be appreciated! var map = new OpenLayers.Map({ div: "map", // QUESTION 1: This likely needs to be EPSG:102671 to match the WMTS stream. Why // does using the correct value completely break things? projection: "EPSG:900913" }); // matrixIds was extracted from capabilities.xml via OpenLayers.Format.WMTSCapabilities() var matrixIds = [{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"0","scaleDenominator":545977.1577380951,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":3,"matrixHeight":2},{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"1","scaleDenominator":272988.5788690476,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":4,"matrixHeight":3},{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"2","scaleDenominator":136493.81696428568,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":6,"matrixHeight":5},{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"3","scaleDenominator":68246.43601190475,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":10,"matrixHeight":9},{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"4","scaleDenominator":34122.74553571428,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":19,"matrixHeight":18},{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"5","scaleDenominator":17060.900297619046,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":38,"matrixHeight":35},{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"6","scaleDenominator":8529.977678571428,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":75,"matrixHeight":68},{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"7","scaleDenominator":4264.516369047618,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":148,"matrixHeight":135},{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"8","scaleDenominator":2131.7857142857138,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":295,"matrixHeight":271},{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"9","scaleDenominator":1065.8928571428569,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":589,"matrixHeight":540},{"supportedCRS":"urn:ogc:def:crs:EPSG::102671","identifier":"10","scaleDenominator":566.9642857142857,"topLeftCorner":{"lon":-17463800,"lat":19496100},"tileWidth":256,"tileHeight":256,"matrixWidth":1107,"matrixHeight":1013}];               var wmts = new OpenLayers.Layer.WMTS({         url: 'http://www.mchenrycountygis.org/arcgis/rest/services/2012aerials/MapServer/tile/{TileMatrix}/{TileRow}/{TileCol}',         isBaseLayer: true,         requestEncoding: 'REST',         layer: '2012aerials',         matrixSet: "EPSG:102671",         matrixIds : matrixIds,         style: 'default',         /* QUESTION 2         Ajax requests fire for incorrect tile URLs such as:       www.mchenrycountygis.org/arcgis/rest/services/2012aerials/MapServer/tile/0/27/29              How do I fix this so we are requesting tiles that actually exist, such as:       www.mchenrycountygis.org/arcgis/rest/services/2012aerials/MapServer/tile/0/135/142              Adding/Removing tileOrigin & tileFullExtent doesn't seem to help...         */         tileOrigin: new OpenLayers.LonLat(-17463800,19496100),         tileFullExtent: new OpenLayers.Bounds(882499.9493310973, 1990000.0923467577, 1029999.9493310973, 2125000.0923467577),         layers: ''     });                          map.addLayers([wmts]);     map.zoomToExtent([807168.3566418022, 2101119.5440438855, 1105331.5420203924, 2131750.0923467577]); // "initial" extent value     //map.zoomToExtent([882499.9493310973, 1990000.0923467577, 1029999.9493310973, 2125000.0923467577]); // "full" extent value     //map.setCenter(new OpenLayers.LonLat(??, ??), 13);


Composite geocoder in ArcGIS 10.2.2 gives different results on second run

I am using a composite geocoder created in ArcGIS 10.2.2 with parcel points, address points, 3 street databases, city and zip code points. When I send the same address database to the geocoder a second time, I end up with a handful of records that get different results. The majority of these are addresses that contain a highway or county road address, but there are some that contain just a street or avenue address. I thought I could use the rematch option to see what the score results were for the various options, but after looking at some of these, I could see that it was re-evaluating on the fly and I was not seeing the original results. Has anyone else noticed this or any idea why this could be occurring? It has happened on small (1000 records) files and large 250,000 ones.


what is the minimum to install of QGIS for it to work?

This may sound like a daft question, but I'm looking to understand what the minimum parts of the install procedure I need to install to get QGIS to work (i.e. can it be installed without GRASS etc?) as I'm looking for the easiest way to roll this product out to a number of users preferably via a packaged installer, so any suggestions/hints would be greatly appreciated. Thanks


How to catch invalid expression errors in Python?

I have a script that creates a new feature using `arcpy.GetParameterAsText`, but I'm having trouble creating an error message for incorrect user input. The input for a parcel # should be something like '1500300000029', but what if 'sdfjfh' is entered? As it stands no error is thrown and a empty feature is created. I've also tried inside PythonWin by assigning a bad entry for 'input' and it runs without error.               import arcpy     import sys          arcpy.env.overwriteOutput = True     try:         Parcels = "..path to data"     arcpy.env.workspace = "C:/Data_new/Temp/Default.gdb"          #input = arcpy.GetParameterAsText(0)         input = "dfdsgdg"                           expression = "PARCEL_ID = " + "'" + str(input) + "'"     arcpy.MakeFeatureLayer_management(Parcels, "parcels_lyr")     arcpy.SelectLayerByAttribute_management("parcels_lyr", "NEW_SELECTION", expression)             arcpy.CopyFeatures_management("parcels_lyr", "output")     except:         print arcpy.GetMessages()      I've also tried               except StandardError, ErrDesc:         arcpy.AddWarning(input)         arcpy.AddWarning(ErrDesc)          except:         arcpy.AddWarning(input)         arcpy.AddWarning(ErrDesc)      What is the appropriate way to handle this?


How do I add a field to a table using the QGIS graphic modeler?

I was wondering if anyone could help as I'm struggling to get anything out of the QGIS graphic modeler. For info I have no programming skills whatsoever. I need to be able to add a new field to an attribute table and update that field with a value that the user will enter. Is there anyway to do that with the graphic modeler? Also is there a comprehensive guide to the graphic modeler? Many thanks


What to do when raster cell size is larger than zonal vector feature in order to collect statistics?

I am trying to extract the mean cell values for a raster dataset based on a survey grid polygon feature using ArcINFO 10.1. The cell size of the raster dataset is about 4 times the size of the survey grid. I am working with the FCLIM climate trend analysis data for equatorial Africa (ftp://www.earlywarning.usgs.gov/fews/pubs/mapping_decadal_variations.pdf) I want the mean value of raster cells/pixels within each polygon grid cell: ![FCLIM MAMJ rain trend raster overlaid by survey grid](http://i.stack.imgur.com/6bdmO.png) I have around 3000 grid cells, using 'zonal statistics as table' tool I get only around 600 results evenly spread about the survey grid. Obviously not the result I am seeking. When I use the 'zonal statistics' tool the result I get the following raster: ![output raster of zonal statistics tool using FCLIM MAMJ rain trend as input raster and survey grid as feature zone data](http://i.stack.imgur.com/2wzEe.png) 2 things have gone wrong with this output (that may not be immediately evident in the image): the cell size of the output raster is still about 4X larger than the feature zone, survey grid polygons; and secondly the raster cells are not coincident with the feature zone survey grid polygons and indeed no longer coincident with the original raster dataset from which it was derived! My question is: How can I calculate the mean value of raster cells/pixels residing within the smaller survey polygon grids? Is it a matter of using a different method of breaking up the raster cells or is there another method I should try?


Getting started with gdal/ogr

I found some batch code using ogr2ogr in order to convert and reproject vector layers. The difficulty for me is to know what to install and which gdal commandline should I use, and how to configure it... The fact is I'm quite lost with all the packages that gdal offers and I don't know :   1. which package am I supposed to install : OSGeo4W, FWTools, MS4W,...?? What is the difference between them ? Can I use any of them ?   2. What is the procedure to set PATH and environment variables ? I work with windows 7.   3. I have already QGIS and GRASS installed on my computer. Is an installation of gdal could damage anything ?


Programmatically export a Qgis project

I have a QGis project that I programmatically modify using the Python xml.etree library. Now, I would like to know if it is possible to add few lines to my script in order to export the modified project to a new jpg image. The python script should be run OUTSIDE the QGis python console. Thanks!


Running Maperitive on a Mac

Really basic question, but I need details and help running Maperitive from my mac. I've tried launching it from the terminal and tried running it through WineBottler. If anyone has experience with this, please help!


How do I calculate Manhattan distance with PostGIS?

I am using the ST_Distance function to calculate the distance between two geometries (a train station and a building) Since I know that all buildings and all train stations are in Chicago, which has an excellent/complete street grid, I would like to use Manhattan (or taxicab) distance. The generic formula for this is the difference in X plus the difference in Y, so Abs(X1-X2)+Abs(Y1-Y2). What PostgreSQL query would make this work?


how to display map center coordinate using openlayers

As far as I know that openlayers can display our mouse coordinate position. However, I do not know how to display map center coordinate? Can anyone help me? thanks for your helping guys!


40k+ points from Shapefile, clustering them and populating on Google Maps/Webapp

I just started a new Python/Django webapp job related to GIS (and i'm totally new to GIS). Here is my task and what i have found so far, and i am stuck as to finding out what's the best way or the global standard way of doing things.   * We are building a webapp where a user can upload a shapefile, and take the records and points from the file then populate them on a map.     1. Upload a shapefile     2. Parsed the shapefile using pyshp library - this seems OK     3. I'm able to read the shapes (coordinates) and the records (value for each field on each shape)     4. Now say for instance i've got about 40,000+ coordinates, each of them with field records in it     5. I'll need to populate those 40k points onto Google Maps - those points should be represented with a marker, and the field records can be represented with a heatmap (depending on which field the user choose to filter it as)     6. Certainly loading 40k points onto Google Maps wouldn't make sense as it will slow down the client side of things. The question is, what can i do to make this fast or make this work? Thinking of clustering nearby coordinates in a grid, and not only that, each coordinate has a field record and we have to somehow cluster them(?) so that the clustered coordinate has an accurate representation of the values in it. A few options that i have looked at are :-   1. Use http(:)//google-maps-utility-library-v3.googlecode.com/svn/trunk/markerclusterer/docs/reference.html, this would speed things up.. but the server will still be sending 40k points, which may not work at a point where we need to load more than 40k points for each different section on the map   2. Was looking at how we can achieve http(:)//www.gdal.org/grid_tutorial.html by using http(:)//www.gdal.org/gdal_grid.html - I've not been able to get the .tiff generation to work yet   3. Also to achieve grid data, i found http(:)//docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.griddata.html - though i don't think this actually does it   4. I just happen to learn about GeoDjango and https(:)//github.com/biodiv/anycluster, and am looking into it to see if this should be something i should use   5. and i'm still continuing to research I don't know if there is any global standard, it looks like a common problem that plenty of people may have come across. I am probably unaware of the lingo terms to search for. A pointer to the right direction, would be appreciated.


How to add constraints in PostGIS?

I have table that containing road network. Table structure is               CREATE TABLE final_all_india_road_network     (       gid serial NOT NULL,       "name" character varying(254),       code character varying(254),       lane character varying(254),       width integer,       "type" character varying(254),       subclass character varying(254),       surface character varying(254),       oneway character varying(254),       bridge integer,       maxspeed integer,       the_geom geometry,       CONSTRAINT final_all_india_road_network_pkey PRIMARY KEY (gid),       CONSTRAINT enforce_dims_the_geom CHECK (st_ndims(the_geom) = 2),       CONSTRAINT enforce_srid_the_geom CHECK (st_srid(the_geom) = 4030)     )      I want to add check geometry constraints to this table. For that I used following Query               "ALTER TABLE final_all_india_road_network       ADD CONSTRAINT enforce_geotype_the_geom CHECK (geometrytype(the_geom) = 'MULTILINESTRING'::text OR the_geom IS NULL);"      But while running this Query following error occured "ERROR: check constraint "enforce_geotype_the_geom" is violated by some row"


Label styling applied in Geoserver are not being applied to WMS layer in Openlayers

I am trying to label a layer of line vectors using SLD in Geoserver. It is simple enough, and it's pretty much copied from the SLD cookbook: http://docs.geoserver.org/latest/en/user/styling/sld-cookbook/lines.html The thing that's getting me is that when I view the Openlayers preview in Geoserver, everything seems to render correctly. Each road is only labeled once and it's not extremely busy looking. However, when I import the layer via WMS to my Openlayers map, it doesn't seem to recognize this line:               <VendorOption name="repeat">400</VendorOption>       Which is supposed to keep it from repeating within 400 pixels. I've also tried:               <VendorOption name="group">yes</VendorOption>       Which should only display a single label on a map per unique text. Any idea what is going on here? I need the labels to look like the Geoserver preview. Map displayed in openlayers after importing layer via WMS: ![WMS Served Map](http://oi57.tinypic.com/ajm5no.jpg) Map displayed in Geoserver Openlayers preview: ![Openlayers preview](http://oi59.tinypic.com/pt1xs.jpg) `<?xml version="1.0" encoding="UTF-8"?> <StyledLayerDescriptor version="1.0.0" xsi:schemaLocation="http://www.opengis.net/sld StyledLayerDescriptor.xsd" xmlns="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"> <NamedLayer> <Name>Roads</Name> <UserStyle> <FeatureTypeStyle> <Rule> <Title>Roads</Title> <LineSymbolizer> <Stroke> <CssParameter name="stroke"> <ogc:Literal>#000000</ogc:Literal> </CssParameter> <CssParameter name="stroke- width"> <ogc:Literal>1</ogc:Literal> </CssParameter> </Stroke> </LineSymbolizer> <TextSymbolizer> <Label> <ogc:PropertyName>name</ogc:PropertyName> </Label> <LabelPlacement> <LinePlacement> <PerpendicularOffset> 10 </PerpendicularOffset> </LinePlacement> </LabelPlacement> <Fill> <CssParameter name="fill">#000000</CssParameter> </Fill> <VendorOption name="repeat">400</VendorOption> </TextSymbolizer> </Rule> </FeatureTypeStyle> </UserStyle> </NamedLayer> </StyledLayerDescriptor>`


Urban heat island

Since there are lots of products in MODIs for Urban heat Island can anyone suggest me which is the best one?and why?I need to start working on Urban heat Island of some city in Asia but before that I need to know which product is best for that.


Euclidean distance tool issue,not obtaining output

I'm trying to perform an exercise but I'm not having the output that I was supposed to have. Which says this Procedure:   1. Open ArcMap and load the Spatial Analyst extension: Customize → Extensions → Spatial Analyst.   2. Add the raster named Sea and the point feature class onept.shp. Enhance the symbol size for the point file.   3. Make the Outputs.gdb your default geodatabase.   4. Enter ArcMap’s Environment Settings via the Geoprocessing → Environments entry. Open the Processing Extent section and set the Extent to Sea. Next, open the Raster Analysis section and set the Cell Size to the same as the raster Sea. Close the environment entry.   5. Open the Viewshed tool and set the Input Raster to Sea and the Input Point to Onepoint. Execute the tool and inspect the result that should look like Figure 2.(I get exactly what it says in the fig 2 which it hasn't been posted, although with different colours) The problem comes in this step. 6.Next, the euclidean distances and directions from the point and outwards should be calculated. Open the Euclidean Distance tool and set the Input raster or feature source data to Onepoint. The output distance raster should be named Distance The Output cell size should already be set to 1852 which is the resolution defined in the Environment. Finally, type Direction for the Output direction raster and start the tool. Inspect the two output rasters that should look like Figures 3 and 4. (They do not look at all as it should) ![Output](http://i.stack.imgur.com/ypT9e.jpg) Figure3 ![Fig3](http://i.stack.imgur.com/77wiA.png) Fig4 ![Fig4](http://i.stack.imgur.com/FveG1.png) I have followed the instructions carefully step by step so I don't know why is returning me this. The only think that I could've set different is the step 4. I set the extent to "Same as layer sea" instead of sea because I didn't have "sea option". Also I enter in geoprocessing cell size as"the same as layer sea" instead of "the same as raster sea". I have no idea what could be the problem, any insights? Editon after Surya comment: This is the result that I have obtained after dropping the layer Direction and Distance. ![enter image description here](http://i.stack.imgur.com/QHkVR.jpg)


Openlayers zooming and panning restrictions don't work

Consider the following example (I tried to keep it as simple as possible).

                  As the `OSM` map is in the Spherical Mercator Projection, the bounds `(-1500000, 4000000, 5600000, 11000000)` should cover most of Europe. I set these bounds both as the `maxExtent` and the `restrictedExtent` of the map. Furthermore, I set `maxScale` to `60 M`. After launching this example, the map really centers to Europe and restricts panning to other parts of the World. But, it absolutely disregards the `maxScale` attribute and allows zooming far away. Now, I try to remove the `restrictedExtent`. I should be able to pan the map freely. This is true. However, now the `zoomToMaxExtent()` doesn't work, although the `maxExtent` is still set! How come neither `maxScale` nor `zoomToMaxExtent()` would work in this example? Am I just gravely misunderstanding the OpenLayers documentation?


Is it necessary/advisable to uninstall QGIS 1.8 first before installing 2.0?

Who has installed QGIS 2.0 on Mac OSX 10.6.8 (Snow Leopard)? Is it necessary/advisable to uninstall QGIS 1.8 first? Have you seen any problems importing QGIS 1.8 projects to 2.0?


GeoEXT using WFS-T draw on point Layer

Hello this is my first post. I'm new with openlayers and geoext and i'm facing a problem while trying to add a new record in a point layer via wfs-t the code is below (it is from http://workshops.opengeo.org/geoext/wfs/editing.html):               function reconfigure(store, url) {                 var fields = [], columns = [], geometryName, geometryType;                 // regular expression to detect the geometry column                 var geomRegex = /gml:(Point|Line|Polygon|Surface|Geometry).*/;                     var types = {                     // mapping of xml schema data types to Ext JS data types                     "xsd:int": "int",                     "xsd:short": "int",                     "xsd:long": "int",                     "xsd:string": "string",                     "xsd:dateTime": "string",                     "xsd:double": "float",                     "xsd:decimal": "float",                     // mapping of geometry types                          **//Here i want to put something for drawing in point layer**                     "Line": "Path",                     "Surface": "Polygon"                 };      I've tried to put "Poin": "Point", but it didn't work! Has anybody an idea about what i have to put there so as to work with point layer too(it already works with line and polygon layers)?? Another point of the code which may cause the problem is the save function:               {         text: "Save",         handler: function() {             app.featureGrid.store.proxy.protocol.commit(                 vectorLayer.features, {                     callback: function() {                         var layers = app.mapPanel.map.layers;                         for (var i=layers.length-1; i>=0; --i) {                             layers[i].redraw(true);                         }                         app.featureGrid.store.reload();                     }             });      Which of the two may cause the problem? Thank's a lot!


how to not print scalebar and how to print the zoom bar that shows on map

I am using javascript api for ArcGIS and need to print the zoom bar on map, but don't want the scale information to be printed. How can I achieve this?


QGIS crash Ubuntu 14.04

Encountering a crash when requesting QGIS 2.2 to access data on desktop machine. Running Ubuntu 14.04 LTS.               $ qgis     Warning: loading of qt translation failed [/usr/share/qt4/translations/qt_en_AU]     qgis.bin(32381)/KSharedDataCache KSharedDataCache::Private::mapSharedMemory: Failed to establish shared memory mapping, will fallback to private memory -- memory usage will increase      kbuildsycoca4 running...     kbuildsycoca4: ERROR creating database '/home/lyndon/.kde/cache-Gerhard-ThinkPad-W530/ksycoca4'! Insufficient permissions in target directory.     Warning: ERROR: Running KSycoca failed.     QGIS died on signal 11ptrace: Operation not permitted.     No thread selected     No stack.     gdb returned 0     Aborted (core dumped)'       Any request, such as opening a file, browsing, creating a new project etc. will cause this crash. I am new to QGIS and will eagerly post additional information, if provided instructions. EDIT: output of:               sudo ls -l /home/lyndon/.kde/cache-Gerhard-ThinkPad-W530/               total 10320      -rw-r--r-- 1 root root 10547304 Apr 25 2013 icon-cache.kcache –       EDIT: running QGIS desktop as root does allow me to open files! But it also pops up with an error message immediately after the program runs (ie after the initial loading/welcome graphic dissapears): **The layer /usr/share/applications/%F is not a valid layer and can not be added to the map** Any idea how I can make QGIS operate from GUI (unity) with root privileges without having to open it via terminal and adding the sudo prefix? I have attempted the solution provided below by edm, but the open window hangs and I recieve errors pertaining to klauncher: **Could not start process Cannot talk to klauncher: Not connected to D-Bus server.** http://askubuntu.com/questions/118822/how-to-launch-application-as-root-from- unity-launcher               $ sudo qgis %F               Warning: loading of qt translation failed [/usr/share/qt4/translations/qt_en_AU]     Error: "/tmp/kde-lyndon" is owned by uid 1000 instead of uid 0.     Error: "/tmp/ksocket-lyndon" is owned by uid 1000 instead of uid 0.     Error: "/tmp/kde-lyndon" is owned by uid 1000 instead of uid 0.     kdeinit4: Shutting down running client.     Connecting to deprecated signal QDBusConnectionInterface::serviceOwnerChanged(QString,QString,QString)     Error: "/tmp/ksocket-lyndon" is owned by uid 1000 instead of uid 0.     Error: "/tmp/kde-lyndon" is owned by uid 1000 instead of uid 0.      QDBusConnection: session D-Bus connection created before QCoreApplication. Application may misbehave.               QDBusConnection: session D-Bus connection created before QCoreApplication. Application may misbehave.     kbuildsycoca4 running...     kbuildsycoca4(8139) VFolderMenu::loadDoc: Parse error in  "/home/lyndon/.config/menus/applications-merged/xdg-desktop-menu-dummy.menu" , line  1 , col  1 :  "unexpected end of file"      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::init: The desktop entry file  "/usr/share/applications/evolution-data-server-uoa.desktop"  has Type= "Application"  but no Exec line           kbuildsycoca4(8139) KBuildServiceFactory::createEntry: Invalid Service :  "/usr/share/applications/evolution-data-server-uoa.desktop"      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S0" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S1" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S2" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S3" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S4" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S5" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S6" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S7" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S8" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S9" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::parseActions: The desktop file "/usr/share/applications/Launchpadlaunchpadnet.desktop" references the action "S10" but doesn't define it      kbuildsycoca4(8139)/kdecore (services) KServicePrivate::init: The desktop entry file  "/usr/share/applications/unity-scope-gdrive.desktop"  has Type= "Application"  but no Exec line           kbuildsycoca4(8139) KBuildServiceFactory::createEntry: Invalid Service :  "/usr/share/applications/unity-scope-gdrive.desktop"       kbuildsycoca4(8139)/kdecore (services) KServicePrivate::init: The desktop entry file "/usr/share/applications/updf.desktop" has Type= " Application" instead of "Application" or "Service"               kbuildsycoca4(8139) KBuildServiceFactory::createEntry: Invalid Service :  "/usr/share/applications/updf.desktop"       kbuildsycoca4(8139)/kdecore (services) KServicePrivate::init: The desktop entry file "/usr/share/applications/unity-lens-photos.desktop" has Type= "Application" but no Exec line               kbuildsycoca4(8139) KBuildServiceFactory::createEntry: Invalid Service :  "/usr/share/applications/unity-lens-photos.desktop"      kbuildsycoca4(8139) kdemain: Emitting notifyDatabaseChanged ("servicetypes", "xdgdata-mime", "services", "apps")     Error: "/tmp/kde-lyndon" is owned by uid 1000 instead of uid 0.     Warning: Object::connect: No such signal org::freedesktop::UPower::DeviceAdded(QDBusObjectPath)     Warning: Object::connect: No such signal org::freedesktop::UPower::DeviceRemoved(QDBusObjectPath)     Error: "/tmp/ksocket-lyndon" is owned by uid 1000 instead of uid 0.     Object::connect: No such signal org::freedesktop::UPower::DeviceAdded(QDBusObjectPath)     Object::connect: No such signal org::freedesktop::UPower::DeviceRemoved(QDBusObjectPath)


How do I skip rows with searchcursor based on an interger value

For example, if I have 100 records, and I want to sample 10% of the records. How do I tell the search cursor in arc10 to read the value of the 9th record, or every 10. I tried to do row[0 + count] and updated a count variable += 9 but it said I can't slice a row. Any suggestions, thanks.


QGIS 2.2 Changeing txt data in theme table not working

A week ago I upgrade from QGIS 2.0.1 to 2.2 and try to edit some simple point themes I made recently. To my surprise, QGIS 2.2 randomly - not allways, does not allow me to change text I type into text cells - like changeing one letter to another and such. My shp file goes into editing mode, I can type text into cells but after exiting cell, old text is displayed. Insertig data with Field calculator does not help, either. I have no problems with same shp files in QGIS 2.0.1. When I try to edit shp again latter, everything works fine ... very frustrating. Does anyone expirience this behavior of QGIS 2.2 as well?


How to Optimize Openstreetmap Tile Rendering?

I am currently trying to render Water (3,198,333 records) and detailed world borders (processed_p) for the world. I am successful in rendering both layers, but I am not happy with performance. My current setup is as follows:   * Amazon EC2 instance (m2.2xlarge)    * 34.2 GB of memory   * 13 EC2 Compute Units (4 virtual cores with 3.25 EC2 Compute Units each)   * 850 GB of instance storage   * 64-bit platform   * I/O Performance: High   * Mapnik, 4 threads   * Postgresql settings based on these benchmarks   * GiST index in place for all of my layers **Some of my Ideas to Increase Performance** :   * Simplifying my features   * Don't render water tiles - when loading with openlayers can set load failure to blue color   * Split world water features into various sections and explicitly tell mapnik its bounds. For example, if I split water features from North America, when Mapnik is generating tiles for Australia it should not query North America to see if any of its features are in the tiles.    * My second idea is a bit more crazy. Since Mapnik performs a spatial query for each tile in order to find features to process, I can perform this query beforehand for all of my layers. I would assume that this would save some processing power.  **My Question:** What else can I do to optimize tile rendering?


How Do I Capture when the Cluster Rendering is Completed

I have a graphics layer that draws clusters using the sample custom cluster from ESRI. http://help.arcgis.com/en/webapi/silverlight/samples/start.htm#CustomClusterer We have a Busy Indicator that we pop up during long calculations to let the user know that something is happening. I would like to pop up this indicator during the clustering, and have it go away once the clusters have done drawing on the screen. I can turn it on, but I am not able to capture a "done rendering" event to turn off the indicator. I've tried capturing the LayoutUpdated event on the Map, but it gets sent several times and am not able to determine which one is the graphics layers saying its done rendering. I found the OnClusteringCompleted event of the GraphicsClusterer and tried to override it by adding it to the custom clusterer but the code never gets there. Any suggestions on how to accomplish what I want?


srtm30 image: some are grayscale some are black & white

I downloaded these DEM files, http://www.dgadv.com/srtm30/ srtm30 globe. I open the DEM files in Multispec. Some of them are black& white, some are automatically grayscale. How do i get them all in grayscale? (attaching image to explain) Thank you!!! ![enter image description here](http://i.stack.imgur.com/a2rUd.jpg) ![enter image description here](http://i.stack.imgur.com/FZYdc.jpg)


Coloring thousands of polygons so each is distinct from its neighbors?

I am using ArcGIS 10, and I have layers of the 2010 Census blocks, block groups, and tracts. There are thousands of polygons in these layers. I want to color them so that each is distinct from its neighbors, but all I can figure out is how to give each one a distinct color. The color of each polygon does not really matter; I just want to easily distinguish each from its neighbor by color, rather than with a border. Or is that something you only see in an elementary school map of the 50 United States?


Trouble with adding WFS layer in ArcMap

I add the following WFS address to ArcCatalog just fine and the different feature classes show up in ArcCatalog. However, when I attempt to drag or Add a feature into the TOC of ArcMap, it gives me a "This dataset appears to be invalid" error and the feature gets added to the TOC but is empty with no feature records. Could I be missing a step here? I'm using ArcGIS 10.1. Here is the WFS address:               http://ogi.state.ok.us/geoserver/wfs?service=WFS&version=1.0.0&request=GetFeature?


How do I draw an EllipticArc (or any other Segment) in a Dynamic Layer?

I am trying to draw an EllipticArc (or any other Segment for that matter) in my Dynamic Layer using the ArcObjects C++ API. For something that seems like it should be straightforward, I am having a very hard time figuring it out. So far I have taken these steps:   1. Create the EllipticArc itself. I have confirmed that I have created one and that its geometry is correct.   2. Cast the EllipticArc to a PointCollection so it may be drawn by DynamicDislay. The problem seems to be that the PointCollection I create contains 0 points, and I am not sure how to solve this problem. Here is what I have so far:               PointCollection4Ptr  finalform(CLSID_Polyline);     IConstructEllipticArcPtr arcCon(CLSID_EllipticArc);          /*Code to construct arc omitted, it was rather lengthy.       I have confirmed that the arc was correctly constructed        and that the object exists*/          IEllipticArcPtr arc = arcCon;             //Cast to a PointCollection so we can draw it     IPointCollection4Ptr    pointcol = arc;               finalform->AddPointCollection(pointcol);              std::cout<<"before pointcount\n";     long pc;     finalform->get_PointCount(&pc);     std::cout<<pc<<std::endl;     std::cout<<"after pointcount\n";          pDynamicDisplay->DrawPolyline(finalform);      The code outputs that the pointcount is 0. Am I completely missing something here? I feel like this should be straightforward, but I have had no luck drawing any Segment on my Dynamic Layer. How is this supposed to be done? Can it be done? Is there a solution for drawing Segments on non-dynamic layers that I will have to use? (I hope not...) Any input is appreciated!


arcpy selectByLocation

EDIT: OK, added two lines to the code below including a projectAs(). When I start my mxd and run this code it selects a polygon. However, if I then pan to a new location - still within the bounds of the target FC - it doesn't work. Is there some update command I need to run after panning? It's as though the dataframe extent is lost after panning. ==== I'm trying to selectByLocation using a pointGeometry. The code below creates a pointGeometry from the center of the current dataframe extent. It works if the dataframe and the target feature class (polygon features) are in the same coordinate system but not when they're in different systems. My usage model is that the target will usually be in a continental projected system or geographic coordinates while the dataframe will be in NAD 1983 UTM 10N to 20N.               import arcpy          mxd = arcpy.mapping.MapDocument("CURRENT")     df = arcpy.mapping.ListDataFrames(mxd)[0]          arcpy.env.outputCoordinateSystem = df.spatialReference     **sr = arcpy.Describe('targetFC').spatialReference**          myPt = arcpy.Point((df.extent.XMax + df.extent.XMin)/2.0, (df.extent.YMax + df.extent.YMin)/2.0)     myPtGeometry = arcpy.PointGeometry(myPt, df.spatialReference)     **myPtGeometry.projectAs(sr)**     arcpy.SelectLayerByLocation_management('targetFC', 'INTERSECT', myPtGeometry)          del df, mxd      Thank you, Scott


How to zoom to a selected feature with pyQGIS

I want to create a function that selects a feature and zooms to it (similar in QGIS). Therefore there is following function:               QgsMapLayerRegistry.instance().addMapLayer(self.vlayer)          def zoomTo(self):         layer = self.vlayer         atable = self.ui.table              selectList=[]         for i in atable.selectionModel().selectedRows():             ID = atable.item(i.row(),0).text()             selectList.append(int(ID))              layer.setSelectedFeatures(selectList)      The selected features are highlighted on the map. But I have no idea how to make a "zoom" to the selected features or some kind of focus them in the middle of the map.


How to work with NREL data in QGIS?

I would like to use QGIS to display solar data from NREL. The file is supposed to provide solar radiation data of Colorado. However, I have some trouble to select this data in QGIS: In QGIS, I select "Add Vector Layer" and select the shp file. QGIS displays the raster for the information (10 km by 10 km data). All squares have the same color. But I can't see the actual data. I was always under the impression that the shapefile would also contain the geographical data. Is that correct? **Do I need to load any other projection file in QGIS to display the solar radiation per latitude and longitude?** **How can I display the NREL data correctly with QGIS?** Thank you for your answer!


Using a CAD template in ArcGIS?

Is there a simple way to create a "CAD template" for map outputs from ArcGIS as in the image below: ![enter image description here](http://i.stack.imgur.com/LqTaq.gif) I'd like to create a standard template and then populate some simple 'fields' with map information (title, scale, date).


How to create a point table from a remote csv file?

The current approach that I am using to load a csv file into a PostGIS enabled database on a server is the following:               -- create table     CREATE TABLE points(      id int,       lat double precision,      lng double precision      );          -- copy data from csv with headers     COPY points (id, lat, lng)     FROM '/home/mydata.csv'     WITH          DELIMITER AS ','          CSV HEADER ;          -- Add point geometry column to table     SELECT AddGeometryColumn('public', 'points', 'geom', 900913, 'POINT', 2);          -- Populate column with point geometries     update points set geom = ST_SetSRID(ST_MakePoint(lng, lat), 900913);          -- Create a spatial index on points     CREATE INDEX idx_points ON points USING GIST (geom);      Although this creates a table with point geometries, there appear to be some restrictions when using the `COPY` command as you need to be a superuser. Also, I can't figure out how to upload a csv file from my laptop to the server using this approach (I need to transfer the file using ftp and then run the script). Is there another strategy that I should be considering for generating a table with point geometries from a csv? I want to perform this entire process in PostGIS/SQL.


Project layer with custom projection

My ultimate goal is to know UTM and UTM Northing and Easting. The data can be anywhere in North america. I calculated UTM using arcpy tool (arcpy.CalculateUTMZone_cartography) and I tried to truncate the value (e.g.'17N') to Project_management as part of Output coordinate string but so far no luck. Please advise. Thanks The code is here:               arcpy.CalculateUTMZone_cartography(outPointSHP, "UTM")          pointUTM= arcpy.SearchCursor(outPointSHP)     for row in pointUTM:         UTM_var= str(row.getValue('UTM'))[30:32]     del pointUTM          ok = UTM_var     arcpy.AddMessage(ok)          out_coordinate_system = os.path.join(connectionPath+'/', r"projections/NAD1983/NAD 1983 UTM Zone "+ok+"N.prj")               arcpy.Project_management(outPointSHP, projPointSHP, out_coordinate_system)


What is the relative popularity of ESRI GDB files vs GDB over RDBMS for custom web apps?

Our organization is interested in knowing the frequency with which people who are using ESRI geodatabases choose GDB files over GDB in relational databases for serving web applications. Opinions and anecdotal evidence are welcome, as well as any hard data on this topic that anyone might have. Thanks!


How to check if wms layer is queryable in openlayers?

I'm running a map with lots of different wms layers published by geoserver. Some of them are queryable, some of them are not. How can I tell that with openlayers api? I'm requesting getfeatureinfo and I need to request only the queryable layers, layers with queryable param set to false always return "Layer is not queryable" when requested - and that's what I want to hide from users. I looked into OpenLayers.Map.layers and found no "queryable" property. UPDATE (to the first commnent): This works fine as long as I use it with localhost:8080/geoserver/wms?service=wms&version=1.1.1&request=GetCapabilities. When I tried to get capabilities from http://www.edpp.cz:8080/geoserver/wms?service=wms&version=1.1.1&request=GetCapabilities, i got empty array of layers as a result. I'm trying it with the following code:               function getCapabilities(offline) {         var url = '';         if (offline) {             url = 'localhost:8080/geoserver/wms?service=wms&version=1.1.1&request=GetCapabilities';         }         else {             url = 'http://www.edpp.cz:8080/geoserver/wms?service=wms&version=1.1.1&request=GetCapabilities'         }         var layers = []             ,result             ,request = OpenLayers.Request.GET({                 url: url                 ,success: function(response) {                     var format = new OpenLayers.Format.XML()                         ,xml = format.read(response.responseText)                         ,text = format.write(xml)                         ,CAPformat = new OpenLayers.Format.WMSCapabilities()                         ,cap = CAPformat.read(xml);                         if (cap.capability) {                             for (var len = cap.capability.layers.length-1; len >= 0; len -= 1) {                                 if (cap.capability.layers[len].queryable) {                                     layers.push(cap.capability.layers[len].name);                                 }                             }                         }                 }             });         return layers; //ARRAY     }


ArcGIS 10.1 Lightweight Runtime

Does anyone know where I can find more information about the new Lightweight Runtime that will be included with ArcGIS 10.1? This was discussed at the recent DevSummit and ESRI Federal User Conference. I am looking for any general information about the product, but more specifically I am trying to find out if there will be any sort of licensing fees to use it. More information can be found here: http://esriaustralia.wordpress.com/2011/02/17/the-new-arcgis-runtime/ Thanks.


URL " http:// <mygisServer>/ArcGIS/rest/service/GPTask/MapServer/jobs/<jobid>" not found in the ArcGIS Server?

I am Using ArcGIS Server 10.03 SP3. I created a geoprocessing service in the server which returns a result map service. I successfully created the gp service and map service with same name. I need to display the result in my application using ArcGIS Javascript API v3.3. The logic works like the sample application http://developers.arcgis.com/en/javascript/samples/gp_resultmapservice/ I tested the gp service in the REST Web interface and successfully executed.The result map is obtaining after completion of the task. But when I use it in the JS API front end , the fire bug tool shows success execution of the gp task , but when it executes               map.addLayer(gpLayer);      the response shows               dojo.io.script.jsonp_dojoIoScript9._jsonpCallback({"currentVersion":10.03,"error":{"code":400,"message":"Unable to complete  operation.","details":["Invalid URL"]}});      When I tried the url http:// /ArcGIS/rest/service/GPTask/MapServer/jobs/ The page not found shows. When I changed the code with               gp.getResultImageLayer(jobinfo.jobId, "output", imageParameters, function (gpLayer) {      map.addLayer(gpLayer)     });      I got the result . Can anybody help to find the exact reason? I need to include the service via ArcGISDynamicMapServiceLayer Here is the My sample code that Did not work               function applyGPTask() {     var gpServiceUrl = "http://localhost/ArcGIS/rest/services/GPTasks/GPServer/GPTask";     gp = new esri.tasks.Geoprocessor(gpServiceUrl);     var params = {         param1: "SUBDIVID=202",         param2:"value2",         param3: "value3"           };           gp.setOutputSpatialReference({ wkid: 4326 });       gp.submitJob(params, gpJobComplete, gpJobStatus, gpJobFailed);      }          function gpJobComplete(jobinfo) {         var mapurl = "http://localhost/ArcGIS/rest/services/GPTasks/MapServer/jobs/" + jobinfo.jobId;     var gpLayer = new esri.layers.ArcGISDynamicMapServiceLayer(mapurl, {         "id": "gpoutput",         "opacity": 0.7     });         map.addLayer(gpLayer);      }


Howto setup and use a remote GIS workstation (QGIS,GRASS,PostGIS)?

I'm using QGIS (and sometimes GRASS and PostGIS) on my 2009 MacBookPro for my thesis project.   I was wondering (during an endless batch generalizing in QGIS) if I can set up a server so I can work from the library/university using the power of a faster home desktop workstation.   So far I understood:   * I can access a remote PostGIS server via QGIS   * I can run GRASS via SSH, even with GUI   * GDAL has PostGIS driver, but it's only for querying and inserting data, the rest of the work is done by the machine that issues the command (or am I wrong?)   * Basic editing of a PostGIS can be done via QGIS, and I can issue SQL commands to make some more advanced editing What I would like to figure out:   * Can I use QGIS plugins remotely? is VNC the only way to accomplish that?   * is there some documentation for PostGIS+QGIS other than creating and visualizing maps (I found only these types of tutorial)? which tools/plugins can I use, and which of them can be run remotely? **UPDATE:**   you all confirmed my thoughts about QGIS, if I want to use the processing power of the server I have to control it via VNC/RDP/NX.   as for my second question, let me rephrase it:   If I use postGIS+QGIS and I want all the processing be done on the server, do I have to (learn and) use only SQL queries with postGIS functions?


Displaying Highrise Delivery Routes

I have a geocoded address table for High Rise buildings in New York. If i'm trying to map delivery routes to a number of these locations, what would be the best and most clear way to display the routes between different floors?


How to make new polygon/feature from point which is NOT centroid?

I've spent all morning searching help sites and digging through the FAQs and worksheets but have failed miserably. I am an ArcGIS novice so I hope this is not too much of a "newbie" question but I'm trying to work out how to make a new polygon/feature from a point which is NOT a centroid?               ----                ----     | . | <--not this   |   |     ----       this-->  .---      I have a series of points where the Easting and Northing listed represents the lower left hand corner of a 200m square. I know how to make a circular buffer which I can convert to a square with the Feature Envelope to Polygon tool but this isnt really what I want. The data is of 3g strengths which is meant to be in grid squares. I cannot upload this file this due to NDA. Any ideas?


How to import ArcGIS 10.1 File geodatabase to MS SQL Server 2012?

I would like to import my File Geodatabase from ArcGIS 10.1 into MS SQL Server 2012. I have never done it before, so I expect to do the import the same way as import data from excel table. Problem is when I open MS SQL Server Import wizard and as a data source I choose Esri geodatabse OLE DB provider - an error occured: _ClassFactory cannot get required class. (Microsoft OLE DB Service Components)_ Please does anybody know what I´m doing wrong, or is there some easiest way, for example with ArcSDE?


QGIS to solve projection mystery in Postgresql

I downloaded some data from local agency. Per QGIS 2.2, I defined the projection to 'EPSG:3587 - NAD83(NSRS2007) / Michigan Central' in postgresql (maybe ArcMap). Data is in Postgresql 9.3. I am putting this data on geoserver and probably thought to store it in the web mercator projection 3587. You can see in the image the black square near Australia is where the features land but should be where the blue swingley circle is. I brought in OpenStreetMap from OpenLayers plugin in QGIS and then brought in the borked layer. The geometry scale doesn't seem to be off...or said differently...it's features seem to measure a reasonable distance. If I had to guess I would say that I projected it correctly but may have defined it incorrectly in postgresql. I could run through the steps again with new data (I did this a few months ago and subsequently got my laptop stolen so no trace of anything I did but my virtual servers humming along). I would like to solve this problem as well as become well versed with QGIS. Any ideas? I was thinking of changing projection in Postgresql but I don't know what I'm doing. Thanks! ![enter image description here](http://i.stack.imgur.com/ZxMXm.png)


Python Hyperlinking Script in Arcmap 10.2

I have the below python script which recursively walks through a directory and finds, then opens the associated parcel pdf file. How would I modify this to not be recursive but to just look in the directory of "\corp.qr.com.au\app\QRGIS\QRData\Leases\Lease_Plans\Active" I can't seem to get it to work. Thanks               import webbrowser, os     def OpenLink([File_Name]):       for r,d,f in os.walk(r"\\corp.qr.com.au\app\QRGIS\QRData\Leases\Lease_Plans\"):         for image in f:           if image.startswith(str([File_Name])):             webbrowser.open(os.path.join(r,image))       return


Is it possible to set the project CRS via a Python command?

Is it possible to set the project CRS of an project in QGIS via a Python command? In my case the user should open QGIS, start my written plugin and this plugin will load spatial data which has the epsg code 31467 (Gauss Krueger Zone 3). Unfortunately QGIS default CRS is WGS84 an so I get into trouble with the units and the scale. My first approach was to use this code:                   my_crs = core.QgsCoordinateReferenceSystem(31467, core.QgsCoordinateReferenceSystem.EpsgCrsId)         self.iface.mapCanvas().mapRenderer().setDestinationCrs(my_crs)      This changes (after the data is already visualized) the CRS to Gauss Krueger Zone 3. At least the espg code in the lower right corner has changed. The units are still in degrees and the scale is totally wrong. 1:5000 in Gauss Krüger is 1:528822376 in WGS 84 (in my case in QGIS). After setting the project CRS manually in the properties the project crs is really GK3 and the scale and the units are right. So my question is, can I set the project CRS in another way than I have done it? I want to avoid the manual way. * * * I solved the problem by adding the following code:               qgis.utils.iface.mapCanvas().setMapUnits(0)     qgis.utils.iface.mapCanvas().refresh()      This sets the units to meter and the scale is updated automatically. The 0 stands for meters, 1 is for foot, 2 is for degree and 3 for unknown.


Pros and Cons of virtualisation technology for database-server and datastorage

Yesterday we had a chat about performance and restorability and I realized how many good things a virtualisation environment can bring to me - but as I'm a little sceptic concerning performance I'm asking here Will a database server suffer from severe performance loss through virtualisation? I don't understand the technology to the last detail, but somehow it's a 'blackbox' more that needs to be processed through hardware. Also will disk access, with all the tricks postgis gives to us be granted (clustering, index etc.) - fragmented clustering is like no clustering! The big advantage is maintenance and scalability. In case of a severe hardware malfunction I can migration within minutes or even realtime to another physical machine. Who has experience and can point me to good websites or literature on this topic. I remember a few things from the last FOSSGIS and a few inhouse benchmarks on ESXi and native Servers and somehow I can't make up my mind whether it's good or not.


How to change Map Scale programmaticaly in C# WinForm Project

I have created form with **axMapControl** and **comboBox** with values of possible scales (e.g. 1:1 000 000, 1: 500 000, etc). How can I change current map scale programmaticaly by using value from selectedItem taken from Combobox?


Is there a way to tell geoserver to send cp1250 http header?

I'm running geoserver with some web pages in data_dir/www. These pages are encoded in cp1250, but when I load them using geoserver (localhost:8080/geoserver/www) they are encoded or they behave as they were encoded in utf-8. As a result i get tons of question marks instead of real letters. I am not able to convert pages into utf-8, so is there a way to tell geoserver not to send utf-8 http header and send cp1250 header instead?


Cause of selected [ArcSDE 10] feature dataset does not contain any feature classes which can participate in topology?

When trying to create a topology for a polygon feature class in ArcSDE 10. I right click the feature dataset that the feature class resides and try to create a new topology. I get an error saying, "The selected feature dataset does not contain any feature classes which can participate in a topology". The feature class I am targeting has a z value. Would this be the culprit?


Finding polygons without right angles

Now we are digitizing some buildings in specified area. The obligatory rule for this work - in most cases buildings should have right angles. We are using QGIS with CAD tools for this work. But sometimes we make mistakes and create polygons with an irregular shape. Does anybody know how can we find such polygons without right angles using open GIS or ArcGIS?


With OpenLayers addLayer(s) Isn't Displayed

I am trying to figure out why the addLayer here isn't displaying the boxes that are present in the example shown below:               function($) {              var box_extents = [             [-10, 50, 5, 60],             [-75, 41, -71, 44],             [-122.6, 37.6, -122.3, 37.9],             [10, 10, 20, 20]             ];         var map;              map = new OpenLayers.Map({             div: 'map',             projection: new OpenLayers.Projection('EPSG:900913'),             'displayProjection': new OpenLayers.Projection('EPSG:4326')         });                   function init() {             var ol_wms = new OpenLayers.Layer.WMS('OpenLayers WMS', 'http://vmap0.tiles.osgeo.org/wms/vmap0?', {                 layers: 'basic'             });                  var boxes = new OpenLayers.Layer.Vector('Boxes');             var i;                  for (i = 0; i < box_extents.length; i++) {                 ext = box_extents[i];                 bounds = OpenLayers.Bounds.fromArray(ext);                      box = new OpenLayers.Feature.Vector(bounds.toGeometry());                 boxes.addFeatures(box);             }                  map.addLayers([ol_wms, boxes]);             map.addControl(new OpenLayers.Control.LayerSwitcher());             var sf = new OpenLayers.Control.SelectFeature(boxes);             map.addControl(sf);             sf.activate();             map.zoomToMaxExtent();         }              init();     })(jQuery);      It works fine with addLayers in another example as displayed here:               (function($) {                  var map;                  function init() {                 map = new OpenLayers.Map('map');                      var layer = new OpenLayers.Layer.WMS( 'OpenLayers WMS',                         'http://vmap0.tiles.osgeo.org/wms/vmap0', {layers: 'basic'} );                 map.addLayer(layer);                      // allow testing of specific renderers via '?renderer=Canvas', etc                 var renderer = OpenLayers.Util.getParameters(window.location.href).renderer;                 renderer = (renderer) ? [renderer] : OpenLayers.Layer.Vector.prototype.renderers;                      var vectorLayer = new OpenLayers.Layer.Vector('Simple Geometry', {                     styleMap: new OpenLayers.StyleMap({                         'default': new OpenLayers.Style(null, {                             rules: [                                 new OpenLayers.Rule({                                     symbolizer: {                                         graphic: false,                                         label: 'Label for invisible point',                                         labelSelect: true,                                         fontStyle: 'italic'                                     },                                     filter: new OpenLayers.Filter.Comparison({                                         type: '==',                                         property: 'topic',                                         value: 'point_invisible'                                     })                                 }),                                 new OpenLayers.Rule({                                     symbolizer: {                                         stroke: true,                                         fill: true,                                         label: 'Polygon with stroke and fill defaults'                                     },                                     filter: new OpenLayers.Filter.Comparison({                                         type: '==',                                         property: 'topic',                                         value: 'polygon_defaults'                                     })                                 }),                                 new OpenLayers.Rule({                                     symbolizer: {                                         stroke: true,                                         fill: false,                                         label: 'Point without fill',                                         labelAlign: 'rb',                                         fontColor: '#ff0000',                                         fontOpacity: 0.4                                     },                                     filter: new OpenLayers.Filter.Comparison({                                         type: '==',                                         property: 'topic',                                         value: 'point_nofill'                                     })                                 })                             ]                         })                     }),                     renderers: renderer                 });                      // create a point feature                 var point = new OpenLayers.Geometry.Point(-111.04, 45.68);                 var pointFeature = new OpenLayers.Feature.Vector(point);                 pointFeature.attributes = {                     topic: 'point_invisible'                 };                      // create a polygon feature from a linear ring of points                 var pointList = [];                 for(var p=0; p<6; ++p) {                     var a = p * (2 * Math.PI) / 7;                     var r = Math.random(1) + 1;                     var newPoint = new OpenLayers.Geometry.Point(point.x + 5 + (r * Math.cos(a)),                                                                  point.y + 5 + (r * Math.sin(a)));                     pointList.push(newPoint);                 }                 pointList.push(pointList[0]);                      var linearRing = new OpenLayers.Geometry.LinearRing(pointList);                 var polygonFeature = new OpenLayers.Feature.Vector(                     new OpenLayers.Geometry.Polygon([linearRing]));                 polygonFeature.attributes = {                     topic: 'polygon_defaults'                 };                      multiFeature = new OpenLayers.Feature.Vector(                     new OpenLayers.Geometry.Collection([                         new OpenLayers.Geometry.LineString([                             new OpenLayers.Geometry.Point(-105,40),                             new OpenLayers.Geometry.Point(-95,45)                         ]),                         new OpenLayers.Geometry.Point(-105, 40)                     ]),                     {                         topic: 'point_nofill'                     });                      map.addLayer(vectorLayer);                 vectorLayer.drawFeature(multiFeature);                 map.setCenter(new OpenLayers.LonLat(point.x, point.y), 4);                 vectorLayer.addFeatures([pointFeature, polygonFeature, multiFeature]);                 var select = new OpenLayers.Control.SelectFeature(vectorLayer, {                     selectStyle: OpenLayers.Util.extend(                         {fill: true, stroke: true},                         OpenLayers.Feature.Vector.style['select'])                 });                 map.addControl(select);                 select.activate();             }              init();     })(jQuery);      They are different enough that it's difficult to compare, but the formatting is pretty much the same: http://dev.openlayers.org/docs/files/OpenLayers/Map- js.html#OpenLayers.Map.addLayer This is a GML layer is a vector layer subclass from what I've gathered: http://docs.openlayers.org/library/layers.html


How to remove overlapping buffer in QGIS?

I wanna know how to do this in QGIS: ![enter image description here](http://i.stack.imgur.com/NkQcG.gif) I've tried many ways, but still can't find the right solution.


OpenLayers dynamic rules defintion & clone function filter

I'm developing on a OpenLayers application which creates dynamic rules for a vector layer using the geostats library. I want to create these rules inside a simple loop through the number of classes. I'm pushing every rule to a array and add this array via stale.addRules(). In my opinion every rule has to be cloned before adding it to the array. The problem is, I am using a OpenLayers.Filter.Function in which the clone() function is not implemented yet (OpenLayers 2.12-rc4). How can I solve this? At the moment, only the last rule gets visualized. It seems like the other rules get lost, due to the reference value issue... My code looks like this:               for (var i = 0; i < numClasses; i++) {         filter_x = new OpenLayers.Filter.Function({                 evaluate: function(attributes) {                     [...]                 }             })         var rule_x = new OpenLayers.Rule({             filter: filter_x,             symbolizer: { fillColor: colors[i],                         fillOpacity: 0.5, strokeColor: "white"}         });         //var clone = filter_x.clone(); //returns null! clone is not implemented         rules.push(rule_x)          }     style.addRules(rules);      Thank you for your help!


How can I add data to a shp file attribute table from arcobjects?

I have a code that selects row data from a listview table in VB.net. I have declared that data as a variable and I can see it if I ask it to dispay a message box (see code below). Now I would like to take that data and place it into either an existing or new shapefile using ArcObjects but I am unclear on the process. Any hints or ideas would be helpful.               Private Sub ListView1_SelectedIndexChanged(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles ListView1.SelectedIndexChanged              'met en mémoire l'item selectionné         Dim i As Integer              If ListView1.SelectedItems.Count = 0 Then             Exit Sub         End If         'Met en mémoire les sous items de l'item séléctionné         Dim x As String = ""         Dim y As String = ""         Dim z As String = ""              'séléction de l'index de l'item         i = Val(ListView1.SelectedIndices.Item(0).ToString)         ' spécifie de prendre toute la rangée sélctionné         ListView1.FullRowSelect = True         ' Set the view to show details.         ListView1.View = View.Details              'Fait en sorte d'extraire les sous items 1, 2, et 3 (les 3 rangée "message, GeoLat, GeoLong"         x = (ListView1.Items(i).SubItems(1).Text)         y = (ListView1.Items(i).SubItems(2).Text)         z = (ListView1.Items(i).SubItems(3).Text)              'Concatenation des trois sous items en une ligne.         Dim a As String = ""         a = i & " " & x & " " & y & " " & z              MsgBox(Err.Number & "  " & Err.Description)         'Pour verifier si la routine fonctionne...         MessageBox.Show(a)          End Sub


COUNTY WebService API lookup from valid street address

I am looking for a solid WebService API that will allow me to submit an http request with a pre-validated US Street/City?state/Zip Address and return COUNTY in xml or json about 200 times per work day. I am aware that Google, MapQuest & Bing do this however the response will go into an internal database, my use does not fit into their free Terms Of Service and I can not afford their large fee for my very small business. It looks like TigerWeb may do it? Is the data good enough? What might a request look like? Many thanks in advance.


ModelBuilder - How to calculate Field using part of path name

![I attempted to include a model parameter for the workspace, rather than use the parse path tool \(show in the image\).](http://i.stack.imgur.com/7BXOM.png)I am trying to calculate a field using part of the feature class path name (specifically, the last three characters of the geodatabase name). Is this possible? I've tried using the parse path tool and including the value as an inline variable subsitution with the VB string right function. Perhaps a more complicated Python script is neccessary?


How to Invert Deleting multiple fields with delete tool in ArcGIS 10.1 SP1

I have a feature (polygon) with nearly 60 fields in it. Now I want to delete all fields except 20(or a list of) and they are not one after one, no they are somewhere in the droplist. But to select for a tons around 40 feature shape files would be a huge searching and clicking time...because he update for every new loaded feature his droplist and i can start searching again and again. Normally it must be go the opposite way...or? I am looking for a tool or script by "invers droping". From the web I didnt get a useful solution. maybe someone here? Thanks for hints...if I forget something to explain, i can draw a scetch...


Sorting points with QGIS

I have a very good idea of what I would like to do but I really do not know if it is possible with QGIS. I have points that each one corresponds to a group size and I would like to use GIS to sort this data. In fact, in order to conserve independence between my group sizes I would like to keep only one group size per day. Thus, is QGIS capable of analysing all my points and telling me how many duplicates I have within 100 meters radius for the same day? Then, if the software finds more than one point per day, can I run an application that will only keep one group size (randomly). I do not know if I am being very clear but I hope you understand my problem. ![enter image description here](http://i.stack.imgur.com/07Ik3.jpg)               Date    Day Month   Year    Week no Season  Time    EASTING NORTHING    DATUM   Groupsize     14/03/2006  14  March       2006    11  WINTER  09:20   599707  5733317 NAD83   1     02/04/2002  2   April       2002    14  SPRING  14:34   602500  5732500 NAD83   500     05/02/2002  5   February    2002    6   WINTER  09:00   600359  5734295 NAD83   250     06/02/2002  6   February    2002    6   WINTER  08:30   599000  5733600 NAD83   600     07/02/2002  7   February    2002    6   WINTER  18:00   599200  5733400 NAD83   650     11/02/2002  11  February    2002    7   WINTER  14:30   599850  5734150 NAD83   565     12/02/2002  12  February    2002    7   WINTER  17:10   599661  5733382 NAD83   29     13/02/2002  13  February    2002    7   WINTER  09:00   599800  5734200 NAD83   650     16/02/2002  16  February    2002    7   WINTER  10:30   599010  5733850 NAD83   100     16/02/2002  16  February    2002    7   WINTER  10:30   601034  5734740 NAD83   150     16/02/2002  16  February    2002    7   WINTER  10:30   600463  5734418 NAD83   75     16/02/2002  16  February    2002    7   WINTER  10:30   600628  5734783 NAD83   150     16/02/2002  16  February    2002    7   WINTER  10:30   600147  5734358 NAD83   100     20/02/2002  20  February    2002    8   WINTER  12:00   602500  5734200 NAD83   325     26/02/2002  26  February    2002    9   WINTER  18:05   599661  5733382 NAD83   155      Thank you,


problem r.walk region GRASS

I'm getting troubles to run r.walk. If I run r.walk in this region:               projection: 3 (Latitud - Longitud.)     zone:       0     datum:      wgs84     ellipsoid:  wgs84     north:      33:39:03N     south:      29:16:03N     west:       33:54:48E     east:       40:34:24E     nsres:      0:00:03     ewres:      0:00:03     rows:       5260     cols:       7992     cells:      42037920      the algorithm runs fine. But isn't the region I need analyze, when I change into this other one:               projection: 3 (Latitud - Longitud.)     zone:       0     datum:      wgs84     ellipsoid:  wgs84     north:      37:21:17.568N     south:      25:06:44.82N     west:       31:46:15.96E     east:       43:23:24.576E     nsres:      0:00:08.999949     ewres:      0:00:08.999272     rows:       4897     cols:       4648     cells:      22761256      The algorthym stucks in the Calculating cost step. It's a widest area, in exchange the number of cells has been reduce a lot. I'm running GRASS 6.4.3 on Ubuntu 12.04. I thougth this could be due to a lack of RAM???. Can anyone help me?? Thanks a lot.


What does this style of line on Topographic maps indicate

If I was reading the purpose of this stackexchange site right, I could also ask general map questions. Not necessarily regarding software. Anyways, I was wondering if someone could inform me what this blue line indicates on this topo map image. I understand the circular line with the ticks indicates a lower elevation, possible sinkhole-like, but I don't know what that blue line is. Also, does that line with the Xs through it normally indicate a fence? ![enter image description here](http://i.stack.imgur.com/AGGVS.png) If it matters at all, this is a map image from an ESRI powered site. EDIT: Here is a link http://www.kgis.org/KGISMaps/map.htm?map=topography&box=2540262:579354:2540837:579771


Rename a SpatialPolygon class object in R

I am completely new to GIS but am quite familar with R. I was able to dissolve some zip codes into a larger sales territory. For reference, my previous question is located here. However, because I will be doing this for 304 sales markets across the US, I will need to be able to give each a name in order to merge on my sales data after the fact. I am assuming this can be done by accessing the ID slot in the object. I am trying to access the `ID` slot. The following code does not work:               tmp2@polygons@Polygons@ID      Results in the following error:               Error: trying to get slot "Polygons" from an object of a basic class ("list") with no slots      Sorry for the long post, but here is my structure. Any help will be greatly appreciated!               An object of class "SpatialPolygons"     Slot "polygons":     [[1]]     An object of class "Polygons"     Slot "Polygons":     [[1]]     An object of class "Polygon"     Slot "labpt":     [1] -71.57027  41.85902          Slot "area":     [1] 0.1317919          Slot "hole":     [1] FALSE          Slot "ringDir":     [1] 1          Slot "coords":                    x        y       [1,] -71.50031 41.68146       [2,] -71.49977 41.67378       [3,] -71.49942 41.66667       [4,] -71.51520 41.66557       [5,] -71.52351 41.66233       [6,] -71.52396 41.66257       [7,] -71.52686 41.66166       [8,] -71.52693 41.66192       [9,] -71.53344 41.66194      [10,] -71.53336 41.66227      [11,] -71.53302 41.66798      [12,] -71.53346 41.66943      [13,] -71.53425 41.66901      [14,] -71.53506 41.66911      [15,] -71.53519 41.66959      [16,] -71.53382 41.67065      [17,] -71.53287 41.67143      [18,] -71.53369 41.67211      [19,] -71.53256 41.67314      [20,] -71.53284 41.67978      [21,] -71.53262 41.68029      [22,] -71.53201 41.68061      [23,] -71.53342 41.68218      [24,] -71.53287 41.68275      [25,] -71.53351 41.68434      [26,] -71.53293 41.68442      [27,] -71.53246 41.68836      [28,] -71.53329 41.68855      [29,] -71.53262 41.69093      [30,] -71.53384 41.69152      [31,] -71.53514 41.69312      [32,] -71.53530 41.69482      [33,] -71.53604 41.69489      [34,] -71.53604 41.69548      [35,] -71.53312 41.69562      [36,] -71.53327 41.69661      [37,] -71.53295 41.69729      [38,] -71.53399 41.69721      [39,] -71.53412 41.69808      [40,] -71.53219 41.69826      [41,] -71.53242 41.70002      [42,] -71.53274 41.69999      [43,] -71.53303 41.70017      [44,] -71.53316 41.70017      [45,] -71.53319 41.70034      [46,] -71.53441 41.70025      [47,] -71.53450 41.70089      [48,] -71.53282 41.70103      [49,] -71.53254 41.70104      [50,] -71.53379 41.70845      [51,] -71.53292 41.70870      [52,] -71.53302 41.71069      [53,] -71.53403 41.71306      [54,] -71.53391 41.71592      [55,] -71.53257 41.71872      [56,] -71.53371 41.71888      [57,] -71.53384 41.72104      [58,] -71.53265 41.72037      [59,] -71.53230 41.72094      [60,] -71.53395 41.72201      [61,] -71.53420 41.72399      [62,] -71.53371 41.72489      [63,] -71.53299 41.72506      [64,] -71.53340 41.72669      [65,] -71.53340 41.72692      [66,] -71.53418 41.72709      [67,] -71.53365 41.72785      [68,] -71.53420 41.72894      [69,] -71.53442 41.72861      [70,] -71.53443 41.73016      [71,] -71.54507 41.73000      [72,] -71.54674 41.73019      [73,] -71.54692 41.73026      [74,] -71.54747 41.73033      [75,] -71.54880 41.72993      [76,] -71.54987 41.72994      [77,] -71.56294 41.72983      [78,] -71.56864 41.72962      [79,] -71.57651 41.72934      [80,] -71.58756 41.71670      [81,] -71.59122 41.71699      [82,] -71.59468 41.72292      [83,] -71.60350 41.72682      [84,] -71.60391 41.72656      [85,] -71.61265 41.72805      [86,] -71.61665 41.72743      [87,] -71.61792 41.72816      [88,] -71.64093 41.72759      [89,] -71.67188 41.72730      [90,] -71.70062 41.72677      [91,] -71.72039 41.72619      [92,] -71.75311 41.72583      [93,] -71.75460 41.72573      [94,] -71.76842 41.72530      [95,] -71.78968 41.72473      [96,] -71.79125 41.77450      [97,] -71.79208 41.79509      [98,] -71.79274 41.80401      [99,] -71.79435 41.84168     [100,] -71.79448 41.84917     [101,] -71.79448 41.84930     [102,] -71.79493 41.85936     [103,] -71.79529 41.86575     [104,] -71.79652 41.89525     [105,] -71.79683 41.90445     [106,] -71.79773 41.92745     [107,] -71.79778 41.92855     [108,] -71.79830 41.95875     [109,] -71.79876 41.98509     [110,] -71.79924 42.00807     [111,] -71.75028 42.00995     [112,] -71.71357 42.01080     [113,] -71.71001 42.01095     [114,] -71.69858 42.01111     [115,] -71.68098 42.01147     [116,] -71.67803 42.01150     [117,] -71.60625 42.01336     [118,] -71.59078 42.01373     [119,] -71.57686 42.01396     [120,] -71.57342 42.01401     [121,] -71.56548 42.01418     [122,] -71.56365 42.01422     [123,] -71.53801 42.01480     [124,] -71.53751 42.01485     [125,] -71.52909 42.01500     [126,] -71.52871 42.01500     [127,] -71.52861 42.01500     [128,] -71.52774 42.01515     [129,] -71.52063 42.01540     [130,] -71.51777 42.01555     [131,] -71.51171 42.01610     [132,] -71.50618 42.01669     [133,] -71.50207 42.01699     [134,] -71.49876 42.01737     [135,] -71.49820 42.01727     [136,] -71.49100 42.01730     [137,] -71.48386 42.01750     [138,] -71.45816 42.01781     [139,] -71.41567 42.01849     [140,] -71.38140 42.01880     [141,] -71.38143 41.99999     [142,] -71.38130 41.98500     [143,] -71.38150 41.96670     [144,] -71.38150 41.95316     [145,] -71.38150 41.95200     [146,] -71.38173 41.92955     [147,] -71.38160 41.91630     [148,] -71.38150 41.91518     [149,] -71.38153 41.91046     [150,] -71.38170 41.89320     [151,] -71.37380 41.89440     [152,] -71.37122 41.89450     [153,] -71.36543 41.89545     [154,] -71.36166 41.89614     [155,] -71.35965 41.89645     [156,] -71.35690 41.89693     [157,] -71.35494 41.89724     [158,] -71.35451 41.89650     [159,] -71.34901 41.89716     [160,] -71.34674 41.89748     [161,] -71.33870 41.89840     [162,] -71.33907 41.89437     [163,] -71.33948 41.89268     [164,] -71.33978 41.89007     [165,] -71.34073 41.88242     [166,] -71.33412 41.86226     [167,] -71.33440 41.86130     [168,] -71.33490 41.86060     [169,] -71.34210 41.84520     [170,] -71.34193 41.84369     [171,] -71.34110 41.84290     [172,] -71.34080 41.84240     [173,] -71.34050 41.84210     [174,] -71.33522 41.83568     [175,] -71.33930 41.83230     [176,] -71.33960 41.83200     [177,] -71.34490 41.82800     [178,] -71.34715 41.82321     [179,] -71.33920 41.80900     [180,] -71.33924 41.80675     [181,] -71.33930 41.80440     [182,] -71.34098 41.79834     [183,] -71.33930 41.79630     [184,] -71.33395 41.79442     [185,] -71.33239 41.79232     [186,] -71.33274 41.79154     [187,] -71.33181 41.78972     [188,] -71.32939 41.78829     [189,] -71.32958 41.78264     [190,] -71.32782 41.78054     [191,] -71.32257 41.77823     [192,] -71.32301 41.77775     [193,] -71.32338 41.77794     [194,] -71.32412 41.77803     [195,] -71.32282 41.77701     [196,] -71.32225 41.77557     [197,] -71.32046 41.77384     [198,] -71.32251 41.77238     [199,] -71.32282 41.77379     [200,] -71.32263 41.77238     [201,] -71.32833 41.76831     [202,] -71.33153 41.76598     [203,] -71.33166 41.76595     [204,] -71.33233 41.76548     [205,] -71.33295 41.76597     [206,] -71.33706 41.76310     [207,] -71.33677 41.76284     [208,] -71.33621 41.76261     [209,] -71.33333 41.76235     [210,] -71.33280 41.76244     [211,] -71.33318 41.76175     [212,] -71.33771 41.76224     [213,] -71.34218 41.75909     [214,] -71.34223 41.75896     [215,] -71.34237 41.75897     [216,] -71.34385 41.75821     [217,] -71.34514 41.75671     [218,] -71.34727 41.75556     [219,] -71.34643 41.75547     [220,] -71.34683 41.75380     [221,] -71.34773 41.75395     [222,] -71.34768 41.75404     [223,] -71.34901 41.75451     [224,] -71.35177 41.75274     [225,] -71.35289 41.75403     [226,] -71.35071 41.75532     [227,] -71.35120 41.75546     [228,] -71.35249 41.75511     [229,] -71.35381 41.75628     [230,] -71.35257 41.75714     [231,] -71.35265 41.75784     [232,] -71.35377 41.75860     [233,] -71.35562 41.75771     [234,] -71.35509 41.75670     [235,] -71.35503 41.75385     [236,] -71.35618 41.75350     [237,] -71.35686 41.75331     [238,] -71.35782 41.75410     [239,] -71.35820 41.75474     [240,] -71.35860 41.75502     [241,] -71.35926 41.75284     [242,] -71.35560 41.75011     [243,] -71.35666 41.74854     [244,] -71.35576 41.74644     [245,] -71.35569 41.74633     [246,] -71.35542 41.74628     [247,] -71.35555 41.74625     [248,] -71.35557 41.74616     [249,] -71.35486 41.74605     [250,] -71.35433 41.74582     [251,] -71.35420 41.74530     [252,] -71.35431 41.74512     [253,] -71.35453 41.74518     [254,] -71.35541 41.74529     [255,] -71.35616 41.74520     [256,] -71.35855 41.74467     [257,] -71.35987 41.74518     [258,] -71.35964 41.74708     [259,] -71.36036 41.74972     [260,] -71.36233 41.75056     [261,] -71.36184 41.75296     [262,] -71.36168 41.75569     [263,] -71.36155 41.75972     [264,] -71.36249 41.76205     [265,] -71.36497 41.76461     [266,] -71.36779 41.76523     [267,] -71.37029 41.76480     [268,] -71.36843 41.76685     [269,] -71.36770 41.76846     [270,] -71.36709 41.76976     [271,] -71.36635 41.77085     [272,] -71.36575 41.77298     [273,] -71.36629 41.77338     [274,] -71.36632 41.77487     [275,] -71.36545 41.77723     [276,] -71.36647 41.77940     [277,] -71.36868 41.78141     [278,] -71.37053 41.78076     [279,] -71.37211 41.78340     [280,] -71.37067 41.78526     [281,] -71.37060 41.78710     [282,] -71.37238 41.79000     [283,] -71.37279 41.78833     [284,] -71.37362 41.79015     [285,] -71.37347 41.79043     [286,] -71.37411 41.79181     [287,] -71.37450 41.79098     [288,] -71.37439 41.79190     [289,] -71.37505 41.79218     [290,] -71.37690 41.79465     [291,] -71.37978 41.79584     [292,] -71.37881 41.79705     [293,] -71.38173 41.80067     [294,] -71.38625 41.80501     [295,] -71.38801 41.80415     [296,] -71.38785 41.80375     [297,] -71.38806 41.80418     [298,] -71.38635 41.80521     [299,] -71.38793 41.80817     [300,] -71.38824 41.80811     [301,] -71.39032 41.80715     [302,] -71.39373 41.81056     [303,] -71.39199 41.81158     [304,] -71.39437 41.81115     [305,] -71.39249 41.81202     [306,] -71.39136 41.81304     [307,] -71.39161 41.81510     [308,] -71.39197 41.81564     [309,] -71.39193 41.81584     [310,] -71.39172 41.81592     [311,] -71.39128 41.81589     [312,] -71.39114 41.81578     [313,] -71.39108 41.81551     [314,] -71.39051 41.81553     [315,] -71.38938 41.81579     [316,] -71.38782 41.81594     [317,] -71.38686 41.81676     [318,] -71.38644 41.81776     [319,] -71.38578 41.81830     [320,] -71.38602 41.81888     [321,] -71.38596 41.81902     [322,] -71.38563 41.81920     [323,] -71.38438 41.82070     [324,] -71.38336 41.82079     [325,] -71.38281 41.82142     [326,] -71.38314 41.82295     [327,] -71.38131 41.82416     [328,] -71.37913 41.82690     [329,] -71.37620 41.82912     [330,] -71.37624 41.82919     [331,] -71.37597 41.82938     [332,] -71.37346 41.83149     [333,] -71.37397 41.83219     [334,] -71.36973 41.83545     [335,] -71.36927 41.83608     [336,] -71.36911 41.83633     [337,] -71.36894 41.83668     [338,] -71.36875 41.83862     [339,] -71.37272 41.83905     [340,] -71.36992 41.84590     [341,] -71.36879 41.84648     [342,] -71.36839 41.84780     [343,] -71.36917 41.84856     [344,] -71.36773 41.84962     [345,] -71.37012 41.84935     [346,] -71.37180 41.85070     [347,] -71.37124 41.85179     [348,] -71.37129 41.85262     [349,] -71.36857 41.85271     [350,] -71.37282 41.85616     [351,] -71.37303 41.85699     [352,] -71.37487 41.85728     [353,] -71.37501 41.85743     [354,] -71.37565 41.86048     [355,] -71.37682 41.86144     [356,] -71.37867 41.86639     [357,] -71.38082 41.87020     [358,] -71.38281 41.87131     [359,] -71.38459 41.87264     [360,] -71.38188 41.86939     [361,] -71.38089 41.86874     [362,] -71.37994 41.86464     [363,] -71.37913 41.86291     [364,] -71.37925 41.86285     [365,] -71.38084 41.86106     [366,] -71.38036 41.85871     [367,] -71.37900 41.85739     [368,] -71.37891 41.85002     [369,] -71.37675 41.84877     [370,] -71.37908 41.84599     [371,] -71.37844 41.84464     [372,] -71.37726 41.84380     [373,] -71.37832 41.84328     [374,] -71.37776 41.84138     [375,] -71.37891 41.84095     [376,] -71.37974 41.83888     [377,] -71.37997 41.83893     [378,] -71.37986 41.83855     [379,] -71.37958 41.83820     [380,] -71.37955 41.83797     [381,] -71.37960 41.83657     [382,] -71.37917 41.83606     [383,] -71.37762 41.83332     [384,] -71.37945 41.83015     [385,] -71.37959 41.82936     [386,] -71.37955 41.82915     [387,] -71.38044 41.82776     [388,] -71.38062 41.82702     [389,] -71.38167 41.82629     [390,] -71.38292 41.82606     [391,] -71.38303 41.82594     [392,] -71.38429 41.82515     [393,] -71.38547 41.82470     [394,] -71.38585 41.82466     [395,] -71.38618 41.82423     [396,] -71.38752 41.82361     [397,] -71.38883 41.82145     [398,] -71.38895 41.82125     [399,] -71.38874 41.82108     [400,] -71.38854 41.81953     [401,] -71.38875 41.81925     [402,] -71.38865 41.81892     [403,] -71.38860 41.81877     [404,] -71.38817 41.81847     [405,] -71.38816 41.81702     [406,] -71.38976 41.81656     [407,] -71.39041 41.81698     [408,] -71.39137 41.81734     [409,] -71.39183 41.81778     [410,] -71.39262 41.81768     [411,] -71.39549 41.81764     [412,] -71.39725 41.81676     [413,] -71.39943 41.81453     [414,] -71.40280 41.81837     [415,] -71.40416 41.81985     [416,] -71.40446 41.82011     [417,] -71.40503 41.82100     [418,] -71.40599 41.82190     [419,] -71.40768 41.82401     [420,] -71.40783 41.82396     [421,] -71.40801 41.82390     [422,] -71.40581 41.82012     [423,] -71.40638 41.81980     [424,] -71.40587 41.81955     [425,] -71.40610 41.81920     [426,] -71.40549 41.81918     [427,] -71.40506 41.81851     [428,] -71.40491 41.81841     [429,] -71.40345 41.81571     [430,] -71.40304 41.81527     [431,] -71.40132 41.81215     [432,] -71.40258 41.81166     [433,] -71.40296 41.81127     [434,] -71.40248 41.81120     [435,] -71.40273 41.81045     [436,] -71.40057 41.81128     [437,] -71.40089 41.81072     [438,] -71.40248 41.81022     [439,] -71.40188 41.80936     [440,] -71.40008 41.80984     [441,] -71.40241 41.80896     [442,] -71.40121 41.80843     [443,] -71.40114 41.80812     [444,] -71.39988 41.80826     [445,] -71.40103 41.80769     [446,] -71.40060 41.80756     [447,] -71.39865 41.80786     [448,] -71.39865 41.80765     [449,] -71.40063 41.80700     [450,] -71.39951 41.80499     [451,] -71.39909 41.80491     [452,] -71.39831 41.80248     [453,] -71.39742 41.80240     [454,] -71.39709 41.80092     [455,] -71.39480 41.80241     [456,] -71.39746 41.80039     [457,] -71.39615 41.79976     [458,] -71.39514 41.79989     [459,] -71.39548 41.79885     [460,] -71.39257 41.80038     [461,] -71.38282 41.79296     [462,] -71.37879 41.78708     [463,] -71.38006 41.78697     [464,] -71.37967 41.78611     [465,] -71.38141 41.78497     [466,] -71.38408 41.78555     [467,] -71.38812 41.78451     [468,] -71.38759 41.78350     [469,] -71.38824 41.78429     [470,] -71.39054 41.78357     [471,] -71.39055 41.78338     [472,] -71.39084 41.78307     [473,] -71.39162 41.78288     [474,] -71.39151 41.78204     [475,] -71.39007 41.78200     [476,] -71.38965 41.78218     [477,] -71.38954 41.78203     [478,] -71.39135 41.78186     [479,] -71.39186 41.78168     [480,] -71.39154 41.78037     [481,] -71.39109 41.77983     [482,] -71.39093 41.77885     [483,] -71.39101 41.77751     [484,] -71.39024 41.77643     [485,] -71.38984 41.77627     [486,] -71.38990 41.77563     [487,] -71.38902 41.77515     [488,] -71.38871 41.77416     [489,] -71.38834 41.77377     [490,] -71.38916 41.76963     [491,] -71.38946 41.76893     [492,] -71.38885 41.76839     [493,] -71.38789 41.76788     [494,] -71.38724 41.76789     [495,] -71.38696 41.76758     [496,] -71.38621 41.76765     [497,] -71.38587 41.76819     [498,] -71.38610 41.76905     [499,] -71.38571 41.76832     [500,] -71.38577 41.76747     [501,] -71.38587 41.76445     [502,] -71.38518 41.76327     [503,] -71.38565 41.76288     [504,] -71.38585 41.76042     [505,] -71.38642 41.76025     [506,] -71.38685 41.76200     [507,] -71.38750 41.76332     [508,] -71.38801 41.76567     [509,] -71.38891 41.76570     [510,] -71.38915 41.76483     [511,] -71.39026 41.76465     [512,] -71.39053 41.76447     [513,] -71.38908 41.76368     [514,] -71.38907 41.76306     [515,] -71.38850 41.76260     [516,] -71.38793 41.76035     [517,] -71.38750 41.75915     [518,] -71.38708 41.75888     [519,] -71.38689 41.75771     [520,] -71.38721 41.75748     [521,] -71.38670 41.75592     [522,] -71.38622 41.75591     [523,] -71.38590 41.75524     [524,] -71.38392 41.75369     [525,] -71.38308 41.75374     [526,] -71.38108 41.75422     [527,] -71.38195 41.75390     [528,] -71.38259 41.75118     [529,] -71.38159 41.74968     [530,] -71.38356 41.74770     [531,] -71.38501 41.74788     [532,] -71.38481 41.74744     [533,] -71.38589 41.74675     [534,] -71.38789 41.74525     [535,] -71.39059 41.74565     [536,] -71.38988 41.74458     [537,] -71.38641 41.74456     [538,] -71.38558 41.74516     [539,] -71.38365 41.74623     [540,] -71.38317 41.74685     [541,] -71.38164 41.74585     [542,] -71.37792 41.74498     [543,] -71.37597 41.74409     [544,] -71.37709 41.74388     [545,] -71.38077 41.73938     [546,] -71.38144 41.73772     [547,] -71.38145 41.73684     [548,] -71.38074 41.73656     [549,] -71.38093 41.73612     [550,] -71.38093 41.73553     [551,] -71.38122 41.73530     [552,] -71.38195 41.73556     [553,] -71.38328 41.73636     [554,] -71.38285 41.73769     [555,] -71.38329 41.73911     [556,] -71.38455 41.73969     [557,] -71.38422 41.73872     [558,] -71.38831 41.73691     [559,] -71.38865 41.73745     [560,] -71.38996 41.73737     [561,] -71.38911 41.73379     [562,] -71.39294 41.73369     [563,] -71.39362 41.73320     [564,] -71.39014 41.73235     [565,] -71.38900 41.73179     [566,] -71.38640 41.73356     [567,] -71.38550 41.73269     [568,] -71.38548 41.73155     [569,] -71.38579 41.73078     [570,] -71.38292 41.72877     [571,] -71.38266 41.72853     [572,] -71.38255 41.72839     [573,] -71.38046 41.72669     [574,] -71.37939 41.72437     [575,] -71.37694 41.72224     [576,] -71.37064 41.71869     [577,] -71.36481 41.71744     [578,] -71.35555 41.71801     [579,] -71.35452 41.71765     [580,] -71.35688 41.71773     [581,] -71.35864 41.71611     [582,] -71.36326 41.71371     [583,] -71.36357 41.71342     [584,] -71.36449 41.71379     [585,] -71.36401 41.71438     [586,] -71.36459 41.71508     [587,] -71.36459 41.71546     [588,] -71.36474 41.71562     [589,] -71.36487 41.71559     [590,] -71.36471 41.71589     [591,] -71.36486 41.71602     [592,] -71.36526 41.71602     [593,] -71.36523 41.71632     [594,] -71.36559 41.71627     [595,] -71.36551 41.71614     [596,] -71.36480 41.71576     [597,] -71.36486 41.71571     [598,] -71.36500 41.71561     [599,] -71.36473 41.71557     [600,] -71.36470 41.71491     [601,] -71.36418 41.71459     [602,] -71.36453 41.71379     [603,] -71.36591 41.71282     [604,] -71.36644 41.71349     [605,] -71.36731 41.71365     [606,] -71.36856 41.71475     [607,] -71.36878 41.71476     [608,] -71.37066 41.71560     [609,] -71.37209 41.71534     [610,] -71.37044 41.71515     [611,] -71.37050 41.71399     [612,] -71.36974 41.71413     [613,] -71.36997 41.71354     [614,] -71.37128 41.71290     [615,] -71.37410 41.71267     [616,] -71.37490 41.71225     [617,] -71.37510 41.71218     [618,] -71.37254 41.71245     [619,] -71.37085 41.71212     [620,] -71.36919 41.71366     [621,] -71.36672 41.71159     [622,] -71.36559 41.70783     [623,] -71.36564 41.70781     [624,] -71.36551 41.70671     [625,] -71.36527 41.70647     [626,] -71.36531 41.70615     [627,] -71.36458 41.70351     [628,] -71.36509 41.70128     [629,] -71.36526 41.70081     [630,] -71.36509 41.70056     [631,] -71.36596 41.69876     [632,] -71.36598 41.69847     [633,] -71.36605 41.69777     [634,] -71.36552 41.69313     [635,] -71.36301 41.69221     [636,] -71.36492 41.68876     [637,] -71.37123 41.68595     [638,] -71.37379 41.67847     [639,] -71.37296 41.67245     [640,] -71.37819 41.66647     [641,] -71.38443 41.66870     [642,] -71.38618 41.67093     [643,] -71.39039 41.67217     [644,] -71.39054 41.68272     [645,] -71.38891 41.68354     [646,] -71.38965 41.68533     [647,] -71.38941 41.68609     [648,] -71.39173 41.68668     [649,] -71.38912 41.68776     [650,] -71.39043 41.68830     [651,] -71.39154 41.68862     [652,] -71.39105 41.69023     [653,] -71.38940 41.69110     [654,] -71.38851 41.69451     [655,] -71.38911 41.69538     [656,] -71.38687 41.69575     [657,] -71.38612 41.69452     [658,] -71.38539 41.69550     [659,] -71.38596 41.69664     [660,] -71.38266 41.69788     [661,] -71.38141 41.69780     [662,] -71.38165 41.69912     [663,] -71.38353 41.69960     [664,] -71.38497 41.69916     [665,] -71.38530 41.69926     [666,] -71.38604 41.69818     [667,] -71.38723 41.69804     [668,] -71.39116 41.69647     [669,] -71.39209 41.69718     [670,] -71.39237 41.69701     [671,] -71.39074 41.69488     [672,] -71.39161 41.69278     [673,] -71.39331 41.69297     [674,] -71.39428 41.69518     [675,] -71.39756 41.69658     [676,] -71.39487 41.69450     [677,] -71.39426 41.69414     [678,] -71.39463 41.69364     [679,] -71.39331 41.69169     [680,] -71.39446 41.69017     [681,] -71.39383 41.68740     [682,] -71.39294 41.68683     [683,] -71.39354 41.68624     [684,] -71.39291 41.68563     [685,] -71.39325 41.68505     [686,] -71.39283 41.68401     [687,] -71.39666 41.68354     [688,] -71.39984 41.68389     [689,] -71.39977 41.68433     [690,] -71.40140 41.68503     [691,] -71.40144 41.68537     [692,] -71.40248 41.68641     [693,] -71.40237 41.68769     [694,] -71.40160 41.68740     [695,] -71.40090 41.68753     [696,] -71.40092 41.68915     [697,] -71.40120 41.68798     [698,] -71.40225 41.68782     [699,] -71.40289 41.68921     [700,] -71.40389 41.68996     [701,] -71.40342 41.69138     [702,] -71.40438 41.69195     [703,] -71.40420 41.69267     [704,] -71.40568 41.69359     [705,] -71.40652 41.69550     [706,] -71.40661 41.69595     [707,] -71.40806 41.69560     [708,] -71.41006 41.69623     [709,] -71.41139 41.69723     [710,] -71.41155 41.69722     [711,] -71.41205 41.69741     [712,] -71.41371 41.69775     [713,] -71.41632 41.69937     [714,] -71.41634 41.69942     [715,] -71.41770 41.70242     [716,] -71.41680 41.69903     [717,] -71.42089 41.69856     [718,] -71.42106 41.69844     [719,] -71.42041 41.69825     [720,] -71.41675 41.69815     [721,] -71.41556 41.69825     [722,] -71.41528 41.69749     [723,] -71.41421 41.69696     [724,] -71.41334 41.69686     [725,] -71.41203 41.69633     [726,] -71.41036 41.69506     [727,] -71.40848 41.69443     [728,] -71.40456 41.68915     [729,] -71.40814 41.68709     [730,] -71.41038 41.68861     [731,] -71.41230 41.68846     [732,] -71.41562 41.68966     [733,] -71.41872 41.68999     [734,] -71.41698 41.68963     [735,] -71.41481 41.68808     [736,] -71.40748 41.68544     [737,] -71.40509 41.68567     [738,] -71.40306 41.68559     [739,] -71.40580 41.68449     [740,] -71.41477 41.68412     [741,] -71.42322 41.68574     [742,] -71.42863 41.68605     [743,] -71.43070 41.68573     [744,] -71.43404 41.68681     [745,] -71.43573 41.68576     [746,] -71.43688 41.68698     [747,] -71.44195 41.68695     [748,] -71.44426 41.68872     [749,] -71.44338 41.68910     [750,] -71.44274 41.68898     [751,] -71.44305 41.69261     [752,] -71.44517 41.69300     [753,] -71.44536 41.69318     [754,] -71.44518 41.69284     [755,] -71.44639 41.69485     [756,] -71.44653 41.69535     [757,] -71.44720 41.69587     [758,] -71.44952 41.69582     [759,] -71.45125 41.69620     [760,] -71.45157 41.69733     [761,] -71.45329 41.69655     [762,] -71.45207 41.69504     [763,] -71.45183 41.69353     [764,] -71.45057 41.69236     [765,] -71.44731 41.69200     [766,] -71.44759 41.69124     [767,] -71.44639 41.69036     [768,] -71.44928 41.68818     [769,] -71.45003 41.68692     [770,] -71.44944 41.68653     [771,] -71.44924 41.68503     [772,] -71.45091 41.68471     [773,] -71.44938 41.68180     [774,] -71.44775 41.67693     [775,] -71.44855 41.67696     [776,] -71.44920 41.67844     [777,] -71.45041 41.67795     [778,] -71.45211 41.67831     [779,] -71.45270 41.67765     [780,] -71.45866 41.67751     [781,] -71.45771 41.67379     [782,] -71.45784 41.67243     [783,] -71.45784 41.67210     [784,] -71.45804 41.67023     [785,] -71.46689 41.66987     [786,] -71.46871 41.67024     [787,] -71.47862 41.67031     [788,] -71.47995 41.67029     [789,] -71.47718 41.67402     [790,] -71.47734 41.68254     [791,] -71.50031 41.68146                    Slot "plotOrder":     [1] 1          Slot "labpt":     [1] -71.57027  41.85902          Slot "ID":     [1] "TRUE"          Slot "area":     [1] 0.1317919                    Slot "plotOrder":     [1] 1          Slot "bbox":             min       max     x -71.79924 -71.32046     y  41.66166  42.01880          Slot "proj4string":     CRS arguments: NA


How can I determine when a Bing map waypoint is going to overlap an existing one, and give the original waypoint a “badge-like” update?

I want to be able to determine when a Bing map waypoint that I want to add is going to overlap an existing one, so that, rather than bunching them all up in one spot like a litter of kittens on a catnip plant, I can instead alter the text/"badge" of the existing waypoint marker so that it displays an indication that it has been "reached" N times. For example, in showing where a person was, and marking the waypoint with a sequential number, I may want to update the text on the waypoint marker from "1" to "1,2" and then "1-3" or possibly just "1+" as there won't be enough room to display something like "1,7,42" etc.