01 July 2020 - Introduction To Global Positioning System by Dr. Ashutosh
01 July 2020 - Introduction To Global Positioning System by Dr. Ashutosh
Ashutosh Bhardwaj
He is Scientist/Engineer-‘SF’ at Photogrammetry and Remote
Sensing Department, Indian Institute of Remote Sensing,
Dehradun. He is a graduate in Civil Engg. from M.R.E.C. (present
MNIT), Jaipur; M.Tech. (Remote Sensing) from BIT, Ranchi & PhD
from Civil Engg. Dept., IIT-Roorkee. He joined ISRO in 2001 after
a brief career as faculty at Dept. of Civil Engg., BITS, Pilani.
• He has been engaged in industry, teaching and research in Surveying,
Photogrammetry, GNSS, Cartography, Remote Sensing and topographic
modelling for the past 20 years. He has published over 50 research papers
in Journals and conferences. He has widely contributed to various
departmental projects on remote sensing and Mapping at NRSC and IIRS.
• He has received appreciation from the Ministry of Planning and National
Development (MPND), Republic of Maldives for rendering expert services
for GPS ground control survey in the National Mapping Project (Maldives).
ORCID iD: https://orcid.org/0000-0002-9241-5427
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Satellite Navigation
Global Navigation Satellite System
&
GPS receivers, surveying
techniques, processing methods,
errors and accuracy
Dr. Ashutosh Bhardwaj
Scientist/Engineer-’SF’, PRSD (GT&OP Group)
[email protected]
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Xiaomi Mi8
Source: https://rntfnd.org/wp-content/uploads/eLoran-Definition-Document-0-1-Released.pdf
Source: https://www.geospatialworld.net/blogs/the-redmi-note-9-pro-with-navic-support-is-
here/#:~:text=The%20Redmi%20Note%209%20Pro%20with%20a%20Qualcomm%C2%AE%20Snapdragon,175)%20on%20the%20Mi%20Store.&text=The%20Navigation%20and%20positioning%20features,%2F%20Galileo%2F%20GLONASS%20%2F%20Beidou.
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
SATELLITE NAVIGATION
• A satellite navigation system is a system that uses satellites to
provide autonomous geo-spatial positioning. Example:
• GLOBAL
• NAVSTAR GPS
• GLONASS
• BEIDOU
Transition: As a direct result of the tragedy killing all 269
• GALILEO aboard Korean Air Lines Flight 007 which strayed into
Russian airspace accidentally & shot down by Soviet air-
to-air missiles on September 1, 1983, near Sakhalin
• REGIONAL Island, Russia, President Ronald Reagan announced on
September 16, 1983 that the GPS system that had
• IRNSS previously been intended for U.S. military use only
• QZSS would now be made available for everyone to use.
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Position
• 1010-20m
- 20 m 10-20m ~10-20m ~10-20m
• Velocity ~ 10cm/s
Velocity 10cm/s 10cm/s currently)
20cm/s 20/40cm/s
• Time ~ 0.1 µs (15ns
• 6/3 0.1
Time orbits
µs (15ns) 20ns 20ns 20/50ns
1 •
st Worldwide
Satellite PRN 4 Coverage
Kosmos1413 *GIOVE-A Beidou-1A
• 24 hour access 12Oct.1982 28Dec,2005 30Oct.2000
22Feb.1978
Ellipsoid WGS84 PZ - 90.11 GTRF CGCS2000* *
• Common Coordinate System
# ##
Space Segment
(Initial Operational Capability(IOC)-1993)
(Full Operational Capability(FOC)-1995)
Block I Block II/IIA
Two Launches: Dec. 23, 2018 and 22 August, 2019
Total IIIA Series: 10 planned; IIIF Series: 11th onwards
Entirely new design
new ground control system (known as OCX $500 million/ $6
billion : FOC by 2021 by Raytheon). OCX passed
cybersecurity tests.
first GPS sent aboard on a SpaceX rocket
M-Code: more powerful GPS 3 signal for military users &
First Launch: 22 Feb 78(78-85)# First Launch: 14 Apr 89(89-97) more secure against jamming or spoofing.
On-Orbit: None, Total=11 Total II Series: 27 (1+14+12+8) S-Band serial telemetry link
Block IIIA/IIIF
GLONASS satellites
Source: https://www.glonass-iac.ru/en/guide/
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
55̊
Equator
Colorado
Springs
Ascension Kwajalein
Hawaii
Islands
Diego
Garcia
GLONASS signals
• The GLONASS satellite signal identifies the satellite and provides:
position, velocity and acceleration vectors at a reference epoch to compute
satellite locations
synchronization bits, data age and satellite health
offset of GLONASS time from UTC (SU) (formerly Soviet Union and now Russia)
almanacs of all other GLONASS satellites
Unlike GPS, all GLONASS satellites transmit the same code at different
frequencies. They derive signal timing and frequencies from one of three on-board
cesium atomic clocks operating at 5 MHz:
For example, L1 = 1602 MHz + (n x 0.5625) MHz where n = the frequency
channel number (n = 0, 1, 2 and so on)
L2= 1246 MHz + n × 0.4375 MHz
• The frequency ratio f2 /f1 is constant for all GLONASS satellites and amounts to
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
The GNSS constellation system’s potential civil applications are many and mirror those of GPS.
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Why GPS?
SA: Selective
availability
Xll
Acc.: 10-100m
Vl
Multi-Satellite Ranging
Timing
• Accuracy of position is only as good as your clock
• To know where you are, you must know when you receive.
• Receiver clock must match SV clock to compute delta-T
• Each Galileo satellite has two master passive hydrogen maser atomic clocks and
two secondary rubidium atomic clocks which are independent of one other
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Interactive Session
• Which is the counterpart system of LORAN: • Which is the counterpart system of Alpha:
• Chyka • Chyka
• GLONASS • Omega
• GPS • GPS
• Alpha • GLONASS
• Galileo is from:
• Beidou is from:
• Russia
• Russia
• European Union
• European Union
• USA
• USA
• Japan
• China
Position Equations
P1 ( X X 1 ) 2 (Y Y1 ) 2 ( Z Z 1 ) 2 b
P2 ( X X 2 ) 2 (Y Y2 ) 2 ( Z Z 2 ) 2 b
P3 ( X X 3 ) 2 (Y Y3 ) 2 ( Z Z 3 ) 2 b
P4 ( X X 4 ) 2 (Y Y4 ) 2 ( Z Z 4 ) 2 b
Where: Pi = Measured PseudoRange (Biased ranges) to the ith SV
Xi , Yi , Zi = Position of the ith SV, Cartesian Coordinates
X , Y , Z = User position, Cartesian Coordinates, to be solved-for
b = User clock bias (in distance units), to be solved-for
• The above nonlinear equations are solved iteratively using an initial estimate of the user
position, XYZ, and b- same for all satellites.
• To solve the user position equations, one must know where the SV is:
• The navigation and time code provides this
• 50 Hz signal modulated on L1 and L2
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
GPS Time
• GPS time is referenced to 6 January 1980, 00:00:00
• Jan 6 = First Sunday in 1980
• GPS satellite clocks are essentially synched to International Atomic Time (TAI)
(and therefore to UTC/zulu time since Jan. 1, 1972 for global civil time)
• TAI, maintained at Lab., France, is the basis for Coordinated Universal Time (UTC), used for
most civil timekeeping
• GPS time = TAI + 13s
• Since 13 leapseconds existed on 1/6/1980
L3(1381.05MHz); L4 (1379.913MHz): used only for a atomic flash detection; Nudet (Nuclear Detection) System
(NDS). L5: 1176.45 MHz (25.5 cm, In-door apps., anti-jamming). L2C: 1227.60MHz (in pre-operational testing
and available on 24 satellites since May, 2017). L1C:1575.42MHz (III series onwards)
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Navigation Message
n In order to solve the user position equations, one must know where the SV is:
n The navigation and time code provides this
u 50 Hz signal modulated on L1 and L2
n The SV’s own position information is transmitted in a 1500-bit data frame
u Pseudo-Keplerian orbital elements
F Determined by control center via ground tracking
u Receiver implements orbit-to-position algorithm
n Also includes clock data and satellite status
n And ionospheric / tropospheric corrections
n International Telecommunication Union (ITU) has reserved 1559-1610MHz band for satellite
based navigation through World Radio Communication (WRC) conferences, held every three
year.
n GPS bands (US Federal Communication Commission): (1215-1240MHz, 1559-1610 MHz, L5-
1164-1188MHz)
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Navigation Message
The Almanac
• In addition to its own nav data, each
SV also broadcasts info about ALL the
other SV’s
• In a reduced-accuracy format
• Known as the Almanac
• Permits receiver to predict, from a
cold start, “where to look” for SV’s
when powered up
• GPS orbits are so predictable, an
almanac may be valid for months
• Almanac data is large
• 12.5 minutes to transfer in entirety
Source:www.glonass-iac.ru/en/GLONASS/ephemeris.php
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Phase Observations
Received Satellite
Wavelength of the signal is 19 cm on L1 and 24 cm
on L2 Phase
Receiver compares self-generated phase with
received phase Generated
Number of wavelengths is not known at the time Phase from
the receiver is switched on (carrier phase ambiguity) Receiver
As long as you track the satellite, the change in
dT
distance can be observed (the carrier phase
ambiguity remains constant)
D = c .dT + N λ
If longer PRN code is used, receiver becomes more resistive to Jamming signal. But, signal
processing is more complex
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Satellite Geometry
Satellite geometry can affect the quality of signals and accuracy of receiver
trilateration.
Positional Dilution of Precision (PDOP) reflects each satellite’s position relative to
the other satellites being accessed by a receiver.
PDOP can be used as an indicator of the quality of a receiver’s triangulated
position.
It’s usually up to the GPS receiver to pick satellites which provide the best
position trilateration.
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Satellite Geometries
Ideal Satellite Geometry Good Satellite Geometry
GNSS
Systems
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Advantages of GNSS
• The use of GLONASS in addition to GPS provides very significant
advantages:
increased availability of satellites & signals.
markedly increased spatial distribution of visible satellites
reduced HDOP and VDOP (DOP) factors
Better atmospheric correction
decreased occupation times means faster RTK results
A larger satellite constellation also improves real-time carrier phase differential
positioning performance.
accurate, robust & reliable services even in bad conditions
Less expensive high-end services
Source: https://www.novatel.com/assets/Documents/Papers/GLONASSOverview.pdf
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
STEPS IN PHOTOGRAMMETRY
Aerial Photo
/ satellite image Input
Aero/Sat. Triangulation/
Control extension and Block
Adjustment
Digital Topographic
maps
Ortho Image
/ Mosaic Generation
Assessment of the MERIT DEM and EarthEnv DEM90 DEMs using GCPs
Table 1: Vertical accuracy measures computed for the experimental sites
Experimetnal ME MAE RMSE St. Dev. LE90
Sites (m) (m) (m) (m) (m)
Jaipur Site 1.37 4.03 3.27 2.97 4.91
Kendrapara -3.33 3.33 4.00 2.21 3.65
site
Dehradun Site 3.17 4.79 7.82 7.15 11.79
Chandigarh 0.27 4.79 5.59 5.58 9.20
Site
Kalka Site -0.19 12.82 16.61 16.58 27.36
Delhi Site -11.19 10.07 13.14 6.90 11.38
FUSED DEMs
Elevation extraction at
GCP locations
Datum Conversion
Figure 1: Location map of the MERIT experimental sites with DEMs
GCPs
Table 2: Comparison of vertical accuracy for the MERIT DEM & EarthEnv DEM90
Experimetnal Sites ME (m) MAE (m) RMSE (m) Statistics generation
DEMs-90m MERIT EarthEnv MERIT EarthEnv MERIT EarthEnv
Jaipur Site 1.37 0.70 4.03 2.28 3.27 3.05
Kendrapara site -3.33 -3.64 3.33 3.64 4.00 4.22 Quality Assessment
Dehradun Site 3.17 0.57 4.79 5.42 7.82 6.55
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
https://www.gim-international.com/files/0812346fc5116deee99fe6f235d260dc.doc
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Services
NAVSTAR GPS:
• Standard Positioning Service(SPS)
• Precise Positioning Service(PPS)
GLONASS:
• Standard Precision Service(SPS)
• High Precision Service(HPS)
Galileo:
• The Open source (OS) ECEF
• The Safety of life (SoL)
• The Commercial Service (CS) by Galileo Operating Concessionaire (GOC)
• The Public Regulated services(PRS): European Police/Antifraud offices.
• The Search & Rescue(SAR-L6) Service:+ to COSPAS-SARSAT system
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Receivers
Single Frequency Dual Frequency
• Baseline Accuracy 1cm/5mm + 2/1ppm (rms) • The “high end” of the GPS Market
• Uses Post process L1 carrier phase • Baseline Accuracy ranging from 5/3/3.5mm +
1/0.4ppm (rms)
• Used for all Surveying tasks with baselines up
to 15Km • Used in all GPS Surveying tasks :-
• Network Densification, Detail Surveys • Geodetic Control Networks, Tectonic Plate
• Real Time Monitoring, Photogrammetric Control,
• Smaller Occupation time Network Densification, Detail Surveys, etc.
• Less expensive alternative to Dual frequency • Real Time
• Most unsophisticated receivers track only L1 • Smaller Occupation time
and use a simplified correction model • New applications are found on a daily basis
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
• Mapping grade
• capable of <3 meters accuracy
• portable, less expensive
• accurate mapping for integration with GIS
• Navigation
• capable of 10 meters accuracy
• light weight, cheap
• basic navigation, limited data storage
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Sources of error
1. [Selective availability]
2. Clock errors
3. Ephemeris errors
4. Atmospheric delays
5. Multipath effects
6. Receiver errors
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Error Budget
User Range Errors (URE) consists of System Errors (Ephemeris Data, & Satellite
Clocks) and atmospheric Errors (Ionospheric/Tropospheric: 90% Dry delay due
to dry gases; 10% Wet delay due to WV & condensed clouds).
User Equipment Errors (UEE) consist of Receiver Noise and Multipath error.
PDOP of 2 means that in the worst case, a 1m URE will result in 2m positional
error.
Typical Error in Meters (per satellite)
Standard GPS Differential GPS
Satellite Clocks 1.5 0
Orbit Errors 2.5 0
Ionosphere 5 0.4
Troposphere 0.5 0.2
Receiver Noise 0.3 0.3
Multipath 0.6 0.6
SA 30 0
30m
• Accomplished by:
• “Dithering” the clock data
• Results in erroneous pseudoranges
• Truncating the nav message data
P
• Erroneous SV positions used to compute user pos.
• Degrades SPS solution by a factor of 4 or more
• Long-term averaging is only effective SA compensator
• Positional accuracy 50m(1sigma), 100m (95%) +/- 100m (95%)
Atmospheric Delay
GPS signals are delayed as they
pass through the atmosphere (a.
layer of charged ions and free
electrons known as ionosphere and
b. the troposphere)
< 10 km > 10 km
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Tropospheric Refraction
• Troposphere is the lowest layer of atmosphere varying from (ground to) 7 to
14kms.
• Pseudo range errors varies from 2m if the satellite is at zenith to 30m for a
5° elevation satellite.
• Hydrostatic delay occurs due to presence of dry gases: about 2 meters in the
zenith direction to 10 meters for lower elevations.
• Wet delay occurs due to water vapor and condensed clouds: only some tens
of centimetres.
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Atmospheric Refraction
Ionsospheric Refraction Depends on:
•Sun’s activity (11 years sunspot cycle),
•Ionospheric thickness & proportions/concentration of ionized particles,
•Season,
•Actual path (i.e. relative position of satellite & receiver) of the wave i.e. signal
•Pseudo range errors vary from 0 to 15m at zenithal incidence to as much as 45m
for low incidence signals.
•delay is proportional to the integral no. of free electrons along the transmission
path and inversely proportional to the square of the transmission frequency.
Multipath Error
SV1
Choke Ring
antenna
SV2
Caused by local reflections of the GPS/GNSS signals that mix with the desired signals
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Signal Obstruction
SV1 SV2
SV3
When something blocks the GPS signal.
Areas of Great Elevation Differences
Canyons
Mountain Obstruction
Urban Environments
Indoors
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Human Error
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Quiz Time
• DGNSS requires minimum…. GNSS receivers: • Beidou use following reference coordinate system:
• 2 • CGCS
• 3 • GTRF
• 4 • ITRF
• 5 • WGS 84
Use
Differential GPS/GNSS
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
High Accuracy
(Receiver position, satellite position, frequency-ionospheric
corrections, time-ambiguity of carrier phase measurements)
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Differential Positioning
• It is possible to determine the position of Rover ‘B’ in relation to
Reference ‘A’ provided
– The coordinates of the Reference
Station (A) are known
DGPS Method
Method to remove errors from GPS measurements
Uses a GPS receiver at a fixed, surveyed location i.e. base station, to measure
error in pseudo range
Base station receives the same GPS signals as the roving receiver and instead
of using timing signal to calculate its position, it uses its known position to
calculate timing.
Pseudorange error for each satellite is subtracted from rover before calculating
position during RTK mode or post processing
Post - Processed
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Real - Time
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
R1 R2
R8
R3
R7
Base
R4
R6
R5
Case1: Survey with Two GNSS Geodetic Receivers
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
R1/3
R1/7
Base
R2/4
R2/6
R1/5
Atmospheric partial
Multipath no
Receiver noise no
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
Point
Line Area
Calendar Contacts
Today Messaging
screen
Directional
MIO Power/USB socket button
INDIAN INSTITUTE OF REMOTE SENSING, DEHRADUN
References
http://www.glonass-ianc.rsa.ru/pls/htmldb/f?p=202:1:15000421459964108253
http://igscb.jpl.nasa.gov/
http://www.navcen.uscg.gov/gps/precise/default.htm
https://rntfnd.org/wp-content/uploads/eLoran-Definition-Document-0-1-Released.pdf
Interface Control Documents:
http://www.navcen.uscg.gov
http://www.Glonass-ianc.ras.ru
http://www.Galileoju.com
Links:
UNAVCO http://archive.unavco.ucar.edu/cgi-bin/dmg/pss
CDDIS http://cddisa.gsfc.nasa.gov/cddis.html
NGS/CORS http://www.ngs.noaa.gov/CORS/ Contact Details of the Faculty:
SOPAC http://sopac.ucsd.edu/
Email- [email protected]
Tel- 0135-2524117