Method to obtain 3D images of a flowing region beneath an object using speckle reflections

Abstract

A method for imaging a flowing media within static regions includes obtaining a plurality of signals using the speckle properties of the flowing media. The plurality of signals are compared to one another such as by subtraction. The static regions are removed from the plurality of signals by the comparison. The remaining signals are combined (such as by summing) to produce an image of the flowing media.

Claims

1. A method for imaging a flowing media behind static regions including the steps of: a) obtaining a plurality of acoustic echo signals containing information about speckle properties of the flowing media; b) comparing each of the plurality of acoustic echo signals to each of the others of the plurality of acoustic echo signals; and c) removing the static regions from the plurality of acoustic echo signals during said step b).

2. The method of claim 1 wherein said step b) includes subtracting the plurality of acoustic echo signals from one another.

3. The method of claim 1 further including the step of enhancing the speckle properties of the plurality of acoustic echo signals by filters or other denoising operations.

4. The method of claim 1 further including the step of enhancing the speckle properties of the plurality of acoustic echo signals by convolution or correlation operations.

5. The method of claim 1 further including the step of enhancing the speckle properties of the plurality of acoustic echo signals by transmission to phase space.

6. The method of claim 5 wherein said step b) includes comparing the differences between the plurality of acoustic echo signals.

7. The method of claim 5 wherein said step b) includes comparing the plurality of acoustic echo signals via subtraction.

8. The method of claim 1 wherein the comparing in step b) is subtracting.

9. The method of claim 1 further including the step of: d) combining the results from each of the comparisons of step b).

10. The method of claim 1 wherein said step b) includes comparing each of the plurality of acoustic echo signals with each of the others of the plurality of acoustic echo signals via principle component analysis or cross comparison subtraction methods.

11. The method of claim 1 further including the step of combining the results from each of the comparisons of step b).

12. The method of claim 1 wherein the flowing media is blood within a skull and within soft tissue, wherein the static regions are the soft tissue within the skull.

13. The method of claim 1 wherein said step a) is performed at a location fixed relative to the static regions such that the plurality of acoustic echo signals pass through the same static regions.

14. The method of claim 1 wherein step a) is performed by a transmitter that remains fixed with respect to the static regions and generates a plurality of acoustic signals that result in the plurality of acoustic echo signals.

15. The method of claim 1 wherein step a) is performed by a transmitter that remains fixed with respect to the static regions and generates a plurality of acoustic signals that result in the plurality of acoustic echo signals.

16. A method for imaging a flowing media behind static regions including the steps of: a) obtaining a plurality of acoustic echo signals containing information about speckle properties of the flowing media; b) comparing the plurality of acoustic echo signals to one another; and c) removing the static regions from the plurality of acoustic echo signals during said step b), wherein the flowing media is blood within a skull and within soft tissue, wherein the static regions are portions of the skull.

17. The method of claim 16 wherein the portions of the skull are of varying thickness.

18. The method of claim 16 wherein said step b) includes subtracting the plurality of acoustic echo signals from one another.

19. The method of claim 16 further including the step of enhancing the speckle properties of the plurality of acoustic echo signals by filters or other denoising operations.

20. The method of claim 16 further including the step of enhancing the speckle properties of the plurality of acoustic echo signals by convolution or correlation operations.

21. The method of claim 16 further including the step of enhancing the speckle properties of the plurality of acoustic echo signals by transmission to phase space.

22. The method of claim 21 wherein said step b) includes comparing the differences between the plurality of acoustic echo signals.

23. The method of claim 16 further including the step of: d) combining the results from each of the comparisons of step b).

24. A method for imaging a flowing media behind static regions including the steps of: a) obtaining a plurality of acoustic echo signals containing information about the speckle properties of the flowing media; b) enhancing the speckle properties of the plurality of acoustic echo signals by transmission to phase space; c) comparing the plurality of acoustic echo signals to one another; and d) removing the static regions from the plurality of acoustic echo signals during said step c), wherein the transmission to phase space generates a plurality of phase signals and wherein said step b) includes comparing dissimilarities between the plurality of phase signals via principle component analysis or cross comparison subtraction.

25. The method of claim 24 wherein step a) is performed by a transmitter that remains fixed with respect to the static regions and generates a plurality of acoustic signals that result in the plurality of acoustic echo signals.

26. The method of claim 25 wherein the flowing media is blood within a skull and within soft tissue, wherein the static regions are the soft tissue within the skull.

27. A method for imaging blood flowing within a skull including the steps of: a) transmitting a plurality of signals toward a portion of the skull and the blood flowing behind the portion of the skull; b) receiving a plurality of echo signals from the portion of the skull and from the blood containing information about speckle properties of the flowing media; c) comparing the plurality of echo signals to one another; and d) based upon said step c), removing portions of the plurality of echo signals caused by the portion of the skull from the plurality of echo signals received during said step b).

28. The method of claim 27 wherein the portion of the skull includes areas of varying thickness.

29. The method of claim 27 wherein step a) is performed by a transmitter that remains fixed with respect to the skull and generates a plurality of acoustic signals that result in the plurality of acoustic echo signals.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 schematically shows a flow imaging system monitoring blood flow in a head.

(2) FIG. 2 shows one embodiment of the flow imaging technique, showing speckles being enhanced, then combined over the special domain to produce an image showing blood flow regions.

(3) FIG. 3 is an example 2D cross section of a flow imaging output from the system of FIG. 1.

(4) FIG. 4 is an example 2D cross section of a flow imaging output from the system of FIG. 1.

(5) FIG. 5 shows multiple cross sections of a flow imaging output from the system of FIG. 1.

(6) FIG. 6 shows an example 3D output of flow imaging system of FIG. 1.

(7) FIG. 7 shows an example 3D output of flow imaging system of FIG. 1.

DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT

(8) A flow imaging system 10 according to one embodiment is shown schematically in FIG. 1. The system 10 includes an imaging device 12, a controller 14 and an image-processing computer 16. The system 10 can be used, for example, as shown in FIG. 1 on a human head 18. The head 18 includes a skull 20 outward of soft tissue 22. A blood vessel 24 is within the soft tissue 22 spaced away from the skull 20. Blood 26 flows through the blood vessel 24.

(9) Generally, images from the imaging device 12 (optionally, as controlled by controller 14) are received and stored by the flow imaging computer 16. Generally, the algorithms described herein are performed by the flow imaging computer 16 on the images stored thereon, as received from the imaging device 12. The flow imaging computer 16 includes a processor, storage (such as memory, mass storage, hard drive, optical drives, or other magnetic, electronic or optical storage), optional graphics processor and user interfaces, such as a display.

(10) Many alternatives for hardware are possible in implementing the method of the present invention. The present techniques discussed here are not restricted to any specific device, and can be implemented in a new device designed for this particular methodology, or an existing device modified to obtain the required signal or data.

(11) 1. Signal Acquisition

(12) The imaging device 12 can be an existing device to obtain the required signals or data. This imaging device 12 (or system) is responsible for signal acquisition.

(13) The scanning device 12 is any device capable of generating, receiving, and storing data created from radiation of the physical factor obtaining an A-Scan (temporal vs. Amplitude data, depth vs. Amplitude, temporal vs. intensity, etc depending on scanning system). This imaging device 12 is able to obtain this data over a 1D, or 2D area in space, in addition to A-Scan data above.

(14) The imaging device 12 can be made up of separate transmitters (source) and receivers, or a single transmitter-receiver. The source, or transmitter, may be any of the following: a single element transmitter, a one-dimensional array, or two-dimensional array. This generates the field to create the radiation for the physical factor. The receiver is preferred to be a one dimensional or two-dimensional array. A motorized single element receiver is possible. This device receives the signal of the radiation from the physical factor.

(15) These signals acquired by the imaging device 12 can be processed in parallel by the methodology, or signals obtained and then processed. It is possible to enhance the results of the methodology by a second data acquisition by the imaging device 12.

(16) The imaging device 12 obtains data in the form of multiple A-Scans (amplitude vs. time), with each set of A-Scans being obtained at a specific spatial region.

(17) Multiple sets of A-Scans can also be taken at different spatial regions form a 2D, or 3D image, this can be done by receiving data at different elements of the receiver, steering of the beam, or controlled physical motion of the receiver. The methodology depends on the scanning system in question.

(18) 2. Methodology & Data Processing

(19) General Overview

(20) The invention uses the set of A-Scans (from scanning system above) to create an image, or sub-set of a larger images indicating the region of blood flow

(21) This is done by comparing the information from the individual acoustic waveforms acquired during the measurement (or comparing processed waveforms). Two examples are provided below:

(22) Comparing pairs of waveforms:

(23) i.e. compare signal 1 with signal 2, signal 2 with 3, and so on for a set number of signals;

(24) i.e. compare signal 1 with signal 2, signal 3 with 4, and so on for a set number of signals;

(25) This comparison can also be done across comparing multiple waveforms:

(26) i.e. signal 1 is compared with signals 2, 3, 4 (and so on), then scan 2 is compared with signals 3, 4 (and so on), the process repeating for a set number of signals.

(27) Alternative patterns are also possible

(28) Techniques such as principle component analysis may be used to find these differences.

(29) The results of these comparisons are combined via non-linear imaging, and produce the pixels to the resulting image, or sub-set of a larger image. The image can then be further enhanced and displayed in two or three dimensions by secondary means. Qualitative or quantitative details on the image can be obtained

(30) General Formulation

(31) The invention's overall algorithm combines multiple difference-based measurements in order to derive information about the blood flow. The algorithm involves the following steps:

(32) Transmitting and receiving multiple signals through the medium containing the static and dynamic components. The position of the transmitter/receiver equipment of the imaging device 12 should remain fixed with respect to the static component of the medium.

(33) Each A-Scan in the received data can be filtered to remove the background noise, and enhance the signal via Low-Pass, Band-Pass, or High-Pass Filters, Convolutions, or Cross Correlations with the known, or estimated signal from the scanning device any other variety of other denoising, and enhancing operations from those knowledgeable in the subject.

(34) For each A-Scan the phase may be measured, as described below (see Measuring Phase subsection).

(35) The differences are extracted using cross comparison of the waveforms as explained in the General overview subsection above. Information about the differences between the data may also be found by techniques such as principle component analysis.

(36) Combine the above differences from all comparisons, corresponding to a specific instance or an interval of time, via summation/integration of their absolute values (possibly of a complex number), squared values; or using any other method that highlights the cumulative changes and suppresses the static component. Information from two separate difference calculations may be combined to produce a new result.

(37) Use obtained number to assign a value to a pixel in the 2D or 3D image. The coordinates of said pixel correspond to spatial-temporal location of the medium where the wave incurred the phase changes measured above. The process above is repeated for different positions of the receiver, or different elements in the receiver, obtaining data at a new spatial position in the object. The formulation can be adjusted to the various receiver types (single element, 1D, and 2D discussed earlier) by those knowledgeable in the subject. The formulation can be adjusted for the different comparison methods discussed above by those knowledgeable in the subject. The formulation can be adjusted for continuous phase information by those knowledgeable in the subject.

(38) Measuring the Phase

(39) Phase information can be obtained via such methods as Fourier Transform (Short Fourier Transforms, Short Shifting Fourier Transforms, and other Fourier methods), and Hilbert Transform. The method need not be restricted to only these transforms. The phase of the signal can be measured at a particular instance of time, such as with a Hilbert transform; or over a finite interval of time, such as with short-Fourier transform. By comparing multiple phase measurements of the signal at various instances in time, or intervals of time, the method obtains information can be obtained about blood flow at various depths through the tissue (corresponding to the time interval). By collecting and processing a number of signals propagating at different spatial positions, it is possible to measure the phase and derive the blood flow information over the entire volume of said medium. From this a variety of 2D and 3D images of the blood flow in the medium can be obtained (3D structure, cross sections, box sections, etc.)

(40) 3. Further Processing

(41) The receiver's position can be altered incrementally to increase the size of the image created. The receiver's position can be altered to enhance a previously imaged area, obtain a finer detailed image, or image the area again. Additional image processing techniques can be used to scale, enhance, smooth, or display the image. The image created can be further refined to deal with irregularities introduced by inhomogenoties in the skull, or tissue.

(42) The present method uses a statistical, or multiple information approach to enhance information lost due to the presence of the skull, or other object.

(43) Multiple signals are obtained (and processed to enhance the signal)

(44) Each signal is compared mathematically to all other signals, this allows the process to turns N signals into roughly N.sup.2/2 measurements

(45) Multiple comparisons are then combined.

(46) The invention is not specifically probing the nonlinear properties of the medium, and as such the invention does not require physical change in signals transmitted into the medium. The invention assumes flowing liquids are flowing and dynamic (non-static), and that this flow of scatterers in the liquid produces different received signals due to this flow producing (a changing speckle response). It is this speckle property of the flowing liquid that is analyzed not nonlinear response.

(47) Possible ways to enhance these speckles from the medium to enhance image contrast have been disclosed above.

(48) The removal of phase inversion (pulse inversion) means the present method requires no specific pulse paring for the pulse comparison. The invention uses N like signals, and allows for the processed signal 1 to be compared to any processed signal from signal 2 through signal N by mathematical means (i.e. subtraction). This multi-comparison option further increases the contrast and detection of the flow.

(49) Comparison of multiple signals by this invention removes static objects, and truly images the flowing region. Pulse inversion techniques assume a flow is essentially static, and there is no differentiation between static and dynamic regions. This method provides options to use, and enhance, A-Scan data (A-Scan data can then used to make 2D B-Scans, 3D C-Scan images, or the 2D and 3D flow images of this invention)

(50) In accordance with the provisions of the patent statutes and jurisprudence, exemplary configurations described above are considered to represent a preferred embodiment of the invention. However, it should be noted that the invention can be practiced otherwise than as specifically illustrated and described without departing from its spirit or scope.