基于快速体积分方程和磁共振的人体电磁特性反演方法.pdf

《基于快速体积分方程和磁共振的人体电磁特性反演方法.pdf》由会员分享，可在线阅读，更多相关《基于快速体积分方程和磁共振的人体电磁特性反演方法.pdf（13页珍藏版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利 (10)授权公告号 (45)授权公告日 (21)申请号 201610190456.8 (22)申请日 2016.03.30 (65)同一申请的已公布的文献号 申请公布号 CN 105877747 A (43)申请公布日 2016.08.24 (73)专利权人 厦门大学 地址 361005 福建省厦门市思明南路422号 (72)发明人 柳清伙 洪荣汉 李胜男 刘娜 张建华 (74)专利代理机构 厦门南强之路专利事务所 (普通合伙) 35200 代理人 马应森 (51)Int.Cl. A61B 5/055(2006.01) 审查员 陈林杰 (。

2、54)发明名称 基于快速体积分方程和磁共振的人体电磁 特性反演方法 (57)摘要 基于快速体积分方程和磁共振的人体电磁 特性反演方法， 涉及核磁共振成像。 结合玻恩迭 代算法或变分玻恩迭代算法或变形玻恩迭代算 法、 快速傅里叶变换、 稳定双共轭梯度算法以及 共轭梯度算法， 进行反复的正演和反演迭代至结 果收敛， 求解人体电磁特性参数， 实现人体电磁 特性参数磁共振成像的方法， 得到人体组织电磁 特性参数分布的二维或三维图像， 呈现组织内各 区域电磁特性参数差异的同时， 也可对其进行定 量研究， 可用于医学上的疾病研究和指导临床诊 断、 治疗。 权利要求书2页 说明书7页 附图3页 CN 105。

3、877747 B 2019.02.01 CN 105877747 B 1.基于快速体积分方程和磁共振的人体电磁特性反演方法， 其特征在于包括以下步 骤： 1)测量实验数据， 利用磁共振的影像技术测量人体， 得到人体的B1+场， 所述B1+场是磁共 振射频场的正旋场； 2)进行迭代计算， 具体步骤为： 假定初始的电特性对比度 E； (0)(r)和磁化率 H； (0)(r),利用方程(1)和(2)， 进行正演 计算， 解出E(0)(r)和H(0)(r)； 所述方程(1)和(2)如下： 其中， Dinv代表反演区域， r和r 为反演区域中的位置， kb为背景介质的波数， 为角频 率， j为虚数符号；。

4、b*，b分别是背景介质的复介电常数和磁导率； Einc(r)和E(r)分别为入射 电场和总电场， Hinc(r)和H(r)分别为入射磁场和总磁场； g(r,r )是标量格林函， E(r)是电 特性对比度， H(r)是磁化率， 定义为以下方程： 其中， 方程(4)和(5)中的 *(r)和 (r)分别是人体复介电常数和磁导率， 而复介电常数 又可用方程(6)表示： 其中，(r)， (r)分别是人体的介电常数， 电导率； 将第n-1次迭代的总电场E(r)和总磁场H(r)， 代入方程(7)和(8)， 进行反演计算， 解 出第n次迭代的 E； (n)(r)和 H； (n)(r)； 所述方程(7)和(8)。

5、如下： 权 利 要 求 书 1/2 页 2 CN 105877747 B 2 其中， Hsca(r)散射磁场； H1+； sca(r)为正旋散射磁场， 定义为H1+； sca(r)Hxsca(r)+jHysca (r)/2； Hxsca(r)和Hysca(r)分别为散射磁场Hsca(r)的x分量与y分量； B1+(r)为磁共振正旋场 的磁通密度， B1+； inc(r)为背景介质下磁共振正旋场的磁通密度； 将步骤解出的 E； (n)(r)和 H； (n)(r)， 代入方程(1)和(2)， 进行正演计算， 解出第n次 迭代的E(n)(r)和H(n)(r)； 重复步骤和， 直到利用 E； (n)(。

6、r)， H； (n)(r)代入方程(7)和(8)算出的B1+(r)与测 量值的相对残差在10以内， 则算法收敛， 结束迭代； 所述迭代的算法采用玻恩迭代算法、 变分玻恩迭代算法或变形玻恩迭代算法； 3)利用方程(4),(5),(6)结合迭代结果中的电特性对比度 E(r)和磁化率 H(r),计算出 (r)， (r)，(r)(或磁化率 H(r)的分布， 输出结果， 计算完毕。 2.如权利要求1所述基于快速体积分方程和磁共振的人体电磁特性反演方法， 其特征 在于在步骤2)第和步骤中， 所述正演计算为己知电特性对比度 E(r)和磁化率 H(r)， 利 用方程(1)和(2)， 结合稳定双共轭梯度算法和快。

7、速傅里叶变换， 计算出总电场E(r)和总磁 场H(r)。 3.如权利要求1所述基于快速体积分方程和磁共振的人体电磁特性反演方法， 其特征 在于在步骤2)第步骤中， 所述反演计算为己知总电场E(r)和总磁场H(r)， 利用方程(7)和 (8)， 结合共轭梯度算法和快速傅里叶变换， 计算出电特性对比度 E(r)和磁化率 H(r)。 权 利 要 求 书 2/2 页 3 CN 105877747 B 3 基于快速体积分方程和磁共振的人体电磁特性反演方法 技术领域 0001 本发明涉及核磁共振成像， 尤其是涉及一种基于快速体积分方程和磁共振的人体 电磁特性反演方法。 背景技术 0002 1991年Haa。

8、cke第一次提出了基于核磁共振的电特性成像(EPT)。 电特性成像与核 磁共振一样都是非侵入式的成像方法， 它利用核磁共振的数据反演出组织的电特性参数分 布。 磁共振检测伴随着热效应， 这可能给人体带来伤害， 尤其是未来要发展的高场磁共振检 测， 热效应与人体的电特性参数分布有关， 了解电特性参数的分布可以提前评估出磁共振 检测给人体带来的热效应及其影响， 这对于核磁共振成像的发展有很大的意义。 癌症由于 其高死亡率和低治愈率而受到人们的广泛关注， 有研究表明， 癌症细胞比正常细胞的电特 性参数高很多， 比如乳癌细胞比正常细胞大200， 而膀胱癌细胞比正常细胞大100， 利用 电特性成像就可以。

9、更清晰明了地分辨出肿瘤和正常组织， 这有利于癌症的早期发现和及时 理疗。 目前， 基于核磁共振的电特性成像方法已经有多种， 比如正交鸟笼线圈法、 多通道传 输接收法、 局部麦克斯韦成像法等等。 但是这些方法都有一定的局限性， 它们都基于麦克斯 韦等式， 大多都假设电特性参数的分布是局部均匀的， 没有将电磁场边界条件考虑在内， 虽 然简化了计算， 但使得电特性参数在组织的边界处会变得不可靠， 且对于噪声也更加敏感。 Balidemaj等人虽提出了利用核磁共振系统测量的数据进行迭代反演的对比源反演方法， 有效地抑制边界出现的错误， 但作者只重建了二维模型， 且没有将射频屏蔽考虑在内。 0003 2。

10、1世纪初定量磁化率成像(QSM)被提出， 利用MRI方法获取磁化率定量图像， 对磁 化率做定量分析， 可以为生物医学研究提供一种非创伤性手段， 而且对组织内部顺磁性铁 含量的测量， 有利于对脑部血管疾病和神经病变疾病的诊断和治疗。 QSM需要对相位信息进 行解缠绕及去除背景场的预处理来获得反映局部磁场变化的场图， 再结合重建算法重建出 磁化率图像。 对于QSM方法， 在病态逆问题的求解情况下， 背景场去除效果不理想的话会对 磁化率的求解产生较大影响， 甚至会导致最终的磁化率成像结果不可用， 目前两种去除背 景场较好的方法是复杂谐波伪影去除法(SHARP)和偶极场投影法(PDF)。 由场图信息重。

11、建磁 化率图像是一个不适定逆问题， 稳定求解和精确定量是关键问题。 常用的QSM重建方法有多 方向采样磁化率计算方法(COSMOS)、 贝叶斯正则化方法、 k空间加权微分法(WKD)等。 由于成 像仪器腔体限制和病人舒适角度问题， 很难获得多个不同摆放方向的脑部成像， 这极大限 制了COSMOS方法的临床应用。 贝叶斯正则化方法， 随着匹配噪声参数增大， 磁化率重建图像 清晰度变好， 显示出很好的去噪能力。 WKD方法成像能显示出脑组织局部细节， 且没有出现 明显的条状伪影， 但成像参数(如场强、 回波时间、 翻转角等)会对相位处理和磁化率分布图 的精度产生影响。 发明内容 0004 本发明的。

12、目的在于针对现有技术中使用的电特性成像(EPT)算法在组织边界处值 说 明 书 1/7 页 4 CN 105877747 B 4 不可靠， 对于噪声更加敏感， 导致成像分辨率较差， 以及磁共振(MR)成像中幅值信息和相位 信息分开成像， 且利用相位信息的磁化率成像需要解缠绕和去除背景场处理等问题， 提供 一种基于快速体积分方程和磁共振的人体电磁特性反演方法。 0005 本发明包括以下步骤： 0006 1)测量实验数据， 利用磁共振的影像技术测量人体， 得到人体的B1+场， 所述B1+场 是磁共振射频场的正旋场； 0007 2)进行迭代计算， 具体步骤为： 0008 假定初始的电特性对比度 E；。

13、 (0)(r)和磁化率 H； (0)(r),利用方程(1)和(2)， 进行 正演计算， 解出E(0)(r)和H(0)(r)； 0009 所述方程(1)和(2)如下： 0010 0011 0012 其中， Dinv代表反演区域， r和r 为反演区域中的位置， kb为背景介质的波数， 为 角频率， j为虚数符号；b*，b分别是背景介质的复介电常数和磁导率； Einc(r)和E(r)分别为 入射电场和总电场， Hinc(r)和H(r)分别为入射磁场和总磁场； g(r,r )是标量格林函， E(r) 是电特性对比度， H(r)是磁化率， 定义为以下方程： 0013 0014 0015 0016 其中，。

14、 方程(4)和(5)中的 *(r)和 (r)分别是人体复介电常数和磁导率， 而复介电 常数又可以用方程(6)表示： 0017 0018 其中，(r)， (r)分别是人体的介电常数， 电导率； 0019 将第n-1次迭代的总电场E(r)和总磁场H(r)， 代入方程(7)和(8)， 进行反演计 算， 解出第n次迭代的 E； (n)(r)和 H； (n)(r)； 0020 所述方程(7)和(8)如下： 0021 说 明 书 2/7 页 5 CN 105877747 B 5 0022 0023 其中， Hsca(r)散射磁场； H1+； sca(r)为正旋散射磁场， 定义为H1+； sca(r)Hxs。

15、ca(r)+ jHysca(r)/2； Hxsca(r)和Hysca(r)分别为散射磁场Hsca(r)的x分量与y分量； B1+(r)为磁共振正 旋场的磁通密度， B1+； inc(r)为背景介质下磁共振正旋场的磁通密度； 0024 将步骤解出的 E； (n)(r)和 H； (n)(r)， 代入方程(1)和(2)， 进行正演计算， 解出 第n次迭代的E(n)(r)和H(n)(r)； 0025 重复步骤和， 直到利用 E； (n)(r)， H； (n)(r)代入方程(7)和(8)算出的B1+(r) 与测量值的相对残差在10以内， 则算法收敛， 结束迭代。 0026 所述迭代的算法可采用玻恩迭代算。

16、法(BIM)、 变分玻恩迭代算法(VBIM)或变形玻 恩迭代算法(DBIM)等。 0027 在步骤2)第和步骤中， 所述正演计算为己知电特性对比度 E(r)和磁化率 H (r)， 利用方程(1)和(2)， 结合稳定双共轭梯度算法(BCGS)和快速傅里叶变换(FFT)， 计算出 总电场E(r)和总磁场H(r)。 0028 在步骤2)第步骤中， 所述反演计算为己知总电场E(r)和总磁场H(r)， 利用方程 (7)和(8)， 结合共轭梯度算法(CG)和快速傅里叶变换(FFT)， 计算出电特性对比度 E(r)和 磁化率 H(r)。 0029 3)利用方程(4),(5),(6)结合迭代结果中的电特性对比。

17、度 E(r)和磁化率 H(r),计 算出 (r)， (r)，(r)(或磁化率 H(r)的分布， 输出结果， 计算完毕。 0030 本发明具体涉及电和磁体积分方程， 结合玻恩迭代算法(BIM)或变分玻恩迭代算 法(VBIM)或变形玻恩迭代算法(DBIM)、 快速傅里叶变换(FFT)、 稳定双共轭梯度算法(BCGS) 以及共轭梯度算法(CG)， 进行反复的正演和反演迭代至结果收敛， 求解人体电磁特性参数 (包括电导率 ， 介电常数 和磁导率 )， 实现人体电磁特性参数磁共振成像的方法， 得到人 体组织电磁特性参数分布的二维或三维图像， 呈现组织内各区域电磁特性参数差异的同 时， 也可对其进行定量研。

18、究， 可用于医学上的疾病研究和指导临床诊断、 治疗。 0031 本发明能有效抑制边界出现的错误、 减小噪声影响、 提高成像分辨率， 结合MR信号 的幅值信息和相位信息， 利用快速傅里叶变换(FFT)加快重建速度、 减少重建时间， 结合玻 恩迭代算法(BIM)或变分玻恩迭代算法(VBIM)或变形玻恩迭代算法(DBIM)， 稳定双共轭梯 度算法(BCGS)， 共轭梯度算法(CG),同时反演电特性参数(电导率 ， 介电常数 )和磁特性参 数(磁导率 )。 0032 本发明的有益技术效果如下： 0033 提出一种基于快速体积分方程算法和磁共振的人体电磁特性反演方法， 对于电特 性成像， 能够有效抑制边。

19、界出现的错误、 减小噪声影响、 提高成像分辨率， 利用快速傅里叶 变换(FFT)加快重建速度、 减少重建时间， 结合变分玻恩迭代算法(BIM)或变分玻恩迭代算 法(VBIM)或变形玻恩迭代算法(DBIM)， 稳定双共轭梯度算法(BCGS)， 共轭梯度算法(CG)， 同 时反演电特性参数(电导率 、 介电常数 )和磁特性参数(磁导率 )。 附图说明 0034 图1是反演出的人脑相对介电常数三维分布。 说 明 书 3/7 页 6 CN 105877747 B 6 0035 图2是反演出的人脑电导率三维分布。 0036 图3是反演出的人脑磁化率三维分布。 0037 图4是反演出的人脑相对介电常数分布。

20、X-Y截面。 0038 图5是反演出的人脑电导率分布X-Y截面。 0039 图6是反演出的人脑磁化率分布X-Y截面。 具体实施方式 0040 本发明利用快速体积分方程算法和磁共振数据进行人体的二维或三维电磁特性 反演， 其中快速体积分方程算法中结合了玻恩迭代算法(BIM)或变分玻恩迭代算法(VBIM) 或变形玻恩迭代算法(DBIM)， 快速傅里叶变换(FFT)， 稳定双共轭梯度算法(BCGS)， 共轭梯 度算法(CG)。 0041 本实施例以变分玻恩迭代算法(VBIM)为例子进行解释。 具体实施方式如下： 0042 1)测量实验数据。 利用磁共振的影像技术测量人体， 得到B1+场。 B1+场是。

21、磁共振射 频场的正旋场。 0043 2)进行迭代计算， 包括以下步骤： 0044 假定初始的电特性对比度 E； (0)(r)和磁化率 H； (0)(r),进行正演计算， 解出E(0) (r)和H(0)(r)。 0045 利用第n-1次迭代的总电场E(r)和总磁场H(r)， 进行反演计算， 解出第n次迭代 的 E； (n)(r)和 H； (n)(r)。 0046 再由步骤解出的 E； (n)(r)和 H； (n)(r)， 进行正演计算， 解出第n次迭代的E(n) (r)和H(n)(r)。 0047 重复步骤和， 直到 E； (n)(r)， H； (n)(r)代入方程(14)和(15)算出的B1+。

22、(r)与 测量值的相对残差在10以内， 则算法收敛， 结束迭代。 0048 步骤和所述正演计算为己知电特性对比度 E(r)和磁化率 H(r)， 利用以下方 程计算总电场E(r)和总磁场H(r): 0049 0050 0051 0052 0053 其中， Dinv代表反演区域， r和r 为反演区域中的位置， kb为背景介质的波数， 为 角频率， j为虚数符号。 b*，b分别是背景介质的复介电常数和磁导率。 Einc(r)和E(r)分别为 入射电场和总电场， Hinc(r)和H(r)分别为入射磁场和总磁场。 A(r)和F(r)分别为磁势矢量 和电势矢量。 g(r,r )是标量格林函数， E(r)是。

23、电特性对比度， H(r)是磁化率， 定义为以下 方程： 说 明 书 4/7 页 7 CN 105877747 B 7 0054 0055 0056 0057 其中， 方程(4)和(5)中的 *(r)和 (r)分别是人体复介电常数和磁导率， 而复介电 常数又可以用方程(6)来表示： 0058 0059 其中，(r)， (r)分别是人体的介电常数， 电导率。 0060 便于计算， 利用冲激函数作为基函数和测试函数， 则方程(1)， (2)， (3)和(4)被离 散化为下列方程： 0061 0062 0063 0064 0065 其中， 离散的单元格中心为： 0066 0067 j1,J,k1,K,。

24、l1,L 0068 x、 y和z分别为单元格在x、 y和z方向的大小。 j， k， l分别为离散后x方向， y方 向和z方向的位置。 J， K， L为x方向， y方向和z方向的离散点数。 0069 而方程(11)和(12)中的V为单元格体积， 定义为Vxyz。 0070 由于， 标量格林函数g(r,r )具有平移不变性， 所以方程(11)和(12)可以利用快速 傅里叶变换(FFT)进行快速计算。 此方法可以节省大量时间以及内存， 使计算更加快速。 0071 离散化后， 可以将方程写成以下形式： 0072 0073 其中， 0074 0075 L代表的是线性操作。 说 明 书 5/7 页 8 C。

25、N 105877747 B 8 0076 直接解方程(13)会花大量时间和内存， 所以在本发明当中利用稳定双共轭梯度算 法(BCGS)和快速傅里叶变换(FFT)计算方程(13)得出总电场E(r)和总磁场H(r)。 0077 步骤所述反演计算为己知总电场E(r)和总磁场H(r)， 利用以下方程计算电特性 对比度 E(r)和磁化率 H(r): 0078 0079 0080 其中， Hsca(r)散射磁场。 H1+； sca(r)为正旋散射磁场， 定义为H1+； sca(r)Hxsca(r)+ jHysca(r)/2， 其中， Hxsca(r)和Hysca(r)分别为散射磁场Hsca(r)的x分量与。

26、y分量。 B1+(r)为磁 共振正旋场的磁通密度， B1+； inc(r)为背景介质下磁共振正旋场的磁通密度。 GbE(r,r )和GbH (r,r )分别为电并矢格林函数和磁并格林函数。 0081 本发明中利用变分玻恩迭代来计算物体的电磁特性分布。 所以方程(14)写成以下 形式： 0082 0083 其中， 0084 Hsca(r)Hsca(r)-Hsca； (n-1)(r) (17) 0085 代表的是测量值与n-1次迭代的散射磁场的差值。 0086 E(r) E； (n)(r)- E； (n-1)(r) (18) 0087 H(r) H； (n)(r)- H； (n-1)(r) (19。

27、) 0088 代表的是第n次迭代与第n-1次迭代计算的电特性对比度 E(r)和磁化率 H(r)的差 值。 0089 为了便于计算， 利用冲激函数作为基函数和测试函数， 则方程(16)写为以下形式： 0090 0091 其中， 电并矢格林函数GbE(r,r )和磁并格林函数GbH(r,r )都具有平移不变性， 所 以也可以利用快速傅里叶变换(FFT)来进行快速计算。 0092 当利用方程(20)计算 E(r)和 H(r)时， 如果直接用矩阵方式进行计算的话， 花费 的时间与内存过于庞大， 所以本发明利用共轭梯度算法(CG)， 结合快速傅里叶变换来计算 E(r)和 H(r)， 然后利用方程(18)。

28、和(19)算出 E； (n)(r)和 H； (n)(r)。 0093 3)利用方程(6),(7),(8)结合迭代结果中的电特性对比度 E(r)和磁化率 H(r),计 算出 (r)， (r)，(r)(或磁化率 H(r)的分布， 输出结果， 计算完毕。 0094 表1 0095 组织相对介电常数电导率(S/m)磁化率 脑干980.5110 说 明 书 6/7 页 9 CN 105877747 B 9 小脑1170.7190 灰质980.511610-8 白质680.291310-8 头骨170.0600 下丘脑980.511210-8 眼睛851.000 舌头750.6520 0096 人脑设置参数如表1,频率为63.87MHz。 说 明 书 7/7 页 10 CN 105877747 B 10 图1 图2 图3 说 明 书 附 图 1/3 页 11 CN 105877747 B 11 图4 图5 说 明 书 附 图 2/3 页 12 CN 105877747 B 12 图6 说 明 书 附 图 3/3 页 13 CN 105877747 B 13 。