Data match-up problem ArcGIS 10.1

I am running into a problem with two data sets matching up, and although it seems like a newbie coordinate system problem, it has me stumped. One data set is a GPS track (Garmin Etrex Legend 30) and the coordinate system is set to WGS 1984. It shows what I believe to be the correct decimal degrees in the bottom right corner of the screen. This is a geographic area in Africa (Malawi) that I'm not familiar with so I'm not sure if it's the correct decimal degrees, but when I move the cursor around there isn't too much variation, suggesting it's correct. ![enter image description here](http://i.stack.imgur.com/Hm32J.png) The other file was obtained from a state agency and I've used the "Define Projection" tool to use the same coordinate system as the GPS file. It's decimal degrees are very different and seem incorrect due to their discrepancy from the other data set. ![enter image description here](http://i.stack.imgur.com/fm3oJ.png) When I open the properties for both files and look at the coordinate system they appear the same. Here is the properties for the GPS file ![enter image description here](http://i.stack.imgur.com/GOEsh.png) and here is the one for the State's contour file. ![enter image description here](http://i.stack.imgur.com/JPQ12.png) Any suggestions as to what is likely an embarrassingly simple coordinate system issue would be greatly appreciated,


Definition Query then Zoom to Layer Automation?

I'm trying to create something in here. I want to have an automated process of doing a DEFINITION QUERY then ZOOM TO LAYER with one step. the concept is like that, a dialog comes up with one form that asks me to put a unique id for the specific parcel that I'm trying to look up, then I type the parcel numder then hit OK then it brings me to that specific parcel directly.   1. Dialog pop-up   2. Empty form to be filled   3. Hit ok   4. It zooms directly to that parcel I'm trying to find a way to do that either with ModelBuilder or Python! An automated way to do the definition query then zoom to layer.


Linebreak in style.Text in Openlayers 3

I have a vector layer from json, and have a style symbolizer:               new ol.style.Text({       text: ol.expr.parse('name')     })      In an vector layer, some of the nodes contains more than one name. I try to find a way to split the names on lines.   * Openlayers seems to ignore '\n' in the json.   * I can't find a way to edit the string, for example with replace(re, "\n")


Convert RT90 into WGS84

Goal:   Gain WGS84 decimal to locate some cities in Google Map. Problem:   I don't know how to convert the value of the RT90 coordinate system into WGS84 in SQL Server.


Convert classified raster map back into vector classes

I have this map (have whole country, this is just a sample, original map also contains color symbology for the value classes) and need to attach the values from it to points. I have no problem to georeference it in QGIS but there is problem with the values. I see two uncomfortable possibilities: > 1) attach the values manually using the map legend. > > 2) manually vectorize the map in QGIS and add the value classes But it would be nice, if I could somehow (semi-)automatically vectorize the map using the color (RGB?) values to build the vector classes in shapefile format. Is this possible? I have here SAGA GIS, QGIS and I have also access to one PC with Arcgis 10 with Spatial and 3D analyst. Is it possible or do I have to use the manual way? thanks PS: it is almost impossible to obtain source data (will try, but I am pessimistic in this way)


OpenLayers - WMSGetFeatureInfo, picking a feature's attribute from a WMS layer created in GIS Cloud

I created a WMS layer in GIS Cloud app. I use GIS Cloud because they provide prepared app to transfer data from mobile to layer. The URL of WMS layer: http://www.giscloud.com/wms/67edc3492aa2a8ec445ea34034b0d5b7?. The name of layer:138865:mobile_fleet. I'm working on thesis and I have to create a web mapping app. I decided to use OpenLayers library but I'm newbie, I have never been programming and I don't have much time to understand all aspects. Nevertheless I tried to follow examples on the web of OpenLayers and decided to use example Feature Info Example as template and then I tried to change parameters of functions and objects. I was successful with creating map viewer and adding WMS layer and moreover I added a WFS layer. I wasn't successful with picking a feature's attributes of WMS layer. The code:                   <!DOCTYPE html>     function init(){     // Map is in mercator this time, so over-ride the default     // options that assume lat/lon.     var options = {       projection: new OpenLayers.Projection("EPSG:900913"),       displayProjection: new OpenLayers.Projection("EPSG:4326"),       units: "m",       numZoomLevels: 10,       maxResolution: 156543.0339,       maxExtent: new OpenLayers.Bounds(-20037508, -20037508,                                               20037508, 20037508.34)};         // Create the map object      map = new OpenLayers.Map('map', options);              // create Google Maps layer     var gmap = new OpenLayers.Layer.Google(           "Google Streets", // the default           {'sphericalMercator': true, numZoomLevels: 20}     );          // create Google Satellite layer     var gsat = new OpenLayers.Layer.Google(         "Google Satellite",         {type: G_SATELLITE_MAP, 'sphericalMercator': true, numZoomLevels: 20}     );          var wmslayer = new OpenLayers.Layer.WMS(         "Canadian Data",         "http://www.giscloud.com/wms/67edc3492aa2a8ec445ea34034b0d5b7?",         {             version: "1.1.1",             layers: "138865:mobile_fleet",             format:"image/png",             srs: "EPSG:4326",             transparent: "TRUE",              },         { isBaseLayer: false }     );                 var wfs = new OpenLayers.Layer.Vector(          "Stavros Features",          {             strategies: [new OpenLayers.Strategy.Fixed()]             , projection: new OpenLayers.Projection("EPSG:4326")             , protocol: new OpenLayers.Protocol.WFS({                 version: "1.1.0",                 url: "http://ec2-46-137-23-35.eu-west-1.compute.amazonaws.com:8080/geoserver/wfs",                 featurePrefix: 'Fleet', //geoserver worspace name                 featureType: "vzorek_WGS84", //geoserver Layer Name                 featureNS: "http://fleet.opengeo.org", // Edit Workspace Namespace URI                 geometryName: "the_geom" // field in Feature Type details with type "Geometry"                 })         });          highlightLayer = new OpenLayers.Layer.Vector("Highlighted Features", {         displayInLayerSwitcher: false,          isBaseLayer: false          }     );              infoControls = {             click: new OpenLayers.Control.WMSGetFeatureInfo({                 url: 'http://www.giscloud.com/wms/67edc3492aa2a8ec445ea34034b0d5b7?',                  title: 'Identify features by clicking',                 infoFormat: "text/plain",                 layers: [wmslayer],                 queryVisible: true             }),             hover: new OpenLayers.Control.WMSGetFeatureInfo({                 url: 'http://www.giscloud.com/wms/67edc3492aa2a8ec445ea34034b0d5b7?',                  title: 'Identify features by clicking',                 layers: [wmslayer],                 hover: true,                 // defining a custom format options here                 formatOptions: {                     typeName: 'water_bodies',                      featureNS: 'http://www.openplans.org/topp'                 },                 queryVisible: true             })         }                   map.addLayers([gmap, gsat,wfs,wmslayer, highlightLayer]);          for (var i in infoControls) {              infoControls[i].events.register("getfeatureinfo", this, showInfo);             map.addControl(infoControls[i]);          }              map.addControl(new OpenLayers.Control.LayerSwitcher());              infoControls.click.activate();         map.zoomToMaxExtent();     }          function showInfo(evt) {         if (evt.features && evt.features.length) {              highlightLayer.destroyFeatures();              highlightLayer.addFeatures(evt.features);              highlightLayer.redraw();         } else {             $('responseText').innerHTML = evt.text;         }     }          function toggleControl(element) {         for (var key in infoControls) {             var control = infoControls[key];             if (element.value == key && element.checked) {                 control.activate();             } else {                 control.deactivate();             }         }     }          function toggleFormat(element) {         for (var key in infoControls) {             var control = infoControls[key];             control.infoFormat = element.value;         }     }          function toggleLayers(element) {         for (var key in infoControls) {             var control = infoControls[key];             if (element.value == 'Specified') {                 control.layers = [water];             } else {                 control.layers = null;             }         }     }          // function toggle(key     </script>      You can see this app: http://ec2-46-137-23-35.eu- west-1.compute.amazonaws.com:8080/geoserver/www/test3.html. When I click on some feature of WMS layer (the red points on the map) I see in Firebug that the response is bad. After s few hours I have found solution. For information: I'm using OpenGeo product deployed on AWS EC2. Geoserver has already implemented proxy so you don't need copy file proxy.cgi from openlayers. All you have to do is open webUI (dashboard), click on proxy and then add a hostname (in my case www.giscloud.com) and in some case mimetype too. That is all.


What data is sent from GPS satellites to receivers?

I know more or less how the GPS system works. But I would like to understand it in more details. Exactly what data is sent from a GPS satellite to a GPS receiver? I think there is multiple types of data, since when I start up a gps receiver, it takes some time before ut becomes useful. As what I know this is because the receiver waits on some data that aren't received very often. I guess that the GPS signal that is used for the position is sent more often, but I don't know. Is there anyone that can explain this in more detail? Or have a useful link?


Simplifying adjacent polygons

I encountered a problem simplifying set of polygons that are adjacent. If I simplify each polygon separately with the Douglas–Peucker algorithm (which is used by many open source tools), the resulting polygons are usually not adjacent anymore. This problem exists, for example, when simplifying borders of countries/provinces. Does anyone has a good solution for it, preferably to be done in postgis?


Simplifying polygons to linestring

I would like to simplify some polygons that represent rivers into linestrings. I assume that holes in polygons can be omitted. Anybody has some good idea how to simply do it? It would be nice, if it was possible to be done directly in postgis or by another open source tool, but brief algorithm description will be sufficient.


Can someone explain in simple terms the licence on Google “satellite” images?

What can I and can I not do with the Google satellite images?


How does a quad-tree cope with very large objects?

If I use a quad-tree to index my map data, how does it cope with a road that goes from one end of the map to the other? What other sort of indexes should I consider to get round this problem?


What are the licence limitations on the free Ordnance Survey (OS) data sets?

The OS has some “free” map data, but what am I allowed to use it for?


Are there any other providers of map data in the UK apart from the Ordnance Survey?

The Ordnance Survey charges a lot for most data, so what are the other options?


In the UK, how do I convert from an address to a lat/long?

Please include where I can buy the dataset and how much it costs?


Projection problem with ESRI JSAPI and ArcGIS map service

I was trying to obtain a gpx file with some coordinates by drawing on an Openlayers map with an ArcGIS baseMap. When I draw the polyline and create the gpx, if I open it on Google Earth, what I see is not what I drawed before, the line is totally different from the original and not positioned where I drawed it. I know it's a projection problem, I've tried trasforming the geometry object from Mercator to Geographic, also getting directly the geographic coordinates from the map coordinates, but nothing. I tried to set "spatialReference" to 4362 and then to 3857, but nothing changes. I'm going to use that .gpx on a gps device (the next week I'll go to the Svalbard islands and I need some gps tracks to go around Longyearbyen by snowmobile, there there aren't any sign of life out the town, so I must be prepared to it), when I'll be there I'll adjust the output right for the device they will rent to me, but now I need to save on the .gpx file almost the right coordinates. I'm getting from the map those coordinates: lat: 61.22582068741976 lon: 4.684820015391338 when I'm expecting instead something around 78. lat and 15. lon. This is some of the code I use to create the map (I'm not pasting the code I know it's not responsible of my problems):                   var initialExtent = new esri.geometry.Extent({"xmin":505615.5801124362,"ymin":8678955.717187276,"xmax":525935.6207525175,"ymax":8689168.654279819,"spatialReference":{"wkid":32633,"latestWkid":32633}});              map = map = new esri.Map("map", {extent: initialExtent, logo : false});         basemapURL = "http://geodata.npolar.no/ArcGIS/rest/services/inspire1/NP_TopoSvalbard_U33_CHL/MapServer/";         map.addLayer(new esri.layers.ArcGISTiledMapServiceLayer(basemapURL));      Here I'm using wkid 32633 that is the default for that map, tried to change with known ones, but nothing happened. And now the code I use to get the coordinates:                       dojo.connect(tb, "onDrawEnd", getCoords);             function getCoords(geo) {                 var r=confirm("Salvare tracciato?");                 if (r==true) {                     geo = esri.geometry.webMercatorToGeographic(geo);                        for ( var path = 0; path < geo.paths.length; path ++ ) {                                            for ( var pt = 0; pt < geo.paths[path].length; pt++ ) {                             tra += '<wpt lat="'+geo.paths[path][pt][1]+'" lon="'+geo.paths[path][pt][0]+'"></wpt>';                         }                        }                 }             }      "tra" is a variable that stores all the code I'll insert into the gpx file with an other function. The "webMercatorToGeographic" function transform the map coordinates to geographic ones.


How to set up parameters in Python toolbox tools?

I build a model builder and export it to python. How can I define parameters for python script in python properties? May you explain about Data type in script proprties? thanks alot


What are the units for Garmin's GPX gpxtpx:speed?

I have a Garmin Nüvi GPS and it's filled with XML records that look like this:               <trkpt lat="38.900207" lon="-77.008990">     <ele>11.67</ele><time>2013-06-11T04:04:49Z</time>     <extensions>        <gpxtpx:TrackPointExtension><gpxtpx:speed>6.86</gpxtpx:speed>         <gpxtpx:course>271.06</gpxtpx:course>        </gpxtpx:TrackPointExtension>        </extensions></trkpt>      (line breaks added for legibility) My question — what are the units for `gpxtpx:speed` ? I've examined all of the trackpoints for my GPS and the `gpxtpx:speed` never gets much above 30, and I KNOW that I'm going faster than 30MPH. I've scanned the documentation and I can't find what the units are.


Get Features From GeoServer Server Using Attribute

Can anyone explain how to query features from GeoServer into an OpenLayers map using an attribute value? Here is what I need to do; 1) A user enters a number or multiple numbers for polygons they want to see. This is a friendly number not an objectid or fid . 2) Query the GeoServer WFS or WMS using the user entered friendly number to get the polygon or polygons they want. For instance the user may have a polygon number of GB2546 and GB2546 is in the STAND attribute of my layer in GeoServer. 3) Add the polygon features to a vector layer. Adding the features (once I get them) into my vector layer is easy but I don't know how to make the request to GeoServer to get the polygon feature information. Is there a control in OpenLayers that will get just the features I need from GeoServer? I know this should be really simple but all the examples I can find demonstrate how to get features from layers that are already loaded into an OpenLayers map. Any help would be appreciated, thanks!


Why does OGR driver.Open() for File Geodatabase give AttributeError: NoneType object has no attribute Open?

There are a number of existing Questions that hint at a possible Answer to this Question, but none seem to be providing me with the solution to my problem:   * How to access feature classes in file geodatabases with Python and GDAL?   * How to add GDAL to existing Python installation?   * Gdal Python Bindings and File Geodatabase API To get **GDAL/OGR 1.10** working on a **Windows 7 SP1 (64bit)** laptop with **Python 2.7.3** from **ArcGIS 10.2 for Desktop** I did the following steps:   1. Installed ArcGIS 10.2 for Desktop which installed Python 2.7.3 - there are no other versions of Python present   2. Downloaded and installed `gdal-110-1600-x64-core.msi` (compiled 2/6/2014 6:11 AM) from http://www.gisinternals.com/sdk/PackageList.aspx?file=release-1600-x64-gdal-mapserver   3. Downloaded and installed `gdal-110-1600-x64-filegdb.msi` (compiled 2/6/2014 6:17 AM) from http://www.gisinternals.com/sdk/PackageList.aspx?file=release-1600-x64-gdal-mapserver. As expected this placed an `ogr_FileGDB.dll` in C:\Program Files\GDAL\gdalplugins.    4. Added a new system environment variable called `GDAL_DRIVER_PATH` and set it to C:\Program Files\GDAL\gdalplugins, and I appended (using ';') the same folder onto the end of my `Path` variable. This was from advice offered in one of the Questions cited above.   5. Downloaded and installed `GDAL-1.10.1.win32-py2.7.exe` from http://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal This enabled me to run some Python scripts from IDLE like the one below to count the number of features in a shapefile.               import os     from osgeo import ogr          daShapefile = r"C:\temp\test.shp"          driver = ogr.GetDriverByName('ESRI Shapefile')          dataSource = driver.Open(daShapefile, 0) # 0 means read-only. 1 means writeable.          # Check to see if shapefile is found.     if dataSource is None:         print 'Could not open %s' % (daShapefile)     else:         print 'Opened %s' % (daShapefile)         layer = dataSource.GetLayer()         featureCount = layer.GetFeatureCount()         print "Number of features in %s: %d" % (os.path.basename(daShapefile),featureCount)      However, when I tried to run the code below:               import os     from osgeo import ogr          daFileGeodatabaseFC = r"C:\temp\test.gdb\testFC"     daFileGeodatabase = r"C:\temp\test.gdb"     featClass = "testFC"          driver = ogr.GetDriverByName('FileGDB')          dataSource = driver.Open(daFileGeodatabase, 0) # 0 means read-only. 1 means writeable.          # Check to see if feature class is found.     if dataSource is None:         print 'Could not open %s' % (daFileGeodatabase)     else:         print 'Opened %s' % (daFileGeodatabase)         layer = dataSource.GetLayer(featClass)         featureCount = layer.GetFeatureCount()         print "Number of features in %s: %d" % (os.path.basename(daFileGeodatabaseFC),featureCount)      I got an error: > Traceback (most recent call last): File "C:\temp\GetFileGDBFeatureCount.py", > line 17, in dataSource = driver.Open(daFileGeodatabase, 0) # 0 means read- > only. 1 means writeable. AttributeError: 'NoneType' object has no attribute > 'Open' After some more reading I decided to use `ogrinfo --formats` in a DOS window to see if FileGDB was recognized and it is! ![enter image description here](http://i.stack.imgur.com/WjL3L.png) **Am I missing a step or two that would complete enabling OGR FileGDB support for my GDAL/Python/ArcGIS configuration?** I do not class myself as a developer so compiling a new version of GDAL, or starting from scratch and installing OSGeo4W instead are options I would prefer to avoid, if at all possible.


How to get Waypoints between two coordinates from Google maps or other map applications?

Suppose I have A (source point) and B (destination point), then I want to get a series of route point or waypoints between A to B from Google maps or other map applications. The result that I need is something like this: (LonA, LatA) (Lon1,Lat1) (Lon2,Lat2) (Lon3, Lat3)................(Lon_n, Lat_n) (LonB, LatB) Are there any available methods or tools to do that? Thx.


Is it possible to cut/filter all layers to an area of interest in one step?

The problem of working with a lot of data (about 100 layers). Very slow computer, rendering long. Vector -> Geoprocessing Tools -> Clip a very good tool, but replace the 100 layers is very high and can go wrong. In addition to work with the same data in different regions, and each time the clip function is complicated work. Is there a faster way, and saves computer resources? On ArcView GIS 3.2 is called the "Area of interest" or Autocad Map 3D "query to filter date".


How can i clip vector layer in qgis and python?

I want to write a python function to create a shape file layer from clipping an input layer using a mask. Does someone know the way how to do it?


Is there a way to execute FME transformers from a PythonCaller script in an FME workspace?

I have a very large and complex Python/arcpy script that does lots of custom calculations, and I am trying to port it to an FME environemnt. The script runs beautifully under ArcGIS. The problem in porting is that some of the ESRI-specific functionality, especially geoprocessing tools, has to be replaced with similar FME functionality. **So, I would like to inquire if there is a way to run FME transformers, e.g. HullReplacer, from a PythonCaller script in an FME workspace.** True, the project could be rewritten completely in Workbench. I have no doubt about that; however, that would take a lot of time and effort ... and to be completely honest, I am doing the port to FME for someone else who actually has the FME license. I do not have access to FME myself. Most of the functions of the existing script have ported almost without change to FME Workbench. Works great. There are a few places, though--right smack dab in the middle of the code--where some FME functionality is needed to replace some ESRI functionality, and "there's a transformer for that." If I could execute a transformer in those places and continue processing the existing code with the output of the transformer, I could easily port the rest of the script. I just don't know if it is possible programmatically to execute a transformer from Python. I can't find any direct mention of this in online FME documentation or in online searches. I have noted that a Python wrapper has been created for the IFMELibrary. That provides a promising hint, but frankly I don't know what the IFMELibrary is or what it entails. There is no documentation for it. Can anyone shed any light on this for me? Thanks, all, for the help. ADDITION: Hmmm. I see in this post some sample code (in the first block) showing a transformer being called through Python. Is this a valid and supported syntax? If so, is there an online documentation set that covers the Python/fmeobjects syntax for calling each of the transformers? CONCLUSION: After a lot of searching, testing, and -- finally!-- being able to submit a question to Safe Software, my conclusion is that one cannot call FME transformers from within a Python script in a PythonCaller transformer. This is disappointing. I must proceed with chopping the scripts into little pieces and redoing the whole process in FME Workbench.


A python tool like GNU's diff for shapefiles?

I often use the GNU `diff` tool to compare two files or even two folders for differences. I was wondering if there is a Python based tool out there that does basically the same for shapefiles? In the first place I am interested in attribute comparison but to find also the differences in geometry would be very handy. Maybe someone has a script at hand I could adopt from? I do not have access to ArcGIS so it needs to be a Python standalone solution.


How to turn all lines into one direction?

How can I turn all features of my polyline (postgis and/or shapefile) into **one** direction - it doesn't matter what direction. I know how to switch the direction of lines but I have no idea how to bring all lines (I have thousands of them) into **one single direction**. ![enter image description here](http://i.stack.imgur.com/kpKVu.png) EDIT: A line can have only two directions, the one or the other like 0 or 1. I want them all to have the direction 0 OR 1. For example the left line in the image is in direction 1 the middle line is in direction 0 and the right one is in direction 1. Now I want the middle line to direction 1 OR the right and left line to direction 0. So that all lines have the same direction 1 OR 0.               geom;given_direction;wanted_direction     line;1;1     line;0;1     line;1;1     line;1;1     line;0;1     line;1;1


How to get Mapnik running on Win?

I have installed Mapnik-2.0.1rc0 in windows 7 32bit machine. I then tested the installation by `import mapnik from *` and it did not return any error. But when I tried to input the config version checking by `mapnik-config -v`, it returns a name error. When I tried to run `rundemo.py`, the following error occurred: > the mapnik library and python bindings must have been compiled and installed > successfully before running this script This probably happened due to improperly set `PATH` and `PYTHONPATH`. **How to add PythonPath?** Thanks


Spatial reference for dummies?

Can anyone recommend a good reference or tutorial that explains the principles of spatial references, projections, and coordinate systems? We are a small team of programmers working with GIS technology, and we often get snagged when working with data in different SRS/CS because we don't understand the concepts.


What is the SRS code for Lambert Conformal Conic projection and GCS Everest Definition 1962?

I have a shapefile with Geographic Coordinate System is "GCS_Everest Definition 1962" and the Projected Coordinate System is "Lambert_Conformal_Conic". I want to add this shapefile in Geoserver. But I don't know which code is used as the SRS for this Projection System. I don't know what to put in "Declared SRS" field. Here is the WKT               PROJCS["Lambert Conformal Conic",        GEOGCS["GCS_Everest Definition 1962",          DATUM["Everest Definition 1962",            SPHEROID["Everest Definition 1962",      6377301.243, 300.8017255654936]],          PRIMEM["Greenwich", 0.0],          UNIT["degree", 0.017453292519943295],          AXIS["Longitude", EAST],          AXIS["Latitude", NORTH]],        PROJECTION["Lambert_Conformal_Conic_2SP"],        PARAMETER["central_meridian", 68.0],        PARAMETER["latitude_of_origin", 32.5],        PARAMETER["standard_parallel_1", 35.31445944444445],        PARAMETER["false_easting", 2743196.4],        PARAMETER["false_northing", 914398.8],        PARAMETER["scale_factor", 1.0],        PARAMETER["standard_parallel_2", 29.655283888888892],        UNIT["m", 1.0],        AXIS["x", EAST],        AXIS["y", NORTH]]


How to find features with NULL values?

I have a layer to analyze and I want to create a thematic map with rule based feature. There is a strange behaviour. I want to find all objects with no data (NULL). All empty objects. The zero value for me is a value. I write: > namefield = '' The result is zero but I have 47 objects with NULL value. This is an error. The strange behaviour is if I write: > namefield != '' The result is 56 objects. This rule works. It is impossible to find objects with no data but qgis find only objects with a value. Is this a bug or a feature? It is a bug for me! I use qgis 1.8 Thanks


How to convert an Erdas training set to .csv?

I am working with Erdas Imagine's Signature Editor to perform maximum likelihood classification. I need to use the same training set created in Erdas for other classification using R. Is there a way to extract the raster band pixel values under the .aoi polygons and export that information to a .txt or .csv? I am looking for the same functionality as the ArcGIS Sample tool where the individual raster band pixel values are extracted underneath points and written to a table (e.g. Figure 2). Alternatively, is it possible to export the signature editor table (i.e. including R,G,B,NIR bands and the class label) to a .csv? Any workarounds on how to integrate Erdas's signature editor with the analytical functionality of R would be very helpful. For example, I need the functionality of the signature editor Figure1 ![enter image description here](http://i.stack.imgur.com/1VeY8.png) Yet I need to convert these values to a format that I can read into R. Note that I need all four bands of the input image along with a class ID (i.e. grid code). This table was generated in ArcGIS. Figure 2 ![enter image description here](http://i.stack.imgur.com/wTtek.png)


OpenLayers - Problems With WPSClient

I am trying to incorporate a WPS running in GeoServer into an OpenLayers client. I am using this simple example as my reference. Using either a local install of OpenLayers.js (v 2.12), or the http://openlayers.org/api/OpenLayers.js, I get an error when the WPSClient is created. In Firefox, I get OpenLayers.WPSClient is not a constructor. Can anyone help me out with this error?


arcpy move attributes from one column to another on a joined layer

Sort of new to Python. In Arcmap you can calculate a field on a joined table to just move the attributes from one column to another with just an = sign. Since my layer has feature linked anno on 7 different layers the calculation has to be done in an edit session. I tried this and it ran over three hours and not sure if it was going to finish. This won't work since I have to do this on about 30 fields. I am trying to do this in arcpy. I was using arcpy.CalculateField_management. I don't think that works based on the error I received. The error is invalid property. I think this means that it wants to make a calculation not just move attributes. Here is my sample code:               import arcpy     from arcpy import env          anno_conn = "\\\\any database"     env.workspace = anno_conn          #Variables     Table_A = "data_1"     Table_B = "data_r3"     Layer_A = Table_A + "_lyr"     Layer_B = Table_B + "_lyr"          try:         #join tables         arcpy.MakeFeatureLayer_management(Table_A, Layer_A)         arcpy.MakeFeatureLayer_management(Table_B, Layer_B)              arcpy.AddJoin_management(Layer_A, "visi_link", Layer_B, "prior_visi")              # Using an edit session - feature linked layer can't be edited outside an edit session               # inside try edits will be saved if successful, will roll back if an exception is thrown         with arcpy.da.Editor(env.workspace) as edit:             arcpy.CalculateField_management(Layer_A, Table_A + ".link_id", "!" + Table_B + ".link_id!", "PYTHON")          except Exception, err:         traceback.print_exc()      It seems to work up until the Calc field line. I've searched for answers to this for a day. Any help is greatly appreciated.


Synergis Weboffice vs Geocortex Essentials

I have one simple question what are the main differences between Synergis Weboffice and Geocortex Essentials? I know they are both used for creating web-based mapping application for end users and they depend on ArcGIS server for working but are there any differences between the tools and features they provide?


Using customized Coordinate System for Archaeological site data

I don't know so much about coordinate systems... In my office we use to deal with spatial data coming from archaeological sites. Each site has its own x-y-z coordinate system (GCS). **Three simple ortogonal cartesian axis**. In the last years we have been managing this spatial data through GIS software (ArcGIS), without using specific coordinate system (just leave it as "undefined") I'd like to know if there exisits any GCS designed to deal with such datasets using simple cartesian orthogonal axis, without grid distortions of the typical GCS. In addition, I'd like to know if this system is suitable for using it in an online mapping application. By the way, we manage 2D (ArcMap) and 3D (ArcScene) environments and work with "mm" as length base unit. If such a thing doesn't exists, maybe someone knows how to create it.


Arcgis not seeing my column header to join tables

I am trying to join a data table to the attribute table of a shapefile. I have done all the necessary steps to makes sure it works: removed spaces, not use numbers in column headers, remove illegal characters etc. BUT still when I try to join the table, Arcmap is not reading the GEOID field that I am trying to use to join the two tables. I cannot see it listed in the fields to join. What could be the problem?


popup when new feature is created

I am trying to create a code where a feature is created, and then automatically a popup window will open. Is there anyone who can assist me please? Here is a part of the code and I want the popup to open right after that one.               var draw = new OpenLayers.Control.DrawFeature(     wfs_layer, OpenLayers.Handler.Point,         {             title: "Create a feature",             displayClass: "olControlDrawFeaturePolygon",             multi: true         });      Later on I have a code for the popup:               //---POPUP fyrir Attribute---//     var select = new OpenLayers.Control.SelectFeature(wfs_layer,{hover: true, highlightOnly: true});       map.addControl(select);       select.activate();              var selectClick = new OpenLayers.Control.SelectFeature(wfs_layer,{});       map.addControl(selectClick);       selectClick.activate();            function onPopupClose(evt) {       select.unselectAll();     }          wfs_layer.events.on({     featureselected: function(event) {         var feature = event.feature;             feature.popup = new OpenLayers.Popup.FramedCloud                 ("pop",           feature.geometry.getBounds().getCenterLonLat(),           null,                     '<form name="inputform">'+                     '

NAME:	\                     <input name="name" type="text" size="30" maxlength="30" value="'+feature.attributes.nafnfitju+'" >

'+                     '
Heimild:	\                     <input name="source" type="text" size="30" maxlength="30" value="'+feature.attributes.heimild+'" >

'+                     '
Code1:	\                     <input name="addrattara" type="text" size="30" maxlength="30" value="'+feature.attributes.addrattara+'" >

'+                     '
Code2:	\                     <input name="Source1" type="text" size="30" maxlength="30" value="'+feature.attributes.source1+'" >

'+                     '
Open/closed:	\                     <input name="open_closed" type="text" size="30" maxlength="30" value="'+feature.attributes.open_closed+'" >

</form>'+                     '
<input type="submit" value="Save" onClick="javascript:ourfunction(\''+feature.fid+'\')">
'+                     '
<input type="submit" value="Delete" onClick="javascript:delete(\''+feature.fid+'\')">
',           null,           true,                     onPopupClose         );         map.addPopup(feature.popup);       },          //destroy popup when feature is no longer selected. Prevents showing 2 Popups at the same time     featureunselected: function(event) {     var feature = event.feature;             map.removePopup(feature.popup);       feature.popup.destroy();       feature.popup = null;     }     });      Hopefully someone can assist me, thanks in advance. Sjonni, Iceland


Is there any way to add an area column to new shapefile created from raster with gdal's polygonize in the same step?

I've been going through a fairly lengthy process of polygonizing very intricate raster images, then creating a new layer with the extra area attribute, so that I can then remove the millions of tiny dots that turn what should be a 5 megabyte shapefile into a 40 megabyte one. When I run the Raster|Conversion|Polygonize with qgis, it shows the gdal script that is run in the command line to complete the task. I don't know many gdal commands other than ogr2ogr stuff, so was wondering if any of you gdal pro's out there knew how to add an area attribute to the outputted shapefile (populated with calculated areas) without a second step by adding something to this script:               gdal_polygonize.py "/Users/douglas/Desktop/polygonized_rasters/zones.tif"      I've clicked a couple of options in this polygonize menu and saw code added to the script, so that gave me the idea that there might be some code I could add to this line to add that area column during the same process. Maybe not, I'm not sure. Thanks though anyway


How to export xy-coordinates from a line segment

I need an export of xy-coordinates of line segments (Quantum GIS). I found an answer to a similar question here but that concerned only xy of the first and last vertex of a line segment. But I need the coordinates of all the vertexes. Can someone help?


MapBasic For Loop ignoring every other row

I have made a loop in MapBasic but is it only updating every other row in the table: ![enter image description here](http://i.stack.imgur.com/unw3Z.png) My code works from the bottom of the table upwards to generate start and end coordinates for polyline features:                   Dim nrows as integer          Dim USX, USY, UEX, UEY as float          For nrows = TableInfo(WorTab, TAB_INFO_NROWS) to 1 Step - 1          Select * from WorTab where rowid = nrows         Dim ObjTYPE as integer         ObjType = ObjectInfo(Selection.obj, OBJ_INFO_TYPE)         If ObjType = 4 then         nPoints = ObjectInfo(Selection.obj, OBJ_INFO_NPNTS)         USX = ObjectNodeX(Selection.Obj,1,1)          USY = ObjectNodeY(Selection.Obj,1,1)         UEX = ObjectNodeX(Selection.Obj,1,nPoints)         UEY = ObjectNodeY(Selection.Obj,1,nPoints)         Update Selection Set SX = USX, SY = USY, EX = UEX, EY = UEY         nrows = nrows -1         Close Table Selection         Else         End if      I'd be very grateful if somebody could point out what I'm doing wrong and why nrows -1 is the wrong choice.


import gpx file to shape file

i am using arcobject .net10.1,arcmap10.1,arcgis10.1. i want a shape file which is created from .gpx file(GPS file) .i have search every where but cound'nt find any information related to my query. i am creating some code as per my knowlegde that code is not working properly. ESRI.ArcGIS.Geoprocessor.Geoprocessor gpq = new ESRI.ArcGIS.Geoprocessor.Geoprocessor(); gpq.OverwriteOutput = true;                       ESRI.ArcGIS.ConversionTools.GPXtoFeatures toolq = new ESRI.ArcGIS.ConversionTools.GPXtoFeatures();             string qfpath = @"F:\Aniket Babar\esri project data\sample waypoint files\fells_loop.gpx";             string qspath = @"C:\Documents and Settings\Administrator.REV3\My Documents\ArcGIS\Default.gdb\fells_loop_GPXtoFeatures1";             toolq.Input_GPX_File = qfpath;             toolq.Output_Feature_class = qspath;                  gpq.Execute(toolq, null);      it gives me error like"Attempting managed execution inside OS Loader lock. Do not attempt to run managed code inside a DllMain or image initialization function since doing so can cause the application to hang." if any one having a solution regarding this please send me this link or any related information.


v.clean in QGIS Python console

I have tried to use v.clean and v.clean.advanced (by GRASS GIS) within a script in QGIS, and was not able. I followed instructions here:, but was not able to do so. I got message: "No module named grass.script". Any hint or reference to any source about how to do this?


ArcGIS raster data - Group statistics for classification

I'm trying to figure out how to classify a group of rasters. These can't be joined as they have overlaps which must be kept. I was wondering whether there's a way of obtaining grouped statistics - eg histograms - which might help me visually classify either groups, or all of the data. Basically, an efficient way to classify them all based on aggregate data. I think I'll probably be able to add symbology to more than one once I have the class limits established. They're all separate ADF files at the moment - I'm wondering whether I should get these into a geodatabase first, somehow.


Stream network doesn't match the reality?

I am trying to create flow accumulation grid in ArcMap using Arc Hydro. I have the DEM (25mX25m) and the stream network shapefile. I have done the following 1)Resampled the DEM to 5m; 2)Burnt the stream network in the DEM using DEM reconditioning; 3)Filled the sinks; 4)Created flow direction grid; and 5)Created flow accumulation grid. The flow accumulation grid I have created was matching well with the reality in most of the catchment. However, in the middle part of the catchment, two streams are connected while they are not connected in reality. In addition, the catchment outlet created by Arc Hydro was not matching with the reality as illustrated below. ![enter image description here](http://i.stack.imgur.com/Ey1HI.jpg) Any recommendations how can I make the flow accumulation grid created by Arc Hydro match the reality would be appreciated? Thanks in advance for your help and time.


GeoServer Web Feature Service (WFS) in Open Layers Problem

I am running GeoServer 2.1.1, and have a datastore from ArcSDE populated with features projected in EPSG:26918. I can successfully add a WMS to an OpenLayers map, but am having issues with WFS. I have the map set to EPSG:900913 in preparation for a mashup with Google base layers. I am under the impression that vector layers are reprojected on the fly, and that absent of any user specification, features will be rendered using the default style. Here is the main chunk of my script:                   map = new OpenLayers.Map('map', {             maxExtent: new OpenLayers.Bounds(                 -128 * 156543.0339,                 -128 * 156543.0339,                 128 * 156543.0339,                 128 * 156543.0339),             maxResolution: 156543.0039,             projection: "EPSG:900913",             units: "m"         });              var bridgeLayer = new OpenLayers.Layer.WMS(             "Bridge",             "http://localhost/geoserver/geodatabase/wms", {             layers: "geodatabase:Bridge",             format: "image/png",             transparent: true         },{             isBaseLayer: true         });              var wfsLayer = new OpenLayers.Layer.Vector(             "Photo Points", {                 strategies: [new OpenLayers.Strategy.Fixed()],                 protocol: new OpenLayers.Protocol.WFS({                     version: "1.0.0",                     url: "http://localhost/geoserver/wfs",                     featureType: "Photo Point",                     srsName:"EPSG:26918",                     featureNS: "http://dev.geodatabase.org",                     geometryName: "SHAPE"                 }),                 visibility: true         });              map.addLayers([bridgeLayer, wfsLayer]);      And so, the WMS displays correctly (and when I add a Google base layer, it lines up correctly). However, the WFS layer fails to appear. According to Firebug, however, my POST request is successful, and I receive an response. Furthermore, the layer object I defined looks to be a fully formed layer, although it does not appear to have any features in its "features" array. Not sure what that means. Screenshots appear below. Something also that I don't understand, when I switch versions to "1.1.0", my POST request fails with a "NoApplicableCode" exception. Additionally, I have tried reprojecting the data in GeoServer to spherical mercator EPSG:900918, and also redefining the srsName parameter, but both experiments have failed. Thanks so much. ![](http://i53.tinypic.com/iwi4ol.jpg) ![](http://i54.tinypic.com/xo45sh.jpg)


Looking for maps with high complexity

For an presentation I need an example of a map which contains very much information. The map should :   * respect basic cartographic principles.   * contain as much quantitive and/or qualitive information as possible without beeing overloaded by too much cartographic details.   * look beautiful as I want to use it in a presentation. Concerning the complexity of the map, I want to use an example by Eduard Imhof who tried to describe the course of a river by words instead of using a map and ended up with a wall of text, showing the data efficiency of the map. I want to use this to show how efficient maps are concerning the volume of data which can be stored on it (compared to the same amount of data stored in a database for example). Any ideas? :)


3D subsurface modelling

I am a student and initiating a thesis project for MS specialize degree in geophysics. I am working on tectono - stratigraphy of some area and want to do the 3D subsurface modelling of the same area but have lesser idea, besides am unable to get some proper guide on this too but since i love doing this I want help plzz... As far as i have studied 3D analyst is not very much suitable for this job yet "may be" used in thesis.


Can I view ArcScene 3D models in a free viewer such as ArcReader or another 3rd party viewer?

I would like to build a complex 3D model in ArcScene, export it and send it to someone else that does not have ArcGIS with 3D Analyst. The idea is to interactively navigate the model during a presentation. I've never used ArcReader or ArcGIS publisher before so I have no idea if a complex 3D model can be published and exported. If not what are my alternatives for viewing the completed model outside if ArcGIS?


How do I include libfcgi files when installing QGIS on SUSE Linux?

I'm using SUSE Linux Enterprise Desktop 11 SP2. When I try to install QGIS from YAST, it seems that the package "libfcgi.so.0" and "libfcgi.so.1" are missing. Is there anyone else that has experienced the same problem? Where can I get those files? Thank you very much!


How to label pie charts in ArcGIS?

How can I add to pie charts in ArcGIS (10.1) labels of the values the charts are created with?


Expand raster extent and set expanded values using nearest neighbor sampling

I have run a hydrodynamic model which has given me a grid of surface elevations throughout a bay. I am trying to process the results in ArcMap (or QGIS if needed). My goal is to take the current grid and buffer/expand (not sure of the correct terminology) the extent of the grid by 2 cells the whole way around, and have the new buffered cells inherit the value (water surface elevation value) from the neighboring cells from which they were expanded. From what I have worked out so far there is an Expand tool in the Spatial Analyst toolbox, however, it only works on integer rasters, and I am dealing with a float raster. I converted my grid to an esri raster grid (displayed in red below) and then used the expand tool to expand the raster by 2 cells (expanded area displayed in green below). I then used the raster calculator and a Con statement to recopy the float values back to the expanded raster. However, this does not copy any values over to the expanded cells. So my question is, is there any way to somehow interpolate or use a nearest neighbor type tool to set the values of the expanded cells to the nearby cells from which they were expanded? I have out together an example image to help me explain my problem: ![Example Image](http://i.stack.imgur.com/nTnoB.jpg) Any ideas/help would be great! Thanks,


Cannot Edit Features of PostGIS Layers in OpenGeo GeoExplorer

I admit I originally posted the question here: http://getsatisfaction.com/opengeo/topics/cannot_edit_features_of_postgis_layers But due to no response and not finding a solution, I thought I'd try to ask here. I have imported some layers from my PostGIS DB into the OpenGeo Suite Geoserver. For the most part they work fine, yet I am unable to edit their features. Every time I try to do this I get the following in FireBug: feature is undefined layer.drawFeature(feature,feature.styl...stFeature=layer.selectedFeatures[0];} OpenLayers.js (line 815) However, if I import a layer that came bundled with the suite such as World - Borders, everything works fine and I am able to do feature editing as expected. I have also seen this thread: http://getsatisfaction.com/opengeo/to... But still I am not sure how to go about it. It is indeed true that we had PostgreSQL and PostGIS installed already, and when I installed the OpenGeo suite all the data dumped itself onto my already existing PostGIS (at least I think so), as everything seems to be running on the same port 5432. What can I do to have my PostGIS layers allow feature editing? This was the main reason we decided to go with GeoExplorer ;). Geoserver error logs report:               INFO [geoserver.wfs] -     Request: getFeature     handle = null     service = WFS     version = 1.1.0     baseUrl = http://10.0.104.14:8081/geoserver/     providedVersion = null     extendedProperties = {}     query = [net.opengis.wfs.impl.QueryTypeImpl@431e34c7 (group: null, propertyName: null, function: null, filter: [ tib_villages.fid-5a164ff3_134c8d54c27_-6fb2 ], sortBy: null, featureVersion: null, handle: null, srsName: EPSG:900913, typeName: [{http://opengeo.org/#tibet}tib_villages])]     maxFeatures = 20     outputFormat = text/xml; subtype=gml/3.1.1     resultType = results     traverseXlinkDepth = null     traverseXlinkExpiry = null     formatOptions = {}     metadata = null       Thanks for your time and I am looking forward to hearing any feedback, elshae


How exactly is an SLD_body object (when extracted from its XML href) structured when using a cElementTree to parse it?

I'm trying to export data from a PostGIS database into a geoserver and I've gotten it to where I can get the PostGIS database to point to the geoserver (I'm doing this as a python script for a GIS program) with gsconfig. So essentially I need to establish the services so that I can see the data in the geoserver data store/namespace that I have created. Therefore I'm using the xml.etree from cElementTree package in Python to parse the data from my database using gsconfig's catalog object. However, I don't know that much about XML and how its structured, but I have the gist of the element tree structure.


How to delete multiple feature classes simultaneously?

I created a number of feature datasets to temporarily store the output of some ModelBuilder models. Each dataset contains several thousand feature classes, and the idea was that I could just delete the feature dataset when I was finished with the ModelBuilder output, instead of deleting each feature class individually (which would take forever!). However, I've been using the Delete (Data Management) tool to delete the feature datasets, and have found that all the feature classes within the dataset, rather than being deleted, are just turning up in the folder the feature dataset was stored in. I assumed that deleting the dataset would remove all that dataset's content also, but this doesn't appear to be the case. I Have tried to delete the feature classes using Windows Explorer (where I could highlight them all simultaneously and then press delete) but they don't show up. So, could someone please tell me if there is a way to delete multiple feature classes (or all content within a feature dataset) in one go? Either in ArcGIS (using v10) or through Windows Explorer?


How can I, in a general way, rescale a raster in ArcGIS Model Builder

What I want to do is rescale a few rasters from whatever they are now (i.e., values between -99 and +1000 OR -256 - +256 OR 0 - 3.23) to a normalised value (i.e. values between 0 and 1). I would prefer to do it in model builder with standard tools, but I am not averse to writing a new one in Python. Any info appreciated.


Why is node tool deleting unselected points in lines?

I am freehand editing lines to a shapefile on Mac QGIS 2.2.0, and have my snap settings at vertex and segment @ 20 pixels. The lines are roadways and I have merged them all into one feature and merged the attributes too. When attempting to delete a handful of points together (by clicking and dragging), I noticed that points in lines way up north somewhere are also deleted. So I tried exiting, re-starting, and zooming out to see if it happened again. The extra nodes that are deleted are not selected. Also, when trying to delete the target points one by one, _different_ points are deleted up north too, also one by one. It is as if they were somehow linked. What could account for this?


Import GRASS and SAGA modules to python script

I'm new at this. I am "building" a python plugin for qgis and I want to import modules from GRASS and SAGA but I don't understand how to do it. The plugin I want to import is Kriging interpolation from SAGA or v.surf.idw from GRASS or interpolation plugin from Raster qgis. I want an interpolation method to convert points in grid. Can you help me? * * * I have another doubt. I run these lines in my code:               vectorlayer_vector =  ftools_utils.getVectorLayerByName(inputLayer)     Elevation = self.lineAttrib.currentText()     cellSize = int(self.linePix.value())     outPath = self.inputLayerCombo3.text()     output = ftools_utils.getRasterLayerByName(outPath)           Processing.runAlgorithm("saga:shapestogrid", vectorlayer_vector, Elevation, 0, 0, 4, None, cellSize, output)      When I run, I get an tiff image but with nothing. The .aux file is not created and when I try to open the tiff image that I created appears this message: c:...teste3.tif is not a supported raster data source So, I guess that the output is wrong. What I have to do to get a correct output? What I missing?


Is it possible to move a visualization to another table? e.g. when merging two tables, it creates a 3rd. OR to completely replace a table and maintain the viz?? I don't want to lose the original table's visualization nor recreate it (i've already done it once this morning when moving between accounts). But I've added 5 columns to integrate into the viz. Anything? or should i just recreate it? thanks, chrism


What is the relationship between the SOC and the DBMS?

I have read lots of documentation about ArcSDE 9.3 recently but there are still are few things I don´t quite understand because I have not installed the infrastructure myself. If the SOM manages the connections to services and the SOCs host the services, what is the relationship between the SOC and the DBMS? Is this Workflow correct?   1. Client application requests features from AGS   2. AGS communicates to the SOM that a service has been requested   3. SOM finds the SOC which is runs the service   4. SOC pulls data out of the DBMS and relays the data to SOM and AGS   5. AGS provides visualization and analysis of data


regionating vector data for use in Google Earth

Does anybody know of free/open source tools to regionate vector (kml) data, like the regionate tool in Google Earth Pro (see http://support.google.com/earth/bin/answer.py?hl=en&answer=175007)? I have only found one, which is Regionator (http://code.google.com/p/regionator/). It is a set of Python scripts and I am not sure how to get it to work so perhaps if anybody has a link to tutorials or the like?.


GRASS algorithm r.los does not produce raster

After running the r.los algorithm through QGIS 2.0, I'm getting a 'Problem loading output layers' box which gives the following message: Oooops! The following output layers could not be open Output raster layer: C:/Users/ffinnm/Documents/CGE/LVIA/Southdown/gis files/viewsheds/hhh.grd The above files could not be opened, which probably indicates that they were not correctly produced by the executed algorithm Checking the log information might help you see why those layers were not created as expected This algorithm requires GRASS to be run. A test to check if GRASS is correctly installed and configured in your system has been performed, with the following result: GRASS seems to be correctly installed and configured I have checked the log information and there is no error logged here. Interestingly, in the log information, the algorithms section states that the 'r.voronoi' algorithm ran 2 seconds after the 'r.los' algorithm. I didn't instruct the r.voronoi algorithm to run (or know what it is!). I have also tried changing both input and output raster formats. Any help here would be much appreciated, thanks


Why does ArcGIS ModelBuilder tool Iterate less than expected Number of Times?

I have an ArcGIS ModelBuilder tool (with an iterator in the workflow) and the tool iterates through approximately 30 features then stops. I have checked that:   1. It is not instructed to iterate only a certain number of times; and   2. That the input feature class does not have a selection on it.  Has anyone else experienced this behavior? Additional details : The iterator churns through a dataset and creates a new FC by unique ID. The next steps are complex and include building terrains. The behavior is usually corrected when I shut down ArcMap and open up a new MXD.


Does OpenLayers mergenewparams repeatedly hit a WMS server when you change the layer name?

Does mergenewparams repeatedly hit a WMS server if you change the layer name? For example, you call the server with layer 'a' then you use mergewnewparams and call layer 'b', then you use mergewnewparams and call layer 'a' again. Does it cache layer 'a' or request a new layer 'a' ? I'm actually using a tilecache, but I assume it works the same for either. Thanks!


Merge an entire Geodatabase in ESRI ModelBuilder

How can I merge an entire Geodatabase in the ESRI (ArcView) ModelBuilder? I have several FeatureClasses created into the GDB and in the model I plan to merge them all. Right now I have imported all the FCs by hand into the Model. This is not the solution I am aiming for because the names and the amount of FCs could vary.


UpdateCursor to insert DDM values from a layer in WGS84?

I have a layer in WGS84 (Geographic, DD) that I need to label in DDM (Degrees Decimal Minutes). **I need to calculate within an arcpy UpDateCursor, a label field with DDM in it.** I saw that you can run the Convert Coordinate Notation, but I need to be able to calculate these labels on the fly within my cursor (this script is also moving points around). After trapsing through the Esri forums I have not seen any snippets with how to achieve this. This is likely a noob question, but this is very early steps for me any arcpy.


How to smooth / "reinterpolate" a raster with GDAL?

Short:   I want to change the resolution of a raster and smooth the gray ramp like shown in the images bellow. The preference is to use GDAL, PIL or Numpy. Description:   I'm kriging points into rasters with an output resolution of 20 meters with the High Performance Geostatistical Library. I don't want to change the output resolution because the interpolation time increase exponentially.   With this resolution the output image is ugly (pixelated and aliased). I don't know if it is conceptually correct but I want the image to be smoother like in the example bellow. It's something like 'reinterpolating' the image into a better resolution one. I'm using python so my preferences are GDAL, Python Imaging Library or Numpy. The answer could be theoretical, like pointing out the algorithm name or the concept of this kind of operation. Source:   ![enter image description here](http://i.stack.imgur.com/PP0wg.png) Destination:   ![enter image description here](http://i.stack.imgur.com/T5g3w.png) EDIT Results with gdalwarp cubic spline:   ![enter image description here](http://i.stack.imgur.com/S7ZDn.png)


How can I properly add PostGIS data into MapGuide Maestro?

I am having trouble connecting a PostGIS database to MapGuide Maestro. Actually, I have added the PostGIS data into Maestro, but when I click to validate the layer it shows an error message: `Validating resource: Session:9754d026-d86b-11e1-8000-00270e208411_en_7F0000010AFC0AFB0AFA//259f852d-fc68-4475-9de1-a5e446fb6e82.FeatureSource Warning - OSGeo.MapGuide.ObjectModels.FeatureSource_1_0_0.FeatureSourceType: Spatial context extent appears to be invalid (or is the default extent returned by FDO provider. Meaning it very likely does not truly represent the extent of your data) Information - OSGeo.MapGuide.ObjectModels.FeatureSource_1_0_0.FeatureSourceType: No primary key defined for class: public:geography_columns, features will not be selectable and feature joins on this class may produce unexpected results` Can anyone please suggest where I went wrong? **Is there any documentation for how to add PostGIS data into MapGuide Maestro?** Also, if I click to see the preview of the layer it gives me an error report of `Fusion Error`. I am using PostgreSQL 8.4, PostGIS 1.5, MapGuide 2.2 and MapGuide Maestro 4.01 Thanks in advance..


Overlaping graphics not visible in featureLayer

I have constructed featureCollection from Json and feature collection has 44 features. I have given this feature collection as input to the Featurelayer which is showing only 32 graphics objects of the featurecollection.I found out that overlapping features are not rendered at all in the feature layer. How could I resolve this issue and what is the possible solution to this problem.? Thank you in advance. Thank you in advance.


What are the equivalents of ESRI products in open source for Web Mapping?

I have been using the following licensed products to publish my GIS data on web.   1. web server: windows 2008 R2,   2. GIS desktop: ArcGIS desktop,   3. GIS web server: ArcGIS server,   4. API: ArcGIS API for JS,   5. Database: Oracle with ArcSDE,   6. IDE: Eclipse. But now I want to switch altogether to open source world mainly because of freedom. Fortunately there is plethora information about open source on web but unfortunately are confusing me a lot. There are soo many products like geoserver, mapserver, openlayer, geotool, vaadin, geomajas, uDig, GRASS etc etc. I am not able to related these products with the one I am already using. Can someone please relate the equivalents?


InfoWindow vs. Popup

Can somebody provide a good explanation/comparison between Info Window vs. Popup Window. It seems both can be used with Feature layers. What are the pros and cons of either and when is more appropriate to use one than the other?


How to setup GeoWebCache to produce same grid as GoogleMaps?

Is it possible to create gridSet in GeoWebCache which has same number of zoomlevels with scaleDenominators compatible with google maps? There is EPSG:4326 default, but no mention of number of zoomlevels/scaleDenominators in documentation.


Use of GIS in fixed telecommunication network management

I want to get an idea how GIS can be used in fixed line telecommunication management? there are open source GIS tools as well as Proprietary tools as well. what is the best tool can be used in above ? thanks. Chandima


Activity Dashboard for LocalHost

ESRI recently put out a new tool called the 'Activity Dashboard' that checks for new maps created, map uses, etc. Sadly, this tool only works for those maps created through the ArcGIS online system. One of the problems I have as the new 'Web Map' guy for my organization is that multiple maps have been created, all of them as (basically) one-offs to fix a specific problem, and no one has a master list of them or the services they are running on. This left me with a very annoying situation. I 'thought' I knew about them all, but I find that I am regularly learning about old maps as they break. The recent address locator that ESRI killed off brought SEVERAL maps to my attention that I wasn't even aware of. I am wondering if there are any tools that have been developed that track maps and/or services that have been created to give someone in my position a better chance of managing this system.


Gaussian Plume Model

I'm looking for code to implement a gaussian plume model in 3D. Any suggestions on where to start? Are there any python scripts out there? Thanks


applying SLD to WMS raster in OpenLayers

I am attempting to adjust the symbology of a loaded WMS raster layer (img file through MapServer). It is the NLCD dataset. The WMS layer loads find with               nlcd = new OpenLayers.Layer.WMS(         'nlcd',         'http://server /fswms/html/rlayers',         {             layers: 'nlcd-2006', transparent: true,             projection: new OpenLayers.Projection("EPSG:4326")         },         {             isBaseLayer: false         }     );      And I have a SLD file that I have created in which I am attempting to make all of the pixels with the value of 11 opacity of 0.               <?xml version="1.0" encoding="ISO-8859-1"?>     <StyledLayerDescriptor version="1.0.0"         xsi:schemaLocation="http://www.opengis.net/sld StyledLayerDescriptor.xsd"         xmlns="http://www.opengis.net/sld"         xmlns:ogc="http://www.opengis.net/ogc"         xmlns:xlink="http://www.w3.org/1999/xlink"         xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">     <NamedLayer>         <Name>nlcd</Name>         <UserStyle>           <Title>ncld test</Title>           <FeatureTypeStyle>             <Rule>               <RasterSymbolizer>                 <ColorMap>                   <ColorMapEntry opacity="0" color="#000000" quantity="11" />                 </ColorMap>               </RasterSymbolizer>             </Rule>           </FeatureTypeStyle>         </UserStyle>       </NamedLayer>     </StyledLayerDescriptor>      The WML layer ends up being map.layers[3] in my test application so I attempting to apply the SLD by the following code but it is not working.               //this is in the init:        OpenLayers.Request.GET({         url: "sldcode.sld",         success: complete     });            //function that attempts to apply SLD     function complete(req) {         mySLD = format.read(req.responseXML || req.responseText);         map.layers[3].mergeNewParams({sld_body:mySLD});     }         Here a link to the dev site: http://rain.nemac.org/~derek/ol/sldTest.html Any help would be greatly appreciated. Thanks, Derek


Visualizing KMZ network link in QGIS

Is it possible to open and visualize and kmz network link in QGIS. If so what are the setps involved? Thanks in advance Cheers


GIS mapping software for visualizing and publishing maps

Currently I am writing my Bachelor thesis and for testing purposes I need software where I can visualize and publish maps and put vector data (point, line, polygon) on them. Can anybody recommend such software? I have already founded GIS Cloud. I need something like this.


Where can DEM data for all of British Columbia be found?

Is there a website or server from which a DEM covering all of BC can be downloaded - preferably at 1:50K scale? GeoBase has tiles for 1:50K and 1:250K, but they do not provide it as a mosaic. It's not very time efficient to download 1000+ tiles and mosaic them.


How to get an extent for Raster using python

I have to get an extent of a given raster so that I can save it and use the same to clip another raster with that extracted extent. As of now I am stuck in the extent part.               import arcpy     from arcpy import env               desc= arcpy.env.extent("D:\GIS @ UTD\Sem 1\GISC 6317\Lab\Lab 10\temp\Tahoe\Tahoe\Emer\erelev.grid")          print desc      The following gives an output 0 0 0 0 NaN NaN NaN NaN. Any suggestions.


How to improve cartographic visualization of overlapping polygons?

I need some advice about how to improve the visualization of cartographic information. User can select different species and the webmapping app shows its geographical distribution (polygonal degree cells), each specie with a range of color (e.g darker orange color where we find more info, lighter orange where less info). The problem is when more than one specie overlaps. What I am currently doing is just to calculate the additive color mix of two colors using http://www.xarg.org/project/jquery-color-plugin-xcolor/ As you can see in this image, the resulting color where two species overlap (mixed blue and yellow) is not intuitive at all. ![enter image description here](http://i.stack.imgur.com/2Zmhu.jpg) Someone has any idea or knows similar tools where to get inspiration? for creating the polygons I use d3.js, so if more complex SVG features have to be created I can give a try. Some ideas I had are... 1) The more data on a polygon, the thicker the border (or each part of the border with its corresponding color) 2) add a label at the center of polygon saying how many species overlap. 3) Divide polygon in different parts, each one with corresponding species color. thanks in advance, Pere


QGIS Advanced String Filter Expression returns too many results

`"GNR"` is a string column and I expect the filter expression `"GNR"='.11'` to only return the feature with the value `'0.11'`. Why are the other features with values of `'.110'`and `'.1100'` included? Is there something wrong with my query? ![string filter expression returns too many results](http://i.stack.imgur.com/enfUf.png)


Conversion of diagonally adjacent pixels into one polygon?

I am trying to vectorize a raster whose digital numbers indicate different classes they belong to. ArcGIS and QGIS yield similar results: separate shapes out of pixels sharing the class that are adjacent only diagonally. In this sense, this question is related with this one and also this one. I tried the suggestions posted in those questions, but none of them worked for me. I can not use the buffering-approach because my polygons refer to networks of areas; the connectivity/topology of those networks must remain intact, and I have nodes very close to each other. Dissolve is not useful either; it apparently does not work on polygons sharing a single vertex. ## ILWIS could solve the issue: I found that the vectorization tool of ILWIS (see explanation here) has an option for scanning for the 8 neighboring pixels during the targeting of each raster pixel to a polygon. ILWIS calls it "Output map 8-connected with smoothing". This might solve the issue by avoiding the polygon breakage during the vectorization when input pixels are placed diagonally. I downloaded and installed the ILWIS version available here. However, after importing the raster I want to vectorize (I used the built-in GDAL option), the drop-down menu to select the input layer in the 'Raster to Polygon' tool does not allow me to select anything and thus I can not use the tool. Does somebody know how to use this ILWIS' tool, or how to solve this problem?


QGIS snapping not working for me

I am new to QGIS, and this is a very basic question. I made a new polygon layer, to trace cadastral boundaries off a raster map. I want the polygon nodes on adjoining acreages to snap together. Some do because I have 'avoid intersections' on. But if I have to manually move a node on one polygon onto a node on the adjoining polygon, I cannot get it to snap to it. I've read the manual, tried to sort it out, but I cannot make it work. How do I make it work? Refer to the screenshot. ![Cadastral boundaries won't snap together](http://i.stack.imgur.com/wgoA2.jpg) I have used MapInfo for 14 years. It does not have an inbuilt way to avoid polygon overlaps, but its snapping is very straightforward.


Why I see different data from PostGIS and ArcSDE

I have installed ArcSDE 10 on PostgreSQL 8.4.3 with PostGIS 1.4.2. I did some tests.   1. I imported one feature layer with 8500 features into ArcSDE through ArcCatalog. and I deleted 8499 features in ArcMap and just kept one. I saved the edits. I checked the data using PostgresSQL client tool, like pgAdmin III. I still see 8500 rows in the table.   2. I imported one feature layer with 3 point features into ArcSDE through ArcCatalog. and I add 2 new points in ArcMap. I saved the edits. So I can see 5 points on this layer. But when checked the data pgAdmin III. I still only 3 rows in the table. I closed and opened ArcMap and pgAdmin several times. and still see they have different data. Is anything I did wrong?


Error when trying to import a raster file into a geodatabase

I am using ArcGIS 9.3.1. For some time I have had no problem importing a raster file into my geodatabase. That's recently changed and I have absolutely no idea why. Every time I try to import a raster file, I receive an error message that says: "Problem generating Output Raster Name (ERROR 999999: Error executing function.)" I am at a loss...any help would be greatly appreciated.


Openlayers error after removing selectFeature control

I am using a selectFeature control for a point layer in OpenLayers 2.12. I have a problem after removing the selectFeature control from the map. The reason why I want to remove the control is because the target layer of the control is destroyed and a new one is added every time the map is moved. In this process I also remove the selectFeature control from the map and add a new one for the updated layer. But that does not work properly. After the first move/zoom event of the map get the following error: _Error: getFeatureFromEvent called on layer with no renderer. This usually means you destroyed a layer, but not some handler which is associated with it._ Here is a version of the function `refresh`, which is called on any pan/zoom event of the map:               var mylayer;     var mycontrol;     refresh = function(){               // Remove previous layers and controls         if(mylayer != null){             mylayer.destroy();             map.removeControl(mycontrol);         }              // Get mylayer         mylayer = new OpenLayers.Layer.Vector("GeoJSON",              ...         );                    // Add layer and selectFeature control         map.addLayer(mylayer);            mycontrol = new OpenLayers.Control.SelectFeature(mylayer,{               onSelect: function(e) {             $("#tmpOutput1").html(e.attributes.name + "\n" + e.attributes.id);             }         });         map.addControl(mycontrol);           mycontrol.activate();             }      **UPDATE** To solve the problem I had to change the layer/control removal to this:               if(mylayer != null){         map.removeControl(mycontrol);         mycontrol.destroy();              map.removeLayer(mylayer);         mylayer.destroy();     }      Note that here the order matters (at least in my example): first remove the control, then remove the layer.


How to retrieve transport infrastructure data from google map?

I imported into postgis osm data using osm2pgsql. My objective was to get a database with roads elements for a specific region of Italy (i.e. traffic lights, driving directions and turning direction at junction). However there are a lot of missing data, because almost all the cells are empty. I would like to know if I did something wrong and if there is an alternative way to obtain these data. For instance I would like to know if there exists a way to retrieve data from google map (transport infrastructure or even real time traffic information)? Thank you Selini


How is a geographic heat map different from a colored statistical thematic map?

From **Wikipedia** , > A heat map is a graphical representation of data where the individual values > contained in a matrix are represented as colors. Apparently, a _geographic_ heat map is only one kind, and this is the only kind relevant to me. Cartography is about geographic maps. In cartography, I often distinguish between topographic and thematic maps. From **Wikipedia** , > A thematic map is a type of map or chart especially designed to show a > particular theme connected with a specific geographic area. In thematic cartography, we can at least distinguish between types of thematic maps: statistical, color, choropleth, isopleth, dasymetric, point-symbol, and so on. **How is a geographic heat map any different from a colored statistical thematic map?**


Multiple geometry columns for points, lines, and polygons in spatialite?

I'm rather new to the GIS world and even newer to Spatialite and trying to create a new DB for our application. I would like to create spatial index(es) on the geometry column(s), but our geometry BLOB's have always existed in a single column as WKT. After doing some research, it looks like you need a separate column for points, lines, polygons, multipoint, etc. I'm assuming this from Recipe #6 in the Spatialite Cookbook which states that you need to create your table first, then use AddGeometryColumn to actually add the column which will contain the geometry. In creating this column, one of the items you specific is the GeometryType. I would much rather keep a single geometry column that could handle all of the shape types, but if the only way to create spatial indexes is to have separate columns, I'll go that route. Any advice on this would be greatly appreciated. Editing this to add more info. Folks before me on this project created their own RTree tables and index to pull data that was within an MBR. My goal is to get rid of the home-grown stuff and use Spatialite for our queries. We have a number of shape files that contain points, lines, multilines, and polygons and I would rather insert their BLOB (WKT) into a single geometry column vs. 1 for each type. Our main queries are to find all shapes that are within an MBR. I hope this helps explain a bit better the end goal. Thanks.


Extend Polyline in ArcGIS without intersecting line

I am trying to extend the ends of each line I have in a seismic survey 2 miles. I have tried using the extend line (Editor) tool but it needs an intersecting line to move to. Does anyone know how to extend a previously created line a specified distance on both ends without intersecting another line?


How to create Race Route with arrow directions and distance markers using ArcGIS for Desktop?

Anyone have any suggestions regarding how I would approach creating a race route with arrow directions and distance markers in ArcGIS? I could do it manually but I thought I recalled seeing something about a plug in that aided in this type of work.


Portal Sharing API - "created" search parameter

I was wondering if anyone had an example of a working REST string searching using the "created" parameter. I have read **the docs** three or four times and every time I try to use it, it returns 0 results. Even though I know there are items in there that were created in the time frame I am providing. Example:               http://<srv address>/sharing/search?f=pjson&q=created:[0000001104537600000%20TO%200000001388534400000]      That is searching everything from January 1, 2005 January 1, 2014. The output from the above:               {       "query": "created:[0000001104537600000 TO 0000001388534400000]",       "total": 0,       "start": 1,       "num": 10,       "nextStart": -1,       "results": []     }      What am I doing wrong?


How to make lines styles with offset join properly?

I want to draw the lanes of a two-lane street in different colors. I tried to apply an offset to a symbol layer which works fine for a single line but the ends of two connected lines do not join properly. Depending on the angle, the lines overlap or there is a gap between the lines. Is there any way to solve my use case?


PostGIS - improving sql statement

I've managed to create the below statement that creates a table based on a point in polygon query. However it creates a column that gives a true/false value which I then delete out all the f values from. Is there a better way of producing this statement?               CREATE TABLE crime_test AS      SELECT a.*, b.town, ST_within(a.the_geom, b.the_geom) "within"     FROM crime_data a, nh_areas b;     DELETE FROM crime_test WHERE within = 'f'


Qgis - Labels automatically resize to fit within the polygon

The question of my problem is simple, I want to make labels automatically reduce to fit within the polignos. Type the biggest size is 12 and goat decrease the size to be restricted to the limits of Poligno. Exists based on the calculator using polignos option field area, but this is time consuming due to the amount of resources. Can I use plugins labels manualemnte easy, but that would take too long and would be basically a CAD but without assistance. ![enter image description here](http://i.stack.imgur.com/Xroab.png)


encoding problem with QGIS 1.8 and CSV import

I have following problem with QGIS 1.8. Lisboa (standalone installer, Win 7 Pro 64bit) - I import CSV - text in CP1250 separated with ";" and using "," as decimal separator. I correctly set the separated text importer, get correct point layer which has all text entries displayed correctly. I save this as shapefile, set the CP1250 encoding but when I load this layer in QGIS (choose CP1250 again during SHP loading) but get point layer with damaged text entries. I use this tool for a long time and did not have such problems. Did you have encounter such problems? thanks


Couldn't load SIP module (QGIS doesn't load Python)

I recently installed QGIS 2.4, but Python does not load. I am running on a Windows 64 bit machine. I receive the error message "Couldn't load SIP module Python support will be disabled". Traceback (most recent call last): File "", line 1, in ImportError: DLL load failed: %1 is not a valid Win32 application. I have another installation of Python at C:/Python27. It is the 32-bit version. When I remove C:/Python27 folder, QGIS works fine, so there is some linkage between the new QGIS installation and the C:/Python27 python interpreter. I can't permanently delete C:/Python because I need it for other applications, so I am wondering if there are any suggestion as to how to fix this. Thanks for any advice.


QGIS error when add WFS

I'm running QGIS 2.4.0-Chugiak on win7 64bit and am getting the following error while trying to add a wfs layer (using the simple Add WFS layer command): The layer http://copadesvr02.copade.neuquen.gov.ar/geoserver/wfs?request=getCapabilities&SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=cite:Comunidades_Indigenas&SRSNAME=EPSG:22182 is not a valid layer and can not be added to the map I managed to add some layers just a while ago but now can't even add those same layers, and I just tried the same WFS service on gvSIG and was successful. Can anyone help me out? The WFS url is: http://copadesvr02.copade.neuquen.gov.ar/geoserver/wfs?request=getCapabilities


Recently I've been designing maps without navigation tools like the Zoom In/Out and Pan tools, or arrows to pan north, pan west, etc. Instead, most modern mapping APIs support the direct use of the map itself to handle navigation - drag the map to pan, double-click to zoom in, scroll the mouse wheel, shift-drag to zoom, etc. **Is there a set of standard accessibility guidelines for online maps?** Are there valid reasons for needing a "click to pan" control, or any of the other old-school navigation controls? For example, are they used by screen readers or other accessibility aids? I note that Google Maps has a Pan control, even though it's not strictly necessary. ![enter image description here](http://i.stack.imgur.com/nww0p.png) Thanks.


How to export MXD layout view to GeoTiff?

It my be old issue, but I can't find right solution. How to produce GeoTiff from layout view not from data view? Any simple solution? Thanks


Changing symbology as per scale reference in ArcMap?

I am using ArcMap 10 and I want to display my village parcel as per given scalewise symbology. I have one idea that copy and paste same layer. But my question is that in a single layer, is it possible to give different symbology as per the scale? e.g               Layer Name    symbol   scale     xyz            -----   100000     xyz            *****   200000


Integrating Maximo with Mapinfo?

I am software developer and have been working on database and Java as major platform for development. I am using Mapinfo for GIS analysis. Now, the system needs to show trouble tickets in Mapinfo that are achieved from Maximo software. How is it possible to connect these two application?


How to test if ArcGis server 10.1 is running correctly?

I've recently installed arcgis server 10.1, but I'm not sure if it installed correctly, or if the proper ports are open. Here is my reasoning and I may be off: 1) When publishing a new service(i.e newMap), the folder "arcgisserver\directories\arcgisoutput\newMap" remains empty, I'm assuming their should be a copy of the service in there now. 2) When trying to access the service through arcgis javascript api, I get a blank page. I copied and pasted the source code from the tab or window that pops up when you press "view" in the services thumbnail into a .html file. Then I replaced the url with the REST URL given to me under the services' capabilities options. I've also tried replacing localhost with my machine name and I still get a blank page.               http://localhost/myorg/rest/services/MyMapService/MapServer      If anyone can provide me with some working test code, that would be great, as I'm learning to use javascript api and server at the moment and I want to make sure its configured correctly. Thank you.


MapInfo - programmatically find if layer is hidden

I am working on an integrated mapping project. I would like a button that toggles a layer. I could hold the layers status in a variable and then show/hide it as appropriate but this seems a little hacky. Is there a way to find out whether a layer's "display" value is set to on or off?


SharpMap UK OS data

I've downloaded the UK raster and the corresponding tab file from the ordnance survey website and have displayed the image using the GDARaster class. however I am now stuck as I have no idea how to use the .tab file that came with the map and how I can display points on the map using long and lat. I can't seem to find a working example of how to do this can anyone help? Thanks.


Looking for QGIS introduction course in Los Angeles

To Whom It May Concern: I am in the process of launching a commercial real estate business and I wish to incorporate QGIS as a tool for keeping and analyzing demographic and market research information. I am new to QGIS and GIS in general and I need to find an instructor, seminar or class in the Los Angeles, California area to learn how to use QGIS. Also, I am an Apple Mac user and compatibility with Mac is important to me. Please advise. L. Carlos Simental


Access Trimble GPS on ArcGIS Engine App

Hopefully this isn't too obscure. I have a Trimble Yuma tablet running Windows 7. I need to build an application that allows viewing of points on a map and editing of a geodatabase (feature layers and tables). I'm hoping to use ArcGIS Engine SDK. BUT...how do I leverage the GPS of the Trimble within my Arc Engine app? Anybody done this with a Trimble before? Any code samples out there? Help!!!


Add or Import Metadata to a shapefile

Is there a way to add metadata to a shapefile in Qgis? I have the shape file and an xml metadata file and I want to import in shape file.


How to combine a vector and a raster file?

I want to combine a raster and a vector file. After combining, can I query on the resulting map to obtain a specific area? How?


Using zip file in html

I have my map files vdcar.js (about 27.79 MB) and vdcar.js.zip (about 9.79 MB) in same folder in server. Also have TestPopden.html (http://ngiip.gov.np/THEMATIC_MAPS/TestPopden.html) and .htaccess (empty file) in same location. I found that it takes maximum time (more than 30 minutes) to open the unzipped map file. What I need to do (with TestPopden.html and .htaccess files) to use zipped version of my map file?               <!DOCTYPE html>                 <title>Population Density</title>       <meta charset="utf-8" />              <link rel="stylesheet" href="leaflet.css" />       <link rel="stylesheet" href="POPCEN/Popcen.css" />


Changing from portrait to landscape layout?

I am trying to switch to from a portrait to landscape layout. I've already switched up the buttons in page and print set up and for some reason, the print layout icon shows my page still having portrait features in a non- completely filled landscape icon. I'm a graduate student in urban planning who still has his GIS training wheels on.


How to geocode US addresses using gisgraphy (400k)

Would someone please direct me in geocoding around 400k addresses using gisgraphy or openstreetmap? Thanks in advance, Aji.


Sextante in QGIS - downloaded, but does not show up

I downloaded Sextante, and I put it into the directory C:\users[my_user_directory].qgis\python\plugins\sextante-1.0.9\sextant I re-started the computer, but there is still no sign of Sextante in QGIS, neither in plug-ins manager, nor menu items. Do I have one directory layer too many in the file path above? Is that why it doesn't work? Thanks.


Converting dataset between coordinate systems

I'm working with NOAA's PRECipitation REConstruction over Land dataset, and am looking for an efficient way to sample the data by pixels on NASA's MODIS grid, using open source tools. The PRECL data uses 0.5 degree blocks, so I'll have quite a few repeated values. I'm currently using a mapping generated by ArcGIS that specifies the array index within the PRECL data for every MODIS pixel. This works (and it's fast, as I'm just doing a series of array lookups), but it'd be nice to get away from that extra data set, especially if there were some way to do the transformation in code, without Arc. Thanks in advance, all!


How to know contents of column names in GRASS?

I am using GRASS to extract a specific area from a vector map. Here is the script that I am running: v.extract PHL_map out=makati where="AREA_NAME = 'Makati" The problem is, I don't know what are the values of the columnnames. I can only view the columnnames but I cant view the values under the columnnames


CQl Filter Don't work with WFS layer in Open Layers

I used a WFS layer from geoserver in open Layers it's shown on MAP but when I want to use a filter it doesn't work I don't know where the error this is the code :               <script type="text/javascript">     OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";          function init() {     //set up a save strategy     var saveStrategy = new OpenLayers.Strategy.Save();          //set up projections     // World Geodetic System 1984 projection     var WGS84 = new OpenLayers.Projection("EPSG:4326");        // WGS84 Google Mercator projection     var WGS84_google_mercator = new OpenLayers.Projection("EPSG:900913");             //Initialize the map     //creates a new openlayers map in the <div> html element with id="map"     var map = new OpenLayers.Map ("map", {        controls:[     //allows the user pan ability     new OpenLayers.Control.Navigation(),     //displays the pan/zoom tools                        new OpenLayers.Control.PanZoom(),     //displays a layer switcher     new OpenLayers.Control.LayerSwitcher()        ],     //projection: WGS84_google_mercator,     displayProjection: WGS84     } );          // map extent     var mapextent =     new OpenLayers.Bounds(-7.70912,33.65164,-7.39876,33.52465).transform(WGS84, WGS84_google_mercator);          //base layers     var openstreetmap = new OpenLayers.Layer.OSM();     var google_maps = new OpenLayers.Layer.Google(     "Google Maps", {     numZoomLevels: 20     }     );     var google_satellite = new OpenLayers.Layer.Google(     "Google Satellite", {     type: google.maps.MapTypeId.SATELLITE,     numZoomLevels: 20     }     );          //wfs-t editable overlay     var wfs_layer = new OpenLayers.Layer.Vector("Editable Features", {     strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],     //projection: new OpenLayers.Projection("EPSG:26910"),     protocol: new OpenLayers.Protocol.WFS({     version: "1.1.0",     // loading data through localhost url path     url: "http:// localhost/geoserver/wfs",     featureNS :  "http://www.opengeospatial.net/cite",     maxExtent: mapextent,     // layer name     featureType: "patrimoine",     // geometry column name     geometryName: "geometry",     schema: "http:// localhost/geoserver/wfs/DescribeFeatureType?version=1.1.0&;typename=cite:patrimoine"     })     });          map.addLayers([openstreetmap, google_maps, google_satellite, wfs_layer]);     map.zoomToExtent(mapextent);          }          function filter(){                      var filterType ="cql";                 //var filter = document.getElementById('zone').value;                 // by default, reset all filters                 var filterParams = {cql_filter: null};                 if (OpenLayers.String.trim(filter) != "") {                     if (filterType == "cql")                          filterParams["cql_filter"] = "batiment='Bastion de la sqala'";          // In this line my parametre it' batiment I made a test on the batiment wich has a value "Bastion de la sqala"                                     }                 // merge the new filter definitions                 mergeNewParams(filterParams);                  }                  function mergeNewParams(params){                 Region.mergeNewParams(params);                  }      and whene I type in console map.layers[2].features it show me an error **Cannot read property '2' of undefined** And Thanks for the help


How do I label my map with attributes?

help is always appreciated. I have a state map in QGIS 2.0 with all of the county names and some other info listed in the attributes table. I would like to have the name of the counties and one of the attributes visually displayed on the map in each respective county. This should be relatively simple...but I can't seem to figure out how to do this. Thanks very much.


Not all features available in OSM exports in TileMill

I'm using TileMill and OSM Bright to style an OpenStreetMap for hiking. The problem is that OSM Bright does not import all the OSM features for styling. I believe it's possible to add features in the osm-bright.osm2pgsql.mml file, but I'm not certain how. Does anyone know how to add, for example, path labels? Thank you.


finding the index of a field with its name using arcpy

I wan to extract the index of a field with its name using arcpy. it seems what I have written is so dumb. is there any simpler method with out using for- loop.                   import arcpy              fieldList = arcpy.ListFields(self.address_ref_layer1)         self.field_name = []         for field in fieldList:             self.field_name.append(field.name)         #self.remove_pseudo_nodes()         print self.field_name         print self.field_name.index(u"NAMN1")      it could also be written:                   fieldList = arcpy.ListFields(self.address_ref_layer1         name_index = [f.name for f in fieldList].index(u"NAMN1")      the ListFiled function returns list of field objects with many properties as well as name but no function to return the index of the field which has a specific name.


How do you set sea level elevation to white in terrain analysis?

I have ASTER files I am working with to create a shaded relief background for some other vector shapefiles. In experimenting with different color-elevation options, I am unable to have the Terrain Analysis tool make the surrounding ocean white in color. I have tried manually adding a white color for elevation values 0.0 to 0.0, but no dice. Going to the Table of Contents->Properties and selecting a different color scheme just replaces the Terrain Analysis colors and I have to re-do Terrain Analysis (although that is likely because I don't know how to 'undo' the color change when the 'undo' option is not available). Does anyone have a solution?


Import data, analyze spatially, and export list of points using Google Earth

I am trying to:   1. import data into Google Earth from an outside data source.   2. Then overlay a map and isolate which datapoints fall into that area.   3. I would like to then export a list of the data points to give to a person that works with our company. Can I do this in Google Earth Pro? Can you direct me to a resource that would walk me through this? Is there a better solution for doing this type of work if it can't?


Normalization of Kernel Density Estimate rasters in ArcGIS 10

I've been told by several GIS instructors that the pixel values of a Kernel Density Estimate raster output by the point KDE tool in ArcGIS 10.1 represents the probability of encountering a point, not the density of points at that location. I've been told that the KDE surface is normalized to a probability density surface by dividing the surface by the total number of points used to make it. My problem is that there is no reference to this normalization procedure in the ArcGIS help files at all, with the possible exception to a reference to the paper by Silverman (1986) from which the quadratic kernel function used in ArcGIS 10.1 is derived. Can any one tell me if this is the case? In addition, if this is the case, if you select a field value instead of just using the point location to build the KDE surface, does this normalization procedure still occur, and how would you interpret the surface?


Convert latitude and longitude into northings and eastings

Hi I have the following UK postcode **dy8 3xt** and know that the latitude and longitude is:-               54.452772     -2.156082      I also know that the Eastings, Northings for the postcode is:-               389490     283880      However I am struggling to find the equation that converts lat/long to northings and Eastings, I would prefer to have the equation in both in jScript and c# (I am being greedy)! Can anyone help? **EDIT** Thanks for your help so far guys, I am starting to learn something here esp. the terminology... Some more info, if you click on this link you can see the results I am looking for. The postcode I entered projects to lat/lng using WG S84 and the grid ref projects to OSGB. So my question is how is this done? **WHAT I LEARNT** Thanks to all that answered, I finally got led to here which I can confirm works great