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

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

【技术实现步骤摘要】

本专利技术涉及核磁共振成像，尤其是涉及一种基于快速体积分方程和磁共振的人体电磁特性反演方法。
技术介绍
1991年Haacke第一次提出了基于核磁共振的电特性成像(EPT)。电特性成像与核磁共振一样都是非侵入式的成像方法，它利用核磁共振的数据反演出组织的电特性参数分布。磁共振检测伴随着热效应，这可能给人体带来伤害，尤其是未来要发展的高场磁共振检测，热效应与人体的电特性参数分布有关，了解电特性参数的分布可以提前评估出磁共振检测给人体带来的热效应及其影响，这对于核磁共振成像的发展有很大的意义。癌症由于其高死亡率和低治愈率而受到人们的广泛关注，有研究表明，癌症细胞比正常细胞的电特性参数高很多，比如乳癌细胞比正常细胞大200％，而膀胱癌细胞比正常细胞大100％，利用电特性成像就可以更清晰明了地分辨出肿瘤和正常组织，这有利于癌症的早期发现和及时理疗。目前，基于核磁共振的电特性成像方法已经有多种，比如正交鸟笼线圈法、多通道传输接收法、局部麦克斯韦成像法等等。但是这些方法都有一定的局限性，它们都基于麦克斯韦等式，大多都假设电特性参数的分布是局部均匀的，没有将电磁场边界条件考虑在内，虽然简化了计算，但使得电特性参数在组织的边界处会变得不可靠，且对于噪声也更加敏感。Balidemaj等人虽提出了利用核磁共振系统测量的数据进行迭代反演的对比源反演方法，有效地抑制边界出现的错误，但作者只重建了二维模型，且没有将射频屏蔽考虑在内。21世纪初定量磁化率成像(QSM)被提出，利用MRI方法获取磁化率定量图像，对磁化率做定量分析，可以为生物医学研究提供一种非创伤性手段，而且对组织内部顺磁性铁含量的测量，有利于对脑部血管疾病和神经病变疾病的诊断和治疗。QSM需要对相位信息进行解缠绕及去除背景场的预处理来获得反映局部磁场变化的场图，再结合重建算法重建出磁化率图像。对于QSM方法，在病态逆问题的求解情况下，背景场去除效果不理想的话会对磁化率的求解产生较大影响，甚至会导致最终的磁化率成像结果不可用，目前两种去除背景场较好的方法是复杂谐波伪影去除法(SHARP)和偶极场投影法(PDF)。由场图信息重建磁化率图像
是一个不适定逆问题，稳定求解和精确定量是关键问题。常用的QSM重建方法有多方向采样磁化率计算方法(COSMOS)、贝叶斯正则化方法、k空间加权微分法(WKD)等。由于成像仪器腔体限制和病人舒适角度问题，很难获得多个不同摆放方向的脑部成像，这极大限制了COSMOS方法的临床应用。贝叶斯正则化方法，随着匹配噪声参数增大，磁化率重建图像清晰度变好，显示出很好的去噪能力。WKD方法成像能显示出脑组织局部细节，且没有出现明显的条状伪影，但成像参数(如场强、回波时间、翻转角等)会对相位处理和磁化率分布图的精度产生影响。
技术实现思路
本专利技术的目的在于针对现有技术中使用的电特性成像(EPT)算法在组织边界处值不可靠，对于噪声更加敏感，导致成像分辨率较差，以及磁共振(MR)成像中幅值信息和相位信息分开成像，且利用相位信息的磁化率成像需要解缠绕和去除背景场处理等问题，提供一种基于快速体积分方程和磁共振的人体电磁特性反演方法。本专利技术包括以下步骤：1)测量实验数据，利用磁共振的影像技术测量人体，得到人体的B1+场，所述B1+场是磁共振射频场的正旋场；2)进行迭代计算，具体步骤为：①假定初始的电特性对比度χE；(0)(r)和磁化率χH；(0)(r),利用方程(1)和(2)，进行正演计算，解出E(0)(r)和H(0)(r)；所述方程(1)和(2)如下：Einc(r)=E(r)-[kb2+▿▿·]∫Dinvg(r,r′)χE(r′)E(r′)dr′+jωμb▿×∫Dinvg(r,r′)χH(r′)H(r′)dr′---(1)]]>Hinc(r)=H(r)-[kb2+▿▿·]∫Dinvg(r,r′)χH(r′)H(r′)dr′+jωϵb*▿×∫Dinvg(r,r′)χE(r′)E(r′)dr′---(2)]]>其中，Dinv代表反演区域，r和r’为反演区域中的位置，kb为背景介质的波数，ω为角频率，j为虚数符号；εb*，μb分别是背景介质的复介电常数和磁导率；Einc(r)和E(r)分别为入射电场和总电场，Hinc(r)和H(r)分别为入射磁场和总磁场；g(r,r’)是标量格林函，χE(r)
是电特性对比度，χH(r)是磁化率，定义为以下方程：g(r,r′)=e-jkb|r-r′|4π|r-r′|---(3)]]>χE(r)=ϵ*(r)-ϵb*ϵb*---(4)]]>χH(r)=μ(r)-μbμb---(5)]]>其中，方程(4)和(5)中的ε*(r)和μ(r)分别是人体复介电常数和磁导率，而复介电常数又可以用方程(6)表示：ϵ*(r)=ϵ(r)-jσ(r)ω---(6)]]>其中，ε(r)，σ(r)分别是人体的介电常数，电导率；②将第n-1次迭代的总电场E(r)和总磁场H(r)，代入方程(7)和(8)，进行反演计算，解出第n次迭代的χE；(n)(r)和χH；(n)(r)；所述方程(7)和(8)如下：Hsca(r)=jωϵb*▿×∫Dinvg(r,r′)χE(r′)E(r′)dr′+[kb2+▿▿·]∫Dinvg(r,r′)χH(r′)H(r′)dr′---(7)]]>H1+;sca(r)=B1+(r)μ(r)-B1+;inc(r)μb---(8)]]>其中，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)(r)，χH；(n)(r)代入方程(7)和(8)算出的B1+(r)与测量值的相对残差在本文档来自技高网...

【技术保护点】
基于快速体积分方程和磁共振的人体电磁特性反演方法，其特征在于包括以下步骤：1)测量实验数据，利用磁共振的影像技术测量人体，得到人体的B1+场，所述B1+场是磁共振射频场的正旋场；2)进行迭代计算，具体步骤为：①假定初始的电特性对比度χE；(0)(r)和磁化率χH；(0)(r),利用方程(1)和(2)，进行正演计算，解出E(0)(r)和H(0)(r)；所述方程(1)和(2)如下：Einc(r)=E(r)-[kb2+▿▿.]∫Dinvg(r,r′)χE(r′)E(r′)dr′+jωμb▿×∫Dinvg(r,r′)χH(r′)H(r′)dr′---(1)]]>Hinc(r)=H(r)-[kb2+▿▿.]∫Dinvg(r,r′)χH(r′)H(r′)dr′-jωϵb*▿×∫Dinvg(r,r′)χE(r′)E(r′)dr′---(2)]]>其中，Dinv代表反演区域，r和r’为反演区域中的位置，kb为背景介质的波数，ω为角频率，j为虚数符号；εb*，μb分别是背景介质的复介电常数和磁导率；Einc(r)和E(r)分别为入射电场和总电场，Hinc(r)和H(r)分别为入射磁场和总磁场；g(r,r’)是标量格林函，χE(r)是电特性对比度，χH(r)是磁化率，定义为以下方程：g(r,r′)=e-jkb|r-r′|4π|r-r′|---(3)]]>χE(r)ϵ*(r)-ϵb*ϵb*---(4)]]>χH(r)=μ(r)-μbμb---(5)]]>其中，方程(4)和(5)中的ε*(r)和μ(r)分别是人体复介电常数和磁导率，而复介电常数又可用方程(6)表示：ϵ*(r)=ϵ(r)-jσ(r)ω---(6)]]>其中，ε(r)，σ(r)分别是人体的介电常数，电导率；②将第n‑1次迭代的总电场E(r)和总磁场H(r)，代入方程(7)和(8)，进行反演计算，解出第n次迭代的χE；(n)(r)和χH；(n)(r)；所述方程(7)和(8)如下：Hsca(r)=jωϵb*▿×∫Dinvg(r,r′)χE(r′)E(r′)dr′+[kb2+▿▿.]∫Dinvg(r,r′)χH(r′)H(r′)dr′---(7)]]>H1+;sca(r)=B1+(r)μ(r)-B1+;inc(r)μb---(8)]]>其中，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)(r)，χH；(n)(r)代入方程(7)和(8)算出的B1+(r)与测量值的相对残差在10％以内，则算法收敛，结束迭代；所述迭代的算法采用玻恩迭代算法、变分玻恩迭代算法或变形玻恩迭代算法；3)利用方程(4),(5),(6)结合迭代结果中的电特性对比度χE(r)和磁化率χH(r),计算出ε(r)，σ(r)，μ(r)(或磁化率χH(r))的分布，输出结果，计算完毕。...

【技术特征摘要】
1.基于快速体积分方程和磁共振的人体电磁特性反演方法，其特征在于包括以下步骤：1)测量实验数据，利用磁共振的影像技术测量人体，得到人体的B1+场，所述B1+场是磁共振射频场的正旋场；2)进行迭代计算，具体步骤为：①假定初始的电特性对比度χE；(0)(r)和磁化率χH；(0)(r),利用方程(1)和(2)，进行正演计算，解出E(0)(r)和H(0)(r)；所述方程(1)和(2)如下：Einc(r)=E(r)-[kb2+▿▿.]∫Dinvg(r,r′)χE(r′)E(r′)dr′+jωμb▿×∫Dinvg(r,r′)χH(r′)H(r′)dr′---(1)]]>Hinc(r)=H(r)-[kb2+▿▿.]∫Dinvg(r,r′)χH(r′)H(r′)dr′-jωϵb*▿×∫Dinvg(r,r′)χE(r′)E(r′)dr′---(2)]]>其中，Dinv代表反演区域，r和r’为反演区域中的位置，kb为背景介质的波数，ω为角频率，j为虚数符号；εb*，μb分别是背景介质的复介电常数和磁导率；Einc(r)和E(r)分别为入射电场和总电场，Hinc(r)和H(r)分别为入射磁场和总磁场；g(r,r’)是标量格林函，χE(r)是电特性对比度，χH(r)是磁化率，定义为以下方程：g(r,r′)=e-jkb|r-r′|4π|r-r′|---(3)]]>χE(r)ϵ*(r)-ϵb*ϵb*---(4)]]>χH(r)=μ(r)-μbμb---(5)]]>其中，方程(4)和(5)中的ε*(r)和μ(r)分别是人体复介电常数和磁导率，而复介电常数又可用方程(6)表示：ϵ*(r)=ϵ(r)-jσ(r)ω---(6)]]>其中，ε(r)，σ(r)分别是人体的介...

【专利技术属性】
技术研发人员：柳清伙，洪荣汉，李胜男，刘娜，张建华，
申请(专利权)人：厦门大学，
类型：发明
国别省市：福建;35