title: Automated/geoprocesing tool to slice, clip, or cut polygons using polylines using ArcGIS Desktop?

text: Looking for a simpler method here. I am trying to split up / cut / slice existing polygons by using existing poly-lines. An example would be splitting a water body or land parcel at the point where a bridge/road crosses it. But the road network may not necessarily result in a closed loop. Since the poly-lines are not necessarily interconnected or continuous, creating a polygon from them is not an option (which eliminates using the split tool). Also, I have tried using a topology with the geometry, but it keeps failing, probably due to the large/complex geometry. Current Workflow: I have accomplished this with the Feature to Polygon tool, combining the lines and polygons together, but it creates additional unwanted polygons anywhere a loop closes. I created a mask of the original polygons, and used this with the overlay->erase tool to remove the unwanted polygons. This still leaves some unwanted geometry (mostly slivers), but is somewhat workable. This seems like an extremely convoluted and round-a-bout way to accomplish (what seems like it should be) a very simple task. Other than manual edits, or using a topology, is there is tool that could accomplish this in a single step? Using: ArcMap \ ArcInfo Desktop 10 SP5 Edit 1: In my case, it is not actually roads as mentioned above. I have water polygons for a coastal area, and need to split the polygons where rock dams or levees have been put across the waterways. Which are typically not interconnected. The water polygons have been simplified and repaired, to the point where I would no longer call the data "dirty" just complex and large. But I've got the workaround mentioned above working for this case. I'm looking more "in general" for a tool that can simply divide polygons using polylines. Edit 2: Mapperz:Thanks for the Model Builder suggestion. I am going to use that as a stop gap solution, for now. Jakub: Thanks for the suggestion. I'm not opposed to a programmatic solution or developing a custom tool, although I've never made one before. I have programming experience, but not in conjunction with Arc. I'd prefer something that directly cuts the geometry though, instead of following the roundabout logistics above. In theory that should reduce resultant gaps, since it wouldn't be subject to multiple iterations of cracking/clustering. Although, I'm not sure that would be as easy or even possible. Edit 3: I'm looking for something to function like on the left in the image. ![Cut Polygon vs Feature to Polygon](http://i.stack.imgur.com/t7zeX.png)


title: Trilateration 3D - Vincenty's Formula

text: with reference to Trilateration using 3 latitude and longitude points, and 3 distances. The accepted answer was votes as shown in wikipedia at: http://en.wikipedia.org/wiki/Trilateration in conjunction with a conversion to Radians before and after. Is there a way to incorporate Vincenty's formulae http://en.wikipedia.org/wiki/Vincenty%27s_formulae to this process to bring a more accurate answer? Kind Regards


title: Holding and using Romanised names with Diacritics

text: I have two True type font files (Permanent Committee on Geographical Names) PCGNA.TTF AND PCGNTN.TTF. The romanised diacritics are well read in MS excel. But when I convert this excel in shp or pgdb feature class, it becomes unreadable. Does anyone know how to overcome this problem?


title: Write string to status bar that doesn't go away

text: I'm attempting to write the number of items matching a given query to the statusbar. when i debug through the code, i am able to step through and see the message hit the statusbar, but then it immediately disappears. Below is the function in question.               public class CollisionViewExtension : Extension     {          private static IQueryable<Collision> AllCollisions         {             get { return _allCollisions ?? (_allCollisions = Database.Query<Collision>(Collision.Query).ToList().AsQueryable()); }         }         ...         public static void DrawCollisionPoints(string whereClause)         {             var matches = new List<Collision>();             try             {                 matches = AllCollisions.Where(string.IsNullOrWhiteSpace(whereClause) ? "1=1" : whereClause, null).ToList();                      if (CurrentSpatialQueryShape != null)                 {                     matches = matches.Where(x => IsPointInPolygon(CurrentSpatialQueryShape, (IPoint)x.Point)).ToList();                 }             }             catch (Exception)             {                 System.Windows.Forms.MessageBox.Show("Invalid query");             }                  FeatureInserter.InsertFeatures(matches);                  var layer = new FeatureLayerLocator(ArcMap.Document.FocusMap).GetFeatureLayerByModelName("CollisionView-Collisions");             if (layer == null)             {                 ArcMap.Document.FocusMap.AddLayer(FeatureLayer);             }                  ArcMap.Application.Caption = string.Format("{0} collisions", matches.Count);             Application.DoEvents();             ArcMap.Application.StatusBar.set_Message((int)esriStatusBarPanes.esriStatusMain, string.Format("{0} collisions", matches.Count));                  ArcMap.Document.ActiveView.PartialRefresh(esriViewDrawPhase.esriViewGeography, FeatureLayer, null);             ArcMap.Application.StatusBar.set_Message((int)esriStatusBarPanes.esriStatusMain, string.Format("{0} collisions", matches.Count));                  CurrentQuery = whereClause;         }         ...     }      I feel like this is pretty much how I've seen it used in all the examples, but for some reason my message is not displayed. Asked another, probably simpler way, is there any way to show the user a count of the features in a particular layer in the statusbar?


title: View PostGIS geometry in QGIS

text: I am a newbie to GIS systems so please be patient. I have a PostGIS database set up with a spatial data table. One of the columns is of Geometry type and contains Points with SRID 4326. I want to view these points overlaid on a map. I can connect to PostGIS from QGIS, load the points as a PostGIS layer and using the OpenLayers plugin also overlay these on a map. The trouble is, it appears that the latitude and longitudes are reversed when displaying them in QGIS! When I connect to the database from command line and use the SELECT query with ST_ASTEXT on the geometry column, I can see that the lat longs are in correct order. Any suggestions or pointers on what might be the problem?


title: How does Flow Mapper Plugin work?

text: I downloaded the Flowmapper plugin for Qgis and I can't seem to find any information on how it is supposed to work... Also the text seems to cut-off and not fully showing for the input fields so I'm not quite sure what the fields should be anyway... Here is an example of what I mean: ![enter image description here](http://i.stack.imgur.com/vGwBC.png) I tried to select a file storing my node coordinates but am not allowed to (greyed-out) pick files that are .csv, .dbf, .shp or seemingly any of the filetypes on my computer... Which is why I started looking for some information on how to use the Flow Mapper plug-in, after finding nothing descriptive I thought it would be easier to ask someone who has maybe used it successfully before. For reference I have a point file of centroids of polygons and I have a collection of lat long points which contain an attribute that should connect them to one of the above mentioned centroids. Maybe there is a better way to do this? QGIS 1.8.0 on Mac Thanks


title: Create Shapefile of all overlapping features from single Shapefile

text: I need to create a Shapefile out of another Spahefile that consists all overlaps of features within the initial Shapefile. I'm sure there is a simple solution using ArcMap or QGIS, but at the moment I am just lost and somehow I struggle to find something via search engines although I am aware that this question must have been asked a thousand times. The screenshot shows the desired result. My new Shapefile should only consist the darker parts independent on the number of overlapping features. ![enter image description here](http://i.stack.imgur.com/4QE3b.jpg)


title: Problem with encoding of CSV and SHP files in Linux

text: I am new to using QGIS. I came up against the encoding issue of shp files. This solution to solved my problem in Windows (for source I had to put windows-1253 instead of utf-8) My question is what solution would I use in Linux systems where I have the very same issue? Where can I "insert" a similar choice?


title: python - Calculate field with part of file name

text: So I previously successfully used the sample code here: Model Builder - How to calculate Field using part of filename               # Import standard library modules     import arcpy, os, sys     from arcpy import env          # Allow for file overwrite     arcpy.env.overwriteOutput = True          # Set the workspace directory      env.workspace = r"C:\temp.gdb"           # Get the list of the featureclasses to process     fc_tables = arcpy.ListFeatureClasses()          # Loop through each file and perform the processing     for fc in fc_tables:         print str("processing " + fc)              # Define field name and expression         field = "SID"         expression = str(fc[:5]) #subsets first 5 characters of fc name              # Create a new field with a new name         arcpy.AddField_management(fc,field,"TEXT")              # Calculate field here         arcpy.CalculateField_management(fc, field, "expression", "PYTHON")       But now I'm running it again, with no changes, even with the same datasets I used before, and it no longer works. I cannot see how or why it errors out. Error message: > Runtime error : ERROR 000539: Error running expression: expression : name > 'expression' is not defined Failed to execute (CalculateField). Help!?!


title: Need help assigning values to field using OGR in python

text: I've got a similar question to this one How to programatically populate a shapefile field with polygon areas in OGR?, but the solution doesn't seem to be working for me. I am trying to import a shapefile from PostGIS into python using OGR, add a field, then populate the field with data from a csv based on an ID and the date. When I print the fields for the feature during the loop, the field gets the correct data but it isn't writing to the layer o3_proj (last two print statements). Am I using the SetFeature function incorrectly? I'm new to PostGIS and ogr, so please speak slowly!               import os, ogr, gdal, csv     driver = ogr.GetDriverByName('PostgreSQL')     datasource = driver.Open("PG: host='localhost' port='54321' dbname='chicago' user='postgres'", update = 1)          o3_proj = datasource.GetLayerByName('o3_monitors(geom3528)')          date = []     maxozone = []     aqi = []     aqs_site_id = []          f = csv.reader(open('/ozone/ozone_data.csv','rU'))     x = list(f)          for iii in range(1,len(x)):         date.append(x[iii][0])         aqs_site_id.append(x[iii][1])         maxozone.append(float(x[iii][2]))         aqi.append(int(x[iii][3]))          daily_field = ogr.FieldDefn('temp_o3',ogr.OFTReal)     o3_proj.CreateField(daily_field)          for feat in o3_proj:         for jjj in range(0,len(date)):             if aqs_site_id[jjj] == feat.GetField(0) and date[jjj] == '03/06/2005':                 o3_proj.SetFeature(feat)                 feat.SetField2(6,maxozone[jjj])                 o3_proj.SetFeature(feat)                             print str(feat.GetField(0))+','+str(feat.GetField(1))+','+str(feat.GetField(2))+','+str(feat.GetField(3))+','+str(feat.GetField(4))+','+str(feat.GetField(5))+','+str(feat.GetField(6))                 print o3_proj.GetFeature(feat.GetFID()).GetField(6)      Result of the print statements:               17-031-0064,41.79079,-87.60165,2001,2012,None,0.036     None     17-031-0076,41.7514,-87.71349,2004,2012,None,0.038     None     17-031-4007,42.06028,-87.86323,2002,2012,None,0.035     None     17-031-4201,42.14,-87.79922,2001,2012,None,0.031     None


title: How to get the coastline from OSM?

text: I've successfully imported all data for Sweden to my PostGIS database. Right now I'm designing my map in Tilemill with OSM Bright as base. But I've got a major problem: where is the coastline? Right now all the roads and buildings are on the water. I've looked through all three `landusages` tables, and the only relevant geometry types are `land` and `island`, but they aren't what I expected. I'm a GIS rookie, so any help is appreciated. :-)


title: Locating 0-points multipart polygons using arcpy

text: By Running this piece of code:               import arcpy     arcpy.env.workspace = "C:\working\GEODATABASES\\tests.gdb"     fcList = arcpy.ListFeatureClasses()     for fc in fcList:         with arcpy.da.SearchCursor(fc, ("OID@", "SHAPE@")) as cursor:         for row in cursor:             print("Feature {0}".format(row[0]))             print "   " + str(row[1].partCount) + " parts for THIS PARTICULAR FEATURE"             for part in row[1]:                 poly = arcpy.Polygon(part)                 print("   Points {0}".format(poly.pointCount))        In theory I am able to list all polygons for a particular feature class and list whether they are multi-part polygons or not. In addition, it shows the number of points per feature part. HOWEVER, it does not do the job properly as I am dealing with a complicate feature class that actually have more feature parts than the script can identify. Even in some of them, it does not obtain the right number of points per part either. By using the OLD reading cursor system in arcpy (not the above one), I get in some polygons a different number of points, but still there are some of them missing, therefore the script cannot access to them. Why is that this script stops counting either points or parts sometimes when entering a multipart-polygon within my feature class? Is it my code not well tuned for complicated multi-part polygons? Thanks,


title: How to install Proj4js.js in my localhost directory?

text: I'm developing a GIS website for vineyard management. I am using Openlayers with the WMTS IGN (French institute for geography) datas. It works fine. I also would like to add the possibility for the user to switch to a geology map provided by the BRGM (geological french institute) geoservice. The issue is that IGN datas avec in Webmercator projection (EPSG:3857) and BRGM datas avec in EPSG:27572 (NTF (Paris)) and i would like to keep map center and bounds while switching from one to the other one. After reading a lot on the internet, i found i needed to use Proj4js for projection transformation. I also have vestor shapes (vineyard fields) and i would need to visualize in both projections by transforming them. I am quite used to implement js librairies (by copying js files in my localhost and include it in the header of my php/html file) but in the case of proj4js, it looks different. The gituhub page only gives this information for installation :               npm install proj4     bower install proj4     jam install proj4     component install proj4js/proj4js      I have no idea how it works and i can't get any clear information from what i read. In the github files, there is no proj4js.js file... Is there somewhere a clear tutorial for installation ? Does someone have experience on this installation ? I tried a "classical" installation downloading 1.1.0 from this link http://download.osgeo.org/proj4js/proj4js-1.1.0.zip but it seems to be an old version and i'm getting errors with the test code below : "lcc:forward: llInputOutOfRange: NaN : NaN" Here is my test code :               <head>          <script src="js/proj4js/lib/proj4js-compressed.js"></script>     <script src="js/proj4js/lib/defs/EPSG3857.js"></script>     <script src="js/proj4js/lib/defs/EPSG27572.js"></script>     <script src="http://openlayers.org/api/2.13.1/OpenLayers.js"></script>          </head>          <body>          <script type="text/javascript">               var source = new Proj4js.Proj("EPSG:3857");     var dest = new Proj4js.Proj("EPSG:27572");          Proj4js.reportError = function (msg) {         throw new Error('Whoops: ' + msg);     };          var position1 = new OpenLayers.LonLat(0, 5700000);              console.log(source);         console.log(dest);         console.log(position1);              console.log(Proj4js.transform(source, dest, position1));               </script>               </body>      Thanks a lot, Vincent


title: Does OSM data contain the direction of travel?

text: I was looking for **Direction of Travel (DOT)** information inside OpenStreetMap shapefiles, but only I've just found **"ONEWAY = 1 or 0"** , but nothing about its direction. If I work on an OSM editor, I can see its arrows which declare its direction. But as you can see, once downloaded as Shapefile type, there is no DOT information about. thanks for your time, ![OSM Metadata](http://i.stack.imgur.com/k0osF.png)


title: Is there a way to use the old command line in Arcgis10?

text: I know that the command line has been replaced by the python window, but I was just wondering if there is any way to use command line within Arc 10? Thanks


title: Is there a way to configure ArcMap to persist the HTML popup window size?

text: I'm specifically asking for when the "Display HTML formatting" is configured as "As a table of the visible fields.", but I think all options behave the same way. My users click the HTML Popup tool![alt text](http://i.stack.imgur.com/5KWDh.png), click the feature on the map to display the HTML popup window, and resize the window until the get the desired size. If they close that window, and repeat the steps...the windows reverts back to the original size and they must resize each time. They feel, and I agree with them, that the window should open up to be the same size as when it was closed. Is this an ArcMap configuration? Or is there a registry value? Is there a way to affect window size if you format using an XSL template?


title: Visualization to inference workflow

text: I provide statistical support for a public health department. As you might imagine, we put together a lot of maps on a regular basis. For me, maps are just another kind of data visualization - useful for getting a feel for the data, for generating and checking hypotheses, etc. But we don't often follow through on actual modeling and hypothesis **testing**. How do you/your organization go about this? What does a workflow that includes inference look like? Who's involved? What tools do you use? What would it _ideally_ look like, if you had your way? Thanks! **EDIT** To be clear, I'm curious about different strategies for going from spatial data to formal, statistical tests of hypotheses about what's going on in the world. For example, let's say I'm trying to target an educational campaign to increase tuberculosis testing. I (personally) would map out the cases of TB against covariates of interest (say, median income or percent foreign-born residents) and try to see if there were any patterns. I might or might not find any; but I would ultimately build a model to estimate the association between those covariates and the number of demographics. This is a critical step because of how good humans are at finding patterns where none exist, or finding uninteresting ones. I know how to do this on my own, but I'm curious about how different organizations institutionalize it (if at all).


title: Convert between ESRI Geometry and WKT

text: It seems like there is no means of converting an ArcObjects geometry to the Well-Known Text representation (and vice versa) in ArcGIS ArcObjects API. The only thing I was able to find is conversion to WKB (the IWkb interface). Is there a way to perform conversion between geometry objects and WKT or do I have to implement it myself? I primarily aim for .NET implementation without much external dependencies. Thanks.


title: Creating visualizations for geometry points with geotools

text: I am a beginner in everything that involves GIS coding and so on. After doing some research i chose Geotools as the tool to help me with my project. Basically i will need to create a new routing algorithm that will get as an input the coordinates of two points : start and destination and a shape file. I have followed the tutorial from geotools, creating a graph etc but when i want to do a simple visualization of lets say the nodes that i have after the graph creation, for the life of me, i do not find a way to do so. Would it be possible for you to point me in a direction that could solve this mystery for me ?


title: Read an INSPIRE GML in QGIS2.0 : gml:xlink resolving

text: I'm trying to read an INSPIRE GML, which countains internal and external links between the elements. I saw on the GDAL site that we can set environnement value in order to load this link in QGIS (http://www.gdal.org/ogr/drv_gml.html) but it seems it doesn't work on this QGIS version. The impacteds lines look like this one: Do someone have some experience in integration of this type of model? Thank you!


title: Free land use data for Manhattan

text: Does anyone know where I can get free land use data for Manhattan? I'm specifically looking for fraction of residential, commercial and other landuse data for each census tract for the New York County. I did search the NYC - Bytes of Big Apple website and found that they charge an appalling 300$ per borough for the land use at tax lot level. Any help is greatly appreciated...


title: Anyone try generating C# classes from http://schemas.opengis.net?

text: I'd like to generate gml/wfs/wms classes from the OGC schemas. Using xsd2code I receive numerous errors, one of them being "Invalid particle derivation by restriction". I'm curious is anyone has successfully generated .NET classes from the OGC schemas and if so, how?


title: Java library for rendering

text: We have some shapefile data and we tried to render them according to scope (boundingbox) from user. For example, the client may send request like this: layer=river&width=256&height=256&bbox=.... Then I will generate an image accordingly. I know this is simiary to a wms server, however I do not want to use geosever or something else. They are heavy. I just have 2 shapfile files with a fixed style. And I found that geoserver use geotools as the render library, so I tried to use geotools myself, however I found that geotools do not provide a accept performance. So I wonder if there are any alternative solutions?


title: In QGIS are there similare construction tools to ArcGIS?

text: For a project we are using ArcGIS to reshape autogenerated features. For that we often use the construction tools polygon, rectangel, circle, etc in ArcGIS. The right angle is also an often used tool. Are there some similar construction tools in QGIS? I found the CADTool-plugin and the Freehand- plugin, but that seems not to be the best solution.


title: How does IP geocoding work?

text: I'm trying to understand how using each part of an IP address influences the accuracy of geocoding. For example, if `15.201.xxx.xxx` is used instead of `15.201.123.xxx` what is the approximate distance in KM? I'm using the IP2Location Lite database to examine who has visited a website. The `README` states the following: `1\. IP2Location Lite Edition is free package with accuracy up to Class C (192.168.1.X) only. It is restricted for non-commercial use.` However, I'm having difficulty finding documentation on what the accuracy of this actually means. Does this accuracy vary over time also?


title: How to create svg line style in QGIS?

text: I need to draw an telephone or energy lines as is done on the polish-core map. Thelephone line was drown by the polyline. I tried the SVG symbols, when I have a lot of single-lines, for each I 'put' symbol svg on the first and the last vertex. For polylines, unfortunately, line markers located on the central vertexs do not rotate to the next segment, but are suspended 'in the air'. If someone has already tested this, I'd like to know the rules. I enclose a snapshot of my unsuccessful endeavors. ![The effect of mine line style](http://i.stack.imgur.com/s8edC.png) And below is the effect that I need, ( _now its generated by Inkscape_ ) ![The effect that I need](http://i.stack.imgur.com/i9iJ9.png) Thanks for the reply.


title: How to use PostgreSQL/PostGIS with ArcGIS for Desktop?

text: Has anyone successfully used PostgreSQL with PostGIS successfully? I am a semi-experienced ArcGIS user and want to build out a database and connect it to ArcCatalog. I've scoured the internet for resources but nobody has really given a good walkthrough. Any help would be much appreciated!


title: How to upgrade php version for mapguide server

text: I am using mapguide server 2.0 . It has PHP version 5.0.I want to use sqlite database3 but it doesn't support sqlite3. Does any one know how to upgarde php version manually? or way to convert sqlite3 database to sqlite2?


title: Hiding layers in mapfile (Mapserver) from WMS requests

text: We have a Mapserver mapfile used in a web app with a lot of layers. We want to publish only some of them with WMS, keeping the others hidden to any WMS request. It has to be the same mapfile and we need each layer has its NAME attribute specified Is there any way to do this?


title: What is the GeoCloud that everyone is talking about

text: Directions Magazine, Pitney Bowes (MapInfo), and many blogs are talking about 'CLOUD's in relation to GIS. I have read about Cloud-Computing but still don't understand the concept.


title: How to create a smallest quadrilateral concave polygon around GPS points?

text: I have this GPS track in a metric grid (Symbology by elevation, if anyone wonders): ![GPS track](http://i.stack.imgur.com/QmhiZ.png) I need to define it's general area - for registrations and permits, with just 4 points. Therefore, I assume that the corners of a circumscribed quadrilateral polygon (4 lines) will do, as that will be the smallest possible polygon area-wise. I don't know which Arcmap tools I can use to do it, though I imagine the result should look something like this: ![enter image description here](http://i.stack.imgur.com/t81qi.png) With a polygon I can easily know the corners' coordinates, but if there's another way to extract them, I'm all ears. Thanks in advance. **Edit:** This is the resulting convex hull tool: ![enter image description here](http://i.stack.imgur.com/OnTTI.png)


title: Handling XML from within ArcGIS

text: In my company, there is a web service that serves geographic data stored in a remote dababase through an URL. This service makes queries to the database and encodes the required data in XML format. So, to query data from a table (for example **User** ), I do something like: `http://myCompanyName/service/QueryUserTable` So the service returns back to the client's browser a XML with all users stored in that table. Now, I would like to do the same from ArcGIS software: The thing would be connect this desktop application to that remote service, retrieve that XML file and parse it **but all inside ArcGIS environment** , not using other scripts or programs outside it. Is it possible? :S Just in case it is not possible to do what I said, what is the best approach to handle and manage traffic in XML using ArcGIS? I am sorry if I posted in the wrong place or this was answered before, the problem is that I wasnt sure what to ask to Google/Stack Exchange to obtain a good answer. Thanks! :D * * * WOW! It is a good idea and I like the way of retrieving the parameters like this long link you pasted. Didn't know about that module since I am pretty new with ESRI software. But there is **another constraint** I did not mention in my first post. As the company is big, there is another "handicap" or problem: At this beginning stage I will have no control over the remote service, **another developers group will "decide" how to encode the data**. I mean, I think I cannot decide if the remote service is SOAP or not, REST or not and therefore, I don't know the XML tags they will use for the encoding. Maybe I can use the module you mention to connect to the service but **I am not sure how data will be served** on the other counterpart... I guess I will be told the general schema of the files, and then I have to try to connect ArcGIS to that remote service and prepare some scripts to parse the incoming XML and allow non-IT users to use it in their layers... Does it make sense? What do you think? It is a bit difficult to explain! :D Thanks for your response and ideas!


title: Are there any Geospatial OData Feeds?

text: I see that OData V3 has support for geospatial types. Has anyone published a feed that uses the new geospatial data types? ![enter image description here](http://i.stack.imgur.com/nIVW3.jpg)


title: Get polygons inside a polygon

text: I have developed a tool in arcobjects on Visual studio 2010, the tool works fine but once installed my ArcMap drawing the map (refreshing) and printing slows down.. I have used `ITopologicalOperator` in the program. Should I release any `com objects` here if yes what object should I release Please guide.


title: Displaying OSM file with markers

text: I'm trying to show OSM file which I exported from openstreetmap api. However, layer I add with OSM file is displaying like just polygons, I cannot show markers (street names, names of buildings etc.) in the map. Here is the map I'm seeing. ![showing like polygones](http://i.stack.imgur.com/uCXaF.png) And here is my openlayers init function.                       var lat=18.9003124; //changed due to company policies.             var lon=30.7738127; //changed due to company policies.             var zoom=15;             var map;                  function init(){                 map = new OpenLayers.Map ("map", {                 controls:[                     new OpenLayers.Control.Navigation(),                     new OpenLayers.Control.PanZoomBar(),                     new OpenLayers.Control.LayerSwitcher(),                     new OpenLayers.Control.Attribution()],                     maxExtent: new OpenLayers.Bounds(-20037508.34,-20037508.34,20037508.34,20037508.34),                     maxResolution: 156543.0399,                     numZoomLevels: 19,                     units: 'm',                     projection: new OpenLayers.Projection("EPSG:900913"),                     displayProjection: new OpenLayers.Projection("EPSG:4326")                 } );                           map.addLayer(new OpenLayers.Layer.OSM());                 map.layers[0].opacity = 0;                      var lonLat = new OpenLayers.LonLat(lon, lat).transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:900913"));                      map.setCenter (lonLat, zoom);                      //Initialise the vector layer using OpenLayers.Format.OSM                 var layer = new OpenLayers.Layer.Vector("Polygon", {                     strategies: [new OpenLayers.Strategy.Fixed()],                     protocol: new OpenLayers.Protocol.HTTP({                         url: "mymap.osm",   //<-- relative or absolute URL to your .osm file                         format: new OpenLayers.Format.OSM()                     }),                     projection: new OpenLayers.Projection("EPSG:4326")                 });                      map.addLayer(layer);      What am I doing wrong?


title: R: distance from one point to many polygons

text: OK so I have a spatial polygons data frame, and I have a point of interest. I want to add a row to the SPDF that contains the distance from the centroid of each polygon to this given point. How do I do this ? Here is a toy dataset:                   library(sp)         library(spdep)         library(rgeos)              # Create toy SPDF         fn <- system.file("etc/shapes/eire.shp", package = "spdep")[1]         prj <- CRS("+proj=utm +zone=30 +units=km")         eire <- readShapeSpatial(fn, ID = "names", proj4string = prj)              # Create a point of interest         pt.coords <-readWKT("POINT(260 5900)")         eire.proj<-CRS(proj4string(eire))         proj4string(pt.coords) <- eire.proj              # check point is on the map         plot(eire)         plot(pt.coords, col="Red", add=TRUE)      Ok so the SPDF data has following columns:                   > names(eire)        [1] "A"       "towns"   "pale"    "size"    "ROADACC" "OWNCONS" "POPCHG"  "RETSALE" "INCOME"  "names"        I want to add a extra column showing the distance from each polygon **_centroid_** to the point we defined. How do I do this ?


title: How to make marker line thicker in QGIS?

text: I'm mapping railways, find that the usual cross marker line is so thin it doesn't show up in Print Composer at web resolution, and can't find any way of making it thicker, only the main line. Is there a way in the just released 2.0? If not, how do I do it in 1.8? My System is QGIS 1.8.0 on OSX 10.6.8


title: How to mark a particular point in MapPanel?

text: There is a component called MapPanel, which displays OpenStreetMaps. I want to mark some places on the displayed map, without changing the underlying map. Those marked places are identified by latitude and longitude. In order to implement this, I wrote following class, which should convert latitude and longitude to pixel coordinates.               @Override     public Point convert(double aLatitude, double aLongitude, Rectangle bounds) {              final Point pointInMapUnits = mapPanel.computePosition(new Point.Double(aLatitude, aLongitude));              Point mapCenter = mapPanel.getCenterPosition();         Dimension mapSize = mapPanel.getMapSize();              final Rectangle panelSize = bounds;         final Point2D.Double panelCenter = new Point.Double(panelSize.x + panelSize.getWidth()/2,                 panelSize.y + panelSize.getHeight()/2);              final double ourObjectXFromMapCenterInMapUnits = pointInMapUnits.getX() - mapCenter.getX();         final double ourObjectYFromMapCenterInMapUnits = pointInMapUnits.getY() - mapCenter.getY();              final double ourObjectXFromMapCenterInPixels = convertToPixelsX(ourObjectXFromMapCenterInMapUnits, panelSize,                 mapSize, mapPanel.getZoom());         final double ourObjectYFromMapCenterInPixels = convertToPixelsY(ourObjectYFromMapCenterInMapUnits, panelSize,                 mapSize,mapPanel.getZoom());              final Point result = new Point((int)(panelCenter.getX() + ourObjectXFromMapCenterInPixels),                 (int)(panelCenter.getY() + ourObjectYFromMapCenterInPixels));              return result;     }      But it doesn't work as following screenshots illustrate. Step 1: Map is centered at the point in question (132.18381, 43.35514). That place is marked by a red dot. ![Image 1](http://i.stack.imgur.com/xmuNP.png) Step 2: Now we scroll the map so that the marked place is located to the left of the position in step 1. The blue dot shows where the red dot should be. ![Image 2](http://i.stack.imgur.com/77i0y.png) Step 3: Now we scroll the map so that the marked place is at the top. Again the blue dot shows the expected position of the red dot. ![Image 3](http://i.stack.imgur.com/NdiPe.png) Here's a video of the process. You can find the code here. How can I change it, in order for the dot to always be displayed at the point (132.18381, 43.35514) ?


title: Python: How to create a Polygon Shapefile from a list of X,Y coordinates?

text: I have a list of X,Y coordinates from a csv file, that represent a polygon. I´m trying to create a Polygon-Shapefile from this list. I´ve been trying around and found a possbility to write the list to a Point-Shapefile. Unfortunately that is not enough for me. Is there any way to get the coordinates in a Polygon-shapefile straight away? Following Brad´s suggestions I tried the following code:               for i in list:         w = shapefile.Writer(shapefile.POLYGON)         w.poly(parts=[list])         w.field('F_FLD','C','40')         w.field('S_FLD','C','40')         w.record('First','Polygon')         w.save('C:/Users/.../Desktop/Shape')      Unfortunately I´m getting an error message:               ShapefileException: Failed to write shapefile bounding box. Floats required.      Looks like there is a problem saving the shapefile to disk. This seems to be a similar question, but I couldn´t work out, how to get it going in my case.


title: Custom ArcGIS GraphicsLayer (JavaScript)

text: I am working on an application that will retrieve data from various sources and construct ESRI GraphicsLayer objects from the data and display it on a map. I have created custom FeatureLayers before but this project requires the use of GraphicsLayers because I need to be able to toggle the visibility of layers. The code below fetches the data from the host and puts it into a GraphicsLayer.               define(["dojo/_base/declare", "dojo/_base/array", "dojo/request", "esri/graphic", "esri/geometry/Geometry", "esri/InfoTemplate"],            function (declare, array, request, Graphic, Geometry, InfoTemplate) {         return declare(null, {           getAllCurrentReadings: function () {             var rtn = [];             var stations = ["S", "SN", "AN", "UP", "GR", "PL", "SR", "J", "N", "FL"];             array.forEach(stations, function (item, i) {               request.post("includes/buoybay_proxy.php", {                 data: {                   "method": "RetrieveCurrentReadings",                   "params": "CBIBS," + item + ",113f8b...f27e0a0bb" // NOTE: id: 1 is necessary as well but is added manually by jsonRPCClient                 },                 sync: true,                 handleAs: "json"               })                 .then(                   function (response) {                     var gfx, attr, t;                     //console.log(response);                     // Now build the Graphic Object and push it into rtn                     gfx = new Graphic();                     gfx.spatialReference = {                       wkid: 102100                     };                          // Define attribute object                     attr = {};                     attr["station"] = response.station;                     attr["title"] = translateStationID(response.station);                     for (var j = 0; j < response.measurement.length; j++) {                       attr[String(response.measurement[j])] = response.value[j];                     }                     gfx.attributes = attr;                          // Define geometry object                     gfx.geometry = new Geometry(gfx.spatialReference, "point");                     gfx.geometry.spatialReference = {                       wkid: 102100                     };                     gfx.geometry.type = "point";                     t = esri.geometry.geographicToWebMercator(new esri.geometry.Point(attr["longitude"], attr["latitude"], gfx.spatialReference));                     gfx.geometry.x = t.x;                     gfx.geometry.y = t.y;                          // Define infoTemplate object                     gfx.infoTemplate = new esri.InfoTemplate();                     gfx.infoTemplate.setTitle(attr["title"]);                     gfx.infoTemplate.setContent("${*}");                          // Define symbol                     gfx.symbol = new esri.symbol.PictureMarkerSymbol("../images/marker.png", 15, 15);                          //console.log(gfx);                     rtn.push(gfx);                   },                   function (error) {                     console.log("Error: " + error + "\n");                   }               )             });             //console.log(rtn);             return rtn;           }         });       }     );      This code seems to construct the GraphicsLayers properly but when I add them to the map object no points are displayed on the map. The code I use to add them to the map object is below.               require(["dojo/parser", "dojo/_base/array", "dijit/layout/BorderContainer", "dijit/layout/ContentPane", "dojo/ready", "esri/map", "esri/layers/ArcGISTiledMapServiceLayer", "js/cbibsGfxModule", "dojo/domReady!"],         function(parser, array, BorderContainer, ContentPane, ready, map, ArcGISTiledMapServiceLayer, cbibsGfxModule) {             var Map, cbibs, gfxLayer, t = [];             function init() {                 Map = new map("mapDiv", {                   basemap: "oceans",                   center: [-77.0357, 38.7877],                   zoom: 7                 });                 dojo.connect(Map, "onLoad", displayData); // Map didn't load until 3rd arg was a function name; why?                      function displayData() {                     cbibs = new cbibsGfxModule();                     t = cbibs.getAllCurrentReadings();                     gfxLayer = new esri.layers.GraphicsLayer();                     array.forEach(t, function (item) {                         gfxLayer.add(item);                         Map.graphics.add(item);                     });                     gfxLayer.spatialReference = {                         wkid: 102100                     };                     //Map.addLayer(gfxLayer); // Add GraphicsLayer to Map object                     console.log(Map); // Custom GraphicLayers are under _layers                 };             };             dojo.ready(init);         }     );      I realize that `gfxLayer.add(item)` and `Map.graphics.add(item)` are somewhat redundant but even with the data in two locations in the Map object the points are still not displayed on the map. I've been working on this for some time now and I really am fresh out of ideas. Any help anyone could offer would be greatly appreciated. Thank you.


title: Problem creating a popup in mobile openlayer

text: I'm trying to use the mobile version of openlayer, unfortunately I could not use framedcloud in it, I could only create a square popup does not look very aesthetic, besides the main need is to identify from a cloud of the associated data points so the popup square is not useful and less use jquery or another library.. The way we achieve using the square was import popup popup openlayers library, but the trick did not work with framedcloud (firebug says "OpenLayers.Popup.FramedCloud is not a constructor")               <link rel="stylesheet" href="../scripts/OpenLayers-2.13.1/theme/default/style.mobile.css" type="text/css" >     ....          <script type="text/javascript" language="javascript" src="../scripts/OpenLayers-2.13.1/OpenLayers.mobile.js"></script>     <script type="text/javascript" la language="javascript" src="../Scripts/OpenLayers-2.13.1/lib/OpenLayers/Popup.js"></script>          .....          function putPopup(x, y, html,map) {             // alert("popup1"+x+" "+y+" "+html)             map.setCenter(new OpenLayers.LonLat(x, y));                  if (popup) {                 map.removePopup(popup);             }             // this.popup = new OpenLayers.Popup("chicken",                     // new OpenLayers.LonLat(x, y), null, html, null, true);               //-------------------------Fail-------------------------------------                                  //popup = new OpenLayers.Popup("chicken",     //                                 new OpenLayers.LonLat(x, y),     //                                 new OpenLayers.Size(120,50),     //                                 html,     //                                 true);               //-------------------------Fail-------------------------------------                                                              popup = new OpenLayers.Popup.FramedCloud("Popup",                  new OpenLayers.LonLat(x, y), null,                 html, null,                 true // <-- true if we want a close (X) button, false otherwise                  );                        //-------------------------Fail-------------------------------------                                            //var myPopupClass = new OpenLayers.Class     //              (     //                  OpenLayers.Popup.FramedCloud,     //                  {     //                      'autoSize': true,     //                      'maxSize': new OpenLayers.Size(500,500)     //                  }            //              );     //                        // var myPopup = new myPopupClass     //                              (     //                                  "chicken",                                           //                                  new OpenLayers.LonLat(x, y),     //                                  null,     //                                  html,     //                                  null,     //                                  true     //                              );     //                              map.addPopup(myPopup);                                   map.addPopup(popup);              };


title: OpenLayers - Get LayerName from GML Response object

text: I am using the WMGetFeatureInfo control in OpenLayers to select items from my layers and return GML. I use this to get an attribute from the GML;                var val = evt.features[0].attributes.UID;     //Where evt is the object returned from "getfeatureinfo" event.      But the GML also contains a whole bunch of other good stuff in between the " wfs:FeatureCollection", and "gml:featureMembers" sections. How can I iterate through the wfs:FeatureCollection or gml:featureMembers because I cannot figure out the syntax. Also, would it be easier to parse the info I need from the XML that's returned? Thanks!


title: Sextante in Qgis 1.9 alpha on Ubuntu can't find saga gis 2.1

text: I am using Qgis 1.9 on Ubuntu 13.04, Sextante is a core part of qgis now. I am trying to access the SAGA 2.1 libraries through qgis and I keep getting a popup window that says SAGA is not installed. I can run SAGA and it works. Normally, I would go into the sextante configuration menu and make sure sextante is looking the right directory (/usr/lib/saga) However, that part of the configuration menu is gone. I went to the qgis options and checked the environment variables and I do not see one for saga, nor is it listed in the PATH variable. I tried adding in a SAGA path variable with no luck. Does anyone here know how to update or add a system variable to qgis 1.9 alpha so that qgis will know where saga is? The end goal will be the ability to have sextante properly access the saga libraries. Many of my Sextante modeller scripts rely HEAVILY on those libraries. I love qgis 1.9 and would hate to have to go back to 1.8 if i do not have to. Thanks EDIT: I was able to install the nightly version of 1.9 alpha on my windows work terminal using the OSGeo4w installer. I went into the advanced install options and made sure I installed SAGA (2.0.8.3) at the same time. Same problem seems to be happening. Everything seems to be working just fine when. Based on these findings I purged everything off of my Ubuntu laptop and tried again ( this time installing SAGA first). Still my problem persist. EDIT 2: I installed the official release of Qgis 2.0.1 from the Ubuntugis unstable PPA today. In this instillation some of the dependencies are dated and can not be used with SAGA gis 2.1. Here is the error message containing the offending package names. Some packages could not be installed. This may mean that you have requested an impossible situation or if you are using the unstable distribution that some required packages have not yet been created or been moved out of Incoming. The following information may help to resolve the situation:               The following packages have unmet dependencies:      saga : Depends: libgdal1 (>= 1.8.0) but it is not going to be installed             Depends: liblas1 but it is not going to be installed     E: Unable to correct problems, you have held broken packages.      The version of `Libgdal1` that gets installed with qgis is 1.7.4 ( something less then 1.8 ) when I try to update these dependencies, it removes half of the dependencies used by qgis effectively removing it.


title: Geocoding using Excel spreadsheet

text: Can I add a spreadsheet with addresses to Qgis and then geocode it?


title: Using multiple fields for a Reclass calculation in ArcGIS Field Calculator

text: I am trying to create a dummy variable from an existing dataset using ArcGIS's Field Calculator tool with a Python script. My conditions require the use of multiple fields to determine the output dummy variable. I have previously run this same script using a single variable condition without error. My code looks like this:               expression= "Reclass(!avg_incpsf!, !strip_dum!)"     codeblock= """def Reclass(avg_incpsf, strip_dum):        if (avg_incpsf >= 0 and avg_incpsf <= 26.1795309) and (strip_dum = 1):          return 1        elif (avg_incpsf >= 0 and avg_incpsf <= 26.1795309) and (strip_dum = 0):          return 2        elif (avg_incpsf > 26.1795309 and avg_incpsf <= 33.9938753) and (strip_dum = 1):          return 3        elif (avg_incpsf > 26.1795309 and avg_incpsf <= 33.9938753) and (strip_dum = 0):          return 4      etc... When I added the strip_dum variable condition I get back the following error: ExecuteError: Failed to execute. Parameters are not valid. ERROR 000989: Python syntax error: Parsing error SyntaxError: invalid syntax (line 2) Failed to execute (CalculateField). I can't seem to figure out how my parameters are invalid. I've also tried to no avail:               expression= "Reclass(!avg_incpsf!, !strip_dum!)"     codeblock= """def Reclass(avg_incpsf, strip_dum):        if (avg_incpsf >= 0 and avg_incpsf <= 26.1795309 and strip_dum = 1):          return 1        elif (avg_incpsf >= 0 and avg_incpsf <= 26.1795309 and strip_dum = 0):          return 2      Suggestions are much appreciated. thanks in advance, mike


title: How to add any WMS layer on Google Earth Plugin in my web application?

text: I have a GIS application which allows user to draw features on Google Earth plugin included on my web application. Now i want to give an option of adding any wms layer on the google earth plugin if needed. I searched on the internet for answers but the only solution was through **Image Overlay** option which is not available on google earth plugin.


title: How to export QGIS elevation data to Mysql WAMP

text: is there a way to export elevation data from QGIS to wamp server(mysql)? please help me. I need to store all the data to mysql wamp.. thanks!


title: Is there any help file to understand resource.xml of network Analyst?

text: I want to translate the direction result of network analyst to turkish. I surf in the interent and found i have to change resource.xml file . the file located in arcgis install directory (...\ArcGIS\Desktop10.0\NetworkAnalyst\Directions) .I want to find a help file or dictionary file to translate each tag and attribute of the resource.xml file.For example what does monobranch means in the xml file? Is there any help file to understand this file ? Thanks for advance


title: Adding a GRASS module manually in QGIS (Mac OS X)?

text: I am trying to add a precompiled non-standard GRASS module to QGIS 2.0.1 (MAC OS 10.9.1, using the binaries by Kyngchaos). Adapting the instructions here (http://hub.qgis.org/projects/quantum- gis/wiki/Adding_New_Tools_to_the_GRASS_Toolbox) I have copied the GRASS module to `/Applications/QGIS.app/Contents/MacOS/grass/bin` and the module description together with two thumbnails to `/Applications/QGIS.app/Contents/Resources/grass/module`. Finally, I have added a line with the module name to the file `default.qgc` residing in `/Applications/QGIS.app/Contents/Resources/grass/config/`. Still, however, the module does not appear in the list of available modules in the GRASS toolbox in QGIS (restarted, of course). Any hints towards a possible solution are very much appreciated.


title: SelectFeature and ModifyFeature in OpenLayers 2.13

text: After lot of experimenting, I have to give up and ask: has anyone had success in using both SelectFeature and ModifyFeature in OL 2.13 (NB, there were changes to ModifyFeature in 2.13, compared to 2.12)? What I want to do: ModifyFeature should work on it's regular way (user can drag points; and click on lines will bring "handles" visible). In addition, I want to catch select and hover events.


title: Map Info: Thematic Mapping Symbols part printing

text: I have produced a series of thematic maps for graduated housing figures, using a symbol which I customized, in three different tables displayed in a layout. This looks fine in the map window, although the symbols are sized much larger than in the layout and do overlap, in order to be a reasonably large size in the layout window. All looks fine in the layout window, no overlapping. However, when printed, a couple of the symbols have parts missing. One has the top half and the other the right hand half. The symbols for these tables are small dots and don't overlap. Any ideas please?


title: circle becomes an oval after saving

text: Using esri desktopv10 i am editing a "tank" layer.   It is selected as the only editable layer It has three features all of which have circle as the diagrammatic property. I use the circle tool and construct a new circle. It draws nicely.   Then I save edit Then I stop editing the layer. Visually it turns into an oval! Thoughts? ![circle after adding feature to layer](http://i.stack.imgur.com/kRbC7.jpg) ![the oval appears after saving and exit editing](http://i.stack.imgur.com/Zl7UW.jpg) Thanks to all for answering!


title: osm_id in a shapefile vs. the uid and id attributes in an OSM file

text: I'm parsing an **OSM file** using the Python package `osmread`. The entities parsed have the following attributes `id`, `version`, `changeset`, `timestamp`, `uid`, and `tags`. On the other hand, I have a **shapefile** , and I use the Python package `pyshp` to read it. I then get the following fields from this (road) shapefile:                   [('DeletionFlag', 'C', 1, 0),         ['osm_id', 'N', 11, 0],         ['name', 'C', 48, 0],         ['ref', 'C', 16, 0],         ['type', 'C', 16, 0],         ['oneway', 'N', 1, 0],         ['bridge', 'N', 1, 0],         ['tunnel', 'N', 1, 0],         ['maxspeed', 'N', 3, 0]]      My question is, if the road shapefile was obtained from the OSM file, which of the two attributes, `id` or `uid`, in the OSM file corresponds to the `osm_id` field in the shapefile?


title: Why mobile routing is so fast?

text: Why my mobile GPS routing calculates routes almost instantly and using pgRouting in my robust PC is so time costly ? **EDIT:** My test: In my cellphone (Android 4.2.2 - 1G RAM - Dual core - Arm v9) I have iGo Primo. In my PC ( 8G RAM - Core i5 ) I have Postgre + PostGIS + pgRouting Table OSM Planet "planet_osm_line" with previous work:               alter table planet_osm_line add column source integer;     alter table planet_osm_line add column target integer;     select pgr_createTopology('planet_osm_line', 0.0001, 'way', 'osm_id');          select * from planet_osm_line where name like '%BR-040%' // (source : 2751784)     select * from planet_osm_line where name like '%BR-101%' // (target : 2040999)          SELECT seq, id1 AS node, id2 AS edge, cost, name, way       FROM pgr_dijkstra('SELECT osm_id as id, source, target, st_length(way) as cost FROM planet_osm_line',  2751784, 2040999, false, false  ) as di        JOIN planet_osm_line pt       ON di.id2 = pt.osm_id ;      The route is (both Brazil Highway): Source: BR-040 Target: BR-101 My cellphone do the job in less than two seconds and my PC takes more than 3 minutes (only the SELECT route - pgr_dijkstra - part).


title: Joining data points to matching polygons

text: I have polygon data on coral reefs at great barrier reef and separate point data for areas of coral bleaching. I have created a clip that only includes points that exist with in the reef polygons. However I need to be able to identify polygons that have bleaching points in them from those that do not. So for example I would like to be able to click on a polygon that has a bleaching point in it and it automatically select all the other polygons that do?


title: Accurate local map projections

text: Is there a map projection that can accurately represent reality at a local area, i.e. when transforming the map could I choose an area that I care about accuracy and allow distortions elsewhere? Like this map: http://upload.wikimedia.org/wikipedia/commons/thumb/8/83/Tissot_indicatrix_world_map_equirectangular_proj.svg/450px- Tissot_indicatrix_world_map_equirectangular_proj.svg.png I want to choose where the distortion happens.


title: C# - Access content from ArcGIS Online

text: I have got an existing .NET C# desktop application which already reads various GIS formats. I would now like to add support to read data from ArcGIS Online. What options have I got to access information such as services, groups, search data etc from ArcGIS Online? Can I use ArcObjects, REST endpoints?


title: Openlayers OverviewMap control with vector layer data display

text: I am creating an OverviewMap control in the map and want to see the content of a Vector layer in the overview map too.                       var mapOverviewOptions = {                 mapOptions: {                     projection: new OpenLayers.Projection("EPSG:900913"),                     displayProjection: new OpenLayers.Projection("EPSG:4326"),                 },                 layers: [baseOverviewLayer, myVectorLayer]                  };             var ovMap = new OpenLayers.Control.OverviewMap(mapOverviewOptions);      This does not work. It only displays the world map background and does not display the vector data that gets rendered on the map. Is it possible to display the vector layer data in the overview map? I tried adding the vector layer to the map object that Overview Control creates,               ovMap.ovmap.addLayer(myVectorLayer)      but that does not help either. Any help would be greatly appreciated.


title: How to avoid setting env.workspace and other paths/settings all the time

text: I have some scripts that are constantly setting env.workspace of arcpy. Usually, I have a caller script( **script1.py** ) that may call some other script ( **script2.py)**. I've created a file **globals.py** with the paths that never change:               import arcpy.env as env          rootPath="S:/workspace/MSc_git/ArcGis/"     layersPath=rootPath+'Layers/'     # Set the workspace environment to local file geodatabase     env.workspace=rootPath+'estagio_db.gdb'      In script1.py I have something like this:               import arcpy     import globals as g     from arcpy import da     arcpy.AddMessage("Importing arcpy outside ArcMap took %s seconds." % (time.clock() - t_start))     from cProfile import run               def main():       ....      And I need to use env.workspace and sometimes the layerspath inscript2.py How can I do this (should I import gloabl.py as well?)? should I pass it by argument when calling a function from script2.py in script1.py or is there any other way to accomplish this? I know we should avoid global vars, but in this case, and since I need to centralize its definition I don't know how to do it well...


title: ArcGIS - Help Using Oracle Function in Definition Query - No ArcSDE

text: I have set up a Query Layer pointing at an Oracle table with a SDO.Geometry column, and that is working and behaving fine. However, when I am attempting to set up a definition query that calls a Oracle function, I get an error message saying "There was an error with the expression. Failed to parse the where clause". I am fairly confident the SQL expression I am using is valid: it works fine inside of SqlDeveloper. Here is what the sql looks like:               select * from gis_meters where gis_getmeterstatus(METERNAME,'2013/08/08 10:25:00','2013/08/09 10:25:00') = 0      In the Query Builder of ArcMap, I am only adding, as the query layer alread has the select *:               gis_getmeterstatus(METERNAME,'2013/08/08 10:25:00','2013/08/09 10:25:00') = 0      After some research, I found some ESRI documentation stating "Coverages, shapefiles, and other nongeodatabase file-based data sources do not support subqueries". I am guessing that the call to the Oracle function (gis_getmeterstatus) is a subquery - and thats why its not working for me? Since I don't have SDE installed in the database, its not considered a geodatabase. Can anyone help me find a way to use this function in my definition query?


title: qgis --- style edges of layer shapefiles to have thicker borders

text: I have two layers I'm working with here, both of North Carolina. One is "townships" the other is "counties". I want the counties shapes to have thick edges with no fill-in color. I couldn't find anything on this after a google search so decided to ask. Thanks


title: How to export aerial imagery as layers in layered PDF?

text: I want to show aerial imagery as layers in my pdf, so people can turn them on and off (e.g historical images). However according to ESRI images become merged when exported as a layered pdf. What is the point of this? The only way around it would be to merge the tif (or ecw in this case) to a single band image but this would result in it being grayscale I think? I tried everything...


title: Setting proxy url for larger request

text: When I query for larger features I get `esri.config.defaults.io.proxyUrl is not set` error. Can someone please explain why proxyurl is needed, how I should set it up, what should it point to, and any other relevant information about it? I am running Websphere and IBM http server.


title: ArcSDE - Is there an automated way of releasing all layers locks?

text: We have an issue where users do not close ArcGIS and layers get locked, and sometimes this affects our nightly updates for some of the layers. I know you can individually remove a lock using the ArcSDE command (which we do), but does anyone know of a way I can create a script that can be run every night to get rid of the locks? We are running ArcSDE 9.2 and Desktop 9.3.1. Thanks


title: Consuming ArcGIS Server map services from WCF

text: I am consuming some map services in my WCF services and for each service that I have I get the client object via this code. Every new service that I add, I get the client object this way. I am curious if I can consume the map service using it's URL so that I dont have to add the map service as a service reference. I would like to use the URL of map service to create client object. Appreciate any feedback. MapServerPortClient client = new WellMapService.MapServerPortClient(); OR MapServerPortClient client = new PipelineMapService.MapServerPortClient();


title: Accessing GPS data with JavaScript in real time

text: Is it possible to access GPS data in real time on a smartphone with JavaScript (without an android or ios application I mean)? Thanks


title: How can i set layers name from map to combobox?

text: I'd like to list map layers from map into combobox. I tried to set itemsource of combobox to: `ItemsSource = "{Binding ElementName=Map, Path=Layers}`. It does not list the layer names. It just lists "Esri.ArcGIS.Client.FeatureLayer". How can I list layers with their names in my combobox?


title: pgRouting 2.0 slow execution of A star algorithm

text: Recently, I wanted to test pgRouting ver. 2.0. I created the database structure with osm2po and imported UK map from http://download.geofabrik.de/europe.html. The data was imported correctly. When I execute a sample query in order to find the shortest path with dijkstra algorithm I receive the result in ~70 sec. When I do the same with astar algorigthm I get the result in over 145 sec. Theoritically, the result of astar should be faster than dijkstra's one but it's not. Queries I execute: --Dijkstra `SELECT seq, id1 AS node, id2 AS edge, cost FROM pgr_dijkstra('SELECT id::integer, source::integer, target::integer, cost::double precision AS cost FROM hh_2po_4pgr'::text, 60, 80, false, false);` --Astar `SELECT seq, id1 AS node, id2 AS edge, cost FROM pgr_astar('SELECT id AS id, source::integer, target::integer, cost::double precision AS cost, x1, y1, x2, y2 FROM hh_2po_4pgr', 60, 80, false, false);` I checked the Astar implementation in pgRouting and it seems that it uses this heuristic: `(|dx|+|dy|)/2`. Am I doing something wrong in my queries or there is something wrong with pgr_astar() implementation in pgRouting?


title: Is it possible to have an infowindow pre-opened in when you load a Cartodb.js map?

text: Can I code my Cartodb.js map somehow so that when someone first opens the map, there is already an infowindow open for them to see?


title: Interpretation of PostGIS geocode function ratings

text: Edit: I am using Postgres 9.2 and recently installed PostGIS along with the Tiger Geocoder from the extras folder and downloaded all the latest Tiger data. I have been testing the Tiger Geocoder function (http://postgis.net/docs/Extras.html) to process a big database of addresses in order to normalize/validate them and get lat-long coordinates. I know that the function returns a rating between 0-100 (0 being the best), but I have not been able to find any documentation stating how to interpret the score. How can I interpret the scores between 0 and 100?


title: Web mapping service with GEBCO bathymetry data

text: Is there any open web mapping service provide map tiles with GEBCO bathymetry data. As far I know Only Google uses GEBCO bathymetrydata for ocean floor model. Others uses ocean as plain blue. I need a service provider of web mapping tiles which provide sea floor with GEBCO data. Thanks


title: Embeded video won't load in Cartodb.js

text: I have custom infowindows in my map that have embed youtube videos. When I'm working directly in cartodb online, everything loads perfectly. Within the content section of my infowindow code, I have the following: < iframe width="640" height="390" src="{{embed_link}}" frameborder="0" allowfullscreen> < /iframe> I am trying to create a similar map using Cartodb.js (using the api from the Cartodb map I have created with the infowindows that are able to load the video), and the videos were loading fine for a while, but now they have stopped loading. Where the embeded video should be, there's just a blank grey box that says: "The file or directory could not be found." Is there any reason why the videos might just stop loading when I create a map using Cartodb.js? Two examples of the embed_link are: "//www.youtube.com/embed/HyCBk9XC2Eo?feature=player_detailpage" and "//www.youtube.com/embed/ZZAFYROf2mY?feature=player_detailpage"


title: Arcpy MosaicToNewRaster problem

text: I can't seem to get the mosaic created for this.                import arcpy           arcpy.env.workspace = r"C:\projectTemp\FloodAnalysis\output\2yr\raster"      arcpy.env.overwriteOutput = True      srName = r"C:\projectTemp\FloodAnalysis\output\2yr\raster"      outPut = r"C:\projectTemp\FloodAnalysis\output\2yr"      rasterName = "2ft.img"      bands = 1      method = "First"      clMap = "First"           arcpy.MosaicToNewRaster_management(srName, outPut, rasterName, "", "", "", bands,      "First", "First")      Thanks for the help.


title: need help with L.mapbox.legendcontrol

text: I am using mapbox.js 1.6.1 to serve .mbtiles created in tilemill (via mbtiles- server.php) on a private host. my knowledge of html and js is very (very) basic, I can't figure out how to use L.mapbox.legendControl properly to add a legend to a L.mapbox.map map. I am using mbtiles-server.php to serve .mbtiles created in tilemill. Here is what works:               //create map, centered on papua       var map = L.mapbox.map('map', 'server/mal.tilejson').setView([-6, 140], 6);      this creates a map with the tile- and gridlayer (hover tooltips) as well as the legend. I have not been successful with creating an empty map and adding layers with their legends showing:               //create map, centered on papua       var map = L.mapbox.map('map').setView([-6, 140], 6);          //baselayer       var srtm = L.mapbox.tileLayer('server/neuguinea_zl3-14.tilejson', {                             maxZoom: 14,                     }).addTo(map);          //overlay      var malgroup = new L.LayerGroup();      var malariatile = L.mapbox.tileLayer('server/mal.tilejson').addTo(malgroup);      var malariagrid = L.mapbox.gridLayer('server/mal.tilejson').addTo(malgroup);       map.addControl(L.mapbox.gridControl(malariagrid));          //create baselayers and overlays     var baseLayers = {       "SRTM v3": srtm       };     var overlay = {       "Malaria": malgroup       };          //add layer control      L.control.layers(baseLayers, overlay).addTo(map);          // create legend       map.addControl(L.mapbox.legendControl());      I was hoping that map.addControl(L.mapbox.legendControl()); would activate the overlay's legend automatically, but it doesn't. Thanks for looking into this!


title: Why can't I reproject from HARN State Plane Ft Intl to State Plane Ft?

text: I'm attempting to reproject a feature class that is:               NAD_1983_HARN_StatePlane_Oregon_North_FIPS_3601_Feet_Intl       to our standard, which is:               NAD_1983_StatePlane_Oregon_North_FIPS_3601_Feet      When I use the Project tool it automatically selects this transformation:               NAD_1983_To_HARN_WA_OR      But when it finishes, the two feature classes are not aligned propery. ![Misalignment](http://i.stack.imgur.com/kdhAs.png) Now there are a gob of other transformations, but how do I really know which one to choose? I've always been a little confused by this. ![Transformations](http://i.stack.imgur.com/IamcB.png) Is the 3 ft misalignment to be expected? Thanks for any guidance. * * * EDIT I tried this transformation:               NAD_1983_HARN_To_WGS_1984_2 + WGS_1984_(ITRF00)_To_NAD_1983      And now the misalignment is well under a foot. But is this acceptable? To just choose a transformation that works best? And is it a true statement that after the reprojection there may be no transformation that aligns the reprojected data perfectly? ![Closer](http://i.stack.imgur.com/AypGj.png)


title: ArcGIS labels with expression not appearing

text: I have a layer which should be labelled based on two different fields. I've chang the Expression in Layer Properties from **[FIELD1]** to **[FIELD1] + '\n' + [FIELD2]**. However, only Field1 is showing in my map, never Field2. I have fixed the scale range and I am using the Maplex Label Engine. Why aren't both fields showing in my label?


title: Automatically insert file path in QGIS table

text: I'm using R to create some plots that refer to specific values of some point in QGIS. Each plot is saved into a directory and for each point in QGIS I added a column with the path to the file. I know 2 different method to display these plots in QGIS:   * using the Action tool (Open file)   * using the Field tab of the Properties menu of the layer changing the Widget from Modify value to Photo What if I have to do tons of plots or if I have a large point dataset? Is there a way to add automatically the path to the file in the attribute table? I tried to use the field calculator but I cannot figure out the problem. Thanks to all BTW: QGIS 2.1.0 master on Debian 7.3


title: How to copy selected features from one sde to another using ArcObjects?

text: Assuming I have two SDEs or other datasources with the same datamodel. One contains features the other one is empty. In ArcMap I have layers from both datasources in my table of contents. Now I would like to copy some features that are already selected in the map view from one sde to the other using C# and ArcObjects. What would be the best way to do this?


title: What would cause Map Tips to flash in ArcMap with Shapefile?

text: Using ArcGIS for Desktop Advanced 10.1 SP 1 I have come across this several times but I have never found a solution for it. The problem I am having is that when I set Map Tips to display for certain shapefiles, the map tips flash on the screen. I've made a YouTube video to demonstrate this problem. There is nothing particularly special about this shapefile... it has about 1800 features from a Geocoding result. As long as I have this shapefile turned on, it causes all of the Map Tips to flash. If I turn this shapefile off, then map tips for other layers are ok (demonstrated in video...) I've tried rebuilding the spatial index, that does not fix it. If I convert the "problem" shapefile to a feature class in a file geodatabase, I no longer have the problem. What in a shapefile might cause the map tips to flash?


title: Crash when editing a layer and when saving in QGIS

text: When I edit the layers, and its small changes like adding a name (or make it appear in the map), QGIS very often crashes (like for every third or fourth or fifth edit, seldom gets past ten edits before crashing. But it can also be the first small change or edit, just afther opening QGIS and the project). Now I use QGIS Dufour 2.0.1 but the same happened with QGIS Lisboa 1.8.0. It also can crash when I try to save the project Im working on. From this question "Why does writing large shapefiles crash QGIS?" it may seeme like its due to filling up the memory? Is there any way to fix this, or get around it? If it's a common problem it must be a major bug or error that need to be fixed. When it happens QGIS seems to work endlesly while standing still, it sometimes comes a (not responing), or simular message in the window titlebar after the file name. If you then try to do anyting like open a menu or wathever, the whole window turn whiteish (or gray if you like), and waiting to see if it comes back dont seeme to work.


title: How to convert data in NAD83 (NSRS 2007) to NAD83 (CORS 96)?

text: I'm quite a ways outside my box on this one, but I would ceratinly appreciate any perspectives on this issue. I have some control points from NGS that are in NAD83 (NSRS 2007) that I would like to merge with the rest of our control point dataset, but our existing data are in NAD83(CORS 96) ... I think ... or it could be NAD83 (2001). Perhaps I should ask: Which NAD83 is inherent in ESRI's NAD_1983_StatePlane_Georgia_West_FIPS_1002_Feet coordinate system (Geographic: GCS_North_American_1983 Datum: D_North_American_1983)? In any case, does anyone know of any transformations for converting NAD83 (NSRS 2007) to one of the "regular" NAD83 datums? NGS will not be publishing any, and ESRI does not have a transformation filter for it. Any help would be appreciated. BTW, I did read the GIS-SE posts from last year regarding NAD83<-->HARN conversions, but I am wondering if anyone has dealt successfully with this issue since then.


title: Visualisation using simple plotting routines or advanced GIS packages?

text: When I have geographic data to be displayed, I usually use a plotting library to visualise the data, for example, combining `matplotlib` and `basemap`, or using the Matlab Mapping toolbox. For example, I have data gridded on a _lat_ - _lon_ grid, and to get them on a map, all I need is a projection and a bit of background map data (usually only coastlines). I am quite satisfied with the results of my visualisation. What does a full-fledged _GIS_ package add that I'm currently missing?


title: Which Windows GIS editors support the creation of 3D shapefiles besides ArcGIS with 3D Analyst?

text: I need to find and **easy and affordable** way to generate 3D linear shapefiles in Windows so one of our users can inspect our tool. We incorrectly assumed that he had 3D shapefile source data available or had a way to produce it. I found an ArcGIS Technical Article on how to convert a 2D shapefile to a 3D shapefile; however, it requires ArcGIS with 3D Analyst. I have ArcGIS available but do not have a 3D Analyst license, so I cannot replicate the steps from the article. The user also has ArcGIS available but I do not know if he has a 3D Analyst license. Regardless, I would like to test the process before sending him the article. I have previously created 3D shapefiles using the OGR/GDAL C++ API and the OGR/GDAL Python API included in FWTools and OSGeo4W. I considered writing a simple program or Python script that can parse an input text file. The Python solution would require the user to install a Python runtime with GDAL support, so I am hesitant to take that approach. MapWindow 4.7.5 has a _3D vector to 2D vector_ option in the _GIS Tools | Vector_ menu, but not the other way around. I just checked and MapWindow is at version 4.8.3, but could not find any reference to 3D shapefiles or vector features in its online documentation. **Is there a Windows GIS editor (besides ArcGIS with 3D Analyst) that supports the creation of 3D shapefiles?**


title: Advice on which version of AutoCAD to purchase

text: I am starting a consulting agengy and I need to determine what my best interests are regarding the purchase of AutoCAD software. I have great faith in the more reputable users on this site, and I'm seeking advisement; particularly from those of you who use **both** ESRI software and map capable versions of CAD. **A little background:** Due to variuos positions I've held in my locale, I have strong expertise in relation to Mining Land Tunure, and the Aggregate Resources Act (in Ontario, Canada). For those of you who are unfamiliar, Aggregate material refers to Pit & Quarry operations (sand, gravel & crushed stone). Also notewothy, there are a number of Environmental studies associated with these strands for which I provide AutoCAD, GIS, and legislative support. I leverage GIS and AutoCAD technology interchangably. I am extremely comfortable going back-and-forth between these two software packages, which allows me to bring a variety of effiiciencies to my worflows. Lastly, I am completely self-taught in AutoCAD, but I have exceptional strength with regard to using it for **my needs**. **My dilemma** I moved from LDD to Civil 3D a few years ago, and I simply can't imagine going back to a less capable version od AutoCAD now that I've experienced the Cadillac of the AotoCAD realm. The difference between map 3D and Civil 3D is $5000 +tax per seat. This is a lot of money at start-up. I only need a few tools (to my knowledge) from the Civil 3D environment (surface from survey, from .shp, etc). I don't need most of its functionality (corridors, alignments, etc, etc ................>. The reseller has indicated that he can give me a few of the Civil tools as part of the Map 3D package, but I have never used it "stand alone". Otherwise, I need spatial data in CAD for site plans via _mapconnect and mapimport_. For the most part, it will be shapefiles. If anyone could advise (who has knowledge of both Civil & map versions) it would be greatly appreciated. I have done more reading than my schedule can afford. Will Map 3D meet my needs? Are there marked differences and shortcomings when measuring Map 3D against the superior product. I appreciate any advisement.


title: Extending ArcGIS

text: I've been tasked with a thesis project where i have to extend the features of ArcGis. I've been asked to create a model written in Python that can run out of ArcGIS 10. This model will have a simple user interface where the user can drag/drop a variety of shapefiles and enter the values for particular variables in order for the model to run effectively. Once the model has finished running, a new shapefile is created that lays out the most cost effective Collector Cable route for a wind turbine from point A to point B. I'd like to know if such a functionality/ extension already exists in ArcGIS so i don't have to re-invent the wheel. If not then what is the best programming language to learn to extend ArcGIS for this (Python vs Visual basic vs Java). My background is Java, PHP, Jquery and Javascript. Also any pointers in the right direction i.e documentation, resources etc would be hugely appreciated


title: How to display coordinates under cursor in web map?

text: I'd like to be able to display the X/Y coordinates of the user's current mouse cursor location such as with this web map. I'm using the .NET Web ADF, C# and ArcGIS Server. Can anyone suggest any code samples for this?


title: Tutorial on how to model flooding using data from gauging stations along rivers in QGIS?

text: Can someone please direct me to a tutorial where I can learn more about how to model flooding using data from gauging stations along rivers using QGIS?


title: How to determine in which grid cell each coordinate of an x,y series is located?

text: I have a circular arena which I have divided into a number of equal-sized 'cells' using a grid-like arrangement of concentric circles and line segments (in grey below) and also a series of x,y coordinates (plotted as a line in black below). ![sampleFig](http://i.stack.imgur.com/RXKiG.png) I was wondering what would be the best method (preferably using R) to find out which 'cell' each point is in ? I guess I could test to see which circle each point lies in by its distance from the centrepoint, and then which sector the point lies in by its angle from the centrepoint, but this seems rather laborious and inelegant.. This question bears some similarity to this one, except that I'm using circular lines to demarcate grid cells rather than linear lines. Thanks!


title: How to install a patch in PostGIS

text: Can anybody tell me how to install this patch in PostGIS? http://trac.osgeo.org/postgis/ticket/413 http://trac.osgeo.org/postgis/changeset/7534 I need to use the ST_OffsetCurve function. but it seems it's not available in the PostGIS installation(I installed the version "postgis- pg90-setup-1.5.3-2"). Can anybody tell me how to make it work?


title: Driving Zone Computation using Web Processing Services

text: I am currently evaluating different OGC Services(WFS,WMS and WPS) for driving zone computation. As I am new to WPS services, I would like some suggestions regarding which OGC services to implement. I am thinking through these issues right now: If I use WPS to calculate driving zones, is it better to integrate this computation as Web Processing Service or to directly compute it in the App Server? Would the client be overwhelmed or exhausted with delays in computation? Is there any way to deal with driving zone computations that exceed some predefined computation time limit, like couple of seconds? Also is there any issue if my background map layer is in different road network than the road network where driving zone computation is being performed? Please correct me if my concept is misleading. Thanks


title: Silence user warnings for GRASS-PYTHON program

text: I've created a program that runs externally from GRASS v6.4.1, using python 2.6. It imports an ESRI shapefile, and exports a seperate KML for every polygon in this vector file. Polygons with islands produce a warning message using "v.out.ogr" saying that features without a category were skipped. Looking into this further, these "skipped" features (islands) are processed later. So the program does what it was designed to do, but I would like it so that this warning message is not printed in the terminal. VERBOSITY setting were set to "0" using:               os.environ['GRASS_VERBOSE']='0' #runs modules silently      But this doesn't take care of warning messages. I found that there is a gisenv variable called DEBUG, can be set to "0" as well so that all warning messages are silenced too. Problem is, I can't seem to get the proper syntax for my GRASS-python script. I've tried:               grass.gisenv() ['set']='DEBUG=0'      Using this line of code atleast allows the program to run, but warning messages are still printed. Other options I've tried terminate the program. Has anyone successfully set the DEBUG level to "0" in a GRASS-Python script? Cheers.


title: Maintain Feature Data without Creating It

text: I work in City Gov. We utilize County Assessor data (Parcels) as they are updated from the Assessors office. Currently we store then in our Enterprise SDE system. As we get new updates, we have to shut down all users of Parcel data in order to delete the old set and overwrite with the new update. Is there a way to update this data without having to shut down users? Some kind of versioning concept for maintainers vs. creators of gis data? We are currently looking into TRUNCATE AND APPEND as an option.


title: Qgis web client search

text: I install Qgis web client on windows 7 and it's working fine! Now I wont to setup search feature!   1. First I added line in Apache config file:   -LoadModule wsgi_module modules/mod_wsgi.so    2. In GlobalOption.js I set var useGeoNamesSearchBox = false; var searchBoxQueryURL   ="http://localhost:89/qgis-web-client-master/wsgi/search.wsgi?query="; var searchBoxGetGeomURL   ="http://localhost:89/qgis-web-client-master/wsgi/getSearchGeom.wsgi"; When i try to search, i do not get anything! Is there any example for searching setup?


title: How to write a python script to give multiple buffers of different sizes in ArcGIS?

text: I want to write a loop that will cause the output to be multiple buffer shapefiles. For example I want a circle buffer of 50ft, 100ft, 150ft..... Can anyone help? I'm new to this and am trying to figure it out myself.


title: Where can I find free Archeology datasets?

text: Are there any sites that gather together free datasets? I am especially interested in published data relating to archeology.


title: Data attributes in Illinois legislative and senate TIGER shapefile

text: > **Possible Duplicate:**   >  Illinois Legislative Zip to District Match Data Set I have to download illinois legislative and Senate shapefile from the census bureau. Do I download State Legislative district from the layer type? Will this file have senator, legislatures name, address and district?


title: layer mask does not mask labels, only vector data

text: I am using ArcGIS 10. I made a layer mask (polygon) that blocks out everyting in a number of vector layers that is outside the area of interest (outside the boundary of two districts). I successfully applied the mask, however I note that the labels outside the districts of interest remain (i.e. the labels inside the mask area). Some of the labels are in an annotation feature class in a geodatabase. Some are not (i.e. I just have the lable features checkbox ticked). All labels are not masked. I have even tried switching it from the Maplex labelling engine to the standard engine, to no avail.


title: Trouble with r.out.gdal in GRASS GIS

text: I'm processing some landsat 8 files. I'm using this scene LC82260982013325LGN00 and I have to apply i.landsat.toar and export the result as a GTiff. However when I use r.out.gdal it takes a very long time (more than 20 minutes per band) and the resulting file is blank and various GB in size, even when I use compression. I tried it GRASS GIS 6.4.4 and 7 beta 2. I have Win8 and 8 gb of ram. What I am doing wrong? **UPDATE:** The region seems to be well set up               g.region -p                                                                          projection: 1 (UTM)     zone:       19     datum:      wgs84     ellipsoid:  wgs84     north:      -53.95464176     south:      -6163823.95464176     west:       -67.95813741     east:       715222.04186259     nsres:      30     ewres:      30     rows:       205459     cols:       23843     cells:      4898758937          r.info map=TOAR_LC82260982013325LGN00_B1@Luciano                                      +----------------------------------------------------------------------------+      | Layer:    TOAR_LC82260982013325LGN00_B1  Date: Fri Jul 25 12:06:01 2014    |      | Mapset:   Luciano                        Login of Creator: Luciano         |      | Location: UTMZone19                                                        |      | DataBase: C:\Users\Luciano\Documents\grassdata                             |      | Title:     ( TOAR_LC82260982013325LGN00_B1 )                               |      | Timestamp: none                                                            |      |----------------------------------------------------------------------------|      |                                                                            |      |   Type of Map:  raster               Number of Categories: -0              |      |   Data Type:    DCELL                                                      |      |   Data Units:   unitless             Vertical datum: (none)                |      |   Rows:         8171                                                       |      |   Columns:      8181                                                       |      |   Total Cells:  66846951                                                   |      |        Projection: UTM (zone 19)                                           |      |            N:   -5918685    S:   -6163815   Res:    30                     |      |            E:     715215    W:     469785   Res:    30                     |      |   Range of data:    min = 0  max = 2.70260743795714                        |      |                                                                            |      |   Data Description:                                                        |      |    generated by i.landsat.toar                                             |      |                                                                            |      |   Comments:                                                                |      |     Reflectance of Landsat-8 OLI/TIRS (method dos4)                        |      |    -----------------------------------------------------------------       |      |     Acquisition date (and time) ........... 2013-11-21 (0.0000 h)          |      |     Production date ....................... 2014-02-17                     |      |                                                                            |      |     Earth-sun distance (d) ................ -0.0000000                     |      |     Sun elevation (and azimuth) ........... -0.00000 (0.00000)             |      |     Digital number (DN) range ............. 0 to 0                         |      |     Calibration constants (Lmin to Lmax) .. -0.00000 to +0.00000           |      |     DN to Radiance (gain and bias) ........ -0.00000 and +0.00000          |      |     Mean solar irradiance (ESUN) .......... -0.000                         |      |     Radiance to Reflectance (divide by) ... -0.00000                       |      |                                                                            |      |     Dark object (1000 pixels) DN = ........ 8739                           |      |     Mode in reflectance histogram ......... -0.00000                       |      |    ------------------------------------------------------------------      |      |                                                                            |      |    i.landsat.toar input_prefix="LC82260982013325LGN00_B" output_prefix=\   |      |    "TOAR_LC82260982013325LGN00_B" metfile="C:\Users\Luciano\Documents\G\   |      |    IS DataBase\Landsat\Landsat 8\LC82260982013325LGN00\LC82260982013325\   |      |    LGN00_MTL.txt" method="dos4" percent=0.01 pixel=1000 rayleigh=0.0       |      |                                                                            |      +----------------------------------------------------------------------------+


title: Shrinking a line to fit within a raster surface?

text: I have a single polyline feature that extends beyond my raster surface. I would like to shrink this line so that the end-points are within the last pixel at the edges of the raster. I have tried converting the raster to a polygon and clipping but it sometimes puts the endpoint just outside of the raster in no data land. Is there a simple method that I could use to intersect/shrink the line? ArcObjects is fine.


title: ArcPy Multiprocessing 000714 - Working in Debug Mode

text: In an attempt to render an interactive window via matplotlib I am using an out of process script in a toolbox. If I debug the script, it is possible to step through and render the window as expected. Running it through ArcMap (not in debug mode) causes a **000714 error**.   1. Arc10.2   2. Running the script out of process   3. Tested with background processing on and off (Geoprocessing Options)   4. pythonExe is being properly found Any thoughts?               import multiprocessing     import sys, os     import arcpy          import matplotlib.pyplot as plt          def show_form():         fig = plt.figure()         x = range(10)         y = range(10)         plt.plot(x,y)         plt.show()          def main():         pythonExe = os.path.join(sys.exec_prefix, 'python.exe')         multiprocessing.set_executable(pythonExe)         multiprocessing.freeze_support = True         jobs = []         pool = multiprocessing.Pool(1)         jobs.append(pool.apply_async(show_form,[]))         pool.close()         pool.join()         del pool         del jobs         arcpy.SetParameterAsText(0, True)                 if __name__ == "__main__":         main()      ## Update 1 As per Luke's comment, I tested `cmd.exe /C test.py`. Other than the long import time for arcpy the script executed as expected with the Tk plot window appearing. ## Update 2 By opening a command window and quickly snapping a screen capture, I can grab an error from the script before it dies. This looks like an issue loading the libraries (32 vs. 64bit?) ![enter image description here](http://i.stack.imgur.com/CF7I9.png) ## Update 3 Arcmap is using the 64bit python by default when running out of process. The next step is to figure out how to disable this, as it appears to be causing a Tk issue as @Luke suggested.


title: Raster dataset doesn't exist error

text: Dataset exist and I am able to get statistics from ArcGIS python window. But it failed as a script with error: ERROR 000865: Input raster: C:\temp\ras does not exist. I have tried following suggestions: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//002z0000000r000000 but to no avail.Please suggest how to handle this               import arcpy          stat = arcpy.GetRasterProperties_management(r'C:\temp\ras', "MEAN")     avg = stat.getOutput(0)     print avg


title: How is GIS used in landscape architecture?

text: I was wandering, if out there somewhere are people that use GIS extensively in the field of landscape architecture (or closely related) and how do they use it? May be some examples of work, maps, visualisations? I am landscape architect, but at the time I am using GIS mainly in the field of hydrology. In the past I did some visibility analyses with use of MapInfo and Vertical Mapper extension and also some landscape character analyses/mapping. Today I am working mainly with ESRI and OpenSource and I do know that capabilities of GIS are much more than visibility analyses. In example in landscape ecology (which is not exactly landscape architecture but there is some remote relation) there are Fragstat and other software that help computing landscape indices also we can build 3D models expressing actual situation more accurate then ever before. There are some examples of use GIS in landscape architecture that I know, but how in general others landscape professionals do use it? What software do they use? Landscape architecture is a broad and interdisciplinary field. In different countries understood in different way. I do have this idea, that in order to get better understanding of my own profession it would be great to hear about use GIS in landscape research/analysis/design/planning etc.


title: Select "interior" lines of a polygon-to-line layer

text: I have a polygon featureclass that has adjacent features. Using Arc10, I am turning those polygons into lines using the "Featureclass to Lines" tool. Now I would like to select **ONLY THE INTERIOR** line(s) (color=grey) so that I can run some ADDITIONAL processes on them. The key is that some (but not all) feature may touch the boundary line (color=blue). I am trying to use "Select By Location", but I cannot come up with a way to accomplish this task. I might just be overthinking this one, but I am stuck. Any suggestions? ![enter image description here](http://i.stack.imgur.com/ca0l5.jpg)