祿曉飛，盛 捷，趙 慧

?

紅外輻射測量系統(tǒng)外場標定方法及飛行目標亮度反演方法

祿曉飛，盛 捷，趙 慧

（酒泉衛(wèi)星發(fā)射中心，甘肅 酒泉 732750）

對于外場紅外輻射特性測量設(shè)備，標定建立了輸出灰度和輸入能量之間的對應關(guān)系。通過構(gòu)建中心紅外系統(tǒng)外場標定體系，分析了直接擴展源標定、間接擴展源標定、小面元黑體近距離成像標定的能量傳遞過程、標定試驗過程以及3種標定方法得到的系統(tǒng)響應度和系統(tǒng)底灰度的物理含義，指出了它們的相同和不同之處，進行了一致性分析。給出了采用標定數(shù)據(jù)計算飛行目標亮度的方法并提出了中距離外場紅外輻射測量實驗方法，對大氣透過率軟件和標定方法進行驗證。

外場標定；平行光管；小面源黑體；幅亮度；大氣透過率

0 引言

隨著反導武器技術(shù)的快速發(fā)展，飛行器突防性能越來越受到重視，對飛行器突防性能綜合評估的要求也越來越高[1]。飛行器突防性能綜合評估包括理論分析和仿真計算、地面靜態(tài)測試、彈道靶試驗以及飛行試驗條件下外場動態(tài)目標紅外特性測量等環(huán)節(jié)。其中外場動態(tài)紅外特性測量最接近實戰(zhàn)條件，是飛行器突防性能評估重要的參考依據(jù)。

在進行外場動態(tài)紅外特性測量試驗前后，需要對紅外輻射測量設(shè)備進行標定，完成2項功能：①選擇紅外系統(tǒng)的工作模式，包括衰減片和積分時間，使得既能捕獲目標又使得設(shè)備像素不飽和。②建立紅外測量設(shè)備的輸出和接收到的能量的關(guān)系，以能夠?qū)︼w行試驗時外場動態(tài)紅外特性測量數(shù)據(jù)進行處理，獲取飛行器的紅外輻射特征，包括輻射強度（亮度）、溫度場等。

現(xiàn)有文獻大多研究實驗室標定方式（例如間接擴展源標定方法），然后用實驗室標定參數(shù)對近距離靜態(tài)目標測量數(shù)據(jù)進行處理。例如，李云紅[2-3]等研究了標定與測量條件相同時的精確測溫技術(shù)。陸段軍等[4]研究了采用積分球?qū)δ曅偷厍蜉椛涮綔y儀進行標定的方法。張曉龍等[5]研究了非致冷長波紅外熱像儀的輻射標定方法。楊詞銀等[6]研究了大口徑紅外光電系統(tǒng)輻射定標及誤差分析方法。

然而對于外場試驗來講，實驗室標定數(shù)據(jù)處理方式并不適用，這是由兩方面造成的：①外場拍攝目標為高速飛行器，背景為天空，紅外輻射測量受到太陽、地面輻射、大氣的影響，環(huán)境十分復雜；②外場試驗的紅外輻射測量設(shè)備大多為機動設(shè)備，經(jīng)緯儀方艙較小，一般無法采用實驗室通用的平行光管標定方式（間接擴展源法）；部分固定式紅外輻射測量設(shè)備雖然可采用間接擴展源標定方式，然而實驗室的工作環(huán)境是可以精確控制的，外場標定試驗的工作環(huán)境很難精確控制。這就要求研究與外場標定試驗配套的標定數(shù)據(jù)處理方法，并且回答“外場標定數(shù)據(jù)處理結(jié)果能否直接應用于飛行試驗測量數(shù)據(jù)中？”本文針對外場紅外輻射測量設(shè)備標定方法多，標定環(huán)境和飛行試驗環(huán)境不同的情況，解決了把外場標定數(shù)據(jù)處理結(jié)果應用于飛行試驗測量數(shù)據(jù)中的難題。

1 外場紅外系統(tǒng)標定體系

紅外系統(tǒng)一般工作在線性區(qū)。通過外場標定建立輸入能量和輸出的關(guān)系。圖1顯示了紅外探測器輸出電壓和接收能量存在如下的關(guān)系。從圖中可看出紅外輻射測量系統(tǒng)與能量可分為3個區(qū)：線性區(qū)、非線性區(qū)和飽和區(qū)。

圖1 紅外探測器接受輻射能量與輸出的典型響應曲線

紅外系統(tǒng)的外場標定方式分為3類：直接擴展源法、間接擴展源法、黑體近距離成像方法[7-12]。針對不同類型紅外輻射測量設(shè)備要選擇相應的標定方式及配套設(shè)備。

1）對于固定紅外系統(tǒng)可建立活動板房，在其內(nèi)設(shè)置滑軌式平行光管，試驗前后采用間接擴展源法對其進行標定；試驗期間將板房移開。在試驗間隙將滑軌式平行光管移動到固定庫房中。

2）活動紅外系統(tǒng)根據(jù)飛行試驗試驗需要布設(shè)在不同的站點，為了保證設(shè)備的機動性，標定設(shè)備的尺寸往往不大，不能內(nèi)置小型平行光管，可使用小面元黑體近距離（幾十米）標定，也可使用直接擴展源進行覆蓋式標定。

3）部分活動紅外系統(tǒng)尺寸較大，而且口徑較小，可內(nèi)置小型平行光管，可使用平行光管對紅外探測器部分像素進行標定。

4）目前大面元黑體已經(jīng)開始在紅外系統(tǒng)中應用，然而大面元黑體的升降溫較慢，面元溫度一致性不易控制，且配套的附屬設(shè)備尺寸較大，不適合活動紅外系統(tǒng)。

2 間接擴展源標定方法

在望遠鏡主光學系統(tǒng)入瞳處放置出瞳面積大于并覆蓋主光學系統(tǒng)入瞳面積的大口徑平行光管，在平行光管焦面上放置均勻面源黑體，以此作為標準紅外輻射源進行標定。平行光管一般由3部分組成：離軸拋物鏡、車床。小面元黑體發(fā)出的紅外光線經(jīng)過平行光管后形成平行光，為了保證精度，需要使得平行光管出射光線與紅外系統(tǒng)的主光軸對準，如圖2所示。

平行光管標定光路見圖3。

2.1 平行光管成像關(guān)系

記平行光管的焦距為g，紅外系統(tǒng)清晰成像的焦距為，紅外探測器單像元尺寸為，面積為d＝2，其對應的小面元黑體面積為B。圖3中，小面元黑體處發(fā)射出紅外光線經(jīng)過平行光管后形成平行光，經(jīng)紅外系統(tǒng)在紅外探測器上成像在探測器處。d和B存在如式(1)的關(guān)系：

2.2 相機像素接收能量

黑體的輻射亮度為B，黑體B發(fā)射且進入探測器像元的能量為[8]：

d表示相機探測器像元d接收光線的亮度，則d接受的能量為（見圖4）：

圖2 紅外系統(tǒng)與平行光管標定示意圖

圖3 間接擴展源標定方法

圖4 探測器接受能量示意圖

令g表示紅外能量經(jīng)過平行光管的能量損失，即：

d＝gB

因此，

這里o，d，都是常數(shù)，只有B是可調(diào)的（通過調(diào)整黑體溫度）。在執(zhí)行試驗時，由于飛行器距離較遠，設(shè)備的焦距設(shè)為。當飛行器的亮度為時（不考慮大氣吸收），探測器中心像元接收的能量為：

如果接受的能量（表現(xiàn)為灰度）與標定實驗中接受的能量相等，則可得到飛行器的亮度。

2.3 平行光管透過率tg標校

標校方法為：

1）用小口徑紅外系統(tǒng)（與紅外輻射測量系統(tǒng)波段相同）依次對面源黑體在2個以上溫度點T測量，計算小口徑相機輻射定標系數(shù)。

令為小口徑紅外系統(tǒng)的響應度，則其對面源黑體定標測量計算公式為：

(T)＋

這里小口徑相機接收的面源黑體輻亮度(T)為：

式中：()為小口徑相機的光譜響應函數(shù)；1為第一輻射常數(shù)；2為第二輻射常數(shù)。

2）用小口徑紅外系統(tǒng)（與紅外輻射測量系統(tǒng)波段相同）依次對平行光管在2個以上溫度點T測量，計算小口徑相機輻射定標系數(shù)。

小口徑紅外系統(tǒng)對平行光管定標測量計算公式為：

式中：g為平行光管幅亮度。

根據(jù)小口徑相機對平管、面源黑體多個溫度的測量結(jié)果，以最小二乘法求解得到平管子系統(tǒng)各個波段的透過率g。

2.4 像素輻照度、入瞳輻照度和像素入瞳輻照度

像素輻照度定義為：

進入到單個像素的能量對應的入瞳輻照度為：

2.5 小面元平行光管實驗

考慮大氣對紅外能量的吸收，目前紅外輻射測量往往在3個波段進行測量：短波1～3mm，中波3～5mm，長波8～12mm。紅外測量設(shè)備本身吸收部分紅外能量，吸收特性由光譜響應()表示，其往往由衰減片檔位確定。因此探測器接受的亮度為：

式中：(,)表示溫度為的黑體在波長處發(fā)射的亮度。

由于飛行器在再入段發(fā)出強烈的能量，因此紅外輻射測量設(shè)備往往會設(shè)置多檔位的積分時間和衰減片。當紅外測量設(shè)備波段1～2、積分時間、衰減片檔位設(shè)定后，利用紅外測量系統(tǒng)（短波、中波、長波）對黑體經(jīng)過平行光管出射的能量進行測量，調(diào)整黑體溫度T，＝1,…,，建立輻射輸入與像素灰度G的關(guān)系：

假設(shè)某臺紅外輻射測量設(shè)備可測量3個波段的紅外圖像，積分時間有5檔，衰減片有3檔，因此共有3×5×3＝45種組合，￥、￥也有45組。

2.6 標定參數(shù)B￥物理含義分析

根據(jù)式(7)、￥的含義是：在外場的條件下，黑體為絕對零點時，探測器產(chǎn)生的灰度輸出，包括2部分：①紅外系統(tǒng)沒有任何能量進入時紅外系統(tǒng)的輸出灰度；②平行光管標定時，光學系統(tǒng)的熱輻射能量opt和平行光管的雜散輻射能量g進入探測器產(chǎn)生的輸出。

3 直接擴展源標定方式

用大面元黑體覆蓋在紅外設(shè)備入瞳上，探測器直接測量黑體發(fā)射的能量，可不考慮平行光管標定和小面元近距離標定引入的雜散輻射等不可控因素，見圖5。然而大面元黑體的制作比較困難，面元溫度一致性和高發(fā)射率難以完成。

入射到探測器上某一像素的能量由入瞳上所有點出射某一方向上能量的集合，不考慮紅外測量設(shè)備的損失，像素接受輻射能量的亮度為黑體的輻射亮度，為dBo/2。

圖5 直接擴展源標定原理示意圖

黑體中心與紅外輻射測量設(shè)備主光軸對齊，設(shè)置測量波段、積分時間、衰減片檔位后，調(diào)整黑體溫度，記錄測量圖像，建立黑體幅亮度和系統(tǒng)響應（灰度）的對應關(guān)系：

式(8)中系統(tǒng)響應度與平行光管標定系統(tǒng)響應度相同，但是系統(tǒng)底灰度0只包括和紅外系統(tǒng)熱輻射產(chǎn)生的能量。

4 黑體近距離成像標定

對于沒有平行光管的紅外輻射測量設(shè)備，可采用黑體近距離成像標定方法，這時不考慮大氣對黑體輻射能量的衰減。小面元黑體放置在距離紅外輻射測量設(shè)備處（一般為十幾米至幾十米，見圖6），調(diào)整和焦距f使得黑體在探測器上清晰成像，且至少充滿1個像素。根據(jù)成像關(guān)系，近距離黑體成像時探測器像面與焦平面不在同一個位置，探測器的接收立體角不為常數(shù)，而與黑體與望遠鏡的距離有關(guān)，f＝＋。近距離標定時中心像素接收的能量為：

黑體中心與紅外輻射測量設(shè)備主光軸對齊，設(shè)置測量波段、積分時間、衰減片檔位后，調(diào)整黑體溫度，記錄測量圖像，建立黑體幅亮度和系統(tǒng)響應（灰度）的對應關(guān)系：

4.1 標定參數(shù)BR物理含義分析

根據(jù)式(10)、B包括2部分：①紅外系統(tǒng)沒有任何能量進入時紅外系統(tǒng)的輸出灰度；②近距離標定時，光學系統(tǒng)的熱輻射能量、地面熱輻射和大氣程輻射進入探測器產(chǎn)生的輸出。因此B受外場影響很大，若外場的條件發(fā)生變化，外場測量的B也是不同的。圖7給出了某設(shè)備白天/夜間中波標定結(jié)果。很明顯，白天標定與夜間標定灰度數(shù)據(jù)相比，B值有一個提升。

4.2 K￥和KR的關(guān)系

對于相同的測量波段、積分時間、衰減片檔位，相同的系統(tǒng)狀態(tài)（、相同），近距離標定和平行光管標定的系統(tǒng)響應度存在如下關(guān)系：

由于f＞，故￥＞K。

同一黑體設(shè)置相同的溫度進行兩種標定實驗，平行光管實驗測得的灰度為￥，近距離成像測得灰度為G：

5 中距離紅外輻射測量實驗

將黑體設(shè)置在距離設(shè)備幾公里處，調(diào)整黑體的溫度，用紅外輻射測量設(shè)備對其進行測量。用設(shè)備提供的大氣透過率軟件計算大氣透過率，對測量目標進行反演，進而對標定方法、數(shù)據(jù)處理方法、大氣透過率計算方法進行驗證。

在穩(wěn)定的環(huán)境中，黑體與紅外系統(tǒng)的距離為，黑體清晰成像時紅外系統(tǒng)的焦距為d，設(shè)定多個溫度點T，紅外系統(tǒng)測量得到灰度G：

這里為黑體與紅外系統(tǒng)之間大氣的透過率。(14)式可以轉(zhuǎn)化為：

5.1 Bd值分析

紅外系統(tǒng)與黑體相距為，其間的大氣會輻射能量，因此值不但包括相機固有的紅外輻射和還應包括大氣的程輻射和地面熱輻射及反射太陽能量。

5.2 驗證大氣透過率計算精度

圖6 小面源黑體近距離標定原理示意圖

圖7 白天、夜晚兩次標定的中波輻射出射度與灰度關(guān)系曲線

5.3 驗證紅外數(shù)據(jù)反演方法

6 飛行器亮度反演

熱特征是識別真假彈頭的重要信息。通過目標紅外輻射參數(shù)的測量，根據(jù)黑體輻射定律，可以獲得目標的溫度特征、平均溫度變化率、譜分布特征、發(fā)射率-面積之積、姿態(tài)運動特征、輻射強度方向圖等重要特征。

6.1 飛行試驗時探測器像元接收能量

像素灰度為：

所以飛行器的波段幅亮度為：

飛行試驗時，紅外輻射測量設(shè)備鏡筒瞄準天空，fly值由光學系統(tǒng)的熱輻射、固有灰度、飛行器與相機之間的大氣輻射、飛行器反射太陽輻射能量產(chǎn)生。

6.2 大氣透過率

從飛行器來的輻射功率在到達紅外傳感器前，會被大氣中某些氣體有選擇地吸收，大氣中懸浮微粒能使光線散射，吸收、散射雖然機理不同，其作用結(jié)果均使輻射功率在傳輸過程中發(fā)生了衰減。大氣的衰減包括大氣分子的衰減（吸收和散射）、氣溶膠的衰減（吸收和散射）。吸收系數(shù)、散射系數(shù)均隨波長、路徑和大氣環(huán)境變化而變化。

根據(jù)Beer定律，單色距離大氣透過率可表示為：

為大氣組分的消光系數(shù)，主要包括：

ext＝a＋m＋m(21)

式中：a、m、m分別為氣溶膠粒子消光系數(shù)、大氣分子散射系數(shù)和吸收系數(shù)。

白天用多波段太陽輻射計測量可見光到近紅外多波段的大氣光譜透過率，其中1～2個波段可用來遙測整層大氣水汽總含量；夜晚用激光雷達測量大氣氣溶膠光學厚度和大氣能見度。用自動氣象站測量地表大氣溫度、濕度、氣壓，以及風速、風向。根據(jù)建立的當?shù)卮髿鈪?shù)數(shù)據(jù)庫和這些實時測量數(shù)據(jù)可以根據(jù)大氣透過率計算軟件（LowTran、ModTran、Cart）計算單色大氣透過率a()。

根據(jù)(18)式：

6.3 背景灰度的物理含義

飛行試驗測量得到的紅外圖像，目標占據(jù)較少的區(qū)域，目標周圍的區(qū)域稱為背景區(qū)域，背景區(qū)域的灰度b是由于光學系統(tǒng)的熱輻射、固有灰度、整個大氣層紅外輻射產(chǎn)生。

對于晚上執(zhí)行的試驗，沒有反射的太陽輻射，可認為fly是背景灰度。特別當飛行器的高程較高，由于外層大氣很稀薄且溫度低，目標與紅外系統(tǒng)之間的大氣輻射近似等于整個大氣層輻射。因此式(19)轉(zhuǎn)化為：

6.4 亮度反演

7 結(jié)論

本文研究了間接擴展源標定方法、直接擴展源標定方法、小面元黑體近距離成像標定方法的能量傳遞過程和標定試驗過程，分析了3種標定方法得到的系統(tǒng)響應度和系統(tǒng)底灰度的物理含義。最后又研究了飛行器的亮度反演方法，討論了波段內(nèi)大氣透過率的等效算法，為外場紅外輻射數(shù)據(jù)處理奠定了基礎(chǔ)。

表1 溫度-波段輻亮度查找表

[1] 徐南榮, 卞南華. 紅外輻射與制導[M]. 北京: 國防工業(yè)出版社, 1997.

[2] 李云紅, 基于紅外熱像儀的溫度測量技術(shù)及其應用研究[D]. 哈爾濱: 哈爾濱工業(yè)大學, 2010.

[3] 李云紅, 孫曉剛, 原桂彬. 紅外熱像儀精確測溫技術(shù)[J]. 光學精密工程, 2007, 15(9): 1336-1341.

[4] 陸段軍, 王模昌, 郁蘊健, 等. 凝視型地球輻射探測儀的輻射標定技術(shù)研究[J]. 紅外與毫米波學報, 2009, 28(1):42-45.

[5] 張曉龍, 劉英, 孫強. 高精度非致冷長波紅外熱像儀的輻射標定[J]. 中國光學, 2012, 5(3): 235-241.

[6] 楊詞銀, 曹立華. 大口徑紅外光電系統(tǒng)輻射定標及誤差分析[J]. 紅外與激光工程, 2011(9): 1624-1628.

[7] 唐嘉, 高昕. 靶場紅外輻射特性測量系統(tǒng)標定方法研究[J]. 紅外與激光工程, 2006(S1): 266-270.

[8] 張建奇, 方小平. 紅外物理[M]. 西安: 西安電子科技大學出版社, 2004.

[9] 陸子鳳. 紅外熱像儀的輻射定標和測溫誤差分析[D]. 長春: 長春光學精密機械與物理研究所, 2010.

[10] 王文革. 輻射測溫技術(shù)綜述[J]. 宇航計測技術(shù), 2005, 25(4): 20-32.

[11] 羅易雪, 湯心溢, 王宇, 等. 基于參考輻射源標定的紅外成像非均勻校正技術(shù)[J]. 紅外技術(shù), 2014, 36(4): 281-285.

[12] 范永杰, 金偉其, 朱麗紅. 紅外場景輻射模擬技術(shù)發(fā)展[J]. 紅外技術(shù), 2013, 35(3): 133-138.

Outdoor Calibration System of Infrared Device and Method of Computing Luminance of Aircraft

LU Xiao-fei，SHENG Jie，ZHAO Hui

（,732750,）

For outdoor infrared device, the calibration process establishes the relationship between input energy and output gray. This paper constructs the outdoor calibration system of infrared device, analyzes the energy transporting process of calibration by collimator, calibration by near blackbody, calibration by overlap blackbody, and demonstrates the process of calibration. By analyzing the physical meaning of calibration parameters, this paper points out the same aspects and different aspects of system response，based on which the method of computing luminance of aircraft is proposed. Also this paper proposes the method of measuring blackbody with middle range, which can be applied to validate the atmospheric transmittance and method of calibration.

outdoor calibration，collimator，small surface blackbody，luminance，atmospheric transmittance

TN216

A

1001-8891(2015)02-0154-07

2014-03-04；

2014-10-01.

祿曉飛（1981-），男，工程師，博士，主要從事紅外物理研究。E-mail：Luxf08@163.com。