Low-Cost Single-Frequency DGNSS/DBA Combined Positioning Research and Performance Evaluation

: In recent years, low-cost single-frequency GNSS receivers have been widely used in many ﬁelds such as mass navigation and deformation monitoring; however, due to the poor signal quality of low-cost patch antennae, it is difﬁcult for carrier phase real-time kinematic (RTK) technology to ﬁx the integer ambiguity. Differential GNSS (DGNSS) positioning with pseudorange can effectively meet the high robustness and reliability requirements for the submeter to the meter level positioning accuracy of UVA/vehicle/aerospace users. To improve the DGNSS positioning accuracy and reliability of low-cost single-frequency GNSS receivers in complex environments, we propose a differential barometric altimetry (DBA)-assisted DGNSS positioning algorithm, which solves the DGNSS observation equations jointly and rigorously with the Earth ellipsoidal constraint equations constructed by the DBA altitude. The DBA altitude accuracy at different baseline lengths was evaluated in detail, and the DGNSS positioning performance of the single-frequency low-cost u-blox receiver NEO-M8T with a patch antenna and DGNSS/DBA combined positioning performance with the BMP280 barometer was analyzed by several sets of static and dynamic experiments under different environments. The results show that the single-frequency NEO-M8T receiver with patch antenna DGNSS positioning accuracy is submeter level in the static environment and drops to meter level in the dynamic environment. GPS+BDS dual system has higher positioning accuracy than single GPS or single BDS. DGNSS/DBA combination has higher positioning accuracy than DGNSS, especially the root mean square error (RMSE) can be improved by 30% to 80% in the U direction and slightly improved in the N and E directions. This study can provide an effective solution reference for various applications of low-cost sensor fusion positioning in the mass consumer market. G.L.; X.D. G.L.; acquisition, S.W., X.D. G.L.


Introduction
With the continuous development and improvement of the Global Navigation Satellite System (GNSS), the number of visible satellites has been greatly increased, which effectively improves the positioning accuracy, reliability, and availability. However, in complex environments such as under trees, urban canyons, tunnels, deep mine pits, and indoors, etc., GNSS signal attenuation is severe, observation quality is poor, or the visible satellites are insufficient and there is the serious multipath effect. These drawbacks will lead to decreased positioning accuracy or even unable to locate, which greatly limit the availability and application of GNSS [1,2]. Combining GNSS with other multi-source sensor technologies to realize the complementary advantages of each system and improve the location-based used both indoors and outdoors. Low-cost barometer altimetry-assisted GNSS positioning navigation is also widely used in the field of aviation flights and smartphone navigation, etc. [22][23][24][25]. The basic principle of barometric altimetry is to use the physical phenomenon that the atmospheric pressure on the Earth's surface gradually decreases with increasing height, but, due to the irregular changes of atmospheric pressure, the altitude error directly calculated by a single barometer is as high as tens or even hundreds of meters, which cannot be used as a constraint to improve GNSS positioning accuracy. Similar to DGNSS, the accurate height can be obtained by barometric correction compensation or differential barometric altitude (DBA). That is, using the property that local atmospheric pressure changes are similar, a barometer is placed at the reference station and another barometer is used as a mobile station to determine its high-precision relative altitude by the differential equation. The altitude accuracy and reliability obtained by DBA mode are high, which can effectively constrain other technology to improve positioning and navigation accuracy.
The user altitude obtained by the DBA system was applied as a virtual satellite in the China Area Positioning System (CAPS), and the construction of independent earth ellipsoidal constraint equations by users' altitude can effectively solve the insufficient number of CAPS satellites and improve its 3D positioning accuracy and availability [26][27][28].
The DBA system also be applied to mobile cellular network positioning for accurately determining the user's height, and reduced 3D positioning to planar positioning which can obtain more desirable positioning accuracy [29]. Mobile cellular base stations can be used as DBA reference stations and transmitted the relevant atmosphere pressure data to the user side through mobile communication networks, achieving GNSS SPP/DBA combined 3D positioning with the altitude accuracy better than 1 m [30]. Inertial/barometric altitude can be fused to measure vertical velocity and height with velocity root mean squared error (RMSE) between 0.04 to 0.24 m/s and RMSE in height between 5 to 68 cm [31]. In addition, a barometer installed on a wearable device can measure vital signs such as blood pressure by detecting the position and orientation of the human body, thus providing a better telemedicine solution for precision medicine [32,33].
At present, to our knowledge, little research has been reported on the DGNSS/DBA combined positioning with low-cost GNSS receivers and a patch antenna, and there is also a lack of research and analysis on the theoretical methods and application effects of DBA, which is worth further study. In this study, we firstly propose a DGNSS/DBA combined positioning algorithm. Second, the DBA altitude accuracy at different baseline lengths is evaluated in detail. Then, the DGNSS performance of single-frequency low-cost NEO-M8T receiver and the accuracy and reliability of DGNSS/DBA combined positioning with low-cost BMP280 barometer are fully evaluated through actual measurement data.
This manuscript is organized as follows: Section 1 is the introduction. Section 2 is the mathematical model of DGNSS/DBA combined positioning, which contains Section 2.1 about DGNSS positioning observation equations; Section 2.2 about the principle of the DBA system; and Section 2.3 on the DGNSS/DBA combined positioning algorithm. Section 3 reports the experiment results, containing Section 3.1 that introduces experimental data; Section 3.2 about DBA altitude accuracy evaluation at different baseline lengths; Sections 3.3 and 3.4 on the DGNSS/DBA combined static and kinematic vehicle positioning performance evaluation for the single-frequency low-cost NEO-M8T receiver and BMP280 barometer, respectively. Section 4 presents the discussions. Section 5 is the conclusion.

DGNSS Positioning Observation Equation
GNSS receivers can simultaneously receive observation data such as code pseudorange, carrier phase, Doppler shift, and SNR. The raw observation values contain the receiver geometric position parameters, clock error as well as hardware delays parameters, and various error corrections, such as tropospheric errors and ionospheric errors. Therefore, Remote Sens. 2022, 14, 586 4 of 22 the raw code pseudorange observation equation between satellite s and receiver i can be expressed as: where P s i denotes the code pseudorange observation between satellite s and receiver i (Unit: m); ρ s i = (X s − X) 2 + (Y s − Y) 2 + (Z s − Z) 2 denotes the geometric distance between satellite s and receiver i at the moment of signal transmission, (X s , Y s , Z s ) and (X, Y, Z) are the satellite s and unknown receiver i antenna center position. c is the speed of light in vacuum; dt i denotes the receiver clock error; dt s denotes the satellite clock error; I s i is the ionospheric error; T s i is the tropospheric error; M s i is the multipath delay error; ε s i contains the code pseudorange measurement noise and other uncorrected errors.
For the short and medium baselines, the receiver clock error and satellite clock error are eliminated in the double-difference observation equation, the ionospheric and tropospheric errors can be neglected, and the DGNSS observation equation can be simplified and expressed as [34]: ∇∆P where ∇∆ denotes the double-difference operator; i and j denotes the reference station receiver and mobile receiver; s 1 and s k denote the reference and nonreference satellite, respectively.

Principle of the DBA System
Single-barometer altimetry affected by atmospheric temperature, seasonal changes, and other factors drifts up to several tens of meters within a day, with poor stability and reliability [35]. Due to the Earth's gravitational field, the space atmosphere pressure and height show a certain regular distribution. Except for the local strong convection zone, the trend of atmospheric pressure variation in the local range of several tens of kilometers shows the same physical characteristics, and the atmosphere is basically in hydrostatic equilibrium in the vertical direction. Usually, the atmospheric pressure is distributed more evenly in the horizontal direction, and the pressure difference is about 1 hPa at a distance of 100 km. Thus, the concept of "difference" can be extended to the field of barometric altimetry with the help of differential GNSS positioning, that is, by setting one (or several) barometric reference points, the barometric altimetry value of the mobile station can be corrected by the precise altitude of the reference station to compensate the influence of atmospheric physical environment changes on the altitude measurement results of the mobile station, thus improving the accuracy of the user altitude of the mobile station.
When the barometric reference station and the mobile station are within a few tens of kilometers, their latitude, gravitational acceleration, and water vapor factors have the same effect on the atmospheric pressure, so the above three errors can be neglected to obtain the simplified DBA formula [36]: where h is the altitude of the mobile station to be found, h 0 is the known altitude of the reference station, P 0 is the pressure of the reference station, P is the pressure of the mobile station, and t m is the average Celsius temperature between the reference station and the mobile station.

DGNSS/DBA Combined Positioning Algorithm
In this study, we propose a DGNSS/DBA combined positioning algorithm, in which the user's altitude obtained by the DBA system in Section 2.2 is used to construct the Earth ellipsoid constraint equation and solved rigorously as an independent observation jointly with the DGNSS observation equation, which is equivalent to adding a virtual satellite located at the center of the Earth [26]. Since the geodetic height is an independent variable Remote Sens. 2022, 14, 586 5 of 22 in GNSS coordinates, an approximate ellipsoid with the altitude h from the reference ellipsoid (WGS-84) can be constructed using the user geodetic height, as shown in Figure 1.

DGNSS/DBA Combined Positioning Algorithm
In this study, we propose a DGNSS/DBA combined positioning algorithm, in which the user's altitude obtained by the DBA system in Section 2.2 is used to construct the Earth ellipsoid constraint equation and solved rigorously as an independent observation jointly with the DGNSS observation equation, which is equivalent to adding a virtual satellite located at the center of the Earth [26]. Since the geodetic height is an independent variable in GNSS coordinates, an approximate ellipsoid with the altitude h from the reference ellipsoid (WGS-84) can be constructed using the user geodetic height, as shown in Figure  1. At this time, when the ground user's geodetic height is not very large, the observation equation after DGNSS/DBA combination can be expressed as: the symbols in the DGNSS observation equation in the first line of Equation (4) are the same as Equation (2).
( , , ) P X Y Z is the 3D coordinate of the ground user; a and b are the long and short semiaxes of the WGS-84 Earth reference ellipsoid, respectively. At this time, when the ground user's geodetic height is not very large, the observation equation after DGNSS/DBA combination can be expressed as: the symbols in the DGNSS observation equation in the first line of Equation (4) are the same as Equation (2). P(X, Y, Z) is the 3D coordinate of the ground user; a and b are the long and short semiaxes of the WGS-84 Earth reference ellipsoid, respectively. Since h is much smaller than the long and short semiaxes of the Earth reference ellipsoid, an approximate reference ellipsoid with a long semiaxis a + h and short semiaxis b + h is used instead without causing much bias [37]. To solve the Earth ellipsoid constraint equation in the second expression of Equation (4) by differential processing, the ellipsoid constraint equation is expanded in the user's approximate position (X 0 , Y 0 , Z 0 ) according to the Taylor series, and only the first-order term is retained, where the partial derivative of X is obtained as: after simplification, we get: similarly, taking partial derivatives of Y and Z yields: Remote Sens. 2022, 14, 586 6 of 22 let α = ∂h/∂X 0 , β = ∂h/∂Y 0 , γ = ∂h/∂Z 0 , that is, the Earth ellipsoidal constraint equation in Equation (4) is linearized at the approximate coordinates (X 0 , Y 0 , Z 0 ) to give: where dh = h −ĥ is the altitude residual, h is the altitude obtained by the DBA system, and h is the geodetic height obtained by the users' at the approximate position. dX, dY, dZ are positional corrections of the receiver antenna center. The detailed conversion process can be found in the original literature [38]. Similarly, the DGNSS observation equation in the first expression of Equation (4) is expanded by the Taylor series at the approximate position (X 0 , Y 0 , Z 0 ), omitting higher-order terms above the first order, and combined with Equation (9) to obtain the DGNSS/DBA combined positioning error equation: in Equation (10), the parameter estimated asx = dX dY dZ T contains three approximate position correction values; are the pseudorange double-difference directional cosine, respectively. l = L 1 · · · L i dh T are the observation value vectors. P = P DGNSS 0 0 P DBA is the DGNSS/DBA combined positioning weight matrix. P DGNSS = Q −1 DD is the a priori weight matrix of GNSS pseudorange double-difference observation equation. Q DD is the GNSS pseudorange double-difference observation values covariance matrix, the stochastic model of GNSS nondifferential observations adopts the sine trigonometric function elevation angle fixed-weight model [39]. According to the error propagation law, the GNSS relative positioning variance-covariance matrix can be expressed as Q DD [34]. The a priori weight P DBA of the DBA can be determined based on the results of the empirical evaluation in Section 3.2. We can solve Equation (10) using the single-epoch weighted least squares method, as Equation (11): where Qxx is the posterior covariance matrix of the parameterx. It can be found that the DBA altitude constraint is equivalent to adding a code pseudorange observation, and the condition number of the error equation coefficient matrix H becomes significantly better, the position dilution of precision (PDOP) value can be effectively reduced. However, the accuracy of the positioning solution is influenced by the accuracy of DBA altitude, i.e., if the accuracy of DBA altitude is better than DGNSS altitude, the improvement effect is obvious, otherwise, the positioning accuracy cannot be improved.

The Introduction of Experiment Data
In the experiment, the reference station included a high-precision geodetic GNSS receiver Trimble NET R9, a Trimble Choking 59,800 antenna, and a low-cost BMP280 barometer. The mobile station consisted of a single-frequency low-cost u-blox receiver NEO-M8T, a patch antenna, and two BMP280 barometers. The mobile station was also equipped with a geodetic receiver Trimble NET R9 for data acquisition, and the postprocessed kinematic (PPK) mode of commercial software Inertial Explorer 8.70 was used to process data to obtain high-precision 3D coordinate sequences as reference values. The nominal resolution of the low-cost BMP280 barometer was 0.01 mbar (0.1 m) and the data sampling rate was set at 1 Hz. As shown in Figure 2, to prevent the effect of crosswind on the barometric pressure, the barometer was placed inside a transparent plastic box with several small holes at the top of the box. All experimental data were recorded and postprocessed by a laptop computer. The height deviation of the BMP280 barometer from the GNSS antenna phase center was compensated by data preprocessing, and the DGNSS/DBA combined positioning analysis was performed by the self-written programs.

DBA Altitude Accuracy Evaluation at Different Baseline Lengths
The altitude accuracy and practical range obtained by the BMP280 barometer DBA model were evaluated through different baseline lengths outdoors. The experimental data were collected on 18 September 2020, and Figure 2 shows the hardware equipment of the reference station and mobile stations. The reference station was arranged on the observation pier on the roof of the Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences (APM, CAS). Figure 3 shows five different mobile station locations at baseline lengths of 0 m, 65 m, 2.6 km (shopping mall plaza), 6.0

DBA Altitude Accuracy Evaluation at Different Baseline Lengths
The altitude accuracy and practical range obtained by the BMP280 barometer DBA model were evaluated through different baseline lengths outdoors. The experimental data were collected on 18 September 2020, and Figure 2 shows the hardware equipment of the reference station and mobile stations. The reference station was arranged on the observation pier on the roof of the Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences (APM, CAS). Medical College), with a data acquisition time of about 30 min for each static point. Figure 4 and Table 1 show the time series and RMSE accuracy statistics of the altitude results of the DBA system under five groups of different baseline lengths.      In Table 1, DBA_min and DBA_max denote the minimum and maximum altitude  In Table 1, DBA_min and DBA_max denote the minimum and maximum altitude results obtained by the DBA model. Mean and STD denotes the average value of DBA altitude and standard deviation. As can be seen from Figure 4 and Table 1, the outdoor DBA altitude RMSE increased gradually with an increase of baseline length, and the DBA altitude RMSE was submeter level within the 2 km baseline length and did not exceed 2 m within the 10 km baseline length. This result can provide a priori information to determine the weight matrix P DBA in Equation (11) for the DGNSS/DBA combined positioning.

DGNSS/DBA Combined Static Positioning Results
This section mainly evaluates the static positioning performance of the single-frequency low-cost NEO-M8T receiver and BMP280 barometer DGNSS/DBA combined positioning algorithm. The experimental data were consistent with Section 3.2, and the 65 m and 6.0 km baseline length data were selected for processing and analysis. Two data processing modes, DGNSS and DGNSS/DBA, were set, and each mode was divided into single GPS, single BDS, and GPS+BDS dual systems by the satellite system.

Baseline Length 65 m
This experiment used 65 m short baseline static data and performed statistical analysis by setting the elevation mask angle from 10 to 40 degrees, and the prior error of the DBA system was σ DBA = 1.0 m. Table 2 shows the average number of visible satellites at different elevation mask angles, and Table 3 lists the average PDOP values for DGNSS and DGNSS/DBA modes at different elevation mask angles. As can be seen from Table 2, the average number of visible satellites was seven to eight for GPS and 12 to 14 for BDS at a low elevation mask angle of 10 or 20 degrees. With the increase of the elevation mask angle, the number of available satellites of both GPS and BDS systems decreased significantly, the satellite space geometry distribution became worse, and the PDOP value gradually increased. The number of GPS satellites was only three at the elevation mask angles of 40 degrees, and the user receiver could not be positioned at this time, while the number of visible satellites of BDS in the China region was larger with eight to 10 visible satellites at the elevation mask angles of 30 or 40 degrees. The GPS+BDS dual system significantly increased the number of visible satellites compared to the single system, which significantly improved the satellite geometry and reduced the PDOP value. As shown in Table 3, increasing a DBA observation was equivalent to adding a virtual satellite, which improved the satellite geometry distribution and reduced the PDOP value; and the reduction of PDOP value was more significant in the environment with a higher elevation mask angle. When the elevation mask angle was 40 degrees, three GPS satellites could not complete the positioning, and adding a DBA observation ensured the availability of user receiver positioning. Figures 5-8 show the north (N)/east (E)/up (U) direction deviation sequence of the two data processing modes at the elevation mask angle 10 to 40 degrees. Each mode included single GPS, single BDS, and a GPS+BDS dual system. Table 4 shows the RMSE values in the N/E/U directions for the two data processing modes at different elevation mask angles.   Table 4 shows the RMSE values in the N/E/U directions for t two data processing modes at different elevation mask angles.
(a) (b)          In this experiment, static data with 6.0 km baseline length were processed and analyzed, the elevation mask angle was 10 degrees, and the priori error of the DBA system was σ DBA = 1.5 m. Figure 9 shows the number of GNSS visible satellites and the PDOP value sequence for both DGNSS and DGNSS/DBA modes at the baseline length of 6.0 km. The RMSE accuracy statistics and corresponding deviation sequence in the N/E/U directions for the DGNSS and DGNSS/DBA mode are shown in Table 5 and Figure 10, respectively. 44.3%/12.2%/82.4%, and the RMSE of the GPS+BDS/DBA combination in the N/E/U dir tions was improved by 33.7%/0%/60.9%, respectively.
3.3.2. Baseline Length 6.0 km In this experiment, static data with 6.0 km baseline length were processed and a lyzed, the elevation mask angle was 10 degrees, and the priori error of the DBA syst was DBA σ =1.5 m . Figure 9 shows the number of GNSS visible satellites and the PD value sequence for both DGNSS and DGNSS/DBA modes at the baseline length of 6.0 k The RMSE accuracy statistics and corresponding deviation sequence in the N/E/U dir tions for the DGNSS and DGNSS/DBA mode are shown in Table 5 and Figure 10, resp tively.
(a) (b)     Table 5 and Figures 9 and 10 show that the PDOP value of BDS was smaller than GPS in the China region due to a large number of observable satellites. The RMSE of single BDS in the N/E/U directions was 0.61/2.27/2.64 m, which was better than that of single GPS in DGNSS mode (2.23/1.12/4.41 m). The RMSE of the GPS+BDS dual system in the N/E/U directions was 0.60/1.01/2.22 m with higher positioning accuracy than single system. Compared to DGNSS mode, the RMSE in the U direction of DGNSS/DBA mode reduced by 0.5 m and 0.1 m for single GPS and single BDS, and there was also some improvement in the N and E directions. However, the GPS+BDS dual system did not improve significantly, due to the higher DGNSS accuracy and the lower DBA height accuracy at 6.0 km baseline length did not prove to be an obvious constraint.  Table 5 and Figures 9 and 10 show that the PDOP value of BDS was smaller than G in the China region due to a large number of observable satellites. The RMSE of sin BDS in the N/E/U directions was 0.61/2.27/2.64 m, which was better than that of sin GPS in DGNSS mode (2.23/1.12/4.41 m). The RMSE of the GPS+BDS dual system in N/E/U directions was 0.60/1.01/2.22 m with higher positioning accuracy than single s tem. Compared to DGNSS mode, the RMSE in the U direction of DGNSS/DBA mode duced by 0.5 m and 0.1 m for single GPS and single BDS, and there was also some provement in the N and E directions. However, the GPS+BDS dual system did not prove significantly, due to the higher DGNSS accuracy and the lower DBA height ac racy at 6.0 km baseline length did not prove to be an obvious constraint.

DGNSS/DBA Combined Kinematic Vehicle Positioning Results
In this section, two kinematic vehicle experiments were designed to evaluate the namic positioning performance of the single-frequency low-cost NEO-M8T receiver a the BMP280 barometer DGNSS/DBA combined positioning algorithm.
The site photo of the mobile station equipment of the kinematic vehicle experim is shown in Figure 11. The platform contained four multisystem dual-frequency geode GNSS antennas which were connected to Trimble Net R9 receivers to obtain the high p cision reference value. The single-frequency, low-cost u-blox NEO-M8T receiver, a pa antenna and a BMP280 barometer were laid on the roof of the car. Two sets of kinema vehicle data were collected, the first set was located in an open urban environment a the second set was in a complex urban environment.

DGNSS/DBA Combined Kinematic Vehicle Positioning Results
In this section, two kinematic vehicle experiments were designed to evaluate the dynamic positioning performance of the single-frequency low-cost NEO-M8T receiver and the BMP280 barometer DGNSS/DBA combined positioning algorithm.
The site photo of the mobile station equipment of the kinematic vehicle experiment is shown in Figure 11. The platform contained four multisystem dual-frequency geodetic GNSS antennas which were connected to Trimble Net R9 receivers to obtain the high precision reference value. The single-frequency, low-cost u-blox NEO-M8T receiver, a patch antenna and a BMP280 barometer were laid on the roof of the car. Two sets of kinematic vehicle data were collected, the first set was located in an open urban environment and the second set was in a complex urban environment.  Table 5 and Figures 9 and 10 show that the PDOP value of BDS was smaller than G in the China region due to a large number of observable satellites. The RMSE of sin BDS in the N/E/U directions was 0.61/2.27/2.64 m, which was better than that of sin GPS in DGNSS mode (2.23/1.12/4.41 m). The RMSE of the GPS+BDS dual system in t N/E/U directions was 0.60/1.01/2.22 m with higher positioning accuracy than single s tem. Compared to DGNSS mode, the RMSE in the U direction of DGNSS/DBA mode duced by 0.5 m and 0.1 m for single GPS and single BDS, and there was also some i provement in the N and E directions. However, the GPS+BDS dual system did not i prove significantly, due to the higher DGNSS accuracy and the lower DBA height acc racy at 6.0 km baseline length did not prove to be an obvious constraint.

DGNSS/DBA Combined Kinematic Vehicle Positioning Results
In this section, two kinematic vehicle experiments were designed to evaluate the d namic positioning performance of the single-frequency low-cost NEO-M8T receiver a the BMP280 barometer DGNSS/DBA combined positioning algorithm.
The site photo of the mobile station equipment of the kinematic vehicle experime is shown in Figure 11. The platform contained four multisystem dual-frequency geode GNSS antennas which were connected to Trimble Net R9 receivers to obtain the high p cision reference value. The single-frequency, low-cost u-blox NEO-M8T receiver, a pa antenna and a BMP280 barometer were laid on the roof of the car. Two sets of kinema vehicle data were collected, the first set was located in an open urban environment a the second set was in a complex urban environment.

Open Urban Environment
This experiment was collected near the industrial park in Caidian District, Wuhan, China, on 20 November 2020, with a data duration of about 80 min. The area has an open urban environment with less shading and good data quality. Figure 12 shows the kinematic vehicle experiment test scene and test trajectory in the open urban environment. Figure 13 shows the number of GNSS visible satellites and the sequence of PDOP values during the kinematic vehicle experiment.

Open Urban Environment
This experiment was collected near the industrial park in Caidian District, Wuhan, China, on 20 November 2020, with a data duration of about 80 min. The area has an open urban environment with less shading and good data quality. Figure 12 shows the kinematic vehicle experiment test scene and test trajectory in the open urban environment. Figure 13 shows the number of GNSS visible satellites and the sequence of PDOP values during the kinematic vehicle experiment. As can be seen in Figure 13, the number of observable satellites of GPS and BDS fluctuated greatly for the kinematic vehicle, resulting in a significant increase in PDOP values compared to a static environment. The addition of DBA observations significantly improved the satellite geometric spatial distribution and reduced the PDOP values. In the open urban environment kinematic vehicle experiment, the comparison of the low-cost BMP280 barometer DBA altitude with the Trimble NET R9 reference altitude is shown in Figure 14. It can be seen that the BMP280 barometer DBA altitude had a high consistency with the reference altitude with an RMSE of 2.10 m. This value can be used to set the a priori weight matrix in the DGNSS/DBA combined positioning process. As can be seen in Figure 13, the number of observable satellites of GPS and BDS fluctuated greatly for the kinematic vehicle, resulting in a significant increase in PDOP values compared to a static environment. The addition of DBA observations significantly improved the satellite geometric spatial distribution and reduced the PDOP values. In the open urban environment kinematic vehicle experiment, the comparison of the low-cost BMP280 barometer DBA altitude with the Trimble NET R9 reference altitude is shown in Figure 14. It can be seen that the BMP280 barometer DBA altitude had a high consistency with the reference altitude with an RMSE of 2.10 m. This value can be used to set the a priori weight matrix in the DGNSS/DBA combined positioning process.

Open Urban Environment
This experiment was collected near the industrial park in Caidian District, Wuhan China, on 20 November 2020, with a data duration of about 80 min. The area has an open urban environment with less shading and good data quality. Figure 12 shows the kine matic vehicle experiment test scene and test trajectory in the open urban environment Figure 13 shows the number of GNSS visible satellites and the sequence of PDOP values during the kinematic vehicle experiment. As can be seen in Figure 13, the number of observable satellites of GPS and BDS fluctuated greatly for the kinematic vehicle, resulting in a significant increase in PDOP values compared to a static environment. The addition of DBA observations significantly im proved the satellite geometric spatial distribution and reduced the PDOP values. In the open urban environment kinematic vehicle experiment, the comparison of the low-cos BMP280 barometer DBA altitude with the Trimble NET R9 reference altitude is shown in Figure 14. It can be seen that the BMP280 barometer DBA altitude had a high consistency The RMSE statistical results and deviation sequence in the N/E/U directions for both DGNSS and DGNSS/DBA modes are given in Table 6 and Figure 15, respectively. The RMSE values of single GPS and single BDS in the N/E/U directions were 1.20/1.32/3.18 m and 1.58/2.11/5.02 m for DGNSS mode, respectively. The GPS/DBA, BDS/DBA, and GPS+BDS/DBA modes improved the RMSE by 40% to 60% in the U direction and increased slightly in the N and E directions, and the DBA altitude showed a good constraint effect.   The RMSE statistical results and deviation sequence in the N/E/U directions for both DGNSS and DGNSS/DBA modes are given in Table 6 and Figure 15, respectively. The RMSE values of single GPS and single BDS in the N/E/U directions were 1.20/1.32/3.18 m and 1.58/2.11/5.02 m for DGNSS mode, respectively. The GPS/DBA, BDS/DBA, and GPS+BDS/DBA modes improved the RMSE by 40% to 60% in the U direction and increased slightly in the N and E directions, and the DBA altitude showed a good constraint effect.  The RMSE statistical results and deviation sequence in the N/E/U directions for both DGNSS and DGNSS/DBA modes are given in Table 6 and Figure 15, respectively. The RMSE values of single GPS and single BDS in the N/E/U directions were 1.20/1.32/3.18 m and 1.58/2.11/5.02 m for DGNSS mode, respectively. The GPS/DBA, BDS/DBA, and GPS+BDS/DBA modes improved the RMSE by 40% to 60% in the U direction and increased slightly in the N and E directions, and the DBA altitude showed a good constraint effect.

Complex Urban Environment
This experiment was collected in Wuchang District, Wuhan, China, on November 20, 2020, with a valid data duration of about 1 h. The area is a complex urban environment, and the occlusion is relatively serious. Figure 16 shows kinematic vehicle experiment test scene and test trajectory in the complex urban environment. Figure 17 shows the number of visible GNSS satellites and the sequence of PDOP values during the kinematic vehicle experiment, and it can be seen that compared with the open urban environment, the number of GNSS visible satellites of the kinematic vehicle in the urban environment was significantly lower and the PDOP values became larger.

Complex Urban Environment
This experiment was collected in Wuchang District, Wuhan, China, on November 20, 2020, with a valid data duration of about 1 h. The area is a complex urban environment, and the occlusion is relatively serious. Figure 16 shows kinematic vehicle experiment test scene and test trajectory in the complex urban environment. Figure 17 shows the number of visible GNSS satellites and the sequence of PDOP values during the kinematic vehicle experiment, and it can be seen that compared with the open urban environment, the number of GNSS visible satellites of the kinematic vehicle in the urban environment was significantly lower and the PDOP values became larger.  Figure 18 shows the DBA altitude results of the low-cost BMP280 barometer during the kinematic vehicle experiment in the complex urban environment compared to the Trimble Net R9 reference altitude. It can be seen that the DBA altitude in the complex urban environment was also very consistent with the reference altitude with an RMSE of 2.19 m, which is approximately the same as the RMSE results of the DBA altitude in the open urban environment in Section 3.4.1. It indicates that the DBA altitude accuracy was reliable and stable in different environments, and could assist GNSS to improve the positioning accuracy in the complex urban environment.

Complex Urban Environment
This experiment was collected in Wuchang District, Wuhan, China, on November 2020, with a valid data duration of about 1 h. The area is a complex urban environme and the occlusion is relatively serious. Figure 16 shows kinematic vehicle experiment t scene and test trajectory in the complex urban environment. Figure 17 shows the numb of visible GNSS satellites and the sequence of PDOP values during the kinematic vehi experiment, and it can be seen that compared with the open urban environment, the nu ber of GNSS visible satellites of the kinematic vehicle in the urban environment was s nificantly lower and the PDOP values became larger.  Figure 18 shows the DBA altitude results of the low-cost BMP280 barometer duri the kinematic vehicle experiment in the complex urban environment compared to t Trimble Net R9 reference altitude. It can be seen that the DBA altitude in the comp urban environment was also very consistent with the reference altitude with an RMSE The RMSE statistics and deviation sequence of the single-frequency, low-cost u-blox NEO-M8T and the geodetic Trimble Net R9 receiver for both DGNSS and DGNSS/DBA modes in the N/E/U directions are given in Table 7 and Figures 19 and 20. The RMSE of the low-cost u-blox NEO-M8T receiver with single GPS and single BDS in DGNSS mode were 4.33/4.69/8.35 m and 4.91/6.91/19.48 m in the N/E/U directions, respectively, and the RMSE of the GPS+BDS dual system was 3.28/5.23/8.91 m. The RMSE statistics of the geodetic Trimble Net R9 receiver in N/E direction in the complex urban environment was maintained at submeter level; the U direction was relatively poor, with RMSE not exceeding 2.7 m. The difference of RMSE between the two GNSS receivers in 3D directions was determined by their hardware performance. Compared with DGNSS mode, the RMSE of single GPS, single BDS, and GPS+BDS dual system in DGNSS/DBA mode of low-cost NEO-M8T receiver slightly worsened in the N and E directions, while the RMSE in U direction could be improved by 50% to 80%, and this improvement ratio is higher than that in the open urban environment. Similarly, for the geodetic Trimble Net R9 receiver, the RMSE of the DGNSS/DBA combination remains the same in the N/E directions compared to DGNSS mode, and the RMSE in the U direction could be improved by 30% to 60%, which verifies the advantage of DBA altitude in assisting DGNSS positioning.  The RMSE statistics and deviation sequence of the single-frequency, low-cost u-blox NEO-M8T and the geodetic Trimble Net R9 receiver for both DGNSS and DGNSS/DBA modes in the N/E/U directions are given in Table 7 and Figures 19 and 20. The RMSE of the low-cost u-blox NEO-M8T receiver with single GPS and single BDS in DGNSS mode were 4.33/4.69/8.35 m and 4.91/6.91/19.48 m in the N/E/U directions, respectively, and the RMSE of the GPS+BDS dual system was 3.28/5.23/8.91 m. The RMSE statistics of the geodetic Trimble Net R9 receiver in N/E direction in the complex urban environment was maintained at submeter level; the U direction was relatively poor, with RMSE not exceeding 2.7 m. The difference of RMSE between the two GNSS receivers in 3D directions was determined by their hardware performance. Compared with DGNSS mode, the RMSE of single GPS, single BDS, and GPS+BDS dual system in DGNSS/DBA mode of low-cost NEO-M8T receiver slightly worsened in the N and E directions, while the RMSE in U direction could be improved by 50% to 80%, and this improvement ratio is higher than that in the open urban environment. Similarly, for the geodetic Trimble Net R9 receiver, the RMSE of the DGNSS/DBA combination remains the same in the N/E directions compared to DGNSS mode, and the RMSE in the U direction could be improved by 30% to 60%, which verifies the advantage of DBA altitude in assisting DGNSS positioning.  It can also be seen that when comparing Tables 3-5 for the static experiment with  Tables 6, 7 for the kinematic experiments, a low-cost receiver with only BDS signals provided lower RMSE values than GPS in static experiments, which was mainly due to BDS having more observation satellites and the PDOP value being lower than GPS at this time. In open and complex urban kinematic experiments, as shown in Figures 13 and 17, BDS had very large fluctuations for the number of observable satellites and a larger average PDOP value (3.02 and 4.10) than GPS (1.97 and 3.16) for both cases. This is the main reason that a low-cost receiver with only GPS signals provides lower RMSE values than with BDS in kinematic experiments.

Discussion
As the low-cost single-frequency GNSS receivers dominate most of the GNSS market [40], there is a strong interest in enhancing their accuracy. Low-cost DBA altitude plays a significant constraining role in improving the DGNSS positioning accuracy.
In the DBA altitude accuracy evaluation experiment, BMP280 barometers can achieve better than 2 m altitude accuracy within 10 km baseline lengths in static environments. The DBA altitude consistent with GNSS reference altitude in Figures 14 and 18 implies that it is also reliable and stable in complex environments. Low-cost single-frequency It can also be seen that when comparing Tables 3-5 for the static experiment with  Tables 6, 7 for the kinematic experiments, a low-cost receiver with only BDS signals provided lower RMSE values than GPS in static experiments, which was mainly due to BDS having more observation satellites and the PDOP value being lower than GPS at this time. In open and complex urban kinematic experiments, as shown in Figures 13 and 17, BDS had very large fluctuations for the number of observable satellites and a larger average PDOP value (3.02 and 4.10) than GPS (1.97 and 3.16) for both cases. This is the main reason that a low-cost receiver with only GPS signals provides lower RMSE values than with BDS in kinematic experiments.

Discussion
As the low-cost single-frequency GNSS receivers dominate most of the GNSS market [40], there is a strong interest in enhancing their accuracy. Low-cost DBA altitude plays a significant constraining role in improving the DGNSS positioning accuracy.
In the DBA altitude accuracy evaluation experiment, BMP280 barometers can achieve better than 2 m altitude accuracy within 10 km baseline lengths in static environments. The DBA altitude consistent with GNSS reference altitude in Figures 14 and 18 implies that it is also reliable and stable in complex environments. Low-cost single-frequency It can also be seen that when comparing Tables 3-5 for the static experiment with  Tables 6 and 7 for the kinematic experiments, a low-cost receiver with only BDS signals provided lower RMSE values than GPS in static experiments, which was mainly due to BDS having more observation satellites and the PDOP value being lower than GPS at this time. In open and complex urban kinematic experiments, as shown in Figures 13 and 17, BDS had very large fluctuations for the number of observable satellites and a larger average PDOP value (3.02 and 4.10) than GPS (1.97 and 3.16) for both cases. This is the main reason that a low-cost receiver with only GPS signals provides lower RMSE values than with BDS in kinematic experiments.

Discussion
As the low-cost single-frequency GNSS receivers dominate most of the GNSS market [40], there is a strong interest in enhancing their accuracy. Low-cost DBA altitude plays a significant constraining role in improving the DGNSS positioning accuracy.
In the DBA altitude accuracy evaluation experiment, BMP280 barometers can achieve better than 2 m altitude accuracy within 10 km baseline lengths in static environments.
The DBA altitude consistent with GNSS reference altitude in Figures 14 and 18 implies that it is also reliable and stable in complex environments. Low-cost single-frequency GNSS receivers with a patch antenna have become increasingly popular due to their lower and lower price. The DGNSS positioning accuracy of single-frequency low-cost GNSS receivers can still meet the submeter positioning accuracy needed by the general public in GNSS-friendly environments. However, the RMSE in the N/E/U directions are all at the meter-level in complex urban environments since low-cost GNSS receivers have poor observation quality, and the positioning accuracy of GPS+BDS dual system is significantly improved compared to single system. There is only a single-epoch resolution algorithm rather than a filtering algorithm is used in this study. In the future, with more and more satellites available for low-cost GNSS receivers and the use of multiple filtering algorithms, DGNSS positioning accuracy is expected to be further improved.
The DGNSS/DBA combined positioning can effectively improve the DGNSS positioning accuracy and meet the demand for real-time positioning applications. The Earth ellipsoid constraint equation constructed by the DBA altitude is equivalent to adding a virtual satellite located at the center of the Earth, effectively improving the spatial geometry structure of the observation satellite. The DGNSS/DBA combined positioning improves the positioning accuracy in the U direction by 30% to 80% compared with the DGNSS positioning, while the positioning accuracy in N and E directions also has a certain improvement effect.
Nowadays, most smartphones integrate both an inexpensive GNSS chip and a barometric pressure sensor. WADGNSS services [10] and a large number of meteorological stations [26] can provide correction information to users. We can achieve higher positioning accuracy without increasing hardware costs. The applications of low-cost DGNSS/DBA, such as indoor and outdoor seamless switching positioning, car navigation, emergency mapping, LBS, and rescue, etc. are likely to increase dramatically. Subsequently, based on the combined positioning of low-cost DGNSS/DBA, the positioning performance research by integrating other sensors, such as MEMS IMU and geomagnetic, etc., will be worth further investigation.

Conclusions
In this study, low-cost single-frequency DGNSS/DBA combined positioning research and performance evaluation was carried out. First, a DGNSS/DBA combined positioning model is proposed. The Earth ellipsoid constraint equation act as a virtual satellite observation at the center of the Earth, effectively improving the spatial geometry structure and PDOP value of the observation satellite. The low-cost BMP280 barometer DBA altitude accuracy is evaluated by different baseline lengths, which is better than the submeter level within 2 km and better than 2 m within 10 km baseline length. In both open and complex urban environment kinematic vehicle experiments, the DBA altitude accuracy is better than 2.20 m, which indicates that the DBA system has highly reliable and stability in different environments in local area.
The low-cost single-frequency NEO-M8T receiver with a patch antenna can achieve submeter level positioning accuracy for DGNSS positioning in the N/E directions and better than 1.5 m in the U direction in a short baseline static environment; as the baseline length increases, the DGNSS positioning accuracy gradually decreases. The positioning accuracy in a kinematic vehicle environment is significantly lower than in the static environment, and the RMSE in the N/E/U directions are all at the meter level in the complex urban environment, and the positioning accuracy of both GPS+BDS dual system is significantly improved compared to single system. The DGNSS/DBA combined positioning for low-cost NEO-M8T receiver and BMP280 barometer improves the positioning accuracy in the U direction by 30% to 80% compared with the DGNSS positioning, while the positioning accuracy in the N and E directions also has a certain improvement effect.