A New Strategy for Improving the Accuracy of Aircraft Positioning Using DGPS Technique in Aerial Navigation

: In this paper a new mathematical algorithm is proposed to improve the accuracy of DGPS (Differential GPS) positioning using several GNSS (Global Navigation Satellites System) reference stations. The new mathematical algorithm is based on a weighting scheme for the following three criteria: weighting in function of baseline (vector) length, weighting in function of vector length error and weighting in function of the number of tracked GPS (Global Positioning System) satellites for a single baseline. The algorithm of the test method takes into account the linear combination of the weighting coefﬁcients and relates the position errors determined for single baselines. The calculation uses a weighting scheme for three independent baselines denoted as ( 1 A ,2 A ,3 A ) . The proposed research method makes it possible to determine the resultant position errors for ellipsoidal BLh coordinates of the aircraft and signiﬁcantly improves the accuracy of DGPS positioning. The analysis and evaluation of the new research methodology was checked for data from two ﬂight experiments carried out in Mielec and D˛eblin. Based on the calculations performed, it was found that in the ﬂight experiment in Mielec, due to the application of the new research methodology, DGPS positioning accuracy improved from 55 to 94% for all the BLh components. In turn, in the ﬂight experiment in D˛eblin, the accuracy of DGPS positioning improved by 63–91%. The study shows that the highest DGPS positioning accuracy is seen when using weighting criterion II, the inverse of the square of the vector length error.


Introduction
The GNSS (Global Navigation Satellites System), as a satellite technology, is continuously developed and widely used in many areas of human life, including the entire aviation industry [1]. However, it is important to note that the benefits of the GNSS in aviation are enormous, ranging from improved flight safety to reduced fuel emissions for example for RNAV (Area Navigation) [2] or the economic factor itself through improved airport capacity and an increased number of flight operations [3]. This is especially important as the GNSS, as one of many on-board sensors, is now a fundamental navigation sensor on board every aircraft. This makes the implementation of this sensor in aircraft positioning and the selection of an appropriate positioning method or technique particularly important for aviation science. GNSS satellite technology allows the user to determine the position of the aircraft using absolute and differential methods [4]. Within the presented work, the Differential DGPS (Differential GPS) method [5] in the GPS (Global Positioning System) satellite system will be discussed. In the differential DGPS technique, the location of the base stations relative to the mobile on-board GPS receiver mounted in the aircraft plays a key role [6]. Moreover, the number of GNSS reference stations involved in DGPS positioning is also very important [7]. The calculation of the aircraft position in DGPS technique usually uses code observations at the L1 frequency of the GPS system, to ensure

Related Papers
The use of the differential DGPS technique in air navigation is crucial for at least several reasons. Many applications of the DGPS method in air navigation can be found in the research literature. In particular, the application of the DGPS technology was evident in the following research areas: the use of the DGPS method as the primary method for aircraft positioning [6,12,13], -the use of the DGPS method to determine the reference position of the aircraft trajectory [14,15], -the use of the DGPS method as part of the GBAS (Ground Based Augmentation System) and the PA (Precision Approach) [16][17][18][19][20][21][22][23][24], -the application of the DGPS method in mapping applications, cartography and photogrammetry [25][26][27][28], -the application of the DGPS method in monitoring the ionosphere [29], -the use of the DGPS method in integration with various positioning sensors to determine flight altitude [30,31], -the determination of the accuracy and precision of DGPS positioning in air navigation [8,[32][33][34], -the determination of the integrity of positioning using the DGPS technique [34,35], -the integration of DGPS and DGLONASS (Differential GLONASS) solutions [6,13,36].
As shown in the research literature, the use of the DGPS technique in air navigation has been thematically extensive. For navigation, the most relevant is the development and elaboration of new mathematical algorithms to improve the accuracy of aircraft positioning, also using the DGPS technique. Therefore, the development of a new position navigation solution can significantly contribute to the improvement of the performance of the determined coordinates of the aircraft derived from the DGPS solution. It may also result in the improved safety of flight and flight operations. Consequently, the main idea of the paper is to develop an optimal algorithm to improve the accuracy of DGPS positioning in aviation using GNSS reference stations. The publication shows a computational strategy for determining the accuracy of the resultant aircraft position using a DGPS solution with three GNSS reference stations. The concept of the algorithm is based on the use of a weighted average model for three weighting criteria. Namely, in the research, a weighting scheme was developed as a function of the baseline (vector) length, the vector length error and the number of tracked GPS satellites. Such a weighting strategy is important for several factors. First, it relates the geometry of the distribution of GNSS reference stations and the flight position of the aircraft by calculating the vector length. Secondly, the weighting scheme is also dependent on the number of GPS satellites taken for the DGPS solution for a single baseline. Finally, the weighting process uses the vector length error values, which are actually the determined average coordinate errors of the stochastic model from the DGPS solution for a single baseline. One can notice that the choice of weighting parameters is not random; therefore, indicating which weighting factor is optimal for improving the accuracy of DGPS positioning for several baselines will be the main task to be carried out in the submitted scientific paper. The obtained research results will only confirm which weighting variant is the most optimal and best for DGPS positioning using several GNSS reference stations. This will also allow further research into the development of modern and advanced algorithms optimising the mathematical model of DGPS positioning in air navigation.
In summary, our main research achievements in the paper relate to the following: development of weighted average model algorithms to improve the accuracy of DGPS positioning using several GNSS reference stations, -implementation of a weighting model for the following three criteria: weighting as a function of baseline (vector) length, weighting as a function of vector length error, weighting as a function of the number of tracked GPS satellites, -definition of algorithms to calculate the resultant position errors, -conducting a formal evaluation of the research results obtained in percentage terms.

Research Method
Section 3 demonstrates the basic algorithm of the DGPS positioning method for a single baseline and a new mathematical model to improve the accuracy of DGPS positioning, including measurement weighting.

Basic DGPS Solution for a Single Baseline
The basic observation equation for the differential DGPS technique in the GPS navigation system for air navigation can be presented as follows [37][38][39]: where ∆-the operator of a single difference for GPS code measurements, it allows determining the difference in code measurements for two satellites tracked by one receiver; → 1A-the vector in the space between the GNSS base station ( → 1 ) and the GPS rover receiver ( → A) mounted onboard the aircraft; P 1A -the value of single code difference (expressed in meters) on the vector → 1A at the L1 frequency in the GPS system; ρ 1A -the geometric distance of the vector → 1A for the single code difference (expressed in geocentric coordinates XYZ) in the GPS system; ρ 1A = ρ 1 − ρ A , ρ 1 -the geometric distance between the GNSS base station ( → 1 ) and GPS satellite; (X 1 , Y 1 , Z 1 )-the geocentric coordinates of the GNSS base station ( → 1 ); (X sat , Y sat , Z sat )the geocentric coordinates of GPS satellite; ρ A -the geometric distance between the GPS rover receiver ( → A) mounted onboard the aircraft and GPS satellite, dinates of aircraft vehicle, e.g., GPS rover receiver ( → A) mounted onboard the aircraft; I 1A -the value of ionosphere delay at L1 frequency for single code difference in the GPS system; T 1A -the value of troposphere delay for single code difference in the GPS system; M 1A -the multipath effect and noise measurement at an L1 frequency for the GPS code measurements.
Equation (1) describes the mathematical model for DGPS positioning in air navigation for a single baseline. For this case, a single difference for GPS code observations on a vector ( → 1A) between a GNSS reference station and an on-board GPS receiver is considered. The position of the aircraft in geocentric XYZ coordinates (X A , Y A , Z A ) is estimated using Kalman filtering, as written below [36,40]: (a) process of "prediction": where A 0 -the matrix of coefficients, x 0 -the estimated values of the designated parameters a priori from the previous step, P 0 -the estimated values of covariance a priori from the previous step, x p -a prediction of state value, P p -the predicted covariance values, Q 0 -the variance matrix of the noise of the measurement process.
(b) process of "correction": where K k -the Kalman gain matrix, P k -the covariance matrix of parameters determined a posteriori, H-the matrix of partial derivatives, R-the covariance matrix of measurements, z-the vector of measured values, I-the unit matrix, x k -the parameters determined a posteriori.
For the purposes of air navigation, the position of an aircraft is much better expressed in ellipsoidal BLh coordinates, as recorded below [41]: where R-the radius of the curvature in the prime vertical, R = a √ 1−e 2 ×sin 2 B , (a, b)-the ellipsoid semi-major axis (a) and semi-minor axis (b), e-the first flattening, e = a 2 −b 2 a 2 , Furthermore, for each aircraft position expressed by the BLh components, their mean error values are determined as shown in the scheme below [42]: where m BLh -the variance-covariance matrix of parameters in the BLh ellipsoidal frame, m BLh = R × P k × R T , R-the conversion matrix from the XYZ geocentric frame to the BLh geodetic frame, mB A -the mean error for geodetic latitude B, mL A -the mean error for geodetic longitude L, mh A -the mean error for ellipsoidal height h.

New Solution to Improve Positioning Accuracy of DGPS Technique for Multiple Baselines
The assumption in Equation (1), however, has a geometric limitation, because the position of the aircraft is determined only from a single baseline, without controllability, which results in a lack of degrees of freedom and limits the verification of the positioning accuracy determination. A much better solution is to apply DGPS using several baselines. In the presented case, DGPS positioning accuracy was analysed on the basis of a solution using 3 baselines and a weight model. In the DGPS technique, 3 components of the aircraft position are determined; therefore, the minimum geometry of the system consisting of the GNSS reference station network and GPS receiver on-board should also have a minimum of 3. Additionally, such a configuration allows 2 degrees of freedom from 2 independent baselines, which means that the user always has control over the calculation of the determined coordinates of the aircraft position. This approach is very good, because it will always enable the elimination of coarse errors from navigation calculations.
The improvement of the accuracy of the determined position using the DGPS technique for several baselines should be based on a weighting model. In this paper, a weighting model for determining the positioning accuracy for the DGPS technique based on the solution of 3 independent baselines is presented. The weighting of measurements in the case of the differential DGPS technique was determined on the basis of three mathematical models taking into account the length of the baseline, the vector length error and the number of tracked GPS satellites. The scheme for determining the resultant accuracy of the aircraft position using the DGPS technique based on 3 independent baselines can be written as follows: (a) criterion I-weighting as a function of baseline (vector) length: where (dB, dL, dh)-the resultant position accuracy of the aircraft in ellipsoidal BLh coordinates; (dB 1A , dL 1A , dh 1A )-the errors of the aircraft DGPS position in ellipsoidal BLh coordinates determined for the baseline 1A, relative to the RTK-OTF (Real Time Kinematic-On The Fly) solution [4,43]; d 3A -the distance between the GNSS reference station ( criterion II-weighting as a function of the error ellipse of the point position: where m 1A -the vector length error for baseline 1A, where ns 1A -the number of tracked GPS satellites for the baseline 1A, ns 2A -the number of tracked GPS satellites for the baseline 2A, ns 3A -the number of tracked GPS satellites for the baseline 3A.

Research Test
The proposed calculation strategy for the mathematical Equations (6)- (11) was implemented and tested for the GPS kinematic data collected from an airborne experiment performed at the civil airport in Mielec, Poland. The test flight was performed over EPML  The aim of the test flight was to determine the basic parameters of the GNSS positioning quality for approach and landing procedures at Mielec airport. As part of the flight test, navigational calculations were also performed for the research method proposed in this paper. Navigational and observation data in RINEX format (Receiver Independent Exchange Format) from the on-board GNSS receiver mounted on Seneca Piper PA34-200T aircraft were used in the calculations. A Topcon HiperPro receiver was placed in the cockpit, which recorded GPS satellite data with an interval of 1 s. Virtual GPS reference stations generated in the positioning service POZGEO-D in ASG-EUPOS (Active Geodetic Network EUPOS) [44] were used as GNSS base stations. Table 1 gives the ellipsoidal BLh coordinates of the GNSS base stations. The location of the GNSS reference stations was selected to ensure that the VirA, VirB and VirC stations were evenly distributed along the flight path of the Seneca Piper PA34-200T aircraft.

Marker of GNSS Base Station Latitude (B) Longitude (L) Ellipsoidal Height (h)
VirA (  An illustrative sketch of the location of the VirA, VirB and VirC reference stations is shown in Figure 1. The aim of the test flight was to determine the basic parameters of the GNSS positioning quality for approach and landing procedures at Mielec airport. As part of the flight test, navigational calculations were also performed for the research method proposed in this paper. Navigational and observation data in RINEX format (Receiver Independent Exchange Format) from the on-board GNSS receiver mounted on Seneca Piper PA34-200T aircraft were used in the calculations. A Topcon HiperPro receiver was placed in the cockpit, which recorded GPS satellite data with an interval of 1 s. Virtual GPS reference stations generated in the positioning service POZGEO-D in ASG-EUPOS (Active Geodetic Network EUPOS) [44] were used as GNSS base stations. Table 1 gives the ellipsoidal BLh coordinates of the GNSS base stations. The location of the GNSS reference stations was selected to ensure that the VirA, VirB and VirC stations were evenly distributed along the flight path of the Seneca Piper PA34-200T aircraft. An illustrative sketch of the location of the VirA, VirB and VirC reference stations is shown in Figure 1.
The navigation calculations were divided into two stages. In the first stage, calculations for mathematical Equations (1)-(5) were performed. This stage was implemented in RTKLIB software [45]. RTKPOST application library was used in RTKLIB software, and the following calculation configuration was set up sequentially:  Table 1, -calculation interval: 1 s.
On this basis, three reports were obtained from the RTKPOST calculation mode. The reports contain information on the ellipsoidal coordinates of the plane, the time of measurement, the number of GPS satellites tracked and the average errors of the determined coordinates. In the next stage, all three reports were imported to the Scilab software [46], in which a navigation application for calculations for the scheme of mathematical Equations (6)-(11) was developed. The entire source code with the positioning algorithm for the three weighting criteria was written in Scilab. It should be additionally emphasised that the flight reference position of the Seneca Piper PA34-200T aircraft was also determined in the RTKPOST mode, separately for the three baselines. The results of the calculations are presented in Section 5.

Results
The proposed weighted average model presented in Equations (6)        From the point of view of accuracy analyses, the key parameter is the determination of the geometric coefficient PDOP (Position Dilution of Precision) [47,48]. The determination of the PDOP is of fundamental importance for the determination of position errors, also in the DGPS technique. Figure 5 shows the results of the PDOP parameter for the flight experiment in Mielec. It is worth noting that for the initial measurement epochs, the PDOP values are high, similar to the position errors in Figures 2-4. On the other hand, for the rest of measurement epochs, the PDOP values are equal to or less than two. Thus, it can be said that the higher the PDOP coefficient is, the lower the positioning accuracy is, which, in turn, translates into worse measurement conditions and a lower GNSS measurement quality [49]. This phenomenon is especially true for the h-component; it can be seen that the initial values of the position errors in  Based on the position errors obtained, it can be seen that the ellipsoidal height h has the lowest accuracy, while the L component has the highest accuracy.
From the point of view of accuracy analyses, the key parameter is the determination of the geometric coefficient PDOP (Position Dilution of Precision) [47,48]. The determination of the PDOP is of fundamental importance for the determination of position errors, also in the DGPS technique. Figure 5 shows the results of the PDOP parameter for the flight experiment in Mielec. It is worth noting that for the initial measurement epochs, the PDOP values are high, similar to the position errors in Figures 2-4. On the other hand, for the rest of measurement epochs, the PDOP values are equal to or less than two. Thus, it can be said that the higher the PDOP coefficient is, the lower the positioning accuracy is, which, in turn, translates into worse measurement conditions and a lower GNSS measurement quality [49]. This phenomenon is especially true for the h-component; it can be seen that the initial values of the position errors in The algorithm of the DGPS method (see Equation (1)) is characterised by strong correlation between the coded measurements recorded by two GNSS receivers for the same moment of time. This is of course related to the differentiation of code observations and a comparison of the individual code observations from the same satellite for the same epoch. It can be said that differential corrections in the DGPS method stimulate this correlation, which is of course influenced by the measurement process itself and the geometric construction in the DGPS method. This, in turn, ultimately translates into the determination of the aircraft position coordinates and positioning accuracy for a single baseline. This correlation is of course characteristic for the baselines 1A, 2A and 3A. Thus, comparing with each other the obtained position errors, shown in Figures 2-4, will actually show how strong this relationship is. Therefore, the relationship between the individual position errors for each baseline was determined, as recorded below: The algorithm of the DGPS method (see Equation (1)) is characterised by strong correlation between the coded measurements recorded by two GNSS receivers for the same moment of time. This is of course related to the differentiation of code observations and a comparison of the individual code observations from the same satellite for the same epoch. It can be said that differential corrections in the DGPS method stimulate this correlation, which is of course influenced by the measurement process itself and the geometric construction in the DGPS method. This, in turn, ultimately translates into the determination of the aircraft position coordinates and positioning accuracy for a single baseline.
This correlation is of course characteristic for the baselines 1A , 2A and 3A. Thus, comparing with each other the obtained position errors, shown in Figures 2-4, will actually show how strong this relationship is. Therefore, the relationship between the individual position errors for each baseline was determined, as recorded below:  Table 2 shows the results of determining the parameters of the relationship (rdB, rdL, rdh). On this basis, it can be seen that the relationship between the position errors between baseline 1A and 2A, as well as 1A and 3A is smaller than one. On the other hand, the relationship between baseline 2A and 3A is larger than one. Going back to Figure 1 and the distribution of the GNSS reference stations, it can be seen that the distance between the VirB and VirC stations is the shortest. It can be concluded that the smaller the distance is between stations, the greater the relationship is between position errors. For the developed calculation algorithm (Equations (6)- (11)), the key element is the determination of measurement weights for criteria I-III. The measurement weights were defined as a function of the baseline (vector) length, the vector length error and the number of GPS satellites tracked. The multiplicity of weighting functions is good, as it allows the user to determine the optimal measurement weighting strategy, which, in turn, translates into an improvement of the determined aircraft coordinates. Figures 6-8 show the measurement weighting values for criteria I-III. In the case of weighting in scheme I, the values of measurement weights are from 0.013 to 0.239 for baseline 1A, from 0.013 to 0.169 for baseline 2A and from 0.025 to 0.164 for baseline 3A. It is worth noting a simple relationship that follows directly from Figure 6. Namely, the more the measurement weight increases, the shorter the geometric distance between the on-board receiver and the GNSS reference station will be. This relationship of course occurs in reverse and is alternating.      the measurement weights vary inversely proportional to the squares of the mean errors, according to Equation (9). It can be seen that the highest value of measurement weights for criterion II reaches 0.125. Figure 8 presents the results of the measurement weights for criterion III. In this criterion, the measurement weights range from 0.111 to 0.250 for all the baselines, 1A, 2A and 3A. For weighting criterion III, the measurement weights are, in most cases, equal and the same. The following Figures 9-11 show the results of determining the resultant values of the position errors determined according to Equations (7), (9) and (11). For the B coordinate, the values of the resultant position errors are as follows:  Figure 7 shows the results of the measurement weights for criterion II. For this criterion, the measurement weights range from 0.003 to 0.125 for baseline 1A, from 0.003 to 0.100 for baseline 2A and from 0.010 to 0.096 for baseline 3A. For weighting criterion II, the measurement weights vary inversely proportional to the squares of the mean errors, according to Equation (9). It can be seen that the highest value of measurement weights for criterion II reaches 0.125. Figure 8 presents the results of the measurement weights for criterion III. In this criterion, the measurement weights range from 0.111 to 0.250 for all the baselines, 1A, 2A and 3A. For weighting criterion III, the measurement weights are, in most cases, equal and the same.     Analysing the results in Figures 9-11, it can be said that the positioning accuracy is highest for the L component and lowest for the h coordinate. This only confirms that the dispersion of position errors is the highest for the h component, while the lowest for the L coordinate.
In the discussion of the results of the resultant position errors, a statement was made that the weighting scheme for criteria I-III significantly reduced the parameter values   1  1  1  2  2  2  3  3  3   (  ,  ,  , 1  1  1  2  2  2  3  3  3   (  ,  ,  ,   In the discussion of the results of the resultant position errors, a statement was made that the weighting scheme for criteria I-III significantly reduced the parameter values (dB 1A , dL 1A , dh 1A , dB 2A , dL 2A , dh 2A , dB 3A , dL 3A , dh 3A ). It is, therefore, necessary to determine the percentage reductions for these position errors. For this purpose, the reductions in the resultant position errors (dB, dL, dh) relative to the parameters (dB 1A , dL 1A , dh 1A , dB 2A , dL 2A , dh 2A , dB 3A , dL 3A , dh 3A ), as shown in Equations (13)- (15), were calculated as follows:  Table 3. On the basis of the results presented in Table 3, it can be concluded that the reduction in position errors using Equations (6)-(11) is quite significant. It is only worth emphasising that in the case of the B component, one can speak of a reduction at the level of 55-94%, depending on the weighting criterion and the baseline. For the L component, the improvement in positioning accuracy ranges from 62% to 94%. For the ellipsoidal height h, on the other hand, the improvement in positioning accuracy is between 63 and 93%. As for the individual weighting criteria, criterion II gives the best results in terms of improvement in positioning accuracy. In contrast, the positioning accuracy improvement results are the worst for weighting criterion III compared to the values for criteria I and II. The results in Table 3 show how important the measurement weighting scheme is for the differential DGPS technique. In the analysis of the results, the parameter of the arithmetic mean for the resultant position errors from the individual criteria I-III was also determined, as recorded below:

Discussion
The scheme of the proposed mathematical model (Equations (6)   For the analysed parameters (dBw, dLw, dhw), the highest positioning accuracy is obviously for the L component, while the lowest is for the h coordinate.

Discussion
The scheme of the proposed mathematical model (Equations (6)  RTKPOST software was also used in the calculations to obtain results from the DGPS and RTK-OTF techniques. In addition, the algorithm scheme (Equations (6)-(16)) was made in Scilab software and the calculation results are shown in a graphical and tabular form. RTKPOST software was also used in the calculations to obtain results from the D and RTK-OTF techniques. In addition, the algorithm scheme (Equations (6)-(16)) made in Scilab software and the calculation results are shown in a graphical and tab form.  1  1  1  2  2  2  3  3  3   (  ,  ,  , Table 4 shows the results of determining the parameters of the relationship ( , , ) rdB rdL rdh . In this case, one can notice a regularity that basically all the parameter results are close to one. While in the flight test in Mielec, the parameter results ( , , ) rdB rdL rdh were varied, in the flight experiment in Dęblin, this is no longer the case.
Interestingly, the values of the rdL correlations are always larger than one, on the other hand, the parameters rdh are always smaller than one. This simply shows that the nature of the correlations of the ( , , ) rdB rdL rdh parameters in the flight experiment in Dęblin is higher due to equal distribution of GNSS base stations and shorter distances between them than in the Mielec experiment.

Parameter
Value Relation  Table 4 shows the results of determining the parameters of the relationship (rdB, rdL, rdh). In this case, one can notice a regularity that basically all the parameter results are close to one. While in the flight test in Mielec, the parameter results (rdB, rdL, rdh) were varied, in the flight experiment in Dęblin, this is no longer the case. Interestingly, the values of the rdL correlations are always larger than one, on the other hand, the parameters rdh are always smaller than one. This simply shows that the nature of the correlations of the (rdB, rdL, rdh) parameters in the flight experiment in Dęblin is higher due to equal distribution of GNSS base stations and shorter distances between them than in the Mielec experiment. It can be observed that the measurement weights determined as a function of the vector length have quite large values. In turn, the measurement weights as a function of the vector length error and the number of tracked GPS satellites are at a similar level, as in the Mielec flight experiment.         (7), (9) and (11). The resultant position errors for the flight experiment in Dęblin were as follows: For the B coordinate: The results of the resultant position errors show that the smallest scatter occurs for the L component and the largest occurs for the h component. Moreover, the highest positioning accuracy is observed for the weighting criterion II. Therefore, the results obtained from the flight test in Dęblin only confirm the conclusions drawn for the test results determined in the flight experiment in Mielec (see the descriptions to Figures 9-11).   Table 5, it can be seen that the positioning accuracy for weighting criterion II is the highest. On the other hand, for weighting criterion III, the worst results were obtained in terms of the improvement of positioning accuracy in relation to criteria I and II. For weighting criterion I, the improvement in accuracy ranges from 78% to 83%. For criterion II, the improvement in accuracy is 90-91%, while for criterion III, the improvement in DGPS positioning accuracy is 63-65%. Energies 2021, 21, x FOR PEER REVIEW 30of 35   Table 5 presents the parameter results (  ,  ,  , III  I  II  III  I  II  III PB PB PB PL PL PL Ph Ph Ph are at a similar level, as in the Mielec flight experiment. Based on Table 5, it can be seen that the positioning accuracy for weighting criterion II is the highest. On the other hand, for weighting criterion III, the worst results were obtained in terms of the improvement of positioning accuracy in relation to criteria   Table 5 presents the parameter results (  ,  ,  , III  I  II  III  I  II  III PB PB PB PL PL PL Ph Ph Ph are at a similar level, as in the Mielec flight experiment. Based on Table 5, it can be seen that the positioning accuracy for weighting criterion II is the highest. On the other hand, for weighting criterion III, the worst results were obtained in terms of the improvement of positioning accuracy in relation to criteria The final stage in the development of the test results for the flight experiment in Dęblin is the determination of the parameters (dBw, dLw, dhw), as shown in Figure 24. The values of the parameters (dBw, dLw, dhw) are as follows: -from −1.679 to 1.449 m for the parameter dBw, -from −1.270 to 0.127 m for the parameter dLw, -from −3.095 to 1.368 m for the parameter dhw.  The research problem shown in this paper concerned the accuracy of DGPS positioning in air navigation. The subject of accuracy research was analysed in scientific papers [6,13,[32][33][34]39]. In these works, the accuracy of DGPS or DGLONASS positioning was also determined in relation to the differential RTK technique. The developed test methodology from Equations (6)- (16) can significantly improve the accuracy of DGPS positioning in air As in the Mielec flight experiment and also in Dęblin experiment, the highest positioning accuracy is seen for the L component, while the lowest positioning accuracy is for the h coordinate.
The research problem shown in this paper concerned the accuracy of DGPS positioning in air navigation. The subject of accuracy research was analysed in scientific papers [6,13,[32][33][34]39]. In these works, the accuracy of DGPS or DGLONASS positioning was also determined in relation to the differential RTK technique. The developed test methodology from Equations (6)- (16) can significantly improve the accuracy of DGPS positioning in air navigation, especially for a solution using several GNSS reference stations. Moreover, the developed algorithm shows that it is still possible to improve the accuracy of DGPS positioning in air navigation, especially when using different configuration of weighting functions. The submitted work is part of a research trend of studies on improving the accuracy of DGPS positioning, similar to articles [6,13,[32][33][34]39].
A similar scope of research, as in the presented article, was conducted in works [50,51]. However, in these papers the research was carried out to improve the accuracy of DGPS positioning in static measurements for a network of GNSS referee stations. In contrast, our work shows an improvement of positioning accuracy in the kinematic mode, during flight experiments. Moreover, papers [50,51] present a mathematical model for the improvement of positioning accuracy using the weighting of measurements as a function of elevation angle to the satellite. On the other hand, our paper presents three weighting schemes, independent from each other, thanks to which it was possible to reduce the positioning errors. The approach presented in this paper and in papers [50,51] shows and emphasises how much of an important role in DGPS positioning model and weighting functions play, enabling us to improve the accuracy parameter.

Conclusions
This paper proposes a new research methodology for improving the accuracy of DGPS positioning using several GNSS reference stations. The paper proposes the use of a new weighting scheme for improving the determination of position errors in a differential DGPS technique. Namely, weighting functions in the form of the inverse of the vector length, the inverse of the square of the distance error and the inverse of the number of tracked GPS satellites for a single baseline were used in the calculations. On this basis, the equations of the linear combinations linking the weighting functions and the position errors determined for individual baselines were written. The weighting scheme for three independent baselines was used in the calculations. The proposed algorithm makes it possible to determine the resultant position errors for ellipsoidal BLh coordinates of the aircraft and significantly improve the accuracy of DGPS positioning.
The study used real observations and GPS navigation data from a Topcon HiperPro on-board receiver and the following three GNSS reference stations: VirA, VirB and VirC. Calculations were performed for the flight experiment in Mielec for the test flight carried out using a Seneca Piper PA34-200T aircraft. On the basis of the performed calculations, the following was found: in the case of the B component, there is a 55-94% error reduction; -for the L component, the improvement in positioning accuracy ranges from 62% to 94%; -for the ellipsoidal height h, the improvement in positioning accuracy ranges from 63% to 93%.
In addition, it is shown that weighting criterion II enables an improvement in DGPS positioning accuracy of 91-94%. In contrast, for weighting criterion I, the accuracy improvement is in the range of 83-89%, while for weighting criterion III it is 55-71%.
As part of the research, DGPS data were also processed for an airborne experiment carried out with a Cessna 172 aircraft at Dęblin airport. In the case of this flight experiment, the proposed research method allowed for an accuracy improvement from 78 to 83% for weighting criterion I, 90-91% for weighting criterion II and 63-65% for weighting criterion III.
The research methodology shown in this paper proved to be effective for improving the accuracy of DGPS positioning. Therefore, the assumptions and purpose of the work were fully achieved. Thus, the obtained research results are satisfactory for the performance of a differential DGPS technique in air navigation. Moreover, two independent flight experiments carried out in Mielec and Dęblin confirmed the high efficiency of the proposed research method.
The mathematical model described in this paper was aimed at improving positioning accuracy in the kinematic mode for air navigation. The presented model has an application in positioning with GNSS satellite technology; hence, the presented work is focused on this problem. The research is based on the fact that in air navigation newer and newer navigation solutions are sought to improve the quality of positioning, including the accuracy parameter. The presented navigation solutions only emphasise that the model of weighting coefficients can be an interesting research method in the area of DGPS positioning in air navigation. The model of weighting coefficients may turn out to be an effective solution for improving the accuracy of DGPS positioning and, thus, for improving the quality of the determined coordinates of an aircraft.

Funding: This research was funded by University Warmia and Mazury in Olsztyn.
Institutional Review Board Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.