CC BY
46Journal of Siberian Federal University. Engineering & Technologies, 2019, 12(1), 72-80
yflK 621.396.98:629.783
Criteria for Selecting GLONASS Satellites to Reduce the Systematic Error in Determining the Coordinates in the Horizon Plane
Vladimir B. Pudlovsky*
Russian Metrological Institute of Technical Physics and Radio Engineering Mendeleevo, Moscow region, 141570, Russia
Received 23.06.2018, received in revised form 02.07.2018, accepted 13.07.2018
The analysis of the effect of systematic errors in pseudorange measurements on the accuracy of determination the coordinates in the plane of the horizon is presented. The criteria for selecting navigation satellites GLONASS are proposed to minimize the systematic components of the error of the planned coordinates, determined on the basis of the method of least squares. The results of checking the effect of selection of satellites when processing real pseudorange measurements by signals of the GLONASS system are presented.
Keywords: navigation receiver, systematic error, coordinates in the horizon plane.
Citation: Pudlovsky V.B. Criteria for selecting GLONASS satellites to reduce the systematic error in determining the coordinates in the horizon plane, J. Sib. Fed. Univ. Eng. technol., 2019, 12(1), 72-80. DOI: 10.17516/1999-494X-0077.
Критерии выбора спутников ГЛОНАСС для снижения систематической погрешности определения плановых координат
В.Б. Пудловский
Всероссийский научно-исследовательский институт физико-технических и радиотехнических измерений Россия, 141570, Московская область, Менделеево
Представлен анализ влияния систематических погрешностей измерений псевдодальностей на точность определения плановых координат (в плоскости горизонта потребителя). Предложены критерии выбора навигационных космических аппаратов ГЛОНАСС для
© Siberian Federal University. All rights reserved
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0). Corresponding author E-mail address: [email protected]
минимизации систематических составляющих погрешности плановых координат, определяемых на основе метода наименьших квадратов. Приведены результаты проверки эффекта отбора спутников при обработке реальных измерений псевдодальностей по сигналам системы ГЛОНАСС.
Ключевые слова: навигационная аппаратура, систематическая погрешность, плановые координаты.
The accuracy of coordinates and height estimating of navigation receiver only by global navigation satellite system (GNSS) signals without any additional information i.e. in the absolute mode is one of the most important metrological characteristic.
The accuracy of the coordinates is generally estimated in an assumption that pseudorange biases to navigational space vehicles (NSV) are equally distributed according to zero-mean normal law and are not correlated between each other [1, 2]. Dilution of Precision (DOP) is generally used as a main criterion of prior coordinates error estimation for these conditions. However, such pseudorange bias idealization is not always justified.
Even for measurements by NSV signals of one GNSS, mainly of GLONASS, the most significant budget components ofequivalent range error are distributed in a way that diifers from normal and non-zero meao value. It is confirmentbyeife ctiveness of various differential and relative navigation methodo [h, 2].
Foe higtee predawn of oavigatien feceivst tNR] m fn fOeolotf mode one should use methods decrearing the mfluenaeof stysiemetic pseudetange biaseswhhoui eny additional information.
In reveral ai|porithmr ff pseudorange p-ecersing of NR mainPy used ttf Rucsian engineers the mcenong of DOP it oaf uitd fn all nor NSV seleet ing i"ail in view" maCe) Moteover, ooe often applies otgoritamr bated an the feast square method (LSM) for measurements processing on all visible NSV to define the receivec caordinftes.
itis weir ORtcwn that "Pa aim of naaigaiion iask (NT) sohttifn in ihe raeeivoe is the receivers oenaition vector definition x (ceotainares ¡end heiehi ef the recfCro)- on-enna os weiO af tCe estimation ef flie recelvep timercala chiPt reNatively to the system tiree reaia) accoMing ft the eseudorange meeruremenrre!i■ not i"t s than to 4 Nit (/=P...w, i->ni. Vector form of these measurements after linearization about true value x is represented in the following way [1, 2]
where ys = Rp - pR - dp - pseudorange measurement result to p-NS considering the correction by thi onPue of prior range/f No thi NS end the rum oP drftorent modeling corrections dp;
1. Introduction
2. Problem definition
y = Gx + s,
(1)
G
a bi gi 1
a2 b2 g 2 1
[abgl] - the matrix of partial pseudorange derivatives on condition
an bn gn 1
vector; £- vector of measurement errors the mean value of which is generally assumed as r = E[e] and variance mat rix a s C = E[ete].
Thetypical algorithm ofnavigation task solution usingLSM fo r the number of measurements more then 4 is describe d 1 ike this [2]:
i^G^ fVy, (2)
where x - tln^ estimation of conditio n vector (CV) accordiing LoM algorithm.
Equations (1) and (2) allow specifying the error oO CV eslimation in [ht foliewing
x-x=(gtg) GT£.
(3)
Tlhe hypothesii as o:f^€;n assumed on the abstnce ol" byeeemotic ereoro in mnasunements (t), ¡su|p|:>ossite;; the meano uiue s hs zeeo
1- = E [o] = 0.
Taking ioys aonsidsration thg assumption CV estimation gyV be untoased as nvell that can be traced from (3)
A = E[i-x^^^G^^O, (4)
where1 = (Ai A2 eS3 A4]T - f ectoo of mean vahiee of CV istimrtion errors.
][n topocentric (i<ncitli coordrnates syetem hs componrnts have the following meaning: W = E [[ — ec^ J — mean value of longitude staor, W2 EE- -Xi]- meanvalueoflatitude error, A^ = E [e3 - e3 ] - mean value od Xeiglit eeior, w 4 EEe te4 e4
] - mean vaiue of estimation error of NR timescales shift m rektira to GNSS
timescale.
Note also that, through the dmgsnal eiements os 17° mttriceo (GTG) oftrn de^^snoi^n^ the valuesofthe DOP, aor exsmple s n (Ioe iotlohmg form [2]:
GDOP = dr[(GTG)f ' =(qu+q2( + 233 + 244 )
0,5
PDOP = (qu
+ 222
HDOP = (qq + q22 )0 5 , (5)
where GDOP - geometric dilution of precision NT solution; PDOP - position dilution of precision in three-dimensional space; HDOP - horizontal dilution of precision in the horizon plane.
The analysis of the pseudorange bias budget shows that condition r = 0 is rarely met in practice especiallyforshortperiods ofcoordinates definitionwiththehelp ofNR.
It follows that the fulfillment of (4) for r ^ 0 is possible only by determining the elements of the matrix G. Hence the purpose of the paper can be formulated as follows: to solve the NT, it is necessary to determine such a configuration of the NSV that will ensure a minimum
- 74 -
displacement of the estimated CV components, tf is? some infosmation on the components r 3 0.
( G' GÎ)-'
3. Theory
If solatso n (2) exist sthen
Thus, conditooa (4tfoi rf ()can be formally fulfilled if
GTr = [abgl]Tr = 0.
It is evMsnt thst for r di=erent from zero its projection on anes of 4-dimentional space of boordinates So be dcfinnd sannof bo eimultaneoesly zero .
Due trr tOia reason iiO fo bntier Co bo liimled fdrtfet by oonditions finding under which the oomoonsanion of syademalic pseudorango btaa ib poboiMe onO0 for 1tai<r ooordinates in plane i.e. the OcUswrng mua. bo fulfille0
Aj = A 2 = 0. (6)
To ^u^t C4)as a normal eystem
GaGA = G'r,
dds^ = I^JE^ 15 = , (7)
Cnndition (<5) penults in orthtganality of vector r tv twt first bauis vectors of its resolution that is expressed thteugh scalar product
(ara)e ib,a) e 0. (8)
Then one should understand when conditions (8) are fulfilled for such type of r all the components ofwhichdiffersfromO.
First variant is when r components are close regarding absolute value. Then columns 1 and r of system (7) are linearly dependent. Therefore the solution of (7), for instance, according Kramer's rule will result in A = A2 = A3 = 0, i.e. the re won t be any shift on all the coordinates and heig ht. Systematic error A4 ofNR timescale thift regarding GNCS seole wit 1 be the onty Co increase.
Coneidering the natore ol tae mool conslderabCo errort eoinpoiing ,4 so ut^or^irt^^ ert imation budget ravo mare vosiants ot r oom]yfnteroC corselaiion nre Ito iio described:
Ay oomponenOs iiirc; posioive aitsl timilar for i^SV pairs;
rii oonponenOs one posilivo eond proportional io NSV elevation angels yj.
Thue condition (toy (o to ^f s^îo t^iie thofoeiowing way: for oariant A:
ri( a + a.) + --- + rnn{ an_i+an) = 0, (9)
r(bi+b2) + ... + rn+ (bn-\ bbn) = 0;
for va riant B:
j=o. ao
j =1 j =1
The efiression for portial dorivativos aJr 0r in topocentric (iocat) cooadinates system and in the polar coordinate s lysiem i2] is define. in the following way:
as = sin p) -ros ;o,
ds = cos atj ■ cos /j, (11)
whe re a, p ange 1 of NpV eighrhnv [
Convidertng this way of aprtifying f and Vj one can rasifynollce t^^t; it is enough to form NSV pa iris ac cordi nj;^ fo cenrrar in mmerry in the plane (the difference by azimuth on 180° and equal elevation angels fac NSV pair) lo fulfill (9) (varianc A).
Functional dependence rj(yj) (variant B) is specific for pseudorange bias components description at the expense of refraction in ionosphere and troposphere. For instance, to recalculate the value of zenith delay in troposphere ZTDj for slant path Zj it is generally used so-called "mapping function" m(j [3].
Srs^MY.
In the same way one should consider the influence of the slant path while calculating the delay inionosphere in rhe wa]]-known Klebnchart rPgoaithm [2]. fn lite first approximation function m(ys) will be
m (aj )=Vsie (as).
The conclusion is that all systematic errors of zenith delays calculation in ionosphere and troposphere will only increase for slant paths of NSV signal propagation in the proportion close to sin-1(Ys). However, these systematic bias can be considered close value for NSV with the equal elevation angles ys. Thus, conditions (9) and (10) are possible for the same NSV pairs selected according to centralsymmetry.
The analysis of condition (8) regarding (11) explains that these conditions are met not only for the given NSV pairs. It is evident that condition (8) is fulfilled by any of NSV for which ys~90° Consequently, without influence on coordinate's error one can include several additional NSV observed in zenith relative to NR in working constellation if it is provided by orbital NSV group from one or several GNSS.
4. Experiment results
The approbation of the influence of the compensation of equal systematic pseudorange bias in GLONASS NSV pairs was conducted on base of the coordinates calculation in the plane comparison with the use of software package RTKLIB v.2.4.2 [4], where the algorithm of definition of absolute coordinatesandheightwiththeleastsquaremethod waspartiallyimplemented.
As reference information there were used the measurement data by GLONASS ranging code (open civil code) only in L1 frequency band. The measurements with period 1 s are obtained in RINEX
-It-
format version 3.03 from NR PolaRx5 (Septentrio), with the help of antenna positioned on the territory of FSUE "VNIIFTRI" in the point with given coordinates.
For further processing the data was used that was recorded during 24 hours dt. 15.03.2018.
NSV selection parameters were formed regarding the following acceptable deviation from satellites search conditions: not more than 10° on elevation angle; not more than 30° on azimuth.
For the location of the NR antenna for the day 15.03.2018, the selection conditions for two pairs of NSV at the viewing angles were made for several time intervals (lasting more than 5 minutes each) with a total duration of more than 400 min. In this selection, an additional condition was used: the difference in the azimuth angle of the viewing vector for NSV from different pairs should exceed 40°.
Moreover one can notice 2 and sometimes 3 pairs of satellites centrally-symmetric regarding NR as well as generally one satellite close to zenith (y>75°). The total sum of GLONASS NR satellites amounts from 8 to 10.
To assess the probability of fulfilling the selection condition for "centrally symmetric" NSVs at other geographic locations, the above-mentioned conditions for the mutual sighting of the GLONASS satellites were simulated for an interval of 8 days. The results of calculating the probability of selection of at least two pairs of NSV are presented in Table 1.
The efficiency of only centrally-symmetric NSV choice for plane coordinates definition was compared to errors received during processing of all visible GLONASS satellites with elevation angles y>5°. The constellations from 4 NSV (2 pairs) has been selected as well as 5 NSV (2 pairs and one NSV close to zenith).
Fig. 1 shows the location example of 10 GLONASS satellites in relation to NR at 00:00 (UTC) for selecting two pairs (R02, R08, R16, R23 numbers) and zenith NSV (R01).
Some results of calculations of the error in determining the coordinates in the plan using the RTKLIB program package using the LSM algorithm are presented in Fig. 2 and in Table 2.
Fig. 2 represents the results of coordinates definition with the LSM during the period of 00:0000:05 (UTC). Blue dots mark the coordinates received according to 9 NSV signals considering measurements correction by a model. Red dots mark the coordinates received according to 5 selected
Table 1. Probability of observation of two pairs of satellites GLONASS
Coordinates of the point, degrees latitude, longitude The total observation interval of two pairs, min Probability of observation of at least two pairs in 8 days
56,02 N 37,22 E 2919 0,25
59,93 N 30,30 E 4597 0,40
68,97 N 33,09 E 4851 0,42
44,50 N 34,16 E 3333 0,29
56,02 N 92,87 E 2802 0,24
80,0 N 100,0 E 3053 0,26
Fig. 1. The view of GLONASS satellites at 00:00 (UTS)
i
— 3 IS V
f 1 1 ' ••• * * y ? *** i ■»
^ wv 1 rn t
1- —1
Fig. 2. Coordinates in the horizon plane according to pseudorange measurements by 5 (red dots) and 9 (blue dots) of GLONASS satellites
NSV (two pairs and one NSV in zenith) without using measurements correction by a model. The point coordinates accepted as standard are shown as intersection of grid full lines.
In the first column of Table 2 the number of NSV 4, 5 or 6 means applying the procedure of satellites pair selection by central symmetry criterion plus one in zenith (for 5 selected NSV).
There were used all pseudorange measurements available at the given time to define the coordinates for NSV number 9 or 10.
Symbols c1 and c2 stand for standard deviation relating the mean values A1 and A2 of error of coordinates under definition.
The capabilities of the RTKLIB package allow you to enable or disable models for correcting pseudoranges due to the influence of the troposphere and the ionosphere when processing measurements conducted in the one frequency band, which is also reflected in Table 2 in the column "Correction by models".
Table 2. Coordinate errors
The number of NSV Correction by models DOP Longitude error, m Latitude error, m Time (UTC), h:min
PDOP HDOP A1 01 A2 02
9 On 1,5 1,0 -3,65 0,42 1,41 0,23 00:00
9 Off 1,5 1,0 0,16 0,20 -0,41 0,27 00:00
4 Off 3,5 1,2 -1,84 0,64 0,78 0,27 00:00
5 Off 2,0 1,2 -0,90 0,46 0,01 0,25 00:00
10 On 1,4 0,8 -2,68 0,40 1,24 0,30 00:05
10 Off 1,4 0,8 -2,48 1,62 -0,27 0,32 00:05
4 Off 3,4 1,2 -1,51 0,49 1,41 0,52 00:05
5 Off 2,0 1,2 -1,19 0,66 0,66 0,55 00:05
9 On 1,8 0,8 2,61 0,23 -1,34 0,45 03:40
9 Off 1,8 0,8 2,64 0,23 -2,42 0,44 03:40
6 On 4,9 1,0 1,81 0,28 -1,71 0,44 03:40
6 Off 4,9 1,0 1,65 0,29 -2,27 0,47 03:40
The values of the DOP are calculated only for the working constellation of NSV, whose pseudorange measurements are used in determining the coordinates.
The analysis of received errors of height definition has not been performed in the course of the work.
The results of coordinates definition error in Table 2 are shown as an example and do not cover the whole scope of processed data received within the period of 24 hours.
5. Discussion of results
For the GLONASS orbital constellation of 24 satellites, the probability of observing two pairs of "centrally symmetric" NSV is from 0.24 to 0.4 or about 400 minutes per day.
The analysis of coordinates definition errors described in Table 2 as well as those received after other measurements have been processed verifies the advantage of NSV selection for partial compensation of defining coordinates shift from the real position of NR antenna.
However it is impossible to fully eliminate the shift of coordinates by NSV selecting. One of the reasons is likely to be the significant difference of pseudorange biases for the pair of selected NSV GLONASS.
As one can see from Table 2 reducing the NSV number the measurement of which were used for coordinate definition can possibly increase the coordinate error dispersion. It is particularly notable during comparison of the values of standard deviation longitude error for operating constellations from 4 and 9 NSV within the period of 00:00 - 00:05. Moreover the mean coordinate value received for selected NSV is as a rule closer to the primary standard.
On the other hand elimination of relatively small NSV number from visible NSV improves the coordinate definition accuracy only to a small extent that might be seen from the errors comparison for operating constellations from 6 and 9 NSV within the period of 03:40 - 03:45 (Table 2). It is easy to understand that in this case the main part of the satellites in an operating constellation has not
changed while pseudorange bias for eliminated and selected NSV is likely to be comparable by the value.
A comparative analysis of the DOP values for a small amount of NSV in the working constellation (about half of the observed satellites) shows that the value of PDOP can not serve as a reliable criterion for the accuracy of the determination of coordinates (in the plane of the horizon). The value of HDOP gives a more adequate estimate of the coordinate error in the plane of the horizon. However, this DOP, calculated in accordance with (5), does not allow to take into account the influence of systematic errors in pseudorange measurements.
Based on the results of analysis of the processed data, it should be separately noted that the compensation of the displacement of the determined coordinates in the in the horizon plane due to selection of the NSV appears more noticeably in the absence (shutdown) of correction of pseudorange measurements by the models included in the RTKLIB software package.
6. Conclusions
Thus the research describes the possibility of pseudorange bias compensation in the pair of NSV during coordinates definition in horizontal plane by the least squares method algorithm. The maximum effect of such compensation can be achieved in the following conditions:
coincidence of absolute value of pseudorange bias in each NSV part;
central-symmetric position of two or more NSV pairs in the projection of spherical coordinates system to horizontal plane (the difference according azimuth to 180° and equal elevation angels for NSV pair);
The usage of NSV located in zenith does not influence the shift of defining coordinates in the plane and allows applying NSV pair with the same elevation angles in different pairs.
The influence of compensation effect of systematic bias is notable during the processing of real pseudorange measurements for selected NSV GLONASS. However, it is impossible to eliminate the shift of defining coordinates.
The use of the DOP values as a criterion for choosing the composition of the working constellation of the NSV does not take into account the influence of systematic errors in pseudorange measurements and does not allow to adequately estimate the accuracy of the coordinates in the horizon plane.
The relevant objective for coordinates definition in absolute mode of NR is to develop new algorithms considering pseudorange bias shift and correlation.
References
[1] GLONASS: Design and Operating Principles, Perov A.I. and Kharisov V.N., Moscow: Radiotekhnika, 4th ed., revised, 2010. 800 p. (In Russian)
[2] Sanz Subirana J., Juan Zornoza J.M., Hernández-Pajares M., GNSS Data Processing Book, Vol. I: Fundamentals and Algorithms, TM-23/1. Noordwijk: ESA Communications, 2013. 223 p.
[3] Pudlovskiy V.B. Tropospheric delay: models and estimation error. The series "Radio communication and radio navigation". Release 3. Radio navigation technology. Perov A.I., Vlasov I.B. Moscow: Radiotekhnika, 2013. p. 50-54. (In Russian)
[4] RTKLIB [Electronic resource] RTKLIB: An Open Source Program Package for GNSS Positioning. - Acess: http://www.rtklib.com/ (date of the application: 14.08.2017).
