Application the SBAS/EGNOS Corrections in UAV Positioning

: The paper presents a new concept of determining the resultant position of a UAV (Unmanned Aerial Vehicle) based on individual SBAS (Satellite-Based Augmentation System) determinations from all available EGNOS (European Geostationary Navigation Overlay Service) satellites for the SPP (Single Point Positioning) code method. To achieve this, the authors propose a weighted mean model to integrate EGNOS data. The weighted model was based on the inverse of the square of the mean position error along the component axes of the BLh ellipsoidal frame. The calculations included navigation data from the EGNOS S123, S126, S136 satellites. In turn, the resultant UAV position model was determined using the Scilab v.6.0.0 software. Based on the proposed computational strategy, the mean values of the UAV BLh coordinates’ standard deviation were better than 0.2 m (e.g., 0.0000018 ◦ = 0.01” in angular measurement). Additionally, the numerical solution used made it possible to increase the UAV’s position accuracy by about 29% for Latitude, 46% for Longitude and 72% for ellipsoidal height compared to the standard SPP positioning in the GPS receiver. It is also worth noting that the standard deviation of the UAV position calculated from the weighted mean model improved by about 21 ÷ 50% compared to the arithmetic mean model’s solution. It can be concluded that the proposed research method allows for a signiﬁcant improvement in the accuracy of UAV positioning with the use of EGNOS augmentation systems.


Introduction
The intensive development of unmanned aerial vehicle technologies has led to them being used in many modern technology areas. One such application is low-altitude imaging for photogrammetric and remote sensing studies. Relatively low costs of procurement and operation made this technology extremely popular [1,2]. However, still up to this day, there are UAVs (Unmanned Aerial Vehicles) available on the market that are equipped with single-frequency GPS (Global Positioning System) receivers, which allow the UAV's position to be determined with an accuracy of only just better than 10 m [3,4]. Moreover, in most of the currently available photogrammetric platforms, the UAVs use GPS pseudorange measurements in the Single Point Positioning (SPP) mode to determine the trajectory for photogrammetric studies [5]. The SPP positioning mode is the basic solution in air navigation for low-cost UAVs. The SPP method uses the L1-C/A code measurements in conjunction with data provided in the navigation message broadcast. Based on the onboard ephemeris, it is possible to calculate the satellites' coordinates, the corrections of the satellite clocks, the ionospheric correction, the relativistic correction, etc. [6]. The satellite position was estimated using the model of Keplerian orbit. In addition the clock bias is modelled using 2nd-degree polynomial model based on broadcast data. Finally, ionosphere correction is evaluated using the Klobuchar model for broadcast data [7]. It is the low-cost method of determining the position in real-time with an accuracy limited to a few meters. A significant advantage of this method is that it is used in every GNSS (Global

Related Works
Unmanned aerial vehicles have a wide range of applications both on the civil and military market. The demand for UAV technology is continuously growing, both in the context of performing air operations within visual line of sight (VLOS) and beyond visual line of sight (BVLOS). In Europe, the initiative of unified European airspace is referred to as U-Space. It includes the regulatory and technical framework for the gradual extension of the UAV's operability to shared airspace, also based on the Unmanned Aerial Systems (UTM) Traffic Management system. In this context, EGNOS and Galileo's use is the main chance to improve UAV navigation quality with regard to the integrity and accuracy of the positioning [18].
Currently, about 80% of GNSS receivers on the market, even the low-cost ones (including autopilot systems for UAVs), can receive corrections from SBAS systems [9,19]. The EGNOS corrections transmitted via the signal on the L1 band are beneficial to saving the battery power of the UAV and extending its flight time, as it enables the improvement of the accuracy of its position determination without additional communication channels and receiving corrections from the ground control system. It also increases the simplicity of system configuration.
Information from EGNOS is the basis for separation from other aircraft and will also play an essential role for UAVs in the future when they are integrated into civilian airspace. As mentioned before, during beyond visual line of sight (BVLOS) flights, it will be possible to detect and avoid other aircraft in a timely and appropriate manner. Moreover, thanks to information from the EGNOS system, it will be possible to increase this type of operation's security. Finally, EGNOS can also help to improve the safety functions of pilotage. This concept aims to use geographic data to establish flight boundary zones to prevent unsafe UAVS flights in sensitive areas. That can be used to restrict flights close to airfields or Energies 2021, 14, 739 3 of 22 above certain altitudes. A navigation solution for positioning based on EGNOS provides the means to determine whether the UAV is crossing a defined boundary zone by assessing the calculated position and the level of certainty [20]. Positioning information is expected to play a key role in achieving the ultimate goal, which is the full inclusion of UAVs in the airspace, and SBAS is by far the best, easiest and most realistic way to achieve this [9]. Currently, there are many UAVs on the market equipped with GNSS modules that support SBAS [21,22].
Much research has also been done into incorporating EGNOS information into lowaltitude imaging. In studies on the acquisition of multispectral imagery from low altitudes [23], taking into account corrections from the EGNOS system made it possible to achieve a sub-meter accuracy of registration of the approximate exterior orientation elements (EOE). However, in study [24], a GNSS sensor assisted by EGNOS was used to increase the accuracy of recording thermal imagery's position from a low altitude.
An interesting application of the SBAS system in UAV technology is the support of the SAR (Search And Rescue) system. The works [25][26][27] show the SBAS system's use as part of the SoL (Safety of Life) service. The SoL service in the EGNOS system is designed to offer an air navigation safety service by aircraft. In the context of the use of UAV technology, this service is to improve the safety of air operations.
The work [28] shows the results of research on using the EGNOS SISNET system for precise positioning of UAVs in Poland. The tests were carried out in real-time and in the post-processing mode. The work also uses a network of GNSS reference stations located in Poland to determine the actual UAV flight trajectory.
The papers [29,30] describe the SBAS system's use in UAV positioning during the landing approach procedure. Analyzing in detail, the work [29] shows the use of SBAS corrections for UAV positioning in the absence of data from the GBAS support system. However, in [30] a model of the GNSS / SBAS approach without using ILS (Instrument landing system) was presented.
An interesting navigation solution was shown in [31], where the UAV flight trajectory was determined, and position errors in the horizontal plane of the flight were determined. The GPS / SBAS solution was included in the calculations. It should be emphasized that the research was conducted based on the PBN (Performance-Based Navigation) concept for UAV positioning.
This study aims to propose a new UAV position solution using EGNOS corrections. The solution used is based on a weighted mean model that considers a single GPS/EGNOS position solution from satellites S123, S126 and S136. The numerical solution's construction uses the weighting of the measurements as a function of the inverse of the square of the mean error of the UAV coordinates in the BLh ellipsoidal frame. Moreover, the proposed numerical solution enables the analysis of the accuracy of the determined resultant UAV coordinates. Such a position model is characterized by a uniform influence of the particular GPS/EGNOS solution on the resultant UAV coordinates. Therefore, it can be concluded that the use of the GPS/EGNOS solution from satellites S123, S126 and S136 causes the resultant UAV coordinates to be defined for the number of degrees of freedom equal to 2.
Finally, the research contributions are focused on: first stage, the methodology describing the mathematical model's scheme using EGNOS corrections from one SBAS satellite will be presented. Next, the second diagram will present the methodology for determining the UAV position based on the GPS solution with the use of corrections from several SBAS satellites. In addition, for the second stage of calculations, analysis of the EGNOS positioning accuracy for UAVs will also be presented.

Designation the UAV Position Based on a Single GPS+EGNOS Solution
The methodology of the first stage is based on the implementation of EGNOS corrections from one SBAS satellite to the mathematical model of GPS positioning in aviation, as presented below [32,33]: where: l-code measurements in the GPS receiver, d GPS/EGNOS -geometric distance satellite-receiver, the long-term and fast EGNOS corrections are applied for designation the geometric distance satellite- (X, Y, Z)-geocentric coordinates of UAV vehicle; X GPS/EGNOS = Xs GPS + δx EGNOS , GPS satellite coordinate along the X axis, Y GPS/EGNOS = Ys GPS + δy EGNOS , GPS satellite coordinate along the Y axis, Z GPS/EGNOS = Zs GPS + δz EGNOS , GPS satellite coordinate along the Z axis, (Xs GPS , Ys GPS , Zs GPS )-GPScoordinatesbasedonephemerisdata, (δx EGNOS , δy EGNOS , δ EGNOS )the long-term and fast EGNOS corrections, c-the speed of light, dtr-receiver clock bias, dts GPS/EGNOS -satellite clock bias, the long-term and fast EGNOS corrections are applied for designation the satellite clock bias, dts GPS/EGNOS = dts GPS + δdts EGNOS , dts GPSsatellite clock bias based on ephemeris data, δdts EGNOS -the long-term and fast EGNOS corrections, Ion EGNOS -ionosphere correction, based on SBAS ionosphere model, Ion EGNOS = MF· ∑ w n ·VTEC n , MF-mapping function, w n -weight, the function of the distance between GRID coordinates in each node and current ionosphere pierce point coordinates, VTEC n -Vertical TEC at each node of GRID coordinates, Trop EGNOS -troposphere correction, based on RTCA-MOPS model, , (ZHD 0 , ZWD 0 )-ZHD and ZWD term at sea level, (λ, β)-water vapour lapse rate and temperature lapse rate, (g, R d )constant coefficients, h-ellipsoidal height, T K -temperature, El-elevation angle, Relrelativistic effect in the GPS receiver, TGD-Timing Group Delay for GPS satellite, Mpmultipath effect in the GPS receiver.
Based on Equation (1), the UAV coordinates are determined in the GNSS navigation system based on a single GPS+EGNOS solution using the least-squares method [34,35], such as below: where: Q X -vector with unknown parameters, N = A T ·p·A-matrix of normal equation frame, A-matrix of plan, v-residuals vector, p-matrix of weights, p = 1 m0 2 ·ml 2 , m0-unit mean error a priori, m0 = 1, ml-pseudorange measurement error matrix, ml = ml 0 sin(El) 2 + m 2 SBAS/EGNOS -mean error of pseudorange, ml 0 -mean error of pseudorange in GPS system, ml 0 = 1 m, m SBAS/EGNOS -mean error of SBAS correction model, L = A T ·p·dl-vector of absolute terms, dl-vector with difference between measurements and modeled parameters, m0 post -unit mean error a posteriori, n-number of observations for each measurement epoch, k-number of unknown parameters for each measurement The position of the UAV is determined with a specific time interval, usually every 1 s. In the observation Equation (1), EGNOS corrections for ephemeris data and GPS satellite clocks, ionospheric and tropospheric delay are applied. The obtained position based on Equation (2) is referenced to XYZ coordinates. In addition, the UAV position is better determined in aerial navigation in ellipsoidal coordinates BLh (B-Latitude, L-Longitude, h-ellipsoidal height). The position of the UAV in BLh ellipsoidal coordinates is described as follows [35]: where: (a, b)-ellipsoid semi-major axis and semi-minor axis, e-first flattening, e = a 2 −b 2 a 2 , R-radius of the curvature in the prime vertical, The mean error values m Qx of the calculated geocentric XYZ coordinates for the UAV may also be expressed in the BLh ellipsoidal coordinate frame, as shown below [36,37]: where: m BLh -variance-covariance matrix of parameters in the BLh ellipsoidal frame, m BLh = R·C Qx ·R T , R-conversion matrix from the XYZ geocentric frame to the BLh geodetic frame, mB-mean error for geodetic latitude B, mL-mean error for geodetic longitude L, mh-mean error for ellipsoidal height h. The calculated position of the UAV in the form of BLh coordinates and their mean errors (mB, mL, mh) will be used for the next stage of the research. It should be noted that Equations (1)-(4) denote an algorithm for determining the UAV position for GPS positioning using SBAS corrections from a single EGNOS satellite. Therefore, it can be said that the Equations (1)-(4) is a single GPS + EGNOS solution in UAV positioning.
Mathematical Equation (5) provides a weighted mean model to determine the resultant UAV position based on the EGNOS corrections from several SBAS satellites. The mathematical model in Equation (5) is based on the application of EGNOS corrections from SBAS satellites S123, S126 and S136.
The measurement weights were defined as the inverse of the squared mean error values of the determined coordinates as below: ph GPS/S126 = 1 mh 2 GPS/S126 ph GPS/S136 = 1 mh 2 GPS/S136 (6) where: (mB GPS/S123 , mL GPS/S123 , mh GPS/S123 )-mean error values of the UAV coordinates determined using the EGNOS corrections from satellite S123, (mB GPS/S126 , mL GPS/S126 , mh GPS/S126 ) -mean error values of the UAV coordinates determined using the EGNOS corrections from satellite S126, (mB GPS/S136 , mL GPS/S136 , mh GPS/S136 )-mean error values of the UAV coordinates determined using the EGNOS corrections from satellite S136.

Accuracy Analysis of the GPS + EGNOS Solution
The accuracy analysis of the presented research method is based on a three-step calculation scheme, as shown below: (I) establishing UAV position errors determined as absolute values from the comparison between the resultant GPS+EGNOS solution and the reference position calculated by the RTK differential technique [38]: (II) establishing statistical accuracy in the form of the RMS error [39]: (III) establishing the mean absolute error: In the first stage of the accuracy analysis, the UAV position errors are calculated as the difference between the resultant GPS + EGNOS solution and the RTK differential method. RMS errors are determined sequentially, which are a statistical measure of accuracy along each BLh axis. Additionally, the accuracy analysis also includes the calculation of the absolute error along each BLh axis.
The mathematical scheme (12)- (14) presented the comparison between multiple GPS+EGNOS solutions and the RTK reference trajectory. This approach is used to determine the UAV position's accuracy in reference to the RTK measurement technique. During the flight, the UAV position was estimated using the SPP GPS solution as a basic positioning model of the flight trajectory. It is necessary to calculate the difference between the multiple GPS+EGNOS solutions and the SPP GPS real model of the UAV's trajectory. It should be evaluated as below [40]: where: (rB, rL, rh)-difference between the multiple GPS+EGNOS solution and the real UAV trajectory from the SPP GPS model, (B SPP,real , L SPP,real , h SPP,real )-UAV trajectory from the SPP GPS model. Equation (15) shows the convergence between the multiple GPS+EGNOS solution and the UAV trajectory from the SPP GPS model. It is a good indicator in navigational analyses of the obtained UAV position.

Research Test
As part of the research analysis, the GNSS data recorded by an AsteRx-m2 UAS receiver placed on the Tailsitter unmanned platform was used. The flying platform was equipped with a single-frequency GNSS receiver, recording data at a frequency of 10 Hz. The Tailsitter platform is a type of VTOL aircraft. The general characteristics of the UAV platform include: capacity: 800 g payload, wingspan: 1.25 m, empty weight: 3.6 kg, max takeoff weight: 4.4 kg. The Tailsitter platform model is presented in Figure 1. In the first stage of the accuracy analysis, the UAV position errors are calculated as the difference between the resultant GPS + EGNOS solution and the RTK differential method. RMS errors are determined sequentially, which are a statistical measure of accuracy along each BLh axis. Additionally, the accuracy analysis also includes the calculation of the absolute error along each BLh axis.
The mathematical scheme (12)(13)(14) presented the comparison between multiple GPS+EGNOS solutions and the RTK reference trajectory. This approach is used to determine the UAV position's accuracy in reference to the RTK measurement technique. During the flight, the UAV position was estimated using the SPP GPS solution as a basic positioning model of the flight trajectory. It is necessary to calculate the difference between the multiple GPS+EGNOS solutions and the SPP GPS real model of the UAV's trajectory. It should be evaluated as below [40]: where: ( , , ℎ)-difference between the multiple GPS+EGNOS solution and the real UAV trajectory from the SPP GPS model, ( , , , , ℎ , )-UAV trajectory from the SPP GPS model. Equation (15) shows the convergence between the multiple GPS+EGNOS solution and the UAV trajectory from the SPP GPS model. It is a good indicator in navigational analyses of the obtained UAV position.

Research Test
As part of the research analysis, the GNSS data recorded by an AsteRx-m2 UAS receiver placed on the Tailsitter unmanned platform was used. The flying platform was equipped with a single-frequency GNSS receiver, recording data at a frequency of 10 Hz. The Tailsitter platform is a type of VTOL aircraft. The general characteristics of the UAV platform include: capacity: 800 g payload, wingspan: 1.25 m, empty weight: 3.6 kg, max takeoff weight: 4.4 kg. The Tailsitter platform model is presented in Figure 1. The UAV platform performed a photogrammetric flight in south-eastern Poland. GNSS navigation data was acquired on April 23rd 2020, between 11:00 a.m.-2:00 p.m. local time. The temperature was 17 °C, while the wind speed was around 2 m/s. Figure 2 shows the UAV flight trajectory in the ellipsoidal coordinate frame BL (B-geodetic Latitude, L-geodetic Longitude). The flight trajectory in Figure 2 is the resultant of the UAV position for horizontal coordinates, calculated according to Equation (7). The minimum value of B component equals 50.499204°, whereas the maximum values were about The UAV platform performed a photogrammetric flight in south-eastern Poland. GNSS navigation data was acquired on 23 April 2020, between 11:00 a.m.-2:00 p.m. local time. The temperature was 17 • C, while the wind speed was around 2 m/s. Figure 2 shows the UAV flight trajectory in the ellipsoidal coordinate frame BL (B-geodetic Latitude, Lgeodetic Longitude). The flight trajectory in Figure 2 is the resultant of the UAV position for horizontal coordinates, calculated according to Equation (7)   The maximum flight altitude was approximately 250 m relative to the height of the route starting point. Figure 3 shows the UAV flight's vertical trajectory as a function of time when the position was recorded. The flight trajectory in Figure 3 is the resultant UAV position for the vertical component, calculated by the Equation (7). The h component changed during the flight from 296.87 to 545.12 m.
The experiment was carried out from 13:02:56 to 13:39:35 according to GPS Time. It was, therefore, necessary to determine the state of the ionosphere during the UAV flight. The VTEC parameter was used to determine the ionosphere's state, which was then determined based on the SBAS ionosphere model. The VTEC parameter values are determined for EGNOS corrections from tracked and available SBAS satellites PRN S123, S126 and S136. The free SBAS Mentor v.1.15 software was used to perform the calculations [41]. Table 1 shows the results of the VTEC parameter. Based on the EGNOS corrections, it was found that the VTEC values for satellites S123, S126 and S136 are the same and   The maximum flight altitude was approximately 250 m relative to the height of the route starting point. Figure 3 shows the UAV flight's vertical trajectory as a function of time when the position was recorded. The flight trajectory in Figure 3 is the resultant UAV position for the vertical component, calculated by the Equation (7). The h component changed during the flight from 296.87 to 545.12 m.
The experiment was carried out from 13:02:56 to 13:39:35 according to GPS Time. It was, therefore, necessary to determine the state of the ionosphere during the UAV flight. The VTEC parameter was used to determine the ionosphere's state, which was then determined based on the SBAS ionosphere model. The VTEC parameter values are determined for EGNOS corrections from tracked and available SBAS satellites PRN S123, S126 and S136. The free SBAS Mentor v.1.15 software was used to perform the calculations [41]. Table 1 shows the results of the VTEC parameter. Based on the EGNOS corrections, it was found that the VTEC values for satellites S123, S126 and S136 are the same and The experiment was carried out from 13:02:56 to 13:39:35 according to GPS Time. It was, therefore, necessary to determine the state of the ionosphere during the UAV flight. The VTEC parameter was used to determine the ionosphere's state, which was then determined based on the SBAS ionosphere model. The VTEC parameter values are determined for EGNOS corrections from tracked and available SBAS satellites PRN S123, S126 and S136. The free SBAS Mentor v.1.15 software was used to perform the calculations [41]. Table 1 shows the results of the VTEC parameter. Based on the EGNOS corrections, it was found that the VTEC values for satellites S123, S126 and S136 are the same and amounted to 1.375 m during the research test. Moreover, Table 1 presents the GIVE parameter index values, which was constant during the research test and was equal to 8. In the first phase of the test, the UAV position calculations were performed in the RTKPOST library of the RTKLIB program (see Equation (1)). Therefore, in the RTKPOST application, the configuration was set as follows [42,43]:  [44]. The mathematical analysis of calculations in Scilab applies to Equations (2)- (7). The test and simulation results are presented in Section 5.

Results
The EGNOS positioning concept assumes that the corrections from SBAS satellites should be identical. However, in many research studies, the problem of the quality of EGNOS corrections was noticed and analyzed. For example, in [45], the accuracy of GPS + EGNOS positioning at airfields in Poland was shown. The authors of the work pointed out that the corrections transmitted by the EGNOS satellites are different, which directly impacts the accuracy of the determined ellipsoidal BLh coordinates for the reference stations.
Another work [46] shows the differences in the coordinates of three GNSS reference stations installed on the aircraft flight path. Based on the performed tests, it was found that there are quite significant differences in the XYZ coordinates of the reference stations based on the EGNOS corrections applied from satellites S120, S124 and S126 and the catalogue position of the base stations. In addition, the paper [47] presents the results of research on the application of EGNOS corrections from satellites S120, S124 and S126 in air navigation. The paper describes the differences between individual GPS + EGNOS solutions determining the position of the aircraft. Based on the research results obtained, it can be seen that the quality of EGNOS corrections significantly influenced the position of the aircraft during the test flight. Moreover, the aircraft's determined position from a particular GPS + EGNOS solution was different, which only emphasizes that the problem of the impact of EGNOS corrections on kinematic positioning in air navigation should be addressed. Therefore, this work and the proposed solution of the weighted average model may significantly influence the further development of scientific works in the context of the implementation of SBAS in aviation.
In support of the considerations on EGNOS corrections' quality in the first stage of the research, the results of the difference of the estimated BLh ellipsoidal coordinates of the UAV position between the individual GPS + EGNOS solutions were shown, in accordance with the Equation (1). Figure 4 shows the results of the B coordinate difference for the UAV position between the various GPS + EGNOS solutions. The B coordinate difference between GPS + S123 and GPS + S126 is from −0.42 to +1.04 m, and the arithmetic mean of the difference is 0.04 m. The B coordinate difference between GPS + S123 and GPS + S136 is from −0.09 to +0.14 m, and the arithmetic mean of this difference is 0.01 m. It can be stated that solutions GPS + S123 and GPS + S136 overlap and are quite closely aligned with each other. On the other hand, the difference of the B coordinate between the solutions GPS + S126 and GPS + S136 ranges from −1.04 to +0.41 m, and the arithmetic mean of this difference is −0.03 m.
In support of the considerations on EGNOS corrections' quality in the first stage of the research, the results of the difference of the estimated BLh ellipsoidal coordinates of the UAV position between the individual GPS + EGNOS solutions were shown, in accordance with the Equation (1). Figure 4 shows the results of the B coordinate difference for the UAV position between the various GPS + EGNOS solutions. The B coordinate difference between GPS + S123 and GPS + S126 is from −0.42 to +1.04 m, and the arithmetic mean of the difference is 0.04 m. The B coordinate difference between GPS + S123 and GPS + S136 is from −0.09 to +0.14 m, and the arithmetic mean of this difference is 0.01 m. It can be stated that solutions GPS + S123 and GPS + S136 overlap and are quite closely aligned with each other. On the other hand, the difference of the B coordinate between the solutions GPS + S126 and GPS + S136 ranges from −1.04 to +0.41 m, and the arithmetic mean of this difference is −0.03 m.  Figure 5 shows the results of the difference in the L coordinate for the UAV position between the various GPS + EGNOS solutions. The L-coordinate difference between GPS + S123 and GPS + S126 is from −1.83 to +0.56 m, and the arithmetic mean of this difference is −0.27 m. The L-coordinate difference between GPS + S123 and GPS + S136 is from −0.16 to +0.14 m, and the arithmetic mean of this difference is −0.01 m. In this case, it is worth noting that the solutions GPS + S123 and GPS + S136 coincide and are quite closely aligned. In turn, the difference in the L coordinate between the solutions GPS + S126 and GPS + S136 ranges from −0.65 to +1.80 m, and the arithmetic mean of this difference is +0.26 m.  Figure 5 shows the results of the difference in the L coordinate for the UAV position between the various GPS + EGNOS solutions. The L-coordinate difference between GPS + S123 and GPS + S126 is from −1.83 to +0.56 m, and the arithmetic mean of this difference is −0.27 m. The L-coordinate difference between GPS + S123 and GPS + S136 is from −0.16 to +0.14 m, and the arithmetic mean of this difference is −0.01 m. In this case, it is worth noting that the solutions GPS + S123 and GPS + S136 coincide and are quite closely aligned. In turn, the difference in the L coordinate between the solutions GPS + S126 and GPS + S136 ranges from −0.65 to +1.80 m, and the arithmetic mean of this difference is +0.26 m.   Figure 6 shows the results of the h coordinate difference for the UAV position between the various GPS + EGNOS solutions. The L-coordinate difference between GPS + S123 and GPS + S126 is from −2.92 to +0.60 m, and the arithmetic mean of this difference is -0.41 m. The h-coordinate difference between GPS + S123 and GPS + S136 is from −0.30 to +0.18 m, and the arithmetic mean of this difference is −0.01 m. For the h component, the results of the solutions GPS + S123 and GPS + S136 coincide and are quite convergent with each other. In turn, the difference in the h coordinate between the solutions GPS + S126 and GPS + S136 ranges from −0.75 to +2.92 m, and the arithmetic mean of this difference is +0.40 m.   Figure 6 shows the results of the h coordinate difference for the UAV position between the various GPS + EGNOS solutions. The L-coordinate difference between GPS + S123 and GPS + S126 is from −2.92 to +0.60 m, and the arithmetic mean of this difference is -0.41 m. The h-coordinate difference between GPS + S123 and GPS + S136 is from −0.30 to +0.18 m, and the arithmetic mean of this difference is −0.01 m. For the h component, the results of the solutions GPS + S123 and GPS + S136 coincide and are quite convergent with each other. In turn, the difference in the h coordinate between the solutions GPS + S126 and GPS + S136 ranges from −0.75 to +2.92 m, and the arithmetic mean of this difference is +0.40 m.
From the calculated BLh coordinate differences for the UAV positions, it can be seen that the EGNOS corrections transmitted by the individual satellites S123, S126, S136 are different. In particular, the difference in coordinates is significant when applying the corrections from satellite S126. From the calculated BLh coordinate differences for the UAV positions, it can be seen that the EGNOS corrections transmitted by the individual satellites S123, S126, S136 are different. In particular, the difference in coordinates is significant when applying the corrections from satellite S126.
Therefore, the impact of individual corrections from EGNOS satellites on the determination of the resultant position of the UAV according to the proposed model of Equations (1)-(8) will be visible. The developed model of solving Equations (1)- (8) to determine the resultant position of the UAV is correct, assuming obtaining different values of BLh coordinates from a single GPS + EGNOS solution (see Equation (1)).
As part of the next step, Figure 7 shows the results of the standard deviation parameters (δB, δL, δh) in the BLh ellipsoidal coordinates. The δB parameter values ranged from 0.01 to 0.62 m. Moreover, the arithmetic mean for the parameter δB is equal to 0.11 m. The δL parameter values ranged from 0.01 to 1.14 m. The arithmetic mean for the δL parameter is 0.20 m. The δh parameter values ranged from 0.01 to 0.81 m. Moreover, the arithmetic mean for the δh parameter is 0.13 m. It is worth noting that from the measuring epoch 47,500 s to 48,500 s, the (δB, δL, δh) parameter values continue to increase up to the maximum value of 1.14 m. Ignoring this period, the (δB, δL, δh) parameter values do not exceed 0.25 m. It should be noted that the highest values of the standard deviation are visible for the L component. The values of the standard deviation along the L axis were determined according to Equation (8). Therefore, the distribution of corrections has the most significant impact on the values of the standard deviation along the L axis v L,i , v L,j , v L,k and the measurement weights (pL GPS/S123 , pL GPS/S126 , pL GPS/S136 ). In the stochastic model, the value of the standard deviation is significantly influenced by the measurement weights compared to the model of equally accurate measurements shown in Equation (18).
( / 123 , / 126 , / 136 ). In the stochastic model, the value of the standard deviation is significantly influenced by the measurement weights compared to the model of equally accurate measurements shown in Equation (18).
The standard deviation values along the L axis are greater than for the horizontal component B or the vertical component h.  Table 2 shows the results of the arithmetic mean for the obtained values of the proportionality coefficients ( , ℎ , ℎ ). Based on the obtained results, it can be concluded that the best fit of the standard deviation values is between the h and L components. The values of the proportionality coefficients vary from 2.09 to 2.89. The proportionality coefficient is important because it determines the relations between the values of standard deviations of individual BLh components. The closer its value is to 1, the more the standard deviation's distribution is in a similar trend. The parameters ( , ℎ , ℎ ) were estimated as below:  Table 2 shows the results of the arithmetic mean for the obtained values of the proportionality coefficients (rBL, rhB, rhL). Based on the obtained results, it can be concluded that the best fit of the standard deviation values is between the h and L components. The values of the proportionality coefficients vary from 2.09 to 2.89. The proportionality coefficient is important because it determines the relations between the values of standard deviations of individual BLh components. The closer its value is to 1, the more the standard deviation's distribution is in a similar trend. The parameters (rBL, rhB, rhL) were estimated as below: (16) where: δB, δL, δh -mean value of standard deviations along to BLh axis. The next stage of the research concerns the analysis of the UAV positioning accuracy. Figure 8 presents the position errors of UAV. In the first part of the analysis, the (DB, DL, Dh) parameters were determined in accordance with Equation (12). The position error values of the (DB, DL, Dh) parameters do not exceed the values of 2.01, 2.09 and 3.98 m, respectively. Table 3 shows the results of the obtained accuracy parameters (RMS DB , RMS DL , RMS Dh ), calculated according to Equation (13). Based on the obtained results, it can be concluded that the highest accuracy is noticeable for the B component The next stage of the research concerns the analysis of the UAV positioning accuracy. Figure 8 presents the position errors of UAV. In the first part of the analysis, the ( , , ℎ) parameters were determined in accordance with Equation (12). The position error values of the ( , , ℎ) parameters do not exceed the values of 2.01, 2.09 and 3.98 m, respectively.  Table 3 shows the results of the obtained accuracy parameters ( , , ℎ ), calculated according to Equation (13). Based on the obtained results, it can be concluded that the highest accuracy is noticeable for the B component     Table 4 shows the results of the obtained mean absolute error values DB, DL, Dh , calculated according to the Equation (14).  In the next step, the BLh coordinates from the GPS + EGNOS solution were compared to the actual flight trajectory from the SPP GPS model based on Equation (15). The differences in coordinates (rB, rL, rh) are shown in Figure 9. The coordinate difference ranges from −0.99 to +0.22 m for the B component, −0.46 to +1.21 m for the L component, for the h component +3,90 to 6.64 m. Based on the obtained (rB, rL, rh) results, a significant shift along the vertical axis h can be seen. The mean value of the differences in coordinates are: −0.34 m for the B coordinate, +0.47 m for the L coordinate, and 5.30 m for the h coordinate. It is worth paying attention to the numerous changes in the coordinate difference in Figure 9. The reason for this phenomenon can be found in the number of tracked GPS satellites. Thus, for the measurement epoch 47,500 s, the number of tracked GPS satellites for the SBAS123 and SBAS136 solution was 7. On the other hand, for the same moment in time from the SBAS126 solution, the number of GPS satellites changed from 8 to 9. Another jump is visible for the measurement epoch 48,500 s. Therefore, with SBAS123 and SBAS136, 8 GPS satellites are tracked, while SBAS126 changes the number of tracked GPS satellites from 10 to 9. Successive changes in the coordinate difference can be explained by the change in the number of GPS satellites being tracked, which obviously affects the UAV ellipsoidal coordinates' performance. from 8 to 9. Another jump is visible for the measurement epoch 48500 s. Therefore, with SBAS123 and SBAS136, 8 GPS satellites are tracked, while SBAS126 changes the number of tracked GPS satellites from 10 to 9. Successive changes in the coordinate difference can be explained by the change in the number of GPS satellites being tracked, which obviously affects the UAV ellipsoidal coordinates' performance.

Discussion
The discussion of the article will cover two main issues: analysis of the accuracy of the UAV positioning based on the proposed method against the SPP code solution without the use of EGNOS corrections, Figure 9. The difference of UAV coordinates between GPS+EGNOS solution and SPP.

Discussion
The discussion of the article will cover two main issues: analysis of the accuracy of the UAV positioning based on the proposed method against the SPP code solution without the use of EGNOS corrections, -analysis of the results of the standard deviation based on the proposed method against the arithmetic mean model.
The major highlights in discussion will concerned to: 1. improvement the accuracy of UAV position based EGNOS solution; 2.
improvement of the values of the standard deviation of BLh ellipsoidal coordinates based on the weighted mean model.
The first part of the discussion shows how the presented research method makes it possible to increase UAV positioning accuracy. For this purpose, two comparisons were presented. At first, the UAV coordinates based on GPS + EGNOS solution (Equations (1)-(5)) were compared in relation to the reference position determined by the RTK technique. At second, the GPS SPP solution (without EGNOS corrections) [48] were compared with RTK true trajectory. Figure 10 shows the position error results along the B axis. The position error values along the B axis were calculated according to Equation (12). It should be noted that the DB position error results for the EGNOS solution are presented and analyzed in Figure 8. Figure 10 additionally shows position errors determined without the use of EGNOS corrections. It should be noted that the DB position errors without the use of EGNOS corrections range from +0.11 to +2.32 m. Additionally, the arithmetic mean is +1.21 m. Comparing both graphs, it can be seen that the EGNOS corrections significantly improved the accuracy of the UAV positioning for the B component. This increase in accuracy is as much as 29%. Figure 11 shows the position error results along the L axis. The position error values along the L axis were calculated according to Equation (12). It should be noted that the DL position error results for the EGNOS solution are presented and analyzed in Figure 8. Figure 11 additionally shows position errors determined without the use of EGNOS corrections. It should be noted that the DL position errors without the use of EGNOS corrections range from −0.49 to +2.46 m. Additionally, the arithmetic mean is 1.02 m. Comparing both graphs, it can be seen that the EGNOS corrections significantly improved the accuracy of the UAV positioning for the L component. This increase in accuracy is around 46%.
corrections. It should be noted that the position errors without the use of EGNOS corrections range from +0.11 to +2.32 m. Additionally, the arithmetic mean is +1.21 m. Comparing both graphs, it can be seen that the EGNOS corrections significantly improved the accuracy of the UAV positioning for the B component. This increase in accuracy is as much as 29%.   Figure 11 shows the position error results along the L axis. The position error values along the L axis were calculated according to Equation (12). It should be noted that the position error results for the EGNOS solution are presented and analyzed in Figure 8. Figure 11 additionally shows position errors determined without the use of EGNOS corrections. It should be noted that the position errors without the use of EGNOS corrections range from −0.49 to +2.46 m. Additionally, the arithmetic mean is 1.02 m. Comparing both graphs, it can be seen that the EGNOS corrections significantly improved the accuracy of the UAV positioning for the L component. This increase in accuracy is around 46%.  Figure 12 shows the position error results along the h axis. The position error values along the Z axis were calculated according to Equation (12). It should be noted that the Dh position error results for the EGNOS solution are presented and analyzed in Figure 8. corrections. It should be noted that the position errors without the use of EGNOS corrections range from −0.49 to +2.46 m. Additionally, the arithmetic mean is 1.02 m. Comparing both graphs, it can be seen that the EGNOS corrections significantly improved the accuracy of the UAV positioning for the L component. This increase in accuracy is around 46%.  Figure 12 shows the position error results along the h axis. The position error values along the Z axis were calculated according to Equation (12). It should be noted that the  It should be noted that the ellipsoidal height h and its accuracy Dh are dependent on the geometrical factor, i.e., the satellite-receiver position system. Mainly, the long-term corrections for the satellite position have a decisive influence on the performance of the vertical coordinate. It should not be forgotten that in GNSS measurements, the tropospheric correction, as one of the atmospheric corrections, also affects the ellipsoidal height value's determination. Summing up, the accuracy of ellipsoidal height determination is influenced by the long-term corrections for ephemeris data and the tropospheric correction for RTCA-MOPS model.
In the second part of the discussion, the authors presented the research analysis results concerning the comparison of the standard deviation value for the UAV position determination model based on the weighted mean and arithmetic mean methods. The mathematical scheme of determining the resultant UAV position using the arithmetic mean model can be written as below: where: i + j + k = 3. In turn, the values of the standard deviation from the weighted mean model can be determined as shown below:  Figure 13 shows the results of standard deviations δB determined from Equations (8) and (11). The δB parameter results for the weighted mean model are shown and described  Figure 7. The maximum results of the δB parameter from the arithmetic mean model are equal to 0.64 m. The mean value of the δB parameter is 0.14 m. Comparing the plots it can be seen that the standard deviation values for the B component from the weighted mean model are smaller than for the arithmetic mean model. The improvement of the standard deviation results from the weighted mean model is about 21% compared to the arithmetic mean method.       (8) and (18). The δL parameter results for the weighted mean model are shown and described in Figure 7. The maximum results of the δL parameter from the arithmetic mean model are equal to 1.16 m. The mean value of the δL parameter is 0.27 m. Comparing the plots it can be seen that the standard deviation values for the L component from the weighted mean model are smaller than for the arithmetic mean model. The standard deviation results from the weighted mean model are about 26% higher than the arithmetic mean method. Figure 15 shows the results of standard deviations δh determined from Equations (8) and (11). The δh parameter results for the weighted mean model are shown and described in Figure 7. The maximum results of the δh parameter from the arithmetic mean model are equal to 1.69 m. The mean value of the δh parameter is 0.26 m. Comparing the plots, it can be seen that the standard deviation values for the Z component from the weighted mean model are smaller than for the arithmetic mean model. The standard deviation results from the weighted mean model are about 50% higher than to the arithmetic mean method.
The research results presented in the discussion only show how important an element in the positioning and orientation in the UAV space is selecting an appropriate mathematical model for the description of the flight trajectory. Similar research work has been analyzed in other scientific articles. Thus, the work [50] shows a mathematical model for determining the UAV position in relation to the base station and the orientation of the UAV using the HPR rotation angles (Heading, Pitch, Roll). Additionally, the paper shows the error in determining the UAV position as a function of the distance from the base station and the errors of HPR rotation angles. The mathematical model shown in [50] has unique properties because it determines the positions and orientations of the UAV in relation to the base station. It can be said that the presented solution is a differential model of the position and orientation of the UAV.
This work shows how important the GNSS satellite technique is today. That is why more and more research works in the world concern the implementation of the GNSS sensor for precise positioning of the UAV. Therefore, in [51] a mathematical model for determining the position and orientation of an aircraft with the use of the single-frequency method for several onboard receivers is shown. Moreover, the paper defines the phase uncertainty solution using the LAMBDA and C-LAMBDA methods for single-frequency receivers. Additionally, in work [51], the effectiveness of determining the orientation of the aircraft for relative GNSS positioning with the use of phase observations and ground reference stations was discussed. The accuracy of the aircraft orientation was set at 0.2 • .    (8) and (18). The parameter results for the weighted mean model are shown and described in Figure 7. The maximum results of the parameter from the arithmetic mean model are equal to 1.16 m. The mean value of the parameter is 0.27 m. Comparing the plots it can be seen that the standard deviation values for the L component from the weighted mean model are smaller than for the arithmetic mean model. The standard deviation results from the weighted mean model are about 26% higher than the arithmetic mean method.  Figure 15 shows the results of standard deviations ℎ determined from Equations (8) and (11). The ℎ parameter results for the weighted mean model are shown and described in Figure 7. The maximum results of the ℎ parameter from the arithmetic mean model are equal to 1.69 m. The mean value of the ℎ parameter is 0.26 m. Comparing the plots, it can be seen that the standard deviation values for the Z component from the weighted mean model are smaller than for the arithmetic mean model. The standard deviation results from the weighted mean model are about 50% higher than to the arithmetic mean method. The research results presented in the discussion only show how important an element in the positioning and orientation in the UAV space is selecting an appropriate mathematical model for the description of the flight trajectory. Similar research work has been analyzed in other scientific articles. Thus, the work [50] shows a mathematical model for determining the UAV position in relation to the base station and the orientation of the UAV using the HPR rotation angles (Heading, Pitch, Roll). Additionally, the paper shows the error in determining the UAV position as a function of the distance from the base station and the errors of HPR rotation angles. The mathematical model shown in [50] has unique properties because it determines the positions and orientations of the UAV in relation to the base station. It can be said that the presented solution is a differential model of the position and orientation of the UAV.
This work shows how important the GNSS satellite technique is today. That is why more and more research works in the world concern the implementation of the GNSS sensor for precise positioning of the UAV. Therefore, in [51] a mathematical model for determining the position and orientation of an aircraft with the use of the single-frequency method for several onboard receivers is shown. Moreover, the paper defines the phase uncertainty solution using the LAMBDA and C-LAMBDA methods for A very interesting approach in UAV positioning may be the PPP-RTK measurement technique [52]. The PPP-RTK method includes elements of PPP positioning, i.e., precise orbits and clocks as well as the characteristics of the phase center of the satellite antenna (or receiver) and the parameters of the atmosphere state from the RTK solution, i.e., the ZTD tropospheric delay model and the VTEC ionospheric delay determination model. Of course, in the case of the PPP-RTK solution, it is necessary to determine the Float phase's uncertainty. For UAV positioning, the PPP-RTK method will increase the accuracy of the determined coordinates to the level of cm~dm. It will also reduce orbit errors, tropospheric and ionospheric delays, and satellite clock errors in the PPP-RTK observation model. In UAV positioning, selecting an appropriate ionosphere and troposphere model is particularly important because the ionospheric correction shortens the vector between the reference station and the UAV position. In contrast, the tropospheric correction affects the determination of the ellipsoidal height of the UAV. Thus, both atmospheric corrections are crucial here to improve the UAV's precise positioning performance.

Conclusions
The work concerns the development of a new resultant UAV position model based on the application of EGNOS corrections in GPS measurements. The paper presents the concept of determining the resultant UAV position on the basis of individual GPS/EGNOS solutions for the SPP code method. For this purpose, a weighted mean model was used to determine the ellipsoidal BLh coordinates of the UAV. The weighting of measurements considers the inverse of the square of the mean position error along the component axes of the BLh ellipsoidal frame. Accuracy measures in the form of position errors, the RMS error and the mean absolute error were determined for the computational strategy used. The new research method was tested to determine the position of the Tailsitter UAV.
On the basis of the performed tests, it was found that the mean values of the standard deviation of the BLh coordinates in the weighted mean model are better than 0.2 m. Moreover, the weighted mean model's use made it possible to increase the accuracy of the UAV position by about 29% for Latitude, 46% for Longitude and 72% for ellipsoidal height compared to the standard SPP positioning in the GPS receiver. It should be mentioned that the values of the standard deviation of the UAV position calculated from the weighted mean model improved by about 21 ÷ 50% compared to the solution from the arithmetic mean model. Such a significant improvement in accuracy in the weighted mean model only emphasizes the proposed solution's effectiveness. This research reinforces the thesis that the use of SBAS in UAV technology is necessary to improve the performance of GPS positioning on the L1 frequency. The authors plan to continue their work on the use of other augmentation systems such as SDCM or GAGAN in UAV technology in Poland.