title: Problem circle with zoom level

text: when I zoom level 16 and over I have two circle segments, I don't know if it's normal ( _projection_ ) or a bug. Have you an idea ? Openlayers : **2.13.1**   Problem With Chrome, IE and Firefox formula circle               function CircleMax(radius)     {     var pointsCircle = [];           var side = 160;      var angle;     var projectedLonLat, geometryPoint;     var points = [];         for (var i = 0; i < side; i++)     {         angle = (i * 360 / side) + side;         projectedLonLat = OpenLayers.Util.destinationVincenty(centerCircle, angle, radius);         projectedLonLat.transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:900913"));                     geometryPoint = new OpenLayers.Geometry.Point(projectedLonLat.lon, projectedLonLat.lat);         pointsCircle.push(geometryPoint);     }     var layerCircleStyle = OpenLayers.Util.extend({}, OpenLayers.Feature.Vector.style['default']);           var style_poly = OpenLayers.Util.extend({}, layerCircleStyle);       style_poly.strokeColor = "red";     style_poly.fillColor = "transparent";     var linearRing = new OpenLayers.Geometry.LinearRing(pointsCircle);     var polygonFeature = new OpenLayers.Feature.Vector(new OpenLayers.Geometry.Polygon([linearRing]),{},style_poly);     circleMaxLayer.addFeatures([polygonFeature]);      } Zoom level 6 ![Zoom le 6 ](http://i.stack.imgur.com/fKj0y.png) Zoom level 16 and over ![image description here](http://i.stack.imgur.com/2dVqA.png)


title: Using multiple basemaps with masking/clipping?

text: Is there a way to limit a basemap to within a boundary in arcmap 10? I am working on a project within San Joaquin county in California. I want the areas surrounding the county to be displayed with a light grey basemap with minimal details. Within the county itself, however, I want to use a more detailed basemap that has the information I am trying to portray. Is there a way to add a second basemap but have it be limited to the county boundaries? Thank you so much for any help you can give me.


title: how to solve Arcgis API for Java Script - Asp proxy error

text: I am using ArcGis 10.1 for ArcGIS Api for Javascript. I have faced the problem of {"code":400,"message":"Unable to complete operation.","details":[]}. while running the application in IIS and my Proxy.config file is look like this:               <serverUrl url="http://localhost:6080/arcgis/rest/services/" --------------------myservices                matchAll="true"></serverUrl>      How can i rectify this error..


title: Null values coerced to zero (0) when exporting to csv or txt

text: I extracted values from several rasters, and in some cases, the value is null. However, when I export the attribute table to a csv or a txt file (from ArcMap), the null values are converted to 0. Is there any way to prevent this from happening? These zeros are a real downer, as throw off my analysis, which is done in R.


title: convert XML to PostgreSQL using MATLAB

text: In a weeks time I should be able to transfer the GIS data provided in the XML format to PostgreSQL database. I have prepared the PostgreSQL database (empty) and looking at Matlab scripts. (I am on ubuntu platform). Matlab is our existing platform.The work is passed on to me by another person who used Mapping tools on MATLAB to export the data to viewers like ArcGIS/Google Earth.This was created for a prototype. (However, in the final product I will be exporting PostgreSQL to QGIS to view the data). I am pretty new to the GIS data/table.Could you suggest me the steps in-detail to have the data in PostgreSQL?


title: which framework is good for closed intranet

text: My requirement is simple. By latitude/longitude, draw icon on the map, draw a line between two icons.I may need put a customer picture over the map. It is not difficult by google map javascript or bing map javascript. The problem is that the system is running in a closed intranet environment. No internet access. Which framework shall I choose? Does it provide map? The system only focus on a specific area, such like a national park, it does not need the map for the whole world. Thanks!!!


title: How to apply Log Transform to Flow Accumulation?

text: I have run a flow accumulation for a DEM that I'm working with. I'm trying to apply a log transform because I've heard it could give a better visual representation. How would I go about doing this process.... Can this be done without Python?


title: AVG and GROUP BY in ArcMap

text: I have a Road Network shp file.(attributes include Names and Volumes ) Each road in my network consists of many links with the same name and different volumes. I want to display Road Names and volumes as labels using sql query or python function to group by names and have volume Average. ** In normal SQL query I would type:               SELECT ROAD_NM, AVG(VOL) as VOL     FROM TABLE     GROUP BY ROAD_NM;


title: How can I convert a DEM into areas of "flatness", with an arbitrary fill/cut to area ratio?

text: I want to find areas of relative "flatness" in a DEM. How can I find flat areas? I want to find a neighbourhood of contiguous cells that have an elevation variance that is within a certain tolerance threshold. I can use the ESRI suite of tools, or any open source GIS.


title: Find maximum and minimum value of points within a polygon buffer?

text: I have a series of 500m buffer zones along a polyline. I also have a series of points following the line, each point represents a 10m x 10m cell. I calculated the values of a field in each point known as "F18". I want to find the average, maximum and minimum values in F18 for all points within each buffer zone and add fields to the attribute table of the buffer for each. ArcGIS 10.1, arcinfo license.


title: Custom scalebar units in ArcMap

text: Is it possible to create a scalebar in ArcMap that shows distances in custom units? I understand how to do this manually, but it would be great if I could define a custom unit of measurement (let's say "1 swimming pool") and use that to dynamically generate a scalebar.


title: Django QGIS API, PyQt4 integrated

text: I am developing a web page through Django, where I want to insert images as static files. These images are produced with QGIS commands, which exist in QGIS API, PyQt4, and other API. So, my question is: Is it possible to incorporate the QGIOS libraries in Django? How do I do that? Thanks!


title: Eclipse export as runnable jar file error

text: I am building a destop application , which takes user input for personal geodatabase and finds feature class in it . It works fine till i run it in eclipse , but when i export it as a runnable jar file from eclipse Juno export option , it export it with warning and when i run the jar file from command promt i get following error ---               java.lang.IllegalArgumentException: URI is not hierarchical         at java.io.File.<init>(Unknown Source)         at com.esri.arcgis.interop.NativeLoader.getArcGISHome(Unknown Source)         at com.esri.arcgis.interop.NativeLoader.loadLibrary(Unknown Source)         at com.esri.arcgis.system.EngineInitializer.initializeEngine(Unknown Source)         at default.FeatureClassFinder.main(FeatureClassFinder.java:69)         at default.BrowsePath.actionPerformed(BrowsePath.java:139)         at javax.swing.AbstractButton.fireActionPerformed(Unknown Source)         at javax.swing.AbstractButton$Handler.actionPerformed(Unknown Source)         at javax.swing.DefaultButtonModel.fireActionPerformed(Unknown Source)         at javax.swing.DefaultButtonModel.setPressed(Unknown Source)         at javax.swing.plaf.basic.BasicButtonListener.mouseReleased(Unknown Source)         at java.awt.Component.processMouseEvent(Unknown Source)         at javax.swing.JComponent.processMouseEvent(Unknown Source)         at java.awt.Component.processEvent(Unknown Source)         at java.awt.Container.processEvent(Unknown Source)         at java.awt.Component.dispatchEventImpl(Unknown Source)         at java.awt.Container.dispatchEventImpl(Unknown Source)         at java.awt.Component.dispatchEvent(Unknown Source)         at java.awt.LightweightDispatcher.retargetMouseEvent(Unknown Source)         at java.awt.LightweightDispatcher.processMouseEvent(Unknown Source)         at java.awt.LightweightDispatcher.dispatchEvent(Unknown Source)         at java.awt.Container.dispatchEventImpl(Unknown Source)         at java.awt.Window.dispatchEventImpl(Unknown Source)         at java.awt.Component.dispatchEvent(Unknown Source)         at java.awt.EventQueue.dispatchEventImpl(Unknown Source)         at java.awt.EventQueue.access$200(Unknown Source)         at java.awt.EventQueue$3.run(Unknown Source)         at java.awt.EventQueue$3.run(Unknown Source)         at java.security.AccessController.doPrivileged(Native Method)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.awt.EventQueue$4.run(Unknown Source)         at java.awt.EventQueue$4.run(Unknown Source)         at java.security.AccessController.doPrivileged(Native Method)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.awt.EventQueue.dispatchEvent(Unknown Source)         at java.awt.EventDispatchThread.pumpOneEventForFilters(Unknown Source)         at java.awt.EventDispatchThread.pumpEventsForFilter(Unknown Source)         at java.awt.EventDispatchThread.pumpEventsForHierarchy(Unknown Source)         at java.awt.EventDispatchThread.pumpEvents(Unknown Source)         at java.awt.EventDispatchThread.pumpEvents(Unknown Source)         at java.awt.EventDispatchThread.run(Unknown Source)     Error: null     Sample failed.  Exiting...     java.lang.NullPointerException         at com.esri.arcgis.interop.NativeLoader.getArcGISHome(Unknown Source)         at com.esri.arcgis.interop.NativeLoader.loadLibrary(Unknown Source)         at com.esri.arcgis.system.EngineInitializer.initializeEngine(Unknown Source)         at default.FeatureClassFinder.main(FeatureClassFinder.java:69)         at default.BrowsePath.actionPerformed(BrowsePath.java:139)         at javax.swing.AbstractButton.fireActionPerformed(Unknown Source)         at javax.swing.AbstractButton$Handler.actionPerformed(Unknown Source)         at javax.swing.DefaultButtonModel.fireActionPerformed(Unknown Source)         at javax.swing.DefaultButtonModel.setPressed(Unknown Source)         at javax.swing.plaf.basic.BasicButtonListener.mouseReleased(Unknown Source)         at java.awt.Component.processMouseEvent(Unknown Source)         at javax.swing.JComponent.processMouseEvent(Unknown Source)         at java.awt.Component.processEvent(Unknown Source)         at java.awt.Container.processEvent(Unknown Source)         at java.awt.Component.dispatchEventImpl(Unknown Source)         at java.awt.Container.dispatchEventImpl(Unknown Source)         at java.awt.Component.dispatchEvent(Unknown Source)         at java.awt.LightweightDispatcher.retargetMouseEvent(Unknown Source)         at java.awt.LightweightDispatcher.processMouseEvent(Unknown Source)         at java.awt.LightweightDispatcher.dispatchEvent(Unknown Source)         at java.awt.Container.dispatchEventImpl(Unknown Source)         at java.awt.Window.dispatchEventImpl(Unknown Source)         at java.awt.Component.dispatchEvent(Unknown Source)         at java.awt.EventQueue.dispatchEventImpl(Unknown Source)         at java.awt.EventQueue.access$200(Unknown Source)         at java.awt.EventQueue$3.run(Unknown Source)         at java.awt.EventQueue$3.run(Unknown Source)         at java.security.AccessController.doPrivileged(Native Method)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.awt.EventQueue$4.run(Unknown Source)         at java.awt.EventQueue$4.run(Unknown Source)         at java.security.AccessController.doPrivileged(Native Method)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.awt.EventQueue.dispatchEvent(Unknown Source)         at java.awt.EventDispatchThread.pumpOneEventForFilters(Unknown Source)         at java.awt.EventDispatchThread.pumpEventsForFilter(Unknown Source)         at java.awt.EventDispatchThread.pumpEventsForHierarchy(Unknown Source)         at java.awt.EventDispatchThread.pumpEvents(Unknown Source)         at java.awt.EventDispatchThread.pumpEvents(Unknown Source)         at java.awt.EventDispatchThread.run(Unknown Source)      I don't have any idea whats wrong neing done. Please provide me the coorect way to export my application. Thanks in advance////


title: examples for Text3dElement

text: i am working in arcgis10. in java. can you all provide me with links to some sample program of `Text3dElement`. til now what i have implemented is given below                              Text3DElement td=new Text3DElement();                       td.setText("run away");                    IElement e=new Text3DElement();                    e=(IElement)td;                    e.setGeometry(point);                          iGlbCnt=new GlobeGraphicsLayer();                   iGlbCnt.addElement(e,0);      ` i was originally using `TextElement`, but it was not working properly on 3d globe. hence, i am trying to work with `Text3dElement` but unable to find any sort of example. thnks for ur help.


title: How to install PostGIS and pgrouting on Windows if Stack Builder fails to connect?

text: I know this has probably been covered and that there are plenty of guides out there, but for the life of me I cannot get pgrouting installed. All guides/how-tos are different or assume prior knowledge. I'm a QGIS newbie and please assume I know nothing! **How do install pgrouting? And what do I need first and how do I install whatever is needed? Can some helpful soul offer a step-by-step guide?** I've been trying to do this for days. I can get Postgres installed but the Stack Builder will not connect. I am thoroughly confused. I'm using Windows 7 Enterprise 64bit Service Pack 1, QGIS Lisboa. I'm now having trouble getting PostgreSQL to install.               Loading additional SQL modules...      Executing cscript //NoLogo "C:\Program Files (x86)\PostgreSQL\9.2\installer\server\loadmodules.vbs" "postgres" " _*_ *" "C:\Program Files (x86)\PostgreSQL\9.2" "C:\Program Files (x86)\PostgreSQL\9.2\data" 5432 Script exit code: 2 Script output: Installing the adminpack module in the postgres database... Executing 'C:\Users\ucesga1\AppData\Local\Temp\rad87BB2.bat'... psql: could not connect to server: Connection refused (0x0000274D/10061) Is the server running on host "localhost" (::1) and accepting TCP/IP connections on port 5432? could not connect to server: Connection refused (0x0000274D/10061) Is the server running on host "localhost" (127.0.0.1) and accepting TCP/IP connections on port 5432? Failed to install the 'adminpack' module in the 'postgres' database loadmodules.vbs ran to completion Script stderr: Program ended with an error exit code Error running cscript //NoLogo "C:\Program Files (x86)\PostgreSQL\9.2\installer\server\loadmodules.vbs" "postgres" " _*_ *" "C:\Program Files (x86)\PostgreSQL\9.2" "C:\Program Files (x86)\PostgreSQL\9.2\data" 5432 : Program ended with an error exit code ![enter image description here](http://i.stack.imgur.com/hCSC3.jpg)


title: WMS: point "pointer" symbol with offset - getFeatureInfo search area slight mismatch

text: I'd like to ask what ways there might be to adjust a point's getFeatureInfo search area so that it matches a "pointer" (i.e slightly offset) symbol. We are running GeoServer 2.2.0 (part of an Opengeo Suite) and the symbology(SLD) for some point layers is various forms of "pointers", that is to say that the center of the icon is not at the center of the point, rather it points to the point, like an arrow or a badge with a small pointer at the bottom or a map pin. (SVG graphics.) There appears to be two ways one can get the point symbols offset: I. Modify the SVG graphic to contain a transparent area in one end, see exhibit A. This method provides an empty area for getFeatureInfo corresponding to the transparent area in the SVG symbol as in geoserver's docs on point slds - one can however then click to the side of the visible symbol and still get a hit, which is unwanted. II. Use the offset function in the SLD for the point layer. The effect this has is that while the "pointer", or pin symbol is properly offset the search area remains centered on the point. See B. Many thanks for any suggestions.


title: Split all OSM roads within boundingbox in 20m segments and save to new table

text: I have a PostGIS Database with all streets from osm within a given boundingbox. Now I want to split all roads from my table planet_osm_roads into segments of max 20 meters length and save those segments to a new postgis table. The geometry of the roads is save in "way geometry(linestring,900913)". Can anyone give me a hint on how to accomplish this? Cheers, Martin


title: Arcpad Message Box and Keyboard Button Issue

text: I'm using ArcPad 8 along with ArcPad Studio to create a custom application for use on Windows Mobile 6 devices. There are data entry forms and on those forms are buttons that, when clicked, present the user some information in a message box. For a button's onclick event, I specify an event script (using JScript) such as:               Application.MessageBox("Some information for the user...");      That all works fine. However, when the user dismisses the message box, the button for toggling the input keyboard/transcriber/recognizer that is located in the lower right corner of the application window immediately disappears, which prevents the user from showing and hiding the keyboard as they attempt to enter data in other various controls on the form. Has anyone else ever had this problem and know of any solutions? Thanks!


title: Topological relationship between two polygon feature layers

text: I have two polygon feature layers: one representing land parcels and one representing roads. Some land parcels (polygons) coincide with a part of a road (another polygon feature). How can I check which land parcels coincide in any of their boundaries with a road? Of course help will be great if a code in VBA and ArcObjects is already available. I note also that I have considered the topology rules defined in http://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=Topology_rules and also and the ITopologyRule interface but I could not be sure which rule in both cases to use for solving my problem. Thanks. Demetris


title: Any freely available vector data in cylindrical projection?

text: I'm trying to develop an application in PHP to plot various points in Europe from a database, and serve them dynamically on a map. I think this should be relatively straightforward when using a map in the cylindrical projection - but I just cannot find one. I only require country borders, and pretty low level of detail. The smallest section of the map to be displayed at one time is likely to be the UK, and the largest output size is probably around 1280 x 1024px. Does anyone know of a freely available (or paid, if free isn't an option) vector? Many thanks, Jack


title: Arcmap: save all layers in Table of Contents to layer files

text: Using python, how can one save all layers in the Arcmap (v10) table of contents pane to layer files on disk?


title: Clearing map cache - ArcGIS Engine

text: Is there a programmatic way to clear the cache of a map in ArcGIS Engine 10.0? It is apparent that there is some caching scheme operating under the covers, and I'd like to clear all of its cached data. To be clear it isn't _leaking_ memory, I just don't want it to be caching anything. This is a fairly vanilla (*.mxd). Just a raster dataset, no basemap layers or anything fancy going on. ![graph of memory](http://i.stack.imgur.com/CFpby.png)


title: QGIS Mosaicing with null values in images

text: I have a set of landuse landcover data which needs to be mosaiced. Whenever i try using QGIS Raster>miscellenous>merge tool. giving no datavalue as '0', it returns either **(A)**                ' Traceback (most recent call last):      File "C:\PROGRA~2\QGISVA~1\bin\gdal_merge.py", line 509, in      sys.exit(main())      File "C:\PROGRA~2\QGISVA~1\bin\gdal_merge.py", line 495, in main      fi.copy_into( t_fh, band, band, nodata )      File "C:\PROGRA~2\QGISVA~1\bin\gdal_merge.py", line 241, in copy_into      nodata_arg )      File "C:\PROGRA~2\QGISVA~1\bin\gdal_merge.py", line 61, in raster_copy      nodata )      File "C:\PROGRA~2\QGISVA~1\bin\gdal_merge.py", line 101, in raster_copy_with_nodata      to_write = Numeric.choose( nodata_test, (data_src, data_dst) )      File "C:\PROGRA~2\QGISVA~1\apps\Python27\lib\site-packages\numpy\core\fromnumeric.py", line 296, in choose      return choose(choices, out=out, mode=mode)      MemoryError       or **(B)** ' The process crashed sometime after starting successfully' I am using QGIS 2.2 on Windows 8 PC. I have a set of 30 images in .tiff format which needs to be mosaiced. I also tried doing it with 2-3 images but the same problem persists.


title: Saving in arc ascii grid format

text: This might be a really simple question but I can't figure it out! I need to save a DEM in Arc ascii grid format so that I can load it in to Matlab. If I explain a bit about what I am doing maybe it will help: I have a 3D plane produced from the CFtool in matlab and I need to work out the volume between this plane and my DEM. I was going to import into matlab and try working it out. I have tried the other way around, i.e. exporting my plane from matlab but this doesn't seem to work. Anybody know how to save the DEM (in tif at the mo) as arc ascii grid format? Or have any other ideas on what to to work out my volume. Working in Arcmap10 and Matlab R2012a


title: Trace Downstream task: which ArcObject interface and method?

text: Does anyone know what ArcObjects interface and method the Trace Downstream trace task on Arcmap's Utility Network Analyst toolbar calls?


title: OpenLayers plugin fails to install

text: When I try to install OpenLayers plugin for QGIS 2.2.0 Valmiera I get the error message the plugin is broken; python said: 'QPyNULLVARIANT'object has no attribute '_ _getitem_ ' Does anyone have any ideas; i am a complete beginner, so I hope for a simple fix. Or, where else can I get base maps. I am using mac OS X 10.9.2 Thanks.


title: Problems installing r.roughness

text: I am trying to install r.roughness in GRASS 6.4.3. It starts compiling but at the end I get this error:               ERROR: Module did not install properly      I cannot identify the problem.   Can somebody that has come across the same error help me?


title: Identify corrupt shapefiles or rasters via Python or ArcObjects

text: I'm interested in a way to identify rasters (e.g., ESRI GRID format) or shapefiles that are corrupt that would prevent them from being used in an analysis tool. I have a large number in structured folder tree and would like to be able to script something that will walk the file tree and be able to ID rasters that, for example if you tried to view it in ArcCatalog it would not allow it to be displayed and accordingly would make any geoprocessing tool fail. Same deal for shapefiles or feature classes. My preference would be somehow to do this in Python (and I wouldn't be opposed to a non-ESRI python module or other) but I'm open for anything. I guess you could hack something together by trying to validate geometry or capitalizing on some aspect of the file design, but what would make the most sense? Or is there another way to do it that is specifically designed for such a purpose? My goal is to just be able to run a simple Quality Control mechanism on some large week-long preprocessing files to make sure they're all sound before I start doing the actual batch processing. thanks, Tom


title: Changing the line color of a marker box using OpenLayers

text:
I'm trying to figure out how to change the line color of a marker box in openlayers. I've been learning from the following example: ![enter image description here](http://i.stack.imgur.com/KYWHt.png)               <!DOCTYPE html>                            <meta http-equiv="Content-Type" content="text/html; charset=utf-8">             <meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=0">             <meta name="apple-mobile-web-app-capable" content="yes">             <title>OpenLayers Boxes Vector Example</title>             <link rel="stylesheet" href="http://openlayers.org/dev/theme/default/style.css" type="text/css">             <link rel="stylesheet" href="http://openlayers.org/dev/examples/style.css" type="text/css">
Boxes Example Vector

                 box, vector, annotation

                 Demonstrate marker and box type annotations on a map.

                         What I tried doing was replacing the `boxes` line with the following:                           //var boxes  = new OpenLayers.Layer.Vector( "Boxes" );                 var styleMapBoxes = new OpenLayers.StyleMap({borderColor: '#008000'});             var boxes  = new OpenLayers.Layer.Boxes( 'Boxes', {styleMap: styleMapBoxes});        but it did not work. I've read through the example here, but still do not understand how I can edit the box style. According to the documentation `borderColor` is a valid attribute.


title: Algorithm to match trajectory

text: I've data in the format (lattitude,longitude,timestamp). What I want to do is between timestamp1 and timestamp2 I want to find out which people might have taken the same routes. Is there some kind of trajectory matching algorithm? Let's say we have 100 users each has a set of points of (lattitude,longitude,timestamp) I want to know which set of users might have matching trajectory. Matching trajectory would be the same route taken, hence in a given set of timestamps the lattitude and longitude should be same or close enough as well as the timestamp should be same or close enough. Close enough can be about 30 seconds for timestamp while for space let it be like 200 metres. I can do this via a brute force approach and I'm looking for better solutions. Posted it on SO: http://stackoverflow.com/questions/16725427/design-an- algortihm-to-match-trajectory I hope I would get more thoughts here.


title: Is there an industry-standard official mapping of GLONASS satellites to "PRN" values?

text: In the U.S. GPS system, satellites are uniquely identified by the PRN number broadcast by each satellite. This ID is usually visible on most GPS receiver platforms that provide access to individual GPS satellite signal information to identify the broadcasting satellite (e.g., accessible on Android via the GPSStatus.Listener and the GpsSatellite.getPrn() method (http://goo.gl/QLse3)). However, the Russian GLONASS system does not currently have the concept of a PRN number broadcast by every satellite, as many of the legacy satellites are an FDMA system, not a CDMA system such as the U.S. So, in GLONASS these satellites broadcast the same code on a different frequency. However, an ID is still available for GLONASS satellites for many GPS receiver positioning platforms that support dual GPS/GLONASS positioning, often via the same access point as the PRN value for U.S. GPS satellites (e.g., on Android via the GpsSatellite.getPrn() method, with seemlingly standardized values across devices I've seen from 65 to 88. See the first screenshot listed in the GPSTest app on Google Play (https://play.google.com/store/apps/details?id=com.android.gpstest) for an example list of U.S. and GLONASS satellites obtained from the GPSStatus.Listener on a Samsung Galaxy S3). I haven't yet been able to find an official industry-standard source for how these ID values for GLONASS satellites are determined. However, proprietary documentation for particular products have shown GLONASS satellites numbered from 65 to 88 [1] and 38 to 61 [2]. Is there any industry-standard official mapping between GLONASS satellites and "PRN" values?


title: Unable to reference GPS waypoints correctly in ArcGIS

text: I'm working on a project. Waypoints were taken in the field with a Garmin Dakota 20. I've tried two different ways of putting them into ArcGIS. One way is through converting the GPX file to a shapefile via DNRGPS. This results in the coordinate system being wrong (even if changing it in the program) and the waypoints end up at the wrong place. The other is with the data interoperability tool which results in the waypoinst missing spatial reference. Even after adding such in through the catalogue this doesnt move the Waypoints anywhere. What am I doing wrong? I want them matched to the current data frame with its maps (Sweref99) **EDIT:** After editing once again through the ArcCatalog i tried changing it to Geographic Coordinate System instead of the projected (still Sweref99) and now, instead of them being nowhere, I get them at the same spot I did through dNRGPS (quite a bit off, but at least on the map).


title: ArcObjects Background Process (Thread) Progress Bar?

text: Using a toolbar, written in ArcObjects, it is possible to push processing to a separate thread so that the GUI does not lock. Is it possible to have the thread update a progress bar in the same position (lower right - arcmap) that a geoprocessing tool does? I must be missing some obvious documentation to gain access to that progress bar update location.


title: Why does install of add-on r.clump2 for GRASS 7 give error?

text: I am trying to install the add-on r.clump2 for grass7. What is wrong and how do I fix it?. Thanks The error code im receiving is the following:               > (Fri Jul  5 15:13:42 2013)                                                 > g.extension -s extension=r.clump2     > svnurl=http://svn.osgeo.org/grass/grass-addons/grass7 Fetching     > <r.clump2> from GRASS-Addons SVN (be patient)... Compiling...     > flag.c:1:23: fatal error: grass/gis.h: No such file or directory      > #include <grass/gis.h>     >                        ^ compilation terminated. make: *** [OBJ.x86_64-apple-darwin10.8.0/flag.o] Error 1


title: Parallel Computing and Raster Operations in ArcObjects gives cryptic error

text: This is a problem I have been struggling with for awhile. The code below takes 1-3 input feature layers (pFLayer) that were stored into a list called pFeatureLayerCollection, casts them into IGeoDatasets and executes EucDistanceFull in parallel to get distance rasters from them. The output IGeoDatasets are added to a list. After the Parallel.ForEach ends, if there is more than 1 output raster, I attempt to get the localminimum of the rasters using the method below. It runs fine if I give it one input raster. It also runs fine if I run it with multiple rasters sequentially (not in parallel) However, if I try to execute it with more than one raster in parallel, it fails with this error every time: **Attempted to read or write protected memory. This is often an indication that other memory is corrupt. (Esri.ArcGIS.SpatialAnalyst at ESRI.ArcGIS.SpatialAnalyst.RasterLocalOpClass.LocalStatistics(IGeoDataset LocalData, esriGeoAnalysisStatisticsEnum Type)**                            ESRI.ArcGIS.Geometry.IEnvelope pExt; //create a new envelope to hold the feature class' size information                 pExt = new EnvelopeClass(); //instantiate the envelope as an object                 pExt.XMin = **[Some given known double value]**;                 pExt.YMin = **[Some given known double value]**;                 pExt.XMax = **[Some given known double value]**;                 pExt.YMax = **[Some given known double value]**;                      string StrOutRasPath = **[Some given valid filepath]**;                      List<IGeoDataset> pMinimumGeodatasetCollection = new List<IGeoDataset>(); //store output distance rasters                           Parallel.ForEach(pFeatureLayerCollection, pFLayer =>                 {                     ESRI.ArcGIS.DataSourcesRaster.IRasterWorkspace pTempRWS;                     ESRI.ArcGIS.Geodatabase.IGeoDataset pOutputRaster = null;                     ESRI.ArcGIS.GeoAnalyst.IRasterAnalysisEnvironment pEnv = null; //Create a raster analysis environment. This sets up an area to do raster operations in                     ESRI.ArcGIS.SpatialAnalyst.IDistanceOp pDistanceOP; //Create raster DistanceOp object                     pDistanceOP = new RasterDistanceOpClass(); //define the newly created object as a new raster distance operation object                          string g_tempFolder = null;                     lock (strOutRasPath) //threadsafe creation of folders from method                     {                         g_tempFolder = CreateTempFolder(strOutRasPath);                     }                               lock (pWSF)                     {                         pTempRWS = (IRasterWorkspace)pWSF.OpenFromFile(g_tempFolder, 0);                     }                                         ESRI.ArcGIS.Geodatabase.IWorkspaceFactory pWSFact; //create general workspace factory (this is the directory you will be working in)                     //pWSFact = new ShapefileWorkspaceFactoryClass(); //set general workspace as a shapefile workspace. WARNING: USING THIS CAUSES A STRANGE COM REWRAPPING ERROR. IT IS AN ESRI BUG                          Type t = Type.GetTypeFromProgID("esriDataSourcesFile.ShapefileWorkspaceFactory"); //DO NOT REMOVE THIS LINE OR THE FOLLOWING TWO LINES. Fixes strange casting error from COM rewrapping. (ESRI BUGFIX)                     System.Object obj2 = Activator.CreateInstance(t);                          pWSFact = (ShapefileWorkspaceFactory)obj2;                          ESRI.ArcGIS.Geodatabase.IFeatureWorkspace pFeatureWorkspace; //create a feature workspace to handle features                     pFeatureWorkspace = (IFeatureWorkspace)pWSFact.OpenFromFile(g_tempFolder, 0); //set feature workspace as the shapefile workspace defined earlier, and specify a file location                                    ESRI.ArcGIS.Geodatabase.IGeoDataset pTempMinRaster;                                  ESRI.ArcGIS.Geodatabase.IGeoDataset pTempGeodataset;                     ESRI.ArcGIS.Carto.IFeatureLayer pTempFeatureLayer;                          ////////////                     //Loop through all features.                     ////////////                          ESRI.ArcGIS.GeoAnalyst.IConversionOp pConvOp;                     pConvOp = new RasterConversionOpClass();                          ESRI.ArcGIS.Carto.IFeatureSelection pFeatureSelection;                     pFeatureSelection = (IFeatureSelection)pFLayer; // ***IMPORTANT*** THIS SETS THE CURRENT FEATURE LAYER IN COMBOBOX TO THE FEATURE SELECTION TO BE USED FOR CONVERSION TO RASTER                      ESRI.ArcGIS.Geodatabase.ISelectionSet pSelectionSet;                     pSelectionSet = pFeatureSelection.SelectionSet; //Defines a selection set of features. The pFeatureLayer information is stored here.                          ESRI.ArcGIS.Geodatabase.IFeatureCursor pFCursor2;                     ESRI.ArcGIS.Geodatabase.IFeature pFeature2;                     ESRI.ArcGIS.Geodatabase.IQueryFilter pQueryFilter;                          pQueryFilter = new QueryFilterClass();                     pFCursor2 = pFLayer.Search(null, false);                     pSelectionSet = pFeatureSelection.SelectionSet;                     pFeature2 = pFCursor2.NextFeature();                          pEnv = (IRasterAnalysisEnvironment)pDistanceOP;                     pEnv.OutWorkspace = (IWorkspace)pTempRWS;                     pEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, cellsize);                          lock (pExt) //threadsafe grab of pExt from outside Parallel.Foreach                     {                         pEnv.SetExtent(esriRasterEnvSettingEnum.esriRasterEnvValue, pExt);                     }                                        ESRI.ArcGIS.Geodatabase.IGeoDataset temporaryGeodataset = (IGeoDataset)pFLayer.FeatureClass;                              pOutputRaster = pDistanceOP.EucDistanceFull(temporaryGeodataset, true, false, false);                                    lock (pMinimumGeodatasetCollection) //threadsafe adding to the list of output geodatasets                     {                         pMinimumGeodatasetCollection.Add(pOutputRaster);                     }                          pFeatureSelection.Clear();                      });                      //LIST CREATION COMPLETE. NOW IF TWO OR MORE RASTERS WERE CREATED, GET THE MINIMUM USING LOCALMINIMUM METHOD                      if (pMinimumGeodatasetCollection.Count == 1) //if just one input feature class                 {                     pMinimumRaster = pMinimumGeodatasetCollection[0];                     pFeatureLayer = pFeatureLayerCollection[0];                      }                      else if (pMinimumGeodatasetCollection.Count == 2) //if two input feature classes, get local minimum                 {                          pMinimumRaster = LocalMinimum2Rasters(pMinimumGeodatasetCollection[0], pMinimumGeodatasetCollection[1], strOutRasPath + "\\tmpVoro0004");                      }      There is the bulk of the code, and below is the implemented LocalMinimum method. LocalMinimum3Rasters does the exact same as LocalMinimum2Rasters, but acommodates for one more raster. The final line of this code is where the program is failing.                   public IGeoDataset LocalMinimum3Rasters(IGeoDataset pRaster1, IGeoDataset pRaster2, IGeoDataset pRaster3, string strOutPath)         {                  ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection pRBCollTmp;             ESRI.ArcGIS.DataSourcesRaster.IRasterBand pRBand;                  ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection pRBColl;             pRBColl = new RasterClass();             pRBCollTmp = (IRasterBandCollection)pRaster1;             pRBand = pRBCollTmp.Item(0);             pRBColl.AppendBand(pRBand);             pRBCollTmp = (IRasterBandCollection)pRaster2;             pRBand = pRBCollTmp.Item(0);             pRBColl.AppendBand(pRBand);             pRBCollTmp = (IRasterBandCollection)pRaster3;             pRBand = pRBCollTmp.Item(0);             pRBColl.AppendBand(pRBand);                  //Create a RasterLocalOp operator             ESRI.ArcGIS.SpatialAnalyst.ILocalOp pLocalOp = new RasterLocalOpClass();                  ESRI.ArcGIS.Geodatabase.IWorkspace pWs;             ESRI.ArcGIS.Geodatabase.IWorkspaceFactory pWSF = new RasterWorkspaceFactoryClass();             pWs = pWSF.OpenFromFile(strOutPath, 0);                  return pLocalOp.LocalStatistics((IGeoDataset)pRBColl, esriGeoAnalysisStatisticsEnum.esriGeoAnalysisStatsMinimum);              }      What could be causing this error? It seems to be directly tied to parallelism - I've noted that all operations work perfectly fine if I use a regular foreach instead of Parallel.ForEach. Extra information: The program creates the two distance rasters every time - If I go to the temporary folders I created in this code after the program crashes, I will find the two temporary rasters wholly intact and perfectly created.


title: How to fill selected feature(s) with transparent red fill using ArcGIS Server 10?

text: I'm writing a Server Object Extension (SOE) for ArcGIS Server 10, and I'd like to change the symbology of features selected from my feature layer such that the features are shaded with a transparent red fill. Here's the code I have so far. This code results in selected features being outlined in red, and filled with a solid red fill, where it is supposed to be a transparent red fill.                       featureClassToQuery = (IFeatureClass)mapServerDataAccess.GetDisplayDataSource(mapName, layerIndex);                  parcelSelectionQueryFilter = new QueryFilter()             {                 WhereClause = string.Format("APN='{0}'", apn)             };                  featureSelection = (IFeatureSelection)featureLayer;             featureSelection.SelectionSymbol = (ESRI.ArcGIS.Display.ISymbol)new ESRI.ArcGIS.Display.SimpleFillSymbol()             {                 Color = new ESRI.ArcGIS.Display.RgbColor() { Blue = 0, Green = 0, Red = 255, Transparency = 50 },                 Outline = (ESRI.ArcGIS.Display.ILineSymbol)new ESRI.ArcGIS.Display.SimpleLineSymbol()                 {                     Color = new ESRI.ArcGIS.Display.RgbColor() { Blue = 0, Green = 0, Red = 0 },                     Style = ESRI.ArcGIS.Display.esriSimpleLineStyle.esriSLSSolid,                     Width = 2.00                 },                 Style = ESRI.ArcGIS.Display.esriSimpleFillStyle.esriSFSSolid             };             featureSelection.SetSelectionSymbol = true;             featureSelection.SelectFeatures(parcelSelectionQueryFilter, esriSelectionResultEnum.esriSelectionResultNew, false);


title: What are the challenges to implementing a version control system for geospatial data?

text: Not that I am in any immediate need of a right answer here, but I've lately seen some efforts to introduce the concept of "(distributed) version control systems" for geographic data. Some examples (that I know of) are the three whitepapers from OpenGeo (1, 2 & 3) and the "Geosynkronisering (geosyncronization)" project by Norwegian GIS Software vendors and the Norwegian Mapping Agency. I've also found this question here, which mentions GeoGit (by OpenGeo), and this question about version control in ArcGIS. Being a developer I know (at least enough to be able to use them) how version control systems for source code (like SVN and Git) works, and my background in geomatics tells me that there are some unique challenges with geographical data that makes the approach not completely similar to the way source code (which is basically text) are handled. So, my question is thus: what are the challenges when dealing with (d)VCS'es for geographical data, how would you solve them, do we need them and are there other attempts to solve these issues than the ones I've mentioned? I know that the OpenGeo whitepapers will answer some of my questions, but what I'm really is after is a more "pedagogical" answer, in the style of "tell me like I'm a 10-year-old", so that I can refer people to a great explanation of the challenges and solutions that geographical data brings to the mix. I hope that someone with some insight will take time to provide some thoughts on the matter, as I said I'm not currently looking to solve a particular problem, but this topic is one that interests me. (Hope this is not a too open-ended question)


title: Python Addin toolbar- Environment mask parameter not being recognised?

text: I've created a button on a python addin toolbar which uses the selected features from the selected layer in the TOC, creates a feature layer from that selection and uses that feature layer to set the environment mask. This is so I can get the raster properties (max, mean) for just my area of interest from a larger raster. I had originally coded it as a tool/script in ArcToolbox and that tool works fine. Now after converting it to be used as addin toolbar button, the tool works but does not recognise the environment mask so gives me the max and mean for the whole slope raster. Can anyone see why this would be the case or offer any suggestions? (suggestions are welcome for more efficient code too :)               import os, arcpy     from arcpy import env     from arcpy.sa import *     import pythonaddins          arcpy.CheckOutExtension("Spatial")          class Terrain_Statistics(object):         """Implementation for FISToolBar_addin.Terrain_Statistics (Button)"""         def __init__(self):             self.enabled = True             self.checked = False         def onClick(self):             mxd = arcpy.mapping.MapDocument('current')             #df = arcpy.mapping.ListDataFrames(mxd, '*')[0]             TocLayer = pythonaddins.GetSelectedTOCLayerOrDataFrame()# Use the selected layer in the Table of Contents             outLocation = 'O:\\GeomasterPublicData\\General\\Terrain.mdb' # hardcoded for now             if not isinstance(TocLayer, arcpy.mapping.Layer): # if a layer in TOC hasn't been selected, give error                 pythonaddins.MessageBox('Please select polygons AND the polygon layer in the TOC', 'Error', 0)             else:                 count = 0                 # count the selection                 count = arcpy.GetCount_management(TocLayer)                 if count >= 1:                     cursor = arcpy.UpdateCursor(TocLayer)                     for row in cursor:                         forest = row.getValue("Forest")                     #forest = '903'                     pythonaddins.MessageBox('Selected ' + str(count) + ' features in ' + forest,'INFO', 0)                     arcpy.MakeFeatureLayer_management(TocLayer, "Selected")                     try:                         arcpy.env.mask = "Selected"                         slopeResult = arcpy.GetRasterProperties_management(outLocation +'\\Slope_' + forest, "MEAN")                         maxResult = arcpy.GetRasterProperties_management(outLocation +'\\Slope_' + forest, "MAXIMUM")                         slope = slopeResult.getOutput(0)                         maximum = maxResult.getOutput(0)                         pythonaddins.MessageBox("Selected " + TocLayer.name + ' polygons' + '\n'                                             + "Slope Model: Slope_" + forest + '\n'                                             + "------------------------------" + '\n'                                             + "Mean Slope: " + str(slope) + '\n'                                             + "Max Slope: " + str(maximum)+ '\n', 'Sale Area Terrain Information', 0)                         #arcpy.mapping.RemoveLayer(df, "Selected")                         except:                         pythonaddins.MessageBox("Something went wrong", 'Error',0)                         #arcpy.mapping.RemoveLayer(df, "Selected")                 elif count < 1:                     pythonaddins.MessageBox("No features selected. Please select polygons and try again", 'INFO', 0)


title: Change from layout view to data view using python? (ArcMap 10)

text: Is there a way to change from layout view to data view using python? I would prefer to do it using arcpy.


title: Proxy issues using WFS from Cadcorp GeognoSIS

text: I have web and WMS/WFS servers hosted on the same machine, although on different ports. I'm currently using OpenLayers to consume the WMS feeds which works well, but now I need to allow the user to extract info about objects so am trying to incorporate WFS into the OpenLayers setup. Being on different ports, I've set up proxy.asp to do the WFS requests:               <%     set objHttp = Server.CreateObject("Msxml2.ServerXMLHTTP")     strURL = Request("url")     objHttp.open "GET", strURL, False     objHttp.Send          If objHttp.status = 200 Then         Response.Expires = 90         Response.ContentType = "application/xml;charset=UTF-8"         Response.BinaryWrite objHttp.responseBody         set objHttp = Nothing     End If     %>      I've tested the proxy using this URL: http:///maps/proxy.asp?url=http%3A//:4326/OgcService/WFS.impl?service=wfs?request=getCapabilities This returns the GetCapabilities which is identical to that returned if I access it directly, but without the proxy. However, when using this WFS, via proxy, I can't get OpenLayers to display any content. At this stage it could well be my OpenLayers code, which I can post if necessary, but the same WFS URL via proxy doesn't work in MapInfo or Cadcorp MapModeller (it does with just the base URL). I suspect the proxy is the issue. Any ideas what the issue is likely to be? Cheers!


title: GDAL: Geotransform - rotated rasters (and ArcGis)

text: I am working on a GDAL driver for a proprietary grid-format. In some cases, the grids are rotated in order to get smaller file sizes. In my GDAL geotransform array, I do the following (showing a hardcoded 30 degrees angle):               rot_deg = 30;     rotation = rot_deg * deg2rad;     poDS->dGeoTransform[0] = originX;     poDS->dGeoTransform[1] = cos(rotation) * cellSizeX;     poDS->dGeoTransform[2] = -sin(rotation) * cellSizeX;     poDS->dGeoTransform[3] = originY;     poDS->dGeoTransform[4] = sin(rotation) * cellSizeY;     poDS->dGeoTransform[5] = cos(rotation) * cellSizeY;      where originX and originY is the _center_ of the corner cell. In addition, I got an identical file in TIFF format, with same min/max + cell values.I then apply an ArcGis geotransform to the .tiff, rotating it 30 degrees aswell. When adding both files as layers, the tiff is 100% correctly aligned to the pivot point, while my GDAL format is slightly displaced. The screenshot shows the overlapping grids, with only the corner pixel visible. The cyan point is the pivot point. The dark pixel on top is the (correct) tiff, while the grey pixel beneath is from my own driver. The misalignment is the same over the whole grid. Anybody got an idea what could be wrong ? ![enter image description here](http://i.stack.imgur.com/7HEvc.png) EDIT/update: In the ArcGis geotransform, it is possible to enter a pivot point, which is the case here (center of the corner cell). However, it seems like GDAL's geotransform is not expressive enough (?) to handle this, hence it uses the lower left corner.


title: How do I import OS raster tiles into QGIS at the correct scale?

text: I have some tiles of OS raster mapping and their associated TFW georeferencing files. The GR files are in the same folder as the raster data and have the same names, ie TQ21.tif/TQ21.tfw. In QGIS 1.8, I add a tile as a raster layer and choose the OSGB 1936/British National Grid CRS. When I load the next tile over, it loads in the correct place. Fine so far. However, the problem is that the scale is completely out - to view a tile I have the onscreen scale indicator at around 1:365000000! If I go to Raster > Miscellaneous > Information there is nothing in the Raster Info box, but the fact that the tiles load in their correct place suggests QGIS _is_ reading the GR files. Am I missing something?


title: How to clear python console in QGIS?

text: After running several hundred lines of commands (and getting their subsequent error statements...), it is nice to be able to clear out the console window and start clean. In the Python Console window in QGIS, I would like to clear it out. If I close it and then open it again, the buffer is still retained. If I 'select all' and delete the content, I don't get a prompt and can't do anything to get one back. Any suggestions on how to clear out the scroll back buffer?


title: How to replace arcpy.CalculateField_management with an update cursor?

text: I'd like to replace this expression               arcpy.CalculateField_management(bufferFile, "Shape_Area", "float('!shape.area@squaremeters!')", "PYTHON")      with an arcpy.da.updatecursor. The reason for this is that i wand to avoid the cryptic syntax of the command above to improve code readability. any ideas if this is possible? thanks!


title: Arcpy & ARCGIS10: Make netCDFRasterlayer, error attribute not found

text: I am trying to create a raster layer from a NetCDF File using Python. If I am using the Arcgis 10 Python window it works without issues. But executing the same script in Pythonxy gives the error: > NetCDF: Attribute not found What could be the problem? I'm just wondering because it works within the Arcgis working area but not outside. Here's my python code:               import arcpy          inNetCDFFile = 'C:/Users/Babs/Dropbox/Masterarbeit/mpihm_cnrmb1nat_precip_monthly_1991_2000.nc'          variable = 'Precip'          XDimension = 'lon'     YDimension = 'lat'          outRasterLayer = 'Precip_2011'     bandDimension = ''     dimensionValue = ''     valueselectionmethod = ''          arcpy.MakeNetCDFRasterLayer_md(inNetCDFFile, variable, XDimension, YDimension,     outRasterLayer, bandDimension, dimensionValue,valueselectionmethod)


title: Using Qgis API and python, how can I return latitude and longitude of a point?

text: I have a point layer, that I'd like to return the Longitude, latitude using the qgis api. I assume I need to use QgsPoint


title: OpenLayers Strategies

text: Can we use the strategies `OpenLayers.Strategy.Fixed` and `OpenLayers.Strategy.BBOX` simulataneously on same vector layer?


title: Where can I download public data set, especially shape file?

text: The places I know include:   * http://www.weogeo.com/   * http://geocommons.com/   * http://koordinates.com/ Any other suggestions?


title: how to find vertices for a geometry in OpenLayers

text: I am having a LineString, where I need to find the vertices logically and save those vertices. Can anyone help me out. Or else how below can help me.               OpenLayers.Geometry.LineString = OpenLayers.Class(OpenLayers.Geometry.Curve, { getVertices: function(nodes) {         var vertices;         if(nodes === true) {             vertices = [                 this.components[0],                 this.components[this.components.length-1]             ];         } else if (nodes === false) {             vertices = this.components.slice(1, this.components.length-1);         } else {             vertices = this.components.slice();         }         return vertices;     }                     }       );        Please help me. Thank you


title: Visually Merging ZIP polygons and Extracting by boundary from SQL Server 2008

text: i have a lot of shapefiles. For simplicity, lets say all the zipcodes of the USA. I've imported them into SQL SERVER 2008. I wish to merge lots of the shapes into one shape. I'm sorta struggling to do this with Sql Server 2008 and sql script. Are there any programs out there that can do this? eg. click on these shapes and merge into one? Also, can i cut shapes out? So lets imagine i have a whole pile of zipcodes together .. and then i draw my own shape .. can I cut out my shape from those zipcode leaving a nice whole. Once more I can sorta do this with SQL Script but I'm getting really tired of this and would really love it if this stuff is all visual. Any ideas, folks?


I am using Django 1.3, PostGIS 2.0 , Postgresql 9.0 and GeoDjango. I had a database with a table for places. This table had latitude and longitude for places in the US. I enabled this database with PostGIS. I added PointField to my model. Then, updated this table with               SELECT AddGeometryColumn('table', 'poly', 4326, 'POINT', 2)          CREATE INDEX "table_poly_id" ON "table" USING GIST ( "poly" )          UPDATE table     SET poly = ST_PointFromText('POINT(' || longitude ||' '|| latitude ||')', 4326)      Now, I want to make query to find places within 5 miles of given latitude and longitude of a user's location. So, I did this               ref_pnt = fromstr("POINT(%s %s)" % ("-87.627778", "41.881944"))     records = Table.objects.filter(poly__distance_lte = (ref_pnt,D(mi=5)))      But, I am getting an error               DatabaseError: Coordinate values are out of range [-180 -90, 180 90] for GEOGRAHY type      I also get an error when I do this               SELECT ST_PointFromText(poly) FROM table           ERROR:  parse error - invalid geometry     HINT:  "0101000020E61000007E" <-- parse error at position 20 within geometry     CONTEXT:  SQL function "st_pointfromtext" statement 1      What could be the problem?


title: How to edit my attribute table with VBA

text: before I start you should know that I am a bloody newcomer when it comes down to coding. I am sitting on some university tasks regarding VBA and ArcGIS: My next task is to create a code snippet that enables the user to edit values (rows) in the attribute table. my attribute table consists of points with basic information such as name, address etc. What I would love to do is load a kind of LIST of all rows and be able to click on the specific row I want to edit. Feel free to post easier ways to get there.. as I said I am extremely new to the world of coding ;-) Thanks everyone Bill


title: Area calculation: Albers Equal Area vs Pseudo Mercator

text: I have a question concerning area calculations in South-East Asia for West- East range that goes from India to Eastern Malaysia and a North-South range from Northern India to Indonesia (xMin,yMin 74.2203,-5.89636 : xMax,yMax 118.927,30.1183). I found two projections that might be suitable for this geographical region:   * Albers Equal Area (Conical) for S-Asia (s. http://spatialreference.org/ref/esri/102028/) and   * Pseudo/Web Mercator (EPSG:3857, s. http://spatialreference.org/ref/sr-org/7483/) However, when I compare the results of both projections, I get large differences in total area (I used QGIS field calculator for area calculation):   * 620,589.41 km^2 for Albers Equal and   * 689,941.23 km^2 for Pseudo Mercator That is a difference of ard. 11% Can anyone help me to understand the differences of these two projections and to know which one is more appropriate for my purpose? Thanks so much in advance


title: How to make corners of openlayers with OSM map rounded?

text:
I'm trying to apply style with border-radius:XXpx to my openlayers map:               map = new OpenLayers.Map(name)     map.addLayer(new OpenLayers.Layer.OSM())      But rounded corners are hidden by tiles of the map: screenshot Is there any way to resolve this? **UPD:** There are following styles applied to my div:

      styles:               #map {         border: solid 1px red;         border-radius: 10px;         overflow: hidden;     }      And overflow:hidden does not help **UPD 2** I can't make this code to work in Chrome (several versions)


title: How to add Custom data (like for Indoor floors) and find routes using open source software?

text: I am a newbie to GIS platforms so maybe whatever I am asking is very simple or obvious. I have read documentation/tutorials and as a last resort I am posting here, so any help is appreciated. I am not asking for code per se just point me to the right direction. I am trying to create an Indoor Map application, I have installed PostGIS and pgRouting on the server. Now, I am not sure what is the easy way to add nodes and ways to the database because in every tutorial/documentation that I came across, people are using OSM data. For me, there is no OSM data. Can somebody explain how can we insert node & ways to database so that pgRouting can work. Thanks.


title: How to insert geometries from shapefile into existing table?

text: I'm running into a problem trying to figure out how I should insert a GIS coverage polygon into an existing table. I'm primarily using QGIS and the PostGIS Manager plugin to handle it. After some trial and error, I believe I have my table set up correctly to accept the polygon object that I want to link to my existing table. Below is the schema, pruned down a bit -- the actual table has more traditional SQL columns that are not applicable to this. However, effectively, I want to enable a spatial query against each of the current rows by using a polygon -- i.e., query against the rows by determining if a point is contained within the polygon.               CREATE TABLE "distribution"     (       location character varying(5) NOT NULL,       unit smallint NOT NULL,       "unitName" character varying(15) NOT NULL,       coverage geometry,       CONSTRAINT "distribution_pkey" PRIMARY KEY (location, unit),       CONSTRAINT enforce_dims_coverage CHECK (st_ndims(coverage) = 2),       CONSTRAINT enforce_geotype_coverage CHECK (geometrytype(coverage) = 'POLYGON'::text OR coverage IS NULL),       CONSTRAINT enforce_srid_coverage CHECK (st_srid(coverage) = 4326)     )      Given the above table and a ESRI shapefile, how do I actually go about inserting the shapefile into the appropriate row's coverage field?


title: Open source engine that is similiar to ArcEngine

text: If I was to develop a mapping application to be a standalone windows app, is there an open source engine that is somewhat comparable to ArcEngine?


title: Focus ArcMap application Window VBA

text: Very simple question to ask, I haven't been able to figure it out up to now though... How would you focus your ArcMAP application window using VBA? I have a script that runs and selects polygons and zooms to them, but i would like arcmap to focus and come to the front at that point. Suggestions?               Dim pDoc As IDocument     Dim pApp As IApplication     Dim pAppROT As IAppROT          Set pAppROT = New AppROT          If pAppROT.Count > 0 Then     Dim p As Integer     Dim x As Integer     p = 0     x = 0     Do While p < pAppROT.Count     If TypeOf pAppROT.Item(p) Is IMxApplication And (pAppROT.Item(p).Document.Title Like "Report.mxd") Then          Set pDoc = pAppROT.Item(p).Document     Set pApp = pDoc.Parent     x = 1     pAppROT.Item(p).Visible = True          End Sub          Set focus = pAppROT.Item(p)          Exit Do     Else          End If          p = p + 1          Loop          End If          If x = 0 Then              Set pDoc = New MxDocument         Set pApp = pDoc.Parent              pApp.OpenDocument ("c:\Report.mxd")              pApp.Visible = True          End If'


title: merging block group data that uses 2000 geography with 2010 geography

text: I would like to combine data sets from 2010 which use 2000 geography with 2011 data that use 2010 census geography. The data is from the contiguous US and is either at block or block group level. I am trying to get a sense of:   1. How difficult that would be to do correctly   2. Where a simple merge using the block or BG coordinates or code would go wrong (such as Las Vegas or Phoenix). I've tried to locate this information on the Census homepage as well as this site, but I am not sure if I am using the best terminology in my search.


title: Is it possible to attain a classified raster with several classes in just one QGIS raster calculator operation?

text: I'm using qgis' raster calculator and wanted to make a new raster (oraster) from the values of an input raster (iraster). something like: if iraster>2 then oraster=1; if iraster >4 then oraster=2; if iraster>6 then oraster=3 all in one operation. Anyone can help?


title: How to set polygon opacity on openlayers?

text:
I have some prolems showing polygons in openlayer. I want to divide one city in areas. For the moment I have zona1 and zona2. I took lat and long from google maps so I think they are correct. My problem, is that polygons are draw in black. I want them to have a inner color but not black and not 100% opacity, I don't know how yo change that. I will explain a little my code because I think it is so strange my idea of showing more than one markers. I have all variables in arrays, so I can get any possition in the array to do operations related. So I create variables for both areas and after that I put them on an array. I have 2 arrays for points. One array keep lat and the other one keep long. So I have to take data from each one each iteration. Sorry if there are any comment in spanish, it is my language. I tried to change all I saw.                <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">                   <title></title>

                       I made some changes to the code. I correct one of my question so I delete it from post. SO my only question about it is how to set opacity.


title: Creating a report from features that have been selected

text: I'm in need of some guidance, I have created an application where the users does the following workflow   1. Enter Postcode   2. Zoom to Postcode   3. Create a buffer at specified distance   4. Spatial Query to find points within the buffer   5. MsgBox appears with the Organisation Name that are within the buffer.   6. Create report with the Organisation details. Number 6 is what I need help on, I want to create a pdf report with details from the attribute table, and a map of the buffer and points, any ideas on what are my options/examples of using Arcobjects in creating reports that look good? Cheers


title: Use arcpy.SearchCursor with raster

text: How can I efficiently subdivide one raster into multiple rasters using the arcpy.SearchCursor? For example if I have a slope and I was to subdivide into 3 rasters:               slope1 with values between 0 and 25degress          slope2 with values between 25 and 50degress          slope3 with values > 50degress      Here I should easily solve the problem with **arcpy.AddFieldDelimiters** If for example i have one categorical file and I also want to do the same thing? If the file I have has only one field with the identification of the category and it's a string.               Value    Count    Type         1    1023     Protected Area 2000         2    1000     Protected Area 2001         3    9555     Protected Area 2002         4    9500     Protected Area 2003         5    200     Construction Site         6    (...)    (...)      Should I use **arcpy.AddFieldDelimiters**? Furthermore, if I want to aggregate certain categories in only one raster, other groups of categories in other rasters? When using the **Search Cursor** we are doing the changes directly to the source file ?should I create a copy for each raster I want to produce?


title: First time using ArcGIS JavaScript API - where do I find the JavaScript API to include (.js) into my solution?

text: I'm writing a web service that will interface with ArcGIS. One of the requirements is to put:               esri.config.defaults.io.proxyUrl = "proxy.ashx";      in the JavaScript portion of the page. Of course, it says "esri is undefined" if I don't include the library that defines "esri". Question: Where do I get this library at? I've searched all over the internet and couldn't find anything! Thanks!


title: Arcgis 10.1 server, arcgis manager shows page cannot be displayed

text: I am facing this issue only in one machine. I installed ArcGIS 10.1 SP1 in Windows 2008 R2 machine, installation went fine. But the select Arcgis manager, localhost:6080/arcgis/manager, it shows only page cannot be displayed. I tried these items, a) uninstall arcgis 10.1 sp1, installed arcgis server 10.1, still the same problem b) uninstall arcgis server products, cleaned the registry, renamed the c:/program files/arcgis/server,... keycodes ... ..commonfile/arcgis... rebooted the server installed arcgis 10.1 sp1, still the same problem c) i verified with windows support team, they confirmed no ports are blocked d) when i compared with good arcgis server in other windows 2008R2, there is a javaw.exe in the process, but this is not in this current machine. Any one have any idea/suggestions/work around to resolve this ? thanks, vadivelan. p


title: What are the differences between Append and load data in ArcGIS?

text: I confused about differences between Append tool in Arctoolbox and Load data in Geodatabase. I think both of them are the same in the result feature class.Both of them use a default output feature class and compare the target feature geometry and schema with another feature but there are differences in some aspects such as performance that I don't know .Another question is that where can i use Append instead of Load data ? Thanks for advance


title: How to store selected attributes of vector layer using PyQGIS

text: I'm following this PyQgis workshop: http://www.qgisworkshop.org/html/workshop/python_in_qgis_tutorial2.html#accessing- data-attributes I have a point layer shapefile. I want to be able to retrieve and store x, y and name of each point( they already exist in attribute table and have indexes of 10, 11 and 1), hopefully in a dictionary format ((x,y) as key, name as value). My problem starts here:               while provider.nextFeature(feat):     map = feat.attributeMap()        I can't figure out how to use this 'map' dictionary. iterating through it I either get "PyQt4.QtCore.QVariant object at 0x23C9C3B0" or "built-in method toString of QVariant object at 0x10D38D50". any help is appreciated!


title: Error when calculate field using VBscript or Python (slice a field)

text: I have an expression in VBscript that works well for labeling:               Function FindLabel (field)      if (INT(RIGHT(field, 2)) <> 1) then        FindLabel = INT(MID(field, 8, 3)) &"-"& INT(RIGHT(field,2))     else         FindLabel = INT(MID(field, 8, 3))     end if     End Function      This expression labels a text in this format: **VCC065-00102** to: **1-2** or from **VCC065-00501** to: **5** (using just the last 5 carachters, if the last one is > 1, it shows the number after a "-"). But, when I tried to use the same script for **calculating a field** , I've got a syntax error error: **Code block:** the same as above **Expression:** > FindLabel ([DATA]) Also, I tried to convert this script to Python: **Code Block:**               def FindLabel (Field):       Var1 = Field[-2:]       Var2 = Field[8:-3]       if (int(Var1)  <> 1 ):         FindLabel = int(Var2) &"-"&  int(Var1)       else:         FindLabel = int(Var2)      **Expression:**               FindLabel (!DATA!)      But it returns an error: _The field is not nullable._ (I don't know why, but the expression is returning a null value) Thanks!


title: Setting a renderer that uses a function for a dynamic layer in ArcGIS JavaScript API

text: How can I, if possible, set a renderer for an `ArcGISDynamicMapServiceLayer` which uses a function to determine feature values instead of just attribute values? When using the following code, the function specified never gets called. Instead, the value to be added last is used for every feature, which kind of seems like a bug to me. That is, the following, which I expected to render features using the default gray symbol only, instead renders all features red.               var map, mapService;     require(["esri/map", "esri/layers/ArcGISDynamicMapServiceLayer", "esri/layers/LayerDrawingOptions", "esri/renderers/UniqueValueRenderer", "esri/symbols/SimpleFillSymbol", "dojo/_base/Color", "dojo/domReady!"],         function (Map, ArcGISDynamicMapServiceLayer, LayerDrawingOptions, UniqueValueRenderer, SimpleFillSymbol, Color) {              var drawingOptions, optionsArray = [];              map = new Map("mapDiv", {             basemap: "osm",             center: [0, 0],             zoom: 5         });              mapService = new ArcGISDynamicMapServiceLayer(             "http://x.x.x.x:6080/arcgis/rest/services/Service_Name/MapServer",             {id: "test"});         map.addLayers([mapService]);              // The renderer displays polygons with a GRAY fill per default.         var defaultSymbol = new SimpleFillSymbol().setColor(new Color("#999999"));         // There are only two specified categories, using YELLOW and RED fill.         var yellowSymbol = new SimpleFillSymbol().setColor(new Color("#ffff00"));         var redSymbol = new SimpleFillSymbol().setColor(new Color("#ff0000"));              var renderer = new UniqueValueRenderer(defaultSymbol, function (feature) {             console.log("!");             return "other";         });              renderer.addValue({value: "yellow",             symbol: yellowSymbol,             label: "Category One (Yellow)", description: ""});         renderer.addValue({value: "red",             symbol: redSymbol,             label: "Category Two (Red)", description: ""});              drawingOptions = new LayerDrawingOptions();         drawingOptions.renderer = renderer;         optionsArray[0] = drawingOptions;         map.getLayer("test").setLayerDrawingOptions(optionsArray);     });      Following this sample, using a `FeatureLayer` instead and adding the renderer by calling `layer.setRenderer(renderer);`. produces the expected result but is not the answer to this question.


title: Where can I find a python GDAL 1.8.0 package?

text: I have python 2.7.4 installed on my RHEL 5 machine along with gdal-devel-1.8.0 and gdal-1.8.0. I tried to install the GDAL 1.9.1 package available on python's page (https://pypi.python.org/pypi/GDAL/1.9.1) but got errors that are related to trying to install v1.9.1 when my gdal libraries are 1.8.0 (see Error installing GDAL 1.9.0 python package). I can't find any other versions of python GDAL bindings. Can anyone help me please?


title: Where to get UK city boundaries

text: Currently I am working on a project where I have to draw boundaries of specific region of UK. Exactly like this![enter image description here](http://i.stack.imgur.com/ujAYS.png) I want to know where to retrieve this kind of GIS data. I have tried openstreetmap data for uk, in that some of the data are unavailable for UK cities. Any suggestion would be highly grateful for me. Thank you in advance.


title: Why can't I see any joined string fields for categorization?

text: After the join I can see all my data in the attributes of the card but when I want to select statistics in the tab style to categorize, I don't see them. I see that the real type variables and not those of type string. Are there any parameters to be changed? * * * original: après avoir fait la jointure je vois bien l'ensemble de mes données dans les attributs de la carte mais quand je veux sélectionner les statistiques dans l'onglet style pour catégoriser, je ne les vois pas. Je ne vois que les variables de type real et pas celles de type string. Y a t-il un paramètre à modifier?


title: Huge point cloud laser data in PostGIS - Storing and Processing it

text: I wonder, how it's possible to store huge sets of laser scanned point cloud data in PostGIS, with the time-aspect for processing it in mind. I know, there exists a geometry-object `Point` in PostGIS. But as far as I know it saves each point in a new tupel, which can make searching for any certain point a very slow process, if a few millions or more of them are stored. I found a paper from HSR Universtiy of Applied Sciences Rapperswill, discussing this topic. It suggests three ways to store such data: `Whole data in one tupel`, `Each point in one tupel` or `Splitting Data into Blocks` which are referenced by info-tables, holding the extends of each block. As the third way seems the most useful for locating stored points, I wonder if anyone already has made some experiences with it? The paper can be found here: http://wiki.hsr.ch/Datenbanken/files/pgsql_point_cloud.pdf Last but not least, I stumpled across a project on github, which seems to deal with point cloud manners in PostgeSQL. Unfortunately not much information about it around the net. So the same question here: Has someone already made some experiences with it? Is it usable for such purposes? Project can be found here: https://github.com/pramsey/pointcloud I would also be glad to hear about other suggestions, ideas or experiences, if there are any. But i must admit, that non-commercial solutions are prefered.


title: running geoserver in tomcat7 with eclipse

text: i downloaded geoserver,put it in tomcat7, lunched tomcat7 with powercmd and published my layers and styles.i created a web application, where i add the laters above to a google basemap. Now when i start tomcat7 with eclipse to run my web, in the console i don't see geoserver deployed on it.   Any idea how to solve this problem. thank you


title: What program is used to create this?

text: What program was used to create this type of mapbook? So if I want to get my roads to look similar to the above PDF - how would I go about doing that. ![enter image description here](http://i.stack.imgur.com/EKAWP.png)


title: Zip Code Boundary Maps for Commercial Use?

text: I need images of cities with features, streets, zip code names and zip code boundaries. For example, take this screen shot: ![enter image description here](http://i.stack.imgur.com/AaPAk.jpg) This is very close to what I need, except not every zip code name is listed in each boundary. If I could get city images like this, with complete zip code names and that allows commercial use, that would be perfect. I don't need zooming or tiling. Static images like the above are perfectly fine. I don't even mind capturing the images myself if there's a service that will let me get this data, provided I can use the images commercially. Do you know of any sources? I googled, but what I found was either lacking some features (as in the above) or did not allow commercial use.


title: How to display several points of a KML with the same coordinates in OpenLayers?

text: I've got a KML, where several points have the same coordinate, similar to following code:               <?xml version="1.0" encoding="utf-8"?>     <kml xmlns="http://www.opengis.net/kml/2.2">       <Document>         <name>test.kml</name>         <visibility>0</visibility>         <open>1</open>         <Placemark>           <name>Place 1</name>           <Point>             <coordinates>7.0901548880586,50.723973666434,0</coordinates>           </Point>         </Placemark>         <Placemark>           <name>Place 2</name>           <Point>             <coordinates>7.0901548880586,50.723973666434,0</coordinates>           </Point>         </Placemark>       </Document>     </kml>      If I click on the POI in an OpenLayers-map, only the last point can be selected. Google Maps' solution to the problem are `next`\- and `previous`-buttons, in Google Earth the placemark "splits" into two (or more) placemarks so you can select one of the placemarks. Is it possible to create a similar behaviour in OpenLayers? There is an analog problem if there are two KMLs with placemarks with the same coordinate. I need a solution for this problem, too.               <?xml version="1.0" encoding="utf-8"?>     <kml xmlns="http://www.opengis.net/kml/2.2">       <Document>         <name>test1.kml</name>         <visibility>0</visibility>         <open>1</open>         <Placemark>           <name>Place 1</name>           <Point>             <coordinates>7.0901548880586,50.723973666434,0</coordinates>           </Point>         </Placemark>       </Document>     </kml>          <?xml version="1.0" encoding="utf-8"?>     <kml xmlns="http://www.opengis.net/kml/2.2">       <Document>         <name>test2.kml</name>         <visibility>0</visibility>         <open>1</open>         <Placemark>           <name>Place 2</name>           <Point>             <coordinates>7.0901548880586,50.723973666434,0</coordinates>           </Point>         </Placemark>       </Document>     </kml>


title: How can I create profile cross-sections in ArcGIS 10?

text: Because the default 'Profile Graph' tool was not sufficient, I wrote a custom VBA script that creates scaled and exaggerated profile cross-sections in ArcGIS but I need more features:   * Visualize more than 1 surface   * Visualize vertical and directional boreholes   * Ability to symbolize / edit the output   * Linked with underlying data - changes in terrain (or profile 3D line) reflected in output A few years ago I tested an extension by A-Prime software called 'CrossView For ArcGIS' but I am not sure if it is compatible with ArcGIS 10 (Contacted A-Prime but they have not yet responded) http://www.aprimesoftware.com/ Area there any alternatives? Other 3rd party software? A functioning script? that can do this?


title: Zonal Statistics as Table Tool

text: I have a slope dataset. I have a polygon shapefile with territory within the extent of the slope map. I have to calculate the real surface area of this shapefile using the zonal statistics as table tool, which I have done and selected the Statistic Type "All" for results. Having these results, how do I know what the raster cell in square meter if for this extent? (ArcInfo 10.1)


title: Is there a way to define the relative position of items in print composer?

text: Is there any chance to make the items in print composer relative to the margins. Let's say for example the legend. I normally align it to the corner on the right bottom where the legend is 10 px away from the boarders of the map. If I add a new item to the legend, the legend gets bigger and I have to align it again. To give the legend a relative position independently from the legends size would be very favourable (tell qgis to always align the legend in the corner on the right bottom with a distance of 10 px always). Than I could define one layout and make several maps without readjusting those items all over again...


title: Configure Google Maps to display travel time by bike?

text: I move around by bicycle but didn't find how to configure Google Maps so that it displays travel time by bike: "x min ride home" would be more useful to me than "x min drive home" :-) Does GM supports this, and if yes, how to configure it accordingly? Thank you.


title: geoserver wrong tile generation on leaflet

text: I am creating world map with bathymetry using Geoserver 2.5.1 and leaflet. It works perfactly upto certain zoom level. but after that i get few tiles drawn in wrong color as shows in images. Does anyone have any idea why geoserver behaving like it or any direction to fix it? As i zoom in more on the map tiles are missing colors as shown below. ![enter image description here](http://i.stack.imgur.com/VHTPK.png) ![enter image description here](http://i.stack.imgur.com/sPdPZ.png) ![enter image description here](http://i.stack.imgur.com/bfcNg.png) ![enter image description here](http://i.stack.imgur.com/y556I.png) ![enter image description here](http://i.stack.imgur.com/h3PB6.png)


title: Lengths of lines within individual but overlapping buffers?

text: I have a set of points for which I have created a 12 mile buffer. Eventually, I would like to know how many kilometers of road exist within each buffer zone individually? The problem I see going forward is that the buffers have overlap. I need each buffer to be individually selectable so that I can apply my roads layer and obtain the desired information.


title: Microsoft visual studio run time library error in Qgis 1.8

text: When i open Qgis project, I get an error from Microsoft Visual C++ Runtime Library: Runtime Error! For other machine its working fine.Since I have uninstalled and re-installed the Qgis,but still same problem.What can I do for fixing this? OS: Windows XP SP2 Thank you for your help.


title: CSWcataloguesource.js plugin "add to map" button not working!

text: I got stuck in using CSWcataloguesource plugin in my OpenGeoSuite custom application. It shows all the layers available in geonetwork but when click the "add to map" button it shows the following error. TypeError: c.getLayer is not a function and m not able to add it as a overlay. Does anyone have idea what this error?and does anyone have used this CSWcataloguesource plugin? Thanks in advance.


title: Open street map accuracy problem

text: I want to locate shops with open street map. For that I send a query to nominatim, a webservice which geocodes coordinates from addresses. I draw a marker at the returned coordinates, and it´s inaccurate about 100 m. (~300 ft). To double check, I search the coordinates on google maps, and it is also the same (wrong) position, so the error seems to be at nominatim. Now I type the address I gave to nominatim and search with osm and google maps, where as osm returns the same inaccurate position, google maps returns the correct. Is it possible that OSM and Nominatim back on the same data in the background (like TIGER?). I'm located in a small town in Germany, is it possible that these areas are still too inaccurate? Can somebody tell me whats wrong there, or if I can solve this?


title: Free sources of mapping data for France

text: We are moving to France soon, and are currently investigating a lot of areas. Does anyone know of good online/offline resources for French mapping? I am very interested in mapping, weather, environmental, geological, hydrological data, as we are preparing to create a zero emissions home out there and want to get the best area. Thanks in advance


title: Schema Lock Test in Model Builder for ArcGIS 10.0 Desktop

text: I have made a model (Model Builder) that exports a Polygon Layer from a Parcel Fabric and Exports it as a set of Features Classes into a set of Feature Datasets. However, there are often other people working within the datasets. I would like to be able to check the dataset for a schema lock before trying to save the data, and then if possible keep checking until the schema lock is removed and the model can complete the process. Is this possible within model builder, and how?


title: Converting .kml to different projection in qgis not working

text: I have run into a slight problem it would seem. Yesterday I downloaded qgis in the hope that it could convert a .kml overlay of some administrative boundaries from EPSG:3857 (Google maps) to EPSG:27700 (Ordnance Survey). What I did: I loaded the .kml file into qgis as a vector file and then right click, save as and type in the new Projection I would like as in this photo: ![](http://i.stack.imgur.com/J0rGb.png) However, when I load the newly created .kml in google earth, it is exactly the same. I also converted both the original and the new .kml file to a .svg using an online tool and compared them. They too were exactly the same. Am I missing something here? Thank you ever so much in advance for any help!


title: Openlayers - how can i show scale values on zoom control

text: I remember seeing once zoom control in OL that had the zoom values on the scale control. I'm looking for such solution now to let the user quick pick desired scale. Any help will be appreciated


title: convert vector with multiple attributes to multiple rasters

text: I have a vector polygon feature with multiple attributes/fields. Ultimately, I want to convert it to a raster 'stack,' whereby each of the attributes/fields will be a different band of the raster feature. What is the best way to do this in as few steps as possible? platforms: QGIS 1.8 lisboa.


title: Assigning centerline address ranges

text: I've been given the task of assigning address ranges to a centerline feature class. I have all the data and I am just going street by street. My question: When assigning ranges, should I just enter the values I find for present addresses or would it be better to extrapolate the range over of the segment based on present values? I am attempting to have further geocoding results be somewhere close to where they should be. From the example in the graphic, where the x's are empty space and not unknown addresses:   * Section A would clearly be from 1 to 3.   * Section B: from 5 to 6 or from 4 to 6?   * Section C: from 12 to 13 or from 7 to 13? There won't be too many geocoding projects where the results will need to be spot on, but there are long stretches of road with scattered addresses (section C) and areas with clusters (section B) ![address example](http://i.stack.imgur.com/NKFsz.jpg) Another question, should the ranges be continuous or will breaks in the range and null values (stretch or side of the road with no addresses) overly effect the geocoding process? Thank you For the need of real example of real data: ![enter image description here](http://i.stack.imgur.com/r6TkG.jpg) The read dots are on the street in blue and range from 3115 - 3141 from left to right.


title: Missing Scottish Administrative Data in OS Opendata sets

text: I'm not sure if this forum is really the correct place to ask, but impertinence would lead me to believe that there may be some one on here that may know the answer.   I'm using the Ordnance Survey Open Data bdline_gb shape files to try to pull together some geographic coverage data for my organisation that covers England Scotland and Wales.   I'm working mainly with the Parish data set; because that's how ridiculously specific my organisation can be (I am determinedly prising them away from **part** parishes!). And fortunately England and Wales administrative Boundaries "build-up" from parishes (or communities in Wales). In the bdline_gb data set there is also very nicely the District, Borough, County and Unitary Authority shape files as well; so I can overlay to see where I have coterminous boundaries that no longer need to be down to Parish level. Unfortunately Scotland doesn't [still] have Parishes; but additionally the data sets don't seem to include any ADMINISTRATIVE boundaries for Scotland, at least not in it's proper coded sense; just the filler polygons without the appropriate codes. I can get electoral data till it's coming out of every orifice just not Administrative-   does anyone know where I can source this?


title: How to create a folder in HOME Folder with ModelBuilder in ArcGis 10?

text: I am just getting started with ModelBuilder and I would like to create a folder in the same folder where the respective mxd is located. The tool have to be usable in different projects, so setting Current Workspace to a specific folder doesn't do the trick. Is there any chance to get access to a path of the Home folder which would be the same in any given project?


title: OGR, writing TAB files, and BOUNDS

text: According to this: http://www.gdal.org/ogr/drv_mitab.html the TAB file writer sets default bounds, with the caveat that 'we should fix this'. Bounds causes issues with respect to coordinate precision, see here: http://www.dse.vic.gov.au/property-titles-and-maps/maps-imagery-and- data/data/vicmap/datums-and-coordinate-systems/victorian-bounding-box I would like to know if anyone knows how I can set bounds using OGR, and specifically, with GeoServer's WFS export using the OGR driver.


title: Is it possible to batch clip rasters in Quantum GIS?

text: I have been able to batch clip rasters in one directory using python in ArcGIS. However, I am not sure how to do the same in Quantum using the python console. Is this possible? If there are different ways to accomplish this, I'm open to suggestions. Thank you! MH


title: How do I stop the map at just one world/globe?

text: I am working on a project to put google maps layer and adding our data layers on top of it. using openlayers and google maps layer with out layer of top of it now the google layer has max resolution of 156543.0339 and 22 zoom levels, I can limit the zoom levels and max resolution etc but how do I stop the map at just one world/globe i.e. I don't want it to scrolling off to another wrapped around world! ![enter image description here](http://i.stack.imgur.com/PLm3b.png) Please see this fiddle for example, http://jsfiddle.net/eaJHn/5 Anyone !!! ?


title: Why is Orthorectify button not active in Image Analysis Windows in ArcGIS10?

text: ![Ortho button not active](http://i.stack.imgur.com/dJcCW.png) I wanted to Ortho rectify GeoEye HR imagery using existing RPC and IFSAR DTM in Image Analysis. According to ArcGIS Resource center http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/using_the_Orthorectify_button/009t000000mv000000/ the button is automatically active if the imagery has RPC file. In my case the ortho button was not active although my Imagery has RPC file. Please advice thank you.


title: Does CartoDB support alphabets other than Latin?

text: Does anyone know if CartoDB supports alphabets other than Latin? I'm trying to upload some data with Cyrillic characters but have had no success. Appreciate any feedback you can provide!


title: How to Select by Attributes on features sharing common value?

text: I am attempting to 'select by attributes' a point feature with the lowest numeric value which shares a common value with other point features in another field. For example: I may have 3 point features with address values: 123, 125 and 127. Each of these point features shares a parcel address value of 131 in another field. I simply need to select the lowest address point for each parcel address so that I can use a definition query to display a subset of the data (Only points with these lowest values). Thank you in advance for the help.


title: ERROR 6: EPSG PCS/GCS code error when running QgsGeometryAnalyzer().buffer

text: I am trying to buffer points from a FileGDB in a PyQGIS script. I setup my environment with the following batch file and run the python script below. The problem is that I get ERROR 6: EPSG PCS/GCS code 102184 not found in EPSG support files when QgsGeometryAnalyzer().buffer executes. The EPSG should actually be 3400 for the input points. How can I get this to work? NOTE: I can enter the commands in the QGIS python console and it works, but the standalone script does not. setup.bat               SET OSGEO4W_ROOT=C:\OSGeo4W     SET QGISNAME=qgis     SET QGIS=%OSGEO4W_ROOT%\apps\%QGISNAME%     SET QGIS_PREFIX_PATH=%QGIS%          SET PATH=%OSGEO4W_ROOT%\bin;%QGIS%\bin;%PATH%     SET PYTHONHOME=%OSGEO4W_ROOT%\apps\Python27     SET PYTHONPATH=%QGIS%\python;%PYTHONPATH%     SET GDAL_DRIVER_PATH=%OSGEO4W_ROOT%\bin\gdalplugins     SET GDAL_DATA=%OSGEO4W_ROOT%\share\gdal      buffer.py               import os, sys     from qgis.core import *     from qgis.analysis import QgsGeometryAnalyzer          outputPath = r"D:\buffer"     gdb = os.path.join(outputPath, "data.gdb")          size = 500          QgsApplication.setPrefixPath("C:\\OSGeo4W\\apps\\qgis",True)     QgsApplication.initQgis()     tpts = gdb + "|layername=tpts"     tbuff = outputPath + "\\buff.shp"     tptsQgsLyr = QgsVectorLayer(tpts, "tptsQgsLyr", "ogr")     QgsGeometryAnalyzer().buffer(tptsQgsLyr, tbuff, size,  False, True, -1)     QgsApplication.exitQgis()


title: How to quick get all symbols for feature layer with ArcObjects

text: I have a `FeatureLayer` and I want to enumerate all the symbols used in it. Simple way - it is to call `SymbolByFeature` for each features. Something approximately likes:               Dim l as IFeatureLayer     Dim c as IFeatureCursor     Dim r as IFeatureRender     Dim f as IFeature     Dim s as ISymbol          Set r = IGeoFeatureLayer(l).Renderer     Set c = l.Search(0, True)     Set f = c.NextFeature          Do Until f is Nothing       Set s = r.SymbolByFeature(f)       [...]  ' check whether this symbol is in collection and add it if necessary       Set f = c.NextFeature     Loop      The problem with this method is that it too long, especially when there is a lot of features in the feature layer.   Question: Is there a faster way to get a collection of all the symbols in a feature layer?


title: How do I get the best print quality from QGIS (and is Arc better)?

text: QGIS exceeds my expectations by a country mile. I am impressed. But; I work with people who use Arc and they have found the best resolution output images are made by creating PDFs. They tell me that .png and the rest out of Arc are not satisfactory and they don't bother with them. Typically we are making the output appear in a report and so the PDFs of the maps and words are assembled in a PDF editor package to make the final document. As a result I have been toying around with QGIS on my Win7 machine and notice that the best output comes by printing directly out of the map composer. This output on my Docuprint C2120 printer is as good as any of the Arc PDFs and is truly impressive. QGIS makes very attractive maps in this configuration. However I don't seem to be able to create PDFs or .pngs and then print them which are anywhere near as good. I find the .svg output to Inkscape unrealiable where for instance the grid/graticule may come out oversized. By way of example, the good quality output is where a sub-6 pt font can be read easily off the printed page. Obviously quality linework in cartography is the holy grail.... Any thoughts? ....Also I am using QGIS trunk...


title: Return results to user of python script from ArcToolbox : EDIT 9.3

text: I have a script in a custom toolbox which sums the values in a given field. After the user indicates which field (currently I have them typing it in but would love to know how to get the ArcToolbox to generate a dropdown after an input is identified) the script returns the value in a win32gui box. Using win32gui means I have to install win32com package on the users machines. Is there another easier solution to return the value of sum-ing the values in a field..preferably if the dialogue box lets you copy the text with the mouse? Thanks. EDIT: I am at 9.3..thanks.


title: Cannot load postgis layer from sql string in pyqgis

text: I'm trying to load a layer using the following code in qgis python console:               sql = "(select the_geom, ID from mytable)"     uri = QgsDataSourceURI()     uri.setConnection("localhost","5432","mydb","postgres","XXXX")     uri.setDataSource("myschema",sql,"the_geom","","ID")     vlayer = QgsVectorLayer(uri.uri(),"LayerName","postgres")     QgsMapLayerRegistry.instance().addMapLayer(vlayer)      It doesn't add the layer but also doesn't return any error. I'm using QGIS 2.4.0-Chugiak


title: How to convert these mapping coordinates

text: I have an older database at work and I've been asked to map the recorded GPS coordinates for each record. Unfortunately, I can't figure out the map coordinates that were recorded. As an example, 1301 Calle Oriente in Milpitas California shows as: 06158991 and 01987998 If I look at google maps, 1301 Calle Oriente is 37.446687, -121.886884 Is there a website that I convert 06158991 and 01987998 to a map like google maps? Is this a valid mapping coordinate format? If so, please let me know what format this is. TIA


title: What could I do to make QGIS find all the plugins?

text: I'm starting to use Ubuntu 12.10 and would like to get my QGIS running on it. Installation was fine but when I tried to add a package called Points to Paths the Plugin Manager was unabled to find it. Indeed there appeared to be only four plugins available on the two repositories that where activated by default (The QGIS Oficial Repository and the orphans repository). So after activating third part repositories there seem to be more plugins but the desired one doesnt appear. Is there a way to install these plugins manually on the console? Or what could I do to make QGIS find all the plugins? Has anybody expirienced similar problems?