Efficient Methods of Utilizing Multi-SBAS Corrections in Multi-GNSS Positioning

Various combining methods have been proposed to utilize multi-satellite-based augmentation system (SBAS) correction to provide accurate position in the global navigation satellite system (GNSS) receiver. However, the proposed methods have not been objectively compared and analyzed, making it difficult to know which ones are effective for multi-GNSS positioning. This paper presents efficient methods of combining multi-SBAS corrections in multi-GNSS positioning by comparing three methods: correction domain integration, measurement domain integration, and position domain integration. The performance of the three methods were analyzed through a covariance analysis that was expanded to multi-GNSS and multi-SBAS. Then, the results were verified by experiments using real measurements and corrections. Furthermore, implementation issues, such as computational complexity, availability, and flexibility, are analyzed. As a result, three methods had the same precision, but different complexity, availability, and flexibility. These results will be important guidelines to design, implement, and analyze navigation systems based on multi-GNSS with multi-SBAS corrections.


Introduction
Global navigation satellite systems (GNSS) have become a major navigation and positioning means in many areas, but the horizontal position accuracy of a standalone GNSS is known to be at least 5.0 m with a 95% confidence level [1]. To improve position accuracy, the receiver has been designed to receive multi-GNSS signals from multi-frequencies and multi-constellation GNSS [2][3][4]. In a multi-GNSS receiver, not only is the global positioning system (GPS) used, but also the Russian global navigation satellites system (GLONASS), the Beidou navigation satellite system (BDS), and Galileo are used to increase the accuracy and reliability of the position [2]. In addition, signal availability can be improved because the probability of receiving more than four satellite signals is improved by using multi-GNSS even in an urban canyon [3]. By using multi-frequencies, accuracy can be further improved by eliminating the ionospheric delay contained in the pseudo-range measurement [4].
Another way to improve the accuracy, precision, and integrity of a GNSS is the use of a differential GNSS (DGNSS). The accuracy and precision of a DGNSS is improved through the elimination of common errors, such as ephemeris error, satellite clock error, and ionospheric and tropospheric delay, using the corrections from reference stations located at known positions. Real-time kinematic (RTK) positioning is a form of DGNSS that uses carrier phase measurements and achieves centimeter-level precision [5,6].
A satellite-based augmentation system (SBAS) uses a geostationary satellite to provide the spatial corrections for a wide area. The fast, long-term, and ionospheric corrections are broadcasted to performance comparison and analysis for the various combination methods presented in this paper will be a guideline for the implementation of navigation systems using multiple SBAS and GNSS.
The rest of the paper is organized as follows: Section 2 briefly describes the SBAS operations and provides the covariance analysis. In Section 3, three integration methods (CDI, MDI, and PDI) for single GNSS positioning with multi-SBAS are presented and their covariance analysis is given. Section 4 expands the results of Section 3 to multi-GNSS positioning with multi-SBAS. Section 5 verifies the results using real measurements and corrections. Section 6 summarizes and discusses the properties of each integration method. The final section draws the conclusions.

Single GNSS Positioning and Error Analysis with Single SBAS Correction
A GNSS positioning procedure with SBAS correction was derived and the performance of the obtained position was evaluated using the covariance analysis. GNSS pseudo-range measurement was denoted as Equation (1), where subscript i represents the GNSS type, such as GPS, GLONASS, Galileo, BDS, and so on: If N i satellites are available in a GNSS, Ψ i ∈ R N i ×1 is a pseudo-range measurements vector, and r i ∈ R N i ×1 is the distance vector between the satellites and a receiver. cB i is the product of a receiver clock bias B ∈ R 1×1 and speed of light c, and cb i is the product of a satellite clock bias vector b i ∈ R N i ×1 and speed of light. E i ∈ R N i ×1 , I i ∈ R N i ×1 , and T i ∈ R N i ×1 indicate the ephemeris error vector, ionospheric delay vector, and tropospheric delay vector, respectively. In addition, a receiver measurement noise vector v i ∼ N(0, σ 2 i I N i ) ∈ R N i ×1 is the additive white Gaussian noise (AWGN) with zero mean and (σ i ) 2 variance, where I N i is an identity matrix with dimension N i .
As a wide-area differential GNSS, an SBAS provides three types of corrections using geostationary satellites to complement and supplement GNSS position accuracy and integrity. First, the fast correction can compensate a range measurement using two parameters: pseudo-range correction (PRC) and range rate correction (RRC). The fast correction parameters are included in Message types (MTs) 2-5 and they are transmitted every 6-60 s to the receiver. Second, the long-term correction eliminates the satellite position and satellite clock error. Because MTs 24 and 25 are transmitted in 120 s periods, they require a relatively longer time than fast correction. Third, ionospheric correction provides the amount of ionospheric delays in the service area using MTs 18 and 26. The receiver can compute the ionospheric delays of the ionospheric pierce point (IPP) by using the surrounding ionospheric grid points (IGPs). Figure 1 shows the positioning procedure using three corrections provided by an SBAS. Using the fast, long-term, and ionospheric corrections, cb i + E i + I i terms in Equation (1) can be compensated. In addition, using the tropospheric corrections in [15] (p. 253), which is not seen in Figure 1 because standard procedures are not defined, the T i term can also be corrected. Therefore, the computed SBAS corrections can be denoted using Equation (2), where the superscript j represents the type of SBAS, including WAAS, MSAS, SDCM, EGNOS, and KASS: Besides the three corrections, SBAS provides integrity and degradation information for each correction for the kth satellite in a GNSS variance of fast correction plus long-term correction (σ Besides the three corrections, SBAS provides integrity and degradation information for each correction for the kth satellite in a GNSS variance of fast correction plus long-term correction In this paper, the variance of SBAS correction for the k th satellite in a GNSS, as the sum of the fast, long-term, ionospheric, and tropospheric corrections, as shown in Equation (3), and was adopted as a weighting in the weighted least squares (WLS) positioning algorithm [16,17]: Because each correction of SBAS is uncorrelated, its covariance matrix is denoted as a diagonal matrix and can be written as Equation (4), where diag indicates a diagonal matrix and the diagonal terms 2 ( ) , 1, , can be computed using Equation (3): The corrected pseudo-range measurements were obtained using Equation (5) The covariance of receiver noise can be denoted as a matrix, which consists of the variance of each satellite: because the each channel's measurement noise in receiver can be different based on the elevation angle of satellite. Therefore, the covariance matrix of corrected measurements can be expressed as Equation (6): By linearizing at a linearization point, Equation (5) becomes Equation (7), where  In this paper, the variance of SBAS correction for the kth satellite in a GNSS, (σ j i ) 2 k , was computed as the sum of the fast, long-term, ionospheric, and tropospheric corrections, as shown in Equation (3), and was adopted as a weighting in the weighted least squares (WLS) positioning algorithm [16,17]: ( Because each correction of SBAS is uncorrelated, its covariance matrix is denoted as a diagonal matrix and can be written as Equation (4), where diag indicates a diagonal matrix and the diagonal The corrected pseudo-range measurements were obtained using Equation (5) by applying the SBAS corrections of Equation (2) to the pseudo-range of Equation (1). The error δη j i = η i − η j i + v i is the sum of receiver noise v i and the remaining terms after applying SBAS corrections η i − η j i ; without loss of generality, it is assumed as AWGN with zero mean vector and R j i covariance matrix [16] (p. 301).
The covariance of receiver noise can be denoted as a matrix, which consists of the variance of each satellite: (σ i ) 2 1 , (σ i ) 2 2 , . . . ,(σ i ) 2 N i because the each channel's measurement noise in receiver can be different based on the elevation angle of satellite. Therefore, the covariance matrix of corrected measurements can be expressed as Equation (6): By linearizing at a linearization point, Equation position, a few iterations are required to calculate this position [16]. The linearization point can be initialized by a known position that is close the user's current or previously calculated position.
Applying WLS with the weighting of Equation (6), the estimated navigation solution δx j i and its covariance are obtained as follows: In our analysis, H j i = H j is held because the satellites are very far away from the user and the effect of SBAS corrections to the line of sight vector is negligible. If the quality of the pseudo-range is improved by SBAS corrections, the covariance of the corrected pseudo-range measurement R j i will be decreased as compared to that of the uncorrected pseudo-range measurement R i . This implies that the precision of the estimated navigation solution δx j i is improved by applying SBAS corrections.

Single GNSS Positioning and Error Analysis with Multi-SBAS Correction
In this section, three methods (CDI, MDI, and PDI) to integrate multi-SBAS corrections for a single GNSS are presented and their performances are compared through a covariance analysis.

CDI
The concept of CDI is shown in Figure 2, where Ψ • is a pseudo-range measurement vector, are the weight vectors. In this section, subscript (•), instead of i, is used to indicate an unspecific GNSS system. The corresponding corrections are weighted and then summed to generate the CDI correction, as shown in Equation (10). And Equation (11) denotes the corrected pseudo-range.
Because the corrections from each SBAS are uncorrelated, the covariance of the corrected pseudo-range R CDI • can be written as By applying WLS to Equation (11), the navigation solution and covariance of CDI can be derived as Equations (13) and (14), respectively, where R CDI • is utilized as a weight matrix for WLS: Sensors 2020, 20, 256 6 of 16 In Equations (13) and (14), H CDI • becomes H • , as mentioned in Section 2, since the satellites are very far away from the user. Using an optimum weighting of Equation (15) [18], Equation (12) becomes Equation (16) by using the fact that the covariance matrix is a diagonal symmetric matrix.
, the navigation solution of Equation (13) can be re-written as that of Equation (17). Similarly, the covariance of CDI can be derived as Equation (18).
This result shows that the precision of the navigation solution may depend on the covariance that was expressed as Equation (18). And the covariance will be reduced by the integration of the inverse of each covariance of corrected measurements. As the simplest and best example, if all corrections of SBASs exhibit the same performance, i.e., R 1 /N SBAS , which means that the precision of the navigation solution is N SBAS times better than that of a single SBAS.
Sensors 2020, 20, 256 5 of 16 user's position, a few iterations are required to calculate this position [16]. The linearization point can be initialized by a known position that is close the user's current or previously calculated position.
Applying WLS with the weighting of Equation (6), the estimated navigation solution j i x δ and its covariance are obtained as follows:

Single GNSS Positioning and Error Analysis with Multi-SBAS Correction
In this section, three methods (CDI, MDI, and PDI) to integrate multi-SBAS corrections for a single GNSS are presented and their performances are compared through a covariance analysis.

CDI
The concept of CDI is shown in Figure 2, where • Ψ is a pseudo-range measurement vector,

MDI
The concept of MDI is shown in Figure 3. In MDI, N SBAS corrections from multi-SBAS were applied to a pseudo-range. Equation (19) is a linearized equation for each corrected pseudo-range and is simply denoted as Equation (20).
Sensors 2020, 20, 256 7 of 16 By applying WLS to Equation (20), the navigation solution and its covariance are obtained as Equations (21) and (22) A comparison of Equations (17) and (18) and Equations (21) and (22) clearly shows that the navigation solution and the covariance of MDI and CDI are identical.

PDI
In the PDI, shown in Figure 4, first, NSBAS navigation solutions were computed independently using the NSBAS-corrected pseudo-range. The procedure for finding each navigation solution is explained in Section 2. Then, NSBAS navigation solutions were weighted-summed to obtain the navigation solution of the PDI A comparison of Equations (17) and (18) and Equations (21) and (22) clearly shows that the navigation solution and the covariance of MDI and CDI are identical.
Equations (26) and (27) show that the navigation solution and covariance of PDI are same as those of CDI and MDI. This means that the three combination methods (CDI, MDI, and PDI) provide the same navigation solution with the same quality if the optimum weightings are used in case of a single GNSS and multi-SBAS.

Multi-GNSS Positioning and Error Analysis with Multi-SBAS Corrections
In this section, the results of single GNSS positioning with multi-SBAS corrections, explained in the previous section, are expanded to multi-GNSS positioning, including GPS, GLONASS, BDS, and Galileo. The CDI, MDI, and PDI were expanded to multi-GNSS and their performances were compared using a covariance analysis. Using Equation (25) as an optimum weighting [18], Equations (23) and (24) become Equations (26) and (27), respectively.
Equations (26) and (27) show that the navigation solution and covariance of PDI are same as those of CDI and MDI. This means that the three combination methods (CDI, MDI, and PDI) provide the same navigation solution with the same quality if the optimum weightings are used in case of a single GNSS and multi-SBAS.

Multi-GNSS Positioning and Error Analysis with Multi-SBAS Corrections
In this section, the results of single GNSS positioning with multi-SBAS corrections, explained in the previous section, are expanded to multi-GNSS positioning, including GPS, GLONASS, BDS, and Galileo. The CDI, MDI, and PDI were expanded to multi-GNSS and their performances were compared using a covariance analysis.

CDI
The concept of CDI for N GNSS GNSS and N SBAS SBAS is shown in Figure 5. By expanding Equation (10), the weighted sum of corrections for the ith GNSS is obtained as: Sensors 2020, 20, 256 9 of 16 Figure 5. CDI for multi-GNSS augmented with multi-SBAS.
Equations (29) and (30) are the corrected pseudo-range based on CDI and its covariance, respectively: By applying optimal weighting [18]

of Equation (31), Equation (30) become Equation (32)
By applying WLS to Equation (29) with the weighting of Equation (31), the estimated navigation solution and its covariance are derived as Equations (33) Equation (34) also shows that, by combining the covariance of correlated measurements, the Equations (29) and (30) are the corrected pseudo-range based on CDI and its covariance, respectively: By applying optimal weighting [18] of Equation (31), Equation (30) become Equation (32) By applying WLS to Equation (29) with the weighting of Equation (31), the estimated navigation solution and its covariance are derived as Equations (33) and (34), respectively. During the derivation, Equation (34) also shows that, by combining the covariance of correlated measurements, the overall combined covariance can be reduced; it is similar in form to Equation (18). For example, if the quality of each SBAS is equal (i.e., R 1 i = . . . = R N SBAS i , i = 1, . . . , N GNSS ), the covariance will be reduced to 1/N SBAS and then N SBAS times performance improvement can be expected. In addition, even if not realistic, if the constellation and number of satellites of each GNSS are the same, i.e., H 1 = H 2 = . . . = H N GNSS is satisfied, the covariance is reduced to 1/(N SBAS × N GNSS ). This implies a (N SBAS × N GNSS ) times improvement in the precision of the navigation solution is possible by combining the multi-GNSS and multi-SBAS in the CDI scheme.

MDI
The concept of MDI for N GNSS GNSS and N SBAS SBAS is shown in Figure 6. The augmented form of the linearized pseudo-range is denoted as:

MDI
The concept of MDI for NGNSS GNSS and NSBAS SBAS is shown in Figure 6. The augmented form of the linearized pseudo-range is denoted as: A comparison of Equations (36) and (37) and Equations (33) and (34) clearly shows that the navigation solution and the covariance of MDI and CDI are identical in multi-GNSS as well as in multi-SBAS.

MDI
The concept of MDI for NGNSS GNSS and NSBAS SBAS is shown in Figure 6. The augmented form of the linearized pseudo-range is denoted as: A comparison of Equations (36) and (37) and Equations (33) and (34) clearly shows that the navigation solution and the covariance of MDI and CDI are identical in multi-GNSS as well as in multi-SBAS.

PDI
The concept of PDI for N GNSS GNSS and N SBAS SBAS is shown in Figure 7. By applying WLS to all corrected pseudo-ranges, (N SBAS × N GNSS )-estimated navigation solutions and their covariances are given as Equations (38) and (39), respectively. cov(δx Sensors 2020, 20, 256 11 of 16 The concept of PDI for NGNSS GNSS and NSBAS SBAS is shown in Figure 7. By applying WLS to all corrected pseudo-ranges, (NSBAS × NGNSS)-estimated navigation solutions and their covariances are given as Equations (38) By expanding Equations (23) and (24), Equations (40) and (41) indicate the PDI navigation solution and covariance, respectively.
Sensors 2020, 20, 256 Equations (43) and (44) show that the navigation solution and covariance of PDI are the same as those of CDI and MDI. This means that the three methods (CDI, MDI, and PDI) give the same navigation solution with the same quality if the optimum weightings are used in cases of multi-GNSS and multi-SBAS.

Experiment and Verification Using Real Measurement and Corrections
To verify the analysis results, the navigation solutions of three methods (CDI, MDI, and PDI) that were calculated using real measurements and corrections were compared. In Chungbuk National University in Korea, the real measurements and corrections for 1000 s were collected using a commercial receiver (M8T, u-blox, Thalwil, Switzerland)) that supports receiving GPS L1 C/A, GLONASS G1 C/A, MSAS, and SDCM. Although Korea is not an included service area for MSAS and SDCM, these signals can be received and the given corrections are reasonable. Currently, MSAS provides corrections for only GPS and SDCM provides corrections for both GPS and GLONASS. Therefore, three types of corrected measurements were available, such as GPS corrected by MSAS, GPS corrected by SDCM, and GLONASS corrected by SDCM.
The commercial receiver that was used generated a local time that was estimated by considering the time offset among the different satellite systems. Because the pseudo-range of each GNSS is measured based on the local time, the user does not need to consider the time offset alignment for positioning using multi-GNSS. Therefore, the measurements and corrections from the commercial receiver could be applied to the three methods mentioned above without additional modification. In the collected real measurements, eight possible GPS measurements could be corrected by MSAS. Additionally, six GPS and five GLONASS measurements that could be corrected by SDCM were received.
Using the collected real measurements, navigation solutions were calculated using the three presented methods of CDI, MDI, and PDI. The experiment was performed with two scenarios. The first was where only GPS measurements and multiple corrections were used. The second was where both GPS and GLONASS measurements were used with two types of corrections. By comparing the results of the two scenarios, the efficiency of the expansion to the multi-GNSS that is presented in this paper was demonstrated.
In general, the reasonable covariance of GNSS and SBAS must be estimated to obtain a more integrated, precise, and accurate position. However, because this experiment focused on comparing the performance of the three integration methods in the same situation, each GNSS and SBAS covariance was initialized as 1.0 m. Subsequently, the positioning results of each method were compared as point-by-point to verify the error analysis results that the precision of the three methods would be equal. The positioning results for the vertical and horizontal errors were presented as east-north-up (ENU) coordinates, as depicted in Figures 8 and 9. Figure 8a,b depicts the horizontal and vertical errors for GPS, MSAS, and SDCM integration using CDI, MDI, and PDI for single GNSS and multi-SBAS that were mentioned in Section 3. To present the efficiency of multi-SBAS integration, results of the positioning using corrected GPS by MSAS are presented together as black circles. Table 1 presents the standard deviations for the east, north, and up axes that indicate the precision. As can be seen from the experiment results, when CDI, MDI, and PDI were applied, precision improved by approximately 0.01-0.07 m than in the case of GPS with MSAS. Moreover, in the scatter plot of Figure 8, red stars (CDI), green circles (MDI), and blue points (PDI) occur in the same positions. The cm-level precisions of CDI, MDI, and PDI were equal. These results verified that the integration of MSAS and SDCM is more precise than when only one SBAS was used. Additionally, the three integration methods had equal precisions. The experimental results are consistent with the error analysis results of Section 3. integrated, precise, and accurate position. However, because this experiment focused on comparing the performance of the three integration methods in the same situation, each GNSS and SBAS covariance was initialized as 1.0 m. Subsequently, the positioning results of each method were compared as point-by-point to verify the error analysis results that the precision of the three methods would be equal. The positioning results for the vertical and horizontal errors were presented as eastnorth-up (ENU) coordinates, as depicted in Figures 8 and 9.    Figure 9. Estimated position errors using CDI, MDI, and PDI for multi-GNSS and multi-SBAS Figure 8a,b depicts the horizontal and vertical errors for GPS, MSAS, and SDCM integration using CDI, MDI, and PDI for single GNSS and multi-SBAS that were mentioned in Section 3. To present the efficiency of multi-SBAS integration, results of the positioning using corrected GPS by MSAS are presented together as black circles. Table 1 presents the standard deviations for the east, north, and up axes that indicate the precision. As can be seen from the experiment results, when CDI, MDI, and PDI were applied, precision improved by approximately 0.01-0.07 m than in the case of GPS with MSAS. Moreover, in the scatter plot of Figure 8, red stars (CDI), green circles (MDI), and blue points (PDI) occur in the same positions. The cm-level precisions of CDI, MDI, and PDI were equal. These results verified that the integration of MSAS and SDCM is more precise than when only one SBAS was used. Additionally, the three integration methods had equal precisions. The experimental results are consistent with the error analysis results of Section 3.  Figure 9a,b depicts the horizontal and vertical errors for GPS, GLONASS, MSAS, and SDCM integration using CDI, MDI, and PDI that were discussed in Section 4 to combine multi-GNSS and multi-SBAS. Additionally, the standard deviations for each axis are presented in Table 2. By adding GLONASS and the SDCM corrections, the standard deviations of each axis decreased by 0.09-0.12 m more than in the case of GPS and MSAS. The experimental results also demonstrated that the cmlevel precisions of CDI, MDI, and PDI were measured equally. Therefore, the error analysis result   Figure 9a,b depicts the horizontal and vertical errors for GPS, GLONASS, MSAS, and SDCM integration using CDI, MDI, and PDI that were discussed in Section 4 to combine multi-GNSS and multi-SBAS. Additionally, the standard deviations for each axis are presented in Table 2. By adding GLONASS and the SDCM corrections, the standard deviations of each axis decreased by 0.09-0.12 m more than in the case of GPS and MSAS. The experimental results also demonstrated that the cm-level precisions of CDI, MDI, and PDI were measured equally. Therefore, the error analysis result that the three methods for multi-GNSS and multi-SBAS have same precision is reasonable.

Summary and Discussion
Besides positioning precision, issues such as computational complexity, availability, and flexibility are also important. In this section, a brief analysis of performance, including implementation issues, is given. Table 3 summarizes the results. As shown in Sections 3 and 4, the covariance analysis indicated that the three methods were the same in precision. However, they differed in computational complexity, availability, and flexibility. Table 3. Comparison of CDI, MDI, and PDI.

CDI MDI PDI
Precision

Flexibility Low Low High
To compare computational complexity, a number of operations, such as addition, multiplication, and division, to obtain a position were analyzed and the results are summarized in Table 4. The number of operations were determined by the number of GNSSs (m = N GNSS ), the number of SBASs (n = N SBAS ), and the number of visible satellites. For ease of comparison, the number of visible satellites of each GNSS was assumed to be the same (l = N 1 = N 2 = . . . = N N GNSS ). In CDI, the dimension of the H matrix and the number of inverse operations are independent of the number of SBASs, so just a simple operation was added to correct the weighting. By contrast, MDI complexity depends on the number of GNSSs, SBASs, and visible satellites. When multi-GNSS and multi-SBAS are concerned, MDI requires more operations. The PDI requires m times (4 × 4) matrix inversion, while CDI and MDI require one. However, the dimension of H matrix of MDI is nml and it is larger than that of PDI(l) and CDI(l). The numerical example of 3 GNSSs and 3 SBASs with 8 visible satellites shows that MDI requires more computational loads.
In GNSS, at least four pseudo-range measurements are required to find position and time. In PDI, each GNSS requires at least four visible satellites, while in CDI and MDI, the position can be found with four satellites regardless of GNSS type. This implies that the availability of PDI is worse than the other methods.
As PDI combines the navigation solutions, it can be easily implemented by combining two or more off-the-shelf receivers that support different GNSSs and SBASs. However, MDI and CDI require a modification of off-the-shelf receivers to integrate the corrections and pseudo-range measurements, such as system time difference correction. This means that PDI is more flexible than CDI and MDI in implementation and expansion.
In summary, CDI, MDI, and PDI have the same precision, but each method has different pros and cons of implementation. Although CDI has less computational load, it is less flexible. Therefore, CDI is suitable for a large number of GNSSs and SBASs. MDI is preferred in urban areas where lots of GNSS and SBAS signals are blocked because of high availability, even though it requires more computing resources. PDI is advantageous in implementation since it does not require modification of existing receivers, even though it has a lack of availability. As more GNSSs and more SBASs are accessible in near future, the ease of implementation might be more important in combing multi-SBAS corrections in multi-GNSS.

Conclusions
Currently, multi-SBAS utilization methods are not officially provided. Although a few methods have already been proposed, each method was verified in different environments and their pros and cons were not presented. In this study, to effectively utilize multi-SBAS corrections in multi-GNSS positioning, three methods, namely, CDI, MDI, and PDI, were presented and their performances were compared. The structures and navigation algorithms for single GNSS with multi-SBAS corrections were presented and then the results were expanded to multi-GNSS. The optimum weights were derived from covariance analysis and, when the optimum weightings of the three methods were applied to the WLS, the analysis revealed that the precisions of the three methods are equal. Experimental results using a commercial receiver and real measurements and corrections confirmed the analysis. Besides precision, implementation issues, such as computational complexity, availability, and flexibility, were also analyzed. CDI requires a small computational load, but is less flexible. MDI provides high availability with more computational complexity. PDI is more flexible, but has less availability. The results are expected to be a useful guideline in the design, implementation, and analysis of an efficient multi-GNSS positioning system using multi-SBAS corrections, which will be common in applications in the near future.