Image Enhancement Spatial Domain PDF
Image Enhancement Spatial Domain PDF
Image Enhancement Spatial Domain PDF
Spatial Domain
Principle Objective of Enhancement
3
Good images
◼ For human visual
◼ The visual evaluation of image quality is a highly subjective
process.
◼ It is hard to standardize the definition of a good image.
◼ For machine perception
◼ The evaluation task is easier.
◼ A good image is one which gives the best machine recognition
results.
◼ A certain amount of trial and error usually is required
before a particular image enhancement approach is
selected.
4
Spatial Domain
◼ Procedures that operate
directly on pixels.
g(x,y) = T[f(x,y)]
where
◼ f(x,y) is the input image
◼ T is an operator on f defined
over some neighborhood of
(x,y)
5
Mask/Filter
◼ Neighborhood of a point
(x,y) can be defined by using a
(x,y) square/rectangular (common
used) or circular subimage
• area centered at (x,y)
◼ The center of the subimage is
moved from pixel to pixel
starting at the top of the
corner
6
Point Processing
◼ Neighborhood = 1x1 pixel
◼ g depends on only the value of f at (x,y)
◼ T = gray level (or intensity or mapping) transformation
function
s = T(r)
◼ Where
◼ r = gray level of f(x,y)
◼ s = gray level of g(x,y)
7
Neighbourhood Processing
Mask Processing or Filter
◼ Neighborhood is bigger than 1x1 pixel
◼ Use a function of the values of f in a predefined
neighborhood of (x,y) to determine the value of g at
(x,y)
◼ The value of the mask coefficients determine the
nature of the process
◼ Used in techniques
◼ Image Sharpening
◼ Image Smoothing
8
Transformation
Identity: s=r, no transformation
L=5
9
Point Processing Techniques :
1. Digital Negative.
2. Contrast Stretching
3. Thresholding
4. Gray level Slicing
5. Bit Plane Slicing
6. Dynamic range compression
10
3 basic gray-level transformation
functions
Negative
◼ Linear function
nth root ◼ Negative and identity
Log
transformations
nth power
◼ Logarithm function
◼ Log and inverse-log
transformation
Identity Inverse Log
◼ Power-law function
◼ nth power and nth root
Input gray level, r transformations
11
Image Negatives
◼ An image with gray level in the range [0, L-1]
where L = 2n ; n = 1, 2…
◼ Negative transformation :
s = L – 1 –r
12
Image Negatives
▪ Image in the range [0, L-1].
L =5
S= 4- r
13
Example of Negative Image
14
Contrast Stretching
◼ Produce higher contrast
than the original by
◼ darkening the levels below
m in the original image
◼ Brightening the levels
above m in the original
image
15
Contrast Stretching
◼ increase the dynamic range of
the gray levels in the image
◼ (b) a low-contrast image : result
from poor illumination, lack of
dynamic range in the imaging
sensor, or even wrong setting of
a lens aperture of image
acquisition
◼ (c) result of contrast stretching:
(r1,s1) = (r min,0) and (r2,s2) =
(r max,L-1)
◼ (d) result of thresholding
16
Thresholding
◼ Produce a two-level
(binary) image
17
18
Log Transformations
s = c log (1+r)
Negative
nth root ◼ c is a constant
and r 0
Log
nth power
◼ Log curve maps a narrow
range of low gray-level values
in the input image into a
Identity Inverse Log wider range of output levels.
◼ Used to expand the values of
dark pixels in an image while
Input gray level, r
compressing the higher-level
values.
19
Log Transformations
◼ It compresses the dynamic range of images with large
variations in pixel values
◼ Example of image with dynamic range: Fourier
spectrum image
◼ It can have intensity range from 0 to 106 or higher.
◼ We can’t see the significant degree of detail as it will be
lost in the display.
20
Example of Logarithm Image
22
Power-Law Transformations
s = cr
◼ c and are positive constants
◼ Power-law curves with
fractional values of map a
Output gray level, s
23
Gamma correction
26
Another example a b
c d
(a) image has a washed-out
appearance, it needs a
compression of gray levels
needs > 1
(b) result after power-law
transformation with = 3.0
(suitable)
(c) transformation with = 4.0
(suitable)
(d) transformation with = 5.0
(high contrast, the image has
areas that are too dark, some
detail is lost)
27
Gray-level slicing
◼ Highlighting a specific
range of gray levels in an
image
◼ Display a high value of all
gray levels in the range of
interest and a low value for
all other gray levels
◼ (a) transformation highlights
range [A,B] of gray level and
reduces all others to a constant
level
◼ (b) transformation highlights
range [A,B] but preserves all
other levels
28
29
Applications
30
Bit-plane slicing
◼ Highlighting the contribution
made to total image
appearance by specific bits
◼ Suppose each pixel is
One 8-bit byte Bit-plane 7
(most significant) represented by 8 bits
◼ Higher-order bits contain the
majority of the visually
significant data
◼ Useful for analyzing the
Bit-plane 0 relative importance played by
(least significant)
each bit of the image
31
Example
◼ The (binary) image for bit-
plane 7 can be obtained by
processing the input image
with a thresholding gray-level
transformation.
◼ Map all levels between 0 and
127 to 0
◼ Map all levels between 129 and
An 8-bit fractal image 255 to 255
32
8 bit planes
Bit-plane 7 Bit-plane 6
33
Histogram Processing
◼ Histogram of a digital image with gray levels in
the range [0,L-1] is a discrete function
h(rk) = nk
◼ Where
◼ rk : the kth gray level
◼ nk : the number of pixels in the image having gray
level rk
◼ h(rk) : histogram of a digital image with gray levels
rk
34
Normalized Histogram
◼ dividing each of histogram at gray level rk by the
total number of pixels in the image, n
p(rk) = nk / n
◼ For k = 0,1,…,L-1
◼ p(rk) gives an estimate of the probability of
occurrence of gray level rk
◼ The sum of all components of a normalized
histogram is equal to 1
35
Example
No. of pixels
6
2 3 3 2 5
4 2 4 3 4
3
3 2 3 5
2
2 4 2 4
1
Gray level
4x4 image
0 1 2 3 4 5 6 7 8 9
Gray scale = [0,9]
histogram
36
300
250
200
150
100
50
0
0 50 100 150 200 250 300
37
300
250
200
150
100
50
0
0 10 20 30 40 50 60 70
38
Histogram Processing
◼ Basic for numerous spatial domain processing
techniques
◼ Used effectively for image enhancement
◼ Information inherent in histograms also is useful
in image compression and segmentation
39
h(rk) or p(rk)
Example
rk
Dark image
Components of
histogram are
concentrated on the
low side of the gray
scale.
Bright image
Components of
histogram are
concentrated on the
high side of the gray
scale.
40
Example
Low-contrast image
histogram is narrow
and centered toward
the middle of the
gray scale
High-contrast image
histogram covers broad
range of the gray scale
and the distribution of
pixels is not too far from
uniform, with very few
vertical lines being much
higher than the others
41
Histogram Equalization
◼ As the low-contrast image’s histogram is narrow and
centered toward the middle of the gray scale, if we
distribute the histogram to a wider range the quality of
the image will be improved.
◼ We can do it by adjusting the probability density
function of the original histogram of the image so that
the probability spread equally.
42
Histogram transformation
s = T(r)
s
◼ Where 0 r 1
1
◼ T(r) satisfies
◼ (a). T(r) is single-
valued and
sk= T(rk)
monotonically
T(r) increasingly in the
interval 0 r 1
◼ (b). 0 T(r) 1 for
0r1
0 rk 1 r
43
2 Conditions of T(r)
45
Applied to Image
◼ Let
◼ pr(r) denote the PDF of random variable r
◼ ps (s) denote the PDF of random variable s
47
Transformation function
49
Finding ps(s) from given T(r)
ds dT ( r )
=
dr dr
dr
p s ( s ) = pr ( r )
r
d
= pr ( w )dw ds
dr 0
1
= pr ( r ) = pr ( r )
pr ( r )
51
r
s = T ( r ) = pr ( w )dw
0
yields
Ps(s)
a random variable s 1
characterized by
a uniform probability
function s
0
52
Discrete
transformation function
◼ The probability of occurrence of gray level in an image is approximated by
nk
pr ( rk ) = where k = 0 , 1, ..., L-1
n
54
Example
Histogram
before after equalization
55
Example
Histogram
before after equalization
The quality is
not improved
much because
the original
image already
has a broaden
gray-level scale
56
Example
No. of pixels
6
2 3 3 2 5
4 2 4 3 4
3
3 2 3 5
2
2 4 2 4
1
Gray level
4x4 image
0 1 2 3 4 5 6 7 8 9
Gray scale = [0,9]
histogram
57
Gray
0 1 2 3 4 5 6 7 8 9
Level(j)
No. of
0 0 6 5 4 1 0 0 0 0
pixels
k
n
j =0
j 0 0 6 11 15 16 16 16 16 16
k nj 6 11 15 16 16 16 16 16
s= 0 0 / / / / / / / /
j =0 n
16 16 16 16 16 16 16 16
3.3 6.1 8.4
sx9 0 0 9 9 9 9 9
3 6 8
Example
No. of pixels
6
3 6 6 3 5
8 3 8 6 4
3
6 3 6 9
2
3 8 3 8
1
Output image
0 1 2 3 4 5 6 7 8 9
Gray scale = [0,9] Gray level
Histogram equalization
59
Note
60
Example
Histogram
before after equalization
61
Example
Histogram
before after equalization
The quality is
not improved
much because
the original
image already
has a broaden
gray-level scale
62
Neighborhood Processing
Smoothing Spatial Filters
64
Spatial Filtering Process
◼ simply move the filter mask from point to point
in an image.
◼ at each point (x,y), the response of the filter at
that point is calculated using a predefined
relationship.
R = w1 z1 + w2 z 2 + ... + wmn z mn
mn
= wi zi
i =i
65
Linear Filtering
◼ Linear Filtering of an image f of size MxN filter
mask of size mxn is given by the expression
a b
g ( x, y ) = w(s, t ) f ( x + s, y + t )
t =− a t =−b
67
Smoothing Linear Filters
68
Smoothing Linear Filters
◼ replacing the value of every pixel in an image by the
average of the gray levels in the neighborhood will
reduce the “sharp” transitions in gray levels.
◼ sharp transitions
◼ random noise in the image
◼ edges of objects in the image
◼ thus, smoothing can reduce noises (desirable) and blur
edges (undesirable)
69
3x3 Smoothing Linear Filters
71
a b
Example c d
e f
74
Median Filters
◼ replaces the value of a pixel by the median of
the gray levels in the neighborhood of that pixel
(the original value of the pixel is included in the
computation of the median)
◼ quite popular because for certain types of
random noise (impulse noise salt and pepper
noise) , they provide excellent noise-reduction
capabilities, with considering less blurring than
linear smoothing filters of similar size.
75
order (0, 2, 3, 3, 4, 6, 10, 19, 97) Median : 4
76
Example : Median Filters
77
Sharpening Spatial Filters
◼ to highlight fine detail in an image
◼ or to enhance detail that has been blurred, either
in error or as a natural effect of a particular
method of image acquisition.
78
Blurring vs. Sharpening
79
Derivative operator
◼ the strength of the
response of a derivative
operator is proportional to
the degree of discontinuity
of the image at the point
at which the operator is
applied.
◼ thus, image differentiation
◼ enhances edges and other
discontinuities (noise)
◼ deemphasizes area with
slowly varying gray-level
values.
80
original blurred image High pass filtered
81
Unsharp masking
f s ( x, y ) = f ( x, y ) − f ( x, y )
sharpened image = original image – blurred image
96
High-boost filtering
f hb ( x, y ) = Af ( x, y ) − f ( x, y )
f hb ( x, y ) = ( A − 1) f ( x, y ) − f ( x, y ) f ( x, y )
= ( A − 1) f ( x, y ) − f s ( x, y )
◼ generalized form of Unsharp masking
◼ A1
97
High-boost filtering
f hb ( x, y ) = ( A − 1) f ( x, y ) − f s ( x, y )
◼ if we use Laplacian filter to create sharpen image
fs(x,y) with addition of original image
f ( x, y ) − f ( x, y )
2
f s ( x, y ) =
f ( x, y ) + f ( x, y )
2
98
High-boost filtering
Af ( x, y ) − f ( x, y )
2
f hb ( x, y ) =
Af ( x, y ) + f ( x, y )
2
◼ A1
◼ if A = 1, it becomes “standard” Laplacian
sharpening
100
Example
101
f
Gx x
Gradient Operator f = = f
G y
y
◼ first derivatives are implemented using the
magnitude of the gradient.
1
f = mag (f ) = [G + G ] 2
x
2
y
2
commonly approx.
1
f 2 f 2 2
= +
x y
f G x + G y
the magnitude becomes nonlinear
102
Gradient Mask z1 z2 z3
z4 Z5 z6
◼ simplest approximation, 2x2 z7 Z8 z9
G x = ( z8 − z 5 ) and G y = ( z 6 − z5 )
1 1
f = [G + G ]
2
x
2
y
2
= [( z8 − z5 ) + ( z6 − z5 ) ]
2 2 2
f z8 − z 5 + z 6 − z 5
103
z1 z2 z3
Gradient Mask z4 z5 z6
z7 z8 z9
◼ Roberts cross-gradient operators, 2x2
Gx = ( z9 − z5 ) and G y = ( z8 − z 6 )
1 1
f = [G + G ]
2
x
2
y
2
= [( z9 − z5 ) + ( z8 − z6 ) ]
2 2 2
f z 9 − z 5 + z8 − z 6
104
z1 z2 z3
Gradient Mask z4 z5 z6
z7 z8 z9
106
Example
107
Example of Combining Spatial
Enhancement Methods
◼ want to sharpen the
original image and bring
out more skeletal detail.
◼ problems: narrow
dynamic range of gray
level and high noise
content makes the image
difficult to enhance
108
Example of Combining Spatial
Enhancement Methods
◼ solve :
1. Laplacian to highlight fine detail
2. gradient to enhance prominent edges
3. gray-level transformation to increase
the dynamic range of gray levels
109
110
111