Ultrasound waveform tomography with spatial and edge regularization
View Patent ↗Synthetic-aperture ultrasound tomography systems and methods using scanning arrays and algorithms configured to simultaneously acquire ultrasound transmission and reflection data, and process the data for improved ultrasound tomography imaging, wherein the tomography imaging comprises total-variation regularization, or a modified total variation regularization, particularly with edge-guided or spatially variant regularization.
1. A synthetic aperture ultrasound tomography imaging method for imaging a tissue medium utilizing a synthetic-aperture ultrasound tomography system comprising two or more ultrasound transducer arrays, the two or more ultrasound transducer arrays being spaced apart from each other, the two or more ultrasound transducer arrays comprising a plurality of transducers, the method comprising:
exciting a first transducer of the plurality of transducers to generate an ultrasound field within the tissue medium;
receiving a transmission signal and a reflection signal at a second transducer of the one or more ultrasound transducer arrays; and
generating an ultrasound waveform tomography image reconstruction, the generating the ultrasound waveform tomography image reconstruction comprising:
applying a regularization function;
applying a plurality of spatially-variant parameters to two or more spatial regions within the ultrasound waveform tomography image reconstruction; and
performing spatially-variant regularization as the regularization function to generate sound-speed reconstructions of the acquired reflection and transmission signals,
wherein generating an ultrasound waveform tomography image reconstruction is a function of computing an acoustic wave property of the reflection and transmission signals by calculating a minimum mean square difference between observed and synthetic waveforms,
wherein the synthetic waveform is generated from simulated ultrasound transmission and reflection data,
wherein the observed waveform is calculated using the reflection and transmission signals, and
wherein the image reconstruction is a function of:
E
(
m
)
=
min
m
{
d
-
f
(
m
)
2
2
+
λ
i
R
(
m
i
)
}
where ∥d−ƒ(m)∥ 2 2 comprises a misfit function, d comprises data relating to the acquired reflection signal and transmission signal, λ i is a positive regularization parameter, and R(m i ) is a spatially-variant regularization term.
2. A method as recited in claim 1 , wherein the spatially-variant parameters vary based on a size of the one or more spatial regions.
3. A method as recited in claim 1 , wherein a gradient of the misfit function is obtained using an adjoint state method.
4. A method as recited in claim 1 , wherein said spatially-variant regularization comprises a modified total-variation regularization comprising two minimization equations:
m
(
k
)
=
arg
min
m
d
-
f
(
m
)
2
2
+
∑
i
λ
1
,
i
m
i
-
u
i
(
k
-
1
)
2
2
,
u
(
k
)
=
arg
min
u
m
(
k
)
-
u
2
2
+
∑
i
λ
2
,
i
∇
u
i
1
,
and
where ∥d−f(m)∥ 2 2 comprises a data misfit function, d comprises data relating to the acquired reflection signal and transmission signals, where λ 1 and λ 2 are both positive regularization parameters, and u i is an auxiliary variable.
5. A method as recited in claim 4 , wherein said step of applying a plurality of spatially-variant parameters comprises:
obtaining approximate locations of the two or more spatial regions; and
applying the plurality of spatially-variant parameters to each of the located two or more spatial regions.
6. A method as recited in claim 5 , further comprising using ray tomography migration results to obtain the approximate locations of the two or more spatial regions.
7. A synthetic aperture ultrasound tomography imaging system for imaging a tissue medium, the system comprising:
two or more ultrasound transducer arrays, the two or more ultrasound transducer arrays being spaced apart from each other, the two or more ultrasound transducer arrays comprising a plurality of transducers;
a processor; and
programming executable on said processor and configured for:
exciting a first transducer of the plurality of transducers to generate an ultrasound field within the tissue medium;
receiving a transmission signal and a reflection signal at a second transducer of the one or more ultrasound transducer arrays; and
generating an ultrasound waveform tomography image reconstruction, the generating the ultrasound waveform tomography image reconstruction comprising:
applying a regularization function;
applying a plurality of spatially-variant parameters to two or more spatial regions within the ultrasound waveform tomography image reconstruction; and
performing a spatially-variant regularization to generate sound-speed reconstructions of the acquired reflection and transmission signals,
wherein generating an ultrasound waveform tomography image reconstruction is a function of computing an acoustic wave property of the reflection and transmission signals by calculating a minimum mean square difference between observed and synthetic waveforms,
wherein the synthetic waveform is generated from simulated ultrasound transmission and reflection data,
wherein the observed waveform is calculated using the reflection and transmission signals, and
wherein the image reconstruction is a function of:
E
(
m
)
=
min
m
{
d
-
f
(
m
)
2
2
+
λ
i
R
(
m
i
)
}
where ∥d−ƒ(m)∥ 2 2 comprises a misfit function, d comprises data relating to the acquired reflection signal and transmission signal, λ i is a positive regularization parameter, and R(m i ) is a spatially-variant regularization term.
8. A system as recited in claim 7 , wherein the spatially-variant parameters vary based on a size of the one or more spatial regions.
9. A system as recited in claim 7 , wherein said programming is further configured to obtain a gradient of the misfit function using an adjoint state system.
10. A system as recited in claim 9 , wherein said spatially-variant regularization comprises a modified total-variation regularization comprising two minimization equations:
m
(
k
)
=
arg
min
m
d
-
f
(
m
)
2
2
+
∑
i
λ
1
,
i
m
i
-
u
i
(
k
-
1
)
2
2
,
u
(
k
)
=
arg
min
u
m
(
k
)
-
u
2
2
+
∑
i
λ
2
,
i
∇
u
i
1
,
and
where ∥d−f(m)∥ 2 2 comprises a data misfit function, d comprises data relating to the acquired reflection signal and transmission signals, where λ 1,i and λ 2,i are both positive regularization parameters, and u i is an auxiliary variable.
11. A synthetic aperture ultrasound tomography imaging system for imaging a tissue medium, the system comprising:
one or more ultrasound transducer arrays;
said one or more ultrasound transducer arrays comprising a plurality of transducers;
a processor; and
programming executable on said processor and configured for:
exciting a first transducer of the plurality of transducers to generate an ultrasound field within the tissue medium;
receiving a transmission signal and a reflection signal at a second transducer of the plurality of transducers of the one or more ultrasound transducer arrays; and
generating an ultrasound waveform tomography image reconstruction, the generating the ultrasound waveform tomography image reconstruction comprising:
applying a regularization function;
applying a plurality of spatially-variant parameters to two or more spatial regions within the ultrasound waveform tomography image reconstruction; and
performing a spatially-variant regularization to generate sound-speed reconstructions of the acquired reflection and transmission signals,
wherein generating an ultrasound waveform tomography image reconstruction is a function of computing an acoustic wave property of the reflection and transmission signals by calculating a minimum mean square difference between observed and synthetic waveforms,
wherein the synthetic waveform is generated from simulated ultrasound transmission and reflection data, and
wherein the observed waveform is calculated using the reflection and transmission signals, and
wherein the image reconstruction is a function of:
E
(
m
)
=
min
m
{
d
-
f
(
m
)
2
2
+
λ
i
R
(
m
i
)
}
where ∥d−ƒ(m)∥ 2 2 comprises a misfit function, d comprises data relating to the acquired reflection signal and transmission signal, λ i is a positive regularization parameter, and R(m i ) is a spatially-variant regularization term.
12. A system as recited in claim 11 :
wherein the plurality of transducers are configured such that a first set of two or more transducers are positioned at an opposing spaced-apart orientation from a second set of two or more transducers such that the first set of two or more transducers face the second set of two or more transducers; and
wherein the first and second sets of two or more transducers are positioned at spaced-apart locations so as to allow for the tissue medium to be positioned in between the first and second sets of two or more transducers.