[发明专利]准静态超声弹性成像位移计算方法和装置

权利要求书

1.一种准静态超声弹性成像位移计算方法，其特征在于，包括：

采集得到连续的多帧超声射频信号；

对于两帧超声射频信号上的任一个位移计算中心点，利用先验时移建立基于邻域窗的迭代互相关函数，进行迭代互相关计算；

获取所述迭代互相关计算的互相关可靠性指数值，判断所述互相关可靠性指数值是否小于第一可信阈值；

如果所述互相关可靠性指数值大于或等于第一可信阈值，根据所述先验时移确定所述位移计算中心点的横向时移，根据所述先验时移和所述迭代互相关计算的迭代互相关结果确定所述位移计算中心点的纵向时移；

如果所述互相关可靠性指数值小于第一可信阈值，对所述位移计算中心点使用邻域内的二维综合互相关搜索，根据邻域内的二维综合互相关系数最大时的二维综合互相关结果，得到所述位移计算中心点的纵向时移和横向时移；

根据遍历全部超声射频信号得到的所有位移计算中心点的纵向时移和横向时移，生成被测生物组织的位移场分布。

2.根据权利要求1所述的方法，其特征在于，所述对于两帧超声射频信号上的任一个位移计算中心点，利用先验时移建立基于邻域窗的迭代互相关函数包括：

记先验时移为[t_p,t_q]，建立如下的迭代互相关函数：

其中，x()和y()分别表示所述两帧超声射频信号，t_p为两帧超声射频信号在横向上的先验时移，t_q为两帧超声射频信号在纵向上的先验时移，t_w和t_h分别是第一帧超声射频信号x()在邻域窗内的横向坐标和纵向坐标，t_w-t_p和t_h-t_q分别是第二帧超声射频信号y()在邻域窗内的横向坐标和纵向坐标，[t_w0,t_h0]是所述第一帧超声射频信号x()上的所述位移计算中心点的坐标，[t_w0-t_p,t_h0-t_q]是所述第二帧超声射频信号y()上对应的位移计算中心点的坐标，[w,h]为邻域窗的宽度和高度，R(t_w0,t_h0)表示[t_w0,t_h0]的邻域窗内的数据与[t_w0-t_p,t_h0-t_q]的邻域窗内的数据的互相关函数，“*”表示共轭运算，A表示迭代互相关结果的幅度，为迭代互相关结果的相位。

3.根据权利要求2所述的方法，其特征在于，所述获取所述迭代互相关计算的互相关可靠性指数值包括：

采用公式计算所述互相关可靠性指数值Qc，其中，

4.根据权利要求2所述的方法，其特征在于，所述获取所述迭代互相关计算的互相关可靠性指数值包括：

采用以下公式计算所述互相关可靠性指数值Qc：

其中，

5.根据权利要求2至4中任一所述的方法，其特征在于，所述根据所述先验时移确定所述位移计算中心点的横向时移，根据所述先验时移和所述迭代互相关计算的迭代互相关结果确定所述位移计算中心点的纵向时移包括：

根据所述迭代互相关结果的相位计算则，所述位移计算中心点的时移为：U_h＝Ps+t_h，U_w＝t_w，

其中，U_h为纵向时移，U_w为横向时移，λ是频率为发射脉冲中心频率的超声波在所述被测生物组织中传播的波长。