Processing velocity information in an ultrasonic system

• US 5,628,321 A
• Filed: 12/18/1995
• Issued: 05/13/1997
• Est. Priority Date: 12/18/1995
• Status: Expired due to Term

• Current Assignee: Diasonics Ultrasound, Inc. (IDB Holding Corporation Limited)
• Original Assignee: Diasonics Ultrasound, Inc. (IDB Holding Corporation Limited)
• Inventors: Scheib, John P., Lin, Sheng-Tz
• Primary Examiner(s): MANUEL, GEORGE C
• Application Number:
  US08/573,698
• Time in Patent Office: 512 Days
• Field of Search: 128/661.08, 128/661.09, 128/661.1, 128/662.01, 73/861.18, 73/861.25
• US Class Current: 600/453
• 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: G01S 15/8979 Combined Doppler and pulse-echo imaging systems
• CPC Class Codes: G01S 7/52071 Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging

Abstract

A method and apparatus in an ultrasonic system for determining an optimum cardiac cycle of a living subject, and for performing vascular measurements. Peak velocities of a spectra of echoes over a period of time due to pulses emitted into a subject under examination are determined. Based upon the peak velocities of the spectra over the period of time, a time series of the peak velocities are determined representing an optimum cardiac cycle. In addition, vascular measurements of the peak velocities for the optimum cardiac cycle can be determined from the optimum cardiac cycle. The peak velocities of the spectra can be performed by dividing the spectra in a plurality of regions according to ranges of the velocity information. A threshold magnitude of the peak velocities is determined by adding a predetermined value to a minimum magnitude. The threshold magnitude is used to find a peak velocity by searching for the peak velocity from a Nyquist velocity of the spectra, identifying a first velocity which has an magnitude exceeding the threshold magnitude, and using the first velocity as the peak velocity. Determining the optimum cardiac cycle of the living subject includes setting a threshold velocity identifying each the peak velocities which exceed the threshold to form a set of candidate peak velocities. Discontinuities in the set of candidate peak velocities are determined, and the optimum cardiac cycle can be determined therefrom. Various vascular measurements for the optimum cardiac cycle include the peak systole (PS) end diastole (ED), include a minimum diastolic deflection (MDD), a rise time, an acceleration index, a pulsatility index, or a resistivity index.

Claims

1. A method of making vascular measurements within a living subject in an ultrasonic system comprising the following steps:

a. determining peak velocities of a spectra of echoes over a period of time due to pulses emitted into a subject under examination, determining said peak velocities includes dividing said spectra into a plurality of regions to determine a threshold magnitude of said peak velocities and identifying a first velocity that exceeds said threshold magnitude;

b. based upon said peak velocities of said spectra over said period of time, determining a time series of said peak velocities representing an optimum cardiac cycle for said subject under examination; and

c. based upon said time series of said peak velocities representing said optimum cardiac cycle, performing vascular measurements of said peak velocities for said optimum cardiac cycle.

2. The method of claim 1 wherein said spectra are obtained by emitting said pulses into said subject under examination, receiving said echoes and determining Doppler frequency shifts from said pulses, and translating said Doppler frequency shifts to velocity information.

3. The method of claim 1 wherein said step of determining said peak velocities of said spectra comprises:

a. dividing said spectra in a plurality of regions according to ranges of said velocity information;

b. determining an average magnitude of said velocity information for each of said plurality of regions;

c. determining a minimum magnitude of said average magnitudes for said plurality of regions;

d. determining a threshold magnitude of said peak velocities by adding a predetermined value to said minimum magnitude; and

e. using said threshold magnitude to find a peak velocity by searching for said peak velocity from a Nyquist velocity of said spectra, identifying a first velocity which has an magnitude exceeding said threshold magnitude, and using said first velocity as said peak velocity.

4. The method of claim 3 further comprising the step of repeating said steps a-e on positive and negative velocities of said spectra, determining a positive peak velocity, and a negative peak velocity, determining a greater absolute value of either said positive or said negative peak velocity and using said greater absolute value as said peak velocity.

5. The method of claim 3 wherein said step of determining said spectra comprises determining a Doppler spectrum, performing a wall-filter on said spectrum, and determining a power estimate of said spectrum.

6. The method of claim 3 further comprising the step of modifying a pulse repetition frequency (PRF) of said pulses if said peak velocity is at said Nyquist velocity in said spectra.

7. The method of claim 3 further comprising the step of averaging said peak velocity over several cardiac cycles of said subject under examination using said method.

8. The method of claim 3 wherein said step of determining said spectra comprises using display values for said ultrasonic system, and power values for each of said display values.

9. The method of claim 8 wherein said display values include color Doppler values representing velocities of a moving substance.

10. The method of claim 9 further comprising the step of averaging said peak velocity over several cardiac cycles of said subject under examination using said method.

11. The method of claim 1 wherein said step of determining said optimum cardiac cycle of said living subject is comprised of the steps of:

a. setting a threshold velocity;

b. identifying each said peak velocity which exceeds said threshold to form a set of candidate peak velocities, and associating a time of occurrence of each said peak velocity with each of said set of candidate peak velocities;

c. identifying discontinuities in said set of candidate peak velocities;

d. locating a candidate peak velocity between said each of said discontinuities to form a set of peak systoles;

e. determining a largest velocity of said set of peak systoles to determine a candidate peak systole;

f. determining whether there is a second peak systole forward in time from said candidate peak systole;

g. if there is not said second peak systole forward in time from said candidate peak systole then selecting a next largest velocity of said peak systoles, using said next largest velocity of said peak systoles as said candidate peak systole and repeating step f until there is a peak systole forward in time from said peak systole;

h. if there is said second peak systole forward in time from said candidate peak systole, then locating a first minimum velocity which occurs before said candidate peak systole and identifying said first minimum velocity as an end diastole (ED) for said optimum cardiac cycle;

i. locating a second minimum velocity which occurs before said second peak systole and identifying said second minimum velocity as a second end diastole (ED); and

j. identifying a time period between said end diastole and said second end diastole as said optimum cardiac cycle.

12. The method of claim 11 wherein said step of performing said vascular measurements for said optimum cardiac cycle includes determining a minimum diastolic deflection (MDD) for said optimum cardiac cycle.

13. The method of claim 11 wherein said step of performing said vascular measurements for said optimum cardiac cycle includes determining a rise time for said optimum cardiac cycle.

14. The method of claim 13 wherein said step of performing said vascular measurements for said optimum cardiac cycle includes determining an acceleration index for said optimum cardiac cycle.

15. The method of claim 13 wherein said step of determining said rise time for said optimum cardiac cycle includes successively determining differences between pairs of adjacent peak velocities between said end diastole and said peak systole, and upon detecting a sign change between a first difference and a second difference, identifying a first adjacent peak velocity as a first systolic peak velocity for said optimum cardiac cycle.

16. The method of claim 15 wherein said step of determining said rise time for said optimum cardiac cycle includes determining a difference in time between said first systolic peak velocity and said end diastole.

17. The method of claim 15 further including the step of determining an acceleration index for said optimum cardiac cycle by determining a difference between said first systolic peak velocity and said end diastole divided by a difference in time between said first systolic peak velocity and said end diastole.

18. The method of claim 11 wherein said step of performing said vascular measurements for said optimum cardiac cycle includes determining a pulsatility index for said optimum cardiac cycle.

19. The method of claim 11 wherein said step of performing said vascular measurements for said optimum cardiac cycle includes determining a resistivity index for said optimum cardiac cycle.

20. An apparatus for making vascular measurements within a living subject in an ultrasonic system comprising:

a. a first circuit for determining peak velocities of a spectra of echoes over a period of time due to pulses emitted into a subject under examination and storing said peak velocities into an array of peak velocities, said first circuit is further configured to divide said spectra into a plurality of regions to determine a threshold magnitude of said peak velocities and identify a first velocity that exceeds said threshold magnitude;

b. a second circuit coupled to said first circuit for determining an optimum cardiac cycle for said subject under examination based upon said time series of said peak velocities; and

c. a third circuit coupled to said second circuit for performing vascular measurements of said optimum cardiac cycle.

21. The apparatus of claim 20 further comprising a fourth circuit coupled to said first circuit for emitting said pulses into said subject under examination, receiving said echoes and determining Doppler frequency shifts from said pulses, and translating said Doppler frequency shifts to velocity information to form said spectra of echoes.

22. The apparatus of claim 20 wherein said first circuit includes:

a. a division means for dividing said spectra in a plurality of regions according to ranges of said velocities;

b. an averaging means coupled to said division means for determining an average magnitude of said velocity information for each of said plurality of regions;

c. a minimum determination means coupled to said averaging means for determining a minimum magnitude of said average magnitudes for said plurality of regions;

d. a threshold determination means coupled to said minimum determination means for determining a threshold magnitude of said peak velocities by adding a predetermined value to said minimum magnitude; and

e. a thresholding means coupled to said threshold determination means for using said threshold magnitude to find a peak velocity by searching for said peak velocity from a Nyquist velocity of said spectra, identifying a first velocity which has an magnitude exceeding said threshold magnitude, and using said first velocity as said peak velocity.

23. The apparatus of claim 22 further comprising:

a. an activation circuit for activating said division means, said averaging means, said minimum determination means, said threshold determination means, and said thresholding means on positive and negative velocities of said spectra; and

b. a peak detection circuit coupled to said activation circuit for determining a greater absolute value of either a positive or a negative peak velocity and using said greater absolute value as said peak velocity.

24. The apparatus of claim 22 wherein said first circuit includes a signal processing circuit for determining a Doppler spectrum, performing a wall-filter on said spectrum, and determining a power estimate of said spectrum.

25. The apparatus of claim 22 further comprising an adjustment circuit for modifying a pulse repetition frequency (PRF) of said pulses if said peak velocity is at a Nyquist point in said spectra.

26. The apparatus of claim 22 further comprising a temporal averaging means coupled to said thresholding means for averaging said peak velocity over several cardiac cycles of said subject under examination.

27. The apparatus of claim 22 wherein said first circuit includes a input means for using display values for said ultrasonic system to determine said peak velocities, and power values for each of said display values.

28. The apparatus of claim 27 wherein said display values include color Doppler values representing velocities of a moving substance.

29. The apparatus of claim 28 further comprising a temporal averaging means coupled to said thresholding means for averaging said peak velocity over several cardiac cycles of said subject under examination.

30. The apparatus of claim 20 wherein said second circuit includes:

a. a threshold determination means for setting a threshold velocity;

b. a thresholding means coupled to said threshold determination means for identifying each said peak velocity which exceeds said threshold to form a set of candidate peak velocities which are stored in an array of candidate velocities, and associating a time of occurrence of each said peak velocity with each of said set of candidate peak velocities;

c. an identification means coupled to said array of candidate velocities for identifying discontinuities in said set of candidate peak velocities;

p1 d. a candidate peak detection means coupled to said identification means for locating a candidate peak velocity between said each of said discontinuities to form a set of peak systoles in an array of peak systoles;

e. a largest velocity detection means coupled to said array of peak systoles for determining a largest velocity of said set of peak systoles to determine a candidate peak systole;

f. a determining means coupled to said array of peak systoles for determining whether there is a second peak systole forward in time from said candidate peak systole and if there is not said second peak systole forward in time from said candidate peak systole then selecting a next largest velocity of said peak systoles, activating said largest velocity detection means for said next largest velocity of said peak systoles as said candidate peak systole until there is said second peak systole forward in time from said candidate peak systole;

g. a minimum velocity detection means coupled to said array of peak velocities for locating a first minimum velocity which occurs before said candidate peak systole and identifying said first minimum velocity as an end diastole (ED) for said optimum cardiac cycle;

h. a second minimum velocity detection means coupled to said array of peak velocities for locating a second minimum velocity which occurs before said second peak systole and identifying said second minimum velocity as a second end diastole (ED); and

i. an optimum cardiac cycle identifying means coupled to said minimum velocity detection means and said second minimum velocity detection means for identifying a time period between said end diastole and said second end diastole as said optimum cardiac cycle.

31. The apparatus of claim 30 wherein said third circuit includes a minimum diastolic deflection determination circuit for determining a minimum diastolic deflection (MDD) for said optimum cardiac cycle.

32. The apparatus of claim 30 wherein said third circuit includes a rise time determination circuit for determining a rise time for said optimum cardiac cycle.

33. The apparatus of claim 32 wherein said third circuit includes an acceleration index determination means for determining an acceleration index for said optimum cardiac cycle.

34. The apparatus of claim 32 wherein said rise time determination circuit includes:

a. successive determination means for successively determining differences between pairs of adjacent peak velocities between said end diastole and said peak systole; and

b. a detection means coupled to said successive determination means for detecting a sign change between a first difference and a second difference, and identifying a first adjacent peak velocity as a first systolic peak velocity for said optimum cardiac cycle.

35. The apparatus of claim 34 wherein said rise time determination circuit further includes a difference circuit coupled to said detection means for determining a difference in time between said first systolic peak velocity and said end diastole.

36. The apparatus of claim 35 wherein said rise time determination circuit further includes an acceleration index determination circuit coupled to said detection means for determining an acceleration index for said optimum cardiac cycle by determining a difference between said first systolic peak velocity and said end diastole divided by a difference in time between said first systolic peak velocity and said end diastole.

37. The apparatus of claim 30 wherein said third circuit includes a pulsatility index determination means for determining a pulsatility index for said optimum cardiac cycle.

38. The apparatus of claim 30 wherein said third circuit includes a resistivity index determination means for determining a resistivity index for said optimum cardiac cycle.

39. A method in an ultrasonic system comprising the following steps:

a. determining peak velocities of a spectra of echoes over a period of time due to pulses emitted into a subject under examination, determining said peak velocities further includes dividing said spectra into a plurality of regions to determine a threshold magnitude of said peak velocities and identifying a first velocity that exceeds said threshold magnitude; and

b. based upon said peak velocities of said spectra over said period of time, determining a set of said peak velocities representing and optimum cardiac cycle for said subject under examination.

40. The method of claim 39 wherein said spectra are obtained by emitting said pulses into said subject under examination, receiving said echoes and determining Doppler frequency shifts from said pulses, and translating said Doppler frequency shifts to velocity information.

41. The method of claim 40 wherein said step of determining said peak velocities of said spectra comprises:

a. dividing said spectra in a plurality of regions according to ranges of said velocity information;

b. determining an average magnitude of said velocity information for each of said plurality of regions;

c. determining a minimum magnitude of said average magnitudes for said plurality of regions;

d. determining a threshold magnitude of said peak velocities by adding a predetermined value to said minimum magnitude; and

e. using said threshold magnitude to find a peak velocity by searching for said peak velocity from a Nyquist velocity of said spectra, identifying a first velocity which has an magnitude exceeding said threshold magnitude, and using said first velocity as said peak velocity.

42. The method of claim 41 further comprising the step of repeating said steps a-e on positive and negative velocities of said spectra, determining a positive peak velocity, and a negative peak velocity, determining a greater absolute value of either said positive or said negative peak velocity and using said greater absolute value as said peak velocity.

43. The method of claim 42 wherein said step of determining said spectra comprises determining a Doppler spectrum, performing a wall-filter on said spectrum, and determining a power estimate of said spectrum.

44. The method of claim 43 further comprising the step of modifying a pulse repetition frequency (PRF) of said pulses if said peak velocity is at said Nyquist velocity in said spectra.

45. The method of claim 39 further comprising the step of averaging said peak velocity over several cardiac cycles of said subject under examination using said method.

46. The method of claim 39 wherein said step of determining said spectra comprises using display values for said ultrasonic system, and power values for each of said display values.

47. The method of claim 46 wherein said display values include color Doppler values representing velocities of a moving substance.

48. The method of claim 47 further comprising the step of averaging said peak velocity over several cardiac cycles of said subject under examination using said method.

49. An apparatus in an ultrasonic system comprising:

a. a first circuit for determining peak velocities of a spectra of echoes over a period of time due to pulses emitted into a subject under examination and storing said peak velocities into an array of peak velocities, said first circuit is further configured to divide said spectra into a plurality of regions to determine a threshold magnitude of said peak velocities and identify a first velocity that exceeds said threshold magnitude; and

b. a second circuit coupled to said first circuit for determining an optimum cardiac cycle for said subject under examination based upon said time series of said peak velocities.

50. The apparatus of claim 49 further comprising a third circuit coupled to said first circuit for emitting said pulses into said subject under examination, receiving said echoes and determining Doppler frequency shifts from said pulses, and translating said Doppler frequency shifts to velocity information to form said spectra of echoes.

51. The apparatus of claim 50 further comprising an adjustment circuit for modifying a pulse repetition frequency (PRF) of said pulses if said peak velocity is at a Nyquist point in said spectra.

52. The apparatus of claim 49 wherein said first circuit includes:

a. a division means for dividing said spectra in a plurality of regions according to ranges of said velocities;

b. an averaging means coupled to said division means for determining an average magnitude of said velocity information for each of said plurality of regions;

c. a minimum determination means coupled to said averaging means for determining a minimum magnitude of said average magnitudes for said plurality of regions;

d. a threshold determination means coupled to said minimum determination means for determining a threshold magnitude of said peak velocities by adding a predetermined value to said minimum magnitude; and

e. a thresholding means coupled to said threshold determination means for using said threshold magnitude to find a peak velocity by searching for said peak velocity from a Nyquist velocity of said spectra, identifying a first velocity which has an magnitude exceeding said threshold magnitude, and using said first velocity as said peak velocity.

53. The apparatus of claim 52 further comprising:

a. an activation circuit for activating said division means, said averaging means, said minimum determination means, said threshold determination means, and said thresholding means on positive and negative velocities of said spectra; and

b. a peak detection circuit coupled to said activation circuit for determining a greater absolute value of either a positive or a negative peak velocity and using said greater absolute value as said peak velocity.

54. The apparatus of claim 52 wherein said first circuit includes a signal processing circuit for determining a Doppler spectrum, performing a wall-filter on said spectrum, and determining a power estimate of said spectrum.

55. The apparatus of claim 52 further comprising a temporal averaging means coupled to said thresholding means for averaging said peak velocity over several cardiac cycles of said subject under examination.

56. The apparatus of claim 52 wherein said first circuit includes a input means for using display values for said ultrasonic system to determine said peak velocities, and power values for each of said display values.

57. The apparatus of claim 56 wherein said display values include color Doppler values representing velocities of a moving substance.

58. The apparatus of claim 57 further comprising a temporal averaging means coupled to said thresholding means for averaging said peak velocity over several cardiac cycles of said subject under examination.

59. A method in an ultrasonic system of determining an optimum cardiac cycle comprising the following steps:

a. determining a time series of velocities over a period of time due to pulses emitted into a subject under examination;

b. identifying discontinuities in said time series of peak velocities which exceeds a threshold;

c. determining a candidate peak velocity between said each of said discontinuities to form a set of peak systoles;

d. determining a peak systole in said set of peak systoles which is followed by a second peak systole and which is preceded by a minimum peak velocity which occurs within a predetermined time before said peak systole and identifying said minimum peak velocity as an end diastole (ED);

e. locating a second minimum velocity which occurs before said second peak systole and identifying said second minimum velocity as a second end diastole (ED); and

f. identifying a time period between said end diastole and said second end diastole as said optimum cardiac cycle.

60. The method of claim 59 further comprising the step of determining a minimum diastolic deflection (MDD) for said optimum cardiac cycle.

61. The method of claim 59 further comprising the step of determining a rise time for said optimum cardiac cycle.

62. The method of claim 61 further comprising the step of determining an acceleration index for said optimum cardiac cycle.

63. The method of claim 59 further comprising the step of determining a pulsatility index for said optimum cardiac cycle.

64. The method of claim 59 further comprising the step of determining a resistivity index for said optimum cardiac cycle.

65. The method of claim 64 further comprising the steps of successively determining differences between pairs of adjacent peak velocities between said end diastole and said peak systole, and upon detecting a sign change between a first difference and a second difference, identifying a first adjacent peak velocity as a first systolic peak velocity for said optimum cardiac cycle.

66. The method of claim 65 further comprising the step of determining a difference in time between said first systolic peak velocity and said end diastole.

67. The method of claim 65 further including the step of determining an acceleration index for said optimum cardiac cycle by determining a difference between said first systolic peak frequency and said end diastole divided by a difference in time between said first systolic peak velocity and said end diastole normalized by the transmit frequency.

68. An apparatus for determining an optimum cardiac cycle comprising:

a. an input means for determining a time series of velocities over a period of time due to pulses emitted into a subject under examination;

b. an identification means coupled to said input means for identifying discontinuities in said time series of peak velocities which exceeds a threshold;

c. a peak detection means coupled to said identification means for determining a peak velocity between said each of said discontinuities to form a set of peak systoles;

d. a determination means coupled to said peak detection means for determining a peak systole in said set of peak systoles which is followed by a second peak systole and which is preceded by a minimum peak velocity which occurs within a predetermined time before said peak systole and identifying said minimum peak velocity as an end diastole (ED);

e. a second minimum velocity detection means coupled to said determination means for locating a second minimum velocity which occurs before said second peak systole and identifying said second minimum velocity as a second end diastole (ED); and

f. an optimum cardiac cycle identifying means coupled to said minimum velocity detection means for identifying a time period between said end diastole and said second end diastole as said optimum cardiac cycle.

69. The apparatus of claim 68 further comprising a minimum diastolic deflection determination circuit coupled to said optimum cardiac cycle identifying means for determining a minimum diastolic deflection (MDD) for said optimum cardiac cycle.

70. The apparatus of claim 68 further comprising a rise time determination means coupled to said optimum cardiac cycle identifying means for determining a rise time for said optimum cardiac cycle.

71. The apparatus of claim 70 wherein said rise time determination means includes:

a. successive determination means for successively determining differences between pairs of adjacent peak velocities between said end diastole and said peak systole; and

b. a detection means coupled to said successive determination means for detecting a sign change between a first difference and a second difference, and identifying a first adjacent peak velocity as a first systolic peak velocity for said optimum cardiac cycle.

72. The apparatus of claim 71 wherein said rise time determination circuit further includes a difference circuit coupled to said detection means for determining a difference in time between said first systolic peak velocity and said end diastole.

73. The apparatus of claim 70 wherein said rise time determination means further includes an acceleration index determination circuit coupled to said detection means for determining an acceleration index for said optimum cardiac cycle by determining a difference between said first systolic peak frequency and said end diastole divided by a difference in time between said first systolic peak velocity and said end diastole and normalized by the transmit frequency.

74. The apparatus of claim 68 further comprising an acceleration index determination means coupled to said optimum cardiac cycle identifying means for determining an acceleration index for said optimum cardiac cycle.

75. The apparatus of claim 68 further comprising a pulsatility index determination means coupled to said optimum cardiac cycle identifying means for determining a pulsatility index for said optimum cardiac cycle.

76. The apparatus of claim 68 a resistivity index determination means for determining a resistivity index for said optimum cardiac cycle.