CN107392869B - 一种基于边缘保持滤波器的人脸图像滤波方法 - Google Patents

一种基于边缘保持滤波器的人脸图像滤波方法 Download PDF

Info

Publication number
CN107392869B
CN107392869B CN201710602140.XA CN201710602140A CN107392869B CN 107392869 B CN107392869 B CN 107392869B CN 201710602140 A CN201710602140 A CN 201710602140A CN 107392869 B CN107392869 B CN 107392869B
Authority
CN
China
Prior art keywords
face image
image
filtering
filter
color space
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
CN201710602140.XA
Other languages
English (en)
Other versions
CN107392869A (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.)
Xi'an Huizhi Big Data Technology Co ltd
Original Assignee
Changan University
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 Changan University filed Critical Changan University
Priority to CN201710602140.XA priority Critical patent/CN107392869B/zh
Publication of CN107392869A publication Critical patent/CN107392869A/zh
Application granted granted Critical
Publication of CN107392869B publication Critical patent/CN107392869B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10024Color image
    • 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/20172Image enhancement details
    • G06T2207/20192Edge enhancement; Edge preservation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种基于边缘保持滤波器的人脸滤波算法,首先获取待滤波的人脸图像;接着将获取到的人脸图像进行压缩,减少计算量;然后将RGB颜色空间的人脸图像转换到Lab颜色空间;再对得到的Lab颜色空间的人脸图像按照本发明提出的滤波算法进行滤波处理;而后将处理后的Lab颜色空间的人脸图像转换到RGB颜色空间;最后将图像进行放大并输出处理后的人脸图像。本发明能够不仅能满足人脸图像滤波较高的滤波要求,又有能效地保持图像的边缘信息,同时算法复杂度较低,计算效率较高,更加具有实用性,对后续提取人脸图像中的有用信息有着积极的意义。

Description

一种基于边缘保持滤波器的人脸图像滤波方法
技术领域
本发明属于图像处理技术领域,具体涉及一种基于边缘保持滤波器的人脸图像滤波方法。
背景技术
数字图像处理是指将图像信号转换成数字信号并利用计算机对其进行处理的过程。早期的图像处理的目的是改善图像的质量,以获得图像中的大量信息。图像处理的常用方法有图像增强、复原、编码、压缩等。
近年来,随着图像处理技术的不断发展,其应用范围也越来越广,深入到人们生活的方方面面,如人脸图像处理技术就是其中重要的一部分。人脸图像处理技术能够显著提高图像质量,不仅能够应用到各种美图应用软件,在人脸识别、人脸跟踪等技术中有着重要的应用价值。分析人脸部图像特征可以发现,人脸图像中含有大量的细节信息如眼睛、鼻子的轮廓、发际线的纹理,同时又含有大量的待去除的“噪声”,如脸部的斑点等无关内容。在对人脸图像进行处理时,既要保证良好的滤波效果,又要能够保留住脸部器官的边缘、轮廓等重要信息,可见人脸滤波算法是一项非常具有挑战的工作。
目前,人脸图像滤波处理有很多方法,其中空间域平滑是数字图像处理中较为常用且有效的方法。经过几十年的发展,空域滤波出现了很多比较成熟的算法,如均值滤波、中值滤波、维纳滤波。但是在本质上,空域滤波算法都属于低通滤波,由于低通滤波后的图像丢失了一些有用的高频信息,致使滤波后人脸图像中细节丢失,滤波后的图像虽然去除了噪声,但使得边缘模糊,对于要求较高的人脸图像处理来说这是难以接受的。所以,各种去噪方法的研究实际是在去噪和保留高频信息之间进行的权衡。
对于上述问题,有人提出应用双边滤波。它不仅采取传统滤波方法的策略,考虑了空间位置信息,还加入了值域相似度的影响,因此,在图像变化平缓的区域,邻域内像素亮度值相差不大,双边滤波转化为高斯低通滤波器;在图像变化剧烈的区域,滤波器利用边缘点附近亮度值相近的像素点的亮度值平均替代原亮度值。因此,双边滤波不仅可以消除图像噪声而且保留了图像边缘信息。双边滤波虽然有上述优点,但是它存在效率低下,滤波效果不佳的缺点。考虑到双边滤波效率低下的问题,目前人脸滤波中比较常见的还有应用引导滤波算法。引导滤波可以构建出与滤波半径无关的滤波模版,大幅度提高滤波速率,但是其又存在着边缘保持效果不佳的问题。
发明内容
本发明的目的在于提供一种基于边缘保持滤波器的人脸图像滤波方法,以克服上述现有技术存在的缺陷,本发明既满足了人脸图像对滤波质量较高的要求,又有效地保持了图像中的边缘等细节信息,同时算法的计算复杂度较低,实用性较强,具有极高的应用价值。
为达到上述目的,本发明采用如下技术方案:
一种基于边缘保持滤波器的人脸图像滤波方法,包括以下步骤:
步骤1:获取待处理的人脸图像;
步骤2:将步骤1得到的人脸图像进行压缩;
步骤3:将步骤2得到的压缩人脸图像从RGB颜色空间转换到Lab颜色空间;
步骤4:将步骤3处理后的人脸图像进行滤波处理;
步骤5:将步骤4中所得的滤波后的人脸图像转换到RGB颜色空间;
步骤6:将步骤5中得到图像进行放大,输出最终人脸图像。
进一步地,步骤2采用了双线性插值的方法对步骤1得到的人脸图像进行压缩,其中压缩倍率为原图的1/N。
进一步地,步骤3中将步骤2得到的压缩人脸图像从RGB颜色空间转换到Lab颜色空间,即压缩人脸图像由R、G、B颜色分量转变为L、a、b三个分量,其中L表示图像的亮度,a表示从红色到绿色的范围,b表示从黄色到蓝色的范围。
进一步地,步骤4采用边缘保持滤波器对步骤3处理后的人脸图像进行滤波处理,边缘保持滤波器公式表述为:
Figure BDA0001357408980000031
其中,Ii为特征图像中的某一像素点,且特征图像为步骤3处理后的人脸图像转换得到的灰度图像,qi为滤波器输出图像中相应的像素点,ak和bk为滤波窗口wk中的一组线性系数,k为滤波半径。
进一步地,提出新的线性系数ak的计算公式为:
Figure BDA0001357408980000032
Figure BDA0001357408980000033
其中,ε为调整系数,pi为步骤3处理后的人脸图像中的某一像素点,Gx(pi)和Gy(pi)分别为pi在水平和竖直方向上与临近像素点的亮度差值,Gradk(pi)为pi对应的梯度值。
进一步地,调整系数ε为0.5。
进一步地,线性系数bk的计算公式为:
Figure BDA0001357408980000034
其中,
Figure BDA0001357408980000035
为步骤3处理后的人脸图像在滤波窗口wk内像素点的均值,
Figure BDA0001357408980000036
为特征图像在滤波框wk内像素点的均值。
进一步地,
Figure BDA0001357408980000037
Figure BDA0001357408980000038
均采用包滤波器计算得到。
进一步地,滤波窗口wk的大小为5×5。
进一步地,步骤6中采用双线性插值法将步骤5中得到图像进行放大,放大倍率为N倍。
与现有技术相比,本发明具有以下有益的技术效果:
本发明方法既满足了人脸图像对滤波质量较高的要求,又有效地保持了图像中的边缘等细节信息,同时算法的计算复杂度较低,实用性较强,具有极高的应用价值。
进一步地,本发明采用新的边缘保持滤波器,采用待滤波人脸图像的灰度图作为特征图像,与待滤波人脸图像一并作为输入,新构造的滤波器在滤波前采用双线性插值法对图像进行压缩,减少了图像中的冗余信息。
进一步地,系数ak、bk计算时能够使用包滤波器(box-filter),能够构建滤波复杂度与滤波半径无关滤波器,极大提高了算法的计算效率。
进一步地,本发明中系数ak的计算结合了待滤波图像和特征图像中梯度信息,对图像中的边缘细节更加敏感。与传统的双边滤波算法相比,新的边缘保持滤波方法算法复杂度较低,计算效率更高,滤波和边缘保持效果更好。
附图说明
图1是本发明的流程示意图;
图2(a)、图3(a)、图4(a)和图5(a)均为原始人脸图像;
图2(b)、图3(b)、图4(b)和图5(b)均为采用双边滤波后的人脸图像;
图2(c)、图3(c)、图4(c)和图5(c)均为采用本发明算法滤波后的人脸图像;
图6(a)、图7(a)、图8(a)和图9(a)分别为图2(a)、图3(a)、图4(a)和图5(a)的局部图;
图6(b)、图7(b)、图8(b)和图9(b)分别为图2(b)、图3(b)、图4(b)和图5(b)的局部图;
图6(c)、图7(c)、图8(c)和图9(c)分别为图2(c)、图3(c)、图4(c)和图5(c)的局部图;
图10为图2(a)在第300行像素点灰度值的变化图;
图11为图2(b)在第300行像素点灰度值的变化图;
图12为图2(c)在第300行像素点灰度值的变化图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
参见图1,本发明提供了一种基于新的边缘保持滤波器的人脸滤波方法。
步骤1、获取人脸图像:利用图像采集设备,获取要滤波处理的人脸图像q。
步骤2、计算步骤1中获取的人脸图像q的尺寸m×n,按照双线性插值法,取原始图像临近的四个像素点的值加权平均得到压缩后的图像,即压缩后图像的大小为原图大小的1/4。
步骤3、将步骤2中压缩后的人脸图像转换到Lab颜色空间,图像由R、G、B颜色分量转变为L、a、b三个分量。
步骤4、根据新构建的边缘保持滤波器,按照公式
Figure BDA0001357408980000051
对人脸图像进行滤波处理。取待滤波人脸图像的灰度图作为特征图像,Ii为其中的某一像素点,qi为滤波后输出的像素点,滤波框wk大小为5。公式中滤波系数ak和bk计算方法为:
Figure BDA0001357408980000052
Figure BDA0001357408980000053
Figure BDA0001357408980000054
取ε为0.5,通过计算待滤波人脸图像中某一点像素点pi在水平和竖直方向上与临近像素点的差值Gx(pi)和Gy(pi)得到点pi梯度Gradk(pi),一并作为系数ak计算公式的输入。
Figure BDA0001357408980000055
为待滤波人脸图像在滤波框wk内像素点的均值,
Figure BDA0001357408980000056
为特征图像在滤波框wk内像素点的均值,k为滤波半径,且k=16。
步骤5、将滤波后的Lab颜色空间的图像转换到RGB颜色空间。
步骤6、根据双线性插值法,对步骤5中得到的图像放大4倍后输出。
图2(a)、图2(b)和图2(c),图3(a)、图3(b)和图3(c),图4(a)、图4(b)和图4(c)以及图5(a)、图5(b)和图5(c)分别采用采用四组图像,用双边滤波器作为对照组与本发明中提出的新的边缘保持滤波器的滤波效果进行对比。从实验结果可以看到与原图和双边滤波后的图像相比,本发明方法在较好的滤除了人脸图像中的斑点等噪声后,较好的保持了图像的边缘细节信息,人物面部轮廓比较清晰。
图6(a)、图6(b)和图6(c),图7(a)、图7(b)和图7(c),图8(a)、图8(b)和图8(c)以及图9(a)、图9(b)和图9(c)不同算法滤波前后的局部放大图,从中可以直观地看到处理效果。图6(a)、图7(a)、图8(a)和图9(a)为待滤波的人脸图像,可见图中人物脸部带有大量斑点噪声,分析双边滤波的效果可以发现,经双边滤波后的图6(b)和图9(b)中人物脸部仍有斑点残留,图7(b)和图8(b)中人物头发、眼部比较模糊。对比应用本发明提出的滤波算法可以看出,四组图片经本发明滤波后人脸图像中的斑点噪声基本去除,图像中边缘等细节保持较好,整体滤波效果有较大提升。
分析图10、图11和图12可以看出,滤波前图像中含有较多噪声,双边滤波能够去除部分噪声,但噪声去除还不够彻底,在灰度变换剧烈的区域的边缘保持效果不够好,本发明提出滤波算法,能够滤除绝大多数噪声,同时尽可能的保留图像中的边缘细节信息,可以看出在图中灰度呈直线上升即灰度变化剧烈的区域的特征,滤波后仍然保留。
表1不同算法对不同尺寸的图像处理用时的对比表
Figure BDA0001357408980000061
Figure BDA0001357408980000071
从表1可以看出,本发明提出的算法耗时大大缩短,极大的提高了算法的实用性。因而,本发明提出的基于新的边缘保持滤波器的人脸图像滤波算法优于已有的滤波算法,对进一步处理以及准确获取人脸图像中的特征信息有着重要的意义。

Claims (1)

1.一种基于边缘保持滤波器的人脸图像滤波方法,其特征在于,包括以下步骤:
步骤1:获取待处理的人脸图像;
步骤2:采用双线性插值的方法对步骤1得到的人脸图像进行压缩,其中压缩倍率为原图的1/N,N为正整数;
步骤3:将步骤2得到的压缩人脸图像从RGB颜色空间转换到Lab颜色空间,即压缩人脸图像由R、G、B颜色分量转变为L、a、b三个分量,其中L表示图像的亮度,a表示从红色到绿色的范围,b表示从黄色到蓝色的范围;
步骤4:将步骤3处理后的人脸图像进行滤波处理;具体采用边缘保持滤波器对步骤3处理后的人脸图像进行滤波处理,边缘保持滤波器公式表述为:
Figure FDA0002689281180000011
其中,Ii为特征图像中的某一像素点,且特征图像为步骤3处理后的人脸图像转换得到的灰度图像,qi为滤波器输出图像中相应的像素点,ak和bk为滤波窗口wk中的一组线性系数,k为滤波半径,滤波窗口wk的大小为5×5;
线性系数ak的计算公式为:
Figure FDA0002689281180000012
Figure FDA0002689281180000013
其中,ε为调整系数,取值为0.5,pi为步骤3处理后的人脸图像中的某一像素点,Gx(pi)和Gy(pi)分别为pi在水平和竖直方向上与临近像素点的亮度差值,Gradk(pi)为pi对应的梯度值;
线性系数bk的计算公式为:
Figure FDA0002689281180000021
其中,
Figure FDA0002689281180000022
为步骤3处理后的人脸图像在滤波窗口wk内像素点的均值,
Figure FDA0002689281180000023
为特征图像在滤波框wk内像素点的均值,
Figure FDA0002689281180000024
Figure FDA0002689281180000025
均采用包滤波器计算得到;
步骤5:将步骤4中所得的滤波后的人脸图像转换到RGB颜色空间;
步骤6:采用双线性插值法将步骤5中得到图像进行放大,放大倍率为N倍,输出最终人脸图像。
CN201710602140.XA 2017-07-21 2017-07-21 一种基于边缘保持滤波器的人脸图像滤波方法 Active CN107392869B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710602140.XA CN107392869B (zh) 2017-07-21 2017-07-21 一种基于边缘保持滤波器的人脸图像滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710602140.XA CN107392869B (zh) 2017-07-21 2017-07-21 一种基于边缘保持滤波器的人脸图像滤波方法

Publications (2)

Publication Number Publication Date
CN107392869A CN107392869A (zh) 2017-11-24
CN107392869B true CN107392869B (zh) 2020-12-01

Family

ID=60335847

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710602140.XA Active CN107392869B (zh) 2017-07-21 2017-07-21 一种基于边缘保持滤波器的人脸图像滤波方法

Country Status (1)

Country Link
CN (1) CN107392869B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108765516A (zh) * 2018-03-26 2018-11-06 阳光暖果(北京)科技发展有限公司 一种基于缩略图的五官区域染色方法
CN110991213A (zh) * 2019-05-14 2020-04-10 杨春燕 自动化大数据分配平台

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7191720B2 (en) * 2003-02-03 2007-03-20 Hbi Branded Apparel Enterprises, Llc Method of manufacture for stitchless garment
CN101316321B (zh) * 2007-05-30 2010-04-07 展讯通信(上海)有限公司 一种基于中值滤波器的图像噪声去除方法和装置
CN100508891C (zh) * 2007-09-04 2009-07-08 陈武凡 在pet成像中最大后验优化图像重建方法
CN101625755B (zh) * 2009-08-06 2012-09-05 西安电子科技大学 基于分水岭-量子进化聚类算法的图像分割方法
US8101069B2 (en) * 2009-11-10 2012-01-24 Sky Bleu Martin External filtering and absorbing device for use in a local containment area
CN101968850B (zh) * 2010-10-21 2012-12-12 重庆大学 一种模拟生物视觉机制的人脸特征提取方法
CN101980286B (zh) * 2010-11-12 2012-02-08 西安电子科技大学 结合双树复小波与双变量模型的sar图像降斑方法
CN102509346A (zh) * 2011-09-30 2012-06-20 北京航空航天大学 基于边缘保持的对象光照迁移方法
CN102881011B (zh) * 2012-08-31 2015-04-29 北京航空航天大学 基于区域分割的肖像光照迁移方法
CN103295010B (zh) * 2013-05-30 2016-06-29 西安理工大学 一种处理人脸图像的光照归一化方法
CN103489166A (zh) * 2013-10-12 2014-01-01 大连理工大学 一种基于双边滤波的单幅图像去雾方法
CN104331683B (zh) * 2014-10-17 2017-07-07 南京工程学院 一种具有噪声鲁棒性的人脸表情识别方法
CN105096273A (zh) * 2015-08-20 2015-11-25 北京邮电大学 一种彩色图像亮度的自动调整方法
CN105469407B (zh) * 2015-11-30 2018-06-22 华南理工大学 一种基于改进的引导滤波器的人脸图像图层分解方法
CN105654437B (zh) * 2015-12-24 2019-04-19 广东迅通科技股份有限公司 一种对低照度图像的增强方法
CN105809677B (zh) * 2016-03-03 2020-12-22 武汉Tcl集团工业研究院有限公司 一种基于双边滤波器的图像边缘检测方法及系统

Also Published As

Publication number Publication date
CN107392869A (zh) 2017-11-24

Similar Documents

Publication Publication Date Title
CN108921800B (zh) 基于形状自适应搜索窗口的非局部均值去噪方法
CN105654436B (zh) 一种基于前景背景分离的背光图像增强去噪方法
CN108205804B (zh) 图像处理方法、装置及电子设备
CN106780311B (zh) 一种结合皮肤粗糙度的快速人脸图像美化方法
US9495582B2 (en) Digital makeup
CN102222328B (zh) 一种边缘保持的自然场景图像自适应加权滤波方法
CN108288258B (zh) 一种针对恶劣天气条件下的低质图像增强方法
CN106683080B (zh) 一种视网膜眼底图像预处理方法
CN111968041A (zh) 一种自适应图像增强方法
CN111598791B (zh) 一种基于改进动态大气散射系数函数的图像去雾方法
CN109544583B (zh) 一种提取皮革图像感兴趣区域的方法、装置及设备
CN109636745B (zh) 基于分数阶微分图像增强算法的最佳阶图像增强方法
CN107292834B (zh) 红外图像细节增强方法
CN107146209A (zh) 一种基于梯度域的单幅图像去雾方法
CN107392869B (zh) 一种基于边缘保持滤波器的人脸图像滤波方法
CN117252773A (zh) 基于自适应颜色校正和导向滤波的图像增强方法及系统
Lei et al. Low-light image enhancement using the cell vibration model
CN111429362A (zh) 一种内窥镜荧光图像的血管增强方法
CN111598790A (zh) 基于对偶深度神经网络的高光谱遥感图像去噪方法
CN103886558A (zh) 基于LoG算子改进的自适应阈值小波去噪算法
CN109064413B (zh) 图像对比度增强方法及采用其的图像采集医疗设备
CN110706180B (zh) 一种极暗图像视觉质量提升方法、系统、设备及介质
Lal et al. Automatic method for contrast enhancement of natural color images
CN116109743B (zh) 一种基于ai和图像合成技术的数字人生成方法及系统
CN116630198A (zh) 一种结合自适应伽马校正的多尺度融合水下图像增强方法

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210205

Address after: Room 508, block a, Rongcheng cloud Valley, No.57, Keji Third Road, high tech Zone, Xi'an City, Shaanxi Province, 710075

Patentee after: Xi'an Huizhi Big Data Technology Co.,Ltd.

Address before: 710064 Shaanxi province Xi'an Beilin District South Second Ring Road No. 33 in the middle

Patentee before: CHANG'AN University