[發明專利]聲速不均勻介質熱聲成像的重建算法

無效

• 申請號: 200810032506.5
• 申請日: 2008-01-10
• 公開（公告）號: CN101214156A
• 公開（公告）日: 2008-07-09
• 發明（設計）人: 汪源源;張弛
• 申請（專利權）人: 復旦大學
• 主分類號: A61B8/00
• 分類號: A61B8/00;A61B5/00
• 代理公司: 上海正旦專利代理有限公司
• 代理人: 陸飛;盛志范

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).]]>