title: Why is my add-in toolbar greyed out?

text: I installed an addIn toolbar for the U.S. Geological Survey's National Cooperative Geologic Mapping Program and I can see it, but all the tools are greyed out. Any idea why I can't see it? I am inserting a jpg of what I see. Thanks!!! ![greyedouttoolbar](http://i.stack.imgur.com/bMMgB.jpg)


title: Vertical Scroll CartoDB not working

text:
I am currently developing a site using CartoDB. I have a fixed div and a table located inside the div. The div is only 300px in height, but the information goes way beyond this height. The issue I am having is when you attempt to scroll while hovering over the box, nothing happens. However, if you click in the div and use the up and down arrows on your keyboard it "scrolls" just fine. Why am I having this issue? The div is fixed directly above the bottom layer which is a map. CSS:               .holder{       position:fixed;       top:80px;       left:8px;       width:226px;       height:900px;       margin:0;       padding:7px 0 0 0;     }     .users-wrapper{       width:226px;       height:300px;       overflow-x:hidden;     }     .users{       position:relative;       display:block;       margin:0;       margin-top:30px;       padding:5px;       background:#FFF;       height:auto;       min-height:5px;       max-height:300px;       width:226px;       max-width:226px;       overflow-x:hidden!important;     }      HTML:

         //Table information located here


title: What parts of the ArcGIS 9.3 Geoprocessor Programming Model diagram are not relevant for ArcGIS 10?

text: The geoprocessor programming model was an excellent tool for python scripting using ArcGIS 9.3. There is still no equivalent diagram for ArcGIS 10 - although there is a request on the ArcGIS ideas site. A previous question on this site asked if the 9.3 diagram is good to use for ArcGIS 10. The answer was that some of it is relevant and it is best not to use it. However, I really like the diagram and would like to reference it when scripting ArcGIS 10. Therefore my specific question is: What parts of the ArcGIS 9.3 diagram are not relevant to ArcGIS 10? Correct me if I'm wrong, but from what I understand, if an equivalent diagram was made for ArcGIS 10, there would be new functions added to the diagram; the diagram would only get bigger.


title: PyQGIS: How to send signal if layer is currently in edit mode

text: I'm looking for a way to send a signal if a layer is in editing mode. I tried editingStarted(), but it wasn't working. I'd like to do something along the lines of the what I have below, but can't find it in the API (where I'm looking: http://qgis.org/api/classQgsVectorLayer.html#a476eaf6ec8f2a45e285e963bcacedcc2)               qgis.utils.iface.setActiveLayer(myLayer)     #capslock is in place of the method i'm searching     if myLayer.ISINEDITMODE() == False:         qgis.utils.iface.actionToggleEditing().trigger()


title: How to interactively get the coordinates of a point in QGIS

text: I have a shapefile with a single layer of type point. It is in the google CRS (EPSG:900913). When I use the 'Identify Feature' tool it gives me lots of information about the point clicked on, but I cant see the coordinates. Is there a way to find out the lat/lng of a point interactively?


title: How can I add a new field in a raster attribute table and loop through it?

text: I am working with an existing raster attribute table. I wonder how can I add a new field to it and then loop through the existing records to populate this new field? Of course I have done this many times for a feature class and a database table but I do not know how to do it with a raster table. I work with VBA. Any code please? Thanks.


I'm having some doubts concerning the layer mechanism. Creating a new layer (ie a new connection between a vector and a DB table), I've imported a vector (rivers) through v.in.ogr, then I've cloned the original table:               db.copy from_table=rivers to_table=rivers_clone      I've created a new connection:               v.db.connect map=rivers table=rivers_clone layer=2      At this point I supposed that the last command had created all the needed structures and fields for the new layer (2), but if I print the categories:               v.category input=rivers layer=2 option=print      no categories are displayed (while they are listed if I choose layer 1). **Question** : I thought that a new connection used the same vector categories, but it seems I was wrong. How to make the layer 2 use the same vector categories from layer 1?


title: How to calculate the flood volume using a LiDAR DEM and a HEC-RAS floodplain polygon

text: I’ve got a vexing GIS analysis problem which is conceptually fairly straightforward, but in reality more challenging to figure out…. Basically, I need to calculate the VOLUME of floodwater within a 100 year floodplai, using ArcGIS 10.0 (all extensions available). My input datasets are:   * LiDAR DEM (1m pixel or 3m pixel) of bare ground (raster)   * 100 year floodplain (output of USACE HEC-RAS analysis), Polygon Shapefile, but attribute table has no useful content. So far I have overlaid the floodplain polygons over the LiDAR raster, [using Extract by Mask tool from Spatial Analyst], so that I have the ground elevations are along the outside edge of the floodplain (think of the top of the bathtub). Key concept is that the top of the floodplain gradually falls in elevation as you move from the headwaters down through the watershed. Then, (this seems to be the impossible part) I need to figure out how I can take the elevations along the outside edge pixels of the raster floodplain (output of previous Extract by Mask) and extend them to the equivalent elevation on the other side of the floodplain, and “raise” all the pixels in between the 2 outside edges of the floodplain to the “top level” of this 100-year floodplain. Once I have this “top elevation” of the floodplain water, I should be able to use the CUT AND FILL tools in ArcGIS to simply subtract the LiDAR ground elevations from the “top floodplain elevation” and compute a volume between these two rasters. I’ve been researching this for days and have yet to find a solution. Perhaps I’m thinking of this completely wrong? I would have expected to find others who have dealt with this problem, but while there are similar problems, none appear to offer a solution. Many thanks in advance!


title: Lightning Data Analysis

text: Hi I am new to model builder in GIS. I have two shapefiles, one is a polygon shapefile of all lightning caused fire in Alaska sine 1939 (~1500 polygons) and the second is a point shapefile with all recorded lightning strike since 1986 in Alaska (~10000 points). What I want to do is create a data/attribute table that outputs the following for varying time scale (1 day through 16 weeks): 1)If the fire had a lightning strike within it on the timescale specified and latitude/longitude of strike 2)If not what is the nearest strike and its distance (point (lightning) to centroid (fire polygon)or point to nearest edge) 3)Polarity of strike 4)Multiplicity of strike The attribute table I have for the lightning data has all strike lat/long, polarity, multiplicity. This is a flow chart of the analysis I want to do: select lightning and fire that falls in specified time (I have dates of both strikes and fire) --> determine if a strike was in the fire --> if yes then put the specified value from the lightning table into the new table along with fire ID --> If no find the nearest strike to the fire and its distance to the fire --> add specified values for this strike to same table What would the basic outline for a model that performs this flow look like? Are there any alternate methods to achieve this with such a large data set? I know this is a lot but any help would be greatly appreciated. Thanks


title: territory analysis

text: I want to determine the territories in which customers to a bank has the shortest path. It can be done with desktop software packages or programming if it is needed. Does anybody has idea on how to start doing it? I want to produce a map like this http://www.esri.com/industries/banking/business/~/media/Images/Content/industries/banking/graphics/DesireLineOverlay- lg.jpg and this map only the territories are important for me.


title: Python GDAL: find out spikes and pits in DEM

text: I'm working with on a DEM (Geotiff) in Python GDAL, and I'd like to find out the spikes or pits in my terrain.   I have no idea how to proceed, does anyone have good algorithm to help me out?   Thanks,   eo


title: Failure to display R script result plots in QGIS

text: I'm trying to run R scripts (3.0.2) from QGIS (2.0.1) on Windows 7 x64. There seems to be a known bug whereby the output results do not display in the QGIS GUI, but rather just an html path ending in *.html.png. I've managed to work around this by having my script end in               jpeg("file")     plot(x,y)     dev.off()      which outputs my plot to a jpeg file which I can open in windows. Does anyone know if there's a more appropriate work around which can take place within QGIS?


title: Topography: How to get the x,y,z coordinates of highest or lowest point/pixel of a raster GIS?

text: Given a .tiff topographic raster image with mountains and deep sea, I am looking for the x,y,z specifics of the highest or lowest point. The x,y,z values being:   * altitude (z-coordinate)   * latitude (y-coordinate)   * longitude (x-coordinate) While `gdalinfo raster.tif` would give a summary of the bounding box and the highest / lowest elevation values (see my previous question here, and note comments below), I would like to find the elevation, latitude, longitude of the highest and/or lowest pixels of my raster. Would be better if these values are got as isolated output integer. I think raster processing iterating on all pixels of the raster and keeping values of the z-highest point may be a good way to go.


title: Google maps api + Open layers

text: I need to open google maps in Latitude: -30.0330600 Longitude: -51.2300000. We report that the object map or call the api?                   <!DOCTYPE html>     <html lang='en'>     <head>      <meta charset='utf-8' />      <title>My OpenLayers Map</title>      <script type='text/javascript' src='OpenLayers.js'></script>           <!-- Chamada a Api Google Maps -->      <script src="http://maps.google.com/maps/api/js?sensor=false&v=3.2"></script>           <script type='text/javascript'>           //Variável map      var map;           //Função a ser chamada para criar o mapa na leitura da página      function init() {           // A classe Map espera dois argumentos "map_element" que o ID do elemento HTLML geralmente uma DIV e opções de mapa {key:value}      map = new OpenLayers.Map('map_element',{         maxExtent: new OpenLayers.Bounds(             -128 * 156543.0339,             -128 * 156543.0339,             128 * 156543.0339,             128 * 156543.0339),         maxResolution: 156543.0339,         units: 'm',         projection: new OpenLayers.Projection('EPSG:900913'),         displayProjection: new OpenLayers.Projection("EPSG:4326"),      });           //Função de inicialização - google.maps.MapTypeId.ROADMAP (the default)      var google_streets = new OpenLayers.Layer.Google(         "Google Streets",         {numZoomLevels: 20}         );                 //Adiciona as camadas ao mapa     map.addLayers([google_streets]);                //Camada de controle que vai mostrar as camadas no mapa      map.addControl(new OpenLayers.Control.LayerSwitcher({}));                center = new google.maps.LatLong(-30.0330600,-51.230000);           if(!map.getCenter(center)){      map.zoomToMaxExtent();           }      }           </script>     </head>          <body onload='init();'> <!-- Chama a função js init() -->           <!--Elemento HTML onde o mapa é exibido -->      <div id='map_element' style='width: 50px; height: 500px;'>      </div>          </body>     </html>


title: Strange structures appear inside dry river in Egypt

text: Browsing Google maps site reveal some strange structures inside one of the dry rivers in Egypt near Quesser , what is really interesting is that those structures do not appear on Bing maps . Searching the internet does not shed light on the subject, also there are no photos present on the site. What could those structures be and why they appear on Google maps only? **Edit** If those structures have been built somewhere in between the two years window. What could those structures be? Is it possible to build structure like this in a given time span at a place like this or are there more structures like this in the world? If this structure serves as a bridge, I don't see any roads leading to it.


title: Embed map in a website

text: I am using arcgis server 10.0. I published a mxd file using arcgis manager resources. I see the web url `http://cosmic- dellt410:8399/arcgis/services/DD/MapServer` but it says get method not allowed. I noted that the wsdl get generated but not the rest api. Any idea why? How do I embed the map on my website?


title: Datums in compound CRSs

text: Is there some specific constraint regarding _vertical_ and _geodetic_ **datums** within a compound Coordinate Reference System (CRS)? E.g. if the ellipsoid model should be the same or not, or if the vertical datum can be barometric, etc. As also _engineering_ CRSs can be compound with geodetic 2D CRSs, then I don't believe there are such sort of constraints. Thanks for any comment.


title: From OSM to shape. Parsing problem

text: I'm working at OpenStreetMap project for android at the moment. I want to retrieve OSM data from my app(nodes, markers and paths) and convert it to shape file (.shp, .dbf and .shx). Later, I'll load it into ArcGIS editor. I've been investigating about it, and I have no clear idea if It's possible a direct conversion between this formats, or I have to step first on KML file and parse it in shape file. I actually don't know how to retrieve the OSM data to save a osm file for being treated. I've been looking for a library or something like that but I have no lucky. Can anyone give me a helping hand? Thanks in advance!


title: How to start developing with Jeppesen C-MAP?

text: I've just downloaded all the contents from the FTP of Jeppesen for their product called **C-MAP CM93**. I've been reading the documentation contained in the SDK folders, but I can't figure out how to start developing in C# or Java with it. Does anybody know what can I do? Where I can find a DLL? Any information on the Web? Thank you!


title: Can't change Data Frame display units in ArcMap 10

text: Last week, I went to View -> Data Frame Properties... -> General tab to change the data frame map display units from Decimal Degrees to Degrees Decimal Minutes. (the Map Units were still in Decimal Degrees, in fact it is greyed out and can't be changed) Today, I wanted to change the Display units back to Decimal Degrees, but after I click Apply and OK, the map status bar still shows coordinates in Degrees Decimal Minutes as I move my cursor around the map. I went back to Data Frame Properties and the Display units have reverted back to Degrees Decimal Minutes! Why is this happening, and what should I do? Thanks.


title: Relative hyperlink paths in ArcMap 10 Map Document Properties?

text: Is there a way to use relative paths (instead of absolute paths) in the "Hyperlink base" field in Map Document Properties for an ArcMap 10 map document? (I have "Store relative pathnames to data sources" checked, but somehow I still can't use relative pathnames in Hyperlink base)


title: How to read and write to .txt in layer properties, action tab?

text: I want to write a certain attribute to a .txt file using the action tab in the layer properties (combined with Hotlink). In the QGIS-2.2-UserGuide, on page 97, there is an example on how to do that in Linux (see image). What is the correct alteration of this for it to work in Windows? ![enter image description here](http://i.stack.imgur.com/2HEIy.jpg)


title: Geographic maps in ArcMap

text: I have OSM data in SHP format inside ArcMap 10. But in ArcMap they are stretched a lot East-West compared to what I see at www.openstreetmap.org. How can I change it in ArcMap? I am not looking for transforming the lat/long coordinates. MapInfo also shows the same SHP file like OSM do it.


title: Adding information from excel into ArcGIS

text: I am working in ArcGIS and have a excel sheet with data pertaining to properties. I want to add this data to a map in ArcGIS of which i have property boundaries as a shapefile? I have a common folio number between the two datasets. However I am a novice and all answers have been over my head. Please point me to a step by step process to inserting databases.


title: Export WebMap using PrintTask Failing Silverlight API

text: I am using a custom tool to export a webmap to PDF. The tool works successfully through ArcGIS 10.1. The problem crops up when consuming the tool through my Silverlight application using PrintTask. The tool works in an intermittent fashion i.e. sometimes it works, and sometimes it gives an error. I turned to ArcGIS Manager log and extracted the following error:               Traceback (most recent call last): File     "C:\arcgisserver\directories\arcgissystem\arcgisinput\Toolbox\ExportPDSitePlanPart.GPServer\extracted\v101\printingtoolbox\ExportPDSitePlan.py", line 48, in result = arcpy.mapping.ConvertWebMapToMapDocument(Web_Map_as_JSON, templateMxd) File "c:\program files\arcgis\server\arcpy\arcpy\mapping.py", line 234, in ConvertWebMapToMapDocument return webmap_tuple(*convertArcObjectToPythonObject(arcgisscripting._convertWebMapToMapDocument(webmap_json, template_mxd, notes_gdb, extra_conversion_options))) RuntimeError: Layer "PD Working Layer": Failed to create layer from service at http://10.58.161.70/arcgis/rest/serv...eatureServer/0. Layer "PD Working Layer": Failed to create layer from service at http://10.58.161.70/arcgis/rest/serv...eatureServer/0. Failed to execute (ExportSP). Failed to execute (ExportSP).      I can confirm that I can access the service from the same server where the tool is running. (i.e. 10.58.161.70). Can someone point me in the right direction? Thanks


title: Using Python script as precondition in ArcGIS ModelBuilder?

text: I have a model that has two preconditions, in which I need to figure out if there is a way to prioritize one precondition to 'run' prior to the second. Is anyone aware of any built-in ability to do so with Model Builder or that can suggest a work around? For instance, if I try to 'link' one precondition to another, I get the error "Can only connect a variable to a process", this is because my 'output' has a precondition that a python script needs to be run against it first, and then that output has the second precondition executed (but only once the original output has been altered by the Python script). ![enter image description here](http://i.stack.imgur.com/kfw3O.jpg)


title: Incomplete Data - Missing line segments from Oracle data

text: I have setup GeoServer 2.5.2 on Windows. Everything was going smoothly until I actually had to display my geometry data from an Oracle Database Table. Some lines are rendered but some are not rendered. The issue seems to be based on the type of lines (curves, etc). This issue does not happen if exporting from oracle to SqLite or Global Mapper. In that case, all lines are kept as in the original data from Oracle. The issue seems to be most likely with the oracle plugin (geoserver-2.5.2-oracle-plugin.zip). Has anyone else ran into this issue and if so, how did you manage to render everything coming from Oracle in GeoServer?


title: GRASS: populating attribute table with specific values

text: I have a vector map that includes 100 points. These points were created using `v.extract` to randomly sample a larger map of 10k points. As such, the 'cat' values range from 1 to more than 10k. What I would like to do is to update the attribute table by rescaling these values to 1 to 100. Here's what I've experienced. If I use               v.category input=pts100 output=foo option=del     v.category input=foo output=pts100 option=add --o      Then when the 'cat' values _appear_ rescaled when I print them with               v.category input=pts100 option=print      However, if I look at the dbtable with               db.select pts100      the 'cat' values have not changed (presumably because it is the key column). So I'm thinking I will need to add a new column using `v.db.addcol`, but I don't know how to specify the range of value to populate the new column. From what I can tell, I have to specify a single value. Ideas?


title: Joining table to shapefile with pyQGIS

text: Joining shapefiles is a relatively straightforward proess in QGIS using 'add vector join' in the layer properties menu. I can add a csv file as a vector layer, as well as a shapefile, and join the two using a common field. I would like to do a similar join using pyQGIS, however I can't seem to find any documentation on how to go about this. Is this possible to do programmatically using pyQGIS with 'add vector join' or the older 'Join by Attributes'? Or, alternately, might there be any other Python modules that might be able to merge a csv with a shapefile before loading as a map layer with pyQGIS?


title: How to correct the "Windows Activation Period has Expired" message on an Esri Amazon AMI?

text: I'm only getting this error with instances that were spun up with Esri AMIs. And I did find one post on the AWS forums (with no resolution) regarding this issue. Does anyone know what the troubleshooting/correction steps are for this issue? ~~I have not started any troubleshooting steps yet.~~   1. Restarted instance ( **No effect** ).   2. Checked KMS . ( **No KMS Issues** ).   3. Turned off firewall, modified C:\Program Files\Amazon\Ec2ConfigService\Settings\config.xml by setting `Ec2WindowsActivate` to `Enabled`.   4. Restarted Ec2Config service.   5. Success!    6. Turned firewall back on. ~~on but I still need to test the KMS...I probably still need to open up a port.~~   7. Firewall is back on and no issues. EDIT: AWS Forum Link: https://forums.aws.amazon.com/thread.jspa?threadID=119200 ![enter image description here](http://i.stack.imgur.com/SICQo.png)


title: QGIS 1.8 tries to use UTM in spite of geo projection file

text: I'm wrestling with a problem that I am sure I've created for myself, but cannot resolve.   I received a set of shapefiles that were expressed with UTM, and I want them to use lat/lon. QGIS will do the translation, of course, but there are hundreds of files and I didn't want to do this manually.   So I obtained the code for converting the UTM (N, Zone 14) to lat/lon and re- wrote the .shp files, taking care to change the bounding box values in the main header, the bounding box values in the record headers, and all the point values (these are polygons).   I replaced the .prj file with this, which reflects the constants I used in the conversion:               GEOGCS["GCS_North_American_1983",   DATUM["D_North_American_1983",    SPHEROID["GRS_1980",6378137,298.257222101]],    PRIMEM["Greenwich",0],    UNIT["Degree",0.0174532925199433]]      Yet QGIS treats the file as if it's UTM -- it displays nothing on the screen, and the cursor position is expressed as easting/northing. If I look at the layer properties, in the Metadata tab it shows this:               Extents:     In layer spatial reference system units : xMin,yMin 677384.09,4068257.57 : xMax,yMax   732884.09,4351757.57     Layer Spatial Reference System:     +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs      Am I correct that there appears to be dissonance between the min/max units and the projection? I cannot find anything that is telling QGIS to display this using UTM, but obviously,   I'm overlooking something.   GlobalMapper displays the file without complaint. Any thoughts on what I'm overlooking? Thanks.


title: How do I get and install the latest version of QGIS in Ubuntu?

text: I'd like to get the latest 64-bit dev build of QGIS on my Ubuntu setup. How exactly do I go about it? I assume it's some kind of arcane terminal command...


title: Different results given by PostGIS and MySQL?

text: I'm developing a spatial database of administrative areas. I started out with PostGIS, but as MySQL is being used in our other projects, I thought I would at least try to create a spatial-enabled MySQL database. However, spatial queries give me different results than PostGIS does. Let's say I have a query like this `SELECT nazev FROM orp as o, kraje as k WHERE Intersects(o.geom,k.geom) AND Within(Centroid(o.geom),k.geom) AND k.nazev = 'Nazev'`. It should return a list of administrative areas within a county with defined name (nazev). Instead, it gives me expected areas + areas that touch boundary of a county, but are not contained by this county. When I run the same query with PostGIS, I get the correct results. Have anyone ever faced this issue? I am using a national CRS (EPSG:102067) and have tried also with EPSG:4326 with same result. Isn't the implementation of the OGC spec the same as in PostGIS? Thanks for any hint.


title: WFS-T: changes not saved on an OpenLayers map

text: I have written a script allowing to display a few WMS and WFS layers and I want to make my 2 WFS layers editable. I have all my controls and can draw / modify / drag / delete my features but, on one of my layers, the changes are not saved, despite the saveStrategy set to "auto: true". I give my code here. The layer **adresspt_wfs** works fine, and the layer **bygg_ovr_wfs** can't save my new edits. If anyone has an idea?               var map;     var wmsurl = "http://localhost:8080/geoserver/wms";     var wfsurl = "http://localhost:8080/geoserver/wfs";     var projection = new OpenLayers.Projection("EPSG:900913");     var displayProjection = new OpenLayers.Projection("EPSG:4326");          function init() {       var bounds = new OpenLayers.Bounds(         1450000, 7490000, 1490000, 7510000);            var format = 'image/png';            var options = {         projection: projection,         displayProjection: displayProjection,         maxExtent: bounds,         units: 'm'       };            map = new OpenLayers.Map('map',options);            // All the base layers:            var osm = new OpenLayers.Layer.OSM('OpenStreetMap','',       { zoomOffset: 11, resolutions: [76.4370282714844,38.2185141357422,19.1092570678711,9.55462853393555,4.77731426696777,2.38865713348389,1.19432856674194] });            var gphys = new OpenLayers.Layer.Google(       "Google Physical",       {'sphericalMercator': true, minZoomLevel: 11, maxZoomLevel: 17, type: G_PHYSICAL_MAP});            var gstr = new OpenLayers.Layer.Google(       "Google Streets",       {'sphericalMercator': true, minZoomLevel: 11, maxZoomLevel: 17});            var ghyb = new OpenLayers.Layer.Google(       "Google Hybrid",       {'sphericalMercator': true, minZoomLevel: 11, maxZoomLevel: 17, type: G_HYBRID_MAP});            var gsat = new OpenLayers.Layer.Google(       "Google Satellite",       {'sphericalMercator': true, minZoomLevel: 11, maxZoomLevel: 17, type: G_SATELLITE_MAP});            var dem = new OpenLayers.Layer.WMS(       "ASTER Gdem", wmsurl,       { layers: 'ASTGTM2_N55E013_dem',           format: format,           srs:'EPSG:4326',           transparent:'false'       },       {isBaseLayer: true });            var nobase = new OpenLayers.Layer.WMS(       "No base map", wmsurl,       { layers: 'luma_project:background',         srs: 'EPSG:4326',         transparent: 'true'       },       {isBaseLayer: true, maxExtent: bounds, maxResolution: 'auto'});            // All the non-base layers (WMS):            var transportation = new OpenLayers.Layer.WMS(         "Lund transportation network", wmsurl,         { layers: 'lu_transportation',           format: format,           srs:'EPSG:4326',           transparent:'true',           zoomOffset: 11,         },         {isBaseLayer: false, visibility: true});            var adresspt = new OpenLayers.Layer.WMS(         "Adress points", wmsurl,         { layers: 'adresspt',           format: format,           srs:'EPSG:4326',           transparent:'true',           zoomOffset: 11         },         {isBaseLayer: false, visibility: false});            var bygg_ovr = new OpenLayers.Layer.WMS(         "Buildings", wmsurl,         { layers: 'bygg_ovr',           format: format,           srs:'EPSG:4326',           transparent:'true',           zoomOffset: 11         },         {isBaseLayer: false, visibility: false});            // All the non-base layers (WFS):            var styles = new OpenLayers.StyleMap({         "default": new OpenLayers.Style({         fillColor: '#ee9900'         }),         "select": new OpenLayers.Style({         fillColor: 'red'         })       });            var saveStrategy = new OpenLayers.Strategy.Save({auto:true});            var bygg_ovr_wfs = new OpenLayers.Layer.Vector("Buildings (edition)", {         stylemap: styles,         visibility: false,         strategies: [           new OpenLayers.Strategy.BBOX(),            saveStrategy           ],         protocol: new OpenLayers.Protocol.WFS({           url: wfsurl,           featurePrefix:"luma_project",           featureType: "bygg_ovr",           featureNS: "http://localhost:8080/luma_project",           srsName: "EPSG:900913",           geometryName: "the_geom",           version: "1.1.0"         })       });            var adresspt_wfs = new OpenLayers.Layer.Vector("Adress points (edition)", {         visibility: false,         strategies: [           new OpenLayers.Strategy.BBOX(),            saveStrategy           ],         protocol: new OpenLayers.Protocol.WFS({           url: wfsurl,           featurePrefix:"luma_project",           featureType: "adresspt",           featureNS: "http://localhost:8080/luma_project",           srsName: "EPSG:900913",           geometryName: "the_geom",           version: "1.1.0"         })       });            // Identification and edition of layers            // Edition            editControls = {         building_draw: new OpenLayers.Control.DrawFeature(bygg_ovr_wfs, OpenLayers.Handler.Polygon),         building_drag: new OpenLayers.Control.DragFeature(bygg_ovr_wfs),         building_modify: new OpenLayers.Control.ModifyFeature(bygg_ovr_wfs),         building_delete: new DeleteFeature(bygg_ovr_wfs, {title: 'Delete building'}),         building_click: new OpenLayers.Control.WMSGetFeatureInfo({             url: wmsurl,              title: 'Identify features by clicking',             layers: [bygg_ovr]         }),         building_hover: new OpenLayers.Control.WMSGetFeatureInfo({             url: wmsurl,              title: 'Identify features by hover',             layers: [bygg_ovr],             hover: true,             // defining a custom format options here             formatOptions: {                 typeName: 'bygg_ovr',                  featureNS: 'http://localhost:8080/luma_project'             }         }),         adresspt_draw: new OpenLayers.Control.DrawFeature(adresspt_wfs, OpenLayers.Handler.Point),         adresspt_drag: new OpenLayers.Control.DragFeature(adresspt_wfs),         adresspt_delete: new DeleteFeature(adresspt_wfs, {title: 'Delete adress point'}),         adresspt_click: new OpenLayers.Control.WMSGetFeatureInfo({             url: wmsurl,              title: 'Identify features by clicking',             layers: [adresspt]         }),         adresspt_hover: new OpenLayers.Control.WMSGetFeatureInfo({             url: wmsurl,              title: 'Identify features by hover',             layers: [adresspt],             hover: true,             // defining a custom format options here             formatOptions: {                 typeName: 'adresspt',                  featureNS: 'http://localhost:8080/luma_project'             }         })       };            for(var i in editControls) {         map.addControl(editControls[i]);         editControls[i].events.register("getfeatureinfo", this, showInfo);       }             // Add the layers on the map            map.addLayers([osm,gphys,gstr,ghyb,gsat,dem,nobase,transportation,bygg_ovr,bygg_ovr_wfs,adresspt,adresspt_wfs]);                 // Add the controls on the map            map.addControl(new OpenLayers.Control.Navigation());       map.addControl(new OpenLayers.Control.LayerSwitcher());       map.addControl(new OpenLayers.Control.OverviewMap());       map.addControl(new OpenLayers.Control.Scale($('scale')));             map.addControl(new OpenLayers.Control.MousePosition({element: $('location')}));                  map.zoomToExtent(bounds);               } // End of init()          function showInfo(evt) {            $('responseText').innerHTML = evt.text;     } // End of showInfo()               function toggleControlEdit(element) {       for(key in editControls) {         var control = editControls[key];         if(element.value == key && element.checked) {           control.activate();         }          else {           control.deactivate();         }       }     } // End of toogleControlEdit()          // This control has been found here:     // http://www.peterrobins.co.uk/it/oledit.html          DeleteFeature = OpenLayers.Class(OpenLayers.Control, {      initialize: function(layer, options) {         OpenLayers.Control.prototype.initialize.apply(this, [options]);         this.layer = layer;         this.handler = new OpenLayers.Handler.Feature(             this, layer, {click: this.clickFeature}         );     },     clickFeature: function(feature) {         // if feature doesn't have a fid, destroy it         if(feature.fid == undefined) {             this.layer.destroyFeatures([feature]);         } else {             feature.state = OpenLayers.State.DELETE;             this.layer.events.triggerEvent("afterfeaturemodified", {feature: feature});             feature.renderIntent = "select";             this.layer.drawFeature(feature);         }     },     setMap: function(map) {         this.handler.setMap(map);         OpenLayers.Control.prototype.setMap.apply(this, arguments);     },     CLASS_NAME: "OpenLayers.Control.DeleteFeature"     }) // End of DeleteFeature      Thanks in advance!


title: How to create a custom projection of PRS92 and UTM 51N for the Philippines?

text: > **Possible Duplicate:**   >  How to define new custom projections in QGIS? How do I create a customized CRS parameter in quantum GIS? I'm try to create a projection of PRS92 and UTM 51N. This projection is for the philippines. Thanks!


title: Sensitivity of Garmin eTrex 30

text: I have a very good experience with Garmin eTrex Legend HCx - here and here they write that HCx has a very sensitive chip, that works in forest, deep valleys and in cities - and I can only confirm - **with HCx I was able to get my position in plane!** **How sensitive is eTrex 30, related to HCx? Have you been able to catch satellites in room / plane / in a narrow street?** eTrex 30 doesn't say to be that sensitive at garmin website. (Could HCx then be more precise/usable in these conditions despite eTrex 30 supporting GLONASS?)


title: How do I change my title of app?

text: We have a map with streets as basemap. We have been trying unsuccessfully, for the last two days to change the basemap to tiled layer with no success at all. I am posting a snippet because script is long. Please let me know if you need additional script. On the surface, it looks simple but it has been extremely challenging. Your generous help is greatly appreciated in advance.


title: Finding percent landcover types with points or polylines

text: I'm using ArcMap 10.1. I had an orthophoto of an old study site. I traced two polygon shapefiles from this photo- one of "fields" and one of "forest". I made a point shapefile of the sampling points, and a polyline shapefile of all possible paths between the sample points. If a butterfly were to fly from one sample point to another (assuming linear flight along the polylines), I am looking to determine; a) what percent of that flight distance was through field, b) what percent of that flight distance was through forest and c) how many times that individual had to cross a forest/field edge. I haven't used GIS since my intro class many years ago, and I can't even figure out the correct vocabulary to search for the answer myself. Any and all help is greatly appreciated! Even if you don't know how to do this, suggestions of what to search for are equally appreciated!


title: measuring adjacency between polygons

text: I am working with some land ownership data. I would like to identify parcels the have small adjacency between other parcels. In a nutshell I want to identify where there is a small overlap between sections. In the example below you can see where one section overlaps another by a small amount of about 15 meters. In most situations this value of overlap is zero. ![enter image description here](http://i.stack.imgur.com/3p0Zy.jpg) Any suggestions how I can detect these situations? I am on arc 10.2 with Spatial Analyst.


title: pyproj, proj4 doesn't return same coordinates after inverse conversion

text: I'm trying to interpolate between coordinates, as a step in that process I convert latlong to UTM so I can work in metric. I don't get the expected latlong when converting back, reproduced below:               >>> import pyproj     >>> pyproj.__version__     '1.9.2'     >>> p = pyproj.Proj("+proj=utm +zone=33 +north +ellps=WGS84 +datum=WGS84 +units=m +no_defs")     >>> p(58.4, 15.6)     (5579373.890337878, 2328887.997220567)     >>> p(*p(58.4, 15.6), inverse=True) # Should return (58.4, 15.6)     (58.284316956873305, 15.63501250163002)      Is there a good explanation for this? How do I do what I want without this issue?


title: Updating filter values and update map

text: I'm developing an ExtJS3.4-GeoExt1-OpenLayers2.12 GIS application. For basic functionality I need to filter layers to get specific data from a Geoserver WFS layer. Filter is done through the variable myVar as a filter value as it follows: First Filter:               customFilter = new OpenLayers.Filter.Logical({        filters: [            new OpenLayers.Filter.Comparison({                type: OpenLayers.Filter.Comparison.EQUAL_TO,                   property: 'myCode',                   value: myVar            })         ]      });      Then the WFS Layer:               wfsLayer = new OpenLayers.Layer.Vector("LayerName", {           protocol: new OpenLayers.Protocol.WFS({                url: "http://localhost:8080/geoserver/wfs",            featureType: "LayerName",            featureNS: "http://localhost:8080/project",            featurePrefix: "schema",            geometryName: "the_geom",            outputFormat: "json",            readFormat: new OpenLayers.Format.GeoJSON(),            srsName: 'EPSG:25831',            defaultFilter: customFilter         }),         strategies: [new OpenLayers.Strategy.Fixed(),             myStrategy         ]     });      And add it to the map:               map.addLayer(wfsLayer);      This works with a pre-set myVar variable before running the app. My idea is that when myVar is changed dinamicaly in my app (this is done), wfs layer updates its filter value too so map data is updated. I've tried refreshing layer as myLayer.refresh() when myVar has been changed on my app but it still takes the initial myVar value so data in the map does not change. Anyone knows?


title: Create circles around a point on a map measuring the radius in hours

text: Disclaimer: I am not experienced in GIS and hope this question fits in this category. I was wondering whether it is somehow possible to create a circle* around a city, using Google Maps or the like, whilst the circle is not measured in miles (or other distance measure) but in hours of driving. The background is the following: I am looking for an ideal location to live if I have to commute between two different places. Ideally one place should be within a ~15 minute range and the other place in a ~45 minute range. I can't find an intersection manually. If there is an automated way how to do this, I could experiment with my parameters to find the ideal place. Any ideas? *Obviously it won't result in a circle.


title: WFS filter for Point in polygon (PIP) query

text: I am trying to construct a wfs request. I want to determine if a certain point is within a polygon. The point coordinates I want to send via the WFS request. The polygon is saved as a layer on my geoserver. How does a WFS request like this look like. My request looks like this, but I am getting an error:               http://myserver.com:8080/geoserver/mapicture/ows?service=WFS&version=2.0.0&request=GetFeature&typeName=mystore:dk_ger_nowater&cql_filter=INTERSECT(geom,%20POINT%20(55.689972%10.239258      The error says:               error on line 4 at column 160: xmlParseCharRef: invalid xmlChar value 16


title: Interpreting Moran's I results

text: I have a sample of data taken from regions of a country! i want to test if there are any spatial autocorrelation on my data using Moran Indice test. the null hypothesis: is that there is no spatial autocorrelation here a sample of my data: ![sample data for my variable with locations coordinates](http://i.stack.imgur.com/Xl9np.png) i used this steps to calculate it on r: library(RODBC) > setwd('e:/r/moran') > > channel <- odbcConnectExcel('moran.xls') > > data <- sqlFetch(channel, 'wilaya') > > inf.dists <- as.matrix(dist(cbind(lon=data$Lon, lat=data$Lat))) > > inf.dists.inv <- 1/inf.dists > > diag(inf.dists.inv) <- 0 > > library(ape) > > Moran.I(data$year2009, inf.dists.inv) and i got these results: $observed -0.02229578 $expected -0.02702703 $sd 0.03455708 $p.value 0.8911011 the $observed ~= $expected and they are negative! does that mean there is a little dispersion! but we cannot reject the null hypothesis, so there is no spatial autocoorelation between the regions! the problem is that i'm not sure about my interpretation ? and i want to know if the method i used for calculating the weight matrix for moran.I function is right?


title: Automatically create polygons from top right and bottom left coordinates?

text: I have a CSV file with hundreds of records. Each record contains 2 matching coordinates, 1) top right corner of a rectangle, 2) bottom left corner of a rectangle.               Polygon ID, Lat 1,       Lat 2,       Long 1,       Long 2               1,          27o 28' 28", 27o 33' 30", 86o 06' 04",  86o 15' 00"      Is there a method which I can use to automate the creation of these polygons? I anticipate there being a few steps.


title: UTFGrid from raster source

text: Is it possible to make utfgrids from a raster source (e.g. GeoTIFF)? I am trying to get utfgrids from a .tif with mapnik.xml and TileStache but it seems not to be working. As a workaround I am polygonizing the tif to a shape which works well but is not very efficient. Has anyone experince with this question? Maybe with other tileservers? This would be an example for the TileMill generated mapnik.xml I am working with:               <?xml version="1.0" encoding="utf-8"?>     <!DOCTYPE Map[]>     <Map srs="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over" maximum-extent="-20037508.34,-20037508.34,20037508.34,20037508.34">          <Parameters>       <Parameter name="bounds">5,43,18,50</Parameter>       <Parameter name="center">16.7693,48.4365,9</Parameter>       <Parameter name="format">png8</Parameter>       <Parameter name="minzoom">0</Parameter>       <Parameter name="maxzoom">15</Parameter>       <Parameter name="scale">1</Parameter>       <Parameter name="metatile">2</Parameter>       <Parameter name="id"><![CDATA[test]]></Parameter>       <Parameter name="_updated">1409125380000</Parameter>       <Parameter name="name"><![CDATA[tmp2m_2_grid]]></Parameter>       <Parameter name="tilejson"><![CDATA[2.0.0]]></Parameter>       <Parameter name="scheme"><![CDATA[xyz]]></Parameter>     </Parameters>               <Style name="utf_grid" filter-mode="first">       <Rule>         <RasterSymbolizer opacity="1" scaling="near"/>       </Rule>     </Style>     <Layer name="utf_grid"       srs="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over">         <StyleName>utf_grid</StyleName>         <Datasource>            <Parameter name="file"><![CDATA[/path/to/test.tif]]></Parameter>            <Parameter name="band"><![CDATA[1]]></Parameter>            <Parameter name="type"><![CDATA[gdal]]></Parameter>         </Datasource>       </Layer>          </Map>


title: Can you use ArcObjects in an MVC service at v10.1

text: I'm slighlty confused with the changes at 10.1 and understand this version no longer supports DCOM or remote arcobjects. Can you use ArcObjects in a webservice (in my case a .NET MVC 4 service) that is deployed on the ArcGIS Server host (so not remote) and uses ArcObjects to perform editing on a gdb featureclass? Or is the only way to go through the SOE framework? If you can use ArcObjects in the .NET MVC service are there any ramifications of doing so over using an SOE?


title: ArcGIS10 VBA requires library references

text: I just upgraded from ArcGIS 9.2 to 10.0, and I'm trying to create some tools using VBA. It doesn't recognize any of the ArcObjects Classes, Interfaces, or Constants unless the reference includes its library name. For example:               Dim pMxDoc as IMxDocument     Set pMxDoc = ThisDocument      doesn't work, but               Dim pMxDoc As esriArcMapUI.IMxDocument     Set pMxDoc = Project.ThisDocument      does. Does anyone know how to fix this? (FYI I have no admin capabilities on this computer.)


title: Polygon creation (Clockwise rotation or not)

text: According to "ESRI Shapefile Technical Description An ESRI White Paper—July 1998", Vertices for a single, ringed polygon are, therefore,always in clockwise order. Otherwise,"Dirty" polygons occur when the rings that define holes in the polygon also go clockwise, which causes overlapping interiors. In my understanding, you can create donuts polygon with "avoid intersection of new polygons" is checked on snapping option, even if both polygon(inner and outer) are created clockwise roated. It has not necessary to care about clockwise rotation when you create polygon in QGIS. Is this a right understanding?? Any advice are welcome!


title: How to compute and fill One-way field in a street shapefile according spacial position

text: I think it is something complicate, but I hope someone could bring up any idea about. I'm using ArcGIS10, and I have a "One-way roads" shapefile with a One-way field with Y/N. But **I need to create a new Field for this layer containing F or T if** :   * **F :** The final point of the line is Northern and Western than the origin of the line.   * **T :** The final point of the line is Southern and Eastern than the origin of the line.  An example,![Example](http://i.stack.imgur.com/p9vAR.png) Another Example: ![Example2](http://i.stack.imgur.com/qi1AL.png) Thanks for your time,


title: Geometry valid in PostGIS turns out to be not valid in Oracle

text: Geometry is in PostgreSQL database. It is valid. Geometry is exported as GML and then inserted into Oracle database. Then function SDO_UTIL.RECTIFY_GEOMETRY(geometry, tolerance) is called, which reports geometry as being not valid. More exactly, function               sdo_geom.validate_geometry_with_context(geometry, tolerance)       gives error               13367 [Element <1>] [Ring <1>].      I understand that Oracle uses tolerance, which means that Oracle sees some points as being same, which will cause errors. Using different tolerance makes errors go away. However changing tolerance is not an option, because I need application to work with the given tolerance. What should I do with geometries in PostGIS to ensure that they will be valid in Oracle after importing to Oracle? (provided I know the tolerance that will be used in Oracle beforehand). Examples:               select sdo_geom.validate_geometry_with_context(SDO_UTIL.FROM_GMLGEOMETRY('<gml:Polygon xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:3301"><gml:outerBoundaryIs><gml:LinearRing><gml:coordinates>564588.6,6438154.86 564592.9,6438149.4 564595.88,6438139.8 564604.32,6438069.13 564602,6438066.15 564595.38,6438064 564594.06,6438060.03 564593.56,6438047.95 564591.74,6438043.64 564573.37,6438036.69 564555.16,6438032.56 564549.87,6438033.55 564547.55,6438037.69 564546.23,6438043.98 564544.24,6438046.95 564540.1,6438047.78 564534.48,6438047.62 564531.83,6438046.13 564531.5,6438042.49 564533.65,6438024.45 564537.62,6438009.39 564540.44,6438003.43 564543.41,6438002.1 564546.89,6438001.44 564591.08,6438005.25 564604.32,6438006.24 564609.62,6438005.41 564609.85,6438005.02 564609.86,6438004.96 564611.6,6438002.1 564613.26,6437992.84 564616.9,6437964.04 564617.73,6437958.74 564616.9,6437952.12 564614.06,6437948.97 564582.31,6437943.02 564500.38,6437928.29 564433.52,6437916.54 564333.64,6437898.53 564243.69,6437883.77 564217.87,6437877.64 564209.59,6437878.47 564201.81,6437885.09 564197.51,6437893.53 564193.37,6437904.95 564184.27,6437928.45 564171.69,6437951.46 564164.24,6437960.73 564157.13,6437974.46 564142.56,6437999.79 564127.67,6438030.24 564126.18,6438036.69 564128,6438040.83 564131.47,6438043.31 564136.61,6438043.15 564145.54,6438041.33 564153.65,6438041.66 564175.83,6438043.81 564193.87,6438046.29 564277.45,6438055.89 564300.12,6438061.85 564342.49,6438075.26 564452.88,6438112.66 564509.65,6438131.36 564521.73,6438131.86 564526.2,6438127.89 564528.35,6438117.96 564531.33,6438113.16 564535.47,6438112.49 564549.37,6438116.3 564555.5,6438121.6 564559.63,6438134.34 564564.6,6438141.62 564571.72,6438142.78 564575.52,6438144.93 564578.17,6438150.89 564582.64,6438154.86 564588.6,6438154.86</gml:coordinates></gml:LinearRing></gml:outerBoundaryIs><gml:innerBoundaryIs><gml:LinearRing><gml:coordinates>564609.86,6438004.96 564609.58,6438005.41 564605.53,6438006.05 564607.18,6438005.22 564609.91,6438004.57 564609.86,6438004.96</gml:coordinates></gml:LinearRing></gml:innerBoundaryIs><gml:innerBoundaryIs><gml:LinearRing><gml:coordinates>564607.18,6438005.22 564606.51,6438005.37 564609.03,6438004.27 564607.18,6438005.22</gml:coordinates></gml:LinearRing></gml:innerBoundaryIs><gml:innerBoundaryIs><gml:LinearRing><gml:coordinates>564584.79,6438130.2 564579.49,6438131.2 564576.18,6438131.03 564571.88,6438129.04 564572.21,6438125.4 564575.52,6438124.25 564578.17,6438123.25 564580.32,6438122.59 564583.63,6438122.59 564585.29,6438124.08 564586.45,6438127.06 564584.79,6438130.2</gml:coordinates></gml:LinearRing></gml:innerBoundaryIs><gml:innerBoundaryIs><gml:LinearRing><gml:coordinates>564574.49,6438043.78 564573.81,6438043.92 564565.46,6438041.86 564565.87,6438040.02 564566.08,6438038.24 564566.56,6438038.17 564574.91,6438040.22 564574.91,6438040.43 564574.49,6438043.78</gml:coordinates></gml:LinearRing></gml:innerBoundaryIs><gml:innerBoundaryIs><gml:LinearRing><gml:coordinates>564418.96,6437980.92 564414.32,6437981.75 564410.02,6437980.92 564404.72,6437979.1 564401.41,6437974.96 564400.25,6437970.49 564402.57,6437968.01 564405.55,6437966.02 564411.18,6437963.87 564415.65,6437963.54 564419.62,6437964.7 564422.43,6437968.01 564424.09,6437972.15 564422.93,6437977.28 564418.96,6437980.92</gml:coordinates></gml:LinearRing></gml:innerBoundaryIs></gml:Polygon>'), 0.005) from dual;               select SDO_UTIL.RECTIFY_GEOMETRY(SDO_UTIL.FROM_GMLGEOMETRY('<gml:Polygon xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:3301"><gml:outerBoundaryIs><gml:LinearRing><gml:coordinates>564588.6,6438154.86 564592.9,6438149.4 564595.88,6438139.8 564604.32,6438069.13 564602,6438066.15 564595.38,6438064 564594.06,6438060.03 564593.56,6438047.95 564591.74,6438043.64 564573.37,6438036.69 564555.16,6438032.56 564549.87,6438033.55 564547.55,6438037.69 564546.23,6438043.98 564544.24,6438046.95 564540.1,6438047.78 564534.48,6438047.62 564531.83,6438046.13 564531.5,6438042.49 564533.65,6438024.45 564537.62,6438009.39 564540.44,6438003.43 564543.41,6438002.1 564546.89,6438001.44 564591.08,6438005.25 564604.32,6438006.24 564609.62,6438005.41 564609.85,6438005.02 564609.86,6438004.96 564611.6,6438002.1 564613.26,6437992.84 564616.9,6437964.04 564617.73,6437958.74 564616.9,6437952.12 564614.06,6437948.97 564582.31,6437943.02 564500.38,6437928.29 564433.52,6437916.54 564333.64,6437898.53 564243.69,6437883.77 564217.87,6437877.64 564209.59,6437878.47 564201.81,6437885.09 564197.51,6437893.53 564193.37,6437904.95 564184.27,6437928.45 564171.69,6437951.46 564164.24,6437960.73 564157.13,6437974.46 564142.56,6437999.79 564127.67,6438030.24 564126.18,6438036.69 564128,6438040.83 564131.47,6438043.31 564136.61,6438043.15 564145.54,6438041.33 564153.65,6438041.66 564175.83,6438043.81 564193.87,6438046.29 564277.45,6438055.89 564300.12,6438061.85 564342.49,6438075.26 564452.88,6438112.66 564509.65,6438131.36 564521.73,6438131.86 564526.2,6438127.89 564528.35,6438117.96 564531.33,6438113.16 564535.47,6438112.49 564549.37,6438116.3 564555.5,6438121.6 564559.63,6438134.34 564564.6,6438141.62 564571.72,6438142.78 564575.52,6438144.93 564578.17,6438150.89 564582.64,6438154.86 564588.6,6438154.86</gml:coordinates></gml:LinearRing></gml:outerBoundaryIs><gml:innerBoundaryIs><gml:LinearRing><gml:coordinates>564609.86,6438004.96 564609.58,6438005.41 564605.53,6438006.05 564607.18,6438005.22 564609.91,6438004.57 564609.86,6438004.96</gml:coordinates></gml:LinearRing></gml:innerBoundaryIs><gml:innerBoundaryIs><gml:LinearRing><gml:coordinates>564607.18,6438005.22 564606.51,6438005.37 564609.03,6438004.27 564607.18,6438005.22</gml:coordinates></gml:LinearRing></gml:innerBoundaryIs><gml:innerBoundaryIs><gml:LinearRing><gml:coordinates>564584.79,6438130.2 564579.49,6438131.2 564576.18,6438131.03 564571.88,6438129.04 564572.21,6438125.4 564575.52,6438124.25 564578.17,6438123.25 564580.32,6438122.59 564583.63,6438122.59 564585.29,6438124.08 564586.45,6438127.06 564584.79,6438130.2</gml:coordinates></gml:LinearRing></gml:innerBoundaryIs><gml:innerBoundaryIs><gml:LinearRing><gml:coordinates>564574.49,6438043.78 564573.81,6438043.92 564565.46,6438041.86 564565.87,6438040.02 564566.08,6438038.24 564566.56,6438038.17 564574.91,6438040.22 564574.91,6438040.43 564574.49,6438043.78</gml:coordinates></gml:LinearRing></gml:innerBoundaryIs><gml:innerBoundaryIs><gml:LinearRing><gml:coordinates>564418.96,6437980.92 564414.32,6437981.75 564410.02,6437980.92 564404.72,6437979.1 564401.41,6437974.96 564400.25,6437970.49 564402.57,6437968.01 564405.55,6437966.02 564411.18,6437963.87 564415.65,6437963.54 564419.62,6437964.7 564422.43,6437968.01 564424.09,6437972.15 564422.93,6437977.28 564418.96,6437980.92</gml:coordinates></gml:LinearRing></gml:innerBoundaryIs></gml:Polygon>'), 0.005) from dual;      EDIT: Added error message shown by Oracle.


title: What is the ruler accuracy in Google Earth Pro?

text: I use Google Earth Pro (7.1) on a daily basis to assess the acreage of US agricultural fields using the GE ruler. What is the accuracy we can expect from Google Earth Pro for surface calculation? (the fields don't always have simple shapes and are sometimes in remote areas) I found a few threads that partially cover this topic. But I didn't find any clear answer. Thank you


title: How to create polygons in ArcMap 10.0 using lat/long data from SQL Server 2008 database?

text: Forgive me if this is very basic, as I'm new to GIS work (though I've been a programmer for 10 years) and so I'm not really sure how to approach this, but here goes... Can I take a table with xmin/xmin ymin/ymax coordinates and use it to create a bunch of polygons in ArcMap 10.0? I have a connection set up to SQL Server 2008 wherein there is a table that has several pieces of information regarding various maps that our team keeps track of. Here's what I'm working with: Said table contains the lat/long data in the form of xmin/xmax ymin/ymax, in addition to several other attributes - gmap_id, author, year, scale, etc. My boss asked me if I could create a connection to our SQL Server, use that information to create polygons on a map of the US, and set up a little attribute window with the other information (author, gmap_id, scale, etc) for each polygon. I've got the first part covered, I think, thanks to this question. I set up the Data Query and pulled the columns that I initially wanted from SQLServer. That gave me a table like so: I had a screen cap but the image won't post due to my rep < 10, so the table looks like:               gmapId, scale, usgsIdF, xmin, xmax, ymin, ymax     18015, 250000, 13139, -121.891, -123.071, 43, 45.5     18016, 125000, 13140, -117, -118.936, 34.25, 35.624     ...      And I now have the data listed as a Source under the Table of Contents window. How can I now translate that information into a bunch of squares on the map? And then can I use that to set up an attribute table? I have a lot of python experience, so if that is the way to go I can do it. I am also good with SQL Server 2008, though I have not used it's spatial capabilities yet. So if there is a solution to be had involving TSQL I can also go that route. I'm just scratching around trying to figure a good way to start.


title: How to get raster corner coordinates using Python GDAL bindings?

text: Is there a way to get the corner coordinates (in degrees lat/long) from a raster file using gdal's Python bindings? A few searches online have convinced me that there is not, so I have developed a work around by parsing the gdalinfo output, it's somewhat basic but I thought it might save some time for people who might be less comfortable with python. It also only works if gdalinfo contains the geographic coordinates along with the corner coordinates, which I don't believe is always the case. Here's my workaround, does anyone have any better solutions? gdalinfo on an appropriate raster results in something like this midway through the output:               Corner Coordinates:     Upper Left  (  -18449.521, -256913.934) (137d 7'21.93"E,  4d20'3.46"S)     Lower Left  (  -18449.521, -345509.683) (137d 7'19.32"E,  5d49'44.25"S)     Upper Right (   18407.241, -256913.934) (137d44'46.82"E,  4d20'3.46"S)     Lower Right (   18407.241, -345509.683) (137d44'49.42"E,  5d49'44.25"S)     Center      (     -21.140, -301211.809) (137d26'4.37"E,  5d 4'53.85"S)      This code will work on files who's gdalinfo look like that. I believe sometimes the coordinates will be in degrees and decimals, rather than degrees, minutes and seconds; it ought to be trivial to adjust the code for that situation.               import numpy as np     import subprocess          def GetCornerCoordinates(FileName):         GdalInfo = subprocess.check_output('gdalinfo {}'.format(FileName), shell=True)         GdalInfo = GdalInfo.split('/n') # Creates a line by line list.         CornerLats, CornerLons = np.zeros(5), np.zeros(5)         GotUL, GotUR, GotLL, GotLR, GotC = False, False, False, False, False         for line in GdalInfo:             if line[:10] == 'Upper Left':                 CornerLats[0], CornerLons[0] = GetLatLon(line)                 GotUL = True             if line[:10] == 'Lower Left':                 CornerLats[1], CornerLons[1] = GetLatLon(line)                 GotLL = True             if line[:11] == 'Upper Right':                 CornerLats[2], CornerLons[2] = GetLatLon(line)                 GotUR = True             if line[:11] == 'Lower Right':                 CornerLats[3], CornerLons[3] = GetLatLon(line)                 GotLR = True             if line[:6] == 'Center':                 CornerLats[4], CornerLons[4] = GetLatLon(line)                 GotC = True             if GotUL and GotUR and GotLL and GotLR and GotC:                 break         return CornerLats, CornerLons           def GetLatLon(line):         coords = line.split(') (')[1]         coords = coords[:-1]         LonStr, LatStr = coords.split(',')         # Longitude         LonStr = LonStr.split('d')    # Get the degrees, and the rest         LonD = int(LonStr[0])         LonStr = LonStr[1].split('\'')# Get the arc-m, and the rest         LonM = int(LonStr[0])         LonStr = LonStr[1].split('"') # Get the arc-s, and the rest         LonS = float(LonStr[0])         Lon = LonD + LonM/60. + LonS/3600.         if LonStr[1] in ['W', 'w']:             Lon = -1*Lon         # Same for Latitude         LatStr = LatStr.split('d')         LatD = int(LatStr[0])         LatStr = LatStr[1].split('\'')         LatM = int(LatStr[0])         LatStr = LatStr[1].split('"')         LatS = float(LatStr[0])         Lat = LatD + LatM/60. + LatS/3600.         if LatStr[1] in ['S', 's']:             Lat = -1*Lat         return Lat, Lon          FileName = Image.cub     # Mine's an ISIS3 cube file.     CornerLats, CornerLons = GetCornerCoordinates(FileName)     # UpperLeft, LowerLeft, UpperRight, LowerRight, Centre     print CornerLats     print CornerLons      And that gives me:               [-4.33429444 -5.82895833 -4.33429444 -5.82895833 -5.081625  ]      [ 137.12275833  137.12203333  137.74633889  137.74706111  137.43454722]


title: Adding a new field in ArcMap from data of an existing field

text: I want to create a new field with years data. There is an existing field with the timestamp YYYYMMDDHHMM. But I just want the years from that field to the new field. That way I can display just the years when I create a graph. IS there anyone that can explain the code to do this. I am not a programmer and need the scripts or python code explained as if I were a 5-year old. Thank you.


title: ArcPy / ArcMap 10.2: How to report data to user?

text: I'm developing an arcpy addin with ArcMap 10.2 that outputs a feature layer based on user input. How can I display attributes from this layer in a MessageBox or other reporting mechanism? I don't want to enable data driven pages, so dynamic tables don't seem to be an option. I'm sure there's an easy solution out there. Ideas? Thanks! Katie


title: QGIS python plugin: PyQGIS for spatial join

text: For a QGIS 2.0 and 2.2 Python plugin, I am trying to update one layer's field attribute with another layer's field attribute based on geometry function QgsGeometry.within(). My first layer is a point layer, the second layer is a buffer of a vector polyline layer containing azimuth measurements. I would like to update the point layer to include azimuth information of the buffer polygon it is within (essentially a spatial join). It is to automatize the process described here: How to horizontally align irregularly shaped point symbol to vector layer in QGIS?. Currently, I'm not getting an error, but nothing is happening either.               rotateBUFF = my buffer polygon layer     pointLayer = my point layer to obtain azimuth data          rotate_IDX = rotateBUFF.fieldNameIndex('bearing')     point_IDX = pointLayer.fieldNameIndex('bearing')     rotate_pr = rotateBUFF.dataProvider()     point_pr = pointLayer.dataProvider()     rotate_caps = rotate_pr.capabilities()     point_caps = point_pr.capabilities()     pointFeatures = pointLayer.getFeatures()     rotateFeatures = rotateBUFF.getFeatures()          for rotatefeat in rotateFeatures:         for pointfeat in pointFeatures:             if pointfeat.geometry().within(rotatefeat.geometry()):                 pointID = pointfeat.id()                 if point_caps & QgsVectorDataProvider.ChangeAttributeValues:                     bearing = rotatefeat.attributes()[rotate_IDX]                     attrs = {point_IDX : bearing}                     point_pr.changeAttributesValues({pointID : attrs})


title: postgis st_linemerge based on attribute

text: I have a postgis table full of very small line segments. I would like to merge them according to an attribute in another field. update line_table set the_geom = ST_LineMerge(the_geom); but how do I do this according to the attribute so that all the small segments with name="west" are joined, and those with name="north" are all joined etc? Thanks for any help, Robert


title: Creating a model to run a complex query

text: I'm just a rookie when it comes to creating models in ArcGIS. I have two shapefiles from which I need to calculate an AVG and a MAX erosion rate for individual properties. The two shapefiles are: Coastal_Properties & Coastal_Change Coastal_Properties has a Property ID (PID), and Coastal_Change has a Rate of erosion in metres per year (RChngMYr). I'd like to retrieve an **Average and a Max RChngMYr for each individual PID** in the data-set. I have no problem running this query using SQL select in MapInfo, but am unsure how to go about it in ArcGIS. It will be better in the long run to create a model in Arc as Property boundaries are dynamic, and it would be terrible to have to run the query in MapInfo everyday. I'd like to create a Geoprocess service from this model and use it in a web application where it would return a phrase such as "The Max erosion rate for PID 12345 is 3.1 m/yr". Any help is appreciated! Thanks in advance


title: ToolValidator not doing anything

text: I have created a simple script tool which has a single parameter of type "String". I have written the tool validation to disable the input field (simply trying to learn tool validation). The validator code is as follows:               class ToolValidator:            def __init__(self):         import arcpy         self.params = arcpy.GetParameterInfo()            def initializeParameters(self):         self.params[0].Enabled = 0         return            def updateParameters(self):         return            def updateMessages(self):         return      Unfortunately, when I open the tool, the field is not disabled (as I expect it should be). Any ideas what I'm missing? SOLUTION: In case it's not clear in the threads below. The solution to the issue was to change Enabled to enabled, which is a case change between ArcGIS 9.3 and 10.0. This capitalization works correctly both in the initializeParameters method and the updateMessages method.


title: Python Script that calls the Voronoi Polygons Function?

text: I've been looking to see if there was a way to call Voronoi Polygons function found in Vector/Geometry Tools, directly from the code but I have not found one. Does this currently exist? Is there a way to call these functions within python coding?


title: How to flash a JSAPI map graphic only when contained in map extent

text: Using ESRI JSAPI 3.7 I'm flashing a graphic that works fine when the graphic is located within the bounds of the map extent (is visible). However, when you move the map so the graphic is located out of the map extent and try to flash, it throws the following error: "0x800a139e - JavaScript runtime error: lang.hitch: scope["setStroke"] is null (scope="[object Window]")". Here's the function I'm using. Can someone suggest maybe a way to flash only if the graphic is located onscreen? Thanks in advance!               function animateGraphic(graphicFlash) {                                               setTimeout(lang.partial(function (animateMe) {                            var shape = animateMe.getDojoShape();                            fx.animateStroke({                                shape: shape,                                duration: 500,                                color: { start: "yellow", end: [0,255,255] },                                width: { start: 25, end: 15 }                            }).play();                        }, graphicFlash), 100);                    }


title: Creating clumps of cells in a raster from seeds of 1 cell (in ArcGIS)

text: I need a process similar to that described in JPD's question: Creating randomly shaped clumps of cells in a raster from seeds of 1 cell/pixel? I have a raster with semi-randomly dispersed "seed" cells with values ranging from 1 to 19. I would prefer a function possible in ModelBuilder of ArcGIS, as this will be added to an already existing model.


title: How does QGIS Topology Checker check the topology of SHP and OSM files?

text: 1. Can anybody make me clear if .shp & .osm (OpenStreetMap Data) are holding topological data/relations or not.    2. if .shp files are not holding then how qgis plugin topology checker checks for the topology of these files.


title: Pyshp 1.2 example from documentation fails

text: I must be doing something really stupid because I am trying to run an example straight from the PyShp (shapefile.py) 1.2 documentation with the only change being the input shapefile and it is failing. Note that it is 32bit Python 2.6 on Server 2008 R2 64bit. The code:               import shapefile     e = shapefile.Editor(shapefile="F:/data/IntermediateWorkspace/5day_pgn.shp")     e.poly(parts=[[[5.1,5],[9.9,5],[9.9,1],[7.5,3],[5.1,1]]])     e.record("Appended","Polygon")     e.save('F:/data/IntermediateWorkspace/5day_pgn')      The error: > Traceback (most recent call last): File "F:\scripts\shapetest.py", line 2, > in e = > shapefile.Editor(shapefile="F:/data/IntermediateWorkspace/5day_pgn.shp") > File "C:\Python26\ArcGIS10.0\lib\shapefile.py", line 1048, in **init** > self.records = r.records() File "C:\Python26\ArcGIS10.0\lib\shapefile.py", > line 520, in records self.__dbfHeader() File > "C:\Python26\ArcGIS10.0\lib\shapefile.py", line 443, in __dbfHeader > fieldDesc = list(unpack("<11sc4xBB14x", dbf.read(32))) error: unpack > requires a string argument of length 32 Thanks for any help!!


title: old project on new PC: bad allocation when opening

text: Just got a new laptop (windows7). When opening a qgis project (shape layers on server), I get an error window "bad allocation. Project was made in Wroclaw, new PC has Lisboa. any suggestions?


title: Filling a Polygon with Lines in ArcGIS

text: i am trying to fill a polygon with lines and calculate total distance of these lines. http://i.stack.imgur.com/KxdqN.png like this example. but i have different distances between these lines like 1 meter, 2 meter, 3 meter.. so i need to take this as an array from client side, calculate length for each of these parameters and return a json response, should i use arcgis modal as a geoprocessing service? or Server object extensions? edit: actually i want to do this operation Filling a polygon with lines using Postgis in ArcGIS but not sure how


title: Prediction of coastal change

text: I have vector data in coastal change for two separate years, like 1950 and 2000. I want to use GIS to create a "Prediction of coastal change" in 50 years, 2050. How do I figure out what is the shape of the new coast line?


title: GRASS 6.4.3 Crashes on Startup; libintl-8.dll cited

text: I recently installed QGIS 2. I used the advanced OSGeo4W installer that I downloaded from the QGIS download page. GRASS is included in the stack. When I attempt to open GRASS, I receive the following error message: "The procedure entry point libintl_setlocale could not be located in the dynamic link library libintl-8.dll." I checked the WinGRASS errors page on the GRASS wiki and did not discover any precedent for this problem either there or by a Google search. My operating system is Windows Vista. I've installed all recommended packages as outlined on the GRASS wiki. Has anyone else encountered this problem? And how might I go about resolving it? Here's a running list of **potential solutions that have not corrected the problem**. I'll update as I get the time to try them:   1. Run Microsoft .NET cleanup utility and reinstall .NET 3.5.    2. Run system file checker from command prompt (>sfc /scannow). Thanks!


While reconciling a child version with the parent version, in a dataset containing Geometric network, the following error was noticed. Error during reconcile: "The version could not be reconciled. Invalid network element id : Field is not editable " We work on ArcGIS 10.1 . Tried to repair the Geometric network and do the reconcile , but to same effect. What is the core issue behind the error and what is the workaround. Experts pls advice.


title: Where can I find information about the development history of QGIS?

text: I know very little about early versions of `QGIS` (far far away from now..) and i'm feeling curious :-) I am looking for various information about the past of QGIS, including:   * Do you know in what year the software was created?    * By whom? (Was it a single person or a group? Do you have some names? Nationality?)


title: Custom Overpass Query in QGIS 2.0

text: QGIS 2.0 uses overpass to query osm data. I want to make customized queries like only buildings or only roads or both. And i only want the geometry so i'll go for skel mode in overpass. I may also need to set the timeout options for large data. How can i make these custom queries from qgis?


title: Problem with Georeferencer not scaling TIFcorrectly

text: I am trying to use QGIS and the Georeferencer to digitize paper maps but I am unable to get the raster to rotate afte I pick my control points. Below is a brief outline of the steps I used to transform my raster.   1. Added two control points and picked corresponding points in map canvas.   2. The transformation type is set to Helmert   3. Resampling method: Cubic   4. Compression: Deflate   5. Specified Output file name   6. Checked "Load in Qgis...", then Ok   7. Then clicked run button. The attached image is the results in QGIS. ![enter image description here](http://i.stack.imgur.com/VRbMo.png) The image should be much smaller


title: How to add legend in WFS openlayers

text: How can I publish WFS openlayer with legend, have tried different examples such as http://api.geoext.org/1.0/examples/vector-legend.js. But when I tried using this firebug shows an error=> GeoExt.MapPanel is not a constructor and the polygon for which the legend is defined in the rule(rule_highlight1 and rule_highlight2) has disappeared from the WFS layers.               var india = new OpenLayers.Layer.WMS(                     "cite:india_state - Tiled", "http://localhost:8080/geoserver/cite/wms",                     {                         LAYERS: 'cite:india_state',                         EPSG:4326,                         tiled: true,                         tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                     },                     {                         buffer: 0,                         displayOutsideMaxExtent: true,                         isBaseLayer: true,                         opacity: 0.8,                          }                 );                               var style = new OpenLayers.Style();                 var rule_highlight1 = new OpenLayers.Rule({                 title: "> 2000m",                 maxScaleDenominator: 3000000,                 filter: new OpenLayers.Filter.Comparison({                     type: OpenLayers.Filter.Comparison.EQUAL_TO,                     property: "classification",                     value: "1",                 }),                 symbolizer: {                     graphicName: "star",                     pointRadius: 8,                     fillColor: "#FE2EF7",                     graphicName: "star",                     fillOpacity: 0.6,                     strokeColor: "#FF0000",                     strokeWidth: 2,                     strokeDashstyle: "solid",                     label: " ${classification}",                     labelAlign: "cc",                     fontColor: "#000000",                     fontOpacity: 1,                     fontFamily: "Arial",                     fontSize: 16,                     fontWeight: "600"}                 });                           var rule_highlight2 = new OpenLayers.Rule({                 title: "< 2000m",                 maxScaleDenominator: 3000000,                 filter: new OpenLayers.Filter.Comparison({                     type: OpenLayers.Filter.Comparison.EQUAL_TO,                     property: "classification",                     value: "2",                 }),                 symbolizer: {                     fillColor: "#0040FF",                     graphicName: "star",                     pointRadius: 5,                     fillOpacity: 0.6,                     strokeColor: "#FF0000",                     strokeWidth: 2,                     strokeDashstyle: "solid",                     label: " ${classification}",                     labelAlign: "cc",                     fontColor: "#000000",                     fontOpacity: 1,                     fontFamily: "Arial",                     fontSize: 16,                     fontWeight: "600"}                 });                               var rule_highlight3 = new OpenLayers.Rule({                     filter: new OpenLayers.Filter.Comparison({                     type: OpenLayers.Filter.Comparison.EQUAL_TO,                     property: "classification",                     value: "2",                 }),                 symbolizer: {                     fillColor: "#0040FF",                     fillOpacity: 0.6,                     strokeColor: "#FF0000",                     strokeWidth: 2,                     strokeDashstyle: "solid",                     label: " ${classification}",                     labelAlign: "cc",                     fontColor: "#000000",                     fontOpacity: 1,                     fontFamily: "Arial",                     fontSize: 16,                     fontWeight: "600"}                 });                      style.addRules([rule_highlight1,rule_highlight2,rule_highlight3]);                 var wfs_layer = new OpenLayers.Layer.Vector("fsa",{                 strategies: [new OpenLayers.Strategy.BBOX()],                 styleMap: style,                 visibility: true,                 protocol: new OpenLayers.Protocol.WFS({                 version: "1.1.0",                 srsName: "EPSG:4326",                 url: "http://localhost:8080/geoserver/cite/wms",                 featureNS :  "http://www.opengeospatial.net/cite",                 featureType: "india_state",                 extractAttributes: true,                 transparent:true                 })                 });                      map.addLayers([india,wfs_layer]);                 map.zoomToMaxExtent();               mapPanel = new GeoExt.MapPanel({         renderTo: "map",         border: false,         layers: map,         center: [6.3, 45.6],         height: 256, // IE6 wants this         zoom: 8     });          legendPanel = new GeoExt.LegendPanel({         layerStore: mapPanel.layers,         renderTo: "legend",         border: false     });      }); Thanks!!!


title: Converting avl to sld or qml

text: Is there a simple way to convert a avl file into a sld or qml? Deegree does not work an my laptop (java problems) Or is there a way to open the original avl file in qgis? Thanks


title: Data Frame Property options have changed

text: I have Googled and searched ESRI help files and experimented until I am blue, not to mention reloading ArcGIS Desktop10 for Repair and adding the SP5 to it. I no longer have the options of border or background for data frame in layout view. I used to -- and now don't. I don't know what changed it yet it is impossible for me to combine data frame layers with the data frame borders visible. This is what I see when I right click on the data frame properties Both on the data frame in the layout view and on the data frame in the TOC in layout view: What happened and how do I fix it? Thanks. ![enter image description here](http://i.stack.imgur.com/MaLMc.jpg)


title: I don't seem to be able to save Layer Styles in QGIS

text: I have created and saved layer styles, not recently I admit but have done in the past. For some reason or other I just don't seem to be able to save any new ones in QGIS. No file is created. It does not even save changes to existing ones, even though it says I am replacing it. I surely must be doing something wrong, but don't know what. I can make the changes in the style and that's fine, it just won't save. I'm running on MS Vista.


title: Measuring straightness of a curve segment (represented as a polyline)

text: I'm working on an automatic elevation contour labeling algorithm and one of the factors that I want to take into account when deciding the positions of labels is how "straight" a particular segment of a contour is. The more straight it is, the more likely it is to be used to place the label on that segment. Each contour is represented by a polyline (but with points close together as to look like a curve to a naked eye). I then have a fixed length (width of a label), say, 100 pixels. If I randomly (or otherwise) choose a contour segment with the width of 100 pixels, I want to be able to get a numerical quantitative value of its straightness (say zero for a totally straight contour segment, some value larger than zero for a not so straight segment, and this value increasing as the crookedness increases). I've searched around for answers but I couldn't find anything really useful. I'd be grateful for any pointers.


title: Area of an EPSG 27700 polygon

text: I am working with EPSG 27700 polygon data in PHP. All the polygons are in the UK. The polygons do not have holes inside them. I'm hoping to be able to calculate the area of these polygons. For this I believe I need to take two steps. Firstly to convert the coords into something where x an y are the same distance and then once that is done to calculate the area inside each polygon. My naive solution would be to use some graphics package to draw the polygon onto an image and count the coloured squares.   1. What's a good coordinate system for this? UK Easting & Northing is my best guess.   2. Is there a PHP library or easily(!) implemented algorithm to find the area inside a polygon defined by a list of points? Alternatively, should I be getting out of my comfort zone and learning to use a proper GIS system?


title: How to generate automatically the name of output parameters in a python script tool?

text: I have some Python script tools in ArcGIS for Desktop. I need to reproduce what all geoprocessing tools do when the first input parameter is introduced: autogenerate a workspace and basename for the outputs. How can I do that?


title: How to merge Shapefiles from multiple folders into Geodatabase Featureclass?

text: I have a question regarding the merge of a lot of folders, which has a lot subfolders, which has the shapefiles stored. I would like to merge all those in one geodatabase. So I thought to use a Python script.               import os, arcpy.da          print os.getcwd()          for dirname, dirnames, filenames in os.walk('.'):         for subdirname in dirnames:             current_dir =  os.path.join(dirname, subdirname)             arcpy.env.workspace = current_dir             fcList = arcpy.ListFeatureClasses(".shp")             Merge = r"D:\Test\Merge_multiple.gdb\Totals"             for fc in fcList:                 print fc                 arcpy.Merge_management([fc,destination],Merge)                 arcpy.Delete_management(desination)                 arcpy.Rename_management(Merge, destination)             break      So, I started the script and after a minute it stops and unfortunately, no result. Because there are several different names in the shapefiles, I thought to use only the suffix .shp. Could you help me? [Edit] I get the path of the working directory in the python window. After this message, the window closes without any errors. ![enter image description here](http://i.stack.imgur.com/H6i7M.png) [Edit2] Based on PolyGeo's advice, I've run the script also in the IDLE environment and this is the result: ![enter image description here](http://i.stack.imgur.com/o42jt.png)


title: ArcObjects : How to identify ITool.OnMouseUp method execution

text: In ArcMap, I have a tool that inherits `BaseCommand`. The tool enables the user to draw a graphic on the screen. I activated this tool (set as current tool) from a Windows UserForm by getting its reference using `ICommandBars.find` method. Now how can the windows form get a notification when the `ITool.OnMouseUp` method is complete? I see that there are no events generated by the `ITool` interface.


title: How to create and export complex symbology?

text: I've got punctual geometries that I'd like to represent in a desktop and a web application, using complex symbology, like this : ![enter image description here](http://i.stack.imgur.com/Vy5sF.png) My idea is to use a symbology editor to create a SLD file that defines the symbols. I used QGis to create the symbology and then tried saving it into a SLD file. The problem is that symbols are simplified and instead of the symbol I created, I get this one: ![enter image description here](http://i.stack.imgur.com/yc1lz.png) Do you know about a SLD editor that supports complex symbology? Do you have any other suggestion different from creating a SLD file?


title: length array in OpenLayers.Layer.Vector

text: how to find length array (markers) in OpenLayers.Layer.Vector ?               var data = new OpenLayers.Layer.Vector("Data", {         styleMap: styleMap,         strategies: [new OpenLayers.Strategy.Fixed()],         protocol: new OpenLayers.Protocol.HTTP({             url: "./data/data.geojson",             format: new OpenLayers.Format.GeoJSON({                 extractStyles: true,                  extractAttributes: true,                 maxDepth: 2             })         })     });


title: Display local language for OSM data in TileMill

text: I'm following the instructions here to create a minimalist style map for Japan: https://github.com/aaronlidman/Toner-for-Tilemill Now, I've downloaded the OSM data, used _osm2pgsql_ , created the tilemill project and it now properly displays in TileMill. However, all the displayed city name labels are in the Romanized (english) characters. Is there a way that I can get the Japanese (Kanji, hiragana, etc) displayed in place of thier romainzed dipictions? ( _Street names do appear in the original Japanese, sometimes with Romanized letters in parens_ ) I'm assuming that the original Japanese names are contained in the osm.pbf, but somewhere along the pipeline they are being ignored. Any ideas?


title: Problems with the Clip and Intersect tool

text: I am having a problem with the clip and the intersect tool on QGIS Valmiera. I have a large vector map of sublocations of a country and am trying to clip them to a district level since I do not need the data for the entire country. I have save both to the project CRS, and have done this multiple times. As a project CRS I am using the WRS 84/Pseudo Mercator (I am using the openlayers plugin for geographic reference). To clip I the large map goes into the Input vector area, and I have made a polygon to represent the district which i add in the "clip layer" field. Then when I clip (or use intersect), nothing ever appears. The layer is created, but is not seen anywhere on the map, and the attribute table, although the fields from the large map are present, has no values. How can I fix this problem?


title: Why do all rows go missing from File Geodatabase Feature Class to *.dbf conversion?

text: ![enter image description here](http://i.stack.imgur.com/xKIul.png)I tried to export a File Geodatabase Feature Class to a *.dbf table/shapefile, but there are no records inside. I believe there is no mistake in exporting--it is simple. The doubt is the file geodatabase format--I cannot use it as a shapefile or export it as a shapefile with intact data. Attached is the original data and the exported data. I am using ArcGIS 10. The data are from the **United States census**. I am only interested in one data column/Field from each state. Is there a fast way to open, export, or convert the geodatabase format to *.dbf?


title: GROUP BY clause in the UpdateCursor arcpy.da

text: I need to specify a "group by" clause for UpdateCursor arcpy.da. From here I found out that it's possible to do using an optional pair of SQL prefix and postfix. But there are no examples and I'm stuck with the syntaxis problem. Can, please, smbd help me? PS. I've tried various modifications of `"(None, 'GROUP BY "INPUT_FID"')"`, but still can't get it. Thank you in advance. Here is the part of the code I'm using:               oDSFields = ('INPUT_FID','NEAR_FID','DISTANCE', popTotalField, popBlacksField)          kernel = 0     alpha = 0.005     #, "(None, 'GROUP BY "INPUT_FID"')"     with arcpy.da.UpdateCursor("outDistanceSorted", oDSFields, "(None, 'GROUP BY "INPUT_FID"')") as cursor:              for row in cursor:                  kernel = float(math.exp(-0.5*math.pow(row[2]/alpha, 2)))


title: Translate Excel formula to be used in field calculator

text: I have a series of equations to create corner points to make rotated polygons based on angles in the attribute table, but I have created it in excel and can't figure out the syntax to get it to work properly in the field calculator. I have 8 different trigonometric formulas to create XY coordinates for the 4 corner points. One of them is as follows:               =COS(ATAN(1500/2250)+RADIANS(MOD(360-([Angle]+270),360)))*(SQRT((1500^2)+(2250^2)))+[X]      Once I understand the syntax for one, I should be able to use it for the other 7 formulas, even though they vary slightly.   Can anyone advise the best way to use this formula in field calculator? Thanks. **Edit** , after receiving @Paul's answer: I figured out the correct syntax of some of the other ones, but two of them are still giving me trouble. One is producing a different result than Excel and one is saying that the syntax is incorrect.               =COS(RADIANS(270)-ATAN(750/1500)+RADIANS(MOD(360-([A]+270),360)))*SQRT((1500^2)+(750^2))+[X]          =COS(-ATAN(1500/2250)+RADIANS(MOD(360-([A]+270),360)))*SQRT((1500^2)+(2250^2))+[X]      I think it has to do with the negative signs and me probably putting them in the wrong spot...


title: Loading Image Services to ArcGIS for Server?

text: I am publishing Image Services to our ArcGIS for Server. As a mosaic dataset, the image has no issue. Overviews are created, and the footprints read proper elevation values. The issue I am running into is when loading into ArcGIS for Server as an Image Service. Side note, the mosaic is a referenced mosaic dataset. Among 5 other product derivatives referenced, only two are having this issue, the other 3 referenced mosaics are loading just fine. Even the original mosaic dataset is "clean" If anyone could offer insight it'd be much appreciated ![Mosaic Dataset](http://i.stack.imgur.com/oVgoY.png) ![Issue when published as an Image Service](http://i.stack.imgur.com/n6AsB.png)


title: Check for extraneous bytes that appear before %PDF at the head of a MapInfo printed PDF file?

text: I am using a MapBasic program that automatically prints individual PDF maps of each polygon in my table. After getting an Adobe Reader XI update, the PDFs no longer open. I receive the following error: The document is damaged and cannot be repaired. Adobe Reader could not open because it is either not a supported file type or because the file has been damaged (for example, it was sent as an email attachment and wasn't correctly decoded). The workaround from Adobe suggests, "modify the way the PDF is produced. Make sure that no extraneous bytes appear before %PDF at the head of the file." Is there a way to check this in MapBasic and modify the program so it correctly writes the file?


title: Geoserver/java out of heap space - how do I increase it?

text: I'm new to Geoserver, Ubuntu and Java, but have downloaded a virtual machine from gisvm.com and am getting up to speed. I got as far as configuring it with some fairly large shapefiles from a project I've worked on previously. My question is related to a problem which I see if I show the shapefile using the OpenLayers layer preview option. I see an error 'OpenLayers map preview code="internalError" Rendering process failed. java.lang.OutOfMemoryError: Java heap space'. Googling has led me to plenty of Java command line options to increase heap space, but I have no idea if this should be applied to an environment variable, in a startup script or as a part of the Geoserver config. Can you help me to understand what I need to edit to get this working? I'm also wondering if I should be splitting my shapefile into smaller pieces and so any guidelines on this matter would be helpful too. Thanks!


title: How to create an interactive web mapping site with comment/rating system?

text: I'm trying not to ask something that's already covered by the other 'beginning web mapping' questions, but I want to do something that I don't see covered here yet. Maybe I'm just missing it. Basically, I want to create an interactive web mapping site. Here's the way I picture it: there is a main map, showing points (pretend they're campsites). When a user clicks on a campsite, in another pane info on that campsite pops up, along with the ability to make comments on it. The user can also add other points, and draw rough polygons, lines, and add text boxes, that all other users see. It will be kind of a giant scratchpad, or scribble map that all users can contribute to. I`d also like to add some kind of thumbs-up, thumbs- down rating system for any comments on the individual campsites, and overall rating system for each user (this looks pretty daunting to me after typing it out!) I'm new to web mapping, as you may have guessed, and I think I know how to proceed to get the maps online, but it's the interactive part that I could really use some general direction on. I have mid-level skills with the basic tools: html, css, php, mysql, if that helps. Many thanks!


title: How to change the page layout from portrait to landscape using Python? (ArcMap 10)

text: I need my script to access the Page Orientation property of the current MXD, and switch it from portrait to landscape (or back). I've looked through the arcpy and arcpy.mapping documentation, and haven't had any luck. Does anyone know how to do this?


title: How do I upload a Featureclass to a Geodatabase in C#?

text: All, I have a featureclass with Arcobjects that came from a Shape file and I wish to import it into a geodatabase. I have the geodatabase connection, I just don't know how to load the feature class to a Geodatabase. Any help would be amazing! Im looking for a code example seeing I can't find it by googling. Hopefully in C#, but I don't mind reading through VB if needed. Thanks. **EDIT:** I wanted to add the Answer thanks to George so here is a code snippet with comments.                ESRI.ArcGIS.Geodatabase.IWorkspace pWorkspace;             ESRI.ArcGIS.Geodatabase.IWorkspaceFactory pWorkspaceFactory;             ESRI.ArcGIS.Geodatabase.IFeatureWorkspace SdeFeatureWorkspace;                  pWorkspaceFactory = new FileGDBWorkspaceFactory();             pWorkspace = pWorkspaceFactory.OpenFromFile(@"C:\tests\test.gdb", 0); //Opens the Geodeabase and creates a connection.             SdeFeatureWorkspace = (ESRI.ArcGIS.Geodatabase.IFeatureWorkspace)pWorkspace;                  IWorkspaceFactory workspaceFactoryShape = new ESRI.ArcGIS.DataSourcesFile.ShapefileWorkspaceFactoryClass(); //creates a shape factory to pull the shape and create the feature class.             IFeatureWorkspace featureWorkspace = (IFeatureWorkspace)workspaceFactoryShape.OpenFromFile(@"C:\Users\Desktop\Newfolder", 0); //Opens the Shape file.             IFeatureClass featureClass = featureWorkspace.OpenFeatureClass("al012009.001_5day_lin"); //opens the featurclass from the shapefile                  SdeFeatureWorkspace.CreateFeatureClass("al012009001_5day_linnn", featureClass.Fields, featureClass.CLSID, featureClass.EXTCLSID, featureClass.FeatureType, featureClass.ShapeFieldName, null); //creates a feature class in the geodb.                  Geoprocessor gp = new Geoprocessor(); //starts a GEOprocess.             gp.SetEnvironmentValue("workspace", @"C:\Tests\test.gdb"); //sets the geodb to use                  ESRI.ArcGIS.DataManagementTools.Append  append = new ESRI.ArcGIS.DataManagementTools.Append(); //this tool uploads the feature class             append.inputs = featureClass;             append.target = @"C:\tests\test.gdb\al012009001_5day_linnn";                  gp.Execute(append, null);


title: How to programmatically get a reference to an ICommand inside ArcGIS?

text: I'm having problems trying to get a ICommand using the method ICommandBars.Find. The ICommand is there, but it's returning null to me. Any ideas?                   public static ICommandItem GetCommandItemByName(string progId)         {             UID uid = new UIDClass();             uid.Value = progId;                  Type t = Type.GetTypeFromProgID("esriFramework.AppRef");             object o = Activator.CreateInstance(t);                  IApplication application = o as IApplication;             ICommandBars commandBars = application.Document.CommandBars;             ICommandItem commandItem = commandBars.Find(uid,false,true);                  return commandItem;         }      The result is always null :(. I need to get a reference to this specific command because it contains some fields and properties that I specified (a IWorkspace object, if it matters) Thanks for your help!


title: Work flow for georeferencing imagery using open source tools?

text: Does anyone have a recommendation for a work flow for georeferencing high resolution aerial photography using open source tools? Comparison of user friendliness and robustness between Quantum GIS versus GDAL? Thank you!


title: How to move a Point perpendicular from nearest Line?

text: I want to move Point perpendicular from nearest line or polygon, and the question is how to do that? example, I have points (house number) near roadbed, and I don't want them to intersect roadbed (buffered street line) I have used near analysis to get street line angle to the points for correct labeling, so the angle of the labels is really good but they are to close the road. So I want to move the points perpendicular away from street line. is there any way?


title: Store Geometry in a database with the ESRI Silverlight API

text: Is it possible to store geometry objects in a database using the ESRI Silverlight API? For this particular case, the client may have SQL Server 2005 installed instead of 2008 and therefore I can not use spatial types. Any solution would also have to support clients who use ArcGIS Server 9.3, and therefore we would not be able to use the new editing capabilities in ArcGIS 10.


title: Store global variables in ArcPad that store user information

text: I know that it's possible to store global variables in ArcPad from reading the help files. I have several field staff that are out everyday -- generally we have them select their name from a "Specialist" domain for each data point they create. Is it possible to store a table somewhere in the background of ArcPad that stores the user's information that can be used to auto-populate fields such as Specialist (user-name)?