Multispectral Body Diffusion-Weighted Imaging
Valentina Taviani1, Shreyas S. Vasanawala1, and Brian A. Hargreaves1

1Stanford University, Stanford, CA, United States

Synopsis

A diffusion-weighted (DW) imaging method was developed to mitigate off-resonance-induced distortion and signal loss, which are problematic for body applications. A 2D RF pulse is used in place of the conventional spectral-spatial excitation used for DW spin echo echo-planar imaging. In the presence of off-resonance, a narrow band of frequencies is excited due to the different bandwidths between excitation and refocusing pulses. By progressively shifting the center frequency, the whole range of off-resonance can be excited and a composite image, corrected for off-resonance-induced distortion, can be reconstructed by estimating the field map from the spectral information.

Purpose

B0 inhomogeneity is a major challenge encountered in body diffusion-weighted (DW) echo-planar imaging (EPI), leading to severe signal loss and distortion. The numerous air-tissue interfaces in the torso, the difficulty to accurately shim large volumes and the presence of fat often result in suboptimal image quality. The use of 2D RF pulses in conjunction with parallel imaging (PI) for coronal whole-body DW-EPI has been described1. While the 2D spatial selectivity can both improve PI performance by phase encoding along the direction with the largest variation in coil sensitivities and reduce distortion by limiting the FOV to the anatomy of interest, the low bandwidth of the 2D RF pulse makes this technique extremely sensitive to off-resonance, resulting in severe signal loss2. We present a method that can uniformly excite the whole imaging volume while correcting residual off-resonance-induced distortion in clinically acceptable scan times. Phantom and volunteer experiments were performed to demonstrate feasibility.

Theory

2D RF pulses can be used to control the extent of the excited volume along both the slice-select and phase-encoding direction. Here, the blipped direction during excitation corresponds to the slice-select direction during imaging, so that only the main excitation lobe is refocused by the 180° refocusing pulse3. Due to the low slice-select bandwidth of the 2D RF pulse (220Hz) and the relatively high bandwidth of the refocusing pulse (1.2kHz), off-resonance spins are shifted outside of the volume refocused by the 180º, so that the combination of the 2 pulses effectively produces a 2D spatial/1D spectral excitation of finite bandwidth. By progressively shifting the center frequency, the whole range of off-resonance can be excited, as demonstrated in Figure 1, where a 30Hz/cm shim gradient was used to simulate off-resonance. Root-sum-of-squares (RSOS) combination of the spectral components restores signal over the whole FOV but does not correct off-resonance-induced distortion. However, knowledge of the spectral components on a pixel-by-pixel basis allows estimation of the off-resonance field (i.e. field map), which, in turn, can be used to correct off-resonance induced distortion (d=ΔB0/prBW; d=shift along phase-encode direction in pixels, prBW=pseudo receiver bandwidth per pixel). A center of mass (COM) approach was used to estimate the field map from the frequency “bin” images4 (cfr. Figure 1).

Materials and Methods

Phantom experiments were performed to assess accuracy of the estimated field map and effectiveness of the distortion correction algorithm. A separate IDEAL acquisition was performed to obtain a reference field map (6 TEs; TE1/DTE=0.9/0.7ms; flip angle=3º; BW=111kHz; FOV=44cm; phase FOV=50%; 7mm slice thickness; matrix size=160×160; 16 slices). The known shape of the phantoms was used to evaluate the performance of the distortion correction algorithm. Healthy volunteers were scanned using the upper 20 elements of a 32-channel receive-only phased-array coil (NeoCoil, Pewaukee, WI) and an 8-channel receive-only neurovascular coil (GE Healthcare, Waukesha, WI) to demonstrate feasibility. All imaging was performed at 3T (GE MR750, Waukesha, WI).

Results and Discussion

Figure 2 shows RSOS and multispectral (MS) reconstructions using IDEAL and estimated field maps. Both field maps were able to successfully correct the distortion induced by the shim gradient. Note that distortion in the original acquired images is reflected in the COM field map, while the IDEAL-generated field is undistorted. In vivo experiments showed similar off-resonance distributions across different subjects for a specific anatomy. In all cases, 4-5 frequency bins were sufficient to cover the range of off-resonance, which could be achieved with scans under five minutes. Figure 3 shows on-resonance and composite MS images of 2 axial abdominal slices. Note that in this case the frequency spread was minimal in the first slice but large enough to cause significant signal loss, successfully recovered by the MS reconstruction, in a different slice (arrow). Figure 4 and 5 show coronal and axial images of the brachial plexus, where conventional on-resonance imaging was non-diagnostic. Overall, good suppression of on-resonance fat was achieved in the MS images (cfr. RSOS, Figure 4), although some off-resonance fat (arrow) remains due to errors in the estimated field map, that could be corrected by acquiring more frequency bins while keeping scan time constant by reducing the number of signal averages per bin.

Conclusion

We developed a DW-EPI method to reduce off-resonance-induced signal loss and distortion, which are especially problematic for body applications. The main limitation consists in the longer duration of the acquisition (proportional to the number of acquired bins). While a few bins are usually sufficient to cover the whole off-resonance range, multi-slice capabilities can be added to the proposed method to extend coverage/reduce scan time.

Acknowledgements

Dr. Maggie Fung for useful discussion; GE Healthcare; NIH P41-EB015891-18.

References

1. Obara P, Taviani V, Loening AM, Iagaru A, Hargreaves BA, Vasanawala SS. Comparison between whole-body coronal and axial DWI performed during PET-MR. Proc. Intl. Soc. Mag. Reson. Med. 2015; 23:3996.

2. Fung M, Estkowski L, Xu D, Banerjee S, Worters P, Wu G, Skinner J, Bayram E. Coronal whole body diffusion imaging with 2D spatially selective excitation (FOCUS). Proc. Intl. Soc. Mag. Reson. Med. 2014; 22:2238.

3. Saritas EU, Cunningham CH, Lee JH, Han ET, Nishimura DG. DWI of the spinal cord with reduced FOV single-shot EPI. Magn Reson Med 2008; 60:468-473.

4. Koch KM, Brau AC, Chen W, Gold GE, Hargreaves BA, Koff M, McKinnon GC, Potter HG, King KF. Imaging Near Metal with a MAVRIC-SEMAC Hybrid. Magn Reson Med 2011; 65:71-82.

Figures

In the presence of off-resonance, a narrow range of frequencies contributes signal due to the different bandwidths between excitation and refocusing pulses. By progressively shifting the center frequency, the whole off-resonance range can be excited. The acquired spectral information can be used to correct off-resonance-induced distortion in the composite image.

MS reconstructions using both IDEAL and estimated field maps successfully corrected off-resonance-induced distortion (cfr. RSOS). COM and IDEAL field maps were in good agreement. Acquisition parameters were: TE=59ms; TR = 2s; FOV=44cm; phase FOV=50%; 7mm slice thickness; matrix size=80×160; 16 slices; 11 frequency bins, 200Hz apart.

DW (b=600s/mm2) on-resonance and MS abdominal images (TE=49ms; TR = 2s; FOV=36cm; phase FOV=50%; 7mm slice thickness; matrix size=80×160; 16 axial slices; 5 bins 200Hz apart; PI factor=2; b=50 (8NEX) and b=600s/mm2 (16NEX)). Note recovered signal in proximity to the spine with MS reconstruction.

Frequency “bin” images, RSOS and MS reconstructions in a healthy volunteer. TE=49ms; TR = 2s; FOV=48cm; phase FOV=50%; 8mm slice thickness; matrix size=128×128; 16 coronal slices; 5 bins 200Hz apart; PI factor=2; b=50 (8NEX) and b=600s/mm2 (16NEX). Note residual unsuppressed off-resonance fat due to inaccurate field map estimation.

On resonance and MS images acquired in a healthy volunteer. TE=49ms; TR = 2s; FOV=48; phase FOV=50%; 8mm slice thickness; matrix size=128×128; 16 axial slices; 5 bins 200Hz apart; PI factor=2; b=50 (8NEX) and b=600s/mm2 (16NEX). Note significant signal loss in the “0Hz” images recovered with the MS approach.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0622