1. Introduction
The integrity of the Global Navigation Satellite System (GNSS) is one of the important factors to ensure civil aviation safety. There are three categories of GNSS integrity augmentation systems: Satellite-Based Augmentation System (SBAS), Ground-Based Augmentation System (GBAS) and Aircraft-Based Augmentation System (ABAS). The first two categories are at the system level, and the latter category is at the user level [
1]. ABAS can be implemented with Receiver Autonomous Integrity Monitoring (RAIM), which provides a navigation solution with guaranteed integrity by consistency checking among measurements [
2].
For an onboard application, RAIM is executed in two steps: the RAIM availability assessment and the satellite fault detection [
1]. The former (RAIM availability assessment) is used to assess in advance whether the navigation solution can meet the integrity risk requirement with the fault detection procedure. For decades, RAIM availability assessment has been achieved by calculating the protection level (PL), the upper bound of the position error corresponding to the integrity risk requirement [
3]. The threshold of the PL is the alert limit (AL), the upper bound of the user-allowed position error. If the PL is lower than the AL, RAIM is considered available; otherwise, it is considered unavailable.
Many studies have focused on determining how to obtain a lower PL to improve the availability of RAIM. Some of these studies were devoted to developing new navigation solution calculation methods, for example, the improved Integrity-Optimized RAIM (NIORAIM) [
4] and the optimal weighted average solution (OWAS) [
5] methods used for the snapshot RAIM algorithm. These methods can obviously decrease the PL with a slight increase in nominal position error. In addition, some studies have committed to accurately modeling the stochastic measurement noise, such as the discrete error-distribution (NavDEN) model proposed by Rife and Pervan [
6] and the distribution model considering both elevation angle and orbit type proposed by Fan [
7]. These measurement noise models are all helpful for obtaining a tight PL.
However, in the process of pursuing a lower PL, i.e., higher availability of RAIM, there is a key issue that is ignored by most researchers: whether the PL can accurately assess RAIM availability. Milner and Ochieng noted this issue [
3]. They qualitatively described that the classic PL, the product of the characteristic slope and the Minimum Detection Bias (MDB) proposed by Brown and Chin [
8], was too optimistic for RAIM availability assessment. The slope is a geometric feature-related parameter that qualitatively describes the relationship between positioning error and pseudo-range residual [
9]. The reason is “PL < AL” might not mean that the integrity risk satisfies the requirement for the measurement bias less than MDB. Meanwhile, its enhancement, abbreviated as the enhanced PL in the following, which provides an additional term to protect against the variation in position error proposed by Angus [
10], is too conservative. The reason is “PL ≥ AL” might not mean that the integrity risk exceeds the requirement. Here, measurement bias means the measurement error caused by the satellite fault, which is different from the measurement noise in the nominal mode. Furthermore, they proposed the ideal PL. It is the minimum PL value that guarantees the integrity risk, satisfying its requirement for arbitrary measurement bias. The ideal PL can prevent RAIM availability assessment from being optimistic or conservative. However, it cannot be solved analytically. A numerical search for the ideal PL begins with an improbably large value [
3], which leads to a large amount of calculation, increasing the computational burden of a GNSS receiver or an onboard computer.
In recent years, most researchers have focused on Advanced RAIM (ARAIM), in which PL is calculated after fault detection [
11]. ARAIM supports multi-constellation dual-frequency GNSS integrity monitoring. Multiple hypothesis solution separation (MHSS) algorithm is used in ARAIM. How to solve the accurate PL for MHSS is a research hotspot, including the PL calculation method for each fault mode and the optimization strategy for the ultimate PL [
12,
13,
14,
15,
16]. Jiang and Wang adopted the ideal PL [
3] in ARAIM and verified it was more accurate than other PLs for the ARAIM availability assessment [
17,
18]. ARAIM is still in the theoretical research stage and is not currently being applied in engineering practice.
Compared with ARAIM, RAIM has two deficiencies. The first is that RAIM is designed for a single constellation, monitoring only a single satellite fault [
2]. ARAIM is designed for double constellations, monitoring not only the single satellite fault but also the multiple satellite faults and the constellation fault [
12]. The second is classic, and the enhanced PLs are not rigorous enough for RAIM availability assessment, while the PL of ARAIM is much more rigorous. However, the on-board calculation of RAIM using the classic or the enhanced PLs is much less than that of ARAIM. For a single constellation, RAIM can still be used, but there is a problem needs to be considered, founding a RAIM availability assessment method both satisfying the rigor and maintaining the low on-board computational burden.
In this paper, a slope-based RAIM availability assessment method is proposed to solve the above problem. The characteristic slope is taken as the assessment basis. Using the ideal slope threshold, this method can achieve a consistent RAIM availability assessment with the ideal-PL-based method. The ideal slope threshold can be calculated offline and searched online because it is only related to one geometric parameter.
The remainder of this paper is organized as follows:
Section 2 states some technical backgrounds.
Section 3 reviews PL-based RAIM availability assessment methods, including the classic, enhanced and ideal PLs. The deficiencies of the classic and enhanced-PL-based methods can be analyzed quantitatively using the rates of optimistic or conservative assessment.
Section 4 proposes the slope-based RAIM availability assessment method after the derivation of the ideal slope threshold.
Section 5 gives an overview of the simulation results of the classic-PL-based, enhanced-PL-based, ideal-PL-based and slope-based methods.
Section 6 concludes this work with a brief summary. The discussion of this paper takes vertical integrity as an example, using one kind of classic RAIM snapshot algorithm, i.e., the least squares residuals (LSR) algorithm. In this paper, the measurement noise is assumed to be independent white Gaussian noise (WGN).
2. Technical Background
Before discussing the RAIM availability assessment methods, some technical backgrounds need to be stated, including the derivation of the integrity risk requirement for the single-satellite fault mode, the definition of the vertical characteristic slope and the specific meaning of RAIM being available.
2.1. Integrity Risk Requirement for the Single-Satellite Fault Mode
The integrity risk
[
19] is the probability of undetected faults causing unacceptably large errors in the estimated position [
20]. HMI is short for hazardous misleading information (HMI).
can be divided into three fault modes, the nominal mode, the single-satellite fault mode and the multiple-satellite fault mode, expressed as:
In Equation (1),
,
and
respectively represent the nominal, the single-satellite fault and the multiple-satellite fault modes, where
is the prior probability of the fault mode
and
is the probability of HMI under the
fault mode.
and
, where
means the number of combinations for choosing one element from K elements;
is the prior fault probability of a satellite, and
K is the total visible satellite number. Taking vertical plane for example,
is calculated with the following equation [
21]:
In Equation (2), is the vertical position error; is the vertical AL; and are the fault detection test statistic and threshold, respectively.
To ensure the integrity of the navigation system,
should be less than its requirement, denoted as
.
Given a geometry between the user and all-in-view satellites,
,
and
,
can be calculated following:
In Equation (4),
is the allowable false alarm probability under the nominal fault mode, satisfying
;
is the probability density function (PDF) of the normal-distributed VPE under the nominal fault mode with mean value 0 and standard deviation value
, where
is explained in
Appendix A.
According to GNSS Evolutionary Architecture Study (GEAS) report [
22], the integrity risk requirement allocated on the multiple satellites fault, denoted as
, can be set
[
22].
To ensure Equation (3) to be true, the integrity risk for single fault mode should satisfy:
where
. Thus, the problem of
has evolved into the problem of
.
Moreover, to ensure that
is nonnegative, the threshold for
, denoted as
, can be derived from the inequality
.
where
represents the inverse function of the cumulative distribution function (CDF) for the standard normal distribution.
indicates that only the integrity risk of the nominal and multiple-satellite fault modes have exceeded the total requirement, i.e.,
.
2.2. Vertical Characteristic Slope
The vertical characteristic slope is defined according to these two parameters [
2],
where
and
respectively characterizes
and
change caused by the measurement bias of the
m-th visible satellite, signed as
, is the faulty satellite. The details of
and
can be seen in
Appendix A and
Appendix B, respectively. For a specific bias value, a faulty satellite with a large slope value will present a high
, and a faulty satellite with a small slope value will present a low
. The “characteristic slope” will be abbreviated as “slope” hereafter.
2.3. Specific Meaning of RAIM Being Available
RAIM being available refers to , i.e., for the arbitrary measurement bias value while RAIM being unavailable refers to for at least one measurement bias.
A specific example is used to intuitively explain the meaning of RAIM being available. The 32-satellite GPS constellation is used in this example. The pseudorange measurement is assumed to be the dual-frequency ionosphere-free combination of L1 and L5. The standard deviation of the measurement noise for
, signed as
, is set according to the ARAIM interim report [
11]. For the location of 37° N latitude, 117° E longitude and height 0 m and the epoch of UTC 14 March 2019 17:15:00, there are 9 visible satellites with a masking angle of 10°. Their vertical slope values are recorded in
Table 1.
Figure 1 presents the base-10 logarithm of
for PRN5, PRN6 and PRN13 with the measurement bias
in the interval of
. These three satellites have the top three vertical slope values, as shown in
Table 1. The
curves for these three satellites follow the same order as the slope values, which illustrates that a faulty satellite with a large slope will present a high
. The
curves for the other six visible satellites must be lower than that of PRN 13 because their vertical slopes are smaller. In this example, RAIM would be unavailable if PRN5 was the faulty satellite because
is larger than
for
values in the range from 11 m to 15 m. Therefore, the intersection between the
curve of the faulty satellite and the
line means that RAIM is unavailable. RAIM would be available if the faulty satellite was one of the other visible satellites except for PRN5, because its
smaller than
at an arbitrary measurement bias. Therefore, the separation between the
curve of the faulty satellite and the
line means that RAIM is available.
Because the faulty satellite is unknown in actual situations, RAIM is considered available only if is true for the arbitrary measurement bias in the worst case, i.e., the satellite with the maximum slope being faulty.
3. PL-Based RAIM Availability Assessment
The classic PL, enhanced PL, and ideal PL are reviewed in this section.
3.1. Classic PL and Its Enhancement
The classic PL, denoted as
, is defined as the product of the maximum slope
and the minimum detectable bias
as follows:
where
is the noncentral parameter of the fault detection test statistic
under single-satellite fault mode [
23].
The enhanced PL, denoted as
, has an additional term that protects against variation in the random error of the position solution on the basis of
as follows [
10]:
where
and
is the standard deviation of the position error distribution.
Slope reflects the relationship between the position error and the pseudorange residual. PL is the projection of the pseudorange residual on the position error as shown in Equations (8) and (9). According to
derived in
Appendix B, the measurement bias
for a faulty satellite
, which causes
to obey
, is
. Based on the position error derivation in
Appendix A, this
makes the VPE under single-satellite fault mode obey:
Contrasting Equation (8) with Equation (10), the classic vertical PL, denoted as , is the expectation of at with the maximum slope. Therefore, at in the worst case. Similarly, the enhanced vertical PL, denoted as , is the upper quantile of the of distribution at with the maximum slope. Therefore, at in the worst case.
3.2. Deficiency of Classic- and Enhanced-PL-Based RAIM Availability Assessment
Here defines H0 as “,” and H1 as “,” in the worst case. Therefore “” or “” respectively mean the RAIM availability assessment is optimistic or conservative.
For
and
where
,
Similar to , changes with the value. Moreover, the value of is known at , for and for . A curve can be generated, which is higher than the curve with and lower than the curve with . The positional relationship among the curve, the curve, and the line jointly determine whether the PL-based RAIM availability assessment is conservative or optimistic. The following is a RAIM availability assessment discussion according to the positional relationship between the curve and the line.
3.2.1. Separation
In this situation, there are three kinds of positional relationships among the
curve, the
curve and the
line, as shown in
Figure 2.
The
curve intersects with the
line, and is higher than the blue dashed
curve, shown as the brown
curve in
Figure 2. In this situation, RAIM is unavailable and the VPL is larger than the VAL, i.e.,
, meaning a successful detection of “RAIM being unavailable”.
The
curve is separated from the
line, and is higher than the blue dashed curve
curve, shown as the black curve
in
Figure 2. In this situation, RAIM is available but the VPL is larger than the VAL, i.e.,
, meaning a conservative assessment.
The
curve is separated from the
line, and is lower than the blue dashed
curve, shown as the blue curve in
Figure 2. In this situation, RAIM is available, and the VPL is smaller than the VAL, i.e.,
, meaning a successful detection of “RAIM being available”.
Therefore, for the condition that the curve is separated from the line, the RAIM availability might be conservatively assessed.
3.2.2. Intersection
In this situation, there are also three kinds of positional relationships among the
curve, the
curve and the
line, as shown in
Figure 3.
The
curve intersects with the
line, and is higher than the blue dashed
curve, shown as the brown
curve in
Figure 3 In this situation, RAIM is unavailable and the VPL is larger than the VAL, i.e.,
, meaning a successful detection of “RAIM being unavailable”.
The
curve intersects with the
line, and is lower than the blue dashed
curve, shown as the black
curve in
Figure 3. In this situation, RAIM is unavailable, but the VPL is smaller than the VAL, i.e.,
, meaning an optimistic assessment.
The
curve is separated from the
line, and is lower than the blue dashed
curve, shown as the blue curve
in
Figure 3. In this situation, RAIM is available, and the VPL is smaller than the VAL, i.e.,
, meaning the successful detection of “RAIM being available”.
Therefore, for the condition that the curve intersects with the line, the RAIM availability might be optimistically assessed.
3.2.3. Tangency
In this situation, there are two kinds of positional relationships among the
curve, the
curve and the
line, as shown in
Figure 4.
The
curve intersects with the
line, and is higher than the blue dashed
curve, shown as the brown
curve in
Figure 4. In this situation, RAIM is unavailable, and the VPL is larger than the VAL, i.e.,
, meaning a successful detection of “RAIM being unavailable”.
The
curve is separated from the
line, and is lower than the blue dashed
curve, shown as the black
curve in
Figure 4. In this situation, RAIM is available and the VPL is smaller than the VAL, i.e.,
, meaning a successful detection of “RAIM being available”.
Therefore, for the condition that the curve is tangent to the line, both optimistic and conservative assessments can be prevented.
Based on the above analysis, the accuracy of PL-based RAIM availability assessment depends on the positional relationship between the curve and the line. There is a risk of conservative assessment when the curve is separated from the line and risk of optimistic assessment when the curve intersects with the line. PL-based RAIM availability assessment is accurate only if the curve is tangent to the line. Because only the value at a specific , i.e., , is determined for or , the position relationship between the entire curve and the line is uncertain. Consequently, both optimistic and conservative assessments might happen when using or to assess whether vertical RAIM is available.
Because , the curve for is much higher than that for . Thus the possibility of intersection between the curve and the line for is much higher than that for , which may lead to an optimistic assessment, while the possibility of separation between the curve and the line for is much higher than that for , which may lead to conservative assessment. Consequently, the optimistic assessment risk of using is higher than that of using for vertical RAIM availability; in contrast, the conservative assessment risk of using is higher than that of using .
3.3. Ideal Protection Level
According to the above analysis, the ideal positional relationship between the curve and the line is tangency, which can prevent both optimistic and conservative RAIM availability assessments. The ideal VPL, denoted as , proposed by Milner and Ochieng, satisfies this condition. It matches the exact required integrity risk for the worst-case bias (WCB), the measurement bias presenting the highest integrity risk. Thus forms a curve tangent to the line. If , the curve must be separated from the line, which means that vertical RAIM is available; otherwise the curve must be tangent to or intersect with the line, which means that vertical RAIM is unavailable.
is the solution of:
4. Slope-Based RAIM Availability Assessment
In addition to the ideal PL, there is another ideal test statistic for RAIM availability assessment: the slope. Both optimistic and conservative RAIM availability assessments can be prevented using the slope once an ideal threshold is found. The following is a deviation of this ideal threshold.
4.1. Derivation of the Ideal Threshold for the Slope
The ideal slope threshold derivation begins with searching for a condition satisfying
for all possible measurement bias values of an arbitrary faulty satellite. To ensure
constantly true, the maximum value of
should be less than
,
Where:
Substituting Equation (7) into Equation (15),
Setting
, Equation (16) can be transformed into:
Setting
,
can be taken as a function of
as follows:
When substituting Equation (18) into Equation (14), an ideal threshold must exist for
with the given
,
,
,
and
, denoted as
, which satisfies the limit situation
According to Equations (18) and (19), once , , and are given, is determined by only one parameter .
Figure 5 presents the numerically solved
for
values with a step of 0.001 at
,
,
,
,
and
. As shown in
Figure 5,
decreases with the increase of
. Each
curve exhibits a nearly constant segment at the beginning and a sharply decreasing segment at the end. The
curve ends when
reaches
, and the spacing between two adjacent curves obviously decreases as
increases.
4.2. Practical Meaning of the Ideal Slope Threshold
Here, the specific example in
Section 2.3 is used to intuitively explain the practical meaning of
. The
value for this example is 4.944. Given
,
,
,
and
, the ideal slope threshold can be numerically solved according to Equation (19), and
.
Figure 6 presents the
curves for
(PRN5),
(PRN6) and
.
for
is calculated according to Equation (19). The
curve for
is tangent to the
line. The
curve for
, which is larger than
, is intersected with the
line. In contrast, the
curve for
, which is smaller than
, is separated from the
line.
It can be deduced that the slope value determines the positional relationship between the curve and the line. forms a tangent curve to the line. If the slope value of a faulty satellite is larger than , its curve would be intersected with the line, meaning that can be satisfied for at least one possible measurement bias value. If the slope value of a faulty satellite was smaller than , its curve would be separated from the line, meaning that can be ensured at an arbitrary measurement bias value. From the analysis of this specific example, the slope is an ideal test statistic for RAIM availability assessment, with the ideal threshold calculated according to Equation (19). Therefore, for each observation epoch, the slope of a visible satellite can be considered as “large slope” if it was larger than and considered as “small slope” if it was smaller than .
4.3. Comparison of the Ideal Slope Threshold and the Ideal Protection Level
The functions of
and
are identical, forming a tangent
curve to the
line. Referring to the equation for
, i.e., Equation (18), the equation for
, can be formulated as follows:
Analyzing this equation, is determined by and with the given , and . Therefore, is only related to while is related to both and .
Considering that the PDF of the noncentral distribution is too complicated, both and should be solved numerically, which will sharply increase the computation burden of a GNSS receiver or an onboard computer. Thanks to the one-to-one correspondence between and , can be calculated offline for discrete values in the range from 1 to and saved in a receiver. It can be on-board searched from the presaved data according to the specific value. However, needs to be calculated online after both and obtained.
4.4. Slope-Based RAIM Availability Assessment Method
Slope-based RAIM availability assessment should be implemented in the worst case to fully prevent the integrity risk. Using slope-based method, RAIM is considered available if the maximum slope is less than the ideal threshold
; otherwise, it is considered unavailable. The specific execution process for slope-based RAIM availability assessment is presented in
Figure 7.
As shown in
Figure 7, the inputs are the observation matrix
H and the weighted matrix
P. The first step is an assessment based on the total number of visible satellites
K. RAIM is considered unavailable if
because fault detection cannot be executed with less than 5 visible satellites. For
, the second step is an assessment based on the
value. RAIM is considered to be unavailable if
because only the sum of
and
has exceed
when
exceeds
. For
, the last step is an assessment based on the maximum slope
. RAIM is finally considered available if
.
can be calculated online according to Equation (6) with given
,
,
,
and
.
should be calculated offline with discrete
with a small step size for different numbers of visible satellites and saved in the GNSS receiver or the onboard computer. Assuming that the total number of visible satellites is
K and
is in the interval of
, where
and
are the indexes of two adjacent discrete points
and
presaved for
K visible satellites,
corresponding to
should be assigned as:
This value is calculated according to a linear fit for
in the interval of
. Because the
curve is convex as shown in
Figure 5, the assigned
value is smaller than the real
value, which may cause a small conservative assessment risk. However, if the
step size is small enough, the assigned
value would be nearly equal to the true value, thereby preventing the small conservative assessment risk.
A simulation is designed to find the desirable range
step.
times of
and
for different location and epochs are collected. As shown in
Table 2, the times of conservative assessment increases with step size widen. Taking into account both the amount of calculation and conservative assessment rate, 0.01~0.02 is the desirable range
step because it is the maximum step size with 0 time of conservative assessment.
It should be mentioned that the calculation amount of the ideal-PL and
is exactly the same. The ideal PL procedure begins with an improbably large VPL of 2000 m and halves the search step by checking if the corresponding integrity risk exceeds the requirement [
3]. For each step of the iteration, the integrity risk needs to be calculated for different bias values with a fixed step. This process is computationally intensive. Similarly, the
procedure begins with an improbably large slope value of 15 and halves the search step. Compared with the ideal-PL-based method, the slope-based method separates the process of numerical iteration from on-board RAIM availability assessment, reducing the burden of on-board computing.
5. Simulation
To compare the performance of the PL-based and slope-based methods, the vertical RAIM availability assessment for a 32-satellite GPS constellation is simulated in worldwide (latitude 60° S~60° N and longitude 180° W~180° E) for a whole day (13 March 2019 0:00:00~24:00:00). The simulation area is meshed as the grid of
and the simulation time step is 300 s. The masking angle is set to 10°. The dual-frequency ionosphere-free combination of L1 and L5 is assumed to be the pseudorange measurement. The standard derivation of the measurement noise is set according to the ARAIM interim report [
11].
The vertical RAIM availability assessment is executed for all grid points, i.e.,
grid points, using the classic-PL-based (
), the enhanced-PL-based (
), the ideal-PL-based (
) and the slope-based methods respectively at each simulation epoch. In the simulation, the parameters related to vertical RAIM availability are
,
[
21],
[
22],
[
11],
, and
. The
values are pre-calculated before the simulation with an
step size of 0.01.
5.1. Specific Example Analysis
Three specific examples for a single grid point and a simulation epoch are chosen to show the RAIM availability assessment using the four methods in detail.
Table 3 records the VPLs,
and
values for these specific examples.
Figure 8,
Figure 9,
Figure 10 and
Figure 11 present the real
curve for
, the
curve for
, and the
curves for
,
and
in these examples, respectively. In
Figure 8,
Figure 9,
Figure 10 and
Figure 11, the
curve for
and the
curve for
are tangent to the
line, while the
curve for
is separated from the
line. The positional relationship between the
curve for
and the
line could be either intersection or separation.
As shown in
Figure 8, the real
curve for
intersects with the
line, which means that RAIM is actually unavailable in the first example. The
curve for
is higher than the real
curve, i.e.,
, which means that RAIM is assessed to be available when using
. All of the
curve for
, the
curves for
and
are lower than the real
curve, i.e.,
,
and
, which means that RAIM is assessed to be unavailable when using
,
or
. Therefore, the RAIM availability is optimistically assessed when using
and successfully assessed using when using
,
or
.
As shown in
Figure 9, the real
curve for
is separated from the
line, which means that RAIM is actually available in the second example. The
curve for
, and the
curves for
and
are higher than the real
curve, i.e.,
,
and
, which means that RAIM is assessed to be available when using
,
or
. The
curve for
is lower than the real
curve, i.e.,
, which means that RAIM is assessed to be unavailable when using
. Therefore, RAIM availability is conservatively assessed when using
and successfully assessed when using
,
or
.
As shown in
Figure 10, the real
curve for
is separated from the
line, which means that RAIM is actually available in the third example. Both the
curve for
and the
curve for
are higher than the real
curve, i.e.,
and
, which means that RAIM is assessed to be available when using
or
. Additionally, both the
curves for
and
are lower than the real
curve, i.e.,
and
, which means that RAIM is assessed to be unavailable when using
or
. Therefore, RAIM availability is conservatively assessed when using
or
, and it is successfully assessed when using
or
in this example.
These three specific examples intuitively illustrate that both optimistic and conservative RAIM availability assessments might happen when using the classic-PL-based method; only conservative assessment might happen when using the enhanced-PL method; and both optimistic and conservative assessments can be prevented when using the ideal-PL-based or the slope-base methods.
5.2. Simulation Results Statistical Analysis
Both the rates of optimistic and conservative assessments for RAIM availability are calculated for each grid point using the classic-PL-based, enhanced-PL-based, ideal-PL-based and slope-based methods. The optimistic assessment rate is indicated by the ratio between the count of epochs in which RAIM is assessed to be available but is actually unavailable and the count of epochs in which RAIM is actually unavailable. The conservative assessment rate is indicated by the ratio between the count of epochs in which RAIM is assessed to be unavailable but is actually available and the count of epochs in which RAIM is actually available.
Table 4 records the specific optimistic and conservative assessment data for all 4 kinds of RAIM availability assessment methods. As shown in
Table 4, both the rates of optimistic and conservative assessments are 0 for each grid point using the slope-based and ideal-PL-based methods. These prove that both optimistic and conservative assessments can be prevented. The reason for the 0 conservative assessment rate of the slope-based method is that the
step size is small enough (0.01) for offline calculation of
, leading to the assigned
value being nearly equal to its true value.
For the classic-PL-based method, there are both grid points with nonzero optimistic assessment rate and grid points with nonzero conservative assessment rate.
Figure 11 and
Figure 12 present the optimistic and conservative assessment rates for each grid point using the classic-PL-based method, respectively. As shown in
Table 4 and
Figure 11, there are 2032 grid points with nonzero optimistic assessment rate, 4.65% of the total grid points, and there are 41 grid points with a 100% optimistic assessment rate, meaning that optimistic assessment always happens when RAIM is unavailable at these grid points. As shown in
Table 4 and
Figure 12, there are 16,678 grid points with nonzero conservative assessment rate, 38.18% of the total grid points, and the maximum conservative assessment rate is 2.59%. Comparing
Figure 12 with
Figure 11, the coverage area of conservative assessment is much larger than that of optimistic assessment, but the maximum optimistic assessment rate is much higher than the maximum conservative assessment rate for a single grid point using the classic-PL-based method.
For the enhanced-PL-based method, there are no grid points with a nonzero optimistic assessment rate. This finding illustrates that the value of
is small enough at
for the enhanced PL, which leads to the
curve always being separated from the
line, preventing optimistic assessment. However, there are 43,038 grid points with nonzero conservative assessment rate, as shown in
Table 4 and
Figure 13, representing 98.53% of the total grid points. The maximum conservative assessment rate is 8.15%. Comparing
Figure 13 with
Figure 12, the conservative assessment risk using the enhanced PL is much higher than that using the classic PL, which is represented by the much larger conservative assessment coverage area and the higher conservative assessment rate for a single grid point.
According to the above simulation results, the performance of the classic PL is the worst on RAIM availability assessment for both the risks of optimistic and conservative assessments. In particular, an optimistic assessment might cause HMI, which is intolerable. The performance of the enhanced PL is better than that of classic PL because the optimistic assessment is prevented. However, the risk of conservative assessment is significant, reducing RAIM continuity. The performances of the slope and the ideal PL are optimal, preventing both optimistic assessment and conservative assessment.