CSE4227 Digital Image Processing
Chapter 03 – Sharpening Spatial Filter
Dr. Kazi A Kalpoma
Professor, Department of CSE
Ahsanullah University of Science & Technology (AUST)
Contact: [email protected]
Google Class code: o5sjihh
CSE | AUST Fall 2022
Today’s Contents
❑Sharpening spatial filter
❑Derivatives of Image
• 1st derivative
• 2nd derivative
❑Laplacian Filter
❑Laplacian Image Enhancement
❑Gradient Operators
❑Difference filters
❑Combining filtering techniques
Chapter 3 from R.C. Gonzalez and R.E. Woods, Digital Image Processing
(3rd Edition), Prentice Hall, 2008 [ Section 3.6, 3.7 ]
Sharpening Spatial Filters
Previously we have looked at smoothing filters which
remove fine detail.
Sharpening spatial filters seek to highlight fine detail.
◼ Remove blurring from images
◼ Highlight edges
◼ Useful for emphasizing transitions in image
intensity
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Sharpening Spatial Filters
DIP, Spring 2012 GS & AS, Bahria University, Islambad 4
Spatial Differentiation
◼ Sharpening filters are based on first- and
second-order derivatives of image.
◼ The derivatives of a digital function are
defined in terms of differences.
◼ i.e. Sharpening filters are based on spatial
differentiation
◼ Differentiation measures the rate of change
of a function.
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Spatial Differentiation
Let’s consider a simple 1 dimensional example
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Spatial Differentiation
A B
DIP, Spring 2012 GS & AS, Bahria University, Islambad
1 st Derivative in Digital Form
❑ The 1st derivative of a function is given by:
f
= f ( x + 1) − f ( x)
x
❑ Its just the difference between subsequent
values and measures the rate of change of the
function
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Image Strip
1st Derivative
8
7
6
5
4
3
2
1
0
5 5 4 3 2 1 0 0 0 6 0 0 0 0 1 3 1 0 0 0 0 7 7 7 7
-1 -1 -1 -1 -1 0 0 6 -6 01st 0Derivative
0 1 2 -2 -1 0 0 0 7 0 0 0
8
6
4
2
0
-2
-4
-6
-8
Derivative is nonzero along the entire ramp, zero in flat area,.
DIP, Spring 2012 GS & AS, Bahria University, Islambad
2 nd Derivative in Digital Form
The 2nd derivative of a function is given by:
f 2
= f ( x + 1) + f ( x − 1) − 2 f ( x)
x
2
Simply takes into account the values both before
(backward) and after (forward) the current value
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Image Strip
8
2nd Derivative
7
6
5
4
3
2
1
0
5 5 4 3 2 1 0 0 0 6 0 0 0 0 1 3 1 0 0 0 0 7 7 7 7
-1 0 0 0 0 1 0 6 -12 6 0 0 1 1 -4 1 1 0 0 7 -7 0 0
2nd Derivative
10
-5
-10
-15
Derivative is nonzero at the onset and end of ramp, stronger
DIP, Spring 2012
response
GS &at
AS,and aroundIslambad
Bahria University, the point.
Sharpening Spatial Filters
1. LAPLACIAN
• Use of 2nd Derivative for Image Enhancement
2. SOBEL (Gradient Operators)
• Use of 1st Derivative for Image Enhancement
Use of 2nd Derivative for Image
Enhancement
The 2nd derivative is more useful for image
enhancement than the 1st derivative - Stronger
response to fine detail
The first sharpening filter we will look at is the Laplacian
DIP, Spring 2012 GS & AS, Bahria University, Islambad
- For Two Dimensional
2 nd Derivative in Two Dimension
DIP, Spring 2012 GS & AS, Bahria University, Islambad 15
1. Laplacian Filter
So, the Laplacian can be given as follows:
2 f = [ f ( x + 1, y) + f ( x − 1, y)
+ f ( x, y + 1) + f ( x, y − 1)]
− 4 f ( x, y)
0 1 0
We can implement it 1 -4 1
using this filter.
0 1 0
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Types of Laplacian Kernels
DIP, Spring 2012 GS & AS, Bahria University, Islambad
0 -1 0
Another Example: -1 4 -1
0 -1 0
+ =
Original Image Laplace Sharpened Laplace filtered
image image
Example: apply the following Laplacian filter on the
highlighted and underlined pixel
153 157 156 153 155
0 -1 0 159 156 158 156 159
-1 4 -1 155 158 154 156 160
154 157 158 160 160
0 -1 0
157 157 157 156 155
Step 1:
154*4 – 158- 156-158-158 = -14
So the value after filter = -14
We call the resultant image: sharpened image.
Step 2:
Filtered image=original + sharpened image
The value in the filtered image=154-14 =130
Laplacian Image Enhancement
0 1 0
1 -4 1
0 1 0
- =
Original Laplacian Sharpened
Image Filtered Image Image
In the final sharpened image edges and fine
detail are much more obvious
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Simplified Image Enhancement
• The result of a Laplacian filtering is not an enhanced image.
• The entire enhancement can be combined into a single
filtering operation
g ( x, y) = f ( x, y) − f 2
= f ( x, y) − [ f ( x + 1, y) + f ( x −1, y)
0 1 0 + f ( x, y + 1) + f ( x, y −1)
1 -4 1
w w w
− 4 f ( x, y )]
1
0 1 0 2 3
f ( x, y ) − 2 f , w5 0
w4 w5 w6
g ( x, y ) =
DIP, Spring 2012
w7 w8 w9 f ( x, y ) + 2 f , w5 0
GS & AS, Bahria University, Islambad
Simplified Image Enhancement
The entire enhancement or sharpening can
be done in one PASS.
g ( x, y) = f ( x, y) − f 2
= 5 f ( x, y) − f ( x + 1, y) − f ( x −1, y)
0 -1 0
− f ( x, y + 1) − f ( x, y −1)
-1 5 -1 We can implement it
0 -1 0 using this filter.
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Simplified Image Enhancement
This gives us a new filter which does the
whole job for us in one step
0 -1 0
-1 5 -1
0 -1 0
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Variants On The Simple
Laplacian
There are lots of slightly different versions of the
Laplacian that can be used:
0 1 0 1 1 1
Simple Variant of
1 -4 1 Laplacian 1 -8 1 Laplacian
0 1 0 1 1 1
-1 -1 -1
-1 9 -1
-1 -1 -1
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Use of 1st Derivatives for Image
Enhancement
The another Sharpening Spatial filters is
SOBEL (Gradient Operators).
What is Gradient of a Digital Image?
DIP, Spring 2012 GS & AS, Bahria University, Islambad
The Gradient of a Digital Image
DIP, Spring 2012 GS & AS, Bahria University, Islambad
The Gradient of a Digital Image
DIP, Spring 2012 GS & AS, Bahria University, Islambad
The Gradient – First-order Derivative
There is some debate as to how best to calculate these gradients.
DIP, Spring 2012 GS & AS, Bahria University, Islambad 28
Gradient Operators
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Gradient Operators
DIP, Spring 2012 GS & AS, Bahria University, Islambad 30
Sharpening Spatial filters :
2) SOBEL (Gradient Operator)
▪ The Sobel operator provides differencing and smoothing
effect of an image.
▪ Sobel operator consists of 3x3 convolution kernels. Gx is
a simple kernel and Gy is rotated by 90°
-1 -2 -1 -1 0 1
0 0 0 -2 0 2
1 2 1 -1 0 1
Gx, Extract horizontal edges Gy, Extract vertical edges
DIP, Spring 2012 GS & AS, Bahria University, Islambad 31
Gradient Operators
Sobel Operator
-1 -2 -1 -1 0 1
f f
= 0 0 0 = -2 0 2
x 1 2 1
y -1 0 1
Gx, Extract horizontal edges Gy, Extract vertical edges
f ( z7 + 2 z8 + z9 ) − ( z1 + 2 z2 + z3 )
+ ( z3 + 2 z6 + z9 ) − ( z1 + 2 z4 + z7 )
Pixel Arrangement
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Gradient Operators
Prewitt Operator
❑ is used for detecting edges horizontally and vertically.
f ( z7 + z8 + z9 ) − ( z1 + z2 + z3 )
+ ( z3 + z6 + z9 ) − ( z1 + z4 + z7 )
Pixel Arrangement
-1 -1 -1 -1 0 1
f f
= 0 0 0 = -1 0 1
x y
1 1 1 -1 0 1
Extract horizontal edges Extract vertical edges
DIP, Spring 2012 Gx GS & AS, Bahria University, Islambad Gy
Sobel Operator: Example
An image of a
contact lens
which is
enhanced in
order to make
defects more
obvious
Sobel filters are typically used for edge
detection
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Sharpening with Sobel
Operator
35
DIP, Spring 2012 GS & AS, Bahria University, Islambad 2016/1/27
Source: D. Lowe
Sharpening with Sobel
Operator
1 0 -1
2 0 -2
1 0 -1
Sobel
Vertical Edge
DIP, Spring 2012 GS & AS, Bahria University, Islambad
(absolute value)
2016/1/27
36
Sharpening with Sobel
Operator
1 2 1
0 0 0
-1 -2 -1
Sobel
Horizontal Edge
DIP, Spring 2012 GS & AS, Bahria University, Islambad
(absolute value)
2016/1/27
37
Sharpening Spatial Filters
Laplacian Sobel
DIP, Spring 2012 GS & AS, Bahria University, Islambad 38
Difference Filter
❑Also called as Emboss filters
❑Enhances the details in the direction specific to the
mask selected
❑Four primary difference filter convolution masks,
corresponding to the edges in the vertical, horizontal, and
two diagonal directions are:
Vertical Horizontal Diagonal 1 Diagonal 2
DIP, Spring 2012 GS & AS, Bahria University, Islambad 39
Difference Filter
Original image Difference filtered image Difference filtered image
added to the original image,
with contrast enhanced
DIP, Spring 2012 GS & AS, Bahria University, Islambad 40
Combining Spatial Enhancement
Methods
Successful image enhancement is
typically not achieved using a single
operation
Rather we combine a range of
techniques in order to achieve a final
result
This example will focus on enhancing
the bone scan
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Combining Spatial Enhancement Methods
(a)
Laplacian filter of
bone scan (a)
(b)
Sharpened version of
bone scan achieved (c)
by subtracting (a)
and (b) Sobel filter of bone
scan (a) (d)
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Combining Spatial Enhancement Methods
Result of applying a (h)
power-law trans. to
Sharpened image (g)
which is sum of (a)
and (f) (g)
The product of (c)
and (e) which will be (f)
used as a mask
(e)
Image (d) smoothed with
a 5*5 averaging filter
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Combining Spatial Enhancement Methods
Compare the original and final images
DIP, Spring 2012 GS & AS, Bahria University, Islambad
Class Work
Consider a 3-bit 4x4 image.
0 2 6 7 Laplacian
1 1 6 4 filter
0 -1 0
4 5 2 7 -1 +4 -1
1 2 6 0 0 -1 0
Find the filtered output image using
• this Laplacian filter,
• a 3 × 3 Mean filter
• a 3x3 Median filter and
• a Sobel operator
Ignore the border pixels in calculation and put zero in the border of
the output image.
