3/7/2012
Derivation of the WGS84-PRS92
Transformation Equation:
A review
Focused Group Discussion
Engr .John Louie D. Fabila
Coordinate Transformation
Bursa-Wolfe
2D
Transformation
3D Conformal
Coordinate Transformation
Mathematical process of converting
coordinate values from one coordinate
system to another
Used in various aspects of geodetic
practice, such as geometric and satellite
geodesy, geodetic astronomy,
photogrammetry and remote sensing
3D Conformal Coordinate Transformation
3D
Transformation
3D Affine
Molodensky
, , , s , X , Y , Z
3-Parameter
where
, , = rotation in x,y,z
Transformation
equations
can be expressed in 7
independent parameters;
x
y
Y
Y
s = scale
3D Polynomial
X , Y , Z = translations in x,y,z
3/7/2012
ROTATION MATRICES
To define the angular orientation of a 3D
coordinate system with respect to another
3D coordinate system, we need to define 3
rotation parameters
ROTATION IN A PLANE
ROTATION IN A PLANE
Given point P,
express its new
XY coordinates
as a function of its
old XY
coordinates, and
the rotation angle
3/7/2012
ROTATION IN A PLANE
Rotation matrices are orthogonal
r112 + r212 = 1
r11r12 + r21r22 = 0
ROTATION IN SPACE
,, are the rotation angle we need to
apply to make the 2 coordinate systems
parallel
Rotation matrix has nine elements,
governed by these 3 parameters
Construct xyz coordinate system parallel
to XYZ whose origin coincides with that of
xyz coordinate system
R 1 = RT
ROTATION IN SPACE
0
0 x '
x 1
y = 0 cos sin y '
z 0 sin cos z '
ROTATION IN SPACE
x
x '
y = M y '
z
z '
3/7/2012
ROTATION IN SPACE
x cos
y = 0
z sin
x
y = M
z
0 sin x
1
0 y
0 cos z
ROTATION IN SPACE
x cos
y = sin
z 0
x
y
z
y = M
z
FROM coord. system
x ' x
y ' = y = M M M
z ' z
m11
M = M M M = m21
m31
cos
0
0 x
0 y
1 z
x
y
z
ROTATION IN SPACE
ROTATION IN SPACE
Coord. system parallel to TO coord.
system
sin
m11 = cos cos
m12 = sin sin cos + cos sin
m13 = cos sin cos + sin sin
x
y
z
m12
m22
m32
m21 = cos sin m22 = sin sin sin + cos cos
m23 = cos sin sin + sin cos
m31 = sin
m32 = sin cos
m33 = cos cos
m13
m23
m33
cos cos
M = cos sin
sin
sin sin cos + cos sin
sin sin sin + cos cos
sin cos
cos sin cos + sin sin
cos sin sin + sin cos
cos cos
3/7/2012
3D Coordinate Transformation
3D Coordinate Transformation
After rotating, scale and translate
X
x ' X
Y = s y ' + Y
Z
z ' Z
X
x X
Y = sM y + Y
Z
z Z
X = sMx + T
Special Case: Datum
Transformations
Almost all modern datums are nearly
aligned, leading to small rotation angles
Datums using the same units of measure
have scale factor nearly equal to 1
The difference in size of ellipsoids used
are in the order of 1x10-6
Assumptions used
For small angles (in radians),
sin =
cos = 1
Product of 2 sines are equal to zero
Since scale is almost equal to 1, use (1+s)
instead of s for scale, where s is the
deviation from a value of 1.
For small values of s, we can use binomial
series expansion:
1
= 1 x + (higher order)
1+ x
3/7/2012
Bursa-Wolfe Transformation
Equation
Coordinate Transformation
cos cos
M = cos sin
sin
sin sin cos + cos sin
sin sin sin + cos cos
sin cos
1
M = RZ
RY
RZ
1
RX
X
Y = sM
Z
cos sin cos + sin sin
cos sin sin + sin cos
cos cos
RY
RX
1
X
1
Y
6
=
1
+
s
10
(
)
RZ
Z PRS 92
RY
Reverse Equation
RZ RY X X
x
1
1
y =
1
RX Y + Y
1 + s 106 RZ
z
RY RX
1 Z Z
RZ
RY X
X
1
X
Y
6
= (1 s 10 ) RZ
1
RX Y
+ Y
Z WGS 84
RY RX
1 Z PRS 92 Z
x X
y + Y
z Z
RZ
1
RX
RY X
X
RX Y
+ Y
1 Z WGS 84 Z
NRMDP results(Jones, 1991)
Rx
Ry
Rz
3.06762"
-4.90291"
-1.57790"
X
Y
Z
s
127.62195
67.24478
47.04305
1.06002
29 stations common
to both WGS84 data
set and the Luzon
Datum data set were
used
3/7/2012
NAMRIA WEBSITE