1. Introduction
Nowadays, unmanned aerial vehicles (UAVs) are frequently preferred for generating photogrammetric outcomes. One of the reasons for this is their ability to quickly generate topographic models of large areas with relatively high spatial accuracy [
1,
2]. Their ability to quickly generate topographic models of large areas with relatively high spatial accuracy is one of the reasons why they are preferred [
1,
2]. Topographic models can be categorized as the digital elevation model (DEM), digital terrain model (DTM) and digital surface model (DSM). In addition to these, outcomes prepared by performing classical geodetic surveys (3D drawing, mapping, etc.) can also be generated faster and easier with UAVs [
3,
4,
5]. In addition to the advantages of UAVs, such as the speed, convenience, and the ability to access places that are not accessible to humans from the air, they also have disadvantages, such as being affected by atmospheric conditions and being unable to generate accuracy below the cm level. Improving the positional accuracy of the outcomes generated by UAVs is also an important research topic.
The integrated global navigation satellite system (GNSS) receivers of UAVs make it impossible to obtain a spatially accurate outcome, neither horizontally nor vertically [
6,
7]. The most commonly used technique for bundle block adjustment (BBA) is seen in the outcomes generated by the photogrammetric method. In order for the adjusted block accuracy to be high, points with known geodetic coordinates are needed. Using ground control points (GCPs) for the BBA is traditional. However, the number of GCPs to be established limits the temporality of the process in UAV operations. An important research topic is whether precise point positioning (PPP) can be an alternative for UAVs, mainly when geodetic position measurements cannot be performed with real-time kinematic (RTK), network RTK, and post-process kinematic (PPK) [
8]. Polat and Uysal [
9] compared airborne LIDAR and UAV-based DEM and found that the vertical accuracy of UAV-based DEM was ±15.7 cm. Martínez-Carricondo et al. [
10] investigated the accuracy of UAV-photogrammetric mapping based on the GCP distribution. They reported that the best accuracy was achieved by placing GCPs at the edge of the study area. Manfreda et al. [
11] suggested that the internal accuracy of the UAV-based model can reach 0.2 cm horizontally and 4 cm vertically. Yu et al. [
12] analyzed the model accuracy and claimed that it is approximately 0.90 m. Zimmerman et al. [
13] analyzed the effect of the UAV height and GCPs on the model accuracy in a selected application area along a coastline and achieved accuracies ranging from 3 to 18 cm. Famiglietti et al. [
14], Žabota and Kobal [
15], Liu et al. [
16], Martínez-Carricondo et al. [
17] and Hayamizu and Nakata [
18] investigated the accuracy of models generated based on the RTK/PPK techniques. Recent studies show that the relative model accuracy obtained in RTK studies is between 0.5 m and 2.5 m [
19,
20], while the model accuracy (relative) in PPK studies is between 2.45 cm and 3.75 cm [
19,
21]. For the models obtained using GCPs with the geodetic surveying method, the horizontal positional accuracy of 8–11 cm and vertical positional accuracy of 9–74 cm was obtained in the measure results with network RTK [
19,
22,
23] (
Table 1).
The scientific studies mentioned above show that the RTK/PPK method is mainly preferred, and there are limited studies exploring PPP-AR as an alternative method for UAV photogrammetry, highlighting the need for further research in this area. This study proposes a new approach by investigating the contribution of the PPP technique to the model accuracy in terms of the positioning of GCPs and this is the main difference from the existing literature. Another novel aspect of this study is the determination of the optimal GNSS session duration by testing the PPP technique at different observation times. The PPP technique has the advantages of being implemented with free online or open-source software, requiring only a single GNSS receiver and not requiring a Global System for Mobile Communications (GSM) connection when surveying in the field. The only disadvantage is that it requires a long convergence time to achieve centimeter position accuracy [
24,
25,
26,
27]. Since PPK is available on a small number of UAVs and PPK decoding is performed with paid software, network RTK does not work everywhere due to the need for a GSM network, and users may not have enough GNSS receivers for classical RTK; the PPP technique can be proposed as an alternative. Whether this proposed technique is suitable for UAVs has received limited attention in scientific studies. This paper addresses this gap by demonstrating the viability and benefits of PPP-AR for GCP positioning.
The aim of this study is to suggest a new approach using PPP techniques for GCP coordinate determination for UAV-based mapping. Within this context, different photogrammetric outcomes are generated. The study looks at the accuracy of the orthomosaic and DEM created without using any GCPs, with four GCPs and eight GCPs. The accuracy is assessed by comparing the horizontal and vertical coordinates obtained from 24 h fixed GNSS sessions of six calibration pillars in the field and the horizontal length differences obtained by electronic distance measuring devices (EDMs). Statistically, Bartlett’s test is applied to determine the accuracy of the results. The results show that the coordinates obtained with the two-hour session did not show a relatively significant difference in positional accuracy compared to the coordinates obtained with the half-hour session. In this context, the PPP technique could be used to establish the BBA of the GCPs for UAVs in large-scale map generation.
3. Application Results and Discussion
In this study, a 115 m × 820 m (width × length) study area was selected within the campus area of Konya Selçuk University (
Figure 4). Within the study area, there were six pillar points whose 3D position accuracies were determined with high precision. These pillars constituted the control base used to determine the calibration parameters of the EDMs. Eight GCPs were installed to cover the application area (
Figure 4) homogeneously. Measurement of GCPs has been performed with the Javad Triumph1 GNSS receiver for two hours, with 30 s recording intervals.
Before starting the measurements, the GCPs were installed in red and white colors so that they could be easily seen from the single image, as shown in
Figure 5a. The GNSS receiver was fixed and leveled on the GCPs with a tripod, and static observations were performed.
Figure 5b shows one of the pillars whose coordinates are obtained as a result of the long-term session. By processing the long-term GNSS observations previously made on the batteries, the coordinates are estimated by geodetic network adjustment with the relative GNSS technique. The RMSEs of the estimated coordinates are in the order of mm.
The GNSS measurements performed at the GCPs are processed with the CSRS-PPP software version 3.0 at observation times of 30 min, 60 min and two hours. These durations are chosen to compare the impact of the session duration on the accuracy of the generated topographic models. The outcomes for the study area are generated in three different scenarios: with no GCPs, with four GCPs (3, 9, 11, 21) and with all the GCPs (
Figure 6).
The horizontal and vertical RMSEs for each scenario are calculated using Equation (8) (
Table 6).
The model accuracies based on the point locations are shown in
Table 6. The horizontal accuracy of the model generated with no GCPs is 102.26–105.40 cm, and the vertical accuracy is 335.57–341.46 cm. The model generated with four GCPs has a horizontal accuracy of 1.07–1.32 cm and a vertical accuracy of 14.20–14.86 cm. When all the GCPs are used, the accuracies decrease to 16.39–16.52 cm horizontally and 24.25–26.42 cm vertically. According to
Figure 6 and
Table 6, the three scenarios’ accuracies are significantly different, and the model with the four GCPs selected at the model corners with the smallest RMSE provides the best result. The main difference is due to the presence of grossly inaccurate points in all the GCPs, which have a negative impact on model accuracy. The four selected GCPs are both homogeneously located points within the block and have the lowest positional errors. These choices also affect the overall model accuracy. Due to the fact that the coordinates of the images, camera calibration values, and local coordinates values are used in the BBA, the error is distributed to the entire model (in case of even one or two grossly inaccurate coordinate values) [
35,
36,
37,
38].
Bartlett’s test is used to statistically compare the RMSEs of the models generated with GCPs at different session durations. Three horizontal and three vertical RMSEs are compared for all the models generated using 0-, 4- and 8-GCPs. Using Equations (9)–(12),
test statistics are calculated and shown in
Table 7.
The values in
Table 7 show the Bartlett’s test statistics of the RMSEs calculated by comparing the coordinates taken from the generated models with the PPP-2 h, PPP-1 h and PPP-30 m coordinates. The calculated test statistics are compared with the critical value in the table (5.9914) and are not significant. The conclusion to be drawn from this is that the coordinate accuracies of the models at three different session durations are consistent with each other. In other words, when the PPP technique is used in the GCPs, there is no significant difference between observing for 30 min and observing for 2 h. However, when the number of GCPs changes, the accuracy of the models generated is significantly different, and the best accuracy is achieved with 4 GCPs. This is due to the fact that there are gross inaccuracies in all the GCPs that affect the model accuracy. To analyze the model accuracies in more detail, the control–baseline lengths are calculated using the pillar coordinates from different models and compared with the known lengths.
When
Figure 7 is evaluated, it is seen that the length differences in the model generated without a GCP vary between 185.2 and 16.6 cm. Moreover, 53% of the differences are greater than 91.6 cm. Differences 1–3, 1–4, 1–5, 1–6, 2–4, 2–5, 2–6, 3–5 and 3–6 are relatively large. The magnitude of these differences is thought to be due to the coarse error in the coordinates of the pillars numbered 4, 5 and 6 taken from the model. The differences of 1–2, 2–3, 3–4, 4–5, 4–6 and 5–6 are smaller than 59.8 cm, which are pretty good values for this model, with an average accuracy of 1 m. If
Table 6 and
Figure 7 are evaluated together, it can be said that while the horizontal accuracy of the model is approximately 1 m, almost half of the differences are 1 m and more considerable. This is due to the low accuracy of the UAV’s internal GNSS receiver and the nonmetric nature of the UAV’s integrated camera.
Figure 8 shows that the differences for the four GCPs range from −12.2 to 22.6 cm. In addition, 93% of the differences are 10 cm or less. Differences 1–5, 1–6, 2–5, 2–6, 3–5, 3–6, 4–5, 4–6, 5–6 are relatively large. All the differences between pillars 5 and 6 are large due to the fact that these pillars are on the edge of the block. The differences of 2–3, 2–4 and 3–4 in the middle of the selected GCPs are 0.2–3.4 cm. These findings indicate that when a length in the middle of the GCPs is measured from the model, accurate results comparable to ground measurements can be obtained. If a low-cost UAV and a GNSS receiver are used to generate the model by selecting four GCPs homogeneously distributed in the region, differences below 10 cm are found to be achievable.
For the model generated using all the GCPs, the differences between the known lengths and the measured lengths ranged between −25.8 and 20.2 cm for 30 min, −25.1 and 18.6 cm for 60 min, and −25.1 and 19.1 cm for 2 h (
Figure 9). The differences of 1–2, 1–4, 1–6, 2–5, 3–4, 4–5 and 5–6 are greater than 10 cm in 30, 60 and 120 min. The main reason for this is that pillars 1, 5 and 6 are located in a region outside the GCPs, and the other reason is that coarse errors are spread over the modelled coordinates of pillars 3 and 4. The fact that 53.3% of the differences are smaller than 10 cm shows that the model generated is more accurate than the model with no GCPs and more inaccurate than the model generated with four GCPs. The RMSEs of the differences calculated using all the GCPs for the 30 min measurements are ±13.3 cm, for the 60 min measurements ±13.0 cm, and for the 2 h measurements ±12.8 cm. It can be seen that the RMSEs calculated for different session durations using all the RMSEs are statistically consistent. In addition, these RMSEs show that the models generated can provide outstanding results at three different PPP session durations. With a low-cost UAV with a nonmetric camera, the average ±13.0 cm accuracy shows that topographic models can be obtained from which length measurements can be performed. In addition to the horizontal distances, the pillar heights and known heights from different models are also compared in this study. The differences are presented in
Figure 10,
Figure 11 and
Figure 12.
When
Figure 10 is evaluated, it is seen that the differences between the heights taken from the model generated without a GCP and the known heights vary between −532.5 and394.7 cm. Differences 1–5, 2–4 and 5–6 are less than 76 cm, while all the other differences are greater than 100 cm. This is due to both the inability of the internal GNSS receiver to provide sufficient accuracy and the lower overall accuracy of the GNSS method in height compared to the horizontal. When the vertical accuracies in
Table 6 and
Figure 10 are considered together, 47% of the differences are around 300 cm and the vertical RMSEs are ±335.57–±341.46. These values show that an accurate height difference cannot be obtained from the model generated without using any GCP.
Figure 11 shows that the differences for four GCPs range from −270.7 to 329.5 cm. Differences 1–2, 1–3, 1–4, 1–5, 2–3, 2–4, 2–6, 3–5, 3–6, 4–5, 4–6 are larger than ten dm. The remaining ones account for 27% of all the differences and are between 4.4 and 7.1 dm. These values show that the height differences determined from the model are more accurate when four GCPs are used than when no GCPs are used. However, it cannot provide precise height information as it is in the order of decimeters.
Using all the control points, the differences between the known heights and measured heights range between −39.29 and 32.94 cm at 30, 60 and 120 min. The differences between 1–2, 1–3, 1–4, 2–4, 2–5, 2–6, 3–4, 3–5, 3–6, 4–5 and 4–6 are greater than 10 cm at 30, 60 and 120 min. The differences of 1–5, 1–6, 2–3 and 5–6 are between 1.9 and 6.3 cm. Precise height differences at the cm level could be obtained for only 26% of the differences. This is mainly due to the fact that the BBA propagates the errors of coarse error points to other points. Nevertheless, the most accurate results in terms of the height differences are obtained in the models generated using all the GCPs. While the differences in
Figure 10 and
Figure 11 are in the decimeter range, differences in the centimeter range are achieved in
Figure 12. The RMSEs of the differences calculated using all the GCPs are ±20.5 cm for 30 min, ±20.5 cm for 60 min and ±20.77 cm for 2 h. It is seen that the vertical RMSEs calculated for different session durations using all the RMSEs are statistically consistent (
Figure 12). UAVs are preferred for generating large-scale maps because they provide fast results, and it is usual to obtain an accuracy of dm and above in the outcomes generated without the use of a GCP. As a result, integrated GNSS systems determine their position as ordinary differential GNSS (D-GNSS) receivers. Studies in the literature with D-GNSS also show that the positioning accuracy can be in the order of dm [
44,
45,
46,
47]. The errors obtained as a result of the accuracy analysis allow the generation of large-scale maps of the outcomes planned to be obtained and to obtain relative accuracy that can be considered sufficient.
In addition, in order to make the accuracy analysis more understandable, the differences between the coordinates taken from different models and the control-based lengths and known lengths are shown in box plot graphs.
When
Figure 13 is evaluated, it is seen that the differences between the lengths taken from the models generated without any GCP and the known lengths are typically distributed. There are no outliers, and the average is 92.5 cm. The maximum value of the differences is 185.2 cm, and the minimum is 16.6 cm. Seven of the differences are above the mean, and eight are below the mean. This shows that statistical calculations and inferences based on the differences can be made reliably.
Figure 14 shows that the differences for the four GCPs are again generally distributed, with a mean of 1.8 cm. The maximum value of the differences is 22.63 cm, and the minimum is −12.16 cm. Seven of the differences are above the mean, and eight are below the mean. This shows that statistical calculations and inferences based on the differences can be performed reliably.
It is seen that the differences between the known lengths and the lengths taken from the models generated with 30, 60 and 120 min measurements using all the GCPs have a normal distribution, there are no outliers, and the mean of the differences is 1.2 cm for the 30 min measurements, 0.9 cm for the 60 min measurements, and 0.9 cm for the 2 h measurements. (
Figure 15). In the 30 min model, the maximum difference is 20.2 cm, and the minimum difference is −25.8 cm; in the 60 min model, the maximum difference is 18.6 cm, and the minimum difference is −25.1 cm; and in the 120 min model, the maximum difference is 19.1 cm, and the minimum difference is −25.1 cm. In all three scenarios, the ranges of variation of the differences are almost equal to each other. These findings show that reliable statistical inferences can be drawn using these differences.
The differences between the control–baseline point heights from different models and the known heights are also shown in the box plot graphs in
Figure 16,
Figure 17 and
Figure 18.
When
Figure 16 is evaluated, it is seen that the differences between the heights taken from the model generated without any GCP and the known heights are normally distributed. There are no outliers, and the mean is −11.8 dm. The maximum difference is 39.5 dm, and the minimum is −53.2 dm. In the box plot, seven of the values are above the mean, and eight are below the mean.
The four control points’ differences are typically distributed again, with a mean of 0.2 dm (
Figure 17). The maximum difference is 32.9 dm, and the minimum is −27.1 dm. Of the differences in the box plot, seven are above the mean, and eight are below the mean. The differences are normally distributed, and there are no outliers.
It is seen that the differences between the heights taken from the models created using all the GCPs and the known heights have a normal distribution with 30, 60 and 120 min measurements. There are no outliers, and the average of the differences is −3.4 cm for the 30 min measurements, −3.4 cm for the 60 min measurements, and −3.6 cm for the 2 h measurements (
Figure 18). When
Figure 16,
Figure 17 and
Figure 18 are evaluated together, it is seen that reliable statistical parameters related to height can be calculated from these data and reliable inferences can be drawn.
The photogrammetric BBA is not a fully automated process. Many factors directly affect the model accuracy, such as marking the GCP locations from the image (performed manually by the operator) and ensuring that the selected GCP location accuracy is appropriate and that the GCPs are homogeneously distributed within the block when installed in the field [
10]. In addition, the location information from which the GCP coordinates are taken (such as how many GNSS satellites are connected and at what angles), the date the images were taken (to determine the solar activity known as the Kp index), the type of UAV used, the UAV camera, UAV-integrated systems and UAV flight planning also have an impact on the model accuracy [
13]. For this reason, it is not possible to talk about absolute accuracy like geodetic methods. In previous studies, the vertical accuracy value found in the comparison of the DEM generated from LIDAR and UAV images was ±15.7 cm [
9], and in another study, it was determined that it could reach 0.2 cm horizontally and 4 cm vertically [
11]. Another study determined that the model accuracy varied between 0.9 and 10 cm as a result of the analysis [
12]. In this study, it is assumed that the horizontal positional accuracy of the 3D elevation models or ortho data generated in such a way that the RMSE value does not exceed ±10 cm and the vertical positional accuracy does not exceed ±20 cm are acceptable values for researchers and practitioners. As a result of the determination of the locations of the GCPs established in the field with PPP methods with different session durations and the accuracy analysis of the models generated with UAVs, it is seen that the measurement time of 30 min with the PPP-AR method made it possible to achieve the specified accuracies.