title: How to setup users and user groups on a distributed architecture?

text: I am learning ArcGIS distributed installation. I planned like this picture. ![enter image description here](http://i.stack.imgur.com/Ayeil.png) I read from installation documents about server users group.   1. agsAdmin   2. agsUsers Which machine contains user groups and users? Web server, SOM or SOC PCs? For example I installed Windows 2008 on PC-1, 2, and 3 machines, all of these are default accounts. I installed PC-1 as SOM. PC-2 as SOC machine. Where my usergroups should be located or all pcs has these groups.


title: How to create a visualisation of the world's aeronautical flight paths?

text: I am a sculptor making a world globe art project, and I need a still image of the whole world’s actual air traffic routes to use as reference material. To be clear, I’m not looking for the ‘air traffic flow’ visualisation, which consists of straight line routes interpolated between the airports. I’m more interested in seeing the structure of the road network in the sky, revealed through the actual movements of the planes. The visualisation im describing appears to have been performed in localised areas, but not yet over the whole globe. I found this visualisation by Aaron Koblin of the air traffic routes over the USA. http://www.aaronkoblin.com/work/flightpatterns/index.html ![USA air traffic](http://www.aaronkoblin.com/work/flightpatterns/1_1920.png) This is basically what I’m looking for, except done over the whole world. And I came across this animation by 'ITO World' of the European air traffic coming back to life after being grounded by a volcanic ash cloud. http://www.itoworld.com/static/airspace_rebooted.html ![europe air traffic](http://www.itoworld.com/static//i/visualisation/ash/ash_europe_x2.jpg) I also found these KML air navigation plans drawn by ICAO:               http://192.206.28.84/Website/ANP%20AERONAUTICAL%20INFORMATION.htm      Scroll down three-quarters of the way until you get to this image: ![air navigation plan](http://i.stack.imgur.com/e2eH6.jpg) ...and then click on the links “MID route UPPER operational" for the mid Atlantic, and the 6 similar links below for KMLs of other parts of the world.   They are useful to me, but I find them confusing in places where lines don’t seem to make sense (ending abruptly and not starting again), and they do not cover the open ocean, or show which routes are used most heavily. And here is http://www.flightradar24.com/. This shows quite a lot of the world’s air traffic movement live. Its receivers do not have full worldwide coverage, and many types of plane do not show up, but as far as I know it’s the most complete coverage out there. I imagine the best solution for getting what I want would be to make a time lapse image of the flightradar24.com data over quite a long period of time – a month or maybe even a year to even out seasonal fluctuations, which might be quite pronounced in places like the arctic. However, if a years worth of data is unmanageable, then I'll have to go for less. As I am using the final digital image as reference material for something I am rendering by hand, I will probably be able to infer any patches that are missing - because of incomplete receiver coverage - using the ICAO navigation plans. So, assuming it is possible, I would greatly appreciate it if somebody could **take me through the process of gathering the necessary data from flightradar24.com and creating the visualisation that I need.** Please go into detail, as I am quite new to GIS, and am not particularly knowledgeable of computers. I am running QGIS on a 4GB RAM PC. If this is the wrong program to be using, or if my computer is too small for this task, please advise me on what I would need. I am hoping to make a large GeoTiff (8 tiles of 13500 x 13500 pixels, a similar size to NASA’s Black Marble 2012 image would be good), showing a map of the worlds flight paths. The most heavily used paths would show up darkest and no traffic at all would be transparent. That way I can use it as an overlay on a world map so I can get my barings. Of course, if any of you find the idea of performing this visualisation fascinating, especially as it doesn’t seem to have been done over the whole world before, feel free to do it and publish the results! Thank You   Martin


I have successfully published a Geoprocessing Tool as GP services, I had called it with success, but when I try to call it, now it gives me this error `GPServer does not exist or is inaccessible` for no reason. MapServer and GeocodeServer work fine I am having a problem just with the GPServer. Any help please.


title: Python: Arcpy: Adding data to a shapefile

text: Complete arcpy beginner here... I have a shapefile with the fields `CNTRY_NAME` and `GDP_HEAD`. I have a list in Python `gdppop2002 = [['Country1, 234.1232341], [Country2, 23.2341234],...etc...]` I'm trying to write each float value in the list to the corresponding country row in the shapefile's attribute table. I need to check if each Country in `gdppop2002` matches a `CNTRY_NAME`, and then write the float value in the list to `GDP_HEAP`. From a couple of hours of googling this appears to require creating an insert cursor, but I'm a little lost. Would anyone be able to point me in the right direction?


title: SELECT assign_vertex_id('ways', 0.00001, 'the_geom', 'gid'); not work

text: I'm using using OSgeo 6.5 dvd and studing pgrouting workshop (link) for get shortest path. I created network topology call "pgrouting - workshop". There are three commands to execute in section 5.2 they are :               -- Add "source" and "target" column     ALTER TABLE ways ADD COLUMN "source" integer;     ALTER TABLE ways ADD COLUMN "target" integer;          -- Run topology function     SELECT assign_vertex_id('ways', 0.00001, 'the_geom', 'gid');      Add source and target was executed successfully. topology function give errors. Heres the error :               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 15 at EXECUTE statement     ERROR:  function startpoint(geometry) does not exist     LINE 1: SELECT gid AS id, StartPoint(the_geom) AS source, EndPoint(t...                           ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.     QUERY:  SELECT gid AS id, StartPoint(the_geom) AS source, EndPoint(the_geom) as target FROM ways     CONTEXT:  PL/pgSQL function "assign_vertex_id" line 27 at FOR over EXECUTE statement      How do i fix this error? I loaded the legacy.sql and noted that error has little bit changed :               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 15 at EXECUTE statement     ERROR:  function point_to_id(geometry, double precision) does not exist     LINE 1: SELECT point_to_id(setsrid(_r.source, srid), tolerance)                ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.     QUERY:  SELECT point_to_id(setsrid(_r.source, srid), tolerance)     CONTEXT:  PL/pgSQL function "assign_vertex_id" line 33 at assignment


title: Gathering data for a large number of census block groups

text: I would like to collect some census data for every block group in New England. I'm primarily interested in median income and median age, and ideally I'd like a large .csv file with one row per census block group and one column per variable (e.g. id variables, median_income, etc.). American FactFinder seems almost useless for this kind of thing. Does the census have a "bulk data export" option?


title: Finding AniMove in Qgis 2.x

text: I've just upgraded to Valmiera in QGis, the only way to do spatial analyses appears to be with the experimental AniMove plugin (?) which is not appearing in my GUI. Can anyone tell me where to find the plugin and how to make it work? Many thanks Zeke


title: I'm new to GPS, and I had a quick question about clarity.

text: I'm working on a project for class, and I'm mapping GPS on Google Earth with a KMZ file. I used a handheld GPS 24°6.691'N + 74°27.833'W. How accurate is my coordinates? I know I'm in decimal minutes! That's about it.


title: How to project geographic coordinates to cartesian with the least possible error?

text: I need to convert geographic coordinates (in decimal degrees) to cartesian with the least possible error? I searched a lot on the internet but it seems that the conversion is not exact. Could you please present a **conversion formula** with the least possible error? I am working on a set of data including hospital locations obtained by Google earth. * * * I am working on a mathematical problem and the hospital locations is just an input parameter for the problem. If the error is less that a few meters it is absolutely acceptable. I just need to project the (long,lat) coordinates to a 2-dimensional cartesian coordinates. Is the following formula usually used for the conversion? x = R * cos(lat) * cos(lon) y = R * cos(lat) * sin(lon) I saw that different sources consider different values for R! what value is better? Some other sources say that this formula is not appropriate as the earth is not an exact sphere!


title: OpenLayers 2.13 : How to display a local tiled raster?

text: I am trying to display (on OpenLayers 2.13) a local tiled raster I made with gdal2tiles.py. Unfortunately it's all blank or pink. Here is my layer definition:               var tiles = new OpenLayers.Layer.XYZ("Local tiles", "tiles/${z}/${x}/${y}.png", {         isBaseLayer: true,         displayProjection: new OpenLayers.Projection("EPSG:4326"),         maxExtent: [-75.1729846, 39.9439032, -75.1529085, 39.9599679]     });          map.addLayer(tiles);      What is wrong? Here is a link to my OpenLayers test page : http://popopopopo.voila.net/index.html And here is the original image: http://popopopopo.voila.net/raster.tif (EPSG: 3857)


title: Issues with pgrouting workshop code in postgis 2.11/pgrouting 2.0

text: I'm having issues getting the pgrouting workshop to work on Postgresql 9.2/PostGIS 2.11 and pgrouting 2.0 on windows. How can I rewrite the following php/sql code and make it compatible with my version:               <?php            // Database connection settings       define("PG_DB"  , "routing");       define("PG_HOST", "localhost");        define("PG_USER", "postgres");       define("PG_PORT", "5432");         define("PG_PASSWD",   "*******");        define("TABLE",   "eastlegon");             $counter = $pathlength = 0;            // Retrieve start point       $start = split(' ',$_REQUEST['startpoint']);       $startPoint = array($start[0], $start[1]);            // Retrieve end point       $end = split(' ',$_REQUEST['finalpoint']);       $endPoint = array($end[0], $end[1]);            // Find the nearest edge       $startEdge = findNearestEdge($startPoint);       $endEdge   = findNearestEdge($endPoint);            // FUNCTION findNearestEdge       function findNearestEdge($lonlat) {              // Connect to database         $con = pg_connect("dbname=".PG_DB." host=".PG_HOST." user=".PG_USER." password=".PG_PASSWD);              $sql = "SELECT gid, source, target, the_geom,                   distance(the_geom, GeometryFromText(                       'POINT(".$lonlat[0]." ".$lonlat[1].")', 4326)) AS dist                  FROM ".TABLE."                   WHERE the_geom && setsrid(                       'BOX3D(".($lonlat[0]-200)."                               ".($lonlat[1]-200).",                               ".($lonlat[0]+200)."                               ".($lonlat[1]+200).")'::box3d, 4326)                  ORDER BY dist LIMIT 1";              $query = pg_query($con,$sql);                $edge['gid']      = pg_fetch_result($query, 0, 0);           $edge['source']   = pg_fetch_result($query, 0, 1);           $edge['target']   = pg_fetch_result($query, 0, 2);           $edge['the_geom'] = pg_fetch_result($query, 0, 3);                // Close database connection         pg_close($con);              return $edge;       }            // Select the routing algorithm       switch($_REQUEST['method']) {              case 'SPD' : // Shortest Path Dijkstra                 $sql = "SELECT rt.gid, AsText(rt.the_geom) AS wkt,                         length(rt.the_geom) AS length, ".TABLE.".id                      FROM ".TABLE.",                          (SELECT gid, the_geom                              FROM dijkstra_sp_delta(                                 '".TABLE."',                                 ".$startEdge['source'].",                                 ".$endEdge['target'].",                                 3000)                          ) as rt                      WHERE ".TABLE.".gid=rt.gid;";           break;              case 'SPA' : // Shortest Path A*                 $sql = "SELECT rt.gid, AsText(rt.the_geom) AS wkt,                           length(rt.the_geom) AS length, ".TABLE.".id                        FROM ".TABLE.",                            (SELECT gid, the_geom                                FROM astar_sp_delta(                                   '".TABLE."',                                   ".$startEdge['source'].",                                   ".$endEdge['target'].",                                   3000)                            ) as rt                        WHERE ".TABLE.".gid=rt.gid;";             break;              case 'SPS' : // Shortest Path Shooting*                $sql = "SELECT rt.gid, AsText(rt.the_geom) AS wkt,                           length(rt.the_geom) AS length, ".TABLE.".id                        FROM ".TABLE.",                            (SELECT gid, the_geom                                FROM shootingstar_sp(                                   '".TABLE."',                                   ".$startEdge['gid'].",                                   ".$endEdge['gid'].",                                   3000, 'length', false, false)                            ) as rt                        WHERE ".TABLE.".gid=rt.gid;";           break;              } // close switch            // Database connection and query       $dbcon = pg_connect("dbname=".PG_DB." host=".PG_HOST." user=".PG_USER." password=".PG_PASSWD);            $query = pg_query($dbcon,$sql);             // Return route as XML       $xml  = '<?xml version="1.0" encoding="UTF-8" standalone="yes" ?>'."\n";       $xml .= "<route>\n";            // Add edges to XML file       while($edge=pg_fetch_assoc($query)) {                $pathlength += $edge['length'];              $xml .= "\t<edge id='".++$counter."'>\n";         $xml .= "\t\t<id>".$edge['id']."</id>\n";         $xml .= "\t\t<wkt>".$edge['wkt']."</wkt>\n";         $xml .= "\t\t<length>".round(($pathlength/1000),3)."</length>\n";         $xml .= "\t</edge>\n";       }            $xml .= "</route>\n";            // Close database connection       pg_close($dbcon);            // Return routing result       header('Content-type: text/xml',true);       echo $xml;          ?>      Thanks,


title: ArcPy: Drop fields from list of fields for arcpy.da.SearchCursor

text: I have a feature class with several fields in the attribute table, and I want to export a few but not all of those fields out to a tab delimited text file. The code below sort of works, half of the fields in the uneeded_fields variable are dropped, but three (SHAPE, long, and utm_y) still end up in my text file. I'm not sure why some of the fields are deleted while the other are passed on to the search cursor. I have tried copying the names from the field set up in Catalog and pasting them into my script so it isn't a typo issue. I'm not sure what I'm missing here. As a side note, will the SHAPE field always be the first field returned unless specified otherwise? If so del fields[0] can take care of that, but the other will still be a problem.               fields = [f.name for f in arcpy.ListFields(fc)]          uneeded_fields = ('SHAPE', 'OBJECTID', 'lat', 'long','utm_x', 'utm_y')          for i,f in enumerate(fields):         if f in uneeded_fields:             del fields[i]          with open(outfile, 'w') as f:         f.write('\t'.join(fields)+'\n') #column headers         with arcpy.da.SearchCursor(fc, fields) as cursor:             for row in cursor:                 f.write('\t'.join([str(r) for r in row])+'\n')


title: ogr2ogr converting MS SQL Server spatial view to MapInfo TAB -error in datum and encoding

text: Running this ogr2ogr command:               ogr2ogr -f "MapInfo File" "D:\V_GeoNode.tab" \     "MSSQL:Server=xxx\SPECIFIC;Database=db1;UID=xxx;PWD=xxx" \     -sql "select * from V_GeoNode" -a_srs "EPSG:25832"      It returns a MapInfo TAB from a SQL server view. The problem is that the datum of the created MapInfo table is WGS84 and not ETRS90 as requested. Also the encoding in the TAB files is set to Neutral and not WindowsLatin1 resulting in Danish letters `ÆØÅ` not to be shown correct. Using the GDAL version 1.10.1 distributed with OSGeo4W (QGIS 2.01). Any hints to get the correct datum and encoding with org2ogr?


title: Ignore Transparent Pixels for Bounds Determination

text: I am displaying an `Image` using GWT OpenLayers on a map. My `Image` has transparent pixels and the system which returns my `Image` returns a footprint that is equal to the image plus the transparent pixels. Since the `Image` has been orthorectified, the `Image` does not occupy those transparent pixels. The `Image` is a transparent PNG. I would like to "shrink" the footprint to match the actual image pixels by adjusting the four corners down to the actual image. My footprint is just a `Polygon` with four corners. Since I am not able to change the system that gives me the non-orthorectified four corners, is it possible in OpenLayers to return the corners of the image for the visible pixels so I can "shrink" my footprint polygon to match the image four corners, essentially tricking the map into producing valid output.


title: How to use directional shading on boundaries for cartographic effect?

text: The attached screenshot of a National Geographic map shows some beautiful directional shading of boundaries. I'm trying to reproduce this sort of shading for a burn perimeter map, as I want to highlight features within the burn perimeter and not overshadow them with a transparent polygon or crosshatch. Is there a way to accomplish this sort of boundary shading using ArcGIS 10? Link to larger map titled: "1946 Political Subdivisions Of India By National Geographic.jpg" ![enter image description here](http://i.stack.imgur.com/kfgjW.jpg)


title: arcgis halo inside border style?

text: The following was marked as duplicate of the question here, but I disagree: ## the answers to that question all involve modifying the geometry or making derived geometries, and they deal with shading and transparency rather than fill patterns. The question below clearly says those options have already been explored and do not do the trick. I'm looking for a way that does not involve modifying the geometry or making derived geometries - let's face it, we're talking about a style here, not about ways to modify geometry - and is generic to any fill type or transparency level. If this is still deemed a duplicate question, then the answer must be "no you can't do that as a style". I'm new to arcgis but have been using qgis for a while. Question on polygon styles in arcgis: is there a way to make an 'inside buffer' outline or border style? I know you can use the buffer tool to make an actual feature that is buffered by a number of geographical distance units, but for this case it would be much better to have a style that has a buffer (or halo) of a certain number of pixels regardless of zoom level, something like this quick hackjob in gimp: ![enter image description here](http://i.stack.imgur.com/uJhos.png)   * solid line for the boundary of the polygon   * inside buffer or halo of a crosshatch pattern >>that does not rotate as a function of polygon outline angle at that point<< (I tried a marker line inside halo (offset by a negative number) and it looks fairly dorky due to rotation along the line)   * the rest is transparent Can it be done? Seems like this isn't that unusual of a thing to do. Boundaries of all sorts are drawn this way on paper maps and they look quite nice. Thanks in advance


title: Conversion problem of projection of MODIS LAI from sinusoidal to geographic

text: I am trying to change MODIS data's (LAI) sinusoidal projection to geographic. Actually I did it. But when I reproject it by ERDAS, the output LAI raster moved about 10 km to west side of actual shape file and I tried it by ENVI, the output raster also moved about 20 km but to north side of shape file. I am using the shape file to check accurrancy of the output raster. It was shape of lake. Also I have tried it in R to convert projection using following script               library(sp)     library(rgdal)     library(raster)     Sys.setenv(TZ='GMT')          lai_real <-raster("E:\\LAI\\LAI_hdf\\MOD15A2_LAI_2014\\New Folder\\mod15a2.a2014169.h23_26.img")         projection(lai_real) <- CRS("+init=epsg:4326")           writeRaster(lai_real, filename="E:\\LAI\\LAI_hdf\\MOD15A2_LAI_2014\\New Folder\\mod15a2.a2014169.h23_26_geo.img", datatype='FLT8S', overwrite=TRUE)      But it doesn't work. How can I convert the projection of the LAI from sinusoidal to geographic (lat/long)? If my explanation is unclear, please visit the link that shared below. I hope then you understand what problem I have to convert sinusoidal to Geographic (WGS84) https://www.dropbox.com/s/cucntnt8kh7m1dl/Hovsgol.PNG.


title: QGIS for Android - installation with QGIS-Installer faild on Galaxy Note 10.1 LTE

text: I once installed QGIS for Android. After deinstallation i'm not able to reinstall QGIS on my Galaxy Note 10.1 LTE. I've enabled App installation from unknown sources! Any ideas where the problem might be? Thanks in advance!


title: how to give pop-up for a point type layer in ms4w

text: My aim is to create pin point on map and click event on that pins.I am able to generate map with the pins on the exact locations but i don't know how to give pop-up when clicked on the pin.Help me with this.Thanks in advance below is my code               LAYER     NAME "Siaya_towns"     DATA Siaya_towns     TYPE POINT     STATUS DEFAULT     CLASS       COLOR 255 255 0       OUTLINECOLOR 0 0 0       SYMBOL "pin.png"       SIZE 8       END     END


title: Store global reference to a service/object

text: In developing with ArcMap, is there any way to register a service that can be used during the life-cycle of a document- for example to do something like: C# .NET               mxApplication.RegisterService(myWebService);       then later-               mxApplication.GetService<WebService>().SaySomething("Hi");      So kind of using the Service Locator pattern or similar.


title: How do I project my WMS tilelayer in leaflet if I want it to display as an overlay on my basemap if they have different EPSG?

text: So my basemap is a world Mercator projection and in wfs format and I'm trying to put a WMS overlay in a different projection over it? Leaflet's API didn't give me much info on changing the epsg so is it even possible?


title: How to create an Shapefile on ArcGIS Server side?

text: I am developing a script which gets an JSON array with point coordinates. This array is encoded to a python dictionary. From this I filter all points (latitude and longitude) and an ID. This coordinates are converted into arcpy- points by using arcpy.Point(). Now in a next step I want to create a shapefile of these points by using Code: `arcpy.CopyFeatures_management(ptList, r"pointTemp.shp")` My question is how to store this shapefile into a folder which is located on ArcGIS Server side. Do I have to give the skript a special path or location? One folder that is on ArcGIS Server side in my case is "C:\temp", I have full read and wright rights, but the script doesn´t create the shapefile in this folder. If I run the script local on my pc, all works fine and the shape is created and stored locally on "C:\temp". Any suggestions or tips? Thanks!


title: How can one select only displayed features using the select by rectangle tool?

text: I have a shapefile with many features, only some of which are displayed. I'm using the "Select Features by Rectangle" tool (the one where you click and drag a box over the features you want to select). I want to only select visible features within the rectangle, however the tool also selects features which are not displayed. Does anyone know how to avoid this? I'm using ArcMap 10.0 Thanks


title: How to repair broken Mosaic Dataset data sources?

text: I have a question. My computer has been replaced with a new Xeon Workstation. Before that, I had a temporary workstation and had ArcGIS 10 on it. While I was waiting for the new workstation to be replaced, I was working on mosaic dataset and worked on several clipped rasters to mosaic them. When I got my new workstation back, I had to backup all of my datasets to external hard drive and then load them to the new computer. So when I bring the mosaic dataset up to my ArcGIS 10 on my new computer, all I see is black! Did I do something wrong when I made the backup or should I run some fixes? Originally, what I did was that I had two NAIPs loaded from the ArcGIS server added them to my ArcMap and had a layer called mask to clipped the raster within the mask. I use them to mosaic them to created one mosaic file saved in Mosaic Dataset. HELP !! I hope I do not have to re-do all the work again. Thanks for your feedback. Thanks !


title: How to add a custom north arrow to ArcMap?

text: I created a custom north arrow in Inkscape, now I cannot seem to add it to ArGIS 10.1. I would ideally like for it to be displayed in the north arrow selector window (See attached image). How can I add my custom north arrow to ArcMap 10.1? ![enter image description here](http://i.stack.imgur.com/bCNqm.png)


title: How to clip the contour map with the polygon using R software?

text: Hi I am wondering how can I create the contour lines from the discrete x, y and z dataset and then add the shapefile to cut the contour lines by polygon. I want to clip the contour map by the shape file. I have created the contour map and added the shape file but I am not able to cut the shape file with polygon. I used akima package to create the contour map. I have posted this question in stackoverflow R help but I didn't get any help so far and I desperately need help on this issue. I will post my code and data links here : The script file image.scale.R can be found on the following location "https://www.dropbox.com/s/2f5s7cc02fpozk7/image.scale.R " The code that I have used so far is as follows:               ## Required packages      library(maptools)      library(rgdal)      library(sp)      library(maptools)      library(sm)      require(akima)      require(spplot)      library(raster)      library(rgeos)          ## Set Working Directory      setwd("C:\\Users\\jdbaba\\Documents\\R working folder\\shape")          ## Read Data from a file      age2100 <- read.table("temp.csv",header=TRUE,sep=",")          x <- age2100$x      y <- age2100$y        z <- age2100$z          ####################################     ##Load the shape file      #####################################      shapefile <- readShapePoly("Export_Output_4.shp")          fld <- interp(x,y,z)          par(mar=c(5,5,1,1)) filled.contour(fld)          ###Import the image.scale      source source("image.scale.R")      # http://menugget.blogspot.de/2011/08/adding-scale-to-image-plot.html               x11(width=8, height=7)      layout(matrix(c(1,2), nrow=1, ncol=2), widths=c(6,1), height=6, respect=TRUE)      layout.show(2)          par(mar=c(4,4,1,2))      image(fld,axes=T)      contour(fld, add=TRUE)     #points(age2100$x,age2100$y, pch=".", cex=2,legend=F)      plot(shapefile,add=T,lwd=2)      box()     par(mar=c(4,0,1,4))      image.scale(fld, xlab="Eastings", ylab="Northings", xaxt="n", yaxt="n", horiz=FALSE)          axis(4)     mtext("Salinity", side=4, line=2.5)      The output of the above code is as follows: ![enter image description here][1] Now, I want to get rid of the colored gradients and the contours from the polygon shapefile and only leave the intersection part. Any help is highly appreciated. Thanks.


title: What specifically is an 'attribute' and how does it work with the related chart?

text: I'm a new QGIS user (and GIS technology in general) and I'm trying to learn its fundamental features. I have a big dilemma about vector shapefiles, and points in particular: basically, I don't understand what specifically is an 'attribute', and how it works with the related chart (I guess that the 'ID' is required for a search within that chart) Please help me! Sorry for my bad english.


title: How to break (both) lines at their point(s) of intersection in QGIS?

text: I'm working with a road network where some of the road segments cross each other but no point is inserted at their intersection. Is there a way to break these line to shorter lines at their points of intersection either completely in QGIS? Thanks


title: Error 000539 in raster calculator in ArcMap 10

text: I'm having trouble averaging 12 rasters in the raster calculator. I keep getting a 000539 error, which I assume means that the syntax in the equation is bad. The problem is I dont know whats wrong with the syntax or how to remedy it. Any help/suggestions would be greatly appreciated. fyi. the equation I enter is ("raster 1" + "raster 2" + ... "raster 12") / (12) thanks, jonathan schmitt


title: LANDSAT Image classification and open source solution(s)

text: I have some experience with ERDAS IMAGINE in doing some supervised/unsupervised LANDSAT imagery classifications and now find myself in the position of needing to do the same but without the access to ERDAS IMAGINE. I've been looking at the GRASS literature and found this here and see that it has similar functionality/tools such as i.class/cluster/gensig as well as PCA and Tassled Cap analysis, which also hold some interest to me. I plan on looking at SAGA as well when I get the chance to see what it can do in respect to classification. I also found this thread as well here but much of what is there dates back a few years. My question is: are there other Open Source alternatives similar to GRASS or SAGA that I should look into as well that you may have used? Just doing my homework before I commit to something...and I am really only interested in the image classification abilities, or lack thereof. I will be using LANDSAT data for an AOI of 700 sq/mi and would like to utilize the ASTER spectral library. The class number would probably be fairly low as it is not a very diversified landscape. Operating system is Windows 7, 32bit(!) and no budget to speak of. Time frame will be mid-summer, full leaf on.


title: Configure SoundEx in ArcEngine

text: Is there a way to configure SoundEx properties in ArcEngine? This query: "100 Fairvew Rd returns "100 Fairview Rd" But this query: "100 Farview Rd" returns no results Is there a way to make the SoundEx a bit less rigid? Here are some concrete examples I would like to get to work if possible (ie. these are my specs) 100 cen => central ave 100 w laur ave -> 100 w laurel ave I'm using the IAddressCandidates2 interface to get my candidates.               IArray resultsArray = addressCandidates.FindAddressCandidates(addressProperties);


title: ArcGIS Geoprocessing Service Error 000732

text: I just published a python script as a geoprocessing service without problems. The script works on desktop perfectly. Now i wanted to run the geoprocessing service and i get this error. ![enter image description here](http://i.stack.imgur.com/DTa22.png) It says that the input-layer is not available or not supported. So i just searched for that layer and i found that it exists and i can load it into ArcMap, too. All other files the script uses have also been copied to the server. Does someone know a solution?


title: Where is the QGIS Parallel Line Plugin?

text: Somehow I believe to recall there does exist some ind of "draw paralel line" plugin for QGIS. Could anyone please point me to - I cant find it :-/


title: Where are the radiation models for Japan catastrophe?

text: I'm interested in seeing authoritative model results showing how radioactive contaminants from Japan might circulate through wind and water. ![map lie](http://i.stack.imgur.com/167Nw.jpg) The AGU blog points out this bad model, but in skimming the sites they reference I didn't see any other models. I see no mention of radiation at NOAA's Global Forecast System. Where are the reputable models?


title: issue with QGIS Scale Level

text: I have different layers of shape files. I have imported all shape files into PostGIS. Now i want to show 3 layers Country, State & District with different scale/zoom level I tried to set scale based rendering for maximum & minimum. But its not working. I am confused how to choose the Scale for it. I am using CRS as WGS84 - 4326 Thank You


title: How can I crawl directory tree and change Mxd relativePaths to true?

text: I have a large number of project folders where each project folder can have many subfolders. These subfolders can contain mxds, shapefiles, personal gdb's and fgdb's. In many cases the user forgot to check the "Use Relative Paths" before saving the mxds, so when the folder gets copied, datasource links become broken. I would like to write a standalone arcpy script that crawls through the directories, opens each mxd, sets the relativePaths to true, then saves the mxd. Is there a script that does this somewhere? The relativePaths property on the MapDocument object is readonly. Is there some other way to change this?


title: How can I call a command programmatically and change its defaults?

text: I am quite new to ArcObjects and development, but from what I have researched I can run a command programmatically like so: [VB]               Dim pUID As New UID     Dim pCmdItem As ICommandItem     pUID.Value = "esriArcMapUI.MxFileMenuItem"     pUID.Subtype = 4 'For Save As command     Set pCmdItem = Application.Document.CommandBars.Find(pUID)     pCmdItm.Execute      However I want to extend this example so when it executes it changes defaults to always the same folder location and filename. Is this possible with ICommandItem, or do I need another interface?


title: OpenLayers tooltips

text: Please point me to examples of combining OpenLayers and any external tooltip libraries like this. I try to use native OL tools for building tooltips here, but they not flexible and customizable. Thanks.


title: Attribute table field value question?

text: I have a field in my attribute table that is land ownership. That field is populated with all different types of ownership but there are redundencies that are just labled differently. Ex. Parks & Recreation AND Parks and Recreation. How do I combine these values within the table without creating a new layer or using dissolve?


title: Problem running i.gensigset in Grass. Any ideas?

text: It just returns `ERROR: Unable to create signature file`. It doesn't say why... Anyone had this before and found a solution?


title: Is it possible to use a WHERE clause with Grass v.in.ogr for DGN sources?

text: I am trying to import a .DGN with Grass v.in.ogr in QGIS. I am able to import the DGN perfectly, but I would like to filter the elements using the WHERE option. From searching, reading and testing I think it should be running correct, but it is erroring out. Any help would be appreciated. The attribute column holding the information I would like to Filter on is called Level and it is an Integer. I want to filter it on a value of 1 (see attached image).![v.in.ogr filer](http://i.stack.imgur.com/tzzXA.jpg)


title: Can TileMill be used to create Paper Map from OpenStreetMap source?

text: I would like to know which software is the best for making paper map from osm source. Is tileMill good for that ?


title: changing central meridian in QGIS causes problems

text: I want to change the central meridian of my world map (WGS84) in QGIS from 0 to -80. The problem is that this distorts the map every time I try it. ![QGIS project showing country outline data. Countries at +/-180 degrees from the new central meridian are unclosed polygons.](http://i.stack.imgur.com/5SBqX.jpg) This is what I do:   * go to Settings/Custom CRS and then click on Add new CRS   * copy existing CRS (choosing WGS84 / World Mercator)   * Edit its parameters (changing only the central meridian) to: +proj=merc +lon_0=-80 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs   * go to Project Properties/CRS and choose my custom CRS And then the whole map is distorted. It works fine when I do it in ArcGIS.


title: Get the total distance from a point to the end of a linestring

text: I have a list of GPS coordinates that make up a bus route:               busRoute = [(x1, y1), (x2, y2), ... (xn, yn)]      I also have a point that may not be exactly on that route (due to GPS error):               currentLocation = (x, y)      Using Python, how can I get the total distance from the bus's `currentLocation` to the end of the `busRoute`?


title: Why are Garmin provided maps inaccurate?

text: Collected tracks on my Garmin GPS (Oregon 450) or in BaseCamp seem to be not in the correct position. There's only a little detail with major highways but it looks as though the route is shifted south, or the highway is shifted north. When viewing the route in Google Earth or Strava, the route follows the road almost perfectly. Are the Garmin provided maps just really inaccurate or is there something else going on? Here's a couple pictures so you can see what's going on. In the Strava image, I'm travelling east-west when I'm north of the highway (which is correct), but on BaseCamp, it shows me travelling east-west when I'm south of the highway. This has happened with all the tracks I've recorded so far on this route. ![Strava image](http://i.stack.imgur.com/odRFq.jpg) ![BaseCamp image](http://i.stack.imgur.com/pmDas.jpg)


title: How to restrict user access to map services on ArcGIS for Server?

text: If we have all our map services published on one server, can we configure it such that our users can only see or access a only a subset of the map services? Thanks.


title: How to select the nearest lines to a point layer?

text: I am a beginner with GIS software and I'm working on QGIS. I have a problem joining the spatial layers. I have a layer with objects in one city and another layer with the streets. I need to attach each object to its street. I try to join attributes by location but I have not been able because it does not take into account the objects that are on the edge of the street (the roundabouts and bridges). So how can I make a join between a point and a layer of layer lines using the closest distance by quantum GIS.


title: Clipping buffers to census tracts

text: I have donut buffers generated to protect anonymity of household locations. They need to be clipped so that they are only inside the census tract that the original point was in. The code takes about 8 hours on a very fast computer. Any optimizations to the following code would be greatly appreciated:               ''' Clipbuff.py     '''     # Import arcpy module     import arcpy, sys, os, datetime          print datetime.datetime.now()     # Set up workspace and scratch workspace     workspace = "C:/Users/percy/projects/CEE/PerturbationProcessDataset/jitters.gdb/UTM/"     arcpy.env.workspace = workspace     scratch_workspace = "C:/users/percy/projects/CEE/PerturbationProcessDataset/scratch/"               # Allow overwriting     arcpy.env.overwriteOutput = True          # Local variables:     points = "OHAS_household"     id_Field = "OBJECTID"     polygons = "Geography_Tract"     buffers = "buff10"               # Open a search cursor on the buffers layer     rows = arcpy.SearchCursor(buffers)     print rows     # Create a for-loop to iterate through buffer IDs (via searchcursor),     # make feature layers of both the buffer and census tract layers,     # select containing census tract,     # clip by selected tract       x = 1     for row in rows:         myVal = row.getValue("SAMPN")         mySelection = '"SAMPN" IN ('+ str(row.getValue("SAMPN")) +")"         # Abstract the current row value         print datetime.datetime.now()         # Make a feature layer with the buffers layer, selecting the OBJECTID of the     current row/record         arcpy.MakeFeatureLayer_management(buffers, "buffers_layer", mySelection)              # Make a feature layer with the corresponding points layer         arcpy.MakeFeatureLayer_management(points, "points_layer", mySelection)              # Make a feature layer for the census tracts layer            arcpy.MakeFeatureLayer_management(polygons, "polygons_layer")              # Select by location on the census tracts layer         arcpy.SelectLayerByLocation_management("polygons_layer", "CONTAINS", "points_layer")              # Clip the queried buffer by its spatially coincident census tract         arcpy.Clip_analysis("buffers_layer", "polygons_layer", scratch_workspace + "buffclip" + str(x))         # Add 1 to the file suffix         x += 1         # Delete the feature layers             arcpy.Delete_management("polygons_layer")         arcpy.Delete_management("buffers_layer")         arcpy.Delete_management("points_layer")          # Set the scratch workspace as the current workspace to list the clipped buffer featureclasses     arcpy.env.workspace = scratch_workspace          # Merge all clipped buffers into one layer     mylist = arcpy.ListFeatureClasses("buffclip*")     mytarget = mylist[0]          del mylist[0]          arcpy.Append_management(mylist, mytarget, "NO_TEST")     print datetime.datetime.now()      New version of code:               # Import arcpy module     import arcpy, sys, os, datetime          print datetime.datetime.now()     # Set up workspace and scratch workspace     workspace = "C:/projects/CEE/PerturbationProcessDataset/jitters.gdb/UTM/"     arcpy.env.workspace = workspace     scratch_workspace = "C:/projects/CEE/PerturbationProcessDataset/scratch1.gdb/"               # Allow overwriting     arcpy.env.overwriteOutput = True          # Local variables:     points = "OHAS_household_subset"     id_Field = "OBJECTID"     polygons = "Geography_Tract"     buffers = "buff1"          arcpy.MakeFeatureLayer_management("Geography_Tract", "polygons" )     arcpy.MakeFeatureLayer_management("buff1", "buffers")     arcpy.MakeFeatureLayer_management(points, "points")          # Open a search cursor on the buffers layer     rows = arcpy.SearchCursor(points)     x = 1     #for row in rows:     row = rows.next()     myVal = row.getValue("SAMPN")     mySelection = '"SAMPN" IN ('+ str(row.getValue("SAMPN")) +")"     print mySelection     #print row[0]     print datetime.datetime.now()     arcpy.SelectLayerByAttribute_management("points", "NEW_SELECTION",mySelection)     ptcount = arcpy.GetCount_management("points")     arcpy.SelectLayerByLocation_management("polygons", "CONTAINS", "points")     arcpy.SelectLayerByAttribute_management("buffers", "NEW_SELECTION",mySelection)     polycount = arcpy.GetCount_management("polygons")     buffcount = arcpy.GetCount_management("buffers")     print ptcount, polycount, buffcount     arcpy.Clip_analysis(polygons, buffers, scratch_workspace + "clippy" + str(x))


title: In QGIS map composer, how do I get the Item Properties dialog back?

text: Yes, well if there is button there, one is bound to click it.... I closed both the Item Properties dialogs in Map Composer. And I can't get them to give me any joy no more... QGIS Trunk


title: Address locator style for my database?

text: I'm a student Studying GIS, and I need to create geocoding for my work. How can I create the adrees locator style for fields with my requirements and my specified shape file..... My attributes are               Id     Survey_No     Owner_Name     Cadestre_No     Soil_Type     Cult_Type     Water_Sour     Rights     Bebts


title: openlayers 3 geojson style

text: here's a part of my html/js file. I cant display my geojson and really dont know the reason why. debug only says : "TypeError: f is null" any1 help ? Thx               var map = new ol.Map({     target :'map',     layers: [         new ol.layer.Tile({             source : new ol.source.OSM()             }),         new ol.layer.Vector({             source: new ol.source.GeoJSON({                 projection: 'EPSG:3857',                 url: 'emprise.geojson'}),             style: new ol.style.Style({                 image: new ol.style.Circle({                     radius:3,                     fill: new ol.style.Fill({color: 'red'})                 })             })         })             ],     view: new ol.View({ <!--2D marche plus...2014 06 29 -->         center:ol.proj.transform([1.22, 47.435],'EPSG:4326', 'EPSG:3857'),         zoom:9                         })         });


title: Line style width property not being applied in QGIS 2.0.1

text: I change the line width of my style and click apply, but the property does not seem to be applying to my line data.


title: GeoServer Error while calculating BBOX

text: I made a custom PL/pgSQL function in PostgreSQL (this function is using pg_tsp and pg_makeDistanceMatrix functions) and then made an SQL View in GeoServer to extract the information from the db using this custom function. But when I am trying to calculate the BBOX Min Max limits in GeoServer, it is throwing the following error message. Please help me to fix this out. Thx. Zia > java.lang.RuntimeException: java.io.IOException: Error occured calculating > bounds for testnew3 at > org.geotools.jdbc.JDBCFeatureSource.getBoundsInternal(JDBCFeatureSource.java:528) > > ----- few more lines ------- > > at > org.mortbay.thread.BoundedThreadPool$PoolThread.run(BoundedThreadPool.java:450) > Caused by: java.io.IOException: Error occured calculating bounds for > testnew3 at > org.geotools.jdbc.JDBCDataStore.getBounds(JDBCDataStore.java:1258) at > org.geotools.jdbc.JDBCFeatureSource.getBoundsInternal(JDBCFeatureSource.java:521) > ... 102 more Caused by: org.postgresql.util.PSQLException: ERROR: Error: > matrix[num][num] in its definition. Where: PL/pgSQL function > "pgr_tsp_latlon" line 22 at RETURN QUERY SQL statement "create temp table za > as select * from pgr_tsp_latlon('select id, st_x(the_geom) as x, > st_y(the_geom) as y from az', ''||tb1||'', x1, y1, x2, y2)" PL/pgSQL > function "pgr_tsp_eucleadian" line 14 at SQL statement at > org.postgresql.core.v3.QueryExecutorImpl.receiveErrorResponse(QueryExecutorImpl.java:2062)


title: Postgis 2.0 simplifying multipolygon

text: I have a data set of national boundaries for the African continent. These are in a multipolygon table in postgis. Each country has multiple entries (rows) such that some of the countries in my table have 20 or 30 rows with individual geometries. I want to combine them multiple rows for each country into a single polygon (if possible). The main reason I want to do this is because I'm creating a map in QGIS with this as my base layer and want to label the countries. At the moment however I get 20 to 30 labels per country - I just want one! I have tried ST_Union which I thought would combine all of my individual geometries into one (per Country) But it only produced geometries for the countries where there was only one geom to begin with:( This is my query:                   SELECT        "africa_ALL".name,      ST_Union(geom)      FROM        public."africa_ALL"       GROUP BY name       ORDER BY name;      Can anyone help out at all, my sql leaves a lot to be desired and I'm desperately trying to get this map finished! All help greatly appreciated. Marty


title: How to preserve dgn cad file attributes when exporting from microstation and importing into Arcmap?

text: I'm trying to import water utility lines from Microstation into ArcMap. The water utility lines are symbolized by size in Microstation, but when I try to view the attribute table in ArcMap the attributes seem to have not transfered over, the pipe size attribute is blank. Any advice or a direction to turn? I'm poor with microstation/CAD but have a solid background in ArcGIS. Thanks!


title: Tool to save own color table with CMYK, GRB and HSL

text: Im looking for a tool to save my colors with CMYK, RGB and HSL. So that I can print and compare, write some notes to every color and easy extract. Any idea? Im mainly using Inkscape, Gimp and QGIS on xUbuntu. Its for maping purpose.


title: No visual after loading geotiff in Qgis Dufour

text: I installed the latest version of QGis 2.0 for mapmaking. Before I used the 1.8 version and I had no problems in loading Geotiffs. With the latest version, it loads the same files but it doesn't show any image. I haven't got any clue why in this version it doesn't load/show these files. Does anybody know?


title: How can I filter a MultiPoint Z / xyz points against a 2d polygon

text: I would like to filter point data from LIDAR scans (3d) by parcel plot lines (2d). I began with loading the lidar data into microsoft sql and performing a STIntersects between the MultiPoint and the Polygon. This provides a result, but loses it's 3rd dimension. I then have tried loading the lidar .shp file into the fiona python library, calling .filter() to narrow down the features, and traversing points checking .within() against the plot line polygon. This seems to work but is super slow. Next I loaded the lidar .shp into the rtree python library. However there doesn't seem to be the best 3d support. The methods that were obvious to me were limited by the inability to get a .bound on the shape that included the 3rd dimension. So when I added the shape to the index is complained the bound was 2d and the shape 3d. Finally I attempted to load each point manually into the rtree index so that I could correct the 3rd dimension manually on the insert. I have not validated that rtree will even work once I have a complete index. Has anyone figured out a non brute-force method?


title: Revert back from RGBA raster to original 1 band raster cell value

text: I have a number or single band raster grids that have been converted to RGBA (png). I no longer have access to the orignal single band raster data but have the png files. Is there a way to convert the png files back to the single band grid? Preferably using QGIS or other opensource GIS?


title: Multiline VBScript Label Expressions in Arcpy

text: I am able to set simple label expressions, but when I try to set multiline expressions (which all VBScript expressions are) it does not get copied into the expression box properly. I have tried using Python's triple quote technique for multilines but that only seperates each line with a tab in the expression window. Anyone know of a way to do this? If so, does the label class need to be set to 'advanced' first, or is it possible to automate this aspect as well? Alternatively, loading in precreated expression files would work just as well if not better if this is possible to script. This is with the use of the Maplex extension. **EDIT** Basically this code below needs to be pasted into the expression box with Arcpy. One line works just fine, but because there are multiple lines, I can't get it to paste it properly.               Function FindLabel ( [SUPP_NO], [SUPP_SUFF])     IF [SUPP_SUFF] = "000" THEN      FindLabel = "#" & [SUPP_NO] & VBNewLine & [Fuse_Count]"


title: Trying to instantiate geoprocessor with .Net vs 2010

text: Retrieving the COM class factory for component with CLSID {5374EC4C-1AA2-4829-A811-DE624ECEC23F} failed due to the following error: 80040111 ClassFactory cannot supply requested class (Exception from HRESULT: 0x80040111 (CLASS_E_CLASSNOTAVAILABLE)).               GeoProcessor gp = new GeoProcessor();      gp.SetEnvironmentValue("workspace", PICKSIMPORT);      IGpEnumList gdbs = gp.ListWorkspaces("", "");      string gdb = gdbs.Next();      while (gdb != "") { Console.WriteLine(gdb.ToString());      gp.SetEnvironmentValue("workspace", gdb);      IGpEnumList fcs = gp.ListFeatureClasses("", "", "");      string fc = fcs.Next();      while (fc != "") { Console.WriteLine(fc.ToString());      fc = fcs.Next();     gdb = gdbs.Next();


title: Create and immediately update field using Field Calculator

text: You can create fields and once created, populate them with data using the Field Calculator. But can you do both at the same time using an expression? Or must you do it in 2 steps?


title: Raster Calculator dynamic map algebra expression

text: Is it possible to dynamically construct a Map Algebra expression that is dependent on the output of another tool, either using model builder or arcpy? For example, say I have a model that performs a raster reclassification that has a list of raster inputs and their respective outputs, i.e. rasterOne -> reclassification -> outputRasterOne rasterTwo -> reclassification -> outputRasterTwo **EDIT:** To clarify, the number of input rasters is not known and therefore the number and names of the output rasters are not either. Because of this, I cannot hardcode them into the map algebra expression. I would then want the map algebra expression to be similar to:               %"outputRasterOne"% + %"outputRasterTwo"% + ...      ![Example Model](http://i.stack.imgur.com/Fn4uq.png)


title: How to make the thickness of traffic volume flow symbology proportional to traffic counts?

text: I can scale the bandwidth proportionally to the traffic volumes; however, I cannot show the directionality of the movements. Volumes for directional moves are stored in two different columns of a database. My workaround for this is to import the layer twice and use two different symbologies, one for each direction. **Is it to possible to do this without having to import the same layer twice?** Note that I need to show the bandwidth starting from the link so I cannot use the same offset because it create blank space between the opposite directions (it seems to me that the offset start from the middle of the thckness). **Edit:** in the uploaded images you can see that the offset depends on the width of the links, so I is not possible to find a single offset value matching for all the links ![Not Desired Effect](http://i.stack.imgur.com/wTXhV.png) ![Desired Effect](http://i.stack.imgur.com/JJhvq.png)


title: How to calculate accuracy of latitude/longitude

text: I'm updating a Windows Mobile application to capture location information using the built-in GPS hardware. The Microsoft API provides a method _GPSGetPosition_ that returns the following struct (ref:GPS_POSITION):               typedef struct _GPS_POSITION {         DWORD dwVersion;         DWORD dwSize;         DWORD dwValidFields;         DWORD dwFlags;         SYSTEMTIME stUTCTime;         double dblLatitude;         double dblLongitude;         float  flSpeed;         float  flHeading;         double dblMagneticVariation;         float  flAltitudeWRTSeaLevel;         float  flAltitudeWRTEllipsoid;         GPS_FIX_QUALITY     FixQuality;         GPS_FIX_TYPE        FixType;         GPS_FIX_SELECTION   SelectionType;         float flPositionDilutionOfPrecision;         float flHorizontalDilutionOfPrecision;         float flVerticalDilutionOfPrecision;         DWORD dwSatelliteCount;         DWORD rgdwSatellitesUsedPRNs[GPS_MAX_SATELLITES];         DWORD dwSatellitesInView;         DWORD rgdwSatellitesInViewPRNs[GPS_MAX_SATELLITES];         DWORD rgdwSatellitesInViewElevation[GPS_MAX_SATELLITES];         DWORD rgdwSatellitesInViewAzimuth[GPS_MAX_SATELLITES];         DWORD rgdwSatellitesInViewSignalToNoiseRatio[GPS_MAX_SATELLITES];     } GPS_POSITION, *PGPS_POSITION;      So I can determine the latitude and longitude but I would like to be able to calculate the accuracy of the result. The application runs on a Trimble Juno and the _GPS Controller_ application on the device provides an accuracy value in metres. Is this something I can determine from the values in the struct or is this something specific to the hardware that would require a native function call? Is the accuracy value based on heuristics or is it a common calculation across different hardware?


title: How I can merge three rasters with sum cell

text: I am working by ArcGIS 10 and I have three raster all of raster have value cell and I want output raster have the sum of pixels values


title: Hover effects using CartoCSS: Is this possible?

text: I'm trying to make a simple map of parcels in TileMill. Right now when I hover over a parcel, the information for that parcel appears in a teaser. That's simple enough. However, when working with particularly small parcels it may not be immediately clear to the user over which parcel one is hovering. What I would like is for the outline of a particular parcel to glow when it is being hovered over. I know in CSS there is a :hover selector, but it would seem that either: a) There is not a direct translation in terms for between CSS and CartoCSS for this selector. -or- b)A hover-like selector is not available in CartoCSS. I've included an example below where I attempted to use the hover selector from CSS. In this case it would seem the hover selector is ignored and the parcel outlines are drawn again using the given attributes. Nonetheless, I offer this faulty example as a potential starting point for any suggestions on how to implement hover features using CartoCSS.               #parcels {       line-color:#151515;       line-width:1;       polygon-opacity: 0.6;     }     #parcels::hover {       line-color:#0AFF47;       line-width:1.25;     }      Thanks for any suggestions!


title: How to render iD editor-style streetmaps?

text: Is there a simple way to render OSM data (let's start with just streets and names for a single metro area) client side as nodes and edges like the new iD editor? I'm comfortable with d3 but am having trouble seeing exactly how they're pulling it off.


title: How to zoom in/out selected elements of MapControl?

text: I have a mapcontrol that includes many elements: text, line. I used this code to add elements into mapcontrol:                pGraphicsContainer.AddElement(pElement, 0);      I can use zoom button to zoom in/out the mapcontrol, but how can I zoom in/out only seleted elements ? Is there any buttons for this purpose in Arcgis ? Thanks and regards, Tai


title: Protocol strings for ISO-19139 metadata

text: I'm working with some ISO-19139 (GMD) metadata. In the `gmd:CI_OnlineResource` element, there is a child element called `gmd:protocol` which should describe the protocol of the resource (e.g. WFS, WMS, etc.). The GMD XML schema seems to let you just put any character string in this, but I've seen different metadata tools behave differently depending on what's in this element. Is there a list of well-established protocol strings to use in this element? I've seen, for example, `OGC:WMS` used to describe a WMS service, but sometimes I'll find variations on that string which contain version numbers. Is there an official (or unofficial) list of strings I can reference?


title: How to add Points to a jvectormap in Miller Projection

text: On jvectormap you can pull down a map of map of austria. I want to be able to plot points (using WGS84 coordinates) onto that map. The map has this information inside: > "projection": {"type": "mill", "centralMeridian": 0.0} > > "bbox": [ { "y": -5942433.950285355, "x": 1060387.4860068962 } , { "y": > -5569704.340427333, "x": 1909872.2578980657 } ] That works. But I need to figure out how the x and y coordinates there were calculated, should be meters from meridian. Figured they are map (svg) file dimenions dependend and stumbled to this and this which is the same thing. Looks like this calculates pixel points more or less accurate, but this is not anything near upper "bbox" thing. Any suggestion welcome.


title: What is the proper way to delete an SQL Server database in order to delete an SDE geodatabase?

text: I'm new to this SDE stuff and would like to know the best way to cleanly delete an SQL Server database in order to delete the SDE geodatabase. Is it as easy as right clicking on the database in SQL Server and hitting delete? Or is there a more preferred way to do this? Also, is there a way to rename the SQL server database that the SDE geodatabase is using? Or is it best to just start fresh with a new database? Thanks in advance.


title: Remove a custom contents view from ArcMap TOC in VB .NET?

text: I am using the ESRI sample, TOCLayerFilterVB2010, which implements a custom TOC view for map layers, as a prototype. For my purposes, I need to be able to turn it on/off, meaning add and remove that view from the TOC horizontal icon bar programmatically. Can someone explain how the custom contents view in this sample gets instantiated? I find that if the dll is registered, it appears. The only way to get rid of it in any ArcMap document is to unregister it. Deleting the Normal.mxt template does not have any effect. Ideally, I would like to have a command toggle that instantiates and destroys the view. I cannot figure out how to control its creation.


title: Is the ServerAction Method Required in a ImapServerCommandAction?

text: I am using an ESRI toolbar control buddied with a map control. I am using a command. Below is a simple example. My question is to clarify - having the ServerAction method is a requirement when you are using the toolbar control and server side action?               public class CustomExtent : IMapServerCommandAction //implements the IServerAction interface     {         public CustomExtent()         {         }                                   #region IServerAction Members          public void ServerAction(ToolbarItemInfo info)     {        Map map = info.BuddyControls[0] as Map;        map.Extent = new Envelope(-90, 10, -70, 5);     }      `


title: Saved maps in geoexplorer (geoserver layer security)

text: I am a bit confused with the save maps function in geoexplorer. I have a vps server with the opengeo suite 3.1, I added a table to the vps local postgis db server to the geoserver.public database. I registered this table as layer in geoserver and it's working fine in geoexplorer composer. My problem is when I save the map in geoexplorer and log out and log in back (as admin) I cannot use the query button and the edit button (they are inactive). It is working fine (so I can use the query button on saved maps) if I don't change the layer security in geoserver. But I changed because I don't want to allow to anonymous users to see any layers. with this settings: the anonymous user cannot see the layers preview in geoserver and in geoexplorer. Admin can see every layer after logged in. But If I save the map (not publish), the query button and the edit button will be inactive. _._.r ADMIN,GROUP_ADMIN _._.w ADMIN,GROUP_ADMIN with the default settings: _._.r * _._.w * The anonymous user can see the layers, admin can use the query button on saved maps. I would like to hide the layers from the anonymous users, and I would like to use the query button on saved maps. Is it possible to do both of them at the same time? I'm sure is possible, just I don't know how can I do that. I cannot solve this problem and it's very frustrating. It would be great if somebody can give me some idea to solve this issue. Thanks, encleadus


title: Free website host for embedding tilemill/mapbox web map

text: I am hoping to design a simple interactive web map using TileMill. I am then planning on uploading that map to Mapbox.com, where I understand I will be given an embed code to put on a website. My question is: will this mapbox map seamlessly embed onto a Wordpress or Wix website? Or should I use a different free website platform? Nota bene: I have little HTML editing experience. Thanks!


title: ArcObjects: Perform action upon refresh command in the TOC

text: We are using custom layers in our plugin used in ArcScene/ArcMap. For quite some time we were trying to figure out how to hook into the event occuring when the user selects "Refresh" in the right click menu item of a layer in the TOC. Any advice how to do that? Thanks.


title: How to define ITRF 92 and ITRF 2008 using Proj4 parameters in QGIS?

text: I know how to define a CRS, but this particular one is giving me a headache. I was unable to find the EPSG code at spatial reference. The main reason for asking this is that raster and vector data available from INEGI (Mexico) comes in ITRF 92. Also new cartography from the same source comes in ITRF 2008. When I load a raster or a vector layer from INEGI their CRS properties only said _"unknown datum based on GRS80 ellipsoid"_. INEGI says that ITRF and WGS84 are "equivalent" to make maps in scales smaller than 1:50 000. Nevertheless I usually make maps in scales greater than 1:50 000 Any ideas are welcome, Thanks beforehand.


title: Geoprocessor Select By Location Using Features Instead of Layers

text: I have a simple spatial feature selection function that includes:                public static void SelectByLocation(IFeatureClass sourceFC, ILayer selectLayer)             {                 Geoprocessor GP = new Geoprocessor();                 SelectLayerByLocation selectLoc = new SelectLayerByLocation();                 selectLoc.in_layer = selectLayer;                 selectLoc.select_features = sourceFC;                 RunTool(GP, selectLoc, null);                 ArcMap.Document.ActiveView.PartialRefresh(esriViewDrawPhase.esriViewGeoSelection, null, null);             }      It works correctly passing an IFeatureClass and an ILayer, however I would like to pass a selected subset of records from the IFeatureClass (for example: only the record with the most recent date or count). I've tried changing it to accept an ILayer as the first parameter and applying a selection to that with no result. I've also tried passing an IFeature object of the intended feature dynamically cast as an IFeatureClass. This latter method returns every row in the select features layer. For reference, I'm calling it now like this:                List<string> affectedFeatures = featureSelector(polyFC, layerList["SelectionLayer"], SelectList);      Where featureSelector() calls the spatial select with "polyFC" as IFeatureClass, layerList["selectionLayer"] as a reference to the ILayer selection layer, and SelectList as a list referencing fields to return from the selection layer. So, I suppose what I'm asking is: is it possible to pass an IFeature to the first method and get a proper spatial selection returned? For what it's worth this is always polygons selecting polygons.


title: Error installing/uninstalling GeoMedia WebMap Professional_6 on Windows 7 (64 bit)

text: On installing GeoMedia `WebMap Professional_06.01.03.08` on `Windows 7 Professional 64 bit`. Following Error window pop up at the end of installation process or If I uninstall it. ![enter image description here](http://i.stack.imgur.com/7dl1Y.png) I have .Net Framework 2, 3.5 and 4 installed. Thanks!


title: Webmapping with pgRouting and Mapserver for Beginners

text: I want to use pgRouting for showing shortest path on hiking trails. I created web map so I need to enable users when they click on start and end point somewhere on hiking trails, that it shows shortest path between those two points. I am beginner with it, so I just need to tell me where to start. I have those lines in postgis database and I am working with phpPgAdmin and Mapserver.


title: Python: Combine multiple Shapefiles

text: i have got a lot of Point-Shapefiles (r20010101.shp, r20010202.shp, r20020101.shp...). Each of them has the same amount of points at exactly the same location. In a column that is named after the Date e.g. "D20010101", where i have got unique Values for every Shape file. What i want to do now, is to Combine or Merge or Join all these Shapefiles in a python command.   The Result should be ONE Shapefile with all those point (but please NO duplicate ones!!).   The Attribute Table of this Shapefile should look like this: FID, Shape, Id, [Value of Date 2001/01/01], [Value of Date 2001/02/02], ...... I Hope that somebody does understand me :)   This has to work, i'm frustrated because I can't find a easy way to do this :( Thank you. Greetings, Markus


title: How to store raster data in spatialite?

text: I want to cache my raster data(large geotiff files) on local database. I thought to use spatialite. I found **librasterlite** which provide functionality to do that. Can anybody point to some tutorial on it. Or is there any other alternative. Programming environment - C++ * * * I want to store geotiff data to spatialite database. I am using gdal to read metadata and reading tiles converting them to hexwkb format and storing it in raster table (id (integer) and raster (blob) as columns).Table is created but when i am opening my database and select raster data column it shows nothing. I mean i can not explore data and see the tile image. Can any body tell me what should be the format in which we should store the tiled data. Are there any function which can automate these things. Like I tried to integrate rasterlite_load in my code but there were memory leakage in the code. Is there any function in spatialite or rasterlite which can take char buffer of tiled image and store it into blob data type in table (and corresponding function for reading).


title: Renaming 1,200 rasters in File Geodatabase using ArcPy?

text: I've been tasked with renaming a large amount of raster files within a file geodatabase (ArcGIS style). I'm taking weather data which has a name format as such: "G"mmddhhhh I'm trying to add a year tag in between the G and month: "G"yyyymmddhhhh Here is the function that I've written:               # ----------------------------------------------------------------------------- #     @timeit     def renameRasters(FGDB):         origdir = arcpy.env.workspace         arcpy.env.workspace = FGDB              for rstr in arcpy.ListRasters("G*"):             try:                 newrstr = "G" + "2014" + rstr.replace("g", "")                 arcpy.Rename_management(rstr, newrstr)                  except:                 print "Failed to rename: " + rstr              arcpy.env.workspace = origdir         return None     # ----------------------------------------------------------------------------- #      This function works fine for a small to medium amount of rasters within a FGDB and takes about 1.5 seconds per raster renamed. I run into issues when there are a large number of files to be renamed. 1,200 files takes > 10 hours. At the same rate as the smaller raster sets this should take < 1 hour... Does anyone have an idea on how to make this process go faster? If not I'll set it up to run over the weekend. I should mention that I am using ArcGIS 10.2.2.


title: Adds ArcGIS controls in the Visual Studio Toolbox window for dragging and dropping into a design view surface

text: I want to drag and drop MapControl into a design view, but I have an error: System.Runtime.InteropeServices.COMException(0x80040154):Classe non enregistrée (Exception from HRESULT:80040154(REGDB_E_CLASSNOTREG)). I find that your map control is not registered or Class not registered, but I don't know how to do it in VS 2008 SDK 9.3, please help.


title: SOC Performance matrix

text: I am monitoring performance of one of our SOC machines. Attached below are the stats that I am getting. it does say request are succeded but I am not sure that it is entirely true. It must take at least sometime. I am looking for some help on how to interpret these results. Are there any ESRI resources that describe what these numbers mean. Thanks jay               Service Usage Time:         Total number of requests:    1          Number of requests succeeded:    1          Number of requests timed out:    0          Avg usage time:   11.955999 Seconds         Min usage time:   11.955999 Seconds         Max usage time:   11.955999 Seconds         Sum usage time:   11.955999 Seconds          Service Wait Time:         Total number of requests:    1          Number of requests succeeded:    1          Number of requests failed:    0          Number of requests timed out:    0          Avg wait time:   0.000000 Seconds         Min wait time:   0.000000 Seconds         Max wait time:   0.000000 Seconds         Sum wait time:   0.000000 Seconds          Service Creation Time:         Total number of requests:    0          Number of requests succeeded:    0          Number of requests failed:    0          Avg creation time:   0.000000 Seconds         Min creation time:   0.000000 Seconds         Max creation time:   0.000000 Seconds         Sum creation time:   0.000000 Seconds


title: ArcMap Form/Button VB.NET MxDocument hook

text: \---UPDATE--- see below I have a scaling raster tool that I'm attempting to migrate from ArcGIS 9.3.1 VBA to VB.NET (Visual Studio 2010). I have a BaseTool that is used to put two graphic points defining the left and right extent of the image that I want to scale. After the second point is added to the view a form opens to input the scaling parameters. My issues is when I go to click the scale command button (within the form) I get the "Null Reference Exception was unhandled by user code, Object reference not set to an instance of an object", error. The error occurs on this line:               pMap = m_pMxDoc.FocusMap          I have my ESRI references imported into the project for both pMap and m_pMxDoc. It seems like the form is not hooked to the ArcMap document? I tried to hook the form on load via the code below, but I get the same error.               Private Sub frmScaling_Load(ByVal sender As Object, ByVal e As System.EventArgs) Handles Me.Load                 Dim hook As Object = Nothing                 m_application = hook                 m_pMxDoc = m_application.Document     End Sub      Any suggestions of what I'm missing would be greatly appreciated. Thanks \---UPDATE--- I created the docable window class and memebers. The On Create sub now is:               Public Sub OnCreate(ByVal hook As Object) Implements ESRI.ArcGIS.Framework.IDockableWindowDef.OnCreate                  m_application = hook             m_pMxDoc = m_application.Document          End Sub      Now when I click the command button within my form to scale the image I get a LoaderLock error (Attempting managed execution inside OS Loader lock. Do not attempt to run managed code inside a DllMain or image initialization function since doing so can cause the application to hang). The error is still on this line:               pMap = m_pMxDoc.FocusMap      Is this error present because my form is open while the scaling action is trying to take place? So, the hook is now on the member "On Create", do I also need to hook within the button click on the form? Are there any other good VB.NET form/code examples within ArcMap that contain some sort of button that interacts with the ArcMap view? I have only found this one,


title: Most suited Remote Sensing

text: What will be the most suited remote sensing data to use for mapping highways in North America?


title: ArcGIS Python Tool - Importing Custom Script Into ToolValidator Class

text: I had posted a question last week about customizing a ToolValidator class and got some very good answers. In working with the proposed solutions, I have created a custom module that performs queries on a db, and will be called by both the ToolValidator class (to provide values for the drop-down lists) and later in the geoprocessing script (to get other parameters based on items selected in drop-down lists). However, I can't seem to actually call the custom module in the ToolValidator class. I have been trying to append to path with no luck. When I try to apply these changes to the script, I get a runtime error: [Errno 9] Bad file descriptor. If I comment out import line, no errors.               sys.path.append('my_custom_module_directory')     import my_custom_module      Many of you may be asking why don't I just implement a custom tool with ArcObjects. The reason is that my end users don't have the privleges necessary to register ANY dlls on their computer. UPDATE: This was happening to me in ArcGIS 10. Interestingly enough, I was initially appending to the path inside of the initialiazeParameters function of the ToolValidator class. If I do the append outside (i.e., on top of) of the ToolValidator class, everything works as expected.               sys.path.append('C:/Working/SomeFolder')     import somescript -------->THIS WORKS          class ToolValidator:       """Class for validating a tool's parameter values and controlling       the behavior of the tool's dialog."""            def __init__(self):         """Setup arcpy and the list of tool parameters."""         import arcpy         sys.path.append('C:/Working/SomeFolder')         import somescript -------> THIS DOESNT WORK         self.params = arcpy.GetParameterInfo()      **UPDATE 2:** I think I found the true cause of my problem. In the code snippets in this post, I have been appending what appear to be real paths (ex C:/Working/SomeFolder) to the sys.path. In my actual ToolValidator class, I was constructing a relative path using `os.path.dirname(__file__)` +"\my_special_folder...". I was anticipating that `os.path.dirname(__file__)` would return the path of the toolbox, since that contains the ToolValidator class. I have come to find this is not the case. As far as I can tell, the ToolValidator class is never actually written to a .py file, and I speculate the the this code is passed to the python interpreter in memory, so `__file__` is useless, or some temp script is persisted and then execfile(path_to_script) is called, again rendering `__file__` useless. I'm sure there are other reasons I am missing as well. Long story short, if I use a hardcoded path, the sys.append works anywhere, relative paths don't work so well in the ToolValidator class.


title: Locations of Geodetic Survey Benchmarks

text: Does anyone know where I can find a list of geodetic benchmark locations for Ontario? (If a list even exists that is...) Thanks


title: Land Use/ Land Cover for Hawaii

text: Talking to a friend of mine in Hawaii and I am from New Mexico. My question to you all is to see if there is a land use/land cover dataset for Hawaii. Thanks.


title: The graduated symbology doesn 't take into account the attribute I need in QGIS

text: I m working on a project on QGIS I have a point shapefile vector that represents my store In the attribute table, i have bunch of data such as size, unit per transaction and most importantly store turnover When i want to graduate the symbology, the turnover isn't present in the list i can choose from Does someone know why? and how to fix it? Thanks


title: ShapeFile to GeoJson conversion

text: I am working with `OpenLayers`. I want to convert `ShapeFile` into `GeoJson` format and send it to browser to show with `OpenLayers`. But problem is when the size of my shape file is so big like 60MB or bigger, then it is not convenient to send whole `GeoJson` data to browser. Anyone can help me, how can i minimize the `GeoJson` file size while my `ShapeFile` size is still very large. Is there any tools or techniques by which i can encode `GeoJson` data to minimize size and send it to client side as well as decode this data to show in map. Thanks


title: File size inflation normal with gdalwarp?

text: After using `gdalwarp` to project and align-to-grid (via `-tap`) a number of rasters I noticed that the output rasters were significantly larger than the original rasters. A fairly thorough web search turned up this Trac issue: > Frank Warmerdam explained the reason: > > "On careful review, the difference in the file in question is because > gdal_translate uses the TIFFWriteScanline() interface to write the output > file from within GTiffDataset::CreateCopy?(), and this only writes as much > of the final 'strip' of the file as is required to complete image area. > **But gdalwarp goes through the blockio interface which writes the complete > final strip, even the portion that falls off the end of the file."** This Trac issue is ~7 years old, however, and I know some changes to the GDAL utilities, including `gdalwarp` have been made since. I'd like to know if the above reasoning still holds and if the file size inflation I'm seeing is "normal." The word "normal" here might be taken to mean _unsurprising_ or _expected_ but, more importantly: **is there anything that can be done to mitigate the effects i.e. reduce the output raster file size?** Below is a table of the file size inflation I'm experiencing.               Input File Size (bytes)     Output File Size (bytes)    Inflation     1437380431                  1698334217                   18%     1428001178                  1698334433                   19%       41683165                   137036637                  228%      The input TIFF files were created in ArcGIS and thus have external Worldfiles, XML and DBF files but these do not make up the difference in file size. Here is a sample `gdalwarp` call as I've used it in all of these cases; the actual execution was handled by a Python `subprocess` (`subprocess.Popen`):               $ gdalwarp -tap -tr 30 30 -t_srs "+proj=aea +lat_1=20 +lat_2=60 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs" -co "COMPRESS=LZW" input_file.tif output_file.tif      I understand that in rare cases compression makes a larger file, but the effect is the same without the LZW compression. The ratios in the table are with LZW compression.


title: GPS Coordinate Handling

text: I have a text file with consecutive coordinates of gps data for ongoing and incoming route(on Same road). i want to convert them in line format. how can i normalise those extra coordinate so that i get only single line for a single route.


title: Editing WFS with max feature limit in QGis

text: I'm serving a large layer (500 thousand polygons) using Geoserver WFS and what to be able to edit it using QGis. I've set a feature limit on the layer so I wont overload GeoServer with the request. I can successully load and edit the WFS layer in QGIS but as soon as I reach the WFS feature limit, QGis will not bring any new polygons as I pan through the map canvas, even when I hit the refresh button in QGis. I was under the impression that QGIS would keep on updating the WFS layer whenever I moved around and that the feature limit in GeoServer would be a limit to a single request and not a total request. Am I getting things wrong?


title: Geoprocessing service saved data to scratch gdb

text: I created and published a geoprocessing service from a toolbox. It works with a file geodatabase on my local drive. When run as a tool in ModelBuilder, it worked well and wrote data into the file geodatabase perfectly. However, when I published it as geoprocessing service and triggered it via REST service, it wrote data into the scratch.gdb in the job directory on the ArcGIS Server directories. How can I configure the tool/model, or service or arcgis server to make it work with the original file geodatabase instead of the scratch file geodatabase? Thanks a lot in advance! Cheers, Alex


title: How to deal with the situation when two rasters stacked together in arcgis?

text: I am doing a project about petroleum geology. There are several formations needed to add attributes (e.g. name, thickness and age). I have the raster data of thickness, so I created a raster catalog and loaded every raster in. However, these formations have different ages. The older formations below have a smaller distribution area, and the younger one above is much bigger, so the old formation is shown totally inside of the younger one. That makes the display quite confusing. So is there anything I can do to make the display clearer? Thanks. The following is the the picture: ![enter image description here](http://i.stack.imgur.com/TgS1q.jpg)


title: Problem with conducting a viewshed analysis

text: I am using the freely available dtm data from the OS UK website (Land-form panorama and its in the ASCII format). when you download the data from the website, the dtm is found in different subfolders (labelled according to the national grid squares, e.g. NN, NO, HY etc.). within each folder there are number of ASC files (which looks like tiles when you view on ArcGIS - I am not sure what the correct technical term is). My sites are spread out across an area of about 23 national grid squares - so there is a large number of ASC files that needs to be added. When I try to open the dtm data on ArcGIS 10.1 I encountered a problem, I can only open the ASC files from about 4 national grid squares before the program freezes and shuts down. The primary focus of my study is to conduct viewshed analysis for my sites. when I add the various ASC files and try to conduct a viewshed analysis it looks like this: ![enter image description here](http://i.stack.imgur.com/Lfn6q.png) The problem is my sites are in clusters and some sites are found spread out into the neighbouring tiles. it looks like when I try to do the viewshed analysis it doesn't take into account the sites that are found in the nearby tiles - which is a huge problem for me. when I go to the 'input' place on the viewshed screen I can only choose one asc file at a time. is there a way to choose more than one? ![enter image description here](http://i.stack.imgur.com/39Pgo.png) I was wondering whether anyone out there could tell me how to add all the asc files from all the different grid squares without the program shutting down and whether there is any way to stitch all the tiles together to make it look smooth and complete (once again I am not sure what the technical term might be for this) so that I can run a proper viewshed analysis. Any information or suggestions would be massively appreciated! thanks Fezter I really appreciate it! I tried doing what you suggested but it still comes up looking like separate tiles: ![enter image description here](http://i.stack.imgur.com/JFcTC.png) Am I doing something wrong? I am also having trouble performing the viewshed analysis here (although I can perform a hillshade), I was able to do it before I created the mosaic dataset. The other problem is I still have to do the hillshade and viewshed analysis one asc file at a time. is there a way around this? one other thing - when I added the raster data to create a mosaic dataset - it comes up as empty green squares - the terrain only becomes visible if I calculate the statistics for each individual asc file - it remains blank when I try to do a batch statistics. I am not sure what I am doing wrong, if you could help me it will be greatly appreciated!


title: How can I see the layers in an .OSM file?

text: I have an osm file. I am trying to convert it to geojson with ogr2ogr. I get an error that ogr2ogr can't handle multiple layers. How can I tell what layers are in the file and remove them? I see from How to remove layers of detail from an osm file? that OSM files might not even have layers. So what is the meaning of this error message?


title: Advice on learning Remote Sensing software?

text: Can you please advice me on which Remote Sensing software ? Which one is popular ? Which one is the least ? Which one is the easy to learn curve ? Pros and Cons welcome ?


title: How to build a portfolio for a GIS/Java career path?

text: I have been a GIS Analyst for over 10 years and am currently studying java (50% through the course). My plan is to develop GIS applications mainly focused towards open source solutions. I would like to develop some form of a portfolio to gain more experience, but I am not sure what the best way to do this would be. I have thought about joining some groups/communities like Geotools, but I am concerned with my lack of experience. I feel I would be more of a hindrance than a benefit. The other idea, is that I could look at developing some in house applications or even some Android application to build up some credibility. If anyone could provide some suggestions, or share their similar experiences, it would be greatly be appreciated. **How should I go about building a portfolio from scratch to further my open-source GIS App career aspirations?**


title: Have to use OpenLayers 2.10, can I force the browser to cache the map tiles and re-use them?

text: OpenLayers introduced client side caching in Version 2.11 and 2.12, along side their performance improvements for other map functions. Needless to say, these functions are intertwined and rely on each other so I don't think I can just pull out the caching code and use it with version 2.10. Is there any way to force the client browser to cache the map tiles and re-use them on subsequent requests? I am hoping this can be done for my WMS TileCache tiles retrieved from my ms4w Apache server. Cheers * * * Here is a live example on my host: http://www.spatialanalysis.ca/cachemytiles/ Check out the network panel for my WMS layer and then check the OSM base layer. The OSM base layer returns "from cache" tiles once they've been loaded once. I have a feeling this is an Apache setup option?


title: Recreation data for Ontario (Canada)?

text: I am aware that a similar question has been asked HERE but I am looking for comprehensive sources of recreational data (specifically for the province of Ontario, Canada). Land Information Ontario has some data ( **I have everything they can provide** ), but I am looking for something covering more categories (and fairly current if possible). I'm hoping to find golf courses, campgrounds, hiking/biking trails, kayaking/canoeing routes and other data of this nature. Other categories are welcome if you have any suggestions. Is there an all-in-one site where I can get such data?


title: How to create cross-sections in ArcGIS?

text: I have some problems to find the correct tools for a university project I am working on. Would be really nice if someone could help me to find the right tools. Is there the possibility if I have a polyline (named "Linedivide") which is divided into, let's say, 5 parts to make ArcGIS draw perpendicular lines of special length at the middle of all 5 parts of the Linedivide? Can someone tell me where I can find a tool for this in the toolboxes? Thank you!


title: Is there a way to export QGIS layers to KML?

text: has had anyone had any luck with the new experimental plugin qgis2kml? I keep getting the error msg **layer : format not yet supported** I had a look but the author has not provided any documentation, although I can't access the author's homepage at the moment. I tried different projections, but other parts of the error message refer to symbology.


title: Plotting latitude on to a (Mercator projection) SVG representing a country map

text: I want to convert (google sourced) **longitude and latitude** coordinates to **x and y** to place points on a **SVG map** I can get the longitude converted to X but I can't nail latitude. http://jsfiddle.net/chrisloughnane/an3BZ/17/   * red dots = place holders   * green dot = calculated position from longitude/latitude I have read so much about Mercator projections and other projection system API's I am now completely confused. I also followed Proj4JS library threads but couldn't find an example that was close to my task. I attempted to emulate @Raphael Wichmann and @betamax solution, unfortunately I got something wrong. Could someone please have a look at my jsfiddle and see if it's an obvious mistake? The second test **function secondconvert(latitude, longitude)** (bottom of javascript pane) provides me with a reasonably accurate **x coordinate** i.e. if I manually enter the **y coordinate** `b.ylat` change to `265` it covers my left red dot place holder nearly perfectly. I think I'm close, any help would be really appreciated. **Note:** the longitude/latitude bounds I used for my SVG map are best guesstimates, worked surprisingly well for the longitude calculation. tia. * * * Original SVG from here.


title: Opensource solutions for finishing maps

text: After seeing the $600 price tag for Adobe Illustrator, I am especially interested in learning open source solutions for finishing maps created in a GIS. What open source, or commonly available, software do you use to add the finishing touches (e.g. arrows, textures, colors, graphics, text, etc) to maps created in Arc, QGIS, or your choice of GIS software? Examples welcome.


title: Submodel not executing within main model

text: My main model contains the following:   * Workspace variable `Input GDB`   * Feature class `Input Feature Class`   * SubModel 1 that produces `OUTPUT` feature class ![Submodel](http://i.stack.imgur.com/sisTc.png) `Parse Path` is used to obtain the name from the input feature class and is stored in the variable `Name`. This is used as an inline variable for the various outputs of the processes inside `SubModel 1`. The workspace variable is also used as an inline variable. My problem now is that the submodel is not executing. It accepted the input feature class and executed successfully on its own. When I added it to the main model it now refuses to run, even with exactly the same input.


title: buffer of clusters in a sparse matrix

text: I work with MATLAB. I have a sparse matrix where I identified different clusters. The value within each cluster is equal, while each cluster has its own unique value. I have 0s as background (outside clusters). Here is an example with clusters 1 and 2:               A=(000002002000        110002222000        111000222200        110000022000        111000000000)      I'd like to use each cluster as "a polygon" and study the value of the outside neighbor pixels (a sort of buffer as in vector data). Obviously in the example it would output 0 as mean all the time, but the point is understanding how to do it, as I have to apply this to another matrix (I work with geolocated data, so I would use the buffer area to find mean values in specific rasters). Is there a way to do that? Also, if so, can I specify the width of this buffer (as number of pixels)?


title: Renaming ArcGIS Server Tile Cache to WorldWind Tile Cache

text: I have a bunch of tile cache created by ArcGIS Server - in WGS84. I would like to re-use these tiles in a WorldWind application. I am thinking all i need to do is rename ArcGIS Server tile cache directories and tile files to match the WorldWind specific structure. Has anyone done this before? Any tips/pointers of where I can start?


title: How to get distinct values from multiple columns in GeoTools?

text: Is there a way to get distinct values from multiple columns in geotools? I know I can get unique values from a column like this:               SimpleFeatureSource featureSource = dataStore.getFeatureSource(layer);     Query query = new Query(layer, Filter.INCLUDE, new String[] { "PROPERTY_A",      "PROPERTY_B", "PROPERTY_C" });     UniqueVisitor visitor = new UniqueVisitor("PROPERTY_A");     SimpleFeatureCollection collection = featureSource.getFeatures( query );     collection.accepts(visitor, null);     Set set = visitor.getResult()      But what I need is an equivalent to `select distinct PROPERTY_A, PROPERTY_B, PROPERTY_C`. **Edit:** I've tried to get unique values from the first column, and then create cql filters with the unique values from the first column to query unique values from the second column and so on, an then join it all in the end. But the performance, as expected, is very poor.


title: Sine inverse Calculation in FME

text: Please help to find formula for sine inverse calculation Safe soft fme. Is there any function in ExpressionEvaluator?. Thanks Venu


title: Shapefile to SVG with ArcMap (while retaining country attributes)

text: Does anybody know how to export the Natural Earth 1:110m shapefile to SVG with ArcMap while retaining the country attributes? From SVG I'm gonna try to convert it again over to RaphaelJS to achieve something like: http://backspace.com/mapapp/javascript_world/ http://backspace.com/mapapp/javascript_world/js/world_570.js But I'm new to all of this so I'm not sure how to go about it. Maybe there's no way but to manually add the country codes for the JS?


title: Qgis crashes when moving layer with many features

text: The qgis does not respond and crashes and I can not do what I need. I tried on minimum hardware requirements but have not found. I'm on a core i7 with 4gb memory windows 7, Any tips?


title: Desktop PC spec to run QGIS

text: I am helping a voluntary sector organisation in East London install QGIS to map 2011 Census data. I would be grateful for any advice on the spec of the Desktop PC we should purchase to run QGIS. Many thanks Jamil


title: How to change field length in QGIS

text: I have a shapefile I made in QGIS 1.7.4. One of the fields is a String of length 80. I'd like to make the field accept longer stings, but I can't figure out how to change its length. How do I do this? Thanks.


title: Geoserver 2.4.3 / OpenLayers - setting zoom level dynamically

text: Expected user behavior is that when user selects the County from a drop-down, the map should zoom to its max extents. I am using SQL View (Postgis on 9.3) and since I do not know the bounds, I gave the Max Bounds of the FL State. Not able to create layer without giving bounds. Can the boundaries be inferred from the SQL? The SQL for the view is as shown below:                  select ST_XMin(b.geom) as xmin,ST_YMin(b.geom) as ymin, ST_XMax(b.geom) as xmax,        ST_YMax(b.geom) as ymax, a.zcta5ce10, a.geom        from us_zip a, tl_2013_us_county b        where b.statefp='12'        and lower(b.name)='%COUNTY_NAME%'        and ST_contains(b.geom,a.geom) order by zcta5ce10      If it is not possible, I'm thinking of getting the max bounds for the selected county before invoking the OpenLayers map.


title: Removal of polygon overlaps

text: I'm using ARC View v10 and have a requirement to buffer 1500+ building polygons by 2m. When buffered an overlap is created as shown highlighted in green in the image below. Where as I would normally use a dissolve or merge to leave only the outline, I need to preserve the identity of each individual polygon after removal of the highlighted areas. Can anyone out there demonstrate how this could be achieved? ![Polygon Overlap](http://i.stack.imgur.com/jm3nk.jpg)


title: Insertion of a raster geotiff in PostgreSQL database

text: I am new in PostgreSQL and psql. What is the method to insert a raster geotiff in a postgreSQL database using command psql? What is the best format? Someone indicates to save a raster image as BLOB format.