title: Geoprocessing ArcGIS 10: Sending lat/long/elevation data to a geoprocessing service

text: I have a silverlight application that allows the user to select points to make a polyline. The polyline is densified per the users input (say every two feet). The new points are sent to an esri elevation service to get the elevations at those points. I have the lat/long/elevation in a c# dictionary with a incremented number as the key. What is the best way to send this over to a geoprocessing service and what should I have my gp service take? I want to send something lightweight like json but I haven't fused with json much (except for trying to deserialize it with c# and failing so I just parsed the return string myself using c#). If I call the gp service (similar to how my code is calling the densify and elevation services, what does the gp service need to take? Are there alternatives? What we will be doing is finding the average percent slope across the polyine. Once the average percent slope is found, we need to generate a polyline shapefile or feature class with the associated percent slope attributed and we'll send this over to autocad folks. EDIT: Here's my code..!   //This is defined as a constant string at the top of the Densify class Densify.SUMBITGEOMETRY = "http://gis.mustangeng.com/ArcGIS/rest/services/Geoprocessing/gpServiceTest/GPServer/gpService/submitJob?Input=";               string dog = string.Join(",", colorList.ToArray());     //String.Join(",", dictElevat.Keys.Select(o => o.ToString()).ToArray());     //HtmlPage.Window.Alert(string.Format("Count: {0}", dog));     HtmlPage.Window.Alert(String.Join(",", dictElevat.Keys.Select(o => o.ToString()).ToArray()));     //UriBuilder builder1 = new UriBuilder(Densify.SUBMITGEOMETRY);     //builder1.Query = string.Format(String.Join(",", dictElevat.Keys.Select(o => o.ToString()).ToArray()));     WebClient webClient1 = new WebClient();     webClient1.Headers["Content-Type"] = "text/plain";     Uri serviceUri = new Uri(Densify.SUBMITGEOMETRY);     string s = String.Join(",", dictElevat.Keys.Select(o => o.ToString()).ToArray());     webClient1.UploadStringAsync(serviceUri, "POST", s);                                 HtmlPage.Window.Alert(serviceUri.ToString());


title: How to generate a random but believable digital elevation model?

text: Is there a way to generate a DEM, programatically or otherwise (arcpy or R?), that will feed into ArcGIS for further spatial analysis? Perhaps this needs to be broken down into smaller incremental steps...   1. Generate a grid   2. Fill grid with values where: `0 > value < maxElevation`   3. Neighbouring cells: `(x1-x2) < maxSlope` Strange request? Anyone have any ideas?


title: Define a projection for one point

text: I'm new to ArcGIS and am using ArcMap 10.1. I wanted to find the distance from a point (lat 36N long 12E) to all the polygons in my shape file. People on this forum helped me out and said that I should create an "XY layer" with only that point, and then use the 'spatial join' tool to find the distance from the polygons to the point. When I use 'spatial join' it says that the "XY layer" I created doesn't have a defined projection. So I tried to define the projection, but no matter which Coordinate System I try, it fails to do so. (My polygon shape file is defined under the GCS_WGS_1984 coordinate system). Thanks for your help!


title: Any reason NOT to convert to a single projection when loading GIS data?

text: I have been given the task of:   * Loading a list of properties with lat/long coordinates into a SQL Server 2008 database.   * Loading points and areas (polygons) of interest from a variety of shapefiles from various local, state, and national sources.   * Allowing the user to query the areas of interest within certain distances of a given property.   * Allow another programmer to extract both the property and points/areas of interest to KML files for visualization in Google Earth. My problem is that each source of shapefiles uses their own projection and SQL Server cannot perform comparisons on geographical data using different projections, nor can it convert from one projection to another. I have determined that, using the pyproj library, I am able to convert a point from a given projection to lat / long. I tested this by converting a point location known to be in Pinellas Park, Florida and then displaying the corresponding lat / long in Google maps — and the displayed point is in Pinellas Park (although Google's estimate of the building number is off a little; 5551 XYZ Ave instead of 5541). My question (at last): is it a safe and reasonable strategy to convert all the points in the incoming shapefiles to lat / long and store them that way, using the SQL Server standard SRID of 4326? Is there any danger of loss of precision or shifting in the conversion? If the converted lat / long polygons are displayed in Google Earth will they be identical to a display produced from the original shape files?


title: Error while creating spatial database connection

text: I'm using ArcGIS 10.0, and I've already created my database using "ArcSDE for PostgreSQL Post installation" .. My database is successfully created in PostgreSQL and now that I'm trying to create a spatial database connection from ArcMap, so I went to : ArcCatalog ==> Database connections ==> Add spatial Database connection.. I've filled all the textboxes with the needed information and when I test the connexion before I submit, I get this error : > Failed to connect to the specified server. > > Entry for SDE instance not found in services file. Although I think I've provided the right information as below : > Server : localhost > > Service : sde:postgres:localhost > > Database : my_project


title: Adding controls to OpenLayers maps - is this good code

text: I have multiple OpenLayers maps, each map gets a subset of controls that I use across all maps. In the past I have added controls in the traditional way, for each control.               var control = new OpenLayers.Control.xxx({ options: whatever });     MAP.addControl(control);      I am considering defining an object that would contain most/all of the controls that I might use. That object would be is a .js file would be loaded into each web page. I would then add the controls I need for that page. the control object, that is in the control.js file               var MY_controls = {         scale:  new OpenLayers.Control.Scale(),         mousePosition :new OpenLayers.Control.MousePosition({numDigits:6}),         panZoomBar: new OpenLayers.Control.PanZoomBar(),         ctrl_zoomin: new OpenLayers.Control.ZoomBox({ disable:true }),         ctrl_zoomout: new OpenLayers.Control.ZoomBox({out:true,disable:true})     };      then on each page I would do something like this               array = ['scale','mousePosition','panZoomBar'];     for ( var x=0; x < array.length; ++x ) {         MAP.addControl(MY_controls[array[x]]);       };      Is this good code?


title: How far ESRI Technical Support goes?

text: Given the number of questions about basic/intermediary usage of ESRI products on this site and in a lot of other mail lists, I became curious to know why people don't ask them to the technical support. So, for registered ESRI users:   How the ESRI technical support works and what kind of questions does it answer?


title: How to convert a bing maps quadkey to mercator coordinates

text: How do you convert a bing maps quad key, such as "03200320023" to coordinates in mercator virtual earth (epsg 3857) coordinate system? I've got an algorithm to convert to coordinates in latitude longitude (e.g. the quad key above comes to "-84.19921875,33.7243385314941,-84.0234375,33.8704147338867"). I could then convert from latitude longitude to mercator coordinate, but that smells far too expensive in processing power. The mercator virtual earth result I'm looking for given the above quadkey is: "-9373014,4011415,-9353446,4030983" Also a question on stackoverflow


title: How to avoid dbo.arcsde LOCK issue while exporting to ArcSDE

text: I had a spatial db connection open in arcmap. While exporting to SDE **with a console app** I get an error in the console: > Lock request conflicts with an established lock [DBO]. Does this just mean someone is using the same ArcSDE connection or is it the feature dataset? Of course when I closed arcmap, and fired the console app again, my feature class exported to ArcSDE successfully. Can I kill any locks programmatically if I am DBO (when other users are using the same FD)? If that's not smart, how can I avoid this lock? How does the GIS community deal with this issue, wait till all users are off and run it at 2am? Btw, the console app will become a batch process updating FCs on a regular basis. I appreciate your feedback, thanks! saim


title: Paris multimodal network dataset not giving all forms of transport available

text: After having completed **Exercise 2 in the official Esri Network Analyst Tutorial** , I decided to find a route in the multimodal dataset, to see how it functions. After adding several network locations along an intentionally long route, I pressed **Solve**. I was quite surprised to find that the route chosen took walking along the streets as the only transportation mode available (despite railway stations and lines being included in the dataset). I then set the Impedance to TravelTime also, but still it is not showing the result. Can you give me suggestion, please?


title: Identify functionality using mapserver

text: I am just creating a simple web application using mapserver. It has the basic pan, zoom in/out functionality. Now what I want to try is to get an identify functional added into it so that when I click on a feature it populates the values in a table or a pop up. Apart from using openlayers could anyone tell me how can this be implemented. Also the data is coming from PostgreSQL. Thank You.


title: Problem loading table in ArcMap

text: I am an ArcGis 10. ( ArcEditor licence) user and I have a problem when trying to load data with the Load Objects Tool in ArcMap. I want to load a table, but after I choose the input data (equipment.dbf) - Add- and press next , ArcMap instantly shuts down and the well known window appears "ArcGis Desktop has encountered a serious application error and is unable to continue. One more thing, I am trying to load the table into a file geodatabase with a geometric network and relationship classes (but that shouldn't be a problem as support says). Does anyone know why? and what is the solution? Thank you, Ligia


title: What components does OpenLayers.Strategy.BBOX require for compilation?

text: I recompiled OpenLayers with minimal settings, and now Vector layer, that uses BBOX strategy, stopped adding bbox parameter to the queries. What components do I need to add? This is the current config:               [first]     [last]     [include]     OpenLayers/Map.js     OpenLayers/Kinetic.js     OpenLayers/Projection.js     OpenLayers/Layer/Vector.js     OpenLayers/Layer/OSM.js     OpenLayers/Layer/Bing.js     OpenLayers/Layer/WMS.js     OpenLayers/Layer/Google/v3.js     OpenLayers/Popup/FramedCloud.js     OpenLayers/Control/Navigation.js     OpenLayers/Control/Zoom.js     OpenLayers/Control/Attribution.js     OpenLayers/Control/SelectFeature.js     OpenLayers/Control/Panel.js     OpenLayers/Control/LayerSwitcher.js     OpenLayers/Renderer/SVG.js     OpenLayers/Renderer/VML.js     OpenLayers/Format/GeoJSON.js     OpenLayers/Protocol/HTTP.js     OpenLayers/Strategy/Fixed.js     OpenLayers/Strategy/BBOX.js     OpenLayers/StyleMap.js     OpenLayers/Rule.js     OpenLayers/Filter/Comparison.js     OpenLayers/Filter/Logical.js     [exclude]      This is quite important, because the complete OL is 990K, while the minimal is only 350K big, and my users may be in remote areas with very slow connections. **It was not just Filter/Spatial**. I noticed that even with Spatial it didn't compose a query string. I bisected the whole list of components and it was... **OpenLayers/Format/QueryStringFilter.js**.


title: Able to circumvent QgsCredentials script called from Pre-execution

text: I've got a stored project that, when opened is connecting to my PostgreSQL DB, and is loading several layers. All good. I'm using a script that is called from Pre-execution (Processing/Options and configuration) that is using the QgsCredentialDialog to distinguish between admin and regular user.               def LoginCredentials():       uri = QgsDataSourceURI()       uri.setConnection("localhost", "5432", "Blah", None, None)       connInfo = uri.connectionInfo()            (success, user, passwd) = QgsCredentialDialog.instance().get(connInfo, None, Nones)       if success:         uri.setPassword(passwd)         uri.setUsername(user)         iface.messageBar().pushMessage("Welcome!", "Welcome "+str(user), level=QgsMessageBar.INFO, duration=3)          def run_script(iface):       LoginCredentials()      Issues:   * one of my 7 layers is able to be loaded before the credentials dialog pops up.    * the credentials dialog is printing a string just above the username/password textboxes with dbname, host, port, user, password and sslmode, all filled out! (I'm guessing the values are read from the .qgs project file, but I rather it didn't show it)    * if I keep hitting the Cancel button I can see that the layers are loaded one by one until everything is loaded and the credentials dialog just disappear.  I'm _very_ new to Q-GIS (4 days), so it might be that I've overlooked something very essential. Cheers!


title: How to keep Table Order in Identify Results?

text: Is there a way to keep the default table order when using the "Identify Feature" tool. I see that you can sort the "identify results" by table column name or value. I don't want the data in this view sorted at all. In our database, we've already standardized our column order so it's easier for the users. Thanks in advance!


title: How to convert int type OBJECTID to ObjectID type recognised by ArcSDE?

text: I have an ordinary table with OBJECTID column (type of int) in a MS SQL Server 2008 R2 instance. I need to store it in an ArcSDE geodatabase instance running on top of the SQL Server. The thing is how to make the OBJECTID field to the esri-compliant Object ID type so that feature operations through Feature Server can be applied down the track?


title: How to import qgis.core to python app?

text: I'm new to Python. I want to create GIS app with Python using QGIS library on Window 7 64bit. The way I do: install QGIS 1.7.0   * I set PYTHONPATH=D:\Program Files (x86)\Quantum GIS Wroclaw\apps\qgis\python   * I run python from cmd in D:\Program Files (x86)\Quantum GIS Wroclaw\bin   * I type import qgis.core, a pop up window appear with text > The program can't start because qgis_core.dll is missing from your computer. > Try reinstalling the program to fix this problem and text in cmd is > ImportError: DLL load failed: The specified module could not be found. Can anyone help me. Thank in advance.


title: What causes ERROR 000732 from MosaicToNewRaster?

text: I'm trying to run the MosaicToNewRaster tool as part of a python script. This is my code:               input_rasters = arcpy.GetParameterAsText(5), arcpy.GetParameterAsText(6), arcpy.GetParameterAsText(7), arcpy.GetParameterAsText(8)     output_location = "C:/Birkbeck_Python/Prog_Assign_1.gdb"     arcpy.MosaicToNewRaster_management(input_rasters, output_location, "R_Dataset", "", "",    "", "1", "", "")      I get the following error message, anybody got any ideas? It works fine when I hardcode it. Thanks ExecuteError: Failed to execute. Parameters are not valid. ERROR 000732: Input Rasters: Dataset '';'';'';'' does not exist or is not supported Failed to execute (MosaicToNewRaster).


title: Help loading Raster to PostGIS 2 .0

text: Similar questions have been asked, but since my reputation is not enough I am unable to comment on answers. (The following is from an answer by nickves for this question Loading raster data to postgis) **Example:** raster2sql -s -t 50x50 -f -I -Y myRaster.tif rasterSchema.RasterName > raster.sql` psql -h localhost -U someuser -d mydb -f raster.sql **Question1:** What does the following mean in relation to the above example:   1. rasterSchema.RasterName   2. raster.sql`psql   3. -f raster.sql Question2: Given the following information, what would the syntax look like for loading a raster to a PostGIS 2.0 database?   1. Image location = C:\Users\Desktop\dem.tif   2. SRID = 2048   3. Database name = Rasters   4. Schema = Elevation My os is windows 7 (not sure if this makes a differnce). Regards


AS a crime analyst you need to establish a few maps that would help police to do targeted patrols, I've given data of suburbs, roads, incidents, boundaries. How can I establish patterns, hotspots and clusters?


title: Simulating Random Agent Movements along Roads on Map

text: I need to simulate agent movements on map as follows: Agent starts a known location and moves along a road. At each road intersection, the agent uniformly chooses one possible road and moves along that. The path chosen is to be traced with a poly-line. I know how to choose a path uniformly and mark it with poly-lines once I have the data. But I am looking for how I can retrieve data like road intersection points and the roads leading from there on.


title: For a non-geocentric ellipsoid, how does the Lat/Long/Ell Ht & semi-major axis & flattening implicitly define the centre of the ellipsoid?

text: For example if we compare the AGD84 (Australian Geodetic Datum) with GDA94 (Geocentric Datum Australia), the centre of the AGD Datum is 200 metres different to GDA, because it is a best fit to the geoid. But how is this difference realised (implicitly) given only the Lat/Long/ Ellipsoidal Height of a Point (Johnstone), the semi-major axis (a) and the flattening (f)? I realise that (a) is different by 23 metres, but it is the angle between the 2 datums that makes the most difference, this is what I am trying to get a handle on. Asking another way, how could I calc the ECEF coords for the centre of the AGD, to see how different it is to GDA? Asking yet another way, how to determine the shift & angles between the axes of of AGD & GDA? I already have the GDA manual - http://www.icsm.gov.au/gda/gdatm/gdav2.3.pdf and am working through writing C++ code for this. I also understand how AGD is not very consistent, so a conversion between the 2 datums requires the use of a model, it is just that I am trying to truly understand the mathematical difference between the two. Sorry for the particularly nerdy question, & thanks in advance for any answers.


title: Storm Water flooding

text: Is there a method in Arc Gis for mapping storm water runoff in an urban environment assuming that the sewers are at full capacity. I just want to map roughly where the water will flow over the land over paved surfaces I have a DEM, a land use map and have already calculated the flow direction and streams network. I wish to find out what areas will flood in the study area during an extreme rainfall event.


title: How to request group of layers from Geoserver with turn off/on capabilities? Geoserver layer organizer

text: I would like to organize my layers by topics, and for each of them (if chosen on the client-side), I would like to request a single-tile image, which contains the predefined layers in the correct rendering order (e.g. if one topic has 5 layers, it renders them from the bottom to top like this: water, roads, highways, railways, poi). But, I also want to be capable to request just four of it, or three, etc. in any combination but still in the correct order. My first idea was Geoserver's "Layer group", but as I've learned that is really more suited for group of layers which are not meant to be turned on/off seperately. (Please let me know if I am wrong about this, though.) And I also kinda lose Getfeatureinfo capability with this. After some research I figured out, that in the GetMap request, if I add all my layers' name in the correct order (from bottom to top-->LAYERS=water,roads,highways,railways,poi), I get the same result, and I can also narrow down the list to the required ones. The GetMap solution seems to by right, but in that case I need to fill up the list for each topic on the client side. So my main question: **is there a solution to get a list of layers from Geoserver by topic in the right rendering order**? Or is there an extension for organizing layers this way? (I haven't seen any extension for that purpose) My idea was that if I organize all layers of a topic to a workspace, I can request a getCapabilities for the workspace WMS endpoint, parse and read out the layers, but unfortunately it gives back the layers in a pseudo-random order, so not really good. Maybe a custom getCapabilities XML can solve this? Please share your ideas, what's the best practice to do this? I'm on Geoserver 2.4.2.


title: Create a web application with ArcGIS Server

text: I have a question about ArcGIS server: I want to make an application for users who can search for proper houses, like a house-finding engine, and I published my map to ArcGIS Server but how can I link my map to the ArcGIS JavaScript API and how can I write queries for my users? For example: in my web application, I want the user to choose and select from drop boxes like price, region street or rooms?


title: Matching items from two tables

text: Can I select items in one table by matching items from a duplicate column in another table in QGIS?


title: Databinding ComboBox to GraphicsCollection - ESRI, WPF runtime

text: I'd like to bind a GraphicsCollection (which is a property in a view model) to a combo box in a XAML window. The problem is I can't figure out how to display attribute values in the combo box. I can get it to display the Geometry property like this:               <ComboBox DataContext="{StaticResource MyViewModel}"                     DisplayMemberPath="Geometry"                     ItemsSource="{Binding CollectionOfGraphics}"  />      But I want to be able to do something like this that will display the OBJECTID or some other attribute:               <ComboBox DataContext="{StaticResource MyViewModel}"                     DisplayMemberPath="Attributes["OBJECTID"]"                     ItemsSource="{Binding CollectionOfGraphics}"  />      I also tried this approach - unsuccessfully:               <ComboBox DataContext="{StaticResource MyViewModel}"                     DisplayMemberPath="Attributes[0]"                     ItemsSource="{Binding CollectionOfGraphics}"  />      Platform: ArcGIS 10.1, WPF Runtime 10.1.1 **Solution based on Antti Kajanus's answer below:**               <ComboBox DataContext="{StaticResource MyViewModel}"                     ItemsSource="{Binding CollectionOfGraphics}" >          <ComboBox.ItemTemplate>               <DataTemplate>                   <TextBlock Text="{Binding Attributes[OBJECTID]}"/>               </DataTemplate>          </ComboBox.ItemTemplate>     </ComboBox>      Note that the attributite name (OBJECTID) did not have to be quoted.


title: How do I use use custom fonts for labelling in Geoserver?

text: I want to use a custom font I have downloaded from the web with Geoserver for my labels. How do I specify the location of my font so that the font-family css parameter can read it? Ando Geoserver 2.1.4


title: What's the difference between hypsography and topography?

text: I've recently come across some data that has been labeled "hypsography". I looked it up and it sounds like topography. But this data source has data labeled "topography" separate from its "hypsography". What's the difference? I've read the wiki for 'Elevation' and its starting to sound like all the "topos" I've been reading are really "hypsos"


title: Cannot get gdalwarp to work

text: I am having some trouble getting gdalwarp to work. I want to clip my input raster `in.img` to the boundary of a polygon stored in a shapefile `polygon.shp` and save it out as `out.img`. The data type of `in.img` is 32-bit float. The raster cells external to the polygon delineation have a value of 0.0. I figure gdalwarp ought to be able to clip `in.img` to my polygon boundary, effectively setting the external 0.0 values to NoData (that is what I _want_ to do, anyway). I am using the following command at the Windows command prompt:               gdalwarp -of HFA -cutline polygon.shp in.img out.img      The problem is, this doesn't seem to achieve anything. `out.img` is essentially a copy of `in.img` in that the cells external to the polygon boundary are still 0.0—my understanding was that this would only happen if I used the `-crop_to_cutline` flag. What am I doing wrong?


title: What is required to get OGR to work with SQL Server?

text: All, I've been having heaps of trouble getting OGR to connect to our SQL Server development environment. (read about it here and here) **SQL Versions** : 2008, 2012   **OS's** : Windows 7, Ubuntu 12.04 My questions are: 1) What permissions does the user need in order to establish a connection to a SQL Server database with OGR? 2) Does anything need to be loaded on the SQL Server machine or application server itself? I used to be able to connect even with QGIS, but that isn't working now either. (or here either) Also: if you have a password with a ! in it, it will throw back a 'bash' error...fyi! Ok thanks again, all!


title: Why OSM data has different attribute values from the same element and same time frame?

text: Here is my question: Does anyone know why OSM has different attribute values from the same element and same frame time? The example explains itself, when I try the Shapefile from different mirrors like **Planet.osm.org** , **bbbike.org** and **geofabrik.de** ; I find that they have a mistake with the Direction of Travel on its attribute value. When I check in the website: here the link , I see on its attribute value that there is no mistake at all and has its right value. ![Example](http://i.stack.imgur.com/3dEUJ.png) So, I thought that **those mirrors extract the exact information** from osm servers, but I think that there is a lack of info. Is that right? Thanks for your time,


title: District heating expansion design with ArcGIS

text: I am quite new at GIS but I have some experience of the basic concepts involved in a workflow of them. What I currently want to do is to calculate the cost of expanding an existing district heating network to the houses of a village. More specifically: I have a village of 310 buildings. In the village there is already a network that some houses are connected to. So, basically I have 310 points in one shapefile and another shapefile that represents the existing district heating network. What I want to achieve is that I want to expand the existing network so that all nodes are connected with the minimum distance to it. That can mean that they can either connect to the existing part or new pipes have to be drawn towards areas and houses can connect to those pipes via sub-pipes (i.e. smaller pipes). I know that overall this is not a very realistic representation of such a design since no constraints are taken into account, but it would be a good starting point for my analysis. Thanks a lot! edit 1: I am capable of using both ArcGIS and QGIS and I do have some programming experience with both Matlab and Python. edit 2: So far I have looked at Geometric Networks and the distance tools of the Spatial Analyst toolbox for ArcGIS. Could either of those be the solution?


title: Openlayers tiles source with poor resolution, how zoom out?

text: I am using Openlayers to display a map from Vienna Open Data (Capabilities) in an app for mobiles, which have higher DPI screens. The resolution of the tiles is not so good (I suppose because they are JPEG instead of PNG). Is it possible to use the tiles from the upper zoom level as if they where smaller tiles for the current level, faking the resolution?


First of all, the code which is below gives me the attribute names of a shapefile correctly if I use PyScripter IDE:               from qgis.core import *     from qgis.gui import *     from PyQt4.QtCore import *     from PyQt4.QtGui import *          QgsApplication.setPrefixPath("C:/OSGeo4W/apps/qgis-dev", True)     QgsApplication.initQgis()           layer = QgsVectorLayer("C:/istanbul/istanbul_geo-wgs84.shp", "sero", "ogr")          fields = layer.pendingFields()          for field in fields.itervalues():         print "Name:%s Type:%s, Length:%s" % ( field.name(), field.typeName(), field.length())      If I want to access attribute names of a shapefile, Aptana Studio gives me an error message below:                   for field in fields.itervalues():     AttributeError: 'QgsFields' object has no attribute 'itervalues'     d:\src\qgis\src\core\qgsvectorlayer.cpp(188) : (QgsVectorLayer::~QgsVectorLayer) entered.      According to the console in Aptana Studio, my system can read vector layers:               d:\src\qgis\src\core\qgsvectorlayer.cpp(1628) : (QgsVectorLayer::setDataProvider)      Extent of layer: 28.0019662445739925,40.7878698688807120 : 29.9327258493754194,41.5270420672795026          d:\src\qgis\src\core\qgscoordinatereferencesystem.cpp(449) : (QgsCoordinateReferenceSystem::createFromWkt) authid recognized as EPSG:4326      I used OSGeo4W installer and configured my IDE based on Nathan's expressions: How to add QGIS Python API Intellisense support in Aptana Studio3/Pydev and finally my "aptana.bat" file to initialize the IDE:               @ECHO OFF     set OSGEO4W_ROOT=C:\OSGeo4W     call "%OSGEO4W_ROOT%"\bin\o4w_env.bat     call "%OSGEO4W_ROOT%"\apps\grass\grass-6.4.3\etc\env.bat     set PATH=%OSGEO4W_ROOT%\bin;%OSGEO4W_ROOT%\apps\qgis-dev\bin;%PATH%     set PYTHONHOME=%OSGEO4W_ROOT%\apps\Python27     set PYTHONPATH=%OSGEO4W_ROOT%\apps\qgis-dev\python     set QGISHOME=%OSGEO4W_ROOT%\apps\qgis-dev     call "C:\Users\Geomatics\AppData\Local\Aptana Studio 3\studio3.bat" %*      How are we gonna tackle this problem? Thanks in advance...


title: r.watershed analysis error from raster map being in different zone and unable to be opened

text: I tried running my DEM (geotiff) file through the `r.watershed` program by the GRASS GIS plugin of QGIS. I believe the error message results from a raster map in the different zone than the current region but I set my GRASS GIS mapset in WGS84-Zone 16N coordinate system where it ought to be. Thank you in advance. Here is the error message I get:               r.watershed elevation=dem2@Tejeras threshold=30 accumulation=dem2_accum drainage=dem2_drain stream=dem2_stream basin=dem2_basin          SECTION 1a (of 5): Initiating Memory.          ￼Raster map is in different zone (16) than current region (0)          ￼Unable to open raster map           ￼unable to open elevation map layer          ￼Subprocess failed with exit code 1          ￼category information for [dem2_accum] in [Tejeras] missing or invalid          ￼category information for [dem2_drain] in [Tejeras] missing or invalid          ￼category information for [dem2_basin] in [Tejeras] missing or invalid          ￼category information for [dem2_stream] in [Tejeras] missing or invalid     Finished with error


title: Spatial Aggregation of attributes from line structures in polygons

text: maybe someone has an idea how to deal with it (or directly knows the equivalent functions in PostGIS) for the following problem: I have a transport network (something like railways) and the attributes of this lines are mostly the transport goods, their corresponding weights and the transportatino dates. In a big railway station lots of these lines come together (all lines are connected and "covered" with a polygon "railway_station" which has not attributes). Now I would like to aggregate the weights of some goods in a defined period based on the polygons area. Do any function like "Spatial Aggregation" or "Spatial Join" exist which is able to give me the expected results? Hopefully someone has an idea! :) Best regards, Jann


title: What SRID is this projection?

text: Anyone know what the SRID is for the following projected coordinate reference system:               {        PROJCS   [           "Florida 0901, Eastern Zone [1983, US Survey feet]",           GEOGCS      [              "NAD 83 [Continental US]",              DATUM         [                 "NAD 83 [Continental US]",                 SPHEROID            [                    "GRS 80",                    6378137.0,                    298.257222101                 ]              ],              PRIMEM         [                 "Greenwich",                 0.0              ],              UNIT         [                 "Decimal Degree",                 0.0174532925199433              ]           ],           PROJECTION      [              "Transverse_Mercator"           ],           PARAMETER      [              "False_Easting",              656166.666667           ],           PARAMETER      [              "Central_Meridian",              -81.0           ],           PARAMETER      [              "Scale_Factor",              0.9999411764705882           ],           PARAMETER      [              "Latitude_Of_Origin",              24.33333333333333           ],           UNIT      [              "U.S. Foot",              0.304800609601219           ]        ]     }      My guess is: 102658. I'm assuming I can just make sure all the different values are the same and I'm good. Thank you for your time.


title: how to configure an arcgis server service to enable multiple editing

text: I'm using ArcGIS Server 9.3.1, ArcSDE, VS.NET 2008 and I'm trying to enable a service to be edited multiple times, the service its going to be public and it will be used to report emergencies and other situations that will require attention. I need some pointers on how to accomplish this. Is it a good idea to have an editable map or its better to save the points coordinates to a database and in a report like app generate all the points at once?


title: How do you make 3d features from survey point data in arcGIS 9.3?

text: I have survey point data brought into ArcGIS using Survey Analyst Extension. I have created my 2D shapefiles of the cross-sections, thalwegs, and roadways but I can not get it to populate the zfield from the point data.


title: Automated Importing of Landsat with new USGS MLT metadata file

text: Has anyone had any luck importing Landsat with the updated MLT files? Besides ENVI 5 SP1 I see no software that can create images from these new files. I can do a manual layer stack but then the band data is not embedded in the header and so I can't pan sharpen and manually updating the headers is a pain when I'm trying to import several dozen files. Thanks, Steve


title: How to determine the most accurate elevation algorithm

text: I use an Android GPS cycling app to track my bike rides and I would like to get a reasonably accurate estimate of the total ascent I did on a ride. I understand that GPS is very inaccurate at determining elevation, so I have been importing the resulting GPX files into several sites that use mapping data to determine elevation, and I've found they're ridiculously inconsistent. For example, from a 30-mile ride I get the following climb data from four sites using the same GPX file:               CardioTrainer   980     Strava          1085     MapMyRide       554     RideWithGPS     2131      The results are consistent over multiple rides, and I can add more sites with an equally mixed set of results. I think I can easily toss RideWithGPS as a badly flawed outlier based on subjective experience alone, but MapMyRide is consistently about half the ascent of the others, so it makes me wonder who is off by 2-to-1. My question is, how can I determine whose algorithm is the most accurate without spending hundreds on a device with a barometric altimeter?


title: Moving window regression (MWR) in R

text: I'm looking for some package in R to estimate the value of a variable by moving window regression. I have read in _Local Models for Spatial Analysis_ , which is a good method to predict values at unsampled locations. A very similar case is mentioned in the book,PRISM. I've been reading various R packages, but have not found anything about MWR, GWR only. Anyone know if there is another way of naming or GWR can be adapted to function as MWR? Thanks Text from book -> Local Models for Spatial Analysis


title: Convert GeoJSON form EPSG 3857 to 4326 with JavaScript

text: Given a GeoJSON in EPSG3857:               {       "type": "Feature",       "properties": {},       "geometry": {         "type": "LineString",         "coordinates": [           [             -13720479.997764934,             5688254.733325758           ],           [             -13720479.47473932,             5688241.794717406           ]         ]       }     }      I am trying to convert the GeoJSON to EPSG 4326 for a Leaflet map. The user input in the map comes in EPSG 3857. I want to do the conversion with JavaScript after the user enters the GeoJSON. So I am looking for a JavaScript library that I can use. I am not looking for tools with whom I can convert the GeoJSON outside of my script (like GDAL/OGR) . I tried Proj4 but didn't get it to work. Is it possible to use Proj4 to convert line strings? If not, what other options should I consider?               var epsg4326 = new Proj4js.Proj('EPSG:4326');     var epsg3857 = new Proj4js.Proj('EPSG:3857');          var geoJSONOBJ4 = new Proj4js.Point(geoJSONOBJ3);     Proj4js.transform(epsg4326,epsg3857,geoJSONOBJ4);


title: in QGIS, while using simplify feature, slider does not appear

text: Using Qgis both 1.8.0 and 1.7.4, the slider does not show up after selecting features and click "simplify feature" botton. Any special setting needed to show the slider? Thanks.


title: How can Flocking Algorithms be used in GIS?

text: The first flocking algorithm was written by Craig Reynolds in 1986. A recent browser implementation can be found here.                   ![flock of seagulls](http://i.stack.imgur.com/1vXwx.gif) The notes refer to " _The use of Flocks to drive a Geographic Analysis Machine_ " (1998) by James Macgill and Stan Openshaw which > uses a flock model, with communication between boids, to better search for > clusters in spacial datasets, just as a natural flock provides better > foraging than individual birds could manage. What other real world examples of this algorithm have been used in GIS, and do any GIS applications implement this functionality?


title: Problem using grid with a vector layer

text: I need to establish a model of maritime pollution and I try to use the grid option to work on ocean sectors. I tried to use my grid on a vector file, .shp, with points. Each point is a maritime pollution fact. Tried with the native Grid (vector) option in QGIS: give me an error Tried with the Qmarxan plugin (seems to be the one i need): give me an error Is it a problem with Python, my layer, my data or a bug problem? I tested Qmarxan with a Geotiff, and it worked quite well. But it seems possible to use it with a vector layer (http://datagistips.blogspot.fr/2012/04/le-carroyage- avec-qgis-et-le-plugin.html) Many thanks! The error with Qmarxan :               Traceback (most recent call last):       File "C:/Users/GF/.qgis//python/plugins\Qmarxan\qm_mkgrid.py", line 189, in accept         self.writeFile(boundBox, sideLen, hexagon)       File "C:/Users/GF/.qgis//python/plugins\Qmarxan\qm_mkgrid.py", line 333, in writeFile         self.ui.progressBar.setValue((row*100)/rowCount)     ZeroDivisionError: integer division or modulo by zero          Version de Python :     2.7.3 (default, Apr 10 2012, 23:31:26) [MSC v.1500 32 bit (Intel)]


title: Best way to delete selected mass amounts of old data from a polygon shapefile attribute table?

text: Dumb problem that for some reason is giving me trouble. I have a polygon shapefile that is loaded with old demographic data in its attribute file, I just received the new data today in xls format. What I want to do is just delete the old demo data from the existing shapefile by selecting the fields in the attribute table and right clicking delete. Then I would do a simple join using the ID numbers for each feature to add the new data from the new table then just export it out as a new shapefile. The problem is that for some reason I can not delete entire columns from the attribute table. I ask the people on this forum how would you go about stripping the data from this polygon shapefile? Thank you!


title: How To GetFeatureInfo From a GeoWebCached Layer in Geoserver

text: I'm really really struggling with this. I have created a map using Geoserver and Openlayers. Things have been working out here and there, as I am new to all this so a lot has been just trial and error. I noticed my map started loading slower and slower so I have decided to use GeoWebCache to cache my layers. The problem is that now I cannot get the feature info to display in a pop-up as it once was when I was getting the layers straight from the server/wms. The part that even confuses me further is that when I provide the pop-up with the url 'http://localhost:8080/geoserver/gwc/service/wms' it actually attempts to send a request. However when I provide the right url from the server where the data actually is (http://localhost:8080/geoserver/wms), it's like if no request is being made. (I am seeing this in firebug). I've also created a geowebcache.xml files and have listed several of my layers as true, but no differences there either :( Any suggestions you can offer me will be greatly appreciated. I am truly very stuck. Many thank yous, elshae


title: spatial database vs spatial datastructure?

text: I am building an automated system that has to process several hundred spatial queries in the form of a polygon every day. The queries have to return both spatial and nonspatial data from different database tables. Each table consists of over 500 000 records. What would be a better approach when performance is a top priority? I) execute a spatial query every time i need (spatial) data? II) first bulk-load all data in-memory into spatial datastructure, and then perform those same spatial queries onto the spatial datastructures? If the latter is the better approach, which kind of datastructure would be recommended? A few possible datastructures could be: \- PR-trees \- tries \- quadtrees \- R-trees Notes: It's safe to assume that there will be no data changes during the processing. The spatial data is 2-dimensional (points, linestrings, polygons)


title: Lock overlay on Custom Geoprocessing Tools

text: How do you get rid of the "lock" icon overlay on a custom geoprocessing tool in ArcToolbox? This seems to appear even with the ESRI Sample Projects.


title: How can I specify the CRS of a raster layer in pyqgis?

text: I manipulating netcdf raster data from the QGIS python console. The metadata of my netcdf file does not provide the Reference system. However I know it is 'EPSG:4326' I am writing a python script to perform tasks automatically, I don't want user interaction. I am using the following commands :               fileName = "/path/to/raster/file.tif"     fileInfo = QFileInfo(fileName)     baseName = fileInfo.baseName()     rlayer = QgsRasterLayer(fileName, baseName)`      But the QGIS windows asked me the Coordinate reference system of this layer. I tried to specified the crs while loading with something like               rlayer = QgsRasterLayer(filename=fileName, basename=baseName, crs='EPSG:4326')      But QGIS crashes. Do you have any idea ? Is it possible to specify a qml file where the crs is defined ? Thank you for your help


title: Use of Geoserver raster processes within OpenLayers: WCS/WPS problem.

text: I've been working with some GeoServer processes using OpenLayers as client. I've been successful in those which use WFS as input (like gs:BufferFeatureCollection). My problem comes when I want to use processes that need as input raster data (tiff). I know that OpenLayers is not a WCS client, but you can send a POST request to the server to get back the result, right?. Particularly I want to execute the gs:RasterZontalStatistics, that needs a tiff and a WFS layer as input, and returns some statistics (average, min, max, etc.) for the tiff inside the geometries of the WFS Layer. This is the code I've used so far:               /*GEOPROCESS - Calculate the average/min/max values of the MDE raster according to user drawing selection*/     function obtainElevationProcess(){         var wpsFormat= new OpenLayers.Format.WPSExecute();              var doc= wpsFormat.write({             identifier: "gs:RasterZonalStatistics",              dataInputs:[{                 identifier:'data',                 reference:{                     mimeType:"image/tiff",                     href: "http://geoserver/wcs",                     method: "POST",                     body: {                         wcs: {                             identifier:'forestry:mde200_V4',                             version:"1.1.1",                             domainSubset: {                                 boundingBox: {                                     projection: 'http://www.opengis.net/gml/srs/epsg.xml#3857',                                     bounds: new OpenLayers.Bounds(-200931.8344047, 4681799.186142525,32948.35238090271, 4948113.428142309)                                 }                                 },                             output: {format:'image/tiff'},                             }                             }                             }                             },                 {                     identifier:"zones",                     reference:{                                 mimeType:"text/xml; subtype=wfs-collection/1.0",                                 href: "http://geoserver/wfs",                                 method: 'POST',                                 body:{                                     wfs:{                                         version:'1.0.0',                                         outputFormat: "GML2",                                         featureNS: "http://localhost:8080/thesis/forestry",                                         featureType:'forestry:Lics',                                         }                                 }                             }                 }],                      responseForm:{                         rawDataOutput:{                             mimeType:"text/xml;subtype=wfs-collection/1.0",                             identifier:"statistics"                     }                 }                      });                   OpenLayers.Request.POST({             url:"http://localhost:8080/geoserver/ows/",             data:doc,             headers: {                 "Content-Type": "text/xml;charset=utf-8"             },             success: function(response){                 alert(response.responseText);                 }         });          }      Once executed, OpenLayers returns this message (selection):               Process execution ecacba6c-a676-46b5-9010-bd4a2bb35632 failed     Failed to retrieve value for input data     Unrecognized request type {service=WCS, Output=image/tiff,   Identifier=net.opengis.ows11.impl.CodeTypeImpl@1ad2fc (value: forestry:mde200_V4, codeSpace: null), DomainSubset=net.opengis.ows11.impl.BoundingBoxTypeImpl@111d0b5 (lowerCorner: [-200931.8344047, 4681799.1861425], upperCorner: [32948.352380903, 4948113.4281423], crs: http://www.opengis.net/gml/srs/epsg.xml#3857, dimensions: null), version=1.1.1}</ows:ExceptionText></ows:Exception></ows:ExceptionReport></wps:ProcessFailed></wps:Status></wps:ExecuteResponse>      Does anyone an idea about how this can be done (if possible)? My environment is: Geoserver 2.3.5 / OpenLayers 2.13 / localhost on Apache. Help will be highly appreciated. Thanks. Constan.


title: If I convert a TIN to a raster in ArcGIS will results from the two versions of the surface be comparable?

text: I have converted a TIN to a raster. I found that the top of a bund wall shifts laterally in position between the TIN and the raster into which it has been converted using the TIN to raster tool in ArcGIS. What is the reason for this? Does this type of thing always happen? Will it be sensible to do volume calcs using the raster with this type of problem or should I stay with the TIN?


title: sqlite version error with ogr2ogr

text: I am trying to randomly select some points from a shape file using ogr2ogr and an SQL query via sqlite but it seems to baulk at the sqlite version:               ogr2ogr -sql "SELECT * FROM probe_2004to2012 WHERE D2004 > 0 ORDER BY RANDOM() LIMIT 10" -dialect sqlite ./Probing/probe_2004_r10.shp ./Probing/probe_2004to2012.shp          ERROR 6: The SQLite version is to old to support the SQLite SQL dialect      Given that I tried updating sqlite3 in my /usr/bin/ from the jan 2013 version and it is now the latest (Oct 2013) version and ogr2ogr still throws up this error I suspect it may be something else that's wrong. I know I can select random points in QGIS using the "Research Tools" toolbox but I would really like to be able to script this as I have a lot of random selecting to do. -update: its seems ogr2ogr supports               -dialect sqlite3      as a valid flag and that seems to get a simplified sql query to parse but the random() selection won't get through. After the discussion below it seems that ogr2ogr is not parsing the -dialect flag "properly". The sqlite3 flag is ignored and ogr defaults to its built in sql, which is why the simple queries work. Why the sqlite flag is coming up ERROR 6 I don't know. Yet.


title: Use of Raster Calculator in ModelBuilder

text: I have the following equation:               3000*(1-exp((-0.000664)*input raster))      I want to make a tool to use this equation on raster images. I want to do this with the Raster Calculator with the help of ModelBuilder. Can you please explain the procedure for making this tool?


title: Why does my QGIS plugin MainWindow start with some widgets not drawn?

text: I have a QGIS plugin written in Python which uses a MainWindow to hold a few other widgets. When the window first gets shown, some widgets are drawn and some are not. When I resize or minimise/restore the dialog, everything is drawn correctly. What's going on? Immediately after dialog shown: ![Immediately after dialog shown](http://i.stack.imgur.com/RZQdQ.png) After minimise/restore or resize: ![After minimise/resize](http://i.stack.imgur.com/rEPIL.png)               from PyQt4 import QtGui     from PyQt4 import QtCore     from qgis.core import *     from qgis.gui import *     from ui_townplanning import Ui_TownPlanning     # from townplanning import msgbox     # create the dialog for zoom to point      townplanning.py: class TownPlanning:                   def __init__(self, iface):             # Save reference to the QGIS interface             self.iface = iface             self.canvas = self.iface.mapCanvas()             # initialize plugin directory             self.plugin_dir = os.path.dirname(__file__)             # initialize locale             locale = QtCore.QSettings().value("locale/userLocale")[0:2]             localePath = os.path.join(self.plugin_dir, 'i18n',                                       'townplanning_{}.qm'.format(locale))                  if os.path.exists(localePath):                 self.translator = QTranslator()                 self.translator.load(localePath)                      if qVersion() > '4.3.3':                     QCoreApplication.installTranslator(self.translator)                  self.username = None             self.data = None                   def initGui(self):             # Create action that will start plugin configuration             self.toolbar = self.iface.addToolBar("TownPlanningToolbar")             self.action = QtGui.QAction(                 QtGui.QIcon(":/plugins/townplanning/icon.png"),                 u"Town Planning", self.iface.mainWindow())             # connect the action to the run method             self.action.triggered.connect(self.run)                  # Add toolbar button and menu item             self.toolbar.addAction(self.action)             self.iface.addPluginToMenu(u"&Town Planning", self.action)              def unload(self):             # Remove the plugin menu item and icon             self.iface.removePluginMenu(u"&Town Planning", self.action)             del self.toolbar              # run method that performs all the real work         def run(self):             self.dlg = TownPlanningDialog(self.iface, layeradded)             self.dlg.show()             self.dlg.repaint()      townplanningdialog.py:               class TownPlanningDialog(QtGui.QMainWindow):         def __init__(self, iface, layer):             QtGui.QMainWindow.__init__(self)             # Set up the user interface from Designer.             self.ui = Ui_TownPlanning()             self.ui.setupUi(self)             self.iface = iface             self.canvas = self.iface.mapCanvas()             self.band = None                  self.ui.previewCanvas = QgsMapCanvas(self.ui.widget)             self.ui.previewCanvas.setCanvasColor(QtCore.Qt.white)             self.ui.previewCanvas.enableAntiAliasing(True)             self.ui.previewCanvas.resize(self.ui.widget.size())                  self.ui.previewCanvas.setLayerSet([QgsMapCanvasLayer(layer)])             self.zoom_extents_preview_canvas()                  self.treeWidget = self.ui.treeWidget             self.lineEditFilter = self.ui.lineEdit             self.dataLabel = self.ui.label             self.zoomAllButton = self.ui.pushButton             self.viewBAButton = self.ui.pushButton_2             self.zoomSelectionCheckBox = self.ui.checkBox             self.noRecordsLabel = self.ui.label_3             self.searchMethodCombo = self.ui.comboBox             self.zoomInButton = self.ui.pushButton_3             self.zoomOutButton = self.ui.pushButton_4             self.panButton = self.ui.pushButton_5      ui_townplanning.py:               class Ui_TownPlanning(object):         def setupUi(self, TownPlanning):             TownPlanning.setObjectName(_fromUtf8("TownPlanning"))             TownPlanning.resize(793, 605)             TownPlanning.setCursor(QtCore.Qt.ArrowCursor)             self.centralwidget = QtGui.QWidget(TownPlanning)             self.centralwidget.setObjectName(_fromUtf8("centralwidget"))             self.lineEdit = QtGui.QLineEdit(self.centralwidget)             self.lineEdit.setGeometry(QtCore.QRect(10, 60, 161, 21))             self.lineEdit.setText(_fromUtf8(""))             self.lineEdit.setEchoMode(QtGui.QLineEdit.Normal)             self.lineEdit.setObjectName(_fromUtf8("lineEdit"))             self.treeWidget = QtGui.QTreeWidget(self.centralwidget)             self.treeWidget.setGeometry(QtCore.QRect(10, 80, 161, 192))             self.treeWidget.setProperty(_fromUtf8("showDropIndicator"), False)             self.treeWidget.setDragEnabled(False)             self.treeWidget.setDragDropOverwriteMode(False)             self.treeWidget.setDragDropMode(QtGui.QAbstractItemView.NoDragDrop)             self.treeWidget.setAlternatingRowColors(True)             self.treeWidget.setSelectionMode(QtGui.QAbstractItemView.SingleSelection)             self.treeWidget.setObjectName(_fromUtf8("treeWidget"))             self.treeWidget.header().setVisible(False)             self.treeWidget.header().setSortIndicatorShown(True)             self.groupBox = QtGui.QGroupBox(self.centralwidget)             self.groupBox.setGeometry(QtCore.QRect(380, 20, 401, 101))             self.groupBox.setObjectName(_fromUtf8("groupBox"))             self.label = QtGui.QLabel(self.groupBox)             self.label.setGeometry(QtCore.QRect(10, 20, 381, 71))             self.label.setText(_fromUtf8(""))             self.label.setObjectName(_fromUtf8("label"))             self.widget = QtGui.QWidget(self.centralwidget)             self.widget.setGeometry(QtCore.QRect(380, 160, 400, 400))             self.widget.setObjectName(_fromUtf8("widget"))             self.pushButton_5 = QtGui.QPushButton(self.widget)             self.pushButton_5.setGeometry(QtCore.QRect(50, 10, 24, 24))             self.pushButton_5.setText(_fromUtf8(""))             self.pushButton_5.setObjectName(_fromUtf8("pushButton_5"))             self.pushButton_3 = QtGui.QPushButton(self.widget)             self.pushButton_3.setGeometry(QtCore.QRect(2, 10, 24, 24))             self.pushButton_3.setText(_fromUtf8(""))             self.pushButton_3.setObjectName(_fromUtf8("pushButton_3"))             self.pushButton_4 = QtGui.QPushButton(self.widget)             self.pushButton_4.setGeometry(QtCore.QRect(26, 10, 24, 24))             self.pushButton_4.setText(_fromUtf8(""))             self.pushButton_4.setObjectName(_fromUtf8("pushButton_4"))             self.pushButton = QtGui.QPushButton(self.centralwidget)             self.pushButton.setGeometry(QtCore.QRect(10, 330, 75, 23))             self.pushButton.setObjectName(_fromUtf8("pushButton"))             self.pushButton_2 = QtGui.QPushButton(self.centralwidget)             self.pushButton_2.setGeometry(QtCore.QRect(10, 300, 101, 23))             self.pushButton_2.setObjectName(_fromUtf8("pushButton_2"))             self.checkBox = QtGui.QCheckBox(self.centralwidget)             self.checkBox.setGeometry(QtCore.QRect(10, 280, 111, 18))             self.checkBox.setObjectName(_fromUtf8("checkBox"))             self.comboBox = QtGui.QComboBox(self.centralwidget)             self.comboBox.setEnabled(True)             self.comboBox.setGeometry(QtCore.QRect(10, 30, 161, 22))             self.comboBox.setObjectName(_fromUtf8("comboBox"))             self.comboBox.addItem(_fromUtf8(""))             self.comboBox.addItem(_fromUtf8(""))             self.label_2 = QtGui.QLabel(self.centralwidget)             self.label_2.setGeometry(QtCore.QRect(10, 10, 71, 16))             self.label_2.setObjectName(_fromUtf8("label_2"))             self.label_3 = QtGui.QLabel(self.centralwidget)             self.label_3.setGeometry(QtCore.QRect(20, 90, 101, 16))             self.label_3.setObjectName(_fromUtf8("label_3"))             TownPlanning.setCentralWidget(self.centralwidget)             self.menubar = QtGui.QMenuBar(TownPlanning)             self.menubar.setGeometry(QtCore.QRect(0, 0, 793, 18))             self.menubar.setObjectName(_fromUtf8("menubar"))             self.menuSearch = QtGui.QMenu(self.menubar)             self.menuSearch.setObjectName(_fromUtf8("menuSearch"))             TownPlanning.setMenuBar(self.menubar)             self.statusbar = QtGui.QStatusBar(TownPlanning)             self.statusbar.setObjectName(_fromUtf8("statusbar"))             TownPlanning.setStatusBar(self.statusbar)             self.action1 = QtGui.QAction(TownPlanning)             self.action1.setObjectName(_fromUtf8("action1"))             self.action2 = QtGui.QAction(TownPlanning)             self.action2.setObjectName(_fromUtf8("action2"))             self.menuSearch.addAction(self.action1)             self.menuSearch.addAction(self.action2)             self.menubar.addAction(self.menuSearch.menuAction())                  self.retranslateUi(TownPlanning)             QtCore.QMetaObject.connectSlotsByName(TownPlanning)


title: Exclude any geometry contained within another in geodjango queryset

text: I am trying to get a set of non-intersecting model objects my a GeoDjango project. I have a model with a polygon geometry field, and it represents a collection of buffers. The model definition is as follows:               from django.contrib.gis.db import models     class MyModel(models.Model):         buffer = models.PolygonField()         objects = models.GeoManager()      Now I would like to get a queryset that gives me all instances of `MyModel` where the buffer polygons are not completely contained within another buffer polygon of any other model instance. I managed to get such a set with the following loop:               from django.db.models import Q     nonintersecting = MyModel.objects.all()     for area in MyModel.objects.all():         nonintersecting=nonintersecting.exclude(Q(buffer__within=area.buffer) & ~Q(id=area.id))      I have a large data set and would like to avoid this loop. I tried using F expressions, but those only compare within the same model instance so that does not work.               nonintersecting=MyModel.objects.exclude(Q(buffer_orig__within=F(buffer_orig)) & ~Q(id=area.id))      My question is if there is a way to do this without the explicit loop shown above?


title: SQL Server query, simillar to Symmetrical Difference in ArcMap

text: i have two tables with geometry data in SQL Server. How can I get all geometry differences with spatial query? I tried:               select k.geom.STSymDifference(r.geom).ToString() from table1 k, table2 r      But it returns to many rows. I would like to get same result as I get it if I use Symmetrical Difference function in ArcMap.


title: How to let user only rearrange layers in a group layer?

text: We have a stand-alone application were we let the users add local shape and raster layers. We would like to add functionality so that the user can rearrange the added layers within a specific group layer. We know it's possible to enable drag and drop on the TocControl but it lets the user rearrange all layers and that's not an option for us. Is there a way to enable drag and drop only on a specific group layer or somehow mimic the same behavior?


title: GarminCustomMap: crash without error

text: My GarminCustomMap crashes without an error. Im using qgis 2.2 on Xubuntu Precise. I tryed with different settings but after 2 - 6 seconds QGIS crashes without any error. What could be the problem? Is there any other way to use my QGIS-Map with Garmin?


title: Link Google Earth with QGis

text: I was wondering if it was possible to link Google Earth files (KML, KMZ) to a QGIS project. The main idea would be that if you'd add a point with coordinates in Google Earth, this point would be added to your QGis project or vice versa. Can this be done?


title: formula for GWPCA cross validation on GWmodel package on R software

text: i want to know the formula of for GWPCA cross validation to get the best bandwuth on GWmodel package on R software. thanks


title: Render tiles in the sea

text: I installed a mapserver on my own server with mapnik, apache and mod_tile. I build an application with WorldWind and my map server. My problem is that there are blank tiles instead of blue tiles to represent the sea. Is there a simple way to remplace these tiles? Or do I have to change files like _generate_tiles.py_ for example?


title: QGIS crash when selecting a Select Single Feature in attribute table

text: I'm loading several layers, styles and setting up some layer relations in a plugin that I execute. Once this is done I manually 'Select Singe Feature' on my map, switch to the "non-table view" in the attribute table (Fig 1.). The next step is where the crash occurs. I click the single selected record in Fig. 1 and the QGIS crash message appear (Fig. 2). How can I figure out what goes wrong, and perhaps fix it? E.g., what to look for in the dump file, or where to send the dump file for a possible bug report. ( **NB.** To be explicit, this is not an April's Fools joke) ![Select Single Feature](http://i.stack.imgur.com/kMRzl.png) ![QGIS dump message](http://i.stack.imgur.com/l7qy4.png)


title: Updatecursor returns with 'variable not defined'

text: I am fighting a bit with the cursor function since this morning. Thanks to the help of the people here and some tutorials I found I got this far               import arcpy          arcpy.env.workspace = "C:\\Users\\Documents\\Scoring\\Ecodis_clip_split" #set workspace to folder with values     value = [] #create list for the returning value of searchCursor     eco = arcpy.SearchCursor ("Anagance.shp", "", "", "area_ha", "")      #search for the value in attributetable     for row in eco:     value.append(row) #write the value into the list          arcpy.env.workspace = "C:\\Users\\Documents\\Scoring\\Wetland"      #set workspace to folder with values     diss = arcpy.UpdateCursor ("diss_1.shp", "", "", "area_eco", "")     for row in diss:         row.setValue(area_eco, row.getValue(value)) #identify all the rows and Update die field area_eco with the extracted value from the search cursor out of the list         cursor.updateRow(row)          del eco     del diss      It bounces back with               area_eco is not defined      Which means - if I understood the concept of python correct - that this is not a variable. Ture, its my field which needs to updated. What am I doing wrong?


title: How to extend lines from "station line" plugin in QGIS

text: I'm trying to create cross sections using the station line plugin (and then import them to HEC RAS). But I can't seem to create lines longer than 99.99 meters from each side of my river line. Is this a bug?


title: FME converter from DWG to Openlayers supported format

text: FME converter is best software that i found till now to convert DWG to some open format vector data file supported by OpenLayers , but it got very awful convert from DWG to KML/GML it drawing unreadable lines on the map. that my input DWG and how it looks in FME: ![enter image description here](http://i.stack.imgur.com/Hu45g.png) After the conversion and read it via OpenLayers it looks like: ![enter image description here](http://i.stack.imgur.com/GQ3dU.jpg) Same thing even if i opened the output KML via FME: ![enter image description here](http://i.stack.imgur.com/zAMjI.png) I use the same projections. Any suggestions to solve it?


title: Terraserver replacment - source for US Topo maps?

text: Now that Terraserver has been shut down (or at least that's how it appears), is there an alternative WMS server for US topo maps? Probably too much to hope to find one that provides vector layers rather than rasters. About a year ago this similar question appeared (actually bulk download of topo maps), but back then Terraserver was still in operation and appeared in several of the answers.


title: convex hull over selected features

text: I am trying to achieve a `convex hull over a feature Set` returned by a query. But have been facing the error below which i am not able to figure out In order to debug to see how the convex hull works i had mocked up features and ran the convex hull for the created poly line features as seen in the below code. This works fine, But the above approach of plotting a convexhull for the feature set fails with the below error. Error: 1) `OPTIONS file:///C:/wip/maps/code/test/proxy.ashx?http://karthikpc:8399/arcgis/rest/se...ver/convexHull Resource failed to load` 2) Another is a type error               TypeError     arguments: Array[2]     get message: function () { [native code] }     get stack: function () { [native code] }     set message: function () { [native code] }     set stack: function () { [native code] }     type: "non_object_property_load"     __proto__: Error      Have given my function below.               function bounding_polygon(featureSet){         var geometries = [];         var geometries_1 = [];         var geoms=[];     var symbol =  new              esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID, new           dojo.Color([255,0,0,0.85]), 2);              //getting geometry from feature set         var geoms = dojo.map(featureSet.features, function(graphic) {          return graphic.geometry;          });         map.graphics.clear();             //function that fails           gsvc.convexHull(geoms,function store_value(result){          var symbol;         switch(result.type){               case "point":                 symbol = new esri.symbol.SimpleMarkerSymbol();                 break;               case "polyline":                 symbol = new esri.symbol.SimpleLineSymbol();                 break;               case "polygon":                 symbol = new esri.symbol.SimpleFillSymbol();                 break;             }          map.graphics.add(new esri.Graphic(result,symbol));         },function error_log(error){});          //trying to run the same for manually created geometries          var new_polyline= new esri.geometry.Polyline(map.spatialReference);      var other_polyline= new esri.geometry.Polyline(map.spatialReference);     new_polyline.addPath([[103.76794615402221, 1.4416487206261546], [103.76790721700029, 1.4427548269642734], [103.76789991705209, 1.4429139656905441]]);      other_polyline.addPath([[103.77897404312868, 1.3687951231704805], [103.7789286010334, 1.36900185852713], [103.7788202535522, 1.3694226417416273]]);      geometries_1[0]=new_polyline;      geometries_1[1]=other_polyline;          gsvc.convexHull(geometries_1,function store_value(result){      var symbol;      switch(result.type){               case "point":                 symbol = new esri.symbol.SimpleMarkerSymbol();                 break;               case "polyline":                 symbol = new esri.symbol.SimpleLineSymbol();                 break;               case "polygon":                 symbol = new esri.symbol.SimpleFillSymbol();                 break;             }          map.graphics.add(new esri.Graphic(result,symbol));              },function error_log(error){});                   }


title: Aggregating Polygons To Meet Privacy Requirements

text: I have a point feature class representing office locations of all employers in a certain industry. The feature class has an attribute for storing the number of employees working at each office. Someone has requested to use this data, spatially-joined to the smallest geographic unit possible — Census Blocks, in this case. However, a privacy agreement prevents the data's release as-is. Instead, it must be suppressed to meet two criteria:   1. Any polygon must contain at least 3 employers (points);   2. No more than 80% of the total employment within a polygon may be by a single employer. I have successfully written a script that spatially joins the points to Census Blocks, keeping the sum and maximum employment in each. Each one that does not meet the suppression criteria is flagged. (Polygons containing no points are not flagged, since there is no data to suppress.) I then check each Block _Group_ to see whether any flagged Blocks are contained within it. Block Groups containing only un-flagged Blocks are then replaced with the Blocks. The resulting features class is then checked against the suppression criteria, to check whether the Block Groups have adequately suppressed the data. The same process is repeated for Tracts, leaving me with a dataset consisting of Tracts (some flagged and some not), Block Groups and Blocks (all un- flagged). The next progression in the geographic hierarchy, however, is the county, which is of no use to the person requesting this data. My question, then, is this: **Are there any commonly accepted methods of aggregating polygons into as many groups as possible, so that all meet some minimum criteria?** Here are some rules that I would like to apply to the aggregation:   * Whenever possible, flagged Tracts should only be aggregated with other flagged Tracts;   * For the flagged Tracts that are not contiguous with any others (or isolated groupings that still do not meet the criteria), they can be joined with Tracts that already meet to criteria, although there may be Tracts with no employers in between them that will also need to be included.   * I would like to keep county boundaries intact unless absolutely impossible (and I anticipate doing this by separating the inputs features into their respective counties before processing them).   * The solution must be in Python, with the use of ArcGIS tools or open-source Python libraries. Ideally, someone can point me to an existing means of implementing this aggregation. If not, I am happy to code the algorithm myself, although a list of specific steps/tools would be much appreciated. The problem strikes me as a special case of redistricting (with dis-contiguous polygons), and to this end I have looked into using PySAL's regionalization algorithms, although it is not clear to me how to check the maximum employer's percentage of total employees using these.


title: Retain ESRI symbology in Mapinfo TAB

text: Is there anyway to retain/easily transfer symbology from ESRI to Mapinfo? The only process I found was posted in 2010 - https://groups.google.com/group/mapinfo-l/browse_thread/thread/883223ff4ae56809 and it is rather convoluted and supposed to be unreliable. ArcGIS 10 and Mapinfo 11 are the core programs. This is related to Is it possible to convert an ESRI style sheet to a font set? and Create a summary table of the symbology of all datasets in multiple mxd's as possible workarounds.


title: Publishing oracle spatial table in geoserver

text: I published oracle **view** in geoserver. using this link. But when I see this in layer preview I get following xml response:               <?xml version="1.0" encoding="UTF-8" standalone="no"?><!DOCTYPE ServiceExceptionReport SYSTEM "http://192.168.70.65:80/geoserver/schemas/wms/1.1.1/WMS_exception_1_1_1.dtd"> <ServiceExceptionReport version="1.1.1" >   <ServiceException code="internalError">           Rendering process failed     java.io.IOException     null     ORA-00904: "NEPALCLUSTER": invalid identifier          </ServiceException></ServiceExceptionReport>                     CREATE OR REPLACE FORCE VIEW "myview"."d_DAILY" ("nepalCluster", "nepalclutter", "ZONE", "MI_PRINX", "MI_STYLE", "GEOLOC", "VAL") AS        select        c.nepalCluster,c.nepalclutter,c.zone,c.mi_prinx,c.mi_style,c.geoloc,NVL(coalesce(optz.VAL,opth.VAL),-1) VAL     from      CLUSTER_REGION c     left join      (select      *      from gDaily a) optz     on      optz.CELL_CLUSTER = c.nepalCluster     left join     (select      *      from zDaily b) opth     on      opth.CELL_CLUSTER = c.nepalCluster;


title: Reading OGC geometry data from SQL server using R

text: I would like to write a script in R that creates a map using OGC geometry data that is stored on a Microsoft SQL server. Is there a way to read OGC geometry data using a query passed from R? Maybe the code would look something like this (But not like this because the following code is actually garbage).               library(RODBC)     library(maptools)     library(maps)          png(file="example%02d.png", width=600, height=480)          con <- odbcDriverConnect('driver={SQL Server};server=SERVERNAME;database=DBNAME;trusted_connection=true')          objects_1 <- sqlQuery(con, 'SELECT OBJID, Shape FROM TABLENAME;')          spplot(objects_1, col="#000000FF", sp.layout = list(otherObjects))


title: OpenLayers serving a MapServer layer: data won't display correctly without using WMS

text: I am trying to serve my Mapfile data as a MapServer layer in OL, rather than a WMS layer. When I display it using WMS, everything works fine. The part of the mapfile I have struggled with was SIZEUNITS. Seems easy enough, but I encountered some strange behavior along the way. The issue: my SIZEUNITS units is miles and the size is 100 (so a 50 mile radius). When I use the WMS layer below, everything is fine:               var cust = new OpenLayers.Layer.WMS("POIs",     "http://192.168.2.164/cgi-bin/mapserv.exe?map=/data2/html/html2/david_POIs4.map",     {layers: ['POIs']}, {wrapDateLine: true, singleTile: true});      However, when I try to switch it over to:               var cust = new OpenLayers.Layer.MapServer("POIs",     "http://192.168.2.164/cgi-bin/mapserv.exe",     {map: "/data2/html/html2/david_POIs4.map",     layers: ['POIs']}, {singleTile: true} );       The points (circles rather, with points in the middle) are HUGE. More or less taking up the whole map. So, why would switching from WMS to MapServer show such glaring differences? A MapServer layer should be simpler to correctly configure. When I was stuck on a similar issue last week I posted here (of course hadn't gotten this far yet). Edit: see the images here for a visual. First map is correct (zoomed in), 2nd map is incorrect. Edit2: I've never had an issue with projections before (I have the same projection set in both the `LAYER` and the `MAP`). Here's the layer in question, and if I use `SIZEUNITS MILES` it seems to automatically get blown out of proportion, regardless of the value you pass to `SIZE` within `STYLE`.                   LAYER                 NAME POIs                 DEBUG 5                 TYPE POINT                 CONNECTIONTYPE OGR                 CONNECTION 'PointsOfInterest_Latest.db'                 DATA 'POIs'                 STATUS ON                 SIZEUNITS miles                 LABELITEM "name"                 CLASS                         STYLE                                 SYMBOL 'Airport'                                 COLOR 8 29 118                                 #WIDTH 5                                 SIZE 10                                 #MINWIDTH 1                                 #MAXWIDTH 1                         END                         #STYLE                         #       SYMBOL 'City'                         #       COLOR 8 29 118                         #       SIZE 6                         #END                         LABEL                                 FONT "andalewt"                                 TYPE truetype                                 ENCODING "UTF-8"                                 MINSIZE 9                                 MAXSIZE 9                                 POSITION AUTO                                 PARTIALS FALSE                                 COLOR 8 29 118                                 OFFSET 1 1                                 PRIORITY 10                                 BUFFER 0                         END # LABEL                 END                      PROJECTION                         "proj=latlong"                         "ellps=WGS84"                         "datum=WGS84"                 END         END      Now, you can see the `SIZEUNITS` set to `MILES` and the `SIZE` set to 10. Being the farthest zoom level out, you should barely be able to see the circles, but they are huge, as you can see here. Any input is appreciated.


title: Can't edit osm_roads in udig

text: There is a PostGIS database that has a road network imported from Open Street Map, stored in osm_roads table. For the testing purposes I need to change some road geometry there. I tried to edit it with uDig Version: 1.4.0b. When I select the layer with osm_roads, I don't see any 'edit' buttons; the only buttons I see are 'create arc', but after drawing the arc it disappears, so I can't edit the road network anyway. Please advise what I should configure to make this changes possible. What is the most likely reason why uDig doesn't allow to edit the road layer?


title: Can I use geoserver to bring sql layer as vector so I can edit it?

text: I've updated geoserver with SQL Store and published a table of geometries. When using WMS I get the LineStrings as img so I can't edit them. I've read that WFS bring the data as vector (features) but I wasn't able to find example how to do it (When I preview layer in OpenLayers it comes as WMS only) How can I do it ?


title: Why does ogr2ogr SHP import error out trying to create a database?

text: All, Can anyone tell me what the error is with OGR2OGR when trying to import a .SHP into MSSQLSPATIAL?               ERROR 1: MSSQL Spatial driver doesn't currently support database creation.     Please create database with the Microsoft SQL Server Client Tools.     MSSQLSpatial driver failed to create MSSQL:server=[]\[];database     =[];uid=[];pwd:[]      The import statement I'm using is this:               ogr2ogr -f MSSQLSpatial "MSSQL:server=[]\[];database=[];uid=[];pwd:[]" C:\[path to file]\police_districts_2877.shp -a_srs EPSG:2877 -lco "GEOM_TYPE=geometry" -lco "GEOM_NAME=geom" -nln "CCD_PoliceDistricts" -progress      As you can see, I'm not trying to create a database, and there is one already created. as per the 'creation issue' cited here: http://www.gdal.org/ogr/drv_mssqlspatial.html I have followed several posts, including this: Shapefile to MSSQL with ogr2ogr fails to make connection, to no avail... EDIT: I'm also getting a connection error, but I can't see why...I think my syntax is correct:               ERROR 1: Unable to initialize connection to the server for MSSQL:server=[]\[];database=[];uid=[];pwd:[],[Microsoft][ODBC SQL Server Driver][SQL Server]Login failed for user '[username]'.      Thanks! -m


title: Is SpatiaLite or SQLite Spatial an Accepted Standard?

text: The reason why I ask this is because I have been undertaking some testing and have found that there seems to be many implementation and not all of them are interchangeable. For example:   1. FME can read and write SQLite Spatial (FDO) and SpatiaLite   2. MapInfo can read and write SQLite Spatial (FDO, I think)   3. QGIS can read and write SpatiaLite But,   * MapInfo wont read QGIS' SpatiaLite   * FME will read and write anything (as usual)   * QGIS will read output from FME in either the FDO or SpatiaLite writer, but can't edit (silent fail) and the same for MapInfo written tables. There are probably other bits of software that read/write some variation of the above but implement it differently. So, to ask a simple question. Is it fair to expect to work in some variation of SQLite and be interoperable? Edit: It turns out that you can read/write spatialite from MapInfo and FME in QGIS via the 'open vector' dialog. I was trying using the 'add spatialite table' dialog. It just remains that MapInfo wont read SpatiaLite from QGIS, but this isn't all that surprising... Edit 2: I am thinking that QGIS wont edit mapinfo and qgis SQLite Spatial layers, despite being able to view them, due to the use of the generic 'geometry' spatial type, rather than 'point' or 'polygon' etc.


title: Tilecache and Mapnik problem

text: **My configs** :               # Rendering OpenStreetMap data with Mapnik; should use metaTiling to     # avoid labels across tile boundaries     [osm-map]     type=MapnikLayer     mapfile=/usr/home/sergeev/osm_mapnik/osm.xml     spherical_mercator=true     tms_type=google     metatile=yes     extent_type=loose      and               var map, layer;                  function init(){                 var options = {                         resolutions: [156543.0339, 78271.51695, 39135.758475, 19567.8792375, 9783.93961875, 4891.969809375, 2445.9849046875, 1222.99245234375, 611.496226171875, 305.7481130859375, 152.87405654296876, 76.43702827148438, 38.21851413574219, 19.109257067871095, 9.554628533935547, 4.777314266967774, 2.388657133483887, 1.1943285667419434, 0.5971642833709717, 0.29858214168548586],                             numZoomLevels: 20,                             maxExtent: new OpenLayers.Bounds(0.05917265625,0.1171263671875,0.0593009765625,0.11720234375),                             projection: 'EPSG:900913'                 };                      map = new OpenLayers.Map( $('map'), options);                 layer = new OpenLayers.Layer.WMS( "VMap0",                          "tilecache.cgi?", {layers: 'osm-map', format: 'image/png' } );                 map.addLayer(layer);                 map.addControl(new OpenLayers.Control.Permalink());                 if (!map.getCenter()) map.zoomToMaxExtent();             }      my `osm.xml` is here: http://pastebin.com/f5EWUsTy **My problem** : I see the map, but i cant't move it: each zoom level is fixed. And the center of the map is somewhere in Africa, but i thought my xml is describing north of Russia. The image below is what i see: ![enter image description here](http://i.stack.imgur.com/WRxNI.png) So what should i do? (I'm totaly new to gis, so, please, be patient :)


title: Advice on how to present a cross section

text: So, I'll try to best explain what I want to do. I have a DEM of a volcano (only showing half the DEM in pictures below) and I have identified a lava flow that I would like to make a cross section/profile of. I have a plane (shown in the pink/red) which represents the base of the flow. I would like to get the flow by its self, ideally with profiles taken along the length of the flow as well as 3 or 4 across the flow. I would also like to show heights between the plane and the DEM. I am currently working in ArcScene 10. I have all the extensions available but can not add new ones as I am on a university computer and do not have administrator privileges. The final product must be visually easy to look at as the end product will eventually be used in a presentation. Any ideas are very welcome! I have a couple of days to play around with this :D![view of the flow from above](http://i.stack.imgur.com/OuY1m.jpg)![View of the flow from side. DEM transparency=30%](http://i.stack.imgur.com/qRyS7.jpg)![Crazy profile](http://i.stack.imgur.com/uMPvT.jpg)


title: How to acquire an ArcGIS access token with Windows Authentication (single sign on)

text: Consider the following scenario:   1. The end user is running an in-house build Windows desktop client that uses the ArcgGIS WPF library for rendering ArcGIS maps.   2. The end user is authenticated via Windows authentication.   3. Now I want the client to acquire an access token for the REST services _without_ prompting the end user for his/her username/password. That is, given that the end user is already logged in via Windows AD, the client should be able to acquire an access token with NTLM/Kerberos or what ever it is Windows authentication is using.   4. The acquired token is then send to an in-house web service which will then access the ArcGIS REST services with the token and thus impersonating the end user without knowing either username or password (this web service is _not_ part of the Windows domain). Is that possible? How? The server is "ArcGIS REST API for 10.1 SP1". I am the client/server programmer - not the ArcGIS administrator :-)


title: ArcPy / DDP - alter graticule logically in varying map scales

text: Morning -this question appears on many sites and I was wondering if there was a script that is now available that resolves it? Please also upvote ideas.arcgis.com/ideaView?id=087300000008MmJAAU#comments My workarounds have been to create multiple mxd's that just hold the same scale map sets and then export only the relevant pages to that scale but this adds a lot of time to the process when changes are required as they have to be done on each mxd. The other related issue is with the centering of elements where when a map is 1:10000 and another is 1:100000 one of them will be off centered from other elements such as the "N" Arrow. Best,


title: How to split non-contiguous feature in QGIS?

text: I have a shapefile with a feature that is non-contiguous. I want to split the areas on the left from the area on the right. However, QGIS's Split Features tool requires you to draw a line over the feature. When I draw a line between the two areas of this feature, it doesn't do anything. How should I split this feature? Here's a screen shot. The feature I want to split is in yellow. ![qgis screen shot](http://i.stack.imgur.com/t1vb6.png)


title: How to fill "shapefile name" and "name of attribute"?

text: When I try to create a shapefile two of the boxes to fill in are "name of the attribute" and " shapefile name". this may be very basic but I can't find and answer. I want to fill them in correctly. What are their purposes and what should one put in them? Why are their two names for the file? Thanks in advance.


title: iterating and exporting all combinations of features as shapefiles

text: I am trying to export individual shapefiles of all 2 pair combinations from a multi-feature polygon shapefile. To explain this better as set of FID's 0, 1, 2, 3, 4 would be exported as individual shapefiles of FID pairs: ('0', '1'), ('0', '2'), ('0', '3'), ('0', '4'), ('1', '2'), ('1', '3'), ('1', '4'), ('2', '3'), ('2', '4'), ('3', '4') Im not so good with python so my attempts with intertools.combination etc have not worked out. Any help would be appreciated. Thanks for any advice and help.


title: How to show cell boundaries on a raster?

text: I am trying to recreate a situation where each cell from a raster dataset has a black border as seen in this image. Any ideas? ![illustration of raster with visible cell boundaries](http://i.stack.imgur.com/2GDJM.png)


title: Problems displaying reprojected WMS in QGIS 2.0.1

text: I'm loading a WMS layer in QGIS 2.0.1 using the coordinate reference system SWEREF99 TM EPSG:3006, a Swedish crs. The layer displays in the map but when I zoom in to Sweden using the Zoom in-button I end up somewhere in the Arabian Sea. The server (a GeoServer-instance) supports and serves the layer in EPSG:3006. I've tried different settings for reprojecting but no differnces. In QGIS 1.8 this works perfect on the same layer. URL: http://geodata.havochvatten.se/geoservices/hav-bakgrundskartor/ows Layer: hav-grundkarta


title: Is there a way to change Table to Table import parameters in ArcCatalog?

text: Is there a way to change Table to Table import parameters in ArcCatalog to either import text longer then 255 characters or truncate the input text to 255 characters? I have 350 tables in a CSV format and I am using import Table (multiple) function in ArcCatalog. Majority of the table import fine but some fail because the source field text is longer then 255 characters. Instead of making the text field longer automatically ArcCatalog importer fails with the following error (identical in all cases):               Failed to convert C:\Users\Jakub Sisak\Desktop\LIDS_Jan2013\Processed_TBL\arasurpt_ara_survey_point.csv. ERROR 001156: Failed on input      OID -1, could not write value 'Area upstream of culvert appears to have been ponded recently, possibly due to culvert blockage/beaver      activity. Sampling started upstream of ponded area. Only fish sampling was done. Station UTM is downstream end. Pond at upstream is      man-made dam of earth and debris and is in poor repair.' to output field COMMENTS     Failed to execute (CopyRows).      Is there a way to change the "importer" paramaters to either truncate the string (accaptable) or increase the length of the text field (ideal) I am using File Geodatabase, ArcGIS Standard (ArcEditor). There is no memo field option in FGDB but text field (default length 255 characters) can be more then 255.


title: When was MFD renamed to FD8?

text: According to my understanding, the Multiple Flow Direction (MFD) algorithm was originally defined in _The Prediction of Hillslope Flow Paths for Distributed Hydrological Modelling using Digital Terrain Models_ by Quinn et al. 1991. Since then it seems to have been renamed to FD8 to prevent confusion and I believe some may even refer to it as MD8 but I have failed to discover when this happened. This originally caused some confusion on a project I am working on. I am looking for a reference or an explanation for when MFD became known as FD8 and if MD8 is also a valid name for the algorithm or am I mixing algorithms altogether?


title: GDALWARP and OGR2OGR / CS2CS not agreeing on transformation

text: I am trying to transform some data (both raster and vector) from WGS84 to a Croatian Coordinate system defined by:               +proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +towgs84=550.499,164.116,475.142,5.80967,2.07902,-11.62386,0.99999445824 +units=m +no_defs      This works great in CS2CS (Rel 4.8.0):               cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +to +proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +towgs84=550.499,164.116,475.142,5.80967,2.07902,-11.62386,0.99999445824 +units=m +no_defs -r      And also great in OGR2OGR (version 1.9.2):               ogr2ogr -s_srs "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" -t_srs "+proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +towgs84=550.499,164.116,475.142,5.80967,2.07902,-11.62386,0.99999445824 +units=m +no_defs" out.shp  in.shp       But in GDALWARP (version 1.9.2), the data comes out c.400m offset to the West (I think it's possibly a rotation error based upon it's appearance):               gdalwarp -s_srs "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" -t_srs "+proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +towgs84=550.499,164.116,475.142,5.80967,2.07902,-11.62386,0.99999445824 +units=m +no_defs" in.tif out.tif      An example of this can be seen below, where I digitised a vector layer around existing raster data, transformed both using the above commands, and then loaded into QGIS: ![example of offset error](http://i.stack.imgur.com/c3eEO.png) I can't for the life of me figure out why this could be - it is worthy of note that I have tried this with multiple Raster datasets from a variety of sources, and they always seem to be offset to the West by the same amount. I have tried upgrading to GDAL 1.10.1, and this has no effect upon the output. Has anyone ever seen this before / got any idea why it could be? Thankyou!


title: vb.net arcobjects enable customized context menu item only if a feature is selected

text: The solution may be very simple, but I haven't made it work. I want to enable a data view context menu item only if a feature (in polyline layer) is selected, all other places the item should in grey. May I say what I need to do is check whether a feature is selected? I searched for this but don't find a proper solution. Or any other concerns? I am using ArcObjects 10.2.2 and VB.NET.


title: How do I load offline google map/earth vrt files as a backdrop in qgis

text: I was trying to follow the process : saving-google-earth-location in order to cache google maps/earth offline for digitizing, but I am fairly new to qgis. When I try to add the vrt file as a raster in qgis I get nothing but a blank image. I also noticed that pyramiding gave me an error which I suspect is because of the resolution. In addition if I wanted to load the vrt as a backdrop/base map does it not require a bunch of tiles loaded to create an image mosaic? I am working in Windows so I would prefer to not have to run some gdal commands in the command prompt. I'm also not so sure how to use tile caching and I would prefer the easiest way possible without having to standup any wms type of server. *I know its against the terms of service.


title: How to integrate pascal code into a QGIS model or plugin?

text: I have a console program that is written in pascal (lazarus) that reads in some raster files and parameter values as input. The code use dynamic 2D arrays to represent the rasters. I would like to know how I can integrate this lazarus (pascal) code into QGIS.


title: Check if feature is selected (Openlayers)

text: I wonder if there is a way to check whether a pointfeature (OpenLayers.Feature.Vector) is selected. I read the API, but didn't find any solution. Any idea?


title: Best way to clip a large raster ECW file ?

text: I am trying to clip a large ECW (details below) but the raster file is too big to be processed entirely. Some details of the ECW below Driver: ECW/ERDAS Compressed Wavelets (SDK 5.0) file size: 50gb Size is 450000, 565081 Pixel size: 0.15 0.15 COLORSPACE=RGB COMPRESSION_RATE_TARGET=9 VERSION=2 Number of bands: 4 The area I want to clip is roughly a 1/5 of the original file. Here are the methods I have tried out with no success:   1. Used Arcgis to save the ecw to tiff / other formats... (I quickly gave up)   2. Used Qgis and its clipper tool... the file creation stayed stuck at about 40%.   3. Used gdal_translate out of OSGeo4W with other options than Qgis. (Tried that thinking that MAYBE freeing some memory not using Qgis would the trick)   4. Used gdal_retile thinking that I would cut the image into pieces and grab the one I wanted. The command "gdal_retile -ps 10000 10000 -of ecw -tileIndex tile.shp -targetDir input.ecw This crashed even quicker" Does anyone have an idea? For info I run windows 7 64bits on an i5-3470 3.2Ghz with 16gb of ram.


title: Converting projected geoTiff to WGS84 with GDAL and Python

text: Apologies if the following question is somewhat stupid, but I am only VERY new to this whole GIS thing. I am trying to convert some projected geoTiff images to WGS84 using gdal in python. I have found a post that outlines the process to transform points within projected GeoTiffs using something similar to the following:               from osgeo import osr, gdal          # get the existing coordinate system     ds = gdal.Open('path/to/file')     old_cs= osr.SpatialReference()     old_cs.ImportFromWkt(ds.GetProjectionRef())          # create the new coordinate system     wgs84_wkt = """     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0,             AUTHORITY["EPSG","8901"]],         UNIT["degree",0.01745329251994328,             AUTHORITY["EPSG","9122"]],         AUTHORITY["EPSG","4326"]]"""     new_cs = osr.SpatialReference()     new_cs .ImportFromWkt(wgs84_wkt)          # create a transform object to convert between coordinate systems     transform = osr.CoordinateTransformation(old_cs,new_cs)           #get the point to transform, pixel (0,0) in this case     width = ds.RasterXSize     height = ds.RasterYSize     gt = ds.GetGeoTransform()     minx = gt[0]     miny = gt[3] + width*gt[4] + height*gt[5]           #get the coordinates in lat long     latlong = transform.TransformPoint(x,y)       My question is, if I want to convert these points and create a new WGS84 GeoTiff file, is this the best way to go about it? Is there a function that exists that will do such as task in 1 step? Thanks!


title: Using ArcGIS Online Base Maps without Esri Software?

text: Since the end of January 2011: > "ArcGIS Online map services are available at no cost regardless of use." - > http://www.esri.com/software/arcgis/arcgisonline/standard-maps.html > > However, the use of ArcGIS Web Mapping APIs may require a subscription. - > http://www.esri.com/software/arcgis/arcgisonline/common-questions.html All the basemaps available from ESRI are accessible by REST services which do not require an API. There is already a new layer type in OpenLayers that can access these services. The Terms of Service state: > In particular, you may not (or attempt to) (a) access Services or resources > not made available to you through a standard web browser or Esri products; **Update** \- these TOS _may_ in fact refer to the esri.com site and "services" rather than the map services themselves. So is it legal and within the terms of service to use any of these layers in _any_ web application such as OpenLayers or through the Google or Bing mapping APIs as these are all used in a standard web browser?