机器视觉中均匀照明成像的LED环形阵列设计 下载: 1310次
1 引言
随着图像处理、工业自动化、机器学习等技术的迅速发展[1],机器视觉[2]在这些领域中发挥了重要作用。获取高质量的图像是机器视觉系统的基本目标,因为图像质量将影响系统的准确性和可靠性[3]。Tarabanis等[4]对计算机视觉中用于识别目标特征的传感器系统进行了具体的研究,认为成像的关键部分可以概括为环境信息(如对象、传感器、照明)及视觉系统的任务信息(如对象特征检测、基于模型的对象识别、场景重建)。
光源对机器视觉检测系统的性能起着重要作用[5],精确的光学结构设计可以提高捕获图像的质量,准确地分离目标和背景信息
近年来,关于LED光学设计性能方面的研究有很多报道。Zhu等[11]利用数学模拟算法在三环LED阵列上设计了一个球体的内表面,实现了均匀照明,此外,还提出一种能量守恒的反射自由曲面。Marchand[12]提出了一种通过控制摄像机和光源的位置来改进机器视觉任务的方法。Li等[13]利用一种基于反馈方法的自由曲面光学系统实现了均匀照度和高出光率。这些研究结果表明,不均匀光照会使图像失真,在特定对象的特征识别、机器人的路径检测、对象分析等应用中会出现误差。结合机器视觉成像质量的影响因子,可知,照度均匀性是机器视觉需要解决的问题,而LED环形阵列灯具的优化设计是解决该问题的重要步骤。
本文在分析三维空间成像过程和Taguchi实验设计方法[14-16]的基础上,提出了一种LED环阵列设计方法,其中Taguchi方法是一种在不同领域解决复杂优化问题的实验设计方法。所提方法的目的是通过优化光源相关参数设置来提高机器视觉照明系统的一致性。通过一系列光线跟踪仿真和实验,验证了所提方法在机器视觉光学优化设计方面的有效性和优越性。
2 成像理论推导
在设计实验之前,需要定义用于机器视觉系统的LED环形阵列的光学性能指标和测量规则。
2.1 成像中的反射模型
彩色相机电子成像具有线性关系,三色通道(k=1,2,3)相机传感器在像素处的响应可以表示为
式中:fk(λ)是第k个通道的颜色过滤器的光谱透射率;d(λ)是探测器的光谱灵敏度;R(λ)是拍摄对象的光谱反射率;ls(λ)是成像光源的光谱辐射;nk是反应常数[17];λH和λL分别是高和低的极限波长。对于普通相机,滤光片的透光率、传感器的灵敏度及光源的亮度是未知的。因此,(1)式可以简化为
式中:参数mk(λ)=fk(λ)d(λ)ls(λ)。
对于光谱反射率,研究人员提出了一些基于物理光学的反射模型,本文选择Cook等提出的反射模型[18]。在这个模型中,任意方向的反射强度依赖于入射能量,而不仅仅依赖于入射强度。双极性反射率R的定义是反射光强度(Lr)与给定方向(在一个小立体角dωi)的入射能量(Ei)的比值。反射光线可以分为两部分,即镜面反射Rs和漫反射Rd。表面反射的几何构型如
式中:s(s>0)和d(d>0)均是反射的两个分量的权值,s+d=1;2δ是光和视觉方向的夹角;γ是入射光和视觉方向的夹角的平分线与曲面法线之间的夹角;Li是入射光强度。
基于能量和强度之间的关系,即Ei=Licos σdωi,σ是光传播方向和表面法向之间的角度,η个光源观察方向的反射光强度和反射能量可以分别被描述为
除了外设光源的照明外,环境光也会照亮物体。对于一个特定的观测方向,环境光的反射量非常小,环境反射率(Ra)与观测方向无关,在均匀漫反射的条件下,Rd=
式中:F(δ, n)、D(m, γ)、G(δ, γ)分别是菲涅耳反射系数、微面分布函数、几何衰减系数;n是物体的折射率,m是微面的方均根斜率,该二参数均取决于物体表面的特性。
2.2 成像灰度等级
相机成像后,传感器在像素处的响应用灰度级表示。对于彩色图像,有3×256个不同的灰度级。但在单色图像中,灰度级仅有256个。
在正常曝光条件下,图像灰度(Ggray)与曝光量(H)之间存在关系[20],表达式为
式中:v和b都依赖于成像条件,并且可以通过提取两个测试点来进行确定;曝光量取决于传感器上的照度(I)和曝光时间(T)。
将(8)式代入(7)式,得
式中:E是拍摄对象反射回的能量Er;S是物体表面的区域,对于确定对象,可认为是一个常数;τ是光学镜头的透光率;P是相机的光圈;b'和b″在化简过程中都是常数。然后将(5)式代入(9)式中,由于环境光在各个方向上是均匀的,它不依赖于成像位置,从而得到图像灰度级的影响因子,表达式为
式中:K是环境光反射能量。
基于(10)式,不同目标表面单元的灰度差可以表示为
表面单元的最佳对比度Cmax可以表示为
由(11)、(12)式可知,在非过曝光照条件下,Ii1与Ii2的差异较小(照度均匀度越高),对于一个确定的物体,灰度差可以较好地反映单元表面特性参数的不同,即图像获得较高的对比度。具有较理想对比度的光环境可以使机器视觉的精度保持在较高的水平。若拍摄对象的反射率相同(Rs1=Rs2),即只设置单一背景,无拍摄对象,投射在背景面上的均匀度可以用成像灰度均匀度表示。
2.3 照度分布推导
为了避免检测对象(如不规则圆形目标、小尺寸目标)上的阴影,光源对提高视觉检测性能具有重要意义[21]。假设一个LED光源是一个朗伯体发射,其照度分布是观察角度的余弦函数。照度分布I(θ)可表示为
式中:I0是入射角θ=0时的照度;ψ值取决于LED发射区与球形封装曲面的相对位置[22]。通常情况下,ψ值大于30,观测角度对于强度的影响较为明显[23]。ψ值可以通过角度θ1/2(θ1/2定义为入视角为0°时的一半照度值对应的视角)获得:
LED的中心位置(x0=rsin θcos φ, y0=rsin θsin φ, z0=rcos θ),距平面屏幕的距离是h。每一个点(x, y, z)对应的照度[24]表示为
式中:L0是LED的发光强度;φ为LED位置的方位角;r为环的半径。LED环形阵列照明如
式中:α为环的数量;β为一个环的LED数量。
在成像过程中,根据(15)、(16)式,在均匀光照条件下,得到不失真的清晰图像。LED单元应安装在光基片上,方便照明均匀性设计。在充分理解成像原理后,根据(10)~(12)式、(15)式和(16)式,需要结合成像系统的参数对LED环形阵列结构照明状况进行分析和研究,以获得更高的对比度。对于特定的拍摄对象,其表面参数(n和m)是确定的,所以其他参数是后续实验重点优化的对象。为了获得良好的相机视觉性能,LED环形阵列的参数优化是成像系统光学设计的核心。
3 获得高质量图像的光环境设计方法
为了获得清晰的物体特征,机器视觉的输出光需要从不同的方向发出,以均匀照射光斑。从
3.1 设计方法介绍
Taguchi方法是一种有效的实验设计工具[25]。该方法可以解决产品性能弱、质量低及成本高的工程问题。在此基础上,提出了一种LED环形阵列均匀照明的灯光设计方法,如
式中:Iave和Imax分别是测试平面上的平均照度和最大照度。对于均匀照明,平均照度与最大照度匹配,理想的均匀度比为1。根据相关的组合结果,可以设计出一个LED照明灯具,同时,可以建立一个实验装置用来比较不同光照条件下的机器视觉成像结果。
3.2 影响因子选择
根据成像理论部分的公式推导,影响成像过程的参数有很多。为了简化设计,一些参数应设为常数值。按照成像模型的规律,δ、γ、r'、θ、φ、α、β、r及ρ是所有与光源位置有关的参数,因此,主要优化参数应该从这些参数中选择。
通过TracePro软件进行一些仿真实验,初步探究参数对照度均匀度的影响。LED的尺寸设置为1 mm×1 mm×1 mm,LED平面上设置每个单元的光线数量为50000,平面屏幕的大小设置为100 mm×100 mm×1 mm,屏幕和环形阵列的距离是200 mm。通过控制变量法测试相关参数与照度均匀度的关系,测试中设置仿真的恒定参数是δ、γ、r'、θ、φ、α和ρ。
表 1. ρ=0°时,LED环形阵列照度均匀度分布
Table 1. Illumination uniformity distribution of LED ring array at ρ=0°
|
图 4. ρ=0°时单环不同半径的照度分布。(a)~(e)照度均匀度分别是23.85%, 22.50%, 20.35%, 17.79%, 15.51%
Fig. 4. Illumination distribution of different radius with ρ=0° in single ring. (a)--(e) Illumination uniformity is 23.85%, 22.50%, 20.35%, 17.79%, and 15.51%, respectively
根据之前的分析,可以发现以下结果。当α、β、r'及ρ确定时,其他参数是可以通过计算得出的,所以,可以选择单一环的LED数量β、相邻环之间的距离r'、环数α、辐射角ρ作为LED环形阵列照明设计的主要参数。将各参数的水准值定在一个合适的范围内,如
表 2. 机器视觉光环境影响因子及其水准设置
Table 2. Influence factors and level setting for machine vision light environment
|
4 实验结果和分析
4.1 模拟分析
选取控制因子和水准设置后,在直交表中设计组合仿真实验,如
式中:yi'是与第i'次实验对应的照度均匀度;s'是实验总次数。
表 3. 仿真模拟直交表
Table 3. Simulation of direct intersection table
|
一个特定的参数对照明均匀性的影响可以通过计算该特定参数的三个水准相应的S/N值的平均值来评估。不同层次的S/N值如
表 4. 不同水准值对应的S/N值
Table 4. Corresponding S/N value of different level values
|
图 5. 组合A3-B1-C1-D1对应的归一化后的照度分布和中间线的照度分布。(a)照度分布;(b)中间线的照度分布
Fig. 5. Normalized illumination distribution and illumination distribution of central line corresponding to the combination of A3-B1-C1-D1. (a) Illumination distribution; (b) illumination distribution of central line
进一步,为了估量各参数对照度均匀性影响的贡献度,采用ANOVA方法进行分析。ANOVA方法的主要目的是评估实验误差并以百分比的方式进行表征。在Taguchi方法中,一般将交互作用视为误差的一部分,当参数效应相比实验误差足够大时,才认定该参数效应为“有意义的”或是“重要的”。有了实验误差后,可以评估每一个参数效应相对于实验误差的“重要性”。只有“重要性”达到某一程度的参数效应,才会被用在预测的公式中;其他“重要性”没达到该程度的参数效应则视为可能只是实验误差造成的,不进行进一步讨论。“重要性”的强弱通过参数ρ'表示:
式中:Sd是方差和;Se是错误的方差和。在本次分析中,Se接近0,这是因为实验运算的错误率极低,具有可重复性。方差St由S/N值的方差和计算得到:
式中:q是单一水准的实验数;εj'是每个因子第j'次实验的S/N值;εq是S/N值的平均值。
采用ANOVA方法来估计设计参数与评价目标之间的因果关系,较高水平的显著性说明计算的因子对评价目标有重要影响,参考(19)~(21)式。
表 5. 均匀性的ANOVA结果
Table 5. ANOVA results of uniformity
|
均匀性主要由单环LED数和照射角度决定。增加单环的LED数量可以改善均匀性,然而,这对应着散热器需要更大的尺寸,灯具的成本也会增加。此外,灯具的尺寸会限制LED的总数。由于LED数量的限制,辐照角度是唯一可以修改以达到高均匀性的参数。
图 7. 不同照射角度下的光强值。(a) ρ=0°;(b) ρ=20°;(c) ρ=40°
Fig. 7. Light intensity at different irradiating angles. (a) ρ=0°; (b) ρ=20°; (c) ρ=40°
图 8. 20°照射角下对应的归一化照度。(a) 2D照度分布图;(b)中间线照度分布图;(c) 3D照度分布图
Fig. 8. Normalized illumination corresponding to irradiating angle of 20°. (a) 2D map of illumination distribution; (b) illumination distribution map of central line; (c) 3D map of illumination distribution
4.2 实验验证与对比
为了验证所提方法的有效性,制作了一个如
为了说明在优化光照条件下捕获图像的对比度,利用MATLAB软件将捕获的三维样本图像转换为灰度图,如
表 6. 测试灯具参数
Table 6. Parameters of test luminaire
|
图 9. 实验系统示意图。(a)设计的LED环形阵列;(b)实验装置
Fig. 9. Schematic of experimental system. (a) Designed LED ring array; (b) experimental setup
图 10. 不同照明环境下的3D样品成像效果。(a)样品图像;(b)最优光环境下拍摄的图像;(c)~(f)非均匀光环境下的图像
Fig. 10. 3D sample imaging effect under different illuminating conditions. (a) Sample image; (b) image taken in optimal light environment; (c)--(f) images in non-uniform light environment
图 11. 样品细节图。(a)最优光环境下的3D样品图像;(b)对应的灰度等级图
Fig. 11. Sample details. (a) 3D sample image in optimal light environment; (b) corresponding gray level map
为了测试灯具的照度均匀性,根据(11)式进行分析,采用设计的灯具对纯黑色不透明布(反射率为0.03)进行照明。在相同条件下由Genie Nano M2420单镜头相机拍摄亮度图像,结果如
式中:V和W分别是图像沿垂直方向和水平方向的像素数;iv和jw分别是像素沿水平方向和垂直方向的位置。并利用MATLAB软件对图像的灰度值均匀性进行分析计算。均匀度U由平均灰度值和最大灰度值计算得到。
优化后的照明灯具的灰度分布如
图 12. 最优光环境下对应的光斑和灰度分布。(a)光斑;(b)中心区域的3D灰度等级分布
Fig. 12. Corresponding light spot and gray distribution in optimal light environment. (a) Light spot; (b) 3D gray level distribution in central area
5 结论
提出了一种用于机器视觉的高质量成像LED环阵列结构的优化照明设计方法。首先,利用灯具结构参数与反射模型、图像灰度对比度之间的关系,找出影响图像质量的主要参数。在此基础上,提出了一种基于Taguchi实验设计方法的对LED环形阵列照度均匀性进行优化的方法,该方法很容易得到最优的设计参数组合。对模拟结果进行方差分析,得到各因子对性能的贡献,结果表明,单环LED数(贡献度为53.10%)对均匀性的影响最大,其次是辐照角度(贡献度为32.85%),照度均匀度达89.31%。然后,设计相关实验加以验证,实验结果表明,所提LED环形阵列光优化方法能够生成高均匀度的三维样品图像,同时照度均匀度可达80.35%。通过仿真和实验分析,可知所提方法可以高效快速地解决机器视觉无法获得高质量图像的光学设计问题。机器视觉进行光学检测准确性时需要高质量图像作为依据,所提方法通过优化照明设计,实现了均匀照度,还原了样品表面真实的对比度信息,可有效解决这方面的问题。在未来工作中,所设计的机器视觉光源可以满足更多针对不同特征对象的需求。
[1] 尹仕斌, 任永杰, 刘涛, 等. 机器视觉技术在现代汽车制造中的应用综述[J]. 光学学报, 2018, 38(8): 0815001.
[2] 梁秀满, 刘文涛, 牛福生, 等. 基于机器视觉的浮选气泡体积和表面积测量研究[J]. 光学学报, 2018, 38(12): 1215009.
[3] 屈恩世, 张恒金, 曹剑中, 等. 对光学设计中照度计算公式的讨论[J]. 光学学报, 2008, 28(7): 1364-1368.
[4] Tarabanis K A, Allen P K, Tsai R Y. A survey of sensor planning in computer vision[J]. IEEE Transactions on Robotics and Automation, 1995, 11(1): 86-104.
[5] Yi S K, Haralick R M, Shapiro L G. Optimal sensor and light source positioning for machine vision[J]. Computer Vision and Image Understanding, 1995, 61(1): 122-137.
[6] Zhu Z M, Liu H, Chen S M. The design of diffuse reflective free-form surface for indirect illumination with high efficiency and uniformity[J]. IEEE Photonics Journal, 2015, 7(3): 1-10.
[8] Pimputkar S, Speck J S. DenBaars S P, et al. Prospects for LED lighting[J]. Nature Photonics, 2009, 3(4): 180-182.
[9] 钱可元, 胡晓佳. 同时实现路面照度与亮度均匀性的LED光源光学系统研究[J]. 光学学报, 2015, 35(2): 0208001.
[10] Dong J T, Lu R S, Shi Y Q, et al. Optical design of color light-emitting diode ring light for machine vision inspection[J]. Optical Engineering, 2011, 50(4): 043001.
[12] MarchandE. Control camera and light source positions using image gradient information[C]//Proceedings 2007 IEEE International Conference on Robotics and Automation, April 10-14, 2007, Roma, Italy. New York: IEEE, 2007: 417- 422.
[13] Li H T, Chen S C, Han Y J, et al. A fast feedback method to design easy-molding freeform optical system with uniform illuminance and high light control efficiency[J]. Optics Express, 2013, 21(1): 1258-1269.
[14] TaguchiG, YokoyamaY, WuY. Taguchi methods, design of experiments[EB/OL]. [2020-06-01].https://www.researchgate.net/publication/292446196_Taguchi_Methods_Design_of_experiments.
[16] 史晨阳. 基于Taguchi方法的曲面底板LED阵列光学设计与优化[D]. 广州: 华南理工大学, 2016.
Shi CY. Optical design and optimazation of cured surface LED array based on Taguchi method[D]. Guangzhou: South China University of Technology, 2016.
[17] Shen H L, Xin J H. Spectral characterization of a color scanner by adaptive estimation[J]. Journal of the Optical Society of America. A, Optics, Image Science, and Vision, 2004, 21(7): 1125-1130.
[18] Cook R L, Torrance K E. A reflectance model for computer graphics[J]. ACM SIGGRAPH Computer Graphics, 1981, 15(3): 307-316.
[19] Wu X J, Gao G M. LED light design method for high contrast and uniform illumination imaging in machine vision[J]. Applied Optics, 2018, 57(7): 1694-1704.
[20] 黄珂. 道路照明测量方法研究[D]. 重庆: 重庆大学, 2006.
HuangK. Study on measurement methods of road lighting[D]. Chongqing: Chongqing University, 2006.
[21] Wu F P, Li S P, Zhang X M, et al. A design method for LEDs arrays structure illumination[J]. Journal of Display Technology, 2016, 12(10): 1177-1184.
[22] Moreno I, Muñoz J, Ivanov R. Uniform illumination of distant targets using a spherical light-emitting diode array[J]. Optical Engineering, 2007, 46(3): 033001.
[23] Moreno I, Avendano-Alejo M, Tzonchev R I. Designing light-emitting diode arrays for uniform near-field irradiance[J]. Applied Optics, 2006, 45(10): 2265-2272.
[24] 史晨阳, 文尚胜, 陈颖聪. 基于Taguchi方法的曲面LED阵列照度问题研究[J]. 发光学报, 2015, 36(3): 651-656.
Shi C Y, Wen S S, Chen Y C. Study on curved surface LED array illumination problem based on taguchi method[J]. Chinese Journal of Luminescence, 2015, 36(3): 651-656.
[25] Lin C F, Wu C C, Yang P H, et al. Application of Taguchi method in light-emitting diode backlight design for wide color gamut displays[J]. Journal of Display Technology, 2009, 5(8): 323-330.
[26] Ramane D, Shaligram A. Optimization of multi-element LED source for uniform illumination of plane surface[J]. Optics Express, 2011, 19(14): A639-A648.
Article Outline
史晨阳, 林燕丹. 机器视觉中均匀照明成像的LED环形阵列设计[J]. 光学学报, 2020, 40(20): 2011003. Chenyang Shi, Yandan Lin. LED Ring Array Design for Uniform Illumination Imaging in Machine Vision[J]. Acta Optica Sinica, 2020, 40(20): 2011003.