Map Projection
|
|
Map projection is a mathematical expression using which the
three-dimensional surface of earth is represented in a two
dimensional plane. The process of projection results in distortion of
one or more map properties such as shape, size, area or direction.
|
|
A single projection system can never account for the
correct representation of all map properties for all the regions of
the world. Therefore, hundreds of projection systems have been
defined for accurate representation of a particular map element for a
particular region of the world.
|
|
Classification of Map Projections
|
|
Map projections are classified on the following
criteria:
|
|
- Method of construction
- Development surface used
- Projection properties
- Position of light source
|
|
Method of Construction
|
|
The term map projection implies projecting the graticule
of the earth onto a flat surface with the help of shadow cast.
However, not all of the map projections are developed in this manner.
Some projections are developed using mathematical calculations only.
Given below are the projections that are based on the method of
construction:
|
|
Perspective Projections : These
projections are made with the help of shadow cast from an illuminated
globe on to a developable surface
|
|
Non Perspective Projections :These projections do not use shadow cast from an
illuminated globe on to a developable surface. A developable surface
is only assumed to be covering the globe and the construction of
projections is done using mathematical calculations.
|
|
Development Surface
Projection transforms the coordinates of earth on to a surface that
can be flattened to a plane without distortion (shearing or
stretching). Such a surface is called a developable surface. The
three basic projections are based on the types of developable surface
and are introduced below:
|
|
1. Cylindrical Projection
|
|
- It can be visualized as a cylinder wrapped
around the globe.
- Once the graticule is projected onto the
cylinder, the cylinder is opened to get a grid like pattern of
latitudes and longitudes.
- The longitudes (meridians) and latitudes
(parallels) appear as straight lines
- Length of equator on the cylinder is equal to
the length of the equator therefore is suitable for showing
equatorial regions.
|
|
|
|
Aspects of
cylindrical projection:
|
|
|
|
(a) Normal: when cylinder has line of tangency to the
equator. It includes Equirectangular Projection, the Mercator
projection, Lambert's Cylindrical Equal Area, Gall's Stereographic
Cylindrical, and Miller cylindrical projection.
|
|
(b) Transverse: when cylinder has line of tangency to
the meridian. It includes the Cassini Projection, Transverse
Mercator, Transverse cylindrical Equal Area Projection, and Modified
Transverse Mercator.
|
|
(c) Oblique: when cylinder has line of tangency to
another point on the globe. It only consists of the Oblique Mercator
projection.
|
|
2. Conic Projection
|
|
- It can be visualized as a cone placed on the
globe, tangent to it at some parallel.
- After projecting the graticule on to the
cone, the cone is cut along one of the meridian and unfolded.
Parallels appear as arcs with a pole and meridians as straight
lines that converge to the same point.
- It can represent only one hemisphere, at a
time, northern or southern.
- Suitable for representing middle latitudes.
|
|
|
|
Aspects of conic
projection:
|
|
|
|
(a) Tangent: when the cone is tangent to only one of the
parallel.
|
|
(b) Secant: when the cone is not big enough to cover the
curvature of earth, it intersects the earth twice at two parallels.
|
|
3. Azimuthal/Zenithal Projection
|
|
- It can be visualized as a flat sheet of paper
tangent to any point on the globe
- The sheet will have the tangent point as the
centre of the circular map, where meridians passing through the
centre are straight line and the parallels are seen as
concentric circle.
- Suitable for showing polar areas
|
|
|
|
Aspects of zenithal projection:
|
|
|
|
(a) Equatorial zenithal: When the plane is tangent to a
point on the equator.
|
|
(b) Oblique zenithal: when the plane is tangent to a
point between a pole and the equator.
|
|
(c) Polar zenithal: when the plane is tangent to one of
the poles.
|
|
Projection
Properties
|
|
According to properties map projections can be
classified as:
|
|
Equal area projection: Also known
as homolographic projections. The areas of different parts of earth
are correctly represented by such projections.
|
|
True shape
projection: Also known as
orthomorphic projections. The shapes of different parts of earth are
correctly represented on these projections.
|
|
True scale or
equidistant projections: Projections that
maintain correct scale are called true scale projections. However, no
projection can maintain the correct scale throughout. Correct scale
can only be maintained along some parallel or meridian.
|
|
Position of
light source
|
|
Placing light source illuminating the globe at different
positions results in the development of different projections. These
projections are:
|
|
Gnomonic projection: when
the source of light is placed at the centre of the globe
|
|
Stereographic Projection: when the source of light is placed at the
periphery of the globe, diametrically opposite to the point at which
developable surface touches the globe
|
|
Orthographic Projection: when the
source of light is placed at infinity from the globe opposite to the
point at which developable surface touches the globe
|
|
|
|
Figure 8: Projections and position of light source
|
|
Constructing Map
Projections
( Adapted from Singh, G 2004, Map work and practical
geography )
|
|
Cylindrical
Projection
Let us draw a network of Simple cylindrical Projection for the whole
globe on the scale of 1: 400,000,000 spacing meridians and parallels
at 30º interval
|
|
Calculations:
|
|
|
|
Figure 9 : Simple cylindrical projection graticule
|
|
|
Steps of
constrution:
|
|
- Draw a line AB, 9.975 cm long to represent
the equator. The equator is a circle on the globe and is
subtended by 360º.
- Since the meridians are to be drawn at an
interval of 30º divide AB into 360/30 or 12 equal parts.
- The length of a meridian is equal to half the
length of the equator i.e. 9.975/2 or 4.987 cm.
- To draw meridians, erect perpendiculars on
the points of divisions of AB. Take these perpendiculars equal
to the length specified for a meridian and keep half of their
length on either side of the equator.
- A meridian on a globe is subtended by 180º.
Since the parallels are to be drawn at an interval of 30º,
divide the central meridian into 180/30 i.e. 6 parts.
- Through these points of divisions draw lines
parallel to the equator. These lines will be parallels of
latitude. Mark the equator and the central meridian with 0º and
the parallels and other meridians. EFGH is the required
graticule.
|
|
Conical
Projection
|
|
Let us draw a graticule on simple conical projection
with one standard parallel on the scale of 1: 180,000,000 for the
area extending from the equator to 90º N latitude and from 60º W
longitude to 100º E longitude with parallels spaced at 15º interval,
meridians at 20º, and standard parallel 45º N.
|
|
Calculations:
|
|
|
|
|
|
Steps of
construction:
|
|
- Draw a circle with a
radius of 3.527 cm that represents the globe. Let NS be the
polar diameter and WE be the equatorial diameter which intersect
each other at right angles at O.
- To draw the standard
parallel 45º N, draw OP making an angle of 45º with OE.
- Draw QP tangent to OP
and extend ON to meet PQ at point Q.
- Draw OA making an angle
equal to the parallel interval i.e. 15º with OE.
|
|
Draw line LM, it represents the central meridian
|
|
|
|
With L as the
centre and QP as the radius, draw an arc intersecting LM at n. This
arc describes the standard parallel 45º N.
|
|
The distance between the successive parallels is 15º.
The length of the arc subtended by 15º is calculated as under:
|
|
|
|
From point n, mark off distances nr, rs, st, nu, uv and
vM , each distance being equal to 0.923 cm. With L as centre, draw
arcs passing through the points t, s, r, u, v and M. These arcs
represent the parallels.
|
|
- Draw OB making an angle
of 20° with OW Length of the arc subtended by 20° is calculated
as under:
|
|
- With O as centre and
radius equal to the arc WB (1.231 cm) draw arc abc.
- From point b, drop
perpendicular bd on line ON. Now db is the distance between the
meridians.
|
|
Keeping in view the number of meridians to be drawn,
mark off distances along the standard parallel toward the east and
west of the point n, each distance being equal to db.
|
|
Join point L with the points of divisions marked on the
standard parallel and produce them to meet the equator.
|
|
|
Figure 10 : Simple conic
projection
|
|
|
Azimuthal
Projection
|
|
Let us draw Polar zenithal equal area projection for the
northern hemisphere on the scale of 1: 200,000,000 spacing parallels
at 15° interval and meridians at 30° interval.
|
|
Calculations:
|
|
|
|
|
|
Steps of construction:
|
|
- Draw a circle with
radius equal to 3.175 cm representing a globe. Let NS and WE be
the polar and equatorial diameter respectively which intersect
each other at right angles at O, the centre of the circle.
- Draw radii Oa, Ob, Oc,
Od, and Oe making angles of 15°, 30°, 45°, 60° and 75°
respectively with OE. Join Ne, Nd, Nc, Nb, Na and NE by straight
lines.
- With radius equal to Ne,
and N’ as centre draw a circle. This circle represents 75°
parallel. Similarly with centre N’ and radii equal to Nd, Nc,
Nb, Na and NE draw circles to represent the parallels of 60°,
45°, 30°, 15° and 0° respectively.
- Draw straight lines AB
and CD intersecting each other at the centre i.e. point N.
- Radius N’B represents 0°
meridian, N’A 180° meridian, N’D 90° E meridian and N’C 90° W
meridian.
- Using protractor, draw
other radii at 30° interval to represent other meridians
|
|
|
Figure 11 : Polar
zenithal equal area projection
|
|
|
Selection of Map
Projection
|
|
Choosing a correct map projection for an area depends on
the following:
|
|
Map Purpose
|
|
Considering the purpose of the map is important while
choosing the map projection. If a map has a specific purpose, one may
need to preserve a certain property such as shape, area or direction
|
|
On the basis of the property preserved, maps can be
categorized as following
|
|
a. Maps that preserve shapes.
|
|
Used for showing
local directions and representing the shapes of the features. Such
maps include:
|
|
- Topographic and
cadastral maps.
- Navigation charts (for
plotting course bearings and wind direction).
- Civil engineering maps
and military maps.
- Weather maps (for
showing the local direction in which weather systems are
moving).
|
|
b. Maps that preserve area
|
|
The size of any
area on the map is in true proportion to its size on the earth. Such
projections can be used to show
|
|
- Density of an attribute
e.g. population density with dots
- Spatial extent of a
categorical attribute e.g. land use maps
- Quantitative attributes
by area e.g. Gross Domestic Product by country
- World political maps to
correct popular misconceptions about the relative sizes of
countries.
|
|
c. Maps that preserve scale
|
|
Preserves true
scale from a single point to all other points on the map. The maps
that use this property include:
|
|
- Maps of airline
distances from a single city to several other cities
- Seismic maps showing
distances from the epicenter of an earthquake
- Maps used to calculate
ranges; for example, the cruising ranges of airplanes or the
habitats of animal species
|
|
d. Maps that preserve direction
|
|
On any Azimuthal projection, all azimuths, or
directions, are true from a single specified point to all other
points on the map. On a conformal projection, directions are locally
true, but are distorted with distance.
|
|
General purpose maps
|
|
There are many projections which show the world with a
balanced distortion of shape and area. Few of these are Winkel Tripel, Robinson and
Miller Cylindrical.For larger-scale maps, from continents to large
countries, equidistant projections are good at balancing shape and
area distortion. Depending on the area of interest, one might use
Azimuthal Equidistant,Equidistant Conic and Plate Carrée.
|
|
Study area
|
|
Geographical location
|
|
The line of zero distortion for a cylindrical projection
is equator. For conical projections it is parallels and for Azimuthal
it is one of the poles. If the study area is in tropics use
cylindrical projection, for middle latitudes use conical and for
Polar Regions use Azimuthal projections.
|
|
Shape of the area
|
|
Young in 1920 described a way of selecting the map
projection which is known as Young’s rule. According to this rule, if
the ratio of maximum extent (z) (measured from the centre of the
country to its most distant boundary) and the width (δ) of the
country comes out to be less than 1.41, Azimuthal projection is
preferable. If the ratio is greater than 1.41 a conical or
cylindrical projection should be used.
Z/δ < 1.41 Azimuthal Projection
Z/δ >1.41 Conical or Cylindrical projections
|
|
Projection Systems
|
|
Given below is the description of the projection systems
that are mostly used:
|
|
Cylindrical Projection
|
|
I.
Equirectangular projection
This is a Projection on to a cylinder which is tangent to the
equator. It is believed to be invented by Marinus of Tyre, about A.D.
100.
|
|
|
|
- Poles are straight lines equal in length to
the equator
- Meridians are straight parallel lines,
equally spaced and are half as long as the equator. All
meridians are of same length therefore scale is true along all
meridians.
- Parallels are straight, equally spaced lines
which are perpendicular to the meridians and are equal to the
length of the equator.
- Length of the equator on the map is the same
as that on the globe but the length of other parallels on map is
more than the length of corresponding parallels on the globe. So
the scale is true only along the Equator and not along other parallels.
- Distance between the parallels and meridians
remain same throughout the map.
- Since the projection is neither equal area
nor orthomorphic, maps on this projection are used for general
purposes only.
|
|
II. Lambert's cylindrical equal-area projection
|
|
It is devised by JH Lambert in 1772. It is a normal
perspective projection onto a cylinder tangent at the equator
|
|
|
|
- Parallels and meridians are straight lines
- The meridians intersect parallels at right
angles
- The distance between parallels decrease
toward the poles but meridians are equally spaced
- The length of the equator on this projection
is same as that on globe but other parallels are longer than
corresponding parallels on globe. So, the scale is true along
the equator but is exaggerated along other parallels
- Shape and scale distortions increase near
points 90 degrees from the central line resulting in vertical
exaggeration of Equatorial regions with compression of regions
in middle latitudes
- Despite the shape distortion in some portions
of a world map, this projection is well suited for equal-area
mapping of regions which are predominantly north-south in
extent, which have an oblique central line, or which lie near
the Equator.
|
|
III. Gall's
stereographic cylindrical projection
Invented by James Gall in 1855, this projection is a cylindrical
projection with two standard parallels at 45ºN and 45ºS.
|
|
|
|
- Poles are straight lines.
- Meridians are straight lines and are equally
spaced.
- Parallels are straight lines but the distance
between them increases away from the equator.
- Shapes are true at the standard parallels.
Distortion increases on moving away from these latitudes and is
highest at the poles.
- Scale is true in all directions along 45ºN
and 45ºS.
- Used for world maps in British atlases.
|
|
IV. Mercator projection
|
|
Gerardus Mercator in 1569 invented this projection.
|
|
|
|
- Parallels and meridians are straight lines
- Meridians intersect parallels at right angle
- Distance between the meridians remains the
same but distance between the parallels increases towards the
pole
- The length of equator on the projection is
equal to the length of the equator on the globe whereas other
parallels are drawn longer than what they are on the globe, therefore
the scale along the equator is correct but is incorrect for
other parallels
- As scale varies from parallel to parallel and
is exaggerated towards the pole, the shapes of large sized
countries are distorted more towards pole and less towards
equator. However, shapes of small countries are preserved
- The image of the poles are at infinity
- Commonly used for navigational purposes,
ocean currents and wind direction are shown on this projection
|
|
V.Transverse Mercator
|
|
This projection results from projecting the sphere onto
a cylinder tangent to a central meridian.
|
|
|
|
- Only centre meridian and equator are
projected as straight lines. The other parallels and meridians
are projected as curves.
- The meridians and the parallels intersect at
right angles
- Small shapes are maintained but larger shapes
distort away from the central meridian.
- The area distortion increases with distance
from the central meridian
- Used to portray areas with larger north-south
extent. British National Grid is based on this projection only.
|
|
Pseudo-cylindrical Projections
A pseudo cylindrical projection is that projection in which latitudes
are parallel straight lines but meridians are curved.
|
|
I. Mollweide Projection
|
|
|
|
- The poles are points and the central meridian
is a straight line
- The meridians 90º away from central meridians
are circular arcs and all other meridians are elliptical arcs.
- The parallels are straight but unequally
spaced.
- Scale is true along 40º 44' North and 40º 44'
South.
- Equal –area projection
- Used for preparing world maps
|
|
II. Sinusoidal Projection
|
|
|
|
- The central meridian is a straight line and
all other meridians are equally spaced sinusoidal curves.
- The parallels are straight lines that
intersect centre meridian at right angles.
- Shape and angles are correct along the
central meridian and equator
- The distortion of shape and angles increases
away from the central meridian and is high near the edges
- Equal area projection
- Used for world maps illustrating area
characteristics. Used for continental maps of South America,
Africa, and occasionally other land masses, where each has its
own central meridian.
|
|
III. Eckert VI
|
|
|
|
- Parallels are unequally spaced straight
lines.
- Meridians are equally spaced sinusoidal
curves.
- The poles and the central meridians are
straight lines and half as long as equator.
- It stretches shapes and scale by 29% in the
north-south direction, along the equator. This stretching
reduces to zero at 49º 16' N and 49º 16' S.
- The areas near the poles are compressed in
north-south direction.
- Suitable for thematic mapping of the world.
|
|
Conical Projection
|
|
I. Bonne’s Projection
|
|
|
|
- Pole is represented as a point and parallels
as concentric arcs of circles
- Scale along all the parallels is correct
- Central meridian is a straight line along
which the scale is correct.
- Other meridians are curved and longer than
corresponding meridians on the globe. Scale along meridians
increases away from the central meridian
- Central meridian intersects all parallels at
right angle. Other meridians intersect standard parallel at
right angle but other parallels obliquely. Shape is only
preserved along central meridian and standard parallel
- The distance and scale between two parallels
are correct. Area between projected parallels is equal to the
area between the same parallels on the globe. Therefore, is an
equal area projection
- Maps of European countries are shown in this
projection. It is also used for preparing topographical sheets
of small countries of middle latitudes.
|
|
II. Polyonic Projection
|
|
|
|
- The parallels are arcs of circles with
different centers
- Each parallel is a standard parallel i.e.
each parallel is developed from a different cone
- Equator is represented as a straight line and
the pole as a point
- Parallels are equally spaced along central
meridian but the distance between them increases away from the
central meridian.
- Scale is correct along every parallel.
- Central meridian intersects all parallels at
right angle so the scale along it, is correct. Other meridians
are curved and longer than corresponding meridians on the globe
and so scale along meridians increases away from the central
meridian.
- It is used for preparing topographical sheets
of small areas.
|
|
Azimuthal/Zenithal Projection
|
|
I. Polar Zenithal
Equal area projection
This projection is invented by J.H Lambert in the year 1772. It is
also known as Lambert’s Equal Area Projection.
|
|
|
|
- The pole is a point forming the centre of the
projection and the parallels are concentric circles.
- The meridians are straight lines radiating
from pole having correct angular distance between them.
- The meridians intersect the parallels at
right angles.
- The scale along the parallels increases away
from the centre of the projection.
- The decrease in the scale along meridians is
in the same proportion in which there is an increase in the
scale along the parallels away from the centre of the
projection. Thus the projection is an equal area projection.
- Shapes are distorted away from the centre of
the projection. Scale along the meridians is small and along the
parallels is large so the shapes are compressed along the
meridians but stretched along the parallels.
- Used for preparing political and distribution
maps of polar regions. It can also be used for preparing general
purpose maps of large areas in Northern Hemisphere.
|
|
II. Polar Zenithal Equidistant Projection
|
|
|
|
- The pole is a point forming the centre of the
projection and the parallels are concentric circles.
- The meridians are straight lines radiating
from pole having correct angular distance between them.
- The meridians intersect the parallels at
right angles.
- The spacing between the parallels represent
true distances, therefore the scale along the meridians is
correct.
- The scale along the parallels increases away
from the centre of the projection.
- The exaggeration and distortion of shapes
increases away from the centre of the projection.
- The projection is neither equal area nor
orthomorphic.
- It is used for preparing maps of polar areas
for general purposes.
|
|
III. Gnomonic
Projection
It is also known as great-circle sailing chart.
|
|
|
|
- The pole is a point forming the centre of the
projection and the parallels are concentric circles.
- The meridians are straight lines radiating
from pole having correct angular distance between them.
- The meridians intersect the parallels at
right angles.
- The parallels are unequally spaced. The
distances between the parallels increase rapidly toward the
margin of the projection. This causes exaggeration of the scale
along the meridians.
- The scale along the parallels increases away
from the centre of the projection.
- The exaggeration and distortion of shapes
increases away from the centre of the projection. The
exaggeration in the meridian scale is greater than that in any
other zenithal projection.
- It is neither equal area nor orthomorphic.
- An arc on the globe which is a part of a
great circle is represented as a straight line on this
projection. This is because the radii from the centre of the
globe are produced to meet the plane placed tangentially at the
pole.
- It is used to show great-circle paths as
straight lines and thus to assist navigators and aviators in
determining appropriate courses.
|
|
IV. Sereographic Projection
|
|
|
|
- The pole is a point forming the centre of the
projection and the parallels are concentric circles.
- The meridians are straight lines radiating
from pole having correct angular distance between them.
- The meridians intersect the parallels at
right angles.
- The parallels are unequally spaced. The
distances between the parallels increase toward the margin of
the projection. The exaggeration in the meridian scale is less
than that in the case of Gnomonic projection.
- The scale along the parallels also increases
away from the meridian and in the same proportion in which it
increases along the meridians. At any point scale along the
parallel is equal to the scale along the meridian.
- The areas are exaggerated on this projection
and the exaggeration increases away from the centre of the
projection.
- A circle drawn on the globe is represented by
a circle on this projection.
- It is used to show world in hemispheres. Also
used for preparing aeronautical charts and daily weather maps of
the polar areas.
|
|
V. Orthographic Projection
|
|
|
|
- The pole is a point forming the centre of the
projection and the parallels are concentric circles.
- The meridians are straight lines radiating
from pole having correct angular distance between them.
- The meridians intersect the parallels at
right angles.
- The parallels are not equally spaced. The
distances between them decrease rapidly towards the margin of
the projection. So, the scale along the meridians decreases away
from the centre of the projection.
- The scale along the parallel is correct.
- The distortion of the shapes increases away
from the centre of the projection.
- It is neither equal area nor orthomorphic.
- The projection is used to prepare charts for
showing the celestial bodies such as moon and other planets.
|
|
Universal
Transverse Mercator (UTM)
|
|
UTM projection
divides the surface of the Earth into a number of zones, each zone
having a 6 degree longitudinal extent, Transverse Mercator projection
with a central meridian in the center of the zone. UTM zones extend
from 80 degrees South latitude to 84 degrees North latitude. The
zones are numbered from west to east. The first zone begins at the
International Date Line (180°).
|
|
|
|
|
|
The particular
transverse Mercator map that is used to represent each zone has its
central meridian running north-to-south down the center of the zone.
This means that no portion of any particular zone is very far from
the central meridian of the transverse Mercator map that is used to
depict the zone. Since a Universal Transverse Mercator zone is 6° of
longitude wide, no portion of a UTM zone is more than 3° of longitude
from the zone's central meridian. Since the distortion in a
transverse Mercator map is relatively low near the map's central
meridian, the result of this close proximity to the map's central
meridian is that the transverse Mercator map used to depict each zone
within the coordinate system contains relatively little distortion.
|
|
Adapted from the
report Map Projections of Europe (2001), the table gives
an account of the commonly used projection systems.
|
|
Property
|
Developable surface
|
Aspect
|
Projections
|
Extent of use
|
Conformal
(True shape)
|
Cylinder
|
Normal
|
Mercator
|
Equatorial regions
(east-west extent)
|
Transverse
|
UTM (Universal
Transverse Mercator)
|
Whole world except polar
areas
|
Oblique
|
Rosenmund Oblique
Mercator
|
Small regions, oblique
& east - west extent
|
Cone
|
Normal
|
Lambert Conformal Conic
|
Small regions, oblique
& east - west extent (1 or 2 standard parallels)
|
Plane
|
Any
|
Stereographic
|
Small regions upto the
hemisphere
|
Polar
|
UPS (Universal Polar
Stereographic)
|
Polar regions
|
Homolographic
(Equal area)
|
Cylinder
|
Normal
|
Lambert Equal Area
|
Equatorial areas
(east-west extent)
|
Cone
|
Normal
|
Albers Equal Area
|
Smaller regions &
continents with east-west extent
|
Plane
|
Any
|
Lambert Azimuthal Equal
Area
|
Smaller regions about
same north-south , east-west extent
|
Equatorial
|
Hammer-Aitoff
|
World
|
Equidistant
|
Cylinder
|
Normal
|
Plate Caree
|
World
|
|
Transverse
|
Cassini Soldner
|
Locally used for large
scale mapping
|
Cone
|
Normal
|
Equidistant Conic
|
Smaller regions &
continents with (1 or 2 standard parallels) east-west extent
|
Plane
|
Any
|
Azimuthal Equidistant
|
Smaller regions about
same north-south , east-west extent
|
|
|
|
|
|
|
Georeferencing
|
It is a process of locating an entity in real world
coordinates. It aligns geographic data to a known coordinate system
representing earth defined through projection systems so it can be
viewed, queried, and analyzed with other geographic data.
|
To georeference a geographic data, the positions of
known points, called control points, are determined. The Ground
Control Points (GCPs) are defined as the points with known
geographical location, whose positions on map correspond to their
positions on earth. GCPs are collected from fixed objects and are
marked on the data to be georeferenced that define where the data is
on earth. The whole data adjusts itself according to these GCPs. At
least three control points are required for georeferencing a
data. Additional control points help increasing the
accuracy. Once the data is georeferenced, each point has a
coordinate associated to it which means the location of any object in
the data (map) can now be determined
|
|
Figure 16 :
Georeferencing an image
|
|
In brief, the process of georeferencing establishes
control points; inputs the known geographic coordinates of these
control points, chooses the coordinate system and other projection
parameters and then minimizes residuals. Residuals are the difference
between the actual coordinates of the control points and the
coordinates predicted by the geographic model created using the
control points. The residuals help in determining the level of
accuracy of the process. The quality of the rectification depends
on the number, accuracy, and distribution of the control points and
the choice of transformation model.
|
Georeferencing Raster Data
|
The raster data occupy a raster space which is defined
as a digital image of the arrangement of the pixels in a grid. The
computer reads the header of the data file and determines the
dimension of the raster space. The position of a cell in the raster
space is referenced by row and column (row, column). This system of
referencing the raster cells is called raster coordinate system. The
origin for this system lies at the upper left corner of the monitor
because the computer monitor displays an image from left to right and
from top to bottom. The method of referencing positions in raster
space is different from that on maps. The origin of the map
coordinate system is the lower left corner. To visualize raster data
spatially raster data needs to be transformed into a map coordinate
system i.e. the raster coordinates (row, column) are transformed into
corresponding ground coordinates (East, West).
|
There are two approaches to georeferencing:
- Image to Map
rectification
- Image to Image
registration
|
Image to Map: rectification is the process by which geometry of
an image is made planimetric. It involves the measurement of the
image coordinates of the reference cells (GCPs) and the corresponding
ground coordinates to relate the image with the real world. The two
sets of coordinates are used to solve a set of polynomial equations
whose order depends upon the amount of geometric distortion in the
raster image. Generally affine transformation is used for the
purpose:
|
|
Where X and Y are the ground coordinates, and, x and y
are the image coordinates.
|
A minimum of three GCPs are required to solve these
equations, though greater the GCPs more accurate is the
rectification. On solving the equations we get the values of the six
(a, b, c, d, e and f) coefficients. Any image coordinate can then be
substituted in the equations to get the corresponding ground position
on the used map coordinate system. The positions of the original grid
cells will have to be interpolated in the mapping coordinate system.
|
After the coordinate transformation the raster cells may
have been oriented differently than the way they were originally in
the raster coordinate system. The attribute value is to be
interpolated for the cells oriented to the new coordinate system.
This is called resampling. There are three common methods of
resampling:
|
Nearest neighbor: In this method, the attribute value of
the original cell nearest to a cell in the output raster layer is
assigned to the corresponding cell.
|
- Bilinear interpolation:
It assigns the value to a cell in the output raster layer by
taking weighted average of the surrounding four cells in the
original grid nearest to it.
- Cubic convolution: It
assigns the value to a cell in the output raster layer by taking
weighted average of the surrounding sixteen cells in the
original grid nearest to it.
|
Among these three the nearest neighbor is preferred
because it doesn’t change the values of the original grid cells
assigned to the reoriented grid cells but it produces blocky images.
The cubic convolution on the other side does change the values but it
generates smoother images.
|
The result of an image to map rectification is a geometrically
correct grid of raster cells.
|
Image to Image registration
is a method of georeferencing a raster layer with the help of another
raster, which is already georeferenced by the process of image to map
rectification. The already georeferenced raster is used as a
reference for the raster which is to be rectified. Control points are
selected from the two raster layers, the coordinate transformation
and resampling is then done in the similar manner as it is done for
image to map rectification. Image to Image registration is used to
spatially match several raster layers to a single reference raster
layer.
|
|
|
|
Comments
Post a Comment