一种利用自回归模型计算磁共振图像表观弥散系数的方法.pdf

1、(19)中华人民共和国国家知识产权局 (12)发明专利 (10)授权公告号 (45)授权公告日 (21)申请号 201410361347.9 (22)申请日 2014.07.28 (73)专利权人 上海理工大学 地址 200093 上海市杨浦区军工路516号 (72)发明人 王丽嘉 裴孟超 董芳 李建奇 聂生东 王远军 (74)专利代理机构 上海申汇专利代理有限公司 31001 代理人 吴宝根 (51)Int.Cl. A61B 5/055(2006.01) G06T 11/00(2006.01) (56)对比文件 CN 101953687 A,2011.01.26,全文. JP 特開2009-

2、285120 A,2009.12.10,全文. JP 特開平6-54816 A,1994.03.01,全文. WO 2008/147923 A1,2008.12.04,全文. 张兆岩等.脑磁共振扩散张量成像最佳扫描 参数的优选. 中华保健医学杂志 .2013,第15卷 (第4期),328-331. Lijun Guo et alDual-source parallel RF transmission for diffusion-weighted imaging of the abdomen using different b values: Image quality and apparent

3、 diffusion coefficient comparison with conventional single-source transmission. Journal of Magnetic Resonance Imaging .2013,第37卷(第4期),875885. 张月浪等.正常腹部实质脏器磁共振弥散加 权成像ADC值和b值研究. 中国医学影像学杂 志 .2006,第14卷(第4期),265-268.(续) 审查员 洪虹 (54)发明名称 一种利用自回归模型计算磁共振图像表观 弥散系数的方法 (57)摘要 本发明涉及一种利用自回归模型计算磁共 振图像表观弥散系数的方法， 利用

4、自回归模型计 算生成ADC磁共振图像。 利用多个且等间隔b值 DWI序列可采集得到多幅不同b值的磁共振图像， 对采集得到的所有像素的磁共振信号进行线性 积分， 即： 利用Simpson数值积分法则对信号点做 关于b值的多段等间隔积分， 再建立自回归模型， 最后对自回归模型参数进行最大似然估计得到 ADC。 是一种更加精确快速的计算E指数衰减系数 的算法， 该算法可以用于磁共振成像技术中快速 精确计算ADC图像。 本发明能够显著ADC磁共振图 像计算效率和精确性， 适用于放射科临床上各种 相关疾病的诊断和研究。 转续页 权利要求书1页 说明书4页 附图1页 CN 104095635 B 2016

5、.07.20 CN 104095635 B (56)对比文件 Dariya Malyarenko et alMulti-system repeatability and reproducibility of apparent diffusion coefficient measurement using an ice-water phantom. Journal of Magnetic Resonance Imaging .2013,第37卷(第5期),1238-1246. So Yeon Kim et alReproducibility of measurement of apparent

6、diffusion coefficients of malignant hepatic tumors: Effect of DWI techniques and calculation methods. Journal of Magnetic Resonance Imaging .2012,第36卷(第5期),11311138. 2/2 页 2 接上页 CN 104095635 B 1.一种利用自回归模型计算磁共振图像表观弥散系数的方法， 其特征在于， 具体包括 如下步骤： 1)用多个且等间隔b值DWI成像序列采集得到多幅不同b值的磁共振图像； 2)对采集的多b值DWI图像信号： 相邻三个数据

7、点， 从第一个b值到第三个b值开始， 进行 两个b值间隔长度的Simpson法则数值积分， 然后从第二个b值到第四个b值进行两个b值间 隔长度的Simpson法则数值积分， 依次下去， 从而形成一个DWI积分数值随b值变化的信号序 列； 3)将步骤2） 形成的DWI积分数值随b值变化的信号序列建立成自回归模型， 最后一个b 值DWI信号值可表示成由邻近的前两个b值DWI信号值的线性组合形式； 4)对步骤3） 建立的自回归模型进行最大似然估计从而求解系数得到表观弥散系数ADC 值； 5)对步骤1） 的磁共振图像中被试体内所有像素计算ADC值， 得到磁共振ADC图像。 2.根据权利要求1所述利用自

8、回归模型计算磁共振图像表观弥散系数的方法， 其特征 在于， 所述步骤2） 中Simpson法则数值积分表示成公式： 其中b表示所设的相等的相邻b值间隔值，表示扩散梯度系数为b时采集的图像 信号； 表示为第i个b值， i=1， 2， ， N-2， N表示b值总数。 3.根据权利要求2所述利用自回归模型计算磁共振图像表观弥散系数的方法， 其特征 在于， 所述步骤3)中DWI积分数值随b值变化的信号序列建立成自回归模型表示成公式： 。 4.根据权利要求3所述利用自回归模型计算磁共振图像表观弥散系数的方法， 其特征 在于， 所述步骤4) 中对所述的自回归模型进行最大似然估计表示成公式： ； 最后对AD

9、C求解得到的解析式表示成公式： ， 其中 ， TE表示所设的相等的相邻b值对应的回波时间差值。 权 利 要 求 书 1/1 页 2 CN 104095635 B 3 一种利用自回归模型计算磁共振图像表观弥散系数的方法 技术领域 0001 本发明涉及一种磁共振成像技术， 特别涉及一种利用自回归模型计算磁共振图像 表观弥散系数的方法。 背景技术 0002 磁共振成像具有无损伤、 无电离辐射、 软组织对比度高、 任意方向断层等特点， 目 前已广泛用于医学临床诊断。 0003 磁共振成像中的表观扩散系数ADC是在DWI成像中衡量生物体内水分子扩散程度 的物理量。 DWI成像中通过改变扩散梯度来改变b值

10、， 从而影响DWI图像信号。 DWI图像信号随 不同b值呈E指数衰减形式， ADC为衰减系数， 是衡量评估生物体组织中水分子扩散强弱， 能 够间接反映生物体组织结构。 在临床上医学， ADC图像被广泛用于诊断脑梗死、 肿瘤、 神经损 坏等疾病。 0004 如果对采集到的多个b值的DWI图像做E指数函数拟合， 可以计算得到ADC磁共振图 像。 目前已有的拟合E指数函数计算衰减系数的方法Levenberg-Marquardt （列文伯格-马夸 尔特法） 非线性拟合方法、 Log-Linear线性拟合方法以及DISC数值积分方法。 0005 上述的Levenberg-Marquardt非线性拟合的方

11、法是基于非线性迭代求解最小二乘 方法， 该方法总体上是精确性较高的方法， 但是该方法的缺点时间效率低下， 耗时长， 对于 计算一份三维图像数据， 由于要逐个像素进行拟合计算， 通常需要做几百万次甚至千万次 非线性拟合， 所需时间非常长。 此外该方法在拟合前需要对初值进行估计， 而初值估计准确 性往往会影响拟合结果的收敛， 因此进一步限制了该方法的准确性和实际应用性。 0006 Log-Linear线性拟合方法是对所有采集的E指数衰减信号求自然对数log， 将非线 性数据转变为线性数据， 然后做线性拟合得到衰减系数的方法。 该方法的优点是计算快且 易于实现， 然而该方法对噪声相当敏感， 对于信噪

12、比低或信号衰减快的数据， 拟合结果的误 差相当大。 0007 DISC数值积分方法是对整个采集到的E指数衰减信号数据求数值定积分， 然后除 以首末两个信号的差值来计算得到衰减系数的方法。 该方法的优点是准确、 快速且易于实 现。 但是由于该方法缺乏对误差最小化的数学分析， 在精度上会有一定的局限。 发明内容 0008 本发明是针对现有的几种拟合计算E指数衰减系数的算法各自的时间长、 依赖初 值估计、 噪声敏感性强、 精度不高的问题， 提出了一种利用自回归模型计算磁共振图像表观 弥散系数的方法， 利用对E指数衰减形式的物理信号建立自回归模型， 然后做最大似然参数 估算， 得到衰减系数的拟合方法。

13、 该方法可应用用于计算生成表观弥散系数 （Apparent Diffusion Coefficient, ADC） 磁共振图像。 是一种更加精确快速的计算E指数衰减系数的 算法， 该算法可以用于磁共振成像技术中快速精确计算ADC图像。 0009 本发明的技术方案为： 一种利用自回归模型计算磁共振图像表观弥散系数的方 说 明 书 1/4 页 3 CN 104095635 B 4 法， 具体包括如下步骤： 0010 1) 用多个且等间隔b值DWI成像序列采集得到多幅不同b值的磁共振图像； 0011 2) 对采集的多b值DWI图像信号： 从第一个b值到第三个b值开始， 进行两个b值间 隔长度 （相邻

14、三个数据点） 的Simpson法则数值积分， 然后从第二个b值到第四个b值进行两 个b值间隔长度 （相邻三个数据点） 的Simpson法则数值积分， 依次下去， 从而形成一个DWI积 分数值随b值变化的信号序列； 0012 3) 将所述DWI信号积分的序列建立成自回归模型， 最后一个b值DWI信号值可表示 成由邻近的前两个b值DWI信号值的线性组合形式； 0013 4) 对步骤3） 建立的自回归模型进行最大似然估计从而求解系数得到表观弥散系 数ADC值； 0014 5) 对步骤1） 的磁共振图像中被试体内所有像素计算ADC值， 得到所述磁共振ADC 图像。 0015 所述步骤2） 中Simps

15、on法则数值积分表示成公式： 0016 0017其中b表示所设的相等的相邻b值间隔值，表示扩散梯度系数为b时采集的 图像信号；表示为第i个b值， i=1， 2， ， N-2， N表示b值总数。 0018 所述步骤3)中DWI信号积分的序列建立成自回归模型表示成公式： 0019 0020 所述步骤4) 中对所述的自回归模型进行最大似然估计表示成公式： 0021； 0022 最后对ADC求解得到的解析式表示成公式： 0023 ， 其中 ， TE表示所设的相等的相邻b值对应的回波时间差值。 0024 本发明的有益效果在于： 本发明利用自回归模型计算磁共振图像表观弥散系数的 方法， 与Levenber

16、g-Marquardt方法相比， 计算速度快数百倍以上， 不需要提供任何初值估 计， 且精度略高； 与Log-Linear方法相比， 计算速度快6-10倍， 抗噪能力强， 精度高很多； 与 DISC数值积分方法比， 计算速度相当， 但是精度更高。 因此， 本发明方法能够显著ADC磁共振 图像计算效率和精确性， 适用于放射科临床上各种相关疾病的诊断和研究。 附图说明 0025 图1为本发明方法得到的人体大脑的ADC磁共振图像图。 具体实施方式 0026 本发明方法关于计算生成ADC磁共振图像包括以下步骤： 说 明 书 2/4 页 4 CN 104095635 B 5 0027 1) 用多个且等间

17、隔b值DWI成像序列采集得到多幅不同b值的磁共振图像； 0028 2) 对采集的多b值DWI图像信号： 从第一个b值到第三个b值开始， 进行两个b值间 隔长度 （相邻三个数据点） 的Simpson法则数值积分， 然后从第二个b值到第四个b值进行两 个b值间隔长度 （相邻三个数据点） 的Simpson法则数值积分， 依次下去， 从而形成一个DWI积 分数值随b值变化的信号序列； 0029 例如， 在b = 0,100,200,300时分别分别采集得到四幅图像， 就对b从0到200积分， 这时候积分用到的是对应b = 0,100,200的第一、 二、 三幅图像信号； 然后再从100到300积 分，

18、 这时候积分用到的是对应b =100,200,300第二、 三、 四幅图像信号； 0030 3) 将所述DWI信号积分的序列建立成自回归模型， 最后一个b值DWI信号值可表示 成由邻近的前两个b值DWI信号值的线性组合形式； 0031 4) 对所述的自回归模型进行最大似然估计从而求解系数得到表观弥散系数ADC 值； 0032 5) 对步骤1） 的磁共振图像中被试体内所有像素计算ADC值， 得到所述磁共振ADC 图像。 0033 其中对步骤2)中进行DWI信号数值积分表示成公式： 0034 （1） 0035其中， 表示第i个积分值；表示扩散梯度系数为b时采集的图像信号； 表示 为第i个b值；，

19、表示为第i个和第i-2个b值DWI图像 之间的信号差值， N表示b值总数， 公式 （1） 是积分的连续形式， 一般可得到解析解， 但我们采 集得到的不同b值下的图像是离散的， 所以在解决实际问题中时， 我们需要对公式 （1） 进行 离散表示， 我们用Simpson积分法将公式 （1） 离散表示成下述公式 （2） ： 0036 （2） 0037 其中公式 （2） 中b表示所设的相等的相邻b值间隔值， 公式 （2） 是积分的离散形 式， 可以得到数值解， 本发明运用此方法得到DWI积分数值随b值变化的信号序列。 0038 其中对步骤3)中将DWI积分的序列建立成自回归模型表示成公式： 0039 （

20、3） 0040 其中对步骤4) 中对所述的自回归模型进行最大似然估计表示成公式： 0041 （4） 0042 其中对步骤4) 中最后对ADC求解得到的解析式表示成公式： 0043 （5） 0044 TE是回波时间， 是采集图像信号时所用的参数， TE表示所设的相等的相邻b值对 应的回波时间差值， 比如采集第一幅图像对应一个b值和TE值， 表示为b1和TE1; 第三幅图 像对应一个b值和TE值， 表示为b3和TE3; TE=TE3-TE1; b=b3-b1） 0045 本发明中，“DWI序列” 是磁共振扩散加权成像， DWI技术中把施加的扩散敏感梯度 说 明 书 3/4 页 5 CN 10409

21、5635 B 6 场参数称为b值或称扩散敏感系数。 通过在成像序列中施加多个不同的扩散梯度采集不同b 值的DWI图像， 采集得到的多个的DWI图像信号随b值呈现以ADC为衰减系数的E指数衰减。 0046 本发明利用自回归模型计算生成ADC磁共振图像， 通过采集多个等间隔b值DWI序 列可获取对应的多幅不同b值的磁共振图像； 对这些磁共振图像中被试体内的所有像素的 磁共振信号进行线性利用Simpson数值积分法， 即： 对信号点关于b值做多段等间隔积分， 再 建立自回归模型， 最后对自回归模型参数进行最大似然估计得到ADC。 0047 以下分步介绍本发明方法用多个且等间隔b值DWI序列采集得到多

22、幅不同b值的磁 共振图像进行ADC值计算， 生成得到ADC磁共振图像的具体实施过程。 其中， 用于计算ADC磁 共振图像数据为多个b值人脑DWI图像数据， 数据来源于西门子3.0T磁共振成像系统， 本实 施例中所采用的b值个数为16。 0048 1、 利用磁共振成像系统用多个且等间隔b值DWI序列采集被试人脑磁共振图像数 据， 获得不同b值的DWI磁共振图像。 0049 2、 利用公式 （1） 和 （2） ， 对所述多b值DWI图像信号： 从第一个b值到第三个b值开始， 每次都进行两个b值间隔长度 （相邻三个数据点） 的Simpson法则数值积分， 从而形成一个 DWI积分数值随b值变化的信号

23、序列。 0050 例如， 本实施例采集的人脑磁共振DWI图像的b值为： 0、 100、 200、 、 1500 （单位： s/ mm2） ， 共16个b值， 相邻b值间隔为100 s/mm2。 对应的人脑DWI图像中前额叶白质区域某个像 素的信号值为： 241、 215、 197、 180、 167、 150、 143、 128、 122、 112、 108、 96、 94、 88、 82和76。 代入 公式 （1） 和 （2） 计算得到14个积分数值分别为43266.7、 39433.3、 36133.3、 33266.7、 30333.3、 28333.3、 25900.0、 24266.

24、7、 22600.0、 21333.3、 19533.3、 18666.7、 17600.0和 16400.0。 0051 3、 利用公式 （3） 将所述DWI信号积分的序列建立成自回归模型， 最后一个b值DWI信 号值可表示成由邻近的前两个b值DWI信号值的线性组合形式。 0052 4、 利用公式 （4） 对所述的自回归模型进行最大似然估计从而求解系数并利用ADC 的解析式公式 （5） 对所述磁共振图像中被试体内所有像素计算ADC值， 得到所述磁共振ADC 图像。 0053 例如将步骤3中的例举的成像参数以及像素信号值代入公式 （5） 计算得到的ADC的 值为844.6 mm2/s。 如图1所示为本实施例利用本发明方法得到的人体大脑的ADC磁共振图 像。 说 明 书 4/4 页 6 CN 104095635 B 7 图1 说 明 书 附 图 1/1 页 7 CN 104095635 B 8