Filling the deadtime gap in zero echo time MRI: principles compared
Romain Nicolas Froidevaux1, Markus Weiger1, David Otto Brunner1, Bertram Jakob Wilm1, Benjamin Emanuel Dietrich1, and Klaas Paul Pruessmann1

1Institute for Biomedical Engineering, University and ETH Zurich, Zürich, Switzerland

Synopsis

Often used for the imaging of compounds with short relaxation times, 3D center-out sequences performing broadband excitation on a ramped-up strong gradient provide very short acquisition delays and high bandwidth. In turn, data located around the k-space center is missed and needs to be recovered. Different approaches, like algebraic ZTE, WASPI or PETRA are used today. An important but still open question is how the relation of the deadtime and the T2s involved affect image quality. We show that critical situations occur for similar values of deadtime and T2, frequently encountered when imaging ultra-short-T2s in solid tissues or hardware parts.

Introduction

Often used for the imaging of compounds with short transverse relaxation times (<1 ms), 3D center-out sequences performing broadband excitation on a ramped-up strong gradient provide very short acquisition delays and high bandwidth [1] (Fig. 1). Because of finite transmit-receive deadtimes ($$$\Delta$$$), data located around the k-space center is missed and needs to be recovered. To do so, different approaches, often referred to as “zero echo time” techniques, like algebraic ZTE [2], WASPI [3,4] or PETRA [5] are used today. Typically, deadtimes of 5-60 us are realized to target tissues with T2s of 150-500 us with overall promising results. However, an important but still open question is how the relation of the deadtime and the T2s involved affect image quality. Therefore, we investigate the performance of state-of-the-art zero-TE techniques depending on the deadtime duration. We focus on the case where components are present with ultra-short-T2s of tens of us coming from hardware parts or solid tissues, showing that critical situations can occur for similar values of deadtime and T2.

Methods

Sequences: The three techniques use the same basic pulse sequence (Fig. 1) and differ only in the way they recover the data lost during the deadtime (Fig 2). In ZTE, this issue is addressed by radial oversampling and algebraic reconstruction [2] up to gap sizes of 3-4 Nyquist dwells (limited by ill-conditioning) [6], optionally adding k=0 [7]. In PETRA, the lacking centre is measured with single-point imaging (SPI). In WASPI, a second ZTE acquisition is performed at lower gradient strength with optional T2 correction [3] or linear merging in an overlap region [4]. Note that strictly speaking, any method relying on the direct measurement of k0 will have TE=$$$\Delta$$$. Thus only ZTE and WASPI (without k0) are “true” zero-echo-time sequences.

After recovering the missing data, 3D images were reconstructed iteratively with the conjugate gradient method [8] by solving:

$$(E^{H}WDE)v=E^{H}WDm$$

with E=encoding matrix, v=voxel intensities, m=measured data, D=density compensation matrix, W=weighting matrix used for merging in WASPI.

Experiments: First, a stack of rubber erasers (T2* ~=324 us) was imaged with a hard polymer-based surface coil (T2* ~=10 us). Parameters: bandwidth=500 kHz, isotropic resolution=3.75 mm, max. gradient=24.5 mT/m, minimum $$$\Delta$$$=6 us.

Second, a proton-free PTFE-based surface coil was used to image a piece of bovine tibia. Long-lived signals were removed by releasing the bone marrow, leaving mainly bone matrix (T2* ~=10 us) and bound water (T2* ~=160 us). Parameters: bandwidth=250 kHz, isotropic resolution=1.6 mm, max. gradient=29.4 mT/m, minimum $$$\Delta$$$=6 us.

The effect of deadtime on image quality was then investigated by removing data up to 40 us.

Results

1D PSF simulations were performed to investigate the sequences’ performance for ultra-short-T2 signal (Fig. 3). Pure ZTE (without k0) is the only method able to reconstruct the expected ideal Lorentzian lineshape. All other methods diverge from it when T2 approaches $$$\Delta$$$ . PETRA and ZTE with k0 exhibit a narrower center lobe and two side lobes of low amplitude. In WASPI, the main lobe amplitude is lowered while side lobes of increased intensity appear, with variations depending on overlap and use of k0.

Fig 4: The performance observed for the ultra-short-T2 signal of the coil confirms the findings of the PSF simulations. Pure ZTE works well at $$$\Delta$$$=6 us while at $$$\Delta$$$=8 us aliasing of out-of-band signal in the FOV center can be removed by adding k0. PETRA exhibits only minor ringing artifacts at smaller $$$\Delta$$$ while coil signal is suppressed at $$$\Delta$$$=40 us. In WASPI, strong ringing artifacts are observed, at $$$\Delta$$$=40 us even from the longer T2 signal of rubber.

Fig 5: Apart from blurring due to apodization, the ultra-short-T2 components of the bone (mainly collagen) are well reconstructed by ZTE and PETRA at $$$\Delta$$$=12 us. In WASPI, this signal is lost associated with data overlap. At $$$\Delta$$$=40 us, image intensity becomes dominated by bound water which is depicted well only with PETRA.

Conclusion

Compounds with short transverse relaxation times (hundreds of us) can be imaged with high fidelity with all investigated “zero-TE” techniques. However, for deadtimes comparable to the T2s involved, careful consideration is required. ZTE provides maximum short-T2 sensitivity with almost ideal PSF at moderate deadtimes, but may suffer from artifacts due to out-of-band signal at larger deadtimes [6]. In this case, PETRA offers favorable performance, and also enables selecting T2 sensitivity by variation of the deadtime duration. In contrast, WASPI can lead to oscillating PSFs and thus image artifacts, whereupon T2 correction [3] is limited to comparable T2 values.

Acknowledgements

No acknowledgement found.

References

[1] Hafner, S. "Fast imaging in liquids and solids with the back-projection low angle shot (BLAST) technique." Magnetic resonance imaging 12.7 (1994): 1047-1051.

[2] Weiger, M., and K. P. Pruessmann. "MRI with zero echo time." eMagRes (2012).

[3] Wu, Yaotang, et al. "Water-and fat-suppressed proton projection MRI (WASPI) of rat femur bone." Magnetic Resonance in Medicine 57.3 (2007): 554-567.

[4] Gibiino, Fabio, et al. "Free-breathing, zero-TE MR lung imaging." Magnetic Resonance Materials in Physics, Biology and Medicine 28.3 (2014): 207-215.

[5] Grodzki, David M., Peter M. Jakob, and Bjoern Heismann. "Ultrashort echo time imaging using pointwise encoding time reduction with radial acquisition (PETRA)." Magnetic Resonance in Medicine 67.2 (2012): 510-518.

[6] Weiger, Markus, et al. "Exploring the bandwidth limits of ZTE imaging: Spatial response, out-of-band signals, and noise propagation." Magnetic Resonance in Medicine (2014).

[7] Kuethe, Dean O., Natalie L. Adolphi, and Eiichi Fukushima. "Short data-acquisition times improve projection images of lung tissue." Magnetic Resonance in Medicine 57.6 (2007): 1058-1064.

[8] Pruessmann, Klaas P., et al. "Advances in sensitivity encoding with arbitrary k-space trajectories." Magnetic Resonance in Medicine 46.4 (2001): 638-651.

Figures

Main pulse sequence. After ramping up the projection gradient, a short pulse is applied and signal is measured as soon as possible (black dots = acquired, white dots = missed). Such 1D spokes are collected in a center-out radial fashion in order to fill the 3D k-space volume of interest.

k-Space acquisition timing. Main data (black dots) is acquired at maximum gradient strength (Gmax), missing the first samples because of the deadtime gap. Other markers show the additional data points required for each method. The slope is inversely proportional to the gradient strength and normalized by Gmax.

1D PSF simulations. Assumed object: delta function at pixel 0, matrix size: 128, T2: 5 Nyquist dwells (k*FOV), deadtime: 3 Nyquist dwells. Left: signal amplitude vs acquired k-space points. Black circles: data on Nyquist grid. Colored markers: additional data acquired in each method. Right: reconstructed delta function (point spread function).

Phantom imaging. Depiction of the same slice (perpendicular to the coil plane) of a 3D volume. Image magnitude is normalized to maximum intensity. A logarithmic gray scale is used to improve contrast at low intensities.

Imaging of bone. Depiction of the same two orthogonal slices reconstructed with different techniques and deadtimes. In WASPI, the datasets overlap over 20 us. At a deadtime of 12 us, ultra-short T2 components (bone matrix) are well visible in PETRA and ZTE.



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