Export SRTM30 Dem file as 3d

i have downloaded some SRTM Dem files here: http://dds.cr.usgs.gov/srtm/version2_1/SRTM30//srtm30 They are from the USGS site. I would like to export them as 3d. I have ran through a few tutorials online, eg 1) http://www.brocku.ca/maplibrary/proc...iew/3dview.htm 2) http://prod.library.utoronto.ca/maplib/3D.pdf But still don't get it. Question: 1) so are the SRTM DEM files exportable to be viewed in 3d via Arcscene/ Arcmap? 2) how do you do it? Thank you!


Need advice on hardware/os upgrade

Forester, not a geopspatial or IT professional, currently using QGIS and SAGA for generally county level geoprocessing of multipspectral imagery and DEM's (cluster analysis, zonal statistics, stream channel modeling). Currently exceeding practical processing limits for laptop/windows configuration using HP EliteBook w/4GB Ram, Intel Core i5, 32 bit Win7 OS. Need advice on logical next step in hardware/OS/software. Future ambitions include deployment of PostGIS and derivation of canopy height models from stereo matched imagery. Thinking low end workstation and Ubuntu OS. Guidance/advice/considerations and recommendations much appreciated. I'm a one man shop trying to operate a bit beyond my capabilities. I do not speak fluent IT or GIS, need a balance of horsepower, idiot proof, moderate learning curve and intermediate support. Cost is important but functionality, flexibility and capability are primary. Thanks.


how to display polygons as points until certain zoom level?

can somebody tell me how to display polygons as points until you zoom to some level. I'm speaking about OpenLayers. So, for example, if the zoom level is less then 10 than it should display feature as a point and when you zoom to 10 or further it should display it as polygon with the shape it really has. I've tried something like this:               var my_style = new OpenLayers.Style({     'fillColor': '#9900cc',     'fillOpacity': .8,     'strokeColor': '#330066',     'strokeWidth': 1,     'pointRadius': 4     });     var vector_style_map = new OpenLayers.StyleMap({     'default': my_style     });          map.events.register("zoomend", map, function () {         if (map.getZoom() < 10) {             console.log("zoom level is less then 10");             layer.styleMap = vector_style_map;         }     });      But it just paint my polygon with defined colors, it doesn't change the shape to point. I'm probably missing something obvious but I would appreciate if you could help me. Thanks


How to update polygon longitude if negative?

I have a polygon, with longitude in [-180, 180], with 2 parts around the dateline. I want to make a map centered on the date line. My plane is to: 1./ update the longitude, so that if lon < 0 , lon=lon+360 2./ dissolve the two part to get a unique polygon. My problem is now 1./ I'm trying:               ogr2ogr nodateline\new_polygons.shp org_polygon.shp -dialect sqlite -sql "SELECT ShiftCoords(geometry, 360, 0) FROM org_polygon WHERE X(geometry)<0"      and the command returns the following error:               Warning 6: Normalized/laundered field name: 'ShiftCoords(geometry, 360, 0)' to 'ShiftCoord'      EDIT: I also tried with `ST_X(geometry)` instead of `X(geometry)`: same error. The error only occur with the "WHERE" clause. a./ what does the error mean? b./ what would be the correct syntax?


Is IGraphicTracker leaking memory?

It seems that when geometries are removed from an IGraphicTracker instance, the geometries remain in memory until IAoInitialize.Shutdown() is called. I have tried removing the geometries by calling the IGraphicTracker.RemoveAll() method and by storing the IDs returned by IGraphicTracker.Add() method and calling the IGraphicTracker.Remove() method for each ID. Is there something that I'm missing in the example provided below? I'm currently using ArcGIS 10.1 and Visual Studio 10 targeting .NET Framework 3.5. An example can be found here: http://forums.arcgis.com/attachment.php?attachmentid=20140&d=1356115382 This problem was originally noted on the ESRI Dev Forum: http://forums.arcgis.com/threads/72516-Memory-Leak-with-Graphic-Tracker


OGR: How to save layer from AttributeFilter to a shape filter?

I select a subset of features from my shape file like so:               node_lyr = nodes.GetLayer()     q_str = "%s = %i" % (field, value)     node_lyr.SetAttributeFilter(q_str)      What is the best way to make a shape file from this subset?


What is a non-geographical feature?

When using the term "feature" in GIS, what is often (if not always) meant is a _geographical_ feature which "is a representation of real world phenomenon associated with a location relative to the Earth" (ISO 19125-2:2004). Are there non-geographical features and if so, are there any examples of such? For example, could a compass or a scale bar be considered non-geographical features?


Conversion from NZMG to latitude and longitude

I want to convert `NEW ZEALAND MAP GRID` coordinates i.e. `Northing` and `Easting` into `WGS84` coordinates i.e. `Latitude` and `Longitude`. I have searched on internet but there is no proper explanation of how to do this, or an online calculator to do this. My final goal is to write a program in `C#` or `JAVA` , which will convert `NZMG` coordinates into `WGS84` coordinates.


What is the most "logical" route generating site (for Europe by car)?

What is the most logical route generating site do you think is out there? By logical I mean optimizing the driving time between two points. (Not always is the fastest way to proceed found along the shortest routes, making use of main-ways and highways when needed.) I tried Google Maps. It seems to make use of highways a lot. Bing Maps, Mapquest, and Map24 all have different calculation methods but they seem to not use highways as much as Google Maps. I mentioned "for Europe by car" because I would like to take a trip from _københavn(denmark)_ to _tekirdağ(turkey)_.


How to copy attributes from one polygon layer to another?

I have a problem that I can't seem to get my head around. I have two polygon layers:   * Polygon A - is a subset of polygon B with the same fields and has identical polygons to Polygon B   * Polygon B - has the attribute data I want to be in Polygon A How can this be done? I tried the QGIS tool "Join Attributes by Location" but as some of the polygons are within others, it tends to link to the first intersect it finds (the outer polygon).


How would I convert this point to WKID 4326?

This is a point in Brewster, MA, USA. How would it be converted to WKID 4326? 324299.3348436486, 832511.5939650714


Looking for advice on getting GIS data into existing AutoCAD site survey

I've got a variety of GIS data from Albemarle County in ShapeFile and DWG format. I'm looking to get the topographic layers and property boundaries together with my existing site survey created by a local surveyor. What's the best practice for doing this? Is it "just XREF it all into one master file" and scale/move around or is there a better way to make this data normalized? I'm in one of the corners of the grid squares that the county has created, so I'd probably need to have 4 of them combined. Edit: I have AutoCAD 2013. I have very little true GIS experience, but I have friends who have offered to help with QGIS.


How can you insert a table record without geometry into a feature?

I am working with ArcGIS and Python and trying to append table rows into a featureclass. You cannot use append and you cannot just make a tableview of the featureclass and then try to append it in. Any ideas? Thanks in advance.


How to draw rotated symbol RASTER type layer in Mapserver?

[PARTIALLY SOLVED]: I was on the wrong track I need to use a UVRASTER, see this part of the docs But doing leaves with with another question about how to get the angle correctly, please see this question I have a raster data source (GRIB2 file) which contains wind direction data for the whole world at an 1 degree interval. I use this mapfile:               MAP         EXTENT      -180 -90 180 90         CONFIG      "CPL_DEBUG" "ON"         CONFIG      "MS_ERRORFILE" "/var/www/mapserver/routeguard/poc.log"         SYMBOLSET   "symbols.sym"         FONTSET     "../fonts/fonts.txt"         IMAGETYPE   PNG8         IMAGECOLOR  255 255 255         SIZE        1600 900          # ------------------------------------------------------------------------------     # the Map Projection     # ------------------------------------------------------------------------------         PROJECTION             #which to use????             "init=epsg:4326" #54004 #4326             "units=m"         END          #----------------------------------------------------------------------------     #grib wave direction layer example     #----------------------------------------------------------------------------         LAYER             NAME "grib_wave_direction"             TYPE RASTER                  DATA "/var/www/mapserver/data/grib/routeguard/4nwcurpol99_2013110412.grb2" #Current direction @ Mean sea level             PROCESSING "BANDS=1.1"                  CLASS                 STYLE                     SYMBOL 'current_S'                     ANGLE [pixel]                     MINSIZE 20                     MAXSIZE 30                     COLOR 0 0 180                 END             END  # CLASS                  PROJECTION                 #"init=epsg:4326"                 "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +pm=-360"             END         END     END #end map file      This generates the following output: ![output from mapserver](http://i.stack.imgur.com/Y0oZj.png) Obviously, the symbol is not drawn, and also not rotated. **How can I draw the symbol and rotate it based on the pixel value?** The symbol is defined in my symbols.sym like so:               SYMBOL         NAME "triangle"         TYPE VECTOR         FILLED TRUE         POINTS             0 1             .5 1             1 1             0 1         END     END      Thank you! I'm on MS 6.4


Ogr2ogr does not support file name with non-ASCII character?

I tried to use the following script to convert shapefile to mif:               ogr2ogr -f "MapInfo File" D:\docs\arcgis_data\shp_mif\BOU.mif D:\docs\arcgis_data\shp\BOU.shp      it works. But I found that if the name of the shapefile has non-ASCII characters like this:               D:\docs\arcgis_data\shp\臺北.shp      Then the covert will fail with message:               Unable to open datasource `D:\docs\arcgis_data\shp\臺北.shp' with the following drivers     ......      Any idea to fix it?


OpenLayers: make polygons not zoom

I have a simple OpenStreetMap for which I want to display a series of circles. Each circle represents readership for that area, so its radius will change based on the numbers. However, when I create the circle and zoom in the map, the circle zooms in as well - it seems like OpenLayers assumes the circle is "glued" to the map and represents a piece of it, when really it's more like a programmatically created marker. How can I get the circle to stay the same pixel radius no matter what the zoom level of the map? Or is there some way to get a marker to be created as a Vector instead of an image?               var map = new OpenLayers.Map({       div: 'geomap',       projection: "EPSG:900913"     });     map.addLayer(new OpenLayers.Layer.OSM());          var epsg4326 = new OpenLayers.Projection("EPSG:4326"); // WGS 1984 projection     var projectTo = map.getProjectionObject(); // The map projection (Spherical Mercator)     // test point     var lonLat = new OpenLayers.LonLat(-74.265982, 40.665443).         transform(epsg4326, projectTo);          var vectorLayer = new OpenLayers.Layer.Vector("Overlay");     var point = new OpenLayers.Geometry.Point(lonLat.lon, lonLat.lat);          var mycircle = OpenLayers.Geometry.Polygon.createRegularPolygon         (             point,             50000,             40,             0         );          var featurecircle = new OpenLayers.Feature.Vector(mycircle);     vectorLayer.addFeatures([featurecircle]);          map.addLayer(vectorLayer);     map.zoomToExtent(vectorLayer.getDataExtent());      UPDATE: Here is the final answer (thanks Rene)!               var map = new OpenLayers.Map({         div: 'geomap',         projection: "EPSG:900913"       });       map.addLayer(new OpenLayers.Layer.OSM());            var epsg4326 = new OpenLayers.Projection("EPSG:4326"); // WGS 1984 projection       var projectTo = map.getProjectionObject(); // The map projection (Spherical Mercator)       var lonLat = new OpenLayers.LonLat(-74.265982, 40.665443).           transform(epsg4326, projectTo);            var vectorLayer = new OpenLayers.Layer.Vector('FSA', {         styleMap: new OpenLayers.StyleMap({           "default": new OpenLayers.Style({             pointRadius: "${radius}",             fillColor: '#ff0000',             fillOpacity: 1.0,             strokeColor: '#000000',             strokeWidth: 2,             strokeOpacity: 1.0           })         })       });       map.addLayer(vectorLayer);            var feature = new OpenLayers.Feature.Vector(           new OpenLayers.Geometry.Point(lonLat.lon, lonLat.lat),           {             radius: 10           }       );       vectorLayer.addFeatures([feature]);            map.zoomToExtent(vectorLayer.getDataExtent());


How to generate geodesic linestring from two points in C#

I'm looking for a **_FREE_** C# class, library or set of functions that will allow me to generate a **_geodesic_** linestring (or array of vertices, whatever) from two decimal degree points. I've found javascript functions that do what I want, but no luck with C# ones. I suppose I could translate it from the above javascript. But I'd prefer to use something tried and tested. Anyone know of such functions, classes or libraries?


Finding municipal-level administrative boundary data for India?

Does anyone know how to find a shapefile with boundaries of Indian municipalities? I have tehsil boundaries as for year 2001, but I would like boundaries of municipalities and, ideally, also historic ones. I believe Census of India has the information but I couldn't find any maps or shapefiles...


Create Friendly Labels for ArcGIS Print Service Options with JS API

How can I change the print options received from an ArcGIS print service to be more friendly and descriptive using the Esri sample print service?


How to create Online IDW(Inverse Distance Weighted) Interpolation?

I have some point data which store as JSON.As follow: `[["Lat":"12.12","Long":"23.23","ZValue":"123"],...]`. I tried to interpolate them as contour line and displayed them on my website ? I do know how to interpolate data in ArcInfo,so I tried to publish a geoprocessing task and using ArcGIS JavaScript API to generate them but it failed to published to ArcServer.![enter image description here](http://i.stack.imgur.com/44kp3.jpg) However,I search for a long time if there exits some other solution but found nothing about how to implement this. Thanks for any suggestion!


How to create a new service area analysis?

I need to create and run a new service area analysis. I'll need to do an analysis of how far you are coming by car from one hospital in 30 min, 60 min and 90 min. Is there anyone who can help me with a model of this in ModelBuilder?


Use Python to parse string

OK, I know that this is Python 101, but I can't seem to figure it out. I have a function setup that will return a string, specifically an azimuth in DMS. I need to take the return values and pull out only the value that appears prior to an ":" So I may get a value x = 123:34:12.4 I need x to then equal 123 and use that to updateRow in an attribute table. As usual... THANKS!


What is the most appropriate Projection for the Gulf coast region?

What are the most appropriate datum and projection I could use for the Gulf coast region - spanning Louisiana, Mississippi, Alabama, and Florida? I have to analyze block group data along the coast. Thanks


How can I edit and save feature attributes of a WMS layer in OpenLayers?

I am creating a data editing application in OpenLayers where the user will need to be able to edit the attributes of features (which are stored on a GeoServer installation). I have seen documentation for editing features in a WFS layer, but I haven't found anything detailing how to edit the attributes of a feature in a WMS layer. Is this feasible to do using a WMS layer? I am currently retrieving data from GeoServer using the GetFeatureInfo request. I would like to load this data into a div, modify it, then save it back to the server. I haven't seen any examples of OpenLayers maps doing this, but I feel like it should involve a simple POST request to the server.


Update GEOS to 3.4

I want to use the rgeos function gDelaunayTriangulation (which requires GEOS 3.4) on my Mac, but when I load rgeos it says:               > library(rgeos)     rgeos version: 0.3-2, (SVN revision 413M)     GEOS runtime version: 3.3.3-CAPI-1.7.4       I have installed and uninstalled both rgeos and GEOS (GEOS framework v3.4.2-2 from www.kyngchaos.com) a number of times and in different orders, but no matter how I try do I get rgeos to say the runtime is anything else than 3.3.3. What should I do to get GEOS 3.4?


How to vary the infoWindow size for different feature layers

I have multiple feature layers on the map, and I wish to display a different sized infoWindow for each layer. As far as I can tell, the map has one infoWindow, which can be resized using: `map.infoWindow.resize(width,height)` Each feature layer then provides the infoTemplate which contains the Title and Content, but no sizing information (right?). How can I show a large infoWindow when clicking on FeatureLayer1, and a small infoWindow when clicking on FeatureLayer2?


Title data from KML/SHP to google maps

Im trying to take land title data from a Shape file and convert it so it will display on Google maps. Works fine going to Google earth I can get the plots and be able to select each one and pull up the relevant data. But I am unable to find anything about doing this on Maps, so far all I have managed is to show the plot data, but cannot call up any data on it. I have been hunting for solutions for a week now and Im at my wits end. Any one have any ideas? Thank you


ST_Geometry vs SDE Binary: seeking recommendations, experiences

I've got some questions I've collected as we try to come to a decision about whether to use ST_Geometry or SDE native geometry for an Oracle/SDE datastore on a new project. The organization has quite a few legacy systems that use the SDE Binary type, and this would be the first to use ST_Geometry if we went that route. Apart from the ability to access spatial data in ArcSDE with SQL, what are the benefits and drawbacks of ST_Geometry? What sort of pitfalls have you experienced in implementing with one type or the other? Can more than one data type live in the same schema? Note: for this project the data in SDE will be all points and polygons.


Can I color different layers with the same pseudocolor color map?

I have multiple DEM layers in my map (downloaded from atlas.louisiana.gov/lidar), but I cannot seem to get the colors to line up. You can see in the attached image that the colors are not continuous between the layers, but they should be. (I've checked the values, and they are). I've tried several ways around this, but none of them seem to work. This data is rendered as Single Band Gray. 1) right clicking on one layer, selecting copy style, and then pasting it into the other layer. 2) saving the style from one layer and loading it into the other. 3) changing the color map to grayscle, inputting min and max, and then reverting back to pseudocolor. (that one seemed silly) It seems like the problem is that I can only do a pseudocolor with standard deviations, and not min/max values. Is that correct? There seems like there should be some way around that... Any ideas?


Distort image to display on leaflet.js based on known 16x16 point grid

I have an image that is a satellite picture, and with php gd I want to distort it so I can overlay it on a leaflet.js map. To make it a bit more simple, I have a grid of points, 16x16 which bind the image to the coordinates, like I chopped the image to 256 pieces and found appropriate coordinates for corners of each of those pieces. It's tedious and inefficient to draw 256 pieces simultaneously, hence I'd like to process the image on the server-side. Is it even possible to do it with just PHP's GD library or is it necessary to reinforce the arsenal? The projection I'm using (and would like to stick with) is standard leaflet.js EPSG:3857. How do I even start about such a task ?


ArcEngine draw a polyline or multipoint

Does anyone have a bit of code to plot a multiline? I just want to add some IPoints to a collection and display in a particular color and width, but I can't find a good way to do this. Right now I'm using this after a dynamic draw event: `dynamicDisplay.DrawMultipleLines(pointcollection);`   But there's no clear place to set the graphic attributes, and it plots incorrectly.


Can I Manually Add Wizards to a Map on CartoDB?

Some of my locations are not getting geocoded in spite of my checking spellings etc.. For instance, I'm trying to load some data on Chhattisgarh, an eastern Indian state and Lakshadweep, an island in the Arabian Sea which falls under the Indian territory. Need urgent help with this.


OpenLayers and the legend

Is it a possibility in OL to add a legend (symbology) directly to (under) particular layer (not only switch on/off the layer)?


QGis: Natural breaks (Jenks) - algorithm is buggy?

I don't know how the Natural Breaks (Jenks) algorithm works in QGis graduated symbolism, and I've deliberately not researched it in order to ask this question, as my point is based on what seems to be intuitively wrong, rather than what may be technically correct (but not useful or expected). So! If you have a range of values and constrain them to a few classes, Jenks does a nice job of splitting up the data:![A](http://i.stack.imgur.com/Mi2mz.png) But if you add further classes, Jenks just repeats the lowest value, often: ![B](http://i.stack.imgur.com/2ab2b.png) Since I don't know (yet - I will now research it) how Jenks works, it'd be remiss of me to say the algorithm is coded "wrongly", however the implementation is peculiar and unexpected. Possibly if it auto-constrained itself to the maximum number of independent classes? Or pooled identical- valued colour bands?


Inherited EditorServerCommandAction double posts returnMessage

I have a bunch of inherited EditorServerCommandActions in my WebADF solution and when ever I use the ReturnMessage inside EditorServerAction() it generates two javascript popups at the clientside. Is this a common know error or am I simply missing something?               public class CustomCommandAction : EditorServerCommandAction     {         protected override void EditorServerAction()         {             ReturnMessage("Foo");         }     }


Converting from DMS to DD using Python in Field Calculator?

I need to convert lat/long that is expressed as degrees, minutes, and seconds in the data _into_ decimal degrees. For example, in the data, they are listed as N335042.06 in the Latitude column, and W86031.04 in the Longitude column. I have done this problem before where I created a script which converted DMS to DD, and vice-versa, so I guess I could use bits from that. But the problem I'm having is how to ignore (for lack of a better word) the 'N' and 'W' in the data? Can I skip them? And DMS are listed all together without any symbols or spaces. Can I use `len()`, `range()`, `split()` to specify what part to read from the value? For example, can do the following? N335042.06 where, 33 = degrees 50 = minutes 42.06 = seconds ...? I came across this ESRI article, but it's in VB. Will probably use it as a reference but some of the terminology/syntax is different from Python. Final code that works!               # Pre-logic     def latDD(x):       D = int(x[1:3])       M = int(x[3:5])       S = float(x[5:])       DD = D + float(M)/60 + float(S)/3600       return DD          # Expression     latDD(!Latitude!)


How to convert a National Bridge Inventory(NBI) ASCII file to a shapefile

I have a problem. I have downloaded the NBI ascii file data from http://www.fhwa.dot.gov/bridge/nbi/ascii.cfm. Now I am not able to convert it into a shapefile. The ASCII file has a latitude longitude column. I understand that it can be imported in a GIS software using X,Y values. But I am not able to make out whether the latitude longitude is in DMS or decimal degrees. I am posting a sample here. Any help would be appreciated Sample Latitude: 40454200   Sample Longitude: 122190700 If there is any other way of converting it, I would like to know.


How to export data from DBF to Excel template using ModelBuilder?

I have created a model in ModelBuilder that uses the iterator tool to run through a polygon layer and clip a raster based on the polygon's unique name in the attribute table. The model then copies the rows of the clipped raster's attribute table and creates a dbf table saved with the polygon name. I now need the dbf table data to be exported or copied and pasted into an excel workbook that contains formulas in a separate worksheet to analyze the data. Is there any way to do this in ModelBuilder? Python? I have used the do.TransferSpreadsheet command in VBA to transfer data from MS Access to an Excel Template in past.. The command works great for automated analysis of complex data and would like something similar for ArcGIS.Thanks in advance!


Implementing WFS within .NET

There are several discussion about implementing (or not implementing) WFS with .net from scratch. However, if I do need to, how do I get started? I know it's a broad question. I apologize for that. I need help getting started. Any kind of walk through will be helpful. edited (8/2/13): thank you for all the comments and responses. I am looking for a step by step walkthrough with details assuming I'm just a GIS person with limited knowledge to a lot of programming concepts. Also, if this is not the right place to ask this question, please let me know.


Umlaute not displayed correctly when loading .mdb MS Access into QGIS (1.9.0-Master) with eVis

Good evening, I am loading an .mdb (MS Access 97 or 2003) to QGIS (1.9) via the eVis database conection. Connecting to the mdb, selecting a table and displaying it as a point-layer through the appropriate x/y-fields works fine. However in the attribute table the German Umlaute ä, ö, ü and also ß are displayed as inverted "?". After days of searching I don't really know how to tackle this one. I am grateful for any recommendation on how to fix this. The umlaute in the mdb / attribute table are displayed correctly in ArcMap. I would love to get this running in QGIS. Thank you.


Entity Framework & Sql Server 2012 spatial

I'm struggling to find any concrete performance of Entity Framework paired with SQL server 2012. We're looking at storing nearly a million points in a database an right now the most common scenario would be to find all points within radius of N. I've mocked up a small prototype using EF5 and it seems fine. However, a senior developer has mentioned that in his opinion EF is bloated and I should consider using a micro-orm or even pure ADO.NET as EF is a bloated monster (lol). Performance is crucial on this project and I would like to know if anyone has experience to prove this guy otherwise (or right for that matter)?


How to adjust raster size for usage in raster calculator?

I want to make a visibility analysis of planned wind turbines in my region. I have a DEM raster and a polygon describing the possible area where the wind turbines might be set up. People are complaining, that the view to a prominent mountain in the region will be disturbed. My "simple" plan is, to raise the DEM-values covered by the polygon up to 200 meters and then perform a visibility analysis from the top of that mountain to be able to identify the regions, where the wind turbines might "block" the view. Problem: I'm not able to successfully combine the vector and the DEM. Steps taken so far:   * Polygon was clipped from the DEM   * Raster values from the clipped part were set to 200 What doesn't work: When using the raster calculator, trying to add DEM + clip, the resulting image only covers the clip. As i understand, the extend of the two rasters should be the same, and that NoData values will prevent operations (googling for 3 hours without success). I fail in resizing my clip to the DEM and also do not get rid of the NoData values.(Really tried every tool under Raster to no avail) Unfortunately, the rastercalc tutorials i found do not cover these "basic" things, so i'm completely lost at the moment. Thanx Bernd p.s. using latest qgis master from osgeo4w


Can I make an ESRI checkout replica from a checkout replica?

If I create a checkout from a clients geodatabase, can I then turn around and put that replica on our ArcGIS Server and create a couple replicas from it? So, then, our replicas would be synced to the "replica" on our server which would then be exported back out to a personal geodatabase. Would that personal geodatabase be able to be synced back to the clients database or would doing all this ruin what makes it a replica? I hope this made sense. Basically, what I'm trying to accomplish is getting a replica from a client that multiple people will be working on, but before syncing it back to the client, we would like the replicas we worked on synced back to 1 database so we can check the connectivity before giving it back to the client.


How to compile QGIS 2.4.0 on Ubuntu Precise 12.04

Been trying to compile this for a couple of hours, using this download: http://qgis.org/downloads/qgis-latest.tar.bz2 and these instructions: http://htmlpreview.github.io/?https://raw.github.com/qgis/QGIS/master/doc/INSTALL.html#toc3. I have problems with `libproj-dev` - I can't install it, get this error: > The following packages have unmet dependencies: libproj-dev : Depends: > libproj0 (= 4.8.0-3~precise2) but 4.8.0-3~precise5 is to be installed E: > Unable to correct problems, you have held broken packages. So I installed the headers by hand - compiler accepted them - and have tried to link in /usr/lib/libproj.so as `PROJ_LIBRARY` in ccmake. The compiler doesn't like this, I get: > make[2]: *** No rule to make target`/usr/lib/libproj.so', needed > by`output/lib/libqgis_core.so.2.4.0'. Stop. make[1]: *** > [src/core/CMakeFiles/qgis_core.dir/all] Error 2 I also get 'Could not find GRASS' a little earlier in the build - don't know if this is significant. I'm a bit confused, `libproj.so` shouldn't be a target? Anyway, can anyone help? My ppa is the ubuntugis stable one - maybe I need unstable?


PostGIS 2.0 "No function matches the given name and argument types"

After a lot of research, I figured out how to build PostGIS 2.0 from source on a Turnkey Linux LAPP (Postgres 8.4) server. The main draw of PostGIS 2.0 for me is the ST_ConcaveHull function, aka "shrink wrap". Now that I have a functioning PostGIS 2.0, I'm trying to implement these functions. I'm working with the NWS Shapefiles for States and Forecast Zones. Most if not all of the geometries I'm dealing with are multipolygons. I've figured out how to draw these polygons, and this is the query I use to request data in the GeoJSON format:               SELECT ST_AsGeoJSON(geom) as geom FROM states WHERE id = 5;      Unfortounately, because I'm working with State sized data, some of the requests are HUGE! The GeoJSON for New Jersey is 2mb! I know PostGIS has some great functions to reduce the complexity of the geometries, but I'm having trouble implementing them. The functions I want to use are:   * ST_Simplify    * ST_Union    * ST_ConcaveHull Unfortunately, I'm getting some very unhelpful error messages when I run queries using these functions.               SELECT ST_AsGeoJson(ST_Simplify(geom)) as geom FROM states WHERE id = 5;          Error: SQLSTATE[42883]: Undefined function: 7 ERROR: function st_simplify(geometry)          does not exist LINE 1: SELECT ST_AsGeoJson(ST_Simplify(geom)) as geom FROM states W...      ^ HINT: No function matches the given name and argument types. You might need to add      explicit type casts.      I went through a LOT of time and effort to make sure that I had PostGIS 2.0 instead of 1.5 for the explicit purpose of using things like ST_ConcaveHull and ST_Simplify, and now Postgre is telling me they don't exist! I found another post on here that gave me a query listing all global functions available to Postgre:               SELECT  proname, proargnames, prosrc      FROM    pg_catalog.pg_namespace n     JOIN    pg_catalog.pg_proc p     ON      pronamespace = n.oid     WHERE   nspname = 'public'      Listed in the results are:   * st_union {geom1,geom2} geomunion   * st_union pgis_union_geometry_array   * st_union aggregate_dummy   * st_concavehull {param_geom,param_pctconvex,param_allow_holes}   * st_simplify LWGEOM_simplify2d   * st_simplifypreservetopology topologypreservesimplify So I'm kind of stuck and confused at this point. The query results say the function doesn't exist, but the output of the above query says they DO exist. Can anybody shed some light on this for me? What am I doing wrong?


Add a popup to a Vector LineString in OpenLayers?

Is it possible to add a popup (without a marker icon) to a Vector (LineString)? This is my code, if possible, how to change it to add the popup? I need the popup to display some additional information to the road segement.               <script>     var ytrajectory=[];     ytrajectory.push('LINESTRING (11.1119339 46.0772706, 11.1117177 46.077171, 11.111613 46.0770837, 11.1115376 46.07702)');     ytrajectory.push('LINESTRING (11.1115376 46.07702, 11.1114585 46.0768611, 11.1114219 46.0767673)');          ...     </script>     function showReconstructedTrajectoryy() {     var vectorLayer = new OpenLayers.Layer.Vector(     "Reconstructed trajectoryy", {     style : reconstructedTrajectoryStyleYellow     });          var parser = new OpenLayers.Format.WKT();     var src = new OpenLayers.Projection('EPSG:4326');     var dest = new OpenLayers.Projection('EPSG:900913');          var s, wktLine, features, reprojectedFeatures;     for (s = 0; s < ytrajectory.length; s++) {     wktLine = ytrajectory[s];     features = parser.read(wktLine);     reprojectedFeatures = features.geometry.transform(src, dest);     vectorLayer.addFeatures(new OpenLayers.Feature.Vector(     reprojectedFeatures));     }          osmmap.addLayer(vectorLayer);     osmmap.zoomToExtent(vectorLayer.getDataExtent());     }


Legend in FlexViewer - Only show classes that are visible in the current map extent

I am using ArcGIS Viewer for Flex to create web map. The legend widget can be set so that the legend in web map to only show **layers** that are available at the current map scale. This is great but I want to set the legend to only show **_classes that are visible in the current map extent_**. Does anyone know how to do it? Thank you


How to add logical operators in wfs openlayers

How i can add logical operators AND, OR both togehter in WFS openlayers.                   <script type="text/javascript">         var map;              function init() {                  var bounds = new OpenLayers.Bounds(                     2309872.5, 2399960.25,                     2388258.25, 2450114.75                 );                     var options = {                     controls: [                                 new OpenLayers.Control.Navigation(),                                 new OpenLayers.Control.PanZoom()                          ],                     maxExtent: bounds,                     maxResolution: 306.1943359375,                     projection: "EPSG:32644",                     units: 'm'                 };                 map = new OpenLayers.Map('map', options);                 var  map_wms = new OpenLayers.Layer.WMS(                     "cite:union - Tiled", "http://localhost:8080/geoserver/cite/wms",                     {                         LAYERS: 'cite:union',                         tiled: true,                         EPSG:32644,                         tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                     },                     {                         buffer: 0,                         displayOutsideMaxExtent: true,                         isBaseLayer: true                     }                  );                            var layerPolygon2 = new OpenLayers.Layer.Vector("Layer Information", {                                             strategies: [new OpenLayers.Strategy.BBOX()],                                             visibility: true,                                             protocol: new OpenLayers.Protocol.WFS({                                             version: "1.0.0",                                             srsName: "EPSG:32644", // this is the default                                             url:  "http://localhost:8080/geoserver/cite/wfs",                                             featureType:"union",                                             featureNS: "http://www.opengeospatial.net/cite",                                             extractAttributes: true,                                             transparent:true,                                             geometryName: "the_geom"                                                       })                          filter: new OpenLayers.Filter.Logical({                     type: OpenLayers.Filter.Logical.AND,                     filters: [                         new OpenLayers.Filter.Logical({                         type: OpenLayers.Filter.Logical.OR,                         filters: [                             new OpenLayers.Filter.Comparison({                             type: OpenLayers.Filter.Comparison.EQUAL_TO,                             property: "classes",                             value: "A?Pt"                             }),                             new OpenLayers.Filter.Comparison({                             type: OpenLayers.Filter.Comparison.EQUAL_TO,                             property: "classes",                             value: "P"                             })                             ]                         )}                         new OpenLayers.Filter.Logical({                         type: OpenLayers.Filter.Logical.OR,                         filters: [                             new OpenLayers.Filter.Comparison({                             type: OpenLayers.Filter.Comparison.EQUAL_TO,                             property: "luse_2005",                             value: "BU"                             }),                             new OpenLayers.Filter.Comparison({                             type: OpenLayers.Filter.Comparison.EQUAL_TO,                             property: "luse_2005",                             value: "CL"                             })                             ]                         )}                          ]                      })                           });                                map.addLayers([map_wms,layerPolygon2]);                      map.zoomToMaxExtent();      } when i try to run above code in browser nothing is being displayed and firebug showing ReferenceError: init is not defined init() error. thanks!!


Print XYZ layer with MapFish Print plugin

I have a XYZ layer(i hope i call it right) in GeoExt/gxp application:               ,{source: "ol",group: "background",fixed: true,type: "OpenLayers.Layer.XYZ",args:        ["ЕЭКО","http://maps.rosreestr.ru/ArcGIS/rest/services/BaseMaps/BaseMapAnno/MapServer/tile/$\{z\}/$\{y\}/$\{x\}",{visibility: true}]}      In MapFish documentation says that print plugin can print TMS and Tile Cache. But its print only WMS layers. I gonna set url to XYZ layer or something? **UPDATE** I found working XYZ print in here http://map.geo.admin.ch/. But There is something strange with layer which set to print. **UPDATE2** More interesting that GeoExt.PrintProvider not send XYZ layer to print plugin. But if i give to it gwc layer i can see gwc layer to post parameters: ![enter image description here](http://i.stack.imgur.com/60ras.jpg) But this layer not print. Why PrintProvider not send XYZ layer in post parameters?


MapInfo 10.0 - saving a MapInfo browser to as a table

I geocoded a table in MapInfo 10.0. In the browser window, the lat and long are viewable, however the columns are not saved to the original table. The "Save Table" option is not available. Are there any options of saving a browser as a table?


How to use the MMQGIS plugin through the Python console?

I know how to make Hub Lines using the GUI, I want to know how to do the same using the Python Console, a tutorial would be great. Thank you


Calculate "CO_FIPS" field based on "COUNTY" field

I have several shapefiles with "COUNTY" already populated, I need to populate the "CO_FIPS" field based on the county. BREVARD = 009, BROWARD = 011 etc..


STIntersects - Update field to give me count of points in polygon

I have a points layer (dbo.ptLayer)   * Around 1M points   * Spatial Geometry Type (dbo.ptLayer.geom)   * No spatial index just yet, but will create one once data gathering complete. I have a polygon layer (dbo.polygonLayer)   * Around 500 polygons.   * Spatial Geometry Type (dbo.polygonLayer.geom) Both have fields called ID. **I want to populate an empty integer field in the polygon layer, with a count of the total number of points within each polygon.** Although I have access to other software products, I am interested to learn what can be done purely within SQL and SQL Server. I believe I should be making use of STIntersects but would like to know what is the best way of doing an update to populate this field.


Generating Point from Centroid of Polygon Error

I have a feature class of polygons within a geodatabase and I'm trying to generate points based on the centroid of the polygon using the Feature to Point tool in ArcGIS 9.3. However, I end up with several points that fall just outside of the polygon. So ... is this the wrong tool to be using for this application? Also, am I making the wrong assumption about this tool generating points from the centroids of polygons?


How to move/edit labels?

I need to move/edit few labels in a job. How can I work on it? I am using QGIS 1.7.4 and the label option in the toolbar (not the layer options)


Select Move Label Tool in QGIS

> **Possible Duplicate:**   >  How to move labels in QGIS? I would like to manually move 6 labels but the Move Label Tool is unable to be selected. I have placed labels by selecting the layer and clicking the "ABC" Labeling button. After going into Toggle Editing, I'm able to Change Labels but not Move or Rotate Labels. I'm running QGIS 1.7.4. I've looked at: How to move labels in QGIS? Is there a way to manually place labels in QGIS? It's a very simple problem but I can't seem to crack it. Thanks all, Zach


Table Manager - Shapefile (Labels)

I have a point shapefile which is an output of the plugin table manager. Now I’m labeling the features, but i can´t seem to get the move label bottom to work, and is active. Is there any connection? Thanks in advance


How to add points to a point vector layer - QGIS

I have a point layer with multiple data points in QGIS, but I want to add some extra data points to it and adjust some info that we missed because the ground was wet and we needed to go around it. I can't figure it out. I would think you would go in edit mode and have an option to add points, but can't find it.


How to add a title in a GetLegendGraphic Request

The GetLegendGraphic request works well. I display the legend at a fixed position in my viewer using OpenLayers. But I would like to know if there is any way to add some custom text to the legend. The legend that geoserver produces is like this: http://docs.geoserver.org/latest/en/user/_images/samplelegend.png But I need to produce a legend with a title: http://minus.com/m4AP102Jj Does anyone know if the GetLegendGraphic request in GeoServer, or the SLD style have a way to add this? **UPDATE (23/08/12)** For now I´ll use a litle trick to add a title by adding in my SLD another custom rule above all others that not match any other rules: Example:               <sld:Rule>       <sld:Name>[NAME]</sld:Name>         <sld:Title>[TITLE]</sld:Title>         <ogc:Filter>           <ogc:PropertyIsEqualTo>             <ogc:PropertyName>[PROPERTY_NAME]</ogc:PropertyName>             <ogc:Literal>[SOME_VALUE_THAT_NOT_MATCH_OTHER_RULES]</ogc:Literal>           </ogc:PropertyIsEqualTo>         </ogc:Filter>         <sld:PolygonSymbolizer>           <sld:Fill>             <sld:CssParameter name="fill">[MATCH_BACKGROUND_COLOR]</sld:CssParameter>           </sld:Fill>         </sld:PolygonSymbolizer>           </sld:Rule>


Does the USDA SSURGO database have a web-based API that I can query?

Based on a previous question, Can I download the entire NRCS SSURGO database?, it appears that it is difficult to import the entire USDA SSURGO database. The Web Soil Survey provides the ability to draw a polygon and then download data. But the polygon has to be fairly small, and there is a lot of pointing, clicking, and waiting to get the desired data. There is also a Geospatial data gateway where it is possible download county-level data. However, these do not meet my needs. This is what I want to do:   1. I start with the lat/lon coordinates of the corners of a rectangle   2. find the area of each soil type in the rectangle   3. find % clay, sand, silt, organic matter as a function of depth for each soil type  Is there a way that I can do this programatically? (Solutions in R preferred, but any open-source solution will be useful)


How can I consume map service using static IP address for VPN

I'm current using the url (http://services.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer) for my iPad map service. However, we need to connect to VPN to access the internal network at the same time. And the VPN server blocks the public Internet and only allow me to setup IP Address based public Internet access. So my question is there static IP addresses associated to the above url I could specify? Thanks


Calculate Distance between a Polygon and a Point Along a Path

I am very new to ArcGis and am trying to figure out how to do this. Basically, this is what I am trying to do: I have a map of Canada and in it I have polygons representing areas (each about the size of a municipality). I also have points describing the location of every airport in Canada and the network of roads all across Canada. **What I need to do is find the distance from each area to the nearest airport using a road.** What I tried to do already was find the centroid of each polygon because I thought that I could use that to calculate the distance between two points along a road. But then I realized that the centroid is not necessarily connected to a road (especially in the northern parts of Canada) therefore I cannot use that. I would really appreciate some sort of guidance as to how I would accomplish this as I am very new to ArcGIS and have not found much to help me (without confusing me) online.


Using VBScript to label shapefile in ArcMap

If you're reading this for the first time, it might be best to focus on section beneath **ANOTHER UPDATE**. I'm a beginner with VBScript and I have a question about using an VBScript expression for a label. I want to label street address numbers and the fractional address number, if it exists, in my map. This data is stored in an address points shapefile. I go Properties > Labels > Expression... in Esri software to enter this expression:               [HSE_NO]&" "& [FRAC_NO]       With this expression the labels will be created next to the address points as follows:                412 1/2 412      See how both labels are on the same line? I would like to have the labels look like this:               412     412 1/2      I want the first label to be the address value without a fraction, and the following lines would be address values with fractions. I welcome your suggestions. **UPDATE** I think what I need is more complex than I realized. As I mentioned, I want first label to be address street value (HSE_NO) without a fraction (where FRAC_NO is NULL). I'm working on getting this to occur with what's below. It's not correct! Any ideas?               Function FindLabel ( [HSE_NO], [FRAC_NO] )     if (IsNull[FRAC_NO]) then         FindLabel = [HSE_NO]          End If     End Function      **ANOTHER UPDATE** I have a point shapefile of addresses. Some address points contain the address number and street name. Some points contain address number, address fractional number, and street name. Points containing fractional values will overlap points containing not containing fractional values. I want to label these overlapping points like this:               1159     1159 1/2      My current VBScript labels addresses in this manner.               1159 1/2 1159      Does anyone know how to get the desired result? Here's my current expression:               Function FindLabel ( [HSE_NO], [FRAC_NO] )       if (IsNull([FRAC_NO])) then         FindLabel = [HSE_NO]       else         FindLabel = [HSE_NO] & " " & [FRAC_NO]        End If     End Function      Let me know if there are questions. Thank you.


Trouble exporting from QGIS to Illustrator

I am having trouble getting my files to transfer or export to a vector file. When I open the file in Adobe illustrator there is a lot of information missing. I have tried to "save selection as a vector..." I have all my layer selected. I tried selecting individual layers, all the layers at once, I selected the layers individually. Next I tried to play with the CRS pull down menu. I used all the options and still get no results. I finally got what may be a file but when i open it in illustrator I get an error message that says "file cannot be opened because there is not enough RAM,, I have 12 gigs of RAM, should not be a problem. Am I missing something? A plug-in?


How to change the raster NoData value with ArcObjects?

I am trying to change the a raster's NoData value with ArcObjects and keep coming up short. Here is the code that I thought would work:               IRasterDataset _dataset = ...get the dataset...     object[] noDataValue = ...set the no data values by raster band...     var rasterbandCollection = (IRasterBandCollection) _dataset;     for(int i = 0; i < rasterbandCollection.Count; i++)     {         ((IRasterProps) rasterbandCollection.Item(i)).NoDataValue = noDataValue[i];     }      This value persists during the ArcMap session, but does not save the changes to disk it appears. Any ideas on how to make this permanent (without having to save a second raster dataset)?


ArcSDE: Specifying fields to not allow null values

I have a feature class in SDE that has various fields where the setting AllowNulls is set to FALSE. When creating a feature in ArcMap in this feature class it still allows nulls to be specified in the attributes, albeit an empty text string. I can then save my edits without any error messages. What's going on? Thanks Rob


Adding Google map layer in Openlayers - what about this function?

How to add Google map as a layer (something like " _var gmap = new OpenLayers.Layer.Google("Google Streets_ ", {visibility: false})"?) to this _mymap.js_ function i already have:               var map = new OpenLayers.Map({         div: 'map',     //  allOverlays: true,         theme: null,         projection: sphericalMercatorProj,         displayProjection: geographicProj,         units: 'm',         numZoomLevels: 18,         maxResolution: 156543.0339,         maxExtent: new OpenLayers.Bounds(             -20037508.34, -20037508.34, 20037508.34, 20037508.34         ),         controls: [             new OpenLayers.Control.Attribution(),             new OpenLayers.Control.Navigation(),             new OpenLayers.Control.PanZoom(),             new OpenLayers.Control.LayerSwitcher()              ],         layers: [             new OpenLayers.Layer.OSM('OpenStreetMap', null),             sundialsLayer,             sprintersLayer,             poisLayer         ],         center: new OpenLayers.LonLat(0, 0),         zoom: 2     });      Thanks. Note: suggested answer is fine -                   layers: [             new OpenLayers.Layer.OSM('OpenStreetMap', null),             new OpenLayers.Layer.Google(             "Google Streets", // the default             {numZoomLevels: 20},             )                  sundialsLayer,             sprintersLayer,             poisLayer         ],


Make Table View with field_info Output

I am using Make Table View to copy data from one table to a table view object and in the process hide and/or rename some fields. When I output the data to a dbf it hides and/or renames the fields accordingly. When I output the data to a geodatabase, it hides the specified fields but does not rename the ones that should be renamed. Any ideas why this is happening? EDIT: After some further testing, it looks like it is changing the field name, and not the field alias in the geodatabase. I tried changing field.name to field.aliasName, but that did not work.               import arcpy     from arcpy import env          # Set the enviroment     env.workspace = '//GIS-SDE/ParcelData.mdb'     env.qualifiedFieldNames = False          # Get the fields from the G_AccountInfo table in ParcelData.mdb     fields = arcpy.ListFields("G_AccountInfo")          # Create a fieldinfo object     fieldinfo = arcpy.FieldInfo()          # Iterate through the fields and set them to fieldinfo     for field in fields:         if field.name == "ACCOUNT_NUM":             fieldinfo.addField(field.name, "ACCOUNT", "VISIBLE", "")         elif field.name == "OWNER_NAME1":           fieldinfo.addField(field.name, "OWNERNAME", "VISIBLE", "")         elif field.name == "OWNER_ADDRESS_LINE1":           fieldinfo.addField(field.name, "OWNERADDRESS1", "VISIBLE", "")         else:             fieldinfo.addField(field.name, field.name, "HIDDEN", "")          arcpy.MakeTableView_management ("G_AccountInfo", "AccountInfo", '//GIS-SDE/ParcelData.mdb', fieldinfo)          # Give the output table a unique name     output = '//GIS-SDE/ParcelData.mdb/cleaned_data'          # Add the table back to the gdb     arcpy.CopyRows_management("AccountInfo", output)


Map Info Workspace problem

I have opened a workspace with a file on a CD and now it will only open when the CD is in the CD drive - I can't find a way around it - any help gratefully received , thanks


map of europe in list of lat/long points

I'm looking for a list of lat/long coordinates which will describe all countries of europe. Ideally the list will also describe on which country each closed path belongs to. I'm imagining something along the lines:                lat    long    cc  type     ------ ------- ---- ----     55.05   -7.25   IE  EU27     55.18   -6.97   IE  EU27     55.25   -6.48   IE  EU27     55.22   -6.1    IE  EU27     55.05   -6      IE  EU27     54.83   -5.7    IE  EU27     54.63   -5.92   IE  EU27     ...              55.23   -6.93   IE  EU27     55.05   -7.25   IE  EU27          58.63   -5      UK  EU27     58.55   -4.67   UK  EU27     58.45   -4.77   UK  EU27     58.58   -4.6    UK  EU27     58.45   -4.48   UK  EU27     58.6    -4.02   UK  EU27     58.63   -3.02   UK  EU27     58.37   -3.12   UK  EU27     57.85   -4.2    UK  EU27     57.85   -3.77   UK  EU27     57.6    -4.42   UK  EU27     57.6    -4.1    UK  EU27     ...               Regarding accuracy, I'd be happy with anything over 5.000 points in total. Where should I look at? I will convert those coordinates to x/y using `fromLatLngToContainerPixel` of google maps and then draw the map using javascript.


How can I convert geomagnetic coordinates to geographic coordinates without doing the math myself?

I have some data which are in geomagnetic coordinates; that is, they are latitude and longitude, but in reference to the magnetic north pole in Canada rather than the geographic north pole. As you might guess, I'd like to match them to other data which are in other coordinate systems. (The data originated as a computation of the auroral oval; i.e., they originated as geomagnetic and there is no geographic source I could consult, even in principle.) The data will be stored in PostGIS. The easy way to do this is to find an EPSG code for the geomagnetic coordinate system, tag the data with that code, and then everything is magically reprojected to whatever I need at the time. However, I can't find any EPSG codes on spatialreference.org containing either _geomagnetic_ or _magnetic_. Various Google searches also turn up nothing. For the level of precision I need, any latitude/longitude coordinate system (e.g., WGS84) modified to have a different north pole is good enough. So if there's an easy way to define such a coordinate system and load it into PostGIS, that would be OK too. **Clarification** For example, suppose that instead of converting from geomagnetic coordinates to geographic, I simply wanted to convert from UTM 15N to WGS84. I could save the UTM object into PostGIS unmodified, and tag it with the EPSG code 26915. Then, when I later made a query, I could say to PostGIS, "please return results in EPSG code 4326", and PostGIS would convert to WGS84 automatically. Similarly, I would like to tag the geomagnetic data with some coordinate system and let PostGIS do the reprojection behind the scenes, without me having to call a `reproject()` function manually.


Map Projection - where to begin

First of all, I realize that this question can be considered like "open" or maybe "too wide" (and definitely very basic for all of you). However, I am completely new in this and I am not looking for more then several hits where to begin, which topics I should search for. Also, I search for some similar questions and I found few open questions with similar meaning (Java, GIS knowledge…What next?, The Best books from ESRI, Where can I find QGIS tutorials and resources?) so I hope someone can help. **The Issues:** * * * I have data in geographical coordinates which I would like to project them into geomagnetic coordinates. * * * What do I need for that? Is information on the internet enough to solve this or I should search for a book? What should I search for?


Also Cant see WFS-T layer with OpenLayers and GeoServer (Windows)

I am in the process of learning for doing WFST - OpenLayers in web. I followed this tutorial in the GisTutor. But i stuck with 2 problems   * the WFS-T layer is not appearing   * the feature created in client cannot be saved (no response when save button is clicked, and accordingly there is no additional row in the PostGIS table of WFST_TEST) I followed this GIS.SE question and I also do not configured the Proxypass as instructed here (also because I do not have apache running the GeoServer). In the latter part of the discussion between Janne Mattila with Iant, Iant suggests that: > don't run from file://somefolder/mytestpage.html always use > localhost:28080/geoserver/www/mytestpage.html this should fix your problem. > – iant and then Janne seems to be happy: > Nice and simple, that fixed it! I don't fully understand why it works like > that, though. Can you explain what happens in OpenLayers when I access the > page using file:// and when I use localhost:28080 instead? – Janne Mattila I finally able to call my html by `localhost:8080/geoserver/www/rinotestpage.html` in my browser, by putting my `rinotestpage.html` and the `THEME` folder (for bringing button PNG files to the app) into the `C:\Program Files\GeoServer 2.3.2\data_dir\www`. But my 2 problems in the above are persisting (WFS-T layer not appeared and feature created in client is not saved) My question:   1. Janne and Iant, how did you crack this? How to solve the problem (WFS-T layer not appearing)?   2. Janne, did you put your `testpage.html` in the same folder as i did ? (`C:\Program Files\GeoServer 2.3.2\data_dir\www`). For additional infos:   1. I am using : PostgreSQL 9.1.9.1 Windows 32-bit + PostGIS 2.0.3.1 + GeoServer 2.3.2 + Windows Server 2003 in VMWare Player 5.0.2   2. I can see in my firebug that there is "POST wfs" request to my `http://localhost:8080/geoserver/wfs`, but the request took forever and no response. I am guessing this is why my WFS-T layer is not shown. ![enter image description here](http://i.stack.imgur.com/Bb3De.jpg) many thanks /Rino


Trying to display some data in OpenLayers - it's not erroring, just not loading :(

I'm trying to display some KML data on a map, using OpenLayers. Nothing complex - but it's just not rendering. Here's my javascript code               var map;          function init(element, kml) {         map = new OpenLayers.Map(element, {             projection: 'EPSG:3857',             layers: [                 new OpenLayers.Layer.Google(                     'Google Streets', // the default                     { numZoomLevels: 20 }                 ),                 new OpenLayers.Layer.Vector('Kml data')             ],             center: new OpenLayers.LonLat(153.132, -26.987)                 // Google.v3 uses web mercator as projection, so we have to                 // transform our coordinates                 .transform('EPSG:4326', 'EPSG:3857'),             zoom: 10         });              // Grab the features.         var features = convertKmlToFeatures(kml);              // Display it.         map.layers[1].addFeatures(features);     }          function convertKmlToFeatures(kml) {         var kmlFormat = new OpenLayers.Format.KML();         return kmlFormat.read(kml);     }      Basically,   * Create a map with two layers.   * Convert the Kml to features.   * Display the features on the map.  I just don't see what I've done wrong :( Oh - on a side note, the KML data is fine, because when I tried it to visualise it on this website, it worked. ![enter image description here](http://i.stack.imgur.com/6cvfa.png)


Problem generating Geodatabase Diagrammer

I am using ArcGIS 10.1 and Visio 2010. When I am trying to generate the diagram using the geodatabase diagrammer it has the following to error messages: Source at Microsoft.Office.Interop.Visio.DocumentsClass.Add(StringFileName) at GeodatabaseDiagrammer.VisioDiagram.OpenVisioDrawing(String fileName, String pathName) in C:\agis\Tools\ArcGIS\10x\GeodatabaseDiagrammer\AddInSourceCode\Version1.2\GeodatabaseDiagrammer\Helpers\VisioDiagram.vb:line 102 Msg: File not found. Then it would pop up another dialog box and say Source: at GeodatabaseDiagrammer.VisioDiagram.SectionHeader() in C:\agis\Tools\ArcGIS\10x\GeodatabaseDiagrammer\AddInSourceCode\Version1.2\GeodatabaseDiagrammer\Helpers\VisoDiagram.vb:line 160 Msg: Object reference not set to an instance of an object.


QGIS regular points around each polygon in a layer (vs ALL polygons)

I have a number of buffered points (some where polygons intersect & are dissolved. I'd like to create regular spaced points grid around EACH polygon -- not ALL polygons. When I do Vector -> Research Tools -> Regular Points --- I have to choose a layer as an Input Boundary, but it doesn't seem to loop through each shape to create a GRID just for that shape, and instead creates one large grid around ALL shapes. The end result: ![enter image description here](http://i.stack.imgur.com/2X4iY.jpg) How can I get a point grid for EACH individual shape vs ALL?


How can I add attributes contained in JSON data to an OpenLayers popup?

I am working on a website using jQuery, Javascript, and OpenLayers. I am adding markers to an OpenLayers Vector layer using the following code:               var pointLayer = new OpenLayers.Layer.Vector("POI Markers", {projection: "EPSG:4326"});     $.getJSON('json/DATA.json', function(data) {       $.each(data.transmitters, function() {         var pointFeatures = [];         var px = this.longitude;         var py = this.latitude;         // Create a lonlat instance and transform it to the map projection.         var lonlat = new OpenLayers.LonLat(px, py);         lonlat.transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:900913"));              var pointGeometry = new OpenLayers.Geometry.Point(lonlat.lon, lonlat.lat);         var pointFeature = new OpenLayers.Feature.Vector(pointGeometry, null, {             pointRadius: 16,             fillOpacity: 0.7,             externalGraphic: 'marker.png',         });              pointFeatures.push(pointFeature);         pointLayer.addFeatures(pointFeatures);       });     });      I am using the following code to select markers and create popups:                var selectControl = new OpenLayers.Control.SelectFeature(pointLayer, {       hover: true,       onSelect: function(feature) {         var layer = feature.layer;         feature.style.fillOpacity = 1;         feature.style.pointRadius = 20;         layer.drawFeature(feature);                  var content = "
Feature:
" + feature.id +                   "

Location:
" + feature.geometry +"
";                  var popup = new OpenLayers.Popup.FramedCloud(             feature.id + "_popup",              feature.geometry.getBounds().getCenterLonLat(),             new OpenLayers.Size(250, 100),             content,             null,              false,              null);             feature.popup = popup;             map.addPopup(popup);         },         onUnselect: function(feature) {             var layer = feature.layer;             feature.style.fillOpacity = 0.7;             feature.style.pointRadius = 16;             feature.renderIntent = null;             layer.drawFeature(feature);                  map.removePopup(feature.popup);         }     });     map.addControl(selectControl);     selectControl.activate();      As you can see, the popup contains information regarding the feature id and geometry. What I would like to do is add information that is contained inside the DATA.json file referenced in the first code block to the corresponding popup. I am using a large JSON file. The following is a sample:               "LOCATION_A": {           "latitude": 10.2070,           "longitude": 25.3215,           "altitude": 1.2,           "status": "CLOSED",           "serverAddress": "192.168.0.1"      }      I would like to add the latitude, longitude and status information to the corresponding popup so when the user selects a popup, they view this information. In addition, it would be useful to change the "externalGraphic" file to reflect the status of each LOCATION. Continuing to use the above JSON file, if a POI were closed, the "externalGraphic" file should be one PNG file, while if the POI were open, it should be another. I am using jQuery in this project and find it fascinating! Any assistance with this problem is greatly appreciated. Thank you all so much.


i have an unhandled exception in asp.net

i am going to create a website and load a gis map to it. i am using arcgis 10 and i created a web application in arcgis server 10 i created a new WebSit in visual studio 2010 and add a map resource manager. add a MapResourceItem for this. after that add a map control and set MapResourceManager property. i config ArcGISIdentity too. but unfortunately when i run application occurs this error > Error processing server request > > [SoapException: Error processing server request] > System.Web.Services.Protocols.SoapHttpClientProtocol.ReadResponse(SoapClientMessage > message, WebResponse response, Stream responseStream, Boolean asyncCall) > +4159 System.Web.Services.Protocols.SoapHttpClientProtocol.Invoke(String > methodName, Object[] parameters) +328 > ESRI.ArcGIS.ADF.ArcGISServer.MapServerProxy.ComputeScale(MapDescription > MapDescription, ImageDisplay MapImageDisplay) +117 > ESRI.ArcGIS.ADF.Web.DataSources.ArcGISServer.MapFunctionality.GetScale(Envelope > extent, Int32 mapWidth, Int32 mapHeight) +564 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map.get_Scale() +794 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map.setExtent(Envelope value, Boolean > addSetExtentCallbackResult) +2556 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map.set_Extent(Envelope value) +39 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map.validatePrimaryMapResource() +2586 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map._OnLoad(EventArgs e) +115 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map.OnLoad(EventArgs e) +37 > System.Web.UI.Control.LoadRecursive() +134 > System.Web.UI.Control.LoadRecursive() +232 > System.Web.UI.Control.LoadRecursive() +232 > System.Web.UI.Page.ProcessRequestMain(Boolean includeStagesBeforeAsyncPoint, > Boolean includeStagesAfterAsyncPoint) +3412


Pointstack with different count-1-symbols / SLD / GeoServer

I try to use the pointstack sld to visualize my data at different scales. I used the example from http://suite.opengeo.org/opengeo- docs/cartography/rt/pointstacker.html, and everything is fine. However, I have a dataset of several thousand points, further defined by specific IDs. E.g. police stations (1402) or fire brigades (1401). Using the example from opengeo I can use only one symbol (or in this case an external vector graphic), but **I would like to have the specific symbol of e.g. fire brigade or police if the count is 1**. If I add a filter for the specific ID, the SLD-validation in GeoServer is fine but the filter is ignored. Therefore all points of the dataset (also non- police) are shown at the layer preview. So my idea is to set a filter at first (only use points with the ID 1401) and use the pointstack filter next. But this does not work. If I do not use the ID filter, all points of the dataset with the count = 1 are getting e.g. the police symbol. ![enter image description here](http://i.stack.imgur.com/qD5y5.jpg) The picture shows some stacked points and a lot of police stations (count = 1). However there aren´t so much police stations. Normally I would add a filter to check only points with the requested ID, so other points should be ignored for visualization (or) I can define point stacks for different IDs with different symbols. Maybe such a filter is not working in combination with the pointstack rule? I hope someone can help me :)


How to compare two rasters and reclass the greater value(python script)?

I am working on an assignment that requires us to reclass a raster(via python) from two rasters. We need to compare the power received at one cell(on raster1) and compare it to the same cell on raster 2. If the value of that cell in raster1 is greater than the value in raster2, the cell would be assigned a value of "1" while a value in raster2 greater than raster1 would get a "2". We also need to reclassify certain values as "0" since they are considered out of range. So the final raster would have three values: 0-out of range/no service 1-better value in raster1 2-better value in raster2 I can try to clarify further if needed.


make a selection by attributes of the fourth part of a range of numbers with arcpy

I am using Arcmap 10.1 and I need to create a script (for use on shapefiles) where I have a input field named "sequence" which is numerical data. I need to select the fourth part of a range of numbers. For example, 1-100 depending on the amount of data in the shapefile and I need to make a selection by attributes of "sequence" where I divide 100/4 = 25 and the selection is the first 1-25, the second selection of 26-50, the third 50-75, etc. The number of records can change and be 150/4 = 37 and the selection would be 1-37 the second third etc 38-75![enter image description here](http://i.stack.imgur.com/KrLZL.png) I have the ide of using arcpy.GetCount_management for me to total the records but do not know how to use the resuldato to divide by 4 and then put the result and make the selection with arcpy.SelectLayerByAttribute_management but not how to do this and if this is the right way to do


How to copy or export list filter from one model to another?

I have a model tool with a variable parameter that is filtered by a string list (county names) I manually created this list filter when I created the model tool. I am now creating a new model tool and would like to use the same list of county names is there anyway to import the list filter from my other model so I do not have to manually create the list again. I am an amateur at best with python so if it is possible to keep the solution within the realm of ModelBuilder that would be best but if python is neccessary bear with me. This is the value list I would like to reuse ![enter image description here](http://i.stack.imgur.com/5qiq2.png)


How to get the feature set for geoprocessing service?

I am calling a geoprocessing service in a web app. The geoprocessing service has one parameter (a feature set). But I am not sure how to find the parameter in the javascript code to send it to the geoprocessing since I am using an editor widget to first create the polygon and then run the geoprocessing service. Can you guys take a look at my code and tell me how I reference a feature set from the editor widget results? Thank you!               ], function(         Map, BasemapGallery, arcgisUtils, Edit,          ArcGISTiledMapServiceLayer, FeatureLayer,         Query, QueryTask,         SimpleMarkerSymbol, SimpleLineSymbol,          Editor, TemplatePicker,          esriConfig, jsapiBundle,         arrayUtils, parser, keys,                borderContainer, contentPane, titlePane,         dom, on       ) {         parser.parse();                     // snapping is enabled for this sample - change the tooltip to reflect this         jsapiBundle.toolbars.draw.start = jsapiBundle.toolbars.draw.start +  "
Press ALT to enable snapping";              // refer to "Using the Proxy Page" for more information:  https://developers.arcgis.com/en/javascript/jshelp/ags_proxy.html         esriConfig.defaults.io.proxyUrl = "/proxy";                       esriConfig.defaults.geometryService = new esri.tasks.GeometryService("http://tfsgis-iisd01:6080/arcgis/rest/services/Utilities/Geometry/GeometryServer");              map = new Map("map", {            basemap: "hybrid",                center: [-98.57, 30.98],           zoom: 6,           slider: true          });              map.on("layers-add-result", initEditor);                   var ActivityArea = new FeatureLayer("http://tfsgis-iisd01:6080/arcgis/rest/services/MyMapService2/FeatureServer/1",{           mode: FeatureLayer.MODE_ONDEMAND,            outFields: ['*']         });              var ActivityPoint = new FeatureLayer("http://tfsgis-iisd01:6080/arcgis/rest/services/MyMapService2/FeatureServer/0",{           mode: FeatureLayer.MODE_ONDEMAND,            outFields: ['*']         });              var stewardship = new FeatureLayer("http://tfsgis-iisd01:6080/arcgis/rest/services/MyMapService2/FeatureServer/2",{           mode: FeatureLayer.MODE_ONDEMAND,            outFields: ['*']         });              map.addLayers([ActivityArea,ActivityPoint,stewardship]);                             function initEditor(evt) {                var templateLayers = arrayUtils.map(evt.layers, function(result){             return result.layer;                });                     var templatePicker = new TemplatePicker({             featureLayers: templateLayers,             grouping: true,             rows: "auto",             columns: 3           }, "templateDiv");           templatePicker.startup();                var layers = arrayUtils.map(evt.layers, function(result) {             return { featureLayer: result.layer };           });                var settings = {             map: map,             templatePicker: templatePicker,             layerInfos: layers,             toolbarVisible: true,             createOptions: {               polylineDrawTools:[ Editor.CREATE_TOOL_FREEHAND_POLYLINE ],               polygonDrawTools: [ Editor.CREATE_TOOL_FREEHAND_POLYGON,                 Editor.CREATE_TOOL_CIRCLE,                 Editor.CREATE_TOOL_TRIANGLE,                 Editor.CREATE_TOOL_RECTANGLE               ]             },             toolbarOptions: {               reshapeVisible: true             }                     };                                    var params = {settings: settings};               var myEditor = new Editor(params,'editorDiv');           //define snapping options           var symbol = new SimpleMarkerSymbol(             SimpleMarkerSymbol.STYLE_CROSS,              15,              new SimpleLineSymbol(               SimpleLineSymbol.STYLE_SOLID,                new Color([255, 0, 0, 0.5]),                5             ),              null           );           map.enableSnapping({             snapPointSymbol: symbol,             tolerance: 20,             snapKey: keys.ALT           });                myEditor.startup();                        }                         on(dom.byId("button"),"click", function() {         gp = new esri.tasks.Geoprocessor("http://tfsgis-iisd01:6080/arcgis/rest/services/CalcFeaturesStewardshipOneParameter3/GPServer/CalcFeaturesStewardshipOneParameter")         //gp.outSpatialReference = map.spatialReference                   var features= [];              var selectGraphics = featureLayer.getSelectedFeatures();              features.push(selectGraphics);         var featureSet = new esri.tasks.FeatureSet();              featureSet.features = features;                   var params = {Stewardship:featureSet}              gp.submitJob(params, completeCallback , statusCallback)                        });              function completeCallback(jobInfo){            }            function statusCallback(jobInfo) {         var status = jobInfo.jobStatus;         if(status === "esriJobFailed"){           alert(status);              }         else if (status === "esriJobSucceeded"){           alert(status);         }       }


New to ArcGIS and want to compute distance between two points "as the fish swims"

We have a layer named Hydrography which is the water area for a given assignment. I have a station point and a target point. I need to compute the distance between those two points but inside of the hydrography. (boats cannot go over land) The DefenseService distance computation looks passable but it seems to be in a straight line since it does not take any Raster as input. The Euclidean seems to almost be there but I cannot find anything that is straight forward. I am using ArcGIS Explorer 2500, ArcGIS for Server 10.1 SP1 Enterprise, and we are writing our services and tools in C#. I have the loading of the hydrography, the creation of the IPoints, but at a loss as to what to do next! Any help appreciated, please.


Can't get MapServer WFS to display in OpenLayers

I'm new to OpenLayers and MapServer and I'm trying to serve a shapefile through a WFS but can't seem to get the data to display. When I use Firebug and HttpFox in Firefox I seem to be getting some data coming back, but it just won't get rendered on my web page. It is entirely possible that I have missed something simple but I can't see what. I've got the following in terms of my setup. MapFile:               MAP          NAME "MapServerWFS"     STATUS on     SIZE 600 600     EXTENT 0 0 700000 1300000     UNITS METERS     SHAPEPATH "shp"     IMAGECOLOR 240 240 240      PROJECTION          "init=EPSG:27700"     END          WEB         TEMPLATE "template.html"         IMAGEPATH "C:/DevProjects/fgsv2/output/"         IMAGEURL "../output/"         METADATA           "wfs_title"          "WFSDemo" ## REQUIRED           "wfs_onlineresource" "http://localhost/cgi-bin/mapserv.exe?map=C:/DevProjects/fgsv2/mapfile.map&" ## Recommended           "wfs_srs"            "EPSG:27700" ## Recommended           "ows_srs"            "EPSG:27700" ## Recommended           "wfs_abstract"       "This text describes my WFS service. This should be some counties." ## Recommended           "wfs_enable_request" "*"  # necessary           "wms_enable_request" "*"  # necessary           "wms_title"           "WMSDemo"         END     END          #this connects to a shapefile     LAYER         NAME "County"         DATA "County"              #CONNECTIONTYPE OGR         #CONNECTION "C:/devprojects/fgsv2/shp/County.shp"         STATUS ON         PROJECTION             "init=epsg:27700"         END         METADATA             "title" "County"             "wfs_title"         "County" ##REQUIRED             "wfs_srs"           "EPSG:27700" ## REQUIRED             "gml_include_items" "all" ## Optional (serves all attributes for layer)             "gml_featureid"     "db_id" ## REQUIRED             "wfs_enable_request" "*"         END              EXTENT 0 0 700000 1300000         TYPE POLYGON         CLASS             NAME "County"             STYLE                 color 100 250 250                 outlinecolor 50 50 250             END         END          PROCESSING "CLOSE_CONNECTION=DEFER"         dump true     END #end of layer                    END # end of map file      The data I'm using is English and Welsh Counties, so the projection is set to EPSG:27700. The OpenLayers map and layers bit is as follows:                               var farm_style = OpenLayers.Util.extend({}, OpenLayers.Feature.Vector.style['default']);                 farm_style.strokeWidth = 3;                 farm_style.strokeColor = "#ff0000";                 farm_style.fillOpacity = 0.5;                                 var mymap = new OpenLayers.Map({ div: "map",                     layers: [                     new OpenLayers.Layer.Vector("County", {                         isBaseLayer: true,                         strategies: [new OpenLayers.Strategy.BBOX()],                         styleMap: new OpenLayers.StyleMap({ "default": farm_style}),                          protocol: new OpenLayers.Protocol.WFS({                             url: "http://localhost/cgi-bin/mapserv.exe?map=C:/DevProjects/fgsv2/mapfile.map",                             version: "1.0.0",                             featureType: "County",                             featurePrefix:"ms",                             featureNS: "http://mapserver.gis.umn.edu/mapserver"                         })                     })                 ,                 ],                     center: [-1, 53],                     zoom: 8                 });      When I view the POST and Content through my Firefox addIns I can see POST:               <wfs:GetFeature service="WFS" version="1.0.0" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.0.0/WFS-transaction.xsd">     <wfs:Query typeName="ms:County">     <ogc:Filter>     <ogc:BBOX>     <ogc:PropertyName>the_geom</ogc:PropertyName>     <gml:Box srsName="EPSG:4326">     <gml:coordinates decimal="." cs="," ts=" ">     -2.64794921875,51.35205078125 0.64794921875,54.64794921875     </gml:coordinates>     </gml:Box>     </ogc:BBOX>     </ogc:Filter>     </wfs:Query>     </wfs:GetFeature>      The reference to EPSG:4326 is abit confusing as I've made no mention of it elsewhere. Is it a default setting somewhere? The Content tab contains the following:                   <wfs:FeatureCollection xsi:schemaLocation="http://www.opengis.net/wfs     http://schemas.opengis.net/wfs/1.0.0/WFS-basic.xsd     http://mapserver.gis.umn.edu/mapserver http://localhost/cgi-bin/mapserv.exe?map=C:/DevProjects/fgsv2/mapfile.map&SERVICE=WFS&VERSION=1.0.0&REQUEST=DescribeFeatureType&TYPENAME=ms:County&OUTPUTFORMAT=XMLSCHEMA">     <gml:boundedBy>     <gml:Box srsName="EPSG:27700">     <gml:coordinates>     294170.798523,90299.411200 655604.702242,588517.298279     </gml:coordinates>     </gml:Box>     </gml:boundedBy>     <gml:featureMember>     <ms:County fid="County.1">     <gml:boundedBy>     <gml:Box srsName="EPSG:27700">     <gml:coordinates>     320262.500400,378383.398000 361791.001700,422761.904601     </gml:coordinates>     </gml:Box>     </gml:boundedBy>     <ms:msGeometry>     <gml:MultiPolygon srsName="EPSG:27700">     <gml:polygonMember>     <gml:Polygon>     <gml:outerBoundaryIs>     <gml:LinearRing>     <gml:coordinates>     330996.596900,394516.299901 338224.102200,379434.497000 328536.798900,380475.199200     326182.199000,378383.398000 320262.500400,388522.196000 330996.596900,394516.299901     </gml:coordinates>     </gml:LinearRing>     </gml:outerBoundaryIs>     </gml:Polygon>     </gml:polygonMember>     <gml:polygonMember>     <gml:Polygon>     <gml:outerBoundaryIs>     <gml:LinearRing>     <gml:coordinates>     330871.099200,405417.000200 335362.999800,402517.699900 337285.798200,405717.803700     344062.801400,398931.696000 348231.000800,404146.903200 351662.604000,402905.498700     361791.001700,394518.799500 355128.200600,395322.399401 354161.500000,387836.602400     350551.103500,389713.904500 344678.600001,382585.295100 334583.800700,388194.599500     326895.800900,406184.799700 336963.201600,422761.904601 330871.099200,405417.000200     </gml:coordinates>     </gml:LinearRing>     </gml:outerBoundaryIs>     </gml:Polygon>     </gml:polygonMember>     </gml:MultiPolygon>     </ms:msGeometry>     <ms:OBJECTID>1    </ms:OBJECTID>     <ms:COUNTY_STR>Merseyside    </ms:COUNTY_STR>     <ms:PMS_REGION>7    </ms:PMS_REGION>     <ms:COUNTRY>ENGLAND    </ms:COUNTRY>     <ms:SNAC_GOR>North West    </ms:SNAC_GOR>     <ms:SHAPE_LEN>2.19563981295e+005    </ms:SHAPE_LEN>     <ms:EARegion>3    </ms:EARegion>     <ms:db_id>1    </ms:db_id>     <ms:updateID>0    </ms:updateID>     <ms:ET_ID>0    </ms:ET_ID>     <ms:ET_X>344814.69641    </ms:ET_X>     <ms:ET_Y>389730.08240    </ms:ET_Y>     <ms:Distance>0.00000000    </ms:Distance>     </ms:County>     </gml:featureMember>     <gml:featureMember>     <ms:County fid="County.2">     <gml:boundedBy>     <gml:Box srsName="EPSG:27700">     <gml:coordinates>     411883.903568,379684.895552 475119.498243,418312.098600     </gml:coordinates>     </gml:Box>     </gml:boundedBy>          These co-ordinates are correct for the data I'm trying to serve so I'm happy that things are appearing to be read ok. Panning and zooming in the map div on my test page will result in the content being altered based on a different bounding box, I just can't see anything. Am I missing some styling information somewhere? I've searched for solutions to this but nothing has worked. There's no proxy set up but my understanding is that for a completely local setup, a proxy is not required. Please can someone take a look at the setup I have and let me know if I've missed something obvious.


Azimuthal Equidistant - 40,000 Km range

I've been looking at making an azimuthal equidistant map. There are several web sites that offer to make one, based on my location, but they all have the same problem. They assume my antenna rotates 360 degrees - so they map the entire world, staarting from my location, and finishing at 20,000km (1/2 earth circumference). However, my antenna setup only covers 180 degrees. As such, I'm looking for a "half circle" shaped map, that ranges to 40,000 km. Does such a map exist? Is "azimuthal equidistant" simply the wrong search term, leading me astray?


Custom Symbols are too small in KMZ File from ArcMap

We are working with the local Police department on an iOS and Android app which will have an embedded Google map showing crime data from the past month. This map is generated via a .KMZ file exported from ArcGIS. We created a set of custom .EMF icons for this map and the Crime Analyst at the Police Department imported them fine, but they are extremely small in the KMZ layer that gets added to the Google map. It looks like the icons get saved to the .KMZ file as 16x16 .PNGs. When the Crime Analyst increases the size of the icons in ArcGIS, it makes them pixelated and blurry. I've also tried unarchiving the .KMZ file and inserting larger versions of the custom symbols into the folder and then zipping it up and converting it to a .KMZ file, but then the Google Map wouldn't display any of the custom symbols. Is there a simple way to retain the original size of the imported custom icons or to increase their size without pixelating them? We have the custom icons as .SVG, .PNG, .EMF, and as a custom icon font. Is it possible to use either the .SVG or the icon font as a custom marker in ArcGIS so that scaling will not be an issue? Let me know if you need more information or to see some of the files. I don't know much about the Arc suite of software and am just the app developer, but I can get information from the Crime Analyst. **EDIT** : Here are the answers to your questions.   * What tool are you using to export the kmz? - In arctoolbox, conversion tools, "layer to kml" tool   * What is the projection of your map document? - NAD_1983_StatePlane_South_Carolina_FIPS_3900_Feet Lambert_Conformal_Conic   * What is the unit of the projection? - Feet   * Imported them into what? - Select a generic symbol. Edit symbol in **Symbol Property Editor**. Symbol Type is **Picture Marker Symbol**. Then select the .EMF file.


How do I get ArcGIS tiles to display a single layer when multiple layers are available?

I am using Bing Maps v7 to display layers from an ArcGIS server. The server displays multiple layers as one, but I would like to only display a single layer. Suppose a map service has the layers: Topographic Info (0), Places Info (5), and Citations(10) I write code to get a tile like so:               var source = new Microsoft.Maps.TileSource({             uriConstructor: function( tile ) {                 return baseUrl + "/tile/" + tile.levelOfDetail                              + "/" + tile.y + "/" + tile.x;             }         });      Which would procdure a URL like this: http://server.arcgisonline.com/ArcGIS/rest/services/World_Topo_Map/MapServer/tile/12/1528/1218 Is it possible to have the generated tile only display information regarding a single layer, such as Topographic Info (0)? I tried variations of things like ?layers=0 in the query string, but nothing seems to work.


How to create QGIS theme map?

I need to make theme map using diagrams like in MapInfo, can I do it in QGIS? Or can I create a surface using attributes?


Azimuth and Distance, NumericalDigitize, Numerical Vertex Edit, XY tools and other plugins are not on QGIS Plugin Manager version 2.0.1

Will these plugins no longer be available in QGIS 2.0.1?


QGIS Cloud in QGIS 2.0 missing and 2.0 Crashing

Is QGIS CLoud still missing off QGIS 2.0. I am having a lot of issues with 2.0. on another matter I cannot load print templates without the system crashing. Are templates non compatible. When i am also loading projects into 2.0 from 1.8 some raster layer properties have reset to default values....eg greyscale back to colour.


QGIS 2.0 Plugins missing

At the risk of a series of downvotes but I have just installed QGIS 2.0 after using 1.7.4 and 1.8 previous. It looks great but there are a lot of plugins missing. I understand from the install that I may need to install them again for fear of system conflicts. However some really good pluggins are missing. The following are the main ones:   * Numerical Vertex Edit    * Easy custom Labeling    * Open Layers Plugin    * Open Streetmap    * QGIS Cloud    * Remove Empty Layers    * Spatial Query This list is not exhaustive but I am sure they are built in somewhere and can't find them. Thanks.


Colorify polygons in a pythonic way in Qgis

I want to addd a Layer and give the layer a color for the border and for the filling. With the code below i am possible to achieve a color for the border, but a did not found a solution to give a color for the filling. In addition, is it possible to change the width of the border line?               uri.setDataSource("shema_name", "layer_name", "geom","","layer_name_id")     lyr = QgsVectorLayer(uri.uri(), "layer_name_you_like", "postgres")          props = { 'color_border' : '255,165,0,125', 'style' : 'no', 'style_border' : 'solid' }     s = QgsFillSymbolV2.createSimple(props)     lyr.setRendererV2( QgsSingleSymbolRendererV2( s ) )          QgsMapLayerRegistry.instance().addMapLayer(lyr)


Csv To Shape Error

I have found a usefull tool to convert csv files into shape files with geotools: http : // docs.geotools.org/latest/userguide/tutorial/feature/csv2shp.html#create-a- shapefile-from-a-featurecollection But i have a problem running it, i get following error: Failed to execute goal org.codehaus.mojo:exec-maven-plugin:1.2.1:exec (default-cli) on project tutorial: Command execution failed. Process exited with an error: 1 (Exit value: 1) The error comes from this line:               final SimpleFeatureType TYPE = DataUtilities.createType(                     "Location",                   // <- the name for our feature type                     "location:Point:srid=4326," + // <- the geometry attribute: Point type                     "name:String"         // <- a String attribute             );      My POM file is this:               <project         xmlns="http://maven.apache.org/POM/4.0.0"         xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"         xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">         <modelVersion>4.0.0</modelVersion>         <groupId>org.geotools</groupId>         <artifactId>tutorial</artifactId>         <version>0.0.1-SNAPSHOT</version>         <packaging>jar</packaging>         <name>tutorial</name>         <url>http://maven.apache.org</url>         <properties>             <project.build.sourceEncoding>UTF-8</project.build.sourceEncoding>             <geotools.version>11-SNAPSHOT</geotools.version>         </properties>         <dependencies>             <dependency>                 <groupId>junit</groupId>                 <artifactId>junit</artifactId>                 <version>3.8.1</version>                 <scope>test</scope>             </dependency>             <dependency>                 <groupId>org.geotools</groupId>                 <artifactId>gt-shapefile</artifactId>                 <version>${geotools.version}</version>             </dependency>             <dependency>                 <groupId>org.geotools</groupId>                 <artifactId>gt-swing</artifactId>                 <version>${geotools.version}</version>             </dependency>             <dependency>                 <groupId>org.codehaus.mojo</groupId>                 <artifactId>exec-maven-plugin</artifactId>                 <version>1.2.1</version>             </dependency>         </dependencies>         <repositories>             <repository>                 <id>maven2-repository.dev.java.net</id>                 <name>Java.net repository</name>                 <url>http://download.java.net/maven/2</url>             </repository>             <repository>                 <id>osgeo</id>                 <name>Open Source Geospatial Foundation Repository</name>                 <url>http://download.osgeo.org/webdav/geotools/</url>             </repository>         <repository> <!--Add the snapshot repository here-->             <snapshots>                 <enabled>true</enabled>             </snapshots>             <id>opengeo</id>             <name>OpenGeo Maven Repository</name>             <url>http://repo.opengeo.org</url>         </repository>       </repositories>         </project>      I hope someone can help me!


How to register as versioned using SDELayer?

This is the normal process using ArcCatalog : ![enter image description here](http://i.stack.imgur.com/9Tiu4.png) I want to automate this process using sdeLayer or any tools - does anybody have an idea as to how to do it? **there is an application sdelayer (sdelayer.exe) in "C:\Program Files\ArcGIS\ArcSDE\pgexe\bin" i want to create a command line to register featurelayer**               sdelayer -o register -l <table,column> -e <entity_mask> -t <storage_type>     {[-C <row_id_column>[,{SDE|USER}[,<min_ID>]]]} [Spatial_Index]     [-E <{empty | xmin,ymin,xmax,ymax}>] [{-R <SRID> | [Spatial_Ref_Opts]}]     [-S <layer_description_str>] [-k <config_keyword>]     [-i {<service> | <port#> | <direct connection>}] [-s <server_name>]     [-D <database_name>] -u <DB_user_name> [-p <DB_user_password>] [-q]


Google map showing Overlapping Boundaries

I am working in process of building an online Geocoding and Reverse Geocoding Service .For That I have been preparing a database by creating a centroid point in a shapefile and sending the point to Google API and getting address polygons.After that i will merge similar address polygons. my problem is this: By reverse geocoding process,I took an example of Two locality polygons say : **_Amjikarai,Bala Vinayaga Nagar_** :the obtained output polygons of two area is: ![enter image description here](http://i.stack.imgur.com/vnbe9.jpg)                       where there is disjoints and Intersections .      So I thought to check the their respective boundaries in Google map.I was shocked that Google itself showing overlapping boundaries of the two areas: ![enter image description here](http://i.stack.imgur.com/H9XcV.jpg) So what is the problem or how can I rectify this problem?


Sum many netCDF rasters

For the first time, I am using NetCDF data files containing daily data about soil moisture at a global level. I want to be able to obtain monthly mean value of soil moisture for each grid cell. **There is one NetCDF file for each day, and each file contains four bands. I only need the first band which contains the data I want.** Using QGIS 1.8, I encounter the following problems: _1. Loading the rasters_ When I want to load the 31 rasters related to January (example), I can select the 31 files, but then, I have to select the band I am interested in for each layer, one by one. **Is there a way to load all of them, but automatically select the band I want?** ![enter image description here](http://i.stack.imgur.com/u9EAc.png) _2. Converting the layers_ The data comes in a floating number format. I detected some issues when it comes to sum values. However, when I convert the layer to a .tiff raster, it seems to transform the value into an integer. This is what I want, but GDAL doesn't want to batch convert all the files into .tiff, either because of the file names of because of the multiband issue. **Is there a way to batch convert NetCDF files into .tiff?** ![enter image description here](http://i.stack.imgur.com/CquXN.png) _3. Batch summing the layers_ Finally, even if I have the layer in .tiff format, I wonder if there is an easier way to sum all the values than by using RasterCalculator. **Is there a tool than batch sums up rasters within a folder?** Also, I noticed that I cannot use RasterCalculator directly on the NetCDF layers because of the file names. If the solution is on this trail, **is there a way to batch change the name of a layer within a NetCDF file?** That's a lot of questions, so many thanks in advance for your help! Damien