Open Source Remote Sensing vs eCognition, ERDAS etc Is it worth it?

What, if any, would be the main benefit of using a commercial tool like eCognition (in particular) or ERDAS as opposed to something like Orfeo Toolbox or OpenCV? (Feel free to suggest other open source software for RS if you think it fills a gap also.) Put aside the GUI as I would likely be scripting it anyway. It seems to me that most of the fundamentals of Remote Sensing are available as easily scriptable OSS elsewhere. Things like mosaicing, edge detection, basic band calcs e.g. (Band_1 - Band_2) that can be strung together for things like NDVI, supervised and unsupervised classification, high/low pass filters etc all seem to be available. Furthermore the basic C++ programs written for these are no more complicated than Python given the level at which the methods are exposed to the user. I am assuming the commercialproducts shine when it comes to more advanced stuff like object recognition and the like. Also I assume they would have the benefit of a helpdesk or support if required. I have used the free trial for eCognition but given that I am usually doing pretty basic stuff I never get a chance to see where it may be better. I am also assuming that it takes a while cutting your teeth on the basics in a program to realise the helpful time savers and get the most benefit out it. Time you don't really get when using a trial version. Is there anything I am missing? Anything that completely stands out that you are thinking 'Wow! this guy is an idiot for not knowing about capability X'? Conversely do you think it is true that most of this stuff is basic technology that has been around for years so it is only the more advanced part of the market pushing the technological envlope that requires these tools?


Geoserver getCapabilties shows wrong OnlineResourse: URL

in my geoserver WMS GetCapabilies document all layeres have an onlineResource parameter which ends with "wms?SERVICE=WMS&" instead of just "wms?" . This has the effect that I can´t cascade my Geoserver layers with other systems as cascading is quite sensitive to irregularities. Can anyone tell if this is normal or is it a bug in geoserver 2.5? Does is correspond to the ogc wms standard? Can I change the GetCapabilities in geoserver per hand?


Color region polygons according to cluster label

I want to present data on a map of our country. I have used R for showing the clusters on a country! Here the steps I took:               setwd('D:/r/cluster2')     channel <- odbcConnectExcel('cluster.xls')     data <- sqlFetch(channel, 'clust9')     y9 <- data.frame(inf=data$infest, faible=data$faible, moyen=data$moyen, fort=data$fort, lon=data$Lon, lat=data$Lat)          library(fossil)          d = earth.dist(y9)       km <- kmeans(d,centers=5)     hc <- hclust(d)         clust <- cutree(hc, k=5 )     set.seed(123)     plot(hc)          y9$clust <- cutree(hc,k=5)     map.AL  <- readOGR(dsn="D:/r/cluster/shp", layer="ALG_boundaries")     map.df  <- fortify(map.AL)          ggplot(map.df)+ geom_path(aes(x=long, y=lat, group=group))+       geom_point(data=y9, aes(x=long, y=lat, color=factor(clust)), size=4)+       scale_color_discrete("Cluster") + coord_fixed()      ![enter image description here](http://i.stack.imgur.com/Geh6I.png) On the map, it's showing just points for the regions that have the same characteristics. However, what I want to show is all the region colored not just a point.


How to represent several points with same coordinates on map?

I'm working on Arcgis 10. I have a point shapefile representing events with several points having coinciding coordinates (They have the same coordinates). For representing them in a dynamic way, I found a nice way - Feature layers maptips. Problem is: I was asked to produce a static (paper) map of the events. Is there a tool that can take these events and shift them apart a bit so all will be seen? (I know that spreading them about will alter their coordinates, but I am prepared to do so on a duplicate shapefile for the display on the map).


varying size of symbol depending on number of datasets

i am trying to map archaeological finds with qgis. some of them have the same coordinates since they are from the same site or very close to each other. is it possible to automatically display the number of finds in a close range through a varying size of the symbol? so, the more datapairs i have in the same spot, the bigger the symbol gets? otherwise 10 finds just appear as one in the map because they are kind of on top of each other...


Dashed lines do not print correct

![enter image description here](http://i.stack.imgur.com/QMF1q.jpg)![enter image description here](http://i.stack.imgur.com/VYkQf.jpg)I have a polygon with a dashed linetype; when viewed in ArcMap the linetype appears correct however when I print to pdf the linetype does not look correct. This might be a pdf quality problem however I am inclined to think that this is an issue within ArcMap.


spatial queries by webgis and virtualizing

i have a question regarding this tutorial: http://www.gistutor.com/postgresqlpostgis/10-intermediate-postgresqlpostgis- tutorials/75-understanding-postgis-spatial-queries.html As u see, it demonstrates how spatial sql statements are working. But my question is, how this tutorial with this webgis was built up. I know that openlayers,json,php,wfs-t,ajax and postgis comes into the play but how all these componentes are working together. Is there any good tutorial out there or do somebody can gives workshops? thank you regards ulan


Splitting lines at a specified distance

Does exist some plugin with it is possible to Splitting lines at a specified distance in Qgis?


OGR Shapefile to CSV: Using Python to get WKT for multipart Polys.

Using ogr2ogr and python, I'm running into a wall with trying to create a csv with WKTs for a shapefile that contains multipart polygons. Currently this is this code I'm using (found it on interweb):               #Purpose: To export a shapefile to WKT          from osgeo import ogr     import sys, os          input = ogr.Open(sys.argv[1])               layer_in = input.GetLayer()     layer_in.ResetReading()     feature_in = layer_in.GetNextFeature()          outfile = open(sys.argv[1] + ".wkt", "w")          while feature_in is not None:              geom = feature_in.GetGeometryRef()         geom_name = geom.GetGeometryName()              outfile.write(str(geom)+ '\n')              feature_in = layer_in.GetNextFeature()      So, do I need to use a different geometry container or do I need to evaluation each feature and if it's a multipart use a different set of geometry classes? I've looked around in the gdal/ogr documentation but I'm having a hard time reconciling it. thanks in advance, matt.


Compute Angle of a line in conventional bearing

I am trying to compute angle between two lines (line segment) using ESRI (ArcGIS, JS) in conventional bearing (Eg: N29'E) I am using Add Toolbar to draw a line segment where point A is a fixed point (already drawn) and when the user is hovering over the map looking for point B, it should display projected angle. Has anybody tried anything like this before or any idea? Thanks!!


How to get the geometry of a polyline to a specific precision?

In python, I use `point = geom.asPolyline()` but it looks like I am not getting any precision just whole numbers? How can I get the geometry and a certain precision?


how to customize the marker using leaflet

I am adding markers to the map as below:               var markers = new L.LayerGroup();     L.marker([Latitude], [Longitude]).bindPopup("Hi there").addTo(markers);     map.addLayer(markers);      I have my own marker icons.Can someone tell me how to add those icons rather than default.....


calculate total Miles in a polygon using Google Map

I have an application that allows the users to create any custom polygon layer on top of google maps, i am interested to know if their is a way to select the road layer in this polygon and then calculate the total miles for these roads.


Reducing number of decimal places after X/Y coordinates in a FC

I'm using arcpy.AddXY_management to add coordinates, and I want to reduce each X/Y field by x number. I've tried to reduce it by string slicing and that obviously didn't work. I also tried setting the precision but that doesn't work on FC's. Anyone have a better idea?


What causes Scripting Error 000714 in ArcGIS when using GetParametersAsText?

I am trying to create a script tool in ArcGIS 10.2 in order to incorporate into a model. The script needs to incorporate three inputs which i've built using GetParameterAsText function then setting up the corresponding links within ArcGIS Script Tool Wizard. However, when I run the script I get **ERROR 000714** (see below). Not sure how to handle this error or why it isn't working to begin with.               import arcpy          inFC = arcpy.GetParameterAsText(0)     outFCName = arcpy.GetParameterAsText(1)     outFolder = arcpy.GetParameterAsText(2)     outFC = outFolder + "\\" + outFCName          try:         if arcpy.Exists(outFC):             arcpy.Delete_management(outFC)         arcpy.CreateFeatureclass_management(outFolder,outFCName,"POLYLINE","#","DISABLED","DISABLED",inFC)         arcpy.AddField_management(outFC,"inFID","LONG","#","#","#","#","NULLABLE","NON_REQUIRED","#")              iCursor = arcpy.da.InsertCursor(outFC, ["inFID","SHAPE@"])              print "Starting!"              with arcpy.da.SearchCursor(inFC,["OID@", "SHAPE@"]) as sCursor:             for row in sCursor:                 inFID = row[0]                 # Print the current multipoint's ID                 #                 # print("Feature {0}:".format(row[0]))                 partnum = 0                      # Step through each part of the feature                 #                 for part in row[1]:                     # Print the part number                     #                     # print("Part {0}:".format(partnum))                          # Step through each vertex in the feature                     #                     prevX = None                     prevY = None                     for pnt in part:                         if pnt:                             # Print x,y coordinates of current point                             #                             # print("{0}, {1}".format(pnt.X, pnt.Y))                             if prevX:                                 array = arcpy.Array([arcpy.Point(prevX, prevY),                                                      arcpy.Point(pnt.X, pnt.Y)])                                 polyline = arcpy.Polyline(array)                                 iCursor.insertRow([inFID,polyline])                             prevX = pnt.X                             prevY = pnt.Y                         else:                             # If pnt is None, this represents an interior ring                             #                             print("Interior Ring:")                     partnum += 1              del iCursor          except Exception as e:         print e.message         arcpy.AddError(e.message)          print "Complete"      Here is the error:               Executing: SplitLines "F:\Extra Data\Prop Frontage\test_for_py.shp" tester.shp "F:\Extra Data\Prop Frontage"     Start Time: Thu Jan 23 11:38:46 2014     Running script SplitLines...     ERROR 000714: Error in script SplitLines.     Error in executing: cmd.exe /C F:\MASTER~2  "F:\Extra Data\Prop Frontage\test_for_py.shp" "tester.shp" "F:\Extra Data\Prop Frontage"          Failed to execute (SplitLines).     Failed at Thu Jan 23 11:38:46 2014 (Elapsed Time: 0.16 seconds)


Is it possible to sync the extents of two open mxds?

I've never seen it, and it's not vital, but I'm wondering if anyone knows of a way to sync the extents of two open mxds? I'm pulling an old ims service into one mxd, and have another mxd open that i will use as the basis of an ArcServer service. I'me trying to match up the layers and symbology, and it would be cool if I could change the extent in one mxd and have it sync with the other. Is it possible? Thanks in advance for any help.


Reprojecting Natural Earth dataset to World Eckert I

My question is connected with QGIS 2. Could anybody explain me how I ought to reproject the Natural Earth dataset there in order my world map to look correctly? When I strive to reproject these data I get the result that you can see on the image below. I was using the next parameters for my projection:               +proj=eck1 +lon_0=120 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs      That is a modified "World Eckert I" projection. ![Figure 1. Natural Earth dataset after reprojecting](http://i.stack.imgur.com/MfGhU.jpg)


How to model a stream in 3-D based on known depths across the stream

I have an excel spreadsheet that contains the elevations across the profile of a stream. I was wondering how to model this in 3-d. Basically there are 6 monitoring stations that are within a stream. I know the elevation of each bank, the slope of each bank, and the depth across the stream at each location. The depth was measured every 20cm across the profile of the stream. eg/ 0.35m across/depth 0.4m, 0.55m across/depth 0.95m, etc. I know that because there are only 6 stations, the result will probably be very blocky, but there are ways to fix that. Basically, I am just looking for somewhere to start. Like I said, I am pretty sure I have all of the pertinent information in an excel spreadsheet, just need to know how to get going. Thanks for all your help.


How do I cluster XY point data? (ArcGIS 9.3)

I have a whole bunch of icons displaying on a map, I want these clustered together as was suggested to me for a different question here... ArcGIS 9.3. How do I increase Picture Marker Symbol performance? A link to an example was provided but I'm not really sure where to begin. I'll preface this by saying I'm a programmer way before I'm a gis person. What I want to know, in total noob terms, is how can I cluster my icons for XY events on the map? Thanks! **EDIT:** Forgot to mention, this is for a C# web ADF GIS portal... If it makes any difference. Is this even possible? I can't seem to find anything for doing this with 9.3 adf.


Openlayers CQL DWITHIN returns errors

I'm trying to use CQL but I guess I do it wrong because I get pink layer in return and error in geoserver : > Binary geometry filter, but first expression is not a property name? (it's a > class org.geotools.filter.LiteralExpressionImpl) > java.lang.IllegalArgumentException: Binary geometry filter, but first > expression is not a property name? This is how I try to do the DWITHIN (get all lines at 200 meter from buildings)               wLayer.mergeNewParams({'CQL_FILTER': "DWITHIN('geom',querySingle('Buildings','geom','INCLUDE'),200,meters)"})      What am I doing wrong ?


Questions about mobile cache, normal tables and basemaps

we are piloting a data gathering application for Windows Mobile and have encountered a few challenges as we don't have previous experience with the ArcGIS Mobile. In future after the pilot, we would need to implement synchronization with a bigger system this application will be part of. 1) We have geodatabase, which contains some standalone tables that have no geometry. We need to be able to update those tables in the mobile application as well. Currently we are trying to create feature classes for each of these tables, add the data into those and then generate some fake geometry point for each row. Is there any other workaround or is this viable solution? 2) Is there any way to add background maps to the application without ArcServer? As I have understood, we need to publish a service in ArcServer and get background maps from the service. Is there any way to have, for example just a .tiff raster map in the device and use that as a background map in the application? This would make our lives much easier during this pilot phase... 3) Are there any known limitations with mobile cache? We have now created the feature classes for few tables and added fake geometries for those. We were able to generate mobile cache with those, but we saw data only for one table, which had only 3 rows. Another table had 14000 rows, and the data was not in mobile cache. The geometry creation had worked because we can see it in ArcMap using the identify tool. 4) Is mobile cache the only solution to save the data? It would be nice to have all the features and tables in one place and it feels such a hassle because we haven't found a way to save normal tables to cache. We are also investigating a way to have two databases, mobile cache and then another database for normal tables. But it feels like there must be easier way around this? We are using ArcGIS Desktop 10 and ArcGIS Mobile 10.1.1 Thanks!


ArcGIS 10 - mosaic: creating 4 files (rasters)

using the Mosaic to new raster tool in order to merge two satellite images, ArcGIS 10 created 4 different files: the coloured image and image1, image2 and image3 (all BW), which are, I guess, three different bands. How to merge them all to have only one file (raster)? (for the purpose of storing the file in the database for the other users, to make it very simple)


Why are the features in my memory layer not displayed?

I am trying to create a memory layer, I followed the directions from the page below but no points appeared on my map. First I tried typing each line into the python console and I also create a plugin with the same code and still nothing appeared on the map. At some point while I was typing the code into the console, it asked for the projection, so I specified my typical projection (epsg 2272, a state plane coordinate system) How to create a new empty vector layer programmatically? I also reviewed http://www.qgis.org/pyqgis-cookbook/vector.html#memory- provider I am experimenting with memory layers because I want to include it into a plugin I am developing.


ArcObjects ITextSymbol.GetTextSize in Python: Why are the out params equal to 0?

This is related to the question I asked earlier about the syntax for this method under Python and comtypes. Now that I know the syntax and am trying to use it, the returned out params are equal to 0. Why is this? Here is my code:               #make TextElement     pTextElement = NewObj(esriCarto.TextElement, esriCarto.ITextElement)     pTextFont = NewObj(stdole.StdFont, stdole.IFontDisp)     pTextFont.Name = 'Arial'     pTextFont.Size = 34.0     pTextSymbol = NewObj(esriDisplay.TextSymbol, esriDisplay.IFormattedTextSymbol)     pTextSymbol.Font = pTextFont     pTextSymbol.HorizontalAlignment = esriDisplay.esriTHACenter     pTextSymbol.VerticalAlignment = esriDisplay.esriTVACenter     txt = 'T>'   #((( hard coded for testing )))     pScrDisp = mxDoc.ActiveView.ScreenDisplay     pDisplay = CType(pScrDisp, esriDisplay.IDisplay)     #Get dimensions of text string     pDisplay.StartDrawing(0, -1)     dblX, dblY = pTextSymbol.GetTextSize(pDisplay.hDC, pDisplay.DisplayTransformation, txt)     pDisplay.FinishDrawing      The out parameters dblX and dblY both equal 0. The hDC and DisplayTransformation inputs are valid. I checked. So, there is something I'm missing. Any insights?


Using model builder, how can I extract certain fields for 1 selected feature from a point shapefile file with multiple features?

I have a point shapefile that is filled with many features. If I have 1 feature selected I need to be able to export certain fields of that selected feature (ID, lat, long) into a dbf file. What tools in model builder can accomplish this?


Creating a new layer from query results

Newbie in arcobject programming. Using C#. Developing a desktop application. I need to join 2 tables (POI layer which is a Point layer and a parcel layer which is a Polygon) where POI.HouseNo=Parcel.UNIT_ID and then create a new shape file which has only the 4 fields POIName,POIType,HouseNo ( all from POI layer attributes), ParcelName ( all from Parcel Layer attributes ) . The user selects the POI Layer and the Parcel layer from combo boxes. Here is the code I have written so far...               IFeatureWorkspace featureWorkspace;     IQueryDef queryDef = featureWorkspace.CreateQueryDef();     queryDef.Tables = cboPOI.Text + "," + cboParcel.Text;     queryDef.SubFields = cboPOI.Text + ".Name," + cboPOI.Text + ".Type," + cboPOI.Text +".HouseNo," + cboParcel.Text + ".ParcelName";     queryDef.WhereClause = cboPOI.Text + ".HouseNo =" + cboParcel.Text+ ".UNIT_ID";      Hope I'm on the right track. Now I need to run this query and get the output to a new shape file. How can I do this? please help.


Creating a quality raster map in QGIS?

I've got 2 raster maps (both with the same values: 0-100). After calculating the average I need to know how precise the map is. Is there any tool in QGis? * * * I try to explain my problem in another way: I have 2 or maybe 3 percentaged raster maps. After putting together ((a+b+c)/3) I need to get a map that shows if the result map is statistically firm. For instance: If all 3 maps have the same value in an area the result must be reliable/firm. If 2 maps have the same values and 1 map is totally different the result is a little bit unreliable. If 3 maps are totally different the result is unreliable. I want to create a colored map (green, orange, red). I already have processed the average map!


How do I allow users to download a GeoTiff on my web map

I am hoping to get help with what I need to do in order to allow a user on my OpenLayers map to download a portion of a GeoTiff they are viewing through a WMS. The data is hosted through GeoServer. I would like the user to be able to simply click a button and download the bounds of the original file I am not sure what steps I need to take in order to begin this process. I have worked a lot with serving geospatial data. But never in a way which allows the user to download the data. Any guide in the right direction would be great. Thanks!


How do I get started with GIS Servers on Amazon EC2?

Are there any tutorials, talks or presentations that discuss the use of EC2 or other cloud computing services to host GIS Servers? Has anyone had any experience in doing it with Amazon's EC2? How hard is it to setup? I know that ESRI has some versions of their server products focused especifically to EC2, has anyone had success in using them? Thank you very much.


Displaying Water Flow Direction Arrows on Roads using ArcGIS for Desktop?

I have road layer and spot heights along the roads. I need to make arrows along the road where flowing water accordingly to the elevation I'm working with ArcGIS 10.1 and 10.2 for Desktop. Roads- Polyline Layer Spot Heights- Point (With Elevation Field) Is there a technique or tool to do this?


Using Space Syntax Analysis plugin for QGIS?

I'm trying to run the Space Syntax Analysis plugin for QGIS. I've added the metadata file that was missing, but it does not appear in the Plugin toolbar when I click on it. Can someone help me to figure out what is going on?


Adobe Flex layer is not displayed, but can zoom to it

I am using Adobe Flex to develop a widget to allow the user to turn on/off 2 layers separately. The strange thing is only 1 layer can be displayed on the map when the layer is turned on, the other layer does not appear when turned on (but can still zoom to the layer). I cannot figure out what happens. The difference of the 2 layers is that the layer not appearing is joined to a database table. Will that affect it? I am using Adobe Flash Builder 4.6 and ArcGIS API for Flex 3.6. Any ideas on it?


Display WMSGetFeatureInfo as tooltips geoExt

I load a map with data and collect information in a json file I want a tooltip appear when mousing over a town on the map I wonder if you have devised. thank you                  var vecLayer = new OpenLayers.Layer.Vector("vector", {             styleMap: styleMap,             protocol: new OpenLayers.Protocol.HTTP({                 // url: "app/data/summits.json",                 //url: "app/data/chargeIndicateur.php",                 url: "../web/app_dev.php/json/indicateur/"+gIndicateur+"/"+gEchelle+"/"+gPeriode,                 format: new OpenLayers.Format.GeoJSON(),                 idProperty: 'Issam'                  }),          , {     // dock: 'right',     xtype: 'toolbar',     id    : 'teste2',          items: [{     text: 'Nom',     labelStyle: 'font-weight:bold;',     }]     } ,      Thanks for quick replay Sorin I ask if you can give me your email


How to create graduated Pie Charts from CSV in QGIS 2.0.1 Defour?

I am trying to make a map in QGIS 2.0.1 ( running on my Mac OS X Maveriks) with the proportionally sized (graduated) pie-charts. I have started with the csv conversion to kml just to place my sample sites on the map, and I am wondering if I can just modify this csv adding some extra columns indicating i) sample size and ii) percentage of each sector for the pie chart. Or there is much less complicated way?


Create A Maximum-Value Composite QGIS

I am looking to create a Maximum-Value Composite raster in QGIS I have three raster images: Oct 2012 Nov 2012 Dec 2012 All are in TIFF Format, Float32, EPSG:4326 - WGS 84 I want to stack these images, then on a pixel-by-pixel basis,examine each value, and retain only the highest value for that pixel location to create a Maximum-Value Composite (refer to http://en.wikipedia.org/wiki/Maximum- value_composite_procedure) e.g.(these are made-up numbers for illustrative purposes): Pixel in Row 205 and Column 106 will have three values: 90 75 100 I want to choose 100, then move to the next pixel and do the same. In the end I should have a final raster with maximum values only, a maximum value composite! I am have QGIS 1.8.0 & 2.0.1, GRASS 6.4.3RC2 and SAGA 2.0.8 at my disposal. Any help will be much appreciated Thanks


Popup is cropped out of window

Sometimes clicking on different features close to window boundary makes part of popup out of window's boundary. ![enter image description here](http://i.stack.imgur.com/R1N37.jpg)               layerInfoWindow = new InfoWindowLite(null, domConstruct.create("div", null, null, map.root));      How can I fix this issue? Thanks


Making a Map with a Legend as a .pdf Without a White Background

I have a shapefile of United States counties. I also have an Excel file of the share of the black population in each county. I'd like to make a map of this using GIS. (Before you say that there are ways to do this outside of GIS, yes I realize that, but I would like to do this in GIS, so please let us not go there.) I am not at all sure how to go about this. Starting backwards: First, I have my shapefile loaded. I know I can save the shapefile as say a .jpg or .pdf under File -> Export Map and then change the Save as Type of .jpg or .pdf. However, when I do so, I get some white in the background from GIS. Is there a way to do this without getting the whtie in the background? Second, how would I include a legend in the map? Third, is there a way I can indicate different shades of one color to mean different things? Like a dark blue is more than 2/3 black, light blue is between 1/3 and 2/3 and white is less than 1/3? Fourth, how do I merge the data on the black population in? I suppose worst case, I can type it in manually by hand. But preferably I wouldn't have to. Is there a way to rearrange the Excel file so that I can merge it in? Like can I merge it in using the fips codes if both files have that? I use ArcMap in ArcGIS 10.1. I very much appreciate any help!


Web map graphics and performance optimization

With a web map (in this case, Bing Ajax API v7), what's more important for performance - the number of individual map shapes/graphics drawn on the map, or the total number of vertices drawn on the map? Put another way, does a single map graphic with 100 vertices use the same amount of resources as 50 map graphics with two vertices each? Thanks.


How to manually edit 16 bit GDEM elevation data on a Geotiff file

I'm building a virtual 3D terrain scenery using both satellite imagery (google maps satellite) and elevation data from GDEM Aster elevation data (geoTIFF format files). My problem is that the elevation data does not perfectly match the satellite images. For example, on some small coastal zones, the sea actually reaches the top of some cliffs... So I'm desperately looking for a way to manually correct the GDEM Aster elevation data (geoTIFF format files) for these zones. I tried using different photo editing software (trial versions and free versions) but none of them lets you set signed or unsigned 16 bit values to individual cells. Even though they are perfectly capable of processing 16 bit greyscale images, they only allow you to set 8 bit values (0-255) in their respective color pickers... I tried : Photoshop CS5, Photoline,, Picture window Pro Paint shop pro 10, Gimp 2 Thank you for any info on any software (Windows) that allows to do this. P.S. I already use Geotiff tools to save and restore the geographic metadata. Best regards


What is the intended Behavior of -skipfailures in ogr2ogr?

I am trying to run ogr2ogr to convert some geojson to shapefile. All is good until any of the geojson features have an invalid geometry (for example, the coordinate array of point geometry is empty). I included the -skipfailures flag, but ogr2ogr quits as soon as it hits the first invalid feature. Am I not doing something wrong here? Example ogr2ogr call:               ogr2ogr -skipfailures -f "ESRI Shapefile" myshape.shp mygeojson.geojson


Performing Raster Noise Reduction and Edge Smoothing?

I have a raster image like this![Raster image](http://i.stack.imgur.com/umDg7.png) I want to delete the isolated noise point and smooth the edge (say the white single points, and I want the edges of the entity look smooth). How can I do this in ArcGIS or in R? The cell size is 30*30.


How to compress enterprise geodatabase with stored procedure?

I'm looking for a stored procedure in ArcSDE for SQL Server which performs geodatabase compression to include it in our daily maintenance plan. Is there any such thing in ArcSDE 10.0?


How do I get started with PostGis

I have been reading a lot about it, but I am struggling to find a good beginner's guide. Which tutorials/books do you recommend? My past experience involves python and arcgis. But I have never worked with python with gis data specifically. I have also worked with mysql before for web development. Thank you very much!


Cell values to No data

![enter image description here](http://i.stack.imgur.com/MU7Rf.png) I need to set cell values in one raster to NoData based on another raster. Both rasters have the same extent and some cells overlap. For those cells that overlap I want to set them to NoData in Raster1 while keeping the rest of the cell values. Please see graphic: Raster1 (green); Raster2 (black). I need to replace the cell values in Raster1 that overlap Raster2 with NoData while keeping the rest of the cell values in Raster1. I'm using the SetNull function in Raster Calculator but I'm getting a blank raster so I guess the expression is not quite right.               SetNull("%Maas_Vol_47_com%","%Maas_Vol_47_inc%")     Raster1 = Maas_Vol_47_com     Raster2 = Maas_Vol_47_inc


Dem file in .adf format, while .tiff is needed

I would like to use this standalone tool to get different kind of Dem visualizations (http://iaps.zrc-sazu.si/en/svf#v). I see that the program accept only .tiff DEM. Now, the DEMs I got are not in that format; rather, I see that the files are .adf I am wondering if I need to convert the DEM I have into .tiff, and how it can be accomplished (possibly in ArcGIS 10.1).


Delete single raster from Raster Dataset

Is there any possibility to delete single raster previously loaded to a Raster Dataset?


How to open osmand obf or pbf files in GIS software?

.obf and .pbf are compressed binary file extensions of OSM data. I want to open these in GIS software like ArcGIS. How is the conversion possible? Is there any automated tool?


Determining the area of each polygon within a polygon boundary?

![enter image description here](http://i.stack.imgur.com/Rwmbd.jpg) I have a series of maps representing a 1000m diameter around a point. These maps contain 40 different types of polygons that represent a land use (ex. Corn crop, Pasture). I want to select all polygons that are within the 1000m boundary and find out each of their areas. I've attempted to create unions and delete the exterior polygons, but I'm only limited to performing a union on 2 features which would make that task very tedious. Is there a faster and simpler way of performing this task? My goal is to create an attbiute table, that will have a unique identifer for each polygon and then a column with the SHAPE_AREA that I can export to excel. I created the maps using ArcGIS, I have access to all licenses and other GIS programs. All polygons are shapefiles, including the boundary circles.


How to obtain buildings base (Ground) Height from Esri Terrain Dataset?

How can I join Height data (surface, DTM) from terrain dataset to buildings polygons ?


QGIS EasyPrint alternatives for v2.0.1?

I recently upgraded to QGIS v2.0.1 and noticed that the EasyPrint plugin is not supported for the latest version. Will this plugin be updated to work in this version, or are there any alternatives for sheet mapping?


Creating a circle with points

I want to arrange 96 points as a circle/ring. Is there a way to make this automatically, for example "Buffer", "vertices"...? Kind regards, Stefan


FME - Calculate distance of point along a line

I have created an FME workspace which inputs a pipeline as well as other linear features which cross it (waterways, roads, railways, fences, other pipelines, etc). The output is a point file which contains the crossing type. What I would like to do is calculate the distance along the pipe from the start of the pipeline for each intersection. So, as you travel along the pipeline, the distance attribute will increase. I'm thinking that I'll need to use some of the linear referencing tools like MeasureGenerator, MeasureExtractor, or something like that, but I'm not sure how to do so. Does anyone have an idea of how to do this?


How to relate Non-Spatial table to Spatial table in two different Postgis DBs?

I need to relate a Non-spatial table to Spatial table, but both are in different database.               Spatial_table(dbname:dist)      gid      district_name      district_code      geom                    Non_Spatial(dbname:census)      ID      district_code      population      male_popu      female_popu      Can anyone please suggest me, How to relate the above to tables to get the query result for the population of specific district? Also can anyone tell me about the difference between Joining and Relating of two tables. Thanks in advance..


Data about administrative subdivisions capital cities

I am looking for data about **second or third rank administrative capitals** in **every African country** (the administrative center of a country territorial subdivision), but I can't find anything close to what I'm looking for. So far, I found data on administrative boundaries (GADM, Natural Earth, UN SALB Project, OSM) and about cities (Natural Earth, OSM, Geonames), but nothing that specifically links an administrative subdivision to its administrative capital. Although Geonames is supposed to provide that kind of information, I briefly checked and found that the data is of very poor quality. Can anyone tell me if what I'm looking for exists? Thank you for your help!


How to resize symbols in the QGIS Print Composer legend?

I am trying to map a count data type via graduated symbol against a chloropleth and the legend is too big. The circles of the graduated symbols are the actual size of their corresponding values in the map, which is good for interpreting the legend but is not good because they overlap and in order for space to be created between them, the legend stretches absurdly along the y axis. Is there any way to scale the circles back to half the size of the mapped symbols? I stretched out the legend as you can see below, aesthetically it wasn't ideal. ![this is what my map looked like when i stretched out the legend, it was rather ridiculous looking](http://i.stack.imgur.com/bd94y.png) got the answer! ![enter image description here](http://i.stack.imgur.com/Fb5O6.png)


Parsing raw data from GPS tracker

As a side project I've been hacking together a simple tracking system for my wayward dog using a cheap GPS / GPRS tracker. I've successfully got it transmitting data to my server over GPRS but while the tracker display shows the correct LAT / LNG the data sent via. GPRS is out (by 8° LAT and -0.37° LNG) so I'm guessing I need to do some post processing? Raw data '#35682303151XXXX#uid1#0#9999#AUT#4#V#00402.6244,W,5556.4382,N,000.38,214#110913#165318#V#00402.6248,W,5556.4386,N,000.22,41#110913#165338#V#00402.6259,W,5556.4371,N,000.32,41#110913#165358#V#00402.6266,W,5556.4364,N,000.05,41#110913#165418##' which has 4 fixes at 20s intervals which I've parsed out as; [lat] => -4.026244 [lng] => 55.564382 [date] => 2013-10-15 16:53:18 [lat] => -4.026248 [lng] => 55.564386 [date] => 2013-10-15 16:53:38 [lat] => -4.026259 [lng] => 55.564371 [date] => 2013-10-15 16:53:58 [lat] => -4.026266 [lng] => 55.564364 [date] => 2013-10-15 16:54:18


Leaflet - Change style on click

I am beginer with Leaflet library and I am facing an issue with a geojson. I have an extern geojson file (points) and I would like to display a circle for each point coordinates and get a bigger circle on a mouse click, with a popup. I have no problem to load the geojson and to view the circle markers, but I can't get both bigger cricle and popup. Here is part of my code :               function clickFeature(e) {                         if(clickedFeature) {clickedFeature.setRadius(3)}                         var layer = e.target;                         e.target.setRadius(5);                         clickedFeature = e.target;                         info.update(layer.feature.properties)                         }        function onEachFeature(feature, layer) {                                         layer.bindPopup(feature.properties.PROPERTY).addTo(MyLayer);                                         layer.on({click: clickFeature,})                                         }          MyGeojson = new L.geoJson(mygeojson, {         onEachFeature: onEachFeature,                pointToLayer: function (feature, latlng) {         return L.circleMarker(latlng, {                     radius: 3,                     fillColor: getColor(feature.properties.TYPE),                     color: "#242121",                     opacity: 1,                     weight: 1.5,                     fillOpacity: 1                 });             }         }).addTo(map);      First I had this same kind of functions and code, but using PNG icon instead of circleMarkers, and it worked well, both bigger icon and popup on mouse click. Any idea what is wrong ?


How to set up multiple Content.ftl files in Geoserver?

So today i started to play around with the content.ftl file for geoserver to change around the popup style. In playing i've discovered that if i call in my content.ftl file which contains code for a layer that isn't in my openlayers html then the popup doesn't work. My question is; Can i set up multiple content.ftl files for different openlayers projects which will contain different layers? If so, how? Example - if i Preview a layer in GeoServer for say planning applications but my content.ftl contains details for styling another of my geoserver layers as well as the planning application style, the popup doesn't appear. Remove the styling for the other layer and the pop up works! Any help would be great. Med I worked it out and heres how.... Instead of having a generic content.ftl file located in the template folder in data_dir, i moved the content.ftl into data_dir/workshop/.....to the layer folder and placed it here. The folder also contained featuretype and layer xml files for reference. So now each of my layers has its own content.ftl and popup style.


Set region around multiple polygons

I have combined multiple objects per row in my layer. Is there a way to build another layer that surrounds the combined objects, per row?See example. For Q-GIS, thx.![enter image description here](http://i.stack.imgur.com/ILfZu.gif)


Select only records with region objects in MapBasic

I am writing a program that finds overlapping objects. I can't perform the AreaOverlap() function on objects that don't enclose an area (e.g. text, points) as the result would be meaningless. **Is there a way to select only records that have objects that are regions?** I have a function that checks whether an object is a region (see below) so I could loop through the objects and check them one-by-one but it would be more efficient to do this as part of a select statement. **My example function:**               Function IsObjectRegion(ObjAlias As Object) As Logical       Do Case ObjectInfo(ObjAlias, OBJ_INFO_TYPE)         Case OBJ_TYPE_REGION,               OBJ_TYPE_RECT,               OBJ_TYPE_ROUNDRECT,               OBJ_TYPE_ELLIPSE           IsObjectRegion = True         Case Else           IsObjectRegion = False       End Case     End Function      **Update: my final, working code using Peter's method and my function:**               Include "mapbasic.def"          Declare Sub main     Declare Function IsObjectRegion(ByVal ObjAlias As Object) As Logical          Function IsObjectRegion(ByVal ObjAlias As Object) As Logical       Do Case ObjectInfo(ObjAlias, OBJ_INFO_TYPE)         Case OBJ_TYPE_REGION,               OBJ_TYPE_RECT,               OBJ_TYPE_ROUNDRECT,               OBJ_TYPE_ELLIPSE           IsObjectRegion = True         Case Else           IsObjectRegion = False       End Case     End Function          Sub Main       Select * from TargetTable Where obj into temp1 noselect       Select * from temp where IsObjectRegion(OBJ) into results noselect       Browse * from results       Close table temp     End Sub


How to run an external exe in ArcGIS Python?

I need help with one question. I would like to run an external executable written in fortran which provides one of the inputs inside a python code. Is it possible to run that exe inside python script? Thanks, Ibe


Loading multiple feature types as one layer in OpenLayers

I have a kml file which contains some polygon and some point features in the same file. I am trying to load and view this layer on my basemap in OpenLayers. The problem is that the polygon features will be drawn on map but not the points. Does this mean that a layer cannot containt both polygons and points in it at the same time? What sollution would you suggest to solve this? (I need these features in one layer and not seperate). Thank you. The code for loading the data in OpenLayers:               Grid = new OpenLayers.Layer.Vector("grid layer", { eventListeners: {                 'featureselected': function (evt) {},                 'featureunselected': function (evt) {},                 strategies: [new OpenLayers.Strategy.Fixed()],                 protocol: new OpenLayers.Protocol.HTTP({                     url: "sample.kml",                     format: new OpenLayers.Format.KML({                         extractStyles: false,                         extractAttributes: true,                         maxDepth: 2                     })                 }),                     styleMap: new OpenLayers.StyleMap({                     "default": new OpenLayers.Style(OpenLayers.Util.applyDefaults({                         pointRadius: 6,                                             strokeWidth: 1,                         strokeColor: "Black",                         fillColor: "Gray",                         fillOpacity: 0.4                     }, OpenLayers.Feature.Vector.style["default"])),                     "select": new OpenLayers.Style({                         fillColor: "Blue",                         fillOpacity: 1                     })                 })             });     map.addLayer(Grid);


Simplest way to limit the memory that the raster package uses in R

What's the easiest way to limit the memory the `raster` package uses in `R` to 5GB of RAM? Ideally I want any raster that exceeds this size to be written to the hard drive and analysed from a temporary file. I'm using Windows 7 (64Bit) and R version 3.0.3. Here's some code I've been playing with on a machine with 16GB of RAM but it keeps maxing out as i'm dealing with large rasters:               require(raster)          memory.limit(size = 5000)     rasterOptions(format="GTiff",overwrite=TRUE,datatype="INT1S",                   tmpdir="C:/Research/BIN",tmptime=1.1,progress="text",chunksize=1000,                   maxmemory=1000)     rasterTmpFile("delete_me_")          r <- raster("myraster.tif")     r[r==0] <- NA      It is at this stage in the script that the RAM usage increases to take all 16GB and then stops working and cancels the function. What do I need to change here? Any help much appreciated.


Online Interactive Directory

I have the following project and I'm figuring out how to best implement it. I would like the project to be as highly automated as possible. The suggestion was made that I make and access table, link it to ArcGIS and then share the directory/map created through sharepoint. Objective: Put together a list of facilities (company owned, company leased, government leased) on a map. Criteria:   * Develop a system that can be easily updated by non-GIS users.   * Background data can be stored in an easily updated program (such as Microsoft Access) and dynamically linked to the mapping function.   * Data sets should be accessible to various people in the company (via something like sharepoint).   * Mapping function should allow for layers of data My main question is the use of access and share point the best way of going about the project? Also lets say I wanted to link my map to an online database (Orcale based) that's protected by a login screen is there any way to do that?


How to set a coordinate transformation GRID

Does anyone know how to set a coordinate transformation in qgis so I can measure the distances in meters or kilometers? If anyone has a clue, would be just great. Thanks!


How to integrate osgeo4w to main python installation?

I'm not sure if I can do this but I would like to access methods from osgeo4w from my main python27 installation. How do I accomplish this or are you only able to import osgeo4w from the python console within QGIS?


How to cut/split multipolygon with blade (linestring) into two pieces using PostGIS?

I struggle to do split operation in postgis I tried to build a query but Im always fail. Id like to cut a multipoligon with a blade (a linestring) into two pieces. These two objects have specific SRID of my country. query:               CREATE TABLE Cut_result AS     SELECT ST_SPLIT(poligon,blade)     FROM (SELECT     ST_MakeLine(700264 126376, 780844 122481)::geometry(LINESTRING,23700) AS blade),     SELECT FROM county AS poligon;      ERROR               ERROR:  syntax error at or near "126376"     LINE 5: ST_MakeLine(700264 126376, 780844 122481)::geometry(LINESTRI...                                ^     ********** Error **********


Is there a way to set different levels of transparency (opacity) for a polygon's border and its interior?

I have some polygons which I'd like to display with their interior displayed at a higher transparency than their borders. Right now, I am doing it by creating two separate layers with the same data source, one for the border (20% transparent) and one for the interior (70% transparent). Actually, I usually end up creating 3 different layers in these cases: one each for the polygon's border and fill, and a third to use in the legend, because neither of the other two will work perfectly in the legend. Here's an example, in this project I have floodplains set to 70% transparent in the fill (because many other items I'm interested in are in the floodplain and I need to see them), and 0% transparent on the border so that its boundaries are clearly visible. Is there a way to set the transparency of the polygon's border separately from that of the interior, so I don't have to duplicate layers? ![enter image description here](http://i.stack.imgur.com/A8ecn.png)


Why are Character Marker Symbol Dialog not showing up in the Symbol Property Editor dialog?

I just installed ArcGIS 10.2.1 on a brand new machine and for some reason Character Marker Symbols are not loading into the Symbol Property Editor. Is anyone experiencing the same issue? Is there a solution for this? ![enter image description here](http://i.stack.imgur.com/qtAjj.png)


split features (vector line) problem

Part of a stream is now underground so I'm trying to use split features to break the stream into two pieces. The manual says to select the feature, then the split features icon, then draw a line across the stream where I want the break. All it does is keep drawing red lines, with vertices wherever I left click, green vertices where I right click, until I crash qgis. (Yes I toggled edit) What should-be-obvious thing am I missing?


Where to license non-automotive POIs for France, Spain and Portugal?

qualified POIs are needed for a commercial project - should cover broad out- of-car usecases ( _not_ needed: parking, gas stations). Could be commercially licensable data or derived from national geodata portals and INSPIRE related initiatives. Or something completely different.


Create custom pixel block size in ArcGIS 10?

I'm using the new IRawBlocks in ArcGIS 10 from the adopted sample in ArcGIS 10 SDK help and I want to create a custom sized pixel block that feeds me a block of values column by column instead of a block. So basically I'll be looping through each column in the raster. Any ideas as I couldn't find anything in the documentation?                       IRawBlocks rwBlocks = (IRawBlocks)rb;             IRasterInfo rsInfo = rwBlocks.RasterInfo;             IPixelBlock pb = rwBlocks.CreatePixelBlock();                  // change pixel blocks size here?????                       // iterate through pixel blocks             for (int pbXcursor = 0; pbXcursor < rsProps.Width; pbXcursor++)             {                 string line = "";                 rwBlocks.ReadBlock(pbXcursor, 0, 0, pb);                      // iterate through cells in pixel block, should only be one pixel column                 for (int pbHeight = 0; pbHeight < pb.Height; pbHeight++)                 {                     object val = pb.GetVal(0, 0, pbHeight);//pbWidth, pbHeight);                     // etc                 }             }      Thanks Rob


How to get radiobutton checked after pop-out window?

I am trying to migrate my program from ArcMap 10 VBA to ArcObjects VB.net add- in in Visual studio express 2008. I write a ArcGIS addin component tool to draw a geometry, perform a spatial query and then pop out a windows form with radio buttons to decide based on which layers in TOC to perform the query. Now what I have encountered is in OnMouseDown Event the radiobutton read from current form and return no value. when I execute the code it opens the form as i expected but doesn't allow me to select radiobutton, besides it doesnot show any error or relevant warnings in VS express errorlist. Here are the codes in Tool Class               Public Class Tool1     Inherits ESRI.ArcGIS.Desktop.AddIns.Tool     Dim frm2 As New Form2          Public Sub New()          End Sub          Protected Overrides Sub OnUpdate()         Enabled = My.ArcMap.Application IsNot Nothing          End Sub          Protected Overrides Sub OnMouseDown(ByVal arg As ESRI.ArcGIS.Desktop.AddIns.Tool.MouseEventArgs)                   ' MyBase.OnMouseDown(arg)         Dim pMxDoc As IMxDocument         Dim pMap As IMap         Dim pActiveView As IActiveView         pMxDoc = My.ArcMap.Document         pActiveView = pMxDoc.FocusMap         pMap = pMxDoc.FocusMap                  Call AddDataFromFile.AddShapeFile("\\Hkhkgsv02\GIS\Basemap\Database_Index", "B1k_index")             Call AddDataFromFile.AddShapeFile("\\Hkhkgsv02\GIS\Basemap\Database_Index", "B5k_index")             Call AddDataFromFile.AddShapeFile("\\Hkhkgsv02\GIS\Basemap\Database_Index", "B10k_index")             Call AddDataFromFile.AddShapeFile("\\Hkhkgsv02\GIS\Basemap\Database_Index", "B20k_index")              'Determine Selection Layer based on which scale         Dim pScreenDisplay As IScreenDisplay         Dim pRubberEnv As IRubberBand              pScreenDisplay = pActiveView.ScreenDisplay         pRubberEnv = New RubberEnvelope               Dim pEnvelope As IEnvelope         pEnvelope = pRubberEnv.TrackNew(pScreenDisplay, Nothing)         'Rectangle Display              If pEnvelope.IsEmpty = True Then Exit Sub         With pScreenDisplay             .StartDrawing(pScreenDisplay.hDC, ESRI.ArcGIS.Display.esriScreenCache.esriNoScreenCache)             .FinishDrawing()         End With         Dim pGeometry As IGeometry         pGeometry = pEnvelope              'Build a Spatial Query Filter              Dim pSFilter As ISpatialFilter         pSFilter = New SpatialFilter         pSFilter.Geometry = pGeometry         pSFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelContains             frm2.Show()                       If frm2.RadioButton1.Checked = True Then                pFlayer = pMxDoc.FocusMap.Layer(FindLayerIndex.Find_Layer_Index("B1k"))             ElseIf frm2.RadioButton2.Checked = True Then                 pFlayer = pMxDoc.FocusMap.Layer(FindLayerIndex.Find_Layer_Index("B5k"))             ElseIf frm2.RadioButton3.Checked = True Then                 pFlayer = pMxDoc.FocusMap.Layer(FindLayerIndex.Find_Layer_Index("B10k"))             ElseIf frm2.RadioButton3.Checked = True Then                 pFlayer = pMxDoc.FocusMap.Layer(FindLayerIndex.Find_Layer_Index("B20k"))             End If          Dim pSelectionLayer As IFeatureSelection         pSelectionLayer = pFlayer         If pFlayer Is Nothing Then Exit Sub         Dim pFeatureClass As IFeatureClass         Dim pFCursor As IFeatureCursor         pFeatureClass = pFlayer.FeatureClass         pFCursor = pFeatureClass.Search(pSFilter, True)              ' Set Array to Store Feature Value         Dim NameValueArray() As String         ReDim NameValueArray(pFeatureClass.FeatureCount(pSFilter))                   'Find Fieldindex to read Tile's Map Name              Dim intFieldIndex As Integer         intFieldIndex = pFeatureClass.FindField("TILE_NAME")     ...          End Sub               End Class      Hope I can get some clue.thanks. Jing


How to cast box3d_extent to geometry (Rare problem with Geoserver and Postgis)

I use Geoserver 2.1.3 and "POSTGIS="1.5.3" GEOS="3.2.2-CAPI-1.6.2" PROJ="Rel. 4.7.1, 23 September 2009" LIBXML="2.7.6" USE_STATS" on PostgreSQL 9.1 From time to time (well - today I can't get around it anymore) I get the following error when I try to "calc native bounds from data" while creating a new layer in Geoserver.                   java.lang.RuntimeException: java.io.IOException: Error occured calculating bounds at      org.geotools.jdbc.JDBCFeatureSource.getBoundsInternal(JDBCFeatureSource.java:485) at      org.geotools.jdbc.JDBCFeatureStore.getBoundsInternal(JDBCFeatureStore.java:179) at      org.geotools.data.store.ContentFeatureSource.getBounds(ContentFeatureSource.java:370) at          ...         Caused by: java.io.IOException: Error occured calculating bounds at          org.geotools.jdbc.JDBCDataStore.getBounds(JDBCDataStore.java:1137) at          org.geotools.jdbc.JDBCFeatureSource.getBoundsInternal(JDBCFeatureSource.java:478) ...          94 more          Caused by: org.postgresql.util.PSQLException: ERROR: function st_envelope(box3d_extent) does not exist          No function matches the given name and argument types. You might need to add explicit type casts. Position: 30      The geometry column contains points, the SRID is defined in the *geometry_columns* tables. And it worked before... Anyways.. I will upgrade to PostGIS2 in a few weeks.. For now I need help on how to transform a PostGIS **box3d_extent** to a **geometry**. With that information I want to dirty-quick-fix this with the creation of a local hack-function **st_envelope(box3d_extent)** which would then call **st_envelope(geometry)**. Thanks in advance for any answers helping with the cast- and/or base problem ;-) Steve


Delete first two letters from a string field

I'm looking for a field calculator python statement to delete only the first two letters of string entries. For Ex:               !Name!.strip([0:2])      I know this syntax is incorrect. Thanks a lot


PostGIS Raster in MapServer

I'm having trouble getting a PostGIS raster to load and display in MapServer. I'm not sure, but suspect it's a version problem. My mapfile layer and the error message are shown below.               LAYER         NAME bio12_2_5m         TYPE raster         STATUS ON         DATA "PG:host=localhost port=31600 dbname='wdata' user='postgres' password='*****' schema='worldclim' table='bio12_2_5m'"           PROCESSING "NODATA=-9999"         PROCESSING "SCALE=AUTO"     END      msDrawMap(): Image handling error. Failed to draw layer named 'bio12_2_5m'.   msDrawRasterLayerLow(): Unable to access file. Corrupt, empty or missing file 'PG:host=localhost port=31600 dbname='wdata' user='postgres' password=' _*_ *' schema='worldclim' table='bio12_2_5m'' for layer 'bio12_2_5m'. Sorry, but table bio12_2_5m doesn't seem to have regular blocking arrangement. Only tables with regular blocking arrangement can be read from now This is WorldClim data loaded into PostGIS with raster2pgsql like:               raster2pgsql -s 4326 bio12.bil worldclim.bio12_2_5m > bio12.sql      To get around what seems like a blocking problem I've also tried:               raster2pgsql -s 4326 -C -r bio12.bil worldclim.test > test.sql     raster2pgsql -s 4326 -t 10x10 bio12.bil worldclim.test2 > test2.sql      And with the MapServer layer I've tried setting mode=1 or mode=2, still no success. Versions ... Ubuntu 12.04 LTS, PostgreSQL 9.1.11, GDAL 1.10.0, PostGIS 2.0.4 with raster support, and MapServer 6.0.1. But (always a but) to be fair I started with GDAL 1.7.3 via the default ubuntu repository, and at the end updated GDAL to 1.10.0 via the Ubuntu GIS "unstable" repository since it seems like the newer versions of GDAL are required to address/support the blocking items. Thanks, very close I suspect and not sure if installing everything from scratch would help starting with the 1.10.0 GDAL. Also tried building GDAL from source and no success there.


Where can I find a shapefile of city boundaries in the US?

Where can I find a shapefile that contains boundaries for the major cities in the United States? I have tried the census website but could not find anything.


Zoom to selected feature reducing MapFramesize using ArcGIS Engine and C#?

I am trying use the the following code for zooming to a selected feature using ArcEngine 10.2.1 and c#. It is performing, **but it is reducing the Mapframe to the feature size**. Which was not happening when i have used the same code using ArcEngine 9.3.               IFeatureSelection pFselection;     IQueryFilter pQueryFilter = new QueryFilter();     pFselection = mainfrm.axPageLayoutControl1.ActiveView.FocusMap.get_Layer(0) as IFeatureSelection;                   pQueryFilter.WhereClause = "STATE = '" + cbxStates.Text.ToUpper() + "'";     pFselection.Clear();     pFselection.SelectFeatures(pQueryFilter, ESRI.ArcGIS.Carto.esriSelectionResultEnum.esriSelectionResultNew,false);     IEnumFeature pEnumFeat;     IFeature pFeature;     pEnumFeat = mainfrm.axPageLayoutControl1.ActiveView.FocusMap.FeatureSelection as IEnumFeature;     pEnumFeat.Reset();     pFeature = pEnumFeat.Next();     pFeatureEnvelope = pFeature.Shape.Envelope;pMapFrame.Map = mainfrm.axPageLayoutControl1.ActiveView.FocusMap;     pMapFrame.ExtentType = ESRI.ArcGIS.Carto.esriExtentTypeEnum.esriExtentBounds;     pMapFrame.MapBounds = pFeatureEnvelope;     mainfrm.axPageLayoutControl1.ActiveView.Extent = pFeatureEnvelope;      can any one suggest the solution.


In ArcMap, [how] can you set a layout element's "name" property?

Map Layout elements appear to have a "name" property that can be used to identify them. In ArcMap is it possible to set this property (interactively), and if so how is this done? I'm using ArcGIS 10.1, and the only place I can think to look is on a map layout element's "Properties..." options (when you right click on an element) however nothing obvious is present unless there is some pseudo-code that you can enter, similar to how dynamic text works? From the help page (ListLayoutElements, in case the link doesn't work): "It is the map document author's responsibility to ensure each page element is given a unique name so that elements can be uniquely identified." So how would a map author do this? Please note, I'm not interested in an arcpy/python solutions (that seems straight forward) as I want this to be applicable for non-python users.


Group Nearest neighbor query in SQL Server 2012 spatial database

I am looking for way to perform GNN query over road networks in SQL Server. Any pointers


WMC as a permalink

I'm using the permalink control in openlayers but there are limitations in the way it handles map's context. For example if users load an external WMS layer it won't appear in the permalink. WMC has less limitations, but it is a XML file so you can't bookmark it and it is less easy to share for general public. The best of both worlds would be a transformation of the WMC file to fit a url. I guess the whole transformation would be compression, serialization, maybe XML to JSON... Is it feasible with limitation of URL's length? Has anyone tried this?


Calculating Lag tolerance and Lag width for experimental variogram in kriging

I need to make a formula for calculating the lag width and lag tolerance of the experimental variogram in kriging. I am programming in asp c#:               double Max_X_Point=0,Y=0,Max_distance = 0,Lagdistance;     for (int i=0;i<KnownPoints.Count;i++)     {     if (Max_X_Point < KnownPoints[i].X)     {     Max_X_Point = KnownPoints[i].X;     Y=KnownPoints[i].Y;     }     }     MaxLagDistance = distVincenty(0, 0, Max_X_Point, Y);     No_Lags =Convert.ToDouble(TextBox7.Text);     for (int i = 0; i < n; i++)     {     if (Max_distance < Variograms[i, 0]) Max_distance = Variograms[i, 0];     }      LagWidth = MaxLagDistance/No_Lags;     while (LagWidth * No_Lags < Max_distance)     {     MaxLagDistance += 1;     LagWidth = MaxLagDistance / No_Lags;     }     Lagdistance = LagWidth;      I try to calculate the total distance of the variogram from (0,0) to a point that has the maximum x_axis coordinate. And check if the maximum lag distance is smaller than the maximum distance in my calculations to ensure embedding all pair points that I calculate the lag with. Please, can any one with experience in variogram calculations or in kriging help me to understand the calculations of lag width and lag tolerance to correct my calculations.


Error importing NOAA’s National Geophysical Data Center as delimited text layer in QGIS version 1.02

I've been struggling with trying to import the most recent data from NOAA’s National Geophysical Data Center earthquake information into Quantum GIS 1.02 for weeks now for a class assignment. I'm attaching three screen shots which show some of what I'm running into. The first is where I import the information in a delimited text layer. ![enter image description here](http://i.stack.imgur.com/sX6Ay.jpg) The second is the error I continue to run into every time I load the file which I've re-downloaded and attempted to use multiple times with the same error. ![enter image description here](http://i.stack.imgur.com/63ncv.jpg) Lastly, is where I'm attempting to make use of the Points in Polygone feature which will load to 100% and then freeze the program. ![enter image description here](http://i.stack.imgur.com/4n7PI.jpg) Any advice or suggestions?


Using GRASS plugin tools from OSGEO4W command line

I am trying to do the interpolation of GRASS vector shapefiles using GRASS plugin tools. My objective is to run interpolation tool in GRASS plugin from OSGEO4W command line, which will help me automate the process in my project. I have found an article on how to do this, but it seems to work for QGIS tools and I don't know how to access the GRASS tools in a manner similar to that mentioned in this website. The website URL is as below: http://nathanw.net/2011/09/27/generating-contours-using-gdal-via-shell-or- qgis/ So if anyone has an idea as to how to do this, please help me out with this.


custom time interval for geowebcache

How to assign time interval to take regular cache using geowebche. For instance every afternoon at 13:00 hrs geowebcache should be started automatically and it should start creating tiles. ?


How to use Arcpy UpdateCursor to update a field?

This should be simple but I've spent days trying everything to get it to work. I need to update a field in a layer based on two other fields. I can't even get it to update based on one field. Could someone please tell me whats wrong. So in this script I want to be classifying the outfield as "poor" if the infield is "Sandy". I will also have another field of numbers (and infield2 = "10").               rows = arcpy.UpdateCursor(infc) #this is my feature layer     for row in rows:         if row.getValue(Infield1) == "Sandy":             row.setValue(Outfield, "poor")         else:             row.setValue(Outfield, "none")         rows.updateRow(row)     del rows      Thanks for all the answers. The script only works if written like this: I did not want to name the field names though, but at least it works. The only way I could do a range 5 - 10 was to put > 5 and <10 (see line 23)               for row in rows:         if row.getValue("Texture") == "Sandy" and row.getValue("Humus") < 10:             row.setValue("Richness", "poor")         elif row.getValue("Texture") == "Sandy" and row.getValue("Humus") > 10:             row.setValue("Richness", "medium")         elif row.getValue("Texture") == "Sandy-Clay" and row.getValue("Humus") < 5:             row.setValue("Richness", "poor")         elif row.getValue("Texture") == "Sandy-Clay" and row.getValue("Humus")  >= 5 and row.getValue("Humus")  <= 10:             row.setValue("Richness", "medium")         elif row.getValue("Texture") == "Sandy-Clay" and row.getValue("Humus") > 10:             row.setValue("Richness", "rich")         elif row.getValue("Texture") == "Clay" and row.getValue("Humus") < 5:             row.setValue("Richness", "poor")         elif row.getValue("Texture") == "Clay" and row.getValue("Humus") > 5:             row.setValue("Richness", "rich")         else:             row.setValue("Richness", "unclassified")         rows.updateRow(row)     del rows, row


IDatasetContainer2 CreateDataset method fails after upgrading to ArcObjects 10

I have a Network Analyst routine that worked for 9.3 and does not work with 10. I was closely following the example at: http://help.arcgis.com/en/sdk/10.0/arcobjects_net/conceptualhelp/index.html#/How_to_create_a_network_dataset/0001000000w7000000/ It bombs when it gets to the line that actually creates the Network Dataset:               INetworkDataset networkDataset = (INetworkDataset)datasetContainer2.CreateDataset(deDataset);      Did something change that would cause these objects (used with Network Analyst) to not work the same way? This is the exact error message:               System.Runtime.InteropServices.COMException was unhandled Message=Exception from HRESULT:      0x8004228E Source=ESRI.ArcGIS.Geodatabase ErrorCode=-2147212658 StackTrace: at      ESRI.ArcGIS.Geodatabase.IDatasetContainer2.CreateDataset(IDEDataset DataElement) at      NetAnalystTestApp.MainForm.PutRoadLayerOnMap() in C:\@ALPACAS\AlpacasCT\experimental\network      routing\NetAnalystTestApp\NetAnalystTestApp\MainForm.cs:line 352 at      NetAnalystTestApp.MainForm.MainForm_Load(Object sender, EventArgs e) in      C:\@ALPACAS\AlpacasCT\experimental\network      routing\NetAnalystTestApp\NetAnalystTestApp\MainForm.cs:line 83 at      System.Windows.Forms.Form.OnLoad(EventArgs e) at      System.Windows.Forms.Form.SetVisibleCore(Boolean value) at      System.Windows.Forms.Control.set_Visible(Boolean value) at      System.Windows.Forms.Application.ThreadContext.RunMessageLoopInner(Int32 reason,      ApplicationContext context) at      System.Windows.Forms.Application.ThreadContext.RunMessageLoop(Int32 reason,      ApplicationContext context) at System.Windows.Forms.Application.Run(Form mainForm) at      NetAnalystTestApp.Program.Main() in C:\@ALPACAS\AlpacasCT\experimental\network      routing\NetAnalystTestApp\NetAnalystTestApp\Program.cs:line 39 at      System.AppDomain._nExecuteAssembly(RuntimeAssembly assembly, String[] args) at      System.AppDomain.ExecuteAssembly(String assemblyFile, Evidence assemblySecurity, String[]      args) at Microsoft.VisualStudio.HostingProcess.HostProc.RunUsersAssembly() at      System.Threading.ThreadHelper.ThreadStart_Context(Object state) at      System.Threading.ExecutionContext.Run(ExecutionContext executionContext, ContextCallback      callback, Object state, Boolean ignoreSyncCtx) at      System.Threading.ExecutionContext.Run(ExecutionContext executionContext, ContextCallback      callback, Object state) at System.Threading.ThreadHelper.ThreadStart() InnerException:      Thanks, Renee


Why would Openlayers always render over other divs?

I'm having trouble getting an Openlayers div to allow other divs to appear above it. Having tried the solution proposed here as well as the more classy-looking one here, I haven't had any luck yet. Openlayers always orders itself over other divs. I've also set the z-index for my other div (which I want on top) as high as it goes, still no effect. Is there somewhere it could be getting overridden? Or some setting that needs an extra trick to get rid of? I have asked this question on StackOverflow as well, where it's arguably more relevant. I'm only posting this here because someone here may have specific Openlayers experience on this matter.


Update geometry(polygon) of shape file based on a point layer

I have two shape files.   1. Polygon shape file   2. Point file contains some points. Now i want to delete all the vertex of polygon shape file which are available in point file. I am using ArcGIS(ArcEditor) platform. Any ideas.


Available iOS GIS applications/framworks for mobile data accquisition

I'm starting to develop an application for iPad which should enable the users to go to the field and acquire georeferenced information. The idea is to enable acquisition of points and polygons, and the filling of a set of attributes for the geometry. Also important is the ability to connect to externals geo-referenced databases, such as a PostGIS one. Ideally, I would want to be somewhat similar to the ArcGis for iOS application, which has a good user interface. Can you guys point me to other applications similar to what I am developing? I already know the following ones:   * Geo-jot   * pcMapper   * iGIS   * GISRoam   * field_assets   * Geopaparazzi Possibly you can point me to something with similar functionality to the ArcGis for iOS? OR Interesting programming frameworks to achieve the desired functionality on iOS? Thank you.


Error when pointing Geoserver to PostGIS

I have a PostGIS database that contains all my GIS data along with metadata. I have it published as a Store in GeoServer. They are both installed on the same windows server 2008r2. When i try to add layers from the store resource and publish it, i encounter an error at the bottom of the "Data" tab of the resource configuration page. Any ideas what might be my problem? ![enter image description here](http://i.stack.imgur.com/090IC.png) When i click on the yellow warning icon i get: ![enter image description here](http://i.stack.imgur.com/0sdBe.png) Here is the error:               org.apache.wicket.WicketRuntimeException: Exception in rendering component: [Component id = message] at org.apache.wicket.Component.renderComponent(Component.java:2563) at org.apache.wicket.markup.html.WebComponent.onRender(WebComponent.java:60) at org.apache.wicket.Component.render(Component.java:2361) at org.apache.wicket.MarkupContainer.renderNext(MarkupContainer.java:1298) at org.apache.wicket.MarkupContainer.renderAll(MarkupContainer.java:1428) at org.apache.wicket.Page.onRender(Page.java:1502) at org.apache.wicket.Component.render(Component.java:2361) at org.apache.wicket.Page.renderPage(Page.java:913) at org.apache.wicket.request.target.component.PageRequestTarget.respond(PageRequestTarget.java:63) at org.apache.wicket.request.AbstractRequestCycleProcessor.respond(AbstractRequestCycleProcessor.java:104) at org.apache.wicket.RequestCycle.processEventsAndRespond(RequestCycle.java:1245) at org.apache.wicket.RequestCycle.step(RequestCycle.java:1316) at org.apache.wicket.RequestCycle.steps(RequestCycle.java:1418) at org.apache.wicket.RequestCycle.request(RequestCycle.java:532) at org.apache.wicket.protocol.http.WicketFilter.doGet(WicketFilter.java:356) at org.apache.wicket.protocol.http.WicketServlet.doGet(WicketServlet.java:124) at javax.servlet.http.HttpServlet.service(HttpServlet.java:707) at javax.servlet.http.HttpServlet.service(HttpServlet.java:820) at org.springframework.web.servlet.mvc.ServletWrappingController.handleRequestInternal(ServletWrappingController.java:158) at org.springframework.web.servlet.mvc.AbstractController.handleRequest(AbstractController.java:153) at org.springframework.web.servlet.mvc.SimpleControllerHandlerAdapter.handle(SimpleControllerHandlerAdapter.java:48) at org.springframework.web.servlet.DispatcherServlet.doDispatch(DispatcherServlet.java:875) at org.springframework.web.servlet.DispatcherServlet.doService(DispatcherServlet.java:809) at org.springframework.web.servlet.FrameworkServlet.processRequest(FrameworkServlet.java:571) at org.springframework.web.servlet.FrameworkServlet.doGet(FrameworkServlet.java:501) at javax.servlet.http.HttpServlet.service(HttpServlet.java:707) at javax.servlet.http.HttpServlet.service(HttpServlet.java:820) at org.mortbay.jetty.servlet.ServletHolder.handle(ServletHolder.java:502) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1124) at org.vfny.geoserver.filters.SetCharacterEncodingFilter.doFilter(SetCharacterEncodingFilter.java:108) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1115) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:265) at org.acegisecurity.intercept.web.FilterSecurityInterceptor.invoke(FilterSecurityInterceptor.java:107) at org.acegisecurity.intercept.web.FilterSecurityInterceptor.doFilter(FilterSecurityInterceptor.java:72) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.ui.ExceptionTranslationFilter.doFilter(ExceptionTranslationFilter.java:124) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.providers.anonymous.AnonymousProcessingFilter.doFilter(AnonymousProcessingFilter.java:125) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.ui.rememberme.RememberMeProcessingFilter.doFilter(RememberMeProcessingFilter.java:142) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.wrapper.SecurityContextHolderAwareRequestFilter.doFilter(SecurityContextHolderAwareRequestFilter.java:81) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.ui.AbstractProcessingFilter.doFilter(AbstractProcessingFilter.java:271) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.ui.logout.LogoutFilter.doFilter(LogoutFilter.java:110) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.context.HttpSessionContextIntegrationFilter.doFilter(HttpSessionContextIntegrationFilter.java:249) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.util.FilterChainProxy.doFilter(FilterChainProxy.java:149) at org.acegisecurity.util.FilterToBeanProxy.doFilter(FilterToBeanProxy.java:98) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1115) at org.geoserver.filters.LoggingFilter.doFilter(LoggingFilter.java:73) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1115) at org.geoserver.filters.ReverseProxyFilter.doFilter(ReverseProxyFilter.java:183) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1115) at org.geoserver.filters.GZIPFilter.doFilter(GZIPFilter.java:41) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1115) at org.mortbay.jetty.servlet.ServletHandler.handle(ServletHandler.java:361) at org.mortbay.jetty.security.SecurityHandler.handle(SecurityHandler.java:216) at org.mortbay.jetty.servlet.SessionHandler.handle(SessionHandler.java:181) at org.mortbay.jetty.handler.ContextHandler.handle(ContextHandler.java:766) at org.mortbay.jetty.webapp.WebAppContext.handle(WebAppContext.java:417) at org.mortbay.jetty.handler.ContextHandlerCollection.handle(ContextHandlerCollection.java:230) at org.mortbay.jetty.handler.HandlerCollection.handle(HandlerCollection.java:114) at org.mortbay.jetty.handler.HandlerWrapper.handle(HandlerWrapper.java:152) at org.mortbay.jetty.Server.handle(Server.java:324) at org.mortbay.jetty.HttpConnection.handleRequest(HttpConnection.java:534) at org.mortbay.jetty.HttpConnection$RequestHandler.headerComplete(HttpConnection.java:864) at org.mortbay.jetty.HttpParser.parseNext(HttpParser.java:533) at org.mortbay.jetty.HttpParser.parseAvailable(HttpParser.java:207) at org.mortbay.jetty.HttpConnection.handle(HttpConnection.java:403) at org.mortbay.io.nio.SelectChannelEndPoint.run(SelectChannelEndPoint.java:409) at org.mortbay.thread.QueuedThreadPool$PoolThread.run(QueuedThreadPool.java:522) Caused by: java.util.MissingResourceException: Unable to find resource: featureTypeReloadWarning for component: resource:tabs:panel:theList:1:content [class=org.geoserver.web.data.resource.FeatureResourceConfigurationPanel] at org.apache.wicket.Localizer.getString(Localizer.java:266) at org.apache.wicket.model.StringResourceModel.getString(StringResourceModel.java:491) at org.apache.wicket.model.StringResourceModel.load(StringResourceModel.java:633) at org.apache.wicket.model.LoadableDetachableModel.getObject(LoadableDetachableModel.java:114) at org.apache.wicket.model.StringResourceModel$AssignmentWrapper.getObject(StringResourceModel.java:230) at org.apache.wicket.Component.getModelObject(Component.java:1565) at org.apache.wicket.Component.getModelObjectAsString(Component.java:1587) at org.apache.wicket.markup.html.basic.Label.onComponentTagBody(Label.java:111) at org.apache.wicket.Component.renderComponent(Component.java:2524) ... 73 more


"Thresholding" for Stream Raster

Using ASTER DEM, I have created Flow Accumulation. However, to create stream order, I was asked to input flow accumultation and stream raster. But I cannot find how to get the stream raster. Some forums say that I should threshold flow accumulation. How do I do this? Please guide me in the right direction. I am using ArcGIS 10.


GIS Buzzwords Topic

What are the main buzzwords today in the GIS market? Let's place the meaning of these ones too, so beginners can enjoy the topic as well. Example: **SDI** "A spatial data infrastructure (SDI) is a framework of spatial data, metadata, users and tools that are interactively connected in order to use spatial data in an efficient and flexible way. Another definition is the technology, policies, standards, human resources, and related activities necessary to acquire, process, distribute, use, maintain, and preserve spatial data." wikipedia definition One answer per post :)


Python objects in my ArcMap addIn aren't instantiated

I've got a weird problem. My addin throws problems intermittently where it doesn't recognize buttons or tools when I try to toggle the visibility in my code. It seems as though the objects aren't instantiated as trying to toggle them works after I try to interact with them. What is the problem and how do I begin to debug it? My code is far too long to post but if there's something I can provide to help can you let me know? Here's a scenario: User starts Arc and clicks a button. The button tries to enable another button with:               anotherButton.enabled = False      Arc throws an error:               global variable "anotherButton" not defined.      User clicks on "anotherButton" and then tries the first button again. This time the command executes and "anotherButton" is disabled.


How to use SAGA ShapesToGrid via Processing in a Python plugin?

I run these lines in my code:               vectorlayer_vector =  ftools_utils.getVectorLayerByName(inputLayer)     Elevation = self.lineAttrib.currentText()     cellSize = int(self.linePix.value())     outPath = self.inputLayerCombo3.text()     output = ftools_utils.getRasterLayerByName(outPath)           Processing.runAlgorithm("saga:shapestogrid", vectorlayer_vector, Elevation, 0, 0, 4, None, cellSize, output)      When I run, I get an tiff image but with nothing. The .aux file is not created and when I try to open the tiff image that I created appears this message: c:...teste3.tif is not a supported raster data source So, I guess that the output is wrong. What I have to do to get a correct output? What I missing? * * * I update my code but I have no results yet. Please see if you can help me. My code are:               inputLayer = self.inputLayerCombo.currentText()     Elevation = str(self.lineAttrib.currentText())     cellSize = int(self.linePix.value())     outPath = self.inputLayerCombo3.text()          Processing.initialize()     Processing.runAlgorithm("saga:shapestogrid", inputLayer, Elevation, 0, 0, 4, "10000.0,14000.0,50000.0,54000.0" , cellSize,outPath)      Is a simple code but I don't have any results. I don't know what else to do. It doesn't create a raster.


ArcGIS Desktop can't return concurrent license

We have one concurrent use ArcView license, managed by a license server box in the office. The server reports that one of our desktops has the license 'borrowed', but when looking at the Desktop Admin tool on that box, it doesn't show anything about having a license checked out. It seems the license is somewhere lost in the ether. Is there a way, from the desktop side, to 'force return' or destroy the local license file? There appears to be no way for me to get the license returned from the desktop side. I've filed a report with Esri, but not holding out hope that my issue will be fixed ASAP. Anyone here have experience with this problem? Would a reinstall of the desktop software help?


Openlayers WMS getFeatureInfo get features with geometry attribute

I am displaying map using WMSlayer and implemented WMSGetFeatureInfo control on wms layer. Tap on a location and get the features using WMS getFeatureInfo service. This service is returning me all the features are located at given Point successfully with feature attributes. My Requirement Along with the Feature attributes, i need to have the feature geometry information too, which not geting as part of feature information from service. So is there any way to capture the geometry information from OpenLayers GetFeatureInfo control? Or Is that a limitation in WMS GetFeatureInfo service?


Setting background RGB value to invisible for multiple rasters in ArcGIS Desktop?

I'm using ArcGIS Desktop 10.0. I need to set my background rgb value of 255,255,255 to transparent. In the layer-Properties Windows > Symbology I can do it for one layer, but can I batch set this, or use arcpy? I have to do it for around 300 layers.


How can I create Multiple joins over 5 tables in ArcGIS?

I have a relational database which uses the following relational schema to create a query: Parcel.IsBooked is joined to Bookinglocation.ID Bookinglocation.IsPartOf is joined to Bookingpage.id Bookingpage.id is joined to Namenumber.isPartOf Namenumber.called is joined to Person.id This is the relationship path from the shapefile "parcel" to the table "Person" which contains the name of the parcel owner. My question is....How is it possible with ArcGIS to select parcels from the shapefile and see the selection in the "Person" table? I seem to only be able to join a shp to a table...but not a table to a table. I need to actually remove names out of the "person" table which fulfill a certain criteria...so exporting tables joined to shapefiles won´t do the job. Is this possible with ArcGIS?


How would I create an input mask for ArcGIS Desktop?

Someone asked me if we can create a simple input mask for geodatabase feature classes so that editing of attributes becomes easier. I am however not sure how this works with ArcGIS 10.x Is this possible with VBA or do I need 3rd Party tools to create masks? The last time I did something similar I could only create tools and toolbars...but no window controls...How could this be done now?


Bounding the Google Maps Search area to my Fusion Layer area?

I have a searchable Google Fusion Layer of Glasgow Postcodes that highlights the searched areas on a map. The search works fine if I enter a full postcode (for e.g. G1 1BL) that matches the data in the Fusion Table. However, if I search for just the first part of the postcode to display all the postcodes in the sector ('G1'), instead of highlighting all the areas on my fusion layer with G1, it relocates my map to China where there is a highway named 'G1'. This isn't the case for all sectors, for example if I search G84 then it highlights the corresponding areas on my map. This happens with other searches too where unless its a valid postcode I am relocated to the other side of the world to another place with this name as opposed to the fusion layer I am querying. Is there a way to fix this or will I have to limit the search to only full postcodes? Thanks in advance!


WFS-T demo with local Geoserver not responding

i have installed the geoserver in my local system and i tried to run the example code from openlayers the code as below               <!DOCTYPE html>                   <meta http-equiv="Content-Type" content="text/html; charset=utf-8">         <meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=0">         <meta name="apple-mobile-web-app-capable" content="yes">         <title>OpenLayers WFS Protocol with Filter</title>         <link rel="stylesheet" href="http://localhost:8080/openlayers/theme/default/style.css" type="text/css">

                     When i click on the save the request to geoserver is in pending for more then 10 min ![enter image description here](http://i.stack.imgur.com/3CO5m.png)


Arcpy scripts fail on big datasets - anyone has experienced similar behavior?

Here a question that is not extremely specific but maybe someone has had similar experiences... I have a few Arcpy scripts that work perfectly when used with small datasets but that silently fail when processing big datasets. E.g. I am using Network Analyst tools to determine, all over the Swiss road network, the accessibility of various services (hospital, mall, bakery, school...). I basically create an OD matrix containing the distance between each populated place and the closest service. Apparently nothing rocket- science. When my script is runnung against a small subset, it performs well and spits out a result. When run against the full road network it (usually) doesn't explicitely fail but all the output tables are just empty. I experience a similar behavior when trying to perform a kind of custom filter on a point dataset: when using a few hundreds of points, no problem but when processing the full 4 Million points, the script runs without any warning but outputs an empty table. I am aware that the question is quite vague but it's on purpose! As I seem to experience similar problems in unrelated situations, I would assume that there is a more global issue to address before trying to work on those specific scripts... So, has anyone experienced similar behavior? FYI: ArcGIS 10.1, Python 2.7, Windows Server 2008, 24GB RAM, Intel Xeon 2.67Ghz (2 processors)


How to record a base layer change event in OpenLayers?

The `OpenLayers.Control` class possesses an `eventListener` property which could be used in the layer switcher control. I am not sure how to make use of it to detect a base layer change event so that I can store the variable name of the chosen layer to be used subsequently when the user navigate back to the map page. Can anyone give me some hints?


Find the nearest point on a line

I have a line that exists between two latlon positions, and a point on a certain latlon position. What I want to know is what the nearest point on that line is (in respect to the other point). So a latlon position on that line. I know how to do the basic math, but I just can't get my head around the latlong calculation. Can anyone help me or give me some pointers on how to do it? Thanks


How to create a SHP file in QGIS using the spreadsheet?

I have a spreadsheet which have columns (X, Y, Att1, Att2 ...). I want to create a point shp file, all column in the spreadsheet to be the column in shp file. How can I do this in QGIS? Thanks.


On using R data with ArcPy to make histogram in ArcMap

I am writing a tool in python which passes parameters to a R script, which in turn creates a DBF file. I read this DBF and want to make a histogram out of the data inside ArcMap as one of the results of the tool. The code snippet is as follows:               import arcpy          simloglkhd = "H:\\CIV\\KULLDORFF\\DATA\\PLOTDATA.dbf"     inputTemplate = "C:\\Program Files (x86)\\ArcGIS\\Desktop10.1\\GraphTemplates\\Default.SimplifiedChinese.tee"     output_graph_name = "VerticalBarGraph1"     verticalBarGraph_bmp = "H:\\CIV\\KULLDORFF\\DATA\\plotdata1.bmp"     # Create the graph     graph = arcpy.Graph()          # Add a vertical bar series to the graph     graph.addSeriesHistogram(simloglkhd, "simloglkhd",10)          # Specify the title of the left axis     graph.graphAxis[0] = "Frequency"          # Specify the title of the bottom axis     graph.graphAxis[2] = "log(lambda)"          # Specify the title of the Graph     graph.graphPropsGeneral.title = "Monte Carlo Distribution of Lambda"          # Output a graph, which is created in-memory     arcpy.MakeGraph_management(inputTemplate, graph, output_graph_name)          # Save the graph as an image     arcpy.SaveGraph_management(output_graph_name, verticalBarGraph_bmp, "MAINTAIN_ASPECT_RATIO", "600", "375")      I am having the following difficulties:   1. I am not able to find a template for this histogram. In the GraphTemplates folder there is PopulationPyramid, but it is still not working. Is there somewhere I can find templates. The create template operation gives a text file instead of a .tee file.   2. I get an error in the line where I addSeriesHistogram saying "Field '10' is not in the dbf"(something like this) while it is only the countBin. Please can someone tell me where I can find good examples of the usage of this function. My dbf file is simple ID(1-999) & simloglkhd(floating point numbers). I really do not know how use this MakeGraph and addSeries_ . I have already read the ArcGIS help. Can you guide me to some more examples and documentation?


Create polygon from page size

Is there an easy way to create a polygon that matches the data frame in my page layout? I'm trying to set up data driven pages and want to have it key off of polygons that frame the boundaries for each page. I would like to see them as polygons so that I can re-position as necessary and quickly visualize what will appear on each page.


How to add LiDAR Intensity attribute to a point feature class?

I need to convert LASD (ESRI LAS Dataset) to ASCII text file. I've export the data to point feature class (LAS to Multipoint), converted multipart to single part and calculate X Y and Z.My first question is - How do I add the LiDAR Intensity attribute to the point feature class ? My Second question Is There a better way to convert LASD to ASCII ? i tried the tool las2txt at LasTools but I still looking for an answer for my first question


QGIS GPS position after entering attributes

I would like to be able to fill a form and, before, during or after, be able to cath a GPS point. I already know the GPS extension but I have to catch the point and enter attributes after... Idea how I can do that? Martin F


UK postcodes to census wards

I'm planning a project using counts of data aggregated by UK postcodes. I'd like to relate the counts to local small area census population data. Does anyone know is there lookup data available to relate UK postcodes to census wards ?


Confidence interval from raster data

I have raster data that represents a probability distribution, i.e. each cell has a probability value (in my case the probability that an animal may be found in the cell), and all the cells add up to 100% (I know for certain the animal is within the extents of my raster). I want to be able to generate vector data for confidence values. For example, the 95% line/polygon denotes the boundary wherein I am 95% confident that I will find the animal. Similarly, if I have a kernel density estimate, How do I generate the XX% line/polygon that borders the densest part of the raster containing XX% of the total population? I am willing to use ArcGIS, or open-source software. If there isn't a tool to perform this for me, what is an algorithm I can implement?


Exporting PostgreSQL table to FileGDB: wkbUnknown layer geometry type

I'm attempting to export a PostGIS (1.5.3) database table as a File geodatabase using the new FileGDB OGR / GDAL 1.9.0 driver:               $ ogr2ogr -overwrite -sql "SELECT * from db_pts LIMIT 1" \         -f "FileGDB" /tmp/test.gdb \         PG:"host='HOST' dbname='DB'"      I get this error:               ERROR 1: Error (-1): FGDB layers cannot be created with a wkbUnknown layer geometry type.     ERROR 1: Terminating translation prematurely after failed translation from sql statement.      The geometry column in the table I'm exported is set up properly and is in the `geometry_columns` table. I think the problem has to do with the way the `-sql` option is returning the layer geometry type. For example, accessing the same database table with `ogrinfo`:               $ ogrinfo -ro -sql "SELECT * from db_allobs LIMIT 1" PG:"host=HOST dbname='DB'"     Layer name: sql_statement     Geometry: Unknown (any)     Feature Count: 1     ...     OGRFeature(sql_statement):0     ...       POINT (-151.472835 69.635245)      It's as though the query doesn't look in `geometry_columns` when building the layer, but later, when exporting the rows it knows that the geometry being exported is a `POINT` (in this case). When I use a similar command to produce a shapefile, all is well, including the layer geometry, which is set to point. Anyone have any ideas?


Accessing OpenLayers on "Manage and Install Plugins"

I am a new QGIS user and I am trying to add the OpenLayers plugins. I went into 'Manage and Install Plugins', and made sure I had the 'All' tab selected. However, whenever I try to search for the pluggin, it never shows up as on option. In the dialog box I wrote out 'OpenLayers'. Other plugins show up, but not that one.


Why is my spatial join count 0 in ArcGIS Desktop?

In arcMap, I have a list of locations with lat/long. From a previous post, I learned about the Spatial Join tool to determine the amount of joins that I have. I tried to use the tool with my states layer and it worked perfectly. However, when I tried to do it on the county layer, it is always giving me with 0 joins. What is the cause of this? The county layer that I am using is from http://www.gadm.org/. The one that I used for the states came from arcgis online. There was one for the countries, but it would give me an error when trying to do a spatial join on the data.


Is there a way to retrieve/query country data from a shapefile programmatically using map coordinates in GeoTools?

I am using Geotools to parse and draw a shapefile. ![Map Drawing](http://i.stack.imgur.com/xemPu.png) Now; when I draw that map I receive coordinates with this bit of code for my `Flight Line` button.               void selectFeatures(MapMouseEvent ev) {                       /*              * Construct a 5x5 pixel rectangle centred on the mouse click position              */                       if(currentClick%2 == 0){             startPoint = ev.getMapPosition();             currentClick++;             }else{              endPoint = ev.getMapPosition();             currentClick++;             }             if(!startPoint.equals(null) && !endPoint.equals(null)){             System.out.println("The two positions are: Start - " + startPoint + "\nEnd - " + endPoint);             System.out.println("Starting point is in ");             }         }      Is there anyway to parse the map/shapefile data for the country that my `startPoint` and `endPoint` are in? ![Show Info Panel](http://i.stack.imgur.com/Udv0M.png) I can see by the `Feature Info` button that was built into the `GeoTools` demo there is a way to parse information from a point; but I am unsure of how to do this programmatically as the `API` is not very well documented.


Clipping very large polygon datasets in ArcGIS?

I'm attempting to clip a large polygon dataset (~5gb, thousands of features) by a much smaller polygon dataset (~40 features). Is there a best practice or most efficient route for performing this task? The standard geoprocessor clip runs indefinitely on a dataset of this size. Would some form of spatial selection and export be more efficient? **EDIT** : Some great answers below. I selected what I view to be the most thorough response, but each answer provides unique insight into the issue. Thanks!