Next Article in Journal
Drone-Based Atmospheric Soundings Up to an Altitude of 10 km-Technical Approach towards Operations
Previous Article in Journal
Consensus Control of Large-Scale UAV Swarm Based on Multi-Layer Graph
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Experimental Study on an Innovative UAV-LiDAR Topographic Mapping System for Precision Land Levelling

1
College of Agricultural Equipment Engineering, Henan University of Science and Technology, Luoyang 471003, China
2
Department of Agriculture and Natural Resources, University of Mohaghegh Ardabili, Ardabil 566199, Iran
*
Author to whom correspondence should be addressed.
Drones 2022, 6(12), 403; https://doi.org/10.3390/drones6120403
Submission received: 9 November 2022 / Revised: 1 December 2022 / Accepted: 6 December 2022 / Published: 8 December 2022
(This article belongs to the Section Drones in Agriculture and Forestry)

Abstract

:
Topographic maps provide detailed information on variations in ground elevation, which is essential for precision farmland levelling. This paper reports the development and experimental study on an innovative approach of generating topographic maps at farmland-level with the advantages of high efficiency and simplicity of implementation. The experiment uses a low-altitude Unmanned Aerial Vehicle (UAV) as a platform and integrates Light Detection and Ranging (LiDAR) distance measurements with Post-Processing Kinematic Global Positioning System (PPK-GNSS) coordinates. A topographic mapping experiment was conducted over two fields in Henan Province, China, and primitive errors of the topographic surveying data were evaluated. The Root Mean Square Error (RMSE) between elevation data of the UAV-LiDAR topographic mapping system and ground truth data was calculated as 4.1 cm and 3.6 cm for Field 1 and Field 2, respectively, which proved the feasibility and high accuracy of the topographic mapping system. Furthermore, the accuracies of topographic maps generated using different geo-spatial interpolation models were also evaluated. The results showed that a TIN (Triangulated Irregular Network) interpolation model expressed the best performances for both Field 1 with sparse topographic surveying points, and Field 2 with relatively dense topographic surveying points, when compared with other interpolation models. Moreover, we concluded that as the spatial resolution of topographic surveying points is intensified from 5 m × 0.5 m to 2.5 m × 0.5 m, the accuracy of the topographic map based on the TIN model improves drastically from 7.7 cm to 4.6 cm. Cut-fill analysis was also implemented based on the topographic maps of the TIN interpolation model. The result indicated that the UAV-LiDAR topographic mapping system could be successfully used to generate topographic maps with high accuracy, which could provide instructive information for precision farmland levelling.

1. Introduction

Surface unevenness of farmlands has been pointed out as a major issue that affects agricultural drainage efficiency [1]. Hu et al. (2014) reported that over 20% of irrigation water could be wasted because of the rough land surface of paddy fields. Rickman (2002) concluded that unleveled rice fields suffered about 24% yield loss on average when compared to effectively leveled ones of the same rice variety and the same fertilizer input [2,3]. As precision farmland levelling facilitates agricultural drainage, improves crop establishment, and increases crop yield, many researchers both in academia and industry have proposed different methodologies for removing mounds and puddles in farmland [4,5]. Generally, two major precision farmland levelling techniques can be enumerated: laser farmland levelling and GNSS (Global Navigation Satellite System) farmland levelling [6]. An accurate topographic map is essential to precision farmland levelling; this map type refers to a 2D (two-dimensional) graphic representation of a terrestrial land surface feature using contour lines, hypsometric tints, and relief shading [7,8,9]. As for the laser leveler, the preset height of the laser transmitter directly determines the cutting depth of each point, i.e., cut-fill ratios. It is almost impossible to set an appropriate height efficiently and accurately without referring to an accurate topographic map [10]. On the other hand, the GNSS leveler measures and adjusts the height of the scraper by comparing the current height of the GNSS receiver with the reference height plane to complete the earth-moving operation. It also needs an accurate topographic map to form the reference height plane. Thus, the accuracy, efficiency, and energy consumption of precision farmland levelling are in high accordance with the delicacy of farmland topographic maps [10].
To generate topographic maps for civil construction surveys, urban ecology modeling, forest monitoring, etc., new technologies including terrestrial laser scanning, aerial photogrammetry, and airborne laser scanning were recently utilized [11,12,13,14,15,16,17,18], while theodolite, total station, and handheld RTK (Real-Time Kinematic)-GNSS modules remain the conventional and primary tools used in common topographic surveys [19,20,21]. Topographic survey using a handheld GNSS module or tripod-based total station has low efficiency and cannot be applied to large-area operations [22]. On the other hand, topographic survey via remote sensing covers large areas and is capable of obtaining adequate data, in spite of the high cost and complicated calibration procedures [23]. Corsini et al. (2013) monitored and mapped a slow-moving compound rockslide using an integration of an airborne laser scanner, terrestrial laser scanner, and automated total station, which quantified slope movement in the order of centimeters to a few decimeters [24]. Rodriguez et al. (2017) evaluated a mobile LiDAR (Light Detection and Ranging) system mounted on a car to develop an architectural analysis by generating a 3D point cloud [25]. Terrone et al. (2021) identified manmade landforms and assessed the morphological evolution of the city of Genoa by coupling historical maps and LiDAR data [26]. There are also plenty of research works and products providing a DSM (Digital Surface Model) based on airborne photogrammetry or satellite stereo-imagery [27,28,29,30,31,32,33]. However, processing aerial images to generate photogrammetric DSMs usually needs strict camera calibration and a large number of spatially well-distributed GCPs (Ground Control Points), which in photogrammetry or the computer vision domain refers to such features that are easily recognizable and distinguishable in both the real world and the images. The spatial resolution, as well as the accuracy, of such topographic maps usually reaches several decimeters to tens of meters [34], which cannot meet the high requirements of topographic maps with centimeter-level accuracy for precision farmland levelling.
Therefore, this study developed and experimentally evaluated an innovative topographic mapping system for precision farmland levelling that was based on a low-altitude UAV (Unmanned Aerial Vehicle). The UAV was equipped with a LiDAR device and PPK (Post-Processing Kinematic)-GNSS modules, which are utilized to conduct topographic survey in a simple and totally autonomous manner. This research integrates multi-source remote sensing data on board the UAV platform, including LiDAR distance measurements, attitude information (pitch and roll) of the UAV’s flight controller, and PPK-GNSS positioning data (latitude, longitude, altitude). The ultimate objective is to generate topographic maps of farmlands with centimeter-level accuracy efficiently, both in terms of time and cost, using geo-spatial interpolation models.

2. Materials and Methods

2.1. Experiment Equipment

A farmland topographic mapping experiment was established over two fields located in Kaifeng City, Henan Province, China. Each field accounts for about 1500 m2 and 1200 m2, respectively, shown in Figure 1 as Field 1 and Field 2. The previous crop was maize. Weeds in the fallow fields were manually removed prior to the experiment so that there were no significant foreign attachments on the ground surface.
A hexacopter was used as the UAV platform, shown in Figure 2a. A high-accuracy LiDAR distance measuring device (JENOPTIK, Jena, Germany) was used to measure the distance between the ground surface and the UAV-LiDAR topographic mapping system, which was rigidly fixed under the UAV platform pointing vertically downwards. Two GNSS modules were used to calculate the 3D positioning coordinates at the frequency of 10 Hz of the UAV-LiDAR topographic mapping system in the PPK manner. One GNSS module was installed on the top of the UAV as a rover receiver, while the other was fixed nearby the field as a base receiver, shown in Figure 2b. The parameters of the equipment of the UAV-LiDAR topographic mapping system are listed in Table 1.
Autonomous UAV flight was conducted on 24 October 2020, using flight paths designed beforehand at the speed of about 5 m/s and altitude of about 30 m above ground level. The flight speed was determined with consideration of the efficiency and attitude stability of the UAV platform. Low flight speed results in low efficiency, and when the flight speed increases, the pitch angle of UAV body correspondingly changes violently, which affects the stability of the topographic surveying system and accuracy of laser ranging data. From the flight speed and positioning frequency of the GNSS module, in-track intervals of the UAV-LiDAR topographic surveying points could be determined as 0.5 m. The cross-track interval of Field 1 and Field 2 was set at about 5 m and 2.5 m, respectively, to cover the whole field in each flight, according to the endurance limitation of common civilian UAVs. The influence of spatial resolution of the topographic surveying points on the accuracy of topographic maps by using interpolation models was studied; this is jointly determined by the cross-track interval and the in-track interval. As the in-track interval was already far less than the cross-track interval, the merits of further reducing the flight speed for acquiring denser topographic surveying points was not discussed. Therefore, the spatial resolution of the topographic surveying data for Field 1 and Field 2 was calculated as about 5 m × 0.5 m and 2.5 m × 0.5 m, respectively. The trajectories of two autonomous flights over the experimental fields are shown in Figure 1. Because the altitude of the UAV-LiDAR topographic surveying system affects the accuracy of LiDAR’s distance measurement, it should be set as low as possible. In this study, the altitude was set to 30 m above ground level to conduct accurate topographic surveying and to simultaneously capture aerial images for further study of data integration.

2.2. Acquiring PPK-GNSS Coordinates

Theoretically, a standalone GNSS receiver usually has a positioning accuracy varying from submeter to several meters [35,36]. On the other hand, an RTK-GNSS module uses a network of ground or virtual reference stations to rectify the GNSS rover receiver’s positioning data, and its positioning accuracy could be significantly improved up to about 2~5 cm. However, as the RTK-GNSS device is embedded with communication modules for data transmission via radio signal or cellular network, it is always too large in size and/or weight to be installed upon small-sized UAV. Therefore, in this study two identical lightweight PPK-GNSS modules were used to collect high-accuracy positioning data. After the completion of a UAV topographic mapping flight, GNSS positioning data of the rover receiver and the base receiver were acquired; we could dynamically rectify the GNSS rover receiver’s positioning data in the post-processing style with the reference of the GNSS base receiver’s positioning data in the following manner.
Firstly, both the GNSS base and rover receiver were held stationary before the UAV platform took off for over 30 min consecutively, and therefore 18,000 sets of the effective static positioning data of each GNSS module at the frequency of 10 Hz were obtained. By using an arithmetic averaging method, the reference coordinate ( x r , y r , h r ) for PPK-GNSS algorithm was calculated from the GNSS base receiver’s positioning data ( x b a s e i , y b a s e i , h b a s e i ) as (3,848,750.598 E, 312,414.385 N, 63.754 m), according to Equations (1)–(3).
x r = i = 1 18000 x b a s e i 18000
y r = i = 1 18000 y b a s e i 18000
h r = i = 1 18000 h b a s e i 18000
where x r , y r , and h r refer to easting, northing, and altitude of the reference coordinate for the PPK algorithm, while x b a s e i , y b a s e i , and h b a s e i refer to easting, northing, and altitude of the positioning data of the GNSS base receiver, respectively.
Subsequently, the positioning data of the GNSS rover receiver ( x r o v i , y r o v i , h r o v i ) and the base receiver over the same period were processed in the RTKLIB development environment using the PPK-GNSS algorithm, where x r o v i , y r o v i , and h r o v i refer to easting, northing, and altitude of positioning data of the GNSS rover receiver, respectively. The PPK-GNSS algorithm calculates position deviations between each positioning datum of the GNSS base receiver ( x b a s e i , y b a s e i , h b a s e i ) and the reference coordinate ( x r , y r , h r ) based on carrier phase analysis. It then uses the deviation information to rectify the corresponsive positioning data of GNSS rover receiver ( x r o v i , y r o v i , h r o v i ) . Thus, PPK-GNSS positioning data were acquired as ( x G P S i , y G P S i , h G P S i ) .
Finally, the GNSS base receiver was kept at the same location, and the GNSS rover receiver on board the UAV platform was used to conduct the topographic mapping experiment. The UAV flight lasted for about 12 min for each field, and in total 675 and 1626 sets of PPK-GNSS positioning data were obtained for Field 1 and Field 2, respectively, after removing noise data during take-off, landing, and turn-around. The horizontal positioning coordinates of PPK-GNSS data were previously given in Figure 1, while the altitudes of PPK-GNSS data are shown in Figure 3. In Figure 3, variations in the UAV’s flight altitude can be observed from 90 m to 92.5 m, which is caused by its intrinsic aerodynamic factors and the influence of turbulence. Upon completion of the topographical mapping flight, the GNSS rover receiver was removed from the UAV platform and fixed onto an aluminum plate. Subsequently, it was dragged manually over the field ground to cross the UAV’s flight paths, shown in Figure 1 as green lines, in order to acquire PPK-GNSS data for accuracy evaluation.

2.3. Rectifying LiDAR Distance Measurements Using UAV Attitude Data

LiDAR is based on the time-of-flight principle and has a measuring range up to 300 m and a 1 mm measuring resolution at the frequency of 2000 Hz [37,38,39,40]. The laser beam divergence of the LiDAR distance measurement device was 1.7 milliradians, and the spatial coverage of each measurement was about 60 mm × 20 mm from a distance of about 30 m. An onboard computer was connected with the LiDAR through a universal serial bus cable; it saves each distance measurement on board in real time. As the output frequency of the PPK-GNSS’s positioning data was 10 Hz, a mean filter with a 200 step size was applied to raw LiDAR distance measurements. Thus, the output frequency of the LiDAR distance measurements was altered from 2000 Hz to 10 Hz in order to align multi-sources remote sensing data according to time sequences and also to improve the LiDAR distance measuring accuracy in the meantime. Due to its aerodynamic characteristics, the UAV body is apt to tilt irregularly during cruise flight [41]. Because electronical gimbals cannot meet the real-time stabilizing requirements of this study, the LiDAR was rigidly fixed below the UAV platform and the UAV’s attitude information was used to compensate for such influences as vibration and air disturbance during flight.
Attitude of a UAV is usually described by means of Euler angles defined as pitch (θ), roll (ϕ), and yaw (ψ) [42,43,44]. In this study, we utilized an extended Kalman filter to integrate data from multiple sensors of the MEMS (Micro-Electro-Mechanical System) gyroscope, accelerometer, and magnetometer for attitude estimation. A gyroscope is an inertial sensor for measuring orientation based on angular momentum principles, and tri-axial angles [ θ g y r o   ϕ g y r o   ψ g y r o ] T can be acquired by integral operation, as expressed in Equation (4). However, due to temperature variations in the gyroscope, errors accumulate along with time and the accuracy of the MEMS gyroscope’s attitude data will be compromised. On the other hand, an accelerometer was used to measure the UAV’s orientation based on the trigonometric functions of the acceleration of gravity components in each axial, expressed in Equation (5) and Equation (6). A MEMS magnetometer is also used for precisely calculating change of the UAV’s heading yaw based on magnetic intensity components in each axial [ m x , m y , m z ] T , in combination with measurements from the gyroscope and accelerometer [45,46].
[   θ g y r o   ϕ g y r o   ψ g y r o ] T =   [ ω x   ω y   ω z ] T d t
where θ g y r o ,   ϕ g y r o ,   and   ψ g y r o are pitch, roll, and yaw of the UAV platform from the gyroscope, respectively, while ω x   , ω y ,   and   ω z are the raw rotation rates (degree/s) of each axial from the gyroscope.
θ a c c e = a r c s i n ( g x g x 2 + g y 2 + g z 2 )
ϕ a c c e = a r c t a n ( g y g z )
where θ a c c e and ϕ a c c e are pitch and roll of the UAV from the accelerometer, respectively, while g x , g y , and g z are the acceleration components readings in each axial from the accelerometer.
Therefore, in order to further improve the accuracy of LiDAR distance measurements, attitude data (pitch θ i and roll φ i ) of the UAV platform were used to obtain the nadir distance, according to Equation (7). Nadir is defined as the angle that points directly downward, or 0°, from the luminaire, as shown in Figure 4. Variations in the nadir distance between the LiDAR device and ground surface during the UAV-LiDAR topographic surveying process are shown in Figure 5. The nadir distance between the UAV-LiDAR system and ground surface varied from 25.5 m to 28.5 m for both Field 1 and Field 2. The variations in the nadir distance between the UAV-LiDAR topographic mapping system and ground surface arise from many factors, such as the ever-changing altitude of the UAV platform during flight and the field terrain itself.
              d n i = d i × cos θ i × cos φ i
where d n i and d i are the nadir distance and the actually measured distance between the LiDAR device and ground surface, respectively, while θ i and φ i are the attitude angle of pitch and roll of the UAV flight controller, respectively.

2.4. Acquiring 3D Coordinates of Ground Surveying Points

The scheme for calculating the elevation of ground surveying points ( h i ) is elaborated in Figure 6 and could be calculated according to Equation (8). Subsequently, by replacing h G P S i with h i , 3D coordinates of the ground surveying points could be obtained as ( x G P S i , y G P S i , h i ) .
h i = h G P S i h f i x d n i
where h i , h G P S i , h f i x , and d n i are the elevation of ground surveying points, the altitude of the UAV’s PPK-GNSS coordinates, the installation height difference between the GNSS rover receiver and the LiDAR device (0.66 m), and the nadir distance between the UAV-LiDAR system and ground surface, respectively.

2.5. Generating Topographic Maps Based on Interpolation Models

In this study, 675 and 1626 sets of effective topographic mapping data for Field 1 and Field 2, respectively, were used to generate topographic maps based on different interpolation models. In geo-statistics, interpolation is a method for approximating and assigning new data for locations where no samples have been taken within the range of a discrete set of known data points. There are various interpolation models such as TIN (Triangulated Irregular Network), IDW (Inverse Distance Weighting), Kriging, and natural neighbor, etc. [47,48,49,50]; however, it remains unclear as to which model is the most suitable for farmland terrain research purposes.
TIN is a vector-based interpolation model that forms a network of triangles of irregular size and shape. TIN’s triangular linear interpolation process begins with determining the triangular to which the interpolation point belongs, then calculating the parameters of a, b, and c for the plane equations of each triangular using the 3D coordinates of three vertices, expressed as Equation (9). The elevation data of every interpolation point within the triangular could be acquired according to Equation (10).
    [ z 1 z 2 z 3 ] = [ 1 ,   x 1 , y 1 1 ,   x 2 , y 2 1 ,   x 3 , y 3 ] [ a b c ]
z = a + b x + c y
where a, b, and c are the parameters for determining the plane equation, while x i , y i , and z i (I = 1, 2, and 3) are the 3D coordinates (easting, northing, and elevation) of each of the vertices of the triangular, and x, y, and z are the 3D coordinates of each interpolating point, respectively.
The IDW interpolation model assigns a weighted value of serval neighboring points by using an inverse distance weighted technique. The weight parameters could be determined from the mathematical model expressed as Equations (11)–(13). The natural neighbor interpolation model is similar to the IDW model, except that it determines weight parameters by calculating proportionate areas instead of inversed distances.
w i = d i p i = 1 n i d w d i p
d i = ( x x i ) 2 + ( y y i ) 2
z ^ ( x , y ) = i = 1 n i d w w i z ( x i ,   y i )
where w i , d i , p , and n i d w are the weighting factor of each measured points, distance of each measured point to interpolating point, power parameter, and number of measured points included in the IDW model, respectively, while ( x ,   y )   and   ( x i , y i ) are easting and northing coordinates of the interpolating point and each measured point, respectively; z ^ ( x , y )   and   z ( x i ,   y i ) are the estimated ground elevation of the interpolating point and each measured point’s ground elevation, respectively.
As mentioned above, TIN, IDW, and natural neighbor interpolation models use the surrounding points’ z–values for estimating the z–value for each interpolating point. However, the Kriging interpolation model incorporates geostatistical relationships among the measured points. Thus, in contrast to the IDW model, the weight parameter for each surrounding measured point in the Kriging interpolation model is not only determined by the inversed distance, but also influenced by the overall spatial autocorrelation of the measured points, expressed as Equation (14).
z ^ ( x , y ) = i = 1 n k r i w i z ( x i ,   y i )
where z ^ ( x , y ) , n k r i ,   w i , and z ( x i ,   y i ) are the estimated z-value of an interpolating point, number of measured points included in the Kriging model, weighting factor of each measured points, and each measured point’s z-value, respectively.

3. Results and Discussion

3.1. Evaluating the Accuracy of a UAV-LiDAR Topographic Mapping System

Horizontal positioning data of the stationary GNSS base receiver in standalone and PPK mode is shown in Figure 7a,b, respectively, from which the horizontal positioning accuracy can be determined as about 3 m and 2 cm. The vertical positioning data varied within the range of 66.180 m to and 66.210 m, with an accuracy of about 3 cm, shown in Figure 7c. Therefore, Figure 7 clearly indicates that PPK-GNSS positioning data have high accuracy and are capable of precisely measuring the 3D coordinates of the UAV-LiDAR topographic mapping system.
Based on the 3D coordinates of each topographic surveying point, topographic maps indicating spatial variations in within-field ground elevation were generated by using ArcMap (ESRI Inc., Redlands, AB, Canada), shown in Figure 8. From the graduated symbols representing different levels of ground elevation, concave features (puddles) at each end of Field 1 can be visually specified and the general high-south-low-north terrain of Field 2, with ground elevations varying from 63.398 m to 63.761 m can be understood.
Accurate 3D coordinates (easting, northing, and elevation) of 12 discrete points distributed around each field were acquired for the accuracy evaluation; these overlapped (or fell very near to) the topographic mapping points. The geo-spatial coordinates of these points are listed in Table 2.
The Root Mean Square Error (RMSE) between 12 pairs of ground elevation data of the handheld PPK-GNSS module and UAV-LiDAR topographic mapping system was calculated as 0.041 m and 0.036 m for Field 1 and Field 2, respectively. The RMSE for the two fields under study showed no significant distinction, although the track interval of the UAV-LiDAR topographic mapping system for the two fields was quite different (5 m for Field 1 and 2.5 m for Field 2). Since the ground elevation of the fields under study varied from 63.339 m to 63.829 m and 63.398 m to 63.761 m, respectively, the value of RMSEUAV showed that the topographic data of the UAV-LiDAR topographic mapping system are of very high precision.

3.2. Topographic Maps of Different Interpolation Models

In ArcMap software, topographic maps representing the ground elevation features of Field 1 and Field 2 using the TIN model were generated, shown in Figure 9. From the maps it may be concluded that ground elevation of Field 1 varies from 63.339 m to 63.829 m, and the distinguishing convex features (ridges) that are shown in white intertwine with concave features (puddles) that are shown in black, with no obvious trend of changing terrain. However, we can visualize a clear low-west-high-east terrain trend for Field 2, with the ground elevation changing from 63.398 m to 63.761 m.
The IDW interpolation result is shown in Figure 10, from which it can be seen more clearly that ridges and puddles spread all around Field 1 and the low-west-high-east terrain trend of Field 2 is also more obvious when compared with the topographic maps based on the TIN interpolation model.
Topographic maps using the natural neighbor interpolation model, shown in Figure 11, were also generated in ArcMap software. From these, the similar terrain details of the two fields under study can be better observed in comparison to with the topographic maps generated by using the IDW interpolation model. Moreover, the topographic maps based on the natural neighbor interpolating model have a more drastic change in ground elevation.
In this study, an ordinary Kriging method with a spherical semi-variogram model was used in ArcMap software to generate topographic maps, shown in Figure 12, from which it may be concluded that the topographic maps based on the Kriging interpolation model are quite similar to the topographic maps based on the IDW model, and the former are featured with more smooth edges in terms of terrain change.

3.3. Evaluating Accuracy of Topographic Maps

The accuracy of topographic maps generated using interpolation models depends on two systems of error. One is the accuracy of each measured point, i.e., the primitive error. The other adheres to a mathematical model of each specific interpolation model, i.e., the interpolation error. In this study, the primitive error is indicated by the values of RMSE between the ground elevation data of the handheld PPK-GNSS module and the UAV-LiDAR topographic mapping system, which proved that the topographic data of the UAV-LiDAR topographic mapping system are very accurate. Interpolation error is the main source that contributes to the accuracy of topographic maps as the z–value of each interpolation point is not only determined by the surrounding points’ z–values, but also is influenced by the neighboring points’ spatial distribution.
In order to evaluate the accuracy of topographic maps generated using different interpolation models, accurate 3D coordinates of 10 discrete points distributed around each field from the trajectory of the handheld PPK-GNSS module were obtained. The spatial distribution of these evaluating points was shown as green dots in Figure 1, and the 3D coordinates were listed in Table 3. Furthermore, according to the easting and northing coordinates of the handheld PPK-GNSS data, the ground elevations of each point in the topographic maps of the different interpolation models were extracted in ArcMap software and listed in Table 3.
The RMSE between the ground elevation data of the handheld PPK-GNSS module and topographic maps of different interpolation models for Field 1 and Field 2 was calculated and listed in Table 4. From Table 4, it could be concluded that although there is no significant distinction between the primitive error of the UAV-LiDAR topographic mapping system for Field 1 and Field 2, the difference of interpolating errors among each interpolation model is obvious. The TIN interpolation model expressed the best performances for both Field 1 with sparse topographic surveying points, and Field 2 with relatively dense topographic surveying points, when compared with the other interpolation models including IDW, natural neighbor, and Kriging. Furthermore, as the spatial resolution of the topographic surveying points is improved from 5 m × 0.5 m to 2.5 m × 0.5 m, the interpolating error of the topographic maps based on the TIN model drops drastically from 0.077 m to 0.046 m.
According to the topographic map of the TIN models that were most suitable to depict farmland terrains, cut-fill analysis was conducted for each field. Based on the statistical results of the ground elevation data from the topographic maps, the average value of each field’s ground elevation data was calculated as 63.597 m and 63.564 m. An average value of ground elevation was assigned to new raster models that shared the same geo-spatial coordinates and resolution with the topographic maps of the TIN model as the desired ground elevation. By comparing the value of ground elevation for each point in the topographic map with the new raster model, operations of cut or fill could be determined, as shown in Figure 13; blue represents the areas where ground elevation is larger than the average value and soils could be cut to form a desired plane, while red represents the areas where ground elevation is smaller than the average value and it needs extra soils to be filled. Furthermore, the results of the cut-fill analysis showed that for Field 1 the volume of cut soil is about 33.2 m3, which accounted for about 616.4 m2 (shown in blue in Figure 13a), while the areas to be filled with extra soils accounted for about 454.8 m2 (shown in red in Figure 13a). On the other hand, for Field 2 the volume of cut soil is about 46.3 m3, which accounted for about 603.8 m2 (showed in blue in Figure 13b), and the areas to be filled with extra soils accounted for about 557.6 m2 (shown in red in Figure 13b).

4. Conclusions

In this study, an innovative topographic mapping system integrating LiDAR distance measurements with PPK-GNSS coordinates was developed using a low-altitude UAV as a platform in a simple, efficient, and totally autonomous fashion. The positioning accuracy of the PPK-GNSS module was experimentally validated as 2 cm, which shows a good capability for conducting topographic mapping of farmlands. Topographic mapping experiments in two fields were conducted and topographic maps indicating spatial variations in within-field ground elevation were generated. The main conclusions were summarized as follows.
(1). The RMSE between 12 pairs of ground elevation data of the handheld PPK-GNSS module and UAV-LiDAR topographic mapping system was calculated as 0.041 m and 0.036 m for Field 1 and Field 2, respectively. As the varying range of ground elevation of Field 1 and Field 2 is 0.49 m and 0.363 m, respectively, the RMSE showed that the topographic data of the UAV-LiDAR topographic mapping system are of high accuracy.
(2). TIN proved to be the most suitable interpolation model for applications of farmland topographic mapping. Moreover, as the spatial resolution of topographic mapping points is improved, the interpolating error of topographic maps based on the TIN model drops drastically from 0.077 m to 0.046 m. Cut-fill analysis based on the topographic maps of TIN model suggested volumes and areas of soils to be cut or filled, which also provides detailed information on setting the height of the desired ground plane and path planning of farmland levels for future study.
Therefore, we can conclude that based on the TIN interpolation model, the UAV-LiDAR topographic mapping system could be successfully used to collect topographic data with high accuracy, which is instructive for precision farmland levelling.

Author Contributions

Conceptualization, M.D. and A.R.; methodology, M.D.; investigation, M.D. and H.L.; resources, M.D.; data curation, M.D.; writing—original draft preparation, M.D.; writing—review and editing, M.D., H.L. and A.R.; visualization, M.D.; supervision, M.D.; project administration, M.D.; funding acquisition, M.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was jointly funded by the National Key Research and Development Program of China, grant number 2019YFE0125500 and the Key Scientific Research Projects of Colleges and Universities in Henan Province, grant number 20A416001.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Du, M.; Noguchi, N.; Okamoto, H.; Kobayashi, N. Topographic Mapping of Farmland by Integration of Multiple Sensors on Board Low-Altitude Unmanned Aerial System. Int. J. Comput. Syst. Eng. 2017, 11, 1205–1209. [Google Scholar]
  2. Hu, L.; Luo, X.; Lin, C.; Yang, W.; Xu, Y.; Li, Q. Development of 1PJ-4.0 Laser Leveler Installed on a Wheeled Tractor for Paddy Field. Trans. Chin. Soc. Agric. Mach. 2014, 45, 146–151. [Google Scholar]
  3. Rickman, J.F. Manual for Laser Land Levelling; Rice-Wheat Consortium Technical Bulletin Series 5; Rice-Wheat Consortium for Indo Gangetic Plains: New Delhi, India, 2002; Volume 110, p. 24. [Google Scholar]
  4. Liang, W.; Liu, C.; Xiao, R.; Liu, C. Review of Conceptual and Systematic Progress of Precision Irrigation. Int. J. Agric. Biol. Eng. 2021, 14, 20–31. [Google Scholar] [CrossRef]
  5. Khalid, A.; Abdalhaleem, H.; ElKamil, T.; Ahmed, G.K.; Rangaswamy, M.; Fahad, A.; Mohamed, E.; Ahmed, A.; Haroon, E. Impacts of Center Pivot Irrigation System Uniformity on Growth of Potato Crop and Residual Soil Nitrogen. Int. J. Agric. Biol. Eng. 2019, 12, 126–131. [Google Scholar]
  6. Levidow, L.; Zaccaria, D.; Maia, R.; Vivas, E.; Todorovic, M.; Scardigno, A. Improving Water-Efficient Irrigation: Prospects and Difficulties of Innovative Practices. Agric. Water Manag. 2014, 146, 84–94. [Google Scholar] [CrossRef] [Green Version]
  7. Abbaszadeh, A.; Kheiri, A.; Hamzeh, A. Subsurface Topographic Modeling Using Geospatial and Data Driven Algorithm. ISPRS Int. J. Geo-Inf. 2021, 10, 341. [Google Scholar] [CrossRef]
  8. Brandolini, P.; Faccini, F.; Paliaga, G.; Piana, P. Man-made Landforms Survey and Mapping of an Urban Historical Center in a Coastal Mediterranean Environment. Geogr. Fis. Dinam. Quat. 2018, 41, 23–34. [Google Scholar]
  9. Gościewski, D.; Gerus, M. Adjusting the Regular Network of Squares Resolution to the Digital Terrain Model Surface Shape. ISPRS Int. J. Geo-Inf. 2020, 9, 761. [Google Scholar] [CrossRef]
  10. Raeisi, H.; Shayannejad, M.; Soltani, A.; Arab, M.; Eslamian, S.; Amoushahi, M.; Marani, M.; Ostad, K. A Simple Method for Land Grading Computations and its Comparison with Genetic Algorithm (GA) Method. Int. J. Res. Stud. Agric. Sci. 2017, 3, 26–38. [Google Scholar]
  11. Johnson, M.; Ouimet, B. Rediscovering the Lost Archaeological Landscape of Southern New England Using Airborne Light Detection and Ranging (LiDAR). J. Archaeol. Sci. 2014, 43, 9–20. [Google Scholar] [CrossRef]
  12. Stroner, M.; Urban, R.; Seidl, J.; Reindl, T.; Brouček, J. Photogrammetry Using UAV-mounted GNSS RTK: Georeferencing Strategies without GCPs. Remote Sens. 2021, 13, 1336. [Google Scholar] [CrossRef]
  13. Remondino, F. Heritage Recording and 3D Modeling with Photogrammetry and 3D Scanning. Remote Sens. 2011, 3, 1104–1138. [Google Scholar] [CrossRef] [Green Version]
  14. Roncella, R.; Bruno, N.; Diotri, F.; Thoeni, K.; Giacomini, A. Photogrammetric Digital Surface Model Reconstruction in Extreme Low-light Environments. Remote Sens. 2021, 13, 1261. [Google Scholar] [CrossRef]
  15. Ping, H.; Wang, H.; Ma, H.; Du, Y. Mini-Review of Application of IoT Technology in Monitoring Agricultural Products Quality and Safety. Int. J. Agric. Biol. Eng. 2018, 11, 35–45. [Google Scholar] [CrossRef] [Green Version]
  16. Ham, Y.; Han, K.; Lin, J.; Golparvar, M. Visual Monitoring of Civil Infrastructure Systems via Camera-Equipped Unmanned Aerial Vehicles (UAVs): A Review of Related Works. Vis. Eng. 2016, 4, 1. [Google Scholar] [CrossRef]
  17. Du, M.; Liu, Y.; Ji, J.; Jin, X.; Zhou, H.; Liu, K. Farmland Topographic Mapping Based on UAV and LiDAR Technology. Trans. Chin. Soc. Agric. Eng. 2020, 36, 60–67. [Google Scholar]
  18. Yang, J.; Li, X.; Luo, L.; Zhao, L.; Wei, J.; Ma, T. New Supplementary Photography Methods after the Anomalous of Ground Control Points in UAV Structure-from-Motion Photogrammetry. Drones 2022, 6, 105. [Google Scholar] [CrossRef]
  19. Resop, P.; Hession, C. Terrestrial Laser Scanning for Monitoring Streambank Retreat: Comparison with Traditional Surveying Techniques. J. Hydraul. Eng. 2010, 136, 794–798. [Google Scholar] [CrossRef]
  20. Keay, S.; Earl, P.; Hay, S.; Kay, s.; Ogden, J.; Strutt, K. The Role of Integrated Geophysical Survey Methods in the Assessment of Archaeological Landscapes: The case of Portus. Archaeol. Prospect. 2009, 16, 154–166. [Google Scholar]
  21. Guo, Z.; Zheng, J.; Cui, Y. Surveying and Setting-out Based on Two-Point Reference Line with Total Station. Bull. Surv. Mapp. 2004, 8, 22. [Google Scholar]
  22. Dardanelli, G.; Maltese, A.; Pipitone, C.; Pisciotta, A.; Brutto, M. NRTK, PPP or Static, That Is the Question. Testing Different Positioning Solutions for GNSS Survey. Remote Sens. 2021, 13, 1406. [Google Scholar] [CrossRef]
  23. Keay, J.; Parcak, H.; Strutt, D. High Resolution Space and Ground-Based Remote Sensing and Implications for Landscape Archaeology: The Case from Portus, Italy. JAS 2014, 52, 277–292. [Google Scholar] [CrossRef]
  24. Corsini, A.; Castagnetti, C.; Bertacchini, E.; Riccardo, R. Integrating Airborne and Multi-Temporal Long-Range Terrestrial Laser Scanning with Total Station Measurements for Mapping and Monitoring a Compound Slow Moving Rock Slide. Earth Surf. Process. Landf. 2013, 38, 1330–1338. [Google Scholar] [CrossRef]
  25. Rodríguez, P.; Jiménez, B.; Muñoz, L.; Arias, P.; Gonzalez, D. Mobile LiDAR System: New Possibilities for the Documentation and Dissemination of Large Cultural Heritage Sites. Remote Sens. 2017, 9, 189. [Google Scholar] [CrossRef] [Green Version]
  26. Terrone, M.; Piana, P.; Paliaga, G.; D’Orazi, M.; Faccini, F. Coupling Historical Maps and LiDAR Data to Identify Man-Made Landforms in Urban Areas. ISPRS Int. J. Geo-Inf. 2021, 10, 349. [Google Scholar] [CrossRef]
  27. Hirano, A.; Welch, R.; Lang, H. Mapping from ASTER Stereo Image Data: DEM Validation and Accuracy Assessment. ISPRS J. Photogramm. Remote Sens. 2003, 57, 356–370. [Google Scholar] [CrossRef]
  28. Doneus, M. Openness as Visualization Technique for Interpretative Mapping of Airborne Lidar Derived Digital Terrain Models. Remote Sens. 2013, 5, 6427–6442. [Google Scholar] [CrossRef] [Green Version]
  29. Harwin, S.; Lucieer, A. Assessing the Accuracy of Georeferenced Point Clouds Produced via Multi-View Stereopsis from Un-manned Aerial Vehicle (UAV) Imagery. Remote Sens. 2012, 4, 1573–1599. [Google Scholar] [CrossRef]
  30. Alshawabkeh, Y.; Baik, A.; Miky, Y. Integration of Laser Scanner and Photogrammetry for Heritage BIM Enhancement. ISPRS Int. J. Geo-Inf. 2021, 10, 316. [Google Scholar] [CrossRef]
  31. Lovitt, J.M.; McDermid, J. Assessing the Value of UAV Photogrammetry for Characterizing Terrain in Complex Peatlands. Remote Sens. 2017, 9, 715. [Google Scholar] [CrossRef] [Green Version]
  32. Jiménez, I.; Ojeda, W.; Marcial, J.; Enciso, J. Digital Terrain Models Generated with Low-Cost UAV Photogrammetry: Methodology and Accuracy. ISPRS Int. J. Geo-Inf. 2021, 10, 285. [Google Scholar] [CrossRef]
  33. Cho, J.-W.; Lee, J.-K.; Park, J. Large-Scale Earthwork Progress Digitalization Practices Using Series of 3D Models Generated from UAS Images. Drones 2021, 5, 147. [Google Scholar] [CrossRef]
  34. Crawford, B.; Swanson, E.; Schultz-Fellenz, E.; Collins, A.; Dann, J.; Lathrop, E.; Milazzo, D. A New Method for High Resolution Surface Change Detection: Data Collection and Validation of Measurements from UAS at the Nevada National Security Site, Nevada, USA. Drones 2021, 5, 25. [Google Scholar] [CrossRef]
  35. Angelo, G.; Piersanti, M.; Pignalberi, A.; Coco, I.; De Michelis, P.; Tozzi, R.; Pezzopane, M.; Alfonsi, L.; Cilliers, P.; Ubertini, P. Investigation of the Physical Processes Involved in GNSS Amplitude Scintillations at High Latitude: A Case Study. Remote Sens. 2021, 13, 2493. [Google Scholar]
  36. Liu, X.; Chen, H.; Jiang, W.; Xi, R. Modeling and Assessment of GNSS/Galileo/BDS Precise Point Positioning with Ambiguity Resolution. Remote Sens. 2019, 11, 2693. [Google Scholar] [CrossRef] [Green Version]
  37. Bewley, H.; Crutchley, S.; Shell, C. New Light on an Ancient Landscape: LiDAR Survey in the Stonehenge World Heritage Site. Antiquity 2005, 79, 636–647. [Google Scholar] [CrossRef]
  38. Devereux, B.J.; Amable, G.S.; Crow, P. Visualisation of LiDAR Terrain Models for Archaeological Feature Detection. Antiquity 2008, 82, 470–479. [Google Scholar] [CrossRef]
  39. Bennett, R.; Welham, K.; Hill, R.; Ford, A. A Comparison of Visualization Techniques for Models Created from Airborne Laser Scanned Data. Archaeol. Prospect. 2012, 19, 41–48. [Google Scholar] [CrossRef]
  40. Pricope, N.G.; Minei, A.; Halls, J.N.; Chen, C.; Wang, Y. UAS Hyperspatial LiDAR Data Performance in Delineation and Classification across a Gradient of Wetland Types. Drones 2022, 6, 268. [Google Scholar] [CrossRef]
  41. Shukla, D.; Komerath, N. Multirotor drone aerodynamic interaction investigation. Drones 2018, 2, 43. [Google Scholar] [CrossRef] [Green Version]
  42. Artale, V.; Milazzo, L.; Ricciardello, A. Mathematical Modeling of Hexacopter. Appl. Math. Sci. 2013, 7, 4805–4811. [Google Scholar] [CrossRef]
  43. Nguyen, D.D.; Rohacs, J.; Rohacs, D. Autonomous Flight Trajectory Control System for Drones in Smart City Traffic Management. ISPRS Int. J. Geo-Inf. 2021, 10, 338. [Google Scholar] [CrossRef]
  44. Zhang, J.; Ren, L.; Deng, H.; Ma, M.; Zhong, X.; Wen, P. Measurement of Unmanned Aerial Vehicle Attitude Angles Based on a Single Captured Image. Sensors 2018, 18, 2655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Chiang, K.W.; Duong, T.T.; Liao, J.K.; Lai, Y.C.; Chang, C.C.; Cai, J.M.; Huang, S.C. On-Line Smoothing for an Integrated Navigation System with Low-Cost MEMS Inertial Sensors. Sensors 2012, 12, 17372. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Zhang, J.; Yuan, H. Analysis of Unmanned Aerial Vehicle Navigation and Height Control System Based on GNSS. Syst. Eng. Electron. Technol. 2010, 21, 643–649. [Google Scholar] [CrossRef]
  47. Sun, Y.; Zheng, W.; Li, Z.; Zhou, Z. Improved the Accuracy of Seafloor Topography from Altimetry-Derived Gravity by the Topography Constraint Factor Weight Optimization Method. Remote Sens. 2021, 13, 2277. [Google Scholar] [CrossRef]
  48. Huang, P. Triangular Irregular Network Generation and Topographical Modeling. Comput. Ind. 1989, 12, 203–213. [Google Scholar] [CrossRef]
  49. Tsai, J. Delaunay Triangulations in TIN Creation: An Overview and a Linear-Time Algorithm. Int. J. Geo-Graph. Inf. Sci. 1993, 7, 501–524. [Google Scholar] [CrossRef]
  50. Watson, F.; Philip, G.M. A Refinement of Inverse Distance Weighted Interpolation. Geoprocessing 1985, 2, 315–327. [Google Scholar]
Figure 1. Flight path of UAV-LiDAR topographic surveying experiment and locations of accuracy validating points. Note: The black lines, green dots, and red dots indicate the flight trajectories of the UAV, the sampling points of handheld PPK-GNSS for evaluating the accuracy of ground surveying points of the UAV-LiDAR topographic surveying system, and the sampling points of handheld PPK-GNSS for evaluating the accuracy of topographic maps based on interpolation models.
Figure 1. Flight path of UAV-LiDAR topographic surveying experiment and locations of accuracy validating points. Note: The black lines, green dots, and red dots indicate the flight trajectories of the UAV, the sampling points of handheld PPK-GNSS for evaluating the accuracy of ground surveying points of the UAV-LiDAR topographic surveying system, and the sampling points of handheld PPK-GNSS for evaluating the accuracy of topographic maps based on interpolation models.
Drones 06 00403 g001
Figure 2. Equipment of the topographic surveying experiment. (a) UAV platform; (b) Base receiver of PPK-GNSS module.
Figure 2. Equipment of the topographic surveying experiment. (a) UAV platform; (b) Base receiver of PPK-GNSS module.
Drones 06 00403 g002
Figure 3. Altitudes of PPK-GNSS positioning data on board the UAV platform. (a) Field 1; (b) Field 2.
Figure 3. Altitudes of PPK-GNSS positioning data on board the UAV platform. (a) Field 1; (b) Field 2.
Drones 06 00403 g003
Figure 4. Scheme for rectifying LiDAR distance measurements using UAV attitude information.
Figure 4. Scheme for rectifying LiDAR distance measurements using UAV attitude information.
Drones 06 00403 g004
Figure 5. Nadir distance between the UAV-LiDAR system and ground surface. (a) Field 1; (b) Field 2.
Figure 5. Nadir distance between the UAV-LiDAR system and ground surface. (a) Field 1; (b) Field 2.
Drones 06 00403 g005
Figure 6. Scheme for calculating elevations of ground surveying points.
Figure 6. Scheme for calculating elevations of ground surveying points.
Drones 06 00403 g006
Figure 7. Positioning accuracy of standalone GNSS and PPK-GNSS. (a) Horizontal positioning accuracy of standalone GNSS; (b) Horizontal positioning accuracy of PPK-GNSS; (c) Vertical positioning accuracy of PPK-GNSS.
Figure 7. Positioning accuracy of standalone GNSS and PPK-GNSS. (a) Horizontal positioning accuracy of standalone GNSS; (b) Horizontal positioning accuracy of PPK-GNSS; (c) Vertical positioning accuracy of PPK-GNSS.
Drones 06 00403 g007
Figure 8. Graduated symbols represent elevations of ground surveying points. (a) Field 1; (b) Field 1.
Figure 8. Graduated symbols represent elevations of ground surveying points. (a) Field 1; (b) Field 1.
Drones 06 00403 g008
Figure 9. Topographic maps generated using the TIN interpolating model. (a) Field 1; (b) Field 2.
Figure 9. Topographic maps generated using the TIN interpolating model. (a) Field 1; (b) Field 2.
Drones 06 00403 g009
Figure 10. Topographic maps generated using the IDW interpolation model. (a) Field 1; (b) Field 2.
Figure 10. Topographic maps generated using the IDW interpolation model. (a) Field 1; (b) Field 2.
Drones 06 00403 g010
Figure 11. Topographic maps generated using the natural neighbor interpolation model. (a) Field 1; (b) Field 2.
Figure 11. Topographic maps generated using the natural neighbor interpolation model. (a) Field 1; (b) Field 2.
Drones 06 00403 g011
Figure 12. Topographic maps generated using the Kriging interpolation model. (a) Field 1; (b) Field 2.
Figure 12. Topographic maps generated using the Kriging interpolation model. (a) Field 1; (b) Field 2.
Drones 06 00403 g012
Figure 13. Cut-fill analysis based on topographic maps of the TIN model. (a) Field 1; (b) Field 2.
Figure 13. Cut-fill analysis based on topographic maps of the TIN model. (a) Field 1; (b) Field 2.
Drones 06 00403 g013
Table 1. Parameters of the equipment of the UAV-LiDAR topographic mapping system.
Table 1. Parameters of the equipment of the UAV-LiDAR topographic mapping system.
EquipmentItemsValues
UAV platformOverall size (mm) 1290 × 650
Motor typeTTA6215
Standard take-off weight (kg)23
Flight controllerPixhawk
Flight altitude (m)30
GNSS moduleGNSS receiversU-blox NEO-M8T
Processing unitIntel Edison
Overall size (mm)45.5 × 27 × 9.2
Weight (g)14
Supply voltage (V)4.75 to 5.5 DC
Power consumption (W)<1
SignalsGNSS, GLONASS, BDS
Out frequency (Hz)10
AntennaTallysman TW4721
LiDAROverall size (mm)136 × 104 × 57
Weight (g)800
Measuring resolution (mm)1
Output frequency (Hz)2000
Measuring range (m)300
Beam divergence (milliradians)1.7
Measuring laser wavelength (nm)905
Supply voltage (V)10 to 30 DC
Power consumption (W)<5
Table 2. 3D coordinates of topographic mapping points measured by UAV-LiDAR topographic mapping system and PPK-GNSS module.
Table 2. 3D coordinates of topographic mapping points measured by UAV-LiDAR topographic mapping system and PPK-GNSS module.
Field NO.Point NO.Easting/mNorthing/mGround Elevation/m
UAV-LiDAR SystemHandheld PPK-GNSS
Field 11312,217.073,848,795.9363.65863.616
2312,217.803,848,801.5263.62863.593
3312,219.743,848,803.2463.6463.655
4312,223.893,848,806.0763.65763.679
5312,226.453,848,809.2063.64963.607
6312,229.303,848,811.5163.40463.476
7312,229.853,848,790.9363.60563.565
8312,233.013,848,793.7363.6763.631
9312,234.603,848,796.0563.51463.568
10312,236.793,848,798.7263.63163.664
11312,239.103,848,802.9963.57263.529
12312,241.293,848,805.1863.68563.662
Field 21312,243.613,848,872.2663.42063.369
2312,247.923,848,874.0663.52263.508
3312,251.313,848,877.7763.52663.509
4312,254.223,848,881.9663.53163.596
5312,257.823,848,887.0563.54963.524
6312,263.323,848,890.2663.53263.584
7312,273.223,848,883.3563.60963.558
8312,270.823,848,880.4663.63563.696
9312,265.423,848,876.8663.62963.542
10312,262.923,848,873.2763.64163.657
11312,260.313,848,869.3663.54263.681
12312,253.823,848,864.6663.53563.578
Table 3. 3D coordinates of sampled points acquired using a handheld PPK-GNSS module and different interpolation models.
Table 3. 3D coordinates of sampled points acquired using a handheld PPK-GNSS module and different interpolation models.
Field NO.Point NO.Easting/mNorthing/mGround Elevation/m
PPK-GNSSTINIDWNatural NeighborKriging
Field 11312,227.813,848,810.6263.46163.48963.47663.50163.497
2312,225.173,848,807.2763.58763.65763.65563.65463.654
3312,222.053,848,804.3963.69363.64563.63663.63263.631
4312,240.043,848,804.3963.61863.63763.64463.64463.644
5312,218.933,848,801.7563.58563.63663.64263.64163.643
6312,238.123,848,800.0763.65663.61963.62963.61763.627
7312,217.253,848,798.1563.57863.63163.64263.63863.642
8312,235.483,848,797.4363.46263.57063.57963.57563.577
9312,234.523,848,795.0363.70463.59063.57363.58563.586
10312,231.643,848,792.6363.77963.63963.63663.63363.634
Field 21312,259.813,848,887.7963.58963.54663.54363.54663.546
2312,257.093,848,884.1063.52063.57263.58163.57863.579
3312,272.053,848,881.3863.67663.62663.62363.62563.624
4312,253.593,848,879.4363.55463.50863.50463.50263.502
5312,268.753,848,878.8563.68363.63663.62663.63463.633
6312,250.093,848,876.7163.53863.48163.47463.48563.483
7312,264.473,848,874.1963.67063.65163.65163.65263.652
8312,246.783,848,873.4163.49563.47163.46263.47163.469
9312,261.363,848,871.6663.51963.57863.59463.58663.588
10312,256.893,848,867.3863.59363.54963.54663.54963.551
Table 4. RMSE between ground elevation data of a handheld PPK-GNSS module and topographic maps of different interpolating models.
Table 4. RMSE between ground elevation data of a handheld PPK-GNSS module and topographic maps of different interpolating models.
TINIDWNatural NeighborKriging
Field 10.077 m0.083 m0.082 m0.082 m
Field 20.046 m0.053 m0.048 m0.049 m
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Du, M.; Li, H.; Roshanianfard, A. Design and Experimental Study on an Innovative UAV-LiDAR Topographic Mapping System for Precision Land Levelling. Drones 2022, 6, 403. https://doi.org/10.3390/drones6120403

AMA Style

Du M, Li H, Roshanianfard A. Design and Experimental Study on an Innovative UAV-LiDAR Topographic Mapping System for Precision Land Levelling. Drones. 2022; 6(12):403. https://doi.org/10.3390/drones6120403

Chicago/Turabian Style

Du, Mengmeng, Hanyuan Li, and Ali Roshanianfard. 2022. "Design and Experimental Study on an Innovative UAV-LiDAR Topographic Mapping System for Precision Land Levelling" Drones 6, no. 12: 403. https://doi.org/10.3390/drones6120403

APA Style

Du, M., Li, H., & Roshanianfard, A. (2022). Design and Experimental Study on an Innovative UAV-LiDAR Topographic Mapping System for Precision Land Levelling. Drones, 6(12), 403. https://doi.org/10.3390/drones6120403

Article Metrics

Back to TopTop