System and method for visualization of tissue microvasculature using ultrasound
View Patent ↗Methods and systems for producing a visually-perceived representation of a sub-millimeter-sized blood vessel located at a depth of many centimeters in the biological tissue, in which the background clutter is suppressed (by at least 30 dB using SVT and additional 23 dB using a combination of morphology filtering and vessel enhancement filtering) as compared to an image obtained with the use of a B-mode ultrasound imaging, while at the same time maintaining the morphology of the blood vessel.
1. A method for generating an image depicting tissue microvasculature using an ultrasound system, the steps of the method comprising:
(a) acquiring ultrasound echo data from a region-of-interest in a subject with an ultrasound system, wherein the region-of-interest contains tissue microvasculature;
(b) generating clutter removed data by performing a clutter removal algorithm on the ultrasound echo data, wherein the clutter removal algorithm is implemented with a hardware processor and a memory;
(c) generating spectrum data by performing a temporal Fourier transform on the clutter removed data;
(d) storing negative frequency data in the spectrum data as negative frequency spectrum data;
(e) storing positive frequency data in the spectrum data as positive frequency spectrum data;
(f) reconstructing a first image from the negative frequency spectrum data;
(g) reconstructing a second image from the positive frequency spectrum data; and
(h) generating a third image that depicts tissue microvasculature by combining the first image and the second image in order to reduce artifacts and noise.
2. The method of claim 1 , wherein the third image is generated by computing a difference between the first image and the second image.
3. The method of claim 1 , wherein the clutter removal algorithm implements clutter removal using a frequency selective filtering algorithm.
4. The method of claim 1 , wherein the clutter removal algorithm implements clutter removal using a low rank approximation.
5. The method of claim 1 , wherein the clutter removal algorithm implements clutter removal using a singular value thresholding.
6. The method of claim 5 , wherein the singular value thresholding comprises global singular value thresholding.
7. The method of claim 1 , wherein motion data are estimated from the ultrasound echo data using a non-rigid motion estimation.
8. The method of claim 7 , wherein the ultrasound echo data are processed to compensate for motion using the motion data before performing clutter removal.
9. The method of claim 7 , wherein the clutter removed data are processed to compensate for motion using the motion data before performing the temporal Fourier transform on the clutter removed data.
10. The method of claim 7 , wherein the motion data are estimated based on a global time-delay estimation.
11. The method of claim 7 , wherein the motion data comprise a displacement field.
12. The method of claim 1 , further comprising processing the third image using a morphology-based filter to remove spatially variant background signals having a spatial scale that is different from a spatial scale of the tissue microvasculature, wherein the morphology-based filter comprises a circular morphology structure.
13. The method of claim 1 , further comprising processing the third image using a vessel enhancement filter that equalizes variations of intensity across the third image to compensate for fluctuations caused by relative orientation of a probe of the ultrasound imaging system and a blood vessel in the tissue microvasculature.