低温环境下红外辐射标定技术研究.pdf

第6卷 第1期2 0 2 3年1月飞控与探测F l i g h tC o n t r o l&D e t e c t i o nV o l.6N o.1J a n u a r y2 0 2 3低温环境下红外辐射标定技术研究*赵海涛1,李士刚2(1.海装上海局驻南京地区第一军事代表室南京2 1 0 0 0 0;2.海装驻上海地区第六军事代表室上海2 0 1 1 0 9)摘 要:低温环境下,红外场景生成装置的标定是对低温红外目标特征评估和精准探测的前提。通过标定过程中热辐射-光子-电子转换测量的传递途径,建立基于红外热像仪灰度的低温标定模型,采用最小二乘法进行拟合,得到模型中的辐射出射度响应函数和系统固有偏置。搭建了低温标定实验装置,根据热像仪对低温辐射源不同辐射出射度的测量结果,得到标定方程。分析了理论值和实际测量值的误差,在低温黑体辐射出射度9.7 9W/m2处,热像仪的辐射出射度误差最大,为-0.1 7W/m2,此时的灰度误差为-9.9 1 D N。关键词:低温;标定;红外场景生成装置;误差分析 中图分类号:T N 2 1 5文献标志码:A文章编号:2 0 9 6-5 9 7 4(2 0 2 3)0 1-0 0 7 5-0 5R e s e a r c ho nI n f r a r e dR a d i a t i o nC a l i b r a t i o nT e c h n o l o g y i nL o wT e m p e r a t u r eE n v i r o n m e n tZ HAO H a i t a o1,L IS h i g a n g2(1.T h eF i r s tM i l i t a r yR e p r e s e n t a t i v eO f f i c eo f t h eN a v a lE q u i p m e n tS h a n g h a iB u r e a u i nN a n j i n g,N a n j i n g2 1 0 0 0 0;2.S i x t hM i l i t a r yR e p r e s e n t a t i v eO f f i c eo f t h eN a v a lE q u i p m e n tD e p a r t m e n t i nS h a n g h a i,S h a n g h a i 2 0 1 1 0 9)A b s t r a c t:T h ec a l i b r a t i o no f i n f r a r e ds c e n eg e n e r a t i o nd e v i c e s i n l o wt e m p e r a t u r ee n v i r o n m e n t i s t h ep r e m i s eo f l o wt e m p e r a t u r e i n f r a r e dt a r g e t f e a t u r ee v a l u a t i o na n da c c u r a t ed e t e c t i o n.T h r o u g ht h e t r a n s m i s s i o np a t ho f t h e r m a l r a-d i a t i o n-p h o t o n-e l e c t r o nc o n v e r s i o n m e a s u r e m e n ti nt h ec a l i b r a t i o np r o c e s s,al o w-t e m p e r a t u r ec a l i b r a t i o n m o d e lb a s e do nt h eg r a y s c a l eo f t h e i n f r a r e dt h e r m a l i m a g e r i se s t a b l i s h e d,a n dt h e l e a s t s q u a r e sm e t h o di su s e dt of i t i t,a n dt h er a d i a t i o ne j e c t i o nr e s p o n s ef u n c t i o na n dt h e i n h e r e n tb i a so f t h es y s t e mi nt h em o d e l a r eo b t a i n e d.Al o w-t e m p e r a t u r ec a l i b r a t i o ne x p e r i m e n t a l d e v i c e i s s e t u p,a n d t h e c a l i b r a t i o ne q u a t i o n i so b t a i n e da c c o r d i n g t o t h em e a s-u r e m e n t r e s u l t so fd i f f e r e n t r a d i a t i o ne m i s s i o nd e g r e e so f l o w-t e m p e r a t u r er a d i a t i o ns o u r c e sb yt h e t h e r m a l i m a g e r.T h ee r r o rb e t w e e nt h et h e o r e t i c a lv a l u ea n dt h ea c t u a l m e a s u r e dv a l u ei sa n a l y z e d.A tt h el o w t e m p e r a t u r eb l a c k b o d yr a d i a t i o ne x i td e g r e eo f9.7 9W/m2,t h er a d i a t i o ne m i s s i o ne r r o ro ft h et h e r m a l i m a g e ri st h el a r g e s t(-0.1 7W/m2),a n dt h eg r a y s c a l ee r r o ra t t h i s t i m e i s-9.9 1 D N.K e y w o r d s:l o wt e m p e r a t u r e;c a l i b r a t i o n;i n f r a r e ds c e n eg e n e r a t i o nd e v i c e;e r r o ra n a l y s i s0 引 言随着红外成像探测系统在卫星、深空探测器、空间反卫星武器中的广泛应用,对低温环境下红外高精度标定的需求与日俱增。在深空冷背景下目标红外特性的研究中,为了模拟深空冷背景下的红外目标特性,通常将红外场景生成装置放置在冷舱所产生的低温真空环境中1-2。对于低温环境下红外成像制导半实物仿真系统,红外辐射标定技术是探测器对低温红外目标特征评估和精准探测的关键。目前的红外辐射标定设备主要有红外热像仪、红外光谱辐射计和红外辐射计三种。在测量低温目标时,红外热像仪最大的不足是其测温下限通*收稿日期:2 0 2 2-1 0-2 5;修回日期:2 0 2 3-0 1-0 7作者简介:赵海涛(1 9 7 8),男,硕士,高级工程师。E-m a i l:z h t_0 0 0 01 6 3.c o m常都在2 5 3 K以上,无法对低于测温下限的目标温度进行直接测量3-4。红外光谱辐射计可以对辐射源固定温度在不同波段范围进行标定,但其结构复杂,大量杂散辐射产生的噪声使低温微弱目标信号的提取十分困难5-6。红外辐射计和红外光谱辐射计都不适用于目标温度分布不均匀的面源辐射出射度测量7-8。由此可见,对于低温红外辐射场景的标定,现有的红外辐射测量设备和方法已经不能满足低温环境下红外场景生成装置的标定要求。本文基于低温环境下热像仪对红外辐射的热辐射-光子-电子转换测量的传递途径,采用红外热像仪灰度的低温标定模型和最小二乘拟合方法,在低于测温极限的温度范围内进行标定,获得热像仪2 5 3 K以下输出的灰度和探测到的辐射出射度之间的关系,得到低于热像仪测温极限的灰度-辐射出射度标定曲线。搭建低温标定实验装置,分析理论值和实际测量值的误差。该方法具有实现简单、误差小的优点,可用于低温环境下红外辐射的标定。1 标定原理在低温环境下,红外辐射标定过程中使用的辐射源为标准黑体,探测器为红外热像仪。黑体发出的辐射经过热像仪光学系统后照射到热像仪焦平面上。标定过程中,黑体的热辐射-光子-电子转换测量传递过程如图1所示。图1 黑体辐射源、传输介质和热像仪传递关系示意图F i g.1 B l a c k b o d yr a d i a t i o ns o u r c e,t r a n s m i s s i o nm e d i u ma n dt h e r m a l i m a g i n gc a m e r a t r a n s m i s s i o nr e l a t i o n s h i ps c h e m a t i c黑体发射率为b b,单元面积dA0的出射光通量聚焦在热像仪光学系统的入瞳处,热像仪光学系统单位面积为dA1,2个面的法向矢量和视线间的角度分别为0和1。传输过程中,o p t为各光学元件的透过率乘积。光学系统入瞳孔径处的光通量被转换为电信号,由放大器进行放大,再经过一系列信号处理和数模转换过程,输出数字图像灰度值。热像仪红外焦平面阵列上单个像元所对应的辐射立体角为9-1 0=Adc o s0/(f/c o s0)2(1)式中,Ad为红外焦平面阵列每个像元的面积;f为光学系统的焦距。辐射源辐射到焦平面阵列单个像元的辐射通量可以表示为1 1=o p tD22c o s0 b bM(Tb b)=o p tAdb b(c o s0)44F2M(Tb b)(2)式中,D为红外镜头的直径;F为红外镜头相对孔径(D/f)的倒数;M(Tb b)为黑体温度为Tb b时的辐射出射度。假设热像仪正对黑体进行探测,则0较小,可认为(c o s0)41(3)在单位时间内,热像仪由于接收到黑体辐射通量而产生的光子数量可以表示为ND e t=h c(4)式中,为波长;h为普朗克常量;c为真空中的光速。在热像仪积分时间t内,由光子转化为电子的数目为S=th c(5)式中,为热像仪红外焦平面阵列的光子量子效率。定义辐射标定系数z为积分时间t内,热像仪产生的电子数与输出的灰度值Gd e t之间的比例系数,灰度值单位为D N,B为热像仪固有偏置。Gd e t=z S+B(6)热像仪输出的灰度值和热像仪接收到的黑体辐射功率密度MB B的关系可以表示为Gd e t=RMB B+B(7)式中,R=o p tAdb bt/4F2h c,R定义为热像仪积分时间内,输出的灰度值对入射辐射出射度的响应。对于同一探测系统而言,R为常数,热像仪输出的灰度值和接收到的辐射出射度呈线性关系。在相同积分时间下,测量黑体在m个辐射出射度情况下的灰度数据,到m个关于热像仪的响应函67飞控与探测第6卷第1期赵海涛,等:低温环境下红外辐射标定技术研究数和系统偏置的方程。用最小二乘法对各辐射出射度对应的灰度数据进行线性拟合1 2-1 3,计算出系统的辐射出射度响应函数R和系统固有偏置B。红外热像仪实际测量的灰度数据Gi和拟合计算得到理论上的灰度数值G i之间的差值的平方和可以表示为2=m i nmi=1(G i-Gi)2,(i=1,2,m)(8)将2分别对辐射出射度响应函数R和系统固有偏置B求偏导数,并令偏导数最小,可以得到最终2个函数的拟合结果。2 实验装置在标定热像仪之前,需要进行如下准备工作:调整热像仪的方位,使热像仪光学系统光轴与黑体中心的法线重合。调整黑体的位置,使黑体的图像充满热像仪视场。将热像仪的测量发射率参数设为黑体的发射率。调节热像仪的焦距,使黑体能够清晰成像,在之后的标定过程中,保持热像仪焦距不变。采用表面涂有高发射率涂层的不锈钢腔体作为低温黑体辐射源,腔体通入液氮制冷,低温黑体表面发射率为0.9 2,温度探测设备的工作波段为7.3 9m,在此范围内,根据黑体普朗克定律,得到黑体温度和辐射 出 射 度 的 关 系 如 图2所 示。当 黑 体 温 度 为1 0 0.0 K时,辐射出射度为0.0 0 0 5 W/m2,当黑体温度为2 5 0.0 K时,辐射出射度为1 3.8 5 W/m2。图2 低温黑体温度-辐射出射度关系F i g.2 R e l a t i o n s h i po f l o w-t e m p e r a t u r eb l a c k b o d yt e m p e r a t u r ea n dr a d i a t i o ne m i s s i o n低温黑体的温度由表面的测温铂电阻进行接触式测温取平均值得到,用热像仪对临近测温铂电阻的标定区域进行测量,标定区域内所有像素点的灰度平均值设为热像仪输出的灰度参数。3 标定结果及误差3.1 标定结果在红外热像仪积分时间为4 0s的情况下,在辐射源温度1 0 0.02 5 0.0 K范围内,对入射的辐射出射度数据MB B和热像仪输出的灰度Gd e t进行标定,标定结果如图3所示。拟合结果显示,红外探测系统的响应函数R为5 9.7 9(D Nm2)/W,系统的固有偏置B为2 3 0.9 7 D N,红外热像仪的灰度标定曲线为Gd e t=5 9.7 9MB B+2 3 0.9 7(9)图3 红外热像仪灰度-标准黑体辐射出射度标定曲线F i g.3 C a l i b r a t i o nc u r v eb e t w e e n i n f r a r e dt h e r m a l i m a g i n gc a m e r ag r a y s c a l ea n ds t a n d a r db l a c k b o d yr a d i a t i o ne m i s s i o n3.2 标定误差当辐射源温度在1 0 0.02 5 0.0 K的低温范围内时,无法通过热像仪对辐射源的温度进行直接测量,因此需要利用3.1节中的热像仪的灰度-辐射出射度标定曲线,通过热像仪的灰度来表征辐射源的辐射出射度。热像仪测量得到的灰度值Gd e t经过灰度-辐射出射度标定曲线计算,得到低温黑体辐射出射度的测量值Md e t。根据此时低温黑体的辐射出射度理论值MB B,得到的实验结果如图4所示。MB B和Md e t之间的拟合关系可以写为Md e t=0.9 9 6MB B+0.0 4 1(1 0)选择多个辐射出射度标定点对热像仪的灰度误差进行标定。在第i个辐射出射度标定点处,通过灰度-辐射出射度标定曲线得到该点的测量值Md e ti。同时,根据此时低温黑体的实际温度值,计算出低温黑体的辐射出射度理论值MB Bi。热像仪辐射出射度的误差Mi可以写为Mi=Md e ti-MB Bi(1 1)77图4 辐射出射度标定实验结果F i g.4 R a d i a t i o ne m i s s i o nc a l i b r a t i o ne x p e r i m e n t r e s u l t s根据式(1 1)计算得到热像仪的辐射出射度的误差如图5所示。在低温黑体辐射出射度9.7 9 W/m2处,热像仪的辐射出射度误差最大,为-0.1 7 W/m2。图5 热像仪辐射出射度标定的误差(1 0 0.02 5 0.0 K)F i g.5 C a l i b r a t i o ne r r o r so f t h e r m a l i m a g i n gc a m e r ar a d i a t i o ne m i s s i o nd e g r e e s(1 0 0.02 5 0.0 K)根据灰度-辐射出射度标定曲线,得到1 0 0.02 5 0.0 K的低温范围内,热像仪的灰度误差如图6所示。热像仪的灰度误差最大为-9.9 1 D N,出现在灰度理论值为8 1 5.9 1 D N处。图6 热像仪灰度值标定的误差(1 0 0.02 5 0.0 K)F i g.6 C a l i b r a t i o ne r r o r so f t h e r m a l i m a g i n gc a m e r ag r a y s c a l ev a l u e s(1 0 0.02 5 0.0 K)4 结 论本文建立了基于红外热像仪灰度的低温标定模型,采用最小二乘法拟合,计算出辐射出射度响应函数和系统固有偏置的标定曲线。通过低温环境下红外辐射标定实验装置,得到温度范围低于热像仪2 5 3 K测温极限的灰度和辐射出射度之间的关系。分析了标定结果的理论值和实际测量值的误差,在低温黑体辐射出射度9.7 9W/m2处,热像仪的辐射出射度误差最大,为-0.1 7W/m2,此时的灰度误差为-9.9 1 D N。该技术适用于低温环境下非接触式的高精度红外辐射标定。参考文献(R e f e r e n c e s)1 K L I ONOV S KAK,B UR R IM.H a r d w a r e-i n-t h e-l o o ps i m u-l a t i o n sw i t hu m b r ac o n d i t i o n s f o r s p a c e c r a f t r e n d e z v o u sw i t hPMDv i s u a l s e n s o r sJ.S e n s o r s,2 0 2 1,2 1(4):1 4 5 5.2D E S T E F ANO M,M I S HR A H,G I O R D ANO A M,e ta l.A r e l a t i v ed y n a m i c sf o r m u l a t i o nf o r h a r d w a r e-i n-t h e-l o o ps i m u l a t i o no f o n-o r b i t r o b o t i cm i s s i o n sJ.I E E ER o b o t i c sa n dA u t o m a t i o nL e t t e r s,2 0 2 1,6(2):3 5 6 9-3 5 7 6.3L I E BMANN F.I n f r a r e dc a l i b r a t i o nd e v e l o p m e n ta tF l u k ec o r p o r a t i o n H a r ts c i e n t i f i cd i v i s i o nC/P r o c e e d i n g so fS P I E-T h eI n t e r n a t i o n a lS o c i e t yf o r O p t i c sa n d P h o t o n i c s,2 0 0 8:6 9 3 9 0 6.4U S AM E N T I A G AR,G A R C I ADF,I B A R R A-C A S T A N E D OC,e t a l.H i g h l ya c c u r a t eg e o m e t r i c c a l i b r a t i o n f o r i n f r a r e dc a m-e r a su s i n gi n e x p e n s i v ec a l i b r a t i o nt a r g e t sJ.M e a s u r e m e n t,2 0 1 7,1 1 2:1 0 5-1 1 6.5WAN G X,L IZ.A s p e c t r a l l yt u n a b l ec a l i b r a t i o ns o u r c eu s i n gE b e r t-F a s t i e c o n f i g u r a t i o nJ.M e a s u r e m e n tS c i e n c ea n dT e c h n o l o g y,2 0 1 8,2 9(3):0 3 5 9 0 3.6Y O OYS,K I MGJ,P A R KS,e t a l.S p e c t r a l r e s p o n s i v i t y c a l i-b r a t i o no f t h e r e f e r e n c e r a d i a t i o n t h e r m o m e t e r a tK R I S Sb yu s i n ga s u p e r-c o n t i n u u m l a s e r-b a s e d h i g h-a c c u r a c y m o n o c h r o m a t i cs o u r c eJ.M e t r o l o g i a,2 0 1 6,5 3(6):1 3 5 4-1 3 6 4.7P R E Z-I Z QU I E R D OJ,S E B A S T I N E,MA R T N E Z GM,e t a l.T h et h e r m a l i n f r a r e ds e n s o r(T I R S)o f t h eM a r se n v i r o n m e n t a ld y n a m i c sa n a l y z e r(ME D A)i n s t r u m e n to n-b o a r dM a r s2 0 2 0,ag e n e r a ld e s c r i p t i o na n dp e r f o r m a n c ea-n a l y s i sJ.M e a s u r e m e n t,2 0 1 8,1 2 2:4 3 2-4 4 2.8S E B A S T I NE,MA R T N E ZG,R AMO SM,e t a l.R a d i o-m e t r i ca n da n g u l a r c a l i b r a t i o nt e s t s f o r t h eME D A-T I R Sr a-d i o m e t e ro n b o a r dNA S A sM a r s2 0 2 0m i s s i o nJ.M e a s-u r e m e n t,2 0 2 0,1 6 4(1 0):1 0 7 9 6 8.9S H IQ,GAO Y,L IZ,e ta l.H i g h-d y n a m i c-r a n g e i n f r a r e dr a d i o m e t e rb a s e do nc h a o sd e t e c t i o nm e t h o dJ.I n f r a r e dP h y s i c s&T e c h n o l o g y,2 0 2 1,1 1 6(3):1 0 3 7 8 7.87飞控与探测第6卷第1期赵海涛,等:低温环境下红外辐射标定技术研究1 0S H IQ,GAOY,L IZ,e ta l.C a l i b r a t i o na n da p p l i c a t i o no fr a d i o m e t e rf o rl o w-t e m p e r a t u r ei n f r a r e d r a d i a t i o n s o u r c eC/P r o c e e d i n g so f1 0t hA p p l i e d O p t i c sa n dP h o t o n i c sC h i n a(AO P C2 0 2 1).B e i j i n g:S P I E,2 0 2 1.1 1 杨国庆.地基红外辐射测量系统提高测量精度的关键技术研究 D.长春:中国科学院长春光学精密机械与物理研究所,2 0 2 0.YAN G G Q.K e yt e c h n o l o g i e sr e s e a r c ho ni m p r o v i n gt h em e a s u r e m e n tp r e c i s i o no fg r o u n d-b a s e di n f r a r e dr a d i o m e t r ys y s t e m sD.C h a n g c h u n:C h a n g c h u nI n s t i t u t eo fO p t i c s,F i n eM e c h a n i c sa n dP h y s i c s,C h i n e s eA c a d e m yo fS c i e n c e s,2 0 2 0(i nC h i n e s e).1 2YUANL,HEZ,L YUG,e t a l.O p t i c a l d e s i g n,l a b o r a t o r yt e s t,a n dc a l i b r a t i o no fa i r b o r n el o n gw a v e i n f r a r e di m a g i n gs p e c t r o m e t e rJ.O p t i c sE x p r e s s,2 0 1 7,2 5(1 9):2 2 4 4 0-2 2 4 5 4.1 3L IC,D I N GY,L I NC,e t a l.O p t o m e c h a n i c a l d e s i g na n d s i m u-l a t i o no fa c r y o g e n i ci n f r a r e d s p e c t r o m e t e rJ.A p p l i e dO p t i c s,2 0 2 0,5 9(1 5):4 6 4 2-4 6 4 9.97