title: Developing Carto Mobile App

text: I have the project to develop a mobile application (Android) being able to get offline maps, go on the field, collect datas about the area we're at and then synchronize it with our data server when we're back to the office. This is for my end of study work (sry, don't know a best translation in english), and i looked at apps like gvSIG Mini, qGIS for Android, qmap, .. But none of this really fits my need, as it is hardly adaptable .. Any ideas on api/existing apps/anything ?


title: How long for surveying quality GPS to get accurate 3D position information

text: I recently read an Ordnance Survey Blog article that seemed to suggest that it was necessary to leave the GPS device in use for four hours to gather requisite data. I am aware that "domestic" GPS units benefit from being left in one position to get better positional accuracy but had assumed that survey quality equipment was "accurate" immediately. What am I missing please? Nige


title: Using FME to merge lines in multiple SHP layers

text: I have multiple layers of complex laine data that I need to flatten int a single layer of lines. The issue is that some of the line data overlaps and intersects and is not end node to end node. How in FME can I take this data anf flatten it into a single set of joined lines?


title: How to use a database in Qgis

text: What is the primary use of database in Q-Gis? How to connect to a database in Qgis?


title: Unable to create a circular polygon with ArcGIS Explorer SDK

text: I am having trouble creating a circle with a specified radius of 200 feet. I have tried several ways to translate the points along the circular arc. here is my current code:               var display = ESRI.ArcGISExplorer.Application.Application.ActiveMapDisplay;          var center = display.TrackPoint();     var convertedRadius = Unit.Convert(200, Unit.Feet, center.CoordinateSystem.Unit);         var firstPoint = GeometryOperations.Move(center, convertedRadius, 0);          var points = new List<Points>();     points.Add(firstPoint);          for (int i = 1; i < 360; i++)     {         var point = GeometryOperations.Rotate(firstPoint, center, (Math.PI / 180) * i) as Point;         points.Add(point);     }          var poly = new Polygon(points, display.CurrentCoordinateSystem);     var circleGraphic = new Graphic(poly, Symbol.Fill.Solid.Red);     display.Graphics.Add(circleGraphic);      I only seem to have a problem creating a polygon. Using the above loop I can place 360 markers in a circle around the center point. Alternately, If there is a way to create a buffer from a point using programmatic methods that would suffice. I am trying to query all shapes that are within 200 feet of a target point.


title: Add print/export button to axToolbarControl

text: I have a map that's being rendered inside a custom application, and now I need to add to the ESRI.ArcGIS.AxToolbarControl an export (to PNG) button and a print button, but honestly don't know how. I've already looked at this and this samples, but couldn't make them work.


title: Where is DatabaseSupport directory in the ArcGIS client?

text: May I know where is this folder stated in the guide? DatabaseSupport directory in the ArcGIS client installation location and move it to a location on the Oracle server. http://resources.arcgis.com/en/help/main/10.1/index.html#//019v0000000s000000


title: convert shapefiles to cad files

text: hi i'm having a difficulty in converting shapefiles to cad file. i'm currently using the ArcGIS 9.3.1 version and i tried to convert it using the toolbox>conversion tool>to cad>export to cad....the conversion was partly successful but then it doesn't have any attributes when i tried to open it in ArcMap. When i try to open it in Autocad i can't see a thing... how can i resolve this? thank you.


title: How convert a genuine UK shapefile(other than the one from natural earth) to d3/topojson?

text: I am following Mike's page (http://bost.ocks.org/mike/map/) to make UK's map. I've converted Mr. Crone's UK postcode shapefile into topojson but failed to display it correctly. http://www.sharegeo.ac.uk/handle/10672/51?show=full               ogr2ogr -f GeoJSON postcodes.json gb_pca_2010.shp     topojson -o uk.json --id-property Label postcodes.json      The resulting json is of 2.2M big. I have used the following code to display it               var height, svg, width;     width = 960;     height = 1160;          svg = d3.select("body").append("svg").attr("width", width).attr("height", height);          d3.json("uk.json", function(error, uk) {       var projection, subunits;       if (error) {         return console.error(error);       }       subunits = topojson.feature(uk, uk.objects.postcodes);       projection = d3.geo.mercator().translate([width / 2, height / 2]);       return svg.append("path").datum(subunits).attr("d", d3.geo.path().projection(projection));     });      The result is that firefox freezes and in the end, I get half black and half white display. :( I have also tried mapshaper to simplify it to 1MB but still not working. What else I could have tried?


title: 10.1 beta - register dll-s

text: I just tried to register my extension for ArcGIS 10 in the 10.1 beta and it failed when I tried to register the dll: "Could not load file or assembly 'ESRI.ArcGIS.ADF.Local' Version=10.0.0.0." Do I have to re-compile it against 10.1 dll-s in order to make it work? For the 9.x version all I had to do was to compile them against the base level (9.2) and it worked in 9.3 as well.


title: Examples of geovisualizations of global connectivity

text: I stumbled upon interesting visualization of global connectivity patterns: > 'Major road and rail networks in Africa, along with transmission line and > underwater cable data.' ![enter image description here](http://i.stack.imgur.com/eQEPh.jpg) (Source, more pictures in the gallery. Original.) Can you point me to other examples of visual representations of connectivity in global scale? Connectivity in this case applies to both transportation and communication. And geovisualizations might include static maps as well as interactive web applications.


title: ArcObjects: How to Speedup Spatial Query to Feature Layer?

text: I have ROAD (POLYLINE) Shape File with more than 400,000 Records. I am working on a code which goes like this:        1. Loop through all Feature of the Shape File   2. Take First Feature   3. Find the LAST VERTEX of that feature. I get a POINT   4. Query to FEATURE Layer all the POLYLINES Connected to that POINT and get back result in feature cusros   5. Do some Process on returned feature of RESULT Cursor   6. Go to Step1 until Loop is not completed. When I do STEP-4, it obviously take lot of time because, I create a BUFFER of Say 1 meter on that POINT and Make a spatial FILTER and Run LAYER.SEARCH QUERY on entire featurelayer using that SPATIAL FILTER. Everytime querying 400,000 records making it very slow. Is there a better way to do this thing ?   * I have to use ArcObjects VBA.    * I do not want to Create NETWORK from the layer.


title: Preparing CSV files for use in ArcGIS for Desktop?

text: I wonder if there are any recomendations how to properly prepare CSV files for use in ArcGIS for Desktop. I ask because I have some troubles using CSV files because ArcGIS attributes wrong field types to my columns and also misinterprets special characters such as á or ê. I have read in the ESRI forum that there is a so-called schema.ini file that defines somehow the field types e.g "Col22=V002 Text" see here http://forums.esri.com/Thread.asp?c=93&f=1149&t=64464 That's kind of funny because I have often seen these .ini files on my disc but never actually wondered what they are good for. It is kind of weird that Excel stores such metadata in an extra file since other programs like R don't do so. I already tried to manipulate this .ini file with little success since I didn't find out how to apply for example "string" type. There are some information on MS sites, see here: http://msdn.microsoft.com/en- us/library/windows/desktop/ms709353%28v=vs.85%29.aspx but I couldn't find a solution. Also I didn't really like the idea to work with this .ini file because it is quite a bunch of work to define and type all the fieldnames when I have let's say 50 columns. And these .ini files might get lost, etc. **So, is there any better way to properly prepare CSV files for use in ArcGIS for Desktop?**


title: repeating records in attribute table

text: I have a point shapefile containing records of bird sightings with location, species and land cover type that has about 30,000 records. I noticed today that one record is repeated about 60 times, with the same gps points, same location name, same everything. But when I highlight all of the repeating records, they are actually all different records. I can't seem to get the attribute table to accurately reflect the records; I've tried copying the file, renaming it, exporting the incorrect records to a new shapefile, restarting Arcmap multiple times, etc. Is the file corrupted? What can I do to fix these records?


title: Leaflet with EPSG4326 Tile Layer (using maptiler.org)

text: I have successfully been using Klokan Petr Pridal's MapTiler (http://maptiler.org ) to take GeoTiff files and create TMS tile pyramids for Leaflet (option #1 which projects to EPSG3857, also called Google Maps compatible Spherical Mercator). But my source data is in EPSG 4326, so I am now using option #3 for WGS84 Plate Caree (Geodetic). I believe it is using gdal2tiles under the hood. The issue I have is that I do not see any of these overlay tiles in Leaflet. Even though am being careful to make sure the map is created with the right CRS. I have tried tuning TMS: false, and TMS: true. No luck. One thing I see is that the numbering even for the top levels z/x/y is weird: (e.g. 0/1/2.png ). Should a Krakow, Poland map really be at coordinates: 5/35/24? or 10/1136/795? With the Spherical Mercator I have 10/568/346!!! Also, all the specs for TMS tile numbering that I have seen are only good for EPSG3857 (spherical Mercator) which I believe has bounds of _85/-85 lat. But EPSG4326 would go to 90. So what gives in the numbering?               var map;     var krakowLayer = null;     var mapOptions = {         center: new L.LatLng(0, 0),         zoom: 1,         tms: true,         crs: L.CRS.EPSG4326     };         map = new L.Map('map', mapOptions);         krakowLayer = new L.TileLayer("./tiles4326/{z}/{x}/{y}.png", { tms: false, opacity: 0.8 });          function addLayers() {         var baseLayers = {             "CloudMade": cloudLayer,         };         var overlays = {             "Krakow Layer": krakowLayer         };         L.control.layers(baseLayers, overlays).addTo(map);     }


title: Using osm2pgsql data with pgrouting 2.0

text: I have successfully manage to make the workshop from of pgrouting (wrote for v1) with pgrouting v2. Now I want to use pgrouting 2 on my postgis database created with osm2pgsql. The problem is first to create the segments and this blog post was very helpfull for this. However,I get:               gis=# select pgr_createTopology('network', 0.01, 'the_geom', 'gid');     ERROR:  We can not determine the srid of "network.the_geom"!      Is there a way to set the srid properly so that `pgr_createTopology` finds it. Or is there an alternative way? I can't believe osm2pgrouting is the only importer for pgrouting. Thanks


title: A couple GeoExplorer questions

text:
I just started working with the OpenGeo Suite. I'd really like to use it for a County web map, but have a couple questions... 1) Is there a way to add a lat/lon or other X/Y location to the URL and have the map center and zoom to the point? Something like "?lat=46.265&lon=-90.206&zoom=7" at the end of the URL. 2) Is the Query tool available in a published map? It's not in the Choose Tools window when publishing. Thanks. === UPDATE === I will drop part 2 and stick with part 1. I guess I forgot to mention the crucial bit that I'm working in GeoExplorer. I did have the tag, but don't fault anybody for missing it. The bounding box isn't working in GeoExplorer. I have the map published in an iframe. Here's a screenshot of the test page without the bbox... i.imgur.com/Wc022Ee.jpg. Here's a screenshot with the bbox... i.imgur.com/7MES4cn.jpg. That's just a blank map that hides behind the real map and only shows up when the real map doesn't. I do have a correct bbox... i.imgur.com/MbAOgVD.jpg in the code. (I guess I need a higher rep before I can add more direct links, so you'll need to copy paste the image URLs if you want to see them.) Here's the source behind the page...
Iron County Test Map

This is a test map, version T-001.  A disclaimer might go here.
                  I tried putting the "?bbox=" at the end of the "src=.../1". Thanks for any help. And thanks for the welcome.


title: How can I convert an Integer field to a Real field in a shapefile?

text: How can I convert an Integer field to a Real field in a shapefile ? I'm searching a tool, or a way to do that (may be under QGis ?!). Regards


title: Server side clustering - cluster (feature) geometry

text: I am trying to implement server side clustering and show it using `OpenLayers`. I am using this good material as an example. As a clustering result, a get some list of features, some of them are clusters. Like this (it is an output of previous example):               Array (          [0] => Array (              [0] => Array ( [id] => marker_3 [lat] => 59.431602 [lon] => 24.757563 )              [1] => Array ( [id] => marker_4 [lat] => 59.437843 [lon] => 24.765759 )              [2] => Array ( [id] => marker_6 [lat] => 59.434776 [lon] => 24.756681 )          )          [1] => Array ( [id] => marker_5 [lat] => 59.439644 [lon] => 24.779041 )          [2] => Array ( [id] => marker_2 [lat] => 59.432365 [lon] => 24.742992 )          [3] => Array ( [id] => marker_1 [lat] => 59.441193 [lon] => 24.729494 )      )      where 0 element is cluster containing 3 features. But how can I compute this clusters geometry? Or how can I visualize it using `OpenLayers`?


title: getting the specified geographic transformation, or at least the GCS, of a data frame

text: First, the TLDR digest:   * I'm getting a **bogus value** from `dataFrame.geographicTransformations` on some machines; how do I get the correct value as shown in the properties dialog? and/ **or**   * Given a spatial reference with a PCS, how do I find out (in Python) **the GCS** or datum that it's based on? And now the details... (Reposted from http://forums.arcgis.com/threads/76061-getting-the-specified- geographic-transformation-or-at-least-the-GCS-of-a-data-frame) I'm writing an arcpy script in ArcGIS Desktop 10.0. This script opens a user- specified MXD and processes user-specified layers in it. Part of the processing involves projecting the data from whatever GCS it's in, to whatever PCS the data frame has in its spatial reference. Sometimes, the input GCS and the output PCS are in the same datum (e.g. WGS_1984), and so I don't need to supply a geographic transformation to the `Project_management()` tool. Other times, they are not in the same datum; e.g. the data frame coordinate system might be Europe_Lambert_Conformal_Conic, which uses the datum ED_1950. If I don't supply the right geographic transformation (or lack thereof), the Project tool fails. How do I know what transformation to use? In the case where the datums are different, I can require that the data frame in the MXD file have a transformation method specified in its properties, e.g. ED_1950_To_WGS_1984_10. If that was the end of the story, my job would be simple: if the data frame has a transformation method in its properties, pass that to `Project_management()`, otherwise omit that parameter. The script would just look at the value of `dataFrame.geographicTransformations`. The problem is that on some machines, when I ask for `mydataframe.geographicTransformations`, **I get a spurious value of "NAD_1927_To_NAD_1983_NADCON".** So if the data frame has no transformation methods specified in its properties dialog, `dataFrame.geographicTransformations` returns `["NAD_1927_To_NAD_1983_NADCON"]`. Which causes `Project_management()` to fail, if that transformation doesn't apply to the GCS and PCS actually being used. For example, when the data's GCS is GCS_WGS_1984, and the data frame's PCS is Asia_Lambert_Conformal_Conic, no geo transformation is needed (both have datum D_WGS_1984), and none is specified. So if I pass `dataFrame.geographicTransformations[0]` to `Project_management()` as the geographic transformation parameter, I will cause it to fail because the NAD transformation is inapplicable. If the data frame _does_ have a transformation method specified in its properties dialog, then (on these same machines) `dataFrame.geographicTransformations` returns a list with "NAD_1927_To_NAD_1983_NADCON" as the first element, followed by the one specified in the properties dialogue. Again, this is on some machines. On others, the spurious "NAD_1927_To_NAD_1983_NADCON" doesn't appear, so I can rely on `dataFrame.geographicTransformations` to be set, or not, as appropriate. So my question is, where is this spurious transformation coming from, and how do I get rid of it, so that if I don't need a transformation, I can avoid passing a bad one to Project_management()? There's nothing in the Environment Settings under Output Coordinates / Geographic Transformations. Or, if I can't get rid of it ... maybe I can detect when the data's GCS and the data frame's PCS have the same datum, so I don't need a geographic transformation. In that case I would just not pass any such parameter to `Project_management()`. The trouble with that approach is, I don't know how to get the datum, or the GCS, of a PCS spatial reference. I can look at `dataframe.spatialReference.PCSCode`, but `spatialReference.GCSCode`, `GCSname`, `datumCode` and `datumName` are blank (see http://help.arcgis.com/en/arcgisdesk...00000p6000000/ \- these properties are only available with a GCS). I know that **ArcGIS desktop knows** the GCS being used by the data frame, because when I open data frame properties, the Coordinate System tab tells me not only the PCS but also the GCS (including datum). But the **GCS is not available** thru `dataframe.spatialReference`, that I can tell. So, if someone could tell me   * how to accurately find the specified geographic transformation of the data frame (i.e. get rid of the uninvited NAD geographic transformation) (this would be the preferred solution); and/or   * how to determine the GCS or datum of a PCS (can my script look it up in some standard place, based on the PCS code?) ... I would be very appreciative.


title: How to import a simulation from ParaView to ArcMap?

text: I made a simulation and reworked the output with **ParaView**. I simulated an avalanche on a DEM and what I want to do is to compare the area of the simulated avalanche with that of the real avalanche. To do this I would use ArcMap 10 where I already have the real area as a shapefile. How can I bring what I have on ParaView to ArcMap? Thank you very much, Michele


title: Get road type (clazz) in osm2po

text: Can I get clazz of a segment in Java? For example,               RoutingResultSegment rrs = graph.lookupSegment(path[j])      How can I get `rrs`'s clazz?


title: How to extract zipped shapefile and check its shape type using arcpy module without writing on disk

text: I have a zipped shapefile saved on my disk. I want to extract it in memory and check the shape type. I'm using arcpy module and want to stick to only default modules comes with python installation. I tried something like ` import cStringIO, zipfile, arcpy stringIO = cStringIO.StringIO() file_ = r"C:\Polyline.zip" zipFile = zipfile.ZipFile(file_) for name in zipFile.namelist(): if name.endswith(".shp"): stringIO.writelines(zipFile.read(name)) print arcpy.Describe(stringIO).shapeType ` But it is not working. I got the trace back: ` Traceback (most recent call last): File "", line 8, in File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy_ _init_ _.py", line 1190, in Describe return gp.describe(value) File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\geoprocessing_base.py", line 374, in describe self._gp.Describe(*gp_fixargs(args, True))) RuntimeError: Object: Describe input value is not valid type ` How can I unzip a shapefile in memory and finds its shape type?


title: Does map export using GarminCustomMap plugin for QGIS work on Mac OS?

text: When I try to export a map using GarminCustomMap (1.4) plugin QGIS 2.01 on Mac OS 10.9 quits with no messages. I suspect this might be a Maverick related issue. Is anyone successfully using the plugin with 10.9? I have reported the crashes on the plugin issues page.


title: Opening InMemory Table to execute some SQL gives Runtime Error 445 - Object doesn't support this action?

text: I have a InMemory that gets created, what I then want to do is run INSERT, UPDATE and DELETE queries on the table. The code i have at the moment is below, I create a new InMemory table, add it to the map (which works).                   Dim pInMemTable As ITable         'create Temp DP table         Set pWSF = New InMemoryWorkspaceFactory         Dim PName As esriSystem.IName         Set PName = pWSF.Create("fred", "InMemory", Nothing, 0)         ' Debug.Print ">>" + PName.NameString + "<<"         Set pFWS = PName.Open              Dim pFlds As IFieldsEdit         Set pFlds = New esriGeoDatabase.Fields         pFlds.AddField MakeField("OBJECTID", esriFieldTypeOID)         pFlds.AddField MakeField("DP_REF", esriFieldTypeString, 10)         pFlds.AddField MakeField("WANTED", esriFieldTypeString, 1)              Dim pUID As New UID         pUID.Value = "esriGeodatabase.Object"              Dim ptable As ITable         Set ptable = pFWS.CreateTable("DP_Temp", pFlds, pUID, Nothing, "")        ' Debug.Print ptable.              Dim pSAT As IStandaloneTable         Set pSAT = New StandaloneTable         Set pSAT.Table = ptable         pSAT.Name = "DP_Temp999"              Dim pMxDoc As IMxDocument         Set pMxDoc = ThisDocument              Dim pSATC As IStandaloneTableCollection         Set pSATC = pMxDoc.FocusMap              pSATC.AddStandaloneTable pSAT         pMxDoc.UpdateContents              Set pInMemTable = pSAT              Dim path As String         path = "InMemory/DP_Temp999"         If pWSF Is Nothing Then             MsgBox "pwsf is nothing"         Else            ' MsgBox "pwsf is something"             'MsgBox pWSF.WorkspaceDescription(False)                  'pwsf.Open(             Set pWS = pWSF.OpenFromFile(path, 0) ' **********PWS doesnt get set***********             If pWS Is Nothing Then                 MsgBox "PWS IS NOTHING"              Else                 MsgBox ("Yes!!")             End If              End If      but the code fails where PWS is getting set as it always returns as Nothing. I have tried various ways of setting the path value e.g. "InMemory/DP_Temp999", "In_Memory/DP_Temp" and trying the "\" as opposed to "/" but it still comes back as nothing. Does anyone have any ideas please? or is there a better way of executing SQL on a table that will allow SELECT, DELETE and INSERTS. ******UPDATE**** Right I have been working on this for a few days have have been able to set pWS, by casting by ptable to a idataset object and then setting pWS to the pdataset propertyset properties which works. However I have then tried the following code, to see if SQL can be executed                   If pWS Is Nothing Then             MsgBox "PWS IS NOTHING"         Else             'MsgBox ("Yes!!")             Dim pWorkspaceProperties As IWorkspaceProperties             Set pWorkspaceProperties = pWS                  Dim pWSProperty As IWorkspaceProperty             Set pWSProperty = pWorkspaceProperties.Property(esriWorkspacePropertyGroup, esriWorkspacePropCanExecuteSQL)             If pWSProperty.IsSupported = True Then                 MsgBox "SQL is supported"                 Dim sql As String                 'sql = "select * from 'DP_TEMP';"                 sql = "UPDATE dp_temp set dp_temp.[DP_REF] = 'Test';"                 pWS.ExecuteSQL sql             Else                 MsgBox "SQL is not supported"             End If         End If      When the code is run, the pWSProperty.IsSupported returns as True, however when the .ExecuteSQL method is run, I get the following error Runtime Error 445 - Object doesn't support this action. So I am a bit confused as to whether SQL can be run or not.


title: QGIS 1.8 compiled against PyQt 4.8+?

text: I'm currently running QGIS on Win7 64bit. My installation says it was compiled against and running against PyQt 4.7.1, although according to the Osgeo4W installer and PyQt4.QtCore.PYQT_VERSION_STR, I actually have PyQt 4.8.3 installed under the Osgeo4w directory. I'm running into a bug in my PyQt code (related to QNetworkAccessManager) that, according to the bug tracking log, was fixed in PyQt 4.8.0. It seems, then, that the PyQt 4.7.1 module is being used, and the 4.8.3 module is being ignored. I'd like to resolve this issue. I've seen a few posts on various forums that suggest that the generic distribution of QGIS 1.8 is compiled against PyQt 4.8. I've tried reinstalling QGIS using the Osgeo4w installer after removing all references to PyQt from my system path, but I still end up with QGIS compiled and running against PyQt 4.7.1. Could anyone shed some light on this issue for me? Should QGIS 1.8 installed through Osgeo indeed be a distribution compiled against PyQt 4.8+? That is, am I doing something wrong in the installation process? Or will I need to compile QGIS manually if I want it to be compiled against PyQt 4.8? Thanks very much. Alex


title: Creating a vector based design from a raster database

text: I have a raster table that I've created using Postgis raster. It has 365 records and the table looks like this:               +-----+------------+---------+------------+     | rid |  rast      |filename |date_of data|     +-----+------------+---------+------------+     | 100 |            | file56  |2002-01-09  |     | 101 |            | file57  |2002-01-10  |     +-----+------------+---------+------------+      The second column _rast_ is raster files in WKT format and data type is _raster_. All my raster files have the same extents and cellsizes. I want to turn this raster design into a vector design like the following:               +-----+---------------+----------+     | ID  |pixel_centroid |array_366 |     +-----+---------------+----------+     | 1200|               |{26,20,..}|     | 1201|               |{28,68...}|        +-----+---------------+----------+      Imagine the raster files have 1600 pixels and accordingly I will have 1600 records in my new table. 'Pixel_centroid' is the point geometry of the pixels and could be calculated using               select ST_centroid((ST_PixelAsPolygons(rast, 1)).geom)  from raster_table      and array_366 is pixel values of all the records in the order of date_of_data:               select (ST_PixelAsPolygons(rast, 1)).val  from raster_table      might be used but I don't know how to put all the values in their place and populate the array.


title: Google Earth plug-in: how to get name/description of a model, from lat/long or click

text: In Google Earth, if you move the mouse over a 3D model (building), it turns light blue. If you click on it, you get a description (or sometimes just an approximate address). Google Earth plug-in web apps like Monster Milktruck or First Person Cam are wonderful educational tools for exploring an area, and learning about it spatially. For example you are walking or driving around Rome, and you see buildings like the Coliseum or the Forum, but you don't know what they are; so it would be very useful to be able to click on them as you can in Google Earth, and get a description. Is there any way to do this in the G Earth Plug-in API? I see that a hitTest can tell you whether you clicked on a model, and if so, what the lat/lon of your hit is. But surprisingly, it doesn't tell us _which_ model was hit. Nor can I find any way to look up the model in Google 3D Warehouse using the lat/lon. Can someone tell me a way to do this? Or is it for sure not allowed?


title: Does Bing maps silverlight API work on Windows phone 7?

text: Has anyone seen any good resources on using Bing Maps Silverlight APIs on WP7?


title: Can I get the metadata about the themes in Oracle MapViewer?

text: I am trying to use the MapViewer xml api to generate a dynamic theme (for generic highlighting support), but don't want to expose table names and such to users, I would like then to communicate in terms of themes. Is there some way that I can obtain the mappings from theme to table so that I can generate dynamic themes based on `jdbc_query` element? The `jdbc_query` needs me to give it some SQL which identifies the instances to be included in the dynamic theme. I want the user to be able to specify a theme name and then be able to construct an sql statement which gets the geometry of the nearest neighbour to a particular point on that theme. To do this I need to know which table the data in the theme came from and which column contains the geometry in that table. Can I find that information using an `info_request` perhaps?


title: bicycle / bike map road grade markers - semi-automated?

text: So I'm making a bike commuting map for our area, and one of the features to show (common on a lot of bike maps) is some type of road grade marker. We're planning to show three categories: 'kinda steep, pretty steep, and crazy steep' - we're not planning to use those names - but you get the drift. (Actually, four categories - anything less than the first category just doesn't get a marker.) The markers are one, two, or three arrows - here's a clip: ![enter image description here](http://i.stack.imgur.com/0an1D.png) The problem: since we live in the foothills - not much flat stuff around here - the task of getting all the grade markers right is pretty daunting and tedious. So - it would be nice to automate it. So I've got the DEM for the area from the NED, and got a nice slope raster from it (thanks to folks that helped figure that out). Overlaying that on the bike map (the roads are a vector layer of course), you can see a way to conceptually automate this task. Basically, making a vector layer that is the boolean AND of the roads layer and the slope layer of certain values (whether by categorizing the slope layer, or by making separate slope layers per grade range ahead of time) should get it done. I'd want the output to be a vector layer. But, I'm not quite clear on how to do all that. Suggestions on workflow? I realize this concept doesn't account for direction (which way is uphill) so I'll plan on doing that by hand (flipping the vector selections as needed, or, doing so with an attribute that is used by the style to decide which way to draw the arrows) UPDATE: this will be done in QGIS 2. A thought: make a vector layer (of shapes) out of raster categories, i.e. a vector shape of areas with slope between 5 and 10 degrees (this would cover one of the categories of steepness), then do a boolean AND of that layer with the roads layer; so now you have a bunch of disjoint lines; that's exactly what is needed; just style these lines with the double-arrow (or as appropriate for that category), and flip direction of each segment as needed (which is easy to do from memory or from a map if needed).


title: Force a collection of unknown number of points to be a polygon

text: I am working with point clusters and want to create the rough shape of the cluster in the form of a polygon using `ST_ConvexHull()`. I'm having problems with single outliers though, as they don't convert to a polygon but to a `POINT`:               SELECT ST_AsText(ST_ConvexHull(ST_Collect(point))) ...          POLYGON((13.01741532 52.40118295,13.382100866593 52.516280103729, ...))     POINT(12.5 51.5)      Is there a way to force small geometries (single or double points) into a polygon? I want to simplify the subsequent processing by guaranteeing polygons; accuracy is not a high priority in this case, making a minimal triangle out of a single point is fine, even if not exact.


title: Looking for a web service for change detection

text: I'm looking for service that provides me a change detection service based on a database of remotely sensed data. Has anybody happen to know such a service? The service I am looking for is a geoportal that alarms the user of a specific change (like oil spill or major coastal change or even fire in forest) in a specific area or globally. It can be either freely available service or on- demand service.


title: How to determine intersection between geojson feature and given coordinates?

text: I'm overlaying a leaflet map with features of geojson, fetched with d3 and added as SVG elements to the DOM. My geojson contains too many features to render at once, so my goal is to add to the map only features, that intersect current cursor location, so I can have a smooth hover effect. I can get current map coordinates of the cursor from a mousemove event, but I have to determine if current coordinates intersect a feature in geojson and exactly what feature is intersected. Thinking further, maybe I have to make a small buffer, so it's easier for the user to hit the closest feature. Is there a fast way in leaflet or d3 (or even pure javascript) to get the intersected feature from geojson feature collection?


title: FEMA LOMA Data Availability?

text: Does anybody know where there might be some good/current FEMA data? I am interested specifically in the LOMA data for Florida. I have found a WMS Service that is hosted by FEMA but I'm having trouble actually bringing it into ArcMap.


title: Use Display Expression for KML Popup, SHP to KML ArcGIS

text: In Arcmap I've created a complex label expression in the Display expression for my layer. I'd like to have that expression displayed when I convert the layer to .kmz file for display in Google Earth. I'm using the ESRI Layer TO KML conversion tool. What happens in the popup right now is the display expression jammed, without the correct formatting, into the feature name box and also the feature description box, followed by the complete list of fields. I've been looking at the HTML popup tab, but don't see any way to use the Display Expression I already have in that tab. Is it possible for this to happen? It looks like otherwise I'll have to write a XSL template. EDIT: Thanks to @KHibma for the comment. What I've ended up doing is   1. removing all of the fields from display in the shapefile properties in ArcMap before the export so they are not included in the kml popup   2. Labeling the features in arcmap with one field so that field value is the "title" of the feature in the popup.   3. adding HTML formatting tags to the Display Expression (which I was writing in Python). this expression appeared, formatted, in the blue box below the feature title. I can't add a screenshot of the result right now, but will be able add it later. This worked fine on a small scale, but when trying to export an entire parcel fabric for a county I would get a C++ runtime error. One interesting thing I found was that, even before the export, if there was an error in the FindLabel function and it was not able to create a label, it would come back as empty. So I added a try/except statement within the function like:               def FindLabel([field1],[field2]):         try:             #code for function         except:             return "there was a problem creating this display"      which worked great. However, when I tried to import sys or traceback in the function and set the error handling message as the return in the except statement, it still returned blank (i.e. didn't return the text of the error handling message). In the end, I fixed all the problems with the expression, but still got the runtime error when making the big export.


title: Is there any accuracy difference between standard DTM editing tools like DT Master and automatic removal of artefacts using ERDAS?

text: I would like to know whether there is any accuracy difference while editing DTM points using standard photogrammetry software like Inpho DT Master and automatics removal of artefacts in ERDAS imagine.


title: When I used geowebcache, my map overviewmap can not use it?

text: I know: gwc / service / wms geowebcache use data, overviewmap is use wms original data, I tried modifying openlayers.js,some difficulties. You encountered the same problem and can help me look at it, thank you！


title: QGIS conversion of Microstation .DGN to .SHP

text: I'm attempting to convert a .dgn file to .shp. I'm reading that support is only provided for v7 microstation files, is this still the case? I've been using ogr2ogr -f "ESRI Shapefile" output.shp input.dgn and it throws the following error. Unable to open datasource `input.dgn' with the following drivers. -> ESRI Shapefile -> MapInfo File -> UK .NTF ... The list goes on. ANy help is appreciated.


title: What are the requirements for using Google's GEvent?

text: For a current project, we need to know when a Google map is done loading all of its tiles. I found the following question, which was for the same basic problem: register event "loadend" on layer So now, we're using the code from that answer:               GEvent.addListener(mGoogleHybridLayer.mapObject, "tilesloaded", function() {      console.log("Tiles loaded");     });      However, this gives a reference error for `GEvent`. Are there some extra requirements for using `GEvent`? Do we need to import it from somewhere? I can't post code we're using, because there isn't any yet - even in the browser's debug console, `GEvent` is undefined. Using Openlayers 2.11. EDIT: maybe I should restate the question into one that can actually be answered: what are all the requirements for being able to use GEvent?


title: Can I overlay layers in Dealul Piscului 1970 and WGS84 in QGIS?

text: I have a 10 Gb ECW file (mosaic of a whole county) that is in Oblique Stereographic projection with Datum: Dealul Piscului 1970. The problem was that no matter the Default CRS for new project was, when I wanted to add this raster the coordinates were all wrong, although the input projection for the raster was correct. This problem was solved by deleting the ECW's ".xml" file that was accompanying the ECW file. Once deleted, the problem was solved and the QGIS was showing the right coordinates. Once the above problem was solved. I wanted to add a vector layer of 30 mb that represents the whole waterway of the Danube Delta plus inland towns. This layer is in WGS84. Now when I try to import this layer does not overlay with the ECW file. The ecw file have been defined with the right projection also the vector layer have been defined with the right projection. In settings/options/CRS/ the "on the fly" is checked and the project projection is WGS84 and also tried the vice-verse with Oblique Sterographic Dealul Piscului. With all the permutations available I still can not overlay the layers. Do you have any other suggestions? Thank you


title: Convert a number in attribute table to binary

text: I would like to convert whole integers from decimal to binary in ArcGIS using the field calculator. I am looking for the equivalent of DEC2BIN excel command. All values are 255 or below so 8 bit is fine. It must have 0 padding to the left. So for example 4 must be 00000100. I can do it by exporting the attribute table to Excel and using the Dec2Bin command in Excel and rejoin the data but I would like it in a model so want to use the field calculator. Dec2Bin (Value, 8) works in Excel. Attached is how I would like it to look when complete. ![enter image description here](http://i.stack.imgur.com/DiMme.png) So basically change the value to and 8 digit binary value stored as a text field with leading 0s so all field have 8 values. I found this code in arcscripts but cannot get it to function. Dim valI as integer Dim valB as string Dim x as integer valI = [fieldName] for x = 7 to 0 step -1 if valI >= 2^x then valB = valB & 1 valI = valI - 2^x else valB = valB & 0 end if Next valB = valB & "00000000" fieldName = valB


title: Can I choose specific SDE version to create replica geodatabase from?

text: I'm trying to automate the process of exporting a replica geodatabase using the geoprocessing tools and eventually a Python script. I can't figure out how to create the replica from a specific version. Is this possible?


title: mapCanvas.saveAsImage not functioning as expected in QGIS 1.8

text: I have a plugin for exporting maps as PNG images. I believe this code was working correctly within 1.6 (but I cannot guarantee I haven't made another change that would have broken it there, too). Basically, I load a set of base map features and then loop through a set of point layers that get applied and the whole thing is exported as an image. I'm not getting the full canvas at all, but only a single attribute loaded from the base map and nothing from the point layers. (bad result, expected result) In the process of debugging this, I discovered a curious fact. If I include a QMessageBox.information right before I do the export, it works correctly! Doesn't work:                   self.mapCanvas.resize(QSize(394,350))         self.mapCanvas.zoomToFullExtent()         currentImagePath = theExportPath + folder         currentImageName = attribute.replace(' ', '_') + ".png"         self.mapCanvas.saveAsImage(currentImagePath + "/" + currentImageName)      Does work:                   self.mapCanvas.resize(QSize(394,350))         self.mapCanvas.zoomToFullExtent()         QMessageBox.information(self.iface.mainWindow(),"Debug","Break")         currentImagePath = theExportPath + folder         currentImageName = attribute.replace(' ', '_') + ".png"         self.mapCanvas.saveAsImage(currentImagePath + "/" + currentImageName)       I tried to add a simple self.iface.mainWindow() in place of the message box but that didn't do the trick. I'm thinking there is some kind of refresh I could do, too, but I tried mapCanvas.refresh() and mapCanvas.updateMap() and these in combination with mainWindow() but could get nothing to work. Did this change in 1.8? Any suggestions to make it work? Including the QMessageBox call will get me going but I have to press "OK" for each map to be exported. With all the export combinations, I'd be doing this more than 1000 times! Thanks! B


title: How to create a multipart polygon?

text: I have two polygons that consist of two separate polygons/lines in QGIS). How can I make this into a multipart polygon, which consists of 1 line in QGIS?


title: Cutting features in ArcGIS 10

text: What´s the easiest/smartest way of getting the opposite output of the clipping tool? I just need the area of feature _A1_ which **is not** covered by feature _A2_ in a new Layer/File. Now, I get the result by union´ing both features and deleting the "wrong" parts by hand.


title: how to make transparent layers with Google earth layer?

text: with your help, i solved most of my problems, but here is another one - i guess, a smaller issue... I managed to get Google earth layer, but when I add my GPS gpx files, are all hidden under the Google map... How to put them onto Google Earth layer? thank you in advance (and sorry for being a dumbass) Krepki


title: Why does intersect not calculate polygon areas?

text: this is strange: I intersected 2 shape-files and the resulting shape does _not_ contain the Shape_Area and Shape_Length of the intersected polygons. Instead the field Shape_Area contains the polygon-area of one of the 2 original Layers. Recalc of the geometries with vector->geometry->add/export column seems to add the areas. That is quite a drag and very slow for my 350.000 polygons. So I am hoping for a systematic approach that adds the areas while generating the intersect. Is intersect done by QGIS (qgis.bin) or by a gdal-tool? Either way: is there a "calculate areas"-option hidden somewhere? THX, stn


title: Extract polygon from several FC in one process

text: I've built a polygon and used Clip (Analysis Tools) in order to extract information from several FC that intersect the polygon. Of course, I used the clip tool for each feature class. My boss wants a single tool that clips all these feature clases at once and saves the result in different FC. Is there a tool in ArcGIS that does that? Select by Graphics is cool, but then I have to export the data "manually". I've created this model, but in the case of the model I don't have the liberty of selecting the area I want to clip everytime. Thank you


title: Good place to get Country, State, City, and Zipcode information?

text: I'd like to populate a database with common location information. I'm not entirely certain that addresses across the globe have similar or even the same entities as the Country, State, City, and Zipcode that we are so familiar with here in the United States. Does anyone know of a web service or other method that would be useful for obtaining such data and populating a local database with its contents? (or conversely just contacting the web service to obtain this information as needed)


title: Housenumber(point) angle from line (street)

text: using arcgis 10 with Point data (Housenumber without angle attribute) and line data (street without angle attribute) I want to create housenumber labels from point with correct angle from nearest streetline. is there any toolbox or script available for this, and if there is a script can you guide me step by step beacuse I dont know where to put script, and python would be the best I believe. my thinking of this is calculate AVERAGE streetline angle and then somehow put this values to nearest point? any help is good help :)


title: Problem in filling table of ArcSDE feature class in SQL Server 2008?

text: I have an ArcSDE point feature class and I add it to ArcMap 10. After that I add some points to it and fill every field for all points in the attribute table and Save Editing. However, the SQL Server table that related to the feature class is empty while the attribute table in ArcMap and ArcCatalog is full (as expected). I checked the connection to the spatial database and that was correct. How do I fill the table in SQL Server when the attribute table in ArcGIS is full?


title: How can I visualize a route in QGIS?

text: I have a route in QGIS that visits every road/edge at least once. By using a marker line I can see the direction of the route, but it doesn't show the order of the segments (see the screenshot below). Is there any style/visualization that will allow me to see the order of the roads in the route? Currently the route is a single polyline, but I can split it into segments if that will make it easier (the route is generated with Python code). ![Route example](http://i.stack.imgur.com/reec7.png)


title: How can I programatically create and add features to a memory layer in QGIS 1.9?

text: I had a working plugin in QGIS 1.8 which read data from an MSAccess database and added it to a series of memory layers. There is some processing involved in the interim so I don't think simply using QGIS to read directly from the database is an option. I would like to move from QGIS 1.8 to 1.9 (primarily due to the improved quality of the print composer). The plugin doesn't work with the new API. I've tried a variety of methods that came up in google searches. One, modifying the below code - from http://www.qgis.org/pyqgis- cookbook/vector.html#memory-provider, i.e. adding geometry and attributes to the dataprovider and then updating the layer - to suit the new API worked a little but the attributes were not visible until I manully entered editing mode (similar to http://hub.qgis.org/issues/3713). An alternative approach, detailed in reply #1 of the above link, added the layer and attributes correctly but I was unable to add features to the layer. Given that it should be a pretty simple task I'm hopeful that someone on here can offer a working example of how this should be done. (P.S. I am not a professional programmer and most of my coding is fairly crude - I welcome any guidance but ask that you excuse some ignorance on my part)               # Receivers = a list of lists returned from a database query          # create layer     vl = QgsVectorLayer("Point", item, "memory")     pr = vl.dataProvider()          # add fields     pr.addAttributes( [ QgsField("Rec_No", QVariant.Int), QgsField("Include",  QVariant.String), QgsField("Label",  QVariant.String), QgsField("X", QVariant.Double),                         QgsField("Y", QVariant.Double), QgsField("Z", QVariant.Double), QgsField("Height", QVariant.Double),                         QgsField("Project_Re", QVariant.String), QgsField("NCA", QVariant.String),                         QgsField("DayCrit", QVariant.Int), QgsField("EveCrit", QVariant.Int), QgsField("NightCrit", QVariant.Int) ] )          for i in range(len(Receivers)):                 # add a feature       fet = QgsFeature()       X = Receivers[i][3]       Y = Receivers[i][4]       fet.setGeometry( QgsGeometry.fromPoint(QgsPoint(X,Y)) )            # Details = a list of results returned from a database query specific to each result in 'Receivers'            if Receivers[i][3] != 0:         Include = 'Yes'       else:         Include = 'No'            fet.setAttributeMap( { 0 : QVariant(Receivers[i][0]), 1 : QVariant(Include), 2 : QVariant(Receivers[i][2]),                              3 : QVariant(Receivers[i][3]), 4 : QVariant(Receivers[i][4]), 5 : QVariant(Receivers[i][5]), 6 : QVariant(Receivers[i][6]),                              7 : QVariant(Details[0]), 8 : QVariant(Details[1]), 9 : QVariant(Details[2]), 10 : QVariant(Details[3]), 11 : QVariant(Details[4]) } )       pr.addFeatures( [ fet ] )          # add a style     vl.loadNamedStyle('C:/OSGeo4W/apps/qgis/python/plugins/Gopher2QGIS/styles/Receiver_Style.qml')          # update layer's extent when new features have been added     # because change of extent in provider is not propagated to the layer     vl.commitChanges()     vl.updateExtents()     vl.updateFieldMap()          QgsMapLayerRegistry.instance().addMapLayer(vl)


title: Web mapping application with time slider and raster data

text: I need to create a web mapping application that can show raster layer of historical maps with a time slider. Basically, something like this: http://goo.gl/GYHOX I have different maps of different years of the same are, that can be showed on a time series and moving on with a time slider. The interface can be simple, the important thing is the time slider. I have programming experience only in HTML/CSS, little of PHP. Any suggestions of which tools to use? I had a look at TileMill but the problem is the time-slider that I have no clue on how to build it. Thanks in advance for who will help me.


title: How specify a corridor minwidth?

text: I’m trying to make a corridor between two substations with ArcGIS, a macro corridor and a micro corridor. I have two questions in this respect and would be grateful for your suggestions. a) How can I specify, that my corridor must have a minwidth of 100m? b) How can I make a linestraighting of my costpath within microcorridor? Thank you very much.


title: How to output visually correlated clusters using Morans I Spatial Autocorrelation?

text: ![this is a screenshot of the geoprocessing output message](http://i.stack.imgur.com/Lqt7v.png) ![this is a screenshot of the attribute table result](http://i.stack.imgur.com/aDTkC.png) How do I use **Morans** I to output visually correlated clusters? I can seemingly generate statistical data, but not one that I know how to visually output on a map. I'm trying to run a morans I analysis on a collection of municipal level boundaries. As I understand the Morans I process, the 'input' feature class is whatever layer is going to be auto correlated, and the input field is whatever variable is going to be measured for the auto-correlation. In the case of the attached photos, I tried to auto-correlate on the basis of the variable R_vicMargNum. Approximately 35 out of the 80 rows in my attribute table do not have any values for that variable; I'm still hunting for the data which can properly fill the fields. 1) I don't understand why the fields which had a value for that variable are 'null' for columsn f7-f9. 2) I don't understand the value of the statistical output.


title: Change number of decimal digits for features in ArcGIS 10?

text: I am editing vertices in polygons and points in ArcMap 10, and realised that they only hold 3 decimal digits for the X, Y coordinates, but I would like at least 6 decimal digits. Even if I type in more digits, it gets rounded up to 3. How do I change this? Thanks.


title: How to create a buffer from a click on map on ARCMap

text: I'm searching a way to creating a buffer of a specified radius starting froma a click on the map.


title: Densify a JTS LINESTRING (add points)

text: I have a LineString I would like to "soften" in a way that if two coordinates (vertices) in the LineString are very far apart, I want to **automatically** add some coordinates in between. The line string should stay the same but with more coordinates describing it. Do you know of a way? ![enter image description here](http://i.stack.imgur.com/1FLCj.png) I have a LineString like the one on the left - and I want to add some coordinates in the between the existing coordinates (like on the right) - hopefully by JTS and hopefully automatically Thanks!!


title: Flex Identify tool version 3 1085 error

text: I am building a flex web map app. I am using the Identify tool version 3, here. Below is the syntax I am using, based on the PDF help file. It, however, gives a 1085 error. Does anyone know if I am using the correct syntax?               <layers onlythese="true">         <layer>             <name>Australia offices</name>             <fields>                 <field name="Address" />                                 <field name="City" />                                 <field name="State" />                                 <field name="Postcode" />                         </fields>             <links>                 <link includeinresults="false" field="" alias="" tooltip=""/>                      <linkprefix></linkprefix>                     <linksuffix></linksuffix>                     <iconfield></iconfield>                     <iconprefix></iconprefix>                     <iconsuffix></iconsuffix>                 </link>             </links>             <zoomscale>15000</zoomscale>             <forcescale>true</forcescale>                  </layer>                  <layer>                         <name>QLD rail</name>             <fields>                 <field name="FNAME" alias="Name"/>                                 <field name="ELECTRIFIE" alias="Electrified"/>                         </fields>             <links>                 <link includeinresults="false" field="" alias="" tooltip=""/>                      <linkprefix></linkprefix>                     <linksuffix></linksuffix>                     <iconfield></iconfield>                     <iconprefix></iconprefix>                     <iconsuffix></iconsuffix>                 </link>             </links>             <zoomscale>15000</zoomscale>             <forcescale>true</forcescale>         </layer>     </layers>


title: Openlayers map changelayer event firing twice for params

text: I have registered **changelayer** event of my map.                   var registerEvents = function () {             map.events.register("changelayer", null, onMapLayerChanged, true);         };              var onMapLayerChanged = function (event) {             console.log(event);         };      when I changed a layer's params as following, console writes two times.               layer.mergeNewParams({'CQL_FILTER': "name=wate"});


title: Problem with OpenLayers Tile Stitching Printing - All Black Image

text: I am trying to get an image of my OpenLayer application. I am using this http://www.mapsportal.org/olprint/ and http://trac.osgeo.org/openlayers/wiki/TileStitchingPrinting as a guide. When I click my print button I get a print failed message and when I go to the TMP folder (c:\wamp\www\tmp) there are IMG files and a JPEG but the JPEG is all black.


title: leaflet and geojson tiles

text: New to mapping. Trying to create a heat map with data that updates every hour for US Counties with leaflet.js similar to http://leaflet.cloudmade.com/examples/choropleth- example.html. But, adding all the US Counties in 1 shot makes map jerky in IE. So found this example of using GeoJSON tiles with leaflet.js, http://bl.ocks.org/3725681 Only problem is that many tiles are not rendered completely. I think some of the z/x/y.json files are missing coordinates, but not sure. Can someone point me to the right direction of the problem here.


title: Count the number of symbols in a 'new symbology' renderer from pyqgis

text: Background : I create a lot of symbology using 'saturation' ramps - for want of a better term, e.g. where the hue and value of each symbol stays the same, but the saturation changes. I couldn't figure out how to do this easily in qgis (EDIT: see comment - its easy when you know how), and each time I classify items, the colours reset to a colour ramp rather than a saturation ramp. I figured this should be relatively easy to do using pyqgis, which it turns out has been, using this code (built for 5 symbols):               from PyQt4.QtGui import QColor              # import this namespace as I want to refer to QColor later     layer = qgis.utils.iface.activeLayer()     i = 0                                       # used to increment the saturation value     for sym in layer.rendererV2.symbols():      # loop through each symbol in the layer using 'new symbology'         i += 50                                          sym.setColor(QColor.fromHsv(40,i,255))      Problem: . The above code works fine, but is not dynamic, in that I've set it up to increase the saturation by 50 each time (giving 5 graduations between 0 and 255). I want to set the saturation to graduate dynamically based on the total number of symbols. The mathematics of this are pretty straight-forward, but I'm not able to get the count from the symbols() method. I've tried:               >>> rendererV2.symbols().count()     Traceback (most recent call last):       File "<input>", line 1, in <module>     TypeError: count() takes exactly one argument (0 given)      ... which as you can see returns an error (although the QT documentation seems to list count() as a method of a QList<> object). I've tried:               >>> rendererV2.symbols().count(QgsSymbol)     0      ... which as you can see returns 0, which isn't true. There are 5 items in the list. I'm sure I must be doing something daft, but being new to python, pyqgis and qt, it is beyond my ability to see what at present. Help?


title: How to save QGIS composer template with image - e.g. company logo?

text: I can't figure out how to save QGIS Print composer template with image. In our map compositions we always have company logo in one place. It is easy to add image (logo) to composer and place it to spot I need. But when I save template, it saves without image (logo). Is there some setting for it? Am I misssing something or there is no way so far to save template including image?


title: How to Project shapefile from NAD 1983 UTM Zone 12 to WGS 1984 UTM Zone 13N?

text: I have a question and I have a shapefile that is in NAD UTM Zone 12 and Datum is North American Datum of 1983. What I would liket to reprojection it to NAD 1983 UTM Zone 13N and change the Datum to WGS 84 ? I have tried different methods but cane out the same. it seems like the Datum doesn't change, is it ? Should it be just one step to use the project to reproject to NAD83 UTM Zone 13 ? or is there more than just one step to get it there ? Currently I have 9.3.1 .


title: How can I conduct Geographically Weighted Analysis using an exponential decay function?

text: I have some points (representing 30 study sites) and I want to calculate the weighted mean of several variables (landscape layers) using a negative exponential decay (weights function) to give more weight to the landscape features closer to the points (study sites) than the features (grid cells) further away. (The series of underlying grid layers represent landscape features such as roads; habitat etc.) Also, how can I use spatial statistics to find the optimal distance (bandwidth) for my distance-decay function? In the literature everyone resorts to creating buffers (vector) at different distances (say 1km, 3km & 5km) to quantify the variables and then uses statistics to determine which buffer distance is significant. Other ecologists (Rhodes et al. 2006) have used the “negative exponential distance weighted density” of each variable and provided a scale parameter (L) which controls how rapidly the influence (ie weighting) of the variable declines with distance, eg L = EXP(-0.002*Distance) I’m using ArcGIS spatial statistics; grid (raster) spatial analysis tools; and investigating GWmodel for R-stats from an earlier post.


title: mapping and update data in database

text: We have two feature layers:   * layer A - drainage segment layer (along road outline)    * layer B - road centre line layer. In my SQL Server database I have a table storing the mapping information of segment ID (from layer A) and Road name (from layer B). ![enter image description here](http://i.stack.imgur.com/3sn5F.jpg) The problem is that when I create a new polyline feature (a segment) in layer A, I need to update the database table also. I am stuck here. I do not know how to update the database table at the same time as when a user creates a new feature in layer A. Should I create a new geoprocess (using python script) to do so? Is there any other solution for it? I am using vb.net arcmap addin (vs 2012), arcmap 10.2 and sql server 2008 r2.


title: Error in Conditional Tool of ArcGIS for Desktop?

text: I am trying to extract stream raster from a flow accumulation raster using Conditional Tool in Spatial Analyst. Here's my input:   1. Input Raster: FLOWACC.tif   2. Expression: Value >= 3.5   3. True Raster: 1   4. Output raster: Stream.tif Error gives me invalid expression. Basically my objective is to isolate all values >=3.5 and consider them as stream raster. I am using ArcGIS 10.2 for Desktop.


title: Qgis doufur, query problem

text: i've a question about a query I need to query results that are multiple of 50 like 100,150,200,300 How can i perform the right query?


title: How to create single contour shapefile from a contour shapefile?

text: Is there a way to create single contour shapefiles from a main contour shapefile. The main shapefile has quite a few contours and I need to isolate them and create individual files by elevation. Selecting contour by contour and creating a shapefile would take a while. Anybody knows of a faster way? Thanks!


title: Interacting between Openlayers and Python

text: I have a openlayers website like this one that gives the user the ability to draw polygons and then creates a GeoJSON file. I want to process this GeoJSON with a python script. The python script returns in the end a string that should be returned to the website and displayed as a info box. What kind of frameworks/apps/steps should I look into to accomplish this project?


title: Help getting my flex service publicly accessible

text: My organization is trying to get our web maps up and running. I have created a simple web map service on our in house ArcGIS server, using the Arcgis flexviewer that is visible locally, works great, and is visible on arcgis.com. We have a 2 server set up where Server A is our internal server that everyone in the company stores their data on. this is the machine we have our install of ArcGIS 10.1 for server on. It is not public, and for internal staff use only.Server B on the other hand is off site and set up by a company that hosts our website to specifically host our web maps. It is in a DMZ and set up to get through our company firewall. It is just a simple web server running IIS 7 and has the ArcGIS Web Adaptor installed .They, like me, have no experience in setting something like this up so I'm learning as much as I can along the way.. We have our flex map up and running [here] which is great. What is not so great is the error that is popping up when the flex map is opened. It looks like my data, which consists of a a few poly lines, some points and some labels, are not making it over to Server B from Server A. I guess my real question is what do I need to do to get the data from Server A show up in my flex app on server B. If anyone can shed some light on this that would be great.


title: Openlayer plugin Google Hybrid date

text: Comment connaître la date de capture d'une image Google Hybrid sur QGIS? How to know when capturing an image on Google Hybrid QGIS? **French Speakers - Please help fix the above machine translation**


title: NDS file format conversion to an executable GIS format

text: Could we open .nds files in any of the GIS mapping environment. it stands for Navigation Data Standards, having the navigation features stored in it. There are 2 types of nds file formats. One is executable one (in sqlite interface) and the other non executable (because of its binary format). I want to customize these nds files to view the data stored in it. I require a software technology to open these files and view the stored data on any of the mapping software which is compatible.


title: How to build a query tool for spatialite data in QGIS?

text: I am trying to build some tool for a collegue of mine to query data from a spatialite database while clicking on vector features. My first guess is, that I could do this with an "action" I define for the layer. The example actions provided for layers unfortunately do not give me any hint how to perform, and more detailed documentation on actions is quite absent. I need a query, cause each feature in layer has many corresponding entries in the database table (one-to-many). While performing the action, I would like to pop up simple table which lists the corresponding data sharing the same id as the feature. In the example action, i can query an attribute of the layer with this:               echo "[% "id" %]"      which then looks like this: ![enter image description here](http://i.stack.imgur.com/MxJ3W.jpg) Ok, that's a start. But what i need is a table which lists all items in a spatialite table with the same id-attribute as the feature. My data structure:               Layer attribute:      id          Data table:      id | name | year      The layer itself is stored in the spatialite database too, if thats of any help. Can anyone give me a hint on how to proceed or direct me to some other web ressources on this topic? Thanx Bernd


title: Automaticallly select and replace all occurences of a font face

text: I have around 50 mxd file that is done some times earlier. The text font used in labeling was not licensed. Now I can Change the text in layout element. But it has also some text in map view too. But arcpy.mapping.ListLayoutElements can not find the text document. Is there any solution in arcpy?


title: OpenLayers: how to display the character degrees in the label

text: I would like to display temperature in label. So how can I write DEGREE SIGN? to be display as ° I was trying write: & # 176; or & deg; or \ & deg; and nothing. Can anybody knows? I greet.


title: Field Type Date invalid when run in Python, but not in Model

text: Using ArcGIS Desktop Advanced 10.1 and PyScripter When running the arcpy.FeatureClassToFeatureClass_conversion, I get the error that the field type of Date is not valid. However, when running this same procedure in the model where the python was exported from, it runs successfully. I have verified that the field type is in fact Date for the field and that the field length is sufficient.   This seeems to be an issue after moving to 10.1 Error Message:               ExecuteError: ERROR 000278: Field Mapping error: Field type invalid for output field xxx and input field xxx      Any ideas?


title: FME: Can I add multiple feature types to the SQL Server Spatial reader?

text: The title's it, really. With .TAB file reader, you can add multiple layers to the reader and use the 'merge feature type' check box with an '*' to read all the layers listed in the source files, which can then be separated out again on featuretype name. With the SQL Server Spatial reader, you can specify many layers in the database selection, but it seems that the 'merge feature type' option doesn't work. Is this a bug, or expected behavior, or am I missing something? EDIT: I have found that I can do this sometimes... so the above may be incorrect. I'll need to test some more.


title: How to Create a Vector Field from Separate u and v Data Files?

text: I am using Quantum GIS and trying to create a vector field using data from two separate netcdfs. One contains the u data, one contains the v data. Can this be done, and if so, how? Thanks. Jason


title: Get the junction id for the roads in a network dataset

text: I have a network dataset with roads and junctions. I am wondering what class I can use in arcpy to extract FromJunctionId and the ToJunctionId for each road segment in the network dataset. Thanks.


title: Lost geometry when joining spatial and non-spatial tables

text: I have two tables I am joining: 1 spatial (Table 1) and one not spatial (Table 2). Table 1 looks like this: | ID | geometry | Table 2: | ID | ... | When I perform the query: > CREATE TABLE AS SELECT * FROM t1, t2 WHERE t1.ID = t2.ID; I get a table with the geometry column included. But unfortunately it is no longer recognized as a spatial table (spatialite_gui at least). Also, this mean I can't do what I really want to do, which is dump the table as shapefile. So my question is: Why is the new table not recognised as a spatial table with a geometry column?


title: How to write labels with varying significant figures/decimal places?

text: I have results that range from ex. 113,009 down to 0.00495 - Can someone recommend a way to write/construct a label expression for the following in ArcMap: • If less than 1, use 2 decimal places (or three if absolutely necessary to show the first non-zero value) • From 1 to 10, use 1 decimal place • If 10 or greater, use no decimal places Or a way to quickly generate a new text column/field which adheres to the above?


title: Why does DataGridView data not update for ArcMap add-in?

text: I have created a dockwindow（name:FrmWellConfig) with a datagridview(datagridview1) control in a new ArcMap add-in project. I then bind the data:               datagridview1.datasouce=wellconfig.tables[wellcolunms];      Why does the datagridview not update? Code behind tool：               FrmWellConfig frm = new FrmWellConfig(SettingWorkspace.AppBasePath, wellConfigFileName, frmSave.MapType, SettingWorkspace.workspace1.CnnStr_InfoDB);      Code behind FrmWellconfig(wellConfig1 is a dataset)：                /// <summary>     /// Designer class of the dockable window add-in. It contains user interfaces that     /// make up the dockable window.     /// </summary>     public partial class FrmWellConfig : UserControl     {         private string cnnStr_OLEDB;         private bool nFlag = false;           public static string appBasePath;         private string wellConfigFileName = "";         sylMapType mapType;                   public sylMapType MapType         {             get { return mapType; }         }         public WellConfig WellConfigData         {             get { return wellConfig1; }         }              public FrmWellConfig(object hook)         {             InitializeComponent();             this.Hook = hook;             this.dataGridView1.DataError += delegate(object sender, DataGridViewDataErrorEventArgs e) { };         }                   public FrmWellConfig(string AppBasePath, string WellConfigFileName, sylMapType _mapType, string CnnStr_OLEDB)         {                  InitializeComponent();             nFlag = true;                    frmType = sylChildFormType.DialogForm;             cnnStr_OLEDB = CnnStr_OLEDB;             appBasePath = AppBasePath;                mapType = _mapType;             dataGridView1.Columns[0].Visible = false;             wellConfigFileName = WellConfigFileName;             wellConfig1.ReadXml(wellConfigFileName);         }              /// <summary>         /// Host object of the dockable window         /// </summary>         private object Hook         {             get;             set;         }              /// <summary>         /// Implementation class of the dockable window add-in. It is responsible for          /// creating and disposing the user interface class of the dockable window.         /// </summary>         public class AddinImpl : ESRI.ArcGIS.Desktop.AddIns.DockableWindow         {             private FrmWellConfig m_windowUI;                  public AddinImpl()             {             }                  protected override IntPtr OnCreateChild()             {                 m_windowUI = new FrmWellConfig(this.Hook);                 return m_windowUI.Handle;             }                  protected override void Dispose(bool disposing)             {                 if (m_windowUI != null)                     m_windowUI.Dispose(disposing);                      base.Dispose(disposing);             }              }


title: Is it possible to export boundaries from Google Map marker?

text: I primarily have a doubt that is it possible to export a marked boundaries from google map marker ?if possible kindly tell how to to do this?


title: Add Surface Information in 3d-analyst

text: One of our guys is using the "Functional Surface->Add Surface Information" feature in 3D-analyst tool in ArcMap 10.1 to add elevation data to an XY list. It's the only feature we're using from the 3D-analyst. Our trial period has ended and we're trying to figure out if we need to buy the whole 3D-analyst tool for $2500 just for this feature, or if there's a cheaper way to do this task since we don't use anything else in 3D-analyst. Since it's a pretty simple task, I might just develop a tool myself to do this using the USGS web service given the cost of 3D-analyst. Is there a cheaper / easier way to do this?


title: How to create and re-use my own map template in Arcmap 10.1?

text: I am working on a project which requires 5 different analyses on the same data. I created five map documents (.mxd) and completed my analyses. Now, the problem is that I want to export my results in map format but I want them to be uniform in display i.e. having a same map template. Is there a way to save one of the maps as a template and then use it for creating other 4 maps? I tried saving one document in 'Map Templates' folder in ArcGis folder but a message box appears which says that I don't have permission to save in that location and I must contact administrator for that, even though I am the only user and administrator of the computer. I don't want to manually add all map elements for 5 different maps, thats too tedious!


title: How to reopen a plugin-window in QGIS?

text: I am currently taking my first steps in playing around with programming a plugin in python 2.7 for QGIS 1.8. While my plugin does what it should, I would like to make it more userfriendly. At the moment, the user is expected to choose an input- and an output-file, which are then processed. In case the user forgets to choose one (or both) of the files, the plugin is displaying an information window - and then quits completely after the user hits "OK".               if src_filepath =="" and dst_filepath =="":         QMessageBox.information(None, "DEBUG:","Please specify files!")`      However, I would rather like to reload the main plugin-window with the "choose file"-dialog. I suspect that the solution is rather easy, but I was not able to invoke this window, so any help is appreciated. The interface code derives more or less directly from the "Plugin Builder" in QGIS.


title: Force dimensionality with shp2pgsql

text: After invoking: `shp2pgsql -W "latin1" theshape.shp > thesql.sql` My console responds with:               Shapefile type: PolygonZ     Postgis type: MULTIPOLYGON[2]      Note the shapefile type is `PolygonZ`. The result is, I have to manually edit the `thesql.sql` file and wrap the geometry column in `ST_Force_2D()` before executing the SQL, otherwise, Postgres complains I'm trying to import a 3D shape into a 2D column. How can I force dimensionality to 2 with `shp2pgsql`? Alternatively, can I force it in QGIS?


title: Upload layers in Geoexplorer

text: I have OGS 4.02, OS is Windows Server 2008 R2 Standard.I have questions about Geoexplorer. I know that anonymous users can not save and publish map and also they can not upload layers. If I create users which are in Geoserver admins (access mode is admin for specific workspace) or editors (access mode is write for specific workspace), why they can not upload their layers (zipped shapefiles) in their workspace, creating maybe even new stores on the fly. Why is that option only avaliable for one person, user with admin role? Is it something with folder security in Windows, is it necessary to create also users and specific folders for them. I know that after succesfull upload, two folders are created: one in C:\ProgramData\Boundless\OpenGeo\geoserver\uploads (tmp prefix) and another in workspace where data are uploaded. So who can only write in that folders through geoexplorer and geoserver port 8080? Is it necessary to create postgis store to enable upload for some power users (editors or small admins)?


title: TableToTable gives unexpected character after line continuation character?

text: I'm working on updating a python script from arcmap 9.3.1 to 10.2.1 and im getting an error for "unexpected character after line continuation character" and im unable to diagnose the character.               arcpy.TableToTable_conversion(sde_parcel_base_temp, sde, "parcel_base", "", "pin \"PIN\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",PIN,-1,-1;pwd \"PWD\" true true false 13 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",PWD,-1,-1;add_num \"Address Number\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",add_num,-1,-1;add_prefix \"Address Prefix\" true true false 1 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",add_prefix,-1,-1;add_strt \"Address Street\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",add_strt,-1,-1;add_type \"Address Type\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",add_type,-1,-1;add_sdir \"Address Suffix\" true true false 2 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",add_sdir,-1,-1;add_munic \"Address Municipality\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",prop_city,-1,-1;add_zip \"Address ZIP\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",prop_zip,-1,-1;add_full \"Full Street Address\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",prop_street,-1,-1;ownr_name1 \"Owner Name 1\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",owner1,0,49;ownr_name2 \"Owner Name 2\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",owner2,0,49;ownr_addr \"Owner Address\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",own_street,-1,-1;ownr_city \"Owner City\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",own_city,-1,-1;ownr_state \"Owner State\" true true false 2 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",own_state,-1,-1;ownr_zip \"Owner Zip\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",own_zip,-1,-1;acres "\Legal Acreage\" true true false 50 Double 0 0 ,First,#,"+ sde_parcel_base_temp +",legal_acreage,-1,-1;prop_code \"Property Class Code\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",property_class,-1,-1;nghbr_code \"Neighborhood Code\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",neighborhood,-1,-1;lgl1 \"Legal Line 1\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lgl2 \"Legal Line 2\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lgl3 \"Legal Line 3\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lgl4 \"Legal Line 4\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lgl5 \"Legal Line 5\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lgl6 \"Legal Line 6\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lrsn \"Proval LRSN\" true true false 50 Long 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;updated \"Last Update\" true true false 50 Date 0 0 ,First,#,"+ sde_parcel_base_temp +",last_update,-1,-1", "")      Above is the part of the script giving me problems. Anybody with a trained eye able to point out my problem for me?


title: Create a subset of a dataset with a specific grid size in ArcGIS using Data Management Tool

text: I am wondering how you can create a subset of a dataset that is a specific grid size (in pixels) using the Clip tool within the Data Management Tool box for raster data. I was trying to find a hint where to set the grid size in the tool's Environment Settings, but the result, after running the tool, did not apply the new grid size. Please does anyone handled this before? Thanks.


title: SQL Command to convert pgRouting result to GeoJSON

text: My goal is to produce a Web UI for pgRouting using this tutorial. I have already imported my data from an OSM pBF to PostGIS using osm2po. Using the pgRoutingLayer plugin for QGIS, I was able to figure out the SQL code to create routes using driving time as the cost.               SELECT ways.*, route.cost AS route_cost FROM ways JOIN (         SELECT * FROM shortest_path(         'SELECT gid AS id,                  source::int4 AS source,                  target::int4 AS target,                  cost::float8 AS cost,                 reverse_cost::float8 AS reverse_cost FROM ways',                 62,19169, True, True))          AS route ON ways.gid= route.edge_id;      My next step is to convert this into GeoJSON as done in Part 1 of the tutorial. This is what seemed logical:               SELECT * INTO export     FROM ways.*, route.cost AS route_cost FROM ways JOIN (         SELECT * FROM shortest_path(         'SELECT gid AS id,                  source::int4 AS source,                  target::int4 AS target,                  cost::float8 AS cost,                 reverse_cost::float8 AS reverse_cost FROM ways',                 62,19169, True, True))          AS route ON ways.gid= route.edge_id;     SELECT ST_AsGeoJSON(ST_UNION(the_geom)) AS geom_union     FROM export;          ERROR:  syntax error at or near ","     LINE 2:  FROM ways.*, route.cost AS route_cost FROM ways JOIN (                         ^      but I can't get it to work. Honestly, I'm not very well acquainted with SQL, can someone point out my mistake? edit: This seems to work, but does it do the same thing?               CREATE TEMPORARY TABLE route_result AS      SELECT ways.*, route.cost AS route_cost FROM ways JOIN (         SELECT * FROM shortest_path(         'SELECT gid AS id,                  source::int4 AS source,                  target::int4 AS target,                  cost::float8 AS cost,                 reverse_cost::float8 AS reverse_cost FROM ways',                 62,19169, True, True))          AS route ON ways.gid= route.edge_id;      SELECT ST_AsGeoJSON(ST_UNION(the_geom)) AS geom_union      FROM route_result;


title: How would I get number of tiles loaded, and total number loading?

text: So new to OpenLayers 3, and I'm trying to create a "loading bar" to indicate how much of a layer has loaded. I'm thinking of making it based on the number of tiles loaded so far divided by the total number of tiles that need to be loaded to render that whole layer from the WMS server-- but I've been looking around the API and I don't know how to get either of those numbers in OpenLayers 3. Does anyone know anything about this, or alternatively does anyone know a different approach I could use to create a layer loading bar? Thanks in advance! :]


title: Running a batch script in GRASS command line - v.dissolve and v.db problems

text: I have been writing a number of scripts for automating isochrone production in GRASS just using the command line. I create files with isochrone lines in them and then a buffer around each of the lines. I then want to dissolve the buffers, so use the command: v.dissolve input=AOS_buffer_XXX output=AOS_buffmerge_XXX Each time this completes it's first file, the script stops. I then want to put an attribute into each line that reflects the filename and I thought I could do this with the command: v.db.addcol map=AOS_buffmerge_XXX columns=GTW_ref INT [--verbose then v.db.update map=AOS_buffmerge_XXX column=GTW_ref value=XXX [--verbose] However when I run this on the several hundred shapes, it ends after the first execution - it seems to allow me to do one, but not significant numbers. Can anyone help me out with why these three functions in GRASS command line won't script? (I'm using GRASS 6.4.2 - would an upgrade help?) Example of the script (I have put annotation where the script just stops - no error message, it just stops):               v.in.ogr dsn=C:\0Data\FK_openspace\network\cleanedfalkirknetwork.shp output=fk_net009 --verbose -–overwrite               v.in.ogr dsn=C:\0Data\FK_openspace\OpenSpaceAndGates\Calculations\AmenityOpenSpace.shp output=AOS_gates1 where="Reference = 1" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\OpenSpaceAndGates\Calculations\AmenityOpenSpace.shp output=AOS_gates2 where="Reference = 2" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\OpenSpaceAndGates\Calculations\AmenityOpenSpace.shp output=AOS_gates3 where="Reference = 3" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\OpenSpaceAndGates\Calculations\AmenityOpenSpace.shp output=AOS_gates4 where="Reference = 4" --verbose --overwrite          v.net input=fk_net009 points=AOS_gates1 output=AOS_fixednet1 operation=connect alayer=1 nlayer=2 thresh=500 --verbose --overwrite     v.net input=fk_net009 points=AOS_gates2 output=AOS_fixednet2 operation=connect alayer=1 nlayer=2 thresh=500 --verbose --overwrite     v.net input=fk_net009 points=AOS_gates3 output=AOS_fixednet3 operation=connect alayer=1 nlayer=2 thresh=500 --verbose --overwrite     v.net input=fk_net009 points=AOS_gates4 output=AOS_fixednet4 operation=connect alayer=1 nlayer=2 thresh=500 --verbose --overwrite          v.net.iso input=AOS_fixednet1 output=AOS_output1 ccats=0-99999999 costs=400 --verbose --overwrite     v.net.iso input=AOS_fixednet2 output=AOS_output2 ccats=0-99999999 costs=400 --verbose --overwrite     v.net.iso input=AOS_fixednet3 output=AOS_output3 ccats=0-99999999 costs=400 --verbose --overwrite     v.net.iso input=AOS_fixednet4 output=AOS_output4 ccats=0-99999999 costs=400 --verbose --overwrite          v.out.ogr input=AOS_output1 type=line dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS layer=1 format=ESRI_Shapefile --verbose      v.out.ogr input=AOS_output2 type=line dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS layer=1 format=ESRI_Shapefile --verbose      v.out.ogr input=AOS_output3 type=line dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS layer=1 format=ESRI_Shapefile --verbose      v.out.ogr input=AOS_output4 type=line dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS layer=1 format=ESRI_Shapefile --verbose           v.in.ogr dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\AOS_output1.shp output=AOS_line_1 where="cat = 1" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\AOS_output2.shp output=AOS_line_2 where="cat = 1" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\AOS_output3.shp output=AOS_line_3 where="cat = 1" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\AOS_output4.shp output=AOS_line_4 where="cat = 1" --verbose --overwrite          v.buffer input=AOS_line_1 output=AOS_buffer_1 type=line distance=35 --overwrite --verbose     v.buffer input=AOS_line_2 output=AOS_buffer_2 type=line distance=35 --overwrite --verbose     v.buffer input=AOS_line_3 output=AOS_buffer_3 type=line distance=35 --overwrite --verbose     v.buffer input=AOS_line_4 output=AOS_buffer_4 type=line distance=35 --overwrite --verbose          v.dissolve input=AOS_buffer_1 output=AOS_buffmerge_1 --verbose     v.dissolve input=AOS_buffer_2 output=AOS_buffmerge_2 --verbose     v.dissolve input=AOS_buffer_3 output=AOS_buffmerge_3 --verbose     v.dissolve input=AOS_buffer_4 output=AOS_buffmerge_4 --verbose          v.db.addcol map=AOS_buffmerge_1  columns=GTW_ref INT --verbose      v.db.addcol map=AOS_buffmerge_2  columns=GTW_ref INT --verbose     v.db.addcol map=AOS_buffmerge_3  columns=GTW_ref INT --verbose     v.db.addcol map=AOS_buffmerge_4  columns=GTW_ref INT --verbose          v.db.update map=AOS_buffmerge_1 column=GTW_ref value=1 --verbose     v.db.update map=AOS_buffmerge_2 column=GTW_ref value=2 --verbose     v.db.update map=AOS_buffmerge_3 column=GTW_ref value=3 --verbose     v.db.update map=AOS_buffmerge_4 column=GTW_ref value=4 --verbose          v.out.ogr input=AOS_buffmerge_1 type=area dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\final layer=1 format=ESRI_Shapefile --verbose --overwrite     v.out.ogr input=AOS_buffmerge_2 type=area dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\final layer=1 format=ESRI_Shapefile --verbose --overwrite     v.out.ogr input=AOS_buffmerge_3 type=area dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\final layer=1 format=ESRI_Shapefile --verbose --overwrite     v.out.ogr input=AOS_buffmerge_4 type=area dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\final layer=1 format=ESRI_Shapefile --verbose --overwrite


title: OpenLayers.Format.GML parser problem

text: I'm trying to parse WFS reponse using `OpenLayers.Format.GML` parser. However, an exception thrown in my try-catch block is `TypeError: a is undefined`. Any pointers as to what I may be doing wrong will be appreciated.               protocol = new OpenLayers.Protocol.WFS({         version: "1.1.0",         url: WFS_HOST,         featureType: "pl",         featureNS: "ws",         srsName: "EPSG:900913"     });          wfsLayer = new OpenLayers.Layer.Vector("PE Vector", {         strategies: [new OpenLayers.Strategy.BBOX()],         protocol: protocol     });          var _Callback = function(resp) {         console.log(resp);         console.log(resp.error);              try {             var gmlParser = new OpenLayers.Format.GML();             gmlParser.extractAttributes = true;             var features = gmlParser.read(resp.responseText);                  if(features){                 for(var i = 0; i < features.length; i++) {                     var feature = features[i];                                         // More operations                 }             }         } catch(e) {             alert("Error: " + e);         }     };          response = protocol.read({         maxFeatures: 100,         callback: _Callback     });      My environment: GeoServer 2.1.3, OpenLayers 2.11, PostGIS 1.5.


title: How to create a hole in a polygon in QGIS?

text: I'm trying to use the tool "delete ring" in order to create a hole in a polygon. The geometry of the hole is another feature of the same shapefile I cannot manage to delete the part of the polygone. Does someone have a tutorial or some explanation?


title: ECW on QGIS 1.8 (Ubuntu 11.10)

text: I'm trying to deal with an ECW file with QuantumGis 1.8 Lisboa, on Ubuntu 11.10. Obviously it tells me that .ecw is not a supported kind of file on QuantumGis. I have tried to install all the libraries but it still doesn't work. What sould I do? Thanks, Jacopo


title: Using render_list to render only map extract

text: I am trying to render an area where I have added an extract from the planet file with render_list. The command I am using is               render_list -x 159 -X 179 -y 160 -Y 410 -n 2 --socket=/var/run/renderd/renderd.sock --min-zoom=10 --max-zoom=10      to render zoom level 10. I already rendered 0-9 zoom levels for the whole planet, but want to render only the areas where I imported maps for the bigger zoom levels. Running the command above gives message "Starting 2 rendering threads" and stays for a long time without doing anything. Any ideas why?


title: Postgis nested query trouble

text: I am having difficulty writing a nested postgis query. I have a single table from which I need to create a view in PostgreSQL. I have created the view with the following command               create view rrop08.rrop08_lws_f AS     select * from rrop08.rrop08_gk3_polygon      where planznr='04.01' OR planznr ='04.02';      The problem here is that the table contains another field called "bez". I need to carry out the first query but also with bez=0; so that none of the 04.01 or 04.02 rows have a bez value other than 0. I have tried this...               create view rrop08.rrop08_lws_f AS     select * from rrop08.rrop08_gk3_polygon      where planznr='04.01' OR planznr ='04.02' AND bez=0;      but this retrieves all rows with bez=0 and not just the 04.01 or 04.02 ones. I need to first select 04.01 or 04.02 and then on this query get those with bez=0...ie.nested. I´m a bit lost on how to do this. THanks for any help, Rob


title: How can i resolve problem with bounds for export?

text: How could I set up bounds in TileMill if i don't want to upload the country around? I simply want to upload my created island in `TileMill` to `Mapbox`. Is there a way to do it? I am using OSM Bright for the design of my map.


title: Calculate Field python code has "Expected end of statement" error

text: I am creating a line shapefile from a csv in a Python (2.7) script. I am assigning an ID field from the csv to a variable and then later using it in the **Calculate Field** tool. The ID field originally looks like this '011947-OCT-084' which ArcGIS(10.1) does not seem to like, so I removed the '-'s and when I use that variable in the tool I get an error: > ExecuteError: ERROR 999999: Error executing function.   >  Expected end of statement   >  Failed to execute (CalculateField). Here is my code ("temp" is the shapefile temporarily stored in in_memory):               with open(inCSV, 'rb') as input:         with open(outCSV, 'wb') as output:             csvwriter = csv.DictWriter(output, headers, delimiter = ',')             csvwriter.writeheader()             for row in csv.DictReader(input):                 uniqueID = row['Common_id']                 newID = uniqueID.replace("-", "")                 arcpy.AddField_management("temp", "uniqueID", "TEXT")                 arcpy.CalculateField_management("temp", "uniqueID", newID)                 row['Common_id'] = newID                 csvwriter.writerow(row)      I'm not entirely sure what the cause behind this is, as I've used variables in the Calculate Tool field before without a problem. Thoughts?


title: IL East 1201 coordinate system

text: What is the difference between NAD CORS96 (Illinois East 1201) and the "normal NAD1983? We have data that was GPS'ed in CORS96, but data was imported into ArcMap using the "normal" NAD1983? What effect, if any, will it have on data? For example, will it affect the location and the area calculations of the data?


title: Does arcpy have a method for accessing related records within a cursor?

text: I have a line feature class that represent sampling segments. I have a data table of field notes in the same Geodatabase that is connected to the feature class with a Relationship Class. Is there a built in method for arcpy to access related records? The work around would be to add another cursor for the second table, but you'd think that if you take the time to relate your data there would be a quick and easy method for accessing related records programatically. Software: ArcGIS 10.1