Trying to make QT + QGIS working

I'm mostly a Java programmer, I worked on C++ many years ago, and knowing this my project manager has asked me to create a little app in Qt that can open a shape file, using QGIS. But I'm encountering some difficulties. I have downloaded the QGIS tutorials of Tim Sutton from GitHub (https://github.com/qgis/QGIS-Code-Examples) and I began working on this: 4_adding_rasters_to_canvas. I have downloaded QGIS with all the sources, version 1.8, and Qt4 (first I tried QGIS 2.0 and QT5, but I have read that there is no compatibility, so I have tried to reproduce the original environment of the tutorials). I have open the project, edited the .pro file to include folders and libraries. I have downloaded OSGEO4W64 (I am on Win7) with all the sources. I have configured the kit under Qt4, using as a compiler Microsoft Visual C++ 10.0 (the only one that doesn't make appear the red error icon in the kits list). I have configured CMake with the path to the CMake 2.8 executable. When I build the project, I see the following errors:               main.obj:-1: error: LNK2019: unresolved external symbol "public: virtual __thiscall QgsApplication::~QgsApplication(void)" (??1QgsApplication@@UAE@XZ) referenced in function _main     main.obj:-1: error: LNK2019: unresolved external symbol "public: __thiscall QgsApplication::QgsApplication(int &,char * *,bool,class QString)" (??0QgsApplication@@QAE@AAHPAPAD_NVQString@@@Z) referenced in function _main      I have done many things for solve this issues, but I want to start again from this point to be sure. My (first for now) question is: I need to compile the QGIS sources under Visual C++ before building the project, because Qt can't see the QGIS dll? Thanks for your time. Bye, Ivan A. **UPDATE** Ok, it seems that I've managed to solve the issue: I've recompiled the 1.8 sources with cmake-gui, and linked the libraries in the build folder. In this way QT sees the libraries. But now I have another issue. There is this code in the tutorial               QgsMapLayerRegistry::instance()->addMapLayer(mypLayer, TRUE);      (mypLayer is a QgsRasterLayer*) and when building it causes the following error:               error: LNK2019: unresolved external symbol "public: class QgsMapLayer * __thiscall QgsMapLayerRegistry::addMapLayer(class QgsMapLayer *,bool,bool)"      I have read that the addMapLayer method is deprecated, so I have replaced it with the method addMapLayers (according to Tim Sutton: http://osgeo- org.1560.x6.nabble.com/Display-raster-layer-on-own-Map-Canvas-td5002935.html ). So I used:               QgsMapLayerRegistry::instance()->addMapLayers(QList<QgsMapLayer*>() << mypLayer);      and the error slightly change but doesn't go away...               error: LNK2019: unresolved external symbol "public: class QList<class QgsMapLayer *> __thiscall QgsMapLayerRegistry::addMapLayers(class QList<class QgsMapLayer *>,bool,bool)"      It's the proper syntax (I found examples on the internet)? There is another way to do this?


Super large scale maps with Leaflet, how?

I have a map/architectural drawing of two floors of a huge building with a lot of rooms. Now I'd like to make an interactive slippy map using preferably Leaflet showing all the rooms and use markers/polygons showing what happens in each room and who's sitting where, etc. These markers should quite easily be editable. Does anyone have any ideas of how to make this happen? First I was thinking about just pretending house covers the whole world and then use standard web mercator with `generate_tiles.py`. Another idea is to use zoomify together with Leaflet. The problem with both methods is that it's difficult to get right coordinates. Best would be to use some kind of a local projection.


PGRouting 2.0 for RHEL5

Is there a PGRouting 2.0 rpm package available for Redhat Enterprise Linux 5.8 and has anyone managed to install it successfully?


LINECAP and PATTERN on PHP mapscript

I am trying to create a railroad line pattern by using LINECAP and PATTERN. The manual says it can be done by using both paramaters on STYLE inside mapfile. Unfortunately, i couldn't find any method or option to do that on php mapscript. Closest i can find is by using setpattern(). But that method is inside symbolObj. Any idea how to do LINECAP and PATTERN with php mapscript ? Thanks. :)


osm2pgsql: Import data in custom projection

I have osm-file and I want to import it to the postgresql via osm2pgsql with custom projection EPSG:20022. I look to the help and found out any keys like -E EPSG: but when I try to do with it I get the data in SRID=4326. Is there a way how to fix it?


Geodjango insert statement: psycopg2.IntegrityError: violates check constraint

Why is the following INSERT statement throwing this error?               INSERT INTO domes_place (placeid, structidx, oscode) VALUES ('10', '1', 'POINT(51.121282945772 1.3138398501222)');     psycopg2.IntegrityError: new row for relation "domes_place" violates check constraint "enforce_srid_oscode"      This is the table definition in postgres:                          Table "public.domes_place"        Column   |          Type          | Modifiers      ------------+------------------------+-----------      placeid    | integer                | not null      structidx  | integer                |        oscode     | geometry               |      Check constraints:         "enforce_dims_oscode" CHECK (st_ndims(oscode) = 2)         "enforce_geotype_oscode" CHECK (geometrytype(oscode) = 'POINT'::text OR oscode IS NULL)         "enforce_srid_oscode" CHECK (st_srid(oscode) = 4326)      Isn't this POINT in WGS84 - i.e. in SRID 4326?


How do I create a static popup that displays a local web page?

I am using OpenLayers 2.11 and would like to place a marker on my map which, when clicked on, would open a popup which then loads a local web page. I would like the popup to be located in the lower, left hand corner of the screen regardless of zoom level or map scrolling by the user. I have been unable to locate an example of how to perform this. Any assistance is appreciated. Thank you!


GeoServer raster tiles are individually contrast-balanced: how to turn this off?

I'm running GeoServer locally as a servlet (in JBoss), and I'm serving a GeoTIFF file through the WMS protocol. In a webpage running an OpenLayers map, I have the GeoTIFF-WMS layer overlaid above an OpenStreetMaps base layer. The problem is that it looks like this: ![GeoTIFF over OpenStreetMap](http://i.stack.imgur.com/mnfR5.png) The GeoTIFF is clearly being tiled (good) and each tile is normalized to a given brightness level (bad). When I zoom in, I can see a similar effect taking place on the smaller tiles, and I know that the original GeoTIFF doesn't have this blockiness: it's a diagonal strip of 16-bit greyscale surrounded by zeros that ought to be transparent. How do I turn off this brightness-normalization feature (or what is it called, so that I can Google it)? Even better, how do I tell GeoServer that pixels with zero values ought to be transparent?


Mapping Trails with GPS - Hardware Suggestions

I have been hired to map about 100km of XC ski trails. I am super proficient in GIS but have not done any detailed GPS work since college when I worked with some sort of Trimble explorers. This is a complex network of trails and each segment needs to be a unique feature. I own a Garmin etrex 30. My guess is that this will not be a great unit for this sort of work. I am thinking of buying a used Trimble GeoExplorer and an antenna that I can attach to my hat so that the GPS can be in a hip pack. I do want my polylines to be accurate and if I remember correctly the Trimble points can be nicely corrected after collection. Any other suggestions for hardware? Do I even need a device this sophisticated?


How can I split polygons in the Natural Earth 10m land shapefile?

This problem, that I have also been struggling with, in QGIS 1.7.3 as well. Here are the steps to reproduce the problem: **1)**Download Natural Earth 10m land shapefile from original website. * _2)_ *Open the file as vector shapefile in QGIS. * _3)_ *On layer browser, right click and enable "Toggle Editing" * _4)_ *Press button "Select single feature by click", then click on the piece of land (becomes yellow) * _5)_ *Click "Split Feature" scissors button. At this point the land layer disappears, (the map display goes white, but land layer is still on layer browser though) and the "Toggle Editing" is disabled without me doing anything. As soon as I re-enable editing, the map comes back. Other features like drawing polygons function normally. Thanks in advance.


Using iface.showAttributeTable() on QgsVectorLayer retrieved by getObject()

I'm unsure if i'm properly using the processing.getObject() function. I'm creating a python processing script for QGIS 2.0 and 2.2, but can't seem to get an attribute table to open programmatically. The code below works fine in the Python Console, but causes QGIS to crash when I tried as a Processing Script. The only difference was that the layer is called by QgsVectorLayer() in the Python Console and processing.getObject() in the Processing script, thus my reasoning behind my question.               from qgis.core import *     from PyQt4.QtCore import *     from PyQt4.QtGui import *     from processing.core.VectorWriter import VectorWriter     import numpy     import processing     import csv     import qgis.utils          #Here we define the input and outputs     #====================================     ##[Example scripts]=group     ##Couche_compteurs_individuels=vector     ##Couche_compteurs_generaux=vector     ##Couche_canalisations_sur_compteur=vector          #And here is the body of the algorithm     #=======================================          #input layers values are always a string with its location.     #That string can be converted into a QGIS object (a QgsVectorLayer in this case))     #using the processing.getobject() method     indivLayer = processing.getobject(Couche_compteurs_individuels)     grosLayer = processing.getobject(Couche_compteurs_generaux)     resLayer = processing.getobject(Couche_canalisations_sur_compteur)          CG_IDlist = []          def AEPstats():         Gros_Feat = grosLayer.getFeatures()         grosIDselection = []         grosidNull = False         for feat in Gros_Feat:             ID_idx = grosLayer.fieldNameIndex('ID_CG')             #selection list of null values             if feat.attributes()[ID_idx] == NULL:                 grosIDselection.append(feat.id())                 grosLayer.setSelectedFeatures(grosIDselection)                 grosidNull = True             #Putting existing attributes in list             else:                 CG_IDlist.append(feat.attributes()[ID_idx])         #Open msg box, attribute table with selected items, and end function         if grosidNull == True:             QMessageBox.warning(                 None,                 'Probleme de saisie',                 "Completez identifiant(s) de reference ('CG_ID') de Compteurs Individuel")             qgis.utils.iface.showAttributeTable(grosLayer)             return      Thanks in advance


Convert vector data with irregular points (no grid) to raster with QGIS?

I am aware that how to convert a vector layer to a raster layer is a quite frequent problem and that there are many questions and answers dealing with it. However, none was helpful concerning my specific problem, so I hope it's okay to ask again. I have a huge dataset consisting of fish catch data - different species and the location where they were caught. The original data is a csv file, but I'm working with it in QGIS and have created a vector file from the spreadsheet. Now I want to use another program (Zonation) that only accepts raster files to calculate species richness in certain areas - is there a possibility to create a raster from a vector layer with irregular points? I have tried the Rasterize tool, but it only gives me a solid rectangle and all individual data points are lost. I'm working with CRS WGS 84 and I have tried to set the cell size to 0.4, but still it's only one solid square. I've never worked with raster files before and I'm very grateful for your help - thanks!


Arcgis CopyFeatures tool extremely slow when exporting to SDE

I have about 10 points in a file gdb that I am attempting to copy over to an sde feature set. In arcpy and in model builder it takes about a minute to execute. Do you know of any tips to speed this up? In SSMS I can query the table almost instantaneously, so there must be a bottleneck somewhere that I am not aware of, since the network is close to zero latency. In the progress bar, it reads "Writing Features 1" for the majority of the time, then seems to complete the rest of the task almost immediately. If it is really supposed to be this slow, I would at least like to have some idea of what it is actually doing during this time. The file is only a few kilobytes total, so I assume this cannot be normal. **Update** : I was trying this operation over a site-to-site VPN. It is a fast connection, so there is still no real reason it should take a full minute to upload 10 points. I tried this on a local connection (local to the SQL Server), and the same operation took 2 seconds. It seems that CopyFeatures has some sort of bug when working through a VPN connection. Any ideas on how I could find a workaround (I need to be able to copyfeatures over the network with a VPN connection)?


Windows 7 ArcGIS 9.3 with TOPO!Pro Runtime Error

A new, clean install of ArcGIS on Windows 7 SP1 and National Geographic TOPO! Pro Suite -runtime error "Runtime Errror (UpdateGroupLayer)#-2147217384 Failed to open raster dataset" TOPO! Pro was designed with XP not Win 7. unable to get Tech Support (closed last year - bummer) Help the fireguy please!!


Setting ProductCode for python script (Citrix environment)

I have set up a python script in ArcToolbox for our users. However when running it with ArcView 10.0 licensing the script fails with: > ERROR 000824: The tool is not licensed I think it's an Identity_analysis command triggering this error. I've googled the topic and it seems there was a SetProduct(code) command, but I can't seem to get it to work. Can someone provide the syntax to get this running, or is it a limitation running ArcGIS in a Citrix environment?


Is a proxy required to use OpenLayer's GetFeatureInfo?

My domain is hosted on HostGator and apparently they don't allow the use of proxies. I have a map that has a WMS point feature hosted from my GeoServer using OpenLayers. I want to be able to select a point and get information from it. Can I do this without using a proxy?


What is the relationship between Map Service Instances and SDE connections?

We have a GPS project based on ArcGIS Mobile 10.0. There are two web map services and each were kept as Nonpooled with Minimum instances 15 and maximum instances 60 by original developer. The database is ORACLE and connection is made through SDE. Recently whole enterprise moved to 10.1 except GPS project which has some crashing issues with 10.1. But that led hosting of this GPS project on an old Test server. Two incidents happened in last two weeks when the connections for GPS projects reached alarming high in few minutes (122 and 58). First time it affected the SDE server and caused total chaos- we had to kill the processes and bounce the SDE server back. Second time, we restarted the services before things out of control. I investigated the first time - only 20 GPS devices were connected that time. I am new to ArcGIS mobile and don't know SDE in detail(It is still under DBA control and I can only investigate the appserver). How will the number of SDE connections being used one time transforms to map services? I want to go to management with the reasons and recommendations.


Interpolation outside the contour of a set of points

I have a field contour in a shape layer. I have another layer compound by a set of points which are the locations where soils samples have been taken within the field. I have in the attribute table of this second shape layer the Phosphorus content of the soil in those points. As I have said, the contour of the set of points is smaller than the field contour, being the contour of the set of points enclosed within the field contour. I have interpolated this set of points by creating a TIN in order to have a value of Phosphorus content for any point in the field. But through this procedure, the area between the contour of the set of points and the contour of the field is excluded from the TIN. I have also tried the Kriging method (this method allows to enter the field contour shape as the Processing Extent in the Environments section, being the referred problem solved) but through this interpolation method a new problem arises: the sample points are lying over classes which value ranges don’t include the real values of the sample points. Then, the original data is in some way adulterated. I guess the IDW method isn’t good for this purpose: there isn’t any point in the field contour, and thus the value will go down towards the field contour in an unreal way due to this fact. Any idea or solution for this issue (interpolating the area excluded from the TIN)?


Rasterize difference between QGIS and FWTools

I am wanting to use FWTools to rasterize a large number of Shapefile polygons. I have used the Rasterize tool in QIS and it works great. But when I use the generated code in FWTools the process does not work. Here is the generated code:               gdal_rasterize -a COVTYPE -ts 5000 3705 -l LCC001k_dis     C:/GlobalData/Vectors/LCC/LCCRasters/LCC001k_dis.shp     C:/GlobalData/Vectors/LCC/LCCRasters/LCC001k.tif      Unfortunately I am not a programmer, so I am not able to create my own script to perform this task. Can anyone assist me with why this code doesn't work in FWTools, or what changes I need to make to utilize the Rasterize functionality.


Transparency in ArcMap Graphic Elements

Is it possible to apply transparency to graphic elements in an ArcMap layout? I'm specifically refering to a neatline -- a box containing my legend, source info, scale bar, etc.


Workflow process to style vector data using GeoServer

I'm starting to get a better understanding of map-servers in general, but I am still unclear on how best to use GeoServer. At the moment, I use web-hosting to store the web-browser related parts of my site, and `php` code to connect to a postgresql/postgis database. All of my data is stored in a desktop server, which I chose to use as it is a very low cost way of storing a large amount of data.               1] web-browser <--> database         This works fine also, as I do not need to style vector data; I am only returning numerical results with geojson vector overlays. However, I would like to incorporate GeoServer into this workflow so that I can apply better cartographic styles to data in the future. I understand the work flow should be similar to the following:               2] web-browser <-->  GeoServer <--> database      My difficulty in understanding this piece, is how to connect GeoServer into the process. Can I run GeoServer on my desktop and connect to it to the browser via `php`, similar to the database approach? Or does GeoServer need to be installed on the web-host and managed there?


Fast Offline OSM Routing Library to Replace MoNav (.NET)

I'm looking for a .NET replacement to MoNav. I'm currently calling a C++ MoNav DLL via COM and it works great--it's fast and memory efficient, and can both load itself and route across the country almost instantaneously. A pre- processed graph for a 5GB OSM file is only ~400MB, which is another plus as I do not need to distribute the source OSM files. However, the COM dependency is a slight issue with certain customers' setups. Its turn-by-turn directions leave a lot to be desired as well. So I was wondering if anyone knows of a .NET replacement? One that preferably used Contraction Hierarchies, has a low memory footprint, and can reliably route across the United States in under a second. That may be asking a lot, I know. My searches have only turned up OsmSharp thus far, and while it looked promising initially, in practice the load times were abysmal for even pre- processed CH graphs (though I am still playing with it and tweaking parameters to try and resolve this). Just wanted to know if anyone knew of a suitable replacement before I attempt to write a C++/CLI wrapper for the original MoNav library and spend a few days (weeks?) on a better turn-by-turn system. So far I haven't found one. Thanks!


maintain different style for each feature in same layer

on feature selection, I am trying to change feature style of selected feature.so that assign different style to each selected feature that exist in same layer. I have done this successfully but when I zoomed the map, feature redraw and assign the last symbology to all selected features. Please resolve my issue Here is my code: var style = OpenLayers.Util.extend({}, OpenLayers.Feature.Vector.style['default']);                   style.strokeColor = color;        style.strokeWidth = width;        selectItem.style = style;        polygonLayer.drawFeature(selectItem);


Edit a classified grid file (.grc) to improve accuracy

I have a .grc file (classified grid file) which contains some clutter information, forest, sea, woodland, etc... I need to alter the data in some area to improve accuracy. Basically a larger forest or possibly create a new value like "tall tree forest", etc.... Is there a way to do this? Thanks in advance


Cannot install node-mapnik on Windows 2008 x64

I have been spending the last few days trying to get node-mapnik running but with no success. Initially inspired by Wilhelm Berg's demo (https://www.youtube.com/watch?v=GlTaRpg3CbQ) which renders pure vector maps using node-mapnik I was interested in getting this running on a Windows Server. I have worked through all of the StackExchange issues on these topics including Wilhelm's checklist for getting this working on Windows 7 x64 but I am at a dead end and my knowledge is really limited in being able to further diagnose the problem. My environment is:   * Windows 2008 Server R2 x64   * Visual Studio Express 2012 Desktop Edition (32 & 64 bit)   * Node.js v0.10.26 (64 bit)   * Python27 (32 bit)    I have taken the following steps:   1. Installed Mapnik 2.2.0 SDK downloaded from git and placed in C:\apps\mapnik-v2.2.0      2. Added to PATH: C:\apps\mapnik-v2.2.0\lib C:\apps\mapnik-v2.2.0\bin      3. Installed Python v2.7 (32bit) downloaded and placed in C:\apps\python27      4. Added to PATH: C:\apps\python27 C:\apps\python27\scripts      5. Added to PYTHONPATH: C:\apps\Python27\Lib; C:\apps\Python27\DLLs; C:\apps\Python27\Lib\lib-tk; C:\apps\mapnik-v2.2.0\python\2.7\site-packages;      6. Rebuild node-gyp using VS 2012 node-gyp rebuild --msvs_version=2012      7. Downloaded and installed protobuf (32 bit) and placed in c:\apps\protobuf      8. Added to PATH: c:\apps\protobuf      9. Installed protobuf using NPM npm install protobuf --msvs_version=2012      10. Downloaded Cairo using git git clone git://anongit.freedesktop.org/git/cairo      11. Copy Cairo files to mapnik: Copy cairo\src\cairo-svg-surface.c to C:\apps\mapnik-v2.2.0\include Copy cairo\src\cairo-svg.h to C:\apps\mapnik-v2.2.0\include Copy cairo\src\cairo-svg-surface-private.h to C:\apps\mapnik-v2.2.0\include      12. Make proj4 epsg files available:  SET PROJ_LIB=C:\apps\mapnik-v2.2.0\share\proj      13. Attempt to install node-mapnik npm install node-mapnik --msvs_version=2012    I receive the following errors: > node-pre-gyp ERR! Source compile required: 404 status code downloading > tarball ... > > mapnik_grid.cpp c:\users\administrator\node_modules\node- > mapnik\build\release\obj\global_interm ediate\vector_tile.pb.h(9): fatal > error C1083: Cannot open include file: 'google/protobuf/stubs/common.h': No > such file or directory (C:\Users\Administrator\n ode_modules\node- > mapnik\build\Release\obj\global_intermediate\vector_tile.pb.cc) > [C:\Users\Administrator\node_modules\node-mapnik\build\mapnik.vcxproj] > c:\users\administrator\node_modules\node-mapnik\src\utils.hpp(11): fatal > error C1083: Cannot open include file: 'mapnik/unicode.hpp': No such file or > director y (..\src\mapnik_geometry.cpp) > [C:\Users\Administrator\node_modules\node-mapnik \build\mapnik.vcxproj] > mapnik_layer.cpp mapnik_cairo_surface.cpp > c:\users\administrator\node_modules\node- > mapnik\src\mapnik_featureset.hpp(6): f atal error C1083: Cannot open include > file: 'mapnik/datasource_cache.hpp': No s uch file or directory > (..\src\mapnik_featureset.cpp) [C:\Users\Administrator\no de_modules\node- > mapnik\build\mapnik.vcxproj] mapnik_image.cpp > C:\Users\Administrator\node_modules\node-mapnik\node_modules\mapnik-vector- > tile \src\mapnik3x_compatibility.hpp(4): fatal error C1083: Cannot open > include file : 'mapnik/version.hpp': No such file or directory > (..\src\mapnik_map.cpp) [C:\U sers\Administrator\node_modules\node- > mapnik\build\mapnik.vcxproj]   >  mapnik_datasource.cpp C:\Users\Administrator\node_modules\node- > mapnik\node_modules\mapnik-vector-tile \src\mapnik3x_compatibility.hpp(4): > fatal error C1083: Cannot open include file : 'mapnik/version.hpp': No such > file or directory (..\src\mapnik_color.cpp) [C: > \Users\Administrator\node_modules\node-mapnik\build\mapnik.vcxproj]   >  node_mapnik.cpp ..\src\mapnik_grid_view.cpp(9): fatal error C1083: Cannot > open include file: 'm apnik/grid/grid.hpp': No such file or directory > [C:\Users\Administrator\node_mo dules\node-mapnik\build\mapnik.vcxproj] > mapnik_projection.cpp C:\Users\Administrator\node_modules\node- > mapnik\node_modules\mapnik-vector-tile \src\mapnik3x_compatibility.hpp(4): > fatal error C1083: Cannot open include file : 'mapnik/version.hpp': No such > file or directory (..\src\mapnik_layer.cpp) [C: > \Users\Administrator\node_modules\node-mapnik\build\mapnik.vcxproj]   >  mapnik_feature.cpp ..\src\mapnik_grid.cpp(9): fatal error C1083: Cannot > open include file: 'mapnik /version.hpp': No such file or directory > [C:\Users\Administrator\node_modules\n ode-mapnik\build\mapnik.vcxproj] > mapnik_palette.cpp ..\src\mapnik_image_view.cpp(11): fatal error C1083: > Cannot open include file: 'mapnik/color.hpp': No such file or directory > [C:\Users\Administrator\node_modu les\node-mapnik\build\mapnik.vcxproj] > mapnik_vector_tile.cpp c:\users\administrator\node_modules\node- > mapnik\src\utils.hpp(11): fatal error C1083: Cannot open include file: > 'mapnik/unicode.hpp': No such file or director y > (..\src\mapnik_cairo_surface.cpp) > [C:\Users\Administrator\node_modules\node-m apnik\build\mapnik.vcxproj] > mapnik_expression.cpp   >  mapnik_memory_datasource.cpp ..\src\mapnik_memory_datasource.cpp(3): fatal > error C1083: Cannot open include file: 'mapnik/version.hpp': No such file or > directory [C:\Users\Administrator\n ode_modules\node- > mapnik\build\mapnik.vcxproj] c:\users\administrator\node_modules\node- > mapnik\src\utils.hpp(11): fatal error C1083: Cannot open include file: > 'mapnik/unicode.hpp': No such file or director y (..\src\mapnik_feature.cpp) > [C:\Users\Administrator\node_modules\node-mapnik\ build\mapnik.vcxproj] > ..\src\mapnik_image.cpp(7): fatal error C1083: Cannot open include file: > 'mapni k/color.hpp': No such file or directory > [C:\Users\Administrator\node_modules\no de-mapnik\build\mapnik.vcxproj] > C:\Users\Administrator\node_modules\node-mapnik\node_modules\mapnik-vector- > tile \src\mapnik3x_compatibility.hpp(4): fatal error C1083: Cannot open > include file : 'mapnik/version.hpp': No such file or directory > (..\src\mapnik_datasource.cpp ) [C:\Users\Administrator\node_modules\node- > mapnik\build\mapnik.vcxproj] C:\Users\Administrator\node_modules\node- > mapnik\node_modules\mapnik-vector-tile \src\mapnik3x_compatibility.hpp(4): > fatal error C1083: Cannot open include file : 'mapnik/version.hpp': No such > file or directory (..\src\mapnik_projection.cpp ) > [C:\Users\Administrator\node_modules\node-mapnik\build\mapnik.vcxproj] > C:\Users\Administrator\node_modules\node-mapnik\node_modules\mapnik-vector- > tile \src\mapnik3x_compatibility.hpp(4): fatal error C1083: Cannot open > include file : 'mapnik/version.hpp': No such file or directory > (..\src\mapnik_palette.cpp) [ C:\Users\Administrator\node_modules\node- > mapnik\build\mapnik.vcxproj] C:\Users\Administrator\node_modules\node- > mapnik\node_modules\mapnik-vector-tile \src\mapnik3x_compatibility.hpp(4): > fatal error C1083: Cannot open include file : 'mapnik/version.hpp': No such > file or directory (..\src\node_mapnik.cpp) [C:\ > Users\Administrator\node_modules\node-mapnik\build\mapnik.vcxproj] > c:\users\administrator\node_modules\node-mapnik\src\utils.hpp(11): fatal > error C1083: Cannot open include file: 'mapnik/unicode.hpp': No such file or > director y (..\src\mapnik_vector_tile.cpp) > [C:\Users\Administrator\node_modules\node-map nik\build\mapnik.vcxproj]


how to calculate the average temperature of each band of a .grd file?

I am having a GRD file in which the temperature at different longitude and latitude are given in a 32x35 array for one day(band),i am having 365 such bands and want to calculate the average of all individual bands.


Calculation to convert db data to Lat/Lng degrees

In a database I have the following Lat/Lng data. I wish to convert it to a readable format.               177514245:34726053  =>  S 32 7.3947, E 115 51.4245     177513828:34726948  =>  S 32 7.3052, E 115 51.3828     177508930:34829424  =>  S 31 57.0576, E 115 50.8930     177508084:34829842  =>  S 31 57.0158, E 115 50.8084     177504234:34831118  =>  S 31 56.8882, E 115 50.4234     177504474:34831039  =>  S 31 56.8961, E 115 50.4474     177507746:34829955  =>  S 31 57.0045, E 115 50.7746     177504623:34830989  =>  S 31 56.9011, E 115 50.4623     177507014:34830191  =>  S 31 56.9809, E 115 50.7014     177506702:34830295  =>  S 31 56.9705, E 115 50.6702      I was given the formula for converting the db data o the readable format but lost the information after a hhd failure. The number before the ':' becomes the Lng and the second the Lat. I don't want to through any red herrings into confuse anyone, but this is what I can remember. Apparently this is a common format. The raw Lat and Lng number was subtracted by a constant. I can't remember if they were the same, and then divided to give degrees and minutes. But take that all with a grain of salt. Thanks in advance, Simon   1. I have no idea. The data in the db is the 18 varchar figure. The OEM software represents the coordinates in the format as shown above.   2. Have been Google ing it for 3 days and have found nothing. Can't remember the site I posted the question on about 3 or 4 years ago so I could go back and read it again. At the time I was told it was a pretty common method of storing coordinates.   3. The db is sqlite, I use RAD Studio. The OEM GUI that displays the data doesn't allow me to access the data through it. But I can read/write the db itself and I wish to map the coordinates. Hence the need to be able to make sense of how:               177514245:34726053      becomes               S 32 7.3947, E 115 51.4245


How to get for a given pixel the values from different overlapped rasters?

I have a few rasters representing different soil properties, like soil pH and sand content. The rasters overlap perfectly, have the same dimensions and same coordinate reference system. I would like to be able to click on any pixel of my map and get the pH and sand values for this location. The "Identify Features" tool seems to return only the value of the selected raster, pH or Sand. It doesn't give the values of each rasters displayed. Are they any equivalent of the "Identify Features" tool that work on multiple rasters? I am QGIS 2.2.0 user on a Windows computer.


how to include part of a map extent in layout view

I have a map of Scotland and within the map layout I'd like to include Shetland within a box at the top so that I don't have to minimise the full extent of the map - I was wondering what the best way to do this would be. Previously i created a duplicate layer in the data view and zoomed it to only fit in Shetland in the layout view but this seemed a clunky way of doing this? Anyone have any advice?


Google earth incorrectly rendering coordinates?

My circles are rendering funky on the poles, and along the equator. ### Here are my coordinates:               158.05553767833,-78.151355482613     152.68654993135,-76.924955252783     147.72303124605,-75.713526831152     143.06274497048,-74.524895329406     138.63332012974,-73.365493766502     134.38233416968,-72.240817945152     130.27100948566,-71.155701792104     126.27011840531,-70.11448759429     122.35726277966,-69.121132703464     118.51502708319,-68.179276670913     114.72969830909,-67.292283062004     110.9903610723,-66.463264689784     107.28824575653,-65.695097793865     103.61625023098,-64.990428769143     99.968582435602,-64.351675871499     96.340488248312,-63.781027589248     92.728040190375,-63.280438894499     89.127969904661,-62.851626275065     85.537532303713,-62.496062233228     81.954392666042,-62.214969784592     78.376530286024,-62.009317375003     74.802153890611,-61.879814541074     71.229625140461,-61.826908561053     67.657387272327,-61.85078227176     64.083896401248,-61.951353160458     60.507553237337,-62.128273775806     56.926633008616,-62.380933438     53.339211222379,-62.708461163863     49.743082525748,-63.109729656798     46.135669299876,-63.583360142432     42.513915669083,-64.127727755403     38.87416120926,-64.740967096403     35.211986616615,-65.420977473287     31.522020664857,-66.165427203563     27.797693492691,-66.971756167889     24.030914923935,-67.837175533906     20.211647006764,-68.758663165469     16.327325438369,-69.732952609442     12.362061986503,-70.756512568067     8.2955243178175,-71.825512172637     4.1013320161292,-72.935764736584     -0.25528729388528,-74.082638186513     -4.8209969442254,-75.260912557308     -9.6596917885399,-76.464550923909     -14.860832875137,-77.686324272452     -20.55272215938,-78.917181740079     -26.923940644354,-80.145162642771     -34.257840694438,-81.353462558233     -42.985222847202,-82.516927422837     -53.749206371382,-83.595765562688     -67.414565058046,-84.525376429557     -84.75320323954,-85.205991865779     -105.36439351601,-85.51452525361     -126.6349523656,-85.372691897859     -145.37174704566,-84.818624018911     -160.36597119513,-83.969463910769     -172.11759253603,-82.937731923535     178.47360903516,-81.800519081463     170.67339362215,-80.605603613521     163.97799351745,-79.38273857854     158.05553767833,-78.151355482613      ### And here is the result: ![Result](http://i.stack.imgur.com/QskM9.jpg) **What could be causing this and how can can I fix it?**


How to update a CartoDB table at a certain time of day

I have a table in CartoDB 'inland' that I would like to have automatically update each day at 8pm PST. I know that PostgreSQL has no internal function to run an UPDATE COLUMN function at a certain time of day. Is there a way to use Cron or PgAgent to do this? Would I need to setup CartoDB on my local machine instead of only using the webservice?


arcgis python window VS pywin32

I have problems while using the ZonalHistogram function in Python. It seems to give strange results, while, with the same code run in arcgis' python window everything works fine. I am trying to study the snow cover per elevation ranges, thus including in the function a MODIS image (showing snow cover) as the input raster, and a reclissified DEM as the input value raster. I tried to modify the reclassification many times, but still not solving the problem... Moreover, inverting the input raster with the input value raster (meaning, inverting the function parameters), everything works fine, but it's a mess since I'm working on thousands of images and having the MODIS not always the same number of columns (while being the DEM always the same). Hope someone could help me. Here is my script, so that you can have a look at it in order to better understand my request:               import arcpy, os, sys, shutil          newpath= r'G:\\Chile_2012\\MODIS10A2.V005\\provvisoria\\Zonal_Histogram\\temporary'     if not os.path.exists (newpath):         os.makedirs (newpath)          arcpy.env.workspace= 'G:\\Chile_2012\\MODIS10A2.V005\\Cobertura_nival\\Cuencas_Copia\\campione'     arcpy.env.overwriteOutput= True     from arcpy.sa import *          listRaster= arcpy.ListRasters ()     elev_ranges= 'G:\\Chile_2012\\MODIS10A2.V005\\DEM\\elev_ranges'     out_dir= 'G:\\Chile_2012\\MODIS10A2.V005\\provvisoria\\Zonal_Histogram\\'          try:         for raster in listRaster:             arcpy.CheckOutExtension ('Spatial')             ZonalHistogram (raster, 'Value', elev_ranges, newpath+ '\\'+ raster)              arcpy.ClearEnvironment ('workspace')         arcpy.env.workspace= newpath         listTable = arcpy.ListTables ()         print listTable              arcpy.Merge_management (listTable, out_dir+ 'Rio_Aconcagua_Merge.dbf')              shutil.rmtree (newpath)          except:         print arcpy.GetMessages ()      I would add that my idea would be to work with layers (by means of functions that could transform rasters in layers) since the table returned in arcgis has a "LABEL" field (coming from the classification in symbology) which is the thing I would like to include in the final table and which is not included while working out of arcgis (in pywin), in order to know on what I am working on.


Is there an exhaustive, searchable catalog of all GIS web services (e.g. OGC WMS/WFS, REST, SOAP)?

Many of us have struggled to get a list of GIS services that are available on the internet. Specifically map services that expose OGC based interfaces like WMS, WFS or Web Service interfaces like REST, SOAP etc. For instance, just for oceanic and marine data there are several online catalogs such as   1) FDGC   2) NOAA NowCoast   3) IMO Mest   4) NASA If you start with weather, transportation, demographics, agriculture, this list becomes huge! These sites are not easy to use, don't have the most updated info and don't have a good taxonomy/ classification of services towards easy discovery of such services. Is there a better attempt out there to create an exhaustive, searchable catalog of all such services? I would like to know their usage rights, free or pay etc. Uptime, performance etc. would also be nice to know.


Reprojecting Natural Earth raster into natearth projection with gdalwarp causing black squares?

I'm trying to reproject the Natural Earth 2 raster layer (http://www.naturalearthdata.com/downloads/10m-raster-data/10m-natural- earth-2/) into the natural earth projection (http://en.wikipedia.org/wiki/Natural_Earth_projection) using gdalwarp. The whole thing looks pretty good except I get 4 black squares by the edges of the map as seen below. gdalwarp options used [gdalwarp included in QGIS: GDAL 1.11.0, released 2014/04/16]: `gdalwarp -overwrite -s_srs EPSG:4326 -t_srs "+proj=natearth +wktext" -r lanczos -of GTiff NE2_HR_LC_SR_W_DR.tif reprojected.tif ` Any ideas? ![natearth reprojection with black squares](http://i.stack.imgur.com/p5wTs.png)


How to programmatically work with featureclasses from multiple geodatabases?

I have produced a Python script which creates multiple buffers of specific (listed) featureclasses within a geodatabase. This is achieved by running each feature through a `for` loop. I would like to know if/how it would be possible to list featureclasses from more than one geodatabase? Currently I am using `env.workspace` to specify the single geodatabase. Listing the features and then running them through the 'for loop'. A possible problem might be that many of the featureclasses in the various geodatabases have the same name.               import arcpy     import os.path     import tempfile     from arcpy import env          # sets the workspace to my gdb     arcpy.env.workspace = "C:/GIS Home/Mining Features (MATTHEW).gdb"          # Output workspace and temp location is defined     output_workspace = "C:/GIS Home/project_1"     temp_workspace = tempfile.mkdtemp()          # list of features to buffer in gdb     fcList = ["Adit_Metalliferous", "Shaft_Metalliferous", "Quarry_Metalliferous", "Dump_Metalliferous", "Engine_House", "Goffan_Working", "Lode_Surface_Outcrop", "Lode_Undefined_Elevation", "Mine_Level", "Mine_Raise", "Portal", "Reported_Subsidence", "Secured_Feature", "Tunnel"]          #create empty lists for output     l_250 = []     l_500 = []     l_750 = []          for featureClass in fcList:         b250 = os.path.join(temp_workspace, featureClass + "_b250.shp")  # make temp output file         b500 = os.path.join(temp_workspace, featureClass + "_b500.shp")  # make temp output file         b750 = os.path.join(temp_workspace, featureClass + "_b750.shp")  # make temp output file         arcpy.Buffer_analysis(featureClass, b250, 250, "", "", "ALL")   # carry out buffer         arcpy.Buffer_analysis(featureClass, b500, 500, "", "", "ALL")   # carry out buffer         arcpy.Buffer_analysis(featureClass, b750, 750, "", "", "ALL")   # carry out buffer         l_250.append(b250)  # append the buffer to the l_250 list         l_500.append(b500)  # append the buffer to the l_500 list         l_750.append(b750)  # append the buffer to the l_750 list          merge_250 = os.path.join(temp_workspace, "250_merge.shp")     arcpy.Merge_management(l_250, merge_250)     merge_500 = os.path.join(temp_workspace, "500_merge.shp")     arcpy.Merge_management(l_500, merge_500)     merge_750 = os.path.join(temp_workspace, "750_merge.shp")     arcpy.Merge_management(l_750, merge_750)          dissolve_250 = os.path.join(output_workspace, "250_dissolve.shp") # create 250_dissolve.shp for final output     arcpy.Dissolve_management(merge_250, dissolve_250)     dissolve_500 = os.path.join(output_workspace, "500_dissolve.shp") # create 500_dissolve.shp for final output     arcpy.Dissolve_management(merge_500, dissolve_500)     dissolve_750 = os.path.join(output_workspace, "750_dissolve.shp") # create 750_dissolve.shp for final output     arcpy.Dissolve_management(merge_750, dissolve_750)


Creating Legends for QGIS 2.0 Feature Blend Mode Outputs

I am working on visualizing ~4100 buffered point features. The shading of the buffers is keyed to a field in the data. Because the buffers overlap, I turned on multiply in the new feature blend mode options. I like the visual result, but what seems to be missing is an ability to generate a meaningful legend for the resulting image. Are there any workarounds for this? Is the only alternative to venture down the path of polygon intersects and symmetrical difference operations? ![buffered points with multiply feature blend mode](http://i.stack.imgur.com/ajhrQ.pngb)


How to use the raster calculator in the processing modeler of QGIS

I am trying to create a easy way to do raster analysis by building a model for two way raster calculations. I have two raster layers that I want to combine into a single raster. For now I've done this manually by doing the raster calculation myself but for the future I want to be able to input these layers into a model without having to go through the trouble. What is the proper annotation of this formula? How I've tried it is not working. ![enter image description here](http://i.stack.imgur.com/eZDZE.png)


Arcgis Server failed to create site?

Failed to create the site. Admin URL 'http://my.site.url:6080/arcgis/admin' is not reachable from 'LOCALHOST': 'Connection timed out'. Anyone have any idea why i am facing this error? i have installed arcgis 10 server of Red Hat and now i am trying to create site it is giving me this ? i have tried restarting serve but problem won't go away?


Cleaning up data in Excel for use in QGIS

I am trying to clean up a large .txt dataset in Excel. I think that a macro or another application may be able to clean up the data they way I want but I have limited programming skills. I need to create fields from different sections in the .txt file so I can then geocode the points. I am conducting a crime analysis for my local police department and up until now I have been using wildcards and the Find & Replace tool to do this, but I have nearly 2 years worth of police calls. Fields needed:"Address", "City", "State", "Type", "Location", "Date", "Day of Week", "Time Received", "Time Dispatched", "Time Arrived", "Time Cleared". I put up the data from January 2012 on Github in repo called SLO-Crime under my name, chadbunn. I will be using CartoDB to geocode the data and then QGIS to perform the analysis. I will probably make some maps in TileMill as well if I have enough time to. Any input would be appreciated because it is getting tedious! Thanks! https://github.com/chadbunn/SLO-Crime/tree/master Text sample below showing a single incident (plain text file, each section composed of a dispatch header summary followed by type / address / officer and other information):               ===============================================================================                                  120101001 01/01/12 Received:00:05 Dispatched:00:06 Arrived:00:12 Cleared:00:53                                   ===============================================================================                                  Type: MENTAL SUBJ                                             Location:LZ1                                   As Observed:                                           Mental Health                                       Addr: 1257 RICH, San Luis Obispo, CA             Clearance Code:Report Filed in                                       Responsible Officer: Hyman, J                                    Units: 4235  ,4241                                    Des: incid#=120101001 AP/KAVANAGH,KOLLEEN 060455 166 PC clr:RTF oc:MENH call=1l


What is difference between GeometryService areasAndLengths method and Measurement widget in ArcGIS JavaScript API?

What's the difference between the GeometryService areasAndLengths method and the Measurement widget to make a caculation in the ArcGIS JS API? Accuracy?


Calculate population coverage

I am relatively new to QGIS and hope I explain this right. I am looking at the percentage of an population that has geographical access to healthcare. Therefore I have one shapefile with village boundaries including population number of each village. Another shapefile is all Healthcare facilities in the area buffered with 5km zones. I now want to calculate the percentage of the village area that is covered by the Healthcare facility buffers and multiply this by the my population numbers. I did manage to create an intersect shapefile of the two formal shapefiles that created correct new polygons of intersection per village. I now however do not know how to sum the areas (Area_1) of all these buffer polygons who are coded equally in the attributable table (under P_06_ID) and divide this by the total village area (AREA). And than multiply this number by the population (Population). All help welcome! ![enter image description here](http://i.stack.imgur.com/W1cbF.png) ![enter image description here](http://i.stack.imgur.com/EQVqh.png) * * * Thank you so much for the help. The command does indeed calculate the percentage of the population that overlaps with that specific polygon. However it still does not sum the different polygons. ![output after command](http://i.stack.imgur.com/W1cbF.png) So what I now want to do is to sum those numbers with the same parish ID automatically (I do have around 500 parishes). So in this example sum PopPerc from cell 699- 709. Any suggestions?


How to use pgRouting's dijkstra_sp_delta_directed with WHERE clause?

Right now I am using PgRouting for vehicle routing and I have problem with Dijkstra shortest path. My query is:               SELECT * FROM shortest_path('SELECT gid  AS id,           source::int4 AS source,           target::int4 AS target,           cost::float8 AS cost,           reverse_cost::float8 AS reverse_cost     from network ',                                257027,     276521,     true,     true)      But, my table network is very large, so I should use "dijkstra_sp_delta_directed" wrapper function with bounding box to make query faster.               Select gid from dijkstra_sp_delta_directed('network', 5700, 6733, 0.1,true,true);      But I don't know how to write some where clause like "from network where road_class=1" Also, I don't know where and how to write (road_class * cost)::float8 AS cost Thank you in advance.


Updating the value of a field with pre-defined values

I am trying change the value of a row for a field whose type is "short" in arcmap using a python script.               with arcpy.da.UpdateCursor(FILENAME, FIELD) as cursor:              for row in cursor:             row = 'Center'             cursor.updateRow(row)      with this i receive the error "RuntimeError: The value type is incompatible with the field type." I know I am giving a string type to a "short" type field. Is there a way around this? There are pre-set values in this field. Do I need to run some other kind of code to create a list of possible values for this field, then select the value based on some kind of list interaction eg `row = fieldlist[1]` Please ask me to clarify anything, but this is everything that seems relevant to me at this point.


Switching from ArcDesktop to ArcGIS Online

I have many GIS files that I was using to create some maps for my work. Our desktop version of ArcGIS expired and we decided to move to the basic online version. I am unable to upload any of my databases and all the shape files that I have are said to have over 1,000 features. Is there any way I can get this to upload? Maybe splitting the data in two parts and uploading them separately? Any help would be appreciated as all my data that I've been working on is completely useless right now.


How do I generate a random raster (or probabilistic cellular automaton) in QGIS?

I'm trying to create something like a cellular automaton land cover change model using QGIS, in which cells randomly change land cover based on probabilities determined from various input rasters. I've already have a map of probabilities of conversion, and I need to take the next step. I found this question, How do I randomly alter a raster map of habitat types?, which looks very similar to what I want to do. The suggested approach is to use a Markov chain or probabilistic cellular automaton, and the commenter writes: "The transitions can be implemented using Map Algebra statements, which makes this method practicable in any raster-based GIS, even those without direct or quick access to cell-level data." Can anybody tell me how to do probabilistic map algebra in QGIS or one of its extensions? I don't see any probability tools in the Raster Calculator. Failing that, is there a tool that will simply generate random rasters to be later combined using the Raster Calculator?


gdal_retile - Out of range band requested

I'm using gdal_retile to create ImagePyramid from a set of GeoTIFF images. I'm using following command:               gdal_retile.py -v -r bilinear -levels 8 -ps 2048 2048 -co "TILED=YES" -co "COMPRESS=JPEG" -targetDir F:/tiles --optfile F:/tileslist.txt      After couple of hours(source images are 320GB) I get following error message:               Traceback (most recent call last):       File "C:\PROGRA~2\FWTOOL~1.7\bin\gdal_retile.py", line 940, in ?         sys.exit(main(sys.argv))       File "C:\PROGRA~2\FWTOOL~1.7\bin\gdal_retile.py", line 857, in main         dsCreatedTileIndex = tileImage(minfo,ti)       File "C:\PROGRA~2\FWTOOL~1.7\bin\gdal_retile.py", line 335, in tileImage         createTile(minfo, offsetX, offsetY, width, height,tilename,OGRDS)       File "C:\PROGRA~2\FWTOOL~1.7\bin\gdal_retile.py", line 467, in createTile         dec.uly+offsetY*dec.scaleY)       File "C:\PROGRA~2\FWTOOL~1.7\bin\gdal_retile.py", line 246, in getDataSet         s_band = sourceDS.GetRasterBand( bandNr )       File "C:\PROGRA~2\FWTOOL~1.7\pymod\gdal.py", line 654, in GetRasterBand         raise ValueError, 'Out of range band requested: %d' % i     ValueError: Out of range band requested: 2      I tried to retile separately images that execution stops on, but then I get no error message. Everything works well, so my guess that it's something wrong with the source image doesn't seem to be valid. Does anyone have any other ideas of what could be wrong or how I should troubleshoot the problem?


Fastest way to remove matched points

Hoping someone can see a faster way of doing this, because at the moment I'm averaging 9 hours per run for what I'm trying to do and well basically getting my behind chewed off by the boss :-) **The Scenario** We have a database running under PostgreSQL/postgis that's used to manage some mapping data for management of a transport fleet. Within this database there are several tables but for my question here, we are interested only in "built up areas" (aka limitedzones) and "Movement Speeds" (aka speeds) - I can't go too much into detail because of NDA's etc, but "limitedzones" is a table that contains a number of MULTIPOLYGONS each representing a geographic area, and has an associated speed limit that a vehicle is allowed to travel whilst inside the bounds of that polygon. "speeds" is a table of POINTS each point following the path of an allowed road and sampled with an allowed speed limit at 1 second intervals, representing a road that a vehicle is allowed to travel on. At present some of these "allowed roads" cross through "Limited areas" where both the "allowed road" and "Limited Area" speed limit are the same. I've been tasked with removing these extra points from the "Allowed Roads" table where they match the conditions set out above. **What I've tried so far** My first attempt was to simply select the points, so lets say for 40 Mph and put them into a memory table, then loop over that table comparing them to the 'limitedareas' table using "ST_Intersects" this didn't work as expected, and failed to remove all the points. It did however give me an idea of how many points I was looking at removing, and the count was easily in excess of 2 million. After a bit more fiddling around, and a few more attempts I finally arrived at the following:               DELETE FROM speeds s     USING limitedzones lz     WHERE ST_Within(s.geometry,lz.geometry)     AND s.speed = lz.speed      and in a similar fashion used the following to interrogate & count:               SELECT s.* FROM speeds s     JOIN limitedzones lz ON ST_Within(s.geometry,lz.geometry)     WHERE s.speed = lz.speed     --LIMIT 10      Both of these work exactly the way I expect, but the select (When not limited) takes about 3 hours to run, and the delete over 9 hours!! If at all possible I need to find a faster way of stripping these points out, as there is likely to be further work to add this into a regular update system, used ad-hoc on a frequent basis to keep the database in order. I have good indexes both B-Tree (gist) and regular on all the columns, and when individual tables are used on there own, everything is nice and fast, so I know it's the join on 2 tables that's causing me the slowdown. Any suggestions are welcome, and I'm even open to do this using a custom function / procedure if need be.


A quick-fix solution for geocoding?

I am working on a prototype for a web-based trip planner and need to have a geocoding function (for the start and destination adresses), up and running the "easy" way. My area that the prototype is covering is apprx.2*4 sqkm. Which option is the fastest and easiest (to learn) as I have very limited time? I have been looking into gisgraphy (downloading and geocoding in the database), nominatim/Mapquest etc. I prefer not be constrained to use a specific map provider, like Google or yahoo. However, if necessary it is possible to change geocoder later on... for now I need the "quick fix" Thanx!


How to symbolize more than one attribute simultaneously in QGIS?

Each symbol can show different properties such as temperature and pressure, etc. I want to change the symbol color according to one property and size with another one. How can I do it in QGIS? Thanks!


How do you scan/quantify vegetation using satellite images?

I am relatively new to GIS (about 4-5 months in) I did not receive any formal education about this so I am completely clueless about how to go about measuring/quantifying vegetation using satellite images. Can somebody please help? Is there any specific kind of image file I need to do this, if there is, how/where can I get them? Thanks! Enrique Roa Undergrad biology student


How to use multiple graphics layers in ArcEngine?

My organization is working on an ArcEngine application that displays a GPS cursor for the user's current position. We are using custom-made GPS code since ESRI's is buggy in ArcGIS 10. Whenever the GPS cursor's position changes, the entire graphics layer redraws rather than just the cursor. This results in a map whose labels are constantly flashing as the cursor moves. We are pretty sure that the problem is caused by having the GPS cursor displaying on the same graphics layer as the labels and any other graphics, but we haven't had any luck finding a way to use multiple graphics layers. **Is there a way to have multiple graphics layers in ArcEngine 10?**


OpenLayers wrapDateLine Problems with ArcGIS93Rest

I'm having issues viewing a ArcGIS93Rest layer on both sides of the dateline -- the data **doesn't cross** the dateline. The problem is that i only see data for one side of the dateline until i pan west and vice versa (see attached pics at bottom). I'm trying to rule out that it's an OpenLayers problem. The **wrapDateLine** option should work with the ArcGIS93Rest layer type because it's a subclass of Layer > HttpRequest > Grid same as OpenLayer WMS Layer. I have no wrapping problems with WMS layers created with UMN Mapserver. However, I can't set the wrapDateLine option as a hash extra option on the ArcGIS93 layer but only after intializing a layer like so. I'm not sure this is right. Is it?:               layer = new OpenLayers.Layer.ArcGIS93Rest('L' + cfg.id, cfg.url + '/export',         { layers: 'show:' + cfg.layer, transparent: 'true',             f: 'image', format: 'png', size: '256,256',             bboxSR: '900913', imageSR: '900913'         });     layer.setOpacity(cfgOpacity);     layer.wrapDateLine = true;      Questions:   1. Does anyone know if wrapDateLine works with ArcGISRest93? Have you seen it work?   2. Has anyone had this problem with ArcGIS Server data before? I'm wondering if it has to do with their service. As I said before, I have no problem with WMS layers from UMN Mapserver  ![west of dateline](http://i.stack.imgur.com/a5whx.png) ![east of dateline](http://i.stack.imgur.com/YV5fM.png)


How to handle AutoCAD Map to QGIS workflow?

I’m looking for some assistance and like many here, I am pretty much a GIS novice. I’m working on a project to map our local cemetery (about 40 acres), locating sections, lots, and graves, etc. I was able to obtain an aerial photo of the cemetery and, using a series of about 6 reasonably accurate GPS waypoints, create a georeferenced TIFF using QGIS and GDAL. This seems to have worked fairly well. Imported as a layer into QGIS, it overlays the Google imagery from the Open Layers plugin well enough for my current use. Since I have a little more drawing control, I’ve been trying to use AutoCAD map for creating the features (call this the “master.dwg”). I imported the TIFF file and drew a few features based on the TIFF image. Exported them to shape files, opened in QGIS and this looked pretty good as well (with the previous TIFF and Open Layers shown). My next step was to create a layout of one of the sections in AutoCAD map. I created a drawing using architectural units so I could enter values as 25’ 6”, etc. Once I completed the layout, I tried inserting the section into the master.dwg. For some reason, the scales seem to be massively out of whack. Basically, the inserted section is probably hundreds of times larger than the georeferenced TIFF. The large dimension of the section is only about 200 feet, yet when zoomed to extents, the TIFF file doesn’t even appear. It is necessary to zoom in quite a bit before it even appears as point in the drawing. Using the measure tool in AutoCAD, measuring the distance between two known points on the TIFF literally registers a Zero (I would have expected it to be about 200 ft). In QGIS, it measures about as expected. Any thoughts as to what is happening here? I’m at a loss right now. Matt =================================================== I've made a little progress. Today, I was able to get some shape files from our local GIS dept. of the area around the cemetery. I've tried MAPCONNECT and MAPIMPORT with the same results... I create a new drawing, assign the projection (MI83-SIF) and then attach / import the shape files. Units for the drawing are standard decimal. Using the MEASUREGEOM, I measure between two points and get 204 (which is what I would expect in feet). When I change the units to Architectural, then MEASUREGEOM returns just over 17'. So, it appears that, while "204" in the shape file should be feet, AutoCAD Map is interpreting it as inches.


directions (path find) given kml supplied bounds

Im trying to generate directions using google.maps.DirectionsService and I really need those directions bounded by the paths defined from within a pre- made kml file. I have reviewed the API docs for google maps and doesnt seem they have this functionality. This is for a custom app where areas are already defined inside a kml file that is layered over google maps. Any ideas? Thanks


WMS tile cache that works with time dimension

Is there a **good WMS tile cache** that understands the **WMS time dimension**? And what configuration is necessary? Open source preferred. I am serving some WMS layers with the time dimension, so the map requests from the client include dates. For example `http://…?…request=GETMAP&…&TIME=2010-12-10` returns the map for 10 Dec 2010. In my undestanding the tile cache will need to treat the `TIME` parameter as an identifier for the maps, just like the layer name. I'm aware of Geowebcache, Tilecache, and MapProxy but I can't find any reference to the time dimension in their documentation. The time dimension has existed since WMS 1.1.0 so I would hope it's been considered by someone? Our layers are served by ncWMS and MapServer, if that's relevant. We would like to use a tile cache to increase the speed for our most frequently requested content. Can anyone recommend a tile cache, and preferably point us to the relevant documentation too?


How to Geocode 300,000 addresses on the fly?

I have a database which has 300,000 address, which are to be shown on map. I know if I geo-code all the address it will be too expensive for me. So I was wondering if it is possible to geo-code the address on the fly/ realtime, when a user would select a address (a property address) it would search through the database and then geo-code the address and then map it with other attributes. It would be really nice if you could share a code, concept, or anything. By the way my backend is on mysql supported by Joomla. Thanks.


SDE Versioned Data Editing Problem

ArcMap is telling our users that data from a certain feature dataset needs to be registered as versioned to edit. ArcCatalog is showing us that the feature dataset in the SDE is not registered as versioned. My question : How does ArcGIS know that the data _should_ be versioned?


ArcObjects: Custom GP tools: GPCompositeDataType and DataType

I got an issue with the parameter value's datatype returned when using composite type parameters. In my case, I create a composite of DEWorkspaceType and DEFeatureClassType:               IGPCompositeDataType gpCompositeDataType = new GPCompositeDataTypeClass();     gpCompositeDataType.AddDataType(new DEWorkspaceTypeClass());     gpCompositeDataType.AddDataType(new DEFeatureClassTypeClass());      No matter what I select in my tool, e.g .gdb, .shp etc, I get "DEWorkspaceType" in the unpacked IGPValue.DataType. If I switch the sequence of the types, using DEFeatureClassTypeClass as the first datatype in the composite, I get DEFeatureClass as the parameter value type.. no matter what I select in the dialog. Am I left with parsing/guessing the datatype in the string from IGPValue.GetAsText() ? (version: 10.0 SP5)


arcpy Make Feature Layer using a where clause

From a table of lines and polylines I am I'm building a seperate table of points (nodes) with their X and Y coordinates stored as attributes. I want the nodes to be unique, so I'm checking to see if a node with the same coordinates has already been inserted by making a feature layer from the nodes table with a SQL clause to selecting any rows with the same coordinates. If the feature layer has zero rows then I insert the new node otherwise I do nothing and move to the next line/polyline. The problem I have is that the SQL clause in the Make Feature Layer seems to be selecting rows even where the coordinates do not match. It works fine for the first one but then after that it picks 1 row each time when most or all of the time it should pic zero rows. This is making my head hurt. Can any one help please? Here is the script..               import arcpy     arcpy.env.overwriteOutput = True          fc_Nodes = 'C:\Users\Rob&LP\Documents\GIS General\DRN Configuration\DRN Sample\LinkNodes_point.shp'          arcpy.DeleteRows_management(fc_Nodes)          fc_Links = 'C:\Users\Rob&LP\Documents\GIS General\DRN Configuration\DRN Sample\DRN_SevUp_14Recs_polyline.shp'          shapeName = arcpy.Describe(fc_Links).shapeFieldName          rows = arcpy.SearchCursor(fc_Links)     InCurs = arcpy.InsertCursor(fc_Nodes)          Ins = 0          for row in rows:           feat = row.getValue(shapeName)      ANode = feat.firstPoint               #check to see if node has already been inserted      arcpy.MakeFeatureLayer_management(fc_Nodes, "fl_P1", """"X" = ANode.X""")      if int(arcpy.GetCount_management("fl_P1").getOutput(0)) == 0:            Ins = Ins + 1            print "NumRows: " + str(arcpy.GetCount_management("fl_P1").getOutput(0))       print "Inserted: "+ str(Ins)       print "ANodeX: " + str(ANode.X)       print ""                 newNode = InCurs.newRow()       newNode.SHAPE = ANode       newNode.setValue("CHECK", 1)       newNode.setValue("X", ANode.X)       newNode.setValue("Y", ANode.Y)       newNode.setValue("ID", Ins)            InCurs.insertRow(newNode)       del ANode       arcpy.Delete_management("fl_P1")       else:            print "NumRows: " + str(arcpy.GetCount_management("fl_P1").getOutput(0))       print "NOT Inserted : " + str(ANode.X)       print ""       del ANode       arcpy.Delete_management("fl_P1")           arcpy.Delete_management("fl_P1")           del InCurs     del rows               print "--------------End---------"      `


Polygonize only pixel values of 1 from a raster

I have a raster mask that I have created in QGIS (2.2.0) that consist of 0, 1 and nodata values. I would like to end up with a polygon dataset that consist of only the value of one pixels from this same raster dataset. I have tried the polygonize tool, but it doesn't seem to be producing what I am looking for as it is taking both the 0 and 1 pixel values into the polygons. Any ideas one how I can accomplish this task would be appreciated. Thanks, Derek


Performing an Arcpy code using Python on all shapefiles in a folder or geodatabase?

I have a geodatabase with many polyline feature classes all with different names. I would like to write a code in python ArcGIS 10.1 to use a geoprossing tool on all classes in the gdb. How would this be done? Would the same logic work for a folder of shapefiles? For this question lets say I want to use `Create Random Points` 500 metres apart on all polylines in this gdb.               import arcpy     from arcpy import env     arcpy.env.workspace = "E:\Example.gdb"          arcpy.CreateRandomPoints_Management (      I asked a similar question on Stack Overflow, except using VBA to loop through a folder of excel workbooks. http://stackoverflow.com/questions/22022672/looping-a-code-through-a-folder- of-workbooks-with-vba


How can I show only some classes of a layer in UMN-mapserver?

I visualize a shape-layer in the mapserver. Now I want to make some requests for only some classes of this layer. I want for example one request for only the water-classes and later for only the wood-classes. Can I realise this with the mapserver?


scratch workspace in-line variable, up one level?

So I'm building a model to be published as a geoprocessing service, and I'm trying to use the 'Save to Layer File' tool to create a symobolized result of my raster output. I've set both current and scratch workspaces to separate file gdb's, and am using in-line variables (%workspace%, %scratchworkspace%) to direct my outputs. I want to use the %scratchworkspace% in-line variable for my 'Save to Layer File' output, but to my knowledge, layer files cannot be stored in a geodatabase... therefore, I need to know how I can set the 'Save to Layer File' output to the folder that the scratch gdb is stored in. Is it possible to specify my .lyr output to be saved "up one level" from the %scratchworkspace%??


How to do stereographic projection conversion (w.r.t. changing point of projection)?

Given a stereographic point P1 projected from point P, and given another point Q, I would like to know how to transform P1's coordinates to that of a stereographic point Q1 projected from point Q. The Q I'm using is pretty close to P, and the points on the sphere/earth to be projected are pretty near the antipode of P/Q, possibly clustered together, and having a high accuracy on each point in order to differentiate them is important. I'm wondering if a simple adding of an offset value to P1 will provide a good approximation? I need to preserve conformality and relative distances of projected points (since distance is not preserved in a stereographic projection, but I think that relative distances of points are preserved?). P's and Q's latitude and longitude will be known. ![2D illustration drawn in paint](http://i.stack.imgur.com/ubQ9C.png) Do let me know if this question is better posted at math.stackexchange.com


Postgres/postgis cursors

I have a road network with StreetName and min/max address ranges for each segment (ie- StreetName = 'Main St', Min = 100, Max = 199). I also have individual address points with similar information (ie- StreetName = 'Main St', HouseNumber = 115). I would like to assign the featureid (gid) of road segments to each address point that has matching StreetName and falls within the min/max range of the road segment. In ArcGIS I would probably write a searchcursor for the roads and that would look cycle through each road segment, get the gid, name, min, & max values, and then do a query statement to find all address points with HouseNumbers that fall within the given segment range (ie- select AddressPoints where StreetName = 'Main St' and min >=100 and max <=199), calculate the gid of the road segment into a field for the selected points in the Address Point layer, then move on to the next road segment and repeat. 1) How would I write a similar cursor in postgres/postgis or 2) Is there a more efficient way of doing this? Thanks.


Is there any free resource to get international raster data (other than USGS)?

I am looking for any raster database that allows free downloads of any kind of raster data pertaining to Damascus, Syria and the surrounding Turkey area. I'm trying to avoid USGS since it requires a lengthy sign-up. I am open to any site that would allow me to query for data based on that location and extract it for free without too much of a sign up process. I don't care if its DEM, topology, orthography, lidar, etc.


Problem of "null" values added for fields with no data when appending data or adding new features to file geodatabase

I am working with a bicycle and pedestrian facility file geodatabase where I am continually appending new records to feature class attribute tables. When I append new records, if there is no data for a certain field, ArcMap will automatically populate the records as "null" for that field. Is there any way to prevent the records from being populated as null? I am afraid that it will send the signal that there are errors in the database for database users. I am considering using "Find and Replace" or possibly a script to replace the "null" values but wanted to see if there was a way to prevent this since I don't want to constantly have to delete the null records every time I append data or add features.


GeoExt 1.1 FeatureStore featurefilter

I'm new to GeoExt and ExtJS, and I've started to set up an small project. Within this project I get my features dynamically generated by an MySQL- database as KML. I import the features as an OpenLayers.Layer.Vecor, store them in a GeoExt.data.FeatureStore and output them in a grid-panel. Now I'd like to filter my features according to certain parameters. So far I've found the following solution: How to use filter in GeoExt.FeatureStore? I've tried to implement this solution and Firebug does not ouput any errors in the script. Noneonetheless I do not get a filtered result in my grid-panel. I also tried to implement the filter within the protocol area, but I don't want to send a DB-request any time I set a filter parameter. Here's my codesnippet:               var events = new OpenLayers.Layer.Vector("events", {         preFeatureInsert: function(feature) {             feature.geometry.transform(               new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject()             );           }     });          var store = new GeoExt.data.FeatureStore({         layer: events,         fields: [             {name: "name", type: "string"},             {name: "description", type: "string"}         ],         proxy: new GeoExt.data.ProtocolProxy({             protocol: new OpenLayers.Protocol.HTTP({                 url: "./php/phpsql_events.php",                 format: new OpenLayers.Format.KML({                          extractStyles: true,                          extractAttributes: true                 })             })         }),         featureFilter: new OpenLayers.Filter.Comparison({                            type: OpenLayers.Filter.Comparison.EQUAL_TO,                      property: "name",                      value: "value"         }),         autoLoad: true     });


Why are services directory and arcgis manager not accessible?

I installed ArcGIS Server 10.0 for Java with its post install and its SP 5. When I try to access the services directory, I get a 404 error from the web server. I don't get that error when I try to get to the arcgis manager web interface. After displaying the arcgis manager web interface, when I try to login, I get an error: Login information is incorrect. The login used is arcgismanager that is member of agsadmin. The password used is the one specified in the post install. Do you have any hint?


Is it possible to change the projection used in OSM2PO? I used the following command to generate a `SQL` file which I loaded into my database:               java -jar osm2po-core-4.2.30-signed.jar prefix=at "filename.osm.pbf.pbf"      This produce the file `at_2po_4pgr.sql`. I then loaded this file into my postgres database:               psql -U postgres  -d routing -a -f at_2po_4pgr.sql      As this table is based on SRID:4326, I wondered if there was any way of changing it so that I can use a different projection system. The columns in the table are as follows:                 id integer NOT NULL,       osm_id bigint,       osm_name character varying,       osm_source_id bigint,       osm_target_id bigint,       clazz integer,       flags integer,       source integer,       target integer,       km double precision,       kmh integer,       cost double precision,       reverse_cost double precision,       x1 double precision,       y1 double precision,       x2 double precision,       y2 double precision,       geom_way geometry(LineString,4326),       CONSTRAINT pkey_at_2po_4pgr PRIMARY KEY (id)      I know that I can use `ST_Transform to convert`geom_way` into a different projection but I was unsure how to update any of the other values as I presumed some were based on the initial projection.


QGIS anno labels

Is there a way to copy and paste Anno Labels in QGIS? I have a layer that contains Parking Restriction information e.g. No Waiting or Limited Waiting 1 hour Mon - Fri 9am - 10pm etc. There are lots of variations of these restrictions. What I'm looking to do is label everything in a neat little rectangle (blue back ground, white text) but it seems that I can only do this manually by creating individual Anno Labels which I can't copy and paste! Med


Postgis 2.0 - public.geometry_columns view does not contain correct info for views

We have recently updated our database to Postgis 2.0. **In our system** we used to create views for tables containing a geometry column. After creating the view we also manually updated the geometry_columns table with necessary info. The views are used to create Geoserver layers and also to feed data for ogr2ogr. **The problem** in Postgis 2.0 is that geometry_columns table was replaced with a view that retrieves data from system catalogs. In our case, for the created views, we don't retrieve correct info in 'coord_dimension', 'srid' and 'type' columns.   _Info_ : after some research i've found out that geometry_columns view is using pg_constraint table to fill those columns, but for that view there are no constraints defined. Does anyone have any idea how can we fix this problem? How can we force gemetry_columns view to retrieve necessary info? Is there a workaround or other way in doing this?


What is the procedure to load LAS files in QGIS 2.0.1?

What is the present procedure to load LAS files in QGIS 2.0.1? After some searching I found a paper outlining the procedure; however, it was for 1.8.something, and related posts here are all rather dated.


How to calculate distances between addresses and cities and coast line?

I am totally new to GIS, and so I apologize if my question is not clear enough. My basic problem is the following:   * I have about 160,000 Danish addresses with longitude and latitude for each address.    * I would like to calculate the distance to the nearest big city (preferably by road and not straight line), and the distance to the coast.    * I have access to a topographical map of the roads and cities in Denmark in the scale of 1:1000000.    * The data structure is vector.    * The format is GDB (ArcGIS).    * The coordinate system is UTM32-ETRS89.    * I use windows 7 64 bit. I have the following question: Is there a free GIS to work with the GDB format (and preferably a reference to a tutorial or similar for absolute beginners)? Any comments are welcome, and I gladly provide more info if needed. Thanks a lot!!


How to add a WMS layer on Cesium Globe?

I downloaded Cesium Java Script library for creating 3D Globes from "http://cesium.agi.com/". Also I found that, they are using Bing satellite map as default. But I need to display layer from My Geoserver. I am very new to this Java Script library. Can you please tell me , How to add a WMS Layer from Geoserver on the Cesium globe? Is it possible? Thanks


Problem with ArcGIS annotation

i have a problem with annotation (ok i will explain) : in database i have :   1. parcelle : a geometry table that contain polygons   2. personne : a table contain parcel's owners    3. pll_personne : is the relationship (n,n) between the (1) & (2) ![enter image description here](http://i.stack.imgur.com/KfvbK.png) **what i want is : to display parcels with annotation Name of there owners see picture please** ![enter image description here](http://i.stack.imgur.com/OmSVW.png) **there is somebody that can help i'm very Jam** !! :(


Installing GCC correctly

I'm writing this as an utter Ubuntu newbie. I want to run the script found on Save points to a shapefile with GDAL. But I get the following error, yes this error is similar to the question posted here. Stick with me.               Python 2.7.2+ (default, Oct  4 2011, 20:06:09)      [GCC 4.6.1] on linux2     Type "help", "copyright", "credits" or "license" for more information.     >>> import point_store     Traceback (most recent call last):     File "<stdin>", line 1, in <module>     File "point_store.py", line 3, in <module>     import osgeo.ogr     File "/usr/lib/python2.7/dist-packages/osgeo/__init__.py", line 21, in           <module>     _gdal = swig_import_helper()     File "/usr/lib/python2.7/dist-packages/osgeo/__init__.py", line 17, in          swig_import_helper     _mod = imp.load_module('_gdal', fp, pathname, description)     ImportError: /usr/lib/python2.7/dist-packages/osgeo/_gdal.so: symbol       VSIFTruncateL, version GDAL_1.8 not defined in file libgdal.so.1 with link      time reference      After reading the question I reference above I tried this:               sudo pip install GDAL==1.8.1      And then I got the following error:               cc1plus: warning: command line option ‘-Wstrict-prototypes’ is valid for        Ada/C/ObjC but not for C++ [enabled by default]     extensions/gdal_wrap.cpp:149:20: fatal error: Python.h: No such file or      directory          compilation terminated.          error: command 'gcc' failed with exit status 1      Does anyone have an idea how I can get my gcc to happy fun times? There are lots of forum posts around about this error but I'm still looking for that "aha" moment.


Data from google maps or bing maps available for offline use? Whole cities, I mean

I was developing some graph-based search algorithms and it would be interesting to use them against some real data. I used openstreetmap against them (downloaded some city-wide OML files from some obscure site), and it all worked well until I started noticing that the data wasn't that trustworthy -- the data near my house had wrong highways directions, and there were roads that weren't properly connected. I was wondering whether I could download full-city google or bing (or actually, any found to be trustworthy source) maps data for this same batch processing? What I basically need is access to all the streets (and to know whether they are one or two-way streets) of a given town of my choice (Lisbon, for instance), as a combination of vertex + edges, or any other format that will let me gather that info, for offline use.


Can I resume (continue) digitizing a polygon feature from any vertex?

When I digitize in ArcMap with a Wacom tablet and a pen I often accidentally double click and finish my sketch. More often then not the accidentally finished sketch crosses itself. I would like to be able to continue digitizing the existing sketch from the last vertex or any vertex I choose. The OOTB "Continue Feature" tool does not work in this situation because it just adds a new feature part to the sketch and does not resolve the "crossing segment" artifact. I also tried using the ESRI's "Modify by Continuing Sketch" tool (developer sample) which does what I need it to do but it always resumes digitizing from the "Last" (red) segment. In theory this is correct unfortunately last vertex is often not the last digitized vertex - Editor seems to assign the last (red) vertex randomly. My question is: As a workaround, **I would like to be able to assign the "last vertex" to the actual last or any vertex in the sketch I choose. Can this be done?** _Figure: Example of an accidentally finished sketch (this is a polygon with a transparent fill) in which the last vertex was assigned to polygon boundary crosses itself not where I finished digitizing (which is indicated by an arrow). This also happens with simple polygons - last vertex becomes another random vertex, not the actual last vertex I created in the sketch._ ![enter image description here](http://i.stack.imgur.com/0dC7h.jpg)


Moving Marker in a Google Map on Android

I am trying to create something similar to http://dev.syskall.com/map/ in my android phone application. Basically I have a List of GeoPoints and I need to make the car/marker move on the map according to the values of the geopoints. Putting the mapView code in a loop does not help at all! The map is invalidated till it reaches the last location in the list and then shows up with the latest location only so I cant see any movement of the marker. If anyone has created this on android - please do get back to me! Appreciate your help very much!


Add WMS group layer (only) to ArcMap

I am trying to add a WMS group layer to ArcMap, without having the sublayers being listed in the TOC.   I will give a little background here: We are sourcing images from Nearmap.com, and their WMS data is served as follows: NearMap   |   Australia   * Adelaide    * \-----2009-10-19    * \-----2009-10-29    * Perth    * \-----2010-10-12  As per Nearmap, the group layer link (Adelaide or Perth shown above) always points to the most recent image avaible, while the sub-layers present date specific images (one of them would of course be the latest one).   If this is true, it releives us from updating and re-publishing our maps every time a new image is available from Nearmap (which happens at least once a month).   When I try to add the group layer in Arcmap, it adds all the sub-layers along with it, and I have to choose a specific date image to publish in our web maps. Does anyone here has experienced the same behaviour in other GIS software, open-source or otherwise?   I was informed by the Nearmap guys that "other" open source software have the "ability" to add group layer directly...any takers for that? Many Thanks   Ujjwal


Use VB to calculate an area of a shape

I'm using Arcmap 10.1, running python 2.7 on a Windows 8.1 machine. I have a script that uses the arcpy.CalculateField_management function to calculate the area of a shape, however, for some reason when i use the python expression to calculate the field, python crashes. The only way i can use the arcpy.CalculateField_management function is by using a VB expression. Does anyone one how i can translate this line to VB `arcpy.CalculateField_management("TAZ_un_Buffer_Dissolve1.shp", "totalArea","!shape.area!", "PYTHON", "")` Have in mind that the line below returns an error. '`arcpy.CalculateField_management("TAZ_un_Buffer_Dissolve1.shp", "totalArea","[shape.area]", "VB", "")`' How else can i accomplish this without using the arcpy.CalculateField_management function?


How to populate fields for line segment with value from polygon to left and right of it?

I have a countywide roads dataset that is used by an emergency dispatch system. There are two fields, "Comm_RT" and "Comm_LT", that need to be populated by the name of the community polygon they reside in. I am having trouble finding a way to populate these two fields when the road borders two different community polygons and thus requires that the "Comm_RT" and "Comm_LT" fields be different. I have started by making selections in ArcMap 10.2 of all the roads that intersect two different polygons. Here is an image of a road segment that would have two different values for "Comm_RT" and "Comm_LT": ![Differing Comm Codes](http://slu.s3.amazonaws.com/comcode.PNG) Is there a way to perform a select by location query or a tool in the Network Analyst extension that could be used to solve this?


Openstreetmap of a specific city of USA

I want to get the osm of a specific city of USA. I can only get the whole state now. But how can I get the osm data for a specific city.


Extracting nodes and edges from osm data

I am trying to extract the edges and nodes of a certain city from osm files and save them to a database. How can I accomplish this. I am not sure how much osm2po and osm2pgrouting tools are acceptable. I mean what is the standard procedure to get the graph. I mean when I use osm2po and tools like that they generate their own data structures. I mean they have this column called geometry of type multiline string and stuff. I am not sure if I should use this. Actually, I want to implement shortest path algorithm of my own. I am not sure if osm is the only source to get such graph. Also, I want to plot the generated graph in matlab. How can I accomplish this?


TMS format specification

I want to convert large map file to tms format. how directory structure, zoom levels , tile size, xml file creation have to be dealt with. Please point to the relevant docs. And is there any gdal c++ api to create tiles(not gdal2tiles.py). Thanks


GeoTiff format (tile and overview confusion)

I have two questions :   1. How does tile creation in geotiff improve performance and how will reading of tiled gtiff images be done?   2. Adding an overview to a geotiff image means we are resampling the original image to a small size and storing it with the original image so that when we are trying to read a large image, the overview will be read first and then the original image (so that reading latency is hidden). Am I right on this? What is the point of adding more than one overview? Can somebody explain briefly about the geotiff format? Thanks


Does ArcGIS Silverlight API 2.4 support Microsoft Silverlight 5.0?

Please let me know whether ArcGIS Silverlight API 2.4 version officially supports or CERTIFIES Microsoft Silverlight 5.0.?


After re-open of old project file, several layer selections no longer available

I'm new to QGIS. I created a QGIS project, loaded 1 to 3 shapefiles via Add Vector Layer (from 2008 TIGER/Line Shapefiles http://www2.census.gov/cgi-bin/shapefiles/state- files). The maps look fine and I can interact with them. Under Layer many of the "middle" options are available after loading the shapefile, e.g., Open Attribute Table, Toggle Editing, Save as ..., Properties. I then save the project. After closing QGIS and re-opening the project file, all of the layer options mentioned above (and a few others) are now greyed out and unavailable. Are these options only available right after loading a shapefile? What good is a project file if one cannot interact with data saved with the project? And why are the Edit options always greyed out and not available? What am I missing? Thanks for helping me understand QGIS.


Applying unique processing extent to each iteration of tool in ArcMAP Model Builder

I am trying to create an automated ArcMAP (10.1) tool that will read in a file (ASCII most likely) of maximum and minimum lat/long values along with some additional variables. The tool will iterate over each row of this file, and run the processes on the unique variables each time. A preliminary tool that only runs on one set of variables works perfectly. ![Working Model](http://i.imgur.com/zFI8M0Z.png) However, the 'Extent' Data Type can't accept variables into it (in any way I have found out at least). So I have no idea how to define a unique processing extent for each run of the tool. Below is an example version of what the model might look like (I have tried several different ways, this is just one). ![Batch Model](http://i.imgur.com/UhWKob1.png) If it matters, I have never tried working with python scripts in ArcMAP, so I'm hesitant to go down that route unless it's fairly straight-forward. Unless there is no other option. The base DEM I will be using will cover the entire globe, so it really is necessary that I only process a small part of it each time I run the tool. Any help is greatly appreciated!


Possible to Clip Imagery Services Layers from ArcGIS Online?

The new Landscape layers such as the various Soil type imagery services look like a great resource. Does anyone know if it possible to clip out small areas of the Landscape imagery service since I only need the soil data on a small area. Clipping it will look much cleaner if I can clip to the feature class of my area of interest. When trying to run the Clip tool in ArcMap 10.2, I can't get the imagery service to show as an option to clip. Any help on this would be appreciated.


Stack - Mapnik, TileMill, OpenLayers, GeoExt

Would this work as a stack for serving webmaps with the basic GIS viewer functionality (pan, zoom, layer switching, info popups)?   * Apache Web Server   * Mapnik to serve the geo-data (TileMill to Style the Mapnik data)   * Open Layers to display in browser   * GeoExt or JQuery for added functionality Do I have that right? Am I missing anything? I've read through a bunch of the threads and other webmap info, still not sure exactly where/how Mapnik fits in. It can replace GeoServer/MapServer in a stack? And TileMill, not sure how the styled map gets added to the server. Thoughts and advice very much appreciated!


C++ coordinate class that automates projection?

Can anyone point me to a C++ class for representing coordinates that automatically handles projection transforms? What I'm thinking of is something that would allow me to create vars that incorporate a projection and have appropriate geographic operators such that comparisons make sense even when the coordinates are in different projections. Using template notation (though I can't figure out how to template on a projection) > Coord < SRID:1234> point1 (x1, y1); Coord < SRID:5678> point2 (x2, y2);               if (point1.east(point2)) ... // true if point1 is east of point2      or               point2 = point1 // This operation assigns the coordinates of point1 to point2,      using the coordinate system of point2 (in this case SRID::5678).


Large mosaic process failing - what are the options? software independent only focussed on result

I need to mosaic about 550Gb of tif imagery together and the software I have tried keeps failing. The area has been split into zones so that the smallest has approx 200 tiles. I have used the latest versions of ERDAS (Imagine and Mapper), ArcINFO and Global Mapper on a 3.30 gigahertz Intel Xeon E31245, DELL, 16GB RAM, 64-bit Win 7 Professional. Mullti-core (4 total), Hyper-threaded (8 total) machine. My C has 700GB free and D has 1.5TB. I am looking into using Grass (never have before) but i.image.mosaic only seems to only handle 4 files...some of mine have 600 tiles. Any other options or opensource software to try? Sorry should add that we can't use a mosaic dataset (or equivalent in other software) as we need to create zones with defined no-data areas as ecw's so that they can be opened in any GIS software and combined with lower resolution/older data when new data doesn't exist seamlessly. ![enter image description here](http://i.stack.imgur.com/2pPT4.jpg) An example of how some mosaiced files look in different sofware. Global Mapper/ERDAS are fine but it's not correct in arcgis. \--- OLDER INFO--- ![enter image description here](http://i.stack.imgur.com/plDGq.png) Sorry for the rough drawing. So having the colored areas as 5 zones will minimize the no data areas in the larger AOI. In arcgis the code is as follows (this is run as a model and not in python as I can't get it to take the tifList input).               arcpy.MosaicToNewRaster_management(tifList+";" +mask,RootOutput,"Tile1.tif","PROJCS['GDA_1994_MGA_Zone_55',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',147.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]","16_BIT_UNSIGNED","0.5","3","MAXIMUM","#")          # Replace a layer/table view name with a path to a dataset (which can be a layer file) or create the layer/table view within the script     # The following inputs are layers or table views: "test2"          arcpy.CopyRaster_management(OutputFile,RootOutput+"Tile1b.tif","#","256","256","NONE","NONE","16_BIT_UNSIGNED")      where tifList should be read in from a csv file but this didn't work in python so I am running the above in a model instead... I have 1.5TB+ free space on my drive but the process crashes with a 9999 error. Would even 100 tiles process? -i.e should we look at breaking the zones up further?


OpenLayers Scaleline always wrong in my virtual star map

I got a problem for the scale line when I using openlayer to show my "star map". Basically, I give a set of 2D coodinates (X,Y) for each star, and generate the map using map server. Now the map can be correctly shown, but when I add scaleline control, I found the scaleline always starts at 500000000KM which is not correct. And the correct scale should be 5400KM, which if I directly use mapserver CGI, it show correctly. Hope anyone can help me out! Here is mapfile : `MAP NAME "starmap1" PROJECTION "init=epsg:4326" END UNITS meters IMAGETYPE PNG EXTENT -31230139.51 -8267352.329 36978851.96 9060799.858 SIZE 800 600` Here is my javascript (using Geoext ):                   panel = new GeoExt.MapPanel({         renderTo: "mappanel",         map: { maxExtent: new OpenLayers.Bounds(-31230139.510000,-25170977.043091,             36978851.960000,25964424.572091),               maxResolution: 85367.949274,             numZoomLevels:40} ,             width: 800,             height: 600,             floating: true,             x: 70,             y: 350         });       panel.map.zoomToMaxExtent();        panel.map.addControl(new OpenLayers.Control.ScaleLine());      Thanks zhizhuo