Comparison of non-parametric T2 relaxometry methods for myelin water quantification

被引:23
作者
Jorge Canales-Rodriguez, Erick [1 ,2 ,3 ,4 ,5 ]
Pizzolato, Marco [5 ,6 ]
Franco Piredda, Gian [1 ,2 ,5 ,7 ]
Hilbert, Tom [1 ,2 ,5 ,7 ]
Kunz, Nicolas [8 ]
Pot, Caroline [9 ,10 ,11 ,12 ]
Yu, Thomas [5 ,13 ]
Salvador, Raymond [3 ,4 ]
Pomarol-Clotet, Edith [3 ,4 ]
Kober, Tobias [1 ,2 ,5 ,7 ]
Thiran, Jean-Philippe [1 ,2 ,5 ]
Daducci, Alessandro [14 ]
机构
[1] Lausanne Univ Hosp, Dept Radiol, Lausanne, Switzerland
[2] Univ Lausanne, Lausanne, Switzerland
[3] FIDMAG Germanes Hosp Res Fdn, Barcelona, Spain
[4] Ctr Invest Biomed Red Salud Mental CIBERSAM, Barcelona, Spain
[5] Ecole Polytech Fed Lausanne EPFL, Signal Proc Lab LTS5, Lausanne, Switzerland
[6] Tech Univ Denmark, Dept Appl Math & Comp Sci, Lyngby, Denmark
[7] Siemens Healthcare AG, Adv Clin Imaging Technol, Lausanne, Switzerland
[8] Ecole Polytech Fed Lausanne EPFL, Anim Imaging & Technol Sect, Ctr Biomed Imaging CIBM, Lausanne, Switzerland
[9] Geneva Univ Hosp, Dept Pathol & Immunol, Geneva, Switzerland
[10] Univ Geneva, Geneva, Switzerland
[11] Ctr Hosp Univ Vaudois CHUV, Div Neurol, Dept Clin Neurosci, Lausanne, Switzerland
[12] Ctr Hosp Univ Vaudois CHUV, Neurosci Res Ctr, Dept Clin Neurosci, Lausanne, Switzerland
[13] Univ Lausanne, Ctr Biomed Imaging CIBM, Med Image Anal Lab, Lausanne, Switzerland
[14] Univ Verona, Comp Sci Dept, Verona, Italy
基金
瑞士国家科学基金会; 欧盟地平线“2020”;
关键词
T-2; relaxometry; Myelin water imaging; Tikhonov regularization; Non-negative least squares; Tissue microstructure;
D O I
10.1016/j.media.2021.101959
中图分类号
TP18 [人工智能理论];
学科分类号
081104 ; 0812 ; 0835 ; 1405 ;
摘要
Multi-component T-2 relaxometry allows probing tissue microstructure by assessing compartment-specific T-2 relaxation times and water fractions, including the myelin water fraction. Non-negative least squares (NNLS) with zero-order Tikhonov regularization is the conventional method for estimating smooth T-2 distributions. Despite the improved estimation provided by this method compared to non-regularized NNLS, the solution is still sensitive to the underlying noise and the regularization weight. This is especially relevant for clinically achievable signal-to-noise ratios. In the literature of inverse problems, various well-established approaches to promote smooth solutions, including first-order and second-order Tikhonov regularization, and different criteria for estimating the regularization weight have been proposed, such as L-curve, Generalized Cross-Validation, and Chi-square residual fitting. However, quantitative comparisons between the available reconstruction methods for computing the T-2 distribution, and between different approaches for selecting the optimal regularization weight, are lacking. In this study, we implemented and evaluated ten reconstruction algorithms, resulting from the individual combinations of three penalty terms with three criteria to estimate the regularization weight, plus non-regularized NNLS. Their performance was evaluated both in simulated data and real brain MRI data acquired from healthy volunteers through a scan-rescan repeatability analysis. Our findings demonstrate the need for regularization. As a result of this work, we provide a list of recommendations for selecting the optimal reconstruction algorithms based on the acquired data. Moreover, the implemented methods were packaged in a freely distributed toolbox to promote reproducible research, and to facilitate further research and the use of this promising quantitative technique in clinical practice. (C) 2021 The Author(s). Published by Elsevier B.V.
引用
收藏
页数:13
相关论文
共 68 条
[1]   Multi-gradient-echo myelin water fraction imaging: Comparison to the multi-echo-spin-echo technique [J].
Alonso-Ortiz, Eva ;
Levesque, Ives R. ;
Pike, G. Bruce .
MAGNETIC RESONANCE IN MEDICINE, 2018, 79 (03) :1439-1446
[2]   MRI-Based Myelin Water Imaging: A Technical Review [J].
Alonso-Ortiz, Eva ;
Levesque, Ives R. ;
Pike, G. Bruce .
MAGNETIC RESONANCE IN MEDICINE, 2015, 73 (01) :70-81
[3]   Symmetric diffeomorphic image registration with cross-correlation: Evaluating automated labeling of elderly and neurodegenerative brain [J].
Avants, B. B. ;
Epstein, C. L. ;
Grossman, M. ;
Gee, J. C. .
MEDICAL IMAGE ANALYSIS, 2008, 12 (01) :26-41
[4]   Proof that gmT2 is the Reciprocal of gmR2 [J].
Bjarnason, Thorarin A. .
CONCEPTS IN MAGNETIC RESONANCE PART A, 2011, 38A (03) :128-131
[5]   Quantitative T2 Analysis: The Effects of Noise, Regularization, and Multivoxel Approaches [J].
Bjarnason, Thorarin A. ;
McCreary, Cheryl R. ;
Dunn, Jeff F. ;
Mitchell, J. Ross .
MAGNETIC RESONANCE IN MEDICINE, 2010, 63 (01) :212-217
[6]  
Bochkarev Vladimir, 2018, Journal of Physics: Conference Series, V1141, DOI 10.1088/1742-6596/1141/1/012148
[7]   Use of the NESMA Filter to Improve Myelin Water Fraction Mapping with Brain MRI [J].
Bouhrara, Mustapha ;
Reiter, David A. ;
Maring, Michael C. ;
Bonny, Jean-Marie ;
Spencer, Richard G. .
JOURNAL OF NEUROIMAGING, 2018, 28 (06) :640-649
[8]   ALGORITHM WITH GUARANTEED CONVERGENCE FOR FINDING A ZERO OF A FUNCTION [J].
BRENT, RP .
COMPUTER JOURNAL, 1971, 14 (04) :422-&
[9]   Noise reduction in multiple-echo data sets using singular value decomposition [J].
Bydder, Mark ;
Du, Jiang .
MAGNETIC RESONANCE IMAGING, 2006, 24 (07) :849-856
[10]   L-curve and curvature bounds for Tikhonov regularization [J].
Calvetti, D ;
Reichel, L ;
Shuibi, A .
NUMERICAL ALGORITHMS, 2004, 35 (2-4) :301-314