[发明专利]一种生物磁光声联合内窥成像方法

权利要求书

1.一种生物磁光声联合内窥成像方法，其特征是，所述方法在建立多层腔体组织横截面模型的基础上，对腔道内的超声回波成像、光声成像和感应式磁声成像的过程进行数值仿真，得到腔道横截面上组织反射的超声回波信号以及组织产生的光声和磁声信号，然后对三种超声信号进行最优加权求和，得到融合后的联合成像信号。

2.根据权利要求1所述的一种生物磁光声联合内窥成像方法，其特征是，所述方法包括以下步骤：

a.建立多层腔体组织横截面模型

成像导管位于多层腔体组织横截面模型的中心，超声换能器位于成像导管顶端，模型所在的坐标系是θ-l极坐标系，其中坐标原点是成像导管顶端中心，θ是极角，l是极径，水平向右的方向为l轴正方向，忽略超声换能器的孔径效应，将超声换能器看作理想的点探测器，其扫描轨迹为平行于成像平面的圆形轨迹，以模型的中心为起始点，将模型等角度划分为N份，每一份近似为层与层之间平行的多层腔体组织，导管处的成像角度为

θ_i＝360(i-1)/N

其中i＝1,2,…,N，θ_i对应的成像区域的角度范围为[θ_ia,θ_ib]，其中θ_ia＝θ_i-180/N，θ_ib＝θ_i+180/N；

b.对成像组织反射的超声回波信号以及组织产生的光声信号和磁声信号进行数值仿真；

④超声回波信号：

从腔体中心沿径向发射超声脉冲，根据不同成分组织的声阻抗差异值以及超声探测器的冲激响应，仿真得到组织反射/散射的超声回波信号；

⑤光声信号：

从腔体中心对周围组织沿径向发射激光脉冲，组织由于光声效应产生光声信号，根据不同成分组织的光吸收系数和散射系数，结合蒙特卡罗模拟和光声波动方程，仿真得到组织产生的光声信号；

⑥磁声信号：

沿腔体的轴向施加静磁场和脉冲磁激励，组织由于磁声效应产生磁声信号，根据不同成分组织的电导率，结合磁声波动方程，仿真得到组织产生的磁声信号；

c.对超声回波信号、光声信号和磁声信号进行融合

对于成像组织中的角度θ_i(i＝1,2,…,N)、位置r处的超声回波信号P_i^(U)(r)、光声信号P_i^(P)(r)和磁声信号P_i^(M)(r)，采用下式进行融合：

其中为超声回波、光声和磁声信号的融合信号，即联合成像信号，W_1i、W_2i、W_3i分别为超声回波、光声和磁声信号的加权因子，且有W_1i+W_2i+W_3i＝1。

3.根据权利要求2所述的一种生物磁光声联合内窥成像方法，其特征是，为了使超声回波、光声和磁声信号的融合信号的均方误差最小，超声回波、光声和磁声信号的加权因子W_1i、W_2i、W_3i应取最优值最优值的计算公式为：

W1i*=1σ1i2(1σ1i2+1σ2i2+1σ3i2)W2i*=1σ2i2(1σ1i2+1σ2i2+1σ3i2)W3i*=1σ3i2(1σ1i2+1σ2i2+1σ3i2)]]>

σ1i2=RUU-RUP]]>

σ2i2=RPP-RPM]]>

σ3i2=RMM-RMU]]>

RUU=k-1kRUU(k-1)+1kPi(U)(k)Pi(U)(k)]]>

RUP=k-1kRUP(k-1)+1kPi(U)(k)Pi(P)(k)]]>

RPP=k-1kRPP(k-1)+1kPi(P)(k)Pi(P)(k)]]>

RPM=k-1kRPM(k-1)+1kPi(P)(k)Pi(M)(k)]]>

RMM=k-1kRMM(k-1)+1kPi(M)(k)Pi(M)(k)]]>

RMU=k-1kRMU(k-1)+1kPi(M)(k)Pi(U)(k)]]>

其中k为超声回波、光声和磁声信号的测量次数，为总均方误差最小时超声回波、光声和磁声信号的最优加权因子，为P_i^(U)(r)的方差，为P_i^(P)(r)的方差，为P_i^(M)(r)的方差，R_UU为P_i^(U)(r)的自相关系数，R_UP为P_i^(U)(r)和P_i^(P)(r)的互相关系数，R_PP为P_i^(P)(r)的自相关系数，R_PM为P_i^(P)(r)和P_i^(M)(r)的互相关系数，R_MM为P_i^(M)(r)的自相关系数，R_MU为P_i^(M)(r)和P_i^(U)(r)的互相关系数。