改进的Sobel算子图像清晰度评价算法.pdf

长春理工大学学报（自然科学版）Journal of Changchun University of Science and Technology（Natural Science Edition）Vol.46No.2Apr.2023第46卷第2期2023年4月收稿日期：2022-10-14基金项目：吉林省科技厅重点研发项目（20200401070GX）作者简介：张永超（1996-），男，硕士研究生，E-mail：通讯作者：李英（1979-），女，博士，副教授，E-mail：张永超，等：改进的Sobel算子图像清晰度评价算法改进的 Sobel 算子图像清晰度评价算法张永超，李英，唐智勇，谢康旗（长春理工大学光电工程学院，长春130022）摘要：针对传统图像清晰度评价算法在显微自动对焦系统中抗噪性弱的问题，在传统梯度函数基础上提出了一种改进的 Sobel 算子图像清晰度评价算法。首先利用自适应中值去噪算法对显微图像预处理，减少噪声干扰；其次增加Sobel 算子模板方向，采用四方向模板来计算像素点梯度；最终使用自适应全局阈值确定图像目标边缘，使用强边缘像素点的梯度平方和作为清晰度评价值，在一定程度上有效地解决了传统评价算子方位局限、抗干扰性差等问题。实验结果表明，该算法与传统梯度评价函数相比，图像处理速度提升了 15.6%，灵敏度因子平均提升了 0.711 9，具有良好的单峰性、抗干扰性和灵敏度。关键词：Sobel 算子；图像清晰度评价；自适应中值去噪；自适应全局阈值中图分类号：TP391文献标志码：A文章编号：1672-9870（2023）02-0096-10Improved Sobel Operator Image Sharpness Evaluation AlgorithmZHANG Yongchao，LI Ying，TANG Zhiyong，XIE Kangqi（School of Opto-Electronic Engineering，Changchun University of Science and Technology，Changchun 130022）Abstract：An improved Sobel operator image sharpness evaluation algorithm is proposed based on the traditional gradientfunction to solve the problem that the traditional image sharpness evaluation algorithm has weak noise resistance in the mi-cro autofocus system.First，the adaptive median denoising algorithm is used to preprocess the micro-image to reduce thenoise interference.Secondly，the Sobel operator template direction is increased，and the four-direction template is used tocalculate the pixel gradient.Finally，an adaptive global threshold is used to determine the edge of the image target，and thesum of the gradient squares of the strong edge pixels is used as the clarity evaluation value，which effectively solves theproblems such as the limitation of traditional evaluation operator orientation and poor anti-jamming.The experimental re-sults show that the algorithm improves image processing speed by 15.6%and sensitivity factor by 0.711 9 on average com-pared with the traditional gradient evaluation function.It has good unimodality，anti-jamming and sensitivity.Key words：Sobel operator；image sharpness evaluation；adaptive median denoising；adaptive global threshold随着现代自动化控制技术、微电子以及数字图像处理等技术的快速发展，自动聚焦技术已经被广泛地应用于照相机、显微镜、内窥镜和摄像机等光电成像设备中。目前显微镜自动对焦方法主要分为主动式对焦与被动式对焦1。主动式对焦通过测量镜头与被测物之间的距离，然后移动切片到焦点处以达到对焦的目的；被动式对焦主要通过图像处理技术，通过采集一系列图像进行分析判断离焦程度来实现对焦。由于基于图像处理技术的被动式对焦方法具有低功耗、低成本、高度集成化和智能化等优点，被广泛运用于光电对焦设备中。以数字图像处理为基础的自动对焦方法中，图像清晰度评价函数是决定整个自动对焦系统准确度和速度的关键。图像的清晰度值是可以反映对焦系统的离焦程度的，图像清晰度值越大，图像越清晰，边缘细节多，说明处于近焦区域；图像清晰度值越小，图像越模糊，说明处于离焦区域；在空域上表现为相邻像素点的灰度值变化较大，在频域上表现为高频分量较多。空间域函数计算简单，稳定性好，频域函数复杂计算量大。近年来人们对空域类梯度评价算法进行了一系列研究，以典型的边缘检测算子为例，如 Sobel、Kirsch、Laplace 等为基础对图像清晰度函数进行改进，以提升评价算法性能，其中Sobel 算子耗时短，操作简单。2017 年褚翔等人2提出了一种基于方向 Sobel 算子边缘加权的调焦评价函数，提升了边缘检测准确度，但运算量大。2018 年许乐灵等人3提出了一种引导滤波自适应双阈值并改进 Kirsch 算子的优化边缘检测算法。2020 年王灿芳等人4提出了一种结合了 Variance 函数与 Brenner 函数优点的新的图像清晰度评价算法。2020 年郑欢欢等人5提出了基于均值计算和差分处理角度构造了 55 的梯度权值模板的改进的 Sobel 算子的边缘检测算法。2022 年 邹 雨 彤 等 人6结 合 人 类 视 觉 系 统（HVS），提出了采用灵活性阈值的多方向权重Sobel 算子的改进图像清晰度评价方法。2022 年李雪梅等人7提出了将加权核范数最小化（WNNM）图像去噪算法与 Sobel 边缘检测算法相结合的边缘检测算法，提升了算法的抗噪性。通过以上文献可知，对基于边缘检测算子的图像清晰度函数的改进方法主要包括改进梯度模板、给梯度模板加权、改进图像边缘提取的阈值、结合其他算法等。但大多数算法较为复杂，运算量大，对硬件要求较高，因此需要对传统梯度函数进行改进，使评价算法既简单又具有强的抗噪性。空域类梯度评价函数中 Sobel 算子虽然可以滤除部分噪声，但是仍然对噪声较为敏感，而且存在方位局限性问题，因此选择在 Sobel算子优点的基础上对传统梯度清晰度算法进行改进。可以增加 Sobel 算子梯度方向至 4 方向（考虑到计算效率问题，不宜增加过多模板）改善方位局限性问题，结合自适应中值滤波和改进的最大类间自适应阈值算法实现图像清晰度计算，用以提升评价函数整体抗噪性、准确性和灵敏度。1传统梯度图像清晰度评价函数介绍目前较为常用的图像清晰度评价函数主要有空域函数、频域函数、信息学函数和统计学函数四大类8。其中，频域函数计算量大，不适于快速测量计算；信息学函数计算结果波动较大；统计学函数的计算结果精度较低。相较而言空域函数，即梯度评价函数计算简单，其利用边缘灰度梯度变化来评价图像清晰度，可识别更多灰度梯度变化，使其调焦范围大、稳定性好，应用范围更加广泛。本节先介绍几种常用的传统灰度梯度函数，用来与改进后的清晰度评价函数的性能进行对比。常用的梯度函数有能量梯度函数（EOG）、Roberts 函数、Tenengrad 函数、Brenner 函数、方差Variance 函数、拉普拉斯 Laplace 函数等9。处理图像时可以把图像看作二维的离散矩阵，然后利用梯度函数获取灰度图像信息，用梯度值来评价图像的清晰度。本节中假设图像为I，其尺寸大小为M N，图像I中某像素点（x，y）的灰度值为f（x，y），图像清晰度评价值为F。（1）能量梯度函数（EOG）能量梯度函数是将图像水平和与垂直方向相邻的像素点的灰度值做差，再进行平方运算，然后将逐个平方差运算的和作为评价值。定义如下：F=xyf()x+1,y-f()x,y2+f()x,y+1-f()x,y2（1）张永超，等：改进的Sobel算子图像清晰度评价算法第2期97长春理工大学学报（自然科学版）2023年（2）Roberts 函数Roberts 函数主要是利用对角线方向像素点灰度值差的平方和作为每个像素点新的梯度值，所有像素梯度值的累加和即为清晰度值。定义如下：F=xyf()x+1,y+1-f()x,y2+f()x+1,y-f()x,y+12（2）（3）Sobel 算子函数Sobel 算子清晰度函数是将图像与 Sobel 算子水平、垂直方向模板卷积计算梯度值，再将梯度平方和作为清晰度值。Sobel 算子两方向模板如下：gx=-101-202-101gy=-1-2-1000121（3）式中，gx、gy分别代表 Sobel 算子水平方向、垂直方向的模板。基于 Sobel 算子的图像清晰度函数F定义，公式如下：F=xyG()x,y2（4）G()x,y T（5）式中，T为阈值；G(x,y)为像素点(x,y)处新的梯度值。定义为：G()x,y=G2x()x,y+G2y()x,y（6）式中，Gx(x,y)和Gy(x,y)分别代表像素点(x,y)与水平、垂直方向模板上卷积的梯度值。其定义式如下：Gx()x,y=f()x,y gxGy()x,y=f()x,y gy（7）（4）Brenner 函数Brenner 函数是将相邻两个像素点灰度之差的平方累加作为清晰度值，优点是计算量小、灵敏度高，缺点是对噪声比较敏感。其定义式如下：F=xyf()x+2,y-f()x,y2（8）Brenner梯度算子可以看作是模板T=-1 0 1与对应位置的像素依次卷积，模板T=-1 0 1是一个带通滤波器，Brenner 梯度算子就是通过带通滤波来滤除大比例的低频能量，保留中频能量。（5）方差函数（Variance）方差函数 Variance 可以用来表示图像灰度信息的离散水平。处在焦点位置的图像灰度变化明显，离散水平高，方差大，而其他位置灰度变化小，方差也相对较小10。因此可以用灰度变化的平均程度来评判图像是否正确对焦。图像清晰度评价值如下：F=xyf()x,y-22（9）式中，为图像的平均灰度。=1MNxyf()x,y（10）式中，M、N是图像的尺寸。（6）Laplace 函数Laplace 函 数 是 将 灰 度 图 像 各 个 像 素 点 与Laplace 算子模板进行卷积运算，得到梯度矩阵G()x,y，取G()x,y的平方和作为评价值，公式如下：F=xyG2()x,y（11）G()x,y=f()x,y L（12）Laplace 算子如下所示：L=0101-41010（13）以上就是常用的梯度图像清晰度评价函数，不同的梯度评价函数因为评价方法计算方式不同，而拥有不同的优点和特点。但都普遍存在一个问题，在外部噪声的干扰下，传统梯度评价函数容易出现局部最值，评价曲线波动较大，灵敏度大大降低，不利于对焦。因此在传统梯度函数基础上，提出一种计算简单、灵敏度高、抗噪性强的图像清晰度函数对提高自动对焦精度具有重要意义。982改进的Sobel算子图像清晰度算法传统的基于 Sobel 算子的清晰度函数只使用水平和垂直两个方向的模板，方位局限、边缘定位不够准确并对含有噪声的图像较为敏感，清晰度曲线灵敏度低，不利于后续自动对焦。因此可以采取增加方向模板的方式对算子进行改进，因考虑到计算效率的问题，只增加对角线模板至四方向模板（0、45、90、135），同时结合自适应中值滤波算法和自适应全局阈值增强图像目标边缘强度，提升算法抗噪性，使用强边缘像素点的梯度平方和作为清晰度值，在一定程度上改善了方位局限性问题，并且提升了算法抗噪性和灵敏度。改进的 Sobel 算子图像清晰度算法主要分为三个部分：图像预处理、Sobel 梯度值计算、阈值确定边缘进行清晰度值计算。流程图如图 1 所示。图 1改进算法流程图2.1自适应中值去噪图像传感器 CCD 在获取和传输图像过程中会引入各种噪声，这些噪声对图像清晰度值计算影响特别大，不利于对焦，所以在对图像进行清晰度值计算之前需要对图像进行去噪预处理，减少噪声对评价算法的影响。图像滤波算法在空域和频域上可分为空间域滤波和频域滤波，空间域滤波法大多为点运算与小区域运算，速度快、计算量小；频域滤波需要对图像进行频域时域变换以及逆变换，算法较为复杂，处理速度慢。空间域滤波主要通过分析噪声在图像空间上的分布特性来去除噪声，常用的空间域滤波方法有均值滤波、中值滤波、统计排序滤波和自适应滤波等；中值滤波器的性能优于均值滤波器，故本文采用自适应中值滤波器进行图像去噪。符号定义如表 1 所示。表 1自适应中值滤波符号定义符号xyZmaxZminZmedmaxZxy物理意义滤波器作用区域区域内最大灰度值区域内最小灰度值区域内灰度中值滤波器允许的最大尺寸坐标为（x，y）点的灰度值自适应中值滤波算法主要被分为两个过程：A层和B层，描述如下：A层：若区域内最小灰度值、最大灰度值和灰度中值满足ZminZmedZmax，则转到B层，否则增 大 窗 口 尺 寸；如 果 窗 口 尺 寸 max，则 重 复A层，否则输出Zmed。B层：如果满足关系式ZminZxyZmax，则输出Zxy，否则输出Zmed。自适应中值去噪算法考虑了被处理像素点Zxy灰度值的情况，根据Zxy的大小决定是否需要进行滤波处理，自适应中值滤波器只有在满足条件Zxy=Zmin或Zxy=Zmax时才会改变被处理点的灰度值11，因此可以减少滤波过程中出现图像失真的问题。该算法的处理过程如下：（1）在以被处理点Zxy为中心的一个自定义窗口对像素点进行排序，找到灰度中值Zmed；（2）若灰度中值Zmed满足Zmed=Zmin或Zmed=Zmax，则增大窗口，重复步骤（1）；（3）若灰度中值Zmed满足ZminZmedZmax，则需要查看被处理点Zxy的灰度值；（4）若关系式ZminZmedZmax成立，则保持被处理像素点Zxy灰度值不变；张永超，等：改进的Sobel算子图像清晰度评价算法第2期99长春理工大学学报（自然科学版）2023年（5）若Zxy=Zmin或Zxy=Zmax，说明Zxy是考察窗口内的灰度极值，这时需要将被处理像素点的值替换为灰度中值Zmed。2.2四方向 Sobel 算子梯度计算传统 Sobel 算子只有水平和垂直两方向的模板，算法简单，计算量小，但对除水平和垂直方向以外的其他方向边缘并不敏感，尤其是受噪声干扰纹理复杂的显微图像，检测效果不佳，抗噪性差。本文将传统 Sobel 算子的方向模板由两个增加到了四个，可检测水平、垂直和对角线多个方向的边缘，改善了传统 Sobel 算子方位局限性的问题，使边缘检测和梯度计算结果更加准确。四方向 Sobel 算子如下：h1=-101-202-101h2=-1-2-1000121h3=012-101-2-10h4=-2-10-101012（14）其中h1、h2分别代表水平和垂直方向模板；h3、h4代表两个对角线方向模板。实际上图像的梯度方向可以取 0360这个范围内任何值，但是考虑到计算效率和边缘宽度限制等问题，故只定义了四个方向模板。虽然继续增加模板的数量可以进一步提高检测的准确性，但同时也降低了运算速度。传统 Sobel 算子处理一帧图片需要916 ms，本文改进的 Sobel 算子处理一帧图片需要 773 ms，相比传统 Sobel 算子运算速度提升了15.6%。扩展了方向模板后，需要将图像与大小为33 的 4 方向模板进行邻域卷积计算，重新计算像素点梯度值。图像中每个像素点的卷积计算公式如下：G1(x,y)=h1 f(x,y)G2(x,y)=h2 f(x,y)G3(x,y)=h3 f(x,y)G4(x,y)=h4 f(x,y)（15）其中，G1(x,y)、G2(x,y)、G3(x,y)、G4(x,y)分别代表灰度图像中每个像素点的水平方向（0）、垂直方向（90）和对角线方向（45、135）的卷积梯度值。计算完每个像素点四方向模板的卷积值后，采用范数为 1 来计算像素边缘强度值，即将水平方向、垂直方向和对角线方向卷积的绝对值之和作为中心梯度值，避免了以范数 2 计算中的开方与平方运算，提升了运算速度。最终可以得到新的梯度图像F(x,y)，则新的梯度图像像素点的梯度计算公式如下：F(x,y)=|G1(x,y)|+|G2(x,y)|+|G3(x,y)|+|G4(x,y)|（16）2.3自适应阈值强边缘像素点提取获得新的梯度图像后，需要对图像强边缘像素点进行提取，将强边缘像素点梯度平方和作为图像清晰度值，这时需要选取合适的阈值对图像进行分割。目前主要有两种自适应阈值选取的方法：基于最大熵原理的自适应阈值选取和基于 Otsu 方法的自适应阈值选取12。最大类间方差法不能准确地将目标与背景分开，目标边缘模糊且背景噪声大13。故采用一种改进的最大类间方差法对图像进行自适应阈值提取强边缘像素点，并计算图像清晰度值。改进的最大类间方差自适应阈值选取主要思想为：依据图像灰度范围，把图像分割为C0和C1两个部分。本节中设图像的灰度取值范围为0，L-1，自适应分割阈值为T，则C0=0，1，2，T，C1=T+1，L-1，各个像素点的平均梯度值如下：Hi=(x,y)WiF(x,y)ni（17）式中，F(x,y)为灰度梯度图像像素点的梯度值；ni为灰度值为i的频数；Wi为每一个像素点的梯度值；Hi为 各 个 像 素 点 的 平 均 梯 度 值；Wi=(x,y)I(x,y)=i。像素均值和梯度方差的计算如下：1000=0THiT,1=0THiL-1-T（18）20=0T(Hi-0)2,21=T+1L-1(Hi-1)2（19）式中，T和L分别为分割出两个区域C0和C1的灰度范围；0、1为两类像素各自的均值；20、21为各部分的梯度方差。灰度值梯度均值方差最大的阈值T计算公式如下：T=Arg max 2021（20）由公式（20）可知，当两方差乘积为最大时，对应的阈值T即为最佳阈值。根据自适应阈值T提取图像强边缘像素点，将图像分割成目标和背景两个部分。将强边缘像素点即像素点梯度值大于等于阈值T的梯度平方和作为图像的清晰度评价值，清晰度值计算公式如下：M=F(x,y)2,F(x,y)T（21）式中，M为被处理图像的清晰度值，该结果可用于后续对焦处理。3实验结果与分析3.1实验设备为了验证改进的 Sobel 算子图像清晰度算法的性能，用 HAYAER300 倍变倍显微镜头在环形光源辅助下，借助机械调焦装置等步长调节显微镜与被测物体间距离，通过 CCD 和图像采集卡获取被测物体一系列离焦、聚焦图片序列；然后在 Win10 操作系统下通过 Matlab 软件分别对改进算法的边缘检测效果、抗噪性和灵敏度进行 了 仿 真 实 验 与 计 算。CCD 型 号 为 IMPERX_bobcat-B2520，分辨率为 2 4482 050，像元尺寸为 3.45 mm，感光芯片尺寸为 66.93 mm，环形光源为 LED 可见光。图 2 为实验图像采集设备图。3.2边缘检测效果实验分析为验证改进算法对目标定位的准确性，将改进算法与传统 Sobel 算子边缘检测效果进行实验对比。实验采用直径为 3 mm 螺母在显微镜下离焦、聚焦、离焦图像序列来进行测试。下面选用其中两幅离焦和聚焦图片，如图 3 所示，进行边缘检测效果对比，实验结果如图 4、图 5 所示。图 2图像采集设备图（a）离焦图像（b）聚焦图像图 3部分图像序列（a）离焦图像边缘检测（b）聚焦图像边缘检测图 4传统 Sobel 边缘检测（a）离焦图像边缘检测（b）聚焦图像边缘检测图 5改进算法边缘检测图 4 是传统 Sobel 算子对螺母的边缘检测效果图，从图中可以看出传统 Sobel 算子由于其方位局限性问题，对目标定位不够准确，丢失了很多目标像素点，边缘检测效果不好。图 5 是本文改进算法边缘检测效果图，可以看出目标轮廓清晰，目标边界得到增强，同传统 Sobel 算法相比张永超，等：改进的Sobel算子图像清晰度评价算法第2期101长春理工大学学报（自然科学版）2023年图像边缘定位更加准确，弥补了传统算法上方位局限性的不足。3.3改进算法抗噪性实验分析CCD 在获取图像的过程中，会引入各种类型噪声，比如由于发热或照明引入的高斯噪声，传输过程中引入的椒盐噪声等。为了对改进算法的抗噪性进行验证，对图像序列分别添加椒盐噪 声 和 高 斯 噪 声，并 将 改 进 算 法 分 别 与 传 统Sobel 算法和传统梯度评价函数进行对比。这里选用螺母图像序列和文字图像序列作为实验对象，并对它们添加不同种类噪声进行抗噪性实验分析，两组实验图像序列如图 6、图 7 所示，这里只展示部分图像序列。图 6螺母图像序列图 7文字图像序列对螺母序列和文字序列分别添加椒盐噪声和高斯噪声，然后对加噪后的全部图像序列进行清晰度评价，将改进的 Sobel 算法和传统 Sobel算法比较，并将清晰度值归一化处理，最终得到两种算法的清晰度评价曲线。实验结果如图 8、图 9 所示。（a）添加椒盐噪声后螺母的图像清晰度评价曲线（b）添加高斯噪声后螺母的图像清晰度评价曲线图 8添加不同噪声后螺母的图像清晰度评价曲线（a）添加椒盐噪声后文字的图像清晰度评价曲线（b）添加高斯噪声后文字的图像清晰度评价曲线图 9添加不同噪声后文字的图像清晰度评价曲线图 8、图 9 是分别添加椒盐噪声和高斯噪声后两组图像序列的评价曲线，从中可以看出传统 Sobel 算子清晰度评价函数对噪声敏感，在椒盐噪声影响下灵敏度低，在高斯噪声影响下还会出现局部最大值，稳定性较差；相比之下，改进的 Sobel 算子清晰度评价函数无论在椒盐噪声还是高斯噪声的影响下仍能保持良好的单峰性，不存在局部多个最值的情况，具有良好的抗噪性和灵敏度。接下来分别给两组图像序列添加噪声密度为0.1、0.25的椒盐噪声以及均值为0，方差为0.02、0.05 的高斯噪声，检验本文改进算法与常用的灰度梯度图像清晰度评价函数的抗噪性。归一化的图像清晰度评价曲线如图 10、图 11 所示。从图 10、图 11 中可以看出常用的图像清晰度函数极易受噪声影响，在椒盐噪声影响下容易出现灵敏度低、单峰性差的问题，即评价曲线上升下降缓慢；在高斯噪声的影响下，常用的清晰度评价函数容易出现局部最大值最小值的情况，不能准确捕获到图像清晰度最佳的图像所对应的位置。本文改进的基于 Sobel 算子的图像清晰度函数，不论是在椒盐噪声影响下还是高斯噪声影响下都能保持住良好的单峰性和无偏性，并且灵敏度较高，相比传统的灰度梯度清晰度评价函数具有更加突出的抗噪性。3.4改进算法灵敏度实验分析图像清晰度函数是用来指导调焦机构找到最佳成像位置的评价函数，高灵敏度是一个评价函数的最基本要求。灵敏度在视觉上表现为评价曲线的陡峭程度，越陡则表明灵敏度越高。为了对本文改进算法的灵敏度进行定量分析，对显微图像序列进行灵敏度评估，采用灵敏度102图 8、图 9 是分别添加椒盐噪声和高斯噪声后两组图像序列的评价曲线，从中可以看出传统 Sobel 算子清晰度评价函数对噪声敏感，在椒盐噪声影响下灵敏度低，在高斯噪声影响下还会出现局部最大值，稳定性较差；相比之下，改进的 Sobel 算子清晰度评价函数无论在椒盐噪声还是高斯噪声的影响下仍能保持良好的单峰性，不存在局部多个最值的情况，具有良好的抗噪性和灵敏度。接下来分别给两组图像序列添加噪声密度为0.1、0.25的椒盐噪声以及均值为0，方差为0.02、0.05 的高斯噪声，检验本文改进算法与常用的灰度梯度图像清晰度评价函数的抗噪性。归一化的图像清晰度评价曲线如图 10、图 11 所示。从图 10、图 11 中可以看出常用的图像清晰度函数极易受噪声影响，在椒盐噪声影响下容易出现灵敏度低、单峰性差的问题，即评价曲线上升下降缓慢；在高斯噪声的影响下，常用的清晰度评价函数容易出现局部最大值最小值的情况，不能准确捕获到图像清晰度最佳的图像所对应的位置。本文改进的基于 Sobel 算子的图像清晰度函数，不论是在椒盐噪声影响下还是高斯噪声影响下都能保持住良好的单峰性和无偏性，并且灵敏度较高，相比传统的灰度梯度清晰度评价函数具有更加突出的抗噪性。3.4改进算法灵敏度实验分析图像清晰度函数是用来指导调焦机构找到最佳成像位置的评价函数，高灵敏度是一个评价函数的最基本要求。灵敏度在视觉上表现为评价曲线的陡峭程度，越陡则表明灵敏度越高。为了对本文改进算法的灵敏度进行定量分析，对显微图像序列进行灵敏度评估，采用灵敏度（a）添加噪声密度为 0.1 的椒盐噪声螺母序列评价曲线（b）添加噪声密度为 0.25 的椒盐噪声螺母序列评价曲线（c）添加均值为0方差为0.02的高斯噪声螺母序列评价曲线（d）添加均值为0方差为0.05的高斯噪声螺母序列评价曲线图 10不同椒盐、高斯噪声影响下螺母序列图像清晰度评价曲线张永超，等：改进的Sobel算子图像清晰度评价算法第2期103长春理工大学学报（自然科学版）2023年因子来评估图像的灵敏度，本文d定为偏离正焦位置一帧图。灵敏度因子的计算公式如下：Msen=Mmax-Mmax+Mmax（22）式中，Msen表示评价函数的灵敏度因子；Mmax表示最 佳 成 像 位 置 即 正 焦 位 置 处 图 像 清 晰 度 值；Mmax+表示偏离正焦d位置处图像清晰度值，d表示横坐标偏移量。表 2 为不同清晰度评价函数的灵敏度因子。从表 2 中可以看出改进的基于 Sobel 算子的图像清晰度函数灵敏度因子比传统的灰度梯度函数灵敏度更高，评价曲线更陡即曲线单调的一侧斜率大。本文改进算法相比传统 Sobel 算子函数灵敏度因子提升了 0.578 5，相较于传统灰度梯度函数提升了 0.578 50.829 5，平均灵敏度因子提升了 0.711 9，并且本文改进算法在噪声干扰下还能保持良好的稳定性，左右两侧单调性好，不存在局部最值，评价曲线无波动性。4结论本文介绍了空间域类梯度评价函数，分析了各自优缺点，在传统的梯度函数基础上，针对传统梯度函数抗噪性差的问题，提出了一种改进的 Sobel 算子图像清晰度评价函数，结合了自适应中值滤波器和最大类间自适应阈值法，提升了算法抗噪性和灵敏度。通过实验证明，本文改进的算法在运行速度和灵敏度因子上分别提高了 15.6%和 0.711 9，相比传统的灰度梯度评价（a）添加噪声密度为 0.1 的椒盐噪声文字序列评价曲线（b）添加噪声密度为 0.25 的椒盐噪声文字序列评价曲线（c）添加均值为0方差为0.02的高斯噪声文字序列评价曲线（d）添加均值为0方差为0.02的高斯噪声文字序列评价曲线图 11不同椒盐、高斯噪声影响下文字序列图像清晰度评价曲线表 2不同评价函数灵敏度因子函数名称灵敏度因子EOG0.199 0brenner0.233 2variance0.054 1Laplace0.067 1传统Sobel0.305 1改进算法0.883 6104函数具有单峰性好、灵敏度高、抗噪性好的特点，有利于自动对焦系统的优化。参考文献1 李惠光，王帅，沙晓鹏，等.显微视觉系统中自动聚焦技术的研究 J.光电工程，2014，41（8）：1-9.2褚翔，祝连庆，娄小平，等.基于改进 Sobel 算子的动态自动调焦算法研究J.应用光学，2017，38（2）：237-242.3 许乐灵，胡石.一种引导滤波自适应双阈值优化边缘检测算法 J.南京理工大学学报，2018，42（2）：177-182.4王灿芳，崔良玉，阎兵.基于 Variance-Brenner 函数的显微图像清晰度评价算法研究 J.装备制造技术，2020（10）：78-82.5 郑欢欢，白鱼秀，张雅琼.一种基于 Sobel 算子的边缘检测算法 J.微型电脑应用，2020，36（10）：4-6.6 邹雨彤，谭毅，罗一涵.检焦系统中基于改进 Sobel算子的图像清晰度评价函数研究 J.电子设计工程，2022，30（7）：150-155.7 李雪梅，钟坚.加权核范数的边缘检测在最小化图像去噪中的应用 J.自动化与仪器仪表，2022（4）：16-20.8鲍歌堂，赵辉，陶卫.图像测量技术中几种自动调焦算法的对比分析J.上海交通大学学报，2005（1）：121-124+128.9李敏.显微视觉系统中自动聚焦算法的研究及微动平台的控制 D.秦皇岛：燕山大学，2016.10 梁隆恺，赵晶，何勇军.一种显微镜自动聚焦算法J.哈尔滨理工大学学报，2018，23（2）：46-52.11 刘晓婷.基于自适应中值滤波的指纹脉冲噪声去除方法 J.警察技术，2022（4）：62-63.12张晶晶，范学良，蔡小舒.计算 Canny 算子双阈值的新方法 J.中国粉体技术，2011，17（1）：67-71.13任莹莹.改进最大类别方差法的电力设备视觉识别技术研究 J.中国新技术新产品，2021（8）：20-22.张永超，等：改进的Sobel算子图像清晰度评价算法第2期105