12 August 2016 Developing an algorithm for enhancement of a digital terrain model for a densely vegetated floodplain wetland
Author Affiliations +
J. of Applied Remote Sensing, 10(3), 036013 (2016). doi:10.1117/1.JRS.10.036013
Abstract
Airborne laser scanning survey data were conducted with a scanning density of 4  points/m2 to accurately map the surface of a unique central European complex of wetlands: the lower Biebrza River valley (Poland). A method to correct a degrading effect of vegetation (so-called “vegetation effect”) on digital terrain models (DTMs) was applied utilizing remotely sensed images, real-time kinematic global positioning system elevation measurements, topographical surveys, and vegetation height measurements. Geographic object-based image analysis (GEOBIA) was performed to map vegetation within the study area that was used as categories from which vegetation height information was derived for the DTM correction. The final DTM was compared with a model obtained, where additional correction of the “vegetation effect” was neglected. A comparison between corrected and uncorrected DTMs demonstrated the importance of accurate topography through a simple presentation of the discrepancies arising in features of the flood using various DTM products. An overall map classification accuracy of 80% was attained with the use of GEOBIA. Correction factors developed for various types of the vegetation reached values from 0.08 up to 0.92 m and were dependent on the vegetation type.
Mirosław-Świątek, Szporak-Wasilewska, Michałowski, Kardel, and Grygoruk: Developing an algorithm for enhancement of a digital terrain model for a densely vegetated floodplain wetland

1.

Introduction

Airborne laser scanning (ALS) is widely applied to retrieve data on topography that is used to develop digital terrain models (DTMs). ALS remains as a cutting-edge methodology for retrieving high-resolution data on land relief and has multiple advantages over traditional surveying techniques. It is considered to be the most time-efficient method and advanced technology due to the capability to derive topographic (bare ground elevation) data by filtering vegetation or other objects on the surface and the capability to produce centimeter resolution DTMs over large areas.1 Within near-natural floodplains, dense vegetation can overestimate ALS terrain measurements by masking true terrain elevation. In this regard, the quality of ALS-based DTMs of vegetated and relatively flat areas, such as lowland floodplains, requires special scientific attention.

Floodplain topography and its complexity influence water flow and the pattern of floodplain inundation, creating spatial and temporal patterns related to the processes of erosion, sediment transport, and deposition.23.4.5.6 The (micro)topographic heterogeneity of riparian habitats shapes unique hydrological conditions present in riverine landscapes that determines biodiversity.7,8 These conditions are relative to the type of water supply and can be considered catchment-scale variables. Topography of a floodplain underpins the extent of floods and the degree of hydration in particular riparian habitats, remaining a strongly local but important element of these geosystems. The importance of the exchange of water, sediments, nutrients, organic matter, and biota between the river and the floodplain has been emphasized in many different studies.910.11.12 Thus, it has been proven that analyzing ecohydrological feedbacks in floodplain wetlands with respect to terrain relief and elevation, on top of the broad spectrum of hydrological data, requires appropriate data on floodplain surface topography.12,13

Most of the topographic analyses in hydrological sciences are performed using DTMs created from the data retrieved either from traditional ground-surveyed optical measurements or, more commonly, remote sensing data, ALS in particular. High quality of ALS-based DTMs that remain as digital representation of bare terrain (surface elevation of the ground without any of the structures or vegetation), significantly improved the quality of hydraulic modeling, fluvial geomorphology analyses, and river dynamic studies.1415.16 The quality of topographic data for floodplains is critically important because only small changes in elevation across these low relief landscapes and small local differences in elevations can entail different functions of these ecosystems and, hence, should be accurately represented by topographic data. However, when variations in DTMs are the result of errors introduced from the chosen data processing method and do not represent the real terrain features, such DTM data may introduce significant errors in habitat assessment and, likely, bias conclusions.

The traditional approach has some limitations related to the high costs of field measurements and/or terrain inaccessibility. Topographic contour maps, especially the ones produced in mid-20th century (in our case, these maps were used in the analysis as they are the only available source of data on land relief available in the appropriate scale of 110,000), also face technical limitations. The uncertainty in elevations introduced from large spacing between contours in flat areas (such as floodplains) and ground measurement control points negatively influence the DTM products obtained through interpolation. By contrast, high spatial resolution DTMs derived by ALS have been shown to be more representative for field slope measurements17 and field-determined elevations18 than DTMs created using topographic contour maps.19 They provide a comprehensive and highly accurate source of elevation data from exposed ground and vegetation surfaces. Nevertheless, while the ALS data collection and processing methods are considered as cutting-edge technology, challenges are still presented to achieve an accurate DTM for low relief floodplains with dense vegetation and periodic inundation. Hence, specific efforts are required to obtain an accurate DTM, which is appropriate for hydrological applications. This is because ALS measurements of height or backscatter are returned from dense and sometimes tall wetland vegetation, a mosaic of flooded and dry lands, the structure of trees and shrubs as well as the time of data acquisition, negatively affecting the quality of a DTM. Consequently, estimates of terrain elevation are often invalid as height measurements are exaggerated due to so-called “vegetation effect” or inundation effects. This problem has been pointed out in numerous studies. Gorte et al.20 have shown that the low vegetation (heights below 20 cm, e.g., grass or branches lying on the ground) is causing problems in determining the bare earth elevation. Hodgson et al.21 found high elevation errors in areas covered by shrubs compared to the other types of the vegetation. Ahokas et al.22 investigated the upward shift of the laser points compared to the ground points to be ±11  cm for grass and ±17  cm for forest. Bollweg and de Lange23 found an upward shift of 8 cm for long dense grass. Su and Bork24 indicated that many ALS last-return pulses may originate from the forest canopy or understoring vegetation rather than the true ground, causing overestimations of ground elevations within the extents of shrubs and forests. Hladik and Alber25 found that mean errors (MEs) of DTM-derived ground elevations for different land cover classes ranged from 0.03 to 0.25 m compared to the ground truth data, with larger offsets for taller vegetation. They also developed species-specific correction factors for selected 10 land cover classes and used these correction factors to modify the ALS-derived DTM. Będkowski and Stereńczak26 showed that the time of the year at which ALS data were acquired may considerably influence the quality of a final DTM. They concluded that elevation differences between DTMs developed on the basis of ALS data acquired in two seasons—spring and summer may be as high as a few meters. They also concluded that untypical relations between DTMs (spring–summer) are apparently connected with the presence of dense broadleaved species in understory canopies. In general, the complexity of DTM development for areas with the dense vegetation is not a new issue in the contemporary literature.20,25 However, the development of ALS DTM processing algorithms and methods needs to be revisited in terms of accuracy, especially for floodplains with high and dense vegetation.

The overall goal of the research presented in this paper is to develop an efficient process to correct ALS data for the “vegetation effect” to gain an accurate DTM that represents the actual terrain surfaces of the research area, which is a broad temperate floodplain. We contest that the methodology based on ALS data and extended to geographic object-based image analysis (GEOBIA),27 along with the application of vegetation height correction coefficients, remains an efficient algorithm for the development of accurate DTMs for floodplain wetlands. By comparing spatial analyses of selected hydrological features of a known flood event, depths of water within the floodplain, and the volume of floodplain water storage, we revealed that the DTM developed with the use of proposed methodology provides much better quality elevation data than the DTM developed with the use of standard data processing procedures.

2.

Materials and Methods

2.1.

Study Site

The northern part of the lower Biebrza River valley in Poland (Fig. 1. 53°23′49′′N—22°28′34′′E bottom left and 53°28′12′′N—22°38′9′′E upper right) was chosen as the research site due to its unique environmental features, including natural vegetation, low human pressure, a natural and dynamically changing river bed with multiple natural oxbows, and a nearly nonmodified flow regime.28

Fig. 1

The study site (the northern part of the lower Biebrza River valley).

JARS_10_3_036013_f001.png

Wetlands of the stretch of the valley analyzed are supplied with surface water and groundwater with spring thaw flooding being a main driver of floodplain ecology.29 The significance of the lower Biebrza River valley in catchment-scale water retention processes has been reported as an important ecosystem service.30 The width of the valley in the analyzed stretch reaches from 2 up to 6 km, and the length reaches 12 km. Terrain elevations within stretch of the floodplain analyzed vary from 103 m above sea level (a.s.l.) to some 107.0 m a.s.l. [Kronstadt normal height coordinate system (PL-KRON86-NH) based on GRS80 ellipsoid]. Locally, beyond the river valley, terrain elevations in the area analyzed reach some 118.5 m a.s.l. Such a low variability of terrain elevations requires special efforts to be undertaken in order to make the DTM representative for hydrological and ecological purposes. Due to the numerous oxbows and meanders of the Biebrza River, its length in the analyzed stretch reaches 18.5 km. Most of the northern part of the valley is covered by floodplain vegetation, which is mainly dependent on seasonal flooding. The presence of dense vegetation cover consisting of alder forest, willow shrubs, reeds, and sedges makes the development of an accurate DTM challenging.

2.2.

Procedure

The procedure applied in this study involved various types of data and consisted of the following tasks:

  • Task 1: Rasterization of rivers and streams vector data from the hydrographic map of Poland (MPHP) at a scale 110,000.

  • Task 2: Processing of ALS point cloud data to obtain uncorrected DTM, digital surface model (DSM), normalized DSM (nDSM), and intensity map.

  • Task 3: The object-based image classification using nDSM map to retrieve map with different height vegetation classes.

  • Task 4: Processing of airborne and satellite imagery to retrieve vegetation indices allowing to distinguish different types of land cover.

  • Task 5: GEOBIA classification of remotely sensed imagery, intensity map, and classification map from task 3.

  • Task 6: Creation of the vegetation map with respect of the vegetation height.

  • Task 7: Application of corrections factor within vegetation classes on DTM from task 2 and additional correction for shrubs and reed using extraction and interpolation.

  • Task 8: Accuracy assessment of the final DTM.

The overall scheme of the data processing steps is presented in Fig. 2. This procedure applies to additional correction method used in order to obtain accurate DTM, task 7 in areas covered by very dense vegetation such as shrubberies or reed. Considering the density of point cloud (4  points/m2) in some places the laser beam did not reach the ground. In these cases, the extraction and interpolation procedures were used.

Fig. 2

The procedure of ALS additional data processing with use of GEOBIA (Task 3 and 5—dotted line).

JARS_10_3_036013_f002.png

2.3.

Airborne Laser Scanning Survey Details

The research in the framework of this paper utilized near-infrared (1064-nm Laser Scanner model and type: Leica ALS70) ALS data acquired in autumn 2011, with a scanning density of 4  points/m2 within the nationwide Polish ISOK program (IT System of the Country’s Protection against Extreme Hazards for Poland).31

The data for analysis were gathered after the growing season, during a period when no or slight flooding or inundation within the research area was present. This is important due to the fact that surface water can contribute to inaccurate elevation measurements using ALS data. In the project, 73 tiles (1×1  km) provided by CODGiK (Central Agency for Geodetic and Cartographic Documentation) were used, covering around 7300 ha in the northern part of the lower Biebrza River valley. The point cloud was acquired with an average elevation error of 0.15 m and an average location error of 0.5 m.32 The analysis used in this study generated a 1-m gridded DTM. The data were delivered in the Polish National Spatial Reference System frames, geodetic reference frame: National Geodetic Coordinate System 1992 (ETRS89/Poland CS92, EPSG code: 2180), and vertical reference frame: Kronstadt normal height coordinate system (PL-KRON86-NH) based on GRS80 ellipsoid.

2.4.

Airborne Laser Scanning Data Processing

ALS data were delivered in Log ASCII Standard (LAS) file format by CODGiK. First, the degree of spatial detail in an unprocessed ALS data was quantified with point density and equal to 4  points/m2 and point spacing was 0.389  points/m.

The LASGround module from LASTools software (rapidlasso GmbH) was used. The LasGround module, the Wilderness algorithm, was applied as recommended for areas with natural vegetation. Nonground points (mainly trees, shrubs, reed, and sedges) were excluded from the analysis to obtain primarily bare earth points. This is important because in areas covered by high and dense vegetation the laser pulse is mainly reflected from the vegetation canopy and rarely reaches the ground. The filtering allowed high objects to be excluded from the analyzed data by taking into account only the points from last-return pulse (the lowest points measured by the laser scanner).

During the ALS point cloud data processing, several products were developed: DTM, DSM, and nDSM—Fig. 3, slope and ALS intensity maps with the spatial resolution of 1  m×1  m.

Fig. 3

DSM, DTM, and nDSM.

JARS_10_3_036013_f003.png

The nDSM was calculated by simply subtracting the last-return product (DTM) from the first-return product (DSM) to give the relative heights of vegetation.33 The ALS intensity, defined as the ratio of incoming to outgoing radiation of a laser pulse, is measured during data acquisition and in this study was used to produce the ALS intensity map. However, the resulting DTM did not give us a final reliable product for areas covered by trees and willow shrubs. In these areas, the laser pulse could not reach the ground surface due to the presence of branches, resulting in an overestimation of the ground elevation as the algorithm used in the analysis treated the lowest points as ground. The further use of GEOBIA classification in the processing procedure helped detect such areas and allowed the application of vegetation correction factors. In the end, a final corrected DTM was obtained. Moreover, in the end, vertical accuracy was determined by ground truthing after DTM creation and additional analysis was performed in order to present the importance of DTM accuracy on some crucial hydrological characteristics.

2.5.

Geographic Object-Based Image Analysis Classification

GEOBIA was used in this study to integrate ALS, optical remote sensing (including satellite and airborne products), and thematic vector data to develop an accurate map of vegetation type and height across the northern part of the lower Biebrza Basin. The classification was performed based on remotely sensed optical imagery (airborne and satellite), laser scanning data, and thematic vector data. Data processing was conducted with eCognition software (Trimble GeoSpatial).

Today, this type of image classification is used by the world’s leading research centers involved in the processing of remote sensing images.3435.36 Processing using the GEOBIA approach is based on objects—groups of pixels representing various features. In the first step of the analysis, the pixels are segmented into objects. Then, objects are assigned to land cover classes defined by the user. This process can be expanded upon and can consist of many steps, including reshaping, resegmenting, and reclassifying the initially created objects. This allows the creation of rule sets for processing large amounts of datasets without any user interaction. In this research, GEOBIA was used to process different types of spatial data (Fig. 2). In the first step (task 1), vector thematic maps, including streams, ditches, and oxbows, were rasterized. Application of the existing vector data in analysis allowed considering these objects more accurately in a final DTM. In the second step (task 2), the DTM, DSM, and nDSM were obtained from ALS point cloud data, and the segmentation and classification of the nDSM were performed in order to obtain the vegetation map on the basis of difference in vegetation height. The use of satellite (Landsat 8) and airborne images in GEOBIA (tasks 4 and 5) allowed vegetated classes to be distinguished from nonvegetated classes and helped delineate the river bed and oxbows. Multiresolution segmentation was chosen for image segmentation process. The scale parameter was set to 20. We used the default value of 0.1 for shape factor, 0.2 for smoothness, and the value 0.7 for compactness to obtain better representation of forest and shrub boundaries. These two parameters control the characteristics of similarity and heterogeneity for each image-object.

The use of three-dimensional point clouds from laser scanners, particularly in the aspect of their integration with the multispectral information originating from digital cameras, allowed the process of image interpretation to be significantly enhanced.37 In task 4, several vegetation indices were tested in order to differentiate areas covered by alder forests, willow shrubs, common reed, reed-manna grass, sedges, and grass. The chosen vegetation indices, which are presented in Table 1, included the normalized difference vegetation index (NDVI38), optimized soil-adjusted vegetation index (OSAVI39), SAVI,40 modified simple ratio (MSR41), and the vegetation vitality ratio (VVR42). All vegetation indices, (Table 1) including intensity map, were assigned the same weight during the segmentation process. Obtained image-objects were further classified into 7 land cover classes.

Table 1

Spectral vegetation indices (SVIs) used in the analysis of satellite and aerial images.

SVIEquation
OSAVINIR−REDNIR+RED+Y
MSR(NIRRED−1)/(NIRRED+1)
NDVINIR−REDNIR+RED
VVR1−a·REDNIR−b

Note: Y is the correction factor that equals 0.16; NIR is the near-infrared band; RED is the red band; and a and b are soil line coefficients.

2.6.

Digital Terrain Model Accuracy Enhancement

Developed in task 2, the ALS DTM represented the false terrain caused by the so-called “vegetation effect.” Therefore, additional correction was performed in task 7 of the presented procedure (Fig. 2) in order to decrease the influence of the vegetation on the resulting DTM. The real-time kinematic global positioning system (RTK GPS) measurements taken in the different types of vegetation within the research area allowed a calculation of the differences between the ALS DTM and elevation represented by ground control points. The positive values indicate that the DTM was above the measured ground control points. The mean value of this shift was calculated for each vegetation class and was then subtracted from the DTM in the given vegetation class in order to correct the elevation values (Table 2). The correction factors for areas overgrown by shrubberies higher than some 2.5 m and reeds higher than some 3.0 m did not allow sufficient correction of the ground elevation. An insufficient number and spatial distribution of reference points measured with the RTK GPS within patches of these vegetation types are likely to have led to unsatisfactory correction of these types. These problems occurred due to difficulties in acquiring high accuracy RTK GPS measurements (because the RTK GPS signal in these high and dense vegetation types was poor). Fortunately, these types of vegetation of the floodplain could be accurately mapped using GEOBIA. Areas covered by very dense clusters of high shrubs (represented here by Cornus sericea or willow shrubs) or dense and high reeds, the elevation data were excluded from the DTM, and interpolation from known ground points surrounding the erased vegetation patch was applied to fill in the gaps. GEOBIA classification attained a satisfactory classification accuracy (Table 2) with an overall value of 80%.

Table 2

Elevation correction factors for different types of vegetation in the northern part of the lower Biebrza Basin.

Vegetation typeNumber of GPS points (n)Total GPS and ALS error (m)Correction factor (m)SD (m)Classification accuracy (%)
Mosaic of grass and sedges830.020.080.0380
Alder forest300.030.110.0484
Loose tussock sedges310.030.280.1087
Reed-manna grass540.020.390.1183
Tussock sedges320.030.430.1767
Reed620.020.630.2082
Shrubs400.030.920.3576

2.7.

Verification of the Digital Terrain Model Processing Results and Influence of Digital Terrain Model Accuracy on Habitat Assessment

All topographic products developed from laser scanning data were verified using RTK GPS surveys performed within the research area4344.45 (Fig. 4).

Fig. 4

Location of the RTK GPS survey points within research area presented on the map resulting from GEOBIA classification.

JARS_10_3_036013_f004.png

A regular distribution of the measurement points was not possible due to the very demanding field conditions (flooding) and the lack of roads or paths allowing the penetration of the valley. These measurements were taken during surveys carried out in recent years in Biebrza National Park46 and by field measurements performed within this study. Data from the current study were collected using high-precision Topcon dual-frequency (L1/L2) RTK Global Navigation Satellite System (GNSS) receivers (GR-3 and GRS-1), and archival data were collected using dual-frequency (L1/L2) Topcon Legacy E receivers with a PGA-1 antenna. The RTK GPS surveying method was applied using the real-time NAWGEO service of the multifunctional precise satellite positioning system of ASG-EUPOS. Estimated precision of the real-time NAWGEO service is as high as 0.03 m (horizontally: XY coordinates) and 0.05 m (vertically: Z coordinate).

The RTK GPS measurements were taken using only GPS reference stations as no GLONASS and GALILEO reference stations are present in northeastern Poland.

All measurements were taken using multipath reduction and only in RTK fixed mode with position dilution of precision values lower than 4. Data points for which the GPS was not able to achieve a good GPS solution (RTK fixed solution), due to signal attenuation by the vegetation, were excluded from further use. RTK fixed solution means that the GPS can see at least five satellites in common and is receiving corrections from the base stations. The measurement accuracy in RTK mode is estimated as ±1.5  cm horizontally and ±2  cm vertically. Points measured were not referenced to any geodetic control network reference point because no such network exists in the research area.

Measurements were scattered throughout the river valley depending on topography and included 307 points taken randomly in different types of vegetation within the study area. At each location point, up to four measurements were taken, and then averaged in order to improve RTK GPS accuracy due to the lack of geodetic benchmarks in the research area. This allowed data to be obtained with an accuracy close to the one recommended by the American Society for Photogrammetry and Remote Sensing47 and the European Spatial Data Research Network48 regarding accuracy validation of ALS data. To quantify DTM error, several measures were computed (Table 3).

Table 3

Quality assessment measures used to measure the differences between values predicted by a model (DTM) (Mz) and the values actually observed (Iz).

Quality assessment measureEquation
MEME=1n∑i=1n[Mz(x,y)−Iz(x,y)]
Mean absolute error (MAE)MAE=1n∑i=1n[|Mz(x,y)−Iz(x,y)|]
Root mean square error (RMSE)RMSE=1n∑i=1n[Mz(x,y)−Iz(x,y)]2

Note: n is the number of test points.

To present the relevance of DTM enhancement for quantification of hydrological features of the floodplain, we compared flood extents, depths of water within the flooded plain, and volumes of flood computed with the uncorrected and corrected DTMs applied to the analysis. Analyses were based on the intersection of grid datasets containing elevations of water table in various flood scenarios. (We used data of flooding in 1%, 5%, 10%, 20%, and 50% recurrence intervals computed by Grygoruk et al.49) Differences in assessed values of these parameters are expected to present the relevance of DTM correction for hydrological estimation of floodplains and their role in shaping ecological features of riparian wetlands.

3.

Results and Discussion

3.1.

Accuracy Assessment Results

The accuracy assessment was performed on two developed DTMs (Figs. 5 and 6). A comparison of the accuracy of the analyzed DTMs shows that the best-fit between the ground-measured elevations and DTM-derived elevations was obtained in the case of the processed ALS DTM with “vegetation effect correction” (Table 4).

Fig. 5

The ALS DTM without “vegetation effect” correction: (a) DTM and (b) hillshade map.

JARS_10_3_036013_f005.png

Fig. 6

The ALS DTM with “vegetation effect” correction: (a) DTM and (B) hillshade map.

JARS_10_3_036013_f006.png

Table 4

ME, MAE, RMSE, minimum (Min), maximum (Max), and SD of discrepancies between ground elevations of analyzed DTMs versus field-measured ground elevations. Values are given in meters.

DTMMEMAERMSEMinMaxSD
DTM with additional correction0.130.300.40−1.101.040.38
DTM without additional correction0.600.101.15−1.207.230.98

Table 5

Calculated flood areas results for each DTM. Values are given in hectares.

DTMFlood recurrence time
1%5%10%20%50%
DTM with additional correction50764956487247544525
DTM without additional correction36923535343533023058

Considering the fact that the variance of n independent variables having the same distribution equals 1/n of the variance of each individual variable, the summaric influence of GPS and ALS accuracy on calculated correction factor reach the value of [(0.052+0.152)/n]0.5, and—in our case—never exceeds 0.03 m (Table 2). This value is significantly lower than the assessed and developed correction factors for particular vegetation types analyzed (Table 2).

The ME of the ALS DTM versus the RTK GPS field measurements of ground elevation reached <0.13  m, which we consider to be acceptable for application in hydrological analyses, such as flood depths and extents calculation. Although the maximum error of this DTM reached 1.1 m, the standard deviation (SD) of all the values analyzed allows the conclusion that although errors higher than 0.5 m have occurred [32% of the data has an error of (0.12±0.4  m)], their frequency of occurrence was considerably lower than in the case of the uncorrected DEM (Fig. 7). The errors higher than 5 m were associated with forest class.

Fig. 7

Histograms of elevation errors for ALS-derived DTMs [(a) without and (b) with additional correction] in comparison to ground-measured data from RTK GPS.

JARS_10_3_036013_f007.png

The application of GEOBIA we present in our study remains the first ever use of this methodology for improving DTM on the basis of vegetation classes distinguished. Reaching the classification accuracy at the level of 80%, which we consider more than satisfactory and very promising as the research task to be further developed, our study opens the field for the improvement of ALS-based DTMs of densely vegetated wetlands. Among the main issues we faced in the presented approach, there was unsatisfactory accuracy of terrain elevation measurements done with RTK GPS in patches of high (>2.5  m) bushes and dense and high (>3.0  m) reeds, which were used to calculate the correction factors related to particular types of the vegetation. The next step in the development of the presented herein methodology of DTM improvement would be to make accurate measurements of terrain elevations within extents of those vegetation types, preferably done with the use of standard optical geodetic devices such as the total station. Such measurements, being much more accurate than the ones done with RTK GPS, would allow gathering high-quality elevation data within patches of high and dense vegetation. Although low accessibility to the study area and harsh field conditions did not allow such measurements to be undertaken, we foresee that the application of optical geodetic measurements would allow reaching even higher accordance of performed DTMs comparing to the field-collected data.

3.2.

Flood Mapping Results

To present how different DTMs applied to flood-mapping procedures affected the final results of the flood area and floodplain volume calculations, we analyzed floods of 1%, 5%, 10%, 20%, and 50% recurrence (Table 5). Water levels were modeled using a one-dimensional hydrodynamic model.50

Being aware that microtopography and the presence of dense vegetation may significantly modify flood distribution and should therefore be determined with high accuracy, the ALS-based DTM with additional correction was expected to give the most reliable results of flood analysis. The result shows that the application of coarse quality elevation data in flood mapping can lead to a significant underestimation of flood extent and flood depth. Such a fact can result in inappropriate decision making when hydrological characteristics play important roles in ecosystem management (Fig. 8). As ecological analyses frequently refer to flood depths, we evaluated maximum and mean depths of the flood along with SDs (Table 6).

Fig. 8

The flood depth results for 50% recurrence time for ALS DTM with (a) additional correction and (b) without such correction.

JARS_10_3_036013_f008.png

Table 6

Calculated water depth results for each DTM. Values are given in meters.

DTMFlood depthFlood recurrence time
1%5%10%20%50%
DTM with additional correctionMax2.422.282.202.101.94
Mean0.890.780.730.660.56
SD0.380.360.340.330.30
DTM without additional correctionMax2.182.031.941.671.59
Mean0.790.700.650.590.50
SD0.390.360.350.330.30

Maximum flood depth reached 2.18  m (ALS DTM without additional correction), which is 0.24 m lower than the value of maximum depth obtained in calculations using ALS with additional processing methods. Again, the results of flood depth calculations obtained with a course DTM applied to lowland floodplains do not allow comprehensive conclusions of floodplain hydrology to be derived.

Knowing that the quality of the DTMs affected flood extent and flood depth calculations, we also compared floodplain water retention capacities (Table 7). Respectively, the calculated flood volumes were higher in average by 13.04  M/m3 for the DTM with correction. The flood volume for 1% flood recurrence time was higher by 15.93  M/m3 and for 50% flood recurrence time by 9.95  M/m3. Our analyses confirmed that in addition to the importance of vegetation influence on computed hydrological features of floodplain wetlands (volumes of vegetation versus volumes of flood),51 it is highly important to consider vegetation-type-dependent quality of DTM to obtain reliable results of topographic and ecohydrological analyses of these ecosystems.

Table 7

Calculated flood volumes for each DTM. Values are given in M/m3.

DTMFlood recurrence time
1%5%10%20%50%
DTM with additional correction44.9638.8235.3931.3425.13
DTM without additional correction29.0324.6222.2119.3915.18

4.

Conclusions

The main goal of this research was to efficiently process an ALS point cloud dataset to obtain an accurate DTM of a densely vegetated floodplain, without losing relevant terrain information that may affect future hydrological analyses wherein a DTM is to be applied. The combination of different types of data, including topography and spatial patterns of wetland vegetation, provided a basis to develop an accurate DTM. Methodology for the development of DTMs based on ALS and including GEOBIA classification and plant height correction coefficients was proven to be an effective tool for high-quality DTM development. Comparison of the DTM produced using the standard procedures with the DTM developed with the proposed methodology showed that RMSE values were up to three times lower for the vegetation-height-corrected DTM. Application of GEOBIA for vegetation mapping, which was a novel approach to DTM improvement, achieved an overall classification accuracy of 80%. Correction factors developed for various types of the vegetation reached values from 0.08 up to 0.92 m and were dependent on the vegetation type. The methodology of removing dense vegetation from the ALS point cloud, then interpolating terrain elevations for these spatial gaps in data considerably increased the accuracy of the DTM. Although the values of vegetation height correction coefficients seem universal for the plant communities analyzed, we strongly recommend to derive appropriate values of these coefficients in particular research sites other than the Biebrza Valley. Flood volume assessment discrepancies between the two DTMs applied reached values as high as 66%. The proposed method, although, was proven to be capable of providing a good quality DTM to be used in hydrological analyses of the temperate floodplain wetlands, is more resource- and time-consuming than the standard one, mainly due to the extensive field measurements that have to be done in order to obtain the accurate field elevation data. In spite of this disadvantage, we recommend this methodology for all the cases of riparian wetlands, where accurate representation of floodplain’s topography may play the crucial role in the appropriate management of wetlands.

Acknowledgments

This study was part of the research task realized in the framework of the Polish-Norwegian Research Program (Norway Grants): “Integrated wetland flow model: a tool for management adaptation in protected wetlands—WETFLOD” (Nor/209804/42/2013). The authors wish to acknowledge valuable remarks of two anonymous reviewers who greatly helped to improve the paper.

References

1. 

P. Tarolli, “High-resolution topography for understanding Earth surface processes: opportunities and challenges,” Geomorphology 216, 295–312 (2014).AWREDI0309-1708http://dx.doi.org/10.1016/j.geomorph.2014.03.008Google Scholar

2. 

D. Hughes, “Floodplain inundation-processes and relationships with channel discharge,” Earth Surf. Processes Landforms 5, 297–304 (1980).ESPLDB0197-9337http://dx.doi.org/10.1002/(ISSN)1931-8065Google Scholar

3. 

J. Lewin and D. Hughes, “Welsh floodplain studies: II. Application of a qualitative inundation model,” J. Hydrol. 46, 35–49 (1980).JHYDA70022-1694http://dx.doi.org/10.1016/0022-1694(80)90034-7Google Scholar

4. 

P. D. Bates and A. P. J. De Roo, “A simple raster-based model for flood inundation simulation,” J. Hydrol. 236, 54–77 (2000).JHYDA70022-1694http://dx.doi.org/10.1016/S0022-1694(00)00278-XGoogle Scholar

5. 

A. Nicholas and C. Mitchell, “Numerical simulation of overbank processes in topographically complex floodplain environments,” Hydrol. Processes 17, 727–746 (2003).HYPRE31099-1085http://dx.doi.org/10.1002/(ISSN)1099-1085Google Scholar

6. 

M. S. Horritt and P. D. Bates, “Evaluation of 1D and 2D numerical models for predicting river flood inundation,” J. Hydrol. 268, 87–99 (2002).JHYDA70022-1694http://dx.doi.org/10.1016/S0022-1694(02)00121-XGoogle Scholar

7. 

P. Banaszuk, A. K. Kamocki and R. Zarzecki, “Mowing with invasive machinery can affect chemistry and trophic state of rheophilous mire,” Ecol. Eng. 86, 31–38 (2016).http://dx.doi.org/10.1016/j.ecoleng.2015.10.005Google Scholar

8. 

W. Kotowski, E. Jabłońska and E. Bartoszuk, “Conservation management in fens: do large tracked mowers impact functional plant diversity?” Biol. Conserv. 167, 292–297 (2013).BICOBK0006-3207http://dx.doi.org/10.1016/j.biocon.2013.08.021Google Scholar

9. 

W. J. Junk, P. B. Bayley and R. E. Sparks, “The flood pulse concept in river-floodplain systems,” Can. Spec. Publ. Fish. Aquat. Sci. 106, 110–127 (1989).CSPSDA0706-6481Google Scholar

10. 

M. C. Thoms, “Floodplain-river ecosystems: lateral connections and the implications of human interference,” Geomorphology 56, 335–349 (2003).http://dx.doi.org/10.1016/S0169-555X(03)00160-0Google Scholar

11. 

S. E. Bunn et al., “Flow variability in dryland rivers: boom, bust and the bits in between,” River Res. Appl. 22, 179–186 (2006).RRAIAQ1535-1459http://dx.doi.org/10.1002/(ISSN)1535-1467Google Scholar

12. 

K. A. Tockner, M. S. Lorang and J. A. Stanford, “River flood plains are model ecosystems to test general hydrogeomorphic and ecological concepts,” River Res. Appl. 26, 76–86 (2010).RRAIAQ1535-1459http://dx.doi.org/10.1002/rra.1328Google Scholar

13. 

M. W. Scown, M. C. Thoms and N. R. De Jager, “An index of floodplain surface complexity,” Hydrol. Earth Syst. Sci. 20, 431–441 (2016).http://dx.doi.org/10.5194/hess-20-431-2016Google Scholar

14. 

C. C. Sampson et al., “Use of terrestrial laser data to drive decimetric resolution urban inundation models,” Adv. Water Resour. 41, 1–17 (2012).http://dx.doi.org/10.1016/j.advwatres.2012.02.010Google Scholar

15. 

G. Wierzbicki et al., “Using VHR multispectral remote sensing and LIDAR data to determine the geomorphological effects of overbank flow on a floodplain (the Vistula River, Poland),” Geomorphology 183, 73–81 (2013).http://dx.doi.org/10.1016/j.geomorph.2012.06.020Google Scholar

16. 

D. Shen et al., “Integration of 2-D hydraulic model and high-resolution Lidar-derived DEM for floodplain flow modelling,” Hydrol. Earth Syst. Sci. 19, 3605–3616 (2015).http://dx.doi.org/10.5194/hess-19-3605-2015Google Scholar

17. 

S. Szporak-Wasilewska et al., “Processing of airborne laser scanning data to generate accurate DTM for floodplain wetland,” Proc. SPIE 9637, 963720 (2015).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.2195223Google Scholar

18. 

X. Shi et al., “A comparison of Lidar-based DEMs and USGS-sourced DEMs in terrain analysis for knowledge-based digital soil mapping,” Geoderma 170, 217–226 (2012).GEDMAB0016-7061http://dx.doi.org/10.1016/j.geoderma.2011.11.020Google Scholar

19. 

J. Vaze, J. Teng and G. Spencer, “Impact of DEM accuracy and resolution on topographic indices,” Environ. Modell. Softw. 25, 1086–1098 (2010).http://dx.doi.org/10.1016/j.envsoft.2010.03.014Google Scholar

20. 

B. Gorte, N. Pfeifer and S. O. Elberink, “Height texture of low vegetation in airborne laser scanner data and its potential for DTM correction,” Int. Arch. Photogramm. Remote Sens. Spatial Inf. Sci. 36, 150–155 (2005).1682-1750Google Scholar

21. 

M. E. Hodgson et al., “An evaluation of lidar- and IFSAR-derived digital elevation models in leaf-on conditions with USGS level 1 and level 2 DEMs,” Remote Sens. Environ. 84, 295–308 (2003).http://dx.doi.org/10.1016/S0034-4257(02)00114-1Google Scholar

22. 

E. Ahokas, H. Kaartinen and J. Hyyppä, “A quality assessment of airborne laser scanner data,” Int. Arch. Photogramm. Remote Sens. Spatial Inf. Sci. 34(Part 3), 23 (2003).1682-1750Google Scholar

23. 

A. Bollweg and R. de Lange, “Wat ruist er door het struikgewas… Welke maaiveldhoogte wordt met laseraltimetrie gemeten bij aanwezigheid van vegetation?” Technical Report AGI-GAR-2003-22, AGI, RWS, The Netherlands (2003). (in Dutch)Google Scholar

24. 

J. Su and E. Bork, “Influence of vegetation, slope, and Lidar sampling angle on DEM accuracy,” Photogramm. Eng. Remote Sens. 72, 1265–1274 (2006).http://dx.doi.org/10.14358/PERS.72.11.1265Google Scholar

25. 

C. M. Hladik and M. Alber, “Accuracy assessment and correction of a LIDAR-derived salt marsh digital elevation model,” Remote Sens. Environ. 121, 224–235 (2012).http://dx.doi.org/10.1016/j.rse.2012.01.018Google Scholar

26. 

K. Będkowski and K. Stereńczak, “Porównanie numerycznych modeli terenu obszarów leśnych generowanych z wykorzystaniem danych skaningu laserowego (LIDAR) uzyskanych w okresie wiosennym i letnim,” Roczniki Geomatyki 7(43), 11–20 (2010). (in Polish)Google Scholar

27. 

G. J. Hay, G. Castilla, “Geographic object-based image analysis (GEOBIA): a new name for a new discipline,” in Object-Based Image Analysis—Spatial Concepts for Knowledge-Driven Remote Sensing Applications, , T. Blaschke, S. LangG. J. Hay, Eds., pp. 75–89, Springer-Verlag, Berlin Heidelberg (2008).Google Scholar

28. 

M. J. Wassen et al., “Eco-hydrological functioning of the Biebrza wetlands: lessons for the conservation and restoration of deteriorated wetlands,” in Wetlands: Functioning. Biodiversity Conservation and Restoration, and R. Bobbink et al., Eds., pp. 285–310, Springer-Verlag, Berlin, Germany (2006).Google Scholar

29. 

F. Keizer et al., “A new look at the flood pulse concept: the (ir)relevance of the moving littoral in temperate zone rivers,” Ecol. Eng. 64, 85–99 (2014).http://dx.doi.org/10.1016/j.ecoleng.2013.12.031Google Scholar

30. 

M. Grygoruk et al., “How much for water? Economic assessment and mapping of floodplain water storage as a catchment-scale ecosystem service of wetlands,” Water 5, 1760–1779 (2013).http://dx.doi.org/10.3390/w5041760Google Scholar

31. 

Z. Kurczyński “Lotnicze skanowanie laserowe–podstawy teoretyczne,” in Podręcznik dla Uczestników Szkoleń z Wykorzystania Produktów LiDAR, pp. 31–532, Main Council of Geodesy and Cartography, Warszawa (2015). (in Polish).Google Scholar

32. 

Z. Kurczyński, E. Stojek and U. Cisło-Lesicka, “Zadania GUGiK realizowane w ramach projektu ISOK,” in Podręcznik dla Uczestników Szkoleń z Wykorzystania Produktów LiDAR, pp. 31–532, Main Council of Geodesy and Cartography, Warszawa (2015). (in Polish).Google Scholar

33. 

A. Räsänen et al., “What makes segmentation good? A case study in boreal forest habitat mapping,” Int. J. Remote Sens. 34, 8603–8627 (2013).IJSEDK0143-1161http://dx.doi.org/10.1080/01431161.2013.845318Google Scholar

34. 

G. J. Hay et al., “An automated object-based approach for the multi-scale image segmentation of forest scenes,” Int. J. Appl. Earth Obs. Geoinf. 7, 339–359 (2005).http://dx.doi.org/10.1016/j.jag.2005.06.005Google Scholar

35. 

Y. Gao and J. F. Mas, “A comparison of the performance of pixel based and object based classifications over images with various spatial resolutions,” Online J. Earth Sci. 2, 27–35 (2008).Google Scholar

36. 

T. Blaschke, “Object based image analysis for remote sensing,” ISPRS J. Photogramm. Remote Sens. 65, 2–16 (2010).IRSEE90924-2716http://dx.doi.org/10.1016/j.isprsjprs.2009.06.004Google Scholar

37. 

M. Szostak, P. Wężyk and P. Tompalski, “Aerial orthophoto and airborne laser scanning as monitoring tools for land cover dynamics: a case study from the Milicz Forest District (Poland),” Pure Appl. Geophys. 171, 857–866 (2014).PAGYAV0033-4553http://dx.doi.org/10.1007/s00024-013-0668-8Google Scholar

38. 

J. W. Rouse et al., “Monitoring the vernal advancement and retrogradation (greenwave effect) of natural vegetation,” NASA/GSFC Type III Final Report, Greenbelt (1974).Google Scholar

39. 

G. Rondeaux, M. Steven and F. Baret, “Optimization of soil-adjusted vegetation indices,” Remote Sens. Environ. 55, 95–107 (1996).http://dx.doi.org/10.1016/0034-4257(95)00186-7Google Scholar

40. 

A. Huete, “A soil adjusted vegetation index (SAVI),” Remote Sens. Environ. 23, 213 (1988).http://dx.doi.org/10.1016/0034-4257(87)90038-1Google Scholar

41. 

J. Chen, “Evaluation of vegetation indices and modified simple ratio for boreal applications,” Can. J. Remote Sens. 22, 229–242 (1996).CJRSDP0703-8992http://dx.doi.org/10.1080/07038992.1996.10855178Google Scholar

42. 

J. Imanishi et al., “The independent detection of drought stress and leaf density using hyperspectral resolution data,” Landscape Ecol. Eng. 3(1), 55–65 (2007).http://dx.doi.org/10.1007/s11355-007-0019-yGoogle Scholar

43. 

P. Baldi et al., “Validation and comparison of different techniques for the derivation of digital elevation models and volcanic monitoring (Vulcano Island, Italy),” Int. J. Remote Sens. 23, 4783–4800 (2002).IJSEDK0143-1161http://dx.doi.org/10.1080/01431160110115861Google Scholar

44. 

T. L. Webster et al., “Flood-risk mapping for storm-surge events and sea-level rise using LiDAR for southeast New Brunswick,” Can. J. Remote Sens. 32, 194–211 (2006).CJRSDP0703-8992http://dx.doi.org/10.5589/m06-016Google Scholar

45. 

J. Oksanen and T. Sarjakoski, “Error propagation analysis of DEM-based drainage basin delineation,” Int. J. Remote Sens. 26, 3085–3102 (2005).IJSEDK0143-1161http://dx.doi.org/10.1080/01431160500057947Google Scholar

46. 

J. Chormański, D. Świątek, T. Okruszko, “Remote sensing limitation in flood modeling verification in wetlands,” in Model Application for Wetlands Hydrology and Hydraulics, , J. Kubrak, T. Okruszko and S. Ignar Eds., pp. 51–72, Warsaw Agricultural University Press, Warsaw (1996).Google Scholar

47. 

M. Flood, “ASPRS guidelines—vertical accuracy reporting for LiDAR data,” in American Society of Photogrammetry and Remote Sensing LiDAR Committee (PAD) (2004).Google Scholar

48. 

J. Höhle and M. Putuckova, “The EuroSDR test: checking and improving of digital terrain models,” in EuroSDR European Spatial Data Research, No. 51, pp. 9–141, EuroSDR Publications, Frankfurt, Germany (2006).Google Scholar

49. 

M. Grygoruk et al., “Analysis of past-present hydrological phenomena of the Biebrza Valley,” HABIT-CHANGE Project Report, Biebrza National Park (2011).Google Scholar

50. 

D. Mirosław-Świątek, J. Kubrak and J. Chormański, “Steady 1-D water surface model of natural rivers with vegetated floodplain: an application to the Lower Biebrza,” in Proc. Int. Conf. on Fluvial Hydraulics River Flow, Vol. 1, pp. 545–553 (2006).Google Scholar

51. 

D. Mirosław-Świątek, S. Szporak-Wasilewska and M. Grygoruk, “Assessing floodplain porosity for accurate quantification of water retention capacity of near-natural riparian ecosystems—a case study of the Lower Biebrza Basin, Poland,” Ecol. Eng. 92, 181–189 (2016).http://dx.doi.org/10.1016/j.ecoleng.2016.03.031Google Scholar

Biography

Dorota Mirosław-Świątek received her PhD in environmental engineering from Institute of Meteorology and Water Management, Warsaw, Poland, in 1998. She is an assistant professor in the Department of Hydraulic Engineering of Warsaw University of Life Sciences (WULS-SGGW). Her scientific expertise covers surface water and groundwater flow modelling, floodplain hydraulics, GIS and risk assessment in environmental systems.

Sylwia Szporak-Wasilewska received her PhD in agricultural sciences from WULS-SGGW, Warsaw, Poland, in 2012. Since then she has worked as an assistant professor in the Water Centre Laboratory, Faculty of Civil and Environmental Engineering (WULS-SGGW). Her research interests include remote sensing, geoinformation systems, ecohydrology, and environmental protection.

Robert Michałowski is an engineer technician at SGGW-WULS, Warsaw, Poland. His expertise covers field measurements (hydrology, hydraulics, geodesy), data processing, and database management. He also develops tools in C/C++, C#, Fortran and various database SQL languages, used in data extraction and analysis, especially with respect to the GIS applications. He is a certified database developer.

Ignacy Kardel received his PhD in agricultural sciences in the field of environmental management from WULS-SGGW, Warsaw, Poland, in 2001. Since then he has worked as an assistant professor in the Department of Hydraulic Engineering (WULS-SGGW). His research interests include hydrology and water management, surface water quality modelling, environmental engineering, GIS, databases, and decision support systems.

Mateusz Grygoruk received his PhD in engineering from Vrije Universiteit Brussel, Belgium, and a PhD in agricultural sciences from WULS-SGGW, Warsaw, Poland, in 2013. Since then he has worked as an assistant professor in the Department of Hydraulic Engineering (WULS-SGGW). He also worked as a GIS specialist in Biebrza National Park. His research interests include ecohydrology, environmental management, and GIS.

Dorota Mirosław-Świątek, Sylwia Szporak-Wasilewska, Robert Michałowski, Ignacy Kardel, Mateusz Grygoruk, "Developing an algorithm for enhancement of a digital terrain model for a densely vegetated floodplain wetland," Journal of Applied Remote Sensing 10(3), 036013 (12 August 2016). http://dx.doi.org/10.1117/1.JRS.10.036013
JOURNAL ARTICLE
16 PAGES


SHARE
KEYWORDS
Vegetation

Algorithm development

Floods

Global Positioning System

Data modeling

Clouds

Satellites

Back to Top