Allowing cross-domain access of third-party WMS server from OpenLayers

I have a web site with an OpenLayers client. I also have a GeoServer running on the same domain. My OpenLayers uses services like WMS, GML, etc. from this GeoServer installation and everything runs fine. Now, I want to add another WMS layer to my existing set of layers. This layer is provided by a third-party WMS (not owned by me). How can I add this layer to my OpenLayers by bypassing the cross-domain issue. PS: I am aware that GeoServer allows adding third-party WMS as a new data store and lets me publish the available layers, but I am facing some projection issues there which I will be asking in a different question.


Using Tilestache and Leaflet To Display Rasters In Native Projection

I have a collection of old Idaho state highway maps in GeoTiff format. The projection that these maps have been georeferenced in is a transverse mercator projection, called IDTM83 (Idaho Transverse Mercator - NAD 83). The projection is defined as: > +proj=tmerc +lat_0=42 +lon_0=-114 +k=0.9996 +x_0=2500000 +y_0=1200000 > +ellps=GRS80 +units=m +no_defs The projection is basically a modified UTM projection that compensates for that fact that the state of Idaho is split by UTM zones 11 and 12. My goal is to serve the maps in the native IDTM83 projection using Tilestache and Leaflet. I am able to serve these same rasters that have already been projected into world mercator, but the the projection distorts much of the line work from the original map. There have been some helpful resources on the internet on how to serve projected data using Tilestache and Leaflet, here:   1. https://gist.github.com/skogsmaskin/6452343   2. http://blog.thematicmapping.org/2012/07/using-custom-projections-with-tilecache.html.  The mapnik configuration for one of my raster maps looks like:               <Map background-color="steelblue"      srs="+proj=tmerc +lat_0=42 +lon_0=-114 +k=0.9996 +x_0=2500000 +y_0=1200000 +ellps=GRS80 +units=m +no_defs"     maximum-extent="2205517,1134032,2834671,2040323">     <Style name="raster-style">       <Rule>           <RasterSymbolizer/>       </Rule>     </Style>          <Layer name="idaho-1947" status="on" srs="+proj=tmerc +lat_0=42 +lon_0=-114 +k=0.9996 +x_0=2500000 +y_0=1200000 +ellps=GRS80 +units=m +no_defs">     <StyleName>raster-style</StyleName>     <Datasource>         <Parameter name="type">gdal</Parameter>         <Parameter name="file">GIS/IdahoMaps/1947_geo_idtm.tif</Parameter>         <Parameter name="format">tiff</Parameter>     </Datasource>     </Layer>     </Map>      Using Python, I am able to generate maps using bounding boxes from the IDTM coordinate system with the Mapnik API. So I know that Mapnik is not the problem. My Tilestache config is pretty basic:                   {       "cache": {"name": "Test","verbose":"True"},       "layers": {         "idaho-1947": {             "provider": {"name": "mapnik", "mapfile": "./mapnik.xml"},             "projection": "TileStache.Goodies.Proj4Projection:Proj4Projection('+proj=tmerc +lat_0=42 +lon_0=-114 +k=0.9996 +x_0=2500000 +y_0=1200000 +ellps=GRS80 +units=m +no_defs',[5545984, 2772992, 1386496, 693248, 346624, 173312, 86656, 43328, 21664, 10832, 5416, 2708, 1354, 677, 338.5, 169.25, 84.625, 42.3125, 21.15625, 10.578125, 5.2890625, 1])"         }}      My leaflet viewer is very basic, set up like on Link #2 above:                 var crs = L.CRS.proj4js('SR-ORG:6677', '+proj=tmerc +lat_0=42 +lon_0=-114 +k=0.9996 +x_0=2500000 +y_0=1200000 +ellps=GRS80 +units=m +no_defs', new L.Transformation(1, -432000, -1, 2040323));       crs.scale = function(zoom) {           return 1 / (234.375 / Math.pow(2, zoom));       };            var map = new L.Map('map', {         crs: crs,         layers: [           new L.TileLayer('http://localhost:8080/idaho-1947/{z}/{x}/{y}.png', {             minZoom: 0,             maxZoom: 4,             continuousWorld: true           })         ],         center: new L.LatLng(45, -116),         zoom: 3,         continuousWorld: true       });      I intend to refine this once I can get the maps to display, but the example is using UTM data and the scales should be pretty similar. I am stuck here. Using the JavaScript examples provided in both the links above, I cannot display any of the maps - only the lightblue background. I know what I am trying to do is possible in ArcGIS Server, but I would really like to leverage some of the OpenSource tools.


How to create custom markers for QGIS?

Is it possible to define a custom marker as symbolizer for a point in QGIS? I'm looking into using a variation on the "Vector Field Marker" with a couple of arrows pointing at dynamic angles (the angles being read as point attributes). I imagine it's doable in C++ (within the QGIS code) but I'm looking for a scripting solution (python). I'm adding a representation of what I'm trying to achieve - for two points A and B with attributes x,y,z. ![custom markers](http://i.stack.imgur.com/Af3Al.png)


geodatabase and relationship

I have two tables in a personal geodatabase and one of them has 13 records and 6 fields, the other has 118 records and 125 fields. There is a relationship class between the two tables where the first table is the origin and the other (bigger) table is the destination. What I need is to take the records satisfying a condition from the second table and sum those records and put into a new field in the origin table


what happens to attribute table while performing union, intersection, difference of gml and shp file in QGIS

I am trying to do an analysis of the data conversion from shapefile to GML file using QGIS. For the analysis I have taken a point, line and polygon file. I first took a shapefile and converted it into a GML file; after converting it I compared the feature count of both files, which were same. I then performed geoprocessing analyses (difference, intersection, union, etc) on both of these files. I am doing a geoprocessing analysis because if the shapefile has been successfully converted into GML without errors, then a difference of both these files will give me a null set. Similarly for others: http://infogeoblog.wordpress.com/2013/01/08/geo- processing-in-qgis/#comments I found that for union, it is not showing the number of tuples in the resultant attribute table as sum of tuples in both GML file attribute table and shapefile attribute table. Similarly, for intersection and difference it is not showing the results as expected (by expected I mean, as we do in set theory). For some files it also shows null values. I want to know what happens to tuples in the attribute table while performing union (intersection, difference) of GML and shapefiles. Do they get added up or is some join performed?


PostGIS | Get degrees clockwise from north of a linestring's beginning and end

I have a table full of linestrings in PostGIS. I would like to know the degrees clockwise from north of the starting point of the linestring and also of the end point. I was thinking that maybe creating a circle with both points could help to get this two values (example: X=32º, Y=209º) but i have no clue how to achieve this.... ![enter image description here](http://i.stack.imgur.com/kzFrn.png) Regards


How do you assign classes to percentage data in a data layer?

As a learning project I need to coneptualize a GIS output for assigning classes to data in a raster layer. In this case i have the unemployment rates for relevant magisterial districts in my database. So I need to create a layer where I can assign classes to the unemployment rate (percentage of people unemployed) as low (0-10%), medium (10-30%), or high (40&+) within each district. I've thought that a raster layer containing the unemployment percentages as cell values could work, but how would I go about outputting this through a GIS? I'm sorry if this question is a bit non-specific, but any tips will be greatly appreciated! Thank you


How to check if REST plugin is installed for GeoServer

Hi I am having trouble uploading vector files onto GeoServer 2.4 using 'Import Vector' in the QGIS Processing toolbox. After reading around about the 'REST' plus in I went to: http:// Also when I run the tool I find an error code 404 in the log. GeoServer is running on a company server along witch apache 6.0 Could someone please tell me how to check the REST plug in, and if it is not installed how to get the REST plug in for GeoServer 2.4.0. Thank you, Sean0


How to Change the Default Displaying Order of Esri Basemaps in basemapGallery?

I just started using ArcGIS JS APIs. In my app, I want to combine a few ESRI basemaps with my own basemap service. I am using basemapGallery dijit to change basemaps. By default, basemapGallery shows the imagery map. I want to show OpenStreetMap as the default. How can I make this happen? The code in my app is simple and directly borrowed from ArcGIS samples. So I did not put it here. Thanks!


Don't move Layout

I'm doing a vb.net button to run in ARCGIS. the button is working, the problem is as I'm working in the layot when I press the button, change the layout coordinates and the map coordinates. Is a way that the layout view don't move and keep the Coordinates (0,0)               I leave the code:          MxDoc = My.ArcMap.Document          ActiveView = MxDoc.ActiveView          Dim x As Double = txtEast.Text          Dim y As Double = txtNorth.Text          Map.MapScale = 4000 'zoom to this scale          Envelope = ActiveView.Extent          ,Create and Add the point coordinates     Point = New ESRI.ArcGIS.Geometry.Point     Point.PutCoords(x, y)          'Envelope Position go to the middle of the point     Envelope.CenterAt(Point)               'Refresh View      ActiveView.Extent = Envelope      MxDoc.ActiveView.Refresh()Map = MxDoc.FocusMap      Thank You!


Openlayers map save feature SRID error into PostGis

OpenLayers.ProxyHost = "proxy.ashx?url=";                  //set up the modification tools             var DeleteFeature = OpenLayers.Class(OpenLayers.Control, {                 initialize: function(layer, options) {                 OpenLayers.Control.prototype.initialize.apply(this, [options]);                 this.layer = layer;                 this.handler = new OpenLayers.Handler.Feature(                     this, layer, {click: this.clickFeature}                 );                 },                 clickFeature: function(feature) {                 // if feature doesn't have a fid, destroy it                 if(feature.fid == undefined) {                     this.layer.destroyFeatures([feature]);                 } else {                     feature.state = OpenLayers.State.DELETE;                     this.layer.events.triggerEvent("afterfeaturemodified",                                                 {feature: feature});                     feature.renderIntent = "select";                     this.layer.drawFeature(feature);                 }                 },                 setMap: function(map) {                 this.handler.setMap(map);                 OpenLayers.Control.prototype.setMap.apply(this, arguments);                 },                 CLASS_NAME: "OpenLayers.Control.DeleteFeature"             });                  function showMsg(szMessage) {                 document.getElementById("nodelist").innerHTML = szMessage;                 setTimeout(                 "document.getElementById('nodelist').innerHTML = ''",2000);             }                  function showSuccessMsg(){             alert("Succes");                 showMsg("Transaction successfully completed");             };                  function showFailureMsg(){             alert("NoneSucces");                 showMsg("An error occured while operating the transaction");             };                  $(document).on("pageshow", "#mappage", function( event ) {                          haritacontrol=1;                      //set up a save strategy                 var saveStrategy = new OpenLayers.Strategy.Save();                     saveStrategy.events.register("success", '', showSuccessMsg);                     saveStrategy.events.register("failure", '', showFailureMsg);                        var bounds = new OpenLayers.Bounds(417138.6205791628, 4243634.852681775,486211.1950145211, 4331542.555877601                     );                     var options = {                         controls:[                         //allows the user pan ability                         new OpenLayers.Control.Navigation(),                         //displays the pan/zoom tools                                            new OpenLayers.Control.PanZoom(),                         //displays a layer switcher                         new OpenLayers.Control.LayerSwitcher(),                          //displays the mouse position's coordinates in a                          //<div> html element with id="coordinates"                         new OpenLayers.Control.MousePosition(                         {div:document.getElementById("coordinates")})                        ],                         maxExtent: bounds,                         maxResolution: 343.3894656086959,                         projection: "EPSG:2322",                         units: 'm'                     };                     map = new OpenLayers.Map('map', options);                          // setup tiled layer                     tiled = new OpenLayers.Layer.WMS(                         "Geoserver layers - Tiled", "http://192.168.1.236:8080/geoserver/Kaski/wms",                         {                             LAYERS: 'KaskiGrup',                                       tiled: true,                             tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                         },                         {                             buffer: 0,                             displayOutsideMaxExtent: true,                             isBaseLayer: true,                             yx : {'EPSG:2322' : true}                         }                      );                      //wfs-t editable overlay                 var wfs_layer = new OpenLayers.Layer.Vector("Editable Features", {                     strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],                          protocol: new OpenLayers.Protocol.WFS({                         version: "1.0.0",                         // loading data through localhost url path                         url: "http://192.168.1.236:8080/geoserver/wfs",                         featureNS :  "http://192.168.1.236:8080/KASKIDB",                         maxExtent: bounds,                         // layer name                         featureType: "DBadres_mahalle",                         // geometry column name                              geometryName: "geom",                         schema: "http://192.168.1.236:8080/geoserver/wfs/DescribeFeatureType?version=1.0.0&typename=KASKIDB:DBadres_mahalle"                     })                     });                      map.addLayers([tiled, wfs_layer]);                 map.zoomToExtent(bounds);                      // add the custom editing toolbar                 var panel = new OpenLayers.Control.Panel(                     {'displayClass': 'customEditingToolbar'}                     );                      var navigate = new OpenLayers.Control.Navigation({                     title: "Pan Map"                     });                      var draw = new OpenLayers.Control.DrawFeature(                     wfs_layer, OpenLayers.Handler.Polygon,                     {                         title: "Draw Feature",                         displayClass: "olControlDrawFeaturePolygon",                         multi: true                     }                     );                                var edit = new OpenLayers.Control.ModifyFeature(wfs_layer, {                     title: "Modify Feature",                     displayClass: "olControlModifyFeature"                     });                      var del = new DeleteFeature(wfs_layer, {title: "Delete Feature"});                      var save = new OpenLayers.Control.Button({                     title: "Save Changes",                     trigger: function() {                         if(edit.feature) {                         edit.selectControl.unselectAll();                         }                         saveStrategy.save();                              // alert('saved');                     },                     displayClass: "olControlSaveFeatures"                     });                      panel.addControls([navigate, save, del, edit, draw]);                 panel.defaultControl = navigate;                 map.addControl(panel);                               });      My code is like that, but when I save feature, it is added but there is no information on table. For example fid columns nothing. Also when I click saved feature, i got this error :               Internal error occurred java.io.IOException null ERROR: Operation on mixed SRID geometries      How can I give columns parameters. Also I found another example about insert but it gives success WFS transaction but nothing change. That code :                   var map;         OpenLayers.ProxyHost = "proxy.ashx?url=";                    $(document).on("pageshow", "#mappage", function( event ) {           var bounds = new OpenLayers.Bounds(417138.6205791628, 4243634.852681775,486211.1950145211, 4331542.555877601             );             var options = {                 controls: [ new OpenLayers.Control.Zoom(),                     new OpenLayers.Control.Navigation()],                 maxExtent: bounds,                 maxResolution: 343.3894656086959,                 projection: "EPSG:2322",                 units: 'm'             };             map = new OpenLayers.Map('map', options);                  // setup tiled layer             tiled = new OpenLayers.Layer.WMS(                 "Geoserver layers - Tiled", "http://192.168.1.236:8080/geoserver/Kaski/wms",                 {                     LAYERS: 'KaskiGrup',                               tiled: true,                     tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                 },                 {                     buffer: 0,                     displayOutsideMaxExtent: true,                     isBaseLayer: true,                     yx : {'EPSG:2322' : true}                 }              );         var roads = new OpenLayers.Layer.WFS(             "Roads",              "http://192.168.1.236:8080/geoserver/wfs",             {typename: 'KASKIDB:DBadres_mahalle'},             {                 typename: 'DBadres_mahalle',                  featureNS: 'http://192.168.1.236:8080/KASKIDB',                  geometryName: "geom",                 extractAttributes: false             }         );         roads.style = OpenLayers.Util.applyDefaults({strokeColor: "#ff0000"},                  OpenLayers.Feature.Vector.style["default"]);                        map.addLayers([tiled,roads]);              var panel = new OpenLayers.Control.Panel(             {displayClass: 'olControlEditingToolbar'}         );              var drawPol = new OpenLayers.Control.DrawFeature(             roads, OpenLayers.Handler.Polygon,             {displayClass: 'olControlDrawFeaturePath'}         );         drawPol.featureAdded = function(feature) {             feature.layer.eraseFeatures([feature]);             // cast to multilinestring             feature.geometry = new OpenLayers.Geometry.MultiPolygon(                 feature.geometry             );             feature.style.strokeColor = "#ff0000";             feature.state = OpenLayers.State.INSERT;             feature.layer.drawFeature(feature);         }                   var drawPoint = new OpenLayers.Control.DrawFeature(             roads, OpenLayers.Handler.Point,             {displayClass: 'olControlDrawFeaturePoint'}         );         drawPoint.featureAdded = function(feature) {             feature.layer.eraseFeatures([feature]);             // cast to multipoint             feature.geometry = new OpenLayers.Geometry.MultiPoint(                 feature.geometry             );             feature.style.strokeColor = "#0000ff";             feature.state = OpenLayers.State.INSERT;             feature.layer.drawFeature(feature);         }                   panel.addControls(             [new OpenLayers.Control.Navigation(), drawPol]         );              map.addControl(panel);         map.addControl(new OpenLayers.Control.LayerSwitcher());              map.zoomToExtent(bounds);                                                                                });


QGIS 2.0 WMTS Issue

When using a WMTS within QGIS 2.0 (using the "?SERVICE=WMTS&REQUEST=GetCapabilities" addition) only one tile per layer is loading. The WMTS works in other GIS software and this problem would appear to be unique to QGIS 2.0. Does anyone know why this would he happening and if there is a workaround to this issue?


(How) Can I define my own placement settings for labels?

Is it possible to define my own Placement settings for my labels? ![Placement properties](http://i.stack.imgur.com/kU8da.png) The predefined ones are not bad, but I would need a few more of my own design. How can I accomplish that?


Using ECQL filter on user form using GeoServer and OpenLayers

I know this is more of a JS question (skills of which I am inadequately inept) but it is using OpenLayers and GeoServer. I am trying to apply a filter with a user form on a Layer, but am not sure why it is not working. If I manually write the code into the CQL_FILTER, it works, but it won't work with the UpdateFilter button.               function init_map() {              var layerExtent = new OpenLayers.Bounds( -9816656, 5296800, -9772016, 5345759);              map = new OpenLayers.Map('map',                  {'restrictedExtent': layerExtent});         basemap = new OpenLayers.Layer.OSM("Simple OSM Map");         map.addLayer(basemap);         markers = new OpenLayers.Layer.Markers("Markers");         map.addLayer(markers);         map.setCenter(             new OpenLayers.LonLat(-71.147, 42.472).transform(                 new OpenLayers.Projection("EPSG:4326"),                 map.getProjectionObject()             ), 12         );              var tiled = new OpenLayers.Layer.WMS(                     "Verizon Cell Towers", "http://mpdflash:7070/geoserver/Milwaukee/wms",                     {                         LAYERS: 'Milwaukee:Verizon Cell Towers',                         tiled: true,                         transparent: true                     },                     {                         displayOutsideMaxExtent: true,                         isBaseLayer: false,                         yx : {'EPSG:6948' : false}                     }                  );                map.addLayer(tiled);              tiled.mergeNewParams({'CQL_FILTER': "filter"})     }          function updateFilter(){             var filter = document.getElementById('filter').value;         }          function resetFilter() {             document.getElementById('filter').value = "";             updateFilter();         }         and
                 <input type="text" size="80" id="filter"/>


Where to get Annual Lightning Data (Zipcode) to use with ArcGIS for Desktop and AnyLogic?

I am building a Dynamic Agent Based Model for School and was looking for Annual Lighting Strike Data for zip codes in Ohio. I can get the state count from NOAA but want to get a bit smaller. I was wondering if anyone knew of a source? (I am using the data to analysis the risk of lighting in and around the Cincinnati Area.) I am using ArcGIS for Desktop to do the maps and **AnyLogic** for my Model.


MapInfo 10.5.1 Update and MS Access Database

I am having problems with a database that I have developed. The data uses MapInfo for any of the spatial parts of it. I installed the update to MapInfo 10.5.1 and I now get an error when I am trying to access an ADODB.recordset. object which previously worked. I get the error: Method 'Connection' of object '_CurrentProject' failed when the database tries to run the code:               rs5.ActiveConnection = CurrentProject.Connection       (where rs5 is a ADODB.Recordset): Has anyone else seen anything like this? Is this related to the updated drivers for Office 2010 that were installed with the update? I am using Access 2007. The database still works on machines that are 10.5.


Configuring & debugging a new mapnik / tirex installation

I've a stack OpenLayers > Apache mod_tile > tirex working with the tirex test backend on a Ubuntu 10.04 virtual machine. However, I now have troubles with the mapnik-example backend included with tirex. OpenLayers works fine, mod_tile does not complain about the newly added "AddTileConfig /tiles/example/ example" to my Apache httpd.conf. Also, the tirex-master receives a request from mod_tile. I suspect this problem to lie with the configuration of the tirex Mapnik rendering backend. Using tirex-status command I can see tirex-master currently rendering the tiles for the example map. The rendering process times out after about 300 seconds and no (meta)tiles are being written to the tile directory. I suspect a wrong configuration or problems with user/write access. Does anyone know where I can find debug information from tirex or mapnik that can help me solve this problem? Thanks in advance. * * * EDIT: I am looking for clues about what is going on in some logfiles. I suspect the problem to lie with the tirex-backend-mapnik process, which seems to be managed by the tirex-backend-manager process. In the [OSM Tirex Logging wiki page][1] it says that the tirex-master and tirex-backend-manager will log to the syslog by default, which should be located under /var/log/messages. This logfile has the following entry at the end ` Dec 23 12:50:13 ubuntumapserver kernel: [10967.139797] mapnik[2079]: segfault at 0 ip 000000000041314b sp 00007fffacfa9530 error 4 in mapnik[400000+1f000]   Dec 23 12:50:14 ubuntumapserver kernel: [10968.391401] mapnik[2077]: segfault at 0 ip 000000000041314b sp 00007fffe23608a0 error 4 in mapnik[400000+1f000] ` Does anyone know have any clues about what can be causing this error? * * * EDIT 2: Hi Dane, Thanks for your reply. I installed Mapnik from source using the [Mapnik wiki tutorial][2] (Installation on a VM guest running Ubuntu 10.04 Server x64) ` python scons/scons.py configure DEBUG=True INPUT_PLUGINS=all OPTIMIZATION=3 SYSTEM_FONTS=/usr/share/fonts/ python scons/scons.py sudo python scons/scons.py install ` I have successfully used Mapnik standalone via the Python commmand line to render images. Furthermore I did an strace during a request from OpenLayers. (I saw you recommending this elsewhere). I can see mod_tile.so being called. However, there is no data about tirex or mapnik. [strace][3]   [config.log from Mapnik installation][4]   Therefore, I guess you might be correct about tirex **not working with a compiled version of Mapnik**. I will try to find out if this is indeed the problem. 1_http://wiki.openstreetmap.org/wiki/Tirex/Logging (About Tirex Logging)   2_http://trac.mapnik.org/wiki/UbuntuInstallationOld (Mapnik installation)   3_http://pastebin.com/L8VBzvgM (strace of a tile request)   4_http://pastebin.com/yvLQDr39 (Mapnik scons config.log)


How to create regular point grid inside a polygon in Postgis?

How to create, inside a polygon, a regular grid of point spaced x,y in postgis? Like the example: ![alt text](http://i.stack.imgur.com/ynKBu.png)


Transform wkb to gml geometry using FME

Maybe (and probably) this is stupid question, but I am totaly new to FME and geometry conversion. I need to transform WKB field from SQLServer to GML geometry using FME software. How can I do this?


Count points in a polygon that changes form ("on the fly calculation")

I have tried to find a plugin for this, with no luck. Perhaps someone knows of a plugin that can do the following: In this project I create a lot of polygons that contain points, and the size of the polygons is dependent of how many points there are inside the polygon. Therefore I need to manually change the size of the polygons to get a certain number of points inside it. It is tedious work to run a new "Points in polygon" command as soon as I have changed a polygon and want to know the number of points. Is there any way to get this calculations directly on screen some how (a sort of "on the fly" point calculation)?


Tool in ArcGIS for stream flow analysis?

I'm just wondering if there is any tool in ArcGIS that can be used to analyse stream flow measurements. For example, estimating runoff, baseflow and baseflow index based on stream flow measurements.


How to store esri geometry types in sql server

I have a txt file containing some data of polylines. Actually it is saved json objects. The features in it are esriPolyLine objects. Now i want to store it in my sql server db as geography data types through EF. Is there anywhere some great examples for doing this. How should the string value loop like to be able to enter the values using DbGeography.From(polylinestring, 4326) Or just a simple example of storing esriGeographyPolyline objects in sqlserver


Using Arcgis geometry in OpenLayers

I am sending request to Arcgis Rest service and getting result as JSON object. ![enter image description here](http://i.stack.imgur.com/rbomz.png) I want to zoom to this Ring array.               var extent = new OpenLayers.Bounds(extent.geometry.rings[0]);      I tried but did not work. I am new at OpenLayers.


Open source alternative to ArcGIS geoprocessing service concept?

I am searching for any examples of implementing a logic that Esri has built with geoprocessing services. You author a custom tool or it can be a system tool >> you publish it exposing this functionality as a web service >> you consume the web service with a client. A classical simple example: user can click on the map to create a point > the coordinates are sent to the server > buffering operation is being executed at the back-end > the buffer zone feature (or just an image) is sent back to client and shown on the map canvas. I've been developing Esri GP services for last 3 years and it has been a nice experience specifically since version 10.1. **Are there any examples of the open source system/solution that one could use to build a similar logic?** I've googled and found sextante gis but there was nothing specific on how this works and how much of logic that is available through Esri is available. 52North also seem to be working on that, but again I am missing the clean explanation of the concept.


Concatenate Fields in ArcMap 10

I added a field to an address layer in ArcMap 10 to concatenate separate address fields. I used the field calculator and a VB Script to accomplish this. So, for example, if the separate fields are "10" "Main" and "Street", then the new field would be "10 Main Street". The problem is that the concatenated field needs to update dynamically. So if the user changes "Main" in the database to "East Main" then that field will update in ArcMap. However, the concatenated field will not update to "10 East Main Street" unless I rerun the Field Calculator script. I need a way for this new field to always reflect the current data, but still do the concatenation. Does anyone have any suggestions? Thanks, Seth


Centralizing shapefile data into database

I've got hundreds of shapefiles from various different GIS projects that I want to start consolidating into a single database platform, currently attempting this with Postgres / PostGIS. Hardly any of the data is standardized - meaning that it's a lot of the same data _types_ , but the particular attribute names / types don't match. Where should I start tackling this? Should I develop a standard model to migrate each shapefile into first (e.g. Hydro_line, transport_line, Hydro_poly standards, etc)? An alternative is to just import each shapefile into Postgres individually, so each shp becomes a table in the database, but I'm not sure about this in terms of performance and organization. Feels sort of like delaying the inevitable... Any advice on dealing with this daunting task?


Autodesk MapGuide Enterprise MrSid on Linux

I am running MapGuide Enterprise 2011 on CentOS 5 linux. Using the Autodesk Raster FDO, i loaded a mrsid file (through the "Load file-based data" method). Whenever i try to view the layer with the mrsid file, the server crashes. most of the time it says "Segmentation Fault", but sometimes it just dies without saying anything. Here's what the error log says when i try to view the layer:               <2011-01-14T10:39:25>   53263248        Autodesk MapGuide Studio v2.2.0.5305       198.245.206.237 Administrator     Error: An exception occurred in FDO component.         Only MOSAIC(), CLIP() and RESAMPLE() are supported.     StackTrace:     - MgFeatureServiceHandler.ProcessOperation() line 83 file FeatureServiceHandler.cpp     - MgOpSelectFeaturesSpatial.Execute() line 110 file OpSelectFeaturesSpatial.cpp     - MgServerSelectFeatures.SelectFeatures() line 266 file ServerSelectFeatures.cpp      I've tried a few other mrsid files with the same problem. TIFF files load fine. Anyone has any insight on this? Thank you!


Problem when converting sosi to shapefile with ogr2ogr on Linux

I am trying to convert a SOSI file (Norwegian geodata format) to a shapefile. I have set up gdal, but apparantly there is a lot of things that I dont have support for. When I do a ./configure, this is the result:                 LIBTOOL support:           yes            LIBZ support:              external       LIBLZMA support:           no       GRASS support:             no       CFITSIO support:           no       PCRaster support:          internal       LIBPNG support:            internal       DDS support:               no       GTA support:               no       LIBTIFF support:           internal (BigTIFF=yes)       LIBGEOTIFF support:        internal       LIBJPEG support:           internal       12 bit JPEG:               yes       12 bit JPEG-in-TIFF:       yes       LIBGIF support:            internal       OGDI support:              no       HDF4 support:              no       HDF5 support:              no       NetCDF support:            no       Kakadu support:            no       JasPer support:            no       OpenJPEG support:          no       ECW support:               no       MrSID support:             no       MrSID/MG4 Lidar support:   no       MSG support:               no       GRIB support:              yes       EPSILON support:           no       WebP support:              no       cURL support (wms/wcs/...):no       PostgreSQL support:        no       MySQL support:             no       Ingres support:            no       Xerces-C support:          no       NAS support:               no       Expat support:             no       libxml2 support:           yes       Google libkml support:     no       ODBC support:              no       PGeo support:              no       FGDB support:              no       MDB support:               no       PCIDSK support:            internal       OCI support:               no       GEORASTER support:         no       SDE support:               no       Rasdaman support:          no       DODS support:              no       SQLite support:            no       PCRE support:              yes       SpatiaLite support:        no       DWGdirect support          no       INFORMIX DataBlade support:no       GEOS support:              no       Poppler support:           no       Podofo support:            no       OpenCL support:            no       Armadillo support:         no       FreeXL support:            no       SOSI support:              no                 SWIG Bindings:             no            Statically link PROJ.4:    no       enable OGR building:       yes       enable pthread support:    yes       enable POSIX iconv support:yes       hide internal symbols:     no      And doing the command results in this: ~/gdal$ ogr2ogr -f "ESRI Shapefile" /media//Geodata/testing/sosi.sos /Documents/test.shp FAILURE: Unable to open datasource `/Documents/test.shp' with the following drivers. I am fairly inexperienced with this, so I am trying to learn. What is going on here, what am I doing wrong?


Problems with ESRI JavaScript API

I'm having some problems to get started using the ESRI JavaSript API. I have a test map service published on my notebook but I can't get it to display. I don't have problems with esri hosted services, just the ones I'm trying to host. Here is some code:               <script src="http://serverapi.arcgisonline.com/jsapi/arcgis/?v=2.1" type="text/javascript"></script>     <link href="http://serverapi.arcgisonline.com/jsapi/arcgis/2.1/js/dojo/dijit/themes/claro/claro.css" rel="stylesheet" type="text/css" >     <script type="text/javascript">         dojo.require("esri.map");         dojo.require("esri.geometry");         dojo.require("esri.layers");         var map;         function init() {             esri.config.defaults.io.proxyUrl = "/proxy/proxy.ashx";             map = new esri.Map("mapDiv");             // var basemapURL = "http://server.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer";             var myURL = "http://brpnb010/arcgis/rest/services/teste/MapServer";                  // var basemap = new esri.layers.ArcGISTiledMapServiceLayer(basemapURL);             var myLayer = new esri.layers.ArcGISTiledMapServiceLayer(myURL);                  // map.addLayer(basemap);             map.addLayer(myLayer);                  dojo.connect(map, "onExtentChange", showExtent);         }              function showExtent(extent) {             var s = "";             s = "XMin: " + extent.xmin +             " YMin: " + extent.ymin +             " XMax: " + extent.xmax +             " YMax: " + extent.ymax;                  dojo.byId("onExtentChangeInfo").innerHTML = s;         }              dojo.addOnLoad(init);     </script>      I'm using Visual Studio to develop, and when the page loads, I get a "Could not load cross domain resources" error. **edit** I'm using version 2.1 of the Javascript API. And I'm developing on localhost, which has the brpnb010 name. Any ideas? Thanks, George


What is the difference between ESRI's SDE and Geospatial dbs?

When would you want to use ESRI SDE (I think it is called Server but I forget the level - please feel free to edit) versus a spatially enabled database. What are the trade-offs on either side? What are the benefits on either side?


What are the differences between spatial databases?

What is the difference between Oracle Spatial, PostGIS, SQLServer2008, SpatiaLite, and MySQL in terms of their spatial support. Are there specific use cases where each shines. What might you want to avoid with any of these databases.


Good conferences for Spatial/GIS developers

What are good conferences for GIS and spatial developers to attend. For ESRI desktop users I think the User Conference is a must attend but if you are an ESRI developer I think the ESRI developer conference is a better event. If you can't make the big shows then the regional events are also a really good deal. FOSS4G is great when it is in your neighborhood but that is a big if.


How can I print a bunch of Arcmap mxd's all at once?

I have a folder with a dozen Arcmap compositions (.mxd files) I need to print. They're all ready to go, I don't need to do any editing. Do I really have to open up each one, wait to refresh, and press print?


Can ArcGIS and OSGeo4W share the same python install?

I don't like unnecessary duplication. Is it possible to install python only once and have both ArcGIS and OSGeo4W use it?


Allow multiple users to share personal geodatabases safely?

We're using ArcGIS 9.3 but for reasons beyond the scope of this question we have to stick with personal geodatabases (MS Access .mdb's) and multiple users. How can we set this up so only a limited few have write access and still allow the bulk of the users to use it?


Calculating travel time from residential points to stores via public transit using ArcGIS for Desktop?

I don't know how to summarize it in just one title, but here is the problem: My goal is to calculate the travel time from residential points to stores via public transit. Basically this period of time consists three parts:   1. The walking time from a residential point to the nearest bus stop "t1", which can be calculated by the Nearest Facility tool from Network Analyst. The resulting table is called "Resident_Stop", of which the number of records is m;   2. The bus-riding time "t2" from this bus stop to the nearest bus stop of the target store via bus route, which is a record from Network Analyst's OD Time Matrix table, namely, "Stop_Stop". The origins and destinations are all bus stops on one single bus route so that the number of records is s*s;   3. The walking time from the ending bus stop to the store "t3", which is basically the same with the first step. The table is "Stop_Store", of which the number of records is n. So the total travel time "t" equals to (t1+t2+t3) and here is where problem occurs: I haven't figured out a way to attach the three tables reasonably to add the times on using a join. Currently I have an OD Time Matrix table of the walking time directly from residential points to stores via road network (no transit involved) and its number of records is m*n. I also separatly joined the table "Resident_Stop" in step 1 and the table "Stop_Store" in step 3 to the m*n-sized table based on residential point IDs and store IDs. So the remaining thing now is to assign the bus-riding time in table 2 based on table 1's stop IDs and table 3's stop IDs to each of the m*n records. Then I can get the total time. Should you get any ideas, I'm willing to discuss with you. I'm now thinking about using Python to handle this but I'm a beginner on that. If you have any books or tutorials to recommand please tell me.


provide an alternative "label" when styling vector layer with openlayers

I'm using openlayers 2.12 and I am styling a vector layer using your typical `new OpenLayers.Style({` Now, this layer does two things. First , renders polygons from Geoserver, so in the `style` I have `'label': '${p_name}'`. Works fine. The second thing the layer does is providing tools to the user to edit an already existing polygon. So when a user chooses to edit a polygon, (resize it for example) the familiar little points appear , as handlers. Now, these points/handlers render "undifined", as name because the `'label': '${p_name}'` is actually undifined. Did not came from Geoserver. Is there a workaround? How can I set an alternative name for the handlers? Thanks in advance


Determining projection of Texas county shapefile?

I have a Shapefile of San Jacinto County in Texas. It thinks it is NAD83 UTM Zone 15N, but the dataset isn't jiving with any other dataset. Max X: 1,073,721 Max Y: 11,223,970 Min X: 895,728 Min Y: 11,013,769 I don't usually work with datasets from this part of the world, so I don't know which projections are in common use. Anyone who is familiar with the area know what projection these extents might fit? This is the county I'm talking about: ![enter image description here](http://i.stack.imgur.com/E5iIh.png)


Setting default browser for QGIS actions

I am trying to open a pdf using QGIS action tool with Adobe Reader. Currently, the action defaults to Internet Explorer which doesn't succeed in opening my document. I'm currently able to open AutoCAD files and photos. Does anyone know how to set the default software used to by QGIS actions? I'm using QGIS version 2.0. my pathname in the attribute field is: explorer C:\Documents and Settings\user\Bureau\file.pdf Action Type: Windows. Name: Open file. Action: [%"pathname"%]. I didn't click the Capture checkbox. Thanks in advance!


How do I create a static popup in OpenLayers that contains HTML content?

I would like to create a static popup using openLayers 2.11 which is always in the lower-right corner of the screen and appears on map load. The idea behind this is to have a logo for my school always appear. I already have a marker layer setup. Maybe I am approaching this incorrectly. Any help would be appreciated. I have tried using some of the examples on the OpenLayers website regarding anchored popups with no luck. Thank you in advance from a poor, starving student. :) The contents of the popup are as follows:
         <link rel="stylesheet" type="text/css" href="css/cube.css">

front

back

top

bottom

left

right


How does a GPS/GSM/GPRS Tracker work without a monthly fee?

GPS newbie here. Looking at this product: http://dx.com/p/gps-gsm-gprs-tracker-for-personal- remote-positioning-29284 How does a GPS tracker work without a monthly plan? Thank You.


Does OSM2PO take into consideration turn restrictions?

It is supposed that for using Shooting star function in pgrouting, we need to have RULE and TO_COST fields in our table, but OSM2PO hasn't created those fields. what must i do for OSM2PO takes into consideration the turn restrictions? Thank you! (sorry, my english is so poor...)


Info from the getResultImageLayer missing

I was finally able to get the result layer overlaid on my application but now I am facing another problem. The output of the layer url does not show any field or attribute details. I am not sure why this is happening? The output looks like this, ![enter image description here](http://i.stack.imgur.com/thvkY.png) As you can see that it is only describing the Extent and other Map properties but the attribute details, which are visible in the jobs folder of the newly created final_shp, is missing on the url and so I am not able to fetch anything onto the map. Any Idea how I can resolve this?


Creating tiles in Mapnik

I have installed mapnik and successfully tested it. Now I am facing problem to create my own customized tiles. I have google it but everywhere it was telling about osm tile generation. My vision is simple. I have a huge map data of my target area. I want to create my own customize map in Mapnik/tilemill. Then create tiles and upload them into server. A step by step instruction will be helpful to me as I am newbie in Mapnik. This my main problem. How will I create tiles server of my own data??? Thanks in advance


How do i get x, y coordinate values when click on map?

I am very preliminary in arc gis and in fact i am a web programmer. I created a local serviece and use this to create a web applcation in ArcGIS > ArcGIS Server for the Microsoft .NET Framework > ArcGIS Server Manager.exe . Then i copy website from "c:/inetpup/wwwroot/projectName" and open this to VisualStudio. Now i want when i click on map, display an pushpin on click location and send x, y to server.


How to add point to a polyline?

I have a `polyline`, and i want to add `point (vertex)` on a specific distance, like two meters for example. How can i do that?


How can I make the layout size 36" x 72" in ArcGIS?

I need to make a poster to present at a conference and I'm wondering if there's anyway to do this all in ArcGIS 10.1. Is it possible to change the layout/page size to 36x72 inches? Is it practical or even feasible to do this? Thanks! Andy


Why is arcpy.mapping.ConvertWebMapToMapDocument slow?

We use arcpy.mapping for print tasks. The slowest part in workflow is arcpy.mapping function **ConvertWebMapToMapDocument** ; it takes 10-12 seconds each time, doesn't matter if map is big or not. We use ConvertWebMapToMapDocument(mapJson, templateMxd) form. templateMxd is empty A3 document with empty data frame and north arrow. Function ExportToPNG takes only 2 seconds while ConvertWebMapToMapDocument eats most of the total printing time. What can be tips to improve performance here? **Edit: ------------------**   Looks it's not possible to solve it for myself as ConvertWebMapToMapDocument implementation spends a lot of time in unnecessary web requests.   Simplest example:               import time     import arcpy          templateMxd = 'C:\\Program Files\\ArcGIS\\Server\\Templates\\ExportWebMapTemplates\\A3 Portrait.mxd'     mapJSON = '{"layoutOptions": {"titleText": "Simple WebMap JSON example"},"operationalLayers": [{"url": "http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer","visibility": true}],"exportOptions": {"outputSize": [1500,1500]},"mapOptions": {"extent": {"xmin": -13077000,"ymin": 4031000,"xmax": -13023000,"ymax": 4053000}},"version": "1.4"}'          timer = time.time()     aTuple = arcpy.mapping.ConvertWebMapToMapDocument(mapJSON, templateMxd)     print('ConvertWebMapToMapDocument took %.2f s' % (time.time() - timer))           Output:               ConvertWebMapToMapDocument took 2.10 s      Now where all this time went? I used Fidler2 to track web requests and got following:               http://maps1.arcgisonline.com/ArcGIS/rest/info?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer/layers?&f;=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer/legend?&f;=json     http://maps1.arcgisonline.com/ArcGIS/rest/info?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/info?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer/layers?&f;=json     http://maps1.arcgisonline.com/ArcGIS/rest/info?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/info?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer/layers?&f;=json          Aggregate Session duration: 00:00:01.687      Only first 4 requests are actually needed other requests are duplicates. For these 4 requests               Aggregate Session duration:    00:00:00.683 Such optimization would save 1 second of 2.10 s spent in ConvertWebMapToMapDocument function.


Is there a way to clip a raster using a vector polygon without having Spatial Analyst?

I have ArcView 10 and I'm trying to clip a DEM (TIFF) to a wilderness boundary (polygon). I know ArcView won't let me do this without Spatial Analyst, but I don't have $2500. I tried exporting my raster and vector in QGIS, but for the life of me, I can't get the the two to line up (different coord systems). I've tried to "DEFINE COORD SYSTEM" in Arc Toolbox, but something just isn't right... I'm about to go insane after 5 hours of failures... Can anyone alleviate this pain??? PLEASE


Editing a KML file or converting KML to Editable format with attributes

I'm using QGIS. I'm trying to either edit a .kml file someone else created in Google Earth or convert it to another editable format (like a .shp file). I've read the following threads without really being able to figure out what the solution is: Any luck with qgis2kml? (seems like this plugin doesn't exist? couldn't find it) How to convert KML to shape file without losing attributes? (I tried the zonums converter mentioned in the answers section, but I lose most of the data in the original .kml file) Converting kml to shape file via ogr2ogr including all attributes (I'm enough of a newbie that I didn't really understand this thread) I've also tried the kml to csv converter here http://www.gps-data- team.com/convert.php. Again, it strips most of the data, giving me only the points and names. Edited to add: saving my .kml file as a .shp file does not work. I get the following error: Export to vector file failed. Error: creation of data source failed (OGR error:Failed to open shapefile .shp/.shp. It may be corrupt or read-only file accessed in update mode. ) Here is a snippet from the beginning of the file. I'm guessing by the comments below, it's just really poorly organized. You have these pieces with the pushpin markers:                   <?xml version="1.0" encoding="UTF-8"?>     <kml xmlns="http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2" xmlns:kml="http://www.opengis.net/kml/2.2" xmlns:atom="http://www.w3.org/2005/Atom">     <Document>     <name>Libanius' Letters, Cribiore.kml</name>     <open>1</open>     <Style id="sn_ylw-pushpin00">     </Style>     <Style id="sh_ylw-pushpin2">         <IconStyle>             <scale>1.2</scale>         </IconStyle>     </Style>     <StyleMap id="msn_grn-pushpin">         <Pair>             <key>normal</key>             <styleUrl>#sn_grn-pushpin0</styleUrl>         </Pair>         <Pair>             <key>highlight</key>             <styleUrl>#sh_grn-pushpin</styleUrl>         </Pair>     </StyleMap>     <Style id="sn_wht-pushpin">         <IconStyle>             <scale>1.1</scale>             <Icon>                 <href>http://maps.google.com/mapfiles/kml/pushpin/wht-pushpin.png</href>             </Icon>             <hotSpot x="20" y="2" xunits="pixels" yunits="pixels"/>         </IconStyle>     </Style>     <Style id="sn_grn-pushpin0">         <IconStyle>             <scale>1.1</scale>             <Icon>                 <href>http://maps.google.com/mapfiles/kml/pushpin/grn-pushpin.png</href>             </Icon>             <hotSpot x="20" y="2" xunits="pixels" yunits="pixels"/>         </IconStyle>     </Style>      .... Then about halfway through, there's a folder tag that contains the place marks :                   <Folder>         <name>Libanius' Letters, Cribiore</name>         <description>Exported on Mon Oct 01 2012 16:11:04 GMT-0700 (PDT)</description>         <Placemark id="18">             <name>Antioch</name>             <description>Home</description>             <styleUrl>#msn_star</styleUrl>             <ExtendedData>                 <Data name="searchrename">                     <value>Antiochia ad Orontem</value>                 </Data>                 <Data name="en_name">                     <value>Theoupolis</value>                 </Data>                 <Data name="featuretyp">                     <value>settlement</value>                 </Data>                 <Data name="pid">                     <value>658381</value>                 </Data>                 <Data name="gr_name">                     <value>Ἀντιόχεια ἡ Μεγάλη</value>                 </Data>                 <Data name="la_name">                     <value>Theoupolis</value>                 </Data>                 <Data name="timeperiod">                     <value>AHRL</value>                 </Data>                 <Data name="perseus_li">                     <value>http://www.perseus.tufts.edu/hopper/text?doc=Perseus%3Atext%3A1999.04.0006%3Aalp</value>                 </Data>                 <Data name="wiki_link">                     <value>http://en.wikipedia.org/wiki/Antioch</value>                 </Data>                 <Data name="path">                     <value>/places/658381</value>                 </Data>                 <Data name="x1">                     <value>36.15818</value>                 </Data>                 <Data name="y1">                     <value>36.202097</value>                 </Data>             </ExtendedData>             <Point>                 <coordinates>36.15818,36.202097,0</coordinates>             </Point>      ..... I tried deleting all the pushpin crap in Oxygen and then reopening the layer in QGIS. The new kml file opens, and the points are there with the data I want. Still not editable, and when I try to save as shapefile I get: "Export to vector file failed. Error: creation of layer failed (OGR error:Failed to open Shapefile `LibaniusShapefile.shp'. ) I think I'm not sure what a well-formed KML file should look like. Thanks.


Why does Kernel Density fail in model with ERROR 000865?

I am attempting to build a model using variables for automation. (see below) ![Model](http://i.stack.imgur.com/HkYDk.jpg) The problem I am having is when the model attempts to run the Kernel Density. It seems as if the Kernel Density will not recognize %FilePreFix%_ts_cd even though a shapefile is created. Below is the error I get. ![Error](http://i.stack.imgur.com/R7Lcg.jpg)


TileMill raster-colorizer

I am trying to use the raster-colorizer in Tilemill (v0.10.1) but it is giving me an error. I did set in the advanced section the band to `band="1"`. What else do I have to do in order to use raster-colorizer? ![enter image description here](http://i.stack.imgur.com/gvy3G.png)


ArcMap lagging when I use the measure tool

ArcMap is lagging when I click the measure tool, my data isn't particularly big (biggest .shp is ~150mb) and my PC is brand new. Any idea what may be causing this?


Raster Analysis - Identifying features from Raster

I have an IDW, tension spline and hill-shade rasters that show an archaeological landscape and from these rasters I want to identify new features which are seen in the landscape as depressions. And from this raster I want to calculate the size and shape of these features. I have not yet come across a form of analysis that would do this, or I am just not understanding the different types of analysis I have looked at. I have been using GRASS, QGis, ArcMap/Secne and surfer, but I'm willing to give anything a try. Thanks Christina


How can I select which interpolation technique will be used in a gdal command? (gdal_fillnodata.py)

I work on Landsat 7 slc-off images and I would like to use the gdal_fillnodata.py command to fill-in the missing pixels. The problem that I have is that the gdal_fillnodata.py uses the "inverse distance weighting" (IDW) interpolation technique and I want to use the kriging geostatistical technique. Do you have any idea how can I do it? Thank you!


How to use IDatabaseCompact.Compact() on an used IWorkspace, without permission denied

I am writting an ArcMap-AddIn with vb.net. I get an error when trying to `compact` my IWorkspace(mdb), the error is: > COMException > > You tried to open a database, which was already opened by the user 'ADMIN' > on Computer 'XXXXXX'. Try it again when the database is available. > > On ESRI.ArcGIS.Geodatabase.IDatabaseCompact.Compact() on > MyProject.MyClass.CompactGDB(IWorkspace pWS) How can i access that workspace, which i used before, to compact it? The Code:               ' CompactGDB         Public Sub CompactGDB(ByVal pWS As IWorkspace)                  Dim pDatabaseCompact As IDatabaseCompact             If (TypeOf pWS Is IDatabaseCompact) Then                 pDatabaseCompact = CType(pWS, IDatabaseCompact)                 If (pDatabaseCompact.CanCompact) Then                     Try                         pDatabaseCompact.Compact()                     Catch ex As Exception                         MessageBox.Show(ex.type & ex.Message & ex.StackTrace, "Error", MessageBoxButtons.OK, MessageBoxIcon.Error)                     End Try                      End If             End If              End Sub      I've got 8 other functions, which work with my workspace. Any suggestions?


Localized Linear Histogram Match using python

I would like to use the local linear histogram match for filling the missing pixels of a Landsat 7 image with another older image (SLC Gap-Filled Products). Do you have any idea how can I do it by using the gdal commands? Thank you!


What kind of different results or what is the difference between two URL based WMS requests

I am trying to retrieve some WMS layers from geoserver so that I can display them in leaflet and I need to populate a URL field and I noticed that many of the layers that I opened up in openlayers had "GetMap&layers" vs examples that were successful with "GetFeatureInfo&Layers" when being served via URL.


What is this PostGIS query doing to show great circle connections?

I'm trying to draw a flight map in CartoDB. My starting point is a single lat long point. We'll pretend it is this:               ST_GeomFromText('POINT(-71.060316 48.432044)', 4326)      And I want draw lines between that point and the points in my CartoDB table. I had a hard time even getting started, but someone pointed me to this map: http://bl.ocks.org/andrewxhill/raw/8406976/ Which is fancier than what I need (I don't need to redraw everytime you click someplace) but includes the basic query:               SELECT ST_Transform(ST_Segmentize(ST_MakeLine(        ST_Transform(the_geom, 953027),        ST_Transform(               (SELECT ST_PointOnSurface(the_geom)                FROM lower_48_zips                WHERE cartodb_id = "+data.cartodb_id+"), 953027)), 100000), 4326)         the_geom FROM us_ski_areas WHERE mindist < 1000      Someone also pointed me towards an example which uses a very similar query:               UPDATE experimental.airroutes     SET the_geom =     (SELECT ST_Transform(ST_Segmentize(ST_MakeLine(            ST_Transform(a.the_geom, 953027),            ST_Transform(b.the_geom, 953027)          ), 100000 ), 4326 )      FROM experimental.airports a, experimental.airports b      WHERE a.id = airroutes.source_id         AND b.id = airroutes.dest_id     );      But I'm having trouble pulling either of these apart. I'm decently conversant in MySQL but totally new to PostGIS and I would love help taking apart what is happening here and replacing source_id with a specific point. I thought I could do something like               UPDATE my_table     SET the_geom =     (SELECT ST_Transform(ST_Segmentize(ST_MakeLine(            ST_Transform(my_table.the_geom, 953027),            ST_Transform(ST_GeomFromText('POINT(-71.060316 48.432044)', 4326), 953027)          ), 100000 ), 4326 )      FROM my_table     );      But I know I need a "where" clause in there and I'm not even sure what it should be.


Stitching MapProxy tiles to receive one georeferenced image?

In my struggle to get a decent aerial image for playing around with as a layer in Biosphere3D (and for other purposes as well), I managed to install Mapproxy and then cached the area in my desired zoom level from a 2m-Orthophoto WMS through the proxy within QGIS.. Well, now I am in possession of tiles worth 60MB, but i see no way how to stitch them together to one image. The tiles are stored in a folder-structure like this: ![enter image description here](http://i.stack.imgur.com/eRGsK.png) The following sets of tiles are in 32/000/781/000/25, 32/000/781/000/25 etc. containing 408.png to 431.png I have no clue how to interpret this structure. I wonder if there is any way (script) to stitch those tiles to receive one georeferenced image. The tiles were captured in WGS84. Grateful for any hint on how to proceed.


How can I determine the location of a visitor to my website?

If I want to provide a website with location appropriate content, what ways can I use to find the visitor's rough location using their IP or header information?


pgrouting download link not working

I am trying to download pgrouting but the download links are not working. I have tried many links but no link is working. can you please help me in this regard. will be thankful


Create Materialized View to Flatten ArcGIS JavaScript API Database

I'm doing a project using ArcGIS JavaScript API v3.5 hosted on ArcServer v10.1. While reading the question below, remember that I **do not** have control over the datasource and cannot just create a materialized view in the database to solve the issue. **Problem:** The issue I'm having is that parts of the database I'm querying are very very normalized -- one table has 11 Foreign Keys that point to parent tables that are even more normalized. You can imagine that to join all the information for a single record for all these relationships is cumbersome. It requires many requests and lots of looping and checking which is slow and error prone. **Questions:** I have 2 questions:   1. What other options are there in the JS API to make querying easier from the client? I see layer types like QueryDataSource in the API reference that seem to be what i want. But there are no examples of how that works and no one online seems to be using that solution. Has anyone here used QueryDataSource?    2. I'm open to options on how to request that the datasources be modified. The most simple solution is to make them create a flattened table joined with all the data. What other options are there? Since I don't know the ArcGIS Server aspect as well I'm a little blind in what to request. For example, is there a way in ArcMap to create materialized-like views for a normalized datasource? Another possible solution is that the API should gather all related records up for me and spew out one big JSON object when i request a child table with many Foreign Keys. This is how I create API's in Django. All related records become nested hashes of a record. Is there someway to construct your API like this with ArcServer? I can't wait to hear what thoughts/solutions others have for this problem


Displaying DEM raster in QGIS

Using Raster>Analysis>DEM to display hillshade produces the message "gdaldem.exe had stopped working" so I asked the community. I'm new to QGIS, and I now realise after a little experimentation that I should have used Raster>Terrain Analysis>Hillshade which now works OK - just in case another newbie has encountered the same problem. Why the former action should crash gdaldem.exe remains a slight concern though.


How to convert Shapefile to TopoJSON with grouping neighbor shapes by some criteria

I have a shapefile with countries, which I convert to GeoJSON using `ogr2ogr` and then converting that json to another TopoJSON, using `topojson` command- line tool. Everything works like charm, but now I want to learn to concatenate/merge/join/cluster (just don't know what is the most appropriate term in GIS context) different neighbor shapes to a bigger ones - say, all African countries to one big Africa shape. Let's assume that all Hdata needed (like continent name) is present in the original shapefile. So, what is the simplest way to merge neighbour shapes while converting shapefile to topojson format? Sorry if this question is about something trivial, it's just that this is a new topic to me, though I'm very glad I should deal with such issues.


Separate polygons based on intersection using PostGIS

I have a PostGIS table of polygons where some intersect with one another. This is what I'm trying to do:   * For a given polygon selected by id, give me all of the polygons that intersect. Basically, `select the_geom from the_table where ST_Intersects(the_geom, (select the_geom from the_table where source_id = '123'))`   * From these polygons, I need to create a new polygons such that intersection becomes a new polygon. So if polygon A intersects with polygon B, I will get 3 new polygons: A minus AB, AB, and B minus AB. Any ideas?


The mysteries of the ArcGIS Server Identity - can you help?

Just getting back into ArcGIS for .NET again after about 2 years of doing the training course. I have a question about the "Add ArcGIS Identity" that you choose in Visual Studio - it doesn't seem to matter what combination of username/password I think of, it always says I have the an invalid username/password (by the way I am using accounts that are in the agsadmin group) on the AGS server. Yet if I copy a ciphered identity tag from the web.config of some random AGS application I made via ArcGIS Manager - it works!!, i.e. the web app loads up the map services and I can see a map in my application. Now, I can't decipher this identity tag as you can probably guess - but does anyone know what I might be doing wrong or what I need to look, or check for in order to get the correct credentials right. I just can't rest until I have this figured out - instead of doing this hack. Cheers **UPDATE - Answer** I made sure my dev PC had the ArcGISWebServices user set with the same password as the AGS physical server - so the "impersonation" was consistant between the two machines and accounts. I replaced out the ciphered details in the web.config and replaced with the ArcGISWebServices one and it works fine. I also tried doing it via the "Add ArcGIS Identity" and now it works but make sure you put the machine name (i.e. the name of your dev pc) in the "Domain or Machine Name" textbox.**


how to retrieve ATMs List from OSM

I am using this command to get the ATMs list `osmosis --read-xml my_map.osm --tf accept-nodes amenity=atm --write-xml output-nodes.osm` But it is giving the list of nodes and ways also...I want only nodes and is there any relation for these nodes and ways.


What is the best way to protect data online?

What techniques do people use when protecting data ( _e.g._ , preventing a user downloading the whole dataset) they are overlaying on top of Google/Yahoo/Bing etc.?


PostGIS importing shapefile error w/ official manager

I'm having some trouble when trying to import an shapefile of the Sao Paulo Brazil's state (I downloaded it from the web and I'm using the own PostGis manager to import it). Here's an image! ![enter image description here](http://i.stack.imgur.com/vGZQv.jpg) The connection is already on. In the "Options" menu, i had to change from "UTF8" to LATIN1, cause if I didn't, I'd get this error message: "Unable to convert field name to UTF-8 (iconv reports "Illegal byte sequence"). Current encoding is "UTF-8". Try "LATIN1" (Western European), or one of the values described at http://www.gnu.org/software/libiconv/." I also wanted to now what is SRID and wich value should i put there. When I changed to LATIN1, I get this other error message: "Failed in pgui_exec(): ERRO: sequência de bytes é inválida para codificação "UTF8": 0xe3 0xa7 0xe3". In english, it says the byte sequence is not valid to UTF8 codification. I'm a bit desperate cause I've never used GIS before and I'll need to deliver a project using it on the end of the year. Thank you guys since now.


Bing maps get Distance between two places

I have referred the sample http://help.arcgis.com/en/webapi/silverlight/samples/start.htm#BingRouting and developed a Bing maps routing tool on my Silverlight GIS website.. now my client wants me to also show the distance between two places whos direction is shown.. I took a look at the Routing class but failed to find anything to get the distance.. I used the ‘haversine’ formula, but it does not give me the exact distance along the route.. or Can I calculate the length of the geometry/polyline drawn on the map.. Will this result be accurate? Please guide me so that I can get the distance between my start and destination..


Comparison of SDI frameworks

What are the main differences between these SDI, especially in the level of geospatial features and ease of integration with existing CMS?   * Cartaro   * Geonode   * GeOrchestra   * EasySDI I wonder for the reasons why SDI software is not so popular?


How to restrict portions of road network from the analysis in ArcGIS 10 Network Analyst?

I'm trying to sort out network analysis in ArcGIS 10. Toy example of finding shortest distance between two points looks something like that : ![closest facility output](http://i.stack.imgur.com/xpDPM.jpg) Where squares are two locations, and the green line is an output from Closest Facility tool. What I need to do, is to restrict motorways and slip roads from the analysis. So the new route should avouid the red segments: ![restricted](http://i.stack.imgur.com/gRUas.jpg) What is the best strategy to proceed in such case? Should I edit my network and delete all the segments I do not want to have in analysis? Or perhaps use linear Barriers? Thanks for help!


Is there a st_buffer inverse function, that returns a width estimation?

When you supply a geometry `g1` and a width `w` to the standard OGC function ST_Buffer, it returns a new geometry `g2`,                 g2 = ST_Buffer(g1,w)      This question is about the inverse: a function where you can supply two related geometries -- `g1` and `g2` satisfies ST_Within(g1,g2) -- and returns the width parameter,                 w = Unbuffer(g1,g2)      Now the interpretation of `w` is "average width". So, **where is this function _Unbuffer_ defined? Are there references for it? Is there a standard library for it?** * * * Complementing after bounty. Examples of applications (estimation of): street width, sidewalk width, river width, canal width, internal polygon semelhance, riparian vegetation buffer width, etc. Example of proposals of credible sources: Pavelsky & Smith 2008 ; ArcMap Script Procedure. They are **not** solutions here because are (very) CPU consumer, "river-specialized", and _can not be implemented by OGC standards_. They calculate average of _river widthes_ by each of many perpendicular (of a central line) lines; so, they are "brute force" algorithms. Example of proposal _without_ credible and/or official sources: troubleshooter here. In this case, you can get the bounty finding sources that makes credible the posted solution.


Why can't I install python-qgis on Ubuntu 12?

I am with ubuntu 12 and Qgis 1.8. I have installed QGIS and grass plugin but I can't install python plugin. How to install python plugin with dependeny in Ubuntu?               ~$ sudo apt-get install python-qgis     [sudo] password for sugandha:      Reading package lists... Done     Building dependency tree            Reading state information... Done     Some packages could not be installed. This may mean that you have     requested an impossible situation or if you are using the unstable     distribution that some required packages have not yet been created     or been moved out of Incoming.     The following information may help to resolve the situation:          The following packages have unmet dependencies:      python-qgis : Depends: python-qt4 (>= 4.1.0) but it is not going to be installed     E: Unable to correct problems, you have held broken packages.


tin model in Qgis

I am having DSM and DTM data .My primary objective is to develope a CHM(Canopy Height Model) in Qgis.which i have done already but in my company they are stressing that CHM model which i already achieved is in Raster but they want the result in TIN model how to do this in Qgis?


"XY to Line" not converting all 58,000 records

I have a dbf (and xls) that contains ~58,000 records of latitude and longitude for both start and end points. I would like to connect each XY coordinate pair using the XY to Line tool found inside the Data Management toolbox. I will be using the Great Circle method fyi. The resultant shapefile contains ~16,000 records/polylines. As far as I can tell, there is no difference between processing the table inside or outside a file geodatabase. Anyone have any thoughts as to why it may not be processing the entire table? This issue persists with both the dBaseIV and XLS file. Cheers, Michael


Getting Field value for a given coordinate out of ESRI Shapefile using GDAL/OGR

I have a bunch of ESRI shapefiles that contain Land Cover data across Canada (provided by Geobase). I need to query the files and obtain the value of a certain field (called COVTYPE) in a given location specified either by lat/lon or UTM coordinates. I can open a file in Quantum GIS and find the needed values manually, but this process needs to be done automatically for a large number of values so I'm using C++ and GDAL. I'm very new to the whole GIS area and I'm just starting to learn GDAL. So far I've learned how to loop over layers, features, fields etc, but I can't figure out how to relate coordinates to features. I went through the official GDAL/OGR tutorials, but I'm still stuck. Any pointers would be appreciated! **Update:** Just to clarify. My input data would be coordinates (lat/lon or UTM). So the way I imagine this should work is:   1. I find which feature the given point is located within   2. read the field value in that feature I know how to do step 2, but I'm stuck at step 1. **Final solution** Technically, the update to Ragi's answer comes very close to what I need to do, but I'll post here a quick summary of how I ended up using it. My way is easier for constructing a point geometry object. So leaving all the boilerplate code out, the essential part is:               OGRPoint* point = new OGRPoint(some_longitude, some_latitude);     layer->SetSpatialFilter(point); //getting only the feature intercepting the point     OGRFeature* feature = layer->GetNextFeature();     int value = feature->GetFieldAsInteger(0); //getting the value - DONE!     //in my case it's 1st field, hence the 0th index


My OpenLayers map was working fine before, showing a KML overlay. Suddenly it broke. Now I see this error in my browser error console. > TypeError: 'undefined' is not a constructor (evaluating 'new > OpenLayers.Layer.GML') My code looks like this:               var kmllayer = new OpenLayers.Layer.GML("KML", "./my-kml-file.kml",         {             format: OpenLayers.Format.KML,          formatOptions: {           extractStyles: true,            extractAttributes: true,           maxDepth: 2         }        });     map.addLayer(kmllayer);


Can I split up v.in.ogr into multiple v.clean commands?

I'm trying to load a large set of polygons into GRASS with " _v.in.ogr_ " and am running out of memory at the "Break Boundaries" stage. I'm wondering if I can break into steps all of the different things that _v.in.ogr_ does to a data set when it is run without the " _-c_ " flag? For example, I know that it must run the equivalent of _"v.clean tool=bpol_ " and " _v.clean tool=rmdupl_ ", but I'm not exactly sure what to run after those 2 to complete the polygon cleaning that normally happens. I figure that I might be able to save some memory if I don't have to run all the operations in one go. Does this makes sense? Thanks. * * * Ok, so it was like digging for gold to find this, but in a presentation given by some people who seem to know GRASS pretty well, they said that the following steps were needed to clean polygon geometry after import: **Cleaning procedures for area import**   * break polygons   * remove duplicates   * break boundaries   * remove duplicates   * clean boundaries at nodes   * change dangles to lines   * remove bridges Link to presentation here: http://geostat- course.org/system/files/GRASS_geostat_landau_2011_intro.pdf * * * Wait, no... there's more to this apparently. I tried those steps listed above (or attempted to at least). I think I'm worse off now than when I started.               v.in.ogr -c     Importing map 52 features...     Number of boundaries: 52     Number of centroids: 52     Number of areas: 42     Number of isles: 42     Number of incorrect boundaries: 10     Number of centroids outside area: 10          1. break polygons       v.clean type=boundary tool=bpol     Number of boundaries: 270     Number of centroids: 52     Number of areas: 2     Number of isles: 4     Number of incorrect boundaries: 268     Number of centroids outside area: 50          2. remove duplicates       v.clean type=boundary tool=rmdupl     Number of boundaries: 190     Number of centroids: 52     Number of areas: 40     Number of isles: 4     Number of incorrect boundaries: 126     Number of centroids outside area: 12          3. break boundaries       v.clean type=boundary tool=break     Number of boundaries: 5147     Number of centroids: 52     Number of areas: 151     Number of isles: 0     Number of incorrect boundaries: 5132     Number of centroids outside area: 40     Number of areas without centroid: 139          4. remove duplicates       v.clean type=boundary tool=rmdupl      Number of boundaries: 4803     Number of centroids: 52     Number of areas: 153     Number of isles: 0     Number of incorrect boundaries: 4781     Number of centroids outside area: 38     Number of areas without centroid: 139          5. clean boundaries at nodes       (No idea what-the-flip this means.  Clean with what?  Do it by hand?)          6. change dangles to lines        v.clean type=boundary tool=chdangle thresh=-1     Number of boundaries: 4803     Number of centroids: 52     Number of areas: 153     Number of isles: 0     Number of incorrect boundaries: 4781     Number of centroids outside area: 38     Number of areas without centroid: 139          7. remove bridges       v.clean type=boundary tool=rmbridge in=eq1_6 out=eq1_7     Number of boundaries: 4420     Number of centroids: 52     Number of areas: 153     Number of isles: 22     Number of incorrect boundaries: 4352     Number of centroids outside area: 38     Number of areas without centroid: 139      Guess I'm off to read the source for _v.in.ogr_ now.


ArcSDE 9.3.1 -> Performance

We have a number of mapservices in ArcGIS Server 10.0 that reads data from FGDBs. We are now moving that data in to ArcSDE and the plan is to use this data from AGS. Is it better to connect to SDE using Direct Connect or going through the ArcSDE services performance wise?


Saga Reclassify grid problem

I am using Qgis 2.1 ver. 5f5cd4c (64 bit version on win7 64) OSGeo installed. I was trying to reclassify a DEM in EPSG:32614. I used the SAGA's Reclassify algorithm in the processing toolbox. I defined that I wanted to use a simple table, defined the table and that I wanted the reclassified grid saved to a tif file. I also wanted the resulted reclassified DEM to be opened after the algorithm finished the calculations. This is what I am repeatedly getting: ![enter image description here](http://i.stack.imgur.com/nulBb.png) I tried to open a new project and load the resulting raster manually, but Qgis does not recognize it as a valid raster data source: ![enter image description here](http://i.stack.imgur.com/EWNl0.png) I was able to open the file and convert it to GeoTIFF in SAGA. I wonder if I am missing something in the reclassifying procedure, or something related with the format of the resulting raster. Thanks in advance for any advice. PS. I was able to use r.reclass to do what I wanted. I also tested the saga reclasify algorithm in other two installations and I have the same problem. Gerardo


Geoserver PrintingProvider crushes after font color change I do to config.yaml

I've installed OpenGeo 3.1 and I want to do some changes in config.yaml. The thing is that after certain changes like changing font color or adding font type it cause the print provider to crush. for example : In A4:               - !text        text: "${mapTitle}"        fontSize: 30        spacingAfter: 30        align: center      if I add `font: Helvetica` or add `fontColor:black` I get error from `http://localhost:8080/geoserver/pdf/info.json` **This makes me think that the printing provider i'm using doesn't support these configurations** I've checked in geoserver , These 2 jars are installed :               printing-2.2-SNAPSHOT.jar     print-lib-1.2-20101018133712.jar      What am I missing here ?!               HTTP ERROR: 500          org.springframework.web.util.NestedServletException: Request processing failed; nested exception is org.ho.yaml.exception.YamlException: Error near line 0: End of document expected.          RequestURI=/geoserver/pdf/info.json     Caused by:          java.lang.RuntimeException: org.springframework.web.util.NestedServletException: Request processing failed; nested exception is org.ho.yaml.exception.YamlException: Error near line 0: End of document expected.         at org.geoserver.monitor.MonitorFilter.doFilter(MonitorFilter.java:177)         at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:70)         at org.geoserver.filters.SpringDelegatingFilter.doFilter(SpringDelegatingFilter.java:45)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.geoserver.platform.AdvancedDispatchFilter.doFilter(AdvancedDispatchFilter.java:49)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.vfny.geoserver.filters.SetCharacterEncodingFilter.doFilter(SetCharacterEncodingFilter.java:109)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:311)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.invoke(FilterSecurityInterceptor.java:116)         at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.doFilter(FilterSecurityInterceptor.java:83)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.access.ExceptionTranslationFilter.doFilter(ExceptionTranslationFilter.java:113)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerAnonymousAuthenticationFilter.doFilter(GeoServerAnonymousAuthenticationFilter.java:53)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter      This is the originl config.yaml :               dpis: [75, 150, 300]          scales:         - 524288000         - 262144000         - 131072000         - 65536000         - 32768000         - 16384000         - 8192000         - 4096000         - 2048000         - 1024000         - 512000         - 256000         - 128000         - 64000         - 32000         - 16000         - 8000         - 4000         - 2000         - 1000          hosts:         - !localMatch             dummy: true         - !ipMatch             ip: 0.0.0.0             mask: 0.0.0.0          layouts:              A4:             metaData: &commonMetaData                 title: "${mapTitle}"                 author: "GeoExplorer"                 subject: "Map printed with the OpenGeo Suite"                 keywords: "map,print"                 creator: "OpenGeo Suite"             mainPage:                 pageSize: A4                 rotation: true                 items:                     - !text                         text: "${mapTitle}"                         fontSize: 30                         spacingAfter: 30                         align: center                     - !map                         spacingAfter: 10                         width: 440                         height: 600                     - !scalebar                         type: bar                         maxSize: 100                         barBgColor: white                         fontSize: 8                         align: right                         spacingAfter: 10                     - !text                         text: "${comment}"              Legal:             metaData: *commonMetaData             mainPage:                 pageSize: LEGAL                 rotation: true                 items:                     - !text                         text: "${mapTitle}"                         fontSize: 30                         spacingAfter: 30                         align: center                     - !map                         spacingAfter: 10                         width: 440                         height: 650                     - !scalebar                         type: bar                         maxSize: 100                         barBgColor: white                         fontSize: 8                         align: right                         spacingAfter: 10                     - !text                         text: "${comment}"              Letter:             metaData: *commonMetaData             mainPage:                 pageSize: LETTER                 rotation: true                 items:                     - !text                         text: "${mapTitle}"                         fontSize: 30                         spacingAfter: 30                         align: center                     - !map                         spacingAfter: 10                         width: 440                         height: 550                     - !scalebar                         type: bar                         maxSize: 100                         barBgColor: white                         fontSize: 8                         align: right                         spacingAfter: 10                     - !text                         text: "${comment}"


Export georeferenced image into Google Earth

I have created a georeferenced image as instructed in this fine tutorial: http://qgis.spatialthoughts.com/2012/02/tutorial-georeferencing-topo- sheets.html Now I would like to get this image into Google Earth, as an Image Overlay or such. Any tips?


Does the Bing Maps API have a distance to coast service?

Does the Bing API contain a distance to coast service from an address? All we need to do is calculate the distance to a major coastline (Atlantic or Gulf of Mexico).


how to combine GPS track logs using modelbuilder?

I'm creating a model to import GPS tracks(.GPX files) into ArcMap, this is to create a network of driveable off road vehicle routes. I've manage to do some of the things I'm trying to do, but there are some problems I haven't been able to figure out, and I was hoping that some of you might be able to help! I'm using ArcGIS desktop 10.2. what i want "out" is: 1) generalized lines from the different tracks( Integrate tool) 2) number of tracks combined to make one line segment (with integrate) 3) average speed on the line segments 4) average time(seconds) it takes to drive the line segment 5) last date driven 6) graphic representation of how many times a line segment is driven So far I've made two separate models, simply because I can't figure out how to combine them. model one: ![enter image description here](http://i.stack.imgur.com/i3pVd.png) This model takes GPX files from the folder "Tracklogs", "GPX to feature" imports the gps points with all its data. "Convert time field" converts a dateTime field into a Date field. "GPX speed" creates a polyline containing speed values for each line segment(line segments between gps points), "Split line at vertices" so I can use "join field" to get both time and speed info togehter. "Add field" and "calculate field" to calulate speed and distance to get seconds it takes to drive on line segment. model two: ![enter image description here](http://i.stack.imgur.com/ajCgo.png) this one merges all files from the first model, ideally I would like to have this one connected to the first model, but when i tried that it would not merge all features together, i think "merge" ran once for every file separately instead of once for all the files.? Integrate gives me a feature class with the averaged line from all the lines and "kernal Density" gives me a graphic representation of how much traffic there has been in different areas.   (in this example there is only two GPX-files, but the model is intended to run potentially for several hundred) attribute tables: ![enter image description here](http://i.stack.imgur.com/6oIU8.png) split_0 and split_1 is the result of the first model,and the SplitLine_Merge is the result of the second model. here you can see some of the problem, the attribute speed and seconds is the same as from split_1, even though the Shape_Lenght is different. I need it to average the speed from both tables when it integrates. so to sum up: 1) generalized lines from different tracks - Accomplished using Integrate 2) Number of tracks combined to make on line segment on the Integrate - not accomplished 3) average speed on the line segments - not accomplished in the generalized lines 4) average time(seconds) it takes to drive the line segment - not accomplished in the generalized lines 5) last date driven - not accomplished in the generalized lines 6) graphic representation of how many times a line segment is driven - accomplished using Kernal Density I know that this is a comprehensive question with several sub-questions so I'm happy with every answer/tip that can help me solve parts of my problem. Thanks in advance!


Problem on MapFishPrint in jpeg formate images

i am using MapFishPrint plugin for printing the my gis map .I am printing the map in two format PDF and JPEG format.my problem is that when i am printing the map in JPEG image formate then the map title and Comment is not printing with JPEG images But map is printing fine.the map title and comment is blank in JPEG image. But when i am printing the image PDF formate then all map is printing fine with map title and comment. So i want to know that JPEG formate is supported by JPEG image in MapFishPrint.? This my code on config.yaml file for A4 page size #=========================================================================== A4 portrait: #=========================================================================== # override the default outputFileName defined above # outputFilename: A4_Layout               mainPage:       rotation: true       pageSize: A4       header:         height: 50         items:           - !text             font: Helvetica             fontSize: 20             align: center             text: '${mapTitle}'            items:         - !map           spacingAfter: 30           width: 440           height: 580         - !text           text: '${comment}'           spacingAfter: 30            footer:         height: 30         items:           - !columns             items:               - !image                 maxWidth: 150                 maxHeight: 100                 align: left                 url: http://www.links-map.com/live/pymble/img/logo1.jpg               - !text                 align: right                 text: 'Page ${pageNum}'      Basically how to print Text with image Thanks


ggmap fortify start and end jumble

I have a census tiger line file for counties available from the Census website. I am trying to display this file on a `ggmap`, but `fortify` seems to jumble the regions. I can't find anyone else who has run into this error, but I'm certainly not doing anything out of the ordinary...               library(maptools)     library(ggmap)     library(rgeos)          # read shapefile     WGS84 <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0")     Counties <- readShapeSpatial("tl_2013_us_county.shp", proj4string = WGS84)          # limit to Ohio     Counties <- Counties[ Counties$STATEFP == 39, ]          # fortify     Counties.fort <- fortify(Counties, region = "GEOID")          # get base map and plot     OhioMap <- get_map("Ohio", source = 'google', maptype = 'terrain', zoom=7)     ggmap(OhioMap) + geom_polygon(aes(x=long, y=lat), fill = 'grey', alpha=0.4,                                    color = 'black', data = Counties.fort)      ![enter image description here](http://i.stack.imgur.com/CP0Mq.png)


cad to shapefile conversion

i have use arcgis10.1,arcobjects.net. i did following code                       ESRI.ArcGIS.Geoprocessor.Geoprocessor gp = new ESRI.ArcGIS.Geoprocessor.Geoprocessor();             gp.OverwriteOutput = true;                  ESRI.ArcGIS.ConversionTools.FeatureClassToShapefile tool = new ESRI.ArcGIS.ConversionTools.FeatureClassToShapefile();                       tool.Input_Features = @"F:\test buffer of JAMNAGAR SHP FILES ex3\FOOTPATH.dwg";             tool.Output_Folder = @"C:\Temp1";             gp.Execute(tool, null);      but after execute this code i got error "Error HRESULT E_FAIL has been returned from a call to a COM component. " if any one have any suggestion or any sample for this conversion or any link then send me please ..


finding length of a custom type array in mapbasic

Hello I am just getting started with MapBasic and I have run in to some trouble trying to count the items in a custom type i have tried accomplishing this by running a loop until an error occur. Problem is that the error trapping does not seem to happen i just keep getting a dialog with an out of range error message what am I missing here? My code is as follows:               Type KeyValuePair     Pos As Integer     Key As Integer     StringVal As String     End Type          Type PopupMenu_Items     ControlId As Integer     Values() As KeyValuePair     End Type          Sub CountPopupMenu_Items (countResult As Integer, arrayToCount() As PopupMenu_Items)     Dim i As Integer      Dim controlId As Integer     i = 1     Do While i = i         controlId = arrayToCount(i).ControlId         OnError Goto endOfArray         i = i + 1     Loop     last_exit:         exit sub     endOfArray:         countResult = i - 1         resume last_exit     End Sub


How to implement viewshed calculation?

I'm looking for a way to implement visibility (aka viewshed). I found a previous GIS SE question that suggested a few GIS applications (not what I need - I'm trying to embed the calculations), and also SAGA. SAGA looks more like what I'm looking for (C++, Java, C# classes), but I can't understand enough of the API documentation to figure out what classes I even need to instantiate (or subclass). I'm not tied to the idea of SAGA - I'm just looking for library / engine that can do intervisibility calcs. The constraints are that it needs to be:   * reasonably "light weight" (since I'd like to be able to cover embedded / mobile)   * reasonably open source (since my application will be open source, although I'm flexible on exactly which license as long as its (L)GPL compatible). Can anyone recommend such a library or engine, and provide a description of or link to which bit of the library or engine I need to use to viewshed calculations? Alternatively / in-addition, I'd also appreciate references to papers or tutorials that explain how to do these calculations in an efficient way (as applicable to embedded / mobile devices, so a GPU based implementation may not be as useful as something that is more general, and handles low power devices) My preference is an existing implementation, rather than creating something myself (which is a a backup option). _Edit: C++ isn't a firm requirement - anything reasonably portable (C#, C++, Java) will do. I'm trying for a library or templates rather than embedding a large application._