1304

Iterative Reconstruction of 23Na Multi-Channel Breast Data Using Compressed Sensing Combined with Anatomical 1H Prior Knowledge
Sebastian Lachner1, Olgica Zaric2, Matthias Utzschneider1, Lenka Minarikova2, Stefan Zbyn3, Bernhard Hensel4, Siegfried Trattnig2, Michael Uder1, and Armin M. Nagel1,5

1Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany, 2High Field MR Center, Department of Biomedical Imaging and Image-guided Therapy, Medical University of Vienna, Vienna, Austria, 3Research Unit of Medical Imaging, Physics and Technology, University of Oulu, Oulu, Finland, 4Center for Medical Physics and Engineering, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany, 5Division of Medical Physics in Radiology, German Cancer Research Centre (DKFZ), Heidelberg, Germany

Synopsis

An iterative reconstruction algorithm for sodium magnetic resonance imaging (23Na MRI) with multi-channel receiver coils is implemented and compared to a conventional gridding reconstruction. Based on compressed sensing (CS) it utilizes a total variation (TV(2)), combined with anatomical weighting factors (AnaWeTV(2)) to preserve known tissue boundaries. Simulated and measured 23Na multi-channel data sets of the female breast were reconstructed. The TV(2) and in particular the AnaWeTV(2) lead to an improved image quality, due to effective noise reduction and the highlighting of structure. The presented CS reconstruction is beneficial especially for high undersampling factors.

Introduction

Sodium magnetic resonance imaging (23Na MRI) suffers from an inherently low signal-to-noise ratio, due to low in vivo concentrations, low MR sensitivity and short relaxation times. As a consequence, the spatial resolution is limited and image quality can be affected by partial volume effects or Gibbs ringing artifacts. Advanced image reconstruction methods, such as compressed sensing (CS) [1], can be used to overcome these restrictions and improve image quality. Based on that, we propose an iterative reconstruction algorithm for multi-channel data sets, which utilizes a total variation (TV(2)), adopted by weighting factors (AnaWeTV(2)) [2] to preserve known tissue boundaries. In this work, the AnaWeTV(2) reconstruction is compared with a conventional TV(2) and the gridding reconstruction [3] using simulated and measured 23Na MRI multi-channel data of the female breast.

Methods

Reconstruction algorithm: The image reconstruction is done via a CS scheme [1]. Thereby the objective function

$$f(x)=\frac{1}{2} \left\|\textbf{A}x-y\right\|_2^2+\tau\sum\limits_{\alpha=x,y,z} \left( \lambda \left\|\textbf{W}_\alpha \textbf{D}_\alpha^{(1)}x\right\|_1+(1-\lambda) \left\|\textbf{W}_\alpha \textbf{D}_\alpha^{(2)}x\right\|_1 \right)$$

is minimized through an iterative algorithm that employs the conjugate gradient method. The first term ensures data consistency by calculating the difference between the raw data y and the current image x; the second term represents the second order total variation, which works as a conventional denoising technique, adopted by weighting factors W [2]. They are extracted from a high-resolution 1H image by calculating the first derivation in each direction, and thus reduce the total variation at known tissue boundaries. Via τ parameter, the data consistency and the sparsity term can be weighed against each other; τ is optimized for each reconstructed channel and data set.

Simulations: Dixon 1H measurements (FOV: 246x320x128 mm3, resolution: (1 mm)3, TAQ = 1:46 min) are used to create a fat and water binary mask (BM) (Fig. 1a,b). The masks are added (1$$$\cdot$$$BMfat + 6$$$\cdot$$$BMwater) to generate a high-resolution ground truth with a tissue contrast similar to that in a 23Na image (Fig. 1c). By multiplying the ground truth with calculated sensitivity maps, a corresponding multi-channel ground truth is obtained. To reduce the spatial resolution to (3x3x1 mm3), high frequency k-space data is cut off, followed by a biexponential decay (T2f* = 2 ms, T2s* = 17 ms) and the regridding to a radial trajectory. Finally, synthetic correlated complex Gaussian noise is added.

The simulated data sets are used to test and optimize the algorithm. Results are compared to the existing ground truth (Fig. 1c). For this, the structural similarity (SSIM) [4] inside the object is calculated to assess the image quality; thereby a value of 1 corresponds to complete agreement between the reconstructed image and the ground truth.

Measurements: This study was approved by an institutional review board. In vivo MR measurements of a female healthy subject (age: 27) were carried out at a 7T whole-body system (Magnetom, Siemens Healthcare, Erlangen, Germany), using a dual-tuned proton/ sodium (1H/ 23Na) bilateral phased-array breast coil (Quality Electrodynamics, OH, USA). A density adapted 3D radial projection reconstruction pulse sequence (DA-3DRP) [5] was employed (TR/ TE = 30 ms/ 0.55 ms, FA = 67°, TRO = 9.98 ms, FOV: 320 mm, nominal resolution: (3 mm)3), to carry out four 23Na measurements with different undersampling factors (USFs):

  • 20.000 projections, USF = 1.8
  • 10.000 projections, USF = 3.6
  • 5.000 projections, USF = 7.2
  • 2.500 projections, USF = 14.4

Results

Fig. 2 depicts the results of reconstructed simulated data, combined via a sum-of-squares (SOS) reconstruction.
The use of a simple TV(2) as sparsity transformation leads to better results compared with the conventional gridding reconstruction: the SSIM indicates an improved image quality, especially for high USFs (Fig. 2a,b).
The AnaWeTV(2) compared with the TV(2) as sparsity transformation preserves known tissue boundaries and improves the resolving of fine structures, both indicated by a higher SSIM value (Fig. 2c). This can also be comprehended by calculating the differences between the TV(2) and the AnaWeTV(2) reconstruction, where edges can be tracked systematically (Fig. 2d).
While the SSIM decreases quickly with higher USFs for the gridding reconstruction, it decreases significantly slower in the case of the CS reconstruction.

The reconstruction of measured data is in agreement with the results of simulations (Fig. 3).

Discussion and Conclusion

The presented CS algorithm leads to a significant improvement of image quality compared to the gridding reconstruction, due to the effective reduction of noise and imaging artifacts. Using the AnaWeTV(2), known tissue boundaries are preserved, which leads to a further improvement. Especially for high USFs, the CS reconstruction shows a clear advantage over the gridding reconstruction. Using the presented algorithm the acquisition time can be reduced while image quality is preserved.

Acknowledgements

This work was supported by the Vienna Science and Technology Fund (WWTF, project LS14-096).

References

[1] Lustig M, Donoho D, Pauly JM. Sparse MRI: The application of compressed sensing for rapid MR imaging. Magn Reson Med. 2007;58(6):1182-1195.
[2] Gnahm C, Nagel AM. Anatomically weighted second-order total variation reconstruction of 23 Na MRI using prior information from 1 H MRI. Neuroimage. 2015;105:452-461.
[3] Jackson JI, Meyer CH,et al. Selection of a convolution function for Fourier inversion using gridding (computerised tomography application). IEEE Trans Med Imag. 1991;10(3):473-478.
[4] Wang Z, Bovik AC, et al. Image quality assessment: from error visibility to structural similarity. IEEE Trans Image Process. 2004;13(4):600-612.
[5] Nagel AM, Laun FB, et al. Sodium MRI using a density‐adapted 3D radial acquisition technique. Magn Reson Med. 2009;62(6):1565-1573.

Figures

Fig. 1: a) 1H fat and water Dixon measurements are used to create b) binary masks (BMs). The BMs are added (1$$$\cdot$$$BMfat + 6$$$\cdot$$$BMwater) to generate a c) high-resolution ground truth with a contrast similar to that in a 23Na image. This reference ground truth is multiplied with sensitivity maps; high frequency k-space data is cut off to lower the resolution, followed by a biexponential T2* decay and the regridding to a radial trajectory; finally, synthetic correlated complex Gaussian noise is added. d) Gridding reconstruction (combined via SOS) of the simulated data set with 20.000 projections.

Fig. 2: Reconstruction of simulated multi-channel breast data with different USFs, combined via SOS. a) Gridding reconstruction. b) TV(2) and c) AnaWeTV(2) CS reconstruction. d) Difference between b) and c). The SSIM value top left assesses the image content by evaluating the present structure; a value of 1 corresponds to complete agreement between the reconstructed image and the ground truth. According to that, the TV(2) leads to better results compared with the gridding reconstruction. The AnaWeTV(2) preserves known tissue boundaries and leads to a further improvement. Especially for high USFs the CS reconstruction shows clear advantages over the gridding reconstruction.

Fig. 3: Reconstruction of measured multi-channel breast data with different USFs, combined via SOS. a) Gridding reconstruction. b) TV(2) and c) AnaWeTV(2) CS reconstruction. d) Difference between b) and c). As for simulated data, the TV(2) and in particular the AnaWeTV(2) lead to an improvement of the image quality, due to the reduction of noise and imaging artifacts and the highlighting of structure. Also here, the CS reconstruction shows its advantage over the gridding reconstruction especially for high USFs.

Proc. Intl. Soc. Mag. Reson. Med. 26 (2018)
1304