Design and performance test of self-calibration multichannel infrared radiometer
-
摘要: 为满足遥感器热红外波段自动化观测定标的需求,研制了具有自动化观测能力的自校准多通道红外辐射计。该设备需具有以下特色功能:1)采用电机驱动镀金反射镜的设计,实现0°~90°仰角的大气下行辐射和地表辐亮度的测量,为消除大气下行辐射对反演地表温度的影响提供了技术手段;2)采用滤光轮分光的方法实现了6个光谱通道的自动设置,结合多通道温度与发射率分离算法可以实现场地温度与发射率的分离,为卫星遥感器热红外波段绝对辐射定标提供了2个关键因子;3)采用内置2个控温精度分别优于0.04 K和0.05 K、发射率均高于0.994,稳定性均优于0.0014的黑体实现内部探测器的实时辐射定标,有效地消除了内部背景辐射对辐射测量的影响,定标不确定度小于0.167%。等效测温不确定度为0.2 K(@303 K, 11 μm),为遥感器热红外波段场地自动化定标应用的开展奠定了基础。Abstract: In order to meet the needs of automated observation calibration for remote sensor thermal infrared band, a self-calibration multichannel infrared radiometer with automated observation capability was developed. The equipment has the following features: 1) the design of motor-driven gold-plated reflector was used to achieve the measurement of downward atmospheric radiation and surface radiance at an elevation angle of 0°~90°, which provided the technical means to eliminate the influence of downward atmospheric radiation on inversion surface temperature; 2) the filter wheel spectrometry was used to realize the automatic setting of 6 spectral channels, and combined with the multichannel temperature and emissivity separation algorithm, the separation of site temperature and emissivity could be achieved, which provided 2 key factors for the absolute radiometric calibration of satellite remote sensor thermal infrared band; 3) a black body with two built-in temperature control accuracy better than 0.04 K and 0.05 K, emissivity higher than 0.994, stability better than 0.001 4 was used to achieve the real-time radiometric calibration of the internal detector, which effectively eliminated the influence of internal background radiation on radiation measurement. The calibration uncertainty is less than 0.167%, and the equivalent temperature measurement uncertainty is 0.2 K (@303 K, 11 μm), which lays the foundation for the application of site automated calibration for remote sensor thermal infrared band.
-
Keywords:
- calibration /
- automation /
- radiometer /
- self-calibration /
- uncertainty analysis
-
引言
机场跑道外来物(FOD,foreign object debris)泛指可能危害航空器或系统的某种外来物质,如飞机和发动机连接件、机械工具、飞行物等[1]。机场跑道外来物一直以来都是危害飞机飞行安全的重要因素之一,通过在跑道周边部署FOD光电探测系统,对跑道路面不同距离、不同大小外来物进行实时监控,快速精确地检测与告警成为保障飞机飞行安全的重要措施。作为光电探测系统的“眼睛”,变焦镜头起着至关重要的作用。变焦镜头是指在像面保持不变的前提下,通过改变系统2个或2个以上透镜组之间间隔实现焦距连续变化,达到小倍率、大视场发现物体,大倍率、小视场识别物体的目的。针对机场跑道外来物光电探测系统远距离(≤500 m)识别微小物体(尺寸≥20 mm×20 mm)的特殊需要,研制高像质、长焦距、大变倍比以及小型化的变焦镜头显得尤为重要。
2007年,华中光电研究所的胡际先[2]设计了一款75 mm~900 mm连续变焦系统,系统总长小于1 000 mm,通光口径140 mm,该光学系统结构合理,各透镜的加工工艺性能良好,并通过了像质测试、环境试验及实景成像试验。2009年,胡际先等人[3]又设计了一款连续透雾变焦系统,该系统焦距为20 mm~400 mm,系统总长不大于280 mm,实现了普通彩色和透雾黑白两种模式的观察,保证了系统的齐焦性,并针对系统小型化做了详细的研究。2015年,姚清华等人[4]设计了一款焦距为10 mm~500 mm的大变倍比长焦距变焦镜头,该系统采用2个新型非球面和一组双层衍射元件,有效地校正了光学系统的像差。2018年,白虎冰等人[5]通过Zemax软件仿真设计了一款变焦范围为100 mm~600 mm的光学变焦系统,整个变焦系统采用球面镜片,口径为135.56 mm,总长为563.956 mm,实现了远距离成像的应用要求。随着变焦镜头焦距以及变倍比的不断增大,系统的色差校正成为一个亟需解决的难题[6-7]。为了获得高像质、大变倍、长焦距的变焦系统,本文采用机械补偿方式实现连续变焦,通过合理分配光焦度,引入多片超低色散(extra-low dispersion,ED)镜片,利用其独特的超低色散特性,有效地校正了系统色差和其他高级像差[8],并通过Zemax光学设计软件进行了像质优化,提高了成像质量。
1 系统参数及设计原理
1.1 设计参数
根据中国民用航空局机场司《机场道面外来物探测设备》信息通告要求,结合现有成熟技术,同时为了适应不同照度,拟选用200万像素Sony逐行扫描1/1.8″CMOS芯片,靶面尺寸为7.76 mm×4.36 mm,对角线长度为9.0 mm,摄像机像元尺寸为3.45 μm。根据该类摄像机像元尺寸,选用镜头焦距满足如下公式:
$$ \frac{f}{R}=\frac{(n\times d)R}{a} $$ (1) 式中:
$ f $ 为镜头焦距;$ R $ 为作用距离,取500 m;$ d $ 为摄像机像元尺寸,取3.45 μm;$ n $ 为探测目标所需像元个数;$ a $ 为目标尺寸,取20 mm。视场满足如下公式:$${\rm{\theta }} = 2{\rm{arctan}}\frac{l}{{2f}}$$ (2) 式中:
$ \theta $ 为视场角;$ l $ 为相机像面对角线长度;$ f $ 为镜头焦距。相应技术参数如表1所示。表 1 光学设计参数Table 1. Optical design parameters参数 技术指标 焦距范围/mm 50~1 000 系统总长/mm ≤400 工作波长/nm 486~656 变倍比 20× F数 6~10 视场角/(°) 10.2~0.51 MTF(@100 lp/mm) >0.2 弥散斑半径(RMS)/μm <10 1.2 设计过程
变焦镜头根据其补偿方式不同,主要分为光学补偿变焦镜头、机械补偿变焦镜头、双组联动变焦光学镜头和全动型变焦光学镜头。其中,机械补偿变焦镜头因为其结构简单、成像质量好且容易实现较大的变倍比与较小的像面漂移被广泛采用[9]。机械补偿变焦光学镜头主要由前固定组φ1、变倍组φ2、补偿组φ3和后固定组φ4组成,如图1所示。变倍组φ2与补偿组φ3各自按照不同的运动曲线移动,从而实现焦距连续变化。目前,变焦系统的结构设计方法主要有两种:一种是近轴高斯光学理论计算法,此法计算复杂,错误率高,容易产生无效解,一般应用于结构较为简单的变焦系统中[5];另外一种是光学设计软件(Zemax、Code V、OSLO等)辅助计算法,此法通过查找专利、文献进行等比缩放,再配合评价函数操作数进行仿真优化得到最终结果,一般用于较为复杂的结构设计中。
本文采用机械变焦结构,通过查阅相关文献以及相关专利[9-11]确认初始结构。根据设计要求,对初始结构各组分进行归一化处理,合理分配各组分的焦距和间隔,然后通过软件中相关操作数进行优化设计。在大变倍、长焦距变焦镜头优化设计过程中,色差校正是需要着重考虑的环节之一。根据消色差理论,整个系统应满足:
$$\nabla \delta _k{'} = - \frac{1}{{n_k{'}\mu _k'^{2}}}{{\Sigma }}{h^2}\frac{\phi }{\upsilon } = 0$$ (3) 式中:
$ {n}_{k}' $ 为第$ k $ 个折射面的折射率;$\; {\mu }_{k}{'} $ 为出射光线与光轴的夹角;$ h $ 为光线的入射角度。根据公式(3)可知,系统色差不仅与透镜的阿贝数$ \upsilon $ 和光焦度$ \phi $ 有关,还与光线在透镜表面的入射高度有关,即要控制透镜在光路中的位置,因此,校正色差主要从材料和结构调整来实现[12-15]。超低色散镜片具有高折射率、高阿贝尔系数,其独特的色散特性可以有效解决系统色散问题。而双胶合镜组不仅可以减少系统镜组数量,简化整体结构,而且有利于校正球差和色差。本文通过采用双胶合结构与超低色散镜片搭配来实现系统色差校正[8, 12]。在长焦状态,前固定组是色差的主要来源,通过采用单片ED镜片FCD100以及胶合镜组(TAFD5G+FCD100)搭配来进行色差优化。其中,选用阿贝尔系数相差较大的两种材料TAFD5G(Nd:1.8348,Vd:42.721)与FCD100(Nd:1.437,Vd:95.1)来进行胶合,不仅减小了前固定组的色差贡献,又避免了镜组光焦度过大导致入射角偏折而增加的高级球差[13]。与此同时,为了平衡变焦状态的色差分布,在补偿组中引入ED镜片FCD100进行短焦状态色差校正。最后通过优化操作数加以约束,调整各个透镜的曲率以及镜组间隔来改变光线的偏折角,从而减小系统的色差贡献。
2 设计结果及像质评价
2.1 设计结果
本文采用机械变焦结构设计了一款大变倍比、长焦距、高像质的变焦镜头。该系统对可见光波段成像,焦距变化范围50 mm~1 000 mm,20倍变焦,系统结构如图2所示。图2从上到下分别为短焦(50 mm)、中焦(300 mm、500 mm)、长焦(1 000 mm)4种状态下外形结构图。由图2可知,系统共有28片透镜组成,总长399.9 mm,最大口径为100 mm,后截距22 mm。其中,前固定组由5片球面透镜组成,为正光焦度,其主要作用是为系统提供相应的工作距以及平衡像差,减轻后固定组校正压力。变倍组由6片球面透镜组成,为负光焦度,通过移动位置来实现变焦系统倍率的改变。补偿组由4片球面透镜组成,为正光焦度,按照一定的曲线规律做非线性运动来补偿变倍组移动所引起的像面变化。后固定组由13片球面透镜组成,为正光焦度,主要采用由冕牌玻璃与火石玻璃搭配的双胶合以及三胶合透镜的结构形式来进行残余像差的矫正。
2.2 像质评价
2.2.1 传递函数
调制传递函数(modulation transfer function,MTF)是所有光学系统性能判据中最全面的判据,特别是对于成像系统而言,它反映了光学系统对物体不同频率成分的传递能力。利用Zemax设计软件可以得到变焦镜头在短焦(50 mm)、中焦(300 mm、500 mm)、长焦(1 000 mm)4种状态下MTF曲线,如图3所示。由图3可知,系统在各个焦距位置均呈现良好的成像质量,在20 lp/mm空间频率时MTF>0.8,50 lp/mm空间频率时MTF>0.5,100 lp/mm空间频率下MTF>0.2,均逼近衍射极限。由此可知,物体在成像中信息传递损失少,有利于机场跑道微小外来物的探测识别。
2.2.2 点列图
弥散斑半径是一个重要的成像质量参数,其斑点半径越小,系统成像质量越好,各组态下弥散斑直径如表2所示。从表2中数据可知,整个光学系统的RMS弥散斑半径在中心视场下各个焦距位置处均小于4 μm,整个光学系统的GEO弥散斑半径在中心视场下各个位置处均小于7 μm,可以看出整个系统在各焦距位置都有比较好的成像质量。
表 2 各组态下弥散斑半径Table 2. Dispersed spots diameter of each configuration焦距/mm 50 300 500 1 000 RMS半径/μm 1.013 1.897 3.324 3.050 GEO半径/μm 1.960 3.345 5,570 6.951 2.2.3 畸变、场曲图
系统畸变、场曲曲线如图4所示。图4(a)~4(d)分别表示系统在短焦(50 mm)、中焦(300 mm、500 mm)、长焦(1 000 mm)4种状态下的畸变、场曲情况。由图4可看出,当系统焦距为50 mm时,系统处于广角状态,全视场为20°,畸变相对较大,接近1.0%;当系统处于中焦和长焦位置时,系统畸变均小于0.1%。与此同时,通过场曲曲线我们发现,无论系统处于短焦、中焦还是长焦状态下,系统场曲均小于0.1 mm,均满足系统对变焦镜头成像的要求。
2.2.4 色差曲线
系统在长焦(1 000 mm)状态下垂轴色差以及轴向像差如图5所示。由图5可看出,由F(0.65 mm)光、C(0.48 mm)光在轴上差异可知,系统在长焦状态下位置色差均在0.2 mm范围内,垂轴色差全视场均在艾里斑内,且垂轴色差最大值在2.5 μm内,保证了长焦状态下的成像质量,满足系统对变焦镜头成像的要求。
2.2.5 变焦补偿曲线
机械补偿的变焦系统需要绘制凸轮曲线,以此来检查变倍组和补偿组运动过程的合理性[10]。通过基于光线追迹法的Zemax宏命令[11],计算出系统变焦过程中变倍组与补偿组的位移量,运用MATLAB软件进行编程,代入初始参数,求得变焦系统变倍组与补偿组的运动变化曲线,如图6所示。图6中横坐标表示凸轮转动的角度,纵坐标表示随凸轮转角变化时变倍组和补偿组相对于像面的距离,其中光阑到像面的距离为零。变倍组随凸轮转角作线性运动,使系统焦距发生变化,补偿组作非线性运动,补偿变焦过程中像面移动,变倍组和补偿组的导程分别为55.48 mm和59.69 mm,变焦曲线平滑无断点,有利于凸轮加工。
3 结论
根据中国民用航空局机场司《机场道面外来物探测设备》信息通告,通过选取合理的初始结构,计算光学参数,采用Zemax光学设计软件对变焦系统进行像质优化,在设计过程中利用超低色散镜片解决了长焦镜头色差难以纠正的问题,通过优化镜片间距以及镜片厚度较好地解决了现有国产大变倍、长焦距变焦镜头无法小型化的问题。经过优化设计的变焦镜头具有较好的成像质量,MTF接近衍射极限,此外,通过研究变焦理论与变焦运动方程,运用MATLAB软件绘制了变倍组与补偿组的运动曲线,证明了变焦运动过程的合理性。该设计也为下一步机场跑道外来物光电探测系统的研发设计提供了铺垫,为其他变焦镜头设计提供了依据。
-
表 1 光学系统结构参数
Table 1 Optical system structural parameters
性能 参数 视场光阑直径/mm 18 孔径光阑直径/mm 26 光阑间距/mm 120 镀金反射镜尺寸/mm 40 视场大小/° 4 表 2 CES100-06黑体性能参数
Table 2 CES100-06 black body performance parameters
性能 参数 型号 CES100-06 辐射面积/mm2 250×250 温度范围/℃ −50~100 控温精度/℃ ±0.03,0~100
±0.05,−50~0发射率 0.992±0.001 控制箱尺寸/mm 133.35×482.6×457.2 电源 100~240VAC±10%(50/60Hz) 表 4 MSIR各通道参数与定标系数
Table 4 Parameters and calibration coefficients of each channel
波段范围/μm ai bi ni di Si Ii 8~14 452.483 76 953.390 08 0.924 34 2.436 87 461.558 08 2 423.595 00 8.1~8.3 57.835 96 1 704.534 51 0.998 24 0.965 94 777.557 77 17 218.022 13 8.3~8.5 47.156 27 1 639.871 17 0.996 66 0.918 54 800.284 12 17 234.965 18 8.9~9.3 55.478 94 1 563.678 19 0.999 63 1.011 89 618.901 69 16 824.752 65 10.3~11.3 53.915 70 1 340.639 07 1.000 67 1.034 66 514.255 59 15 363.374 07 11.5~12.5 35.209 79 1 192.883 98 1.002 36 1.036 04 487.120 23 15 180.425 98 表 5 辐射定标不确定度分析
Table 5 Uncertainty analysis of radiometric calibratio
不确定源 符号 不确定度类型 数值/% 等效温度不确定度/K(@303.15 K,11 μm) 定标黑体发射率不确定度 u1 B 0.140 - 控温不确定度 u2 B 0.072 - 定标系数拟合不确定度 u3 A 0.019 - 背景噪声引入的不确定度 u4 A 0.052 - 合成不确定度 U - 0.167 0.2(k=2) 表 3 内置黑体发射率定标结果
Table 3 Calibration results of built-in black body emissivity
次数 环温黑体 加热黑体 拟合斜率 发射率 拟合斜率 发射率 1 −0.001 2 0.993 19 −0.001 9 0.993 9 2 −0.001 3 0.993 27 −0.002 0 0.994 3 −0.002 1 0.994 07 −0.002 5 0.994 5 4 −0.002 0 0.993 91 −0.002 2 0.994 2 5 −0.001 7 0.993 72 −0.001 8 0.993 8 6 −0.001 8 0.993 76 −0.001 9 0.993 9 7 −0.001 7 0.993 68 −0.002 2 0.994 2 8 −0.001 1 0.993 08 −0.002 4 0.994 4 9 −0.002 1 0.994 13 −0.002 1 0.994 1 10 −0.001 1 0.993 12 −0.002 3 0.994 3 平均值 0.993 6 - 0.994 1 标准偏差 - 0.04% - 0.023% -
[1] 张勇, 祁广利, 戎志国, 等. 卫星红外遥感器辐射定标模型与方法[M]. 北京: 科学出版社, 2015: 21-22. ZHANG Yong, QI Guangli , RONG Zhiguo, et al. The model and method of radiometric calibration for satellite infrared remote sensor[M]. Beijing: Science Press, 2015: 21-22.
[2] 邱刚刚, 李新, 韦玮, 等. 基于场地自动化观测技术的遥感器在轨辐射定标试验与分析[J]. 光学学报,2016,36(7):0701001. doi: 10.3788/AOS201636.0701001 QIU Ganggang, LI Xin, WEI Wei, et al. Experiment and analysis of on-orbit radiometric calibration for remote sensors based on in-site automated observation technology[J]. Acta Optica Sinica,2016,36(7):0701001. doi: 10.3788/AOS201636.0701001
[3] 吕佳彦, 何明元, 陈林, 等. 基于敦煌辐射校正场的自动化辐射定标方法[J]. 光学学报,2017,37(8):0801003. doi: 10.3788/AOS201737.0801003 LYU Jiayan, HE Mingyuan, CHEN Lin, et al. Automated radiometric calibration method based on Dunhuang radiometric calibration site[J]. Acta Optica Sinica,2017,37(8):0801003. doi: 10.3788/AOS201737.0801003
[4] 张孟, 韦玮, 张艳娜, 等. 基于场地自动化观测技术的SNPP VIIRS高频次在轨辐射定标[J]. 光学学报,2019,48(4):0428001. ZHANG Meng, WEI Wei, ZHANG Yanna, et al. High-frequency on-orbit radiometric calibration of SNPP VIIRS based on in-site automated observation technology[J]. Acta Optica Sinica,2019,48(4):0428001.
[5] TSUCHIDA R, WELCH E, ZALEWSKI, et al. ASTER preflight and inflightcalibration and the validation of level 2 products[J]. IEEE,1998,36(4):1161-1172.
[6] TONOOKA H, PALLUCONI F, HOOK S, et al. Vicarious calibration of ASTER thermal infrared bands[J]. IEEE,2005,43:2733-2746.
[7] SCHOTTR. Thermal infrared radiometric calibration of the entire landsat 4, 5, and 7 archive (1982–2010)[J]. Remote Sens. Environ.,2012,122:41-49. doi: 10.1016/j.rse.2011.07.022
[8] 童进军, 邱康睦, 李小文, 等. 一种卫星遥感仪器热红外通道在轨绝对辐射定标新方法[J]. 红外与毫米波学报,2005,24(4):277-280. doi: 10.3321/j.issn:1001-9014.2005.04.009 TONG Jinjun, QIU Kangmu, LI Xiaowen, et al. New method of in-flight absolute calibration for thermal infrared channals of satellite sensor[J]. Infrared millim waves,2005,24(4):277-280. doi: 10.3321/j.issn:1001-9014.2005.04.009
[9] 张允祥, 李新, 韦玮, 等. 基于多通道温度与发射率分离算法的敦煌场地红外特性研究[J]. 光学学报,2019,39(10):1028003. doi: 10.3788/AOS201939.1028003 ZHANG Yunxiang, LI Xin, WEI Wei, et al. Study on infrared characteristics of dunhuang site based on multi-channel temperature and emissivity separation algorithm[J]. Acta Optica Sinica,2019,39(10):1028003. doi: 10.3788/AOS201939.1028003
[10] BOREL C C.Iterative of surface emissivity and temperature for a hyperspectral sensor[C]// First JPL workshop on remote sensing of land surface emissivity, Pasadena:Jet Publication Laboratory,1997:1-5.
[11] 张勇, 祁广利, 戎志国, 等. 卫星红外遥感器辐射定标模型与方法[M]. 北京: 科学出版社, 2015: 21-22. ZHANG Yong, QI Guangli, RONG Zhiguo, et al. The model and method of radiometric calibration for satellite infrared remote sensor[M]. Beijing: Science Press, 2015: 21-22.
[12] 田国良, 柳钦火, 陈良富, 等. 热红外遥感[M]. 第2版. 北京: 电子工业出版社, 2014: 208-213. TIAN Guoliang, LIU Qinhuo, CHEN Liangfu, et al. Thermal remote sensing[M]. 2nd ed. Beijing: Electronic Industry Press, 2014: 208-213.
[13] GILLESPIE A R. Lithologic mapping of silicate rocks using TIMS[C]//Kahle A B,Abbott E,eds,The TIMS data user’s workshop,pasadena:Jet Publication Laboratory,1985:29-44.
[14] WATSON K. Spectral ratio method for measuring emissivity[J]. Remote Sens. Environ.,1992,42:113-116. doi: 10.1016/0034-4257(92)90094-Z
[15] KERR P S, HOOK S J. Separating temperature and emissivity in thermal infrared multispectral scanner data: Implications for recovering land surface temperature[J]. Geoscience and Remote sensing,1993,31:1155-1164. doi: 10.1109/36.317447
[16] SALISBURY J W, DARIA D M. Thermal-infrared remote sensing and Kirchhoff’s law: 1: laboratory measurements[J]. Journal of GeophysicalResearch,1994,99:11897-11911.
[17] GILLESPIE A R, MATSUNAGA T, ROKUGAWA S, et al. A temperature and emissivity separation from advanced space-borne thermal emission and reflection radiometer(ASTER) images[J]. IEEE Transactions on Geoscience and Remote Sensing, 1998, 36(4): 1113-1126.
[18] 邱文龙. 黑体辐射源校准技术研究[D]. 西安: 西安工业大学, 2007. QIU Wenlong. Research on Calibration Technology of Blackbody Radiation Source[D]. Xi’An: Xi’An Industrial University, 2007.
-
期刊类型引用(10)
1. 胡长德,海书亮,徐小琴,王晨斌,马盛中. 一种调光调焦变倍训练系统的研究. 光学技术. 2025(01): 13-18 . 百度学术
2. 杨玉霞,王艳云. 基于单运动组元的双视场像方远心显微镜设计. 光学技术. 2025(02): 183-188 . 百度学术
3. 王起,张国芳. 基于同轴柯勒照明的大变倍比紧凑型显微镜设计. 中国光学(中英文). 2024(04): 921-931 . 百度学术
4. 吴羽婷,林志强,王敏. 15 mm~300 mm宽光谱变焦光学系统设计. 应用光学. 2023(03): 491-499 . 本站查看
5. 金鑫,肖作江. 大光斑激光眩目光学系统设计. 长春理工大学学报(自然科学版). 2023(04): 56-60 . 百度学术
6. 张坤,曲正,钟兴,李京宸,张茜. 30倍轻小型变焦光学系统设计. 光学精密工程. 2022(11): 1263-1271 . 百度学术
7. 陈鑫,甄胜来,方健,陶善静,俞本立. 用于偏振成像系统的20 mm~200 mm变焦物镜设计. 应用光学. 2022(05): 846-852 . 本站查看
8. 宋鹏飞,王春艳,孙昊,王志强,刘欢. 10倍变焦安防镜头光学系统设计. 长春理工大学学报(自然科学版). 2022(04): 39-43 . 百度学术
9. 王晓艳,徐高魁. 高隔离度激光通信终端光学系统设计. 红外与激光工程. 2021(07): 228-232 . 百度学术
10. 赵迎,接长伟. 激光通信终端光学系统杂散光抑制技术研究. 激光与红外. 2021(10): 1378-1382 . 百度学术
其他类型引用(5)