title: Using Mosaicking Images in Erdas software

text: Do you know any dissertation which has written about Mosaic briefly?   I did the Mosaic for my Landsat images but i want to explain and write about but i do not know how and what shall i write about it ? is here any dissertation has written the mosaic briefly? Sorry for my bad English so that is why i do not know how to write so i really needs 100 word for this process ...if someone helps me ...i will appreciate your help.


title: QGIS Gazetteer Plugin connection to internal database

text: A few posts on here mention that the Astun QGIS-Gazetteer-Plugin could be adapted to reference a PostGIS database. Has anyone done this adaptation and could share it? I want to use a gazetteer stored in a PostGIS database, but don't want to have to create a webserver.


title: How do I get the centroid of line with OGR

text: Can somebody explain how to get the centroid of a feature with ogr? I have no idea what the `centroid()` function itself does.               >>> linelyr = lineshp.GetLayer()      >>> f = linelyr.GetFeature(0)      >>> g = f.geometry()     >>> g.GetPoint()     (5.85136, 49.974, 0.0) # starting point of the line          c =g.Centroid()     >>> c.GetX()     0.0


title: Export attribute table (feature class) to text file in batch?

text: What would be the most efficient way to export attribute tables from multiple feature class files stored in a file geodatabase? Export to .txt or .csv. Is there any sort of batch process? Or python script? I have hundreds of files within this geodatabase and really need to automate this process.


title: List the projection & datum of all the files in a folder

text: In QGIS/GRASS, is it possible to list the projection details of all the shp,tab,dxf files in a folder ? We are having 100 - 150 files in too many folders. We would like to create a table of Proj details of all the files in those folders. Any idea ?


title: Visualize bounding box to square like polygon

text: I want to add script in my plugin to visualize bounding box into a polygon shape, bounding box i try to visualize have format like this:               (X1,Y1) (X2,Y2) (X1,Y2) (X2,Y1)      Does QGIS/PYTHON have modul/script to do such thing?


title: MYSQL 5.6 Cluster / chain points by distance to any other point without millions of calculations?

text: I have a table of approx 10,000 - 20,000 dynamic points across the globe stored as lon/lat points. I want to cluster any point within (roughly) 3km of any other point into one point (and return the averaged center point of each cluster if possible). Unlike some other solutions I have seen, I do not just want points within 3km of the central cluster point, it is fine to have a point within each cluster many kilometres away from the central cluster point as long as there is a chain of points within 3km of each other joining the points together into a single chain/cluster. The only way I can think of doing this is looping over every point, finding any other points within 3km of that point, and then recursively each of those points etc until i find all in that perticular chain/cluster, and then repeat for all the remaining un-clustered points - this would result in millions of calculations. Is there an easier way of doing this? Preferably with MYSQL 5.6 geospatial abilities, but im open to any viable solution.


title: PostGis geometry name lookup failed for table

text: I have a table name erpAssets which has column name Region of type polygon. it contains entries of points also like (-79.4609576808001,43.9726680183837)). I have written query which convert these points to geometry which is used by mapnick to generate tiles. ( SELECT replace(concat(concat('POINT',concat(trim(TRAILING ')'FROM concat('(',trim(LEADING '(' FROM region::TEXT)::TEXT) ),')')),''),',',' ') as data from "erpAssets" where parent=326815) as foo   ) as points My osm.xml file is as follows:               <?xml version="1.0" encoding="utf-8"?>     <!DOCTYPE Map [     <!ENTITY % entities SYSTEM "inc/entities.xml.inc">     %entities;     ]>     <Map background-color="#bb0000" srs="&srs900913;" minimum-version="2.0.0">       &fontset-settings;      <Style name="theme_park">          <Rule>           <PointSymbolizer file="&symbols;/alpinehut.p.16.png"                 />              </Rule>          </Style>          &layer-shapefiles;     <Layer name="theme_park" status="on" srs="&osm2pgsql_projection;">         <StyleName>theme_park</StyleName>         <Datasource>           <Parameter name="table">             (SELECT ST_GeomFromText(data)     FROM     (     SELECT replace(concat(concat('POINT',concat(trim(TRAILING ')'FROM  concat('(',trim(LEADING '(' FROM region::TEXT)::TEXT) ),')')),''),',',' ') as data      from "erpAssets" where parent=326815) as foo                    ) as points                </Parameter>           &datasource-settings;         </Datasource>     </Layer>     </Map>      Tiles are not generated and i am getting error on console as follows: Postgis Plugin: geometry name lookup failed for table '"erpassets"'. Please manually provide the 'geometry_field' parameter or add an entry in the geometry_columns for '"erpassets"'. I have converted region column into geometry by aliasing then why i am not getting points on map?


title: Using OGR Intersection Method Python

text: I'm currently working to find the intersections between linears in a shapefile. I'm a newbie to the GDAL/OGR library. I've done some googling and came up with what seems to be a useful example(see below). There is a problem though. I don't see the method ExportToWkt(). My intellisense(or whatever it's called in IDLE) for the variable "intersection"doesn't have that option. Can anybody tell me what I'm doing wrong. Or just give me some general guidance to comparing two features in a shapefile and finding the intersection(if they intersect at all). Thanks. ![enter image description here](http://i.stack.imgur.com/xcFyf.png)               from osgeo import ogr          wkt1 = "POLYGON ((1208064.271243039 624154.6783778917, 1208064.271243039 601260.9785661874, 1231345.9998651114 601260.9785661874, 1231345.9998651114 624154.6783778917, 1208064.271243039 624154.6783778917))"     wkt2 = "POLYGON ((1199915.6662253144 633079.3410163528, 1199915.6662253144 614453.958118695, 1219317.1067437078 614453.958118695, 1219317.1067437078 633079.3410163528, 1199915.6662253144 633079.3410163528)))"          poly1 = ogr.CreateGeometryFromWkt(wkt1)     poly2 = ogr.CreateGeometryFromWkt(wkt2)          intersection = poly1.Intersection(poly2)          print intersection.ExportToWkt()


title: Transforming between epsg:3857 and Google Maps Tile coordinates

text: Hio, I used to use a simple python project (https://github.com/hrldcpr/mercator.py) to convert between wgs84 and google maps tile coordinates. This project doesn't have a license, and looks like it may be inaccurate (assuming a spherical earth?), so I decided to use pyproj. I understand that google maps uses a projection called epsg:3857 to draw their tiles. When I convert from wgs84 to epsg:3857, I get numbers that I just don't understand:               >> import pyproj     >> epsg3857 = pyproj.Proj(init='epsg:3857')     >> wgs84 = pyproj.Proj(init='EPSG:4326')     >> latlng = (45.618, -73.604)     >> pyproj.transform(wgs84,epsg3857,latlng[0],latlng[1])     (5078172.5310075525, -12357511.560866801)      I know that the lat,lng coordinate I took is in the tile `http://b.tile.openstreetmap.org/11/605/731.png`, But I just can't figure out how to go from the returned numbers to the tile coordinates. Firstly, the numbers shouldn't be negative (tile coordinates are positive). Trying different scales based on the 11 zoom levels, I just can't figure this out. Searching the internet, I can't figure out what the epsg:3857 coordinates actually denote.


title: How to rasterize SpatialPolygons in R?

text: I'm trying to extract the bathymetry values of my area of interest from a world bathymetry raster layer using the 'rasterize' function in the {sp} package. *Edits: I found the 'extract' function which seems to be more what I'm looking for. This is what I've done so far:               > class(subarea0) #This is my area of interest (Eastern Canadian Arctic Sea)     [1] "SpatialPolygons"     attr(,"package")     [1] "sp"          > extent(subarea0)     class       : Extent      xmin        : -82.21997      xmax        : -57.21667      ymin        : 60.2      ymax        : 78.16666          library(marelac)     data("Bathymetry")#World bathymetric data in library (marelac)     names(Bathymetry);class(Bathymetry);str(Bathymetry)     [1] "x" "y" "z"     [1] "list"     List of 3      $ x: num [1:359] -180 -179 -178 -177 -176 ...      $ y: num [1:180] -89.5 -88.5 -87.5 -86.5 -85.5 ...      $ z: num [1:359, 1:180] 2853 2873 2873 2873 2873 ...            raster_bath<-raster(Bathymetry)#Transformed into a raster layer         extent(raster_bath) <- extent(subarea0)#Transform the extend of my raster to the extend of my SpatialPolygons          >ras_sub0<-rasterize(subarea0,raster_bath)#rasterize my SpatialPolygons (*Edits: not the function that I need here, but I am still interested to learn what results mean)     Found 2 region(s) and 10 polygon(s)     > plot(ras_sub0)     > plot(subarea0, add=TRUE)      ![enter image description here](http://i.stack.imgur.com/nJMrB.jpg)               > ras_sub0     class       : RasterLayer      dimensions  : 180, 359, 64620  (nrow, ncol, ncell)     resolution  : 0.06964709, 0.0998148  (x, y)     extent      : -82.21997, -57.21667, 60.2, 78.16666  (xmin, xmax, ymin, ymax)     coord. ref. : +proj=longlat +ellps=WGS84      values      : in memory     min value   : 1      max value   : 2      layer name  : layer       I don't understand the result. Why am I getting 2 colors for each of my polygons? What do they mean? How can I ultimately get bathymetry depth contour? Is this something to do with my resolution or changing the dimensions? *Edits: ok, I've now done the following:               v <- extract(raster_bath, subarea0)#Extract data from my Raster Layer for the locations of my SpatialPolygons      v is a list and I am not quiet sure how/under what form to rebind this info with my spatial polygon... Thank you!


title: IFeatureLayer.Seach() throws exception HRESULT:0x80041538

text: My code as bellow give error HRESULT:0x80041538 whenever I use GROUP BY or ORDER BY postfix. Trying in Personal(.mdb) and MSSQL database datas - same result. What can be the reason. Asked similar question before, but according to what I found that time, in MSSQL it should be working.               IQueryFilter qf = new QueryFilter();     ((IQueryFilterDefinition2)qf).PostfixClause = "GROUP BY " + groupbyfields;     ILayer layer = Utils.getLayer(map, layer_name);     ICursor ccursor = null;      cursor = ((IFeatureLayer)layer).Search(qf, false);


title: Looking For Free Global Bedrock Geology Data

text: I've searched around for a freely available global bedrock geology data set, but I have yet to find one. Anyone come across such dataset?


title: Can QGIS Read Spatialite Views?

text: I have a spatialite database I created that spans several years of crime data. To facilitate a real-world scenario where I may simplify the interaction for a user by introducing a View that simply limits the time span or what types of crimes are returned, I want to know if it is possible for QGIS to read this View like it would any other spatial table in my database. The point here is that I obviously do not want to create a separate table for every possible View I may use, and I am trying to make it simple by not requiring the user to know the database and the appropriate SQL to subset the data themselves. A similar workflow in the ArcGIS world might be creating separate layer files that the user could open in ArcMap that access the database (say, SQL Server) and present only the limited data set. I created a View in my database, but it does not show up in the list of feature classes I can import when I connect to the spatialite db in QGIS. Maybe Views aren't considered spatial? (This is something I'm going to look into later, but I thought I already did; I just don't recall the conclusion.) In any case, how might one replicate the sort of work flow I want or the similar ArcGIS alternative?


title: Node.js transformation from NAD83 to WGS84

text: I've used the excellent npm shapefile library to import Canadian electoral districts into MongoDB. Here's a sample file: http://ftp2.cits.rncan.gc.ca/pub/geott/electoral/2011/fed308.2011.zip Is anyone familiar with a node.js tool that would allow me to do this? I'm familiar with Django and the GIS tools provided by Geodjango, but am new to GIS development in the Node.js ecosystem. Thanks in advance!


title: Attribute inspector not editable - FLEX or Javascript

text: I am using ArcGIS 10.2. I have an application that is using the editor and the attribute inspector. I need to add a string showing the current date to a field. Everything works but after creating the feature and selecting it again I can see the string in that field in the attribute inspector but the attribute inspector is not editable anymore for that feature. Is it possible to add a string to a field and still have that feature editable in the attribute inspector? Thanks                   var    currDateTime:String = ""         private function drawEnd(event:DrawEvent):void            {             dateFormatter.formatString = 'YYYY/MM/DD HH:NN:SS';           currDateTime = String(dateFormatter.format(dtCurrent) + '.' + dtCurrent.getMilliseconds());                    FeatureLayerstewardship.addEventListener(FeatureLayerEvent.EDI}TS_STARTING, editsStartingHandler);                    }                 private function editsStartingHandler(event:FeatureLayerEvent):void                     {                         if (event.adds)                         {                                  event.adds[0].attributes['DrawDate']  =  currDateTime;                         }                               }


title: How to set up ArcGIS Server 10.1 test environment on Mac laptop with VMware?

text: In my production environment I have ArcGIS Server 10.0 Enterprise, ArcSDE 9.3, and ArcIMS. I want to shut down ArcIMS as soon as possible for obvious reasons, and upgrade ArcSDE 9.3 to 10.1 and ArcGIS for Server 10.0 to 10.1. As I am new to ArcGIS for Server 10.1, I would like to set up a test environment just to get to know the software architecture but will not be able to do this on my network at work. I would therefore like to do this on my laptop with VMware. Could anyone guide me as to how I should do this on my Mac with VMware?


title: openlayers Arcgis WFS mousehover

text: i have a problem using ArcGIS wfs layer in openlayers. I would like to add the layer to the map as vector layer and then using a mousehover to highlight the features. However i cannot bring it to work. Code example:               var buildings_wfs = new OpenLayers.Layer.Vector("Buildings_WFS",                        {strategies: [new OpenLayers.Strategy.Fixed()],                         visibility: true,                         protocol: new OpenLayers.Protocol.WFS.v1_1_0({                         srsName: "EPSG:3857",                         version: "1.1.0",                         featureType: "Buildings",                         url: "http://XXX/MapServer/WFSServer",                         featureNS: "http://XXX/MapServer/WFSServer",                         geometryName: "Shape",})});      any help is much apreciated.


title: How can I move two layers at once in QGIS?

text: I'm trying to move two digitized layers, one containing points, the other polygons to a different location at the same time. QGIS doens't seem to allow this. Even if I select both layers and features it will only move one (the one I last selected, so which ever feature is shown in editor) I know ArcGIS allows multiple layers to be moved at once, just with the moving tool or with the spatial adjusment. I haven't been able to reproduce this in QGIS. Thanks for any help!


title: Linking sld to vector layer

text: I have created a vector layer by reading a geojson file and I have a SLD file created for this layer. I am wondering how do I link the SLD file to the layer. I had a read of this example on Openlayers website: http://openlayers.org/dev/examples/sld.html I am having a bit of trouble understanding the code for this example- which are using radiobuttons to change the style to the selected button. I believe the linking in my situation will be much simple than what is showing in the example- just linking one SLD file to one single vector layer. Any guidance will be greatly appreciated.


title: Coxcomb chart in qgis

text: I would like to use special pie chart for my map, where the angle is not important, but area or radius of the arc carries information. I found, that this type of pie chart is calling Coxcomb/polar/rose (in Czech it is known as comparative) chart (see http://blogs.sas.com/content/jmp/files/2012/03/spending1.png or http://visual.ly/learn/coxcomb-chart), basicaly it is bar chart masked in pie :). Do you have any suggestions how to create this chart automaticaly or semiautomaticaly in gis (I prefer QGIS)? Thanks for help, Daniel


title: Looking for free/paid RSS to GeoRSS converters

text: I am looking for an RSS to GeoRSS converters . I was able to find only one converter i.e. GeoNames (its crashing some-times) and Meta-carta (looks like closed) So does anyone knows free/paid RSS to GeoRSS converters ? Software's : ArcGIS Suit 10.1 & ArcGIS JSAPI 3.3/3.4 Thanks in advance : )


title: ArcGIS Geoprocessing Script - Calling External EXE Works in Desktop - Causes GP Service to Quit Abruptly

text: I have an interesting problem on my hands... I have a generic python script that does a bunch of geoprocessing functions using the standard arcgisscripting module (Yes, I am still at 9.3.1). At one point in this script, I call out to an external exe. This external exe does some heavy duty raster processing using gdal17.dll. After the .exe does its thing, my script continues on doing additional processing using the results from the .exe - still using arcgisscripting functions. The tool runs beautifully in a desktop environment. However, as soon as I publish it as an ArcGIS Server Geoprocessing Service - it exits abruptly in the first calling of arcgisscripting function **AFTER** the gdal-based exe is called. I can't log or trap an error - it just plain 'ol quits. I am very aware that mixing gdal in the arcgis geoprocessing environment is looking for trouble. However, porting the logic from the exe is just not an option at this point - if it is even possible. The fact the script works fine in the desktop environment gives me a glimmer of hope that it may work in a gp service. Any suggestions greatly appreciated.


title: How do you choose the classifications on choropleth maps?

text: I find myself more and more struggling over defining break points when displaying choropleth (aka thematic) maps to view by _others_. Does anyone have any suggested references that help guide, both how to choose the type of scale used _and_ the appropriate number of break points? In particular for the number of bins I have only ever seen arguments for a limiting number (e.g. you should not use any more than 5). * * * To be more specific about what I am looking for, most references I have come across about the subject are similar to the document referenced by julien in this post, and I'm just looking for a more in-depth discussion about the topic. A few specific use cases I come across frequently (for examples of my struggles);   * When displaying data that has a large right skew, I'm typically hesitant to display an exponential scale. I fear (for the audiences I am typically displaying maps to) this would cause a greater amount of cognitive burden reading the scale and mapping actual attribute values to the colors. Are my fears incorrect? Also for these types of distributions I find it difficult to justify any particular number of bins.   * When displaying many small multiple maps, how do I choose an appropriate scale that allows one to visualize relationships effectively both within and between the small multiples? My de-facto standard when the attribute scales vary to a great extent is to use quintiles within each seperate distribution. Are quintiles too many classifications and creating too great a cognitive burden to compare between the panels? I assume people understand quantile classifications are equivalent to rankings (and thus when classed that way aids in interpreting between panels), is this assumption correct? I initially wrote a paragraph trying to describe the goals of such maps, but I suspect my goals are pretty typical so it was unneeded. The only thing to clarify again is that these are for viewing by other people (like in reports, publications) and are not really for my own exploratory data analysis (although I would suspect good advice should translate to either). Perhaps a good reference may describe the potential goals of such maps, and trade-offs associated with using different classification schemes. I would be interested in both specific and general references.


title: Where can I download all 3,141 County US Polygons clipped with high resolution shoreline data?

text: After having done quite a bit of research I have been unable to find a single shapefile download which includes all of the current 3,141 counties with their associating data clipped with the highest resolution shoreline data available. Does anyone know where I can find this data?


title: How to delete a table field in ArcGIS?

text: I am trying to delete a table field: From ArcMap table view, I right click on the field heading and select `Delete field`. Unfortunately, I obtain the following message:               Failed to delete the selected fields.     General function failure [TABLE_NAME]      and the field remains. How should I proceed to delete this damned field?


title: How to query distict row with the same label from different linestring with multiple segment?

text: I have a multilinestring street table in postgis with two columns: label and prefix. This table has 2 different roads with the same name and each road contains multiple segments. Roads were splitted on each intersection during authoring in QGIS so that the table has total 14 rows. How to query the two so that i have two unique street rows with prefix and label column in the result?


title: Splitting raster pixels by a polygon's edge?

text: I have a a polygon shapefile of different regions, and a raster grid of precip data. I want to take the polygon and add up all the raster precip data in each zone. I could use zonal statistics to do this, however the tool does not split pixels exactly where the polygon edge cuts through the a pixel, it instead takes either the whole pixel value or not at all. I could convert the raster data to polygons, but I was wondering if there was a way to do this without going to polygons.


title: DEM sources for arctic regions

text: Are there any 90m DEMs out there for regions north of where the SRTM dataset ends? I looked at ASTER, but it seems to have a lot of unfilled voids, and it's important that the resolution matches the datasets I already have.


title: JTS difference op followed by an intersection op not returning expected results

text: I have two overlapping polygons and I want to get the line segment(s) where the edges cross. ![enter image description here](http://i.stack.imgur.com/z6gE3.png) So I begin by performing a Difference operation. That gives me the area where A is "outside" of B. ![enter image description here](http://i.stack.imgur.com/Wtl0k.png) Then I do an Intersection op. I expect that to give me the line segment where the two polygons "touch", but as you can see I get a smaller than expected result that doesn't extend along the entire edge where they touch. You can see in the above shot where the green result is, but here it is isolated from the two polygons. ![enter image description here](http://i.stack.imgur.com/aTNLw.png) What am I doing wrong? Is there a better way to go about this?


title: Issue selecting feature for IdentifyTask (Javascript API)

text:
I am having an issue selecting features in my web map for an IdentifyTask operation. I have a dynamic map service that contains a few sublayers that I would like to be able to click and have a popup or infowindow display the attribution for the selected feature. Currently onclick of the "map", the IdentifyTask execution fires off. The problem is that the precision of the map click seems to be off. When I click a feature with the cursor, the feature is not selected and the InfoWindow displays 'No Results' way below where I clicked the cursor. The InfoWindow typically is supposed to display where a user clicks the cursor right? I have to click above the intended feature in order to actually select the feature and have the IdentifyTask recognize it. Its weird. I played around with the tolerance values for the IdentifyParameters object but no luck. Has anybody ever had a similar issue? Is there a way to specify the precision for a mapClick event? Any help would be greatly appreciated. Below is the code I'm using to show the results of the IdentifyTask.               //Create IdentifyTask object     userConfig.identifyTask = new esri.tasks.IdentifyTask(currentServiceUrl);          //Create and setup IdentifyParameters     userConfig.identifyParams = new esri.tasks.IdentifyParameters();     userConfig.identifyParams.height = userConfig.map.height;     userConfig.identifyParams.width = userConfig.map.width;     userConfig.identifyParams.mapExtent = userConfig.map.extent;     userConfig.identifyParams.tolerance = 3;     userConfig.identifyParams.returnGeometry = true;          //Define the layers to identify     var layersToId = [];     if (currentUtilId == "barksdale") {         var layers = userConfig.barksdaleService.visibleLayers;         for (var i = 0; i < layers.length; i++) {             if (layers[i] != 27 && layers[i] != 40 && layers[i] != 82) {                 layersToId.push(layers[i]);             }         } //End loop     }     else if (currentUtilId == "malmstrom") {         var layers = userConfig.malmstromService.visibleLayers;         for (var i = 0; i < layers.length; i++) {             if (layers[i] != 27 && layers[i] != 40 && layers[i] != 82) {                 layersToId.push(layers[i]);             }         } //End loop     }     else if (currentUtilId == "minot") {         var layers = userConfig.minotService.visibleLayers;         for (var i = 0; i < layers.length; i++) {             if (layers[i] != 27 && layers[i] != 40 && layers[i] != 82) {                 layersToId.push(layers[i]);             }         } //End loop     }     else if (currentUtilId == "whiteman") {         var layers = userConfig.whitemanService.visibleLayers;         for (var i = 0; i < layers.length; i++) {             if (layers[i] != 27 && layers[i] != 40 && layers[i] != 82) {                 layersToId.push(layers[i]);             }         } //End loop     }     else if (currentUtilId == "warren") {         var layers = userConfig.warrenService.visibleLayers;         for (var i = 0; i < layers.length; i++) {             if (layers[i] != 27 && layers[i] != 40 && layers[i] != 82) {                 layersToId.push(layers[i]);             }         } //End loop     }     userConfig.identifyParams.layerIds = layersToId;     userConfig.identifyParams.geometry = evt.mapPoint;          var deferred = userConfig.identifyTask.execute(userConfig.identifyParams);          deferred.addCallback(function (response) {         // Response is array of identify result objects         // Return an array of features.         return dojo.map(response, function (result) {             var feature = result.feature;             var attribution = feature.attributes;             var infoTemplate = new esri.InfoTemplate();             var templateString = "";             //Parse and print the key-value pairs in the attribution object             for (prop in attribution) {                 templateString = templateString + "" + prop + ": " + attribution[prop] + "
";             }             infoTemplate.setTitle("Identify Results");             infoTemplate.setContent(templateString);             feature.setInfoTemplate(infoTemplate);             return feature;         });     });          userConfig.map.infoWindow.setFeatures([deferred]);     userConfig.map.infoWindow.show(evt.mapPoint);


title: Create Layer editing/Analysis in an ArcGIS web map application using Javascript API

text: I want to make an ArcGIS web map application in which user can select different layers of the published map (data of the layer in Postrgresql) and do some analysis based on layer data.import data to the application. Some kind of Interactive map. 1) Do I need Geodatabase or just registered postgresql database with Arcgis for Server and ArcGIs for Desktop 10.1 is enough? 2) any sample template for javascript API is available? Thank you


title: Tiling/Caching ArcGIS Server

text: A few weeks ago we crashed our server by trying to tile a large image service (600gbs) (multiple times without clearing the cache) which caused us to exceed the servers free disk space. We fixed the issue by deleting the service and clearing the server cache. From this point forward we have not been able to publish and tile an image service, however we are able to publish a regular map service. Does anyone know what might be causing us to not be able to publish/tile image services? Is there a potential break in code somewhere in a config file on the server? We are running Arc 10.1 and we have substantial free disk space. Thank you.


title: Layers not lining up

text: I have three layers I need to try match up at the borders. This is not a projection problem, all layers are projected to Albers NAD83. The base layer (Green) is 30m DEM, with a US boundaries layer on top of the DEM and finally a polygon "layer of interest" on top (Red). As you can see in the figure below, none of these layers match up at the borders. What is the best way to get these layers to match up in ArcMap 10.1? This is mainly for aesthetics, as I have it the red layer is transparent and it just looks bad. EDIT: So figure 2 (bottom) is what I have done thus far. (Note: I have changed the colors around in figure 1 just to make the different layers easier to distinguish for you all.) As you can see, the problem with the layers lining up isn't all that bad at this scale.. but I'm picky like that. Figure 1 is zoomed in to an area on the coast of South Carolina (apologies for no scale bar). @Kevin, the state boundaries were requested because this project spans many states and the hillshade (while not required) was also requested and adds some pizaz. Thanks everyone for comments so far. ![enter image description here](http://i.stack.imgur.com/KjwsJ.png) ![enter image description here](http://i.stack.imgur.com/exDbT.jpg)


title: Map from Mapbox as layer in OpenLayers

text: I want to use this map from mapbox as base layer in OpenLayers: http://mapbox.com/#/tileset/iraq I used this example to make it work. Everything is fine when I keep the code as it is with original layer "world-light", but when I change it to "iraq", it doesn't work. I checked it in Firebug and it looks fine - it asks for "iraq" tiles as the original page. Do you have any ideas how to solve it? Thanks


title: Standardised Spatial Information System Vocabulary

text: I've had some musing about this for some time and thought I'd ask the forum. _`Are there any resources out there that provide a standard, comprehensive dictionary of acronyms, terms and keywords for the Spatial Industry?`_ I've seen the Spatial Industry make leaps and bounds into other sectors outside of it's geology/geography roots, and with this terms such as GIS (geographic information system), Coordinate System, I believe are beginning to lose their potency in definition and are more historical terms due to the introduction of spatial technology to an increasing audience of users, developers and enthusiasts. The most common example I come across, is the term Coordinate System. Typically used to defined a projected and geographic coordinate system to the uninitiated, it can become confusing or add ambiguity to an answer if the merits or differences between coordinate systems is being discussed. As a result I prefer to use the term Spatial Reference System and offer the sub categories of projected reference system and geographic reference system. Essentially, I'd prefer to align myself with a set of agreed terms that can be addressed. If such a reference exists. Cheers, OptimizePrime


title: How do I render data quicker in Arcmap?

text: Is there a way around to speed up the data into ArcMap ? When I converted the raster to polygon and some data would load up quicker and some other load slow. I am wondering if I could add the spatial index to the data that load slow. Would that work ? or Is there anything that I could make the data load quicker ? Thanks


title: How display fields from geodata service c#

text: I want to display my geodata stored in the arcgis server. I did the connection to the server but the problem is the parameter versionName. I will not use arcSDE so I have ignored versionName. Can someone help me to find the solution? thx :) 'ESRI.ArcGIS.GISClient.IAGSServerObject so = GetSelectedServerProxy(); ESRI.ArcGIS.GeoDatabaseDistributed.IGeoDataServer gdso = (ESRI.ArcGIS.GeoDatabaseDistributed.IGeoDataServer)so; //++ Queries for subfields if the Get Fields checkbox is enabled ESRI.ArcGIS.Geodatabase.IQueryFilter qf = new ESRI.ArcGIS.Geodatabase.QueryFilterClass(); qf.WhereClause = ""; //++ based on the start index and count entered. ESRI.ArcGIS.Geodatabase.IResultPortionInfo range = new ESRI.ArcGIS.Geodatabase.ResultPortionInfoClass(); range.StartIndex = 0; range.Count = 500; //++ Run a tablesearch based on the version, table name, fields selected, startindex and count ESRI.ArcGIS.GeoDatabaseDistributed.IGDSQueryResultPortion res = gdso.TableSearch(versionName,"tabl", qf, range);'


title: Examples of geospatial scientific posters

text: Can you point me to examples of interesting scientific posters featuring results of geospatial analysis and/or geovisualizations?


title: Where can I get a world map with boundaries of regions and communities for commercial use?

text: I am looking for a vector file worldmap which should include as many boundaries reaching to a very small level (like ZIP-Codes/small communities/etc) as possible. The data of Natural Earth is quite good, but does not go that much in detail. The GADM data would be nearly perfect, but it is not allowed for company use. Requirements:   * Consistent boundaries for nearby countries (no overlapping!).   * As much detail as possible (refering to boundaries).   * Country boundaries and subdivisions should be separated.   * The data is for **company use** , therefore I should be able to buy it, or it should be free to use for companies. I also thought of extracting it from the OpenStreetMap data. Do you think this is possible? If this forum is not the right place for this question, please give me a hint for a better place.


title: make: pyuic4: command not found

text: I need to develop a QGIS plugin for a project. First, I launched the following command for the installation on Ubuntu 12.04:               sudo apt-get update     sudo apt-get install qgis     sudo apt-get install python-qgis      Then, I installed the **Plugin Builder** plugin. Finally, I tried to follow this simple tutorial. However, I'm stuck at the execution of the **make** command, which give to me the following error: **make: pyuic4: command not found** I did some research, but I'm unable to find a solution. Can someone help me?


title: Is there a way to create a plugin in C++ without compiling all the Qgis source?

text: This is my second post and I apologize if it's post the same day on the first. After several researchs, I found the only way to create a plugin in C++ is to compile the Qgis source with the plugin created thanks to the "plugin_builder.py" which can be found in the "src/plugins" on the git qgis (https://github.com/qgis/Quantum-GIS/tree/master/src/plugins) To create a simple plugin which can be add simply by putting the folder of the plugin in ".qgis/python/plugins" or "C:/OSGeo4W/apps/qgis/python/plugins" the plugin must be written in Python, It is rigth ? So my question is simple : Is there a way to create a plugin in C++ without compiling all the Qgis source ? Thank you very much Regards


title: Calculating a projected distance

text: I am pretty new to this and I am not sure I understand all the concepts, so I might be using words wrongly or obviously not getting some details right. Here goes. We mean to do those calculations within our program using the ArcObjects API. This is an issue related to how to use the ArcObjects API to achieve the desired result We want to do a proximity search. From a given point in projected coordinates, I want to look for objects (I think the word in ArcGIS is 'features') within a given distance given in meters (or some other real-world unit). Keep in mind here that the objects are _not_ contained in an ArcGIS database: we have the objects in projected coordinates in some other, external, list. All I want is to convert the distance from meters to a projected distance, given in projected units, whatever that is.


title: Allowing Add/Delete/Update of features (web application): OpenLayers + WFS and GeoServer?

text: I've developed a WPF application that has a custom SQL Server database and uses the Esri ArcGIS WPF API to allow users to view, edit and query shapes (features) on a custom map. Maps are created from tiles using the Bing (Quadkey) system. Users can also create multiple layers. This application has a great focus on the editing aspect. Altought I use the ArcGIS WPF API from Esri, there are no other services/servers from this company being used for this solution. Layers, features and symbols are stored using custom structures. But now I'm looking to make the transition to a web application. I've studied and tested differente mapping libraries: Leaflet, OpenLayers, Google Maps API, ArcGIS JavaScript API, etc. I know about WMS/WFS and mapping servers like GeoServer and MapServer (but I don't know what would be the best match). So, finally, my question: what would be the "best" (recommended) solution for a web application using custom (tile) maps that should have a great experience to allow users to add/delete/delete features, query, filter and save them? I know that the "best" solution can vary greatly. But what I'm really looking for is some advice from experience users that have done this. **OpenLayers + WFS and GeoServer** would be a good match? Thanks! Any help would be appreciated! :)


title: USGS DEM Easting and Northing metric and unit

text: Greetings to all! First post in here. Firstly, I am a complete noob in geography. But I have to do some programming work with DEM files. I got some DEM files from http://geobase.ca/geobase/en/find.do?produit=cded the 1:50,000 version and read in the data. According to the DEM specifications, there is a quote: > "... a 4x2 array containing the ground coordinates of the quadrangle > boundary for the DEM. The coordinates of the quadrangle corners are ordered > cloclwise direction beginning with the SW corner. The array is stored as > pairs of eastings and northings." Question: 1. Is easting and northing same as UTM metric? 2. Are UTM coordinates supposed to be written as easting-northings or northing-eastings? 3. What is the standard unit for UTM coord? meters? arc seconds? In the DEM file i have, I figured it is in arc seconds. 4. How to I convert, specifically for this DEM (see above) I have, from the "4x2" array to meters? Below is one of the DEM header data i extracted. Which this data, specifically, the SE coord (2nd line) and the 4x2 array (SW, NW, NE, SE coord) is it sufficient to point on the map where this quadrangle is? I tried very hard to map this data to google earth but i cannot seem to work out the maths.... hope someone can enlighten me. Thanks in advance!               DEM Name: 092g01DEMe     SE Coord: -12215 0.0000, 49 0 0.0000     DEM Level Code: 1     Ground Reference System: Geographic     Ground Reference Zone: 0     Ground Unit: arc_seconds     Elevation Unit: meters     Ground Resolution (lat, lng, elev): 0.75, 0.75, 1     Elavation Array Szie: 1201 x 1201     Percentage void: 0     SW Coord: -440100, 176400     NW Coord: -440100, 177300     NE Coord: -439200, 177300     SE Coord: -439200, 176400


title: How to create a generic tool with ModelBuilder in ArcGIS?

text: We are working on the creation of a tool with **model Builder** in ArcGIS 10.1. We want to build a generic tool that can be launched from any computer. In order to do so, we have set the input files as parameters in the model. When launched from the **Toolbox** , the paths on our computer come up in the parameter boxes. The user can change them, but we would like to have blank boxes instead. To do so, we delete all relative path of the parameters in the model. When we test the model in ArcMap, we get the following message “ **Error Warning 000873** ” which says that our workspace does not exist. Also we have a second kind of error: “ **ERROR 001000** ” which means that some value of zone field does not exist. **Is there anyone who has executed this operation with success and could explain us the right way of process?**


title: Altering composer label items in QGIS with Python

text: I'm only new to QGIS and Python so this may exists already, and hopefully be easy to do. What I'm hoping to do sometime in the future is build a map book generator using QGIS as the base. In order to do this I need to be able to select all the labels in the composer and replace some tags, so something like $NAME$ would replace with a variable the user has set at runtime. I have been having a look at the QGIS APIs [and source] but just can't seem to get it to do what I need. The code I have so far:               composerList = self.iface.activeComposers()     if(len(composerList) < 1):           return          composerView = composerList[0]     composition = composerView.composition()     if(composition is None):           return          for item in composition.selectedComposerItems():     //Need to check for QgsComposerLabel type and update the text.      I'm stuck when it comes to checking the type of the current item in Python, `selectedComposerItems()` just returns each item as a `QgsComposerItem` and calling `displayText()` fails because it doesn't exists on the type QgsComposerItem. I have tried the following code but it crashes QGIS when it hits a legend or scale bar:               for item in composition.selectedComposerItems():         try:             item.__class__ = qgis.core.QgsComposerLabel             if item.displayText() == "$NAME$":                 print "I found one"                 item.setText("Hello World")         except:             print "Fail"      The idea was to just try setting `__class__` on every item as `qgis.core.QgsComposerLabel` and if it throws an exception then I know it's not a label, but at the moment if it hits a legend or scale bar it crashes QGIS which is not good.


title: PostGIS + ArcSDE mixed environment

text: I'm trying out a mixed environment here at work with PostGIS and ArcSDE together. I'm having some problems creating feature classes with custom SRIDs. I imagined that all that I needed to do was update spatial_ref_sys table, but it did not worked also. Any tips on how to achive this with custom SRIDS? With normal srids it's working just fiiine. Thanks


title: How do you create a point at the midpoint of a polyline with Python?

text: I have a polyline shapfile (that contains multiple lines) and am trying to create a point at the midpoint of each line using python. I am currently using the ET Geowizards tool 'PolylineToPoints' which does exactly what I want it to...but I am trying to avoid using ET Geowizards (as it is causing other issues). Is there a way to do this in python using Arc tools and functions? Thanks a bunch!


title: Inaccurate result while clipping on gis data

text: I am trying clip polyline data using a polygon geometry. I am using itopologicaloperator to do this task. However, resultant clipped data is off by 1 feet. Meaning that the polylines at the boundaries of the polygon geometry goes beyond the boundary by 1 feet. Here is my case: 1) the polygon geometry is in GCS_North_American_1983 2) the polylines are in GCS_WGS_1984 I tried to generate the result by projecting the polygon geometry data to GCS_WGS_1984 ((I used NAD1983_To_WGS1984_5 for geographic transformation to equidistant conic projection) before clipping the data . The resultant clipped data is not accurate in that polylines cross the polygon by 1 feet or more. I tried to generate result using ESRI tools with accurate results. I do not understand why I could not do this with my code. Here is my code. Any help would be appreciated. Thanks. code:               ITopologicalOperator topo;     IFeatureCursor cursubdiv = null;     IClone geoclone = null;     IFeatureCursor curtrack = null;     IFeature feattrack = null;     IGeometry geotrack = null;     IGeoDataset geosubdiv = null;     ESRI.ArcGIS.esriSystem.IUnitConverter unitConverter = null;          //tracks are in gcs_1984     //Polygon geometry is in GCS_North_American_1983     //so convert track geometry to gcs_1984     //project clipped data to USA_Contiguous_Equidistant_Conic     try     {            geosubdiv = (IGeoDataset) psubdivftly.FeatureClass;            //get the geometry projected       geoclone = (IClone) clipgeo ;       IGeometry clipgrid = (IGeometry) geoclone.Clone();       // this polygon geometry is in GCS_North_American_1983                 ISpatialFilter spfilt = new SpatialFilter();       //project the polygon geometry to GCS_1984       clipgrid = ProjectFromEquidistantToGCS(clipgrid);                 spfilt.geometry = clipgrid       spfilt.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects       curtrack =  pGeoTrack.FeatureClass.Search(spfilt, true);            //create featurebuffer to insert clipped data into the input layer       IFeatureBuffer subdivfeatbuff =  memcliptrackly.FeatureClass.CreateFeatureBuffer();       IFeatureCursor insertcur = memcliptrackly.FeatureClass.Insert(true);            //check each feature is within or crosses and add to the virtual track layer.       while ((feattrack = curtrack.NextFeature ()) != null)       {         geotrack = feattrack.ShapeCopy;         IRelationalOperator relop = (IRelationalOperator)geotrack;         if (relop.Within(clipgrid) == true)         {           geotrack = projectFromGCSToEquidistant2(geotrack); //projects it  to USA_Contiguous_Equidistant_Conic           subdivfeatbuff.Shape = geotrack;           insertcur.InsertFeature(subdivfeatbuff);         }         else         {           if (relop.Crosses(clipgrid) == true)           {                       topo = (ITopologicalOperator)geotrack;             //get the clipped geometry             topo.Clip(clipgrid.Envelope);             // geoclone = (IClone)topo;             geotrack = projectFromGCSToEquidistant2((IGeometry)topo); //projects it  to USA_Contiguous_Equidistant_Conic             subdivfeatbuff.Shape = geotrack;             insertcur.InsertFeature(subdivfeatbuff);             }         }                 }                 insertcur.Flush();


title: Arcgis server not showing map service

text: I have a map service enabled in my server and have kml and wms as well as mapping and feature access capabilities turned on. I can view the wms in qgis, I can veiw the kml in google earth. However when I add the map as an operational layer to my config.xml in my sample application, I get no errors, the other basemap layers come up fine, and the operational layers show up in my widget (with my layer checked on), BUT it doesn't show up. I have changed my map service document to be the same coordinate system as the layer I am trying to show. The data is coming from postgresql db as an sde service. Any help would be great More info: I still have this problem. I have tried to add services that use different cs. I have reprojected my data to be the same as the basemap service in the flex sample app. I have tried to use a config.xml that only contained my service (and changed the xy extent to match. Still I can see it as wms and kml but not as a mapservice.


title: Run a python script every time the page is changed in Data Driven Pages? (ArcGIS 10)

text: I'm setting up a Data Driven Pages mapbook mxd for a series of parking locations. I have dynamic text that updates according to the current map sheet, but I also have text that needs to update using a script. I need the dynamic text to return the LotType attribute, but it is coded using a domain, so it returns an integer rather than the string associated with the integer (such as "Special Permit", "Fac/Staff", "Metered", etc). I have a VBscript from an old MPS Atlas MXD that returns the associated string, but I don't know of any way to use a script with non-MPS Atlas dynamic text. Is it possible to embed a python script into the MXD that runs whenever the user changes pages in the map book? Or is there a way to include scripts within the core Arc10 dynamic text? Thanks


title: How to show attributes of a new feature in Attribute table area when using IEditor3 and IEditTemplate, ArcGIS 10

text: ArcObjects SDK for .NET and ArcGIS 10 Desktop Quick snippet to set the stage:               IEditor3 editor.StartOperation();       IFeature newFeature;     // fill out newFeature.Shape, etc.     newFeature.Store();     editor.StopOperation("Adding new feature");      At this point I'd like to automatically "show" the attributes for this new feature in the Attribute display/table, which could be visible alongside the Create Features template list. Equivalent process through ArcMap itself:   * Start the Editor via the Editor Toolbar.   * Click a feature template from the **Create Features** list.   * Click the **Attributes** button to open the Attributes display.   * Mouse-down in ArcMap. See new/default attributes for the new feature be displayed in the Attributes table area. Is the possible with ArcObjects for .NET?


title: How to filter layer of osm data?

text: I have a table osm_roads that has a column named _geometry_ geometry geometry(LineString,900913). If I want to filter a city to draw districts, I cannot do it by geographic coordinates? What if in openJUMP or QGIS I want to select by bounding box, such a table as that? Thank you in advance.


title: How to convert GDB data to shapefiles without ArcGIS?

text: We've got some data (supposedly it's a simple collection of points--or at least that is all we need to extract) from a 3rd party in a GDB folder (I don't know how to tell what version) but we don't have an ESRI Arc* product. We are not primarily a GIS company, and our tools basically support only shapefiles. How can we convert the GDB to shapefiles? Note: I tried following the instructions to get GDB support in QGIS (after uninstalling QGIS and re-installing via osgeo4w-setup) but it didn't work; QGIS still says "Xyz.GDB is not a valid or recognized data source".


title: references for the colloquial naming of spatial data resolution

text: In our era of data deluge, many folks, including myself, quest for ill-defined "high-resolution" datasets. I suspect that colloquial description of spatial resolution as "high" or "low" is and will continue to be a set of terms relative to the what was available in the recent past. However, I'm wondering: are there peer-reviewed conventions for subjective description of spatial resolution? Perhaps some sort of Beaufort scale for spatial resolution?


What tool do I need to have satellite background images to draw features of a postgis database? I looked about Geomajas but seems very complex for what I'm trying to do. I have osm data imported in a postgis database, as a layer in QGIS. The second layer is what I am drawing, and the third layer that I ask here is how to have satellite imaginery on bottom of the two mentioned ones as a visual aid in drawing features.


title: ArcObjects Runtime Binding 10.0

text: I installed the ArcOjbects SDK (if I remember correctly) and in visual studio I have referenced the needed assemblies. Namely ESRI.ArcGIS.Geoprocessor. When in stepping through in debug mode I hit this routine:               void BindingArcGISRuntime(object sender, EventArgs e)             {                      //ProductCode[] supportedRuntimes = new ProductCode[] { ProductCode.Desktop };                 //foreach (ProductCode c in supportedRuntimes)                 //{                 //    if (RuntimeManager.Bind(c))                 //        return;                 //}                           if (!RuntimeManager.Bind(ProductCode.Desktop))                 {                     // Failed to bind, announce and force exit                     Console.WriteLine("Invalid ArcGIS runtime binding. Application will shut down.");                     System.Environment.Exit(0);                 }             }      But I don't hit the "Console.WriteLine("Invalid ArcGIS runtime binding. Application will shut down.");" piece Then straight to this auto-generated routine:               if (RuntimeManager.ActiveRuntime == null)                 {                     EventHandler temp = ResolveBindingEvent;                     if (temp != null)                     {                         temp(this, new EventArgs());                     }                 }                      return Initialize();      Any ideas? Thanks.


title: How to calculate stream order from vector data in ArcGIS 10

text: I know that NVision Solutions has a tool for calculating stream order from vector hydrography data for ArcGIS 9.3.1. Is there something similar that works with ArcGIS 10?


title: Symbolizing Layers from a PostGIS Database to many Client Computers

text: I was hoping someone could help with a symbolizing issue. I have a PostGIS database on a server with a number of layers. The layers of a series of polygons, lines and points that all share the same symbology. I want to symbolize the layers with the same symbol schema on computers in lab that are connected to the database. I have used the MultiQML tool to symbolize the layers and saved them on the server, which did not load the correct symboled features on the connected lab computers. Do I have to have the QML files on each lab computer, or is the a way to have a QML file on the server that each lab computer will be able to read via the database connection and symbolize the layers correctly without needing to store the QML locally? Thanks, Ryan


title: How to reverse geocode without Google?

text: I have a basic tracking system i've written with OpenLayers using google and openstreetmap and I want to use geocoding to get street names etc However, I'm a bit confused. I was hoping there was a easy/free way with OL and OSM however, i'm not sure there is I've tried using google, but it doesn't work with bulk data. It allows me to get about 3 or 4 back and then doesn't work for the rest..... i think they require a wait period between each request anyway, can anyone recommend a free way to geocode without google? please be patient with my lack on knowledge in this areas Thanks


title: Mosaicking a selection from a raster catalog and returning the mosaicked raster

text: I have pretty much given up on trying to get help from ESRI forums so I thought I would try my luck with the good people here... Here is what I have (in C#)... I have an IFeatureCursor that is giving me a selection from an IRasterCatalog (so basically it is a cursor full of IRasterCatalogItems). What I would like to do is take these rasters and mosaic them inside a method and return the mosaicked raster. I have tried what I thought would be obvious beginning with an IRasterCollection...               IMosaicRaster m_mosaic = (IMosaicRaster)m_rasColl;     m_mosaic.WhereClause = "";          ISaveAs m_saveAS = (ISaveAs)m_mosaic;     IRasterDataset rasDataset =             (IRasterDataset)m_saveAS.SaveAs("tmp_" + selectedFeatureClass.FeatureClassID.ToString(),                                            tmpWS, "TIFF");           IRaster extractedRaster = rasDataset.CreateDefaultRaster();          // Compute the statistics for the raster.     IRasterCalcStatsHistogram m_calcStats = new RasterCalcStatsHistogramClass();     IStatsHistogram m_stHisto = new StatsHistogramClass();          m_calcStats.ComputeFromRaster(extractedRaster, 1, m_stHisto);          return extractedRaster;      This gave me what seemed to be a properly mosaicked raster, and the values seemed to be correct when I used the identify tool in ArcMap, however, the display looked like a raster that has not had statistics calculated and I can't seem to figure out why. So I started looking at less obvious ways to do it but I haven't had much luck so far. Suggestions? Thanks a bunch!


title: How to get rpy2 for QGIS on OSX?

text: I am a new QGIS user, have installed QGIS on a Mac running Snow Leopard OSX 10.6.8 following all the instructions. QGIS seems to work well for everything other than when I boot up I get the message "Unable to load the plugin: Unable to load required package rpy2. Please ensure that both R, and the corresponding version of Rpy2 are correctly installed". I have installed Rpy2 several times and reloaded R. I have R version 2.15 loaded (but also RStudio). Does anyone have any ideas?


title: How to get FGDB support in GDAL 1.9 without compiling?

text: I was hoping that I could use the 'master' (nightly-build) version of QGIS to get fgdb support without compiliing (I'm not a developer and frankly have wasted days of productive time trying to get things to compile with no luck). Since it seems that gdal 1.9 support fgdb, I thought I could just download the developer version. However, the latest version of QGIS (191a229), even with gdal 1.9 does not seem to allow me to add my gdb file. Is the gdal version included with the master build just not compiled with fgdb support? Does anyone know if there's a way to do this without compiling a new version?


title: Georeferencing a whole heap of raster imagery

text: I have a whole heap of satelite imagery coming in that is not georeferenced. They have a grid overlaid so I do know the bounding box/coordinates of the images (and projection). I just need them georeferenced for use in ArcGIS. I have manually done a few by referencing to a point dataset I traced over the top - and it works great. I just need to automate this process. Each image is exactly the same, so the same points overlaid will always match up (even if it doesnt, the application of this does not need to be accurate). Basically I need to apply the 'link table' information from the example I've already done, and apply it to all the images... Possible in python/model builder?


title: Naming your Output feature class?

text: I am working in a work space in python on ArcGIS 10.1 on a remote desktop. Here is my code               import arcpy     from arcpy import env     env.workspace = "E:\Example"          PST = arcpy.ListFeatureClasses ("", "Point")     LST = arcpy.ListFeatureClasses ("", "Polyline")          for point in PST:         for polyline in LST:             arcpy.SplitLineAtPoint_management (polyline, point, point + "500" + ".shp", 1)      I want the output feature class to have the same file name as the point feature class, but with 500 appended to the end. Currently the way the code is now, it will name the output "Example.shp500.shp" which is giving me an error of failing to execute. How can I name my output "Example500.shp"?


title: Need software to create Info-graphic map

text: I want to make an info graphic map for history theme hence I am researching on Google for suitable software but not able to find suitable software and very less information regarding available related with info-graphic maps. We are using Open source platform ,while searching I found this map and looking similar kind of maps. Found following useful links : Making info-graphics in ArcMap Blog :Infographics, GIS, and Mapping So does anyone knows about info-graphic maps ? & its software/reference sites ?


title: How to get the result of a geometry service relation method?

text: I'm using ArcGIS 10.0 flex API 3.0. I'd like to test if two geometries are equal. For that, the relation method is called like follows:                       var relationParameters:RelationParameters = new RelationParameters();             relationParameters.geometries1 = [geom1];             relationParameters.geometries2 = [geom2];             relationParameters.spatialRelationship = RelationParameters.SPATIAL_REL_RELATION;             relationParameters.comparisonString = "IDENTICAL (G1, G2) = TRUE";             _geometryService.addEventListener(GeometryServiceEvent.RELATION_COMPLETE, onGeomCompared);             _geometryService.relation(relationParameters);      The service is supposed to return a boolean that tells if the two geometries are equal or not. Isn't it? The GeometryServiceEvent returned by the service doesn't contain any boolean and the result object of the event is empty. How can I get the response of the service for my request?


title: How to calculate distance in meter between QgsPoint objects in PyQgis?

text: I tried to calculate distance between QgsPoint objects in qgis, i used QgsPoint `StartPoint.sqrDist(EndPoint)`. but i don't find an expected result, i also tried               import math     math.sqrt(StartPoint.sqrDist(EndPoint))      but still getting a very high result. does somebody know a way how to calculate distance in meter between 2 Qgspoint objects.


title: Wrong object size after projecting spatial data from one coordinate system to another

text: I am having problems reprojecting items in ArcMap 10.1 from British National Grid To WGS84. (WKID: 27700 Authority: EPSG to WKID: 102100 Authority: ESRI) All Polygon Layers I project [I use Project (Data Management)] from BNG to WGS are basically twice as big in metres as they were originally. Even the Shape Area and Shape Length inside the attribute table changes. Any ideas on this would be greatly appreciated. Thank you.


title: Extracting random point coordinates using ArcGIS 10 model builder

text: I have a question about using model builder in ArcGIS 10.1. In my model, I am creating a processing extent, and within that I am creating a random point to be used in further processing. This model is iterative, and for each iteration, I am creating a new random point (effectively overwriting the values for the point).   What I would like to be able to do is to write the x-y coordinates of these random points off to a file before the next iteration is performed (either all together or as separate files for each point).   In a way, what I want is the inverse of the "Get Cell Value" tool and to be able to write that value to a file.   I guess it wouldn't be too difficult in python, but I was hoping to be able to do this in model builder.


title: How do I truncate a fieldname in model builder using field iterator

text: I have about 9,000 tables with rather long fieldnames, I need to chop off part of the fieldname, so for instance I want to change 'GP_CodeA82010' to 'A82010'. I am doing this using field iterator and the calculate field, but the calculate field does not seem to like the iteration, so for instance instead on 'A82010' it gives me an output value of '1'. I have tested a few python expression on the calculation field; Right( [FieldName], 6) and slice function- "%Name%"[7:13] Does anyone know a work around this?


title: "car Icon moving in a polyline from start to end destination"

text: I am new in google maps api programming and I want to make a map which is the same as this map example "http://econym.org.uk/gmap/example_cartrip2.htm)" does anybody here can help me to start and build this kind of application for android? tnx... :)


title: Is there a way to find out which pixel values are around the current pixel?

text: I would like to write a cycle which takes a current pixel, look at the surrounding pixels, choose the pixel with max/min value and skip to that pixel...and so on. Is there in arcpy such a function which can work with pixel's neighborhood?


title: Do ArcGIS clients that direct-connect to postgresql/postgis databases override configuration settings?

text: I went in to modify some settings on a Postgresql database via the PGAdmin config editor and some of the settings don't have their default values and they are higher because they've been overridden (e.g. `maintenance_work_mem` ). I know PGAdmin will do this. And that I can do it manually in PSQL. Do the ArcGIS clients ever override postgresql configuration settings? And if so, what triggers it? Is it tool-based? Or when the software starts up? Or is Esri entirely hands-off?


title: GeoTIFF Mosaic/Crop/Translate: Easiest Method in QGIS?

text: I routinely have to mosaic multiple GeoTIFFs together, extract a rectangular area, and save that area as a JPG with a world file. Typically I've been using commercial software for that (Geographic Imager) but would like to switch that task over to QGIS. I can eventually get the task completed using multiple tools, but the process takes much longer and involves creation of intermediate GeoTIFFs along the way. In Geographic Imager I simply load the images into Photoshop, crop, and SaveAs and I'm done. So what's the easiest, quickest workflow to:   1. Load multiple GeoTIFFs   2. Mosaic them together   3. Select an area to extract (preferably by drawing a box)   4. Save selected area out to JPG with accompanying world file Thanks!


title: How to convert a string from GetParameterAsText to uppercase?

text: This is my first question here and I'm sure it's ridiculously simple. But I have Googled and no luck. I am trying to make a tool in ArcToolbox using a Python script so users can input a street and lot number and it will zoom to the correct location/s in the cadastre. It's working but I have to type the name of the street in upper case in the toolbox and I'd like it work for uppercase or lower. I assume I should use something like .upper() but I don't know where to put it. I've tried and failed in a few places. Here's my snippet:               RoadName = "'" + arcpy.GetParameterAsText(0) + "'"


title: What are the parameters for my connection QGIS

text: I am new QGIS, I was wondering what the parameters are for the New MSSQL Connection in the dialog window? Provider/DNS: `? when local ?` tried leaving empty Host: `? when local ?` tried my computer user name/server instance name Database: `full path to file I assume` ![Dialog Window](http://i.stack.imgur.com/djIky.png) Error: ![Error](http://i.stack.imgur.com/j1GdA.png)


title: Is there a way to automatically select a feature when another feature in a different layer shares attribute data?

text: I have a mapset displaying foraging areas, each with a unique ID. I spatially joined the mapbooks layer, so each mapbook sheet identifies the foraging area it is depicting (1:1), so sheet 5 is foraging area 18. So when sheet 5 draws, I want it to select foraging area 18, which has a matching attribute. I want to do this so I can change the symbology of the target feature to make the map stand out more. I'm using ArcGIS 9.3.x with the mapbooks extension.


title: How can I remove dipoles from rasterlayer in QGIS 2.2?

text: I work with geomagnetic survey data. When convertet into a rasterlayer, every object with a high magnetic potential (for exampel an iron-nail on the surface) is shown as very sharp anomalie called dipole (a black and a white circle tied together). These dipoles disturb the further analysis of the rasterdata. Normally I mark them with dots in a second vectorlayer, but this is a very timeconsuming work. Is there any extension, with wich I can filter these dipoles directly from the rasterlayer?


title: How to map factors and their relative importance in ArcGIS 10.1 for Desktop?

text: I want to know how I can calculate a map from a set of factors (in the form of raster layers) and Their relative importance for example: slope 20% lithology 10% . . .


title: exporting segmentation results from ecognition to other software for anaylsis

text: I want to create vegetation indices in ERDAS but based on segmentation derived in eCognition. I have tried different ways of exporting my segmentation results to a TIFF file which i can further use for analysis, however I am not getting undesired results.


title: Bounding Box in CuRL Tool(REST API)

text: I have set my time interval using CuRL tool for automatic Geowebcache.. How to set bounding box in curl command?


title: gdalwarp needs bounding box

text: I try to mask a raster with gdalwarp and I always get the answer: cannot compute bounding box of cutline. what is wrong with my cutline shapefile? or what requirements should it have to work?


title: User defined Polygon disappears with featureadded and selectfeature

text:
I am trying to achieve the following functionality:   1. User draws a polygon.   2. Once complete, fetch features in that polygon( Using WFS).   3. On Select polygon, show the features that intersect that polygon. My issue is: Once user draws a polygon, I use featureadded to fetch the features within the polygon. Once fetched, I can see the details using alert with no issue. But the circle/polygon disappears once I start zooming in /out. I want the polygon to be still visible. Again on select, the polygon disappears as soon as the polygon is selected. I want the polygon to be visible at all times and only be removed when I delete it. I am using this example as reference : http://openlayers.org/dev/examples/draw-feature.html Any inputs on how I can have the polygon to not disappear? Here is the code I am using:               var map, drawControls, polygonLayer;     var url = "http://localhost:8080/geoserver/wms";     var wfsurl = "http://localhost:8080/geoserver/wfs";     var bounds = new OpenLayers.Bounds(-20037508.34, -20037508.34, 20037508.34,20037508.34);     var EPSG4326 = new OpenLayers.Projection("EPSG:4326");     var EPSG900913 = new OpenLayers.Projection("EPSG:900913");          function init() {      var options = {          controls: [],          maxExtent: bounds,          restrictedExtent: bounds,          maxExtent: bounds,          numZoomLevels: 4,          maxResolution: 156543.0339,          projection: EPSG900913,          units: "m",          displayProjection: EPSG4326           };      map = new OpenLayers.Map('map', options);      // Using OSM as base map      var layer = new OpenLayers.Layer.OSM(          "OpenStreet Maps");                states = new OpenLayers.Layer.WMS(          "NA states", url, {          layers: 'topp:states',          tiled: true,          styles: '',          transparent: true,          format: 'image/png'      }, {          isBaseLayer: false,          displayOutsideMaxExtent: false,          ratio: 1,          buffer: 0,          yx: {              'EPSG:4326': true          }      });           polygonLayer = new OpenLayers.Layer.Vector("Polygon Layer", {          displayInLayerSwitcher: false      });      circleLayer = new OpenLayers.Layer.Vector("Circle Layer", {          displayInLayerSwitcher: false      });      map.addLayers([layer, states, polygonLayer, circleLayer]);      map.addControl(new OpenLayers.Control.PanZoomBar({          position: new OpenLayers.Pixel(2, 15)      }));      map.setCenter((new OpenLayers.LonLat(4, 47).transform(EPSG4326, EPSG900913)), 3);      map.addControl(new OpenLayers.Control.Navigation());      map.addControl(new OpenLayers.Control.MousePosition());      map.addControl(new OpenLayers.Control.Scale($('scale')));      map.addControl(new OpenLayers.Control.ScaleLine({          geodesic: true      }));      map.addControl(new OpenLayers.Control.LayerSwitcher());           selectControl = new OpenLayers.Control.SelectFeature([polygonLayer, circleLayer], {          onSelect: onFeatureSelect,          // onUnselect: onFeatureUnselect      });      drawControls = {          polygon: new OpenLayers.Control.DrawFeature(polygonLayer,          OpenLayers.Handler.Polygon),          circle: new OpenLayers.Control.DrawFeature(circleLayer,          OpenLayers.Handler.RegularPolygon, {              handlerOptions: {                  sides: 40                   }          }),          select: selectControl      };      for (var key in drawControls) {          map.addControl(drawControls[key]);          drawControls[key].events.register('featureadded', drawControls[key], function (f) {              onFeatureAdded(f.feature);          });      }      map.events.register('click', map, function (e) {          var bounds = map.getExtent().clone();          bounds = bounds.transform(EPSG900913, EPSG4326);          mouseLoc = map.getLonLatFromPixel(e.xy);               OpenLayers.Request.GET({              url: url,              params: {                  REQUEST: "GetFeatureInfo",                  EXCEPTIONS: "application/vnd.ogc.se_xml",                  BBOX: bounds.toBBOX(),                  SERVICE: "WMS",                  INFO_FORMAT: 'text/html',                  QUERY_LAYERS: [map.layers[1].params.LAYERS],                  FEATURE_COUNT: 50,                  Layers: [map.layers[1].params.LAYERS],                  WIDTH: map.size.w,                  HEIGHT: map.size.h,                  x: Math.floor(e.xy.x),                  y: Math.floor(e.xy.y),                  buffer: 0,                            //format: format,                  styles: map.layers[0].params.STYLES,                  srs: map.layers[0].params.SRS              },              success: setHTML          });               OpenLayers.Event.stop(e);           });     }          function setHTML(response) {           document.getElementById('nodelist').innerHTML = response.responseText;          }            function toggleControl(element) {        for (key in drawControls) {            var control = drawControls[key];            if (element.value == key && element.checked) {                control.activate();            } else {                control.deactivate();            }        }        if (element.value == "delete") {          polygonLayer.destroyFeatures();          polygonLayer.redraw();          circleLayer.destroyFeatures();          circleLayer.redraw();      }           }          function onFeatureAdded(e) {           var wfsProtocol = new OpenLayers.Protocol.WFS({          url: wfsurl,          featureType: 'states',          featureNS: "http://www.openplans.org/topp",          srsName: EPSG4326,               geometryName: "the_geom"      });                var pfilter = new OpenLayers.Filter.Spatial({          type: OpenLayers.Filter.Spatial.INTERSECTS,          value: e.geometry.transform(EPSG900913, EPSG4326)      });           wfsProtocol.read({          filter: pfilter,          callback: processSpatialBboxQuery,          scope: new OpenLayers.Strategy.BBOX()                });                  }        function onFeatureSelect(e) {           var cql = "INTERSECTS(the_geom," + e.geometry + ")"      alert(cql);      states.mergeNewParams({          'CQL_FILTER': cql      });          }     // Process the features fetched. For the time being we are just showing  POS Name of first feature         function processSpatialBboxQuery(e) {      alert("hi");      };

             <input type="radio" name="type" value="none" id="noneToggle"                 onclick="toggleControl(this);" checked="checked" />                 Navigate
                 <input type="radio" name="type" value="polygon" id="polygonToggle"                 onclick="toggleControl(this);" />                 Draw Polygon
                 <input type="radio" name="type" value="circle" id="circle"                 onclick="toggleControl(this);" />                 Draw Circle
                 <input type="radio" name="type" value="select" id="selectToggle"                 onclick="toggleControl(this);" />                 Select feature
                 <input type="radio" name="type" value="delete"                 id="delete" onclick="toggleControl(this);" />                 Delete feature
                <input type="checkbox" name="allow-pan" value="allow-pan" id="allowPanCheckbox" checked=true onclick="allowPan(this);" />                 allow pan while drawing


title: Difficulty in running the script

text: I have developed a script tool for adding a layer in the current map document. It did not work.. The message is "Completed script addlayer...". But there were no layers added in the TOC. Another issue is when i type the same script in the python console window within Arcmap, the layer gets added. After adding the layer through console window, when i execute the script tool, it works and the layer gets added. I don't know what is the problem... pls help.. this is my script..               import arcpy     mxd = arcpy.mapping.MapDocument("CURRENT")     df = arcpy.mapping.ListDataFrames(mxd, "*")[0]     addLayer = arcpy.mapping.Layer(r"D:\Academic\M.E_thesis\data\ascii files\maxent\soil.lyr")     arcpy.mapping.AddLayer(df, addLayer, "BOTTOM")


Can Python toolboxes be used with ArcObjects GeoProcess.AddToolbox? I'm trying to execute a tool in a Python toolbox (.pyt) from .NET. First I call:               GP.AddToolbox(@"C:\path\to\MyToolbox.pyt")       And that call doesn't get an error. Then I build the parameter list and then I execute it like this:               result = (IGeoProcessorResult2)GP.Execute("ToolName", parameters, null);      This immediately throws the "Error HRESULT E_FAIL has been returned from a call to a COM component." error, so nothing is returned to result. My next step is to create a dummy script and build it into a .tbx toolbox and leave the .NET code the same to see if it's something else, but it did make me wonder if this is even possible. I couldn't fine any examples of using AddToolbox with a .pyt file.


title: Randomly shuffle values from a column into a series of new columns in a feature class in ArcGIS for Desktop?

text: I am trying to run a Monte Carlo simulation on a point feature class. In order to do so, I am trying to take the values from a column of observed data and randomly "shuffle" them between the points over n additional "simulated" columns. In my case, I'd like to add 199 new columns, each with the values of the the observed data randomly redistributed among the points. My observed data are binary (1 and 0) I have been able to do this in R previously and bring the table back into ArcGIS, but the tables have become so unwieldly that I am having difficulty bringing them back into ArcGIS, and would like to run the analysis in ArcGIS, rather than moving in and out. Does anyone have any idea how to do this? I have tried using cursor, but have failed.


title: Comparing two maps with nice web interface?

text: Is it possible in Leaflet or OpenLayers to compare two maps with a bar to slide between the two? Examples include:   * http://www.esri.com/services/disaster-response/hurricanes/hurricane-sandy-the-aftermap   * http://www.abc.net.au/news/specials/qld-floods/ (not a map) I know it's possible, I don't know how to do it though. Anyone got any ideas? Examples?


title: Problem in splitting and removing Polyline

text: ![Two polylines intersecting](http://i.stack.imgur.com/rD4Mj.png) I am stuck at a problem for weeks, Please Help!!! I am making program in C# Here is what I need to do --- I have two polylines that are intersecting with each other at some points. I want to delete those portion of one polyline which are going above the other polyline. I am attaching an image explaining the problem I am going through. In this image, I have two Polylines Line1 and Line2. So according to that image, I want to remove the **dashed** portions of Line2. Please Help!!! Thanks a ton in advance Here is what I have done till now.... First I am finding out all the points of intersection of the two lines and saving them in a layer.... Then I am Splitting these lines at these points using this code...               while (pPointF != null)     {         IPoint pPoint = (IPoint)pPointF.Shape;         ISpatialFilter pSF = new SpatialFilterClass();              pSF.Geometry = pPoint;         pSF.GeometryField = "Shape";         pSF.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects;              IFeatureCursor pLineCursor = pFeatureclass1.Search(pSF, true);         IFeature pLineF = pLineCursor.NextFeature();              while (pLineF != null)         {             IPolycurve pPolycurve = (IPolycurve)pLineF.Shape;             IPoint pToPoint = pPolycurve.ToPoint;             IPoint pFromPoint = pPolycurve.FromPoint;                  if (pFromPoint.X != pPoint.X && pFromPoint.Y != pPoint.Y && pToPoint.X != pPoint.X && pToPoint.Y != pPoint.Y)             {                 IFeatureEdit pFeatureEdit = (IFeatureEdit)pLineF;                 pFeatureEdit.Split(pPoint);             }             pLineF = pLineCursor.NextFeature();         }     }      The problem with this code that it is giving different results than expected. As I can see from attribute table, it is generating more than expected parts of line i.e. if there are 5 intersection points then there should be 6 parts of a line, nut it is creating 10-12 parts. Also the number of parts created for both lines are different, which logically should be same. Due to these problems, when I am trying to delete unwanted parts of line, the arcGis crashes..... Please help with this problem. . . . Now I am using a Geoprocessor to do splitting, following is the code:                       SplitLineAtPoint pSplitter = new SplitLineAtPoint();             pSplitter.in_features = pPolylayer1;             pSplitter.point_features = pPointsFeatureLayer;             string outfile = txtDirName.Text + "\\" + "splitted";             pSplitter.out_feature_class = @outfile;                  Geoprocessor GeoPro = new Geoprocessor();             GeoPro.Execute(pSplitter, null);      But on the last line, where Geoprocessor is being executed, the ArcGIS crashes... Please tell me where is the problem...


title: Replicating QGIS style in R with rgdal

text: I want to use rgdal to plot a shapefile graduating the colour of features according to the values of one field (`density`). I couldn't find online detailed indications on how to do it and the package documentation doesn't explain the options for `plot()`. This is what I composed so far               library(rgdal)            shapefile = "shapefile.shp"     map <- readOGR(dsn = shapefile, layer = "shapefile", verbose=FALSE)          plot(map, xlim = c(6.70, 18.32), ylim = c(35.2, 47.6), border=NA, add=TRUE)      This is the style in QGIS I am trying to replicate: ![enter image description here](http://i.stack.imgur.com/hihky.png)


title: Field calculator fails when computing First X coordinate in ArcGIS10

text: I have this line in python to compute First X coordinates into the Field Calculator from each polyline:               !Shape!.firstPoint.X      Then, when I inspected the attributes I see this rare thing: ![Example](http://i.stack.imgur.com/JAiwX.png) You may see that I have two fields ( _StartX_ and _InicioX_ ), both with the same Python code, one is a string field and the other a Float one. But **the computed VALUE is wrong** :( **You can clearly see that "spatially" the features are situated in different positions, but on its computed X Value are the same.** **Could someone explain me that?** Does anybody know any other way to compute this X,Y start and end coordinates from lines? Thanks for your time.


title: Is it possible to connect to FEMA's WMS and WFS using QGIS?

text: I have just found out about the wonderful world of WMS and WFS. I have added the FEMA versions of both to QGIS. WMS and WFS. The WMS seems to connect and adds something to QGIS, yet it is not visible. The WFS connects, but I get an error that it is not a valid layer. So... I guess that my question is - Is anyone willing to test these and maybe give some guidance on how to better use these? I have also tried NYS WMS and the Nationwide Inventory WMS `http://107.20.228.18/ArcGIS/services/FWS_Wetlands_WMS/mapserver/wmsserver?` with no luck. Similar results to FEMA's WMS.


title: Counting road intersections and blocks

text: I have a shp file that consists of the road system of Mexico City. I would like to count the number of road intersections (crossings) and blocks (which I would define as a segment starting from one intersection and ending in another). This within a buffer of a certain size: a circle centered in a polygon within another shp file. I am mainly using ArcGIS but am quite familiar with R, Python and QGIS if the solution is simpler on any of those. I realize their are a few similar questions around here but havent been able to use any of the information provided there.


title: Layers with same coordinates not lining up [QGIS]

text: I have a project where I have a few layers (State of Texas, cities, roads, etc.). They are all resaved and projected to NAD83 / UTM zone 14N EPSG: 26914. The project CRS is also NAD83 / UTM zone 14N EPSG: 26914 and on the fly projection is enabled. I added a text delimited file and cannot get them to line up. I checked the attribute table of the delimited file and the coordinates are correct. I cross-checked this with a point close to where a few of the points from the delimited file are and they are very close so in theory they should be lining up. I tried to save the delimited file and save it with the same projection and add the shp file back in and that did not work. Anyone have any ideas? I can't think of anything else to try.


title: pyscripter out of memory

text: I'm running a loop that create feature classes at a geodatabase. After creating 221 datasets I get the error message "pyscripter out of memory". Is there a way to clear the physical memory of computer in order to resume with the script ? I've alredy tried :                 gc.enable       clear = lambda: os.system('cls')       clear()       gc.collect()      ![Windos Task Manager](http://i.stack.imgur.com/jTEJp.png) ![enter image description here](http://i.stack.imgur.com/5KQhI.png) ![enter image description here](http://i.stack.imgur.com/Eh94I.png) ![Out of Memory in pyscripter](http://i.stack.imgur.com/UBUV9.png)


title: Measuring distance between discharge point, along river, until next intersection with another river?

text: i hope someone can answer me my simple question. I have two sets of data   1. a river network --> line feature   2. discharge points along those rivers --> Point Data The discharge points represent wastewater treatment plants. I try to find those treatment plants which don't have a inflow of another river, for a minimum of 2 km downstream. This means the river from this point on does not intersect with another river for at least 2km. How can i find those plants? Are there any tools for this? What about the Arc Hydro tools?


title: CartoDB table seems to be corrupt

text: Have 218 rows in a table in the free version of CartoDB. Suddenly behaving strangely - not accepting edits to fields in existing rows, only showing 200 rows when sorted Descending on row number (but showing 218 when Ascending), etc. Is there a way to refresh the database or at least find out where the problem is?


title: how to refresh a layer using openlayer

text: in geoserver I configured a layer containing positions that change in time (they are smart phone positions indeed). I'am able to insert fresh position in geoserver datastore (a postgis db), but I'm not able to automatically reload the layer in client side using openLayers. In addition I wish to load from the layer only the information contained in the visualized BBOX (ie what the user will see on its screen). Could anyone provide me a suggestion?


title: Develop 3D buildings with roof structure from raw data

text: I recently received a task to develop 3D buildings with roof structure. All I have now is satellite images and aerial photographs covering the entire urban area. After reading some previous methods, it seems possible to achieve this with building/roof edges delineated and height information (e.g. a digital surface model). Is it possible to perform these tasks in ArcGIS? If so, with which tools? I am pretty new to GIS so I am not familiar with tools involved in such tasks.


title: GeoServer setup for intercepting WFS-T requests and handle them via PHP

text: I'm developing a web app which uses GeoServer side by side with a PHP framework. My goal is to be able to intercept WFS-T requests to update some records whenever a WFS transaction is executed (insert/update/delete). Eg: I want to be able to keep track of modifications to the features saving `user_id` of who made the change and duplicating the affected record before the changes take place to keep an history of changes on features. My actual setup is as follows:   * apache proxies all `/geoserver/*` requests to port 8080 (Tomcat/GeoServer) except for those directed to `/geoserver/wfs/*`   * these last (wfs) are managed by my php framework (they are first checked for session authentication) which reads the request XML body and backs up the interested records accordingly before passing on the request via php cURL I'm quite satisfied of the working for now but...   * I'm not sure it's a good approach (seems a bit complex): I wonder if there's a built in way of managing record backups in GeoServer   * I wonder if there could be security implications in leaving all requests (except for WFS ones) go directly to GeoServer Any suggestion either on php or on GeoServer would be greatly appreciated.


title: Binding licence for using ArcEngine

text: I have just started developing with arcengine and some how I am stick in first steps.. When I want to bind licence to start appplication I get an error which means I haven't any licence..(I think so)..Here is program.cs code                   using System;     using System.Collections.Generic;     using System.Linq;     using System.Windows.Forms;     using ESRI.ArcGIS;     namespace WindowsFormsApplication4     {         static class Program         {             /// <summary>             /// The main entry point for the application.             /// </summary>             [STAThread]             static void Main()             {                // ESRI.ArcGIS.RuntimeManager.BindLicense(ProductCode.EngineOrDesktop);                           if (!RuntimeManager.Bind(ProductCode.EngineOrDesktop))                 {                     MessageBox.Show(                         "Unable to bind to ArcGIS runtime. Application will be shut down.");                     return;                 }                 Application.EnableVisualStyles();                 Application.SetCompatibleTextRenderingDefault(false);                 Application.Run(new Form1());             }         }     }      I checked Arcgis Administrator and here is licences that I have..Can you please help me about binding licence?![ARCGIS Licence ](http://i.stack.imgur.com/EbYEt.png)


title: Parcel ID Mapping (brainstorming)

text: In my project I have a table full of residential properties by Parcel_ID. Table 2 is a TIGER streets shapefile. I would like to build a query in MS Access displaying all properties within a mile radius of my selected property from table 1. ![project structure](http://i.stack.imgur.com/icQDh.jpg) My questions are :   1. Should I use open-source postgis (spatial databases) and build this project in SQL 2008 to query off map location rather than tabular data?   2. Or Should I use a FileGeodatabase to import the shape file into MS ACCESS?  Thanks for any help or reply, PS - I have the shapefile loaded in ArcGIS explorer - but not sure how to export. 6/17/2012 update* _*_ I posted the state parcel tables below: **County Mapping Parcel Table**   OBJECTID PARCEL_TAG **PARCEL_I** PARCEL_C LOTNUM SUBDIVISIO SUBDIVNUM STNUMBER PRE_DIR STREET_NAM SUFFIX SUF_DIR FULL_STNAM CITY COUNTY STATE TOWNSHIP ZIPCODE ESTSQFT BLOCK MODDATE STATUS **SHAPE_AREA** **SHAPE_LEN** | 1 1251683 8011477 8011477 7809 1723 4015 N KENWOOD AVE KENWOOD AVE INDIANAPOLIS MARION IN WASHINGTON 46208 5280 1004 **** _5279.952157_ *** **_343.9990607_**


title: How to use a spectral library

text: I want to classify an Aster image using the ASTER Spectral Library - Version 2.0 Anyone here ever done this? Any suggestions? Note: I'm using opensource GIS and Remote Sensing software Thanks.


title: Why does the File Geodatabase that I created not appear to any other users on the network?

text: I created a File Geodatabase to store files that I use in maps for my coworkers. The folder is in a location where all users should have access and permissions. Strangely enough, after I created the gdb, I'm the only one that can see or access it. It appears in my Catalog view, but not in windows explorer. The map document that accesses the fcs in the gdb works fine for me but shows broken links for my colleagues. Any reason why this might be happening?? ![Where's my database?](http://i.stack.imgur.com/7YPlW.png)


title: How to Export a Shapefile as KML Using R

text: All, I have exported an ESRI shapefile (polygon) to KML using R, but the resulting layer shows only the outer extent of the polygon. How can I fill in the area with a color? Any help is greatly appreciated. The code is as follows:               library(rgdal)     setwd("G:\\GIS_SJR\\GIS Shapefiles\\Boundaries")     polygon <- readOGR(".", "20km_buffer_albers")     polygonWGS <- spTransform(polygon, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))       writeOGR(polygonWGS, dsn="polygonWGS.kml", layer="polygonWGS", driver="KML")


title: How to add a graphic table element in Python without Production Mapping extension?

text: I'd like to add a table in the Map Layout that includes only items visible at the current extent of the dataset. I know that the Production Mapping Extension lets you insert a Graphic Table Element filtered by the dataset current extent (I suppose). But we haven't this extension. Has anyone ever tried to do something like this in Python?


title: Matching Landsat 8 with previous Landsat images

text: I am working with Landsat 8 in ENVI. My image's dimensions are 8921 x 13131. I need to apply a forest mask to this image with dimensions of 8805 x 12321. I try to resize an image, but it doesn't work. (After that I need to calculate forest disturbance index for 2013 and subtract it from DI for 2000, so I really need an 2013 image (Landsat 8) to be resized). I resize it using basic tool - spatial subset - subset using - file - and then pick my forest mask. However, as a result I get a file of 8806 x 12321 instead of 8805. If I adjust the number of samples manually, it gives me 8805. But then I try to apply a mask and that's what I get (see a screenshot). They just don't match. And when I try to overlay my mask (which I made based on Landsat 7 image) with my Landsat 8 image, those two images have different frames. Does anybody have any good advice? Thanks a lot. ![enter image description here](http://i.stack.imgur.com/OOFMu.png)


title: Simplify Multigeometry polygons in QGIS

text: I am really new to QGIS and KML so this might be obvious but I am too ignorant to know what to do. I have some KML files I would like to use in an application but they are much more complicated than what I need. They contain multigeometry with a number of polygons and I would like to merge these regions to simply get the outline of the total region. In addition, I would like to reduce the number of vertices. ![Displaying layer](http://i.stack.imgur.com/ATYxc.png) I saved this as a shapefile so I could access the Merge Selected Features tool and I merged all of the features: ![Merged features](http://i.stack.imgur.com/fiGlh.png) But when I try to simplify the feature I get the error: > This feature cannot be simplified. Check if it has enough vertices to be > simplified In the documentation for Simplify Feature I saw it mentioned: > If a feature cannot be simplified (e.g. MultiPolygons), a message shows up. So I thought this was likely the problem. Is there any way I can simplify this, or convert the multipolygon to a simple polygon? Edit: During my searches I found mention of a plugin called fTools that seems to have some possibly useful functions, such as "Multipart to singleparts" and other vector tools that would be useful, however I do not seem to have it by default in QGIS 2.2, nor can I find it on the qgis plugins page. Is this plugin gone?


title: WFS protocol - Update method

text: I have a question related to the Update method, in case the feature doesn't exists, the update method will generate a new feature? (i.e. will behave as Insert)? Or maybe, in contrast, an insertion using an existing fid (feature id) will perform an update over this feature? Thanks in advance.


title: evt.layers undefined while trying to display Legend widget in Javascript API

text: I cannot get an ArcGIS Javascript API Legend widget to show up. An excerpt of my code is as follows. I believe I'm failing somewhere near the console.log that I call...undefined comes up in the console at this point. So for some reason evt.layers is not getting populated. Have I provided enough code? If you need more snippets I can copy/paste more.                       parser.parse();                  //initialize the map object             map = new Map("map", {                 basemap:"gray",                 center: [-96.53, 38.374],                 zoom: 13             });                  //access the hydro layer as a feature layer             var hydro = new FeatureLayer("http://sampleserver3.arcgisonline.com/ArcGIS/rest/services/Hydrography/Watershed173811/MapServer/0", {                 mode: FeatureLayer.MODE_ONDEMAND,                 outFields: ["*"]             });                  //add hydro to map             map.on("layer-add-result", function (evt) {                 console.log(evt.layers);                 var layerInfo = arrayUtils.map(evt.layers, function (layer, index) {                     return {layer:layer.layer, title:layer.layer.name};                 });                 if (layerInfo.length > 0) {                     var legendDijit = new Legend({                         map: map,                         layerInfos: layerInfo                     }, "legendDiv");                     legendDijit.startup();                 }             });             map.addLayers([hydro]);         });