CN105654445B - 一种基于小波变换边缘检测的手机图像去噪方法 - Google Patents
一种基于小波变换边缘检测的手机图像去噪方法 Download PDFInfo
- Publication number
- CN105654445B CN105654445B CN201610060892.3A CN201610060892A CN105654445B CN 105654445 B CN105654445 B CN 105654445B CN 201610060892 A CN201610060892 A CN 201610060892A CN 105654445 B CN105654445 B CN 105654445B
- Authority
- CN
- China
- Prior art keywords
- image
- edge
- value
- noise
- point
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000003708 edge detection Methods 0.000 title claims abstract description 16
- 230000009466 transformation Effects 0.000 title abstract description 4
- 238000009792 diffusion process Methods 0.000 claims abstract description 47
- 238000009499 grossing Methods 0.000 claims abstract description 20
- 238000001914 filtration Methods 0.000 claims abstract description 14
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 230000011218 segmentation Effects 0.000 claims abstract description 4
- 238000010276 construction Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 4
- 230000001131 transforming effect Effects 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 6
- 238000001514 detection method Methods 0.000 abstract description 2
- 238000003706 image smoothing Methods 0.000 abstract description 2
- 230000000717 retained effect Effects 0.000 abstract 1
- 230000000007 visual effect Effects 0.000 description 3
- 239000006002 Pepper Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 235000016761 Piper aduncum Nutrition 0.000 description 1
- 235000017804 Piper guineense Nutrition 0.000 description 1
- 244000203593 Piper nigrum Species 0.000 description 1
- 235000008184 Piper nigrum Nutrition 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- 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/10—Image acquisition modality
- G06T2207/10024—Color image
-
- 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/20004—Adaptive image processing
-
- 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/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Picture Signal Circuits (AREA)
Abstract
本发明公开了一种基于小波变换边缘检测的手机图像去噪方法,能够在有效平滑噪声的同时保留图像边缘细节,去除彩色噪声,从而得到理想的去噪效果,包括如下步骤:第一步,读入彩色噪声图像,将输入的图像从RGB彩色空间转换到YUV彩色空间;第二步,用高斯滤波器平滑图像,对Y通道图像进行三次低通滤波;第三步,利用小波变换算法检测图像边缘;第四步,采用阈值分割得到边缘的二值化图像;第五步,在图像边缘区域采用各向异性扩散平滑噪声;第六步,在图像平滑区域采用均值滤波去噪;第七步,对平滑区域图像进一步去除彩色噪声;第八步,将图像从YUV彩色空间转换到RGB彩色空间,得到最终的去噪图像。
Description
技术领域
本发明属于图像处理领域,涉及一种基于小波变换边缘检测的手机图像去噪方法。
背景技术
手机拍照由于其便捷性广泛应用于生活实际中,人们对手机图像要求越来越高。在图像处理领域,去除噪声一直是研究热点。手机图像去噪过程中易损害边缘信息,为了提高手机图像的质量,在去除噪声的同时又不能使图像的边缘轮廓变模糊,对图像边缘处和平滑区域采用不同的方法去噪。现有的边缘检测算子主要检测信号的高频部分,受噪声影响严重,并不能检测到图像的最佳边缘。基于小波变换的边缘检测算法通过构造小波系数,分离信号和噪声,用检测模极大值方法可有效检测出图像的边缘。
由于图像边缘往往携带大部分的信息,对图像边缘去噪的同时也要保留图像的细节。现有的去噪算法主要适用于特定类型的噪声,如均值滤波主要用于去除高斯噪声,中值滤波主要用于去除椒盐噪声或称脉冲噪声,基于双变量收缩函数的去噪方法主要用于去除高斯噪声,手机图像所含噪声并不属于某种特定类型,且含有彩色噪声。本发明提供一种基于小波变换边缘检测的手机图像去噪方法,通过设置扩散系数,在各个方向上扩散程度根据梯度模的不同而改变,在图像边缘实行较弱的平滑,受噪声类型影响不大。通过减小图像的色度值进一步平滑彩色噪声,从而使处理后的手机图像视觉效果更好。
发明内容
技术问题:本发明提供一种可以有效去除噪声并保留图像边缘的基于小波变换边缘检测的手机图像去噪方法。
技术方案:本发明的一种基于小波变换边缘检测的手机图像去噪方法,包括如下步骤:
步骤1:读入一帧大小为L1×L2×3的噪声彩色图像I0,其中L1和L2为正整数,分别表示噪声彩色图像I0的总行数和总列数,然后将输入的噪声彩色图像I0从RGB彩色空间转换到YUV彩色空间:
其中,R表示RGB彩色空间中的红色分量,G表示RGB彩色空间中的绿色分量,B表示RGB彩色空间中的蓝色分量,Y表示YUV彩色空间中的亮度值,U表示YUV彩色空间中的饱和度,V表示YUV彩色空间中的色调值,转换后的YUV噪声图像记为f,大小为L1×L2×3;
步骤2:用高斯滤波器平滑YUV噪声图像f,对Y通道进行三次低通滤波,得到平滑图像f1:
f1=((fY*h1)*h2)*h3
其中,fY表示YUV噪声图像f在Y通道的分量,*表示卷积,h1、h2、h3为低通滤波器,分别为:
h1=[0.125 0.375 0.375 0.125],
h2=[0.125 0 0.375 0 0.375 0 0.125],
h3=[0.125 0 0 0 0.375 0 0 0 0.375 0 0 0 0.125]
步骤3:对平滑图像f1,用小波变换算法检测图像边缘,具体方法如下:
由平滑图像f1和离散化后的小波函数μx和μy进行卷积,得到平滑图像f1沿水平和垂直方向的梯度Mxf1和Myf1:
Mxf1=f1*μx,
Myf1=f1*μy
其中,x和y分别表示水平和垂直方向坐标轴,μx和μy分别为离散化后水平和垂直方向的小波函数,构建方法如下:
构建两个二维小波函数μx(x,y)和μy(x,y):
其中,θ(x,y)是二维高斯平滑函数,μx(x,y)和μy(x,y)分别表示二维高斯平滑函数θ(x,y)沿x和y方向的一阶偏导数,离散化后水平和垂直方向的小波函数μx和μy分别为:
根据平滑图像f1沿水平和垂直方向的梯度Mxf1和Myf1构建图像梯度的模Cf1和幅角Af1:
其中,|Mxf1|表示Mxf1的绝对值,|Myf1|表示Myf1的绝对值,表示的反正切函数,
沿幅角Af1方向寻找图像梯度的模Cf1的局部极大值点,具体方法为:
其中,i和j分别表示整幅图像的行坐标和列坐标,f2(i,j)表示边缘图像f2在点(i,j)处的值,Cf1(i+1,j)、Cf1(i+1,j+1)、Cf1(i,j+1)、Cf1(i-1,j+1)、Cf1(i-1,j)、Cf1(i-1,j-1)、Cf1(i,j-1)、Cf1(i+1,j-1)分别表示图像梯度的模Cf1在点(i+1,j)、(i+1,j+1)、(i,j+1)、(i-1,j+1)、(i-1,j)、(i-1,j-1)、(i,j-1)、(i+1,j-1)处的值,Af1(i,j)表示幅角Af1在点(i,j)处的值,并以此找出所有图像梯度的模Cf1的局部极大值点,构成边缘图像f2;
步骤4:用阈值分割得到边缘的二值化图像,将YUV噪声图像f分为边缘区域和平滑区域,其具体方法为:
将小波变换后的边缘图像f2拉伸到0~255区间,获得对比度增强后的图像f3:
其中,f2(i,j)表示边缘图像f2在点(i,j)处的值,f3(i,j)表示对比度增强后的图像f3在点(i,j)处的值,f2max和f2min分别为边缘图像f2的灰度最大值和最小值,
设置阈值T,比T大的像素点为255,比T小的像素点为0,得到二值边缘图像f4:
其中,f4(i,j)表示二值边缘图像f4在点(i,j)处的值,阈值T根据图像灰度的不同选取不同的值:
其中,fave表示YUV噪声图像f所有像素点的灰度均值,由所有二值边缘图像f4中灰度值为255的像素点所对应的YUV噪声图像f中的区域构成边缘区域含噪图像P,由所有二值边缘图像f4中灰度值为0的像素点所对应的YUV噪声图像f中的区域构成平滑区域含噪图像Q;
步骤5:采用各向异性扩散方法对边缘区域含噪图像P去噪,具体算法如下:
根据各向异性扩散模型,按照如下方式,构建拉普拉斯算子8方向离散化的各向异性扩散方程:
Pt+1(i1,j1)=Pt(i1,j1)+λ[rN·▽NP+rS·▽SP+rW·▽WP+rE·▽EP+rNE·▽NEP+rSE·▽SEP+rNW·▽NWP+rSW·▽SWP]t(i1,j1)其中,i1和j1分别表示边缘区域图像的行坐标和列坐标,t表示迭代次数,设置为20次,Pt(i1,j1)表示边缘区域含噪图像P在点(i1,j1)处第t次迭代后的各向异性扩散去噪值,迭代初始值为边缘区域含噪图像P在Y通道的分量在点(i1,j1)处的值,即P0(i1,j1)=PY(i1,j1),λ表示控制扩散总体强度的常数,取λ=1/4,r表示扩散系数,其选取遵循图像梯度越大,扩散系数越小,图像梯度越小,扩散系数越大的原则,故扩散系数r的表达式为:
其中,D表示扩散门限,取值为10,▽表示梯度算子,F(▽P)表示关于▽P的函数,▽P表示边缘区域含噪图像P的梯度,在各个方向的分量为:
▽NP(i1,j1)=P(i1-1,j1)-P(i1,j1),▽SP(i1,j1)=P(i1+1,j1)-P(i1,j1),
▽EP(i1,j1)=P(i1,j1+1)-P(i1,j1),▽WP(i1,j1)=P(i1,j1-1)-P(i1,j1),
▽NEP(i1,j1)=P(i1-1,j1+1)-P(i1,j1),▽SEP(i1,j1)=P(i1+1,j1+1)-P(i1,j1),
▽SWP(i1,j1)=P(i1+1,j1-1)-P(i1,j1),▽NWP(i1,j1)=P(i1-1,j1-1)-P(i1,j1)
其中,N、S、E、W分别表示上、下、右、左方向,NE、SE、NW、SW则分别表示右上、右下、左上、左下方向,P(i1-1,j1)、P(i1+1,j1)、P(i1,j1+1)、P(i1,j1-1)、P(i1-1,j1+1)、P(i1+1,j1+1)、P(i1+1,j1-1)、P(i1-1,j1-1)分别表示边缘区域含噪图像P在点(i1-1,j1、(i1+1,j1)、(i1,j1+1)、(i1,j1-1)、(i1-1,j1+1)、(i1+1,j1+1)、(i1+1,j1-1)、(i1-1,j1-1)处的值,▽NP(i1,j1)、▽SP(i1,j1)、▽EP(i1,j1)、▽WP(i1,j1)、▽NEP(i1,j1)、▽SEP(i1,j1)、▽SWP(i1,j1)、▽NWP(i1,j1)分别表示边缘区域含噪图像P在点(i1,j1)处上、下、右、左、右上、右下、左上、左下方向上的梯度,
扩散系数r在每次迭代中作为边缘区域含噪图像的梯度▽P的函数被更新:
其中, 分别表示扩散系数r在点(i1,j1)处第t次迭代时上、下、右、左、右上、右下、左上、左下方向上的分量,分别对应边缘区域含噪图像P在点(i1,j1)处第t次迭代时上、下、右、左、右上、右下、左上、左下方向上的梯度▽NPt(i1,j1)、▽SPt(i1,j1)、▽EPt(i1,j1)、▽WPt(i1,j1)、▽NEPt(i1,j1)、▽SEPt(i1,j1)、▽SWPt(i1,j1)、▽NWPt(i1,j1)的函数,
经过一定迭代次数,自适应平滑噪声,并保留梯度强的图像边缘;
步骤6:采用均值滤波对平滑区域含噪图像Q进行去噪,得到平滑区域均值滤波后的图像H:
其中,i2和j2分别表示平滑区域图像的行坐标和列坐标,Q(i2,j2)表示平滑区域含噪图像Q在点(i2,j2)处的值,H(i2,j2)表示平滑区域均值滤波后的图像H在点(i2,j2)处的值,Zi2,j2表示中心在点(i2,j2)处、大小为m×n的矩形子图像邻域,取m=n=17个像素点;
步骤7:寻找出平滑区域均值滤波后的图像H的彩色噪声点,再对彩色噪声点进行平滑处理,其中,彩色噪声点的寻找方法为:
如果满足以下条件,则点(i2,j2)为彩色噪声点:
其中,δ表示判断彩色噪声的阈值,取δ=30,HU(i2,j2)和HV(i2,j2)分别表示平滑区域均值滤波后的图像H在点(i2,j2)处的U值和V值,即饱和度和色调值,
对彩色噪声进行平滑处理的具体方法为:缩小HU(i2,j2)和HV(i2,j2)的值,得到平滑区域均值滤波后的图像H在点(i2,j2)处的缩小后的饱和度HU2(i2,j2)和色调值HV2(i2,j2):
HU2(i2,j2)=HU(i2,j2)/k,
HV2(i2,j2)=HV(i2,j2)/k
其中,k表示阈值,根据平滑区域均值滤波后的图像H在点(i2,j2)处亮度值HY(i2,j2)的不同选取不同的值:
步骤8:将经步骤7处理后的彩色图像从YUV彩色空间转换到RGB彩色空间:
并将转换后的图像输出。
有益效果:与现有方法相比,本发明具有如下优点:
1、由于图片受噪声干扰严重,现有的边缘检测算子对噪声较为敏感,易把噪声较强的点误判为边缘,损失细节信息。为了精确定位边缘,本发明利用小波变换检测局域突变的能力,通过构造水平和垂直方向的小波函数,与低通滤波后的图像卷积,取得平滑噪声后的图像梯度;根据图像边缘梯度大,噪声点梯度小,通过检测图像梯度的模极大值检测局域突变,减小噪声点的误判;而弱边缘的模值比较小,为了防止边缘丢失,沿图像梯度的幅角多个方向,在像素点邻域通过判断幅角所在范围检测图像梯度的模极大值,对于弱边缘也能在局部区域根据边缘的方向性被检测出来,且多个方向的划分有利于细化边缘,弱化噪声的影响。针对不同灰度值的图片设置合适的阈值定位边缘,可抑制噪声的干扰,提高了边缘检测的精度,从而有效区分图像边缘处和平滑区域。
2、为了在去噪的同时避免边缘模糊,本发明在图像边缘处采用各向异性扩散去噪,设置自适应扩散系数,在各个方向上扩散程度根据梯度模的不同而改变,噪声点梯度低,扩散系数大,平滑效果越强,而在梯度高的边缘区域,扩散系数小,从而实现在图像边缘实行较弱的平滑,保留图像边缘细节。由于待去噪的图像较大,手机图像去噪实时性要求高,本发明在各向异性扩散基础上增加梯度的方向性,有利于提高单次迭代的扩散效果,从而减少扩散的迭代次数,在保证去噪效果的同时提高了运算效率,便于应用实现。
3、现有的图像去噪方法对于固定类型的噪声去噪效果好,而手机图像受拍摄条件如光线、拍摄对象等因素的影响,图像噪声类型特殊,噪声点呈大片块状,图像亮度值越低噪声点越多,且彩色噪声严重,噪声分布不规律,随机性大,不属于高斯噪声、椒盐噪声等现有的噪声类型。本发明在图像平滑区域采用均值滤波可平滑大片块状噪声,再通过减小图像的饱和度和色调值,且根据亮度值的不同缩减程度不同,进一步平滑彩色噪声,从而使处理后的手机图像去噪效果理想,视觉效果更好。
附图说明
图1为一种基于小波变换边缘检测的手机图像去噪方法流程图。
图2为手机拍摄的彩色噪声图像。
图3为用一种基于小波变换边缘检测的手机图像去噪方法去噪后的图像。
具体实施方式
下面结合说明书附图和实施例,清楚、完整地描述本发明方法的详细过程。
步骤1:读入一帧大小为L1×L2×3的噪声彩色图像I0,其中L1和L2为正整数,分别表示噪声彩色图像I0的总行数和总列数,然后将输入的噪声彩色图像I0从RGB彩色空间转换到YUV彩色空间:
其中,R表示RGB彩色空间中的红色分量,G表示RGB彩色空间中的绿色分量,B表示RGB彩色空间中的蓝色分量,Y表示YUV彩色空间中的亮度值,U表示YUV彩色空间中的饱和度,V表示YUV彩色空间中的色调值,转换后的YUV噪声图像记为f,大小为L1×L2×3;
步骤2:用高斯滤波器平滑YUV噪声图像f,对Y通道进行三次低通滤波,得到平滑图像f1:
f1=((fY*h1)*h2)*h3
其中,fY表示YUV噪声图像f在Y通道的分量,*表示卷积,h1、h2、h3为低通滤波器,分别为:
h1=[0.125 0.375 0.375 0.125],
h2=[0.125 0 0.375 0 0.375 0 0.125],
h3=[0.125 0 0 0 0.375 0 0 0 0.375 0 0 0 0.125]
步骤3:对平滑图像f1,用小波变换算法检测图像边缘,具体方法如下:
由平滑图像f1和离散化后的小波函数μx和μy进行卷积,得到平滑图像f1沿水平和垂直方向的梯度Mxf1和Myf1:
Mxf1=f1*μx,
Myf1=f1*μy
其中,x和y分别表示水平和垂直方向坐标轴,μx和μy分别为离散化后水平和垂直方向的小波函数,构建方法如下:
构建两个二维小波函数μx(x,y)和μy(x,y):
其中,θ(x,y)是二维高斯平滑函数,μx(x,y)和μy(x,y)分别表示二维高斯平滑函数θ(x,y)沿x和y方向的一阶偏导数,取x={-2,-1,0,1,2},y={-2,-1,0,1,2},离散化后水平和垂直方向的小波函数μx和μy分别为:
根据平滑图像f1沿水平和垂直方向的梯度Mxf1和Myf1构建图像梯度的模Cf1和幅角Af1:
其中,|Mxf1|表示Mxf1的绝对值,|Myf1|表示Myf1的绝对值,表示的反正切函数,
沿幅角Af1方向寻找图像梯度的模Cf1的局部极大值点,具体方法为:
其中,i和j分别表示整幅图像的行坐标和列坐标,f2(i,j)表示边缘图像f2在点(i,j)处的值,Cf1(i+1,j)、Cf1(i+1,j+1)、Cf1(i,j+1)、Cf1(i-1,j+1)、Cf1(i-1,j)、Cf1(i-1,j-1)、Cf1(i,j-1)、Cf1(i+1,j-1)分别表示图像梯度的模Cf1在点(i+1,j)、(i+1,j+1)、(i,j+1)、(i-1,j+1)、(i-1,j)、(i-1,j-1)、(i,j-1)、(i+1,j-1)处的值,Af1(i,j)表示幅角Af1在点(i,j)处的值,并以此找出所有图像梯度的模Cf1的局部极大值点,构成边缘图像f2;
步骤4:用阈值分割得到边缘的二值化图像,将YUV噪声图像f分为边缘区域和平滑区域,其具体方法为:
将小波变换后的边缘图像f2拉伸到0~255区间,获得对比度增强后的图像f3:
其中,f2(i,j)表示边缘图像f2在点(i,j)处的值,f3(i,j)表示对比度增强后的图像f3在点(i,j)处的值,f2max和f2min分别为边缘图像f2的灰度最大值和最小值,
设置阈值T,比T大的像素点为255,比T小的像素点为0,得到二值边缘图像f4:
其中,f4(i,j)表示二值边缘图像f4在点(i,j)处的值,阈值T根据图像灰度的不同选取不同的值:
其中,fave表示YUV噪声图像f所有像素点的灰度均值,由于光照比较弱的图片噪声相对较大,设置比较大的阈值能够滤除大部分噪声,将不同图像按灰度均值的不同分为四类,取不同的阈值T,由所有二值边缘图像f4中灰度值为255的像素点所对应的YUV噪声图像f中的区域构成边缘区域含噪图像P,由所有二值边缘图像f4中灰度值为0的像素点所对应的YUV噪声图像f中的区域构成平滑区域含噪图像Q;
步骤5:采用各向异性扩散方法对边缘区域含噪图像P去噪,具体算法如下:
根据各向异性扩散模型,按照如下方式,构建拉普拉斯算子8方向离散化的各向异性扩散方程:
Pt+1(i1,j1)=Pt(i1,j1)+λ[rN·▽NP+rS·▽SP+rW·▽WP+rE·▽EP+rNE·▽NEP+rSE·▽SEP+rNW·▽NWP+rSW·▽SWP]t(i1,j1)其中,i1和j1分别表示边缘区域图像的行坐标和列坐标,t表示迭代次数,设置为20次,Pt(i1,j1)表示边缘区域含噪图像P在点(i1,j1)处第t次迭代后的各向异性扩散去噪值,迭代初始值为边缘区域含噪图像P在Y通道的分量在点(i1,j1)处的值,即P0(i1,j1)=PY(i1,j1),λ表示控制扩散总体强度的常数,取λ=1/4,r表示扩散系数,其选取遵循图像梯度越大,扩散系数越小,图像梯度越小,扩散系数越大的原则,故扩散系数r的表达式为:
其中,D表示扩散门限,取值为10,▽表示梯度算子,F(▽P)表示关于▽P的函数,▽P表示边缘区域含噪图像P的梯度,在各个方向的分量为:
▽NP(i1,j1)=P(i1-1,j1)-P(i1,j1),▽SP(i1,j1)=P(i1+1,j1)-P(i1,j1),
▽EP(i1,j1)=P(i1,j1+1)-P(i1,j1),▽WP(i1,j1)=P(i1,j1-1)-P(i1,j1),
▽NEP(i1,j1)=P(i1-1,j1+1)-P(i1,j1),▽SEP(i1,j1)=P(i1+1,j1+1)-P(i1,j1),
▽SWP(i1,j1)=P(i1+1,j1-1)-P(i1,j1),▽NWP(i1,j1)=P(i1-1,j1-1)-P(i1,j1)
其中,N、S、E、W分别表示上、下、右、左方向,NE、SE、NW、SW则分别表示右上、右下、左上、左下方向,P(i1-1,j1)、P(i1+1,j1)、P(i1,j1+1)、P(i1,j1-1)、P(i1-1,j1+1)、P(i1+1,j1+1)、P(i1+1,j1-1)、P(i1-1,j1-1)分别表示边缘区域含噪图像P在点(i1-1,j1、(i1+1,j1)、(i1,j1+1)、(i1,j1-1)、(i1-1,j1+1)、(i1+1,j1+1)、(i1+1,j1-1)、(i1-1,j1-1)处的值,▽NP(i1,j1)、▽SP(i1,j1)、▽EP(i1,j1)、▽WP(i1,j1)、▽NEP(i1,j1)、▽SEP(i1,j1)、▽SWP(i1,j1)、▽NWP(i1,j1)分别表示边缘区域含噪图像P在点(i1,j1)处上、下、右、左、右上、右下、左上、左下方向上的梯度,
扩散系数r在每次迭代中作为边缘区域含噪图像的梯度▽P的函数被更新:
其中, 分别表示扩散系数r在点(i1,j1)处第t次迭代时上、下、右、左、右上、右下、左上、左下方向上的分量,分别对应边缘区域含噪图像P在点(i1,j1)处第t次迭代时上、下、右、左、右上、右下、左上、左下方向上的梯度▽NPt(i1,j1)、▽SPt(i1,j1)、▽EPt(i1,j1)、▽WPt(i1,j1)、▽NEPt(i1,j1)、▽SEPt(i1,j1)、▽SWPt(i1,j1)、▽NWPt(i1,j1)的函数,
经过一定迭代次数,自适应平滑噪声,并保留梯度强的图像边缘;
步骤6:采用均值滤波对平滑区域含噪图像Q进行去噪,得到平滑区域均值滤波后的图像H:
其中,i2和j2分别表示平滑区域图像的行坐标和列坐标,Q(i2,j2)表示平滑区域含噪图像Q在点(i2,j2)处的值,H(i2,j2)表示平滑区域均值滤波后的图像H在点(i2,j2)处的值,Zi2,j2表示中心在点(i2,j2)处、大小为m×n的矩形子图像邻域,取m=n=17个像素点;
步骤7:寻找出平滑区域均值滤波后的图像H的彩色噪声点,再对彩色噪声点进行平滑处理,其中,彩色噪声点的寻找方法为:
考虑到RGB空间和YUV空间存在彩色空间转换,则RGB空间下,如果满足以下条件,则点(i2,j2)为彩色噪声点:
YUV空间下,如果满足以下条件,则点(i2,j2)为彩色噪声点:
其中,δ表示判断彩色噪声的阈值,取δ=30,HR(i2,j2)、HG(i2,j2)、HB(i2,j2)分别表示平滑区域均值滤波后的图像H在点(i2,j2)处的R、G、B值,Have(i2,j2)表示HB(i2,j2)值的均值,即Have(i2,j2)=[HR(i2,j2)+HG(i2,j2)+HB(i2,j2)]/3,HU(i2,j2)和HV(i2,j2)分别表示平滑区域均值滤波后的图像H在点(i2,j2)处的U值和V值,即饱和度和色调值,
对彩色噪声进行平滑处理的具体方法为:缩小HU(i2,j2)和HV(i2,j2)的值,得到平滑区域均值滤波后的图像H在点(i2,j2)处的缩小后的饱和度HU2(i2,j2)和色调值HV2(i2,j2):
HU2(i2,j2)=HU(i2,j2)/k,
HV2(i2,j2)=HV(i2,j2)/k
其中,k表示阈值,由于亮度暗的区域噪声较大,当平滑区域均值滤波后的图像H在点(i2,j2)处亮度值HY(i2,j2)较小时,对应的k稍大,而HY(i2,j2)较大时对应的k稍小,根据HY(i2,j2)的不同k选取不同的值:
步骤8:将经步骤7处理后的彩色图像从YUV彩色空间转换到RGB彩色空间:
并将转换后的图像输出。
一种基于小波变换边缘检测的手机图像去噪的应用实验:
在应用实验中,对大量实际拍摄的噪声图像进行了测试。图2中上图是一张大小为4128×3096的手机拍摄的彩色噪声图像,下图是由图2上图红框标注区域放大后的图像,图3中上图是利用本发明提出的一种基于小波变换边缘检测的手机图像去噪方法处理后的图像,下图是由图3上图红框标注区域放大后的图像,可以看出,图2中的图像含有大量噪声,彩色噪声十分明显,尤其在边缘区域受噪声干扰严重,利用本发明的方法可以避免把噪声较强的点误判为边缘,并保留梯度小的弱图像边缘,有效区分图像边缘处和平滑区域;有效去除图像中的噪声,在图像边缘实行较弱的平滑,避免边缘模糊;在图像平滑区域去除较大噪声,使平滑区域的图像更加光滑,并且进一步平滑彩色噪声,从而使处理后的手机图像视觉效果更好。
应理解上述实施例仅用于说明本发明技术方案的具体实施方式,而不用于限制本发明的范围。在阅读了本发明之后,本领域技术人员对本发明的各种等同形式的修改和替换均落于本申请权利要求所限定的保护范围。
Claims (1)
1.一种基于小波变换边缘检测的手机图像去噪方法,其特征在于,按照以下步骤进行:
步骤1:读入一帧大小为L1×L2×3的噪声彩色图像I0,其中L1和L2为正整数,分别表示噪声彩色图像I0的总行数和总列数,然后将输入的噪声彩色图像I0从RGB彩色空间转换到YUV彩色空间:
其中,R表示RGB彩色空间中的红色分量,G表示RGB彩色空间中的绿色分量,B表示RGB彩色空间中的蓝色分量,Y表示YUV彩色空间中的亮度值,U表示YUV彩色空间中的饱和度,V表示YUV彩色空间中的色调值,转换后的YUV噪声图像记为f,大小为L1×L2×3;
步骤2:用高斯滤波器平滑YUV噪声图像f,对Y通道进行三次低通滤波,得到平滑图像f1:
f1=((fY*h1)*h2)*h3
其中,fY表示YUV噪声图像f在Y通道的分量,*表示卷积,h1、h2、h3为低通滤波器,分别为:
h1=[0.125 0.375 0.375 0.125],
h2=[0.125 0 0.375 0 0.375 0 0.125],
h3=[0.125 0 0 0 0.375 0 0 0 0.375 0 0 0 0.125]
步骤3:对平滑图像f1,用小波变换算法检测图像边缘,具体方法如下:
由平滑图像f1和离散化后的小波函数μx和μy进行卷积,得到平滑图像f1沿水平和垂直方向的梯度Mxf1和Myf1:
Mxf1=f1*μx,
Myf1=f1*μy
其中,x和y分别表示水平和垂直方向坐标轴,μx和μy分别为离散化后水平和垂直方向的小波函数,构建方法如下:
构建两个二维小波函数μx(x,y)和μy(x,y):
其中,θ(x,y)是二维高斯平滑函数,μx(x,y)和μy(x,y)分别表示二维高斯平滑函数θ(x,y)沿x和y方向的一阶偏导数,离散化后水平和垂直方向的小波函数μx和μy分别为:
根据平滑图像f1沿水平和垂直方向的梯度Mxf1和Myf1构建图像梯度的模Cf1和幅角Af1:
其中,|Mxf1|表示Mxf1的绝对值,|Myf1|表示Myf1的绝对值,表示的反正切函数,
沿幅角Af1方向寻找图像梯度的模Cf1的局部极大值点,具体方法为:
其中,i和j分别表示整幅图像的行坐标和列坐标,f2(i,j)表示边缘图像f2在点(i,j)处的值,Cf1(i+1,j)、Cf1(i+1,j+1)、Cf1(i,j+1)、Cf1(i-1,j+1)、Cf1(i-1,j)、Cf1(i-1,j-1)、Cf1(i,j-1)、Cf1(i+1,j-1)分别表示图像梯度的模Cf1在点(i+1,j)、(i+1,j+1)、(i,j+1)、(i-1,j+1)、(i-1,j)、(i-1,j-1)、(i,j-1)、(i+1,j-1)处的值,Af1(i,j)表示幅角Af1在点(i,j)处的值,并以此找出所有图像梯度的模Cf1的局部极大值点,构成边缘图像f2;
步骤4:用阈值分割得到边缘的二值化图像,将YUV噪声图像f分为边缘区域和平滑区域,其具体方法为:
将小波变换后的边缘图像f2拉伸到0~255区间,获得对比度增强后的图像f3:
其中,f2(i,j)表示边缘图像f2在点(i,j)处的值,f3(i,j)表示对比度增强后的图像f3在点(i,j)处的值,f2max和f2min分别为边缘图像f2的灰度最大值和最小值,
设置阈值T,比T大的像素点为255,比T小的像素点为0,得到二值边缘图像f4:
其中,f4(i,j)表示二值边缘图像f4在点(i,j)处的值,阈值T根据图像灰度的不同选取不同的值:
其中,fave表示YUV噪声图像f所有像素点的灰度均值,由所有二值边缘图像f4中灰度值为255的像素点所对应的YUV噪声图像f中的区域构成边缘区域含噪图像P,由所有二值边缘图像f4中灰度值为0的像素点所对应的YUV噪声图像f中的区域构成平滑区域含噪图像Q;
步骤5:采用各向异性扩散方法对边缘区域含噪图像P去噪,具体算法如下:
根据各向异性扩散模型,按照如下方式,构建拉普拉斯算子8方向离散化的各向异性扩散方程:
其中,i1和j1分别表示边缘区域图像的行坐标和列坐标,t表示迭代次数,设置为20次,Pt(i1,j1)表示边缘区域含噪图像P在点(i1,j1)处第t次迭代后的各向异性扩散去噪值,迭代初始值为边缘区域含噪图像P在Y通道的分量在点(i1,j1)处的值,即P0(i1,j1)=PY(i1,j1),λ表示控制扩散总体强度的常数,取λ=1/4,r表示扩散系数,其选取遵循图像梯度越大,扩散系数越小,图像梯度越小,扩散系数越大的原则,故扩散系数r的表达式为:
其中,D表示扩散门限,取值为10,表示梯度算子,表示关于的函数,表示边缘区域含噪图像P的梯度,在各个方向的分量为:
其中,N、S、E、W分别表示上、下、右、左方向,NE、SE、NW、SW则分别表示右上、右下、左上、左下方向,P(i1-1,j1)、P(i1+1,j1)、P(i1,j1+1)、P(i1,j1-1)、P(i1-1,j1+1)、P(i1+1,j1+1)、P(i1+1,j1-1)、P(i1-1,j1-1)分别表示边缘区域含噪图像P在点(i1-1,j 1)、(i1+1,j1)、(i1,j1+1)、(i1,j1-1)、(i1-1,j1+1)、(i1+1,j1+1)、(i1+1,j1-1)、(i1-1,j1-1)处的值, 分别表示边缘区域含噪图像P在点(i1,j1)处上、下、右、左、右上、右下、左下、左上方向上的梯度,
扩散系数r在每次迭代中作为边缘区域含噪图像的梯度的函数被更新:
其中, 分别表示扩散系数r在点(i1,j1)处第t次迭代时上、下、右、左、右上、右下、左上、左下方向上的分量,分别对应边缘区域含噪图像P在点(i1,j1)处第t次迭代时上、下、右、左、右上、右下、左上、左下方向上的梯度 的函数,
经过一定迭代次数,自适应平滑噪声,并保留梯度强的图像边缘;
步骤6:采用均值滤波对平滑区域含噪图像Q进行去噪,得到平滑区域均值滤波后的图像H:
其中,i2和j2分别表示平滑区域图像的行坐标和列坐标,Q(i2,j2)表示平滑区域含噪图像Q在点(i2,j2)处的值,H(i2,j2)表示平滑区域均值滤波后的图像H在点(i2,j2)处的值,Zi2,j2表示中心在点(i2,j2)处、大小为m×n的矩形子图像邻域,取m=n=17个像素点;
步骤7:寻找出平滑区域均值滤波后的图像H的彩色噪声点,再对彩色噪声点进行平滑处理,其中,彩色噪声点的寻找方法为:
如果满足以下条件,则点(i2,j2)为彩色噪声点:
其中,δ表示判断彩色噪声的阈值,取δ=30,HU(i2,j2)和HV(i2,j2)分别表示平滑区域均值滤波后的图像H在点(i2,j2)处的U值和V值,即饱和度和色调值,
对彩色噪声进行平滑处理的具体方法为:缩小HU(i2,j2)和HV(i2,j2)的值,得到平滑区域均值滤波后的图像H在点(i2,j2)处的缩小后的饱和度HU2(i2,j2)和色调值HV2(i2,j2):
HU2(i2,j2)=HU(i2,j2)/k,
HV2(i2,j2)=HV(i2,j2)/k
其中,k表示阈值,根据平滑区域均值滤波后的图像H在点(i2,j2)处亮度值HY(i2,j2)的不同选取不同的值:
步骤8:将经步骤7处理后的彩色图像从YUV彩色空间转换到RGB彩色空间:
并将转换后的图像输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610060892.3A CN105654445B (zh) | 2016-01-28 | 2016-01-28 | 一种基于小波变换边缘检测的手机图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610060892.3A CN105654445B (zh) | 2016-01-28 | 2016-01-28 | 一种基于小波变换边缘检测的手机图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105654445A CN105654445A (zh) | 2016-06-08 |
CN105654445B true CN105654445B (zh) | 2018-08-21 |
Family
ID=56488978
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610060892.3A Active CN105654445B (zh) | 2016-01-28 | 2016-01-28 | 一种基于小波变换边缘检测的手机图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105654445B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106604057A (zh) * | 2016-12-07 | 2017-04-26 | 乐视控股(北京)有限公司 | 视频处理方法及装置 |
CN108428215A (zh) * | 2017-02-15 | 2018-08-21 | 阿里巴巴集团控股有限公司 | 一种图像处理方法、装置和设备 |
CN107203009B (zh) * | 2017-05-17 | 2019-05-14 | 杰创智能科技股份有限公司 | 一种基于小波特征向量提取的手机检测方法 |
JP7003558B2 (ja) * | 2017-10-12 | 2022-01-20 | カシオ計算機株式会社 | 画像処理装置、画像処理方法、及びプログラム |
CN107784639B (zh) * | 2017-11-02 | 2021-03-09 | 长安大学 | 一种无人机遥感图像改进多边滤波去噪方法 |
CN110246087B (zh) * | 2018-03-07 | 2021-06-04 | 舜宇光学(浙江)研究院有限公司 | 参考多通道的多分辨率去除图像色度噪声的系统和方法 |
CN108513044B (zh) * | 2018-04-16 | 2020-11-13 | 深圳市华星光电技术有限公司 | 图像平滑处理方法、电子装置及计算机可读存储介质 |
CN108985349B (zh) * | 2018-06-26 | 2021-09-28 | 福州大学 | 一种考虑图像集复杂性的石材纹理图像cnn识别方法 |
CN109492653B (zh) * | 2018-11-15 | 2022-04-12 | 深圳市比邻星精密技术有限公司 | 乳腺病灶体积测量方法、装置、计算机设备及存储介质 |
CN109801235B (zh) * | 2018-12-28 | 2023-03-28 | 佛山科学技术学院 | 一种绿萝叶片病害原因检测方法及装置 |
CN112070717B (zh) * | 2020-08-05 | 2024-06-04 | 煜邦数字科技(广东)有限公司 | 基于图像处理的输电线路覆冰厚度检测方法 |
CN112918956A (zh) * | 2021-02-20 | 2021-06-08 | 陆伟凤 | 一种基于图像识别技术的垃圾分类系统 |
CN116124754B (zh) * | 2023-04-17 | 2023-07-25 | 山东省中源联科生物工程集团有限公司 | 通过图像分析检测透明质酸钠凝胶中交联剂残留量的方法 |
CN116993629B (zh) * | 2023-09-27 | 2023-12-19 | 福建晟哲自动化科技有限公司 | 基于图像分解的平滑方法、装置、电子设备及存储介质 |
CN118570708A (zh) * | 2024-08-02 | 2024-08-30 | 淄博思创信息科技有限公司 | 一种视频智能分析方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102222328A (zh) * | 2011-07-01 | 2011-10-19 | 杭州电子科技大学 | 一种边缘保持的自然场景图像自适应加权滤波方法 |
CN104331863A (zh) * | 2014-10-15 | 2015-02-04 | 中国航空工业集团公司洛阳电光设备研究所 | 一种图像滤波去噪方法 |
CN104766278A (zh) * | 2015-03-19 | 2015-07-08 | 天津大学 | 基于自适应平滑系数的各向异性滤波方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7747045B2 (en) * | 2006-06-30 | 2010-06-29 | Fujifilm Corporation | Method and apparatus for diffusion based illumination normalization |
-
2016
- 2016-01-28 CN CN201610060892.3A patent/CN105654445B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102222328A (zh) * | 2011-07-01 | 2011-10-19 | 杭州电子科技大学 | 一种边缘保持的自然场景图像自适应加权滤波方法 |
CN104331863A (zh) * | 2014-10-15 | 2015-02-04 | 中国航空工业集团公司洛阳电光设备研究所 | 一种图像滤波去噪方法 |
CN104766278A (zh) * | 2015-03-19 | 2015-07-08 | 天津大学 | 基于自适应平滑系数的各向异性滤波方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105654445A (zh) | 2016-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105654445B (zh) | 一种基于小波变换边缘检测的手机图像去噪方法 | |
US9692991B2 (en) | Multispectral imaging system | |
KR100872253B1 (ko) | 이미지 센서에 의해 생성된 이미지의 노이즈 제거 방법 | |
WO2013168618A1 (ja) | 画像処理装置及び画像処理方法 | |
CN102222328B (zh) | 一种边缘保持的自然场景图像自适应加权滤波方法 | |
CN111784605B (zh) | 一种基于区域指导的图像降噪方法、计算机装置及计算机可读存储介质 | |
KR20150117646A (ko) | 적어도 하나의 추가적인 이미지를 사용한 이미지 향상과 에지 검증을 위한 방법 및 장치 | |
CN107403435B (zh) | 一种基于rgb颜色空间的彩色图像边缘提取方法 | |
CN103347190B (zh) | 一种边相关的联合彩色去马赛克和放大方法 | |
CN116823686B (zh) | 一种基于图像增强的夜间红外与可见光图像融合方法 | |
CN112583999B (zh) | 相机模块的镜头脏污检测方法 | |
US9153013B2 (en) | Image processing apparatus, image processing method and computer readable medium | |
CN110517206B (zh) | 彩色摩尔纹消除方法及装置 | |
CN111510691A (zh) | 颜色插值方法及装置、设备、存储介质 | |
KR20110005408A (ko) | 영상 잡음 제거 장치 및 방법 | |
Wang et al. | An efficient method for image dehazing | |
CN113706521A (zh) | 碳纤维表面毛团检测方法、装置、存储介质和电子设备 | |
KR20140109801A (ko) | 3d이미지 품질을 향상시키는 방법과 장치 | |
Wang et al. | Single-image dehazing using color attenuation prior based on haze-lines | |
CN104766287A (zh) | 一种基于显著性检测的模糊图像盲复原方法 | |
CN110136085B (zh) | 一种图像的降噪方法及装置 | |
CN115841632A (zh) | 输电线路提取方法、装置以及双目测距方法 | |
Xie et al. | Local average fitting active contour model with thresholding for noisy image segmentation | |
Anandh et al. | Visual Enhancement of Underwater Images Using Transmission Estimation and Multi-Scale Fusion. | |
Liu et al. | Infrared ship target image smoothing based on adaptive mean shift |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |