[发明专利]一种适用于弹性成像的高抗噪局部频率估计方法

权利要求书

1.一种适用于弹性成像的高抗噪局部频率估计方法，其特征在于：该方法包括以下步骤：

步骤1、根据待分析信号S_A设定估计带宽B，顺序划分整个0到π频率范围为N份，其中第i个频率范围的中心频率为ρ_i,ρ_i+1/ρ_i为定值；

步骤2、将S_A分别通过由中心频率ρ_i确定的滤波器R_i(ρ)，获得N个输出结果，其中第i个输出结果为S_Ai；

步骤3、顺序等分0～360度为M个方向范围，其中第k个方向范围内的方向为将每个S_Ai分别输入M个表达式为的方向滤波器，得到M个输出结果；第k个输出结果为其中是空间向量场与第k个方向滤波器方向的夹角；

步骤4、对每个S_Ai对应的输出结果组合计算后得到中心频率为ρ_i时的空域数据共计得到N个空域数据f_Ai；

步骤5、取相邻的ρ_i+1和ρ_i为一组，利用步骤4计算得到的N个f_Ai计算各组窄带估计结果，其中第h个窄带估计结果记为ω_h，h＝i＝1,2,......N-1；

步骤6、由步骤5所得的窄带估计结果ω_h计算宽带估计结果W，γ是用来调节加权系数(||f_Aif_A(i+1)||)的参数；

步骤2所述滤波器的表达式为

ξ＝0.03*SNR+0.15-0.2*B，且ξ∈[-1,0)∪(0,1]，SNR是待分析波动图像的信噪比；

ω_h的表达式为：或者

其中FT^-1为傅立叶逆变换，F_A(i+1)和F_Ai分别对应空域信号f_A(i+1)和f_Ai的频域表示，α＝B²β/8，β＝1/log₂(ρ_i+1/ρ_i)。