Dongyeob Han1, Min-Oh Kim1, Honpyo Lee1, Taehwa Hong1, and Dong-Hyun Kim1
1Yonsei University, Seoul, Korea, Republic of
Synopsis
A simultaneous, free-breathing water/fat separation
and T1, T2 quantification method was proposed. Dual TR (in-phase and out-phase
TR) and varying sinusoidal flip angle was used with FISP acquisition. For
motion robustness, random rotating golden angle trajectories were applied. T1,
T2 and Δφfat of fat were pre-determined using
the fat dominant region mask, then water/fat signal combined dictionary was
generated. The results show that the water/fraction maps from the proposed
method were in good agreement with conventional breath-hold results. Furthermore,
measured T1, T2 values were in good agreement with the values from the previous
research. Purpose
In abdomen imaging, several problems such as motion
and B0 field inhomogeneity hinder accurate diagnosis and cause imaging
artifacts. Often, the fat signal is separated based on techniques originating
from the Dixon method1. The gating methods2 or fast
imaging methods3 with single breath-hold are used to overcome the
motion artifact. Herein, to acquire water/fat fraction and T1, T2
quantification simultaneously without these limitations, a dual TR FISP with
random rotating golden angle radial acquisition while free breathing was
proposed.
Methods
Figure
1 shows the proposed algorithm. In the data acquisition step, dual TR pattern
with alternating in-phase TR (4.6 ms) and out-phase TR (6.9 ms) and varying
sinusoidal flip angle pattern were applied with FISP (unbalanced SSFP)
acquisition scheme which is robust to B0 field inhomogeneity4.
Random rotating golden angle radial trajectories5 were also applied
to grant incoherency of undersampling artifacts and motion robustness. A total
of 600 highly undersampled (16 spokes for each images, reduction factor 8) 2D
images were acquired in 3T Siemens scanner (Tim Trio) with 2x2mm2
resolution, TE=TR/2, 256x256 mm2 FOV, 5 mm slice thickness and 1560
Hz/Px bandwidth for a scan time 54 sec/slices.
To generate water/fat signal combined dictionary, fat's T1, T2 and Δφfat (phase difference between in-phase TR and
out-phase TR due to chemical shift of fat) were pre-determined with an assumption
that the fat's T1, T2 and Δφfat values are uniform
in the subject. Fat dominant region mask was generated from the difference
between the average phase images acquired from in-phase TR and out-phase TR as
shown in Fig. 1. An average fat signal was acquired using the fat dominant
region mask, then T1, T2 and Δφfat of the
subject's fat were pre-determined. Finally, water/fat signal combined
dictionary was generated based on a water signal dictionary applying complex
sum of the fat signal evolution which was pre-determined. The water/fat
dictionary has water fraction : [0:0.05:1], water T1 : [500:25:1980,
2000:200:3000] ms and water T2 : [5:5:145, 150:50:250] ms.
To validate the proposed method, single breath-hold
(exhale) multi-echo GRE was acquired with TE : 2.2, 3.0, 3.8 ms, TR : 6.1 ms
and FA : 18˚ with same FOV and
resolution from the above. For water/fat separation, a graph cut based method6
was applied. Furthermore, an additional breath-hold experiment was performed
using the proposed dual TR FISP pulse sequence to validate the motion
robustness of random rotating golden angle trajectories. The data was collected
during 16 breath-hold (exhale) states (breath-hold duration 600 TRs=3.4 sec)
which were synchronized with the pulse sequence. All experiments were performed
with a healthy volunteer (IRB-approved) in both axial and sagittal slices.
Results
Figure 2 and 3 shows the results of the proposed
method for each axial and sagittal slices. Water/fat fraction maps from the
proposed method are in good agreement with the fraction maps from the
multi-echo gradient echo data in both axial and sagittal slices. Signal
evolutions in point 'x' shows that there are no meaningful difference between
the signal evolution pattern of breath-hold and free breathing. This signal
evolution results and comparison results shown in Figures 2 and 3 between
breath-hold and free-breathing shows that the random rotating golden angle
trajectories can successfully reduce the motion artifacts due to the
respiratory motion. Furthermore, ROI-averaged T1 and T2 values in the liver from the proposed
free-breathing method are T1 : 726 ms and T2 : 52 ms, which are in good
agreement with previous known T1, T2 values in liver3.
Discussion
& Conclusion
The
proposed method based on the random rotating golden angle radial acquisition
and dictionary matching method can successfully suppress motion artifacts arising
from respiration. However, further analysis for non-periodic or hard-breathing respiratory
motion is required, due to all experiments being performed with a healthy
volunteer in this study.
Streaking artifact due to strong flow signal
from the aorta was shown in the results. Additional flow suppress preparation technique
may also required.
We have presented a simultaneous water/fat separation and T1, T2 quantification method using dual TR FISP with random rotating
golden angle radial trajectories while free breathing in the abdomen. The
method would be useful for patients who cannot maintain breath-hold and also
for quantification research in the abdomen.
Acknowledgements
This research was supported by Global PH.D Fellowship Program through the National Research Foundation of Korea(NRF) funded by the Ministry of Education (NRF-2015H1A2A1030698).
References
1.
Dixon WT. Simple proton spectroscopic imaging. Radiology 1984; 153: 189-194.
2.
Ehman RL, McNamara MT, Pallack M, Hricak H, Higgins CB. Magnetic resonance
imaging with respiratory gating: techniques and advantages. American Journal of
Roentgenology 1984;143:1175-1182.
3.
de Bazelaire CM, Duhamel GD, Rofsky NM, Alsop DC. MR imaging relaxation times
of abdominal and pelvic tissues measured in vivo at 3.0 T: preliminary results
1. Radiology 2004;230:652-659.
4.
Jiang Y, Ma D, Seiberlich N, Gulani V, Griswold MA. MR fingerprint-ing using
fast imaging with steady state precession (FISP) with spiralreadout. Magn Reson
Med 2014. doi: 10.1002/mrm.25559.
5.
Winkelmann S, Schaeffter T, Koehler T, Eggers H, Doessel O. An optimal radial
profile order based on the Golden Ratio for time-resolved MRI. Medical Imaging,
IEEE Transactions on 2007;26:68-76.
6.
Hernando D, Kellman P, Haldar J, Liang ZP. Robust water/fat separation in the
presence of large field inhomogeneities using a graph cut algorithm. Magnetic
Resonance in Medicine 2010;63:79-90.