[发明专利]超声弹性成像方法、装置、设备及存储介质

权利要求书

1.一种超声弹性成像方法，其特征在于，包括：

获取目标组织区域形变前与形变后的相邻稀疏信号，所述相邻稀疏信号包括第一稀疏回波信号和与所述第一稀疏回波信号相邻的第二稀疏回波信号，所述第一稀疏回波信号为所述目标组织区域形变前的稀疏信号，所述第二稀疏回波信号为所述目标组织区域形变后的稀疏信号；

对所述第一稀疏回波信号和所述第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值；

对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，确定所述第一稀疏回波信号和所述第二稀疏回波信号的线性图像数据，并根据所述线性图像数据计算第二形变位移值；

对所述第一形变位移值和所述第二形变位移值进行求和运算，得到所述目标组织区域的最终形变位移值，并根据所述最终形变位移值生成所述目标组织区域的目标应变弹性图像；

其中，所述对所述第一稀疏回波信号和所述第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值，包括：

基于预设的零相位迭代算法，对所述第一稀疏回波信号和所述第二稀疏回波信号在同一位置纵向点上进行相位差估算，得到所述第一稀疏回波信号与所述第二稀疏回波信号之间的第一相位差；

基于所述零相位迭代算法，对所述第一稀疏回波信号和所述第二稀疏回波信号在与所述纵向点相邻的多个横向点进行相位差估算，得到多个横向点对应的第二相位差；

对所述第一相位差和多个所述第二相位差进行均值运算，得到所述第一形变位移值；

所述零相位迭代算法的计算公式为：

其中，D_pre(i+u,p)为上一帧信号相对于当前帧信号的形变位移值，Z1_pre(i+u+k,p+h)为所述第一稀疏回波信号，Z1_pos(i+u+k,p+h)为所述第二稀疏回波信号，i为横向方向的信号线号，p为纵向方向的信号点号，u为横向内的计算次数，j为虚数因子，K为Mask窗口的横向尺寸，H为Mask窗口的纵向尺寸，w₀为角频率，cOnj为共轭复数；

所述对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，确定所述第一稀疏回波信号和所述第二稀疏回波信号的线性图像数据，并根据所述线性图像数据计算第二形变位移值，包括：

对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，得到所述第一稀疏回波信号的第一高位宽模值和所述第二稀疏回波信号的第二高位宽模值；

基于预设的数据动态范围，对所述第一高位宽模值和所述第二高位宽模值进行压缩到低位宽，得到所述第一稀疏回波信号对应的第一线性图像数据和所述第二稀疏回波信号对应的第二线性图像数据；

根据所述第一线性图像数据和所述第二线性图像数据，计算所述第二形变位移值；

其中，对所述第一稀疏回波信号和第二稀疏回波信号进行求模运算的计算公式为：

Module_pre(i,p)表示第一稀疏回波信号的第一高位宽模值，I1_pre(i,p)表示第一稀疏回波信号的实部，Q1_pre(i,p)表示第一稀疏回波信号的虚部，Module_pos(i,p)表示第二稀疏回波信号的第二高位宽模值，I1_pos(i,p)表示第二稀疏回波信号的实部，Q1_pos(i,p)表示第二稀疏回波信号的虚部。

2.如权利要求1所述的超声弹性成像方法，其特征在于，所述获取目标组织区域形变前与形变后的相邻稀疏信号，包括：

获取所述目标组织区域形变前与形变后的多帧原始回波信号；

基于预设时间间隔，提取多帧所述原始回波信号中相邻的两帧原始回波信号；

对相邻的两帧所述原始回波信号进行下采样，得到所述第一稀疏回波信号和所述第二稀疏回波信号。

3.如权利要求1所述的超声弹性成像方法，其特征在于，所述根据所述第一线性图像数据和所述第二线性图像数据，计算所述第二形变位移值，包括：

基于预设光流法，对所述第一线性图像数据和所述第二线性图像数据进行运算，得到所述第二形变位移值。