title: distanceTo is giving wrong result in OpenLayers?

text: Hi i am using `distanceTo from OpenLayers` to calculate the distance between two locations these are map details                   Geographic  = new OpenLayers.Projection("EPSG:4326");          Mercator = new OpenLayers.Projection("EPSG:900913");                  map = new OpenLayers.Map({        div: "map_canvas",        projection: Mercator,        displayProjection: Geographic,        center: new OpenLayers.LonLat(0, 0),        minResolution: "auto",        minExtent: new OpenLayers.Bounds(-1, -1, 1, 1),        maxResolution: "auto",        maxExtent: new OpenLayers.Bounds(-180, -90, 180, 90),            units : 'km'      }); these are my lat and lon values               lat = 47.162494, lon = 19.503304     lat = 46.2525684, lon = 20.147060499999952      google maps function result `111861.5962110549`               function distanceBetweenMarkers(latlng1, latlng2){         return google.maps.geometry.spherical.computeDistanceBetween(latlng1, latlng2);     }      openlayers function result `6237980.800397318`               function distanceBetweenMarkers(latlng1, latlng2){             var point1 = new OpenLayers.Geometry.Point(latlng1.lon, latlng1.lat);             var point2 = new OpenLayers.Geometry.Point(latlng2.lon, latlng2.lat);             return point1.distanceTo(point2);         }      why is the difference between `google maps computeDistanceBetween and openlayers distanceTo?`


title: Determine direction of rivers, which villages lie upstream and which downstream

text: I want to determine which villages lie upstream and which villages lie downstream big cities built on the side of rivers in Indonesia. I have the .shp file of villages and cities, the .shp file of all rivers in the country and some .tif files showing elevation (standard SRTM rasters, although I cannot find the exact resolution right now). I know how to identify villages located close to rivers (using spatial join with the rivers .shp file). However, I do not know how to determine the direction of the rivers. Should I discard my rivers .shp file and try to create a Digital Elevation Model from the SRTM .tif data? Using the .tif data I just did: FLOW DIRECTION + FLOW ACCUMULATION, but the computer is stuck since 40 minutes ago and even if it works I would not know how to continue.. Any suggestion? Thank you. /michele


title: QGIS WMS ver 1.3.0 Exception with GeoServer

text: I'm trying to load WMS layers from a local GeoServer into QGIS so that I can edit them. I can access and edit the same layers if I access them directly from a PostGIS database. However, trying to load them as WMS layers in QGIS fails with an XML Parsing Exception:                <<?xml version="1.0" encoding="UTF-8"?><ServiceExceptionReport version="1.3.0"    xmlns="http://www.opengis.net/ogc" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"   xsi:schemaLocation="http://www.opengis.net/ogc http://localhost:8080/geoserver/schemas /wms/1.3.0/exceptions_1_3_0.xsd">   <ServiceException>       javax.xml.transform.TransformerException: Translator error     Translator error     Error occurred trying to write out metadata for layer: Pk50095     Unable to acquire a reader for this coverage with format: WorldImage     </ServiceException></ServiceExceptionReport>     Tried URL: http://localhost:8080/geoserver/wms?SERVICE=WMS&REQUEST=GetCapabilities      In addition, when I try to access the `GetCapabilites` URL directly, an XML Parsing error is generated for WMS ver 1.3.0. However, version 1.1.1 GetCapabilities URL executes fine. Any idea on what may be wrong with my server installation?


title: How can I use R-Tree to find points within a distance in Spatialite?

text: Create a table with 3D points:               create table points (coords pointz);     select CreateSpatialIndex('points', 'coords');      Insert 13k rows:               insert into points values (geometryfromtext('pointz(...)'))      Do               SELECT         p.rowid, count(*)     FROM         points p, points v     WHERE         distance(v.coords, p.coords) < 1     GROUP BY         p.rowid;      Will see an output of several (1..10) points per second (will take a few hours to process 13K points). `EXPLAIN` shows that the index is not used. A query with `PTDistanceWithin` function works the same way. How can I use the spatial index? I checked if it's possible to use an MBR, but those functions are for 2D only. **added:** here's a function reference, and all the functions I could use in usual case, work in 2D only!


title: Modifying the LocateWidget.mxml Locate and Clear functions

text: I am using Robert's compiled TOC widget and attempting to modify the LocateWidget.mxml for the flexviewer. The workflow I have in mind is that user performs an address search clicks the locate button and all my defined operational layers turn on. The clear button would then turn them off. My attempt at the code is posted here: http://forums.arcgis.com/threads/57851-Modifying-the-LocateWidget.mxml-Locate- and-Clear-functions


title: How to get spatialite_osm_net from spatialite

text: I am using windows 7. I downloaded spatialite-tools (precompiled) from the following site. But there is no "spatialite_osm_net" tool available. For one of my project I was trying to convert the OSM road network to spatialite road network. http://www.gaia-gis.it/spatialite-2.3.1/binaries.html It looks like without this tool It is not possible to do it. For your information, I tried to make this tools from the following site according to instruction using MinGW. But It does not work when I come to $make command. I do not know is there anything missing? or how Should I do?               $ cd c:/spatialite-tool-2.4.0          $ export "CFLAGS=-I/usr/local/include"          $ export "LDFLAGS=-L/usr/local/lib"          $ export "PKG_CONFIG_PATH=/usr/local/lib/pkgconfig"          $ ./configure --target=mingw32          $ make          $ make install-strip      So now how can I get spatialite_osm_net???


title: Disable layer checkbox if no geometry on that zoomlevel (source=wms)?

text: I have a map application using GXP which displays different data coming from different WMS sources. Some geometry is only available on certain zoomlevels (available via MaxScaleDenominator on the GetCapabilities). It should be only possible to select an overlay layer when you are at that minimum zoom level, otherwise it should be disabled. Any idea on how to do this?


title: Trying to make xml data readable for Arcgis

text: I am trying to convert xml data into tabular form for the process of geocoding. We have the xml file with info that we are going to run through an address locator. Do I need a script for this? Can we import directly into arc without coverting to xls, maybe using Python?


title: How to automatically reduce polygon size?

text: I am working with two polygon layers (counties & states) which I will 'intersect' to find which counties are within which states. Due to data inaccuracies I have to change all county polygons and make them a bit smaller. At this stage it doesn't really matter by how much, as long as I can do this change automatically. I understand that polygons which are currently perfectly aligned now, will afterwards have a 'gap' between them, but that's fine. I've tried 'Simpliyfy Geometries' but this only 'smoothens' edges and doesn't help in this case. Is there any option in QGIS to reduce the size of polygons automatically? (i.e. move all polygon points 200 meters closer to the 'polygon mid-point'). Cheers !


title: GeometryServer > Project with Datum Transformation

text: I am facing the problem that when i use the extent in the url and i go from Amersfoort(WKID 4289) to WGS84(WKID 4326) that without the tranformation the point from longitude and latitude is placed about 100m away from the original point. In our company this isn't what we want. In ArcGIS 10.1 it is possible in the rest interface to apply the transformation, but i get a repeated error, when i remove the transformation i get no error. I use the Transformation: esriSRGeoTransformation_Amersfoort_To_WGS1984 Can you tell me what i am doing wrong, i put a image with this post and a sample url: https://ourserver.com/arcgis/rest/se...d=false&f=html ![REST Interface GeometryServer](http://i.stack.imgur.com/cDXFY.png)


title: Why am I unable to publish a postgis layer to geoserver?

text: I left arcgis and microsoft sql and now I am trying to connect Qgis+PostgreSQL (PostGIS) + GeoServer. I've brought my postgis database to the server, created a "store". And now I'm trying to publish a layer. So I went to Layers- Choose my layer- Publish- Edit Layer (Data/Publishing)- configure everything there- SAVE and nothing happens... I can't see my layer in "Layer Preview" So maybe configuring was wrong. I chose the native SRS, point style, authority name: admin, authorithy URL (is it the server's URL?? ). What is layer identifier, as well? Thank your patience and answers:)


title: GeoJSON is valid according to JSONLint but jquery gives error

text: I have the following GeoJSON (from webservice)               {     "Features": [         {             "Id": "1",             "Geometry": {                 "Coordinates": [                     {                         "Latitude": 31.2812,                         "Longitude": 34.8027,                         "Altitude": 0                     }                 ],                 "Type": 0,                 "CRS": null,                 "BoundingBoxes": null             },             "Properties": {                 "Color": "#00FF55"             },             "Type": 7,             "CRS": null,             "BoundingBoxes": null         },         {             "Id": "2",             "Geometry": {                 "Coordinates": [                     {                         "Latitude": 31.7825,                         "Longitude": 35.1966,                         "Altitude": 0                     }                 ],                 "Type": 0,                 "CRS": null,                 "BoundingBoxes": null             },             "Properties": {                 "Color": "#FF5500"             },             "Type": 7,             "CRS": null,             "BoundingBoxes": null         }     ],     "Type": 8,     "CRS": null,     "BoundingBoxes": null      } This is the webservice code :               public string GetItemsFromDBGeoJSON()     {         List<GeoJSON.Net.Feature.Feature> fList = new List<GeoJSON.Net.Feature.Feature>();         using (SqlConnection conn = new SqlConnection(connectionString))         {             using (SqlCommand cmd = new SqlCommand("dbo.GetGPSData", conn))             {                 conn.Open();                 SqlDataReader reader = cmd.ExecuteReader();                      while (reader.Read())                 {                     string name = (string)reader["GPSUserName"];                     double lon = Convert.ToDouble(((decimal)reader["Lon"]));                     double lat = Convert.ToDouble(((decimal)reader["Lat"]));                     string clr = (string)reader["GPSUserColor"];                          GeoJSON.Net.Geometry.GeographicPosition latlon = new GeoJSON.Net.Geometry.GeographicPosition(lat,lon,0);                     GeoJSON.Net.Geometry.Point point = new GeoJSON.Net.Geometry.Point(latlon);                                              Dictionary<string,object> properties = new Dictionary<string,object>();                     properties.Add("Color", clr);                     GeoJSON.Net.Feature.Feature feature = new GeoJSON.Net.Feature.Feature(point, properties);                     feature.Id = name;                                              fList.Add(feature);                 }                 conn.Close();             }         }         GeoJSON.Net.Feature.FeatureCollection fCol = new GeoJSON.Net.Feature.FeatureCollection(fList);         JavaScriptSerializer serializer = new JavaScriptSerializer();              return serializer.Serialize(fCol);     }      The Ajax call gives me error :               SyntaxError: JSON.parse: unexpected character      Here is the request using JQuery               $.ajax({         type: "POST",         url: "http://localhost:36361/Service.asmx/GetItemsFromDBGeoJSON",         dataType: 'json',         success: function (response) {             console.log('success');         },         error: function (req, status, error) {             console.log(error);         }     });      Also , Using Firebug Console , I see this response from the request:               <?xml version="1.0" encoding="utf-8"?><string xmlns="http://tempuri.org/">{"Features":[{"Id":"1","Geometry":{"Coordinates":[{"Latitude":31.2812,"Longitude":34.8027,"Altitude":0}],"Type":0,"CRS":null,"BoundingBoxes":null},"Properties":{"Color":"#00FF55"},"Type":7,"CRS":null,"BoundingBoxes":null},{"Id":"2","Geometry":{"Coordinates":[{"Latitude":31.7825,"Longitude":35.1966,"Altitude":0}],"Type":0,"CRS":null,"BoundingBoxes":null},"Properties":{"Color":"#FF5500"},"Type":7,"CRS":null,"BoundingBoxes":null}],"Type":8,"CRS":null,"BoundingBoxes":null}</string>      Any idea why ? I want to show it on leaflet layer.


title: How can I map a matrix of origins and destinations in QGIS?

text: I have a central point on a map, and I have other points in a separate layer (the location of a venue, and the origin cities of the venue visitors). Is there a plugin which would generate output similar to airline destination maps like this; ![a example](http://www.r-bloggers.com/wp-content/uploads/2010/11/map_web.jpg) They don't have to be "great circle" type of lines, just straight lines will do. This is for visual representation of point of origin for visitors of certain events. Any tips or plugin recommendations?


title: How display tiles by using earth plugin in browser?

text: I was using MAp Tiler to generate tiles for my maps, And everything is great. But I have problems with displaying those tiles in 3D version by using earth plugin. The data which gave me Map tiler was with files but for API 2, but I want to use now API 3 version. All works fine with using maps, bu I cannot dislay tiles by using earth plugin. Please give me soem egzamples, if anyone had to do something with this subject. Regards


title: Get coordinates in both decimal degrees and meters in an attribute table

text: I have an attribute table of points and the coordinates are stated in decimal degrees. I would like to add two fields where the coordinates are also stated in meters or kilometers. I am using ArcGIS 10.2.


title: Looking for guidance with beginner use cases

text: **Background** I'm more or less completely new to GIS. At work we use a program called Microvision which does stuff like work out drive distance (in addition to straight line distance) from which you can then create distance bands. Plus has some other mapping features like creating layered ranges. However it has got to be one of the worst programs I have ever had the displeasure of using (not in terms of features, but just in the way it has been coded). So I'm looking at open source alternatives instead. My problem is that I really have no idea where to start. I downloaded SAGA but quickly realised how far out of my depth I am. The features SAGA boasts are cool, but in some respects it is may more powerful than I need. I also don't understand how to import a base map into SAGA. I found some demo maps which are cool, but I have no idea how to create those maps in the first place. Someone on another forum suggested Quantum GIS, so I've going to give that a go very soon. **My main needs/desires are below,** if you can suggest something to match them, I'd be really grateful:   1. A prebuilt map of the world (or least Australia), or some process to follow to create one. This is so I can then overlay our own maps and views over top of it.   2. Something with a lot of customisation and preferably a good API. An ability to manipulate and transform data within the program would be great (even if it requires grunt work coding to do). Microvision is one of the most rigid pieces of software I've ever used and doesn't let you customise anything.   3. Typically my data will be at a postcode level, so the ability to map postcodes to latitudes and longitudes would be awesome. Something that comes with postcode catchments, or even census collection districts for Australia would be AMAZING.   4. Drive time, aka googlemaps style. This isn't essential, but it would really help us for some of the work we do. Of course I can do approximations using straight line distance, but Australia isn't the flattest of countries so something that can look at drive distances would be fantastic. Help appreciated :)


title: Batch model stores incorrect model parameter values

text: I am trying to batch process about ten files using a model I created. There are five model parameters, which I can fill in for the first feature class without a problem. However, when I try to fill in the output file information for the subsequent feature classes, the model "remembers" the values for the first feature class and changes it. For example, for the output for the intersect of `fc1`, I specify `fc1_intersect`. Likewise with `fc2`, I specify `fc2_intersect`. As soon as I click away to fill in the next values, the output of `fc2`'s intersect will be changed to the same as `fc1`'s intersect. This causes the model to fail. I have not specified any default output locations or anything. I used to have a similar problem in ArcMap, where the tools would "remember" the previous layers used, but closing and reopening it would solve it. This is not working now with ArcCatalog, and the time I was hoping to save by batch processing is now being wasted trying to figure out why it keeps changing my input.


title: Get a pixel value

text: How can i do this operation: I want to **get a pixel value** fro the satellite image `(Landat, GeoTIFF, 7 bands)` using qgis-python. What operators should i use? (for example in qgis's Python Console)


title: How to vectorize a Raster-Classification? (qgis, grass)

text: I classified my raster-DHM into elevation levels using colormap. Its not a "real" classification, just a colormap. So I cant vectorize the classes to get polygones for every class. Whats the best way, to classify and vectorize a raster? The Rastercalculator? How to use it, if i want to classify my DHM in elevation gradiation like: 1000-1200m; 1200-1400m? Thanks! Martin


title: How do I start to debug COM error "HRESULT E_FAIL" for TableToGeodatabase geoprocessing tool?

text: I've got some code that creates a table using OLEDB connections to an access table. I want to convert this table to a Geodatabase table, so I ask a question How do I convert an access table to a proper geodatabase table? I tried to use the code there, which threw HRESULT E_FAIL, and then after reading Blah's link about using the geoprocessor I came up with this code:               Geoprocessor GP = new Geoprocessor();          ESRI.ArcGIS.ConversionTools.TableToGeodatabase tableToGeodatabaseTool =         new ESRI.ArcGIS.ConversionTools.TableToGeodatabase();               tableToGeodatabaseTool.Input_Table = _cSM.Geodatabase_FilePath + "\\" + TableName;     tableToGeodatabaseTool.Output_Geodatabase = _cSM.Geodatabase_FilePath;          GP.Execute(tableToGeodatabaseTool, null);      ... Which also throws HRESULT E_FAIL. I'm not sure why, and I'm not sure how to go about debugging it. This code is being run from a custom toolbar in ArcGIS. It seems like ArcGIS is locking the access database, but I assume because I'm using ArcObjects that's not the problem. Any ideas? **Edit** : Here are the details to the exception that is being thrown. You will note that "inner exception" is empty; it is returned as null:               System.Runtime.InteropServices.COMException was unhandled by user code            Message="Error HRESULT E_FAIL has been returned from a call to a COM component."            Source="ESRI.ArcGIS.Geoprocessing"            ErrorCode=-2147467259            StackTrace:             at ESRI.ArcGIS.Geoprocessing.GeoProcessorClass.Execute(String Name, IVariantArray ipValues, ITrackCancel pTrackCancel)             at ESRI.ArcGIS.Geoprocessor.Geoprocessor.ExecuteInner(IGPProcess process, ITrackCancel trackCancel, IGeoProcessor igp, IVariantArray iva)             at ESRI.ArcGIS.Geoprocessor.Geoprocessor.Execute(IGPProcess process, ITrackCancel trackCancel)             at Reports.Validation.validate_StumpageValues() in C:\Users\Public\Documents\SilvAssist\Reports\Reports\Validation.cs:line 131             at Reports.frm_Options.run_Validations() in C:\Users\Public\Documents\SilvAssist\Reports\Reports\frm_Options.cs:line 160             at Reports.frm_Options.bgWorker_Validations_DoWork(Object sender, DoWorkEventArgs e) in C:\Users\Public\Documents\SilvAssist\Reports\Reports\frm_Options.cs:line 190             at System.ComponentModel.BackgroundWorker.OnDoWork(DoWorkEventArgs e)             at System.ComponentModel.BackgroundWorker.WorkerThreadStart(Object argument)      InnerException:


title: Solution for web mapping

text: I have already built a desktop application for the management of archaeological settlements. Currently the data are stored in shapefiles (spatial data) and MYSQL (descriptive information). The application is based on MapWindow and MapWindow plugins. In the database there are some relationships between tables (one-to-many, many-to-many). I need to convert my application to a web mapping application. I intent to share my spatial data with geoserver.The web application should have forms to search the descriptive data of mysql and tools to select and zoom to the corresponting spatial data. Almost all the existed web mapping solutions have query capabilites but only on the spatial data (i.e. fields of a shapefile). Unfortunately my data are stored separately in mysql. What's the best solution to implement my web mapping application? Do I have to make custom tools? I intend to use Geomoose and some special plugins or Mapbender. Are they the appropriate platforms? I need a platform to join the spatial data and my mysql data in a easy and usefull way. I try to avoid extra programming for the forms and the connection between spatial and descriptive data. Thank you, Leonidas


title: Join Zip Codes Table to Nearest City from Second Table

text: I have two tables, siteareas and ziplookup. The former contains large cities in the US; the latter all US zips with lat/lngs. I have a foreign key into the former table in the latter: sitearea_id. I have made the join for all cities in the zip table that area present in the siteareas table. I wish to do so for the remaining zip codes by finding the nearest city from siteareas and updating ziplookup with that city's sitearea_id. Here's what I have so far:               update ziplookup set sitearea_id = siteareas.id      FROM siteareas           inner join ziplookup           on ST_DWithin(siteareas.lnglat, ziplookup.lnglat, 0.1)      where ziplookup.sitearea_id is null      order by siteareas.lnglat <-> ziplookup.lnglat      limit 1 ;      This is obviously not working, but I need a way to order by the closest result and limit to one result. This is probably more of a SQL question, but any help is appreciated.


title: Wrong coordinates after resizing the map

text: I catch a click event on openlayers map and show a popup with featureinfo (if any) with position set as `map.getLonLatFromPixel(new OpenLayers.Pixel(e.xy))`. The user might set the map to show in fullscreen with a button, a little of JS code and some css styling:               #geomap_f, #mapbox_f {         background: white;         position: fixed !important;         left: 0;         top: 0;         width: 100% !important;         height: 100% !important;         z-index: 49;     }      The problem arises when the user clicks the fullscreened map. The geoserver throws `org.geoserver.platform.ServiceException: 979, 8330 not in dimensions of image: 1017, 943` and e.xy object holds some really huge numbers (979, 8330 in this case). I know I need to tell the map or the geoserver (or both) that the mapsize was changed and that I want them to recalculate the relationship between xy and latlon. I've tried using some of OpenLayers.Map methods with no success so far though. How can it be done? Do I need to resize the map differently? I already call map.updateSize() after resizing, that doesn't seem to influence the coordinates.


title: GDAL 1.9.1 ogr2ogr KML to PostGIS not incuding ExtendedData

text: I have a KML file utilizing a Schema in the document and ExtendedData in each placemark. The ExtendedData does not seem to pass in conversion to PostGIS or shapefile. Any assistance would be much appreciated.


title: Correlation between point feature (disease) and elevation

text: I have a problem! How can I calculate the correlation between point feature and Dem raster? Hypothesis: disease prevalence in high altitude is more than low altitude. I must do this in ArcGIS 10 and classify the altitude to 15 or 20 classes. For example 0-200 meter above sea-level, 200-400 m and ... . I will join these altitude classes with disease prevalence shapefiles and then calculate the correlation between two entities. When I look at the layers, I see the correlation between altitude and disease.


title: ST_Equals Postgis problems

text: I am trying to find overlapping lines using ST_EQUALS on two tables This query returns 0 rows:               SELECT * FROM geom1 c, geom2 h WHERE st_equals(c.simple_geom, h.simple_geom);      This query runs correctly and returns overlapping lines:               SELECT * FROM geom1 c, geom2 h WHERE st_astext(c.simple_geom)=st_astext(h.simple_geom);      I am not sure why ST_equals fails. Both colums are indexed, I thought it was due to geom being multilinestring, but after converting to simple_geom with just LineString (using St_geometryN(geom,1)) it is still giving unexpected results. So it is not indexes, and neither MultiLineString column. What are other possible sources of this behaviour?


title: import parameters for r.watershed - (watershed basin analysis program) tool from grass using Qgis 2.2

text: I want to run this tool from Qgis 2.2 but i don't know how to use the PARAMETERS to run this tool. ![enter image description here](http://i.stack.imgur.com/38UIg.png) The red frame are the PARAMETERS I don't know. Can someone to tell me what those PARAMETERS are and what they mean?


title: How to add another layer origin in my installed OpenLayers plugin?

text: Can I modify my installed Openlayers plugin files to add a new layers origin? If yes, what files?


title: What are the existing datasets of world maritime routes?

text: I have a list of cargo ports of all over the world and need to show some port- to-port trips on a map, and also to compute lengths of these trips. For that, I need a dataset of actual maritime routes (even approximative) between important ports preferably free. Is there something available?


title: Algorithm to merge adjacent raster datasets

text: I am trying to figure out how to merge adjacent raster datasets. Is there any clear algorithms (may be a pseudo code)?


title: How to analyze the signatures for a Landsat clasification with Sextante in Quantum Gis

text: I’m defining my training data for a Landsat classification. I´m digitalizing the polygons with the field data in QGis and I want to analyze the signatures of the Landsat images to avoid errors. In ERDAS use de Signature editor, but I think I can do it with the toolbox Sextante in QGIS. I can’t find in the web an example of this. Anyone knows?


title: Agent Analyst tool interface will not open on ArcGIS 10.1

text: I am having the same sort of problem as others seem to have had (for me specifically the agent analyst tool interface won't open). I have ArcGIS 10.1 sp1 installed on my computer. I have followed these steps so far:   1. added the refresh.config.exe text file (ArcGIS 10.1 using Repast's Agent Analyst Extension - refresh.exe does not work)   2. changed values 10.0 to 10.1 in files agent_analyst.bat and agent_analyst_run.bat in agent analyst installation folder 3- added the ArcGIS 10.1 installation path "for example C:\Program Files (x86)\ArcGIS\Desktop10.1" to environment variables ARC_JAVA & ARC_JAVAW & ARCOBJECTS_JAR (Which GIS software or plugins can perform a Cellular Automata Simulation?) I restarted and still have had no luck (just a black box flashes up when I create a new agent analyst tool or try to edit the tool)... I don't know if there is some other step I need to do or if I have done one of the above steps incorrectly. Any advice would be greatly appreciated!


title: Can I add an ArcToolbox tool to QGIS?

text: I have been using a publicly available tool in ArcGIS called SPreAD-GIS. It is a series of python scripts that are "implemented as an ArcGIS toolbox" (with tbx file). I am new to QGIS. Is there a way to incorporate and use SPreAD-GIS in QGIS? Any help would be appreciated. Thank you.


title: Qgis Spatial Join: problem with a column!

text: I run the function "Join attributes by location" in Qgis 2.2 between a point vector layer and a polygon vector layer. When I open the resulting table of contents of point layer there is a column of numbers (fields from 10 to 13) that fail, because there is the same value in every row. For example "-214", a value that does not exist in the overlay polygon shape. However, this column is the one I need for doing further process. Does anyone know what the problem could be? Any tricks? Perhaps problem of "encoding" systems is causing this? Additional information: I tried with Saga process but in Qgis 2.2 there is no option for Saga folder in the configuration menu. I tried to change encoding (UTF-8, System) but the results is the same! Carlo Z


title: Can Hyperlink be based on two fields for path and file name separately?

text: I need to create a hyperlink field. I have one field with the file name, and another field that contains only part of the file path (the real path is unknown). The table contains 14,000 records with 14,000 different DWG file names. Is it possible to create a Hyperlink field in a table or feature class that will start to search automatically with windows browser for the file at a specific folder according to the field with the estimated Path? i.e. one field with the file name to look for, and another field with part of the path to look at.


title: Is there a downloadable data source for mapping out rivers with names for the whole world?

text: I'd like to be able to extract rivers with their names and the flow size (not direction) for the whole world. The scale I'll be drawing it at is about 2 centimeters per 10 real world km - nothing too precise! At this point I've looked at the hydrosheds dataset available for download and it seems great except for missing the river names...do you guys know of any other that is available for download? I've also looked at the HYDRO1K (unable to download) and the Global Runoff World Rivers...I can't get a copy (asked already). Or - do you know how I could automate the naming of rivers after I have drawn them from hydrosheds? In other words - some online DB or API that I can use lat-long to ask for river name at this or that lat-long coordinate? I use Java and am expecting to have to code this out by using lat-long pairs... Thanks!


title: QGIS Web Client: Where to upload my .shp files?

text: I've installed QGIS Server and QGIS Web Client on an Ubuntu 12.04 Server according to these instructions: https://github.com/qgis/qgis-web-client When I load my own project with the Web Client, it shows the correct table of contents but no map, just the box that says "loading map". A WMS service of my project is running fine and the two sample projects (helloworld, naturalearth) are displayed by the Web Client as well. Looks like the client doesn't find the data for the project. On the server my folder structure looks like this: /qgis-web-client/projects/myprojectfile.qgs /qgis-web-client/projects/shapefiles/myfirstlayer.shp /qgis-web-client/projects/shapefiles/mysecondlayer.shp /qgis-web-client/projects/shapefiles/... The project was created with QGIS 1.8 for Windows using relative paths and this folder structure: /myproject/myprojectfile.qgs /myproject/shapefiles/myfirstlayer.shp ... Where do I have to put the "shapefiles" folder? Do I need to specify that location somewhere?


title: Somebody knows a test to my application using WMS 1.1 (as a client) with OGC Standard?

text: I need a test to certificate my application with the OGC Standard. Please if somebody knows a site or something in which i can do this please contact me or answer this topic. Thanks. My e-mail is cyberchino23@gmail.com Greetings...


title: is the GPS network free or available for any device?

text: I have a question (I'm really sorry if is kind of dumb question). But if I have a device with GPS (smarthphone,pda,pocket pc, etc) Does the device need any additional service to enable GPS data? I mean is the GPS network free or available for any device with the appropriate hardware ? I'm kind of lost and probably I need to use GPS for one application running on windows mobile 6.5


title: Is GPS more accurate on specific hours of the day?

text: I have heard that the GPS system is more accurate, and less accurate on some specific hours on the day. If this is true, how can I get the information on when GPS is more or less accurate at my location? And how big is the difference? I guess this depends on where in the world you are, I'm living very north, so it's maybe has a bigger effect here. I also think it has to do with how many GPS satellites you get signals from at the moment.


title: Clip raster to multi raster

text: I have a raster aerial photo. I have county wide parcel data. **I want to clip/screen shot the raster at each individual parcel**. Think of this application - a post card with an image of your house --here is your house-- I want to have each parcel attribute to contain information + hyperlink of the aerial photo. Any other solutions welcome, but this is the only method I could think of. Although I am unsure of how to achieve it efficiently. Thanks! Using Arcinfo. 10.


title: Is there a decent open source map-editor?

text: I would like to have an open source map editor for drawing a city map? I don't want to use OpenStreetMap. It will be a custom map project. I don't need an advanced GIS system, but a simple map-editor, where I can use my own custom symbols. It must support Transverse Mercator projection and be able to calculated distances between objects on the map.


title: How to compute the distance between centroids and polygons?

text: Suppose I have an Excel spreadsheet with the latitude and longitude of the centroid of every county in the U.S. I want to know the shortest distance from each county to any part of Oklahoma (presumably, somewhere along the border). I then want this distance in this Excel spreadsheet as a new column. Is this possible to do and how? I know next to nothing about GIS


title: Is it possible to sort by Layer Name when adding vector data from file geodatabase in QGIS?

text: Is it possible to order the data alphabetically when adding data from a file geodatabase in the _Select vector layers to add..._ window in QGIS? It is ordered by the Layer ID number, with no way to sort by Layer Name, Number of Features or Geometry Type. Shapefiles are sorted alphabetically which makes it easy to find the data I am looking for, but the Layer ID from the file geodatabase makes it very difficult to add the data quickly as many of the names are very similar (ie wpt2013_03_12_R22, wpt2013_03_13_R2) but have very different Layer ID numbers. It would be much more efficient to keep the majority of the data in file geodatabases so that it can be used effectively in both ArcGIS and QGIS


title: How to Make QGIS on Windows use Spatialite 3.0 instead of 2.4?

text: I spent all afternoon trying to figure out why my complex triggers are working fine in Spatialite-GUI but QGIS refuses to write to the layer. Finally remembered I ran into a similar issue a few months back. Problem is, there's no workaround this time. I'm using triggers to force a point to snap to the nearest line on another layer, using the midpoint if no geometry is supplied, or the nearest point if it is. Works great in 3.0, but sadly, almost all of the geo functions I'm using (Interpolate Line, Shortest Distance, etc) are not on the list of functions in 2.4. So... I know there's been discussion about removing the built-in Spatialite from QGIS. That seems to imply there'd be a way to point QGIS at Spatialite installed somewhere else. Does anyone know how to make that happen? Running Windows 7, QGIS 1.7.4. Sorry if it sounds stupid, but I've never compiled anything so I have pretty much no idea what I'm talking about.


title: What kind of spatial index does arcgis use?

text: My first master's research topic was about spatial indexes. I learned about R-Trees and R+-tree, R*-trees and other various but very similar topics. This may not be published information, but what spatial index type does ESRI use? As a secondary question can this type be changed and how? There are some situations where a different type of index would boost performance. In Oracle's case it can use the ST_Geometry spatial data type which can be indexed to either an R-tree or a quad tree or both: http://docs.oracle.com/html/A88805_01/sdo_intr.htm Other databases use binary storage mechanisms which are non-spatial by default and do not seem to use the spatial indexing available. In SQL 2005 and lower spatial datatype didn't exist, but ESRI was still able to use SQL databases. According to this case study SQL 2012: http://www.microsoft.com/casestudies/Case_Study_Detail.aspx?CaseStudyID=710000000547 it uses the automated spatial performance type, but doesn't list what kind of R-Tree it is. Bonus points for reading material answering this question. Thanks.


title: How to convert LAS file to numpy array?

text: I have begun learning how to manipulate LAS data in python and wanted to see how others handle LAS files. I would like to read the points (I am using a numpy array), and filter out classes 1 and 2 (unclassified and ground) to a separate array. I have the following code but cannot seem to get the points filtered.               # Import modules     from liblas import file     import numpy as np          if __name__=="__main__":         '''Read LAS file and create an array to hold X, Y, Z values'''         # Get file         las_file = r"E:\Testing\ground_filtered.las"         # Read file         f = file.File(las_file, mode='r')         # Get number of points from header         num_points = int(f.__len__())         # Create empty numpy array         PointsXYZIC = np.empty(shape=(num_points, 5))         # Load all LAS points into numpy array         counter = 0         for p in f:             newrow = [p.x, p.y, p.z, p.intensity, p.classification]             PointsXYZIC[counter] = newrow             counter += 1      I have seen arcpy.da.featureClassToNumpyArray, but I did not want to import arcpy nor have to convert to shapefile. Any suggestions on filtering or easier way to read LAS data into a numpy array are welcome.


title: Orthographic projection in Python Basemap

text: I'm trying to plot a GeoTIFF in Python using Basemap. I can do so correctly with `cyl` projection like this:               m = Basemap(projection='cyl',                 llcrnrlat=LL[1], urcrnrlat=UR[1],                 llcrnrlon=LL[0], urcrnrlon=UR[0])      But if I try to do so in the `ortho` projection, it does not work. I'm using this function call:               m = Basemap(projection='ortho',                 lon_0 = origin[0],                 lat_0 = origin[1],                 llcrnrlat=LL[1], urcrnrlat=UR[1],                 llcrnrlon=LL[0], urcrnrlon=UR[0])      Basemap / matplotlib crashes with the error:               AttributeError: 'Basemap' object has no attribute '_height'      If I remove the llcrnX arguments, the `ortho` projection works, but I see the whole globe. I guess this is OK, if I can then somehow plot my image using the corner coordinate. Basemap documentation states that the corners may be specified.


title: How to calculate the intersected area from a layer with overlapping buffers in qgis

text: I need to calculate the area of landusetypes in 50m around various dots (all in one shapfile). My problem is, the buffers around the dots are overlapping each over in some places. I tried to intersect/clip the layer (e.g. grasland) with the layer "buffer", but the overlapping parts are only counted once to the overlapping buffers not twice. Is there any tool in qgis which might help me? Or is there any other way to get this output? Thanks for your help.![enter image description here](http://i.stack.imgur.com/Mw5ww.jpg)


title: Create random points, only where it overlaps specific values in a raster

text: I need to create a specific number of random points THAT FALL ON A RASTER with values between 0 and 1. The trick is the raster has a lot of holes - it is a grid of linear features where there are height values on the linear features, but everything in between is a NA. I have selected all my used locations that fall on the linear features in the raster, and intersected them to get the height value. Now I need to create the same number of random points, restricted to being where there are raster values (i.e no points are generated in the holes, where there are NA values). Any advice? I've considered converting my raster to a polygon, however the raster is over a large extent (60 x 90 km) and is very high resolution (1m pixels), so this would probably take a long time.


title: Adding transportation schedules to a network dataset

text: How do I add transportation schedules to a multimodal network dataset? For example, Bus 1 runs between 08:00 and 16:00, and Bus 2 between 14:00 and 23:00. If a passenger wants to solve for a route starting at 16:30, Bus 1 cannot be included in the result as it is not in service at that time. After some research, I came across this article about managing a public transit system in GIS using ArcInfo 8 (the article is from 2004). The transit data model is essentially the same as the one that ESRI currently has available, except they added a table called `ScheduleTimetable`, which contains the bus schedules which are related to routes. This table can't be added to the network dataset though, as file gdb tables can't be stored inside a feature dataset (and therefore can't be used in a network dataset). How would I be able to access the time fields so they can be used as restrictions for the routes?


title: EPSG:4038 Gridset does not work with caching on GeoServer 2.2

text: I am trying to setup a cache for a GeoTIFF base layer. I set it up and I went to the built in GeoWebCache page where it generates test map for the caching. I saw the built in projection - EPSG:900913 - EPSG:4326 I tried these two links and it shows open layers map and it generates cached tiles properly. So I set to work to create a gridset for ESPG:4038 using this manual http://docs.geoserver.org/latest/en/user/webadmin/tilecache/gridsets.html#webadmin- tilecaching-gridsets (I added 20 levels of zoom) I added this gridset to the base map layer cache and now I see these three items in the GeoWebCache page.   * EPSG:900913   * EPSG:4326   * ESPG:4038 I tried the first two links and it works properly (showing the maps and I see the indicators that certain amount of disk space are being used by the cache layer). But when I try to open ESPG:4038, the generated open layer map is blank and there's no activity in the disk space indicator for the cache. When I check on Google Chrome, it return 256x256 transparent png tiles. I'd appreciate if anyone can point out which part I do wrong


title: Migrate python scripts from 9.3 to Arcgis10

text: What changes are needed to migrate scripts from the ArcGis 9.x environment/gp- object to the ArcGis 10/arcpy environment? Can anyone provide a good overview what the changes are? Can someone share pitfalls etc from their own experience?


title: Arcgis Visualbasic Templates

text: ![enter image description here](http://i.stack.imgur.com/rzc5x.png)I am using ArcGis 10 Template windows application(Desktop) in VisualStudio 2010. i have added an ArcGis Engine Map Control and engine TOC Control in that. when I am Adding my digitized map to the map control it is only showing my map not the digitized features and layers in map control. whereas in the TOC Control it is showing the name of all the layers. please help me.


title: how to create custom widget in mapguide fusion template?

text: I am using mapguide maestro 5.0.1. I want to add custom fusion widget which have option to search layer in multiple ways. I am not able to find any tutorial to add custom widget useful. can you please suggest any tutorial which explain everything in detail. I have working module of this option in ajaxviewer can I use it in fusion.


title: How to export and import ST_GEOMETRY from Oracle environment to another Oracle environment?

text: I have a geographic layer stored in an Oracle Geodatabase with ST_GEOMETRY type. I'd like to export that layer and its data into a SQL INSERT format in order to import it into another Oracle database which is in another server machine. After the export I get lines like this: `Insert into myTable(OBJECTID,NAME,SHAPE) values ('2989','MyName','SDE.ST_GEOMETRY(8,1120,672587.528700002,6824864.9106,676521.2799,6829005.2801,NULL,NULL,NULL,NULL,8002536.39765853,16898.2791219904,3,'oracle.sql.BLOB@5e105008'));` I don't know how to deal with the BLOB type. I suppose `oracle.sql.BLOB@5e105008`references a memory address where some data is stored. How can I export that belob and import it in the new Oracle DB?


title: Generate a kml file from Google results

text: Is it possible to generate a kml file of searched items in Google Maps. For example:   * Search for highschools in Virginia   * Return the results as a kml file with the schoolname address and location


title: Offline map(google/openlayers) for HTML5(Angular) based hybrid application

text: I have a requirement where the customer has explicitly asked for offline support for the map. The rest of the application works well offline, but i am not sure how to get this done using **Google** or **OpenLayers**. I have gone through here and i couldn't find any answer to what i was looking for. Please do note that i am only using `angular.js` and `phonegap.js` and nothing else.I have worked with google maps as well as openlayers and im inclined to using these as i feel its much light weight and community support is awesome. **Q. Has anyone implemented or come across offline usage of google maps or openlayers. Other reliable map solutions are also welcome, but just need to be fast and mobile friendly.** What i have done : I have downloaded the entire source repo from github(openlayers) and included in my application. Although it takes hell of a time loading in browser, the application when installed had no performance issue. The issue was even though i included all files in my app, openlayers still fetches the map images(tiles) from another server and it is like dynamic(I had asked related here : Openlayer maps issue with ssl).               http://c.tile.openstreetmap.org/14/8745/5622.png     http://c.tile.openstreetmap.org/14/8745/5625.png     http://b.tile.openstreetmap.org/14/8745/5623.png      and so on. I have not gone with google maps as everything i got is direct api call to google server. Any good leads will be greatly appreciated.


title: Can I transfer custom QGIS settings from one computer to another?

text: Can I export custom general settings from my QGIS installation to a file so that I can transfer the settings to QGIS on a different computer? Does QGIS have them already saved to a file located somewhere on hard disk? Thanks!


title: How Do I configure PostGIS 2.0 on windows to use ST_GeomFromGeoJSON?

text: I have installed PostGIS 2.0 on a Postgres 9.1 database on a Windows 2012 Server. I need to write an insert query, which uses ST_GeomFromGeoJSON(). The Docs for this query say that I need - JSON-C >= 0.9, and very helpfully give a link on how to configure PostGIS for that. Unfortunately, the article which talks about configuration, only shows how you configure PostGIS on a Linux system. What steps do I need, to configure PostGIS 2.0 on windows to use ST_GeomFromGeoJSON?


title: drupal cartaro ask me postgis 1.5

text: I would install drupal cartaro, I have a host service with postgis 2.0 and a service with geoserver installed. Cartaro recommend me to use postgis 1.5. Do you know some solution for to use pg 2.0? I wouldn't change my database installation.


title: Examples of Beautiful Maps (at a county level or closer)

text: The Examples of Beautiful Maps question was outstanding and it gave us some great results. I was wondering if anyone else had similar examples, but at a county or parcel level. At those levels, I wonder if people will react the same as they did to the Examples of Beautiful maps, or if there will be other facts to consider at that resolution. Thanks to everyone in advance.


title: OpenStreetMap shp layer is shifting

text: I am dealing with this issue more than 2 month and i am hoping that some one have a lead for me. I am using Openlayers to show some vectors layers like buildings, points, gardens and more layers that i got from the government data, and i got shp file of all the streets from OpenStreetMap. I want to show everything together so the streets will suite to all the other vector layers. I am using Geoserver to show all the layers with WMS, my layers projection is ITM (Israel) and the OpenStreetMap layers is the google projection. when i show everything together the street layers is shifting little bit and its look like this: ![It is not suite](http://i.stack.imgur.com/oEePk.png) as you can see the layers are not as it should be... and its happening also with the others layers from openstreetmap like buildings etc. I searched all over the internet and i didnt found solution... hope some one will help me...


title: How to configure Geonetwork to be accessed from remote machines?

text: I just installed Geonetwork and have been playing around with it on our server (localhost:8080), now I want to access this from other machines on the network but when I type in the address I get a page not found error. I found this link through googling and altered the jetty.xml file (I couldn't find a mapServices.xml file). As a result I can now get to the page but none of the mapping components work (they all still refer to geoserver on localhost). Is there a location where I can make this change once, instead of going through heaps of file and individual changes? Thanks Andrew


title: Open a location in Bing Maps 3D

text: Does anyone know how I can link to a location in Bing Maps and have it open in 3D mode? A link formatted like this: http://www.bing.com/maps/default.aspx?lvl=15&cp=37.9609257054894~-87.3786394334287 will open that location using their standard 2D map. Is there anything I can add in to have it open in 3D instead? Thanks, Seth


title: How to catch when layers visibility change in ArcGIS?

text: I'm trying to regenerate a Legend in a Pagelayout whenever the layers visibility are changed in the TocControl. What is the best way to catch when a layer visibility is changed?


title: Compare fields of two tables and calculate field

text: Can you help me to compare fields in two tables? I have two tables, which contain ID and Measure fields. I need to take stroke by stroke values from field1(table1) and find all values in field1(table2) which are less than sought-for value. And then take largest of them and put it on stroke field2(table1). Example showed on picture below. It's will be great, if it would be python script. ![enter image description here](http://i.stack.imgur.com/ECQbX.jpg)


title: Google earth balloon not fully visible

text: When i load my KML on google earth, the balloon appears whole but when I embed on a webpage using the embedKML tool, the data in the balloon (In this case a table) appears to be cut. The close button is not showing. I have tried increasing width and height in the tool but there is no change. ![Image of my google earth after embedding it on a webpage ](http://i.stack.imgur.com/Yx7S8.png)


title: Center (and add points too) a map with LonLat

text: I'm new to OpenLayers, so I'm sure this will be an easy win for someone :) I'm trying to center a map on a point (and place markers) using lon/lat but it appears to take everything as meters?                   map = new OpenLayers.Map("map",{         projection: new OpenLayers.Projection("EPSG:3857"),         units: "dd",         numZoomLevels: 18,         maxResolution: 156543.0339,         layers: [             new OpenLayers.Layer.OSM("OpenStreetMap", null, {                 transitionEffect: 'resize'             })         ],         center: new OpenLayers.LonLat(-2.626944, 57.563222),         zoom: 7     });     var markers = new OpenLayers.Layer.Markers( "Markers" );     map.addLayer(markers);      N.B. this code isn't designed to add markers, I will need to but the problem is clearly with this chunk as I'm being centered (almost) at the 0,0 I've tryed the following code               var size = new OpenLayers.Size(21,25);     var offset = new OpenLayers.Pixel(-(size.w/2), -size.h);     var icon = new OpenLayers.Icon('http://www.openlayers.org/dev/img/marker.png', size, offset);     var i=0;     function Go() {         x = document.getElementById("x");         i=i+10;         x.innerHTML = i;         markers.addMarker(new OpenLayers.Marker(new OpenLayers.LonLat(i,10),icon));     }      Calling go on a mouse click, and indeed my marker moves a small amount (I'm guessing 10 meters) to the right every click.


title: Select features from attribute table and then create a new field with the selected attribute data

text: I have a shape file which contains land use typology for my administrative area. I’m investigating deficiency in certain land uses and there are 2 ‘standards’ that we have to use. These standards use different land uses. The Shapefile has a field in the attribute table that contains all the land uses, for example each polygon can have a land use of either A1 – 12, B1 – 12, C1-12 etc. One standard requires land uses A1 – 6 and C5 -9 and the other A1-12 and B1-9. The Land uses A, B & C have separate fields. I have created an SQL query in the advanced search that selects these land uses. What I now want to do is create 2 new fields, one which populates the selected land use for each standard. I know I can select the attributes and save as a new selection but ideally I would like to populate 2 new fields with the attributes from the existing land use field.


title: Fastest way to add shapely geometries to a PostGis database using geoalchemy2

text: I am looking for the absolute fastest approach to adding +10,000 shapely geometry objects to a PostGIS database using geoalchemy2. I'm having no speed issues when uploading to a localhost postgres server, but when I upload directly to my remote server it takes 35x as long! I realize this might be a server question and not appropriate for this site, but curious if anyone has a successful workflow for: **Quickly uploading shapely geometries to a PostGIS database using geoalchemy2**


title: Multiple ArcGISDynamicServiceLayer in 1 map instance

text: Hope you can help me with this. We are trying to add the functionality to load multiple ArcGISDynamicServiceLayer to map instance as sub maps. Currently, it is loading fine but there are times that the layers are null. Am i doing it wrong? The process is generally creating a new instance of ArcGISDynamicMapServiceLayer with a specific map url and then add it to the current map instance(ESRI.ArcGIS.Client.Map). I'm not sure if this is the most efficient way of doing this. Is there another way in implementing this? We are using ArcGIS 10.1 in a Silverlight Viewer. Thanks.


title: New Zealand Projection System

text: I have a database with Easting and Northing points of New Zealand. The data structure is as follows Easting | Northing 2595425 | 6622565 2594585 | 6622595 2595000 | 6622085 2594500 | 6622005 2595300 | 6622915 I need to create shapefile from this data. My Final goal is to convert this data to WGS84 system.


title: How hidden/not include particular layer getting added in Tree panel?

text: I am trying to show tree panel using GeoExt js 2,Ext js 4 and Openlayer js 2.12. Now i want some/one/particular layer in map/mappanel hidden/ not to be include in tree panel. my code for store and tree panel is as follow:               var store = Ext.create('Ext.data.TreeStore', {                 model: 'GeoExt.data.LayerTreeModel',                 root: {                     expanded: true,                     children: [                          {                             plugins: ['gx_baselayercontainer'],                             expanded: true,                             text: "Base Maps"                         }, {                             plugins: ['gx_overlaylayercontainer'],                             expanded: true,                     text: "Overlay Layer",                     component: {                     xtype: "gx_wmslegend",                     layerRecord: mapPanel.layers,                     showTitle: false,                     cls: "legend"                 ]                 }             });                       // create the tree with the configuration from above             tree = Ext.create('GeoExt.tree.Panel', {                 border: true,                 region: "west",                 title: "Layers",                 width: 150,                 split: true,                 collapsible: true,                 collapseMode: "mini",                 autoScroll: true,                 store: store,             enableDD: true,                 rootVisible: false,             lines: false             });       I have used " **displayInLayerSwitcher:false** " but i thing it is for layerswitcher control in openlayer. Please help me for this.


title: GeoServer - Show Maps in a Application Java

text: How can I display a GeoServer's Map inside a Java's Desktop Application? I would like to display it inside a Java App instead of a browser. Is this possible?


title: Creating square buffer in qgis

text: How can I create a buffer with square ends along a line segment in qgis? I've buffered lines using v.buffer in the GRASS plugin, but the -c option for square ends is not exposed in that interface. Is there another way to do this with GRASS, or another tool? Thanks! update: I ended up using the grass qgis plugin with the "shell" module and just used the grass commands. If you want to see the results in qgis you can use the "add grass vector layer" tool.


title: How can i add a north arrow in QgsComposition?

text: I'm trying to add an svg image of north arrow in QgsComposition using this code:               mapRenderer = iface.mapCanvas().mapRenderer()     theComposition = QgsComposition(mapRenderer)     theComposition.setPlotStyle(QgsComposition.Print)     northarrowIcon = QgsComposerPicture(theComposition)     northarrowIcon.setPos(QPointF(50,110))     northarrowIcon.setPictureFile("./svg/north_arrows/north-arrow_3_simple_symmetric_triangular.svg")     theComposition.addItem(northarrowIcon)      i don't have an execution problem but the arrow is not shown. Can someone check and give me an idea what can be the cause of that.


title: how do you manage columns in an ArcSDE/Oracle feature class?

text: Our enterprise has a long history with Oracle and ArcSDE. And lots of tension between the Oracle people and GIS people. Oracle folk often don't like the requirement of using expensive 'nonstandard' tools (e.g. anything other than TOAD or SQLPLUS) to manipulate business tables. But I digress already! So for example today I wanted to do something simple - a column in a 15-year old non-versioned low-precision sdebinary feature class was marked as non- nullable. Recent editing clients don't use that column anymore, but for various reasons we don't want to drop it. **So I wanted to change the column to nullable**. 1) do in in TOAD easy peasy. Even if the feature class is locked by ArcMap sessions or ArcGIS Server map services. That worked. BUT - this does not update the geodatabase system tables. **SDETABLE -o describe** shows consistent with with what TOAD shows, but ArcGIS still sees the column as non- nullable. The **column_registry.object_flags** column (see Appendix E, Administering ArcSDE for Oracle) shows '4' for this column. Whatever that means. 2) do it in the SDE commandline tools - **SDETABLE -o alter_column**. That worked too, but first required me to identify locks on the feature class, track down which ArcGIS Server services include layers connecting as those users, stop them, then run the command. Lots of literal running around. AND - this also does not update the geodatabase system tables, so ArcGIS still sees the column as non-nullable. 3) use ArcCatalog 10.1 feature class properties dialog? In my experience everything it tells you about an ArcSDE column is read-only. So I guess not. 4) a Toolbox tool? I didn't find a handy 'change column metadata' tool in 10.1. 4) anything else I missed? What will you do when the commandline tools go away? Are things any better in 10.2? thanks for feedback, MC


title: How to get data from TileMill tileLayer properties in Mapbox for click/hover

text: I've successfully used Leaflet to load a geoJSON polygon-shapefile with data, and added click and hover interactivity to the layer. For a large dataset, I switched to mapbox.js 1.6.1, used TileMill to export a tiles, used infostream's mbtiles-php (https://github.com/infostreams/mbtiles-php) to serve the data. I'm very confused, now, how to use mapbox.js to access show the polygon attributes in the tile layer? I have successfully loaded the tiles onto a page using mapbox.js - but don't know where to start to pull in the attribute table. I did not find an appropriate Mapbox example on the example page, and the Mapbox support links are all directed to their contact page. Did I misunderstand how the purpose of Mapbox/TileMill?


title: Project twitter geo data into wgs 84 map

text: I'm currently working on displaying twitter geo data (currently just long/lat -> cartesian without any analysis). This is the world map (source: naturalearthdata): http://dl.dropbox.com/u/14747311/world.png This is the rendered data: http://dl.dropbox.com/u/14747311/overlay.png. The mapping is done via mercator projection (with R = 1):               x = ((long + 180)/360)*x_dim     y = (y_dim/2) - ((arsinh(tan(lat))*(y_dim/2) (lat converted to radiant of course)      Now, the shape is apparent, but of course, the projections don't match. How can I fix that? (programming experience is plentiful, so it's more the theoretical side I'm concerned with).


title: Qspatialite plugin problems with queries using inner join

text: I'm trying to create a layer in QGIS from a sqlite database using a query from 2 tables: one has spatial data and the other doesn't. When I build the query using Qspatialite using the 'Advance SQL Editor' it throws an error saying that the first column in the query has an ambiguous column name. Copying and pasting from an MSAccess query, it finds the matching records without a mutter. Here are segments of the two queries: MSAccess version (working)               SELECT Inspections.ID, Inspections.InspectionDate, Inspections.InspectionBy, Inspections.PathID, Inspections.Section, Inspections.GISLink, Inspections.Observations, Inspections.Notes, Inspections.ActionPriority, Sections.Geometry AS Geometry     FROM Inspections INNER JOIN Sections ON Inspections.GISLink = Sections.GISLink;      Qspatialite version (not working)               SELECT 'Inspections'.'InspectionDate',     'Inspections'.'InspectionBy',     'Inspections'.'PathID',     'Inspections'.'Section',     'Inspections'.'GISLink',     'Inspections'.'Observations',     'Inspections'.'Notes',     'Inspections'.'ActionPriority',     'Sections'.'GISLink',     'Sections'.'Geometry' AS Geometry     FROM 'Inspections', 'Sections'     inner join 'Inspections'     on 'Inspections'.'GISLink' = 'Sections'.'GISLink'      Could anyone give me an idea of what wrong's here - when I remove the INNER JOIN the query runs, but obviously gives me the duplicate values you'd expect! Is this a problem with QSpatialite in the way it constructs the queries with '' all over the place? I've just tested both queries in QGIS DB manager and same result the copy- paste version works, but the 'Advance SQL editor' version throws the same error... Many thanks - and hoping that I can save someone else frustrating hours of SQL - I was assuming the differences in MSAccess would prevent the query from running in Qspatialite and as a last ditch I did the copy-paste. Doh! I'm using QGIS 1.8.0 and QSpatialite 6.0.4


title: How do I know the vertical/horizontal units of my DEM? I'm creating a slope using gdaldem

text: I'm using gdaldem to create a slope map. I can apply a scale. The default is 1. The notes say:               Ratio of vertical units to horizontal.     If horizontal unit of source DEM is degrees (e.g WGS84 projection), use scale=111120     if the vertical units are meters (or scale=370400 if they are in feet)      Can someone explain this to me in layman's terms. My DEM is in NAD27 UTM 18N and in the projection properties it says linear units Meters and angular units degrees. How do I know what the vertical units are and what the horizontal units are? Thanks.


title: Open attribute table without highlight selected feature in ArcObjects

text: I am using ArcGIS 10.2.1, I have a need to display the attribute table only for selected features, but I don't want to highlight the selected features. The purpose of open the attribute table is to view feature information and after that allow user to zoom to the feature. Is it possible in ArcObjects? Or is there any alternative solution instead of open attribute table?


title: Tile server with dynamic mapnik rendering, caching, and pre-emptive rendering

text: I've been making some pretty (but computationally expensive) tiles with TileMill and I'd like to serve them up. First, I served them straight from TileMill. Every tile gets rendered dynamically, which is simple, but slow. Next, I exported a region as MBTiles then served it with MBTiles-server. Every tile is served statically, which is fast, but exporting all the tiles is too slow (multiple days). I'm looking at TileStache, which combines dynamic rendering with caching. I imagine that this means the first time a region is looked at, it will be extremely slow, then quick thereafter. I have plenty of spare CPU time (8 cores, not being used for anything else) and plenty of disk space (terabytes). I'd really like a server that can:               a) serve cached tiles if available     b) dynamically render tiles if necessary     c) continuously render tiles in the background to maximise the cache hit rate      Is mod_tile the only option? I'd prefer to avoid Apache, and mod_tile's documentation looks messy. For context, I want to serve tiles for the whole of Australia, down to zoom 15.


title: How to group and count attribute data?

text: I have a problem with group same data to count it. for example, i want to group same code to one and count them. i can do it when i use MapInfo (Pic 1) but i can't in QGIS. (Pic 2: i can classify data and count it but can't get them in attribute table.) It's possible to do it in QGis? please tell me. Thank you. ![enter image description here](http://i.stack.imgur.com/yb6Kr.jpg)


title: how to create a Heat-Map in ArcMap 10.1

text: Is there any tutorial about creating a heat-map in ArcMap 10.1?


title: Export a large raster to a comma or tab delimited text file in ArcGIS 10.2 for Desktop

text: I have a large raster with multiple attributes in 200-meter cells, and I would like to export this to a tab or comma delimited file which contains the X,Y coordinates and all of the attributes. If it's possible I'd also like to have the file arranged so it goes East to West across the whole raster like a typewriter. Thanks for looking and thanks in advance for any help you can provide. I'm currently in the process of trying raster to point, assigning X,Y, to the points and exporting the table.


title: Problem when converting Shape file to KML files in QGIS

text: Problem when creating KML files in QGIS? I select Layer CRS as WGS 84 (EPSG:4326) When I open the kml layer in Google earth it is wrong. In Google earth, appears some segments every where on the globe. Can anyone explain to me what is the wrong?


title: Add points to shapefile and display over map interactively using Python Add-Ins

text: I want to add points to a shapefile interactively using pythonaddins and display them over map. Points are added to shapefile but not to display. Please suggest how to fix this discrepancy.               class AddPoints(object):        """Implementation for points_addin.tool1 (Tool)"""        def __init__(self):            self.enabled = False            self.shape = 3             def onMouseDownMap(self, x, y, button, shift):            xy = (x,y)            cursor = arcpy.da.InsertCursor(r'C:/temp/points.shp', ("SHAPE@XY"))            cursor.insertRow([xy])                 arcpy.RefreshTOC()            arcpy.RefreshActiveView()


title: Crowdsourcing: Alternatives to Ushahidi

text: We all know that Ushahidi offers a versatile solution for crowdsourcing information and its effective depiction in maps. I am interested to know if there is any application(free and open source) on par with or superior to Ushahidi! I would be using the application to collect data from users via tweets, web forms, emails, SMS and map them using the location information provided by them; all of which can be accomplished by Ushahidi. My deployment would be a permanent online deployment to collect ongoing information and not the one which would be live for few days/weeks. Hence, one of the major limitation I foresee with Ushahidi is the requirement of manual intervention for approval of reports, which may need more personnel in the event of large number of responses. Is there any possibility to build a custom solution for such purpose using different open source tools?


title: postgis : is the operator <-> works for Geography also?

text: With recent postgres//postgis version, you can now use the operator <-> to find nearest rows to a point : > Indexed Nearest Neighbour Search in PostGIS Correct me if I'm wrong, but according to my tests :               SELECT name, gid     FROM geonames     ORDER BY geom <-> st_setsrid(st_makepoint(-90,40),4326)     LIMIT 10;      is the same precision and a bit slower as if I use the Postgres Point type. Now I would like to use this new operator with **Geography** type to do better calculus on a spheroid. Is it possible ?


title: Import MySQL Spatial vector layer in QGIS

text: I have a spatial table called 'places' in MySQL, with a spatial column 'pt' (point). I tried to import it in QGIS and no points appear. I created a view 'places_as_text' which reproduces the table but shows the points as WKT (used astext(pt) in the definition of the view). Still nothing appears in QGIS. Any help appreciated


title: Excluding layers from snapping using IPointSnapper

text: I'm working with ArcObjects in C#, and am trying to exclude all map layers but one (the only one I want to snap to) from the snapping environment. I'm using the ISnappingEnvironment and IPointSnapper interfaces. I've put the following code in the OnClick method of the tool I'm using:                m_Snapper = m_SnappingEnvironment.PointSnapper;      ISet anExclLyrSet = new SetClass();      m_Snapper.ExcludedLayers(ref anExclLyrSet);           // ...loop through all map layers, adding all but the one of interest to anExclLyrSet           m_Snapper.ClearCache();      Basically, I'm getting poor snapping performance due to the presence of a basemap layer that contains multiple feature layers. If I turn that layer off (make invisible), the performance is good. I've confirmed that all the layers (for composite type layers, I add each sub- layer to the exclusion as well) are being added to the exclusion set, but I can't seem to get the performance to improve without turning off that basemap layer. I'm not sure what I'm doing wrong - does the cache need to be cleared after each individual layer is excluded? I'd appreciate insights anyone might have regarding this. Thanks! Mark C.


title: Identifying my coordinates from external datasource

text: I'm working with data received from an agency outside of my office and two columns in the dataset are labeled 'CompX' and 'CompY'. Here is a sample set pulled from the spreadsheet:               +-----------+------------+     |   CompX   |   CompY    |     +-----------+------------+     | 584490.53 | 4513000.14 |     | 584565.52 | 4513137.77 |     | 583596.45 | 4509846.2  |     +-----------+------------+      This data should all be in New York City, however, using these as decimal lat longs (like 58.449053 & 45.1300014) shows up in eastern Europe. Anyone familiar with this type of coordinate system? I'm trying to create a map based on this data. Thanks!


title: How do I implement outbound arcobjects interfaces in C#?

text: I would like to write a class in C# that implements one or more **outbound** interfaces, for example ILayerEvents and IActiveViewEvents. Does anyone know how to do this with C#? **Update** Following @Chaz's response, I have tried this out.               class MyCustomLayer : BaseCustomLayer, ILayerEvents_Event     {         public override bool Visible         {             get { return base.Visible;}             set             {                 base.Visible = value;                 if (this.VisibilityChanged != null)                     this.VisibilityChanged(value);             }         }         public override void Draw(esriDrawPhase drawPhase, IDisplay Display, ITrackCancel trackCancel)         {             //throw new NotImplementedException();         }         public event ILayerEvents_VisibilityChangedEventHandler VisibilityChanged;     }      I expected that when the map sees that a layer being added implements ILayerEvents_Event, that it would wire events accordingly. However, when I add my custom layer to the map this doesn't happen. I wrote a test command that toggles `ILayer.Visible` but the checkbox in the TOC doesn't change ... `this.VisibilityChanged` is always null, indicating that there are no subscribers. Or maybe I'm wrong in assuming the TOC somehow subscribes to this interface?               protected override void OnClick()     {         MyCustomLayer l = null;         for (int i = 0; i < ArcMap.Document.FocusMap.LayerCount; i++)         {             var layer = ArcMap.Document.FocusMap.get_Layer(i);             if (layer is MyCustomLayer)                 l = layer as MyCustomLayer;         }         if (l != null)             l.Visible = !l.Visible;         else         {             l = new MyCustomLayer() { Name = "a custom layer" };             ArcMap.Document.FocusMap.AddLayer(l);         }     }


title: Does OpenLayers have the ability to mask a raster with another layer?

text: Does OpenLayers have the ability to mask (in the GIS sense) one raster layer by another vector (or raster) layer? For instance this would be the ability to show only the parts of one layer that are spatially coincident with another.


title: Does ArcMap have a Pan To Layer command?

text: In an Attribute Table, I am able to choose Pan to Feature (rather than Zoom to Feature) to refocus the map without changing the zoom level. Is there a similar Pan to Layer option (and/or shortcut key) when in the Table of Contents, to match the Zoom to Layer? (The workaround is, obviously, Zoom to Layer and then reset the scale -- I'm just looking for a simpler way.)


title: Pull out multiple elevation points from an ISurface

text: I have an `ISurface` object and wish to pull out several elevations with one call; instead of using the `GetElevation(IPoint)` method. The ISurface is representing a dem with 30m cell size. Is there an `ISurface` method that pulls out several elevations with one call? Is there another way of pulling elevations from a dem? A query?


title: stacked colours according attributes

text: I have a hierachical List of attributes e.g. column (1) S (2) SD (3) SD0 or (1) A (2) AN (3) An1 etc. So there is one S with several SD s and several SD0 in the attribute table. When I display column (3) I wonder if there is an automatic way to set the colours in a stacked way i.e. all S different grey colours all A different blue etc. All I could find was the rule based way but then it is still a lot of work because there are so many types in column 3. Is there a way to use the other columns, in my case for example column 1 to classify? Thank you for help!!


title: How to create polyline feature from two input points by following/tracing another polyline?

text: I have a list of things each have two grid references, start and stop points. I can turn these into points and join them with polylines. But what I need to do is have the lines that join the points the follow the route of another polyline, in this case it is a river. I have 38 lists (for various rivers), each with 10-20 items on them, so hopefully I can find a quick soloution.


title: How can I enter metadata into a qgiscloud project

text: I published a project into qgiscloud: https://qgiscloud.com/tragewegen/IDEAproject_DonkeyMobileBlackCarbon. I want metadata be published in the metadatatab next to the legendtab. Can somebody tell me how to do that? Thanks in advance


title: Fixing the OpenLayers MBTiles Example

text: I've been trying to implement the OpenLayers Example for reading MBTiles: http://projects.bryanmcbride.com/ol_mbtiles/ Right off, the link to the necessary mbtiles.php script was broken, but I found what I believe to be a copy of it, here: https://gist.github.com/1818011 After exporting my MBTiles file from TileMill, I saw in the SQLite that the mbtiles.php file's SQL query was using the wrong columns and tables for pulling the appropriate tile images. The original query:               $sql = "SELECT * FROM tiles WHERE zoom_level = $zoom AND tile_column = $column AND tile_row = $row";      What it should be based off what TileMill Exported:               $sql = "SELECT images.tile_data FROM map, images WHERE map.zoom_level = $zoom AND map.tile_column = $column AND map.tile_row = $row AND map.tile_id = images.tile_id";      Once I did that and stopped doing bindColumns on SQLite columns I wasn't using, I could pull the tile images from my MBTiles file using the php script. That's where I get stuck. I followed the Javascript code from the OpenLayers example linked to above, but the TMS layer keeps calling for tile rows and columns that don't exist. Here's the relevant bits of my JS code. What I'm not showing is my creation of an OSM base layer. I think the problem lies in the math used to calculate the x and y values in mbtilesURL():               function load_mbtiles(){                 var mbtilesLayer = new OpenLayers.Layer.TMS("MBTiles Overlay", "mbtiles.php", {                     getURL: mbtilesURL,                     isBaseLayer: false                 });                      map.addLayer(mbtilesLayer);             }                  function mbtilesURL (bounds) {                 var db = "filename.mbtiles";                 var res = this.map.getResolution();                 var x = Math.round ((bounds.left - this.maxExtent.left) / (res * this.tileSize.w));                 var y = Math.round ((this.maxExtent.top - bounds.top) / (res * this.tileSize.h));                 var z = this.map.getZoom();                      return this.url+"?db="+db+"&z="+z+"&x="+x+"&y="+((1 << z) - y - 1);                  }      Can anyone think of why this function would be calling the completely wrong tiles?


title: Auto fill column in attribute table

text: Is there a way to fill a column in the attribute table with the same value, similar to a paste . Is there something in the calculator that can be used to fill selected records/attributes.


title: How to show point labels with rotation and offset

text: > In ArcGIS 10.2, I have a point layer with street names and rotation field. I > need to use the rotation field AND show labels with an offset and centered > on top of the points (I want all names to show centered and over the > points). > > Any way to do this ? > > It's very easy in QGIS, but cannot find a way in ArcGIS... only option is > either rotation OR offset but not both together ?


title: Problems with QGIS plugin traceDigitize on Mac?

text: I'm trying to get the traceDigitize plugin to work with QGIS 2.0.1 on my Mac while following this tutorial but I get this error as soon as I ctrl-click:               An error has occured while executing Python code:          Traceback (most recent call last):       File "/Users/oskarlin/.qgis2/python/plugins/traceDigitize/traceDigitize.py", line 135, in createFeature         if (geom.isGeosValid()):     AttributeError: 'NoneType' object has no attribute 'isGeosValid'          Python version:     2.7.5 (default, Aug 25 2013, 00:04:04)      [GCC 4.2.1 Compatible Apple LLVM 5.0 (clang-500.0.68)]          QGIS version:     2.0.1-Dufour Dufour, f738351          Python path: ['/Applications/QGIS.app/Contents/MacOS/../Resources/python/plugins/processing', '/Applications/QGIS.app/Contents/MacOS/../Resources/python', '/Users/oskarlin/.qgis2/python', '/Users/oskarlin/.qgis2/python/plugins', '/Applications/QGIS.app/Contents/MacOS/../Resources/python/plugins', '/Library/Frameworks/SQLite3.framework/Versions/B/Python/2.7', '/Library/Frameworks/GEOS.framework/Versions/3/Python/2.7', '/Library/Frameworks/GDAL.framework/Versions/1.10/Python/2.7/site-packages', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python27.zip', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat-darwin', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat-mac', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat-mac/lib-scriptpackages', '/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib-tk', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib-old', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib-dynload', '/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/PyObjC', '/Library/Python/2.7/site-packages', '/Applications/QGIS.app/Contents/Resources/python/plugins/fTools/tools']      Could this be a Mac specific problem or am I doing something wrong? I'm using OSX version 10.9 (Mavericks) with QGIS 2.0.1 and GDAL Complete framework.


title: How to solve error "insufficient privileges to query iis metabase" in ArcGIS server 10?

text: I want to create a web application in ArcGIS server 10 .I have Windows 2008 R2 and Log in as Administrator . I created a service in ArcGIS server successfully but when i want to create a web application using Web applications the error **Insufficient Privileges to query IIS metabase** does not allow to create a web .I surf in the internet such as this solution and another solutins and found some recommended solutions but they don't work. Is there any solution to solve the problem? the user Administrator is in agsadmin and agsuser groups. Thanks for Advance


title: QGIS save as SLD style for graduated styled layer different from what it's shown in QGIS

text: I've found that when applying a graduated style to a layer in QGIS, if I save the style as SLD to use it in my Geoserver, then the SLD rendered layer looks different. So I've made a small example: ![QGIS styled layer](http://i.stack.imgur.com/YRMZE.png) The layer contains 4 elemtents, with a feature value that is 1,2,3,4 for each.   * First color style class represents values from 1 (included) to 2 (included)    * 2nd style class represents values > 2 and <= 3 (right?)   * 3td style class represents values > 3 and <= 4 but the Geoserver WMS rendered result with the same QGIS style saved as SLD is: ![enter image description here](http://i.stack.imgur.com/Jj0zS.png) where the polygon with feature value = 1 is missing. In fact the exported SLD contains instruction to style values > 1 and <= 2 for the first style class:                     <se:Rule>           <se:Name>1.0000 - 2.0000</se:Name>           <se:Description>             <se:Title>1.0000 - 2.0000</se:Title>           </se:Description>           <ogc:Filter>             <ogc:And>               <ogc:PropertyIsGreaterThan>                 <ogc:PropertyName>somedata</ogc:PropertyName>                 <ogc:Literal>1</ogc:Literal>               </ogc:PropertyIsGreaterThan>               <ogc:PropertyIsLessThanOrEqualTo>                 <ogc:PropertyName>somedata</ogc:PropertyName>                 <ogc:Literal>2</ogc:Literal>               </ogc:PropertyIsLessThanOrEqualTo>             </ogc:And>           </ogc:Filter>           <se:PolygonSymbolizer>             <se:Fill>               <se:SvgParameter name="fill">#edf8fb</se:SvgParameter>             </se:Fill>             <se:Stroke>               <se:SvgParameter name="stroke">#000000</se:SvgParameter>               <se:SvgParameter name="stroke-width">0.26</se:SvgParameter>             </se:Stroke>           </se:PolygonSymbolizer>         </se:Rule>      and this can be a problem, as I need to hand modify all the SLD I'm producing with QGIS now. Also, the label shown in QGIS is not exported in the SLD as you can see the style do not render any layer label. QGIS in use is 2.01 on Linux. Is there a way to solve this?


title: QGIS 1.8 Project in 2.0.1 Bad Layers

text: I have a confusing problem with a project generated in QGIS 1.8 when trying to open this on a new windows pc with Darfour (same OS, Win7) The hard drive including the project has been moved from old machine to new and as such no paths etc should have changed. ![http://picpaste.com/qgisBadLayers.jpg](http://picpaste.com/extpics/qgisBadLayers.jpg) http://picpaste.com/qgisBadLayers.jpg I'd obviously rather not have to manually correct all 40, and want to understand exactly what has changed / what the issue is. The layers are UK OS data. Many thanks * * * Many thanks for your reply. The data in this case are in shapefiles from the OS Digimap service The issue was that Darfour expects the following               <datasource>./OS VectorMap District (Vector)SU/data/SU/SU_AdministrativeBoundary.shp</datasource>      and the original project are formatted thusly:               <datasource>./OS VectorMap District (Vector) SU/data/SU|layername=SU_Airport</datasource>      a couple of quick find and replace to take the layername out and add .shp did the trick


title: how to make public transportation map?

text: I have collected the public transport routes and i want to map a map that looks like this. My problem is i am not able to show the shared segments as parallel lines. I have not been able to find good resources on this. I have good experiences with making maps and GIS softwares. If this requires special tools , then i can learn those too. Although i am looking for a static map, any possibility of making an interactive map would be better.


title: Open all popups OpenLayers

text: I have some WMS layers with attributes and I want to show some info when I click on some of them. Currently I'm using:               info = new OpenLayers.Control.WMSGetFeatureInfo({url: 'http://localhost:8080/geoserver/SIMI/wms',                      title: 'Identify features by clicking',                     queryVisible: true,                     eventListeners: {                         getfeatureinfo: function(event) {                             map.addPopup(new OpenLayers.Popup.FramedCloud(                                 "chicken",                                  map.getLonLatFromPixel(event.xy),                                 null,                                 event.text,                                 null,                                 true                             ));                         }                     }                 });      map.addControl(info);           info.activate();      Inside the init() function and works perfectly. When I click something it shows the info. But I want the popups to be opened as default, when I load the page containing the map. Is this possible? Also, I'd like to be able to open the popup only for a layer (I'm using a layer group and currently using content.ftl to filter out stuff but it's not viable since it still opens a popup with no content). Thank you!


title: Is a master's degree needed for a GIS career?

text: Currently I'm on the 5th year of education at the university (Dept. of geography, specialisation - "GIS"). In my country (Belarus) we have a system of education that differs from the majority of the other states. So after these five years "bachelor degree" is not written in a diploma (but it is implied). And there is a possibility to have a master's degree (one year more) that is (as far as I know) treated in other countries. I have an intention to leave the country as there isn't any possibility to find a GIS-related job here. Is the point of view, that in this case, it's much better to be a master than a bachelor?


title: Splitting strings into new rows, retaining geometry

text: **UPDATED BELOW** I have 9 shapefiles, each shapefile has 700,000+ rows and I have roughly 12 attribute fields in each shapefile. There are three attribute fields, common to every shapefile, that I am interested in - Type, Code, Area. In the majority of records, there is one entry for each field (| B | X2 | 3.45 |) but in some cases you have more than one entry in a field (up to 8) separated by comma; | A, B, B | X2, ZZ9, GGR | 1.2, 3.1, 0.23 | It is nearly always the case that if there is 3 records in one field, the other two fields concerned will also have 3 records. If this does not apply, then the data is incomplete and can be ignored. It is also the case that the 1st entry in one field is associated with the 1st entry in another field. I have found how to ArcGIS 10 Python split string and ArcMap 10.1 Python split inconsistent string but that results in deconcatenating the string by the comma delimiter and creating new attribute fields for every comma separated entry. The thing i need is for the deconcatenated data to create new records NOT new fields. In the instance above;   1. | A, B, B | X2, ZZ9, GGR | 1.2, 3.1, 0.23 | would become   1. | A | X2 | 1.2 |    2. | B | ZZ9 | 3.1 |    3. | B | GGR | 0.23 | the top one is the original record with the superfluous data removed into two new records below it. I need the new records to retain all the other information from the other 9 fields that dont need to be touched and also the geometry as they refer to the same polygon space. thanks very much, i'm stuck! (Answers for ArcGIS 10.1, FME, R, Python all gladly accepted) **EDIT 9th APR:** i used the following code to split concatenated entries (by a comma) in fields 'TYPE' and 'AREA' into new rows (there are between 1 and 8 concatenated entries in TYPE and AREA):               import arcpy     searchC =  arcpy.da.SearchCursor("FeatureClass", ("TYPE","AREA", "SHAPE@"))          insertC = arcpy.da.InsertCursor("newFeatureClass", ("TYPE","AREA", "SHAPE@"))           for row in searchC:        A = [x.split(",") for x in row[:-1] ]             for i in range(min(len(A[0]),len(A[1]))):            newRow = [A[0][i] , A[1][i], row[2]]            insertC.insertRow(newRow)     del insertC      this produces the correct amount of new rows based upon the content of the original concatenated rows. but when doing this, i only get the 2 fields TYPE and AREA next to a new FID. each entry also has a uniqueID that has to come across as well, so i tried (amateurishly) this;               import arcpy     searchC =  arcpy.da.SearchCursor("FeatureClass", ("UNIQUEID", "TYPE", "AREA", "SHAPE@"))          insertC = arcpy.da.InsertCursor("newFeatureClass", ("UNIQUEID", "TYPE","AREA", "SHAPE@"))           for row in searchC:         A = [x.split(",") for x in row[:-1] ]              for i in range(min(len(A[0]),len(A[1]),len(A[2]))):             newRow = [A[0][i] , A[1][i], A[2][i], row[3]]             insertC.insertRow(newRow)     del insertC      However, this now leaves the original concatenated row plus one new row that takes only the first entry from TYpe and AREA. UNIQUEID does populate both rows. Any help? i want to split out the concatenated entries but make sure the UNIQUEID that applies to all is carried across, thanks very much


title: How to write sderaster commands

text: I am trying to receive information about a Raster image in an sde Database. I am using the sderaster -o list command which gives me the error "ERROR: No [opotions]". I have tried to enter all the possible combinations of raster name, but I can´t get further with it. I am guessing that table refers to the dataset name. but what about the connection paramaters or database name? Does sderaster only work with sde services, or also with direct connect? Version: Arcsde 9.3


title: Viewing .las format LIDAR data in QGIS?

text: Is there an easy path to visualizing this LIDAR data? I have QGIS, and some USGS LIDAR data in (.las) format downloaded from http://lidar.cr.usgs.gov/ . This means I have both the (.las) and metadata in (.xml) format. I am aware liblas, but not how to apply it to this task. I am running on Ubuntu 11.04 with QGIS 1.7.0-Wroclaw. Presumably GRASS is available to me as well. A similar question for ArcGIS is How to convert LAS to GRID/Raster for ArcGIS input? . I just need to get a sense of this data and the registration... conversion to a DEM would be OK if I can visualize it.


title: How to import LAS to QGIS?

text: I need to import LAS files to Qgis. I´ve never worked with Qgis till now. Could anybody explain me this problem to the details? I´ve installed Quantum GIS v. 1.8.0. Thanks


Is it possible to load an .las file directly into QGIS? I would like to be able to load LiDAR data and to be able to do analysis based various feature attributes; e.g., elevation, classifications, etc. While LASTools (which was suggested in What is the procedure to load LAS files in QGIS 2.0.1?) is a nice set of tools, it still didn't produce an .las point file that could be read by QGIS


title: Display errors with QGIS Python Console

text: I'm using qgis 2.0.1 and I'm just starting to try to use the python console. From what I can tell the functionality of the console is fine (I'm not getting any unexpected errors) but when I'm typing in the command line the entire line goes blank with every other keystroke. This makes it pretty frustrating to edit my code. Has anyone else experienced this issue and found a solution for it? Also, since I'm just starting with pyqgis are there any IDE's that people find easier to use than the qgis python console? Thanks!


title: Query to create polygone from lines that touches in POSTGIS

text: I am trying to create polygons from a line layer with a spatial query in postgis ![lines](http://i.stack.imgur.com/BdOyr.png) the probleme is i dont have a lot of information on the lines only id of each line and geometry thanks for your help


title: osm2pgsql: modifing the default.style - newly added columns are empty

text: I would like to import a osm-db (from planet) with the "name_de" and "source" column. So I added/changed those lines in the default.style:               node,way   source       text         linear         node,way   name_de      text         linear      After import, both columns are created but empty. What I am doing wrong?


title: adding custom tools to a toolbar

text: I have created 3 custom tools, using vb.net, to run in ArcGIS. They all runs fine.They are inside a same project, in three different tool classes . Now i want them all to be added in a single toolbar, so that it will be easier to access them. To do so, i am using the basetoolbar Item, inside Extending ArcObjects here is little bit of whats inside the class               Public Sub New()          '     'TODO: Define your toolbar here by adding items     '      AddItem("CustomTool_Openlayer")      AddItem("CustomTool_selectfeature")      AddItem("CustomTool_splitline")               'BeginGroup() 'Separator     'AddItem("{FBF8C3FB-0480-11D2-8D21-080009EE4E51}", 1) 'undo command     'AddItem(New Guid("FBF8C3FB-0480-11D2-8D21-080009EE4E51"), 2) 'redo command     End Sub          Public Overrides ReadOnly Property Caption() As String         Get             'TODO: Replace bar caption             Return "My VB.Net Toolbar"         End Get     End Property      after i run the program, i get the MY VB.net Tool bar in the arcmap, but it looks something like this, and none of the tools work. ![enter image description here](http://i.stack.imgur.com/fL3e8.png) this is the small description of the tool               MyBase.m_category = "mytools" 'Localizable text.             MyBase.m_caption = "Open layer And Select Feature" 'Localizable text.             MyBase.m_message = "Open layer And Select Feature" 'Localizable text.             MyBase.m_toolTip = "Open layer And Select Feature" 'Localizable text.             MyBase.m_name = "CustomTool_Openlayer"       someone tell , what to write in the `additem()` function, i aslo tried with the GUIDs of the tools               Public Const ClassId As String = "4f76ddec-f1b4-472e-90e2-97257c9c23a2"         Public Const InterfaceId As String = "abdd6064-2ed1-47d4-8af5-fe4836515927"         Public Const EventsId As String = "ac56be15-367c-4082-966f-129871b28438"        1. should i use this basetoolbar item, or is there some other way. i just read that in basetoolbar commands are added only not tools. ?? please correct me.


title: How many levels for administration divide?

text: Excuse my ignorance, I just want to know how many levels for boundaries divide in a global level? What's Admin 4 btw?


title: OSM-based military map

text: Is there an openstreetmap based system that makes it possible to draw something like military symbols? It should be possible to draw frontline, zones of control, battle events, troops, vehicles and so on. I think, there should be such project, because it would useful for analyzing modern conflicts or even planning historic reconstruction events.