Distortion and Displacement

Distortion and Displacement
Distortion – “…any shift in the position of an image on a photograph that alters the perspective
characteristics of the image.”
Displacement – “…any shift in the position of an image on a photograph that does not alter the
perspective characteristics of the photograph”
Types of Distortion
Types of Displacement
Film & print shrinkage
Curvature of the Earth
Atmospheric refraction of light Image motion (incl. object ht.)
Lens distortions
Tilt
Topography or relief
Lens distortion radiates from the PP
• Causes an image to appear closer to or farther from the PP
• most serious near the edges of the photo
Usually
negligible
Distortion and Displacement
I.
Principal Points (PP) & Conjugate Principal Points (CPP) – Not the only points of interest on
an aerial photo
A)
Nadir & Isocenter – also important, especially on oblique photos.
B)
PP – physical or optical center of an uncropped photo (intersection of diagonal fiducials)
1) If a photo is truly vertical, PP is directly beneath the aircraft
2) Not so on and obliques
3) Distortions due to lens imperfections radiate from the PP, but are typically very
slight….so they are ignored.
C)
Nadir is the gravitational center of the photo, and is directly beneath the plane
regardless of flight attitude (pitch & roll) at the time of film exposure.
1) If tilt is not too great, nadir will be on the photo
2) Displacement due to elevation differences are radial from Nadir
3) Displacement do to elevation differences can cause serious problems when trying to
measure distance and directions accurately
4) Displacement do to elevation differences is what allows us to see photos in stereo
and will provide us with 2 of the 3 techniques for heights from aerial photos.
Distortion and Displacement
Extreme effects of topographic displacement and distortion…
•
•
Can we accurately measure scale, distance, or area using this photo?
Is average PRS helpful for area calculation?
Displacement
However, Terrain Variation…
• displaces both the top & bottom of
vertical objects.
• How much & direction (toward/away
from Nadir) depends on magnitude of
elev. ∆ (+/-) w/respect to Nadir elev.
On a vertical photograph over flat terrain…
• The farther any vertical object is from
nadir…the greater the object’s top is
displaced away from nadir
• The object’s bottom position remains
stationary (i.e., no displacement)
Distortion and Displacement
D) Isocenter – is always on a line
through the PP & Nadir
called the Principal Meridian,
and is located half way
between these two points.
1) Displacement due to tilt is radial from Isocenter
2) Displacement due to tilt is zero on a truly vertical photos, but increases proportionally as
tilt angle increases.
3) Tilt displacement on slightly tilted photos is usually less in magnitude than displacement
from elevation differences, but tilt is much more difficult to detect, calculate, and correct.
(easier to detect in Chicago vs. rural Iowa…why?)
•
Images on “up side” of a tilted photo are
displaced toward photo center
•
Images on “down side” of titled photo
displace away from photo center
Topographic Displacement
For photo scale problems, we have used these symbols:
D = Ground distance  GD
d = Photo distance  PD
Hg = Flying height above ground  H
From now on used these symbols
GD = Ground distance
PD = Photo distance
H = Flying height above ground/datum
Why? Due to symbology overlap with displacement problems.
H
h
d
r
A
E
=
=
=
=
=
=
Flying height above ground
Height (+/–) of an object on the ground (tree, building, etc.)
Displacement distance on photo (inches or mm)
Radial distance on photo (inches or mm) from nadir to the displaced point
Altitude of aircraft above MSL
Elevation of ground level or datum
Displacement:
All about Similar Triangles
In this Example…
r = photo dist.: nadir to building top
r’ = photo dist.: nadir to building base
d = displacement dist.: top – base  r – r’
Displacement ( d )
d = r – r’
Consider following relationships:
f  r
as H-h  R
f  r’ as H  R
Since d = r – r’
We have to solve for
r and r’ separately…
MSL
Topographic Displacement:
1
𝑓
𝐻−ℎ
=
𝑟
𝑅
∴
𝑓(𝑅)
𝑟=
𝐻−ℎ
All about Similar Triangles
𝑑 = 𝑟 − 𝑟′
𝑓𝑅
𝑑 = 𝐻−ℎ −
𝑓𝑅
𝐻
𝐻(𝑓𝑅)
𝑑 = 𝐻(𝐻−ℎ) −
2
𝑓
𝐻
=
𝑟′
𝑅
∴
𝑓(𝑅)
𝑟′ =
𝐻
𝐻−ℎ 𝑓𝑅
𝐻 𝐻−ℎ
𝑑=
𝐻(𝑓𝑅)
𝐻(𝐻−ℎ)
𝑑=
𝐻𝑓𝑅 − 𝐻𝑓𝑅 + ℎ𝑓𝑅
𝐻(𝐻−ℎ)
ℎ𝑓𝑅
𝑑 = 𝐻(𝐻−ℎ)
𝑑=
𝑟ℎ
𝐻
−
𝐻𝑓𝑅−ℎ𝑓𝑅
𝐻 ℎ−𝐻
r
𝑑 = 𝑟 − 𝑟′
Topographic Displacement at ISU
Measurements
based on…
Photo # 30
Scale 1:24100
H = 12000’
NADIR
Suburb
E = 1035.07’
hsub = +119.09’
River
𝑫𝒊𝒔𝒑𝒍𝒂𝒄𝒆𝒎𝒆𝒏𝒕: 𝒅 =
Actual
Actual
𝑑𝑟𝑖𝑣
E = 876.05’
hriv = -39.93’
𝒓𝒉
𝑯
rriv  nadir to river  13.15 * 0.2 = 2.63”
rsub  nadir to suburb  18.65 * 0.2 = 3.73”
𝑟𝑟𝑖𝑣 ∗ ℎ𝑟𝑖𝑣
2.63" ∗ (−39.93′)
=
=
= −0.00875"
𝐻𝑛𝑎𝑑𝑖𝑟
12000′
𝑑𝑠𝑢𝑏
Displacement
NADIR
E = 915.98’
Correct map
position 0.00875”
farther away
from nadir
𝑟𝑠𝑢𝑏 ∗ ℎ𝑠𝑢𝑏
3.73" ∗ (+119.09′)
=
=
= +0.03701"
𝐻𝑛𝑎𝑑𝑖𝑟
12000′
Terrain displacement is
extremely slight on 1:24100
photos of ISU Campus
Correct map
position 0.037”
closer to nadir
Topographic Displacement:
All about Similar Triangles
𝒅=
𝒓𝒉
𝑯
We get: 𝒉 =
𝒅𝑯
𝒓
From:
Measurements
based on…
Photo # 24
𝒓
𝒓′
𝒅
H = 3000’
NADIR
𝑟 ′ = 22.45 * 0.2 = 4.49”
𝑟 = 22.95 * 0.2 = 4.59”
Carver Hall
𝑑 = 𝑟 − 𝑟′
𝑑 = 4.59" − 4.49 = 0.1"
NADIR
ℎ=
Ground
Note: inches cancel
out leaving h in feet.
0.1" ∗ 3000′
300′
=
= 65.36′ ≈ 65′
4.59"
4.59
Tilt Displacement: Radial from Isocenter
• Is this a truly “vertical” photograph?
• How can one tell?
• Where are the PP & Nadir points?
Displacement due to Tilt: 𝒅𝒕
If you take a photo using a 6” CFL with 3o of Tilt (𝒕), how much displacement is there in an
image if Y is 4”?
Principal Meridian
“Nose-up” side of tilt (horizon side)
in radians
(4")2
𝑑𝑡 =
6"
− 4"
0.0532
dt
“Isoline” or
tilt axis
PP
I
Y
𝑑𝑡 = 0.14" ⇒ Correct map position
actually 0.14” farther
from nadir on photo
N
displacement
Tilt Displacement
On a tilted photograph, the nadir point was determined by intersecting lines passing
through perfectly tall and clearly visible features on the photograph.
The distance between the nadir & the principal point was measured to be 0.5 inches.
What was the angle of tilt of the camera at the time of exposure if a 6 inch CFL lens was
used?
|𝑃𝑃 − 𝑛𝑎𝑑𝑖𝑟|
𝑇𝑖𝑙𝑡 𝐴𝑛𝑔𝑙𝑒 𝜃 = 𝑡𝑎𝑛−1
𝐶𝐹𝐿
𝑇𝑖𝑙𝑡 𝐴𝑛𝑔𝑙𝑒 𝜃 = 𝑡𝑎𝑛−1
0.5 𝑖𝑛.
6 𝑖𝑛.
= 4.764𝑜
Nothing in radians!
Another Topographic Displacement Problem
What can be done to make displacement as small as possible?
• The smaller the flying height or CFL, the greater the displacement.
• The taller the object (or deeper the hole) the more serious the displacement.
• Displacement is more serious the further we move from the nadir.
Suppose we fly at 13,750ft. and photograph a level terrain with a single 200ft. hill
(above datum elev.)…
How serious is photo displacement if 𝒓 = 𝟐. 𝟐𝟓“ and scale is 1:20,000?
𝒅=
𝒓𝒉
𝑯
=
2.25" ∗ 200𝑓𝑡.
13,750𝑓𝑡.
= 0.033“ …on the photo. What’s the ground distance?
1′
𝐺𝐷 = 0.033" ∗ 20,000 ∗
= 55′
12"
Another Height Problem
Rearrange displacement formula to get height formula as before…
𝑑=
𝑟ℎ
𝐻
⇒
ℎ=
𝑑𝐻
𝑟
NOTE: 𝒉 formula only works when 𝒅 (i.e., 𝒓 & 𝒓′) is actually visible on the photo
Suppose we are using photos at a scale of 1:10,000.
The flying height for the photos was 10,000 feet. The base of a building can be
seen 3.16 inches (𝑟′)from the nadir. The top of the same building is 3.18 inches (𝑟)
from the nadir.
What is the height of the building?
𝒉=
𝒅𝑯
3.18"−3.16" ∗ 10000𝑓𝑡.
=
= 63𝑓𝑡.
𝒓
3.8"
Area Determination Error Due to TILT
A
B
Above Datum…
Area is over estimated
Below Datum…
Area under estimated
Solid lines around A & B = true map positions of 2000-feet square areas at MSR = 12,000 on a single
photo…
Dashed lines show photo position & shape of same areas on PSR = 12,000 photo
•
Area A is 600’ higher than Nadir elevation (1000’)
• Here, PSR can be re-calculated & correct areas can be obtained
•
Left extreme of area B is 600’ above nadir elev. & right extreme is 600’ below Nadir elev.
• No mathematical correction is possible for average scale difference.
• In B, can only estimate area accurately near the axis of tilt
• If one of the photos in a stereo pair captures “your area” near Nadir, use that photo.
Effects of Topo Displacement on Area Determination
A more subtle example….
Week 4
Wednesday
Height Determination: The Shadow Method
2 methods for measuring heights on a single aerial photo:
1) Topographic displacement method
•
•
•
•
Object being measured must be vertical (top  bottom)
Object must be far enough from nadir to see displacement
Photo scale must be large enough to reliably measure displacement
Must be able to see both top & bottom of the object
2) Sun-Angle Shadow method
•
•
•
•
•
Must be able to measure full length of object’s shadow
Object’s top must be sharply pointed/distinct and transfers to shadow
Must know sun elevation angle above horizon
Must know the precise time the photo was taken
The rest is simple trigonometry to get object height (ℎ)
Actually,
somewhat
rare!
Height Determination: The Shadow Method
Not a trivial task to determine sun elevation angle!
1) Need the specific Lat/Lon of the object you wish to measure…
•
•
Can be interpolated from 1:24000, 7.5 minute Topoquad maps
Can also use Google Earth or other GIS sources
2) Need the exact time of photo exposure
•
•
Photos we use in class have this info clipped off
The original still has it…
3) Need to plug time and location (Lat/Lon) into a Sun Elevation Calculator
or use solar altitude tables is info above is known.
•
http://www.esrl.noaa.gov/gmd/grad/solcalc/azel.html
Height Determination: The Shadow Method
Perfect!
1) Vertical
2) Pointy top
Shadow correct
Shadow long
Shadow long
Instances that lead to error
in height estimation using
tree shadows
Shadow short
Shadow short
Shadow long
Shadow top
not resolved
on photo
Brush or snow
Shadow short
Shadow short
Height Determination: The Shadow Method
So, if all object conditions are met and we have all the photo ephemeris info…
Then height by shadow measurement is a simple process:
Sun Elevation Angle (𝜃)
Object Height
Shadow Length
𝐹𝑖𝑟𝑠𝑡, 𝑐𝑜𝑛𝑣𝑒𝑟𝑡 𝜽 𝑓𝑟𝑜𝑚 𝑑𝑒𝑔𝑟𝑒𝑒𝑠 𝑡𝑜 𝑟𝑎𝑑𝑖𝑎𝑛𝑠
𝑂𝑏𝑗𝑒𝑐𝑡 𝐻𝑒𝑖𝑔ℎ𝑡
𝑇𝑎𝑛 𝜃 =
𝑆ℎ𝑎𝑑𝑜𝑤 𝐿𝑒𝑛𝑔𝑡ℎ
⇒
𝑑𝑒𝑔𝑟𝑒𝑒𝑠 ∗ 𝜋
180
𝐻𝑒𝑖𝑔ℎ𝑡 = 𝑇𝑎𝑛 𝜃 ∗ (𝑆ℎ𝑎𝑑𝑜𝑤 𝐿𝑒𝑛𝑔𝑡ℎ)
Shadow Method: Determining Object Height
1) Begin by measuring the photo length (PD) of the shadow and converting that
to a ground distance (GD) using the formula….
1
𝑃𝐷"
=
𝑃𝑆𝑅
𝐺𝐷"
⇒
𝑃𝑆𝑅 ∗ 𝑃𝐷" = GD“
1′
...GD" ∗
= 𝐺𝐷′
12"
2) Calculate the tangent of sun elevation angle & multiply by shadow length
Example:
1) You determine the shadow length of a radio tower is 200 ft.
2) You find that the sun elevation angle was 43o 
3) you find… 𝑇𝑎𝑛 0.750492 = 0.932515
4) With this, Height = 0.932515 ∗ 200 ft. = 𝟏𝟖𝟔. 𝟓 𝐟𝐭.
43 ∗ 𝜋
180
= 0.750492 radians
Shadow Method: Determining Sun Elevation
A more difficult problem is determining the sun’s elevation angle when an
object’s height and its shadow length are known.
For this, you can find sun's elevation angle by solving the height formula for
tangent of sun's elevation angle.
𝐻𝑒𝑖𝑔ℎ𝑡 = 𝑻𝒂𝒏 𝜽 ∗ (𝑆ℎ𝑎𝑑𝑜𝑤 𝐿𝑒𝑛𝑔𝑡ℎ)
For Example:
1) Suppose you know an object is 50 ft. tall and the shadow length = 38 ft.
2) 𝑇𝑎𝑛 𝑆𝑢𝑛′ 𝑠 𝑒𝑙𝑒𝑣𝑎𝑡𝑖𝑜𝑛 𝑎𝑛𝑔𝑙𝑒 =
3) This
=
50 𝑓𝑡.
38 𝑓𝑡.
𝐻𝑒𝑖𝑔ℎ𝑡
𝑆ℎ𝑎𝑑𝑜𝑤 𝐿𝑒𝑛𝑔𝑡ℎ
= 1.3158
4) Arctan(1.3158) = 0.92093 radians

0.92093 ∗180
𝜋
= 52.8𝑜
Determining Photo-Specific Solar Position
Most often, you won't know the height of objects casting shadows, so…
Options:
1) Can use a solar altitude table which is based on latitude (LAT) of the
photography & time the photo was taken.
•
Will give an approximate value for sun's altitude depending on how
closely LAT & time match those available in the table
•
It is possible to interpolate the table for time, date, & LAT to get an exact
answer for sun's altitude
2) Or, can use solar ephemeris to determine Sun’s altitude. To use ephemeris
must know the following:
•
•
•
•
Angle X ….Sun’s declination (LAT) on the day of the photo (look up in
ephemeris
Angle Y ….LAT of photography (find this on a map)
Angle Z ….is the hour angle, or difference in longitude (LON) between
position of the sun and locality of the photography (can be calculated)
𝑆𝑖𝑛 𝐴𝑙𝑡. 𝐴𝑛𝑔𝑙𝑒 = 𝐶𝑜𝑠𝑋 𝐶𝑜𝑠𝑌 𝐶𝑜𝑠𝑍 + 𝒐𝒓 − 𝑆𝑖𝑛 𝑋 𝑆𝑖𝑛𝑌
+ from March 21 – Sept. 23 in Northern Hemisphere
– from Sept. 24 – March 20 in Southern Hemisphere
FYI only
Determining Photo-Specific Solar Position (Option 1)
Sun Position Calculator: http://aa.usno.navy.mil/data/docs/AltAz.php
output
Determining Photo-Specific Solar Position (Option 1)
Height Determination: The Shadow Method
Very hard to read clock!
Requires use of the dreaded
Magnifying Comparator
The Shadow Method: In Presence of Slopes
When shadows fall on either an uphill or downhill slope, the method
𝑯𝒆𝒊𝒈𝒉𝒕 = 𝑻𝒂𝒏 𝜽 ∗ (𝑺𝒉𝒂𝒅𝒐𝒘 𝑳𝒆𝒏𝒈𝒕𝒉)
…will give incorrect height estimates unless adjustments are made.
An adjustment formula can be used if degree of slope is known …(via clinometer etc.)
𝑯𝒆𝒊𝒈𝒉𝒕 = 𝑺𝒉𝒂𝒅𝒐𝒘𝑳𝒆𝒏 ∗ 𝑪𝒐𝒔 𝑺𝒍𝒐𝒑𝒆 ∗ [𝑻𝒂𝒏 𝑺𝒖𝒏𝑬𝒍𝒆𝒗 + 𝒐𝒓 − 𝑻𝒂𝒏 𝑺𝒍𝒐𝒑𝒆 ]
+ if shadow falls uphill
− if shadow falls downhill
Angles in Radians 
𝜃∗𝜋
180
Example:
1:15840 photo is taken at 11:30 CST on August 20st at 40o N LAT . A tree casts a shadow on a
uphill face of a 35o slope. On the photo the shadow measures 0.03”, what is the tree’s height?
ℎ = 0.03" ∗ 𝐶𝑜𝑠 0.610851 ∗ 𝑇𝑎𝑛 1.064651 + 𝑇𝑎𝑛 0.610851
ℎ = 0.03" ∗ 0.81916 ∗ (1.804048 + 0.700186)
ℎ = 0.03" ∗ 2.051369
ℎ = 0.061541" ∗ (15840/12) = 𝟖𝟏. 𝟐′
≈ 𝟖𝟐′
Get sun elev. from SOLAR TABLE
Hotspot
If the sun's altitude is greater than 90 minus 1/2 the lens angle of the camera being used
to take the photos, rays from the sun can be reflected directly onto the film causing a sun
spot to appear on the photo
> 60o
60o to 52.5o
52.5o to 40o
< 40o
Hotspot
>
Sun's Image or Hot Spot
will cause loss of detail in a
small area. It appears as a
white circular spot on the
photo.
=
<
No Hotspot
𝜃 = 90 −
𝜃
1
∗ Lens Angle
2
Hotspot & No Shadow Spot
No Shadow Spot: same distance from the Principal
point as the Sun's Image along a line drawn from the
sun's image through the PP. No Shadow point will be
on the opposite side of PP from the sun's image
PP
No
Shadow
spot
Hotspot
No Shadow Spot
1)
No Shadow Spot: opposite side of PP from the sun's image (hotspot)
on a line running through the PP.
2)
Large area where no shadows will appear on the photo
•
Caused by the fact that the displaced images of objects fall
directly on top of the shadow being cast by the object.
•
Result is loss of detail over a large area & reduction in ability to
interpret aerial photos.
•
Ideally, would like to obtain photos in such a way that both sun
spot & no shadow spot are avoided.
VS.