[发明专利]一种剪切波波速计算方法及弹性成像设备

摘要

本发明主要提出一种剪切波波速计算方法及弹性成像设备，该方法先通过对不同检测位置的“形变‑时间”曲线进行频谱分析，根据有效频率范围对“形变‑时间”曲线进行滤波、插值。利用特定窗函数与插值后的“形变‑时间”曲线相乘增加互相关匹配运算的精度。最后对不同检测位置的“形变‑时间”曲线进行互相关匹配运算；根据检测位置距离和偏移大小计算出剪切波速度。该弹性成像设备中具有这样的计算剪切波波速的模块。这种方法通过对检测位置的“形变‑时间”曲线进行预处理，增加“形变‑时间”曲线的信噪比，然后利用互相关匹配方法寻找不同检测位置的“形变‑时间”最大偏移位置，从而计算出剪切波速度值。

主权项

1.一种剪切波波速计算方法，其特征在于：包括以下步骤：步骤1、确定待观察深度位置的步骤；该步骤中选取两个横向检测位置，根据形变估计结果矩阵作出所选的横向检测位置对应得“形变‑时间”曲线；步骤2、确定有效频率范围；设计低通滤波器的步骤；该步骤中对各横向检测位置对应得“形变‑时间”的曲线进行FFT频谱分析，选取有效频率范围，并设计FIR低通滤波器系数；步骤3、滤波的步骤；该步骤中使用步骤2中设计的低通滤波器对各个检测位置的“形变‑时间”曲线进行低通滤波；步骤4、插值的步骤；对滤波后的“形变‑时间”曲线进行样条插值；步骤5、对插值后的“形变‑时间”曲线加窗的步骤；该步骤中，选取窗函数类型、参数，根据“形变‑时间”曲线长度设计窗函数系数，并与“形变‑时间”曲线相乘；步骤6、相关匹配运算步骤；该步骤中，对各位置加窗后的“形变‑时间”曲线进行互相关匹配运算；步骤7、计算剪切波速度值的步骤；该步骤中根据互相关匹配结果找到最大偏移值，结合检测位置间隔计算出剪切波速度值。