title: Reclassifying a raster based on values in another raster

text: I am attempting to create a new raster based on the values of two rasters I have already created. I am using ArcGIS 10.1   * Raster1 is an area with my true values.   * Raster2 was created by reclassifying Raster1 using focal statistics. This reclassified the raster1 pixel values according to the maximum pixel value within the kernel. The goal for my reclassification is to compare the values in Raster1 to the values in Raster2. If the Raster1 value is < 50% of the value in Raster 2 then the pixel in the new raster is set to 0, else the value is kept as the value in Raster1. I am hoping that I can work something out using the ArcGIS toolbox. I am not an advanced user so I don't understand programming.


title: How to permanently rearrange the order of features in a Shapfile with QGIS?

text: I have created a trimmed/clipped map of Europe, the Middle East and North Africa in QGIS 1.7.3 by deleting all the other countries from a world map shapefile. I would like to permanently rearrange the countries that are left in the shapefile in a particular sequence in the attribute table. I searched for other similar questions on this website and one answer suggested that the .dbf file be modified using a spreadsheet program, e.g., MS Excel 2003. I tried doing that and it seemed to work; but then I realized that when I pointed at Austria in the shapefile, the name that came up was Egypt. Hence, on the map itself, the names of the countries are all mixed up now. Can you please suggest a way to arrange the countries in the attribute table according to a specific sequence without its affecting the countries on the map itself? Thank you.


title: How do I create measured line segments along existing polylines from a starting point?

text: I realize my problem may be solved using Linear Referencing, and if that's the case I'm looking for help figuring out the workflow involved there. I'm working in ArcGIS 10.2 Desktop Basic, and this is what I started with. I had a bunch of points where there were animal sightings, and those point were mostly in between two other points called "camera locations" that are start & end points for these animal sighting points. These are separate point features. I was supposed to connect the animal sighting points into a generalized line from the start to the ending camera locations, to have a line. I have this for about a dozen different camera start & end locations, a connect-the-dots line in between them now. Now, I'm supposed to start at the "starting" camera point side of the connect-the-dots line, and measure distances from 0-400', 400-800', 800-1200', and 1200-1600' but if the line ends before 1600' just cut it off at the ending camera location point on the line. To be clear, these are like 12 different polylines in a feature (shapefile), and they each need to be broken up into the measures above. I'm wondering if someone can take my problem and help me make sense of how to get those measured points along the lines I drew from connecting the previous dots. While I've been using GIS for years, somehow I've never used linear referencing and I'm not sure which of its tools I need to use and when. Thanks, I appreciate any help!


title: OpenLayers.Popup.FramedCloud won't autosize

text: I don't see why this shouldn't work, but it doesn't.. :P > var popuphtml=""; $.each(persons.persons,function(i,item){   >  if(item.in_loc == feature.attributes['locid'].substring(8)){ popuphtml > +=""+ "" + item.lastName +" " + item.firstName + "   > "; } }); >   >   >               popup = new OpenLayers.Popup.FramedCloud("chicken", >               feature.geometry.getBounds().getCenterLonLat(), >               null, >               popuphtml,   >               null, >               true, >               onPopupClose); >               popup.autoSize = true; >               panMapIfOutOfView=true; >               feature.popup = popup; >               map.addPopup(popup); >


title: symbol created from emf still acts like a picture and blocks out surrounding vector data

text: I have a lot of CAD blocks that I want to convert to symbol markers to use in ArcGIS 10. Everything I read says the emf is converted to vector data, but that does not seem to be the case. Procedure: open Style Manager>Open the Style I created for the blocks>Click on Marker Symbol folder>Right click for new marker symbol>The only option that allows emf is the Picture Marker option. When I bring the emf picture in it's the block on what is effectively a piece of paper and when I use it as my point symbol, the surrounding area is blanked out. Changing the tier level in the TOC is not a feasible option, and I think that is only a bandaid and not a solution. I think what I really need is a character marker but I don't see an option for creating them from an emf file. What am I doing wrong?


title: Returned BBOX of a WMS with varying aspect ratio

text: WMS GetMap queries take a BBOX parameter for the map and a WIDTH and HEIGHT parameter for the image. If the aspect ratio of the image (WIDTH/HEIGHT) doesn't match the aspect ratio of the BBOX (dx/dy) does the spec say it will distort the image into that WIDTH/HEIGHT aspect or does it say it will increase the area of the map returned in the image. The server I'm playing with is doing the latter, but I see mentions on the web of the former behaviour. In the case of increasing the area of the map, does the spec say that either the vertical or horizontal extent of the image will match the requested BBOX (depending on if its too tall or too wide)? I suppose the safest thing is to always call with matching aspect ratios... Here's some test grabs with three different aspect ratio WIDTH/HEIGHT settings and the same BBOX for your comparison:               wget -O test.png "http://www.osmgb.org.uk/OSM-GB/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=OSM-GB:__all__&SRS=     EPSG:27700&BBOX=532685,176236,542182,183918&FORMAT=image/png&WIDTH=400&HEIGHT=400"     wget -O test2.png "http://www.osmgb.org.uk/OSM-GB/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=OSM-GB:__all__&SRS     =EPSG:27700&BBOX=532685,176236,542182,183918&FORMAT=image/png&WIDTH=400&HEIGHT=800"     wget -O test3.png "http://www.osmgb.org.uk/OSM-GB/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=OSM-GB:__all__&SRS     =EPSG:27700&BBOX=532685,176236,542182,183918&FORMAT=image/png&WIDTH=800&HEIGHT=400"


title: Is it possible to add rings using other polygons?

text: I would like to "Add Rings" but using other polygons. Example: https://www.dropbox.com/s/ue1k3b646ivqedo/holes.JPG I would like to delete the two selected polygons and create holes under them in the big polygon where they are. Is that possible? Thanks!


title: how to get the projection from a vectorlayer in qgis?

text: Hi I am trying to set the projection on a raster to match that of a vector point layer. Thus I need to findout what is the projection of a give layer, to use it in the GDAL.Dataset.SetProjection() so that I can create the GEOTIFF with the appropriate projection. How do I do That? (in Python)


title: What is the width of the beach of a little island polygon?

text: This is a "classic" problem, **that have no practical utility** (see note-1 for applications), but illustrate another GIS problems... ![little island figure](http://i.stack.imgur.com/xj7i3.png) We can do some approximations to modeling the situation,   * the real center (coconut) is not at the ST_Centroid point, but can be approximated to;   * the real shape (a island) is not a circle, neither a square, they are optional references. So, **translating the question to an equation** , and using OGC standard functions... " **What is the value of _w_?**" at the equation               ST_Area(ST_Buffer(ST_Centroid(littleIsland), w, quad_segs)) = ST_Area(littleIsland)      When `quad_segs=1`, the reference is a square, when `quad_segs>=8` the reference is near to a circle (see examples of generated polygons). * * * NOTE-1: For understand this problem into a context of applications, see "Is there a st_buffer inverse function, that returns a width estimation?" and see "How to fit a polygon to a “best buffer” about a reference geometry in it?". NOTE-2: the "inverse function" is so general and complex, because can operates with point, line and polygon inputs. Here we can discuss the simplest case, when the input is like a point, with no area or length.


title: Miller cylindrical projection - central meridian

text: I'm new to projections topics and I've got a question that may be silly, but I hadn't found answer so far. I'm trying to use Miller's cylindrical projection for converting osm lat/lon to x,y values. I don't know how should I choose central meridian. It should be calculated for every transformation I make or I choose it once? What should be value of it? I use this code to calculate miller's projection:               import sys     import math          lon = float(sys.argv[1])     lat = float(sys.argv[2])               x = lon # - CENTRAL_MERIDIAN?     y = 1.25 * math.log(math.tan(math.radians((math.pi/4.0) + (2.0/5.0*lat))))


title: How to make RASTER from irregular point data without interpolation

text: I was trying to make a raster image from an irregularly spaced point database. The data looks like-               > head(s100_ras)              x       y         z     1 267573.9 2633781 213.29545     2 262224.4 2633781  69.78261     3 263742.7 2633781  51.21951     4 259328.4 2633781 301.98413     5 264109.8 2633781 141.72414     6 255094.8 2633781  88.90244      I want these 'z' values within a mesh which I created by               # Create a fine mesh grid     my_mesh=expand.grid(seq(min(s100_ras$Y),max(s100_ras$Y),l=100),                         seq(min(s100_ras$X),max(s100_ras$X),l=100))      I also want the z-values to be assigned as 'NA' for those mesh points that are outside the data points. The points over the mesh looks like this: https://drive.google.com/file/d/0B6GUNg-8d30vYzlwTkhvaHBFTnc/edit?usp=sharing when I plot               plot(my_mesh)     points(s100_ras$Y, s100_ras$X, pch="*", col='blue')      The problem is that I'm not sure how to build on this, the following steps don't work because my mesh grid and data points are not of the same scale!!               library(rgdal)     library(raster)     xyz<-cbind(my_mesh, s100_ras)     r <- rasterFromXYZ(xyz)     image(r)      If I try to make a raster just by using the data points (without any mesh), R throws an error since my data is irregularly spaced!               library(sp)     s100_ras <- data.frame(expand.grid(x = s100_ras$Y, y = s100_ras$X),                             z = as.vector(s100_ras$mean))     coordinates(s100_ras) <- ~x+y     proj4string(s100_ras) <- CRS("+proj=utm +zone=46 +datum=WGS84")     gridded(s100_ras) = TRUE          suggested tolerance minimum: 0.916421      Error in points2grid(points, tolerance, round) :        dimension 1 : coordinate intervals are not constant      Moreover, I was trying to play with 'rasterize' function (for irregular grids) of 'raster package', but couldn't get a way with it :(. I know how to interpolate and make a regular grid, but for the sake of originality, I want to AVOID interpolation. Is it possible to make a raster of irregularly spaced data points without idw or kriging methods? Thanks in advance.


title: Extraction of stream gradient for branched river networks

text: At the moment I am working on the extraction of stream gradients for a large basin, as a consequence, rivers are strongly branched.So far, I read through the relevant posts on this site, but so far I am not able to apply the proposed methodologies to my problem. I started with the common procedure of extracting the river network from the 30-by-30 m AsterDEM. What I want to do now, is to measure the gradient of the resulting streams in fixed intervals. The interval for a river reach depends on its mean gradient, which I determined using HEC GeoHMS (the steeper the mean slope, the more often I want to sample). Thus I divided the streams in four categories. Now, I need to divide the resulting streams in equal intervals. I wanted to make use the technique as proposed here: 1. Thus converting the stream network into a route (using the “linear referencing toolbox”) and then use route events to divide the route in reaches of equal lengths. Unfortunately I have no idea how to prepare an event table in order to divide the routes in equal intervals (e.g. one event every 100 m), especially in a branched network. I would highly appreciate if you could guide me through the workflow, as I have no experience with the linear referencing tools! Rafael


title: How to use Collect Values and Tabulate Areas in ModelBuilder?

text: I have problems with the **Tabulate Area** tool. The polygons in my shapefile overlap and because of this ArcMap does not calculate proper values. I try to solve this using ModelBuilder where it is possible to do the **Tabulate Area** for each polygon separately but at the end I would like to combine rows back to one table. I thought that the solution may be the model only tool called **Collect Values**? Can someone help me how to add rows to one table or to save each **Tabulate Area** result under different names? I hope that the question is understandable. If not, I will try to rewrite it. Here is also a picture. ![enter image description here](http://i.stack.imgur.com/0bNmm.jpg)


title: How to choose an area (of mxn metres) within a raster that have values above a certain range

text: I've got a raster of concentrations of iron on Mars and I want to be able to choose several areas of either m by n dimensions in metres, or just any area in m^2, that encapsulate regions that have high concentrations. For example: I want all regions with concentrations over 5 for an area of 50x50 metres, or an area of 2500 m^2. What process would you use? Would I create a grid with 50x50 metre cells? Then perform some sort of calculation that tabulates the average cell concentration and output concentrations over some value? I don't have any idea how to do that, but even if I did the areas don't need to be perfect squares so it wouldn't be the best solution just better than not do any analysis. Either way if anyone has any ideas how to choose areas like this either way I'd really appreciate any help. Thanks


title: Creating a feature of a part of a state boundary from a shapefile of U.S. States

text: I have a shapefile of U.S. states. I want to create a feature layer from it so that I can do a Near analysis with it. Suppose that I am interested in the Georgia boundary, but only the part of Georgia that touches South Carolina? How could I extract just that segment? Or, suppose I am additionally interested in the part of Georgia that touches both North Carolina and South Carolina but not at all Alabama, Tennessee, or Florida? Thank you very much for any help.


title: Create 6 band IMG file with GDAL out of two JPG files, one with WorldFile coordinates

text: I have two JPG files # 1 ortomap: http://imageshack.us/photo/my-images/703/ortofoto1.jpg/ with WorldFile:               0.683772574816354     0.000000000000000     0.000000000000000     -0.683739485024190     395057.114562241190000     741149.458682749890000      # 2 aerial photo in form of JPG: http://imageshack.us/photo/my- images/688/aerialphoto.jpg/ **How can I create a 6 band IMG out of those two JPGs but according to the coordinates from the first JPG WorldFile using GDAL?** I am after some exemplary commadline commads


title: formatting geojson from postgis data using node.js

text: I am trying to build a GeoJSON object from a SQL query to some GIS point data in a postgis postgresql database. A snippet of my node.js app.js is below. As it stands I understand building the type and features, but don't know how to attach a properties properties to each GeoJSON record (in the below, it all renders at the end, separate from (not collated with) the features). THE QUESTION: What do I need to do so that the properties attach (collate) for each record in the loop that builds the GeoJSON So that it looks more like this http://www.geojson.org/geojson-spec.html#examples?               function GrabData(bounds, res){            pg.connect(conn, function(err, client){            var moisql = 'SELECT ttl, (ST_AsGeoJSON(the_geom)) as locale from cpag;'                 client.query(moisql, function(err, result){         var featureCollection = new FeatureCollection();              for(i=0; i<result.rows.length; i++){           featureCollection.features[i] = JSON.parse(result.rows[i].locale);           featureCollection.properties[i] = JSON.parse(result.rows[i].ttl); //this is wrong        }             res.send(featureCollection);        });          });     }           function FeatureCollection(){        this.type = 'FeatureCollection';        this.features = new Array();        this.properties = new Object;  //this is wrong      }


title: Convert X Y Pixels to Easting Northing

text: Hopefully someone can help. I'm trying to find out what the formula would be to convert the X, Y of an image to the Easting and Northing coordinates. I've downloaded some of the raster maps from OS UK web site and I'm looking at adding the miniscale tif file to one of my applications. I've been given a number of coordinates in the easting northing format and I need to convert this to an x y of the image so I can place a label or image on top of the mini scale tif file. I've been searching online quite a bit but everything seems to point towards applications like ArcGIS and other GIS software. However I'd quite like to make this a simple project, so I'd like to avoid embedding 3rd party sources within the application. Does any body have any good starter guides or examples on this, or what the forumla would be to work out the x y from the easting and northing I have? Thanks in advance.


title: Search for geojson polygon in mongodb implementing a custom vector tileserver

text: I am building a vector tileserver in node.js using mongo 2.4.7 as database. My database contains geojson data from multiple sources and they are in many layers. It is important to render the layers in the correct sort order. Since I have little, or no control of the client rendering the individual tiles I have to clip the geojson data before returning it to the client so one “lower lying” polygon does not hide another higher lying tile in an adjacent tile. I have created an simple rectangle ->polygon clipper that works for now - but this should be improved later on. If anyone knows an efficient clipper library that is usable from node.js please let me know. This is, however, not the issue. My real challenge is the finding the correct polygons within the database. I am currently using both $geoIntersects and $geoWithin in my query. Let say that we have a polygon A that in zoom level x just fits within one tile bounding box (tile A). Here the Polygon will be returned by the $geoWithin clause and all is well. Going several zoom levels in, the area from the tile A now fits within 9 tiles in a square. How do I find Polygon A when the server needs to render tile number 5 (the middle one)? And yes - in this situation I would like to render a square covering the complete tile due to the clipping. I have created a drawing to illustrate my challenge. ![enter image description here](http://i.stack.imgur.com/aK8zS.png)


title: How to split one modelbuilder process in several GUI windows?

text: ![enter image description here](http://i.stack.imgur.com/YYzp3.jpg)I am working with ArcGis 10.2. I would like to know how I can have the different steps of a modelbuilder's process displayed in successive windows once the previous main step is done? A bit like in a AV script. Let's say, I have a process with 4 main steps (each has at least 6 to 10 Tools/steps) requiring each time different Inputs, Parameters and Outputs. What I will get as "GUI" tool will be a very long list of parameters to inform. I would like to avoid this in order to have an easier tool to handle by "basic" ArcGis users. I am quite Ok with the modelbuilder but sadly don't have skills in coding. Many thanks I have edited my post with a screenshot. I'd like to have the possibility to ask the user if the first Step (in the middle) needs to be done or not - regarding the answer either run the second step (last box), after execution of the first step or directly after the question " do you need to reclass your data?" I hope this is helpful. Many thanks


title: How to read Oracle data with PyQGIS? ("identifier is not long" error)

text: Do you have any idea how load Oracle spatial data to QGIS with using the PyQGIS? Is it possible? I try:               uri = QgsDataSourceURI()     uri.setConnection("*host*", "1521", "*database*", "*username*", "*password*")            sql = "(SELECT fid, geom FROM target.ms_ft_bh095_v)"     uri.setDataSource("",sql,"GEOM","","fid")     layer = QgsVectorLayer(uri.uri(), "test", "oracle")      The error message is: ORA-00972: identifier is not long. But FID max value is cca 20 000. Is Oracle driver correct with this syntax? Ms_ft_bh095_v is view in this moment.


title: Maintaining times with points2one in QGIS?

text: I have a problem that might be trivial, but I really cannot find a solution! I have a csv files with several points (x,y). These points represent the position of different animals at different times, and thus can be joined in several tracks (each point belongs to a single track, but a track can have from 2 to several points). Now these points have clearly also a time, but if I join them using Points2one I loose this information that I need to calculate the speed of the animals. Is there a way to have a multime time column (or at least first-last time) or to calculate the speed of an animal if I have such dataset?


title: Can i use CreateMapSDDraft in the arcpy.mapping library to set the scale range on tiled services at 10.1?

text: I know that esri has exposed many properties to the acrpy lib for publishing services. I need to published a tiled/cached map service but i need to set the min/max scale when i do that. Can this be accomplished via python code? I don't see it as a property on the help page of the CreateMapSDDraft script. Is it just not supported or am i doing something wrong?


title: IF-Statement in Field Calculator

text: I want to use an IF/ELSE-Statement to define a field in the attribute table of my raster. Therefore I have a field "VALUE". The condition should be: If "VALUE" < 4 ...then use die field "A".....if "VALUE" > 4 then set 0.5. I was searching for answers but couldn't find a proper one.


title: Joining by location - DIVA-GIS .dbf layer not included

text: I want to get place names from a .dbf file displayed on my map. I've imported the .dbf file as a vector layer. The gazetteer table does not have any ID fields so I was trying to use the "Join by location" option. I've made sure that the same CRS is set on the gazetteer as on all the other layers. However, when I open the "join by location" dialog the gazetteer layer is not included into the list of the available layers. Why is it not displayed? Maybe this has to do with how the fields containing the coordinates are named in the table?


title: QGIS relations widget in custom form

text: I've just begun exploring QGIS relations option, now in version QGIS 2.2, details here: http://blog.vitu.ch/10112013-1201/qgis-relations This is a great implementation of something I've been attempting to manually organize for a while now. Now, i'm wondering how to integrate the relations widget into a custom form (currently, it automatically adds only if using the auto generated form). I'm new to PyQt Designer, so i'm not sure where to start. Thanks in advance.


title: RemapValue Usage

text: I successfully produce and empty cost raster at the end of this script. I am fairly confident it has something to do with RemapValue parameters? Does anyone know why my out put cost distance raster is empty?               #-------------------------------------------------------------------------------     import arcpy     import os     import time     from arcpy import env     from arcpy.sa import *     from decimal import Decimal     arcpy.CheckOutExtension("spatial")     env.overwriteOutput = 1          #Begin code     #To be ommitted ***"Hardcoded Testing"***     env.workspace = r'\\Scratch.gdb'     landuseRaster = r'\\Landuse_Raster_Clip'     subs = r'\\Subs_w_dist'     Field = 'Name'     destValue = 'Coulter'     originValue = 'Amarillo 34th St'          #Where clause sql expression     where_clause_destination = "\"{}\" = '{}'".format(Field,destValue) #produces: "Name" = 'Coulter'     where_clause_origin = "\"{}\" = '{}'".format(Field,originValue) #produces: "Name" = 'Amarillo 34th St'          #Make feature layer     arcpy.MakeFeatureLayer_management (subs, "origin_Selection", where_clause_destination)     arcpy.MakeFeatureLayer_management(subs, "destination_Selection", where_clause_destination)          # Process: Slope     arcpy.AddMessage("Computing slope...")     arcpy.gp.Slope_sa(raster_Clip, 'slope_Raster', "DEGREE", "3.2808333")          # Process: Slice     arcpy.gp.Slice_sa("slope_Raster", 'slice_Slope', "10", "NATURAL_BREAKS", "1")     slice_Slope = "slice_Slope"          remapSlope = RemapValue([[1, 1], [2, 2],[3, 3], [4, 4], [5, 5], [6, 6], [7, 7], [8, 8], [9, 9],[10,10],["NODATA", "NODATA"]])          remapLanduse = RemapValue([['Cropland and pasture',2],['Herbaceous Rangeland',1],['Nonforested wetland',10],['Other agricultural land',1],['Strip mines, quarries and gravel pits',2], ['Mixed rangeland',1],['Lakes',10],['Residential',3],['Commercial and Services',3],['Confined feeding operations',3],['Transportation, communications and services',10], ['Transitional areas',3],['Shrub-brushland rangeland',1],['Other urban or built-up land',3],['Bare exposed rock',3],['Industrial',2],['Reservoirs',10],['Forested wetland',10],['Mixed urban or built-up land',3],["NODATA", "NODATA"]])          myWOTable = WOTable([[landuseRaster, 80, "LANDUSE", remapLanduse],                          [slice_Slope, 20, "Value", remapSlope],                          ], [1, 10, 1])          # Process Weighted Overlay     arcpy.AddMessage("Creating weighted raster...")     outWeightedOverlay = WeightedOverlay(myWOTable)     outWeightedOverlay.save("weighted_Raster")     weighted_Raster = "weighted_Raster"          # Process: Cost Distance     arcpy.gp.CostDistance_sa("destination_Selection", weighted_Raster, "distance_Raster", "", 'backlink_Raster')


title: English-taught bachelor programs in cartography and GIS science?

text: Is there any English-taught bachelor program in Cartography and GIS science in Belgium, the Netherlands, Austria, Germany, Denmark, Sweden, Finland? I've done massive homework about that, but haven't succeeded. If you know any course or program like that, please let me know. As the only international language that I know is English, therefore I can only enroll on English-taught programs. Moreover, programs in USA, UK and Canada are unaffordable to me, so I want to do my study in the above countries.


title: r.stream.* commands do not work in qgis environment

text: I have been working for a while with some grass commands in QGIS environment (r.watershed, r.terraflow, r.thin, r.to.vect etc) with out any problem. I haven't been able to get any results at all from r.strem.xxxxx commands. I use for input the results from r.watershed but I keep receiving the Ooops.... message. Do I need to install extra files for the r.stream extension?? Any help highly appreciated. The message received is the following: > Oooops! The following output layers could not be open Output basin map > (Strahler): > C:\DOCUME~1\aant\LOCALS~1\Temp\processing\e3cc4eb40cfa41c09e638f569ef4b6ac\strahler.tif > Output basin map (Shreve): > C:\DOCUME~1\aant\LOCALS~1\Temp\processing\a921888aca254949823497f5e2f85056\shreve.tif > Output basin map (Horton): > C:\DOCUME~1\aant\LOCALS~1\Temp\processing\dfdd7cf3acf841858803e77019a0a0b3\horton.tif > Output basin map (Hack): > C:\DOCUME~1\aant\LOCALS~1\Temp\processing\04cabe9871454d57a38083575feaac64\hack.tif > Output basin map (Top): > C:\DOCUME~1\aant\LOCALS~1\Temp\processing\d1dc1d81119e472fbcaf0a468b7bd8d6\top.tif > The above files could not be opened, which probably indicates that they were > not correctly produced by the executed algorithm Checking the log > information might help you see why those layers were not created as expected > This algorithm requires GRASS to be run. A test to check if GRASS is > correctly installed and configured in your system has been performed, with > the following result: GRASS seems to be correctly installed and configured I run QGIS ver.2.4 on WinXP and Win7.


title: Identify duplicate polygons to enable removal after reviewing attributes

text: I wish to label overlapping polygons so that I can manually locate each overlap then review the attributes of each overlapping polygon and decide which suite of attributes to maintain before dissolving the polygons to ensure a dataset free of overlap. I know the manual review is time consuming but I cannot see another way of doing this if you don't know which attributes should be maintained before looking at the overlapping polygons. I wish to be able to calculate the relational area of specific attributes within the dataset without biasing the analysis due to overlap. Ideally I would like a way to add a field that added a code (e.g. number) that would link each polygon with other polygons that it either entirely or partially overlaps with. Then I can determine which of these polygon's attributes to retain. I have very large datasets with numerous different attributes. I use Arc10 with win7 64 bit processor and ArcInfo licence. Any ideas and help is greatly appreciated. Thanks


title: How to export 1.5 million records from File Geodatabase to Microsoft Access table

text: I have a table that I want to export from arcmap that has around 1.5 million records. I need to import it into access. Now every time I have tried to export it into a format that access will read, like a csv or dbase, the records become corrupted and instead of displaying there values i they have contents such as "3.186345345345-3". Does anyone know of a way to export to access that will work perfectly? Thanks.


title: Some help installing pgrouting_dd on mac os x

text: I hope someone may be able to help me but i've been grapling for hours trying to install pgrouting with Driving Time to no success. I have been able to build pg routing no problem. However if i trying to cmake with the -DWITH_DD=ON flag set when performing the make I get the following and it fails:               [ 10%] Building C object         extra/driving_distance/src/CMakeFiles/routing_dd.dir/alpha.o     [ 20%] Building CXX object extra/driving_distance/src/CMakeFiles/routing_dd.dir/alpha_drivedist.o     [ 30%] Building CXX object extra/driving_distance/src/CMakeFiles/routing_dd.dir/boost_drivedist.o     [ 40%] Building C object  extra/driving_distance/src/CMakeFiles/routing_dd.dir/drivedist.o     Linking CXX shared library ../../../lib/librouting_dd.dylib      Undefined symbols:.... Does anyone have any ideas on how to overcome this? I should also point out that I am using os x 10.6.x C


title: Is it possible to access individual pages within an existing pdf using arcpy.mapping?

text: I've used the arcpy.mapping module for some time now to create mapbooks and other pdf documents, using a combination of Data Driven Pages and Python. I would like to take an existing pdf and save each individual page as separate pdfs. I have read through the online documentation for `PDFDocument`, and while I can see the usefulness of the `pageCount` property, I don't see how I can access an individual page. The only reference I see to a page number is when using the `insertPages` method to insert a page before an existing page. I know that using something like `pypdf` enables one to iterate over the pages in a pdf and other useful things, but I'd like to know if I can do this using `arcpy` before downloading it.


title: Extraction of DEM from netCDF file into ASCII XYZ format

text: I need to extract a DEM of a waterbody from a netCDF file so I can ultimately convert the DEM into xyz format. I'm not sure how to go about doing this, and so far, I have tried to create raster layers of the DEM in ArcGis. However, I'm not sure if this is the best approach to obtain the data in xyz format. If anyone could provide some ideas about how I would go about handling this conversion I would appreciate it. Thanks, Michael


title: How to design markers - raster or vector?

text: I have several layers that are point geometries. I am defining them custom markers symbols and such symbols have colors (red, yellow and green) to represent values. Since there are several of them in an area, I started with marker clusters, but quickly realized that clusters are good only if all markers mean the same (like redfin home markers). I would like to stick to markers but to avoid visual clutter - I would like to resize them based on the zoom level. If I design these markers in SVG, I could achieve this - as its easy to size SVG than raster formats, similar to this link http://mapbox.com/tilemill/docs/guides/advanced-map- design/#images_as_icons Please send me any examples or sample sites you have developed or come across so that I can see the visual appeal of having several custom markers in maps. Thanks, Ramesh


title: Parcel Fabric vs. Polygon feature class (for land parcel management)

text: I'm wondering if I should migrate my existing parcel polygons (stored as a feature class in a geodatabase) to a parcel fabric in ArcMap. What are the advantages of doing this? I work for a county so the parcels are used mainly to store attributes (who owns which parcels in a graphic representation etc.) and general information - NOT for being survey accurate. Any info is greatly appreciated. Thanks!


title: add location search to openlayers map

text: I am trying to figure out how to add a search functionality to my openlayers map. I want the user to type in an address and the map to zoom into the location. I am looking at different geocoders and can't figure out a simple solution to this. Nominatim seems liek a possible solution. But I am trying to find something that I can develop on the front end. I have OSM data for my background and would prefer to use this. Any tip or ideas to point me in the right direction would be appreciated! Thanks


title: OpenLayers WPS function fails in IE

text: I am making a WPS call in OpenLayers which works successfully in Chrome and Firefox but not in IE9. The response from IE tells me "Process failed during execution\nCould not convert to target type com.vividsolutions.jts.geom.Geometry" The same code returns intersecting polygons in both other browsers so I think it's down to how IE handles my GML/XML. Can anyone help out? Code as follows: function wpsExecute(featureA, featureB, identifier) {               var polygonElement;     if (BrowserInfo.browser == "Chrome") {         polygonElement = "Polygon"     } else {         polygonElement = "gml:Polygon"     }          // set up the process specific parameters     var options = {         // set the wps process         identifier : identifier,         dataInputs : [ {             identifier : "a",             data : {                 complexData : {                     attributes : {                         mimeType : "text/xml",                         subtype : "gml/3.1.1"                     },                     // get the gml:Polygon node from the xml DOM                     value : featureA.getElementsByTagName(polygonElement)[0]                 }             }         }, {             identifier : "b",             data : {                 complexData : {                     attributes : {                         mimeType : "text/xml",                         subtype : "gml/3.1.1"                     },                     // get the gml:Polygon node from the xml DOM                     value : featureB.getElementsByTagName(polygonElement)[0]                 }             }         } ],         responseForm : {}     };          var format = new OpenLayers.Format.WPSExecute();     // console.log(format.write(options));          try {              OpenLayers.Request.POST({             url : WPS_URL,             data : format.write(options),             async : false,             success : parseExecute,         });      WHERE featureA and featureB are as follows:               var gmlOptions = {         featureType : "feature",         featureNS : "http://example.com/feature"     };          // new gml 3 format     var gml3 = new OpenLayers.Format.GML.v3(gmlOptions);               compareFeatureGml = gml3.write(compareFeature);          var xmlFormat = new OpenLayers.Format.XML();     var featureB = xmlFormat.read(compareFeatureGml);      Thank you all!


title: Does ArcMap 10 have a command line switch that is roughly equivalent to Visual Studio's /SafeMode?

text: Anytime I had an ArcMap 9.x related incident open with ESRI, they always had me do the usual steps of uninstalling extensions and tools, renaming the ESRI profile folder to ESRI_old, and renaming the ESRI registry key under HKCU to ESRI_old. I am unable to find a list of arcmap.exe command line switches, but I feel like a /SafeMode type of switch would really speed up troubleshooting. The only command line switch I know of is `/log`. Does anyone know of a /SafeMode equivalent or have a link to arcmap.exe command line switches?


title: Cycle through / add name of clipping feature to output of clips in Model Builder?

text: I have a feature class named Study Areas with multiple polygons, named in an attribute field as Study Area A, B, C, etc. I will be using these polygons to clip various polygon feature classes that represent different landscape attributes, like Land Cover, Soils, etc. The end goal is an analysis of how much of each type of Land Cover acreage is within Study Area A, B, and so on. I would like to do this through Model Builder (I don't know Python) but need to find a way that I could create outputs that have the Study Area name appended onto them so that the outputs can be easily organized (there's going to be hundreds of outputs). So running through the model, I would need 1) a component that cycles through the Study Area feature class so that only one study area polygon is used for clipping at a time and 2) a clipped output named something like Study_Area_A_Land_Cover Thanks


title: Reverse Geocoding - calculate every time, or store result in database?

text: Several questions, such as this one, discuss whether to store only lat/long in the database & calculate the street address every time it is needed, but they seem to miss a few points & seem to be opinion, rather than authoritative. Hence this question, seeking a final, definitive clarification, with justification. 1) reverse geocode every time   pro - smaller database size (why should everyone duplicate Google's database?)   cons - might occasionally fail whereas storing the result will not 2) store address in the database   pro - less likely to hit a Google imposed limit on the number of requests   con - could be a _large_ amount of data. A few hundred vehicles storing a new location every 5 minutes or so, 16 hours per day... probably more coding If I take approach 2), I plan a new MySql table for addresses, indexed by lat/long and storing the street address. Every time a vehicle reports position, look for a matching address. If none, calculate & add. In addition: later, when a client requests movement history, if the address is blank, attempt to calculate & store if found. Please note that I am restricted to a MySql database. Thanks in advance for any help


title: OpenGeo Suite start up failure

text: I started the openGeo Suite 3.0.1 in Windows 8 and after a while I get this error: `Giving up on starting the Suite as it is taking longer than expected` Any idea?


title: How to iterate all input folders (containing coverages) using Python with ArcGIS 9.3?

text: I'm trying to write python scripts to batch convert thousands of coverages to shapefiles in arcgis 9.3. The coverages files are in different folders:               d:/area_A/building     d:/area_B/building     d:/area_C/building     ...      and the output file will be merged or unioned into a single shapefile:               d:/output/building_merged.shp      I'm stuck in iterating all the input folders..here is what I've got so far:               import arcgisscripting     gp = arcgisscripting.create()          try:             gp.Workspace = "d:/area_A/building"            gp.FeatureClassToShapefile("polygon","d:/output")       except:              print gp.GetMessages()      Thank you in advance to point me into the right direction.


title: "GIS" Mapping on Monument - Scaling

text: I have a Baroque tombstone for which I have to map the damage phenoma. Normally one can use a vector drawing programm but one has always two boundaries between the different polygons. Questions: Is it possible to add a layer of a photograph of the tombstone and to scale it afterwards in Quantum GIS? If one adds a layer, QGIS asks for the CRS and, of course, tombstones has no CRS. Besides the „boundary question“ it should then be possible to calculate the percentage of the different damage phenomena which would be very nice. Thanks in advance, Jörn


title: Advice on ideal server specifications

text: I am looking for your thoughts on ideal server specifications for a moderately sized habitat mapping project. My current server is about to reach its storage capacity, so I am researching buying a new one that I can dedicate to the project. The following is a list of my basic server requirements:   * 10-12 TB Storage Capacity   * RAID 5   * $5k limit   * Primarily used to support ArcGIS, Erdas, and ENVI applications I would like to get your thoughts on system configurations that work well for you, given the previously stated constraints. Any advice on specific hardware and server software would be greatly appreciated. Thanks for your input.


title: Pyscripter hangs on spatial join?

text: I've got a script which loops through a series of FC layers (tracker points). Using a FOR loop, they're reprojected then joined to a streets layer. Sometimes it works fine, then other times it never completes the spatial join and stays hanging in debugging mode & I have to manually stop it. It does the reproject OK, just not the join. Any idea why it might be doing this? Cheers               for tracker in trackerPoints:              trackerLoc = workingFolder +"\\" + tracker              gpsReproject = gpsID + "_Reproject"         arcpy.Project_management(trackerLoc,gpsReproject,"PROJCS['TM65_Irish_Grid',GEOGCS['GCS_TM65',DATUM['D_TM65',SPHEROID['Airy_Modified',6377340.189,299.3249646]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',200000.0],PARAMETER['False_Northing',250000.0],PARAMETER['Central_Meridian',-8.0],PARAMETER['Scale_Factor',1.000035],PARAMETER['Latitude_Of_Origin',53.5],UNIT['Meter',1.0]]","TM65_To_WGS_1984_2","GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433],METADATA['World',-180.0,-90.0,180.0,90.0,0.0,0.0174532925199433,0.0,1262]]")              gpsJoin = gpsID + "_Join"         arcpy.SpatialJoin_analysis(gpsReproject,streets,gpsJoin,"JOIN_ONE_TO_ONE","KEEP_ALL","","CLOSEST","#","Distance")  # Spatial Join streets & tracker points


title: How to set up Python 2.7 to work in both QGIS and ArcGIS?

text: I have both ArcGIS 10.1 and QGIS 1.8 installed and both include Python 2.7 in slightly different setups and in different folders (c:\Python27\ArcGIS10.1\ and c:\Program Files (x86)\Quantum GIS Lisboa\apps\Python27). It is only ArcGIS which include python.exe, while it is really QGIS I would like to start experimenting with for some scripting, but preferably just running my scripts from the command console. Are there any recommendations on how to do this smartest? I'm thinking of settings environment variables etc. I am new on both QGIS and python.


title: The relation between the Feature's Geometry and its externalGraphic

text: I wonder what is the relation between the `Geometry.Point` one sets for a `Feature.Vector` and the image provided as `externalGraphic` for the style, in OpenLayers. Is the `Geometry.Point` of a feature for the image:   1. The centre point of the image,   2. or the left top most corner of the image. In case it is 2 then please advice how would I set the values of offsets if, for example, the image is 32x32 (WxH). Is it simply 16 (W/2) for the value of `graphicXOffset`, and accordingly for the `graphicYOffset`? I need to know this information to set the image centre to be the point of a feature. Sadly the documentation doesn't say a word about it, or I am unable to find it.


title: How to get server info from WMSCapabilitiesStore

text: I am using the following code to retrieve a list of data layers hosted on a WMS Server:               var store = new GeoExt.data.WMSCapabilitiesStore({ url: fullUrl });          // Add some callbacks to handle various situations       store.on('load',      successFunction);     store.on('exception', failureFunction);          store.load();                function successFunction(store, records)      {         // Generate a nice table or something     }           That works well enough for retrieving layer data, but while I am connecting to the server, I also want to get the service/title attribute for the server itself. That data is contained in the server's response to the GetCapabilities request, and it must be stored somewhere, but I can't figure out how to retrieve it from the data returned to successFunction.


title: Applications and methods for TerraSAR X

text: People who already used or are familiar with TerraSAR X, can you point out:   * Main applications and methods   * Opensource software tools Thanks


title: Failed to create reader using geoserver Image Pyramid plug-in

text: I am following instructions on this link, http://docs.geoserver.org/stable/en/user/tutorials/imagepyramid/imagepyramid.html I am pointing the URL to the folder normally but I get the following error, "Could not list layers for this store, an error occurred retrieving them: Failed to create reader from file:///mnt/MapData/bmpyramid and hints null" The difference is that in my case, I used a mosaic.vrt file to create the image pyramids. If that was the problem, what would the solution be?


title: Printable Esri symbol sets, a.k.a. style references

text: Does anyone know of good printable documents of the as-shipped Esri symbol sets, also called Style References? The new with Arcgis 10 ability to search symbol sets is a great boon, but one can still only find things for which you know the keywords for, and only few so many at a time on screen. I find hard copy of things like this much easier and faster to look through and locate the right one for the map at hand. For example this morning I was looking for a good glacial moraine symbol and searching the style sets didn't locate anything relevant. So I decided to bite the bullet and build my own. In the process I went browsing through the _...\Desktop10\Styles_ directory and discovered several pictures I've never noticed represented in any styles that would be perfect for the project I'm working on. Since they're here they very likely already exist in a symbol set and I'd be wasting my time creating them anew. Thus the desire for a printed set to help me find them. In previous versions there is a developer sample called Style Dump, but it doesn't work in Arcgis 10 (and the results were pretty crude and wasted paper anyway).


Our GIS Server setup includes 1SOM and 2SOC servers. I am getting this error in Server log. I have tried several settings but none seem to fix the issue, looking for some feedback on what is the actual fix. Error1: Server Object instance creation failed on machine Error2: Server has no associated SOC machines or all SOC machines are at capacity. Add a SOC machine or change the capacity of the available SOC machines Thanks Jay


title: Where can I find a list of non deliverable zip codes

text: Here is an interesting problem I am facing... There are certain areas within the US where the USPS does not offer postal delivery to those addresses. One example here is my Zipcode 92067 which is "Rancho Santa Fe". All homes within the 92067 zip code can only accept regular postal mail when such mail is sent to ones PO box. So, my question here is: where can I find a list of all zipcodes or for that matter zip4s which don't have rural delivery. Thanks in advance! * * * Based upon the data from Jeffery I created the following dataset, mapped and plotted the results based on those zipcodes. Here are the results of that work: https://www.google.com/fusiontables/DataSource?snapid=S441148K1ew


title: Arcgis desktop connection to oracle

text: I am trying to create a database connection from ArcGIS desktop 10.1 to SDE Geodatabase on Oracle server with direct connection using Oracle easy connect syntax So, from ArcCatalog I use "Add database connection" then I chose Oracle in database platform then in the instance I write down myserver: port/serviceid so I write dbserer:1539/devDatabase but it gives me ora-12154 TNS:could not resolve the connect identifier I thought by using the easy connect I don't have to worry about TNS names


title: Are there any decent free GIS projects built in C# or with a good C# API?

text: I'm looking for a free mapping / gis project with a **C#** interface, preferably written in **C#** and open source. Does such a project exist? If so, are there any that do not require connecting to a map server and will allow me to distribute map tiles on disc?


title: Displaying cross-sections on a web page

text: I'm looking for inspiration and possible tools on how to display cross- sections on a (dynamic) webpage. More specifically I have to show on these cross-sections how underground pipelines are located compared to the surface. Apart from the 2D information on the surface (high resolution) we will also have access to lidar data of the surface around these pipelines. I'm basically considering two paths. One is to make a true cross-section, like eg this example: ![enter image description here](http://i.stack.imgur.com/KeGCR.jpg) Another option would be to create a small strip of land, a bit like this image (but with a much more boring look, but and showing the depth of the pipelines also. ![enter image description here](http://i.stack.imgur.com/MO1Md.jpg). This has the advantage that the surrounding terrain is also shown, which is useful as the position of the pipelines is not always very precise. Although I have some ideas on how to create the first one (basically making a projection along the transect), I wonder if any tools exist that could create images like the second one. And in general I'd be interested to know if anyone has already knows of (or made) similar applications.


title: How can I create a WFS from netCDF file(s)?

text: Currently I'm working on a project that uses OpenLayers and Geoserver to view data from a postGIS server (both on a map through a WMS, and as data through a WFS). Now, I'd like to be able to compare the output that I currently have with data from netCDF file(s). So far I've tested ncWMS and used it as a WMS server, but as far as I know, if I want to get the raw data WMS supports getfeatureinfo, but here I can only display data from 1 pixel at a time, so to compare a lot of data would resolve in a huge amount of requests to the WMS server. I'd like to find an alternative, where it would be easy for a user to compare his/her data from netCDF files with the data in the database. Of course users could download the raw data as is now, and then compare it somehow on their own, but that's not what we're interested in. What Im thinking: the user gets an option to upload the file(s) to the application and then the application extracts whatever data it wants. I was perhaps thinking it would send the file to a server that would supply a WFS to the application. I know it is possible to upload the netCDF to a postGIS database, but that seems like overkill? Or is it perhaps the most efficient way? I would also like an option that is reliable and not dependent on third-party software that might now be sufficiently tested or up to date. I'm quite new to all of this and would really appreciate all help I could get :) **Edit:** I haven't figured this out yet. What I've moved onto now is to test the netCDF API for Java from uninet, see http://www.unidata.ucar.edu/software/thredds/current/netcdf- java/v4.2/javadoc/ucar/nc2/package-summary.html However, I'm having difficulties understanding how to extract point data from the arrays with raster data that I seemingly get out of it.


title: Can a custom javascript web application be uploaded and run on ArcGIS online?

text: I'm researching arcgis online and have a subscription account. From what I can make out from the documentation and some experimentation, I've found the following:   * You can create a feature service or tiled map service and have ArcGIS online host it (but not a dynamic map service). You can reference from anywhere using the REST service.   * You can create a "web map" that gives you a viewer and simple widgets, scale bar, zoom, etc. Using the ID of the web map, you can reference anywhere either arcgis online or your own site.   * You can create a "web application" from one of a couple dozen template that gives you a few more options for customization such as popups and even editing widgets. You can run it on ArcGIS online. You can also download the template, customize if desired and run on your server.   * You can upload a reference to a "web application", basically a bit of publicity and add it to your to your gallery. However, the web application only runs on your own webserver. Question: Can you upload your custom web mapping application to ArcGIS online and have them host it? I'm understanding that this is not possible. Can someone confirm or deny? Thanks!


title: Generalize lines with specific points

text: It's been a while since I was here. I have this problem now. I have electricity line feature class and point feature class intersecting line features. What I need to do is simplify this line feature class but effects must not effect those vertices which intersect with points. Using ArcGIS 10.2.2. Any ideas?


title: ArcGIS 10.2 Copy Rows error

text: I'm very confused by this one. I have run a polygon neighbors analysis in ArcGIS 10.2 and appended some data from another table to help me whittle down the neighbors list before my next action. So I have the following:               bn = "BLPUs_PolygonNeighbors"     arcpy.MakeTableView_management(bn, "Neighborscut")     arcpy.SelectLayerByAttribute_management("Neighborscut", "NEW_SELECTION", where2)     arcpy.SelectLayerByAttribute_management("Neighborscut", "SUBSET_SELECTION", where1)     arcpy.SelectLayerByAttribute_management("Neighborscut", "SUBSET_SELECTION", where3)     arcpy.Delete_management(bn)     arcpy.CopyRows_management("Neighborscut", "BLPUs_PolygonNeighbors")      But when I run this, it says ERROR 000733: Output Table: Same as input C:\Users\Ant\Documents\ArcGIS\ITN.gdb\BLPUs_PolygonNeighbors Failed to execute (CopyRows). Which seems to suggest it thinks "Neighborscut" is actually also called "BLPUs_PolygonNeighbors" Any ideas??


title: postgis, st_within, and a join

text: So I'm in the middle of learning postgis/postgres. I've a polygon table. I've a point table. The polygon layer contains a foreign key that points to another table. For the example the Foreign key points to the type of polygon. So TypeFK in the polygon table references TypeID in the type_tbl; I can do a: select type from polygon as a, type_tbl as b where a.typefk = b.typeid; Now I know the polygon is a type of tree. I can do a: select typefk from polygon join points on ST_Contains(polygon.geom, points.geom); I know what typefk they fall in. I'm ultimately trying to return the type of polygon the points are located within. So - gurus of postgis/SQL - how do I combine the two? I've been sitting here too long without a break unfortunately and am just not getting it. Thanks much for any and all help.


title: How to update new field from existing field?

text: I need to update my new field/column (Property Type) with either "detached", "semi", or "terrace" depending on an existing column (Property Count). So... where the count = 1 I need to populate Property Type with "detached" where the count = 2 I need to populate Property Type with "semi" where the count >=3 I need to populate Property Type with "terrace" How would I write this in the field calculator? I've tried various methods without success.


title: What is the best conversion tools?

text: I need to obtain a lot of base map in the following formats from a PostgreSQL DB:   * *.shp   * *.tab   * *.js (GEOJSON) I know that ogr2ogr is able to obtain these formats type. Do you know of other tools to achieve this ? (a tool usable with some comand lines) Thanks for all your responses,


title: Set ZIndex for Google Fusion Tables in Google Maps API

text: Is there any way to set the ZIndex of a Google Fusion Table layer in the Google Maps API? I have several Fusion layers added as overlays to a map with the ability to turn them on and off. I want certain layers to always be ontop if they are on, however Fusion Tables/Google Maps just renders them in the order that they are drawn. For example: Layer 1 starts On. Layer 2 starts Off. Layer 2 should Always be on top if it is being shown. Click Layer 2 On, and it is on top.   Now, click Layer 1 Off and the On again, and now it is on top as it is the most recent layer drawn. I'd ideally like to set the ZIndex of Layer 2 to higher than Layer 1 so that it is always on top, but there doesn't seem to be a way to do this. Any ideas?


title: How to iterate through features in a shapefile?

text: I am interested in iterating through a point shapefile to perform tasks on a point by point basis using R. In Python, this can be accomplished using the following workflow:               import arcpy          shp = r'C:\path\to\shapefile.shp'          with arcpy.da.SearchCursor(shp, ["Field", "SHAPE@"]) as cursor:         for row in cursor:             print row[0]             # Perform some action with row[1] (i.e. the point geometry)      How can I iterate through a point shapefile to access individual point geometry using R?


title: How to write Shapely geometries to shapefiles?

text: Can someone demonstrate a simple way to write geometry data-structures from shapely into shapefiles? I am particularly interested in polygons with holes and linestrings. It would also be beneficial to stay away from arcpy (so osgeo, pyshp, etc. would all be better).


title: How do I symbol numbers into groups?

text: I am working on creating soils map and need to group symbolize into greater numbers and within between numbers.. Here is an example of what I am talking about ... ![enter image description here](http://i.stack.imgur.com/hVdSR.jpg)


title: Global Geo Api places

text: I'm looking for a global geoapi for the places. I see Facebook places, google places, simplegeo, factual etc... Why not create an API that includes all places? Free and easily accessible as OpenStreetMap for the webmapping ? If people are interested, write me! Thanks


title: How to georeference using GDAL tools?

text: I have a .tiff image that I am trying to overlay in google earth. For that I require to georeference it first. The image projection is UTM with NAD83 as datum. I know the lat/long of the four corners of the image. I am doing the following way:   1. gdalinfo to get the four coordinates.   2. gdal_translate to specify lat/long info   3. gdal2tiles.py to tile the image. The only thing is bugging me is that I never specified UTM and NAD83, then how come it can overlay my image? Does anybody know how to specify this projection information?


title: Some Help to select some features from a vector layer

text: I have a vector layer (kml) that has some point with scientific name of plants and animals, some description of the place where the people saw them and the coordinates. Basically the table associated with the kml is something like this: scientific name//description//lat//long I'm working with OpenLayers, and I really appreciate if you can help to figure out about how to write the correct code: application needed: When I select a feature (from vector 3), All other features that share the same scientific name are highlighted too. Here is my html code: ///////////////////////////////////////////////                                    init();


title: What are the differences between the various ArcGIS Server options?

text: This document makes me want to stab my eyes with a fork. Can someone provide a simple summary of the differences between Enterprise vs. Workgroup, and Basic vs. Standard vs. Advanced? As a bonus, what can I expect to spend on the different options?


title: Access Feature Classes on Multiple Database Schemas with ArcSDE

text: I am using ArcObjects to access a GeoDatabase that was set up with ArcSDE 10 and SQL Server 2008 R2. In this particular case, I have some tables set up on alternate schemas (this is to replicate a clients' environment). The OpenFeatureClass method on the IFeatureWorkspace is working correctly when I access feature classes that are in the default schema (SDE). That same method will fail, however, when I try to access feature classes on the other schemas. The error I'm getting is "Error HRESULT E_FAIL has been returned from a call to a COM component." Has anyone had a similar issue? Is there something I should do to ensure the other feature classes are accessible to ArcObjects? _**UPDATE_** Brad Nesom made a good point in the comments about whether the alternate schema tables are spatial. The table WERE spatial, but for some reason after I changed the schema, they were changed to non-spatial. It's possible that ArcSDE does not allow you to modify schemas for spatial tables.


title: Unable to view cached map service unless logged in as admin in GIS Server

text: I recently re-created a cache for a map service due to SDE server changes. When the cache finished I was unable to view it unless I manually went to the cache folder and viewed each tile one by one in photoshop. After a little research I discovered the virtual directory in the cache directory was incorrect. After changing the virtual directory I can now view the cached service in ArcCatalog through my admin login, but cannot view it through the "Use GIS Services". This happens on both the internal and external sites I have. I've checked security/permissions on everything I can think of, folders, groups, etc. I'm almost certain this is an easy fix, but I'm coming up with nothing. Any help would relieve some stress! Thank you in advance.


title: Creating a Map with Spatial filtered layer on the fly

text: I have the following requirement. I created a web application with User accounts. based on the User type I need to filter certain point layer with a selected cities region. Ie , when user log in, his accessible cities (city ids are retrieved from a db) will get. With these city IDs , select corresponding city geometries regions from the map. After doing an Union operation on all these polygons (city layer) , give resulting final polygon to a spatial Query filter , so that certain point layers need to be shown filtered by this polygon. So that user can only see the his accessible points features filtered by cities regions. I am using ArcGIS JS API V.3.3 with ArcGIS Server 10 Kindly please help me how can achieve this goal. Is a solution possible via REST API/ JS API? Geoprocessing solution also accepted. Thanks in advance Shabeerali


title: How to replicate QGIS vector smoothing effect in ArcGIS for Desktop?

text: You can see from the screenshots that QGIS has, in my opinion, much more pleasant on-screen rendering of vector features. I assume there is some sort of smoothing algorithm used in QGIS. Is there a way to replicate this on- screen smoothing effect in ArcGIS? (click images for full effect) **ArcGIS for Desktop ArcMap 10.2** ![enter image description here](http://i.stack.imgur.com/H1ZYZ.jpg) **QGIS 2.2 Valmiera** ![enter image description here](http://i.stack.imgur.com/y2eGh.jpg)


title: How can I select vector features that overlay raster in arcpy?

text: I have two datasets: point feature class and raster (stored in file geodatabase). The spatial extent of the raster is smaller than the point dataset. Is there a way I could select (or create new, temporary layer/feature class) only the points that fall into the extent of the raster? Basically I'd like to have an output similar to `Clip` or `Intersect` tools, but they don't seem to accept rasters as inputs. Is there readily available arcpy tool(s) that could be used to achieve that? Or would it require some custom solution?


title: I cannot draw a line with a function

text: I'm working with GeoExt and Openlayers and HTML, and I need to draw a line when the user clicks an HTML button. I created a Function Prueba, but nothing happens. If I want to add a control with my function, it works perfectly. Thanks for your help. Here is my js code:               Ext.Loader.setConfig({         enabled: true,         disableCaching: false,         paths: {             GeoExt:"GeoExt"         }     });          Ext.require([        'GeoExt.panel.Map'     ]);                Ext.onReady(function() {         var me = this,layers = [];           var apiKey =   "AqTGBsziZHIJYYxgivLBf0hVdrAk9mWO5cQcb8Yux8sW5M8c8opEC2lZqKR1ZZXf";                  //Mapa OpenLayer             map = new OpenLayers.Map('map');                  //Mapas OpenLayer WMS:Web Map Service////             var Wms = new OpenLayers.Layer.WMS("OpenLayers WMS", "http://vmap0.tiles.osgeo.org/wms/vmap0",{layers: 'basic'});             var OpnSt= new OpenLayers.Layer.OSM.Mapnik("Mapa OpenStreet");                  //Mapas de Google/////             var gphy = new OpenLayers.Layer.Google(                 "Google Fisico",                 {type: google.maps.MapTypeId.TERRAIN}             );             var gmap = new OpenLayers.Layer.Google(                 "Google Calles",                 {numZoomLevels: 20}             );             var ghyb = new OpenLayers.Layer.Google(                 "Google Hybrido",                 {type: google.maps.MapTypeId.HYBRID, numZoomLevels: 20}             );             var gsat = new OpenLayers.Layer.Google(                 "Google Satelite",                 {type: google.maps.MapTypeId.SATELLITE, numZoomLevels: 22}             );             //Mapa Bing             var bing =new OpenLayers.Layer.Bing({                 key: apiKey,                 type: "AerialWithLabels",                 name: "Bing Aereo + Texto",                 transitionEffect: 'resize'             });                  var bingRoad = new OpenLayers.Layer.Bing({                 key: apiKey,                 type: "Road",                 // custom metadata parameter to request the new map style - only useful                 // before May 1st, 2011                 metadataParams: {                     mapVersion: "v1"                 },                 name: "Bing Caminos",                 transitionEffect: 'resize'             });                    //Creamos el Panel Switcher          map.addControl(new OpenLayers.Control.LayerSwitcher());         map.addControl(new OpenLayers.Control.ScaleLine());         map.addControl(new OpenLayers.Control.MousePosition({prefix:'Longitud | Latitud: ', separator: ' | '}));         //Agregamos los tipos de Mapas         map.addLayers([OpnSt,gmap,gphy,ghyb,gsat,bing,bingRoad]);              var mapPanel = Ext.create("GeoExt.panel.Map", {             renderTo: "DivMap",             layout: 'fit',             border:false,             width: "100%",             height: "100%",             map: map,             zoom: 2              });     });                    function Prueba(){         var vector = new OpenLayers.Layer.Vector("Lineas");         map.addLayer(vector);         var line= new OpenLayers.Control.DrawFeature(vector, OpenLayers.Handler.Path);         map.addControl(line);     }         and here is my HTML code when user clicks:               <input type="checkbox"  value="line" onclick="Prueba()" name="type"/>


title: How to make OpenLayers.Popup draggable?

text: I've used jquery-ui draggable option as following but the moving popup can not be ended once dragging has started. Is there any another method to make OpenLayers.Popup draggable?               popup = new OpenLayers.Popup.AnchoredBubble(                     "JInfo",                                            // id                      jFeature.geometry.getBounds().getCenterLonLat(),  // lonlat                      new OpenLayers.Size(ISA.Size.POPUP_WIDTH, ISA.Size.POPUP_HEIGHT),                      html,                                           // contentHTML                      null,                                           // anchor                      true,                                           // closeBox                      function(){                                    // closeBoxCallback                               jSelectControl.unselect(jFeature);                      });             jFeature.popup = popup;         map.addPopup(popup);         $("#JInfo").draggable();      Thanks advance for your helps, Yasemin


title: how to stop ArcGIS 10.1 automatically creating mapped drives?

text: I have a MXD which I am using to make a PMF. I'm using ArcMap 10.1. I want my layers only to have UNC paths (starting with \\server...) as otherwise some users will get broken links if we have inconsistent mapped drives. Unfortunately, sometimes ArcMap 10.1 switches some layers to my mapped drive equivalent of the file path. I don't know why. When I try to change them all back to UNC paths in a wholesale fashion by disconnecting the relevant mapped drive so I can then switch one layer to the UNC path and let that propogate through all, ArcMap 10.1 decides instead to re establish the relevant mapped drive when starting up. Anyone know a fix?


title: Best Practices for Tracking Feature Changes and Supporting Documentation

text: In my organization, it is very important that we have documentation from an authoritative source to support every feature that we add to our data and every edit that we make to each feature. Currently, our organization keeps an excel file on which we track every change we make and include links and file paths to the documentation that supports the change. Lately, these excel spreadsheets have become a bit cumbersome. I am curious if anyone else deals with this level of documentation, and how they do it. I know about metadata and the track edits tool, but neither seems to provide the level of detail that we need to maintain data lineage and documentation for every single feature. Is there a tool or a preferred method for dynamically linking each feature in a feature class to detailed information about it's lineage and documentation? Thanks!


title: How to save (Commit) a temporary table in Mapinfo?

text: I tried the following:               Commit Table Query1 As "C:\Temp\XYZ.TAB" TYPE Query Interactive      It does not work.


title: How to: Conditional Symbology in ArcMap

text: I'm in charge of a database that contains property information. I would like to symbolize parcels of land that we have a 'written agreement' with the owner differently than parcels that we do not have written agreements with. My solution was to create a calculated field (1 or 0) based on if there was any data in the "Agreements" field. > I'm looking for a more dynamic solution (not having to re-run the field > calculation) so that when an agreement is entered into the db, the symbology > automatically changes. This is important because our field staff take the db out into the field on hand-held devices; they need to know which **properties they are allowed to access (green)** vs. which **properties they cannot access (red).**


title: Literature on bias in spatial weights matrix

text: I have a question with regard to a possible bias as a result of how one constructs the spatial weights matrix. Consider a case where one works with areal data, for instance we're interested in whether the crime rate in one neighbourhood depends on the rate in other states. The spatial weights matrix is this case can be constructed using contiguity,e.g. queen contiguity based on 1st order neighbours, or using a distance cut off. Working with areal data the polygons used as unit of analysis are rarely equally size. So it seems to me that if one uses the distance based spatial weights matrix the results might be biased due to this variability in size among the polygons. I wondered whether anyone knows of any good references from the literature on this subject? I didn't have much success with my internet search so far. Thanks in advance.


title: Pictometry Image Control?

text: Does anyone have any references for this? (Is it OK to ask about it here?) Specifically the ImageControl API. I'm tasked with writing a program using it.


title: How do I serve MBTiles with TileStache?

text: I'm am looking for a way to serve tiles while offline. I'm using Windows 7. Here's what I've done so far:               -Uploaded a portion of OSM data into TileMill     -Converted the OSM data into a file called osm.mbtiles     -Installed Python and Mapnik     -Installed TileStache (to C:\TileStache-1.49.10) and its dependencies (Pillow,   ModestMaps,SimpleJSON, Werkzeug)     -Placed osm.mbtiles into the TileStache folder (to C:\TileStache-1.49.10\osm.mbtiles)     -Written test.cfg to look like this:          {       "cache":       {         "name": "Test"       },       "layers":        {         "osm":         {             "provider": {"name": "mbtiles", "tileset": "osm.mbtiles"}              }       }     }          -Place the config file here: C:\TileStache-1.49.10\test.cfg      Then, after all of that is done, I run:               C:\TileStache-1.49.10>tilestache-server.py -c test.cfg      After I run that command, I see:                * Running on http:<IP_ADDR>:8080      So, I open up my web browser and go to:               localhost:8080/osm/preview.html      which is then routed to:               http://localhost:8080/roads/preview.html#10/37.8000/-122.2600      I see nothing on this page, but the server in the terminal is spitting out this: ![enter image description here](http://i.stack.imgur.com/TRgBL.png) Also, I've tried all the steps above using world map converted to .mbtiles using zoom level 0-8 and then manually changing the URL of the tilestache server page to be:               http://localhost:8080/roads/preview.html#3/50/50      But, I still get the same errors and nothing shows up. What should I try to fix this?


title: OpenLayers Enable KML feature based on user form

text: I have a KML layer with over 300 features, but I only want to display features based on a user form. Each feature has a unique ID and Sector and I have added the entire KML to the map but have the visibility set to false. How would I go about adding only the features to the map that contain the unique ID and Sector values specified after a user hits submit on a form?


title: How to automatically merge Texas survey block data

text: I have Data from the Texas rail road commission. I'm trying to get the block data to display with a thick line. The problem is the data is in surveys. I exported a copy of the shapefile. If you take a look at the L2BLOCK colunm you see that I started merging them together, but it take forever if you have to do so for a lot of counties. ![enter image description here](http://i.stack.imgur.com/tPWCF.png) I want the blocks to have a thicker line to seperate them from the surveys. Here's what the table above looks like. ![enter image description here](http://i.stack.imgur.com/OVhbJ.png) You can see that some of them are merged and some not. Is there way to automate merges based on a same block numbers in the attribute table?


title: Redraw info window after centering map on click (ESRI Javascript API)

text: When I click on one of the features on my map the map gets re-centered using the following event handler:               var pointClick = projects.on("click", function(e){         // Center map on selected feature         map.centerAt(e.mapPoint);              // A bunch of other stuff     });      Unfortunately the info window describing that point does not move to the new location unless the map gets redrawn (by moving to another zoom level, for example). Is there a method to force the info window to refresh? The information in the info window reflects the point that was most recently clicked on, but the location of the info window does not change. Screencast: http://www.screencast.com/t/fo1Nec0CBzk6 I can provide the code if necessary.


title: QGIS is not able to find GRASS. Choose the path of installation folder of GRASS

text: I've downloaded and then correctly installed `QGIS 1.7.4` for windows (Win 7 64bit is my version) in `C:\Programs\Quantum GIS`. But everytime I start QGIS this popup appear: `"QGIS is not able to find the installation of GRASS. Choose the path of the GRASS installation folder".` I've chosen every installation folder in Quantum GIS but the program not recognise it! So if I don't choose the folder, I can't use GRASS. In the installation folder there are a file called `"grass64.bat"`, but not `"grass64.exe"`. Could it be the problem? I've installed QGIS also in other folder and `D:\Programs\\...` but without results.


title: Different level of detail(layers) on different zoom level at  OpenLayers map

text: I want display different detail(layers) on different zoom level on OpenLayers map. I mean when zoom in should display lot of layers and when zoom out those unsuitable layers should hide automatically. Different layers should display on different zoom level ,not using layer switcher. Can this do using OpenLayers? If not what is the open source solution? If yes, pls give some example or links to get some idea Thanks for reading and answering


title: Is it possible to use assign_vertex_id() POLYGONs?

text: I'm trying to use pgrouting in order to create a graphe and I have two questions : **Is it possible to use the function assign_vertex_id() with a POLYGON geometry type?** I always have the following error when I try to use this fonction :               psql (9.0.3)     Type "help" for help.          bpi=# SELECT assign_vertex_id('rangee', 0.00001, 'geometrie', 'id');     NOTICE:  CREATE TABLE will create implicit sequence "vertices_tmp_id_seq" for serial column "vertices_tmp.id"          CONTEXT:  SQL statement "CREATE TABLE vertices_tmp (id serial)"          PL/pgSQL function "assign_vertex_id" line 14 at EXECUTE statement          ERROR:  query string argument of EXECUTE is null     CONTEXT:  PL/pgSQL function "assign_vertex_id" line 31 at EXECUTE statement      the table rangee is :               bpi=# \d rangee                    Table "public.rangee"              Column    |         Type          | Modifiers      -------------+-----------------------+-----------       id          | integer               | not null       zone_id     | integer               |       geometrie   | geometry              |       description | character varying(50) |       source      | integer               |       target      | integer               |          Indexes:      "rangee_pkey" PRIMARY KEY, btree (id)      "rangee_spatial" gist (geometrie)          Triggers:          rangee_geom_srid BEFORE INSERT OR UPDATE ON rangee FOR EACH ROW EXECUTE PROCEDURE geometrie_srid()      I don't understand this mistake. Could you help me please? My second question is about pgrouting. **Is it possible to create a dual graph with this tool?** Thanks a lot and I'm sorry for my english.


title: Reprojecting a shape to UTM33

text: I still have trouble reprojecting shape files :( I would like to reproject it to UTM33N / WGS84 (EPSG 32633) I think the projection of the shape might be Pulkovo 1984/ Gauss-Krüger Zone 4 (EPSG 2398) or ETRS89 / UTM33N (EPSG 25833) http://projektidee.org/images/extern/ffh.zip I failed to reproject it though... Could anyone here help me with that? best regards and thanks desputin **[this scackexchange is getting on my nerves!!!! I can't even answer to your questions so I have to edit this fist thread!!]** Yes sorry. I'm trying to do this with Qgis. What I do is right click on layer an save as - then I choose the Projection system UTM33N/WGS84 and save the file to a new shape. The other thing was to first assign the other projection (pulkovo) to the shape layer and then save it to a new shape file. With both ways of dealing with it I didn't succeed...


title: Extracting point data from Digital Raster Graphic (DRG) in ArcMap?

text: I am trying to run an IDW analysis based off stream and springs data. I have the hydrography.eoo file from GeoComm but am now wanting to create a point layer of the springs in the area which show up on the DRG. Is there any way to extract point data from a Digital Raster Graphic (DRG) in ArcMap without manually creating a spring layer? I'm using ArcGIS 10.1 for Desktop.


title: Opening ArcGIS layer package in open source software?

text: It seems that ESRI introduced the layer package in ArcGIS 9.3.1. I have been provided with layer packages containing essential data, but do not have access to this software. Ideally, I'd like to be able to open them using some open source tool. I mainly use R, but Quantum GIS or GRASS would be okay at a stretch. Suggestions about the best way to access this data are appreciated.


title: PostGIS - 1.5 to 2.0 easy update

text: I have a PostGIS 1.5 database and want to update to 2.0. I read that because this is major release I have to follow the 'Hard update' procedure. However the spatial data I have is only test data so I can drop it. If this is the case, is there easier procedure to update to 2.0 if I don't care about spatial data (I do care about non-spatial tables)?


title: Adding arbitrary code to ArcGIS ModelBuilder?

text: Is there any way to add arbitrary code snippets to an ArcGIS model? I envision something similar to the "show code block" the expression builder in advanced labelling or field calculator. For example my project of the moment could benefit greatly from fetching the username of the person running the tool. I know in python this is a simple `import os; os.environ.get('USERNAME')`. I think I could write that snippet to `get_username.py`, add it to a toolbox in Arccatalog, then add that tool to the model and use it as per-requisite model parameter for my current tool. That's a lot of work for such a small thing, and next time I might want something else, like inline variable substitution or a compound variable. I think a code snippet tool or similar could be address all of these needs in a generic, broadly helpful manner.


title: How to start an ArcGIS Server 10.0 cached MapService without database connections?

text: I have a cached ArcGIS 10.0 MapService that I need to start but the DB which was used as the datasource for creating the tiles is no longer available. Is there a way to start the service without datasources?


title: Draw graphic on screen without requiring a double click to finish sketch

text: I am trying to get the user to draw a line(only 2 vertices) on the screen(just like the georeferencing link/spatial adjustment link). I got the below snippet from ESRI website which does give me a polyline but it requires the user to double click to finish the sketch.               public IPolyline GetPolylineFromMouseClicks(IActiveView activeView)             {                    IScreenDisplay screenDisplay = activeView.ScreenDisplay;                    IRubberBand rubberBand = new RubberLineClass();               IGeometry geometry = rubberBand.TrackNew(screenDisplay, null);                    IPolyline polyline = (IPolyline)geometry;                    return polyline;                  }      I would like the tool to finish the drawing as soon as two points are clicked, thus avoiding a double click. Can somebody help me on this? I seem to be making a mess of a simple situation!


title: How to display indoor map using OpenLayers?

text: I created my osm file in JOSM, and i imported it to my database by osm2pgsql tool. I created new layer in my osm.xml file, which selecting data with query:               select way from planet_osm_line,     (select parts from planet_osm_rels WHERE (((hstore(tags)->'level'))  = '0')) p     where osm_id = any (p.parts))       when i run my python script to generate png file i got an error: > RuntimeError: Postgis Plugin: PostGIS: unable to query the layer extent of > table 'planet_osm_rels' because we cannot determine the geometry field name. > Please provide either an 'extent' parameter to skip this query, a > 'geometry_field' and/or 'geometry_table' parameter, or add a record to the > 'geometry_columns' for your table. is any one can tell me what i'm doing wrong and how i can fix that?


title: How to create tiles with points from MongoDB?

text: I have a (potentially) large data set stored in MongoDB. Each record is just `lat`, `lon` and some float-point `value`. Each pair of `lat` and `lon` is unique, but the whole set is not covering the Earth, there may be regions with no data at all (so they're not stored). I need to draw my own tiles based on this data, where mentioned `value` displayed with respective color (say, from black for minimal values to white for maximal). My dataset is changing in time, so I have to re-render tiles periodically. **I've installed GDAL but I can't find a way to fetch the data.** I've tried to export points to XYZ format, but got no luck: `gdal2tiles.py` doesn't want to work with it (I think the cause is that my coords are not sequential: I found a sample XYZ data set in Internet, and it was converted OK, and the only difference is that file had a value for each node in grid). Maybe it will work if I populate my database with some initial values for each possible coords pair, but it will cost much time, disk space and DB performance, so it's the last resort. How can I make this stuff work? What am I doing wrong? Is there something to read that will help me?


title: Create a generic shape and draw it as a polygon with center (or base) and rotation?

text: In Java (jts), is there a way to draw a generic shape, eg. bounded in meters, and apply this to a coordinate, with rotation around a center or base? Something like:               Box box = new Box(10, 100, METERS);     Rectangle rect =          GeometricShapeFactory.createRectangle(box, 65.100, 7.123, Math.toRadians(45));


title: Syntax error: keyword can't be an expression with the function "Con"?

text: From How can I reclassify an integer raster to a float raster?, I would like to reclassify an integer raster "cpraster3.tif" to a float raster. For example, I would like to reclassify "Value" = 1 into 0.27, "Value" = 2 into 1.00 ... Here is my code:               test_reclass = Float(Con(("cpraster3.tif", "Value" = 1, 0.27, Con("cpraster3.tif", "Value" = 2, 1.00, Con("cpraster3.tif", "Value" = 3, 2.9, Con("cpraster3.tif", "Value" = 4, 2.94, Con("cpraster3.tif", "Value" = 5, 5.37, Con("cpraster3.tif", "Value" = 6, 10.35, Con("cpraster3.tif", "Value" = 7, 16.49, Con("cpraster3.tif","Value" = 7, 22.08, 0.0)))))))))      But I have this error message: > Failed to run script - syntax error - keyword can't be an expression


title: add a Custom search to an openlayers map

text: I have looked at this thread: add location search to openlayers map But it is not the result that I am trying to achieve. Probably best if i show you: http://www.potholestudios.com/Content/pid=46.html I have made this Interactive Game Map using Open Layers. What i would like to do is add a search Box, so when say someone types into a Search Box "Deep Space 9" the map would fly to coordinates (-33.87, -8.58) Deep Space 9 in the Beta Ursae sector block. Now, I'm guessing I'd have to list the X,Y coordinates of every system in a list with Names. So that when a Name is Entered in a Search Box it moves the map to the listed coordinates. But i have no idea even where to start with this. any help would be appreciated.


title: Seeking free nautical route planner

text: Not actually a route planner, but it amounts to the same thing. I want to know the distance a vessel at sea must travel between two points. Obviously this isn't just a straight line thing, if there is land on that straight line. I would prefer a solution with google maps, but anything will do really, just so long as it tells me how far the vessel must travel to get from point A to point B (it doesn't need to tell me how to get there, that would be a bonus). Thanks in advance for any help. @Geoist +1 "You want something that takes into account most efficient routes given currents, tides, etc?". Yes, please, if that is possible. If it costs extra, then I would prefer to start with a free solution and add that later.


title: How to get The character spacing from a ITextElement arcgis engine 9.3 C#

text: i have a problem.I can set character spacing to a IElement as ITexElement like :               TextSymbolClass pTextSymbol = new TextSymbolClass();     pTextSymbol.CharacterSpacing = 10 ; //etc     ((ITextElement)firstElement).Symbol= pTextSymbol;     //firstElement is a Element in map      Now i want to get value of CharacterSpacing when i focus to a ITextElement in Map.But when i used this code below, some thing happen and it break from function and have no exception.               TextSymbolClass pTextSymbol2 = new TextSymbolClass();         pTextSymbol2 = (TextSymbolClass)((ITextElement)firstElement).Symbol; //break from here         value=pTextSymbol2.CharacterSpacing ;      So, how can I get the CharacterSpacing value of a Element,can you show me another way to get it.I'm using Arcgis engine 9.3, C# .Thank you so much.