CN109785323A - 一种基于中频滤波的图像聚焦测度实现方法 - Google Patents
一种基于中频滤波的图像聚焦测度实现方法 Download PDFInfo
- Publication number
- CN109785323A CN109785323A CN201910101912.0A CN201910101912A CN109785323A CN 109785323 A CN109785323 A CN 109785323A CN 201910101912 A CN201910101912 A CN 201910101912A CN 109785323 A CN109785323 A CN 109785323A
- Authority
- CN
- China
- Prior art keywords
- image
- field emission
- value
- emission array
- focus field
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000001914 filtration Methods 0.000 title claims abstract description 28
- 230000009466 transformation Effects 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000003384 imaging method Methods 0.000 abstract description 8
- 239000000284 extract Substances 0.000 abstract description 2
- 230000001737 promoting effect Effects 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000006002 Pepper Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于中频滤波的图像聚焦测度实现方法,属于被动成像技术领域。该方法首先对图像进行傅里叶变换,得到频域系数,然后构造带通滤波器并对图像进行带通滤波。其次,对带通滤波后的图像进行分块处理并计算每一幅子图像像素灰度值的方差,并将其作为该子图像的清晰度信息。最后,求出所有子图像清晰度信息的均值,并把该均值作为整幅图像的聚焦测度值。本发明采用中频滤波和图像分块的方式来提取图像细节信息,具有原理简单、计算复杂度低的优点,同时通过上述步骤的实施,降低了噪声对图像细节信息的干扰,尤其在低对比度成像条件下的噪声鲁棒性强。本方法适合于相机的被动成像系统,便于推广使用。
Description
技术领域
本发明属于被动成像技术领域,具体涉及一种基于中频滤波的图像聚焦测度实现方法。
背景技术
日常生活中的拍照设备,例如单反相机、具有照相功能的手机以及路口的监控抓拍系统等都可以获得清晰的图像。然而清晰图像的获取是依赖于拍照设备的自动聚焦性能实现的。目前,市场上的拍照设备主要采用被动成像的自动聚焦技术。其核心是设计一个用于评价图像清晰度的聚焦测度,通过聚焦测度来选择最清晰的图像并保存。因此,一个性能优良的聚焦测度实现方法直接影响到所拍摄图像的质量。
目前使用较为广泛的图像聚焦测度方法是基于图像细节信息来构建的,例如基于图像边缘提取的聚焦测度。典型的方法有图像一阶高斯导数方法、二阶导数方法、一阶偏导数方法、梯度求和方法以及拉普拉斯求和方法等。此类方法的本质是构造一个尺寸为3×3或者5×5卷积模板,用该模板与整幅图像进行卷积运算。卷积运算的结果是提取出图像的边缘信息,然后采用取绝对值或者平方求和的形式来构造整幅图像的聚焦测度。此类构造聚焦测度的方法存在如下两个主要缺点。首先是卷积运算复杂度高,需要对整幅图像的所有像素点进行遍历操作,目前还没有较为成熟的快速算法,尤其是硬件设备上的快速算法,使得此类聚焦测度方法的聚焦实时性指标较差。其次是图像的噪声与边缘信息均属于高频信息,经过卷积运算后会增强噪声信息,也就是说此类聚焦测度易受噪声影响,最终导致错误的聚焦。
另外一类聚焦测度方法是基于图像变换的方法,即在变换域(频域)提取图像的高频信息来构建聚焦测度。典型的方法有细尺度小波系数的和、高频与低频小波系数的比率、基于离散余弦变换的聚焦测度和基于曲波变换的聚焦测度。这些基于变换的方法有一个共同的特点,就是对图像进行变换后提取高频信息,以此来作为图像聚焦测度值。这类方法和前面基于边缘提取的方法在思想方法上是一致的,均强调高频信息。只不过前者是在空域采用卷积的方法来构造聚焦测度,而后者是在频域采用变换的方式来构造聚焦测度。基于图像变换的方法也容易受噪声影响,并且有些变换的计算复杂度比较大,例如小波变换和曲波变换,均没有成熟的硬件快速算法。
上述两类方法还有一个共同的特点是以整幅图像的全局信息来构建聚焦测度。如果图像的背景比较均匀或者比较平滑,此时极易受噪声影响,使得相应的聚焦测度无法反映出图像的清晰度信息。例如用手机的拍照功能在光照条件比较弱的室内或者夜间场景下进行拍照,我们会发现(手机)相机的自动对焦功能并不好用,所拍摄出的图像存在模糊现象并且有颗粒感。这就是聚焦测度算法失效的一个体现。因此,如何构建具有噪声鲁棒性的聚焦测度具有重要的研究意义与实用价值。
发明内容
针对现有技术中存在的问题,本发明提供了一种具有强噪声鲁棒性的图像聚焦测度实现方法,本发明所采取的技术方案如下:
一种基于中频滤波的图像聚焦测度实现方法,包括如下步骤:
步骤S1:将原始图像像素的行数和列数均调整为2n的整数倍,其中n为正整数,得到图像f(x,y),图像f(x,y)的行数和列数分别用M和N表示,即x=0,1,...,M-1,y=0,1,...,N-1,然后计算图像的傅里叶变换,得到频域系数T(u,v),其中u=0,1,...,M-1,v=0,1,...,N-1;
步骤S2:构造频域的带通滤波器H(u,v),并与图像f(x,y)的频域系数T(u,v)进行乘法运算得到G(u,v),即G(u,v)=T(u,v)·H(u,v),以此来实现频域的带通滤波;
步骤S3:对G(u,v)进行逆傅里叶变换得到带通滤波后的图像g(x,y);
步骤S4:对g(x,y)进行分块处理,得到大小为2n×2n像素的子图像Si(x,y),其中i=1,2,...,M×N/22n;子图像Si(x,y)的变量取值为:x=0,1,...,2n-1,y=0,1,...,2n-1;
步骤S5:计算每一幅子图像Si(x,y)像素(灰度)值的方差σi,并作为该子图像的清晰度信息;
步骤S6:计算所有子图像清晰度信息σi的均值,并把该均值作为整幅图像的聚焦测度值。
优选地,所述步骤S1中的傅里叶变换公式为:
其中“i”为虚数单位。
优选地,所述步骤S2中的带通滤波器H(u,v)定义如下:
其中,T0为带通滤波的中心频率,取值为H(u,v)最大值的二分之一;W为带通滤波器的频带宽度,取值为(M+N)/4。
优选地,所述步骤S3中对G(u,v)进行逆傅里叶变换的公式为:
与现有技术相比,本发明具有的有益效果是:本发明原理简单,采用中频滤波和图像分块的方式来提取图像细节信息,具有计算复杂度低的优点,同时通过步骤S2~S6的实施,很大程度上降低了噪声对图像细节信息的干扰,使得本方法所得到的聚焦测度有较高的噪声鲁棒性,适合于相机的被动成像系统,尤其是低对比度成像条件下的抗噪声能力较强,适合于推广使用。
附图说明
图1为本发明的实施步骤框图。
具体实施方式
为了便于技术人员理解本发明的技术方案,现结合说明书附图和实施例对本发明的技术方案作进一步详细描述。
本发明提出了一种基于中频滤波的图像聚焦测度实现方法,其实施步骤框图如图1所示,在本实施例中,子图像分块参数选择n=3,则该方法的具体步骤细化为:
步骤S1:将原始图像像素的行数和列数调整为8的整数倍,可以通过对图像进行裁剪或者插值来实现,这样便得到图像f(x,y)。这里,f(x,y)的行数和列数分别用M和N表示。之所以对图像行数和列数进行调整是因为本发明所提出的方法是基于图像分块实现的,当n=3时,需要把图像划分为大小为8×8像素的若干幅子图像。计算图像f(x,y)的傅里叶变换,得到频域系数T(u,v),其中u=0,1,...,M-1,v=0,1,...,N-1。傅里叶变换公式为:
为了降低图像噪声(主要是高斯噪声、椒盐噪声和乘性冲击噪声)的影响,接下来的步骤S2~S6起到了关键作用。
步骤S2:构造频域的带通滤波器H(u,v),并与图像f(x,y)的频域系数T(u,v)进行乘法运算得到G(u,v),即G(u,v)=T(u,v)·H(u,v),以此来实现频域的带通滤波。带通滤波器H(u,v)定义如下:
其中,T0为带通滤波的中心频率,取值为H(u,v)最大值的二分之一。W为带通滤波器的频带宽度,取值为(M+N)/4。由于图像噪声属于高频信息,通过对频域的带通滤波可以滤除了大部分高频噪声信息,提高了聚焦测度的噪声鲁棒性,也为后续步骤的量化处理奠定了基础。
步骤S3:对G(u,v)进行逆傅里叶变换得到带通滤波后的图像g(x,y)。逆傅里叶变换的公式为:
步骤S4:对步骤S3所得到的图像g(x,y)进行分块处理,,子图像分块参数选择n=3,得到若干幅大小为8×8的子图像Si(x,y),其中i=1,2,...,M×N/64。之所以对图像进行分块处理是出于降低计算复杂度的考虑。这与传统的空域滤波计算过程中采用逐像素点进行处理的方式相比,大大降低了计算复杂度。此外,基于分块的计算方法还有平滑滤波的效果,能够进一步降低噪声对图像清晰度量化结果的影响。
步骤S5:计算每一幅子图像Si(x,y)像素(灰度)值的方差σi,并作为该子图像的清晰度信息。对于一幅聚焦良好的图像,其所包含的细节信息较多,具体有边缘信息或者区域边界信息,并且这些细节信息均出现在图像的局部区域中。而这部分信息的提取对于聚焦测度的计算至关重要。这也是本发明实施步骤S5中求子图像Si(x,y)像素值方差σi的原因。实际上,图像越清晰,图像的亮度变化越明显,从图像像素值的角度来考虑就是清晰图像的像素值具有较大的离散度。在统计学上通常使用方差来衡量这一离散特性。因此,本发明在步骤S5中用子图像Si(x,y)像素值的方差σi来作为子图像的清晰度信息。
步骤S6:计算所有子图像清晰度信息σi的均值,并把该均值作为整幅图像的聚焦测度值。
在其他施例中,可以根据待处理图像的像素值以及对聚焦方法的实时性需求等因素选择n值,n值越大,分块数量越少,方法的实时性越好,但是聚焦效果越差;反之,n值越小,分块数量越多,聚焦效果越好,但是方法的实时性变差,所以选择一个折中的n值非常重要。
应当说明的是,上述实施例均可根据需要自由组合。以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (4)
1.一种基于中频滤波的图像聚焦测度实现方法,其特征在于:具体步骤如下:
步骤S1:将原始图像像素的行数和列数均调整为2n的整数倍,其中n为正整数,得到图像f(x,y),图像f(x,y)的行数和列数分别用M和N表示,即x=0,1,...,M-1,y=0,1,...,N-1,然后计算图像的傅里叶变换,得到频域系数T(u,v),其中u=0,1,...,M-1,v=0,1,...,N-1;
步骤S2:构造频域的带通滤波器H(u,v),并与图像f(x,y)的频域系数T(u,v)进行乘法运算得到G(u,v),即G(u,v)=T(u,v)·H(u,v),以此来实现频域的带通滤波;
步骤S3:对G(u,v)进行逆傅里叶变换得到带通滤波后的图像g(x,y);
步骤S4:对g(x,y)进行分块处理,得到大小为2n×2n像素的子图像Si(x,y),其中i=1,2,...,M×N/22n;子图像Si(x,y)的变量取值为:x=0,1,...,2n-1,y=0,1,...,2n-1;
步骤S5:计算每一幅子图像Si(x,y)像素值的方差σi,并作为该子图像的清晰度信息;
步骤S6:计算所有子图像清晰度信息σi的均值,并把该均值作为整幅图像的聚焦测度值。
2.如权利要求1所述的基于中频滤波的图像聚焦测度实现方法,其特征在于:所述步骤S1中的傅里叶变换公式为:
其中“i”为虚数单位。
3.如权利要求1所述的基于中频滤波的图像聚焦测度实现方法,其特征在于:所述步骤S2中的带通滤波器H(u,v)定义如下:
其中,T0为带通滤波的中心频率,取值为H(u,v)最大值的二分之一;W为带通滤波器的频带宽度,取值为(M+N)/4。
4.如权利要求1所述的基于中频滤波的图像聚焦测度实现方法,其特征在于:所述步骤S3中对G(u,v)进行逆傅里叶变换的公式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910101912.0A CN109785323B (zh) | 2019-01-25 | 2019-01-25 | 一种基于中频滤波的图像聚焦测度实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910101912.0A CN109785323B (zh) | 2019-01-25 | 2019-01-25 | 一种基于中频滤波的图像聚焦测度实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109785323A true CN109785323A (zh) | 2019-05-21 |
CN109785323B CN109785323B (zh) | 2024-01-30 |
Family
ID=66504209
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910101912.0A Active CN109785323B (zh) | 2019-01-25 | 2019-01-25 | 一种基于中频滤波的图像聚焦测度实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109785323B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113960778A (zh) * | 2021-09-29 | 2022-01-21 | 成都西图科技有限公司 | 一种基于中频滤波的动态步长聚焦方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6087577A (ja) * | 1983-10-19 | 1985-05-17 | Matsushita Electric Ind Co Ltd | ビデオカメラの自動焦点調節装置 |
US6201899B1 (en) * | 1998-10-09 | 2001-03-13 | Sarnoff Corporation | Method and apparatus for extended depth of field imaging |
US6370279B1 (en) * | 1997-04-10 | 2002-04-09 | Samsung Electronics Co., Ltd. | Block-based image processing method and apparatus therefor |
US20070092242A1 (en) * | 2005-10-26 | 2007-04-26 | Canon Kabushiki Kaisha | Optical apparatus and focus adjustment method |
KR20070074293A (ko) * | 2006-01-09 | 2007-07-12 | 주식회사 팬택앤큐리텔 | 이동 통신 단말기 및 그 카메라 자동 초점 제어방법 |
CN101943839A (zh) * | 2010-07-06 | 2011-01-12 | 浙江大学 | 一体化的自动聚焦摄像机装置及清晰度评价方法 |
CN102129694A (zh) * | 2010-01-18 | 2011-07-20 | 中国科学院研究生院 | 一种图像显著区域检测方法 |
JP2011177373A (ja) * | 2010-03-02 | 2011-09-15 | Shimadzu Corp | X線撮影装置 |
WO2013126000A1 (en) * | 2012-02-21 | 2013-08-29 | Flir Systems Ab | Image processing method with detail-enhancing filter with adaptive filter core |
CN104637064A (zh) * | 2015-02-28 | 2015-05-20 | 中国科学院光电技术研究所 | 一种基于边缘强度权重的离焦模糊图像清晰度检测方法 |
CN105354817A (zh) * | 2015-09-25 | 2016-02-24 | 济南中维世纪科技有限公司 | 一种噪声图像自动聚焦方法 |
WO2017089736A1 (fr) * | 2015-11-27 | 2017-06-01 | Kerquest | Procede d'authentification et/ou de controle d'integrite d'un sujet |
CN107240092A (zh) * | 2017-05-05 | 2017-10-10 | 浙江大华技术股份有限公司 | 一种图像模糊度检测方法及装置 |
-
2019
- 2019-01-25 CN CN201910101912.0A patent/CN109785323B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6087577A (ja) * | 1983-10-19 | 1985-05-17 | Matsushita Electric Ind Co Ltd | ビデオカメラの自動焦点調節装置 |
US6370279B1 (en) * | 1997-04-10 | 2002-04-09 | Samsung Electronics Co., Ltd. | Block-based image processing method and apparatus therefor |
US6201899B1 (en) * | 1998-10-09 | 2001-03-13 | Sarnoff Corporation | Method and apparatus for extended depth of field imaging |
US20070092242A1 (en) * | 2005-10-26 | 2007-04-26 | Canon Kabushiki Kaisha | Optical apparatus and focus adjustment method |
KR20070074293A (ko) * | 2006-01-09 | 2007-07-12 | 주식회사 팬택앤큐리텔 | 이동 통신 단말기 및 그 카메라 자동 초점 제어방법 |
CN102129694A (zh) * | 2010-01-18 | 2011-07-20 | 中国科学院研究生院 | 一种图像显著区域检测方法 |
JP2011177373A (ja) * | 2010-03-02 | 2011-09-15 | Shimadzu Corp | X線撮影装置 |
CN101943839A (zh) * | 2010-07-06 | 2011-01-12 | 浙江大学 | 一体化的自动聚焦摄像机装置及清晰度评价方法 |
WO2013126000A1 (en) * | 2012-02-21 | 2013-08-29 | Flir Systems Ab | Image processing method with detail-enhancing filter with adaptive filter core |
CN104637064A (zh) * | 2015-02-28 | 2015-05-20 | 中国科学院光电技术研究所 | 一种基于边缘强度权重的离焦模糊图像清晰度检测方法 |
CN105354817A (zh) * | 2015-09-25 | 2016-02-24 | 济南中维世纪科技有限公司 | 一种噪声图像自动聚焦方法 |
WO2017089736A1 (fr) * | 2015-11-27 | 2017-06-01 | Kerquest | Procede d'authentification et/ou de controle d'integrite d'un sujet |
CN107240092A (zh) * | 2017-05-05 | 2017-10-10 | 浙江大华技术股份有限公司 | 一种图像模糊度检测方法及装置 |
Non-Patent Citations (7)
Title |
---|
LIQIANG GUO,LIAN LIU,HAIJIANG SUN: "Focus Measure Based on the Image Moments", INTERNATIONAL CONFERENCE ON MECHATRONICS AND AUTOMATION, pages 1151 - 1156 * |
SWATI DEWANGAN, ET.AL: "Image Smoothening and Sharpening using Frequency Domain Filtering Technique", RESEARCHGATE, vol. 5, no. 4, pages 120 * |
卢东兴,莫家庆: "边缘特征的光学图像清晰度判定与分析", 激光杂志, vol. 37, no. 7 * |
王勋,金剑秋,章志勇: "《图形图像数字水印方法》", 30 November 2011, pages: 35 * |
王岑,黄顺吉: "快速多项式变换(FPT)算法计算二维离散傅里叶变换(DFT)的一种新的改进方法", 《信号处理》, 2 April 1990 (1990-04-02) * |
邹建成,牛少彰: "数学及其在图像处理中的应用", 北京邮电大学出版社, pages: 120 * |
韩九强,杨磊: "《数字图像处理——基于XAVIS组态软件》", 31 August 2018, pages: 161 - 162 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113960778A (zh) * | 2021-09-29 | 2022-01-21 | 成都西图科技有限公司 | 一种基于中频滤波的动态步长聚焦方法 |
CN113960778B (zh) * | 2021-09-29 | 2024-07-30 | 成都西图科技有限公司 | 一种基于中频滤波的动态步长聚焦方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109785323B (zh) | 2024-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101901602B1 (ko) | 디지털 사진에서 노이즈를 제거하는 장치 및 방법 | |
Kim et al. | Wavelet-domain color image enhancement using filtered directional bases and frequency-adaptive shrinkage | |
KR101112139B1 (ko) | 부호화된 영상의 확대비 및 노이즈 강도 추정장치 및 방법 | |
WO2014114098A9 (zh) | 终端侧时间域视频质量评价方法及装置 | |
Akbari et al. | Adaptive saliency-based compressive sensing image reconstruction | |
CN109447903B (zh) | 一种半参考型超分辨率重构图像质量评价模型的建立方法 | |
CN110351453A (zh) | 一种计算机视频数据处理方法 | |
Zhang et al. | Single image dehazing based on fast wavelet transform with weighted image fusion | |
CN108614998B (zh) | 一种单像素红外目标检测方法 | |
CN109859196A (zh) | 一种基于分块pca的图像聚焦测度实现方法 | |
Rabie | Adaptive hybrid mean and median filtering of high-ISO long-exposure sensor noise for digital photography | |
CN109801289A (zh) | 一种基于Gabor变换的图像聚焦测度实现方法 | |
CN109801288A (zh) | 一种基于方向统计特性的图像聚焦测度实现方法 | |
CN106296591B (zh) | 基于马氏距离的非局部均匀数字图像去噪方法 | |
CN109785323A (zh) | 一种基于中频滤波的图像聚焦测度实现方法 | |
CN106375675A (zh) | 一种航空相机多曝光图像融合方法 | |
CN109886943A (zh) | 一种基于数论变换的图像聚焦测度实现方法 | |
CN109859194A (zh) | 一种基于局部边缘检测的图像聚焦测度实现方法 | |
Joy et al. | RECENT DEVELOPMENTS IN IMAGE QUALITY ASSESSMENT ALGORITHMS: A REVIEW. | |
CN109934876A (zh) | 一种基于二阶矩函数的图像聚焦测度实现方法 | |
CN109886942A (zh) | 一种基于图像峭度的聚焦测度实现方法 | |
Wang et al. | Spatio-temporal ssim index for video quality assessment | |
Kerouh et al. | Wavelet-based blind blur reduction | |
CN109934801A (zh) | 一种基于分块哈达玛变换的图像聚焦测度实现方法 | |
CN109859195A (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 |