Erpeng Dai1, Xiaodong Ma1, Zhe Zhang1, Chun Yuan1,2, and Hua Guo1
1Center for Biomedical Imaging Research, Department of Biomedical Engineering, School of Medicine, Tsinghua University, Beijing, China, People's Republic of, 2Vascular Imaging Laboratory, Department of Radiology, University of Washington, Seattle, WA, United States
Synopsis
One of the challenges for interleaved EPI (iEPI)
DWI is the phase inconsistency among different shots. Several methods, performed
either in the image or k-space domain, have been proposed to solve this problem
with extra acquired navigator data. However, the navigator is usually acquired
with a lower bandwidth in the phase encoding direction than the image echo,
which can cause different distortion levels. In this study, the effects of such
distortion for the image or k-space based reconstruction are investigated. It
has been shown that the k-space based method is more tolerant to the navigator
distortion.Purpose
Compared with single-shot EPI, interleaved EPI (iEPI)
DWI can provide higher spatial resolution and less geometric distortions.
However, physiological motion can cause phase inconsistency among different
shots, thus inducing ghost artifacts in DWI images. Several methods, performed either
in the image
1 or k-space
2-4 domain, have been proposed
to solve this problem with extra acquired navigator data. However, the
navigator is usually acquired with a lower bandwidth in the phase encoding (PE)
direction than the image echo, which can cause different distortion levels. In
this study, the effects of such distortion on the performance of the image or
k-space based reconstruction are investigated.
Methods
Sequence The sequence diagram for navigated spin echo EPI DWI
acquisitions is shown in Fig. 1(a). The echo spacing of the navigator (Tnav)
is first set to the same as the image echo. Then it is reduced to 1/Rnav
of the echo spacing of the image echo, where Rnav is the reduction
factor. Noticeably, to ensure enough sampling points for a navigator, when Rnav>1,
the number of PE lines for the navigator is increased. The
corresponding sampling trajectories for different Rnav (Rnav
= 1, 2, 3) are shown in Fig. 1(b). Meanwhile, the sampling bandwidth in the
readout (RO) direction and the slope time of the RO gradient for the navigator
are unchanged.
Reconstruction The iEPI DWI reconstruction can be
performed either in the image or k-space domain. In the image domain, the low
resolution navigator image is first reconstructed, and then the shot-to-shot
phase variations are calculated and corrected 1. While in the k-space
domain, data from different shots can be taken as encoded by both phase variations
and coil sensitivities, and the missing data of each shot are recovered through
convolution, like GRAPPA 2-5. In this study, image reconstruction
using image-space sampling function (IRIS) 1 is used as the image
based reconstruction, while synergistic image reconstruction with phase
variation and sensitivity (SYMPHONY, previously named SEPARATE) 2, 4
is used as the k-space based reconstruction. It should be noted that for IRIS
used here, no image registration between navigators and DWI images is
performed. After the reconstruction, the SNR is calculated using the pseudo-multiple
replica method 6 with 128 repetitions.
Experiments All scans were
performed on a Philips 3.0T Achieva TX MRI scanner (Philips Healthcare, Best,
The Netherlands) using a 32-channel head coil. All human studies were performed
under IRB approval from our institution. The imaging parameters were as
follows: FOV = 240 × 240 mm2, acquisition
voxel size = 1 × 1 × 4 mm3, 4 slices
with gap = 16mm, 8 shots with 29 echoes per shot, no partial Fourier, TE = 77ms,
TR = 3s, scan time = 1.75min, diffusion preparation was applied in 3 orthogonal
directions with b = 800 s/mm2, NSA = 1. For different Rnav,
the acquisition matrix sizes of the navigator are summarized
in Table 1.
Results and
Discussion
Fig.
2(a) shows the navigator images for different Rnav. For a
better illustration, all navigator images are zero-padded to 240 × 240. The
slice consisting of frontal lobe of the brain (indicated by green line in Fig. 2(b)) is chosen for display, where distortions are severe at the air-tissue interfaces (indicated
by red arrow head). A T2 weighted TSE image with the same resolution 1 × 1 × 4
mm3 is shown in Fig. 2(c) as reference. As shown, with Rnav
increasing, the navigator distortion is largely reduced.
Fig.
3 shows a comparison between DWI images reconstructed with the image (a) or
k-space (b) based methods for different navigator distortion levels. With a severely
distorted navigator, the DWI image reconstructed from the image based method is
largely corrupted (yellow arrow heads in Fig. 3(a)). When the navigator
distortion is minor (Rnav=3), the image based method can achieve
comparable performance as the k-space based method and even show a slightly
higher SNR, as shown in Fig. 3(c) and (d). In general, the k-space based method
works well for different navigator distortion levels.
Conclusion
In interleaved EPI DWI, the k-space based
reconstruction is more tolerant to navigator distortions, compared with the
image based method. This indicates that the k-space method can be a better
choice when the image registration between image and navigator is not
easily implemented.
Acknowledgements
This
work was supported by National Natural Science Foundation of China (61271132,
61571258) and Beijing Natural Science Foundation (7142091).References
[1] Jeong HK, et al. Magn Reson
Med 2013.
[2] Ma X et al. ISMRM
2015; p2799.
[3] Liu W et al. Magn Reson Med.
2015;00.
[4] Guo H at al. ISMRM Workshop
on SMS Imaging 2015.
[5] Griswold M a et al.
Magn Reson Med. 2002;47(6).
[6] Robson PM, et al. Magn Reson Med 2008;60.