High spatial and temporal resolution dynamic contrast-enhanced magnetic resonance imaging
View Patent ↗A method for high spatial and temporal resolution dynamic contrast enhanced magnetic resonance imaging using a random subsampled Cartesian k-space using a Poisson-disk random pattern acquisition strategy and a compressed sensing reconstruction algorithm incorporating magnitude image subtraction is presented. One reconstruction uses a split-Bregman minimization of the sum of the L1 norm of the pixel-wise magnitude difference between two successive temporal frames, a fidelity term and a total variation (TV) sparsity term.
1. A method of contrast enhanced magnetic resonance imaging, the method comprising:
(a) acquiring a plurality of magnetic resonance imaging data sets of a subject from a scanner by randomly under-sampling Cartesian k-space; and
(b) reconstructing images from the acquired data sets with a split-Bregman minimization of the sum of the L 1 norm of the pixel-wise magnitude difference between two successive temporal frames, a fidelity term and a total variation (TV) sparsity term.
2. A method as recited in claim 1 , further comprising optimizing weighting parameters by retrospectively under-sampling the k-space and minimizing root mean square error between original k-space images and reconstructed images.
3. A method as recited in claim 2 , further comprising normalizing acquired data by maximum image intensity prior to image reconstruction.
4. A method as recited in claim 1 , wherein said random under-sampling of Cartesian k-space uses variable density Poisson-disk random patterns.
5. A method as recited in claim 1 , further comprising:
plotting a sampling trajectory on a ky-kz plane in k-space; and
sampling k-space at the beginning of each temporal frame to avoid motion sensitivity and eddy currents that can affect imaging.
6. A method as recited in claim 1 , wherein said reconstruction further comprises:
minimizing cost function for a volume by dividing a cost function into multiple Bregman distances and updating the distances according to a norm and a threshold.
7. A method as recited in claim 1 , further comprising:
acquiring dynamic contrast enhanced magnetic resonance imaging under-sampled k-space datasets for a plurality of temporal frames; and
reconstructing images from the datasets by solving:
(
I
1
,
I
2
)
=
argmin
{
U
1
F
(
I
1
)
-
K
1
2
2
+
λ
TV
(
I
1
)
+
μ
I
2
e
i
ϕ
1
-
I
1
1
U
2
F
(
I
2
)
-
K
2
2
2
+
λ
TV
(
I
2
)
+
μ
I
2
-
I
1
e
1
i
ϕ
2
1
}
;
wherein fidelity terms ∥U 1 F(I 1 )−K 1 ∥ 2 2 and ∥U 2 F(I 2 )−K 2 ∥ 2 2 compare the Fourier transform (F) of the estimated images (I 1 and I 2 ) to measured samples (K 1 and K 2 ) using sampling masks (U 1 and U 2 );
wherein the total variation terms λTV(I 1 ) and λTV(I 2 ) are the L1-norm of the discrete gradients of I 1 and I 2 ; and
wherein the ∥I 2 |e iϕ 1 −I 1 | 1 and |I 2 −|I 1 |e iϕ 2 | 1 terms are the L1 norm of the pixel-wise magnitude subtracted image: ∥I 2 |−|I 1 ∥ 1 =∥I 2 |e iϕ 1 −I 1 | 1 =|I 2 −|I 1 |e iϕ 2 | 1 .