Ultrasound clutter filtering with iterative high pass filter selection

• US 6,997,876 B2
• Filed: 04/16/2004
• Issued: 02/14/2006
• 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/825,719
• Publication Number:
  US 20040199078A1
• Time in Patent Office: 669 Days
• Field of Search: 600453-456, 738/612.5
• US Class Current: 600/455
• 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

A system and method for ultrasound clutter filtering is provided. A processor is configured to iteratively select an optimal high pass filter for the progressive, ordered filtering of clutter from ultrasound color flow imaging data. The high pass filter input for each iterative selection and ordered set of high pass filters is the same original ultrasound color flow imaging data. The high pass filters have different cutoff frequencies whereby each high pass filter can be implemented using different structures. The system and method allow for filtering of clutter from ultrasound color flow imaging data until the clutter is substantially removed.

Claims

1. An ultrasound clutter filter system, comprising:

a processor configured to iteratively select an optimal high pass filter for progressive, ordered filtering of clutter from ultrasound color flow imaging data wherein a determination of whether a high pass filter is optimal is made by comparing an index to a threshold for that index, the index being computed using a mathematical formula including a mean frequency and a magnitude of a filtered signal, wherein a high pass filter input for each iterative selection is the original ultrasound color flow imaging data.

2. The system of claim 1 wherein the processor is configured for selecting said optimal high pass filter for performing filtering over each packet of flow data that corresponds to an acoustic point in a color flow region of interest.

3. The system of claim 1 wherein the processor is configured to select the optimal high pass filter is if the filtered signal is less than a threshold.

4. The system of claim 1 wherein the processor is configured to select the optimal high pass filter is if a magnitude of a portion of the filtered signal that is within a predetermined frequency range is less than a threshold.

5. The system of claim 1 wherein the processor is configured to select the optimal high pass filter is if a total energy of the filtered signal is less than an energy threshold.

6. The system of claim 1 wherein the processor is configured for selecting said optimal high pass filter for performing filtering in real time while imaging.

7. The system of claim 1 further comprising a finite number of high pass filters from which the optimal high pass filter is selected.

8. The system of claim 7 in which the finite number of high pass filters is two.

9. The system of claim 7, wherein the finite number of high pass filters comprises a finite impulse response filter.

10. The system of claim 7, wherein the finite number of high pass filters comprises an infinite impulse response filter.

11. An ultrasound clutter filter system, comprising:

a processor configured to iteratively select an optimal high pass filter for progressive, ordered filtering of clutter from ultrasound color flow imaging data wherein a criterion for selecting the optimal high pass filter is if a mean frequency of filtered signal data is less than a clutter frequency threshold wherein if the mean frequency is less than the clutter frequency threshold is determined by whether an absolute value of an imaginary part of a first order autocorrelation of the filtered signal data is less than a constant times a real part of the autocorrelation, where the constant is determined by the clutter frequency threshold, wherein a high pass filter input for each iterative selection is the original ultrasound color flow imaging data.

12. The system of claim 11 wherein the processor is further configured to select the optimal high pass filter further comprises if the filtered signal data is less than a threshold.

13. The system of claim 11 wherein the processor is further configured to select the optimal high pass filter further comprises if a magnitude of a portion of the filtered signal that is within a predetermined frequency range is less than a threshold.

14. The system of claim 11 wherein the processor is further configured to select the optimal high pass filter further comprises if a total energy of the filtered signal is less than an energy threshold.

15. The system of claim 11 wherein the filtering high pass filter is configured to perform over each packet of flow data that corresponds to an acoustic point in a color flow region of interest.

16. The system of claim 11 wherein the filtering high pass filter is configured to perform in real time while imaging.

17. The system of claim 11 further comprising a finite number of high pass filters from which the optimal high pass filter is selected.

18. The system of claim 17 in which the finite number of high pass filters is two.

19. An ultrasound clutter filtering method, comprising:

iteratively selecting an optimal high pass filter for progressive, ordered filtering wherein the iteratively selecting comprises computing an index using a mathematical formula including a mean frequency and a magnitude of a filtered signal and determining whether a filter is optimal by comparing the index to a threshold, wherein a high pass filter input for each iterative selection is original ultrasound color flow data; and

filtering clutter from the ultrasound color flow data until the clutter is substantially removed.

20. The method of claim 19 wherein the filtering comprises filtering an individual packet of the ultrasound color flow data that corresponds to an acoustic point in a region of interest.

21. The method of claim 19 wherein the filtering is performed in real time while collecting the ultrasound color flow data.

22. The method of claim 19 wherein the iteratively selecting further comprises determining if the filtered signal is less than a threshold.

23. The method of claim 19 wherein the iteratively selecting further comprises determining if a magnitude of the filtered signal is less than a than a frequency threshold.

24. The method of claim 19 wherein the iteratively selecting further comprises selecting a high pass filter from a finite number of high pass filters.

25. The method of claim 19 wherein the iteratively selecting further comprises selecting one of two high pass filters.

26. The method of claim 25, wherein one of the two high pass filters comprises a finite impulse response filter.

27. The method of claim 25, wherein one of the two high pass filters comprises an infinite impulse response filter.

28. The method of claim 19 wherein the iteratively selecting further comprises determining if a magnitude of a color flow signal in a preselected frequency range is less than a color flow signal threshold wherein a phase shift is less than a phase shift threshold.

29. An ultrasound clutter filtering method, comprising:

iteratively selecting an optimal high pass filter for progressive, ordered filtering of ultrasound color flow data wherein the iteratively selecting comprises determining if a magnitude of a color flow signal in a preselected frequency range is less than a color flow signal threshold, wherein the determining includes determining whether an absolute value of an imaginary part of a first order autocorrelation of the color flow data is less than a constant times a real part of the autocorrelation, where the constant is determined by a frequency threshold, wherein a high pass filter input for each iterative selection is the original ultrasound color flow data, and;

filtering clutter from the ultrasound color flow data until the clutter is substantially removed.

30. The method of claim 29 wherein the filtering comprises filtering an individual packet of the ultrasound color flow data that corresponds to an acoustic point in a region of interest.

31. The method of claim 29 wherein the filtering is performed in real time while collecting the ultrasound color flow data.

32. The method of claim 29 wherein the iteratively selecting further comprises determining if the filtered signal is less than a threshold.

33. The method of claim 29 wherein the iteratively selecting further comprises determining if a magnitude of frequency of the filtered signal is less than a frequency threshold.

34. The method of claim 29 wherein determining further comprises determining whether a phase shift is less than a phase shift threshold.

35. The method of claim 29 wherein the iteratively selecting further comprises selecting a high pass filter from a finite number of high pass filters.

36. The method of claim 29 wherein the iteratively selecting further comprises selecting one of two high pass filters.