CN103886553A - 图像非局部均值去噪方法和系统 - Google Patents
图像非局部均值去噪方法和系统 Download PDFInfo
- Publication number
- CN103886553A CN103886553A CN201410086580.0A CN201410086580A CN103886553A CN 103886553 A CN103886553 A CN 103886553A CN 201410086580 A CN201410086580 A CN 201410086580A CN 103886553 A CN103886553 A CN 103886553A
- Authority
- CN
- China
- Prior art keywords
- image
- mrow
- search
- level
- msub
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000013598 vector Substances 0.000 claims abstract description 98
- 239000011159 matrix material Substances 0.000 claims abstract description 75
- 238000000513 principal component analysis Methods 0.000 claims abstract description 33
- 238000001914 filtration Methods 0.000 claims abstract description 15
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- NAXKFVIRJICPAO-LHNWDKRHSA-N [(1R,3S,4R,6R,7R,9S,10S,12R,13S,15S,16R,18S,19S,21S,22S,24S,25S,27S,28R,30R,31R,33S,34S,36R,37R,39R,40S,42R,44R,46S,48S,50R,52S,54S,56S)-46,48,50,52,54,56-hexakis(hydroxymethyl)-2,8,14,20,26,32,38,43,45,47,49,51,53,55-tetradecaoxa-5,11,17,23,29,35,41-heptathiapentadecacyclo[37.3.2.23,7.29,13.215,19.221,25.227,31.233,37.04,6.010,12.016,18.022,24.028,30.034,36.040,42]hexapentacontan-44-yl]methanol Chemical compound OC[C@H]1O[C@H]2O[C@H]3[C@H](CO)O[C@H](O[C@H]4[C@H](CO)O[C@H](O[C@@H]5[C@@H](CO)O[C@H](O[C@H]6[C@H](CO)O[C@H](O[C@H]7[C@H](CO)O[C@@H](O[C@H]8[C@H](CO)O[C@@H](O[C@@H]1[C@@H]1S[C@@H]21)[C@@H]1S[C@H]81)[C@H]1S[C@@H]71)[C@H]1S[C@H]61)[C@H]1S[C@@H]51)[C@H]1S[C@@H]41)[C@H]1S[C@H]31 NAXKFVIRJICPAO-LHNWDKRHSA-N 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明提供一种图像非局部均值去噪方法,包括:获取待去噪的图像,抽取预设大小的图像块,读取图像块的RGB三个颜色分量的空域数据组成矢量矩阵,并经过主成分分析和降维,生成降维后的残差矢量矩阵;对图像块,根据预设的S级搜索区域,在搜索区域内根据预设的每级搜索区域的搜索步长,以图像块的中心像素点为中心从上到下,从左到右移动预设大小的搜索窗口,得到搜索块;其中,各级搜索区域的搜索步长依次递增;计算搜索块对应的残差矢量矩阵与像素点对应的矢量的欧几里得距离,再根据欧几里得距离计算图像块的权重值,最后根据权重值得到图像块中心像素点的滤波值,完成图像的去噪处理。本发明还提供对应的系统,时间复杂度较低,去噪速度快。
Description
技术领域
本发明涉及图像处理技术领域,特别是涉及一种图像非局部均值去噪方法,以及一种图像非局部均值去噪系统。
背景技术
数字图像与视频在成像过程中,往往引入多种噪声,包括大气噪声,摄像器材的传感器噪声,尤其是光照度较低的成像图像,其引入的噪声往往比较严重。图像去噪技术包括空域去噪与频域去噪技术。频域去噪由于需要经过空间域到频率域的转换,往往复杂度较高,而且传统的维纳滤波去噪效果并不好。传统的空域滤波去噪方法,如均值滤波,中值滤波,统计排序滤波器等滤波去噪后往往导致图像模糊,去噪效果较差。目前效果较好的空域滤波方法是非局部均值去噪方法(NLM,Non-Local Means)。这种方法计算当前图像块与该块周围的图像块的相似性,并由相似性计算权重值,令相似块的中心点的加权平均值作为当前图像块的中心像素的滤波值。这种方法的算法复杂度简单,且能在保留图像细节的同时较好地去除图像噪声,但其时间复杂度较高,复杂度为O(M×N×(2r+1)2×L2),其中M与N分别为图像的长度与宽度,r为图像当前块的搜索半径,L为图像块的长宽大小(宽度与长度相同)。随着搜索半径r的增大,图像去噪效果更好,但时间复杂度迅速增高。在实际应用中,往往需要快速的图像/视频去噪技术,实现实时处理,较高的时间复杂度限制了NLM算法的应用。
当图像噪声较严重时,图像细节被噪声掩盖,会对传统的NLM算法中计算图像块的相似性造成影响,无法获得准确的相似性块。同时,研究表明在大多数情况下图像块的相似性会随着搜索半径的增大而下降,较远距离的搜索块的权重较低,而传统的NLM算法的搜索步长时间复杂度较高。
发明内容
基于此,本发明提供一种图像非局部均值去噪方法和系统,其时间复杂度较低,图像去噪速度快。
一种图像非局部均值去噪方法,包括如下步骤:
获取待去噪的图像,以所述图像中每个像素点为中心抽取预设大小的图像块,读取每个图像块的RGB三个颜色分量的空域数据组成矢量矩阵,并经过主成分分析和降维,生成降维后的残差矢量矩阵;
对所述图像块,根据预设的S级搜索区域,在各级搜索区域内根据预设的每级搜索区域的搜索步长,以所述图像块的中心像素点为中心从上到下,从左到右移动预设大小的搜索窗口,得到搜索块;其中,S为大于或等于3的整数;第1至第S-1级的搜索区域都为所述图像中以In i+x,j+y为中心,长宽为(2rs+1)的方形区域,且不包括低于当前级的搜索区域,In为所述待去噪的图像,rs为第S级搜索区域的搜索半径,第S级搜索区域为所述图像中第1~S-1级搜索区域以外的区域;x和y分别为垂直和水平方向相对于像素点(i,j)的位移,各级搜索区域的搜索步长依次递增;
计算所述搜索块对应的所述残差矢量矩阵与所述像素点对应的矢量的欧几里得距离,再根据所述欧几里得距离计算所述图像块的权重值,最后根据所述权重值得到所述图像块中心像素点的滤波值,完成所述图像的去噪处理。
一种图像非局部均值去噪系统,包括:
降维模块,用于获取待去噪的图像,以所述图像中每个像素点为中心抽取预设大小的图像块,读取每个图像块的RGB三个颜色分量的空域数据组成矢量矩阵,并经过主成分分析和降维,生成降维后的残差矢量矩阵;
搜索模块,用于对所述图像块,根据预设的S级搜索区域,在各级搜索区域内根据预设的每级搜索区域的搜索步长,以所述图像块的中心像素点为中心从上到下,从左到右移动预设大小的搜索窗口,得到搜索块;其中,S为大于或等于3的整数;第1至第S-1级的搜索区域都为所述图像中以In i+x,j+y为中心,长宽为(2rs+1)的方形区域,且不包括低于当前级的搜索区域,In为所述待去噪的图像,rs为第S级搜索区域的搜索半径,第S级搜索区域为所述图像中第1~S-1级搜索区域以外的区域;x和y分别为垂直和水平方向相对于像素点(i,j)的位移,各级搜索区域的搜索步长依次递增;
去噪模块,用于计算所述搜索块对应的所述残差矢量矩阵与所述像素点对应的矢量的欧几里得距离,再根据所述欧几里得距离计算所述图像块的权重值,最后根据所述权重值得到所述图像块中心像素点的滤波值,完成所述图像的去噪处理。
上述图像非局部均值去噪方法和系统,在噪声图像中抽取图像块,读取空域数据构成矢量矩阵,再经过主成分分析和降维,生成降维后的残差矢量矩阵,采用降维的方式,可在图像噪声较严重时,获得更准确的图像相似性块;接着在设置搜索区域时采用变步长的搜索方式对图像进行去噪,可显著降低图像去噪的时间复杂度,提高图像的去噪速度。
附图说明
图1为本发明图像非局部均值去噪方法在一实施例中的流程示意图。
图2为本发明图像非局部均值去噪方法在一实施例中搜索区域的示意图。
图3为本发明图像非局部均值去噪方法在一实施例中搜索区域的另一示意图。
图4为本发明图像非局部均值去噪系统在一实施例中的结构示意图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细说明,但本发明的实施方式不限于此。
如图1所示,是本发明一种图像非局部均值去噪方法的流程示意图,包括如下步骤:
S11、获取待去噪的图像,以所述图像中每个像素点为中心抽取预设大小的图像块,读取每个图像块的RGB三个颜色分量的空域数据组成矢量矩阵,并经过主成分分析和降维,生成降维后的残差矢量矩阵;
本步骤对图像进行PCA(主成分分析,Principal Component Analysis)预处理,先依次读取图像中的像素点,以每个像素点为中心抽取图像块,通过每个图像块包含的RGB三个颜色分量的空域数据,每个图像块用矢量表示,所有矢量组成矢量矩阵H;
具体地,以所述图像中每个像素点为中心抽取预设大小的图像块的步骤可为:
对所述图像进行边界对称拓展,得到拓展后的图像;其中,所述图像的大小为M×N,M为图像长度,N为图像宽度,拓展后的图像大小为(M+l)×(N+l);
所述经过主成分分析和降维,生成降维后的残差矢量矩阵的步骤为:
从所述矢量矩阵中选取R个矢量构成一随机矢量矩阵进行主成分分析,获得PCA投影矩阵;
将矢量矩阵中的每个矢量减去所述选取的R个矢量的均值,得到残差矢量矩阵;
将PCA投影矩阵与残差矢量矩阵相乘,得到所述降维后的残差矢量矩阵;
本实施例中进行主成分分析和降维,可在图像噪声较严重时,获得更准确的图像相似性块,获得更好的图像去噪的效果,提高方法的鲁棒性,同时由于维度的降低,可提高图像的去噪速度。
S12、对所述图像块,根据预设的S级搜索区域,在各级搜索区域内根据预设的每级搜索区域的搜索步长,以所述图像块的中心像素点为中心从上到下,从左到右移动搜索窗口,得到搜索块;其中,S为大于或等于3的整数;第1至第S-1级的搜索区域都为所述图像中以In i+x,j+y为中心,长宽为(2rs+1)的方形区域,且不包括低于当前级的搜索区域(也即是,当前一级的搜索区域,其是以In i+x,j+y为中心,长宽为(2rs+1)构成一方形区域,在该方形区域中去除低于当前级的其他级搜索区域后,剩余的区域为当前级的搜索区域),In为所述待去噪的图像,rs为第S级搜索区域的搜索半径,第S级搜索区域为所述图像中第1~S-1级搜索区域以外的区域;x和y分别为垂直和水平方向相对于像素点(i,j)的位移,各级搜索区域的搜索步长依次递增;
本步骤对图像进行变步长搜索,对于每个图像块,设置S级搜索区域,S的个数可为大于或等于3的整数,对于第1~S-1级的搜索区域,都是以像素点(i,j)为中心,长宽为(2rs+1)的方形区域(且不包括低于当前级的搜索区域),最后的S级搜索区域则为图像中第1~S-1级搜索区域以外的剩余区域;rs为第s级搜索区域的搜索半径;其中,如果该方形区域的某一侧边界越过图像In本身的边界,则以图像的边界为该方形区域的该侧边界;
在每级搜索区域内,根据该搜索区域预设的搜索步长,从上到下,从左到右移动搜索窗口,获得搜索块Pi+x,j+y,x和y分别为垂直和水平方向相对于(i,j)的位移,其中,各级搜索区域的搜索步长依次递增,能有效提高后续图像去噪的速度。
S13、计算所述搜索块对应的所述残差矢量矩阵与所述像素点对应的矢量的欧几里得距离,再根据所述欧几里得距离计算所述图像块的权重值,最后根据所述权重值得到所述图像块中心像素点的滤波值,完成所述图像的去噪处理;
具体地,可根据下式计算所述欧几里得距离:
其中,di,j,k为所述欧几里得距离,k表示矢量vi,j,k的第k维度,q为所述残差矢量矩阵维度。
根据下式计算所述图像块的权重值:
其中,wi+x,j+y为所述图像块的权重值,σ为噪声标准差。
根据下式计算所述像素点的滤波值:
其中,z表示图像块的RGB色度空间,Ω表示所述搜索区域。
接下来再通过一具体实施例详细阐述本发明。
以M×N大小的噪声图像In为例子进行说明,M为图像长度,N为图像宽度。本实施例包括PCA预处理和变步长NLM(非局部均值,Non-local Mean)滤波去噪两个步骤。
(1)PCA预处理
①对M×N大小的噪声图像In进行边界对称拓展,得到(M+3)×(N+3)大小的图像I'n,然后在图像I'n中从上到下从左到右,以像素(i,j)为中心,抽取7×7大小的图像块,其中3≤i≤M-4,3≤j≤N-4。每个块包含RGB三个颜色分量的空域数据,每个块用矢量表示,所有矢量组成矢量矩阵H。
②从矢量矩阵H中选择R个矢量并组成新的矢量矩阵,称为随机矢量矩阵,R的计算方式可为R=min(M×N,5000)。对这个随机矢量矩阵进行PCA(主成分分析)分析,获得PCA投影矩阵,其维度为q(一般可为30),以及这R个矢量的均值。
(2)变步长NLM去噪
①对于噪声图像In中以(i,j)为中心的7×7大小的图像块Pi,j,其中0≤i≤M-1,0≤j≤N-1;如图2中矩形虚线内的7×7大小的块,设置3级搜索区域,第s(s=1,2)级的搜索区域为以In i+x,j+y为中心,长宽为(2rs+1)的方形区域,rs为第s级搜索区域的搜索半径,如果该方形区域的某一侧边界越过图像In本身的边界,则以图像的边界为该方形区域的该侧边界。第3级搜索区域为第1级与第2级搜索区域外的区域。图2中r1=5,r2=11,搜索窗口大小为7×7。第s级区域内的搜索窗口的搜索步长为s,即第1,2,3级区域的搜索窗口的搜索步长依次递增,分别为1,2,3。
例如图3中,灰色7×7方形块Pi-3,j-5位于第1级搜索区域,带斜线的7×7方形块Pi-11,j-11,位于第2级搜索区域。在图3中,在第1级搜索区域内,矩形搜索窗口从Pi-3,j-5块移动到Pi-2,j-5块,搜索步长(移动步长)为1,而在第2级搜索内,矩形窗口从Pi-11,j-11块移动到Pi-9,j-11,搜索步长(移动步长)为2。
②在搜索区域内根据该搜索区域的搜索步长,从上到下,从左到右移动搜索窗口,获得搜索块Pi+x,j+y,x和y分别为垂直和水平方向相对于(i,j)的位移,计算该搜索块所对应的残差矢量vi+x,j+y与vi,j的欧几里得距离:
其中k表示矢量vi,j,k的第k维度。
③根据欧几里得距离di,j,k计算图像块Pi+x,j+y的权重值:
其中σ为噪声标准差。
④根据以下公式计算像素Ii,j的滤波值:
上式中z表示图像的RGB色度空间,Ω表示搜索区域;变步长NLM滤波去噪结束。
对整幅图像的每个像素采用以上①-④步骤进行变步长NLM去噪,即完成整幅图像的去噪处理。
本实施例采用PCA分析与降维的方式,可在图像噪声较严重时,获得更准确的图像相似性块,获得更好的图像去噪的效果,提高算法的鲁棒性。同时由于维度的降低,可提高图像的去噪速度。
传统的非局部均值去噪方法采用固定步长(步长为1)的搜索方式,在搜索区域较大时,搜索次数较多,计算量较大,而本实施例采用变步长的搜索方式,在搜索半径较大时,搜索步长增大,从而减少了搜索的次数,可降低图像的时间复杂度,提高图像的去噪速度。如对于7×7的搜索窗口,PCA降维后维度q=30。当第一级搜索区域的半径为r1=5,第2级搜索区域的半径为r2=11,对于M×N大小的图像,上述实施例与传统NLM方法的时间复杂度比较如下表1所示。随着搜索半径的增大,本实施例的时间复杂度较传统方法下降的比例更高,从而提高图像去噪的速度。
表1本实施例与传统方法的时间复杂度比较
搜索半径r | 传统NLM方法 | 本发明方法 | 降低百分比 |
3 | O(2401×M×N) | O(1470×M×N) | 38.8% |
5 | O(1470×M×N) | O(3630×M×N) | 38.8% |
7 | O(11025×M×N) | O(4470×M×N) | 59.5% |
9 | O(17689×M×N) | O(5550×M×N) | 68.6% |
11 | O(25921×M×N) | O(6870×M×N) | 73.5% |
14 | O(41209×M×N) | O(7950×M×N) | 80.7% |
本发明还提供一种图像非局部均值去噪系统,如图4所示,包括:
降维模块41,用于获取待去噪的图像,以所述图像中每个像素点为中心抽取预设大小的图像块,读取每个图像块的RGB三个颜色分量的空域数据组成矢量矩阵,并经过主成分分析和降维,生成降维后的残差矢量矩阵;
本模块对图像进行PCA(主成分分析,Principal Component Analysis)预处理,先依次读取图像中的像素点,以每个像素点为中心抽取图像块,通过每个图像块包含的RGB三个颜色分量的空域数据,每个图像块用矢量表示,所有矢量组成矢量矩阵H;
具体地,所述降维模块还可用于:
对所述图像进行边界对称拓展,得到拓展后的图像;其中,所述图像的大小为M×N,M为图像长度,N为图像宽度,拓展后的图像大小为(M+l)×(N+l);
所述降维模块还可用于:从所述矢量矩阵中选取R个矢量构成一随机矢量矩阵进行主成分分析,获得PCA投影矩阵;将矢量矩阵中的每个矢量减去所述选取的R个矢量的均值,得到残差矢量矩阵;将PCA投影矩阵与残差矢量矩阵相乘,得到所述降维后的残差矢量矩阵。
本实施例中进行主成分分析和降维,可在图像噪声较严重时,获得更准确的图像相似性块,获得更好的图像去噪的效果,提高方法的鲁棒性,同时由于维度的降低,可提高图像的去噪速度。
搜索模块42,用于对所述图像块,根据预设的S级搜索区域,在各级搜索区域内根据预设的每级搜索区域的搜索步长,以所述图像块的中心像素点为中心从上到下,从左到右移动预设大小的搜索窗口,得到搜索块;其中,S为大于或等于3的整数;第1至第S-1级的搜索区域都为所述图像中以In i+x,j+y为中心,长宽为(2rs+1)的方形区域(且不包括低于当前级的搜索区域),In为所述待去噪的图像,rs为第S级搜索区域的搜索半径,第S级搜索区域为所述图像中第1~S-1级搜索区域以外的区域;x和y分别为垂直和水平方向相对于像素点(i,j)的位移,各级搜索区域的搜索步长依次递增;
本模块对图像进行变步长搜索,对于每个图像块,设置S级搜索区域,S的个数可为大于或等于3的整数,对于第1~S-1级的搜索区域,都是以像素点(i,j)为中心,长宽为(2rs+1)的方形区域(且不包括低于当前级的搜索区域),最后的S级搜索区域则为图像中第1~S-1级搜索区域以外的剩余区域;rs为第s级搜索区域的搜索半径;其中,如果该方形区域的某一侧边界越过图像In本身的边界,则以图像的边界为该方形区域的该侧边界;
在每级搜索区域内,根据该搜索区域预设的搜索步长,从上到下,从左到右移动搜索窗口,获得搜索块Pi+x,j+y,x和y分别为垂直和水平方向相对于(i,j)的位移,其中,各级搜索区域的搜索步长依次递增,能有效提高后续图像去噪的速度。
去噪模块43,用于计算所述搜索块对应的所述残差矢量矩阵与所述像素点对应的矢量的欧几里得距离,再根据所述欧几里得距离计算所述图像块的权重值,最后根据所述权重值得到所述图像块中心像素点的滤波值,完成所述图像的去噪处理。
具体地,可根据下式计算所述欧几里得距离:
其中,di,j,k为所述欧几里得距离,k表示矢量vi,j,k的第k维度,q为所述残差矢量矩阵维度。
根据下式计算所述图像块的权重值:
其中,wi+x,j+y为所述图像块的权重值,σ为噪声标准差。
根据下式计算所述像素点的滤波值:
其中,z表示图像块的RGB色度空间,Ω表示所述搜索区域。
本发明图像非局部均值去噪方法和系统,在噪声图像中抽取图像块,读取空域数据构成矢量矩阵,再经过主成分分析和降维,生成降维后的残差矢量矩阵,采用降维的方式,可在图像噪声较严重时,获得更准确的图像相似性块;接着在设置搜索区域时采用变步长的搜索方式对图像进行去噪,可显著降低图像的时间复杂度,提高图像的去噪速度。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种图像非局部均值去噪方法,其特征在于,包括如下步骤:
获取待去噪的图像,以所述图像中每个像素点为中心抽取预设大小的图像块,读取每个图像块的RGB三个颜色分量的空域数据组成矢量矩阵,并经过主成分分析和降维,生成降维后的残差矢量矩阵;
对所述图像块,根据预设的S级搜索区域,在各级搜索区域内根据预设的每级搜索区域的搜索步长,以所述图像块的中心像素点为中心从上到下,从左到右移动预设大小的搜索窗口,得到搜索块;其中,S为大于或等于3的整数;第1至第S-1级的搜索区域都为所述图像中以Ini+x,j+y为中心,长宽为2rs+1的方形区域,且不包括低于当前级的搜索区域,In为所述待去噪的图像,rs为第S级搜索区域的搜索半径,第S级搜索区域为所述图像中第1~S-1级搜索区域以外的区域;x和y分别为垂直和水平方向相对于像素点(i,j)的位移,各级搜索区域的搜索步长依次递增;
计算所述搜索块对应的所述残差矢量矩阵与所述像素点对应的矢量的欧几里得距离,再根据所述欧几里得距离计算所述图像块的权重值,最后根据所述权重值得到所述图像块中心像素点的滤波值,完成所述图像的去噪处理。
2.根据权利要求1所述的图像非局部均值去噪方法,其特征在于,以所述图像中每个像素点为中心抽取预设大小的图像块的步骤为:
对所述图像进行边界对称拓展,得到拓展后的图像;其中,所述图像的大小为M×N,M为图像长度,N为图像宽度,拓展后的图像大小为(M+l)×(N+l);
在所述拓展后的图像,从上到下从左到右,依次以每个像素点(i,j)为中心,抽取t×t大小的图像块,其中l≤i≤M-l-1,l≤j≤N-l-1,t为预设的图像块长度,并且t为奇数,l为预设的拓展长度,
3.根据权利要求1所述的图像非局部均值去噪方法,其特征在于,所述经过主成分分析和降维,生成降维后的残差矢量矩阵的步骤为:
从所述矢量矩阵中选取R个矢量构成一随机矢量矩阵进行主成分分析,获得PCA投影矩阵;
将矢量矩阵中的每个矢量减去所述选取的R个矢量的均值,得到残差矢量矩阵;
将PCA投影矩阵与残差矢量矩阵相乘,得到所述降维后的残差矢量矩阵。
4.根据权利要求1所述的图像非局部均值去噪方法,其特征在于,根据下式计算所述欧几里得距离:
其中,di,j,k为所述欧几里得距离,k表示矢量vi,j,k的第k维度,q为所述残差矢量矩阵维度。
5.根据权利要求4所述的图像非局部均值去噪方法,其特征在于,根据下式计算所述图像块的权重值:
其中,wi+x,j+y为所述图像块的权重值,σ为噪声标准差。
6.根据权利要求4所述的图像非局部均值去噪方法,其特征在于,根据下式计算所述像素点的滤波值:
其中,z表示图像块的RGB色度空间,Ω表示所述搜索区域。
7.根据权利要求1所述的图像非局部均值去噪方法,其特征在于,如果该方形区域的某一侧边界越过所述图像本身的边界,则以所述图像的边界为所述方形区域的侧边界。
8.一种图像非局部均值去噪系统,其特征在于,包括:
降维模块,用于获取待去噪的图像,以所述图像中每个像素点为中心抽取预设大小的图像块,读取每个图像块的RGB三个颜色分量的空域数据组成矢量矩阵,并经过主成分分析和降维,生成降维后的残差矢量矩阵;
搜索模块,用于对所述图像块,根据预设的S级搜索区域,在各级搜索区域内根据预设的每级搜索区域的搜索步长,以所述图像块的中心像素点为中心从上到下,从左到右移动预设大小的搜索窗口,得到搜索块;其中,S为大于或等于3的整数;第1至第S-1级的搜索区域都为所述图像中以In i+x,j+y为中心,长宽为(2rs+1)的方形区域,且不包括低于当前级的搜索区域,In为所述待去噪的图像,rs为第S级搜索区域的搜索半径,第S级搜索区域为所述图像中第1~S-1级搜索区域以外的区域;x和y分别为垂直和水平方向相对于像素点(i,j)的位移,各级搜索区域的搜索步长依次递增;
去噪模块,用于计算所述搜索块对应的所述残差矢量矩阵与所述像素点对应的矢量的欧几里得距离,再根据所述欧几里得距离计算所述图像块的权重值,最后根据所述权重值得到所述图像块中心像素点的滤波值,完成所述图像的去噪处理。
10.根据权利要求8所述的图像非局部均值去噪系统,其特征在于,所述降维模块还用于:从所述矢量矩阵中选取R个矢量构成一随机矢量矩阵进行主成分分析,获得PCA投影矩阵;将矢量矩阵中的每个矢量减去所述选取的R个矢量的均值,得到残差矢量矩阵;将PCA投影矩阵与残差矢量矩阵相乘,得到所述降维后的残差矢量矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410086580.0A CN103886553B (zh) | 2014-03-10 | 2014-03-10 | 图像非局部均值去噪方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410086580.0A CN103886553B (zh) | 2014-03-10 | 2014-03-10 | 图像非局部均值去噪方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103886553A true CN103886553A (zh) | 2014-06-25 |
CN103886553B CN103886553B (zh) | 2017-02-01 |
Family
ID=50955428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410086580.0A Expired - Fee Related CN103886553B (zh) | 2014-03-10 | 2014-03-10 | 图像非局部均值去噪方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103886553B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104700375A (zh) * | 2015-03-27 | 2015-06-10 | 麦克奥迪(厦门)医疗诊断系统有限公司 | 一种基于主成分分析的病理图像视觉效果改善方法 |
CN106023093A (zh) * | 2016-05-09 | 2016-10-12 | 吉林大学 | 一种基于改进图像块匹配的非局部均值图像去噪方法 |
CN106204462A (zh) * | 2015-05-04 | 2016-12-07 | 南京邮电大学 | 基于图像多特征融合的非局部均值去噪方法 |
CN107786780A (zh) * | 2017-11-03 | 2018-03-09 | 深圳Tcl新技术有限公司 | 视频图像降噪方法、装置及计算机可读存储介质 |
CN108154481A (zh) * | 2016-12-06 | 2018-06-12 | 展讯通信(上海)有限公司 | 图像去噪方法、装置及电子设备 |
CN110097511A (zh) * | 2019-04-11 | 2019-08-06 | 成都微光集电科技有限公司 | 一种图像降噪方法 |
CN111247554A (zh) * | 2017-11-06 | 2020-06-05 | Eizo株式会社 | 图像处理装置、图像处理方法以及图像处理程序 |
CN111260580A (zh) * | 2020-01-17 | 2020-06-09 | 珠海全志科技股份有限公司 | 一种基于图像金字塔的图像去噪方法、计算机装置及计算机可读存储介质 |
CN111724325A (zh) * | 2020-06-24 | 2020-09-29 | 湖南国科微电子股份有限公司 | 一种三边滤波图像处理方法及装置 |
CN113034387A (zh) * | 2021-03-05 | 2021-06-25 | 成都国科微电子有限公司 | 一种图像去噪方法、装置、设备及介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100215282A1 (en) * | 2009-02-23 | 2010-08-26 | Van Beek Petrus J L | Methods and Systems for Imaging Processing |
CN102298773A (zh) * | 2011-09-19 | 2011-12-28 | 西安电子科技大学 | 一种形状自适应的非局部均值去噪方法 |
CN102436646A (zh) * | 2011-11-07 | 2012-05-02 | 天津大学 | 基于压缩感知的ccd噪声估计方法 |
US20120121162A1 (en) * | 2010-11-11 | 2012-05-17 | Samsung Electronics Co., Ltd. | Filtering apparatus and method for high precision restoration of depth image |
CN102609914A (zh) * | 2012-01-17 | 2012-07-25 | 天津大学 | 图像传感器的信号相关噪声估计方法 |
CN102930508A (zh) * | 2012-08-30 | 2013-02-13 | 西安电子科技大学 | 基于图像残余信号的非局部均值图像去噪方法 |
CN103020918A (zh) * | 2013-01-09 | 2013-04-03 | 西安电子科技大学 | 基于形状自适应邻域均值的非局部均值去噪方法 |
CN103093434A (zh) * | 2013-01-27 | 2013-05-08 | 西安电子科技大学 | 基于奇异值分解的非局部维纳滤波图像去噪方法 |
-
2014
- 2014-03-10 CN CN201410086580.0A patent/CN103886553B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100215282A1 (en) * | 2009-02-23 | 2010-08-26 | Van Beek Petrus J L | Methods and Systems for Imaging Processing |
US20120121162A1 (en) * | 2010-11-11 | 2012-05-17 | Samsung Electronics Co., Ltd. | Filtering apparatus and method for high precision restoration of depth image |
CN102298773A (zh) * | 2011-09-19 | 2011-12-28 | 西安电子科技大学 | 一种形状自适应的非局部均值去噪方法 |
CN102436646A (zh) * | 2011-11-07 | 2012-05-02 | 天津大学 | 基于压缩感知的ccd噪声估计方法 |
CN102609914A (zh) * | 2012-01-17 | 2012-07-25 | 天津大学 | 图像传感器的信号相关噪声估计方法 |
CN102930508A (zh) * | 2012-08-30 | 2013-02-13 | 西安电子科技大学 | 基于图像残余信号的非局部均值图像去噪方法 |
CN103020918A (zh) * | 2013-01-09 | 2013-04-03 | 西安电子科技大学 | 基于形状自适应邻域均值的非局部均值去噪方法 |
CN103093434A (zh) * | 2013-01-27 | 2013-05-08 | 西安电子科技大学 | 基于奇异值分解的非局部维纳滤波图像去噪方法 |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104700375A (zh) * | 2015-03-27 | 2015-06-10 | 麦克奥迪(厦门)医疗诊断系统有限公司 | 一种基于主成分分析的病理图像视觉效果改善方法 |
CN106204462A (zh) * | 2015-05-04 | 2016-12-07 | 南京邮电大学 | 基于图像多特征融合的非局部均值去噪方法 |
CN106023093A (zh) * | 2016-05-09 | 2016-10-12 | 吉林大学 | 一种基于改进图像块匹配的非局部均值图像去噪方法 |
CN108154481B (zh) * | 2016-12-06 | 2021-09-21 | 展讯通信(上海)有限公司 | 图像去噪方法、装置及电子设备 |
CN108154481A (zh) * | 2016-12-06 | 2018-06-12 | 展讯通信(上海)有限公司 | 图像去噪方法、装置及电子设备 |
CN107786780B (zh) * | 2017-11-03 | 2021-04-16 | 深圳Tcl新技术有限公司 | 视频图像降噪方法、装置及计算机可读存储介质 |
CN107786780A (zh) * | 2017-11-03 | 2018-03-09 | 深圳Tcl新技术有限公司 | 视频图像降噪方法、装置及计算机可读存储介质 |
CN111247554A (zh) * | 2017-11-06 | 2020-06-05 | Eizo株式会社 | 图像处理装置、图像处理方法以及图像处理程序 |
CN110097511A (zh) * | 2019-04-11 | 2019-08-06 | 成都微光集电科技有限公司 | 一种图像降噪方法 |
CN110097511B (zh) * | 2019-04-11 | 2021-06-15 | 成都微光集电科技有限公司 | 一种图像降噪方法 |
CN111260580A (zh) * | 2020-01-17 | 2020-06-09 | 珠海全志科技股份有限公司 | 一种基于图像金字塔的图像去噪方法、计算机装置及计算机可读存储介质 |
CN111724325A (zh) * | 2020-06-24 | 2020-09-29 | 湖南国科微电子股份有限公司 | 一种三边滤波图像处理方法及装置 |
CN111724325B (zh) * | 2020-06-24 | 2023-10-31 | 湖南国科微电子股份有限公司 | 一种三边滤波图像处理方法及装置 |
CN113034387A (zh) * | 2021-03-05 | 2021-06-25 | 成都国科微电子有限公司 | 一种图像去噪方法、装置、设备及介质 |
CN113034387B (zh) * | 2021-03-05 | 2023-07-14 | 成都国科微电子有限公司 | 一种图像去噪方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103886553B (zh) | 2017-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103886553A (zh) | 图像非局部均值去噪方法和系统 | |
CN107945125B (zh) | 一种融合频谱估计法和卷积神经网络的模糊图像处理方法 | |
US9471964B2 (en) | Non-local mean-based video denoising method and apparatus | |
CN105096280B (zh) | 处理图像噪声的方法及装置 | |
CN103761710B (zh) | 基于边缘自适应的高效图像盲去模糊方法 | |
CN103533214B (zh) | 一种基于卡尔曼滤波和双边滤波的视频实时去噪方法 | |
CN108564597B (zh) | 一种融合高斯混合模型和h-s光流法的视频前景目标提取方法 | |
CN106204447A (zh) | 基于总变差分和卷积神经网络的超分辨率重建方法 | |
CN107730536B (zh) | 一种基于深度特征的高速相关滤波物体追踪方法 | |
CN107749987B (zh) | 一种基于块运动估计的数字视频稳像方法 | |
CN109215053B (zh) | 一种无人机航拍视频中含暂停状态的运动车辆检测方法 | |
CN104796582B (zh) | 基于随机喷射retinex的视频图像去噪与增强方法及装置 | |
CN103841298B (zh) | 一种基于颜色恒量和几何不变特征的视频稳像方法 | |
CN102572223B (zh) | 一种用于视频降噪的相似块搜索方法 | |
CN106097256B (zh) | 一种基于图像盲去模糊的视频图像模糊度检测方法 | |
CN103268482B (zh) | 一种低复杂度的手势提取和手势深度获取方法 | |
Kryjak et al. | Real-time implementation of the ViBe foreground object segmentation algorithm | |
CN102222327A (zh) | 基于Treelet变换和最小均方误差估计的图像去噪方法 | |
CN106254720A (zh) | 一种基于联合正则化的视频超分辨率重建方法 | |
CN109741258B (zh) | 基于重建的图像超分辨率方法 | |
Alam et al. | Space-variant blur kernel estimation and image deblurring through kernel clustering | |
CN103745442B (zh) | 基于非局部小波系数收缩的图像去噪方法 | |
CN111260590B (zh) | 图像降噪方法及相关产品 | |
CN112750089A (zh) | 基于局部块最大和最小像素先验的光学遥感影像去雾方法 | |
CN105069764B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: Kezhu road high tech Industrial Development Zone, Guangzhou city of Guangdong Province, No. 233 510670 Patentee after: VTRON GROUP Co.,Ltd. Address before: 510670 Guangdong city of Guangzhou province Kezhu Guangzhou high tech Industrial Development Zone, Road No. 233 Patentee before: VTRON TECHNOLOGIES Ltd. |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170201 |