fatal error LNK1104 qgis build

I try to build qgis1.8.0 using Microsoft Visual c++ 2010 in a Windows 7 64bit. all the dependency are installed and the solution qgis1.8.0.sln was generated succesfully. but my problem is the appearence of those when i build the solution: fatal error LNK1104: cannot open file '....\core\RelWithDebInfo\qgis_core.lib fatal error LNK1104: cannot open file '....\providers\grass\RelWithDebInfo\qgisgrass.lib fatal error LNK1104: cannot open file '....\analysis\RelWithDebInfo\qgis_analysis.lib does anybody know how to fix this issue ? thank you


Have Legends been implemented for ESRI's REST API?

ESRI's 10.0 documentation for legend indicates that this url should return a legend: http://map.floridadisaster.org/GIS/rest/services/Facilities/Critical_Facilities/MapServer/legend But it doesn't work. Is there something special I need to do when publishing a mapservice to enable legends? If not, does anyone know of a url that returns a legend?


Why doesn't this simple OpenLayers code work?

The problem is that when I create a new shape and add it to the layer and programmatically select the new shape, I cannot unselect that shape using the mouse unless I either click the shape or click a different shape. How do I get the new shape to be selected by default and unselect it by just clicking on the map?               function test3() {          //Get the WKT and Attributes from Shape where Shape_ID = 26786     var feature26786 = layerUserShapes.getFeaturesByAttribute('Shape_ID', 26786)[0];     var formatter = new OpenLayers.Format.WKT();     var originalWKT = formatter.write(feature26786);     var originalAttributes = feature26786.attributes;          //Unselect all the Features     controlPolygonSelecter.unselectAll();          //delete all the Features     layerUserShapes.removeAllFeatures();          //Build a new Feature with the same WKT and same Atributes, and add it to the Layer     var newGeo = OpenLayers.Geometry.fromWKT(originalWKT);     var polygonFeature = new OpenLayers.Feature.Vector(newGeo, originalAttributes); //build a Feature from Geometry & Attributes     layerUserShapes.addFeatures([polygonFeature]); //add this feature to the map          //Select the new Feature     var NewFeatureToSelect = layerUserShapes.getFeaturesByAttribute('Shape_ID', 26786)[0];     controlPolygonSelecter.select(NewFeatureToSelect);          //Now, try to click on the map to UNselect this new Feature. It won't work. The Feature will remain highlighted/selected unless you actually click on it and then click off. Then, the shape will return to its UNhighlighted status and the event for Unselect will fire.      }


Data Driven Pages - Locator Map Zoom

I need to create a map series that highlights a single tax parcel on a 800 scale map. On the same map, I have a separate data frame that zooms into the highlighted parcel. I'm trying to get it to where DDP will zoom into each individual parcel in the second data frame. However it needs to zoom to the parcel, and adjust scale to always have a tight zoom on the given parcel. The main frame MUST stay at 800 scale, and the smaller should always be a tight zoom on the parcel. I'm not having any trouble on the main frame, but how using DDP can I get it to auto zoom in on the parcel in the smaller data frame? I've attached a map example of What it should look like. ![enter image description here](http://i.stack.imgur.com/R5uvH.png)


Getting started with 'shell scripting loops' in GRASS

I'd like to work through the MODIS-Aqua example **here**. But typically the hardest bit is getting off the ground. What specifically does it mean by 'shell scripting'? is it referring to the command console? Python console? To get started I tried the 'Decompress' example `bzip2 -d A20023532002360.L3m_8D_SST_4.bz2` in both command and python consoles - neither worked - I expect I need to get the file path in there, but not sure of the syntax. The file currently resides in the Project>Mapset folder where I thought the default path would lead - is that the best place for it? The next big hurdle will be creating the loop to process hundreds of hdf files, but i can't get off the mark until I work out some of the very basics. Sorry, this is my first foray into GRASS scripting. Any help would be most appreciated!


Are there any sources of UK building footprint data available?

I noticed that there's another question related to this topic, but the question was related to US data. I'm looking at the possibility of producing an interactive map for a UK based city, which would benefit from an overlay describing the building footprints in the city centre. Does anyone know of a good source for this information? Ideally the source would be free, but commercial suppliers could also be considered. Many thanks in advance.


Cached Map Services concept and feature layer

I am having performance issue in my application. My feature layer has 220 polygons, but when I load the map with the feature layer, the browser hangs. On the advice of forum experts, I have now developed cached services, but I don't know what to do with it. As far I know, the cache service loads PNG file. So for each extent or zooming, we have to load different levels of map. Is this right? Secondly, I want to add interactive feature to the map where I want the user to be able to select some countries on the map. So when they click on the map, I can do some javascript action. For this to happen, I have to load the feature layer and then add my onclick functionality. Can I achieve this with Cache service? Do i have to add Feature Layer to the basemap? Please explain. Below is Rest service               Map Name: Layers          Legend          All Layers and Tables          Dynamic Legend          Dynamic All Layers          Layers:              WorldCountries_50K (0)         WorldCountries_100K (1)         WorldCountries_250K (2)         WorldCountries_500K (3)          Description:          Copyright Text:          Spatial Reference: 102100  (3857)               Single Fused Map Cache: true          Tile Info:              Height: 256         Width: 256         DPI: 96         Levels of Detail: 20             Level ID: 0 [ Start Tile, End Tile ]                 Resolution: 156543.03392800014                 Scale: 5.91657527591555E8             Level ID: 1 [ Start Tile, End Tile ]                 Resolution: 78271.51696399994                 Scale: 2.95828763795777E8             Level ID: 2 [ Start Tile, End Tile ]                 Resolution: 39135.75848200009                 Scale: 1.47914381897889E8             Level ID: 3 [ Start Tile, End Tile ]                 Resolution: 19567.87924099992                 Scale: 7.3957190948944E7             Level ID: 4 [ Start Tile, End Tile ]                 Resolution: 9783.93962049996                 Scale: 3.6978595474472E7             Level ID: 5 [ Start Tile, End Tile ]                 Resolution: 4891.96981024998                 Scale: 1.8489297737236E7             Level ID: 6 [ Start Tile, End Tile ]                 Resolution: 2445.98490512499                 Scale: 9244648.868618             Level ID: 7 [ Start Tile, End Tile ]                 Resolution: 1222.992452562495                 Scale: 4622324.434309             Level ID: 8 [ Start Tile, End Tile ]                 Resolution: 611.4962262813797                 Scale: 2311162.217155             Level ID: 9 [ Start Tile, End Tile ]                 Resolution: 305.74811314055756                 Scale: 1155581.108577             Level ID: 10 [ Start Tile, End Tile ]                 Resolution: 152.87405657041106                 Scale: 577790.554289             Level ID: 11 [ Start Tile, End Tile ]                 Resolution: 76.43702828507324                 Scale: 288895.277144             Level ID: 12 [ Start Tile, End Tile ]                 Resolution: 38.21851414253662                 Scale: 144447.638572             Level ID: 13 [ Start Tile, End Tile ]                 Resolution: 19.10925707126831                 Scale: 72223.819286             Level ID: 14 [ Start Tile, End Tile ]                 Resolution: 9.554628535634155                 Scale: 36111.909643             Level ID: 15 [ Start Tile, End Tile ]                 Resolution: 4.77731426794937                 Scale: 18055.954822             Level ID: 16 [ Start Tile, End Tile ]                 Resolution: 2.388657133974685                 Scale: 9027.977411             Level ID: 17 [ Start Tile, End Tile ]                 Resolution: 1.1943285668550503                 Scale: 4513.988705             Level ID: 18 [ Start Tile, End Tile ]                 Resolution: 0.5971642835598172                 Scale: 2256.994353             Level ID: 19 [ Start Tile, End Tile ]                 Resolution: 0.29858214164761665                 Scale: 1128.497176         Format: PNG32         Compression Quality: 0.0         Origin: X: -2.0037508342787E7         Y: 2.0037508342787E7         Spatial Reference: 102100  (3857)          Initial Extent:              XMin: -2.6548096424178407E7         YMin: -1.4376721731989654E7         XMax: 2.204125777392E7         YMax: 2.400182614911706E7         Spatial Reference: 102100  (3857)          Full Extent:              XMin: -2.2041257773920003E7         YMin: -3.578420532071707E7         XMax: 2.204125777392E7         YMax: 2.400182614911706E7         Spatial Reference: 102100  (3857)          Units: esriMeters          Supported Image Format Types: PNG32,PNG24,PNG,JPG,DIB,TIFF,EMF,PS,PDF,GIF,SVG,SVGZ,BMP          Document Info:              Title:         Author:         Comments:         Subject:         Category:         Keywords:         AntialiasingMode: Fast         TextAntialiasingMode: Force          Supports Dynamic Layers: true          MaxRecordCount: 1000          MaxImageHeight: 2048          MaxImageWidth: 2048          Supported Query Formats: JSON, AMF          Min Scale: 5.91657527591555E8          Max Scale: 1.8489297737236E7      HTML code Display Different Map Services Based on Zoom Level               <link rel="stylesheet" href="http://js.arcgis.com/3.6/js/dojo/dijit/themes/claro/claro.css">     <link rel="stylesheet" href="http://js.arcgis.com/3.6/js/esri/css/esri.css">


Select Features by Polygon - Not Displaying the Shape Drawn

I am a relative newbie to QGIS 1.8, have looked through the previous answers, but have found nothing that resolves my issue. I have a number of datazones as polygons on a layer. One of the datazones contains two islands that I need to split. With the layer selected and visible, editing turned on, I choose Select Features by Polygon. I click on a point in the sea around the island, move my mouse to another point, click again, a line appears, repat this until I get to the point where I want to put the last dot and right click. What I understand should happen at this point is that a filled shape should appear, but it doesn't, although the islands become selected. Is this a bug or am I missing something obvious? The same is true of all of the vector layers from other shapefiles. Thanks a lot Andy


Poor rendering performance in ArcGIS Server tiled map service and crashes logged

I have a number of map services in ArcGIS Server (10.1, running on 64-bit Windows Server 2008 R2) and only one of them is cached. Caching is configured on-demand as there's no way I can store the full tile set on disk. Also many of the tiles would be empty as my dataset is quite sparse. Features are only shown when zoomed quite far into the map, which means the service can go some time without being hit (requests are made from ArcGIS JavaScript API v3.2 which doesn't request a tile unless the user is within the visible range). After a period of inactivity (I'm not sure how much) the first hit to the tile service takes a long time to respond. I can't say for sure if this only happens when generating a new tile on demand, when loading an existing tile, or both. The delay itself is annoying but it's also accompanied by an error in the AGS logs:               [Severe] Services containing process crashed. MapServiceName.MapServer instance has crashed      The logs provide no more information, and performance picks up after the initial crash. The map service's data source is a pair of File GeoDatabases, both are registered with the server, both are in EPSG:4326 and output in EPSG:3857 in line with the ESRI base map. Features are very simple (all polygons, usually four-sided shapes) and don't seem to have any geometry problems. The cache disk is nowhere near full and I haven't previously encountered any permissions issues with any of the AGS accounts. Can anyone suggest why the service might be crashing, why performance might be slow (if not connected to the crash), and where I might want to take my investigation? I have intentionally avoided the AGS forums as I would actually like someone to see my question before the new year ;) Thanks.


Osm2pgsql crashed on import. Can I continue with what has been imported so far?

I recently imported the osm planet file using osm2pgsql. The process crashed (hung my whole box) when creating the geometry indices and other indices on the tables. Scouring through the source code, it looks fairly trivial to add the indices by hand, however I was wondering if anybody knew if the database is consistent at this point (ie, all row data comitted). Don't want to have to start over with the whole import really.


QGIS crashes when adding vector layer from mysql or testing connection to mysql

For a while I have been adding layers from my mySql database and all went well. Then I got crashes when adding new vector layers. Then I only got crashes and it is not possible to add layers from the mySql database any more (crash) or create a connection, then test it (crash). Sporadically I can get to the list of the tables in the mySql database after creating a new connection but usually QGIS crashes. I tried deinstalling and reinstalling QGIS. Did not help. I am using the menu "layer" > "add vector layer" > "database" (in German, so maybe a little different). The table I am trying to add is named "tblTeilPopMassnPolygon", contains an id "TPopMassnPolygonId" with Primary Key and Auto Increment. Plus a polygon field "TPopMassnPolygonPolygon" and a field named "TPopMassnGuid", varchar 40, to save GUID's in. There are 6000 datasets in the table. I had tested them with the polygon data before and it worked. So there is probably nothing wrong with the data itself. I tried using shorter names for table and fields but that didn't help. How can I get QGIS to connect to my mySql tables again?


Create TIN from Elevation Points in Python

I want create a TIN from point features with elevation using a python script. This is my script:               import arcpy     from arcpy import env     import json as simplejson     import urllib     import sys,os          arcpy.CheckOutExtension("3D")          env.workspace = r"D:\test\Pythontest"          OutputTin = "PointFTin"     Delaunay = "constrained_delaunay"          arcpy.CreateTin_3d(OutputTin,{},"D:\test\Pythontest\PlotGrid_label_prj.shp"+ "Shape.Ele masspoints", Delaunay)      When I run this script I get this error:               RuntimeError: Object: Error in executing tool      Where is it going wrong?


MapInfo TAB - Integer fields and null values

Can MapInfo TAB integer fields hold null values? Is there a spec document on MapInfo TAB attributes available for download?


Conversion from N,E to Lat/Long And Lat/Long to N,E will always be same?

I have an excel tool for the conversion of following lat/long: Lat: 23° 11’ 24.67503’’ S Long: 118° 46’ 6.15742’’ E , which gives me following N,E: Northing: 7,434,317.43 Easting: 680,991.69 But I am not clear that why conversion of the N,E back to Lat/Long is not giving me same results i.e. Lat: 23° 11’ 24.67503’’ S Long: 118° 46’ 6.15742’’ E Can anyone help me please?


How to Modify Polygons in GeoDjango?

Just getting acquainted with GeoDjango. Have Basic Tutorial app created with sample world data. In trying out the admin, I appreciate seeing the polygon on the map and a couple basic tools for adding a new polygon and a tool for moving existing points. Question: What options are there for actually DELETING existing points? Are there other tools for easily creating/editing map polygons. I am targeting Google Maps API. My GeoDjango App will be the back end for generating Google Maps. Looking for options for creating/modifying boundaries/polygons to be maintained by GeoDjango. Thank you :)


What is the least painful way to style ENC data for the Web?

Currently working on a project that among other things requires displaying NOAA ENC for a web based intranet portal. We downloaded S-57 ENC data from NOAA site and successfully imported them into Postgres DB using "ogr2ogr". We now need to serve them up with the right colors and symbols that must follow S-52 standard. I have looked at OpenCPN and ArcGIS-ArcMap and other desktop tools which have beautiful renditions of S-57 data with colors, daylight/ night and symbols. I would like to do the same for a web app that is served by GeoServer. The styling is hopefully accomplished using Styled Layer Descriptor (SLD) files. I've looked at Atlas Styler, OpenGeo Styler etc. Nice tools but i'd hate to hand-craft them as there are a ton of layers and features, not to mention validating all of them to get it right. It would be nice if there is a way to download or import SLD files from good sources that correspond to S-57 layers. What is the least painful way to render ENC data on Web? Thanks


Modify feature: moving or change the shape

When the user make some edit on features on ArcMap, the user can make some changes like moving the feature ,edit the attribute or change the shape. how could i know if the user move the feature or changed the shape. I want to put different conditions on the editing according to what the user did.


Stats problem, QGIS heatmap

I am new to QGIS 1.8 and have been trying to make a heatmap. From what I can tell, the reason I am not getting the results I want is that the max value of the raster is erroneously measured at 2. Max should be 15-20 from what I can tell. I have tried to run the operation multiple times with different parameters and the max value is 2 every time. I have used the Identify tool to verify that the problem is not the cells themselves. Any cell with value above 2 is symbolized the same. Any ideas on how to fix this?


GeoServer 2.4.1: How to use JSONP in a WFS request?

I've ready that I can enable JSONP by setting ENABLE_JSONP in web.xml:                 <context-param>         <param-name>ENABLE_JSONP</param-name>         <param-value>true</param-value>       </context-param>      And I've read you can set outputFormat=json&format_options=callback:myCallback but I only get json back. Is this functionality broken/reverted? I'm on Windows 8.1.


How to determine aspect areas for different elevation bands in a basin in QGIS?

I am working in QGIS and have the dem rater, the vector basin boundary, the hypsometric curve and the aspect raster. I need to determine the north, south and east-west facing areas in each of the 15 elevations bands already determined. How can I do this in QGIS?


Query area without "AREA" field

In ArcGis 10.2 I'd like to show only polygons over 1m2, using something like a query. Can I do this without adding an AREA field and calculating/populating it in my table ?


How to activate a layer in QGIS layers legend

How can I activate a layer in qgis map legend automatically using python code . I am using               qgis.utils.iface.setActiveLayer(self.theLayerToEdit)      but the layer is still not checked in the maplegend.


What are problems & biases stemming from mathematics behind using Monte Carlo Markov Chain to make heatmap?

My buddy and I are using Monte Carlo Markov Chains to make a heat map - in the interest of rigour, and out of curiosity: "What are the problems & biases that stem from the mathematics behind using Monte Carlo Markov Chain to make a heatmap, and what alternatives or corrective measures would you use to remedy this?" I'm finding it difficult to find a comparison of different predictive algorithms ! What I want to know is what biases are induced by the mathematics behind this algorithm, and which simpler/more complex proceedures would one use to correct for these.


Is there an alternative to the Integrate tool that will not change one of the inputs

I have two polygon layers, and would like to snap the outlines of one of them (colored poly's) to the other (grey) without changing the grey in anyway (so Integrate seems out). It should only snap if the edge of the colored poly's is within a distance of the grey poly edges already, something like 200ft. There are several thousand of them that need snapped so it can't be a by hand method. Am I forgetting some tool that could do that since Integrate changes both inputs? I have ArcInfo 10.1 most if not all extensions and Xtools Pro. Thanks. ![example](http://i.stack.imgur.com/ABlLT.jpg)


Is there Open Source equivalent to Collector for ArcGIS to use with QGIS/PostGIS?

From my reading of the **Esri web pages** the **Collector for ArcGIS** is tightly tied to that platform. I would be delighted to learn that I am wrong about this. I have a strong interest in generic mobile data collection apps for iOS and Android devices with GPSs. Particularly ones that can collect data while out of coverage, and then sync when they get back in to coverage. Ideally, they would also be able to cache portions of the data locally for use when off line. What I envision is a web based application that allows one to specify your data layout, constraints, etc. and which will then generate the the backend database and allow access via a restful API. What I need to do is enter data from a form on an mobile device, preferably both iOS and android, which will update tables in a backend database, potentially something with GIS extensions such as PostGIS. That would enable me to easily incorporate the data into QGIS. It must be able to work offline and many of the areas we work have not 3G coverage. Ideally you would be able to also cache selected data on the device for offline use. Map integrations is not important. I have found some good cloud based solutions but I really want something where I have control of the back end data collection as some of the data is potentially sensitive. e.g. reporting sighting of rare reptiles which are a smuggling threat. What alternatives are the for mobile apps to get data into QGIS.


Legend Problem with arcGIS javascript API

I am new to arcGIS javascript API.I tried to display a simple dynamic map service with legend but its not displaying legend. I am not able to find error in my code. Please help me in resolving the problem.               <!DOCTYPE html>                     <meta http-equiv="Content-Type" content="text/html; charset=utf-8">         <meta name="viewport" content="initial-scale=1, maximum-scale=1, user-scalable=no">         <title></title>         <link rel="stylesheet" href="http://{Myserver}/arcgis_js_api/library/3.7/3.7/js/esri/claro/claro.css">             <link rel="stylesheet" href="http://{Myserver}/arcgis_js_api/library/3.7/3.7/js/esri/css/esri.css">

Zoom In

Zoom Out

Full Extent

Prev Extent

Next Extent

Pan

Deactivate

                               enter code here


Can QGIS merge CSV files with different colums?

I am trying to merge a month's worth of daily CSV files without using ArcGIS. Each file has ~25,000 entries so this is a high data load project. Sample data can be found here. Note that we are using the Prerun product. Arc can handle it (using simple merge geoprocess tool), but we have to make this open source as it will be used to teach students who don't have access to ESRI software. We tried a simple c# code but ran into the problem that **some of the days have extra columns** , which snafu'd our code. Before we spend the time to make our code much more robust to handle the disparity in the products, I wanted to ask if there was a way to merge tables (CSV files) within QGIS? This would be the simplest methodology for the people we're trying to benefit.


How to map flood areas from a set of points containing surface and height of water?

I am working on the damages made by the typhoon Haiyan. I use QGIS 2.0 with extensions and toolbox, and MapInfo. I have vectors, rasters of the place. I want to map the flood area along the coastline with analyze of the height of water. The picture below shows a result approaching what i would like to obtain ![area of inundation and height of water](http://i.stack.imgur.com/gez46.jpg) Beside i have a set of 70 points, with long/lat, area inundation, height of water ![table of data](http://i.stack.imgur.com/biNwm.jpg) that are represented on this map ![map of data](http://i.stack.imgur.com/O2AmQ.jpg) I suppose that a Satellite image has been used to get the points. My questions are:   1. is it possible to represent the areas of inundation with the attributes in the table?   2. If no, what is the method they could have used ? (may be not easy to guess it!)   3. If yes, what is the method to proceed?  Anyone has an idea ? Thanks a lot for your help:)


Loading Tiger & ACS Data Into CartoDB

Does anyone have experience loading the pre-joined TIGER shape files and ACS 5 year data files by block group into CartoDB? I have downloaded the data from http://www.census.gov/geo/maps- data/data/tiger-data.html by the 2007-2011 block group data, but the zip file isn't recognized in CartoDB.


How to allow more polyline vertices for a ST_GEOMETRY in Oracle/SDE?

I'm receiving an error writing geometry (in ArcObjects) that has a large number of vertices (around 900). So big, but not ridiculously big. I would expect this to be no problem and I don't recall this issue in SQL Server or PostGreSQL. The geometry storage is ST_GEOMETRY. The error is:   `ORA-24816: Expanded non LONG bind data supplied after actual LONG or LOB column` ..and it appears to be related to the number of vertices that the polyline has on it. Lots of vertices: I get the above-mentioned error. Less vertices: no problems. As a test, I generalized the geometry down to two points: success. I've started to work my way up and gotten up to 120 vertices with success. Is there a better way to determine that limit? Is there a way to increase the number of vertices that I can have on a line? Is there a DBTUNE parameter that manages this? I've googled the actual ORA error, but I'm uncertain as to how to apply this to an ArcSDE repository. A couple of notes:   1. I don't own or administer the ArcSDE repository and can't modify it.   2. Breaking the long polyline feature into several smaller ones is not an option ( _for this project_ ).   3. Write using `feature.Shape = geometry; feature.Store();`   4. Code executes as expected for most of the features--I have 2 that have high number of points that cause the ORA error. **Oracle version:** ` Oracle Database 11g Enterprise Edition Release 11.2.0.2.0 - 64bit Production   PL/SQL Release 11.2.0.2.0 - Production   CORE 11.2.0.2.0 Production   TNS for Linux: Version 11.2.0.2.0 - Production   NLSRTL Version 11.2.0.2.0 - Production`


How do I have logic in ESRI's popup configs for the Flex api/viewer?

The popup window config functionality I'm referencing is this... http://help.arcgis.com/en/webapps/flexviewer/help/index.html#//01m30000002q000000 This is my config...               <?xml version="1.0" ?>     <configuration>         <title>Wind</title>         <description>             <![CDATA[                 Forecast is valid for {DATA_DATE} UTC                 Wind Direction: {WIND_DIRECTION} degrees                 Wind Speed: {WIND_SPEED} knots                 Wind Chill: {WIND_CHILL} F             ]]>         </description>         <fields>             <field name="DATA_DATE">                 <format useutc="true" />             </field>              <field name="WIND_DIRECTION"/>                      <field name="WIND_SPEED"/>             <field name="WIND_CHILL"/>         </fields>     </configuration>      The problem is that sometimes I want to display N/A for the Wind Chill value. This only occurs when wind chill is -999. Is it possible to display something different than what is in the layer data in this config driven popup? If so, how? I'd rather not have the column for wind chill be a varchar and have the string "N/A" in it.


FME Desktop Pro 2014 : CityGML to KML Collada

Hi I'm pretty new at using FME (just started this morning to tell the truth) In a few words I'd like to export a CityGML model to KML/Collada files in order to load these in GoogleEarth and other 3D clients. In details : **Input** : A CityGML containing 300 buildings referencing textures AND DTM data I want to get rid of (ReliefFeature and TINRelief). The general structure of my CityGML is : `<CityModel> <Building> <BuildingPart> <Appearance> <SurfaceDataMember> ... </Appearance> <RoofSurface> ... <WallSurface> ... </BuildingPart> </Building> ... </CityModel>` **Output** : Idealy I'd like to generate a collection of KMZ files containing KML + Collada + Textures. Each KMZ file would represent a portion (a tile) of the original CityGML file. **My question** : Where to start ? :) **What I've already done** :   1. I've been able to load my CityGML file in FME. A certain number of FeatureTypes were imported (Building, BuildingPart, WallSurface, RoofSurface,CityModel).   2. I've been able to define a Collada writer and translate my data. As a result I obtained a .dae file for each featureType ... which is not really what I was trying to achieve => I thought of getting a Collada file containing each building with complete geometry (roofs+walls) as a first step. **So** :   * Are there any example FME project or tutorials that might help me to achieve this kind of advanced transformation ?    * At least, is there a way to iterate through the building list, get roofs+ walls and concatenate the geometry to produce a Collada containing the list of all building geometries ? Hope it's clear enough ;) Thanx in advance Emmanuel


How to get the nodes on the absolute boundary of a TIN with ArcObjects

I am trying to get only the nodes on the absolute boundary of a Tin. But GetNeighbor or LeftTriangle methods return edges whose nodes have NaN Z and very weird X-Y coordinates, so the row is inserted in the feature class without a feature. Is there any other method or am I doing anything wrong? Thanks, Mar Dim i As Integer For i = 1 To pTinAdvC1.TriangleCount               Set pTinTriangleC1 = pTinAdvC1.GetTriangle(i)              Dim j As Integer         For j = 0 To 2                  If pTinTriangleC1.Edge(j).GetNeighbor Is Nothing Then                  Dim pTinEdge As ITinEdge             Set pTinEdge = pTinTriangleC1.Edge(j)                  Dim pBorderP As IPoint             Set pBorderP = New Point             Set pZAwarePoint = pBorderP             pZAwarePoint.ZAware = True                  pTinEdge.FromNode.QueryAsPoint pBorderP                  Set pPointLayer = pScene.Layer(2)             Set pPointFC = pPointLayer.FeatureClass             Set pFeat = pPointFC.CreateFeature             On Error Resume Next             Set pFeat.Shape = pBorderP             pFeat.Store                  End If              Next          Next


QGIS Georeferencer and Coordinate Systems

I'm having an issue with Ubuntu and Windows versions of QGIS where the georeferencer does not properly assign the coordinate system to the output raster. When I assign the output CRS to 'EPSG:900913 - Google Mercator', the output is set to "EPSG:3395 - WGS 84 / World Mercator" Its not until i 'set layer CRS' of the newly georeferenced raster to 'EPSG:900913 - Google Mercator' that the raster draws correctly. What am I doing wrong? Thanks!


How to convert dbf to raster the short way?

I have got a dbf file and i want make it to raster file but shor way, Look at my long way, 1- First i add my dbf file 2- and i create tin http://i1211.photobucket.com/albums/cc421/mtnkou/2.jpg 3- at the and my raster file (tin to raster) http://i1211.photobucket.com/albums/cc421/mtnkou/3.jpg How can i do dbf to raster with short way, i dont want create tin file, Just dbf to raster, How can i do ? i am student and i apologize to everyone for this stuation. Thank you very much every answer.


Grayscale Bing Basemap - Is it possible in ArcMap?

Is it possible to achieve a grayscale look for the Bing basemap in ArcMap 10? I noticed that its possible with the maps hosted by ESRI since we can access the server that is hosting them. I haven't been able to find anything on Bing. Even using the OpenStreet Maps would be a great alternative, but again there doesn't seem to be an open server to use the trick posted on ESRI's mapping blog. **Edit** It looks like this is currently not possible, so I'm accepting the "its not possible answer". If this does become possible in the future please submit an answer!


Is there a simple way to detect what generation a MrSID file is?

I just found out there are a few types of MrSID files \- generations MG2,MG3,MG4. As I have a user which uses Bentley Microstation V8 and apparently can see only some of the generations (and thus I need to convert between the different generations)- I wanted to know if there is a quick way to detect what generation a SID file is (I have AutoCAD 2007, ArcGIS 9.3.1)? Thanks in advance.


Georeferencing an image with GDAL, based on QGIS Georeferencer Plugin

I would like to georeference an image with GDAL. The image already has a spatial reference set (EPSG:31276). For the georeferencing, I have a list of 4 ground control points (also in EPSG:31276). Using the QGIS Georeferencer Plugin, I can correctly georeference the image. I then select "Generate GDAL Script", which results in the following GDAL commands:               gdal_translate -of GTiff -gcp 6.55253e+06 -5.07402e+06 6.55273e+06 5.07395e+06 -gcp 6.55246e+06 -5.0711e+06 6.55263e+06 5.07106e+06 -gcp 6.55477e+06 -5.07382e+06 6.55511e+06 5.07368e+06 -gcp 6.55483e+06 -5.07197e+06 6.55525e+06 5.07197e+06 "C:/gis_sample_data/sample_image/image1.tif" "C:/Users/XYZ/AppData/Local/Temp/image1.tif"          gdalwarp -r near -order 1 -co COMPRESS=NONE  "C:/Users/XYZ/AppData/Local/Temp/image1.tif" "C:/gis_sample_data/sample_image/image1_georef.tif"      Those commands run without an error, but the resulting image is not correctly projected. Here is a small part from `gdalinfo "C:/gis_sample_data/sample_image/image1.tif"` (the original image):               Corner Coordinates:     Upper Left  ( 6551897.000, 5074277.000) ( 18d40' 4.21"E, 45d48'39.73"N)      whereas `gdalinfo "C:/gis_sample_data/sample_image/image1_georef.tif"` (result from the GDAL call) shows               Corner Coordinates:     Upper Left  ( -534350.794,  156938.268) ( 34d51'24.74"W,  0d52'47.63"N)      How would I have to modify the GDAL commands to correctly georeference the image? The sample image is part of the Baranja Hill data set, which can be downloaded from here. The ground control points are those:               mapX,mapY,pixelX,pixelY,enable     6552733.64613180514425039,5073946.04477077350020409,6552533.84482758585363626,5074021.48275862075388432,1     6552633.32449856773018837,5071059.01110315229743719,6552461.44827586226165295,5071104.32758620660752058,1     6555111.64040114637464285,5073678.52041547279804945,6554765.36206896603107452,5073821.32758620660752058,1     6555249.11819484271109104,5071965.62141833826899529,6554833.5,5071973.08620689623057842,1      When copied into a text editor and saved with the .points extension, they can be imported in the QGIS Georeference Plugin. EDIT: Changed negative to positive values in `gdal_translate`. Results of `gdalinfo`:               Corner Coordinates:     Upper Left  ( -591887.923,  223739.838) ( 35d 8'28.93"W,  1d14'58.74"N)      By the way, I also tried adding `-t_srs EPSG:31276`, `-s_srs EPSG:31276` and `-a_srs EPSG:31276` and various combinations thereof to the gdal commands, but it didn't fix the problem.


Numbering rows and columns in polygon grid (fishnet)

I seem to be a bit stuck on what i thought would be quite straightforward. I've created a fishnet of polygons, all 25 x 25 metres, 100 polygons by 100 polygons. My main task is to create a relative co-ordinate system based on the distance from a vertex somewhere near the middle of the grid (but NOT the middle), which will be the origin 0,0 Basically, the polygon square - where the origin sits at the very bottom left of it - will be X=0 and Y=0. The next square to its right will be X=25 Y=0 and so on and so on to the extents in all directions. The increasing blocks of 25 just represent the actual distance in metres from the origin. I thought I could maybe number the columns and rows of the polygon grid then work from there, assigning an X and Y value based on column and row number but I cant even seem to do that. Any ideas? thanks a lot


Can GeoRSS feed be used for raster image?

GeoRSS can be used to show vector data (point, line, polygon, etc). Can it be used to show raster image, say GeoTIFF? Is it part of the specification? Any idea? Thanks. Cheers, Alex


Edit imported osm maps on GeoServer with iD or Potlatch 2

I have installed geoserver on Ubuntu and imported osm map data into it. Now I want to attach a map editor(iD or potlatch 2) to this suite to edit maps but I have not a clue. Any help or documentation? In the near future I will need to get change sets and also upload my changes to osm. Is it possible to do this?


How to use GeoJSON from Geoserver?

Please help me on how can get my map from Geoserver. I want create my country web but I can nott see it in my browser but I can see in Geoserver. Please help me I am very new openlayers and Geoserver.                           <title>OpenLayers Example</title>


Android offline app to visualize GeoTIFF maps (or any well-known raster format)

I like to go to mountain, and got lots of GeoTIFF maps, which I can convert to other formats with GlobalMapper (JPG+JGW, ECW, ...). However, I haven't found an Android app to: \- show these raster images \- show my position, at the same time I tested QGIS Mobile, but it did not show my position. Any other alternative? Thank you very much and best regards.


How can I find a point inside a polygon in PostGIS?

**How can I find a point that is guaranteed to be within a given polygon in PostGIS?** I am aware of the `ST_Centroid` function. However, the centroid is not always within a polygon, see below: ![a centroid lying outside a polygon](http://postgis.refractions.net/docs/images/st_centroid03.png) Furthermore, I would like to avoid using a point that is on the polygon boundary, but rather want one that is inside the boundary (and not inside a hole in donut shaped polygons).


ST_HasBBOX does not exists

I have the lastest PostGIS version. When I query: `select ST_HasBBOX(the_geom) from table_geom limit 1` I get this error: `ERROR: function st_hasbbox(geometry) does not exist` Why this function does not exists


PostGIS & Geoserver: One large vector layer or many small layers?

I have a geoserver that displays some basic country data for Brazil at a course spatial scale (1:1 milion). Now the geography institute (IBGE) is making available the topo sheets at 1:250,000 scale, one shape file for each of the 548 tiles. My question is: In order to store those data in PostGIS and serve via Geoserver, would I be better of merging all topo sheets into one? Or is it better to have each sheet as a separate table in PostGIS? This, of course, considering that there will be no topology issues with the sheets. Thanks Daniel


Is there a FOSS solution for legend parsing in a web mapping environment?

Our web mapping requirements are fairly standard - some basic querying and the ability to print out the resulting map (or save the image). Typically we've used ESRI products in the past (right now we're using GeoCortex) but we still haven't been able to solve the problem of how to parse the legend information to match what is actually being shown on the map. For example, if there are four polygons showing up on the map, there should only be four legend items for the polygons. Is there an open source solution for this problem? **EDIT** : Let's see if I can clarify this more. Currently, we have a geology layer with 33,000+ polygons. The legend for this layer is a single image and contains several hundred legend items (ie. it's the legend for the entire geology layer, stored in a JPEG/PNG). Typically, a user will zoom in on an area (such as a single NTS map sheet, or even smaller) and so a subset of geology polygons will be shown on screen. They will then save the map as a print document (PDF or image), which is generated by GeoCortex/ArcGIS Server. The print document will contain a header, an image of the map and a legend. However, the legend for the geology layer will be the legend for the _entire_ geology layer, not what is actually shown on the map (ie. a much smaller subset). So, I'm looking to see if there is a **FOSS solution** to enable me to clip out or generate on the fly an subset of the geology legend which reflects the actual items on-screen/on the map, to go with the print document. I hope that clarifies things; I apologize if it doesn't!


In ArcMap 10.0, how do you vertical dock the Tools toolbar to the left of the map and right of the TOC?

In general, I really like the docking environment in ArcGIS 10. As I was showing off ArcGIS 10 to some of my end-users, I couldn't get the toolbar docked to the left of the map and to the right of the table of contents (TOC). Am I missing something here? I want to configure like ArcMap 9.3.1: screenshot I think as a workaround, end users can dock the toolbar to the left of the TOC, and auto-hide the TOC(personally, I will do this). Or they can dock the TOC to the right of the map (and maybe also auto hide the TOC).


ogrinfo returns a negative coordinate for extent of mif data file

I'm trying to get the extent of a MIF file using ogrinfo. The MIF file is supposed to be in epsg:23030 projection. This is:               PROJCS["ED50 / UTM zone 30N",         GEOGCS["ED50",             DATUM["European_Datum_1950",                 SPHEROID["International 1924",6378388,297,                     AUTHORITY["EPSG","7022"]],                 AUTHORITY["EPSG","6230"]],             PRIMEM["Greenwich",0,                 AUTHORITY["EPSG","8901"]],             UNIT["degree",0.01745329251994328,                 AUTHORITY["EPSG","9122"]],             AUTHORITY["EPSG","4230"]],         UNIT["metre",1,             AUTHORITY["EPSG","9001"]],         PROJECTION["Transverse_Mercator"],         PARAMETER["latitude_of_origin",0],         PARAMETER["central_meridian",-3],         PARAMETER["scale_factor",0.9996],         PARAMETER["false_easting",500000],         PARAMETER["false_northing",0],         AUTHORITY["EPSG","23030"],         AXIS["Easting",EAST],         AXIS["Northing",NORTH]]      But with ogrinfo, I get this as SRS WKT               PROJCS["unnamed",         GEOGCS["unnamed",             DATUM["European_Datum_1950",                 SPHEROID["International 1924",6378388,297],                 TOWGS84[-87,-98,-121,-0,-0,-0,0]],             PRIMEM["Greenwich",0],             UNIT["degree",0.0174532925199433]],         PROJECTION["Transverse_Mercator"],         PARAMETER["latitude_of_origin",0],         PARAMETER["central_meridian",-3],         PARAMETER["scale_factor",0.9996],         PARAMETER["false_easting",500000],         PARAMETER["false_northing",0],         UNIT["Meter",1]]      And this as extent               Extent: (-28050.000000, 3896708.000000) - (1139100.000000, 4874250.000000)      Outside the limits of epsg:23030 projection, and with negative coordinates. I don't know if it's related, but the data file covers all Spain, and epsg:23030 doesn't cover all the country. There's a small part covered by epsg:23029. So, I guess the problem is related with the way the file was exported or georeferenced in the origin (MapInfo). So, what could be the problem with my data?


geoserver wms layer is very slow

I am using open-layers to display a map and I have vectors layer with some features that I'm loading with ajax from MySql, and I'm using geoserver to load all the background from shp files(like buildings streets etc. ) with default style, like this:                       var generalLayers = new OpenLayers.Layer.WMS(                     "general", "http://XXX/geoserver/city/wms",                     { "transparent": "true",                         LAYERS: 'city:Fence, city:Buildings, city:streets',                         format: format,                         tiled: true,                         tilesOrigin: map.maxExtent.left + ',' + map.maxExtent.bottom                     },                     {                         buffer: 0,                         isBaseLayer: false,                         displayOutsideMaxExtent: true,                         yx: { 'EPSG:2039': false },                         visibility: true                     }                 );      the features from MySql are loading after 1.7 seconds, but the WMS layers is completely loading after 19.6 seconds!! (12 seconds after the features), all this time i have only the features without the background. I have VPS so its now because the server. why it taking so much time? how can i improve it? its because im loading from shape files? if i will import everything(the shp files) to the db and I will load the background with MySql plugin for geoserver as WMS layer its will be faster? thanks.


After import some data from OSM file into Postgre, the coordinates numbers are too big and having trouble to cast to Geography:               Coordinate values are out of range [-180 -90, 180 90] for GEOGRAPHY type      In my table:               select ST_AsText( the_geom )  from routes limit 4          "LINESTRING(-4902130.86 -2265840.91,-4902083.04 -2265907.16)"     "LINESTRING(-4903296.22 -2266295.18,-4903363.69 -2266333.39)"     "LINESTRING(-5705594.67 -3518944.1,-5705523.22 -3518937.37)"     "LINESTRING(-5484056.84 -1884629.16,-5484093.47 -1884510.24)"      I think the numbers should be between (-90 / 90, -180 / 180) ... Info: 1: select ST_SRID( the_geom ) returns 900913 for all lines. 2: select UpdateGeometrySRID('public', 'routes', 'the_geom', 4326) does not converts the numbers, but SRID now shows 4326.


Weird Error when setting the source path in the "PICTURE_ELEMENT"

I am trying to set the `sourceImage` to a tiff file using this code:                       m=1             pic = arcpy.mapping.ListLayoutElements(mxd,"PICTURE_ELEMENT" , str(m))[0]             pic.sourceImage = r"D:\PDF\tables\{0}.tiff".format(str(m)) ##error here      I get this error: `IOError: Cannot set path based on current settings.`


Using Calculate Field (Data Management) tool to calculate 3D Length in ArcGIS 10.1

I'm relatively new to the world of model builder and scripting, so it's likely there is a very simple solution to the question I've been trying to figure out. I'm trying to create a tool that takes a given shapefile, interpolates it over a given surface, adds a "Length_3D" field, then calculates that length (3D) in the attribute table using the Calculate Field tool in the Data Management Toolbox. I've done the same for 2D length using !shape.length@miles! as the expression and I have no issues there. What I can't seem to figure out is how to write an expression that will allow me to calculate the 3D length from my already interpolated line. Right now the model runs all the way through until this point then throws an error. I'm able to open the newly created ZM shape and calculate the 3D length by using the "Calculate Geometry" tool in the attribute table then selecting "3D Length". I'd really like to be able to perform this in the model though. Has anyone else had a similar issue or know how to write the expression? I've tried "!shape.3DLength@miles!, !shape.3D_Length@miles!, and !shape.3D Length@miles!" with no luck. Thanks!


Bathymetry dataset of Italy

I am looking for a dataset of bathymetrics data of Italy. I shall use them into an application for collision avoidance of vessels into the port. I need to know how deep the sea bottom is in each point, and that's why data should be as much as possible precise and graided. So I consider only data more resoluted than 3 arc-second. Can someone help me to find open source or not too much expensive dataset with this constraints?


Georeference image is mirrored

![enter image description here](http://i.stack.imgur.com/NbgAd.jpg) I'm working with QGIS and I'm trying to geo-reference a tobin lease map. The map datum: NAD27, projection: UTM Zone 14N. The process to georeference is fairly straight forward but the results I'm getting is off. The program is basically mirroring the image. I have tried it three times with the same results. I'm not sure what I am doing wrong.


Load LAS files in MapInfo

Is it possible to open LAS files within MapInfo? If so, how is accomplished in MapInfo 10.0?


ArcGIS flow accumulation

I am trying to run ArcGIS 10.1 flow accumulation on a 10000 x 10000 cell unsigned 8Bit integer flow model (but 80% of that is no data). The data has been through a "fill" procedure once. The only unusual aspect is that I have inverted the DEM, so the water flows from 0 into negative values.. could be an issue?? The model has been running for 24 hours, but while ArcGIS claims it is processing, the GPU shows no activity and the output file is empty. My machine is a 2Ghz quadcore, 64 bit Win7 OS. How long might I expect this to run, and is there anyway to tell if it is working at all? Thanks


How to build pgRouting with CGAL?

When running cmake . for CGAL 3.9, i get the following output:               root@postgresql local/CGAL-3.9# cmake -DWITH_CGAL_Qt3=OFF -DWITH_CGAL_Qt4=OFF -DWITH_CGAL_ImageIO=OFF .     -- The C compiler identification is GNU     -- The CXX compiler identification is GNU     -- Check for working C compiler: /usr/bin/gcc     -- Check for working C compiler: /usr/bin/gcc -- works     -- Detecting C compiler ABI info     -- Detecting C compiler ABI info - done     -- Check for working CXX compiler: /usr/bin/c++     -- Check for working CXX compiler: /usr/bin/c++ -- works     -- Detecting CXX compiler ABI info     -- Detecting CXX compiler ABI info - done     -- Build CGAL from release: CGAL-3.9     -- Packagenames: CGAL-3.9     -- CGAL_MAJOR_VERSION=3     -- CGAL_MINOR_VERSION=9     -- CGAL_BUGFIX_VERSION=     -- CGAL_SONAME_VERSION=8     -- CGAL_SOVERSION     =8.0.0     -- CGAL_REFERENCE_CACHE_DIR=     -- Building shared libraries     -- Targetting Unix Makefiles     -- Using /usr/bin/c++ compiler.     -- USING CMake version: 2.8.0     -- System: Linux     -- USING GCC_VERSION = '4.4.3'     -- Using gcc version 4 or later. Adding -frounding-math     -- Build type: Release     -- USING CXXFLAGS = ' -frounding-math -O3 -DNDEBUG'     -- USING LDFLAGS = ' '     -- Could NOT find GMP     -- Could NOT find MPFR     -- CGAL_Core needs GMP, cannot be configured.     CMake Error at cmake/modules/FindBoost.cmake:917 (message):       Unable to find the requested Boost libraries.            Boost version: 1.40.0            Boost include path: /usr/include            The following Boost libraries could not be found:                    boost_thread            No Boost libraries were found.  You may need to set Boost_LIBRARYDIR to the       directory containing Boost libraries or BOOST_ROOT to the location of       Boost.     Call Stack (most recent call first):       cmake/modules/CGAL_SetupBoost.cmake:12 (find_package)       cmake/modules/CGAL_SetupDependencies.cmake:32 (include)       CMakeLists.txt:422 (include)               -- Boost include:     /usr/include     -- Boost libraries:     -- Boost definitions:     -- USING BOOST_VERSION = '1.40.0'     ...      ive tried various env variables. I am wondering if i need to install / compile some other boost thing? I have some libboost* in /usr/lib Any suggestions on things to try would be appreciated.


Georeference to x y coordinates with one lat long origin?

We are running parallel transect lines for a survey, which will roughly cover a rectangular area. When making the transects, we will periodically capture spatial locations of features and events we encounter. Due to technical issues, we can only record the latitude/longitude coordinates of the very beginning of the first transect (i.e. one corner of the rectangle). The locations of all our logging of other features/events will be recorded relative to the origin, with x y coordinates (running along the east- west north-south axes, respectively). I would like to georeference all our data points from the survey into a map, where I can click on any data point and get its coordinates in latitude and longitude. How should I do this? I prefer a solution that could be done in QGIS 1.7+, but ArcGIS 10 is OK, too. Thanks!


Spatialite android trigger problem

I am developing an android application that uses the database spatialite, I use this project to integrate the database on android https://code.google.com/p/spatialite-android/ It works normally, except when I use the clauses' UPDATE 'or' INSERT' it there's the trigger that prevents the request to be executed, i get the following message:'geometry violates Geometry constraint [geom-type or SRID not allowed]. whatsoever geometries that I want to update or just attribute data, but when I execute the same clause in spatialite_gui for example, it works normaly. I tried to delete the trigger, the problem is solved, but after it is impossible to execute spatial queries, the table is corrupted. I can't understand or locate the problem,to try to solve it. thank you for helping me. Edit: my question is: why the trigger prevents 'insert' and 'update' on android and not on saptialite_gui?


Any plans to change GP toolbox format (.tbx file) to text-based format?

Have you heard if esri has any plans to change GP toolbox format (.tbx file) to text-based format in future? The issue is that we can't track changes in binary tbx file using source control system. With text based format it would be easy. Most of our tools are .NET (COM) based. Some used for printing are Python based. Therefore we use a custom toolbox (.tbx), not Python toolbox (.pyt). Or alternatively (maybe this would be even better approach) esri can provide a tool buildTbx which takes some text-based toolbox descriptor file, finds what tools should be in toolbox, what are inputs and outputs, also other files needed to be included in toolbox (.py files) and creates that binary tbx file as an output.


How to Create Polygon shapefile based on change in point value?

Is there a function in ArcGIS that would take analyze the adjacent point or polygon based on percent difference. If the adjacent value is within less than the percent it would dissolve otherwise create a new area. Maybe I need to use raster data or stay with point data? Or complete data analysis outside ArcGIS before creating a polygon shapefile based on change in value.


Using ArcMap 10.1, is it possible to create a mapbook with face to face maps showing before and after digitizations?

I want to create a mapbook that displays changes in river flows over time. I have one shapefile of my recent digitizations, as well as ones done prior. I was hoping to display them, new on the left and old on the right. Does anyone know if this is possible?


How create geographic column in postgis table?

I merged several point shape files in postgis, but when I load the new merged shape in QGIS it is not possible to see points in the main window. To create the merged table I used this piece of code found in the forum: CREATE TABLE merged (id serial primary key, attrib1 integer, attrib2 varchar(15),....); SELECT AddGeometryColumn('merged','geom',,','XY'); INSERT INTO merged (attrib1, attrib2, ...., geom) SELECT attribA, attribB,...,geom FROM table_1; INSERT INTO merged (attrib1, attrib2, ...., geom) SELECT attribA, attribB,...,geom FROM table_2; I think that the problem is that in the merged table the geographic column is empty. I tried to populate this column, but my attempt filed. How can I populate this column? Any help will be appreciated. Thanks in advance


Background map as open data from geoserver

I'm in the process of opening up as much data to the public as possible. I've been told to open up for the background map (tiles?) which is used from a geoserver(.org), so others can use the maps/tiles? for their own projects. Which way would be the best way to do it? Is it a collection of images which should be exported, is it some sort of vector-graphics or documents which can be exported or...? I really do hope such a question is OK, I'm in no way a GIS expert (or GIS- user for that matter).


How to get a list of map scales programmatically in ArcMap?

I need to get all map scales, which are shown in the ArcMap scalebar. I need list of value, something like:               1:1500     1:3000     1:6000     1:12000     1:24000     1:50000     1:100000      Does anyone know how to get it?


nearest neighbor within table using PostGIS

using PostGIS 2.0 and PostgreSQL 9.1 I have a table of photographic observations, sco_photos. The location column contains a PostGIS geometry point specifying where the observation was made. The table is not enormous (about 7,000 rows). I am trying to write a query that will return, for each observation, the nearest other observation in the same table, along with the distance between the two observations. I can easily do this for a fixed point, e.g., how far is each observation from some central point, but the trick is locating the nearest other point in the table and measuring the distance to it. This is the closest I've come thus far:               SELECT a.photo_id, b.photo_id, st_distance(a.location, b.location)      FROM sco_photos a, sco_photos b      WHERE st_distance(a.location, b.location) =           (select min(st_distance(a.location, c.location))             from sco_photos c where c.photo_id <> a.photo_id);      ...but this seems just to hang. Performance for this application is not terribly important, as I need this for analysis not for real-time querying. Thanks for any advice.


QGIS Profile from Line tool requiring shapely.wkd?

I am trying to install the Profile from Line tool for QGIS, and the system informed me I need shapely.wkd. A search on the QGIS site for shapely takes me to the Contour Plugin...can anyone help me out with this?


OpenLayers: Show a styled title similar to pop-up but always visible?

I need to basically create boxes with titles inside, set to different styles (different coloured box, different FontAwesome icon etc.) set to specific points on the map, ideally replacing the label titles for LinePoint features. What would be the best/most ideal solution? I presume creating pop-ups that are always "popped up/visible" - but I don't want the little 'arrow' directly below (nor the white box outline) - I suppose those can be disabled by CSS ? Except I'm using the canvas rendererer (I can change to dom or webgl, I figured canvas looked nicer but unsure which works best for which situations)


How to export coordinates of selected nodes with ArcMap?

In ArcMap, I imported a .shp file with .dbf file. In the .dbf file, there is an attribute that denotes the node type. In ArcMap, by the menu "Selection > Select by Attributes", I chose the nodes that I want. The nodes chosen became light-blue. Now, how can I export the location coordinates of the chosen nodes with their node names? The node names are attributes recorded in the .dbf file. This is my first time to use any GIS software, so please be patient.


openlayers 3 function resolution

I'd like to make visible my var text between two resolutions so to speak between 150 and 50 for instance. but I didn't find how to do this yet.               var communes = new ol.layer.Vector({             title:'communes',             source: new ol.source.GeoJSON({                      url: 'communes.geojson'}),             style: function(feature, resolution) {     var text = resolution < 150 ? feature.get('nom_comm') : '';     if (!styleCache[text]) {       styleCache[text] = [new ol.style.Style({         fill: new ol.style.Fill({           color: 'rgba(255, 255, 255, 0.2)'         })      Any help ? Thx


Feature "Add to wishlist" in Info Window in Google Maps API 3?

I have Google Map and a lot of markers. Is there possibility to add feature when user clicks on marker and Info Window shows up, that user has option to add that location to wishlist and associate that to either a user ID/cookie?


Open NetCDF files in ArcGIS using Multidimension Tools

as in the title mentioned, iam looking for a way to get *.nc data displayed and managed in ArcGIS. The downloadsource of my data is NCEP/NCAR Reanalysis 1: Surface Flux To be a bit more accurate i want climate information, like wind speed and solar radiation for a given timespan. Let's say i want the wind speed information from the 01/01/2013 - 31/12/2013. My first step was to download the file representing year 2013 from Wind speed 2013 about 50 megabyte If im using the "Make NetCDF-Raster Layer" from the "Multidimension Tools" i get follwing error: > ERROR 000276: One or both dimensions have variable spacing in their > coordinate values Failed to execute (MakeNetCDFRasterLayer). When i set the dimension value to time, shown in figure 1 ![Make NetCDF with Dimension = time](http://i.stack.imgur.com/84HVA.jpg), i get a different error: > ERROR 999998: Unexpected Error. Failed to execute (MakeNetCDFRasterLayer). Anyway, i think this is due to the extent of the file. But iam not exactly sure about the error messages. Although iam able to create a Subset of the specific file Subset Link, iam not able to get the data over the year 2013, because he will average the time span iam entering into one file. What iam looking for are multiple tables, each of it representing the wind speed values at a specific time. The goal is to get the information about the windspeed for one year for 100 points located in europe. So if i have 100 points i have 100 tables, row numbers = timesteps and column would be the windspeed value. # Update 1: 2. July 2014 14:51               Dimensions:     uwnd                1457 19 8     time                1457     lat             19     lon             8      Variables: uwnd               parent_stat         Other     long_name           4xDaily u-wind at 10 m     valid_range         -32765 -12325     least_significant_digit     1     level_desc          10 m     GRIB_id             33     unpacked_valid_range        -102.199996948242 102.199996948242     dataset             NMC Reanalysis     missing_value           32766     statistic           Individual Obs     GRIB_name           U GRD     actual_range            -23.5 28.2999992370605     var_desc            u-wind     add_offset          225.449996948242     precision           2     units               m/s     scale_factor            0.00999999977648258      time               long_name           Time     axis                T     standard_name           time     delta_t             0000-00-00 06:00:00     actual_range            2013-01-01 00:00:00 2013-12-31 00:00:00     units               hours since 1-1-1 00:00:0.0      lat               long_name           Latitude     axis                Y     standard_name           latitude     actual_range            69.5216979980469 35.2374992370605     units               degrees_north      lon               long_name           Longitude     axis                X     standard_name           longitude     actual_range            345 358.125     units               degrees_east      NC0                4.2.6      Is there anyway i can get the data i want to? It would be sufficient to have the data in table format, but it would be nice to display them in ArcGIS as well. Any suggestions on this topic are highly appreciated.


Python error message after download (QGIS 1.8.0/ Windows 7)

Hello to all members of the forum, I downloaded from the OSM website a small area in the neighborhood of Strasbourg, France: X (longitude): 7.6808 7.8216 Y (latitude): 48.5075 48.52 I got the two regular .osm and .osm.db files, but I also got an error message: > "Une erreur est apparue lors de l'exécution du code python : > > Traceback (most recent call last): File "C:/PROGRA~1/Quantum GIS > Lisboa/apps/qgis/./python/plugins\osm\OsmPlugin.py", line 289, in > downloadOsmData > self.dbm.addDatabase(self.dbFileName,self.dlgLoad.pointLayer,self.dlgLoad.lineLayer,self.dlgLoad.polygonLayer) > File "C:/PROGRA~1/Quantum GIS > Lisboa/apps/qgis/./python/plugins\osm\OsmDatabaseManager.py", line 212, in > addDatabase self.dbConns[key]=sqlite3.connect(key) OperationalError: unable > to open database file > > Version de Python : 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 > bit (Intel)] > > Version de QGIS : 1.8.0-Lisboa Lisboa, f440b60 > > Chemin vers Python : ['C:/PROGRA~1/Quantum GIS Lisboa/apps/qgis/./python', > 'C:/Users/cfermanian/.qgis//python', > 'C:/Users/cfermanian/.qgis//python/plugins', 'C:/PROGRA~1/Quantum GIS > Lisboa/apps/qgis/./python/plugins', 'C:\PROGRA~1\Quantum GIS > Lisboa\bin\python27.zip', 'C:\PROGRA~1\Quantum GIS > Lisboa\apps\Python27\DLLs', 'C:\PROGRA~1\Quantum GIS > Lisboa\apps\Python27\lib', 'C:\PROGRA~1\Quantum GIS > Lisboa\apps\Python27\lib\plat-win', 'C:\PROGRA~1\Quantum GIS > Lisboa\apps\Python27\lib\lib-tk', 'C:\PROGRA~1\Quantum GIS > Lisboa\apps\qgis\bin', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27', > 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages', > 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\PIL', > 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32', > 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32\lib', > 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\Pythonwin', > 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\wx-2.8-msw- > unicode', 'C:\PROGRA~1\Quantum GIS > Lisboa\apps\qgis\python\plugins\fTools\tools']"   1. Is it a default error message from QGIs ?   2. the names of streets & roads do not appear in the QGis main viewer, is it normal ?   3. could it be a side-effect of the python error above ? Great if somebody has a clue...


Export WMS-layer to dwg or dxf

Is it possible to export a WMS-layer to dwg or dxf. I would like to use lines, circles etc. in Autocad, but i can't find out if it is possible?


What are data visualization and analyses pros/cons of QGIS over ArcGIS for Desktop?

Many of my colleagues use GIS on a very sporadic basis for simple data visualization and basic analysis (clipping, buffering, etc.) and I am hoping to transition these non-core GIS people to QGIS so that we aren't relying solely on ArcGIS and it's associated licensing costs. I've used QGIS a bit in the past, but was hoping that some of you could explain the pros/cons of using it over ArcGIS. What if someone wanted to do more advanced raster analysis? What is the support like for data formats? How about interfacing with databases? I'm just looking for a high-level overview and any hints/tips/pitfalls you may be able to provide information on.


Create FeatureClass in Personal Geodatabase and store data

I have to create a new featureClass in my existing personal GeoDatabase. And save some polylines into this featureClass. How can i do it? I need c# code to do this. Edit: I have a problem during Geometry type featureClass creation. Here the code i'm using:               IFields fields = objectClassDescription.RequiredFields;     ESRI.ArcGIS.Geodatabase.IFieldsEdit fieldsEdit = ESRI.ArcGIS.Geodatabase.IFieldsEdit)fields; // Explicit Cast     ESRI.ArcGIS.Geodatabase.IField field = new ESRI.ArcGIS.Geodatabase.Field();     ESRI.ArcGIS.Geodatabase.IFieldEdit fieldEdit = (ESRI.ArcGIS.Geodatabase.IFieldEdit)field;     fieldEdit.Name_2 = "Shape";     fieldEdit.Type_2 = ESRI.ArcGIS.Geodatabase.esriFieldType.esriFieldTypeGeometry;     fieldEdit.IsNullable_2 = true;     fieldEdit.AliasName_2 = "SampleShap";     fieldEdit.DefaultValue_2 = "test";     fieldsEdit.AddField(field);     fields = (ESRI.ArcGIS.Geodatabase.IFields)fieldsEdit;     ESRI.ArcGIS.Geodatabase.IFieldChecker fieldChecker = new ESRI.ArcGIS.Geodatabase.FieldChecker();     ESRI.ArcGIS.Geodatabase.IEnumFieldError enumFieldError = null;     ESRI.ArcGIS.Geodatabase.IFields validatedFields = null;     fieldChecker.ValidateWorkspace = (ESRI.ArcGIS.Geodatabase.IWorkspace)workspace;     fieldChecker.Validate(fields, out enumFieldError, out validatedFields);     IFeatureClass featureClass = featureWorkspace.CreateFeatureClass("Test01", validatedFields, null, null, ESRI.ArcGIS.Geodatabase.esriFeatureType.esriFTSimple, strShapeField, "");      I get this error:               The module of database management can't find the table or the input query 'GDB_DataChanges'      How can i solve it?


QGIS users - What plugins do you consider to be essential?

The title says it all - vector, raster, geoprocessing plugins?


Is there an ArcGIS tool equivalent to "Dissolve" that can be used on tables?

I have a data table that I need to trim down, merging records which have the same value in an ID field and performing a couple math operations on various other fields (summing, averaging). The Dissolve tool in ArcGIS does this for feature classes, but it doesn't work on tables. I realize there's probably a tool in ArcGIS that does this (it's basically a Dissolve with no geometry) but I haven't found it and am in a bit of a hurry. Is there one? I have an ArcInfo license available if necessary.


Problem - ST_AsGeoJSON return text with accolade

I am saving a result of ST_AsGeoJSON(..) into csv file. But the results I got can't be interrupted as an object later on. How can I polish and delete the accolades that I don't need. "{ **""* _type** "" **:** "" **MultiPolygon** "" **,** "" **coordinates** ""_*:[[[[-74.0758273156424,40.5922324042084],[-74.0761549604646,40.5921150798826],[-74.076200629116,40.5921869978165],[-74.0760676410383,40.5922346191496],[-74.0758729840191,40.5923043204707],[-74.0758273156424,40.5922324042084]]]]}"


Topo To Raster executing but not creating the raster

I am using **Topo to Raster by File** (ArcGIS version 10.1) via arcpy on a large stream and point dataset. Arcpy runs fine and executes without issuing any error message, however it does not create the raster as it normally does. When I look in the project workspace I can see that the S File that Topo to Raster creates temporarely (before deleting it when execution is terminated) have been created and its size is slightly over 2GB. Not sure what is the S File used for as it is a binary and can't inspect it. Is there a limit on the size of the S File created by Topo to Raster? What can I do to get around this issue in addition to reduce the size of the point dataset? Note that: There is enough space in the hard disk thereore the issue is not space; and The tool runs ok when I reduce the size of the dataset (option which I would like to avoid).


Does India Have a CORS (Continuously Operating Reference Stations) Network?

Does India have a CORS (Continuously Operating Reference Stations) network? If yes, can anybody provide more details about it?


OSX: Installing GRASS addons in qGIS?

I would like to install r.stream.basins in GRASS through qGIS (2.0.1-Dufour). I tried using `g.extension ext=r.stream.basins` and got the error "the 'make' program is required. Please install first". Is there are better way to install addons from within qGIS? I'm really having trouble figuring this out from the available documentation but it seems like it should be fairly simple. Any help would be greatly appreciated, thanks.


How to display an elevation profile and store the DEM elevation back to the KML line?

When I draw a track in Google Earth, for example along a road, I can right- click it to "Show elevation profile", and it displays a 2D plot of distance vs elevation. However, if I save the KML, the elevation coordinate is always zero. So let's suppose I import the KML LineString in QGIS, and I have a SRTM DEM layer, how can I:   1. Display an elevation profile (distance along road vs DEM elevation at point);   2. "Replace" the (zeroed) elevation profiles of KML by DEM elevation values, so that I could export back to kml and have elevation "hardcoded" in the LineString.


how to stack landsat 8 layers correctly?

i downloaded a landsat8 satellite image from USGS, it was separated into 12 layers. So, what i would like to do is to stack some of them. To do that, i am not sure which layer should be excluded from the stacking. Any help?


How to determine the accuracy of getFeatureInfo on a WMS layer?

I have added a WMS layer to the ArcGIS Server JSAPI and am attempting to query it using `getFeatureInfo` when the user clicks on the map. A reproducible sample is available here. This is mostly working, except that the accuracy of the results is very suspicious. If I click within a few hundred metres of a polygon boundary, I often get back the neighbouring polygon's results: ![enter image description here](http://i.stack.imgur.com/BoGCy.png) In this case I have clicked within Staffordshire Moorlands 010B (as shown by the infoWindow's anchor) but the result shows the neighbouring polygon 011C. In order to run the `getFeatureInfo`, I am calculating the map extent in its native web mercator, converting this into the units of the WMS layer, and using this as the bounding box. Is this a valid approach? The question ArcGIS Server WMS GetFeatureInfo tolerance implies that there may be some trickery required to get the tolerance right? Edit: here is a new sample which doesn't use any reprojection, so both requests use the same bbox, SRS, etc. Here is the GetMap request and the GetFeatureInfo request (the location of the click is about the same as the screenshot, ie within the 010B polygon, though the GetFeatureInfo result says 011C)


Openlayers panel control editing after initialization

How do I reference controls that have already been added to a panel. I am attempting to destroy or modify the default controls associated with the editingToolbar. These changes would be done outside of the initialization function. EDIT:My intent was to change the layer that a drawFeature controls was adding features too.


How to identify misaligned corners in PLSS sections

Sometimes PLSS section corners do not align ![enter image description here](http://i.stack.imgur.com/rIvIR.jpg) How could I automate the detection of these corners and generate a new shapefile showing where this occurs? I have a shapefile of 147k section polygons that I would like to analyze. I am on ArcGIS 10.1 and have Spatial Analyst and 3D analyst licenses. Based on comments below, I know I could count the sides of polygons to detect some of these occurrences. This would not work in every case where corners are misaligned. Any other suggestions? On a side note, why do the corners in PLSS sections not always align?


Is there any working example for pgr_trsp ?

I'm trying to use pgrouting pgr_trsp function but result is useless. Exactly same code with pgr_dijkstra working OK and find shortest path between start and stop but with pgr_trsp it is look like random jumps around map. I'm using this code `SELECT a.seq, a.id1 AS node, a.id2 AS to, a.cost, cr_2po_4pgr.osm_name, ST_AsText(ST_GeomFromEWKB(decode(cr_2po_4pgr.the_geom,'hex'))) as geodata FROM pgr_dijkstra ('SELECT gid as id, source, target, cost FROM cr_2po_4pgr', $start,$end,false, false) as a left join cr_2po_4pgr on a.id2 = cr_2po_4pgr.gid` Based on documentation trsp should be almost same as dijkstra but i can use turn restriction code. But even basic is not working. What i'm doing wrong ?


Is there an open source alternative to ArcGIS for Desktop's Integrate tool?

I need to snap about 1000 lines in one layer to the lines in another layer (there are about 5000 lines in the second layer). The lines to be snapped are about 15 map units away from the base layer lines. It looks like ArcGIS could do this using the **Integrate** tool in the Data Manager toolbox. Is there a QGIS/GRASS/SAGA/other open source alternative to Integrate or will I need to crack open my Python skills? Thanks for any help/advice on this


How can I merge postcode sector polygons into larger territories in QGIS?

I'm new to GIS so patience might be required :-) I'm building a basic web-based mapping tool for business use, to be used as a demo to show that we need one but managers don't know it yet! So far I've built a HTML page displaying an OpenLayers OSM base layer. I've been supplied with some postcode sector polygons in shapefile format that I'm wanting to merge into our business territories (wards, districts etc.). Aims were to plot these as layers on top of OSM that the user can switch between etc. Using QGIS, i don't seem to be able to import or link to a table of our geography hierarchy and dissolve the polygons by these fields. Is this best achieved by importing CSV or linking to db etc. I'm a bit lost :-) our geo hierarchy:               Disrict, Ward, Neighbourhood, PC Sector     11, 11A, 11A15, AA1 1 etc.      Any advice would be gratefully received :-)


What are some of the OSX GIS tools

Most Opensource tools are cross platform; however, some GIS tools such as ArcGIS are platform specific (windows only). What are some of the GIS tools that are usable on the mac? **EDIT** Although thank you for putting the basic apps in your answers, my question specifies that Most if not all of opensource tools are cross platform, I was really hoping to get more obscure libraries, packages, and programs that I didn't know about on OSX similar to dotspatial and sharpmap for windows


Elevation graph from GPX file

I've got many GPX files in the following format:               <?xml version="1.0" encoding="utf-8"?>     <gpx xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xsd="http://www.w3.org/2001/XMLSchema" version="1.1" creator="c174a4ce-d611-43ef-b4d9-e4122bc8f2a4" xmlns="http://www.topografix.com/GPX/1/1">       <rte>         <name>Withens Clough, Cragg Vale and Stoodley Pike</name>         <rtept lat="53.705918666306154" lon="-2.02137437495311">           <ele>268</ele>         </rtept>         <rtept lat="53.705457119907116" lon="-2.0253609418514311">           <ele>293.8</ele>         </rtept>         ....       </rte>     </gpx>      I use this is a website to generate a route using OS Openspace maps. However, this doesn't create an elevation profile. **Is there a web service or something I can run on my web server that can generate GPX elevation graphs for embedding into webpages?** My web server runs PHP. I've found plenty of Google map ones but don't need the map, just the graph, but in the implementations I've looked at it seems like one depends on the other and don't want to compromise loading times by including two maps or remove the OS maps which are a lot better, in my view.


What is the best web editor and server for postgis?

I'm looking for a web editor (similar to openstreetmap) that works with an map Server (postGIS). This should permit edit and create new structures on map. What is the best option? Can you help me? Thank you for your time.


Geocortex using MS WF

I am going to be looking at the new Geocortex software soon, and wondered if anyone else was using it, what they think of it and if they had any online resources they could share. Thanks in advance


how to convert UTM northing and longitude to latitude

I am doing some work with the Military Grid Reference System (MGRS), which is based on Universal Transverse Mercator (UTM). I find that I need to convert from UTM northing and longitude (actually, delta angle between desired longitude and appropriate UTM central meridian) to latitude. I've searched for days and cannot find this conversion. I need an algorithm. Can anyone point me to a reference to such an algorithm, or to open source software that does this conversion?


How to use CQL_FILTER for the highest value in a layer field?

I need to query/filter my wms in openlayers to only see the polygons classified with the highest values in a given field. Would be something like:               wms_layer.mergeNewParams({'CQL_FILTER': "ND = max(ND)"});      But this is not working. Any ideas?


How to convert OSM polyline data to polygons?

I just now downloaded OSM data .I am having the data as polylines. I want that polylines as polygons is it possible in Qgis?below I also posted the screen shot of the polyline also.![enter image description here](http://i.stack.imgur.com/MPkJK.jpg)


Splitting polygons in OpenLayers is not working

I have some WFS layers in OpenLayers, containing Polygons. I try to make a tool to split these polygnos and save them to db.               //The layer which I want spliting work on     var contract_layer = map.getLayersByName("contracts");     //The spliting control     var split_control = new OpenLayers.Control.Split({            layer: contract_layer,            listeners: {                'beforesplit' : function(){                    alert("before");                }            }     });       and then I add this control to my map and make it active:               toolbarItems.push(new GeoExt.Action({            control: split_control,            map: map,            toggleGroup: 'map'        }));             map.addControl(split_control);        split_control.activate();      And now when I chose the tool, I can draw a line which ends with double-click, but that's it. nothing else.Nothing splits and nothing is sent to geoserver. What do you think I'm missing?


Popup for ArcGISDynamicMapServiceLayer point data

I was going through this example from ESRI Resource popup. I was able to integrate it to my polygon layer but it is failing for point and line layer.I am not sure but               query.spatialRelationship = esri.tasks.Query.SPATIAL_REL_INTERSECTS;      could be the reason.I tried the other combination of the task but was not successful if anyone has another example or if anyone could suggest something in existing example then it will be very helpful.I am using ArcGIS server 10.2,Javascript API 3.1.


GDAL - projection from pixels to geo coordinates

I am using GDAL library in C++. On my dataset, I have obtained geo transform with `GetGeoTransform` to GT array (GCP is NULL). Than I use this values for this equation (Source):               Xgeo = GT(0) + Xpixel*GT(1) + Yline*GT(2)     Ygeo = GT(3) + Xpixel*GT(4) + Yline*GT(5)      I also add 0.5 offset to `Xpixel` and `Yline`, to calculate coordinate in the "center" of pixel. If I put [0,0], I got correct result in lat/lon. However, if I change `Yline` to something other than 0, my longitude is messed up - instead of 44.0, I got -349.2... (latitude is still correct). What am I doing wrong ? **My data informations** \- origin, pixel size and top left are taken from GT array used in equation above; projection is `GetProjectionRef`               Projection is `GEOGCS["Coordinate System imported from GRIB file",DATUM["unknown",SPHEROID["Sphere",6371200,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]]'     Origin = (0,0) - GT(2), GT(4)     Pixel Size = (0.5,-33.268) - GT(1), GT(5)     Top Left = (5.75,66.634) - GT(0), GT(3)


Using a value from "iterate feature selection' to calculate a field

I'm trying to streamline a model that identifies all zones that intersect each parcel in our Local Government Area (eg. Parcel x has Zone B1, B2 and E4). I am trying to make use of iterators in ArcGIS10. The problem is in bold italics below; I have a Zones layer and a Parcels layer. I select all polygons of Zone A (using iterate feature selection) then do a select by location to get parcels that intersect those selected polygons. I now want to check if Field1 is empty, if it is, **_use calculate field to populate Field1 with the "Zone A" value used to select my zone in the first place._** I then continue with some simple if-then-else logic and the model then iterates with Zone B. ![Screen Capture of my model](http://i.stack.imgur.com/uMeNC.jpg) I built a working version of this in ArcGIS 9.3.1 that works but it is extremely cumbersome and making small changes to the process is very time-consuming. Any suggestions are much appreciated. P.S. I have no python experience and only a little VB knowledge but if writing an arcpy script is the only solution to my problem then I'm more than happy to start learning!


TINs: What advantage does unequal sampling of data points have over regular sampling?

Triangulated Irregular Networks. Under what conditions is one sampling method more appropriate over the other?


Split edges at intersections in PostGIS

I'm attempting to do a walking distance analysis using OpenStreetMap data (osm2pgrouting) and PostGIS and following underdark's blog posts, however I need to exclude residential and service roads that cross primary and secondary roads, as they are unsafe to cross. To do this I think my first step is to split edges at intersections. I'm guessing I'll want to use the ST_Split function, but I'm not sure how to put it all together in one SELECT statement. (A later question will be how to remove such crossing edges based on the highway=* tag, but I need to figure this one out first.)


QGIS 1.7 + MacOS 10.6 : georeferencing from map canvas bugs

I am experiencing issues with the gdal plugin to georeference a raster. It works fine when I enter the coordinates manually but it does not work when I try to link the plugin with the map canvas. Does anyone had this issue and was able to solve it? Thanks Guillaume QGIS 1.7 MacOS 10.6.7


Newbie PostGIS Geometry and Multilinestring Clarification

I'm having some trouble figuring out if I've got the right data loaded for PgRouting. Using the bash script included in PostGIS 2.0, I loaded Tiger2010 data for U.S. State California.   1. The edges table contains the_geom column, data type geometry. Using Underdark's example, it seems I need the roads in a Multilinestring format in order to begin generating routes. Can the data in the edges table be converted into type multilinestring? The table query is below.   2. I altered my geocoded location data to data type geography. PgRouting's shortest_path function needs the data in integer format. How do I convert the geography type point into an integer that shortest_path can use? Thanks for your patience.               CREATE TABLE tiger.edges     (       gid integer NOT NULL DEFAULT nextval('edges_gid_seq'::regclass),       statefp character varying(2),       countyfp character varying(3),       tlid numeric(10,0),       tfidl numeric(10,0),       tfidr numeric(10,0),       mtfcc character varying(5),       fullname character varying(100),       smid character varying(22),       lfromadd character varying(12),       ltoadd character varying(12),       rfromadd character varying(12),       rtoadd character varying(12),       zipl character varying(5),       zipr character varying(5),       featcat character varying(1),       hydroflg character varying(1),       railflg character varying(1),       roadflg character varying(1),       olfflg character varying(1),       passflg character varying(1),       divroad character varying(1),       exttyp character varying(1),       ttyp character varying(1),       deckedroad character varying(1),       artpath character varying(1),       persist character varying(1),       gcseflg character varying(1),       offsetl character varying(1),       offsetr character varying(1),       tnidf numeric(10,0),       tnidt numeric(10,0),       the_geom geometry,       CONSTRAINT edges_pkey PRIMARY KEY (gid),       CONSTRAINT enforce_dims_the_geom CHECK (st_ndims(the_geom) = 2),       CONSTRAINT enforce_geotype_the_geom CHECK (geometrytype(the_geom) = 'MULTILINESTRING'::text OR the_geom IS NULL),       CONSTRAINT enforce_srid_the_geom CHECK (st_srid(the_geom) = 4269)     )     WITH (       OIDS=FALSE     );     ALTER TABLE tiger.edges OWNER TO postgres;     GRANT ALL ON TABLE tiger.edges TO postgres;     GRANT ALL ON TABLE tiger.edges TO gis_group;          -- Index: tiger.idx_edges_tlid          -- DROP INDEX tiger.idx_edges_tlid;          CREATE INDEX idx_edges_tlid       ON tiger.edges       USING btree       (tlid);          -- Index: tiger.idx_tiger_edges_countyfp          -- DROP INDEX tiger.idx_tiger_edges_countyfp;          CREATE INDEX idx_tiger_edges_countyfp       ON tiger.edges       USING btree       (countyfp);          -- Index: tiger.idx_tiger_edges_tfidl          -- DROP INDEX tiger.idx_tiger_edges_tfidl;          CREATE INDEX idx_tiger_edges_tfidl       ON tiger.edges       USING btree       (tfidl);          -- Index: tiger.idx_tiger_edges_tfidr          -- DROP INDEX tiger.idx_tiger_edges_tfidr;          CREATE INDEX idx_tiger_edges_tfidr       ON tiger.edges       USING btree       (tfidr);          -- Index: tiger.tiger_edges_the_geom_gist          -- DROP INDEX tiger.tiger_edges_the_geom_gist;          CREATE INDEX tiger_edges_the_geom_gist       ON tiger.edges       USING gist       (the_geom);


OpenLayers: how to save added point in server side

* I'm new in OpenLayers.   * I have simple application, which use a picture as map,   * on this picture users can add, move, or delte point, and then save.   * I have implemented adding, moving, deleting functionalities, only except save.   * codes are in below : codes here   * anyone who knows how to save change on server side? save to database, file?


How to validate Personal Geodatabase against XML schema?

Is their any tool in arcgis model builder which can validate a personal geodatabase and a XML schema to give out an report saying if the gdb is validated or not ...


How to use SLD Select control?

I have a WMS layer and i want to get fid of polygon in it by click it. How can i do it?   I know that i can create WFS layer from this WMS and get feature attributes from layer or store. But its possible do this use just WMS? I dont want use WFS Layer couse it may contain ~30k features. But if i will use WFS Layer with visible=false selection control will be works? I find SLDSelect control its can help me?   Maybe its possible to do like in gxp library?


Invalid Geometry Error after ST_ConvexHull

I've a Postgis table of around 80,000 points and want to be able to populate a second table summarizing the data with convex hulls based on points with a common attribute. For example, the data contains the telephone exchange which each point is associated with - I want to create a convex hull for all points in each telephone exchange. It ought to create around 1000 overlapping polygons. I tried using the normally excellent QGis ftools convex hull plugin with the option to create hulls based on an input field but my computer spent 6-hours maximising 3 of the four cores before crashing QGis. And I've tried doing it directly in sql: > INSERT INTO new_table(exchange_name, the_geom)   >  SELECT exchange_name, ST_ConvexHull(ST_Collect(the_geom)) As the_geom   >  FROM first_table   >  GROUP BY exchange_name But this complains about invalid geometry - I've tried geometry constraints of "MULTIPOLYGON" and generic "GEOMETRY" so not sure why this should happen - the source table loads perfectly in QGis. The Grass "v.hull" tool only creates a convex hull for all points and not based on an attribute. Any tips or suggestions would be very welcome!


Why is shapefile layer not visible on WMS layer in QGIS?

In my QGIS project I have added a WMS layer for which the CRS is WGS84 EPSG:4326. Now I want to show my shapefile layer which has the same CRS. QGIS shows the layer separately but does not show with the above WMS layer. Can any one help me to understand why? I am using Quantum GIS 1.7.0 I made up shape file for my project for a city area in which point locations are there.with some attribute data. I have directly drag and drop shape file in q GIS. is it ok or i made some wrong thing?please guide me..


Batch Conversion software for latitude and longitude into shapefile format

I have a huge list of latitude and longitude and I want them to be converted into a shapefile format using some batch processing.


Editing QGIS point layer: retrieve attributes of intersecting line layer

I'm creating a custom editing form using QGIS 2.0 and would like to modify the mouse event when creating a point. I've gotten a start by reading this thread: How to programatically check for a mouse click in QGIS In fact, the code below is copied and modified directly from it, so thank you if you helped me with this code snippet. My goal: When creating a new point, I would like to retreive feature attributes of another (line) layer that intersects with the mouse click, using a rectangle filter. The mouse click opens the form when creating the point, while automatically retrieving data from another line layer intersecting with the rectangle buffer. This is in the interest of comparing data in these two layers that are supposed to correspond. Following the lead of the mentioned thread, i'm attempting to make use of the following code without creating a separate Map Canvas. To put it to use, what would I provide, if anything, as the Python Init function under Layer Properties > Fields? Currently I have mypyfile.formOpen as I've added python logic as detailed here: http://nathanw.net/2011/09/05/qgis-tips-custom- feature-forms-with-python-logic/ I may be going about this in the wrong way, so i'm open to suggestions to direct me towards creating this data collection tool.               class PointTool(QgsMapTool):         def __init__(self, canvas):             QgsMapTool.__init__(self, canvas)             self.canvas = self.iface.mapCanvas()             self.iface = iface             self.pointTool = QgsMapToolEmitPoint(self.canvas)              def run(self):             self.canvas.setMapTool(self.pointTool)             QObject.connect(self.pointTool, SIGNAL("canvasClicked(const QgsPoint, Qt::MouseButton)", self.selectFeature))              def selectFeature(self, point):             #resLayer previously defined             if not resLayer or resLayer.type() != QgsMapLayer.VectorLayer():                 QMessageBox.warning(None, "No!", "Select a vectorlayer")                 return                  width = self.iface.mapCanvas().mapUnitsPerPixel()*2             rect = QgsRectangle(point.x() - width, point.y() - width, point.x() + width, point.y() + width)             rect = self.iface.mapCanvas().mapRenderer().mapToLayerCoordinates(layer, rect)                  resLayer.select([], rect)             feat = QgsFeature()             ids = []             while resLayer.nextFeature(feat):                 ids.append(feat.id())             #selecting the feature or features that intersect with mouse click buffer             resLayer.setSelectedFeatures(ids)          #code continues to guide user data entry using python logic


what is the relevance of cosine in determining map resolution?

"At any latitude and zoom level, you can determine the scale by using the following equation"               Map resolution = 156543.04 meters/pixel * cos(latitude) / (2 ^ zoomlevel)      is a quote from (http://msdn.microsoft.com/en-us/library/aa940990.aspx) Can anyone explain to me the relevance of the cos(latitude) part? I am having difficulties visualising this.


IToolbar uncheck items - ArcGIS

I currently have an app with two menu bars. One being the AxToolbarControl in a UserControl and one in a separate menu bar containing in-house stuff. I currently have the issue where if I select something in the AxToolBarControl and then select something in the other menu-bar the item in the AxControl is still checked. This obviously results in two commands being executed on the mousedown event handler. I've looked at the IControl and IToolbar but it looks like .Checked is a ReadOnly. Any ideas?


How to collect data via WFS and Popup

I have a map with a WFS-Layer and I am able to add Polygons to my data stored in a PostGIS-database, served with Geoserver (Openlayers, Geoserver 2.1.4, PostGIS 2.0 for PostgreSQL 9.1).) How can I add a description to the polygon? E.g., after drawing a polygon I want to specify the crop that are grown there, i.e. insert information into the Postgis-table-column "crop". I found this example: http://gis.ibbeck.de/ginfo/apps/OLExamples/OL26/examples/styles_unique_with_group_wfs.html. The Popup type "edit name and description" would be perfect for my means, unfortunately I can't work out how it works as i have a really hard time understanding its code, greenhorn that I am...   * How is it done?   * Or is there an easier example/ way of how to use a popup to edit your attributes? Thanks so much, I appreciate any help I can get! * * * Here's my WFS-Layer and the tools used so far. The desired column to edit/add Information in the postgis-table "test" is "crop".               function init() {     Save-strategy     var saveStrategy = new OpenLayers.Strategy.Save();          //empty map, bounds are test-layer bounds (EPSG:32647)       map = new OpenLayers.Map({         div: "map",         allOverlays: true,         maxExtent: new OpenLayers.Bounds(         653237.69439077,1519879.063165,655229.57939001,1520825.6733868         )         });          //WFS-Layer Test= editable data     var test = new OpenLayers.Layer.Vector("Editable Features", {         strategies: [new OpenLayers.Strategy.Fixed(), saveStrategy],             protocol: new OpenLayers.Protocol.WFS({         url: "http://..../wfs",         featurePrefix: 'testkf',         featureNS: "http://.../testkf",         extractAttributes: true,         featureType: "test",         geometryName: "geom",         })     });          map.addLayer(test);          //Toolbar:         var panel = new OpenLayers.Control.Panel(         {'displayClass': 'customEditingToolbar'}         );              var draw = new OpenLayers.Control.DrawFeature(         test, OpenLayers.Handler.Polygon,         {         title: "Draw Feature",         displayClass: "olControlDrawFeaturePolygon",         multi: true         });              var save = new OpenLayers.Control.Button({         title: "Save Changes",         trigger: function() {         if(edit.feature) {         edit.selectControl.unselectAll();         }         saveStrategy.save();         },         displayClass: "olControlSaveFeatures"         });               panel.addControls([save, draw]);     map.addControl(panel);       map.addControl(new OpenLayers.Control.LayerSwitcher());     map.zoomToMaxExtent();     }      EDIT: Hi, so nobody knows anything about this topic? Or is there another problem with my question? I found another example, however, this is way too complex for me to understand yet: http://dev4.mapgears.com/bdga/bdgaWFS-T.html This is exactly what I am looking for. Can anyone break it down to me what the steps are to achieve this? Or is there any way to improve my question? Thanks for your help!