紅外焦平面探測器響應率不穩定像元快速篩選方法

崔 坤，陳凡勝，蘇曉鋒，梁清華，唐玉俊

?

紅外焦平面探測器響應率不穩定像元快速篩選方法

崔 坤1,2，陳凡勝2，蘇曉鋒2，梁清華1,2，唐玉俊1,2

（1. 中國科學院大學，北京 100049；2. 中國科學院上海技術物理研究所，上海 200083）

隨著紅外探測技術的發展，制冷型紅外焦平面探測器在航空航天和民用領域的應用越來越廣泛。探測器像元響應率的不穩定性對紅外焦平面探測器定量化應用產生非常大的影響。目前檢測響應率不穩定像元的方法主要通過定量化的黑體標定的方法來實現，該類方法需要改變黑體溫度從而得到像元的響應率，操作步驟和計算過程比較復雜，不易于高頻次的工程測試實現。本文提出了一種通過計算像元隨積分時間的相對響應率變化來快速篩選不穩定像元的方法，該方法每次測試時僅需要對著一個固定溫度的黑體，改變積分時間獲取一組探測器成像數據，通過計算每個像元的相對響應率，然后對比像元在不同測試時得到的相對響應率變化，通過與設定閾值的比較就可以實現響應不穩定性像元的篩選，自定義閾值還可以實現對響應不穩定像元進行分級。實驗證明，該方法能夠有效篩選出響應率不穩定的像元，操作簡便，便于工程實現。

紅外焦平面探測器；響應率不穩定；相對響應率；黑體標定

0 引言

紅外探測器像元響應率指在一定的幀周期或行周期條件下，紅外焦平面各像元對單位輻射功率產生的輸出信號電壓[1]。因此，像元量子效率和暗電流的改變都會引起像元響應率改變，另外像元的響應特性并非一成不變，而是隨著使用過程發生一定的變化[2]。經過長期測試發現，紅外焦平面中僅有少部分像元存在響應率變化較大的情況，該類像元的存在會影響探測器的定量化應用，尤其是弱小目標探測這類典型的應用環境，因此，根據像元響應率的變化對不穩定性像元進行篩選具有重要的意義。

在探測器像元響應穩定性研究方面，國外有研究機構做了大量的研究工作，其中斯必澤（Spitzer）科學中心的S. Carey、J. Surace和加州理工學院的S. Willner等人對斯必澤太空望遠鏡中的紅外探測組相機（IRAC）在軌運行的最初4.5年獲取的數據進行研究，他們發現其中一個8mm的探測器組件中過熱像元的數量以每年60個的速度在增長[3]。

2010年4月，JWST（James Webb Space Telescope，詹姆斯·韋伯太空望遠鏡）的NIR Spec研究小組發現原來測試合格的5mm截止波長的H2RG探測器經過約18個月的存儲，發生功能退化，通過像元暗電流測試的方法發現探測器內出現大量暗電流偏大的像元，暗電流從原來的小于0.01e－s－1，增大為0.1e－s－1～60 e－s－1。同年年底，NIR Cam研究小組發現同樣的問題出現在NIR Cam所用的4個5mm截止波長的H2RG探測器中[4]。

目前，國內外研究像元穩定性的方法主要為黑體標定法，暗電流測試法[4]，NETD[2,5]等。黑體標定法需要改變黑體溫度，并嚴格計算每個像元接收到的輻射通量，計算比較復雜，并且受實驗環境的影響較大。暗電流測試和NETD方法是一種間接評價像元響應穩定性的方法，操作步驟比較復雜。以上兩種方法都不利于高頻次的工程測試實現。相對改變黑體溫度，改變成像時的積分時間則變得非常容易，無論是像元量子效率的改變還是暗電流的改變都會影響像元隨積分時間變化的響應率，因此本文使用變積分時間的相對響應率變化來篩選不穩定性像元。文中首先對探測器內像元的響應模型進行分析，然后給出變積分時間像元響應率變化的定義，之后給出文中提出方法的原理和步驟，最后通過實驗對該方法的可行性進行驗證。

1 像元篩選方法原理

1.1 單像元響應模型

CMOS型碲鎘汞紅外焦平面探測器像元主要由光敏元、單元讀出電路、列選通電路及輸出緩沖器組成。其中，CTIA單元讀出電路結構是一種使用廣泛的讀出電路結構，該結構具有很多的優點：能提供穩定的探測器偏置，且有高的光電流注入效率、高增益和低噪聲[6]。因此，文中以CTIA結構為例分析探測器像元的輸出響應模型。單像元電路模型如圖1所示。

圖1 紅外探測器單像元模型

圖1中，光敏元完成光電轉換，積分電路實現電流電壓的轉換，為了便于分析，將這兩部分單獨出來分析。另外，光電二極管采用光生電流ph、電阻d和電容d來等效。該兩部分的等效電路如圖2所示。

圖2 像元光電轉換電路

圖2中int為積分電容，rest是積分復位MOS管，ref是放大器正向輸入端的參考電壓，det是光電二極管的端電壓，運算放大器的增益為c。

從(2)式可以看出，當運放近似理想運放的時候， CTIA結構像元輸出與積分時間呈線性關系。不過在實際設計時，并非必須將運放的放大器倍數設計的極大，通過一些電路結構優化方法同樣可以有效降低響應輸出的線性度，文獻[8]使用電路結構優化的方法使得CTIA型單元讀出電路的非線性度小于1%。目前，工程應用的CTIA結構的輸出非線性度都非常小。因此可以近似采用(2)式進行像元響應特性分析。

圖2中，暗電流可以表示為：

因此，根據(2)和(3)式，再將列選通電路和輸出緩沖器及噪聲考慮進去，像元最終的響應輸出電壓變為：

式中：1、2、3分別是源跟隨器、列緩沖器和輸出緩沖器的增益。根據光電轉換的原理，光生電流正比于像元接收到的波段輻射通量，因此光生電流可以表示為：

ph＝(5)

目前紅外焦平面探測器像元的讀出電路結構有多種形式，如自積分型（SI）、直接注入型（DI）等，雖然這些電路結構形式不一致，不過對于電壓型輸出的探測器而言，一般都是利用光生電流（包括暗電流）在電容上積分實現電流電壓的轉換，因此電壓輸出形式和CTIA結構類似，可用如下模型來表示：

1.2 像元篩選方法的原理

由單像元響應模型可知，當像元的量子效率、暗電流改變時，即使像元接收同樣的波段輻射通量，像元輸出電壓隨積分時間的斜率也會改變。因此可以通過分析每個像元隨積分時間的響應率來篩選響應不穩定像元。文獻[1]中給出了像元隨輻射功率的響應率的計算公式如式(8)所示：

由于像元接收波段輻射通量的計算受多個因素的影響，嚴格計算比較困難，因此采用計算相對響應率的方式來替代直接計算絕對響應率，即：

那么，由(10)和(11)式可以得到相對響應率的計算公式如下：

對于正常使用環境中的正常響應的像元而言，光生電流一般會遠大于暗電流，因此可以得到：

在不考慮輻射非均勻性的情況下，由(12)～(14)式可以得到：

從(15)式可以看出，在一定的波段輻射通量的范圍內，像元隨積分時間的相對響應率不再受波段輻射通量的影響，因此在實際測試中不需要嚴格計算像元接收到的輻射通量，只需要在測試過程中黑體的溫度穩定即可。另外，由于像元的光子噪聲跟像元接收到的輻射功率有關，因此在兩次不同測試時，最好采用同一黑體溫度來降低噪聲對測試精度的影響。

相對響應率僅能夠反應一次測試中像元的響應特性，為了方便對多次測試結果進行對比，文中采用隨積分時間的相對響應率變化來定量化分析像元穩定性的變化情況。像元的相對響應率變化定義為：對于同一像元，兩次實驗測得的相對響應率的差值除以第一次實驗測得的相對響應率，然后取百分比，具體計算公式如下：

另外，可以通過像元相對響應率變化來對不穩定像元進行分級，從而實現紅外面陣探測器不穩定像元的快速分級篩選。

2 像元篩選方法步驟

步驟1：設置面源黑體的溫度略高于測試環境溫度，在實驗室條件下建議黑體溫度采用25℃，調節黑體與探測器的相對位置，使黑體覆蓋探測器全視場，記下黑體與探測器的相對位置，下次測試需要保持相對位置不變（主要消除輻射非均勻性的影響）。待黑體溫度穩定以后，紅外面陣探測器測試系統上電，逐漸增大成像時的積分時間使探測器所有像元的輸出電壓均值逐漸增大到飽和輸出電壓，并記下此時的積分時間。

步驟2：細分積分時間序列，從0.1～0.9線性得到15個積分時間，然后逐漸改變積分時間，得到一組變積分時間的成像數據，每個積分時間下獲取100幀連續圖像數據。

步驟3：將每個積分時間下獲取的多幀數據通過累加平均的方式得到一幀圖像，降低測試系統噪聲。

3 實驗驗證

實驗采用320×256的紅外中波面陣探測器，其像元采用CTIA讀出電路結構。實驗測試方案示意圖如圖3所示。

圖3 實驗測試方案示意圖

通過表1可以看出，當探測器輸出在30%～75%范圍內變化時，有99.92%的像元的相對響應率變化小于5%，說明在一定的輻射通量范圍內，像元隨積分時間的相對響應率基本是不變的，和接收輻射通量近似無關，在每次實驗時可以用來分析像元響應穩定性的變化。

另外，為了驗證本文提出的方法能夠有效地篩選出響應率不穩定的像元，根據本文第3部分給出的像元篩選步驟，對該探測器進行2次測試，測試間隔大于10天，取出相對響應率變化的絕對值大于5%的5個像元。作為對比，在每次測試時，固定積分時間為2ms，改變黑體溫度獲取一組變黑體溫度數據用于計算像元隨輻射通量的響應率和響應截距。由于像元隨積分時間和輻射通量的響應線性度較好，在10%～90%動態范圍內通過線性擬合得到每個像元隨積分時間或者輻射通量的相對響應率和截距（數據采集電路采用16bit AD轉換器），如表2、表3所示。

表1 不同輸出條件下滿足閾值條件的像元比例

為了更清晰地顯示像元響應特性的變化，下文給出4號像元在兩次試驗中，隨積分時間和隨黑體溫度變化的響應曲線如圖4、圖5所示。圖4、圖5中，兩次試驗的均值響應曲線近似一致。

通過表2、表3可以看出篩選出的像元隨積分時間的相對響應率變化和隨輻射通量的相對響應率變化有相同的變化趨勢，并且幅值接近，這與式(7)相符，同時也說明這些像元的暗電流變化不大。由于計算像元隨積分時間的相對響應率變化比較容易，可以利用該方法作為一種預篩選的方法，對有漂移趨勢的像元做進一步跟蹤和監測。由式(7)可以看出，像元隨積分時間的響應率還能反映像元暗電流的變化。與此不同，暗電流的變化恰恰體現在像元隨輻射通量輸出的截距上。

另外，本文提出的方法是在像元輸出為線性的情況下得到的，對于部分單元讀出電路結構本身具有一定的非線性，此外對于像元信號輸出通路上的源跟隨器和輸出緩沖器而言，當輸入電壓較小時或者近似等于供電電壓時，輸出都會呈現一定的非線性。這也是目前多數焦平面探測器像元響應輸出呈現“S”型的原因之一。因此，本文提出的方法一般使用10%～90%動態范圍計算像元響應率。

4 結論

本文提出了一種紅外焦平面探測器響應不穩定像元的快速篩選方法，該方法以像元隨積分時間的相對響應率變化作為篩選依據。文中在CTIA結構的單像元響應輸出模型分析的基礎上給出了像元篩選方法的原理和具體實施步驟。最后，通過實驗驗證本文提出的方法能夠有效篩選出響應率變化的像元，并且該方法操作和計算方便，易于工程實現。

表2 隨積分時間相對響應率和截距

表3 隨輻射通量的相對響應率和截距

圖4 像元隨積分時間的響應曲線

圖5 像元隨輻射通量的響應曲線

[1] 丁瑞軍, 梁平治, 唐紅蘭, 等. GB/T 17444-2013. 紅外焦平面陣列參數測試方法[S].北京:中國標準出版社,2014. DING Ruijun, LIANG Pingzhi, TANG Honglan, et al. GB/T 17444-2013. Measuring methods for parameters of infrared focal plane arrays[S]. Beijing: Standards Press of China, 2014.

[2] Greiner M E, Smith R L, Timlin H A. Uniformity and stability in 2-dimensional infrared focal plane arrays [C]//, 1994, 2225:176-184.

[3] Carey S, Surace J, Lacy M, et al. Stability of the infrared array camera for the Spitzer space telescope[C]//, 2008, 7010:10-17.

[4] Bernard J Rauscher, Carl Stahle, Robert J Hill. Commentary: JWST near-infrared detector degradation-finding the problem, fixing the problem, and moving forward[J]., 2012, 2 (2): 021901.

[5] State-of-the-art MCT IR-modules with enhanced long term and cycle stability[C]//, 2012, 8353: 83532L-1.

[6] 周世椿. 高級紅外光電工程導論[M]. 北京: 科學出版社, 2014: 227-228.

ZHOU Shichun.[M]. Beijing: Science Press, 2014: 227-228.

[7] 文勇, 劉思超, 金友山, 等. CTIA型讀出電路非線性的研究[J]. 激光與紅外, 2009, 39(9): 978-981. WEN Yong, LIU Sichao, JIN Youshan, et al. Investigation on the nonlinearity of CTIA readout circuit[J]., 2009, 39(9): 978-981.

[8] 孫偉鋒, 宋文星, 謝亮, 等. CTIA型單元讀出電路非線性分析及優化設計[J]. 航空兵器, 2010, 2010(6): 48-52. SUN Weifeng, SONG Wenxing, XIE Liang, et al. Analysis and optimization on nonlinearity of CTIA readout circuit[J]., 2010, 2010(6): 48-52.

A Fast Screening Method of Pixels with Unstable Response Rate in IRFPA

CUI Kun1,2，CHEN Fansheng2，SU Xiaofeng2，LIANG Qinghua1,2，TANG Yujun2

（1. University of Chinese Academy of Science, Beijing 100049, China 2. Shanghai Institute of Technical Physics, Chinese Academy of Science, Shanghai 200083, China;）

With the development of infrared technology, cooled IRFPA detectors are widely used in aerospace and civil fields. The instability of pixels’ response rate has a great influence on the quantitative application of IRFPA detectors. The methods of detecting pixels with unstable response rate are mainly through quantitative blackbody calibration currently, which needs to change the blackbody temperature so as to obtain the pixels’ response rate, and the steps and calculations are very complicated. It is not easy for high-frequency engineering applications. This paper presents a fast screening method for pixels with unstable response rate by calculating the change of the pixels’ relative response rate with the integration time. In each test, it requires changing the integration time to obtain a set of imaging data with detector in front of a fixed temperature black body and calculating each pixel’s relative response rate. The change of each pixel’s relative response rate will be calculated. Then pixels with unstable response rate will be screened out through comparison with a setting threshold, and pixels can be classified by user-defined threshold. The experimental results show that the method can effectively screen out the pixels with unstable response rate, and the operation is simple and convenient, and it is easy to be realized.

IRFPA detectors，unstable response rate，relative response rate，blackbody calibration

TN215

A

1001-8891(2017)02-0130-06

2016-06-18；

2016-10-27.

崔坤（1989-），男，博士生，主要從事空間遙感圖像信息處理方面的研究。E-mail：cuikun_1989@163.com。

陳凡勝（1978-），男，博士，研究員，博士生導師，主要從事空間遙感方面的研究。E-mail：cfs@mail.sitp.ac.cn。

中國科學院上海技術物理研究所創新基金（CX-60）。