← Back to Patent Explorer
US5909515A Expired / Lapsed US ⏱️ 0.0 years remaining

Method for the temporal filtering of the noise in an image of a sequence of digital images, and device for carrying out the method

Assignee: US Philips Corp
•
Filing: 1996-07-01
•
Granted: 1999-06-01
•
Expires: 2016-07-01
📄 Original PDF 🌐 Google Patents 🔄 Sync Now

Abstract

Abstract A method for the temporal filtering of an image of a sequence of noisy digitized images includes the evaluation of an integration equation which produces the filtered intensity (P t C ) of a given pixel A t (x,y)! in the present image (J t P ) by evaluating the sum of the noisy intensity (I t P ) of the present (t) and the product of a gain factor (K t P ) and the difference of said noisy intensity of the present (I t P ) and the causal filtered intensity (P t-1 C ). According to this evaluation, the gain factor (K t C ) is a function of the causal gain factor (K t-1 C ) and of a continuity coefficient (α t C ) which measures the probability of intensity continuity between the present noisy intensity (I t P ) and the causal filtered intensity (P t-1 C ). A device for carrying out this method includes a sub-assembly whose input (100) receives the present noisy sample (I t P ) and which includes a memory (10) which supplies a causal filtered intensity (P t-1 C ), a memory (17) which supplies a causal gain factor (K t-1 C ), and means for calculating the integration equation.

Inventors

  • 👤 Sherif Makram-Ebeid

Bibliographic Details

Application Number —
Jurisdiction US
Data Source google_patents
Last Synchronized 2026-11-09 11:21
CPC Classifications
G; G06; G06T; G06T5/00; G06T5/70; G06V; G06V10/00; G06V10/20; G06V10/30; G06T2207/00; G06T2207/10; G06T2207/10016; G06T2207/10116; G06T2207/10121; G06T2207/20; G06T2207/20172; G06T2207/20182; G06T2207/30; G06T2207/30004

Claims Analysis

Total Claims: 21
Claim 1
Independent Claim 171 words
A method for temporal filtering of noise in an image (J t P ), which is referred to as the present image and forms part of a sequence of images in the form of two-dimensional matrices of pixels which have digitized, noisy intensity values which are referred to as samples (I j (x,y)), the method comprising: evaluating a present filtered causal sample (P t C ), sample which is the temporal filtering of a present noisy sample (I t P ) corresponding to a pixel in a given location (x,y) in the present image, by forming a linear combination of the present noisy sample (I t P ) and earlier samples in the sequence, referred to as noisy causal samples (I j C ), multiplied by weights (b j C ), which constitute coefficients between 0 and 1 of a probability of intensity continuity between said noisy causal sample, whereto the weight is attached, and the present noisy sample, the weight attached to the present noisy sample having the value 1.
Claim 2
Independent Claim 171 words
A method for temporal filtering of noise in an image (J t P ), which is referred to as the present image and forms part of a sequence of images in the form of two-dimensional matrices of pixels which have digitized, noisy intensity values which are referred to as samples (I j (x,y)), the method comprising: evaluating a present filtered causal sample (P t C ), sample which is the temporal filtering of a present noisy sample (I t P ) corresponding to a pixel in a given location (x,y) in the present image, by forming a linear combination of the present noisy sample (I t P ) and earlier samples in the sequence, referred to as noisy causal samples (I j C ), multiplied by weights (b j C ), which constitute coefficients between 0 and 1 of a probability of intensity continuity between said noisy causal sample, whereto the weight is attached, and the present noisy sample, the weight attached to the present noisy sample having the value 1.
Claim 3
Dependent Claim 71 words
A method as claimed in claim 1, in which the weight attached to a given causal noisy sample is the product of coefficients of probability of continuity between the successive noisy causal samples, from the given causal noisy sample until the present noisy sample, and in which the linear combination of the causal noisy samples and present noisy sample is normalized by the sum of the weights relating to said samples.
Claim 4
Dependent Claim 113 words
A method as claimed in claim 2, in which for the evaluation of the present filtered causal sample (P t C ) of the present instant (t) the linear combination is formed by performing the evaluation of a recursive relation, referred to as a causal integration relation, comprising the determination of a sum of the filtered causal sample (P t C ) of the preceding instant (t-1) and the product of a factor, referred to as a causal gain factor (K t C ) of the present instant, and a difference between the present noisy sample (I t P ) and the filtered causal sample (P t-1 C ) of the preceding instant.
Claim 5
Dependent Claim 70 words
A method as claimed in claim 3, in which said causal gain factor (K t C ) of the present instant (t) is evaluated in a recursive manner by the ratio of a causal gain factor (K t-1 C ) of the preceding instant (t-1) to the sum of said causal gain factor (K t-1 C ) of the preceding instant and a causal continuity coefficient (α t C ).
Claim 6
Dependent Claim 83 words
A method as claimed in claim 4, in which the evaluation of the causal integration relation comprises: the determination of the absolute value of a causal difference (Δ C ) between the present noisy sample (I t P ) of the present instant (t) and the filtered causal sample (P t-1 C ) of the preceding instant, and the determination of the causal continuity coefficient (α t C ) as a decreasing function (F C ) of said causal difference (Δ C ).
Claim 7
Dependent Claim 70 words
A method as claimed in claim 5, in which for the evaluation of the causal continuity coefficient (α t C ) the decreasing function is a causal function (F C ) whose argument (z C ) is given by the causal difference (Δ C ), normalized by a factor (S t C ) that takes into account the variance of the noise relative to each sample of the causal difference.
Claim 8
Dependent Claim 57 words
A method as claimed in claim 6, in which the causal decreasing function (F C ) has a constant maximum value smaller than or equal to 1 when its argument (Z C ) is between 0 and 1 and decreases to a minimum value equal to or larger than zero when its argument is larger than 1.
Claim 9
Dependent Claim 53 words
A method as claimed in claim 7, in which the filtering power is controlled by the maximum and minimum values of the causal decreasing function (F C ) which are chosen so as to smooth the noise peaks and to smooth the noise teeth which follow the discontinuities due to a spatial movement.
Claim 10
Independent Claim 48 words
The method of claims 6 wherein said normalization factor (S t C ) is evaluated by multiplying a standard deviation of the noise (σ B ) by the square root of two divided by two minus said causal gain factor (K t-1 C ) of the preceding instant.
Claim 11
Dependent Claim 63 words
A method as claimed in claim 4, in which said causal gain factor (K t C ) is the inverse of the sum of the weights attached to the causal samples of the linear combination, and the causal continuity coefficient (α t C ) is the weight of the causal sample (I t-1 C ) preceding the present sample (I t P ).
Claim 12
Dependent Claim 83 words
A method as claimed in claim 10, in which the evaluation of the causal integration relation comprises: the determination of the absolute value of a causal difference (Δ C ) between the present noisy sample (I t P ) of the present instant (t) and the filtered causal sample (P t-1 C ) of the preceding instant, and the determination of the causal continuity coefficient (α t C ) as a decreasing function (F C ) of said causal difference (Δ C ).
Claim 13
Dependent Claim 70 words
A method as claimed in claim 11, in which for the evaluation of the causal continuity coefficient (α t C ) the decreasing function is a causal function (F C ) whose argument (z C ) is given by the causal difference (Δ C ), normalized by a factor (S t C ) that takes into account the variance of the noise relative to each sample of the causal difference.
Claim 14
Dependent Claim 57 words
A method as claimed in claim 12, in which the causal decreasing function (F C ) has a constant maximum value smaller than or equal to 1 when its argument (z C ) is between 0 and 1 and decreases to a minimum value equal to or larger than zero when its argument is larger than 1.
Claim 15
Dependent Claim 53 words
A method as claimed in claim 13, in which the filtering power is controlled by the maximum and minimum values of the causal decreasing function (F C ) which are chosen so as to smooth the noise peaks and to smooth the noise teeth which follow the discontinuities due to a spatial movement.
Claim 16
Independent Claim 48 words
The method of claims 12 wherein said normalization factor (S t C ) is evaluated by multiplying a standard deviation of the noise (α B ) by the square root of two divided by two minus said causal gain factor (K t-1 C ) of the preceding instant.
Claim 17
Independent Claim 271 words
A device for temporal filtering of noise in images which are part of a sequence of images, each image being in the form of a matrix of pixels, comprising: an image processing system for supplying a present noisy digitized sample (I t P ) of a pixel having a given location (x,y) in an image (J t P ) in the form of a matrix of pixels arriving at a present instant (t), and a sub-assembly, referred to as a causal sub-assembly, whose input receives the present noisy sample (I t P ) and which comprises a memory for supplying a digitized filtered causal sample (P t-1 C ) of the preceding pixel, in the same location (x,y) in the matrix of the instant at the given location and a digitized noisy causal sample (I t-1 C ) of the preceding instant at the given location, a memory for supplying a causal gain factor (K t-1 C ) of the preceding instant, and calculation means for evaluating a causal integration equation and for delivering on its output the value of a present filtered causal sample (P t C ), which is the temporal filtering of the noisy present sample (I t P ), wherein the causal integration equation is a recursive relation comprising a sum of the filtered causal sample (P t-1 C ) of the preceding instant (t-1) and the product of a causal gain factor (K t C ) of the present instant (t), and a difference between the present noisy sample (I t P ) and the filtered causal sample (P t-1 C ) of the preceding instant.
Claim 18
Independent Claim 271 words
A device for temporal filtering of noise in images which are part of a sequence of images, each image being in the form of a matrix of pixels, comprising: an image processing system for supplying a present noisy digitized sample (I t P ) of a pixel having a given location (x,y) in an image (J t P ) in the form of a matrix of pixels arriving at a present instant (t), and a sub-assembly, referred to as a causal sub-assembly, whose input receives the present noisy sample (I t P ) and which comprises a memory for supplying a digitized filtered causal sample (P t-1 C ) of the preceding pixel, in the same location (x,y) in the matrix of the instant at the given location and a digitized noisy causal sample (I t-1 C ) of the preceding instant at the given location, a memory for supplying a causal gain factor (K t-1 C ) of the preceding instant, and calculation means for evaluating a causal integration equation and for delivering on its output the value of a present filtered causal sample (P t C ), which is the temporal filtering of the noisy present sample (I t P ), wherein the causal integration equation is a recursive relation comprising a sum of the filtered causal sample (P t-1 C ) of the preceding instant (t-1) and the product of a causal gain factor (K t C ) of the present instant (t), and a difference between the present noisy sample (I t P ) and the filtered causal sample (P t-1 C ) of the preceding instant.
Claim 19
Dependent Claim 260 words
A device as claimed in claim 16, in which the calculation means for evaluating the causal integration equation and for delivering the filtered causal sample (P t C ) comprise: an adder for determining an absolute value of a causal difference (Δ C ) between the present noisy sample (I t P ) of the present instant (t) and the stored filtered causal sample (P t-1 C ) of the preceding instant (t-1), and tables for determining a causal continuity cofficient (α t C ) as a decreasing function (F C ) of the causal difference (Δ C ) normalized by a factor (S t C ) of the variance of the noise (σ B ), a table for determining a factor (K t C ), referred to as the causal gain factor of the present instant (t) and equal to the quotient of the stored causal gain factor (K t-1 C ) of the preceding instant (t-1) and the sum of said causal gain factor (K t-1 C ) of the preceding instant and the causal continuity coefficient (α t C ), adders and a multiplier for determining the sum of the filtered causal sample (P t-1 C ) of the preceding instant (t-1) and the product of said causal gain factor (K t C ) of the present instant and a difference between the noisy present sample (I t P ) and the filtered causal sample (P t-1 C ) of the preceding instant, thus producing the value of the present filtered causal sample (P t C ).
Claim 20
Dependent Claim 72 words
A device for observation of a medical image, comprising; a system for supplying data of a sequence of medical images, each image being in the form of a two-dimensional matrix of pixels having digitized intensity values, an image display system, and an image processing system having access to the data of the sequence of images and to the display system, said image processing system including a device as claimed in claim 10.
Claim 21
Dependent Claim 73 words
A device for observation of a medical image, comprising; a system for supplying data of a sequence of medical images of medical images, each in the form of a two-dimensional matrix of pixels having digitized intensity values, an image display system, and an image processing system having access to the data of the sequence of images and to the display system, said image processing system including a device as claimed in claim 11.

References Cited (Backward Citations)

US4367490A Cited
US4663773A Cited
US4639784A Cited
US5018179A Cited
US5296937A Cited
US5309237A Cited
US5467380A Cited
US5684720A Cited

USPTO Annuity & Maintenance Windows

Annuity Stage Window Opens Payment Deadline Surcharge Period Until Estimated Fee Status
3.5 Years (Year 4) 2002-05-31 2002-11-30 2003-06-01 $2000.00 Window Closed / Lapsed
7.5 Years (Year 8) 2006-05-31 2006-11-30 2007-06-01 $3760.00 Window Closed / Lapsed
11.5 Years (Year 12) 2010-05-31 2010-11-30 2011-06-01 $7700.00 Window Closed / Lapsed

Assigned Categories & Tags

No categories assigned to this patent yet.

Assign Category Tag