Reverberant shear wave field estimation of body properties
A reverberant shear wave field in an object such as a patient's body or organ causes deformations in one or more selected directions measured with an imaging modality such as ultrasound or MR equipment or other imaging equipment, to estimate displacements in one or more selected directions over time increments and then viscoelastic properties such as stiffness or other parameters of the ROI.
1 . A method of using a reverberant shear wave field in a body to estimate viscoelastic properties of hidden regions of interest in the body, comprising:
actuating a vibration source that produces plural shear waves at selected frequencies that interact with each other and with structures in the body to produce a reverberant shear wave field in the body;
imaging a region of interest in the body in the presence of said reverberant shear wave field with an imaging device and producing an estimate of measured motion of the region of interest in a selected first direction;
receiving in an image processor as inputs the selected frequencies of the shear waves and said estimate of measured motion, and to process the inputs with computer algorithms to provide an estimate of one or more viscoelastic properties of said region of interest in the body; displaying said estimate of one or more viscoelastic properties of the region of interest on a computer display; and
controlling said steps of producing shear waves, imaging, computer processing, and computer display with a controller operatively coupled to said vibration source, said imaging device, said computer display, and said image processor.
2 . The method of claim 1 , in which the vibration source produces shear waves at one or more frequencies in the range of 1600-2400 Hz.
3 . The method of claim 1 , in which the vibration source produces shear waves at one or more frequencies in the range of 1000-4000 Hz.
4 . The method of claim 1 in which the vibration source comprises 3-7 individual vibration sources.
5 . The method of claim 1 in which the vibration source comprises more than 7 individual vibration sources.
6 . The method of claim 1 in which the vibration source is an extended vibration surface.
7 . The method of claim 1 , in which said estimate of measured motion includes a position in said first selected direction of the region of interest at plural respective times and speed of the change in position in said first selected direction, and the image processor is configured to take images of the region of interest and calculate said one or more estimates as a function of auto-correlation of said images.
8 . The method of claim 1 , in which the imaging device is further configured to measure motion in a second selected direction and provide an estimate of the measured motion in the second direction.
9 . The method of claim 1 , in which the image processor is configured to provide a map of said one or more viscoelastic properties of each of plural points within the region of interest.
10 . The method of claim 1 , in which the image processor is configured to provide a two-dimensional map of said one or more viscoelastic properties of each of plural points within the region of interest.
11 . The method of claim 1 , in which the image processor is configured to provide a three-dimensional map of said one or more viscoelastic properties of each of plural points within the region of interest.
12 . The method of claim 1 , in which the step of using an imaging device comprises using an ultrasound scanner to produce said time sequence of ultrasound images of the region of interest and said estimate of motion.
13 . The method of claim 1 , in which the step of using an imaging device comprises using an MRI scanner to produce said time sequence of ultrasound images of the region of interest and said estimate of motion.
14 . The method of claim 1 , in which the step of using an imaging device comprises using an Optical Coherence Tomography (OTC) system configured to provide a time sequence of OTC images of the region of interest and said estimate of motion.
15 . The method of claim 1 , in which the step of using an imaging device comprises using an x-ray imaging system configured to provide a time sequence of x-ray images of the region of interest and said estimate of motion.