对阿尔茨海默病的扩散磁共振成像处理方法
阅读说明：本技术 对阿尔茨海默病的扩散磁共振成像处理方法 (Diffusion magnetic resonance imaging processing method for Alzheimer disease ) 是由 马国林 杜雷 高文文 刘秀秀 王依格 于 2020-05-25 设计创作，主要内容包括：对阿尔茨海默病的扩散磁共振成像处理方法,首先在3D T1图像上勾画出感兴趣区ROI,将勾画好的T1图像配准到dMRI图像中,随后计算出双侧海马的α,H,μ和ADC的平均值,然后使用FSL工具对获得的图像进行头部运动和涡流畸变校正,使用b值为0和954s/mm～(2)时获得的图像计算ADC值,采用FM模型对图像进行分析以获取异常扩散参数,根据基于FM模型的dMRI理论获得由扩散引起的信号衰减,再进行统计学分析获得结果,因此能够较好地反映生物活细胞扩散过程,达到真正信号衰减曲线和拟合曲线之间的最佳一致性。(A diffusion magnetic resonance imaging processing method for Alzheimer disease comprises the steps of firstly delineating a region of interest ROI on a 3D T1 image, registering the delineated T1 image into a dMRI image, then calculating the average values of alpha, H, mu and ADC of bilateral hippocampus, then using an FSL tool to correct head movement and eddy current distortion of the obtained image, and using b values of 0 and 954s/mm 2 The ADC value of the obtained image is calculated, the FM model is adopted to analyze the image so as to obtain abnormal diffusion parameters, signal attenuation caused by diffusion is obtained according to a dMRI theory based on the FM model, and then statistical analysis is carried out to obtain a result, so that the diffusion process of the living biological cells can be better reflected, and the optimal consistency between a real signal attenuation curve and a fitting curve is achieved.)
技术领域
本发明涉及医学图像处理的技术领域，尤其涉及一种对阿尔茨海默病的扩散磁共振成像处理方法。
背景技术
阿尔茨海默病(Alzheimer's disease,AD)是一种复杂而常见的神经退行性疾病，其特征是记忆丧失和认知功能下降。根据世界卫生组织(World Health Organization,WHO)的数据，痴呆影响了全世界近4750万人，并且每年仍会增加大约770万新病例。阿尔茨海默病是最常见的痴呆类型，在这些病例中可能占60-70％，阿尔茨海默病对患者的日常生活造成了重大影响。AD的发病机制非常复杂，主要包括β淀粉样蛋白(Aβ)和过度磷酸化tau蛋白的形成分别导致Aβ斑块和神经原纤维缠结形成，进而导致神经元死亡。现在AD的诊断较复杂并且准确率不高，因此在临床研究中开发一种有效的AD诊断方法十分重要。
扩散磁共振成像(diffusion magnetic resonance imaging,dMRI)是一种可以描述生物组织中水分子随机运动的强大工具，它能够提供关于组织微观特性的独特信息，并对检测大脑灰质和白质的变化高度敏感。水分子在神经系统中的扩散过程是有方向性的。这种方向依赖，即各向异性，主要是由于固有的轴突膜阻碍水分子的扩散和轴突的密集堆积导致。在微米水平上测量水分子扩散的各向异性能够反映微观结构的变化。
与其他磁共振成像方法相比，dMRI研究细胞水平水的扩散过程，远优于传统的毫米级分辨率图像。目前临床上应用最广泛的dMRI技术之一是表观扩散系数(apparentdiffusion coefficient,ADC)，dMRI获得的b值在0～1000s/mm2之间。ADC的临床应用广泛，可用于鉴别脑肿瘤和对肿瘤进行分级，也可用于鉴别轻度认知障碍(mild cognitiveimpairment,MCI)和AD患者。扩散张量成像(diffusion tensor imaging,DTI)是研究中广泛应用的dMRI技术之一。从DTI中获得的分数各向异性(fractional anisotropy,FA)和平均扩散率(mean diffusivity,MD)是MCI和AD患者脑异常的潜在生物标志物。然而，ADC和DTI都有一定局限性。首先，传统的dMRI采用单指数模型，它假设生物组织中存在规则的扩散过程。然而许多研究发现观测到的dMRI信号衰减曲线偏离了单指数形式。而ADC和DTI均基于常规的单指数扩散模型，观测到的扩散时间反映了扩散的非高斯性质。二是DTI张量模型过于简单，这意味着其各项指标受微观结构的许多特征的影响。为了解决这一问题，许多基于不同理论和反常扩散过程的模型被开发出来，以找到真正信号衰减曲线和拟合曲线之间的最佳一致性。这些模型能对疾病类型和分级进行更详细的检测，但其信号衰减曲线仍不太匹配。
发明内容
为克服现有技术的缺陷，本发明要解决的技术问题是提供了一种对阿尔茨海默病的扩散磁共振成像处理方法，其能够较好地反映生物活细胞扩散过程，达到真正信号衰减曲线和拟合曲线之间的最佳一致性。
本发明的技术方案是：这种对阿尔茨海默病的扩散磁共振成像处理方法，其包括以下步骤：
(1)图像获取：在磁共振扫描仪上，使用8通道头线圈，使用Stejskal-Tanner单次自旋回波平面成像序列获得所有受试者的dMRI图像；
(2)左右侧海马被选为感兴趣区ROI，首先在3D T1图像上勾画出ROI，然后将勾画好的T1图像配准到dMRI图像中，随后计算出双侧海马的α,H,μ和ADC的平均值，α指诺亚指数，量化随机过程的波动，当α＝2时增量是高斯分布，而当0<α<2时增量是Lévy分配；Η是赫斯特指数，描述分子轨道的自相似性；μ是记忆参数，且μ＝Η-1/α，当μ>0时，增量呈正相关且为长程依赖，而当μ<0，增量是负相关且为短程依赖；
(3)使用FSL工具对获得的图像进行头部运动和涡流畸变校正，使用b值为0和954s/mm2时获得的图像计算ADC值，采用FM模型对图像进行分析以获取异常扩散参数，根据基于FM模型的dMRI理论，由扩散引起的信号衰减通过公式(1)计算：
Dα,Η指反常扩散的扩散系数，γ指旋磁比，G0指扩散梯度幅值，Δ指梯度分离时间，η指一个无量纲数，是由α,H,δ和Δ确定的；
(4)采用单因素方差分析比较轻度AD患者、中度AD患者和健康对照者这三组在年龄、受教育程度和MMSE评分方面的差异，数据用均数±标准差来表示，性别采用卡方检验进行比较，P值<0.05为差异有统计学意义。
本发明首先在3D T1图像上勾画出感兴趣区ROI，将勾画好的T1图像配准到dMRI图像中，随后计算出双侧海马的α,H,μ和ADC的平均值，然后使用FSL工具对获得的图像进行头部运动和涡流畸变校正，使用b值为0和954s/mm2时获得的图像计算ADC值，采用FM模型对图像进行分析以获取异常扩散参数，根据基于FM模型的dMRI理论获得由扩散引起的信号衰减，再进行统计学分析获得结果，因此能够较好地反映生物活细胞扩散过程，达到真正信号衰减曲线和拟合曲线之间的最佳一致性。
附图说明
图1示出了左右侧海马被选为感兴趣区。
图2示出了一个AD患者(第一行，58岁男性)和一个健康对照组(第二行，60岁男性)的轴位T1加权图像、α、H、ADC图，肉眼所见双侧海马在AD患者和健康对照者之间未见明显差异。
图3示出了AD患者与健康对照组之间的比较。散点图(A)和表(B)显示α和ADC值可以鉴别AD患者和健康对照组，而H和μ值在鉴别二组时未见明显统计学差异。AD组共24人，健康对照组共11人。采用双样本t检验，P<0.05为差异有统计学意义。
图4示出了轻度AD和中度AD患者之间的比较。散点图(A)和表(B)显示α和ADC值可以鉴别二组，而H和μ值在鉴别二组时未见明显统计学差异。轻度AD患者共12人，中度AD患者共12人。采用双样本t检验，P<0.05为差异有统计学意义。
图5示出了由α,ADC,α+ADC组成的ROC曲线图以鉴别AD患者(n＝24)和健康对照者(n＝11)。如图所示α+ADC的曲线下面积(AUC＝0.848,左侧ROI；AUC＝0.856,右侧ROI)高于单独α和单独ADC。
图6示出了由α,ADC,α+ADC组成的ROC曲线图以鉴别轻度AD患者(n＝12)和中度AD患者(n＝12)。如图所示α+ADC的曲线下面积(AUC＝0.861,左侧ROI；AUC＝0.868,右侧ROI)高于单独α和单独ADC。
图7示出了AD患者右侧海马的α值和MMSE分值的相关性。进行皮尔逊偏相关分析。n＝24。*,p<0.05。
图8示出了AD患者双侧海马的α值和MoCA分值的相关性。进行皮尔逊偏相关分析。n＝24。*,p<0.05。
图9示出了根据本发明的对阿尔茨海默病的扩散磁共振成像处理方法的流程图。
具体实施方式
如图1所示，这种对阿尔茨海默病的扩散磁共振成像处理方法，其包括以下步骤：
(1)图像获取：在磁共振扫描仪(3.0T磁共振扫描仪，GE Healthcare,DiscoveryMR750,USA)上，使用8通道头线圈，使用Stejskal-Tanner单次自旋回波平面成像序列获得所有受试者的dMRI图像；
(2)左右侧海马被选为感兴趣区ROI，首先在3D T1图像上勾画出ROI，然后将勾画好的T1图像配准到dMRI图像中，目的是更准确的确定海马的边界，排除不明确的体素，随后计算出双侧海马的α,H,μ和ADC的平均值，α指诺亚指数，量化随机过程的波动，当α＝2时增量是高斯分布，而当0<α<2时增量是
Lévy分配；Η是赫斯特指数，描述分子轨道的自相似性；μ是记忆参数，且μ＝Η-1/α，当μ>0时，增量呈正相关且为长程依赖，而当μ<0，增量是负相关且为短程依赖；
(3)使用FSL工具对获得的图像进行头部运动和涡流畸变校正，使用b值为0和954s/mm2时获得的图像计算ADC值，采用FM模型对图像进行分析以获取异常扩散参数，根据基于FM模型的dMRI理论，由扩散引起的信号衰减通过公式(1)计算：
Dα,Η指反常扩散的扩散系数，γ指旋磁比，G0指扩散梯度幅值，Δ指梯度分离时间，η指一个无量纲数，是由α,H,δ和Δ确定的；
(4)采用单因素方差分析比较轻度AD患者、中度AD患者和健康对照者这三组在年龄、受教育程度和MMSE评分方面的差异，数据用均数±标准差来表示，性别采用卡方检验进行比较，P值<0.05为差异有统计学意义。
本发明首先在3D T1图像上勾画出感兴趣区ROI，将勾画好的T1图像配准到dMRI图像中，随后计算出双侧海马的α,H,μ和ADC的平均值，然后使用FSL工具对获得的图像进行头部运动和涡流畸变校正，使用b值为0和954s/mm2时获得的图像计算ADC值，采用FM模型对图像进行分析以获取异常扩散参数，根据基于FM模型的dMRI理论获得由扩散引起的信号衰减，再进行统计学分析获得结果，因此能够较好地反映生物活细胞扩散过程，达到真正信号衰减曲线和拟合曲线之间的最佳一致性。
优选地，所述步骤(1)中，在常规dMRI序列扫描过程中，为了拟合FM模型没有固定扩散梯度分离时间，扩散梯度分离时间Δ依次为27.060ms，39.560ms，52.060ms，每个Δ值的扩散梯度幅值G0依次为15.67、19.68、24.73、31.06、39.01和49.00mT/m，目的是在对数轴上近似均匀间隔，梯度时间常数δ固定为20.676ms；因此，在每个梯度方向获得了18个非零b值：151、239、377、595、939、1481、245、387、611、964、1521、2399、339、535、845、1333、2103和3317s/mm2；为了减小扩散各向异性的影响，在三个正交方向x轴、y轴和z轴上应用了扩散梯度，共获得12张b＝0的图像。
优选地，所述步骤(1)中，dMRI序列参数：重复时间TR＝3800ms；回波时间TE＝110ms；翻转角FA＝90°；层厚＝5.0mm；视场FOV＝240mm×240mm；矩阵大小＝128×128；加速因子＝2；激发次数＝2；总扫描时间为8分33秒；T1成像参数：TR＝6.7ms；TE＝Min Full；FA＝12°；层厚＝1.0mm；视场FOV＝256mm×256mm；矩阵大小＝256×256，扫描时间＝4分10秒。
优选地，所述步骤(2)中，由一名有经验的放射科医生在3D T1图像上使用MRICRON软件手动勾画出ROI。
优选地，所述步骤(3)中，所有拟合过程均采用MATLAB中的trust-region-reflection非线性拟合算法进行。
优选地，所述步骤(4)中，采用双样本t检验比较α,Η,μ和ADC值在AD组和健康对照组之间的差异及在轻度AD组和中度AD组之间的差异。
优选地，所述步骤(4)中，为了量化α,Η,μ和ADC值在鉴别AD组和健康对照组和鉴别轻度AD组和中度AD组的敏感度和特异度，采用受试者工作特征曲线ROC(receiveroperating characteristic)进行评估，并计算了曲线下面积AUC(area under thecurve)。
优选地，所述步骤(5)中，采用皮尔逊偏相关分析获得AD患者的α,Η,μ和ADC值与MMSE和MoCA量表分值之间的相关性。
以下更详细地说明本发明的实施例。
所有受试者的人口统计学信息及临床信息汇总到表1。本发明共纳入24例AD患者(男性9例，女性15例，平均年龄68.96岁，范围50-79岁)和11例健康对照者(男性2例，女性9例，平均年龄65.27岁，范围54-78岁)。AD组又分为轻度AD组和中度AD组。这三组患者在年龄、受教育程度方面无统计学差异(P>0.05)。三组患者的MMSE分值有统计学差异(P<0.05)。
表1
图1显示了轴位T1图像中的双侧海马。图2显示了AD组和健康对照组的α,Η,ADC图像和T1图像，肉眼所见双侧海马在AD患者和健康对照者之间未见明显差异。
表2
所有受试者双侧海马的FM相关参数值及ADC值均汇总到表2，用均数±标准差表示。图3和图4分别显示了鉴别AD患者与健康对照者的散点图以及鉴别轻度AD和中度AD患者的散点图。如图3所示，α值(左侧ROI，P＝0.011；右侧ROI，P＝0.001)和ADC值(左侧ROI，P＝0.001；右侧ROI，P＝0.001)可以鉴别AD患者和健康对照者。此外，图4也显示α值(右侧ROI，P＝0.015)和ADC值(左侧ROI，P＝0.011；右侧ROI，P＝0.022)可以鉴别轻度AD和中度AD患者。而Η和μ值未能区分AD组及健康对照者、轻度AD和中度AD患者(P>0.05)。
ROC曲线分析也显示出鉴别AD患者与健康对照者的能力。图5为鉴别AD患者与健康对照者的ROC曲线，如图所见，相比于α(AUC＝0.78左侧ROI，AUC＝0.811右侧ROI)和ADC(AUC＝0.847左侧ROI，AUC＝0.833右侧ROI)，α+ADC(AUC＝0.848左侧ROI，AUC＝0.856右侧ROI)显示出更大的曲线下面积，当α的阈值为1.5939(左侧ROI)和1.58415(右侧ROI)时，显示出最高的约登指数(敏感性+特异性-1)，此时的敏感度和特异度分别为0.917和0.636(左侧ROI)、0.917和0.727(右侧ROI)。同样地，对于鉴别轻度AD和中度AD患者(图6)，相比于α(AUC＝0.813右侧ROI)与ADC(AUC＝0.792左侧ROI，AUC＝0.806右侧ROI)，α+ADC(AUC＝0.861左侧ROI，AUC＝0.868右侧ROI)能获得更大的曲线下面积；当α＝1.55115(右侧ROI)时能获得最高的约登指数，此时敏感度和特异度分别为0.917和0.75。
在AD患者中FM相关参数α值与MMSE分值(P<0.05；图7)和MoCA分值(P<0.05；图8)呈正相关。其他FM相关参数值和ADC值与MMSE分值和MoCA分值无显著相关性。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何形式上的限制，凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属本发明技术方案的保护范围。
- 上一篇：一种医用注射器针头装配设备
- 下一篇：磁共振摄像装置