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

摘要

本发明公开了一种适用于弹性成像的高抗噪局部频率估计算法，涉及数字信号处理技术领域，本发明中提出了一种适用于弹性成像的高抗噪局部频率估计算法，其特点在于将传统对数正态滤波器乘以频率指数加权的系数，获得新表达式ξ＝0.03*SNR+0.15‑0.2*B，ξ∈[‑1,0)∪(0,1]；并提出一种新的窄带局部频率估计算法或者替换先前的局部频率估计算法；本发明中新的窄带估计表达式在相同噪声条件下可以进一步提高小肿瘤的辨识能力，也就是可以进一步提高灵敏度；实验表明本发明提出的新算法可以实现0dB极端恶劣条件下的弹性成像，这使得低场磁共振弹性成像成为可能；在高性噪比情况下可以实现更高的灵敏度，可以提高高场条件下高质量高分辨率的成像效果。

主权项

1.一种适用于弹性成像的高抗噪局部频率估计算法，该方法包括以下步骤：步骤1、根据待分析信号SA设定估计带宽B，顺序划分整个0到π频率范围为N份，其中第i个频率范围的中心频率为ρi,ρi+1/ρi为定值步骤2、将SA分别通过由中心频率ρi确定的滤波器Ri(ρ)，获得N个输出结果，其中第i个输出结果为SAi步骤3、顺序等分0～360度为M个方向范围，其中第k个方向范围内的方向为将每个S_Ai分别输入M个表达式为的方向滤波器，得到M个输出结果；第k个输出结果为k＝1,2,3,...,M，其中是空间向量场u与第k个方向滤波器方向的夹角步骤4、对每个S_Ai对应的输出结果组合计算后得到中心频率为ρ_i时的空域数据共计得到N个空域数据f_Ai步骤5、取相邻的ρi+1和ρi为一组，利用步骤4计算得到的N个fAi计算各组窄带估计结果，其中第h个窄带估计结果记为ωh，h＝i＝1,2,......N‑1步骤6、由步骤5所得的窄带估计结果ωh计算宽带估计结果W其特征在于：步骤2所述滤波器的表达式为ξ＝0.03*SNR+0.15‑0.2*B，且ξ∈[‑1,0)∪(0,1]，SNR是待分析波动图像的信噪比。