title: Spatial overlay data from many lines into a single line

text: Hello I have joined lines in my data set into longer lines (based on similar attributes) to reduce the number of rows in my data set. Now I want to take the information from the original smaller lines and update it into the larger lines (two different tables). I would like to concatenate the attributes coming from many lines into the columns for the larger single line. I'm using postGIS for this. thanks


title: Report Designer (ArcGIS 10.1) - Using Element Properties

text: I've finally got a table generated that will work for a street index and I'm attempting to use the ‘Create Report’ tool with ArcGIS 10.1 to generate what I need for a map. Essentially, I have it setup and working but the gap that is generated between each street name is to large. I've been playing around with the ‘Element Properties’ and realized that the 'Can Grow' behavior is where my problem may lie. If I set this to 'True' the street names longer then the space available will be given more space horizontally. A nice ability to have but by having this setup it seems to cause every other street name to take up more space horizontally. If I switch the ‘Can Grow’ to ‘False’ it fixes the issue but doesn’t allow longer street names enough space. This is a problem because I have limited space on my map. Has anyone had this issue before? Should I be using excel and if so is there documentation on the process? Am I missing something in the ‘Element Properties’? Appreciate any ideas. Thanks, Brian ![enter image description here](http://i.stack.imgur.com/7Jhto.jpg) ![enter image description here](http://i.stack.imgur.com/oSR4V.jpg) ![enter image description here](http://i.stack.imgur.com/fxPky.jpg) ![enter image description here](http://i.stack.imgur.com/Vjnbh.jpg)


title: How to pull lat/lng from a map canvas and update a database

text: So if i want to pull out lat/lngs after user clicks on a map canvas (lat/lngs are pulled using any javascript API (leaflet, wax etc.) and add this to a postgresql database using a query, how do i do this? I came across: https://github.com/commandprompt/postgres-js But how do you generally add values to a db on the fly?


title: How to insert data into fusion table using php?

text: I am using fusion tables to generate pins over the map. now I want to insert more coordinates data into the table using php but couldn't get any exact idea how to do it? Can anyone help me for that? Thanks in advance


title: Record Date and User for an edited feature

text: The client has a Feature Class in SDE 9.3.1. When a feature in the feature class is edited the client wants to automatically record the date/time and who edited it in the feature class attributes. I have found a few references to putting database triggers on the "A" table. Also for creating a class extension. Is there an plusses or minuses for doing either one? Are there any other options that I am missing? Thanks


title: All latitude and longitude within 1 mile circle area

text: I am trying to generate pins over the map of my database coordinates by using fusion table and I want to generate only within 1 mile area circle around the specific location. To get 1 mile area's coordinates similar to my db coordinates I've to put condition of < and > to fusion table but I don't know how to get all the coordinates within 1 mile area circle or what should be the value for condition in query. Can any one please help me to figure out, it'll really appreciable.


title: How to become a GIS Mobile Developer?

text: I am mid level GIS Analyst\GPS Coordinator who would like to become a GIS mobile developer at some point. At my workplace we work in the ESRI Platform,ArcGIS Server etc Being a GIS Analyst without any programming experience under my belt.I have started taking java programming classes.At some point in time i would like to create native apps. However,right now I am not sure how to get there. Should I start off by creating apps using ArcGIS server and then branch off in to native apps etc? So I guess my question is where do I start? Formal training, languages etc.There is so much information but i am having a hard time figuring out where to go. Hopefully some of you can shed some light from the prospects of someone who is experienced professional in the field. Thank you for your time and my apologies for the vagueness of my question.


title: Web Feature Service Question

text: If I set up a Web Feature Service that provides point (location) data, can my WFS also provide a bunch of additional data for events that occurred at those points? Is there a better way to do that? Would a GIS tool such as an ESRI product even know how to deal with the event data? Should the WFS only include the point data? What am I trying to do? I have a database full of data that has samples/observations that occurred at monitoring locations over time. Our database users want to be able to pull the data into their GIS tools over the web (not through database connectivity) and then do all sorts of cool stuff with the data in their GIS tools. Is WFS the rights type of service for what I'm trying to do?


title: How to make a private OSM (Save personal information to) a map using OpenLayer or Tilemill or mapbox

text: OSM (OpenStreetMap) is great for making a real time public map for general purposes but what if I want to add private info to a map such as the location and date of my first kiss or a note about every location I posted fliers etc. I'd like to know how to add private information to OSM. I have heard that Tilemill or OpenLayer are best for this but cannot find a video tutorial on adding private map data to a private layer and being able to back that info up and update it or transfer it to others to edit or update etc... As I research more I'm seeing that Tilemill/Mapbox are what I am looking for. Which makes no sense. It is expected that OSM has this functionality and what does the name TileMill have to do with OSM or maps at all? I think a lot is lost in new possible community members coming to OSM for this functionality and being totally confused and sent to something called TileMill. Still haven't found any vid tutorials for tilemill/mapbox or any clear info that these are infact what I am searching for.


title: drag and drop csv file and display map using leaflet js

text: I am new in php and leaflet js. I want to drag and drop .csv file in php and using the values of csv file, I want to display map using leaflet js. Can any one help me.


title: What steps does SAGA GIS take when calculating Topographic Wetness, specifically when handling depression filled (by ArcMap) DEM

text: I am a grad student involved with hydrological modeling. In the model I work on, we are calculating the **Topographic Index** ( **TI** ; without transmissivity - or set to 1) using Quinn's MFD algorithm. We do not use the exponent factor p, which is the same as having it set to 1. I am interested in determining how the TI calculation in SAGA handles large flat areas that might exist on a watershed DEM that has been delineated in ArcMap. The watershed delineation process in ArcMap involves filling depressions according to Jensen and Dominique, which creates these large flat areas. Our model uses a TI calculation algorithm that handles the flat areas reasonable well, but results in some striping in the direction of the processing (left -> right, then down, etc.). This results from the fact that we sort the elevation values from highest to lowest, moving from the top left corner to the bottom right during the sort, and then proceed with calculating the TI values from highest to lowest. When a large depression or flat area is reached, we temporally elevate the current cell to force flow into a neighbor or neighbors. This, due to the sorting having occurred left to right, then down, results in striping. The TI values are all close, and high (near the extreme of the distribution) as the should be, but not equal. But the SAGA TI algorithm seems to be able to give them a uniform value. I was wondering if anyone could explain what steps the SAGA TI algorithm takes when handling such large flat areas. See the two embedded pics below for more info - pay particular attention to the large flat area in the middle of the watershed. The SAGA output is on the top, my algorithm is on the bottom. I'm thinking that some sort of recursive tree search/dependence graph algorithm is needed. Any comments or ideas are greatly appreciated! I haven't been able to determine this from the help text. ![SAGA - TI](http://i.stack.imgur.com/lznNT.png) ![Model - TI](http://i.stack.imgur.com/qrCjI.png)


title: Leaflet marker cluster by Administrative boundary

text: Is it possible to perform clustering based on administrative boundary using the leaflet.markercluster plugin? How would I do that? I already have the boundary layer.


title: show .bil file with openlayers

text: There is a WMS service http://www.nasa.network.com/elev?REQUEST=GetCapabilities It can return files with the next formats: image/bil application/bil application/bil16 application/bil32 Then we want to draw them with OpenLayers. How can we do that? Is there any library (java preferable) to work with .bil format? Cheers, Dmitry


title: Convert several SpatialPointsDataFrames to raster

text: There are 100 SpatialPointsDataFrames in my Workspace. I loaded them on this way:               filenames <- list.files(path="",                             pattern="XYhectareTravelTimes_ez+.*shp")          for(i in filenames){                   filepath <- file.path("/",i)                   assign(i, readShapePoints(filepath))      They are called like this:               XYhectareTravelTimes_ez10.*shp      XYhectareTravelTimes_ez11.*shp      XYhectareTravelTimes_ez12.*shp       etc. How can I convert them to rasters looping through the Workspace? I am very new R user and hoping to find some help. thanks a lot.


title: Openlayers pop-up FeatureSelected doesn't invoke even FeatureClicked works well

text:
My code snippet is: //Vector Layer with listeners:               var vectorLayer = new OpenLayers.Layer.Vector("Kuzzular");              vectorLayer.events.on({             'featureselected': function(evt) {                 var feature = evt.feature;                 var popup = new OpenLayers.Popup.FramedCloud("popup",                         OpenLayers.LonLat.fromString(feature.geometry.toShortString()),                         null,                                 feature.attributes.name+"
"+feature.attributes.address+"
"+feature.attributes.speed + " km/h",                         null,                         true,                         null                 );                 popup.autoSize = true;                 popup.maxSize = new OpenLayers.Size(1600,1600);                 popup.fixedRelativePosition = true;                 popup.relativePosition = "tl";                 popup.panMapIfOutOfView = true;                 popup.calculateRelativePosition = function () {                     return 'tr';                 };                      feature.popup = popup;                 map.addPopup(popup);             },             'featureunselected': function(evt) {                 var feature = evt.feature;                 map.removePopup(feature.popup);                 feature.popup.destroy();                 feature.popup = null;             }         });                    map.addLayer(vectorLayer);     //create selector         var selector = new OpenLayers.Control.SelectFeature(vectorLayer,{             hover:true,             autoActivate:true         });              map.addControl(selector);      //create vector features               var feature =           new OpenLayers.Feature.Vector(                                 new OpenLayers.Geometry.Point( unit.lon,unit.lat ).transform(epsg4326, projectTo),                                 {name:unit.name, address:unit.address, speed:unit.speed},                                 {externalGraphic: '/assets/admin/layout/img/avatar/' + unit.id + '.png', graphicHeight: 74, graphicWidth: 74, graphicXOffset:-18, graphicYOffset:-71  }                         );                         feature.id = unit.id;                         vectorLayer.addFeatures(feature);      If I change `featureselected` with `featureclicked`, popups are shown well but they don't disappear as I wish. What is the problem here? Additionally, can I give an offset to the popup? And how?


title: Adding polygons from an existing layer to a new layer (shapefile)?

text: I have an existing shapefile that contains parcel data. I would like to add some of the existing parcels to a new shapefile. I am using ArcMap 10.1.


Is there any way to have the legend text above the symbol and beside it, not just above? What i do at present is to cut the end of the text off and insert a label then group them, however i want to resize the box and possibly have the text go into two lines, or preferably wrap around the symbol. Excuse the sloppy work below, a rush job for illustrative purposes. ![enter image description here](http://i.stack.imgur.com/8l3OM.jpg) ![enter image description here](http://i.stack.imgur.com/26MsX.jpg)


title: Add minutes contained in a column to a date

text: In CartoDB I want to add minutes contained in a column of data type number to a date. e.g. SELECT my_initial_date + 'xx minute'::INTERVAL as my_date FROM table What is working is:               SELECT my_initial_date + '7 minute'::INTERVAL as my_date FROM table      This adds 7 minutes to each my_initial_date. What I want to do is the same but with a variable number of minutes stored in another column, e.g.:               SELECT my_initial_date + 'column_name minute'::INTERVAL as my_date FROM table      This doesn't work in CartoDB! How to do this?


title: Should I use OpenLayers + GeoNetwork + GeoServer together?

text: At the moment I'm doing an internship and my assignment is as following:   * design (and build) a (prototype) web application    * with a map on which the users can show selected data.    * The server collects metadata and a direct download url from the source.    * Users can download the datasets directly from the source.   * A system like this that already exists and it should be able to use layers from this system through a WMS protocol. Now there has already been research done, but it wasn't done internally and I'm not sure if she understands what GIS is. She refers to GIS as a GIS server. GIS is the whole package though, isn't it? The research report came up with using OpenLayers, GeoServer and GeoNetwork together. OpenLayers is referred to as the webviewer, GeoServer as the GIS server and GeoNetwork as a geoportal. As I understand it though, GeoNetwork can display maps just fine. Should I disregard the research done and just use GeoServer + GeoNetwork? Or maybe some other combination? As long as it's free OpenSource (part of the assignment) the language doesn't matter. I've tried (and am still researching) to find the best solution myself, but a little nudge to the correct tools would be awesome. Before this I hadn't even heard of GIS, though I do find it pretty interesting.


title: ArcGIS JS API: Convert ESRI Date Type into String

text: I am following this ESRI tutorial to display information from a related table for the selected features. The example in itself for my project works fine, however, **my data include two Date fields** and when I display them on the right side grid, I see something like "1237507200000" as starting date and "1237161600000" as ending date. I understand I have to format this `esriFieldTypeDate` so I can display it as a text field like "YYYY-MM-DD", but I don't know how. It is **important to highlight** that this related records are stored in a `dojo.data.ItemFileReadStore` as in the example. Please, does anyone know how to do this? Thanks!


title: Is there a way to visualize a 3D full/solid body (like a full layer) in ArcScene

text: I am trying to present layer and simulated layer (how it decreased in a year, for example) in ArcScene, but showing two layers one on each other with a slight difference (like max 30 cm.) does not look really great.   I tried to show the difference, but taking the heights from properties it again makes it just a layer surface, not a full body. Is there a possibility to visualize a layer in ArcScene as a body for which you can choose the upper heights (like elevation) and the inner heights (like for ex. mean see level or so..)? _[edit]   _ So, the talk goes about Extruding, just stupidly I cannot find such option.. For illustrating. Two colorful layers (impossible to see two..) represent lets say 2010 and 2011 and the difference in between them I just made a separate layer (in greyscale). And the properties box in which there is no tab 'Extrude'.. ![enter image description here](http://i.stack.imgur.com/CHiug.jpg)


title: How to uncheck a layer in LayerSwitcher on map load?

text: I have four layers in the LayerSwitcher. When the map is loaded all the layers are checked by default. I do not want to show all Layers as checked.How do I uncheck a layer when the map is loaded?


title: Startup python script to re-join layers

text: In QGIS 2.2, there is a bug which involves using the **JOINS** properties on a vector layer (that is already being **joined** with several other vector layers) whereby the attributes are not fully loaded. Example:               Result_Layer > Middle_Layer > Layer_1                                 > Layer_2                                 > Layer_3      In the above case, the _Result_Layer_ will only show the values of its own columns and their respective fields and about a dozen columns coming from the _Middle_Layer_ and their respective values before showing **NULL** for the rest of the attribute table. The number of column and field attributes I am using in total is extensive ( ~200 columns; ~20,000 fields) so I haven't tested this on a much smaller scale. The workaround at the moment is to **re-join** the _Result_Layer_ to the _Middle_Layer_. * * * UPDATE: I've been looking at refreshing the `mapCanvas()` to see if that would somehow 'rejoin' the layers and bring back the values. I've tried the following in the Python console:               from qgis.utils import iface     iface.mapCanvas().refresh()      It probably refreshes but doesn't automatically rejoin. I tried also using `self.iface.mapCanvas().refresh()` but I get an error. This is for a client who wants to stick with QGIS 2.2 and does not want to use PostGIS/SpatiaLite as suggested by @mapBaker. Therefore, the only alternative I can think of is to write a script that automatically rejoins all layers when loading up the project. I would greatly appreciate any help/guidance with this problem (if it is possible!) Also, as a matter of curiosity, is this problem fixed in QGIS 2.4?


title: ESRI Map Cache Errors at Release 10.2.2

text: I recently posted this issue with esri. Has anyone experienced this?: We are experiencing brand new caching issues at 10.2.2. I am unable to either DELETE_TILES or RECREATE_ALL_TILES from an existing cache. The error we keep receiving is: Line 30Output failure, error string = Error moving bundle Failed to cache extent: -9223211.962860 3087885.815688 -9079003.170544 3271103.291998 at scale 577790.55428899999. The error occures when running RECREATE_ALL_TILES or DELETE_TILES operations of the Manage Map Cache Tiles (run as a tool or script) when run from a dedicated server gp cluster machine. Our clustered setup consists of a network file server storing our config, security and data-stores, as well as all server and system directories. Two servers serve as a mapCluster for map and image services, a single server in a gpCluster for caching and user-defined gp services, and a virtualized web server for the WebAdaptor. This has occurred as of lastweek whether reading from the GIS SDE database OR from a data-store file gdb. All permissions for directory shares and folder and file security have been accounted for. When the error occurs, the entire cache becomes corrupt, making any other caching operations in-operable (delete cache, delete tiles, etc). During the recreate, the bundlx file is being removed but the bundle file is not, thus not allowing the tool to overwrite thebundle from the admin-defined D:/arcgistemp folder(s) on the GP server. Thus far, the only solution has been to kill the service, stop the site, and then remove the entire cache directory, restart the site, recreate both the service and the cache. Has anyone encountered this?


title: SDE 10.1 Layers Not Loading in ArcGIS Engine Application

text: I have an ArcGIS application that works great when loading layers from SDE 10.0 but fails when trying to load the layers when they're pointing to a SDE 10.1 server. The only difference is the layer source in the MXD. When the layer points to the SDE 10.0 server, everything loads as expected. If I change the layer source to point to the exact same layer (name, rows, etc.; it was merely copied over from 10.0 to 10.1) in the SDE 10.1 server, a red exclamation point appears next to the layer name.   * The SDE 10.1 layers are visible in ArcMap so this shouldn't be a port or instance issue   * This isn't a licensing issue; I have the app set to exit if the license is invalid/unavailable Is there something special I need to do to get an Engine application to communicate with an SDE 10.1 server? Is this a service file issue? I tried adding a reference to the service name and port for the new server in the services file but it didn't seem to help, but 10.0 and 10.1 seem to do services/instances a little differently; perhaps I'm missing something. UPDATE: The more I think about this, the more it seems to be a "direct connect" issue. Is it possible that I'm not connecting to SDE because things aren't set up right on the client? If so, how do I set up direct connect on an ArcGIS Engine client?


title: Password protecting tables in cartodb

text: I am trying to find an easy way to have private tables, yet exposing them to certain clients if they are authenticated. This authentication can be achieved by a token system that allows the client to perform queries only if the token or the credentials are valid. The approach here was to have a function with the authentication logic and whatever queries or operations you need to expose, and call that function from the client along with his credentials. Problem with this is that you end up with a pg function for every operation, need to precisely define the input and output of the function, no polimorphy, so every time anything changes (like a data type) all functions need to be recreated. Not very flexible. So another approach would be to have just one function that acts as a wrapper for whatever SQL you send it, checks credentials, sanitizes the sql and sends back the response.               CREATE OR REPLACE FUNCTION runSQL(thesql text,client text, hash text) RETURNS SETOF RECORD as $$     DECLARE     sql text;     expiry time;     client_info RECORD;     m   varchar;     --here go the blacklisted strings     blacklist varchar[] := array['UPDATE ','INSERT ', 'ALTER ', 'DROP ', 'TRUNCATE', 'DELETE '];     BEGIN         expiry:=INTERVAL '1 hour';         sql:= 'SELECT * FROM auth WHERE client_id = $1 AND client_hash = $2 AND created_at > (NOW() - $3)';             EXECUTE sql using client, hash, expiry INTO client_info ;                 -- check if it's a valid user and has an up to date token             IF client_info IS NULL THEN                 RAISE EXCEPTION 'Authorization failed for client %', client;             END IF;              sql:= thesql;                 -- sanitize the sql         FOREACH m  IN ARRAY blacklist                 LOOP                 IF position(m in upper(thesql)) = 1 OR position(concat(' ', m) in upper(thesql))>0 THEN                         RAISE EXCEPTION 'Not enough privilleges';                 END IF;                 END LOOP;                 --send back results             RETURN QUERY EXECUTE sql ;     END;     $$ LANGUAGE plpgsql     SECURITY DEFINER;      And usage would be:               SELECT * FROM runSQL('SELECT cartodb_id, the_geom FROM sometable', 'client_id', 'token') as f(cartodb_id int, geom geometry)      My concern here is not opening a huge security hole in my database. Do you think it would be enough just to check for potentially dangerous strings in the query (like INSERT, DELETE, DROP, ALTER, GRANT, REVOKE...) or are there other security issues with this approach? Any other ideas on how to have a SECURITY DEFINER function that only allows certain operations like read, or access to just some tables? Thanks! David


title: Setting Script Tool Parameters Using Python and ArcGIS 9.3.1

text: Can't quite figure this out and i'm sure it's a simple fix. I have a script that utilizes the AGSSOM tool. It first stops all published services, allows for databases to be tidied up by Compressing/Uncompresing and Compacting, then resarts the services. The script is a tool in a tool box. I'm currently trying to set a parameter where the user picks which geodatabases to clean up (our services use several geodatabases and they are constantly being manipulated). Within the Parameters tab of the scripts properties, I've added a parameter and assinged it to Data Element as the Data Type (anyone have a better suggestion?). Additionally, I've set the MultiValue to Yes. Here is an example of some of the code.               inGDB = gp.GetParameterAsText(0)     gp.CompressFileGeodatabaseData_management(inGDB)      When I load multiple GDBs It crashes whith an error:                <class 'arcgisscripting.ExecuteError'>: Failed to execute. Parameters are not valid.     Input file geodatabase data: Dataset 'Z:\Data\Reference.gdb;Z:\Data\As-Built.gdb' does not exist     Failed to execute (CompressFileGeodatabaseData).          Failed to execute (ArcgisServerUpdater).       If I load a single GDB It works. Hope I'm not covering a similar question.


title: Detect Table of Contents Changes in ArcMap using ArcObjects

text: What is the preferred and/or most efficient method for detecting a change in the table of contents in ArcMap using ArcObjects? My current implementation relies on a dictionary with the names and indices of the layers.               public static Dictionary<string, int> layerEnumerator()     {         IEnumLayer enumLayer = ArcMap.Document.FocusMap.get_Layers(null, false);         enumLayer.Reset();         ILayer layer = enumLayer.Next();         Dictionary<string, int> layerDictionary = new Dictionary<string, int>();         int lyrNdx = 0;         while (!(layer == null))         {             if (!(layerDictionary.ContainsKey(layer.Name)))             {               layerDictionary.Add(layer.Name, lyrNdx);             }             lyrNdx++;             layer = enumLayer.Next();         }                  return layerDictionary;         }      Generally I will call this when I need to find the position of a specific layer in the ToC. Is there a more efficient means of gathering what's in the ToC at any given time (such as an event for ToC Contents Changed)? The solution I've provided above seems inadequate and inefficient at best.


title: QGIS Crashing in Ubuntu / 'Save As' browsing

text: When trying to create a new shapefile in QGIS on Ubuntu, I get a crash whenever I browse to a folder to save the file. Has anyone else experienced this?


title: How to mosaic Landsat 7 images correctly?

text: This is a result of mosaicing of two Landsat 7 images (layer stacked) in TIFF format. ![enter image description here](http://i.stack.imgur.com/ScycV.jpg) I was aming to put these images together with no overlaping visible line between them. The process was: Added raster layers (two tiff files) then Raster- Miscellaneous-Merge then input files and then named the output file in Tiff format. Where is my mistake? I use QGIS 1.8.0.


title: cad drawing overlay

text: I want to overlay a cad drawing on a google map using geoserver. Any tips on how i can do this. Do i need to convert to cad file to a shapefile. Im presuming i do. If so what tools are available to do this? Thanks in advance


title: Generating full planet tiles from OpenStreetMap data?

text: i am tested using two scenarios,      1) First, i am used switch2osm for tile rendering. its good it handles many things. But it takes long days to pre-render all tiles and it takes huge memory in TB. 2) Second, i am used tilemill to generate MBTiles and served using https://github.com/chelm/mbtiles-server But i dont know how much efficient for larger using. Other than this, Is any other way serving Tiles for whole planet in efficient way. Machine config : 1.5 TB and 128 GB main memory.


title: How to scale with radius attribute

text: new user of QGIS. I would like to style vector data based on an attribute (radius). The actual unit is in nautical miles. I don't want to use size-scale field or buffer to style the point rather I would like to make sure that the actual circle created would be drawn to scale based on its attribute. Thanks in advance.               LONGITUDE,LATITUDE,WIND,PRESSURE,CATEGORY,DATE,RADIUS(m)     145.9,6.5,55,985,4,11/5/2013 0:00,101860     144.6,6.5,60,980,4,11/5/2013 6:00,111120     142.9,6.9,65,975,5,11/5/2013 12:00,120380     141.3,7.1,75,965,5,11/5/2013 18:00,138900     139.7,7.3,85,950,5,11/6/2013 0:00,157420     138,7.6,95,930,5,11/6/2013 6:00,175940     136.2,7.9,105,920,5,11/6/2013 12:00,194460     134.4,8.2,115,905,5,11/6/2013 18:00,212980     132.8,8.7,115,905,5,11/7/2013 0:00,212980     131.1,9.3,115,905,5,11/7/2013 6:00,212980     129.1,10.2,125,895,5,11/7/2013 12:00,231500     126.9,10.6,125,895,5,11/7/2013 18:00,231500     124.8,11,110,910,5,11/8/2013 0:00,203720     122.5,11.4,90,940,5,11/8/2013 6:00,166680     120.5,11.9,90,940,5,11/8/2013 12:00,166680     118,12.2,90,940,5,11/8/2013 18:00,166680     116.6,12.3,90,940,5,11/9/2013 0:00,166680     114.8,13.5,90,940,5,11/9/2013 6:00,166680     113.1,14.4,85,945,5,11/9/2013 12:00,157420     111.4,15.4,80,950,5,11/9/2013 18:00,148160     110.3,16.5,75,955,5,11/10/2013 0:00,138900     109,17.9,70,960,5,11/10/2013 6:00,129640     108,19.4,65,965,5,11/10/2013 12:00,120380     107.5,20.4,60,975,4,11/10/2013 18:00,111120     107.1,21.5,40,990,3,11/11/2013 0:00,74080


I'm developing a map application (Several shapefiles with OpenstreetMap or BingMap on background and some dynamically generated features on the map). I'vs created some part of it with Telerik RadMap for silverlight. but I'm far more comfortable with Javascript and HTML (I'm mainly a web developer). So, is there any advises on how to choose the right technologies for a web based map application? is there Any shortcomings for Javascript components that support GeoJSON (leafletjs, KendoMap, ...)? Any comparisons between `GeoJSON and HTML/JavaScript vs Shapefile and Silverlight`?


title: Symbolizing File Geodatabase feature class with Layer File created from AutoCAD drawing?

text: I received an autocad drawing file which I converted into feature classes with the CAD to Geodatabase tool. **I would like to symbolize the exported data in ArcMap the same way it was symbolized in the input drawing file I received.** I visualized the drawing file in ArcMap and saved the "Drawing layers" Layer files as per usual but I am unable to import this symbology into the exported layer properties. When I import the saved layer files, the resulting "Import Symbology Matching Dialog" correctly matches all 3 fields (Linetype, Colo, LineWt) but clicking OK doesn't do anything - dialog disappears and symbology remains unchanged. Alternatively, when I load the layer file itself, it is linked with the DWG and source can only be changed for the CAD Feature Class type. How can i symbolize my exported data the same way as the original DWG is shown in ArcMap.


title: Accessing smb share on Trimble YUMA with python

text: I need to write a script that would import a .shp on the YUMA to a gdb on the script/host pc for manipulation with arcpy. I can't seem to find a smb python module. So Im not entirely sure on how to do this. Have I just not found it? surly its came up before.


title: Why is feature.geometry() None in my user-defined expression function?

text: I'm trying the user-defined expression functions for QGIS, as listed on Nathans blog: http://nathanw.net/2012/11/10/user-defined-expression-functions- for-qgis/ I always get the feature.geometry() as None, when running the user expression in the expression builder updating a text column? This happens for all geometries like point, lines and polygons. The data format is ESRI Shape. The code is in an userfunctions.py referenced from startup.py. Simplified example:               from qgis.utils import qgsfunction     from qgis.core import QGis          @qgsfunction(0, "Python")     def test(values, feature, parent):         geom = feature.geometry()           if geom is None:              return 'geom is None'              return 'geom is NOT None'      My text column is updated with the text: 'geom is None'. Why can I not reference the geometry of a feature?


title: How to configure a write only layer on Geoserver?

text: This is probably a wrong way of doing it, but I was exploring this option only because I do not know how to implement the right solution. We have a layer to which features are added using WFS-T. We have configured Geoserver to authenticate and authorize via LDAP. While querying for features, we would like Geoserver to return features based on the user/role. Since I do not know how to set feature based security (row level security), my thought was to see if we can make the layer write only and not allow any read operation. The read will be done through a SQL Parametric View layer will add a WHERE cause to filter by an unique value. For doing that, I have this setting in layers.properties workspace.layer.w=ROLE_USER workspace.layer.r=SUPERUSER However this doesn’t seem to work and I am not able to do any WFS-T on the layer although this user has the correct role. Reading s What would be the right strategy to implement this? Thanks in advance.


title: Can spatial data in Oracle be imported into SQL server?

text: Oracle supports spatial data and so does SQL Server. What are the things to consider if we have to migrate spatial data from Oracle to SQL Server? Are there any features in Oracle that are not there in SQL Server for supporting spatial data? We're talking about data maintained by ESRI ArcGIS if that affect answer of the question.


title: Cannot create new spatial database connection in ArcCatalog?

text: I have searched a lot in the net, but i just found Add Spatial Database Connection fails? here. However, I cannot see the confirmed answer. I want to know if anyone has found the solution? my ArcGIS desktop version and ArcSDE version is 10.0 and Sqlserver R2 2008. The problem is that when i click test in spatial database connection wizard in ArcCatalog tree can see the message "process succeed" but in fact after clicking ok nothing is created.


title: Do large vertex counts in PostGIS cause performance issues in WMS?

text: I am wondering if large geometry (geometry with large vertex count) in PostGIS might cause performance issues for GeoServer (or other product) in serving the WMS (not WFS)? This question also ask about effect of large vertex count but more in PostGIS process. Many thanks in advance.


title: Import arcpy using PyScripter or IDLE gives No module named RuntimeError?

text: I'm trying to figure out this problem. Yesterday I installed PyScripter and since then, scripts doesn't work. When I run any script (in PyScripter or IDLE) which uses "import arcpy", I get this error message:               import arcpy          Traceback (most recent call last):       File "<pyshell#0>", line 1, in <module>         import arcpy       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\__init__.py", line 17, in     <module>         from geoprocessing import gp       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\__init__.py", line 14, in <module>     from _base import *       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\_base.py", line 14, in <module>         import arcgisscripting       File "C:\Python26\ArcGIS10.0\lib\arcgisscripting.py", line 81, in <module>         from RuntimeError import RuntimeError     ImportError: No module named RuntimeError      Are there any suggestions on a potential to fix?


title: QGIS custom labeling by clicking on polyline (and polygon)

text: QGIS has now few really good options for custom labels but I can't find a way how easily manage labels along polylines. Is there any tool which helps me to place labels by clicking on polyline so the label will be anchored to that point on line - **to automatically set position above/on the line and proper rotation or even better curved with the line**? Similarly with polygon - click on polygon to exact spot where I want the label to be. I know there is workaround to create columns for labels x/y/rotation and move and rotate labels manually but that is often very time consuming and clunky way compared to some other GIS sw.


title: Geoserver WPS process list enabling /disabling through functionality code provided by geoserver WEBINF/lib

text: How can i enable or disable WPS processes through functionality code, i haven't found that related functionality code in geoserver WPS extension jar files those are in WEBINF/lib. Please guide me the specific location of WPS process functionality and enabling and disabling process


title: How I can sum 360 raster files in a folder with the Arcgis raster calculator?

text: In a folder I have 360 precipitation raster files and I want to sum them all to the raster calculator


title: Apache2 is missing conf.d file

text: I am following this tutorial about how to create my own tiles and got stuck on slide 20. I am trying to tell my web server about my mod_tile installation. However the folder `/etc/apache2` doesn't have a `conf.d`. How can I fix that?


title: Sensor Observation Service Filter Request

text: I am a student and currently developing a web application using Sensor Observation Service. In my Web-Client, I am sending request to SOS Server per sensor and getting back the observations. The main functionality of Web-Client is to perform Spatial Selection on map and see the Observations of the selected region. In this scenario, the region selected gives street/sensor Id's(which is fetched using WFS), and the client sends request to SOS for each street/sensor Id. So, it means if I select 1000 street Id's, I have to send 1000 request to SOS, which will obviously slow down my application. My question: Is there any way, I can send filter request to SOS ? Like, I send all my sensor Id's to SOS and get the observations of all at once? Or any other suggestion to improve my application behavior a part from SOS Filter Request? Thank you very much in Advance! Regards, Rajat Arora


title: FME error - failed to validate FME worker licence key

text: The FME licence I have access to is a floating one. The issue seems to be easy to recognize (see attached log) however there are some doubts I would like to ask for. While the process is really machine consuming I have some suspicious that the error of the FME is not due to "Failed to validate FME worker licence key" but it's more related to the task the machine has to deal with. Can it be ? The error always occurs in the same point of proccessing time. How does FME check for the licence key ? Just at the start of the workbench or there are some "milestones" during the processes where the licence has to be confirmed again ? Please find below a part of the log. Thanks in advance for the help. Gregory * * * _Failed to validate FME worker license key. Program Terminating Translation FAILED. Failed to validate FME worker license key. Program Terminating Translation FAILED. Failed to validate FME worker license key. Program Terminating Translation FAILED. ... Last line repeated 2 times ... FME Session Duration: 4 hours 56 minutes 38.5 seconds. (CPU: 220.4s user, 11.6s system) END - ProcessID: 1324, peak process memory usage: 1074456 kB, current process memory usage: 355076 kB Error in child process '33', exited with code '2' Program Terminating Translation FAILED. Warning: not all FMESessions that were created were destroyed before shutdown. This may cause instability_


title: Exporting data from a WFS service to csv file using qgis

text: I have a service that contains a list of points with some meta data and I would like to export this to a csv file. I remember doing this some months ago and it worked well, but now I can't find the method (or plugin?) that I used. Has anyone done something similar and know how to export data from a WFS service to a CSV file using QGIS?


title: Why does a tool not call its constructor or OnActivate() method?

text: I'm working on an add-in for ArcMap 10.0 in .Net version 3.5. It was finely added in arcgis desktop 10.0 as addin tool and working finely. I added some code in the addin project and build that project again. But At this point, when I run ArcMap and select the tool, the tool is deactivated. In fact, when I put break point on my project and debug it, tool's any method are not called. I wonder, the OnActivate() method is not called either. Even the constructor of the Addin Tool class is not called. ### Question: How do I activate my tool now? Why does the tool not call it's constructor or OnActivate() method?


title: Issues with installing pgrouting in mac

text: I installed both postgresql and postgis using brew in mac. How can I now install pgrouting. I couldn't install it using brew. Any suggestions?


title: What is difference between map layer and spatial dataset (e.g. feature class, shapefile)?

text: This may be very trivial question. To me they seems very similar in concept. What is the difference between feature class and feature layer?


title: How much is 1 meter in decimal degrees?

text: I supposed to know this but unfortunately I do not know. I am using Tester to select objects that are bigger than 1 m. My data are in EPSG:4326. How much is 1 meter in decimal degrees?


title: I'm looking for an idea to calculate shortest path in a line Vector layer using a starting point

text: I'm looking for an idea to calculate shortest path time in a line Vector layer using a starting point in a qgis python plugin based on distance and speed. is there a solution where i can get ideas? I'm thinking to use Dijkstra's algorithm for calculation and i found a link with portion of a python code, but i need to know how do i use parameters as distance and speed, and manage intersections.


title: QGIS 2.2 freezes when I select a Grass command via Processing toolbox

text: I'm on Ubuntu 13.10, and I have never been able to get the Grass plugin to work in QGIS 2.2 via the Processing toolbox. As soon as I select any Grass command, QGIS locks up, no menu appears. I then have to close QGIS. I've tried Grass 6.4.3-2 and 6.4.3-3 from the stable and unstable repos, on a fresh install of Ubuntu on my Desktop and Laptop (6.4.3-3 is supposed to fix a possibly identical problem, but the "fix" doesn’t work for me). Is anyone else having trouble getting Grass to work under Ubuntu 13.10 via the Processing toolbox? Any solutions? Thanks.


title: Is it possible to connect Maptitude 5 with Oracle Spatial?

text: I have a version of Maptitude 5, and I need to connect it to Oracle Spatial. I have the client installed on the machine, including the connection to the database using Qgis and gvSIG happens normally. Using Maptitude: 1-Open File Menu 2-Files of type: Oracle Spatial Layer 3-In the window that opened've User, Password and Service 4-E gave an error that the connection was not performed with these parameters, but these are the parameters used in Qgis and gvSIG. I know this is an old version of Maptitude, but I wonder if anyone uses or has used this version connected to Oracle Spatial.


title: How to split apart overlapping features

text: I have a set of 20 features that are the exact same shape (a rectangle) that are stacked on top of one another in qGIS. Is there a way to programmatically separate them? Ideally I'd like to separate them into smaller chunks in the same space. I.e, if originally they are all 10 miles across and stacked on top of each other, have it separate them so they are each .5 miles across and staggered so as one ends, the next one begins. Is that possible or do i need to do it by hand?


title: Stratified Random sampling in ArcGIS

text: I am trying to generate random points for a survey, but I want the points to fall based on 2 polygons. I have created a polygon fishnet, and have a layer of forested habitats. I am trying to generate points so that there is one point per grid, and so that the point falls within a forested habitat in the grid. Is there a straightforward way to do this?


title: How do you calculate the area of polygons which overlap a buffer around a point?

text: I'm new to QGIS (Dufour 2.0.1) and have plotted nests as points on a googlemaps base-map and have created buffer zones around the nests. I've also plotted ponds as polygons on the same base-map. I want to know how to calculate the area of the ponds which fall within the buffer zone PER NEST.


title: Globe draped layer updating - ArcObjects vs Desktop

text: I have a globe document that includes several layers that are draped on the globe's surface, but no elevation layers. If I open the document in ArcGlobe (desktop) and add an elevation layer, the appearanceof the draped layers changes - they are draped on the elevation. However, if I use this same globe document in a Visual Studio application and add the elevation layer programmatically, the appearance of the draped layers doesn't change. Does anyone know why this is the case? Here's the code I'm using to add the elevation layer programmatically: IGlobe theGlobe = axGlobeControl1.Globe; IRasterLayer mergedElevLyr = new RasterLayerClass(); mergedElevLyr.CreateFromFilePath("D:\@ALPACAS\AlpacasCT\experimental\3D\data\MergedElev_111711.img"); mergedElevLyr.Visible = false; theGlobe.AddLayerType((ILayer)mergedElevLyr, esriGlobeLayerType.esriGlobeLayerTypeElevation, true); theGlobe.GlobeDisplay.RefreshViewers(); Any thoughts will be appreciated! Mark


title: VBA ArcGIS 9.3 Field Calculator: 'countif' across several fields

text: I have an attribute table with 15 columns that contain either a '0' or a taxonomy number. I would like a summarized column that shows the number of different species.               location  type1  type2  type3  etc etc...  nr_of_species         1      123    234    0      .........      2         2      456     0     0      .........      1         3       0      0     0      .........      0      Is there a way to do this with VBA and field calculator? I have no experience with VBA, so i might need a complete dummy-guide. thanks!


title: Setting transparency in background colors of neatline fame?

text: Is there an easy way to set transparency levels in background colors of a neatline frame or any other frame for that matter in ArcGIS?


title: How can I disable and enable an ArcGIS 10 Add-In button?

text: Is there a way to enable / disable Add-In buttons / tools?


title: Set visible/ratio for make feature layer based on a list

text: So, I have a feature class with three integer fields. I need to create a feature layer based on this feature class to use it later for running some GP tools where the values of the integer fields will be used for some calculations (with ratio).               fieldA fieldB fieldC       25     60     70      However, I don't know what field should be used, this is something that I read from a look-up table. So, when I have read the look-up table, I get a list that looks like:                [fieldA,fieldB]      Great, now I know which fields I want to use ratio for, i.e., which one will be taken into analysis. However, I am not sure how I should approach this, that is, how can I specify in the fieldinfo variable what field should have VISIBLE RATIO based on the input list with their names. The whole fieldinfo variable is just a string and I can't see any smart way of "getting VISIBLE RATIO right after the fieldA and the fieldB" (as in the example scenario).               fieldinfo = "fieldA fieldA VISIBLE RATIO; \                  fieldB fieldB VISIBLE RATIO; \                  fieldC fieldC HIDDEN NONE"          arcpy.MakeFeatureLayer_management(input_fc,inputfeaturelayer,"#","#",fieldinfo)      I looked at the `fieldinfo` arcpy class, but not really sure if this will help me here and if yes how I should use it. **EDIT:** before using the `field_info` I was forced to use the str.replace() to update the `fieldinfo` object.               for field in fields_list:              updated_makefeaturelayer_fieldmapping = updated_makefeaturelayer_fieldmapping.replace('{0} {0} HIDDEN NONE'.format(field),'{0} {0} VISIBLE RATIO'.format(field))


title: How can I get responsive pinch-zoom when using OpenLayers 3 and Google Maps

text: Following the OpenLayers 3 Google Maps integration example, I'm able to put my own custom tile layer over google maps imagery. However, on mobile devices and tablets, pinching-to-zoom only resizes the custom tile layer during the duration of the pinching. The google map doesn't zoom until the user lets go of the pinch. The above OL3 example behaves the same way on those devices; the custom vector layer is the only thing that resizes on pinch until the user lets go. I'm assuming this is due to how the google maps object is so decoupled from OpenLayers in version 3. Per the example's code, the apparent prescribed way of getting the google map to zoom is to set up an event handler for change:resolution in the OpenLayers view which will in turn set a zoom level on the google maps object:               view.on('change:resolution', function() {       gmap.setZoom(view.getZoom());     });      I'm assuming that this event is only firing when the user lets go of the pinch, which would explain this behavior. Aside from switching to a proper OL3-managed Bing Maps layer, are there any good workarounds for this?


title: OpenLayers - possible to display detailed GIS data?

text: I have an OpenLayers app to plot point features, but need to read in a very detailed multipolygon vector layer to plot an accurate location. I'm currently loading the polygon layer via TinyOWS. I've set maxZoom to 2000 in OpenLayers and am using BBOX strategy. When I zoom in the first load is ok, but the whole layer gets loaded and I can see it when I zoom out despite the maxZoom. The application now runs too slowly. I'm new to this and wondered if someone could point me in the right direction. I got som pointers here: Different level of detail(layers) on different zoom level at OpenLayers map and here: OpenLayers GML Polygon Performance Here's the WFS call for the polygon layer:               pos_skotsel = new OpenLayers.Layer.Vector("Pos_Skötsel", {         strategies : [new OpenLayers.Strategy.BBOX()],         projection: new OpenLayers.Projection("EPSG:4326"),         protocol: new OpenLayers.Protocol.WFS({             version: "1.1.0",             srsName: "EPSG:4326",             url: "http://XX.XXX.XXX.XXX/cgi-bin/tinyows",             featurePrefix: "tows",             featureNS :  "http://www.tinyows.org/",             featureType: "XXXX",             maxScale: 2500,             geometryName: "geom"         })     });           map.addLayers([gphy, wfs, pos_skotsel]);      Thanks in advance, Mini


title: AGS 10.1 geoprocessing service returning the wrong URL

text: I have an image export geoprocessing service that runs on ArcGIS server 10.1. It works by writing an image out to the scratch folder, and returning the location of that image. Basically just this:               Output_File = os.path.join(arcpy.env.scratchFolder, output)     #Export the file ...     arcpy.SetParameterAsText(8, Output_File)      This works fine in most situations, and returns a url that looks like this:               http://www.oursite.com/arcgis/rest/directories/arcgisjobs/...      but on one client's network, when they connect to our hosted server, the service returns something that looks a lot more like this, despite the fact that the same instance of the service returns normally otherwise:               http://(local server name).(internal domain name).local:6080/arcgis/rest/directories/arcgisjobs/...      Is this the result of a misconfiguration on the server? And is there a reason why it would return the local network name of the server in some contexts?


title: Default value in WPS input

text: I'm using GeoServer and want to add a WPS with default values. How can I define an optional parameter as a WPS input? I think, that I have to define this parameter with a multiplicity of 0, but how can I use a default value? Do I have to add something in the ApplicationContext.xml? Thanks for helping.


title: Label all features of shapefile with PyQGIS

text: I have a problem when trying to label all features of a shapefile with a Python routine. If I open my shapefile in QGIS "by hand", and select a .qml with pre-defined label selection, the labels are correctly displayed. If I load the .qml file via a Python routine (see line below), all style features (colors etc.) are there, EXCEPT for the labels. No idea why. `myLayer.loadNamedStyle("someStyle.qml")` Since the line above did not give me the desired labels, I tried to create a more extensive routine, see below.               myProvider = myLayer.dataProvider()     myLayer.startEditing()     for myFeat in myLayer.getFeatures():        #Some calculations in here...        #The variable "Some_Value" is what I'd like to use as label             Index_myLabelField=myLayer.fieldNameIndex("myLabelField")        Value_myField={myFeat.id():{Index_myLabelField:QVariant(Some_Value)}}        myProvider.changeAttributeValues(Value_myField)        myLayer.commitChanges()             myLayer.label().setLabelField(0,Index_myLabelField)        myLayer.label().labelField(0)        print myLayer.label().fieldValue(0,myFeat)   #Returns the desired label string             #Create map renderer        mapRenderer = QgsMapRenderer()        rect = myLayer.extent()        mapRenderer.setExtent(rect)        mapRenderer.setLayerSet(QStringList(myLayer.id()))        myRenderContext = mapRenderer.rendererContext()             myLayer.select(myFeat.id())        print myLayer.isValid()  #Returns true        myLayer.label().renderLabel(myRenderContext,myFeat,1)   #Error: Segmentation fault (core dumped)      Until the last line, where I try to render the label, everything seems to work fine. So what is going wrong in the last line? There are two things I am unsure about:   1. In the QGIS API, the function to tell QGIS the field I want to use as label is  setLabelField (int attr, int fieldIndex) What is the "attr" in this context? What should I put here? I used "0" so far...   2. In the renderLabel() function, there is a "bool selected" needed in third position of the call. So what I am doing is selecting the feature (myFeat), and telling the function a "1" here. Is that correct? This post Rendering and Labeling Shapefile with PyQGIS does not fully help me, as the OP there does not seem to label all features of his shapefile. I am able to add one label to the shapefile (kind of like a title) with the described method there, but I am looking to have a label next to every feature (e.g. like street names on a map). I would appreciate any suggestions here. If you do know, why the qml-loading in the beginning of my post does not work, please also let me know. I really don't care HOW I get my labels in there, I just need them somehow. Thanks for any help!


title: Digitizing with OpenLayers Plugin seems incorrect

text: I am attempting to digitize some buildings based on floor plans using QGIS. To do this, I have taken the blueprints, georeferenced them against openlayers bing aerial plugin. I then digitized the floor plans in EPSG: 3857. That all seems to work (Also, much faster to "save as image" and add the bing image back to the map). My issue... I added the digitized shapefile to ArcMap with the ArgGIS imagery basemap, and it's not anywhere near correct... I believe the projection used for that is web mercator, which seems to be the same as the Pseudo Mercator that I've used (I projected the file in ArcToolbox from Psuedo Merc to Web Merc and they lined up perfectly). Does this mean that the imagery from the OpenLayers plugin is off? I know, I guess since I have ArcGIS access, I could just do it there, but I really want to be able to use Quantum in order to scale this project up a bit without costly ArcGIS licenses. Thanks for the help! Here are pictures that show what I'm talking about: QGIS with Bing aerial with labels in openlayers. Upgraded to 2.3 after comment below, but the issue still exists: ![QGIS with Bing aerial with labels in openlayers. Upgraded to 2.3 after comment below, but the issue still exists:](http://i.stack.imgur.com/SDkLd.jpg) ArcMap with Imagery with labels turned on... shouldn't be in the desert, which is in the correct state. And if I project the shape to the 102100 that is web mercator, there is no change... ![ArcMap with Imagery with labels turned on... shouldn't be in the desert, which is in the correct state. And if I project the shape to the 102100 that is web mercator, there is no change...](http://i.stack.imgur.com/wZZGz.jpg)


title: Get Coordinates of selected item on map using ArcGIS Engine

text: I have `IFeatureSelection`. Then select from it items by `QueryFilter`. How can i get coordinates of selected item? (If there are more then 1 object, then select first)? Selected item is Point in Attribute Table. I am using C# to write code. * * * Using samples and Links from @artwork21 and @GISKid i think that i can get coordinates of selected item using this code:               var fClass = featureLayer.FeatureClass;     var fCursor = fClass.Update(null, false);     var aFeature = fCursor.NextFeature();     var tPoint = aFeature.Shape as IPoint;     // then i get Coordinates from tPoint.      But i get new problem - method `Update` isn't implemented. Have SomeBody any advices? Thanks.


title: Using SWAT model with QGIS

text: I am trying to use SWAT (Soil and Water Assessment Tool) for watershed management activities. It is a set of commandline tools and needs a complex set of input data some of which can be generated from DEMs. There is a plugin available for ArcGIS and MapWindow which makes life easier. I want to know if there is any project/work going on to integrate SWAT with QGIS. Or, is there any other soil and watershed modelling tool available for QGIS which can generate/estimate runoffs, sediment load, erosion vulnerability etc. from relevant input data like rainfall, slopes, stream network, soil type, land use etc.


title: Finding Polygons which do not touch any other polygon in a single oracle database table

text: I currently have an oracle table of approximately 100,000 records - all of which are polygons. I want to find all the polygons which do not touch any other polygon within the same table. There are a few previous posts which use the ST_TOUCHES query but this relies on you having two geometries to compare against from different tables. Unless my understanding is incorrect? Where as what I want is a query which looks at the one table and one Geometry column and looks to see if a polygon touches another.


title: Error 999998 in Arcpy 10.1

text: I found an RuntimeError: ERROR 999998: Unexpected Error when lunching Arcpy in my ArcGis 10.1 such as this:               Runtime error      Traceback (most recent call last):       File "", line 50, in        File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\sa\Functions.py", line 3840, in Minus         in_raster_or_constant2)       File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\sa\Utils.py", line 47, in swapper         result = wrapper(*args, **kwargs)       File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\sa\Functions.py", line 3837, in wrapper         return _wrapLocalFunctionRaster(u"Minus_sa", ["Minus", in_raster_or_constant1, in_raster_or_constant2])     RuntimeError: ERROR 999998: Unexpected Error.      I've extensively gone through forums and found similar problems, but nothing proposed has resolved this problem. I haven't tried to reinstall ArcGIS, though. My code is: Also, when I execute the code, it gives OK until the `Print "Ok1"`, which means that `transmissionCapacity = Minus(tc1, losses)` is not executed. Error 999998 worries me a lot. Any wise advice from you guys? Thanks.


title: Countries not loaded in Nominatim base installation?

text: I just installed the Nominatim geocoder and loaded the Andorra.pbf dataset. I saw that in the Nominatim project, there are multiple .sql file in the data/ folder with country names and polygon that were injected in the database during the installation process. Files: - country_name.sql - country_naturalearthdata.sql - country_osm_grid.sql - us_state.sql - us_country.sql I was assuming that with this base installation I would be able to query Nominatim with different country names and that it would provide me results... but it does not. No result found each time. Did I miss a step in the installation process to be able to use the content of those file in the Nominatim query? Is there a way to use this data? Thanks


title: Georeferencing Tools in ArcGIS Engine 10.1

text: I'm new with ArcGIS or even GIS it self. My question might be a newbie question but I need help on building a Georeferencing tool set like the one in available ArcMap. I found that it is not available to be added to toolbar control. BTW I'm build the program using Windows Form and C#.


title: Implement a least cost path in a graph

text: I am trying to produce a routing functionality in a web-GIS that, given a start and an end point, will retrieve the least cost path taking into account points along edges representing major costs, thus avoiding them. I would like to ask what I could do to start producing this function. I have a graph with edges and nodes, I have points representing costs along edges, I have quite a good knowledge of Python and I was thinking to use PostGIS to store my data. I had a look at the pgRouting extension of PostGIS and it seems to be the solution I am looking for. Anyway I would like to ask someone more experienced to provide advice on what I need, what I need to know and how to get started. Thanks!


title: Add Oracle spatial view Layer

text: Is it possible to add an Oracle Spatial View on QGIS MAP? My view is registered on MDSYS metadata but, when I connect with Oracle DB, it is dissabled on the list. QGIS Desktop 2.2.0 Oracle Database 11g Enterprise Edition Release 11.2.0.4.0 - 64bit Production PL/SQL Release 11.2.0.4.0 - Production CORE 11.2.0.4.0 Production TNS for 64-bit Windows: Version 11.2.0.4.0 - Production NLSRTL Version 11.2.0.4.0 - Production tks!


title: storing date range in PostgreSQL/PostGIS for time parameter in WMS/WFS

text: Using GeoServer I want to deliver some objects/features with WFS and WMS. The objects and their attributes have a date range during which they are valid. Say, we have street lamps and each street lamp has an `ID`, `POINT (geom)`, `start_date`, `end_date`. The WMS and WFS request can use the `time` parameter to filter the resulting objects (see WFS GetFeature with time, WMS GetMap with time). Using the time filter, only the street lamps that a within `start_date` and `end_date` would be returned by GeoServer . However, I now saw that for `time` the respective time column must be stored with a `date` type. That is, in geoserver one cannot configure start or end date columns for a WMS or WFS layer. It is only possible to configure a single date column (same for MapServer). What's the best way to store my street lamps in PostgreSQL/PostGIS? For a lamp that was created on Jan 01, 2014 and destroyed on Jan 20, 2014 would I need to have 20 data sets in my database:               ID GEOM DATE     01 xy   2014-01-01     01 xy   2014-01-02     ...     01 xy   2014-01-19     01 xy   2014-01-20      Is there a way to store date ranges in PostgreSQL/PostGIS that can be used by GeoServer? Or is this impossible?


title: How to save and load vector layer attribute value map, editable state in .qml style file?

text: I just upgraded to Qgis 2.4 and there were problems with style file:   * When I load style file for Qgis <= 2.2, the color of features is right but when I edit them, the value map (value - description pair) is lost.   * The value map I created for layers in Qgis 2.4 is not kept in .qml style file when I choose 'save style' as 'Qgis layer style file'. Is this a bug in Qgis 2.4? **ADDED INFORMATION:** Sample files: http://www.mediafire.com/download/bbk02sj3jl9iwtc/test_data.7z For line layer _line.shp_ with 3 attributes: id, name, type - I have style file _2_2_line_style.qml_ generated from Qgis 2.2 which have following lines:                 <edittypes>         <edittype labelontop="0" editable="0" type="0" name="id"/>         <edittype labelontop="0" editable="1" type="0" name="name"/>         <edittype labelontop="0" editable="1" type="3" name="type">           <valuepair key="first type" value="1"/>           <valuepair key="second type" value="2"/>           <valuepair key="third type" value="3"/>         </edittype>       </edittypes>      In Qgis 2.2, when I load that line layer with that style and add more lines, there will be a pop up window that allows me to enter values with a drop down list: ![enter image description here](http://i.stack.imgur.com/0zwJ9.png) But in Qgis 2.4, when I load the same layer with the style file and add (create new) features, **the drop down list is gone** : ![enter image description here](http://i.stack.imgur.com/nlXZF.png) More over, in Qgis 2.4, when I set the value map for that line layer manually, I can get the drop down list but when I save the style as .qml and inspect the file, there isn't any key - value line like the 2.2 style file. ![enter image description here](http://i.stack.imgur.com/NsjDk.png)


title: json issue: GetFeatureInfo with raster images

text: I have succesfully typed a GetFeatureInfo request with the following structure directly into my browser:               http://192.168.1.20/geoserver/workspace/ows?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetFeatureInfo&CRS=CRS:84&BBOX=39.32295698174142,-8.4236378760969881,39.363077494358599,-8.3933494101086144&WIDTH=1114&HEIGHT=841&LAYERS=layer&QUERY_LAYERS=layer&FORMAT=image/png&INFO_FORMAT=application/json&FEATURE_COUNT=50&I=396&J=579         It returns me onto the screen:               {"type":"FeatureCollection","features":[{"type":"Feature","id":"","geometry":null,"properties":{"GRAY_INDEX":162}}]}          This image was a 8-bit jpeg. Now, when I try to do the same over a 16-bit raster image, it does not return any value:               {"type":"FeatureCollection","features":[]}          I was wondering whether the raster type is a key-fact0r on this issue. Do I have to run the query in a different way?


title: Converting Arcobjects in VB.Net to C#

text: I am attempting to learn ArcObjects using a tutorial that uses VB.NET and I am trying to convert it to C#. I have Visual C# 2008 Express Edition. This is an ArcMap Add-In project that creates a button that displays a message box with the name of the top most layer in the data frame. Here is the VB.NET code:               Dim pMxDoc As IMxDocument     pMxDoc = My.ArcMap.Application.Document     Dim pMap As IMap     pMap = pMxDoc.FocusMap      Dim pFLayer As IFeatureLayer2     pFLayer = pMap.Layer(0)          Dim strDisplayField As String     strDisplayField = pFLayer.DisplayField          MsgBox(strDisplayField, vbOKOnly, "Test")      My conversion in C# so far:               using ESRI.ArcGIS.ArcMapUI:     using ESRI.ArcGIS.Carto;          IMxDocument pMxDoc = ArcMap.Application.Document as IMxDocument;     IMap pMap = pMxDoc.FocusMap;     IFeatureLayer2 pFLayer = pMap.      There is no 'Layer' property. Only 'LayerCount' to choose from. I even tried this:               IFeatureLayer2 pFLayer = ArcMap.Application.Document as IFeatureLayer2;     pFLayer = pMap.      this didn't work either. I have looked all over the ArcMap Object Diagram and I can't find the FeatureLayer2 class anywhere. According to ArcObjects SDK for .NET there is a Layer member for the IMap interface, but I don't have access to that. Any help would be grealty appreciated.


title: When digitizing, objects are not saved well

text: I've some problems with digitizing. If I draw a line or a polygon, the line or polygon is saved in the attribute list, but not visible and its location is not saved either (when I zoom to the object it goes to coordinates 0,0). Both the project CRS and layer CRS are the same. Sometimes it happens, sometimes not, and I don't have any idea what I might be doing differently. If it happens, most of the time QGIS crashes when I try to save the changes. Thanks, Silke I'm using QGIS, Lisboa version


title: MAP Director for ArcGIS

text: I have installed the Cartegraph 8, where I would like to know if it is possible to add the MAPDirector for ArcGIS Extension to ArcMap. As per the following documentation, it seems like the extension should be available from the Customize > Extensions > Map Director for ArcGIS. Unfortunately I was not able to get the extension displayed in ArcMap. Cartegraph Documentation: http://ua.cartegraph.com/ua/Content/MapSolutions/MapArc/ArcMapExtension_Task.htm I could also find the following ESRI link, however I was not able to get the how to steps to add this extension to ArcMap. ESRI Documentation: http://www.esri.com/news/arcnews/spring03articles/business-partners.html I am suspecting if this is something to do with the registration of the DLL. However when I went to the C:\ProgramFiles (x86)\Cartegraph\, I could only see a system folder, where I can see many folders such as bin, Utils, etc. However I am not sure which DLL that I would need to register to make the extension available in ArcGIS. ArcGIS Version: 10.1 Cartegraph Version: 8 OS: Win 7 64-bit Please suggest me on this issue.


title: VRML export issue ArcScene

text: when I try to export simple Tins (or documents which contain Tins) from my ArcScene to .wrl the export fails. Please see my example: I have the following table (xlsx):               X   Y   hoehe          28  -2  0     -2  -2  0     28  48  0     -2  48  0      I import it into ArcScene through add xy data. Than I build a Tin from it using hoehe as height field (There is no difference between using python or going manually):               arcpy.CreateTin_3d(r'C:\Users\aasen\Documents\ArcGIS\Default.gdb\floor_tin', '#',     'Tabelle1$ Events' hoehe Mass_Points <None>", 'DELAUNAY')      Doing this I get a nice flat Tin. Now I'm exporting the Tin with 3D -> VRML where I select the options: Indent text [check]   visible layers [check]   GeoVRML [no check]   Center [no check]   ElevationGrid [check]   JPEG Quality = high The resulting VRML (attached as unbenannt.rar) is incomplete. freeWRL is giving me the following error code:               FreeWRL expected library version freeWRL version 1.22.12, got libfreeWRL version      1.22.12...     windir path=[C:\Windows]     get_current_dir returns[C:/Users/aa/Desktop/]     argv[0]=C:\Program Files (x86)\freeWRL\freeWRL\freeWRL.exe     argv[1]=C:\Users\aa\Desktop\Unbenannt.wrl     posixlycorrect=[(null)]     starting createWindow32     wm_create     starting createcontext32     got hdc     Depth Bits = 24     created context     made a window     resize_GL     showed window     updated window - leaving createwindow     opengl version=3.3.0     ERROR: Expected a closing brace after fields of a node; at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."error in parser_fieldValue by call 2     Expected sfnode Value for a fieldtype! at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."ERROR: Expected a closing brace after fields of a node; at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."error in parser_fieldValue by call 2     Expected sfnode Value for a fieldtype! at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."ERROR: Expected a closing brace after fields of a node; at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."ERROR:Expected "]" before end of MF-Value at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."ERROR: Expected a closing brace after fields of a node; at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."ERROR:Expected "]" before end of MF-Value at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."ERROR: Expected a closing brace after fields of a node; at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."ERROR:Expected "]" before end of MF-Value at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."ERROR: Expected a closing brace after fields of a node; at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."ERROR:Expected "]" before end of MF-Value at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."ERROR: Expected a closing brace after fields of a node; at: "400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000     ..."Parsing complete, but have unrecognized data at end of input:"400                  diffuseColor     1,000 1,000 1,000                  emissiveColor    0,000 0,000 0,000                  shininess        1,000                  specularColor    0,000 0,000 0,000      ....)?¤"Parser Unsuccessfulresize_GL      (This post can also be found on http://forums.arcgis.com/threads/77817-VRML- export-issue) Thank you so much for your help!


title: plotting a map of north greenland in R

text: My knowledge in GIS and R is again tested when I want to map a small part of northern Greenland (called Washington Land) in R (and plot analysis data on top later). I've managed to plot it with               library(maps)     library(mapdata)     library(rworldmap)     map("worldHires", "Greenland", xlim=c(-67,-63), ylim=c(78,82))     newmap <- getMap(resolution = "high")     plot(newmap, xlim=c(-67, -63), ylim=c(78,82), asp=1)      But resolution and projection is horrible for this area and my need. I have a good map in MapInfo, but how should I output a map, so it can load in R and R can plot data accurately on top of it?


title: Enlarge region/features in shape file

text: I'm completely new to QGIS. I have a shape file over the counties in Denmark (see here for what the shape file looks like). The counties around Copenhagen are rather small and I'd like to add an enlarged/"zoomed in" mini map to make it easier to distinguish between the counties. Thus, both the original counties and the "zoomed in" counties could be linked with the same name to make it easy to see the information in the counties. I know this will produce a map that is not geographically correct, but it will be used to convey information. I have installed the plugin that allows me to scale a set of selected features but then the remaining regions/features of the map disappeared. Can anyone point me in the right direction.


title: ArcMap - Opening the attribute table gives 'Could not load data from the data source' error

text: ArcMap wont display the attribute table of a dataset. I get an error 'Could not load data from the data source' and it says 'An invalid SQL statement was used'. It then opens a table with the correct number of fields, but no field headers and no data. It does however display the attribute table if I make a selection and ask it to show only the selection, and I can select every record in the table and it works. It also displays the values if I use the identify button. The datasource is a table in a Geodatabase and was created from FME. FME reports no errors in loading the data, and loads all of the records it is given. Using ArcGIS 10.1 SP1. The table has no special characters in the table name, or the field names. Edit 10/05/2013 The Geodatabase will open in QGIS. If I export to another features class in a file geodatabase using ArcCatalog it the works. I have tried using the FILEGDB and GEODATABASE_FILE writers in FME and the attribute table will not open in ArcMap. Both work in QGIS 1.8


title: xy tolerance from c#.net

text: For some graph work, I need the value of X,y tolerance. I found ISpatialReferenceTolerance.XYTolerance Property. But I couldn't find any option that how can I get `ISpatialReferenceTolerance` interface from the map? Can anybody give me some pointer? Thanks.


title: Difference in vector / raster area measurements

text: I'm comparing raster and vector distance zones. So I created a line, calculated the Euclidean distance, and reclassified the raster to two zones: from 0-10m and 10-20m. To compare it with vector distance, I created a Multiple Ring Buffer with the same line and the same two zones, 0-10m and 10-20m. When I'm calculating the area, for example for the 0-10m zone, the (vector) buffer has an area of 1475 m² while the Euclidean distance raster for the same zone has an area of 1493 m²! Can someone please explain why those two values differ?


title: OGR UnionCascaded gives strange lines (flat polygons)

text: I'm using OGR to perform a union using UnionCascaded. All my polygons have been created by GDAL polygonize which ensures that all polygons follow pixel edges and should therefore align nicely. But I get these strange flat polygons looking like lines: ![Example of union result with strange flat polygons](http://i.stack.imgur.com/8G4eu.jpg) They are positioned at locations where two polygons met before they where unioned. **Can you recognise this pattern and what lead to it?**


title: Why does new feature created not snap to existing feature vertices in ArcGIS 10?

text: I have created a polygon feature earlier then I wanted to create points on the polygon vertices but it won't snap on top of the vertex. I wanted to know the reason why it does not snap and is there any ways to convert the vertices into points instead of manually doing it? Thanks


title: How to limit GEOjson data to specific area?

text: I am currently working with a shapefile from the United States Drought Monitor that countains all drought data from all of the states in the US. I am converting it to geojson with ogr2ogr and would like to filter the geojson output so that it only contains data within the bounds of California. How can I do this? I am wondering if there are optional arguments I can pass in to something like `ogr2ogr -f GeoJSON drought.json USDM_20140617.shp` Or would I only be able to limit the results by doing a join between shapefiles (for example, a join with a shapefile that contained data for California's boundaries)? Shapefiles are here Thanks!


title: FindTask result returns an array of features but only displays the last item in the popupwindow

text:
executing the complete findtask url generated by my application from the address bar I get back a good json object that returns 112 features I was trying to implement functionality similar to this sample http://developers.arcgis.com/javascript/samples/popup_sidepanel/. If you click right on the border or two states a pager appears in the upper left corner. I can't see why the handlers for SelectionChange and SetFeatures in the InitializeSidebar Method aren't performing like the sample.               var map, myPopup, UtilTemplate;          //modules     // @formatter:off     require(["esri/map", //map             "esri/layers/ArcGISDynamicMapServiceLayer",  //map             "esri/graphic", //map             "esri/geometry/Extent", //map             "esri/InfoTemplate", //map             "esri/layers/FeatureLayer",  //map             "esri/layers/Field",  //map             "esri/layers/LabelLayer", //map             "esri/layers/LayerInfo", //map             "esri/layers/ImageParameters", //map                  "esri/tasks/query", //map             "esri/tasks/QueryTask", //map             "esri/tasks/FeatureSet", //map             "esri/tasks/FindTask",             "esri/tasks/FindParameters",                  "esri/tasks/IdentifyParameters", //map             "esri/tasks/IdentifyTask", //map             "esri/tasks/IdentifyResult", //map                  "esri/symbols/SimpleFillSymbol", //map             "esri/symbols/SimpleLineSymbol", //map             "esri/symbols/SimpleMarkerSymbol", //map             "esri/symbols/TextSymbol", //map             //"esri/renderers/Renderer",                     "esri/renderers/SimpleRenderer", //map                       "esri/symbols/Font", //map                       "esri/dijit/PopupTemplate", //map             "esri/dijit/HomeButton", //map             "esri/dijit/Basemap",              "esri/dijit/BasemapLayer",             "esri/dijit/BasemapGallery",             //"esri/dijit/BasemapToggle", //map             "esri/dijit/Geocoder", //map             "esri/dijit/Scalebar", //map             "esri/dijit/Legend", //map             "esri/dijit/OverviewMap",             "esri/dijit/Measurement",             "esri/dijit/Popup",                  "agsjs/dijit/TOC",  //map                  "dojo/ready", "dojo/parser", "dojo/on", "dojo/dom", //all             "dojo/dom-construct", "dojo/number", "dojo/query", "dojo/fx", //"dojo/Deferred", "dojo/promise/Promise", "dojo/promise/all",                  "dojo/_base/Color", //all             "dojo/_base/declare", //all             "dojo/_base/array", //all             "dojo/_base/lang", //all                       "dijit/layout/BorderContainer", "dijit/layout/ContentPane",                  "dijit/registry", "dijit/form/Button",             "esri/graphicsUtils", "dojo/_base/connect", "esri/domUtils"                                 ],              function(Map, ArcGISDynamicMapServiceLayer, Graphic, Extent, InfoTemplate,             FeatureLayer, Field, LabelLayer, LayerInfo, ImageParameters, Query, QueryTask, FeatureSet, FindTask, FindParameters, IdentifyParameters, IdentifyTask, IdentifyResult,             SimpleFillSymbol, SimpleLineSymbol, SimpleMarkerSymbol,TextSymbol, SimpleRenderer, Font, PopupTemplate, HomeButton, Basemap, BasemapLayer, BasemapGallery, Geocoder, Scalebar, Legend, OverviewMap, Measurement, Popup,             TOC, ready, parser, on, dom, domConstruct, number, query, fx, Color, declare, arrayUtils, lang,             BorderContainer, ContentPane, registry, Button, graphicsUtils, connect, domUtils) {          // @formatter:on              ready(function() {     //Setup Popup             var myPopup = new Popup({fillSymbol: new SimpleFillSymbol(SimpleFillSymbol.STYLE_NULL, new SimpleLineSymbol(SimpleLineSymbol.STYLE_SOLID, new Color([102,0,204]), 5), new Color([255,255,0,0.25])),                  lineSymbol: new SimpleLineSymbol(SimpleLineSymbol.STYLE_SOLID, new Color([255,0,0]), 5),                 titleInBody: false, offsetX: 100, offsetY: 100,                 pagingInfo: false                 },  domConstruct.create("div")); //1st color fill is for line, 2nd for fill                             var extentInitial = new Extent({             "xmin" : -9230904.490170034,             "ymin" : 3111585.686814823,             "xmax" : -9086285.632654676,             "ymax" : 3173576.116754018,             "spatialReference" : {                 "wkid" : 3857             } // center and zoom no longer needed         });         var myLods = [{"level" : 0, "resolution" : 4891.96981024998, "scale" : 18489297.737236},  //levels 5 - 19 of world street map; 0-14 index levels for center and zoom                 {"level" : 1, "resolution" : 2445.98490512499, "scale" : 9244648.868618},                 {"level" : 2, "resolution" : 1222.99245256249, "scale" : 4622324.434309},                 {"level" : 3, "resolution" : 611.49622628138, "scale" : 2311162.217155},                 {"level" : 4, "resolution" : 305.748113140558, "scale" : 1155581.108577},                 {"level" : 5, "resolution" : 152.874056570411, "scale" : 577790.554289},                 {"level" : 6, "resolution" : 76.4370282850732, "scale" : 288895.277144}, //initial extent zoom                 {"level" : 7, "resolution" : 38.2185141425366, "scale" : 144447.638572},                 {"level" : 8, "resolution" : 19.1092570712683, "scale" : 72223.819286},                 {"level" : 9, "resolution" : 9.55462853563415, "scale" : 36111.909643},                 {"level" : 10, "resolution" : 4.77731426794937, "scale" : 18055.954822},                 {"level" : 11, "resolution" : 2.38865713397468, "scale" : 9027.977411},                 {"level" : 12, "resolution" : 1.19432856685505, "scale" : 4513.988705},                 {"level" : 13, "resolution" : 0.597164283559817, "scale" : 2256.994353},                 {"level" : 14, "resolution" : 0.298582141647617, "scale" : 1128.497176}                 ];                                //setup event handlers for the next/previous buttons                 on(dom.byId("previous"), "click", selectPrevious);                 on(dom.byId("next"), "click", selectNext);                                 map = new Map("cpCenter", {                 basemap : "hybrid",                 extent : extentInitial,                 lods : myLods,                 showAttribution: false,                 infoWindow: myPopup,                 displayGraphicsOnPan: false,                      //nav: true             });                  map.on("load", function () {               //Create the find parameters              findParams = new FindParameters();               findParams.returnGeometry = true;               findParams.layerIds = [0];               findParams.searchFields = ["ACCOUNT","Name1"];               findParams.outSpatialReference = map.spatialReference;               console.log("find sr: ", findParams.outSpatialReference);              });                   //map.infoWindow.set("popupWindow",false);              map.infoWindow.popupWindow = false;               initializeSidebar(map);               function initializeSidebar(m){                     var popup = m.infoWindow;                          //when the selection changes update the side panel to display the popup info for the                      //currently selected feature.                      connect.connect(popup, "onSelectionChange", function(){                         displayPopupContent(popup.getSelectedFeature());                     });                          //when the selection is cleared remove the popup content from the side panel.                      connect.connect(popup, "onClearFeatures", function(){                         //dom.byId replaces dojo.byId                         dom.byId("featureCount").innerHTML = "Click to select feature(s)";                         //registry.byId replaces dijit.byId                         registry.byId("leftPane").set("content", "");                         domUtils.hide(dom.byId("pager"));                     });                          //When features are associated with the  map's info window update the sidebar with the new content.                      connect.connect(popup, "onSetFeatures", function(){                         displayPopupContent(popup.getSelectedFeature());                         dom.byId("featureCount").innerHTML = popup.features.length + " feature(s) selected";                              //enable navigation if more than one feature is selected                          popup.features.length > 1 ? domUtils.show(dom.byId("pager")) : domUtils.hide(dom.byId("pager"));                     });                 }                function displayPopupContent(feature){                     if(feature){                         var content = feature.getContent();                         registry.byId("leftPane").set("content", content);                     }                 }                      function selectPrevious(){                     map.infoWindow.selectPrevious();                 }                      function selectNext(){                     map.infoWindow.selectNext();                 }                   var home = new HomeButton({                 map: map             }, "HomeButton");             home.startup();              var overview = new OverviewMap({map: map, opacity: 0.5, attachTo: "bottom-left"});             overview.startup();                   var measurement = new Measurement ({map: map},"measurementDiv");             measurement.startup();                  connect.connect(measurement.area, "onClick", measureButtonClickHandler);             connect.connect(measurement.distance, "onClick", measureButtonClickHandler);             connect.connect(measurement.location, "onClick", measureButtonClickHandler);              function measureButtonClickHandler() {                  if (this.checked){                  myPopup.clearFeatures();             lyrParcels.disableMouseEvents();                      }                   else if (this.checked == false) {                      lyrParcels.enableMouseEvents();              }              }                        var scalebar = new Scalebar({             map: map,             scalebarUnit : "dual",             //attachTo: "botttom-left"         });               /*InfoTemplates*/ //Info Templates don't access domain descripts, popup templates do//              var wwwParFlds = ["ACCOUNT","Name1","Land_Val","Improvemt","Just"]; //"ACCOUNT","PID"                        var plyParcels = new InfoTemplate();             plyParcels.setTitle("Parcel Info");             plyParcels.setContent("Account Number: ${ACCOUNT}
Land Value:  ${Land_Val:NumberFormat}
Improvement Value:  ${Improvemt:NumberFormat}
Just Value:  ${Just:NumberFormat}
Market Area:
Owner Name:  ${Name1}");                                              //"
 View Property Detail");               /*InfoTemplates*/                                //STYLE PARCELS         var renderer = new SimpleRenderer(                new SimpleFillSymbol(SimpleFillSymbol.STYLE_SOLID,         new SimpleLineSymbol(SimpleLineSymbol.STYLE_SOLID,         new Color([255,0,0]), 2),new Color([255,255,0,0.1])              ));                                  var lyrParcels = new FeatureLayer("https://ags2.scgov.net/arcgis/rest/services/ScpaInternal/AppraiserWM/MapServer/0", {             mode : FeatureLayer.MODE_ONDEMAND,             opacity : 0.9,             infoTemplate : plyParcels,                  maxAllowableOffset: calcOffset(),             outFields : wwwParFlds             });                  lyrParcels.setMinScale(2256.994353);             lyrParcels.setRenderer(renderer);              var lyrBuildClass = new FeatureLayer("https://ags2.scgov.net/arcgis/rest/services/ScpaInternal/scpaAppraiserGIS_WM/MapServer/26", {             mode : FeatureLayer.MODE_ONDEMAND,             visible : false,             opacity : 0.8,             maxAllowableOffset: calcOffset(),                  });                             //  lyrParcels.setVisibility(true);                         /*calc feature layer offsets*/         on(map, 'zoom-end', function() {                      lyrParcels.maxOffset = calcOffset();                      lyrParcels.setMaxAllowableOffset(lyrParcels.maxOffset);                  });                  function calcOffset() {                 return (map.extent.getWidth() / map.width);             }                   /*add the layers!*/              map.addLayers([lyrParcels,lyrBuildClass]);          /*find by pid*/              registry.byId("search").on("click", doFind);              findTask = new FindTask("https://ags2.scgov.net/arcgis/rest/services/ScpaInternal/AppraiserWM/MapServer/");                         function doFind() {               //Set the search text to the value in the box               findParams.searchText = dom.byId("searchText").value;               findTask.execute(findParams, showResults);             }                  function showResults(results) {               //This function works with an array of FindResult that the task returns               map.graphics.clear();               var symbol = new SimpleFillSymbol(                 SimpleFillSymbol.STYLE_SOLID,                  new SimpleLineSymbol(SimpleLineSymbol.STYLE_SOLID, new Color([98, 194, 204]), 2),                  new Color([98, 194, 204, 0.5])               );                    //create array of attributes               var items = arrayUtils.map(results, function (result) {                 var graphic = result.feature;                 graphic.setSymbol(symbol);                 map.graphics.add(graphic);                           return result.feature.attributes;               });                        map.setExtent(graphicsUtils.graphicsExtent(map.graphics.graphics),true);                                                                                }          /*Find by PID*/               /*GeoCoder*/         var locatorUrl = "https://ags2.scgov.net/arcgis/rest/services/CompositeLocator/GeocodeServer";           var scgGeocoder = [{             url: locatorUrl,             name: "CompositeLocator",             singleLineFieldName: "SingleLine",             placeholder: "Find an Address"             }];              var geocoder =  new Geocoder({             map: map,             autoComplete: true,             autoNavigate: false,             arcgisGeocoder: false,             geocoders: scgGeocoder,                  },dom.byId("divSearch")); //dom.byId: if not present, will not clear text from text box, or selection graphics             geocoder.startup();                  geocoder.on("select", showLocation);             geocoder.on("clear", removeInfoGraphics);     /*GeoCoder*/              function showLocation(evt) {                  map.graphics.clear();             var point = evt.result.feature.geometry;                   var ptSymbol =  new SimpleMarkerSymbol(SimpleMarkerSymbol.STYLE_SQUARE, 10,         new SimpleLineSymbol(SimpleLineSymbol.STYLE_SOLID,         new Color([255,0,0]), 1),         new Color([0,255,0,0.25]));                            var graphicPoint = new Graphic(point, ptSymbol);             map.graphics.add(graphicPoint);                       var queryParcels = new Query;             queryParcels.geometry = evt.result.feature.geometry;                       var parcelsQuery = lyrParcels.selectFeatures(queryParcels,FeatureLayer.SELECTION_NEW); //deferred evacQuery                  map.infoWindow.setFeatures([parcelsQuery]);  //map.infoWindow.setContent(evt.result.name);             map.infoWindow.show(evt.result.feature.geometry);                  // Center and zoom the map on the result             if (point !== undefined) {                 map.centerAndZoom(point, 13); //7 for sph             }          //      lyrFranZones.on("selection-complete", returnNull); // call getDomainVals to get the domain values from the attribute fields for data grid or an infoTemplate              }         function returnNull() {             var retVal = null;              }               /* legend */         map.on("layers-add-result", function(results) { // sets up an array to store layer index positions and returns the layer legend values and layer title name         var toc = TOC({             map: map,             layerInfos: [{layer: lyrParcels, title: "Parcels"},{layer:lyrBuildClass, title: "Building Class", slider: "true", collapsed: "true"}]         }, "tocDiv");         toc.startup();         });          /* legend*/          /* CleanUp*/         function removeInfoGraphics() {             map.infoWindow.hide();             map.graphics.clear();             //var clearStore = new ItemFileReadStore({data: {identifier: "",  items: []}});             //gridEvac.setStore(clearStore); //clears the grid, kudos Michael!         };     /* CleanUp*/              }); //brace for on ready function     }); // brace for require and Aliasees function


title: ArcMap 10.1 Spatial Join: 'Contains' and 'Within' question

text: I have two polygon datasets: electoral districts and postal code areas. I'm performing a spatial join to see which postal code areas are completely within a district (the area can touch the borders, but doesn't cross them; i.e. not situated in 2+ districts). The way I see it I can do it two ways: (1) Districts as target features, postal code areas as join features with a Contains Clementini match option (2) Postal code areas as target features, districts as join features with a Within Clementini match option I assume these are the same thing. Yet I have ran both and they yield different results. Theoretically: why am I getting different results? And which way is deemed "more correct"? Thanks for your input in advance!


title: How to change delimiter in attribute table of QGIS?

text: I have imported a .kml file into QGIS, and the data is all shoved together into a single field like this: ![delimeter](http://i.stack.imgur.com/uaF3a.png) Fine, it's pipe delimited - easy to change the delimiter like in LibreOffice calc imports I thought. But could not find a solution within QGIS. Next attempt: change the delimiter in the kml file to the common default comma delimited using Sed: `sed -e 's/|/,/g' in.kml > out.kml` Still no joy - it looks exactly the same as above, only with commas instead of pipes. So someone please help me, and I'm sure it will be of use to others. (P.s. a colleague and I have tried the same thing in ArcGIS and the same problem is encountered).