Department of Mathematics, University of Pavia, Via Ferrata 5, 27100 Pavia, Italy; Advanced Imaging and Artificial Intelligence Center, Department of Neuroradiology, IRCCS Mondino, Foundation, Via Mondino 2, 27100 Pavia, Italy; INFN, Istituto Nazionale di Fisica Nucleare, Pavia Unit, Via Bassi 6, 27100, Pavia, Italy.
INFN, Istituto Nazionale di Fisica Nucleare, Pavia Unit, Via Bassi 6, 27100, Pavia, Italy; Department of Physics, University of Pavia, Via Bassi 6, 27100, Pavia, Italy.
Comput Methods Programs Biomed. 2024 Nov;256:108399. doi: 10.1016/j.cmpb.2024.108399. Epub 2024 Aug 28.
Magnetic Resonance (MR) parameters mapping in muscle Magnetic Resonance Imaging (mMRI) is predominantly performed using pattern recognition-based algorithms, which are characterised by high computational costs and scalability issues in the context of multi-parametric mapping. Deep Learning (DL) has been demonstrated to be a robust and efficient method for rapid MR parameters mapping. However, its application in mMRI domain to investigate Neuromuscular Disorders (NMDs) has not yet been explored. In addition, data-driven DL models suffered in interpretation and explainability of the learning process. We developed a Physics Informed Neural Network called Myo-Regressor Deep Informed Neural NetwOrk (Myo-DINO) for efficient and explainable Fat Fraction (FF), water-T (wT) and B1 mapping from a cohort of NMDs.A total of 2165 slices (232 subjects) from Multi-Echo Spin Echo (MESE) images were selected as the input dataset for which FF, wT,B1 ground truth maps were computed using the MyoQMRI toolbox. This toolbox exploits the Extended Phase Graph (EPG) theory with a two-component model (water and fat signal) and slice profile to simulate the signal evolution in the MESE framework. A customized U-Net architecture was implemented as the Myo-DINO architecture. The squared L norm loss was complemented by two distinct physics models to define two 'Physics-Informed' loss functions: Cycling Loss 1 embedded a mono-exponential model to describe the relaxation of water protons, while Cycling Loss 2 incorporated the EPG theory with slice profile to model the magnetization dephasing under the effect of gradients and RF pulses. The Myo-DINO was trained with the hyperparameter value of the 'Physics-Informed' component held constant, i.e. λ = 1, while different hyperparameter values (λ) were applied to the squared L norm component in both the cycling loss. In particular, hard (λ=10), normal (λ=1) and self-supervised (λ=0) constraints were applied to gradually decrease the impact of the squared L norm component on the 'Physics Informed' term during the Myo-DINO training process. Myo-DINO achieved higher performance with Cycling Loss 2 for FF, wT and B1 prediction. In particular, high reconstruction similarity and quality (Structural Similarity Index > 0.92, Peak Signal to Noise ratio > 30.0 db) and small reconstruction error (Normalized Root Mean Squared Error < 0.038) to the reference maps were shown with self-supervised weighting of the Cycling Loss 2. In addition muscle-wise FF, wT and B1 predicted values showed good agreement with the reference values. The Myo-DINO has been demonstrated to be a robust and efficient workflow for MR parameters mapping in the context of mMRI. This provides preliminary evidence that it can be an effective alternative to the reference post-processing algorithm. In addition, our results demonstrate that Cycling Loss 2, which incorporates the Extended Phase Graph (EPG) model, provides the most robust and relevant physical constraints for Myo-DINO in this multi-parameter regression task. The use of Cycling Loss 2 with self-supervised constraint improved the explainability of the learning process because the network acquired domain knowledge solely in accordance with the assumptions of the EPG model.
磁共振(MR)参数映射在肌肉磁共振成像(mmMRI)中主要使用基于模式识别的算法进行,这些算法的特点是在多参数映射的背景下计算成本高且可扩展性差。深度学习(DL)已被证明是一种快速进行 MR 参数映射的强大且高效的方法。然而,它在 mMRI 领域中的应用来研究神经肌肉疾病(NMDs)尚未得到探索。此外,数据驱动的 DL 模型在学习过程的解释和可解释性方面存在困难。我们开发了一种基于物理的神经网络,称为 Myo-Regressor Deep Informed Neural NetwOrk(Myo-DINO),用于从 NMD 队列中高效且可解释地进行脂肪分数(FF)、水-T(wT)和 B1 映射。总共选择了 2165 个切片(232 个受试者)来自多回波自旋回波(MESE)图像作为输入数据集,其中使用 MyoQMRI 工具箱计算了 FF、wT、B1 地面真实图。该工具箱利用扩展相位图(EPG)理论和双组件模型(水和脂肪信号)以及切片轮廓来模拟 MESE 框架中的信号演化。作为 Myo-DINO 架构,实现了一个定制的 U-Net 架构。平方 L 范数损失由两个不同的物理模型补充,以定义两个“物理信息”损失函数:循环损失 1 嵌入单指数模型来描述水质子的弛豫,而循环损失 2 结合了 EPG 理论和切片轮廓来模拟梯度和 RF 脉冲作用下的磁化退相。Myo-DINO 是在保持“物理信息”组件的超参数值不变(即 λ=1)的情况下进行训练的,而在两个循环损失中,对平方 L 范数组件应用了不同的超参数值(λ)。特别是,应用了硬(λ=10)、正常(λ=1)和自我监督(λ=0)约束,以在 Myo-DINO 训练过程中逐渐降低平方 L 范数组件对“物理信息”项的影响。Myo-DINO 在使用循环损失 2 进行 FF、wT 和 B1 预测时表现出更高的性能。特别是,自我监督加权循环损失 2 显示出与参考图具有很高的重建相似性和质量(结构相似性指数>0.92,峰值信噪比>30.0 dB)和较小的重建误差(归一化均方根误差<0.038)。此外,肌肉-wise FF、wT 和 B1 预测值与参考值显示出良好的一致性。Myo-DINO 已被证明是一种在 mMRI 环境中进行 MR 参数映射的强大且高效的工作流程。这提供了初步证据,表明它可以成为参考后处理算法的有效替代方法。此外,我们的结果表明,循环损失 2,其中包含扩展相位图(EPG)模型,为 Myo-DINO 在这个多参数回归任务中提供了最稳健和相关的物理约束。使用具有自我监督约束的循环损失 2 提高了学习过程的可解释性,因为网络仅根据 EPG 模型的假设获取领域知识。