Surveying II AAiT

ammanuel bekele
DEPARTMENT OF CIVIL ENGINEERING, AAIT 1 of 74
Table of Contents
1 Curves ............................................................................................................................................... 3
1.1 General ....................................................................................................................................... 3
1.2 Types of curves and their use ..................................................................................................... 3
1.3 Circular Curve ............................................................................................................................ 4
1.4 Setting Out of Circular Curves................................................................................................... 6
1.5 Transition Curves ....................................................................................................................... 9
1.6 Setting Out of Transition Curves from the Tangent Point ....................................................... 12
1.7 Vertical Curves ........................................................................................................................ 13
2 Control Survey ................................................................................................................................ 17
2.1 Introduction .............................................................................................................................. 17
2.2 Triangulation ............................................................................................................................ 18
2.3 Classification of Triangulation Systems .................................................................................. 18
2.4 Triangulation figures and choices ............................................................................................ 19
2.5 Angle and side conditions ........................................................................................................ 20
2.6 Strength of Figures ................................................................................................................... 21
2.7 Adjustment of Triangulation networks .................................................................................... 25
2.8 Triangulation Computation ...................................................................................................... 34
2.9 Trilateration .............................................................................................................................. 35
2.10 Adjustment of Trilateration .................................................................................................. 37
2.11 Computation of lengths and positions .................................................................................. 38
3 Contour Lines and digital terrain model ......................................................................................... 39
3.1 Contours and Contour Lines .................................................................................................... 39
3.1.1 General Rules for Contour ................................................................................................ 40
3.1.2 Methods of Contouring ..................................................................................................... 40
3.1.3 Uses of Contour Lines (Contour Map) ............................................................................. 41
3.2 Digital Terrain Model (DTM) .................................................................................................. 41
3.2.1 Overview........................................................................................................................... 41
3.2.2 DTM Creation ................................................................................................................... 42
3.2.3 Derivation from DTM....................................................................................................... 43
3.2.4 Application of DTM ......................................................................................................... 43
3.2.5 Uncertainty and Errors in DTMs ..................................................................................... 43
4 Photogrammetry ............................................................................................................................. 44
4.1 Introduction .............................................................................................................................. 44
4.2 Types of aerial photographs ..................................................................................................... 44
4.3 Scale of Aerial Photographs ..................................................................................................... 44

4.4 Relief Displacement ................................................................................................................. 47
4.5 Stereoscopy .............................................................................................................................. 48
4.6 Parallax..................................................................................................................................... 50
4.7 Use and Products of Aerial photographs.................................................................................. 55
4.8 Applications of Aerial photographs ......................................................................................... 56
5 Introduction to Global Positioning System (GPS) Surveying ........................................................ 57
5.1 Over of GPS ............................................................................................................................. 57
5.2 GPS Signal ............................................................................................................................... 58
5.3 Reference coordinate system for GPS ..................................................................................... 59
5.4 Fundamentals of GPS positioning............................................................................................ 61
5.5 Errors in GPS observation........................................................................................................ 62
5.6 Linear combinations of GPS observables ................................................................................ 63
5.7 GPS Positioning Modes ........................................................................................................... 64
5.8 Other Satellite navigation systems ........................................................................................... 66
5.9 Application of GPS .................................................................................................................. 66
6 Introduction to Geographic Information System (GIS) .................................................................. 67
6.1 Introduction .............................................................................................................................. 67
6.2 Definitions of GIS .................................................................................................................... 67
6.3 Key elements of GIS ................................................................................................................ 67
6.4 History of GIS .......................................................................................................................... 67
6.5 Components of GIS .................................................................................................................. 68
6.6 Data Model in GIS ................................................................................................................... 69
6.6.1 Vector model .................................................................................................................... 70
6.6.2 Raster model ..................................................................................................................... 71
6.6.3 Comparison of Vector and Raster..................................................................................... 72
6.6.4 Triangulated Irregular Network (TIN).............................................................................. 73
6.7 Data Capturing and Entry techniques in GIS ........................................................................... 73
6.8 Application areas of GIS .......................................................................................................... 74
1 Curves
1.1 General
In the geometric design of motorways, railways, pipelines, etc., the design and setting out of curves is
an important aspect of the engineer’s work. The initial design is usually based on a series of straight
sections whose positions are defined largely by the topography of the area. The intersections of pairs of
straights are then connected by horizontal curves. In the vertical design, intersecting gradients are
connected by curves in the vertical plane.
1.2 Types of curves and their use

In general curves can be listed under three main headings, as follows:
(1) Circular curves of constant radius.
(2) Transition curves of varying radius (spirals).
(3) Vertical curves of parabolic form.
Curves used in horizontal planes to connect two straight tangent sections are called horizontal curves.
Two types are used: circular arcs and spirals. Both are readily laid out in the field with surveying
equipment. A simple curve [Figure 1.1(a)] is a circular arc connecting two tangents. It is the type most
often used. A compound curve [Figure 1.1(b)] is composed of two or more circular arcs of different
radii tangent to each other, with their centers on the same side of the alignment. The combination of a
short length of tangent (less than 100 m) connecting two circular arcs that have centers on the same
side [Figure 24.1(c)] is called a broken-back curve. A reverse curve [Figure 1.1(d)] consists of two
circular arcs tangent to each other, with their centers on opposite sides of the alignment.
Figure 1.1 Types of Circular Curve
Compound and reverse curves are unsuitable for modern high-speed highway, rapid transit, and
railroad traffic and should be avoided if possible. However, they are sometimes necessary in
mountainous terrain to avoid excessive grades or very deep cuts and fills. Compound curves are often
used on exit and entrance ramps of interstate highways and expressways, although easement curves are
generally a better choice for these situations.
Easement curves are desirable, especially for railroads and rapid transit systems, to lessen the sudden
change in curvature at the junction of a tangent and a circular curve. A spiral makes an excellent
easement curve because its radius decreases uniformly from infinity at the tangent to that of the curve
it meets. Spirals are used to connect a tangent with a circular curve, a tangent with a tangent (double
spiral), and a circular curve with a circular curve. Figure 1.2 illustrates these arrangements.
Figure 1.2 Use of spiral transition curves
The effect of centrifugal force on a vehicle passing around a curve can be balanced by superelevation,
which raises the outer rail of a track or outer edge of a highway pavement. Correct transition into
superelevation on a spiral increases uniformly with the distance from the beginning of the spiral and is
in inverse proportion to the radius at any point. Properly superelevated spirals ensure smooth and safe
riding with less wear on equipment. As noted, spirals are used for railroads and rapid-transit systems.
This is because trains are constrained to follow the tracks, and thus a smooth, safe, and comfortable
ride can only be assured with properly constructed alignments that include easement curves. On
highways, spirals are less frequently used because drivers are able to overcome abrupt directional
changes at circular curves by steering a spiraled path as they enter and exit the curves.
1.3 Circular Curve

A simple circular curve shown in Fig. 1.3, consists of simple arc of a circle of radius R connecting two
straights AI and IB at tangent points T1 called the point of commencement (P.C.) and T2 called the
point of tangency (P.T.), intersecting at I, called the point of intersection (P.I.), having a deflection
angle  or angle of intersection . The distance E of the midpoint of the curve from I is called the
external distance. The arc length from T1 to T2 is the length of curve, and the chord T1T2 is called the
long chord. The distance M between the midpoints of the curve and the long chord, is called the mid-
ordinate. The distance T1I which is equal to the distance IT2, is called the tangent length. The tangent
AI is called the back tangent and the tangent IB is the forward tangent.
Figure 1.3 Elements of circular curve
Formula to calculate the various elements of a circular curve for use in design and setting out, are as
under.
Tangent length (T) =
Length of curve ( )
Long chord (L) =
External distance (E) = ( )
Mid-ordinate (M) = ( )
Chainage of T1 = chainage of P.I. – T
Chainage of T2 = chainage of T1 +
Degree of the circular curve

The rate of curvature of circular curves can be designated either by their radius (e.g., a 1500-m curve)
or by their degree of curve. There are two different designations for degree of curve, the arc definition
and the chord definition. By the arc definition, degree of curve is the central angle subtended by a
circular arc of 100 m [see Figure 1.4]. This definition is preferred for highway work. By the chord
definition, degree of curve is the angle at the center of a circular arc subtended by a chord of 100 m
.This definition is convenient for very gentle curves and hence is preferred for railroads. The formulas
relating radius R and degree D of curves for both definitions are shown by
Figure 1.4 Degree of a curve using arc definition
1.4 Setting Out of Circular Curves
One approach for setting out circular curves is to locate the tangent points and intersection points that
define the construction lines and then proceed with the setting out of the pegs along the curve by taking
measurements from these positions.
The pegs that mark the center line of the curve can be set out from the tangent points:
1. By tape and off set
2. By angle and distance measurement
3. By angle measurement only
4. By Coordinates
1. Methods Using Tape
I. Offset from Long Chord: the method is suitable for curves of short radius.
AB = AO – OB
= AO  OU 2  UB 2
= R R 2  (L / 2) 2
Draw CE parallel to TU, then
y = EB = EO – BO
EO2 = CO2 – CE2 EO = R2  x2
y  R 2  x 2  R 2  (L / 2) 2
Fig 1.5 Offset from long chord

II. Offset from Tangent: the method is suitable for short curve
R 2  AB 2  AO 2  y 2  ( R  x) 2
Therefore
1
x2 2
y  R R 2  x 2  R  R(1  )
R2
Expanding by using Binomial theorem,
1 x2
y  R  R(1   ....)
2 R2
x2
 Approximately
2R
Fig 1.6 Offset from tangent
2. Setting out using one thodolite and tape by deflection angle method
Arc TA = 2R ( in radians)

 2 R *
180
= chord TA (if chord not greater
than (R/20))
TA *180
Therefore  deg
2R
TA
  1718.9 min
R
Fig. 1.7 Setting out by deflection angle method
3. Setting out using two theodolite
Fig. 1.8 Setting out by two theodolite method

4. Setting out using Coordinate

Today, because of the availability of total station and GPS/GNSS instruments with data collectors,
circular curves are often staked using the coordinate method. For this procedure, coordinates of the
points on the curve to be staked must first be determined in some reference coordinate system. In
Figure bellow, assume that the azimuth of the back-tangent going from A to V is known, the
coordinates of the PI (point V) are known, and that the defining parts of the curve have been
computed using the following equations. Using the tangent distance and azimuth of the back
tangent, after the departure and latitude are computed, the coordinates of the point of
commencement (A) calculated by
XA = XV + T sin AzVA
YA = YV + T cos AzVA
where AzVA is the back azimuth of line AV.
With the coordinates of the PC known, coordinates of points on the curve can be computed using the
same deflection angles and subchords used to stake out the curve by the total chord method. Deflection
angles are added to the azimuth of AV to get azimuths of the chords to all stations to be set. Using the
total chord length and chord azimuth for each station, departures and latitudes are calculated, and
added to the coordinates of A (the PC) to get the station coordinates.
With coordinates known for all curve points, they can be staked with the total station occupying any
convenient point whose coordinates are also known in the same system. The PC, PT, PI, or curve
midpoints are points that are often used.
AzAO = AzAV + 90°

XO = XA + R sin AzAO
YO = YA + R cos AzAO
XP = XO + R sin AzOP
YP = YO + R cos AzOP
Fig 1.9 Setting out by coordinate method

To stake the curve points, the total station instrument is set on the curve’s center point, a backsight
taken on point A, and the azimuth of line OA indexed on the horizontal circle. To stake any point such
as P, the azimuth of OP is placed on the instrument’s circle and the stake placed on the line of sight at
a distance R from the instrument.
5. Setting out with inaccessible intersection point
Fig 1.10 Setting out with inaccessible intersection point
1.5 Transition Curves

The transition curve is a curve of constantly changing radius. It used to connect a straight to a curve of
radius R, then the commencing radius of the transition will be the same as the straight (∞), and the final
radius will be that of the curve R.
The centrifugal force (i.e. P=WV2/gR) increases as R decreases, and thus, since the radius of the
straight is infinity, the centrifugal force would increase instantaneously from zero to its maximum
value (assuming no change in V), as the vehicle moved from the straight to curve. Passengers in the
vehicle would thus experience a lateral shock as the tangent was passed. To avoid this curve variable
radius is inserted between the straight and the circular curve in order that the centrifugal force may
build up in a gradual and uniform manner. This curve is called a transition curve.
Fig 1.11 Super elevation

Derivation of Transition Curve Equations

 P increases uniformly with distance l from the beginning of the transition curve, i.e. P  l
 At the point where the radius of the transition curve = r
WV 2 1
P Hence P for constant
gr r
velocity
And so l  1 / r since P  l
 lr = K = LR
Fig 1.12 Derivation of transition curve

Length of Transition Curve
Let L = total length of the transition curve, R= radius of circular curve entered and V= uniform
velocity of the vehicle.
2
The radial acceleration at the circular cure = f  V
R
Therefore if the time taken to travel along the transition curve be t seconds,
t=L/V
V2 L V3
Rate of change radial acceleration = f /( L / V )    / 
3
R V LR
Note: A rate of change of radial acceleration of 1/3 m/s is a COMFORT LIMIT above which side-
throw will be noticed
K l
dl  rd  d  d  dl
l K
2
l
On integrating,    A A=0, since =0
2K
when l = 0
l2 l2 l2
  
2 K 2 LR 2lr
Now we can write
dy  dl cos(  d / 2)
dx  dl sin(  d / 2)
l3
Integrating, x  and y  l
6K
Consider points X and Y which is near together
on the transition curve, X being l from T.
Relation Between  and 
x l 2 (1   2 / 14   2 / 440) 
tan    
y 6K (1   2 / 10   2 / 216) 3
For small angle  tan  =  therefore, the maximum  = 1 = L/2R
Shift
Where transition curves are introduced between the tangents and a circular curve of radius R, the circular curve
is “shifted” inwards from its original position by an amount BP = S so that the curves can meet tangentially.
This is equivalent to having a circular curve of radius (R+S) connecting the tangents replaced by two transition
curves and a circular curve of Radius R, although the tangents are not the same, being B and T respectively.
Referring the figure below <NMO = <RTO = 90
BM = NT1 = Maximum offset on transition curve
Also shift = S = BP  BM  MO  NT1  ( PO  MO)
L3 L3  1 1
2 4
 L3 1
2
=  ( R  R cos1 )    R  R(1    ...)  R

6LR 6LR  2! 4!  6 LR 2
L3 R L2 2
=  ( )
6LR 2 2 LR
L2
=
24 R
RL2 L
Also Q1T1  PT1  R1 = 
2 LR 2
Fig 1.13 Derivation of shift
Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

1.6 Setting Out of Transition Curves from the Tangent Point
Step 1
L2
Calculate the shift S from the expression S 
24 R
Step 2

Calculate IB = ( R  S ) tan
2
Step 3
 L
BT = L/2 then IT = ( R  S ) tan 
2 2
Step 4
Either calculate offset from
l3 y3
x or x 
6K 6K
OR
l2
Calculate the deflection angles  from particular distances l from T using the fact that    /3  rad
6 RL
1800 l 2 1800 L
= min  max  L / 6R rad OR =
 RL R
Step 5
Setting out the circular curve
Calculate the deflection angles for the circular curve from
c
 Circular  1718.9 min
R
Where c is the chord length for the circular arc
The angle subtended by the circular arc T1T2 at the center of that arc is ( - 21)

Length of the arc T1T2 = R(  21 )
180
Step 6
Set up the theodolite at T1 and sight back on T; then , transit the telescope and locate T1T3 by setting off an
angle 2/3 1. Set out the circular curve by the deflection angle method, from this tangent.

Step 7
Setting the other transition curve
Point U is located from I using the relationship IU = IT. The transition is then set out from the tangent point U
and tangent UI by either of method given in step 4.
1.7 Vertical Curves

Curves are needed to provide smooth transitions between straight segments (tangents) of grade lines for
highways and railroads. Because these curves exist in vertical planes, they are called vertical curves. An
example is illustrated in Figure below, which shows the profile view of a proposed section of highway to be
constructed from A to B. A grade line consisting of three tangent sections has been designed to fit the ground
profile. Two vertical curves are needed: curve a to join tangents 1 and 2, and curve b to connect tangents 2 and
3. The function of each curve is to provide a gradual change in grade from the initial (back) tangent to the grade
of the second (forward) tangent. Because parabolas provide a constant rate of change of grade, they are ideal
and almost always applied for vertical alignments used by vehicular traffic.
Two basic types of vertical curves exist, crest and sag. Curve a is a crest type, which by definition undergoes a
negative change in grade; that is, the curve turns downward. Curve b is a sag type, in which the change in grade
is positive and the curve turns upward. There are several factors that must be taken into account when designing
a grade line of tangents and curves on any highway or railroad project. They include (1) providing a good fit
with the existing ground profile, thereby minimizing the depths of cuts and fills, (2) balancing the volume of cut
material against fill, (3) maintaining adequate drainage, (4) not exceeding maximum specified grades, and (5)
meeting fixed elevations such as intersections with other roads. In addition, the curves must be designed to (a)
fit the grade lines they connect, (b) have lengths sufficient to meet specifications covering a maximum rate of
change of grade (which affects the comfort of vehicle occupants), and (c) provide sufficient sight distance for
safe vehicle operation.

Fig 1.14 Vertical curve: summit and sag
Elevations at selected points (e.g., at 20, 30, or 40 m) along vertical parabolic curves are usually computed by
the tangent-offset method. It is simple, straightforward, conveniently performed with calculators and computers,
and self-checking. After the elevations of curve points have been computed, they are staked in the field to guide
construction operations so the route can be built according to plan.
Curve Design
For flat gradient it is sufficiently accurate to treat the length along the tangents, the length along the curve, a
chord AC as equal to length 2l so as to be in accordance with the survey information.
The basic equation for a simple parabola is:
Y=C.Xs2
Where Y is the vertical offset from gradient to curve, Xs is the distance from the start of the curve, and C is a
constant.
Vetical curve
+ q%
+ p% -q%
+ p%
Difference in gradients (p - q)
Difference in gradients (p+ q)
-p%
-p%
+ q%
-q%
Difference in gradients (p+ q) Difference in gradients (p - q)
Fig. (b)
Fig. (a)
where p and q are in %

The properties of the parabola give: -

 The vertical through the intersection of the tangents B is a diameter and bisects AC.
 BD=DE. D is the vertex of the parabola.
 The offsets from the tangent AB are proportional to the square of distance from A.
 Offsets from tangent at D are proportional to the square of distance from D.
B'
( p* L / 100)
(p/100)
B
+ p% (p+ q)%
Xs y h1
D ( q* L / 100)
-q%

H1
A Horizontal Horizontal C
L L
Fig. (c)
B
+ p% -q%
D
h h
A
S E C
Fig. (d)
From Fig (c)

y  xs
2
y  C * xs
2
y  C ( x sec ) 2
y  Kx 2
When x=2L
Y=B’C=PL/100 + qL/100
( p  q) L
y  Kx 2  K (2 L) 2 
100
pq
K
400 L
( p  q) 2
y  Kx 2  x
400 L

( p  q) L2 ( p  q) L
When x = L, Y= BD  
400 L 400
pL
Height of B above A =
100
pL qL
Height of C above A = 
100 100
The angle between AC and the horizontal
L ( p  q) p  q
  radians
100 2 L 200
pL ( p  q) ( p  q)
BE   L L  2 BD
100 200 200
 BD  DE
For other case where the difference in gradients = p-q
( p  q) x 2
y
400 L
Height of a point
xp ( p  q) x 2
H1  
100 400 L
For maximum value of H1,
dH1
0
dx
p 2( p  q ) x
 0
100 400 L
2p
x L
Pq

2 Control Survey
2.1 Introduction
The determination of the precise position of a number of stations, usually spread over a large area, is referred to
as control surveying. Control surveys can be horizontal or vertical. The objective of horizontal control surveys
is to establish a network of control stations whose position is specified in terms of latitude and longitude, where
as the objective of vertical control survey is to provide the elevations of fixed benchmarks with respect to the
mean sea level datum. These control surveys are located where other surveys can be conveniently and
accurately tied in to them.
The results of the control surveys are used as a basis from which surveys of smaller extent can be originated.
Boundary surveys, construction, route surveys, topographic and hydrographic surveys and other may be
involved. Horizontal control can be carried out by precise traversing, by triangulation, by trilateration, and
perhaps by same combination of these methods. Intersection, resection, and satellite positioning are also the
other methods in horizontal control surveys. The exact method used depends of the terrain, equipment available,
information needed and economic factors.
With traversing, a series of horizontal distances and angles are measured. This is generally cheaper due to
shorter sights and convenience to carry it out under less favorable weather conditions than the other methods. Its
disadvantage is that, there are fewer checks available for locating mistakes in the work and the whole system
can rather easily sway or bend. To check a traverse it is necessary to form a loop returning to its starting point
or to tie it in to previously established control points.
A triangulation consists of a series of joined or overlapping triangles in which an occasional line (called the
base line) is measured and all other sides of the triangles are calculated from angles measured at the vertices of
the triangles. The lines of a triangulation system form a network that ties together all the triangulation stations at
the vertices of the triangles. A triangulation has the following advantages:
1. More redundancies or checks are available i.e. more than one route can be followed to compute the
length of a line.
2. There is little tendency for the system to sway or bend i.e. azimuths can be easily and accurately carried
or established throughout system.
3. Outstanding landmarks such as steeples, water tanks, etc can be located by establishing directions from
different stations.
Its disadvantages are it needs long-range intervisibility, which in turn requires the erection of special towers and
signals, making the system the most expensive. Moreover, a good weather is required to attain intervisibility.
With trilateration, the lengths of the sides of a series of joined or overlapping triangles are measured (usually
with the EDM equipment) and the angles are computed from the lengths. It has the following advantages:
1. It is more accurate than the other two types due to the fact that distances can be measured more
accurately than angles.
2. It is generally less expensive than triangulation.
3. More checks are available.

Unlike triangulation, it is not easy to position transmission towers, steeples, water tanks, etc. by the EDM
because to do so requires reflectors on these landmarks. These landmarks can be located if angular
measurements are made to them.
In combined triangulation and trilateration systems all sides and angles in the joined or over lapping triangles
are measured. This method provides the strongest control network.
2.2 Triangulation
Triangulation as a form of horizontal control is applied when a large area is to be surveyed and when the
methods of traversing would not be expected to maintain a uniformly high accuracy over the entire area. The
methods of triangulation require a maximum number of precise angle measurements and a minimum number of
distance measurements. The triangles are developed in to a net of interconnected figures, and certain lines,
called base lines, must be measured in order to compute the other sides in the net.
Triangulation is necessary to control the location of large bridge structures, state and federal highways, dams,
canals, and other engineering works of a massive nature. A large project, such as boundary location, power
development, water resources development, flood control, irrigation, or reclamation, requires triangulation in
order to maintain the necessary accuracy throughout the system.
Once a triangulation system, whether large or small, has been developed, measured and adjusted, the points in
the system then furnish control for subsequent traversing, minor triangulation, trilateration, intersection, and
resection needed for day-to-day engineering operations.
2.3 Classification of Triangulation Systems
Survey or accuracy standards are generally defined as being the minimum accuracies deemed necessary to
obtain certain specific objectives. Survey specifications can be defined, as the field operations or “recipes”
needed to achieve the particular standards desired.
Accuracies required for horizontal control depend on the type of survey and the ultimate use of the control
points. There are three orders of triangulation based on their uses: first order (primary), second order
(secondary) and third order (tertiary) triangulation.
First order is the highest accuracy and is required for developing the national network of horizontal control, for
the study of small crustal movements in areas of seismic activity, and for large metropolitan control expansion.
Since it covers a large area, the effect of earth’s curvature is to be taken in to account-geodetic triangulation. It
may cover the whole country primary grids may be provided.
Second order (secondary) triangulation provides point at greater density than first order triangulation. This
network is adjusted to fit its parent primary triangle or its surrounding primary control. Less refinement is
needed as the network is surrounded by the primary control. It is recommended for controlling extensive land
subdivision and construction.
Third order (tertiary) triangulations used to establish control for local developments and improvements,
topographic and hydrographic surveys or other such projects for which they provide sufficient accuracy. They
are not carried out for rural areas. They may not be adjusted to fit the national network.

2.4 Triangulation figures and choices

Although triangles are the basic figures in triangulation networks, they are not allowed to exist alone because
they do not provide sufficient checks on measurements and subsequent computations. Triangles are combined
to form other geometrical entities. The most common figure used is the braced quadrilateral (fig1.1b).It is best
suited to long narrow systems. Braced polygons (fig1.1c) can also be used in the case of wide systems.
(a) Chain of triangles (b) chain of quadrilaterals (c) chain of polygons.
Fig 2.1 triangulation figures
Triangulation procedure
1. Reconnaissance, select the locations of stations;

2. Evaluation of the strength of figures;
3. Erection of signals, and in some cases, tower for elevating the signals and /or instruments;
4. Observation of directions or angles;
5. Measurement of base lines;
6. Astronomic observation at one or more locations,
7. Computations including reduction to sea level, calculation of the lengths of all sides and coordinates for
all stations, and adjustment of the triangulation network to provide the best estimates of co-ordinates of
all points.
Reconnaissance
The success of any triangulation depends on the reconnaissance, which is the most difficult and exacting task of
the extensive survey.
Reconnaissance consists of selection of stations, determination of the size and shape of the resulting triangles,
the number of station to be occupied and the number of angle or direction to be observed. The intervisibility and
accessibility of stations, the usefulness of stations in latter work, the cost of the necessary signals, and the
convenience of base-line measurements are considered.
In order to acquire the information, available maps, survey information and aerial photographs of the area are
used. Field inspection and rough measurements may also be carried out to supplement the study, select
favorable location for station and check suitability.
Heavily wooded country is the most difficult in which to carry triangulation. A valley of proper width, with
peaks on either side, is the most favorable. The stations are located on higher points, provided that their
locations will give the best-shaped triangles. The ideal condition is one in which no or low towers are used and
no clearing is required.

A most important and difficult part is to determine the height of towers necessary to make the sight between any
two of the stations clear of obstruction.
A mistake on the part of the person may lead to a delay of the work, especially in the first and second order
triangulation work where observation is done after dark.
2.5 Angle and side conditions
To ensure homogenous results from the computations of the sides of the triangles in a triangulation net, the
network must have adequate geometric strength. Two important factors that affect geometric strength are the
angle and side conditions in the network, and the magnitude of the angles observed.
Angle condition equations in a figure express the following:
1. The sum of the interior angles in a polygon equals some multiple of 1800.
2. A station equation exists if an angle  1 can be expressed in terms of an angle 1at a station.
3. If a horizon is closed, a center-point equation which states that the sum of these angles is equal to 3600 is
required.
To see how angle conditions are determined, consider fig 1.2 below. Start with one line a and no angle. Add
another line b at an angle 1from a. Next, add another line c to yield 2. This process continues adding line n to
yield (n-1) angles, so that n lines result in (n-1) angles to produce a determinate figure.
b a
1 1
a a a c
2
Fig. 2.2 Angle Condition
Any angles measured in excess of these required for the determinate case are redundant and need a condition
equation. For example, for a plane triangle, if three angles are measured, there is one redundancy for which the
angle condition is
1+2+3=1800
It is necessary to evaluate the number of angle conditions in a given figure . Let
CA=total number of angle conditions (including center-point equations)

L= number of lines in a polygon;
A= number of angles measured in the polygon.
Then CA=A-(L-1)=A-L+1
Angle condition equations can be satisfied without having consistent lengths in sides therefore, a side condition
equation is necessary when lengths of a side of a triangle can be computed in more than one route using the law
of sines.
To locate any third station in a triangle three lines a needed, and two lines are needed to locate any additional
station. All other lines in excess of there are redundant and need a condition equation.

Fig 2.3 (a) Inconsistent Sides (b) Side Condition
The number of extra lines or the number of side conditions is

CS = n’ – 3 – 2 (s-3) = n’ – 2s + 3
Where CS = number of side conditions;
n’= number of side in a figure;
S = number of stations in the figure.
The total number of side and angle conditions is there fore
C = CA + CS = (A – L + 1) + (n1 – 2S + 3 ) or
C = n - no where n-no of angles measured,
no - minimum no of angles required to plot the figure.
2.6 Strength of Figures
When a triangulation project is being evaluated in the preliminary stages of the work, it is necessary to
determine the strength of figure for the network. This step is required in order to ensure uniform accuracy
throughout the network. The strength of figure is a function of:
1. The geometric strength of the triangles that make up the network. Ideally, the triangles should be
equilateral.
2. The number of stations occupied for angle or direction measurements lines occupied at only one end
should be avoided whenever possible.
3. The number of angle and side conditions used in adjusting the network. this number should be large in
proportion to the numbers of observations.
The following equation is used to evaluate the strength of figure:
R
D C
D
 
 A  A  B  B
2 2

Where R = Strength of figure
D = the number of directions observed in the figure.
C = total number of angle and side conditions
δ A and δB are the respective differences of the sines.
The subscripts A and B refer to the two distance angles in each of the triangles under consideration.
For a triangle with all stations occupied
C = CA + CS = (A – L + 1) + (n1 – 2s + 3)
where n1 is the total number of lines including the known side in a figure, A is the number of lines observed
in both direction including the known side, s is the total number of stations and L is the number of station
occupied.

In the above equation, the two factors (D – C)/D and ( A   A  B   B ) are related only to the number of
2 2
conditions and observations and the geometry of the triangles. Thus the value for R is independent of the
precision of the measurements. Consequently this procedure is useful primarily as a means of comparing
various network configurations so as to obtain optimum geometric conditions and a desirable number of
conditions versus number of measurements. It is also used as a means for determining the most favouable
route for calculating through a network of triangulation.
Note:
1.The stronger the figure, the lower the value of R becomes.
2.The A angle is opposite the side to be computed, whereas the B angle is opposite the known side
3.The third angle in each of the triangles is referred to as the azimuth angle the strength of the figure does
not depend on the size of the azimuth angle, only on the A and B.
Example:
Compute the first and second-best strength for the following figure AB is the starting side.
C D
B A Fig 2.4
Solution;
There are 10 directions as the line AB is fixed
Therefore, D = 10.
Also A = 8; L = 6; n1 = 6; s = 4
The total number of side and angle conditions is
C = A-L+n1 –2S + 3
=8 – 6 + 1 + 4 – 2* 4+3 =4
The values of  A   A  B   B are obtained from table-1.2 and are tabulated below for the distance angles
2 2
A and B. the smaller of the two angles is read from the top row.
R
DC
D
  
 A2   A  B   B 2
D  C 10  4
  0.6
D 10
Triangle Known Computed Distance Angles  A2   a b   b 2

Side Side
ABD AB BD 1120 390 5.4
BCD BD CD 74 80 1.0
6.4*0.6=3.84=R2
ABC AB AC 105 47 3.2
ACD AC CD 85 63 1.4

4.6*0.6=2.76=R1
ABD AB AD 31 39 28.6
ACD AD CD 85 33 11.6
40.2 *0.6=24.12
ABC AB BC 27 47 29
BCD BC CD 74 24 28.4
54.4*0.6=32.64
DC
As the number of conditions increases versus the number of measurements, the term decreases.
D
Also as the distance angles A and B approach the ideal value of 450 , the term in the summation
become smaller. The stronger the figure, the smaller the value of R.
Triangulation Stations
The points forming the triangulation stations are selected on the basis of visibility as for example on the top of
hills or church steeples, or radio towers, or water tanks. As a result, the points are not uniformly spaced. Some
of the points are obviously inaccessible and it is necessary to establish eccentric stations from them and to
determine the distance and directions from the main stations to the eccentric stations.
Sometimes it is necessary to build special towers for making the observations. These towers contain one tower
built inside another so that the towers supporting the instrument and the instrument operator are independent.
For triangulation of lower accuracy, a small pole signal about 2m high guyed in place, or an object already in
place might be used for sighting. As most triangulation works of higher accuracy are carried out between late
afternoon and night, an electric lamp or an automobile headlight may be used as a signal. Signal should be free
from phase.
In general, the type of signal used depends on the length of the line and the accuracy required whereas its form
depends on the locality and available materials.
Triangulation stations should be marked and referenced very carefully for use at later days. The signals used
must permit centering the instrument if the station is to be occupied (for instance, an iron pipe set vertically so
that a pole can be inserted in to it for sighting and this pole can be removed when centering the instrument over
the station)
When towers are required due to flat terrain, heavy timber or other factors, it is necessary to determine the
minimum height of towers required from the equation for earth’s curvature and atmospheric refraction.
Example:
Determine the height of a tower to permit seeing the surface of the earth at a distance of 10 km.
H = 0.0675d2 where d is the distance in km and h is the height of the tower in m.
If towers of equal height are placed on both ends of the line, each tower is 5km from the midpoint and the
height of each tower is h = 0.675* 25 = 1.69m.

To minimize the effects of refraction the line of sight should clear the ground by at least 3.5m, so each tower is
5.2 m high.
Baseline and Angle Measurement
Since the computed sides of a triangulation system can be no more accurate than the baselines, every precaution
to ensure accuracy is taken in measuring these lines. The length of any baseline is determined primarily by the
desirability of strong figures in the net. Ordinarily, the longer the baseline, the easier it will be found to secure
strong figures. Short baselines are connected to the triangulation system through a basenet.
Because of errors in angular measurements, the accuracy of the computed lengths will decrease as the distance
from the baseline increases. On extensive triangulation the accuracy is maintained by measuring additional
based lines. The number of baselines required will depend on the strength of the figures through which
computations are performed. When a maximum allowable cumulated value for R has been reached, a new
baseline is required. The maximum values to R for the various orders of triangulate are given in table 1.1
Baselines can be measured either by tapes or by EDM equipment. For measurement of higher precision, invar
tape is always used. Some of the activity required while taping are:
 Clearing obstructions along the baseline:
 Setting stakes on line at distances apart equal to the length of the tape to be used;
 Fastening copper strips to the tops of these stakes for marking the tape lengths;
 providing support for the tape;
 Running levels to determine difference in elevation where required;
 Recording field temperature and pull and applying systematic corrections etc.
Once systematic corrections have been applied to the measured lengths, the remaining step is to reduce all the
triangle lengths to equivalent sea level distances in all extensive triangulation surveys.
In the figure, A and B are two points on the earth’s surface at an
B average elevation of h above sea level.
D h Since the lengths of the arcs are proportional to their radii,
MSLD S D DR
 or s 
R R D h R h
A
Where S = sea level length, m
D = measured length at average elevation of
H above sea level .
Fig 2.5 Reduction to sea level R = 6,372,200m
Example:
The average elevation measured at the site of a measured baseline is 540m above mean sea level. The following
additional data are given the tape has a standard length of 50.002m at 20 o c when supported at the 0 and 50 m
points with a tension of 9kg; the coefficient of thermal expansion of the tape is 3.6* 10-7 / oC and the tape
weighs 1.5kg, the recorded length of the baseline is 1050.241m; the average temperature was 25.5oC ; the cross-

sectional area of the tape is 0.06cm2; the stakes were set on a 3% slope at 50-m intervals; the standard interval
and tension were used throughout the operation. Compute the length of the baseline.
Solution
The corrections to be applied are: corrections for absolute length, temperature, slope and reduction to sea level
distance.
1. Correction for absolute length:
c (50.002  50.000)
C a  L *1050.214   0.042m
l 50
2. Correction for temperature:
Ct   L(T Ts )  3.6 *107 (25.5  20)   0.002m
Total slope distance = 10.50.241 + 0.002 = 1050.285m
3. Correction for slop:
Cslope= s(1 – cos a ) = 1050.285(1 – cos (tan –1 3/100)) = -0.472m
Reduced horizontal length of the base = 1050.285 – 0.472 = 1049.813m
4. Reduction to sea level distance:
DR 1049 .813 * 6372200
Sea level distance, S   1049 .724 m
Rh 6372200  540
The use of EDM equipment for base line measurement only needs a clear line of sight between two ends of a
base. Provided that such a base is available in the main triangle, there is no need for a basenet. The operation
with EDM equipment has the advantage that the time for observation is independent of the length of the base.
Moreover, measurement is entirely independent of the type of terrain, making the possibility of selecting
baselines at more frequent intervals to increase the strength of the figure.
Instruments used for measuring angles depend on the desired accuracy of the triangulation network For first
order work, direction theodolites which can be read directly to 0.2 should be used .For second-order work, a
repeating theodolite reading to 20’’ to 30’’ can be used if advantage is taken of the added precision gained by
repeating angles.For more information on thiodolite surveying the student is referred to surveying I.
2.7 Adjustment of Triangulation networks
The primary objective of triangulation is to determine the co-ordinates of the triangulation stations for further
surveys. Before computations of the positions carried out, the net should be adjusted as follows.
1. All baselines have to be corrected for systematic errors and reduced to sea level distances;
2. Observed directions and angles should be reduced to average values and triangles checked for closure;
3. Necessary reductions to center should be made where eccentric stations have been occupied;
4. Spherical excess should be computed where lines and triangles are of sufficient length.In the adjustment
process, the following conditions have to be satisfied:
(1)Angle conditions (2) side conditions (3) length conditions (4) azimuth conditions and (5) position
conditions.
There as to methods of adjustment: the approximate method and the least square method.

Reduction to center
In measuring triangulation angles, it may not be possible to occupy a station directly below an excellent target,
such as a church steeple or a lighthouse. Targets may also be blown out of position and angles read on them
have to be corrected to the true position of the triangulation station. There are thus two types of problems: when
the instrument position is not the true position and the measured angles should be corrected to what they would
be at the true station: and when he target is blown out of position.
In fig 1.6 shown below, the angles were measured at an eccentric station A’ instead of the true station A
B
D A
A'
E Fig 2.6 C
Steps to calculate the correction to the measured angles:

1. Calculate the directions of all lines with respect to the line A’A’ as a meridian.
2. Find the corrections that refer these directions to the true station by solving the triangles such as ABA’,
ACA’, ADA’, and AEA’,
l L l
In triangle ABA’   sin e  sin a
sin e sin a L
l can be measured directly, and made as small as possible. L is not known and is taken as A’B which is
estimated from ABA’, taking A’ as the triangulation station. Since the values of the angles at B, C, D and E in
the above figure are very small, their values can be obtained in seconds of angle by dividing both sides of the
above equation by sin 1’’ Thus, if e denotes the angle in seconds,
l sin a
e' '  *
sin 1 L
Example:
A theodolite is set at an eccentric station E and the following observations are made to four triangulation
stations. Compute the directions from A to B, C, D and G. prepare an abstract of directions, taking the
direction from A to B as 00 00’00.0’’.
AE = 1.5 m A B
 C
E
G D

To station Mean direction

A 2470 19’32.8’’
B 18 o52’ 53.2
C 37002’44.6’’
D 1200 45’16.6
G 1540 13’26.0
Solution
1. Calculate the directions of all lines with respect to EA as a meridian.
All directions have to be increased by 3600 - 2470 19’328’’ = 1120 40’27.2’’
Line Mean direction (new)
EA 00 00’00.0’’
EB 1310 33’20.4’’
EC 1490 43’11.8’’
ED 2330 25’43.8’’
EG 2660 53’ 53.2’’
2. Since the distances L1, L2, L3 and L4 from station A to stations B, C, D and G are not known, take them to be
approximately equal to distances from station E to stations B, C, D and G respectively.
3. Calculate the corrections that refer these lines to station A as follows:
Stations B C D G
Distance 8670m 7950m 4755m 4650m
0 0
Angle a 131 33’20.4 149 43’11.8’’ 233’25’43.8’’ 266013’53.2’’
Sin a + 0.74831 + 0.50423 -0.80312 -0.99784
6 6 6
1/sin1” 0.3094*10 0.3094*10 0.3094*10 0.3094*106
e’’ + 26.8’’ +19.6’’ -52.2’’ -66.4’’
A (corrected) 131033’47.2’’ 149043’31.4’’ 233024’51.6’’ 266012’46.8’’
Spherical Excess
When a theodolite is set up and leveled at a triangulation, the vertical axis of the instrument assumes a plumb
position and the lime of sight sweeps through a vertical plane. The vertical plane containing a line of sight to a
distant station intersect the sea level surface along an arc.
A B’
C’
Cc B
Fig 2.7 Spherical excess
In fig above, the vertical plane containing the line of sight from A to B intersects the sphere along arc AB; the
direction reading on the theodolite’s horizontal circle is represented by line AB’. The horizontal angle obtained

from the two direction readings is equal to the spherical angle i.e. <C’AB’ = <CAB .The same situation would
exist at stations B and C. Thus the theodolite measures the spherical angles at the triangulation stations forming
the spherical triangle. The amount by which the sum of three angles of a spherical triangle exceeds 180 0 is
called the spherical excess. Before the closing error in the angles of a large triangle can be determined, spherical
excess must be known.
The amount of this excess is dependent on the area of the triangle and the latitudes of the vertexes. It is
approximately equal to 1’’ for every 75.5 mile 2 of area.
A more exact value for the excess in seconds of angle is
bc sin A 1
e' '  2
 bc sin A where m  2
2R 2 sin 1' ' 2 R sin 1' '
Where b, c and A are the two sides and the included angle of the triangle, and R is the radius of
curvature of the earth. The values of m are available from tables.
The excess is divided equally between the three angles of the triangle, the corrections being subtracted from the
observed values.
Triangulation Adjustment by the Approximate Method
Connections in triangulation nets can be either by means of triangles or quadrilaterals, or by a combination of

different kinds of figures. When triangles are used the only adjustment to be made is to satisfy the condition that
the sum of the three angles of each triangle should equal 1800 plus any spherical excess. The adjustment is made
ordinarily by correcting each angle by one – third of the closure in that triangle. Although the work involved is
simple, the only check on the computed lengths is afforded by measuring additional base lines.
When quadrilaterals are used, additional checks on angular measurements are provided by the added geometric
conditions in each figure. Checks on the computed lengths are afforded by following four different
combinations of triangles in each quadrilateral. Unless the measured angles have been adjusted, it is unlikely
that the four computed lengths will agree exactly. The purpose of the adjustment is to correct the observed
angles so that the various conditions relative to the sums of the angles will be satisfied, and in addition, so that
the computed lengths of a side will have the same numerical value regardless of the triangles used in the
computation.
D C
h g f e
a d
c
A b B
Fig 2.8
In the quadrilateral shown above, the sum of the eight angles should equal 360 0, and the sum of the three angles
of any one triangle should equal 1800, ignoring the spherical excess. In addition, b + c = g + f, and a + h = d + e.
The side condition can be developed as follows:
sin b
In triangle ABC, BC = AB
sin e

sin d sin b sin d

In triangle BCD, CD = BC  AB
sin e sin e sin g
sin f sin b sin d sin f
In triangle CDA, DA = CD  AB
sin a sin a sin e sin g
sin h sin b sin d sin f sin .h

In triangle DAB, AB  DA  AB
sin c sin a sin c sin e sin g
sin b sin d sin f sin h

Hence, 1
sin a sin c sin e sin g
Taking the logarithms of both sides of the equation,

(log sin b + log sin d + log f + log sin h) -
(log sin a + log sin c + log sin e + log sin g) = 0
Unless this condition is satisfied, the value of the computed length will depend on the triangles used in the
computation.
Example:
Adjust the quadrilateral shown below by the approximate method. All angles are observed with equal precision.
Angle Value
D 1 270 02’12.6’’
C 5 6 2 53041’ 09.3’’
4 7
3 550 25’25.1’’
4 43051’17.0’’
5 56011’57.3’’
6 24031’25.2’’
7 39053’01.4’’
3
2 1 8 8 59023’35.7’’
B A
Fig 2.9
Solution:
Sum of the eight angles in the quadrilateral = 3600 00’03.6”
Correction to each angle for 3600= 3.6’’/ 8 = 0.45’’
Sum of opposite angles (1 + 2) = (5 + 6)
1+2 = 80043’21’’ 5+6=80043’21.6’’
Correction to be applied to these angles = 0.6’’ / 4 =0.15’’(to be added to 1 and 2, and deducted from 5 and
6).
Sum of opposite angles (3+4) = (7+8)
3+4 = 99016’41.20’’; 7+8 = 99016’36.20’’

Correction to be applied to these angles = 5’’ /4 = 0.125’’(to be added to 7 and 8, and subtracted from 3 and
4).
The side condition equation is developed as
(log sin 1+log sin 3+log sin 5+log sin 7)-(log sin 2+log sin 4+log sin 6+log sin 8) = 0
The solution is effected by using the method indicated in the table below. The sums of the log sins of alternate
angles adjusted for angular closure are found as in column 4.
Difference between sums of log sins of alternate angles = 0.00001205
The sums of the alternate log sines fail to agree by 120.5 in seven place values, and hence the angles must be
corrected by the amount that will eliminate this difference. The change must be effected by correcting each by
equal amount in order to keep the sum at 3600 and to keep the sums of opposite angles equal. The value of the
correction in seconds is found by dividing the total change in log sines by the total difference in log sines for 1’’
of angle.
The correction is to be added to those angles whose log sines are to be raised and subtracted from those angles
whose log sines are to be brought down. If the corrections are accurate, the sums of the log sines of the
corrected angles will agree within the last one or two in seven decimal place.
Triangulation Adjustment by the Least squares Method: condition Equations
In a completed quadrilateral with six sides and eight angles observed, the conditions to be satisfied are (1) the
angles in any three of the four triangles must add up to 1800 plus spherical excess, and (2) the angles of the
closing side of the quadrilateral must be the same no matter which way it is computed. The method of condition
(1) (2) (3) (4) (5) (6) (7) (8)

Stat Ober. angles For For Logsine of Difference of CorrectedCorrected logsines
ion 3600 opp. adj. angles logsines for 1” angles
(ang angles
les)
1 27002’12.6” 12.15” 12.30” -0.342406963 0.000004125 12.93” -0.3424040363
0.000001547 08.37” -0.093783478
2 53041’9.3” 08.85” 09.00” -0.093782503
3 55025’25.1” 24.65” 23.40” -0.084407123 0.000001450 24.03” -0.084406209
4 43051’17.0” 16.55” 15.30” -0.159375662 0.000002191 14.67” -0.159377043
5 56011’57.3” 56.85” 56.70” -0.080411744 0.000001409 57.33” -0.080410856
6 24031’25.2” 42.75” 24.60” -0.381882389 0.000004614 23.97” -0.381885297

0.000002518 02.83” -0.192981444
7 39053 ’01.4 00.95” 02.20” -0.192983031
0.000001245 35.87” -0.065157042
8 59023’35.7’’ 35.25” 36.50” -0.065156257
-0.700208862  00.00” -0.700208
3600 0’ 03.6” 00.00” 00.00” -0.700196812 =0.000019099 73
equations for the adjustment of the quadrilateral is explained by a specimen problem.

Example
Adjust the quadrilateral in the example given previously by the method of condition equations.
Solution:
Let the residuals of the measured angles 1 through 8 be v1, v2 v3 v4 v5 v6 v7 and v8 respectively.
The resulting angles including the correction are
1. 27002’12.6’ + v1 2. 53041’09.3’’ + v2 3. 55025’25..1’’+v3 4. 43051;17.0’’ + v4

5. 56011’57.3’’ +v5 6. 24031’25.2’’+v6
0
7. 39 53’01.4’’+ v7 8. 59023’35.7’’+v8
Considering angle conditions of three of the four triangles
Triangle ABC, A 27002’12.6’’+ v1

B 53041’09.3’’+ v2 + 55025’25.1’’ + v3
C 43051’17.0 + v4
180000’04.0’’ + v1 +v2+ v3 + v4 = 180000’00.0’’
Triangle BCD, B 55025’25.1+v3

C 43051’17.0’’+ v4 + 56011’57.3’’ + v5
D 24031’25.2’’ + v6
180000’04.6’’ + v3 + v4 + v5+v6 = 180000’00.0’’
Triangle CDA, C 56011’57.3’’ + V5

D 24031’25.2 + V6 + 39053’01.4’’ + V7
A 59023’35.7’’ + V8
179059’59.6’’ + V5 + V6 + V7 + V8 = 180000’00.0’’
From the three conditions we have,

V1 + V2 + V3 + V4 + 4,0 = 0 (1)
V3 + V4 + V5 + V6 + 4.6 = 0 (2)
V5 + V6 + V7 + V8 - 0.4 = 0 (3)
The side condition can be developed in a number of ways:
1. Considering all the four triangles ABC, BCD, CDA and DAB in turn
sin 1sin 3 sin 5 sin 7

1
sin 2 sin 4 sin 6 sin 8
AB Sin 4 AB sin 7 AC sin 6  7 

2.  ;  ; 
AC Sin2  3 AD sin 2 AD sin 5

sin 2 sin 4 sin 6  7 

 1= Angles 1 and 8 are omitted here.
sin 2  3sin 5 sin 7
DA sin 2 CD sin 8 DB sin 4  5

3.  ;  ; 
BD sin 1  8 DA sin 5 CD sin 3
sin 2 sin 4  5sin 8

 1= Angles 6 and 7 ar omitted here.
sin 3 sin 5 sin 1  8
CD Sin8 BC Sin6 CA Sin2  3
4.  ;  ; 
CA Sin6  7  CD Sin3 BC Sin1
sin 2  3sin 6 sin 8
1  Angle s4 and 5 are omitted here.
sin 1sin 3 sin 6  7 
BC sin 1 BA sin 7 BD sin 4  5 sin 1sin 4  5sin 7

5.  ;  ;  1 Angles 2 and 3 are
AB sin 4 BD sin 1  8 BC sin 6 sin 4 sin 6 sin 1  8
omitted here.
The selection of one of the five sets is in such a manner as to make the constant term in the side condition
equation as large as possible. For this, the set that contains the smallest combination of angles is chosen. In this
specific problem, set 5 is selected because the sum of the omitted angles 2 and 3 is larger than the sum of the
two angles any other station
Sin 1 sin(4+5) sin 7=sin 4 sin 6 (1+8)
Taking the logarithm of both sides
log sin 1 + log sin (4+5) + log sin 7 = log sin 4 +log sin 6 + log sin (1+8)
In order to make a condition equation, the corrections to the angles are included.
log sin (1 + v1) + log sin (4 + 5 + v4 + v5) + log sin (7+v7)=
log sin (4 + v4) + log sin (6 +V6) + log sin (1 + 8 + V1 + V8)
If the corrections are in seconds, each logarithm can be expressed in a more convenient form, as for instance,
log sin (1 +V1) = log sin 1 +(change in value of log sin 1 per second)*V1
log sin (1+V1) = -0.342405725 + 4. 126*10-6V1

log sin (4+5 +V4 + V7) = -0.00672087551 –0.373*10-6V4 -0.373*10-6V5
log sin (7 + V7) = -0.192985047 + 2.520* 10-6 V7
-0.542111647 + 10-6 (4.126v1 – 0.373v4 – 0.373v5+2.520v7)
log sin (4 + V4)= -0.159371936 + 2.191* 10-6V4

log sin (6+V6) = -0.38187962 + 4.615*10-6 V6
log sin (1+8 + V1 + V8) = -0.000843514 + 0.131*10-6v1+0.131*10-6 v8
-0.54209507 + 10-6(2.191v4 + 4.615v6 +0. 131v1 + 0.131v8)
Transposing all terms to the left and simplifying, we have

3.995v1 – 2.564v4 – 0.373v5 – 4.615v6 + 2.52v7 – 0.138v8 – 16.577 = 0 (4)

Equations 1 through 4 form the four condition equations.
If we assume that all angles are measured with equal precision, the weight given to each angle is unity and
hence ∑V2 must be a minimum in order to arrive at the most probable values for the angles in the solution
.Hence
V12 +V22+V32 +V42 +V52 +V62 +V72 +V8 2= minimum (5)
The four condition equations contain eight unknowns resulting in no unique solution. The condition of least
squares allows us to find a unique set of V’s that satisfy the four condition equations. Let the four condition
equations be multiplied by constants C1, C2, C3, and C4, respectively, known as Lagrange multipliers.
C1 (V1 + V2 + V3 + V4 + 4.0) = 0 (1a)

C2 (V3 + V4 + V5 + V6 + 4.6) = 0 (2a)
C3 (V5 + V6 + V7 + V8 - 0.4) = 0 (3a)
C4 (3.995V1 – 2.564V4 – 0.373V5 – 4.615V6 + 2.52V7 – 0.138V8 – 16.577) = 0(4a)
Adding the left-hand side of the four equations and differentiating with respect to the variables we get,
dv1(C1 + 3.995C4) + dv2(C1) + dv3 (C1 + C2) + dv4 (C1 + C2 – 2.564C4) + dv5 (C2 + C3 –0.373C4) + dv6 (C2 +
C3 – 4.615C4) + dv7(C3 + 2.52C4 )+ dv8 (C3 – 0.131C4) – 16.577 = 0
Differentiating equation 5 and comparing to that above result, we get
V1 = C1 + 3.995C4
V2 = C1
V3 = C1 + C2
V4 = C1 + C2 – 2.564C4
V5 = C2 + C3 – 0.373C4
V6 = C2 + C3- 4.615CC
V7 = C3 + 2.52C4
V8 = C1 – 0.131C4
Substituting for V’s back into equations 1 through 4, we get
4C1 + 2C2 + 1.431C4 + 4 = 0

4C1 + 4C2 – 7,552C4 + 4.6 =0
2C2 + 4C3 – 2606C4 – 0.4 =0
1.431 C1 –7.552C2 – 2.606C3 + 50.341C4 – 16.577 = 0
Solving the above equation simultaneously,

C1 = - 1.33 C1 = 0.226
C2 = 0.299 C4 = 0.423
Solving for V’s,

V1=0.39 V2 =-1.30 V3 =-1.00’’ V4 =-2.09’’
V5 =0.37’’ V6 =-1.42’’ V7 =1.29’’ V8 =0.17’’

2.8 Triangulation Computation
Computation of lengths
Two sides of each triangle are computed by the law of sines, since one side of the triangle is always known and
the three angles have been measured and adjusted. In computing the sides of the triangles in a quadrilateral, the
solution of the two triangles is sufficient to compute the position of the forward triangulation stations. The two
triangles chosen must be the strongest route through the quadrilateral.
To check the accuracy of the fieldwork or the consistency of the figure after adjustment, computations of the
forward triangulation stations can be repeated through the second-best strong routes. However, only the results
for the strongest route will be used in further computations.
Example:
Compute the lengths of sides of the triangles in order to determine the positions of the forward stations given
that the length of AB is 1192.938m in the example above.
Solution
Here the angles adjusted by the approximate method are used to compute the lengths of the sides of the
triangles.
1. Compute the strength of figures for the four possible combinations of triangles following the method
discussed. The strongest chain is through triangles CBA and DCA (R1 = 4.536), followed by the chain
through triangles DAB and CDB (R2 =5.4).
2. Write the formula used to compute the desired sides. Computation is carried out through the strongest
route (CBA and DCA), but to provide a check on the computations, it is repeated through the second-
strongest route (DAB and CDB).
Side station angle Adjusted angles Length, m
BA 1192.938
C 4 43051’14.67’’
B 2+3 109006’32.40’’
A 1 27002’12.93’’
180000’00.00.
CB=BA(sin 1/sin4) 782.691
AC=BA(sin(2+3)/sin4) 1626.970
AC 1626.970
D 6+7 64024’26.80’’
C 5 56011’57.33’’
A 8 59023’35.87’’

180000’00.00’’
DA=AC (sin5 / sin(6+7)) 1499.049
CD=CA (sin8 /sin(6+7)) 1552.636
BA
1192.938
D 7 39053’02.83’’
A 1+8 86025’48.80’’
B 2 53041’08.37’’
180000’00.00’’
DA=AB(sin2 / sin7) 1499.049
BD=AB(sin(1+8)/sin7) 1856.760
DB 1856.760
0
C 4+5 100 03’12.00’’
D 6 24031’23.97’’
B 3 55025’24.03’’
180000’00.00
BC=DB(sin6/sin(4+5)) 782.691
CD=DB (sin3/ sin(4+5)) 1552.636
Note: Inconsistencies in computed lengths along both chains should not be greater than 0.003m.
Computation of positions
The objective of triangulations is to establish the horizontal positions of the triangulation stations in the network
relative to one another and with respect to the horizontal datum. The network must include at least one point,
and preferably two or more points, the horizontal position of which is known with respect to a national datum,
or co-ordinate system.
If the co-ordinates of two stations A and B are known, the length of the line AB, LAB, calculated from
L AB  ( X B  X A ) 2  YB  Y A 
2
is used as a starting line in the computations of the lengths of the quadrilateral. The azimuth from the north of
the line AB, AABN, is computed from,
XB  XA
AABN  tan 1
YB  YA
Where X A ,X B ,Y A ,and Y B are the X- and Y- co-ordinates of A and B ,respectively.
2.9 Trilateration

A triangulation system can be converted to a pure trilateration system by measuring the lengths of the lines
directly, using the EDM’ s, without measuring any horizontal angles in the network .The net can be adjusted to
obtain the co-ordinates of the stations. However, in order to maintain accuracy of the azimuths of the lines in
the net, astronomical observations are made to selected stations. These measured azimuths impose conditions
that must be satisfied in the adjustment process.
Trilateration is frequently combined with triangulation in order to strengthen a net that may have serious
deficiencies in geometric conditioning, or trilaterated lines are used as check bases in a chain of figures to meet
the criteria of strength of figures. A hybrid net of triangulation and tailateration imposes several rigid geometric
conditions that must be satisfied in the adjustment of the net. In certain instances, the rigidity is necessary for
example; the measurement of very small displacements due to earthquake fault movement must be made over
fairly large areas. Structural deformation or displacement of dams is also made with the highest possible
accuracy. These measurements must be duplicated periodically with the same degree of reliability in order to
reflect the actual condition and to avoid measuring errors, as errors of smaller magnitude may cause a large
destruction. The internal accuracy and reliability of such a network is greatly enhanced by a hybrid
triangulation-trilateration system.
Before lengths of lines in the trilateration net can be used in any subsequent computations, their slope lengths
determined by the EDM instrument and corrected for atmospheric conditions must be reduced to corresponding
sea level distances. Auxiliary data required for determining slope corrections are either the vertical or zenith
angles measured at the two ends of each line or elevations of the two ends of each line. The use of the later will
be discussed in the following section.
Reduction of slope distance by station elevations
When the elevations of the two ends of a trilateration line AB are given, the reduced horizontal distance AB’
can be calculated from,
AB’=(G2-BB’2)1/2-BB’ sin C”:
Where c’ = 16.19” /km (angular error of curvature)
In which G is the slope distance, BB’ is the elevation difference between the two ends of the line AB and AB’ is
the reduced horizontal distance refer to fig below.
G
B
/2
A
B’ C D
M N
.
/2 /2

Fig 2.10 Reduction of slope distance by station elevations

The correction for reduction to sea level, c1, is
c1=AB’ hA/R
Where hA is the elevation of station A and R is the radius of curvature of the earth.
Example:-
Trilateration station A is at an elevation of 576..055m above sea level, and station B is at an elevation of
1693.740m. The slope distance measured between the two stations by EDM is 20, 454.50m. Reduce the
distance to sea level.
Solution;
Given G= 20, 454.50m; hA = 576.055m, hB= 1693.740m.
Then BB= dhAB = 1693.740-576.055 = 1117.685m
The reduced horizontal distance is
AB’=(G2-BB’2)1/2-BB’ sin c”
c”=16.19”/km*20.4545=331.158355”
Substituting the values for the quantities in the equation
AB = 20, 422.15 m
The solution is iterative, and this distance is used to compute c”
C” = 16.19” /km* 20.42215km = 330.6346085”

and computing AB’ using the new value of c”
AB’ = 20,422.149m
The correction for reduction to sea level is,
c1=AB’ hA/R =1.846m
Reduced sea leveldistance=20,420.303m
If length is greater than 20 miles, the chord has to be corrected by
c2=MN3/(24*R2)
To get the sea level length of the line. The result from the above equation is added to the chord MN.
2.10 Adjustment of Trilateration
When extreme accuracy is not required, the angles in each triangle are computed from the lengths of the sides,
and then these angles are adjusted by the approximate methods explained for triangulation. As a check on the
computed angles, the three angles in any one triangle must add up to 1800 exactly whether spherical excess
exists or not. If the spherical angles in a triangle are needed for computation of spherical coordinates, one-third
the spherical excess is added to each computed angle.
Any approximate adjustment of trilateration will produce inconsistencies in subsequent computations.

Therefore, if the work is on a large project with high accuracy, the least squares method is recommended.

2.11 Computation of lengths and positions
If an approximate method of adjustment of the computed angles is used first, the lengths of the sides must be
recomputed to be consistent with the adjusted angles. Before this recalculation can be performed, the strongest
route of triangles in each quadrilateral must be selected, and one side of the triangle must be held fixed and
equal to the reduced sea level distance. If, however, the trilateration begins from a line joining two points with
fixed positions then the fixed length of that line will be used to compute the sides.
The positions of the trilateration stations are computed following the same procedure as that of triangulation
when the approximate method of adjustment is used.
Exercise
In the quadrilateral shown in the accompanying figure, all angles are observed with equal precision with the
following results:
D
C
5 6
4
7 Angle Value
1 28 000, 29.9”
2 390 39’ 44.9”
3 59059” 54.2”
3 8
2 1 4 52019’ 47.0”
B 5. 390 17’30.2”
A
6. 280 22’ 44.0”
7. 530 10’ 43.1”
8. 590 09’ 03.1
Fig 2.11
Adjust the quadrilateral by both the approximate and condition equation methods. Also establish the positions
of the forward stations given that XA = 20.450.768m; YA = 21,000.659m; XB = 21, 450. 762m and YB = 18,
191.470m.

3 Contour Lines and digital terrain model

3.1 Contours and Contour Lines
A contour is an imaginary line that connects points of equal elevations on the ground surface. A line joining
several closely spaced ground points of equal elevation on a drawing is called a contour line. Thus contours on
the ground may be represented be contour lines on the map. On a given map successive contour lines represent
elevations differing by a fixed vertical distance called the contour interval. The choice of this contour interval
depends upon the following factors.
i. The nature of the ground: In a flat and uniformly sloping country, the contour interval is small, but in
mountainous region the contour interval should be large otherwise the contours will come too close to
each other.
ii. The purpose and extent of the survey: Contour interval is small if the area to be surveyed is small and
the maps are required to be used for the design work or for determining the quantities of earth work etc.
Wider interval shall be used for large area and comparatively less important works.
iii. Scale of the maps: The contour interval should be in the reverse ratio to the scale of the map i.e. the
smaller the scale the greater the contour interval.
iv. Time and expense of field and office work: The smaller the contour interval the greater is the amount of
field work and plotting work.
The choice of map scale also depends upon;
i. The clarity with which features can be shown.

ii. The cost: - the larger the scale the higher the cost.
iii. The contour interval
Typical map scales, map uses ad corresponding contour interval are shown in the table below.
Map Scale Typical Use Contour Interval for

mountainous terrain
1:1000 Design 0.25m
1:2000 Design 0.50m
1:5000 Planning 1m
1:10000 Planning 2m
1:25000 Regional Planning 2.5m
1:50000 Regional Planning 5m
1:100000 Regional Planning 5 or 10m
1:250000 Regional Planning 10m
1:500000 State/National 20 r 50m

3.1.1 General Rules for Contour
The principal characteristics of contour lines are as follows:
i. All points on contour lines have same elevation.
ii. Flat ground is indicated where the contour lines are widely separated and steep ground is indicated where
contour lines run close together.
iii. The contour lines connecting points on uniform slopes are spaced uniformly.
iv. Contours crossing a human-made horizontal surface (roads, railroads) will be straight parallel lines as they
cross the facility.
v. A series of closed contour lines on the map represent a hill if the higher values are inside and a depression
if the higher values are outside.
vi. Contour lines across ridge or valley lines are at right angles. If the higher values are inside the bend or
loop, it indicates a ridge and if the higher values are outside the bend, it represents a valley.
vii. Contour line cannot end anywhere but close on themselves either within or outside the limits of the map.
viii. Contour lines cannot merge or cross one another.
3.1.2 Methods of Contouring
There are mainly two methods of contouring:

a. Direct Method
- In this method the contours to be located are directly traced out in the field by locating and
making a number of points on each contour.
- These method is very slow and tedious as a lot of time is wasted in searching points on the
same contour. Compared to the other, it is very accurate and suitable for small areas where
great accuracy is required.
b. Indirect Method
- In this method, the points located and surveyed are not necessary on the contour lines. Spot
elevations will be taken along a series of lines laid out over the area, e.g. grid or cross section.
The contours are then determined by interpolation based on the spot elevations.
- This method is less tedious, cheaper & quicker than the direct method but they are less
accurate when compared to the direct method.

Interpolation of Contour
The process of spacing the contours proportionally between the plotted ground points established by indirect
method is termed as interpolation of contour.
Measuring Slope from Contour
By measuring the distance between two points (∆H) on two contour lines, whose elevations are, of course,
known, the slope of the line connecting these two points can be calculated.
To lay out a highway, railway, canal or any other communication line at a constant gradient, the alignment can
easily plotted on a map.
3.1.3 Uses of Contour Lines (Contour Map)
Enable an engineer to approximately select the most economical and suitable site for an engineering
work such as reservoir, dam, runways, highways etc.
Drawing of a section/profile.
Help in computation of quantities of earthwork.
Determination of catchment areas
Determination of Reservoir Capacity
3.2 Digital Terrain Model (DTM)
3.2.1 Overview
 A digital terrain model represents the relief of the earth’s surface by a set of x, y, z date points where x is
the horizontal x-position [e.g. x longitude, y is the horizontal y-position (e.g. latitude) and z is the
altitude of the earth’s surface [e.g. elevations above see level). It is collection of terrain data as a
sequence of discrete [x, y, z] data points.
 The points are usually horizontally regular spaced in the form of a square grid.
 Just as there are engineering design criteria for selecting a contour interval to represent terrain for a
given application, so to similar criteria are used to select point spacing so that the DTM adequately
represents the terrain. These criteria depend on the potential uses for the data, accuracy requirements, the

terrain character and other factors. The advantages of the regular grid layout a simplified data collection
routine and ease of data access by subsequent programs.
 As the scale dose for maps the grid cell size determines the resolution and degree of generalization of
the DTM
 The grid call size of DTM varies from a few meters (high resolution
DTM, mostly available only for small areas) to medium resolution [20 to 100m) and low resolution
DTM such as the 30 arc sec DTM with 925mm grid cell size which covers the entire world
 Digital terrain model should not be confused with digital elevation model (DEM).
 DEM is only the more general expression for digital surface data but
one must define the kind of surface the elevation date are for :e.g. DEM of the vegetation surface are
for: e.g. DEM of the vegetation surface, DEM of groundwater surface or DEM of the relief of the earth
surface which is also called digital terrain model.
3.2.2 DTM Creation

There are two steps necessary for the creation of a DTM.
1. Collecting original x, y, z coordinates.
2. Interpolations of a regular grid DTM
1. Collecting original x, y, z co-ordinates the altitude data of the earth’s surface can be collected from the
following data source.
I. Digitalization of contour lines from topographic maps.
II. Stereoscopic measurements from aerial photos
III. Stereoscopic measurements from (optical) satellite data
 DTM of regulation surface gives already regular spaced data
IV. Radar satellite data
 gives already regular spaced data
V. Laser scanning measurements
VI. Field measurements or ground surveys
2. Interpolation of a regular grid DTM.
1. The x, y, z co-ordinates collected from the above sources are usually irregularly spaced.

2. The data have to be transformed into regular space grid data using different interpolation
algorithms.
3. The chosen grid cell size usually depends mostly on the density of the collected irregular spaced
data.
4. The quality of the DTM is mainly determined by the density and accuracy of the collected
original altitude data.
3.2.3 Derivation from DTM

i. Calculation of local morphometric values for each grid cell of DTM such as slope angle and
curvature.
ii. Calculation of complex morphometric values for each grid cell of the DTM such as measures of
drainage basins, distance to flow lines and water sheds etc.
iii. Derivation of morphographic features: linear features such as crest lines, flow lines and edges; area
features such as valley grounds summit areas subdividing of slopes etc.
iv. Visualization of the relief, computing contour lines, analytic hill shading and 3d-views.
3.2.4 Application of DTM

 DTM were first used for engineering purposes such as planning of rods, railroads, channels and
water dams or for calculation of mass (soil) movements for bigger hooking areas.
 DTM is a base component of all modern geographical information systems or spatial
information systems.
3.2.5 Uncertainty and Errors in DTMs

1. The accuracy of a DTM depends on the accuracy of its sources data and on the model
resolution.
2. Two DTMs produced from the same data will not contain the same information if their
resolution and sampling strategies are different.
3. A DTM may contain badly formed links or spot height errors which must be corrected by
hand or automatically prior to use.

4 Photogrammetry
4.1 Introduction
Photogrammetry is the science of making measurements on photographs. Terrestrial photogrammetry deals
with photographs taken from a known ground position. Aerial photogrammetry deals with those photographs
taken from the air. Measurements made from photographs can be used to obtain horizontal distances between
points, elevation of points, compilation of topographic and planimetric maps, preparation of profiles and
cross-sections, construction of mosaics, production of orthophotos, etc. Photographs can also be interpreted
for geological, agricultural and engineering investigation, forestry, ecology and environmental analysis,
resource inventory, military intelligence, etc.
4.2 Types of aerial photographs

Photographs taken with the optical axis of the camera held essentially vertical or perpendicular to the ground
surface are called vertical aerial photographs. These have simple geometry and hence they are preferred in
photogranmmetric applications. Some tilt angle (up to 50) is tolerated. Oblique aerial photographs are taken
with the camera axis tilted intentionally. A low oblique is that photograph which bears no image of the
horizon due to small tilt angle. A high oblique shows the image of the horizon with the ground image.
Most planimetric and topographic mapping, mosaic construction and orthophoto production are done by using
vertical aerial photographs although high obloquies are sometimes used in the preparation of small-scale
planimetric maps and charts due to large ground coverage.
4.3 Scale of Aerial Photographs

The scale of a photograph is the ratio of the distance on the photograph to the corresponding distance on the
ground. A vertical aerial photograph resembles a planimetric map in that it shows both the planimetric and
cultural features on the ground in their relative positions. It is different from planimetric maps, however, in
two respects: (1) the photograph does not contain the standard map symbols, which are essential to the
map;(2) the scale of the photograph varies in different portions of the photograph due to the ground relief.
Consider fig 4.1
a b Photo plane a b Photo plane

o o
A B
Ground A O B O
A’ O’ B’ Datum
Datum A’ O’ B’
(a) (b)
Fig 4.1 Scale of photograph on flat and undulating ground
In fig 4.1 (a) points A, O and B all lie at the same elevation. The horizontal distances AO and OB are equal to
A’O’ and O’B’ on the reference datum. A truly vertical photograph taken with the camera at L would show

the images of A, O and B at a, o and b. The ratio ao/ A’O’ equals the ratio ob/O’B’, and the scale of the
photograph is uniform across the photograph.
In fig 4.1 (b) points A, O and B are at different elevations and the horizontal distances between A and O, and
O and B are A’O’ and O’B’ on the datum. If a vertical photograph is taken with the camera axis at L, the
points A, O and B are shown on the photograph as a, o and b. The ratio ao/ A’O’ does not equal the ratio
ob/O’B’, and the scale of the photograph varies across the photograph due to variation in elevations of the
ground.
Referring to fig 4.1(b)
Dis tan ce on photo ao ob

Scale   1 1 1 1
Dis tan ce on ground A O OB
ao f ob f
Scale   ; Scale  
A1O1 H  h A O1 B1 H  hB
f
 Scale, S E  ( 4.1)
H  hE
Where SE is the scale of a vertical photograph for a given elevation E;

f is the focal length of the camera in inches or millimeters;
H is the flying height above datum, in feet or meters; and
h is the elevation of the point, line or area above datum, in feet or meters.
Provided that the relief is not extremely variable, the average scale of a single photograph or a set of
photographs may be desired to be able to measure distances in any are of the photograph(s). It is given by
f
Sav  (4.2)
H  hav
in which Sav is the average scale of the photograph (s) and hav is the average elevation of the area covered by
the photograph, in ft or m, above datum.
When the scale of the photograph is not given it can be determined by measuring a distance between known
landmarks on the photograph and comparing to the corresponding distance on the ground, as for instance the
distance between two read intersections.
The scale of he photograph can also be determined by comparing a distance measured on the photograph with
the corresponding distance measured on a map with a known scale. The photo scale is then found from the
following relation.
photo scale photo dis tan ce
 (4.3)
Map scale Map dis tan ce

Solved Examples:
1.A camera with a focal length of 152.35mm and a picture size 230* 230mm is used to photograph an area
from an altitude of 2400m above sea level. The average ground elevation is 420m above sea level. What is
the average scale of the photography?
Solution
1) From eq. (4.2)
f 152.35mm 1mm
Sav   
H  hav 2400mm  420m 13m
 Sav  1 : 13,000
2) The actual ground area covered by a single photograph
A= 230*13000*230*13000mm2
A= 8,940,100m2
2. Determine the amount of ground area covered by a single photograph of example 1.

Solution
1mm on photo= 13m on ground

1mm2 on photo= 169m2 of area on ground
230*230mm2 of photo area= 230*mm2*169m2/mm2
=8,940,100m2of ground area is covered by a single photograph.
3. Two points lying at an elevation of 268.5m appear on the photograph taken in example1 above. The
distance between these points scales 22.78mm on the photograph. What is the ground distance?
Solution
At the given elevation of 268.5m,
f 152.35mm
S   1 : 13,990
H  h 2400m  268.5m
Dis tan ce on photo 22.78 1

S  
and also Dis tan ce on ground G.D 13.99
 Ground dis tan ce  138 .71m
Exercise: The distance between two road intersections on an aerial photograph is 6.50cm, and the same
distance measures 2.30cm on a map with a scale of 1:50,000. What is the approximate scale of the
photograph?

4.4 Relief Displacement
Relief displacement is the displacement of the image of a ground point on the photograph from the position
the image would have if the point were on the datum, or it is the condition in which a point on the ground is
displaced from its true map position.
In fig 4.2 points B and C are located at elevations of hB and hC above datum. On a vertical photograph their
images appear at the same point. The datum positions B and C would have images at b and c on the
photograph. Point p is the principal focus of the photograph formed by joining opposite fiducial marks at the
comers or the edges of the photograph. Because of the elevations of B and C above datum, both b and c have
been displace outward along the radial lines pb and pc respectively.
The relief displacement dr for a point is give by
rh
dr  (4.4)
H
where r is the radial distance from the principal point to the image point (pb);
h is the elevation of the point above datum; and
H is the flying height above datum, in same units as h.
p, v c b
O
B
E
H
h
Datum
C P, V
Fig 4.2 Relief displacement
From eq.4.4 if can be see that the relief displacement of a point depends on the position of the point on the
photograph. It increases outward toward the edge of the photograph. At the principal point it is zero. Also the
relief displacement increases as the elevation of the point above datum increases and decreases with the flying
height.
Relief displacement can be used to calculate the height of objects in certain situations where the relief
displacement can be measured from the photograph, for instance, for features such as flagpoles, utility poles,
some transmission towers, comers of buildings, water tanks, etc. where both the top and bottom of the object
can be observed on the photograph and the top occurs vertically above the bottom.
Relief displacement can also be used to determine the actual length of lines, angles between lines, areas and
volumes from photographs after reducing all points of the photograph to a common datum.

Solved Example: On an aerial photograph taken at an altitude of 1500m, the relief displacement of the of the
image of a flagpole was measured as 1.6 cm. The distance from the center of the photograph to the image of
the top of the flagpole is 11.0 cm. If the base of the flagpole is at an elevation of 200. 0 m A.O.D, what is the
height of the flagpole?
Solution:
From eq.44
d r  rh / H  1.6cm  11.0cm * h / 1500

h  218.0m
The height of the flag pole above ground is therefore 218.0-200.0=18.0m
4.5 Stereoscopy
Stereoscopy is the ability of the individual to perceive the object space in three dimensions through using both
eyes. Each human eye represents a single camera, and monocular viewing results in flat perspective and the
person’s ability to perceive depth is hampered. Binocular viewing allows the person to view the object from
two different locations due to the separation between the eyes. The pair of camera stations when taking
overlapping photographs is therefore similar to viewing with both eyes.
Viewing with one eye fixes one direction from the eye to the object only, which is in sufficient for fixing the
object’s distance from the viewer. On the other hand viewing with two eyes fixes two directions from both
eyes, hence giving the location of the object at their intersection. The two directions from both eyes form a
convergence angle (parallactic angle) . The closer the point to the eyes, the larger the convergence angle
between the two directions (refer fig 4.3 a). Since1>2, the observer perceives P1as closer than P2, and h is a
function of 1- 2.
Fig 4.3 Stereoscopy

In fig 4.3(b) an idealized tower is photographed from two camera stations. If it were possible to see the left
photograph with the left eye and the right photograph with the right eye only, the observation will perceive the
tower in three dimensions as shown. The lines joining the eyes and the two images of the top of the tower will
intersect at an angle 1, while those passing through the images of the tower base make an angle 2. Since
1>2, the top of the tower will appear closer to the viewer than its base and the tower will then be perceived
in three dimensions.
Stereoscopic Viewing of Aerial Photographs

When aerial photographs are viewed stereoscopically, they produce the three dimensional image of the terrain
in the photograph. The three-dimensional view helps in the study of the terrains as it relates to engineering
projects and also provides the basis for topographic mapping by photogrammetric techniques. There are three
requirements for stereoscopic viewing: (1) the two aerial photographs must provided two views of the terrain
taken from two different camera positions;(2) the two aerial photographs must be oriented properly for
viewing; and (3) the viewer must have normal binocular vision.
1. Photograph overlap: When aerial photographs for project areas are to be taken, flight lines are aid out on a
flight map with a spacing that will cause photographs to cover an overlapping strip of ground both laterally
and consecutively. The lateral overlap between photographs amounts to about 25% of the width of the
photograph of each side. The actual spacing on the flight map may be determined from the scale at which the
photographs are to be taken (eq.4.3). The photographic distance 75% of the width of the photograph if a 25%
overlap is to be maintained between strips (fig 4.4)
Fig. 4.4 Photograph overlap
Along the flight line each photograph overlaps an area that overlaps the area taken by the previous photograph
by about 60%. Therefore, the time interval between exposures is selected to achieve so during the flight. The
large overlap between successive photographs helps to:
(1) Provide coverage of the entire ground area from two view points, such coverage being important for
stereoscopic viewing; (2) allow only the central portion of the photograph to be used for mosaic construction,
thus eliminating to a great extent errors due to relief displacement.
2. Stereoscopic Viewing Devices: The primary purpose of a stereoscopic viewing device is to ensure that the
right eye sees only the photograph on the right and the left eye sees only the photograph on the left. Two
common devices are used: the pocket stereoscope and the mirror stereoscope.

3. Orientation of Photographs for Stereoscopic Viewing: The photographs being viewed stereoscopically are
oriented in the same relative position as the camera existed in the field at the time of photography. The principal
point of each photograph is located by joining opposite fiducial marks at the four comers of each photograph,
and the principal point on the right-hand photograph is identified and marked on the left-hand photograph and
the principal point of the left-hand photograph identified and marked on the right-hand photograph. The
photographs are then overlapped along the flight line with the overlapping images superimposed. Then the two
overlapping images are separated along the flight line with the overlap inwards until a satisfactory fusion is
obtained through the stereoscope. Vision in both eyes is necessary so that the different images of the terrain
contained in the separate photographs can be transferred to the brain by the two eyes of the viewer. The stereo
model of the terrain is then perceived by the brain based on the slight differences in the two views of the terrain
recorded on the photographs.
4.6 Parallax
The algebraic difference of the distances of two images of a ground point from their perspective principal point,
measured parallel to the air base is called parallax.
Referring to fig 4.5 below, a point A at an elevation of h above datum appears at a and a’ on two consecutive
exposures. It can be seen from the figure that the image of A has moved through a total algebraic distance of x-
x’ on the principal plane. This apparent movement of the image of an object between successive exposures is
called parallax of the object.
The parallax p= x-x’ where x and x’ are the coordinates on the left and right photographs respectively, with
respect to the base line. For convenience the position of the image a’ is shown on the left photograph. From the
figure, triangles L1a’’a and AL2L1are similar. Hence,
f p fB
  p (4.5)
H h B H h
Where p= parallax of the object (A), mm
f= the focal length of the lens of the camera, mm
B= the distance between the two camera stations (known as the air base), m.
H= the flying height above datum, m
h=elevation of the object above datum, m.
B
L1 L2
f
a’’ a a’
p x x’
H
A
h
Fig 4.5 Parallax

f b b( H  hav )
We know Sav   B
H  hav B f
Introducing a new quantity Ho=H-hav
bH o
B
f
Substituting for B in equation 4.5
bH o
p (4.6)
H h
Where b= the average distance between the principal points of the overlapping photographs (also called
photo base), mm
Ho= the flying height above mean terrain of the area, m.
Parallax is a direct indication of elevation and can be measured on a pair of overlapping aerial photographs by
means of a parallax bar.
Fig 4.6 Parallax bar
As shown in fig 4.6 above, two measuring marks are located on the bottom surfaces of two pieces of glass
attached to the bar. The left mark is fixed to the bar while the right mark is moved to the left or to the right by
tuning the micrometer screw. A movement to the left is an increase in parallax whereas a movement to the
right is a decrease. If the bar is placed under a mirror stereoscope, the two measuring marks will appear as a
single mark. This is called the floating mark. By turning the micrometer screw one way or the other, this
single mark appears to move up or down according to the increase or decrease in the parallax of the marks.
Fig 4.7 Adjusting photographs for stereoscopic viewing
Suppose that it is required to measure the difference in parallax between points a and b shown on the stereo
pairs in fig. 4.7b. The photographs are first superimposed as shown in fig 4.7a and then separated by a suitable
distance under the stereoscope until satisfactory fusion is obtained (fig 4.7b). The parallax bar is moved over

the photographs to point a, and the micrometer screw is turned until the floating mark appears to touch the
stereoscopic image at point a. The micrometer is read, say as 12.10mm. Then the bar is moved to point b, the
floating mark is again brought to the surface of the stereoscopic image at b. The micrometer is read as, say,
9.65mm. The difference in parallax, dpab, from a to b is thus 9.65-12.10=-2.45mm. In this manner, the
difference in parallax between any two points, such as the top and bottom of a tree or a building, or between
the top and bottom of a mountain can be measured. The distance o1o2is the average photo bade b.
All parallax measurements are made parallel to the base line.
If the parallax of point a is given as pa , the parallax of b will be Pb=pa+dpab.
If only elevation difference between two points, for which no parallax is given, is required, it can be estimated
from
dpH o
dh  (4.7)
dp  b
in which dh is the elevation difference(m), dp and b are in mm, Ho is the flying height above the ground(m).
Coordinates in the ground system can also be obtained in terms of parallax.
From equation 4.1
f
SE 
H h
Substitution in to equation 4.5 gives

p
p  SE B  SE  (4.8)
B
in which p is the parallax of the point in question and B is the air base. If the coordinate in the terrain is
established with the X-axis parallel to the photographic axis of flight x-axis, then
x B
X   x and
SE p
(4.9)
y B
Y y
SE p
In which X and Y are the ground coordinates in the specific system described above and x and y are the
coordinates of the image point in the left photograph with respect to the axis of flight coordinate system. On
the basis of equation 4.9, the distance between two points could be determined if their elevation are known, or
at least the elevation of one point is known in which case the measured parallax is used together with the
image coordinates from equation 4.8. Then, the distance between the two points is given by
D  ( X1  X 2 )2  (Y1  Y2 )2 (4.10)
Solved Examples
1. A rectangular area 130kmx120km is to be mapped from aerial photographs taken to a scale of 1/20000.
The focal length of the lens of the camera to be used is 152mm and each print is to be 230mm square.

Provision is to be made for a 60% overlap between successive exposures and a 25% lateral overlap. Find
(a) the average height above ground at which the plane must operate; (b) the time interval between
exposures in any one strip if the operating speed of the plane is 200km/h; and (c) the minimum number of
photographs required.
Solution
a) From equation 4.2,

f
S av  , but hav  0 at ground level
H  hav
f 152 mm
Hence, H    152 * 20000 * 10 3  3040 m
S av 1 / 20000
b) Let the flight line be parallel to the 130km length. Since there is 60% overlap between successive
exposures, the effective length of each photograph is 40% of 230mm.
0.4*230=92mm
The ground distance covered by this photo length is
92 mm*20000*10-3=1840m.
130,000
Number of photograph per strip=  70.65  71 photos
1840
The operating speed of the plane is 200km/h. To cover the length of 130km, the plane needs
130/200=0.65 hour.
Since the exposures are at regular intervals
0.65
Time interval between exposures= hour  33.12 sec
70.65
c) First find the number of strips required. The width of the area to be photographed is 120km. A 25% lateral
overlap results in an effective photo distance of 0.75*230=172.5mm. The ground distance covered by this
width is 172.5*20000*10-3=3450m.
120,000
Number of strips =  35 strips
3,450
Minimum number of photograph required = 71*35=2485
2. On the overlap of a pair of vertical aerial photographs taken at a height of 2500m AOD with a 152mm
focal length camera are shown two points A and B. Point B is a point on a pass through a range of hills,
while point A is the center of a bridge in a valley. In order to estimate the amount of rise between these
points, parallax measurements were made as follows: Point A: mean reading 5.90mm
Point B: mean reading 11.43mm
The mean level of the valley containing the principal points of the two photographs is 82.00m AOD,
while a BM on the bridge near A was 74.55m AOD. If the respective photo bases are 89.1mm and
91.4mm, calculate the height of B above A.

Solution
From equations 4.5 and 4.6
fB bH o
p 
H h H h
Ho= 2500-82.0=2418 m above mean terrain.
Since the elevation of A can be sufficiently estimated to be 74.55m AOD,
bH o b1  b2 89.1  91.4
Pa  b   90.25 mm
H  ha 2 2
90.25 * 2418
 pa   89.97 mm
2500  74.55
From the parallax bar readings, the difference in parallax is found as
dpab=pb- pa=11.43-5.90=5.53mm.
pb=pa+dpab=89.97+5.53=95.50mm.
bH o bH o 90.25 * 2418
 pb   hb  H   2500   214 .93mAOD
H  hb pb 95.50
Therefore, the height of B above A is 140.40m
3. Given the photo coordinates of A and B as A (40mm, 75mm) and B (-10mm, 80mm) of the left
photograph of example 2 above, determine the ground distance between these points.
Solution
From equation 4.2
f 152mm
Sav    1mm : 15.908m
H  hav 2500  82.0
The air base B= b/Sav= 1435.69m

Then from equation 4.9
B 1435.69
XA  x  40  638.30m
pa 89.97
B 1435.69
Y A y  75  1196.81m
pa 89.97
B 1435.69
XB  x  10  150.33m
pb 95.50
B 1435.69
YB  y  80  1202.67m
pb 95.50
 D  ( X A  X B ) 2  (YA  YB ) 2  788m

4.7 Use and Products of Aerial photographs
Products of aerial photographs

A variety of products are derived from aerial photographs when they are utilized for various operations.
Some of thee products are listed below:
1.Paper- print copies of aerial photographs themselves, which may be used for general planning and for field
operations.
2.Mosaics, which are composed of segments of photographs assembled to give the appearance of a continuous
picture of the terrain. Since each photograph is a perspective, the mosaic contains distortions and should not
be considered as accurate map.
3.Orthophotos specially prepared photographic representations of the terrain without distortions (relief and
tilt) which can be used as a plan metric map. In some cases, contour lines are superimposed on an
orthophoto, resulting in orthophoto maps, which are the equivalent of topographic maps.
4.Planimetric maps, photogrammetrically prepared, contain only the horizontal positions information of
terrain features.
5.Topographic maps, which show both the planimetric detail and contour lines, complied
photogrammetrically by a stereo plotter.
6.Digital data, to produce a system of ground coordinates of a limited number of control points or a dense
network representing the total terrain surface- digital terrain model (DTM)- in the form of the XY- grid
system.
Mosaics
A mosaic is an assembly of overlapping aerial photographs to form one continuos picture of the terrain. In an
uncontrolled mosaic, overlapping photographs are placed down consecutively and stapled to a board. No
horizontal control is used during the assembly. All but the central portion of each photograph is trimmed away,
leaving a small amount of overlap. Use of the central portion of each photograph decreases distortions due to
relief and tilt displacements. A controlled mosaic is constructed to give both high pictorial quality and good
accuracy. This is horizontally controlled by picture points whose horizontal positions on the ground are
established by ground surveys. Each control point has to be identified on the photographic segment to be used.
The accuracy of such a mosaic depends on minimizing the effects of relief and tilt displacements. As for relief,
the smaller the area around the center of each photograph to be used in the mosaic, the less the relief
displacement due to lesser value of r. Therefore, long focal length higher – altitude flying and increased side
and end overleaps are preferred. With regard to tilt the regular aerial photographs are usually replaced by
rectified prints produced by rectifiers. Rectificación produces an equivalent truly vertical photograph from the
same exposure station. The rectification process also allows for bringing all rectified prints to the same scale.
Once the assembly is finished, the mosaic is then photographed with a large copy camera and reproductions are
made at the desired scale.
Orthophotos
An orthophoto is an improvement over a mosaic where the images on the aerial photograph are manipulated
optically mechanically or electronically to remove the perspective aspect and the tilt displacement. The

resulting photograph is an orthophoto. A series of orthophotos are placed together to form orthomosaics to
cover a large area.
Con tour lines can also be plotted during the stereo plotting process of an orthophoto and these can be combined
with the imagery of the area covered. The result is an orthophoto map.
Topographic Information
Topographic information from aerial photographs is obtained by measuring parallax.,
f *B
p
H h
In addition to determining vertical information on a point to basis, it is possible to plot contour lines, in
perspective, using a combination mirror stereoscope and parallax bar. If f, B and H are known in the above
formula, the parallax value p for each contour h (i.e. the plane at elevation h in the terrain) can be calculated.
The corresponding reading on the parallax bar is then set, fixing the floating mark in space where the floating
mark remains in contact with the perceived terrain surface, this motion would describe a contour line. A
drawing pencil attached to the assembly would draw such a contour line. These contours are not orthographic
since they are drawn in relation to one of the photographs, which is perspective. This principle is applied in the
case of large-scale topographic mapping from aerial photographs using stereo plotter. Contour lines produced
photogrammetrically are superior to those produced by ground methods due to the avoidance of interpolation
between spot heights in the former case.
4.8 Applications of Aerial photographs
The most extensive use of aerial photographs is in topographic mapping and orthophotographs. Other
application includes.
 Design, location construction and maintenance of modern highways at all phases, eg. Earthwork
quantities are calculated from photogrammetric models, pavement condition and possible erosion of
embankment can be interpreted from aerial photographs.
 Dam settlement, wave action structural deformation and channel sedimentation can be analyzed from
photographic products.
 In hydrology, analysis of slopes, ground coverage, watershed areas and snow depth can be determined
from photographs these are all important in run off estimation.
 In agriculture, crop inventory and crop diseases analyses are performed by using aerial photographs.
Colored prints are prepared to distinguish among diseases.
 In land surveying photographs are used for identification and location of boundary lines and corners, and
determination of soil and vegetal covers.

5 Introduction to Global Positioning System (GPS) Surveying

5.1 Over of GPS
Global Positioning System (GPS) is a satellite based navigation system that was developed by the U.S.
Department of Defense (DoD) in the early 1970s.Initially, GPS was developed as a military system to fulfill
U.S. military needs. However, it was later made available to civilians, and is now a dual-use system.
GPS provides continuous positioning and timing information, anywhere in the world under any weather
conditions and it is a one-way- ranging system.
GPS consists of 24 operational satellites and some spares. GPS orbits are ellipses with maximum eccentricity is
about 0.01, The six equally spaced orbital planes are inclined at 55° to the equator, resulting in five hours above
the horizon. The system therefore guarantees that at least four satellites will always be in view. The semi major
axis of a GPS orbit is about 26,560 Km .The corresponding GPS orbital period is about 11 hours, 58 minutes.
GPS segments:
GPS consists of three segments-
1. Space segment- consists of 24 satellites, each transmits signal which has number of components. Each
satellite has a fundamental frequency of 10.23 MHz and transmits two L-band radio signals. Signal L1 has
a frequency of 1575.42 MHz (10.23 × 154) and L2 a frequency of 1227.60 MHz (10.23 x 120). The carrier
frequencies and the codes mainly used to determine the distance from the user’s receiver to the GPS
satellites. Navigation message contains along with other information, the coordinates of the satellites as a
function of time.
2. Control segment-The control segment has the task of supervising the satellite timing system, the orbits
and the mechanical condition of the individual satellites. Neither the timing system nor the orbits are
sufficiently stable to be left unchecked for any great period of time. The satellites are currently tracked by
five monitor stations, situated in Kwajalein, Hawaii, Ascension and Diego Garcia, with the master control
in Colorado Springs.
3. User segment – Includes all military and civilian users. With a GPS receiver connected to a GPS antenna,
a user can receive the GPS signals, which can be used to determine his or her position anywhere in the
world. GPS is currently available to all users worldwide at no direct charge.

5.2 GPS Signal
All GPS satellites have atomic clocks and transmit a microwave radio signal derived from such clock,
composed of two carrier frequencies (L1= 1575.42 MHz, L2 = 1227.6 MHz) and modulated by two digital
codes and a navigation message. Even though all of the GPS satellites transmit the same L1 and L2 carrier
frequencies, the code modulation is different for each satellite, which significantly minimizes the signal
interference. Each part of the signal is different from every other, called pseudo-random code. This sequence is
repeated continuously. All GPS receivers know this sequence and repeat it internally. The receiver picks up the
satellite's transmission and compares the incoming signal to its own internal signal. By comparing how much
the satellite signal is lagging, the travel time becomes known. There are three types of code on the carrier
signals, C/A code, P-code and Navigation Message.
The characteristics of these codes and navigation message are shown in the following table.

C/A code P- code Navigation Message

Chipping 1.023MHz 10.23MHz 50 bit per second
rate
Length per 293 m 29.3m 5950 Km
chip
Repetition 1ms 1 week N/A
Code type Gold PRN N/A
Carriers on L1 L1,L2 L1,L2
Code nature Courser code appropriate for 10 times finer that C/A Very course
initially locking onto the code
signal.
Included Time according to the Time according to the Ephemeris, Satellite clock
information satellite clock when the satellite clock when the corrections, Almanac data,
signal was transmitted. signal was transmitted. ionospheric information,
and satellite health data.
Application Moderate Accuracy, Benign High Accuracy, ECM For all the cases
Environment, P-Code Environment, Secure
Acquisition
The GPS signal contains three different types of information: pseudo random code, almanac data and ephemeris data.
The pseudo random code is simply an I.D. code that identifies which satellite is transmitting information.
Almanac data is data that describes the orbital courses of the satellites. Every satellite will broadcast almanac
data for every satellite. Your GPS receiver uses this data to determine which satellites it expects to see in the
local sky. It can then determine which satellites it should track. Almanac data is not precise it can be valid for
many months.
Ephemeris data is data that tells the GPS receiver where each GPS satellite should be at any time throughout
the day. Each satellite will broadcast its own ephemeris data showing the orbital information for that satellite
only.
5.3 Reference coordinate system for GPS

A reference system is the complete conceptual definition of how a coordinate system is formed. It defines the
origin and the orientation of fundamental planes or axes of the system. It also includes the underlying
fundamental mathematical and physical models. A conventional reference system is a reference system where
all models, numerical constants and algorithms are explicitly specified. Time system is also included in
reference systems.
A number of time systems are used in the world. From all of them, UTC and GPS time are the most important
time system for GPS users.
The fact that the topographic surface of the Earth is highly irregular makes it difficult for the geodetic
calculations. To overcome these problem geodesists adopted a smooth mathematical surface, called the
reference surface to approximate the irregular shape of the Earth, Which is sphere for low accuracy positioning,
or ellipsoid for high accuracy positioning. The reference ellipsoid used of GPS work is the WGS84 ellipsoid
(geocentric). With semi major axis (a) = 6378137m and f= 1/298.257223563.

Geodetic datum is a mathematical surface, or a reference ellipsoid, with a well defined origin and orientation.
A geocentric datum is geodetic datum with its origin coinciding with the center of the earth.
A geodetic datum is uniquely determined by specifying eight parameters: two parameters to define the
dimension of the reference ellipsoid; three parameters to define the position of the origin; and three parameters
to define the orientation of the three axes with respect to the earth. In addition to the geodetic datum, the so
called vertical datum is used in practice as a reference surface to which the heights (elevations) of points are
referred. The vertical datum is often selected to be the geoid; the surface that best approximates the mean sea
level on a global basis.
Because the objective of GPS surveys is to locate points on the surface of the earth, it is necessary to have a so
called terrestrial frame of reference, which enables relating points physically to the Earth. Geocentric
coordinate system is a coordinate system in which its the origin coincides with the center of the Earth. Its axis
passes through the Greenwich meridian in the plane of equator. The Z axis coincides with the conventional
Terrestrial pole and Y axis held right hand rule.
Conversion from Satellite reference system to geocentric system four angular parameters is required:
Inclination angle, the argument of perigee, right ascension and Greenwich hour angle.
Of particular importance to GPS users is the 3-Dgeodetic coordinate system. In this system the coordinates of a
point are identified by the geodetic latitude (ϕ), geodetic longitude (λ) and height above the reference surface.

The satellites coordinates as given in the broad cast ephemeris will refer to the WGS 84 reference system.
Therefore, a GPS user who employs the broadcast ephemeris in the adjustment process will obtain his or her
coordinates in WGS 84 system as well.
Ethiopia is using local Adindan datum, with Clark 1880 ellipsoid, so that the coordinates directly found from
GPS have to be transformed into this local datum before using it.
5.4 Fundamentals of GPS positioning
The fundamental technique of GPS is to measure the ranges (distances) between the receiver and a few
simultaneously observed satellites. GPS positioning is based on trilateration, which is the method of
determining position by measuring distances to points at known coordinates. At a minimum, trilateration
requires 3 ranges to 3 known points. But in case of GPS, since there is an offset in clocks of Satellite and
receiver, clock error will be treated as unknown in addition to the coordinates of a point. Such that one
additional measurement is required, i.e. minimum 4 simultaneous “Pseudo ranges” to 4 satellites are necessary.
When a GPS receiver is switched on, it will pick up the GPS signal through the receiver antenna. Once the
Receiver acquires the signal, Time that the signal is transmitted from the satellite is encoded on the signal, using
the time according to an atomic clock onboard the satellite. Time of signal reception is recorded by receiver
using an atomic clock. A receiver measures difference in these times and it will process it using its built-in
software to calculate the range.
Distance = Travel time x Speed of light
The partial outcome of the signal processing consists of the distance to the GPS satellites through the digital
codes (Known as the pseudo ranges) and the satellite coordinates through the navigation messages.
GPS receivers in determining distances to satellites employ two fundamental methods: code ranging or
pseudorange measurement and carrier phase- shift measurements.
Pseudorange measurement is a measure of the range or distance between the GPS receiver and the GPS
satellite which is needed for the position computation. Either the P-code or the C/A code can be used for
measuring th pseudorange. When the PRN code is transmitted from the satellite, the receiver generates an exact
replica of that code. After some time, equivalent to the signal travel time in space, the transmitted code will be
picked up by the receiver. By comparing the transmitted code and its replica the receiver can compute the signal
travel time. Multiplying the ravel time by the speed of light gives the range between the satellite and the
receiver.
The receiver and satellite clocks are not perfectly synchronized and the measured range is contaminated, along
with other errors and biases, by the synchronization error between the satellite and receiver clocks. For this
reason, this quantity is referred to as the pseudorange, not the range.

Carrier- Phase measurement is another way of measuring the ranges. The range would simply be the sum of
the total number of full carrier cycles plus fractional cycles at the receiver and the satellite, multiplied by the
carrier wavelength. The ranges determined with the carriers are far more accurate than those obtained with the
codes. This is due to the fact that the wavelength of the carrier phase, 19 cm in the case of L1 frequency, is
much smaller than those of the codes.
The carriers are just pure sinusoidal waves. This means that all cycles look the same. Therefore, a GPS receiver
has no means to differentiate one cycle from another. It can only measure a fraction of a cycle very accurately,
while the initial number of complete cycles remains unknown or ambiguous. This is therefore, commonly
known as the initial cycle ambiguity, or the ambiguity bias. Fortunately, the receiver has the capability to keep
track of the phase changes after being switched on. This means that the initial cycle ambiguity remains
unchanged over time as long as no signal loss.
It is clear that if the initial cycle ambiguity parameters are resolved accurate range measurements can be
obtained, which lead to accurate position determination.
Cycle Slips is a discontinuity or a jump in the GPS carrier- phase measurements, but an integer number of
cycles, caused by temporary signal loss. This signal loss can be due to obstruction of the GPS satellite signal by
building, bridges, trees, and other objects. This is mainly because the GPS signal is a weak and noisy signal.
Cycle slips may occur briefly or remain for several minutes or even more. Cycle slips must be identified and
corrected to avoid large errors into computed coordinates. This can be done using several methods. Examining
the so called triple difference observable, which is formed by combining the GPS observables in a certain way,
is the most popular in practice.
5.5 Errors in GPS observation
Even though GPS is the most accurate positioning method, it is not free from errors. GPS pseudo code and carrier phase
measurements are subject to different types of errors. These errors can be Satellite dependent, propagation dependent or
receiver dependent.
Satellite dependent errors

GPS carrier and pseudorange code measurements are subjected to errors due to satellite erroneous data transmition to the
receiver.
a) Satellite clock error
b) Satellite ephemeris errors
c) Selective Availability
d) Satellite Geometry
Signal propagation dependent errors

Satellite signal propagates in the atmosphere and this travel media is composed of different layers, namely
ionosphere and troposphere. Due to the presence of these layers the signal propagation is characterised by delay
and bending of the light ray.

a. Ionospheric delay
b. Tropospheric delay
c. Multipath error
Receiver dependent errors

The space segment and the signal propagation media is not the only source of error and biases in GPS
measurement, but also there are errors that depend up on the receiver clocks and the orientation of the antenna.
a. Receiver clock error
b. Antena- phase-center variation
5.6 Linear combinations of GPS observables

GPS measurements are corrupted by a number of errors and biases as discussed in previous section, which are
difficult to model fully. The unmodeled errors and biases limit the positioning accuracy of the stand alone GPS
receiver. Fortunately, GPS receivers in close proximity will share to a high degree of similarity the same errors
and biases. As such, for those receivers, a major part of the GPS errors budget can simply be removed by
combining their GPS observables.
In principle, there are three groups of GPS errors and biases: satellite related errors, receiver related errors, and
atmosphere related errors and biases. The measurements of two GPS receivers simultaneously tracking a
particular satellite contain more or less the same satellite related errors and atmospheric errors. The shorter the
separation between the two receivers, the more similar the errors and biases are. Therefore, if we take the
difference between the measurements collected at the two receivers, the satellite related errors and the
atmospheric related errors will be reduced significantly. In fact, as shown in the next chapter, the satellite clock
error is effectively removed with this linear combination. This linear combination is known as between
receiver single difference.
Similarly, the two measurements of a single receiver tracking two satellites contain same clock errors.
Therefore, taking the difference between these two measurements removes the receiver clock errors. This
difference is known as between satellite single difference.
When two receivers track two satellites simultaneously, two between receiver single difference observables
could be formed. Subtracting these two single difference observables from each other generates the so called
double differences. This linear combination removes the satellite and receiver clock errors. The other errors are
greatly reduced. In addition, this observable preserves the integer nature of the ambiguity parameters. It is
therefore used for precise carrier phase based GPS positioning.

Another important linear combination is known as the triple difference, which results from differencing two
double difference observables over two epochs of time. As explained in the previous sections, the ambiguity
parameters remain constant over time, as long as there are no cycle slips as such, when forming the triple
difference, the constant ambiguity parameters disappear. If, however, there is a cycle slips in the data, it will
affect one triple difference observable only, and therefore will appear as a spike in the triple differences data
series. It is for this reason that triple difference linear combination is used for detecting cycle slips.
All these linear combinations can be formed with a single frequency data, whether it is the carrier phase or the
pseudorange observables. If dual frequency data is available, other useful linear combinations could be formed.
One such linear combination is known as the ionosphere free linear combination. As shown in the next chapter,
ionospheric delay is inversely proportional to the square of the carrier frequency. Based on this characteristic,
the ionosphere-free observable combined to form the so called wide lane observable an artificial signal with an
effective wavelength of about 86cm. this long wavelength helps in resolving the integer ambiguity parameters.
5.7 GPS Positioning Modes
Positioning with GPS can be performed by either of two ways: point positioning or relative positioning.
Point positioning employs one GPS receiver that measures the code pseudoranges to determine the user’s
position instantaneously as long as four or more satellites are visible at the receiver. This positioning method is
used mainly when a relatively low accuracy is required. To determine the receiver’s point positioning at any
time, the satellite coordinate as well as a minimum of four ranges to four satellites are required.
Relative Positioning or differential positioning employs two GPS receivers simultaneously tracking the same
satellites to determine their relative coordinates. Of the two receivers, one is selected as a reference, or base,
which remains stationary at a site with precisely known coordinates. If both receivers track at least four
common satellites, a position accuracy level of the order of a sub centimeter to a few meter can be obtained.
Tracking more than four common satellites would improve the precision of the GPS position solution.

Static GPS surveying is a relative positioning technique that depends on the carrier phase measurements. It
employs two or more receivers simultaneously tacking the same satellites. One receiver, the base receiver, is
setup over a point with precisely known coordinates such as a survey monument. The other receiver, the remote
receiver, is set up over a point whose coordinates are sought. The observation time varies from about 20
minutes to few hours depending on the distance between the base and the remote receivers, the number of
visible satellites, and the satellite geometry.
Fast (Rapid) static surveying is similar to static surveying, but in this case only the base receiver remains
stationary over the known point during the entire observation session. The rover receiver remains stationary
over the unknown point for a short period of time only and then move to another point whose coordinates are
sought. This method is suitable when the survey involves a number of unknown points located in the vicinity of
known points.
Stop- and -Go GPS surveying is another carrier phase based relative positioning technique. In such case, a
base receiver that remains stationary over the known point and one or more rover receivers. The rover receiver
travels between the unknown points, and makes a brief stop at each point to collect the GPS data.
RTK GPS is a carrier phase based relative positioning technique that like previous methods employs two (or
more) receivers simultaneously tracking the same satellites. This method is suitable when:
 The survey involves a large number unknown point located in the vicinity.
 The coordinates of the unknown points are required in real time.
 The line of sight, the propagation path, is relatively unobstructed.
Because of its ease of use as well as its capability to determine the coordinates in real time, this method is the
preferred method by many users. The base receiver remains stationary over the known point and is attached to a
radio transmitter. The rover is normally carried in a backpack and is attached to a radio receiver.
Real time differential GPS is a code based relative positioning technique that employs two or more receivers
simultaneously tracking the same satellites. It is used when a real time meter level accuracy is enough. This
method is based on the fact that the GPS errors in the measured pseudoranges are essentially the same at both
the base and the rover as long as the baseline length is within a few hundred kilometers.

5.8 Other Satellite navigation systems
GPS has never been the only satellite positioning system; other systems such as Transit preceded it, and Russia
has a system called Glonass, which has remained mainly a military system. In addition, the European Union has
planned an entirely independent satellite system called Galileo, which will further enhance navigational
accuracy as well as provide a number of other services (e.g. for search and rescue). In particular, it is intended
that there should be no common mode of failure between GPS and Galileo. In the future the China’s
COMPASS will come.
5.9 Application of GPS

GPS has been available for civil and military use for more than two decades. That period of time has witnessed
the creation of numerous new GPS applications. Because it provides high-accuracy positioning in a cost
effective manner, GPS has found its way into many industrial applications, replacing conventional methods in
most cases.
GPS for civil engineering applications
Civil engineering works are often done in a complex and unfriendly environment, making it difficult for
personnel to operate efficiently. The ability of GPS to provide real-time submeter- and centimeter-level
accuracy in a cost-effective manner has significantly changed the civil engineering industry. Construction firms
are using GPS in many applications such as road construction, Earth moving, and fleet management. In road
construction and Earth moving, GPS, combined with wireless communication and computer systems, is
installed onboard the Earthmoving machine. Designed surface information, in a digital format, is uploaded into
the system. With the help of the computer display and the real-time GPS position information, the operator can
view whether the correct grade has been reached. In situations in which millimeter-level elevation is needed,
GPS can be integrated with rotated beam lasers.
The same technology (i.e., combined GPS, wireless communications, and computers) is also used for
foundation works (e.g., pile positioning) and precise structural placement (e.g., prefabricated bridge sections
and coastal structures). In these applications, the operators are guided through the onboard computer displays,
eliminating the need for conventional methods. GPS is also used to track the location and usage of equipment at
different sites. By sending this information to a central location, GPS enables contractors to deploy their
equipment more efficiently. Moreover, vehicle operators can be efficiently guided to their destinations.

6 Introduction to Geographic Information System (GIS)

6.1 Introduction
Almost everything happens at somewhere. Largely we humans are confined in our activities to the surface and
near surface of the Earth. We travel over it and in the lower levels of the atmosphere, and through tunnel dig
just below the surface. We dig ditches and bury pipelines and cables, construct mines to get at mineral deposit,
and drill wells to access oil and gas. Keeping track of all this activity is important, and knowing where it occurs
can be the most convenient basis for tracking. Knowing where something happens is of critically important if
we want to go there ourselves or send someone there, to find other information about the same place, or to
inform people who live nearby. In addition, most (perhaps all) decisions have geographic consequences.
Therefore, geographic location is an important attribute of activities, policies, strategies, and plans. Geographic
information systems are a special class of information systems that keep the track not only of events, activities,
and things, but also of where these events, activities, and things happen or exist. Because location is so
important, it is an issue in many of the problems society must solve. Some of these are so usual that we almost
fail to notice them- the daily questions of which route to take to and from work, for example. Others are quite
extraordinary occurrence, and require rapid, concerted, and coordinated response by a wide range of individuals
and organizations. Problems that involve an aspect of location, either in the information used to solve them, or
in the solutions themselves, are termed as geographic problems.
6.2 Definitions of GIS

Different types of definition can be give to GIS depends on the profession and the purpose it is using. But
common to all definition is the spatial location. For our purpose we can define GIS as “A GIS is a system,
consisting of hardware, software, data, procedures and proper organizational context, which compiles, stores,
manipulates, analyses, models and visualizes spatial data, to solve planning and management problems”
(Christiansen, 1998). GIS as a system and science can be defined as follow: Geographic Information system
(GIS) is a computer-based system used to store, analyze and display spatial information. Geographic
Information Science (GISc) is a set of principles, procedures and technologies providing scientific foundation to
spatial analysis or spatial science.
6.3 Key elements of GIS

GIS has three key elements: the Geographic, Information and the System. Geographic: It shows the real world
or the spatial realities (Example: the location of a city, location of a school etc) Information: It is about the data
and their meanings (Example: The name of a city, the area of a city, the population density of a city etc)
System: The system is also about the computer technology (Example: Computer hardware &Software).
6.4 History of GIS

The history of GIS is split in to three eras: The era of innovation, era of commercialization, and era of
exploitation.
1) Era of Innovation (1957- 1977): It is the era where GIS is introduced to the world. It was created by
Harvard researchers in the Harvard Laboratory for computer graphics and spatial analysis. The most
important events in the era of innovation were the foundation of ESRI (Environmental System and
Research Institute) and the launch of Landsat I.

2) Era of Commercialization (1981-1999): It is the era where GIS is used to make a business. Hence, a
number of Government and private organizations were established to make GIS a worldwide profit
making industry. The main events of this era were the launch of ArcInfo, introduction of GPS operation
(It is used for navigation, surveying and mapping), and Internet GIS products.
3) Era of Exploitation (1999 -): It is the era where we are now. It is distinct by a high number of GIS users.
The prominent activities of this era are the availability of more than one million users, Launch of
IKONOS and QUICKBIRD satellites, and the introduction of Google earth and Mobile mapping.
6.5 Components of GIS

The people, software, data, hardware and procedures are the five components of a GIS:
Fig.6.1 GIS Components
Hardware: It consists of the computer systems on which the GIS software will run.
Fig.6.2 Hardware component of GIS
Software: It provides tools to manage analyze and effectively display and disseminate spatial data and spatial
information. The main functions of GIS software’s are analytical functions that provide means for driving new
geo-information from the existing spatial and attribute data. The key components of GIS softwares are:
 A database management system (DBMS)
 Tools for the input and manipulation of geographic information
 Tools that support geographic query, analysis, and visualization

 A graphical user interface (GUI) for easy access to tools

GIS software’s can be divided as public domain and commercially available software’s. Example: ArcView,
ArcGIS, ENVI, ERDAS, IDRISI, ILWIS, and many more others.
Data: It is the most important component of a GIS. Geographic data and related tabular data can be collected in
house, compiled to custom specifications and requirements, or purchased from a commercial data provider. A
GIS can integrate spatial data with other existing data resources, often stored in a DBMS. The integration of
spatial and tabular data stored in a DBMS is a key functionality afforded by GIS.
People: GIS technology has limited value without the people who manage and develop plans for applying it to
real world problems. GIS user ranges from technical specialists, who design and maintain the system, to
customers who use it to help them perform their everyday work. The identification of GIS specialist's vs. end
users is often critical to the proper implementation of GIS technology.
Method: Method include how the data will be retrieved, input into the system, stored, managed, transformed,
analyzed, and finally presented in a final output. The procedures are the steps taken to answer the question
needed to be resolved. The ability of a GIS to perform spatial analysis and answer these questions is what
differentiates this type of system from any other information systems. A successful GIS operates according to a
well designed implementation plan and business rules, which are the models and operating practices unique to
each organization.
6.6 Data Model in GIS

The heart of any GIS is the data model, which is a set of constructs for description and representation of objects
and processes in the digital environment of the computer. People (GIS users) interact with operational GIS in
order to undertake tasks like making maps, querying databases, and performing analysis and presentation of out
puts. Because the types of analyses that can be undertaken are strongly influenced by the way the real-world is
modelled, decisions about the type of data model to be adopted are vital to the success of a GIS project.
Generally there are three types of data which are involved in GIS based analysis and mapping, namely spatial
data, attribute data and Meta data.
Spatial data is also known as geospatial (Coordinate) data or geographic information/data. It is the data or
information that identifies the geographic location of features and boundaries on earth, such as natural or
constructed features, parcels, roads, buildings.
Attribute data is non geographic, tabular and descriptive data that GIS links to spatial data. Attribute data is
collected and compiled for specific areas like states, cities, and so on and often comes packaged with map data.
Metadata is data about the data itself, which includes, the spatial reference system attached to the spatial data,
accuracy, description about every attribute fields, and so on.
Spatial data models begin with a conceptualization and also there are rules which govern the view of real world
phenomena or entities. GIS store information about the world as a collection of thematic layers that can be
linked together by geography. This simple but extremely powerful and versatile concept has proven invaluable
for solving many real-world problems from tracking delivery vehicles, to recording details of planning
applications, to modeling global atmospheric circulation. The thematic layer approach allows us to organize the

complexity of the real world into a simple representation to help facilitate our understanding of natural
relationships.
Consider a road map suitable for use at a statewide or provincial level. This map is based on a conceptualization
that defines roads as lines. These lines connect cities and towns that are shown as discrete points or polygons on
the map. Road properties may include only the road type. Examples may include a limited access interstate,
state highway, country road, or some other type of road. The roads have a width represented by the drawing
symbol on the map. However this width, when scaled, may not represent the true road width. This
conceptualization identifies each road as a linear feature that fits into a small number of categories. All state
highways are represented by the same type of line, even though the state highways may vary. Some may be
paved with concrete, others with bitumen. Some may have wide shoulders, others not, or dividing barriers of
concrete, versus a broad vegetated medium. We realize these differences can exist within this conceptualization.
Therefore, GIS represents real world objects (buildings, roads, land use, elevation) with digital data. Real world
objects can be divided into two abstractions: discrete objects (a house) and continuous fields (elevation).
Two methods are used to reduce geographic phenomena to forms that can be coded in computer databases, and
we call these raster and vector.
Fig.6.3 Thematic Laysers
In principle, both can be used to code both fields and discrete objects, but in practice there is a strong
association between raster and fields, and between vector and discrete objects.
6.6.1 Vector model

The first conceptualization defines discrete entities that may be represented by discrete objects, a vector data
model. A farm field, road, wetland, cities and census tracts are examples of discrete entities that may be
represented by discrete objects. A vector data model uses coordinates to store the shape of a spatial entity and
associated attribute data to define discrete objects.

In vector data models, there are there vector object used to represent the geometry of real world entities, namely
point, line and polygon. Groups of coordinates define the location and boundaries of these objects, and these
coordinate data plus their associated attributes are used to represent the entities. In the vector world, the point is
the building block from which all spatial entities are constructed. The smallest spatial entity, the point, is
represented by a single (x, y) coordinate pair.
The vector representation using point, line or polygon depends up on the map scale and the extent of the work,
since a city can be represented by a point in small scale map but polygon in large scale map which covers small
area.
Fig.6.4 Relationships between points,lines and polygons
6.6.2 Raster model

Raster data models define the world as a regular set of cells in a grid pattern. The grid based square or
rectangular cell is called pixel is organized in rows(from top to bottom) and columns(from left to right ) where
in each cell a single value of attribute which represents the phenomena or entities of interest is stored. These
cells are usually evenly spaced in the x and y directions and the organizations are quite easy to code in an array
structure in most computer languages.
Thus, individual cell is used as the building block for creating images of points, lines and areas in the raster
world.
These might also represent photographic or scanned images. A raster cell stores a single value, it can be
extended by using raster bands to represent RGB (Red, Green, Blue) colors, color maps (a mapping between a
thematic code and RGB value), or an extended attribute table with one row for each unique cell value. A point
is indicated with a single cell, a line by several cells with the same value forming a linear grouping, and an area
by a clump of cells all having the same value.
Raster data models are the natural means to represent “continuous” spatial features or phenomena. Elevation,
precipitation, slope, and pollutant concentration are examples of continuous spatial variables. These variables
characteristically show significant changes in value over broad areas. The gradients can be quite steep (e.g., at
cliffs), gentle (long, sloping ridges), or quite variable (rolling hills). Because raster data may be a dense
sampling of points in two dimensions, they easily represent all variations in the changing surface. Raster data
models depict these gradients by changes in the values associated with each cell.

Raster data sets have a cell dimension, defining the size of the cell. The resolution of the raster dataset is its cell
width in ground units. For example, in a Landsat TM raster image, each cell may be a pixel that represents an
area of 30 meters by 30 meters.
The figure below shows representation of the given reality both in raster and vector data models.
Fig.6.5 Raster and Vector representation of real world
6.6.3 Comparison of Vector and Raster

The question often arises, “which are better, raster or vector data models?” the answer is neither of both.
Neither of the two classes of data models is better in all conditions or for all data. Both have advantages and
disadvantages relative to each other and to additional, more complex data models. As an example, elevation
may be represented as sets of contour lines in a vector data model or as a set of elevations in a raster grid. The
choice often depends on a number of factors including the predominant type of data (discrete or continuous), the
expected type of analysis, available storage; the main sources input data or the expertise of human operators.
Raster data models exhibit several advantages relative to vector data models. First raster data models are
particularly suitable for representing themes or phenomena that change frequently in space. Each raster cell may
contain a value different than its neighbors. Thus, trends as well as more rapid variability may be represented.
Raster data structures are generally simpler, particularly when a fixed cell size is used.
Raster data models also facilitate easy overlays, at least relative to vector models. Each raster cell in a layer
occupies a given position corresponding to a give location on the earth surface. Data in different locations align
cell-to-cell over this position. Thus overlay involves locating the desired grid cell in each data layer comparing
the values found for the given cell location. This operation is quite rapid in raster data structure.
Finally raster data structures are the most practical method for storing, displaying, and manipulating digital
image data such as aerial photographs and satellite imagery.
Vector data models provide some advantages relative to raster data models. First, vector models generally lead
to more compact data storage, particularly for discrete objects. Large homogenous regions are recorded by the
coordinate boundaries in a vector data model, the same regions are recorded as a set of cells in case of raster.
Vector data are more natural means for representing networks and other connected features. Vector data by their
nature store information on intersections (nodes) and the linkage between them (lines).

Vector data models are easily presented in a preferred map format. Humans are familiar with continuous line
and rounded curve representations in hand or machine drawn maps and vector based maps show these curves
whereas raster data often show a “stair step” edge for curved boundaries, particularly when the cell resolution is
large relative to the resolution at which the raster is displayed. Vector data models facilitate the calculation and
storage of topological information. Topological information aids in performing adjacency, connectivity and
other analysis in an efficient manner. Topological information also allows some forms of automated error and
ambiguity detection, leading to improved data quality.
6.6.4 Triangulated Irregular Network (TIN)

TIN is a data model commonly used to represent terrain heights. It store GIS data for 3D surface model.
Typically the x, y and z locations for measured points are entered into the TIN model. These points are
distributed in space and the points may be connected in such a manner that the smallest triangle formed from
any three points may be constructed. The TIN forms the connected network of triangles and therefore the basic
unit is a triangle. Because a triangle consists of three lines connecting three nodes, each triangle will have three
neighbours (except those on the side or periphery). The triangle is represented by a sequence of three nodes.
Each triangle may have other associated attributes such as population density, crime rate, etc. in another table.
The TIN model usually uses some forms of indexing to connect neighbouring points. Each edge of a triangle
connects to two points, which in turn each connect to other edges. These connections continue until the entire
network is spanned. Thus, the TIN is a rather more complicated data model than the simplest raster grid when
the objective is terrain representation
Fig.6.6 Triangulated irregular network data model
6.7 Data Capturing and Entry techniques in GIS

The first step of using a GIS is the provision of data. Much of the success of a GIS project depends on the
quality of the data that is entered into the system, and thus this phase of a GIS project is critical and must be
taken seriously. Data collection is one of the most time-consuming and expensive yet important, of GIS tasks.
Two main types of data capture are:

 Primary data sources are those collected in digital format specifically for use in a GIS project. Raster
data capture includes remote sensing; aerial photograph, DEM, The advantage of such a GIS data
capture is that there is a consistency in the data generated, and the whole process can be regularly and
systematically assorted to get accuracy of the data at a very cost effective manner. Vector data capture
refers ground surveying and GPS. Ground survey is a very time-consuming and expensive activity, but it
is still the best way to obtain highly accurate point locations.

 Secondary data sources are digital and analog datasets that were originally captured for another
purpose and need to be converted into a suitable digital format for use in a GIS project. Paper maps,
aerial photographs, and images are scanned (converting analogue map into a computer readable form)
and georeferenced so that they provide geographic context for raster data source. Secondary vector data
capture involves digitizing (converting raster data into vector data and the process called as
Vectorization) vector objects from maps and other geographic data sources. Entry of coordinates using
Coordinate geometry and conversion existing digital data are also another data capturing techniques
Data transfer: It involves importing digital data from other sources (E.g.: Internet, WAN, LAN and from
physical media CD ROMS).
Attribute data can be entered by direct data loggers or manual keyboard entry. Data has to be entered and stored
in a proper digital format in the computer to apply any GIS analysis.
6.8 Application areas of GIS

GIS can answer number of questions as long as it is based on spatial location. Condition can be accessed by
answering what is....? Type of question, Where is....? Question is addressed by specifying the required event.
Trends and patterns are also studied by analyzing temporal data.
Generally GIS is applicable in number of fields since everything is related to spatial location.
 Water Resource:- GIS provides calculations and methods for watershed characteristics, flow statistics,
flood zoning and forecasting, flow direction, ground and surface water management, flow accumulation
and facilitates the watershed delineation by using Digital Elevation Models (DEMs)
 Agriculture: - GIS is used in a variety of such as managing crop yields, monitoring crop rotation
techniques, and projecting soil loss for individual farms or entire agricultural regions.
 Business: -GIS is a tool for managing business information by using the location. Business sites can be
selected, sale territories can be optimized and customers location can be tracked
 Environmental and natural resource management: -GIS can be used to produce maps, inventory
species, measure environmental impact, or trace pollutants, study geological features, analyze soil,
access seismic information.
 Engineering: - GIS is a tool for traffic analysis and planning, topographic map preparation, location,
decision and suitability analysis, pavement management, route selection, landuse planning, urban
planning, infrastructural mapping.
 Risk management: - GIS can help in risk management and analysis by showing which area will be
subjected to natural and manmade disaster and prevention measures can be developed. e.g. soil erosion,
flood, draught and diseases
Thematic cartography: - GIS is used to show different even distributions in forms of maps, tables
and charts.

Surveying II AAiT

Uploaded by

Copyright:

Available Formats

Surveying II AAiT

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Surveying II AAiT

Uploaded by

Copyright:

Available Formats

ammanuel bekele

DEPARTMENT OF CIVIL ENGINEERING, AAIT 1 of 74

4.3 Scale of Aerial Photographs ..................................................................................................... 44

1.2 Types of curves and their use

Figure 1.1 Types of Circular Curve

Figure 1.2 Use of spiral transition curves

1.3 Circular Curve

Figure 1.3 Elements of circular curve

Degree of the circular curve

Figure 1.4 Degree of a curve using arc definition

1.4 Setting Out of Circular Curves

1. Methods Using Tape

Fig 1.5 Offset from long chord

Fig 1.6 Offset from tangent

3. Setting out using two theodolite

Fig. 1.8 Setting out by two theodolite method

4. Setting out using Coordinate

where AzVA is the back azimuth of line AV.

AzAO = AzAV + 90°

Fig 1.9 Setting out by coordinate method

5. Setting out with inaccessible intersection point

Fig 1.10 Setting out with inaccessible intersection point

1.5 Transition Curves

Fig 1.11 Super elevation

Derivation of Transition Curve Equations

Fig 1.12 Derivation of transition curve

Relation Between  and 

=  ( R  R cos1 )    R  R(1    ...)  R

Fig 1.13 Derivation of shift

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

1.6 Setting Out of Transition Curves from the Tangent Point

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

1.7 Vertical Curves

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

Fig 1.14 Vertical curve: summit and sag

where p and q are in %

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

The properties of the parabola give: -

From Fig (c)

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

2.3 Classification of Triangulation Systems

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

2.4 Triangulation figures and choices

(a) Chain of triangles (b) chain of quadrilaterals (c) chain of polygons.

Fig 2.1 triangulation figures

1. Reconnaissance, select the locations of stations;

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

2.5 Angle and side conditions

CA=total number of angle conditions (including center-point equations)

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

Fig 2.3 (a) Inconsistent Sides (b) Side Condition

The number of extra lines or the number of side conditions is

2.6 Strength of Figures

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

Triangle Known Computed Distance Angles  A2   a b   b 2

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

Lecture Note: - CENG 2803: SURVEYING II 2012/13 AY

Baseline and Angle Measurement