[发明专利]磁声耦合成像的微弱声信号频域检测处理方法

权利要求书

1.一种磁声耦合成像的微弱声信号频域检测处理方法，其特征在于，在磁声耦合成像系统中进行如下步骤：

(1)确定传感器(5)空间运动初始角度θ_A和终止角度θ_B，以及空间的步进角θ_step；

(2)由扫描驱动装置(6)驱动传感器(5)至空间运动初始角度θ_A位置处；

(3)根据激励信号，以及磁声耦合成像系统函数和传感器频响特征，分析确定待测信号的频段分布特征，待测磁声信号频段上下限ω_a和ω_b，及设定频率间隔ω_step；

(4)设锁定频率ω_i＝ω_a，对待测磁声信号进行锁定放大，测得锁定频率ω_i下的幅值A_i，相位φ_i，并存储；

(5)判断当前锁定频率ω_i：如果ω_i＜ω_b，则令ω_i＝ω_a+ω_step，再重复第(4)步骤继续进行在锁定频率ω_i下锁定放大求得幅值相位；否则转入第(6)步骤；

(6)调出所有频率对应的幅值A和相位φ信息；

(7)利用傅立叶逆变换，求得时域磁声信号波形y_j(t)；

(8)判断当前传感器(5)空间运动角度θ_j是否为终止角度位置，若未达到终止角度位置，则返回步骤(3)；若达到终止角度位置时停止扫描，转入第(9)步骤。

(9)调用存储的各个测量位置的时域波形，代入图像重建算法；

(10)重建组织电导率分布特性图像。

2.根据权利要求1所述的磁声耦合成像的微弱声信号频域检测处理方法，其特征在于，进行步骤(2)的同时，由激励源(4)对待成像目标(1)进行激励，产生磁声耦合信号。

3.根据权利要求1所述的磁声耦合成像的微弱声信号频域检测处理方法，其特征在于，步骤(3)所述的分析是采用傅立叶变换进行的。

4.根据权利要求1所述的磁声耦合成像的微弱声信号频域检测处理方法，其特征在于，步骤(8)所述的判断当前传感器(5)空间运动角度θ_j是否为终止位置是判断：是否为θ_j＜θ_B。

5.根据权利要求1所述的磁声耦合成像的微弱声信号频域检测处理方法，其特征在于，步骤(7)还包括存储域磁声信号波形y_j(t)。

6.根据权利要求1所述的磁声耦合成像的微弱声信号频域检测处理方法，其特征在于，步骤(9)所述的图像重建算法是：再利用傅里叶逆变换，结合格林函数法，建立三维柱坐标系，利用洛伦兹力散度法，推出声源重建公式为

▿·(J×B0)=-12πc3∫∫S0dS01-(z0-z)2|r-r0|21t×∂2p(r0,t)∂t2|t=|r-ro|/c]]>

其中，π为圆周率，c为真空中电磁波传播速度，S₀为积分曲面，r为待成像目标空间矢量，r₀为观测点位置，z₀为观测点柱坐标中z方向分量，z为待成像目标空间矢量z方向分量，t为时间，B₀为静磁场，。

根据麦克斯韦方程组，以及电磁场边界条件，进一步推出待测组织的电导率

σ(r)=|▿·(J×B0)||B1×B0|]]>

其中，B₁为交变磁场，

利用上述声源重建公式和待测组织的电导率公式即求得待成像组织电导率分布图像。