光学显微成像系统图像清晰度评价函数的对比李雪;江旻珊【摘要】图像清晰度评价函数是评价各类成像系统成像质量的一个关键函数,为找到合适的图像清晰度评价算法,采用MATLAB软件对16种适用于光学显微成像系统的清晰度评价函数进行仿真,定量分析了不同算法的灵敏度、单峰性、无偏性以及运算速度.实验表明:Laplacian函数具有较高的单峰性、无偏性和灵敏度;存在高斯噪声时,Brenner函数、Tenengrad函数和基于Prewitt算子的函数以及中值滤波-离散余弦函数稳定性好;而存在椒盐噪声时,Roberts函数综合性能最优.【期刊名称】《光学仪器》【年(卷),期】2018(040)001【总页数】11页(P28-38)【关键词】清晰度评价函数;高斯噪声;椒盐噪声;Laplacian函数【作者】李雪;江旻珊【作者单位】上海理工大学光电信息与计算机工程学院,上海200093;上海理工大学光电信息与计算机工程学院,上海200093【正文语种】中文【中图分类】TG115.21引言图像清晰度评价函数是评价光学显微成像系统成像质量的关键函数,而一个良好的图像清晰度评价函数应满足单峰性、无偏性、灵敏度高、信噪比大及计算量大的要求[1]。

目前国内外学者提出了很多图像清晰度评价函数,每种评价函数各有优缺点,并且在不同的环境下有不同的评价效果,但大多只能在某些特定情况下才具有较好的效果。

由于自动对焦的系统十分复杂和多变,一种函数无法适应所有的对焦场合,因此本文将对目前比较常见的十几种图像清晰度评价函数进行对比分析,研究其主要性能和特点,以便为这些函数的正确使用提供一定的依据。

1 图像清晰度评价函数(1) Brenner梯度函数(Brenner)[2]Brenner梯度函数仅仅考察被判断点与其中一个相邻像素点之间的灰度差值,简洁实用、计算量少,其表达式为(1)式中I(x,y)为像素点(x,y)处的灰度值。

Brenner梯度算子可以看作是模板T=[-1 0 1]和对应位置的图像像素[I(x,y)I(x+1,y) I(x+2,y)]依次进行卷积,模板T=[-1 0 1]是一个带通滤波器,Brenner梯度算子正是通过带通滤波来滤除比例较大的低频能量,保留图像中的中频部分能量。

(2) 改进的Brenner梯度函数(ImprovedBre)改进算法采用两个模板滤波器来克服阈值对传统Brenner算法评价结果的影响[3]。

两个滤波器模板分别是带通滤波器T=[-1 0 1]和高通滤波器G=[1 -1]。

用两个滤波器模板分别对图像进行滤波,通过计算低频部分能量和高频部分能量来评价图像的清晰度。

改进的清晰度评价算法可定义为(2)改进的算法在滤除比例较大的低频成分的同时,保留了图像细节丰富的中高频成分。

(3) 绝对方差函数(AbsVar)[4]绝对方差函数是一个较为简单的算法,它仅仅是将像素点与其相邻位置的特定一个像素进行比较,其表达式为(3)因为正焦图像的边缘有较大的灰度差值,所以正焦图像经过这个算法计算之后会得到相对于离焦图像更大的值,且图像越模糊,其值越小。

(4) Roberts梯度函数(Roberts)[5]该函数对于绝对方差函数来说,考察的不仅仅是像素点和其中一个邻域像素的灰度关系,而是将像素点和其右下角的邻域像素进行比较,同时再比较其垂直方向下的两个邻域像素之间的灰度差值,两组差值的和作为函数的最终结果。

函数表达式为(4)该函数使用了被判断点及其外沿三个像素点的灰度信息,实际上是以某一点为中心的连续梯度的近似,处理边缘特性的时候比绝对方差和Brenner梯度函数要好。

(5) Laplacian函数(Laplacian)[6]Laplacian函数就是利用了边缘检测的Laplacian算子,其表达式为(5)式中:T为阈值;G(x,y)为图像与Laplacian算子的卷积。

Laplacian算子为(6)该评价函数使用了被判断点及其周围四个像素点的灰度信息。

Laplacian函数是二级微分算子,对孤立噪声点的响应是对阶跃边缘响应的4倍,对单像素线条的响应是对阶跃边缘响应的2倍,对线端及斜向边缘的响应大于垂直及水平走向边缘的响应,所以它不及Roberts梯度函数和梯度向量模方函数。

(6) Tenengrad函数(Tenengrad)该函数和Laplacian函数类似,同样运用边缘检测的思想[7]。

与之不同的是Tenengrad函数采用的是Sobel算子,利用Sobel算子来估算图像在水平方向和垂直方向的梯度,为使图像边缘的梯度放大,对梯度进行平方运算。

其定义为(7)其中(8)式中Gx(x,y)和Gy(x,y)分别是图像与gx和gy的卷积,其中gx和gy分别是水平方向和垂直方向的Sobel算子,可表示为(9)Tenengrad函数是先加权平均再微分,分别对水平和垂直方向进行模板运算,对噪声有一定的抑制能力。

(7) 梯度向量模方函数(SGVM)梯度向量模方函数是一个灰度变化梯度和的表达式,只选取了梯度标量数值信息作为灰度变化量描述[8],其函数定义为(10)(8) 自相关函数(Autocorrection)当图像细节越多,即高频能量越丰富时,自相关函数的曲线就越尖锐[9]。

该方法兼顾空间域函数的优点并且克服了空间域函数对噪声干扰敏感和对对焦区域选择要求高的缺点,具有较好的实时响应性能。

该函数表达式为(11)该函数能较好地克服图像中噪声和高亮目标等因素的影响,但是计算量较大。

(9) 熵函数法(Entropy)熵函数法是根据香浓信息论提出来的,香浓认为熵越大的时候信息量就越多,然后将此原理应用到图像清晰度评价中,因此可以认为在图像能量一定的情况下,图像熵越大图像越清晰[10]。

定义为(12)(10) 全频段积分函数(Integral)[11]在整个频域段进行积分,表达式为(13)(11) 相邻灰度差分算子绝对值之和(Sum)[12]用差分绝对值代替乘方和开方,即对被测点及其邻近点的灰度作差分运算,提取该点灰度值的变化大小,图像灰度差分绝对值之和算子为(14)式中M、N分别为图像长度和宽度的像素数。

(12) 中值滤波-离散余弦函数(MDCT)[13]函数表达式为(x,y)⊗g)2(15)其中该函数可以将图像信息由空域转换至频域,图像越模糊高频分量越少,模糊图像中相邻像素间的相互作用,表现为频域分布中高频成分的缺失,以高频成分的大小作为图像清晰度的评价依据。

(13) 图像能量函数(ImageEnergy)[14]使用图像能量函数作为清晰度评价函数的原因,是由于图像在离焦成像时,离焦量越大,弥散损失的能量就越大。

该函数表达式为(16)(14) 平面微分平方和(PlanarDiff)对图像进行微分运算提取图像中景物的边缘和轮廓,可以评判图像中高频分量的大小,并判断对焦正确与否[15]。

其表达式为(17)该算法最大的特点在于计算方法简单而且计算量小,适用于需要实时自动对焦的场合。

(15) 基于Prewitt边缘检测算子的清晰度评价函数(Prewitt)[16]Prewitt算子与Sobel算子非常类似,只是模板系数不同。

Prewitt算子的sx和sy 分别用卷积模板表示为和图像清晰度评价函数为(18)式中:和Sy(x,y)分别是图像与sx和sy的卷积。

Prewitt边缘算子是一种一阶微分算子,是在图像空间利用两个方向模板与图像进行邻域卷积来完成检测。

(16) 基于LOG边缘检测算子的清晰度评价函数[17](LOG)LOG算子对图像进行边缘检测时,输出的LOG(x,y)是通过卷积运算得到的。

LOG(x,y)常用模板为应用LOG边缘检测算子评价清晰度时,最后整幅图像的清晰度评价公式为(19)这种方法是先将图像和高斯滤波器进行卷积运算,这种运算既平滑了图像又降低了噪声,由于孤立的噪声点和较小的组织结构将同时被滤除,会不利于清晰度评价的精度,但此方法会增强调焦函数的单调性。

2 算法比较本文在多个环境下对16种算法进行了比较和验证,所有算法均由MATLAB编程,并且在相同的计算机环境下测试。

为了验证这些算法的性能,我们对同一个电路板拍摄了一组在纯白背景下的图像。

我们选择了在清晰成像位置附近,每隔0.05 mm采集一幅图像,焦面前和焦面后各10幅,一共采集到21幅图,将所得到的21幅图像从1到21依次排序,得到一组序列号,并在不同环境条件下对各个算法的单峰性,无偏性、灵敏度以及计算量进行比较。

2.1 单峰性、无偏性和灵敏度图1为各个函数运算并平滑归一化之后的结果。

图1 理想情况下各清晰度评价函数运算结果Fig.1 The results of all sharpness functions由图1可以看出,熵函数、全频段积分函数以及图像能量函数表现极差,甚至不能找出最清晰的那幅图像。

而Tenengrad函数和自相关函数虽然能分辨出最清晰的图像,但是仍然能清楚地看出存在局部极值点。

下面我们来对比一下其他算法的性能,理想情况下,运算结果如图2所示。

图2 理想情况下部分函数的运算结果Fig.2 The results of part of sharpness functions由图2可见:在单峰性方面,每个函数都体现了极好的单峰性,并且函数在第一幅图像到最后一幅图像的调焦范围内都只有一个极值点;在无偏性方面,每个函数的极值点都在横轴为10的位置,即每个函数的最佳焦距位置都指向第10幅图像的位置,函数之间的结果得到相互的证明;在灵敏度方面,Laplacian函数和中值滤波-离散余弦函数的灵敏度最高,基于Prewitt边缘检测算子的函数的灵敏度是这几个函数中最低的。

2.2 计算量图3为各个函数运行一次所耗费的时间。

图3 理想情况下各函数运行耗时Fig.3 Time- consumption of each algorithm 从图3可以看出,熵函数耗时最短,而相邻灰度差分算子绝对值之和函数耗时最长。

但是在耗时方面,由于受计算机因素的限制以及函数运行环境的影响,在此仅仅作为参考。

为了验证以上结果的可靠性,我们用有自动对焦功能的显微镜拍摄了南瓜茎细胞在显微镜离焦-合焦-离焦状态下的一段视频,从中间截取了25张图像对以上结果进行验证。

图4的3张图分别是显微镜离焦、合焦、再离焦时拍摄到的南瓜茎细胞图。

图4 南瓜茎细胞图像Fig.4 Images of cell of pumpkin stem理想情况下,各函数运算结果如图5所示。

图5 理想情况下各清晰度评价函数运算结果Fig.5 The results of all sharpness functions从图5 中可以看出:熵函数、图像能量函数以及全频段积分函数均不能识别出最清晰的那幅图像;其次,仍旧可以看到Tenengrad函数和自相关函数拥有几个不可忽视的局部峰值;其余的函数均表现出来较好的无偏性,均在第13幅图的位置达到了最大值,而且都只有一个极值点,证明了这些函数优秀的单峰性;而在灵敏度方面,Lapacian函数的灵敏度仍旧是最高的。