发明名称 一种基于改进型阈值函数的小波变换图像去噪方法
摘要 本发明涉及一种基于改进型阈值函数的小波变换图像去噪方法。包括以下步骤:首先选择基小波,确定小波分解层数,对含噪图像进行多尺度小波变换;然后对每一层的细节系数分别确定一个阈值;再采用改进型阈值函数进行阈值处理;最后将低频系数和阈值处理后的高频系数进行小波重构,即可得到去噪后的图像。本发明提供的去噪方法,克服了硬阈值函数的不连续性的缺陷,同时减小了软阈值函数中存在的恒定偏差。采用本发明后能够有效的把信号和噪声区别开来,并且的在剔除噪声的同时保护图像信息边缘,改善了图像的峰值信噪比,具有更好的图像去噪效果。
申请公布号 CN103854264B 申请公布日期 2016.05.11
申请号 CN201410120642.5 申请日期 2014.03.28
申请人 中国石油大学(华东) 发明人 鄢志丹;杨春梅;郑金吾;卢洋
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 青岛联信知识产权代理事务所 37227 代理人 段秀瑛;王月玲
主权项 一种基于改进型阈值函数的小波变换图像去噪方法,其特征在于:所述方法包括以下步骤:步骤1:选择基小波,确定小波分解层数k,对含噪图像f(x,y)进行多尺度小波变换,分别获得不同的子带系数:最高层低频系数LL<sub>J</sub>,以及水平细节系数LH<sub>k</sub>,垂直细节系数HL<sub>k</sub>和对角细节系数HH<sub>k</sub>,k=1,2,…,J;保持LL<sub>J</sub>不变,对每一层的细节系数LH<sub>k</sub>、HL<sub>k</sub>、HH<sub>k</sub>分别确定一个阈值;步骤2:分别采用改进的阈值函数对每一层的细节系数LH<sub>k</sub>、HL<sub>k</sub>、HH<sub>k</sub>进行阈值处理,得到处理后的小波高频系数<img file="FDA0000924792740000011.GIF" wi="439" he="102" />步骤3:将最高层低频系数LL<sub>J</sub>和阈值处理后的高频系数<img file="FDA0000924792740000012.GIF" wi="407" he="103" />k=1,2,…,J进行小波重构,即可得到去噪后的图像<img file="FDA0000924792740000013.GIF" wi="176" he="94" />其中,在所述的步骤1中,阈值的确定方式如下:<maths num="0001"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><msub><mi>&sigma;</mi><mi>n</mi></msub><msqrt><mrow><mn>2</mn><msub><mi>lnN</mi><mi>k</mi></msub></mrow></msqrt><mo>*</mo><mi>m</mi></mrow>]]></math><img file="FDA0000924792740000014.GIF" wi="396" he="86" /></maths>其中:m为调节因子,当处理水平细节系数LH<sub>k</sub>和垂直细节系数HL<sub>k</sub>时,m=1;当处理对角细节系数HH<sub>k</sub>时,<img file="FDA0000924792740000015.GIF" wi="174" he="69" />N<sub>k</sub>表示第k层相应的水平、垂直、对角方向的小波系数个数;σ<sub>n</sub>为噪声标准差;在所述步骤2中,采用的改进型阈值函数数学表达式如下:<maths num="0002"><math><![CDATA[<mrow><mover><mi>w</mi><mi>&Lambda;</mi></mover><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>sgn</mi><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow><mo>(</mo><mo>|</mo><mi>w</mi><mo>|</mo><mo>-</mo><mfrac><msup><mi>T</mi><mn>2</mn></msup><mrow><mrow><mo>|</mo><mi>w</mi><mo>|</mo></mrow><mo>+</mo><msup><mi>T</mi><mn>3</mn></msup></mrow></mfrac><mo>,</mo><mo>|</mo><mi>w</mi><mo>|</mo><mo>&gt;</mo><mi>T</mi></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><msup><mrow><mo>|</mo><mi>w</mi><mo>|</mo></mrow><mn>3</mn></msup><mrow><mn>1</mn><mo>+</mo><msup><mi>w</mi><mn>2</mn></msup></mrow></mfrac><msup><mrow><mo>(</mo><mfrac><mi>w</mi><mi>T</mi></mfrac><mo>)</mo></mrow><mn>4</mn></msup><mo>,</mo><mrow><mo>|</mo><mi>w</mi><mo>|</mo></mrow><mo>&le;</mo><mi>T</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000924792740000016.GIF" wi="713" he="318" /></maths>式中,w为原始的小波细节系数,<img file="FDA0000924792740000017.GIF" wi="50" he="87" />为估计的小波细节系数,T为阈值。
地址 266000 山东省青岛市经济技术开发区长江西路66号
您可能感兴趣的专利