基于非下采样轮廓波变换的图像去噪方法及装置.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN106204495A(43)申请公布日2016.12.07(21)申请号201610566766.5(22)申请日2016.07.18(71)申请人北京邮电大学地址100876北京市海淀区西土城路10号(72)发明人杜军平梁美玉刘红刚(74)专利代理机构北京超凡志成知识产权代理事务所(普通合伙)11371代理人朱文杰(51)Int.Cl.G06T5/00(2006.01)权利要求书2页说明书8页附图3页(54)发明名称基于非下采样轮廓波变换的图像去噪方法及装置(57)摘要本发明提供了一种基于非下采样轮廓波变换的图像去噪方法及装置，所述方法包括：利用非下采样轮廓波变换将待去噪的图像进行多尺度分解，获得每个尺度的低频分量的第一系数以及每个尺度的高频分量在各个方向上的第二系数；根据每个尺度的第二系数中与该尺度的上一尺度的第一系数的相关性，对各第二系数进行阈值收缩处理或保留，获得各尺度的去噪系数；根据所述去噪系数进行非下采样轮廓波逆变换，获得去噪后的图像。CN106204495ACN106204495A权利要求书1/2页1.一种基于非下采样轮廓波变换的图像去噪方法，其特征在于，所述方法包括：利用非下采样轮廓波变换将待去噪的图像进行多尺度分解，获得每个尺度的低频分量的第一系数以及每个尺度的高频分量在各个方向上的第二系数；根据每个尺度的第二系数中与该尺度的上一尺度的第一系数的相关性，对各第二系数进行阈值收缩处理或保留，获得各尺度的去噪系数；根据所述去噪系数进行非下采样轮廓波逆变换，获得去噪后的图像。2.根据权利要求1所述的方法，其特征在于，所述对各第二系数进行阈值收缩处理或保留包括：当所述第二系数为第一尺度的系数，对该第二系数进行阈值收缩处理，获得该第二系数对应的去噪系数，所述第一尺度为所述图像进行非下采样轮廓波变换后的第一个尺度；当所述第二系数不为第一尺度的系数且该第二系数的父系数为噪声信号，对该第二系数进行阈值收缩处理，获得该第二系数对应的去噪系数，所述父系数为该第二系数所在尺度的上一尺度的第一系数；当所述第二系数不为第一尺度的系数且该第二系数的父系数为有效信号，保留该第二系数作为去噪系数。3.根据权利要求2所述的方法，其特征在于，所述方法还包括：计算所述第二系数对应的自适应性阈值；当所述第二系数的父系数大于该第二系数对应的自适应阈值，判定该第二系数的父系数为有效信号，当所述第二系数的父系数小于该第二系数对应的自适应阈值，判定该第二系数的父系数为噪声信号。4.根据权利要求3所述的方法，其特征在于，所述计算所述第二系数对应的自适应性阈值包括：通过公式进行计算，其中，w1(j，d，p)为在尺度j的方向d上的该第二系数p的值，Thr(j，d，p)为该第二系数p对应的自适应性阈值，Num为该第二系数局部邻域窗口中包含的第二系数的数量，w1(j，d，l)是以w1(j，d，p)为中心的局部邻域窗口中的第二系数，是噪声方差估计函数，是位于尺度j的各个方向的第二系数，τ为预定值。5.根据权利要求1所述的方法，其特征在于，对第二系数进行阈值收缩处理包括：通过最小化无偏风险估计SURE算法计算线性参数ad,1以及ad,2，构造与该第二系数相关的线性函数θd(ω)＝ad，1t1(ω)+ad，2t2(ω)，其中，θd(ω)表示在方向d上的该第二系数对应的去噪系数，t1(ω)＝ω，表示该第二系数的值。6.一种基于非下采样轮廓波变换的图像去噪装置，其特征在于，所述装置包括：2CN106204495A权利要求书2/2页变换模块，用于利用非下采样轮廓波变换将待去噪的图像进行多尺度分解，获得每个尺度的低频分量的第一系数以及每个尺度的高频分量在各个方向上的第二系数；处理模块，用于根据每个尺度的第二系数中与该尺度的上一尺度的第一系数的相关性，对各第二系数进行阈值收缩处理或保留，获得各尺度的去噪系数；逆变换模块，用于根据所述去噪系数进行非下采样轮廓波逆变换，获得去噪后的图像。7.根据权利要求6所述的装置，其特征在于，所述处理模块还包括：收缩单元，用于当所述第二系数为第一尺度的系数，对该第二系数进行阈值收缩处理，获得该第二系数对应的去噪系数，所述第一尺度为所述图像进行非下采样轮廓波变换后的第一个尺度，或者当所述第二系数不为第一尺度的系数且该第二系数的父系数为噪声信号，对该第二系数进行阈值收缩处理，获得该第二系数对应的去噪系数，所述父系数为该第二系数所在尺度的上一尺度的第一系数；保留单元，用于当所述第二系数不为第一尺度的系数且该第二系数的父系数为有效信号，保留该第二系数作为去噪系数。8.根据权利要求7所述的装置，其特征在于，所述装置还包括：计算模块，用于计算所述第二系数对应的自适应性阈值；判定模块，用