Nonuniformity analysis and correction of DMD imaging in infrared target simulation system
-
摘要: 基于数字微镜阵列(DMD)的红外目标模拟系统具有高分辨率、高帧率和大温度范围等优点,能够经济、有效地检测和评估红外探测系统的工作性能。为了研究基于DMD的红外目标模拟系统中照明系统的结构设计对成像均匀性的影响并减小像面照度不均匀性,提出了一种消除成像不均匀的方法,对DMD照明系统的结构参数进行了分析和设计,以成像光束截面积之比作为相对照度指标,建立了数学模型,用软件进行数值仿真,得到了能消除因渐晕导致的照度不均匀现象的结构参数和像面各点相对照度分布在88.75%~1;并用TracePro验证了理论分析的结果,得到像面各点相对照度分布在90.91%~1,与理论分析结果相符,证明了在DMD照明系统结构参数设计合理的情况下,该方法可以有效降低因渐晕导致的像面照度不均匀性,为红外目标模拟系统结构参数的设计和像面照度分析提供了参考。Abstract: The infrared target simulation system based on the digital micromirror device (DMD) has the advantages of high resolution, high frame rate and large temperature range, which can detect and evaluate the working performance of the infrared detection system economically and effectively. In order to study the influence of the structure design for illumination system on the imaging uniformity and reduce the illumination uniformity of the image surface in the infrared target simulation system based on the DMD, a method to eliminate the imaging nonuniformity was proposed. In which the structure parameters of the DMD illumination system were analyzed and designed, the ratio of the imaging beam cross-sectional area was taken as the relative illumination index, the mathematical model was established with the software for numerical simulation, and the structure parameters which could eliminate the illuminance nonuniformity caused by the vignetting was obtained. The relative illuminance distribution of each point on the image plane was ranged from 8.75%~1; the results of the theoretical analysis were verified by the TracePro, and the relative illuminance distribution of each point on the image plane was ranged from 90.91%~1, which was consistent with the theoretical analysis results. It was proved that in the case of the rational design of the structure parameters for DMD illumination system, the image plane illuminance nonuniformity caused by the vignetting could be reduced effectively, which provided a reference for the design of structural parameters for infrared target simulation system and image plane illuminance analysis.
-
引言
激光三角法作为一种非接触式测量方法,由于其具有材料适应性广、不易损伤表面、结构简单、测量范围大、抗干扰性强、测量准确度高、可用于实时在线测量等特点,已经在三维轮廓、厚度、料位、位移及振动测量中得到广泛应用[1-4]。
激光三角法是通过将光源发射出的激光束投射到被测表面,在另一方向通过成像观察反射光点的位置,从而计算出物点的位移或被测表面高低变化量[5]。由于入射光和反射光构成一个三角形,所以这种方法被称为三角测量法,又可按入射光线与被测工件表面法线的关系分为直射式和斜射式。并且由于直射式的优点是光斑较小,光强集中,不会因被测面不垂直而扩大光斑,一般体积较小[6]。
作者在前期工作中,以传统的直射式激光三角法为基本光学原理,以一字线激光器为入射光源,建立了一种表面粗糙度非接触测量系统,并结合数字图像处理技术实现粗糙度的测量。然而对于表面粗糙度这种微观几何量的测量方法而言,其特点是测量范围较小而测量精度应该较高。因此本文重点针对这种光学测量系统的测量精度进行探讨。
研究人员从光源与探测器等器件、测量对象、光学结构形式等方面进行考虑并采用了多种措施方法提高激光三角法的测量精度,例如, 用光学双稳装置来稳定激光光源的光强[7];考虑被测表面倾斜对测量方法的影响并加以校正[8];位移传感器结构、光斑定位不确定度及定标方法对测量精度的影响[9];物体表面颜色及粗糙度对测量结果的影响及补偿[10]。采用双三角光路修正投影光束轴线偏摆和投影反射光斑重心偏移的方法[11]等。
事实上,还可通过提高仪器的测量灵敏度来改善系统的测量精度。对于直射式三角测量系统,它的测量灵敏度除了与放大倍率[12]有关外,还与仪器结构参数有着密切的联系,其中主要涉及物距和工作角。另外,物距变化导致放大倍率发生变化,因此物距与放大倍率相联系。文献[12]对此讨论后得出结论:物距应尽可能小,这既有利于提高成像系统放大倍率,同时也可缩小传感器的体积;若为平面被测物,在实现精密测量时,像距与物距比一般取1.2~2.5。而对于工作角这一结构参数对系统灵敏度影响的分析报道甚少。于是本文通过建立表面粗糙度激光三角测量系统,采用理论推导与实验方式分析了工作角与灵敏度之间的关系。
1 测量系统物像关系
直射式三角法测量光路如图 1所示。图中l和l′分别为参考面的物距和像距;α为激光光轴与成像光轴之间的夹角,定义为工作角;β是光接收器基线与成像光轴间的夹角,定义为成像角;z是被测物体相对于基准位置的位移量;x是像点相对于成像基准点之间的距离。
图 2为激光束、成像透镜及光接收器CCD三者间的位置关系图,图中O是被测物体的基准位置;P点是成像光点的基准点。为了保证测量的精度,α和β必须满足Scheimpflug条件[13]。
$$ l\text{tan}\alpha =l\prime \text{tan}\beta $$ (1) 即激光束轴线、成像透镜主面及焦平面三者的延长线交于一点,或者三者相互平行,这是激光三角测量传感器实现精密测量的前提条件。因此要求光接收器的受光面与成像光轴成一夹角β,此时一定景深范围内的被测点都能成像在探测器上。
若确定出像点在成像面上的位移x时,利用相似三角形比例关系,按下式求出被测面的位移:
$$ z=\frac{lx\text{sin}\beta }{l\prime \text{sin}\alpha -x\text{sin}(\alpha +\beta )} $$ (2) 对于本文的表面粗糙度三角法测量系统而言,输入量用表面起伏量z来表示,输出量则为光接收器上的像点的位移量x。如较小的表面起伏能转化为较大的输出偏移,便可提高系统测量精度。于是将(2)式重写为x与z之间的关系式:
$$ x=\frac{zl\prime \text{sin}\alpha }{z\text{sin}(\alpha +\beta )+l\text{sin}\beta } $$ (3) 2 测量灵敏度分析
灵敏度是传感器静态特性的一个重要指标,定义为输出量的增量与引起该增量的相应输入量增量之比,也就是传感器在稳态下的输出变化与输入变化的比值:
$$ ~S=\underset{\Delta z\to 0}{\mathop{\text{lim}}}\, \frac{\Delta x}{\Delta z}=\frac{\text{d}x}{\text{d}z} $$ (4) 对(3)式微分得到灵敏度表达式:
$$ S=\frac{\text{d}x}{\text{d}z}=\frac{ll\prime \text{sin}\alpha \text{sin}\beta }{{{[l\text{sin}\beta +z\text{sin}(\alpha +\beta )]}^{2}}} $$ (5) 显然,测量灵敏度与系统结构参数有着密切联系。
当成像透镜选择好以后,焦距便已知,并且根据系统结构确定出物距后便可得出像距;再当工作角给定后便可求得成像角β。因此β仅为导出参数,通过将(1)式代入(5)式而消去,由此灵敏度公式进一步简化为
$$ ~S=\frac{ll\prime \text{sin}\alpha }{{{(l+z\text{cos}\alpha )}^{2}}} $$ (6) 在表面粗糙度测量系统中,l远大于z。此外,对于图 1直射式结构,α为0°~90°,于是有0≤cosα≤1。这时(6)式的分母可近似表示成:
$$ l+z\text{cos}\alpha \approx l $$ (7) 于是重写(6)式为
$$ S=(l\prime /l)\text{sin}\alpha $$ (8) 由(8)式可知,测量灵敏度与系统放大倍率M=l′/l有关,并且测量灵敏度与放大倍率成线性正比关系。由此验证了文献[12]中“应尽可能减小物距以有利于提高成像系统的放大倍率,从而改善测量灵敏度”这一结论。
而灵敏度与α成非线性正弦关系,且工作角越大,灵敏度越高。但是,到达光接收器上的光功率将随α的增加而减小,因此应当在确保接收的光功率大于其灵敏度的前提下再尽量增大α,这将有利于提高仪器的灵敏度以改善测量精度。以下将讨论最佳工作角的确定。
根据近似前的(6)式可知,灵敏度与表面起伏量也有关系,在z=0点处的灵敏度是最高的;而当z越大时,测量灵敏度则降低。考虑这一点,三角法在平面物体粗糙度测量中具有更好的测量效果。
3 工作角的确定
将放大倍率l′/l设为常数,根据(8)式仿真得到图 3所示的关系曲线。当α较小时,灵敏度随α的增大而迅速提高;但当α增大到一定角度时,灵敏度提高,速率则明显变缓。并且为了产生较大的α,还需调整物距l和成像透镜的焦距f,从而导致放大倍率发生变化。因而在α较大的情况下,需根据(8)式分析放大倍率与工作角对测量灵敏度的综合影响。
选择图 4所示的表面光洁度等于$\nabla $3的刨床量具作为被测对象,建立图 5所示的实验系统,来验证工作角与灵敏度之间的关系、以及放大倍率与工作角之间的相互制约关系。测量系统建立在JX13C图像处理万能工具显微镜上,选择放大倍率为1倍的显微物镜;一字线激光器作为光源,输出波长为635 nm,直线度≥1/5 000。
考虑到图像采集系统与激光器的物理尺寸,工作角不能太小,实验中能实现的最小α约为30°。在保证系统放大倍率不变情况下,α可从30°增大到70°附近,同时拍摄不同工作角所对应的激光光带图,如图 6(a)~(c),3幅图对应的α分别为30°、50°和70°。为了进一步取得更大的工作角并保证清晰成像,需增大物距并相应地调整透镜焦距,图 6(d)为α=80°时的激光光带。
测量过程中使用同一标准样块,因此在一定长度范围内输入量的变化量Δz是不变的。根据灵敏度定义式(4)可知,输出量的变化量Δx越大,表明测量灵敏度越高。输出量x表现为光斑像点在探测器上的位移,于是增量Δx则对应于激光光带的起伏程度。图 6(a)~(c)表明,工作角从30°增加至70°时,激光光带起伏程度逐渐增大,由此说明灵敏度随工作角的增加而提高。而当工作角进一步增加至80°时,光带图像的起伏程度又逐渐减小,如图 6(d)所示,说明灵敏度出现下降趋势。究其原因,则是由于为取得更大的α并保证成像清晰,需通过增大物距并适当调整焦距从而导致系统放大倍率降低;这种情况下,(8)式中的sinα项尽管是增大的,但它对改善灵敏度的作用在减缓;另外由于(8)式中l′/l项在减小,于是在二者综合影响下导致灵敏度呈降低趋势。
4 结论
由三角法灵敏度公式分析出灵敏度包含3个影响因素,即被测表面起伏量、光学放大倍率与工作角,并给出了每一参数对测量灵敏度的独立影响结果。结合工作角的研究结果得出,提高测量灵敏度除了通过提高成像系统放大倍率、选用高分辨率探测器、尽可能在基准平面附近进行测量之外,还需增大工作角。但这不是任意的,将受两方面条件的限制:其一,工作角增大到较大角度时,为保证清晰成像,需动态调整物距和物镜焦距反而导致系统放大倍率降低。其二,当工作角逐渐增大时,光接收器接收到的散射光的强度将逐渐减小。因此,必须保证在接收器件所接收到的光强大于其灵敏度的情况下增大工作角,只有这样才能真正提高测量灵敏度。论文的工作可为仪器设计中测量灵敏度分析、结构参数设计提供借鉴。
-
表 1 符合假设的结构参数表
Table 1 Structure parameters according to hypothesis
mm d r ${y_0}$ ${{\textit{z}}_0}$ 180 40 60 140 180 40 64 140 160 42 60 140 142 54 74 166 180 32 46 100 180 35 50 110 170 38 55 128 -
[1] 王海鹏. 基于DMD的红外目标模拟光学系统研究[D]. 哈尔滨: 哈尔滨工业大学, 2016: 1-26. WANG Haipeng. Research on optical system of infrared scene projector based on digital mirror device[D]. Harbin: Harbin Institute of Technology, 2016: 1-26.
[2] 孙永雪. 基于DMD的红外双波段目标仿真器研究[D]. 哈尔滨: 哈尔滨工业大学, 2013: 3-25. SUN Yongxue. Research of dual-band infrared target simulator based on digital-micromirror device[D]. Harbin: Harbin Institute of Technology, 2013: 3-25.
[3] 姚园, 王德江, 徐正平, 等. 基于DMD的红外成像制导目标模拟器研究[J]. 激光与光电子学进展,2013(50):072302-1-072302-2. YAO Yuan, WANG Dejiang, XU Zhengping, et al. Design of target simulator based on DMD for infrared imaging guided system[J]. Laser & Optoelectronics Progress,2013(50):072302-1-072302-2.
[4] 钟兴, 张元, 金光. 大视场光学系统像面照度均匀性优化[J]. 光学学报,2012,32(3):0322004-1-0322004-2. ZHONG Xing, ZHANG Yuan, JIN Guang. Optimization of image illumination uniformity in large field of view optical system[J]. Acta Optic Sinica,2012,32(3):0322004-1-0322004-2.
[5] 姚园. 红外成像制导目标仿真器的设计[D]. 哈尔滨: 哈尔滨工业大学, 2011: 1-22. YAO Yuan. Optical design of target simulator for infrared imagine guided system[D]. Harbin: Harbin Institute of Technology, 2011: 1-22.
[6] 庞广宁. 基于DMD的红外目标模拟系统关键技术研究[D]. 长春: 长春理工大学, 2018: 1-27. PANG Guangning. Key technology on DMD based infrared scene simulation[D]. Changchun: Changchun University of Science and Technology, 2018: 1-27.
[7] 李丹, 薛芸芸, 曹雯, 等. 基于微透镜阵列的DMD芯片投影系统照明优化[J]. 光学学报,2013,33(1):0122002-1-0122002-4. LI Dan, XUE Yunyun, CAO Wen, et al. Optimization of DMD illumination system with microlens Array[J]. Acta Optica Sinica,2013,33(1):0122002-1-0122002-4.
[8] 韩庆, 王健, 熊峥, 等. 用于长波红外目标模拟器的DMD衍射特性分析[J]. 红外与激光工程,2017,46(5):93-99. HAN Qing, WANG Jian, XIONG Zheng, et al. Diffraction characteristics analysis for DMD-based scene projectors in the long-wave infrared[J]. Infrared and Laser Engineering,2017,46(5):93-99.
[9] 李雄杰. 面向AMOLED屏亮度均匀性的光学检测系统[D]. 哈尔滨: 哈尔滨工业大学, 2017: 1-10. LI Xiongjie. An optical detection system for AMOLED screen brightness uniformity[D]. Harbin: Harbin Institute of Technology, 2017: 1-10.
[10] 杨洪飞, 夏晖, 马贝, 等. 红外数字TDI成像非均匀性校正[J]. 中国图像图形学报,2018,23(3):459-465. YANG Hongfei, XIA Hui, MA Bei, et al. Non-uniformity correction of infrared digital TDI imaging[J]. Journal of Image and Graphics,2018,23(3):459-465.
[11] 刘兰娇. 红外辐射源能量光谱分布测试研究[J]. 科技资讯,2014(33):42-43. LIU Lanjiao. Study on the measurement of energy spectrum distribution of infrared radiation source[J]. Science & Technology Information,2014(33):42-43.
[12] 赵亿坤, 王景辉, 原遵东, 等. 低温标准黑体辐射源的性能研究[J]. 计量学报,2017(5):589-592. ZHAO Yikun, WANG Jinghui, YUAN Zundong, et al. Study on the performance of the low-temperature blackbody radiation source[J]. Acta Metrologica Sinica,2017(5):589-592.
[13] 舒心, 郝小鹏, 宋健, 等. 100~400 K真空红外亮温标准黑体辐射源研制[J]. 计量学报,2019,40(1):13-16. doi: 10.3969/j.issn.1000-1158.2019.01.03 SHU Xin, HAO Xiaopeng, SONG Jian, et al. Research of 100~400 K vacuum infrared radiance temperature standard blackbody source[J]. Acta Metrologica Sinica,2019,40(1):13-16. doi: 10.3969/j.issn.1000-1158.2019.01.03
[14] 李卓, 牟达, 唐志方, 等. 高均匀性长波红外远心照明系统设计[J]. 长春理工大学学报,2017,40(5):9-12. LI Zhuo, MU Da, TANG Zhifang, et al. Design of high uniformity LWIr illumination system based on telecentric light path[J]. Journal of Changchun University of Science and Technology,2017,40(5):9-12.
[15] 付秀华, 马国俊, 吕景文, 等. 基于TracePro的开放型高反射回归器件的研究[J]. 光学学报,2019,39(11):1131003-1-1131003-1. FU Xiuhua, MA Guojun, LYU Jingwen, et al. Research on open high reflection regression device based on TracePro[J]. Acta Optica Sinica,2019,39(11):1131003-1-1131003-1.
-
期刊类型引用(2)
1. 杨丹丹,高礼华,沈悦,胡张斌,孟刚. 基于多学科交叉的国家重点实验室数字展厅设计创新. 包装工程. 2024(20): 500-512+517 . 百度学术
2. 宋涛,祁继辉,侯培国,赵明宇,李坤. 基于B样条曲线的投影图像边缘融合方法. 应用光学. 2021(06): 1040-1047 . 本站查看
其他类型引用(1)