張發(fā)強(qiáng)，張玉發(fā)，鄧 強(qiáng)，程正東，朱 斌

(1.脈沖功率激光技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，電子工程學(xué)院，安徽 合肥230037;2.電子工程學(xué)院基礎(chǔ)部，安徽 合肥230037)

1 引言

建立合理準(zhǔn)確的模型，是光電成像系統(tǒng)研究的一個(gè)重要方面，可以為光電成像系統(tǒng)的研制和評(píng)價(jià)提供理論依據(jù)和參考，具有重要的研究?jī)r(jià)值。光電成像系統(tǒng)建模是一種數(shù)學(xué)模型的構(gòu)建，需要綜合考慮目標(biāo)和背景特性、光電成像系統(tǒng)性能、大氣環(huán)境及多種其他因素。目前光電成像系統(tǒng)建模所采用的方法主要有以下四種:調(diào)制傳遞函數(shù)法、光線追跡法、點(diǎn)擴(kuò)散函數(shù)法和基于圖像像素處理法［1］。

以上幾種建模方法各有優(yōu)缺點(diǎn)，相對(duì)于其他幾種建模方法，調(diào)制傳遞函數(shù)法把光電成像系統(tǒng)中的光學(xué)系統(tǒng)、探測(cè)器、信號(hào)處理電路和顯示裝置等組成部分有機(jī)地結(jié)合起來(lái)，使模型具有普遍適用性，同時(shí)又可以將各個(gè)部分分開，獨(dú)立地對(duì)光電成像系統(tǒng)各組成部分進(jìn)行分析，具有開放、靈活的特點(diǎn)，該方法深入物理過(guò)程，囊括了大多數(shù)物理現(xiàn)象，所建立的模型能夠較準(zhǔn)確地描述真實(shí)系統(tǒng)。

本文以調(diào)制傳遞函數(shù)為理論依據(jù)，分別計(jì)算了光電成像系統(tǒng)各組成部分的調(diào)制傳遞函數(shù)，求得了系統(tǒng)總的調(diào)制傳遞函數(shù)，建立了基于MTF的光電成像系統(tǒng)模型。

2 調(diào)制傳遞函數(shù)的理論基礎(chǔ)

在利用MTF建模過(guò)程中，假設(shè)光電成像系統(tǒng)是線性、空間不變系統(tǒng)［2］。調(diào)制傳遞函數(shù)是光學(xué)傳遞函數(shù)(OTF)的模，描述的是系統(tǒng)對(duì)輸入各種頻率圖像的對(duì)比度傳遞能力，表征系統(tǒng)的頻率通過(guò)量，圖1說(shuō)明了調(diào)制傳遞函數(shù)的物理意義。

圖1 調(diào)制傳遞函數(shù)示意圖

從圖1可以看出，物體的對(duì)比度經(jīng)過(guò)光電成像系統(tǒng)后變小，圖像變得模糊，不同的空間頻率對(duì)應(yīng)的對(duì)比度變化曲線就是調(diào)制傳遞函數(shù)。

由信號(hào)處理理論可知，MTF是光電成像系統(tǒng)點(diǎn)擴(kuò)展函數(shù)(psf)在頻域的傅里葉變換，以二維圖像為例，利用psf和MTF表示光電成像系統(tǒng)的輸入和輸出之間的關(guān)系:

式中，x，y為空間坐標(biāo);νx，νy為空間頻率;式(1)和式(2)分別表示了輸出圖像與輸入圖像在空域及頻域上的關(guān)系?？梢钥闯?，通過(guò)調(diào)制傳遞函數(shù)可以將空域上的卷積轉(zhuǎn)化為頻域上的乘積，簡(jiǎn)化了運(yùn)算。

3 光電成像系統(tǒng)建模分析

如果一個(gè)系統(tǒng)是分立的，并且是線性平行不變(LSI)的，則整個(gè)系統(tǒng)的調(diào)制傳遞函數(shù)是各個(gè)部分調(diào)制傳遞函數(shù)的乘積。

根據(jù)以上的思想，光電成像系統(tǒng)可以作為L(zhǎng)SI系統(tǒng)進(jìn)行研究，求出各個(gè)組成部分的MTF，系統(tǒng)總的MTF為各個(gè)組成部分MTF的乘積，將系統(tǒng)總的MTF與輸入的圖像的傅里葉變換相乘，最后再做反傅里葉變換，并疊加合理的噪聲，即可得到系統(tǒng)輸出圖像，這就是建立基于MTF的光電成像系統(tǒng)模型的核心思想。

所以在利用調(diào)制傳遞函數(shù)對(duì)光電成像系統(tǒng)進(jìn)行建模之前，首先要分析光電成像系統(tǒng)的組成部分，并計(jì)算各部分的調(diào)制傳遞函數(shù)，通常光電成像系統(tǒng)包含以下幾個(gè)組成部分:目標(biāo)和背景、大氣傳輸、光學(xué)系統(tǒng)、探測(cè)器、信號(hào)處理電路、顯示裝置，其模型如圖2所示。

圖2 光電成像系統(tǒng)MTF模型

從圖2可以看出，光電成像系統(tǒng)的輸出圖像和輸入圖像之間的關(guān)系如下:

式中，σ(x，y)為疊加的噪聲，所以求出光電成像系統(tǒng)總的MTF是建模的關(guān)鍵問(wèn)題，下面對(duì)光電成像系統(tǒng)各組成部分的MTF進(jìn)行分析。

3.1 大氣的MTF

大氣對(duì)光電成像系統(tǒng)的影響主要來(lái)源于大氣分子散射、氣體吸收以及氣溶膠散射和吸收三個(gè)部分對(duì)目標(biāo)圖像傳輸?shù)挠绊?，吸收和散射引起到達(dá)傳感器上的光能量減少，而散射和湍流會(huì)引起圖像的模糊和細(xì)節(jié)損失。根據(jù)大氣對(duì)光傳輸?shù)牟煌绊?，將大氣調(diào)制傳遞函數(shù)分為湍流MTF和氣溶膠MTF［3－4］。其中湍流MTF又被分為長(zhǎng)曝光時(shí)間MTF和短曝光時(shí)間MTF，其表達(dá)式分別為:

式中，C2n為大氣折射率常數(shù)，通常采用經(jīng)驗(yàn)?zāi)Ｐ陀?jì)算;r為光傳輸?shù)穆窂介L(zhǎng)度;λ為波長(zhǎng);μ在用于近場(chǎng)時(shí)取1，用于遠(yuǎn)場(chǎng)時(shí)取0.5;D為光學(xué)系統(tǒng)的入瞳直徑，當(dāng)曝光時(shí)間t〈〈10 ms時(shí)，認(rèn)為其滿足短曝光條件［5］。

大氣的氣溶膠MTF分為低空間頻率MTF和高空間頻率MTF，其表達(dá)式如下［3］:

式中，νc是近似為a/λ的截止頻率;a是微粒半徑;Aa是大氣的吸收系數(shù);Sa是大氣散射系數(shù);從上式可以看出，大氣的高空間頻率MTF其實(shí)就是透過(guò)率，只有在低空間頻率時(shí)MTF才是空間頻率的函數(shù)。

通過(guò)上面的分析可以看出，光在大氣中傳輸時(shí)，會(huì)受到氣溶膠微粒的吸收和散射，以及湍流的影響，在工程計(jì)算中，假設(shè)氣溶膠微粒對(duì)輻射光傳輸?shù)挠绊懞屯牧鲗?duì)光傳輸?shù)挠绊懯仟?dú)立的，根據(jù)MTF的特點(diǎn)，大氣總的MTF為湍流MTF與氣溶膠MTF的乘積:

上式即為光電成像系統(tǒng)建模中大氣部分的調(diào)制傳遞函數(shù)表達(dá)式。

3.2 光學(xué)系統(tǒng)的MTF

光學(xué)系統(tǒng)是光電成像系統(tǒng)的重要組成部分，接收來(lái)自景物的光輻射，并把景物圖像投射到系統(tǒng)的探測(cè)器陣列上。

在實(shí)際分析光學(xué)系統(tǒng)的MTF時(shí)，應(yīng)該分析引起光學(xué)系統(tǒng)成像質(zhì)量變化的因素。不同的因素對(duì)應(yīng)不同的MTF，光學(xué)系統(tǒng)的調(diào)制傳遞函數(shù)包括衍射極限的調(diào)制傳遞函數(shù)MTFdif和像差產(chǎn)生的調(diào)制傳遞函數(shù)MTFabe兩部分組成［6］。

Gaskill給出了圓形孔徑的光學(xué)系統(tǒng)衍射限MTF:

式中，νc=D/λ是光學(xué)系統(tǒng)的截止頻率，單位為c/mrad;D是光瞳直徑，只有當(dāng)空間頻率小于截止頻率時(shí)，上式才成立，當(dāng)空間頻率大于截止頻率時(shí)，MTF為零。

由像差引起的彌散斑能量分布近似為高斯函數(shù)，具有圓對(duì)稱形式，其傳遞函數(shù)為:

式中，σ=σc/f是角度度量的標(biāo)準(zhǔn)方差，單位為c/mrad;f為光學(xué)系統(tǒng)的焦距。σ與彌散斑所占能量的百分比有直接的關(guān)系，假定所有的能量集中在彌散斑內(nèi)，則標(biāo)準(zhǔn)方差σc取彌散斑直徑1/4。

則光學(xué)系統(tǒng)總的MTF為:

3.3 探測(cè)器的MTF

探測(cè)器是光電成像系統(tǒng)的核心部件，它既完成視場(chǎng)內(nèi)圖像的掃描取樣，同時(shí)通過(guò)光電效應(yīng)將輻射信號(hào)轉(zhuǎn)化成電信號(hào)。探測(cè)器對(duì)輸入光信號(hào)的作用主要有三種:采樣、傳輸和擴(kuò)散，探測(cè)器總的調(diào)制傳遞函數(shù)等于這三種調(diào)制傳遞函數(shù)的乘積。探測(cè)器種類眾多，這里以常用的凝視型CCD探測(cè)器為例進(jìn)行分析。

這里只給1－D探測(cè)器陣列的具體推導(dǎo)情況，2－D陣列的推導(dǎo)是完全相似的，設(shè)其部分尺寸為a，采樣間隔d，根據(jù)采樣定律和探測(cè)器的規(guī)格可以得出探測(cè)器的采樣MTF［3］:

從上式可以看出，探測(cè)器的采樣要同時(shí)受到感光區(qū)尺寸a和采樣間隔d的影響。

CCD在進(jìn)行信號(hào)采樣的同時(shí)，將光信號(hào)轉(zhuǎn)換成電信號(hào)并在像元之間轉(zhuǎn)移傳輸，在傳輸過(guò)程中由于各種原因，必然存在損耗，由此而導(dǎo)致的MTF稱之為傳輸MTF。

設(shè)CCD像元的數(shù)目為n，像元間的電荷轉(zhuǎn)移損失率為ε，CCD的耐奎斯特頻率為νN，則CCD的傳輸調(diào)制傳遞函數(shù)為［7］:

從上式可以看出，nε越大，MTFt就下降的越厲害，所以為了提高CCD的傳輸MTF，減少電荷轉(zhuǎn)移損失率是至關(guān)重要的。

CCD中的光信號(hào)在進(jìn)行傳輸時(shí)，其產(chǎn)生的電荷會(huì)擴(kuò)散到相鄰的像元，引起MTF的下降，由于擴(kuò)散產(chǎn)生的MTF為［8］:

式中，α為材料對(duì)光的吸收系數(shù);d為表面耗盡區(qū)寬度;L0為光生電荷擴(kuò)散長(zhǎng)度;L是與空間頻率有關(guān)的量:

從上面兩式可知，隨著空間頻率的增加，MTFd(ν)的值下降。

通過(guò)上面的分析，可得探測(cè)器總的調(diào)制傳遞函數(shù)為:

上式即為光電成像系統(tǒng)建模中探測(cè)器部分的調(diào)制傳遞函數(shù)。

3.4 信號(hào)處理電路的MTF

信號(hào)處理電路是光電成像系統(tǒng)的重要組成部分，不同的光電成像系統(tǒng)有不同的信號(hào)處理電路，但都可以看作是各種濾波器的集合，在這里主要討論高通和低通濾波器的調(diào)制傳遞函數(shù)。在討論信號(hào)處理電路的MTF的時(shí)候，特別要注意電路所傳遞的是時(shí)域信號(hào)，其對(duì)應(yīng)的傳遞函數(shù)是時(shí)間頻率域，所以在實(shí)際的應(yīng)用中要將其轉(zhuǎn)換成空間頻率域，這樣才能和整個(gè)光電系統(tǒng)的MTF統(tǒng)一，其MTF表示如下［9］:

式中，MTFe1(ν)為低通濾波器的調(diào)制傳遞函數(shù);MTFe2(ν)為高通濾波器的調(diào)制傳遞函數(shù);νt0為濾波器的3dB頻率;νt0=1/2πRC;R和C為濾波器的電阻和電容;ν0為νt0對(duì)應(yīng)的空間頻率;ω為掃描角頻率，對(duì)不同的成像體制，ω有不同的表示方式。

整個(gè)信號(hào)處理電路可以看作是高通濾波器和低通濾波器的組合，一個(gè)濾波器的輸出是另外一個(gè)濾波器的輸入，從系統(tǒng)響應(yīng)的角度，認(rèn)為整個(gè)信號(hào)處理電路的MTF是高通濾波器的MTF和低通濾波器的MTF的乘積:

上式是從普遍意義上對(duì)光電成像系統(tǒng)的信號(hào)處理電路MTF進(jìn)行了分析，當(dāng)針對(duì)具體的系統(tǒng)時(shí)，要分析其電路構(gòu)成，求出其MTF值。

3.5 顯示裝置的MTF

顯示裝置的種類很多，其中最常用的是監(jiān)視器是陰極射線管CRT，其傳遞函數(shù)如下［10］:

式中，σm為光斑尺寸參數(shù)，可以采用收縮光柵技術(shù)測(cè)出。

3.6 系統(tǒng)總的MTF

光電系統(tǒng)總的調(diào)制傳遞函數(shù)(MTF)是系統(tǒng)各個(gè)組成部分調(diào)制傳遞函數(shù)的乘積，根據(jù)上面的分析，可得系統(tǒng)總的MTF為:

根據(jù)上式可以建立基于MTF的光電成像的模型，分析輸入和輸出圖像的關(guān)系，進(jìn)行系統(tǒng)仿真實(shí)驗(yàn)。

4 非線性噪聲影響分析

基于MTF光電成像系統(tǒng)模型的理論基礎(chǔ)是假定光電成像系統(tǒng)是線性的。然而在實(shí)際中，光電成像系統(tǒng)各組成部分都不同程度的有非線性效應(yīng)，為了提高模型的精度，常采用加入合理的噪聲模型模擬非線性效應(yīng)的作用。非線性效應(yīng)在光電成像系統(tǒng)中主要表現(xiàn)在以下幾個(gè)方面［11］:

(1)大氣中氣溶膠分布密度和湍流狀態(tài)隨機(jī)的非線性變化;

(2)光學(xué)系統(tǒng)透過(guò)率隨入射光光譜成分的非線性變化;

(3)探測(cè)器量子效率隨入射光非線性變化;

(4)信號(hào)處理電路中的非線性效應(yīng);

(5)顯示裝置隨光強(qiáng)的非線性變化。

在實(shí)際模型的建立中，為了反映這種非線性效應(yīng)對(duì)光電成像系統(tǒng)模型的影響，通過(guò)在MTF模型上疊加適當(dāng)?shù)脑肼暷Ｐ?。通常采用疊加三維噪聲模型的方法，三維噪聲模型能夠全面表征所有的噪聲源特征，具有普遍意義。通過(guò)把噪聲分為可以管理的數(shù)據(jù)集，使得復(fù)雜問(wèn)題簡(jiǎn)單化，能夠較好的反映光電成像系統(tǒng)的非線性效應(yīng)對(duì)輸出圖像的影響［12］，三維噪聲模型的數(shù)學(xué)表達(dá)式如下:

式中，σsys為系統(tǒng)總噪聲，其他7種噪聲含義如表1所示。

表1 三維噪聲模型中的噪聲分量描述

從表1可以看出，系統(tǒng)總的噪聲是由7種產(chǎn)生機(jī)理不同的噪聲組成，在具體的光電成像系統(tǒng)中并不是存在所有的噪聲，可能是其中幾種噪聲占主導(dǎo)地位，這和系統(tǒng)所采用的成像技術(shù)有密切的關(guān)系［13］。

5 仿真分析

根據(jù)上面調(diào)制傳遞函數(shù)的計(jì)算方法和光電成像系統(tǒng)的使用要求及性能指標(biāo)，在系統(tǒng)研制之前就可以建立其MTF模型進(jìn)行仿真模型，以評(píng)價(jià)和改進(jìn)光學(xué)成像系統(tǒng)，采用圖1所示的仿真思路，再疊加合理的噪聲，就可以建立光電系統(tǒng)的輸入輸出模型。

選擇標(biāo)準(zhǔn)圓形標(biāo)準(zhǔn)靶板為輸入圖像進(jìn)行模擬，采用Matlab編程將上述模型程序化，針對(duì)具體光電成像系統(tǒng)求得各組成部分的MTF，進(jìn)而求得系統(tǒng)總的MTF，疊加合理的噪聲模型，可得到系統(tǒng)的輸入和輸出的圖像，如圖3和圖4所示。下面以前期設(shè)計(jì)的一個(gè)紅外成像系統(tǒng)為例，對(duì)其進(jìn)行仿真分析，評(píng)價(jià)其系統(tǒng)質(zhì)量。

圖3 光電成像系統(tǒng)的輸入圖像

圖4 光電成像系統(tǒng)的輸出圖像

從上圖可以看出輸入圖像經(jīng)過(guò)光電成像系統(tǒng)以后，輸出圖像對(duì)比度和分辨率有所降低，并且變的有些模糊，但總體上對(duì)目標(biāo)進(jìn)行了良好的成像，能夠滿足發(fā)現(xiàn)和識(shí)別的要求。

6 結(jié)論

本文在分析了調(diào)制傳遞函數(shù)的基本理論之后，提出了建立基于MTF的光電成像系統(tǒng)的模型，通過(guò)分析光電成像系統(tǒng)中的各個(gè)環(huán)節(jié)的傳遞函數(shù)，最終計(jì)算得到了系統(tǒng)總的MTF，通過(guò)對(duì)輸入和輸出圖像進(jìn)行仿真分析，可以為光電系統(tǒng)的設(shè)計(jì)提供理論依據(jù)。

［1］ WANG Xuewei，LI Ke，WANG Shili．Research on modeling and simulation methods of infrared imaging system［J］．Computer＆Digital Engineering，2012，40(6):124－126．(in Chinese)王學(xué)偉，李珂，王世立．紅外成像系統(tǒng)建模仿真方法研究［J］．計(jì)算機(jī)與數(shù)字工程，2012，40(6):124－126．

［2］ LNaiguang．Fourier optics［M］．Beijing:CHINA MACHINE PRESS，2006:144－163．(in Chinese)呂乃光．傅里葉光學(xué)［M］．北京:機(jī)械工業(yè)出版社，2006:144－163．

［3］ Ronald G Driggers，Melvin H Friedman，Jonathan Nichols．Introduction to infrared and electro-optical systems［M］．ARTECH HOUSE，1999:163－262．

［4］ WU Pengfei，F(xiàn)ANG Shuai，XU Qingshan，et al．Restoration of blurred image based on atmospheric MTF［J］．Journal of Atmospheric and Environmental Optics，2011，6(3):196－202．(in Chinese)武鵬飛，方帥，徐青山，等．基于大氣MTF的退化圖像復(fù)原［J］．大氣與環(huán)境光學(xué)學(xué)報(bào)．2011，6(3):196－202．

［5］ YAN Juncheng，WANG Defei，QI Wenzong．Research on the electro-optic performance of infrared detectors influenced by the atmosphere turbulence［J］．Laser＆Infrared，2009，39(2):162－165．(in Chinese)閆俊成．大氣湍流對(duì)紅外探測(cè)器光電性能影響的研究［J］．激光與紅外，2009，39(2):162－165．

［6］ Gaskill，J，Linear Systems．Fourier transforms and optics［M］．New York:Wiley，1978:305－307．

［7］ GUO Yue，LIU Xiaopeng，YANG Hua．MTF analysis for CCD［J］．Spacecraft Recovery＆Remote Sensing，2004，25(3):25－28．(in Chinese)郭悅，劉曉鵬，楊樺．CCD的調(diào)制傳遞函(MTF)分析［J］．航天返回與遙感，2004，25(3):25－28．

［8］ REN Hang，ZHANG Tao．Effect of carrier diffusion on modulation transfer function of CCD［J］．Journal of Applied Optics，2009，30(2):257－262．(in Chinese)任航，張濤．載流子擴(kuò)散對(duì)CCD調(diào)制傳遞函數(shù)的影響［J］．應(yīng)用光學(xué)，2009，30(2):257－262．

［9］ HU Fangming．The technical research on modeling and performance evaluation of Electro-optical imaging system［D］．Xi'an:Xidian University，2005．(in Chinese)胡方明．光電成像系統(tǒng)建模及性能評(píng)估技術(shù)研究［D］．西安:西安電子科技大學(xué)，2005．

［10］Vollmerhausen，R．Impact of display modulation transfer function on the quality of sampled imagery［C］．SPIE Proceeding of Aerospace/Defense Sensing and Controls，Orlando FL，Apr．1996．

［11］LIU Huitong，WANG Qi，ZHENG Zhiwei，et al．Analysis and calculation of the nonlinear response in infrared focal plane arrays［J］．J．Infrared Millim．Waves，2001，20(4):253－258．(in Chinese)劉會(huì)通，王琪，鄭志偉，張明，紅外焦平面陣列非線性響應(yīng)的分析和計(jì)算［J］．紅外與毫米波學(xué)報(bào)，2001，20(4):253－258．

［12］Gerald C．Holst．Testing and evaluation of infrared imaging systems［M］．Second Editon．America:JCD Publishing and SPIE Optical Engineering Press，1998．171－249．

［13］ZHANG Zhenzhong，LI Anan．Physical effect modeling of electro-optical imaging system［J］．Journal of Xianyang Normal University，2011，26(2):27－31．(in Chinese)張振中．光電成像系統(tǒng)物理效應(yīng)建模［J］．咸陽(yáng)師范學(xué)院學(xué)報(bào)，2011，26(2):27－31．