[发明专利]声速不均匀介质热声成像的重建算法

权利要求书

1.一种声速不均匀介质热声成像的图像重建方法，其特征在于：对声速均匀介质热声成像的时域重建法进行修正，使声波在介质中两点间的传播时间不再与距离成正比，而是一个待定值；估计出待测介质中的声速分布，通过微扰理论计算出声波在介质中两点间的传播时间，然后进行图像重建；其中：

所述的修正的时域重建法为：

A(r)=≈-η∫Ω0∫dΩ01t∂p(r0,t)∂t|t=T(r,r0)---(3)]]>

这里，A(r)是介质的电磁波吸收分布，r₀是以原点为圆心、半径为r₀的圆球面轨迹，p(r₀，t)是在圆球面上探测到的热声波，t是时间，T(r，r₀)是两点r与r₀间的声波传播时间，η是一个常数；

当待测介质的声速c(r)不均匀时，有：

u(r)＝u₀+ε·u₁(r)         (4)

这里，u(r)＝1/c(r)为介质的慢度，u₀是常数，ε是一个很小的值；由微扰理论计算T(r，r₀)的方法为：

T(r，r₀)＝T₀(r，r₀)+ε·T₁(r，r₀)+ε²·T₂(r，r₀)       (5)

T0=∫lu0ds,T1=∫lu1ds,T2=∫l12u0(u12-|▿T1|2)ds---(6)]]>

其中l是r和r₀间的线段；

估计介质的慢度分布u_e(r)由下式估计：

u_e(r)＝T(-r/|r|·r₀，r/|r|·r₀)/2r₀                    (12)

其中计算T(-r₀，r₀)的方法为：

T(-r0,r0)=argmax(Rr0(t))---(10)]]>

Rr0(t)=∫-∞∞S(r0,τ)S(-r0,-(τ-t))dτ---(5)]]>

R_r0(t)是S(r₀，t)和S(-r₀，-t)的相关函数，S(r₀，t)定义为：

S(r0,t)=[∫0tp(r0,t)dt]·t---(8).]]>