The present invention concerns a method for characterizing a punctual location inside a medium with ultrasound, and for determining an image of a region inside said medium, in using ultrasound waves propagation inside said medium. The present invention also concerns the apparatus implementing said methods.
Ultrasound imaging of a medium is a very common technique for imaging a medium, and in particular for imaging a human body. Lots of ultrasound imaging techniques are known.
For example, the U.S. Pat. No. 6,095,980 discloses an imaging system and method for detecting linear and nonlinear scatterers inside the medium, called “pulse inversion Doppler”. This system measures the ultrasound response under multiple excitation levels, so as to improve sensitivity in the detection of non-linear responses. Additionally, this method often uses the injection of contrast agents inside the medium for again improving the sensitivity of ultrasound imaging.
Unfortunately, with such an imaging system, the user still gets a weak contrast and in spite of the use of contrast agents, some tissue signals (i.e. without contrast agents) remain in the image and then the bloods vessels contours are insufficiently defined in the image.
One object of the present invention is firstly to provide a method for determining a physical characteristic on a punctual location inside the medium.
To this effect, the method comprises the following steps:
Thanks to these features, the physical characteristic on the punctual location of the medium is determined with an improved accuracy.
In various embodiments of the method, one and/or other of the following features may optionally be incorporated.
In an aspect of the method, the step c) of calculating a phase difference comprises the following sub-steps:
In an aspect of the method, the step c) of calculating a phase difference comprises the following sub-steps:
In an aspect of the method, the physical characteristic is proportional to the phase difference.
In an aspect of the method, the succession of steps a) to c) are repeated several times for providing a plurality of phase differences, and wherein:
In an aspect of the method, the emitted sequence comprises a number greater than two emitted pulses, the received sequence comprises at least the same number of received pulses than the number of emitted pulses, and wherein:
In an aspect of the method:
In an aspect of the method, the second limit is higher than the first limit.
In an aspect of the method, the first limit is 0.3 radian, and preferably 0.1 radian.
In an aspect of the method, the second limit is 0.3 radian, and preferably 0.5 radian.
Another object of the invention is secondly to provide an apparatus for determining physical characteristic on a punctual location inside a medium. The apparatus comprises:
The processing unit implements the method for determining a physical characteristic on a punctual location above specified.
Another object of the invention is thirdly to provide a method for determining an image of a region inside a medium, the image being composed of a plurality of pixels, and wherein the method comprises for a plurality of locations inside the region.
The method comprises the following steps:
Thanks to these features, the image of the medium is determined with an improved sensitivity and accuracy. The image contrast of the image is improved compared to known ultrasound techniques.
In various embodiments of the method, one and/or other of the following features may optionally be incorporated.
In an aspect of the method, the step c) of calculating a phase difference comprises the following sub-steps:
In an aspect of the method, the step c) of calculating a phase difference comprises the following sub-steps:
In an aspect of the method: during step d), the pixel is proportional to the phase difference.
In an aspect of the method, the succession of steps a) to c) are repeated several times for providing a plurality of phase differences for each location inside the region, and:
In an aspect of the method, the emitted sequence comprises a number greater than two emitted pulses, the received sequence comprises at least the same number of received pulses, and:
In an aspect of the method, each pixel in the image corresponds to a location inside the region.
In an aspect of the method:
In an aspect of the method:
In an aspect of the method, the second limit is higher than the first limit.
In an aspect of the method, the first limit is 0.3 radian, and preferably 0.1 radian.
In an aspect of the method, the second limit is 0.3 radian, and preferably 0.5 radian.
In an aspect of the method, two emitted pulses of the emitted sequence have an amplitude ratio of greater or equal to two.
In an aspect of the method, before step a), a contrast agent is introduced into the medium.
In an aspect of the method, the contrast agents comprises microbubbles.
Another object of the invention is fourthly to provide an apparatus for determining an image of a region inside a medium, said apparatus comprising
The processing unit implements the method for determining an image of a region inside a medium above specified.
Other features and advantages of the invention will be apparent from the following detailed description of its embodiments given by way of non-limiting example, with reference to the accompanying drawings.
In the drawings:
In the figures, the same references denote identical or similar elements.
The medium 11 is for instance a living body and in particular human or animal bodies, or can be any other biological or physic-chemical medium (e.g. in vitro medium). The volume of medium comprises variations in its physical properties. For example, the medium may comprise tissues and blood vessels, each one having various physical properties. For example, the tissue may comprise an area suffering from an illness (e.g. cancerous cells), or any other singular area, having various physical properties in comparison to other area of the medium. Some portions of the medium 11 may include some added contrast agent (e.g. micro bubbles) for improving the contrast of physical properties of these portions. When insonified by ultrasound waves, such contrast agents generate non-linear echoes. Therefore, a well-known use of such contrast agents is the injection of a fluid containing a predetermined quantity of contrast agents inside blood vessels. Then, the physical characteristic of such blood vessels can be more easily detected in comparison to a physical characteristic of a tissue that does not comprise the contrast agent, said contrast agent only flowing in the vessels.
The physical characteristics, that can be detected by the method that senses the medium via ultrasound waves, are mechanical properties of the medium, like stiffness, or else. The method distinguishes values and/or variations of said physical properties. For example, the method may detect mechanical interface between two materials in the medium: e.g. it can detect bubbles shell.
The apparatus 10 may include:
In a variant, a single electronic device could fulfil all the functionalities of the electronic unit 13 and of the processing unit 14. The processing unit 14 may be a computer.
The probe 12 can comprise a curved transducer so as to perform an ultrasound focusing to a predetermined position in front of the probe. The probe 12 can comprise a linear array of transducers, few tens of transducers (for instance 100 to 300) juxtaposed along an axis X so as to perform ultrasound focusing into a bi-dimensional (2D) plane. The probe 12 can comprise a bi-dimensional array so as to perform ultrasound focusing into a tri-dimensional (3D) volume.
A first configuration of the method represented on
The above processing unit 14 controls the electronic unit 13 and the probe 12 for:
The ultrasound waves We, Wr toward and from the location can be a focused wave (beam) or a non focused beam.
A well-known beamforming method can be used:
The at least two different amplitudes of emitted pulses can be produced by varying the transmit voltage or preferably by varying the aperture size (i.e. by varying the number of transducers elements contributing to emit the emitted ultrasound wave). The aperture may also be divided into two or more groups of elements.
In the first configuration of the method, the processing unit 14 further operates the following steps:
Each emitted sequence ES of ultrasound waves We comprises several emitted pulses (at least two). Also, each received sequence RS of ultrasound waves Wr comprises several received pulses.
Each pulse in these sequences may be a more complex signal, for example being a combination of a plurality of sine signals. However, we will consider in the present description that the pulses are only sine signals for the sake of simplicity.
A sine signal or pulse comprises a phase Φ and amplitude X. The phase of such sine signal is the angle at its origin in a time scale representation of the signal. The amplitude is the maximum value of the sine signal over time. Such phase can be determined or estimated by various methods: time shift relative to a reference clock signal, inter-correlation with said reference clock signal, calculation of Arctan(Im/Re) in case the signal is represented by complex values (for example, in case of IQ demodulation), or any other known method.
A phase difference between two sine signals or pulses can be also determined by any one of the above-cited methods, and e.g. simply by differences of predetermined phases of each signal.
Then, the phase difference of step c) can be determined by various methods. Only two of them will be presented in the following description.
According to a first method for determining phases, the phase difference δ between two sequences of two sine signals (as used in step c) of the method) can be defined as the difference of:
The defined phase difference δ is δ=δr−δe.
This phase difference can also be calculated by the time shifts of the signals relative one to another: in this case, the phase difference δ is given by:
δ=2.π.f.(Tr−Te)=2.π.f.(Δt)
where
If the time difference (Tr−Te) is greater than 2.π, it should be reduced to the value modulo 2.π.
According to a second method for determining phases, the two received sine signals may be scaled (weighted) so as to have the same amplitude, and then the resulted weighted signals are subtracted one to the other, i.e.:
δr=Φ(a.Sr2−b.Sr1).
Eventually, the first weight coefficient a may be unity, and only δ is determined so as b.Sr1 has the same amplitude as Sr2.
The inventors have found that the above defined phase difference δ of the received signals (received pulses) relative to the emitted signals (emitted pulses) surprisingly depends on the nature of the medium 11 at the focused point P0. Therefore, the physical characteristic of said medium at the location P0 can be determined.
This effect is due to the difference of amplitude of first and second emitted pulses (Se1, Se2), and to the non-linear behaviour of the medium at point P0. This non-linear behaviour of medium is usually analysed in ultrasound imaging via amplitude harmonics analysis (second harmonic, third harmonic . . . ). This is an amplitude and frequency analysis. In present case, the method uses phase difference analysis. Of course, such phase technique can also be combined with harmonic analysis.
Advantageously, the first and second emitted signals ultrasound pulses (signals Se1, Se2) have a ratio of amplitudes greater or equal to two.
Additionally, a sum of the emitted signals in the emitted sequence (ES) is null. Therefore, a simple method of summing the received signals cancels the first harmonic and this sum can be directly used to estimate the phase and/or harmonic behaviour of the material in the medium 11.
More specifically, the emit sequence ES can comprise only two emit pulses of inverted amplitudes, denoted as a +1, −1 sequence, such as in pulse inversion method as described in U.S. Pat. No. 6,095,980. In that case, the weighting factors a, b in the above second method for calculating the phase difference should also be preferably (+1, +1) so that the sum of the signals is null.
More specifically, the emit sequence ES can comprise three emit pulses of +½, −1, +½ amplitudes, such as in power modulated method as described in U.S. Pat. No. 6,638,228. In that case, the above second method for calculating the phase difference should for example use the following formula:
δr=Φ(Sr1+Sr2+Sr3),
i.e. having the following weighting factors (+1, +1, +1).
Such method may be generalized as follow:
In the method, the amplitudes and the weighting factors are chosen so that e1.r1+e2.r2+e3.r3=0. Then, the combined signal CS should be a null signal if the medium at the location P0 behaves as a linear material.
The inventors have discovered that looking at the phase of the resulting signal (the phase of the combined signal), that is the phase difference δ, is representative of the nature of the medium 11 at the location P0.
Thus, the phase difference δ measured at the given punctual location P0 can be used for determining the type or nature of medium at this punctual location P0.
In a first variant, the physical characteristic PC of the material of medium 11 at or near the punctual location P0 is proportional to the value of the previously determined phase difference δ.
In a second variant, the succession of steps a) to c) are repeated several times, for example M iterations, for providing a plurality of phase differences δ, said plurality having then M values of phase difference: δ1, . . . , δM.
Then, during the step d), the physical characteristic PC is determined by:
The M iterations can be proceeded with the same first and second amplitudes for the first and second emitted pulses (signals Se1, Se2) in the emitted sequence ES.
The M iterations can be proceeded with predetermined varying values of amplitudes for the first and second emitted signal Se1, Se2 in the emitted sequence ES, so as to test various variations of amplitudes and to test various possibilities of non-linearities in the medium. Therefore, thanks to the above variant, a non-linear behaviour of the medium 11 at location P0 can be more easily distinguished from a linear behaviour.
In a third variant, the emitted sequence ES may comprise more than two signals of various amplitudes: it may comprise a number N of pulses (signals): three, four, five or even more. The received sequence RS then comprises at least the same number N of received pulses (signals), and preferably the same number N of received pulses (signals). Each received signal Sr1 of index i must be correctly associated to the corresponding emitted signal Sei of same index i.
Then, during the step c), a plurality of couples or pairs (index i and index j, i=/j, and i, j<=N) of emitted-received signals can be extracted, each couple of emitted signals Sei, Sej having a different amplitude. For each couple, a phase difference δij is determined as explained above, so as to provide a plurality of phase differences.
Then, during the step d), the physical characteristic PC is determined by (ij being all the predetermined couples)):
The above mathematical functions are usually well known defined or equivalent of them.
The N signals may all have different amplitudes or some of them can be identical. Then, various variations of amplitudes can be tested or experienced during the process. Therefore, a non-linear behaviour inside the medium 11 at location P0 can be more easily distinguished from a linear behaviour.
This distribution curve 20 comprises a mean value (abscissa of the maximum of the curve) near zero radian, and a small standard deviation or variance (the width S of the curve).
The absolute value of the mean value A is for instance lower than 0.1 radian, and the standard deviation is higher than 0.1 radian.
Such measurement corresponds to a location in the medium 11 having a substantially linear behaviour: It may be a location or region of tissue, or not including any contrast agent (e.g. micro-bubbles).
This distribution curve 20 comprises a mean value A (abscissa of the maximum of the curve) that is non zero radian, and a standard deviation S or variance much higher (the width of the curve) than in the case of the
The absolute value of the mean value A is for instance higher than 0.3 radian, and the standard deviation is higher than 0.3 radian. Both values can therefore be used to distinguish a first type of material (physical characteristic) inside the medium from a second type of material (physical characteristic).
Therefore, all the above variants of the method can then include the following steps for classifying or segmenting the nature or type of physical characteristic at the location P0 inside the medium 11:
Eventually, the second limit L2 is higher than the first limit L1, and if the physical characteristic is inside the interval of these two limits (L1, L2), the type of material is not determined.
Eventually, the second limit L2 is equal to the first limit L1.
Thanks, to these classification steps, the type of physical characteristic is determined.
For example, the first limit L1 is 0.3 radian, or 0.1 radian, so as to determine if the physical characteristic is low (linear behaviour).
For example, the first second limit L2 is 0.3 radian, or 0.5 radian, so as to determine if the physical characteristic is high (non-linear behaviour).
These limits can be adapted to each application of the method: distinguishing a region including a contrast agent from another region not including the contrast agent, or distinguishing tissue to blood vessels, or distinguishing illness cells to healthy cells.
A second configuration of the method is, illustrated on
The method uses identical or similar elements of the above disclosed apparatus 10 of
The image produced by the method is composed of a plurality of pixels (for example, a number K of pixels), each pixel corresponding to a different location (Pk) inside the region R, k being an index to identify each pixel in the image or each location in the region R. Eventually, the image may be composed of only one pixel. However, the image may preferably comprise more than one ten thousand pixels (100×100 image).
The second configuration of the method (imaging method) mainly differs from the previous first configuration of the method by the scanning of a plurality of locations inside a region R so as to generate an image of said region.
At each location Pk inside the region R, the processing unit 14 controls the electronic unit 13 and the probe 12 for:
Similarly as for the first configuration, the ultrasound waves We, Wr can be focused or non focused waves, according to known techniques.
The emitted and received signals (representing the pulses) are also similar or identical to those as represented on
In the method of the second configuration of the invention, the processing unit 14 operates then the following steps for each location Pk inside the region R:
The definitions given for a signal phase and a phase difference δ during the first configuration of the method are still valid. There are various methods for determining a phase or phase difference.
The phase difference δ at step c) can be determined by various methods, and least one of the following two methods.
According to a first method, the phase difference δ is a difference of two phase differences: the one concerning the received signals (δr: measured) and the one concerning the emitted signals (δe: known). The phase difference is then for example calculated by the following formula (explained in the first configuration of the method):
δ=δr−δe=Φ(Sr2)−Φ(Sr1)−δe.
Then, the method comprises the following sub-steps:
According to a second method, the phase difference δ is a phase of a combined signal, said combined signal being a weighted sum of the received signals.
For example, in case of two pulses (signals), the phase difference can be calculated by the following formula:
δr=Φ(a.Sr2−b.Sr1),
The weighting factors a, b being determined to have cancelation of the weighted signals (a.Sr2, b.Sr1) if the medium has a linear behaviour.
More generally, the phase difference δ is:
δ=Φ(CS)=Φ(r1.Sr1+r2.Sr2+r3.Sr3).
The amplitudes e1, e2, e3 for the emit pulses and the weighting factors r1, r2, r3 for the received signals are chosen so that e1.r1+e2.r2+e3.r3=0.
Then, the method comprises the following sub-steps:
In a first variant, the pixel value corresponding to each location Pk inside the region R is directly the value of the previously determined phase difference δ.
In a second variant, the succession of steps a) to c) are repeated several times, for example M iterations, for providing a plurality of phase differences δ, said plurality having then M values of phase difference: δ1, . . . , δM.
Then, during the step d), the pixel value is determined by a mean value or standard deviation value or a variance value of the plurality of phase differences.
The M iterations can be proceeded with the same first and second amplitudes for the first and second emitted pulses (signals Se1, Se2) in the emitted sequence ES.
The M iterations can be proceeded with predetermined varying values of amplitudes for the first and second emitted pulses (signals Se1, Se2) in the emitted sequence ES, so as to test various variations of amplitudes and to test various possibilities of non-linearities of the medium at each location Pk. Therefore, a non-linear behaviour of each location Pk in the region can be more easily distinguished from a linear behaviour.
In a third variant, the emitted sequence ES may comprise more than two pulses of various amplitudes: it may comprise a number N of pulses: three, four, five or more pulses. The received sequence RS then comprises at least the same number N of received pulses, and preferably the same number N of received pulses. Each received signal Sr1 of index i (corresponding to a pulse i) must be correctly associated to the corresponding emitted pulse (i.e. emitted signal Sei) of same index i.
Then, during the step c), a plurality of couples or pairs (index i and index j, i=/j, and i, j<=N) of emitted-received signals can be extracted, each couple of emitted signals Sei, Sej having a different amplitude. For each couple, a phase difference δij is determined as explained above, so as to provide a plurality of phase differences.
Then, during the step d), the pixel value is determined by a mean value e or standard deviation value or a variance value of the plurality of phase differences
The N emitted pulses (signals) may all have different amplitudes or some of them can be identical. Then, various variations of amplitudes can be tested during the process. Therefore, a non-linear behaviour of the medium 11 at each location Pk inside the region R can be more easily distinguished from a linear behaviour.
In a fourth variant, the region R is subdivided in a plurality of sub-regions, corresponding to sub-area in the image.
Between step c) and step d) of the method, a sub-region phase difference (δsr) is calculated, said sub-region phase difference (δsr) being a mean value of the plurality of phase differences (δ) of the locations belonging to said sub-region calculated during step c).
Then, a pixel in the image corresponds to a sub-region.
Additionally, any of the above variant can be combined to another one to have a method that provides a more sensitive image of the medium: non-linear behaviour are more easily distinguished from a linear behaviour.
Similarly to the first configuration of the invention, all the variants of the method can then include following steps for classifying or segmenting the image:
Eventually, the second limit L2 is higher than the first limit L1, and if the pixel value is inside the interval of these two limits (L1, L2), the type the pixel is not determined.
Eventually, the second limit L2 is equal to the first limit L1.
Thanks, to these classification steps, the image can be segmented into sub-areas having different properties.
For example, the first limit L1 is 0.3 radian, or 0.1 radian, so as to determine a sub-area of a first type (linear behaviour).
For example, the second limit L2 is 0.3 radian, or 0.5 radian, so as to determine a sub-area of a second type (non-linear behaviour).
These limits can be adapted to each application of the method: distinguishing tissue to blood vessels, or distinguishing illness cells to healthy cells.
For example, the
Then, in all the present methods can take advantage in using contrast agents inside the medium 11.
For example, contrast agents penetrate blood vessels but not into the tissues. Then, the contrast between the different area of the medium 11 with and without contrast agents will be increased, as illustrated on
Ultrasound contrast agents could be small gas-filled micro-bubbles contrast agents. Small gas-filled micro-bubbles could be for example seeded into a blood stream or from encapsulated micro-bubbles. Contrast agents could consist of gas-filled microscopically small bubbles encapsulated by an elastic shell. There are a variety of micro-bubbles contrast agents. Micro-bubbles could be heavy gas/air mixture (PFC, SF8), encapsulated (shell—lipid, albumin, polymer). Micro-bubbles differ in their shell makeup, gas core makeup, and whether or not they are targeted. Micro-bubble size is fairly uniform. They lie within a range of 1 to 10 micro-meters in diameter. That makes them smaller than red blood cells, which allows them to flow easily through the circulation as well as the microcirculation. Micro-bubbles could be administrated by intravenous injection to the systemic circulation. They could be removed dissolved in circulation, filtered by liver and cleared 15 minutes.
Micro-bubbles have a high degree of ability to reflect the ultrasound waves. When micro-bubbles are caught in an ultrasonic frequency field, they compress, oscillate, and reflect a characteristic echo, generating a non-linear ultrasound backscatter. The ultrasound pulses are thus reflected in a different way by the gas in the medium containing micro-bubbles and by the soft tissue surroundings.
This is a Continuation application under 35 U.S.C. § 120 claiming priority to U.S. patent application Ser. No. 15/111,214 filed on Jul. 13, 2016; which is a 35 U.S.C. § 371 National Phase filing of PCT/EP2015/051369 filed on Jan. 23, 2015; which claims priority under 35 U.S.C. § 119 to U.S. Provisional Application No. 61/930,805 filed on Jan. 23, 2014.
Number | Name | Date | Kind |
---|---|---|---|
5577505 | Brock-Fisher et al. | Nov 1996 | A |
5706819 | Hwang et al. | Jan 1998 | A |
6095980 | Burns et al. | Aug 2000 | A |
6398732 | Brock-Fisher et al. | Jun 2002 | B1 |
6497665 | Hunt et al. | Dec 2002 | B1 |
6506158 | Kawagishi et al. | Jan 2003 | B2 |
6638228 | Brock-Fisher et al. | Oct 2003 | B1 |
11058401 | Bruce | Jul 2021 | B2 |
Number | Date | Country |
---|---|---|
2001286472 | Oct 2001 | JP |
5041258 | Oct 2012 | JP |
2008016992 | Feb 2008 | WO |
2010121265 | Oct 2010 | WO |
Entry |
---|
International Search Report related to Application No. PCT/EP2015/051369; reported dated May 11, 2015. |
Number | Date | Country | |
---|---|---|---|
20210267576 A1 | Sep 2021 | US |
Number | Date | Country | |
---|---|---|---|
61930805 | Jan 2014 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 15111214 | US | |
Child | 17326000 | US |