[发明专利]一种剪切波波速计算方法及弹性成像设备在审
| 申请号： | 201811418565.6 | 申请日： | 2018-11-26 |
| 公开（公告）号： | CN109589138A | 公开（公告）日： | 2019-04-09 |
| 发明（设计）人： | 翁嘉淳;董永刚 | 申请（专利权）人： | 深圳中科乐普医疗技术有限公司 |
| 主分类号： | A61B8/08 | 分类号： | A61B8/08 |
| 代理公司： | 深圳市智科友专利商标事务所 44241 | 代理人： | 周小年 |
| 地址： | 518000 广东省深圳市宝安区石*** | 国省代码： | 广东;44 |
| 权利要求书： | 查看更多 | 说明书： | 查看更多 |
| 摘要： | |||
| 搜索关键词： | 形变 检测位置 弹性成像设备 剪切波波速 互相关 匹配运算 剪切波 预处理 相乘 有效频率范围 大小计算 频谱分析 最大偏移 窗函数 信噪比 偏移 滤波 匹配 | ||
1.一种剪切波波速计算方法，其特征在于：包括以下步骤：
步骤1、确定待观察深度位置的步骤；该步骤中选取两个横向检测位置，根据形变估计结果矩阵作出所选的横向检测位置对应得“形变-时间”曲线；
步骤2、确定有效频率范围；设计低通滤波器的步骤；该步骤中对各横向检测位置对应得“形变-时间”的曲线进行FFT频谱分析，选取有效频率范围，并设计FIR低通滤波器系数；
步骤3、滤波的步骤；该步骤中使用步骤2中设计的低通滤波器对各个检测位置的“形变-时间”曲线进行低通滤波；
步骤4、插值的步骤；对滤波后的“形变-时间”曲线进行样条插值；
步骤5、对插值后的“形变-时间”曲线加窗的步骤；该步骤中，选取窗函数类型、参数，根据“形变-时间”曲线长度设计窗函数系数，并与“形变-时间”曲线相乘；
步骤6、相关匹配运算步骤；该步骤中，对各位置加窗后的“形变-时间”曲线进行互相关匹配运算；
步骤7、计算剪切波速度值的步骤；该步骤中根据互相关匹配结果找到最大偏移值，结合检测位置间隔计算出剪切波速度值。
2.根据权利要求1所述的剪切波波速计算方法，其特征在于：所述的步骤2中，选取有效频率范围时，频率衰减在-3dB范围内的。
3.根据权利要求1所述的剪切波波速计算方法，其特征在于：所述的步骤2中，设计FIR低通滤波器系数时，采用用窗函数法设计，窗函数选择切比雪夫窗函数。
4.根据权利要求1所述的剪切波波速计算方法，其特征在于：所述的步骤5中，窗口类型选择Tukey窗函数，参数选取0.5。
5.一种弹性成像设备，包括发射声辐射冲击组织施加激励的装置，接收待观察深度位置的装置，对接收到的剪切波进行处理的装置，其特征在于：所述的对接收到的剪切波进行处理的装置中包括剪切波波速计算装置；其特征在于：所述的剪切波波速计算装置包括：
确定待观察深度位置的装置；该装置选取两个横向检测位置，根据形变估计结果矩阵作出所选的横向检测位置对应得“形变-时间”曲线；
对各个检测位置的“形变-时间”曲线进行低通滤波的低通滤波器；
插值的装置；对滤波后的“形变-时间”曲线进行样条插值；
对插值后的“形变-时间”曲线加窗的装置；该装置选取窗函数类型、参数，根据“形变-时间”曲线长度设计窗函数系数，并与“形变-时间”曲线相乘；
相关匹配运算装置；该装置对各位置加窗后的“形变-时间”曲线进行互相关匹配运算；
计算装置；计算装置根据互相关匹配结果找到最大偏移值，结合检测位置间隔计算出剪切波速度值。
6.根据权利要求5所述的弹性成像设备，其特征在于：还包括对“形变-时间”曲线进行FFT分析的FFT频谱分析装置。
7.根据权利要求6所述的弹性成像设备，其特征在于：所述的低通滤波器为FIR低通滤波器，采用用窗函数法设计，窗函数选择切比雪夫窗函数。
8.根据权利要求6所述的弹性成像设备，其特征在于：所述的相关匹配运算装置的窗口类型选择Tukey窗函数，参数选取0.5。
该专利技术资料仅供研究查看技术是否侵权等信息，商用须获得专利权人授权。该专利全部权利属于深圳中科乐普医疗技术有限公司，未经深圳中科乐普医疗技术有限公司许可，擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作，请联系【客服】
本文链接：http://www.vipzhuanli.com/pat/books/201811418565.6/1.html，转载请声明来源钻瓜专利网。
