Lecture UPC 6.v0.0
Lecture UPC 6.v0.0
Lecture 6
Differential positioning with
Code pseudoranges
Contact: [email protected]
gAGE/UPC
Authorship statement
The authorship of this material and the Intellectual Property Rights are owned by
J. Sanz Subirana and J.M. Juan Zornoza.
This authorship statement must be kept intact and unchanged at all times.
gAGE/UPC
24 April 2014
and these errors vary ‘slowly’ with time. That is, the errors are
correlated on space and time.
Error
gAGE/UPC
Short-baselines
Long-baselines
Pi j i j c ti t j Ti j I i j Ki K j Mi j P ij
Lij i j c ti t j Ti j Ii j i j Ni j bi b j mij L ij
Satellite coordinates r j
deviation r j r j r0j
When approximate values of both
receiver and satellite APC positions
are taken, a linearization around
them yields:
i j 0ij
i j 0ij ρˆ 0ij ri ρˆ 0ij r j
j
ρ̂ 0 i r0j roi
ρˆ ij ri ρˆ 0ij j
User coordinates r0 roi
deviation ri ri r0i
gAGE/UPC
Pi j i j c ti t j Ti j I i j Ki K j Mi j P ij
Lij i j c ti t j Ti j Ii j i j Ni j bi b j mij L ij
Pu j uj c tu t j Tu j Iuj Ku K j P uj
Pr j rj c tr t j Tr j I rj K r K j P rj
rru ru rr
Pi j i j c ti t j Ti j I i j Ki K j Mi j P ij
Lij i j c ti t j Ti j Ii j i j Ni j bi b j mij L ij
j rru ru rr
• Ionosphere ( I i )
j
• Troposphere (Ti )
ru rr
• Wind-up (i )
j
PRN06
PRN06
PLAN-GARR: 15km
PRN30
( L1 ) L1ru
sat
rusat PRN30 ( P1 ) L1ru
sat
rusat
Dif. Instrumental
Dif. Receiver clock: delays and carrier
Main variations Common Dif. Tropo. and Iono. : ambiguities:
for all satellites Small variations constant
Master of Science in GNSS @ J. Sanz & J.M. Juan
8
gAGE
Single-Difference of measurements
(corrected by geometric range!!)
research group of Astronomy and Geomatics
PRN06
PRN06
PLAN-GARR: 15km
PRN30 PRN30
( L1 ) L1ru
sat
rusat ( P1 ) L1ru
sat
rusat
PRN06 PRN06
IND1-IND2: 7m IND1-IND2: 7m
gAGE/UPC
PRN30 PRN30
( L1 ) L1ru
sat
rusat ( P1 ) L1ru
sat
rusat
Master of Science in GNSS @ J. Sanz & J.M. Juan
9
gAGE
Single-Difference of measurements
(corrected by geometric range!!)
research group of Astronomy and Geomatics
PRN06 PRN06
IND1-IND2: 7m IND1-IND2: 7m
PRN30 PRN30
( L1 ) L1ru
sat
rusat ( P1 ) L1ru
sat
rusat
Dif. Instrumental
Dif. Receiver clock: delays and carrier
Main variations Common Dif. Tropo. and Iono. : ambiguities:
for all satellites Small variations constant
Master of Science in GNSS @ J. Sanz & J.M. Juan
10
gAGE
Linear model for Differential Positioning
Code and carrier measurements
research group of Astronomy and Geomatics
Pi j i j c ti t j Ti j Ii j Ki K j P ij
Lij i j c ti t j Ti j Ii j i j Ni j bi b j L ij
where: j j ρ r j
i 0i
ˆ 0ij ri ρˆ 0ij r j
j
ru
j
0 ru ρˆ ru ρˆ rr ρˆ r ρˆ r
j
0u
j
0r
j
0u
j j
0r
j receiver Ref.
Station
gAGE/UPC
Baseline
0 ruj ρˆ 0uj rru ρˆ 0 ru
j
rr ρˆ 0 ru
j
r j rru ru rr
ru rr
gAGE/UPC
gAGE
Linear model for Differential Positioning
and where the rru estimate will be affected in
research group of Astronomy and Geomatics
ε site
Range error due to Range error due to User rru ru rr
reference station Sat. coordinates receiver Ref.
Sta.
coordinates uncertainty uncertainty
u
j
u
j To Earth Centre
b
ρˆ 0 ru
j
ε site ε site
Hint:
ρˆ 0 ru
j
ε site ρˆ 0uj ρˆ 0 rj ε site rj
ρ0j u ρ0j r b
ε site ε site
ρˆ uj
b r rru ρˆ rj
User ru
Refe
receiver Stati
Note: the following approaches have been taken: rru ru rr
j ρ j
ρ j
ρˆ 0 r
0r 0r
j
j
0jr 0ju 0r
ru rr
ρˆ j ρ0u ρ
j j
gAGE/UPC
0u
0u 0uj j
broadcast
Satellite
ephemeris
location error
True position ε
ρuser
user
Differential range error due to
satellite obit error
ρ ref
ε
ref
user ρuser ρref
ε ε
user ref
user
A conservative bound:
b
User
with a baseline b 20km
gAGE/UPC
Reference Station 20 1
20000 1000
CREU and ε
EBRE user
Position
from
Satellite
broadcast
Differential range error from
location error
ephemeris between CREU and EBRE
True position ε
ρuser
user
ρ ref
ε
ref
user ρuser ρref
ε ε
user ref
user
gAGE/UPC
288 km of baseline
User
Reference Station
CREU-EBRE
CREU Differential positioning
Absolute positioning
gAGE/UPC
288 km of baseline
eph
1 2
2
1
1 1
gAGE/UPC
1 2 1 2 constant 0
eph
2 2 eph
eph
1 2
2
1
1 1
gAGE/UPC
1 2 1 2 constant 0
eph
2 2 eph
1
where a b cos
u ρˆ bˆ ρˆ bˆ ρˆ T ρˆ - ρˆ ρˆ T bˆ
I bˆ ρˆ ρˆ T bˆ I ρˆ ρˆ T bˆ
2
Note: in this 3D problem is NOT the elevation of ray.
Note: u sin uˆ
( B A ) 2 a b sin b
a a ε T uˆ ε T sin uˆ
2 2
b sin
ε T I ρˆ ρˆ T
1 b Where: b = b bˆ
b sin
Master of Science in GNSS
Note: ε ρ @ J. Sanz
is the baseline
& J.M. Juan vector
gAGE
ORBIT TEST :
Broadcast orbits
Along-track Error (PRN17)
research group of Astronomy and Geomatics
PRN17:
Doy=077, Transm. time: 64818 sec
rov2
11 km
15.2 km
gAGE/UPC
εT I ρˆ ρˆ T
19.7 km b
rov3
Exercise:
so that the residual error ionospheric, tropospheric and ephemeris are small
compared to the typical errors due to receiver noise and multipath, it can be
assumed: j
Tru I ru 0 ;
j
ˆ j j ρ0ru εeph 0
Working with smoothed code, a residual error of about 0.5 metres could be
tolerable, but for carrier based positioning it should be less than 1 cm to
allow the carrier ambiguity fixing.
• The typical spatial gradient of the ionosphere (STEC) is 1-2 mm/km (i.e.
gAGE/UPC
0.1-0.2 m in 100km), but it can be more than one order of magnitude higher
when the ionosphere is active.
• The correlation radio of the troposphere is lower than for the ionosphere. At
10km of separation the residual error can be up to 0.1-0.2 m. Nevertheless,
90% of the tropospheric delay can be modelled and the remaining 10% can
be estimated together with the coordinates (for high precision applications).
For distances beyond a ten of kilometres or significant altitude difference it
would be preferable to correct for the tropospheric delay at the reference
station and user receiver.
Carrier-smoothed code:
Pseudorange code measurement errors due to receiver noise and multipath
can be reduced smoothing the code with carrier measurements.
Smoothed codes of 0.5m (RMS) can be obtained with 100 seconds smoothing.
On the other hand, the ionospheric error is substantially eliminated in
gAGE/UPC
differential mode and the filter can be allowed for lager time smoothing
windows.
The left hand side of previous equation can be spitted in two terms:
one associated to the reference station and the other to the user:
gAGE/UPC
Pruj 0 ruj Pu j 0 uj Pr j 0 rj
Master of Science in GNSS
rru ru rr ru ε site @ J. Sanz & J.M. Juan
gAGE Differential code based positioning
Pruj 0 ruj ρˆ 0 uj ru c tru Kru P ruj
research group of Astronomy and Geomatics
The left hand side of previous equation can be spitted in two terms:
one associated to the reference station and the other to the user:
Pruj 0 ruj Pu j 0 uj Pr j 0 rj
• The term Pr 0 r is the error in range measured by the reference
j j
Position Position
– The user receiver applies the PRC and RRC to correct its own
gAGE/UPC
DGNSS with code ranges: users within a hundred of kilometres can obtain
one-meter-level positioning accuracy using such pseudorange
Master of Science in GNSS
corrections.
@ J. Sanz & J.M. Juan
32
gAGE Differential code based positioning
The user applies this differential correction to remove/mitigate common
research group of Astronomy and Geomatics
errors: P j j PRC j ρ ˆ j r c t j
u 0u 0u u ru P ru
The previous system for navigation equations is written in matrix notation as:
1
T
ρˆ 01u
Pref
1
2
ρˆ
T
Pref
2
0u
1 ru
c tru
n
Pref T
ρˆ 0nu 1
gAGE/UPC
where
Pref j Pu j 0 uj PRC j
Note: we have considered here only code measurements. The carrier based
positioning will be treated next, using double differences of measurements and
targeting the ambiguity fixing.
Master of Science in GNSS @ J. Sanz & J.M. Juan
35
gAGE
Range Differential Correction Calculation
Broadcast Actual SV
research group of Astronomy and Geomatics
SV Position Position
– The user receiver applies the PRC and RRC to correct its own
gAGE/UPC
DGNSS with code ranges : users within a hundred of kilometres can obtain
one-meter-level positioning accuracy using such pseudorange
Master of Science in GNSS
corrections.36
@ J. Sanz & J.M. Juan
gAGE
research group of Astronomy and Geomatics
GODN
USN3 76 m
GODS
ftp://cddis.gsfc.nasa.gov/highrate/2013/
gAGE/UPC
GODS
USN3
S/A=off
DGPS DGPS
GODN
GODN
GODS: reference
gAGE/UPC
GODS: reference
USN3 station
GODS: reference
USN3 station
PERFORMANCE CATEGORY I
Requirements
?
Time to alarm 6 s
80
60
40
20
0
-20
-40
-60
-80 ANALYSIS NOTES
-100
-120 - Data taken from Overlook PAN Monitor Station,
equipped with Trimble SVeeSix Receiver
-140 - Single Frequency Civil Receiver
gAGE/UPC
GPS/GLONASS Satellites:
Time to alarm is from minutes to hours
No indication of quality of service
Health Messages:
GPS up to 2 hours late
GLONASS up to 16 hours late
gAGE/UPC
Continuity:
Less than 10-5 Chance of Aborting
a Procedure Once it is Initiated.
Availability:
gAGE/UPC
NSE
Confidence bound
gAGE/UPC
SBAS
gAGE/UPC
GBAS
DQM: Checks the validity of the GPS ephemeris and clock data for each
satellite that rises in view of the LGF and at each time new navigation
data messages are broadcast.
MQM: Confirms the consistency of the pseudorange and carrier-phase
measurements over the last few epochs to detect sudden step and any
other rapid errors.
MRCC: Examines the consistency of corrections for each satellite across all
reference receivers.
sm-monitor: Helps ensure a Gaussian distribution for the correction error
with zero mean and that the broadcast spr _ gnd overbounds the actual
errors in the broadcast differential corrections.
MFRT: Verifies that the computed averaged pseudorange corrections and
correction rates fit within the message field bounds. This is from [RD-5]
Receiver Noise
70
Artemis
IOR
AOR-E
60
50
Latitude (°)
40
30
gAGE/UPC
20
-40 -30 -20 -10 0 10 20 30 40
Longitude (°)
RIMS
(x 37)
NLES
(x 6)
EWAN
gAGE/UPC
http://egnos-user-support.essp-sas.eu/egnos_ops/egnos_system/system_description/current_architecture
Do-229A, B, C, D
DIRECTION OF DATA FLOW FROM SATELLITE; MOST SIGNIFICANT BIT (MSB) TRANSMITTED FIRST
IODP
Clock IODP
PRN Integrity
GLONASS
Mask Information
DATA
(1) (6)
IODS IODP
IODP
Service Acceleration
Message Information
(7)
gAGE/UPC
(27)
Ionospheric Ionospheric
Mask Corrections
(18) IODI
IODI
(26)
Master of Science in GNSS @ J. Sanz & J.M. Juan
70
gAGE
Message Time-Outs:
Users can operate even when missing Messages
research group of Astronomy and Geomatics
System Latency
master station
Correction time-Out
Corrections
UTC Timing Data 12 300 None None
Almanac Data 17 300 None None
(*) Fast Correction Time-Out intervals are given in MT7 [between 12 to 120 sec]
Master of Science in GNSS @ J. Sanz & J.M. Juan
72
gAGE
research group of Astronomy and Geomatics
PRN MASK (MT01)
Bit No 1 2 3 4 5 6 . 38 . 120 . 210
Value 0 1 0 1 1 0 1 1 0
PRN GPS GPS GPS GLONASS AORE
PRN 2 PRN 4 PRN 5 Slot 1 PRN 120
PRN mask 1 2 3 21 29
Number
Each MT01 contains its associated IODP
Assignment
Up to 51 satellites in 210 slots. PRN
Slot
120-138 GEO/SBAS
number, between 1 to 51. 139-210 Future GNSS/GEO/SBAS/Pseudolites
• Primarily Removed SA
– Common to ALL users
– Up to 13 Satellites Per Message
– Pseudorange Correction /confidence Bound
– Range Rate Formed by Differencing
– UDRE degrades Over Time
• Acceleration Term in MT 7
gAGE/UPC
(t ) PRC n RRC n (t t n )
PRC (t Y C1 PRC * t sat dt sat
PRC n PRC o
RRC n TGD IONO TROP
t n to
DIRECTION OF DATA FLOW FROM SATELLITE; MOST SIGNIFICANT BIT (MSB) TRANSMITTED FIRST
( RSSUDRE 0 MT10)
N75
N65
N55
N50
S50
0 1 2 3 4 5 6 7 8
S55
S65
gAGE/UPC
S75
S85
MT10
IONO Corrections
+ MT18
IONO GIVE
+ MT26
Degradation Param.
s s s
2 2
flt
2
UIRE s 2
air s 2
tropo
T
x
X N , E, U , cdt
d N2 d NE d NV d NT
d NE d E2 d EV d ET
G T WG
2
d d
1
Px dN d
2 2 2 2
d V2 HPL 6.00 N d NE
E E 2
d NV d EV dVT
d 2
NT d ET dVT d T2 2
W PY1
VPL 5.33 dV
s 12 0
Py
0 2
sN s i2 s i2, flt s i2,UIRE s i2,air s i2,tropo
-Ks-s s Ks
gAGE/UPC
X ~ N (0,1)
P( X 5.33) 107
s
fc rrc ltc er , if
2
RSSUDRE 0 ( MT10)
research group of Astronomy and Geomatics
t tu tlat
2
MT25
fc a 2
MT2-6,24 tltc , v0 or v1
IODFcurrent , IODFprevious 3
0 , ( IODF IODFprevious ) mod 2 1 Cer , RSSUDRE ,
rrc
current
a I fc Brrc
4 t t t of , ( IODF current
IODFprevious ) mod 2 1
Ciono _ ramp , Ciono _ step ,
Iiono , RSSiono
Master of Science in GNSS @ J. Sanz & J.M. Juan
80
gAGE
Degradation of Ionospheric Corrections
1
R cos E 2
s F s
research group of Astronomy and Geomatics
2
2 2 2 Fpp 1 e
UIRE pp UIVE Re hI
iono
tiono,GIVEi
Protection
Level
True
Error
N
s i2 s i2, flt s i2,UIRE s i2,air s i2,tropo
sV2i s i2
gAGE/UPC
VPL KV
i 1
GEOMETRY
EGNOS
Availability
Integrity
gAGE/UPC
Accuracy
EGNOS
Availability
Integrity
Continuity
gAGE/UPC
Accuracy