[发明专利]一种基于对数正态滤波器的局部频率估计方法

说明书

本发明公开了一种基于对数正态滤波器的局部频率估计方法，涉及数字信号处理技术领域，本发明中提出了新的窄带估计结果的表达式：和新的计算宽带估计结果公式原有技术只能在α＝1的情况下进行窄带估计，使得当滤波器带宽设定后，用来覆盖0～π频段的滤波器组中心频率比值关系也就不能改变，使得滤波器估计带宽和滤波器个数相互制约，本发明从理论上突破了这一限制，使得可以保持滤波器估计带宽不变的情况下，采用更多滤波器个数，从而可以最大限度的发挥对数正态滤波器的优越性，从而实现了同时提高基于对数正态滤波器的局部频率估计算法的灵敏度和抗噪声性能的有益效果。

技术领域

本发明主要涉及数字信号处理技术在核磁共振弹性成像领域的应用，具体为一种基于对数正态滤波器的局部频率估计方法。

背景技术

在临床上医生采用手工触诊判定组织病变历史非常悠久，组织穿刺也是一种获取组织特性的重要手段，但是其非常不利于病人的健康。核磁共振弹性成像是一种非常有用的无损组织弹性测量技术，但是目前主要还是科研使用，没有形成正式临床诊断标准及规范。核磁共振弹性成像需要通过在人体组织输入机械波，让组织产生振动，通过核磁共振弹性成像专用相位敏感序列采集机械波在人体组织中的位移分布，然后通过核磁共振弹性成像反演算法获得组织弹性。众所周知，机械波的衰减随着频率的升高而变得迅速，意味着频率越高，能传入身体组织的深度越浅。为了获得完整的组织特性，我们需要尽可能低的机械波驱动频率，然而为了获得更精确更多细节的信息，需要用更高的机械波驱动频率。然而这是一对矛盾，不可能同时满足。其中重要的一个环节就在于核磁共振弹性成像的反演，反演算法目前最成熟常用的算法有两种，局部频率估计算法(LFE)和基于波动方程的代数反演算法(AIDE)，都还存在居多问题。特别是在低频率驱动的情况下，AIDE算法虽然探测灵敏度高，但是存在反演伪影严重、抗噪声性能差等问题，LFE算法虽然反演伪影少，但对噪声干扰也非常敏感，并且小肿瘤探测灵敏度不高。为了进一步提高核磁共振弹性成像反演算法性能，本发明从基于对数正态滤波器组的局部频率估计算法入手，深入分析频率估计灵敏度影响因素以及宽带频率估计误差产生的原因，从而提出新的见解。

发明内容

本发明的目的在于克服现有技术的不足，提供了一种基于对数正态滤波器的局部频率估计方法，以解决现有技术中AIDE算法虽然探测灵敏度高，但是存在反演伪影严重、抗噪声性能差等问题，LFE算法虽然反演伪影少，但对噪声干扰也非常敏感，并且小肿瘤探测灵敏度不高等技术问题。

发明了一种基于对数正态滤波器的局部频率估计方法，该方法包括以下步骤：

步骤1、获得需要分析局部频率性质的待分析信号S_A，根据灵敏度需求设定估计带宽B，将整个0到π频率范围用N个中心频率顺序分割，第i个中心频率为ρ_i，i＝1,2,..........,N，且使得相邻两个频率范围内的中心频率ρ_i+1/ρ_i为定值，确定参数，β＝1/log₂(ρ_i+1/ρ_i)，β为相邻中心频率的对数倒数比例关系系数；

步骤2、将S_A分别通过由中心频率ρ_i确定的传递函数为的对数正态滤波器，获得N个输出结果，其中第i个输出结果为S_Ai；

步骤3、将0～360度方向顺序等分为M块，其中第k个块的方向表示为且该方向为表达式为的方向滤波器的方向，将步骤2所得到的每一个S_Ai分别送给M个表达式为的方向滤波器，得到M个输出结果；其中第k个输出结果为k＝1,2,3,...,M，其中是空间向量场与第k个方向滤波器方向的夹角；

步骤4、将步骤3中每个S_Ai对应所得到的M个输出结果组合计算，得到中心频率为ρ_i时的输出结果