View source for Exercise: Measuring point coordinates using WMS Google maps

Coordinates are often needed for ground control points (GCP) with the objective to georeference remote sensing data or scanned maps without graticule. Google maps can be used as one source for registering point coordinates.

== Import highly accurate DGNSS (Differential Global Navigation Satellite System) points stored in a ASCII CSV file ==
# Click [[file:QGIS_2.0_Add_delimited.png|25px]] {{button|text=Add delimeted Text File }}. Select the input file ''GBData\GNSS\Festpunkte_Forstbotanischer_Garten.csv''. Define X field: ''East'' and Y field ''North''. {{button|text=OK}}. (Fig. A) [[Image:importascii.png|thumb|right|250px|'''Figure A:''' Import of point coordinates from a comma separated ASCII file]]
# Specify the coordinate reference system (CRS)  ETRS / UTM zone 32N. Search for the EPSG code: ''25832''. Three precise geodetic survey points are loaded into the viewer. The current active project CRS is displayed on the lower left frame of the viewer as EPSG:25832.

== Use Google maps as background layer ==
# Activate {{mitem|text=Web --> Openlayers plugin --> Google Maps --> Google Satellite}}. A high spatial resolution aerial othophoto originally aqcuired in 2009 by GeoBasis-DE/BKG (Bundesamt für Kartographie und Geodäsie, Frankfurt and provided as web map service (WMS) by Google is displayed in the viewer. Note, the active CRS is automatically changed to EPSG:3857 (WGS84 / Pseudo Mercator).
# Label the geodetic points. {{mitem|text=Layer --> Properties --> Layer --> Labels}}. Switch to {{button|text=Show labels for this layer}} in the drop list and Label with the field ''ID''. Change text size to ''12'' and text color to white.
# Zoom close in to a geodetic point. What do you notice? The Google layer is misaligned when zooming in to about 1:1000 scale! Find the maximum zoom level where Google layers are still correctly displayed (Fig. B). [[Image:vp_gm.png|250px|thumb|right| '''Figure B:''' Geodetic survey points overlaid on Goggle Satellite layer using Openlayers plugin.]]

== Display GNSS points of a long term measurement before and after differential postprocessing ==
# Add uncorrected points of a long term measurement using GNSS receiver Trimble GeoExplorer 6000. {{mitem|text=Layer --> Add layer --> Add vector layer}} open the shapefile ''\GBData\GNSS\PosPkt.shp''. Zoom in to point ''VP1'' and observe the scattered point cloud. Measure the extent of the GNSS point using the ''Measure line Tool''.
# Now open the GNSS points as a result of a differential postprocessing: {{mitem|text=Layer --> Add layer --> Add vector layer}} and open ''\GBData\GNSS\PosPkt_DGPS.shp''. Measure the extent of the postprocessed GNSS point using the ''Measure line Tool''. Which horizontal accuracy can be expected? 
== Register point coordinates based on Google Maps ==
# Select {{mitem|text=Layer --> New --> New Shapefile Layer}}. 
# As layer type, select {{button|text=Point}}. Click the {{button|text=Specify CRS}} button and select ETRS89 / UTM zone 32N (EPSG:25832). 
# Add a new attibute: type ({{typed|text=vpid}} into the {{button|text=Name}} field of the {{button|text=New attribute}} section.
# Select {{button|text=Text data}} as data type.
# Confirm with {{button|text=Add to attributes list}}. 
# Confirm with {{button|text=OK}} and enter path and  file name (e.g. {{typed|text=VP_gm}}) in the following menu.
# Select the new shapefile in the Layers Panel. Start the edit mode by clicking the {{button|text=Toggle editing}} button [[file:QGIS_2.0_Edit.png|25px]].
# Zoom in to VP1 (Fixed point marked with "Meisselzeichen") with maximum zoom level where Google layers are still aligned. Click {{button|text=Add feature}} [[file:QGIS_2.0_Addpoint.png|25px]] to start digitizing and mark a point with left click and entering the attributes id = {{typed|text=1}} and vpid = {{typed|text=VP1}} in the appearing window.
# Switch off the edit mode of the point layer by clicking the {{button|text=Toggle editing}} button [[file:QGIS_2.0_Edit.png|25px]].
# {{mitem|text=Vector --> Geometry Tools --> Export/Add geometry columns}}. Select the vector layer {{typed|text=VP_gm}} and click {{button|text=OK}} and {{button|text=Yes}}. Close the dialog with {{button|text=OK}} and {{button|text=Close}}.
# Select the {{typed|text=VP_gm}} in the layers panel. Right click {{button|text=Open Attribute table}}. Report the coordinates ''XCOORD'' and ''YCOORD'' and calculate the Euclidean distance [meters] to ''East'' and ''North'' of the precise survey point VP1 in the attribute table {{typed|text=Festpunkte_Forstbotanischer_Garten}}.


[[category:Forestry Application]]