高精度激光定位散斑机理分析及抑制
1 引 言
激光具有良好的单一方向性、高能量性和高相干性,常作为高精度视觉定位测量系统中的成像载体[1-2]。激光照射成像介质时,成像介质内分子接收到激光能量后会发生受激运动,导致激光投影成像时存在随机复合散斑[3-4],这会影响成像光斑的稳定性,进而增加光斑质心的定位误差。因此,抑制随机散斑是从根源上提高视觉定位测量系统中光斑质心提取精度的关键。
近年来,许多学者针对散斑抑制方法开展了研究。王锐[5]等提出建立多光束分光照明系统对激光散斑效应进行抑制。张圣涛[6]使用多个互不相干光源,利用半导体激光二极管阵列的合理排布,引进角度多样性降低激光散斑。这类方法需要改进系统光源的功率或数量,会导致光斑过饱和,增加光斑质心的提取难度。梁传样[7]等使用衍射光学元件实现激光整形和散斑抑制。该方法改变了激光成像的高斯型,很难完成光斑质心的高精度提取。Deng[8]等采用旋转球透镜来抑制激光投影散斑。付瀚毅[9]提出转动随机透镜阵列抑制散斑的结构。孙鸣捷[10]等提出对激光照明系统中混光棒施加连续作用力,使混光棒连续形变从而实现激光散斑抑制。这种方法虽抑制了散斑,但也会让成像光斑位置发生物理偏移,增大成像光斑质心的提取误差。上述方法虽能抑制散斑,但在视觉定位测量系统中,光源的光强与波长、成像靶面的位置都需动态改变。这种动态改变会降低光斑质心的提取精度,所以不适用于以光斑质心定位为基础的视觉定位测量系统。鉴于此,本文主要针对视觉定位测量中系统保持相对静止时的随机散斑进行抑制,提高成像光斑的稳定性及光斑质心的定位精度。
本文理论分析了视觉定位测量系统中随机散斑产生的原因,推导出空间中任意一点总波扰表达式。定量分析成像介质密度、粗糙度、吸光性、激光功率与成像光斑图像信噪比之间的关系,并在此基础上设计验证实验。提取不同物理特性成像介质上的光斑质心,获得光斑质心位置稳定性随成像介质物理特性变化的趋势,验证理论分析的正确性。
2 原 理
2.1 随机激光散斑的成因
图1是以激光为载体的高精度视觉定位测量系统。理论上,当相机和成像靶面相对静止时,成像系统记录的激光光斑应处于较稳定的状态。但实际使用这类视觉定位测量系统时,在相机、成像介质和激光器位置保持不变时,相机获得的光斑成像图稳定性却较差,严重影响光斑质心位置的提取精度。
稳定激光器输出的激光模式为高斯光束,其光场分布可用麦克斯韦方程组来表示[11]。结合麦克斯韦方程组与介质方程,得到自由空间稳态传输的电磁场所满足的亥姆霍兹方程:
式中:分别为坐标轴的单位矢量,是波数,是振幅。在近轴条件下,激光器产生的高斯光束是方程(1)的一个特解,可表示为:
式中:为最大幅值, 为空间中任意一点,可以看成是一个沿轴缓慢变化的复函数,为虚部。将式(2)代入式(1),忽略项得:
式(3)为抛物线型微分方程,用试解法求解,最后得到高斯光束的表达式:
式中:表示高斯光束的等相面曲率半径, 表示光斑半径,表示附加相移。式(4)含振幅和相位两部分,它反映高斯光束场分布及光束在传播过程中的变化规律。
将式(4)代入惠更斯-菲涅尔公式中,则空间中位置矢量为的散射微元所产生的散射子波扰为:
当空间中子波叠加时,空间中任意一点总波扰可表示为:
式中:为倾斜因子,主要用来调整子波朝着各个方向传播波幅;为空间中某位置矢量;为子波源到点的位置矢量;是散射积分曲面;为某位置矢量总波扰。
图2为折射光在成像介质内部碰撞的过程。由分子运动论、光的折射定律和光的波粒二象性定理可知,用激光照射成像介质,折射光会进入到成像介质内部,内部部分分子吸收其他分子传递的能量后,不只在平衡位置振动,从而导致成像介质内部分子始终在随机运动[12-13]。因此,散射面所产生的散射子波在方向和强度上随机。由式(6)可知,如果散射面射出的子波在强度和方向上发生随机变化,那么总波扰的强度也会发生变化。这些随机子波在空间中相互干涉,经过光学系统后,在像面上形成许多随机分布的明暗动态变化图样,出现随机激光散斑现象。
2.2 成像介质表面出射光信噪比分析
由2.1节分析可知,成像介质表面随机散射光导致随机激光散斑的出现。因反射光不会进入成像介质内部,所以在不考虑空气扰动的情况下,反射光本身是稳定的。以激光为光源,从成像介质表面进入到空间中的反射光和散射光可表示为:
式中:表示纯反射光,表示反射光中与散射光发生干涉的部分,表示纯散射光,和成正比。可见光斑的不稳定性主要来自散射光自干涉,以及散射光与反射光的交叉干涉。图3表示了3种光的分布情况,3部分扇形内线条越凌乱说明其稳定性越差。
图 3. 纯反射光、混合光和纯散射光示意图
Fig. 3. Schematic diagram of pure reflected light, mixed light and pure scattered light
因相机视场范围有限,并非所有从成像介质表面出来的光都能进入到相机中。根据式(7)和式(8),像面上记录的能量可表示为:
式中:为视场范围,是相机像面中记录的激光总信号。决定最终像面中光斑的整体亮度。若将定义为稳定的信号,将定义为噪声,则信噪比可表示为:
由式(10)可知,为提高信噪比,可以增大,但这会使值过大,成像光斑过亮,不利于提取质心。因此,只能在减小的基础上提高信噪比。
2.3 成像介质粗糙度对光照的影响
根据微平面理论,在微观尺度下,任何表面都是由若干个细小平面组成[14]。微表面粗糙程度可表示为:
式中:为采样点个数,为采样点与微平面平均高度的差值。
由双向反射分布函数(Bidirectional Reflected Distribution Function, BRDF)的定义可知,只有当微平面法线与的半程向量方向一致时[15],相机才能接收到微平面反射光。如图4所示(彩图见期刊电子版),虽然A,B两点都能被相机和光照到,但只有B点的光能反射到相机中。此外,微表面法线分布的不规则性也会导致散射光从成像介质内部出来的方向更加随机。
Trowbridge等在式(11)基础上,使用高斯分布来近似统计微表面的法线分布与半程向量重合的比率[16],即:
式中:为粗糙度,为微平面法线,为光位置矢量,为视线位置矢量,为半程向量。
式(12)是通过统计微平面法线分布与半程向量重合度来估计进入到相机像面上的能量,但这种关系并非绝对成立。因为粗糙微表面会导致微平面间相互遮蔽,这种遮蔽会损耗光线能量[17]。图4蓝色入射光展示了微平面对入射光和出射光的遮蔽现象,从图4可以看出,微平面遮挡关系依赖于相机视线方向和入射光方向。
为量化微平面这种不规则性排列对光能量的影响,Smith利用统计学理论,近似求得微平面之间未被遮挡部分的比例,公式[18]如下:
式中:为射线位置矢量,为粗糙度,为微平面法向量。
考虑到粗糙度会影响反射和入射光线能量,由式(13)可知,任意微平面未被遮挡的概率可表示为:
式中:为入射光位置矢量,为视线位置矢量。
结合式(12)和式(15)可知,成像介质粗糙度对入射光、散射光和反射光的影响可表示为:
式中:为单束入射光或散射光能量,为最后进入到相机像面的能量。
综上所述,增加粗糙度会导致散射光和镜面反射光的自损耗,并提高超出相机视场的概率,引起减少,降低,进而提高光斑稳定性,减弱像面上记录的光斑亮度。虽然增加粗糙度能抑制散射光在空间中的传播,但反射光也会在一定程度上被抑制。此外,过大的粗糙度导致像面上的能量过小,增加光斑质心提取的难度,因此粗糙度要适中。
2.4 成像介质吸光性对光照的影响
根据朗伯比尔定律,物质本身可以吸收光子,物质的颜色是对不同波长可见光选择性吸收的结果。某物质吸收了大部分自然光,那么该物质呈现黑色,反之则表现为白色[19]。物质对光的吸收度可表示为:
式中:为摩尔吸光系数,为吸收层厚度,为吸光物质浓度。
由式(17)可知,相比于同浓度、同厚度的白色物质,暗色物质有更强的吸光性,所以用加厚的暗色物质作为成像介质,能让更多光子在成像介质内部消耗,使成像介质表面射出的散射子波数量更少、能量更小,即抑制的产生,提高光斑稳定性,减小成像平面的光斑亮度。结合2.2节的分析,粗糙度和吸光性均会一定程度地抑制散射光,选择合适粗糙度和吸光性的成像介质,能提高成像介质表面出射光的信噪比,抑制视觉定位测量系统中随机激光散斑的产生。
3 实 验
本文从成像介质密度、表面粗糙度、吸光性和激光输出功率4个方面开展实验。在成像介质实验中,为排除背景光对漏孔靶面光斑提取的影响,实验在暗室进行。其他实验在自然光下完成。实验装置包括相机、全站仪和成像靶面。选取具有不同物理特性的成像介质作为激光照射的靶面,来验证靶面与散射光强度之间的关系。以全站仪的激光器作为激光光源射向靶面,相机仰视看向靶面并将采集的数据实时显示在电脑屏幕上。实验过程中保持相机、激光器和成像靶面的位置不变,激光照射方向可变。其中,相机分辨率为1 2801 024,帧率为250 frame/s,激光器功率在5~200 mW之间,激光器距离靶面10 m。本文使用一种改进的高斯拟合法[20]对靶面中的激光光斑进行质心提取,以坐标稳定性来评判不同材质的成像靶面对激光散斑的抑制能力。为了防止激光光斑能量过低时光斑提取受到外界环境干扰,或光斑能量过高出现饱和现象,本文将激光器的输出功率设定为50 mW,并采集连续4 s内的1 000张照片进行质心提取。为评价光斑稳定性,使用贝塞尔公式计算光斑质心标准差,然后统计光斑质心的均值(Mean)、极差值(Range)和标准差(STD)。
3.1 材料密度数据对比实验
由2.1节分析可知,激光照射处靶面内部分子的随机运动导致激光随机散斑现象。为验证该结论,本实验选取材质相同、漏孔尺寸分别为0.1,1,2 mm的白色尼龙布(WNC)作为靶面。图5中方框内是成像靶面上的激光光斑。图5左下角为靶面材料。为展示成像光斑的细节,将实验中靶面上的光斑放大12倍,图5右上角为靶面成像光斑放大图。后续的靶面材料粗糙度实验和吸光性实验中,对靶面成像光斑进行相同的放大处理。随着靶面密度的减小,会有更多的光透过靶面,即单位面积内产生的散射子波更少,在成像靶面上记录的光斑稳定性增强,亮度降低。具体实验数据如表1所示。
表 1. 三种密度靶面激光质心坐标的均值,极差值和标准差
Table 1. Mean, range, and STD of laser centroid pixel coordinates on three kinds of density target surfaces
|
从表1可知,漏孔越大,靶面密度越小,光斑质心坐标的极差值和标准差越小,光斑越稳定。2 mmWNC靶面上光斑质心坐标的极差值和标准差最小,光斑质心坐标提取的标准差在0.009 pixel内。靶面密度与光斑稳定性成反比,与理论分析结果一致。但漏孔靶面不适用于实际测量,因为大多数实验在多光源环境下进行,背景光会影响靶面光斑提取。
3.2 粗糙度数据对比实验
当光照环境、相机和激光器位置以及激光功率不变时,粗糙度增加导致散射光和镜面反射光自损耗的增大,超出相机视场的概率也会相应增加。为验证粗糙度对激光散斑的影响,本实验选取黑色羊绒布(BCC)、黑色玻纤幕布(BFC)、黑色光面亚克力板(BAC)作为靶面,如图6所示,3种材料的表面粗糙度依次下降。图6中方形框内是成像靶面上的激光光斑,右上角为靶面成像光斑放大图,右下角为靶面材料。
由2.3节的理论分析可知,在一定范围内,粗糙度越大,光斑稳定性越高,亮度越弱。具体数据如表2所示。
表 2. 三种粗糙度靶面激光质心坐标的均值,极差值和标准差
Table 2. Mean, Range, and STD of laser centroid pixel coordinates on three kinds of roughness target surfaces
|
由表2可知,粗糙度最大的BCC靶面的光斑坐标极差值和标准差最小,光斑稳定性最高,光斑质心坐标的标准差在0.084 pixel内,所得到的实验结果和理论分析一致。
3.3 吸光性数据对比实验
因物质对光的吸收能力与物质本身的摩尔吸光系数以及厚度成正比。为验证吸光性对散射光的影响,本实验选取白色羊绒布(WCC)、黑色羊绒布(BCC)和加厚版黑色羊绒布(TBCC)作为靶面。3种靶面均为同种材质,且密度相同。如图7所示,方形框内是成像靶面上的激光光斑,右上角为靶面成像光斑放大图,右下角为靶面材料。
图 7. 三种吸光性材质靶面光斑图像
Fig. 7. Images of light spot on target surface of three light absorbing materials
根据式(16)可知,3种靶面吸光性从左至右逐渐增强,所以加厚版黑色羊绒布上的光斑应最稳定且光斑亮度最低,实验数据如表3所示。
表 3. 三种吸光性材质靶面激光质心坐标的均值,极差值和标准差
Table 3. Mean, Range, and STD of laser centroid pixel coordinates on three light absorbing target surfaces
|
由表3可知,TBCC靶面的光斑稳定性最高,光斑质心坐标的标准差在0.031 pixel内,这个结果与理论分析的结果一致。
3.4 激光功率数据对比实验
由表2和表3中的实验数据可知,加厚黑色羊绒布的标准差最小。为测试该材质对动态散斑的抑制范围,本实验将激光器的输出功率分别调为10,50,70,90,110,130,170和200 mW,并计算了质心坐标的均值、极差值和标准差,实验数据如表4所示。图8为8种功率下靶面成像光斑放大13倍后的细节图。图9为8种功率下光斑质心STD的变化曲线。
表 4. 八种功率激光的质心坐标均值,极差值和标准差
Table 4. Mean, range, and STD of laser centroid pixel coordinates on eight kinds of power lasers
|
由表4和图9可知,随着功率的增加,光斑稳定性逐渐降低。激光功率在10~70 mW时,光斑质心的STD较为稳定,而且功率的增加并未明显影响其稳定性。然而,当激光功率大于70 mW时,光斑逐渐呈现饱和状态,光斑质心计算呈现较大的随机性。因此,实际应用时应在激光光斑不出现饱和的状态下选择最佳材质。随着入射光能量的增加,随机散射光的数量和能量也会增加。由2.2节对信噪比的分析可知,增大,光斑不稳定性会增强,并且也会提升,使光斑出现饱和状态,理论分析结论和实验结果一致。
4 结 论
本文针对高精度视觉定位测量系统中激光随机散斑所引起的光斑不稳定性,进而导致定位误差的问题开展研究,结合激光发光原理和分子运动论,分析出随机散斑是成像介质内分子受激后进行无规则运动,进而导致散射子波在空间中随机干涉所产生的现象,这种现象与散射子波的强度以及方向有关。从图像信噪比角度分析散射光和反射光,得出只有抑制散射光才能提高光斑稳定性和质心提取准确性的结论。在此基础上,使用粗糙度适当和吸光性的成像介质能有效抑制散射光的出现。为验证理论分析的正确性,从成像介质表面密度、粗糙度、吸光性和激光器输出功率4个方面展开实验,得到的实验结果与理论分析的结论完全一致。实验数据表明,激光功率在50 mW左右时,选择STD最小的加厚黑色羊绒布材质作为成像靶面,光斑稳定性在0.04 pixel内,能较好地抑制随机散斑,提升光斑质心的定位精度。
[1]
[3]
[4]
[5]
[7]
[8] DENG L X, DONG T H, FANG Y W, et al. Speckle reduction in laser projection based on a rotating ball lens[J]. Optics & Laser Technology, 2021, 135: 106686.
[9]
[12] PAUL B, KAMAL S, QURESHI T. Multipath wave-particle duality in classical optics[J]. Optics Letters, 2020, 45(12): 3204.
[14] CHENG J, CHENG X L, MENG X C, et al. A Monte Carlo emissivity model for wind-roughened sea surface[J]. Sensors, 2019, 19(9): 2166.
[15] SUN L L, ZHAO F. Bidirectional reflectance distribution function algorithm based on the Poynting vector analysis[J]. Optical Engineering, 2021, 60(6).
[16] TROWBRIDGE T S, REITZ K P. Average irregularity representation of a rough surface for ray reflection[J]. Journal of the Optical Society of America, 1975, 65(5): 531-536.
[17] BUTLER S D, NAUYOKS S E, MARCINIAK M A. Comparison of microfacet BRDF model to modified Beckmann-Kirchhoff BRDF model for rough and smooth surfaces[J]. Optics Express, 2015, 23(22): 29100.
[18] IGOUCHKINE O, ZHANG Y B, MA K L. Multi-material volume rendering with a physically-based surface reflection model[J]. IEEE Transactions on Visualization and Computer Graphics, 2018, 24(12): 3147-3159.
[19] MAYERHÖFER T G, PAHLOW S, POPP J. The bouguer-beer-lambert law: shining light on the obscure[J]. Chemphyschem: a European Journal of Chemical Physics and Physical Chemistry, 2020, 21(18): 2029-2046.
[20]
Article Outline
徐晶, 王钎潮, 任正玮, 付飞蚺, 方明. 高精度激光定位散斑机理分析及抑制[J]. 光学 精密工程, 2023, 31(13): 1900. Jing XU, Qianchao WANG, Zhengwei REN, Feiran FU, Ming FANG. Mechanism analysis and suppression of speckle in high-precision laser location[J]. Optics and Precision Engineering, 2023, 31(13): 1900.