Method and apparatus for real-time calculation and display of strain in ultrasound imaging

• US 6,099,471 A
• Filed: 10/07/1998
• Issued: 08/08/2000
• Est. Priority Date: 10/07/1997
• Status: Expired due to Term

• Current Assignee: GE Vingmed Ultrasound AS (GE Healthcare Technologies, Inc.)
• Original Assignee: General Electric Company
• Inventors: Torp, Hans, Olstad, Bjorn, Heimdal, Andreas, Kristofferson, Kjell
• Primary Examiner(s): Jaworski, Francis J.
• Application Number:
  US09/167,896
• Time in Patent Office: 671 Days
• Field of Search: 600/438, 600/440, 600/441-443, 600/454-456
• US Class Current: 600/438
• 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: B23P 19/04 for assembling or disassembling parts
• CPC Class Codes: B23Q 17/20 for indicating or measuring workpiece characteristics, e.g. contour, dimension, hardness
• CPC Class Codes: G01S 15/8979 Combined Doppler and pulse-echo imaging systems
• CPC Class Codes: G01S 7/52042 determining elastic properties of the propagation medium or of the reflective target diagnostic techniques involving the measurement of strain A61B8/485

Abstract

An ultrasound system and method for calculation and display of strain velocity in real time is disclosed. Strain velocity may be determined from tissue velocity. Tissue velocity determined by measuring the pulse-to-pulse Doppler shift at range positions along an ultrasound beam and calculating tissue velocity based on the Doppler shift. The strain velocity is then calculated as a gradient of tissue velocity. Alternatively, linear regression methods may be used to calculate strain velocity from tissue velocity. Alternatively, strain velocity may be determined directly from the difference in Doppler shift between range positions. The result of the strain velocity determination may be displayed in a number of manners such as M-mode, color-coded video images or time-variation curves, and may be displayed in combination, or as a mixture, with a B-mode image. A reliability index may be calculated and used to modify the display of strain velocity information.

Claims

1. A method for providing real-time strain velocity images of a spatial region comprising:

firing ultrasound pulses into an area of interest of the spatial region;

acquiring echo signals for a plurality of range positions along ultrasonic beams in an area of interest using a beam interleaving technique to acquire echo information along different beam directions during a time interval between two consecutive ultrasound pulse firings in a common direction;

estimating, in real-time, strain velocities from said acquired echo signals for said range positions inside said spatial region; and

displaying estimated strain velocities for each range position at spatial coordinates on a display unit associated with said spatial region to provide a real-time image of said strain velocities for said spatial region.

2. The method according to claim 1, wherein the step of estimating strain velocities comprises:

estimating tissue velocity for range positions along the ultrasonic beam based on the echo signals; and

calculating the strain velocity as a spatial derivative of the tissue velocity.

3. The method according to claim 2 wherein the spatial derivative is found with a linear regression of the tissue velocity for range positions along the ultrasonic beam.

4. The method according to claim 3, wherein the step of estimating strain velocities comprises:

estimating tissue velocity for range positions along the ultrasonic beam based on the echo signals; and

calculating the strain velocity by determining a velocity difference between estimated tissue velocities associated with at least a first and second range positions and dividing the velocity difference by a distance between the first and second range positions.

5. The method according to claim 1, wherein the step of estimating the strain velocity comprises:

estimating a complex pulse-to-pulse correlation R(r) for a number of range positions along the ultrasonic beam based on the echo signals;

determining a strain correlation function, S(r), over a radial distance dr according to an equation S(r)=conj(R(r))*R(r+dr); and

calculating the strain velocity according to an equation SV(r)=c/(4π

drTfo) phase (S(r)).

6. The method according to claim 1, wherein the step of estimating the strain velocity comprises:

estimating a complex pulse-to-pulse correlation for a number of range positions along the ultrasonic beam, based on the echo signals;

calculating a strain correlation function from at least two range positions separated by a given radial distance; and

calculating the strain velocity based on the phase of the strain correlation function.

7. The method according to claim 6 wherein the strain correlation function is given by multiplying the conjugate of the complex pulse-to-pulse correlation for a first range position by the complex pulse-to-pulse correlation for a second range position where said second range position is located the given radial distance from said first range position.

8. The method according to claim 7 wherein the strain velocity is given by dividing a numerator defined as the product of the phase angle of the strain correlation function and the speed of sound by a denominator defined as the product of 4, π

, the given radial distance, the ultrasound frequency and the time between consecutive pulses of said multiple of pulses.

9. The method according to claim 6 further comprising the step of temporally averaging the strain correlation function before calculating the strain velocity.

10. The method according to claim 6 further comprising the step of spatially averaging the strain correlation function before calculating the strain velocity.

11. The method according to claim 1 wherein the step of displaying the estimated strain velocity comprises combining a B-mode tissue image and the strain velocities at corresponding spatial coordinates in a color encoding.

12. The method according to claim 1 wherein the step of displaying the estimated strain velocity comprises combining an M-mode tissue image and the strain velocities at corresponding spatial coordinates in a color encoding.

13. The method according to claims 1 wherein the step of displaying the estimated strain velocities includes displaying the estimated strain velocity for a complete two-dimensional region.

14. The method according to claim 1 wherein the estimated strain velocities are estimated for a three-dimensional region and displayed with a three-dimensional visualization technique.

15. The method according to claim 1 wherein multiple line acquisition techniques are applied to further increase the resulting frame-rate.

16. The method according to claim 1 further comprising the steps of:

calculating a strain reliability index as a function of signal power for each sample point inside the said spatial region; and

modifying the display of the strain velocity according to the strain reliability index.

17. The method according to claim 1 wherein the strain reliability index given by dividing the absolute value of the product of the conjugate of the complex pulse-to-pulse correlation for a first range position and the complex pulse-to-pulse correlation for a second range position by the signal power at one of said first or second range positions.

18. The method according to claim 1 wherein a harmonic of the echo signals is used to estimate the strain velocities.

19. The method according to claim 1 wherein the strain velocity images are generated simultaneously with the application of an external force to produce a pressure gradient in the spatial region.

20. The method according to claim 1 wherein the strain velocity images are generated while simultaneously applying changing pressure against the imaged spatial region.

21. The method according to claim 1 wherein the echo signals are acquired using intravascular ultrasound techniques to assess regional changes in the vessel wall properties.

22. A method for generating real-time strain velocity images comprising:

acquiring harmonic signals for echo signals received from a plurality of range positions along an ultrasonic beam in an area of interest to cover a spatial region;

estimating in real-time strain velocities for said range positions inside said spatial region based on the harmonic signals for the acquired echo signals; and

displaying estimated strain velocities for each range position at spatial coordinates on a display unit associated with said spatial region to provide a real-time image of said strain velocities for said spatial region.

23. The method according to claim 22, wherein the step of estimating strain velocities comprises:

estimating tissue velocity for range positions along the ultrasonic beam based on the echo signals; and

calculating the strain velocity as a spatial derivative of the tissue velocity.

24. The method according to claim 23 wherein the spatial derivative is found with a linear regression of the tissue velocity for range positions along the ultrasonic beam.

25. The method according to claim 22, wherein the step of estimating strain velocities comprises:

estimating tissue velocity for range positions along the ultrasonic beam based on the echo signals; and

calculating the strain velocity by determining a velocity difference between estimated tissue velocities associated with at least a first and second range positions and dividing the velocity difference by a distance between the first and second range positions.

26. The method according to claim 22, wherein the step of estimating the strain velocity comprises:

estimating a complex pulse-to-pulse correlation R(r) for a number of range positions along the ultrasonic beam based on the echo signals;

determining a strain correlation function, S(r), over a radial distance dr according to an equation S(r)=conj(R(r))*R(r+dr); and

calculating the strain velocity according to an equation SV(r)=c/(4π

drTfo) phase (S(r)).

27. The method according to claim 22, wherein the step of estimating the strain velocity comprises:

estimating a complex pulse-to-pulse correlation for a number of range positions along the ultrasonic beam, based on the echo signals;

calculating a strain correlation function from at least two range positions separated by a given radial distance; and

calculating the strain velocity based on the phase of the strain correlation function.

28. The method according to claim 27 wherein the strain correlation function is given by multiplying the conjugate of the complex pulse-to-pulse correlation for a first range position by the complex pulse-to-pulse correlation for a second range position where said second range position is located the given radial distance from said first range position.

29. The method according to claim 28 wherein the strain velocity is given by dividing a numerator defined as the product of the phase angle of the strain correlation function and the speed of sound by a denominator defined as the product of 4, π

, the given radial distance, the ultrasound frequency and the time between consecutive pulses of said multiple of pulses.

30. The method according to claim 27 further comprising the step of temporally averaging the strain correlation function before calculating the strain velocity.

31. The method according to claim 27 further comprising the step of spatially averaging the strain correlation function before calculating the strain velocity.

32. The method according to claim 22 wherein the step of displaying the estimated strain velocity comprises combining a B-mode tissue image and the strain velocities at corresponding spatial coordinates in a color encoding.

33. The method according to claim 22 wherein the step of displaying the estimated strain velocity comprises combining the M-mode tissue image and the strain velocities at corresponding spatial coordinates in a color encoding.

34. The method according to claims 22 wherein the step of displaying the estimated strain velocities includes displaying the estimated strain velocity for a complete two-dimensional region.

35. The method according to claim 22 wherein the estimated strain velocities are estimated for a three-dimensional region and displayed with a three-dimensional visualization technique.

36. The method according to claim 22 wherein said spatial region is limited to a single point and the strain velocities are displayed as a spectral or strain velocity curve versus time.

37. The method according to claim 22 wherein a said acquiring echo signals is performed with a low pulse repetition frequency.

38. The method according to claim 37 wherein said low pulse repetition frequency is combined with a beam interleaving technique for increasing the frame-rate by acquiring information from different beam directions in the time interval between two consecutive pulses in a given direction.

39. The method according to claim 38 wherein multiple line acquisition techniques are applied to further increase the resulting frame-rate.

40. The method according to claim 22 further comprising the steps of:

calculating a strain reliability index as a function of signal power for each sample point inside the said spatial region; and

modifying the display of the strain velocity according to the strain reliability index.

41. The method according to claim 40 wherein the strain reliability index given by dividing the absolute value of the product of the conjugate of the complex pulse-to-pulse correlation for a first range position and the complex pulse-to-pulse correlation for a second range position by the signal power at one of said first or second range positions.

42. The method according to claim 22 wherein the strain velocity images are generated simultaneously with the application of an external means to produce a pressure gradient in the spatial region.

43. The method according to claim 22 wherein the strain velocity images are generated while simultaneously applying changing pressure against the imaged spatial region.

44. The method according to claim 22 wherein the echo signals are acquired using intravascular ultrasound techniques to assess regional changes in the vessel wall properties.