Generating high-resolution digital elevation models for wetland research using Google EarthTM imagery – an example from South Africa E Hoffmann and F Winde*
North-West University, Potchefstroom Campus, School of Environmental Science and Development, Private Bag X6001, Potchefstroom 2520, South Africa
Abstract Digital elevation models (DEM) generated in geographical information systems (GIS) have proven to be useful tools in hydrological research, aiding, amongst others, the delineation of catchment areas, identification of drainage patterns and flow pathways as well as for runoff determinations. They are of particular value in areas of comparatively flat topography, where such tasks are often difficult to perform. However, owing to the fact that elevation differences in wetlands typically fall below or just into the range of contour intervals of standard topographic maps being generally 20 m, and 5m for some areas, the latter fail to give enough detail. This means that sufficiently detailed relief information is often difficult to obtain for wetland research. Sitespecific, high-resolution relief surveys are too expensive, relative to many research budgets, to constitute a viable alternative. Based on an approximately 12 km² study area surrounding a karst-related peatland, this paper presents a method to retrieve the required high-resolution elevation data at 1 m intervals, at low cost, from satellite imagery in Google EarthTM. The paper describes procedures used to capture and process the data using GIS ArcDesktop™ to produce a high-resolution contour map and DEM. For quality assurance purposes the generated map is visually compared to 5 m and 20 m contour intervals of standard topographic maps (1:50000) issued by the Chief Directorate Surveys and Mapping (CDSM). After correcting an off-set of 5 m it was found that the deviation of the generated contour map based on Google data from CDSM contours was in the same order as the deviation between the 2 CDSM contour sets. Finally, all 3 contour maps were compared to a contour map of 0.5 m-interval resolution specifically generated for the study area using aerial photography from an airborne survey. This too confirmed the overall good reliability of the generated, Google EarthTM-based, contour map. Although Google EarthTM’s contour models are based on data of the Shuttle Radar Topography Mission (SRTM), the direct use of freely available SRTM data for localised, high-resolution DEM yielded unsatisfactory results. This may be due to (unspecified) procedures, or unknown data sources employed by Google EarthTM that enhance the quality of relief data. With updating intervals of 2 to 4 years, satellite imagery in Google EarthTM offers the additional advantage of containing much more recent information on relevant hydrological features than the outdated topographic maps available for the study area. It is concluded that the presented method allows the generation of high-resolution DEMs especially for areas of flat topography where adequate relief information is either not available or too costly to generate. These DEMs are useful for further wetland research.
Keywords: digital elevation model, DEM, Google Earth, GIS, wetland, hydrology, relief, contour intervals, Gerhard Minnebron peatland, SRTM, topographic maps
Introduction
27°0'0"E
27°15'0"E
27°30'0"E
27°45'0"E
28°0'0"E
28°15'0"E
Wonderfonteinspruit catchment
GMB
Randfontein
Kerkskraal Dam
0
r
Wonderfonteinspruit
10
Westonaria
Khutsong
20 Km
1:500,000
Carletonville Welverdiend Blyvooruitzicht
Deelkraal 26°30'0"S
Gerhard Minnebron
26°15'0"S
Mohlakeng
Mooi Rive
Legend Du
To it S
Gerhard Minnebron catchment
Canals near Gerhard Minnebron Non-perennial stream
pr uit
perennial stream
1
Boskop Dam
R 50
Owing to potential threats, associated with continued and poorly-regulated peat mining activities, to the ecological integrity of a peatland in the semi-arid karst region of the North-West Province (South Africa), the Department of Water Affairs and Forestry (DWAF) (now the Department of Water Affairs (DWA)), in conjunction with the Department of Agriculture (DoA) initiated a research project to assess impacts of peat extraction on the hydrological functioning of the Gerhard Minnebron (GMB) peatland (see Fig. 1).
Dams near Gerhard Minnebron Potchefstroom
Wetlands
* To whom all correspondence should be addressed. +2718 299 1582; fax: +2718 299 1582; e-mail:
[email protected] Received 20 April 2009; accepted in revised form 12 November 2009.
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
Gerhard Minnebron wetland 0
2,500
26°45'0"S
Settlements
Dolomite
5,000 m
Figure 1 Location of study area and upstream hydrography
53
E 28°
27° E
96 km
it inspru erfonte Wond
91 k
m
Mo oi R i
ver
26°S
GMB wetland
Height difference: 220 m
N
Figure 2 Example of a (triangular irregular network-based) DEM for an area of approximately 8 700 km² surrounding the studied wetland at Gerhard Minnebron (GMB), generated from 20 m contour intervals extracted from 16 topographical 1:50000 scale maps, comprising all of the sheets of quadrants 2627 A through D, as issued by the CDSM
27°S
In the course of the study a GIS (geographic information system) was implemented to spatially represent data acquired for this wetland. This paper focuses on how high-resolution elevation data acquisition and integration into a GIS can be done cost-effectively by using virtual globe systems such as Google Earth™ to model areas with little change in relief. Hydrological modelling and GIS development were historically used in parallel but as separate avenues of specialisation. This continued until the late 1980s, when efforts to improve the analytical abilities of GIS began to blossom and researchers began to draw these 2 fields closer together (Raper and Kelk, 1991; Sui and Maggio, 1999). This fusion has resulted in the integration of hydrological techniques by GIS program suites (Sui and Maggio, 1999). ESRI’s ArcDesktop™ suite is one such a program suite and was selected to set up a GIS for this study. Its 3D Analyst extension was of particular use and others, for example, ArcHydro, show considerable value regarding the visualisation and modelling of hydrological processes and features. In spite of the challenges faced regarding the use and integration of GIS in hydrological studies, there have been numerous examples of such endeavours, as is noted by Sui and Maggio (1999) and Raper and Kelk (1991). Wetland modelling-related integration examples can also be seen in studies regarding runoff processes (Drogue et al., 2002; Ko and Cheng, 2004; Soulsby et al., 2006), flood predictions (Al-Sabhan et al., 2003), water quality concerns (Janauer, 1996; Bobba et al., 2000; Janssen et al., 2005; Jewitt et al., 2004), and water balance questions (Alemaw and Chaoka, 2003; Fraser et al., 2001; Lacroix et al., 2002; Portoghese et al., 2005), to name but a few cases. A crucial component in modelling of hydrological processes such as surface runoff, catchment delineation and the determination of the nature of surface-ground water interactions is a DEM (digital elevation model) that can accurately capture the terrain in sufficient detail. DEMs also allow surface flow pathways to be identified without resorting to complicated automated hydrological modelling techniques. Relief reflected by standard topographic contours at a 20 m interval is acceptable for DEMs at regional scale, even for terrains with comparatively flat topography, such as the Highveld (see Fig. 2). However, such an interval is commonly insufficient to capture the required degree of detail in smaller areas, especially
54
wetlands with their characteristically small changes in relief. The hydrology of the larger study area is controlled to a large extent by the presence of highly karstified dolomite which hosts some of SA’s largest aquifers. The Gerhard Minnebron wetland itself owes its existence to the continuous outflow of dolomitic groundwater from the Gerhard Minnebron karst spring, yielding some 60 to 80 Mℓ/d. While falling outside the surface catchment of the Wonderfonteinspruit (WFS), which is the single largest tributary of the Mooi River, the Gerhard Minnebron Spring is possibly hydraulically linked to the WFS system via a network of underground karst channels present in what is known as the Boskop-Turffontein compartment (Winde, 2006). The location of the GMB wetland in relation to outcropping dolomite is shown in Fig. 1. At the Gerhard Minnebron wetland the maximum vertical difference between the highest and lowest elevation points does not exceed 15 m. This means that in standard topographic maps, such as the 1:50 000 series (which is the highest resolution map series covering the whole of South Africa) displaying 20 m contour intervals, this vertical difference in topography will be lost. For areas covered by (aerial) orthophotographs at 1:10 000 scale the Chief Directorate: Surveys and Mapping (CDSM) is able to provide derived contour information at 5 m intervals in digital format. While this greatly improves the vertical resolution of DEMs at sub-regional scale it is still insufficient for adequately capturing hydrologically-relevant topographic detail in wetlands, where relief changes well below 5 m are often crucial for determining flow patterns. In South Africa professional consulting firms specialising in aerial surveys are often the only alternative source for generating sufficiently detailed relief data. Depending on the location of the site (distance from airport) this option is commonly too costly for many research budgets. In the particular case of the 10.8 km² study area at GMB, costs for the aerial survey at 0.5 m contour resolution amounted to ZAR100 000 (Minnaar, 2009) equating to ZAR9 259/km 2. For ’metropolitan and growth areas the CDSM provides higher resolution DEMs (not contour maps) based on 50 m x 50 m (0.25 ha) grid cells at the cost of shipping only. For selected areas in more remote regions DEMs of somewhat lower resolution (200 m x 200 m – 4 ha and 400 m x 400 m grid size – 16 ha) are available (CDSM, 2007).
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
N
Gerhard Minnebron Karst spring
Gerhard Minnebron wetland
Elevation of cursor position [m above mean sea level]: 1401m
Coordinates of cursor position [degree minute second]
Cursor position
Date of image [MM DD YYYY]
Figure 3 A screenshot of Google Earth™ showing the study area in 3-D fashion with status bar at the bottom indicating elevation and coordinate reading for the cursor position on the image
However, for the vast majority of areas such information is not available, as was the case with the study area. This paper presents a cost-efficient alternative for creating the required level of vertical resolution by using satellite imagery provided, at a relatively low cost via the internet, by Google’s Google Earth™ application. The core of the method is based on elevation readings given by Google Earth™ for any position of a cursor which is moved across a chosen satellite image. Depending on user preference this elevation is provided in feet or meters above mean sea level (m a. m. s. l.) at intervals of 1 ft (approximately 30.5 cm) or 1 m, respectively, and displayed in a control bar at the bottom of the image. Together with the elevation, the exact coordinates for the locations indicated by the cursor position are provided (Fig. 3). In the specific case of this study it means that Google Earth Google EarthTM’s elevation data are at a resolution 5 to 20 times higher than available datasets provided by the CDSM. In Google Earth Google EarthTM this dataset supports an additional function which allows users to view the relief in their area of interest in 3-D fashion with a maximum vertical exaggeration factor of 3. For some areas in the GMB wetland this exaggeration was insufficient to visualise height differences important for determining pathways of surface flow during flood events, for example. However, once these data have been imported into a GIS the exaggeration factor can be chosen at will. In order to use this detailed relief information provided by Google EarthTM in a DEM a method had to be designed to import Google elevation data into a usable GIS format for use in ArcDesktop™. This method is presented in this paper, accompanied by a critical evaluation of the quality and reliabi lity of the resulting GIS-based DEM.
Origin and accuracy of Google EarthTM elevation data Google Earth™ provides high-resolution elevation data using the virtual globe system, which started in June 2005 and used Shuttle Radar Topography Mission (SRTM) data for its elevation baseline.
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
The SRTM was an initiative, by the National Aeronautics and Space Administration of the USA (NASA), the USA’s National Geospatial-Intelligence Agency, the German Space Agency and the Italian Space Agency, to produce a standardised, complete and high-resolution DEM of planet Earth using Interferometric Synthetic Aperture Radar (InSAR) technology (Farr et al., 2007). Taking 10 d or 149 orbits the mission covered an area between 60° N latitude and 56° S latitude leaving the largely uninhabited polar regions out. Each orbit produced a strip of images in 225 km wide swaths (Farr et al., 2007). Almost the entire targeted area (99.96%) was mapped at least once, with only a few patches over North America having been missed. Data errors of this mission are similar to those of the USA’s National Elevation Dataset (Farr et al., 2007). With InSAR cells generated at 1 arc-second by 1-arc-second (~ 30 m x 30 m) the vertical height error for absolute error (i.e. measured against ground truth) is less than 16 m and the relative (i.e. measured against generated data) linear vertical height error is less than 10 m (Farr et al., 2007). Regarding the circular and horizontal geo-location, SRTM data have an absolute error of less than 20 m and a relative error of less than 15 m. It should also be noted that this method is not able to penetrate the vegetation canopy and thus, in densely-vegetated areas, overestimates ground elevation by the height of trees (Farr et al., 2007). Since the semi-arid Highveld, into which the study area falls, is only sparsely vegetated with bushes and scrubs this error is of marginal, if any, importance. Since SRTM data underlie the Google EarthTM DEM, attempts were made to use these data directly for generating a high-resolution relief model of the wetland. However, the contours derived directly from SRTM data (as downloaded free of charge from the website of the Consortium for Spatial Information – an initiative from the Consultative Group for International Agriculture: http://srtm.csi.cgiar.org/), by sub jugating them to high resolution interpolation, yielded unsatisfactory results with many contours showing an unnatural, blocky appearance and frequently being open-ended (Fig. 4, next page). A possible reason for the fact that contours based on the original SRTM data are less satisfactory than those generated via Google EarthTM may be related to the continuous refinement of elevation data in Google EarthTM through successive addition of high resolution data from other sources as they become available for different areas (Taylor, 2009). This is implemented by not only utilising satellites but also other (unspecified) technologies (Aubin, 2009). Another reason may relate to the fact that SRTM-data had to be cut out from a larger set covering a 5° latitude by 5° longitude grid cell since the computational power required to generate 1 m contours for the entire grid cell exceeded the capacities of available computers. However, extracting the data relevant for the study area results in some deviation from the main dataset owing to interpolation-related errors associated with the newly defined margins of the cut-out area. This resulting misfit between original data from the entire grid cell and cut-out data for the study area is illustrated in Fig. 5 using 20 m intervals as an example (instead of 1 m intervals to reduce the computational requirements). In view of the above-described difficulties the extraction of improved elevation data as used in the final Google EarthTM product was attempted. In addition to improved accuracy compared to SRTM data this had the additional advantage of the relief information being geographically linked to continuously-updated land use data. Since the inception of the
55
27°7'30"E
tS
pr u
it
27°8'0"E
´
Non-perennial stream Perennial stream Wetlands
0
1.000 Meters 27°8'5"E
56
62,5
125
250 Meters 27°8'15"E
26°29'35"S
1400 m
26°29'25"S
26°29'55"S
26°29'15"S
1400 m
26°28'55"S
26°29'5"S
Deviations between original and cut-out SRTM contours
500
500
(For the elevation ≤ 1405 m a.m.s.l)
1.000 Meters
26°29'30"S
Legend 20 m SRTM contours for an entire tile Google Earth derived contours
250
250
5 m SRTM contours for the study area 1 m SRTM contours for the study area
26°30'0"S 27°8'30"E
20 m SRTM contours for the larger study area
0
Non-perennial stream Perennial stream Wetlands
26°29'0"S
Figure 5 (bottom) Deviation between the 1400 m contour interval generated based on original SRTM data as contained in 5° grid cells and SRTM data extracted from the dataset for the study area (indicated by double pointed arrows). Both contour sets not only display a blocky and unnatural appearance but also differ from the corresponding contour generated using elevation data extracted from Google Earth
27°7'30"E
27°9'30"E
26°28'30"S
Mo oi R
Figure 4 (right) 1 m interval contours of GMB (the wetland) as derived from SRTM data displaying an overall unnatural and blocky appearance
27°9'0"E
Legend
0
r
To i
27°8'30"E
i ve
Du
27°8'0"E
peat project in late 2005 the imagery covering the wetland was updated 3 times, which is in line with the average update interval indicated by Google EarthTM (Fig. 6, next page). Compared to available topographic information on the area, which dates back to 1975 (1:10 000 black and white orthophotographs) and 1995 respectively (1:50 000 map series, sheet Rysmierbult 2627 AC), Google EarthTM imagery is significantly more recent. This is of particular importance with regard to possible quantification of impacts associated with peat mining, which results in a visible reduction of wetland areas and increasing extension of open water surface. With most of the peat mining activities commencing on a large scale after 2003, these impacts on the area’s hydrography are not captured in any of the standard topographic sources. However, owing to a rather low resolution of the first set of Google imagery for the study area, many of the required details were still obscured. Resolution improved significantly with the 1st update covering the whole area of interest by images dated August - September 2004 (Fig. 6). At this resolution, details such as the size and shape of the several ponds created by peat mining within the wetland, the course of the drainage line, and the flow status of nearby streams and dams were discernable. The situation improved further with the 2nd update of imagery, after which the majority of the study area was now covered by satellite imagery as recent as June 2007 (Fig. 6, next page).
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
August 2004
WF S
River Mooi
June 2007
WF S
River Mooi
WF S
River Mooi
Pre-2004
open water
open water
open water
Gerhard Minnebron
Gerhard Minnebron
Gerhard Minnebron
Figure 6 Sequence of Google EarthTM images of the wider study area, stretching from approximately 7 km upstream of the confluence of the Mooi River and the Wonderfonteinspruit (WFS), as updated between 2006 and 2008, displaying dates when underlying satellite imagery was captured. (Owing to increasing resolution in the more recent images a higher eye height had to be chosen in order to keep the size of the displayed area constant)
Generating high-resolution DEM using Google Earth™ elevation data
SRTM
elevation data Elevation data from other sources
Google intren processing
The generation of a high-resolution DEM consisted of 3 general steps: • In the 1st step, a procedure was developed to capture Google Earth™ elevation point data and subsequently import them into the GIS ArcDesktop™. • The 2nd step comprised of the generation of a high-resolution contour map and DEM in GIS ArcDesktop™ using the imported elevation data from Google Earth™. • In the final step, the quality and reliability of the generated
contour map underlying the DEM were assessed by comparison with standard CDSM contour data as well as high-resolution relief data generated from airborne survey data. With an elevation interval of 0.5 m the latter provides the highest resolution of all datasets and was therefore used as a benchmark. However, having been generated by Computer Aided Design (CAD) software, the data first had to be converted into a GIS-compatible format. A schematic overview of the process of generating a highresolution DEM, in the form of a flow diagram, is depicted in Fig. 7.
Google Earth image of study area with elevation data at 1m/ 1ft resolution Grid-based capturing of Google Earth elevation data at 1m intervals
Import of unprocessed SRTM data of 5°-grid cell into GIS Arc Map
Export of captured data from Excel to GIS (Arc Map) Generation of contour map for study area of 1 m resolution (‚1m Google‘)
Quality assessment: Inter-comparison of generated ‚1m Google‘ map with other contour sets
Cutting out of
SRTM data
relevant for study area
Airborne survey of study area (CAD – format)
CDSM contours in topographic map 1:50000 of study area
Conversion of CAD data into GIS format
Generation of 1 m – interval contour map of study area
5 m interval (digital format only)
20 m interval (‚20m CDSM‘)
(‚5m CDSM‘)
(‚1m SRTM‘)
Hardcopy (jpg format)
Digital (vector format)
Generation of 0,5 m – interval contour map of study area (‚0,5m Survey‘)
Figure 7 Flow path diagram depicting the data sources and different steps of generating a high resolution DEM/contour map based on Google EarthTM elevation data and its subsequent comparison to other contour datasets as a form of quality assessment (shaded blocks: contour sets that were compared with each other). For results of the inter-comparisons of the 6 different contour datasets see Table 2
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
57
27°9’27“
3460 m
26°27‘48“
27°7’22“
8“ 6“
ive r
N
GMB eye
3590 m
Mo oi R
Grid point: E10
26°29’45“
GMB wetland
3470 m
Figure 8 Google Earth™ exported image displaying the extent of target area around the Gerhard Minnebron (GMB) wetland and grid for which elevation data at 1 m interval were extracted
Import of Google elevation data As a 1st step the extent of the target area had to be determined, taking possible hydraulic links between the wetland and the adjacent stream systems, especially overland flow during high flow conditions, into account. This resulted in a target area exceeding the actual size of the Gerhard Minnebron wetland (Fig. 8). Using the ‘display grid option’ of Google Earth™ this target area was overlain by a grid of markers, with 6″ longitudinal and 8″ latitudinal difference between points resulting in a total of 16 columns and 20 rows (Fig. 8). The columns were named with
capital letters in alphabetical order from west to east, starting at the column intersecting the confluence between the Gerhard Minnebron stream and the Mooi River. Later 4 columns were added to the west named W to Z. Together with the rows numbered 0 to 19 from south to north each point of the grid received thus a unique alpha-numerical ID-code (Fig. 8). Each grid point was marked by a so-called ‘placemarker’, a Google EarthTM function in which additional information can be recorded and displayed on screen. Using this ‘placemarker’ function for each grid point, the elevation was displayed when the cursor was held over the point in question and manually recorded. In addition to the elevation of the grid points, every point at which elevation changed (by 1 m) was logged and captured together with its coordinates in a Microsoft Excel™ spreadsheet. Changes in elevation between grid points were established by slowly moving the cursor along the longitudinal and the latitudinal grid lines until the elevation reading changed. This was frequently done through zooming in to obtain a more accurate spatial reference of the point in question. All elevation points were placed between 2 adjacent grid points (horizontally, east-west, or vertically, north-south) and named by combining the code of the 2 points separated by an arrow indicating the direction of cursor movement. Requiring an estimated 30 h, a total of 2080 elevation points were captured covering an area of some 12.44 km² (equalling some 2.5 man-hours per km² at a rate of some 70 points/h). At an hourly rate of ZAR300 this would result in ZAR700/km². After inclusion of proportional costs for the internet connection and the licence fees for the GIS software, this is over 10 times more cost-efficient, for this area, than the high-resolution, airborne survey. Points where the elevation changed by 1 m in a W-E direction are grouped as ‘horizontals’ while points changing in a N-S direction are termed ‘verticals’ (Table 1). Cases where elevation changes of more than 1 m between 2 neighbouring grid points occurred resulted in 2 or more points having the same point name. However, since x and y coordinate values are provided for a point they cannot be confused. The coordinates obtained from Google EarthTM had to be converted from the degree-minute-second format into decimal
Table 1 Example of Microsoft Excel™ table used for recording information associated with elevation points for the GMB study area as retrieved from Google Earth (total length: 2 080 rows) Elevation [m a.m.s.l.]
1399 1400 1400 1401 1409 1410 1411 1412 1408 1407 1409 1417 1420 1415 1400 1398
58
South [DDMMSS]
East [DDMMSS]
Point Name
Point Type
DeciDegSouth
DeciDegEast
262945 262939 262933 262927 262842 262836 262830 262824 262818 262812 262806 262800 262754 262748 262945 262945
270752 270752 270752 270752 270722 270722 270722 270722 270722 270722 270722 270722 270722 270722 270726 270745
A_0 A_1 A_2 A_3 W_10 W_11 W_12 W_13 W_14 W_15 W_16 W_17 W_18 W_19 W_0 -> X_0 Z_0 -> A_0
Grid points Grid points Grid points Grid points Grid points Grid points Grid points Grid points Grid points Grid points Grid points Grid points Grid points Grid points Horizontals Horizontals
-26.49583333 -26.49416667 -26.49250000 -26.49083333 -26.47833333 -26.47666667 -26.47500000 -26.47333333 -26.47166667 -26.47000000 -26.46833333 -26.46666667 -26.46500000 -26.46333333 -26.49583333 -26.49583333
27.1311111 27.1311111 27.1311111 27.1311111 27.1227778 27.1227778 27.1227778 27.1227778 27.1227778 27.1227778 27.1227778 27.1227778 27.1227778 27.1227778 27.1238694 27.1291111
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
degrees to be usable in ArcView GIS. This conversion was executed as follows:
(IDW) technique of ArcDesktop™ to interpolate between these points a 1 m interval contour map of the study area was generated (Fig. 9). Based on the retrieved set of elevation data from Google EarthTM a DEM was generated in ArcDesktop™ allowing for a 3-D-view of the study area. Enhancing elevation data through colour-coding aided the identification of potential pathways for water from the Mooi River possibly entering the peatland during flood events as well as additional outflow (seep) points (Fig. 10, next page).
a.. DeciDegEast = deg + (minute/ 60) + (second/ 3600) b.. DeciDegSouth = - deg - (minute/ 60) - (second/ 3600) where: DeciDegEast and DeciDegSouth are the decimal degree coordinates for the latitudes and longitudes, respectively. Newer versions of Google Earth™ tools/options show how latitude/longitude can be changed, sparing the user from manually calculating these coordinates.
Assessing quality of the generated contour map
It should be noted that the resolution could be increased by a factor of 3, to 0.3 m contour intervals, by choosing in Google EarthTM to display elevation in feet instead of meters. Since this would have resulted in approximately triple the amount of time spent on capturing the elevation points, only changes at 1 m intervals were used in this pilot study.
In order to assess the reliability of the generated contour map and the associated DEM of the study area, the contour information was compared against 2 sets of standard topographic data (20 m and 5 m contour intervals from CDSM). The latter 2 sets were also compared to each other to assess how much same-source contour data may deviate from each other. This, in turn, provides some baseline as to what degree of deviance is regarded as acceptable by the CDSM, and was thus used to assess the quality of the generated 1 m Google contour data. In addition, elevation data from an aerial photography survey later became available to the project. This dataset provided the study with, among other features, 0.5 m contoured elevation data (acquired digitally from DWAF, 2008). Having the highest resolution this dataset was used as a benchmark against which all other contour datasets were compared. Generated in a computer-aided design (CAD) program, digital elevation data underlying this map had to be converted to be used in the ArcView GIS. This procedure, as well as pertinent technical details of the airborne survey, is briefly outlined below.
Generating contour maps and DEM in GIS ArcDesktop™ Out of the 2 080 points logged, only points where the elevation displayed in Google EarthTM changed, e.g. from 1 410 m to 1 409 m (in total 1 760 points, termed ‘horizontals’ or ’verticals’) were used to generate the contour map, resulting in a density of 1 elevation point per 7 069 m². The remainder are grid points that were omitted since their elevation could not be determined as accurately, as most of them fall somewhere in between points where elevation changes. The use of the 320 grid points would therefore have introduced an error margin of up to 1 m per point. Using the Inverse Distance Weight 27°7'20"E
27°7'40"E
27°8'0"E
27°8'20"E
27°8'40"E
27°9'0"E
27°9'20"E
1405
26°28'20"S
1420
1415
5 143
1410
1410
1410
5
10
1400
1400
Canal
1:20,000
26°29'20"S
1435
1440
Legend
143 0
14 05
0 14 1
Original image in Google Earth displaying the grid from which a total of 1 760 elevation points were retrived underlying the map displayed on the right hand side
26°29'0"S
14 0
0
14 0
26°28'40"S
14
05 14
1m-Google contour
250
500
1,000 m
0 14 0
Perennial stream 142 5
0
Wetlands
26°29'40"S
Non-perennial stream
Figure 9 Map of the study area displaying 1m interval contours generated through IDW (inverse distance weight) interpolation in GIS ArcDesktop™ based on 1 760 point elevation data retrieved from Google Earth™
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
59
N Possible inflow pathways during flood events in the Mooi River
Figure 10 Colour-coded DEM (based on IDW interpolation) of the wider Gerhard Minnebron wetland area based on 1 m interval Google Earth™ elevation data (from the SRTM) showing possible eventdriven flow paths between the Mooi River and the wetland
Mooi River
Du Toit Spruit
GMB wetland
X +2 929 200
Legend
X +2 929 200 14 04
1407.5
Y - 14 800
14
Y - 14 700
Y - 14 600
Y - 14 500 1403
Y - 14 400
1405
Y - 14 300
.5
Y - 14 200
Y - 14 100
5 .0 140
140 3
.0
5 1404. 14
05 .5
14
14
04
.0
14
03
.5
05 .0
3 .0 140
Fallow Land 14
07
Cad Renderer
1404.5
Rocky
.0
14
02
.5 1402.0 14 01
Ruin
14
.5
Res
01
14
14
00
Ruin
.5 1407.5
1404.0
1399
.0
1399.0
.5
1398
.5
98 .0
Fallow Land
1:400,000,000
2 0065
00
13
W/ Tank
Ruin
Solid, 6, 0
X +2 929 300
.0
X +2 929 300
d6.dgn Polyline
1406.5 1406. 0
1406.5
Rocky
Ruin
07
Y - 14 900
1404.5
.0
1405.5
1405.0
1406.0
Solid, 4, 2
1405
Ruin
Ruin 140 6
Solid, 2, 0
140 6
Solid, 0, 0
139 6
1404
.0
14
X +2 929 400
.5 03
X +2 929 400
Medium-Dashed, 0, 0
.0 14
04
1404.5
1405.0
5 .5 140
.0
07
1406
14 07
.0
08
.5
14
06
14
14
.5
.0
.0 1397 8 .5 140
139 8
1396
13
Solid, 1, 1
96
1408
5 1404.
.5
140 2
08
1402
1405.5
.0
14
07
.5
14
14
06
Land
.0
13
14 6.5 140
06
5 .5 140 .0
14
05
.5 03
96
.0
.0 14 04.
14
5
14
04
.0
14
06
LongDashed, 0, 1
.5
03 .0
1402.5
140 5
14
06
Short-dashed, 56, 0
.0
1405.5
14
03
14 .5
03
.014
02
.514
02
14 .0
01 .5
14
01
.0
14
.5 14
00
.0
13 99
99 .5
2 0066
00
13
140 5.0
14
01
Solid, 0, 1
1404.5
.5
14
04
.0
.5
140 4
X +2 929 600
1401.0
Rocky
0.5 140
1397.5
00
8.5 139 1398.0
1397.0
1396
.5
14
Rocky
1399
140 2
.5
1396
.0 5.5 139
140 2
13 98. 5
139 9 .0
1403.0
14 03 .5
.5
1398
1396.0
Y - 14 300
Y - 14 200 R.L. FROM
STATION
G.T.S. BEA CONS
Y - 14 400
R.L.
o
SHEET LAYOUT
C LAR KE27 1880 1492.80
Y
Y - 14 500
CO-ORDINAT E SYSTEM LO.
ME TRES ABO VE M.S.L.
X
X
Y - 14 600
o
CLAR KE 1880 27 1492.80
Y
C5
METRES ABO VE M.S.L.
D5
R.L.
Y - 14 100
F3
E5
E6
E7
F5
F6
F7
G5
G6
G3
G4
H3
H4
Gauteng Region
SURVEYED
H5
COMPILED
J5
Siyazi DTM Services
PHOTOGRAPHY
10
20
40
60 SCALE 1/1000
80
100
120
140
m
No.
01/'08
TRACED
04/'08
TOPOGRAPHY CHECKED
W611 03/'08 REVISION
Siyaz i DTM Services
MOOI RIVER G.W.S
04/'08
CONTOUR SURVEY OF GERHARDMINNEBRON PEAT STUDY
CADASTRAL CHECKED DATE
DISTRICT AREA
POTCHEFSTR OOM 40
Ha
CODE
SHEET SAR
CON
DSM
Aerial photography survey elevation map dated March 2008 Generated by Computer Aided Design (CAD) software (.dgn-file type) Source: DWAF (2008)
Apart from the visual comparison employed in this paper there are other, more formalised and quantifiable, procedures such as those used by Gorokhovich and Voustianiouk (2006) and Sun et al. (2003). However, in view of the associated debate regarding the best method, and the fact that visual comparison proved to be sufficient for the purpose intended, no further efforts were made to quantify deviations between the different contour sets in a formalised manner.
60
D6
Short-dashed, 2, 0
Figure 11 An example for 1 out of 22 separate CAD design files covering an 800 m long and 500 m wide part of the surveyed study area
Dotted, 0, 0
DEPARTMENT OF WATER AFFAIRS & FORESTRY
PRETO RIA 0001
D7
E4 F4
J4
0
RE PUBLIC OF SOUTH AFRICA
DIRECTORATE SPATIAL AND LAND INFORMATION MANAGEMENT PRIVATE BAG X313
C6 D6
Y - 14 700
G.T.S. BEACONS
Y - 14 800
STATION
X +2 929 700
Y - 14 900
CO-OR DINAT E SYSTEM LO. R. L. FROM
Short-dashed, 0, 0 Medium-Dashed, 6, 0
1398 .0 1397.5 1397.0 .5 1396
1396
X +2 929 700
Solid, 1, 2
140 3
.0
1399.0 139 9
95
14
140 2 .0
Land
140 6
Reeds
13
X +2 929 500
140 6
W/Tank
Concr.
140 0
X +2 929 600
Short-dashed, 1, 0
.0
07
X +2 929 500
Solid, 1, 0
1405.0
14
139 6
LOC. NO.
C230/41
REG. N O.
155335/08
Dotted, 0, 1 Solid, 0, 2 Dotted, 0, 2 Solid, 3, 2 Solid, 8, 0 Solid, 8, 2
Airborne relief mapping (0.5 m contour interval) The mapping is based on a stereoscopic evaluation of a series of aerial photographs taken from an aircraft at a height of 4 000 m above ground level, accurate enough for use in 1:2 000-scale plans (Smith, 2009). Twelve ground-control points of known elevation and coordinates were employed to geo-reference the images acquired from the aircraft. Geo-referenced stereoscope images were produced and, by using these images and a stereoplotter, contours and other features were drawn in (Smith, 2009).
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
27°7'40"E
26°28'20"S
0
250
500
27°8'0"E
27°8'20"E
27°8'40"E
27°9'0"E
Figure 12 Map with 0.5 m contours generated from the aerial photography survey data displaying reed cover and open water areas
1,000 m
Legend 0.5m-Survey contours in 1m intervals Area covered by reeds Area covered by open water
26°28'40"S
Non-perennial stream
Perennial stream
26°29'0"S
1:15,000
27°8'30"E
27°8'50"E
1:2,500
14 0
140 5
7
6
0 14
26°28'40"S
26°29'20"S
27°8'40"E
1397
04 14
26°29'40"S
14 03
1403
1402
1401
140 0
9 139 1398
1397
1396
1396
26°30'0"S
6 139
1396
26°28'50"S
1395
Legend 0.5m-Survey contours in 1m intervals
1397
Area covered by reeds Area covered by open water
0
62.5
125
250 m
The data from the airborne survey were then digitised data and mapped topographic features was detected; these using a CAD program locking all feature information in a misfits are possibly related to inaccuracies in geo-referencing single CAD file (design files, .dgn). Within these files colour and projection. However, for the purpose of the project, the and cartographic symbols are used to identify different feamisfits were regarded as tolerable. The DEM produced from tures. An example of such a file for a segment of the study area the CAD dataset differs from the Google-based DEM in that it is displayed in Fig. 11. was based on CAD-generated contours and not on interpolated The conversion of CAD data into GIS-compatible shape files elevation points. Figure 13 shows the CAD-based DEM of the was complicated by the fact that for symbols consisting of more GMB wetland also displaying some of the features captured by than 1 part (such as a dashed line which is made up of many the aerial survey. small lines), each separate part of that symbol is identified in GIS as a topographic feature. While N these single parts are still automatically grouped into one larger shape file (e.g. a catchment boundStructures ary), such shape files cannot be used for determinGMB Rock outcrop Wetland ing enclosed areas. This is only possible once Bushes each part of the shape file is connected to others Contour lines forming one larger entity (using among other techniques the ’merge function’ in the ‘Editor’ tool of ArcDesktop™). Where CAD lines were enclosing areas such as mining ponds, wetlands, etc., of which the area size had to be determined, this connection of separate symbol parts had to be perCa na formed, rendering the import of CAD data rather l Mooi R time-consuming. Before the area enclosed by such iver lines could be determined, all these lines had to be linked together to form an uninterrupted line Farmland around the area which could then be transformed into a polygon using the ‘buffer’ and ‘union’ functions in ArcDesktop™. Only after the poly500 m gons are created, can the size of enclosed areas be determined. Examples for such areas (wetlands and mining ponds) are displayed in Fig. 12. Figure 13 Compared to the 1:50 000 topographic map, DEM based on aerial photography survey contour data (0.5 m intervals) which spatial deviation between the GIS-imported CAD include features such as rock outcrops, structures, vegetation, rivers and canals
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
61
27°9'0"E
27°9'30"E
0
27°10'0"E
250
500
1,000 m
0
250
500 m
1:20,000
26°28'20"S
26°28'30"S
1440
26°28'0"S
1:5,000
26°28'10"S
27°8'30"E
26°29'0"S
1440
27°9'40"E
27°9'50"E
Legend
Figure 14 Mismatch between same-source data: ‘Original contours’ being the contours before mismatch was reset and ‘CDSM fitted contours’ being the contours set to the form and location of the contours found on the national topographic maps (1:50 000)
26°29'30"S
Fitted countours Original contours Non-perennial stream Perennial stream
27°7'30"E
27°8'0"E
27°8'30"E
27°9'0"E
27°9'30"E
ins r
nte
i ve
rfo
oi R
26°28'0"S
Wo n
de
Mo
pru
it
27°7'0"E
To it S
438 m
Du
Max. horizontal deviation
were subsequently taken as the base dataset to adjust the digital (vector) contour lines (Fig. 14). The adjusted set of the digital 20 m interval contour data was subsequently used for the following comparisons.
Deviations between samesource data: 20 m CDSM vs. 5 m CDSM contours Owing to the generally flat topo graphy of the wetland only 3 x 20 m interval contours occur within the larger study area and can thus be compared to the 5 m interval contours. Figure 15 depicts the deviations between 2 sets of contour lines both originating from the CDSM. Legend In assessing the deviation Non-perennial stream between contour lines, it is imporPerennial stream tant to note that identical degrees of 5m-CDSM contours deviation in steep terrain result in 20m-CDSM contours generally smaller horizontal disWetlands 0 250 500 1,000 m tances between contour lines (i.e. ‘misfits’) than in flat terrain where the lower topographic gradient (the Figure 15 height difference kept constant) Deviation areas between 20 m and 5 m CDSM contours illustrated using the 1 400, 1 420 and 1 440 m a. m. s. l. contour lines of both sets results in a larger horizontal distance between contour lines. With some 438 m the largest horizontal misfit occurs at the 1400 m Comparing Google and CDSM data contour in the (flat) confluence area where the Du Toit Spruit Before comparing 20 m interval CDSM contour data with joins the Mooi River (Fig. 15). the 5 m CDSM contours it was noted that the contour data in vector format (i.e. discrete .shp files) showed some deviation 1 m Google contours vs. 20 m CDSM contour data from identical contours as displayed in scanned images of To compare the 1 m contour map (‘Google contours’), only the 1:50 000 topographic maps imported into ArcDesktop™. The 2 x 20 m CDSM contour lines (‘CDSM contours’) which occur topographic image’s contour lines (as the original data source)
62
1420
1440
26°29'30"S
1440
1400
1420
1400
1400
26°29'0"S
1400
26°28'30"S
pr ui t
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
27°8'30"E
27°9'0"E
sp ru
it
27°8'0"E
26°28'0"S
Wo n
de rfo nte in
R iv er
Max. h ru it
orizo n
ta 540 m l deviatio n
14 05
Non-perennial stream 20m-CDSM Countours
140 0
00 14
Perennial stream
14 20
1m-Google contours Wetlands
0
142 Go 0 fits ogle 14 25 CD 142 SM 5 c o 14 20 ntou co r nt o ur
5 140
Legend 26°29'30"S
00 14 140 5
Sp
140 0
To it
Go fits ogle CD 140 SM 5 c 14 ont 00 o co ur nt o 1405 ur
26°29'0"S
26°28'30"S
Du
250
500
1,000 m
Figure 16 1 m Google contours compared with 20 m CDSM contours illustrating an overall good fit and a +5 m off-set in Google data using the 1 400 m and 1 420 m a. m. s. l. contour lines indicating a maximum horizontal misfit of 540 m between the 2 sets of contour lines 27°7'30"E
27°8'0"E
27°8'30"E
27°9'0"E
27°9'30"E
t ei ns
ive
r
on
oi R
26°28'0"S
Wo n
de rf
Mo
pr u
it
27°7'0"E
Max. horizontal deviation 490 m it
5 140 00 14
00 14
Perennial stream
10 14
1m-Google contour 5m-CDSM contours Wetlands
0
250
500
1,000 m
20 14
30 14
143 5
Non-perennial stream
14 25
14 05
Legend
140 5
pr u
14 05
To it S
26°28'30"S
Du
14 15
Comparison of all contour data to a high-resolution airborne survey: 0.5-m-contour interval The survey was done at 2 different levels of accuracy with the GMB wetland itself mapped at a scale of
27°7'30"E
Mo oi
26°29'0"S
1-m-Google contours vs. 5-m-CDSM contour data Comparing 1 m Google contours with 5-m-CDSM contours once again showed an overall good fit between the 2 datasets (Fig. 17). As observed in the 2 previous comparisons the largest horizontal deviation occurs again on the 1 400 m contour line next to the confluence of the Mooi River and the Du Toit Spruit (490 m) (Fig. 17). As demonstrated earlier, this is most likely not due to an above average error in either of the datasets but rather due to the fact that this area is particularly flat resulting in larger horizontal misfits than identical misfits in steeper terrain.
27°7'0"E
26°29'30"S
in the area of interest could be used, namely the 1 400 m and the 1 420 m a. m. s. l. lines. Comparing these with the corresponding contour lines, from the generated Google map, suggest an overall satisfactory fit between the 2 sets of contours (Fig. 16). However, Google contour lines matched best with CDSM contours which generally displayed a 5 m lower elevation (Fig. 16). The reason for this 5 m off-set could not be determined. It is, however, relatively easy to adjust by subtracting 5 m from all contour heights in the generated Google map. As was the case with the 5 m vs. 20 m CDSM contour data, the largest horizontal deviation (540 m) occurs at the 1 400 m contour at the confluence area of WFS and Mooi River, this time, however, extending into the Du Toit Spruit area (Fig. 15). Based on the maximum horizontal deviation it appears that the generated 1 m-Google map fits well to the 20 m CDMS contours (with a maximum horizontal deviation of 540 m) while the 5-m-contour set fits to the 20 m CDMS contours (also created by the CDSM) with a maximum horizontal deviation of 438 m. The fact that the generated map based on Google data displays an even better fit to the standard 20 m CDSM contour than the 5 m contour generated by the CDSM themselves indicates a better-than-expected reliability of the generated Google map.
Figure 17 Selected 1 m Google contours (bold dashed line) compared against corresponding 5 m CDSM contours (5 m off-set to be considered)
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
63
1: 4 000 while the wider study area extending to the confluence of the WFS and the Mooi River was captured at a scale of 1:8 000. However, the contour interval resolution for both areas was kept at 0.5 m. The location of the 2 areas mapped using airborne survey data is indicated in Fig. 18.
For the purpose of comparing contour lines only the area mapped at higher accuracy (1: 4 000) covering the actual wetland is used. The 0.5 m contour interval map produced for the study area is displayed in Fig. 19. The noticeable lack of contour lines in the actual wetland illustrates the flatness of the area, where differences in the micro-relief driving the development of flow patterns are frequently well below the 0.5 m interval. It should also be noted that the confluence area of the Du Toit Spruit and the Mooi River falls outside the 0.5 m part of the survey (Fig. 19). For the subsequent comparison with contour data of alternative origin, only those contour lines of the 0.5 m set which correspond with the contours of the different datasets falling into the surveyed area are displayed. 20 m CDSM contours: Measured against 2 contours from the 20 m CDSM contour dataset (1 420 m and 1 440 m a. m. s. l.) the aerial photography survey contours require no off-set correction displaying an overall good degree of congruence between the 2 contour sets (Fig. 20). Owing to the limited area for which the survey at this resolution was done, only parts of the 20 m contours can be compared against corresponding 0.5 m survey contours. For the area covered, a comparatively small horizontal deviation of a maximum 150 m at the 1 400 m contour line indicates an overall good fit between the 2 datasets (Fig. 20, next page).
Figure 18 Topographic map of the study area depicting the areas surveyed by airborne photography using 2 different scales of mapping (1:4 000 for the wetland – Area A and 1:8 000 for the wider study area, Area B) (DWAF, 2008)
5 m CDSM contours: As Fig. 21 (next page) indicates, the 0.5 m survey contours and the 3 relevant 5 m CDSM contours (i.e. 1 400 m, 1 405 m, 1 410 m) match well. 27°7'30"E
27°8'0"E
27°9'0"E
ru sp ein
Ri v er
26°28'0"S
1401
W on
de
rfo
26°28'40"S
Mo oi
nt
14 12
4 140
27°8'30"E
it
27°7'0"E
1398
26°28'50"S
26°28'30"S
Du
To it S
pr u
it
141 6
141 2
140 8
26°29'0"S
1404
142 0
1401
Non-perennial stream
Canal
Perennial stream
Wetlands
1:2,500
0.5m-Survey contours 0
25 50
100 m
Figure 19 0.5 m survey contours for the wider study area as well as the wetland part
64
27°9'10"E
26°29'30"S
27°9'0"E
26°30'0"S
27°8'50"E
Legend
Legend
0
250
Non-perennial stream Perennial stream 0.5m-Survey contours 1m-Google contour Wetlands
500
1,000 m
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
27°7'30"E
27°8'0"E
27°8'30"E
27°9'0"E
27°9'30"E
ins r
nte
i ve
rfo
oi R
26°28'0"S
Wo n
de
Mo
pru
it
27°7'0"E
To it S
Max. horizontal deviation
pr ui t
150 m
26°28'30"S
Du
Figure 20 20 m CDSM contours compared against 0.5 m survey contours
1420 1400
1420
1400
0
26°29'0"S
1400
Legend
14 0
1400
1400
Non-perennial stream 26°29'30"S
Perennial stream 20m-CDSM contours 0.5m-Survey contours Wetlands
0
250
500
1,000 m
27°7'30"E
27°8'0"E
27°8'30"E
27°9'0"E
27°9'30"E
t ei ns on
r
26°28'0"S
Wo n
de rf
Mo oi Riv e
pr u
it
27°7'0"E
Du To it S
it
26°28'30"S
1405 1405
1400
1405
0.5m-Survey contours
14 1
0
Perennial stream
14 1
14 00
140 m
Max. horizontal deviation
Non-perennial stream
14 05
0
Legend 26°29'30"S
1400
14 0
1400
0
1405
140 0
26°29'0"S
Figure 21 5 m CDSM contours as solid lines and the 0.5 m survey contours as broken lines with each contour identified by a unique signature
pr u
5m-CDSM contours Wetlands
0
250
500
1,000 m
Unlike deviations in previous comparisons, the maximum deviation was not found to be associated with the lowest-lying contour line in the study area (1 400 m a. m. s. l.) but occurred at the 1 405 m contour where it amounts to some 140 m (Fig. 21). 1 m Google contours: Before comparing the 2 datasets the Google contours had to be corrected by reducing the displayed
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
elevation for each contour line by 5 m. Owing to the relatively large number of contours only the following 5 pairs of corresponding contour lines are compared with each other: 1 400 m, 1 405 m, 1 410 m, 1 415 m and 1 420 m a. m. s. l. (Fig. 22). The deviation between the 0.5 m survey and 1 m Google contours appears to be more complicated than a simple misfit, with a maximum horizontal deviation of 221 m (Fig. 22).
65
In addition there is also a sporadic difference between the shape of 1 m and 0.5 m contour lines covering the study area. The results of the various inter-comparisons of different contour sets discussed above are summarised in Table 2.
Conclusions and recommendations A common problem for hydrological research in topographically flat areas such as wetlands is the lack of sufficiently
1405
1410
26°28'30"S
27°9'0"E
27°7'0"E
27°7'30"E
27°8'0"E
27°8'30"E
27°9'30"E
27°9'0"E
26°28'0"S
221 m
Max. horizontal deviation
pru it
26°28'30"S
Du To it S
14 10
14 00 m
14 05
14 0
Legend
1400 1420
140 0
26°29'0"S
140 0
05 14
Non-perennial stream 0.5m-Survey contours
Legend
1m-Google contour Wetlands
Non-perennial stream
1,000 m
Perennial stream 1m-Google contours
Figure 22 1 m Google contours (dashed) and 0.5 m survey contours (solid) compared for selected contours
1400
26°29'0"S
62.5
125
250 m
142
1405
Wetlands 0
5
0.5m-Survey contours
1410
500
141 5
250
14 10
0
1405
26°29'30"S
Perennial stream
Table 2 Overview of results of the inter-comparison of different contour datasets available for the study area with the generated 1m Google contour map
20 m CDSM (lvector - shp)
5 m CDSM (vector – shp)
1 m Google
1 m SRTM
66
Figs. 5 and 6: - SRTM contours are blocky and often open-ended - Horizontal misfit of approx. 50 m between original 5° x 5° grid cell data and those cut out for study area - 1 400 m contour is split in 3 parts and deviates from corresponding Google contour
-
-
5 m CDSM (vector shp) Fig. 17: - A vertical off-set of 5 m was detected - After correcting this off-set an overall good fit appeared - Max. horizontal deviation: 490 m (1 400 m contour line)
-
-
20 m CDSM (image – jpg) Fig. 16: - 2 x contour lines evaluated (1 400/1 420 m) - A vertical off-set of 5 m was detected - After correcting this off-set an overall good fit appeared - Maximal horizontal deviation 540 m (1 400 m contour line) - 1 m Google contours appear somewhat smoother than 20 m CDSM contours Fig. 15: - 3 x contour lines evaluated (1 400/ 1 420/1 440 m) - Maximal horizontal deviation 438 m (1 400 m contour line) - Positive and negative deviations occur - 5 m contour line appears to be smoother than 20 m line Fig. 14: - At 1:50 000 map for study area a horizontal misfit of approx. 30 m was detected (1 400 m contour line) - Positive and negative deviations occur
0.5 m survey Fig. 22: - A vertical off-set of 5 m was detected - More complex pattern of deviation: horizontal misfit max. 221 m (1 405 m contour line) shape of contours - Deviations mainly confined to flat areas; in steep terrain generally good fits
Fig. 21: - 3 x contour lines compared (1 400/ 1 405/1 410 m) - An overall good fit - Max. horizontal deviation: 140 m (1 405 m contour line)
Fig. 20: - only 2 x 20 m contour lines could be compared - Overall good fit - Maximal horizontal deviation: 150 m (1 400 m contour line)
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
detailed elevation data required to reliably delineate catchment areas, understand flow patterns within the wetland and to identify potential hydraulic interactions between wetland and adjacent streams, to name but a few examples. GIS-based DEMs prove to be useful tools in addressing most questions. Compared to standard topographic maps displaying 20 m contour intervals, and 5 m intervals for selected areas, elevation data provided by Google EarthTM, an internet-based global visualisation application, are provided at a high enough resolution to create 1 m intervals, at a minimal cost. Although SRTM data underlie the Google Earth™ relief information, the direct use of freely-available SRTM data downloaded from the internet yielded unsatisfactory results that were found to be inferior to contours generated by the procedure described in this paper. Apart from errors associated with the interpolation of elevation data cut-out from larger sets, this may be caused by the continuous improvement of the terrain database in Google Earth™ which in addition to satellites also includes other technologies and alternative data sources. This may partly explain why the vertical accuracy of Google elevation data is higher than the 16 m indicated for the original SRTM data. Even after introducing potential errors associated with the extraction of elevation data described in this paper, the generated contours and DEMs at 1 m contour resolution were found to compare favourably with not only relief data of standard topographic maps but also results of a high resolution aerial survey. However, with little or no information available on the exact procedures and interpolation techniques used in Google EarthTM that so obviously improve the original SRTM dataset, such good agreement may be area-specific and can therefore not be unconditionally extrapolated to other regions. In addition to higher resolution relief data, Google Earth™ also provides more up-to-date hydrographical information, as a result of the close updating intervals – 2 to 4 yr compared to the approximately 20 yr intervals for the 1:50 000 topographic maps. This is of particular importance in environments where anthropogenic impacts such as peat mining result in fastchanging hydrographical conditions. The method developed to capture elevation data from Google Earth™ and import these into ArcDesktop GIS proved to be successful in creating a reliable, high-resolution DEM for the study area. After the off-set corrections of 5 m, deviations of the generated 1 m contour lines from 20 m contours contained in the 1:50 000 topographical map of the CDSM, the deviations were found to be less than those between the 2 CDSM contour datasets (i.e. 20 m and 5 m interval contours). As sources for potential errors it should be noted that, apart from some deviations between the 20 m and 5 m CDSM contour lines, we also detected discrepancies between printed contours in scanned images and contours in digital format (vector data), both provided by the CDSM. Compared to contour lines at 0.5 m intervals from a photogrammetric airborne survey, more complex patterns of deviation to the generated 1 m Google contours were found, while the (significantly less dense) CDSM contour data continued to show an overall satisfactory degree of congruence. However, with approximately ZAR10 000 per km 2 for the airborne survey the presented method of sourcing elevation data from Google EarthTM is over 10 times more cost-efficient, provided that internet access, appropriate GIS software and entry-level GIS operating skills are available. The generated 1m resolution DEM was found to significantly enhance the visual understanding of the micro-relief
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010
in the studied wetland which despite numerous field visits remained invisible to the naked eye. This, in turn, aids the understanding of how landforms may influence hydrological processes such as direction of water flow, development of drainage patterns, occurrence of groundwater discharge areas, etc. Applying the computer-generated DEM to the study area also aided the delineation of catchment boundaries and identification of presently unmapped surface depressions along which flood waters enter the wetland. Moreover, maximum peat accumulation levels could also be determined as they are controlled by the elevation of the water table at the outflow point. In summary it is concluded that the presented methodology allows the retrieval, conversion and application of elevation data from Google Earth™ in order to produce reliable GISbased DEMs that have significantly higher resolutions than the standard topographic maps. These DEMs and contour maps have proven to be particularly useful tools for hydrological research in low-relief energy environments such as wetlands for which adequate contour data are either not available, or too costly to generate.
Acknowledgements Funding of the project (Project No. 2006-231) and logistical support by DWAF (now DWA) is gratefully acknowledged. The researchers would also like to give heartfelt thanks to Ms S Rudolph and the DWA agents at the Boskop Dam and Aquapark offices for their patience, help, data and other saving graces, as well as to T de Klerk and D Cilliers for their input. We also wish to thank an anonymous reviewer, in particular, for her/his encouraging and constructive comments.
References ALEMAW BF and CHAOKA TR (2003) A continental scale water balance model: a GIS-approach for Southern Africa. Phys. Chem. Earth 28 957-966. AL-SABHAN W, MULLIGAN M and BLACKBURN GA (2003) A real-time hydrological model for flood prediction using GIS and the WWW. Comput. Environ. Urban Syst. 27 9-32. AUBIN M (2009) Google Earth: From Space to Your Face…and Beyond. URL: http://www.google.com (Accessed 17 August 2009). BOBBA AG, SINGH VP and BENGTSSON L (2000) Application of environmental models to different hydrological systems. Ecol. Modell. 125 15-49. CDSM (2007) Personal communication. 28 June 2007. Enquiry 4253/07. E-mail from
[email protected] to E Hoffmann. DROGUE G, PFISTER L, LEVIANDIER T, HUMBERT J, HOFFMANN L, EL IDRISSI A and IFFLY JF (2002) Using 3D dynamic cartography and hydrological modelling for linear streamflow mapping. Comput. Geosci. 28 928-994. DWAF (2008) Topographic map 1:50 000 sheet 2627 AC (Rysmierbult), 1995 indicating the flight direction and number of strips of the airborne photography survey. Department of Water Affairs and Forestry, Pretoria, South Africa. FARR TG, ROSEN PA, CARO E, CRIPPEN R, DUREN R, HENSLEY S, KOBRICK M, PALLER M, RODRIGUEZ E, ROTH L, SEAL D, SHAFFER S, SHIMADA J and UMLAND J (2007) The Shuttle Radar Topography Mission. URL: http://www2.jpl. nasa.gov (Accessed 30 June 2008). DOI:10.1029/2005RG000183. FRASER CJD, ROULET NT, and LAFLEUR M (2001) Groundwater flow patterns in a large peatland. J. Hydrol. 246 142-154. GOROKHOVICH Y and VOUSTIANIOUK A (2006) Accuracy assessment of the processed SRTM-based elevation data by CGIAR using field data from USA and Thailand and its relation to the terrain characteristics. Remote Sensing Environ. 104 409-415.
67
JANAUER GA (1996) Macrophytes, hydrology, and aquatic ecotones: a GIS-supported ecological survey. Aquat. Bot. 58 379-391. JANSSEN R, GOOSEN H, VERHOEVEN ML, VERHOEVEN JTA, OMTZIGT AQA and MALTBY E (2005) Decision support for integrated wetland management. Environ. Modell. Software 20 215-229. JEWITT GPW, GARATT, JA, CALDER, IR and FULLER L (2004) Water resources planning and modelling tools for the assessment of land use change in the Luvuvhu Catchment, South Africa. Phys. Chem. Earth 29 1233-1241. KO C and CHENG Q (2004) GIS spatial modeling of river flow and precipitation in the Oak Ridges Moraine area, Ontario. Comput. Geosci. 30 379-389. LACROIX MP, MARTZ LW, KITE GW and GARBECHT J (2002) Using digital terrain analysis modeling techniques for the parameterization of a hydrologic model. Environ. Modell. Software 17 127-136. MINNAAR A (2009) Personal communication. Department of Water Affairs and Forestry, Sub Directorate: Remote Sensing, Pretoria. PORTOGHESE I, URRICCHIO V and VURRO M (2005) A GIS tool for hydrogeological water balance evaluation on a regional scale in semi-arid environments. Comput. Geosci. 23 15-27.
68
RAPER JF and KELK B (1991) Three-dimensional GIS. In: Maguire DJ, Goodchild MF and Thind DV (eds.) Geographic Information Systems: Principles and Applications. Longman: London. 1 299-316. SMITH N (2009) Personal communication. Department of Water Affairs and Forestry, Potchefstroom, South Africa. SOULSBY C, TETZLAFF D, RODGERS P, DUNN S and WALDRON S (2006) Runoff processes, stream water residence times and controlling landscape characteristics in a meso-scale catchment: An initial evaluation. J. Hydrol. 325 197-221. SUI DZ and MAGGIO RC (1999) Integrating GIS with hydrological modelling: practices, problems, and prospects. Comput. Environ. Urban Syst. 23 33-51. SUN G, RANSON KJ, KHARUK VI and KOVACS K (2003) Vali dation of surface height from shuttle radar topography mission using shuttle radar altimeter. Remote Sensing Environ. 88 401-411. TAYLOR F (2009) NASA/Japan Release Terrain Data – Google Earth’s is Better URL: http://www.gearthblog.com (Accessed 17 August 2009). WINDE F (2006) Challenges for sustainable water use in dolomitic mining regions of South Africa – a case study of uranium pollution. Part I: Sources and pathways. Phys. Geogr. 27 333-347.
Available on website http://www.wrc.org.za ISSN 0378-4738 (Print) = Water SA Vol. 36 No. 1 January 2010 ISSN 1816-7950 (On-line) = Water SA Vol. 36 No. 1 January 2010