CN105678774A - 一种基于主成分分析的图像噪声水平估计方法 - Google Patents
一种基于主成分分析的图像噪声水平估计方法 Download PDFInfo
- Publication number
- CN105678774A CN105678774A CN201610016782.7A CN201610016782A CN105678774A CN 105678774 A CN105678774 A CN 105678774A CN 201610016782 A CN201610016782 A CN 201610016782A CN 105678774 A CN105678774 A CN 105678774A
- Authority
- CN
- China
- Prior art keywords
- image
- noise level
- matrix
- principal component
- component analysis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000000513 principal component analysis Methods 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims abstract description 49
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 4
- 230000017105 transposition Effects 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 description 5
- 239000013598 vector Substances 0.000 description 4
- 238000001914 filtration Methods 0.000 description 2
- NAPPWIFDUAHTRY-XYDRQXHOSA-N (8r,9s,10r,13s,14s,17r)-17-ethynyl-17-hydroxy-13-methyl-1,2,6,7,8,9,10,11,12,14,15,16-dodecahydrocyclopenta[a]phenanthren-3-one;(8r,9s,13s,14s,17r)-17-ethynyl-13-methyl-7,8,9,11,12,14,15,16-octahydro-6h-cyclopenta[a]phenanthrene-3,17-diol Chemical compound O=C1CC[C@@H]2[C@H]3CC[C@](C)([C@](CC4)(O)C#C)[C@@H]4[C@@H]3CCC2=C1.OC1=CC=C2[C@H]3CC[C@](C)([C@](CC4)(O)C#C)[C@@H]4[C@@H]3CCC2=C1 NAPPWIFDUAHTRY-XYDRQXHOSA-N 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000012847 principal component analysis method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于主成分分析的图像噪声水平估计方法,包括:步骤1,将原始图像划分为若干大小相同的图像块;步骤2,将选取的图像块组成图像数据矩阵,并求取图像数据矩阵的协方差矩阵;步骤3,对协方差矩阵做特征分解,将所得的最小特征值作为噪声水平的估计值,并利用该噪声水平的估计值计算阈值;步骤4,计算每个图像块对应的梯度矩阵,以及梯度矩阵对应的协方差矩阵的迹;步骤5,选取迹小于等于阈值的图像块,重复进行步骤2~步骤4,直至连续两次计算得到的噪声水平的估计值保持不变或迭代达到预定次数后停止。本发明利用PCA更准确地找到图像中平滑的块,具有准确性和鲁棒性,适用于大范围内噪声水平和各种场景。
Description
技术领域
本发明涉及计算机图像处理技术领域,具体涉及一种基于主成分分析的图像噪声水平估计方法。
背景技术
噪声水平是许多图像处理程序中的重要参数,比如,图像去噪,图像分割等等。当下大多数去噪算法的性能很大程度地依赖于对噪声水平的估计,而解决办法也基本上是假设噪声水平是已知的,但是显然不适用于真实的场景。
最常见的噪声模型是加性高斯白噪声,以此为例,噪声水平估计算法的目标就是在给定单张噪声图像的条件下,估计出未知的高斯噪声的标准差σn。关于这个问题,已经相继提出了许多算法。一般,可以将这些方法分为,基于滤波的方法,基于图像块的方法和基于统计的方法。
基于滤波的方法,首先使用一个高通滤波器对图像进行滤波以抑制图像结构,然后通过计算噪声图像与滤波后图像的差值数据的方差得到最终的噪声水平估计值。这种类型的方法通常假设这个差值图像就是噪声,但这个假设并不是始终成立,尤其在图像有比较复杂的结构或者说良好的细节的时候。
基于图像块的方法中,图像被分成许多个图像块。每一个图像块可以认为是图像中一个有特定大小的,比如N*N的矩形窗口。有最小方差的图像块,其亮度变化最小,而一个均匀的图像块的亮度变化主要是由噪声引起的。
最近,Zoran和Weiss提出了一种基于统计的方法来分析DCT域滤波后的图像,其峰值的变化主要由噪声引起,通过分析这些变化来计算噪声水平。
到目前为止,如何准确地估计输入图像的噪声水平值依然是一个挑战,尤其是当图像纹理比较丰富的时候,所以,一个鲁棒的噪声水平估计算法是迫切需要的。
发明内容
本发明提供了一种基于主成分分析的图像噪声水平估计方法,利用主成分分析(PCA)的方法更准确地找到图像中平滑的块,具有准确性和鲁棒性,适用于大范围内噪声水平和各种场景。
一种基于主成分分析的图像噪声水平估计方法,包括:
步骤1,将原始图像划分为若干大小相同的图像块,相邻图像块之间具有重叠部分;
步骤2,将选取的图像块组成图像数据矩阵,并求取图像数据矩阵的协方差矩阵;
步骤3,对协方差矩阵做特征分解,将所得的最小特征值作为噪声水平的估计值,并利用该噪声水平的估计值计算阈值;
步骤4,计算每个图像块对应的梯度矩阵,以及梯度矩阵对应的协方差矩阵的迹;
步骤5,选取迹小于等于阈值的图像块,重复进行步骤2~步骤4,直至连续两次计算得到的噪声水平的估计值保持不变或迭代达到预定次数后停止。
本发明从单张图像中选择低秩的没有高频成分的图像块,选择的过程利用了图像块的梯度信息与统计特征,然后利用主成分分析的方法选择图像块,对图像的噪声水平进行估计。
步骤1中,对原始图像进行划分时,首先确定一个固定尺寸的窗口,将窗口按照从左往右,从上往下的顺序依次移动固定的步长,每次移动后位于窗口内的图像即为一个图像块。图像块有重叠区域,但大部分的图像块是独立的,各图像块的大小相等,每个图像块以列向量的形式表示,所有选取的图像块的列向量构成图像数据矩阵。
初次计算时,选取所有的图像块进行步骤2的操作。
步骤5中所述的噪声水平的估计值保持不变,是指连续两次计算得到的噪声水平的估计值的差值在设定的范围之内。
图像块的梯度矩阵对应的协方差矩阵的迹表征了图像块的纹理强度,将迹小于等于阈值的图像块作为选定的图像块,重复进行步骤2~步骤4。
最后一次计算得到的噪声水平的估计值即为所需要求的噪声水平的估计值。
作为优选,所述阈值的计算公式如下:
其中,F-1(δ,α,β)为Gamma分布的累积分布函数的逆函数;α为形状参数,β为尺度参数,δ为置信度;σn为高斯噪声的标准差;N2是图像块中的像素数;Dh为水平方向的求导算子;Dv为垂直方向的求导算子;tr表示求迹操作,T表示矩阵转置操作。
作为优选,每个图像块对应的梯度矩阵的协方差矩阵如下:
其中,为图像块yi对应的梯度矩阵;Dh为水平方向的求导算子;Dv为垂直方向的求导算子;T表示矩阵转置操作。
作为优选,协方差矩阵的迹的计算公式如下:
其中,n为高斯噪声,Dh为水平方向的求导算子;Dv为垂直方向的求导算子;T表示矩阵转置操作。
与现有技术相比,本发明的有益效果为:
(1)本发明基于图像纹理来选择平滑、低秩的图像块,具体是利用图像块的梯度矩阵,利用PCA的方法来计算,能够更准确地找到图像中平滑的块。
(2)本发明提出的噪声水平估计算法是独立于场景的,算法具有准确性和鲁棒性,适用于大范围内噪声水平和各种场景。
附图说明
图1为本发明基于主成分分析的图像噪声水平估计方法的流程图。
具体实施方式
下面结合附图,对本发明做详细描述。
如图1所示,一种基于主成分分析的图像噪声水平估计方法,包括以下步骤:
(1)将输入的噪声图像分为大小为9×9的图像块,具体方式为:设定窗口的大小为9×9,然后依次由左往右,从上往下移动窗口扫描噪声图像,每次移动后位于窗口内的图像部分即为一个图像块,移动步长均为1个像素。
(2)将所有的图像块写成列向量的形式,记为yi,yi是第i个的图像块,然后将所有的图像块组合成一个大的数据矩阵,记为Y。
(3)按照如下公式,计算数据矩阵Y对应的协方差矩阵Σy:
其中:M是图像块的个数,yi为列向量。
(4)对协方差矩阵做特征分解,得到其所有的特征值,并将最小的特征值作为噪声水平的估计值,数学表达式如下:
其中,为估计的噪声水平,Σy为图像块组成的协方差矩阵,λmin表示Σy最小的特征值。
将初始噪声水平的估计值记为对应于选取所有的图像块进行计算。
(5)利用计算出的(初始即为),按照如下公式计算阈值threshold:
其中,F-1(δ,α,β)是Gamma分布的累积分布函数的逆函数,α为形状参数;β为尺度参数;δ为置信度;σn为高斯噪声的标准差;N2是图像块中的像素数;Dh和Dv,分别表示水平和垂直方向的求导算子;tr表示求迹操作。
具体实施中,N=9; Dh=[-0.500.5], T表示转置操作; δ=0.99。
(6)计算每个图像块对应的梯度矩阵
其中,Dh和Dv,分别表示水平和垂直方向的求导算子。
(7)计算每个图像块梯度矩阵对应的协方差矩阵的迹,定义为:
其中,为图像块yi对应的梯度矩阵;Dh和Dv,分别表示水平和垂直方向的求导算子。
的迹定义为:
其中,n为高斯噪声,T表示转置操作;Dh和Dv,分别表示水平和垂直方向的求导算子。
(8)上述求迹的过程对每个图像块都要算一次,将求得的每个迹的值,与threshold比较,如果小于等于threshold,则将该图像块进行标记,即记录其下标i。
(9)转到步骤(2),将步骤(8)中标记的图像块作为新的选中的图像块,并再顺序进行步骤(2)~步骤(8)。
具体实施时,当第(4)步计算出第5个噪声水平值的时候,整个过程停止,并以此结果作为最终的噪声水平的估计值。
Claims (4)
1.一种基于主成分分析的图像噪声水平估计方法,其特征在于,包括:
步骤1,将原始图像划分为若干大小相同的图像块,相邻图像块之间具有重叠部分;
步骤2,将选取的图像块组成图像数据矩阵,并求取图像数据矩阵的协方差矩阵;
步骤3,对协方差矩阵做特征分解,将所得的最小特征值作为噪声水平的估计值,并利用该噪声水平的估计值计算阈值;
步骤4,计算每个图像块对应的梯度矩阵,以及梯度矩阵对应的协方差矩阵的迹;
步骤5,选取迹小于等于阈值的图像块,重复进行步骤2~步骤4,直至连续两次计算得到的噪声水平的估计值保持不变或迭代达到预定次数后停止。
2.如权利要求1所述的基于主成分分析的图像噪声水平估计方法,其特征在于,所述阈值的计算公式如下:
其中,F-1(δ,α,β)为Gamma分布的累积分布函数的逆函数;α为形状参数,β为尺度参数,δ为置信度;σn为高斯噪声的标准差;N2是图像块中的像素数;Dh为水平方向的求导算子;Dv为垂直方向的求导算子;tr表示求迹操作,T表示矩阵转置操作。
3.如权利要求2所述的基于主成分分析的图像噪声水平估计方法,其特征在于,每个图像块对应的梯度矩阵的协方差矩阵如下:
其中,为图像块yi对应的梯度矩阵;Dh为水平方向的求导算子;Dv为垂直方向的求导算子;T表示矩阵转置操作。
4.如权利要求3所述的基于主成分分析的图像噪声水平估计方法,其特征在于,协方差矩阵的迹的计算公式如下:
其中,n为高斯噪声,Dh为水平方向的求导算子;Dv为垂直方向的求导算子;T表示矩阵转置操作。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610016782.7A CN105678774A (zh) | 2016-01-11 | 2016-01-11 | 一种基于主成分分析的图像噪声水平估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610016782.7A CN105678774A (zh) | 2016-01-11 | 2016-01-11 | 一种基于主成分分析的图像噪声水平估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105678774A true CN105678774A (zh) | 2016-06-15 |
Family
ID=56300109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610016782.7A Pending CN105678774A (zh) | 2016-01-11 | 2016-01-11 | 一种基于主成分分析的图像噪声水平估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105678774A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107295217A (zh) * | 2017-06-30 | 2017-10-24 | 中原智慧城市设计研究院有限公司 | 一种基于主成分分析的视频噪声估计方法 |
CN107451981A (zh) * | 2017-08-15 | 2017-12-08 | 四川长虹电器股份有限公司 | 基于dct和梯度协方差矩阵的图像噪声水平估计方法 |
CN107909586A (zh) * | 2017-12-11 | 2018-04-13 | 厦门美图之家科技有限公司 | 图像噪声计算方法及装置 |
CN108629264A (zh) * | 2017-03-18 | 2018-10-09 | 上海荆虹电子科技有限公司 | 用于图像处理的方法和装置 |
CN113298770A (zh) * | 2021-05-20 | 2021-08-24 | 武汉工程大学 | 图像噪声水平估计方法、装置及计算机存储介质 |
CN113658118A (zh) * | 2021-08-02 | 2021-11-16 | 维沃移动通信有限公司 | 图像噪声程度估计方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537624A (zh) * | 2015-01-05 | 2015-04-22 | 西安电子科技大学 | 基于ssim校正聚类稀疏表示的sar图像降斑方法 |
CN104637060A (zh) * | 2015-02-13 | 2015-05-20 | 武汉工程大学 | 一种基于邻域主成分分析-拉普拉斯的图像分割方法 |
-
2016
- 2016-01-11 CN CN201610016782.7A patent/CN105678774A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537624A (zh) * | 2015-01-05 | 2015-04-22 | 西安电子科技大学 | 基于ssim校正聚类稀疏表示的sar图像降斑方法 |
CN104637060A (zh) * | 2015-02-13 | 2015-05-20 | 武汉工程大学 | 一种基于邻域主成分分析-拉普拉斯的图像分割方法 |
Non-Patent Citations (2)
Title |
---|
STANISLAV PYATYKH等: "Image noise level estimation by principal component analysis", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 * |
XINHAO LIU等: "Single-Image Noise Level Estimation for Blind Denoising", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108629264A (zh) * | 2017-03-18 | 2018-10-09 | 上海荆虹电子科技有限公司 | 用于图像处理的方法和装置 |
CN108629264B (zh) * | 2017-03-18 | 2022-09-27 | 上海荆虹电子科技有限公司 | 用于图像处理的方法和装置 |
CN107295217A (zh) * | 2017-06-30 | 2017-10-24 | 中原智慧城市设计研究院有限公司 | 一种基于主成分分析的视频噪声估计方法 |
CN107295217B (zh) * | 2017-06-30 | 2020-06-12 | 中原智慧城市设计研究院有限公司 | 一种基于主成分分析的视频噪声估计方法 |
CN107451981A (zh) * | 2017-08-15 | 2017-12-08 | 四川长虹电器股份有限公司 | 基于dct和梯度协方差矩阵的图像噪声水平估计方法 |
CN107909586A (zh) * | 2017-12-11 | 2018-04-13 | 厦门美图之家科技有限公司 | 图像噪声计算方法及装置 |
CN107909586B (zh) * | 2017-12-11 | 2020-07-03 | 厦门美图之家科技有限公司 | 图像噪声计算方法及装置 |
CN113298770A (zh) * | 2021-05-20 | 2021-08-24 | 武汉工程大学 | 图像噪声水平估计方法、装置及计算机存储介质 |
CN113658118A (zh) * | 2021-08-02 | 2021-11-16 | 维沃移动通信有限公司 | 图像噪声程度估计方法、装置、电子设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105678774A (zh) | 一种基于主成分分析的图像噪声水平估计方法 | |
CN103049892B (zh) | 基于相似块矩阵秩最小化的非局部图像去噪方法 | |
CN106204447A (zh) | 基于总变差分和卷积神经网络的超分辨率重建方法 | |
CN106127688B (zh) | 一种超分辨率图像重建方法及其系统 | |
CN107358589B (zh) | 一种结合梯度直方图与低秩约束的去噪方法 | |
CN104408742B (zh) | 一种基于空间‑时间频谱联合分析的运动目标检测方法 | |
CN107705265B (zh) | 一种基于总曲率的sar图像变分去噪方法 | |
CN109118446B (zh) | 一种水下图像复原及去噪方法 | |
CN105335947A (zh) | 图像去噪方法和图像去噪装置 | |
CN107644217A (zh) | 基于卷积神经网络和相关滤波器的目标跟踪方法 | |
CN106097256B (zh) | 一种基于图像盲去模糊的视频图像模糊度检测方法 | |
CN102298774B (zh) | 基于联合相似性的非局部均值去噪方法 | |
CN104200434A (zh) | 一种基于噪声方差估计的非局部均值图像去噪方法 | |
CN107295217B (zh) | 一种基于主成分分析的视频噪声估计方法 | |
CN108932699A (zh) | 基于变换域的三维匹配调和滤波图像去噪方法 | |
CN106296649B (zh) | 一种基于水平集模型的纹理图像分割方法 | |
CN103475897A (zh) | 一种基于失真类型判断的自适应图像质量评价方法 | |
CN103150725B (zh) | 基于非局部均值的susan边缘检测方法及系统 | |
CN111768349A (zh) | 一种基于深度学习的espi图像降噪方法及系统 | |
CN109741258B (zh) | 基于重建的图像超分辨率方法 | |
CN106204502A (zh) | 基于混合阶l0正则化模糊核估计方法 | |
CN110660077A (zh) | 一种融合多个特征的多尺度目标跟踪方法 | |
CN111191680A (zh) | 一种基于非凸运动辅助的目标检测方法 | |
CN104637060A (zh) | 一种基于邻域主成分分析-拉普拉斯的图像分割方法 | |
CN103745442B (zh) | 基于非局部小波系数收缩的图像去噪方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160615 |