Reference systems, datums and coordinate systems

TOWARDS IMPLEMENTING UNIVERSAL TRANSVERSE MERCATOR
COORDINATES IN GHANA.
John Ayer
Department of Geomatic Engineering
KNUST
Kumasi - Ghana
Abstract.
Over the years, different Nations have defined and used different Spatial Reference
Systems for mapping and spatial data gathering and archiving. In recent times, the
need to share information and also map and manage transboundary resources on
Regional, Continental and Global basis is leading to the adoption of common
reference frame initiatives such as the African Reference Frame (AFREF).
The Universal Transverse Mercator (UTM) coordinate system has been proposed for
some years now for global mapping but its adoption has been slow because it
required new mapping to replace existing maps and map data held by different
Nations over the years. However, its adoption has advantages for data integration.
Besides, recent technological advancements such as the development of Global
Navigation Satellite System (GNSS) which computes coordinates based on WGS84
ellipsoid directly and also provides a program for projecting results in UTM has
provided a further impetus for adopting UTM coordinates.
This paper discusses the different scenarios that could result from the adoption of
UTM and proposes what steps are required for its implementation in Ghana.
Keywords.
Spatial Reference Systems, Universal Transverse Mercator, GNSS, Horizontal
Datums, LWOST.
1
Introduction.
1.1 Reference systems, Datums and Coordinate systems
Spatial reference systems help to integrate data from different sources together such
as images, maps and ground surveys. They also allow one to refer to positions on and
near the Earth’s surface uniquely and unambiguously by allowing the assignment of
unique values to each position in space.
Examples of reference systems include the International Terrestrial Reference System
(ITRS), the Ghana Coordinate System, the Geodetic Reference System of latitude and
longitude definitions using various ellipsoids, Metes and Bounds, and Postal codes
etc.
Datums refer to a set of parameters (values) and control points used to define the
shape of the earth, and the origin and orientation of the coordinate system used by a
country. Ordinarily, the word means facts which are known from which other
conclusions can be drawn. When used with reference systems, it means origin points
with assigned values which could then be used for deriving values of other points.
One may therefore define Horizontal datums for horizontal measurements, Vertical
Datums for height measurements and combined Datums that may combine both
horizontal and vertical measurements.
To understand the concept of datum properly, it must be noted that Location will only
make sense when it is compared to some frame of reference. For example one may
speak of a Building with house number H215 as being opposite the General Post
Office in Accra. This description is relative but in absolute terms, the same building
could be stated as having positional coordinates “05° 35’ 06” N, 00 ° 01’ 09” W” (of
course value depends on ellipsoid!!), or “333382.141 N, 1256469.743 E in Ghana
National Datum (GND)” or “X=6348703.172, Y= -2124.826, Z=616629.997
Cartesian values” or else “618475.6234 N, 830306.7194 E, in UTM Z 30 N”. Thus
the same building has different positional coordinate values depending on the chosen
reference frame.
2
1.2 Horizontal Datum
A Horizontal Datum comprises a base reference for horizontal coordinates.
It
includes the position of an initial point of origin and an ellipsoid that models the
surface of the earth in the region of interest. The complete Datum definition consists
of the origin point of the X,Y,Z, axis (three parameters), the orientation of the axes
(three parameters), the ellipsoid used in terms of its shape and size (two parameters)
[Steed, J., 1990].
Why an ellipsoid?
The picture below (figure 1) is that of earth from ERS-1
Figure 1
It could be seen that:
•
The shape is complex and not suitable for any mathematical calculations.
•
The simplified mathematical surface that can approximate to portions of its
surface is the ellipsoid.
•
No single ellipsoid will accurately fit all parts so each region may have
separate best fitting ellipsoid.
•
A global average ellipsoid may be defined.
The reason an ellipsoid is used and not a perfect sphere is because the Earth spins
making it bulge out a little bit at the equator. Since the actual Earth is irregularly
shaped, one can come up with a number of different values for the major radius (a),
3
the minor radius (b), and the flattening (f), depending on where they are measured.
[Snyder, J.P.]
1.3 Traditional Horizontal Datums
Many Nations established their own national datums using various national standards
and procedures by calculating ellipsoids that fit well locally. They also established
initial point location and orientations usually with astronomic observations. The
results of the different datums are that, transformations are necessary to convert
between them because ellipsoid shape and size differences result in origin shifts and
rotations. [EPSG, 2005].
To illustrate some limitations to the traditional approach, consider the difference in a
position as defined on the global average datum WGS84, (ellipsoid values Semimajor Axis: a = 6378137 m, Semi-minor Axis: b = 6356752.3142m, Flattening Ratio:
f=(a-b)/a = 1/298.257223563) and the Ghana Horizontal Datum ( a= 6378299.996 m.
b=6356751.69m. and f=1/296)
Positional coordinate values on WGS84, ( N 06° 04' 49.84347" E 0° 03 ' 0.85954" )
Positional coordinate values on Ghana Datum (N 06 ° 04 ' 39.9216 "
E 0 ° 02 '
59.7823 ")
Separation difference between the points on the different datum definitions is 306.6m.
which is more than 1000 feet.
1.4 Vertical Datums
Like horizontal measurements, elevation only has meaning when it is referenced to
some start point. For this reason, different levels of starting references have variously
been used such as the low tide level, the high tide level and the mean sea level which
coincides with a local geoid definition. Fig 2. The Global Positioning System (GPS)
however uses a height above the global ellipsoid.
4
Figure 2: Height references
1.5 Coordinate systems
A Coordinate system is a spatial reference system that allows one to assign unique
values (coordinates) to each location and establish relationships between sets of such
coordinates. It enables the interpretation of a set of coordinates as a representation of
a unique position in a real world space. [Ayer, J. and C. Fosu, 2008]
In Ghana, we have been associated with five (5) different types of coordinates from
Colonial days until now namely:
i. Geodetic or Geographic Coordinate System
This is a Global Coordinate System based on angles defined on the ellipsoidal
Earth surface. It is expressed in longitudes and latitudes. [Survey Records, 1936]
ii. Cartesian coordinate system
This measure position of a point from a defined geocentric origin using three
perpendicular axes in three-dimensional space denoted by X, Y, Z. The Z axis
aligns with the axis of the rotation of the earth. Conversion between geodetic
coordinates and equivalent Cartesian coordinates use standard mathematical
formulations. [Steed, J., 1990]
In Ghana, Cartesian coordinates are not directly used for mapping but these
are also given when a GPS is used.
iii. Projected coordinates
These are planar two dimensional Cartesian-like coordinates that result from
performing a mathematical mapping from a point on the Earth's surface to a
5
plane. They relate the coordinates of points on earth’s curved surface with
equivalent values of the same points on a plane or flat surface (MAP).
Many such mathematical mappings are used for particular purposes but the type
used in Ghana and many other parts of the World is known as a Transverse
Mercator Mapping. [Survey Records, 1936]
The Transverse Mercator projection is a cylindrical projection. A scale
distortion is often applied along the central meridian to minimize distortions.
The map projection is itself defined by a set of equations for converting
geographic latitudes and longitudes into projected Eastings and Northings and
the reverse calculation. [Snyder, J. P. 1987]
iv. Local coordinates
In Ghana, local coordinates used to be defined for each town based on a
Cassini projection. Usually, a point close to the centre of the town was chosen
with assigned false origin values of X=10000ft. and Y=10000ft. This together
with azimuth bearing formed the basis for plan mapping and drawing of town
maps. [Survey Records, 1936]
v. The Universal Transverse Mercator Coordinate System
This is a specific implementation of the transverse Mercator projection
designed for world wide mapping in which the world is divided into zones
with each zone projected individually. There are 60 zones, each 6° longitude
wide starting from 180° W eastwards. [Toms R.M, 1996]
The UTM is a standard means of stating any position on earth, by giving two
numbers — the easting or how far East, and the northing or how far north. To
use UTM to specify a position anywhere, one must also state the zone. The
origin for each zone is the intersection of the zone central meridian and the
equator.
In the Northern hemisphere, the northing is zero on the Equator and increases
the farther one moves north from the Equator. The false easting on the zone
central meridian is 500000m. Also, the Scale along the central meridian for
each zone is 0.9996
6
In the Southern hemisphere, the false Northing is 10000000m. and the false
Easting is 500000m. [D.M.A., 1990]
1.6 How do we deal with areas of interest that cross UTM zones?
This problem is certain to occur with areas that are wider than 6° of longitude or on
some occasions a narrower area of interest may be bisected by the edge of a UTM
zone
One has to understand that; all features need to be encoded using a consistent
coordinate system so the usual method is to assign the coordinates to the more
predominant zone that contains the majority of the area of interest.
It is also possible to encounter an instance when the area of interest crosses the
Equator. Note that, for Northern zones, the Equator’s Northing is at 0m but for
Southern zones, the Equator’s Northing is at 10,000,000m. Coordinates must be
consistent. But since 0 ≠ 10,000,000, the usual solution is to assign all coordinates to
the Southern zone. [Defence Mapping Agency. (1987)]
2.0 Historical setting- Ghana National Datums for Plan and Height
2.1 National Geodetic framework
Ghana National framework consists of monuments planted at points whose
coordinates are well known. National geodetic controls in Ghana dates back to June
1904 when observations for latitude were done by Captain P. G. Guggisberg from a
pillar established in the house of the then Secretary of Native affairs. This pillar was
connected by traverse to GCS 547 in Accra and longitude for GCS 547 was
determined by telegraphic signal exchange with Cape Town in November 1904 to
establish both latitude and longitude values for GCS 547. The initial pillars, GCS 547
and GCS 121 at Legon were connected by triangulation to obtain their latitudes and
longitudes. [Survey Records, 1936]
7
Other series of triangulation nets were established to obtain latitude and longitude of
trigonometric points. The controls in Northern territories and low lying regions were
established by traversing.
The Main Horizontal framework therefore consists of five chains of triangulations and
a number of primary and secondary traverses.
The Calculated best fitting local ellipsoid that fitted the measurements done in the
then Gold Coast was suggested by the British war office using initial measurements.
This ellipsoid has subsequently come to be known as the War Office ellipsoid.
Ghana Horizontal Datum therefore comprises the following:
• Ellipsoid: Ghana War Office.
– Parameters of ellipsoid :
•
a=6378299.996 m.
•
b=6356751.69 m.
• inverse flattening 296
– feet
log
to
meter
conversion
factor
=
0.304799706846.
f −
= 1.484014544967 (Departmental tables, pg. 48)
m
Easting and northing coordinates are derived using a Transverse Mercator
projection. The origin of coordinates is the intersection of latitude 4 degrees 40 (4o
40’ N) minutes North and longitude 1 degree west (1o W). This origin is given
values of X =0.00 ft., Y = 900000.00ft. To avoid distortion at the eastern and
western edges, true distance along the central meridian was decreased by 1/4000.
This is equivalent to a scale of 0.99975 along the central meridian and this secant
implementation of the transverse Mercator meant that scale is correct at
Y=434000 and Y=1366000. ie 466000 feet from the central meridian on either
side. (W 002o 16’ 50” and E 000o 16’ 50”) [Technical Instructions to Cadastral
Branch]
For Ghana then, horizontal datum is defined by:
• Ellipsoid = War office
• Semi major axis a = 6378299.996m.
• Inverse flattening 1/f = 296
8
• Latitude of origin 4o 40’ N
• Longitude of origin = 1oW
• False Northing = 0.000ft.
• False Easting = 900000.000 ft.
• Scale along central meridian = 0.99975.
This datum definition is referred to in some publications as Accra Datum.
2.2 Local Town Frameworks
Most Town surveys used to be surveyed based on Local Frameworks determined and
executed to fit the whole town before the establishment of the National geodetic
framework.
All towns have now been connected to the main national framework and formulae for
conversion worked out at the Survey Department (SD) Head office for transformation
between such local values and National values.
Because these cover limited regions of within 100km2, affine transformations were
sufficient. Thus for Accra, the conversion formulae are:
X N = 302406.75 + 0.999829 X L + 0.001355YL
YN = 1154438.16 + 0.999829YL − 0.001355 X L
Local to National conversions for other Towns were derived from the equations:
X N = X o + X L + F0 X L + 2 RX L YL (Y0 − 900000)
YN = Y0 + YL + F0YL + R(YL − X L )(Y0 − 900000)
2
2
Where F0 = scale factor defined for each point and R is a constant taken to
be 1.149 x 10-15
In addition, Cartographic charts are designed for use to avoid manual calculations
by the Survey Department. [Survey Records]
2.3 Ghana Vertical Datum
Level heights in Ghana are based on the Mean Sea Level (MSL) determined in Accra
in 1923 as height datum. Very old lines were based on the Low Water of Ordinary
Spring Tide (LWOST) at Accra or LWOST at Sekondi as Datums.
9
Levels based on the LWOST of Sekondi can be converted to MSL by subtracting 3.00
ft. whereas those values based on LWOST of Accra can be converted to MSL by
subtracting 2.85 ft. [Survey Records]
2.4
Leigon Datum
In 1977, the triangulation and traverse framework was re-adjusted together using the
method of variation of Geographical coordinates.
In this new adjustment, it was agreed to hold the existing latitude and longitude values
for GCS 121 (5o 38’ 52.2700” N, 0o 11’ 46.08” W) at Legon as fixed (error free).
Following this adjustment, projected coordinates were computed on the original
Transverse Mercator projection using the same origin of the intersection of latitude 4o
40’N and longitude 1oW to which again false origin values of 0.00 ft. for Northing
and 900000.00 ft. for Easting were adopted.
However ellipsoid proposed for projection was Clarke 1880 because other African
Nations around were already using that ellipsoid and it was thought that this would
enhance future regional mapping. [Survey Records, Aug. 1978]
Although this Re-computation was completed in August 1978, the results were not
implemented or made public. However, it has come to be severally known as “Leigon
datum”, or Ghana metric grid in some publications.
The positional shift between the two different ellipsoidal projected coordinates is on
the average of 25 feet perhaps why it was never used.
3.0 Transformation Procedures for Implementing UTM Coordinates in Ghana
Transformation is the conversion of coordinates from one coordinate system to
another. Transformation involves the conversion of geodetic coordinates from one
geodetic datum to another geodetic datum, usually involving changes in the shape,
orientation, and center position of the reference ellipsoid. [Ayer, J. 2008].
3.1 The Need for Geodetic Transformations
Maps defined on different spheroids, projections, and origin points [datums] can
not be compared or integrated without converting them first into a common datum
10
frame. This conversion is often difficult because there are no immediate tools for
computing them. Performing transformations between planimetric coordinates of
different geodetic datums is a complex mathematically involving process. It
requires the selection of an appropriate transformation model which is
computationally demanding so the availability of user-friendly software tools
would be an advantage. Usually transformations are defined between a
local/National datum and WGS84.
Figure 3 below shows various scenarios of such transformations between the
WGS84 Datum and the Ghana datum.
Figure 3 Transformation Scenarios
In Scenario one, a GPS is used to collect data corrected against a reference
Continuous Reference Station (CORS) data in International Terrestrial Reference
Frame (WGS84). Since the result is also required in WGS84, there is no need for
the application of any transformation.
For Scenario two, the GPS data is corrected against a reference CORS data in
WGS84. If the result is required in the Ghana National Datum (GND), then a
transformation is needed from WGS84 to GND.
In Scenario three, measurements are done in GND. If the results are however
required in WGS84, then again a transformation from GND to WGS84 is
required.
11
In Scenario four, measurements are done in GND and the results are also required
in GND so here again no transformation is required.
Conclusions and Recommendations
As indicative from above, if any new survey is done using a GPS it is possible to
project results directly to UTM once it is properly referenced to a CORS and there
is no need for any transformation.
Alternatively, UTM controls must be established with GPS and thereafter classical
means may be used and results processed in UTM.
A third option is to use classical means of survey based on existing SD controls
and then transform results to UTM. This third option requires that SD must design
software for doing the transformations. In addition Cartographic charts prepared
by SD can also be used to deduce the transformations.
A forth option is when a GPS is used for Surveying but the results are required in
the GND. This is the only option which may require the use of published SD
transformation parameters as input for GPS processing.
In either way, the major challenge is the conversion of all existing data into the
new UTM grid and this requires the software for direct conversion and the
cartographic charts recommended for option three.
References
1. Ayer, J. 2008, Transformation Models and Procedures for Framework
Integration of the Ghana National Geodetic Network, The Ghana Surveyor,
Journal of the Ghana Institution of Surveyors.
2. Ayer, J. and C. Fosu,2008, Map Coordinate Referencing and the Use of GPS
Datasets in Ghana, Journal of Science and Technology, Vol. 28, No. 1
3. D.M.A., Defense Mapping Agency (1990). Datums, Ellipsoids, Grids and
Grid Reference Systems, DMA Technical Manual 83538.1, Fairfax, Virginia,
USA.
4. Defence Mapping Agency. (1987) Department of defence World Geodetic
System 1984: its definition and relationship with local geodetic systems
12
(second edition). Technical Report no. 8350.2, Defence Mapping Agency,
Washington.
5. EPSG 2005. Coordinate Conversion and Transformation including Formulas.
European Survey and Petroleum Group Guidance Note number 7 part 2.
6. Gold Coast Survey Department Records 1936. Report on the
Chains of
Triangulation surveyed in the Colony
7. Primary and secondary adjustment of TM coordinates for Ghana Aug. 1978.
8. Report on Survey and adjustment of a Network of Precise Traverses measured
between 1924 and 1930 in the Colony.
9. Report on the Secondary Triangulation Surveyed in the Colony prior to 1931.
10. Snyder, J. P. 1987. Map Projections – A Working Manual. USGS Prof. Paper
1395
11. Snyder, J.P., Map Projections used by the U.S. Geological Survey, Geological
Survey Bulletin; 1532, U.S. Government Printing Office, Washington D.C.,
1982.
12. Steed, J., (1990). A practical approach to transformation between commonly
used reference systems. Aust. Surveyor, 35(3), 248-264.
13. Survey Records, Aug. 1978 Primary and secondary adjustment of TM
coordinates for Ghana.
14. Tables for Use in the Survey Department.
15. The Technical Instructions for the Cadastral Branch of the Gold Coast Survey,
1936.
16. Toms R.M, Efficient Transformations from Geodetic to UTM Coordinate
Systems, Fifteenth Workshop on Standards for Interoperability of Distributed
Simulations, Orlando, FL, September 16-20, 1996, August 7, 1996
13