System and method for adaptive clutter filtering in ultrasound color flow imaging

• US 6,733,455 B2
• Filed: 06/11/2002
• Issued: 05/11/2004
• Est. Priority Date: 08/20/1999
• Status: Expired due to Term

• Current Assignee: Shenzhen Mindray Bio-Medical Electronics Co., Ltd. (Excelsior Union Ltd.)
• Original Assignee: Zonare Medical Systems Incorporated (Excelsior Union Ltd.)
• Inventors: Ji, Ting-Lan, Mo, Larry Y. L., McLaughlin, Glen W., Chou, Ching-Hua
• Primary Examiner(s): Jaworski, Francis J.
• Application Number:
  US10/167,606
• Publication Number:
  US 20020169378A1
• Time in Patent Office: 700 Days
• Field of Search: 600/453-456, 738/612.5, 361/135, 382/162, 382/261, 382/263, 382/265, 382/272
• US Class Current: 600/454
• CPC Class Codes: A61B 8/00 Diagnosis using ultrasonic, sonic or infrasonic waves imaging of objects using sonar G01S15/00
• CPC Class Codes: A61B 8/06 Measuring blood flow measuring volume flow in general G01F, e.g. G01F1/66, G01F1/72; measuring speed of fluids in general G01P5/00
• CPC Class Codes: A61B 8/08 Detecting organic movements or changes, e.g. tumours, cysts, swellings A61B8/02 - A61B8/06 take precedence
• CPC Class Codes: A61B 8/0833 involving detecting or locating foreign bodies or organic structures
• CPC Class Codes: A61B 8/0866 involving foetal diagnosis; pre-natal or peri-natal diagnosis of the baby
• CPC Class Codes: A61B 8/13 Tomography A61B8/10, A61B8/12 take precedence; tomography for radiation diagnosis A61B6/02
• CPC Class Codes: A61B 8/4433 involving a docking unit
• CPC Class Codes: A61B 8/4438 Means for identifying the diagnostic device, e.g. barcodes
• CPC Class Codes: A61B 8/4455 Features of the external shape of the probe, e.g. ergonomic aspects
• CPC Class Codes: A61B 8/462 characterised by constructional features of the display
• CPC Class Codes: A61B 8/467 characterised by special input means
• CPC Class Codes: A61B 8/488 involving Doppler signals
• CPC Class Codes: G01S 15/8981 Discriminating between fixed and moving objects or between objects moving at different speeds, e.g. wall clutter filter
• CPC Class Codes: G01S 7/52084 related to particular user interfaces special user input means for ultrasonic medical diagnostic devices A61B8/467

Abstract

An adaptive clutter filtering for ultrasound color flow imaging is provided including an iterative algorithm that is used to select the best clutter filter for each packet of color flow data. If significant clutter motion is present, a high pass filter cutoff frequency is automatically set to suppress the clutter and associated flash artifacts. The cutoff frequency is chosen according to the frequency of the clutter—the lower the clutter frequency, the lower the cutoff frequency can be. If clutter frequencies are low, lower filter cutoffs allow for maximum low flow detection. In this manner, the filter cutoff frequency can be optimized based on the data for each pixel in the color flow image.

Claims

1. A system comprising:

a processor configured to iteratively select an optimal filter from an ordered set of filters for filtering out clutter from ultrasound image data while collecting ultrasound image data, using as criteria for selecting the optimal filter

whether a filtered signal is less than a signal threshold, and

whether a magnitude of frequency of the filtered signal is less than a frequency threshold, which is determined by whether a phase shift is less than a phase shift threshold, which in turn is determined by whether a real part of a first order autocorrelation of the filtered signal is greater than a constant times an imaginary part of the autocorrelation, where the constant is determined by the frequency threshold.

2. An adaptive clutter filtering method for color flow imaging that selects an optimal set of High Pass Filter (HPF) parameters based on an In phase component and Quadrature component of a data packet (I/Q data) for each pixel location in an image, the adaptive clutter filtering method having an iterative loop comprising:

applying an initial set of HPF parameters to the I/Q data to produce filtered I/Q data;

computing a magnitude related index and a phase change related index of the filtered I/Q data;

testing the magnitude related index and phase change related index using corresponding threshold tests, based on predefined threshold constants, to determine if clutter is present in the filtered I/Q data; and

deciding whether to

repeat the applying, the computing, the testing and the deciding using another set of HPF parameters, or

exit the iterative loop and compute flow related parameters for color flow image formation.

3. The adaptive clutter filtering method of claim 2 wherein the optimal set of HPF parameters is selected from a group of sets of HPF parameters ordered in a sequence according to a measure of an amount of clutter rejected.

4. The adaptive clutter filtering method of claim 3, wherein the iterative loop sequentially searches the group of sets of HPF parameters for the optimal set of HPF parameters according to the measure, and exits when either a clutter status change occurs, or when a last set of HPF parameters in the sequence has been reached.

5. The adaptive clutter filtering method of claim 4, wherein the iterative loop performs a binary search for the optimal set of HPF parameters from the group of sets of HPF parameters, starting from a middle of the sequence, and stopping when either a clutter status change occurs or when an extreme set of HPF parameters is reached.

6. The adaptive clutter filtering method of claim 4, wherein selecting a set of HPF parameters corresponds to changing a passband frequency cutoff value and sharpness of a frequency transition band.

7. The adaptive clutter filtering method of claim 4, wherein computing the magnitude related index includes computing a total power of the filtered I/Q data.

8. The adaptive clutter filtering method of claim 4, wherein computing the magnitude related index includes computing a sum of absolute values of the In phase component and of the Quadrature component of the filtered I/Q data.

9. The adaptive clutter filtering method of claim 4, wherein computing the phase change related index includes computing a real part and an imaginary part of an autocorrelation of the filtered I/Q data.

10. The adaptive clutter filtering method of claim 4, wherein the magnitude threshold test includes comparing the magnitude related index to a predefined clutter magnitude threshold that is a function of an iteration number.

11. The adaptive clutter filtering method of claim 4, wherein the magnitude threshold test includes comparing the magnitude related index to a predefined clutter magnitude threshold that is a function of a currently selected set of HPF parameters.

12. The adaptive clutter filtering method of claim 4, wherein the testing includes comparing |y | with C*x, where x and y denote the real and imaginary parts of an autocorrelation of the In phase and Quadrature components of the filtered I/Q data, and C is a predefined constant that corresponds to the largest possible phase change per pulse repetition period for clutter.

13. The adaptive clutter filtering method of claim 4, further comprising:

processing a predetermined percentage of pixels in a color flow region of interest;

after the processing, adjusting the initial set of HPF parameters based on statistics of previously selected sets of HPF parameters; and

processing pixels of the color flow region of interest not previously processed based on the adjusting.

14. The adaptive clutter filtering method of claim 4, wherein a maximum number of iterations allowed for selecting a set of HPF parameters is larger in a predefined central portion of a color flow region of interest, and smaller in other portions.

15. A method comprising:

collecting ultrasound data;

iteratively selecting an optimal filter from an ordered set of filters for filtering out clutter from an ultrasound image signal while collecting the ultrasound data, including

filtering the ultrasound image signal to produce a filtered ultrasound image signal, and

determining

whether the filtered ultrasound image signal is less than a signal threshold, and

whether a magnitude of frequency is less than a frequency threshold, by determining whether a phase shift is less than a phase shift threshold, which in turn is determined by whether a real part of a first order autocorrelation of the filtered ultrasound signal is greater than a constant times an imaginary part of the autocorrelation, where the constant is determined by the frequency threshold.

16. A system comprising:

means for iteratively selecting an optimal filter from an ordered set of filters for filtering out clutter from ultrasound image signal while collecting the ultrasound data, including

means for performing an initial filtering of the ultrasound image signal to produce a filtered ultrasound image signal, and

means for determining

whether the ultrasound image signal is less than a signal threshold, and

whether a magnitude of frequency is less than a frequency threshold, having means for determining whether a phase shift is less than a phase shift threshold having means for determining whether a real part of a first order autocorrelation of the ultrasound image signal is greater than a constant times an imaginary part of the autocorrelation, where the constant is determined by the frequency threshold.