Next Article in Journal
Laser-based Thickness Control in a Double-Side Polishing System for Silicon Wafers
Previous Article in Journal
Event-Based Feature Extraction Using Adaptive Selection Thresholds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rectification of Bowl-Shape Deformation of Tidal Flat DEM derived from UAV Imaging

1
Department of Civil Engineering, Sunchon National University, 255 Jungangro, Suncheon, Jellanamdo 57922, Korea
2
Department of Civil Engineering, Chonnam National University, 77 Yongbongro, Bukgu, Gwangju 61186, Korea
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(6), 1602; https://doi.org/10.3390/s20061602
Submission received: 15 January 2020 / Revised: 6 March 2020 / Accepted: 11 March 2020 / Published: 13 March 2020
(This article belongs to the Section Remote Sensors)

Abstract

:
It is necessary to periodically obtain topographic maps of the geographical and environmental characteristics of tidal flats to systemically manage and monitor them. Accurate digital elevation models (DEMs) of the tidal flats are produced while using ground control points (GCPs); however, it is both complicated and difficult to conduct GPS surveys and readings of image coordinates that correspond to these because tidal flat areas are not easy to access. The position and distribution of GCPs affect DEMs, because the entire working area cannot be covered during a survey. In this study, a least-squares height-difference (LHD) DEM matching method with a polynomial model is proposed to increase the number of DEM grids while using a presecured precise DEM to rectify the distortion and bowl effect produced by unmanned aerial vehicle (UAV) images. The most appropriate result was obtained when the translation parameters were quadratic curve polynomials with an increasing number of grids and the rotation parameters were constant. The experimental results indicated that the proposed method reduced the distortion and eliminated the error caused by the bowl effect while only using a reference DEM.

1. Introduction

Coastal areas are vulnerable regions that are easily affected by natural and artificial environmental changes. There is growing interest in these areas in a number of associated fields, as various forms of land (sediment) have either emerged or disappeared, owing to these changes [1,2]; a unique ecological environment among these is the intertidal flat. The intertidal flat is submerged in water during high tide and its surface is fully exposed during low tide. It also possesses significant features in terms of habitat for shellfish, water purification, maintenance of diverse species, flood control, and recreational and scenic resources [3]. The Korean intertidal flat is one of the largest tidal flats in the world. It is also an internationally protected shelter for endangered migratory birds and is considered to be one of the prosperous ecosystems inhabited by various species living in clusters [4].
It is necessary to periodically obtain topographic maps that understand the geographical and environmental characteristics of tidal flats to systemically manage and monitor them. The bottom topography of an intertidal flat in a time series manner can be used, particularly to estimate the annual or seasonal variations in the sediment budget and coastal changes [5,6,7]. Very few investigations have been carried out so far despite the importance of intertidal flats, as access to such areas is restricted. An intertidal flat consists of inflow of seawater, soil with no solidification, and deposits of mud and silt. Therefore, non-intrusive and temporally regular measurements over intertidal flats are essential.
Until now, the topography of tidal flats was obtained by using either airborne light detection and ranging (LiDAR) or echo sounders [8,9]. In 1999, Irish et al. introduced a scanning hydrographic operational airborne lidar survey (SHOALS) that could simultaneously measure the depth of water and terrain of the coastal area [10]. Measurements using the airborne LiDAR system incurred high costs and involved a number of processing tasks. Moreover, the system did not work adequately in the case of tidal flats, as they contain remnants of water and puddles of sea water [5,11]. The use of echo sounders was also not suitable, as the ship carrying the hull-mounted sonar could not conduct the survey in conditions where the water-level was very shallow, such as in intertidal flats. Owing to these limitations, the waterline method was widely utilized to generate the topography of tidal flats [12,13]. It is known that the waterlines extracted from optical and SAR satellite images can be stacked to construct contour lines corresponding to the height of the tide at the time of passage of the satellite [6]. However, this method requires a large amount of satellite data that are acquired in a short period of time and at different tidal heights, as well as very accurate records of tidal gauge data. Moreover, in 2007, Kim et al. reported that the waterlines extracted from SAR images may be different from the actual waterlines, depending on the radar frequency used [14]. Consequently, the above-mentioned methods were not suitable for generating the bottom topography of intertidal flats.
In recent times, a number of studies were conducted to assess the possibility of generation of the topography of tidal flats using unmanned aerial vehicle (UAV) images and photogrammetric software, such as Agisoft PhotoScan or Pix4D. According to these studies, on-site GCPs were required for the accurate development of digital elevation model (DEM) [2,7,15]. However, a survey of the GCPs was difficult, as it was not easy to approach the tidal flats in the Southwest coast of the Korean peninsula, where the proportion of silt and clay content were high (http://www.ecosea.go.kr) [16].
Without GCPs, the DEM generated from the UAV images and photogrammetric software probably have a nonlinear distortion that overestimated the elevation of the terrain. The systematic overvaluation of the terrain increased as the distance from the center increased, due to the ‘bowl or doming effect’ [17,18], which affected the external part of the area that is covered by the GCPs [19]. Recently, UAVs with an integrated on-board real-time kinematics (RTK) have been used for mapping terrain. In an area of about 550 × 330 m2 containing buildings, roads, and meadows, the average DSM accuracy had a ground sampling distance that was obtained while using SenseFly eBee-RTK was approximately 3.7 without the GCPs [20]. However, an assessment of the accuracy of the DEM shows an average root mean square error (RMSE) of approximately 2 m without GCPs when a DJI phantom 4 RTK is flown on a north-to-south trajectory over a 2 km coastal site [21]. The base station must be placed even if no GCP has to be surveyed, especially in case of inaccessible areas.
This study proposes an least-squares height-difference (LHD) DEM matching technique with a polynomial model to correct the distortion in the DEM produced by the UAV, while using the existing reference DEM such as LiDAR or pre-compensated dense DEM instead of the GCPs. The LHD technique is based on the height difference in the same-plane positions between the two DEMs for the calculation of seven three-dimensional transformation parameters (translations, rotation angles, and scale) [22,23,24]. However, the distortion due to the bowl effect could not be removed, even if the DEM produced by the UAV images can use the reference DEM to correct the these parameters. Therefore, an LHD matching technique with a polynomial model was developed for refinement of the distorted DEM. This approach determines the polynomial coefficients that depend on the DEM grid lines for transformation parameters.
The remainder of this paper is organized, as follows: The fundamentals of the proposed LHD matching are introduced in Section 2, with an explanation of the various input parameters. Section 3 details the fields of experiment, data processing, and instrumentation. Section 3 also discusses the experimental results, where the analysis focuses on the enhancement of accuracy of the corrected DEM. Finally, Section 4 presents the conclusions.

2. Proposed Approach of LHD Matching

A UAV DEM built without GCPs usually contains a deformation, owing to the bowl effect. The bowl effect is probably affected by location errors in the tie-point that is used by the software for self-calibration and aerial triangulation. Therefore, a non-metric camera with inaccurate interior orientation parameters (IOPs, i.e., focal length, principal point displacement, and lens distortion) is mounted on the UAV and the exterior orientation parameters (EOPs, i.e., position and rotation of the principal point) are determined by errors that are propagated by the bundle block adjustment process between overlapping images when a camera model with initially inaccurate values is used [17,18,19].
The purpose of our study was to rectify the DEM that is produced by the overlapping UAV images and photogrammetric software using LHD matching, because it is difficult to collect GCPs via field surveys in an inaccessible area, such as a tidal flat. We used only the ground coordinates of the pre-established reference DEM, in particular, to apply the LHD. The heights of the DEM were extracted by determining same-plane positions within a grid interval in both the reference and produced DEMs, as shown in Figure 1.
Without the GCPs, a bowl-shaped deformation cannot be rectified, even by applying the LHD matching method while using the seven transform parameters and reference DEM (Figure 1). We proposed a new LHD matching technique that was applied by using polynomial models of the transformation parameters as functions of the number of grid spaces between the two DEMs to overcome this problem (Equation (1)).
Figure 2 shows a flowchart of the steps used in the study. LHD matching is based on the three-dimensional (3-D) similarity equation (Equation (1)) that specifies the relationship between the reference and UAV DEMs (translations, rotation angles, and scale). The polynomial orders were selected after testing these functions.
A seven-parameter 3-D similarity transformation was used to express the geometric relationship between the reference DEM and UAV DEM points:
[ X R E F Y R E F Z R E F ] = s R [ X U A V Y U A V Z U A V ] + [ X T Y T Z T ]
where X R E F , Y R E F , and Z R E F are the 3-D coordinates of the reference DEM; and, X U A V , Y U A V , and Z U A V are the 3-D coordinates of the UAV DEM. X R E F , Y R E F and X U A V , Y U A V are the same-plane coordinates, whereas Z R E F and Z U A V are the different heights.   X T , Y T , and Z T are the translations, s   i s a scale factor, and R =   [ r 11 r 12 r 13 r 21 r 22 r 23 r 31 r 32 r 33 ] is an orthogonal rotation matrix that contains the rotational angles (ω, φ, κ) with respect to the X, Y, and Z axes, respectively. X T = X 0 + X 1 · l + X 2 · l 2 + X 3 · l 3 , Y T = Y 0 + Y 1 · l + Y 2 · l 2 + Y 3 · l 3 , Z T = Z 0 + Z 1 · l + Z 2 · l 2 + Z 3 · l 3 , ω = ω 0 + ω 1 · l + ω 2 · l 2 + ω 3 · l 3 , φ = φ 0 + φ 1 · l + φ 2 · l 2 + φ 3 · l 3 , κ = κ 0 + κ 1 · l + κ 2 · l 2 + κ 3 · l 3 . l   is the line as grid space of the DEM, the direction of which is the same as the flight path. X 0 , X 1 , X 2 , X 3 , Y 0 , Y 1 , Y 2 , Y 3 , Z 0 , Z 1 , Z 2 , Z 3 , and ω 0 , ω 1 , ω 2 , ω 3 , φ 0 , φ 1 , φ 2 , φ 3 , κ 0 , κ 1 , κ 2 , κ 3 are the polynomial coefficients of the grid lines.
Equation (1) can be linearized by Taylor series expansion, as shown in Equation (2), based on the theory of Rosenholm and Torlegård (1988):
X R E F + Δ X = s o ( r 11 o X U A V + r 12 o Y U A V + r 13 o Z U A V ) + X T o Δ κ Y U A V + Δ φ Z U A V + Δ s X U A V + Δ X T Y R E F + Δ Y = s o ( r 21 o X U A V + r 22 o Y U A V + r 23 o Z U A V ) + Y T o + Δ κ X U A V Δ ω Z U A V + Δ s Y U A V + Δ Y T Z R E F + Δ Z = s o ( r 31 o X U A V + r 32 o Y U A V + r 33 o Z U A V ) + Z T o Δ φ X U A V + Δ ω Y U A V + Δ s Z U A V + Δ Z T
where the index o indicates an initial value. Δ X T = Δ X 0 + Δ X 1 · l + Δ X 2 · l 2 + Δ X 3 · l 3 , Δ Y T = Δ Y 0 + Δ Y 1 · l + Δ Y 2 · l 2 + Δ Y 3 · l 3 , Δ Z T = Δ Z 0 + Δ Z 1 · l + Δ Z 2 · l 2 + Δ Z 3 · l 3 , Δ ω = Δ ω 0 + Δ ω 1 · l + Δ ω 2 · l 2 + Δ ω 3 · l 3 , Δ φ = Δ φ 0 + Δ φ 1 · l + Δ φ 2 · l 2 + Δ φ 3 · l 3 , Δ κ = Δ κ 0 + Δ κ 1 · l + Δ κ 2 · l 2 + Δ κ 3 · l 3 .
Equation (2) can be rewritten in terms of the Δ X   and Δ Y   portions, as in Equation (3):
Δ X = Δ X T Δ κ Y U A V + Δ φ Z U A V + Δ s X U A V
Δ Y = Δ Y T + Δ κ X U A V Δ ω Z U A V + Δ s Y U A V
If the Z-axis of the reference DEM is the same as that of the UAV DEM, then Z R E F = f ( X U A V , Y U A V ) . The equation can be linearized as in Equation (4).
Δ Z = f X U A V Δ X + f Y U A V Δ Y
Δ X     and   Δ Y from Equation (3) are substituted in Equation (4), which is then expressed, as follows:
Δ Z = f X U A V ( Δ X T Δ κ Y U A V + Δ φ Z U A V + Δ s X U A V ) + f Y U A V ( Δ Y T + Δ κ X U A V Δ ω Z U A V + Δ s Y U A V )
Δ Z of Equation (5) is substituted into the Z term in Equation (2), and the transformation parameters are finally determined by Equation (6):
Z R E F i Z R E F i o = f i X U A V i Δ X T f i Y U A V i   Δ Y T + Δ Z T
+ [ f i X U A V i X U A V i f i Y U A V i   Y U A V i + Z U A V i ] Δ s
+ (   Y U A V i + f i Y U A V i   Z U A V i ) Δ ω
+ ( X U A V i f i X U A V i   Z U A V i ) Δ φ
+ (   f i X U A V i Y U A V i f i Y U A V i   X U A V i ) Δ κ
where i = 1,2, …, n (n = number of points), and Z R E F i o = s o ( r 31 o X U A V i + r 32 o   Y U A V i + r 33 o   Z U A V i ) + Z T o , which is the initial height calculated while using the initial polynomial parameters. In the Z term in Equation (1), f i X U A V i     and f i Y U A V i   are the slopes with respect to the X and Y directions of the two grid spaces in the produced DEM. If the value of the slope is close to zero, the accuracy of the parameters is low, as shown in Equation (6).
Δ   is computed by the least squares method and iteration ( Z T 1 = Z T o + Δ Z T ; Δ Z T = Δ Z 0 + Δ Z 1 · l + Δ Z 2 · l 2 + Δ Z 3 · l 3 ) to calculate the transformation parameters. The transformation parameters are then determined until the Δ terms are nearly zero. LHD matching uses the height difference at the same-plane locations between the two DEMs. It is important that the height difference ( d i ) of any point ( i ) should not be largely different from those of other points. Therefore, we first determined the same-plane positions between the two DEMs within a grid interval and then calculated the height differences at those points before applying the LHD matching. Only the point where d i is within 2 σ (σ is the standard deviation) of d m   can participate in the LHD matching operation to limit the gross error point:
d m = d 1 + + d n n
where d m   is the average height difference between all of the points in the search.

3. Experiment and Results

The test site chosen to apply the above-mentioned method was a tidal flat area in Hampyung Bay, which is located on the west coast of the Korean Peninsula (Figure 3). The site is an embayment situated in a deep inland from the west coast and it is a unique geographic feature that is formed without large rivers. The UAV images were acquired while using a DJI Phantom 4, as shown in Table 1, at 2:20 pm on August 12, 2016, at the onset of high tide from low tide. The photogrammetric process was implemented by the PhotoScan to obtain the DEM of the site. The camera was self-calibrated during the alignment of the image without GCPs. At this stage, corresponding points were identified on the images and matched; in addition, the position of the camera for each image from onboard global navigation satellite system data and camera calibration parameters were automatically refined and stored in the EXIF metadata of the images [25]. With the dense stereo matching of the PhotoScan, a 10 cm resolution DEM was generated by the UAV images (Figure 4). Figure 5 shows the produced UAV DEM and LiDAR DEM (resolution 1 m; obtained in 2011) used as the reference DEM. The UAV and LiDAR DEMs are both digital surface models that represent the mean sea level elevations of the reflective surfaces of all features. The difference between the two DEMs due to vegetation is not high, because most of the test area is covered by tidal flats and the data includes roads, paddy fields, and wild land.
The IOPs and EOPs of the camera influence the accuracy of the DEM derived from UAV photogrammetry [26]. The EOPs were calculated by the position of the GPS equipped on the DJI Phantom 4 and bundle adjustment of the PhotoScan; however, these were not accurate. GCPs must be installed and surveyed on the tidal flat area to obtain the accurate DEM by correction of the IOPs and EOPs. However, the surveying of GCPs is almost impossible, as access to the muddy tidal flat of the Southwest coast is very difficult. Therefore, in this study, the produced DEM was corrected while using a pre-secured LiDAR DEM with a new LHD DEM matching method. In the experiment, we applied the diagonal direction (Line A in Figure 5) as the flight path of the UAV DEM to determine the polynomial models of the transformation. Next, we used the new LHD matching method with the following eleven cases to obtain the transformation parameters:
  • Case 1: 7 parameters ( X 0 ,   Y 0 , Z 0 ,     ω 0 ,   φ 0 ,   κ 0 , s )
  • Case 2: 10 parameters ( X 0 ,   Y 0 , Z 0 , X 1 ,   Y 1 , Z 1 , ω 0 , φ 0 , κ 0 , s )
  • Case 3: 13 parameters ( X 0 ,   Y 0 , Z 0 , X 1 ,   Y 1 , Z 1 , ω 0 , φ 0 , κ 0 , ω 1 , φ 1 , κ 1 , s )
  • Case 4: 13 parameters ( X 0 ,   Y 0 , Z 0 , X 1 ,   Y 1 , Z 1 , X 2 ,   Y 2 , Z 2 , ω 0 , φ 0 , κ 0 , s )
  • Case 5: 16 parameters ( X 0 ,   Y 0 , Z 0 , X 1 ,   Y 1 , Z 1 , X 2 ,   Y 2 , Z 2 , ω 0 , φ 0 , κ 0 ,     ω 1 , φ 1 , κ 1 , s )
  • Case 6: 19 parameters ( X 0 ,   Y 0 , Z 0 , X 1 ,   Y 1 , Z 1 , X 2 ,   Y 2 , Z 2 , ω 0 , φ 0 , κ 0 , ω 1 , φ 1 , κ 1 , ω 2 , φ 2 , κ 2 , s )
  • Case 7: 22 parameters ( X 0 ,   Y 0 , Z 0 , X 1 ,   Y 1 , Z 1 ,   X 2 ,   Y 2 , Z 2 ,   ω 0 , φ 0 , κ 0 ,   ω 1 , φ 1 , κ 1 , ω 2 , φ 2 , κ 2 , ω 3 , φ 3 , κ 3 , s )
  • Case 8: 16 parameters ( X 0 ,   Y 0 , Z 0 , X 1 ,   Y 1 , Z 1 , X 2 ,   Y 2 , Z 2 , X 3 , Y 3 , Z 3 ,   ω 0 , φ 0 ,   κ 0 , s )
  • Case 9: 19 parameters ( X 0 ,   Y 0 , Z 0 , X 1 ,   Y 1 , Z 1 , X 2 ,   Y 2 , Z 2 , X 3 , Y 3 , Z 3 , ω 0 , φ 0 , κ 0 ,   ω 1 , φ 1 ,   κ 1 , s )
  • Case 10: 22 parameters ( X 0 ,   Y 0 , Z 0 ,   X 1 ,   Y 1 , Z 1 , X 2 ,   Y 2 , Z 2 , X 3 , Y 3 , Z 3 ,   ω 0 , φ 0 , κ 0 , ω 1 , φ 1 ,   κ 1 , ω 2 , φ 2 , κ 2 , s )
  • Case 11: 25 parameters ( X 0 ,   Y 0 , Z 0 ,   X 1 ,   Y 1 , Z 1 , X 2 ,   Y 2 , Z 2 ,   X 3 , Y 3 , Z 3 , ω 0 , φ 0 , κ 0 , ω 1 , φ 1 , κ 1 , ω 2 , φ 2 , κ 2 , ω 3 , φ 3 , κ 3 , s )
Table 2 shows the estimated constant terms (six parameters) of the polynomial model in all cases. The results of scale factor are almost 1.0 in all cases; therefore, we have omitted the scale representation. Standard deviations of shift and rotation parameters that are determined by the least squares solution do not exceed 0.3 m and 0.1°, respectively, except Cases 10 and 11. In particular, the six parameters of Case 4 and Case 8 are similar; however, Case 4 has a higher accuracy than Case 8. Overall, Cases 10 and 11 are the least reliable. The degree of the polynomial model does not necessarily have to be high or overparameterized for using it with this method; this shows that the reliability of the parameter can be decreased.
The produced DEM (Figure 6a) was modified by the transformation parameters through the proposed LHD matching that was based on the eleven cases using LiDAR DEM as the reference DEM (Figure 6b). Table 3 shows a quantitative comparison of each modified DEM with LiDAR DEM. The error before DEM correction is 1.4 m with a maximum error of 7.8 m. After DEM correction, the RMSE and maximum error fell to 1.2 m and 6.8 m, respectively, for Case 1. The accuracy slightly improved from Case 1 to Case 3; the accuracy was not improved for Case 4. The error in the DEM is relatively large when compared to DEM resolution (10 cm) due to the inflow water from the tidal flat. There is a possibility that tie-point, self-calibration, and block-adjustment errors in PhotoScan due to the reflection of the underwater area during the process of stereo matching affect the bowl effect.
Figure 7 shows a visual comparison before and after correction for Cases 1, 2, 3, 10, and 11. The left side of Figure 7 shows the point clouds from the UAV images on the mesh of the LiDAR DEM in the lateral view of the 3-D surface, and the right side shows a comparison of the deviations in height between the UAV DEMs and LiDAR DEM.
On the left side of Figure 7, it can be seen that the original UAV DEMs (Figure 7a), and the corrected Cases 1 and 2 DEMs (Figure 7b,c), are almost parabolic in shape (dashed lines), and point clouds are inconsistent with the LiDAR mesh (elliptical shape), whereas the corrected DEMs in Cases 3, 10, and 11 (Figure 7d–f) closely match the LiDAR mesh. In the color maps on the right side of Figure 7, Cases 1 and 2 DEMs are lower than the LiDAR DEM in the central part and higher on the outer part, whereas the heights of the DEMS for Cases 3, 10, and 11 are closer to that of the LiDAR DEM. As the corrected DEMs for Cases 3–11 were almost similar, these comparisons were omitted.
For a detailed comparison of the corrected DEMs, Figure 8 presents the height deviation in each of the cases between the LiDAR DEM and UAV DEM before and after the correction in the Line B profile, as shown in Figure 5. The analysis of accuracy was performed while using the height deviation between the two profiles of the corrected DEMs and LiDAR DEM.
In Figure 8, the errors in Cases 1 and 2 were reduced when compared to the original UAV DEM. However, the error due to the bowl effect still existed even after the correction by applying parameters 7 and 10. Yet, it was observed that this error was removed by applying parameter 13 to Case 3. Therefore, if the transformation parameter is applied by the LHD method, which does not take the polynomial as a flight path of the DEM without the GCP, the bowl effect cannot be eliminated. The optimal parameters were selected from among the 11 cases by calculating the RMSE in the height deviation between the two DEMs after correction (Table 3 and Table 4). In particular, the error sharply decreases from Case 3 to approach the spatial resolution of the UAV image (10 cm) as compared to that before the correction, and the difference in accuracy from Case 3 to Case 11 is less than cm, as shown in Table 4. Applying more parameters showed no significant effect. Table 2 shows that the parameters tended to become slightly inferior as higher-order polynomials were used (Case 10 and Case 11). Therefore, in this study, it was surmised that the optimal parameter to correct the UAV-made DEM while using the reference DEM and proposed LHD matching method was Case 4, which was similar to the bowl effect geometry. It is the secondary curved polynomial with constant rotation having the principal axis in the direction of line, as shown in Figure 5. The experimental results verified that the proposed technique reduced the error of the DEM that was produced by the UAV using only the reference DEM and also eliminated the error due to the bowl effect without requiring GCPs.

4. Conclusions

It is important to obtain a DEM of a tidal flat periodically in order to monitor its environmental condition. It is difficult to acquire GCPs in the tidal flat, although UAV images for accurate and economic DEM production are appropriate. Therefore, we proposed a polynomial LHD DEM matching method by which the position and rotation distortion and bowl effect error were removed, without the requirement of GCPs, in a tidal flat DEM produced from UAV images using a photogrammetry software. The test site chosen to apply the proposed method was a tidal flat area in Hampyung Bay that was located on the west coast of the Korean Peninsula. The experiments were conducted in 11 cases with different correction models to determine the optimal LHD polynomial model. In the experimental investigations, we proposed an optimal model, in which the translation parameters were second-order polynomials and rotation parameters were constant as the number of DEM grids increased in the flight direction. The RMSE in the height of this case was nearly 15 cm. This result was close to the DEM resolution (10 cm) from the acquired UAV image. The experimental results showed that the proposed method could reduce the error in the tidal flat DEM while using only the reference DEM and eliminated the bowl effect error. This method is very effective for correction in the DEM produced by UAV images obtained in inaccessible areas where GCP survey is difficult.
The limitation of the method that is proposed in this paper is that the accuracy of the corrected DEM is affected by the accuracy of the reference DEM. As the height accuracy of LiDAR DEM used as the reference DEM is 15–20 cm, the UAV DEM can be regarded as having similar or lower accuracy, and it is necessary to use an independent reference point to evaluate the actual accuracy of the corrected UAV DEM.

Author Contributions

H.L. conceived the idea, designed the framework; H.L. and D.H. carried out the experiments and analyzed the data; H.L. and D.H. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Ministry of Education of the Republic of Korea and the National Research Foundation of Korea (NRF-2018R1D1A1B06049484).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Marani, M.; D’Alpaos, A.; Lanzoni, S.; Carniello, L.; Rinaldo, A. Biologically-controlled multiple equilibria of tidal landforms and the fate of the Venice lagoon. Geophys. Res. Lett. 2017, 34, 1–5. [Google Scholar] [CrossRef]
  2. Kim, K.L.; Kim, B.J.; Lee, Y.K.; Ryu, J.H. Generation of a large-scale surface sediment classification map using unmanned aerial vehicle (UAV) data: A case study at the Hwang-do tidal flat, Korea. Remote Sens. 2019, 11, 229. [Google Scholar] [CrossRef] [Green Version]
  3. Baker, C.; Lawrence, R.L.; Montagne, C.; Patten, D. Change detection of wetland Ecosystems using Landsat imagery and change vector analysis. Wetlands 2007, 27, 610–619. [Google Scholar] [CrossRef]
  4. Koh, C.H.; Khim, J.S. The Korean tidal flat of the Yellow Sea: Physical setting, ecosystem and management. Ocean Coast. Manag. 2014, 102, 398–414. [Google Scholar] [CrossRef]
  5. Won, J.S.; Na, Y.H.; Kim, S.W. Tidal flat DEM generation by satellite remote sensing. In Proceedings of the 2003 International Geoscience and Remote Sensing Symposium, Toulouse, France, 21–25 July 2003. [Google Scholar]
  6. Ryu, J.H.; Kim, C.H.; Lee, Y.K.; Won, J.S.; Chun, S.S.; Lee, S. Detecting the intertidal morphologic change using satellite data. Estuar. Coast Shelf Sci. 2008, 78, 623–632. [Google Scholar] [CrossRef]
  7. Jaud, M.; Grasso, F.; Dantec, N.L.; Verney, R.; Delacourt, C.; Ammann, J.; Deloffre, J.; Grandjean, P. Potential of UAVs for monitoring Mudflat Morphodynamics (Application to the Seine Estuary, France). ISPRS Int. J. Geo Inf. 2016, 5, 50. [Google Scholar] [CrossRef] [Green Version]
  8. Wimmer, C.; Siegmund, R.; Schwäbisch, M.; Moreira, J. Generation of high precision DEMs of the Wadden Sea with airborne interferometric SAR. IEEE Geosci. Remote Sens. 2000, 38, 2234–2245. [Google Scholar] [CrossRef]
  9. Ross, P.H.; Ustin, S.L.; Hastings, A. Use of lidar to study changes associated with Spartina invasion in San Francisco Bay Marshes. ISPRS J. Photogramm. 2006, 100, 295–306. [Google Scholar] [CrossRef]
  10. Irish, J.L.; Lillycrop, W.J. Scanning laser mapping of the coastal zone: The SHOALS system. ISPRS J. Photogramm. 1999, 54, 123–129. [Google Scholar] [CrossRef]
  11. Kim, D.J.; Moon, W.M.; Kim, G.B.; Park, S.E.; Lee, H.S. Submarine groundwater discharge in tidal flats revealed by space-borne synthetic aperture radar. Remote Sens. Environ. 2011, 115, 793–800. [Google Scholar] [CrossRef]
  12. Mason, D.; Hill, D.; Davenport, I.; Flather, R.; Robinson, G. Improving inter-tidal digital elevation models constructed by the waterline technique. In Proceedings of the Third ERS Symposium, European Space Agency Publications Division, Florence, Italy, 14–21 March 1997; pp. 1079–1082. [Google Scholar]
  13. Choi, J.K.; Ryu, J.H.; Lee, Y.K.; Yoo, H.R.; Woo, H.J.; Kim, C.H. Quantitative estimation of intertidal sediment characteristics using remote sensing and GIS. Estuar. Coast Shelf Sci. 2010, 88, 125–134. [Google Scholar] [CrossRef]
  14. Kim, D.J.; Moon, W.M.; Park, S.E.; Kim, J.E.; Lee, H.S. Dependence of water-line mapping on radar frequency used for SAR images in intertidal areas. IEEE Geosci. Remote Sens. Lett. 2007, 4, 269–273. [Google Scholar] [CrossRef]
  15. Long, N.; Millescamps, B.; Guillot, B.; Pouget, F.; Bertin, X. Monitoring the topography of a dynamic tidal inlet using UAV imagery. Remote Sens. 2016, 8, 387. [Google Scholar] [CrossRef] [Green Version]
  16. EcoSea. Available online: http://www.ecosea.go.kr (accessed on 3 December 2019).
  17. Ouédraogo, M.M.; Degré, A.; Debouche, C.; Lisein, J. The evaluation of unmanned aerial system-based photogrammetry and terrestrial laser scanning to generate DEMs of agricultural watersheds. Geomorphology 2014, 214, 339–355. [Google Scholar] [CrossRef]
  18. James, M.R.; Robson, S. Mitigating systematic error in topographic models derived from UAV and ground-based image networks. Earth Surf. Process. Landf. 2014, 39, 1413–1420. [Google Scholar] [CrossRef] [Green Version]
  19. Jaud, M.; Passot, S.; Bivic, R.L.; Delacourt, C.; Grandjean, P.; Dantec, N.L. Assessing the accuracy of high resolution digital surface models computed by PhotoScan® and MicMac® in sub-optimal survey conditions. Remote Sens. 2016, 8, 465. [Google Scholar] [CrossRef] [Green Version]
  20. Forlani, G.; Dall’Asta, E.; Diotri, F.; Cella, U.M.D.; Roncella, R.; Santise, M. Quality assessment of DSMs produced from UAV flights georeferenced with on-board RTK positioning. Remote Sens. 2018, 10, 311. [Google Scholar] [CrossRef] [Green Version]
  21. Taddia, Y.; Stecchi, F.; Pellegrinelli, A. Using DJI PHANTOM 4 RTK drone for topographic mapping of coastal areas. In Proceedings of the International Archives of Photogrammetry, Remote Sensing & Spatial Information Sciences, Enschede, The Netherlands, 10–14 June 2019. [Google Scholar]
  22. Rosenholm, D.; Torlegård, K. Three-dimensional absolute orientation of stereo models using digital elevation models. Photogramm. Eng. Remote Sens. 1988, 54, 1385–1389. [Google Scholar]
  23. Ebner, H.; Strunz, G. Combined point determination using digital terrain models as control information. In Proceedings of the International Archives of Photogrammetry and Remote Sensing, Kyoto, Japan, 1–10 July 1988. [Google Scholar]
  24. Karras, G.; Petsa, E. DEM matching and detection of deformation in close-range photogrammetry without control. Photogramm. Eng. Remote Sens. 1993, 59, 1419–1424. [Google Scholar]
  25. Agisoft, L.L.C. Agisoft PhotoScan User Manual. Available online: https://www.agisoft.com/pdf/photoscan-pro14en.pdf (accessed on 30 December 2019).
  26. Gindraux, S.; Boesch, R.; Farinotti, D. Accuracy assessment of digital surface models from unmanned aerial vehicles’ imagery on glaciers. Remote Sens. 2017, 9, 186. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Least-squares height-difference (LHD) matching between a reference digital elevation model (DEM) and a unmanned aerial vehicle (UAV) DEM from an UAV (a) and a common result from the UAV DEM, as affected by the bowl effect (b).
Figure 1. Least-squares height-difference (LHD) matching between a reference digital elevation model (DEM) and a unmanned aerial vehicle (UAV) DEM from an UAV (a) and a common result from the UAV DEM, as affected by the bowl effect (b).
Sensors 20 01602 g001
Figure 2. Flow-chart for proposed LHD matching and correction of DEM.
Figure 2. Flow-chart for proposed LHD matching and correction of DEM.
Sensors 20 01602 g002
Figure 3. Test location (a) and orthomosaic image taken by the camera onboard the UAV (b).
Figure 3. Test location (a) and orthomosaic image taken by the camera onboard the UAV (b).
Sensors 20 01602 g003
Figure 4. PhotoScan process for DEM (a); camera locations and image overlap in PhotoScan processing report (b).
Figure 4. PhotoScan process for DEM (a); camera locations and image overlap in PhotoScan processing report (b).
Sensors 20 01602 g004
Figure 5. UAV DEM (a), Line A is the direction for the polynomial models of the transformation parameters, Line B is the profile for error check, and light detection and ranging (LiDAR) DEM (b).
Figure 5. UAV DEM (a), Line A is the direction for the polynomial models of the transformation parameters, Line B is the profile for error check, and light detection and ranging (LiDAR) DEM (b).
Sensors 20 01602 g005
Figure 6. (a) UAV DEM and (b) LiDAR DEM.
Figure 6. (a) UAV DEM and (b) LiDAR DEM.
Sensors 20 01602 g006
Figure 7. DEM (a) before correction and after correction for (b) Case 1, (c) Case 2, (d) Case 3, (e) Case 10, and (f) Case 11. The figures on the left are overlaid with UAV point clouds and LiDAR mesh; the figures on the right are color maps of height deviations between UAV DEMs and LiDAR DEM (in meters).
Figure 7. DEM (a) before correction and after correction for (b) Case 1, (c) Case 2, (d) Case 3, (e) Case 10, and (f) Case 11. The figures on the left are overlaid with UAV point clouds and LiDAR mesh; the figures on the right are color maps of height deviations between UAV DEMs and LiDAR DEM (in meters).
Sensors 20 01602 g007aSensors 20 01602 g007b
Figure 8. (a) DEM profiles before and after correction using the UAV images and LiDAR DEM for (b) Case 1, (c) Case 2, (d) Case 3, (e) Case 10, and (f) Case 11. The circled section shows the gross error in the UAV DEM due to tide inflows.
Figure 8. (a) DEM profiles before and after correction using the UAV images and LiDAR DEM for (b) Case 1, (c) Case 2, (d) Case 3, (e) Case 10, and (f) Case 11. The circled section shows the gross error in the UAV DEM due to tide inflows.
Sensors 20 01602 g008
Table 1. Technical data of Phantom 4 and camera.
Table 1. Technical data of Phantom 4 and camera.
ItemsParameters
Phantom 4Weight1380 g
Max. speed20 m/s
Flight time28 min
GNSSGPS, GLONASS
CameraFocal length3.61 mm
FOV94°
Sensor size12.4 M (4000 × 3000)
Table 2. Six transformation parameters and their standard deviations.
Table 2. Six transformation parameters and their standard deviations.
Case #Parameters
X 0 / σ X 0 ( m ) Y 0 / σ Y 0 ( m ) Z 0 / σ Z 0 ( m ) ω 0 / σ ω 0   ( ° ) φ 0 / σ φ 0 ( ° ) κ 0 / σ κ 0 ( ° )
1−0.59/0.0120.33/0.013−0.03/0.0040.34/0.0020.32/0.002−0.02/0.003
2−0.18/0.0700.66/0.07410.85/0.1311.19/0.008−0.51/0.008−0.02/0.015
3−0.59/0.0120.33/0.01311.27/0.0040.34/0.0020.32/0.001−0.02/0.003
40.42/0.0571.78/0.0618.88/0.0741.34/0.0110.64/0.0120.01/0.008
51.86/0.0850.19/0.0802.97/0.2250.81/0.018−0.18/0.017−0.21/0.060
62.22/0.0870.25/0.0802.94/0.2270.84/0.019−0.15/0.017−0.46/0.060
71.95/0.0980.41/0.0853.65/0.2331.21/0.021−0.05/0.019−0.29/0.061
82.27/0.1152.67/0.1208.87/0.0761.35/0.0420.64/0.0460.03/0.012
9−0.86/0.1301.16/0.1272.93/0.2260.81/0.044−0.17/0.043−0.20/0.060
10−0.48/0.2331.63/0.19815.76/0.6001.85/0.083−1.18/0.072−0.19/0.241
11−1.07/0.2371.86/0.19817.79/0.6102.32/0.083−1.18/0.0710.05/0.243
Table 3. Height deviations between the LiDAR DEM and UAV DEM before and after correction.
Table 3. Height deviations between the LiDAR DEM and UAV DEM before and after correction.
Correction BeforeAfter: Case #
1234567891011
RMSE (m)1.371.241.211.121.111.111.101.111.101.101.101.10
Max (m)7.846.766.285.996.005.985.965.955.985.955.955.94
Table 4. Height deviations between the two DEMs before and after correction on the line B profile in Figure 5.
Table 4. Height deviations between the two DEMs before and after correction on the line B profile in Figure 5.
CorrectionBeforeAfter: Case #
1234567891011
RMSE (m)0.9560.5170.4750.1500.1490.1530.1560.1660.1520.1570.1490.163
Max (m)2.2582.2182.0601.4681.4401.3851.3931.4261.4571.4031.4141.451

Share and Cite

MDPI and ACS Style

Lee, H.; Han, D. Rectification of Bowl-Shape Deformation of Tidal Flat DEM derived from UAV Imaging. Sensors 2020, 20, 1602. https://doi.org/10.3390/s20061602

AMA Style

Lee H, Han D. Rectification of Bowl-Shape Deformation of Tidal Flat DEM derived from UAV Imaging. Sensors. 2020; 20(6):1602. https://doi.org/10.3390/s20061602

Chicago/Turabian Style

Lee, Hyoseong, and Dongyeob Han. 2020. "Rectification of Bowl-Shape Deformation of Tidal Flat DEM derived from UAV Imaging" Sensors 20, no. 6: 1602. https://doi.org/10.3390/s20061602

APA Style

Lee, H., & Han, D. (2020). Rectification of Bowl-Shape Deformation of Tidal Flat DEM derived from UAV Imaging. Sensors, 20(6), 1602. https://doi.org/10.3390/s20061602

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop