一種紅外成像系統(tǒng)的畸變測試方法

王 棪,郭哲民，劉國平

一種紅外成像系統(tǒng)的畸變測試方法

王 棪,郭哲民，劉國平

（昆明物理研究所，云南 昆明 6550223）

紅外成像系統(tǒng)圖像畸變控制的良好與否對其能否在應(yīng)用平臺上發(fā)揮應(yīng)有作用極為關(guān)鍵，因此在實(shí)驗(yàn)室內(nèi)對紅外成像系統(tǒng)的圖像畸變進(jìn)行測試、分析是極其必要的，本文提出了一種結(jié)合質(zhì)心亞像素識別和精密測角的局部畸變測試方法，讀取點(diǎn)靶在各個局部視場的質(zhì)心變化并進(jìn)一步計(jì)算即可得到畸變。采用該方法可以較好地解決點(diǎn)列質(zhì)心法測試大視場系統(tǒng)操作性不高以及局部畸變法定位不夠精確的問題，采用本方法對某型號紅外成像系統(tǒng)進(jìn)行了局部畸變測試，取得了與理論較為相符的測試結(jié)果，其相對畸變測試誤差不超過0.02%，可以很好地滿足紅外成像系統(tǒng)的畸變測試需求，并對成像系統(tǒng)性能進(jìn)行評估，有效反饋圖像畸變對紅外成像系統(tǒng)探測能力的影響。

畸變；紅外成像系統(tǒng)；質(zhì)心

0 引言

當(dāng)前，紅外成像系統(tǒng)作為輔助作戰(zhàn)的重要武器裝備，正被越來越多地應(yīng)用到各類作戰(zhàn)場景，其成像質(zhì)量的優(yōu)劣對能否實(shí)現(xiàn)精確打擊至關(guān)重要。其中，圖像畸變作為一個關(guān)鍵指標(biāo)，它雖不直接影響圖像清晰度，但會改變視場內(nèi)目標(biāo)的形狀，使圖像發(fā)生明顯的拉伸或壓縮，實(shí)戰(zhàn)中，這類目標(biāo)形狀的改變將會干擾觀察者對目標(biāo)性質(zhì)的判定。參考機(jī)載光電系統(tǒng)，視場可達(dá)20°～30°，當(dāng)在極限觀測距離上搜尋、定位目標(biāo)時，目標(biāo)形狀的改變將影響對目標(biāo)的識別概率，降低目標(biāo)被探測和識別的可能性，換言之，可能導(dǎo)致部分使用者產(chǎn)生誤判，進(jìn)而給作戰(zhàn)代入干擾，因此，如何對其在實(shí)驗(yàn)室進(jìn)行精確檢測極為關(guān)鍵。

在對紅外成像系統(tǒng)的畸變進(jìn)行檢測時，常需要考察若干離散的局部視場區(qū)域的畸變情況，要分析不同視場范圍的畸變大小，有兩種可行的思路，第一種，是通過對大量離散點(diǎn)數(shù)據(jù)進(jìn)行采集，擬合出適用于全視場范圍的畸變函數(shù)，再將目標(biāo)視場值代入其中，解算出離散的若干個局部畸變的大小，此方法可用于相機(jī)的畸變標(biāo)定和修正[1-3]。第二種，則是分別對各個離散的局部視場進(jìn)行數(shù)據(jù)采集，解算得到對應(yīng)的局部視場區(qū)域畸變大小[4-5]。

對國內(nèi)外常用的紅外成像系統(tǒng)畸變檢測技術(shù)進(jìn)行調(diào)研發(fā)現(xiàn)，第一種思路有一定的應(yīng)用，如美國的OPTIKOS公司，采用若干個點(diǎn)組成的陣列作為目標(biāo)，將其置于準(zhǔn)直儀焦面，為被測相機(jī)輸出目標(biāo)，通過分析局部視場區(qū)域內(nèi)的點(diǎn)列的質(zhì)心極矩[4]，得到畸變大小。采用第一種思路的缺點(diǎn)在于，不同視場大小的被測紅外系統(tǒng)，需要采用不同規(guī)格的點(diǎn)陣去匹配，效率上有欠缺，另外，當(dāng)被測相機(jī)視場較大時，采用點(diǎn)陣往往只能覆蓋相機(jī)的局部區(qū)域，若采用焦距更短的大視場平行光管，則其像質(zhì)可能無法滿足測試精度要求，實(shí)際中性能檢測所選用的準(zhǔn)直儀焦距至少是被測相機(jī)的3～5倍[4]。因此采用準(zhǔn)直儀輸出目標(biāo)覆蓋被測相機(jī)視場的方式往往不可行。參考可見光相機(jī)廣角鏡頭的畸變檢測方法，在有限距離的豎直平面，如光滑墻面上，制作大面積的棋盤目標(biāo)[5-6]，由于紅外成像系統(tǒng)的調(diào)焦范圍有限，最小清晰距離往往在數(shù)十米以上，采用這類方法將需要巨大的空間，客觀條件較為苛刻，且紅外成像系統(tǒng)的靶標(biāo)是通過能量強(qiáng)弱來區(qū)分，同時，畸變測試在靶標(biāo)的尺寸精度上要求較高，大面積的高精度熱靶在加工上有難度。

根據(jù)第二種思路設(shè)計(jì)的檢測方法更常見，國內(nèi)外的一些測試系統(tǒng)研發(fā)機(jī)構(gòu)，多采用國家軍用標(biāo)準(zhǔn)中描述的方式，將十字靶標(biāo)或矩形靶標(biāo)作為觀察目標(biāo)，其大小約占整個視場的5%，將其中心與光學(xué)系統(tǒng)主點(diǎn)重合時的圖形外徑作為基準(zhǔn)，將目標(biāo)圖像分別轉(zhuǎn)動至各個局部視場的邊緣，得到邊緣區(qū)域目標(biāo)外徑并與中心基準(zhǔn)進(jìn)行比較，得到其在水平和垂直方向的失真量從而解算出畸變。第二種思路對短焦大視場的紅外系統(tǒng)很適用，但隨著紅外成像技術(shù)的更新和迭代，一些長焦的紅外成像系統(tǒng)正越來越多地被應(yīng)用，對于此類紅外成像系統(tǒng)，上述的常用方法存在不足：采用具有一定外徑尺寸的目標(biāo)將其置于不同局部視場區(qū)域下，由于目標(biāo)本身占據(jù)了一定的視場范圍，因此解算得到的數(shù)據(jù)并不能嚴(yán)格對應(yīng)所考察的局部視場，視場定位的誤差帶來的影響將會被放大。

本文提出了一種基于質(zhì)心識別的局部畸變測試方法，測試目標(biāo)由均勻黑體輻射準(zhǔn)直儀靶標(biāo)輸出，通過測試軟件對小目標(biāo)進(jìn)行質(zhì)心亞像素識別，結(jié)合精密轉(zhuǎn)臺對角度的高精度控制，從而實(shí)現(xiàn)畸變測試，對于不同的被測系統(tǒng)，只需在特別定制的同一靶標(biāo)上，選取不同規(guī)格的圓孔即可，第一種思路在測試大視場時，用點(diǎn)列覆蓋大視場不易操作，采用本方法，只需通過轉(zhuǎn)動轉(zhuǎn)臺將所用靶標(biāo)圖像轉(zhuǎn)動至所需測試的視場方位并讀取目標(biāo)質(zhì)心即可，也無需考慮準(zhǔn)直儀視場與被測系統(tǒng)視場的匹配問題。第二種思路在測試長焦小視場系統(tǒng)時，容易因所用靶標(biāo)本身大小而影響到定位精度，采用本方法，只需將目標(biāo)質(zhì)心坐標(biāo)與需測試的局部視場位置重合即可，靶標(biāo)尺寸大小不影響局部視場的定位。測試結(jié)果的不確定度只依賴于測試設(shè)備的不確定度和質(zhì)心識別不確定度。通過這種方式，即保有質(zhì)心法的測試精度，也解決了局部畸變定位可能不準(zhǔn)確的問題。

1 測試原理

從畸變的原理入手，建立局部視場的畸變檢測模型，假設(shè)一個焦距為的紅外成像系統(tǒng)對準(zhǔn)直儀輸出的無限遠(yuǎn)目標(biāo)成像，目標(biāo)張角為，經(jīng)紅外成像系統(tǒng)成像的理想像高為￠，實(shí)際像高為。實(shí)際像高可能大于或小于理想像高，因此，畸變的符號將隨像高的增大或壓縮而改變，將畸變進(jìn)行泰勒展開[7]，將只包含奇數(shù)項(xiàng)，且展開通常取3項(xiàng)：

實(shí)際檢測中，普遍采用的是相對畸變的形式，則光學(xué)系統(tǒng)的相對畸變可表示為：

式(2)代入式(1)，相對畸變的級數(shù)展開形式為：

另外，像高與焦距和視場角相關(guān)，其表達(dá)式為：

＝×tan(4)

假設(shè)光學(xué)系統(tǒng)對一個9×9的點(diǎn)列目標(biāo)成像，其圖像如圖1所示，建立坐標(biāo)系，光學(xué)系統(tǒng)主點(diǎn)作為坐標(biāo)原點(diǎn)，點(diǎn)列目標(biāo)中心與主點(diǎn)重合，點(diǎn)列目標(biāo)的各個圓孔的理想質(zhì)心坐標(biāo)已知，采用紅外系統(tǒng)對點(diǎn)列成像，各個圓孔的實(shí)際質(zhì)心坐標(biāo)通過亞像素質(zhì)心識別算法得到，由于光學(xué)系統(tǒng)畸變所造成的質(zhì)心偏移即為對應(yīng)局部區(qū)域的相對畸變，綜合式(3)和式(4)，建立形如式(5)的方程組，通過求解該方程組，即可得到相對畸變的級數(shù)表達(dá)式。利用該表達(dá)式可以對全視場畸變進(jìn)行標(biāo)定：

圖1 畸變測試模型

通過上述模型，可以對全視場的畸變進(jìn)行檢測和標(biāo)定，考慮局部視場的檢測，如上文所述，全視場范圍的點(diǎn)列靶對一些視場較大的被測相機(jī)很難做到匹配，并且我們只針對離散的部分區(qū)域進(jìn)行分析，因此，無需建立分析連續(xù)數(shù)據(jù)的方程，從檢測的角度，我們也無需校正圖像畸變，因此，可對上述畸變模型進(jìn)行簡化：根據(jù)式(3)，相對畸變的大小隨視場增大而增大，對于同一局部視場區(qū)域，視場范圍的4個邊角位置畸變最大，在進(jìn)行畸變檢測時，同一局部視場，取最大值作為該局部視場的畸變值，依此作為判斷被測件畸變是否滿足設(shè)計(jì)要求的依據(jù)。

因此，可將局部視場畸變測試的模型簡化：采用單一圓孔取代多個圓孔組成的點(diǎn)列作為被測相機(jī)的目標(biāo)，如圖2所示，建立坐標(biāo)系，假設(shè)虛線矩形框標(biāo)注的是待測量的局部視場區(qū)域的邊界，假設(shè)需檢測3個局部視場的相對畸變，則可將圓孔分別置于3個區(qū)域的共12個邊界位置，邊界位置的理想坐標(biāo)根據(jù)被測相機(jī)的視場求得，利用軟件讀取圓孔在該位置的質(zhì)心坐標(biāo)，將實(shí)際質(zhì)心坐標(biāo)與理想坐標(biāo)相減，得到質(zhì)心偏移(,)，系統(tǒng)圖像的水平分辨率為，垂直分辨率為，系統(tǒng)水平視場為垂直視場為，根據(jù)式(6)，得到相對畸變：

在對全部12個邊角區(qū)域執(zhí)行檢測后，針對每個區(qū)域得到的4組數(shù)據(jù)，取其最大值作為該局部視場的畸變值。

采用這種辦法，畸變解算的基礎(chǔ)是質(zhì)心偏移量，與質(zhì)心分析的目標(biāo)外徑無直接關(guān)聯(lián)，因此，采用這種方法，可以嚴(yán)格地對應(yīng)局部視場區(qū)域，解算結(jié)果即為該局部視場的畸變值。

圖2 局部視場測試模型

根據(jù)上述模型，設(shè)計(jì)了紅外成像系統(tǒng)畸變檢測光路，如圖3所示，目標(biāo)在均勻輻射源照射下，經(jīng)過平面反射鏡和離軸拋物鏡放大后，以平行光輸入被測相機(jī)。

圖3 紅外成像系統(tǒng)畸變測試光路

2 測試過程

根據(jù)圖3所示的光路，配置相應(yīng)組件，搭建了測試平臺，如圖4所示，測試平臺包括目標(biāo)模擬、角度控制和被測件3個模塊，選擇某紅外成像系統(tǒng)對目標(biāo)模擬模塊提供的輸出成像，并輸入至高分辨率圖像監(jiān)視器和測試軟件。

圖4 測試平臺

根據(jù)上文描述的測試方法，制定測試步驟并執(zhí)行：1）調(diào)整轉(zhuǎn)臺，將圓孔目標(biāo)像置于成像系統(tǒng)圖像的左邊緣，清零轉(zhuǎn)臺讀數(shù)，再轉(zhuǎn)置右邊緣，讀取轉(zhuǎn)臺數(shù)據(jù)作為相機(jī)視場數(shù)據(jù)；2）對各個局部視場的待測區(qū)域進(jìn)行編號，如圖5，根據(jù)視場數(shù)據(jù)計(jì)算各編號位置對應(yīng)的方位和俯仰角度，以度分秒作為單位，數(shù)據(jù)見表1。3）調(diào)整轉(zhuǎn)臺使成像系統(tǒng)與平行光管光軸平行，打開成像系統(tǒng)的十字分化，將圓孔目標(biāo)與其中心保持重合，清零轉(zhuǎn)臺；4）按照計(jì)算出的局部視場角度轉(zhuǎn)動轉(zhuǎn)臺，并一一采集圖像存儲，目標(biāo)的移動軌跡見圖6；4）采集完圖像后由軟件解算出各個局部視場的質(zhì)心數(shù)據(jù)并代入公式(7)解算得到相對畸變。

圖5 局部視場編號

表1 局部視場對應(yīng)角度

圖6 目標(biāo)移動軌跡

3 測試結(jié)果

依據(jù)上述測試流程，測試得到了該紅外成像系統(tǒng)的畸變數(shù)值，為分析測試結(jié)果的合理性，利用軟件對光學(xué)鏡頭畸變進(jìn)行了診斷，取出3個局部視場的畸變最大值轉(zhuǎn)換為百分比，擬合曲線，該曲線和實(shí)測結(jié)果對比如圖7所示，可看到，兩者的變化趨勢是相吻合的，隨著局部視場區(qū)域的增大，畸變值逐漸增大，利用本測試方法得到的畸變測試結(jié)果其規(guī)律符合畸變指標(biāo)的規(guī)律。

圖7 仿真實(shí)測數(shù)據(jù)對比

為進(jìn)一步評估畸變測試結(jié)果對紅外成像系統(tǒng)成像質(zhì)量判斷所起到的作用，我們進(jìn)行了一組性能測試，選取一定空間頻率的四桿靶作為目標(biāo)，在不同局部視場、不同畸變區(qū)域?qū)ζ溥M(jìn)行MRTD（minimum resolvable temperature difference）測試，通過測試驗(yàn)證：MRTD測試結(jié)果的變化趨勢是否與畸變測試結(jié)果變化趨勢相符，即畸變的測試結(jié)果是否能一定程度衡量系統(tǒng)對目標(biāo)的探測能力；實(shí)際結(jié)果如圖8所示，測量結(jié)果與畸變測試數(shù)據(jù)變化趨勢吻合，畸變越大，MRTD測試結(jié)果越大，也即目標(biāo)觀測效果雖畸變的增大而下降，同時，畸變測試結(jié)果和MRTD測試值在3個局部視場均小于MRTD預(yù)估值和標(biāo)準(zhǔn)值，由此可知，畸變測試結(jié)果既符合目標(biāo)觀測的客觀變化趨勢，也能一定程度地反應(yīng)系統(tǒng)的探測能力變化，畸變結(jié)果可用于評估圖像失真是否影響到系統(tǒng)的探測能力，有效反饋系統(tǒng)圖像畸變對探測能力的影響。

圖8 MRTD測試值對比

4 誤差分析

本方法的測試誤差源主要包含：1）亞像素級質(zhì)心識別誤差；2）準(zhǔn)直儀輸出平行度誤差；3）高精度轉(zhuǎn)臺角度誤差。

準(zhǔn)直儀平行度誤差主要源于平行光管的焦面離焦和平行光管系統(tǒng)波像差，我們采用的準(zhǔn)直儀光學(xué)系統(tǒng)，其鏡頭組件包含離軸拋物鏡和平面反射鏡，該系統(tǒng)的平行性誤差主要來源有兩類：1）拋物鏡和平面鏡鏡面組件的面形誤差；2）拋物面和其焦面的相對位置的失調(diào)誤差。這兩類誤差的存在將會導(dǎo)致平行光管出射波前與理想平面波存在差異，從而輸出光線存在角像差，這類誤差將會對角度測量結(jié)果造成影響。此類偏差記為c，采用五棱鏡法對平行光管進(jìn)行標(biāo)定，平行光管的平行度極限誤差為＝32，包含因子＝2，則平行光管的平行度誤差為：

1＝/＝1.5002

高精度轉(zhuǎn)臺的誤差主要源于轉(zhuǎn)臺轉(zhuǎn)動的絕對角度與真值的偏離，此類偏離記為r，采用23面棱體和自準(zhǔn)直管對高精度轉(zhuǎn)臺進(jìn)行標(biāo)定，高精度轉(zhuǎn)臺的測試半寬為12，假設(shè)其服從均勻分布，其誤差為：

亞像素級質(zhì)心識別，誤差主要源于噪聲對算法處理的影響，噪聲的不穩(wěn)定性和隨機(jī)性將會使圖像目標(biāo)的灰度值產(chǎn)生波動，目標(biāo)灰度值矩陣作為質(zhì)心解算的輸入，其數(shù)值的變化將影響質(zhì)心解算的結(jié)果，此類影響記為z。

根據(jù)圖3所示的測試平臺，執(zhí)行分析質(zhì)心識別誤差分析的實(shí)際測試和圖像采集，紅外成像系統(tǒng)對平行光管無限遠(yuǎn)圓孔目標(biāo)成像，保持被測紅外成像系統(tǒng)與平行光管相對位置固定，將紅外成像系統(tǒng)增益調(diào)至最高級，使噪聲干擾最大化，連續(xù)采集10幀圖像，采用質(zhì)心亞像素識別算法對其進(jìn)行識別并解算，其質(zhì)心識別的向和向標(biāo)準(zhǔn)差均不超過0.02pixel，優(yōu)于0.1個像素。

綜合以上分析，在進(jìn)行相對畸變檢測時，平行光管的平行度和高精度轉(zhuǎn)臺精度將會影響視場測試精度和局部視場的定位精度，質(zhì)心識別精度將會影響質(zhì)心的識別，每次測試，兩者均會造成影響，視場角的角度誤差用式(8)計(jì)算，角度測試的誤差為1.802。

則相對畸變的誤差為：

將上文數(shù)據(jù)代入式(9)，得到各個局部視場的畸變測試誤差，數(shù)據(jù)見表2，可看到，越接近中心視場對畸變要求越為嚴(yán)格，在30%視場，畸變測試的誤差為0.014%，滿足使用要求。

表2 畸變測試誤差

5 結(jié)論

本文提出了一種局部視場相對畸變的測試方法，該方法提出了一種解算畸變的數(shù)學(xué)模型，結(jié)合該數(shù)學(xué)模型和實(shí)際測試條件，搭建了畸變檢測系統(tǒng)，該系統(tǒng)能夠結(jié)合精密角度測量和質(zhì)心亞像素識別，較為可靠地輸出被測紅外成像系統(tǒng)的局部視場畸變，采用該方法，能夠有效地規(guī)避采用點(diǎn)列目標(biāo)、十字目標(biāo)和方孔目標(biāo)所出現(xiàn)的問題，其測量的誤差不超過0.02%，滿足常見紅外成像系統(tǒng)的指標(biāo)要求，能夠較為可靠地在實(shí)驗(yàn)室對被測相機(jī)的畸變和進(jìn)行檢測和評估，并能作為系統(tǒng)探測能力評估的有力參考。

[1] 安德里安·凱勒, 加里·布拉德斯基. 學(xué)習(xí)OpenCV3[M]. 北京: 清華大學(xué)出版社, 2018.

Adrian Kaehler, Gary Bradski.3:C++[M]. Beijing: Qsinghua University Press, 2018.

[2] ZHANG Haofeng, ZHAO Chunxia, LIU Jianfeng. Straight line feature based image distortion correction[J]., 2007, 6(2): 83-86.

[3] 吳國棟, 韓冰, 何煦. 精密測角法的線陣CCD相機(jī)幾何參數(shù)實(shí)驗(yàn)室標(biāo)定方法[J]. 光學(xué)精密工程, 2007, 15(10): 1628-1631.

WU Guodong, HAN Bing, HE Xu. Calibration of geometric parameters of line array CCD camera based on exact measuring angle in lab[J]., 2007, 15(10): 1628-1631.

[4] Gerald Holst.[M]. Washington: SPIE Press, 1998.

[5] Fryer J G, Brown D C. Lens distortion for close-range photo- grammetry[J]., 1986, 52(1): 51-58.

[6] Brown D C. Close-range camera calibration[J]., 1971, 37(8): 855-866.

[7] 郁道銀, 談恒英. 工程光學(xué)[M]. 北京: 機(jī)械工業(yè)出版社, 2011.

YU Daoyin, TAN Hengying,[M]. Beijing: China Machine Press, 2011.

A Method for Testing Distortion of an Infrared Imaging System

WANG Yan，GUO Zhemin，LIU Guoping

(Kunming Institute of Physics, Kunming 650223, China)

It is necessary to test and analyze the distortion of an infrared imaging system in a laboratory because the ability of the system to play a role on the platform is dependent on the distortion control of an infrared imaging system. We have proposed a method based on a test of sub-pixel precision and precise angle measurement to calculate distortion using centroid difference data. It is laborious to test the distortion of an infrared system with a large field of view using a method which requires the computation of a large number of coordinates of points. Furthermore, the location of the field view could be a problem for the partial distortion method. By using the proposed method, the aforementioned problems can be addressed. We applied the method to evaluate a thermal system and the results were in agreement with those of the simulation software. The error control was less than 0.02%,which can meet the requirement. Moreover, the performance of the system was evaluated to analyze how the distortion affects the ability of the system.

distortion, infrared imaging system, centroid

TN215

A

1001-8891(2021)11-1061-05

2021-06-30；

2021-08-03.

王棪（1991-），男，工程師，現(xiàn)從事紅外成像系統(tǒng)性能評估技術(shù)研究，E-mail：wywang_y@163.com