WPSExecute Format Issue

I have:               function createWPSExecuteRequest(namespace, layerName, popFilter) {         var request = OpenLayers.Format.XML.prototype.write(new OpenLayers.Format.WPSExecute().writeNode('wps:Execute', {             identifier: 'gs:Aggregate',             dataInputs: [{                 identifier: 'features',                 reference: {                     mimeType: 'text/xml',                     href: WFS_HOST,                     method: 'POST',                     body: {                         wfs: {                             featureType: layerName,                             version: '1.1.0',                             featureNS: namespace,                             srsName: 'EPSG:4326',                             defaultFilter: popFilter                         }                     }                 },                  input: {                     identifier: 'aggregationAttribute',                     data: {                         literalData: 'grid_code'                     }                 },                 input: {                     identifier: 'function',                     data: {                         literalData: 'Sum'                     }                 },                 input: {                     identifier: 'singlePass',                     data: {                         literalData: 'Yes'                     }                 }                  }],             responseForm: {                 rawDataOutput: {                     identifier: 'result'                 }             }         }));              return request;     }      but the generated GML is shown below:               <wps:Execute xmlns:wps="http://www.opengis.net/wps/1.0.0" version="1.0.0" service="WPS">         <ows:Identifier xmlns:ows="http://www.opengis.net/ows/1.1">gs:Aggregate</ows:Identifier>         <wps:DataInputs><wps:Input>             <ows:Identifier xmlns:ows="http://www.opengis.net/ows/1.1">features</ows:Identifier>             <wps:Reference mimeType="text/xml" xlink:href="/geoserver/wfs" xmlns:xlink="http://www.w3.org/1999/xlink" method="POST">                 <wps:Body>                     <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.0.0">                         <wfs:Query typeName="p_layer" srsName="EPSG:4326" xmlns:undefined="http://pe.example.com/nig"/>                     </wfs:GetFeature>                 </wps:Body>             </wps:Reference>             </wps:Input>         </wps:DataInputs>         <wps:ResponseForm>         <wps:RawDataOutput>         <ows:Identifier xmlns:ows="http://www.opengis.net/ows/1.1">result</ows:Identifier>         </wps:RawDataOutput>         </wps:ResponseForm>     </wps:Execute>      Any idea what I may be doing wrong? **UPDATE:** I figured part of this out. Its the braces. A revised version that generates everything except the filter is shown below:               function createWPSExecuteRequest(namespace, layerName, popFilter) {         var request = OpenLayers.Format.XML.prototype.write(new OpenLayers.Format.WPSExecute().writeNode('wps:Execute', {             identifier: 'gs:Aggregate',             dataInputs: [{                 identifier: 'features',                 reference: {                     mimeType: 'text/xml',                     href: WFS_HOST,                     method: 'POST',                     body: {                         wfs: {                             featureType: layerName,                             version: '1.0.0',                             featureNS: namespace,                             srsName: 'EPSG:4326',                             defaultFilter: popFilter                         }                     }                 }},                 { identifier: 'aggregationAttribute',                 data: {                     literalData: {                         value: 'grid_code'                     }                  }},                 {identifier: 'function',                 data: {                     literalData: {                         value: 'Sum'                     }                 }},                 {identifier: 'singlePass',                 data: {                     literalData: {                         value: 'Yes'                     }                 }}             ],             responseForm: {                 rawDataOutput: {                     identifier: 'result'                 }             }         }));              return request;     }      The filter is not generated yet. Any idea will be appreciated.


Using PostgreSQL with QGIS and ArcGIS?

Is it possible to use PostgreSQL with QGIS and ArcGIS? i.e. one database for different clients. Is it anything that I have worry or any problems to use the DBMS from two different software? I know that I need PostGIS and ArcSDE and I can't manipulate or store raster data with PostGIS 1.5 in PostgreSQL from QGIS. What would you recommend or what to avoid?


Convert spatial Geographic data from sql server 2008 to Json

Is there any example to convert spatial data (point,line,polygon) stored as geographic in sql server 2008 into Json and display data using Arcgis Javascript Api 3.0.I have tried converting geographic data into Json given Geometry to Json but as my data stored in geographic data format it failed.Thank you in advance.


Analysis of lines in a point in QGIS

I have 33 layers, 1 with points and 32 with lines. I want to extract/highlight the points where more than **one** line is intersecting. I have tried with buffers but I'm not getting anywhere.


Underwater terrain modeling

My task is to model the surface of a lake bed, calculate volume of it and visualize it as 3D a surface. I have acquired depth measurements in a small lake using echosounder. The precision should be pretty high (<5 cm for depths and 3-5 cm for RTK GPS position). Points are placed regularly with 2.5 m intervals between them. I also have digitized shoreline from high resolution orthomosaic. Which would be the best interpolation algorithm for such tasks? (QGIS, GRASS, SAGA, maybe some other free software does such things?). The biggest challenge has been to find an algorithm, in which I can define the shoreline as 0 value which has to be taken in account.


Create a Leaflet custom checkbox control

I want to create a custom checkbox control which will simply set a flag in jquery/javascript: if checked the flag = 'clustered' or if unchecked flag = 'unclustered'. So far I have a control on the map but its not a checkbox and how do i get the state of the checkbox - if its checked/unchecked. code:               function MapShowCommand() {       alert("checked / unchecked"); //set flag     }          L.Control.Command = L.Control.extend({         options: {             position: 'topleft',         },              onAdd: function (map) {             var controlDiv = L.DomUtil.create('div', 'leaflet-control-command');             L.DomEvent                 .addListener(controlDiv, 'click', L.DomEvent.stopPropagation)                 .addListener(controlDiv, 'click', L.DomEvent.preventDefault)             .addListener(controlDiv, 'click', function () { MapShowCommand(); });                  var controlUI = L.DomUtil.create('div', 'leaflet-control-command-interior', controlDiv);             controlUI.title = 'Map Commands';             return controlDiv;                 }         });     var test = new L.Control.Command();     map.addControl(test);


Caching layers on the fly and seeding combination

I have created a map in geoserver Im am using geoserver 2.5.1 with the intergrated geowebcache. The map seems to be working fine and the seeding also seems to be working fine provided I don't use a bounding box for some strange reason. But also want it to cache on the fly as well so I can seed the main sections of the map that would be used and have on the fly caching in the areas that arn't being used as much. Am I mistaken in thinking that it saves those certain tiles like the seeding does or does it do it virtually or save in a temp tile folder and then goes away instantly tht person either closes the map program. I have set everything that I could set but for some reason on the fly caching does not seem to work. Am i missing something in the settings? Or am misunder standing something with regards to it as in the geoserver documention states that it can combine both on the fly caching as well as seeding the layers.


Update row with another not None row

I'm working with an attribute table in ArcMap 10.2 where there are some 30 fields mainly consisting of NULL values. I don't want to display to users a table full of NULL values as it'll be difficult to read. See below image: ![enter image description here](http://i.stack.imgur.com/x3lvQ.png) So the result I'd like to create is a Field called **SUMMARY** that includes the " _field name_ : _Value_. _field name_ : _Value_ " of any row that is NOT NULL **e.g. A: T. D: 10/01/2011** In this way a user can select a feature with the identify tool and see what values there are rather than scrolling through loads of NULL values. Does anyone have suggestion on the general construct of the script? Including the da.SearchCursor, da.UpdateCursor, and I guess by way of row index if statements to query if not none add to SUMMARY field with delimators ":" after field name and "." after row value.


How to programmatically add an algorithm from external plugin to QGIS processing toolbox?

I'm writing a python plugin for QGIS and I'd like to know if it's possible to provide an interface to it from the processing toolbox. I can see how the existing algorithms are loaded into processing via providers, but I don't see any obvious way of loading an algorithm without modifying the module. Perhaps it isn't possible and I need to submit a feature request, or attempt at coding it myself. I was expecting something like:               import processing     from myplugin import mysweetalgorithm     processing.load_algorithm(mysweetalgorithm)


Software for opening .img files

I have some cpc precipitation data stored in IMG format.Could anyone suggest a software beside ERDAS imagine or Arcgis which can be freely downloaded?


D3.js with dynamic geojson data (SQL server/.NET)

I have a $.ajax function which returns the geojson data from the asp.net mvc server method. I want to know how I could send that as a parameter to d3.json() function. Better if anyone shared code snippet to proceed. I have already tried as below but no success. d3.json(returnedGeojsonObject, function(error, json) {} Exact code snippet $(document).ready(function () { $('#btnGetMap').click(function () { getShopsFrom(1, 1); }); });               function getShopsFrom(mallID, floorID) {         var params = { MallID: mallID, FloorID: floorID };              $.ajax(         {             url: '<%:Url.Action("GetShops") %>',             dataType: 'html',             data: params,             type: "POST",             success: function (result) {                 alert("Got the details..voila");                 parseResultShopDetails(result);             },             error: function (req, status, error) {                 alert("Unable to get shop details : " + error);             }         });     }          function parseResultShopDetails(data) {         var width = 960, height = 1160;              var projection = d3.geo.albers()                         .center([0, 55.4])                         .rotate([4.4, 0])                         .parallels([50, 60])                         .scale(1200 * 5)                         .translate([width / 2, height / 2]);              var path = d3.geo.path()                     .projection(projection);              var svg = d3.select("body").append("svg")                         .attr("width", width)                         .attr("height", height);                 d3.json(data, function (error, uk) {             svg.append("path")                     .datum(topojson.feature(uk, uk.objects.subunits))                     .attr("d", path);                     });             }


Where is the Merge selected features button in QGIS 2.2?

Ever since I updated to QGIS 2.2 I've been a bit lost with the new menu layout (I was running on Wroclaw until then). I can't seem to find the "Merge selected features" button. I need to merge several polygons inside a shapefile and can't find how to do it in this version. Is the Merge selected features button still included on QGIS 2.2? Where can I find it? What's the current icon for it?


How to generate TFW files for a list of images?

I have a few digital images that are not georeferenced. E.g.:               Img01.tif       Img02.tif       ...     Img99.tif        On the other hand, I have a text file, which contains georeferencing information for these images. E.g., I have a file `c:\tfws\list.txt` which contains the following information:               Img01.tfw 0.25 0 0 - 0.25 456258,125 4569852,125       Img02.tfw 0.25 0 0 - 0.25 456586,125 4570001,125        Img03.tfw 0.25 0 0 - 0.25 456952,125 4570300,125       ...     Img99.tfw 0.25 0 0 - 0.25 458412,125 4575123,125        I need to generate a TFW file for each image in the text file. E.g. create a file named `C:\tfws\img??.tfw` containing the following information:               0.25       0       0       -0.25       456258,125       4569852,125        Does anyone know of a tool able to generate these TFW files automatically?


What are open data sources about current sea state?

I wonder where (links/APIs) can be obtained data that is related to current (alike realtime - no past) sea state data (worldwide, but regional is good two, salt\waves length etc)? I am only interested in dynamically changing data as close to realtime as possible.


How to increase resolution of 3D image in GRASS NVIZ

I am trying to make 3D visualization of terrain data and orthophoto using GRASS NVIZ via QGIS. I have been able to import both the DTM and the orthophoto in NVIZ window and successfully draped the image over DTM. But the resolution looks very poor. Is there any way to increase the resolution ? It seems like NVIZ GUI is not so user-friendly. For example, in Arcscene, there is a option to optimize the rendering that can increase the resolution of the 3D draped image. Here is the screeshot from GRASS NVIZ:![enter image description here](http://i.stack.imgur.com/QFU07.jpg)


Creating batch file from ERDAS IMAGINE Spatial Modeler?

How to create batch file from the ERDAS IMAGINE Spatial Modeler? I want to do the same process to different input images, something similar which we can do by batch processing of the existing programs in Erdas. If there is no inbuilt function then one can use model maker to make *.gmd (graphical model) file but how to batch it?


SQL query to have a complete geojson feature from PostGIS

I would like to get a geojson feature with properties from PostGIS. I have found an example to have a feature collection but I can't make it works for just a feature. so far I tryed to modify the feature collection query of the example. but the output is not valid.


QGIS Intersection problem

The QGIS Intersection tool is not working for two shapefiles in Version 2.2 or 2.4 (haven't tried earlier versions). I get a new shapefile with all the attributes, but no features. However, when I use ArcMap with the two shapefiles, the intersection happens correctly. The shapefiles are both in the same coordinate system. I have used the QGIS tool with success with other shapefiles that I got from a QGIS tutorial. The shapefiles that I'm using are all from the Data BC website. I've uploaded two samples to http://vulnerablevancouver.ca/qgis_intersection/ if anyone wants to have a look. I put the zip files there that include the metadata.


Making 'persistent' ArcMap tool without creating Python add-in?

I want to make an addon that, while editing in ArcMap, will track where I pan and move the viewport. I have the pieces together, but I don't like debugging with the Python Add-In wizard, because I have to recompile and restart ArcMap each time I modify the code. Is there another way to create this type of 'persistent' script without using an add-in? If I make a toolbox tool and run it, will it run in the background, and how would I stop it? All the toolbox tools I've ever created have a certain task they perform and, once completed, the tool stops running. Same with ModelBuilder.


QGIS WEB Cliente para windows 7

Bom dia, estou com alguma dúvidas a respeito do QGIS WEB Cliente que não consigo achar uma solução. Eu sei que tenho que fazer a instalação do Qgis Serve como esta escrito neste tutorial Anita Graser QGIS Server on Windows 7 I also know that I have low ZIP QGIS Web Client, but where should I put it in Windows? After that should follow this documentation In my attempt put the contents of the ZIP file into C: \ OSGeo4W \ apache \ htdocs \ The map model HelloWorld works But I have no idea where to go, my project is stopped.


How to create a rainfall raster from rainfall contours?

I downloaded a free rainfall data from the Philippine GIS Data Clearinghouse (PhilGIS) which are vector lines. How can I extract the data points from this? I want to extract it because i will be using these points to make a raster file for the rainfall data that i downloaded. I'm using QGIS and global mapper for this.


Is there an equivalent to arcpy.ImportToolbox() in 9.3.1?

Is there an equivalent to `arcpy.ImportToolbox()` using `ArcGIS 9.3.1`? I'm assuming it would use `arcgisscripting` instead of `arcpy` but I can't find any documentation.


QGIS 2.2 crashes with custom SVG fill

This may be a bug, but I am not sure how or if I should submit it. I wasn't able to find much with a Google search. I developed a custom SVG fill and implemented it. When I go to change the style settings again QGIS crashes. I had this happen on two different machines, both using Windows 7 and QGIS 2.2. You can see the custom fill and the error in the following screenshot:![QGIS Crash using custom SVG fill.](http://i.stack.imgur.com/EkiU5.png) I also have the .dmp file, not sure how to view these, didn't want to open with Notepad and I don't see a way to attach the file here. **Has anyone experienced this type of crash?** **EDIT:** I have gotten this a number of times doing various task in different map projects.


Can I set individual scale levels for an ArcGISDynamicMapServiceLayer?

I have a map service that contains 6 layers. I turned one layer into a feature layer. I used this layer in a query and show the result. I then would like to just show the other 5 layers as a ArcGISDynamicMapServiceLayer, but I didn't set any min/max scale levels for the layers when I published the MXD. I figured it'd be best to do this programmatically, but I can't figure out how. Can I set the min/max scale levels independently for each layer in the ArcGISDynamicMapServiceLayer? If not, what's the best way to approach this? Thanks! EDIT Can someone explain to me how this question isn't clear? I don't believe it could get any more clear.


How to remove ArcGISTiledMapServiceLayer completely

I have a problem related to "ArcGISTiledMapServiceLayer". I have an ArcGISTiledMapServiceLayer map service, which I have added to Map object. This service work as BaseMap. Now I want to add another ArcGISTiledMapServiceLayer map service layer to the same Map object. Before adding the new layer, I just removed the previous tiled map service layer by using method `objMap.removeAllLayers()`. But this method only removes the layer not the extent and it create problem when I add a new layer. So I want to remove completely the ArcGISTiledMapServiceLayer (layer, extent, scale etc.). I have tried all the below methods   1. `objMap.removeAllLayers();`   2. `objMap.removeLayer(objMap.getLayer('myLayerID'));` Any solution?


Correct WFS URL Format (Gaia WFS-T Extender)

I am trying to add a GeoServer WFS layer to Gaia WFS-T Extender for editing. The available layers are available alright but I am unable to load them into Gaia. Using the same URL I was able to add and edit the layer in both Quantum GIS and UDig. Any hints on where I could be missing it with Gaia? [ why Gaia: It is light-weight and packs functionality I'm looking for]


Draw polygon with gradient color in Geoserver

I have polygon like sector. I want to fill it with gradient color. Is it possible to do this in Geoserver? Thanks. **UPDATE** I found some code with filling image               <FeatureTypeStyle>      <Rule>        <PolygonSymbolizer>          <Fill>            <GraphicFill>              <Graphic>                <ExternalGraphic>                  <OnlineResource                    xlink:type="simple"                    xlink:href="colorblocks.png" />                  <Format>image/png</Format>                </ExternalGraphic>              <Size>93</Size>              </Graphic>            </GraphicFill>          </Fill>        </PolygonSymbolizer>      </Rule>      i can put here big circle which is gradient, but does it has rotation or position of image?


CQL Filter in Geoserver very slow?

I am experimenting with CQL filtering in Geoserver. I have two layers, one (A) has 8000 features, the other one (B) has 3000 features or so. I would like to search of layer A for intersection with layer B's features. Not all of Layer B's features, just a few based on their feature ID. I've been searching around and I can't find a way yet to do the **feature ID + geometry filter** , I hope someone can suggest me how to do that. So in the meantime, I have tried with another attribute I've put in, though, and the CQL Filter is currently like this with WMS:               layers=A     CQL_FILTER=INTERSECTS(the_geom,querySingle('B','the_geom','OBJECTID=6139'))      When I run this in my browser , the query goes forever and timed out. There is an error in the jetty : **org.geoserver.platform.ServiceException: This requested used more time than allowed and has been forcefully stopped. Max rendering time is 60.0s** The example from Geoserver (the one using restricted layer vs bug sites) doesn't have this problem. I'm wondering if I'm writing the query wrong, whether my PC's just not up to it, or something else. Should I be doing this in WFS instead?


Uploading unprojected shape files using GeoServer REST API

I have a number of shapefiles that lack the `.prj` file. I'm able to upload them to geoserver using the rest api. However, I'm unable to serve them after upload. Which brings me to the question as to whether I need to reproject before uploading to geoserver. If I need to reproject, then which options are available? Thanks all in advance. I'm using GeoServer 2.1.3. UPDATE: Layers created from the shapefile lack projection information. When I update the `Declared SRS` field, the layers are served without issues.


Geocode adress with appartments

I have a table with addresses which I try to geocode in ArcGIS 10.1 with a single field locator. The problem is that some of the addresses includes an apartment number which I do not have in my locators address list. Logically that would be easy to solve when comparing the two address-fields, one could just compare the table addresses with the locators where the address 'BEGINS WITH' but ArcGIS seems to match the two rather with 'IS' thus giving all the addresses with appartment numbers a score of 0. Is it possible to get around this in any other way than removing all apartment numbers from the table I want to code?


Error when opening attribute table of a feature layer with a programmatically generated join

I'm creating a query table with IQueryDef and then using MemoryRelationshipClassFactory to create an inner join between a feature layer and the query table. Both the query table and the join are created without incident and work correctly in ArcMap, (only the features with matching records in the query table are displayed on the map). However, when I attempt to open the attribute table for the feature layer in ArcMap I get the following error: _Could not load data from the data source. If you can correct the problem, press the refresh button to reload data. Possible problems can include bad network connection, invalid field, etc. Field is not editable._ The table then opens with empty rows. I'm not sure what the "Field is not editable" message is telling me. I thought perhaps there was something wrong with the query table, but I can create the table and add it to ArcMap as a stand alone table and it opens correctly. It seems there is something in the way I am doing the join that is causing trouble with the feature layers Display Table. Below is the relevant code:                   Dim qd As IQueryDef = ftrWksp.CreateQueryDef         qd.Tables = "swCBCleaning"         qd.SubFields = "*"         qd.WhereClause = pre & "CompletedDate" & suf & " >= date '2012-11-08 00:00:00'"              Dim qn2 As IQueryName2 = CType(New TableQueryName, IQueryName2)         qn2.QueryDef = qd         qn2.PrimaryKey = "swCBCleaning.OBJECTID"         qn2.CopyLocally = False         Dim dsName As IDatasetName = CType(qn2, IDatasetName)         Dim dataSet As IDataset = CType(ds.Workspace, IDataset)         dsName.WorkspaceName = CType(dataSet.FullName, IWorkspaceName)         dsName.Name = "swCBCleaning"         Dim name As IName = CType(qn2, IName)         Dim tblJoin As ITable = CType(name.Open, ITable)              Dim memRelClassFactoryType As Type = Type.GetTypeFromProgID("esriGeodatabase.MemoryRelationshipClassFactory")         Dim memRelClassFactory As IMemoryRelationshipClassFactory = CType(Activator.CreateInstance(memRelClassFactoryType), IMemoryRelationshipClassFactory)              Dim relClass As IRelationshipClass         Dim dispRelClass As IDisplayRelationshipClass              relClass = memRelClassFactory.Open("Join", geoFtrLayer.DisplayFeatureClass, "GlobalID", CType(tblJoin, IObjectClass), "FeatureGlobalID", "forward", "backward", esriRelCardinality.esriRelCardinalityOneToMany)         dispRelClass = CType(geoFtrLayer, IDisplayRelationshipClass)         dispRelClass.DisplayRelationshipClass(relClass, esriJoinType.esriLeftInnerJoin)


How to enable tilestache-seed logging?

We're migrating from a working 'tilecache' server to a new 'tilestache' server. I have been chasing my tail for 2 days now trying to figure out how to get tilestache-seed to cache some tiles. I'm failing miserably and neither the tilestache docs nor google are helping at all. It must be me. I have a configuration that looks like this:               {         "cache": {             "name": "Disk",             "path": "/opt/tiles",             "umask": "0000",             "dirs": "safe"         },         "logging": "debug",         "layers": {             "staging": {                 "provider": {                     "name": "url template",                     "template": "http://example.com/common/services/imagegen/?mapservice=cities&cs=102113&crs=102113&layers=6&format=png&bbox=$xmin,$ymin,$xmax,$ymax"                 }             }         }     }      A typical tilestache-seed command line looks like this               ./seed -c tilestache.json -l staging -b -65.167510000 47.589390000 -65.167490000 47.589410000      Note: `seed` is a symlink to `tilecache-seed.py` \- PIA to type the full path so much :) Note also that 'imagegen' is a proprietary application that sits in front of ArcIMS to manage other things about the image generation. Nothing magic, just additional capability that is irrelevant in this particular case. Running the command line and _nothing_ at all happens. No output at all from tilestache-seed. It is not making a call to the URL in the config file - I can see that because imagegen writes the generated image to an output directory on the web server, and there's nothing being written in response to seed being run. If I manually massage the seed command line and substituting the -b parameters into the proper spots in the URL and paste that in a browser, I generate an image just fine. So why isn't seed doing what I think it should be doing? I figure there must be a way to see what's going on inside seed. Well, no command line option to turn on "verbose" mode. Ah, I can use the Test cache and turn on verbose there. No, that changes nothing - no output of any sort. Note that the example in the tilestache docs for the Test cache appears to be in error:               "verbose": True      throws an error when running seed.               "verbose": true      fixes the error, but doesn't generate any additional output. Can someone provide me guidance, please, on a) what might be wrong with my configuration and/or b) how to turn on some sort of diagnostic mode so that I can debug this further? I'm not a python guy, though that's a skill I probably ought to acquire. Thanks! Bo PS: Any reference to 'tilecache' should read 'tilestache' - been typing the former for too many years :)


To Use or Not to Use Relationship Classes?

I am building a government database from a massive amount of data. Now I am following the standard approach of the ESRI Local Government Data Model, but modified to fits the needs and wants of a specific government. I've noticed that in the Local Government Data Model ESRI uses a lot of relationship classes. I've used relationship classes in the past, but tend to steer clear of them. This is because I usually don't see a reason why I would want attributes about a shape in a separate table from the feature class. For almost all of the data I am assembling I still cannot see a reason for this. For Parcel Fabric I can see some reasons where it would be helpful, but still I believe it would be best to still keep the attributes directly in the feature class. Now to be clear I am going to use relationship classes to reference feature classes to feature classes, this question is more specified to why I would store data in a non-spatial table that could just as easily be store in a feature class. Anyways thanks for any input folks!


Is it possible to do a conversion and add layers with command prompt in ArcGIS 9.3?

Is it possible to operate ArcGIS 9.3 from command prompt? I want to convert an ascii file to raster format and adding layers to the obtained raster file. How can this be done ? Thanks .


onSelect sometimes firing before onUnselect

I have a map that displays satellite coverage zones. Some of these coverages overlap. I am using the SelectFeature functionality to fire off **onSelect** and **onUnselect** events. However, when moving from one feature to another (without moving on the map itself first), the onSelect method sometimes fires before the onUnselect method. As the onUnselect method is being used to clear a popup created as part of the onSelect method, this means that sometimes the popup is created and then destroyed, rather than destroyed and then created. Any ideas? Below is a snippet of the relevant portion of code (in CoffeeScript):               beams = JSON.parse localStorage.getItem("beams")     $.each beams, (i, item) ->         beamTrendData = getTrendData4Beam item.ShortName         beamStatus = getBeamStatus beamTrendData         switch beamStatus             when "green"                 style = styleMapGreen             when "orange"                 style = styleMapOrange             when "red"                 style = styleMapRed         #if (item.ShortName is "G17") or (item.ShortName is "T14R") or (item.ShortName is "A5") or (item.ShortName is "YAM") or (item.ShortName is "IS14") or (item.ShortName is "W7")         #   visible = true         #else         visible = true         coverageLayer = new OpenLayers.Layer.Vector(item.ShortName,             strategies: [ new OpenLayers.Strategy.Fixed() ]             visibility: visible             protocol: new OpenLayers.Protocol.HTTP(                 url: "resources/satelliteFootprints/" + item.OuterBeamEIRPKML                 format: new OpenLayers.Format.KML(                     extractStyles: false                     extractAttributes: false                     maxDepth: 2                 )             )             styleMap: style         )         Application.map.addLayer coverageLayer         Application.hoverControl = new OpenLayers.Control.SelectFeature coverageLayer,             multiple: false             hover: true             toggle: true             clickout: false             onSelect: beamSelect             onUnselect: beamUnselect             callbacks:                 click: (feature) -> beamClick feature         Application.map.addControl Application.hoverControl         Application.hoverControl.activate()      I'm basically adding several KML layers, and for each one associating a SelectFeature object that in theory triggers events for both hover over and click (as I want to add different functionality to the hover and click events). In IE9 (other than unselect and select sometimes firing in the wrong order) it all seems to work just fine. Strangely though for any other browser the events only fire on the topmost layer. EDIT: I fixed the cross browser issue, by taking the adding of the hoverControl out of the loop and simply adding it an array containing the vector layers all at once. So:               beams = JSON.parse localStorage.getItem("beams")     beamFeatures = [];    # <------------     $.each beams, (i, item) ->         beamTrendData = getTrendData4Beam item.ShortName         beamStatus = getBeamStatus beamTrendData         switch beamStatus             when "green"                 style = styleMapGreen             when "orange"                 style = styleMapOrange             when "red"                 style = styleMapRed         #if (item.ShortName is "G17") or (item.ShortName is "T14R") or (item.ShortName is "A5") or (item.ShortName is "YAM") or (item.ShortName is "IS14") or (item.ShortName is "W7")         #   visible = true         #else         visible = true         coverageLayer = new OpenLayers.Layer.Vector(item.ShortName,             strategies: [ new OpenLayers.Strategy.Fixed() ]             visibility: visible             protocol: new OpenLayers.Protocol.HTTP(                 url: "resources/satelliteFootprints/" + item.OuterBeamEIRPKML                 format: new OpenLayers.Format.KML(                     extractStyles: false                     extractAttributes: false                     maxDepth: 2                 )             )             styleMap: style         )         Application.map.addLayer coverageLayer         beamFeatures.push coverageLayer    # <------------          Application.hoverControl = new OpenLayers.Control.SelectFeature beamFeatures,         multiple: false         hover: true         toggle: true         clickout: false         onSelect: beamSelect         onUnselect: beamUnselect         callbacks:             click: (feature) -> beamClick feature     Application.map.addControl Application.hoverControl     Application.hoverControl.activate()      I haven't checked yet whether the select and unselect events fire in the right order all the time yet. Thanks, Kurt


Extract Values to Point from Raster Issue

I am trying to use the tool Extract Values to Points on a File Geodatabase raster, but I am getting some weird errors. I've used this tool multiple times before with no issue, but not since I upgraded to 10.1, so I'm wondering if there is an issue with it. The errors I got back so far have been:   1. This error happened when I specified a shapefile as the output: _The specified geodataset is invalid for the conversion. The syntax help for the tool you are using should indicate the supported types of input data. For example, for the particular conversion, a rasterlayer may be an invalid type of geodataset._   2. Next I have an error about it needed at least one band in the raster, the raster does has one band. This occurred when I specified a File Geodatabase as the output.   3. Finally I had the good ole Error: 999999 for no apparent reason. The raster is a DEM derived from LiDAR that was not created by a different company. It's too large to be made available for testing, but here are its specs:               Raster_Information          ColumnsandRows  5119,8241     NumberofBands   1     CellSizeXY      1, 1     UncompressedSize    160.93 MB     Format      FGDBR     Source_Type     Generic     Pixel_Type  floating point     Pixel_Depth     32 Bit     Colormap    absent


Colour matching satellite imagery from different time periods

Is it possible to colour match two separate satellite images taken at two different times? I am working with two high resolution, multispectral images that were taken approximately one year apart and I want to do some change analysis with them. Unfortunately one of the images is extremely washed out in urban areas and I'm trying to correct it (see attached samples). **Image From Time 1** ![Sample Time 1](http://i.stack.imgur.com/cENsF.jpg) **Image From Time 2** ![Sample Time 2](http://i.stack.imgur.com/RuCq4.jpg) Ideally, I would like to have both of the images to appear as similar as possible, preferably like the first one, so I can identify changes between the two time periods. I am using ArcGIS 10.2 Professional, but if there are other options to deal with this I would like to hear about those as well.


Using reader Feature Type Name as Writer Feature Type Name in FME

I've search and searched for a way to do this, but no method seems to work. Basically I want to use the input feature class name (feature type in FME) as the output name. This normally works, however the tool I'm using (ChangeDetector) adopts the "revised" input as the output Feature name instead of the "original" input name. Even linking the "original" input name to the output name attribute does not work. It keeps the "revised" input name. I've tried using various transformers like the SubstringExtractor to pull the name and set the writer to fan out by this attribute, but even that fails. While I would prefer to get this working within the workbench, I do plan on using this workbench in a python script, so it would be handy to get it working through that as well. I setup the workbench to run using               subprocess.Popen(["pathToFme", "pathToWorkbench"])      This works just fine. Subsequently I tried adding parameters at the end, to input the name based on code samples I found while searching this...               subprocess.Popen(["pathToFme", "pathToWorkbench"] --fileName: "nameOfOutput")      with "fileName" being the name of the parameter within the workbench, but I can't seem to get it to run, let along pick up the parameter input from the script. I couldn't get os.system to properly execute the workbench, so I want to stay away from that. Basically what I am trying to achieve is the ability to take a large number of inputs, run them through the changeDectector, and output with the same name as the "original" input. Any help or suggestions would be appreciated. Thanks.


What are resources to learn ArcGIS for someone who already knew Python

I know that there are several resources recommendation questions, however, they are all from GIS people trying to learn Python (for example, here). My situation is the opposite: I have worked extensively in Python, but know nothing about ArcGIS. What is the best way to learn? To be more specific, I want to learn the following:   * Create maps from shapefiles   * Incorporate other social science-y data on map. For example, conflict site, social network, economic activity, etc.   * Most importantly, statistical modeling with a spatial component (network, hotspot analysis, etc.)


Graphically depict statistical difference between thematic rasters

The graph shows a regression between two raster datasets with classified tree canopy cover: class 1 = 0-10%, class 2 = 10 - 20% ...class 6 = >60% (Figure 1, Figure 2). The spatial resolution for raster 1 is 1m while the spatial resolution for raster 2 is 30m. Note that raster 2 was resampled to 1m for analysis purposes in the script (Appendix A). I am looking for an elegant, graphical approach to statistically describe the difference between two thematic rasters at different spatial scales. My first attempt was to 1) create a raster stack, 2) randomly sample the stack, 3) run a regression and descriptive statistics, 4) graph the results (Appendix A). What other graphical methods exist for statistically comparing thematic raster data? I usually use R, Python and MATLAB, although I would be happy with generic solutions too. * * * **Figure 1** ![enter image description here](http://i.stack.imgur.com/m7Ceq.jpg) * * * **Figure 2** ![enter image description here](http://i.stack.imgur.com/ITbCt.png) * * * **Appendix A**               require(raster)     require(spatstat)     require(ggplot2)     require(gridExtra)     require(hexbin)     require(rgdal)          # Read the TIF data     file = 'C:/path/raster1.tif'       file2 = 'C:/path/raster2.tif'                    # Create raster objects from the single band canopy products     r = raster(file)     r2 = raster(file2)          # Resample r2 to 1m for analysis purposes     r2 = resample(r2, r, method = "ngb")          # Stack the two raster layers for analysis     raster <- stack(r, r2)          # sample random locations in raster stack and report values in a dataframe     df = data.frame(sampleRandom(raster, size=1000, cells=TRUE, sp=TRUE))          ## Do the regression and plot the results     new = df$X3711203_ne     old = df$X3711203_ne_30m          # Calculate RMSE and other values     fit <- lm(old ~ new)     rmse <- round(sqrt(mean(resid(fit)^2)), 2)     coefs <- coef(fit)     b0 <- round(coefs[1], 2)     b1 <- round(coefs[2],2)     r <- round(sqrt(summary(fit)$r.squared), 2)          # Build equation see ?plotmath     eqn <- bquote(italic(y) == .(b0) + .(b1)*                     italic(x) * "," ~~ r == .(r) * "," ~~ RMSE == .(rmse))          eqn_text = as.character(as.expression(eqn))          # Get max range to automatically update plot size     li = c(range(new)[2], range(old)[2])     range = tail(sort(li),1) + 0.05          # Plot results     p1 = ggplot(df, aes(new, old))  +        geom_point() +       geom_smooth(method = lm) +       stat_sum(aes(size = ..n..)) +       xlab("Cover Class (1m)") +       ylab("Cover Class (30m)") +       annotate("text", x= 0, y= range, label=eqn_text, hjust=0, size=4, face="italic", parse=TRUE) +       xlim(0, range) +        ylim(0, range) +       ggtitle("Tree Canopy Cover Class")          p1


How to export query result on raster table in PostGIS

I did success query on a raster table in postGIS(raster data clipped by a vector polygon), now I want to view the query result in QGIS. But how to export the query result to another table. In vector data, it's simple, just "CREATE TABLE schema.someTable AS (...sql in here...); How about raster data? Any suggestion? Thanks.


I want to get rid of unnecessary map and data from shp file, how?

I have loaded the base map from esri's defualt provided maps but its for the whole world, and i have gotten the roads data too, but its for whole continent. I only need for one city i am going to be working on to save my RAM. I am using arcgis 10.2


Unable to load PostGIS vector layer into QGIS using python cookbook recipe

I'm having trouble loading a PostGIS layer using Dufour (from osgeo4w), using the guide at http://www.qgis.org/en/docs/pyqgis_developer_cookbook/loadlayer.html. I just get the 'invalid layer' message. My parameters appear to be perfect.               from PyQt4 import QtGui     from PyQt4 import QtCore     from qgis.core import *     from qgis.gui import *     from ui_townplanning import Ui_TownPlanning     # create the dialog for zoom to point               class TownPlanningDialog(QtGui.QDialog):         def __init__(self, iface):             QtGui.QDialog.__init__(self)             # Set up the user interface from Designer.             self.ui = Ui_TownPlanning()             self.ui.setupUi(self)             self.iface = iface             self.canvas = self.iface.mapCanvas()             self.addPostGISLayer('myhost', '5432', 'glis', 'user', 'pass',                                  'tpadmin', 'planning_application_boundary_validation',                                  'wkb_geometry')              def addPostGISLayer(self, host, port, dbname, username, password,                             schema, table, geom_col):             '''Adds a PostGIS table to the map'''             uri = QgsDataSourceURI()             uri.setConnection(str(host), str(port), str(dbname), str(username), str(password))             uri.setDataSource(str(schema), str(table), str(geom_col), '', str(key_col))             vlayer = QgsVectorLayer(uri.uri(), str('test'), 'postgres')             # QgsMapLayerRegistry.instance().addMapLayer(vlayer)             if not vlayer.isValid():                 msgbox('Layer not loaded', uri.host(),                        uri.database(), uri.port(), uri.username(), uri.password(),                        uri.schema(), uri.table(), uri.geometryColumn())          def msgbox(*args):         a = [str(i) for i in args]         text = '\n'.join(a)         qmsg = QtGui.QMessageBox()         qmsg.setText(str(text))         qmsg.exec_()      **Code now works after correction**


How to return geo-tagged social media posts?

I have seen quite a few ESRI web apps that utilize Geo-Tagged social media posts as a point data layer. I'm wondering how one goes about querying things like Twitter, Flickr and YouTube to return posts related to a specific location. I'm also hoping to be able to filter these results based on hash- tags. Twitter is my main interest, but I know when I get that working I'll have people clamoring for the others. My hope is to add this to our public web viewer where citizens view pictures, videos and tweets about our public lands. Also, are these Geo-Tagged posts limited by your own personal followers / friends? Here is an ESRI example of what I am looking for: ESRI Australia Wildfires Map


Can FME Desktop (node-locked) be installed on an Amazon EC2 or VPC instance?

I did do some searching in the FME help and the legal notices, but I'm not getting any hits (probably due to vocabulary). Are there any legal restrictions for FME Desktop that would prevent me from installing onto an EC2 or VPC Amazon instance? I'm not sure if (or how much) that would compete against the FME Cloud offering. I'm mostly interested for test and evaluation, but I am also curious if anyone is using an EC2/VPC as a permanent home for FME Desktop. This another question that is inspired by this closely related question:   Is there an optimum RAM configuration to use for FME Desktop?


ArcMap 10.0 won't load

ArcMap 10.0 on a Windows 7, 32-bit computer stopped loading two days ago. It worked perfectly before then. ArcCatalog, etc., all continue to open fine. Task Manager indicates ArcMap is trying to load, but very, very slowly. I deleted Normal.mxt, and tried to restart ArcMap, but it has not yet successfully opened up and has not re-created Normal.mxt.


Error opening feature class from gdb

I have a file geodatabase that a person gave to me for work. I try to open the feature classes in Arcmap 10.2 and some of it opens perfectly but with the others appears a message that says: Error opening feature class. Sometimes when I try to open the properties of the layers I get "Failed to edit the selected object(s)." What can I do in this case? ![enter image description here](http://i.stack.imgur.com/uOdzG.png)


How to start a local GIS Users Group

My immediate geographic area is lacking an organized local GIS Users Group. I belong to a group that meets a couple times a year, but its about 1.5 hour drive. I'm thinking about starting up a Users Group for our County, we have about 50 users just in my organization that I think could benefit, not to mention the other users from different organizations. Problem is, I've never done anything like this and don't exactly know where to begin. Does anybody have tips or pointers for how to start a group or how to make a GIS group successful?


How to export Symbology from ArcGIS to Microstation or QGIS?

I am using ArcGIS 10.1 and trying to figure out how to export my symbols of elements in ArcGIS to CAD or QGis readable format. Problem is, that I can only export symbology to Style, or export map to DGN, but there will be export just only graphics without symobology. Please, do you know how to solve this? Thanks a lot!


Where is Google Maps extension for ArcGIS?

I need to import Google Maps into ArcGIS for Desktop and apparently I need to activate a Google Maps extension. I'm working on ArcGIS 10.0 and there are no such extension. Do I have to install a new extension?


Utility GIS - Type of data cleanup needed & Effort

I am trying to get into Utility GIS and I have been asked by my boss to size up the effort of what it takes to clean up Pipeline, Utility asset GIS data and figure out if there is a way think of a factor like a per linear mile man days of effort. The data is in the form of ESRI File G Geodatabase. Typical clean up related tasks I have are : 1\. Data consistency by using domains (Manual Effort) - Per Record Multiplier 2\. Topology checks (ArcToolbox) Per Featureclass Multiplier 3\. CAD to GIS conversion (Per Feature Class Multiplier) 4\. Visual QCs Can the hive mind suggest any other tasks typically for cleaning pipeline data and any assessment of what you would think of a man-day multiplier. Appreciate the help!


what raster format for r.los in qgis

What raster format should I use to create viewsheds in QGIS using the grass r.los plugin? I am using the Ordnance Survey Opendata DTM (tried both ascii and shapefile formats) and srtm geotiffs but when I run r.los I can't select an input raster - does that mean the raster is the wrong format?


Clean way to install gdaltools plugin for QGIS on OSX

I have qgis 1.8.0 on osx mountain lion and am installing the gdaltools 1.2.29. I have run into a problem that a number of people are seeing. The plugin cannot find my gdal python (see below). Anyway, I can fix the problem hackily (if that is a word) by symlinking                 cd /Library/Python/2.7       sudo mv site-packages/ site-packages-old       sudo ln -s /Library/Frameworks/GDAL.framework/Versions/1.9/Python/2.7/site-packages /Library/Python/2.7      So that fixes the prob because the python path list that the plugin knows about includes /Library/Python/2.7/site-packages But the real solution is to add a path to the list of paths that the python plugin[s] know about. How do you do that? \----start of error message---- Couldn't load plugin GdalTools due an error when calling its classFactory() method Traceback (most recent call last): File "/Applications/QGIS.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 164, in startPlugin plugins[packageName] = package.classFactory(iface) File "/Users/aberezin/.qgis//python/plugins/GdalTools/ **init**.py", line 32, in classFactory from GdalTools import GdalTools File "/Applications/QGIS.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 309, in _import mod = _builtin_import(name, globals, locals, fromlist, level) File "/Users/aberezin/.qgis//python/plugins/GdalTools/GdalTools.py", line 39, in raise ImportError( error_str ) ImportError: No module named osgeo [python-gdal] Python version: 2.7.2 (default, Jun 20 2012, 16:23:33) [GCC 4.2.1 Compatible Apple Clang 4.0 (tags/Apple/clang-418.0.60)] QGIS version: 1.8.0-Lisboa Lisboa, Python path: ['/Applications/QGIS.app/Contents/MacOS/../Resources/python', '/Users/aberezin/.qgis//python', '/Users/aberezin/.qgis//python/plugins', '/Applications/QGIS.app/Contents/MacOS/../Resources/python/plugins', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python27.zip', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat- darwin', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat- mac', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat- mac/lib-scriptpackages', '/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib- tk', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib- old', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib- dynload', '/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/PyObjC']


How to edit Attributes on WFS Layer through Popup?

I am working on a project for university. So far I have displayed WfsLayers which can be select from a panel, also I can click on my points whereupon a Popup will show up with some information which is derived from my Geoserver which is connected to my Post GIS table. So far so good. I now would love to be able to edit my attributes in the popup box. I found this example (http://gis.ibbeck.de/ginfo/apps/OLExamples/OL26/examples/styles_unique_with_group_wfs.html), which looks just great and is exactly what I am looking for but being a coding newbie in general I have terrible difficulties understanding the code. I tried to use parts but just can't get it to work at all. Does someone have an easier example I could try to work with?               <!DOCTYPE HTML>               <title>OpenLayers Demo</title>


Draw middle ordinate perpendicular to chord, out to curve, with arcpy

I've been assigned to find the radius of road curves. My boss found this article describing a pretty straightforward method of doing so:   1. find the start/end points of the curve;   2. create a chord;   3. create the middle ordinate (from chord midpoint, perpendicular line that goes up to the curve)   4. do a bit of math with the chord and ordinate lengths to find the radius. Manually, this works very well. The downside is that I'm working with thousands of miles of road, so we want to build an automated tool to be able to quickly analyze our entire existing network, and also be able to check roads digitized in future. I have adapted this script to find transition points where the curve begins/ends -- comparing line segments to see if the road direction is changing (e.g. is it still turning right, or has it started curving left instead). I therefore have chords for all the curves, but now I'm struggling with how to make middle ordinates. It needs to:   1. start at the midpoint of the chord   2. be perpendicular to the chord (I'm not positive how to do this yet, but found a script that will probably work and/or be adaptable)   3. be long enough to reach the road, but no longer (???) So, long question short: how can I find the (approximate) distance between my chord midpoint and the road? (Another StackExchange question describes working with the Near_analysis function which looks like it could work, but I currently have all of those chord midpoints just stored as a python list, not a feature class. I'd like to avoid generating a bunch of new features if possible.)


Merge point features within a given radius

I have a file in QGIS with 5000 points. Some of the points are really close to each other (5m). How can I cluster the point features to one point feature if they are within 5m? The new points should be centroids from the given points.


remove whitespace Python code results in odd character - Field Calculator

I need to remove whitespace from the end of a string in a field in ArcMap. I tried first the VBA function RTrim, but whitespace still remained. I thought maybe the whitespaces were tabs, so I tried 'Replace([field], vbTab, ""), but no luck. Then I tried the Python function rstrip, as so: !FIELDNAME!.rstrip() ..but this resulted in the whitespace being replaced by this character: Â There are several of these characters in each of the fields cells trailing the string. I have looked online for some meaning of this symbol, but I haven't found any. Doe anyone know what this means? I have still not successfully removed the trailing whitespace. here is a screenshot of the attribute column after I run the Python code given above: ![enter image description here](http://i.stack.imgur.com/kPe3V.png)


Compiling QGIS 2.0.1 on Centos (SIP problems)

Trying to compile QGIS 2.0.1 on CentOS 6.4 (64 bit):               /usr/share/sip/PyQt4/QtGui/qsyntaxhighlighter.sip: In function ‘PyObject* meth_QgsExpressionHighlighter_setCurrentBlockUserData(PyObject*, PyObject*, PyObject*)’:     /usr/share/sip/PyQt4/QtGui/qsyntaxhighlighter.sip:78: error: ‘qtgui_wrap_ancestors’ was not declared in this scope     make[2]: *** [python/CMakeFiles/python_module_qgis_gui.dir/gui/sipguipart3.cpp.o] Error 1     make[1]: *** [python/CMakeFiles/python_module_qgis_gui.dir/all] Error 2     make: *** [all] Error 2      This is with SIP 4.15.2, Python 2.6.6, and PyQt-x11-gpl-4.10.3 Any help please?


Can I create a table showing all the Srid values for all the tables in my DB?

Hi I can identify the current SRID of one table (at a time) in my database, but I was wondering if there was a way to bring back a single table that would tell me all the datasets that do not presently have a SRID against the geometry? This is in SQL Server 2008. Any help would be greatly appreciated.


Create new geodatabase in modelbuilder and place queried data in new database

![Query existing database and place queried data into newly created geodatabase](http://i.stack.imgur.com/eMPs8.jpg) I have built a model that will take an existing database and query all feature classes in that dataset using a SQL Expression and I then want all the queried data from each feature class to be exported by feature class to a newly created database. I know I can manually create a new geodatabase and then set the output to that the geodatabase in the "Select" tool but I want the model to create the database automatically. The model here shows how I have tried to do it but by doing it the way illustrated here it creates a new GDB for each iteration, exports the feature to the new GDB and then moves on the next iteration, overwriting the previous feature class. How can I create the geodatabase first (once and only once), then run the iteration and place each each iterated feature class into the newly created database? Any advice?


Python Field Calc for Document Hyperlink

I am trying to create a script that will hyperlink to a .tiff or .pdf. These files are stored in one folder on a shared server. I want the script to look in a field: 'PLAT_IMAGE', which has the plat number (ex: 12048-001-01) and find it's match in the folder on the shared server and then write for example: \\\server\engineering\Asbuilts\Plats\12048-001-01.tif Making the document hyperlink active in ArcMap when using Identify My problem is that the calculation is going through and not giving me any errors but yet it's not calculating the field, so I don't know what I'm doing wrong!               import os, arcpy, fnmatch          def  x(PLAT_IMAGE):         path = "J:\\engineering\\Asbuilts\\Plats"         plat = PLAT_IMAGE     for file in os.listdir(path):         if os.path.isdir(os.path.join(path, file)):             for file in os.listdir(os.path.join(path, file)):                 if fnmatch.fnmatch(file,plat + '.*'):                     return "\\" + "\\mercury\gisdata$\engineering\Asbuilts\Plats" + "\\" + file          x( !PLAT_IMAGE! )      Any help is much appreciated


Georeferencer plugin will not load in QGIS

I have downloaded and installed QGIS. I would like to georeference an image using the core georeference plugin. When loading it no toolbar appears and the 'update', 'reinstall' etc boxes in plugin menu are also greyed out. Any clues?


How to get neighborhood values count using a geoprocessing tool in ArcGIS?

I have a raster (black=0, red=1, and white=NoData). I am wondering if there is any tool in arcgis which can help me only count number of 0s and 1s separately in a 3*3 neighborhood. I need an attribute table with that information. Any hint of relevant tool would be appreciative.


How to force ArcGIS server to create KML?

Attempting to serve arcgis server (AGS) 10 sde data to Google Earth (GE).   I don't have any problem opening the service or the network link in GE on my desktop.   When I try to open in ipad device nothing happens.   Finally today I got a message "kmz not supported on iphone". My network links generated by AGS are all kmz.   **Is there a way to force AGS to use kml instead?**   ...or has anyone been successful with AGS service on an ipad? **EDIT:** I do have specific reasons for using GE. I hope to negate those with some time. However everything I have built to this point has been flex. (I was informed that we now have ipad devices).   So either rebuild all my widgets in java, or find a solution that allows kml on the ipad.   Update: I get an error if I try to open kmz,   if I unzip it and put the kml in my network link folder it just doesn't work.   if I try to publish a network link with the kml ags says no kmz found.   I guess I am on to rebuild all my widgets in java.


Alternative to Copy Features in ModelBuilder to display a result

I have created a model that takes in a raster, converts it to a point shapefile, then extracts a set of points based on location, then further selects a set of points based on a SQL expression. I tried to run it once but it was taking absolutely forever. I would have thought this step to be akin to 'Export Data' in the regular window. Is this the right method? Current model creation below. The goal of this procedure is to find a region of high density crops in a selected country. I will eventually need to determine the coordinate footprint of the area so I can export that information to fine tune a forecast model to focus on that particular area. ![Model](http://i.stack.imgur.com/dH0Ou.jpg)


What are some lightweight python modules to convert LAS data formats to text formats?

I searched on the net and found pyLAS, but it seems to require quite a few other modules and I couldn't get the install right. Any pointers will be much appreciated. Thanks!


In Opengeo 4.0.2 suite, enabling labels causes map to show pink screen. The map shows up properly if labels are switched off. I am using Opengeo 4.0.2 suite on Ubuntu 13.10 server. How do I resolve this problem?


It is possible put highcharts chart in openlayers popup?

As title. I build a openlayers map on my website.   and there are many feature vector on the map. now I want to put a highcharts charts in popup.   to show the temperature in every city   something Like this   Is it possible ? How? Thank you! I used contentHTML but it is failed highcharts code               $('#chart').highcharts({         chart: {             type: 'area'         },         title: {             text: 'Area chart with negative values'         },         xAxis: {             categories: ['Apples', 'Oranges', 'Pears', 'Grapes', 'Bananas']         },         credits: {             enabled: false         },         series: [{             name: 'John',             data: [5, 3, 4, 7, 2]         }]     });      part of openlayers code               var content = "

" ;             var popup = new OpenLayers.Popup.FramedCloud("popup",                 OpenLayers.LonLat.fromString(feature.geometry.toShortString()),                 null,                        content,                 null,                 true             );


Error 'GDAL extension not installed' while installing cartaro

I am trying to install Cartaro on my Ubuntu 12.04 LTS. I have followed every step from this guide: http://cartaro.org/documentation/installation But on the 'Verify Requirements' page I get this warning which says ''GDAL extension not installed'. Also when I click on proceed with installation, on Setup database page I see mysql and SQLite as options type not PostgreSQL as shown here Here is what I did: After using               sudo unzip geoserver-2.3.5-war.zip geoserver.war -d /var/lib/tomcat7/webapps/      I opened /etc/default/tomcat7 and added the following code in the end               JAVA_OPTS="${JAVA_OPTS} -DENABLE_JSONP=true"      Then I edited geoserver/WEB-INF/web.xml and added the following code just after web-app starting tag.                 <session-config>         <session-timeout>33333</session-timeout>       </session-config>      Please tell me what did I do wrong and what is the correct way to install the GDAL extension. Thanks


Declination Diagram

1. What point on a UTM map does the declination diagram refer to? Is it the centre? Is it the edges? Is there a standard?   2. The central meridian on a UTM map points to true north. Are grid north and true north arrows identical only at this central meridian? If they are identical then what point on the UTM map does the grid north arrow on the declination diagram refer to?   3. If all grid lines on a UTM map point to true north (causing increased distortion with distance from the central meridian) then why is there a convergence angle (grid north to true north) on the declination diagram?


UTM convertor from one zone to next zone

I want to convert the utm cordinate systenm WGS 84 "E,N " from 38 N ZONE to 39 N or to 37 N. pleas send me the formula and programs that I can do it. thank you


How to modify existing grass gis tools and run them in QGIS?

I would like to take an existing tool of the grass tool box and modify the code to the specific needs of my task. As I´m struggling to start some kick- off would be great. How to access code? What do I have to consider? How to access/run the code in QGIS?


Can I upload an LPA file (Long Profile) to QGIS

I have a .lpa file that I'd like to upload to QGIS - Can this be done and if so, how? Thanks!


How to export qgis as a kml file?

I have a .qgis contour plot.   I would like to convert this to a .kml file.   How specifically would I go about doing this?


gdal_proximity distance to meters

I followed the answer given in this topic for searching the distance to the nearest coastline given a lat/lon coordinate. The distance output obtained is in degree (I think) but I need the distance expressed in meters. By citing the answer given: > There will be a units problem if you create the distance raster in > geographic coordinates. In order to get kilimetres you will need to > rasterize the sea in a suitable projected coordinate system and then project > the distance raster to geograpic coordinates to do the lookup. This will > give you the fastest response time as you're not projecting the query > coordinate at the time of lookup. At the moment my python code is the following:               from osgeo import gdal     ds = gdal.Open("/home/user/distance.tif")          gt = ds.GetGeoTransform()          topLeftX = gt[0]     dimX = gt[1]     topLeftY = gt[3]     dimY = gt[5]          band = ds.GetRasterBand(1)     #Numpy array of the distances in the whole area     data = band.ReadAsArray()     #Indexing by longitude/latitude position            def getDistance(coord):         x = int((coord[0] - topLeftX)/dimX)         y = int((coord[1] - topLeftY)/dimY)         distance = data[y,x]         return distance      The coordinate system is WGS84. I think I have to rasterize the shape file in a Projected Coordinate System but then I don't know what to do. Any suggestions?


Surface distance on DEM between two coordinates

Given two coordinates on a digital elevation model (DEM), how do I compute the actual distance traveled between the two locations, assuming a straight line route? How will changing the resolution of my DEM affect the results?


openlayers image layer stretch with no blur effect

I use an Image Layer to draw some data at specific locations on my web application using gwt-openlayers. Assume the original image is 400x400 pixels. Whem I zoom in the map, the image is stretched. The problem is that the stretch effect is 'blurry'. How can I avoid this blurry effect and instruct openlayers to apply a normal 'raw' stretch on the image? Is there some options somewhere? When I say blurry, I mean this: ![enter image description here](http://i.stack.imgur.com/YVNNg.png) The image above demonstrate the concept. I have a 400x400 pixels image. The image is formed by 4 sub squares 200x200 pixels each. Top left square is blue, top right is red, bottom left is yellow and bottom right is red. When I draw this image as an image layer in openlayers and I zoom the map, the center of the image gets blurry and the four colors are blended. Is there a way to avoid this blur effect and have a 'raw' stretch? Thanks. Here's the code (in gwt-openlayers) I use to create the image layer:               org.gwtopenmaps.openlayers.client.layer.ImageOptions imgOptions = new org.gwtopenmaps.openlayers.client.layer.ImageOptions();     imgOptions.setIsBaseLayer(false);     imgOptions.setTransitionEffect(TransitionEffect.RESIZE);     imgOptions.setAlwaysInRange(true);     imgOptions.setLayerOpacity(0.7);     imgOptions.setProjection("EPSG:900913");          imgLayer = new org.gwtopenmaps.openlayers.client.layer.Image("my_name", null, bounds, new Size(1, 1), imgOptions);     map.addLayer(imgLayer);


How can I relate the osm2po osm_source_id to data in OpenStreetMap?

I created and imported OpenStreetMap (OSM) data into Postgres using osm2po: e.g. `java -Xmx8g -jar osm2po-4.6.9/osm2po-core-4.6.9-signed.jar config=osm2po.config.eu_xc prefix=eu tileSize=x,c europe-latest.osm.bz2` and `psql -d osm2po_eu_xc -q -f "eu_2po_4pgr.sql"` which creates the table _eu_2po_4pgr_ and has datasets with columns _"osm_source_id_ ", _"osm_target_id"_. I created and imported the exact same osm data (using the exact same bz2 file) using osm2pgsql (0.69) e.g. `time osm2pgsql -smuc europe-latest.osm.bz2 -d osm_eu` which creates table like _planet_osm_point_ which has a _"column osm_id"_. I can't find the _"osm_source_id"_ in say _planet_osm_point_. What I would like is to join the osm2po table with the original OSM Data, in order to have the complete information. Is that possible or am I missing something fundamentally?


How to Resize the Qgis Main Window from python

I would like to write a plugin that simply resizes the qgis main window to make it small. How can I access the qgis main window through the qgis api? Use Case: We are having problems with qgis taking over the screen, so that the minimization toolbar is inaccessible, and we can't use anything else on the desktop until we have exited qgis. This seems to be a result of incapatability between Qt and our very old windowing system. I know how to fix this problem using the QGIS.conf file or the registry, but we have many users of varying ability, and many are switching between windows and linux (tunneled through Solaris) servers, which can cause the problem to recur. Fixing the problem with repeated edits of the config file (or the windows registry) is really too much to ask of them. So, I want to make a button that they can just press, and fix it.


GeoWebCache: misalignment of tiles

I'm trying to test simple web map based on PostGIS/Geoserver/OpenLayers stack on Windows Server 2008. When I load my WMS layer by:               var gwctest = new OpenLayers.Layer.WMS(         "gwctest",         "http://localhost:1979/geoserver/wms",         {         'layers': 'Index2000:index_vect_s100',         'format': 'image/png',         'srs': 'EPSG:3857',          'units': 'm',         'transparent': true         },         {'opacity': 0.8, 'isBaseLayer': false, 'visibility': false}          );      Everything seems to be fine. ![alt text](http://i.stack.imgur.com/pdHAh.jpg) However, when I try to add WMS layer to my map using GeoWebCache, replacing               "http://localhost:1979/geoserver/wms"      with:               "http://localhost:1979/geoserver/gwc/service/wms"      My data suddenly move from correct position: ![alt text](http://i.stack.imgur.com/p6H7y.jpg) Any ideas on what may be the cause? EDIT: some additional info on GWC mailing list post.


Error in performing transect in ambur R?

I'm trying to filter a transect, window size is 5 and filter individual baselines is set to 0. I select the transect shapefile, then the inner baseline shapefile and it gives me this error.               Error in SpatialLinesDataFrame(shape.prep2, trandata) : row.names of data and      Lines IDs do not match.      I did everything correct up to this point, add fields to shapefiles, construct transects, stuck on filter transects. Anyone know how to fix this error or knows what it means? The syntax for SpatialLinesDataFrame is               function (sl, data, match.ID = TRUE)     {     if (is.character(match.ID)) {         row.names(data) = data[, match.ID[1]]         match.ID = TRUE     }     if (match.ID) {         Sl_IDs <- sapply(slot(sl, "lines"), function(x) slot(x,              "ID"))         data_IDs <- row.names(data)         mtch <- match(Sl_IDs, data_IDs)         if (any(is.na(mtch)))              stop("row.names of data and Lines IDs do not match")         if (length(unique(mtch)) != length(Sl_IDs))              stop("row.names of data and Lines IDs do not match")         data <- data[mtch, , drop = FALSE]     }     if (nrow(data) != length(sl@lines))          stop("length of data.frame does not match number of Lines elements")     new("SpatialLinesDataFrame", sl, data = data)     }     <environment: namespace:sp>      not too familiar with code sorry. The error is happening because one of those if statements. Looking at the documentation on the R website I see that shape.prep2 is the object in SpatialLines class and trandata is an object in the data.frame class. "object of class data.frame; the number of rows in data should equal the number of Lines elements in sl." It must be something wrong with one of the shapefiles.


How can to merge DEMs in ArcGIS for Desktop 10.2?

I'm trying to merge some dems downloaded from USGS. When I used arcgis 9.2 there was an expression to be used in raster calculator to do that. I used it again but unfortunately it doesn't work in arcgis 10.2. What can I do?? I do not have the toolbar Production Contouring


Generating new shp files according to attribute in QGIS

I know that it is possible to color a vector map by attribute. I separated a map into 4 colors based on an attribute and now I want to generate a new map for each color. Is there a way to do this in QGIS?


Looking for web and desktop tools with GML support

I am a graduate student and a novice in GIS field. I have few GML files that contain vector polygon data and I want to see them in a map format. The problem is that I do not have access to ArcGIS/Grass GIS software installed in my machine. I am wondering if there are online GIS website where I can upload my file and visualize it. Also is there any simple and small size free/open source software that I can install on my computer?


How to import excel point and line data into ArcGIS 10.1 as shapefile?

I have two separate excel documents that I'd like to import into ArcMap, but I don't exactly know how to create a shapefile. The two files are as follows: 1) Point data. This spreadsheet has a long list of points (road intersections) with unique IDs and XY coordinates. This one shouldn't be too hard to import. 2) Line data. This spreadsheet has "From" and "To" columns with the unique IDs from the point data spreadsheet. It also has the XY coordinates of the points. Is there any way to import the line data spreadsheet so that ArcMap draws a line between each of the points? Another thing to consider is that many of these roads are not simply straight lines, so I don't know how to account for that. I'm using ArcGIS 10.1. Thanks!


Agent Analyst - How to reduce memory use during many, many loops?

I am currently using the Agent Analyst adaptation of Repast within the ArcGIS framework. The model I've built is capable of doing what I'd like (ie. randomly generating biologically-sensible migration pathways in order to find NDVI- maximized routes). However, in the loop I'm doing where I make an ArrayList for each randomization ('Positions'), and then iterate this into another, sorted ArrayList ('Movements', I find that after about 25-35k loops (ie. TickCount) the Java program crashes. Upon inspection it looks like the RAM used by Java during this process rapidly grows to take over all available memory. I'm looking for any input into what might be going on that is overloading the RAM (e.g., retained ArrayLists in the cache, 'new' rasters loaded rather than reloaded during each Week iteration). This Noob is very grateful for any guidance. Here is the relevant part of the NqPy code created within the Vector Agent part of the model (where I suspect the problem is occurring):               def step():       #self.model.loadTracker()       tc = self.model.getTickCount()       self.PABUid = int(tc)       self.totalNDVI = 1000000.0       point = self.the_geom.coordinate       point.x = self.model.startX       point.y = self.model.startY       self.setPositions(ArrayList())       positions = self.getPositions()       positions.add("PABU"+self.PABUid)       self.model.iterationWeek = self.model.startWeek        wk = self.model.iterationWeek       eW = self.model.endWeek       while wk <= eW:         self.rotate()         self.distance()         #movement output from the 'distance' action is evaluated to see if the agent would         #move into an inhospitable enviroment (eg. ocean)         self.propMove()         #print "CheckingNDVI"         nX = self.model.newX         nY = self.model.newY         self.model.rasterName = "NDVI"+wk+".tif"         rasterName = self.model.rasterName         raster = self.model.getRaster(rasterName)         val = raster.getPixelValueAtMap(nX,nY,0)         #print "NDVI value: ",val         self.PixelValue = val         if (val > self.model.minNDVI):           self.model.isSuitable = true         else:           self.model.isSuitable = false         self.calcMove()              positions.add("Week: "+wk)         point = self.the_geom.coordinate         positions.add(Double(point.x))         positions.add(Double(point.y))         positions.add(Double(self.PixelValue))         self.totalNDVI = self.totalNDVI+self.PixelValue         #self.model.trackPABU()         wk = wk+1         self.model.iterationWeek = wk         #self.model.writeAgents()         #self.model.updateDisplay()            positions.add(0, Double(self.totalNDVI))       #print positions            movements = self.getMovements()       iterator = movements.listIterator()       iterator.add(str(positions))            #print movements       Collections.sort(movements)       if tc>self.model.SimsToKeep:         movements.remove(0)                 if tc % 5000 == 0:         for strAgent in movements:           file = File(self.model.outputDir+tc+"output.txt")           fw = FileWriter(file.getAbsoluteFile(),true)           bw = BufferedWriter(fw)           bw.newLine()           bw.write(str(strAgent))           bw.close()         print tc


What is most suitable DEM to delineate watershed (from spot height, ASTER and SRTM)?

I have three DEMs of our country which are spot height DEM, SRTM DEM and ASTER DEM. SRTM has 90m resolution and ASTER has 30m resolution. My supervisor told me to delineate watershed using these three DEMs separately, and also using merged DEMs. What is the main differences between the three types of DEM that I mentioned above? How do I merge these three DEMs? What is the most suitable DEM to delineate watershed?


Alternative to ST_ClosestPoint, one which returns a vertex of the geometry

ST_ClosestPoint returns the closest point on a geometry, as in the example at http://postgis.net/docs/ST_ClosestPoint.html What I need is to get the closest point which is _already_ a part of the geometry. Example:               geometry = [1,2,3,4]     Point = any position      `ST_ClosestPoint(geometry,point)` returns something like `[2.5]` which is _between_ existing vertexes. How can I get a result without interpolation (so it should be 1,2,3 or 4 and nothing else possible)? Thanks!


Lost relation when importing from OSM to PostgreSQL with ogr2ogr

I am currently trying to import all administrative boundaries (countries, states, cities, etc) from OSM. I downloaded a small file from Geofabrik (Corsica), in order to run my tests. I understand that I could run the same operation worldwide using planet- latest.osm. First, I extracted all boundaries using the following command like               osmosis --read-xml file="corse-latest.osm" --tf accept-ways boundary=administrative --used-node --write-xml out1.osm      Then, I import everything into PostgreSQL using               ogr2ogr -f PostgreSQL PG:"[connection details]" out1.osm -lco COLUMN_TYPES=other_tags=hstore      Which fills up the multipolygon table with the list of cities, and that is exactly what I wanted. However, if I look for "Ajaccio", there is a lot of useful tags associated with the node 51341115, that do not show up in the multipolygon. Instead, they are associated to a point linked to that node ID. In the XML file, I have been able to link that point to the relation 73283, inside of which are standing all the ways of the city, and the tags that I actually find in the multipolygon table. The question is then the following: how can I link the multipolygon to the point in order to extract the tags, using solely the database? I'm not specifically tied to those tools, I'll use anything that works and that spares me writing code. Many thanks :)


How to display images as cross sections in 3d

I want to display a base map and several vertical cross sections in 3D. The base map image is easy to rectify as a geotiff. How can I also display the cross sections? The cross sections are at present just rectangular images. I can work out the depth/pixel conversion and the surface xy locations for each column of pixles. I have no DTM so the base map is flat and the cross sections simply reside directly underneath. My preferred options are (1) anything written in Python (2) open source GIS software. eg QGIS, uDig etc. (3) perhaps point cloud software like meshlab? Regards


Workflow in ArcGIS 10

Our agency has collabrated with six townships and one city for preparing an inventory of stormwater network. We are using ArcGIS 10 (2 ArcView licenses) and ArcPad 7.1 (stuck with it since our Mobilemapper CX is running on Windows CE). We are using ESRI File based Geodatabase (chose it because it allows multiple readers and one editor per feature datasets and it has one feature dataset for every political subdivision). We have finished 25% of the project. My primary concerns are: a) We cannot edit same feature dataset at a time. Mostly, my intern and I collect and edit data for same area/region at a time. b) Collecting data in the field is not efficient enough since we cannot see each others edits. (using ArcPad Data Manager) c) Sharing data with client is not a smooth task (converting the data into KML format so that it can be viewed on ArcGIS Explorer). I understand that the solution to my problem is ArcSDE and ArcGIS server but we don't have funding. We want to finish the project ASAP. Please recommend a workflow and open source software where we can work and share the data with clients more efficiently. Thank you.


How to use the output file of a script as an input file for the rest of the script in arcpy?

If my question does not make sense or needs clarification please let me know. I am struggling to solve it so I was struggling to phrase it into a question. I have this code here:               import arcpy,os     from arcpy import env          #Name .gdb and set location. This will create the .gdb     out_folder_path = ("C:/Users..../Documents/test")     out_name= "Test.gdb"     arcpy.CreateFileGDB_management(out_folder_path, out_name)          # This will load the .xml into the .gdb created above     target_gdb = "C:/Users/..../Documents/test/Test.gdb"     in_file = "C:/Users..../MASTER.xml"     import_type = "DATA"     config_keyword = "DEFAULTS"          # Execute ImportXMLWorkspaceDocument     arcpy.ImportXMLWorkspaceDocument_management(target_gdb, in_file, import_type)      I am struggling to convert this into a tool. The problem I am encountering is I want the user to just set the .env settings, then input the database name and then location. Then the tool will take over and create the file gdb and then load the .xml workspace into it. The .xml workspace has been hardcoded to a certain location. I am just struggling with the target_gdb for the .xml import. It does work as is with hardcoded locations and for one time use. But for those who do are not python saavy having to change code everytime is not feasible. Again I want this to become a tool. I did find this formum: Automating script to set geodatabase? but this just confused me even more. Updated code:02/4/14               import arcpy,os     from arcpy import env      # Name .gdb and set location. This will create the .gdb               out_folder_path = arcpy.GetParameterAsText(0)     out_name= arcpy.GetParameterAsText(1)     arcpy.CreateFileGDB_management(out_folder_path,out_name)      # This will load the .xml into the .gdb created above               target_gdb = out_folder_path + os.path.basename(out_name)     in_file = "C:/Users/.../MASTER.xml"     import_type = "DATA"     config_keyword = "DEFAULTS"      # Execute ImportXMLWorkspaceDocument               arcpy.ImportXMLWorkspaceDocument_management(target_gdb, in_file, import_type)


Is there a Deployment build of GeoExt 1.1?

The latest version of GeoExt that is available for download at http://geoext.org/downloads.html (version 1.1) appears to be a development version. They are individual files with whitespace and comments. Usualy there is a minified, compress version of javascript libraries available, to be used in production environment. Nowhere could I find such a version for GeoExt 1.1. Is there a deployment build available?


Removing duplicate rows from a shapefile

I would like to write some c# code to delete duplicate rows from a shapefile stored in a folder. Can anyone provide me with an example or pointers on the best way to attempt this, please?


Identifying Features With a Relationship Class

I have a feature class which represents points where vegetation clearing has occurred. I also have a table which has a list of vegetation species which were cleared at each location and the amount of vegetation which was in the sampled area. I have a relationship class in the geodatabase which links the two based on a `Clearing_ID` field. When the user clicks on a point with the Identify tool, they can drill down to see the list of plant species present. However, the list appears as numbers which is the `OBJECTID`. Ideally, I'd like the user to see the list of species present as the default field instead of the object ID. ![enter image description here](http://i.stack.imgur.com/v9YZt.jpg) Is there a way to do this?


Google Maps Markers inconsistent as they're drawn?

I've playing around with the google maps api and am puzzled at the following behavior. If I use `mPoint` as the LatLng for my marker, the marker is rendered on a different point on the map as opposed to putting the same value directly into the properties of the marker. Code chunk is as follows:                       var mPoint = [new google.maps.LatLng(38.991300,-76.936165)];     //        var GPoint = new google.maps.Marker({map:map,position:AVPoint});       var marker = new google.maps.Marker({         map: map,         position: new google.maps.LatLng(38.991300,-76.936165),         //position: new google.maps.LatLng(mPoint),         draggable: true       });      Why would this happen?


How can I develop a vehicle tracking application on a tight budget?

I need to develop and GIS application in .net with following Requirement, Can you please suggest the options how can i start with. I am new to GIS World, but I have a good knowledge. I neither developed any GIS application before! 1) We need to track vehicles on the map, the data recieved from Multiple Controllers connected Over Ethernet. 2) Display of Movements of Multiple Vehicles to be Provided for simultaneous operation. 3) Displaying Latitude, Longtitude, Altiitude, Roll, yaw, Pitch, Velocity with Every Vechicle. 4) Unique Icons for Multiple vehichles (icon should be configurable) 5) Provision for the modes to facilitate selection of vehicles for simultaneous display on Map. 6) Automatically move the Vehicle icon representing the respective vehicle, to the correct geographic location on map when vehicle change the position. 7) Automatic PAN movement of the area concerned on the digital map during Online Vehicle Tracking. 8) Calculate Actual disptance traveled by vehicle along with Aerial Distance, 9) Calculate relative distance between two vehicles. 10) Digital Map retrieval and display as per the area of operation. 11) Provision for additional icons creation on the map. (These icons should be displayed for futher use) 12) Provision of Information Layer (Peliminary/Detailed information about Terrain etc), 3D ANALYSIS 13) Geo-Referencing Facility. It should be able to digitize it by georeferencing. 14) Import digital masp irrespective of their original/native format. 15) OS Platform independency.   1. It should support Military grid coodrinate Data. 17) Provision for Defininf Waypoints with positioning Information 18) Multiple Access to various Digital map,image map like geotiff,.jpeg,.dted image 19) Integration of scanned maps, satellite images . 20) Zoom in Zoom out, Zoom in Max Scal1 1:1 and zoom out 1:50000 21) User Specified scale and Dynamic display of scale in the Mapview 22) Exchange of GIS Data for generation of Reports. 23) Import/Export Maps using .dgn (Microstation DGN) format.


How do I get the path of a layer file (.lyr)?

In C#, how do I get the path of a layer file? I can get the datasource and name of the layer but I need to get the path to the .lyr file. I want to record the path so my extension can copy and load the .lyr file later. This gives the path to the datasource (.mdb geodatabase in my case) that the .lyr file references:               var dataset = (IDataset) map.Layer[0];     var name = dataset.Workspace.PathName + "\\" + dataset.Name      Using that info plus the FeatureClass name I can open the layer and get the layer's polygons on the map but it is not the same as opening the .lyr file.


Why does field calculator expression sometimes return a binary value (0 or 1)?

I am running simple field calculator expressions and experiencing some of the silliest results I have ever seen. I am working on census data (living as a shape file) and trying to perform a simple [fieldName]/[fieldName] calculation on a newly created field. ArcGIS happens to enjoy playing around with me... and returns values 0 and 1 in no particular order/reason that I can really justify. Both input fields from the field calculator expression are of double data types. I have created a new field of double data type as well. However, the problem persists with other data types, including long/short integer. Is this a bug? Has anyone else experienced this? Shortest route to create + update a new field other than moving between SPSS/Excel and ArcGIS all day long?


How can I use different config.yaml files for application using mapfish print in geoserver?

System: Ubuntu 10.04 server / Tomcat6 / geoserver2.1.2 with integrated GeoWebCache I have various webmapping applications running which utilize the mapfish print module.Some are in tomcat6/webapps some are in my apache2 DocRoot. I need to have various mapfish print layouts to show different logos, legends, scalebars etc depending on which webmapping application is using it....but I only have one config.yaml file to style the pdf layout. Qustion: How can Have I various versions of config.yaml and link them to the different webapps? I´d be grateful for any replies, yours, Rob


How to create a maplayer from mapfile?

I am new in Mapserver, now I want to create maplayer , so its required shapefile for this, and if yes how can i used it. Map file is given below               eg2.map     MAP      NAME ENG1      STATUS ON      SIZE 800 500      IMAGETYPE GIF      IMAGECOLOR 240 240 240      SHAPEPATH "/www/geog5870/geog5870test/data"      EXTENT 87000 8000 660000 660000      UNITS METERS      WEB      IMAGEPATH "/www/tmp/"      IMAGEURL "/tmp/"      END      LAYER      NAME "EnglishGOR"      TYPE POLYGON      STATUS DEFAULT      DATA "England_gor_2001_area"      CLASS      NAME "English Government Office Regions"      OUTLINECOLOR 60 60 60      COLOR 255 255 0      SYMBOL 0      END      END      END


Transform an Irish grid reference (29903) into UK grid reference (27700) on insert

I have an insert script to input several hundred points, however some of the points use the Irish Grid Reference system. I have separated these out however was unsure how to convert or transform? these into 27700 on insert. I have the line below that grabs the UK reference, and the POST_GIS column is set to 27700. `st_geomfromtext('POINT(" . addslashes($data[6]) . ")',27700)`


Customise the splash image of QGIS

Where can I find the splash image which is shown while starting QGIS, and how to customise it?


What projection appeared in the NYTimes this morning?

![enter image description here](http://i.stack.imgur.com/GmSG6.png) Perhaps not as GIS related as some, but definitely about maps. Can anyone identify the map projection that appeared in the New York Times this morning? It's not one I've seen before, and a quick scan through _Map Projections: A Working Manual_ did not yield the result I was looking for.


Programmatically start and stop map services with Python

Is it possible to programmatically stop and start ArcGIS Server map services using Python?


Reusing Layer Selection on Toolbar

I'm working on a custom toolbar for an extension which will have operations that work upon a specific, user-selectable layer. Is it possible to directly use the layer selector used by ESRI's toolbars in your own toolbar? For example, the 3D Effects toolbar contains a "Layer:.." combo box for selecting the layer on which to operate. It includes a full tree of the scene, filtered to certain types of layers. Is this control a reusable object?


How to generate a grid of horizontal equal-area rectangles in a Mercator projection?

I'm new to GIS so pardon my ignorance and bear with me. I'm generating a grid relative to a LonLat point. To make the cells all equal area, I'm creating it in 2163 projection. The problem is that when I reproject to 900913 and stick it on a Google Maps baselayer in OpenLayers, the whole grid is rotated clockwise. I understand this is a purely cosmetic issue, but it's a project requirement. Is it possible to change the map's projection to display it vertically? If not, is there another map source I can use that uses 2163 or can be projected to it? If all of this is out of a question, is there another way to generate equal area grids for Mercator 900913 which will line up vertically?


Does PostGIS not work with psql?

I have a sql script that creates several tables, a trigger, and a trigger function. All the commands work in pgAdmin when I execute them in the "Query" section. My goal is to put all these commands in a script so that I can set up the database on any machine relatively easily and quickly. When I attempt to run the script through psql, however, I get the following errors on all of the lines where I attempt to use `Geography` types:               ERROR:  type "geography" does not exist      It doesn't make sense because I can copy and paste the exact same queries into the query editor of pgAdmin and it works just fine. Does psql not support PostGIS? Or maybe just not the most recent version of PostGIS, 1.5.3, which includes the new `Geography` types? Is there something I can do to fix this? An example of one of my queries:               CREATE TABLE source_imagery (         id SERIAL PRIMARY KEY,         image_type VARCHAR(1000),         image_path VARCHAR(1000),         boundary GEOGRAPHY(POLYGON, 4326),         image_time TIMESTAMP,         catalog_time TIMESTAMP DEFAULT CURRENT_TIMESTAMP     );      Error message:               ERROR:  type "geography" does not exist     LINE 5:     boundary GEOGRAPHY(POLYGON, 4326),                          ^


ArcPad 10.0.2 - Removing shortcuts from ArcPad today

I have created a shortcut that goes to ArcPad today menu, but now I want to remove it. How does one do this?