Photogrammetry Handout
Photogrammetry Handout
Principal Points A, B, C,
and D, which lie at
elevations above datum
of hA, hB, hC, and, hD
respectively, are imaged
on the photograph at a, b,
c, and d.
The scale at any point
can be expressed in terms
of its elevation, the
camera focal length, and
the flying height above
datum.
From the given Figure, from similar triangles Lab and LAB, the following expression can be
written:
It is apparent from these relationships that photo scale increases at higher elevations
and decreases at lower ones.
This concept is seen graphically in the above Figure, where ground lengths AB and
CD are equal, but photo distances ab and cd are not, cd being longer and at larger
scale than ab because of the higher elevation of CD.
In general, by dropping subscripts, the scale S at any point whose elevation above datum is
h may be expressed as
Here S is the scale at any point on a vertical photo, f is the camera focal length, H the flying
height above datum, and h the elevation of the point.
Use of an average photographic scale is frequently desirable, but must be
accepted with caution as an approximation.
For any vertical photographs taken of terrain whose average elevation above datum is the
average scale is
Example
1.Suppose that highest terrain h1, average terrain havg, and lowest terrain h2 are
610, 460, and 310 m above mean sea level, respectively. Calculate the maximum
scale, minimum scale, and average scale if the flying height above mean sea level is
m and the camera focal length is 152.4 mm.
Solution by maximum scale occurs at maximum elevation
Equation (27.8) can be used to locate the datum photographic positions of images on a vertical photograph.
can also be applied in computing heights of vertical objects such as buildings, church steeples, radio
towers, trees, and power poles. To determine heights using the equation, images of both the top and
bottom of an object must be visible.
Example. In Figure the above figure, radial distance to the image of the
base of the pole is 75.23 mm, and radial distance to the image of its top is
76.45 mm. The flying height H is 4000 ft above datum, and the elevation
of B is 450 ft. What is the height of the pole?
Where AB is the horizontal length of ground line AB, H the flying height above
datum, and the elevations of the control points above datum, and x and y
the measured photo coordinates of the control points.
In the above Equation all variables except H are known. Hence a direct solution
can be found for the unknown flying height.
The equation is quadratic, so there are two solutions, but the incorrect one will be obvious and
can be discarded.
Error Evaluation
Some of the more significant sources of errors are:
I. Errors in photographic measurements, e.g., line lengths or photo
coordinates
II. Errors in ground control
Sources I&II can be minimized precise, properly calibrated equipment and
suitable caution are used in making the measurements.
III. Shrinkage and expansion of film and paper
Source III can be practically eliminated by making corrections.
IV. Tilted photographs where vertical photographs were assumed
Magnitudes of error introduced by source IV depend upon the severity of
the tilt.
Example: Assume that a vertical photograph was taken with a camera having a focal
length of 152.4 mm. Assume also that a ground distance AB on flat terrain has a
length of 1524 m and that its corresponding photo distance ab measures 127.0
mm. Flying height above ground may be calculated, as follow:
Interior Orientation
Three parameters plus distortion
principle point H, (x0,Y0)
Principal distance C (C=f (focal length)
Geometric Definition
Principal points: intersection point of optical axis through image plane
Principal distance: distance of projection center from image plane
Calibration
H, C are calibrated in laboratory and are assumed to be known for image measurements
Alternatively in-situ calibration (requires GCPs)
Exterior Orientation
Defines the orientation (position and attitude) of the image, i.e. the image coordinate
system in relation to an exterior-(world-, object-) coordinate system.
2D image coordinate system
Defined by interior orientation
In image plane (2D- Cartesian)
Origin in principal point
3D image coordinate system
Defined by interior orientation
3D- Cartesian
Origin in projection center
3D/2D coordinate systems are fixed to camera/image plane
Suppose that a pair of aerial photography is taken from exposure stations L1 and L2 so
that the building appears on both photos.
Flying height above ground is H′, and the distance between the two exposures is B, the
air base.
Object points A and B at the top and bottom of the building are imaged at a1 and b1 on
the left photo and at a2 and b2 on the right photo.
placing photos laid on table so left eye sees left image and right eye sees the right image
Stereoscopic model (called stereo model) appears below table top, overlapping pair
photographs is stereopair
Brain judges height of building from differences in parallactic angle
4.3.1Stereoscopes
It is quite difficult to view photographs stereoscopically without the aid of optical devices. These
difficulties in stereoscopic viewing overcome by use of stereoscopes instruments. Stereoscopy is
a three dimensional perception obtained with the aid of a stereoscope equipment. Enables you to
view an object from 2 different camera positions to obtain a 3-dimensional view.
In order to get a stereoscopic vision of a certain area, most photogrammetrists use either:
i. Lens stereoscopes/pocket stereoscope
ii. Mirror stereoscopes
I. Pocket stereoscope
It is one of the most popular lens stereoscope commonly known as pocket stereoscope.It has plane-
convex lens, upper side flat with a focal length of 100mm.The rays entering the eyes are parallel
and converge at infinity and have been accommodated (focused) at 100mm distance .The pocket
stereoscope is cheep, transportable, and can be used in the field. It has 2-4 x lens magnification .
Disadvantages of pocket stereoscope
I. Limited magnification
II. The distance between corresponding points on the photos must be equal to or smaller
than the eye base. With normal size photographs this becomes difficult or impossible
without bending or folding the photos.
Page | 14
Mirror stereoscope LAS DU/2023
Photogrammetry hand out
Use of stereoscopes
In stereoscopic viewing, it is important to orient the photos so that the left and right eyes see the
left and right photos, respectively.
If the photos are viewed in reverse, a pseudoscopic view results in which ups and downs
are reversed; e.g., valleys appear as ridges and hills appear as depressions.
Accurate and comfortable stereoscopic viewing requires that the eye base
The line joining the centers of the stereoscope lenses, and the flight line all be parallel.
For vertical photographs, the flight line is the line from the center of the left photo to the
center of the right photo.
In below figure the left photo was exposed from a lower flying height than the right photo,
and consequently its scale is larger than the scale of the right photo. Even though the photos
are truly vertical and properly oriented, y parallax exists at both points a and b due to
variation in flying heights.
In below figure the left photo was exposed from a lower flying height than the right photo, and
consequently its scale is larger than the scale of the right photo. Even though the photos are
truly vertical and properly oriented, y parallax exists at both points a and b due to variation in
flying heights.
The effect of tilted photos is illustrated in the following figure. The left photo is truly vertical and
shows positions of images a through d of a square parcel of property on flat terrainThe right photo
was tilted such that the same parcel appears as a trapezoid. In this case, y parallax exists
throughout the stereoscopic model as a result of the tilt. Most serious y parallaxes usually occur
from improper orientation of the photos, a condition which can be easily corrected.
Photographic scale
Altitude above terrain
Air base terrain relief
Figures a and b represent, respectively, the taking of a pair of vertical overlapping photographs
and the stereoscopic viewing of those photos. In figure a, the camera focal length is f, the air base
is B, the flying height above ground is H′, the height of ground object AC is Z, and the horizontal
ground distance KC is D. In Figure a, assume that Z is equal to D. In Fig. b, i is the image distance
from the eyes to the photos, be is the eye base, h is the distance from the eyes to the perceived
stereomodel, z is the stereomodel height of object A′C′, and d is the horizontal stereomodel
distance K′C′. Note that while the ratio Z/D is equal to 1, the ratio z/d is greater than 1 due to
vertical exaggeration. An equation for calculating vertical exaggeration can be developed with
reference to these figures. From similar triangles of Figure a.
Example1: Estimate the stereoviewing base-height ratio if the height of the stereoscope above the photos is 10.00
cm, and the photos are placed 5.0 cm apart.
Example 2 :Using the results of Example 1, calculate the approximate vertical exaggeration for vertical aerial
photos taken with a 152.4-mm-focal-length camera having a 23-cm-square format if the photos were taken with 60
percent end lap.
Note: If a 305-mm-focal-length camera had been used, the B/H′ ratio would have been 0.30, and vertical
exaggeration would have been reduced to 2.
The parallax of any point is directly related to the elevation of the point.
Parallax is greater for high points than for low points.
Shows the two photographs of in superposition. Parallaxes of object points A and B are pa and
pb, respectively. Stereoscopic parallax for any point such as A whose images appear on two
photos of a stereopair, expressed in terms of flight-line photographic coordinates,
By equating similar triangles of the above fig , formulas for calculating hA, XA, and YA may be
derived. From similar triangles L1oay and L1AoAy,
Now substituting the above into each of Eqs. (b) and (a) and reducing gives
Where as hA is the elevation of point A above datum, H is the flying height above datum, B is the
air base, f is the focal length of the camera, pa is the parallax of point A, XA and YA are the ground
coordinates of point A in the previously defined unique arbitrary coordinate system, and xa and ya
the photo coordinates of point a measured with respect to the flight-line axes on the left photo.
Example .A pair of overlapping vertical photographs was taken from a flying height of 1233 m
above sea level with a 152.4-mm-focal-length camera. The air base was 390 m. With the photos
properly oriented, flight-line coordinates for points a and b were measured as xa = 53.4 mm, ya =
50.8 mm, xa′ = –38.3 mm, ya′ = 50.9 mm, xb = 88.9 mm, yb = –46.7 mm, xb′ = –7.1 mm, yb′ = –
46.7 mm. Calculate the elevations of points A and B and the horizontal length of line AB. Solution
Example: flight-line axis x and x′ coordinates for the images of a vertical control point C were
measured as xc = 14.3 mm and xc′ = –78.3 mm. If the elevation of point C is 591 m above sea
level, calculate the elevations of points A and B of that example, using parallax difference
Other errors of lesser consequence such as camera lens distortion and atmospheric
refraction distortion
Chapter: Five
Analytical Photogrammetry: part one
5.1 Introduction
Analytical photogrammetry is a term used to describe the rigorous mathematical calculation of
coordinates of points in object space based upon camera parameters, measured photo
coordinates, and ground control. Analytical photogrammetry generally involves the solution of
large, complex systems of redundant equations by the method of least squares.The evolution of
computer technology and photogrammetric software, however, made analytical photogrammetry
a commonplace technique.
Analytical photogrammetry forms the basis of many modern hardware and software
systems, including:
stereoplotters (analytical and softcopy),
digital terrain model generation,
orthophoto production,
Whereas x′a and y′a are corrected photo coordinates and xa and ya are measured coordinates. The
ratios xc/xm and yc/ym are simply scale factors in the x and y directions, respectively.
Example .For a particular photograph, the measured x and y fiducial distances were 233.8 and
233.5 mm, respectively. The corresponding x and y calibrated fiducial distances were 232.604 and
232.621 mm, respectively. Compute the corrected values for the measured photo coordinates
which are listed in columns (b) and (c) in the table below.
The radial distortion value was the radial displacement from the ideal location to the actual image
of the collimator cross, with positive values indicating outward displacements.
The form of the polynomial, based on lens design theory, is
In above Eq, Δr is the amount of radial lens distortion, r is the radial distance from the principal
point, and k1, k2, k3, and k4 are coefficients of the polynomial. The coefficients of the polynomial
are solved by least squares using the distortion values from the calibration report. To correct the x,
y position of an image point, the distance r from the image point to the principal point is computed
and used to compute the value of Δr. This is done by first converting the fiducial coordinates x and
y, to coordinates and , relative to the principal point & r can be compute with
Example: An older USGS camera calibration report specifies the calibrated focal length f =
153.206 mm and coordinates of the calibrated principal point as xp = 0.008 mm and yp = –0.001
mm. The report also lists mean radial lens distortion values given in columns (a) and (b) of the
table below. Using these calibration values, compute the corrected coordinates for an image point
having coordinates x =62.579 mm, y = –80.916 mm relative to the fiducial axes.
Solution Compute r values (in meters) in column (c) by the following equation
Using the least squares method the following k values were computed.
Compute the distance from the principal point to the image point
that light rays travel in straight paths, and to compensate for the known refracted paths, corrections
are applied to the image coordinates.
In this equation, α is the angle between the vertical and the ray of light,
and K is a value which depends upon the flying height above mean sea level and the elevation of
the object point. A convenient method, adapted from the Manual of Photogrammetry, is to compute
K by
Example: A vertical photograph taken from a flying height of 3500 m above mean sea level
contains the image a of object point A at coordinates (with respect to the fiducial system) xa =
73.287 mm and ya = –101.307 mm. If the elevation of object point A is 120 m above mean sea
level and the camera had a focal length of 153.099 mm, compute the x′ and y′ coordinates of the
point, corrected for atmospheric refraction.
Solution Compute r
Generally, processing image for the purpose of ortho-rectification, the minimum number
of GCPs required is three.
In establishing the mathematical relationship between ground and image space, some
parameters must be determined. The parameters include a scale factor (describing the scale
difference between image space and ground space); X, Y, Z coordinates (defining the
positional differences between image space and ground space); and three rotation angles
(omega, phi, and kappa) that define the rotational relationship between image space and ground
space. It is highly recommended that a greater number of GCPs be available than are actually
used in the block triangulation. Additional GCPs can be used as check points to independently
verify the overall quality and accuracy of the block triangulation solution. The result of the
analysis is an RMSE that defines the degree of correspondence between the computed values
and the original values. Lower RMSE values indicate better results
Tie Points
A tie point is a point that has unknown ground coordinates, but is visually recognizable in the
overlap area between two or more images and computed during block triangulation. It can be
measured both manually and automatically and should be visually well-defined in all images.
Tie points should also be well distributed over the area of the block and typically, nine tie points
in each image are sufficient for block triangulation. The following figure shows the placement of
tie points.
In a block of images with 60% overlap and 25-30% side lap, nine points are sufficient to tie
together the block as well as individual strips as shown in the following figure
where 𝑥� and 𝑦� are the photo coordinates of image point 𝑎�: 𝑋�𝐴�, 𝑌�𝐴�, and 𝑍�𝐴� are object space
coordinates of point A. 𝑋�𝐿�, 𝑌�𝐿� , and 𝑍�𝐿� are object space coordinates of the exposure station 𝑓�
is the camera focal length 𝑥�0 and 𝑦�0 are the coordinates of the principal point (usually known
from camera calibration) the 𝑚�′𝑠� are functions of three rotation angles, and most often
omega, phi, and kappa are the angles employed .
In the above Eq. Subscripts 1 and 2 affixed to terms D, E, and F indicate that the terms
apply to either photo 1 or photo 2. The m’s again are functions of the three rotation angles
omega, phi, and kappa. One coplanarity equation may be written for each object point
whose images appear on both photos of the stereopair. The coplanarity equations do not
contain object space coordinates as unknowns; rather, they contain only the elements of
exterior orientation of the two photos of the stereopair.
Like collinearity equations, the coplanarity equation is nonlinear and must be linearized by using
Taylor’s theorem and solved iteratively for corrections to approximations of the orientation
parameters.
In Eqs. (11-6) and (11-7), the terms are as defined in Secs. 11-4 and D-5. Two equations are formed
for each control point, which gives six equations if the minimum of three control points is used. In
this case a unique solution results for the six unknowns, and the residual terms on the right sides
of Eqs. (11-6) and (11-7) will be zero. If four or more control points are used, more than six
equations can be formed, allowing a least squares solution. Since the collinearity equations are
nonlinear, and have been linearized using Taylor’s theorem, initial approximations are required
for the unknown orientation parameters.
The linearized forms of the space intersection equations for point A are
The process of forming an analytical stereomodel involves three primary steps: interior orientation,
relative orientation, and absolute orientation. After these three steps are achieved, points in the
analytical stereomodel will have object coordinates in the ground coordinate system. These points
can then be used for many purposes, such as digital mapping, serving as control for orthophoto
production, or DEM generation. The three orientation steps can be performed as distinct
mathematical operations, or it is possible to combine them in a simultaneous solution.
Photo mosaicking
It is a combination of air photographs, cut and fitted in proper layout and sequence forming
one composite picture of the terrain with continuity of detail. It is prepared for qualitative
study of the terrain before any work is started on the project. For all kinds of preliminary
studies, without difficulty of maps; mosaics are found to be of great help both from speed
and cost of view. A mosaic differs from a map in the following ways: In a map scale is
uniform throughout and ground objects appear in their azimuth or directions while in a
mosaic this is not so because of distortions by tilt and relief; A mosaic is more useful to a
layman who is not conversant with the map symbols and interpretations as the detail in a
mosaic appears in its natural form; A mosaic can be prepared in a much shorter times as
compared to map making which is a long and costly affair; andA mosaic can be prepared
only of flat terrain and of vertical photos. It is extremely difficult to make mosaic of an
oblique photo or a vertical photo of highly mountainous terrain because of large scale
difference.
CHAPTER: SIX
Part one Flight line planning
6.1 Introduction
Why flight planning?
Photography available for a particular area could be outdated for applications such as land
use mapping.
10yrs vs 5yrs Rural vs Urban respectively
Available photography may have been taken in the wrong season and time.
Flights are usually scheduled between 10 a.m. and 2 p.m. for maximum illumination and
minimum shadow. When new photography of an area is required, overall, a great deal of
time, effort, and expense go into the planning and implementation of a photographic
mission. Below, we discuss the geometric aspects of the task of flight planning parameters.
photographs in the direction of flight line is known as forward overlap and is kept at about 60%.
The overlap between the adjacent flight lines is known as the side overlap and is
kept at about 20% - 30%.
EXAMPLE
A study area is 10 km long in the east-west direction and 16 km in the north-south direction. A
camera having a 152-mm-focal-length lens and a 230-mm format is to be used. The desired photo
scale is 1: 25,000 and the nominal end-lap and side-lap are to be 60 and 30 percent respectively.
Beginning and ending flight lines are to be positioned along the boundaries of the study area and
the average Elevation is 300 m above datum. Perform the computations necessary to develop a
flight plan.
Solution
1. Determine the flight direction of the aircraft.
Use north-south flight lines. Note that using north-south flight lines minimizes the number of
lines required and consequently the number of aircraft turns and realignments necessary.
2. Find the flying height above terrain:
3. Determine ground coverage per image from film format size and photo scale:
4. Determine ground separation between photos on a line for 40 percent advance per photo (i.e.,
60 percent end-lap):
= uncover-end lap* ground coverage per image from film format size
5. Find the time between exposures if the aircraft speed is 160 km/hr, the time between
exposures is
7) Compute the number of photos per 16-km line by dividing this length by the photo advance.
: Add one photo to each end and round the number up to ensure coverage
Remember that (1+1) is a total of one photo added to each end of the flight line to ensure
complete coverage of the area of interest.
8. Find the distance between flight lines? (Side lap=30%)
9) Find the number of flight lines required to cover the 10-km study area width by dividing this
width by distance between flight lines
Figure of plus