一种基于低秩恢复的非局部图像去噪方法技术

本发明专利技术公开了一种基于低秩恢复的非局部图像去噪方法，包括以下步骤：图像采集并对其进行灰度变换处理；对灰度变换后的图像，通过全局寻找相似像素块建立三维相似性矩阵，然后对三维矩阵的离散余弦变换和Hadamard变换系数进行硬阈值滤波，得到去除部分噪声的相似像素块的初始估计，去除大部分噪声的同时提高低秩恢复去噪的相似像素块匹配准确度；（4）以变换域滤波作为先验知识，对初始去噪的参考像素块在搜索区域寻找相似像素块，然后利用原始图像对应位置的相似块构成相似性矩阵，对相似性矩阵进行低秩矩阵分解，从而实现有效的分离噪声和信号，得到最终去噪的图像。本发明专利技术的算法简单快速，信噪比高、一致性好，特别适用于大规模图像高质量降噪的需求。

【技术实现步骤摘要】
一种基于低秩恢复的非局部图像去噪方法
本专利技术涉及一种非局部图像噪声去除方法，属于图像处理

技术介绍
图像噪声是一种普遍存在的现象，它很大程度上是由于成像系统、传输介质以及图像在获取、传输与记录过程中不可避免地受到各种干扰造成的，容易导致图像质量降低从而阻碍图像的直观形象表达，甚至直接影响到后面的特征提取和目标识别等关键技术的成败。所以对图像进行去噪处理是很必要的。一方面减小噪声，同时保护边界和细节，对后续的图像分析有很好的效果。研究表明绝大部分旳噪声都可以用高斯白噪声、脉冲噪声或者两种噪声的混合噪声来表示。目前，图像去噪的方法主要包括局部方法(1971年Tukey提出的中值滤波；1982年Serra提出的数学形态学滤波；1992年Rudin提出的总变分、变换域滤波等)、非局部方法(非局部均值；2006年Aharon提出的基于字典学习等；2009年Wright提出的低秩恢复去噪)等。局部方法主要利用像素点的邻域来估计像素点的真实值，这种方法因为只考虑到局部的信息，对于单个像素的相似性度量易受噪声影响，很难准确地表征像素间的实际相似性，所以容易造成边缘模糊、细节丢失，去噪效果不够理想。Buades2005提出了一种与局部方法截然不同的图像去噪策略，即非局部均值图像去噪方法。这一方法充分利用了图像中存在自相似性的性质，通过加权平均一些相似像素能较好地估计像素的真实值，取得了前所未有的图像去噪效果。非局部思想通过全局寻找一些相似像素块能较好地估计参考像素块的真实值，能很好地保留纹理细节信息，并且视觉效果更好。但是目前全局寻找相似像素块的匹配准确率不高，并且没有考虑整体的相关性以及局部信息，容易引入虚假信息。近年来，基于低秩矩阵恢复的去噪模型因其考虑了整体的相关性逐渐成为图像处理领域的研究热点。主要思想是从较大的但稀疏的误差中恢复出本质上低秩的数据矩阵，一副图像如果没有噪声的话是低秩的，但噪声的加入会增加矩阵的秩。所以对含噪声的图像进行低秩恢复可以将原来的数据分解为两个矩阵，一个是低秩恢复的去噪图像(由于元素间存在某种内在的联系)，另一个是稀疏噪声。传统的低秩表示利用奇异值硬阈值实现核范数最小化，没有考虑像素块相似度的差异，容易产生条纹失真，去噪效果不够理想。
技术实现思路
针对上述现有技术存在的问题，本专利技术的目的是借助非局部的方法，提供一种基于低秩恢复的非局部图像去噪方法，以适应图像的大规模高质量降噪的需求，达到高质量图像去噪的目的。本专利技术通过以下具体步骤实现：(1)图像采集：获取图像；(2)预处理：对采集的图像进行灰度变换；(3)变换域滤波：对灰度变换后的图像进行初始的变换域去噪，去除大部分的噪声的同时，提高后续去噪过程的块匹配准确度；a)参数设置：图像灰度化后，设定参考像素块和相似像素块的搜索区域尺寸，假定当前参考像素块为P，设定在搜索区域内寻找P的相似像素块Q的数目；b)相似像素块Q的样本集获取：通过块匹配方法在以P为中心的搜索区域内寻找相似像素块Q，相似像素块集合可以表示为：c(P)＝{Q:d(P,Q)≤λ1}其中，d表示利用像素灰度特征计算两个像素块之间的归一化欧式距离，λ1是两个像素块假定相似的距离阈值，小于规定阈值的则认为搜索块和参考像素块相似，划为同一个样本集；c)三维相似性矩阵构成：将寻找的相似像素块集合c(P)按照相似度大小排序，显然第一个像素块是参考块P本身，然后聚集构成三维相似性矩阵，其中每一层代表一个相似性像素块；d)联合滤波：对三维相似性矩阵进行三维变换，其中包含进行一次二维离散余弦变换和一次一维Hadamard变换。设定阈值γ，对变换后的系数矩阵进行硬阈值滤波，只保留大于阈值的对应系数，然后反变换得到去噪后的三维相似性矩阵；因为图像的大部分真实信息都集中在三维相似性矩阵能量顶点处，而噪声往往集中在三维相似性矩阵能量的底部，因此通过变换域的硬阈值滤波可以在保证能量不损耗的情况下，有效去除噪声的同时保留绝大部分图像信息。e)初始估计：当完成了三维变换域滤波后，得到每个相似像素块和相应的每个像素点的估计值，通过计算重叠像素块的加权平均估计值来得到初始去噪后的图像；(4)低秩恢复：以变换域滤波为先验知识，进一步对初始去噪后的图像进行低秩矩阵分解，实现有效分离噪声和信号，得到最终去噪的图像。实验证明通过低秩恢复不仅能够去除图像中的噪声，还能够很好地保持图像中的真实边缘结构与重要的纹理细节，无论从客观的性能指标还是从主观的视觉感受来评价都优于同类方法。a)相似像素块的样本集获取：与初始去噪搜寻相似像素块的思想一致，通过块匹配方法在以P为中心的搜索区域内寻找相似像素块Q，相似像素块集合可以表示为：c(P)＝{Q:d(P,Q)≤λ2}其中，d表示利用像素灰度特征计算两个像素块之间的归一化欧式距离，λ2是两个像素块假定相似的距离阈值，小于规定阈值的则认为搜索块和参考像素块相似，划为同一个样本集；b)相似性矩阵构成：与初始去噪不同的是利用原始图像对应位置的相似像素块Q'代替Q构成相似性矩阵，对于每一个相似像素块Q'的列，按照从左到右的顺序首尾连接，得到相似像素块的列向量，将匹配的相似像素块对应的列向量按照相似度大小排列，构成相似性矩阵s(P)；s(P)＝o(P)+n(P)，其中o(P)是待恢复的相似性矩阵，而n(P)是噪声矩阵；c)低秩恢复：假设干净的图像是低秩的，而噪声的加入会增加矩阵的秩，所以利用加权核范数最小化从相似性矩阵s(P)中恢复出干净的相似性矩阵o(P)。具体步骤如下：对相似性矩阵进行奇异值分解，根据奇异值估计相应的权值，并根据权值对奇异值收缩，最后根据奇异值更新相似性矩阵，得到去噪后的相似性矩阵，即从s(P)中恢复得到o(P)；d)最终估计：将去噪后的相似块矩阵对噪声图像中同一个像素点的多个估计值进行加权平均，获得去噪后的每个像素点的最终估计值。e)遍历图像中的所有像素点，重复上述步骤，得到最终去噪后的图像。由以上操作步骤，即可实现图像的高质量去噪，将噪声与干净图像分开，从而提高图像的后续处理效果。本专利技术基于非局部去噪技术，通过变换域去噪和低秩恢复方法，达到高质量图像去噪效果。本专利技术不仅能够去除图像中的噪声，还能够很好地保持图像中的真实边缘结构与重要的纹理细节，无论从客观的性能指标还是从主观的视觉感受来评价都优于同类方法。由于本专利技术的算法简单快速，信噪比高、一致性好，特别适用于大规模图像高质量降噪的需求。附图说明图1为本专利技术的整体流程图。图2为本专利技术具体实施例的去噪结果图。其中，a)水下原始图像；b)水下噪声图像；c)去噪效果图像。图3为本专利技术具体实施例中不同方法的去噪效果图像。其中，a)原始水下噪声图像b)三维块匹配BM3D去噪效果图、c)非局部均值NLM去噪效果图、d)高斯尺度混合GSM去噪效果图、e)字典学习K-SVD去噪效果图、f)Ours五种去噪方法。图4为本专利技术具体实施例中不同去噪方法的对比结果图。具体实施方式下面以水下光学图像的图像去噪为例，结合附图详细说明本专利技术的具体实施过程。实验图像由水下机器人(AutomaticUnderwaterVehicle，AUV)在海底采集。本专利技术的整体流程如图1所示，具体详细过程如下：(1)水下图像获取：通过AUV在本文档来自技高网...

【技术保护点】
一种基于低秩恢复的非局部图像去噪方法，其特征在于它包括以下具体步骤：(1)图像采集；(2)对采集的图像进行灰度变换处理；(3)对灰度变换后的图像，通过全局寻找相似像素块建立三维相似性矩阵，然后对三维矩阵的离散余弦变换和Hadamard变换系数进行硬阈值滤波，得到去除部分噪声的相似像素块的初始估计，去除大部分噪声的同时提高低秩恢复去噪的相似像素块匹配准确度；(4)以变换域滤波作为先验知识，对初始去噪的参考像素块在搜索区域寻找相似像素块，然后利用原始图像对应位置的相似块构成相似性矩阵，对相似性矩阵进行低秩矩阵分解，从而实现有效的分离噪声和信号，得到最终去噪的图像。

【技术特征摘要】
1.一种基于低秩恢复的非局部图像去噪方法，其特征在于它包括以下具体步骤：(1)图像采集；(2)对采集的图像进行灰度变换处理；(3)对灰度变换后的图像，通过全局寻找相似像素块建立三维相似性矩阵，然后对三维矩阵的离散余弦变换和Hadamard变换系数进行硬阈值滤波，得到去除部分噪声的相似像素块的初始估计，去除大部分噪声的同时提高低秩恢复去噪的相似像素块匹配准确度；(4)以变换域滤波作为先验知识，对初始去噪的参考像素块在搜索区域寻找相似像素块，然后利用原始图像对应位置的相似块构成相似性矩阵，对相似性矩阵进行低秩矩阵分解，从而实现有效的分离噪声和信号，得到最终去噪的图像。2.根据权利要求1所述的非局部图像去噪方法，其特征在于步骤(3)中所述的全局寻找相似像素块的方法是利用快速的块匹配，满足c(P)＝{Q:d(P,Q)≤λ}；其中d表示利用像素灰度特征计算两个像素块之间的归一化欧式距离，λ是两个像素块假定相似的距离阈值，P是参考像素块，Q是搜索区域内任意像素块，c(P)是相似性像素块样本集。3.根据权利要求1所述的非局部图像去噪方法，其特征在于步骤(3)中...

【专利技术属性】
技术研发人员：年睿，刘雪飞，王志远，肖玫，
申请(专利权)人：中国海洋大学，
类型：发明
国别省市：山东,37