CN104637064B - 一种基于边缘强度权重的离焦模糊图像清晰度检测方法 - Google Patents
一种基于边缘强度权重的离焦模糊图像清晰度检测方法 Download PDFInfo
- Publication number
- CN104637064B CN104637064B CN201510091228.0A CN201510091228A CN104637064B CN 104637064 B CN104637064 B CN 104637064B CN 201510091228 A CN201510091228 A CN 201510091228A CN 104637064 B CN104637064 B CN 104637064B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- pixel
- gradient
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- 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/20024—Filtering details
-
- 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/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Studio Devices (AREA)
Abstract
本发明提供一种基于边缘强度权重的离焦模糊图像清晰度检测方法,首先对输入图像进行预处理,直方图均衡化方法修正亮度和对比度,对数码相机系统中存在的恒定功率加性噪声采用维纳滤波处理,使用中值滤波器和高斯滤波器分别对椒盐噪声和高斯噪声处理。然后采用四个方向边缘梯度算子对每个像素点的梯度进行检测,根据检测得到的梯度大小排除局部亮暗点和孤立噪声点的干扰,对剩余像素作进一步处理。接着对剩余像素的方向梯度与设定阈值比较,区分出边界上梯度值相对较大的强边缘像素和边界附近梯度相对较小的弱边缘像素并分别赋予不同的权重。最后所有像素进行最大梯度平方累加,得到整幅图像的清晰度检测值。
Description
技术领域
本发明涉及一种模糊图像清晰度检测方法,特别是一种基于边缘强度权重的离焦模糊图像清晰度检测方法,属于数字图像处理和光电跟踪测量领域。该方法实现基于数字图像处理技术的自动聚焦检测,对镜头和光电传感器采集到的每一帧数字图像给出清晰度评价并判断出聚焦最准确的图像,为后续控制镜头运动有效提供反馈和指导,在保持已有算法低复杂度的基础上,增加检测梯度方向,区分强弱边缘的作用,客观反映了图像的清晰程度,有效提高了算法的灵敏度、抗噪性和抗亮度变化等。可以广泛应用于各种光电测量设备和自动聚焦系统中。
背景技术
随着相机、摄像机等数码产品的流行,电子扫描仪、医学显微镜等精密仪器的发展,以及卫星导航、计算机视觉等高端技术的需要,对获取图像的质量要求越来越高。图像质量的好坏与清晰度直接相关,清晰度即模糊程度,是指人眼能感受出的相邻影像间的明显程度。在图像的采集、传输、压缩和滤波等各种处理过程中都可能产生模糊,比如曝光期间成像系统与被摄物体相对运动产生的运动模糊,成像系统聚焦不良或部分景物处于散焦状态产生的离焦模糊,光的衍射、以及压缩之后高频丢失等产生的各类模糊等。因此,为快速准确获取高质量清晰图像,需要对模糊图像进行清晰度评价并根据评价结果展开相应的反馈调节,即是自动聚焦的过程。
近年来,随着计算机硬件水平以及数字图像处理技术的发展,一种新兴的、基于数字图像处理的自动聚焦技术开始蓬勃发展起来。目前国内外对这一领域投入的大量的研究,很多成果已经在现有数字成像系统中广泛应用。基于图像技术的清晰度评价方法主要分为时域法、频域法、统计学法和信息熵法等四类,其中时域法主要包括Brenner函数、Tenengrad函数、Robert函数、Laplace函数、方差函数和平方梯度函数等灰度梯度法;频域法主要有功率谱法、小波变换法等;而统计学法有变化率函数、标准化变化率函数、Vollath自相关函数、基于梯度的标准相关函数等;信息熵法主要是利用图像的信息熵进行评价处理。上述方法中时域法算子简单,运算速度快,但精度不高,频域法灵敏度有所提高,但复杂度很高,实用受限,而基于统计理论和信息熵的清晰度评价函数对背景噪声、光照条件等比较敏感,这类评价函数在复杂场景下可能会失效。因此,目前各种基于数字图像处理的自动聚焦评价算法各有局限,存在各种缺点,需要进行进一步完善。
发明内容
本发明技术解决问题:针对现有技术的不足,提供一种基于边缘强度权重的离焦模糊图像清晰度检测方法,在保持较低复杂度的基础上,增加检测梯度方向,区分强弱边缘的作用,客观反映了图像的清晰程度,有效提高了算法的灵敏度、抗噪性和抗亮度变化等,适用于实时自动聚焦系统中。
为实现这样的目的,本发明的技术方案:一种基于边缘强度权重的离焦模糊图像清晰度检测方法,包括如下步骤:
步骤一、图像预处理:
基于数字图像处理技术的聚焦检测方法,主要是对镜头和CCD采集到的每一帧数字图像判断聚焦是否准确,成像系统是否清晰,并给出反馈信号控制镜头的运动,直到采集到的图像符合使用要求,即完成自动调焦。实际应用中,随着镜头位置和目标的不断变化,所获得的图像的平均亮度和视场会有所改变,成像器件暗电流背景噪声、图像采集过程中的椒盐噪声及电路处理过程中产生的固定模式噪声等都会对获取图像的清晰与否造成影响,因此,在进行模糊图像清晰度评价之前有必要进行预处理操作。图像预处理主要包括三个部分:直方图均衡化、维纳滤波、均值滤波和高斯滤波。首先对输入图像使用直方图均衡化方法修正亮度和对比度,修正输出图像的灰度动态范围为0~255(256个灰度级数),以达到符合人眼观察和后续处理的目的。然后对图像采用维纳滤波器处理,滤波器窗口选取邻域3×3,减少系统中存在的恒定功率加性噪声,尤其是固定模式噪声对图像质量的影响。经过维纳滤波后的图像会有轻微“振铃”效应,接着使用中值滤波器对可能存在的椒盐噪声进行有效滤除。最后,针对可能存在的高斯噪声,选取高斯低通滤波器处理,滤波器窗口大小为7×7,标准差为0.5。自此完成对输入图像的预处理操作。
步骤二、像素梯度判断:
对经过步骤一得到的图像进行像素梯度计算和判断处理。在进行像素梯度判断之前,选取一定的聚焦区域进行处理,好的聚焦区域可以减少参与运算处理的数据量,有利于提高聚焦速度。通常情况下,被摄物体大都位于场景的中央,我们选取中心区域作为聚焦窗口。中央窗口定义w为:
其中,m和n分别对应图像大小的行和列,f(x,y)是第x行,第y列的灰度值。
对中央窗口区域内的像素进行梯度计算和判断处理。梯度计算和判断的方向主要有水平方向、垂直方向和对角方向1以及对角方向2,对中心位置位于第x行,第y列的像素f(x,y),所述梯度检测算子如下:
水平方向梯度检测算子Ix:
垂直方向梯度检测算子Iy:
对角方向1梯度检测算子I1:
对角方向2梯度检测算子I2:
其中,水平方向是平面直角坐标系中0°或180°方向,垂直方向是平面直角坐标系中90°或-90°方向,对角方向1平面直角坐标系中45°或-135°方向,对角方向2是平面直角坐标系中135°或-45°方向。
根据公式(2)~(5)计算得到图像中每个像素点处的4个方向梯度,然后结合孤立噪声点和局部过亮过暗点的特性设定一个阈值T,利用如下公式(6)筛选出非孤立噪声和非局部过亮过暗点:
其中,筛选出非孤立噪声和非局部过亮过暗点的阈值T取值0.4。
步骤三、算子权重赋值:
对经过步骤二排除孤立噪声点和局部过亮过暗点后的像素进行算子权重赋值,筛选出的像素已经排除了噪声等因素对结果的影响,利用公式(7)、(8)对像素的方向梯度与设定阈值T1、T2进行比较,区分出边界上梯度值相对较大的强边缘像素和边界附近梯度相对较小的弱边缘像素,并对强边缘赋值大权重a,弱边缘赋值小权重值b。对中心位置位于第x行,第y列的像素f(x,y)的四个方向梯度分别为Ix、Iy、I1、I2,图像中的强弱边缘像素分别定义如下:
强边缘像素:max(Ix,Iy,I1,I2)-min(Ix,Iy,I1,I2)≥T2 (7)
弱边缘像素:T1≤max(Ix,Iy,I1,I2)-min(Ix,Iy,I1,I2)<T2 (8)
其中,判断强弱边缘的阈值T1=0.2,T2=1.4。
步骤四、计算全部像素点的梯度和值得到整幅图像的清晰度评价值:
经过步骤三处理,即可获得图像中非孤立噪声和非局部过亮过暗点的强弱边缘像素及其相应的权重值,利用如下像素梯度和算子计算位于第x行,第y列的像素f(x,y)的梯度值s(x,y):
计算得到每个像素点处的梯度和值s(x,y)后,最后累加聚焦区域内的全部像素点的梯度和值即可得到整幅图像的清晰度评价值,公式如下(10)所示:
本发明与现有技术相比的有益效果在于:
(1)本发明采用基于边缘强度权重的离焦模糊图像清晰度检测方法,与传统仅仅计算图像中某些方向的灰度梯度方法相比,增加到检测水平、垂直、对角方向的灰度梯度检测,选取了3×3模板算子,在定位精度和计算效率之间取得了很好的平衡。
(2)本发明采用基于边缘强度权重的离焦模糊图像清晰度检测方法,与常规自动聚焦方法相比,增加了直方图均衡化、维纳滤波、均值滤波和高斯滤波处理等图像预处理,在清晰度评价算子计算之前很好地减弱了噪声的影响,对图像的亮度和对比度作了修正,提高了抗噪声能力,拓宽了评价算法适用性。
(3)本发明采用基于边缘强度权重的离焦模糊图像清晰度检测方法,选取了中央窗口聚焦区域进行像素灰度梯度计算,与传统采用全窗口计算灰度梯度方法相比,减少了背景参与运算处理的数据量,有利于提高聚焦速度。
(4)本发明采用基于边缘强度权重的离焦模糊图像清晰度检测方法,不同于经典的聚焦算法把背景局部过亮过暗点和孤立噪声引入到聚焦算法中,使其灵敏度和精确度下降,本发明结合边缘和噪声特性,引入设定阈值,去除局部亮暗点和孤立噪声对算法灵敏度的影响,有效地抑制了噪声和局部点的干扰,提高了灵敏度和抗噪声能力。
(5)本发明采用基于边缘强度权重的离焦模糊图像清晰度检测方法,区别于现有大部分算法利用图像的全局性数据进行清晰度评价,本发明强调边界上梯度值相对较大的强边缘像素和边界附近梯度相对较小的弱边缘像素对人眼视觉清晰度贡献的差异性,并对强、弱边缘像素赋予不同的权重值,最终得到的清晰度评价结果更符合人眼视觉特性,也提高了聚焦图像的准确性和有效性。
总之:本发明在在保持较低复杂度的基础上,增加检测梯度方向,区分强、弱边缘的作用,客观反映了图像的清晰程度,有效提高了算法的灵敏度、抗噪性和抗亮度变化等,适用于实时自动聚焦系统中。
附图说明
图1为本发明方法的整体框架实现流程图;
图2为本发明选取的其中三幅不同离焦程度的标准输入图像和三幅不同离焦程度的真实序列图像,图2(a)、图2(b)和图2(c)分别是模糊圆半径为R=15、0、9时的不同离焦的标准彩色序列图像,图像分辨率为768×512像素,图2(d)、图2(e)和图2(f)分别是不同离焦的真实场景彩色序列图像,图像分辨率为768×576像素。为方便显示和记录,后续所有彩色图像均以黑白图像显示;
图3(a)、图3(b)分别是分辨率为768×512像素的标准过亮图像和处理后的正常图像,图3(c)、图3(d)分别是分辨率为768×576像素的真实噪声图像和处理后的正常图像;
图4(a)是上述方法对标准正常图像的清晰度评价值归一化曲线图;
图4(b)是上述方法对过暗图像的清晰度评价值归一化曲线图,其中过暗图像的灰度均值为23;
图4(c)是上述方法对过亮图像的清晰度评价值归一化曲线图,过亮图像的灰度均值为175;
图4(d)是上述方法对过弱噪声图像的清晰度评价值归一化曲线图,弱噪声图像的高斯噪声均值为0,方差为25;
图4(e)是上述方法对过强噪声图像的清晰度评价值归一化曲线图,强噪声图像的高斯噪声均值为0,方差为57;
图5(a)是上述方法对标准正常图像的清晰度评价值归一化曲线图;
图5(b)是上述方法对过暗图像的清晰度评价值归一化曲线图,其中过暗图像的灰度均值为32;
图5(c)是上述方法对过亮图像的清晰度评价值归一化曲线图,过亮图像的灰度均值为174;
图5(d)是上述方法对过弱噪声图像的清晰度评价值归一化曲线图,弱噪声图像的高斯噪声均值为0,方差为25;
图5(e)是上述方法对过强噪声图像的清晰度评价值归一化曲线图,强噪声图像的高斯噪声均值为0,方差为57。
具体实施方式
下面结合附图对本发明的实施例作详细说明。本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于以下的实施例。
如图1所示,本实施例的算法流程分为四个步骤:图像预处理、像素梯度判断、算子权重赋值、梯度求和。
本实例提供了一种基于边缘强度权重的离焦模糊图像清晰度评价方法,具体包括如下步骤:
步骤一:图像预处理。本实施例选取31幅分辨率为768×512像素,不同离焦程度的标准彩色图像序列和实际拍摄的60幅分辨率768×576像素,不同离焦程度的真实场景图像序列分别作为输入图像。不同离焦程度的标准彩色图像序列是通过成像系统的高斯离焦退化模型处理,即理想的点扩散函数用二维高斯函数替代并与图像卷积,二维高斯函数如公式(11)所示:
其中
σ为扩散函数,R为模糊圆半径,序列号从1~31分别是R=15...1,0,1,...15的不同离焦的图像。
实际应用中,随着镜头位置和目标的不断变化,所获得的图像的平均亮度和视场会有所改变,成像器件暗电流背景噪声、图像采集过程中的椒盐噪声及电路处理过程中产生的固定模式噪声等都会对获取图像的清晰与否造成影响,因此,在进行模糊图像清晰度评价之前有必要进行预处理操作。图像预处理主要包括三个部分:直方图均衡化、维纳滤波、均值滤波和高斯滤波。首先对输入图像使用直方图均衡化方法修正亮度和对比度,修正输出图像的灰度动态范围为0~255(256个灰度级数),以达到符合人眼观察和后续处理的目的。然后对图像采用维纳滤波器处理,滤波器窗口选取邻域3×3,减少系统中存在的恒定功率加性噪声,尤其是固定模式噪声对图像质量的影响。经过维纳滤波后的图像会有轻微“振铃”效应,接着使用中值滤波器对可能存在的椒盐噪声进行有效滤除。最后,针对可能存在的高斯噪声,选取高斯低通滤波器处理,滤波器窗口大小为7×7,标准差为0.5。自此完成对输入图像的预处理操作。
步骤二:像素梯度判断。对经过步骤一得到的图像进行像素梯度计算和判断处理。在进行像素梯度判断之前,选取一定的聚焦区域进行处理,好的聚焦区域可以减少参与运算处理的数据量,有利于提高聚焦速度。通常情况下,被摄物体大都位于场景的中央,按照公式(1)中央窗口定义选取中心区域作为聚焦窗口。对中央窗口区域内的像素进行梯度计算和判断处理。对中心位置位于第x行,第y列的像素f(x,y),如公式(2)~(5),计算和判断梯度的方向主要有平面直角坐标系中0°或180°方向的水平方向、90°或-90°方向的垂直方向和45°或-135°方向以及135°或-45°方向的两个对角方向共四个梯度方向。计算得到图像中每个像素点处的4个方向梯度后,然后结合孤立噪声点和局部过亮过暗点的特性设定一个阈值T=0.4,利用公式(6)筛选出非孤立噪声和非局部过亮过暗点,作进一步处理。
步骤三:算子权重赋值。对经过步骤二排除孤立噪声点和局部过亮过暗点后的像素进行算子权重赋值,筛选出的像素已经排除了噪声等因素对结果的影响,利用公式(7)~(8)对像素的方向梯度与设定阈值T1=0.2,T2=1.4进行比较,区分出梯度大的边界区域强边缘像素和模糊区域数量巨大的梯度相对较小的弱边缘像素,并对强边缘赋值大权重a=3,弱边缘赋值小权重值b=1。
步骤四:算子权重赋值。经过步骤三处理,获得图像中非孤立噪声和非局部过亮过暗点的强弱边缘像素点及其相应的权重值,利用公式(9)即可计算得到位于第x行,第y列的像素f(x,y)的梯度值s(x,y)。完成对每个像素点处的梯度和值s(x,y)计算后,最后利用公式(10)累加聚焦区域内的全部像素点的梯度和值即可得到整幅图像的清晰度评价值。依次对图像序列中的每一帧图像进行清晰度评价,可以得到全部图像序列的清晰度值,找出其中最清晰的图像序列号即可完成全部聚焦过程。
为了验证本发明方法的有效性,选取常用时域清晰度评价算子Brenner函数、平方梯度函数、Tenengrad梯度函数、Laplacian函数、Krish函数、Robert梯度函数、Variance方差函数以及本发明提供的方法分别对不同情况下的不同离焦程度的31幅标准序列图像和60幅真实场景图像进行聚焦实验,图4和图5分别是各种评价算子对不同程度离焦的标准图像序列和真实场景图像得出的清晰度评价值归一化曲线图。
图4为常用时域清晰度评价算子Brenner函数、平方梯度函数、Tenengrad梯度函数、Laplacian函数、Krish函数、Robert梯度函数、Variance方差函数以及本发明提供的方法对31幅分辨率为768×512像素不同程度离焦的标准图像序列实验结果图,横坐标代表图像序列号从1~31,纵坐标代表清晰度归一化评价值。
图5为常用时域清晰度评价算子和本发明提供的方法对60幅图像分辨率为768×576像素不同程度离焦的真实场景图像序列实验结果图,横坐标代表图像序列号从1~60,纵坐标代表清晰度归一化评价值。
从图4(a)和图5(a)可以看出,除Robert和Laplacian函数外,常用时域清晰度评价算子和提出的方法都能从不同模糊程度的正常序列中准确判断出标准图像第16幅和真实场景图像第26幅最清晰,满足调焦函数要求的单峰性、无偏性、单调性等特点;从图4(b)、(c)和图5(b)、(c)中可以看到,经过预处理后的提出算法对图像的亮度具有较好的适应能力,其它评价算子会随着亮度的变化灵敏度有所下降,部分算子甚至会失效;图4(d)、(e)和图5(d)、(e)可以看到,常用时域算子会随着噪声的干扰加大,灵敏度下降,当噪声大到一定程度,大部分评价算子失效,无法正确聚焦,而本发明提出的算法仍然具有较宽的陡峭宽度、较大的清晰度比率、较高的灵敏度和较小的局部极值因子,能够正确判断出最清晰图像的准确位置,这表明提出方法具有较好的抗噪性能。
为了评估检测方法的复杂度,在相同配置的计算机上使用不同评价算子对不同离焦程度的31幅标准序列图像进行实验,统计得到对全部图像处理的运行时间(单位s)如表1所示。
表1为在相同配置的计算机上使用不同评价算子对不同离焦程度的31幅标准序列图像进行实验,图像分辨率为768×512像素,统计得到对全部图像处理的运行时间(单位s)。
评价算子 | 运行时间(s) | 评价算子 | 运行时间(s) |
Brenner | 11.236451 | Laplacian | 14.469344 |
平方梯度 | 9.289950 | Robert | 9.870406 |
Tenengrad | 22.827378 | Variance | 89.195970 |
Krish | 16.944005 | Proposed | 23.911242 |
可以看出,本发明提出的清晰度检测算法的复杂度并不高,平均每幅图像处理时间仅0.7s,在所有常用检测算子中处于中等水平,但灵敏度和抗噪性能提升很大。因此,总的来说,本发明在在保持较低复杂度的基础上,增加检测梯度方向,区分强、弱边缘的作用,客观反映了图像的清晰程度,有效提高了算法的灵敏度、抗噪性和抗亮度变化等,为后续用于实时高精度自动聚焦系统提供了基础。
本发明未详细阐述部分属于本领域技术人员的公知技术。
本技术领域中的普通技术人员应当认识到,以上的实施例仅是用来说明本发明,而并非用作为对本发明的限定,只要在本发明的实质精神范围内,对以上所述实施例变化,变型都将落在本发明权利要求书的范围内。
Claims (5)
1.一种基于边缘强度权重的离焦模糊图像清晰度检测方法,其特征在于,包括如下步骤:
步骤一、图像预处理:
图像预处理主要包括:
(1.1)直方图均衡化,首先对输入图像使用直方图均衡化方法修正亮度和对比度,以达到符合人眼观察和后续处理的目的;
(1.2)维纳滤波和均值滤波:对图像采用维纳滤波器处理,减少系统中存在的恒定功率加性噪声,固定模式噪声对图像质量的影响;经过维纳滤波后的图像会有轻微“振铃”效应,接着使用中值滤波器对存在的椒盐噪声进行有效滤除;
(1.3)高斯滤波,针对可能存在的高斯噪声,选取高斯低通滤波器处理;自此完成对输入图像的预处理操作;
步骤二、像素梯度判断:
对经过步骤一得到的图像进行像素梯度计算和判断处理;在进行像素梯度判断之前,选取中心区域作为聚焦窗口,中央窗口定义w为:
<mrow>
<mi>w</mi>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mfrac>
<mn>3</mn>
<mn>8</mn>
</mfrac>
<mi>m</mi>
<mo>&le;</mo>
<mi>x</mi>
<mfrac>
<mn>5</mn>
<mn>8</mn>
</mfrac>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mfrac>
<mn>3</mn>
<mn>8</mn>
</mfrac>
<mi>n</mi>
<mo>&le;</mo>
<mi>x</mi>
<mfrac>
<mn>5</mn>
<mn>8</mn>
</mfrac>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munder>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,m和n分别对应图像大小的行和列,f(x,y)是第x行,第y列的灰度值;
对中央窗口区域内的像素进行梯度计算和判断处理;梯度计算和判断的方向有水平方向、垂直方向和对角方向1以及对角方向2,对中心位置位于第x行,第y列的像素f(x,y),梯度检测算子如下:
水平方向梯度检测算子Ix:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mi>x</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mo>*</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mo>*</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
垂直方向梯度检测算子Iy:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mi>y</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mo>*</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mo>*</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
对角方向1梯度检测算子I1:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mn>1</mn>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>*</mo>
<mi>f</mi>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mo>*</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
对角方向2梯度检测算子I2:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mn>2</mn>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>*</mo>
<mi>f</mi>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mo>*</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
根据公式(2)~(5)计算得到图像中每个像素点处的4个方向梯度,然后结合孤立噪声点和局部过亮过暗点的特性设定一个阈值T,利用如下公式(6)筛选出非孤立噪声和非局部过亮过暗点:
<mrow>
<mo>|</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>x</mi>
<mo>+</mo>
<mi>I</mi>
<mi>y</mi>
<mo>+</mo>
<mi>I</mi>
<mn>1</mn>
<mo>+</mo>
<mi>I</mi>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mn>4</mn>
</mfrac>
<mo>-</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>x</mi>
<mo>,</mo>
<mi>I</mi>
<mi>y</mi>
<mo>,</mo>
<mi>I</mi>
<mn>1</mn>
<mo>,</mo>
<mi>I</mi>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>></mo>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤三、算子权重赋值:
对经过步骤二排除孤立噪声点和局部过亮过暗点后的像素进行算子权重赋值,筛选出的像素已经排除了噪声因素对结果的影响,利用公式(7)对像素的方向梯度与设定阈值T1、T2进行比较,区分出边界上梯度值相对较大的强边缘像素和边界附近梯度相对较小的弱边缘像素,并对强边缘像素赋值大权重a,弱边缘像素赋值小权重值b;对中心位置位于第x行,第y列的像素f(x,y)的四个方向梯度分别为Ix、Iy、I1、I2,图像中的强弱边缘像素分别定义如下:
强边缘像素:max(Ix,Iy,I1,I2)-min(Ix,Iy,I1,I2)≥T2 (7)
弱边缘像素:T1≤max(Ix,Iy,I1,I2)-min(Ix,Iy,I1,I2)<T2 (8)
步骤四、计算全部像素点的梯度和值得到整幅图像的清晰度评价值
经过步骤三处理,即可获得图像中非孤立噪声和非局部过亮过暗点的强弱边缘像素及其相应的权重值,利用如下像素梯度和算子计算位于第x行,第y列的像素f(x,y)的梯度值s(x,y):
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>a</mi>
<mo>*</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>max</mi>
<mo>(</mo>
<mi>I</mi>
<mi>x</mi>
<mo>,</mo>
<mi>I</mi>
<mi>y</mi>
<mo>,</mo>
<mi>I</mi>
<mn>1</mn>
<mo>,</mo>
<mi>I</mi>
<mn>2</mn>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
<mi> </mi>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>x</mi>
<mo>,</mo>
<mi>I</mi>
<mi>y</mi>
<mo>,</mo>
<mi>I</mi>
<mn>1</mn>
<mo>,</mo>
<mi>I</mi>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>x</mi>
<mo>,</mo>
<mi>I</mi>
<mi>y</mi>
<mo>,</mo>
<mi>I</mi>
<mn>1</mn>
<mo>,</mo>
<mi>I</mi>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mi>T</mi>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>b</mi>
<mo>*</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>max</mi>
<mo>(</mo>
<mi>I</mi>
<mi>x</mi>
<mo>,</mo>
<mi>I</mi>
<mi>y</mi>
<mo>,</mo>
<mi>I</mi>
<mn>1</mn>
<mo>,</mo>
<mi>I</mi>
<mn>2</mn>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
<mi> </mi>
<mi>T</mi>
<mn>2</mn>
<mo>></mo>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>x</mi>
<mo>,</mo>
<mi>I</mi>
<mi>y</mi>
<mo>,</mo>
<mi>I</mi>
<mn>1</mn>
<mo>,</mo>
<mi>I</mi>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>x</mi>
<mo>,</mo>
<mi>I</mi>
<mi>y</mi>
<mo>,</mo>
<mi>I</mi>
<mn>1</mn>
<mo>,</mo>
<mi>I</mi>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mi>T</mi>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
<mi> </mi>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>x</mi>
<mo>,</mo>
<mi>I</mi>
<mi>y</mi>
<mo>,</mo>
<mi>I</mi>
<mn>1</mn>
<mo>,</mo>
<mi>I</mi>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>x</mi>
<mo>,</mo>
<mi>I</mi>
<mi>y</mi>
<mo>,</mo>
<mi>I</mi>
<mn>1</mn>
<mo>,</mo>
<mi>I</mi>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo><</mo>
<mi>T</mi>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
权重a=3,权重值b=1;
计算得到每个像素点处的梯度值s(x,y)后,最后累加聚焦区域内的全部像素点的梯度和值即得到整幅图像的清晰度评价值,公式如下(10)所示:
<mrow>
<mi>S</mi>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mfrac>
<mn>3</mn>
<mn>8</mn>
</mfrac>
<mi>m</mi>
<mo>&le;</mo>
<mi>x</mi>
<mfrac>
<mn>5</mn>
<mn>8</mn>
</mfrac>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mfrac>
<mn>3</mn>
<mn>8</mn>
</mfrac>
<mi>n</mi>
<mo>&le;</mo>
<mi>x</mi>
<mfrac>
<mn>5</mn>
<mn>8</mn>
</mfrac>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<mi>s</mi>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
2.根据权利要求1所述的基于边缘强度权重的离焦模糊图像清晰度检测方法,其特征在于,所述步骤一中直方图均衡化修正输出图像的灰度动态范围为0~255,256个灰度级数,维纳滤波器的处理窗口选取邻域3×3,高斯滤波器的窗口大小为7×7,标准差为0.5。
3.根据权利要求1所述的基于边缘强度权重的离焦模糊图像清晰度检测方法,其特征在于,所述步骤二中水平方向是平面直角坐标系中0°或180°方向,垂直方向是平面直角坐标系中90°或-90°方向,对角方向1平面直角坐标系中45°或-135°方向,对角方向2是平面直角坐标系中135°或-45°方向。
4.根据权利要求1所述的基于边缘强度权重的离焦模糊图像清晰度检测方法,其特征在于,所述步骤二中筛选出非孤立噪声和非局部过亮过暗点的阈值T设定为0.4。
5.根据权利要求1所述的基于边缘强度权重的离焦模糊图像清晰度检测方法,其特征在于,所述步骤三中判断强弱边缘的阈值T1=0.2,T2=1.4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510091228.0A CN104637064B (zh) | 2015-02-28 | 2015-02-28 | 一种基于边缘强度权重的离焦模糊图像清晰度检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510091228.0A CN104637064B (zh) | 2015-02-28 | 2015-02-28 | 一种基于边缘强度权重的离焦模糊图像清晰度检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104637064A CN104637064A (zh) | 2015-05-20 |
CN104637064B true CN104637064B (zh) | 2017-09-12 |
Family
ID=53215770
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510091228.0A Active CN104637064B (zh) | 2015-02-28 | 2015-02-28 | 一种基于边缘强度权重的离焦模糊图像清晰度检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104637064B (zh) |
Families Citing this family (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794702B (zh) * | 2014-10-08 | 2016-02-17 | 江苏省电力公司电力科学研究院 | 一种基于图像复杂度分析方法的特定场景下电力视频图像模糊检测方法 |
CN104917969B (zh) * | 2015-05-30 | 2018-01-19 | 广东欧珀移动通信有限公司 | 一种图像处理的方法及移动终端 |
CN105184747B (zh) * | 2015-09-09 | 2019-05-10 | 天津光电高斯通信工程技术股份有限公司 | 低照度图像对比度的提升方法 |
CN105631854A (zh) * | 2015-12-16 | 2016-06-01 | 天津天地伟业数码科技有限公司 | 一种基于fpga平台的自适应图像清晰度评价算法 |
US10325351B2 (en) * | 2016-03-11 | 2019-06-18 | Qualcomm Technologies, Inc. | Systems and methods for normalizing an image |
CN106097339A (zh) * | 2016-06-12 | 2016-11-09 | 成都甄识科技有限公司 | 一种基于模糊集的图像边缘检测方法 |
CN106295509B (zh) * | 2016-07-27 | 2019-11-08 | 浙江工业大学 | 一种面向非均匀退化视频中对象的结构化跟踪方法 |
CN107578373A (zh) * | 2017-05-27 | 2018-01-12 | 深圳先进技术研究院 | 全景图像拼接方法、终端设备及计算机可读存储介质 |
CN107240078A (zh) * | 2017-06-06 | 2017-10-10 | 广州优创电子有限公司 | 镜头清晰度验算方法、装置及电子设备 |
CN109131076A (zh) * | 2017-06-27 | 2019-01-04 | 上海蔚兰动力科技有限公司 | 驾驶辅助系统及驾驶辅助方法 |
CN107507173B (zh) * | 2017-08-15 | 2021-07-27 | 上海交通大学 | 一种全切片图像的无参考清晰度评估方法及系统 |
CN107688222A (zh) * | 2017-09-18 | 2018-02-13 | 济南中维世纪科技有限公司 | 一种基于虚焦报警算法的用于定焦ipc上的自动对焦方法 |
CN108269240B (zh) * | 2017-12-29 | 2021-03-23 | 杭州昕龙医疗科技有限公司 | 一种图像自动聚焦方法、电子设备及存储介质 |
CN108288267B (zh) * | 2018-01-17 | 2022-04-05 | 中国矿业大学 | 一种基于暗通道的扫描电镜图像清晰度无参考评价方法 |
CN109078342B (zh) * | 2018-08-17 | 2020-03-27 | 浙江大丰实业股份有限公司 | 舞台吊杆高度维护系统 |
CN110148147B (zh) * | 2018-11-07 | 2024-02-09 | 腾讯大地通途(北京)科技有限公司 | 图像检测方法、装置、存储介质和电子装置 |
KR102592605B1 (ko) * | 2018-12-06 | 2023-10-24 | 삼성전자주식회사 | 이미지 신호 처리기, 이미지 신호 처리기의 동작 방법, 및 이미지 신호 처리기를 포함하는 전자 장치 |
CN109785323B (zh) * | 2019-01-25 | 2024-01-30 | 淮阴师范学院 | 一种基于中频滤波的图像聚焦测度实现方法 |
CN109859151B (zh) * | 2019-01-25 | 2023-07-14 | 淮阴师范学院 | 一种基于局部直方图的聚焦测度实现方法 |
CN111860057A (zh) * | 2019-04-29 | 2020-10-30 | 北京眼神智能科技有限公司 | 人脸图像模糊和活体检测方法、装置、存储介质及设备 |
CN110378893B (zh) * | 2019-07-24 | 2021-11-16 | 北京市博汇科技股份有限公司 | 图像质量评价方法、装置和电子设备 |
CN110610470B (zh) * | 2019-09-18 | 2022-12-09 | 西安汇智信息科技有限公司 | 基于多方位梯度比较的相机多焦点清晰图像提取方法 |
CN110807745B (zh) * | 2019-10-25 | 2022-09-16 | 北京小米智能科技有限公司 | 图像处理方法及装置、电子设备 |
CN110852999B (zh) * | 2019-10-29 | 2023-03-10 | 北京临近空间飞行器系统工程研究所 | 图像扫描系统及图像扫描方法 |
CN110930363A (zh) * | 2019-10-29 | 2020-03-27 | 北京临近空间飞行器系统工程研究所 | 曲面模糊图像清晰度评价值确定方法、设备及存储介质 |
CN111161211B (zh) * | 2019-12-04 | 2023-11-03 | 成都华为技术有限公司 | 一种图像检测方法及装置 |
US11276156B2 (en) * | 2020-01-07 | 2022-03-15 | GM Global Technology Operations LLC | Gaussian image quality analysis tool and method for operation |
CN111242861B (zh) * | 2020-01-09 | 2023-09-12 | 浙江光珀智能科技有限公司 | 去除tof相机杂散光方法及装置、电子设备及存储介质 |
CN111340722B (zh) * | 2020-02-20 | 2023-05-26 | Oppo广东移动通信有限公司 | 图像处理方法、处理装置、终端设备及可读存储介质 |
CN113326722B (zh) * | 2020-02-29 | 2023-06-02 | 湖南超能机器人技术有限公司 | 基于序列模式的图像模糊检测方法及设备 |
CN113822835B (zh) * | 2020-06-04 | 2023-12-29 | 合肥君正科技有限公司 | 一种易于硬件实现的图像清晰度评价方法 |
CN112818737B (zh) * | 2020-12-18 | 2024-02-02 | 广州视源电子科技股份有限公司 | 视频识别方法、装置、存储介质及终端 |
CN112914727A (zh) * | 2021-03-19 | 2021-06-08 | 联仁健康医疗大数据科技股份有限公司 | 非目标障碍物分离方法、系统、医疗机器人及存储介质 |
CN112967265A (zh) * | 2021-03-22 | 2021-06-15 | 联仁健康医疗大数据科技股份有限公司 | 医用夹取方法、系统、医疗机器人及存储介质 |
CN113096014B (zh) * | 2021-03-31 | 2023-12-08 | 咪咕视讯科技有限公司 | 视频超分处理方法、电子设备及存储介质 |
CN113554628A (zh) * | 2021-07-27 | 2021-10-26 | 苏州微景医学科技有限公司 | 图像处理方法、装置及计算机可读存储介质 |
CN113808038A (zh) * | 2021-09-08 | 2021-12-17 | 瑞芯微电子股份有限公司 | 图像处理方法、介质及电子设备 |
CN114972141A (zh) * | 2022-05-13 | 2022-08-30 | 华侨大学 | 一种再模糊理论的双模态聚焦分析方法 |
CN115631171B (zh) * | 2022-10-28 | 2023-09-15 | 上海为旌科技有限公司 | 一种图片清晰度评价方法、系统及存储介质 |
CN115410135B (zh) * | 2022-11-01 | 2023-04-07 | 中国民航大学 | 自主式搬运的航空行李特征感知重构方法、系统及其应用 |
CN116152261B (zh) * | 2023-04-24 | 2023-06-27 | 济南奥盛包装科技有限公司 | 一种印刷制品质量的视觉检测系统 |
CN116527922B (zh) * | 2023-07-03 | 2023-10-27 | 浙江大华技术股份有限公司 | 图像编码方法及其相关装置 |
CN117152135B (zh) * | 2023-10-30 | 2024-01-23 | 济宁市市政园林养护中心 | 一种道路施工裂缝缺陷评估检测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101660946A (zh) * | 2009-09-28 | 2010-03-03 | 山东神戎电子股份有限公司 | 热像仪自动聚焦方法及装置 |
CN102404602A (zh) * | 2011-09-23 | 2012-04-04 | 浙江工业大学 | 一种基于清晰度测试卡的摄像机清晰度检测方法 |
CN102831392A (zh) * | 2012-07-09 | 2012-12-19 | 哈尔滨工业大学 | 一种远距离虹膜跟踪与采集装置及其方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8724009B2 (en) * | 2010-05-05 | 2014-05-13 | Hiok Nam Tay | Auto-focus image system |
US8366273B2 (en) * | 2011-01-31 | 2013-02-05 | National Chiao Tung University | Iris image definition estimation system using the astigmatism of the corneal reflection of a non-coaxial light source |
-
2015
- 2015-02-28 CN CN201510091228.0A patent/CN104637064B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101660946A (zh) * | 2009-09-28 | 2010-03-03 | 山东神戎电子股份有限公司 | 热像仪自动聚焦方法及装置 |
CN102404602A (zh) * | 2011-09-23 | 2012-04-04 | 浙江工业大学 | 一种基于清晰度测试卡的摄像机清晰度检测方法 |
CN102831392A (zh) * | 2012-07-09 | 2012-12-19 | 哈尔滨工业大学 | 一种远距离虹膜跟踪与采集装置及其方法 |
Non-Patent Citations (5)
Title |
---|
A Spectral and Spatial Measure of Local Perceived Sharpness in Natural Images;Guong T.Vu et al;《IEEE Transactions on Image Processing》;20110929;第21卷(第3期);第934-945页 * |
Efficient auto-focus algorithm for optical measurement system;Ni Jun et al;《Proceedings of the SPIE - The International Society for Optical Engineering》;20081119;第7283卷;第728344-1——728344-5页 * |
一种图像增强新方法;闫河 等;《计算机工程与应用》;20101210;第46卷(第34期);第195-198页 * |
基于区域加权的视频质量评价方法;朱宏 等;《宁波大学学报( 理工版)》;20130131;第26卷(第1期);第33-38页 * |
基于边缘差异的虚拟视图像质量评价方法;张艳 等;《电子与信息学报》;20130831;第35卷(第8期);第1894-1900页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104637064A (zh) | 2015-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104637064B (zh) | 一种基于边缘强度权重的离焦模糊图像清晰度检测方法 | |
Liu et al. | Estimating generalized gaussian blur kernels for out-of-focus image deblurring | |
CN110599433B (zh) | 一种基于动态场景的双曝光图像融合方法 | |
US7999867B2 (en) | Image edge detection apparatus and method, image sharpness emphasizing apparatus and method, recorded meduim recorded the program performing it | |
Li et al. | Single image dehazing using the change of detail prior | |
CN109714519B (zh) | 一种自动调整图像画面的方法和系统 | |
CN109410207A (zh) | 一种基于ncc特征的无人机巡线图像输电线路检测方法 | |
CN108596878A (zh) | 图像清晰度评价方法 | |
DE102010001520A1 (de) | Durch einen Flugzeugsensor unterstütztes Iriserfassungssystem und -Verfahren | |
CN104077759A (zh) | 一种基于色觉感知及全局质量因子的多曝光度图像融合方法 | |
CN106679581B (zh) | 一种高温环境下物体变形测量方法 | |
CN107726990B (zh) | 一种金属板料成形应变测量中点阵网格图像采集与识别方法 | |
CN103679672B (zh) | 基于边缘垂直距离匹配的全景图像拼接方法 | |
CN105139391B (zh) | 一种雾霾天气交通图像边缘检测方法 | |
CN103973990A (zh) | 宽动态融合方法及装置 | |
DE102014117120A1 (de) | Bildverarbeitungsvorrichtung | |
CN110189375A (zh) | 一种基于单目视觉测量的图像目标识别方法 | |
CN112019751B (zh) | 基于标定信息的自动聚焦方法 | |
CN108305232A (zh) | 一种单帧高动态范围图像生成方法 | |
TW202224404A (zh) | 相機模組的鏡頭髒污偵測方法 | |
CN106548471A (zh) | 粗精调焦的医学显微图像清晰度评价方法 | |
CN102542535A (zh) | 虹膜图像的去模糊方法 | |
Habeeb et al. | Contrast enhancement for visible-infrared image using image fusion and sharpen filters | |
Paul et al. | Filters used in x-ray chest images for initial stage tuberculosis detection | |
CN103778615B (zh) | 基于区域相似性的多聚焦图像融合方法 |
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 |