Geoserver does not see new layers

I have an existing GeoServer instance running (version 2.1.3). This GeoServer has multiple Workspaces and Stores serving a number of layers as WFS to an OpenLayers map. I want to update my Stores every quarter with the most recent layers. The process goes like this: 1) Delete the GeoServer Store (for example Store: "3rdQ2012") under Workspace "A" that contains the previous quarters shapefile. 2) Create a new GeoServer Store (for example Store: "4thQ2012") under Workspace "A" for the new layer. 3) Edit my OpenLayers map project to point to the new layer. For some reason when I do this GeoServer is not serving the new layer. When I look at the GeoServer log file I see this message for the affected layer: "2012-12-21 16:15:33,965 WARN [geotools.gml] - Could find type for _nameOfLayer_ in the schema, generating type from feature." When I go to the GeoServer "Layer Preview" page I can successfully see and select "OpenLayers" under the common formats column GeoServer successfully opens a map window and displays the layer. But, when I select GML under the common formats column GeoServer opens a new window and displays this : " This XML file does not appear to have any style information associated with it. The document tree is shown below. http://schemas.opengis.net/wfs/1.0.0/OGC-exception.xsd"> java.lang.NullPointerException null " It appears that GeoServer is not seeing my new layer or perhaps something is not updating correctly after updating my stores and layers. I have tried restarting my GeoServer service as well as doing restarts on the server. I have also deleted the workspace, stores, and layers multiple times in an attempt to force GeoServer to update but all to no avail. Any suggestions would be appreciated. Thanks!


display a selected WMS Info

I am using the WMSGetFeatureInfo and I am also obtaining the Info of the base raster and i do not need that. Is it possible to obtain the FeatureInfo of only selected WMS features? How can I modify this code to accommodate only selected WMS layer?               info = new OpenLayers.Control.WMSGetFeatureInfo({                 url : 'http://niles.iplantcollaborative.org:8080/geoserver/wms',                 title : 'Identify features by clicking',                 queryVisible : true,                 eventListeners : {                 getfeatureinfo : function(event) {                 map.addPopup(new OpenLayers.Popup.FramedCloud("chicken", map.getLonLatFromPixel(event.xy), null, event.text, null, true));                 }                 }                 });                 map.addControl(info);                 info.activate();      **UPDATE 2.0** I have added the code,               gfLayers = [];                 var lyr = map.layers;                      for(var a=0; a < lyr.length; a++){                     if(lyr[a].CLASS_NAME == "OpenLayers.Layer.WMS" && lyr[a].visibility == true){                         if(lyr[a].isBaseLayer == false){                             gfLayers.push(lyr[a].name);                         }                          }                  };                      info = {                     hover : new OpenLayers.Control.WMSGetFeatureInfo({                         url : 'http://niles.iplantcollaborative.org:8080/geoserver/wms',                         title : 'Identify features by clicking',                         //layers : gfLayers,                         layers : [pima, carib],                         queryVisible : true,                         hover: true,                         eventListeners : {                             getfeatureinfo : function(event) {                                 deletePopUp();                                 popup = new OpenLayers.Popup.FramedCloud("popInfo", map.getLonLatFromPixel(event.xy), null, event.text, null, true);                                 map.addPopup(popup);                             }                         }                     })                 };                 map.addControl(info.hover);                 info.hover.activate();      Now I am getting Empty PopUps. not sure why?


TF.Net versus the NetTopologySuite

In a recent post, a responder made reference to the Topology Framework .Net. I'd never heard of it before, having gently used the NetTopolgySuite (and DotSpatial) before with satisfactory results. Why would one use TF.Net over NetTopologySuite, or vice versa? Why does TF.Net exist if NetTopologySuite has been around for quite a while?


Reprojection - Learning the ins and outs

I would like to expand my knowledge coordinate systems, projections and reprojecting data, with a view of writing a reprojection engine. Could you provide me with some resources (books, websites, scrolls) on the maths behind coordinate systems and reprojecting data?


How do I Capture GDAL projection and create matching output (OGR or GDAL) shape files?

Working in C++, using GDAL and OGR: Depending on what the user has supplied for image input, I need to output point (or other) shape files in the same projection as the input raster images. The output shape files must be readable by ESRI as defined to be in the same projection as the input image. We are working with many input images and output shape files, so the projection process must be automated. How do I best transfer the input GDAL projection information to create a corresponding output shape file of the same projection? Do I have to go to OGR and use well-known text protocols to output the shape file? Can I instead use the GDAL projection information to write out the shape file using GDAL somehow?


Spring and OpenLayers conflicts on label rendering?

I need to use Jquery to handle other elements on the webpage. I tried to put labels on a vector layer with the following style:                 var vector_style = new OpenLayers.Style({                 'fillColor': 'red',                 'fillOpacity': .8,                 'strokeColor': '#aaee77',                 'strokeWidth': 3,                 'pointRadius': 8,                 'label': "${name}",                 'labelAlign': "cc",                 });      But the labels do not show. I am just wondering if there is a conflict on this one "${name}". I tested it again with a simple page without Jquery. The labels do show. Please let me how to fix this. Thanks. **Update:** I just figured out the conflict may come from Spring or JSTL. I use Firebug tool to check the javaScript source code of the style. It appears as this:                   var vector_style = new OpenLayers.Style({                 fillColor: 'red',                 fillOpacity: .8,                 strokeColor: '#aaee77',                 strokeWidth: 3,                 pointRadius: 8,                 label: "",                 labelAlign: "cc",                 });      I returned an variable from the server side:                   model.addAttribute("name", "test");      then all features show the label "test". It approved that ${name} is recognized as a spring variable.


Flex Geoprocessing Widget output parameter invalid

I want to publish a ArcToolbox with a script on a Flex application as a widget. I have ArcGIS server from which I will publish it as a geoprocessing service. The tool has an input parameter which is a text file and the output parameter has to be a folder which will route the user to choose a folder on computer. When I try to publish this tool, it says that the output parameter - folder is invalid. Is it possible that the folder parameter is not supported by flex output data? Any way to do that?


How to find a point on given line that its distance from another point is given also?

I have this input:   1. (lat1,lon1), (lat2,lon2) - Two coordinates that create the line   2. (lat0, lat0) - Point on this line    3. d - distance I want to find a point (lat, lon), where   1. (lat, lon) is on the line (mentioned above)    2. The distance between (lat, lon) and (lat0, lon0) is d The direction of (lat0, lon0)-to-(lat, lon) is not importent right now. [and of course, the lat/lon is given as geographic coordinates in degrees]. Can anyone knows any algorithm for that? Many thanks!


pyQGIS : Way in order to know modifications?

I 'm writing a pyqgis script. I succeed in accessing gis vector data (lists) from my remote server (thanks to urllib) . I succed in displaying geometries and attributs. 1/ Now, i would like manage layer 's editions and modifications. I must find a way to know the polygons that were been modified , deleted or added. Could you throw light for me? 2/ In a 2nd time i would like send these editing, these modifications toward the remote server. Could you tell me how to do that please? in order to be more explicit i paste the usefull github link: github.com/picardie-nature/clicnat-qgis/blob/master/plugin.py Thanks


convert my large svg image into tiles

This is probably going to annoy the more experienced map people here, but I'm sort of new to all this, or at least I don't have a lot of experience. I know this is a bit strange and not ideal, but I have a huge SVG file of geographical shapes and regions. I want to be able to plot points and draw shapes in Openlayers. The problem is that, of course, this file is over 50 MB, so I can't exactly expect users to load the whole thing then navigate around.   I have a less detailed map for when users are zoomed out, so I'd like to use this map as they zoom in, but then tiling would of course be best. I don't need much detail beyond what I have so using just these two maps as starting points is fine. Like I said, I was planning on using OpenLayers since I've used it in the past but I'm not tied to it. How is the best way for me to proceed?   Thanks.


PostGIS query with input from point layer

I have a table of rivers in PostGIS, each river linking to its outflow. I also have a QGIS layer of some selected sources, and would like to select the entire river downwards from each source (i.e. find the closest segment to each point, then recurse through the outflow links until hitting the ocean). Is there a way to do this interactively? I'd like to just update the point layer and have the query use these points. I don't think QGIS offers any help with that (parametrized queries?), but I could store the point layer as a shapefile, load that as a temporary table and use that in my query.   * is there a way to load a shapefile in a subquery? I don't want to use `shp2pgsql` manually, but load the current version of the file each time the query is run.   * can I maybe somehow trigger actions in QGIS when a layer has been edited? i.e. update the river layer when the point layer has been saved with the new query output, or, alternatively, run a script that loads the point data, runs the query, dumps its data into a shapefile, loads that in QGIS and refreshes the view?


Conversion from .shp to.swf?

Does anyone know which free software can I use to convert a .shp file to a .swf?! I have created a shapefile on QGIS which I would project on STAT PLANET PLUS. In addition, I would like to display some other maps when the user selects a specific indicator. In order to be able to project additional maps, the files have to be in .swf format, so my shapfiles have to be converted from .shp to .swf. I think maybe the trial version of Adobe Illustrator can make the conversion (but I dont know how), any other ideas?! Any help would be appreciated. Thanks.


Population density map

I have California population density data from from the U.S. Census in an Excel file. The spatial information is provided in ZIP Code Tabulation Areas (ZCTAs). How do I plot a population density map with that data using R? This is what I have been doing: dataset               library(maptools)     library(RColorBrewer)     library(classInt)     library(maps)          ## set the working directory.     setwd("../Downloads")          ## load the file     pop <- read.table("mre.txt", header=TRUE, row.names=1, sep=",")          #select color palette and the number colors (levels of income) to represent on the map     colors <- brewer.pal(9, "YlOrRd")          #set breaks for the 9 colors      brks<<-classIntervals(pop$POP10, n=9, style="quantile")     brks<- brks$brks          #plot the map     plot(pop, col=colors[findInterval(pop$POP10, brks,all.inside=TRUE)], axes=F)      I also tried this other way:               CaliforniaZTACs<-readShapePoly("[tl_2010_06_zcta510.shp][2]")     plot(CaliforniaZTACs)     CaliforniaZTACs      but could not associate a population density value from my dataset to the location ID of the tl_2010-06-zcta510.shp file.


From where I can get 4-band satellite imagery of Sweden?

I'm looking for good quality, fine resolution 4-band satellite photos of Sweden for academic use (and I'm willing to pay for it). What I'm after is as detailed as possible photos of certain part of Sweden where I would like to classify vegetation with commonly used indices. In order to do so I need NIR band, which is missing in the data provided by Swedish authorities. They instead offer readymade vegetation maps, which are a bit outdated for my purposes. Photos should be taken after 2005, preferably for a couple of years in different seasons. In addition I'm looking for good quality weather data for Europe and Africa, as accurate as possible in order to evaluate the conditions birds face during migration from wintering grounds. Could anyone help?


How to export offsets from Trimble Rangefinder?

I already asked this question at http://www.rpls.com/forums/viewtopic/8/1276880, but didn't get any answer so far. Maybe here someone can help... I have a question about the use of Trimble Products. I didn't find a specific forum for that, so I hope that I might get some help here. If I'm in the wrong forum, maybe anybody could point out where I could find help? We have the following devices: GPS Pathfinder Pro, Handheld with TerraSync and a TruPulse 360b Laser Rangefinder. In a (dense) forest, we want to map all trees with 13 hectares (~3000 trees). Some gaps in this forest allow measurements with the dGPS, but mostly there is no reception. By measuring distances, azimuths and slopes between neighboring trees, we could calculate coordinates for all trees, based on the few GPS measurements we take in the gaps. For this, I can create a feature (without position --> no reception) for each measurement and take an offset measurement from one tree to another with the Rangefinder. I can than export the data file to GPS Pathfinder Office software (V4.2). In the software, I can see the offset for each feature (distance, azimuth, slope) by clicking on it. But then there's a problem: In order to use the offset for the calculations, I'd have to export it. In this isn't working. It is possible to export a lot of attributes for each feature, but apparently not the offset. So here are my questions: \- Is it somehow possible to export the offset measurements with GPS Pathfinder Office (V4.2)? \- If not: Is this possible with the newest version of this software? \- If not: Could any other software accomplish this? Any help would be greatly appreciated! I really wouldn't like to copy thousands of measurements by hand.


How to convert gps receiver native format to google lat long or XY format

I have received some GPS co-ordinates in the following format like - "17267690N 078279340E".. How to convert the co-ordinates to show it on Google map.


How show some fields only when using WMSGetFeatureInfo?

How to show some important field only when using WMSGetFeatureInfo.I use this example. When I click on maps it show all the fields of attribute table.What part should I modify in above mentioned code. Pls help me? Any idea wellcome


Problem uninstalling License Manager

If i try to unnistall current version of License Manager (ver 9.3.1) i get the following message: ![alt text](http://i.stack.imgur.com/PptPF.jpg) This a problem because now i can't install the new version of License Manager for ArcGis 10. Do you know the solution?


SDE views not showing values in joined table

I created an SDE view between a point feature class and a table (stored in SDE) - the join is done on a text field called RFRNC. I changed the view so that it is now a LEFT OUTER JOIN in line with this doc: http://www.jamesrichards.com/post/2009/07/16/How-To-Create-an-ArcSDE-Spatial- View-With-an-Outer-Join.aspx Problem is - when create a row in my flat table in the RFRNC column with a value I know exists in the featureclass - it never shows up or updates in the view - what is going on here? Any ideas greatly appreciated.


Use lat / long to create polyline that can be passed into DrawPolyline

public void DrawPolyline(ESRI.ArcGIS.Carto.IActiveView activeView)     {       if(activeView == null) return;       ESRI.ArcGIS.Display.IScreenDisplay screenDisplay = activeView.ScreenDisplay;            screenDisplay.StartDrawing(screenDisplay.hDC, (System.Int16)ESRI.ArcGIS.Display.esriScreenCache.esriNoScreenCache); // Explicit Cast       ESRI.ArcGIS.Display.IRgbColor rgbColor = new ESRI.ArcGIS.Display.RgbColorClass();       rgbColor.Red = 255;            ESRI.ArcGIS.Display.IColor color = rgbColor; // Implicit Cast       ESRI.ArcGIS.Display.ISimpleLineSymbol simpleLineSymbol = new ESRI.ArcGIS.Display.SimpleLineSymbolClass();       simpleLineSymbol.Color = color;            ESRI.ArcGIS.Display.ISymbol symbol = (ESRI.ArcGIS.Display.ISymbol)simpleLineSymbol; // Explicit Cast       ESRI.ArcGIS.Display.IRubberBand rubberBand = new ESRI.ArcGIS.Display.RubberLineClass();       ESRI.ArcGIS.Geometry.IGeometry geometry = rubberBand.TrackNew(screenDisplay, symbol);       screenDisplay.SetSymbol(symbol);       screenDisplay.DrawPolyline(geometry); //I WANT TO CREATE THIS FROM LAT / LONGs!       screenDisplay.FinishDrawing();     }      The above code snippet works but the problem is I want to create the IGeometry from lat \ longs not from clicking the screen. How can I create an polyline IGeometry object for the drawPolyline function using standard lat \ long coordinates?


GRASS : r.resamp.rst error : "Not enough disk space--cannot write files"

Today, I am confronted to a difficulty I cannot manage. When I try to use the tool "r.resamp.rst", it gives me as answer:               *r.resamp.rst input=map@PERMANENT ew_res=11 ns_res=11 pcurv=curv                Processing all selected output files will require     261320 bytes of disk space for temp files.     Temporarily changing the region to desired resolution ...     Changing back to the original region ...     Percent complete:     Not enough disk space--cannot write files     Not enough disk space--cannot write files     Not enough disk space--cannot write files     Not enough disk space--cannot write files     .....     interpolate() failed     dnorm in mainc after grid before out1= 12.961481     ERROR: split_and_interpolate() failed*      I have 14GB left on my hard disk so I don't understand why they are saying me that :/ I tried to find solution, and I found some answers (e.g. here : http://lists.osgeo.org/pipermail/grass-dev/2005-April/018152.html) but I don't understand ^^ They speak about a "ulimit -f" function, but where should I write it ? It doesn't work in GRASS. I am on Windows, maybe is that the problem. SO in conclusion, I am more than lost so if someone could help me it would be perfect! Thank you, b.raoul


ArcMap extension 9.3 vs 10

I have updated our ArcMap extension to work with ArcMap 10 - mainly just change reference from ESRI.ArcGIS.ADF to ESRI.ArcGIS.ADF.Local and set value of 'Specific Version' property for all references to false, and also use the new ESRIRegAsm.exe to register it. So I am wondering now if I can build a single extension that will work with both 9.3.1 and 10.0 or do we need separate binaries for each product?


Perform a spatial join in C# with ArcObjects

I feel like this has to be on the internet somewhere, but I haven't been able to find it. **I'm attempting** to perform a spatial join on a polygon and point layer. My goal is to push the polygon's ID to each point that it contains. In PostGIS, it might look close to this: `SELECT a.*, b.* FROM a, b WHERE ST_CONTAINS(a, b)`. **I have:**   1. A polygon IFeatureClass, from a geodatabase   2. A point IFeatureClass, from a geodatabase   3. The name of the polygon field I would like to turn this into a point IFeatureClass with an extra field. **My question:** How do I do this? Bonus points for links to relevant documentation.


Which training courses in Geodetics can you recommend?

Can anyone recommend any training courses in geodesy which would cover the following objectives:   * Case studies on coordinate issues   * Overview of geodetic datums, coordinate reference systems and transformations   * Use of Global navigation satellite systems ie. GPS    * Methods in projection methods and methods in measuring the Earths surface   * Vertical datums, geoidal models, vertical CRS and their transformations   * Orthorectification, georeferencing and other data referencing methods Thanks Rob


QGIS, Postgis: Is it possible to create a new postgis layer?

Situation: I have a postgis vector lay that has lots of line segments (the features). I want to create a set of new layers each of which contains a subset of the features of the original layer and have these layers also stored back in postgis. Ideally with the actual features stored just once. The only way I can see to create a new layer as opposed to load one from some source is via Layer->new which gives an option of shapefile or spatialite. Choosing shape file I first tried to specify gid and fid attributes for the new layer but SPIT complains about duplicate field gid when I try to load it into postgis. If I do not create any attributes for the new layer SPIT happily loads it, presumably duplicating all the features. I am guessing this is a 'feature' rather than an oversight because if it were possible to share features between layers then then what would happen if a feature were split? If things were to work properly then all layers with referenced that feature would need updating with the new gid of the new portion. Is this interpretation correct? Russell


Can county level info be brought into a map?

I have a csv file of Clovis projectile points that were recorded only with the county in which they were found, no coordinates at all, no shapefiles. Is there a way to bring these into a state map with counties to show county level finds? * * * I have a geocoded county outline set for VA. Imported with no problems. I have geocoded site locations for VA. Imported with no problems. Now I have 920 Clovis spear points that were recorded to the county level only and which are not geocoded in anyway, except that the county is listed. I asked the question here and was directed to MMGIS and got the same error message that 2 other folks got so that didn't work. Then tried to use the vector join. That brought in only some of the info from the Clovis csv file I wanted and then only up to the number of counties in VA (105). So how do I do this. The ArcGIS folks can do it but can QGIS do so and if so, can someone give this newbie step by step instructions as I have tried dozens of combinations to no effect that is useful. Thanks


How do I create a drop down box / combo box of bookmarks from an ArcGIS Online webmap?

Without using the clunky ArcGIS bookmark dijit, how do I create a drop down box of bookmarks from an ArcGIS Online webmap (or JSON) that zoom to the extents when clicked?               //FYI, JSON format     "bookmarks": [         {           "extent": {             "spatialReference": {"wkid": 102100},             "xmax": -1.2251063434141435E7,             "xmin": -1.2561092020865962E7,             "ymax": 4977986.7113767695,             "ymin": 4811659.737828305           },           "name": "Utah County"         },         {           "extent": {             "spatialReference": {"wkid": 102100},             "xmax": -1.2615515185004937E7,             "xmin": -1.3235572358454255E7,             "ymax": 4310538.580390846,             "ymin": 3977884.633293776           },           "name": "San Bernardino County"         }     ]


How to list all of the contents of a geodatabase with arcpy?

This seems like a simple problem, but I have not found a way to accomplish it. I am writing a script that needs to record the paths of everything inside of a geodatabase to a text file, and later copy those files to a new location and zip them up. I have been able to use `arcpy.ListFeatureClasses` to grab feature classes, and `arcpy.ListRasters`, and `arcpy.ListDatasets`, but how do I obtain a list of, say, Address Locators in my gdb? I haven't found a way to do it with arcpy. **There should be a simple function that lists all of the contents of a geodatabase. Is there one I'm missing?** Or is there any other way to accomplish this in python?


How to install GDAL with TIFF support on Linux?

I am trying to install GDAL on a linux machine and tested some sample commands like `gdalinfo --version` This worked fine, but when I started working with images like `.tif` I got the following error .               Warning 1: TIFFReadDirectoryCheckOrder:Invalid TIFF directory; tags are not sorted in ascending order     ERROR 1: MissingRequired:TIFF directory is missing required "ImageLength" field     gdalinfo failed - unable to open '1.tif'.      I reinstalled `libtiff` but I am still getting the same error. Iam working on CentOS release 5.7 (Final) .Is there any one who worked on this. Any help is appreciated.


How to compile liblwgeom independently?

For using liblwgeom function in spatialite, I need to compile the library. I found that it is included in the PostGIS code. My question is, how to compile the library without compiling the whole PostGIS?


How is a "Map Package" created programmatically?

Given a loaded map document (mxd) how do I find the geodatabase objects related to the contained layers? I want to copy them with the map document so the map can be moved onto each user's local machine from a network share. Is there an API to do this? Can someone point to source code that does this? To clarify (after all the great feedback): How is a "Map Package" created programmatically using the ArcEngine API? Is this even possible?


Retrieve tile path of superoverlay in Google Earth Plugin using Javascript

I've created a superoverlay using MapTiler and am displaying that in a webpage with the Google Earth plugin. I'd like to identify the current tile (level, x, y) being displayed with a mousemove event. I'm able to get the current lat/lon using getGlobe(), e.g.               google.earth.addEventListener(ge.getGlobe(), 'mousemove', function(e) {         var lat = e.getLatitude();         var lon = e.getLongitude();     }      The biggest hurdle seems to be getting the correct superoverlay level based on camera altitude as the transition between levels is dictated by the minLodPixels parameter in the KML files. I'm not sure how best to determine when a tile is, say, 128 pixels wide on the screen, especially when this will be based on how large the map window div is. Once the level is determined, it should be trivial to figure out the X, Y offsets of the tile. I'm also assuming that the minLodPixels/maxLodPixels parameters will stay consistent across all tiles in the superoverlay (MapTiler has these at 128/-1). Any thoughts? **Update 2012.02.22** I've delved a bit further into this and found some good documentation on minLod in the KML documentation. From http://code.google.com/apis/kml/documentation/kmlreference.html#region, > Google Earth calculates the size of the Region when projected onto screen > space. Then it computes the square root of the Region's area (if, for > example, the Region is square and the viewpoint is directly above the > Region, and the Region is not tilted, this measurement is equal to the width > of the projected Region). If this measurement falls within the limits > defined by and (and if the is in view), the Region is active. So if I'm able to find out the screen real estate of each hierarchical tile at a particular LookAt, I should be able to find the level I'm at based on what tile minimally exceeds the minLod threshold. As a test, I've read the SuperOverlay's root KML ('doc.kml') and retrieved the coarsest tile bounds:               var INFO = '/path/to/doc.kml';     $.get(INFO, function(data) {         // Get the corners of the coarsest resolution         var box = $(data).find("LatLonAltBox");         north = parseFloat($(box).find("north").text());         south = parseFloat($(box).find("south").text());         east = parseFloat($(box).find("east").text());         west = parseFloat($(box).find("west").text());     });      From there, I calculate area on the mousemove event of GEWindow.               google.earth.addEventListener(ge.getWindow(), 'mousemove', function(e) {              // Track the screen coordinates of the coarsest tile         var view = ge.getView();         var alt_mode = ge.ALTITUDE_RELATIVE_TO_GROUND;         var ne = view.project(north, east, 0.0, alt_mode);         var se = view.project(south, east, 0.0, alt_mode);         var sw = view.project(south, west, 0.0, alt_mode);         var nw = view.project(north, west, 0.0, alt_mode);              // Get the square-root of area (using equation 3 from          // http://mathworld.wolfram.com/Quadrilateral.html)         var qX = (se.getX() - nw.getX());         var qY = (se.getY() - nw.getY());         var pX = (sw.getX() - ne.getX());         var pY = (sw.getY() - ne.getY());         var area = Math.sqrt(0.5 * Math.abs(qX * pY - pX * qY));     });      When I track the tile's visibility in GE though, it disappears at a value higher than 128 (the minLod) when viewed without tilt, and appears at a value lower than 128 when viewed with tilt. I don't have enough rep to post pictures, but I've posted example screenshots here. In these pictures, the white polyline shows the tile boundary. I'm obviously still doing something incorrectly.


How do I solve a "No PROJ.4 translation" error with gdal?

I've installed `gdal` using the `brew` utility on OS X so the `proj` and `geos` dependencies were also installed. I'm now trying to warp a GeoTIFF so that I can import it into Google Earth. I am getting a `proj` error, and cannot figure out how to solve it. I use the following command:               gdalwarp -t_srs "EPSG:4326" in.tif out.tif      and get the following error message:               Processing input file in.tif.      ERROR 1: No PROJ.4 translation for source SRS, coordinate transformation initialization has failed.     0...10...20...30...40...50...60...70...80...90...100 - done.      There is very little on Google about this error. I have set my `GDAL_DATA` and `PROJ_LIB` environment variables to `/usr/local/shar/gdal` and `/usr/local/share/proj` respectively. I do not think this is a problem with the input file since it appears to have a coordinate system defined already, checked with `gdalinfo`:               Driver: GTiff/GeoTIFF     Files: in.tif            in.tfw     Size is 341, 558     Coordinate System is:     LOCAL_CS["WGS 84 / Antarctic Polar Stereographic",         UNIT["metre",1,             AUTHORITY["EPSG","9001"]]]     Origin = (-3143782.158253818750381,3005869.481321755331010)     Pixel Size = (6268.758042377240599,-6268.758042377240599)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=BAND     Corner Coordinates:     Upper Left  (-3143782.158, 3005869.481)      Lower Left  (-3143782.158, -492097.506)      Upper Right (-1006135.666, 3005869.481)      Lower Right (-1006135.666, -492097.506)      Center      (-2074958.912, 1256885.987)


How to visually sort the order of Points with QGIS, GDAL, ...?

The following Points layer is the basis of a rendering in a Android application. ![Points layer in QGIS](http://i.stack.imgur.com/qA99a.png) Android allows to draw a shape from the provided **points**. I **cannot** use LineStrings. ![Points layer in Android](http://i.stack.imgur.com/kxk2J.png) However, the order of points exported from QGIS seems to be odd. Is there any way to **visually sort** the order of points in QGIS or using the GDAL tools? The **order of feature ids** in the final export format (such as CSV) needs to be kind of a **roundtrip** as visualized in the following picture. If there are multiple paths or a "hole" it is perfectly okay to have multiple point layers. ![Points layer visualizing the roundtrip](http://i.stack.imgur.com/qVo0p.png) ## Export experiments I tried different steps to get hold of the latitude/longitude coordinates:   1. QGIS: Export points layer as CSV, add layer create options: `GEOMETRY=AS_YX`   2. Convert into Android specific format JSON format using regular expressions. This has the disadvantage that the feature id is still included. The order of points however is **not** as desired.   1. `ogr2ogr -f CSV output.csv input.shp -lco GEOMETRY=AS_YX`   2. Convert into Android specific format JSON format using regular expressions. Here the feature id is already omitted. The order of points however is **not** as desired.   1. QGIS: Export points layer as GML   2. `ogr2ogr -f "CSV" -select "" output.csv input.gml -lco "GEOMETRY=AS_YX"`   3. Convert into Android specific format JSON format using regular expressions. The order of points however is **not** as desired. The order of points is **different** to what the former exports/conversions output. * * * The original layer is also available as a LineString layer.


ArcMap 10.1 - UTM Zone Display Mismatch

I've recently noticed that the UTM zone displayed in the Identify box and the "Goto XY" box display a different UTM zone than the one I'm in. I'm working in WGS84 NAD83, Zone 13N but when I use the Identify tool and look at the location of the point it shows the location as '13S X,Y'. Also, if I type in a set of coordinates manually using the "Go To XY" tool, I type it in as 13N X Y, and hit enter, after it flashes the location the zone I had input is changed to 13S. _Can someone explain why the zone is being displayed as 13S and not 13N?_ I'm assuming this is one of those weird geometric things that I never picked up in school. Hopefully someone can enlighten me! ![enter image description here](http://i.stack.imgur.com/wTfML.png)


Stripping spatial reference

I am creating an exercise for remote sensing students and would like to take an existing raster and remove the spatial reference from the file so that the students can practice using another image to orthorectify the first. I am using ArcGIS 10. Any hints?


How can I flatten a multi-band raster to a single-band raster in ArcMap 10.2?

I need to convert a multi-band PNG image into a single-band raster, for analysis and visualisation. How can I do this in ArcGIS 10.2?


url completion in cartoDB info windows

By far the most of my data is taken up a column of long strings pointing to the source of the data for each point. I need to keep that. However, it would be neat if CartoDB could use this information (i.e. build the URL) without me providing the full URL, only the relevant ending (and save on the size of the table by not repeating the common part as a long string). Is this doable in their info window?


Google teardrop patented

We want to indicate a point on a map and like the idea of using a teardrop but it seems that google has the teardrop idea patented (http://www.google.com/patents/USD645052 http://www.google.com/patents/USD661704 http://www.google.com/patents/USD620950) Does this mean that even if we create our own, unique teardrop shape that we would be violating the patents that google has on the tear drop shape? Obviously it wouldn't look like googles, but it would look similiar since it would be a teardrop shape. Any other ideas on a good point marker? We are using Esri's tools, not google maps. Thanks!


Is it possible to replaceDataSource of layer with dataset stored in_memory using ArcPy?

I use in_memory workspaces to speed up some Python script tools that get published as Geoprocessing Services significantly. The tools open MXDs and perform, among other things, an ExportToPDF. In the MXD there is a layer that references a dataset stored on disk. At runtime a new dataset gets created with the same schema as that above, so that it can replace it as the data source but, given a choice, I would like that new dataset to ONLY be written to the in_memory workspace. My testing suggests that using in_memory workspaces like this is not supported by either Layer.replaceDataSource or Layer.findAndReplaceWorkspacePath. Is there a way to replace the data source of a layer with a feature class stored in the in_memory workspace so that it can be displayed in a map that gets exported to PDF before the in_memory workspace is deleted? It looks like this would possible using ArcObjects and an InMemoryWorkspaceFactory but here I need an ArcPy solution. If there is no way to do this using ArcPy at ArcGIS 10.0 (and nothing coming at 10.1) then what I will do is to propose as an ArcGIS Idea that Layer.replaceDataSource be enhanced to support a workspace_type=INMEMORY_WORKSPACE option. Note: I have also posted this to the ArcGIS Discussion Forums but have had no replies there so far. An **ArcGIS Idea** has also been submitted for this to be enhanced.


How to run multiple iterations of a model?

I am working with a model that is build in modelbuilder and I have to test with with over 150 different use cases. Is there a way I could upload a table with use cases as it attribute and then have the model run one time on all the 150 cases, one by one. For instance, have it run the first time one case #1 and after is has finished run again at case #2 and so on. I know I have to use iteration but unfortunately it only runs ones. Below is an image of the model: ![http://i.imgur.com/dH7tuhA.png](http://i.stack.imgur.com/DWHuB.png) Thank you E


How to create FeatureLayer without Legend in Table Of Contents (TOC) using ArcGIS Engine?

So what I am trying to accomplish is having a IGroupLayer that contains an IFeatureLayer and in the table of contents (TOC) there would be no + next to the feature layer (i.e. there is no legend). So far what I tried to do is this, which is not working:               ILegendInfo legendInfo = geoFeatureLayer.Renderer as ILegendInfo;     ILegendGroup legendGroup = legendInfo.get_LegendGroup(0);     legendGroup.Visible = false;


What is gid in OpenStreetMap data model?

I'm using an old code written by an ex-fellow and I wonder what is GID in OSM? The gid is basically used everywhere. The project uses PostGIS and PgRouting. Any pointers to that? I know this is a naive question but I need a start. Thanks. As an example, this method comes from this interface.               /**          * @param gid the gid of the openstreetmap entity we want to retrieve          * @return the OpenstreetMap entity or null if not found          */         public OpenStreetMap getByGid(final Long gid) ;


How to Install ArcGIS 10.2

I found out that ArcGIS 10.2 was released yesterday, which is great since I am looking forward to it. Although so far I haven't been able to find any information on how the update will be installed. I was hoping that it would be more of a patch update, but alas nothing to be found here. So it's probably a full uninstall and re-install similar to 10.0 to 10.1. Does anyone have any concrete information (with s source) on this?


How to copy an SDE geodatabase to a file geodatabase using Python?

I would like to make an exact copy (domains, feature datasets, feature classes, etc.) of an SDE database to a file geodatabase. I have tried several possibilities, including:   1. using the Copy (Data Management) process   2. creating a new GDB and manually copying each feature dataset from the SDE   3. exporting an xml workspace document from the SDE and importing it into the GDB The `Copy_management` process does not seem like it would work for copying an SDE to a GDB, since the input and output data types must match. The process of importing each feature dataset into a new GDB could probably be automated using Copy_management as well by iterating through each feature dataset, though it seems this could cause problems of an incomplete copy if there was an error with one of the processes. Exporting and importing xml workspaces seems to work, though this process creates incredibly large files when the process is used on large geodatabases. Is there a more straightforward way to copy the contents and schema of an SDE to a GDB than the ways mentioned, in a way that can be automated? If not, are there any reasons that the above possibilities should not be used in this process?


Why is ArcGIS for Desktop not joining records when values have dot (.) in text string?

ArcGIS is not joining records from a table that have a dot (.) in the text string. Has anyone else come across this problem? Is it a known issue?


How can I convert/transform coordinates from ETRS89 to WGS84 Utm 36N in QGis

I am a student and just getting to know qgis. At the momment I have to combine data which had been measured in ETRS89 with data in WGS84 UTM 36N. I created several vektorlayers for the data in ETRS89. Then I tried to transform the data into wgs84 by using "on the fly transformation" but it seemed to me that there was no change visible. The coordinates for one of our fixpoints are in WGS84 279000.000/4207000.000/1616.875 but Qgis shows me the coordinates 542000.00/4206000.00 (the height I cant even see). How can I transform the coordinates? I dont mean a change in projektion I want to have a proper translation/conversion of the coordinates? Could somebody please explain it step by step to me? Thank you very much! Chrissi Fist of all let me thank you for your quick reply. The more western points are the correct ones in the open street map. I think we still not have solved the problem. May be I should explain it a little bit better... I have got some data of turkish students. They contain meassurementpoints which had been aquired with a totalstation. We also meassurred with a total station in the same area but we meassured in WGS84 UTM 36N relying on fixpoints wich had been set witch GPS. If you now look at the data, there has to be a difference because the coordinates of our tukish colleaugues begin with 542.../420... and ours 279../420... and also the height is divergent. But im not quite shure whether the information of the turkish student was correckt that it is ETRS. I really dont know how I can bring this data together


open QGIS project with parameters

I am looking at a piece of work to deploy QGIS as the standard within an organisation. There is an existing piece of software that I need to link to QGIS so that when a button is pressed within the app, it opens QGIS, opens a given project and pans to the correct location (using coordinates in the application). If I were to call QGIS from commandline, what would I need to a) open a given project and b) how would I pass coordinates from one system to automate the map window pan to location. I'm guessing that it is going to be some python code? That being the case, are there any examples of similar ways of starting QGIS? Thanks


Converting from ArcGIS TIN to LandXML?

I am working on a project that involves some back and forth between ArcGIS 10.1 and Bentley Microstation and Geopak. Microstation/Geopak has the ability to import/export LandXML files for 3D surfaces. I have been bringing LandXML files into ArcGIS through the "LandXML to TIN" (3D Analyst) tool. Now that I've done some processing in ArcGIS I'd like to send the surfaces back to MicroStation. Using LandXML as an intermediate again would make sense except that ArcGIS seems to lack a "TIN to LandXML" tool. Does anyone know how to convert a TIN to Land XML? Or, alternatively, any format that can be read by MicroStation/Geopak?


How get Map extents for draw polygon when final click is done?

I am using openLayer js to draw a polygon and now i want map extents of the drawn polygon one final click (i.e. user finish drawing polygon or double click is done). For drawing the polygon i am using following code:               action = Ext.create('GeoExt.Action', {             text: "Draw Poly",             icon: 'Images/poly.png',             control: new OpenLayers.Control.DrawFeature(                 vector, OpenLayers.Handler.Polygon             ),             map: map,             // button options             toggleGroup: "draw",             allowDepress: false,             tooltip: "Draw Polygon",             // check item options             group: "draw"         });      this code is written using openlayer js, extjs and geoextjs. Any please tell how to get map extends for the drawn polygon.


QGIS plugins - how to use a custom behaviour on right click and actionPan on left click?

I have been struggling with this problem for quite a long time. Here we go : I implemented a maptool that draws a red cross on the point that is being right-clicked on. There are good how-tos, such as this topic : How to programatically check for a mouse click in QGIS But I happen to completely lose the left-click behaviour. This means that I cannot navigate through my layer with left click. I used a little trick ("iface.actionPan().trigger()") to get the pan action anyway, but it sets the mapTool to another object, an instance of mapToolPan. This later on causes various bugs, because I am calling methods that belong to my homemade PointTool class on a mapToolPan instance. If the user right clicks on QGIS, the toolPan is set to my homemade PointTool class again. Here is the code where I need help (I only implemented the canvasReleaseEvent because I only use this event) :               class PointTool(QgsMapTool):          def __init__(self, canvas):         QgsMapTool.__init__(self, canvas)         self.canvas = canvas          def canvasReleaseEvent(self, event):         """ Checks which click was used : if left, triggers the pan action, else, adds point. """              if event.button() == QtCore.Qt.RightButton:             # The right click triggers the point selection             x = event.pos().x()             y = event.pos().y()             self.clickedPoint = self.canvas.getCoordinateTransform().toMapCoordinates(x, y)             self.showPoint(self.clickedPoint)             self.canvas.refresh()              elif event.button() == QtCore.Qt.LeftButton:             print("release clicked left")             # The left click triggers the pan action             iface.actionPan().trigger()          def showPoint(self, point):         """ Adds point in qgis main window """              # Using a vertex marker to show the selected point as a red cross         m = QgsVertexMarker(self.canvas)         m.setCenter(point)         m.setColor(QColor(255,0,0))         m.setIconSize(10)         m.setIconType(QgsVertexMarker.ICON_X)              self.canvas.refresh()      So the question is : how can I keep my right-click behaviour AND still have normal left-click behaviour? Questions and remarks are welcome, and of course any help would be greatly appreciated! :)


How to install the QGIS Raster Lang plugin?

I´m trying to install the **rasterlang** plugin for QGIS. I tried several things on this website but it seems to be quite complicated because I need to install binary files and other files. The idea would be to make it the way that when I type in my console: **from rasterlang.layers import layerAsArray** then It would work nicely. Thanks


How to center the map based on some (x,y) in QGIS?

I want to center and recenter the map loaded on QGIS using constantly changing (x,Y) Point. How should I do this using some code? My GPS device will write (x,y) to a text file. From there I get (x,y) and plot two points at a time over a map on QGIS using below Code.                   QUrl url =  QUrl::fromLocalFile("/home/mit/Desktop/Bagfiles/file1.txt");         url.addQueryItem("delimiter",",");         url.addQueryItem("xField", "x");         url.addQueryItem("yField", "y");         url.addQueryItem("crs","epsg:4723");         while (true)         {         emit drawVectorLayer(QString::fromAscii( url.toEncoded() ), "layer", "delimitedtext");        usleep(10000000);         }      Now I want to center the map based on these values.


How CONVERT a lat/lon point TO a x/y point ON a world mercator proj image in JS?

Thx to read and think about my problem! :) I have a map world image like that (mercator proj): http://www.learnnc.org/lp/multimedia/14596 And I have a lat/lon point to place on it in javascript, I am looking for a solution since long time but impossible to solve it.. I have found this example (first of the page) which do this: http://www.netzgesta.de/mapper/ (mouse lat/lon position is write at bottom/left map) but impossible to retrieve the function from the source.. So is there someone who got an idea? And I will need to do the opposite, pointer will be draggable and I will need to convet to lat/lon. thanks Bertrand from france


Getting exceptions with simple hello worldwind java app

Let me preface this by saying I'm very new to Java, eclipse and worldwind. So, I've got a bare bones java app. I imported worldwind.jar into the build path. Created a simple JFrame, code below, and added the following to the Project > properties > Run/Debug settings > Arguments > VM arguments : > -Djava.library.path=Path\To\My\Worldwind\Directory Here is the JFrame source               public class TestMap extends JFrame {          private JPanel contentPane;          public static void main(String[] args) {         EventQueue.invokeLater(new Runnable() {             public void run() {                 try {                     TestMap frame = new TestMap();                     frame.setVisible(true);                 } catch (Exception e) {                     e.printStackTrace();                 }             }         });     }          public TestMap() {         setDefaultCloseOperation(JFrame.EXIT_ON_CLOSE);         setBounds(100, 100, 450, 300);         contentPane = new JPanel();         contentPane.setBorder(new EmptyBorder(5, 5, 5, 5));         contentPane.setLayout(new BorderLayout(0, 0));         setContentPane(contentPane);              JSplitPane splitPane = new JSplitPane();         contentPane.add(splitPane, BorderLayout.CENTER);              WorldWindowGLCanvas worldWindowGLCanvas = new WorldWindowGLCanvas();         splitPane.setRightComponent(worldWindowGLCanvas);              Model m = (Model)WorldWind.createConfigurationComponent(AVKey.MODEL_CLASS_NAME);     }          }      When I run the program, the worldwind panel is black, and i get exceptions streaming through the console output. Apr 18, 2012 2:44:18 PM gov.nasa.worldwind.WorldWindowGLAutoDrawable display SEVERE: Exception while attempting to repaint WorldWindow java.lang.NullPointerException at gov.nasa.worldwind.render.DrawContextImpl.getGlobe(Unknown Source) at gov.nasa.worldwind.view.BasicView.apply(Unknown Source) at gov.nasa.worldwind.AbstractSceneController.applyView(Unknown Source) at gov.nasa.worldwind.BasicSceneController.doRepaint(Unknown Source) at gov.nasa.worldwind.AbstractSceneController.repaint(Unknown Source) at gov.nasa.worldwind.WorldWindowGLAutoDrawable.doDisplay(Unknown Source) at gov.nasa.worldwind.WorldWindowGLAutoDrawable.display(Unknown Source) at com.sun.opengl.impl.GLDrawableHelper.display(GLDrawableHelper.java:78) at javax.media.opengl.GLCanvas$DisplayAction.run(GLCanvas.java:288) at com.sun.opengl.impl.GLDrawableHelper.invokeGL(GLDrawableHelper.java:194) at javax.media.opengl.GLCanvas.maybeDoSingleThreadedWorkaround(GLCanvas.java:265) at javax.media.opengl.GLCanvas.display(GLCanvas.java:130) at javax.media.opengl.GLCanvas.paint(GLCanvas.java:142) at sun.awt.RepaintArea.paintComponent(Unknown Source) at sun.awt.RepaintArea.paint(Unknown Source) at sun.awt.windows.WComponentPeer.handleEvent(Unknown Source) at java.awt.Component.dispatchEventImpl(Unknown Source) at java.awt.Component.dispatchEvent(Unknown Source) at java.awt.EventQueue.dispatchEventImpl(Unknown Source) at java.awt.EventQueue.access$000(Unknown Source) at java.awt.EventQueue$3.run(Unknown Source) at java.awt.EventQueue$3.run(Unknown Source) at java.security.AccessController.doPrivileged(Native Method) at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source) at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source) at java.awt.EventQueue$4.run(Unknown Source) at java.awt.EventQueue$4.run(Unknown Source) at java.security.AccessController.doPrivileged(Native Method) at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source) at java.awt.EventQueue.dispatchEvent(Unknown Source) at java.awt.EventDispatchThread.pumpOneEventForFilters(Unknown Source) at java.awt.EventDispatchThread.pumpEventsForFilter(Unknown Source) at java.awt.EventDispatchThread.pumpEventsForHierarchy(Unknown Source) at java.awt.EventDispatchThread.pumpEvents(Unknown Source) at java.awt.EventDispatchThread.pumpEvents(Unknown Source) at java.awt.EventDispatchThread.run(Unknown Source) Any idea what I'm doing wrong?? EDIT: Forgot to mention, I'm on a 64bit machine, I downloaded the AMD64 version of jogl and put the .dll's in the worldwind folder. That worked to get the demo app running.


Definitions of arcgis script parameter Data Type options

When adding a python script to a toolbox in Arcgis one of the steps is to set parameters and specify what Data Type each requires. The list is long, perhaps a hundred, and many are not immediately understandable. For example I just tried to require data type "Database Connections", but when I run the tool and select "Database Connections\Connection to foobar.sde" I'm told the item is invalid. Aww nuts. Now I could just ask _"what is the "Database Connection" data type?"_ here on GIS.se and get on with my work, but then tomorrow I'll have to do it again for _Disk Connection_ and _Decimate_. (might do wonders for my reputation level though...) So lets just get 'em all in one go, I hope: **Where do I find the definitions of all the Data Types available for tool parameters?** ![enter image description here](http://i.stack.imgur.com/fcLMj.png)


GIS learning resources for a .NET developer

There are couple of questions which are related to learning GIS or GIS skill sets here but I would like to ask very specific question. My current skills: _I don't claim that I am a good programmer_ :-))   * .Net, programming Language C# (4+ years exp) and sybase PowerBuilder (3 years exp)   * Experience working in web and desktop applications.   * SQL   * DBMS (SQL-Server, ASA)   * Javascript, Ajax, HTML, Webservices Can you please help me to find resources which I can use to start GIS development using .net technologies! Regards,


Delete SDE feature class features with ArcPy

I'm working on a script that will empty my SDE database of all data. How would I go about deleting features from feature classes using ArcPy? I have an SDE connection created with a user that should have delete privileges. So far I've tried using the arcpy.DeleteFeatures_management from code and from the Delete Features toolbox in data management. My DB has SDE feature datasets with SDE feature classes inside. The Delete Features tool box errors out with: > ERROR 000732: Input Features: "Dataset Database Connections\Connection to > sde.sde\User.MyFeatureDataset\User.MyFeatureClass" does not exist or is not > supported. The documentation is not really helping me. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//001700000036000000


How to extract raster DEM elevations to vector attribute tables in ArcGIS 10?

I have an ArcMap 10 project with multiple layer groups. One layer is a raster DEM, another layer group is a bunch of vector shapefiles (including points, lines and polygons). The client wants all vector data to contain Z values. Populating the attribute tables for the vector data with X and Y values is easy enough. But how do I populate the Z values? Is there a way that I can use Field Calculator to "grab" the height values from the DEM layer for each vector feature? I should add that we don't have a licence for the 3D Analyst or Spatial Analyst extensions. Many thanks for helping a newbie out :)


how to build a routeable network for pgRouting with srid 3785

I'm following the pgrouting-workshop http://workshop.pgrouting.org/chapters/topology.html and It works, but in my project I'm working with the mercator projection(srid 3785) and there is not a routed network like the Open Street Map to work in, so I need to build my own.. so for testing purposes with the help of google maps and Proj4js I build a simple graph like this ![enter image description here](http://i.stack.imgur.com/92Hbq.png) point1 = "POINT (-7363703.972517611 -1966750.4952863816)"   point2 = "POINT (-7363530.794875402 -1966719.4427436402)"   point3 = "POINT (-7363414.945004408 -1966694.3618437336)"   point4 = "POINT (-7363365.977533163 -1966681.2242294976)"   point5 = "POINT (-7363419.722318676 -1966658.5319867264)"   point6 = "POINT (-7363344.479618957 -1966793.4911147922)"   // it's in lon / lat format Then the geom column in ways table   way_red = "MULTILINESTRING((point1, point2, point3, point4))"   way_green = "MULTILINESTRING((point3, point5))"   way_blue = "MULTILINESTRING((point4, point6))" **Questions:**   1. Is it correct the format I builded the ways ??   2. Exist a correct or pre-established format to construct a routeable network?   3. the workshop uses "MULTILINESTRING", is it ok if I use just "LINESTRING"?? Thanks in advance, and sorry for my english.


Multi Region/Location Use Case on EC2

Context: ArcGIS applications for one of our customers are currently running on multiple sites in Americas, Europe, Australia and Africa geos (15 different countries). They are considering a cloud model. As you all are aware AWS provides 5 regions and support few edge locations for each of the regions. Challenge: The coverage is not adequate for the current scenario. Considering that the applications are currently hosted at 15 different locations. Running intensive analysis on a ArcGIS cloud based solution results into a substandard user experience is my understanding. Question: My question to the experts on this forum is cloud a suitable model to be recommended for my client. I have read about caching and elastic caching but are they really going to be effective considering the geo spread. My major concern is around latency with EC2s limited coverage and hence user experience will be taking a big hit. Case Studies: I am keen to hear about case studies of similar nature use cases if any. Thanks in advance. Regards, Sameer


How to find most frequent values?

**Entry layers:** _polygons_ \- no attributes (columns/values), _points_ \- needed attributes (columns/values) **How to get this result?: polygon layer with attribute (column/value) of most often (most frequent - MODUS) of values from points values within polygons** _Probably best for me with SQL select in Mapinfo but I can convert layers and do it in Qgis or ArcGIS if you know the simple way in these SW._ Additionally: How to get count of points within polygons and add it as column (attribute) to polygon layer in Mapinfo? - **solved - this not help to main question** Example of points within polygons - polygons represented by outlines (a lot of overlapping - it is import from "nonGIS" software through dxf - nobody cared about topology etc.): ![Example of points within polygons](http://i.stack.imgur.com/BO6ZS.png) Here is example how to get **maximum/minimum/mean**... there is a simple solution in ArcMap - right click on layer - _Joins and Relates_ / _Join..._ / _Join data from another layer based on spatial location_ (as described here: In ArcGIS how can I append the number of points that are located inside each individual block groups of the entire US?) - but I need **MODUS**.


How to get WKT from SHAPE column in the ArcSDE geodatabase?

I have an ArcSDE geodatabase running on top of SQL Server 2008. A table for a feature layer was added from ArcCatalog. The SHAPE column is a type of int which was automatically converted from ArcCatalog. **My question is:** In SQL Server Management Studio, what is the best SQL syntax to output WKT for each of the records in the table? p.s. `[SHAPE].STAsText()` doesn't work in               SELECT [OBJECTID] ,[SHAPE], [SHAPE].STAsText() FROM [sde].[sdl_admin].[POINT_LOCATION]          Cannot find either column "SHAPE" or the user-defined function or aggregate      SHAPE.STAsText", or the name is ambiguous.


issues importing arcpy in Python IDLE

Hi I would want to use arcpy in my python code. I am able to import it through the ARCGIS desktop python console. But I am not able to import it into Python IDLE. I get the following error               Traceback (most recent call last):          File "C:\Python26\lib\site-packages\win32com\client\dynamic.py", line 246, in      __getitem__raise TypeError("This object does not support enumeration")       I mentioned the following folders as a part of my python path.               C:\Program Files\ArcGIS\Desktop10.0\arcpy;     C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy;     C:\Python26\ArcGIS10.0\Tools\Scripts      I found the same issue being mentioned. import arcpy yields "TypeError: This object does not support enumeration" The user says it is due to the following: > I wrote for Arcgis 9 and 10 in the same directory. Something about the > initialization logic in import arcpy finds the custom arcgisscripting.py > used for running python 2.6 with arcgis 9.3 when they are located together. I have only ARCGIS 10 and i am not aware of how arcpy looks up for its modules. Please do help me resolve this issue.


Is there a zonal statistics tool to calculate median and mode raster values for each polygon?

I'd like to calculate zonal statistics for a raster where the zones are defined by polygons in a vector layer within QGIS. Specifically I'm interested in calculating the median and mode raster values for each feature. QGIS includes a 'zonal statistics' plugin, but it only seems to calculate mean and sum statistics. Is there another plugin that could do this task for me? I've looked at the source code for the existing statistics plugin and it looks like it would be fairly easy to modify, although my skills in C++ are fairly lacking. QgsZonalStatistics source: http://www.qgis.org/api/qgszonalstatistics_8cpp_source.html How to use QGIS Zonal Stats Plugin from Python Console? Raster statistics available in ArcGIS: http://resources.arcgis.com/en/help/main/10.1/index.html#//009z000000w7000000


arcpy: can I use SelectLayerByAttribute to compare fields from joined tables?

I can't figure out the syntax to select table records where one field's value exceeds that of another. The fields are from a pair of table views connected w/AddJoin(). One view is of a fgdb table and the other is of a point feature class in a separate fgdb. My ArcGIS10.0 syntax looks like:               arcpy.MakeTableView_management(r'C:\Workspace\wk_NHDPlus\NHDPlus\NHDPlus_national.gdb\POIs\GagesII', 'gagesTV', '"REGION" = \'02\' AND "COMID" = 0')     arcpy.AddJoin_management('msrTV', 'ComID_Measures_REACHCODE', 'gagesTV', 'REACHCODE', 'KEEP_ALL')        arcpy.SelectLayerByAttribute_management('msrTV', 'NEW_SELECTION', '"GagesII.MEASURE" >= "ComID_Measures02.ComID_Measures_BOTMEAS" ')      I've also tried:               newName = arcpy.AddFieldDelimiters(r'C:\Workspace\wk_NHDPlus\NHDPlus\NHDPlus_national.gdb\ComID_Measures02', "ComID_Measures02.ComID_Measures_BOTMEAS")     arcpy.SelectLayerByAttribute_management('msrTV', 'NEW_SELECTION', '"GagesII.MEASURE" > '+newName)      The result in both cases is the entirety of the original 'msrTV', which I know is incorrect. Is what I'm trying to do possible? Hopefully I'm just fouling up the syntax. _*_ update: also tried:               newName = arcpy.AddFieldDelimiters(r'C:\Workspace\wk_NHDPlus\NHDPlus\NHDPlus_national.gdb\ComID_Measures02', "ComID_Measures02.ComID_Measures_BOTMEAS")     oldName = arcpy.AddFieldDelimiters(r'C:\Workspace\wk_NHDPlus\NHDPlus\NHDPlus_national.gdb\ComID_Measures02', "GagesII.MEASURE")     sqlquery=oldName+'>='+newName     arcpy.SelectLayerByAttribute_management('msrTV', 'NEW_SELECTION', sqlquery)


Where are the Block and Lot fields in the ESRI Local Government Information Model Geodatabase Schema?

I am looking at the ESRI Local Government Information Model (LGIM) geodatabase schema and I cannot find the Feature Class that holds the parcels' Block and Lot information. Do you know which Feature Class in LGIM schema holds the Block&Lot information? Thanks


how to replace GTileLayer in google maps API v3

I need to port code of one application written in Google Mapa API v.2 to API v.3. There is several layers taken via TMS, that are displayed using GTileLayer class. It is as simple as:               var tilelayer_tematyczne = new GTileLayer(null, null, null, {         tileUrlTemplate: 'http://domain/geoserver/gwc/service/gmaps?layers=layer&zoom={Z}&x={X}&y={Y}&format=image/png',         isPng:true,         opacity:1      });          tilesTematyczne = new GTileLayerOverlay(tilelayer_tematyczne);     map.addOverlay(tilesTematyczne);      I do not see any similar class in API v3. How do you cope with that problem? How to display TMS on top of Google Maps?


Creating non-spatial Many to One Joins/Relationships between tables

I may have missed a trick here but I can't seem to find a solution. I have a table of ~1.6 million postcodes, also containing co-ordinates, and a separate file detailing ~6000 a population with postcodes attached. I am trying to attach the co-ordinates from the postcode file to the population file. There is, but not always, two or more individuals associated to one postcode. The theory, in my head, is simple. Join the postcode table to the population table on the postcode attribute. To start with this idea in ArcGIS 10.1 seemed simple, using a basic join, but what seems to be happening is that a 1:1 relationship is being enforced. After researching further I found the Make Query Table tool which claimed to solve this M:1 problem. Though that doesn't seem to be working. I am currently in a two-step process to see if I can solve this. First is looking through some, don't particularly want to look through all files in either table, to see if there is a formatting error that would explain this error. The second is looking into using MySQL to see if that can be used to join the tables. Or at least create a query of the two tables that will produce the desired output table. First off if you have any advice on making this work in ArcGIS that'd be brilliant, I'm sure it can do it, just not sure how to make it work. Alternatively if you know of another method that would easily join these tables, be that QGIS, MySQL or something else, then please do say. Hope you can help, and if there's anything I've missed out that would help you understand my problem, let me know.


Difference Between GLC2000 Versions

What are the differences between the various versions (1 through 5) of the Global Land Cover (GLC) 2000 product? I am going to use northern North America, Europe and Asia (north of 50°N). Is it better to mosaic the continental products (even though they are different versions) or use the global product (v 1.0). I want the best data source that I can get for land classification, so if anyone has other recommendations, I would appreciate it.


QGIS and Ubuntu 'browse' dialogue causing crashing

Below is an image of 'the suspect dialogue' in Ubuntu that when you encounter it in QGIS you are most definitely going to crash the application. http://ow.ly/i/ZGNc/original You find this particular dialogue - which doesn't appear with every menu - when you use: \- polygon centroids \- raster extraction / projection \- create new shapefile (from desktop when you browse to location of new file) ubuntu 12.04 qgis 1.8


How can I split a segment into two segments?

The problem is how to split segments or in other words how to insert new points into a polygon so that the segment where I created the new point is now two segments connected by said point. Every time I create a new point within a segment, the segment stays the same length as before and is not being split by the new point. (I don't know why the following snapshots have double columns on the left side but just ignore that. was not like this in the program). First I have the segment I want to split: ![segment is dark red between two points](http://i.stack.imgur.com/sDfi6.jpg). Then I go to edit mode for the segment layer, click on the knot tool and select the segment: ![location where the segment should be split](http://i.stack.imgur.com/7ictB.jpg). I double-click on the segment where i want a new point. The "new place for the point" appears: ![enter image description here](http://i.stack.imgur.com/u1idA.jpg). I then go to pointlayer, activate the edit mode, click on "add new element" and create a new point near the segment: ![new point is shown near the "x"](http://i.stack.imgur.com/kJYZl.jpg). After I activate the snap-function and define the radius, I place the point on the "x" with the knot-tool: ![point is now in the right place](http://i.stack.imgur.com/N74vp.jpg). Now if I go on the info-pointer and click on the segment, the segment is still as one piece and is not split by the new point: ![still one segment](http://i.stack.imgur.com/5qQbW.jpg). Why?


Technique using GDAL to overlay / fit one PNG map to directly overlay on another?

I have one .PNG basemap with known North edge of Latitude = 52.48, a South edge of Latitude = 21.94, an east edge of Longitude = -61.87 and a west edge of Longitude = -129.37. And I have another PNG image overlay with a known North edge of Latitude = 50.03, a South edge of Latitude = 22.02, an east edge of Longitude = -67.24 and a west edge of Longitude = -124.99. Is there a way using GDAL or other simple scripting tools to process the second PNG image to be in an output that directly overlays (same x/y dimensions) the basemap? I have no trouble generating map tiles that overlay each other beautifully - but I seem to be having no luck figuring out the syntax of the seemingly simpler case. Help in the form of a simple example is much appreciated! To be clear: the goal is an output .png file that has the exact same x/y dimensions as the basemap, and which can directly overlay on top of it.


Selection and Popups on Overlapping/Coincident Vector Layer Polygons in OpenLayers

I have a GeoJSON file of polygons that I am bringing into my OpenLayers session. Some of the polygons occupy the same space on the map. I am using a select control and popup control like:-               selectControl = new OpenLayers.Control.SelectFeature(geojson_layer,     {         onSelect: onFeatureSelect, onUnselect: onFeatureUnselect,         clickout: true, toggle: true,         multiple: false, hover: false,         toggleKey: "ctrlKey", // ctrl key removes from selection         multipleKey: "shiftKey" // shift key adds to selection     });     map.addControl(selectControl);     selectControl.activate();     geojson_layer.events.on(     {         "featureselected": function(e)          {             showStatus("selected feature "+e.feature.attributes.APPREF+" on  Applications");         },         "featureunselected": function(e)              {             showStatus("");             }         });                function onFeatureSelect(feature)              {                 selectedFeature = feature;                 popup = new OpenLayers.Popup.AnchoredBubble                 (                        "chicken",                      feature.geometry.getBounds().getCenterLonLat(),                     null,                     "
App Reference: " + feature.attributes.APPREF +"
Proposal:
" + feature.attributes.PROPOSAL+ '
Further Details'+"
",                     null,                      true,                      onPopupClose                 );                 feature.popup = popup;                 popup.autoSize = true;                 popup.minSize = new OpenLayers.Size(150,50);                 popup.maxSize = new OpenLayers.Size(250,150);                 //popup.panMapIfOutOfView = true;                 popup.closeOnMove = true;                 map.addPopup(popup);             }             function onFeatureUnselect(feature)              {                 map.removePopup(feature.popup);                 feature.popup.destroy();                 feature.popup = null;             }      Where more than 1 polygon occupies the same space - can anyone advise how I allow for drilling down and/or making the other features at the location available for a popup (and therefore the link out for further information). I suspect the answer lies somewhere in Z-indexing, but can't wrap my head around how I would integrate that from the example thats provided. Presumably there is a way of handling this type of scenario - if someone could enlighten me that would be very much appreciated Thanks


Extent Intersection when WrapAround180 enabled

Allowing WrapAround180 allows map extents to go beyond the true Web Mercator extent. Once this happens I am not able to use map.extent.intersects(point). My work around was to   1. Normalize extents in lat/long   2. Convert point to lat/long   3. Split extents to two parts if crossing date line.   4. Test intersection with extent(s) Is there an easier way? Here's my code: var extent = map.extent; var pointToTest = someGraphic;                    var normMin = wmUtils.xyToLngLat(extent.xmin, extent.ymin);                var normMax = wmUtils.xyToLngLat(extent.xmax, extent.ymax);                    //test if minx is greated than maxx          if (normMin[0]> normMax[0]) {                var pointMin1 = new Point([normMin[0],normMin[1]],new SR({ wkid:4326 }));           var pointMin2 = new Point([180,normMax[1]],new SR({ wkid:4326 }));           var wmMin1 = wmUtils.geographicToWebMercator(pointMin1);           var wmMin2 = wmUtils.geographicToWebMercator(pointMin2);                var pointMax1 = new Point([-180,normMin[1]],new SR({ wkid:4326 }));           var pointMax2 = new Point([normMax[0],normMax[1]],new SR({ wkid:4326 }));           var wmMax1= wmUtils.geographicToWebMercator(pointMax1);           var wmMax2= wmUtils.geographicToWebMercator(pointMax2);                     var normExtent1 = new Extent(wmMin1.x,wmMin1.y,wmMin2.x,wmMin2.y, new SR({ wkid:102100 }));           var normExtent2 = new Extent(wmMax1.x,wmMax1.y,wmMax2.x,wmMax2.y, new SR({ wkid:102100 }));              } else {               var pointMin = new Point(normMin,new SR({ wkid:4326 }));          var wmMin = wmUtils.geographicToWebMercator(pointMin);          var pointMax = new Point(normMax,new SR({ wkid:4326 }));          var wmMax = wmUtils.geographicToWebMercator(pointMax);          var normalizedExtent = new Extent(wmMin.x,wmMin.y,wmMax.x,wmMax.y, new SR({ wkid:102100 }));             }             if (normMin[0]> normMax[0]) {             if ((!normExtent1.intersects(pointToTest.geometry) && !normExtent2.intersects(pointToTest.geometry))) return;           } else {                         if (!normalizedExtent.intersects(pointToTest.geometry)) return;        }


QGis: Difference between Mac OSx and Windows

i try to import a KML-Placemark file with QGIS (--> Add Vectorlayer, UTF-8). When i import the File with my Mac (OS 10.6) i don´t have any problems viewing all the attributes stored in the KML-File (they are stored in the Tag "ExtendedData") ![enter image description here](http://i.stack.imgur.com/CIDao.png) but when i import the same KML-File in the same QGis Version on a Windows 8 OS i only see the attributes "Name" and "Description" (The points are shown correctly) ![enter image description here](http://i.stack.imgur.com/btjyo.png) this is referring to my question here: KML in QGis with additional Data i thought the problem is solved but it is only to 50 % :) here is my updated kml-file:               <?xml version="1.0" encoding="UTF-8"?><kml xmlns="http://www.opengis.net/kml/2.2">     <Document>     <Placemark id="PM1">     <TimeStamp>             <when>2007-01-14T21:05:02Z</when>           </TimeStamp>     <name>Data1</name>     <description>AA</description>     <ExtendedData>     <Data name="Orientierung">     <displayName>Orientierung</displayName>     <value>1</value></Data>     <Data name="Brennweite">     <displayName>Brennweite</displayName>     <value>279/100</value>     </Data>     <Data name="Datum">     <displayName>Datum</displayName>     <value>2013:04:04 13:20:49</value>     </Data>     <Data name="Pfad">     <displayName>Pfad</displayName>     <value>file:///C:\Temp\20130404_132050.JPG</value>     </Data>     <Data name="Kamera">     <displayName>Kamera</displayName>     <value>GT-P3110</value></Data>     </ExtendedData>     <Point>     <coordinates>15.448561,47.11721,0</coordinates>     </Point>     </Placemark>          <Placemark id="PM2">     <TimeStamp>             <when>2007-01-14T21:05:02Z</when>           </TimeStamp>     <name>Data2</name>     <description>AA</description>     <ExtendedData>     <Data name="Orientierung">     <value>5</value></Data>     <Data name="Brennweite">     <value>400/400</value>     </Data>     <Data name="Datum">     <value>2013:04:04 13:20:50</value>     </Data>     <Data name="Pfad">     <value>file:///C:\Temp\HUGO.JPG</value>     </Data>     <Data name="Kamera">     <value>GT-P3110</value></Data>     </ExtendedData>     <Point>     <coordinates>16.458561,48.15721,0</coordinates>     </Point>     </Placemark>     </Document>          </kml>      does anyone know why this problem is appearing?? hope you can help me!   thanks michik


Geoserver won't cache wms server with style in workspace

I am using Geoserver 2.5 on ubuntu and am trying to cache a layer that has it styles in workspaces. When I try to start the cache with seed/truncate nothing happens and I get things like this in the logs:               2014-04-07 10:56:03,939 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - Converted URL to lowercase, from: '/gwc/rest/seed/wis_gent:maaiveld'; to: '/gwc/rest/seed/wis_gent:maaiveld'  and httpMethod= GET     2014-04-07 10:56:03,939 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - ~~~~~~~~~~ antPath= /** methodList= [GET]     2014-04-07 10:56:03,939 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - Candidate is: '/gwc/rest/seed/wis_gent:maaiveld'; antPath is /**; matchedPath=true; matchedMethods=true     2014-04-07 10:56:03,939 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - returning ADMIN     2014-04-07 10:56:03,940 TRACE [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/gwc/rest/seed/wis_gent:maaiveld]     2014-04-07 10:56:03,941 DEBUG [org.geoserver.filters] - Compressing output for mimetype: text/html;charset=ISO-8859-1     2014-04-07 10:56:03,942 DEBUG [org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1] - SecurityContextHolder now cleared, as request processing completed     2014-04-07 10:56:04,028 TRACE [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/gwc/rest/web/geowebcache_logo.png]     2014-04-07 10:56:04,029 DEBUG [org.geoserver.filters] - Not compressing output for mimetype: image/png     2014-04-07 10:56:04,034 DEBUG [org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1] - SecurityContextHolder now cleared, as request processing completed     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - Converted URL to lowercase, from: '/gwc/rest/seed/wis_gent:maaiveld'; to: '/gwc/rest/seed/wis_gent:maaiveld'  and httpMethod= POST     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - ~~~~~~~~~~ antPath= /** methodList= [GET]     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - Candidate is: '/gwc/rest/seed/wis_gent:maaiveld'; antPath is /**; matchedPath=true; matchedMethods=false     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - ~~~~~~~~~~ antPath= /** methodList= [POST, DELETE, PUT]     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - Candidate is: '/gwc/rest/seed/wis_gent:maaiveld'; antPath is /**; matchedPath=true; matchedMethods=true     2014-04-07 10:56:05,610 DEBUG [org.geoserver.security.RESTfulPathBasedFilterInvocationDefinitionMap] - returning ADMIN     2014-04-07 10:56:05,611 TRACE [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/gwc/rest/seed/wis_gent:maaiveld]     2014-04-07 10:56:05,613 DEBUG [org.geoserver.gwc.layer] - --> GWC Seeder Thread-6 submitting getMap request for meta grid location [0, 0, 0] on org.geoserver.gwc.layer.GeoServerMetaTile@27897fc     2014-04-07 10:56:05,614 INFO [org.geoserver.wms] -      Request: getServiceInfo     2014-04-07 10:56:05,616 DEBUG [org.geoserver.ows] - Getting layers and styles from LAYERS and STYLES     2014-04-07 10:56:05,616 ERROR [org.geoserver.ows] -      org.geoserver.platform.ServiceException: No such style: wis_bestemming         at org.geoserver.wms.map.GetMapKvpRequestReader.parseStyles(GetMapKvpRequestReader.java:1386)         at org.geoserver.wms.map.GetMapKvpRequestReader.read(GetMapKvpRequestReader.java:318)         at org.geoserver.wms.map.GetMapKvpRequestReader.read(GetMapKvpRequestReader.java:80)         at org.geoserver.ows.Dispatcher.parseRequestKVP(Dispatcher.java:1412)         at org.geoserver.ows.Dispatcher.dispatch(Dispatcher.java:622)         at org.geoserver.ows.Dispatcher.handleRequestInternal(Dispatcher.java:263)         at org.springframework.web.servlet.mvc.AbstractController.handleRequest(AbstractController.java:153)         at org.geoserver.gwc.GWC.dispatchOwsRequest(GWC.java:1162)         at org.geoserver.gwc.layer.GeoServerTileLayer.dispatchGetMap(GeoServerTileLayer.java:572)         at org.geoserver.gwc.layer.GeoServerTileLayer.getMetatilingReponse(GeoServerTileLayer.java:517)         at org.geoserver.gwc.layer.GeoServerTileLayer.seedTile(GeoServerTileLayer.java:711)         at org.geowebcache.seed.SeedTask.doActionInternal(SeedTask.java:136)         at org.geowebcache.seed.GWCTask.doAction(GWCTask.java:76)         at org.geowebcache.seed.MTSeeder.call(MTSeeder.java:36)         at org.geowebcache.seed.MTSeeder.call(MTSeeder.java:25)         at java.util.concurrent.FutureTask.run(FutureTask.java:262)         at java.util.concurrent.ThreadPoolExecutor.runWorker(ThreadPoolExecutor.java:1145)         at java.util.concurrent.ThreadPoolExecutor$Worker.run(ThreadPoolExecutor.java:615)         at java.lang.Thread.run(Thread.java:744)     2014-04-07 10:56:05,617 ERROR [org.geowebcache.seed.MTSeeder] - Problem communicating with GeoServer     org.geowebcache.GeoWebCacheException: Problem communicating with GeoServer         at org.geoserver.gwc.layer.GeoServerTileLayer.getMetatilingReponse(GeoServerTileLayer.java:522)         at org.geoserver.gwc.layer.GeoServerTileLayer.seedTile(GeoServerTileLayer.java:711)         at org.geowebcache.seed.SeedTask.doActionInternal(SeedTask.java:136)         at org.geowebcache.seed.GWCTask.doAction(GWCTask.java:76)         at org.geowebcache.seed.MTSeeder.call(MTSeeder.java:36)         at org.geowebcache.seed.MTSeeder.call(MTSeeder.java:25)         at java.util.concurrent.FutureTask.run(FutureTask.java:262)         at java.util.concurrent.ThreadPoolExecutor.runWorker(ThreadPoolExecutor.java:1145)         at java.util.concurrent.ThreadPoolExecutor$Worker.run(ThreadPoolExecutor.java:615)         at java.lang.Thread.run(Thread.java:744)     Caused by: java.lang.IllegalStateException: Expected: RenderedImageMap, got null         at org.geoserver.gwc.layer.GeoServerTileLayer.dispatchGetMap(GeoServerTileLayer.java:575)         at org.geoserver.gwc.layer.GeoServerTileLayer.getMetatilingReponse(GeoServerTileLayer.java:517)         ... 9 more      Using a style without a workspace works great. Does anyone have any idea what could be wrong and if I can use a style in the workspace?


python.multiprocessing and “FATAL ERROR (INFADI) MISSING DIRECTORY”

I have a question which is closely related to this old question here: python.multiprocessing and "FATAL ERROR (INFADI) MISSING DIRECTORY" I have also looked at this, which seems to relate to something similar: Multiprocessing Errors - ArcGIS implementation just like in the linked question, I run into the error: "FATAL ERROR (INFADI) MISSING DIRECTORY" at random times, but never the BND-error I have multiple folders with thousands of shapefiles in each, where I want to convert two of the shapefile fields to two different rasterlayers. the shapefiles were converted from csvs via multiprocessing, with no issues whatsoever. I have tried most of the workarounds from the previous question, such as creating temporary and unique scratchworkspaces and workspaces, check that the infile exists and create a string in memory and flush it, but to no avail. I have also tried splitting the workload, so that it will first write all rasters for one field, then do all the rasters for the second field. so my question is, would anybody know of another potential workaround? I was thinking that it might be possible to somehow make sure that when a new raster is written to the folder, that no other processes try to access the space, but I'm unsure how to code it. relevant code excerpts below:               import arcpy     import os     import sys     import arcgisscripting     import multiprocessing     from functools import partial          def multiproc(infile,path,savepath1,savepath2):           try:                 arcpy.Exists(infile)                  # Set local variables             inLayer = path + "sample.xml"             fieldList = arcpy.ListFields(infile[:-4]+".dbf", "field1")             inData = infile + " field1"             outLayer = infile[:-4]                  # Check out the ArcGIS Geostatistical Analyst extension license             arcpy.CheckOutExtension("GeoStats")             # Execute CreateGeostatisticalLayer             arcpy.GACreateGeostatisticalLayer_ga(inLayer, inData, outLayer)                arcpy.GALayerToGrid_ga(outLayer, savepath1 + "/" + infile[:-4])             arcpy.Delete_management(outLayer)                  #starting second raster creation process             arcpy.Exists(infile)             fieldList = arcpy.ListFields(infile[:-4]+".dbf", "field2")             inData = infile + " field2"             outLayer = infile[:-4]                  # Check out the ArcGIS Geostatistical Analyst extension license             arcpy.CheckOutExtension("GeoStats")             # Execute CreateGeostatisticalLayer             arcpy.GACreateGeostatisticalLayer_ga(inLayer, inData, outLayer)             arcpy.GALayerToGrid_ga(outLayer, savepath2 + "/" + infile[:-4])             arcpy.Delete_management(outLayer)                   except:             import logging             LOG_FILENAME = 'c:/error.txt'             logging.basicConfig(filename=LOG_FILENAME,level=logging.DEBUG,)                  logging.debug(infile)                  logging.exception('Got exception on main handler')             pass               if __name__=='__main__':     # Set environment settings     #setting paths for the directories           ##Finding paths and list of files                #multiprocessing part         Pool=multiprocessing.Pool         workers=multiprocessing.cpu_count()         pool=Pool(processes=workers)              partial_multiproc=partial(multiproc,path=path,savepath1=savepath1,savepath2=savepath2)         pool.map(partial_multiproc,filelist)         pool.close()         pool.join()      EDIT1: after having moved the data and tried on an internal harddrive (the previous attempts have been on an external HD), I also got the BND-error on rare occassions, whereas the INFADI is much more prevalent. EDIT2: I have been unable to find a full solution so far, but a potential workaround lies in creating temp-folders that are deleted with arcpy at the end of the multiproc-function. even though the error persists, the faulty folders are - as far as I can tell - not deleted, meaning that at least it is possible to know which rasters were not completed and redo these later.


Coarse heat map placement using gdal

I have generated a coarse heatmap of where data is available by 1 degree x 1 degree cells that I want to overlay on top of a world image. Using a python script I merged two images (heat map and a world map -- generated with gdal_rasterize). Reviewing the picture I have a vertical registration issue (I am guessing a flattening issue). The picture can be sent if needed but I cannot post it at this time (I do not have a reputation >= 10). How can I use gdal to merge these images correctly? Thanks


Off by about 85 meters when converting NAD27 UTM to WGS84 lat/lon

Looking through other answers on converting to or from NAD27 CONUS, I don't see an answer that seems to fit my problem. I am creating paper topographic maps by layering information derived from USGS DEM data (NAD83) for contours and slope shading. I am using OpenStreetMap (OSM) data (WGS84) for trails, roads, water features, etc. Finally I am overlaying a UTM grid. All this is being formatted using Mapnik. Everything is aligning well (trails and road appear correct with respect to the contours. The waterways from OSM also nicely match up with the drainages as shown by contours derived from USGS DEM data, etc.) And if I use WGS84 for my UTM grid it seems to line up well too, at least the positions of map features are the same within my measuring capability as a USGS topo map. But if I overlay a NAD27 UTM grid my easting grid lines appear to be off by about 85 meters. The northing grid lines appear to be correct. I could live with 5 or 10 meters as that is the accuracy of my GPS but not 85 meters. I create the UTM grid by determining the UTM measurements for the edges of my map and then going to points the nearest even 1000 meters off paper and putting ways/lines between the points so located. The position of the points is converted to from UTM to lat,lon and the points and ways are written to a OSM XML file which is used as the source for the UTM grid layer in Mapnik. Conversion of UTM to lat,lon is using pyproj in a Python script. My instantiation of the projection objects is as follows:               wgs84utm = pyproj.Proj(proj='utm', zone=v['zone'], datum='WGS84', ellps='WGS84')     nad27utm = pyproj.Proj(proj='utm', zone=v['zone'], datum='NAD27', ellps='clrk66')      When I convert a point from UTM to lat,lon I use the following:               if (proj == 'wgs84'):         point = wgs84utm(easting, northing, inverse=True)     elif (proj == 'nad27'):         point = nad27utm(easting, northing, inverse=True)      One thing that confuses me is I only see one NAD27 UTM setup in the pyproj espg list. Based on my GPS and other mapping programs I was expecting to see several (CONUS, Alaska, Canada, etc.). There is a NAD27 CONUS but it is an Albers projection, not a UTM projection. Based on values for that I tried               nad27utm = pyproj.Proj(proj='utm', zone=v['zone'], datum='NAD27', ellps='clrk66', lat_1=29.5, lat_2=45.5, lat_0=23, lon_0=-96, x_0=0, y_0=0)      But the results were the identical. I am a newbie with python and map projection conversions, so I'm pretty sure that I've made some mistake in my usage of the projection library but don't know what. I'd appreciate any suggestions on where I have gone wrong. Thanks!


Building real-time point/polyline connection in GIS

In the first instance I need to create the network by constructing the basic elements, the nodes and links. This is relatively easy work as SharpMap has provided geometry tools that allow users to directly create point/polyline via simple mouse clicks e.g. capturing coordinate to generate point and two coordinate to construct line. These functions have been encapsulated and can be directly used. For this task, however, every time I construct a new polyline, I need to keep it connected with two nodes that I manually created, so then the network is connected otherwise the node and link will be separated. I really want to know the programming procedure on how to do this in real- time. That is, every time I create a line, I hope it will search the existing points and display if there are existing points that can be connected. If it can be connected, e.g. by displaying different colour, and then I can snap one of nodes (starting node or ending node) of the line on that point. If I do this way, I think I can construct the connected network as required. I am wondering if anyone has such experience to build such function, or you can provide me some guidance in this aspect? Thanks for the helps.


GeoServer map layers/tiles Loaded in the screen or Not

I'm using OpenLayers - 2.11,GeoServer 2.13, ExtJS 3.4, and Java for developing GIS based web applications. Please guide me to check if all my GeoServer map layers/tiles Loaded in the screen or not. While dragging map is taking time to load next tiles. So that time i want to mask all the map body and unmask after loading all the map tiles. ![enter image description here](http://i.stack.imgur.com/lOVyq.png)


how to remove and change qgsvertexmarker in QGIS using c++ code?

In QGIS following code is used to set marker at point A.                 QgsVertexMarker * marker = new QgsVertexMarker(mPlugin->iface()->mapCanvas());       marker->setColor(Qt::black);       marker->setIconType(2);       marker->setCenter(A);      How to delete this marker using C++ code? In python it is                canvas.scene().removeItem(marker)      I am using this marker in a loop where point A will change continuously so i would like to delete marker at previous position and insert in new. when i use below code marker is deleted permanently.                 mPlugin->iface()->mapCanvas()->scene()->removeItem(marker);      so how to use this marker to plot moving point?


How to work with the STRtree Index when object in the index change it's position?

I'm work with the NTS(NetTopologySuite), I get large amounts of Point geometries, for fast querying the point object, I make the STRtree Index like this, _STRtree tree = new STRtree(); foreach (var p in allPoint) { tree.Insert(p.EnvelopeInternal, p); }_ While some point objects in the "allPoint" change it's position. And I find I can't query the point which change position anymore.. How to work with it? Thanks.


Access Stamen tiles flavors with OpenLayers

i'm using OpenLayers (v2.13.1) and would like to use Stamen Toner tiles as background. The example is simple to understand : http://maps.stamen.com/test/openlayers.html But i'd like to use the "lite" variant of the toner tiles, how is it possible with OpenLayer ? Thanks


CalculateField_management

How to put the field from arcpy.GetParameterAsText() into the expression here: CalculateField_management(in_table, field, **expression** , {expression_type}, {code_block}) Now it is like **"!Name!.replace(\"о\",\"r\")"** but here the field _!Name!_ is already firmly defined and isn't depend on arcpy.GetParameterAsText() thank you


Can we have a local BingMap server or Google Earth server?

Is it possible to have all the map data come from local server instead of Google Earth or Bing Map but still use their tools as front end. For example Bing Map control (from their SDK) would be talking to some local server in the network. I want this kind of setup in a place where there is no internet and the client wants to use Bing Map toolset. Same question applies to Google Earth. OR is it that its only possible with something like ArcGIS like server software but in that case Google Earth or Bing Map control won't be its front end. Please excuse me for my lack of understanding of GIS in general


Field Calculator If/Then in ArcGis 10.1

I've been looking through these threads, i found something similar and thought it could be the answer, but it doen't seem to work for me... I could do this the slow way with Select By Attributes, but It's time consuming and I will have to repeat it several times for other works. I have a field called "Kultur" (String) with different codes (DG, A, FF, ...). I want to change this information into a numeric information, so i added a new field called "Kultur_Cl" (Short Integer). In my case if "Kultur" is DG, "Kultur_Cl" should be 1. If "Kultur" is A, "Kultur_Cl" should be 2, ... Here is what i tried : ![enter image description here](http://i.stack.imgur.com/JFKae.png) And here is what it looks like in ArcMap: ![enter image description here](http://i.stack.imgur.com/d1nDV.png) Once it's done i would like to do : If "Kultur_Cl" = 1 then "Kultur_Pkt" = 15, if "Kultur_Cl" = 2 then "Kultur_Pkt" = 18, ... Would this go the same way, or is it different because there would only be numeric data? Thanks for your help !!


Geometry gets modified when updating only attribute data

Using ArcGIS Runtime for Windows Mobile 10.1. We are trying to update attribute data of features and somehow the geometry gets also updated when we save the changes to the mobile cache and synchronize them back to the database. If I update multiple features, all previously updated features get the geometry of the feature that was last updated. Here is an example code which reproduces the issue. The firstRow gets the geometry of secondRow when the data is synchronized back to the database. In mobile cache everything seems to be correct, atleast I can see the features correctly on the map and select them. And also the attribute data is correctly synchronized.                       int firstId = 2105650708;             int secondId = 2105650709;                  MobileCache mobileCache = new DatabaseService().mobileCache;                  FeatureSource source = mobileCache.FeatureSources["TABLE"];                  FeatureDataTable firstDt = source.GetDataTable(new QueryFilter(                 String.Format("{0} = {1}", Columns.OBJECTID, firstId)));             FeatureDataRow firstRow = firstDt[0];             firstRow[Columns.DESCRIPTION] = "FIRST";             firstDt.SaveInFeatureSource();                  FeatureDataTable secondDt = source.GetDataTable(new QueryFilter(                 String.Format("{0} = {1}", Columns.OBJECTID, secondId)));             FeatureDataRow secondRow = secondDt[0];             secondRow[Columns.DESCRIPTION] = "SECOND";             secondDt.SaveInFeatureSource();      Has anyone experienced anything similar?


How to get turn angle, line orientation and line length in ArcMap 9.3.1?

I am mapping and analysing animal movements using arcgis 9.3.1 along with Hawthorne Beyer's GME, and am trying to utilise the correlated random walk commands, but in order to do this I need to find out the turn angles, step length and line orientation. Is there a way to easily do this in arcmap for an entire vector layer? Such as is done with calculating total or portion line length by adding a new field and calculating geometry? I am not adept at using scripts, so if the solution requires this, can you also include instruction on how to run it? Very many thanks in advance,