title: Join .csv file and .shp of georeferenced points

text: I have georeferenced a map and saved a selection of some points on the map. These points are cities of departure and arrival and I want to join the .shp with a and excel .csv with the information of the migration of my individuals to represent their way. -If I join on my .shp (the selected points and the origin cities on .csv) I only get 15 out of 40 lines as many individuals started their trip from the same city -If I joint on my .csv I have a complete table attributes. I save a selection and want to open it as .shp but I only have three extensions (.dbf, .prj, .qpj) How can I get a complete joint table attributes between a .cvs of 40 items and a .shp of 15 items and make sure the items are georeferenced on my map? I hope this makes sense, I'm still a beginner ;) Thank you for your help


title: Problem using shortest_path with find_nearest_link_within_distance function

text: I'm tried to work on the scenario that avoiding go through some path if there is a barrier near that path. It can run well with the query below:               SELECT * FROM shortest_path('SELECT gid AS id, source::int4, target::int4, shp_leng::double precision AS cost         FROM network1         where gid <> 35'         ,18, 36, false, false);      Which means that the the path with gid 35 in the network1 table won't be taken into the calculation of shortest_path function. But now, if I want to use find_nearest_link_within_distance function to fine the nearest path by a specific point, I got some problem, my code is like this:               SELECT * FROM shortest_path('SELECT gid AS id, source::int4, target::int4, shp_leng::double precision AS cost         FROM network1         where gid <>          (find_nearest_link_within_distance(     'POINT(120.218325 22.998469)',     0.5,     'network1'))'         ,18, 36, false, false);      The error information:               ERROR:  syntax error at or near "POINT"      LINE 5: 'POINT(120.218325 22.998469)', Seems like find_nearest_link_within_distance can't be put between ' '. I can't find the solution, can someone help me out with that?


title: Create polygon of buildings from Quickbird images in QGIS

text: I have Quickbird/Ikonos/GeoEye image and want to extract all the buildings from the image as an independent layer in QGIS.   How can i achieve that?


title: combine between WFS and WMS

text: we developing a project which contains a lot of data ,   using WMS layers for displaying the points,   Also we used WMS to return points information in a popup window ,   however to highlight selected point we used WFS layer. **The problem now is** :   sometimes the popup displayed without highlighting the selected point and vice versa ..! how could i combine between these tow features (highlight selected point and show the popup for it)? notes:   i think one reason is the clicking event defined different,   WFS using:   `selectControl.events.register("featureselected", this, function(e) { select.addFeatures([e.feature]);}` and WMS using:   `popupwindow1=new OpenLayers.Control.WMSGetFeatureInfo({   layers: [layer1],   autoActivate: true,   maxFeatures: 3,   infoFormat: "application/vnd.ogc.gml",   eventListeners: {....}` **The second problem is** :   how could i return features from WFS layer without adding the vector layer to the map?


title: Remove unused references in ArcGIS 9.3.1 vba

text: Can anyone suggest me how will I find and remove unused references used in a project done in ArcGIS 9.3.1 vb6. Thanks Khaled


title: GeoJSON to Shape file conversion with geotools

text: Can someone please provide me an example of how to convert GeoJSON to Shape file using geotools?


title: Comparing two Digital Elevation Models (DEMs) from LAS files?

text: I have two Lidar las files, one is original let's say with X points. And the other is copy of the first las file but with Y points, where Y is less than X. Now, I want to compare how the Digital Elevation Models of these two las files vary... I want to get information like RMSE, standard deviation, etc... I would appreciate, if anyone could tell me what softwares, or ways to get the comparison info...


title: Shapefile import with Oracle + GeoServer produces points instead of vector shapes

text: I have implemented the example at http://www.gistutor.com/openlayers/22-advanced-openlayers- tutorials/47-openlayers-wfs-t-using-a-geoserver-hosted-postgis-layer.html and got it working properly. Now I am trying to get the same example running, but with **Oracle Spatial** (or Locator). I imported the test data shapefile to my Oracle DB using Oracle Map Builder (https://forums.oracle.com/forums/thread.jspa?threadID=1030707&start=0&tstart=0). The data imported without errors, and I was able to configure GeoServer layer to use Oracle datasource for the data. However, when I look at the data with GeoServer's "layer preview", it's just points when it should be vector shapes. This is how my layer looks like with Postgis (and how it should be): ![enter image description here](http://i.stack.imgur.com/aWpL9.png) And this is how it looks like with Oracle: ![enter image description here](http://i.stack.imgur.com/i49Es.png) Any ideas on what went wrong? Loading the shapefile, configuring Oracle, configuring GeoServer, or what...?


title: Openlayers - Vector layer polygon doesn't move in-synch with Google V3 baselayers

text: When a polygon is drawn on the map and the map is dragged, the position of the polygon relative to the GM V3 base changes.               function init() {         var map = new OpenLayers.Map('map');         var gsat = new OpenLayers.Layer.Google(             "Google Satellite",             {type: google.maps.MapTypeId.SATELLITE, numZoomLevels: 22}         );         drawVctr = new OpenLayers.Layer.Vector("Editable Vector");         map.addLayers([gsat, drawVctr]);              // Google.v3 uses EPSG:900913 as projection, so we have to         // transform our coordinates             var lng = "-71.118004";          var lat =  "42.389118";          map.setCenter(new OpenLayers.LonLat(lng, lat).transform(             new OpenLayers.Projection("EPSG:4326"),             map.getProjectionObject()         ), 19);         drawPlygn = new OpenLayers.Control.DrawFeature(             drawVctr,             OpenLayers.Handler.Polygon         );               map.addControl(drawPlygn);     }


title: How to use raw_input() in QGIS python console?

text: I'm hoping to use raw_input() inside a loop to prompt the user when s/he wants to continue running the script. However, I'm getting an error:               Traceback (most recent call last):        File "<input>", line 1, in <module>     EOFError: EOF when reading a line      Is this function supported in qgis' python console? If not, do you have any suggestion? I don't need to read the user's input. I just need the loop to pause for a while to let the user do some other stuff, then when ready, to resume working again.


title: Import xyz data from csv into gdal

text: I have many CSV file in the form of lon, lat, Z that in would like to import into gdal for further processing. I need to do this on the fly in a python script so any step involving gdal_grid or anything else from the command line is not an option. It seemed that creating a vrt xml would solve the problem but I am unable to get gdal to load the data. Loading data this way seems to be such a common task I can't believe there is not a way to do this.


title: Exporting Leica Geo Office points to view/edit in ArcGIS for Desktop?

text: I've got some points recorded with a leica 1200 rover, which I now have imported into **Leica Geo Office**. I want to export these to view and edit them in ArcGIS for Desktop. I recorded their properties as annotations. From Geo Office I can export them as either CSV or shapefile (plus some options that I think aren't viable). However the annotations do not seem to come along when using the shapefile export. Also the location of the points seem to be closer to Brazil rather than somewhere in East England. No prj file is generated along with the shp. If I export the data to csv I do get the annotations, but get coordinates which are not familiar at all (as I thought it would just be WGS84 lat lon decimal degrees): X 3851050 Y 83323 I have hardly any experience at all with Geo Office, this is basically my first touch... As I have ~60 copying the annotations isn't that much of a strain but getting points on the right location would be nice. Anyone out there using Geo Office and able to help me out?


title: Extract coordinates for series of georeferenced files in ArcGIS 10?

text: I have a series of 1900 spatially continuous images that are georeferenced. Each image is a tif file paired with a tfw file. These images are easily open and shown in ArcMap 10. However, I need a way to extract and tabulate all of their associated coordinates (in this case, the centerpoints of each image) into a spreadsheet. How do I do this? Please note that I have no scripting abilities. Thanks!


title: ArcPy: DDP Error in exporting only occurs with "SELECTED"

text: I'm trying to use DDP (Data Driven Pages) and am running into a problem. I'm using arcpy.SelectLayerByAttribute_management to select a series of records, and then trying to use the ddp.exportToPDF(pdfddp, "SELECTED") but keep running into an error (below). It shouldn't be the select by attribute function as I run it in ArcMap and it runs correctly and I can see the selected features highlighted. Also, I can successfully use the ddp.exportToPDF(pdfddp, "CURRENT") but what I need is to be able to have a user select a series of features and have those be passed.                mxd = arcpy.mapping.MapDocument("C:\\temp\\Graffiti_DM4.mxd")     df = arcpy.mapping.ListDataFrames(mxd, "*")[0]     pdfddp="C:\\temp\\pdfddp1.pdf"     arcpy.SelectLayerByAttribute_management("C:\\temp\\test.lyr","NEW_SELECTION", '"Incident_Z" ='  + "'10002'")     ddp.exportToPDF(pdfddp, "SELECTED")      And here is the error               Runtime error      Traceback (most recent call last):       File "<string>", line 1, in <module>       File "c:\program files\arcgis\desktop10.1\arcpy\arcpy\utils.py", line 181, in fn_         return fn(*args, **kw)       File "c:\program files\arcgis\desktop10.1\arcpy\arcpy\_mapping.py", line 446, in exportToPDF         return convertArcObjectToPythonObject(self._arc_object.exportToPDF(*gp_fixargs((out_pdf, page_range_type, page_range_string, multiple_files, resolution, image_quality, colorspace, compress_vectors, image_compression, picture_symbol, convert_markers, embed_fonts, layers_attributes, georef_info, jpeg_compression_quality, show_selection_symbology), True)))     AttributeError: PageLayoutObject: Error in exporting pages      As always, any help is EXTREMELY appreciated!


title: OSM Nominatim Server Installation Projection Error

text: Environment: Ubuntu 12.4 Postgresql 9.1.4 PostGIS 2.0.1 Latest Nominatim installation instructions posted at http://wiki.openstreetmap.org/wiki/Nominatim/Installation The installation is failing with (command built by the Nominatim php build script): /osm2pgsql/osm2pgsql -lsc -O gazetteer --hstore -C 1100 -d nominatim /home/user/osm_data/central-america.osm The -l tag is causing the following error: Projection code failed to initialise osm2pgsql SVN version 0.80.0 (64bit id space) ERROR: Error executing external command: /home/user/Nominatim/osm2pgsql/osm2pgsql -lsc -O gazetteer --hstore -C 1059 -d nominatim /home/user/osm_data/central-america.osm.pbf In trying to troubleshoot, where is the -l tag pulling the projection code from? I can get the import to run without the -l but the import later fails with an SRID error which I can only suppose results from this initial error. Can anyone enlighten me on what I should be looking at to fix the issue? Thanks in advance.


title: How can I make a polyline wrap around the world?

text: I am using leaflet maps to create a representation of a round-the-world challenge. I would like to add a polyline which heads east from tokyo and then appears west of south america on the map - but instead I get a line which crosses the map in the opposite direction (see yellow line). ![enter image description here](http://i.stack.imgur.com/2QuR9.png) I think this is probably related to datelines and /or coordinate systems but am a bit sketchy on the detail. Can anyone explain the theory behind what I need to do to get this to work? I am using the bluemarble projection from Nasa:               var bluemarble = new L.TileLayer.WMS("http://demo.opengeo.org/geoserver/wms", {     layers: 'bluemarble',     attribution: "Data © NASA Blue Marble, image service by OpenGeo",     minZoom: 2,     maxZoom: 5     });


title: arcpy.mapping Graduated Colors Symbology

text: I'm trying to figure out how to apply graduated color symbology to a feature class by using arcpy.mapping in ArcMap 10.1. To start, a script runs and a feature class is added to a map document, no problem. This feature class needs the graduated colors. The value field used for the graduated colors is created during the script. In order to apply the graduated colors, I have to create a layer file from the feature class. This is where I start to get confused. I know UpdateLayer has to be used, and I have to base the update on a source layer? I tried this with no luck. If someone can just guide me on what order the steps are to apply a graduated color symbology to a feature class using arcpy.mapping, it will be appreciated. I should be able to figure out the rest of the scripting. Thanks!


title: Postgres/Postgis spatial index - no speed up

text: I have a spatial table in a postgres/postgis database. Every row in it represents a Polygon. It is of following form:               +----+--------+     |gid |   way  |     +----+--------+     |241 | 01030..|      The geometric column is "way" which contains the geometry for a polygon. In WKT it is: POLYGON(('....')). I am doing a lot of ST_Contains queries on this table, to test whether two polygons are contained in each other e.g.:               Select ST_Contains(a.way, b.way) From table AS a, table AS b Where a.gid = 15 And b.gid = 16      I was wondering how to speed up this query and added a spatial index on the table:               CREATE INDEX table_way_gist ON table USING gist(way);      But actually I dont see a speed up. I create the index AFTER I filled the table with all polygons BEFORE I do the ST_Contains queries. Should the index be added before filling a table? Are there special requirements on the table to work with the index? The projection (srid) of the geometric column way is set to 900913. I am using: psql (PostgreSQL) 9.1.4 / POSTGIS="1.5.3" Thx for help.


title: How/Where can I find vector data for building footprints in Texas?

text: I've done a lot of research on this topic over the past 8 or 9 months, but I'm a GIS newb, so I may have overlooked some possibilities for finding/extracting this data. A short list of resources I have researched (an exhaustive list would be very long):   1. Tax Net USA - Text data on properties is free. You have to pay to get footprint images. I don't mind paying, but their search interface is cumbersome, and the images they provide are small raster PNGs with 'cruft' that I don't need, such as dimensions. Additionally, each image is 'divorced' from the real world - there is no lat/lon or compass orientation information in the images. On the plus side, the images are detailed and are often accurate to within a foot or so.   2. ArcGIS - The Texas maps seem to have acceptable coverage and reasonably accurate building footprint polygons, but I have no idea how I would go from a street address to a vector graphic of the building(s) at the address using ArcGIS (or whether that's even feasible).   3. OpenStreetMap - Very sparse coverage of buildings in the central Texas area, the area I'm most interested in. I'm not frightened by a highly technical solution. I'm a software developer, so I'm willing to install, write and/or interface with software to do the job. Requirements in order of importance: Vector representation, accuracy, "queryability" (publicly available API would be great), lat/lon & compass orientation info. I'm OK with querying each address individually, on an as-needed basis, but a bulk dataset would also be interesting.


title: Select By Attribute using geoprocessing without using layers?

text: I am selecting features based on attributes within an ArcGIS Add-In. The selections are pretty simple(e.g. "DESCRIPTION = left"). Is there a tool or method of doing this that doesn't require passing layers around? Is there a simple ArcObjects method I may have looked over? After the selection I will need to pass the selected features into another geoprocessing step - so I need to get the output selection as a path (string). Currently I am: 1) Creating a feature layer from a given file. 2) Selecting features from the layer. 3) Saving the layer using Copy Features to get the output as a new shapefile.


title: ArcGIS 10.0/arcpy: Can't retrieve the points in a polyline geometry obect: any idea why?

text: My arcpy script builds a Python dictionary, dictCLPnts, that holds the geometry objects of a polyline feature class. The following code snippet shows the attempt to access the geometry object with the key "rwy":               feat = dictCLPnts[rwy]     arcpy.AddMessage('feat is type  %s' % (type(feat)))     arcpy.AddMessage('partcount = %s' % (str(feat.partCount)))     arcpy.AddMessage('pointcount = %s' % (str(feat.pointCount)))     part0 = feat.getPart(0)     arcpy.AddMessage('part0 is type  %s' % (type(part0)))     part0.reset()     pnt1 = part0.next()     if not pnt1:  #in case there is a Null pnt at the beginning of the part         pnt1 = part0.next()     arcpy.AddMessage('pnt1 is type  %s' % (type(pnt1)))     X1 = pnt1.X     Y1 = pnt1.Y      Here is the output:               feat is type  <class 'arcpy.arcobjects.geometries.Polyline'>     partcount = 1     pointcount = 1002     part0 is type  <class 'arcpy.arcobjects.arcobjects.Array'>     pnt1 is type  <type 'NoneType'>     <type 'exceptions.AttributeError'>: 'NoneType' object has no attribute 'X'      Looks to me like the points are there (and there isn't any reason for them not to be). I checked the original data, and the points are all there, z-values and all. It draws fine, and Check Geometry reports nothing amiss. I am just baffled by this. Yes, I could refactor and use a cursor (yet again), but that would be hideous at this point in the code. Dictionaries are also a lot faster. If anyone has any insight on this one, please let me know. Thank you, all.


title: schools: how are they tagged - with the same method all over the world?

text: are all schools tagged in the same way? or do some countries do use a specal tag-(system)? **update: thanks** to the responders: what does this mean https://taginfo.openstreetmap.org/tags/amenity=school Typ countl Objekte all: 592 443 0.02% [Node] Node 325 974 0.36% [Way] Way 260 384 0.11% [Relation] Relation 6 085 0.23% end of **update: thanks** we have the central ressource for the tag-system http://wiki.openstreetmap.org/wiki/Tag:amenity%3Dschool > how can we evaluate that? > > This area should envelop the full grounds of the school including all the > buildings, sports facilities and grounds. For schools with multiple sites > the multipolygon relation can be used. Tag the element with amenity=school > and name=* for the name of school. Suggested additional tags include: >   >   >     operator=* - Name of operator, often the local education authority. >     addr=* - Address details >     capacity=*, for the number of pupils taught at the school >     isced:level=*, for the educational level (proposed tags) >     fee=yes if the the school makes a direct charge for core services. >     religion=*, if the school is associated with a particular religion (also > denomination=*) >     wikipedia=*, for a link to a Wikipedia article about the school >     website=*, for a link to the school's own website >   >      Public-images-osm logo.svg building = school Convent School's > Pippal.JPG Description A generic school building. Used on these >   > > elements may be used on nodes should not be used on ways may be used on > areas may be used on relations Useful combination >   >   >     name=* >     amenity=school >     amenity=college >     amenity=kindergarten >     See also >     building=university >   How to map > there is just one school building: Map the outline of the area Area of the > building and add the tags amenity=school and building=school. Of course you > can add name=* and address as well. Well - i wonder how we can get a concrete and concise result if we do a request for all the schools in the world               <!--     This query looks for nodes, ways and relations      with the given key/value combination.     Choose your region and hit the Run button above!     -->     <osm-script output="json" timeout="25">       <!-- gather results -->       <union>         <!-- query part for: “building=school” -->         <query type="node">           <has-kv k="building" v="school"/>           <bbox-query {{bbox}}/>         </query>         <query type="way">           <has-kv k="building" v="school"/>           <bbox-query {{bbox}}/>         </query>         <query type="relation">           <has-kv k="building" v="school"/>           <bbox-query {{bbox}}/>         </query>       </union>       <!-- print results -->       <print mode="body"/>       <recurse type="down"/>       <print mode="skeleton" order="quadtile"/>     </osm-script>        btw - this above code does not function - and this below also not. what is wrong here?               This query looks for nodes, ways and relations      with the given key/value combination.     Choose your region and hit the Run button above!     -->     <osm-script output="json" timeout="25">       <!-- gather results -->       <union>         <!-- query part for: “building=school” -->         <query type="node">           <has-kv k="building" v="school"/>           <bbox-query {51.5557914,0.2118915,51.5673083,0.2369398}/>         </query>       <!-- print results -->       <print mode="body"/>       <recurse type="down"/>       <print mode="skeleton" order="quadtile"/>     </osm-script>      can you help out!?


title: Infinite plane map using Polymaps library

text: I'm working on a web project for which I need a simple square-tile map which is infinite in all directions emanating from (0,0). From my research the Polymaps library looks like a good fit, and I've been basing my experiment on the Tile Grid example (http://polymaps.org/ex/tiles.html). That example code is very close to what I'm looking for, except the map is only infinite on the longitudinal axis; the latitude range is fixed at +/- 85 degrees. Is there a way I can remove this restriction (eg: by turning off it's use of the spherical mercator projection and just use the tile indexes?), or will I have to modify the source code for Polymaps to achieve that? Any guidance or suggestions would be greatly appreciated.


title: multiuser testing in web gis

text: We have build a web based application using aspmap and postgresql database (spatial database) While we were checking our app for multiuser testing, we found that the results are overlapped for the parameters selected by both users. The method we have used for rendering is - we have one column in database and we are updating the same and based on that field we are doing color rendering on map. Now if one user selects product x and another product y and both of them fires same query at same time. both of them gets mix results on map. as the same field in spatial table is getting update immediately one after another. IS this right approach for any web application or are there any alternatives which can control such situation?


title: problem in scale of projected raster

text: I have a shapefile (my basin limit) and one raster file that I projected it to the shapefile coordinate system and set its cell size to 1000m(area=1Km2). my basin area is 5400 Km2 and it should contain almost 5000 cell of raster file at least. but just one and a half cell put on this range. what's the problem? ![enter image description here](http://i.stack.imgur.com/iTSwa.jpg) ![enter image description here](http://i.stack.imgur.com/LA6Nm.jpg) ![enter image description here](http://i.stack.imgur.com/ac0Sa.jpg) ![enter image description here](http://i.stack.imgur.com/ijY1N.jpg)


title: Capturing GPS coordinates on the ground in Boliva?

text: Assuming no cellular or internet connection, what is the best way to capture GPS coordinates on the ground in Boliva? Preferably a free approach. For example, are their any android apps that can do this without a cellular connection?


title: How to build SQL Expressions in ArcGIS 10.1 Python script tool dialog?

text: I've made a python tool that requires the user to establish which state to select within an input's table (the "Prepared Footprint"). However, once the user drops in the layer and tries to write the SQL statement to select the state which is contained in the state field of the Prepared Footprint, there are no available fields to select from. ![enter image description here](http://i.stack.imgur.com/58nLz.png) Is there some option within the properties of the python script tool that will allow the Query Builder to populate the fields and attributes from the input layer like the native Select by Attributes in ArcMap? If not, how might I code this into my source script to feed fields and their attributes to the Query Builder of the "State" parameter? (experience level with building tools: Beginner) * * * EDIT: Here's a screenshot of the tool's parameters. When I select "State" and try to change the "Obtained from" paramter property nothing happens. There's no drop down, no cursor, nothing but a highlighted row for "Obtained from". Has anyone ever had this problem before? ![enter image description here](http://i.stack.imgur.com/AoSZ0.png) * * * EDIT: Found the problem. A "Layer" data type cannot be used as the input for "Obtained from". Feature Layer works instead for this case. ![enter image description here](http://i.stack.imgur.com/JF1bC.png)


title: Map services error in Firefox bug but don't understand

text: I am currently harvesting many web map services on an online map view. The services are configured in Oracle (Toad). Most of the maps seem to be working apart from 10 or so. I can see them in the ToC but no map layer or corresponding legend is showing. In firefox I get the following error: TypeError: mapService is undefined I've opened all the services in ArcGIS without any problem, checked I'm using the right xml source files and if they all correspond correctly with the Oracle database and it's fine. Can anyone help me please???


title: Use Gdal2tiles to generate 512px tiles

text: I'm a newbie in gdal2tiles & python but I would like to modify gdal2tiles to create 512px tiles for retina display (screen). Can someone suggest how I can do this?


title: Connecting QGIS 2.x to Firebird Database via ODBC?

text: I tried to connect a Firebird Database in QGIS 2.2 and 2.4 via ODBC to a database and get following error message: > Unable to initialize ODBC connection to DSN for SYSDBA/masterkey@Akdees > [Microsoft][ODBC Driver Manager] The specified DSN contains an architecture > mismatch between the Driver and Application **Funny, that it does work with QGIS 1.8!** but fails with QGIS 2.2 and 2.4. Does anyone have an idea as to what the reason is and how to make it work again?


title: TypeError with ArcGIS Flex: cannot convert views

text: I am new to Flex, and Flash builder, I followed the tutorial in this link: https://developers.arcgis.com/flex/guide/tutorial-query-your-data.htm to write a mobile mapping application that includes a query task. After writing the code, no errors appeared in the Problems tab but there were 6 warnings. And when running the app, the following errors appeared (also shown in screenshots):               An ActionScript error has occurred:     TypeError: Error #1034: Type Coercion failed: cannot convert views::AGISHomeView@6396ae1 to spark.components.View.     at spark.components::ViewNavigator/createViewInstance()[E:\dev\4.y\frameworks\projects\mobilecomponents\src\spark\components\ViewNavigator.as:1999]     at spark.components::ViewNavigator/commitNavigatorAction()[E:\dev\4.y\frameworks\projects\mobilecomponents\src\spark\components\ViewNavigator.as:1932]     at spark.components::ViewNavigator/commitProperties()[E:\dev\4.y\frameworks\projects\mobilecomponents\src\spark\components\ViewNavigator.as:1301]     at mx.core::UIComponent/validateProperties()[E:\dev\4.y\frameworks\projects\framework\src\mx\core\UIComponent.as:8219]     at mx.managers::LayoutManager/validateProperties()[E:\dev\4.y\frameworks\projects\framework\src\mx\managers\LayoutManager.as:597]     at mx.managers::LayoutManager/doPhasedInstantiation()[E:\dev\4.y\frameworks\projects\framework\src\mx\managers\LayoutManager.as:783]     at mx.managers::LayoutManager/doPhasedInstantiationCallback()[E:\dev\4.y\frameworks\projects\framework\src\mx\managers\LayoutManager.as:1180]”      * * * Here is the code I used:               <?xml version="1.0" encoding="utf-8"?>     <s:Application xmlns:fx="http://ns.adobe.com/mxml/2009"                          xmlns:s="library://ns.adobe.com/flex/spark"                          xmlns:mx="library://ns.adobe.com/flex/mx"                          xmlns="http://ns.adobe.com/air/application/4.6"                          xmlns:esri="http://www.esri.com/2008/ags">                <s:layout>                 <s:VerticalLayout gap="10"                                           horizontalAlign="center"                                           paddingBottom="20"                                           paddingLeft="25"                                           paddingRight="25"                                           paddingTop="20"/>           </s:layout>           <fx:Script>                 <![CDATA[                       import com.esri.ags.FeatureSet;                       import mx.controls.Alert;                       import mx.rpc.AsyncResponder;                            private function doQuery():void                       {                             queryTask.execute(query, new AsyncResponder(onResult, onFault));                                  function onResult(featureSet:FeatureSet, token:Object = null):void                             {                                   // No code needed in this simple sample, since the                                   // graphics layer is bound to the query result using                                   // graphicProvider="{queryTask.executeLastResult.features}"                             }                                  function onFault(info:Object, token:Object = null):void                             {                                   Alert.show(info.toString(), "Query Problem");                             }                       }                 ]]>           </fx:Script>           <fx:Declarations>                 <!-- Layer with US States -->                 <esri:QueryTask id="queryTask"                                         showBusyCursor="true"                                         url="http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/5"                                         useAMF="false"/>                 <esri:Query id="query"                                   outSpatialReference="{myMap.spatialReference}"                                   returnGeometry="true"                                   text="{qText.text}">                            <esri:outFields>                             <fx:String>MED_AGE</fx:String>                             <fx:String>POP2007</fx:String>                       </esri:outFields>                 </esri:Query>                </fx:Declarations>             <s:Panel height="60"                        backgroundColor="0xB2BFC6"                        title="Query a layer (search for a state)">                      <s:layout>                       <s:HorizontalLayout/>                 </s:layout>                      <s:TextInput id="qText"                                    width="100%"                                    enter="doQuery()"                                    text="California"/>                      <s:Button click="doQuery()" label="Do Query"/>           </s:Panel>           <esri:Map id="myMap">                 <esri:extent>                       <esri:Extent xmin="-14298000" ymin="2748000" xmax="-6815000" ymax="7117000">                             <esri:SpatialReference wkid="102100"/>                       </esri:Extent>                 </esri:extent>                      <esri:ArcGISTiledMapServiceLayer                       url="http://server.arcgisonline.com/ArcGIS/rest/services/World_Physical_Map/MapServer"/>                </esri:Map>     </s:Application>


title: Cant find Define Line Portion tool

text: I am following the procedure on this help page: http://resources.arcgis.com/en/help/main/10.1/index.html#//00390000002z000000 However when at step 3, I can't seem to find the _Define Line Portion_ tool. Can someone please help locate that tool? I've also posted the picture of how my interface looks like.![enter image description here](http://i.stack.imgur.com/b1Z1u.png)


title: How to save MakeNetCDFTableView output?

text: I was trying to extract the time-series data from the netcdf using grid points (lat/long) by `MakeNetCDFTableView`. I have more than 100 grid points and I want to extract the data from the netcdf file. I was trying the following pyscript, but couldn't save the output table. Could you please help me fix my problem?               #Import system modules     import arcpy          # Set local variables     inNetCDFFile = "F:/VIC_FORCINGS/7071.nc"     variable = "prec"     outTableView = "F:/VIC_FORCINGS/precipm"     rowDimension = "time"     dimensionValue = "64;-146"     valueSelectionMethod = ""          # Execute MakeNetCDFTableView     arcpy.MakeNetCDFTableView_md(inNetCDFFile, variable, outTableView, rowDimension, dimensionValue,valueSelectionMethod)          outTable = arcpy.CopyRows_management(outTableView)     print


title: Colour issue when reprojecting image mosaic

text: We have a couple of instances of GeoServer running, one of which is the latest version and the other is version 2.1.3. The older version has Ordnance Survey base mapping which is loaded as a number of EPSG:27700 image-mosaics. These are then reprojected by GeoServer into the Google projection EPSG:900913 and displayed as part of a layer group where the style switches the layers on / off depending on the scale. We believe this is a fairly standard use of GeoServer to provide a base map coverage and it works fine. In the newest versions of GeoServer the same base map layer creates a horrible orange coloured map as it seems to be breaking the colour banding of the layers. Seems to be turning the large swaths of blue sea into orange, as well as screwing the colours on the rest of the map. I can't upload a screenshot as our organisations internet policies restrict it. Please note, you need to reproject two image-mosaics and combine them in a layer group or GetMap to reproduce the error. We wondered if anyone else has experienced the bug in the latest versions, and whether there are any solutions or workarounds anyone has come across.


title: How to use DotSpatial to project/reproject shapefiles on the fly

text: I would like to use DotSpatial to reproject shapefiles on the fly before loading them to GeoServer. I've been looking for tutorials on how to get started with this in vain. Any pointers on where I can start will be highly appreciated.


title: How to break a single polygon into multiple polys?

text: I created voronoi polygons after which I dissolved them on the basis of one attribute. Now there are cases where in for one single attribute "abc" I have3 polygons at different places with same object id. I would like to break these into 3 polygons as they are at different locations. Any tool in ArcGIS or QGIS which can help me in doing that. Thank You.


title: How to fix Error 14 -- Unable to open ~/.qgis//qgis.db?

text: I get this error: Error 14 -- Unable to open ~/.qgis//qgis.db -- QGIS version 1.8.0-Lisboa -- Ubuntu 12.04 It was discovered when I tried to define a custom CRS and save it under: settings/custom CRS in QGIS 1.80.-Lisboa under Ubuntu 12.04 and saving the custom CRS indicated this error. I notice that the qgis.db file is located at: /usr/share/qgis/resources/qgis.db I created a link but this did not work either. Can I tel QGIS to look for the qgis database at the above location rather that ~/.qgis//? Looks like permissions are read and write for all users. Thanks


title: OpenLayers on iPhone with navigation: marker-click event won't work

text:
I'm making a map with openlayers. For the navigation on the iPhone I added OpenLayers.Control.Navigation(), so zoom and pan work fine. On the desktop, i can click a marker and a popup appears. But on the phone nothing happens? Any suggestions?               map = new OpenLayers.Map("map", {projection: mercator,         controls: [                            new OpenLayers.Control.LayerSwitcher(),             new OpenLayers.Control.Navigation(),             new OpenLayers.Control.ZoomPanel(),             new OpenLayers.Control.MouseDefaults()                  ],             units: 'meters',     });                    temp = tCoord(lon, lat);             var feature = new OpenLayers.Feature(markerLayer, temp);             feature.closeBox = true;             feature.popupClass = OpenLayers.Class(OpenLayers.Popup.Anchored, {minSize: new OpenLayers.Size(320, 100)} );             feature.data.popupContentHTML = "
" + name +"
" + desc                                              +"
Route hierher
";             feature.data.overflow = "hidden";                  var tempMarker = new OpenLayers.Marker(temp, icon.clone());             tempMarker.feature = feature;              tempMarker.events.register("mousedown", feature, markerClick);             allPoints.push(new Array(zeilenCounter, tempMarker));      (tCoord transforms from wgs84 to spherical mercator)


title: How to query the duplicated field of QGIS?

text: I have a column which have many duplicated field record. How can I find those duplicated features? Please suggest.


title: List of Virtual Globes

text: I kept a listing of virtual globes in 2005 and 2006, but it's now woefully out of date http://geography2.blogspot.com/2005/11/listing-of-virtual-globes.html. Today, the main four seem to be: Google Earth (proprietary and free), Microsoft Bing Maps 3D (proprietary and free), ESRI ArcGIS Explorer (proprietary and free), and NASA World Wind (free and open source). Many, if not most, are being integrated into the web browser. What is out there? Any special or particularly innovative capabilities?


title: OpenLayers Point Feature Vanish on Zoom-In

text: I use OpenLayers to animate a point moving along a line. To do this I use several map layers and one vector layer to contain one LineString feature and one Point feature. Currently, the animation is working but I have one issue where every time I zoom in the map, the Point feature vanish even though the LineString is still shown. When I zoom out, the Point feature is shown again but when I try to zoom in again, the Point vanished. These are screenshots of what happen ![Default Zoom](http://i.stack.imgur.com/zIhzO.jpg) ![Try to zoom in, point vanished](http://i.stack.imgur.com/MJJIz.jpg) ![Zoom out again, point returns](http://i.stack.imgur.com/Kiokm.jpg) This is a portion of source code I use to animate the Point:               /* == Fitur Simulasi Gerak di Rute */     var themarker = new OpenLayers.Feature.Vector(         new OpenLayers.Geometry.Point(map.getCenter().lon, map.getCenter().lat),         null,         {externalGraphic: imgUrl+"symbol/perusak kawal.png", pointRadius: 10, rotation: -90}     );     var linePath;     var indexPoint = 0;     var bergeraklah = function() {         var pindah = function(dari, ke, asal) {             var jumlah = 10; // hardcoded             var xmove = Math.abs(asal.x - ke.x)/jumlah;             var ymove = Math.abs(asal.y - ke.y)/jumlah;             // x-axis             if (dari.x < ke.x) {                 if (dari.x + xmove < ke.x) {                     dari.x += xmove;                 } else {                     dari.x = ke.x;                 }             } else {                 if (dari.x - xmove > ke.x) {                     dari.x -= xmove;                 } else {                     dari.x = ke.x;                 }             }             // y-axis             if (dari.y < ke.y) {                 if (dari.y + ymove < ke.y) {                     dari.y += ymove;                 } else {                     dari.y = ke.y;                 }             } else {                 if (dari.y - ymove > ke.y) {                     dari.y -= ymove;                 } else {                     dari.y = ke.y;                 }             }         };              var interval = 100; // hardcoded         if (themarker.isMoving) {             var tujuan = linePath.geometry.components[indexPoint];             if (themarker.geometry.x != tujuan.x || themarker.geometry.y != tujuan.y) {                 pindah(themarker.geometry, tujuan, linePath.geometry.components[indexPoint-1]);                 g.ruteLayer.drawFeature(themarker);                 window.setTimeout(bergeraklah, interval);             } else {                 themarker.isMoving = false;                 bergeraklah();             }         } else if (indexPoint < (linePath.geometry.components.length-1)) {             var prev = linePath.geometry.components[indexPoint];             // fetch next point             var tujuan = linePath.geometry.components[++indexPoint];             // hitung kemiringan             var bearing = tc1(prev.x, prev.y, tujuan.x, tujuan.y);             themarker.isMoving = true;             themarker.style.rotation = bearing - 90; // gambarnya hadap kanan 90 derajat             pindah(themarker.geometry, tujuan, linePath.geometry.components[indexPoint-1]);             g.ruteLayer.drawFeature(themarker);             window.setTimeout(bergeraklah, interval);         } else {             // selesai             g.simulating = false;             g.ruteLayer.removeFeatures([themarker]);         }     };               $('#button-simulasi').bind('click', function(e) {         if (!g.simulating) {             g.simulating = true;             indexPoint = 0;             linePath = g.ruteVector.clone();             themarker.geometry.x = linePath.geometry.components[0].x;             themarker.geometry.y = linePath.geometry.components[0].y;             g.ruteLayer.addFeatures([themarker]);             themarker.isMoving = false;             bergeraklah();         }     });     /* == End of Fitur Simulasi Gerak di Rute */      Does anybody know what cause the vanishing Point? Because I tried to google it but still don't find any solution.


title: How can I retrieve an attribute's domain using ArcIMS ActiveX Connector?

text: I have an ArcIMS Image Service that's consuming data stored in an ArcSDE database (both v9.3.1). Many of the layers have attribute columns that use coded value domains defined in the SDE database, and I'd like to be able to display the values defined in the domains. Is it possible to retrieve the all values in the domain associated with an attribute column via an ArcXML request to the ActiveX connector?


title: What's the fastest way to do a bounding box query in PostGIS

text: I have a Postgres table, with almost 2 million rows, with a long-lat `coordinates` field in the form `POINT(-73.4938 33.2405)`. Supposing there's a geospatial index on that field, what's the most efficient, fastest way to select all the rows within an arbitrary bounding box? The box is like `SW long-lat: -74.0042 40.7688`, `NE long-lat: -73.8809 40.7984`.


title: Finding measurements of unknown points relative to reference points with known lat, long and measurements?

text: I am trying to measurements of the distance of multiple points relative to reference points of known latitude, longitude and measurements along a river. Is there tool that can do this? Would dynamic segmentation or linear reference work to able to see the measurements in the attribute table?


title: How to use Openlayers TMS with another SpatialReference?

text: I am trying to use OpenLayers with my TMS (tilecache). It fetches tiles well but tiles are not in right order. I use EPSG:102067 - Krovak East-North, which has negative axis. Maybe because of this it is flipped. ![look at image](http://i42.tinypic.com/2mcy51l.jpg) Is there some solution, how make it work?!                   <script type="text/javascript">         var lon = -651799;         var lat = -1127290;         var zoom = 5;         var map, layer;              function init(){             map = new OpenLayers.Map( 'map', {maxResolution:1.40625/2} );             layer = new OpenLayers.Layer.TMS( "TMS",                     "http://localhost/cgi-bin/tilecache-2.11/tilecache.cgi/", {layername: 'osm', type:'png', serviceVersion: '1.5.0'} );             map.addLayer(layer);             map.addControl(new OpenLayers.Control.LayerSwitcher());             map.setCenter(new OpenLayers.LonLat(lon, lat), zoom);         }          </script>      Thank you


title: consecutive generate Sql query using script arcpy

text: I am using arcgis 10.1 I am currently doing a lot of sql query feature class and I am thinking of automating this process so far my escrip is as follows               import arcpy     arcpy.SelectLayerByAttribute_management("casa","NEW_SELECTION"," [color] = 'blanco' ")      then I get a count of the selected               contador = int(arcpy.GetCount_management("casa").getOutput(0))     print contador      > 88 then do another query on the same feature class and I get him the count and repeat this process as many times on the same feature class conducted many consultations my first question is how could I modify this scrip to use and able to put in modelbuider as input parameter one feature class and make all of these query as defined in escrip my second question if I could list all the results it gives me the count of each query and export them to a table


title: how to set width & height of a polygon in QGIS

text: I can't setting width and height of a polygon How to do that in QGIS? Thank you


title: MapBasic error when creating button with conditional functions: "Commandinfo: Argument 2 out of range"

text: I am new to MapBasic and this particular forum so please excuse my somewhat simple questions and potential violations of the GIS SE protocol - not intentional :) Also, thank you in advance for any help. I ran into this old thread regarding buttons with conditional functions and tested the code suggested by Peter Horsbøll Møller there: Making a Button in mapinfo with conditional function However I ran into the following error message "Commandinfo: Argument 2 out of range". I've seen the same issue reported on other discussion forums with no outcome. Has anybody seen/resolved this issue before? I'm using MapBasic 10.0/MapBasic IDE/MI Professional 12.0. Thank you!


title: Export MXDs files to PDFs using Python

text: I'm trying to write a python script that will export multiple MXDs to individual PDFs while using a feature class field to rename the pdfs (this can be an optional field, if not check it can copy the mxd file name). I copied a script from ESRI, just need the renaming of pdfs part to work. import arcpy, os, string               #Read input parameters from script tool     mxdList = string.split(arcpy.GetParameterAsText(0), ";")     outPDFpath = arcpy.GetParameterAsText(1)          #Create a new PDF object to store the results     outputPDF = arcpy.mapping.PDFDocumentCreate(outPDFpath)


title: Creating a polyline from point data that represents multiple linear features in ArcGIS 10.0

text: I have a problem where I have to somewhow create a polyline feature class in ArcGIS 10.0 from points that were collected in the field using a Trimble GeoXT that represent the start and end of multiple linear features such as stone walls, wire fences, etc... The Trimble GeoXT collects the date and the exact time that the point for the start and end of the linear feature was collected. The idea is I am walking a forest boundary line where sometimes there is a stone wall that starts and follows the boundary, but then stops and there is nothing for 20 feet or so and then a wire fence starts and follows the boundary and may end or transition into a stone wall, something else or nothing at all along this forest boundary line. I need a python script that can draw a line from the start of one linear feature, such as a stone wall, to its end point and stop and then do the same thing for all the other types of linear features (stone walls, wire fences, barb wire fences and virginia rail to be exact). Is there any way this is possible? I am not very skilled using Python scripting. The Points To Lines tool in ArcGIS 10.0 only draws a single line through every single point regardless of what type of linear feature the point represents and continues this single line where there is no fence. The Points To Lines tool also does not cary over any of the attribute data that is associated with the original point layer. I hope I have provided enough information. Please ask for more information if needed.


title: change xslt to match all OSM tags

text:
Hello and good day, I am brand new to GIS and to this forum so please bear with me for asking a newbie question. I have an xslt that works on a split osm file[germany.osm.bz2] What I want to do is parse all elements of `'amenity'` and `@v='restaurant`. In other words i want to gather all "available" elements (nodes, ways and relations) to match the elements I need to rework the given filter below so that we can get all elements.... See the current XSLT -               <xsl:stylesheet version = '1.0'             xmlns="http://www.w3.org/1999/xhtml"             xmlns:xml_split="http://xmltwig.com/xml_split"             xmlns:xsl='http://www.w3.org/1999/XSL/Transform'>              <xsl:output method="text" encoding="UTF-8"/>         <xsl:template match="/">                      <xsl:for-each select="xml_split:root/node/tag[@k='amenity' and @v='restaurant']">                 <xsl:value-of select="../@id"/>                 <xsl:text>	</xsl:text>                 <xsl:value-of select="../@lat"/>                 <xsl:text>	</xsl:text>                 <xsl:value-of select="../@lon"/>                 <xsl:text>	</xsl:text>                 <xsl:for-each select="../tag[@k='name']">                     <xsl:value-of select="@v"/>                 </xsl:for-each>                 <xsl:text>
</xsl:text>             </xsl:for-each>         </xsl:template>          </xsl:stylesheet>      it gives back the following results               "27376049   49.6499423  8.5543447   Clubhaus SC Olympia Lorsch"     "27407910   49.7045590  8.5115993   Forsthaus Jägersburg"     "27433397   51.4588114  7.0073003   Unperfekthaus"     "27433643   50.0081188  9.0691993   Schützenhaus"     "27453754   49.6358087  8.4547972   Feldschlössel"     "27474067   50.1880211  9.1742722   Pizzeria A66"     "27475403   49.0062148  8.3369945   Marty's Restaurant"     "27505752   48.4081456  11.7959129  Gasthof Nagerl"     "27530181   49.5054702  8.7089191   Zum Pflug"      I want to have a more in depth results with all tags. What modifications do I need to make to the xslt-file? I hope that I was able to ask the question correctly. If you have any questions do not hesitate to ask.. **Update:** moved the question - thx to Stev_K. Thx for all your help here. Guess that the question has a far more appropiate place at Stackoverflow!


title: How is the "place" determined for this Tiger 2010 shapefile address record?

text: In most cases, an ADDR.DBF record links to an EDGES.DBF record that links to a FACES.DBF record that links to a PLACE.DBF record that contains the place. Sometimes, the place key in FACES.DBF is blank, the COUSUB key can be used instead, and COUSUB.DBF contains the place. But what do you do in the following case?               TL_2010_2047_ADDR.DBF contains an address range with     primary key: ARID = 4003948456882     TLID = 85556021     SIDE = R      which links to edge:               TL_2010_22047_EDGES.DBF     primary key: TLID = 85556021     TFIDR = 201079140      which links to face:               TL_2010_22047_FACES.DBF     primary key: TFID = 201079140     COUSUBFP10 = 94267     PLACEFP10 = blanks      which has only a COUSUB key, not a PLACE key, but unfortunately the NAME10 field in the TL_2010_22_COUSUB.DBF file for key 94267 contains "2", which is not a viable place name (the correct place is probably something like "White Castle"). What are the rules by which COUSUBFP10 and/or PLACEFP10 are either used or ignored, how is the correct place name supposed to be determined, and where are those linkage rules documented?


title: How to let ArcGIS Server see related records in a geodatabase

text: I have an ArcSDE geodatabase which contains multiple relationship classes. If I add this layer to ArcMap, the relationship is automatically applied (eg, the Identify dialog allows me to drill-down to the related features). However, if I publish this map document to ArcGIS Server, the REST endpoints for the layer don't show the relationship classes. Similarly, if I create a new ArcMap document and create a Relate within ArcMap (using > Joins and Relates in the TOC) this also fails to appear in ArcGIS Server's REST directory. This sample REST directory does show the related data. How has this been set up?


title: How to use a field-to-RGB mapping for symbology in QGIS?

text: Using QGIS version 1.7. I have a plain text file that lists a set of rgb values against a code. I want to use this colour table to colour a polygon layer by mapping one of its attribute fields ('map_symb') to a code in the text file. the colour table is very long, and looks like this:               $ head gsv1Msymbology.txt     MAPCODE RED GREEN   BLUE     Oc  143 255 255     WAT 255 255 255     Qa  244 250 202     Qdl 195 239 218     Na  248 255 238     Qd2 227 255 190     Qxw 248 255 238     Qns 255 148 83     Qn  255 202 190     ....       So I want to match my 'map_symb' attribute to a value in MAPCODE, and use the corresponding RGB values to colour the polygons. Is there a gui way to do this? Any way? Ben.


title: Can use of GML Simple Features profile by Web Feature Service (WFS) be advertised?

text: The ArcGIS 10.1 documentation on Working with WFS services says: > The GML Simple Features profile is recommended for use with ArcGIS and is > the basis for WFS support in ArcGIS. Is there an easy way to determine whether a particular WFS is encoding using that or an alternative profile - perhaps by looking at its GetCapabilities? It's interesting to note that not a single service URL was suggested when I asked whether anyone knew of a public WFS using the Simple Features profile. UPDATE: I'm also trying to uncover an answer to this via the ArcGIS Discussion Forums now.


title: How can I estimate kernel density, taking barriers (i.e. a priori zero-probability regions) into account?

text: How can I estimate kernel density, taking barriers (i.e. a priori zero- probability regions) into account? At the moment, regions where the population cannot occur are clipped out for display. However, I imagine this is biasing the results, since observations near the zero-probability regions don't have as much impact as they might. I have several thousand observations, so a tool or programmatic solution is a necessity. Bonus points for an implementation in ArcGIS or R.


title: Failed to open raster - ArcGIS 10 SP2, Python 2.6.5, Windows 7

text: I have an arcpy / python script that is working fine with ArcMap 10 no Service pack (Build 2414) on windows XP, but consistently failing on using ArcMap 10 SP2 (Build 3200). Getting the ESRI 99999 error. In a nutshell, this script is supposed to reclassify and convert four raster files to polygon (shp). Here is the code. Any help or ideas would be appreciated.               import arcpy, os, string     arcpy.CheckOutExtension("Spatial")     from arcpy.sa import *          #set output location     outLoc = "O:/gis/!KEYHOLE/DAILY_SHP"          #get file date, set data directory     fileDate = arcpy.GetParameterAsText(0)     dataFolder = "O:/gis/PRED_SVC/MET_PROJ/FOR_OBS_LTG_DATA"          #add file date message     arcpy.AddMessage("...using file date " + fileDate + "...")          #set precip amount file based on date input     prefixList = ["erc", "thz", "prpdur_", "prpamt_"]     for prefix in prefixList:     metName = prefix + fileDate     metFile = os.path.join(dataFolder, metName)     print metName     print metFile          #output workspace, turn overwrite output on     arcpy.env.workspace = "c:/tmp"     arcpy.env.overwriteOutput = 1          #reclassify     if prefix == "erc":         remap = RemapRange([[-99, -1.282, 1], [-1.282, -0.842, 2], [-0.842, -0.524, 3], [-0.524, -0.253, 4], \                  [-0.253, 0, 5], [0, 0.253, 6], [0.253, 0.524, 7], [0.524, 0.842, 8], [0.842, 1.282, 9], \                  [1.282, 1.645, 10], [1.645, 1.881, 11], [1.881, 99, 12]])         arcpy.AddMessage("...reclassifying ERC...")          elif prefix == "thz":         remap = RemapRange([[-99, -1.645, 1], [-1.645, -1.282, 2], [-1.282, -0.842, 3], [-0.842, -0.524, 4], \                  [-0.524, -0.253, 5], [-0.253, 0.253, 6], [0.253, 0.524, 7], [0.524, 0.842, 8], \                  [0.842, 1.282, 9], [1.282, 99, 10]])         arcpy.AddMessage("...reclassifying 1000HR fuels...")          elif prefix == "prpdur_":         remap = RemapRange([[0, 0.009, 1], [0.009, 1.9999, 2], [1.9999, 4.9999, 3], \                             [4.9999, 9.9999, 4], [9.9999, 25, 5]])         arcpy.AddMessage("...reclassifying precipitation duration...")          else:         remap = RemapRange([[0, 0.0009, 1], [0.0009, 0.09999, 2], [0.09999, 0.249999, 3], \                             [0.249999, 0.499999, 4], [0.499999, 0.999999, 5], \                             [0.999999, 1.999999, 6], [1.999999, 1000000, 7]])         arcpy.AddMessage("...reclassifying precipitation amount...")          rcls = Reclassify(metFile, "Value", remap)           #make raster layer, convert to poly     metLyr = arcpy.MakeRasterLayer_management(rcls, "metLyr")     arcpy.env.workspace = outLoc     arcpy.env.overwriteOutput = 1     outPolyName = string.replace(prefix, "_", "") + "_poly.shp"     outPolyFile = outLoc + os.sep + outPolyName     arcpy.AddMessage("...converting " + metFile + " to poly...")     print outLoc     print outPolyName     print outPolyFile     arcpy.RasterToPolygon_conversion(metLyr, outPolyFile)     del rcls          arcpy.CheckInExtension("Spatial")


title: Connecting ArcSDE to Quantum GIS for dummies?

text: I really have a need to be able to connect to and edit ArcSDE data in Quantum GIS. This question has been asked before and apparently is possible, but the instructions are just not specific enough for me to follow... Is it possible to connect to ArcSDE using QuantumGIS? https://hub.qgis.org/wiki/17/Connecting_to_ArcSDE_databases I was able to get the sde-9.3/ogr_SDE.dll from the MapServer zip download and place it in C:\OSGeo4W64\bin\gdalplugins, but then what? I've read I need to compile it, but how and with what? I've read I need to create a vrt file. Is there a tool that creates a vrt file or should I just open notepad and write one. Then how do I add the vrt file to Quantum GIS? Thanks for any help.


title: QGIS - exercise data

text: I´m working currently through the QGIS training manual. However, the data for the spatial analysis exercise in chapter 9.4 (http://docs.qgis.org/2.2/en/docs/training_manual/complete_analysis/analysis_exercise.html) is not available under the link provided (can´t find the archive "exercise_data/more_analysis"). Somebody has the data at hand or can tell me where to download the archive. thank you & best regards, Flo


title: count attribute in raster

text: I'm trying to calculate the total number of counts a raster contains using python script. When I ran the script, I didn't get any error. However, I also didn't get the number of counts for each row in the attribute table. Any suggestion/thought would be greatly appreciated! Code:               import arcpy          #To overwrite output     arcpy.env.overwriteOutput = True          #Set environment settings     arcpy.env.workspace = "C:/Subhasis/Test/raster-read"          #Build Raster Attribute Table     arcpy.BuildRasterAttributeTable_management("sti-9", "Overwrite")          # loop through attribute table     Rows = arcpy.SearchCursor( "sti-9" )           for row in Rows:         count = row.getValue( "COUNT" )          print count


title: Linking data to road colour scheme

text: I'm getting started with a qgis map showing the quality of roads in my town. I have various data which create an index of 1-5. How do I link this index data to road data in a map, showing the bad roads (5) in red, and good roads (1) in green? What is the correct term for this process? Any pointers on this would be greatly appreciated. Regards, Case.


title: How to add map layers from Oracle database VIEWS in MapServer mapfile?

text: I am using MS4w v 3.0.3 (with Java mapscript) with Oracle Spatial 11g. Can we add layers to MapServer mapfile from **_oracle database views_**? If yes, is there any specific way to create views with spatial column for mapserver mapfile? I am getting the following error when I try to add a layer from a view : `Exception in thread "http-bio-8080-exec-1" java.lang.UnknownError: msDrawMap(): Image handling error. Failed to draw layer named 'amplifier'.;msOracleSpatialLayerWhichShapes(): OracleSpatial error. Error: ORA-13226: interface not supported without a spatial index`. Regards, Shiva


title: arcpy Tool works in debug mode but not in normal model

text: When I try to run the following simple script from arcmap it fails but when run using right--click debug and then F9 (debug) in python scripter it works fine. I don't understand why??? ![enter image description here](http://i.stack.imgur.com/AVkON.png)               # Modified by: Georgec Corea (ATGIS) coreagc@gmail.com          # Author:  ESRI     # Date:    July 5, 2010     # Version: ArcGIS 10.0     # Purpose: This script will print one or more map documents to a local printer.     #          The script is intended to run within a script tool.  There are two     #          parameters:     #               1) Select Map Documents to Print,     #               2) Select Output Printer (auto populated using a validation script)     #     #Notes: The print order of the MXDs is based on how they are entered.  The MXD     #       at the top of the list is first followed by those below it.          import arcpy, string, datetime, shutil, os     import arcpy.mapping as MAP          #Read input parameters from script tool     MXDList = string.split(arcpy.GetParameterAsText(0), ";")     #printer = arcpy.GetParameterAsText(1)          count=0     #Loop through each MXD and print     for MXDPath in MXDList:         count=count+1         arcpy.AddMessage(str(count)+'. Working on '+str(MXDPath)+' @ '+str(datetime.datetime.now())+"\n")     ##    try:         mxd=arcpy.mapping.MapDocument(MXDPath)         outPDFpath = mxd.filePath[:mxd.filePath.rfind('.')]+'.pdf'         ddp = mxd.dataDrivenPages         #MXD = MAP.MapDocument(mxd)             #arcpy.AddMessage(ddp, MXD)         if str(os.path.isfile(outPDFpath))=='False':                 outputPDF = arcpy.mapping.PDFDocumentCreate(outPDFpath)                 mxd.dataDrivenPages.exportToPDF(outPDFpath)             #arcpy.mapping.ExportToPDF(mxd, outPDFpath)                 name='c:\\temp\\'+outPDFpath[outPDFpath.rfind('\\'):]                 shutil.copy2(outPDFpath,name)                 arcpy.AddMessage('Printed '+str(outPDFpath)+' @ '+str(datetime.datetime.now())+"\n")         else:                 arcpy.AddMessage('File Exists '+str(outPDFpath)+' skipped pdf creation @ '+str(datetime.datetime.now())+"\n")     ##    except:     ##        errorm=arcpy.GetMessages()     ##        arcpy.AddMessage(' Error creating file...' +str(outPDFpath)+'\n'+str(errorm)+'\n'+' continuing...')              #MAP.PrintMap(MXD, printer)         #Remove variable reference to file     del mxd      **When run in debug** ![enter image description here](http://i.stack.imgur.com/mHbI5.png)


title: Is there a public Web Map Tiling Service (WMTS) available?

text: Does anyone know of an "almost always available" service URL for a public facing Web Map Tiling Service (WMTS)? I'm keen to test WMTS client capabilities without having to configure a WMTS server.


title: How do I intersect DEM and a polyline?

text: I'm working on Qgis, I need the longitudinal section of a river. Thanks


title: Clustering a number of points with the same coordinate around a central point

text: Is there a plugin in QGIS to randomly distribute points with the same coordinates around the coordinate? Similar to the "shift point plugin", but instead of producing a circle of points the points are distributed around the coordinate randomly within a set distance or polygon. Is there a plugin that does this or is there a combination of tools I could use to do this? I have looked online but can't seem to find anything.


title: Retrieve a point coordinates from map canvas using QGIS API

text: I'm trying to build a plugin that firstly captures the map coordinates. I have tried to adapt the map tool explained in this post: How to programatically check for a mouse click in QGIS The code I used for the main function is below:               def canvasReleaseEvent(self, event):         #Get the click         x = event.pos().x()         y = event.pos().y()              canvas = self.canvas         canvas.mapRenderer().setProjectionsEnabled(True)         canvas.mapRenderer().setDestinationCrs(QgsCoordinateReferenceSystem(27700))         espg = canvas.mapRenderer().destinationCrs().authid()              point = self.canvas.getCoordinateTransform().toMapCoordinates(x, y)      Unfortunately the when I print point it always returns as (0,0). I have tried various ways to set the crs of the points but none work. Any ideas anyone?


title: Problem with r.drain in Dufour and 2.1

text: I am trying to perform to define the shortest route between two points using r. walk, r.cost and r. drain. I have successfully done it using grass under Lisboa standallone living in win7. When I try to replicate my procedure under dufour or 2.1 (using the same data), I execute successfully r.walk and r. cost., nevertheless I get I big black raster as a result of r.drain. Data is the same, same DEM same friction cost and same points. The problem happens with dufour standallone 32, 64 and OsGeo 32 and 64 under win 7 64 and Vista 64. I would appreciate any hint on the matter. Thanks in advance Gerardo


title: How to open .bds file in qgis

text: I have a Land Registry boundary file which is in .bds (Base Dataset) format, a native file format in Cadcorp SIS. I want to open this file in QGIS. Is it possible, or do I have to do any conversions?


title: Has anyone else been unable to locate and install the GRASS plugin in QGIS 2.0.1?

text: After uploading and installing QGIS 2.0.1, I am unable to locate and install the GRASS plugin. The plugin is not listed in the 'Plugin Manager'. The 'Plugin_Installer' is also shown in red, indicating that there is a problem with the plugin. Has anyone else experienced the same problem? The look and feel of 2.0.1 is great though, and there are some nice new features.


title: Are there any web services providing country border data?

text: I am an Android developer (working in Java) and want to write a geographical quiz game. The quiz will need a list of countries and will randomly choose one to highlight on a map (probably using Google Maps). Then the user can type in the name of the country. What web services are available that will give me the polygon data of the border of a given country?


title: GRASS GIS - r.water.outlet - What is my mistake?

text: I'm new to GIS - more so in GRASS GIS. I installed GRASS GIS 6.4.3svn in my Windows 7 machine. No problem in installation.   * Downloaded SRTM-Geotiff file   * Location, Mapset all o.k   * Used _r.in.gdal_ to import SRTM file   * Filled the sinks using _r.terraflow_   * _Masked_ the filled raster   * Zoomed the required region & set the computational region   * Ran _r.watershed_ using _threshold=10000_   * Got drainge, basin, stream & visual rasters   * Visual rendering of rasters (drainge, basin, stream & visual) - o.k   * Different Watersheds (basins) displayed for entire area   * Zoomed stream raster to find lat/lon of mouth of river (cross checked with Google earth)   * Converted the lat/lon as decimal   * Ran _r.water.outlet_ using _stream_ as input & gave the lat/lon as easting & northing   * **Sadly the output display is blank**   * Rerun _r.water.outlet_ - only a small dot in the lat/lon point. No delineation of watershed for the given easting & northing   * All the above procedures as per manuals   * Self learning of GRASS GIS - no formal training **Since I'm new here, I'm not permitted to upload the result images - Sorry...** What is the mistake I'm doing? How to get the watershed (basin) for the given lat/lon? Kindly help. Thanks in advance...


title: Multiple polygons and rasters for Zonal Statistics as table in ModelBuilder

text: I have a large number (>17,000) of raster image files that I am creating by iteration. I also have clipped out corresponding polygon files with identical names in a separate folder (maybe not necessary, but my approach so far). I would like to calculate zonal statistics for these corresponding files. I have tried to set up a model that iterates through the rasters with a submodel that iterates through the polygon files with varying preconditions. When I run the model I am only getting statistics for the first file in the polygon folder. I think I understand why this is happening, but not an alternative way to set it up. I have included a picture of the main model and sub model. I am not very familiar with coding, so I would like to avoid that if possible, but am open to any suggestions! Let me know if there is any other information that would be helpful. Thanks! ![SubModel](http://i.stack.imgur.com/xlHVc.jpg) ![Main Model](http://i.stack.imgur.com/3tMdq.jpg)


title: Converting ms-access database into easy to use spatial database

text: I've recently come into ownership of an MS-Access database (2007 accdb), which I want to be able to use in a GIS. I know a little bit about GIS (I've used QGIS and ArcGIS quite a bit) but very little about databases. There is currently no spatial data in it but I do have various shapefiles which I want to link it to and ideally merge all these files/databases into one spatial database. If I had access to ArcMAP I'd probably just set up an ESRI Geodatabase, but that's not available. From what I've read PostGIS might be something I could use but as I say I know very little about databases beyond MS-Access. Is there a user friendly system someone could point me in the direction of? It would need to be compatible with MapInfo 7 (don't ask!) or QGIS. Any help/advice would be much appreciated.


title: how to use custom fonts in QGIS mapserver on Debian Linux

text: I've got a QGIS setup running, except for the custom fonts. They work just fine on my desktop machine (OS X), but they don't work on the server (Debian Linux). I've got the .ttf files installed in the system, so that should work. I also found something online about a fonts.txt file, so I created that:               Merriweather    /usr/share/fonts/truetype/Merriweather-Regular.ttf     Kingthings Calligraphica    /usr/share/fonts/truetype/Kingthings_Calligraphica_2-webfont.ttf      but it doesn't seem to do anything. Unfortunately, documentation on this seems to be lacking. Has anyone got this working and can tell me how to correctly setup my fonts?


title: How to increase symbol size of polygons when zoomed out?

text: I have many small polygon scattered. How do I increase size of polygons when zoomed out in qgis? As seen in the image the symbol appears to be lost when I zoom out. **Note:** I have created a thematic layer on polygon. **Background :** All layers are vector layer from database. ![enter image description here](http://i.stack.imgur.com/NWWDF.png)


title: Importing Bing Maps into application using ArcGIS Engine?

text: I'd like to create an application using ArcGIS Engine 10. I want to import Bing Maps into my application as base map. I surf the internet but couldn't find any tools or API to import these maps to my application. Is there any way to import the maps to application using ArcGIS Engine?


title: Unable to change the layer unit in QGIS

text: I have been working with QGIS for 2 weeks and I got the hang of it now. The problem now is I cannot calculate the area of a created polygon. Because the layer unit is not changing to meters. When changing from degrees to meters, it just jumps back. How to fix this problem? Thanks


title: Duplicate ObjectID in ArcSDE Table

text: After making a check in from an arcpad AXF file, i remark that there's duplicate object ID records , this Arcsde pgsql table is versioned and edited by Axf checkins and a .net engine application. the weird behaviour is that the objectid's from checkin start from the last objectID checked in from AXF in the last time, but between these times i added several records using the .net app so the objectids produced here are duplicated by the last check in, the .net app adds record directly through pgsql, so i set the objectID to (Max(ObjectID) + 1). i'm using ArcSDE 10.1 and PGsql 9.0.8. i'm not sure but i suppose that the problem come from generating ObjectID in .net app level, How i could generate ObjectID and inform Arcsde that this OID it has been used ?


title: Add 3rd party repositories button not in version 1.8.0 Lisboa

text: The available manuals and other help sources all refer to the 'Add 3rd Party Repositories button' in Plugins>Fetch Python Plugins and then the Repositories Tab. My problem is that the installation of QBIS 1.8.0 Lisboa does not have the Add 3rd Party Repositories button where stated. How do I then add 3rd party repositories please


title: Viewshed over Multiple DEM Models Possible?

text: I finally got the viewshed working correctly over one DEM model. However, I am wondering is it possible to do the viewshed over multiple DEMs? E.g. A Multiple County Viewshed based off one point. The dataset source does not allow downloads over 2gb. So I'm forced to download each county separately. I think a vaguely remember a way to combine rasters. If that is possible, which tool do I use ?


title: Making ArcGIS for Desktop play sound on speaker when model completes?

text: I would like to be able to sleep while Models are running for much of the night in ArcGIS for Desktop . I would like to have ArcGIS for Desktop play a sound when it finishes, to awaken me. Is this possible?


title: What is the max latitude and longitude values possible?

text: What is the max value of latitude and longitude values possible available that fit on with in the world coordinate system.


I hear people complain: > my biggest gripe with google spatial apis is that they hide the fact that > latitude/longitude is not a universal coord system My question is **Is latitude / longitude a universal coordinate system?** What are the implications if it is not?


title: Reducing the Size of a 3D model created in Sketchup

text: I created a building in sketchup. It is of higher size(25 mb), so I exported it as Collada file(6 mb). Here I reduced the size of the images to the maximum possible level. Still the collada file (3 mb) is of higher size. Is there any way to automatically reduce the size of collada file except image compression. The file is exported to kmz format, there the size is just 600 kb. But I want the file in collada format so that I can import it easily as an multipatch file in Arcgis. Any help towards 3D model Size reduction.


title: ST_AsLatLonText returns wrong data

text: I want to list all cities in `lubelskie vovoidenship`, which is relation 130919. I also want to get `latLong` data for each city. To perform this I do:               SELECT          name, ST_AsLatLonText(way), place      FROM          planet_osm_point     WHERE          place IN ('city', 'town', 'village', 'hamlet')         AND ST_Contains((SELECT way FROM planet_osm_polygon WHERE osm_id = -130919), way)      Query runs ok, but problem lays within the results. For example:               ...     "Lublin";"15°1'48.000"S 45°23'24.000"E";"city"     ...      The coordinates for **every** city are **wrong**. What should I do? **UPDATE** Data is imported using `osm2pgsql`:               osm2pgsql -U postgres -W -d gis2 --create ~/Pobrane/poland-latest.osm      Data downloaded from: http://download.geofabrik.de/europe/poland.html Result using `ST_AsEWKT()`:               "Lublin";"SRID=900913;POINT(2512485.39 6665744.97)";"city"


title: Query on the Google map engine

text: i have created a data store in the google maps engine and added to the google maps using V3.. which is using Visualization library, i need to Set query on it, for example Get only some features which are having name=chiru How can we do this? because google map engine gives only REST API..


title: How to make osm2pgsql import a tag?

text: I loaded OSM data into PostgeSQL with Osmosis&Osm2PgSQL, and see many tags are missing. I need to make it import `lanes` tag. How can I make it save this tag too? **edit:** my solution was to copy /usr/share/osm2pgsql/default.style into a new file, edit it adding a line:               way   lanes         int4         linear      then call osm2pgsql with parameter:               osm2pgsql ... --style my.style


title: Cannot display Oracle raster images using QGis

text: I am dealing with QGIS 2.2 and the plug-in Oracle spatial georaster. Apparently the connection to the database can be made, as the raster images are listed once you press 'Connect', but when try to display them, either I see a black box or anything at all. In theory my QGIs 2.2.0 has the GDAL version 1.10.1. I downloaded it using the OSGEO4w package for Windows. Have similar errors been reported? Is there any way of testing where the problem might be?


title: How to move a layer in QGIS Toc above an existing group?

text: If in QGIS TOC a group already exists and I add a new layer via               QgsVectorLayer* lyr= mQGisIface->addVectorLayer(uri->uri(),Title, myProviderName);      the new layer will be inserted within the group. How can I move it outside/above the group?


title: Raster Calculator Blurring Image Results

text: I am trying to run a NDVI analysis using raster calculator on Arg Map 10.2.1 and the results create blurry or bleeding image, almost as if the map has been dipped in water. I am using the default environment settings in the raster calculator. Both the bands being entered are very clear and have high resolution, I can see individual plants and rows in a farm field, but when I run the NDVI the resolution seems to get blurred. It almost seems like the raster calculator seems to automatically create a Guassian Blur for the image. Has anyone else encountered this and know of ways to overcome this blurring? Also I have also noticed that when running a raster calculator, the blurring only seems to occur when subtracting or dividing image bands, but not adding or multiplying, is there a reason for this? ![This is an example of the Bands resolution](http://i.stack.imgur.com/iLpoq.png) ![An example of the bad raster output](http://i.stack.imgur.com/T24vD.png) This is an example of the Bands resolution An example of the bad raster output


title: Replacing double quotes in filename using

text: I am working with data in an SDE, when I create a new feature class within the SDE I get a file name which contains my user name in double quotes. I do not have administrative permissions to the SDE and therefore cannot get the file name to write without having my user name in double quotes. I was wondering if there is a way to replace the double quotation marks in the file name so I can write an output to a shapefile. When I try to output to a shapefile I get an error stating 'invalid characters'. I've tried .replace('"','') but that doesn't seem to work. The only replace function I can get to work is replacing full stops in the file name by stating .replace(".","") using double quotes. Is there a way I can replace these troublesome double quotes. Thank you for for any light you can shed on this subject.


title: How can I set map extent in arcgis silverlight api?

text: Simple question but I don't know the answer. I want to set my map extent such that users see the map only in map extent and they don't allow to see greater extent. How can I set the extent map? Thanks in advance


title: Saving vector layer with specific symbology

text: I am creating a polygon feature layer with soil classes. I need to share this layer with some colleagues and I need the layer to maintain the color ramp I have created. Is there any specific format where the color ramp is saved together with the layer? Since I don´t know how to do that, I was thinking on saving the color ramp that I created as a style. How can I export this style? thank you very much for your help


title: Multipoint vs point PostGIS

text: Is it normal to use the geometry type 'multipoint' instead of 'point' knowing that I use it to store only one point. How can I convert from multipoint to point


title: Where can I find USGS SRTM 3 index grid or what procedure should I employ to select tiles?

text: I would like to download all USGS SRTM 3 tiles for Ontario, Canada from the USGS site: http://dds.cr.usgs.gov/srtm/version2_1/SRTM3/North_America/ Where can I obtain the index that will assist me in querying the tiles I need to download. Alternatively, is there another easier way to identify the tiles? When I was downloading the tiles for Manitoba, Canada I found a range of files then just downloaded the individual files from the site. Is there a similar range of files for Ontario? Thanks


title: Transferring from Windows to Linux OS for ESRI ArcReader

text: I have a .pmf file and I would like to open it under ubuntu. However all the layers are directing to a C:/ source. As Linux OS does not have such drive, how can I change the layer source path? Thanks! Edit: I only have ArcReader with me..... Is it a must that I need to use a editor to output a new .pmf again for linux OS?


title: How to use a SQL 2005 table as attribute data in QGIS

text: We have financial data in a SQL 2005 database which we would like to link to a cadastral shape file in QGIS. The idea is to be able to display up-to-date outstanding balances for each plot. How can we link QGIS to the SQL table. We are using QGIS 2.0.


title: Cannot upload an IMG file into PostGIS using raster2pgsql

text: I cannot upload into PostGIS a raster with the following code:               raster2pgsql -s 4326 -I -C -M "C:\lyr\gebco.img" -F public.gebco ¦ psql -d test -U postgres -W      It says ERROR: Unable to read raster file: public.gebco Is it because of the sentence, or because of the configuration at my database? I don't know what to change.


title: esri.dijit.Legend.startup not working in IE10

text: In Internet Explorer (IE) 10, legend.startup in the below code is not working. It just gets struck. The same is working in IE8.                   function generateLegend(results) {     legend = new esri.dijit.Legend({         map : map,         layerInfos : legendLayers     }, "legendDiv");     legend.startup();     }


title: Width and distance when using line patterns in QGIS

text: I am trying to make a cross hatch polygon fill (example link attached) using the values (defined by law) of 0.1mm for line width and 0.8mm for distance between the lines. But I don't get/see any hatch with such a small values, so I'd like to know if the values less than 1mm are supported anyway or I have to solve this differently? Thx very much on any information regards, dejan example link: http://i.stack.imgur.com/wG5Mo.png


title: maptool doesn't receive key pressed events after QMenu.popup() function

text: I am developing a qgis python plugin which uses a maptool receiving key pressed events by keyPressEvent(self, event) function, mouse movements events by canvasMoveEvent(self, event) function and mouse press events by canvasPressEvent(self, event) function. When the user right clicks on the canvas I want to display a contextual menu creating a QMenu object and running its popup function. When the control returns to canvas my maptool continue to receive only mouse movements and mouse press events and NOT key pressed events. Why canvas doesn't send key pressed events anymore? Is that any workaround? Thank you


title: Where can i find Shapefiles of city along with the world file (.jgw) and the jpg image?

text: Are there any links where i can download a shapefile of a certain city along with aerial view jpeg image and the corresponding world file (.jgw) ? I need it to test an algorithm. Thanks,


title: How do I embed a KML map using the "new" Google Maps?

text: I have a KML link that I can use to generate an embedded Google Map using "old" Google Maps and would like to use it instead to create a "new" Google Map embedding. To generate the an "old" embedding, all I need to do is click on the link icon on the page for my KML, which generates an `iframe` I can use for embedding:               <iframe width="425" height="350" frameborder="0" scrolling="no" marginheight="0" marginwidth="0" src="https://maps.google.com/?q=http:%2F%2Fshare.abvio.com%2F9d18%2Fde58%2F4d38%2Fa340%2FCyclemeter-Cycle-20140629-0905.kml&ie=UTF8&t=m&ll=42.429738,-71.31901&spn=0.135997,0.21338&output=embed"></iframe>      But this is an "old style" embedding, and I would like to generate a "new style" one from the same KML. Is there a simple mapping from one to the other or from my KML URL directly to the form of the "new `iframe`?


title: Registering SQL Spatial View in ArcSDE and updating Row_ID?

text: I am having issues registering an inner join SQL view with a spatial column in ArcSDE. Any help would be much appreciated. Here is the background info: ArcSDE 10, SQL Server Enterprise 2008 R2, SQL 'Geometry' Spatial Type, Non- Versioned database, DBO schema, My Windows login is a sysadmin user Here is the goal: The view I am trying to register is an inner join of a polygon feature class to a table. The left of the join is the table, which in the database has a ONE_TO_MANY relationship with the feature class, where the feature class is the primary key as ONE and the table is the foreign key as MANY. I need to display in ArcMap the polygon geometries of all the feature class records as joined to the table. Since this is a ONE_TO_MANY relationship, it will result in coincident/duplicate polygon geometries, but that's OK. Initially I am adding the view to SDE using sdetable -o create_view:               sdetable -o create_view -T Int_Veg_TEST -t "Interaction,Vegetation" -c     "Interaction.ObjectID AS INT_OID, Interaction.INT_UUID, Interaction.EVENT_UUID,     Vegetation.ObjectID AS VEG_OID, Vegetation.VEG_UUID, Vegetation.Shape" -w     "Interaction.EVENT_UUID = Vegetation.VEG_UUID" -i -D      This finishes OK, but the resulting view only displays the count and geometries from the right of the join (feature class) and I need it to display the feature class geometry and count as joined to the left (table). I would love to set the Row_ID when creating the view but that option is not available in -o create_view. When using the sdetable -o describe option, it shows that ArcSDE has pulled the feature class ObjectID field as the Row_ID:               Table Int_Veg_TEST:     Column name    Attribute    type    Null?  Length,DPs   RowID Column?     -------------------------------------------------------------------------------     INT_OID                  SE_INT32  NOT NULL   10     INT_UUID                 SE_UUID   NULL       38     EVENT_UUID               SE_UUID   NULL       38     VEG_OID                  SE_INT32  NOT NULL   10     VEG_UUID                 SE_UUID   NULL       38     SHAPE                    SE_SHAPE  NULL        0     OBJECTID                 SE_INT32  NOT NULL   10          SDE Set      From here I am having trouble using command line sdetable or sdelayer to update the Row_ID of an existing spatial view. Anyone know if this is possible using command line? I can update the Row_ID manually in the sde_table_registry rowid_column, and it works. All polygon geometries for the join are displaying in ArcMap. 2 problems here though: A: SDE still thinks it is maintaining the Row_ID. And I can't find where in the sde registry to manually update this to USER.               Table Int_Veg_TEST:     Column name   Attribute   type Null?   Length,DPs   RowID Column?          INT_OID       SE_INT32     NOT NULL        10          SDE Set     INT_UUID      SE_UUID      NULL            38     EVENT_UUID    SE_UUID      NULL            38     VEG_UUID      SE_UUID      NULL            38     VEG_OID       SE_INT32     NOT NULL        10     SHAPE         SE_SHAPE     NULL            0      B. I'd rather not update the SQL tables directly. I've also tried creating a new view using SQL Server Management Studio and registering it in SDE using sdelayer -o register, but this also does not work: SQL:               SET QUOTED_IDENTIFIER ON     GO           CREATE VIEW [Int_Veg_TEST]          WITH SCHEMABINDING          AS           SELECT dbo.INTERACTION.OBJECTID AS INT_OID, dbo.INTERACTION.INT_UUID,          dbo.INTERACTION.EVENT_UUID, dbo.VEGETATION.OBJECTID AS VEG_OID,      dbo.VEGETATION.VEG_UUID, dbo.VEGETATION.SHAPE          FROM dbo.INTERACTION     INNER JOIN dbo.VEGETATION ON dbo.INTERACTION.EVENT_UUID = dbo.VEGETATION.VEG_UUID          GO          CREATE UNIQUE CLUSTERED INDEX Ix_INT_OID_UNIQ_CLUSTERED     ON Int_Veg_TEST (INT_OID)          GO      command line:               sdelayer -o register -l Int_Veg_TEST,Shape -e a -t Geometry -C INT_OID,USER -i -D           Error: Underlying DBMS error (-51).     Error: Cannot Create Layer.     Microsoft SQL Server Native Client 10.0: Cannot alter 'INT_VEG_TEST' because it is not        a table.      More information is available.


title: leaflet.js ImageOverlay projection bug

text: We are trying to display a image over the leaflet map using "ImageOverlay" layer, but it does not work according to expected. In Google Maps works perfectly. Apparently Leaflet just stretches the image to fit the specified bounds while Google Maps reshapes the image. We are generating the image using matplotlib with longitude,latitude coordinates as x,y axis. **Leaflet Projection** ![Leaflet Projection](http://i.stack.imgur.com/ieBYJ.png) **Google Maps Projection** ![Google Maps Projection](http://i.stack.imgur.com/jM0aF.png)


title: Can conversion of MapInfo WOR to ArcGIS MXD preserve styles and coordinate systems?

text: Is there any ways to migrate from MapInfo's data with WOR project file to ArcGIS MXD project file without leaving all styles and coordinate system? I know that I can convert from TAB files to SHP files via ogr2ogr but what's about the styles information? P.S. I have MapInfo Professional 12 and ArcGIS Desktop 9.3.1. Also, I can't use the Tab Reader Tool now.


title: Why doesn't the mosaic to new raster tool or the mosaic tool work when I try to use them?

text: I am trying to mosaic 4 raster datasets of elevation together using the 'Mosaic to new raster' tool in ArcGIS 10.1. I have even tried using the 'mosaic' tool and it still doesn't work. The errors that appear state that the data sets must be in the same projection (Which they should already be, I checked in the properties), that they must have the same number of bands (which they do), or just a generic error message. Normally this tool works for me quite well. I'm not sure what the problem is this time, could it be that the rasters are in BAG format? I would appreciate any information that could help and I would be glad to offer any more necessary information in order to clarify my situation.


title: Help finding Postcode Boundaries for The UK, Spain, Canada, Ireland, Brazil?

text: I would be forever grateful for anyone's help finding current Postal Code shapefiles for any or all of these countries. Spain Brazil UK Ireland Canada Free is obviously better but because I may use them commercially I might need to pay for licensing.   Thanks for your time!


title: How to get the color of a vector layer in pyQGIS?

text: I have see this topic (How to change the color of a vector layer in pyQGIS?) but I'm looking for get the color not to change. How I can resolve my problem ? Thanks


title: Basic If/Then in ArcGIS 10.2 Field Calculator (Python or VB)

text: I've been looking through these threads and can't seem to figure out my answer yet. Very new to VBA and Python, but lots of years with Arc.I know I can do this the slow way with Select By Attributes, but It's time consuming. I am attempting a spatial join between cases of a disease (points layer) and US Census Tracts (polygon layer). This requires count data. For each point/case, I have a field called YEAR with date ranges 2001 to 2012 depending on what year the case happened. I need a count column for each year. For example, the first one I'm calling COUNT01. If the case date in YEAR is 2001, then COUNT01 will have a 1 in it. If it is any other year (2002-2012), then there will need to be a 0. I will have COUNT02, COUNT03...COUNT12 columns. Can't have a "Null" value. Here's what I've tried so far ![enter image description here](http://i.stack.imgur.com/O98HG.jpg) And here is what the data looks like in Arc. ![enter image description here](http://i.stack.imgur.com/pLkNo.jpg)


title: COGOing Legal Descriptions

text: I successfully used the COGO traverse tool to create a polygon for a piece of property described in a legal description. My issue arose because it did not line up with data I already have of the area. I did not expect it to line up perfectly, but I did expect it to be closer than it was. The projections for all of the data in the map are the same so I'm not sure if it is because there is a specific projection used for legal descriptions or what. If anyone has experience with COGO and legal descriptions I would appreciate the help. I am using ArcMap 10.2. The data I already have for the area is vast, but I am using a buildings feature class (polygon features), and a quarter section feature class (polygon features) from the datasets in our database. I also have a georeferenced image in my map that shows the property boundary surrounding the particular apartments I am looking at. This image is based on the legal description and comes from the engineers who worked on the apartment complex. Again, all data is set to the same projection. My steps include:   1. Creating a polyline shapfile and a polygon shapefile.   2. In my editor toolbar I went to editing options and changed the "Direction Type" to Quadrant Bearing, and the "Direction Units" to Degrees Minutes Seconds. I checked the "Ground to Grid Connection" box, then I made the "Distance Factor" 1 since my projection is in State Plane US Feet.   3. I started an editing session and began using the COGO traverse tool to create a polyline feature that according to the legal description commences at the corner of section 28 (remember I have quarter section map, which has the referenced section) and ends at the "True Point of Beginning."   4. From the True Point of Beginning I use the COGO traverse tool again to create a polygon for the property boundary based on the legal description. If more details about my steps are still needed please let me know, but for now I am trying to keep it simple to avoid boring or dissuading the reader from reading. Also, for information purposes, this project was given to me to see if it was possible to use GIS to map Legal Descriptions. So far I can't report yay or nay because yes I successfully mapped using a legal description, but no it did not end up where it should. It is very very close, but the boundary I COGO'ed does not cover all the apartments in my building feature class that it needs to. It just needs to move to the south about 100 ft. I am having a hard time believing the legal description is right and our building feature class is wrong because that would mean all the data we have in our database is wrong. I am also not convinced that the legal description is wrong. So that leaves me with thinking I made a technical user error. I can't make any recommendations until I figure out what the issue is and whether it can be fixed. This is a screenshot of what I have been working on. The polygon is supposed to surround the smaller buildings in the middle. ![map screenshot](http://i.stack.imgur.com/8P5Xr.png) This is a screenshot of the legal description I have been using. ![metes & bounds description](http://i.stack.imgur.com/3ERLP.png)


title: How can I style many layers of vector grid once in QGIS?

text: After I created many of vector grid layers, I now need to style them in a batch system. Is there any way to do this? If that is not possible, how can I use the layer duplicate option to inherit the style from the original grid; and then change the extents of the duplicated layer?


title: How to create a standard deviation map from three rasters in QGIS?

text: I have 3 normalised (0-100) raster maps. Is it possible to get a map which shows the standard deviation? I'd like to see the areas with smaller and higher std deviation. Thank you!


title: Lidar Elevations Data

text: I am trying to acquire a method of detecting sub cm elevation without manually measuring from the field. Can LIDAR data produce this for me? Example: A 1500 foot long pipe conveying water from one point to another. I need to know the rise and falls of the terrain where the pipe will be lying. Can LIDAR produce this for me? I am needing coverage for agricultural land throughout the Delta. Thanks.


title: Using MAX in SQL Expression of ArcGIS ModelBuilder?

text: I am trying to use a SQL expression in order to select values that are >= 25% of its max value in that column. For instance, `VALUE >= MAX(VALUE)*0.25`, and those will ultimately be the selected points it pulls out. Is there syntax to make that happen, or am I going to need to create a new column and use field calculator on it. Biggest issue with that is I will have to do this for about 250 files. I'm also trying to create this in modelbuilder. When I try it as above it just returns there is an error with the expression.