[发明专利]基于小波变换模极大值法的磁共振信号工频噪声消减方法

权利要求书

1.一种基于小波变换模极大值法的磁共振信号工频噪声消减方法，其特征在于，包括以下步骤：

a、向理想的磁共振观测信号中加入工频噪声，构建含工频噪声的全波磁共振信号，并对含噪信号进行n层小波分解，n取7；

b、提取出各层细节系数d₁，…，d_n以及最后一层近似系数a_n；

c、用小波变换模极大值去噪算法更新每一层细节系数；

d、利用重建的小波系数重构信号；

步骤c包括：

1)求出每一尺度上的小波系数模极大值点及其相应的模极大值；

2)在最大尺度上设定阈值，更新最大尺度上信号产生的模极大值点；

3)根据尺度J的模极大值点及其邻域范围，更新尺度J-1上的模极大值点；

4)在尺度1上保留与尺度2位置相对应的模极大值点的位置；

5)利用各个尺度上剩余的模极大值点进行小波系数的重建；

在最大尺度上设定阈值，更新最大尺度上信号产生的模极大值点包括:

对最高层细节系数的模极大值进行硬阈值处理，当最高层细节系数的模极大值大于阈值λ时，保留其对应的模极大值点，反之剔除其对应的模极大值点，更新最大尺度上细节系数的模极大值点；

根据尺度J的模极大值点及其邻域范围，更新尺度J-1上的模极大值点包括：利用最大尺度上的模极大值点位置构造搜索邻域U(x_Jn)，其中x_Jn表示尺度J上的第n个模极大值点，在邻域U(x_Jn)内搜寻尺度J-1上的细节系数模极大值点，保留尺度J-1上位于邻域U(x_Jn)内的模极大值点，清除其余模极大值点，从而更新尺度J-1上的细节系数模极大值点；

在尺度1上保留与尺度2位置相对应的模极大值点的位置包括：将其他位置的模极大值点置为0；

求出每一尺度上的小波系数模极大值点及其相应的模极大值包括：

求取工频噪声的均方误差σ；

将工频噪声的均方差σ带入式(1)求取阈值λ，式(1)其中λ是实数，N是最高层细节系数的长度。

2.按照权利要求1所述的方法，其特征在于，步骤3)结束后令J＝J-1，判断J是否大于2，若J2，循环步骤3)，直到尺度数为2时，停止循环。