发明名称 基于K-SVD和局部线性嵌套的图像降噪系统和方法
摘要 本发明公开了一种基于K-SVD和局部线性嵌套的图像降噪系统和方法,具体涉及到基于字典学的信号稀疏表示和重构及流形学的图像降噪技术。它采用K-SVD方法为框架的基于字典学的信号稀疏表示和重构技术,在求解信号稀疏表示时将局部线性嵌套作为约束条件加入目标函数,以此加强分解出的稀疏系数之间的联系,克服随机噪声对稀疏系数的影响,从而获得比原K-SVD方法更优的图像降噪效果。
申请公布号 CN102789633B 申请公布日期 2015.09.02
申请号 CN201210228638.1 申请日期 2012.07.02
申请人 河海大学常州校区 发明人 汤一彬;单鸣雷;朱昌平;韩庆邦;高远;殷澄
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 董建林
主权项 一种基于K‑SVD和局部线性嵌套的图像降噪方法,其特征在于利用以下模块:采样模块,计算拉普拉斯矩阵L模块,目标函数构造及字典、稀疏系数优化模块,估计图像块获取模块,整体估计图像块获取模块;所述方法包含以下步骤:(1)、在采样模块中输入含噪图像,对该图像采样出N个M<sub>1</sub>×M<sub>2</sub>像素大小的图像块,并记录各图像块在含噪图像中的位置,对第k个图像块像素矩阵B<sub>k</sub>按列进行顺序堆砌,k=1,2,...,N,构成一(M<sub>1</sub>M<sub>2</sub>)×1的列矢量Y<sub>k</sub>;(2)、在计算拉普拉斯矩阵L模块中对全体图像块所对应的列矢量按照局部线性嵌套方法,根据各矢量点在空间几何中的结构,计算距离的权值矩阵W,并进而构造拉普拉斯矩阵L;(3)、在目标函数构造及字典、稀疏系数优化模块中利用K‑SVD算法构造目标函数为:<img file="FDA0000743193210000011.GIF" wi="675" he="113" />约束条件为||Y<sub>k</sub>‑DX<sub>k</sub>||<sub>2</sub>≤ε进行参数D和X的优化求解,其中,Y=[Y<sub>1</sub> Y<sub>2</sub> ... Y<sub>N</sub>],D为字典,X<sub>k</sub>为Y<sub>k</sub>在D上投影的稀疏系数,k=1,2,...,N,X=[X<sub>1</sub> X<sub>2</sub> ... X<sub>N</sub>],||.||<sub>F</sub>为Frobenius范数,||.||<sub>2</sub>为2范数,Tr(.)为矩阵求迹运算,(.)<sup>T</sup>为矩阵转秩运算,β为加权系数,ε为分解误差门限;Y<sub>k</sub>为一列向量,是矩阵Y的一列;(4)、在估计图像块获取模块中利用字典D和稀疏系数X<sub>k</sub>,对列矢量Y<sub>k</sub>进行估计,其求解公式为<img file="FDA0000743193210000012.GIF" wi="231" he="85" />k=1,2,...,N,并将估计量<img file="FDA0000743193210000013.GIF" wi="59" he="84" />按行进行顺序分解,构造对应大小为M<sub>1</sub>×M<sub>2</sub>的估计图像块像素矩阵<img file="FDA0000743193210000021.GIF" wi="96" he="83" />(5)、在整体估计图像获取模块中将步骤(4)估计出的图像块,根据步骤(1)中的图像块在含噪图像中的位置信息,覆盖至含噪图像的对应位置,并记录含噪图像中每个像素点上需要覆盖的估计图像块个数w<sub>i,j</sub>,i,j为像素位置坐标,及对应的估计图像块的像素值<img file="FDA0000743193210000029.GIF" wi="95" he="89" />k=1,2,...,w<sub>i,j</sub>;(6)、利用含噪图像中每个像素点上需要覆盖的估计图像块个数w<sub>i,j</sub>,对<img file="FDA0000743193210000022.GIF" wi="74" he="91" />进行权值相加,最终获得降噪后图像对应位置的像素值<img file="FDA0000743193210000023.GIF" wi="111" he="93" /><img file="FDA0000743193210000024.GIF" wi="692" he="180" /><img file="FDA0000743193210000025.GIF" wi="77" he="92" />和α分别为含噪图像的像素值和对应的加权系数。
地址 213022 江苏省常州市新北区晋陵北路200号