title: .NET Components triggering ITools on ArcEngine/Objects

text: I'm building a simple map interface on a .net component that I need to redistribute around. Well, my initial thoughts was to create the tools just the way I'm used to do it in ArcMap, inheriting BaseTool and then set the MapControl.CurrentTool = tool in my onclick form events. What is the best way to accomplish this? This one? Also, how can I pass the hook to those tools? I'm using this class:               public class CommandPoolManager     {         public ICommandPool _CommandPool;              public ICommandPoolEdit GetCommandPoolEdit         {             get { return (ICommandPoolEdit)_CommandPool; }         }              public ICommandPool GetCommandPool         {             get { return this._CommandPool; }         }              public CommandPoolManager(object hook)         {             this._CommandPool = new CommandPoolClass();             ICommandPoolEdit commandPoolEdit = (ICommandPoolEdit)this._CommandPool;             commandPoolEdit.SetHook(hook);         }              public void AddCommandToPool(ICommand command)         {             ICommandPoolEdit commandPoolEdit = GetCommandPoolEdit;             commandPoolEdit.AddCommand(command,null);             commandPoolEdit.CallOnCreate(command);         }              public void AddCommandToPool(ICommand command,UID uid)         {             ICommandPoolEdit commandPoolEdit = GetCommandPoolEdit;             commandPoolEdit.AddCommand(command, uid);             commandPoolEdit.CallOnCreate(command);         }     }      Would that be the correct way to use CommandPool? I'm without options for testing this, since it's a component to be distributed. The real question is: Am I on the right track?


title: Import PDF to qgis

text: i have a PDF containing a region of Baltimore city that contains lots of polygons in it. you can download it from here: http://www.baltimorecity.gov/Portals/0/agencies/planning/public%20downloads/Typology2011%20Final.pdf I also have a kmz file, containing each part of Baltimore city, you can download it here: https://data.baltimorecity.gov/api/file_data/-gbDlGngMlcHrfsfZepc1QzhJbN42zwUOcIhJ8ZJOVY?filename=Housing%2520Typology%25202008.kmz All i want to do is to open the kmz file in Qgis, then import the pdf file into qgis as background and then select the city regions with mouse and based on the color of the region in the PDF file, edit the attributes of the kmz file and after all i want to export it as an shape file. How am i able to import the pdf file as a background image and fits it to city regions in order to be able to select the regions regarding the background image?


title: How to get the layers of my project to show while tracking with GPS?

text: I've been using/testing several versions of various GIS software and run into occasional issues. Right now, my biggest with QGIS is with GPS connection/downloading from/tracking. These are all critical issues for me. 1.) I think I did finally figure out the need to set data transfer to NMEA in/out on my Garmin GPSMap60CSX. QGIS repeatedly crashed on clicking connect, I hope this was my only error and is now corrected. 2.) Downloading from/uploading to GPS unit... clicked to download tracks and waited 10 minutes before giving up and canceling with no success. Any tips/pointers? 3.) I've posted in the old forum about this one... I am unable to see any layers of the open project when live tracking. Last time (I think using Copiapo) I was eventually able to zoom until the .shp file was visible (basically stumbled onto that one), but not the .sid file that formed the base map. Now, using Wroclaw there is nothing at all visible except the gps track. Just creating a track doesn't help, I need to be able to see my location in relation to various physical features and property boundaries in real-time.


title: OpenLayers + GeoServer + Oracle: WFS-Transaction update error

text: I am trying to implement wfs-t in a map application. tried using ArcGIS server wfs-t, and didn't make it since there is a lock mechanism while trying to edit a feature.ArcGIS WFS-T + OpenLayers; then decided to move on using GeoServer... unfortunately, stuck again. sincerely wish somebody in this site can give me some advice. thanks a lot! **what I did in JS code is:** update the value of an attribute, and set the state of that feature as OpenLayers.State.Update;               feat.attributes['MATERIAL'] = 'CMPtest';     feat.data['MATERIAL'] = 'CMPtest';     feat.state = OpenLayers.State.UPDATE;      then: using `SaveStrategy.save()`; **the error response text is:**                <?xml version="1.0" ?>     <wfs:WFS_TransactionResponse version="1.0.0" xmlns:ogc="http://www.opengis.net/ogc" xmlns:wfs="http://www.opengis.net/wfs" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.opengis.net/wfs http://XXXX/geoserver/schemas/wfs/1.0.0/WFS-transaction.xsd">       <wfs:InsertResult>         <ogc:FeatureId fid="none"/>       </wfs:InsertResult>            <wfs:TransactionResult>              <wfs:Status>                <wfs:FAILED/>              </wfs:Status>              <wfs:Message>           update error         </wfs:Message>            </wfs:TransactionResult>     </wfs:WFS_TransactionResponse>      **the request xml is:(from Chrome console)**               <?xml version="1.0" ?>         <wfs:Transaction service="WFS" version="1.0.0" xmlns:wfs="http://www.opengis.net/wfs" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.0.0/WFS-transaction.xsd">           <wfs:Update typeName="feature:SAMPLE_FEATURE" xmlns:feature="http://www.opengeospatial.net/cite">             <wfs:Property>               <wfs:Name>                 SHAPE               </wfs:Name>               <wfs:Value>                 <gml:Point srsName="EPSG:26912" xmlns:gml="http://www.opengis.net/gml">                   <gml:coordinates cs="," decimal="." ts="">                     401550.55,4647703.82                   </gml:coordinates>                 </gml:Point>               </wfs:Value>             </wfs:Property>             <wfs:Property>               <wfs:Name>                 MATERIAL               </wfs:Name>               <wfs:Value>                 CMPtest               </wfs:Value>             </wfs:Property>             <wfs:Property>               <wfs:Name>                 CENTER_LON               </wfs:Name>               <wfs:Value>                 -112.188498069               </wfs:Value>             </wfs:Property>             <wfs:Property>               <wfs:Name>                 CENTER_LAT               </wfs:Name>               <wfs:Value>                 41.9749785934               </wfs:Value>             </wfs:Property>             <wfs:Property>               <wfs:Name>                 CENTER_Z               </wfs:Name>               <wfs:Value>                 1353.59340947               </wfs:Value>             </wfs:Property>             <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">               <ogc:FeatureId fid="SAMPLE_FEATURE.1"/>             </ogc:Filter>           </wfs:Update>         </wfs:Transaction>


title: Where can I find the highest resolution shapefiles for counties or states

text: I am looking for a download source which has the highest resolution (greatest detail/accuracy) for all counties in the US. More specifically I run into this problem mainly when dealing with the shorelines. I realize that one can download high resolution shoreline data, overlay the tiger shapefiles for the counties and then clip them but I figured someone might have already done this.


title: I would like to know about Careers in GIS

text: I am currently doing my MS software engineering. I am very much interested in GIS. I would like to know about the various opportunities available in this field and also the skills that i will have to acquire to enter into this field. i would like it to be software related. Thank you


title: How to get started on a GIS Programming Career?

text: I am already working as GIS Consultant in UK. Being an engineering consultancy (not IT or GIS based) company there wasn't any need of GIS development in the organisation. But due to some recent project requirements my company is trying to promote us towards GIS development. So my question is that, which route should I progress my career, like: Python - C# - .Net   OR Java? Quicker advice will be much appreciated as I am in India for my 3 weeks and have an opportunity to join some basic programming courses. Recent visit to the local computer institute, C# and .Net was the main stress.


title: Is there a commandline tool for intersecting Shapefiles?

text: in QGIS vectorlayers can be intersected. That is apparently done by ftools and seems to be integrated into the QGIS-program itself. That is different from many other functions that call GDAL-tools. I would like to intersect esri-shapefiles in a script and without having to load anything into QGIS. Is there a GDAL-tool or some other commandline-suite that will do this? It would be nice to have a tool that handles the shapefile natively (GDAL- tools do that) without having to import it or convert it to some internal format. An alternate solution would be a commandline-interface to ftools. THX stn


title: Converting TAB to MIF / MIF to TAB

text: How I can export (import) data from TAB format to MIF and vice versa?


title: Research papers on SWE (Sensor Web Enablement) and SOS (Sensor Observation Service) with security

text: I am looking for research papers (conference/journals) on the topic of SWE (Sensor Web Enablement) and SOS (Sensor Observation Service) also dealing/including security aspects. I have already found some sources (after cleaning up, will post them here too) but i would appreciate it, if you could post any expert material or websources where to search.


title: How to style the WFS layer in openscales?

text: > **Possible Duplicate:**   >  Styling the WFS point layer (openscales 1.1 & Flex 3) I am very new to openscales... I want to change the styles of the WFS point layer dynamically by giving rules. I want to show some jpeg/iamge as the point in the openscales. Any idea how to define the rules to change the styles of default point style. thanks in advance..


title: Does MySQL has function which will return the points inside polygon?

text: I have MySQL database with polygons, and also the database with points (column longitude and column latitude), is it possible to return all points inside one polygon? Thanks!


title: Openlayers Vector Layers eventListeners: how to specify different data without defining multiple layers

text: I have a map. It has two origins, and two destinations. The code for one origin is here:               vectors_origin = new OpenLayers.Layer.Vector("Origin", {                 eventListeners:{                         'featureselected':function(evt){                         var feature = evt.feature;                         var popup = new OpenLayers.Popup.FramedCloud("popup",                         OpenLayers.LonLat.fromString(feature.geometry.toShortString()),                         null,                         json.orig1_iata_code.iata+"<br>",                         null,                         true,                         null                         );                         popup.autoSize = true;                         popup.maxSize = new OpenLayers.Size(400,800);                         popup.fixedRelativePosition = true;                         feature.popup = popup;                         map.addPopup(popup);                 },                 'featureunselected':function(evt){                 var feature = evt.feature;                 map.removePopup(feature.popup);                 feature.popup.destroy();                 feature.popup = null;                 }         }         });      Note this line: `json.orig1_iata_code.iata+"<br>",` This takes care of one origin. Now since all 4 points (2 origins, 2 destinations), this HTML is going to differ a bit. There will be an `orig1_iata_code`, `orig2_iata_code`, `dest1_iata_code`, and `dest2_iata_code` I really don't want to have to define 4 separate vector layers. Sure it will work, but that seems pretty tacky (correct me if I'm wrong). So where that html is called above, `json.orig1_iata_code.iata+"<br>",`, is there any way to specify 'multiple' optional parameters? Multiple I've done before, but all of them will be called every time.


title: ArcGIS 10.1 crashes on use of python window

text: When I open the python window in either ArcMap or ArcCatalog (ArcGIS 10.1) as soon as I press a key it freezes the entire program and it crashes. It doesn't produce any warning, error, or pop up with info afterwards it just crashes. I tried searching for similar problems through the forum but none of their solutions worked. I installed Service Pack 1 and that didn't solve the problem. I can use IDLE and/or PythonWin to import Arcpy but some custom python tools won’t work. Some of those tools crash it in the same way. Has anyone experienced anything like this or have any idea on how I can get python window to work for Arc 10.1? (I have also tried uninstalling and re-installing ArcGIS 10.1)


title: Spatial Analysis - Calculate a cells elevation (height) above the nearest stream cell

text: I have searched for some ideas on how to formulate this analysis, but have come up short. I imagine there is a pretty straight forward answer; any help would be greatly appreciates. Data: DEM, 1/9th arc second (~10-meter), 239,891,470 cells, projected; Streams, National Hydrology Data set, high resolutions, 9,470 features vector; Streams, Same as above, but converted to a raster Problem: Need to create a raster of the same resolution and extent as the DEM where each cells's value is its height (in meters) above the nearest stream (line or cell). I am working with ArcGIS 10.1 with spatial analyst. I would prefer to run this analysis in Python, but can do it in R if it is easier. Edit to clarify: each cell in my stream raster has the value of the elevation at that point. All cells that are not within a stream have a value of zero. For every cell in the extent of the DEM, I am looking to find the nearest stream cell and subtract the value of that cell (the elevation) from the elevation of the cell being analyzed. Calculate the difference in elevation between each cell and the nearest cell representing a stream. Thank you for your consideration!


title: Why does ArcGIS 10.1 Upgrade Geodatabase tool throw ERROR 999999: Error executing function?

text: I am trying to upgrade File Geodatabases using a Python script that runs fine until the step where it runs the **Upgrade Geodatabase** tool. The error is: > "ERROR 999999: Error executing function.\nFailed to execute" The error can also be seen when the **Upgrade Geodatabase** tool is run manually so I have not re-included the code which can be found in an **earlier question** which led to this one. Does anyone know what would cause this error?


title: Checking if an attribute matches a lookup table

text: I have a point file representing plants (Species 1, Species 2, etc.) and a polygon file of ecozones (Zone A, Zone B, etc.). I also have a lookup table for which zones each plant species is known to grow in, which looks like: Species A, Zone 1 Species A, Zone 2 Species A, Zone 3 Species B, Zone 2 Species B, Zone 3 Species B, Zone 4 etc. My goal is to determine whether or not each plant point is inside a zone it is known to grow in. By using a spatial join, I am able to apply the zone code to each point, but how can I relate this value to the lookup table to determine whether the zone assigned is in the lookup table? Is there something I should be doing before the spatial join?


title: Loading Raster and vector data In MapGuide using Maestro

text: I am new to Mapguide , working on it for my summer Internship.I have Installed it perfectly. Now the task is to load vector data and Raster data in map Guide (mainly have to load *.shp files and *.tif files). I have find not much resource that provide exact steps. Also I don't know where to put data files ? In Repository or somewhere else ? How to link these file ? Can any one help me in this?


title: Communication Between Browser and QGIS Server/QGIS API

text: Suppose I have QGIS server up and running with some vector data. Now a user accesses this data from a webpage. Would it be possible in some ways for the user to send a request to the QGIS Server and make it do some Vector Geo- processing at the server side and return the results? (For example user sends a select query on a point layer. The Server will calculate the Convex Hull and return the WKT.) As per the QGIS website, > QGIS Server provides a web map service (WMS) using the same libraries as the > Quantum GIS (QGIS) desktop application. If it is not possible at present, would it be possible to create something like a REST service(i got the term from ArcGIS REST API) which can hold communication between user and QGIS server? What are the difficulties in implementing them?


title: ArcGIS Autoincrement field for file geodatabase using domains

text: I have read this thread "Autoincrement field when New Feature added" which suggests a python script to autoincrement a field. Can anyone say if there is anyway to do this using domains or any other simpler way in a file geodatabase? I have to edit streets and would like to define a field which simply counts upwards. For ArcGIS users that are not well versed in Python, scripts are pretty daunting, and having to run a tool every time an editing session is finished is also pretty inefficient. It surely can´t be toooo difficult to have a new sequence-like data type in ArcGIS, can it?


title: How to use v.to.3d in GRASS GIS?

text: I am trying to use GRASS GIS for the first time. After quite a lot of trouble I managed to create a location and a new mapset, and my shapefile is now displayed in the window "GRASS GIS Map Display" (for some reason only in 3d view). My next step would be to use the v.to.3d module to add a value to the z-coordinate. I don't seem to get it right, first I tried to use "Search module", and when I found v.to.3d I doubleclicked and filled out the appearing form. That didn't work... Then I tried writing it manually in the "Command console", but being a total GRASS GIS-beginner I'm not sure how to write? I get, or I think I get, that I need the command 'v.to.3d', an input, output and in some way add a height value. Can someone please give me a hand? Exactly how can I do this?


title: Convert GML to GeoJSON

text: I have downloaded a GML file from Statistics Canada and would like to convert it into GeoJSON to use parts of it in a D3 visualization on the web. However, I did not find a converter that would do that for me (I tried http://converter.mygeodata.eu/ , but the GeoJSON it returned did not contain any actual geometry). What is the best approach to convert GML into GeoJSON?


title: How can i display markers from a database?

text: I'm an absolute GIS beginner and here is what i am trying to achieve. I am developing a web application that allows users to :   1. type an address, and be displayed the corresponding slipping map   2. select a point on this map, and create a marker with custom info   3. this marker and associated info will be stored in a database   4. every user requesting an address will be presented a slipping map, showing all the markers other users have defined. So far, i have learnt how to get lat/lon from an address, using nominatim, and display a map via openlayers. I think i know how to store markers in a database, but i don't have a clue how to display an openlayers slipping map, with my markers (several thousands) applied on it ? Is it possible without installing my own mapserver/geoserver ? Thanks for your hints. Christine


title: Is there a way to access ps.map in python?

text: i am wondering whether its possible using python interface to load layers and edit them in Print/Cartographic Composer without displaying, and finally export to PDF? something that can be done via ArcPy (ArcGIS python module). thanks for any inputs!


title: Better way to handle HDF and GeoTIFF files

text: I am using GDAL and MODIS reprojection tool in my python project and want to know better way for math new indices of HDF data (MODIS data). Also I want to visualize them in web. Which format is more flexible for math and visualize my results?


title: How to add new raster store and publish for ECW fromat in GeoServer 2.4-RC2?

text: I am using GeoServer 2.4 RC2 and Tomcat7 on windows xp OS. I have extracted the gdal-data.zip (2.4.0) and geoserver-2.4-SNAPSHOT-gdal-plugin into Tomcat 7.0\webapps\geoserver\WEB-INF\lib. Now Geoserver shows ECW in raster data source list. ![enter image description here](http://i.stack.imgur.com/uq7i0.png) Then I did the following; ![enter image description here](http://i.stack.imgur.com/p7AdF.png) My ECW file path is G:\ECW\DB010_7500.ECW. ![enter image description here](http://i.stack.imgur.com/EUN2P.png) I gave Data Source Name as ECW (folder name). Then I clicked the save button but nothing happened. I navigated to Layers->Add a new resource, I can see site: ECW in dropdown list. But there is no effect when I select it from dropdown list. ![enter image description here](http://i.stack.imgur.com/B3RPd.png) So I am not able to publish it. I don't know these steps are right way to publish ECW image on Geoserver. Please suggest me to do this.


title: How to add coded value to domain?

text: I wish to add a new value in my domain and I get this message : An invalid SQL statement was used [UPDATE GDB_Domains SET DomainName = 'Statut_Pygargue',Description = Statut de l'activite du nid', MergePolicy = 3, SplitPolicy = 3 WHEREID = 20] I am the owner of the file geodatabase.


title: How do I generate a correlogram from a GeoTiff using R?

text: I'm trying to evaluate the scale of spatial autocorrelation for a bunch of rasters I have as GeoTIFFs. I'd like to be able to do this entirely using R, in order to maintain smooth workflow. What's the best way of going from GeoTIFF to correlogram using R?


title: Applying SBAS corrections properly

text: The SBAS standard, RTCA-MOPS, contains lots of unclear information. I am trying to implement SBAS corrections by the help of EGNOS SDK which is an open source software by ESA. When I apply long corrections, results get worse. Fast and ionospheric corrections make it better or worse. Can you suggest any more detailed documents and source code?


title: How to convert contour lines to LAS points?

text: i have downloaded free DEM form the NET and i generated a contour line with 10 meter for interval, so i just want to convert the contour lines to LAS point file to create a new Terrain dataset,and is this a correct way to create terrain dataset or not, if not what is the best way to create terrain dataset with a free data.


title: Reading environment variables in arcpy

text: I have a script that tries to read from the arc environment using               import arcpy     tolerance = arcpy.env.XYTolerance      as recommended here: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//002z0000000w000000 Alas `tolerance` always ends up as `None`. Even when I set it explicity in the `Environments...` part of the script dialog box. Is this a bug, or am I barking up the wrong tree?


title: ArcPy: Using lblclass expression to change font style and color

text: I´m trying to label features of a dataset using ArcPy and the **lblclass** function. While commands like "VBNewline" and "Round" work pretty well, I´m not able to change colors or font style. An expression that works in ArcMap doesn´t work in the script. Could anyone help me out please? I´m using ArcGIS 10.0 Code that works:                #Adding Labels         layer = arcpy.mapping.ListLayers(mxd, "")[0]          if layer.supports("LABELCLASSES"):             for lblclass in layer.labelClasses:                 lblclass.className = "Flaeche"                 lblclass.expression = "[FID] & VBNewLine & Round([AREA_mm2],2)"             lblclass.showClassLabels = True         layer.showLabels = True         arcpy.RefreshActiveView()       Code that works in ArcMap but not in ArcPy               #Adding Labels     layer = arcpy.mapping.ListLayers(mxd, "")[0]      if layer.supports("LABELCLASSES"):         for lblclass in layer.labelClasses:             lblclass.className = "Flaeche"             lblclass.expression = ""<CLR red='255' green='0' blue='0'>" & [FID] & VBNewLine & Round([AREA_mm2],2) & "</CLR>""         lblclass.showClassLabels = True     layer.showLabels = True     arcpy.RefreshActiveView()


title: How to equalise colour, brightness and contrast across several scanned paper maps

text: I am a vectors guy mostly, but for a current project I need to work with a bunch of old scanned paper maps (ww2 bomb damage maps for London, if your interested!) We scanned and georeferenced the maps, and now want to create a tiled composite layer for serving on a website. We will cut out the borders obviously, thats not the issue here. Problem is, that there are unsightly visual colour and brightness differences between map sheets. I am a bit at a loss as to how to equalise them so to give a coherent visual look. I read up on histogram equalising, but my current toolbox (Manifold GIS, GDAL, GeoServer) doesnt seem to have the required functionality. Example of 4 scans already georeferenced: ![Example of 4 scans already georeferenced](http://i.stack.imgur.com/JesBO.jpg)


title: How to fix QGIS error "Entry Point could not be located"?

text: Today, after returning from vacations, I started my QGIS 1.7 and got this terrible message: ![enter image description here](http://i.stack.imgur.com/8GHei.png) Help, anyone?


title: Launching QGis 1.8.0 Lisboa on Windows 8 getting an error "the procedure entry point could not be located in the dynamic link library" [Solved]

text: > **Possible Duplicate:**   >  How to fix QGIS error “Entry Point could not be located”? I have a problem with QGIS 1.8.0 Lisboa installed on a desktop with Windows 8 Pro. When I launch qgis.bat the message "The procedure entry point ?minimum@QAbstractSlider@@QBEHXZ could not be located in the dynamic link library C:\Progra~2\Quantum GIS Lisboa\bin\phonon4.dll" appears. I would like to precise that when I installed QGIS 1.8 for the first time it ran correctly. I think that some Windows update or new software installations caused the problem. If I try to re-install QGIS nothing changes. Is there a possibility to restore the situation? Thank you, have a nice day


title: libmysql.dll-procedure entry point error

text: > **Possible Duplicate:**   >  How to fix QGIS error "Entry Point could not be located"? I am working with qgis,pyqt4 and python2.5. I was getting error qtcore4.dll could not be located.so i set PATH variable as C:\Python25\Lib\site- packages\PyQt4 so i get rid off it. But now m getting the error: procedure entry error point free_defaults could not be located in the dynamic link library LIBMYSQL.dll. what path to set? is it of mysql?


title: Running QGIS on PC using XT (Desk top)

text: > **Possible Duplicate:**   >  How to fix QGIS error "Entry Point could not be located"? The error shown below occurs each time I try to run QGIS: > The procedure entry point gtCPLError could not bee located in the dynamic > link library geotiff.dll Incidentally I can run QGIS on my note book


title: How to fix Qtcore4.dll can't find procedure entry point?

text: > **Possible Duplicate:**   >  How to fix QGIS error "Entry Point could not be located"? I am working with QGIS, python2.5 and pyqt4 for py2.5. My system variables are as follows:               path:=C:\Program Files\Quantum GIS Wroclaw\apps\qgis\bin;C:\Python25;C:\Program Files\Quantum GIS Wroclaw\apps\qgis;          pythonpath:=C:\Program Files\Quantum GIS Wroclaw\apps\qgis\python;      But when i try to run a .py file through the command shell it gives me this error:               python.exe-entry point not found.      The procedure entry point: ??0QDataStream@@QAE@PAVQByteArray@@H@Z could not be located in the dynamic link library QtCore4.dll      Qtcore4.dll is present in Qgis folder, why does the error occur? and it is not present in system 32 so what do i do then? it was present then i deleted it it still says that.


title: Why can I not start Quantum GIS?

text: > **Possible Duplicate:**   >  How to fix QGIS error "Entry Point could not be located"? Why can I not start Quantum GIS? What do I do about this following error message? "The procedure entry point TIFFMergeFieldInfo could not be located in the dynamic link librarylibtiff.dll".


title: What are arc's guarantees on data precision, XYTolerance (i) when displaying (ii) in arcpy?

text: Interesting behaviour: Two line segments in a dataset have almost, _but not quite_ , identical endpoints. The difference between them is less than the XYTolerance for the data set. Reading the coordinates into arcpy, they come out as different to one another. But if I zoom in and in and in, far beyond the XYTolerance, ArcMap still displays them as identical - no gap appears between the lines. Is this to be expected? Does arc make any guarantees on precision of data drawn on screen, or data passed to arcpy? What will happen if the points are separated by the XYTolerance exactly? Update: I tried to test the boundary case, but the plot thickens. The following script creates a feature class that displays correctly in ArcMap despite having features smaller than the cluster tolerance. Still I have other data that doesn't display correctly under the same circumstances.               import arcpy          out_polyline_feature_class = arcpy.GetParameterAsText(0)          # t is tolerance     t=0.001          ic = arcpy.InsertCursor(out_polyline_feature_class)          edges = [ [ (0,0), (0,t/2) ],               [ (0,0), (t,0) ],               [ (0,0), (-t,-t) ],               [ (0,0), (-2*t,0) ]             ]          line = arcpy.Array()     pnt = arcpy.Point()     for edge in edges:         line.removeAll()         for point in edge:             pnt.X, pnt.Y = point             line.add(pnt)              row = ic.newRow()              row.shape = line         ic.insertRow(row)         del row          del ic


title: Do postgres 9.1, postgis 2.0 from stackbuilder and pgrouting 1.0.3. work together?

text: I am new to pgrouting and I would like to check if the version I have installed works fine with postgres and postgis. I have installed postgres 9.1, postgis 2.0 from stackbuilder and pgrouting 1.03. Any suggestion how I can check if pgrouting works fine?


title: Is WGS84 itself a Coordinate Reference System?

text: I was wondering the difference between **WGS84** and **EPSG:4326** (see this). Can then WGS84 on its own be defined as a _Coordinate Reference System_? Or is it just a _datum_ so that a _Coordinate System_ must be attached to it to fully define a _CRS_? (see e.g. this document, figure below Sec. 6.4.1) Thanks for helping.


title: extract shapefile attributes from raster colors

text: I have two layers:   * a shapefile with states and counties (shp)   * a raster file (with colors and some text and borders, but not a perfect map) (tiff)  I want to use the colors in the raster file to associate values to one of the attributes of the shapefile. I've started to do so manually with qgis. I overlayed the shapefile and the raster, and I read the color and fill in the value in the attribute editor: ![enter image description here](http://i.stack.imgur.com/iinQ5.jpg) I was wondering if this can be done programmatically. The pseudocode would be something like:               for county in polygons_in_shapefile:         center_of_polygon = county.center()         color = raster_color(center_of_polygon.x, center_of_polygon.y)         if color == blue:            county.value = 1         if color == red:            county.value = 2         ... etc ...      I'm a total noob when it comes to GIS, so maybe it's very obvious (and the code may be totally the wrong strategy). Any ideas, tutorials or links would be appreciated.


title: How to Dock ArcMap's Table of Contents

text: Question.. I have ArcView 10.1 and am, without luck, searching for/how to "dock the TOC" on the L side of the screen. When I grab the TOC, I see the buttons light up with arrows for up/down, L/R etc... How do I dock the TOC??? Thanks


title: Can I also import road type using osm2po?

text: I have imported the entire UK OSM dataset using the excellent osm2po. The metatdata of the imported table is shown below and my question is whether any of the columns map to ROAD TYPE. If not may I instruct osm2po to include road type as one additional column? Of course, I wish to do this in order to help estimate accurate costs for different types of journey. Many thanks,                       Table "public.result_2po_4pgr"         Column     |       Type        | Modifiers      ---------------+-------------------+-----------      id            | integer           | not null      osm_id        | bigint            |       osm_name      | character varying |       osm_source_id | bigint            |       osm_target_id | bigint            |       clazz         | integer           |       flags         | integer           |       source        | integer           |       target        | integer           |       km            | double precision  |       kmh           | integer           |       cost          | double precision  |       reverse_cost  | double precision  |       x1            | double precision  |       y1            | double precision  |       x2            | double precision  |       y2            | double precision  |       geom_way      | geometry          |


title: How can I import georeferenced photos (jpg, kml, kmz) using QGIS or ArcGIS for Desktop?

text: I've got several georeferenced digital photos as jpg, kml, and kmz. How can I import them into my GIS projects (QGIS or ArcGIS for Desktop) so I can open them as a popup?


title: SLD Editor with GUI?

text: I am setting up a GeoServer instance to serve up a few different PostGIS databases. I have GeoServer up and running however, it takes me forever to enter the SLD into GeoServer by hand. Is there a better way to generate SLD for GeoServer via some GUI based application? Thanks, Tyler


title: How to select island polygons witin the same polygon layer?

text: i have a dataset containing many polygons. In this dataset some polygons have other polygons (as a donut) inside of them. So those island polgons are contained by only one polygon but it originated from the same layer and has the same attributes. I want to select and merge them but i dont know how to select/find them. It is possible to select polygons based on the number of donuts inside of them? ![enter image description here](http://i.stack.imgur.com/bS4Gk.png)


title: Problem registering database to ArcGIS for server 10.1

text: I want to register an SDE geodatabase to ArcGIS for server 10.1 . I use SQl server 2008 R2 , ArcGIS 10.1 . I backup my database from server "a" and restor to server "B" then I used registering database in ArcGIS publisher .I successfully registered database (fig1), but when i want to publish my map ,the warning raised **" datasource is not registerd with the server and data will be copied to the server"** (fig2) . How do i get rid of the warning ? ![fig1](http://i.stack.imgur.com/z6bNZ.png) ![fig2](http://i.stack.imgur.com/txUvr.png)


title: How to specify one layer among many with getWMSInfo?

text: I am using `getWMSFeatureInfo` with multiple point layers. When I click on a point on map, how to specify which layer that point belongs to?


title: How to handle "Mount: Access Denied" when installing QGIS on Win7?

text: I tried to install qgis on my laptop with Windows 7 Professional and it is a 64 bit machine. after going through the procedure described below I got the message "Mount Access Denied". What is going on? Any help would be greatly appreciated. Thank you. I downloaded osgeo4w-setup clicked on install using "Advanced Install"; Selected "Install from local directory"; Selected "C:\OSGEO4W" as the root directory for the install; Expanded "Desktop" Selected "qgis" Selected "1.8.0-5" Expanded "Libs" Selected "GDAL" Selected "1.9.2-2" Clicked "NEXT" Message: "Mount: Access Denied"


title: update line Exception(openlayers+geoserver+oracle )

text: I have table in Oracle. I publish with geoserver. When I wft-t this layer the insert operation is ok, but update Exception I don't know what's wrong. Please help me. table is : create table ROAD ( ID NUMBER, TYPE NUMBER, SUBTYPE NUMBER, NAME VARCHAR2(50), GEOM SDO_GEOMETRY ) exception : `2013-03-28 08:42:35,733 INFO [geoserver.wfs] - Request: getServiceInfo 2013-03-28 08:42:35,863 ERROR [geoserver.ows] - org.geoserver.wfs.WFSException: java.lang.StringIndexOutOfBoundsException: String index out of range: -1 at org.geoserver.wfs.Transaction.transaction(Transaction.java:110) at org.geoserver.wfs.DefaultWebFeatureService.transaction(DefaultWebFeatureService.java:170) at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method) at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:39) at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:25) at java.lang.reflect.Method.invoke(Method.java:597) at org.springframework.aop.support.AopUtils.invokeJoinpointUsingReflection(AopUtils.java:318) at org.springframework.aop.framework.ReflectiveMethodInvocation.invokeJoinpoint(ReflectiveMethodInvocation.java:183) at org.springframework.aop.framework.ReflectiveMethodInvocation.proceed(ReflectiveMethodInvocation.java:150) at org.geoserver.ows.util.RequestObjectLogger.invoke(RequestObjectLogger.java:54) at org.springframework.aop.framework.ReflectiveMethodInvocation.proceed(ReflectiveMethodInvocation.java:172) at org.springframework.aop.framework.JdkDynamicAopProxy.invoke(JdkDynamicAopProxy.java:202) at $Proxy25.transaction(Unknown Source) at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method) at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:39) at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:25) at java.lang.reflect.Method.invoke(Method.java:597) at org.geoserver.ows.Dispatcher.execute(Dispatcher.java:772) at org.geoserver.ows.Dispatcher.handleRequestInternal(Dispatcher.java:272) at org.springframework.web.servlet.mvc.AbstractController.handleRequest(AbstractController.java:153) at org.springframework.web.servlet.mvc.SimpleControllerHandlerAdapter.handle(SimpleControllerHandlerAdapter.java:48) at org.springframework.web.servlet.DispatcherServlet.doDispatch(DispatcherServlet.java:923) at org.springframework.web.servlet.DispatcherServlet.doService(DispatcherServlet.java:852) at org.springframework.web.servlet.FrameworkServlet.processRequest(FrameworkServlet.java:882) at org.springframework.web.servlet.FrameworkServlet.doPost(FrameworkServlet.java:789) at javax.servlet.http.HttpServlet.service(HttpServlet.java:727) at javax.servlet.http.HttpServlet.service(HttpServlet.java:820) at org.mortbay.jetty.servlet.ServletHolder.handle(ServletHolder.java:487) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1093) at org.geoserver.filters.ThreadLocalsCleanupFilter.doFilter(ThreadLocalsCleanupFilter.java:23) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:74) at org.geoserver.filters.SpringDelegatingFilter.doFilter(SpringDelegatingFilter.java:45) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.platform.AdvancedDispatchFilter.doFilter(AdvancedDispatchFilter.java:49) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.vfny.geoserver.filters.SetCharacterEncodingFilter.doFilter(SetCharacterEncodingFilter.java:109) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:311) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.invoke(FilterSecurityInterceptor.java:116) at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.doFilter(FilterSecurityInterceptor.java:83) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.access.ExceptionTranslationFilter.doFilter(ExceptionTranslationFilter.java:113) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerAnonymousAuthenticationFilter.doFilter(GeoServerAnonymousAuthenticationFilter.java:53) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.authentication.www.BasicAuthenticationFilter.doFilter(BasicAuthenticationFilter.java:150) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.geoserver.security.filter.GeoServerBasicAuthenticationFilter.doFilter(GeoServerBasicAuthenticationFilter.java:82) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.context.SecurityContextPersistenceFilter.doFilter(SecurityContextPersistenceFilter.java:87) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.springframework.security.web.FilterChainProxy.doFilter(FilterChainProxy.java:173) at org.geoserver.security.GeoServerSecurityFilterChainProxy.doFilter(GeoServerSecurityFilterChainProxy.java:97) at org.springframework.web.filter.DelegatingFilterProxy.invokeDelegate(DelegatingFilterProxy.java:346) at org.springframework.web.filter.DelegatingFilterProxy.doFilter(DelegatingFilterProxy.java:259) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.filters.LoggingFilter.doFilter(LoggingFilter.java:71) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.filters.GZIPFilter.doFilter(GZIPFilter.java:41) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.filters.SessionDebugFilter.doFilter(SessionDebugFilter.java:46) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.mortbay.jetty.servlet.ServletHandler.handle(ServletHandler.java:360) at org.mortbay.jetty.security.SecurityHandler.handle(SecurityHandler.java:216) at org.mortbay.jetty.servlet.SessionHandler.handle(SessionHandler.java:181) at org.mortbay.jetty.handler.ContextHandler.handle(ContextHandler.java:726) at org.mortbay.jetty.webapp.WebAppContext.handle(WebAppContext.java:405) at org.mortbay.jetty.handler.ContextHandlerCollection.handle(ContextHandlerCollection.java:206) at org.mortbay.jetty.handler.HandlerCollection.handle(HandlerCollection.java:114) at org.mortbay.jetty.handler.HandlerWrapper.handle(HandlerWrapper.java:152) at org.mortbay.jetty.Server.handle(Server.java:324) at org.mortbay.jetty.HttpConnection.handleRequest(HttpConnection.java:505) at org.mortbay.jetty.HttpConnection$RequestHandler.content(HttpConnection.java:842) at org.mortbay.jetty.HttpParser.parseNext(HttpParser.java:648) at org.mortbay.jetty.HttpParser.parseAvailable(HttpParser.java:211) at org.mortbay.jetty.HttpConnection.handle(HttpConnection.java:380) at org.mortbay.io.nio.SelectChannelEndPoint.run(SelectChannelEndPoint.java:395) at org.mortbay.thread.BoundedThreadPool$PoolThread.run(BoundedThreadPool.java:450) Caused by: java.lang.StringIndexOutOfBoundsException: String index out of range: -1 at java.lang.String.substring(String.java:1937) at org.geoserver.wfs.WFSException.init(WFSException.java:103) at org.geoserver.wfs.WFSException.<init>(WFSException.java:23) at org.geoserver.wfs.UpdateElementHandler.checkValidity(UpdateElementHandler.java:127) at org.geoserver.wfs.Transaction.execute(Transaction.java:225) at org.geoserver.wfs.Transaction.transaction(Transaction.java:104) ... 89 more `


title: How do you parse an HTML formatted attribute value to seperate columns in a database table?

text:
All, I have an attribute field formatted using HTML all ready to go for google earth markers. The KML in question is here: http://dl.dropbox.com/u/95231813/state_pts_MI.kml Here is the value:
         Facility type:

         Geoid: 261450002001020

         County: Saginaw County

         2010 Correctional Population: 18
      ![enter image description here](http://i.stack.imgur.com/KJk6T.png) (i tried to paste the actual HTML but couldn't figure out how to display it as text rather than formatted html!!!) I've love to have seperate attribute fields for the columns within that string: Facility Type   County   2010 Correctional Population   How can this be done? QGIS? LibreOffice? Python? Thanks!!! -mb


title: How to make a spatial view in PostGIS and add it as a layer in QGIS?

text: In PostGIS 2.0, what is the most efficient way to create a spatially enabled view that is viewable in QGIS? I followed the context of Part 2 of @Mike Toews answer to this question, and the view shows up as a point layer (in my case) in the "Add PostGIS Layers" dialog. But after creating the view and spatially enabling it, I received the following error when trying to add the view to QGIS: > You must select a table in order to add a layer Can someone please explain why the spatial view will not load into QGIS, and what I need to do in order to make it loadable as a layer in QGIS? **UPDATE 1:** ( _I'm currently using QGIS v1.8 from the OSGEO installer_ ) I view is built from the geonames database, loaded into PostGIS. The geoname table/spatial data can be selected/viewed in QGIS just fine. When I try to select the " **mt_view** ", QGIS sees the view in the dialog, but won't let me select it. It gives a dotted box around the area, but does not highlight and select the view. ![QGIS cannot select the mt_view as a layer](http://i.stack.imgur.com/r3xJd.png) **UPDATE 2:** I just noticed that when I run my first query:               CREATE OR REPLACE VIEW data.test AS     SELECT * FROM data.geoname     WHERE admin1='MT'      the view shows up in the QGIS Add Layer Dialog, but is unselectable. Then when I run:               SELECT Populate_Geometry_Columns('data.test'::regclass);      I get the result of "0", where according to the documentation, I think it should be returning "1". Any thoughts?


title: Is there a way to show pop-up infowindow for every new point in torque map without clicking or hovering on it?

text: Trying to work with torque map, but I want to show a small story behind every point... Thank you for your help


title: What kind of map is this

text: I would like to have/create a map of Kazakhstan in the same way as in the attached picture (grey map of UK and Irland) using QGIS? I don't know what kind of map is this, what its name and where I can get it? And how to have it in this small size with a point on it? Thanks in advance. ![Grey map](http://i.stack.imgur.com/M5CnY.jpg)


title: Postgis shapefile loader 2.1 error: dbf file can not be opened

text: The .shp, .dbf and the rest of files are fine, I open them in QGIS, AGOL with no problem. In fact, I've tried 3 different layers, and get always the same error. Anyone can please help, as I really need to migrate my .shp's to postgis, and have no idea of why it's not able to open .dbf? Postgres v9.3 x86 Postgis shapefile loader v2.1 Thanks a lot in advance! PS I've tried changing the default character encoding from the default 'UTF-8' to 'ISO 8859-14' but same error.


title: Converting map coordinates using AutoCad to Excel sheet

text: I have some maps of pieces of lands and I want to convert the coordinates to excel sheet I want to arrange the data like this in the excel sheet: land number( the pink number in the pic) accompanied by its coordinates in each row I will be really thankful if you can help me ![enter image description here](http://i.stack.imgur.com/hBQw3.png) link to the map file above: ## update I want to have 3 columns in the excel sheet 1) representing the land number text value 2) representing the the coordinates of that land separated by "_" like this `23.87718162405882_57.10554858573619_23.877357092058385_57.105595188322354_23.877292855536815_57.10588146117277` each two representing a point/corner 3) the block number where the land is if the format of the coordinates can not be like above then there is no harm having 3 columns x and y or easting and northing plus point number. therefore each land will have rows of the same points number


title: How to combine several datapoints in one location to a single point varying in size accordingly?

text: This might be a very simple thing, but I just can't figure it out how to do it in QGIS (alternatively I'm also grateful for explanations for ArcGIS). I have a fisheries dataset with multiple catches in one location. I'd like to display these catches as a dot varying in size according to the number of catches - like dots displaying population size in a city. However, the excel sheet I have is one catch each row, with the same coordinates but different attributes (length, weight, etc). Is there a tool in QGIS that can merge all points in one location (other attributes not needed then)? Thanks a lot for your help!


title: Output PostGIS geometry data as centroid or polygon depending on zoom level (GeoServer)

text: This is what I want: when I zoom out, cities of the country should appear as points and at a certain level when I zoom in, I want them to appear as polygons (showing the area of the city and boundaries). How can I do that? I am thinking of creating two tables "cities_points" and "cities_polygons" in my PostGIS database. Any idea? (actually not only the cities, but it's the same for other components of my map).


title: Qgis-Web-Client - Windows - Red Screen

text: Well, I'm a Windows user and a newbie with Web mapping technology, but I've had success with qgis server and the "getmap" and "getcapabilities" work well. Now I'm trying to have the same success with Qgis Web Client, but appeared a little problem: The metadata, legend and layers appear, but the map is not showing up (instead appears a red broken layer). What could be the problem? ![windows - Qgis web client](http://i.stack.imgur.com/VKZrM.jpg) I don't know, but maybe I could have done something wrong in files like qgis- web-client.conf, or globaloptions.js, so i'll show you screenshots of it. Qgis-web-client.conf: ![Qgis-web- client.conf](http://i.stack.imgur.com/vQUQF.jpg) GlobalOptions.Js: http://jsfiddle.net/nYfLZ/ ![GlobalOptions.Js](http://i.stack.imgur.com/MRDap.jpg) ![GlobalOptions.Js](http://i.stack.imgur.com/WMSiK.jpg) ![GlobalOptions.Js](http://i.stack.imgur.com/plIcQ.jpg) ![GlobalOptions.Js](http://i.stack.imgur.com/39oOw.jpg) Or, do I need to edit GISProjectListing.js ? http://jsfiddle.net/rWV7v/ GISProjectListing.js Screenshot: ![GISProjectListing.js](http://i.stack.imgur.com/ZKMdk.jpg) Or it could be something else? ![Qgis](http://i.stack.imgur.com/Uu0YF.jpg)


title: Where are TIGER's "Consolidated Cities"?

text: The TIGER reference docs make numerous references to Consolidated Cities. It's also one of the checkboxes on TIGERWeb. But there are only 5 files in the Consolidate cities directory and available for download (http://www2.census.gov/geo/tiger/TIGER2011/CONCITY/). And I can't find any in the Places file (i.e. no lines have MTFCC='G4120'). So can someone please explain how to find which "Consolidated City" (if any) any given (lat,long) belongs too in the US ? Thanks


title: World map picture to globe

text: i have 26 pictures of the globe and i want to bound them. How can i proceed to abtain a flatworld map? it looks circa like this but in 26 pieces. map Can somebody help me?


title: Which is the best Android tool to visualize JPG+JGW images?

text: I've got some georreferenced JPG+JGW images, and would like to visualize them on my Android device. Which is, in your opinion, the best tool? I've tried OSMdroid (after converting JPG+JGW with MapTiler), but I can't create points to mark some places. Are other alternatives? Thank you very much in advance.


title: How to make interpolation with SLD?

text: I want to make color interpolation with SLD on Geoserver. Interpolation works fine. But i need to know min, middle and max property values to make it correct. So i want to know this values, i says:                 <LineSymbolizer>              <Stroke>                <CssParameter name="stroke">                  <ogc:Function name="Interpolate">                    <!-- Property to transform -->                    <ogc:PropertyName>AverageDailyIntensCar</ogc:PropertyName>                         <!-- Mapping curve definition pairs (input, output) -->                    <ogc:Literal>0</ogc:Literal>                    <ogc:Literal>#fefeee</ogc:Literal>                         <!--<ogc:Literal>15000</ogc:Literal>-->                    <ogc:Div>                      <ogc:Function name="max">                         <ogc:Literal>0</ogc:Literal>                         <ogc:PropertyName>AverageDailyIntensCar</ogc:PropertyName>                      </ogc:Function>                      <ogc:Literal>2</ogc:Literal>                    </ogc:Div>                    <ogc:Literal>#00ff00</ogc:Literal>                         <!-- <ogc:Literal>-->                       <ogc:Function name="max">                         <ogc:Literal>0</ogc:Literal>                         <ogc:PropertyName>AverageDailyIntensCar</ogc:PropertyName>                      </ogc:Function>                   <!-- </ogc:Literal>-->                    <ogc:Literal>#ff0000</ogc:Literal>                         <!-- Interpolation method -->                    <ogc:Literal>color</ogc:Literal>                    <!-- Interpolation mode - defaults to linear -->                  </ogc:Function>                </CssParameter>              </Stroke>           </LineSymbolizer>      Have no validation errors, but in layer preview i get all features with same colors.   My main problem is a middle value:                              <ogc:Div>                  <ogc:Function name="max">                     <ogc:Literal>0</ogc:Literal>                     <ogc:PropertyName>AverageDailyIntensCar</ogc:PropertyName>                  </ogc:Function>                  <ogc:Literal>2</ogc:Literal>                </ogc:Div>                <ogc:Literal>#00ff00</ogc:Literal>      Here something wrong becouse when i use Litteral for middle value interpolation works fine and i get features with differend colors.   Whats wrong here?


title: Write WKT/WKB to File Geodatabase using FGDB API

text: Is there a way to set the geometry of a row using either WKT or WKB with the File Geodatabase API? Is there any other way to write spatial data to the file gdb using the fgdb api without having to go through the highly onerous task of defining the number of parts, defining the number of points per part, building point arrays and part collections all before you can finally set the geometry?


title: Is there any efficient way to convert Geojson to WKT?

text: I want to use ST_GeomFromText() to load vector data to PostGIS. However, I don't know how to convert from Geojson type to WKT. N.B. I know I can use OGR2OGR to upload the geojson directly, but I just want to confirm if there is any way to convert it or not. Thanks.


title: Why does AddField give ERROR: 000732 Input Table: does not exist or is not supported?

text: I was handed off a a bunch of state shape files and wrote a script to add a few new fields and calculate some attributes. The first line in my script is -               arcpy.AddField_management('C:\WB_prj\cityPts\ghosts.gdb\vietnam_adm3_citypt_albers_ghost', 'City_cnt', 'SHORT', '#', '#', '#', '#', 'NULLABLE', 'NON_REQUIRED', '#')      I have all of the state shapefiles in a geodatabase. The script runs perfect for some of the files but for others I get an error that is shown below. My pathways are correct and I can manually add the fields through ArcMap but when I try and run the script for some of them I get the error. Any ideas on how to solve this? I'm running ArcGIS 10.0 > File "C:\WB_prj\scripts\ghostRegions_final.py", line 79, in add_fields() > File "C:\WB_prj\scripts\ghostRegions_final.py", line 14, in add_fields > arcpy.AddField_management('C:\WB_prj\cityPts\ghost.gdb\vietnam_adm3_citypt_albers_ghost', > 'City_cnt', 'SHORT', '#', '#', '#', '#', 'NULLABLE', 'NON_REQUIRED', '#') > File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\management.py", > line 2625, in AddField raise e ExecuteError: Failed to execute. Parameters > are not valid. ERROR 000732: Input Table: Dataset > C:\WB_prj\cityPts\ghost.gdb ietnam_adm3_citypt_albers_ghost does not exist > or is not supported Failed to execute (AddField).


title: What is the tiling scheme of TileCache?

text: I am using TileCache from MetaCarta to create tiles for a WMS service. It creates a tile with an URL like http://example.com/tiles/Admin/06/000/000/357/000/000/158.png. I am trying to understand how to start from Map Coordinates, and get the appropriate tile. I am looking for some documentation which explains the Folder structure/hierarchy of the server cache directory


title: Adding table from Postgis

text: While trying to import a table containing raster data in QGIS , Raster column is not being imported ... I am getting following messages:   1. Database connection was successful , but the accessible tables could not be determined.   2. The view 'public.raster_columns' has no column suitable for use as a unique key. Quantum GIS requires that the view has a column that can be used as a unique key. Such a column should be derived from a table column of type int4 and be a primary key, have a unique constraint on it, or be a PostgreSQL oid column. To improve performance the column should also be indexed. **Please help.** Table Used:               CREATE TABLE everest     (       rid serial NOT NULL,       rast raster,       filename text,     )      * * * After adding primary key constraint also , I am getting same error message. I use raster2pgsql command to create a sql file corresponding to the raster tif file. And thereafter pgadmin to create table in postgis from the sql file , Here , rid is primary key and of type serial .


title: Equivalent of Simple Data Loader in Modelbuilder

text: I have a model which processes a user specified polygon feature class `InputFC`, with the final result a dissolved (based on a user specified dissolve field `D_FIELD`) polygon feature class `OutputFC`. The model also creates a file geodatabase and an empty feature class called `FinalFC` using `InputFC` as a template (for the attribute schema). If I were doing this manually, I would now load the data from `OutputFC` into `FinalFC` using the `Simple Data Loader` in ArcCatalog. Since they have the same table structure, ArcCatalog would automatically match `D_FIELD` in `OutputFC` to its matching `D_FIELD` in `FinalFC` and load the data into it. The rest of the attributes in `FinalFC` would be empty. Is this functionality available as a tool for ModelBuilder, or is there a better way to get the dissolved data into the correct table structure?


title: How to draw lines, dashed lines, routes and add images to the QGIS mapcanvas?

text: I need to draw a rectangle in a map canvas and inside it draw routed dashed lines and add small images to it. These actions are to be performed using the QGIS API. Is there any place with examples that I can start with? All that I need is a direction of implementing this.


title: Adding shapefile directly to geoserver

text: I have py script which is creating shapefiles using shapefile.py. Once the shapefile is created I have to add them manually to the geoserver. Is there a way to add it dynamically to the geoserver? through postgis or something? **UPDATE** Ok so seems my project manager has changed his plan, he wants me to add the shapefiles to the PostGIS Database which will update the newly added shapefiles to the geoserver where there will be a connection to the PostGIS Database. Is it possible to update the newly added shapefiles to reflect in the Geoserver?


title: Comparison of WMS and WFS services

text: I am writing my thesis about WMS and WFS services. I would like to compare these two services. Can you suggest factors on which I can do this comparison? (In my mind I have: functionality, fastest data transmission.) Maybe you can also advise me about any good literature where I can find appropriate information?


title: Using Python in Arcmap to partition strings

text: I have a title "Text that remains Proposed Thing to replace" > var_a = New Text I want the output to be > "Text that remains Proposed New text" Currently outputs this > "(u'Text that remains ', u'Proposed', u' Thing to replace')" Trying to add the var_a               b = a.partition('Proposed') + var_a      But get this error > : can only concatenate tuple (not "str") to tuple Failed to execute > (LegendText).               # Loop through each text element in the map document     for textElement in arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT"):          # Replace Title text     if textElement.name == "Title":          # With this text     a = textElement.text     b = a.partition('Proposed')     textElement.text = str(b)      Any ideas?


title: What are my options in terms of implementing a WPS server?

text: I know about PyWPS, are there any good alternatives?


title: ArcGIS: Place all labels on a map non overlapped

text: I was wondering if there is a way to place all labels on a map non overlapping, but by extending an arrow/line from the label to the point. I've tried the call out box, but I couldn't find a way to enable all labels, and second extend the lines


title: How much Water on Earth(USGS) : A Wrong Map representation?

text: I happened to read a USGS article, How much water is there on, in, and above the Earth? showing a representation of all the water(which includes water from Oceans, Seas,Bays,Ice caps, Glaciers, Permanent Snow and ground water) displayed as a small sphere as shown below. ![USGS water on earth representation](http://i.stack.imgur.com/eIvFF.jpg) I find it hard to believe although the site has given good statistics for it. Is this a correct representation or is it a map gone wrong?


title: Topographic tile cache needed to replace OSM street map data in an app

text: I have developed and marketed an application that successfully uses a pre loaded tile cache of Open Street Map .png files for a small city area about 20km x20km to support the presentation of tourist information on the web, iPad and iPhone. It is written in HTML5 and uses MoSync to create the wrapper to deliver it on the mobile devices. This approach works well for cities but I would like to extend the app to work in countryside areas (in the UK initially) where if I could access it, Ordnance Survey mapping 25000,50000 and 250000 would deliver a much better user experience. My specific questions are: How do I get hold of topographic tiles? Can I get tiles that have the same slippy map structure as OSM. Can I use Ordnance Survey as a source for the UK and if so is it a one off license fee? What zoom levels are available? I expect 25000, 50000 and 250000. Is 100000 available? What formula do I use to identify the specific tile from a lat and long pair? Will all the tiles at a specified zoom level be the same delta lat and delta long? If not how do I tile them?


title: How to select features with NULL Geometry from Shapefiles?

text: I have been given a Shapefile that contains records with NULL geometry, along with records with non-NULL geometry. This situation appears to cause problems in QGIS (1.8 and 2.0) when editing the geometries. I have been able to remove the NULL records by using the Marquee Select tool in the Map pane to select all the features in the dataset and saving the selection as a new shapefile, what I would like to know is how to identify and delete the NULL geometry records in situ in QGIS. I have been unable to do this using the 'Select By Expression' functionality in QGIS 2.0 with expressions such as the following:                $area = 0           $geometry IS NULL           geomToWKT(  $geometry ) =  ''      Edit: I have uploaded a sample TAB file containing five records, two of which (ID 3 and 4) have null geometries. The format is a MapInfo TAB file as I am having problems creating a Shapefile that preserves the NULL geometries! (I am not sure of the provenance of the Shapefile that sparked the question to begin with)


title: Select more than one place at a time in Google Earth?

text: I've got a Google Earth .kmz file opened in my copy of Google Earth 6.0.1.2032 running on Mac OS X 10.6.7. There are about 5000 points/"places" in the file and they show up fine. What I would like to do is create folders in Google Earth for those points to categorise them. Unfortunately, I noticed that the Google Earth interface does not let me select more than one place at a time! Is there another solution for me to manage and categorise those 5000+ places in the .kmz file other than drag and dropping them individually into folders in Google Earth?! (I really don't want to have to do that!) Thanks.


title: How can I get the elevation of a point in openlayers?

text: I want to know if it is possible to obtain the elevation of a point with some type of control of others. I am painting a linestring with a series of coordinates but I need to know the elevation of every point in order to print a chart. I know that Google has an API for elevation, but I need to use only openlayers. Note: My baselayer is a google map. thanks!


title: Matplotlib Basemap Resolution with PyQt makes Python Quit Unexpectedly

text: I'm having difficulties using Matplotlibs Basemap in combination with a PyQt4 interface. Basically I create a Pyqt4 Canvas in which I plot a Matplotlib Basemap and some markers. I call the function in               self.button3 = QtGui.QPushButton('plot', self)      self.button3.clicked.connect(lambda: self.plot(geom0))      And pass a geometry to my plot(geom0) function.               def plot(self, geom0):                     x1, y1 = zip(*((geom0.GetX(i), geom0.GetY(i)) for i in range(geom0.GetPointCount())))                     print geom0                          self.m = Basemap(projection='merc', resolution = None, area_thresh = 0.1,llcrnrlon=6.66321, llcrnrlat=50.48198, urcrnrlon=7.54761, urcrnrlat=51.46085)                     self.m.drawlsmask(land_color='lightgreen',ocean_color='aqua',lakes=True)                          x,y = self.m(x1, y1)                     self.m.plot(x, y, 'D-', markersize=4, linewidth=2, color='k', markerfacecolor='b')                          self.canvas.draw()      This works fine and the canvas is filled and a map without any kind of resolution is shown with my Linestring from `geom0` but when I change the resolution to low `'l'` or high `'h'` Python quits unexpectedly. Changing the resolution to low or high takes a while to calculate and I guess that is where the problem lies. Do you have any ideas how I can tell PyQt to hold his horses until the Basemap is calculated? Thanks a bunch //Dalton


title: Importing USGS digital line graphs (OPT format) to QGIS?

text: How do I get USGS digital line graphs downloaded as .opt format into a format I can view on QGIS?


title: Why can't I connect to my new local PostGIS database?

text: I'm trying to connect to local PostGIS instance (PostgreSQL.app on mac) with QGIS 2.0, and it gives general connection errors. Connection with same parameters with other clients (psql, PGAdmin) is ok, with both socket and port, so it can't be server configuration.


title: Qgis 1.8 is not recognizing sdw file

text: I'm running qgis 1.8, and trying to load a sid file into my map. It imports fine, but seems to ignore the sdw file that comes with it, as the raster shows up without any apparent georeferencing. Is there a way to force the info in the sdw file to be applied to the sid? Thank you for any insight.


title: How to document GIS analysis?

text: New to QGIS & GIS-meta While I know I can annotate within a map, what is being used for documenting your work:sources of data, choices of color/per layer with each map made ( thinking along the lines of templateing, so as not have to deconstruct a map after the fact ). laundry list of what makes up the map, font selections, color choices. Terence


title: ArcMap10 using Feature to Line to simplify dual carriageways

text: I need to clean a roadmap shapefile to replace all dual carriageways (DCs) with single lines. The `Feature to Line` tool nearly works, but is destroying some of the road links. My current workflow (since the tool doesn't discriminate between lines of different attribute) is based on a strategy of splitting the shapefile into two - one with the DCs and roads they link to to process, and another with everything rest to be rejoined at the end. Specifically:   1. use a buffer `Select Layer by Location` to identify all roads within 1 meter of a DC    2. create new DC layer from selected features   3. apply `Feature to Line` tool to the DC layer with an approx 25m `XY tolerance` (max distance between carriageway lanes) However I'm finding that the new simplified DC line objects (typically the roads that linked to DCs) do not properly match up with their continuations in the untouched layer. To illustrate: ![enter image description here](http://i.stack.imgur.com/wbJ7C.png) Very grateful for any suggestions on how to approach this so the links with the grey roads are maintained. Thanks in advance.


title: Cannot input new data into FeatureDataForm

text: I am wanting to update the fields on my server using a FeatureDataForm. I have followed the documentation, but it still does not work. It also does not have the buttons in the form that say 'OK' and 'Delete' like the documentation does. Any help will be appreciated. XAML (FeatureDataForm):               <Border x:Name="FeatureDataFormBorder" Visibility="Collapsed"                      HorizontalAlignment="Right" VerticalAlignment="Top" Margin="10,10,10,0" Width="300" Height="300" >                     <Border.Effect>                         <DropShadowEffect Color="Black" Direction="-45" BlurRadius="20" Opacity=".75" />                     </Border.Effect>                     <esri:FeatureDataForm x:Name="MyFeatureDataForm"  Foreground="Black" BorderBrush="{StaticResource PanelGradient}" BorderThickness="12"                                              IsReadOnly="False" LabelPosition="Left" />                 </Border>      XAML(MyMap):               <esri:Map x:Name="MyMap" Background="#FFE3E3E3" WrapAround="True"  MouseRightButtonDown="MyMap_MouseRightButtonDown" Extent="-9834972.92753924,4441899.425293319,-9833977.88119163, 4442762.485358352" Loaded="MyMap_Loaded">                     <esri:Map.Layers>                         <esri:ArcGISTiledMapServiceLayer ID="MyLayer" Url="http://services.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer" />                         <esri:FeatureLayer ID="WaterValves" Url="http://localhost:6080/arcgis/rest/services/CalvertCity_Test_2/MapServer/0" Renderer="{StaticResource MySimpleRenderer}"                                                       OutFields="*" MouseLeftButtonUp="FeatureLayer_MouseLeftButtonUp" MouseLeftButtonDown="FeatureLayer_MouseLeftButtonDown_1" DisableClientCaching="True"  Mode="OnDemand"/>                         <esri:GraphicsLayer ID="CandidateGraphicsLayer">                             <esri:GraphicsLayer.MapTip>                                 <Grid>                                     <Rectangle Stroke="Gray"  RadiusX="10" RadiusY="10" Fill="#77FF0000" Margin="0,0,0,5" >                                         <Rectangle.Effect>                                             <DropShadowEffect/>                                         </Rectangle.Effect>                                     </Rectangle>                                     <Rectangle Fill="#DDFFFFFF" Stroke="DarkGray" RadiusX="5" RadiusY="5" Margin="10,10,10,15" />                                     <StackPanel Orientation="Vertical" HorizontalAlignment="Center" Margin="30,20,30,30">                                         <TextBlock Text="{Binding [Address]}" HorizontalAlignment="Left" Foreground="Black" />                                     </StackPanel>                                 </Grid>                             </esri:GraphicsLayer.MapTip>                         </esri:GraphicsLayer>                         <esri:ArcGISDynamicMapServiceLayer ID="CalvertCity" Url="http://localhost:6080/arcgis/rest/services/CalvertCity_Test_2/MapServer" />                     </esri:Map.Layers>                                 </esri:Map>      C#:               private void MyMap_Loaded(object sender, RoutedEventArgs e)             {                 MyFeatureDataForm.FeatureLayer = MyMap.Layers[1] as ESRI.ArcGIS.Client.FeatureLayer;             }                       private void FeatureLayer_MouseLeftButtonDown_1(object sender, GraphicMouseButtonEventArgs e)             {                 e.Graphic.Selected = !e.Graphic.Selected;                 if (e.Graphic.Selected)                     MyDataGrid.ScrollIntoView(e.Graphic, null);             }      I believe this tells me that the map is editable: ![enter image description here](http://i.stack.imgur.com/95k91.png)


title: Change Tile url in Openlayer-3

text: I am trying to change my Tile Url in Openlayer-3. I made this Fiddle. With the code:               var map;     var icons = [];          var source = new ol.source.OSM();     map = new ol.Map({         layers: [         new ol.layer.Tile({source: source}),         new ol.layer.Vector({                   source: new ol.source.Vector({features:icons})         })],         renderer: "canvas",         target: 'map',         view: new ol.View2D({             zoom: 11         })     });          // Change map source url dynamicly     // This line gives error in console:     // Resource interpreted as Image but transferred with MIME type text/html:     source.setUrl('http://a.tile.openstreetmap.org/');          map.getView().setCenter(transform(5, 52));          function transform(lng, lat) {         return ol.proj.transform([lng, lat], 'EPSG:4326', 'EPSG:3857');     }      The code breaks on the line:               source.setUrl('http://a.tile.openstreetmap.org/');      It gives a message in the console (Chrome: Ctrl + shift + j):               Resource interpreted as Image but transferred with MIME type text/html:      Any solutions?


title: Hardware specifications for QGIS

text: New to QGIS, running on a Mac with adequate RAM & HD space. Looking to future I find no real hardware specifications. I want to find more information, can someone point me towards QGIS forum that might have discussions about rendering speed/time in relation to size of documents/maps being built/organized. Terence


title: Wetlands Symbol for QGIS point pattern fill

text: Is there a marker available anywhere for point pattern fill that resembles the standard wetlands symbol used in USGS topographic and other maps? Any good standard symbol that will be intuitively recognized as wetlands?


title: Computing altitude from contour table - slow ST_Dwithin query

text: I am using PostgreSQL 9.2/PostGIS 2.0. I have a table of house addresses and a table of contour lines which are as follows:               address_list          house_id | geom (point, 4283)     ---------+---------         1    | Point         2    | Point          sa_elevation             gid   | altitude  |geom (MultiLineString, 4283)     ---------+-----------+-----------------         1    |  100m     |          2    |  550m     |       I am trying to compute the altitude for each house in the address_list table (approx 6,000 houses) using a 'nearest neighbour' approach ie. find the nearest contour line to the home and tell me what the altitude is for that contour line. My syntax is included below. I use a left join because I don't know if all houses will have a contour line within 3km (following PostGIS in Action p. 137). Note that 3112 is in metres. I've tested the syntax on 5 houses and it works. But the problem is that if I run it on the whole table and leave it running all night it never seems to finish. I have run vacuum analyse and haves spatial indices on each table. Could someone suggest a faster procedure? Or is there anything else I should do? Many thanks.               SELECT DISTINCT ON(a.house_id)          a.house_id,          e.altitude AS altitude     FROM address_list AS a     LEFT JOIN sa_elevation AS e      ON ST_Dwithin (ST_Transform(a.wkb_geometry, 3112),ST_Transform(e.wkb_geometry, 3112), 3000)     ORDER BY a.house_id, ST_Distance(ST_Transform(a.wkb_geometry,3112), ST_Transform(e.wkb_geometry,3112));      The explain plan is:                  Unique  (cost=232208649.81..232208684.17 rows=6871 width=4109)       ->  Sort  (cost=232208649.81..232208666.99 rows=6871 width=4109)             Sort Key: a.house_id, (st_distance(st_transform(a.wkb_geometry, 3112), st_transform(e.wkb_geometry, 3112)))             ->  Nested Loop Left Join  (cost=0.00..232183921.91 rows=6871 width=4109)                   Join Filter: ((st_transform(a.wkb_geometry, 3112) && st_expand(st_transform(e.wkb_geometry, 3112), 3000::double precision)) AND (st_transform(e.wkb_geometry, 3112) && st_expand(st_transform(a.wkb_geometry, 3112), 3000::double precision)) AND _st_dwithin(st_transform(a.wkb_geometry, 3112), st_transform(e.wkb_geometry, 3112), 3000::double precision))                   ->  Seq Scan on address_list a  (cost=0.00..212.71 rows=6871 width=36)                   ->  Materialize  (cost=0.00..46527.78 rows=42852 width=4073)                         ->  Seq Scan on sa_elevation e  (cost=0.00..24845.52 rows=42852 width=4073)


title: QGIS Invalid object name 'geometry_columns' QODBC3:Unable to execute statement

text: I was using Lisboa 1.8 with SQL Express and tables with geometry. Connections were working with 1.8 and are working with ArcGIS, so I know my tables and geometry are fine. I installed 2.4 Chugiak, mapped my SQL Server, test connection works fine, but now I get [Microsoft][ODBC SQL Server Driver][SQL Server] Invalid object name 'geometry_columns'. QODBC3:Unable to execute statement. Can someone tell me what is going wrong all of a sudden. Thanks.


title: Does Terrain Navigator Pro have Attribute Tables?

text: I learned about mapping with ESRI ArcGIS, I am now trying to use **Terrain Navigator Pro (TNP)** and finding that there are things I can't figure out because they are named differently. Is there a different name for an attribute table in TNP? Or does a similar table even exist in TNP?


title: Change layer border in arcmap through arcobject

text: I wish to change the border of layer in arcmap through arcobject vb.net program .Through my program i have added layer to arcmap now i want to change its property .Can anyone please tell me where to start with?


title: Dissolve large dataset in ArcGIS and QGIS

text: I have more than 80k point data. I have created few buffer from 1 to 5 km of them. Now I want to get the covered area by 1 km to 5 km buffered range. For this I have tried to dissolve those individual buffered file. In ArcGIS 10 it took 6 hours and then replied an error. Also in progress bar it was showing 33% completed. I think this 33% progress is some kind of bug as it stay 6 hours. In QGIS it took 8 hours to complete 61% then stopped. Is there any other way or tools to dissolve those large dataset in a reasonable time?


title: USA MSA/City/County with SImplest Border in JSON

text: I deal with a lot of poloygon-related JSON representing US places regularly. One thing that could make reading JSON with my eyeballs easier is having a short list of US places whose borders have short(er) JSON representation (usually as ESRI Javascript API's concept of Geometry-Rings, but other formats are used as well). As a secondary consideration, the more population within the borders of a US place, the more useful it will be to shortlist and use. With those criteria stated, anybody have a favorite US city/MSA/county that has a short JSON-represented border? Or have an empirical-defined list of such? For example, the US city of Granville, North Dakota has a reasonably short JSON border vector/geometry: `"rings":[[[-11225173.4082,6151922.0798999965],[-11225175.3006,6151872.743600003],[-11225176.4138,6151845.650600001],[-11225176.7477,6151750.1567],[-11225176.8591,6151705.1696000025],[-11225181.7571,6151705.1696000025],[-11225173.6308,6151696.138800003],[-11225167.2856,6151684.0977],[-11225162.6101,6151635.933399998],[-11225159.2706,6151554.991300002],[-11225158.4913,6151536.261],[-11225165.0592,6151483.080600001],[-11225163.612,6151407.993100002],[-11225160.6063,6151248.622400001],[-11225158.4913,6151135.576399997],[-11225156.9328,6151108.652900003],[-11225302.984,6151107.983999997],[-11225372.2247,6151107.6494999975],[-11225476.0858,6151105.810000002],[-11225608.89,6151103.636100002],[-11225753.3828,6151101.2949],[-11225751.2677,6150804.641000003],[-11225750.7111,6150721.365699999],[-11225927.4865,6150721.867299996],[-11225927.4864,6150672.370800003],[-11226098.028,6150672.872400001],[-11226097.9166,6150609.497100003],[-11226513.1384,6150610.667599998],[-11226796.6693,6150612.506899998],[-11226874.4816,6150611.6708],[-11226882.6079,6150611.001999997],[-11226882.6079,6150758.321000002],[-11226882.3853,6151024.036600001],[-11226779.9714,6151024.371100001],[-11226779.5261,6151375.382700004],[-11226761.0471,6151380.566799998],[-11226756.817,6151381.737499997],[-11226663.0859,6151390.099100001],[-11226353.5063,6151393.443700001],[-11226362.078,6151574.223200001],[-11226359.8516,6151909.3693],[-11226336.0292,6151910.0383],[-11226238.6247,6151923.4177],[-11226032.127,6151918.734899998],[-11226014.9838,6151919.069399998],[-11225757.947,6151924.254000001],[-11225626.8126,6151923.752300002],[-11225173.4082,6151922.0798999965]]]`


title: Why does image resolution decrease when I import an image into the georeferencer?

text: I'm using QGIS 1.8.0 on Windows 7. Usually the georeferencer functions very well. However, when I try to import a particular aerial photograph (in either .tif or .jpg format) the resolution of the image is reduced within the georeferencer itself.   ![resolution before and after](http://i.stack.imgur.com/NqaQZ.jpg) This generally does not happen with the other photos. Why would this particular photo lose resolution when imported into the georeferencer? Something about the format of the photo? Thank you for your time!


title: QGIS Metatools error

text: I haven't been able to successfully use the QGIS MetaTools plugin for describing metadata. Has anyone used it before? When I try to create metadata for a shapefile, the plugin creates the xml file but then I get a python error:               An error has occured while executing Python code:          Traceback (most recent call last):     File "C:/Users/daniel/.qgis2/python/plugins\metatools\metatools.py", line 267, in doEdit     standard = MetaInfoStandard.tryDetermineStandard(self.metaProvider)     File "C:/Users/daniel/.qgis2/python/plugins\metatools\standard.py", line 35, in tryDetermineStandard     print "TryDetStd ", text, " Text"     IOError: [Errno 9] Bad file descriptor      I'm using QGIS 2.4.0-Chugiak Chugiak, installed from the OSGeo4W installer with python version 2.7.4. Has anyone seen this error?


title: ArcMap "Could not open the specified file"

text: Has anyone ever encountered this message? I am getting this when I try to open certain shapefiles. I am working on a problem that generates shapefiles from spreadsheets using ThinkGeo's .NET tools. When I try to open ones where I have approximately 142 columns (of mostly string, some float, and a couple datetime fields) I get this message. I don't get anything else besides this message. I can open the file in ArcPad 8.0 and QGIS 2.0.1, but it fails every time with ArcMap. I have run the shapefiles through Andrew Williamson's ShapeChecker, which found no problems. Going through the Check Geometry tool gives me an ERROR 000229 saying it cannot open the file. It doesn't have a lock on it, so there is something else happening. Any tools or suggestions that I should look in to? **Edit:** I am using ArcMap 10.2.2. I can view the files in ArcCatalog, but when I go to preview them, I get a message saying "Open Failed" as the title and "Error opening feature class". I tried loading this up in ArcCatalog 9.3.1 and ArcMap 9.3.1. I am able to see the shapes drawn out, but anytime I try to view the attribute data, it crashes without explanation. This is point data generated from x/y coordinates that were generated from ArcMap originally. What I don't fully understand about these errors is how to tell exactly what it is. I have up to 255 columns of data that I can include. I am able to include all the data if I break up the columns into smaller chunks. It seems like for whatever reason my 'no crash' limit is around 142 columns. As for record count. ArcMap throws an _ERROR 000229 Cannot Open_ when I try to run it on the shapefile. ShapeChecker verified that the row count and feature count match though. I have also checked the field name length as I generate them. None go over 10 characters. I haven't found anything definitive yet as to what causes my errors or how to solve them.


title: Connecting isolines (i.e.points of equal value) in 3D scene using ArcGIS for Desktop?

text: I have a point shp theme with fields Lat, Long, Altitude (m), and Number. The Number field – is a numerical parameter of each point (Lat, Long). The Altitude field – is the parameter of elevation of each point (Lat, Long) above the surface. I made a 3d map in ArcScene of ArcGIS10 where I represented each point as a symbol of different size (according to the “Number” field), and then – elevated each point above the surface (according to the “Altitude” field). Is it possible to connect by isolines (or something like that) all points of the equal value of the “Number” field (i.e. all points of the same size) in this 3d map?


title: Is the filename of a shapefile the same as its Feature Class name?

text: This is a very basic question, but I'm not finding the answer on Google or here... I have a bundle of files I have to submit in an incredibly specific format. Our sub-consultants gave us a set of shapefiles, but the shapefiles have a different file name than the required feature class. All the required attributes are correct. For example, they named the file "boundary.shp" (plus all of it's accompanying files) whereas the feature class needs to be named "gis_boundary_limit". Can you have the feature class name be different than the shape file name? If I just rename the files have I renamed the feature classes? Is there a way in ArcMap to rename a feature class within a shapefile? Forgive how basic this question is, I had been operating on the assumption that shapefiles and feature classes were interchangeable.


title: Print function causes ArcMap to hang

text: ArcMap hangs after clicking on print. I get the dreaded spining globe icon as if it is working, but does nothing. I eventually go to task manager to end the task. Any ideas are greatly appreciated.


title: How to improve positioning of images taken with GPS and Gyro info?

text: I'm sort of 'recreating' google street view with MIT right now. We have a car with a camera on top and we're driving at 10mph down a road. As we drive we're shooting pictures at 30fps and each photo is tagged with a timestamp. I can interpolate between images and GPS route so I can tell you the "estimated" position at which every image was taken. However, we have some GPS errors, especially in downtown areas with multipath, etc. To take a first stab at this, we bought a better GPS unit from Trimble the Aardvark GPS + DR which is a decent GPS unit as well as a Gyro. We also tied in the wheel speed pulse to the unit so it can estimate the position of the car even if GPS drops out for a second or two. Over time though, drift occurs and it can tell us that we're 10-15 meters off the road until we get a decent GPS fix again. We've done some basic metrics for the position, filtering out positions that have a faster than expected speed, given the previous points, and that can give me a decent guess at what parts of the route are bad. That said, is there any way I can improve things further? I was going to do some map matching but I wonder if I should only do that in areas I decide are "bad" or if map matching is good for everything. The real purpose of this is to be able to type in an address and pull up all images of that address. Right now we're close, but not close enough. P.S.: Our current stack is Python, Google Maps and OSM mostly.


title: Qgis problem clipping a raster .asc file, the output is unreadable

text: I am currently using QGIS (1.8.0) to prepare files for use in MaxEnt to predict species distribution in a new territory in Argentina. I need to a) cut the raster files (bioclim, and others) to a single province, and b) create a .asc file from a .shp file (not sure if this is possible in QGIS) My problem is with point 'a'. I have read the various questions on this site regarding raster clipping which were very useful. I have set the CRS for each file (WGS84 EPSG4326), and the Raster/Extractions/Clipper tool using a mask does create a file. However, when I open this file using TextPad (with which I can view the original larger file) the text is in an unknown language/script and does not open in MaxEnt. I was worried the problem my lay in the fact that my .shp mask did not fully fall within my original raster, so I tried clipping by extent, so I clipped only within the extent of the raster, but got the same results. Has anyone any suggestions as to what might be the problem?


title: How to read Geotiff tags and pixel per pixel values using GDAL C#

text: I'm new here and I'm new at programing with GDAL C#. My question is how I can get the value of SRTM GEOTIFF pixel to compare with it's neighbors, and how can I get the values of Geotiff tags.


title: Calculation amount or number of certain pixel-values in a raster

text: It seems pretty simple but until now I could not find a proper solution: In a given raster i want to know the number or ratio of e.g. red pixels. Edit: to be more precise: I have a number a arial photos and want to estimate vegetation cover, so I tried to estimate/count the number or amount of green pixels. If there is an other solution to estimate vegetation cover, i am open to any other way. First Attemps were raster statistics and inspecting histogramms, but the results were not clear to me.


title: Convert multiple TIFFs into one?

text: I have downloaded the Blue Marble dataset, and it consists of eight GeoTIFF images, labeled like this: A1 B1 C1 D1 A2 B2 C2 D2 I would like to use these tiles for GeoServer to be served as WMS. Unfortunately, I don't know how to use multiple images as data source. 1) Can this be somehow done in GeoServer? 2) If not, what tool can be used to merge such images into single huge GeoTIFF file? IrfanView is out of question, it freezes trying to open a single image...


title: Calculating spaces between roads from PostGIS database

text: Is it possible to calculate the spaces between the roads using a select statement in a PostGIS/PostGre Database? I tried using some combinations of the predefined functions as ST_Intersects & ST_Difference but did not succeed yet. I would appreciate every help.


title: Python scripts execution in ArcGIS Runtime WPF app

text: DEAR PEOPLE FROM THE FUTURE: Here's what we've figured out so far… I am facing a quite interesting requirement: > "run python scripts within an ArcGIS Runtime app". checking the concepts in the ArcGIS Resources for ArcGIS Runtime SDK for WPF i can read: > Geoprocessing in the runtime provides you with a framework for answering > your GIS questions. The API works with geoprocessing though geoprocessing > service tasks. These tasks are either hosted online, through an ArcGIS > Server, or can be hosted locally, by providing the local server with a > runtime enabled geoprocessing package (GPK). > > Geoprocessing services can host one or more tasks, which come from a tool, > script, or model, authored in ArcGIS for Desktop and then published to > ArcGIS for Server or packaged as a geoprocessing package but also I can read the following statement > Not all geoprocessing tools or functionality are supported in the runtime > local server. Please see the supported geoprocessing tools list. a few lines after is the next statement: > Python scripts Python is an open-source, general-purpose programming > language used as a scripting language in ArcGIS geoprocessing. Python > scripts can be shared for use with the runtime as geoprocessing packages. > Scripts can be turned into script tools, and then packaged into a > geoprocessing package. The ArcGIS Runtime contains its own imbedded Python I understand that I can run python scripts using the geoprocessing packages, but not all the functionality that I can get from ArcPy library is available, Am I wrong? is there a known workaround that allows to run python scripts? (with out using ArcObjects, just ArcGIS Runtime)


title: How do I explain what GIS is in a business context

text: This is hypothetical so don't ask for details: I've heard of this GIS thing with computers and maps and stuff and it sounds like it might apply to this project I'm working on. Could you explain what it is in enough detail for me to decide if it's worth me spending time/money looking in to? In short, I'm asking for a (long) elevator pitch on what GIS offers. * * * While the answers so far have said a few things that GIS does, nobody has given a picture of what GIS is that differentiates it, that makes it more than just a minor specialty in IT (like C++ or Python). As a starting point, _is the following wrong and if so how_? > GIS uses computers to deal with and analyse multiple data sets that are > _interconnected/related by location_ information. > > GIS can help answer questions like: > >   * What is common to these locations? >   * Where do these conditions exist? >


title: Edit Sketch options for digitizing polygons

text: Trying my hand at building a methodology on QGIS as opposed to ESRI, and I'm having some issues with simple digitizing of irregular polygons. The edit sketch graphic that appears when adding a new polygon feature is only transparent in the area where no coverage exists of previously added vertices. The problem arises when trying to digitize a polygon that has internal angles. It gets to a point in the digitizing that you can no longer see your base map to digitize from. I have looked at QGIS digitizing options and how found nothing to adjust the transparency of the edit sketch polygon. so I might be able to see the underlying map as I digitize. Screenshots below to illustrate my point. Is there a way to combat this, other than to complete the polygon and then go in and modify individual nodes to position correctly? ![Boundary to be digitized](http://i.stack.imgur.com/NSe0J.png) ![Overlap issue and inability to see inner corner](http://i.stack.imgur.com/Bk9IG.png) [SOLVED] - By Jake below. I missed the simple ability to modify the Alpha channel of the rubberband setting.


title: Display labels for all point features in ArcGIS

text: I might be missing something silly, but is there a way to guarantee that each point label in my point layer gets displayed?


title: How to best determine the future location of a Fire Station?

text: I have been tasked with creating a map and doing some analysis that shows the future growth (10 years) of our county. They would like to see where the future growth is likely to occur and also show where fire station(s) should be located to support the new growth. My first thought is to look at the census. They give me the 2010 numbers and I can also get a current estimate, however, finding anything beyond that is tricky. In short, I'm not exactly sure the best methodology to use for this analysis. Can anybody give me some pointers?


title: Can't connect QGIS 2.0 to postgresql 9.3 remote server with postgis 2.1

text: I am trying to connect qgis to my postgresql server using the add postgis layer option in qgis dufour. I am not getting any details from qgis other then 'connection failed check settings.' However, when I check the logs from postgresql I am getting that there isn't a function 'select postgis_version' despite the fact that I can use that function within psql and see that it is version 2.1. I enabled the local_addresses=* in postgresql.conf and I'm trying to connect using tcp/ip. So my entry for my local machine in the pg_hba is: hostnossl all all myipaddress 255.255.255.0 trust


title: Can someone recognize the CRS?

text: I'm using QGIS 1.8.0 and have some shapefiles with a "User Defined CRS" which I don't recognize. Is anyone able to recognize it? Here is an image of the settings I see on QGIS. ![enter image description here](http://i.stack.imgur.com/cqSFm.png) Thanks in advance.