主权项 |
一种基于线性泊松方程的磁热声成像电导率重建方法,其特征在于:所述的磁热声成像电导率重建方法包括以下步骤:第一步:获取磁热声信号由具有线性调频电流激励源(1)产生线性调频电流激励信号经激励线圈作用在导电物体(4)上,导电物体(4)由于感应电流的作用产生焦耳热,进而产生电磁热声信号;利用超声换能器(6)接收产生的电磁热声信号,然后通过超声信号处理、采集子系统(7)对信号进行采集和处理;控制电路实现对电流激励源、超声换能器(6)和超声信号处理、采集子系统(7)的控制和同步;利用超声换能器(6)对产生的电磁热声信号进行断层圆周扫描,获取每个断层圆周上的电磁热声信号,将获得电磁热声信号进行存储;第二步:获取热声源分布根据电磁热声的产生原理,得出电磁热声成像的声压波动方程:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mo>▿</mo><mn>2</mn></msup><mi>p</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mn>1</mn><msup><msub><mi>c</mi><mi>s</mi></msub><mn>2</mn></msup></mfrac><mfrac><msup><mo>∂</mo><mn>2</mn></msup><msup><mrow><mo>∂</mo><mi>t</mi></mrow><mn>2</mn></msup></mfrac><mi>p</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mfrac><mi>β</mi><msub><mi>C</mi><mi>P</mi></msub></mfrac><mi>S</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><msup><mi>δ</mi><mo>′</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000633077010000011.GIF" wi="1407" he="162" /></maths>其中r为超声换能器与热声源之间的距离,p(r,t)是声压,cx为热声源在介质中的传播声速,C<sub>P</sub>为导电物体的比热容,β为导电物体的热膨胀系数,δ(t)是狄拉克函数,S(r)是热声源分布,t是时间项;根据声压波动方程,利用时间反演法得到热函数的重建公式:<img file="FDA0000633077010000012.GIF" wi="1574" he="145" />其中R=|r′‑r|,<img file="FDA0000633077010000013.GIF" wi="192" he="131" />r'是超声换能器的位置,S<sub>d</sub>是超声换能器所在的平面,p′是声压对时间的一阶导数,n是r′位置S<sub>d</sub>的法线矢量;选取导电物体的某一断层面z=z<sub>0</sub>,利用公式(2)获取断层面上的热声源S(x,y,z<sub>0</sub>);导电物体上所有热声源分布通过断层数据S(x,y,z<sub>0</sub>)在z方向的插值得到,或者通过移动超声换能器通过测量的方式获取;第三步:获取导电物体标量电位空间分量导电物体的热声源直接与导电物体的电场强度和电导率相关:S=σE<sup>2</sup>=σE·E (3)其中E是磁场强度;因生物组织的电导率较低,电场强度的空间分量E表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>E</mi><mo>≈</mo><mo>-</mo><mo>▿</mo><mi>φ</mi><mo>-</mo><msub><mi>A</mi><mn>1</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000633077010000021.GIF" wi="1075" he="77" /></maths>其中A<sub>1</sub>为一次磁矢位的空间分量,一次磁矢位与导电物体无关,是电流激励源在真空中产生的磁矢位,φ是标量电位的空间分量,<img file="FDA0000633077010000022.GIF" wi="47" he="51" />是哈密顿算符;根据电流连续性定理,有:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mo>▿</mo><mo>·</mo><mi>σ</mi><mrow><mo>(</mo><mo>▿</mo><mi>φ</mi><mo>+</mo><msub><mi>A</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000633077010000023.GIF" wi="1140" he="83" /></maths>其中,σ为导电物体的电导率,对导电物体进行空间离散,给出导电物体电导率的初值[σ]<sup>0</sup>,将[σ]<sup>0</sup>和一次磁矢位的空间分量A<sub>1</sub>代入公式(5),利用线性有限元方法求解得到标量电位的空间分量[φ]<sup>1</sup>;第四步:定义目标函数由公式(3)和公式(4)得到:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>σ</mi><mo>=</mo><mfrac><mi>S</mi><mrow><mi>E</mi><mo>·</mo><mi>E</mi></mrow></mfrac><mo>=</mo><mfrac><mi>S</mi><mrow><mrow><mo>(</mo><mo>▿</mo><mi>φ</mi><mo>+</mo><msub><mi>A</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>·</mo><mrow><mo>(</mo><mo>▿</mo><mi>φ</mi><mo>+</mo><msub><mi>A</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000633077010000024.GIF" wi="1254" he="146" /></maths>定义目标函数:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>σ</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><mi>S</mi><mo>-</mo><mi>σ</mi><mrow><mo>(</mo><mo>▿</mo><mi>φ</mi><mo>+</mo><msub><mi>A</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>·</mo><mrow><mo>(</mo><mo>▿</mo><mi>φ</mi><mo>+</mo><msub><mi>A</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000633077010000025.GIF" wi="1264" he="101" /></maths>第五步:最优化迭代求解导电物体电导率:将标量电位的空间分量[φ]<sup>1</sup>和一次磁矢位的空间分量A<sub>1</sub>代入目标函数,通过最优化算法进行迭代,直到误差满足终止条件,所得的电导率即为最终导电物体的最终电导率。 |