Reliability Analysis of Network Real-Time Kinematic

Themultistate reliability theory was applied to the network real-time kinematic (NRTK) data processing chain, where the qualities of the network corrections, baseline residuals, and the associated variance-covariance matrices are considered as the system state vectors. The state vectors have direct influence on the rover receiver position accuracy. The penalized honored stochastic averaged standard deviation (PHSASD) is used to map the NRTK sensitive data, represented by the states vectors to different levels of performance. The study shows that the improvement is possible by identification of critical components in the NRTK system and implementation of some parallelism that makes the system more robust.


Introduction
High accuracy positioning with GNSS is carried out using both code and carrier phase data from the GNSS satellites.To obtain position accuracies at the cm or mm level using carrier phase data, an important part of the data processing is to estimate the initial oscillator phase offset, the socalled ambiguity, for each receiver-satellite pair.Resolution of ambiguities requires that the influence of most errors sources in the positioning process is reduced to the cm-level, and high accuracy positioning is therefore often done in a relative mode where the position of a GNSS receiver located in an unknown position (the rover) is determined relative to one or more reference stations located in positions known on beforehand [1].
With relative carrier phase based GNSS positioning the effects of the distance-dependent error sources such as uncertainties in satellite positions and atmospheric effects on the satellite signals induced by the ionosphere and troposphere are reduced.Also the effects of satellite and receiver clock errors in the positioning process are reduced by relative positioning, and all this in combination makes it possible to resolve the ambiguities and thereafter obtain positions for the rover at the cm or mm level.
For high accuracy GNSS positioning in real time, the real time kinematic (RTK) technique has been developed.Traditionally this is based on a reference station transmitting data to the rover where the data is used in estimation of the position of the rover in a relative or differential mode [2].
Using a network of reference stations for RTK, the socalled NRTK technique, provides the opportunity for applying more advanced algorithms for estimation of the distancedependent errors within the network and thereby possibilities for providing a more robust service.Such operational NRTK services exist in many regions and countries today and have become an indispensable tool in high accuracy navigation and surveying.
A brief description of the NRTK functionality is as follows: the first step is collection of raw observations from the network of reference stations, solving the ambiguities within the reference network, and generating error estimates.Then an interpolation/smoothing scheme is applied to generate the NRTK corrections for the user location.For information on how to avoid a loss of information under interpolation of NRTK data, the interested reader is referred to [3].
The NRTK corrections are then transmitted to rover receivers.Several NRTK techniques exist and the most commonly used at present are, for instance, the Master Auxiliary Concept (MAC) [4,5], the Virtual Reference Station (VRS) concept [6], and the FKP techniques [7], as well as the Network Adjustment (NetAdjust) concept developed by Raquet and Lachapelle [8,9].Multistate system reliability theory has been a research topic for many years, for instance, extension of the system from two-state to multistate reliability [10,11] to compute the mean performance level at any given time  and stochastic evaluation and bound computation of multistate coherent systems [12], further, studies on application of reliability analysis to GNSS data processing [2], a comparative GNSS reliability analysis [13], reliability analysis under GNSS weak signals [14], accuracy and reliability of multi-GNSS real-time precise positioning [15], and robust reliability testing in case of signal degradation environment [16,17].
In order to avoid confusion between the reliability theory definition used in different fields, we refer to the traditional reliability theory to the context of statistical testing based on the theory developed by Baarda [18], while component, system, binary, and multistate reliability terms will be used in case of component and system reliability computation.
Our aim is to provide the user in the field with continuously high quality corrections with the ability to identify the periods for which the reliability of the network RTK performance is reduced in terms of accuracy and availability.Therefore, solution quality indicators describing the reliability of the network RTK are needed to transfer the status of the network to the user in the field.Intensive research has been conducted recently in this field to derive these quality indicators and can be classified into two main classes; (i) spatially correlated (ionosphere, troposphere, and orbital) error indicators; (ii) residuals errors indicators.Most network RTK used quality indicators are residual integrity monitoring (RIM) and irregularity parameters (IP) quality indices [19], residual interpolation uncertainty (RIU) [20], geometrybased quality indicator (GBI) [21], and ionospheric index I95 [20].An elegant presentation summarizing the network RTK quality indicators can be found in [22].
In order to apply reliability analysis to a NRTK system, the starting point is the decomposition of the block diagrams of NRTK processing chain into simple components and computation of the system reliability.Figure 1 shows three levels of data processing modules, the network, the baseline, and the rover receiver modules, where  denotes rover,  denotes reference station, and (, , ) denote errors in the position , ,  coordinates, respectively.Based on these levels, we can build the reliability block diagrams for the NRTK processing chain and compute the reliability for the entire system.
The rest of the paper is organized as follows: Section 2: brief introduction to component, system, and binary state reliability theory and deterministic and stochastic reliability.
Section 3: NRTK blocks diagram determination and module reliability computations.
Section 4: multistate reliability theory applied to NRTK processing chain.
Section 5: some test results.
Section 6: procedures used to validate the NRTK system reliability.
Section 7: discussions and conclusion.
Test data used in this investigation is described in Appendix A.

Reliability Analysis
The aim of this section is to introduce the basic of the component, system, and binary state reliability theory, including deterministic and stochastic reliability, followed by decomposing the block diagrams of the NRTK processing chain into simple components, and compute the entire system reliability.The description of the structural relationship between the components and the system must be defined.Figure 2 illustrates the concept and each main block will be treated separately in the coming sections.

Structure Functions.
In order to construct the entire reliability block diagram of the NRTK data integrity, it is necessary to define the elementary building blocks as Bernoulli indicator function [23, pp. 28-29].The function will signal if a unit or system is functioning or not.
A random variable  is said to be a Bernoulli random variable if its probability function is given by (1).The indicator and structure functions are given below by following the notation given by Natvig [24]: 1, success with probability , 0, failure with probability (1 − ) . (1) where ∐ is read "ip" and denotes the parallel coupling operator.
Journal of Electrical and Computer Engineering Note that the structure functions given by ( 4) and ( 5) are an n-out-of-n structure and an 1-out-of-n structure, respectively.Other structure functions exist, for instance, the bridge Natvig [24, p. 12].

System Reliability Computation.
The structure functions are defined; now is time to compute the system reliability.
We move from the deterministic model to the stochastic one by introducing the random variables.Some notations are needed to represent the state vector, structure function, and the reliability.We follow the notations given by Hoyland and Rausand [25,.We denote the state variables of the  independent component at time  by  1 () ,  2 () , . . .,   () .
Let  denote the system reliability, then we have To avoid confusion,   = (  () = 1) is the probability of functioning of the th component and referred to as the component reliability.

Reliability of Series Structures.
The reliability function   (p) of the series system of  independent components is given by the expression: If all components have the same (), (12) becomes {()}  .For  = 5 and  = 0.99, then the reliability   (p) = 0.951.An important remark is that the reliability of a series structure is at most as reliable as the least reliable component, that is,   (p) ≤ min  (  ()).

Reliability of Parallel Structures.
The reliability function   (p) of the parallel system of  independent components is given by the expression: If all components have the same (), (13) becomes 1 − {1 − ()}  .For  = 5 and  = 0.99, then the reliability   (p) = 1.
An important remark is that the reliability of a parallel structure is at least as reliable as the most reliable component, that is, (p) ≥ max  (  ()).
Details on how to compute the reliability of parallel structure in general are given in Appendix C.

Reliability of k-out-of-n Structures.
The reliability function  (−) (p) of the k-out-of-n system of  independent components with equal probability   () =  is given by the expression:

Reliability in the NRTK Processing Chain
The aim of this section is to determine the structure functions, NRTK module's reliability, and the corresponding block diagrams.

Considerations around NRTK Data
Processing.Some considerations around the NRTK data integrity are introduced below as a background for the design process and to ease the discussions in the following sections.For more information on GNSS data processing, the reader is referred to [1,2,26,27].The key to precise positioning is the correct ambiguity resolution and validation.With ambiguities resolved to wrong integer numbers, there will be offsets in the position solution, and with float ambiguities (ambiguities that are not fixed to integer values) the position solution is inaccurate and also very unstable and sensitive to changes in satellite geometry.Good satellite-receiver geometry, as, for instance, expressed by the so-called DOP factor (dilution of precision), is important to perform successful ambiguity resolution and achieve centimeter level accuracy in real time.
Spatiotemporal models that describe well the variations of the spatially correlated errors in the corrections field are also an important key for reliable NRTK positioning.
Robust estimation algorithms to handle large data sets are also a key factor becoming more important in the future as observations from several GNSS systems to a larger degree will be combined in one processing loop.Today, most NRTK systems operate with data from the American GPS and the Russian GLONASS system.Including data from the European Galileo as well as the Chinese Beidou systems in NRTK operations will soon be the norm for most NRTK services.With satellites from more GNSS systems being available the satellite-receiver geometry on the rover side is improved.This is especially important when the user is operating in constricted environments such as narrow street canyons or forest areas.

NRTK Corrections Reliability Analysis.
The main function of the NRTK is to provide the rover in the field with high quality corrections on an epoch-by-epoch basis.
From Figure 2, five modules  , are defined for the network and the reliability of each module will be evaluated. , corresponds to level  ∈ {, , } and module  ∈ {1, 2, . . ., 5}, where  stands for network,  for baseline, and  for rover receiver.For instance, the module  ,1 corresponds to the first network module, the Station Data Integrity, as shown in Figure 2, and  ,1 denotes the first of the rover modules, i.e., the Prediction Function as shown in Figure 2.

Reference Receivers Data Integrity.
Generation of high quality raw observations at reference receivers requires reliable hardware (HWC) and software (SWC) components.Expensive hardware and sophisticated algorithms are keys to achieve this goal.Figure 3 shows the concept and each component will be treated separately in the next sections.
Let software and hardware components be represented by the modules  ,1, and  ,1, , respectively.Software Component of Module  ,1 Definition.SWC requires an ensemble of sequential checks on raw observations.This includes (i) Let  1,1 denote the satellite data integrity algorithm.This algorithm will discard the measurements from unhealthy satellite(s) or from satellite(s) for which we do not have the orbital data.
(ii) Let  1,2 denote controlled cycle-slip algorithm.This task requires investigation of carrier phase discontinuities by examination of loss of lock (LLI) indicator and signal-to-noise ratio (SNR) flags.
Figure 5 shows the block diagram of the network reweighting algorithm component  1,10 .

Block Diagram and Reliability of Software Component.
The structure function of the software component is well described by the 5-out-of-10 structure function.This means that the more algorithms check is, the more reliable raw observations become.
In order to produce reliable raw observations of high quality, it is necessary to perform at least five checks from a total of ten.With ten algorithms' check, we can generate high quality raw observations while five algorithms' check x 1,7 x 1,8 x 1,9 x 1,10 x 1,2 x 1,2 x 1,3 x 1,3 x 1,4 x 1,4 x 1,5 x 1,6 x N,10,1 x N,10,2 x N,10,3 will produce acceptable level raw observations.The selection of the algorithms is independent of the order.The structure function Φ(x) of Figure 4 is given by the expression Assuming that the individual algorithms are independent with equal probability (), then the reliability of software component is Due to the fact that y() ∼ binom (, ()).The block diagram of the  ,1,10 is given by Figure 5.
The corresponding reliability reads  The elements of the hardware component of the network module  ,1 are as follows: (i) Let  1,1 denote the GNSS receiver type.
(ii) Let  1,2 denote the rover software known as firmware to decode the GNSS signals.
(iv) Let  1,4 denote the choke ring that allows better reception of low elevation angle GPS satellites and improved multipath rejection.
(v) Let  1,5 denote duplicated system (as discussed in the next section).
The structure function Φ(y) is given by the expression The reliability of the hardware of the module  ,1 reads Hardware Component Improvement.In order to ensure continuous raw data delivery at the station, duplicated hardware components are recommended.This task is accomplished by a parallel coupling of the module  ,1 , defined by Figure 6.
Applying the definition of reliability computation of parallel coupling (Equation (C.2)), the duplicated HWC reliability reads where   () is given by the Equation (19).The drawback of a duplicated system is the financial issues.

Network Ambiguity Resolution.
As mentioned in Section 3.1.The key for precise positioning is the correct ambiguity resolution and validation.The module  ,2 is composed of two main parts, the first part is the ambiguity processing algorithms performance, and the second part is the statistical ambiguity quality indicators derived from the processing chain of the ambiguity.The components of the module  ,2 are as follows: (i) Let  2,1 denote the float solution of the ambiguity obtained via least square or Kalman filter.
(iii) Let  2,3 denote the validation procedures to validate the final solution.
(iv) Let  2,4 denote the administration of the ambiguities.
(v) Let  2,5 denote the success rate of the ambiguities resolution.
(vi) Let  2,6 denote the ambiguity dilution of precision (ADOP) [30].ADOP measures the precision of the ambiguities and can be viewed as a quality indicator.
The block diagram of the  ,2 is given by Figure 7. Structure Function of  ,2 .The structure function of the module  ,2 reads Reliability of  ,2 .Reliability of the module  ,2 reads dropping the index 2 in the expression of   () and compute

Network Corrections
Quality.The quality of the network corrections depends on various parameters, for instance, the estimation algorithms, network status (sparse/dense), the covariance functions, and the smoothing/interpolation algorithms.The components of the module  ,3 are as follows: (i) Let  3,1 denote the network reference receivers' separation.Dense network is attractive because the network corrections are better estimated with short distances between reference receivers.
(ii) Let  3,2 denote the quality of the estimation algorithm used to estimate the parameter vector  ∈ Θ.
(iii) Let  3,3 denote the quality of the covariance function used to model the network correlation errors.
(iv) Let  3,4 denote the quality of the interpolation algorithm used to generate the user corrections.Parallel interpolation algorithms will enhance the quality of the user corrections generation and avoid the information loss.
Block Diagram of  ,3 .The block diagram of the module  ,3 is given by Figure 8.
Structure Function of  ,3 .The structure function of the module  ,3 reads Reliability of  ,3 .Dropping the index 3, the reliability of  ,3 reads Amelioration Potential of  ,3 .Our aim is to provide the user in the field with high quality corrections on an epochby-epoch basis.The interpolation/smoothing algorithm plays a central role.We can implement different parallel interpolation/smoothing algorithms that compete about the quality of service parameters.The corrections will be sent from the algorithm with higher score.For more information on this topic, the interested reader is referred to [3].

Network Data Integrity.
The module  ,4 responsibility is to carry out the quality control on the corrections field and the corresponding variance-covariance matrices.This includes the following: (i) Let  4,1 denote global test statistics to detect any extremal events that can bias the rover position, corrections field investigation.
(iii) Let  4,3 denote the application of the imputation algorithm to compute the statistics.
Serial coupling is the appropriate choice for the module  ,4 and the block diagram is given by Figure 10.
3.2.5.Network Probability Score.The module  ,5 computes the network quality indicators in terms of the successfully ambiguities resolution and the quality of the network corrections.This is the first state vector of the system under investigation.
(i) Let  5,1 denote quality indicator for the corrections field.
(ii) Let  5,2 denote quality indicator for the uncertainty of corrections field.
(iii) Let  5,3 denote quality indicator for the ambiguities expressed by ADOP from Section 3.2.2.
(iv) Let  5,4 denote the number of common satellites used in the computation.
(v) Let  5,5 denote the number of rejected satellites from computation.
Note that  5,1 and  5,2 are actually the network RTK quality indicators (RIM, IP, RIU, GBI, and ionospheric status indicator), defined in Section 1.
Serial coupling is the appropriate choice for the module  ,5 and the block diagram is represented by Figure 11.

NRTK Baseline Reliability
Analysis.This is the second level of the NRTK data processing.The corrections generated by the network are involved to generate the computation point (CP), then the unknown rover coordinates are determined relative to the computation point.This method is known as the relative positioning technique.Similarity between the baseline and the network data processing exists as we see in the coming subsections.

Rover Receiver Data Integrity.
The module  ,1 is similar to the module  ,1 defined in Section 3.2.1.The quality of the raw observations collected by the rover receiver depends strongly on the statistics methods used to check for anomalies.All variables defined in the module  ,1 are applicable to the module  ,1 .

Baseline Ambiguity Resolution.
The module  ,2 is similar to the module  ,2 defined in Section 3.2.2.All variables defined for the module  ,2 are applicable to the module  ,2 .

Baseline Data Integrity.
The module  ,3 is similar to the module  ,4 defined in Section 3.2.4.All variables defined for the module  ,4 are applicable to the module  ,3 .

Baseline Probability
Score.The module  ,4 is similar to the module  ,5 defined in Section 3.2.5.All variables defined for the module  ,5 are applicable to the module  ,4 .

NRTK Rover Reliability Analysis.
The investigation of the rover position error is the final check.The quality is measured in terms of standard deviations of the topocentric coordinates (, , ).

Rover Prediction Function.
This module  ,1 uses information from the double-difference variance-covariance matrix to compute the prediction of the position error.The number of satellites used in the computation is considered as a parameter.
(i) Let  3,1 denote inspection of the main diagonal of variance-covariance matrix.

Rover Position
Error.This module  ,2 computes the standard deviations of the rover position error (, , )  and assigns a final score to determine the state of the rover receiver accuracy.

Multistate Reliability Analysis
Since the rover position accuracy cannot be represented by a binary system with two performance level states as functioning or failed, the multistate system (MSS) approach is chosen to deal with situations where more than two levels of performance are considered.The material used to construct this section is from Natvig [32].

Definition of NRTK Performance Levels.
Based on the values computed from the score modules  ,5 (Section 3.2.5)and  ,4 (Section 3.3.4),respectively, a single judging number is assigned to determine the performance level of the rover position accuracy.
The states represent level of performance ranging from the perfect functioning level perfect down to the complete failure level catastrophic.Five states are defined for a NRTK system, namely, start, perfect, acceptable, rejected, and catastrophic states.
The rover position accuracy is well described by the state diagrams given by Figures 12 and 13, respectively.
The user requirements regarding the position accuracy are satisfied.
(3) State 2: rejectable position accuracy.Major complications due to the atmosphere, multipath, or algorithms failure.The user requirements are not satisfied.
(4) State 1: catastrophic state.The NRTK system is down and not delivering the corrections to the user in the field.
(5) State 0: start state.The process always starts at this state and can reach any other states.
Note that the states 0, 4, 3, 2 are transient states while the state 1 is absorbing state.Note that the probabilities  | must be computed from real data.

Penalized Honored Stochastic Averaged Variance.
Based on state vectors data, our aim is to construct the NRTK state diagram from the network and baseline.The variancecovariance matrix (VCM) is considered as the state vector and the average variance with respect to the number of observed satellites ( sat ) is used to compute the quality indicator on an epoch-by-epoch basis.In addition, the number of rejected satellites ( rej ) and the geometry factor (DOP) are used to penalize/honor the average variance.

Penalized Average Variance Component.
The total average variance computed from the VCM shall be penalized in case of rejection of satellite(s) with bad data and causes the increase of DOP indicator.The penalized function shall look similar to the following: (1) Penalized least square (PLS) proposed by Green and Silverman [33, p. 5] where   is the observations, (  ) is the curve we fit to the data,   (⋅) is the second derivative of the function (⋅), and  is the smoothing parameter and defines the rate of change between the residuals and local variations.Anyway, minimizing () gives the best compromise between smoothness and goodness-offit.A large value of  will make the penalty term more in action, while with a small value the first term will be the main contribution.
(2) Information criteria type penalizing the model complexity.Denote by  the model to be investigated and dim() is the length of its parameter vector .
Akaike's information criterion (AIC) [34, Chap.2] The Bayesian information criterion (BIC) of Schwarz (1978) takes the form of a penalized log-likelihood function where the penalty is equal to the logarithm of the sample size times the number of estimated parameters in the model [34,Chap. 3] 4.2.2.Honored Average Variance Component.Detection and rejection of satellite(s) with bad data are a good thing.The check algorithms shall be honored as long as the DOP values remain in the acceptable region.The value of the horizontal dilution of precision (HDOP) is expected to be less or equal 2.0, that is, HDOP ≤ 2.0.
Forming a new stochastic variable   = HDOP  −  HDOP , then we can monitor the values of   over time.Note that  HDOP corresponds to the mean value of the HDOP in time span Δ.

Balanced Average Variance.
Our aim is to put together the pieces defined in Sections 4.2.1 and 4.2.2, respectively, and find a way to balance between the satellite(s) rejection  rj and the HDOP value.
The exponential reweighting algorithm type is an option.The algorithm places more importance to more recent data by discounting older data in an exponential manner.For epoch , let   = (( obs, −  rej, ) −  const )/ sys , where  const ≥ 8 is the user defined parameter and corresponds to the minimum number of satellites required to compute a reliable solution and preserve a good HDOP value, and  sys = 31 is the total satellites in GPS constellation.
A suitable stabilization factor  ∈ [0, 1] is chosen such that the penalized honored average variance ( phav ) reads where  , are diagonal elements of the covariance matrix of the baseline,  obs, is the number of satellites with valid data used in the computation,  rej, is the number of rejected satellites by the algorithms, and   is the window size which is user defined.

Penalized Honored Average Variance Validation.
The parameter vector of (35) is  = (, , ).Our aim is to study the variation of the second term of (35) and try to get some valid answers.
Full details of the penalized honored average standard deviation algorithm are given in Appendix B.

NRTK Residuals Contribution.
NRTK residuals generated by the network and baseline data processing are considered as state vectors and will be to construct the state diagram.The procedure is defined as follows: (i) choose the time window Δ = 10 seconds; (ii) compute the standard deviation of the residuals  res .
Figure 18 shows the concept; (iii) choose a suitable strategy to map the computed values of  res .

NRTK Reliability Results
The aim of this section is to present the results from the analysis.The level of performance ranging from the perfect functioning level perfect down to the complete failure level catastrophic shall be determined from the data.

Horizontal Dilution of Precision (HDOP).
The geometry of the visible satellites is considered as an important factor in achieving high quality results especially for point positioning and kinematic surveying.Anyway, the geometry changes with time due to the relative motion of the user and satellites.
A measure of the instantaneous geometry is the dilution of precision (DOP) factor.The DOP values are computed from the variancecovariance matrix in the ECEF coordinate system and converted to the topocentric local coordinate system with its axes along the local north, east, and up (i.e., vertical) by rotational matrix  by applying the law of covariance propagation.
The DOP value can be defined in various ways; PDOP value in the local system is identical to the value in the global system.In addition to the PDOP, two further DOP definitions are used; HDOP, the dilution of precision in the horizontal position, and VDOP, denoting the corresponding value for the vertical component.
Acceptable horizontal DOP value is HDOP ≤ 2.0. Figure 14 shows the computed HDOP for the analyzed data set.In addition, Figure 15 shows the viewed satellites.Clearly, the number varies between 8 and 10 satellites.

Rover Level of Performance
Prediction.We will predict the rover level of performance ranging from the perfect functioning level perfect down to the complete failure level catastrophic.Three classification lines are chosen in order to separate the computed average standard deviations using (35) into four decisions regions based on the values of  pred ∈ {0.2, .03,.4}. Figure 16 shows the concept.
The computation of  pred using (35) proceeds as follows: (i) Based on the sliding window size Δ, form a data matrix from the baseline residuals.Compute the  variance-covariance matrix VCM res and get and sort in ascending order the diagonal elements .
(ii) Compute the averaged standard deviation, the first term of ( 35)  (iv) Compute the second term of ( 35) where   is the HDOP computed from the solution variance-covariance matrix.The computed HDOP values are shown in Figure 14.

NRTK Reliability Validation
The aim of this section is to introduce the procedures used to validate the NRTK system reliability.The level of performance ranging from the perfect functioning level perfect down to the complete failure level catastrophic shall be determined from data.

NRTK State Diagram Definition.
As we mentioned in the introduction, the key of the NRTK method is the measurement of the distance-dependent errors.The variations in the ionospheric and tropospheric fields are assumed constant in a period of time Δ < 10 s.Since we have access to the rover position error in the topocentric coordinates system (  ,   ,   ), the state diagram is computed as follows:     /  , where   is the number of times we are visiting the state  and   is the total number of events.We see that each row of P (50) has almost identical entries; this confirms that (41) converges to some values as  → ∞.It seems that the existence of a limiting probability that the process  will be in state  after a large number of transitions , and the value is independent of the initial state.

States Transition
For this data set, the probability of being in state 4 is 33.12%, state 3 is 54.59%, and state 2 is 8.0%.

State Diagram.
Information from the network (Section 3.2) and the baseline (Section 3.3) is combined in such way that a single judging number is mapped to the performance levels with states  = {0, 1, 2, 3, 4} and the transitions probabilities  | are obtained.This task is accomplished by counting the frequencies and computing the associated probability, that is,   =   /  , where   is the number of times we are visiting the state  and   is the total number of events.
In order to define a single judging number, various schemes are considered, for instance, an averaging scheme or to assign different weights to each component.

Discussions and Conclusion
An improvement of the rover position accuracy can be achieved by applying procedures for multistate reliability analysis at the system and user level in NRTK.More concretely, the network corrections, baseline residuals, and the associated variance-covariance matrices are considered as the system states and have a direct influence on the rover position accuracy.
The use of the multistate reliability analysis will help us to get some concrete answers to the following problems (i) can we trust the corrections provided by NRTK to the user?
(ii) at which level?
(iii) what are the amelioration potentials?
The weaknesses and the strengths of the system have to be identified and the amelioration potential can be achieved by modifying the serial critical components coupling into paralleled one with a cost effectiveness.
The methods tested make it possible to identify the NRTK critical component with bad data so this can be eliminated or downweighted in the positioning process leading to an improvement in the rover position from epoch to epoch.
It is expected that the suggested approach will reduce the number of wrong or inaccurate rover positions encountered by NRTK users in field, which subsequently will lead to a more efficient work flow for NRTK users.(i)  obs, is total number of satellites used in the estimation process.(ii)  rej, is total number of satellites rejected by the algorithms.(iii)  const is user defined parameter.Default is set to 8. (iv)  sys is total number of satellites in the GNSS constellation.For the GPS, the value is set to 31.
( where   is the window size and is user defined.Default is set to 10. (4) Operation level: the second term in (B.1) shall operate on the same level as the first term and shall have the same unit.This is accomplished by adjustment of parameters ,   , and   .
(5) Stabilization factor  can be implemented by the Danish method.

C. Reliability Computation Technique
In this section we present how the computation of the reliability of parallel structure is carried out step by step.In case of more than two components, the computation of the reliability function (Φ()) is straightforward.The first step is to divide the whole system components into two main components and applying the formula given by (C.2).The next step is to substitute each individual reliability function with the corresponding terms.The last step is simple calculations.

2 JournalFigure 1 :
Figure 1: NRTK processing chain.The left panel shows NRTK module that produces the network corrections.The curved line indicates the output generated by the network module.Middle panel shows the local and the baseline processing module.Right panel shows the rover position solution module.

Figure 6 :
Figure 6: Block diagram of hardware component of module  ,1 .The left block describes the GNSS receiver and the right is for the GNSS antenna.

Figure 7 :
Figure 7: Network ambiguity resolution module  ,2 .The first block of this figure represents the ambiguity processing algorithms and the second block is the ambiguity quality indicators.

Figure 12 :
Figure 12: Start state of the rover position accuracy.

Figure 13 :
Figure 13: State diagram of the rover position accuracy.

Figure 15 :
Figure 15: Number of satellites used in the computation.Baseline of ∼41 km, year: 2014, DOY: 85. Plotted as red dots connected by blue lines.
(i) Choose the time window Δ = 10 seconds.(ii) For each component, compute the standard deviation of the rover position error   = (  ,   ,   ).
Figure 17   shows the concept.(iii)Based on the user requirements and the computed values of   , map this value to the performance defined levels, namely, the states  = {0, 1, 2, 3, 4}.(iv)Compute transitions probabilities  | from data.This task is accomplished by counting the frequencies and computing the associated probability, that is,   = Standard deviations of rover position error

Figure 17 :
Figure 17: Standard deviations of the rover position errors (  ,   ,   ) in the topocentric coordinate system.Horizontal lines are used for classification.Baseline of ∼1 km, year: 2013, DOY: 152.

1 )Figure 19 :
Figure 19: Test area used in this investigation, from Rogaland region.Composed of 6 reference receivers.

x 1 x 2 Figure 20 :
Figure 20: Block diagram of two parallel components.
1,3 denote uncontrolled cycle-slips in the observations.The algorithm uses the observation combinations for this purpose.The interested reader is referred to [27, pp.95-101].

Table 3 :
Reference receiver coordinates, Euref89 .Trimble NetR9 receivers, tracking GPS, and GLONASS satellite signals.Baselines vary between 35 and 112 km and the height difference between the sites is about 225 m.Tables1, 2, and 3 give a full description of subnetwork while Figure19shows the location of reference receivers.The data used for testing was collected on day of year (doy) 152 in 2013 and doy 85 in 2014, respectively.
) Position domain quality indicator   which is defined as the HDOP Stabilization factor  combines and balances between   and   .A reasonable combination is to binomial/exponential trial: