Geodesy Lec 1
Geodesy Lec 1
Essam Ghanem
0/27
Webster: a branch of applied
mathematics concerned with the
determination of the size and shape of the
earth and the exact positions of points on
its surface and with the description of
variations of its gravity field
1/27
Vertical Datums
The Geoid
Gravity: Local Attraction
Unfortunately, the density of the earth‟s crust is not uniformly
the same. Heavy rock, such as an iron ore deposit, will have a
stronger attraction than lighter materials. Therefore, the geoid
(or any equipotential surface) will not be a simple mathematical
surface.
2/27
Vertical Datums
The Geoid
What is the GEOID?
• “The equipotential surface of the
Earth‟s gravity field which best fits, in
the least squares sense, global mean
sea level.”
• Can‟t see the surface or measure it
directly.
• Modeled from gravity data.
3/27
Vertical Datums
The Geoid
Equipotential Surfaces
Topography
4/27
An ellipsoid of So
revolution
we squash is the
the figure
spherewhich
to would be
obtained by rotating an ellipse
fit better at the about
poles. its shorter axis.
The GRS80 ellipsoid is used for the Globe.
Close Fit At The Equator This creates a spheroid
b = 6,356,752.31414 m
But The Poles Are Out
a = 6,378,137.00000 m
7/27
Geodetic Datum
• Geodetic datum – The
information that ties a geoid
model to an ellipsoid model,
and to a real point on (or in)
the Earth by a coordinate
system
What are Geodetic Datum?
`Geodetic Datum„ is the process of fixing an ellipsoid
of a particular orientation and position with respect
to the Earth.
To fix the ellipsoid to the mean earth, basically, we
need to know the relationship between the chosen
ellipsoid and the earth in terms of position and
orientation. It is defined by seven parameters.
GRS80-WGS84 HELMERT1906
GEOID
9/27
What are Geodetic Datum?
The 7 parameters that define the Geodetic Datum:
• Two parameters for the size and shape of the ellipsoid
• Three parameters for the distances of the ellipsoid
center from the center of the earth
• One parameters for the directions of the rotation axis
with the mean rotation of the earth
• One parameters for the direction of 0 longitude with
the earth‟s international 0 longitude
10/27
Geodetic Datum
Two types: horizontal and vertical
datum
horizontaldatum - forms the basis
for the computations of horizontal
control surveys in which the
curvature of the earth is considered
vertical
datum – elevations are
referred
11/27
Local Horizontal Datum
Components
• Parameters of the
ellipsoid
– axis length
– flattening value
• Parameters that tie
the ellipsoid to the
origin point (known
place on the Earth)
12/27
The Ellipsoid
The 7 parameters that define the Geodetic Datum:
• The ellipsoid is an ellipse of revolution rotating around
its minor axis. Its size and shape are chosen to best fit
the geoid in the area of interest thus allowing as many
ellipsoids as you care to define.
• It is this figure however that terrestrial observations
must be reduced and to which geodetic coordinates
refer. Z
• The ellipsoid is an ellipse of
revolution rotating around its
minor axis. Its size and shape b
are chosen to best fit the geoid
in the area of interest thus a Y
allowing as many ellipsoids as a O
you care to define.
X
Rotational axis
13/27
The Geometry of Ellipsoid
Reference ellipsoids are identified by a name and often by a year
like Helmert1906 and defined by either:
•semi-major and semi-minor axes, or the relationship
between the semi-major axis and the flattening of
the ellipsoid (expressed as its eccentricity).
N. America Europe
O1
O2
S. America Africa
15/27
Basic Geodesy
Local vs. Global Reference Ellipsoid
HELMERT1906
GRS80-WGS84
GEOID
16/27
Basic Geodesy
EGYPT
ELLIPSOID DEFINITIONS
HELMERT1906 :
a = 6,378,200 m 1/f = 298.3
17/27
The World
18/27
The Greenwich Meridian
19/27
The Greenwich Meridian
Meridians are great circles that pass through the Polar
axis of the Earth. In 1884, the International Meridian
Conference established the meridian passing through the
transit instrument at the Royal Observatory in
Greenwich, England as the prime meridian. The prime
meridian is also known as the Greenwich Meridian, and
has a longitude of 0°.
• Meridians to the west of
Greenwich are assigned
values of 0° to 180° West
longitude. Meridians to
the east of Greenwich are
assigned East longitudes.
• positive in the eastern
hemisphere
• negative (-) in the western
hemisphere
• longitudes are
traditionally designated
by the Greek letter λ. 20/27
The Equator
21/27
The Equator
22/27
The Equator & Parallel of Latitude
Parallels of Latitude are approximately same distance apart
and run in a north-south direction. The are defined as the
arc distance in the plane of the meridian from the
equatorial plane to the prime vertical.
Due to earth flattening, 1° of latitude is approx. 69.4 mi
near poles and 68.7 mi. near
equator.
Latitude are considered
• positive in the northern
hemisphere
• negative (-) in the
southern hemisphere
• latitudes are traditionally
designated by the Greek
letter φ.
23/27
The Centre of the World
24/27
P
25/27
The Coordinate Systems
The coordinate system, in general, is a reference
frame in which any point can be defined either in
two or three dimensions.
The coordinates of any point are necessary to be
known in order to:
• Define the horizontal and/or the vertical position.
• Deal with different points on the earth’s surface.
• Obtain any geometric relation between different points.
• Plot easily and accurately any point on the map plane
26/27
The 3D Geodetic Curvilinear Coordinate
Geodetic latitude : is the angle measured form the equator to the point on the
meridian plane of that point. varies from zero to 90 north or south.
Geodetic longitude : is the angle measured from Greenwich meridian to the
point on the geodetic equator anti-clockwise direction. varies from 0 to 360
east or form 0º to 180º east and form 0º to 180º west.
Geodetic height h: is the distance between the point and the sphere measured on
the ellipsoid normal direction.
27/27
Complications
• Ellipsoid does not possess a constant radius of
curvature.
• Radii of curvature are dependent on latitude only.
• Two radii are of interest
– Meridian
– Prime Vertical
28/27
Complications
The radius used for the longitude is called the
"Radius of Curvature in the prime vertical".
The radius used for the latitude change to North
distance is called the "RM Radius of Curvature in
the meridian.
“RN is the radius of a
circle that is tangent to
the ellipsoid at the
latitude and has the
same curvature as the
ellipsoid in the north-
south direction there.
29/27
The Radii of Ellipsoid
Merdians : Meridians are great circles on the circumference
of the ellipsoid which pass through the north and south
poles. Any plane containing a meridian and the polar axis is
a meridian plane.
Radii of Ellipsoid : Because
the Earth is approximated by
an ellipsoid and not a sphere,
the great circle that defines
the prime vertical at P has a
radius RN that is different
than the radius RM in the
meridian at P. The lengths of
these two radii, which are
collinear at any point, are
used in many geodetic
computations. 30/27
The Radii of Ellipsoid
31/27
Example:
Using the GRS80 ellipsoidal parameters, what are the radii for the
meridian and prime vertical for a point of latitude 41°18’15.0132” N
What is the radius of the great circle that is at an azimuth of 142°14’36”
at this point?
Solution
For the GRS80 ellipsoid, yields:
32/27
Meridian Convergence - ()
Angular convergence dependent on the distance between the meridians,
the latitude and the dimensions of the earth
33/27
Different Types of Latitudes
34/27
Geocentric Latitude
35/27
Reduced Latitude
36/27
The 3D Cartesian Coordinate
ECEF X, Y, and Z Cartesian coordinates define three dimensional
positions with respect to the center of mass of the reference ellipsoid.
The Z-axis points from the center toward the North Pole.
The X-axis is the line at the intersection of the plane defined by the
prime meridian and the equatorial plane.
The Y-axis is defined by the intersection of a plane rotated 90o east
of the prime meridian and the equatorial plane.
.
37/27
The Geodetic to Cartesian Coordinates
It should be mentioned here; that any of the above two sets
of coordinates can be mathematically transformed to the
other set and vice versa.
Z
h
g
Where: N is the Prime Vertical Y
radius and is equal to: j
X
Y
X
38/27
The Cartesian to Geodetic Coordinates
39/27
Greenwich Z Meridian of longitude
N
meridian Parallel of latitude
=0°
P - Geographic longitude
•
j - Geographic latitude
R - Mean earth radius
W O j •
E
Y
R
•
Equator =0°
•
X O - Geocenter
40/27
Length on Meridians and Parallels
Length on a Meridian:
AB = Re R
R D
(same for all latitudes)
C
Re B
Length on a Parallel: Re
A
CD = R = Re Cos
(varies with latitude)
41/27
Example:
What is the length of a 1º increment along
on a meridian and on a parallel at 30N, 90W?
Radius of the earth = 6370 km.
Solution:
A 1º angle has first to be converted to radians
radians = 180 º,
so 1º = /180 = 3.1416/180 = 0.0175 radians
43/27
Calculations on Sphere
Calculations on Sphere
Calculations on Sphere
Napier's rules
The Geodetic Height
P
Ellipsoid Heights h
Heights determined
using GPS are
Ellipsoid
referenced to the
WGS 84 Ellipsoid
49/27
Vertical Datums
Ellipsoid vs. Geoid
• Ellipsoid
– Simple Mathematical Definition
– Described by Two Parameters
– Cannot Be 'Sensed' by
Instruments
• Geoid
– Complicated Physical
Definition
– Described by Infinite
Number of Parameters
– Can Be 'Sensed' by Instruments
50/27
Vertical Datums
Ellipsoid vs. Geoid
High Density
51/27
52/27
Geoid Model
53/27
The Relation between the geoid
and the ellipsoid
h is the geodetic height, which is
defined as the distance from the
point to the ellipsoid surface
measured on the perpendicular line
to the ellipsoid at that point
H is the orthometric height, which is
defined as the distance from the
point to the geoid surface measured
on the perpendicular line to the geoid
at that point.
Orthometric heights H are ref. to a
Datum which is M.S.L.
M.S.L approximates the Geoid
N is the geodetic undulation, which is
the difference between the geodetic
height and the orthometric height.
54/27
The Relation between the geoid
and the ellipsoid
55/27
Vertical Datums
h =H+N
H is measured traditionally
h is measured with GPS Observations
N is modeled using Gravity Models
57/27
58/27