基于MCMC采样和阈值低秩逼近的图像去噪方法技术

本发明专利技术公开了一种基于MCMC和阈值低秩逼近的图像去噪方法。本方法在去噪过程中，利用蒙特卡罗采样产生图像块，并根据直方图多个统计特征，得到满足马尔科夫链的相似度判决函数；对各类图像块簇进行奇异值分解，并依据图像相应的先验信息，对奇异值进行自适应阈值估计；基于分解后的低秩结构，借助低秩逼近算法进行图像重构，从而达到去噪目的。本方法充分利用图像非局部几何结构信息相似的特性以及低秩矩阵能较好地处理高维数据的优势，克服了以往非局部均值遍历搜索方法选块复杂度高的缺点，增强了选块的鲁棒性，并在一定程度上减少了算法的运行时间。

【技术实现步骤摘要】

本专利技术涉及一种基于MCMC采样和阈值低秩逼近的图像去噪方法，属于视频图像 数据处理的

技术介绍
图像是人类活动中最常用的信息载体，已成为人们获取外界原始信息的主要途 径。然而，复杂的监测场景、恶劣的天气条件、光照变化以及采集设备的精度等易使图像受 噪声污染严重、视觉效果模糊、目标特征难以分辨，从而对图像的后续处理产生影响，无法 满足应用环境的需求。因此，作为图像预处理中不可或缺的环节，图像去噪对后续图像特征 提取、分割、压缩等起到很大作用，能够使图像更加真实地被还原，并再现目标场景。 近年来，稀疏表示得到广泛关注，Rubinstein等人在文献"Double sparsity: Iea rning sparsed dictionaries for sparse signal approximation'，提出 了稀疏 K-SVD 字 典学习方法。该方法同时兼顾了 Wavelet、Curvelet、DCT等字典的简单性，以及K-SVD、在 线学习等方法训练获得的字典具有的灵活性，解决了 K-SVD字典结构性差、高维数据处理 难的问题。然而，此方法未利用各图像块之间存在的几何结构相似特性。 非局部均值（Non-Local Means, NLM)去噪方法是Buades等在2005年提出的，此 方法提出了非局部思想的概念，充分利用了图像内部的自相似性。该方法在保持纹理和边 缘信息方面有不错的效果，但是在图像的匀质区域会出现伪纹理的现象。同时，该方法是从 像素点到像素点逐个进行处理的，导致计算复杂度很高，且未考虑图像块的能量集中性。 Dabov等人在非局部相似的理论基础上，于2007年提出了 BM3D算法。BM3D算法 是一种空域和频域相结合的去噪方式，不仅包含了非局部去噪的思想，而且利用了变换域 滤波方法。该方法去噪后的图像具有较高的信噪比，且纹理信息保留较好，特别是对变化比 较缓慢的边缘保持效果较好，是目前最理想的图像去噪算法之一。BM3D算法利用了图像变 换的能量集中性，但相似块群的构造仅仅考虑了图像块的整体信息而并未利用图像的结构 信息。 为克服BM3D算法的不足，基于图像局部结构相似特性，一些结构聚类的图像去 噪算法被相继提出。典型代表是基于聚类的稀疏表示算法（Clustering-based Sparse Representation, CSR)，CSR结合了 NLM和BM3D的思想，通过减小稀疏表示误差（Sparse Coding Noise，SCN)以提高图像的重构质量。该方法的去噪效果优于BM3D，具有较高的信 噪比，对纹理和边缘的保持也取得了非常好的效果。然而，由于通过不断的迭代更新完备主 成分分析（Principal Component Analysis, PCA)字典和稀疏表示系数，对图像块进行稀疏 表示，导致算法的时间复杂度大大增加。
技术实现思路
本专利技术目的在于结合上述去噪方法的优势并针对其不足，提出一种基于马尔科夫 蒙特卡罗采样和自适应阈值低秩逼近的图像去噪方法，结合非局部相似的特性以及低秩逼 近处理在处理高维数据时的优势，克服了字典训练的复杂度高的问题。本专利技术旨在通过蒙 特卡罗采样产生图像块，并根据直方图多个统计特征，采用满足马尔科夫链的相似度判决 函数选择图像相似块，以充分利用图像的结构信息。此外，在低秩条件下利用块之间的相似 性对相似性矩阵进行秩最小化约束逼近，以减小相似块融合不准确的现象，有效地保留图 像的细节、纹理等结构信息，从而提高去噪效果。 为实现上述目的，本专利技术的实现方法，包括如下步骤： (1)输入一幅N行M列的含噪图像Y : ⑵设置参数： 观测图像Y e RmXn，样本点s]，s]的集合S ,，样本点数m。，空间方差〇s，期望图像 1,加性噪声N，相似度评价函数ψ ; (3)采用块直方图多个统计特征相似度判据进行马尔科夫蒙特卡罗采样，即对观 测图像Y中的每个点 S]进行迭代操作，以获得相似匹配块簇； 3a)初始化参数：虚，Sj = Φ，s j = s j。，k = 1 ; 3b)对整幅图像进行采点： 按照高斯建议性分布从位置s]D开 始，随机取一点s' jk; 3c)通过s' ]k的接受概率决定是否接受s' #为序列集Ω中的点，接受概率为： 其中，Ψ (s' ]k/S]。）是用来度量采样点s' ]k与估计点s ,。近似程度的一个客观评价 函数； 3d)本方法采用直方图多个统计特征的组合来度量图像块之间的相似度，其结构 相似性评价函数为： 3e)若Ψ (s' _jk/s_j())值越大，则表示s' _jk与s 。结构信息越相似，可通过判定来决定 是否放入序列集Ω中； 3f)若 k < 111。+1，则转到 3b)，否则转到（4); (4)基于序列集Ω，获得每一个图像块的相似性矩阵Y(s])，C](Y( S]))= {s.j0, Sjl，…，S.jk}，令 C j (Y (Sj) ) = Yj; (5)依据图像的先验信息估计阈值矢量w，对序列集Ω的每一个图像块相似矩阵 Yj进行SVD分解（U 〃 Sj, Vj) = Svd(Yj)，其中Uj是左奇异矩阵，S j是奇异值对角矩阵，V』是 右奇异矩阵；从而获得估计僅 (6)对分解后的低秩结构进行秩最小化逼近，按如下公式进行秩最小化逼近： 其中，乂^是待求的变量，I I. I 表示加权核范数，if是等式右边的最优值。该等 式的最优解是通过求解以下等式： 其中，Sw(.)表示阈值算子。因此，经秩最小化逼近后得到去噪后的样本集Λ\ ; (7)若含噪图像中的每个像素点S]都处理完毕，转到（8);否则，转到（3); (8)聚集去噪后的图像样本集JTi以形成整幅干净图像X。 本专利技术与现有的技术相比具有以下优点： 1.本专利技术结合了局部、非局部和低秩逼近技术，在去除噪声的同时更能实现和抑 制为纹理，并能很好的保留待去噪图像中的纹理细节信息。 2.本专利技术由于对搜寻区域内的图像块利用直方图统计特性评价函数，并通过马尔 科夫蒙特卡罗采样获取相似块簇，增强了相似块的准确性。 3.本专利技术由于使用相似性矩阵秩最小化的逼近问题，避免了非局部均值中权值趋 同的问题，以及字典学习稀疏表示不够准确的问题。 4.本专利技术基于先验信息的特点，针对不同的奇异值采用不同的奇异值阈值计算方 式，使得相似性矩阵秩逼近更准确，最大限度的平滑光滑区域，并保持纹理信息。 5.本专利技术由于是在空域中进行的，实现过程简单，且去噪效果良好。【附图说明】 图1为本专利技术的总体流程图。 图2为不同去噪方法对图像Road的去噪效果图。 图3为不同去噪方法对图像Rail的去噪效果图。 具体的实施方式 参照附图1，本专利技术的实现步骤如下： 步骤1，输入一幅N行M列的含噪图像Y : 步骤2,设置参数： 观测图像Y e RmXn，s]表示样本点，S ,表示s ,的集合，m。表示样本点数，σ s表示 空间方差，表示期望图像，N表示加性噪声，Ψ表示相似度评价函数； 步骤3,采用块直方图多个统计特征相似度判据进行马尔科夫蒙特卡罗采样，即对 观测图像Y中的每个点 s]进行迭代操作，以获得相似匹配块簇； 3a)初始化参数：i , Sj= Φ，s .j= s .j0, 本文档来自技高网...

【技术保护点】
一种基于马尔科夫蒙特卡罗(Markov‑chain Monte Carlo,MCMC)采样和阈值低秩逼近的图像去噪方法，包括如下步骤：(1)输入一幅含噪图像Y：Y＝X+N，其中，Y＝{Y(s),s∈S}表示含噪图像，X＝{X(s),s∈S}表示不含噪图像，N＝{N(s),s∈S}表示加性噪声，S是离散格点s的集合；(2)设置参数：观测图像Y∈Rm×n，sj表示样本点，Sj表示sj的集合，m0表示样本点数，σs表示空间方差，表示期望图像，N表示加性噪声，Ψ表示相似度评价函数；(3)采用块直方图多个统计特征相似度判据进行马尔科夫蒙特卡罗采样，即对观测图像Y中的每个点sj进行迭代操作，以获得相似匹配块簇；3a)初始化参数：Sj＝φ，sj＝sj0，k＝1；3b)对整幅图像进行采点：按照高斯建议性分布即第k个候选点必须在以第k‑1个样本点为中心点的高斯分布意义下取得，从位置sj0开始，随机取一点s'jk；3c)通过s'jk的接受概率决定是否接受s'jk为序列集Ω中的点，接受概率为：其中，Ψ(s'jk/sj0)是度量采样点s'jk与估计点sj0近似程度的客观评价函数；3d)本方法采用直方图多个统计特征的组合来度量图像块之间的相似度，其结构相似性评价函数为：上式中，ρ＝{m,σ,R}表示像素点所对应图像块的直方图多个统计特征集合，其中：均值m表示纹理平均亮度，标准差s表示纹理平均对比度，平滑度R表示纹理亮度的相对平滑度；对于灰度一致的区域，其平滑度R为1；对于灰度级的值有较大差异的区域，R为0，起到保护边缘和细节信息的作用；m、s、R都对应上式中的αi；表示s'jk、sj0点 对应图像块在直方图统计特征ρ下，统计特征分量αi对应的值；3e)若Ψ(s'jk/sj0)值越大，则表示s'jk与sj0结构信息越相似，可通过判定来决定是否放入序列集Ω中；3f)若k＜m0+1，则转到3b)，否则转到(4)；(4)基于序列集Ω，获得每一个图像块的相似性矩阵Y(sj)，cj(Y(sj))＝{sj0,sj1,…,sjk}，令cj(Y(sj))＝Yj；(5)依据图像的先验信息估计阈值矢量w，对序列集Ω的每一个图像块相似矩阵Yj进行SVD分解(Uj,Sj,Vj)＝svd(Yj)，其中Uj是左奇异矩阵，Sj是奇异值对角矩阵，Vj是右奇异矩阵；从而获得估计值(6)对分解后的低秩结构进行秩最小化逼近，按如下公式进行秩最小化逼近：其中，Xj是待求的变量，表示加权核范数，是等式右边的最优值；该等式的最优解是通过求解以下等式：Xj＝Sw(Yj)＝UjSw(Σj)VjT其中，Sw(.)表示阈值算子，因此，经秩最小化逼近后得到去噪后的样本集(7)若含噪图像中的每个像素点sj都处理完毕，转到(8)；否则，转到(3)；(8)聚集去噪后的图像样本集以形成整幅干净图像X。...

【技术特征摘要】

【专利技术属性】
技术研发人员：罗晖，徐瑶，王玮，王培东，
申请(专利权)人：华东交通大学，
类型：发明
国别省市：江西;36