Chrominance Assisted Sharpening of Images

Similar documents
1.Discuss the frequency domain techniques of image enhancement in detail.

Frequency Domain Enhancement

Practical Image and Video Processing Using MATLAB

Image Filtering Josef Pelikán & Alexander Wilkie CGG MFF UK Praha

CoE4TN4 Image Processing. Chapter 3: Intensity Transformation and Spatial Filtering

Image Filtering. Median Filtering

CSE 564: Visualization. Image Operations. Motivation. Provide the user (scientist, t doctor, ) with some means to: Global operations:

Digital Image Processing

IMAGE PROCESSING: AREA OPERATIONS (FILTERING)

Image Enhancement. DD2423 Image Analysis and Computer Vision. Computational Vision and Active Perception School of Computer Science and Communication

Chapter 3. Study and Analysis of Different Noise Reduction Filters

CSE 564: Scientific Visualization

EE482: Digital Signal Processing Applications

Digital Image Processing

A Spatial Mean and Median Filter For Noise Removal in Digital Images

Using the Advanced Sharpen Transformation

Prof. Feng Liu. Winter /10/2019

Filip Malmberg 1TD396 fall 2018 Today s lecture

Table of contents. Vision industrielle 2002/2003. Local and semi-local smoothing. Linear noise filtering: example. Convolution: introduction

Image analysis. CS/CME/BIOPHYS/BMI 279 Fall 2015 Ron Dror

Image Enhancement in spatial domain. Digital Image Processing GW Chapter 3 from Section (pag 110) Part 2: Filtering in spatial domain

Image Enhancement. Image Enhancement

CoE4TN4 Image Processing. Chapter 4 Filtering in the Frequency Domain

Lecture 3: Linear Filters

Image Processing. Adam Finkelstein Princeton University COS 426, Spring 2019

IMAGE ENHANCEMENT IN SPATIAL DOMAIN

Spatial Domain Processing and Image Enhancement

DIGITAL IMAGE PROCESSING UNIT III

Prof. Vidya Manian Dept. of Electrical and Comptuer Engineering

A Division of Sun Chemical Corporation. Unsharp Masking How to Make Your Images Pop!

AN EFFICIENT IMAGE ENHANCEMENT ALGORITHM FOR SONAR DATA

Fourier Transforms and the Frequency Domain

TIRF, geometric operators

On the evaluation of edge preserving smoothing filter

Head, IICT, Indus University, India

Filtering in the spatial domain (Spatial Filtering)

Digital Image Processing

Image analysis. CS/CME/BioE/Biophys/BMI 279 Oct. 31 and Nov. 2, 2017 Ron Dror

Images and Filters. EE/CSE 576 Linda Shapiro

Enhancement Techniques for True Color Images in Spatial Domain

DIGITAL IMAGE DE-NOISING FILTERS A COMPREHENSIVE STUDY

Image analysis. CS/CME/BioE/Biophys/BMI 279 Oct. 31 and Nov. 2, 2017 Ron Dror

PARAMETRIC ANALYSIS OF IMAGE ENHANCEMENT TECHNIQUES

Lecture No Image Filtering (course: Computer Vision)

An Efficient Noise Removing Technique Using Mdbut Filter in Images

Image Processing COS 426

SURVEILLANCE SYSTEMS WITH AUTOMATIC RESTORATION OF LINEAR MOTION AND OUT-OF-FOCUS BLURRED IMAGES. Received August 2008; accepted October 2008

Motivation: Image denoising. How can we reduce noise in a photograph?

Image Enhancement in the Spatial Domain

Computing for Engineers in Python

Chapter 2 Image Enhancement in the Spatial Domain

ABSTRACT I. INTRODUCTION

[Kaur, 2(8): August, 2013] ISSN: Impact Factor: INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY

LINEAR AND NONLINEAR FILTER FOR IMAGE PROCESSING USING MATLAB S IMAGE PROCESSING TOOLBOX

Image filtering, image operations. Jana Kosecka

Motivation: Image denoising. How can we reduce noise in a photograph?

Vision Review: Image Processing. Course web page:

CS534 Introduction to Computer Vision. Linear Filters. Ahmed Elgammal Dept. of Computer Science Rutgers University

Achim J. Lilienthal Mobile Robotics and Olfaction Lab, AASS, Örebro University

June 30 th, 2008 Lesson notes taken from professor Hongmei Zhu class.

Filtering. Image Enhancement Spatial and Frequency Based

Filters. Materials from Prof. Klaus Mueller

Transforms and Frequency Filtering

Digital Image Processing

Adaptive Fingerprint Binarization by Frequency Domain Analysis

FPGA IMPLEMENTATION OF RSEPD TECHNIQUE BASED IMPULSE NOISE REMOVAL

TDI2131 Digital Image Processing (Week 4) Tutorial 3

Matlab (see Homework 1: Intro to Matlab) Linear Filters (Reading: 7.1, ) Correlation. Convolution. Linear Filtering (warm-up slide) R ij

Midterm Examination CS 534: Computational Photography

Image Processing for feature extraction

Chapter 6. [6]Preprocessing

Image Enhancement contd. An example of low pass filters is:

A Histogram based Algorithm for Denoising Images Corrupted with Impulse Noise

Keywords Fuzzy Logic, ANN, Histogram Equalization, Spatial Averaging, High Boost filtering, MSE, RMSE, SNR, PSNR.

CEE598 - Visual Sensing for Civil Infrastructure Eng. & Mgmt.

8.2 IMAGE PROCESSING VERSUS IMAGE ANALYSIS Image processing: The collection of routines and

קורס גרפיקה ממוחשבת 2008 סמסטר ב' Image Processing 1 חלק מהשקפים מעובדים משקפים של פרדו דוראנד, טומס פנקהאוסר ודניאל כהן-אור

UM-Based Image Enhancement in Low-Light Situations

Last Lecture. photomatix.com

New Spatial Filters for Image Enhancement and Noise Removal

Index Terms: edge-preserving filter, Bilateral filter, exploratory data model, Image Enhancement, Unsharp Masking

Digital Image Processing. Digital Image Fundamentals II 12 th June, 2017

Computer Vision for HCI. Noise Removal. Noise in Images

Image Processing. What is an image? קורס גרפיקה ממוחשבת 2008 סמסטר ב' Converting to digital form. Sampling and Reconstruction.

Image preprocessing in spatial domain

02/02/10. Image Filtering. Computer Vision CS 543 / ECE 549 University of Illinois. Derek Hoiem

Review and Analysis of Image Enhancement Techniques

Image Manipulation: Filters and Convolutions

OFFSET AND NOISE COMPENSATION

Computer Vision, Lecture 3

Multimedia Systems Giorgio Leonardi A.A Lectures 14-16: Raster images processing and filters

ENEE408G Multimedia Signal Processing

Preparing Remote Sensing Data for Natural Resources Mapping (image enhancement, rectifications )

Image Smoothening and Sharpening using Frequency Domain Filtering Technique

Performance Analysis of Average and Median Filters for De noising Of Digital Images.

Design of Practical Color Filter Array Interpolation Algorithms for Cameras, Part 2

Maine Day in May. 54 Chapter 2: Painterly Techniques for Non-Painters

TDI2131 Digital Image Processing

8. Lecture. Image restoration: Fourier domain

Algorithm for Image Processing Using Improved Median Filter and Comparison of Mean, Median and Improved Median Filter

Transcription:

Blekinge Institute of Technology Research Report 2004:08 Chrominance Assisted Sharpening of Images Andreas Nilsson Department of Signal Processing School of Engineering Blekinge Institute of Technology

Chrominance Assisted Sharpening of Images Andreas Nilsson (730325-3335) Department of Signal Processing School of Engineering Blekinge Institute of Technology Ronneby, Sweden 2004-08-17

Abstract This paper presents a method to sharpen an image with less amplification of the noise than conventional methods. The proposed design method uses unsharp filtering for the sharpening part and differences in the chrominance data to decide if sharpening should be performed or not. The basic idea is to only do the sharpening if there are any changes in the chrominance information. Examples are included of sharpened images with less amplification of the noise. 1

Introduction Different image sharpening methods have been known for a long time and can be studied in [1], [2] and [3]. In the literature the sharpening is often performed with either a high pass filter or by an unsharp filter or mask. One problem with these methods is that the sharpening is applied to the entire image. This can result in that not only contures but also unwanted noise gets amplified. In an image there is often areas which are not desired to sharpen e.g. clothes, blue-sky, buildings etc. Generally textures like these have the same color over the whole region and even if there is a change it is of a low-frequency characteristic. This makes it possible to use the chrominance components in the decision step. To be able to perform this the image, which is commonly represented in RGB (red-green-blue), must be represented in a format which separate the luminance and chrominance components, e.g. Y C b C r where Y represents the luminance, C b the chrominance blue and C r the chrominance red. In this format the sharpening can be performed on the luminance, Y, and the decision to sharp or not can be performed using the chrominance, C b and C r, where sharpening is excluded when there is little variation in the chrominance components. In Section 1, a description of different sharpening methods is given. Section 2 presents the method using the chrominance data to decide if sharpening should be performed or not. Section 3 presents the results including image examples. The report ends with a conclution in section 4. 2

1 Sharpening To perform sharpening there are two basic methods: highpass filtering and unsharp filtering which are presented in [1], [2] and [3]. These can be used in different ways presented in this section. 1.1 Highpass Filtering The simplest way to sharpen an image is to use a high-pass filter, which emphasizes the high frequencies components and reduces the low frequencies. Some typical examples of highpass masks used for contrast enhancements are shown in Eq. 1 h HP = h HP = h HP = 0 1 0 1 5 1 0 1 0 1 1 1 1 9 1 1 1 1 1 2 1 2 5 2 1 2 1 (1) A drawback of this approach is that the results from areas with constant or slow variation will be very small; this will result in an image with enhanced edges over a rather dark area. To preserve the low-frequencies parts of the image it is significantly better to us an unsharp filter. 1.2 Unsharp Filtering Unsharp filtering or unsharp masking has been used by photographic artists and printing industries for a long time. The method is a simple way to enhance the edges and other high frequency components but with the possibility to keep the basic low frequencies in the image. This is performed via a procedure where the original image first is blurred, unsharped, and then subtracted from the original to produce an edge image g(x, y). This can be expressed as g(x, y) = f(x, y) f blurr (x, y) (2) where f(x, y) is the original image and f blurr (x, y) is the blurred or unsharped version of f(x, y). This image can then be used for sharpening if it is added back into the original image f sharp (x, y) = f(x, y) + m g(x, y) (3) 3

where m is a scaling factor that decides the amounts of sharpening. The blurred version f blurr (x, y) of the original image can be generated by convolution with a 2-D Gaussian filter where the discrete Fourier transform of it is expressed as G(n 1, n 2 ) = 1 n 2 1 +n2 2 2πσ 2 e 2σ 2 (4) and σ represents the Gaussian filters spread factor. In Eq. 5 a suitable integervalued convolution kernel that approximates a Gaussian with a σ of 1.0 is shown. 1 4 7 4 1 h gauss = 1 4 16 26 16 4 273 7 26 41 26 7 (5) 4 16 26 16 4 1 4 7 4 1 The blurred version f blurr (x, y) can also be generated by using a common lowpass (LP) filter as illustrated in Eq 6. h LP = 1 9 h LP = 1 10 h LP = 1 16 1 2 1 1 2 1 2 4 2 1 2 1 (6) which are common for image enhancement. A more common way to perform the unsharp filtering is to generate g(x, y) by using the negative Laplacian, the second-order derivative, operator. The Laplacian L(x, y) of an image with pixel intensity values I(x, y) is given by: L(x, y) = 2 I 2 x 2 + 2 I 2 y 2 (7) Three commonly used discrete negative Laplacian kernels are shown in Eq. 8. 0 1 0 h Laplacian = 1 1 4 1 0 1 0 h Laplacian = 1 4 1 8 1 (8)

h Laplacian = 1 1 2 1 2 4 2 1 2 1 By using these negative Laplacian filters g(x, y) can be generated direct by convolution and Eq. 3 can then be used. 2 Chrominance Assisted Sharpening In a sharpened image the high frequencies of the image are amplified. This results in an image where contours and edges will be more prominent. The problem is that this, most likely, also will amplify unwanted high-frequency noise in the image. In this section the sharpening is only applied on the luminance component for pixels, where there are difference in the chrominance component compared to the surrounding pixels. 2.1 Chrominance Variations To decide if there is a variation in the chrominance some measure must be used. If the two chrominance components C b and C r are seen as two axes in a 2-D plane the chrominance distance between two pixels can be calculated as d = (C b1 C b2 ) 2 + (C r1 C r2 ) 2 (9) The distance d shall then be compared with some threshold k. This distance calculation is performed for the nearby neighbors for each pixel in the image and the maximum is chosen to be compared with k. To decide the threshold k some empirical tests must be performed on the camera which generates the images, due to the different characteristics between cameras. Thus, different cameras generate different amount of noise in the luminance and chrominance components. 2.2 Pre-filtering Depending on the characteristics of different cameras the amount of noise in luminance and chrominance components will differ. For cheaper cameras the amount of noise is essential higher in the chrominance components and since the decision is based on variation in the chrominance components it can be profitable to do some pre-filtering before the decision is taken. This will not influence the subjective result since the filtered components are only used for the decisionmaking. 5

The pre-filtering can be done by both a common LP-filter like the mean filter; the filter mask is shown in Eq. 10, or a median filter. h LPmean = 1 9 (10) The mean filter attenuates high frequencies present in the image, but if the noise includes single very unrepresentative value this can significantly affect the mean value of all the pixels in its neighborhood. The filter will also blur edges, which can affect the result of the decision to sharpen or not. An alternative is to use a median filter, which is more robust for single very unrepresentative values of the components. The median filter sorts the components in the surrounding of the chosen component and then replaces it with the middle component value. As a result of this the filter will not create new unrealistic values when an edge is filtered. What to choose depends on the type of noise the image includes. 2.3 Implementation This section shows a flow chart of the method, see Figure 1. 3 Results In this section the results on images from three different cameras are shown. The unsharp filtering method is used for all three cases. The different cameras demand different setups of pre-filtering strength on the sharpening, f, and threshold, k. Also, both common LP filtering as negative Laplacian filters have been used to generate g(x, y). The LP filter used is shown in Eq. 11 h LP = 1 25 1 1 1 1 1 1 1 1 1 1 and the negative Laplacian filter usedd is shown in Eq. 12. h Laplacian = 1 1 8 1 (12) (11) Image one, Vacation, is taken with a Canon Ixus 400, the second, Flowergirl, with a Canon S10 and the third, House, with a SonyEricsson 610 mobile phone. It has not been possible to calculate any PSNR so just subjective evaluation has been used to tune in the different parameters in the method. 6

START Read new pixels luminance and chrominance components (Y, Cb and Cr) Pre-filtering chrominance NO YES Pre-filtering chrominance Calculate d d > k NO YES Sharpening is applied on the luminance component and the original chrominanse is used. The original pixel used. Figure 1: Flowchart describing the selective sharpening filter. 3.1 Vacation The original image is shown in Figure 2. The sharpened image with both LP filtering and Laplacian filtering is shown in Figure 3. In Figure 4 the mask that illustrates where the sharpened pixels are used is shown. The final results are shown in Figure 5. No pre-filtering on chrominance is used. 3.2 Flowergirl The original image is shown in Figure 6. The sharpened image with both LP filtering and Laplacian filtering is shown in Figure 7. In Figure 8 the mask that illustrate where the sharpened pixels are used is shown. The final results are 7

shown in Figure 9. No pre-filtering on chrominance is used. 3.3 House The original image is shown in Figure 10. The sharpened image with both LP filtering and Laplacian filtering is shown in Figure 11. In Figure 12 the mask that illustrate where the sharpended pixels are used is shown. The final results are shown in Figure 13. LP pre-filtering has been performed on the chrominance components with the same filter, h LP, used to generate g(x, y). 4 Conclusion This report presents a method where chrominance components are used to decide if sharpening is required or not. The sharpening has been performed in a general way using unmask filtering with lowpass filters or negative Laplacian filters. The decision-making has been performed by calculating the variation in chrominance component. No sharpening has been performed if there is little variation in the chrominance components. In section 3 the results from three different cameras are presented. In the two first examples, Vacation and Flowergirl no pre-filtering of the chrominance components has been performed, depending on the low noise contribution from the cameras. In the third example, House, the result from a much cheaper camera is shown where it is a much higher noise contribution. Therefore, pre-filtering has been used in this example. The threshold k is different for all three cases to avoid amplifying of the noise. If k would be reduced a much higher amount of noise would be included in the image which would lower the subjective result. In all examples it can be seen that with the chrominance assisted sharpening method a smoother sharpening is performed and no or little noise is amplified. The evaluation and parameter tuning has been performed based on subjective opinions. References [1] William K. Pratt, Digital Image Processing: PIKS Inside, Third Edition, 2001, John Wiley & Sons, Inc. [2] Jae S. Lim, Two-Dimensional SIGNAL and IMAGE Processoing, 1990, Prentice Hall PTR. [3] Rafael C Gonzales, Richard E. Woods, Digital Image Processing, 1993, Addisio-Wesley Publishing Company 8

Vacation Figure 2: Original image. Figure 3: Left Sharpening with LP filter. Right Sharpening with Laplacian filter. 9

Figure 4: The resulting Mask, showing where d > k. Figure 5: Left Assisted Sharpening with LP filter. Right Assisted Sharpening with Laplacian filter. 10

Flowergirl Figure 6: Original image. Figure 7: Left Sharpening with LP filter. Right Sharpening with Laplacian filter. 11

Figure 8: The resulting Mask, showing where d > k. Figure 9: Left Assisted Sharpening with LP filter. Right Assisted Sharpening with Laplacian filter. 12

House Figure 10: Original image. Figure 11: Left Sharpening with LP filter. Right Sharpening with Laplacian filter. 13

Figure 12: The resulting Mask, showing where d > k. Figure 13: Left Assisted Sharpening with LP filter. Right Assisted Sharpening with Laplacian filter. 14

Chrominance Assisted Sharpening of Images Andreas Nilsson ISSN 1103-1581 ISRN BTH-RES--08/04--SE Copyright 2004 by the author All rights reserved Printed by Kaserntryckeriet AB, Karlskrona 2004