title: How do I create a non-spatial web map

text: This question pertains to mapping arbitrary data using web map tile technology. Most of the time the data source for a web map is spatial in nature (OSM). I'd like to create a site map using the cartographic sense of the word. Rather than a boring list of links I'd like to model a site map after cartography from the Metro in Washington DC. This could be used as a sort of alternative to a standard navigation bar. Is it possible to churn any old image/vector into a set of tiles at different zoom levels, if I were to say, draw vectors manually with inkscape?


title: Which UTM zone 32N with prefix should I use in Rhineland-Palatinate (Germany)?

text: I would like to store and use GIS-Data in Rhineland-Palatinate (Germany). The coordinate System should be UTM 32N **with 32 as prefix** (false easting = 32,500,000; not 500,000). I looked at **EPSG 4647** [ETRS89 / UTM zone 32N (zE-N)].   But I'm confused because EPSG says: > Scope: Used in Schleswig-Holstein. > > Information Source: Landesvermessungsamt Schleswig-Holstein. Which UTM 32N coordinate system (EPSG-code) should use in Rhineland- Palatinate? EDIT: I will convert a large database from Gauss-Krueger to UTM. The database (and the new coordinate system) will be used for many years. We will exchange data with other (UTM with and without prefix, Gauss-Kueger). Since 4647 to refers to Schleswig-Holstein I am not sure whether 4647 is good for us. But there are no other EPSG codes with zone prefix.


title: How do I install ArcGIS 10.1 and use my current Python 2.7.2 install?

text: I just did a clean OS install and installed ActivePython 2.7.2 I want to use this Python install for use with ArcGIS 10.1 I just installed ArcGIS 10.1 and it did a new Python install. What do I need to do to get Python for ArcGIS to install on my current Python?


title: How best to fix a non-noded intersection problem in PostGIS?

text: I'm using a `PL/R` function and `PostGIS` to generate voronoi polygons around a set of points. The function that I am using is defined here. When I use this function on a particular dataset I get the following error message:               Error : ERROR:  R interpreter expression evaluation error     DETAIL:  Error in pg.spi.exec(sprintf("SELECT %3$s AS id,        st_intersection('SRID='||st_srid(%2$s)||';%4$s'::text,'%5$s')      AS polygon FROM %1$s WHERE st_intersects(%2$s::text,'SRID='||st_srid(%2$s)||';%4$s');",       :error in SQL statement : Error performing intersection: TopologyException: found non-noded      intersection between LINESTRING (571304 310990, 568465 264611) and LINESTRING (568465      264611, 594406 286813) at 568465.05533706467 264610.82749605528     CONTEXT:  In R support function pg.spi.exec In PL/R function r_voronoi      From examining this part of the error message:               Error performing intersection: TopologyException: found non-noded intersection between     LINESTRING (571304 310990, 568465 264611) and LINESTRING (568465 264611, 594406 286813)      at 568465.05533706467 264610.82749605528      This is what the problem listed above looks like: ![enter image description here](http://i.stack.imgur.com/MSxmL.png) I initially thought that this message might be caused by the existence of identical points, and tried to solve this using the `st_translate()` function, used in the following way:               ST_Translate(geom, random()*20, random()*20) as geom       This does fix the problem, but my concern is that I am now translating all points up to ~20m in the x/y direction. I also can't tell what an appropriate translation amount is need. For example, in this dataset through trial and error a `20m * random number` is ok, but how can I tell if this needs to be bigger? Based on the image above I think the problem is that the point is intersecting with the line while the algorithm is trying to intersect the point with a polygon. I'm not sure what I should be doing to ensure that the point is within a polygon, rather than intersecting with a line. The error is occurring on this line:               "SELECT        %3$s AS id,        st_intersection(''SRID=''||st_srid(%2$s)||'';%4$s''::text,''%5$s'') AS polygon      FROM        %1$s      WHERE        st_intersects(%2$s::text,''SRID=''||st_srid(%2$s)||'';%4$s'');"      I have read through this previous question, What is a "non-noded intersection"? to try to better understand this problem, and would appreciate any advice on how best to solve it.


title: Importing matplotlib basemap and shapely

text: I'm not sure if this is the right place to ask this, but here we go. I have a python script that uses matplotlib's basemap and another part that uses shapely to do an intersection of 2 polygons. If basemap is imported before shapely and I run the intersection I get this exception:                 intersect_poly = grid_poly.intersection(data_poly)     File "/sw/lib/python2.7/site-packages/shapely/geometry/base.py", line 334, in intersection       return geom_factory(self.impl['intersection'](self, other))     File "/sw/lib/python2.7/site-packages/shapely/topology.py", line 53, in __call__       "This operation produced a null geometry. Reason: unknown")     shapely.geos.TopologicalError: This operation produced a null geometry. Reason: unknown      If I import shapely first everything works fine. I would assume this is because of some "funkiness" in the way they are accessing the GEOS library. I've checked that in both situations the same library file is loaded in shapely ("print shapely.geos._lgeos"). Does anyone have an idea as to why this is happening and if there is a right way of doing this? Does this happen for anyone else? In the mean time I can just make sure to import shapely first (not sure if that affects basemap yet). Otherwise maybe I'll skim through the basemap source and ask this question on one of the mailing lists. Any suggestions on which mailing list I should post this to? Thanks. I'm using OSX(10.7) with a fink install that has "libgeos3.3.3-shlibs", "libgeos3.3.1-shlibs", "libgeos3.3.1", "libgeos3.3.0-shlibs", "libgeos3.3.0", and "shapely-py27 (1.2.16-1)" installed. And here's a simple test script that reproduces the problem (flip the imports and it works):               from mpl_toolkits import basemap     from shapely import geometry          g_ring = [(-88.462425, 26.992203), (-57.847187, 26.992203), (-57.847187, 17.599869), (-88.462425, 17.599869), (-88.462425, 26.992203)]     grid_g_ring = [(-123.044, 59.844000000000001), (-49.384999999999998, 57.289000000000001), (-65.090999999999994, 14.335000000000001), (-113.133, 16.369), (-123.044, 59.844000000000001)]          data_poly = geometry.Polygon(g_ring)     grid_poly = geometry.Polygon(grid_g_ring)          print grid_poly.intersection(data_poly).area


title: Layer visibility using maxscale minscale not working

text: I'm using Openlayer and Geoserver. I have a series of layers that have a scale parameter like: 1:50000-1:25000 where 50000 is the max scale while 25000 is the minscale. I have a map like in this example:               map = new OpenLayers.Map(     {             allOverlays: false,             projection: "EPSG:32632",             maxExtent: bounds,             scales:[500000,350000,250000,100000,25000,20000,15000,10000,5000],                     units: 'm'         });      And I have this layer WMS:                var layer = new OpenLayers.Layer.WMS("layer1",wms_address,             {                 layers:"ws:layer1",                 transparent:true                                         },             {             //        minResolution: 14.035643,     //        maxResolution : 7.01782148                     maxScale:50000,             minScale:25000         });       and I would expect that the layer within that scale it dosen't show up in the map. I need this kind of behaviour. However if I set the scale 1:25000, the layer is present inside the map (I saw it with Firebug) but it is not showed up in the map! I also tried with max and min Resolution but nothing, the layer is in the map but I cannot see it. What I'm doing wrong? Please help me, I looked everywhere in every forum I don't know how to solve this problem! Thank you very much. Laura


title: How to decrypt/decompile old toolbox tbx-file in ArcGIS 10

text: I have inherited an arcgis 10 toolbox tbx file for a web application. Are there any way to open the fil in for example Model Builder to see its interiors?


title: Leaflet: Linking a marker click event to an image

text: I'm breaking down my question into three parts. I am a beginner to Leaflet and to writing code. Currently, I have a JSON file with all of my markers linked to my map page. The idea is that every marker will have a different image so when the user clicks on the marker the thumbnail for the image will appear in the map, which will link to an image in the website. Here is an example of a marker in the JSON file: var mapshs4 = { "type": "FeatureCollection", "features": [ { "type": "Feature", "id": 0, "properties": { "idcode": "A", "popupContent": "Here!" + '< img src="thumbnail">' }, "geometry": { "type": "Point", "coordinates": [ -89.2, 22.4 ] } }   1. Is there a way to resize the popup box to the image in the popup of a marker.   2. At the same time, I would like the click event to connect to a larger sized image embedded within the site. My guess is that I would have to attach an ID to each image and to each marker to link them together.   3. Conversely, how do you make it so that when the user clicks on the image embedded in the site, it will go to that marker on the map? (For an example that uses a label: http://www.theglobeandmail.com/life/home-and-garden/real-estate/exclusive-how-much-are-homes-in-your-toronto-neighbourhood/article4520994/) Any help would be much appreciated.


title: How to represent two fields of the same states or district layer using different color borders and fills in QGIS?

text: If I had a set of data in QGIS that was gradually shading a state or district the exact way I wanted -- but I also need to portray more information about that state or district with additional colors, is that possible to do using colored borders? If so, can somebody walk me through that? And if not, has anybody else thought of another way to represent multiple things about one location? I guess the ideal would be something like shading the state lame yellow to bright yellow pending on how high their statistic was -- and if the state was Democratic-controlled, put a blue border around it, and if it was Republican- controlled, put a red border around it. Something like that. Any suggestions? Thanks! -- *Edit: I think I've potentially figured it out by adding another layer and creating my own border-only style and classify one as Democrat (blue) and one as Republican (red). If the polygon fills are 100% transparent and this layer was on top of the other one, I'd imagine that this should work, right? But every time I get to hitting "ok" to where the map should change to reflect the new borders, it always loads forever to the point where I have to 'Force Quit.' Any thoughts on that? Thanks again!


title: Fetch Google Map tile, without using Google Map API

text: Does anyone here know is there any other way to fetch Google map tile, without have to use Google Map API? Microsoft Bing Map does provide Rest API. wonder if possible to with Google Map, Cheer.


title: QGIS plugins with python dependencies

text: Problem: I have a Windows XP machine with Python 2.7.3 installed at ‘C:\Python27' and the scipy module installed at ‘C:\Python27\Lib\site- packages'. I installed QGIS 1.8 recently and started exploring plugins. One of the plugins I wanted to use required the scipy module, which, as mentioned above, was already installed. However, the plugin couldn't find scipy because it was not in the ‘C:\Program Files\Quantum GIS Lisboa\apps\Python27\Lib\site- packages’ folder. I considered installing a separate copy of scipy in ‘C:\Program Files\Quantum GIS Lisboa\apps\Python27\Lib\site-packages’ but backed away from using the scipy Windows installer when it couldn't see the QGIS python installation and didn’t give me the option to install scipy under QGIS. Work-around: So I went into the ‘C:\Program Files\Quantum GIS Lisboa\apps\Python27\Lib’ folder and found the 'site.py' file. I opened it in a text editor and, right below the import statements, I found the text: PREFIXES = [sys.prefix, sys.exec_prefix] and I changed it to read: PREFIXES = [sys.prefix, sys.exec_prefix, 'C:\Python27\Lib\site-packages'] Now the QGIS python installation is able to find my original scipy module and it seems to be working fine. Presumably, QGIS plugins will be able to take advantage of the many other [python(x,y)] modules I have installed in this directory. Hope this is useful for others.


title: Symbology Stretch Multiple Raster

text: I need to change stretch of a raster using properties - Symbology to **Minimum-Maximum** like image bellow :   ![enter image description here](http://i.stack.imgur.com/NhuoM.png) Problem is I have about 172 raster, is there any way to make this automaticaly change? I have arcinfo license


title: Model Builder is not intersecting erased and selected files

text: I'm having trouble debugging my model that I'm creating in ArcMap 10.1 ModelBuilder. The model is just a simple string of geoprocesses, but there's a point at which the model isn't exporting the expected features. Specifically, I've narrowed it down to a point at which the model takes a polyline file which had been erased (interim_erased) and intersects it with a selected polygon file (Inactive_ESA), and at this point the output file (interim_intersect) is empty. The only thing I see that could be happening is that the "ESA" isn't getting selected at the "Select by Attributes" step (The output being "Inactive_ESA"), because the erased file "interim_erased" has the expected features in its table. Also, I've double-checked the selection expression - it's correct and works fine outside of ModelBuilder. I also have the license to execute Erase. Maybe, if someone knows how to look into ModelBuilder at what's actually getting selected at the "Select Layer By Attributes" point in the model? ![part of the model process which has output a blank attribute table](http://i.stack.imgur.com/s5X2S.png) ++++++++++++ ADDED 2014-07-31 +++++++++++++++ The ESA polygons are broken down by active(green) and inactive(red), the polylines as well, active(green) and inactive(red). I'm selecting the green polygons in the ESA, and selecting the red polylines. Then I'm Erasing the red polylines with the green ESA. (This part works, and my intermediate data for this process outputs the expected result - red polylines only on the red ESA with no red polylines overlapping onto green ESA). However, when this erased file of red polylines is intersected with my selected Inactive_ESA (which is where I'm selecting the red polygons), the table is empty. So this suggests to me that the ESA isn't actually getting selected by the "Select Layer By Attribute (9)" process. How might I see what's actually getting selected at this step? (because I've double checked my expression and it's correct and works fine outside of the model) ![ESA and polylines included in the model](http://i.stack.imgur.com/PpXp0.png)


title: How to use Ordinary Least Square (OLS) & Geographically Weighted Regression (GWR)?

text: I have managed to join the data set(xl. file) with a shapefile. I want to run these models, OLS and GWR in Arcmap, my problem is, I'm bit confuse about these fields, "Input Feature Class", "Output Feature Class", and "Input Features"?, what do I need to do in order to achieve this task?


title: How to mask the images in GRASS

text: I have MODIS NDVI images, which needs to be masked with VI_Quality. I am using GRASS software and executing GRASS commands using shell scripts. I need to perform raster calculation (multiplication) for NDVI images with the VI_Quality. I want to multiply ith NDVI image with the ith image of VI_Quality. I have used the following codes :               for file in $NDVIFILES;     do     r.mapcalc "$file = $NDVIFILES * $QAFILES;     done          i=1     for file in $NDVIFILES;     do     r.mapcalc "file $((i++)) = NDVIFILES $((i++)) * QAFILES $((i++))";     done      I am unable to do this raster multiplication in a loop. Please let me know the correct syntax to do this operation. Thanks


title: Adding KmlLayer using local file (File API)

text: I am trying to load a local KML file to the map for this I am using FileReader object (File API). To use google.maps.KmlLayer I need to use readAsDataURL to transform to a URL the file, the problem is that KML File is not loading to the map and I don't get any error. I think may be url attribute in KmlLayer doesn't accept data url. This is the code:               gp.addEventListeners = function(){       $("#file-upload").change(function (evt){           var file  = evt.target.files[0]; // FileList object           var reader = new FileReader();                reader.onloadend = function () {              var ctaLayer = new google.maps.KmlLayer({               url: reader.result             });             ctaLayer.setMap(gp.map);           }                if (file) {             reader.readAsDataURL(file);           }            });     }      In OpenLayes the same code is working. What I can do? I need help :( Sorry for my English.


title: Anyone who has the ESRI Census 2000 Tiger/Line data?

text: The ESRI Census 2000 TIGER/Line Data were free. http://www.esri.com/data/download/census2000-tigerline http://arcdata.esri.com/data/tiger2000/tiger_download.cfm But they vanished. The download link is broken. Is there anyone out there who has the data?


title: ArcGIS 10 - Text background not showing up when exporting

text: I have highways labeled with the standard shield and route number on it. For my main data frame the labels come up correctly. I also have 4 inset maps where **only the number will show up, not the shield** when exporting to a file (any type). Here are some things that I tried that did not work:   * Exporting the labels to annotation   * Creating replacement graphics (in each inset's data frame view)   * For each of the above, using .png background instead of .emf The only solution that seems to work is creating graphics as part of the document itself and overlaying those on top of the inset. This is very bad practise however and the graphics will not pan with the inset. Just to clarify, the shields (text background) shows up in the layout screen when editing the map but not in the exported product after exporting the map to a .pdf or other extension. The version of ArcGIS being used is 10.1 SP1. I cannot upgrade or change the installation. Suggestions tried:   * Deleting the main dataframe, then exporting - same result   * Moving the inset dataframe to a new document **UPDATE:** It appears that ArcGIS will hide (or not draw) background symbols that are smaller than a certain size compared to the map (% or ratio). When I made the print size smaller (it was originally 34x40 inches) or the symbols larger, the shields showed up. Does anyone know how to change this option on ArcGIS and possibly show all symbols regardless of size? Following comment suggestions, the links provided in this forum post did not help.   * http://support.esri.com/en/knowledgebase/techarticles/detail/28548   * http://support.esri.com/en/knowledgebase/techarticles/detail/30753   * http://forums.esri.com/Thread.asp?c=93&f=1731&t=235317   * http://forums.esri.com/Thread.asp?c=93&f=1731&t=178684   * http://support.esri.com/en/knowledgebase/techarticles/detail/17783 However, following the last link showed that the exported EMF file was not correct, signaling that "the problem exists within ArcMap, usually in converting the elements drawn in the map to the EMF format."


title: WMS GetFeatureInfo and CQL on Multiple Layers

text: I am having issues with querying CQL filtered layers with OpenLayers and GeoServer, using the GetFeatureInfo control. It is ignoring the CQL filter and just querying the entire layer (as evidenced by both the results and the request it sends to GeoServer). I've seen examples where the CQL filter is passed to GeoServer along with the GetFeatureInfo request, but I have several layers visible and each of them has a different CQL filter; so I'm not sure how/if I can handle this by appending my filter to the request URL. Does anyone have any ideas how I can solve or work around this issue?


title: How to handle a UnicodeEncodeError with sampling tool?

text: I'm using qgis and I keep getting an error when using the point sampling tool, its keeping me from getting the sampling done. Does anyone has an idea what is it about? Thanks! ![error](http://i.stack.imgur.com/KquJW.png)


title: Looking for geodata of Europe's major roads

text: As stated in the title, I'm currently looking for a shapefile/other geodata source of Europe's major roads (highways). If anybody knows a free datasource (or a possibility to extract that data from Open Street Map without processing the world file), it would be greatly appreciated.


title: Merging many to one data rows in arcpy

text: The dataset 'left_right_merge' at the bottom of my screen shot shows the dataset that I'm currently working with. The NODE values indicate a single point feature, but is duplicated with many attributes. I want to collapse the dataset so that it looks like the 'test_data_nodes' at the top of the screen shot. How would I do this in python? ![enter image description here](http://i.stack.imgur.com/eFsrZ.jpg)


title: ArcGIS SP4 fails to open file geodatabases created with ogr2ogr

text: I've been using the new FileGDB support in gdal-1.9.0 to create file geodatabases that can be opened in ArcGIS. For example (I'm doing this all on one line, so I'm not certain this will run as indented here):               ogr2ogr -overwrite -f 'FileGDB' \         -sql "SELECT                    CAST(gml_id AS character(255)), id,                    CAST(scene_id AS character(255)),                    CAST(license_type AS character(255)),                    CAST(url AS character(255))               FROM SDMI_Inventory" \         sdmi_inv.gdb \         WFS:"http://wfs.alaskamapped.org/inventory"      This worked great until 10.0 Service Pack 4, at which point ArcMap / ArcCatalog started reporting errors similar to this when you try to add / preview the layer:               Error opening feature class     The item does not have a definition. [SDMI_Inventory]     The item does not have a definition.     The XML being loaded could not be parsed. [physical name: SDMI_INVENTORY]     The item does not have a definition.     The XML being loaded could not be parsed. [physical name: SDMI_INVENTORY]      I've got the latest (1.2) API and the latest GDAL (1.9.0), x86_64 GNU/Linux. Anyone have any ideas how to get this to work again?


title: Geoserver SQL Parameter View

text: I'm trying to use Geoserver Parameter SQL views, but I'm having trouble getting geoserver to recognize the parameters (and I'm not able to follow the documentation - it's not very clear to me). Here's the code I'm working with so far, does anyone have any suggestions?               SELECT *     FROM Table1 INNER JOIN Table2     ON Table1.ID = Table2.ID     WHERE Table1.Date = %Date%      I'm trying to get geoserver to recognize each "Date" as a parameter - there are about 25 total dates, and I want the user to be able to select a date (or multiple dates) for display. The dates are in yyyy-mm-dd format. Thanks for your help.


title: Hibernate Spatial 4 and Postgis 2.0

text: I'm having some issues integratings this technologies:   * Hibernate Spatial 4.0-M1   * Postgis 2.0.2 (with compiled jdbc 2.0.2)   * Hibernate 4.1.1 The specific error is:               Caused by: org.postgresql.util.PSQLException: Can't infer the SQL type to use for an instance of org.postgis.PGgeometry. Use setObject() with an explicit Types value to specify the type to use.      The entity anotation is:               @NotNull     @Column(columnDefinition="Geometry")     @Type(type="org.hibernate.spatial.GeometryType")     private Point geom;      And the Point creation example is:               Location location = new Location();     WKTReader fromText = new WKTReader();     Point geom = null;     try {         geom = (Point) fromText.read("POINT(-56.2564083434446 -34.8982159791812)");     } catch (ParseException e) {         throw new RuntimeException("Not a WKT string:" + "SRID=4326;POINT(-56.2564083434446 -34.8982159791812)");     }     if (!geom.getGeometryType().equals("Point")) {         throw new RuntimeException("Geometry must be a point. Got a " + geom.getGeometryType());     }     location.setGeom(geom);     locationDAO.insert(location);      Thanks!


title: Is there any way to use mbutil on windows?

text: I have used mbutil on ubuntu. Its a very handy tools. I am curious that is there any way to use it on windows 64 bit machine. If there is a way please help me to do so. Thanks in advance.


title: Ogr2Ogr SQLite Geometry Functions

text: I believe title is all clear.I'm using Ogr2ogr utilities & binaries to read geometry. I used something like this:               ogr2ogr -f "SQLite" test2.sqlite test.sqlite -sql "SELECT ST_X(GEOMETRY),ST_Y(GEOMETRY),COL1,COL2 FROM TABLE"      Error says that there is no function ST_X,where I also tried X(..) that doesnt work either. What am I missing ? I checked the documentation for SpatialLite.There is no other inline function that export X,Y portions of geometry to double.I also need to create geometry from LAT,LONG columns from selection. How can i do these ?


title: ogrListLayers() in rgdal returns character(0) for a kml file

text: I try to read a kml file in R. In the past, there was a post (link). I checked stack overflow and found similar posts. By following the examples, I was trying to read a kml file. I could read the kml file using `getKMLcoordinates()` in the maptool package and extract geo data (i.e., longitude and latitude). I also used the `rgdal` package in order to see layer names. But, `ogrListLayers(dsn = "name.kml")` surprisingly returned character(0). My questions are:   1. Is is possible for kml files not to have layer names?   2. If there is no layer name, is there still any hope to find the names? Or do you have to simply give up? P.S. I originally posted this message in stack overflow, but I have not received any suggestions. So I decided to post this message here.


title: What does "Error: Layer not found: -1" from ArcGIS Silverlight API 2.4 mean?

text: We are using ESRI SL 2.4 API- and our application has one dynamic service and one feature layer. Only few of the users getting below error, not sure whether this is issue with Microsoft sliverlight version. Any quick help on this. Error: Layer not found: -1 at ESRI.ArcGIS.Client.Layer.OnInitializationFailed(EventArgs e) at ESRI.ArcGIS.Client.Layer.Initialize() at ESRI.ArcGIS.Client.FeatureLayer.info_Error(Object sender, TaskFailedEventArgs e) at ESRI.ArcGIS.Client.FeatureService.FeatureLayerInfo.OnError(Exception err) at ESRI.ArcGIS.Client.FeatureService.FeatureLayerInfo.Execute_Completed(Object sender, DownloadStringCompletedEventArgs e) at System.Net.WebClient.OnDownloadStringCompleted(DownloadStringCompletedEventArgs e) at System.Net.WebClient.DownloadStringOperationCompleted(Object arg)


title: Migrating Lines to ArcGIS polyline

text: I'm trying to import geometry lines from Oracle Spatial database to ArcGIS enterprise geodatabase. and some data will cause the error. Its says not enough points for feature class, but when I look at offending object. it clearly has 2 or more coordinates stored, but I notice the start and end point is the same. I have check the following a) Using geometry filter to take out lines only b) Use FME inspector to view the lines. b) Used geometry validator 4) Used duplicate coordinate -> Valid dataset will cause the same error How can I solve this problem? May i know is it possible to set FME to skip those data


title: Easy Custom Labeling plugin does allow to save?

text: I use the **Easy Custom Labeling plugin** chose positions of the labels and safe the layer and project for continuing later. When opening again the project, the layer opens but without any objects inside. Is there a reason why the layer is saved, but not its content?


title: How to select behind feature which is overlapped with the same layer on QGIS 1.8.0

text: My issue: I received data which was created without "Avoid intersection of new polygon". There a few overlapped feature with the same layer. I could fine overlapped feature by the spatial query. But then, I could not select behind small feature. At present, I select target feature from attribute table and delete that feature. I'm hoping I could be able to select feature without attribute table. Grateful for any help!!


title: Neighborhood analysis based on occurrence of value

text: How can I do, either in ArcGIS, grass or another tool, the following: for each cell (pixel), if there is at least one cell with value '4' in a neighborhood of 5 cells around it, reclassify that cell to 0, otherwise don't change the value? Thanks.


title: Weather conditions forecasting / interpolation

text: I'm trying to carry out weather conditions interpolations. I have information of 92 weather stations for 800,000 km2. As you can imagine the extension of land is huge and I need some kind or weather forecasting models to predict wind velocity and direction and solar irradiation for points in between the stations via interpolations. This means if I have a weather station in A and another one in C, I can know what is happening in B with weather forecasting models. Do you know if QGIS does something alike? (maybe there is a plugin to download ?!) Are you familiar with any models that could help me out? Thanks in advance!!


title: Should I move from ArcGIS modelbuilder to python?

text: I have been learning how to use ModelBuilder with ArcGIS 10.2, and I like that it shows me visually how to process data using ArcGIS's functions. I think visually. I see, however, many replies by folks that users should drop MB in favour of Python. I am not a programmer, and don't foresee dropping everything to learn a new language. It's not that I can't; it's more that I don't have the time to learn a new language. Why should I, a non-programmer, move from ModelBuilder to Python, and if I really should, is there a relatively painless way to make the transition? Is there a visual Python equivalent to MB, for example? This trend of pushing people towards Python and away from MB troubles me. I have noted that development for MB seems at a standstill, ditto bug-swatting for MB. Can anyone shed some light on a murky topic for me? Thx, Russell in Canada


title: How to search and zoom to result using WMS?

text: I am making a web gis application (openlayers/geoserver) to display a layer of properties, I need to add a functionality of search in it based on the name of the owner. A list of similar names owners with address should be displayed and when i click a particular item the map should be zoomed on the particular property location. The problem is that i want this functionality using WMS. (I do not want it in WFS as the vector layer can be copied using QGIS or other tools, it is not secure). Please let me know how will this be possible using WMS.


title: error performing wfs-t insert in postgres 9.1.4

text: I am getting **error performing insert** exception while performing wfs-t insert operation when i use **postgres 9.1.4 and postgis 2.0**. but i get **successful insert** when using **postgres 9.0.1 and postgis 1.5** (on another computer). is this some kind of bug in postgres/postgis or i am doing something wrong ? my script is below               OpenLayers.ProxyHost= "/cgi-bin/proxy.cgi?url=";     //variables     //var host = "http://192.168.10.11";     var host = "http://localhost";     var port = ':8888';     var centerX = 85.3403;//491213.721224323//-123.1684986291807;//9497800;     var centerY = 27.70774;//5456645.24607268//49.245339757767844;//3212000;     var center = new OpenLayers.LonLat(centerX, centerY);     var range = 1.5;//10000000//.5;//10000;     var map_bound = [centerX-range,centerY-range,centerX+range,centerY+range];     var map_wd = 490;     var map_ht = 360;     var legend_wd = 100;     var legend_ht = map_ht;     var tbar_wd = map_wd + legend_wd;     var tbar_ht = 100;     var zoom = 15;          var proj_sp = new OpenLayers.Projection("EPSG:4326");          var DeleteFeature = OpenLayers.Class(OpenLayers.Control, {         initialize: function(layer, options) {             OpenLayers.Control.prototype.initialize.apply(this, [options]);             this.layer = layer;             this.handler = new OpenLayers.Handler.Feature(                 this, layer, {click: this.clickFeature}             );         },         clickFeature: function(feature) {             // if feature doesn't have a fid, destroy it             if(feature.fid == undefined) {                 this.layer.destroyFeatures([feature]);             } else {                 feature.state = OpenLayers.State.DELETE;                 this.layer.events.triggerEvent("afterfeaturemodified",                      {feature: feature});                 feature.renderIntent = "select";                 this.layer.drawFeature(feature);             }         },         setMap: function(map) {             this.handler.setMap(map);             OpenLayers.Control.prototype.setMap.apply(this, arguments);         },         CLASS_NAME: "OpenLayers.Control.DeleteFeature"     });          function oldBrowser() {         var browser = navigator.appName;         var b_version = navigator.appVersion;         var version = parseFloat(b_version);         var useragent = navigator.userAgent;         switch (browser){             case 'Microsoft Internet Explorer':                 browser = "MSIE";                 version = useragent.substr(useragent.lastIndexOf('MSIE') + 5, 3);                 if(version<10)                     return true;                 break;             case 'Netscape':                 if (useragent.lastIndexOf('Chrome/') > 0) {                     browser = "Chrome";                     version = useragent.substr(useragent.lastIndexOf('Chrome/') + 7, 10);                 }                 else if (useragent.lastIndexOf('Firefox/') > 0) {                     browser = "Firefox";                     version = useragent.substr(useragent.lastIndexOf('Firefox/') + 8, 5);                 }                 else if (useragent.lastIndexOf('Safari/') > 0) {                     browser = "Safari";                     version = useragent.substr(useragent.lastIndexOf('Safari/') + 7, 7);                 }                 else                 { alert("undefined browser"); }                 break;             case 'Opera':                 version = useragent.substr(useragent.lastIndexOf('Version/') + 8, 5);                 break;         }     }          Ext.onReady(function() {             if(oldBrowser()){             alert("Your browser version is pretty old."                  +"please update to latest version"                  +"performance may be reduced.");         }              var saveStrategy = new OpenLayers.Strategy.Save();         /*saveStrategy.events.register("start", '', alert('saving'));         saveStrategy.events.register("success", '', alert('save success'));         saveStrategy.events.register("failure", '', alert('save failure'));         */         //map configuration         var map_options = {             allOverlays:true,             maxExtent:new OpenLayers.Bounds(map_bound[0],map_bound[1],map_bound[2],map_bound[3]),             numZoomLevels: zoom,             controls: [new OpenLayers.Control.PanZoomBar(),                 new OpenLayers.Control.MousePosition(),                 new OpenLayers.Control.ScaleLine(),                 new OpenLayers.Control.Scale(),                 new OpenLayers.Control.Navigation()             ],             minScale: 50000,             projection: "EPSG:4326"         }         var map = new OpenLayers.Map(map_options);         var navi = new OpenLayers.Control.NavigationHistory();         map.addControl(navi);         center.transform(proj_sp, map.getProjectionObject());         //layers configuration         var world = new OpenLayers.Layer.WMS(             "OpenLayers WMS"             ,"http://localhost:8888/geoserver/wms"             ,{layers: 'nurc:Arc_Sample'}             ,{displayInLayerSwitcher:false}         );              var roadswgs = new OpenLayers.Layer.WMS(             "roadswgs"             ,host+port+"/geoserver/wms"             ,{  layers: "postgis:trans_ln_8"                 ,transparent: true             }             //,{'reproject': true}         );              var portwgs = new OpenLayers.Layer.WMS(             "airportwgs"             ,host+port+"/geoserver/wms"             ,{  layers: "postgis:airport_runway"                 ,transparent : true             }             //,{'reproject': true}         );              //var vector = new OpenLayers.Layer.Vector('editable');              var vector = new OpenLayers.Layer.Vector('editable'             ,{  strategies: [new OpenLayers.Strategy.BBOX(),saveStrategy]                 ,protocol: new OpenLayers.Protocol.WFS({                     url:  host+port+"/geoserver/wfs",                     featureType: "airport_runway",                     featureNS: "http://geoserver/postgis",                     version: "1.1.0",                     geometryName: "geom"                 })                 /*,styleMap: new OpenLayers.StyleMap({                     pointRadius: 5,                     strokeColor: "#ff005b",                     strokeWidth: 1,                     fillOpacity: 0.4,                     fillColor: '#ffffff'                 })*/             }         );         alert('vector.protocol.url='+vector.protocol.url);         var wfsroads = new OpenLayers.Layer.Vector("PostgisWFS",              {   strategies: [new OpenLayers.Strategy.BBOX()],                 protocol: new OpenLayers.Protocol.WFS({                     url:  host+port+"/geoserver/wfs",                     featureType: "trans_ln_8",                     featureNS: "http://geoserver/postgis",                     version: "1.1.0"                                     })                 /*,styleMap: new OpenLayers.StyleMap({                     pointRadius: 6,                     strokeColor: "#ff005b",                     strokeWidth: 1,                     fillOpacity: 0                 })*/             }         );         map.addLayers([wfsroads,roadswgs,vector,portwgs/*world,,gsat*/]);                //adding OSM layer         //map.addLayer(new OpenLayers.Layer.OSM("OSM"));         //adding google satellite layer         map.addLayer(new OpenLayers.Layer.Google("Google Satellite",{type:G_SATELLITE_MAP, numZoomLevels: 22}));                        //mappanel configuration         var mapPanel = new GeoExt.MapPanel({             renderTo: 'gxmap',             height: map_ht,             width: map_wd,             map: map,             title: 'A Simple GeoExt Map',             extent: map_bound,             center: center,             //zoom: zoom         });         //layer tree configuration         var layerList = new GeoExt.tree.LayerContainer({             text: 'All Layers',             layerStore: mapPanel.layers,             leaf: true,             expanded: true         });              var layerTree = new Ext.tree.TreePanel({             title: 'Layers',             renderTo: 'layerTree',             root: layerList,             width: legend_wd,             //height: '*',             autoheight: true         });              //legend configuration         var legendPanel = new GeoExt.LegendPanel({             title:'Legend',             layerStore: mapPanel.layers,             renderTo: "legend",             border: true,             width: legend_wd,             //height: legend_ht,             autoheight: true         });              //editor tools         var editor = "editor tools";         { //snap for drawing          var snaproad = new OpenLayers.Control.Snapping({             layer: vector,             targets: [{                 layer: wfsroads,                 tolerance: 15,                 greedy: true             }]         });         //snaproad.activate();         var snapvector = new OpenLayers.Control.Snapping({             layer: vector,             targets: [{                 layer: vector,                 tolerance: 15,                 greedy: true             }]         });         snapvector.activate();}         { //Draw Feature control         var draw = new OpenLayers.Control.DrawFeature(             vector,     //layer             OpenLayers.Handler.Path, //Handler             {   multi:true,                 type: OpenLayers.Control.TYPE_TOGGLE,                 //handlerOptions : {maxVertices: 5}             }//options         );              /*draw.featureAdded = function(feature) {             //feature.layer.eraseFeatures([feature]);             // cast to multilinestring             //alert(feature.geometry);             //feature.geometry = new OpenLayers.Geometry.MultiLineString(feature.geometry);             //alert(feature.geometry);//feature.style.strokeColor = "#ff0000";             feature.state = OpenLayers.State.INSERT;             feature.layer.drawFeature(feature);         }*/         map.addControl(draw);         // DrawFeature control, a "toggle" control         var draw_ac = new GeoExt.Action({             text: "draw",             control: draw,             map: map,             // button options             enableToggle: true,             tooltip: "draw feature",             toggleGroup: editor         });}         { //Select Feature control         var select = new OpenLayers.Control.SelectFeature(             vector,             {multiple: false, hover: true,             }         );         var select_ac = new GeoExt.Action({             text: "select",             control: select,             map: map,             // button options             enableToggle: true,             tooltip: "modify feature",             toggleGroup: editor         });         }         { //Modify Feature control         var modify = new OpenLayers.Control.ModifyFeature(vector,{type: OpenLayers.Control.TYPE_TOGGLE});         map.addControl(modify);         var modify_ac = new GeoExt.Action({             text: "modify",             control: modify,             map: map,             // button options             enableToggle: true,             tooltip: "modify feature",             toggleGroup: editor         });         }         { //Delete Feature control         var del = new DeleteFeature(vector, {title: "Delete Feature"});         map.addControl(del);         var del_ac = new GeoExt.Action({             text: "delete",             control: del,             map: map,             // button options             enableToggle: true,             tooltip: "delete feature",             toggleGroup: editor         });         }         { //Save Feature control         var save_ac = new GeoExt.Action({             text: "save",             //icon: 'icons/icon_zoomnext.png',             //control: saveStrategy.save(),             onClick: function(){                 //var s=confirm('save?');                 //if(s==true){                     saveStrategy.save();                 //  alert('saved');                 //}                  //else                 //  alert('not saved');             },             disabled: false,             //deactivateOnDisable: true,             //activateOnEnable: true,             tooltip: "save road"         });}         { //Measure Feature control         var length = new OpenLayers.Control.Measure(OpenLayers.Handler.Path, {             eventListeners: {                 measure: function(evt) {                     alert("The length was " + evt.measure + ' ' + evt.units);                 }             }         });         var area = new OpenLayers.Control.Measure(OpenLayers.Handler.Polygon, {             eventListeners: {                 measure: function(evt) {                     alert("The area was " + evt.measure + ' sq ' + evt.units);                 }             }         });         var measu_length = new GeoExt.Action({             //text: "measu_length",             icon: 'icons/icon_ruler.gif',             control: length,             map: map,             // button options             enableToggle: true,             tooltip: "measu_length",             toggleGroup: editor         });         var measu_area = new GeoExt.Action({             //text: "measu_area",             icon: 'icons/tool_measurearea.png',             control: area,             map: map,             // button options             enableToggle: true,             tooltip: "measu_area",             toggleGroup: editor         });}         { //identify Feature control         var info = new OpenLayers.Control.WMSGetFeatureInfo({             url: host+port+'/geoserver/wms',              eventListeners: {                 getfeatureinfo: function(event) {                     map.addPopup(new OpenLayers.Popup.FramedCloud(                         "chicken",                          map.getLonLatFromPixel(event.xy),                         null,                         event.text,                         null,                         true                     ));                 }             }         });         var identify_ac = new GeoExt.Action({             //text: "identify",             icon: 'icons/icon_query.gif',             control: info,             map: map,             // button options             enableToggle: true,             tooltip: "identify",             toggleGroup: editor         });}         { //pan Feature control         var pan = new OpenLayers.Control.DragPan();         var pan_ac = new GeoExt.Action({             //text: "pan",             icon: 'icons/tool_pan.gif',control: pan,             map: map,             // button options             enableToggle: true,             tooltip: "pan map",             toggleGroup: editor         });         }         {//tools simple         var toolsbar = "tools"          //tools toolbar         //zoom to max extent tool         var maxext = new GeoExt.Action({             //text: "Max Extent",                    icon: 'icons/icon_fullextent.gif',             control: new OpenLayers.Control.ZoomToMaxExtent(),             map: map,             deactivateOnDisable: true,             activateOnEnable: true,             tooltip: "zoom to max extent"         });              var zoomin = new GeoExt.Action({             //text: "Zoom In",             icon: 'icons/icon_zoomin.gif',             control: new OpenLayers.Control.ZoomIn(),             map: map,             deactivateOnDisable: true,             activateOnEnable: true,             tooltip: "zoom inside",         });              var zoomout = new GeoExt.Action({             //text: "Zoom Out",             icon: 'icons/icon_zoomout.gif',             control: new OpenLayers.Control.ZoomOut(),             map: map,             deactivateOnDisable: true,             activateOnEnable:true,             tooltip: "zoom outside",         });              var pan = new GeoExt.Action({             text: "Pan",             control: new OpenLayers.Control.Pan(),             map: map,             deactivateOnDisable: true,             activateOnEnable:true,             tooltip: "drag to move",         });              var prev_ext = new GeoExt.Action({             //text: "previous",             icon: 'icons/icon_zoomprevious.png',             control: navi.previous,             disabled: false,             deactivateOnDisable: true,             activateOnEnable: true,             tooltip: "previous in history",         });              var next_ext = new GeoExt.Action({             //text: "next",             icon: 'icons/icon_zoomnext.png',             control: navi.next,             disabled: false,             deactivateOnDisable: true,             activateOnEnable: true,             tooltip: "next in history"         });              }         var toolbar = new Ext.Toolbar({                 //embed tool in toolbar             renderTo: 'tools',             width: tbar_wd,             height: tbar_ht,             autoWidth : true,             autoHeight : true,             items: [draw_ac,'-',modify_ac,'-',del_ac,'-',save_ac,'-',measu_length,'-',measu_area,'-',select_ac,'-',identify_ac,'-',pan_ac,maxext,'',zoomin,'',zoomout,'',prev_ext,'',next_ext]         });              // the part below is just to make the bounds show up on the page         function updateBounds(){             var bbox = map.getExtent().toArray();             code = "Bounds      lower left: "+bbox[0]+','+bbox[1]+"     upper right"+bbox[2]+','+bbox[3];             document.getElementById('bounds').innerHTML = code;         }         // update the bounds with each map move         map.events.register('moveend', map, updateBounds);         // and update the bounds on first load         updateBounds();          });


title: How to pan (rather than zoom) when using a 2-finger drag on a Mac

text: Using Mac OSX, dragging the trackpad with 2 fingers performs a scroll of the web page, equivalent to using the mouse wheel. However, if the cursor is located within an OpenLayers map, the 2-finger drag performs a zoom-in on the map. Is it possible to change this action to perform a pan instead of a zoom? Eg, see the ArcGIS Server JS API's behaviour on a Mac with a trackpad. Thanks


title: I am trying to put autocomplete of Google API V3 for places, but giving some error

text: I am trying to put textbox within the infowindow which is working fine in order to get direction and route. And now i want to make textbox autocomplete for google api places. I have ref this link "autocomplete in google api for places":https://developers.google.com/maps/documentation/javascript/places#places_autocomplete I am trying to put autocomplete in google API V3 but some error of **a is null {main,places}.js** following are my code: Find way #directions-panel { height: 100%; float: right; width: 390px; overflow: auto; }                 #map-canvas {         margin-right: 400px;       }            #control {         background: #fff;         padding: 5px;         font-size: 14px;         font-family: Arial;         border: 1px solid #ccc;         box-shadow: 0 2px 2px rgba(33, 33, 33, 0.4);         display: none;       }            @media print {         #map-canvas {           height: 500px;           margin: 0;         }              #directions-panel {           float: none;           width: auto;         }       }     </style>                Start: End: --> can any one please tell me weather i am going right. and what correction i need to do for same.


title: How to insert non-spatial data through WFS-T insert?

text: When i perform wfs-t insert operation using postgis database only spatial data and primary key gets inserted. Can anyone point me to some examples in which attribute data are also inserted along with spatial data using openlayers?


title: Data defined label alignment option greyed out on new Quantum GIS Dufour

text: I am trying to customise the positioning of labels in relation to points in a vector file layer in the new QGIS 2.0. In QGIS 1.8 I was able to enter two customised X/Y columns in the attribute table of a shapefile to precisely align the labels around the points. However, from the new QGIS 2.0 menu 'Layer label options'>'Placement'>'Data defined', the X and Y 'Alignment' options are greyed out and I am unable to enter the columns. Please see a screen dump of the issue below Can anybody explain why this might be the case? Thanks for your consideration. ![enter image description here](http://i.stack.imgur.com/0eiNe.jpg)


title: how do you identify ridges on an elevation raster?

text: I was wondering whether anyone might be able to tell me in layman's terms how to extract or differentiate between ridges and hilltops from 10m elevation data using ArcGIS 10.1? Any advice would be greatly appreciated


title: Is it possible to convert a Geoserver SLD to Quantum QML format?

text: I have a web application which runs on Geoserver. I am migrating my application to desktop and am struck at applying the same color schemes and symbology. I am not able to find a way to convert SLDs to QML for implementing in desktop application. I am aware that Quantum GIS 1.8 has the SLD input in styling. I specifically need a solution for Quantum GIS 1.6, as it is more preferred in my organization. Any suggestion is highly appreciated.


title: SAGA grid VRT, file how to apply a LUT?

text: I've created a heap of saga grid files for some hydro survey work. I'm looking at a way to read these into a mapguide raster/grid layer, I'm able to load up the grids using a vrt file, I can't figure out how to apply a color table to them though, values should bein the order of 0 to 50 float ... So I guess I need help on stylizing the source using a vrt file?? And perhaps some other program external to mapguide to check the display of the vrt dataset? Any help much appreciated, just ask if theres more info I can supply. an example of the vrt source I'm playing with               <VRTDataset rasterXSize="2957" rasterYSize="4598">               <GeoTransform> 5.3048100000000000e+005, 5.0000000000000000e-001, 0.0000000000000000e+000, 6.9916320000000000e+006, 0.0000000000000000e+000,-5.0000000000000000e-001</GeoTransform>       <VRTRasterBand dataType="Float32" band="1">         <NoDataValue>-9.99990000000000E+004</NoDataValue>          <ComplexSource>       <SourceFilename relativeToVRT="1">15251007.sdat</SourceFilename>       <SourceBand>1</SourceBand>       <SourceProperties RasterXSize="2957" RasterYSize="4598" DataType="Float32" BlockXSize="2957" BlockYSize="1" />       <SrcRect xOff="0" yOff="0" xSize="2957" ySize="4598" />       <DstRect xOff="0" yOff="0" xSize="2957" ySize="4598" />       <LUT>0:50,15:256</LUT>             <NODATA>-99999</NODATA>     </ComplexSource>


**One of my WMS layer is not loading the image tiles.** See the code:                   <script type="text/javascript">         var map, ia_wms;         function init(){             map = new OpenLayers.Map('map', {zoomDuration: 1,projection: 'EPSG:3857'});         var esr_all = new OpenLayers.Layer.WMS("ESR_Boundary","http://localhost/cgi-bin/mapserv?map=/var/www/mstest/plant_sr.map",{'layers':"ESR_outline",'transparent':true,'format':'image/png'}                ,{                 isBaseLayer: false,opacity: 100                 }                 );         var mslayer = new OpenLayers.Layer.WMS( "MapServer Layer",                                               "http://localhost/cgi-bin/mapserv?map=/var/www/osm-demo/mapserver-utils-svn/osm-google_test.map&",                                               {'layers': "land0 borders0 places0 land1 borders1 places1 land2 borders2 places2 land3 borders3 places3 land4 landuse4 transport_areas4 waterarea4 borders4 places4 land5 landuse5 transport_areas5 waterarea5 roads5 borders5 places5 land6 landuse6 transport_areas6 waterarea6 waterways6 roads6 borders6 places6 land7 landuse7 transport_areas7 waterarea7 waterways7 roads7 borders7 places7 land8 landuse8 transport_areas8 waterarea8 waterways8 railways8 roads8 borders8 places8 land9 landuse9 transport_areas9 waterarea9 waterways9 railways9 roads9 borders9 places9 land10 landuse10 transport_areas10 waterarea10 waterways10 railways10 roads10 aeroways10 borders10 places10 land11 landuse11 transport_areas11 waterarea11 waterways11 railways11 roads11 aeroways11 borders11 places11 land12 landuse12 transport_areas12 waterarea12 waterways12 railways12 roads12 aeroways12 borders12 places12 land13 landuse13 transport_areas13 waterarea13 waterways13 railways13 roads13 aeroways13 borders13 places13 land14 landuse14 transport_areas14 waterarea14 waterways14 railways14 roads14 aeroways14 borders14 places14 land15 landuse15 transport_areas15 waterarea15 waterways15 railways15 roads15 aeroways15 borders15 places15 land16 landuse16 transport_areas16 waterarea16 waterways16 railways16 roads16 aeroways16 borders16 places16 land17 landuse17 transport_areas17 waterarea17 waterways17 railways17 roads17 aeroways17 borders17 places17 land18 landuse18 transport_areas18 waterarea18 waterways18 railways18 roads18 aeroways18 borders18 places18"}                                                );             map.addLayers([esr_all, mslayer]);         map.addControl(new OpenLayers.Control.LayerSwitcher());         map.zoomToExtent(new OpenLayers.Bounds(69.00,21.00,70.00,22.00));         map.setCenter(new OpenLayers.LonLat(69.57,21.55).transform(         new OpenLayers.Projection('EPSG:4326', 'EPSG:3857'),         map.getProjectionObject()         ), 14);         }     </script>      See the image: ![esr_all layer is showing](http://i.stack.imgur.com/jxGfh.png) I get an error: > msWMSLoadGetMapParams(): WMS server error. Invalid layer(s) given in the > LAYERS parameter. A layer might be disabled for this request. Check > wms/ows_enable_request settings. Instead if change it as a mapserver layer,                   var mslayer = new OpenLayers.Layer.MapServer ( "MapServer Layer",                                               "http://localhost/cgi-bin/mapserv?map=/var/www/osm-demo/mapserver-utils-svn/osm-google_test.map&",                                               {'layers': "land0 borders0 places0 land1 borders1 places1 land2 borders2 places2 land3 borders3 places3 land4 landuse4 transport_areas4 waterarea4 borders4 places4 land5 landuse5 transport_areas5 waterarea5 roads5 borders5 places5 land6 landuse6 transport_areas6 waterarea6 waterways6 roads6 borders6 places6 land7 landuse7 transport_areas7 waterarea7 waterways7 roads7 borders7 places7 land8 landuse8 transport_areas8 waterarea8 waterways8 railways8 roads8 borders8 places8 land9 landuse9 transport_areas9 waterarea9 waterways9 railways9 roads9 borders9 places9 land10 landuse10 transport_areas10 waterarea10 waterways10 railways10 roads10 aeroways10 borders10 places10 land11 landuse11 transport_areas11 waterarea11 waterways11 railways11 roads11 aeroways11 borders11 places11 land12 landuse12 transport_areas12 waterarea12 waterways12 railways12 roads12 aeroways12 borders12 places12 land13 landuse13 transport_areas13 waterarea13 waterways13 railways13 roads13 aeroways13 borders13 places13 land14 landuse14 transport_areas14 waterarea14 waterways14 railways14 roads14 aeroways14 borders14 places14 land15 landuse15 transport_areas15 waterarea15 waterways15 railways15 roads15 aeroways15 borders15 places15 land16 landuse16 transport_areas16 waterarea16 waterways16 railways16 roads16 aeroways16 borders16 places16 land17 landuse17 transport_areas17 waterarea17 waterways17 railways17 roads17 aeroways17 borders17 places17 land18 landuse18 transport_areas18 waterarea18 waterways18 railways18 roads18 aeroways18 borders18 places18"}                                                );      it works and I get tiles for both the layers. Why is it not when both are WMS layers? In my **osm-google_test.map** file the **WEB** portion is like this:               WEB     #IMAGEPATH "/tmp/ms_tmp/"     #IMAGEURL "/ms_tmp/"         IMAGEPATH '/opt/fgs/tmp/ms_tmp/'          IMAGEURL  '/ms_tmp/'     TEMPLATE  '/var/www/mstest/templatesample.html'     #labelcache_map_edge_buffer "10"         METADATA             "wms_enable_request"    "*"             #"ows_enable_request" "*"             "wms_srs" "EPSG:900913 EPSG:4326 EPSG:3857 EPSG:2154 EPSG:310642901 EPSG:4171 EPSG:310024802 EPSG:310915814 EPSG:310486805 EPSG:310702807 EPSG:310700806 EPSG:310547809 EPSG:310706808 EPSG:310642810 EPSG:310642801 EPSG:310642812 EPSG:310032811 EPSG:310642813 EPSG:2986"                 "labelcache_map_edge_buffer" "-10"                 "wms_title" "osm france"         END      END


title: Defining a WHERE clause from a subquery for a query to create a handy view

text: It is about my implementation from an answer of this question. It is more an postgresql issue than a GIS issue, whereas it is about geometries and spatial queries. I want to create a view from the query below, whereby I can select the result by an ID (that defines a linestring). My problem is, that I had no success with coding the query in that way, that for every linestring only one `newpoint` is created that is half the length away from the linestring (see last subquery 'newdistance'). I tried to `SELECT DISTINCT ON (id)...` but that failed. The linestring table:               id smallint NOT NULL,     geom geometry(linestring, 31468)      I have the same problem in another query, to "extract" a WHERE clause (here in the first part of the query) from an subquery. My goal is to have a view `newpoint` (from the query below) whereby I can query like this               SELECT * FROM newpoint WHERE id=1      The query:               WITH line AS         (SELECT             geom         FROM linestring             WHERE id=1),     points AS         (SELECT             ST_StartPoint(geom) AS A,             ST_EndPoint(geom) AS B         FROM line),     azimuth AS         (SELECT             ST_Azimuth(A,B) AS azimuth         FROM points),     distance AS         (SELECT             ST_Distance(A,B) AS length         FROM points),     newdistance AS         (SELECT             length + (length * 1/2) AS newlength         FROM distance)          SELECT         ST_SetSRID(ST_Translate(A, sin(azimuth) * newlength, cos(azimuth) * newlength), 31468) AS newpoint     FROM newdistance, azimuth, points      Setting the WHERE clause id=1 at the end of the query didn't succeed. **EDIT** When I use               SELECT DISTINCT ON (geom) ...      in the last part of the query it works ONLY with 2 objects (linestrings). With more than 2 linestrings it doesn't work. Here is a picture, how does it has to look like. ![enter image description here](http://i.stack.imgur.com/Zk006.png)


title: Admin Sample GXP Viewer

text: I have Geospatial application which is using GXP, GeoServer, GeoExt and PostgisSQL. We would like to enhance the application by adding a user friendly Admin system, similar to what GeoNode has. This sytem would let us create and manage users and create projects specific to users. Can someone give us a head start on how we can do this? Or any place we can download it and configure it according to our requirements. I would really appreciate your help. Thanks


title: How to figure out the SRID of a point by comparing it to another point with known SRID?

text: I have a set of store locations defined in the coordinate system of Bottincarto, which I need to display on Google Maps iOS. Here is an example of coordinates for Notre-Dame-de-Paris:   * 600926 2428323 (Notre-Dame-de-Paris on http://maps.bottincarto.com/...)   * 48.852968 2.349902 (Notre-Dame-de-Paris on http://www.google.fr/maps/...) Here is an example of coordinates for Trafalgar Square:   * 428285 2726837 (Trafalgar Square on http://maps.bottincarto.com/...)   * 51.508039 -0.128069 (Trafalgar Square on http://www.google.fr/maps/...) I already tried a linear regression, but it is not even good as an approximation. How can I reverse engineer BottinCarto projection system?


title: How can I zoomToExtent and then zoom out?

text: I use openlayers 2.12 I use HTML5's geolocation API to get the coords and then I insert them in a layer as feature (point). After that, I zoom in , since the layer has no other feature.               map.zoomToExtent(geoloc.getDataExtent());      `geoloc` is the name of the layer... My problem is, how do I zoom-out automatically after the `zoomToExtent`? Because the map zooms in too much and I can see just the point and very little details. It would be more user-friendly to zoom out a bit Thanks


title: Convert Shape Length Area Data to Latitude Longitude

text: I am attempting to import publicly available shapefile data from http://gis.utah.gov/data/boundaries/zip-codes/ into Tableau for a data visualization project. Unfortunately, the .shp file and .gdb database that are available for download include only the shape length / area for each of the Utah zip codes. My goal is to convert these from the current NAD 1983 UTM 12N (EPSG: 26912) format into latitude / longitude coordinates that can be read into Tableau. I am using Quantum GIS (QGIS 2.0.1) and would appreciate any and all help.


title: Migrating 10.0 map services to ArcGIS 10.2 for Server?

text: I have a few map services that I would like to move to ArcGIS for Server 10.2. Each server is on a separate machine. I have come across instances where people migrate the entire arcgisserver folder, (arcgisinput, arcgisoutput, cache, etc.) however my 10.0 services only maintain an input and an output. What's the most time efficient manner in migrating single map services and not an entire instance from 10.0-10.2?


title: Star Markers - unable to get tie points in Inpho Match-AT

text: What do the star markers mean in Trimble's Inpho Match-AT software, specifically in the multi-photo measurement feature? I have been unable to generate tie points in an aerial triangulation in photos which have no star in the center icon. ![enter image description here](http://i.stack.imgur.com/3XqGJ.jpg) I don't know what they mean or why they don't have stars. The photos are present and connected, they have initialized external orientations and correct camera parameters, and each file does actually have data. I haven't been able to find anything in any help documentation.


title: Qgis installation fails on Ubuntu 14.04

text: I am having problems installing Qgis on Ubuntu 14.04 (Trusty). I tried several approaches: installing from the Ubuntu standard repos, installing from http://qgis.org/debian and installing from http://qgis.org/debian-nightly, using both apt-get and aptitude. Everytime, the installation fails because libnetcdf.so.6 cannot be found:               Setting up qgis-providers (2.3.0+git20140530+ca6d4d9~trusty1) ...     /usr/lib/qgis/crssync: error while loading shared libraries: libnetcdf.so.6: cannot open shared object file: No such file or directory      ls /usr/lib/libnet* gives:               /usr/lib/libnetcdf.a       /usr/lib/libnetcdf_c++.so.4      /usr/lib/libnetcdff.so        /usr/lib/libnetcdf.so        /usr/lib/libnetpbm.so.10     /usr/lib/libnetcdf_c++.a   /usr/lib/libnetcdf_c++.so.4.1.0  /usr/lib/libnetcdff.so.5      /usr/lib/libnetcdf.so.7      /usr/lib/libnetpbm.so.10.0     /usr/lib/libnetcdf_c++.so  /usr/lib/libnetcdff.a            /usr/lib/libnetcdff.so.5.1.0  /usr/lib/libnetcdf.so.7.1.1      So the library is indeed missing. Any help is highly appreciated!


title: Using a shapefile in QGIS makes it unusable in ArcGIS?

text: We currently use Arcgis 9.2, and we are looking at moving forward with other software, and we have been testing QGIS recently. Basically we are trying to copy the existing projects we have, and re-create them in QGIS. But I have came across quite a big problem, re-created a project which has a couple of polygon shapefiles. They came into QGIS fine as always. I then opened up the .mxd project in ArcGIS, and noticed a couple of deadlinks... they were the two shapefiles I had just brought into QGIS! Tried a couple of other random shapefiles and the same thing happens... seems a single shapefile cant be used in both QGIS and ArcGIS at the same time? Looking at the shapefile in catalog and it shows as a green question mark. If i close QGIS the shapefile opens fine again in ArcGIS. I could of course create a copy of each shapefile and point QGIS to the copies, but it means duplicating A LOT of shapefiles, which is far from ideal. Anyone else come across this? And if so, is there a solution/work around?


title: Is it possible to add dynamic controls to a dockable window?

text: I am writing a small application that uses the dockable window. What will happen is that for each layer in the TOC a label and checkboxes will appear within the dockable window. However what do i need to do to add the control onto the Dock Window. As an example i have                               IDataset dataset = enumDataset.Next();                 while (dataset != null)                 {                     MessageBox.Show(dataset.Name);                     //create a label object                      Label label = new Label();                     label.Name = "lbl_" + dataset.BrowseName;                     label.Text = dataset.Name;                     label.Top = top;                     label.Height = height;      But there is no               this.Controls.Add(label);      Does anyone know what the code should be to add the control? Thanks


title: ArcGIS Fill Editing Tool

text: Is there a way to fill while editing in ArcGIS? You know similar to the popular paint tool, where an enclosed area is filled with a colour. I'd like to fill an enclosed area (enclosed by poly lines or polygons or both). I'd like to be able to do this in an edit session, a batch process or tool would not provide the functionality that I am looking for. Here is an image to better explain the situation: ![ArcGIS Fill Example](http://i.stack.imgur.com/HW5AC.png) I need to fill the area between the High Water Mark (the dark blue polyline) and the lake (the blue polygon). I've highlighted the area I need to fill with a red scratch. Now there is a lot of gaps and inaccuracies in the data, so I need to manually check each one as I do it. It takes forever to trace first the watermark and the lake. I don't see why there isn't an ArcGIS editing tool to fill in gaps. If Microsoft can create one anyone can.


title: how to geowebcache my own group layer in geoserver?

text: I am unable to geowebcache my layer group in geoserver. But in Geowebcache demo page, default layer groups are available. How to add my group layer into geowebcache? please guide me.


title: QGIS 1.8.0 ask me to install scipy

text: I need help in a step by step to install scipy via QGIS or via python, I'm not familiar with python language so the guide of Nathan's blog wasn't useful to me due that OSGeo4W don't recognize the commands there described. I have win8 and win7 OS. I'd like to have some help with this issue becuase I can't download LecoS properly that way. TVM!


title: GIS Certification

text: I have a question regarding GIS certification. My main goal is to develop GIS web applications. Currently I am a developer with over 20 years building software and for the past 10 building web application using mostly Java on the server and HTML/Javascript/CSS on the client. I also have a lot of experience with data integration. I am currently enrolled in the Penn State's Postbaccalaureate Graphical Information System program. I just completed the first course this week and there are an 4 additional classes after that. The first course was very good and informative and I feel I learned a lot of background knowledge in the GIS field. I am trying to decide if I should continue on with the additional classes and complete the program. Two reasons why I am reluctant to continue: 1) Not sure if this is the best way for me to get into the GIS industry given I want to develop in the field 2) The classes are expensive Please offer any suggestion on this subject area. Thanks in advance.


title: How to create and visualize groundwater plume models?

text: My water department has asked me about visualization of groundwater plume models. I have no idea where to start. All the links refer me to Cameo/Aloha but these are for Atmospheric modeling as far as I can see. Has anyone any experience with point source plumes in groundwater? I guess the models are similar to those from atmospheric contamination models but I also need to understand which variables are necessary for a successful model.


title: Creating custom size images from WMS cache

text: I am currently working on my own WMS service (using PHP). The service works currently if the requested image is in 256x256px size but currently does not support any other size. My service basically returns cache tiles that were generated by Arcgis. My question is if it is possible to generate correct tiles from these arcgis 256x256px tiles that are larger or smaller then 256x256px. For example if the request is asking for 512x512px tile could I generate the tile on the fly using the 256x256px tiles? One idea on how to do this is by finding all the 256x256px images that are part of the requested larger/smaller image, merge them and cut out the request image.


title: Getting data of an area covered by a polygon in ArcGIS

text: We have a layer that contains land parcel information such as area for each parcel, and a unique parcel ID. We also have multiple polygons that cover some parcels. We're trying to figure out how to determine which parcels are covered by the polygon and which parcels are not covered by the polygon. We really have no idea how to proceed so any help is appreciated. Thanks everybody.


title: How to define which columns in a Spatialite database are the latitude and longitude?

text: I have imported a CSV containing points with coordinates specified in columns named Latitude and Longitude. How do tell spatialite that these columns are the coordinates of of my points?


title: NVIS module not found in QGis

text: When I try to use the nvis module in Grass, QGis, pops-up a warning window saying "Module nviz not found". how can I solve this issue and use the nviz?


title: What lens and camera settings are recomended for UAV (drone) aerial imagery capture?

text: I have been toying with the idea of purchasing/bowrrowing a UAV and attempting to capture my own aerial imagery. I have pretty good DSLR Canon gear and several lenses but I am wondering what lens, focal length and camera settings, including capture intervals, vehicle speed, etc. are recommended.


title: Labeling dilemma

text: How would you label this map logically? I will preempt the question, with the statement it is for a fire truck, so we really have to worry about seeing the address numbers at 3 am with one contact in. I have tried a few different options but any suggestions would be greatly appreciated!![enter image description here](http://i.stack.imgur.com/wyteb.png)


title: How to expedite process of creating features

text: I am creating TIN out of each two cross sections and in the end I merge them all together and create a seamless TIN using python scripting. First I create all the TIN and save them somewhere in my local hard drive and finally I merge them. Someone told me there should be a way to do this process in the memory to save time and them create the final TIN and copy that in my hard drive. I already use feature layer for cross sections to expedite the process. I was wondering if there is still any trick I can use to lower the process time?


title: How to deal with political boundaries in grid cell analysis?

text: I have point data for Africa and I want to aggregate these data to grid cells which I create by a fishnet. Some grid cells cover parts of 2 or 3 countries. Normally, such cells are by different methods completely assigned to a single country. I would like to ask whether there is a method to create grid cells that is more sensitive to national borders? Can I split such cells with country borders and assign them to different countries or is it possible to duplicate such cells and assign them partially to one country? I use ArcGIS 10.2. Thank you for your help!


I use Quantum GIS 2.2.0 Valmiera I have a polyline layer. I can drag&drop polyline nodes by using "Node Tool". Is it possible to drag&drop line precisely 6 meters to right or left?


title: Are the ESRI .lyr and .mxd data format specs published?

text: Are the binary formats of .lyr and .mxd files published?


title: Insert image in attribute table(unable to save edits operation attempted on empty geometry)

text: I have been trying to associate/insert images in an attribute table. I have created a field called pictures and data type is Raster. ( I have used **this resource video** to do the task) Once I load the picture in the attribute table, every thing is ok. I can see the picture in the attribute table. But once I want to save the edits, I am getting this error: unable to save edits operation attempted on empty geometry Thank you in advance for your help. ArcGIS 10.1 (Advanced), Windows 8


title: why readOGR() is not reading projection information from .spj - r

text: 0 down vote favorite I successfully import a polygon shapefile with readOGR, however the information in the .prj file is not being imported.               test class : SpatialPolygonsDataFrame nfeatures : 19407 extent : 35551.4, 1585917, 6318047, 9408727 (xmin, xmax, ymin, ymax) coord. ref. : NA          summary(test) Object of class SpatialPolygonsDataFrame Coordinates: min max x 35551.4 1585917 y 6318047.3 9408727 Is projected: NA proj4string : [NA]      the projection should be utm21n complex zone


title: Snapping doesn't work when many layers loaded!

text: I am using Arcmap 10, when I have only limited number of layers loaded, the snapping experience is perfect but when I have many layers loaded specially polygons, snapping process is painful. My cursor turns to "Busy" shape (the sand clock) and can't snap to the feature. I have no idea why it happens. Any solution?


title: List of actual GPS position

text: I have a question - where can I find a list of actual GPS position, for example - ships, airplanes etc. I need it for my project for College - I will create mobile/desktop/web app, which will show on map current positions of this elements, their last road and save info about that elements in database. Also, is there any website with informations about weather problems with lat/lng info? Thanks in advance, I've tried to find something about it, but without result..


title: ArcGIS Server REST get polygon centroid

text: I have a Site reference feature service data in polygons served up by ArcGIS Server. I need to get the centroid coordinate of a polygon by the Id but can’t find a method out of the box. I’m working on a GeoEvent Processor Service that will use this Site service to grab a centroid from a site Id. So far I believe my options are:   * Custom GP service, python – takes in polygon rings, returns point   * Server Object Extension – same as above   * Sql view, query layer served up as a feature service that uses STCentroid (The database uses a SQL GEOGRAPHY type) on the Site layer. Would query the layer by ID to return the geometry.  The first two options keep me on track with a service orientated architecture however it certainly seems overkill to be writing an SOE for one function. The final option could be a ‘goer’ but it means my logic is coupled to SQL Server.


title: repeat skobbler maps horizontally with OpenLayers?

text: I'm trying to get skobbler maps to repeat horizontally if you load this page http://developer.skobbler.com/demo/skobbler-tiles-openlayers.html and zoom out, the map will not repeat horizontally However if you use this link and zoom out http://developer.skobbler.com/web- tiles/embedded it works as I want (i.e. it repeats across the screen) However, I'm using OpenLayers and not an Iframe Is what I want to do possible with OpenLayers? Thanks Mark


title: Applying same color ramp over multiple maps with different values in ArcGIS10?

text: I'm doing a rainfall analysis for multiple events in ArcGIS10.1, and I want to apply the same color ramp (with a constant range) for each map. For example, all of the events I have range from zero to roughly 6.5 inches. The color ramp I created is a simple rainbow scheme, with red representing 0 and purpleish representing 6.5. My problem is is that obviously, not all of the events range from 0-6.5. Some go from 0.3-0.7, others from 0.1-1.3, etc. How do I go about applying the 0-6.5 color ramp on each map so that the values are represented with the proper color? Let me know if this question needs any clarification. Thanks! EDIT -- So perhaps I'm doing something wrong. ![](http://i.imgur.com/TfEydve.jpg) I made the color ramp as shown (saved it as a layer file and transformed it into a raster) and when I go to import a symbology to one of my interpolated layers ![](http://i.imgur.com/CuTq03g.jpg) as shown, it stretches the entire color ramp from the lowest value to the highest value, not just the parts of the color ramp it "should". Any advice on what I'm missing?


title: Is it real to install PostgreSQL 9 on the virtual machine Windows Server R2?

text: I have Virtual Box with GuestAddition. And virtual machine "Windows Server 2008 R2" When I use real computer with Windows Server 2008 R2 I can install postgresql-9.2.3-2-windows. But when I install postgresql-9.2.3-2-windows on virtual machine I have error "Problem running post-install step. Installation may not complete correctly" Please, help me.


title: Why are the results different for the same point using GEOS Point.transform and Postgis ST_Transform?

text: I'm trying to define a default point setting for OSM layers in my application and I want users to be able to define the default latitude and longitude in _EPSG:4326_. What I'm doing right now is to use GEOS `Point.transform` to turn the EPSG:4326 input into _EPSG:3857_ :               _default_point = Point(         getattr(settings, 'DEFAULT_LON', None) or 0,         getattr(settings, 'DEFAULT_LAT', None) or 0,         srid=4326     )     _default_point.transform(3857)     default_lat, default_lon = _default_point.get_coords()      Then, i.e. for `(-2.41544723477347,43.2140883725249)` I'm getting `(-268966.0382412352, 4810570.312724996)`. That, as for what I can see in the map, is wrong. Now, when trying over psql               SELECT ST_AsText(ST_Transform(ST_SetSRID(ST_Point(-2.41544723477347,43.2140883725249),4326),3857));                      st_astext                      -------------------------------------------      POINT(-268886.356213004 5344615.26048709)      Which is right. What is going on then? Why are the results different? EDIT: `select srid, proj4text from spatial_ref_sys where srid = 3857 or srid = 900913 or srid = 4326` yields:                  4326 | +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs         3857 | +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +units=m +k=1.0 +nadgrids=@null +no_defs      900913 | +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs


title: Open Arabic Shapefile in Qgis

text: I have a shape-file all it's fields are in Arabic. When I opens it in qgis 2.0 all the data is replaced by "?". I changed the encoding to Arabic or UTF-8 but with the same result. How can I Open an Arabic Shapefile in Qgis.


title: How to use GDAL with PyPy?

text: I am trying to optimize my python code using pypy. While I am able to import most libraries and run succesfully, I am not able to import GDAL libraries. The same libraries when imported in Python 2.7 run fine. I have tried adding GDAL location to PYTHONPATH, and even tried downloading GDAL Python bindings to PyPy directory, somehow pypy does not want to recognize GDAL, and comes with the following error.               Traceback (most recent call last):            File "app_main.py", line 51, in run_toplevel            File "ImportENC-P.py", line 50, in <module>              from osgeo import gdal, ogr, osr            File "B:\PYPY\site-packages\osgeo\__init__.py", line 21, in <module>              _gdal = swig_import_helper()            File "B:\PYPY\site-packages\osgeo\__init__.py", line 13, in swig_import_helper              import _gdal          ImportError: No module named _gdal


title: Geoserver WFS GetFeature returns no features when reprojecting PostGIS layer

text: I may be misunderstanding the GeoServer documents, but when accessing a GeoServer 2.2 WFS 1.1 layer, should I be able to reproject to an arbitrary SRS via the SRSNAME parameter? I'm generating a bounding box from the location of a click on a WMS map where I know there is a feature. The map is projected into UTM zone 15 (EPSG:32615). So I generate the following URL: http://server.domain.com/geoserver/MyWorkspace/wms?SERVICE=WFS&VERSION=1.1.0&REQUEST=GetFeature&TYPENAME=MyWorkspace%3AMyLayer&MAXFEATURES=9&SRSNAME=EPSG%3a32615&BBOX=382138.16774714%2c4587006.11489167%2c382138.16774714%2c4587006.11489167 This query returns 0 features, and I don't know why. It's essentially the same query I was using for another layer, also served from GeoServer, that was sourced from a shapefile in EPSG:32615, where the query worked. So here's what's different about this layer:   1. I'm using PostGIS for a backend on the new layer; the old one was a shapefile   2. The shapefile was in EPSG:32615; the PostGIS table is in EPSG:4326   3. The PostGIS table is using the "geography" data type rather than the "geometry" data type, if that's relvant. Any ideas?


I can't seem to access the plugin repository. Is this an issue with my system or is there something going on at the server end?! I get the following error as wells sometimes a window saying: "Error reading repository: QGIS Official Repository HTTP request failed" It's snowing like crazy here and that often messes with the internet but every thing else seems to be working fine.


title: openlayers load local gml

text: I'm trying to load a simple gml file with openlayers, but it doesn't load. Pretty sure it's a simple issue, but I figure out what's wrong. Thanks in advance, Pedro                           var meu = new OpenLayers.Layer.Vector("States", {                 protocol: new OpenLayers.Protocol.HTTP({                     url: "teste.GML",                     format: new OpenLayers.Format.GML()                 }),                 strategies: [new OpenLayers.Strategy.Fixed()],             });                  map.addLayer(meu)      Tried firefox and safari in mac os and firefox in linux


title: How to programmatically add a new plugin repository in QGIS?

text: How can I add the name of a repository and URL to the third-party plugin repository using python in QGIS 2.01?


title: For loop within SDE

text: I have a list of SDE connection files in a folder in my computer. I can use arcpy.listWorkspaces("*", "SDE") to find the workspaces. But I then want to dig deeper into each of the SDE files and list the feature classes. I tried:               sdeFiles=[]     listWorkSpaces=arcpy.ListWorkspaces("*", "SDE")     for sde in listWorkSpaces:         arcpy.env.workspace=sde         getFeatureClasses=arcpy.ListFeatureClasses()         getFeatureClasses.sort()         sdeFiles.extend(getFeatureClasses)     print sdeFiles      So when I print the sdeFiles list, it will list all the SDE feature classes in all the SDE files, but sometimes it will list them 6+ times, so it's duplicating it. I imported the list to excel and removed duplicates. It created 1514 duplicates and only 650 unique values. Frustrating! Can anyone tell me why?


title: How to display information in a div outside the map?

text: I tried a bit around and created a simple map with two polygons(most of the code I got from the Interactive Choropleth Map Tutorial). I added a `"L.control()"` to show information about the polygons in a `div`. Contrary to the tutorial I would like to have this `div` outside the map. To get a clue what I mean please look here: http://jsfiddle.net/Tbr9p/


title: postgis and osm2pgsql: authentication failed

text: Id like to import osm into my postgis to use it with qgis. I followed exactly these instructions: http://wiki.openstreetmap.org/wiki/PostGIS/Installation At step #Activate_PostGIS:               psql --username=gisuser --dbname=gis --command="\d"      This error occured:               psql: FATAL:  Peer authentication failed for user "gisuser"      (Same after importing my *.osm by doing this:               osm2pgsql -s -U gisuser -d gis path-to-file.osm      Error:               Error: Connection to database failed: FATAL:  Peer authentication failed for user "gisuser"      I have no idea what to do, so im asking here. Using xubuntu 12.10, postgres 9.1.11


title: How to use MOLUSCE plugin in QGIS?

text: I am trying to explore MOLUSCE ver 3.0.4 in QGIS Dufour. I am trying to simulate the land use change using this tool. I have following querries:   1. It asks for Initial and final year at the first step. Where to enter the exact year for which user is interested for prediction?   2. It gives output file of simulation. For which year this simulation is? Since we already have given Initial and Final year.


title: Open source/free software for generating a 3D view of mapping?

text: Is there any software that will generate a 3D view given a DTM/height data and raster mapping, and output an image in a printable resolution? I am looking for something similar to what you'd find in a commercial hike planning application, but free or open source preferably.


title: How to utilise single 'randomForest' model for several rasters in R?

text: I am using randomForest package in R to derive species distribution models. When I predict the model to the RasterStack containing the predictors in which the values were extracted to build my model, everything goes fine. The problem is that I would like to predict this model I generated in the present scenario to the same set of predictors in the past and future. However, R gives me an error message. Past and present predictors RasterStack are the same extension, resolution, etc. This is not the problem. Reading the help of the 'predict' function {raster package} in R, I found the following: "The names in the Raster object should exactly match those expected by the model. This will be the case if the same Raster object was used (via extract) to obtain the values to fit the model". So, according to this, I can't predict my model to another RasterStack that was not used to obtain the values for the model. So, I would like to know another way that I could predict values from a present scenario to past or future scenarios. Thanks a lot, Mariana


title: How to get hypsometric curve from elevation data?

text: I have my elevation data but cannot make an attribute table from it and so I'm not sure how to get any of that data into a table in order to make a hypsometric curve.


title: Missing level for cache map service

text: I am having this issue where caching is missing one of the levels (level 13). I am trying to find out how to resolve this issue. L9 - 1,155,581 - OK L10 - 577,790 - OK L11 - 288895 - OK L12 - 144447, OK L13 - 72223, MISSING, don’t show at all for some areas L14 - 36111, OK L15 - 18055, OK L16 - 9027, OK L17 - OK L18 - OK Thanks Jay


title: Disappearing Point Features in Geoserver using WMS

text: I have a shapefile with roughly 6,500 points worldwide in it that I am trying to serve with Geoserver 2.2.1 using WMS. All has been apparently fine with it until I implemented a filtering feature on my client application which uses leaflet. When I append a CQL_FILTER (attribute filter, not spatial) into the WMS query, I noticed missing features when zoomed out. When I zoomed back in, they would sometimes reappear, but not always. See the image below - ![Side by side comparison](http://i.stack.imgur.com/IoA2B.png) At the zoom level on the left, Atlanta is not displayed. When I zoom in, it is. However, sometimes not even the point in Tampa is displayed at the zoom level on the left. If I zoom out 3 more levels, no points at all will be displayed. I'm not sure the issue is the CQL_FILTER param, as with 6,500 points it's hard to notice a few missing points at a global scale, but the particular filter I'm showing here as an example filters down to only 3 features, and when 1 to 3 of them are missing depending on the zoom level it's particularly noticeable. I can recreate this same behavior in geoserver's own layer preview using the CQL Filter in the advanced map options, so I'm fairly certain it is not my client application. As for my own steps at debugging this issue, I have tried the default Points style/sld file to rule out my own layer styling. I have turned off all the caching that I am aware of. I double checked that my projections are correct - I created the shapefile in ArcGIS 10 using WGS_1984_Web_Mercator_Auxiliary_Sphere as the projection, and the layer is set to EPSG:3857 in geoserver, which I think is equivalent. I also upgraded from geoserver 2.2 to 2.2.1 and have had the same issue in both. I have also deleted geoserver's spatial index file (.qix) and let it get recreated since I've seen similar issues in Arc with corrupt spatial indexes, but obviously that did not work either. Here's a snapshot from Geoserver's own layer preview with the CQL Filter on and zoomed into the same area as shown above. The red circle is approximately where I should be seeing another point (Atlanta). ![Openlayers Example](http://i.stack.imgur.com/Zutqm.png) I've tried tweaking every other setting I can think of, but am having no luck. I've also looked up and down the geoserver logs and enabled verbose logging, and do not see any errors/exceptions. I also do not see any bad requests in chrome's developer tools. Thanks in advance for any help. If I'm missing any critical info, I'll provide what I can, but this is an internal/non-public application.


title: Multiple queries in a row

text: I have osm data imported into PostgreSQL with osm2po tool (wich makes the data routable). Also I have imported data relationed with nodes and node_tags with osmosis tool. I need to obtain a set of nearest nodes from a given location (maybe 30 or 40). Then I need to calculate the shortest path from the given node from each others. It means to execute shortes_path query multiple times in a row, wich takes about 2 seconds (1 query 2 seconds: 30 queries 60 -> too long).               SELECT * FROM shortest_path('        SELECT id AS id,               source::int4 AS source,               target::int4 AS target,               cost::float8 AS cost        FROM tt_2po_4pgr',     275742,     274494,     false,     false)      What would be the best way to do this in a reasonable time? Thank in advance! Pd: excuse my english **Edited:** In accordance to _dkastl_ answer, I tried with wrapper functions and effectively I obtain better results 2 or 3 times better, but isn't acceptable for me yet. With the following query:               explain analyze SELECT *         FROM dijkstra_sp_delta('dr_2po_4pgr', 278174, 289786, 0.3);      I obtain this results: explain


title: Left click then right click is treated as double click in Mac QGIS Layers

text: I am having an issue with QGIS 1.8.0 on my Mac where I left click then right click which gets treated as a double click if the timing is too close together and brings up the Layer Properties box. This usually happens when I am focused on a different program, left click to gain QGIS' focus and then right click to get a context menu for a layer. This is force of habit I guess.   1. Is there a way to disable the Layer Properties from appearing with a double click? I know that you can choose between the Attribute Table and the Properties Box.   2. Should this issue be added as a bug report to QGIS?


title: Enabling ArcSDE Archiving with Python and C API using ArcGIS 10.0?

text: I have a script that does some register as versioned and updating privileges. How would I call into the c api for ArcSDE to enable archiving with python? http://help.arcgis.com/en/geodatabase/10.0/sdk/arcsde/api/capi/versioning/archcreate.htm Any direction would be helpful.


title: why args.ScreenPoint.IsEmpty is true?

text: I use a map control and set a MapClick event for it. And i have a tool that ClientAction set to point. Event fired when i clicked on the map. But args.ScreenPoint was empty. What's wrong ? I want get x and y in map click.


title: How to create a line with start and end points which are in another layer in QGIS?

text: In my map, there are lots of points in layer **A**. I want to make a line layer **B** and I want to make start & end points of lines snapping to points in layer A manually. Can I use points of layer A as a coordinate attribute selection index for my line layer B? Thank you.


title: Very slow connections to ArcSDE with over 20,000 datasets

text: I have a problem with ArcSDE that I have never come accross before and I was wondering if anyone could shed any light on this problem. **Problem ArcSDE 9.3.1:** > Connnecting to an ArcSDE instance (Oracle Spatial) with over 20,000 datasets > (raster) causes ArcGIS/ArcObjects code to connect very slowly - minutes in > fact before a successful connection is made to the database. However if a > native SQL connection is made, the connection is almost instant. Has anyone else had this problem and if so what can be done to remedy this problem.


title: Command doesn't show up on ICommandBar

text: I'm trying to put a custom command into the Help menu on startup if it isn't already there. This is my code to do so:               private void OnCreate()     {         // other stuff in OnCreate()...         // registering my method for the right event         m_appStatusEvents = m_application as IApplicationStatusEvents_Event;         m_appStatusEvents.Initialized += new IApplicationStatusEvents_InitializedEventHandler(AddCmd);     }          private void AddCmd()     {         UID myCmdId = new UIDClass();         myCmdId.Value = myCmdIdValue;  // myCmdIdValue is my command's UID              ICommandBars commandBars = m_application.Document.CommandBars;              UID helpMenuId = new UIDClass();         helpMenuId.Value = helpMenuIdValue; // helpMenuIdValue is esriArcMapUI.MxHelpMenu         ICommandBar helpMenu = null;         helpMenu = (ICommandBar)commandBars.Find(helpMenuId, false, false);         // here, helpMenu.Count is 3         if (helpMenu != null)         {             ICommandItem myCmd = helpMenu.Find(myCmdId, false);             if (myCmd == null)             {                 Object index = Type.Missing;                 helpMenu.Add(myCmdId, ref index);                 // here helpMenu.Count is 4, so my command was added             }         }     }      I think it's pretty straight forward: I'm looking for the Help menu and add my command if it isn't already there. In the debugger, I can see that an item is added to the Help menu, as its `Count` goes from 3 to 4. Problem is: The command doesn't show up. I can add it manually, so I know it's loaded and correct, it just doesn't show up when I try to add it automatically. Edit: I just tried to refresh the Help menu with `(helpMenu as ICommandItem).Refresh()` but that didn't help either. Edit2: I just added this to debug:               for (int i = 0; i < helpMenu.Count; i++)     {         Debug.WriteLine(helpMenu[i].Caption);     }      And it shows   * &ArcGIS Desktop-Help   * &ArcGIS Resource Center   * &About ArcMap...   * myCommand Caption So it's definitely there, it just won't show!