Next Article in Journal
Polarizer-Free AOTF-Based SWIR Hyperspectral Imaging for Biomedical Applications
Previous Article in Journal
Fault Detection and Classification in MMC-HVDC Systems Using Learning Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Accuracy Positioning Based on Pseudo-Ranges: Integrated Difference and Performance Analysis of the Loran System

National Time Service Center, Chinese Academy of Sciences, Xi’an 710600, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(16), 4436; https://doi.org/10.3390/s20164436
Submission received: 7 July 2020 / Revised: 28 July 2020 / Accepted: 5 August 2020 / Published: 8 August 2020
(This article belongs to the Section Remote Sensors)

Abstract

:
The Long Range Navigation (Loran) system as a backup of the Global Navigation Satellite System (GNSS) is a good choice. The dominant deterioration factors of position accuracy are the pseudo-range measurement errors and the geometric dilution of precision (GDOP). This paper focuses on the algorithm integrated difference with pseudo-ranges to improve the position accuracy. Firstly, the theoretical prediction of propagation delay and raw measurement are compared. The results show that the measured pseudo-range consists of a constant term and a temporal term, which reflect the propagation situation along the true path. Secondly, a position solution algorithm based on a pseudo-range and difference is presented, exceeding the limit of a single chain. Finally, some simulation tests are implemented utilizing the new proposed position algorithm to verify the differential performance. This method can reduce the GDOP conveniently through increasing the number of transmitters. In view of the amplitude and characteristics of errors in measurement, systematic error and random noise are distinguished and discussed. The absolute accuracy responds to the pseudo-range bias that is different from geometric distance and repeatable accuracy is mainly influenced by random noise. The difference method can improve the absolute accuracy via the correction degree without changing the geometry of the transmitters.

1. Introduction

Accurate position is a significant topic that should be intensively studied. Generally, it is easy to get the required accuracy using the Global Navigation Satellite System (GNSS). However, the vulnerability of GNSS to unintentional and intentional interference signals can be found frequently nowadays [1,2,3]. For national security and economic effectiveness, a reliable and complementary navigation system is needed desperately [4]. The suitability of the Loran for a backup navigation system has been evaluated and reported [5]. Since then, as a leading role, the United States Coast Guard (USCG) has tried to improve the performance of the Loran system dramatically by the modernization of equipment [6], which mainly included the improvement of transmitters and reference clocks. The improved effect and application are also very obvious [7]. Besides the US, the United Kingdom (UK), South Korea and other countries have also introduced much research on alternatives to the GNSS backup system [8,9,10,11]. They all identified Loran as the terrestrial radio navigation system that has the potential to fulfil those requirements. China has also carried out relevant research [12,13,14]. Meanwhile, China’s major national infrastructure project for developing terrestrial Loran began recently. Three new transmitters using the Loran-C mechanism will be built in the west region of China, which will construct a full coverage timing and positioning network. In order to meet the stated performance [15] as a backup, many researchers and organizations mostly focus on the correction of the Additional Secondary Factor (ASF), which still plays an importance role in position and timing [16,17,18,19,20,21]. The results show that ASF varies with time and can be modified by the difference method, which will be used to improve position accuracy. Nevertheless, the effect of difference is not involved or discussed [22]. So, it is necessary to research the different position solutions and performance to upgrade the position accuracy.
As we all know, geometric distance applied to position can be derived by the product of propagation delay and velocity when the propagation path is smooth and homogeneous. When the propagation distance between the transmitter and receiver is known, the Loran signal propagates along the surface of the Earth, whose velocity is slowed by the atmosphere and the surface dielectric properties of the Earth [19,21]. Based on this, a model predicting propagation delay based on relative dielectric and equivalent conductivity is developed, which contains three main components: Primary Factor (PF), Secondary Factor (SF) and ASF [16,23,24]. However, weather experienced along the propagation path influences the refractive of the atmosphere and the ground conductivity and then cause variations in the propagation delay [24,25,26,27,28]. In other words, propagation delay is not constant and changes with time, so it is difficult to calculate the delay of the Loran signal instantaneously. In some work, researchers have divided propagation delay into spatial and temporal components based on actual measurements, which correspond to a pseudo-range via velocity. Some studies about the impact of the temporal component on Loran position, especially on aviation Non-Precision Approaches (NPAs) and maritime Harbor Entrance and Approach (HEA), have also been done [29]. In fact, systematic error and random noise are involved in the space and temporal components. The influence of observation error on position and the influence of processing on error mitigation need to be studied in detail.
This paper focuses on a new position method of the Loran system, based on measurement pseudo-range combined difference. The performance of this way via position accuracy is discussed. Before this, the distinction between predicted propagation delay and the measurement one is presented. A fit formula and statistical calculation about the measured propagation delay are illustrated, which reflect the characteristics of delay and pseudo-range variations over time. Accounting for the features of pseudo-ranges, some simple tests using new algorithms are simulated. Firstly, a comparison of the Geometric Dilution of Precision (GDOP) between three and four sites is considered, showing that the new method can easily decrease the GDOP. Then, position accuracy is presented in order to estimate the difference performance with the impact of observation error, including SF, ASF and random noise.

2. Pseudo-range and Position Solution Method

2.1. Comparion between Predicted Propagation Delay and Measurement

2.1.1. Predicted Propagation Delay

Generally, it is important to convert the propagation delay to the ground truth range in positioning [16]. If the truth range is known, the propagation delay is derived from the distance divided by the speed of light in a vacuum. However, the Loran signal propagates out radially from the transmitter by groundwave, traveling parallel to the surface of the Earth. So, it does not travel at the free space speed of light, rather, it is slowed by the atmosphere and the surface of the Earth. The groundwave accumulates a delay compared to the expected propagation time over the same distance in free space. This delay is the accumulation of three theoretical components or factors [16,17,18,19,23,24,25]. One part is the P F , which is the time delay of the long wave signal propagation through the atmosphere. The second part is the seawater S F , which is the additional delay due to the signal traveling over the sea. The last one is the A S F , the additional delay due to the signal traveling over the land. The specific formula can be seen in Equation (1).
T p = P F + S F + A S F
where P F = R d n s C , the light speed C in the vacuum is 299,792,458 m/s and R d is the distance of signal propagation between the transmitter and the receiver. The refractive index of the atmosphere n s , assumed as a constant, means that the speed of the signal is a fraction lower than the speed of light in a vacuum. At present, it is difficult to distinguish the S F and the A S F , so S F is often incorporated with A S F as a function of distance R d , carrier frequency f , relative dielectric ε r and equivalent conductivity σ e , as shown in Equation (2):
S F + A S F = a r g W ( f , R d , ε r , σ e ) ω
where ω = 2 π f corresponds to the Loran carrier frequency of 100 kHz and a r g W is the phase of signal attenuation function W about the Loran. It is known that the value of S F + A S F is a smaller modified term comparing to P F [30]. According to the parameters in [30], S F + A S F via distance for each type of ground is given in Figure 1. When the propagation distance increases the S F + A S F also increases. At a constant distance, S F + A S F is the least for average sea water.
In the above model, the total propagation delay is calculated using some simplification, similarly to the result of Brunavs’ formula [22]. The result of the calculation is constant if the ground truth distance is known and the dielectric constant with equivalent conductivity is determined. At this level, the effect of climate and weather is not taken into account, or of ground elevation on the index of refraction of the atmosphere on the surface of the ground. Meanwhile, it is not enough to consider the dielectric and equivalent conductivity changes due to freezing rain and snowy weather. These variations are often lumped together for the user. Real-time factors, such as climate and weather, will result in a distinct difference between predicted value and true propagation delay. That is to say, the total propagation delay is not accurate, although it is possible to predict it from Equations (1) and (2).

2.1.2. Measured Propagation Delay and Assessment

In order to describe the true value, a modified formula is presented:
T p ( t ) = P F + S F + A S F + Δ A S F ( t )
where Δ A S F ( t ) accounts for all of the time varying aspects, reflecting the impact of various real-time factors. In other words,   Δ A S F ( t ) depends on the weather and ground conductivity along the propagation path over time. Any changes in climate and the conductivity would change Δ A S F ( t ) . If the receiver position is not accurate or the propagation path is irregular, the fluctuations would also be contained in Δ A S F ( t ) . Similarly, we also divide the propagation delay into two parts: a spatial component T p , which is the constant part of T p ( t ) , and a temporal component Δ A S F ( t ) .
A better way to get propagation delay is to measure it directly. For this purpose, a Loran receiver is used. The measured propagation delay N is the difference between the Time of Arrival ( T O A ) and Time of Transmission ( T O T ) of Loran pulses coherently modulated on the carrier ground wave. The measurement principle is shown in Figure 2.
τ in T O T ( τ ) represents the transmitter’s frame and r in T O A ( r ) is the receiver time frame. The propagation delay N between receiver and transmitter is defined by Equation (4) below:
N = T O A ( r ) T O T ( τ ) = T p + Δ A S F ( t ) + δ T r + P r + ξ ( t )
where δ T r and P r are the clock bias and the internal processing delays of the receiver, respectively, and ξ ( t ) represents the random residual error due to the unavoidable noise in artificial measurement. In the Loran frequency band, the dominant source of noise is atmospheric noise, which is caused by lighting discharges and man-made noise and local interference from, for example, switch-mode power supplies. Other sources of noise may include transmitter pulse timing jitter or receiver-related noise [31]. It is likely the accuracy of N itself is a function of t , consisting of a constant part N (i.e., N = T p + δ T r + P r ) different from T p and a temporal part Δ A S F ( t ) + ξ ( t ) .
Figure 3 presents a raw measurement of N in a fixed position (109.7503° E, 38.2503° N) for four days, receiving the signal of a BPL transmitter (109.5431° E, 34.9486° N) in Pu Cheng (PC, a city in Shaanxi Province, China). Here, abnormal measurements are removed from the results in advance. The propagation environment is rather homogeneous. Obviously, N varies via time, as seen from the blue line in Figure 3a. A mean value is also plotted by the red line, giving the constant part of N . The average value of the observation is 1253.145 µs, while the predicted constant T p is 1223.860 µs based on the parameters of average land [30] using Equations (1) and (2). Depending upon the navigation application, the difference originating from the effect of δ T r + P r is negligible, because each transmitter is synchronized and the user receiver clock bias is common. Due to the influence of random noise on raw measurements, polynomial fitting is used to derive the trend of observation, seen as the green curve in Figure 3a.
In order to compare the amplitude of each component, Figure 3b gives the decomposition of N , Δ A S F ( t ) and ξ ( t ) . It is clear that N , the mean value over four days, is the determinant of position and Δ A S F ( t ) is much smaller than the former. However, its variation as time advanced is no more than 30 ns no matter whether it increased or decreased. The amplitude of ξ ( t ) is similar to Δ A S F ( t ) , whose histogram of each day is plotted in Figure 4. The standard deviations (STDs) of calculation for ξ ( t ) are 10.343 ns, 13.521 ns, 11.463 ns and 9.647 ns, respectively. Besides, normal distribution fitting curves with the same STDs are also presented, which are used to explain the characteristic of ξ ( t ) .

2.2. Position Solution Method Integrated Difference

The difference method is usually used to correct common errors between the user and referent station. A scheme of the difference correction system in this paper is shown in Figure 5.
Based on the above illustration, the procedure using difference is:
  • Consistency test of reference station receiver and user receiver, to determine the processing delay deviation between two receivers.
  • Identify the accuracy position of a reference station and get the reference station’s pseudo-range, then derive the difference correction based on the accurate location.
  • Obtain pseudo-range observations from at least three transmitters.
  • Calculate the correction of each pseudo-range, then use the least squares solution to find the location.
The details are listed in the following part. The measured propagation delays between the transmitting stations and receiver are converted into pseudo-ranges by multiplying the speed. A pseudo-range’s equation for referent station is given, based on the delay relationship, which is used to calculate the difference correction:
N i ( t ) C = n s R d i + C ( S F i + A S F i + Δ A S F i ( t ) ) + C ( δ T r + P r ) + C ξ i ( t ) .
where i corresponds to the i th transmitter and i = 1 , 2 , 3 Other parameters have the same meaning as above. Here, it is necessary that the local clock of the reference station is consistent with the time system of the transmitter or that the time deviation is known and all the transmitters in one chain are synchronized or the inherent difference is known in advance. The relative difference correction, removing ξ i ( t ) , is expressed as:
C ( S F i + A S F i + Δ A S F i ( t ) ) = N i ( t ) C n s R d i C ( δ T r + P r )
In the previous work, the database of differential correction is necessary, which cost a large number of measurements [19]. In this paper, according to the analysis of S F i + A S F i in Section 2.1.1 and the amplitude of Δ A S F i ( t ) in Section 2.1.2, differential correction for user position is derived depended on distance relation. A simple data map around the referent station is calculated. The pseudo-range integrated difference equation for the user is given as:
N i ( t ) C C ( S F i + A S F i + Δ A S F i ( t ) ) = n s R d i + C ( δ T r + P r ) + C ξ i ( t )
where C ( S F i + A S F i + Δ A S F i ( t ) ) on the left side of Equation (7) is the correction corresponding to different transmitters. It is noted that the environmental features of the path from the transmitter to the reference station are similar to those of the path to the user’s receiver. If the paths are very different, the difference method is not suitable.
In the geodetic coordinate system, the Loran system transmitter’s position ( φ i , λ i ) is known. The receiver position ( φ , λ ) , R d i and δ T r + P r are three unknown parameters to be sought. It is necessary to measure three pseudo-ranges used to at least find a position. The algorithm is similar to the GNSS pseudo-range solution based on all available measurements using least squares. Differing from the former algorithm [22], another processing method for the distance between the transmitter and the user is presented in this paper, and details can be seen in Appendix A. With the help of matrix inversion, the main linearized observation equation is thus given by:
A · X = B
where A is the design matrix, X is the unknown quantity and B is the observation vector eliminating the difference correction. The least squares solution is derived as:
X = ( A T · A ) 1 · A T · B
The assumed position P 0 ( φ 0 , λ 0 , T u 0 ) is updated by the increment ( Δ φ , Δ λ , Δ T u ) and the process is repeated. If the new update is less than 1mm, then no further iterations are considered necessary and we arrive at a “best” Loran fix.

3. Results and Discussion

Raw measurements of propagation delay or pseudo-range contain a number of biases and errors. In order to illustrate the practicability of the pseudo-range method, some detailed simulation tests are discussed in the following parts. The simulation environment is a seawater path, which is partly homogeneous. Another propagation type is considered in the Discussion.

3.1. Geometric Dilution of Precision

In the Loran positioning system, the GDOP between the receiver and the transmitters is a dominant deterioration factor of positioning accuracy. The impact of geometry on the accuracy performance of a ranging system is well understood [32]. An effective method for improving position accuracy is to optimize the geometrical placement of transmitters, thereby decreasing the GDOP. On the other hand, another way to improve performance is by increasing number of stations involved in the position solution when the transmitter station’s location has been confirmed, However, it is difficult to diminish the GDOP using traditional hyperbolic navigation position fixing due to chain restrictions. The pseudo-range method proposed in this paper is convenient to progress two chains or more, not limiting the number of transmitters. As an example, two chains on the east coast of China are used to analyze the GDOP. The details are listed in Table 1.
The GDOP is calculated for three and four sites, respectively, using the geodesic and its azimuth, shown in Figure 6. An area with GDOP ≤ 20 is presented. The left one is for three sites (M, X, Y) and the right one is for four sites (M, X, Y, Z). It is obvious that the coverage area of the three sites is smaller than for the four sites. That is to say, the value of the GDOP for the four sites is smaller than for the three sites on fixed points. Table 2 lists three points within coverage, arbitrarily including their coordinates and corresponding results, which are for later analysis.

3.2. Position Accuracy

When referring to the accuracy of a positioning system, it is necessary to distinguish its absolute accuracy and repeatable accuracy. According the Loran-C User Handbook [33], the absolute accuracy is defined as the accuracy of a position with respect to the geographic or geodetic coordinates of the Earth. The repeatable accuracy, then, is the accuracy with which a user can return to a position whose coordinates have been measured at a previous time with the same navigational system [31,33]. The true test of the position method is the position accuracy between solutions in a receiver and true fix, where the absolute accuracy and repeatable accuracy are contained.
Before discussion, a specification of some simulation parameters is illustrated in Table 3. Three preferential user sites within coverage range are presented, with coordinates the same as in Table 2. The distances between the test site and each transmitter using Equations (A2) and (A3) in the Appendix A, as well as P F , are also given. Similarly, the predicted value of S F and A S F , assuming all-seawater paths according to Equation (2), are calculated and are presented in Table 3. Based on the former two items, the pseudo-ranges not containing time variations and clock biases are simulated, as seen in Table 3. Additionally, four stations are considered regardless of the number of participants. Here, the different correction progress is ignored and the performance of difference is stressed.

3.2.1. Position without Difference Correction

The method without difference is tested. Based on the observation characteristics of the propagation delay referred to in Section 2.1.2, the pseudo-range consists of the distance between the user and transmitter, S F and A S F . Here, the temporal Δ A S F is ignored due to a small amplitude. Meanwhile, random noise with an STD of 50ns (i.e., about 15m, which is smaller than the typical value) is added to sampling point 1,000 of each pseudo-range. Actually, the range of noise determined by the signal-to-noise ratio (SNR) of the received signal is temporarily not considered here. All terms in the pseudo-range are considered as measurement biases besides distance related to position. Figure 7 gives the scatter position error using the algorithm proposed in Section 2.2. The longitude and latitude error are represented by the horizontal and vertical axis, respectively. The upper three subgraphs, with respect to the user’s position A, B, C, are the results of the three sites composed of M, X, Y and the bottom three subgraphs are the results of the four sites composed of M, X, Y, Z. The red star is the statistical average corresponding to the case without random noise. A significant absolute offset of position emerges due to the impact of S F and A S F . It implys that the pseudo-range suffers from large measurement biases resulting in absolute accuracy in the order of hundreds of meters or more, seen in the 95% error radius. The error’s directional components, along longitude and latitude in the form of a statistical average, are listed in Table 4, where the STDs of each component represents the repeatable accuracy. Comparing the upper three subgraphs and the bottom one, the repeatable accuracy is improved because the GDOP for the four sites is smaller than for the three sites, that is, a smaller GDOP suppresses the influence of random noise. However, absolute accuracy does not always improve because it also depends on observation error.

3.2.2. Position with Difference

Based on the results of Section 3.2.1, the difference correction is used to mitigate the influence of S F and A S F to obtain the actual distance between transmitter and receiver. For example, the amplitude of S F and A S F is reduced to 70%, that is to say, 30% of the bias in the pseudo-range is corrected. The results are presented in Table 4. Similarly, 70% of biases are corrected and are listed in Table 4, too. It shows that the latitude and longitude direction’s positioning average error decreased proportionally and gradually with the increase in correction. Compared with no difference, the latitude average error is reduced to 873.035 m and 370.236 m, respectively, corresponding to 30% and 70% correction for the A site. However, the STD is almost unchanged because the random noise is not corrected. Ideally, difference would correct all the system errors, implying that there is no S F or A S F in the observations. The absolute error is almost nearly zero for both three and four transmitters, as seen in Table 4. Figure 8 presents the scatter plot of position error by the influence of random noise only, using the algorithm proposed in Section 2.2 with difference. The longitude and latitude errors are represented by the horizontal and vertical axis, respectively. The upper three subgraphs are the results of three sites, composed of M, X, Y. The bottom three subgraphs are the results of four sites, composed of M, X, Y, Z. There is almost no position bias from the statistical mean results, which means a good absolute accuracy performance. The 95% error radius for the B site is smaller than for the A and C sites due to its smaller GDOP. Comparing the upper and bottom subgraphs, the 95% error radius decreases for the bottom case, which also responds to the decreased GDOP for the four transmitters. In a practical application, the difference correction is not always or necessarily proportional, so the correction effect based on arbitrary correction is not improved at all times.

3.2.3. Discussion

In Section 3.2.1 and Section 3.2.2 of the paper, the main influence of S F + A S F and random noise on position accuracy is discussed. The effect of Δ A S F with a smaller amplitude is omitted because its influence is similar to S F + A S F , which is system bias corresponding to position. It is shown that difference correction can calibrate the system error only, improving the absolute accuracy. However, if the correction scale of the system error in the observed pseudo-range is inconsistent, i.e., all the observation is non-weighted, then a process method considering unequal precision is very necessary. The least squares solution equation with weight could be expressed as follows:
Z = ( A T · W · A ) 1 · A T · W · b
Naturally, it is important to determine the weighting matrix. One method is based on the SNR of each signal, seen reference [19] for details. As for random noise, it is eliminated significantly by filtering or smoothing the pseudo-range. The repeatable accuracy is improved using a pseudo-range after filtering, such as the Kalman filtering progress. When the proposed pseudo-range method is used, the user receiver needed to obtain three transmitter signals at least. The position of the new transmitter and geometric placement must be considered before using this new method in order to ensure signal reception.
In the above discussion, the simulation path is seawater, where S F + A S F is the smallest one of all the ground types shown in Figure 1. Clearly, in another propagation type, such as average land, the absolute position bias will be much larger than that for seawater. The difference method is necessary in order to get a higher absolute accuracy. Compared with the actual difference, the degree of correction depends on the correlation of the two paths. One is the path between the transmitter and the reference station. Another one is the path between the same transmitter and the user receiver. A more complicated situation is that the degree of correction is different for all transmitters. A more detailed analysis should be considered in future work.

4. Conclusions

In this paper, a theoretical predicted propagation delay of the Loran signal and influential factors are discussed. Moreover, the direct measurement of propagation delay between user and transmitter are illustrated. Both theoretical analysis and measured results showed that propagation delay could be divided into a spatial constant term and temporal term. In particular, we presented a new position solution algorithm combining difference using the Loran system conveniently based on pseudo-range measurement. Some conclusions about this method are:
The new algorithm, not limited by the chain, is better than the traditional hyperbolic method. Besides, it is capable of decreasing the GDOP because of greater participation of the transmitters in the position solution. The repeatable accuracy could be improved without changing the amplitude of random noise.
In addition, the absolute accuracy corresponded to the measurement’s bias in the pseudo-range, such as S F and A S F . So, in the new method using difference, it is convenient to eliminate the measurement bias of the pseudo-range to improve the absolute position accuracy due to the clear error sources. In detail, absolute positioning offset decreased proportionally with the increase in correction degree about the pseudo-range.
The work in this paper could be applied to the Loran receiver with a filter restraining random noise, thereby improving the repeatable accuracy of the Loran system further. It is able to promote the safety of HEA and the development of the Loran receiver. In the future, the Loran system may create greater value within the navigation field.

Author Contributions

B.Y. gave the position solution algorithm with integrated difference and carried out the simulation test and drafted the manuscript. Y.L. contributed to the discussion of the propagation delay, especially the influence of climate. W.G. provided the scheme of actual measurement and corresponding date. Y.H. supervised the work, arranged the architecture and contributed to the writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the Natural Science Foundation of China for Distinguished Young Scholars, grant number 11803040 and Frontier Science of the Chinese Academy of Sciences, grant number QYZDJ-SSW-JSC034.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Position Solution Algorithm

The Loran is a navigation system that operates at 100 kHz. The Loran uses a hyperbolic positioning technique to estimate one’s location [34]. Then, pseudo-ranges like GNSS are developed [21].
Different from the former algorithms in [22], another processing method is presented in this paper. The pseudo-range equation is expressed as:
N i ( t ) C = n s R d i + C ( S F i + A S F i + Δ A S F i ( t ) ) + C ( δ T r + P r ) + C ξ i ( t )
The algorithm is based on a GNSS-like pseudo-range solution of all available measurements using least squares. In order to describe it clearly, C ( S F i + A S F i + Δ A S F i ( t ) ) is neglected due to its smaller amplitude, however, its influence is reflected in the measurement.
In the geodetic coordinate system, the Loran system transmitter’s position ( φ i , λ i ) is known, and the receiver position ( φ , λ ) , R d i and δ T r + P r are three unknown parameters to be sought. The distance between the transmitter and receiver on an ellipsoid could be calculated by many formulas. A universal formula used in the field of navigation is the Andoyer–Lambert formula [35]:
R d i = S 0 i + Δ S i
where:
{ S 0 i = a δ 0 i c o s δ 0 i = sin φ i s i n φ + c o s φ i c o s φ c o s ( λ λ i ) Δ S i = a f 4 [ s i n δ 0 i δ 0 i 1 + c o s δ 0 i ( s i n φ + s i n φ i ) 2 s i n δ 0 i + δ 0 i 1 c o s δ 0 i ( s i n φ s i n φ i ) 2 ]
S 0 i is the arc distance of a sphere, a is the radius of the Earth, δ 0 i is the geocentric corner between the transmitter and receiver, Δ S i is a correction due to the sphere to ellipsoidal surface, f = ( a b ) / a is the oblateness of the ellipsoid and a and b refer to semi-major and semi-minor axes, respectively, of the WGS-84 ellipsoid.
It is necessary to measure three pseudo-ranges at least used to find a position. So, the process is started by an assumed position P 0 ( φ 0 , λ 0 , T u 0 ) and Equation (A1) had to be linearized using a Taylor expansion:
{ ρ 1 n s R d 1 ( φ 0 , λ 0 , T u 0 ) R d 1 φ Δ φ | P 0 + R d 1 λ Δ λ | P 0 + R d 1 T u Δ T u | P 0 ρ 2 n s R d 2 ( φ 0 , λ 0 , T u 0 ) R d 2 φ Δ φ | P 0 + R d 2 λ Δ λ | P 0 + R d 2 T u Δ T u | P 0 ρ 3 n s R d 3 ( φ 0 , λ 0 , T u 0 ) R d 3 φ Δ φ | P 0 + R d 3 λ Δ λ | P 0 + R d 3 T u Δ T u | P 0
where ρ i = N i C and R d i represents the distance between the transmitter and the assumed position. With the help of matrix inversion, the linearized observation equation is thus given by:
A · X = B
where
A = [ R d 1 φ R d 1 λ R d 1 T u R d 2 φ R d 2 λ R d 2 T u R d 3 φ R d 3 λ R d 3 T u ] ,   X = [ Δ φ Δ λ Δ T u ] , B = [ ρ 1 n s R d 1 ( φ 0 , λ 0 , T u 0 ) ρ 2 n s R d 2 ( φ 0 , λ 0 , T u 0 ) ρ 3 n s R d 3 ( φ 0 , λ 0 , T u 0 ) ]
The most important process is to get the partial derivatives in matrix A .
R d i φ = a δ 0 i φ + a Δ δ i φ
where
δ 0 i φ = cos φ i sin φ cos ( λ λ i ) sin φ i cos φ 1 [ sin φ i sin φ + cos φ i cos φ cos ( λ λ i ) ] 2
However, Δ δ i is very complex, which is a function of δ 0 i ,   φ and λ . The partial differential of φ is given by:
Δ δ i φ = Δ δ i δ 0 i · δ 0 i φ + Δ δ i φ
The right part is expressed:
{ Δ δ i δ 0 i = f 4 [ [ ( c o s δ 0 i 1 ) ( 1 + c o s δ 0 i ) + ( s i n δ 0 i δ 0 i ) s i n δ 0 i ] · ( s i n φ + s i n φ i ) 2 ( 1 + c o s δ 0 i ) 2 ] f 4 [ [ ( c o s δ 0 i + 1 ) ( 1 c o s δ 0 i ) ( s i n δ 0 i + δ 0 i ) s i n δ 0 i ] · ( s i n φ s i n φ i ) 2 ( 1 c o s δ 0 i ) 2 ] Δ δ i φ = f 4 [ ( s i n δ 0 i δ 0 i ) [ 2 ( s i n φ + s i n φ i ) ] c o s φ 1 + c o s δ 0 i ] f 4 [ ( s i n δ 0 i + δ 0 i ) [ 2 ( s i n φ s i n φ i ) ] c o s φ 1 c o s δ 0 i ]
Besides,
R d i λ = a δ 0 i λ + a Δ δ i λ
where
δ 0 i λ = cos φ i cos φ sin ( λ λ i ) 1 [ sin φ i sin φ + cos φ i cos φ cos ( λ λ i ) ] 2
Similarly, the partial differential of Δ δ i   a b o u t λ is given by:
Δ δ i λ = Δ δ i δ 0 i · δ 0 i λ + Δ δ i λ
Now, all components in matrix A are listed in Equations (A5) and (A6). Based on the expression of the matrix, the least squares solution is derived as:
X = ( A T · A ) 1 · A T · B
The assumed position P 0 ( φ 0 , λ 0 , T u 0 ) is updated by the increment ( Δ φ , Δ λ , Δ T u ) and the process is repeated. If the new update is less than 1mm, then no further iterations are considered necessary and arrive at a “best” Loran fix. The whole process of position solution based on pseudo-ranges is clear. The accuracy of the derived position solution can be calculated by comparing the ground truth fixes.

References

  1. Grant, A.; Williams, P.; Ward, N.; Basker, S. GPS Jamming and the Impact on Maritime Navigation. J. Navig. 2009, 62, 173–187. [Google Scholar] [CrossRef] [Green Version]
  2. Last, D. GNSS: The Present Imperfect. Inside GNSS 2010, 5, 60–64. [Google Scholar]
  3. Zhao, X.; Zhan, X.Q.; Liu, X.; Li, S.J. GNSS Vulnerability Analysis and Assessment. J. Aeronaut. Astronaut. Aviat. 2014, 46, 12–20. [Google Scholar]
  4. Kim, H.; Lee, J.; Oh, S.H.; So, H.; Hwang, D.H. Multi-Radio Integrated Navigation System M&S Software Design for GNSS Backup under Navigation Warfare. Electronics 2019, 8, 188. [Google Scholar] [CrossRef] [Green Version]
  5. Volpe, J.A. Vulnerability Assessment of the Transportation Infrastructure Relying on the Global Positioning System. National Transportation System Center. August 2001. Available online: https://rntfnd.org/wp-content/uploads/Vople_vulnerability_assess_2001.pdf (accessed on 5 August 2020).
  6. FAA report to FAA Vice President for Technical Operations Navigation Services Directorate. Loran’s Capability to Mitigate the Impact of a GPS Outage on GPS Position, Navigation, and Time Applications. March 2004. Available online: http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.126.7557 (accessed on 5 August 2020).
  7. Offermans, G.; Bartlett, S.; Schue, C. Providing a Resilient Timing and UTC Service Using eLoran in the United States. Navig. J. Inst. Navig. 2017, 64, 339–349. [Google Scholar] [CrossRef]
  8. Johnson, G.W.; Siszek, P.F.; Hartnett, R.J.; Shalaev, R.; Wiggins, M. An Evaluation of eLoran as a Backup to GPS. In Proceedings of the 2007 IEEE Conference Technologies for Homeland Security, Woburn, MA, USA, 16–17 May 2007; pp. 95–100. [Google Scholar]
  9. Basker, S. Development of a Radio Navigation Plan for General Lighthouse Authorities of the United Kingdom and Ireland. J. Navig. 2007, 60, 11–16. [Google Scholar] [CrossRef]
  10. Yu, D.-H. Study of Alternative Navigation systems for GNSS in South Korea. J. Inf. Commun. Converg. Eng. 2010, 8, 524–527. [Google Scholar] [CrossRef]
  11. Griffioen, J.W.; Oonincx, P.J. Suitability of Low-Frequency Navigation Systems for Artillery Positioning in a GNSS Denied Environment. J. Navig. 2013, 66, 35–48. [Google Scholar] [CrossRef] [Green Version]
  12. Hu, A.P.; Gong, T. Research Status and Progress on the Enhance Loran-C Navigation Technology. Mod. Navig. 2016, 1, 74–78. [Google Scholar]
  13. Yan, B.R.; Li, Y.; Guo, W.; Hua, Y. A Method of Long-wave Timing Based on Difference. Acta Astronomica Sin. 2018, 59, 54. [Google Scholar]
  14. Li, Y.; Hua, Y.; Yan, B.R.; Guo, W. Experimental Study on a Modified Method for Propagation Delay of Long Wave Signal. IEEE Antennas Wirel. Propag. Lett. 2019, 18, 1719. [Google Scholar] [CrossRef]
  15. Williams, P.; Basker, S.; Ward, N. E-Navigation and the Case for eLoran. J. Navig. 2008, 61, 473–475. [Google Scholar] [CrossRef]
  16. Williams, P.; Last, D. Mapping the ASFs of the Northwest European LORAN C System. J. Navig. 2000, 53, 225–235. [Google Scholar] [CrossRef]
  17. Siszek, P.; Johnson, G.; Oates, C.; Wiggins, M.; Hartnett, R. Methods for Developing ASF Grids for Harbor Entrance and Approach. In Proceedings of the 35th Annual Technical Symposium, Groton, MA, USA, 24–25 October 2006. [Google Scholar]
  18. Johnson, G.W.; Siszek, P.F.; Hartnett, R.J.; Dykstra, K.; Shalaev, R. Temporal Variation of Loran ASFs and their Effects on HEA Navigation. In Proceedings of the European Navigation Conference ENC2008, Toulouse, France, 22–28 April 2008. [Google Scholar]
  19. Hargreaves, C. ASF Measurement and Processing Techniques to Allow Harbour Navigation at High Accuracy with eLoran. Master’s Thesis, The University of Nottingham, Nottingham, UK, 2010. [Google Scholar]
  20. Yun, S.C.; Sang, W.H.; Sang, R.Y.; Chansik, P.; Sung, H.Y.; Sang, J.L. A Study on the ASF Correction Age and Error for Effective eLoran Data Channel Utilization in Korea. J. Korean GNSS Soc. 2013, 2, 109–114. [Google Scholar]
  21. Wang, X.Y.; Zhang, S.F.; Sun, X.W. The Additional Secondary Phase Correction System for AIS Signals. Sensors 2017, 17, 736. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Minimum Performance Standards for Marine eLORAN Receiving Equipment. Radio Technical Commission for Maritime Services. Available online: https://rtcm.myshopify.com/collections/maritime-navigation-equipment-standardsproducts/copy-of-differential-gnss-package-both-of-the-current-standards-10402-3-and-10403-3 (accessed on 11 February 2020).
  23. Wnong, X.; Jinlin, Z.; Wenkui, L.; Li, G. Prediction and Measurement of Loran-C Additional Secondary Factor (ASF). In Future Wireless Networks and Information Systems 2012; Springer: Berlin/Heidelberg, Germany, 2012; pp. 225–232. [Google Scholar]
  24. Sung, H.Y.; Chang, B.L.; Young, K.L.; Jong, K.L.; Young, J.K.; Sang, J.L. Accuracy Improvement Technique for Timing Application of LORAN-C Signal. IEEE Trans. Instrum. Meas. 2011, 60, 2648–2654. [Google Scholar]
  25. Hwang, S.W.; Choi, Y.S.; Lee, S.J. An analysis on ASF variations and temperature. In Proceedings of the IEEE/ION Position, Location and Navigation Symposium—PLANS 2014, Monterey, CA, USA, 5–8 May 2014. [Google Scholar]
  26. Campbell, L.W.; Doherty, R.H.; Johler, J.R. Loran-C System Dynamic Model: Temporal Propagation Variation Study. Electromagnetic Wave Propagation. 1979. Available online: https://xueshu.baidu.com/usercenter/paper/show?paperid=9b630e5dfc60760f6819f512108e9f8a&site=xueshu_se&hitarticle=1 (accessed on 5 August 2020).
  27. Mungall, A.G.; Costain, C.C.; Ekholm, W.A. Influence of Temperature-Correlated Loran C Signal Propagation Delays on International-Scale Comparisons. Metrologia 1981, 17, 91–96. [Google Scholar] [CrossRef]
  28. Ning, L.; Gang, M.; Lachapelle, G.; Elizabeth, C. ASF Effect Analysis Using an Integrated GPS/eLORAN Positioning System. In Proceedings of the ION NTM, Monterey, CA, USA, 18–20 January 2006. [Google Scholar]
  29. Johnson, G.W.; Shalaev, R.; Oates, C.; Siszek, P.F.; Hartnett, R.; Lown, D.; Shmihuk, K.A. Precedure for Creating Optimal ASF Grids for Harbor Entrance & Approach; Institute of Navigation, GNSS 2006: Fort Worth, TX, USA, 2006. [Google Scholar]
  30. ITU-R P. Ground-Wave Propagation Curves for Frequencies between 10 kHz and 30 mHz. 2007, pp. 368–369. Available online: https://www.itu.int/rec/R-REC-P.368/en (accessed on 5 August 2020).
  31. Safar, J.; Vejrazka, F.; Williams, P. Assessing the Limits of Eloran Positioning Accuracy. Int. J. Mar. Navig. Saf. Sea Transp. 2011, 5, 93–101. [Google Scholar]
  32. Kaplan, E.D. Understanding GPS Principles and Application; Artech House: Norwood, MA, USA, 2006; pp. 322–328. [Google Scholar]
  33. U.S. Coast Guard and the U.S. Coast Guard Auxiliary. Loran-C User Handbook. Available online: https://www.loran.org/otherarchives/-1992%20-LoranC%20User%20Handbook%20-%20USCG.pdf (accessed on 30 January 2020).
  34. Qiu, D.; Boneh, D.; Lo, S.; Enge, P. Reliable Location-Based Services from Radio Navigation Systems. Sensors 2010, 10, 11369–11385. [Google Scholar] [CrossRef] [PubMed]
  35. Lambert, W.D. The Distance between Two Widely Separated Points on the Surface of the Earth. J. Ish. Acad. Sci. 1942, 32, 125–130. [Google Scholar]
Figure 1. The trend of S F + A S F via distance.
Figure 1. The trend of S F + A S F via distance.
Sensors 20 04436 g001
Figure 2. The schematic diagram of propagation delay.
Figure 2. The schematic diagram of propagation delay.
Sensors 20 04436 g002
Figure 3. The schematic diagram of measured propagation delay and decomposition. (a) Raw measurement; (b) Decomposition of propagation delay.
Figure 3. The schematic diagram of measured propagation delay and decomposition. (a) Raw measurement; (b) Decomposition of propagation delay.
Sensors 20 04436 g003
Figure 4. Character of observation noise in propagation delay.
Figure 4. Character of observation noise in propagation delay.
Sensors 20 04436 g004
Figure 5. A scheme of the difference correction system.
Figure 5. A scheme of the difference correction system.
Sensors 20 04436 g005
Figure 6. The GDOP of user based on azimuth, the left one is for three sites and the right one is for four sites.
Figure 6. The GDOP of user based on azimuth, the left one is for three sites and the right one is for four sites.
Sensors 20 04436 g006
Figure 7. Scatter plot of position accuracy effected by SF + ASF and random noise. The upper three subgraphs are for three transmitters and the bottom three subgraphs are for four sites.
Figure 7. Scatter plot of position accuracy effected by SF + ASF and random noise. The upper three subgraphs are for three transmitters and the bottom three subgraphs are for four sites.
Sensors 20 04436 g007
Figure 8. Scatter plot of position accuracy with complete difference, only including random noise in the pseudo-range. The upper three subgraphs are for three transmitters and the bottom three subgraphs are for four cases.
Figure 8. Scatter plot of position accuracy with complete difference, only including random noise in the pseudo-range. The upper three subgraphs are for three transmitters and the bottom three subgraphs are for four cases.
Sensors 20 04436 g008
Table 1. Location of sites used in the calculation.
Table 1. Location of sites used in the calculation.
FlagLocationLatitude (N)Longitude (E)
MRongCheng37.0644122.3228
XHeLong42.7199129.1075
YXuanCheng31.0689118.8860
ZRaoping23.7239116.8958
Table 2. The GDOP of user sites.
Table 2. The GDOP of user sites.
Test SiteLatitude(N)Longitude(E)GDOP
Three SitesFour Sites
A25.000125.00018.4763.721
B35.000130.0004.2952.774
C25.000135.00018.2926.946
Table 3. Specification of simulation parameters, ns = 1.000315.
Table 3. Specification of simulation parameters, ns = 1.000315.
Test SiteTransmitters Distances   ( R d / m ) P F / μ s S F + A S F / μ s Pseudo-Range ( ρ / m )
AM1,361,597.3124,543.2302.7191,362,412.499
X2,001,273.4226,677.6324.1782,002,526.074
Y901,698.6333,008.6901.673902,200.037
Z834,208.6942,783.4971.521834,664.528
BM728,502.8402,430.7891.285728,887.962
X860,486.0812,871.1771.580860,959.626
Y1,125,186.9223,754.4022.1801,125,840.485
Z1,780,238.8645,940.1083.6741,781,340.376
CM1,800,246.826,142.8643.7201,801,362.021
X2,038,381.5886,801.4514.2632,039,659.613
Y1,718,993.1715,735.7503.5351,720,052.800
Z1,841,004.3836,142.8643.8131,842,147.447
Table 4. The position error’s directional components, random noise with STD 50ns.
Table 4. The position error’s directional components, random noise with STD 50ns.
SF + ASFNumberTest SiteLat Error
Average/m
Lat Error
STD/m
Log Error
Average/m
Log Error
STD/m
YesThreeA1240.696191.300−1353.64865.061
B381.3929.958−474.16849.089
C270.754154.478−595.915134.361
FourA371.54430.192−1097.38435.255
B576.35515.359−737.31234.628
C304.61448.458−622.78468.172
Yes 1ThreeA873.035191.346−948.42365.974
B267.89330.632−331.79249.735
C190.103159.826−418.229140.051
FourA260.04430.605−768.59034.547
B403.78216.262−514.94436.179
C212.64046.631−434.99264.4859
Yes 2ThreeA370.236185.145−407.63163.306
B115.49830.396−144.48749.337
C81.238162.273−176.940142.134
FourA111.06429.137−328.44433.520
B173.12216.335−222.15336.548
C88.51747.753−184.27767.344
NoThreeA0.920197.696−2.28867.840
B−0.53829.8991.44050.042
C−10.626150.57610.775133.709
FourA−0.62631.4270.60737.680
B1.26216.771−1.19937.655
C−0.31146.7410.47464.937
1 S F + A S F is decreased 30% in pseudo-range. 2 S F + A S F is decreased 70% in pseudo-range.

Share and Cite

MDPI and ACS Style

Yan, B.; Li, Y.; Guo, W.; Hua, Y. High-Accuracy Positioning Based on Pseudo-Ranges: Integrated Difference and Performance Analysis of the Loran System. Sensors 2020, 20, 4436. https://doi.org/10.3390/s20164436

AMA Style

Yan B, Li Y, Guo W, Hua Y. High-Accuracy Positioning Based on Pseudo-Ranges: Integrated Difference and Performance Analysis of the Loran System. Sensors. 2020; 20(16):4436. https://doi.org/10.3390/s20164436

Chicago/Turabian Style

Yan, Baorong, Yun Li, Wei Guo, and Yu Hua. 2020. "High-Accuracy Positioning Based on Pseudo-Ranges: Integrated Difference and Performance Analysis of the Loran System" Sensors 20, no. 16: 4436. https://doi.org/10.3390/s20164436

APA Style

Yan, B., Li, Y., Guo, W., & Hua, Y. (2020). High-Accuracy Positioning Based on Pseudo-Ranges: Integrated Difference and Performance Analysis of the Loran System. Sensors, 20(16), 4436. https://doi.org/10.3390/s20164436

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

Article Metrics

Back to TopTop