An Infrared Micro Scanner Measurement and Calibration Method Based on Image Processing
-
摘要: 在以红外焦平面为核心的红外成像系统中,微扫描器件可以有效提高整个系统的空间分辨率。针对微扫描器件的检测,本文提出了一种基于图像处理的测量与校准方法,并搭建了一套检测系统用于对微扫描器件进行检测与校准。以某型微扫描器件为测试对象,实验结果表明本文所提方法在测量精度、重复精度以及不确定度方面均达到了良好的效果,可以为微扫描器件的设计、生产提供基础支撑。Abstract: In infrared imaging systems, in which the core is an infrared focal plane array, a microscanner can enhance the spatial resolution of the entire system. To test microscanners, this study developed a measurement and calibration method based on image processing and built a system to measure and calibrate microscanners. Using a microscanner as a test subject, the test results indicate that the proposed method has a significant effect on the measurement accuracy, repetition accuracy, and uncertainty. The method can provide technical support for the design and manufacture of microscanners.
-
Keywords:
- infrared focal plane array /
- micro scanner /
- measurement method /
- image processing
-
0. 引言
对于航天产品,无论是卫星平台还是有效载荷,都需要进行充分而适度的环境模拟试验,以达到发现隐患、优化设计,提高在轨可靠性的目的。在热控分系统中,包含一项极为重要、且贯穿整个研制阶段的环境模拟试验,热平衡试验。通过该试验,达到验证热设计的正确性和合理性、修正热分析数学模型的目的。在热平衡试验中,能否准确地模拟外热流直接决定试验的成败,是影响试验结果的关键热边界条件。对于近地轨道航天器,空间外热流主要由太阳直射、地球反照和地球红外3部分辐射能量组成[1]。
在热平衡试验中,模拟外热流的方式主要包括太阳模拟器、表面粘贴电阻加热片、红外加热笼、红外灯阵等。其中,太阳模拟器主要模拟与太阳光相接近的光谱分布、辐照度等,但是这种模拟方式成本高、技术复杂,需配套复杂的运动模拟器模拟航天器在轨的姿态变化。红外灯阵模拟的外热流辐照强度可以根据电压来调节,对瞬态的外热流模拟较容易实现,但是红外灯的辐照强度具有明显的方向性,不容易控制模拟表面的热流均匀程度,并且红外灯阵模拟外热流,需要考察模拟表面对红外灯的辐射光谱吸收特性。表面粘贴电阻加热片模拟外热流,模拟的是航天器表面的吸收外热流,要求加热器粘贴后不能对模拟表面造成损伤,并且要在加热器表面喷涂与模拟表面相同的涂层。针对红外加热笼模拟外热流,其以结构灵活、模拟简单等特点而被广泛应用在热试验中。红外加热笼模拟的目标外热流即为空间相机在热分析过程中所得到的轨道周期平均外热流。通常配合热流计,计算该热流密度下热流计敏感片的温度,并以热流计的温度作为反馈来控制红外加热笼的功率[2-5]。
杨晓宁[6-7]等通过优化红外灯阵不同的参数来降低模拟的热流密度的不均匀性。同时,还研究了红外加热笼带条的覆盖系数对到达卫星表面外热流均匀性的影响,仿真结果显示加热笼距离卫星表面越远,覆盖系数越大模拟出的外热流均匀性增长的越快。季琨等[8]针对热试验中红外加热笼之间的干扰进行了优化设计,通过防辐射屏和补偿边设计,可以达到减小热干扰的目的。韩继广等[9]研制出一种运动装置,可以改变红外加热笼与航天器表面之间的相对距离,来实现瞬态外热流的模拟。孙玉玮等[10]对红外加热笼的边缘效应进行了研究,通过增添边缘加热笼以及将边缘加热笼翻转一定角度的方式来改善航天器表面到达热流密度的均匀程度。此外,对圆台型航天器表面和与之相匹配的红外加热笼进行了研究,得出了圆台型加热笼的优化设计参数[11]。谢吉慧等[12]针对红外加热笼与航天器表面结构匹配问题,建立一套数字化结构匹配方法,可以使红外加热笼与航天器表面的匹配成功率达到100%。房红军等[13]对加热笼的布点进行了优化设计,通过对热流计布置位置的优化使外热流的误差均方根由9.2%降低至3.2%。但是调整布点位置改变的仅为实际到达卫星表面的平均外热流的大小,并不能改善到达卫星表面上外热流的均匀程度。
而对于特征尺寸较小的表面,采用红外加热笼模拟外热流时,可能存在总到达能量偏低、热流密度均匀性差等问题。因此,本文将基于某微型空间相机,以对外热流最为敏感的散热面为研究对象,开展红外加热笼仿真分析与优化设计的研究。采用有限元法建立红外加热笼-热流计的系统仿真模型,分析加热笼尺寸和热流计粘贴位置对模拟表面总到达能量和热流密度均匀性的影响。
1. 计算模型
在某微型空间相机的热平衡试验中,通过红外加热笼模拟散热面的到达外热流,在一种工况下,模拟出的外热流是恒定的。使用热流计粘贴在散热面的中心来测量模拟的到达外热流。实际试验外热流模拟装置示意图如图 1所示,相机整机包络尺寸为980 mm×200 mm×253 mm,散热面宽和高的尺寸为160 mm×167 mm;加热笼由框架、弹簧、钢带、底座组成,其中钢带通过弹簧与框架相连,形成闭合回路。加热笼的宽和高的尺寸为233 mm×255 mm,分别为加热笼左右边框的最大宽度和上边框至底板之间的距离,加热笼钢带的带条宽度为6 mm,厚度0.1 mm,带条间隔为9 mm;黑片热流计的外径为40 mm;红外加热笼与散热面之间的距离为50 mm。
为了分析该装置模拟出的外热流情况,采用UG/SST模块建立红外加热笼-黑片热流计系统的有限元模型,如图 2所示。其中,壳单元为15055个、节点18408个,共计33463个。对模型进行不同程度的简化,由于加热笼电阻丝与加热笼框架之间采用弹簧相连接,实际热阻很大,所以在建立有限元模型过程中忽略其耦合关系,朝向被加热表面的一侧喷涂黑漆,红外半球发射率为0.9,背向被加热表面的一侧为材料表面自然状态。结合散热面、红外加热笼的轮廓尺寸和黑片热流计的外径,将散热面划分为9个区域,分别建立9个黑片热流计模型,均匀分布在散热面的相应区域,经测量黑片热流计的敏感片及补偿片的红外半球发射率为0.9,并对对应区域的热流计进行编号;热沉温度设置为100 K;散热面喷涂白漆,红外半球发射率和太阳吸收率分别为0.92和0.12。
根据热分析结果得出的散热面高温工况的轨道周期平均外热流,以及散热面参数可以计算出热流计的目标温度:
$$T_{\mathrm{s}}=\sqrt[4]{\frac{\alpha_{\mathrm{s}}}{\varepsilon \sigma}\left(Q_1+Q_2\right)+\frac{Q_3}{\sigma}}$$ (1) 式中:αs为散热面的太阳吸收率;ε为散热面的红外半球发射率;Q1、Q2、Q3分别为太阳直射、地球反照、地球红外的辐射能量;σ为斯特藩-玻尔兹曼常数。
在热试验中,红外加热笼和被加热表面均为平面,并且加热笼的尺寸与被加热表面尺寸基本一致,加热笼与被加热表面之间存在一定距离,这就会出现被加热表面呈现中间受到的热流多,周围受到的热流少的现象[3]。仅采用单一点处的热流计温度值反映出的整个面的热流密度值与实际到达外热流存在误差,因此需要对试验中的误差进行分析。
目前对于均匀程度的定义还不完善,本文采用9个热流计测得的热流密度的均方偏差S来表征到达散热面的均匀程度。
$$S=\sqrt{\frac{\left(t_1-t\right)^2+\left(t_2-t\right)^2+\cdots+\left(t_9-t\right)^2}{9}}$$ (2) 式中:t1~t9为9个敏感片达到平衡时的温度;t为目标热流计温度。
2. 仿真分析
由于红外加热笼主要模拟高温工况下的外热流,这里仅对高温工况的模拟情况进行分析。经热分析得到的轨道周期平均外热流为423.9 W/m2,并以此作为红外加热笼模拟的目标外热流,代入公式(1)计算出高温工况下的目标热流计温度为20.9℃,以此作为控制条件,改变加热笼的加热功率,反复进行稳态计算,当目标热流计达到目标温度且加热笼的占空比为75%时,以此时的功率作为热试验中加热笼需要调节的功率。
9个热流计敏感片达到平衡时的温度如表 1所示,温度云图如图 3所示。
表 1 传统加热笼控制方法下敏感片平衡温度Table 1. Equilibrium temperature of sensitive sheet under conventional heating cage control methodSensitive sheet code Steady state average temperature /℃ RLJ-1 -14.1 RLJ-2 4.5 RLJ-3 -13.1 RLJ-4 2.2 RLJ-5 20.9 RLJ-6 2.7 RLJ-7 18.8 RLJ-8 0.2 RLJ-9 18.9 由表 1、图 3可以看到,在热试验中利用红外加热笼模拟外热流,并使用中心位置的RLJ-5热流计温度作为反馈,来控制红外加热笼的开关,达到平衡时的平均到达热流密度为341.2 W/m2,较目标到达外热流偏低19.6%,不满足总到达能量的保守设计原则。所引起的各点温度偏差最大值可达-34.0℃,9个区域的热流均方偏差为102.0。由图 3的仿真结果可知,该种传统加热笼设计带来的外热流模拟与热分析所得的外热流存在较大误差,由于中心的RLJ-5敏感片温度与周围的温度相比较高,所以利用中心RLJ-5热流计的温度来控制加热笼,会导致模拟的平均外热流明显低于在轨期间的外热流,这会使试验工况无法覆盖在轨期间的极端高温工况。并且加热笼模拟出的外热流的不均匀程度很高,因此若要改善这一现象,需对加热笼进行优化设计。
3. 优化设计
本文针对上述缺点提出改进方法:
① 针对模拟的平均到达热流密度较在轨期间偏低的问题,可以采用周围温度较低的热流计对加热笼进行控制,使得模拟的平均热流密度略大于在轨期间平均热流密度;由于9个热流计具有对称性,分别采用RLJ-1、RLJ-2和RLJ-4热流计的温度作为反馈控制加热笼,进行仿真分析,仿真结果如下,9个热流计敏感片达到平衡时的温度如表 2所示。
表 2 不同热流计控制加热笼时敏感片的平衡温度Table 2. Equilibrium temperature of the sensitive sheet when controlling the heating cage with different heat flow metersSensitive sheet code Steady state average temperature/℃ RLJ-1 to control RLJ-2 to control RLJ-4 to control RLJ-1 19.6 0.4 2.2 RLJ-2 40.7 20.0 21.9 RLJ-3 20.6 1.3 3.1 RLJ-4 37.9 17.5 19.4 RLJ-5 59.2 37.3 39.2 RLJ-6 38.5 18.0 19.9 RLJ-7 14.4 -4.5 -2.7 RLJ-8 35.8 15.5 17.3 RLJ-9 14.5 -4.5 -2.7 经计算,采用RLJ-2和RLJ-4热流计的温度来控制加热笼得到的模拟平均热流密度分别为375.7 W/m2和385.5 W/m2,依然小于在轨期间的平均热流密度,而采用RLJ-1热流计作为控制加热笼的温度数据时,平均热流密度为493.1 W/m2,高于在轨期间的平均热流密度。因此,在试验中可以使用RLJ-1热流计来测量到达外热流并采用其温度数据控制加热笼。
② 针对模拟的平均外热流均匀性差的问题,可以在空间条件允许的条件下,尽可能扩大加热笼尺寸,并使散热面处在加热笼的中心位置,以提高均匀性。另外,扩大加热笼的尺寸势必会导致除被模拟表面的其他表面受到影响,可以采用防辐射屏来遮挡热流。改进后结构示意图如图 4所示。
在原有的加热笼尺寸的基础上,分别将宽和高依次增加30 cm和20 cm,并利用RLJ-1热流计的温度反馈对加热笼进行控制,进行仿真分析,得到的结果如下,9个热流计敏感片达到平衡时的温度如表 3所示。各个尺寸下,平均热流密度曲线如图 5所示,热流密度均方偏差曲线如图 6所示。
表 3 改进结构的敏感片的平衡温度Table 3. Equilibrium temperature of sensitive sheet with improved structure/℃ Sensitive sheet code Heat cage size/mm 233×255 263×275 293×295 323×315 353×335 383×355 413×375 443×395 RLJ-1 19.6 19.4 20.1 20.1 19.7 20 20.2 20.7 RLJ-2 40.7 33 29.5 27.1 24.8 24.3 23.9 24.1 RLJ-3 20.6 20.2 20.9 20.9 20.3 20.8 20.9 21.1 RLJ-4 37.9 32.5 29.7 27.6 25.6 24.8 24.4 23.9 RLJ-5 59.2 45.6 39.1 35.2 31.1 29.2 28 27.3 RLJ-6 38.5 32.8 30 27.9 25.9 25.1 24.6 24 RLJ-7 14.4 16.9 18 18.7 18.6 19.6 20.2 20.4 RLJ-8 35.8 29.4 26.4 24.8 23.1 23.4 23.4 23.6 RLJ-9 14.5 16.8 17.7 18.5 18.4 19.4 20 20.4 由表 3和图 5、图 6可以看出利用RLJ-1热流计温度作为反馈,可以提升模拟的平均外热流至高温工况下的轨道周期平均外热流以上;并且扩大加热笼可以有效改善模拟外热流的均匀程度。但是当尺寸扩大到353 mm×335 mm时,模拟的平均外热流与在轨期间平均外热流的差值为69.2 W/m2,热流密度均方偏差下降至27.0,继续扩大尺寸不能明显减小模拟的平均外热流与在轨期间平均外热流的差值,且无法有效地降低热流密度的均方偏差,也就是说继续扩大尺寸已经无法有效地改善到达散热面的热流均匀程度,由于扩大加热笼的尺寸势必会提升加热笼所需的功率,因此在实际的热试验中可以选用尺寸为353 mm×335 mm的加热笼模拟外热流,既可以满足模拟的平均外热流的均匀性要求,又可以做到空间和功率方面的平衡。
4. 结论
本文根据某微型空间相机在热试验中的外热流模拟装置,采用有限元法建立了系统仿真模型,分析了在传统加热笼控制方法下,到达模拟表面的总能量及热流密度均匀性,并对其结构进行了改进,通过分析改进装置的仿真结果,得出以下结论:
1)对于采用传统加热笼控制方法出现的到达模拟表面的总能量偏低的问题,通过调整热流计粘贴位置,可以使到达模拟表面的热流密度高于在轨期间到达散热面的热流密度,进而满足总到达能量的保守设计原则。
2)对于红外加热笼在该热试验中模拟外热流时出现的热流密度不均匀的问题,通过扩大加热笼能够实现提升均匀度的效果,热流密度的均方偏差从最初的102.0降低至27.0,并且平均热流密度与目标热流密度的差值从69.2 W/m2下降至13 W/m2,从而使试验工况更接近空间相机在轨期间的工作状态。
-
表 1 数学模型参数p计算结果
Table 1 Calculation results of parameter p
No. The X direction displacement (pixel) The Y direction displacement (pixel) The piezo positioner displacement/μm 1 -0.002 -0.004 0 2 0.621 0.639 5 3 1.398 2.206 10 4 2.875 3.327 15 5 3.655 4.573 20 6 5.111 5.878 25 7 6.474 6.993 30 8 7.533 8.109 35 9 8.661 9.221 40 10 9.81 10.321 45 11 10.616 11.579 50 12 12.042 12.871 55 The parameter p of the X direction: 0.2260 pixel/μm,
the parameter p of the Y direction: 0.2364 pixel/μm表 2 某型微扫描器件的测试结果
Table 2 The test results of a micro scanner
The design displacement No. X direction displacement test Y direction displacement test Real value Difference Real value Difference 12.5 μm 1 12.573 0.073 12.089 -0.411 2 12.303 -0.197 12.221 -0.279 3 12.567 0.067 12.121 -0.379 4 12.485 -0.015 12.231 -0.269 5 12.506 0.006 12.075 -0.425 6 12.504 0.004 12.053 -0.447 7 12.576 0.076 12.231 -0.269 8 12.208 -0.292 12.221 -0.279 9 12.507 0.007 12.113 -0.387 10 12.510 0.01 12.113 -0.387 11 12.450 -0.05 12.072 -0.428 12 12.372 -0.128 12.157 -0.343 13 12.451 -0.049 12.174 -0.326 14 12.455 -0.045 12.178 -0.322 15 12.448 -0.052 12.173 -0.327 16 12.498 -0.002 12.149 -0.351 17 12.486 -0.014 12.083 -0.417 18 12.366 -0.134 12.078 -0.422 19 12.454 -0.046 12.158 -0.342 20 12.457 -0.043 12.133 -0.367 Mean 12.4588 -0.0412 12.1412 -0.3589 表 3 本文实验不确定度汇总
Table 3 The summary of the test uncertainty in this paper
Uncertainty component Ui Uncertainty source Uncertainty Sensitivity coefficient ci $ \left| {{c_i}} \right| \cdot {U_i} $ U(m) of the X direction The repeat test of the X direction 0.02014 μm c1=4.425 0.08912 μm U(m) of the Y direction The repeat test of the Y direction 0.01278 μm c1=4.230 0.05406 μm U(m1) of the X direction The error of the piezo positioner at the X direction -0.0041 μm c2=2.791 -0.01144 μm U(m1) of the Y direction The error of the piezo positioner at the Y direction -0.0045 μm c2=2.699 -0.01215 μm The synthetic uncertainty of the X direction: 0.08985 μm; The expand uncertainty of the X direction: 0.17970 μm
The synthetic uncertainty of the Y direction:0.05541 μm; The expand uncertainty of the Y direction: 0.11082 μm -
[1] Bagavathiappan S, Lahiri B B, Saravanan T, et al. Infrared thermography for condition monitoring–a review[J]. Infrared Physics & Technology, 2013, 60: 35-55.
[2] Kogure S, Inoue K, Ohmori T, et al. Infrared imaging of an A549 cultured cell by a vibrational sum-frequency generation detected infrared super resolution microscope[J]. Optics Express, 2010, 18(13): 13402-13406. DOI: 10.1364/OE.18.013402
[3] Lanfrey D B, Trinolet P, Pistone F, et al. New IR detectors with small pixel pitch and high operating temperature[C]//Proc. of SPIE, 2010, 7854: 78540M.
[4] 吕侃, 王世勇. 超分辨率技术在红外微扫描中的应用[J]. 电子设计工程, 2011, 19(13): 166-169. DOI: 10.3969/j.issn.1674-6236.2011.13.050 LV Kan, WANG Shiyong. Application of super-resolution techniques in infrared micro-scanning[J]. Electronic Design Engineering, 2011, 19(13): 166-169. DOI: 10.3969/j.issn.1674-6236.2011.13.050
[5] 张良, 仇振安, 杨小儒, 等. 红外系统微扫描技术研究[J]. 激光与光电子学进展, 2012, 49(4): 1-5. https://www.cnki.com.cn/Article/CJFDTOTAL-JGDJ201204024.htm ZHANG Liang, QIU Zhen'an, YANG Xiaoru, et al. Research of infrared micro-scanning technology[J]. Laser & Optoelectronics Progress, 2012, 49(4): 042302 https://www.cnki.com.cn/Article/CJFDTOTAL-JGDJ201204024.htm
[6] 吴新社, 邓芳轶, 陈敏, 等. 旋转式红外微扫描器研制[J]. 红外与毫米波学报, 2011, 30(3): 263-267. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYH201103017.htm WU Xinshe, DENG Fangyi, CHEN Min, et al. Development of rotary infrared micro-scanner[J]. Journal of Infrared and Millimeter Waves, 2011, 30(3): 263-267. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYH201103017.htm
[7] 王学伟, 李珂, 王世立. 红外成像系统微扫描成像重建算法研究[J]. 光电工程, 2012, 39(12): 122-126. https://www.cnki.com.cn/Article/CJFDTOTAL-GDGC201212023.htm WANG Xuewei, LI Ke, WANG Shili. Microscanning reconstruction algorithm for IR imaging system[J]. Opto-Electronic Engineering, 2012, 39(12): 122-126. https://www.cnki.com.cn/Article/CJFDTOTAL-GDGC201212023.htm
[8] 代少升, 张德州, 崔俊杰, 等. 基于微扫描的红外超分辨率成像系统的设计[J]. 半导体光电, 2017, 38(1): 103-106, 112. https://www.cnki.com.cn/Article/CJFDTOTAL-BDTG201701026.htm DAI Shaosheng, ZHANG Dezhou, CUI Junjie, et al. Design of infrared super-resolution imaging system based on micro-scanning[J]. Semiconductor Optoelectronics, 2017, 38(1): 103-106, 112. https://www.cnki.com.cn/Article/CJFDTOTAL-BDTG201701026.htm
[9] 黄燕, 沈飞, 黄整章, 等. 压电式高精度位移微扫描控制系统设计[J]. 光学精密工程, 2016, 24(10s): 454-460. HUANG Yan, SHEN Fei, HUANG Zhengzhang, et al. Micro-scanning control system design for piezoelectric high-precision displacement[J]. Editorial Office of Optics and Precision Engineering, 2016, 24(10s): 454-460.
[10] 王忆锋, 侯辉, 冯雪艳. 红外焦平面器件微扫描技术的发展[J]. 红外技术, 2013, 35(12): 751-758. http://hwjs.nvir.cn/article/id/hwjs201312002 WANG Yifeng, HOU Hui, FENG Xueyan, Development of microscan techniques in infrared focal plane array[J]. Infrared Technology, 2013, 35(12): 751-758. http://hwjs.nvir.cn/article/id/hwjs201312002
[11] 王林波, 王延杰, 邸男, 等. 基于几何特征的圆形标志点亚像素中心定位[J]. 液晶与显示, 2014, 29(6): 1003-1009. https://www.cnki.com.cn/Article/CJFDTOTAL-YJYS201406024.htm WANG Linbo, WANG Yanjie, DI Nan, et al. Subpixel location of circle target center based on geometric features[J]. Chinese Journal of Liquid Crystals and Displays, 2014, 29(6): 1003-1009. https://www.cnki.com.cn/Article/CJFDTOTAL-YJYS201406024.htm
[12] 梁智滨, 吴鹏飞, 李灵巧, 等. 基于改进Zernike矩和均值漂移的插针位置检测方法[J]. 桂林电子科技大学学报, 2021, 41(4): 305-311. https://www.cnki.com.cn/Article/CJFDTOTAL-GLDZ202104008.htm LIANG Zhibin, WU Pengfei, LI Lingqiao, et al. Pin position detection based on improved Zernike moment and mean shift[J]. Journal of Guilin University of Electronic Technology, 2021, 41(4): 305-311. https://www.cnki.com.cn/Article/CJFDTOTAL-GLDZ202104008.htm
[13] 田光宝, 王见, 王博文. 单目相机非合作目标提取及位姿检测[J]. 红外与激光工程, 2021, 50(12): 20210166. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202112050.htm TIAN Guangbao, WANG Jian, WANG Bowen. Monocular camera non-cooperative target extraction and pose detection[J]. Infrared and Laser Engineering, 2021, 50(12): 20210166. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202112050.htm
[14] Ghosal S, Mecrotra R. Orthogonal moment operator for subpixel edge detection[J]. Pattern Recognition, 1993, 26(2): 295-306.
[15] 卢达, 白静芬, 林繁涛, 等. 基于映射常数的动态量值不确定度评定方法[J]. 电测与仪表, 2022, 59(6): 53-59. https://www.cnki.com.cn/Article/CJFDTOTAL-DCYQ202206008.htm LU Da, BAI Jingfen, LIN Fantao, et al. Evaluation of uncertainty for dynamic values based on mapping constants[J]. Electrical Measurement & Instrumentation, 2022, 59(6): 53-59. https://www.cnki.com.cn/Article/CJFDTOTAL-DCYQ202206008.htm
[16] 中国国家标准化管理委员会. 测量不确定度评定和表示: GB/T 27418-2017[S]. 北京: 中国标准出版社, 2018. Standard Administration. Guide to Evaluation and Expression of Uncertainty in Measurement: GB/T 27418-2017[S]. Beijing: Standards Press of China, 2018.