Next Article in Journal
Agnostic Envelope Linearization of Dynamically Supplied Power Amplifiers for Mobile Terminals
Next Article in Special Issue
Optimizing the Use of RTKLIB for Smartphone-Based GNSS Measurements
Previous Article in Journal
Comparison of Three Prototypes of PPG Sensors for Continual Real-Time Measurement in Weak Magnetic Field
Previous Article in Special Issue
A Combined Elevation Angle and C/N0 Weighting Method for GNSS PPP on Xiaomi MI8 Smartphones
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Instantaneous Best Integer Equivariant Position Estimation Using Google Pixel 4 Smartphones for Single- and Dual-Frequency, Multi-GNSS Short-Baseline RTK

1
National School of Surveying, University of Otago, 310 Castle Street, Dunedin 9016, New Zealand
2
Geomatic Innovation Research Group, Faculty of Built Environment and Surveying, Universiti Teknologi Malaysia, Johor Bahru 81310, Johor, Malaysia
3
Geoscience Australia, GPO Box 378, Canberra, ACT 2601, Australia
4
FrontierSI, Goods Shed, Village Street, Docklands, VIC 3008, Australia
5
Research School of Earth Sciences, Australian National University, Canberra, ACT 2601, Australia
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(10), 3772; https://doi.org/10.3390/s22103772
Submission received: 10 February 2022 / Revised: 9 May 2022 / Accepted: 9 May 2022 / Published: 16 May 2022
(This article belongs to the Special Issue Precise Positioning with Smartphones)

Abstract

:
High-precision global navigation satellite system (GNSS) positioning and navigation can be achieved with carrier-phase ambiguity resolution when the integer least squares (ILS) success rate (SR) is high. The users typically prefer the float solution under the scenario of having a low SR, and the ILS solution when the SR is high. The best integer equivariant (BIE) estimator is an alternative solution since it minimizes the mean squared errors (MSEs); hence, it will always be superior to both its float and ILS counterparts. There has been a recent development of GNSSs consisting of the Global Positioning System (GPS), Galileo, Quasi-Zenith Satellite System (QZSS), and the BeiDou Navigation Satellite System (BDS), which has made precise positioning with Android smartphones possible. Since smartphone tracking of GNSS signals is generally of poorer quality than with geodetic grade receivers and antennas, the ILS SR is typically less than one, resulting in the BIE estimator being the preferred carrier phase ambiguity resolution option. Therefore, in this contribution, we compare, for the first time, the BIE estimator to the ILS and float contenders while using GNSS data collected by Google Pixel 4 (GP4) smartphones for short-baseline real-time kinematic (RTK) positioning. It is demonstrated that the BIE estimator will always give a better RTK positioning performance than that of the ILS and float solutions while using both single- and dual-frequency smartphone GNSS observations. Lastly, with the same smartphone data, we show that BIE will always be superior to the float and ILS solutions in terms of the MSEs, regardless of whether the SR is at high, medium, or low levels.

1. Introduction

The key to high-precision global navigation satellite system (GNSS) positioning and navigation is carrier-phase integer ambiguity resolution. The ambiguity-fixed GNSS baseline, as obtained by integer least squares (ILS), is commonly expected to be superior to its float counterpart if the integer ambiguity success rate, i.e., the probability of correct integer estimation, is close to the maximum value of one. On the other hand, the float solution usually becomes the positioning preference when the success rate is too low. Alternatively, the best integer equivariant (BIE) estimator, as introduced by Teunissen [1], can be used since it always provides the optimal solution in terms of the mean squared errors (MSEs). Several studies have already investigated the extremely poor multipath suppression and linearly polarized patch antenna in the smartphone that is the foremost hurdle to achieving centimeter-level accurate positioning [2,3,4,5,6,7,8,9]. Since smartphones track GNSS signals with a poorer quality than geodetic-grade receivers and antennas, the ILS SR is typically different from one, and the BIE estimator would then be the preferred option. Laurichesse et al. [10] used undifferenced ambiguity resolution in precise point positioning (PPP) mode with smartphone measurements. Darugna et al. [9] and Warnant et al. [11] compared the positioning performance of different multi-GNSS positioning techniques (i.e., single-point positioning (SPP), differential GPS (DGPS), and real-time kinematic (RTK)) with different smartphone models. A much more recent study from Paziewski et al. [12] assessed the performance of several Android smartphones in relative positioning, whereby centimeter-level precision was achievable with fixed ambiguities. In this contribution, we study, for the first time, the single-baseline RTK positioning performance of the BIE estimator using smartphone GNSS data.
Verhagen and Teunissen [13] proved that this estimator is always optimal in terms of the MSE, while Wen et al. [14] demonstrated the use of the BIE estimator for GNSS precise point positioning (PPP). In Brack et al. [15] and Brack [16], a sequential BIE approach was developed. Subsequently, Teunissen [17] extended the theory of integer equivariant estimation by developing the principle of BIE estimation for the class of elliptically contoured distributions, while Odolinski and Teunissen [18] analyzed the BIE performance for low-cost, single- and dual-frequency, short- to long-baseline multi-GNSS RTK positioning, and they found that the BIE positions reveal a ‘star-like’ pattern when the ILS SRs are high. Odolinski and Teunissen [19] recently compared the RTK positioning performance of BIE for the multivariate normal and multivariate t-distribution.
The recent development of smartphone GNSS chipsets, such as Broadcom BCM47755 embedded, makes precise code-based positioning [20], PPP [21], and instantaneous, centimeter-level RTK positioning possible with Android-based smartphones [2,3,12,22,23,24]. In this contribution, we assess the BIE positioning performance using Google Pixel 4 (GP4) smartphones and compare the performance to that of the ILS and float estimators. The performance of the BIE estimator is assessed from the MSE perspective while using multi-GNSS smartphone data for an instantaneous (single-epoch) and single-baseline RTK model, while employing single- and dual-frequency observations.
This contribution is organized as follows: in Section 2, we describe the functional model that is employed for instantaneous single-baseline RTK positioning. In this section, we emphasize the differences among the float, ILS, and BIE estimators. In Section 3, we present the smartphone GNSS data and stochastic model settings. The setup configuration deployed in this study is deemed to be the best configuration by having the smartphones placed in an upright position, as examined and proven in Yong et al. [22]. Then, in Section 4, we provide an analysis of the single- and dual-frequency RTK positioning performance under different model strengths. We further analyze the performance of BIE and compare it to the ILS and float contenders from the MSE perspective. Lastly, a summary with conclusions is given in Section 5.

2. Instantaneous, Single- and Dual-Frequency, Multi-GNSS RTK Using the Float, Integer Least Squares, and Best Integer Equivariant Estimators

In this section, we describe the functional model employed for the single-baseline RTK model while tracking single- and dual-frequency multi-GNSS observations using GP4 smartphones. We also introduce the float, ILS, and BIE estimators.

2.1. Functional Model

We assume that the two GP4 smartphones track GPS, Galileo, QZSS, and BDS code and carrier-phase frequencies on two frequencies. We make use of broadcast ephemerides for satellite orbits and clocks. The relative ionospheric, tropospheric delays and satellite orbit errors can be assumed negligible, since we employ short baselines. The single-epoch (instantaneous) and linearized double-differenced (DD) system of observation equations reads as follows:
y = A a + B b ,
where y is the vector of DD carrier-phase and code observations, A is the design matrix of the DD integer ambiguities in vector a , and B corresponds to the design matrix of the real-valued baseline components b . We employ system-specific reference satellites when performing the between-satellite single-differences. We note that, if we would take a common reference satellite on the overlapping frequencies between the systems, it could further strengthen the model [25,26]. For the stochastic model, we use an elevation weighting sine function as employed in RTKlib v2.4.3 [27].

2.2. Float Estimation

To obtain the float solution, denoted with a ‘hat’ symbol, we estimate the ambiguities and baseline components as real valued parameters in a standard least-squares solution, obtaining
[ a , ^ b ^ ] T ,   [ Q a ^ a ^ Q a ^ b ^ Q b ^ a ^ Q b ^ b ^ ] ,
where a ^ , b ^ are the vectors of the float ambiguities and baseline components with dimensions n and p , respectively, Q a ^ a ^ , Q b ^ b ^ denote the corresponding variance covariance (VCV) matrices, and Q a ^ b ^ = Q b ^ a ^ T are the corresponding covariance matrices between the ambiguities and baseline components.

2.3. Integer Least-Squares Estimation

By using the float ambiguities a ^ in Equation (2) we can find the integer least-squares solution of the ambiguities, denoted with a ‘check’ symbol, by solving the following problem:
a ˇ = arg min | | a ^ z | | Q a ^ a ^ 2 ,
where | | . | | Q a ^ a ^ 2 = ( . ) T Q a ^ a ^ 1 ( . ) . This ILS problem is efficiently solved using the LAMBDA (least-squares ambiguity decorrelation adjustment) method [28], finally yielding the following fixed baseline solution:
b ˇ = b ^ Q b ^ a ^ Q a ^ a ^ 1 ( a ^ a ˇ ) .
The corresponding VCV matrix, provided that the uncertainty in a ˇ can be neglected, reads
Q b ˇ b ˇ = Q b ^ b ^ Q b ^ a ^ Q a ^ a ^ 1 Q a ^ b ^ .
The precision of the fixed ILS baseline b ˇ in Equation (4) is driven by the very precise phase data provided that the ILS SR is sufficiently high, while in the single-epoch case the float solution b ^ in Equation (2) is driven by the precision of the code data. This means, in the case that the ILS SR is sufficiently high, that the ILS solution is expected to have at least a two-order of magnitude positioning precision improvement compared to that of the float solution.

2.4. Best Integer Equivariant Estimation

When the ILS SR is low, the user will normally prefer the float solution b ^ in Equation (2) rather than the ILS solution b ˇ in Equation (4). The alternative is to use the BIE estimator to solve for the ambiguities [1]. Assuming normally distributed GNSS data, the BIE estimator we use is denoted with an ‘overline’ symbol and reads
a ¯ λ = z ϵ Θ a ^ λ z exp ( 1 2 | | a ^ z | | Q a ^ a ^ 2 ) z ϵ Θ a ^ λ exp ( 1 2 | | a ^ z | | Q a ^ a ^ 2 ) ,
Note in Equation (6) that the BIE solution is not always necessarily an integer as it is a weighted sum of integers. We also remark that, for the BIE estimator, no ratio test is needed [1].
The BIE baseline solution can then be derived as
b ˇ = b ^ Q b ^ a ^ Q a ^ a ^ 1 ( a ^ a ¯ λ ) ,
where a ˇ in Equation (4) is replaced by a ¯ λ in Equation (6). Note that the true BIE estimator involves an infinite weighted sum over the whole space of integers for z, which is computationally impractical. Hence, in Equation (6) we make use of a finite integer set Θ a ^ λ instead [29], which can be defined as follows:
|| a ^ z || Q a ^ a ^ 2 < λ 2 ,
where the threshold λ 2 can be determined from a central chi-squared distribution χ 2 with n degrees of freedom and a small significance level α = 10 9 . Note that, for very weak models, the number of candidates within this threshold in Equation (8) can reach several tens of thousands of candidates, whereas, for stronger models, at most a few candidates are usually obtained.

3. Google Pixel 4 Smartphone GNSS Data Collection

This section describes the short-baseline setup configurations of the GP4 smartphones while using (1) external antennas and (2) internal smartphone antennas. In this section, we also outline the stochastic model settings.

3.1. Setup Configuration with External and Internal Smartphone Antennas

Figure 1 depicts the external and internal antenna setups to evaluate the positioning performance of the GNSS data observed by GP4 smartphones. The GP4 smartphones are capable of tracking dual-frequency GPS L1 + L5, Galileo E1 + E5a, QZSS L1 + L5, and BDS B1 code and carrier-phase observations. The GP4 smartphones logged the GNSS measurements at a 1-s measurement interval via the Geo++ RINEX Logger vers. 2.1.6.
When evaluating the performance using external antennas, the GP4 smartphones are placed in two separate radiofrequency (RF) shielding boxes to avoid them from receiving the GNSS signals other than from the dedicated reradiating antenna (see Figure 1a–d). The GNSS signals are collected from two distinct active low-cost antennas, Swift GPS500, and then reradiated via a passive antenna inside the RF shielding boxes. The signal amplifier is connected between the rooftop antenna and reradiating antenna to mitigate the effect of signal attenuation over a 30 m LMR-400 flexible low loss communication coaxial cable. A similar setup was validated in Yong et al. [22] that benchmarked the smartphones with survey-grade receivers, to assure that no GNSS signal leakage was experienced in the RF shielding box. The duty-cycling settings of the GP4 smartphones were disabled during the experiment to obtain continuous carrier-phase observations [3].
In addition to the short baseline with external antennas, we assessed the positioning performance of the short baseline while using the smartphone internal antennas (see Figure 1e). Note that the built-in antennas of the GP4 smartphones have been found to be sensitive to poor quality GNSS signals and the surrounding environment [7,22], which means that we can expect larger multipath errors to be present for this setup.

3.2. Stochastic Model Settings

The stochastic model was determined by fitting the empirical 95% confidence interval/ellipses to the formal counterparts, as derived from the corresponding VCV matrices of the positions. The empirical VCV matrix was estimated from the positioning errors obtained by comparing the estimated positions to very precise benchmark coordinates, whereas the formal VCV matrix was obtained by the average of all single-epoch formal VCV matrices of the entire observation time span [30]. We used independent datasets to analyze and to formulate the stochastic model for the subsequent sections, and the stochastic model settings were determined for different elevation cut-off angles to formulate the most realistic stochastic models possible. By using realistic stochastic models, we could assure that we obtained the best possible ambiguity resolution and positioning performance results. The different elevation cut-off angles were used to mimic situations in urban canyon environments or when low-elevation multipath is present.
Table 1 depicts the range of the undifferenced and zenith-referenced standard deviations (STDs) utilized in the stochastic models, together with the observation span of the external and internal antenna setup configurations in Figure 1. Note that each GNSS and/or frequency had equal weighting in this article, similar to the RTKlib implementation [27].
Table 1 shows that the code STDs improved by a factor of approximately five when external antennas were used instead of the internal smartphone antennas that were more sensitive to multipath. For example, the code STD increased from a maximum of 1.4 m when using external antennas to approximately 6.0 m when internal antennas were used. The corresponding phase STDs increased from a maximum of 2 mm to 4 mm when using internal antennas. Note that the ambiguity resolution performance in the subsequent sections is driven by the number of satellites and frequencies used, as well as the stochastic model.

4. Instantaneous, Short-Baseline, Single- and Dual-Frequency RTK and BIE Positioning with Google Pixel 4 Smartphones

In this section, we investigate the BIE estimator and compare the performance to that of the commonly used ILS and float estimators. The benefit of using the BIE estimator is that the MSEs are always smaller than or at least as good as the float and ILS counterparts. For instance, in practice, when the ILS success rate is lower than the desired 100%, the user usually opts for the float solution when, in fact, BIE would be the preferred option.
The number of correctly fixed epochs, used below, was determined by the number of epochs where the estimated local east, north, and up coordinate errors were all below or equal to 0.05 m. The ILS success rate was then computed as follows:
P S E = #   of   correctly   fixed   epochs total   #   of   epochs × 100 % .
In the results below, we investigate the BIE performance for GP4 smartphones when using both external and internal smartphone antennas while collecting single- and dual-frequency multi-GNSS data.

4.1. BIE with External Antennas for Single-Frequency RTK

Figure 2 depicts the float (black), ILS (magenta), and BIE (green dots) horizontal RTK positioning errors using GP4 smartphones in a short-baseline RTK setup, while using low-cost external antennas. The positioning errors were determined by comparing the estimated positions to very precise benchmark coordinates. These benchmark coordinates were determined using geodetic GNSS receivers and antennas, a Kalman filter, and a multi-epoch model while assuming the ambiguities to be time-constant. Any phase center offsets and variations of the smartphones were neglected in this process. From top to bottom rows and left to right columns, we depict the results for various elevation cut-off angles resulting in ILS SRs of 11.5%, 54.9%, 79.7%, 94.3%, 99.8%, and 99.9%, respectively. We depict L1+ E1 + L1 + B1 GPS + Galileo + QZSS + BDS results, where the zoom-in windows show at least a two-order-of-magnitude improvement when going from ambiguity float and incorrectly fixed ILS solutions to that of the correctly fixed ILS positioning errors. Note that the float solutions, as depicted by black dots, become more precise as the model strength increases.
Figure 2 (top row and left column) shows that many of the incorrectly fixed ILS solutions (magenta dots) are at the meter level, and that the BIE solutions (green dots) resemble the float solutions (black dots underneath the green dots). When the ILS SRs increased, however, such as in the right column and second row, we could see BIE solutions starting to outperform their float counterparts, with a much larger density of BIE solutions with millimeter- to centimeter-level positioning precisions, as shown in the zoom-in windows. In the second row and right column, as well as in the third row, we can further see ILS solutions with larger positioning errors than BIE even though the ILS SRs ranged from 94.3% to 99.9%.
To further illustrate the optimal performance of the BIE estimator in terms of the positioning MSEs, we depict in Figure 3 the MSE ratios, with respect to the float MSEs, as a function of the ILS SRs. Note that the MSE is here the sum of the variances of the east, north, and up errors, since our estimated positions are unbiased. The float MSE ratio is equal to one and is depicted as a full blue line, whereas the ILS and BIE counterparts are depicted as dashed magenta and full green lines, respectively. Note that these MSE ratio results resemble those of Odolinski and Teunissen [19], albeit based on completely different datasets (smartphone vs. low-cost RTK receiver data).
Figure 3 shows, as expected, that the BIE MSE ratio is equal to that of the float solutions when the ILS SR is close to 0%, and that BIE is equal to ILS when the ILS SR is close to 100%. Most importantly and for all other cases, we can see that the BIE MSE ratios are smaller than those of the float and ILS solutions, respectively. This shows that using the BIE estimator on smartphone data for RTK positioning will give the optimal positioning performance, as measured by the MSEs.
Table 2 shows the percentage of the 3D position errors within a range of 0.05 m, 1 m, 2 m, and 4 m. The given percentages can provide a practical understanding of the distribution of the RTK positioning errors. The largest percentages for each scenario are shown in bold to distinguish the values from each other when rounded to two decimal places. The percentage of the ILS for the weakest model with the lowest ILS SR of 11.52% shows that the BIE solutions are similar to the float solutions. For the strongest model with the highest ILS SR of 99.92%, the results show that BIE resembled the ILS solutions. For all other scenarios, it becomes clear that BIE outperformed both the ILS and the float solutions in terms of not obtaining very large positioning errors, while also having a smaller likelihood than ILS of very small positioning errors (unless the ILS SR is very high). These results are similar to the cumulative distribution functions (CDFs) as discussed in Odolinski and Teunissen ([18]; Figure 4) and Verhagen and Teunissen ([13]; Figure 1).
To also illustrate the corresponding positioning precisions of the different estimators, Table 2 depicts the positioning standard deviations (STDs), the mean number (#) of satellites, and the employed elevation cut-off angles. We can see in Table 2 that the BIE and float solutions have similar STDs for the east, north, and up components when the ILS SR is 11.5% (with a slightly better performance for the BIE estimator), and that both solutions have STDs that are much better than their ILS counterparts. When the ILS SR increased to 54.9%, we can further see that the BIE estimator start to significantly outperform both the float and the ILS solutions, with STDs in east, north, and up of 2.344 m, 2.217 m, and 7.797 m, respectively. The corresponding STDs for the float and ILS solutions are 2.791 m, 2.605 m, and 9.455 m, and 2.661 m, 2.642 m, and 8.990 m in east, north and up, respectively. In other words, the BIE solutions have STDs that are up to more than 1 m smaller (in the up component) than the ILS counterparts.
Lastly, we can see in Table 2 that the BIE solutions have a better performance than their ILS counterparts when the ILS SRs reached values of 99.8% and 99.9%, with smaller positioning STDs by up to several centimeters in east, north, and up, as well as better performance by even several meters than the float solutions. This implies that when the ILS SR is different from the desirable 100%, the BIE estimator will indeed outperform the float and ILS estimators, and this is true even when smartphone GNSS data are used.

4.2. BIE with Internal Antennas for Dual-Frequency RTK

In this section, we investigate the corresponding BIE performance when the internal antennas of the smartphones are used. Since, with the smartphone internal antennas, the multipath errors are more significant than when external antennas are used [22], in this section, we use dual-frequency L1 + L5 GPS, E1 + E5a Galileo, L1 + L5 QZSS, and B1 BDS observations to further strengthen the model.
Figure 4 depicts, as in Figure 2, the float (black), ILS (magenta), and BIE (green dots) horizontal RTK positioning errors using GP4 smartphones in a short-baseline RTK setup, but while using the internal antennas of the smartphones. From top to bottom rows and left to right columns, we depict the results for various elevation cut-off angles resulting in ILS SRs of 9.6%, 53.9%, 72.1%, 84.0%, 91.9%, and 95.4%, respectively. The zoom-in windows show at least a two-order-of-magnitude improvement when going from ambiguity float and incorrectly fixed ILS solutions to that of the correctly fixed ILS positioning errors. Note again that, as the model become stronger, the float solutions, as depicted by black dots, become more precise. We also remark here that it is evident that the float and incorrectly fixed ILS solutions have a much poorer precision than in Figure 2. This degradation in precision when internal antennas are used is indeed due to their sensitive to multipath effects, where the code observations, which dominate the precision of the single-epoch float solutions, are more affected [31].
Figure 4 shows, similar to Figure 2, that many of the incorrectly fixed ILS solutions (magenta dots) have errors at the meter level, and that the BIE solutions (green dots) more or less resembles the float solutions (black dots underneath the green dots) when the models are weak (at the top row). When the ILS SRs increases, however, such as in the right column and second row, we can again see that BIE solutions start to outperform their float counterparts, with a much larger density of BIE solutions with millimeter- to centimeter-level positioning precisions as shown in the zoom-in windows. In the second row and right column, as well as in the third row, we can again see ILS solutions with larger positioning errors than BIE despite the fact that the ILS SRs ranged from 84.0% to 95.4%.
To show the above superior performance of the BIE estimator in a different way, Figure 5 illustrates the 95.4% ILS SR scenario (see Figure 4f), but now with each solution in a separate subplot. The figure shows that the BIE solutions are indeed superior to their float and ILS counterparts, with a better precision than both estimators and fewer large positioning errors than the ILS estimator.
To again illustrate the optimal BIE performance, we depict in Figure 6 the MSE ratios, with respect to the float MSEs, as a function of the ILS SRs. The float MSE ratio is depicted as a full blue line, whereas the ILS and BIE counterparts are depicted as dashed magenta and full green lines, respectively. Figure 6 shows, similar to Figure 3, that the BIE MSE ratio is equal to that of the float solutions when the ILS SR is close to 0%, and that BIE is close to ILS when the ILS SR is also close to 100% (i.e., 95.4%). Most importantly and for all other cases, we can again see that the BIE MSE ratios are smaller than those of the float and ILS solutions. This shows that using the BIE estimator on smartphone data for RTK positioning, even when the internal smartphone antennas are used, will give the optimal positioning performance.
Table 3 depicts the corresponding percentage of the 3D position errors within a range of 0.05 m, 1.0 m, 2.0 m, and 4.0 m, the positioning standard deviations (STDs), the mean number (#) of satellites, and the employed elevation cut-off angles when the internal smartphone antennas are used. We can again see that the percentages of the position errors are consistent with the CDFs in, e.g., Figure 4 of Odolinski and Teunissen [18]. Similarly, the BIE solutions have always smaller STDs for the east, north, and up components when the ILS SR is between 9.6% and 95.4%, with better STDs by more than 1 m to several tens of centimeters than the ILS solutions in east, north, and up, and better performance by several meters than the float solutions when the ILS SR is high. This implies again that when the ILS SR is different from the desirable 100%, the BIE estimator will indeed outperform the float and ILS estimators, and this is true even when smartphone GNSS data with internal antennas are used.

5. Discussion

In this contribution, we analyzed the best integer equivariant (BIE) estimator for real GNSS data collected by Google Pixel 4 (GP4) smartphones and antennas. We compared the instantaneous (single-epoch) positioning performance of BIE to the float and integer least squares (ILS) estimators that are commonly used when the ILS success rate (SR) is different from one and close to one, respectively. Radiofrequency (RF) shielding boxes and reradiating antennas were used to track GNSS signals from external low-cost antennas, consisting of L1 + L5 GPS, E1 + E5a Galileo, L1 + L5 QZSS, and B1 BDS code and carrier-phase observations. The short-baseline real-time kinematic (RTK) performance was also evaluated while using the GP4 internal smartphone antennas. We investigated the BIE performance both when single-frequency and dual-frequency measurements were employed for the combination of GPS + Galileo + QZSS + BDS. We showed that the BIE positioning performance was superior to that of the ILS and float estimators when the ILS SR is different from one. This was demonstrated to be true on the basis of real multi-GNSS data collected by the GP4 smartphones and antennas.
Our BIE performance evaluation consisted of comparing the estimated positions to very precise benchmark coordinates, and the optimality of the BIE estimator was further evaluated through its position mean squared errors (MSEs) and standard deviations (STDs). It was shown that the BIE performance resembles that of the float estimator when the ILS SR is very low and was similar to that of the ILS when the ILS SR is very high. For all other cases, we demonstrated that BIE outperformed both the float and the ILS estimators even when on the basis of real GP4 smartphone data while using external and internal smartphone antennas. Future studies could involve evaluating the GP4 smartphone BIE RTK positioning performance for longer baselines, when the relative atmospheric delays need to be estimated [19].

Author Contributions

C.Z.Y., conceptualization, methodology, formal analysis, visualization, and writing the manuscript; K.H., software development with BIE implementation into RTKlib; E.R. and S.M., conceptualization; R.O., conceptualization, methodology, formal analysis, investigation, and writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Geoscience Australia for the project 1002B: “Precise GNSS positioning with smartphones”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The broadcast ephemerides were used for satellite orbits and clocks. The Google Pixel 4 smartphone observation data are stored at the School of Surveying data facilities, the University of Otago, and can be made available upon request by contacting the corresponding author C.Z.Y. via email.

Acknowledgments

The fifth author was the recipient of a Geoscience Australia research grant on the project 1002B: “Precise GNSS positioning with smartphones”. We logged the GNSS data using a Geo++ RINEX logger from Android-based smartphones. All this support is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Teunissen, P.J.G. Theory of integer equivariant estimation with application to GNSS. J. Geod. 2003, 77, 402–410. [Google Scholar] [CrossRef] [Green Version]
  2. Pesyna, K.M.; Heath, R.W.; Humphreys, T.E. Centimeter positioning with a smartphone-quality GNSS antenna. In Proceedings of the 27th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS+ 2014), Tampa, FL, USA, 8–12 September 2014; pp. 1568–1577. [Google Scholar]
  3. Humphreys, T.E.; Murrian, M.; Van Diggelen, F.; Podshivalov, S.; Pesyna, K.M. On the feasibility of cm-accurate positioning via a smartphone’s antenna and GNSS chip. In Proceedings of the 2016 IEEE/ION Position, Location and Navigation Symposium, Savannah, GA, USA, 11–16 April 2016; pp. 232–242. [Google Scholar]
  4. Håkansson, M. Characterization of GNSS observations from a Nexus 9 Android tablet. GPS Solut. 2019, 23, 21. [Google Scholar] [CrossRef] [Green Version]
  5. Tomaštík, J.; Chudá, J.; Tunák, D.; Chudý, F.; Kardoš, M. Advances in smartphone positioning in forests: Dual-frequency receivers and raw GNSS data. For. Int. J. For. Res. 2021, 94, 292–310. [Google Scholar] [CrossRef]
  6. Aggrey, J.; Bisnath, S.; Naciri, N.; Shinghal, G.; Yang, S. Multi-GNSS precise point positioning with next-generation smartphone measurements. J. Spat. Sci. 2019, 65, 79–98. [Google Scholar] [CrossRef]
  7. Zhang, X.; Tao, X.; Zhu, F.; Shi, X.; Wang, F. Quality assessment of GNSS observations from an Android N smartphone and positioning performance analysis using time-differenced filtering approach. GPS Solut. 2018, 22, 70. [Google Scholar] [CrossRef]
  8. Guo, L.; Wang, F.; Sang, J.; Lin, X.; Gong, X.; Zhang, W. Characteristics Analysis of Raw Multi-GNSS Measurement from Xiaomi Mi 8 and Positioning Performance Improvement with L5/E5 Frequency in an Urban Environment. Remote Sens. 2020, 12, 744. [Google Scholar] [CrossRef] [Green Version]
  9. Darugna, F.; Wübbena, J.; Ito, A.; Wübbena, T.; Wübbena, G.; Schmitz, M. RTK and PPP-RTK Using Smartphones: From Short-Baseline to Long-Baseline Applications. In Proceedings of the 32nd International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2019), Miami, FL, USA, 16–20 September 2019; pp. 3932–3945. [Google Scholar]
  10. Laurichesse, D.; Rouch, C.; Marmet, F.X.; Pascaud, M. Smartphone applications for precise point positioning. In Proceedings of the 30th International Technical Meeting of the Satellite Division of The Institute of Navigation, Portland, OR, USA, 25–29 September 2017; pp. 171–187. [Google Scholar]
  11. Warnant, R.; Vyvere, D.; Van, L.; Warnant, Q. Positioning with single and dual frequency smartphones running Android 7 or later. In Proceedings of the 31st International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2018), Miami, FL, USA, 24–28 September 2018; pp. 284–303. [Google Scholar]
  12. Paziewski, J.; Fortunato, M.; Mazzoni, A.; Odolinski, R. An analysis of multi-GNSS observations tracked by recent Android smartphones and smartphone-only relative positioning results. Measurement 2021, 175, 109162. [Google Scholar] [CrossRef]
  13. Verhagen, S.; Teunissen, P.J.G. Performance comparison of the BIE estimator with the float and fixed GNSS ambiguity estimators. In A Window Future Geod; Springer: Berlin, Heidelberg, 2005; pp. 428–433. [Google Scholar]
  14. Wen, Z.; Henkel, P.; Brack, A.; Günther, C. Best integer equivariant estimation for precise point positioning. In Proceedings of the ELMAR-2012, Zadar, Croatia, 12–14 September 2012; pp. 279–282. [Google Scholar]
  15. Brack, A.; Henkel, P.; Günther, C. Sequential best integer-equivariant estimation for GNSS. Navigation 2014, 61, 149–158. [Google Scholar] [CrossRef]
  16. Brack, A. Partial Carrier-Phase Integer Ambiguity Resolution for High Accuracy GNSS Positioning. Ph.D. Thesis, Lehrstuhl fur Kommunikation und Navigation, Technische Universitat Munchen, München, Germany, 2019. [Google Scholar]
  17. Teunissen, P.J.G. Best integer equivariant estimation for elliptically contoured distributions. J. Geod. 2020, 94, 82. [Google Scholar] [CrossRef]
  18. Odolinski, R.; Teunissen, P.J.G. Best integer equivariant estimation: Performance analysis using real data collected by low-cost, single- and dual-frequency, multi-GNSS receivers for short- to long-baseline RTK positioning. J. Geod. 2020, 94, 91. [Google Scholar] [CrossRef]
  19. Odolinski, R.; Teunissen, P.J.G. Best integer equivariant position estimation for multi-GNSS RTK: A multivariate normal and t-distributed performance comparison. J. Geod. 2022, 96, 3. [Google Scholar] [CrossRef]
  20. Paziewski, J.; Sieradzki, R.; Baryla, R. Signal characterization and assessment of code GNSS positioning with low-power consumption smartphones. GPS Solut. 2019, 23, 98. [Google Scholar] [CrossRef] [Green Version]
  21. Purfürst, T. Evaluation of Static Autonomous GNSS Positioning Accuracy Using Single-, Dual-, and Tri-Frequency Smartphones in Forest Canopy Environments. Sensors 2022, 22, 1289. [Google Scholar] [CrossRef] [PubMed]
  22. Yong, C.Z.; Odolinski, R.; Zaminpardaz, S.; Moore, M.; Rubinov, E.; Er, J.; Denham, M. Instantaneous, Dual-Frequency, Multi-GNSS Precise RTK Positioning Using Google Pixel 4 and Samsung Galaxy S20 Smartphones for Zero and Short Baselines. Sensors 2021, 21, 8318. [Google Scholar] [CrossRef] [PubMed]
  23. Paziewski, J. Recent advances and perspectives for positioning and applications with smartphone GNSS observations. Meas. Sci. Technol. 2020, 31, 091001. [Google Scholar] [CrossRef]
  24. Hesselbarth, A.; Wanninger, L. Towards centimeter accurate positioning with smartphones. In Proceedings of the 2020 European Navigation Conference (ENC), Dresden, Germany, 23–24 November 2020; pp. 1–8. [Google Scholar]
  25. Odijk, D.; Teunissen, P.J.G. Characterization of between-receiver GPS-Galileo inter-system biases and their effect on mixed ambiguity resolution. GPS Solut. 2013, 17, 521–533. [Google Scholar] [CrossRef]
  26. Odolinski, R.; Teunissen, P.J.G.; Odijk, D. Combined BDS, Galileo, QZSS and GPS single-frequency RTK. GPS Solut. 2015, 19, 151–163. [Google Scholar] [CrossRef]
  27. Takasu, T.; Yasuda, A. Development of the low-cost RTK-GPS receiver with an open source program package RTKLIB. In Proceedings of the International Symposium on GPS/GNSS, Jeju, Korea, 4 November 2009; pp. 1–6. [Google Scholar]
  28. Teunissen, P.J.G. The least-squares ambiguity decorrelation adjustment: A method for fast GPS integer ambiguity estimation. J. Geod. 1995, 70, 65–82. [Google Scholar] [CrossRef]
  29. Teunissen, P.J.G. On the computation of the best integer equivariant estimator. Artif. Satell. 2005, 40, 161–171. [Google Scholar]
  30. Odolinski, R.; Teunissen, P.J.G.; Odijk, D. First combined COMPASS/BeiDou-2 and GPS positioning results in Australia. Part II: Single- and multiple-frequency single-baseline RTK positioning. J. Spat. Sci. 2014, 59, 25–46. [Google Scholar] [CrossRef]
  31. Khodabandeh, A.; Zaminpardaz, S.; Nadarajah, N. A study on multi-GNSS phase-only positioning. Meas. Sci. Technol. 2021, 32, 095005. [Google Scholar] [CrossRef]
Figure 1. GP4 smartphones in a short-baseline setup configuration with external antennas (a,b) inside an RF shielding box (c,d), and with the smartphone internal antennas (e) on the rooftop of the building of the School of Surveying in Dunedin, New Zealand.
Figure 1. GP4 smartphones in a short-baseline setup configuration with external antennas (a,b) inside an RF shielding box (c,d), and with the smartphone internal antennas (e) on the rooftop of the building of the School of Surveying in Dunedin, New Zealand.
Sensors 22 03772 g001
Figure 2. Horizontal (north/east) scatter plot of the multi-GNSS, single-frequency (L1 GPS + E1 Galileo + L1 QZSS + B1 BDS) GP4 data derived with the BIE (green dots), ILS (magenta dots), and ambiguity float (black dots) estimators for instantaneous RTK positioning with external antenna for a short baseline in Dunedin, New Zealand, based on 8 h of data (1 s measurement interval): (a) 11.5% ILS SR (38° cut-off); (b) 54.9% ILS SR (35° cut-off); (c) 79.7% ILS SR (32° cut-off); (d) 94.3% ILS SR (30° cut-off); (e) 99.8% ILS SR (20° cut-off); (f) 99.9% ILS SR (10° cut-off).
Figure 2. Horizontal (north/east) scatter plot of the multi-GNSS, single-frequency (L1 GPS + E1 Galileo + L1 QZSS + B1 BDS) GP4 data derived with the BIE (green dots), ILS (magenta dots), and ambiguity float (black dots) estimators for instantaneous RTK positioning with external antenna for a short baseline in Dunedin, New Zealand, based on 8 h of data (1 s measurement interval): (a) 11.5% ILS SR (38° cut-off); (b) 54.9% ILS SR (35° cut-off); (c) 79.7% ILS SR (32° cut-off); (d) 94.3% ILS SR (30° cut-off); (e) 99.8% ILS SR (20° cut-off); (f) 99.9% ILS SR (10° cut-off).
Sensors 22 03772 g002
Figure 3. MSE ratios of GP4 single-frequency, short-baseline, and instantaneous RTK positioning errors using external antennas, with BIE (green line), ILS (dashed magenta line), and ambiguity float (blue line), all versus the float solutions.
Figure 3. MSE ratios of GP4 single-frequency, short-baseline, and instantaneous RTK positioning errors using external antennas, with BIE (green line), ILS (dashed magenta line), and ambiguity float (blue line), all versus the float solutions.
Sensors 22 03772 g003
Figure 4. Horizontal (north/east) scatter plot of the multi-GNSS, dual-frequency (L1 + L5 GPS, E1 + E5a Galileo, L1 + L5 QZSS, and B1 BDS) GP4 data derived with the BIE (green dots), ILS (magenta dots), and ambiguity float (black dots) estimators for instantaneous RTK positioning with smartphone internal antennas for a short baseline in Dunedin, New Zealand, based on 6 h of data (1 s measurement interval): (a) 9.6% ILS SR (40° cut-off); (b) 53.9% ILS SR (30° cut-off); (c) 72.1% ILS SR (25° cut-off); (d) 84.0% ILS SR (20° cut-off); (e) 91.9% ILS SR (15° cut-off); (f) 95.4% ILS SR (10° cut-off).
Figure 4. Horizontal (north/east) scatter plot of the multi-GNSS, dual-frequency (L1 + L5 GPS, E1 + E5a Galileo, L1 + L5 QZSS, and B1 BDS) GP4 data derived with the BIE (green dots), ILS (magenta dots), and ambiguity float (black dots) estimators for instantaneous RTK positioning with smartphone internal antennas for a short baseline in Dunedin, New Zealand, based on 6 h of data (1 s measurement interval): (a) 9.6% ILS SR (40° cut-off); (b) 53.9% ILS SR (30° cut-off); (c) 72.1% ILS SR (25° cut-off); (d) 84.0% ILS SR (20° cut-off); (e) 91.9% ILS SR (15° cut-off); (f) 95.4% ILS SR (10° cut-off).
Sensors 22 03772 g004
Figure 5. The 95.4% ILS SR scenario (see Figure 4f): horizontal (north/east) scatter plot of the multi-GNSS, dual-frequency (L1 + L5 GPS, E1 + E5a Galileo, L1 + L5 QZSS, and B1 BDS) GP4 data derived with the ILS (left column), BIE (middle column), and ambiguity float (right column) estimators for instantaneous RTK positioning with smartphone internal antennas.
Figure 5. The 95.4% ILS SR scenario (see Figure 4f): horizontal (north/east) scatter plot of the multi-GNSS, dual-frequency (L1 + L5 GPS, E1 + E5a Galileo, L1 + L5 QZSS, and B1 BDS) GP4 data derived with the ILS (left column), BIE (middle column), and ambiguity float (right column) estimators for instantaneous RTK positioning with smartphone internal antennas.
Sensors 22 03772 g005
Figure 6. MSE ratios of the GP4 dual-frequency, short-baseline, and instantaneous RTK positioning errors using internal antennas, with BIE (green line), ILS (dashed magenta line), and ambiguity float (blue line), all versus the float solutions.
Figure 6. MSE ratios of the GP4 dual-frequency, short-baseline, and instantaneous RTK positioning errors using internal antennas, with BIE (green line), ILS (dashed magenta line), and ambiguity float (blue line), all versus the float solutions.
Sensors 22 03772 g006
Table 1. Stochastic model settings in terms of undifferenced and zenith-referenced STD range (minimum to maximum) of the code and phase observations. In the last column, the time span of the data to be analyzed is also depicted (in universal coordinated time, UTC).
Table 1. Stochastic model settings in terms of undifferenced and zenith-referenced STD range (minimum to maximum) of the code and phase observations. In the last column, the time span of the data to be analyzed is also depicted (in universal coordinated time, UTC).
Antennas UsedRange of Phase STDs (m)Range of Code STDs (m)DoY (Hours of Day), hh:mm:ss UTC
External antenna0.001–0.0021.198–1.413228–2021 (8 h), 13:35:00–21:34:59
Internal antenna0.003–0.0045.985–5.998344 and 345–2020 (6 h), 21:03:00–02:57:59
Table 2. Empirical STDs of the ILS, BIE, and float solutions for single-frequency and instantaneous RTK using GP4 smartphones, based on 8 h observations for a short-baseline external antenna data experiment (see Figure 2). Comparisons of the percentage of the 3D position errors within a range of 0.05 m, 1.0 m, 2.0 m, and 4.0 m are also given. Bold values show the largest percentages for each scenario and estimator.
Table 2. Empirical STDs of the ILS, BIE, and float solutions for single-frequency and instantaneous RTK using GP4 smartphones, based on 8 h observations for a short-baseline external antenna data experiment (see Figure 2). Comparisons of the percentage of the 3D position errors within a range of 0.05 m, 1.0 m, 2.0 m, and 4.0 m are also given. Bold values show the largest percentages for each scenario and estimator.
ModeILS Success Rate (%)Percentage of the 3D Position Errors (%)STD (m)Mean # of SatellitesElevation Cut-Off Angle
≤0.05 m<1.0 m<2.0 m<4.0 mENU
ILS11.5211.5246.1269.2791.183.5095.62415.0448.6
BIE0.0946.9072.7193.473.2385.28314.17038
Float0.0044.7071.3493.343.2825.33214.358
ILS54.9254.9272.7985.0896.162.6612.6428.99010.6
BIE5.4970.9288.2597.892.3442.2177.79735
Float0.0055.6582.1097.552.7912.6059.455
ILS79.7579.7587.3292.8698.011.8491.8836.03511.7
BIE29.0585.9794.5498.921.6341.5965.17232
Float0.0058.4584.8298.172.5072.4628.096
ILS94.2994.2996.3897.8599.351.1261.0243.07313.5
BIE67.8195.8698.2799.561.0140.8872.66730
Float0.0060.5986.4198.642.3702.2587.147
ILS99.8599.8599.9099.9499.990.1830.1920.42517.9
BIE99.6599.9099.9599.990.1700.1790.40020
Float0.0070.9893.4199.781.8461.8654.762
ILS99.9299.9299.9799.9899.990.0950.1470.44423.6
BIE99.9199.9799.9899.990.0930.1430.44010
Float0.0078.0396.8199.981.5561.6333.411
Table 3. Empirical STDs of the ILS, BIE, and float solutions for dual-frequency and instantaneous RTK using GP4 smartphones, based on 6 h observations for a short-baseline internal antenna experiment (see Figure 4). Comparisons of the percentage of the 3D position errors within a range of 0.05 m, 1.0 m, 2.0 m, and 4.0 m are also given. Bold values show the largest percentages for each scenario and estimator.
Table 3. Empirical STDs of the ILS, BIE, and float solutions for dual-frequency and instantaneous RTK using GP4 smartphones, based on 6 h observations for a short-baseline internal antenna experiment (see Figure 4). Comparisons of the percentage of the 3D position errors within a range of 0.05 m, 1.0 m, 2.0 m, and 4.0 m are also given. Bold values show the largest percentages for each scenario and estimator.
ModeILS Success Rate (%)Percentage of the 3D Position Errors (%)STD (m)Mean # of SatellitesElevation Cut-Off Angle
≤0.05 m<1.0 m<2.0 m<4.0 mENU
ILS9.629.6229.2844.4465.8610.993812.097229.70687.79
BIE0.5227.6645.2667.9510.672311.723128.562040
Float0.0024.1941.9866.4610.838711.852228.8412
ILS53.8953.8964.8972.6784.466.18216.109615.940211.87
BIE24.9860.4772.5086.155.87955.596214.844330
Float0.0031.8454.6179.306.44366.469316.6158
ILS72.1572.1578.8783.6390.864.60114.937210.836413.74
BIE46.2273.8982.5391.384.35704.53329.918625
Float0.0036.3260.8984.865.26555.710612.0913
ILS84.0184.0188.1091.3195.323.03213.05107.363015.49
BIE65.9386.1890.9595.722.86522.71736.808820
Float0.0040.1065.2988.754.26464.59099.8794
ILS91.9091.9094.2996.0198.051.99982.06444.117717.56
BIE80.0393.5295.9998.291.88191.89003.787715
Float0.0045.5771.1792.703.64344.07137.1847
ILS95.3995.3996.6497.6898.761.59671.38882.727819.81
BIE88.3796.3597.7599.071.48311.23862.448310
Float0.0047.8174.5594.413.38563.69806.1203
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yong, C.Z.; Harima, K.; Rubinov, E.; McClusky, S.; Odolinski, R. Instantaneous Best Integer Equivariant Position Estimation Using Google Pixel 4 Smartphones for Single- and Dual-Frequency, Multi-GNSS Short-Baseline RTK. Sensors 2022, 22, 3772. https://doi.org/10.3390/s22103772

AMA Style

Yong CZ, Harima K, Rubinov E, McClusky S, Odolinski R. Instantaneous Best Integer Equivariant Position Estimation Using Google Pixel 4 Smartphones for Single- and Dual-Frequency, Multi-GNSS Short-Baseline RTK. Sensors. 2022; 22(10):3772. https://doi.org/10.3390/s22103772

Chicago/Turabian Style

Yong, Chien Zheng, Ken Harima, Eldar Rubinov, Simon McClusky, and Robert Odolinski. 2022. "Instantaneous Best Integer Equivariant Position Estimation Using Google Pixel 4 Smartphones for Single- and Dual-Frequency, Multi-GNSS Short-Baseline RTK" Sensors 22, no. 10: 3772. https://doi.org/10.3390/s22103772

APA Style

Yong, C. Z., Harima, K., Rubinov, E., McClusky, S., & Odolinski, R. (2022). Instantaneous Best Integer Equivariant Position Estimation Using Google Pixel 4 Smartphones for Single- and Dual-Frequency, Multi-GNSS Short-Baseline RTK. Sensors, 22(10), 3772. https://doi.org/10.3390/s22103772

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop