Apparatus and method for adaptively reducing noise in a noisy input image signal
-
0Associated
Cases -
0Associated
Defendants -
0Accused
Products -
90Forward
Citations -
0
Petitions -
8
Assignments
First Claim
1. An adaptive apparatus for spatially reducing noise in a noisy input image signal comprising:
- a low-pass filter receiving said noisy input image signal to generate a noisy low-spatial frequency image signal;
a first pixel-based serial-to-parallel converter receiving the noisy low-spatial frequency image signal to generate a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to predetermined pixels-window characteristics;
a pixel-based local window segmentation processor comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels;
a counter generating a selected-pixels count signal;
a second pixel-based serial-to-parallel converter receiving the noisy input image signal to generate a group of noisy input image parallel signals associated with the locally considered pixel and according to the predetermined pixels-window characteristics;
a mean estimator combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel;
a minimum-mean-square-error filter receiving the noisy input image parallel signals, the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a noise-filtered output image signal according to an input noise statistic signal.
8 Assignments
0 Petitions

Accused Products

Abstract
The basic configuration of Single local Adaptive Window Spatial Noise Reducer (SAW-SNR) is based on a preliminary de-noising low-pass filter followed by homogenous region segmentation to the considered pixel in a given local window. The configuration is composed also of an adaptive local mean estimator, an adaptive local statistic estimator which is preferably an economic standard deviation (SD) estimator and finally, a minimum-mean-square-error (MMSE) based de-noising technique. The proposed segmentation configuration outperforms existing spatial noise reducers in term of subjective and objective performances, in term of edge preservation, noise reduction in both homogenous regions or picture edges and Peak Signal to noise Ratio (PSNR). A second configuration in the form of a Parallel Multiple local Adaptive Window Spatial Noise Reducer (Parallel M-AW-SNR), is a combination of several basic configurations which implements different segmented windows. The M-AW-SNR, which is the less complex configuration for multiple spatial noise reducers, reduces further residual noise as compared to the basic configuration. A third configuration combines the basic configuration of SAW-SNR with a controllable noise variance estimator. This generic configuration allows an adaptive local control of noise reduction level, which can be useful for some correlated noise such as ringing noise in DCT-based decompressed images or cross-luminance noise in composite decoded images.
103 Citations
View as Search Results
Using statistics to locate signals in noise | ||
Patent #
US 7,869,980 B2
Filed 11/03/2005
|
Current Assignee
International Business Machines Corporation
|
Sponsoring Entity
International Business Machines Corporation
|
IMAGE DENOISING METHOD | ||
Patent #
US 20110069902A1
Filed 07/02/2010
|
Current Assignee
PixArt Imaging Incorporated
|
Sponsoring Entity
PixArt Imaging Incorporated
|
VIDEO SIGNAL PROCESSING APPARATUS, ENHANCEMENT GAIN GENERATING METHOD, AND ENHANCED GAIN GENERATING PROGRAM | ||
Patent #
US 20110037901A1
Filed 07/29/2010
|
Current Assignee
Sony Corporation
|
Sponsoring Entity
Sony Corporation
|
IMAGE PROCESSING DEVICE AND METHOD | ||
Patent #
US 20110026837A1
Filed 07/29/2010
|
Current Assignee
Casio Computer Company
|
Sponsoring Entity
Casio Computer Company
|
VISUAL PROCESSING DEVICE, VISUAL PROCESSING METHOD, PROGRAM, DISPLAY DEVICE, AND INTEGRATED CIRCUIT | ||
Patent #
US 20110109765A1
Filed 01/13/2011
|
Current Assignee
Panasonic Intellectual Property Corporation of America
|
Sponsoring Entity
Panasonic Intellectual Property Corporation of America
|
Method and apparatus for short range motion adaptive noise reduction | ||
Patent #
US 7,961,972 B2
Filed 05/01/2007
|
Current Assignee
Intel Corporation
|
Sponsoring Entity
Intel Corporation
|
Video processing system and method with recursive and non-recursive filtering in separate high frequency bands | ||
Patent #
US 7,969,510 B2
Filed 07/14/2003
|
Current Assignee
Snell Limited
|
Sponsoring Entity
Snell Limited
|
VISUAL PROCESSING APPARATUS, VISUAL PROCESSING METHOD, PROGRAM, RECORDING MEDIUM, DISPLAY DEVICE, AND INTEGRATED CIRCUIT | ||
Patent #
US 20110122148A1
Filed 01/31/2011
|
Current Assignee
Panasonic Intellectual Property Corporation of America
|
Sponsoring Entity
Panasonic Intellectual Property Corporation of America
|
Techniques for Reducing Noise While Preserving Contrast in an Image | ||
Patent #
US 20110090240A1
Filed 06/05/2009
|
Current Assignee
FotoNation Limited
|
Sponsoring Entity
FotoNation Limited
|
Digitized image stabilization using energy analysis method | ||
Patent #
US 7,961,966 B2
Filed 01/04/2005
|
Current Assignee
Etron Technology Inc
|
Sponsoring Entity
Etron Technology Inc
|
VISUAL PROCESSING DEVICE, DISPLAY DEVICE, VISUAL PROCESSING METHOD, PROGRAM, AND INTEGRATED CIRCUIT | ||
Patent #
US 20110085085A1
Filed 12/20/2010
|
Current Assignee
Panasonic Intellectual Property Corporation of America
|
Sponsoring Entity
Panasonic Intellectual Property Corporation of America
|
Method And System For Increasing Signal-To-Noise Ratio | ||
Patent #
US 20110170796A1
Filed 01/29/2007
|
Current Assignee
Canadian Space Agency
|
Sponsoring Entity
Canadian Space Agency
|
Using statistics to locate signals in noise | ||
Patent #
US 8,024,156 B2
Filed 06/29/2010
|
Current Assignee
International Business Machines Corporation
|
Sponsoring Entity
International Business Machines Corporation
|
Information processing apparatus and information processing method | ||
Patent #
US 7,734,095 B2
Filed 10/18/2005
|
Current Assignee
Canon Ayutthaya Limited
|
Sponsoring Entity
Canon Ayutthaya Limited
|
Method and system for noise reduction in an image | ||
Patent #
US 7,668,396 B2
Filed 09/29/2003
|
Current Assignee
VIXS Systems Incorporated
|
Sponsoring Entity
VIXS Systems Incorporated
|
Multi-touch input discrimination | ||
Patent #
US 7,855,718 B2
Filed 01/03/2007
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
Apple Inc.
|
IMAGE PROCESSING APPARATUS AND METHOD OF PROCESSING IMAGE | ||
Patent #
US 20100066868A1
Filed 09/11/2009
|
Current Assignee
Canon Ayutthaya Limited
|
Sponsoring Entity
Canon Ayutthaya Limited
|
SYSTEM AND METHOD FOR MULTI-SCALE SIGMA FILTERING USING QUADRATURE MIRROR FILTERS | ||
Patent #
US 20090214131A1
Filed 02/27/2008
|
Current Assignee
Xerox Corporation
|
Sponsoring Entity
Xerox Corporation
|
Image enhancement method using local illumination correction | ||
Patent #
US 7,590,303 B2
Filed 09/29/2005
|
Current Assignee
Samsung Electronics Co. Ltd.
|
Sponsoring Entity
Samsung Electronics Co. Ltd.
|
METHOD AND SYSTEM FOR NOISE LEVEL DETECTION IN IMAGE DATA | ||
Patent #
US 20090257674A1
Filed 04/11/2008
|
Current Assignee
Xerox Corporation
|
Sponsoring Entity
Xerox Corporation
|
Method, apparatus, and computer product for eliminating noise | ||
Patent #
US 20080024668A1
Filed 12/28/2006
|
Current Assignee
Fujitsu Limited
|
Sponsoring Entity
Fujitsu Limited
|
Method of repairing scratches in digital images | ||
Patent #
US 7,324,704 B2
Filed 09/11/2003
|
Current Assignee
Primax Electronics Limited
|
Sponsoring Entity
Primax Electronics Limited
|
Methods and apparatuses providing noise reduction while preserving edges for imagers | ||
Patent #
US 20080043124A1
Filed 09/17/2007
|
Current Assignee
Aptina Imaging Corporation
|
Sponsoring Entity
Aptina Imaging Corporation
|
Multi-Touch Input Discrimination | ||
Patent #
US 20080158185A1
Filed 05/31/2007
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
-
|
Apparatus and method for adaptive 3D noise reduction | ||
Patent #
US 7,418,149 B2
Filed 07/29/2005
|
Current Assignee
SENSIO TECHNOLOGIES INC.
|
Sponsoring Entity
ALGOLITH INC.
|
Enhancing text-like edges in digital images | ||
Patent #
US 7,433,535 B2
Filed 09/30/2003
|
Current Assignee
Hewlett-Packard Development Company L.P.
|
Sponsoring Entity
Hewlett-Packard Development Company L.P.
|
Method and apparatus for short range motion adaptive noise reduction | ||
Patent #
US 20080273119A1
Filed 05/01/2007
|
Current Assignee
Intel Corporation
|
Sponsoring Entity
Intel Corporation
|
DECODING WITH EMBEDDED DENOISING | ||
Patent #
US 20080285655A1
Filed 05/16/2008
|
Current Assignee
Hong Kong University of Science and Technology
|
Sponsoring Entity
Hong Kong University of Science and Technology
|
Apparatus and methods for image restoration | ||
Patent #
US 20080317375A1
Filed 06/21/2007
|
Current Assignee
SOUTHERN MISSISSIPPI THE UNIVERSITY OF
|
Sponsoring Entity
SOUTHERN MISSISSIPPI THE UNIVERSITY OF
|
Image enhancement method using local illumination correction | ||
Patent #
US 20070071350A1
Filed 09/29/2005
|
Current Assignee
Samsung Electronics Co. Ltd.
|
Sponsoring Entity
Samsung Electronics Co. Ltd.
|
Method and apparatus for image noise reduction | ||
Patent #
US 20070133893A1
Filed 12/14/2005
|
Current Assignee
Aptina Imaging Corporation
|
Sponsoring Entity
Aptina Imaging Corporation
|
Fast edge-oriented image interpolation algorithm | ||
Patent #
US 7,286,721 B2
Filed 09/11/2003
|
Current Assignee
CHEVRON TECHNOLOGY VENTURES LLC
|
Sponsoring Entity
CHEVRON TECHNOLOGY VENTURES LLC
|
Method and system for image processing for structured light profiling of a part | ||
Patent #
US 7,302,109 B2
Filed 08/28/2003
|
Current Assignee
General Electric Company
|
Sponsoring Entity
General Electric Company
|
Predicting ringing artifacts in digital images | ||
Patent #
US 6,845,180 B2
Filed 03/16/2001
|
Current Assignee
Sharp Laboratories of America Incorporated
|
Sponsoring Entity
Sharp Laboratories of America Incorporated
|
Video deblocking method and apparatus | ||
Patent #
US 6,922,492 B2
Filed 12/27/2002
|
Current Assignee
Google Technology Holdings LLC
|
Sponsoring Entity
Motorola Inc.
|
Apparatus and method for adaptive 3D noise reduction | ||
Patent #
US 20060056724A1
Filed 07/29/2005
|
Current Assignee
SENSIO TECHNOLOGIES INC.
|
Sponsoring Entity
SENSIO TECHNOLOGIES INC.
|
Information processing apparatus and information processing method | ||
Patent #
US 20060083443A1
Filed 10/18/2005
|
Current Assignee
Canon Ayutthaya Limited
|
Sponsoring Entity
Canon Ayutthaya Limited
|
Spatio-temporal filter unit and image display apparatus comprising such a spatio-temporal filter unit | ||
Patent #
US 7,034,892 B2
Filed 01/23/2002
|
Current Assignee
Koninklijke Philips N.V.
|
Sponsoring Entity
Koninklijke Philips N.V.
|
Video processing | ||
Patent #
US 20060109383A1
Filed 07/14/2003
|
Current Assignee
Snell Limited
|
Sponsoring Entity
Snell Limited
|
Digitized image stabilization using energy analysis method | ||
Patent #
US 20060146139A1
Filed 01/04/2005
|
Current Assignee
Etron Technology Inc
|
Sponsoring Entity
Etron Technology Inc
|
System and method for sharpening vector-valued digital images | ||
Patent #
US 20060182364A1
Filed 02/16/2005
|
Current Assignee
SOZOTEK INC
|
Sponsoring Entity
SOZOTEK INC
|
Video signal encoding device and video signal encoding method | ||
Patent #
US 20060209949A1
Filed 12/25/2003
|
Current Assignee
Mitsubishi Electric Corporation
|
Sponsoring Entity
Mitsubishi Electric Corporation
|
Method of reparing scratches in digital images | ||
Patent #
US 20050058366A1
Filed 09/11/2003
|
Current Assignee
Primax Electronics Limited
|
Sponsoring Entity
Primax Electronics Limited
|
Method and system for image processing for structured light profiling of a part | ||
Patent #
US 20050046872A1
Filed 08/28/2003
|
Current Assignee
General Electric Company
|
Sponsoring Entity
General Electric Company
|
Method and system for noise reduction in an image | ||
Patent #
US 20050069221A1
Filed 09/29/2003
|
Current Assignee
VIXS Systems Incorporated
|
Sponsoring Entity
VIXS Systems Incorporated
|
Enhancing text-like edges in digital images | ||
Patent #
US 20050069217A1
Filed 09/30/2003
|
Current Assignee
Hewlett-Packard Development Company L.P.
|
Sponsoring Entity
Hewlett-Packard Development Company L.P.
|
Fast edge-oriented image interpolation algorithm | ||
Patent #
US 20050058371A1
Filed 09/11/2003
|
Current Assignee
CHEVRON TECHNOLOGY VENTURES LLC
|
Sponsoring Entity
CHEVRON TECHNOLOGY VENTURES LLC
|
Method and system for noise estimation from video sequence | ||
Patent #
US 20050107982A1
Filed 11/17/2003
|
Current Assignee
Eastman Kodak Company
|
Sponsoring Entity
Eastman Kodak Company
|
Video deblocking method and apparatus | ||
Patent #
US 20040126034A1
Filed 12/27/2002
|
Current Assignee
Google Technology Holdings LLC
|
Sponsoring Entity
Google Technology Holdings LLC
|
Spatio-temporal filter unit and image display apparatus comprising such a spatio-temporal filter unit | ||
Patent #
US 20020101543A1
Filed 01/23/2002
|
Current Assignee
Koninklijke Philips N.V.
|
Sponsoring Entity
Koninklijke Philips N.V.
|
Predicting ringing artifacts in digital images | ||
Patent #
US 20020131647A1
Filed 03/16/2001
|
Current Assignee
Sharp Laboratories of America Incorporated
|
Sponsoring Entity
Sharp Laboratories of America Incorporated
|
Noise reduction method, medium, and system | ||
Patent #
US 20070296871A1
Filed 06/04/2007
|
Current Assignee
Samsung Electronics Co. Ltd.
|
Sponsoring Entity
Samsung Electronics Co. Ltd.
|
Video signal encoding device and video signal encoding method | ||
Patent #
US 8,094,715 B2
Filed 12/25/2003
|
Current Assignee
Mitsubishi Electric Corporation
|
Sponsoring Entity
Mitsubishi Electric Corporation
|
SYSTEM AND METHOD FOR REDUCING NOISE IN AN IMAGE | ||
Patent #
US 20120002896A1
Filed 06/30/2010
|
Current Assignee
Samsung Electronics Co. Ltd.
|
Sponsoring Entity
Samsung Electronics Co. Ltd.
|
Apparatus and methods for image restoration | ||
Patent #
US 8,073,277 B2
Filed 06/21/2007
|
Current Assignee
SOUTHERN MISSISSIPPI THE UNIVERSITY OF
|
Sponsoring Entity
SOUTHERN MISSISSIPPI THE UNIVERSITY OF
|
Multi-touch input discrimination | ||
Patent #
US 8,130,203 B2
Filed 05/31/2007
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
Apple Inc.
|
System and method for multi-scale sigma filtering using quadrature mirror filters | ||
Patent #
US 8,180,169 B2
Filed 02/27/2008
|
Current Assignee
Xerox Corporation
|
Sponsoring Entity
Xerox Corporation
|
Multi-touch input discrimination | ||
Patent #
US 8,243,041 B2
Filed 01/18/2012
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
Apple Inc.
|
ANISOTROPIC DENOISING METHOD | ||
Patent #
US 20120224784A1
Filed 03/01/2011
|
Current Assignee
FotoNation Limited
|
Sponsoring Entity
FotoNation Limited
|
Image processing apparatus and method of processing image for reducing noise of the image | ||
Patent #
US 8,300,120 B2
Filed 09/11/2009
|
Current Assignee
Canon Ayutthaya Limited
|
Sponsoring Entity
Canon Ayutthaya Limited
|
Visual processing device, display device, visual processing method, program, and integrated circuit | ||
Patent #
US 8,311,357 B2
Filed 12/20/2010
|
Current Assignee
Panasonic Intellectual Property Corporation of America
|
Sponsoring Entity
Panasonic Corporation
|
Image processing apparatus and method of processing image for reducing noise of the image | ||
Patent #
US 8,345,127 B2
Filed 06/20/2012
|
Current Assignee
Canon Ayutthaya Limited
|
Sponsoring Entity
Canon Ayutthaya Limited
|
Method and system for increasing signal-to-noise ratio | ||
Patent #
US 8,358,866 B2
Filed 01/29/2007
|
Current Assignee
Canadian Space Agency
|
Sponsoring Entity
Canadian Space Agency
|
Optimal denoising for video coding | ||
Patent #
US 8,369,417 B2
Filed 05/18/2007
|
Current Assignee
Hong Kong University of Science and Technology
|
Sponsoring Entity
Hong Kong University of Science and Technology
|
Multi-touch input discrimination | ||
Patent #
US 8,384,684 B2
Filed 12/10/2010
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
Apple Inc.
|
Visual processing apparatus, visual processing method, program, recording medium, display device, and integrated circuit | ||
Patent #
US 8,401,324 B2
Filed 01/31/2011
|
Current Assignee
Panasonic Intellectual Property Corporation of America
|
Sponsoring Entity
Panasonic Corporation
|
Visual processing device, visual processing method, program, display device, and integrated circuit | ||
Patent #
US 8,406,547 B2
Filed 01/13/2011
|
Current Assignee
Panasonic Intellectual Property Corporation of America
|
Sponsoring Entity
Panasonic Corporation
|
GENERAL BANDING AND CODEC BANDING ARTIFACT REMOVAL | ||
Patent #
US 20130128122A1
Filed 11/23/2011
|
Current Assignee
Stmicroelectronics Asia Pacific Pte Limited
|
Sponsoring Entity
Stmicroelectronics Asia Pacific Pte Limited
|
Noise reduction method, medium, and system | ||
Patent #
US 8,467,003 B2
Filed 06/04/2007
|
Current Assignee
Samsung Electronics Co. Ltd.
|
Sponsoring Entity
Samsung Electronics Co. Ltd.
|
System and method for reducing noise in an image | ||
Patent #
US 8,472,724 B2
Filed 06/30/2010
|
Current Assignee
Samsung Electronics Co. Ltd.
|
Sponsoring Entity
Samsung Electronics Co. Ltd.
|
SYSTEMS, METHODS, AND APPARATUS FOR DETECTING IRREGULAR SENSOR SIGNAL NOISE | ||
Patent #
US 20130197849A1
Filed 10/11/2010
|
Current Assignee
General Electric Company
|
Sponsoring Entity
General Electric Company
|
Multi-touch input discrimination | ||
Patent #
US 8,531,425 B2
Filed 07/27/2012
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
Apple Inc.
|
Method and system for noise level detection in image data | ||
Patent #
US 8,208,753 B2
Filed 04/11/2008
|
Current Assignee
Xerox Corporation
|
Sponsoring Entity
Xerox Corporation
|
Multi-touch input discrimination | ||
Patent #
US 8,542,210 B2
Filed 02/15/2013
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
Apple Inc.
|
Video signal processing apparatus, enhancement gain generating method, and enhanced gain generating program | ||
Patent #
US 8,634,027 B2
Filed 07/29/2010
|
Current Assignee
Sony Corporation
|
Sponsoring Entity
Sony Corporation
|
METHOD AND SYSTEM FOR SUBSTANTIALLY REMOVING DOT NOISE | ||
Patent #
US 20140056536A1
Filed 08/27/2012
|
Current Assignee
Toshiba Medical Systems Corporation, Toshiba Corporation
|
Sponsoring Entity
Toshiba Medical Systems Corporation, Toshiba Corporation
|
Techniques for reducing noise while preserving contrast in an image | ||
Patent #
US 8,723,879 B2
Filed 06/05/2009
|
Current Assignee
FotoNation Limited
|
Sponsoring Entity
FotoNation Limited
|
General banding and codec banding artifact removal | ||
Patent #
US 8,731,323 B2
Filed 11/23/2011
|
Current Assignee
Stmicroelectronics Asia Pacific Pte Limited
|
Sponsoring Entity
Stmicroelectronics Asia Pacific Pte Limited
|
Multi-touch input discrimination | ||
Patent #
US 8,791,921 B2
Filed 08/19/2013
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
Apple Inc.
|
Decoding with embedded denoising | ||
Patent #
US 8,831,111 B2
Filed 05/16/2008
|
Current Assignee
Hong Kong University of Science and Technology
|
Sponsoring Entity
Hong Kong University of Science and Technology
|
Anisotropic denoising method | ||
Patent #
US 8,879,841 B2
Filed 03/01/2011
|
Current Assignee
FotoNation Limited
|
Sponsoring Entity
FotoNation Limited
|
Image denoising method | ||
Patent #
US 8,885,965 B2
Filed 07/02/2010
|
Current Assignee
PixArt Imaging Incorporated
|
Sponsoring Entity
PixArt Imaging Incorporated
|
Multi-touch input discrimination | ||
Patent #
US 9,024,906 B2
Filed 07/28/2014
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
Apple Inc.
|
IMAGE PROCESSING SYSTEM AND METHOD | ||
Patent #
US 20150235100A1
Filed 08/07/2013
|
Current Assignee
Foundation of Soongsil University-Industry Cooperation
|
Sponsoring Entity
Foundation of Soongsil University-Industry Cooperation
|
Multi-touch input discrimination | ||
Patent #
US 9,256,322 B2
Filed 03/25/2015
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
Apple Inc.
|
Image noise level detection and removal system and method | ||
Patent #
US 9,607,238 B2
Filed 08/07/2013
|
Current Assignee
Foundation of Soongsil University-Industry Cooperation
|
Sponsoring Entity
Foundation of Soongsil University-Industry Cooperation
|
SYSTEM AND METHOD FOR STORING HIERARCHICAL DATA RELATED TO AN IMAGE | ||
Patent #
US 20170249305A1
Filed 05/15/2017
|
Current Assignee
Pixia Corp.
|
Sponsoring Entity
Pixia Corp.
|
Multi-touch input discrimination | ||
Patent #
US 9,778,807 B2
Filed 12/30/2015
|
Current Assignee
Apple Inc.
|
Sponsoring Entity
Apple Inc.
|
System and method for storing hierarchical data related to an image | ||
Patent #
US 10,049,123 B2
Filed 05/15/2017
|
Current Assignee
Pixia Corp.
|
Sponsoring Entity
Pixia Corp.
|
System and method for storing points of polygons related to an image | ||
Patent #
US 10,311,098 B2
Filed 07/11/2018
|
Current Assignee
Pixia Corp.
|
Sponsoring Entity
Pixia Corp.
|
Image signal processing apparatus | ||
Patent #
US 5,499,111 A
Filed 08/04/1992
|
Current Assignee
Panasonic Communications Company Limited
|
Sponsoring Entity
Matsushita Graphic Communications Systems Inc.
|
Apparatus for producing a noise-reducded image signal from an input image signal | ||
Patent #
US 5,499,057 A
Filed 08/22/1994
|
Current Assignee
Sony Corporation
|
Sponsoring Entity
Sony Corporation
|
Noise reduction apparatus for video image | ||
Patent #
US 5,404,179 A
Filed 09/11/1992
|
Current Assignee
Matsushita Electric Industrial Company Limited
|
Sponsoring Entity
Matsushita Electric Industrial Company Limited
|
Apparatus and method for smoothing images | ||
Patent #
US 5,442,462 A
Filed 06/10/1992
|
Current Assignee
D.V.P. TECHNOLOGIES LTD. 18 HATZOUK STREET
|
Sponsoring Entity
D.V.P. TECHNOLOGIES LTD. 18 HATZOUK STREET
|
Video noise reducer with detection of both total and temporal picture changes | ||
Patent #
US 5,191,419 A
Filed 04/24/1991
|
Current Assignee
BTS BROADCAST TELEVISION SYSTEMS GMBH
|
Sponsoring Entity
BTS BROADCAST TELEVISION SYSTEMS GMBH
|
Time domain television noise reduction system | ||
Patent #
US 5,260,775 A
Filed 04/11/1991
|
Current Assignee
Quivanus Processing LLC
|
Sponsoring Entity
Farouda Yves C.
|
Image control apparatus for a television receiver | ||
Patent #
US 5,272,533 A
Filed 01/14/1992
|
Current Assignee
HITACHI VIDEO INFORMATION SYSTEM INC. A CORPORATION OF JAPAN, Hitachi Ltd.
|
Sponsoring Entity
HITACHI VIDEO INFORMATION SYSTEM INC. A CORPORATION OF JAPAN, Hitachi Ltd.
|
Method and apparatus for improving the signal to noise ratio of an image formed of an object hidden in or behind a semi-opaque random media | ||
Patent #
US 5,140,463 A
Filed 03/08/1990
|
Current Assignee
RESEARCH FOUNDATION OF CITY OF COLLEGE OF NEW YORK
|
Sponsoring Entity
Robert R. Alfano, Kwong M. Yoo
|
Noise reduction system for color television signal | ||
Patent #
US 5,161,018 A
Filed 04/29/1991
|
Current Assignee
Sony Corporation
|
Sponsoring Entity
Sony Corporation
|
Multiple measurement noise reducing system using space-variant filters | ||
Patent #
US 4,503,461 A
Filed 02/22/1983
|
Current Assignee
Board of Trustees of the Leland Stanford Junior University
|
Sponsoring Entity
Board of Trustees of the Leland Stanford Junior University
|
System for coring an image-representing signal | ||
Patent #
US 4,523,230 A
Filed 10/22/1984
|
Current Assignee
RCA Corporation
|
Sponsoring Entity
RCA Corporation
|
Adaptive image enhancement system | ||
Patent #
US 4,334,244 A
Filed 07/28/1980
|
Current Assignee
Hughes Electronics Corporation
|
Sponsoring Entity
MAGNAVOX GOVERNMENT AND INDUSTRIAL ELECTRONICS COMPANY
|
Noise reduction system | ||
Patent #
US 4,163,258 A
Filed 12/21/1976
|
Current Assignee
Sony Corporation
|
Sponsoring Entity
Sony Corporation
|
28 Claims
- 1. An adaptive apparatus for spatially reducing noise in a noisy input image signal comprising:
a low-pass filter receiving said noisy input image signal to generate a noisy low-spatial frequency image signal;
a first pixel-based serial-to-parallel converter receiving the noisy low-spatial frequency image signal to generate a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to predetermined pixels-window characteristics;
a pixel-based local window segmentation processor comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels;
a counter generating a selected-pixels count signal;
a second pixel-based serial-to-parallel converter receiving the noisy input image signal to generate a group of noisy input image parallel signals associated with the locally considered pixel and according to the predetermined pixels-window characteristics;
a mean estimator combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel;
a minimum-mean-square-error filter receiving the noisy input image parallel signals, the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a noise-filtered output image signal according to an input noise statistic signal. - View Dependent Claims (2, 3, 4, 5, 6)
wherein; g*(i,j;
x,y) is said values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with i=−
k, . . . ,+l;
j=−
m, . . . ,+n, said window having dimensions of (k+l)×
(m+n), k,l,m,n being appropriate positive integers;
g*(x,y) is said locally considered pixel value; and
T is a predetermined threshold value.
- 3. An apparatus according to
claim 2 , wherein a value of said mean pixel value signal is given by: - 4. An apparatus according to
claim 1 , wherein said input noise statistic signal is an input noise variance signal, said minimum-mean-square-error filter including:a variance calculator for generating an estimated variance signal from said noisy input image parallel signals, said segmented local window parallel signals, said selected-pixels count signal and said mean pixel value signal;
a weight calculator for generating a weight signal from the estimated variance signal and the input noise variance signal;
a first adder having an inverting input receiving said mean pixel value signal and a positive input receiving the noisy input signal associated with the locally considered pixel to generate a difference signal;
a multiplier receiving the weight signal and the difference signal to generate a weighted difference signal;
a second adder having a first positive input receiving the mean pixel value signal and a second positive input receiving the weighted difference signal to generate said noise-filtered output image signal.
- 5. An apparatus according to
claim 1 , wherein said input noise statistic signal is an input noise standard deviation signal, said minimum-mean-square-error filter including:a standard deviation calculator for generating an estimated standard deviation signal from said noisy input image parallel signals, said segmented local window parallel signals, said selected-pixels count signal and said mean pixel value signal;
a weight calculator for generating a weight signal from the estimated standard deviation signal and the input noise standard deviation signal;
a first adder having an inverting input receiving said mean pixel value signal and a positive input receiving the noisy input signal associated with the locally considered pixel to generate a difference signal;
a multiplier receiving the weight signal and the difference signal to generate a weighted difference signal;
a second adder having a first positive input receiving the mean pixel value signal and a second positive input receiving the weighted difference signal to generate said noise-filtered output image signal.
- 6. An apparatus according to
claim 5 , wherein values of said segmented local window parallel signals are given by:wherein; g*(i,j;
x,y) is said values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with i=−
k, . . . ,+l;
j=−
m, . . . ,+n, said window having dimensions of (k+l)×
(m+n), k,l,m,n being appropriate positive integers;
g*(x,y) is said locally considered pixel value;
T is a predetermined threshold value;
wherein said mean pixel value signal is defined by;
- 7. An adaptive apparatus for spatially reducing noise in a noisy input image signal comprising:
a plurality of parallel-connected adaptive spatial noise reducers each presenting a distinct set of predetermined pixels-window characteristics, each said noise reducer receiving said noisy input image signal to generate a corresponding pre-filtered output image signal and an averaging unit receiving each said pre-filtered output image signal at a corresponding positive input thereof to generate a noise-filtered output image signal, wherein each said noise reducer comprises;
a low-pass filter receiving said noisy input image signal to generate a noisy low-spatial frequency image signal;
a first pixel-based serial-to-parallel converter receiving the noisy low-spatial frequency image signal to generate a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to the set of predetermined pixels-window characteristics;
a pixel-based local window segmentation processor comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels;
a counter generating a selected-pixels count signal;
a second pixel-based serial-to-parallel converter receiving the noisy input image signal to generate a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics;
a mean estimator combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel;
a minimum-mean-square-error filter receiving the noisy input image parallel signals, the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate said pre-filtered output image signal according to an input noise statistic signal. - View Dependent Claims (9)
wherein; g*(i,j;
x,y) is said values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with i=−
k, . . . ,+l;
j=−
m, . . . ,+n, said window having dimensions of (k+l)×
(m+n), k,l,m,n being appropriate positive integers characterizing said distinct set of predetermined pixels-window characteristics;
g*(x,y) is said locally considered pixel value;
T is a predetermined threshold value; and
wherein a value of said mean pixel value signal is given by;
- 8. An adaptive apparatus for spatially reducing noise in a noisy input image signal comprising:
a low-pass filter receiving said noisy input image signal to generate a noisy low-spatial frequency image signal;
a plurality of parallel-connected adaptive spatial noise reducers each presenting a distinct set of predetermined pixels-window characteristics, each said noise reducer receiving said noisy low-spatial frequency image signal to generate a corresponding pre-filtered output image signal and an averaging unit receiving each said pre-filtered output image signal at a corresponding positive input thereof to generate a noise-filtered output image signal, wherein each said noise reducer comprises;
a first pixel-based serial-to-parallel converter receiving the noisy low-spatial frequency image signal to generate a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to said set of predetermined pixels-window characteristics;
a pixel-based local window segmentation processor comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels;
a counter generating a selected-pixels count signal;
a second pixel-based serial-to-parallel converter receiving the noisy input image signal to generate a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics;
a mean estimator combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel;
a minimum-mean-square-error filter receiving the noisy input image parallel signals, the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate said pre-filtered output image signal according to an input noise statistic signal.
- 10. An adaptive apparatus for spatially reducing noise in a noisy input image signal comprising:
an adaptive spatial noise reducer including;
a low-pass filter receiving said noisy input image signal to generate a noisy low-spatial frequency image signal;
a first pixel-based serial-to-parallel converter receiving the noisy low-spatial frequency image signal to generate a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to a set of predetermined pixels-window characteristics;
a pixel-based local window segmentation processor comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels;
a counter generating a selected-pixels count signal;
a second pixel-based serial-to-parallel converter receiving the noisy input image signal to generate a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics;
a mean estimator combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel;
a minimum-mean-square-error filter receiving the noisy input image parallel signals, the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a noise-filtered output image signal according to an input noise statistic signal; and
a controllable noise statistic estimator including;
a high-pass two-dimensional filter receiving the noisy input image signal to generate a noisy horizontal/vertical high-spatial frequency image signal;
a third pixel-based serial-to-parallel converter receiving the noisy horizontal/vertical high-spatial frequency image signal to generate a group of noisy horizontal/vertical high-spatial frequency image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics;
a statistic calculator combining the noisy horizontal/vertical high-spatial frequency image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a resulting noise statistic signal associated with the locally considered pixel;
a noise statistic estimator unit generating said input noise statistic signal from the resulting noise statistic signal. - View Dependent Claims (11, 12, 13, 14, 15, 16, 17)
wherein; g*(i,j;
x,y) is said values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with i=−
k, . . . ,+l;
j=−
m, . . . ,+n, said window having dimensions of (k+l)×
(m+n), k,l,m,n being appropriate positive integers;
g* (x,y) is said locally considered pixel value; and
T is a predetermined threshold value.
- 12. An apparatus according to
claim 11 , wherein a value of said mean pixel value signal is given by: - 13. An apparatus according to
claim 10 , wherein said input noise statistic signal is an input noise variance signal, said controllable noise statistic estimator is a controllable noise variance estimator, said statistic calculator is a variance calculator, said resulting noise statistic signal is a resulting noise variance signal, said noise statistic estimator unit is a noise variance estimator unit. - 14. An apparatus according to
claim 13 , wherein said noise variance estimator includes:a noise mean variance estimator receiving said input resulting noise variance signal to generate a mean variance signal;
a condition detector receiving said noisy input image signal to detect a spatial characteristic of said locally considered pixel;
a weighting device operatively connected to the condition detector for combining the mean variance signal with a weight value corresponding to a detected spatial characteristic to generate said input noise variance signal.
- 15. An apparatus according to
claim 10 , wherein said input noise statistic signal is an input noise standard deviation signal, said controllable noise statistic estimator is a controllable noise standard deviation estimator, said statistic calculator is a standard deviation calculator, said resulting noise statistic signal is a resulting noise standard deviation signal, said noise statistic estimator unit is a noise standard deviation estimator unit. - 16. An apparatus according to
claim 15 , wherein values of said segmented local window parallel signals are given by:wherein; g*(i,j;
x,y) is said values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with i=−
k, . . . ,+l;
j=−
m, . . . ,+n, said window having dimensions of (k+l)×
(m+n), k,l,m,n being appropriate positive integers;
g*(x,y) is said locally considered pixel value;
T is a predetermined threshold value;
wherein said mean pixel value signal is defined by;
- 17. An apparatus according to
claim 15 , wherein said noise standard deviation estimator includes:a noise mean standard deviation estimator receiving said input resulting noise standard deviation signal to generate a mean standard deviation signal;
a condition detector receiving said noisy input image signal to detect a spatial characteristic of said locally considered pixel;
a weighting device operatively connected to the condition detector for combining the mean standard deviation signal with a weight value corresponding to a detected spatial characteristic to generate said input noise standard deviation signal.
- 18. An adaptive method for spatially reducing noise in a noisy input image signal comprising the steps of:
i) filtering said noisy input image signal to generate a noisy low-spatial frequency image signal;
ii) converting the noisy low-spatial frequency image signal to a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to a set of predetermined pixels-window characteristics;
iii) comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels;
iv) generating a selected-pixels count signal;
v) converting the noisy input image signal to a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics;
vi) combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel; and
vii) processing the noisy input image parallel signals with a minimum-mean-square-error filter using the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a noise-filtered output image signal according to an input noise statistic signal. - View Dependent Claims (19, 20, 21, 22, 23)
wherein; g*(i,j;
x,y) is said values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with i=−
k, . . . ,+l;
j=−
m, . . . ,+n, said window having dimensions of (k+l)x (m+n), k,l,m,n being appropriate positive integers;
g*(x,y) is said locally considered pixel value; and
T is a predetermined threshold value.
- 20. A method according to
claim 19 , wherein a value of said mean pixel value signal is given by: - 21. A method according to
claim 18 , wherein said input noise statistic signal is an input noise variance signal, said minimum-mean-square-error algorithm including the steps of:a) generating an estimated variance signal from said noisy input image parallel signals, said segmented local window parallel signals, said selected-pixels count signal and said mean pixel value signal;
b) generating a weight signal from the estimated variance signal and the input noise variance signal;
c) subtracting said mean pixel value signal from the noisy input signal associated with the locally considered pixel to generate a difference signal;
d) multiplying the weight signal with the difference signal to generate a weighted difference signal;
e) adding the mean pixel value signal and the weighted difference signal to generate said noise-filtered output image signal.
- 22. A method according to
claim 18 , wherein said input noise statistic signal is an input noise standard deviation signal, said minimum-mean-square-error algorithm including the steps of:a) generating an estimated standard deviation signal from said noisy input image parallel signals, said segmented local window parallel signals, said selected-pixels count signal and said mean pixel value signal;
b) generating a weight signal from the estimated standard deviation signal and the input noise variance signal;
c) subtracting said mean pixel value signal from the noisy input signal associated with the locally considered pixel to generate a difference signal;
d) multiplying the weight signal with the difference signal to generate a weighted difference signal;
e) adding the mean pixel value signal and the weighted difference signal to generate said noise-filtered output image signal.
- 23. A method according to
claim 18 , wherein values of said segmented local window parallel signals are given by:wherein; g*(i,j;
x,y) is said values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with i=−
k, . . . ,+l;
j=−
m, . . . ,+n, said window having dimensions of (k+l)×
(m+n), k,l,m,n being appropriate positive integers;
g*(x,y) is said locally considered pixel value;
T is a predetermined threshold value;
wherein said mean pixel value signal is defined by;
- 24. An adaptive method for spatially reducing noise in a noisy input image signal comprising the steps of:
i) filtering said noisy input image signal to generate a noisy low-spatial frequency image signal;
ii) converting the noisy low-spatial frequency image signal to a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to a set of predetermined pixels-window characteristics;
iii) comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels;
iv) generating a selected-pixels count signal;
v) converting the noisy input image signal to a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics;
vi) combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel;
vii) processing the noisy input image parallel signals with a minimum-mean-square-error filter using the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a pre-filtered output image signal according to an input noise statistic signal;
viii) repeating said steps ii) to vii) according to at least one further complementary set of pixels-window characteristics to generate a further pre-filtered output image signal according to the input noise statistic signal; and
ix) averaging said output image signals to generate a noise-filtered output image signal. - View Dependent Claims (25)
wherein; g*(i,j;
x,y) is said values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with i=−
k, . . . ,+l;
j=−
m, . . . ,+n, said window having dimensions of (k+l)×
(m+n), k,l,m,n being appropriate positive integers characterizing said distinct set of predetermined pixels-window characteristics;
g*(x,y) is said locally considered pixel value;
T is a predetermined threshold value; and
wherein a value of said mean pixel value signal is given by;
- 26. An adaptive method for spatially reducing noise in a noisy input image signal comprising:
i) filtering said noisy input image signal to generate a noisy low-spatial frequency image signal;
ii) converting the noisy low-spatial frequency image signal to a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to a set of predetermined pixels-window characteristics;
iii) comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels;
iv) generating a selected-pixels count signal;
v) converting the noisy input image signal to a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics;
vi) combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel;
vii) filtering the noisy input image signal to generate a noisy horizontal/vertical high-spatial frequency image signal;
viii) converting the noisy horizontal/vertical high-spatial frequency image signal to a group of noisy horizontal/vertical high-spatial frequency image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics;
ix) combining the noisy horizontal/vertical high-spatial frequency image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a resulting noise statistic signal associated with the locally considered pixel;
x) generating an input noise statistic signal from the resulting noise statistic signal; and
xi) processing the noisy input image parallel signals with a minimum-mean-square-error filter using the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a noise-filtered output image signal according to the input noise statistic signal. - View Dependent Claims (27, 28)
wherein; g*(i,j;
x,y) is said values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with i=−
k, . . . ,+l;
j=−
m, . . . ,+n, said window having dimensions of (k+l)×
(m+n), k,l,m,n being appropriate positive integers;
g*(x,y) is said locally considered pixel value; and
T is a predetermined threshold value.
- 28. A method according to
claim 27 , wherein a value of said mean pixel value signal is given by:
1 Specification
The present invention relates to the field of image noise reducing, and more particularly to apparatus and method for adaptively reducing noise in a noisy input image or a sequence of images, wherein the noise is of an additive Gaussian type.
While existing prior patents and other publications on noise reducing techniques are abundant, they generally concern coring techniques applied in high frequency part of a considered image. Image de-noising techniques can be classified as spatial or temporal ones. Of course, a series combination of spatial and temporal techniques is possible and generally beneficial.
Temporal filters are generally applied for a sequence of images in which the noise component between two successive images is supposed to be non-correlated. The temporal filtering techniques are essentially based either on motion detection or motion estimation. The filter structure can be either infinite impulse response (IIR) or finite impulse response (FIR) with frame delay elements. In general, the temporal techniques are more performing than the spatial ones. The system cost is essentially due to the frame memory and the motion estimation hardware. Temporal de-noising techniques can be found, for example, in U.S. Pat. No. 5,161,018 to Matsugana; U.S. Pat. No. 5,191,419 to Wicshermann; U.S. Pat. No. 5,260,775 to Farouda; U.S. Pat. No. 5,404,179 to Hamasaki and in a recent publication entitled “The Digital Wetgate: A Third-Generation Noise Reducer”, Wischerman; G., SMPTE Journal, February 1996, pp. 95-100.
Spatial noise reducing techniques can be applied for either still pictures or sequence of images. These techniques are described in many available textbooks such as: “Fundamentals of Electronic Image Processing”, Weeks; A. R. Jr., SPIE Optical Engineering Press, Bellingham, Wash., 1996; “Two-Dimensional Signal and Image Processing”, Lim; J. S., Prentice-Hall, Englewood Cliffs, N.J., 1990; and “Nonlinear Digital Filters: Principles and Applications”, Pitas, J. and al., Kluwer Academic Publishers, Boston, 1990. In general, spatial noise reducing techniques can be divided further into three categories.
In the first category, the spatial nonlinear filters are based on local order statistics. Utilizing a local window around a considered pixel, these filters are working on this set of pixels ordered from their minimum to their maximum values. For example, the median filter, the min/max filter, the alpha-trimmed mean filter, and their respective variants can be classified in this category. These filters work well for removing impulse like salt-and-pepper noise. However, for the small amplitude noise these filters can blur some details or small edges.
In the second category the coring techniques are applied in another domain different from the original image spatial domain. The chosen domain partly depends on noise nature. U.S. Pat. No. 4,163,258 to Ebihara teaches the use of the Walsh-Hadamard transform domain, while U.S. Pat. No. 4,523,230 to Carlson et al. discloses some sub-band decomposition. Finally, the homomorphic filter, working in the logarithmic domain, is the classical one for removing multiplicative noise and shading from an image.
In the third category, the filters are locally adaptive and the noise removing capacity is varying from homogenous regions to edge regions. These filters give good results for additive Gaussian noise. A well-known filter in this category is the minimum-mean-square-error (MMSE) filter as originally published in “Digital image enhancement and noise filtering by use of local statistics”, Lee; J. S., IEEE Trans. on PAMI-2, March 1980, pp. 165-168. Referring FIG. 1, a general block diagram of the prior art Lee'"'"'s MMSE noise reducer is illustrated. Let a fixed dimension window centered on the considered or current pixel. The filtered pixel output f*(x,y) is additively composed of the local mean value obtained at output 12 of a mean estimator 10 and a weighted difference of the noisy pixel g(x,y) and the local mean intensity values. The optimum weight K determined by MMSE at an output 14, which corresponds to a kind of coring technique, is equal to the local variance ratio of the true image and the noisy one. The Lee'"'"'s MMSE filter efficiently removes noise in homogenous image regions while reserving the image edges. However, the noise essentially remains in edge or near-edge regions. Moreover, the required variance calculation is expensive for hardware implementation.
In “One-dimensional processing for adaptive image restoration”, Chan; P. et al., IEEE Trans. on ASSP-33, February 1985, pp. 117-126, there is presented a method for noise reducing in edge regions. The authors propose the use, in series, of four (4) one-dimensional MMSE filters respectively along 0°, 45°, 90° and 135° directions. The obtained results are impressive for large variance noise. However, for small noise, the filter can blur some image edges. Moreover, the noise variance output estimation at each filter stage require costly hardware.
For a same purpose, in “Digital image smoothing and the Sigma filter”, Lee; J. S., Computer Vision, Graphics, and Image Processing-24, 1983, pp. 255-269, there is proposed a said Sigma filter as illustrated in FIG. 2. For noise removing, this filter calculates with a segmentation processor 16 using a local window of 5×5 dimensions combined to a mean estimator 18, the mean value of similar pixel intensities to that of a central considered pixel g(x,y), to obtain f*(x,y). A pixel in the window is said similar to the considered pixel if the intensity difference between these two pixels is smaller than a given threshold value. Usually, the threshold value is set equal to twice the noise standard deviation. For small noise, the Sigma filter works well, except for some pixels with sharp spot noise. For the latter case, J. S. Lee has suggested also, in a heuristic manner, the use of immediate neighbor average at the expense of some eventually blurred picture edges. Generally, the MMSE filters yield better objective results in term of peak signal to noise ratio (PSNR), than the Sigma filter.
Independently to the Lee'"'"'s contribution, U.S. Pat. No. 4,573,070 to Cooper essentially discloses a Sigma filter for a 3×3 window. Moreover, in the same Patent, Cooper finally combines in a single configuration the said Sigma filter, an order statistic filter and a strong impulse noise reduction filter.
In the above-cited publication of A. R. Weeks Jr., there is described the adaptive double-window-modified-trimmed mean (DW-MTM) filter. This filter yields as output the trimmed mean value of similar pixel intensities to the median value in a local window. The adaptive DW-MTM filter is able to eliminate both salt-and-pepper noise and Gaussian noise but at the expense of filtered image blurring.
Finally, for correlated noise such as ringing/quantified noise in Discrete-Cosine-Transform (DCT) based decompressed image or cross-luminance noise in NTSC/PAL decoded image, there is still a need for a generic configuration or technique for implementing a real time noise reducer.
It is therefore an object of the present invention to provide an apparatus and a method for adaptively reducing noise in a noisy input image signal which provides spatial noise reduction in both homogenous and edge regions using a robust adaptive local segmented window, while preserving picture edges.
Another object of an aspect of the present invention is to provide, via a local segmented window, a technique that can adaptively estimate local mean and local standard deviation or variance in a non-stationary environment of a picture.
Yet another object of an aspect of the present invention is to provide an economic MMSE-based spatial noise reducer.
Yet another object of an aspect of the present invention is to provide a combined spatial noise reducer in which multiple local segmented windows are considered.
Yet another object of an aspect of the present invention is to provide a generic configuration for correlated noise variance estimator and for controllable mechanism of local noise reduction level.
According to one or more of the above objects, from a broad aspect of the present invention, there is provided an adaptive apparatus for spatially reducing noise in a noisy input image signal comprising a low-pass filter receiving the noisy input image signal to generate a noisy low-spatial frequency image signal. The apparatus further comprises a first pixel-based serial-to-parallel converter receiving the noisy low-spatial frequency image signal to generate a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to predetermined pixels-window characteristics, a pixel-based local window segmentation processor comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels, and a counter generating a selected-pixels count signal. The apparatus further comprises a second pixel-based serial-to-parallel converter receiving the noisy input image signal to generate a group of noisy input image parallel signals associated with the locally considered pixel and according to the predetermined pixels-window characteristics. The apparatus also comprises a mean estimator combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel, and a minimum-mean-square-error filter receiving the noisy input image parallel signals, the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a noise-filtered output image signal according to an input noise statistic signal.
According to a further broad aspect of the invention, there is provided a an adaptive apparatus for spatially reducing noise in a noisy input image signal comprising a plurality of parallel-connected adaptive spatial noise reducers each presenting a distinct set of predetermined pixels-window characteristics, each said noise reducer receiving the noisy input image signal to generate a corresponding pre-filtered output image signal and an averaging unit receiving each pre-filtered output image signal at a corresponding positive input thereof to generate a noise-filtered output image signal. Each noise reducer comprises a low-pass filter receiving the noisy input image signal to generate a noisy low-spatial frequency image signal and a first pixel-based serial-to-parallel converter receiving the noisy low-spatial frequency image signal to generate a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to the set of predetermined pixels-window characteristics. Each noise reducer further comprises a pixel-based local window segmentation processor comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels, a counter generating a selected-pixels count signal, and a second pixel-based serial-to-parallel converter receiving the noisy input image signal to generate a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics. Each noise reducer also comprises a mean estimator combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel, and a minimum-mean-square-error filter receiving the noisy input image parallel signals, the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate the pre-filtered output image signal according to an input noise statistic signal.
According to another broad aspect of the invention, there is provided an adaptive apparatus for spatially reducing noise in a noisy input image signal comprising an adaptive spatial noise reducer including a low-pass filter receiving the noisy input image signal to generate a noisy low-spatial frequency image signal, and a first pixel-based serial-to-parallel converter receiving the noisy low-spatial frequency image signal to generate a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to a set of predetermined pixels-window characteristics. The adaptive spatial noise reducer further includes a pixel-based local window segmentation processor comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels, a counter generating a selected-pixels count signal, and a second pixel-based serial-to-parallel converter receiving the noisy input image signal to generate a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics. The adaptive spatial noise reducer also includes a mean estimator combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel, and a minimum-mean-square-error filter receiving the noisy input image parallel signals, the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a noise-filtered output image signal according to an input noise statistic signal. The apparatus further comprises a controllable noise statistic estimator including a high-pass two-dimensional filter receiving the noisy input image signal to generate a noisy horizontal/vertical high-spatial frequency image signal, and a third pixel-based serial-to-parallel converter receiving the noisy horizontal/vertical high-spatial frequency image signal to generate a group of noisy horizontal/vertical high-spatial frequency image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics. The controllable noise statistic estimator further includes a statistic calculator combining the noisy horizontal/vertical high-spatial frequency image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a resulting noise statistic signal associated with the locally considered pixel, and a noise statistic estimator unit generating the input noise statistic signal from the resulting noise statistic signal.
According to another broad aspect of the invention, there is provided an adaptive method for spatially reducing noise in a noisy input image signal comprising the steps of: i) filtering the noisy input image signal to generate a noisy low-spatial frequency image signal; ii) converting the noisy low-spatial frequency image signal to a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to a set of predetermined pixels-window characteristics; iii) comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels; iv) generating a selected-pixels count signal; v) converting the noisy input image signal to a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics; vi) combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel; and vii) processing the noisy input image parallel signals with a minimum-mean-square-error filter using the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a noise-filtered output image signal according to an input noise statistic signal.
According to another broad aspect of the invention, there is provided an adaptive method for spatially reducing noise in a noisy input image signal comprising the steps of: i) filtering the noisy input image signal to generate a noisy low-spatial frequency image signal; ii) converting the noisy low-spatial frequency image signal to a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to a set of predetermined pixels-window characteristics; iii) comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels; iv) generating a selected-pixels count signal; v) converting the noisy input image signal to a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics; vi) combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel; vii) processing the noisy input image parallel signals with a minimum-mean-square-error filter using the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a pre-filtered output image signal according to an input noise statistic signal; viii) repeating said steps ii) to vii) according to at least one further complementary set of pixels-window characteristics to generate a further pre-filtered output image signal according to the input noise statistic signal; and ix) averaging said output image signals to generate a noise-filtered output image signal.
According to another broad aspect of the invention, there is provided an adaptive method for spatially reducing noise in a noisy input image signal comprising: i) filtering the noisy input image signal to generate a noisy low-spatial frequency image signal; ii) converting the noisy low-spatial frequency image signal to a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to a set predetermined pixels-window characteristics; iii) comparing values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with the locally considered pixel value to generate segmented local window parallel signals associated with selected pixels; iv) generating a selected-pixels count signal; v) converting the noisy input image signal to a group of noisy input image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics; vi) combining the noisy input image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a mean pixel value signal associated with the locally considered pixel; vii) filtering the noisy input image signal to generate a noisy horizontal/vertical high-spatial frequency image signal; viii) converting the noisy horizontal/vertical high-spatial frequency image signal to a group of noisy horizontal/vertical high-spatial frequency image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics; ix) combining the noisy horizontal/vertical high-spatial frequency image parallel signals with the segmented local window parallel signals and the selected-pixels count signal to generate a resulting noise statistic signal associated with the locally considered pixel; x) generating an input noise statistic signal from the resulting noise statistic signal; and xi) processing the noisy input image parallel signals with a minimum-mean-square-error filter using the segmented local window parallel signals, the selected-pixels count signal and the mean pixel value signal to generate a noise-filtered output image signal according to the input noise statistic signal.
The preferred embodiments of the present invention will now be described with reference to the accompanying drawings, in which:
FIG. 1 is a general block diagram of the prior art Lee'"'"'s MMSE noise reducer;
FIG. 2 is a general block diagram of the prior art Lee'"'"'s Sigma filter;
FIG. 3 is a block diagram illustrating a first preferred embodiment of a spatial noise reducer according to the invention, in which a single local adaptive window is utilized;
FIG. 4 illustrates some examples of window characteristics versus current pixel;
FIG. 5 is a block diagram representing a pixel-based sequential-to-parallel converter;
FIG. 6a is a graphic notation of a scalar product device;
FIG. 6b is a specific example of a scalar product device implementing two vectors each having three (3) components;
FIG. 7 illustrates a functional detailed diagram of the first basic embodiment of FIG. 1 using scalar product devices;
FIG. 8 is a general block diagram representing a second preferred embodiment of a spatial noise reducer according to the present invention, in the form of a parallel multiple local adaptive window spatial noise reducer;
FIG. 9 is a block diagram representing a third preferred embodiment of s spatial noise reducer according to the invention, which combines a single local adaptive window spatial noise reducer with a controllable noise variance estimator.
Referring now to FIG. 3 a first embodiment of the present invention in the form of a Single local Adaptive Window Spatial Noise Reducer (SAW-SNR) will be described. This basic configuration provides a segmentation-based MMSE spatial noise reducer for efficient edge preservation, and noise reduction in both homogenous or edges picture regions. The configuration can be also considered as a suitable combination of the Sigma and the MMSE filters. In the block diagram shown in FIG. 3, appropriate delays allowing signal synchronization required by the various operations of the signal processing are not represented, implementation of such delays being well known in the art. An example of such an implementation will be described later in reference to FIG. 7. A noisy image input signal g(x,y) at 301, which is typically a luminance component of the image, is applied to the SAW-SNR generally designated at 318. The couple (x,y) represents the current coordinates of a considered pixel. The input signal g(x,y) is a noisy version of an original image f(x,y). A noise component n(x,y) is supposed to be additive, zero mean and independent to the original image, i.e.:
The purpose of a noise reducer is to provide a cleaned output signal f*(x,y) which should be as close as possible to the original and unknown signal f(x,y). A usual criterion is the minimum of mean square error:
in which the image output signal f*(x,y) is a linear combination of the noisy observed image input g(x,y) and a local constant value:
Substituting equation (3) in equation (2) and minimizing the mean square error, the expression of f*(x,y) in equation (3) becomes:
wherein μ(x,y) is the local image mean value of f(x,y) or that of g(x,y)and the weight K(x,y) is determined by:
In the latter equation, the symbol σg(x,y) denotes the local standard deviation (SD) of the input image g(x,y), and σn corresponds in turn to that of the noise signal n(x,y). Substituting these results into equation (2) yields the minimum value of the performance indices:
In equation (6), the symbol σf2(x,y) denotes the local variance of the clean image f(x,y). The two equations 4 and 5 imply the knowledge of the local mean μ(x,y) and the local SD σg(x,y). Unfortunately, the latter values are generally unknown. They should be estimated by some suitable manner.
In the Lee'"'"'s MMSE filter, these values are calculated in a local fixed dimension window around the central pixel (x,y) as follows
and
wherein the summation are done over all i and j in the given local window, g(i,j;x,y) being values of a group of noisy input image parallel signals, NT denoting the pixel total number in the window. In homogenous regions where pixels are similar, equation (7) yields generally good results. However, in near edge or edge regions, the estimation given by equation (7), which does not take care of the non-stationary signal, becomes erroneous. In a similar manner, equation (8) can overestimate the signal variance in edge or near edge regions, the weight K(x,y) in equation (5) becomes unity (1) and the system given by equation (4) does not reduce any more noise in such cases. Moreover, the square values in equation (8) imply many expensive multipliers in an hardware implementation.
Meanwhile, in the Lee'"'"'s Sigma filter, the weight K(x,y) is set to be fixed and equal to zero (0), which imply that the Sigma filter is always sub-optimum, and the mean value estimation is calculated only with selected similar intensity pixels in the window, that is:
wherein the summation are done again over all i and j in the given local window, N bein a value of said selected-pixels count signal. However, in order to consider both homogenous and non-homogenous regions, the pixel intensities g(i,j;x,y) are weighted by window binary signals w(i,j;x, y) defined as follows:
wherein g(x,y)=g(0,0;x,y) is the intensity of the window central and considered pixel, T being a predetermined threshold value. The suggested threshold value in Sigma filters is equal to 2σn. The higher the threshold value T, the more noise is reduced, but at the expense of blurring true small image edges. If the threshold value T is small, equation (11) will leave unchanged isolated strong spot noise.
In order to combine the main advantages of Lee'"'"'s MMSE and Sigma filters, the noisy input signal g(x,y) shown at 301 in FIG. 3, is firstly fed to a preliminary low pass (LP) filter 302. The low pass filter 302 plays a very important role by partly reducing strong spot noise. Its output g*(x,y) allows robust segmentation in the presence of noise. For additive noise, it can be a simple filter with impulse response function lp(x,y) described by the following matrices:
The filter output g*(x,y) 303 is sent to a first pixel-based serial-to-parallel converter (PB-SPC) 306 receiving the noisy low-spatial frequency image signal to generate a group of noisy low-spatial frequency image parallel signals associated with a locally considered pixel and according to predetermined pixels-window characteristics, which parallel signals generated at 303′ are fed to a pixel-based local window segmentation processor 304 that selects pixels in a local window which are similar to the considered current pixel. An implementation example of a PB-SPC will be described later with reference to FIG. 5. The segmentation working with filtered image, yields at its output 320 a set of window binary signals w(i,j;x,y) defined now as:
wherein:
g*(i,j;x,y) is the values of the noisy low-spatial frequency image parallel signals associated with pixels included within the window with i=−k, . . . ,+l; j=−m, . . . ,+n, the window having dimensions of (k+l)×(m+n), k,l,m,n being appropriate positive integers;
g*(x,y) is the locally considered pixel value; and
T is a predetermined threshold value.
Since the noise component in the filtered image g*(x,y) is reduced, the threshold value T fed at an input 305 can be also reduced. Practically, the threshold value T can be set to σn or usually to a fixed value selected from 8 to 16 in an eight (8) bits video system, if noise is weak enough. The original noisy image signal g (x,y) at 301 is sent to a second PB-SPC 321 receiving the noisy input image signal to generate a group of noisy input image parallel signals g(i,j;x,y) associated with the locally considered pixel and according to the predetermined pixels-window characteristics. The parallel signals g(i,j;x,y) generated at the PBSPC output 322 connected to a line 323 and the window binary signals w(i,j;x,y) at 320′ are sent together to the mean value μ estimator 307 having an output 308 described by the above equations (9a) and (9b). In these equations, the term N is simply a count of selected pixels in the local window which are similar to the considered pixel (x,y) within the tolerance of threshold value T. A selected-pixels counter 327 with signal input w(i,j;x,y) 320 can be done by appropriate adders to realize equation (9b) above. The counter output signal at 328 represented in a few bits is sent to the mean value μ estimator 307. The SAW-SNR further comprises a minimum-mean-square-error filter generally designated at 329 receiving the noisy input image parallel signals g(i,j;x,y) at 322, the segmented local window parallel signals w(i,j;x,y) at 330, the selected-pixels count signal N at 331 and the mean pixel value signal μ(x,y) at 332 to generate a noise-filtered output image signal f*(x,y) according to an input noise statistic signal at 314. More specifically, the filter 329 comprises a statistical calculator 312 receiving the above signals g(i,j;x,y), w(i,j;x,y) and N to either generate a standard deviation (SD) signal or a variance signal at 313. In order to reduce the inherent complexity of equation (8) above, the proposed calculator 312 is preferably a standard deviation calculator yielding a SD value σg(x,y) at 313 estimated by absolute deviation mean as follows:
In the latter equation, N is again the number of selected similar pixels in the local window. Moreover, depending on the noise distribution, an appropriate value for constant C can be chosen equal to about 1.25 for additive Gaussian noise, or to about 1.15 for additive uniform noise. The filter 329 further comprises a weight calculator 319. The current SD σg(x,y) at 313 as given by equation (14) and the noise input SD σn at 314 are preferably applied to the weight calculator 319 generating weight K(x,y) as already described by equation (5) above. Although less efficient, the statistical calculator 312 may also be a variance calculator implementing equation (8) to generate a variance value σg2(x,y) in which case the noise input variance σn2 is fed to the weight calculator 319 to generate weight K(x,y) using same equation (5). In the case of white noise of constant variance σg2, instead of a weight calculator 319, a lookup table (LUT) is preferably used, which is built by pre-calculating necessary weights K(x,y). In the latter case, the value of constant C in equation (14) can be chosen fixed and the input noise variance σn2 represented by a few bits can be entered by an end-user to obtain the desired noise reduction level.
After proper delaying, the noisy image signal at 322 is applied in turn to an adder 309 included in the filter 329 through an input line 324. The adder output signal at 325 representing the current difference (g(x,y)−μ(x,y)) is weighted by the current weight K(x,y) signal at 315 by a multiplier 310 included in the filter 329. The multiplier output 316 and the mean value μ(x,y) are applied together to an adder 311 included in the filter 329 to yield the filtered image signal f*(x,y) 316 of the desired Single local Adaptive Window Spatial Noise Reducer (SAW-SNR) 318.
It is worthwhile to note that, in addition to the usual noisy input g(x,y) at 301 and the filtered output f*(x,y) at 316, the SAW-SNR 318 has also an input 300 for a window characteristics data (WIC) signal. This signal is used to select some window size and its relative position to the considered current pixel at the coordinates (x,y). In another words, the selected window is usually, but not absolutely necessary centered to the current pixel. Some pixel offsets or different window sizes can be useful for a more sophisticated spatial noise reducer, which will be described later according to a second preferred embodiment of the present invention.
Referring now to FIG. 4, there is illustrated a few examples of window characteristics relatively to the current pixel position. The window 401 is a 5×5 square one centered to the current pixel 400, a centered square window being suitable for progressive image. The dotted window 402 is also a centered window of dimension 3×7, i.e., 3 lines and 7 columns, an horizontal rectangular window being suitable for interlaced signal. For the case of a single adaptive window spatial noise reducer, centered windows are recommended. Moreover, equations (7), (9) and (10) above show that, at least in a homogeneous region, the larger the window the better the result. FIG. 4 illustrates also two windows 403 and 404 that are not centered on the current pixel. The window 403 has one pixel being offset to the left, while the window 404 has one pixel being offset to the right. Horizontal offset windows are suitable for economic purpose in a Parallel Multiple local Adaptive Window Spatial Noise Reducer (Parallel M-AW-SNR), which will be described later in detail with reference to FIG. 8.
Before referring to FIG. 7, which illustrates a functional implementation block diagram of the first embodiment in detail, in order to simplify and facilitate the presentation, let'"'"'s consider now FIG. 6 that represents an example of a pixel-based sequential-to-parallel converter (PB-SPC) of 3×5 dimension. In order to create a local window for filtering or for segmentation, it is necessary to use a PB-SPC that converts a scanned image signal to multiple parallel signals corresponding to all pixel intensity values in the given window. The scanned image input signal g(x,y) at 501 is simply applied to appropriate line and pixel delays in series. For a window of 3×5 dimension, it takes at least two line delays 502, 504 and three (3) series of four (4) pixel delays 510-516, 520-526 and 530-536, respectively. In this example, if the window is centered, the current pixel is located at an output 523. The other output signals 501, 503, 505, 511, . . . ,513, 533, . . . ,537, correspond to the current pixel neighbors in the window. Turning now to FIG. 6a, there is illustrated a graphic notation of a scalar product (SP) from two vectors {right arrow over (a)}=(a1, a2, . . . ,an) at an input 611 and {right arrow over (b)}=(b1, b2, . . . ,bn) at an input 612, according to the following relation:
wherein the summing is done over i for i=1,2, . . . ,n. In the notation, n is the component number of the vectors {right arrow over (a)} and {right arrow over (b)}. The scalar product (SP) at 616 is done by a scalar multiplier 613. Some implementation details of the scalar multiplier are shown in FIG. 6b. It consists of n (illustrated with n=3) parallel multipliers 613-1, 613-2, . . . , 613-n and an adder 615 of n multiplier results 614-1, 614-2, . . . ,614-n. The inputs of the ith multiplier 613-i are the components ai, 611-i and bi, 611-i of the respective vectors {right arrow over (a)} and {right arrow over (b)}.
Referring now to FIG. 7, an implementation for the first embodiment of the proposed SAW-SNR basically described above with reference to in FIG. 3 will be now explained in detail. The noisy input image g(x,y) at 701 is applied to an input PB-SPC 702 of dimension 3×3 as part of the low-pass filter represented at 302 in FIG. 3. The PB-SPC vector output at 717 is applied in turn to a SP 741 that also receives a vector of low pass filter coefficients at 740. Assuming that equation (12a) above describes the low pass filter, all component values of the coefficient vector at 740 are equal to unity (1). The proportional term {fraction (1/9)} in equation (12a) can be further absorbed by the threshold value in equation (13) above. Meanwhile, the output 718 of the PB-SPC 702, corresponding to the current pixel, is an appropriate delayed version of the input noisy image signal g(x,y). The SP output 703, corresponding to the low pass output signal 303 in FIG. 3, is sent to a PB-SPC 704, corresponding to the PB-SPC 306 in FIG. 3, of desired window dimension in order to get the pixel intensity vector at 705 of the given local window. A single input 706 represents the current filtered pixel intensity output signal g*(x,y). Let NT being the window total pixel number and also the component number of the vector at 705 in which the signal component at 706 corresponds to the current considered pixel intensity. The vector at 705 and the current component at 706 are then applied to a vector adder 707 as illustrated by a double circle, which in fact includes NT parallel adders, to calculate the NT differences (g*(i,j;x,y)−g*(x,y)) as required by equation (13) above. The vector adder output 708 is applied to a vector full-wave rectifier 739 necessary now for calculating |g*(i,j;x,y)−g(x,y)| for all i and j in the window. For segmentation purpose, the rectifier output vector at 740 and the threshold vector at 709 are sent to a vector adder 710, from which a resulting vector at output 711 is fed to a comparison detector 712. The detector output vector {right arrow over (w)} at 713 composed of NT binary components {w(i,j;x,y)} is sent now to three (3) SP 715, 724 and 728. The SP 715 is simply provided as part of the selected similar pixel counter, and is provided with a second input 714 fed by a constant and unitary (1) vector. The SP output 716 represents the local selected similar pixel number N, corresponding to the current pixel intensity. Meanwhile, the delayed version of the input noisy signal g (x,y) at 718 is sent to a PB-SPC 719, corresponding to the PB-SPC 321 in FIG. 3, of the same desired window dimension than that of PB-SPC 704. The parallel output vector at 720 and the segmented window vector {right arrow over (w)} at 713 are applied together to the SP 724 that yields at its output 735 the summing value of similar pixel intensities in the local window. The total and local value signal at 735 and the local number N signal at 716 are applied to a ROM or LUT 737 to get a local mean μ at 721. It is pointed out that a few bits can represent the value N. The ROM or LUT 737 is proposed to avoid an expensive division. In order to determine the local SD σg of the noisy input image, the vector of NT parallel windowed noisy pixels at 720 and the obtained local mean μ at 721′ are applied also to a vector adder 722. The vector adder output 723 is sent in turn to a vector full-wave rectifier 725 generating at 726 an output difference vector {|g(i,j;x,y)−μ(x,y)|} of NT components. The difference vector signals and the segmented window vector {right arrow over (w)} signals are sent to the SP 728 that yields at the output 729 the weighted sum (ΣΣw(i,j;x,y)·(|g(i,j;x,y)−μ(x,y)|)) as necessary for equation (14) above. The multiplication by a constant C equal to about 1.25 for Gaussian noise and the division by N selected similar pixel number at 716 can be integrated in a ROM or LUT 731 generating at an output 743 an estimated local SD σg of the noisy input image, which is fed with the SD σn of the input noise at 730 to a weight calculator 742 implementing equation (5) above. The current calculated weight output K(x,y) at 732 is used to weigh the current difference (g(x,y)−μ(x,y)) at 727 via a multiplier 733. The multiplier weighted output signal at 736 is sent to an adder 734. The above weighted current difference signal at 736 and the local mean μ(x,y) at 721, via an adder 734, give the noise-filtered output image signal f*(x,y) at 738.
Referring now to FIG. 8, a second embodiment of the present invention in the form of a Parallel Multiple local Adaptive Window Spatial Noise Reducer (M-AW-SNR) will be now described. Let g(x,y) be the noisy image input signal at 801, which corresponds to the same input signal at 301 in FIG. 3. The noisy signal is applied simultaneously to L parallel SAW-SNRs, respectively designated at 318-1, 318-2, . . . and 318-L. Each SAW-SNR 318, which was basically described above with reference to FIG. 3, is working now with a distinct set of window characteristics WIC fed at inputs 300-1 to 300-L. The SAW-SNR outputs, 316-1 to 316-L, are added together via an adder 820 as part of an averaging unit 824. The adder output 821 is divided by L and then rounded by a calculator 822 included in the averaging unit 824. The result output 823 is also the desired final signal of the Parallel M-AW-SNR. It is to be understood that the L low pass filters included in the L parallel SAW-SNRs can be replaced by a single low pass filter 802 indicated in dotted lines in FIG. 8. The proposed Parallel M-AW-SNR configuration can be justified by its performance. Let hi(x,y) be the ith SAW-SNR output, 318-i, defined as follows:
The residual noisy component ri(x,y) is supposed to be non-correlated to the clean image f(x,y). This supposition is valid at least in edge regions where noise can remain after applying a SAW-SNR. The correlation between two components ri(x,y) and rj(x,y) resulted from two different but overlapped windows with similar size, is supposed equal to followings:
The term σr2 denotes the variance of the residual SAW-SNR noise, and ρij is the normalized cross-correlation between noisy components ri(x,y) and rj(x,y). Let JPM and J denoting respectively the performance indices of a Parallel M-AW-SNR and a SAW-SNR. It can be shown that the improvement J/JPM is equal to:
J/JPM=L2/ΣΣρij>1 (18)
The improvement has been verified experimentally even in the case of a SAW-SNR, which utilizes, as window, the union of all individual ones in a given Parallel M-AW-SNR. The proposed Parallel M-AW-SNR is an averaging configuration that can also be justified by its simplicity. In fact, it would be much more complicated to obtain an optimum linear combination with L SAW-SNRs output images. The Parallel M-AW-SNR configuration is composed of L parallel SAW-SNRs as previously described. It could be also possible to utilize SAW-SNRs in series (not shown). If the residual noise after each SAW-SNR is supposed to be still non-correlated to the clean image f(x,y), it can be theoretically shown that the performance indicia JSM associated with a Serial M-AW-SNR is given by the following equation:
In equation (19), the terms σf2(x,y) and σn2 denote respectively the local clean image variance and the original input noise variance. Hence, for the Serial M-AW-SNR, a theoretical improvement of:
over the single SAW-SNR is possible. However, this expression shows that, in edge regions, the improvement is not necessary appreciable. Moreover, serial configuration requires estimation of each individual residual noise, which imply a costly hardware implementation. Therefore, the Parallel M-AW-SNR is recommended over the Serial M-AW-SNR.
Turning now to FIG. 9, a third preferred embodiment of the present invention in the form of a Single local Adaptive Window Spatial Noise Reducer with controllable noise statistic estimator will be now described. The noisy image signal g(x,y) is applied at 301, as usual, to the basic SAW-SNR 318 generating the noise-filtered output image signal f*(x,y) at 316. At the same time, the noisy image is sent also to a High-Pass Two-Dimensional Filter HP2DF 901, which extracts only very high vertical and horizontal frequency components in the input image. The filter band-pass bandwidth is chosen such that the filter output 902 contains mainly noise components, but not clean image. Justification for the choice of a HP2DF can be explained generally by the wide-band nature of noise. The utilize HP2DF is preferably a diamond filter, an example of impulse response function thereof being as follows:
The noisy horizontal/vertical high-spatial frequency image signal at 902 is fed to a PBSPC 911 to generate at output 913 a group of noisy horizontal/vertical high-spatial frequency image parallel signals associated with the locally considered pixel and according to the set of predetermined pixels-window characteristics. The BPSPC output 913 with the segmented local window parallel signals w(i,j;x,y) at 320 and the selected-pixels count signal Nat 328 available from the SAW-SNR 318 are sent to a statistic calculator 912 to generate a resulting noise statistic signal associated with the locally considered pixel. As explained before regarding the first preferred embodiment, while the statistic calculator 912 is preferably a SD calculator generating a SD resulting signal σR(x,y) at 904, it could be a variance calculator directly generating a variance resulting signal σR2(x,y) at 904. The calculator 912 is similar to the calculator 312 described before with reference to FIG. 3, but provided with a built-in mean value μ estimator. The calculator resulting signal at 904 is applied at the input of a look-up table LUT 905 that estimates, in turn, at its output 907, a mean value of local noise input SD σm(x,y) (or variance σm2(x,y) ). The LUT input-output relationship between the two local standard deviations σr(x,y) (or variance σr2(x,y)) and σm(x,y) (or variance σm2(x,y)) can be described by the following method. Let consider the linear portion in Eq.5 expressed by:
wherein the unknown additive noise variance σn2(x,y) is supposed now to be varying. It is thus necessary to pre-estimate this variance value for each pixel located at (x,y).
In many situations where the processing is well defined, such as for NTSC or PAL encoding/decoding and DCT based compression/decompression, an available original and clean test signal f(x,y) can be used for noise pre-learning. FIG. 9 illustrates partly a proposed configuration at 960 used for performing an off-line noise variance pre-estimation. The original test signal f(x,y) at 950 is applied to the above-mentioned processing at 951 which gives a test noisy image signal g(x,y) at 952. The additive test noise signal n(x,y) at 954 is then obtained by the difference (g(x,y)−f(x,y)) provided by an adder 953 and is sent in turn to a statistic calculator 955 similar to the calculator 912. The test noise SD σm(x,y) (or the test noise variance σm2(x,y)) estimation is done in the same context of that of the signal g(x,y) in the proposed SAW-SNR 318, with the segmented window parallel signals w(i,j;x,y) at 320 and the selected-pixels count signal N at 328. In words, for a considered pixel at (x,y), one obtains a pair of SD values (σr(x,y),σn(x,y)) (or a pair of variance values (σr2(x,y),σn2(x,y))). For the whole test picture or set of test pictures, a given value of σr(orσr2) can have many resulted values of σm(or σm2). In order to obtain a unique input-output relationship for the LUT 905, it is necessary, for a given σr(or σr2), to define a single value σm representing all possible values of σm. For the preferred SD calculation, proposed estimations for σm are as follows:
or
The estimation (22) or (23) can be done then on an off-line basis by a data storage and estimation device 957. The input-output result (σr,σm) at 904 and 958 respectively, permits the establishment of a pre-calculated LUT 905 for real time processing involving an unknown image. If the memory LUT 905 is large enough, some controllable bits can be fed at input 903 representing a learning or functional condition, for example for NTSC, PAL or 12 Mbit MPEG. The main requirement of the method is the prior knowledge of the processing to create the noisy image g(x,y) from the clean image f(x,y). The LUT output signal σm(x,y) at 907 could be sent directly (not shown) to the input 314 of the SAW-SNR 318 as local SD of the input noise. However, the estimated value σm(x,y) in equations 22 or 23 is only a representative one and is not necessary always good for all spatial varying image conditions. In such cases, the controllable noise statistic estimator preferably comprises a conditional weighting unit 909 which could be another LUT that receives the signal σm(x,y) at 907 and which is controlled, in turn, by the output 910 of a condition detector 908. The condition detector is driven by the input noisy image signal g(x,y) through line 301′. Depending on the application nature, The condition detector 908 classifies the current pixel as belonging to some specific regions such as edge regions, flat regions, near edge near flat regions, texture regions etc., and accordingly applies an appropriate weight or function at 910 to the mean value σm(x,y) fed at 907 for each region. As an example, for the case of ringing noise, the conditional weighting output can be (½)·σm(x,y) in contour or texture regions, 2·σm(x,y) in near edge near flat regions, (¼)·σm(x,y) in flat regions (to protect small details that does not necessarily create the ringing noise) and σm(x,y) elsewhere. The conditional weighting output is then sent to the SAW-SNR input 314 as local SD (or variance) of the input noise. It is to be understood that, for specific noise such as the coding block effect, the proposed method may be used provided appropriate specific conditional techniques and data are applied.