Lec 15 Multiscale Re
Lec 15 Multiscale Re
ANALYSIS
Image Representation
• Wavelet Pyramids
•Gaussian pyramids
• Laplacian Pyramids
• Applications
What is a good representation
for image analysis?
• Pixel domain representation tells you
“where” (pixel location), but not
“what”.
– In space, this representation is too localized
• Fourier transform domain tells you
“what” (textural properties), but not
“where”.
– In space, this representation is too spread out.
• Want an image representation that gives
you a local description of image events
—what is happening where.
– That representation might be “just right”.
Image representation
Fourier Analysis
Breaks down a signal into constituent
sinusoids of different frequencies
Spatial
Freq.
Fourier basis
Spatial
Freq.
Wavelet basis
Spatial
Short Time Fourier Analysis
In order to analyze small section of a
signal, Denis Gabor (1946), developed a
technique, based on the FT and using
windowing : STFT
STFT (or: Gabor Transform)
Spatial
Freq.
Fourier basis
Spatial
Freq.
Wavelet basis
Spatial
Wavelet's properties
.Flexibility •
Scaling
Wavelet analysis produces a time-
scale view of the signal.
Scaling means stretching or
compressing of the signal.
scale factor (a) for sine waves:
f ( t ) sin( t ) ; a 1
f ( t ) sin( 2t ) ; a 1 2
f ( t ) sin( 4t ) ; a 1 4
Scaling (Cont’d)
f ( t ) (t ) ; a 1
f ( t ) ( 2 t ) ; a 1 2
f ( t ) ( 4 t ) ; a 1 4
Wavelet function
b – shift •
a , b x 1
a
x b
a
coefficient
a – scale •
coefficient
a , bx , by x , y 1a
x bx
a , y by
a 2D function •
a , bx , by x , y
a , bx , by x , y
1
1a x bx
x abx ,
,
y by
y aby
a a a
Various Wavelet basis
Wavelet - Frequency domain
Wavelet bands are split recursively
image
H L
H L
H L
Discrete Wavelet Transform
LENA
LH
HL HH
Wavelet - Frequency domain
Wavelet decomposition - 2D
Frequency domain
Edges (derivatives):
f(x)
f (x)
Space Required for Pyramids
Image Pyramids
No details in image -
Low (blurred image)
Resolution low frequencies
Low resolution
High resolution
Image Pyramid
Frequency Domain
Low resolution
High resolution
Image Blurring = low pass filtering
*
=
~
=
*
=
*
=
Image Pyramid
Low resolution
High resolution
Gaussian pyramid construction
filter mask
Repeat
• Filter
• Subsample
Until minimum resolution reached
• can specify desired number of levels (e.g., 3-level pyramid)
The whole pyramid is only 4/3 the size of the original image!
Gaussian Pyramid
Level n
1X1
Level 1
2n-1 X 2n-1
Level 0
2n X 2n
Image Pyramids
High resolution
Gaussian Pyramid
Image pyramids
• Gaussian Pyramids
• Laplacian Pyramids
What does blurring take
?away
smoothed – original
Band-pass filtering
Need this!
Original
image
Gl+1
Expand
Reduce
Gl G’l
Ll = Gl - G’l
original
?What does blurring take away
smoothed – original
Laplacian Pyramid
Gaussian Laplacian
Pyramid Pyramid
expa
nd
ex
p
- =
an
d
ex
p
- =
an
d
- =
The Laplacian Pyramid
Li Gi expand( Gi 1 )
Gaussian Pyramid Gi Li expand( Gi 1 ) Laplacian Pyramid
Gn Ln Gn
ex pa
nd
G2
an-
L2
exp
d
=
G1 L1
ex
pa -
nd =
G0 L0
- =
Image Pyramids - Comparison
Gaussian Pyramid
Laplacian Pyramid
Wavelet Pyramid
Gaussian
Frequency
Laplacian
Pyramid Domain
Pyramid
Laplace Pyramid -
No scaling
50
100
50
150
100
200
150
50
250
50 100 150 200 250
200
100
250
150
50 100 150 200 250
200
250
50 100 150 200 250
50
100
50
150
100
200
150
250
50 100 150 200 250
200
250
50 100 150 200 250
from: B.Freeman
Reconstruction of the original image
from the Laplacian Pyramid
Laplacian Gl = Ll + G’l
Pyramid
expand
+ =
expand
+ =
expand
+ =
Original
= Image
Multiresolution Spline - Example
Glued Splined
laplacian level 4
laplacian level 2
laplacian level 0
High Frequencies
Middle Frequencies
Low Frequencies
Multiresolution Spline
High Frequencies
Middle Frequencies
Low Frequencies
Multiresolution Spline - Example
Wide Transition
© prof. dmartin
Pyramids in Frequency
Domain
Gaussian Pyramid Laplacian Pyramid
Wavelet Decomposition
Fourier Space
Image Fusion