Detection of motion in vibro-acoustography
View Patent ↗Harmonic motion is produced in a subject using vibro-acoustography. An ultrasonic imaging system repetitively interrogates the subject and the Doppler shift in the reflected echo signals is analyzed to measure the phase and amplitude of harmonic motion produced in the subject at different prescribed frequencies. Shear wave propagation through the subject is determined from this information and mechanical properties related to “stiffness” of the subject are determined. A Kalman filter is employed in the phase and amplitude measurement to extract the harmonic motion information from background noise.
1. A vibro-acoustic system which comprises:
a) a force generator for imparting harmonic motion to a subject at a prescribed frequency, the harmonic motion propagating over a given distance;
b) an ultrasonic system including:
an ultrasonic transducer;
a transmitter connected to the ultrasonic transducer and being operable therewith to produce an ultrasonic beam which repeatedly applies ultrasonic energy to each of a plurality of spaced apart prescribed points on the subject through which the imparted harmonic motion propagates;
a receiver connected to the ultrasonic transducer and being operable therewith to receive from each of the spaced apart prescribed points a corresponding plurality of ultrasonic echoes; and
c) a processor connected to receive the ultrasonic echoes and including a Kalman filter which for each of the spaced apart prescribed points extracts from the corresponding plurality of ultrasonic echoes an amplitude and a phase of the harmonic motion at that prescribed point, and wherein the processor is operable to calculate a mechanical property using the extracted information.
2. The system as recited in claim 1 in which the receiver operates to produce quadrature I and Q signals from each plurality of ultrasound echoes and the Kalman filter receives the quadrature I and Q signals.
3. The system as recited in claim 1 in which the subject is tissue and the mechanical property is complex stiffness.
4. A method for measuring a mechanical property of a subject, the steps comprising:
a) applying a force to the subject which imparts harmonic motion thereto at a prescribed frequency, the harmonic motion propagating over a given distance;
b) sensing the harmonic motion at a plurality of spaced apart prescribed points on the subject, through which the harmonic motion propagates, with an ultrasonic system that applies a pulsed ultrasonic beam to each prescribed point and receives a corresponding plurality of ultrasonic echo signals from each prescribed point;
c) processing the corresponding plurality of ultrasonic echo signals to detect the harmonic motion at each of the prescribed points and to extract a phase of the harmonic motion at each of the prescribed points;
d) calculating the mechanical property using the extracted phases.
5. The method as recited in claim 4 , wherein the subject is tissue.
6. The method as recited in claim 4 , wherein the mechanical property is one of shear elasticity and viscosity.
7. The method as recited in claim 4 in which step c) includes:
i) demodulating the received ultrasonic echo signals to produce quadrature I and Q signals;
ii) calculating the arctangent of the ratio of said I and Q signals; and
iii) filtering the signal which results from step ii) with a Kalman filter to extract the phase of the harmonic motion therein at the prescribed frequency.
8. A method for measuring a mechanical property of a subject, the steps comprising:
a) applying a force to the subject which imparts harmonic motion thereto at a prescribed frequency;
b) sensing the harmonic motion at a plurality of prescribed points on the subject with an ultrasonic system that repeatedly applies a pulsed ultrasonic beam to each prescribed point and receives a corresponding plurality of ultrasonic echo signals from each prescribed point;
c) processing the corresponding plurality of ultrasonic echo signals to detect the harmonic motion at each of the prescribed points and to extract a phase of the harmonic motion at each of the prescribed points;
d) calculating a first shear wave speed using the extracted phases at the prescribed points;
e) repeating steps a) through c) at a different prescribed frequency and calculating a second shear wave speed using the extracted phases at the prescribed points at the different prescribed frequency; and
f) calculating the mechanical property using the first wave speed and the second wave speed.
9. The method as recited in claim 8 in which the subject is tissue.
10. The method as recited in claim 8 wherein step e) is repeated a plurality of times to calculate a plurality of corresponding shear wave speeds and the mechanical property is calculated using all the calculated shear wave speeds.
11. The method as recited in claim 8 in which the mechanical property is one of shear elasticity and viscosity.
12. The method as recited in claim 8 in which both shear elasticity and viscosity are calculated using the first wave speed and the second wave speed.
13. The method as recited in claim 8 in which steps b) and c) include:
i) demodulating received ultrasonic echo signals to produce quadrature I and Q signals;
ii) calculating the arctangent of the ratio of said I and Q signals; and
iii) filtering the signal which results from step ii) using a Kalman filter to extract a phase of the harmonic motion therein at the prescribed frequency.
14. The method as recited in claim 13 , further including modeling a desired harmonic signal representative of the harmonic motion by a differential equation and recursively estimating the amplitude and the phase of the harmonic motion at each of the prescribed points using the Kalman filter.
15. The method as recited in claim 8 in which step a) includes applying an acoustic force to the subject.