title: How to get arcpy.AddMessage() to display a message from a script tool called from another script tool?

text: I am using a python script tool within an ArcToolbox in ArcGIS 10.2. The script calls another python script tool while running but does not display the messages from the called script when using the arcpy.AddMessage() method. Example: scribtbox (ArcToolbox) two script tools   * script1   * script2 script1:               import arcpy          # output to the screen     arcpy.AddMessage("Doing routine 1...")     routine1()          # call script2     arcpy.script2_scriptbox      script2:               import arcpy          # output to the screen     arcpy.AddMessage("Doing routine 2...")     routine2()      The output model run dialog box in ArcGIS says:               > Doing routine 1...      and then finishes without saying the message from script2. I'd like it to say:               > Doing routine 1...     > Doing routine 2...      and then finish. Is this functionality not built in, am I doing something wrong, or is there some other issue at hand. Thanks.


title: SAGA-GIS Produces different results in command line than GUI

text: I have been working on a lidar processing workflow using SAGA-GIS. I determined my workflow on a single tile in the SAGA GUI and now I need to batch the operations on 100 more tiles. I wrote a batch file to call the SAGA Command Line to run all the las files in a folder. I noticed, though, that the out DSM created by the command line process looks much more blurred and interpolated than the DSM I created in the GUI for the same tile. Then, I started to investigate, so I started running each individual process of my workflow in the GUI and in the command line with the same data. The very first step is to import the LAS file. In the GUI, this doesn't save to a .spc file, it just seems to do it in memory. With the command line, it will save a .spc file from the LAS file. These two files are different: ![Pictures of Results and Min/Max Values](http://i.stack.imgur.com/IZVWv.png) The picture in the upper left is the point cloud created when I imported the LAS file in the SAGA GUI. Notice the values circled in the lower left image. Those are the Min/Max values of the upper left point cloud. The picture in the upper right is the SPC file that was created when running the Import LAS module via the command line. Notice the Min/Max values circled in the bottom right. As far as I can tell, all the parameters I used were identical. In fact, I ran the GUI version first and then right-clicked on the module to create the batch file with all the same parameters. Can anyone help me with this? The DSM's created by the batch file are not nearly as clear as the one created by the GUI because they require more interpolation and cannot be used for identifying tree canopy.


title: How to export QGIS contour attributes in a DXF file?

text: I'm new to this. I have a .tiff data file that generates contours in QGis. That's super. I save it as .dxf. The contours don't have any attributes or data. That's bad, because I want to make a 3D landscape in Revit using the information. How do a save/export the contour lines from QGis so that they have a Z-elevation attribute or at least some information... Thankful for help!!! * * * I got an answer, but the post was changed so i don't know where it went... Anyhow... I ticked the attributes box when extracting the contours, but the attributes doesn't show in the dxf...


title: "Error HRESULT E_FAIL has been returned from a call to a COM component." during using geoprocessor by VB.net

text: I am also get a error shows like "Error HRESULT E_FAIL has been returned from a call to a COM component." And the Exception snapshot shows the "stacktrace" is at "ESRI.ArcGIS.Geoprocessing.GeoProcessorClass.Execute". Can you help me to solve it? The thing is that I can using the geoprocessor by C#, but can't using it by vb.net. Is it something wrong with the setting or .NET framework? Also, the Com Exception is -2147467259! Private Sub RunBuffer()                   Dim gp As IGeoProcessor2 = New GeoProcessor()         gp.OverwriteOutput = True         Dim sev As Object = Nothing              Try             Dim parameters As IVariantArray = New VarArray()             gp.AddToolbox("C:\Program Files (x86)\ArcGIS\Desktop10.0\ArcToolbox\Toolboxes\Data Management Tools.tbx")             gp.SetEnvironmentValue("workspace", "C:\gp\AirportsAndGolf\golf.gdb")                  parameters.Add("C:\gp\AirportsAndGolf\golf.gdb\Airports")             parameters.Add("C:\gp\AirportsAndGolf\golf.gdb\Airports_buff2")             parameters.Add("1000 meters")                  gp.Execute("Buffer_analysis", parameters, Nothing)                 Console.WriteLine(gp.GetMessages(sev))              Catch ex As Exception             ' Print generic exception messages.             Console.WriteLine(ex.Message)             ' Print geoprocessing execution error messages.             Console.WriteLine(gp.GetMessages(sev))         End Try      End Sub


title: Can't insert AI CS6 files as objects into ArcMap 10.1

text: In the past I have inserted Adobe Illustrator files into ArcMap as objects. My old version of AI was 8.0 and ArcMap was 9.2. Now I have arcmap 10.1 and AI CS6. When the process doesn't shut down arcmap, it gives the message.."failed to create object. Make sure object is entered in system registry. ???????


title: geoserver/openlayers: multiple layer geotiff

text: I have a geotiff file with multiple layers (46 layers), each layer corresponding to a different date (the geotiff layers corresponding to sea surface temperature, 1 layer per year). I imported this geotiff in geoserver, but I can't render it with the geoserver visualisation tool (based on openlayers). 1./ does geoserver support multiple layers > 3 bands? 2./ if yes, when I query the WMS in openlayers, how can I select a given layer?


title: DEM error when trying to do topographic correction in GRASS

text: I'm new to GRASS gis and I'm trying to do a topographic correction (i.topo.corr) of 9 ASTER bands using SRTM dem. I keep getting this error message: > ERROR: Illumination model is of CELL type Is this error referring to my elevation raster, and if so how do I change the type to the appropriate one?


title: Export Web Map Task: Configure HTTP response headers

text: Is there any way to configure the http headers (response headers) of a geoprossessing task hosted by ArcGIS for Server? I would like to make the export web map task, return an image which would be directly downloaded by my browser.


title: Using Geoserver

text: I am new to GIS and interested in using Geoserver. I have had a look at tutorials and read some books on this but still feeling clueless. I intend to use C# to accomplish some tasks eg, plot points on a map, draw polygons from points of interest and show heat maps from such points. Looking forward to getting some pointers. Thank you.


title: How do I not specify a CustomerWeightField when scripting CustomerDerivedTA_ba?

text: When you run the Business Analysis toolbox wizard for creating a Customer Derived Trade area, you have the option of weighting each customer by count or for specifying a field (such as sales or number of visits) to use as a weight. You select this with radio buttons and if you select "count", then you don't have to specify a field. However, when scripting this tool, the CustomerWeightField is a required field. Leaving it blank yields a "This field was not found in the customer layer" type error message. In my case, I just want to use the count as the customers I'm analyzing do not have any other valid fields for weighting the results. Is there a magic value to enter for CustomerWeightField that will force it to use count? The documentation says: > The field used to calculate the trade areas. This is based on either the > number of customers (count) or the calculated weighted value assigned to > each customer. A fallback, I guess, would be to add a weight field and set it to 1, but I don't want to have to rely on modifying the input data.


title: co-registering a vector file to a georeferenced raster image

text: i would like to co-register a vector file to a raster image by a simple translation in envi. Is it possible to do this in envi? I have tried using the georeferencing tool in arcgis but the vector layer cannot be selected for georeferencing!


title: Convert String to Integer in QGIS tables

text: I want it to know how to convert string field values to integer values in QGis Dufour... there are numbers and strings in the same field, when I use toreal(string) it stops and it won´t create the field. Thanks in advance for your help! CIR


title: What do you think is happening to my QGIS Georeferencer

text: I'm pretty new with QGis, but have used tools as Openlayers and Geoserver for some years now without never being asked to diplay a georeferenced raster image. Which gets me to use QGis in order to build a Geotiff from a simple scanned PNG or JPEG. I've followed many tutorials, but I'm somehow still struggling to import any image file in my Georeferencer plugin window. Here's what I get when I import a jpeg image: ![enter image description here](http://i.stack.imgur.com/G24ZN.png) It looks like a badly read jpeg. Any idea ? Thanks for your help ! **EDIT** I've tried with a PNG by converting the original jpeg with an image tool (like IrfanView) and the image is now diplayed. It's better, but doesn't answer exactly why the jpegs aren't working ...


title: OpenLayers LineString: Not displaying the line (array of points)

text: I'm perplexed as to why tail_path is not drawing the line (can't see anything on the map) but the alert seems to be returning the correct data.               var path_travelled = new Array();     var path = json_coords.path_data;     for (i in path)       {                                    var lat = path[i].LINE_LAT;         var long = path[i].LINE_LONG;         path_travelled.push(new OpenLayers.Geometry.Point(long, lat));         alert(path_travelled);       }     var tail_path = new OpenLayers.Geometry.LineString(path_travelled);     vectors.addFeatures([tail_path]);      The alert for the path_travelled variable returns this (it adds a new POINT every time you click the OK button): ![enter image description here](http://i.stack.imgur.com/cQz1N.jpg) I'm using vectors.addFeatures (I have another variable I also use here but doesn't apply to this question so just removed it to avoid confusion). Any input on the matter is appreciated.


title: Hit test for dynamic object

text: I am drawing a dynamic moving vehicle like this:               m_dynamicSymbolProperties.set_DynamicGlyph(esriDynamicSymbolType.esriDSymbolMarker, this.bluePoliceCarGlyph);     m_dynamicSymbolProperties.SetScale(esriDynamicSymbolType.esriDSymbolMarker, .3f, .3f);     currentUnitPosition = getNextPosition();     dynamicDisplay.DrawMarker(currentUnitPosition);      Does anyone know an easy way to do a hit test on my vehicle? I am having trouble with this code example, http://help.arcgis.com/en/sdk/10.0/arcobjects_net/conceptualhelp/index.html#/How_to_identify_objects_in_dynamic_display/0001000004s2000000/ * * * One of the issues I have with the example is that ESRI states: _ArcGIS versions 9.2 and later have the following three display systems that use OpenGL as their underlying drawing pipe—ArcGlobe/GlobeControl, ArcScene/SceneControl, and MapControl's dynamic display._ So I am trying to get a hit test for a dynamic display object, but according to the tutorial, part 1 (link above), I have to use the OpenGL display system to give my dynamic display object a watermark in order for the hit test to work. This makes no sense at all.               uint ID = Convert.ToUInt32(r[0]);     OpenGL.GL.glLoadName(ID);      * * * Another issue I have is when this code is called               IArray dynamicHitArray = m_dynamicDisplay.Locate(x, y,         esriDynamicSelectionMode.esriDSMLayers, this as IDynamicLayer);      Another event is fired off to DrawDynamicLayer and I will get a stackoverflow unless I do this               if(!someBool)     {         someBool = true;          IArray dynamicHitArray = m_dynamicDisplay.Locate(x, y,             esriDynamicSelectionMode.esriDSMLayers, this as IDynamicLayer);            someBool = false;      }


title: Is there a function or tool that will add vertices to polylines where endpoints are snapped to edges? (ArcMap)

text: My friend and I are creating a walking network, and we have been having issues getting the routing to work. We just realized that the reason why our routes are messed up is because we snapped our routes to edges as well as endpoints. The routing works fine where the lines are snapped at endpoints because they share a common vertex, but when snapped to edges only one line has a vertex at that location, and the network doesn't see them as connected. Is there a tool or function that will automatically create a matching vertex on the other line?


title: Top-of-atmosphere reflectance conversion with WorldView-2 images

text: I recently received a set of radiometrically-corrected dynamic range adjusted 8-byte WorldView-2 images. I want to convert to top-of-atmosphere reflectance and would like to use the FLAASH in ENVI 4.8 to do so. However, in reviewing the WorldView-2 radiometric correction information provided by Digital Globe I noted that the report states that "Caution is advised when applying the equations provided here to pan-sharpened products, dynamic range adjusted (DRA) products, or WorldView-2 mosaics with radiometric balancing because the generation of these products may apply non-linear transformations to the pixel DN values." As the images I received were dynamic range adjusted radiometric balanced, my question is what care exactly should be taken when converting to TOA reflectance? Some of my photos are mosaiced with images from separate dates so I assume that I use the spatial subset function in FLAASH to limit the atmospheric correction to each date. What am I missing?


title: How do I identify and gracefully kill hung Gp service processes?

text: I've published a gp tool as an asynchronous non-pooled gp service. When the gp tool hangs I would like to allow the user to kill it, without impacting other services. If I put a webservice on the SOM that kills a process based on an ID, will it harm other services? Related questions: Is there a way for a webservice running on the SOM machine to kill a process on a (different) SOC machine? The first thing the gp tool does is write out the hostID and processID to a log file, but I'd like a system to work for any gp service, not just ones that log this info. Is there a request I can make to the gp server that returns a hostID & processID for an arbitrary jobID? **Update:** I've decided to first implement a method to Cancel GP jobs, and then only pursue the process killing strategy if canceling proves ineffective. Since the REST API doesn't currently have a Cancel method, I will be doing this via SOAP. I've also been shown how to log the process IDs, so that way I won't need to rely on my own logging.


title: Recoloring geotiff for visual presentation

text: I have a rather big Geotiff. It is the mosaic of four RapidEye satellite images. They are all a bit bland/grayish. When combined I would like to freshen up their colors so that they are more visually appealing. The imagery will be used as a backdrop for other data. I have tried right-clicking on the raster in the layers list and selecting Properties. Using Symbology allows for detailed adjustments. But there is no easy way to auto-adjust the image color balance. Some image software like IrfanView has "auto adjust colors". IrfanView can not open my geotiff image, but this kind of functionality represents what I am looking for. To the extent I am doing analysis on the imagery I will of course not mess with the color balance.


title: How to download Leaflet vectors in GPX format?

text: I would like to use the leaflet.draw plugin to allow users to digitize vectors on the map. At the end the data should be downloadable as a GPX-file. Is there any sample that shows how this can be done?


title: Generating GeoJSON with Python

text: Hi I want to programatically create a GeoJSON file using polygons from a shapefile but adding attributes from my own application. This is easily done for a shapefile:               def create_data_dayer(self,varlist, data):         """         Creates a new shape to contain data about nodes.         varlist is the list of fields names associated with         the nodes.         data is a list of lists whose first element is the geocode         and the remaining elements are values of the fields, in the         same order as they appear in varlist.         """         if os.path.exists(os.path.join(self.outdir,'Data.shp')):             os.remove(os.path.join(self.outdir,'Data.shp'))             os.remove(os.path.join(self.outdir,'Data.shx'))             os.remove(os.path.join(self.outdir,'Data.dbf'))         # Creates a new shape file to hold the data         if not self.datasource:             dsd = self.driver.CreateDataSource(os.path.join(self.outdir,'Data.shp'))             self.datasource = dsd             dl = dsd.CreateLayer("sim_results",geom_type=ogr.wkbPolygon)         #Create the fields         fi1 = ogr.FieldDefn("geocode",field_type=ogr.OFTInteger)         dl.CreateField(fi1)         for v in varlist:             #print "creating data fields"             fi = ogr.FieldDefn(v,field_type=ogr.OFTString)             fi.SetPrecision(12)             dl.CreateField(fi)              #Add the features (points)         for n,l in enumerate(data):             #Iterate over the lines of the data matrix.             gc = l[0]             try:                 geom = self.geomdict[gc]                 if geom.GetGeometryType() != 3: continue                 #print geom.GetGeometryCount()                 fe = ogr.Feature(dl.GetLayerDefn())                 fe.SetField('geocode',gc)                 for v,d in zip (varlist,l[1:]):                     #print v,d                     fe.SetField(v,str(d))                 #Add the geometry                 #print "cloning geometry"                 clone = geom.Clone()                 #print geom                 #print "setting geometry"                 fe.SetGeometry(clone)                 #print "creating geom"                 dl.CreateFeature(fe)             except: #Geocode not in polygon dictionary                 pass             dl.SyncToDisk()      since I have all the geometries on a dictionary by geocode (self.geomdict) I simply create the features, set the fields and clone the geometries from pre- existing layer (code loading that layer ommited for simplicity). All I need now is a way to generate the GeoJSON from the combination of fields and geometries, naturally with the help of OGR to get the rest of the file right (CRS, etc. as from the source map) How do export the feature collection generated as above?


title: How to antialiase tiles when seeding a layer from an GeoTiff in GeoServer

text: I have an GeoServer running that uses GeoWebCache. As a data source I use a 1.5 GB big map for a layer. The outcome of the tiling and scaling are very rough looking tiles. Is there a way to antialiase the tiles, that they are more pleasant to look at?


title: QGIS: rotatable markers with offset from point

text: Im wondering, if its possible (in QGIS), to set a rotatable marker with offset from point, like in my example below. Any suggestion? ![enter image description here](http://i.stack.imgur.com/fvOF4.png)


title: I am looking for a list of mapping libraries?

text: I am looking for the top APIs that I can use to display maps from a number of map servers/sources. (like openlayers and openscales). So are there any other popular APIs for **other languages than javascript and flex** (C++,Java,C#... ) ? (it doesn't matter if it is free or not)


title: rpy2 installation ERROR "Tried to guess R's HOME but no R command in the PATH"?

text: I tried to install rpy2 (pip install rpy2). I want to use the SDAPP4 extension with QGIS (QGIS shows up the message missing rpy2...). The system tells me always the same error: `Tried to guess R's HOME but no R command in the PATH`. What can I do to solve the problem? I'm using:               Windows XP     Python 2.7     R 2.11.0 (directory: C:\Programme\R\R-2.11.0)     QGIS 1.8      Sorry, I'm a beginner. Thanx!


title: What is the latest version of R that supported gpclib?

text: Does anyone know what the latest version of R that supported the gpclib package? I was using the latest version but when I tried to install gpclib, it said that it was unable due to my version of R. The gpclib page gives the minimum version and requirements. Do I have to downgrade my version of R?


title: 10.1 Marker Symbol Rotation Problem

text: I am trying to create a marker symbol (and then save the symbol as a style), but the origin of the items makes them rotate improperly (see pictures). Is there any way to merge all of the layers of the symbol? Or modify the origin of rotation? ![Before rotation](http://i.stack.imgur.com/e4whK.png) ![After Rotation - Incorrect](http://i.stack.imgur.com/2Dmg7.png) EDIT FOR COMMENTS: I am using the simple marker with multiple layers. The circle is the default circle with a second smaller white one added on top to create the label box. The main problem is the triangle should rotate around the origin of the circle, but I don't know how to modify this is the style editor. Thoughts?


title: Converting geometry_dump to Point geometry vertex

text: If this was a Point geometry `010100000001DE02098A1B52C090A0F831E62E4640`, is it possible to represent the above geometry_dump Point Vertex. I am trying with `St_collection`...,


title: Authentication in Geoserver 2.4.x from a web application

text:
I am trying to build a web application using html and javascript from which I want to login to geoserver hosted at: https://weaveprod.ucdp.utah.edu/geoserver/web/ I have written the following code for authentication using javascript.But it is always giving login failed.Kindly help me out.               HTML code:                    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />     <title>Untitled Document</title>          <link rel="stylesheet" type="text/css" href="login.css">
     <form>         Username:         <input type="text" id="username" name="username">         Password:         <input type="password" id="password" name="password">
             <input type="checkbox" id="checkbox" name="checkbox">Keep me logged             in             <input type="submit" value="Login" onclick="login(option)">
              </form>
                           Javascript code:          function login (options) {     // url del servlet del geoserver     var url = options.server;     // parametros para el login     params = "username=" + options["user"] + "&password=" + options["password"];          var contentType = "application/x-www-form-urlencoded";     //se inicializa la petición ajax     var ajax = $.ajax({         data : params,         type : "GET",         crossDomain: true,         dataType: "jsonp",         contentType : contentType,         url : url     });     // se ejecuta cuando la peticion finaliza     ajax.done(function() {              if ($.cookie("JSESSIONID") != null && options && options.success) {             options.success();         }     });     // si ocurrio un error al realizar la peticion     ajax.fail(function() {         if (options && options.failure) {             options.failure();         }     });     // se ejecuta siempre al final de la petición, sin importar que esta     // haya fallado     ajax.always(function() {         if (options && options.always) {             options.always();         }     });     };      NOTE: Earlier in the javascript code, I used a POST request.But it is still giving a login failed message.The webpage and the server are in different domains and I do not have access to the source code of the server. Thanks, Sayan


title: How to reproject an EVF?

text: I have an .evf file and it is not opening in QGIS. How can we reproject the vector layer (specifically if we want to convert it from lat/lon to UTM projection)?


title: How to simply label raster cells?

text: I think the title says everything but: I want to generate a small sample raster (eg. euclidean distance or cost distance) to explain the method of cost-distance analysis, and I'd like to label the pixels with the corresponding cost-value. I already tried to lay a fish-net over the raster in arcgis, but I didn't find a good way to copy the values into the attribute table. I'd appreciate every hint! :)


title: Why the modifications I do in QGIS doesn't appear when I save an OSM file I created?

text: Through QGIS I loaded an OSM file I created and I'm trying to visualize what happens with the file in terms of codes after some modifications like adding points, lines or polygons. But the problem is that when I save the file with the OpenStreetMaps plugin the modifications don't appear as I was expecting... Here you can see the file I created and if you want, feel free to do the same test adding points or whatever, saving and then see if something changes in the new saved file.               <?xml version="1.0" encoding="UTF-8"?>     <osm version="0.6" generator="OpenStreetMap server">     <bounds minlat="52.2615973629" minlon="10.5262132592" maxlat="52.2653605629" maxlon="10.5313793802"/>     <node id="100" lat="52.264859" lon="10.528042" version="1" user="ribe_lu" visible="true" timestamp="2012-09-10T13:47:15Z"/>     <node id="101" lat="52.264859" lon="10.529796" version="1" user="ribe_lu" visible="true" timestamp="2012-09-10T13:47:15Z"/>     <node id="102" lat="52.262619" lon="10.528042" version="1" user="ribe_lu" visible="true" timestamp="2012-09-10T13:47:15Z"/>     <node id="103" lat="52.262619" lon="10.529796" version="1" user="ribe_lu" visible="true" timestamp="2012-09-10T13:47:15Z">         <tag k="made_by" v="Lucas"/>     </node>     <node id="104" lat="52.264186" lon="10.529078" version="1" user="ribe_lu" visible="true" timestamp="2012-09-10T14:20:17Z"/>     <node id="105" lat="52.26326" lon="10.529078" version="1" user="ribe_lu" visible="true" timestamp="2012-09-10T14:20:17Z"/>     <node id="106" lat="52.26326" lon="10.528801" version="1" user="ribe_lu" visible="true" timestamp="2012-09-10T14:20:17Z">         <tag k="made_by" v="Lucas"/>     </node>     <way id="25" visible="true" timestamp="2011-01-09T14:05:48Z" version="1" user="ribe_lu">         <nd ref="100"/>         <nd ref="101"/>         <nd ref="103"/>         <nd ref="102"/>         <nd ref="100"/>         <tag k="area" v="yes"/>         <tag k="Schloss" v="Arkaden"/>         <tag k="Postal_Code" v="38100"/>     </way>     <way id="26" visible="true" timestamp="2011-01-09T14:05:48Z" version="1" user="ribe_lu">         <nd ref="104"/>         <nd ref="105"/>         <nd ref="106"/>         <nd ref="104"/>         <tag k="Test_by" v="Lucas"/>     </way>     <relation id="75" visible="true" timestamp="2012-09-10T14:50:10Z" version="1" user="ribe_lu">         <member type="way" ref="25" role="outer"/>         <member type="way" ref="26" role="inner"/>         <tag k="type" v="multipolygon"/>     </relation>     </osm>      Thank you very much, Best regards, Lucas.


title: How to exclude outliers from raster using ArcGIS for Desktop?

text: I calculated terrain curvature from a DEM. The DEM I use is interpolated from point elevation data. Generally, the DEM seems to be good, but some artifacts from interpolation exist. Those artifacts cause extreme values for curvature. Around 90% of the curvature values are between -25 and 25, but the other 10% are between -500 and 500. I can use histogram equalize to properly visualize the curvature raster. For further analysis however I'd like to replace this extreme values. So my question is: How can I exclude this extreme values from the raster? I imagine something like this would work (pseudo-code):               m = median("curvature")     IF("curvature.value") > m * 10 THEN ("curvature.value") = m * 10     IF("curvature.value") < m * 10 THEN ("curvature.value") = m * 10      What tool would I use for this? I tried Raster Calculator, but I couldn't get the syntax right. Software: ArcGIS 10 with spatial and 3D analyst


title: Which is the correct projection in ArcGIS for EPSG:3857/900913?

text: I want to publish a mapservice in ArcGIS server, which will be overlayed on Google Maps or Bing Maps. For this, the Map Service needs to be in EPSG:3857/900913. The Google Maps projection is supposed to assume a Spherical Earth Projection. So I went to > Coordinate Systems\Projected Coordinate Systems\World (Sphere-based) There seems to be a _Mercator (sphere).prj_ , but it does not seem to have the proper Projection parameters. I then went to > Coordinate Systems\Projected Coordinate Systems\World & there are 3 WGS84 based mercator projections: ![Coordinate Sytem](http://i.stack.imgur.com/XKfa5.png) I see that in a lot of ESRI literature, they use WGS 1984 Web Mercator (Auxiliary Sphere). Is that the Google Maps Projection? What is the _WGS 1984 Web Mercator.prj_ then? (As a bonus question, this image is from ArcGIS desktop 10.0, but does the answer again change at ArcGIS 10.1?)


title: GPS and GIS for real time monitoring of an object

text: How is it possible to store coordinate of object in server in real time using feasible hardware ( can be used by student) GPS into database and retrived in real time??   1. What are the names of equipment ?   2. **Any link to good resource?** The purpose of this experiment is to provide real time location of object through internet(web).


title: How to show an infowindow on click of an Image symbol marked at specific coordinates in ArcGIS map.

text: markers(Pushpins) are loaded with values from XML and is displayed on the map.How to enable click events on the markers to display InfoWindow.                       <esri:Map x:Name="MyMap" Style="{StaticResource MyMapStyle}" >             <esri:Map.Layers>                 <esri:ArcGISTiledMapServiceLayer ID="StreetMapLayer"                                              Url="http://server.arcgisonline.com/ArcGIS/rest/services/ESRI_StreetMap_World_2D/MapServer" />             <esri:GraphicsLayer ID="GraphicLayerMarkers">                      </esri:GraphicsLayer>             </esri:Map.Layers>         </esri:Map>


title: Error point sampling tool

text: I'm usig QGIS 1.8.0 and I did some interpolations. Now I want to do some sampling with the point sampling tool but I keep getting this error (I attach a picture of the error message). Does anybody knows why? Any help would be really appreciated. ![enter image description here](http://i.stack.imgur.com/GLZbm.png)


title: Can i install opengeo suite in VPS windows?

text: I want to quickly set up a gis server which i can use via internet .i guess opengeo suite is best for newcomer and also it has community edition . but if you think there is option let me know. I need to work with WMS, WFS, WFS-t Services . I am developer i am also looking to build some custom web app in future . so is VPS (windows) good ? shall i concerned with windows server version ? eg 2008 , 2012 ? I think company is good option https://www.arvixe.com/ ? has anyone done this before ? I just want to make sure that my money is not wasted .. thanks for your feedback. Is ubuntu or centOS good option if i don't know about linux ?


title: Rotating Point feature problem

text: I am trying to create new symbols for our Stormwater Stream Restoration dept.. Utilizing some of the ESRI came symbols, I have adjusted their size and color and added them to the map. When I click the Rotate tool, it shows the feature rotating in the preview as I turn it, but when I let go of the mouse, the screen blinks, but the symbol never rotated. Also, when I'm in Symbology/Advance/Rotation, the only options in the drop down are: *None, Random, and expression. Btw... I don't know how to code or write any scripts yet. Thanks


title: How to split a line at point positions?

text: I am a rooky with QGIS and i am stuck from a basic problems. I have a polyline and i have some points (2) on the polyline. I would like to know if there is a tool or script that i could use to split my original polyline into 3 polylines? (Starting point to point A, Point A to B and B to ending point of the polyline)? I am looking for an automatic process rather splitting the line manually. Thanks a lot for your help Baptwo


title: Is there a way to properly have GDAL raise exceptions in Python?

text: I am working with GDAL in Python and have wanted to get GDAL to raise exceptions so that I can do something useful with them. For instance if a file cannot be opened with gdal.Open() I would like to throw an exception instead of gdal returning None. I would just like a way to capture the errors instead of them being dumped to sys.stdout. Currently I have tried using the `gdal.UseExceptions()` function but came across an issue. When `gdal.UseExceptions()` is set it seems to hide certain errors and not raise an exception on them. For instance `ERROR 5: Access window out of range in RasterIO() Requested (1,15) of size 25x3 on raster of 26x17` disappears and no exception is raised. This seems pretty dangerous. Any ideas are appreciated, thanks.


title: Floodplain generation given river section polylines, flood level elevation points and a DEM of the land surface

text: I am looking to generate flood plain polygons to represent 20 and 100 year flood zones, in the most accurate and efficient way possible, given the following datasets: 1) River cross section polylines 2) River flow polylines (center lines between opposing shorelines) 3) 20 year and 100 year elevation values for the vertices of each cross section 4) A digital elevation model (raster with an elevation value for each pixel) of the land surface Here is my proposed methodology so far: 1) Form polygons delimiting each section, the river center line and inland areas well beyond the flood plain 2) For each polygon formed in the previous step, extract the intersecting raster surface (Extract by Mask) 3) Isolate the pixels of the DEM that are less than or equal to the flood elevation (20 and 100 year) value for the section using the Reclassify tool 4) Perform a raster (reclassified raster from previous step) to vector (polygons) to form floodplains between each section. 5) Merge the polygons together to form a continuous floodplain I am looking for insight and suggestions to produce the most accurate results possible given the available data. I have access to ArcGIS and the Spatial Analyst and 3D Analyst extensions. Cheers!


title: ArcObjects 10.1

text: just updated to 10.1 and I'm having issues moving from 10.1. I get the errors with my ADF reference when trying to create custom tool bars examples:               'BaseToolbar' is ambiguous in the namespace 'ESRI.ArcGIS.ADF.BaseClasses'.        Any suggestions on how to convert? I have the following references imported:               Imports Microsoft.VisualBasic     Imports System     Imports Microsoft.Office.Interop     Imports ESRI.ArcGIS.Geodatabase     Imports ESRI.ArcGIS.ADF     Imports ESRI.ArcGIS.esriSystem     Imports System.Windows.Forms     Imports ESRI.ArcGIS.Carto     Imports ESRI.ArcGIS.ADF.Connection.Local


title: OpenLayers ModifyFeature and styling

text: In OpenLayers, I have a Vector layer with a StyleMap and SelectFeature and ModifyFeature controls, like this:               vector_layer = new OpenLayers.Layer.Vector('My Vector Layer', {extractAttributes: true,                 styleMap: new OpenLayers.StyleMap({'default':{                     strokeColor: "${customColor}",                     strokeOpacity: 1,                     strokeWidth: 2,                     fillColor: "${customColor}",                     fillOpacity: 0.5                 }}) });      map.addLayer(vector_layer);               selectControl = new OpenLayers.Control.SelectFeature(vector_layer,                     {onSelect: onFeatureSelect, onUnselect: onFeatureUnselect});     map.addControl(selectControl);     selectControl.activate();     modify_control = new OpenLayers.Control.ModifyFeature(vector_layer);     map.addControl(modify_control);      When I tick a checkbox on the page, the ModifyFeature control is activated:               modify_control.activate();     modify_control.mode |= OpenLayers.Control.ModifyFeature.ROTATE;     modify_control.mode |= OpenLayers.Control.ModifyFeature.DRAG;      The problem is this: when I remove the StyleMap, ModifyFeature shows the handles to drag and rotate the feature. _With_ the StyleMap, they are not visible. Any suggestions how to fix this? Thanks!


title: pyqgis: How to get features of selected geometries

text: I wonder if anybody could provide an example how to get the attributes of selected Features? I tried the folowing Code in the python Console - but I'm stuck at the point where I'd like to get the attributes:               qgis.utils.iface     canvas = qgis.utils.iface.mapCanvas()     cLayer = canvas.currentLayer()     selectList = []     if cLayer:         count = cLayer.selectedFeatureCount()         print count         selectedList = layer.selectedFeaturesIds()         for f in selectedList:             # This is where I'm stuck             # As I don't know how to get the Attributes of the features      I would be really glad if anybody could help. Regards


title: Best Landsat-5 TM band combination for detecting fire scars

text: I am looking for the best band combination with Landsat 5 imagery to detect fire scars.


title: Access to the nearest line by clicking on OpenLayers

text: i need to access to the nearest line by clicking on my Layer how can i do it? for example take a look at this photo.


title: Matching two point clouds without loss (nearest neighbor, but more complicated)

text: I have two point clouds, one with demographic and one with electoral data. Points in both clouds represent the same empirical entity (say, a village, or an urban neighbourhood), but are not in the exact same lat/long (different datasources, GPS inaccuracies, etc). Moreover, one village can have several points in each cloud. I am looking for a cluster or nearest neighbor algorithm which allows me to merge both point clouds. I have thought about two strategies: a) do a nearest neighbor matching, starting from the first point cloud and adding a variable to each point that directs me to the respective closest point in the second point cloud. The problem is that no method for nearest neighbor matching which I found (GRASS, SAGA, R) ensures that all points in the second cloud are picked at least once as this "respective closest point" for the first cloud. A neverending circle which could probably approximately be optimized - but I have no idea how. Any suggestions welcome. b) the second strategy sees it as a clustering problem, in which I attempt to recreate the underlying empirical entity "village", which might contain multiple electoral and demographic points, but must contain at least one of each type. In other words, I would need a bottom-up hierarchical cluster algorithm which groups points together based on their spatial proximity, but simultaneously (!) ensures that each cluster contains at least one point from the electoral point cloud and one point from the demographic point cloud. Again: any suggestions welcome. c) there are probably other solutions that I am unaware of - I am truly a beginner here... It would be marvellous if somebody could point me in the right direction. Thank you in advance!


title: Is the ArcMap style database format available somewhere?

text: The *.style files are MS Access databases where each style type is a table. The table has one binary column where the actual style is written. Is this binary format available? I would like to import some custom color ramps, but I can't see a way to do it... and I hoped in the style database!


title: Algorithm to find geohashes which satisfy query condition

text: I have an application where I need to query geohashes by specifying the bounding rectangle. Since geohashes map 2D to 1D is there a method by which I may find the prefixes of all geohashes which satisfy the query condition (bounding rectangle). I tried to google but I did not find any algorithm. It will be great if there is some library which does that for me efficiently in C/c++. Even if there is no library then too I can can work with an algorithmic solution.


title: How do I make KML overlays appear in OpenLayers at all zoom levels?

text:
I need to overlay KML polygons on basic map in OpenLayers. I need that KML is visible on every zoom level. I did some coding, but my KML is not showing up. What am I doing wrong? Thanks.               <!DOCTYPE html>                     <title>OpenLayers Tutorial - Basic Map Setup</title>
OpenLayers Testing


title: VBscript Regular Expression in ArcMap Field Calculator

text: I'm looking for a regex in VBscript that I can use on a ArcMap table. I know what it would look like in python: `pat =r'[0-9]{3,4}'` This filters all numeric characters from the string which have the length 3 to 4. So from AO123 you'd get 123 and from VO2-324C, 324. What would be the equivalent in VBScript? How do I implement it with the RexExp Object? ArcGIS has a Field calculator that let's you execute VBScript code. But even the most simple regex (see screenshot) fails with the error message object needed: 're'. ![enter image description here](http://i.stack.imgur.com/Ay8bJ.jpg)


title: Check if field value for selected features is null

text: My feature class contains two fields: an identifying field `U_ID`, and a descriptive text field `DESC`. This field currently has `<Null>` for all features. I would like to select all the features with the same ID and populate this field with sequential values (using code I have previously used here), only if there are `<Null>` values in `DESC`. My model receives the feature class as input into a `Feature Selection Iterator`, grouping by `U_ID`. I then want to check for null values in the `DESC` field for those `Selected Features`. If a `<Null>` is found in `DESC` for any of the `Selected Features`, return true and populate the field with values. If there are values for all the features, the condition is false. This is what the model looks like: ![enter image description here](http://i.stack.imgur.com/qmhes.png) For example: the iterator selects `U_ID` = AA01. There are ten features in the feature class with that ID, so `Selected Features` contains those 10 features. The `Calculate Value` tool contains the following:               checkIfNull("%DESC%")      * * *               def checkIfNull(f):       return "true" if not f else "false"      Seeing as `DESC` is `<Null>` for those ten features (or even just one of the features), the precondition is true, and `Calculate Field DESC` should be executed. However, the precondition returns false and the model ends. Changing the code to               return "true" if f is None else "false"      yields the same result. I don't think I'm accessing the values in `DESC` correctly, because my test for `<Null>` seems to be right, judging by the answers to other questions here.


title: Mapping x, y pixel coordinates to latitude and longitudes on the WGS84 datum

text: _Apologies in advance for the long question, I am fairly confused at this point and I figure if I show my thought process any flaws early on in the piece can be more easily pointed out._ I have a grid of 1000 by 1000 (origin is top left), on which I'm drawing a variety of polygons. I'm trying to convert the x, y locations of each of the vertices of these polygons programmatically to longitude and latitude for use in geojson. Effectively, I want to map my grid around the globe. I have run into a few problems, I suspect largely because I'm fairly inexperienced with GIS. I spent some time reading up on what I could, this explanation of map projections and coordinate systems was particularly helpful, and after reading it and other sources, I had the following plan. With the bounds of -180 degrees to 180 degrees longitude, and 84 to -80 latitude, I should be able to convert arbitrary geographic coordinates to geocentric (x, y) and then scale those values to exist between 0 and 1000, and reverse the process to go the other way. For the tools to do this, I decided upon pyproj, a wrapper around proj4 as **1)** We use python   **2)** It seems fairly specific to the task at hand which is usually a good sign From my reading, the bounds of longitude and latitude (in the utm projection, which is what I plan to use) are -180 to 180, and 84 to -80               from pyproj import proj          p = Proj(proj='utm', zone='56H', ellps='WGS84')          def x_and_y(lon, lat):         x, y = p(lon, lat, errcheck=True)         return (x, y)      I thought i'd use the following code to work out the x, y values of the lat and long boundaries, and use when scaling from 0-1000 to whatever the min and max x and ys would be. This is where I run into a brick wall, effectively. **1)** Using 0, 0 as lat and long gives me `RuntimeError: latitude or longitude exceeded limits`, is this expected?   **2)** Actually, a lot of values give me errors, i.e., if I do the following:               longs = range(-180, 181, 12)     lats = range(84, -81, -4)          #print longs     #print lats          for lon in longs:         for lat in lats:             print '%s  %s' % (lon, lat)             print '%s  %s' % x_and_y(lon, lat)      I get errors from and including -108 degrees longitude through to and including 60 degrees longitude. At this point I realised my understanding must be fundamentally flawed and decided to turn to you. I guess I have the following questions: **1)** Am I travelling down the right path here? If my goal is to map (x, y) coordinates onto a global map, is this the way to go?   **1a)** If I'm not, recommendations?   **2)** Clearly my understanding is flawed, can anyone see and explain where and why?   **3)** Are there any tutorials or code samples around for this kind of thing? The pyproj documentation is fairly limited. Finally, if anyone has recommended reading for this kind of work in general, I'd love a pointer or two on where to start. Thanks for your time! PS: If this is better suited for stack overflow, please let me know.


title: Add Spatial Database Connection fails?

text: I created a Database Server (SQL Express 2008.) I created a "New GeoDatabase..." I double-click on "Add Spatial Database Connection." I fill out the dialog and change the transactional version. I press the "Test Connection" button. The "Test Connection" message box says, "Connection succeeded." I press the "OK" button twice and... NO connection is created. I cannot create a new spatial database connection. (I searched for this on the ESRI forum and, of course, the question is there but as usual there is no answer.) Has anyone else come across this problem?


title: Collapse layers in ArcMap Table of Contents using ArcObjects

text: Can someone offer up a Christmas present here? I'm stuck on how to collapse all layers in my ArcMap projects (9.3). I've tried a few bits of sample code, but they all seem to fail for some reason or other. Some of the layers might be grouped, and some are imagery or WMS layers if that makes a difference. Here is some code I have tinkered with, but does not work:               Dim pMxdoc As IMxDocument     Set pMxdoc = ThisDocument     Dim pFeatLayer As IFeatureLayer     Dim pEditLayers As IEditLayers     Set pEditLayers = pEditor     Dim pMap As IMap     Set pMap = pMxdoc.FocusMap     Dim pLayer As ILayer          'Collapse any expanded grouplayers     Dim pGroupLayer As IGroupLayer     For i = 0 To pMxdoc.FocusMap.LayerCount - 1     Set pLayer = pMxdoc.FocusMap.Layer(i)     'Check if the reached layer is a group     If TypeOf pLayer Is IGroupLayer Then         Set pGroupLayer = pLayer         if pGroupLayer.Expanded = true Then             pGroupLayer.Expanded = false         End if     End If     Next i     pMxdoc.UpdateContents     pMxdoc.ActiveView.Refresh      ![enter image description here](http://i.stack.imgur.com/rbiys.jpg)


title: Is there a way to get the map document name/path from a published map service?

text: I have several map services published to ArcGIS Server 10.0 SP2 and I have layer metadata files exported to HTML that I want to link to the layers in the map service. Currently the metadata HTML files are named based on the data source (eg 'Zip_Code_Poly.htm') but I need to correlate the layer name/ID from the map service to that file. I was thinking that the easiest way to do it would be to run an arcpy script to get all the layer data sources. I've run into an issue though that the client (this is a web app) needs to look up these HTML files by service name and layer name/ID but without hard coding a mapping from service -> MXD I can't figure out what map document to use in arcpy. Is there a way using arcpy or something else that I can look at a published map service and get the map document it was created from?


How do I make geoserver update layers and stores? I've updated the tables in PostGIS used by a data store and layer, but I cannot seem to get it to update. I have the new data running, but I had to create a new Store and Layer just to get it to work. Is there an easier way to refresh stores and/or layers? Thanks


title: Geoserver image generation error

text: I'm using Geoserver + PostGIS and I got a image generation error when i test my layer in url: http://localhost:8089/geoserver/pmc/wms?service=WMS&version=1.1.0&request=GetMap&layers=pmc:utbs&styles=&bbox=272206.156,7448903.5,303656.5,7484798.0&width=448&height=512&srs=EPSG:4326&format=image/jpeg this is the image error: ![enter image description here](http://i.stack.imgur.com/kam9A.jpg) http://img695.imageshack.us/img695/9963/layerimage.jpg.


title: How do I find out the resolution in (centimeters per pixel) of a TIF?

text: I tried using GDAL and got:               # gdalinfo nicholas_map_2.tif      Driver: GTiff/GeoTIFF     Files: nicholas_map_2.tif            nicholas_map_2.tfw     Size is 28799, 27875     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (-122.325695999999994,37.876846000000000)     Pixel Size = (0.000000441883000,-0.000000350276000)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       COMPRESSION=LZW       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  (-122.3256960,  37.8768460) (122d19'32.51"W, 37d52'36.65"N)     Lower Left  (-122.3256960,  37.8670821) (122d19'32.51"W, 37d52' 1.50"N)     Upper Right (-122.3129702,  37.8768460) (122d18'46.69"W, 37d52'36.65"N)     Lower Right (-122.3129702,  37.8670821) (122d18'46.69"W, 37d52' 1.50"N)     Center      (-122.3193331,  37.8719640) (122d19' 9.60"W, 37d52'19.07"N)     Band 1 Block=28799x32 Type=Byte, ColorInterp=Red       Mask Flags: PER_DATASET ALPHA      Band 2 Block=28799x32 Type=Byte, ColorInterp=Green       Mask Flags: PER_DATASET ALPHA      Band 3 Block=28799x32 Type=Byte, ColorInterp=Blue       Mask Flags: PER_DATASET ALPHA      Band 4 Block=28799x32 Type=Byte, ColorInterp=Alpha      I assume pixel size is what I'm looking for? But is that in degrees?


title: How to change table name of command raster2pgsql

text: I am completely newbie with GIS and PostGIS as well I'm trying to load raster into postgreSQL with raster2pgsql. It runs well. However, the table name of the raster is generated automatically. For example, the raster named: landsat.tiff, then the table name is landsat. My question is: Can we change that name of table? and how we can do it? Thanks in advance! **EDIT:** Here's what I typed in Windows Command line:               \raster2pgsql -s 4326 *.TIFF -t 100x100 > demo.sql      and insert into database               \psql -U postgres -d demo -f demo.sql


title: Is it possible to create MBtiles from a TIF?

text: I have a geo-referenced tif (it's really big around 1GB) and I'd like to convert it to mbtiles so I can serve them and view the map. What is the best tool to do this conversion?


title: How to use InsertCursor to add points to a point feature class?

text: I am struggling with a very simple task of inserting "SHAPE@XY" into a point feature class. Following is what I have:               xy = [(361087.15436001675, 198747.07296730898), (361530.56426259346, 198976.62248019245),           (361692.3015462744, 197815.78617074457), (361845.0959187896, 198297.07296730898),           (362135.3345548634, 197552.7531621556)]          fc = "C:/temp/xy.shp"     cursor = arcpy.da.InsertCursor(fc, ["SHAPE@XY"])     for row in xy:         cursor.insertRow(row)          del cursor      it throws back error:               TypeError: sequence size must match size of the row      I used field calculator to add x and y field into the same shapefile form where I retrieved these points, and it works fine to add above set of co- ordinates to a multiline feature class, but above set of lines fail. Any suggestions?


title: Videos in CartoDB Infowindows

text: I need to know how to embed a video (YouTube or Vimeo) into an infowindow in CartoDB in the same way that an image appears as an image header. Potentially more complicated, when there is no video link I'd like an image to appear instead. Can anyone assist me (as a non-programmer) with the code for this and WHERE to insert it in the HTML? I have viewed the "How do you embed iFrame content into a CartoDB infoWindow?" question, but I cannot make that work for either YouTube or Vimeo videos. Also what I am asking for is for it to be a part of the "image header" in CartoDB which is not what that question answers. Thank you!


title: Can I upload a video file to embed in a custom infowindow in Cartodb?

text: I want to embed video clips into custom infowindows in Cartodb, and I know how to do this via a youtube link, but is there a way I can embed a video from a file I have saved on my computer? Or does the video need to come from a website?


title: Need help installing QGIS 2.0 Need GDAL then SQLite3 3.7.16

text: OSX 10.7.5 I need help installing QGIS 2.0 It says I then need GDAL. When I try to install GDAL it says SQLite3 3.7.16 or newer framework is required. I can't find SQLite3 3.7.16.


title: Field Calculator extract three last digits from another field

text: How everyone! I have two fields the on is the cadastral number which is for example 051020213001 and I am trying to fill the other field (label) with the last 3 digits from the cadastral number, in this case 001. Is there any way I can do that automatically and not manually because I have 2800 polygons in which I have to do that. I am using QGIS 2.2.0


title: Reproduce raster with the same extent, and the same number of columns and rows

text: I have two raster A and B. A and B have the same grid size but different extent. How can I make A has the same number of columns and rows, as well as the spatial extent of B?


title: How can a leaflet SVG overlay be placed outside of the leaflet-overlay-pane?

text: I have SVG on top of GeoJSON in my map. How can a leaflet SVG overlay be placed outside of the leaflet-overlay-pane while keeping GeoJSON in place in the leaflet-overlay-pane? (For example, how could I put just the SVG into the leaflet-shadow-pane and then turn pointer-events off on the leaflet-shadow- pane... thus making it possible to interact with the GeoJSON?) My apprehension is this: Leaflet map layers are stored in panes: http://leafletjs.com/reference.html#map-panes However, it is not clear to me just exactly how these panes are ordered and what each of them contains. according the Leaflet's css (http://cdn.leafletjs.com/leaflet-0.7.2/leaflet.css), the order of these panes makes more sense:   * **.leaflet-map-pane { z-index: <"?">; }** //Pane that contains all other map panes. (presumably all the below map panes)     1. **.leaflet-tile-pane { z-index: 2; }** //Pane for tile layers.     2. **.leaflet-objects-pane { z-index: 3; }** //Pane that contains all the panes except tile pane. (presumably all the below panes)       1. **.leaflet-overlay-pane { z-index: 4; }** //Pane for overlays like polylines and polygons.         * **.leaflet-overlay-pane svg { z-index: <"?">; }**         * **.leaflet-overlay-pane <"?"> { z-index: <"?">; }**       2. **.leaflet-shadow-pane { z-index: 5; }** //Pane for overlay shadows (e.g. marker shadows).       3. **.leaflet-marker-pane { z-index: 6; }** //Pane for marker icons.       4. **.leaflet-popup-pane { z-index: 7; }** //Pane for popups. I am specifically interested in **.leaflet-overlay-pane svg** and any subsequent panes that reside inside **.leaflet-overlay-pane** (i.e. **GeoJSON pane** ). Thanks


title: Re-registering a GEOTIFF in qgis

text: I'm working on some images taken from a UAV over the course of a year. They are OSGB 1936 projection as is the whole project. Annoyingly they don't quite line up so I am trying to a few of the images using the georeferencer tool. ![The control points](http://i.stack.imgur.com/0tQ1G.png) ![The settings](http://i.stack.imgur.com/x0Vm0.png) The problem is, is that when I run it, the resulting alignment is miles away. ![Where it ends up](http://i.stack.imgur.com/SlrvT.png) You can see from the georeferencing screen that it doesn't need to move much. Every tutorial I've seen assumes that it's a completely new image, i.e. jpg or something, that needs referencing rather than a previously referenced image. Am I going about this the wrong way?


title: What is the best method to snap features to road network?

text: I have a road network with some of my tract centroids and facilities that fall off of it, however I would like to execute my analyses so that all features are part of the road network? My facilities are hospitals and incidents are census tract centroids. I am executing a closest facility analysis for drive time and distance. I would like my output to adjust these features to the network with an output in an attribute table which shows the distance the feature was adjusted to the network. Much like the near analysis, but actually adjusting features so that they fall within the network.


title: How do I calculate the new shapefile areas after clipping?

text: I'm working with the WWF ecoregions, which use the GCS_WGS_1984 geographic coordinate system. I have clipped the ecoregions with a species range (consisting of two polygons), but the new attribute table doesn't tell me the new 'clipped' ecoregion areas. From reading around it seems I need a projected coordinate system to calculate these new areas. However the attribute table of the ecoregions contain area (km2) as well as shape area. Would these have been calculated using a projected coordinate system? If so, is there a way of finding out which one they used? Complete beginner here, thanks for any help in advance.


title: How do I open images georeferenced in ArcMap in QGIS? (.aux.xml)

text: We have several ArcMap georeferenced images (raster datasets) that we are attempting to open in QGIS without any luck. These are not inside of a geodatabase and are not Arcinfo Binary Grid format or Arcinfo Binary Coverages. QGIS supports GDAL rasters and the ArcMap created .aux.xml file contains the georeference points, but QGIS apparently does not recognize this file as it places my Texas image in a very small area of the Atlantic ocean off the coast of Africa (with a google maps basemap added). (THe only other files in the raster dataset are .tfwx and .ovr) Is there something special I need to do within QGIS to apply the coordinates in the aux.xml created by ArcMap, or something else I'm missing? Thanks in advance, Justin Here are the contents of one such aux.xml file:               <?xml version="1.0"?>          -<PAMDataset>          <Metadata/>               -<Metadata domain="Esri">          <MDI key="PyramidResamplingType">NEAREST</MDI>          </Metadata>               -<Metadata domain="xml:ESRI" format="xml">               -<GeodataXform xmlns:typens="http://www.esri.com/schemas/ArcGIS/10.1" xmlns:xs="http://www.w3.org/2001/XMLSchema" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:type="typens:PolynomialXform">          <PolynomialOrder>1</PolynomialOrder>               -<SpatialReference xsi:type="typens:ProjectedCoordinateSystem">          <WKT>PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0],AUTHORITY["EPSG",3857]]</WKT>          <XOrigin>-22041259.177065704</XOrigin>          <YOrigin>-33265069.154224794</YOrigin>          <XYScale>135385247.70505452</XYScale>          <ZOrigin>-100000</ZOrigin>          <ZScale>10000</ZScale>          <MOrigin>-100000</MOrigin>          <MScale>10000</MScale>          <XYTolerance>0.001</XYTolerance>          <ZTolerance>0.001</ZTolerance>          <MTolerance>0.001</MTolerance>          <HighPrecision>true</HighPrecision>          <WKID>102100</WKID>          <LatestWKID>3857</LatestWKID>          </SpatialReference>               -<SourceGCPs xsi:type="typens:ArrayOfDouble">          <Double>5.994281757832141</Double>          <Double>5.6742238617780458</Double>          <Double>2.0592992159490677</Double>          <Double>4.5229583592636118</Double>          <Double>2.014907522695573</Double>          <Double>8.8668953443025202</Double>          <Double>2.8272115205686532</Double>          <Double>7.7051733542912189</Double>          </SourceGCPs>               -<TargetGCPs xsi:type="typens:ArrayOfDouble">          <Double>-10752195.345278211</Double>          <Double>3857531.2694506948</Double>          <Double>-10756947.906142347</Double>          <Double>3856205.4206410511</Double>          <Double>-10756935.616529588</Double>          <Double>3861470.3939956115</Double>          <Double>-10755980.997953685</Double>          <Double>3860058.8390891682</Double>          </TargetGCPs>          <Name/>          </GeodataXform>          </Metadata>          </PAMDataset>


title: How to 'reverse clip' (erase) in R

text: A reverse clip saves only the part of your spatial object that is _outside_ the bounds of another object, as opposed to a regular clip which saves the parts that are _inside_ the other object. Here's how to do it in ArcMap. The question is: how to do this in R? Reproducible example (on Linux machines):               system("wget 'https://github.com/Robinlovelace/Creating-maps-in-R/archive/master.zip' -P /tmp/")     unzip("/tmp/master.zip", exdir = "/tmp/master")     uk <- readOGR("/tmp/master/Creating-maps-in-R-master/data/", "ukbord")     lnd <- readOGR("/tmp/master/Creating-maps-in-R-master/data/", "LondonBoroughs")     plot(uk)     plot(lnd, add = T, col = "black")      What I want here to do is to save all of the UK _except_ for London. Visually, I want the black shape in the resulting image to be a hole. ![enter image description here](http://i.stack.imgur.com/I2Peb.png)


title: Show layers TOC in QGIS 2.0.1

text: I somehow closed the layer TOC in QGIS 2.0.1. How do I get it back? Feeling a bit lost here.


title: How to filter wrong GPS readings?

text: I have absolutely no experience with GIS, it just so happens that one of the hobby projects I am doing at the local hackerspace includes tracking a model boat with a GPS... I hope this is the right place to ask! [If not, I would appreciate suggestions on where to turn for help]. Here's my problem: today we went testing the new model boat, but when I processed the log, I noticed that for some strange reason, some of the GPS readings are totally off (see picture below: we never actually sailed on the left of the trees). ![Google Earth screenshot](http://i.stack.imgur.com/V9KXd.png)   1. I was wandering if there is a **standard way to filter a GPS log to spot and eliminate such noisy signals**. The only obvious way I could think of was to tell the filter the a maximum speed for the boat that was tracked, and simply exclude any signal that would imply the boat to have travelled fastest than that maximum speed. Is there a better way?   2. Since our boat is a model one and most of the times it travels at speeds below the knot, I was wondering if there is a **clever way to "rectify" points that are probably wrong**. An example: suppose the boat is travelling east to west and the track looks something like: ...⋅... it is quite possible that the central dot is a glitch in the GPS, but rather than dropping the log signal altogether, I would like to "correct" it [each log record contains information about the electronics too, and this data is always correct]. Again: is there a standard way to "smooth" a GPS track? I thought some sort of solution could pass through the standard deviation, but I am clueless on how to proceed. Additional information: we are using python for processing the log, but while code snippets are always appreciated, I'd be more than pleased with a general explanation of the algorithm that I would implement myself. Many thanks in advance for your time and expertise.


title: How to select and copy a group of polygons in QGIS?

text: I'm new to QGIS. I opened a few shapefiles that cover Europe. Now I want to select a smaller area and if possible copy the polygons in this area so I can work with a smaller set. How should I do this in QGIS? I started to choose "toogle to editing", then I choosed the "node tool" and started to make rectangle. But then I get a popup (see screenshot below) telling me to set a "tolerance". I looked in (as suggested in the popup) "Settings" > "Project properties" > "General", but couldn't find and any tolerance value here. From the User Guide (version 1.7.0) in chapter 4.5.1 it was suggested that I should change the tolerance in "Settings" > "Options" > "Digitizing" and I changed from 0.00000 to 1.00000. It was also suggested to edit "Settings" > "Snapping options" so I changed the tolerance to the same value there. Now I tried with the "Node tool" again, but I got the same popup and couldn't select any polygons. How can I select an area of polygons in QGIS? that I later can copy to a new project in QGIS to get a smaller set to work with. ![enter image description here](http://i.stack.imgur.com/lVva9.png)


title: Connect to ArcSDE geodatabase over LAN?

text: I am running ArcGIS Server 10, with ArcSDE 10. I have a geodatabase that I would like my colleague to connect to over the LAN at our office. How do I go about doing this? I would prefer a list of simple steps than a link to the ESRI Web Help URL. Thanks in advance! :)


title: How to create a map of only one state in QGIS?

text: This is a very simple question but a Google search did not provide an answer. Basically suppose you have a map of the U.S. with states and county layers, as well as other layers like roads, etc. I want to work with a map of only one state, with all layers (e.g. counties, roads, etc) in that state. I remember in ArcGIS it had something to do with taking an intersection, not sure how it works in QGIS.


title: How to call a DLL file from VBA?

text: I have a DLL file that I have loaded and registered it in ArcGIS as a command icon. This icon launches a window where the user should do the following selections: an input file, a field name, select between 2 option buttons and an output file name. However, I need to run this DLL file via a VBA and ArcObjects code so as the user will not see the window of the DLL tool and the above selections will be given automatically by the code within a loop so as to take various output files. Is that possible? Please help. Thanks Demetris Eventually I found a function that may be is approprite to call the DLL file. I declare it below and I try to call it via the Calldll() procedure. However, I receive the following message "Run-time error 13- Type mismatch". I think that the arguments have the appropriate data type. Could you please have a look in the code below to find what is going on? Thanks again Demetris                Public Declare Function CreateParcelsShape Lib "C:\Program    Files\ArcGIS\CreatePoly.dll" (Centroids As IFeatureClass, Block As  IPolygon, pFDS As IFeatureDataset)          Public Sub CallDll()          Dim pMxDoc As IMxDocument     Set pMxDoc = ThisDocument          'Get the active map (data frame)     Dim pMap As IMap     Set pMap = pMxDoc.FocusMap          'Get the layers     Dim pLayers As IEnumLayer     Set pLayers = pMap.Layers          Dim pLayer1 As ILayer     Set pLayer1 = pLayers.Next               Do Until pLayer1 Is Nothing            If pLayer1.Name = "Block12" Then         Exit Do            End If          Set pLayer1 = pLayers.Next     Loop               Dim pBlocks As IFeatureLayer     Set pBlocks = pLayer1          Dim pBlocksFC As IFeatureClass     Set pBlocksFC = pBlocks.FeatureClass               Dim pBlocksCursor As IFeatureCursor     Set pBlocksCursor = pBlocksFC.Search(Nothing, False)          Dim pBlockFeature As IFeature     Set pBlockFeature = pBlocksCursor.NextFeature          'Get the shape of the polygon (Block)     Dim polygon As IPolygon     Set polygon = pBlockFeature.Shape          Dim pLayers2 As IEnumLayer     Set pLayers2 = pMap.Layers          Dim pLayer2 As ILayer     Set pLayer2 = pLayers2.Next          Do Until pLayer2 Is Nothing               If pLayer2.Name = "CentroidsBlock12" Then     Exit Do          End If          Set pLayer2 = pLayers2.Next     Loop          Dim pCentroids As IFeatureLayer     Set pCentroids = pLayer2          Get the dataset     Dim pFDataSet As IFeatureDataset     Dim pFWorkspace As IFeatureWorkspace     Dim pAWFactory As IWorkspaceFactory          Set pAWFactory = New AccessWorkspaceFactory     Set pFWorkspace = pAWFactory.OpenFromFile("C:\LACONISS\GAPopulation.mdb", 0)     Set pFDataSet = pFWorkspace.OpenFeatureDataset("Polygons")          'Call the DLL function     CreateParcelsShape pCentroids, polygon, pFDataSet          End Sub


title: Detecting map changes using Openlayers

text: Is it possible to delete a vector layer when a user pans or zooms a map? I saw in the documentation that there is a function called onMapResize. Can this also be used to detect a pan change? I want to delete a point if the map is panned or zoomed and tried using it in the following way (which does not work):               map.OpenLayers.vectorLayer.onMapResize(console.log('map was moved'))     //use removeAllFeatures() to remove vector points


title: User Input in Openlayers Problem

text:
I thought you could give me some advice regarding my problem. I have a working WFS Layer and I want to be able to type in some text which then should be used as the filter value. Afterwards the map should be refreshed, however only with the features selected through my input text. For Example I want to see Restaurants, so I type it in and it refreshes my map. This is my code so far. In advance - thanks for your help.                                 <title>OpenLayers Demo</title>
         <form name="form" id="form">             <input type="text"  name="Search" id="searchText1" value="" />         <input type="submit"  />           </form>


title: How to import shp2pgsql?

text: I want to import shape file into PostGIS. There is 1 database `postgres` in it but no table added yet. I am using this command:               shp2pgsql -I -s 4326 districts.shp districts | psql -d postgres -U postgres       and it's giving me some error: > Unable to convert data value to UTF-8 (iconv reports "Invalid or incomplete > multibyte or wide character"). Current encoding is "UTF-8". Try "LATIN1" > (Western European), or one of the values described at > postgresql.org/docs/current/static/multibyte.html Is my command wrong ??


title: How to offset points along line in QGIS?

text: I used the "Profile from Line" plugin to create evenly spaced points along a line in QGIS, but I'd like to offset them to either side of the line. I spose I could draw two lines parallel to the existing one and create evenly spaced points along each of them then delete the new lines, but I'm curious if there is a different way to offset points the way I've shown in my MS Paint illustration. ![offset points along line](http://i.stack.imgur.com/C4hoh.jpg)


title: How to overlay a Shapefile on a WMS layer using Geotools?

text: I am trying to overlay a shp file on wms layer using the example http://docs.geotools.org/latest/userguide/tutorial/raster/image.html. Each time I want to read the format of the file to get a reader I get this error :               Exception in thread "main" java.lang.UnsupportedOperationException: Trying to get a reader from an unknown format.     at org.geotools.coverage.grid.io.UnknownFormat.getReader(UnknownFormat.java:62)     at com.qedrix.map.maplotr.Demo1.displayLayers(Demo1.java:127)     at com.qedrix.map.maplotr.Demo1.main(Demo1.java:260)      My code follows:               public class Demo1 {          private AbstractGridCoverage2DReader reader = null;          private StyleFactory sf = CommonFactoryFinder.getStyleFactory();     private FilterFactory2 ff = CommonFactoryFinder.getFilterFactory2();          /**      * This method examines the names of the sample dimensions in the provided      * coverage looking for "red...", "green..." and "blue..." (case insensitive      * match). If these names are not found it uses bands 1, 2, and 3 for the      * red, green and blue channels. It then sets up a raster symbolizer and      * returns this wrapped in a Style.      *       * @return a new Style object containing a raster symbolizer set up for RGB      *         image      */     private Style createRGBStyle() {         GridCoverage2D cov = null;         try {             cov = reader.read(null);         } catch (IOException giveUp) {             throw new RuntimeException(giveUp);         }         // We need at least three bands to create an RGB style         int numBands = cov.getNumSampleDimensions();         if (numBands < 3) {             return null;         }         // Get the names of the bands         String[] sampleDimensionNames = new String[numBands];         for (int i = 0; i < numBands; i++) {             GridSampleDimension dim = cov.getSampleDimension(i);             sampleDimensionNames[i] = dim.getDescription().toString();         }         final int RED = 0, GREEN = 1, BLUE = 2;         int[] channelNum = { -1, -1, -1 };         // We examine the band names looking for "red...", "green...",         // "blue...".         // Note that the channel numbers we record are indexed from 1, not 0.         for (int i = 0; i < numBands; i++) {             String name = sampleDimensionNames[i].toLowerCase();             if (name != null) {                 if (name.matches("red.*")) {                     channelNum[RED] = i + 1;                 } else if (name.matches("green.*")) {                     channelNum[GREEN] = i + 1;                 } else if (name.matches("blue.*")) {                     channelNum[BLUE] = i + 1;                 }             }         }         // If we didn't find named bands "red...", "green...", "blue..."         // we fall back to using the first three bands in order         if (channelNum[RED] < 0 || channelNum[GREEN] < 0 || channelNum[BLUE] < 0) {             channelNum[RED] = 1;             channelNum[GREEN] = 2;             channelNum[BLUE] = 3;         }         // Now we create a RasterSymbolizer using the selected channels         SelectedChannelType[] sct = new SelectedChannelType[cov.getNumSampleDimensions()];         ContrastEnhancement ce = sf.contrastEnhancement(ff.literal(1.0), ContrastMethod.NORMALIZE);         for (int i = 0; i < 3; i++) {             sct[i] = sf.createSelectedChannelType(String.valueOf(channelNum[i]), ce);         }         RasterSymbolizer sym = sf.getDefaultRasterSymbolizer();         ChannelSelection sel = sf.channelSelection(sct[RED], sct[GREEN], sct[BLUE]);         sym.setChannelSelection(sel);              return SLD.wrapSymbolizers(sym);     }          public void displayLayers() {              File rasterFile = fetchWmsImage();              AbstractGridFormat format = GridFormatFinder.findFormat(rasterFile);              this.reader = format.getReader(rasterFile);              // Initially display the raster in greyscale using the         // data from the first image band         Style rasterStyle = createRGBStyle();              // Create a basic style with yellow lines and no fill         Style shpStyle = SLD.createPointStyle("point", Color.YELLOW, Color.GRAY, 0.0f, 1.5f);              MapContent map = new MapContent();         map.setTitle("ImageLab");              MapViewport vp = new MapViewport();              org.geotools.map.Layer rasterLayer = new GridReaderLayer(reader, rasterStyle);         map.addLayer(rasterLayer);              saveImage(map, "final.jpeg", 583);          }          public File fetchWmsImage() {              URL url = null;         try {             url = new URL("http://184.106.187.247:8080/geoserver/rg/wms?version=1.1.0");         } catch (MalformedURLException e) {             // will not happen         }              WebMapServer wms = null;         try {             wms = new WebMapServer(url);                  WMSCapabilities capabilities = wms.getCapabilities();             Layer[] layers = WMSUtils.getNamedLayers(capabilities);                  GetMapRequest request = wms.createGetMapRequest();                  request.setFormat("image/png");             request.setDimensions("583", "420");             request.setTransparent(true);             request.setSRS("EPSG:900913");             request.setBBox("-13019428.542822,3922163.1648461,-13013051.407366,3929863.8567165");             //request.setProperty("isBaseLayer", "false");             //request.setProperty("opacity", ".2");                  for (Layer layer : WMSUtils.getNamedLayers(capabilities)) {                 if (layer.getName().equals("rg:parcels"))                     request.addLayer(layer);             }                  System.out.println(request.getFinalURL());                  GetMapResponse response = (GetMapResponse) wms.issueRequest(request);             BufferedImage image = ImageIO.read(response.getInputStream());                  File rasterFile = new File("C:\\Users\\samabhik\\Workspace\\MAP\\data\\out.png");             ImageIO.write(image, "png", rasterFile);                  try{                 Thread.sleep(2000);             }catch(InterruptedException e){}                  return rasterFile;              } catch (ServiceException e) {             // TODO Auto-generated catch block             e.printStackTrace();         } catch (IOException e) {             // TODO Auto-generated catch block             e.printStackTrace();         } finally {              }              return null;     }          public void saveImage(final MapContent map, final String file, final int imageWidth) {              GTRenderer renderer = new StreamingRenderer();         renderer.setMapContent(map);              Rectangle imageBounds = null;         ReferencedEnvelope mapBounds = null;         try {             mapBounds = map.getMaxBounds();             double heightToWidth = mapBounds.getSpan(1) / mapBounds.getSpan(0);             imageBounds = new Rectangle(0, 0, imageWidth, (int) Math.round(imageWidth * heightToWidth));              } catch (Exception e) {             // failed to access map layers             throw new RuntimeException(e);         }              BufferedImage image = new BufferedImage(imageBounds.width, imageBounds.height, BufferedImage.TYPE_INT_RGB);              Graphics2D gr = image.createGraphics();         gr.setPaint(Color.WHITE);         gr.fill(imageBounds);              try {             renderer.paint(gr, imageBounds, mapBounds);             File fileToSave = new File(file);             ImageIO.write(image, "jpeg", fileToSave);              } catch (IOException e) {             throw new RuntimeException(e);         }     }          public static void main(String[] args) {              Demo1 demo = new Demo1();         demo.displayLayers();     }      } My pom dependency looks like this:               <dependency>             <groupId>org.geotools</groupId>             <artifactId>gt-shapefile</artifactId>             <version>${geotools.version}</version>         </dependency>         <dependency>             <groupId>org.geotools</groupId>             <artifactId>gt-swing</artifactId>             <version>${geotools.version}</version>         </dependency>         <dependency>             <groupId>org.geotools</groupId>             <artifactId>gt-epsg-hsql</artifactId>             <version>${geotools.version}</version>         </dependency>         <dependency>             <groupId>org.geotools</groupId>             <artifactId>gt-geotiff</artifactId>             <version>${geotools.version}</version>         </dependency>         <dependency>             <groupId>org.geotools</groupId>             <artifactId>gt-image</artifactId>             <version>${geotools.version}</version>         </dependency>         <dependency>             <groupId>org.geotools</groupId>             <artifactId>gt-wms</artifactId>             <version>${geotools.version}</version>         </dependency>         <dependency>             <groupId>org.geotools</groupId>             <artifactId>gt-coverage</artifactId>             <version>${geotools.version}</version>         </dependency>         <dependency>             <groupId>org.geotools</groupId>             <artifactId>gt-imageio-ext-gdal</artifactId>             <version>${geotools.version}</version>         </dependency>      I am using geotools version 9-SNAPSHOT I have added the native libraries from gdal in java classpath in eclipse. I am totally stuck and don't know what's wrong. I even tested the geoserver WMS URL and it works:               http://184.106.187.247:8080/geoserver/rg/wms?SERVICE=WMS&LAYERS=rg:parcels&ISBASELAYER=false&FORMAT=image/png&OPACITY=.2&HEIGHT=420&TRANSPARENT=TRUE&REQUEST=GetMap&BBOX=-13019428.542822,3922163.1648461,-13013051.407366,3929863.8567165&WIDTH=583&STYLES=&SRS=EPSG:900913&VERSION=1.1.1      Just ran gdalinfo.exe to read the image content and it follows:               Driver: PNG/Portable Network Graphics     Files: ..\..\Workspace\MAP\data\out2.png     Size is 583, 420     Coordinate System is `'     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  (    0.0,    0.0)     Lower Left  (    0.0,  420.0)     Upper Right (  583.0,    0.0)     Lower Right (  583.0,  420.0)     Center      (  291.5,  210.0)     Band 1 Block=583x1 Type=Byte, ColorInterp=Red       Mask Flags: PER_DATASET ALPHA     Band 2 Block=583x1 Type=Byte, ColorInterp=Green       Mask Flags: PER_DATASET ALPHA     Band 3 Block=583x1 Type=Byte, ColorInterp=Blue       Mask Flags: PER_DATASET ALPHA     Band 4 Block=583x1 Type=Byte, ColorInterp=Alpha


title: programmatically create and publish geotiff layer in geoserver

text: We need to automatize process of creating datastore and publishing geotiff layer. Also, before publishing we need to change Declared SRS and to save Native Bounding Box to our database. Geotiff file will be uploaded to the file system by our customer. Yes, I`ve seen Create a Layer in GeoServer using REST, also I see that there is publishGeoTIFF(...) methods in Geoserver-manager in GeoServerRESTPublisher.java. Another solution that I see - is to make all the operations with **bot written with selenium or htmlUnit**. Could you please suggest us, what approach is more elegant and suitable for us. Is there any actions in our workflow that we will not be able to do with Geoserver-manager. P.S.: we also have seen at gsrcj project, but it can`t create programmatically geotiff layer. Thank you in advance, Dmitry.


title: Instructions for building a Flex App using open source PostGIS DB and Geoserver

text: Anyone have a link for a good guide on how to build a Flex app to work with Open source REST services? I am familiar with how to do it in ESRI and Arc Server Manager, but don't want to pay for a ArcServer license. Thanks! -Adam


title: Imposm full database schema import template

text: I noticed that the default imposm template for importing OSM data into postgis lacks some types of features (powerlines for example). I was searching the web to find a complete schema import template (in the form of python classes, like the default template) but couldn't find any, does it exist? It would be much easier to have a complete template, then remove what you dont need. Thanks


title: Importing XML file in ArcGIS failed

text: I have created a model in Enterprise Architect (in ArcGIS technology). Inside I have several point, polygon, line feature classes and a couple of domains which correspond to them. From the level of EA I export my model to ArcGIS XML Workspace. Then when I try to import it in ArcGIS 10 I get an error as below: ![enter image description here](http://i.stack.imgur.com/IDq1N.jpg)


title: Python module to delete SHP features (without Desktop GIS installed)

text: I have a piece of software (not Arc) that runs nightly on a non-production machine which updates shapefile features from an external database. Periodically I would like to delete all of the features in the shapefile (not the file itself, which must remain) and let the software "rebuild" the shapefile from scratch. I would like to automate this process. I don't have any GIS software currently installed on that machine. I was hoping that I could script a routine in Python that would automatically delete the features, much like Arc's Delete Features geoprocessing tool. Are there any Python modules that would allow me to do this? Preferrably open- source? I looked at Shapely and PyShp but didn't see anything that would allow me to delete the features in mass or that matched a WHERE clause. They can write features and analyze them, but haven't seen DELETE FEATURES functions anywhere. I must certainly be missing something simple...


title: UniqueValueRenderer and ClassBreaksRenderer configuration

text: So I can create one of these renderers based on parameters I let the user set. Some parameter combinations either create invalid conditions or no conditions at all - missing data, too many breaks or too few, and so on. This obviously unacceptable. I've considered doing what the desktop product does and display the list of breaks and let the user fiddle with the parameters until they do something useful. Has anyone done that? Is there a sample somewhere on how to do that? If not can someone offer advice on either doing that or suggest another way to ensure that a renderer "works?"


title: IUniqueValueRenderer on Polyline and Points made my Features all invisible

text: I am currently using IUnique value renderer like this:                foreach (string nodeType in nodes)           yield return new KeyValuePair<string, ISymbol>(nodeType, CreateNewSymbol() );               public static ISimpleLineSymbol CreateSimpleLineSymbol(IRgbColor rgbColor, Double width, esriSimpleLineStyle style)     {         SimpleLineSymbolClass simpleLineSymbol = new SimpleLineSymbolClass();              simpleLineSymbol.Color = rgbColor;         simpleLineSymbol.Style = style;         simpleLineSymbol.Width = width;              return simpleLineSymbol;     }      finally this is called during creation of the layer:                 /// <summary>             /// Set's the unique symbol types on the feature layer             /// </summary>             /// <param name="featureLayer"></param>             /// <param name="layerType"></param>             public static void SetFeatureLayerSymbols(IFeatureLayer featureLayer, IEnumerable<KeyValuePair<string, ISymbol>> layerSymbols)             {                 IFields fields = featureLayer.FeatureClass.Fields;                      //Get indexes of certain fields                 int id = fields.FindField("TYPE");                      IUniqueValueRenderer uniqueValueRenderer = new UniqueValueRenderer();                      //The fields which are used to map to the symbol                 uniqueValueRenderer.FieldCount = 1;                 uniqueValueRenderer.set_Field(0, "TYPE");                      foreach (KeyValuePair<string, ISymbol> values in layerSymbols)                     uniqueValueRenderer.AddValue(values.Key, null, values.Value);                      if (uniqueValueRenderer.ValueCount > 0)                     (featureLayer as IGeoFeatureLayer).Renderer = uniqueValueRenderer as IFeatureRenderer;             }      I can see the index in the TOC showing the proper unique values and the colors, but on the map it is invisible, I CAN HOWEVER HIGHLIGHT IT. Thanks and Regards,


title: Suggestions for Geocoding Messy, Inconsistent, Single-Field Addresses?

text: I have nine spreadsheets containing information about a total of 33,401 unique events in the greater Chicago area. I have been asked to geocode all of these, if possible, and I am certainly no stranger to geocoding. However, the location information for each of them is the worst I have ever seen, written in a single field, with no particular convention. I have no ZIP codes or city names, but I do have county names in almost every case. When street names are included, they are often missing their suffix ("Ave", "St", "Rd"). State and US highways are both frequently coded as the indiscriminate "Rte/Rt". The majority of locations have been written as intersections (often with extra irrelevant information), such as:               SB Pulaski & 162nd St.     I-55 @ Rt.30     Devon and Cicero (Il 50) NW corner TS     NB Rt.41 @ Half Day Rd. Exit Ramp.      In the case of Interstates, these "intersections" often do not actually exist -- it's just referencing a street the Interstate passes over. A fair number have (relatively) proper addresses:               1800 s Wolf rd. south of Oakton, north of Touhy.     1010 S. Rt. 14 - in front of Thunderbird Country C     Grayslake Maintenance Yard, 217 N. Baron, Grayslak      Some are more vague, but still generally well-specified:               South bound Busse rd south of Oakton and Higgins     EB Elgin-O'Hare W of Rohlwing Rd     NB IL-59,  1 mile north of IL-132      And some are almost certainly impossible to locate without additional context:               EB Elgin  O'Hare expresway     Prairie View Rest Area     Comm Center/Stevenson Yard      My question is, given the wide range of formats and anything-goes approach for specifying addresses in these datasets, are there any suggested methods for parsing at least some of this stuff into a reasonably clean set of geocodable addresses? I've so far been stumped and have been going through the painful process of making sense of individual records in Google Maps. I want to trim off as much manual work as possible, as I'd prefer not to spend the next three years working on this.


title: Proxy configuration problem with Glassfish and GeoServer

text: I have Geoserver running on 8080 and Glassfish on 9999.So for wfs-t i had to use some proxy script.But at the response part of firebug i got following code: http://localhost:9999/schemas/ows/1.0.0/owsExceptionReport.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ows="http://www.opengis.net/ows"> Could not determine geoserver request from http request org.geoserver.platform.AdvancedDispatchFilter$AdvancedDispatchHttpRequest@eba9b0


title: How to add a background map to a java project?

text: I'm new at GIS. I am trying to make a java application that can display a map with certain points. I have achieved this using the geotools library. However, I want to put the map from OpenStreetMap or Google Streets in the background also. I already know that PostGIS can be used with GeoTools. However I don't know where to start. How is the best way to start once I have installed PostGIS? Hello, I have also come across KML. I'm sorry I'm a beginner. What KML parsing mean, and how can I use this to show a map in my Java Application?


title: How can I add dynamic text to a layout from an attribute table?

text: I was trying to figure out how to use dynamic text from an attribute table in a collection of Data Driven Pages in ArcGIS 10. I've looked in the help files, but can't find anything on dynamic text from attribute tables.


title: How to convert a vector boundary to a raster map?

text: I need to output a PNG file of a raster map overlayed with boundaries present in a vector map. The problem is, the r.out.png script only exports a raster map. My boundaries are in a vector format. When I tried to convert the vector boundaries into a raster, in order for me to combine them, the converted raster from vector contained the filled vector boundaries. I only need the lines to be present in the raster map. Is there a way to do it? I cannot use the d.rast and d.vect since I am calling the scripts in a shell script so there would not be any GRASS windows open.


title: Use buffer tool to draw different radii?

text: ![enter image description here](http://i.stack.imgur.com/XnP63.png) I am new to Arcmap-Arceditor 10 and could really use your help. I have been mapping ethanol plants (which utilize corn to produce ethanol) and corn production amounts in the counties of the USA. My goal is to draw buffer zones (circles) around the ethanol plants (dots). However, the tricky part is that the circles must reflect what region an ethanol plant buys corn from in order to sustain its ethanol production. I know how much corn each plant needs and I know the average corn production/square mile in each county. How do I get Arcmap to calculate the radiuses of the buffer zones it needs to draw corn from? This gets especially tricky when I have an ethanol plant lying on the border of two counties and each county has a different corn production amount/square mile. Do you know if it's possible to do this in arcmap? If so, how should I proceed? Thank you so much! In the image attached, the dots are ethanol plants. The red regions show the county is self-sufficient in producing enough corn for its ethanol plants. The orange regions show that corn must be imported from other counties to support its ethanol plants.


title: Create metadata for multiple datasets simultaneously in ArcGIS 10

text: I am trying to find a way to create metadata through Arc Catalog for multiple datasets simultaneously. I'm working within a file geodatabase, which composes of multiple feature classes and rasters. The data has common theme so I'd like to create template of metadata for one item and then populate other datasets within the same metadata information. I'm planning to use ISO standard metadata template. I've done some research and it appears that there are three options (but neither of them too successful so far): - use 'Metadata Importer' tool within Conversion Toolbox (however I will be able to copy metadata on one-to-one basis only) - X-tools pro have 'Batch Metadata Editing' option available but after setting the parameters within ArcMap I can't see these canges through Arc Catalog - third option is adding 'Batch Importer' (http://edndoc.esri.com/arcobjects/9.0/Samples/Metadata/Importers/Batch_Importer/Batch_Importer.htm). I managed to add it to ArcCatalog but I'm not able to run it (I get error message 'ICommand_OnClick() ). I'm working with ArcGIS 10. Any help/suggestions would be much appreciated. Thanks in advance,


title: Can JQuery Mobile UI provide a user-created attribute table for a user-created point feature?

text: I'm looking into using the JQuery Mobile UI, and I would like the user of the app to be able to create their own point on the map and add a few short attributes. Any idea if this is possible? Thanks in advance.


title: Cannot establish connection to port 6080

text: From what I understand, the default port for ArcGIS for Server is port 6080. I have installed ArcGIS for Server 10.2 and now I cannot connect to the Server Manager. I am using Windows 8. My firewall is completely turned off. I have restarted the service about 100 times and I have reinstalled ArcGIS for Server about 10 times at this point. The port is not shown when I do a netstat -an in command prompt as well. Any help would be appreciated. Thank you in advance.               URL I used: http://localhost:6080/arcgis/manager      Picture of "error": ![In Chrome](http://i.stack.imgur.com/jiPAY.png) ![In FireFox](http://i.stack.imgur.com/EDWBq.png) ![In IE](http://i.stack.imgur.com/6Ol1K.png)


title: Why does QGIS 2.0 print composer not render the legend names?

text: I made a map using the print composer in QGIS 2.0. However, in the legend it does not display the name of all the items listed. I have updated the legend items, removed the items and added them again but does not work. Also, if I add an image to the map it is not displayed (item send to front does not make a difference). Am I missing a setting or is it a bug? ![enter image description here](http://i.stack.imgur.com/dAVob.gif)


title: What's the best way to add Tilemaps to QGIS without having to use Webmercator?

text: I'm looking for a good way to add tilemaps in qgis that are either in a custom projection or that I can reproject in qgis. I have used the openlayers plugin to add basemap layers, but it forces me to use webmercator. I have also spent some time playing with adding layers as a gdal wms as described here http://www.gdal.org/frmt_wms.html, but layers added in this manner are very slow even in their native projection. My end goal would be to have ortho photos in a fast format in my local projection, and to be able to make some custom basemaps.


title: How can I maintain MapInfo GRC file attributes in QGIS/GDAL?

text: I am trying to manage a classified MapInfo grid file in QGIS - I'm wondering if anyone has any experience with maintaining the attributes of the classified grid? QGIS seems to assign them just a number (Grid value 51, 52, 53 instead of site1234_sector1, site1234_sector3). It does seem to recognize that the different attribute values exist, and lists them in the metadata, but I don't see any simple way to pull that table linking the new grid values to the old grid attributes. Any advice would be appreciated. Thanks! EDIT: Current system that works involves copying the attribute list out of the metadata tab and pasting it into excel. This works, but clearly isn't ideal.