CN107392869B - 一种基于边缘保持滤波器的人脸图像滤波方法 - Google Patents
一种基于边缘保持滤波器的人脸图像滤波方法 Download PDFInfo
- 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
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 230000003321 amplification Effects 0.000 claims 1
- 238000003199 nucleic acid amplification method Methods 0.000 claims 1
- 230000002146 bilateral effect Effects 0.000 description 13
- 230000000694 effects Effects 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000001815 facial effect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
Images
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
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
-
- 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/20172—Image enhancement details
- G06T2207/20192—Edge 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处理后的人脸图像进行滤波处理,边缘保持滤波器公式表述为:
其中,Ii为特征图像中的某一像素点,且特征图像为步骤3处理后的人脸图像转换得到的灰度图像,qi为滤波器输出图像中相应的像素点,ak和bk为滤波窗口wk中的一组线性系数,k为滤波半径。
进一步地,提出新的线性系数ak的计算公式为:
其中,ε为调整系数,pi为步骤3处理后的人脸图像中的某一像素点,Gx(pi)和Gy(pi)分别为pi在水平和竖直方向上与临近像素点的亮度差值,Gradk(pi)为pi对应的梯度值。
进一步地,调整系数ε为0.5。
进一步地,线性系数bk的计算公式为:
进一步地,滤波窗口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、根据新构建的边缘保持滤波器,按照公式对人脸图像进行滤波处理。取待滤波人脸图像的灰度图作为特征图像,Ii为其中的某一像素点,qi为滤波后输出的像素点,滤波框wk大小为5。公式中滤波系数ak和bk计算方法为:
取ε为0.5,通过计算待滤波人脸图像中某一点像素点pi在水平和竖直方向上与临近像素点的差值Gx(pi)和Gy(pi)得到点pi梯度Gradk(pi),一并作为系数ak计算公式的输入。为待滤波人脸图像在滤波框wk内像素点的均值,为特征图像在滤波框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不同算法对不同尺寸的图像处理用时的对比表
从表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处理后的人脸图像进行滤波处理,边缘保持滤波器公式表述为:
其中,Ii为特征图像中的某一像素点,且特征图像为步骤3处理后的人脸图像转换得到的灰度图像,qi为滤波器输出图像中相应的像素点,ak和bk为滤波窗口wk中的一组线性系数,k为滤波半径,滤波窗口wk的大小为5×5;
线性系数ak的计算公式为:
其中,ε为调整系数,取值为0.5,pi为步骤3处理后的人脸图像中的某一像素点,Gx(pi)和Gy(pi)分别为pi在水平和竖直方向上与临近像素点的亮度差值,Gradk(pi)为pi对应的梯度值;
线性系数bk的计算公式为:
步骤5:将步骤4中所得的滤波后的人脸图像转换到RGB颜色空间;
步骤6:采用双线性插值法将步骤5中得到图像进行放大,放大倍率为N倍,输出最终人脸图像。
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)
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)
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集团工业研究院有限公司 | 一种基于双边滤波器的图像边缘检测方法及系统 |
-
2017
- 2017-07-21 CN CN201710602140.XA patent/CN107392869B/zh active Active
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 |