一种利用自回归模型计算磁共振图像表观弥散系数的方法
【专利摘要】本发明涉及一种利用自回归模型计算磁共振图像表观弥散系数的方法，利用自回归模型计算生成ADC磁共振图像。利用多个且等间隔b值DWI序列可采集得到多幅不同b值的磁共振图像，对采集得到的所有像素的磁共振信号进行线性积分，即：利用Simpson数值积分法则对信号点做关于b值的多段等间隔积分，再建立自回归模型，最后对自回归模型参数进行最大似然估计得到ADC。是一种更加精确快速的计算E指数衰减系数的算法，该算法可以用于磁共振成像技术中快速精确计算ADC图像。本发明能够显著ADC磁共振图像计算效率和精确性，适用于放射科临床上各种相关疾病的诊断和研究。
【专利说明】一种利用自回归模型计算磁共振图像表观弥散系数的方法
【技术领域】
[0001] 本发明涉及一种磁共振成像技术，特别涉及一种利用自回归模型计算磁共振图像 表观弥散系数的方法。
【背景技术】
[0002] 磁共振成像具有无损伤、无电离辐射、软组织对比度高、任意方向断层等特点，目 前已广泛用于医学临床诊断。
[0003] 磁共振成像中的表观扩散系数ADC是在DWI成像中衡量生物体内水分子扩散程度 的物理量。DWI成像中通过改变扩散梯度来改变b值，从而影响DWI图像信号。DWI图像信 号随不同b值呈E指数衰减形式，ADC为衰减系数，是衡量评估生物体组织中水分子扩散强 弱，能够间接反映生物体组织结构。在临床上医学，ADC图像被广泛用于诊断脑梗死、肿瘤、 神经损坏等疾病。
[0004] 如果对采集到的多个b值的DWI图像做E指数函数拟合，可以计算得到ADC磁共 振图像。目前已有的拟合E指数函数计算衰减系数的方法Levenberg-Marquardt (列文伯 格-马夸尔特法）非线性拟合方法、Log-Linear线性拟合方法以及DISC数值积分方法。
[0005] 上述的Levenberg-Marquardt非线性拟合的方法是基于非线性迭代求解最小二 乘方法，该方法总体上是精确性较高的方法，但是该方法的缺点时间效率低下，耗时长，对 于计算一份三维图像数据，由于要逐个像素进行拟合计算，通常需要做几百万次甚至千万 次非线性拟合，所需时间非常长。此外该方法在拟合前需要对初值进行估计，而初值估计准 确性往往会影响拟合结果的收敛，因此进一步限制了该方法的准确性和实际应用性。
[0006] Log-Linear线性拟合方法是对所有采集的E指数衰减信号求自然对数log，将非 线性数据转变为线性数据，然后做线性拟合得到衰减系数的方法。该方法的优点是计算快 且易于实现，然而该方法对噪声相当敏感，对于信噪比低或信号衰减快的数据，拟合结果的 误差相当大。
[0007] DISC数值积分方法是对整个采集到的E指数衰减信号数据求数值定积分，然后除 以首末两个信号的差值来计算得到衰减系数的方法。该方法的优点是准确、快速且易于实 现。但是由于该方法缺乏对误差最小化的数学分析，在精度上会有一定的局限。
【发明内容】
[0008] 本发明是针对现有的几种拟合计算E指数衰减系数的算法各自的时间长、依赖初 值估计、噪声敏感性强、精度不高的问题，提出了一种利用自回归模型计算磁共振图像表观 弥散系数的方法，利用对E指数衰减形式的物理信号建立自回归模型，然后做最大似然参 数估算，得到衰减系数的拟合方法。该方法可应用用于计算生成表观弥散系数(Apparent Diffusion Coefficient, ADC)磁共振图像。是一种更加精确快速的计算E指数衰减系数的 算法，该算法可以用于磁共振成像技术中快速精确计算ADC图像。
[0009] 本发明的技术方案为：一种利用自回归模型计算磁共振图像表观弥散系数的方 法，具体包括如下步骤： 1) 用多个且等间隔b值DWI成像序列采集得到多幅不同b值的磁共振图像； 2) 对采集的多b值DWI图像信号：从第一个b值到第三个b值开始，进行两个b值间 隔长度(相邻三个数据点）的Simpson法则数值积分，然后从第二个b值到第四个b值进行 两个b值间隔长度(相邻三个数据点）的Simpson法则数值积分，依次下去，从而形成一个 DWI积分数值随b值变化的信号序列； 3) 将所述DWI信号积分的序列建立成自回归模型，最后一个b值DWI信号值可表不成 由邻近的前两个b值DWI信号值的线性组合形式； 4) 对步骤3)建立的自回归模型进行最大似然估计从而求解系数得到表观弥散系数 ADC 值； 5) 对步骤1)的磁共振图像中被试体内所有像素计算ADC值，得到所述磁共振ADC图 像。
[0010] 所述步骤2)中Simpson法则数值积分表示成公式：
【权利要求】
1. 一种利用自回归模型计算磁共振图像表观弥散系数的方法，其特征在于，具体包括 如下步骤： 1) 用多个且等间隔b值DWI成像序列采集得到多幅不同b值的磁共振图像； 2) 对采集的多b值DWI图像信号：从第一个b值到第三个b值开始，进行两个b值间 隔长度(相邻三个数据点）的Simpson法则数值积分，然后从第二个b值到第四个b值进行 两个b值间隔长度(相邻三个数据点）的Simpson法则数值积分，依次下去，从而形成一个 DWI积分数值随b值变化的信号序列； 3) 将所述DWI信号积分的序列建立成自回归模型，最后一个b值DWI信号值可表不成 由邻近的前两个b值DWI信号值的线性组合形式； 4) 对步骤3)建立的自回归模型进行最大似然估计从而求解系数得到表观弥散系数 ADC 值； 5) 对步骤1)的磁共振图像中被试体内所有像素计算ADC值，得到所述磁共振ADC图 像。
2. 根据权利要求1所述利用自回归模型计算磁共振图像表观弥散系数的方法，其特征 在于，所述步骤2)中Simpson法则数值积分表示成公式：
其中W表示所设的相等的相邻b值间隔值，表示扩散梯度系数为b时采集的图 像信号Λ表示为第i个b值，iU，……，N-2, N表示b值总数。
3. 根据权利要求2所述利用自回归模型计算磁共振图像表观弥散系数的方法，其特征 在于，所述步骤3)中DWI信号积分的序列建立成自回归模型表示成公式：
4. 根据权利要求3所述利用自回归模型计算磁共振图像表观弥散系数的方法，其特征 在于，所述步骤4)中对所述的自回归模型进行最大似然估计表示成公式：
最后对ADC求解得到的解析式表示成公式：
」抵表示所设的相等的相邻b值对应的回波时间差值。
【文档编号】G06T11/00GK104095635SQ201410361347
【公开日】2014年10月15日 申请日期:2014年7月28日 优先权日:2014年7月28日
【发明者】王丽嘉, 裴孟超, 董芳, 李建奇, 聂生东, 王远军 申请人:上海理工大学