一种权值窗口自适应非局部均值图像去噪方法技术

本发明专利技术公开了一种权值窗口自适应非局部均值图像去噪方法，能够根据图像局部结构特性自适应控制权值窗口的大小，在保护边缘结构的同时抑制噪声，从而使得图像质量显著提高，它包括如下步骤：第一步，初始化读入一帧噪声图像；第二步，构建结构张量矩阵；第三步，根据所构建的结构张量矩阵，构建边缘结构指示符，定位像素点所在区域特征；第四步，利用边缘结构指示符对图像进行区域分类；第五步，根据每个像素点所属区域的类别，确定每个像素点的邻域大小；第六步，根据所确定的邻域大小，构建邻域间的相似性度量函数；第七步，筛选相似性最高的S个点；第八步，构建去噪模型，得到去噪图像。

【技术实现步骤摘要】

本 专利技术属于图像处理领域，涉及。
技术介绍
图像去噪技术属于图像处理和底层视觉中的关键性问题，是后续模式识别和高层理解的基础，有着广泛的应用需求，该技术可以应用到交通监控、军事、医学等多个领域。因此，通过图像去噪技术提高图像质量，具有重要的理论意义和现实价值，收到国内外学术界和商业界的极大关注。非局部均值方法是当前新兴而又非常有效的一种图像去噪方法，它利用大多数图像内存在很多相似性图像块的特点，在图像中搜索灰度相似的块进行匹配来估计噪声点的值。非局部均值方法权值窗口大小的选择至关重要，选择较大的权值窗口容易可以很好地去除噪声，但同时会导致边缘结构的模糊化，反之，选择较小的权值窗口，可以很好地保护边缘结构，但噪声不能很好地得到抑制。由于图像边缘结构是图像中最重要的视觉特征，保护图像边沿结构尤为重要。为了达到这比目标，非局部均值去噪方法中的权值窗口应该具有结构自适应性，即当待处理的像素点为边缘结构上的点时，采用较小的权值窗口函数，以更好地保护边缘结构；当待处理的像素点属于平坦区域或噪声点是，采用较大的权值窗口函数，以抑制噪声。为此，本专利技术提出了一种基于图像区域自适应的非局部均值去噪方法，该方法能够根据图像的局部结构特性自适应选取权值窗口的大小，从而可以更好地保护边缘结构和抑制噪声。
技术实现思路
技术问题:本专利技术提供一种可以更好地保护图像边缘结构和抑制噪声的权值窗口自适应非局部均值图像去噪方法。技术方案:本专利技术的权值窗口自适应非局部均值图像去噪方法，包括如下步骤:步骤1:读入一帧大小为M1XM2X 3的噪声彩色图像uQ，其中M1和M2为正整数，分别表示图像矩阵的行数和列数，然后将输入的噪声彩色图像从RGB彩色空间转换到YCbCr彩色空间，转换后的噪声图像记为f，大小为M1XM2X 3 ；步骤2:构建噪声图像f所有像素点(i，j)的结构张量矩阵T。(i, j)本文档来自技高网...

【技术保护点】
一种权值窗口自适应非局部均值图像去噪方法，其特征在于，该方法包括以下步骤：步骤1：读入一帧大小为M1×M2×3的噪声彩色图像u0，其中M1和M2为正整数，分别表示图像矩阵的行数和列数，然后将输入的噪声彩色图像从RGB彩色空间转换到YCbCr彩色空间，转换后的噪声图像记为f，大小为M1×M2×3；步骤2：构建噪声图像f所有像素点(i,j)的结构张量矩阵Tσ(i,j)Tσ(i,j)=Gσ(i,j)*(fx(i,j))2Gσ(i,j)*fx(i,j)fy(i,j)Gσ(i,j)*fx(i,j)fy(i,j)Gσ(i,j)*(fy(i,j))2,]]>其中，fx(i,j)表示噪声图像f在点(i,j)处沿坐标轴x方向的导数，fy(i,j)表示噪声图像f在点(i,j)处沿坐标轴y方向的导数，Gσ(i,j)表示标准差为σ的高斯核函数，其表达式为：G(i,j)=12πσ2exp[-(i2+j2)/2σ2]]]>其中，exp表示公知的以e为底的指数函数，标准差σ取值范围为0.5；步骤3：构建边缘结构指示符，具体方法为：利用雅克比方法求得结构张量矩阵Tσ(i,j)在噪声图像f中所有像素点(i,j)处的主特征值λ1(i,j)和次特征值λ2(i,j)分别为：λ1(i,j)=12{Gσ(i,j)*(fx(i,j))2+Gσ(i,j)*(fy(i,j))2+(Gσ(i,j)*(fx(i,j))2-Gσ(i,j)*(fy(i,j))2)2+4(Gσ(i,j)*fx(i,j)fy(i,j))2},]]>λ2(i,j)=12{Gσ(i,j)*(fx(i,j))2+Gσ(i,j)*(fy(i,j))2-(Gσ(i,j)*(fx(i,j))2-Gσ(i,j)*(fy(i,j))2)2+4(Gσ(i,j)*fx(i,j)fy(i,j))2},]]>然后根据结构张量矩阵Tσ(i,j)的主特征值λ1(i,j)和次特征值λ2(i,j)，构建如下边缘结构指示符ρ(i,j)：ρ(i,j)＝0.3×ln(|λ1(i,j)‑λ2(i,j)|)其中，|·|表示绝对值，ln是公知的以e为底的对数函数；步骤4：利用边缘结构指示符ρ(i,j)，按照如下方法对噪声图像f的所有像素点(i,j)进行区域分类：其中，V表示区域分类的总类数，且为正整数，取V＝6，Cp(p＝1,2,…,V)表示第p类像素点集，p为区域分类的序号，λmin表示整个图像区域中最小的次特征值，λmax表示整个图像区域中最大的主特征值；步骤5：根据每个像素点所属区域的类别，按照如下方式选取邻域Ni,j的大小R(i,j)：步骤6：构建像素点(i,j)的邻域Ni,j和任意像素点(k,l)的邻域Nk,l的相似性度量函数w(i,j,k,l)为：w(i,j,k,l)=exp{-||Nk,l-Ni,j||22h2}]]>其中，(k,l)表示一个与(i,j)不同的像素点的坐标，表示Nk,l‑Ni,j的L2范数，h表示权值函数的衰减因子，取h＝1；步骤7：根据相似性度量函数w(i,j,k,l)，按照如下方式，从噪声图像f中筛选出与像素点(i,j)相似度最高的S个点，取S＝20：Bi,j‾={(k1,l1),(k2,l2)···,(kS,lS)|w(i,j,k1,l1)≥w(i,j,k2,l2)≥···≥w(i,j,kS,lS)}]]>其中表示与像素点(i,j)相似度最高的S个像素点的集合；步骤8：对噪声图像f去噪，具体方法为：对于噪声图像f中的每一个像素点(i,j)，根据下式去噪，得到像素点(i,j)去噪后的像素值f‾(i,j)=1Q(i,j)Σ(k,l)∈Bi,j‾w(i,j,k,l)f(k,l)]]>其中，Q(i,j)为(i,j)处的归一化常数，Q(i,j)的表达式为Q(i,j)=Σ(k,l)∈Bi,j‾w(i,j,k,l)]]>得到所有像素点去噪后的像素值后，即构成噪声图像f去噪后的图像...

【技术特征摘要】
1.一种权值窗口自适应非局部均值图像去噪方法，其特征在于，该方法包括以下步骤: 步骤1:读入一帧大小SM1XM2X 3的噪声彩色图像Utl，其中M1和M2为正整数，分别表示图像矩阵的行数和列数，然后将输入的噪声彩色图像从RGB彩色空间转换到YCbCr彩色空间，转换后的噪...

【专利技术属性】
技术研发人员：曾维理，季赛平，李聪，路小波，费树岷，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏;32