CN107248148B - 图像降噪方法及系统 - Google Patents

图像降噪方法及系统 Download PDF

Info

Publication number
CN107248148B
CN107248148B CN201710450067.9A CN201710450067A CN107248148B CN 107248148 B CN107248148 B CN 107248148B CN 201710450067 A CN201710450067 A CN 201710450067A CN 107248148 B CN107248148 B CN 107248148B
Authority
CN
China
Prior art keywords
image
edge
contour
points
original image
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
Application number
CN201710450067.9A
Other languages
English (en)
Other versions
CN107248148A (zh
Inventor
邵科
马伟剑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SmartSens Technology Shanghai Co Ltd
Original Assignee
SmartSens Technology Shanghai Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by SmartSens Technology Shanghai Co Ltd filed Critical SmartSens Technology Shanghai Co Ltd
Priority to CN201710450067.9A priority Critical patent/CN107248148B/zh
Publication of CN107248148A publication Critical patent/CN107248148A/zh
Application granted granted Critical
Publication of CN107248148B publication Critical patent/CN107248148B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G06T5/70
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details
    • G06T2207/20032Median filtering
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Abstract

本发明提供了一种图像降噪方法及系统,在对原始图像进行L次滤波,分别获得L个平滑滤波图像,再分别对L个平滑滤波图像进行边缘提取,从而分别获得原始图像的L个边缘图像,并且在L个边缘图像的基础上获得原始图像的轮廓图像,在以上的过程中,可以在较好的提取原始图像的轮廓细节的同时排除噪声的干扰,从而获得质量较佳的轮廓图像。其中L为大于等于2的自然数,L次平滑滤波的平滑系数各不相同,且第L次平滑滤波的平滑系数最大。对第L平滑滤波图像进行中值滤波,获得取值第L平滑滤波图像,根据原始图像、去噪第L平滑滤波图像以及轮廓图像即可获得原始图像降噪的图像。上述方法和系统不仅适用于静态图像处理,还适用于实时图像处理。

Description

图像降噪方法及系统
技术领域
本发明涉及图像处理领域,尤其是一种图像降噪方法及系统。
背景技术
随着电子技术的不断发展,数字相机或者具有数码摄像头的手机等已经深入千家万户,成为大家记录美丽自然风光、铭刻生活以及精彩瞬间的必备工具。
噪点,主要是指数字相机设备的感光组件将光线作为接收信号接收并输出的过程中所产生的图像中的粗糙部分,也指图像中不该出现的外来像素,通常由电子干扰产生,看起来就像图像被弄脏了,布满一些细小的噪点。
目前,为了进一步提高拍照图像的画质,各类拍照设备均会使用一定的降噪算法以滤除不必要的拍照噪点。当用数字相机或手机拍摄光线不是特别充足的场景,比如“夜景”、或者“路灯”时,因为硬件本身的原因,使用较高的增益后图像会出现明显的噪点,这时就需要先用硬件或者软件降噪算法做降噪处理,然后再生成照片。而采用不同的硬件、软件算法、以及不同的参数来降噪,结果会有很大差异。但是在无论采用何种方法对其进行降噪处理,在去除噪点和保留图像的细节二者之间总是无法兼得。
发明内容
本发明的目的在于提供一种图像降噪方法及系统,以实现在平滑噪声的同时尽可能的保留细节。
为了达到上述目的,本发明提供了一种图像降噪方法,包括以下步骤:
对一原始图像进行L次平滑滤波,分别获得L个平滑滤波图像,并分别提取所述L个平滑滤波图像的边缘,获得L个边缘图像,其中L为大于等于2的自然数,所述L次平滑滤波的平滑系数各不相同,且第L次平滑滤波的平滑系数最大;
去除第L平滑滤波图像中的孤立噪点,获得一去噪第L平滑滤波图像;
根据所述L个边缘图像获得所述原始图像的轮廓图像;以及
根据所述原始图像、去噪第L平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像。
优选的,在上述的图像降噪方法中,根据所述L个边缘图像获得所述原始图像的轮廓图像的步骤包括:
将所述L个边缘图像二值化,获得二值化的L个边缘图像;以及
基于二值化的第L边缘图像进行轮廓搜索。
优选的,在上述的图像降噪方法中,根据所述L个边缘图像获得所述原始图像的轮廓图像的步骤还包括:
当二值化的第L边缘图像中断时,查找其他(L-1)个边缘图像中对应位置是否连续,当其他(L-1)个边缘图像中对应位置连续时,则在所述二值化的第L边缘图像中的中断处添加一点,然后对所述二值化的第L边缘图像继续进行轮廓搜索;否则停止所述二值化的第L边缘图像在该处的轮廓搜索,并从下一个有效点开始对所述二值化的第L边缘图像继续进行轮廓搜索。
优选的,在上述的图像降噪方法中,所述有效点是指属于所述二值化的第L边缘图像的边缘轮廓的点。
优选的,在上述的图像降噪方法中,通过中值滤波去除第L平滑滤波图像中的孤立噪点。
优选的,在上述的图像降噪方法中,根据所述原始图像、去噪第L平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像的步骤包括:
对于所述原始图像中的非边缘区域中的点的像素值,采用所述去噪第L平滑滤波图像中对应点的像素值;
对于所述原始图像中的边缘区域中的点的像素值,采用所述原始图像中对应点的像素值;
对于所述原始图像中的部分边缘区域的点的像素值,由以下公式获得:
像素值=原始图像中对应位置的点的像素值*a+去噪第L平滑滤波图像中对应位置的点的像素值*(1-a);
其中,a为权重,其取值为0-1。
优选的,在上述的图像降噪方法中,对于所述原始图像中的任意一M*M区域,如果该M*M区域中不包含任何所述轮廓图像中的轮廓点,则该M*M区域为非边缘区域;
如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以上,则该M*M区域为边缘区域;
如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以内,则该M*M区域为部分边缘区域;
M的取值为大于等于3的自然数,k的取值为30-50。
优选的,在上述的图像降噪方法中,所述L次平滑滤波的平滑系数的取值范围为0-1。
本发明还提供了一种图像降噪系统,包括:L个平滑滤波器、L个边缘提取器、孤立噪声滤波器、轮廓提取器以及加权融合器;
其中,所述L个平滑滤波器用于对一原始图像进行L次平滑滤波,分别获得L个平滑滤波图像;
所述L个边缘提取器分别提取所述L个平滑滤波图像的边缘,获得L个边缘图像,其中L为大于等于2的自然数,所述L次平滑滤波的平滑系数各不相同,且第L次平滑滤波的平滑系数最大;
所述孤立噪声滤波器用于去除第L平滑滤波图像中的孤立噪点,获得一去噪第L平滑滤波图像;
所述轮廓提取器用于根据所述L个边缘图像获得所述原始图像的轮廓图像;
所述加权融合器用于根据所述原始图像、去噪第L平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像。
优选的,在上述的图像降噪系统中,所述轮廓提取器根据所述L个边缘图像获得所述原始图像的轮廓图像的步骤包括:
将所述L个边缘图像二值化,获得二值化的L个边缘图像;以及
基于二值化的第L边缘图像进行轮廓搜索;
当二值化的第L边缘图像中断时,查找其他(L-1)个边缘图像中对应位置是否连续,当其他(L-1)个边缘图像中对应位置连续时,则在所述二值化的第L边缘图像中的中断处添加一点,然后对所述二值化的第L边缘图像继续进行轮廓搜索;否则停止所述二值化的第L边缘图像在该处的轮廓搜索,并从下一个有效点开始对所述二值化的第L边缘图像继续进行轮廓搜索。
优选的,在上述的图像降噪系统中,所述加权融合器用于根据所述原始图像、去噪第L平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像的步骤包括:
对于所述原始图像中的非边缘区域中的点的像素值,采用所述去噪第L平滑滤波图像中对应点的像素值;
对于所述原始图像中的边缘区域中的点的像素值,采用所述原始图像中对应点的像素值;
对于所述原始图像中的部分边缘区域的点的像素值,由以下公式获得:
像素值=原始图像中对应位置的点的像素值*a+去噪第L平滑滤波图像中对应位置的点的像素值*(1-a);
其中,a为权重,其取值为0-1。
优选的,在上述的图像降噪系统中,对于所述原始图像中的任意一M*M区域,如果该M*M区域中不包含任何所述轮廓图像中的轮廓点,则该M*M区域为非边缘区域;
如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以上,则该M*M区域为边缘区域;
如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以内,则该M*M区域为部分边缘区域;
M的取值为大于等于3的自然数,k的取值为30-50。
在本发明提供的图像降噪方法及系统中,在对原始图像进行L次滤波,分别获得L个平滑滤波图像,再分别对所述L个平滑滤波图像进行边缘提取,从而分别获得所述原始图像的L个边缘图像,并且在所述L个边缘图像的基础上获得所述原始图像的轮廓图像,在以上的过程中,可以在较好的提取所述原始图像的轮廓细节的同时排除噪声的干扰,从而获得质量较佳的轮廓图像。其中L为大于等于2的自然数,所述L次平滑滤波的平滑系数各不相同,且第L次平滑滤波的平滑系数最大。对第L平滑滤波图像进行中值滤波,获得取值第L平滑滤波图像,根据原始图像、去噪第L平滑滤波图像以及轮廓图像即可获得所述原始图像降噪的图像。上述方法和系统不仅适用于静态图像处理,还适用于实时图像处理。
附图说明
图1为本发明一实施例中图像降噪方法的示意图;
图2为图1中步骤S3的流程示意图;
图3为本发明一实施例中图像降噪系统的示意图。
具体实施方式
下面将结合示意图对本发明的具体实施方式进行更详细的描述。根据下列描述和权利要求书,本发明的优点和特征将更清楚。需说明的是,附图均采用非常简化的形式且均使用非精准的比例,仅用以方便、明晰地辅助说明本发明实施例的目的。
本发明一实施例提供了一种图像降噪方法,包括以下步骤:对一原始图像进行L次平滑滤波,分别获得L个平滑滤波图像,并分别提取所述L个平滑滤波图像的边缘,获得L个边缘图像,其中L为大于等于2的自然数,所述L次平滑滤波的平滑系数各不相同,且第L次平滑滤波的平滑系数最大;去除第L平滑滤波图像中的孤立噪点,获得一去噪第L平滑滤波图像;根据所述L个边缘图像获得所述原始图像的轮廓图像;以及根据所述原始图像、去噪第L平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像。
在本发明一实施例中,L取值为2。具体的,如图1所示,图1为本发明一实施例中图像降噪方法的示意图。所述图像降噪方法包括以下步骤:对原始图像进行第一平滑系数的滤波,获得第一平滑滤波图像,然后再提取第一平滑滤波的边缘,从而得到第一边缘图像,如图1中的步骤S1。对所述原始图像进行第二平滑系数的滤波,获得第二平滑滤波图像,并提取所述第二平滑滤波图像的边缘,从而获得第二边缘图像,如图1中的步骤S2。基于所述第一边缘图像和第二边缘图像获取所述原始图像的轮廓图像,如图1中的步骤S3。对所述第二平滑滤波图像进行滤波以消除所述第二平滑滤波图像中未被平滑的孤立噪点,从而获得一去噪第二平滑滤波图像,如图1中的步骤S4。最后,再将所述原始图像与所述去噪第二平滑滤波图像基于所述轮廓图像进行加权融合,从而获得降噪后的图像,如图1中的步骤S5。
具体的,在上述步骤S1中,所述第一平滑系数取值范围为大于等于0且小于等于1,其中,当所述第一平滑系数取值为0时,则对所述原始图像的平滑滤波效果最弱,当所述第一平滑系数取值为1时,对所述原始图像的平滑率效果最强。
较佳的,所述第一平滑系数的取值范围为0-0.5,即所述第一平滑系数的取值范围为大于等于0且小于等于0.5,可以为0至0.5之间的任意值,例如,可以为0.1,0.2,0.3,0.4以及0.5。
在对所述原始图像进行第一平滑系数的平滑滤波后,获得所述第一平滑滤波图像,再提取所述第一平滑滤波图像的边缘,从而获得所述第一边缘图像。由于所述第一平滑系数的取值范围为0-0.5,使得所述第一边缘图像在消除所述原始图像中一部分噪点的同时,尽量多的保留所述原始图像的边缘细节。
进一步的,在上述步骤S2中,所述第二平滑系数的取值范围为大于等于0且小于等于1,其中,当所述第二平滑系数取值为0时,则对所述原始图像的平滑滤波效果最弱,当所述第二平滑系数取值为1时,对所述原始图像的平滑率效果最强。
较佳的,所述第二平滑系数的取值范围为0.5-1,即所述第二平滑系数的取值范围为大于等于0.5且小于等于1,可以为0.5至1之间的任意值,例如,可以为0.5,0.6,0.7,0.8,0.9以及1。
在对所述原始图像进行第二平滑系数的平滑滤波后,获得所述第二平滑滤波图像,然后提取所述第二平滑滤波图像的边缘,从而获得所述第二边缘图像。由于所述第二平滑系数的取值范围为0.5-1,从而使得所述第二边缘图像较所述第一边缘图像而言,消除所述原始图像中更多的噪点,使得所述第二边缘图像的非边缘区域更平滑,同时所述第二边缘图像也将损失一些边缘细节。
需要说明的是,在上述步骤S1和步骤S2中,所述第一平滑系数的取值小于所述第二平滑系数的取值,并且,在本发明的实施例中,上述步骤S1和步骤S2的顺序并不做任何限定,可以根据需要任意设定,只需要保证所述第一平滑系数小于所述第二平滑系数即可。
接下来,如图1中的步骤S3,在根据上述步骤S1和步骤S2获取所述第一边缘图像和第二边缘图像后,基于所述第一边缘图像和第二边缘图像提取所述原始图像的轮廓,从而获得所述原始图像的轮廓图像。
具体的,如图2所示,图2为图1中步骤S3的流程示意图。图1中的步骤S3包括以下步骤,首先将所述第一边缘图像和第二边缘图像分别进行二值化,进而获得二值化的第一边缘图像和二值化的第二边缘图像,如图2中的步骤S31。
基于所述二值化的第二边缘图像进行轮廓搜索,如图2中的步骤S32。当所述二值化的第二边缘图像中断时,查找所述二值化的第一边缘图像中的对应位置是否连续,如图2中的步骤S33。
当所述二值化的第一边缘图像中的对应位置连续时,执行步骤S34,否则执行步骤S35。
步骤S34,在所述二值化的第二边缘图像中添加该点,使得所述二值化的第二边缘图像在该中断处连续,并重复上述步骤S32至步骤S33,对所述二值化的第二边缘图像继续进行轮廓搜索。直至获得所述原始图像的轮廓图像。
步骤S35,当所述二值化的第一边缘图像中的对应位置中断时,则停止所述二值化的第二边缘图像在该处的轮廓搜索,并从下一个有效点开始,重复上述步骤S32至步骤S33,对所述二值化的第二边缘图像继续进行轮廓搜索。直至获得所述原始图像的轮廓图像。
所述有效点是指属于所述二值化的第二边缘图像的边缘轮廓的点,具体的判断方法为,从一点开始,可以找到连续的n个点,则表明该点为有效点,否则,则认为该点为噪点。具体的,所述连续的n个点是指任意方向连续的n个点,并不仅限于沿着轮廓连续的n个点,其中,n的取值为大于等于1的自然数,较佳的,n的取值范围为1-5,即可以取1至5中的任意数字,例如可以为2,3,4以及5。当然,在本发明的其他实施例中,所述n的取值还可以为大于5的其他值,本发明并不以此为限。
在上述步骤S4中,对上述步骤S2获得的所述第二平滑滤波图像进行中值滤波,获得所述去噪第二平滑滤波图像,以去除所述第二平滑滤波图像中的孤立噪点。所述孤立噪点是指与四周点的值有明显差异的点。对其该点及其四周点进行中值排序后,判断该点与中间点的差值,当该差值大于一定值时则用中间值将其替换,由于该点一般为最大值点或最小值,在对其进行中值滤波后,可以将其有效的滤除,对于平滑部分的点,当前点与中间点的差值不大,不做替换。
在上述步骤S5中,根据所述原始图像、所述去噪第二平滑滤波图像以及所述原始图像的轮廓图像获得降噪后的图像。
具体的,将原始图像划分为三种区域,第一种为非边缘区域,第二种为边缘区域,第三种为部分边缘区域。具体的,对于所述原始图像中的任意一M*M区域,如果该M*M区域中不包含任何所述轮廓图像中的轮廓点,则该M*M区域为非边缘区域。如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以上,则该M*M区域为边缘区域。如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以内,则该M*M区域为部分边缘区域。
进一步的,对所述非边缘区域,采用所述去噪第二平滑滤波图像中对应的点的像素值。对于所述边缘区域,采用所述原始图像中对应的点的像素值。对于部分边缘区域,根据所述原始图像以及所述去噪第二平滑滤波图像中对应的点的像素按照一定的比例加权获得。
具体而言,对于所述部分边缘区域的点,其像素值由以下公式获得:
像素值=原始图像中对应位置的点的像素值*a+去噪第二平滑滤波图像中对应位置的点的像素值*(1-a),
其中,a为权重,其取值为0-1。
在本发明一实施例中,M的取值为大于等于3的自然数,较佳的,M的取值为3-10,即M可以取3至10中的任意值,例如,可以取值为3、4、5、6、7、8、9以及10中的任意值。k的取值为30-50,在本发明一实施例中,k的取值可以为30至50中的任意一个值。以M取值为5,k的取值30为例,即当所述轮廓图像中的轮廓点占该5*5区域中总点数的30%以上,则该5*5区域为边缘区域。如果该5*5区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该5*5区域中总点数的30%以内,则该5*5区域为部分边缘区域。
通过上述步骤S5的加权融合,即可获得对所述原始图像进行降噪后的图像。在获取所述降噪后的图像的过程中,基于所述原始图像的轮廓图像,将所述原始图像和去噪第二平滑滤波图像进行融合获得,可以在平滑噪声的同时保留细节。并且,上述图像降噪方法不仅适用于静态图像处理,还适用于实时图像处理。
需要说明的是,在本发明的其他实施例中,上述L的取值还可以是其他大于2的值,例如为3、4或者其他值,本发明并不以此为限。当L的取值大于2时,对所述原始图像可以采用更多的平滑滤波系数进行滤波,由由于各个L的值均不相同,从而可以获得更多平滑等级的平滑滤波图像,进而可以更进一步的提高所述原始图像降噪的效果。
本发明实施例还提供了一种图像降噪系统,具体的,所述系统包括:L个平滑滤波器、L个边缘提取器、孤立噪声滤波器、轮廓提取器以及加权融合器;其中,所述L个平滑滤波器用于对一原始图像进行L次平滑滤波,分别获得L个平滑滤波图像;所述L个边缘提取器分别提取所述L个平滑滤波图像的边缘,获得L个边缘图像,其中L为大于等于2的自然数,所述L次平滑滤波的平滑系数各不相同,且第L次平滑滤波的平滑系数最大;所述孤立噪声滤波器用于去除第L平滑滤波图像中的孤立噪点,获得一去噪第L平滑滤波图像;所述轮廓提取器用于根据所述L个边缘图像获得所述原始图像的轮廓图像;所述加权融合器用于根据所述原始图像、去噪第L平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像。
在本发明一实施例中,L的取值为2,具体的,如图3所示,图3为本发明一实施例中图像降噪系统的示意图。所述系统包括:第一平滑滤波器、第一边缘提取器、第二平滑滤波器、第二边缘提取器、孤立噪声滤波器、轮廓提取器以及加权融合器。
具体的,所述第一平滑滤波器用于对一原始图像进行第一平滑系数的滤波,从而获得一第一平滑滤波图像;所述第一边缘提取器用于提取所述第一平滑滤波图像的边缘,从而获得所述原始图像的第一边缘图像。所述第二平滑滤波器用于对所述原始图像进行第二平滑系数的滤波,从而获得所述第二平滑滤波图像。所述第二边缘提取器用于提取所述第二平滑滤波图像的边缘,从而获得所述原始图像的第二边缘图像。所述孤立噪声滤波器用于去除所述第二平滑滤波图像中的孤立噪点,从而获得一去噪第二平滑滤波图像,具体的,所述孤立噪声滤波器包括但不限于中值滤波器。所述轮廓提取器用于基于所述第一边缘图像和第二边缘图像提取所述原始图像的轮廓图像。所述加权融合器用于根据所述原始图像、去噪第二平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像。
所述第一平滑滤波器的平滑系数小于所述第二平滑滤波器的平滑系数。具体的,所述第一平滑滤波器的平滑系数为第一平滑系数,所述第二平滑滤波器的平滑系数为第二平滑系数.所述第一平滑系数取值范围为大于等于0且小于等于1,其中,当所述第一平滑系数取值为0时,则对所述原始图像的平滑滤波效果最弱,当所述第一平滑系数取值为1时,对所述原始图像的平滑率效果最强。较佳的,所述第一平滑系数的取值范围为0-0.5,即所述第一平滑系数的取值范围为大于等于0且小于等于0.5,可以为0至0.5之间的任意值,例如,可以为0.1,0.2,0.3,0.4以及0.5。
在对所述原始图像进行第一平滑系数的平滑滤波后,获得所述第一平滑滤波图像,再提取所述第一平滑滤波图像的边缘,从而获得所述第一边缘图像。由于所述第一平滑系数的取值范围为0-0.5,使得所述第一边缘图像在消除所述原始图像中一部分噪点的同时,尽量多的保留所述原始图像的边缘细节。
进一步的,所述第二平滑系数的取值范围为大于等于0且小于等于1,其中,当所述第二平滑系数取值为0时,则对所述原始图像的平滑滤波效果最弱,当所述第二平滑系数取值为1时,对所述原始图像的平滑率效果最强。较佳的,所述第二平滑系数的取值范围为0.5-1,即所述第二平滑系数的取值范围为大于等于0.5且小于等于1,可以为0.5至1之间的任意值,例如,可以为0.5,0.6,0.7,0.8,0.9以及1。
在对所述原始图像进行第二平滑系数的平滑滤波后,获得所述第二平滑滤波图像,然后提取所述第二平滑滤波图像的边缘,从而获得所述第二边缘图像。由于所述第二平滑系数的取值范围为0.5-1,从而使得所述第二边缘图像较所述第一边缘图像而言,消除所述原始图像中更多的噪点,使得所述第二边缘图像的非边缘区域更平滑,同时所述第二边缘图像也将损失一些边缘细节。
所述轮廓提取器基于所述第一边缘图像和第二边缘图像提取所述原始图像的轮廓图像的过程包括:首先将所述第一边缘图像和第二边缘图像分别进行二值化,进而获得二值化的第一边缘图像和二值化的第二边缘图像,如图2中的步骤S31。接着基于所述二值化的第二边缘图像进行轮廓搜索,如图2中的步骤S32。当所述二值化的第二边缘图像中断时,查找所述二值化的第一边缘图像中的对应位置是否连续,如图2中的步骤S33。
当所述二值化的第一边缘图像中的对应位置连续时,执行图2中的步骤S34,否则执行图2中的步骤S35。
步骤S34,在所述二值化的第二边缘图像中添加该点,使得所述二值化的第二边缘图像在该中断处连续,并重复上述步骤S32至步骤S33,对所述二值化的第二边缘图像继续进行轮廓搜索。直至获得所述原始图像的轮廓图像。
步骤S35,当所述二值化的第一边缘图像中的对应位置中断时,则停止所述二值化的第二边缘图像在该处的轮廓搜索,并从下一个有效点开始,重复上述步骤S32至步骤S33,对所述二值化的第二边缘图像继续进行轮廓搜索。直至获得所述原始图像的轮廓图像。
上述步骤S35中的所述有效点是指属于所述二值化的第二边缘图像的边缘轮廓的点,具体的判断方法为,从一点开始,可以找到连续的n个点,则表明该点为有效点,否则,则认为该点为噪点。具体的,所述连续的n个点是指任意方向连续的n个点,并不仅限于沿着轮廓连续的n个点,其中,n的取值为大于等于1的自然数,较佳的,n的取值范围为1-5,即可以取1至5中的任意数字,例如可以为2,3,4以及5。当然,在本发明的其他实施例中,所述n的取值还可以为大于5的其他值,本发明并不以此为限。
具体的,所述孤立噪声滤波器去除所述第二平滑滤波图像中的孤立噪点,从而获得一去噪第二平滑滤波图像的过程包括:去除所述第二平滑滤波图像中的孤立噪点。所述孤立噪点是指与四周点的值有明显差异的点。对其该点及其四周点进行中值排序后,判断该点与中间点的差值,当该差值大于一定值时则用中间值将其替换,由于该点一般为最大值点或最小值,在对其进行中值滤波后,可以将其有效的滤除,对于平滑部分的点,当前点与中间点的差值不大,不做替换。
进一步的,所述加权融合器根据所述原始图像、去噪第二平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像的过程包括:
将原始图像划分为三种区域,第一种为非边缘区域,第二种为边缘区域,第三种为部分边缘区域。具体的,对于所述原始图像中的任意一M*M区域,如果该M*M区域中不包含任何所述轮廓图像中的轮廓点,则该M*M区域为非边缘区域。如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以上,则该M*M区域为边缘区域。如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以内,则该M*M区域为部分边缘区域。
进一步的,对所述非边缘区域,采用所述去噪第二平滑滤波图像中对应的点的像素值。对于所述边缘区域,采用所述原始图像中对应的点的像素值。对于部分边缘区域,根据所述原始图像以及所述去噪第二平滑滤波图像中对应的点的像素按照一定的比例加权获得。
具体而言,对于所述部分边缘区域的点,其像素值由以下公式获得:
像素值=原始图像中对应位置的点的像素值*a+去噪第二平滑滤波图像中对应位置的点的像素值*(1-a),
其中,a为权重,其取值为0-1。
在本发明一实施例中,M的取值为大于等于3的自然数,较佳的,M的取值为3-10,即M可以取3至10中的任意值,例如,可以取值为3、4、5、6、7、8、9以及10中的任意值。k的取值为30-50,在本发明一实施例中,k的取值可以为30至50中的任意一个值。以M取值为5,k的取值30为例,即当所述轮廓图像中的轮廓点占该5*5区域中总点数的30%以上,则该5*5区域为边缘区域。如果该5*5区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该5*5区域中总点数的30%以内,则该5*5区域为部分边缘区域。
在获取所述降噪后的图像的过程中,基于所述原始图像的轮廓图像,将所述原始图像和去噪第二平滑滤波图像进行融合获得,可以在平滑噪声的同时保留细节。并且,上述图像降噪方法不仅适用于静态图像处理,还适用于实时图像处理。
需要说明的是,在本发明的其他实施例中,上述L的取值还可以是其他大于2的值,例如为3、4或者其他值,本发明并不以此为限。当L的取值大于2时,对所述原始图像可以采用更多的平滑滤波系数进行滤波,由由于各个L的值均不相同,从而可以获得更多平滑等级的平滑滤波图像,进而可以更进一步的提高所述原始图像降噪的效果。
综上,在本发明实施例提供的图像降噪方法及系统中,在对原始图像进行L次滤波,分别获得L个平滑滤波图像,再分别对所述L个平滑滤波图像进行边缘提取,从而分别获得所述原始图像的L个边缘图像,并且在所述L个边缘图像的基础上获得所述原始图像的轮廓图像,在以上的过程中,可以在较好的提取所述原始图像的轮廓细节的同时排除噪声的干扰,从而获得质量较佳的轮廓图像。其中L为大于等于2的自然数,所述L次平滑滤波的平滑系数各不相同,且第L次平滑滤波的平滑系数最大。对第L平滑滤波图像进行中值滤波,获得取值第L平滑滤波图像,根据原始图像、去噪第L平滑滤波图像以及轮廓图像即可获得所述原始图像降噪的图像。上述方法和系统不仅适用于静态图像处理,还适用于实时图像处理。
上述仅为本发明的优选实施例而已,并不对本发明起到任何限制作用。任何所属技术领域的技术人员,在不脱离本发明的技术方案的范围内,对本发明揭露的技术方案和技术内容做任何形式的等同替换或修改等变动,均属未脱离本发明的技术方案的内容,仍属于本发明的保护范围之内。

Claims (9)

1.一种图像降噪方法,其特征在于,包括以下步骤:
对一原始图像进行L次平滑滤波,分别获得L个平滑滤波图像,并分别提取所述L个平滑滤波图像的边缘,获得L个边缘图像,其中L为大于等于2的自然数,所述L次平滑滤波的平滑系数各不相同,且第L次平滑滤波的平滑系数最大;
去除第L平滑滤波图像中的孤立噪点,获得一去噪第L平滑滤波图像;
根据所述L个边缘图像获得所述原始图像的轮廓图像;以及
根据所述原始图像、去噪第L平滑滤波图像基于所述原始图像的轮廓图像进行加权融合获得所述原始图像降噪后的图像;
其中,根据所述L个边缘图像获得所述原始图像的轮廓图像的步骤包括:
将所述L个边缘图像二值化,获得二值化的L个边缘图像;
基于二值化的第L边缘图像进行轮廓搜索;以及
当二值化的第L边缘图像中断时,查找其他(L-1)个边缘图像中对应位置是否连续,当其他(L-1)个边缘图像中对应位置连续时,则在所述二值化的第L边缘图像中的中断处添加一点,然后对所述二值化的第L边缘图像继续进行轮廓搜索;否则停止所述二值化的第L边缘图像在该处的轮廓搜索,并从下一个有效点开始对所述二值化的第L边缘图像继续进行轮廓搜索。
2.如权利要求1所述的图像降噪方法,其特征在于,所述有效点是指属于所述二值化的第L边缘图像的边缘轮廓的点。
3.如权利要求1所述的图像降噪方法,其特征在于,通过中值滤波去除第L平滑滤波图像中的孤立噪点。
4.如权利要求1所述的图像降噪方法,其特征在于,根据所述原始图像、去噪第L平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像的步骤包括:
对于所述原始图像中的非边缘区域中的点的像素值,采用所述去噪第L平滑滤波图像中对应点的像素值;
对于所述原始图像中的边缘区域中的点的像素值,采用所述原始图像中对应点的像素值;
对于所述原始图像中的部分边缘区域的点的像素值,由以下公式获得:
像素值=原始图像中对应位置的点的像素值*a+去噪第L平滑滤波图像中对应位置的点的像素值*(1-a);
其中,a为权重,其取值为0-1。
5.如权利要求4所述的图像降噪方法,其特征在于,对于所述原始图像中的任意一M*M区域,如果该M*M区域中不包含任何所述轮廓图像中的轮廓点,则该M*M区域为非边缘区域;
如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以上,则该M*M区域为边缘区域;
如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以内,则该M*M区域为部分边缘区域;
其中,M的取值为大于等于3的自然数,k的取值为30-50。
6.如权利要求1所述的图像降噪方法,其特征在于,所述L次平滑滤波的平滑系数的取值范围为0-1。
7.一种图像降噪系统,其特征在于,包括:L个平滑滤波器、L个边缘提取器、孤立噪声滤波器、轮廓提取器以及加权融合器;
其中,所述L个平滑滤波器用于对一原始图像进行L次平滑滤波,分别获得L个平滑滤波图像;
所述L个边缘提取器分别提取所述L个平滑滤波图像的边缘,获得L个边缘图像,其中L为大于等于2的自然数,所述L次平滑滤波的平滑系数各不相同,且第L次平滑滤波的平滑系数最大;
所述孤立噪声滤波器用于去除第L平滑滤波图像中的孤立噪点,获得一去噪第L平滑滤波图像;
所述轮廓提取器用于根据所述L个边缘图像获得所述原始图像的轮廓图像;
所述加权融合器用于根据所述原始图像、去噪第L平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像;
所述轮廓提取器根据所述L个边缘图像获得所述原始图像的轮廓图像的步骤包括:
将所述L个边缘图像二值化,获得二值化的L个边缘图像;
基于二值化的第L边缘图像进行轮廓搜索;以及
当二值化的第L边缘图像中断时,查找其他(L-1)个边缘图像中对应位置是否连续,当其他(L-1)个边缘图像中对应位置连续时,则在所述二值化的第L边缘图像中的中断处添加一点,然后对所述二值化的第L边缘图像继续进行轮廓搜索;否则停止所述二值化的第L边缘图像在该处的轮廓搜索,并从下一个有效点开始对所述二值化的第L边缘图像继续进行轮廓搜索。
8.如权利要求7所述的图像降噪系统,其特征在于,所述加权融合器用于根据所述原始图像、去噪第L平滑滤波图像以及所述原始图像的轮廓图像获得所述原始图像降噪后的图像的步骤包括:
对于所述原始图像中的非边缘区域中的点的像素值,采用所述去噪第L平滑滤波图像中对应点的像素值;
对于所述原始图像中的边缘区域中的点的像素值,采用所述原始图像中对应点的像素值;
对于所述原始图像中的部分边缘区域的点的像素值,由以下公式获得:
像素值=原始图像中对应位置的点的像素值*a+去噪第L平滑滤波图像中对应位置的点的像素值*(1-a);
其中,a为权重,其取值为0-1。
9.如权利要求8所述的图像降噪系统,其特征在于,对于所述原始图像中的任意一M*M区域,如果该M*M区域中不包含任何所述轮廓图像中的轮廓点,则该M*M区域为非边缘区域;
如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以上,则该M*M区域为边缘区域;
如果该M*M区域中包含所述轮廓图像中的轮廓点,且所述轮廓图像中的轮廓点占该M*M区域中总点数的k%以内,则该M*M区域为部分边缘区域;
其中,M的取值为大于等于3的自然数,k的取值为30-50。
CN201710450067.9A 2017-06-14 2017-06-14 图像降噪方法及系统 Active CN107248148B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710450067.9A CN107248148B (zh) 2017-06-14 2017-06-14 图像降噪方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710450067.9A CN107248148B (zh) 2017-06-14 2017-06-14 图像降噪方法及系统

Publications (2)

Publication Number Publication Date
CN107248148A CN107248148A (zh) 2017-10-13
CN107248148B true CN107248148B (zh) 2023-08-25

Family

ID=60019361

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710450067.9A Active CN107248148B (zh) 2017-06-14 2017-06-14 图像降噪方法及系统

Country Status (1)

Country Link
CN (1) CN107248148B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110264480B (zh) * 2018-03-12 2023-07-07 株式会社理光 一种图像分割的方法、装置及计算机可读存储介质
CN108391067B (zh) * 2018-05-10 2023-09-29 浙江芯劢微电子股份有限公司 一种基于rgb-ir传感器的去噪增强装置及方法
CN109639929B (zh) * 2019-01-11 2020-10-30 珠海全志科技股份有限公司 图像降噪方法、计算机装置及计算机可读存储介质
CN110189348B (zh) * 2019-05-29 2020-12-25 北京达佳互联信息技术有限公司 头像处理方法、装置、计算机设备及存储介质
CN110880165A (zh) * 2019-10-15 2020-03-13 杭州电子科技大学 一种基于轮廓和颜色特征融合编码的图像去雾方法
CN112188178B (zh) * 2020-09-30 2022-02-22 Tcl华星光电技术有限公司 图像显示方法和图像显示装置
CN113850800B (zh) * 2021-10-15 2024-04-30 郑州磨料磨具磨削研究所有限公司 一种硬脆材料划切缝崩边检测方法
CN116596802B (zh) * 2023-06-08 2023-10-24 常州市宏发纵横新材料科技股份有限公司 一种复杂背景下的捆绑纱提取方法、装置、设备及介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7292733B2 (en) * 2002-10-16 2007-11-06 Matsushita Electric Industrial Co., Ltd. Image processing apparatus and image processing method
JP2007306501A (ja) * 2006-05-15 2007-11-22 Fujifilm Corp 画像処理方法、画像処理装置、および画像処理プログラム
CN102184523A (zh) * 2011-04-13 2011-09-14 展讯通信(上海)有限公司 基于多分辨率和过完备变换的数字图像降噪方法及系统
CN102456221A (zh) * 2010-10-25 2012-05-16 新奥特(北京)视频技术有限公司 一种图像噪声快速去除方法
CN102479383A (zh) * 2010-11-30 2012-05-30 上海银晨智能识别科技有限公司 椒盐噪声去除方法及装置
CN102999916A (zh) * 2012-12-12 2013-03-27 清华大学深圳研究生院 一种彩色图像的边缘提取方法
CN103208104A (zh) * 2013-04-16 2013-07-17 浙江工业大学 一种基于非局部理论的图像去噪方法
CN104077767A (zh) * 2013-03-29 2014-10-01 北京千橡网景科技发展有限公司 检测名片边缘的方法和装置
CN104299191A (zh) * 2014-09-28 2015-01-21 南通大学 一种超声图像去噪方法及中值滤波方法
CN104835127A (zh) * 2015-05-19 2015-08-12 中国农业科学院农业信息研究所 一种自适应平滑滤波方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8849057B2 (en) * 2011-05-19 2014-09-30 Foveon, Inc. Methods for digital image sharpening with noise amplification avoidance

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7292733B2 (en) * 2002-10-16 2007-11-06 Matsushita Electric Industrial Co., Ltd. Image processing apparatus and image processing method
JP2007306501A (ja) * 2006-05-15 2007-11-22 Fujifilm Corp 画像処理方法、画像処理装置、および画像処理プログラム
CN102456221A (zh) * 2010-10-25 2012-05-16 新奥特(北京)视频技术有限公司 一种图像噪声快速去除方法
CN102479383A (zh) * 2010-11-30 2012-05-30 上海银晨智能识别科技有限公司 椒盐噪声去除方法及装置
CN102184523A (zh) * 2011-04-13 2011-09-14 展讯通信(上海)有限公司 基于多分辨率和过完备变换的数字图像降噪方法及系统
CN102999916A (zh) * 2012-12-12 2013-03-27 清华大学深圳研究生院 一种彩色图像的边缘提取方法
CN104077767A (zh) * 2013-03-29 2014-10-01 北京千橡网景科技发展有限公司 检测名片边缘的方法和装置
CN103208104A (zh) * 2013-04-16 2013-07-17 浙江工业大学 一种基于非局部理论的图像去噪方法
CN104299191A (zh) * 2014-09-28 2015-01-21 南通大学 一种超声图像去噪方法及中值滤波方法
CN104835127A (zh) * 2015-05-19 2015-08-12 中国农业科学院农业信息研究所 一种自适应平滑滤波方法

Also Published As

Publication number Publication date
CN107248148A (zh) 2017-10-13

Similar Documents

Publication Publication Date Title
CN107248148B (zh) 图像降噪方法及系统
Zhang et al. Enhancement and noise reduction of very low light level images
EP2852152B1 (en) Image processing method, apparatus and shooting terminal
CN106251318B (zh) 一种序列图像的去噪装置及方法
WO2016206087A1 (zh) 一种低照度图像处理方法和装置
WO2015184408A1 (en) Scene motion correction in fused image systems
CN105894460A (zh) 图像滤波方法及装置
CN103606132A (zh) 基于空域和时域联合滤波的多帧数字图像去噪方法
CN109389560B (zh) 一种自适应加权滤波图像降噪方法、装置及图像处理设备
KR20130077726A (ko) 디지털 사진에서 노이즈를 제거하는 장치 및 방법
KR20100016278A (ko) 신호에 대한 희망하는 효과를 달성하기 위하여 신호에 커널을 적용하는 효과를 조절하기 위한 기술
CN109584198B (zh) 一种人脸图像质量评价方法、装置及计算机可读存储介质
CN111161167B (zh) 基于中通道补偿和自适应大气光估计的单幅图像去雾方法
KR20130031573A (ko) 영상 처리 방법 및 영상 처리 장치
WO2023273868A1 (zh) 图像降噪方法、装置、终端、存储介质
Yu et al. Image and video dehazing using view-based cluster segmentation
CN101562692A (zh) 噪声降低装置及方法
CN105338221B (zh) 一种图像处理方法及电子设备
CN110264415B (zh) 一种消除抖动模糊的图像处理方法
JP2008112349A (ja) 移動物体検知装置
KR101651889B1 (ko) 극저조도 영상의 화질 개선 장치 및 방법
KR101527962B1 (ko) 비디오 영상의 움직임 객체 추출 방법
CN110705492A (zh) 一种舞台移动机器人障碍目标检测方法
Xiong et al. Improved haze removal algorithm using dark channel prior
KR100647955B1 (ko) 디포커스 현상을 보정하는 영상처리장치

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 200233 A1101, 900 Yishan Road, Xuhui District, Shanghai

Applicant after: Siteway (Shanghai) Electronic Technology Co.,Ltd.

Address before: 200233 A1101, 900 Yishan Road, Xuhui District, Shanghai

Applicant before: SHANGHAI YE XIN ELECTRONIC TECHNOLOGY CO.,LTD.

CB02 Change of applicant information
CB02 Change of applicant information

Address after: Room 612, 6th floor, building 3, 111 Xiangke Road, China (Shanghai) pilot Free Trade Zone, Pudong New Area, Shanghai 200131

Applicant after: Starway (Shanghai) Electronic Technology Co.,Ltd.

Address before: 200233 A1101, 900 Yishan Road, Xuhui District, Shanghai

Applicant before: Siteway (Shanghai) Electronic Technology Co.,Ltd.

GR01 Patent grant
GR01 Patent grant