[发明专利]结合回波时移技术和2D热传导方程的超声波温度成像方法

权利要求书

1.一种结合回波时移技术和2D热传导方程的超声波温度成像方法，用于对离体的猪肝组织进行成像，包括组织升温阶段和组织降温阶段：

在组织升温阶段：

A.在烧灼之前，获取超声波射频信号r(x,z,t₀)，此时的超声波声速取c₀＝1540m/s；

B.开始烧灼；

C.在时间点t＝t_i，获取超声波射频信号r(x,z,t_i)，对两个时间点t_i-1和t_i的射频信号做一维互相关分析，得到射频信号在两个时间点的时间偏移累加的回波信号时间偏移图通过下式获得

D.为了降低温度估计中出现的纹波效应，使用8阶低通巴特沃斯滤波器Butterworth对获得的回波信号时间偏移图进行轴向和侧向滤波；

E.根据下列公式：

其中

其中，T(z)为温度，c₀(z)为依赖于组织深度z的超声波声速，为时间偏移对轴向的微分，α(z)为热膨胀系数，β(z)代表超声波声速变化与组织温度变化之间的线性关系；将D步骤的计算结果对时间偏移图进行轴向微分，再乘上系数k和即可得到组织温度分布图，大小为m×n；

在组织降温阶段：

A.首先获得烧灼停止时刻的原始射频信号，记为r(x,z,t_h)；

B.使用基于回波时移的超声波温度成像技术获得时刻t_h的组织温度分布图；

C.根据以下公式：

其中代表位置(x_i,y_j)在时刻t_h的组织温度，l＝α²(t_h+1-t_h)，为热扩散系数，q为热传导系数，c为介质比热容，ρ为介质密度，取q＝0.586,单位：W/m·K,对于位置为(x_i,y_j)的像素点在时刻t_h+1的温度可以通过该点以及该点周围的点在时刻t_h的温度u^h估计出：设计一个滑动窗使用该滑动窗扫描整个图像矩阵，估计时刻t_h+1的温度图像中每一个像素点的温度值，另外，由于图像边缘的像素点使用滑动窗的时候会超出图像矩阵的索引范围，使用室内温度值T₀对图像矩阵进行扩充，使图像的大小变为(m+2)×(n+2)；

D.经过组织降温阶段步骤C即可得到时刻t_h+1的组织温度分布；然后对该图像重复组织降温阶段步骤B即可得到下一时刻的组织温度分布图像。