基于雙向反射分布函數(shù)的海面紅外偏振特性表征模型

張 弛，吳 鑫，謝 建

(西安電子科技大學(xué) 物理與光電工程學(xué)院，陜西 西安 710071)

1 引 言

隨著CCD技術(shù)和探測(cè)技術(shù)的發(fā)展，針對(duì)目標(biāo)與背景的光學(xué)散射特性的研究被廣泛地應(yīng)用到目標(biāo)識(shí)別[1]、醫(yī)學(xué)檢測(cè)[2]和計(jì)算機(jī)圖形學(xué)[3]等領(lǐng)域中。海面的光學(xué)散射特性是目標(biāo)與背景光學(xué)特性研究的重要內(nèi)容，受到了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注。受到海面背景的影響，目標(biāo)的光強(qiáng)信息常常淹沒(méi)于背景中，降低了海上目標(biāo)探測(cè)與識(shí)別的正確率與效率。例如，海面對(duì)太陽(yáng)或天空輻射的反射會(huì)產(chǎn)生強(qiáng)烈的耀斑，當(dāng)耀斑進(jìn)入光學(xué)探測(cè)系統(tǒng)中時(shí)，產(chǎn)生的強(qiáng)烈雜波會(huì)嚴(yán)重影響成像效果，造成目標(biāo)丟失。與光強(qiáng)探測(cè)相比，基于偏振信息的探測(cè)技術(shù)能夠獲得更多關(guān)于目標(biāo)和背景的輪廓與細(xì)節(jié)信息。研究海面的偏振特性，可為海背景下的目標(biāo)探測(cè)與識(shí)別、目標(biāo)隱身與反隱身以及海上搜救等應(yīng)用提供新的手段。

在紅外波段，海面對(duì)太陽(yáng)和天空散射所產(chǎn)生的太陽(yáng)耀光以及海面的自發(fā)輻射存在著特定的偏振特性。太陽(yáng)耀光和海面自發(fā)輻射的偏振方向和偏振強(qiáng)度有明顯差異。Cooper[4]和Gregoris[5]通過(guò)Cox-Munk 波浪斜率模型分析了太陽(yáng)耀光的偏振特性，并通過(guò)實(shí)驗(yàn)測(cè)量發(fā)現(xiàn)在紅外偏振成像中，在中波紅外波段，探測(cè)器接收的輻射主要來(lái)自海面耀光；在長(zhǎng)波紅外波段，接收的輻射主要來(lái)自海面的自發(fā)輻射。Shaw[6]等的研究表明，海面的偏振發(fā)射率會(huì)隨著觀測(cè)角度、風(fēng)速和波長(zhǎng)等因素的變化而發(fā)生改變，但文中未涉及海面散射光的偏振特性研究。Chang[7]等認(rèn)為海面自發(fā)輻射的偏振度與觀測(cè)角度以及粗糙度等因素有關(guān)，海面越粗糙，偏振度越低。但同樣缺乏海面散射光的偏振特性建模和分析。

以往的研究表明，海面的紅外偏振輻射特性與風(fēng)速、海面的幾何分布情況以及海水的介電特性等多種因素有關(guān)。由于實(shí)時(shí)監(jiān)測(cè)與現(xiàn)場(chǎng)實(shí)驗(yàn)測(cè)量等研究手段會(huì)耗費(fèi)大量的人力物力，為了減少研究成本和縮短研究周期，紅外偏振成像建模與仿真技術(shù)是一種研究海面紅外偏振輻射特性的有效手段。國(guó)外早在2002年就開(kāi)始了對(duì)紅外偏振成像仿真的研究[8-10]，國(guó)內(nèi)在偏振雙向反射函數(shù)(polarized Bidirectional Reflectance Distribution Function，pBRDF)理論[11-12]與材料的偏振特性方面[13-14]取得了一定的進(jìn)展，但對(duì)海面的紅外偏振特性研究較少，尤其是紅外偏振成像仿真領(lǐng)域的研究則鮮有報(bào)道。

在研究海面光學(xué)性質(zhì)的過(guò)程中，通常采用Cox-Munk模型[10,15]對(duì)海面的斜率分布進(jìn)行預(yù)測(cè)。Cox-Munk模型中假設(shè)海面為高斯粗糙面，斜率方差、峰度和偏度等參數(shù)是基于海拔高度為12.5 m處的風(fēng)速計(jì)算得出，所以當(dāng)風(fēng)速過(guò)大或過(guò)小時(shí)，Cox-Munk模型中參數(shù)的不確定度增加[16]。因此，在海面的紅外偏振成像仿真研究中，建立一種符合實(shí)際海情的海面模型具有重要意義。與Cox-Munk模型相比，基于海浪譜和快速傅里葉變換的海面生成方法不僅可以獲得海面的斜率信息，還可以獲得海面的高度信息。目前，海浪譜主要有Neumann譜、PM譜、Jonswap譜和Elfouhaily譜[17]等。其中，Elfouhaily譜是基于實(shí)測(cè)數(shù)據(jù)得到的，并且是一個(gè)全波數(shù)范圍內(nèi)的海浪譜，因此本文選用Elfouhaily譜。在紅外波段，與海面幾何尺寸相比，波長(zhǎng)較小，需將海面分成許多微小面元。在偏振仿真成像過(guò)程中，大多數(shù)文獻(xiàn)[15,18]均將小面元近似看成光滑平面，僅考慮了鏡面反射，忽略了其他方向上的散射，在相同面元數(shù)的情況下，這種方法較為粗略。

本文提出將pBRDF模型和偏振自發(fā)輻射模型結(jié)合起來(lái)，建立了一種紅外偏振輻射模型。首先，在海面光學(xué)特性表征時(shí)，本文將每個(gè)微小面元近似為有一定斜率分布的微粗糙面，這樣更加符合實(shí)際情況，也提高了偏振仿真成像的效率，計(jì)算更加便捷。然后，通過(guò)數(shù)值計(jì)算得到了不同風(fēng)速下偏振度隨觀測(cè)角的變化曲線(xiàn)，對(duì)比分析了不同風(fēng)速對(duì)計(jì)算結(jié)果的影響。最后，建立中波紅外偏振成像的仿真模型，生成了海面和船艦的紅外偏振圖像，并結(jié)合提出的紅外偏振輻射模型對(duì)仿真圖像進(jìn)行了理論分析。該研究結(jié)果可以為海上目標(biāo)探測(cè)、識(shí)別與分類(lèi)等應(yīng)用提供理論依據(jù)和數(shù)據(jù)支撐。

2 海面的生成

本文利用Monte Carlo方法模擬生成海面的高度場(chǎng)。Monte Carlo法的基本思想是對(duì)白噪聲進(jìn)行傅里葉變換，在頻域?qū)Ｃ娴墓β首V進(jìn)行濾波，再進(jìn)行傅里葉逆變換得到海面的高度場(chǎng)[20]。設(shè)海面在x和y方向上的長(zhǎng)度分別為L(zhǎng)x和Ly，等間隔離散點(diǎn)數(shù)分別為Nx和Ny，相鄰點(diǎn)的距離分別為Δx=Lx/Nx和Δy=Ly/Ny。海面上每點(diǎn)的坐標(biāo)可表示為xm=m·Δx和yn=n·Δy，其中m=-Nx/2+1,…,Nx/2且n=-Ny/2+1,…,Ny/2。每點(diǎn)的海面高度z(xm,yn)可表示為[23]:

(1)

(2)

式中：kmk=2πmk/Lx，knk=2πnk/Ly，N(0,1)表示均值為0，方差為1的高斯分布。方向海浪譜Ψ(kmk,knk)的表達(dá)式為[17]：

Ψ(kmk,knk)=

(3)

斜率方差σ2是描述海面粗糙程度的一個(gè)重要參數(shù)，具體計(jì)算公式如下[24]：

(4)

其中：

(5)

(6)

風(fēng)速分別為3，5，10 m/s時(shí)，通過(guò)FFT海面生成方法以及Cox-Munk模型計(jì)算得到的斜率方差σ2以及絕對(duì)誤差δ如表1所示。從表1中可以看出，本文FFT海面生成方法計(jì)算得到的σ2與經(jīng)驗(yàn)值吻合得較好，當(dāng)風(fēng)速為5 m/s時(shí)，絕對(duì)誤差僅為0.000 2。

表1 不同風(fēng)速下的斜率方差σ2和絕對(duì)誤差δ

Tab.1 Slope varianceσ2and absolute errorδwith different wind speeds

Wind speed/(m·s-1)σ2Cox-MunkFFTδ30.018 360.0210.002 6450.028 60.028 80.000 2100.054 20.0520.002 2

3 海面的紅外偏振輻射模型

3.1 基于pBRDF的反射偏振模型

圖1 海面坐標(biāo)系和微面元坐標(biāo)系Fig.1 Coordinate systems of sea surface and single microfacet

BRDF可以表示為散射輻亮度與入射輻照度的比值[25]：

(7)

式中：θi和φi分別是入射光的天頂角和方位角；θr和φr分別是散射光的天頂角和方位角；Lr是散射輻亮度；Ei是入射輻照度。真實(shí)的海面是大量微小面元的集合體，受到風(fēng)速的影響，海表面存在著高低起伏不定的波浪，因此這些微小面元的法線(xiàn)方向各異。菲涅爾定律和斯涅爾定律不能直接應(yīng)用到海面的BRDF計(jì)算中，必須考慮不同微面元對(duì)散射光的不同貢獻(xiàn)。圖1為海面坐標(biāo)系和微面元坐標(biāo)系。其中，n是粗糙面的法線(xiàn)，nμ是微面元的法線(xiàn)，nx是粗糙面本地坐標(biāo)系的x軸；θ是微面元法線(xiàn)與粗糙面表面法線(xiàn)的夾角；β是微面元上法線(xiàn)與入射光的夾角；θ和β可由球面三角幾何得到[21]：

cos 2β=cosθicosθr+sinθisinθrcos(φi-φr),

(8)

(9)

偏振的pBRDF在標(biāo)量的BRDF基礎(chǔ)上發(fā)展而來(lái)，描述了材質(zhì)對(duì)不同偏振方向的光的反射特性，可以表示為[25]:

(10)

式中:Lr和Ei分別是反射輻亮度和入射輻照度的矢量形式。根據(jù)微面元理論，通過(guò)將標(biāo)量BRDF函數(shù)與4×4 Mueller矩陣M作用可以得到F的具體表達(dá)式[25]：

Sh(θi,θr,σ)·M(θi,φi;θr,φr;λ),

(11)

式中Mueller矩陣的各分量Mij可由瓊斯矩陣各分量推出，文獻(xiàn)[21]給出了Mij的表達(dá)式。Sh為遮擋函數(shù)，其表達(dá)式[26]為：

(12)

(13)

其中：θi和θr分別為入射光和反射光的天頂角，erfc(x)是互補(bǔ)誤差函數(shù)?？紤]到紅外波段海水對(duì)輻射的衰減很大，經(jīng)過(guò)多次反射并被探測(cè)器接收到的漫反射光所占的比重極少，故F中未包含漫反射部分。

3.2 偏振輻射模型

來(lái)自海面的紅外輻射由兩部分組成：海面對(duì)太陽(yáng)以及天空的反射輻射Lr和海面的自發(fā)輻射Le，不考慮大氣對(duì)輻射傳輸?shù)挠绊懀綔y(cè)器接收到的紅外輻射Ld表示[15]為：

Ld=Lr+Le,

(14)

Lr=Lsun,r+Lsky,r,

(15)

其中，Lr由海面對(duì)太陽(yáng)的反射輻射Lsun,r和海面對(duì)天空的反射輻射Lsky,r兩部分組成。不同波段下，海面對(duì)太陽(yáng)的反射和對(duì)天空的反射所占的比重不同。在中波波段，主要是對(duì)太陽(yáng)的反射[22]。因此，式(14)可以簡(jiǎn)化為：

Ld=Lsun,r+Le,

(16)

其中[18-19]：

Lsun,r=F(θi,φi;θr,φr;λ)·E(θin,φin,λ)·Sin,

(17)

Le=εsea·LBB(λ,T),

(18)

式中：Lsun,r是海面反射的太陽(yáng)光譜輻亮度，Sin是入射的Stokes矢量，E(θin,φin,λ)為入射到海面上的太陽(yáng)光譜輻照度，LBB(λ,T) 是與海面溫度T相同的黑體輻射亮度，其表達(dá)式如下[27]：

(19)

式中：第一輻射常數(shù)C1=3.741 6×108(W·μm4·m-2) ，第二輻射常數(shù)C2=1.438 8×104(μm·K)。

由于海水對(duì)紅外輻射的衰減很大，可以將海水看成是不透明的，因此根據(jù)基爾霍夫定律，發(fā)射與反射是相關(guān)的。當(dāng)海水處于局部熱平衡狀態(tài)時(shí)，能量達(dá)到守恒狀態(tài)，吸收和發(fā)射相等。海水的發(fā)射率εsea與理想黑體發(fā)射率εbb的關(guān)系[19，28]如下：

εsea=

(20)

式中：dΩi是單位立體角。理想黑體輻射是非偏振的，黑體發(fā)射率可以表示為εbb=[1,0,0,0]T。因此可得出海面的紅外發(fā)射率的表達(dá)式[19,28]為：

(21)

式中：F00，F(xiàn)10，F(xiàn)20和F30分別為矩陣F中對(duì)應(yīng)分量。

在傳統(tǒng)的光線(xiàn)跟蹤方法中，需將二維隨機(jī)粗糙海面劃分為不同的單個(gè)微小面元，并將每個(gè)微面元近似當(dāng)作鏡面。為了有效地描述海面的紋理特征，需要采用更精確的幾何劃分方式，這勢(shì)必會(huì)增大計(jì)算時(shí)間。為了提高計(jì)算效率，本文對(duì)每個(gè)微小面元采用pBRDF模型進(jìn)行計(jì)算，將單個(gè)面元散射貢獻(xiàn)進(jìn)行疊加，得出總的接收輻射為：

(22)

式中：A為生成的二維隨機(jī)海面的面積，M和N分別為在x軸和y軸上的三角面元個(gè)數(shù)。

4 仿真結(jié)果

由式(14)可知，在紅外波段下，探測(cè)器接收到的偏振信息受反射的輻射和自發(fā)輻射共同影響。當(dāng)入射方位和探測(cè)方位確定后，可分別計(jì)算出特定入射幾何和觀測(cè)幾何條件下的反射輻射亮度向量Lsun,r和自發(fā)輻射亮度向量Le，進(jìn)一步可計(jì)算海面自發(fā)輻射的偏振度DOLPe和海面反射輻射的偏振度DOLPr。為了驗(yàn)證理論模型中自發(fā)輻射模型和pBRDF模型的有效性，分別計(jì)算了海面的DOLPe和DOLPr，結(jié)果如圖2～圖4所示。

4.1 偏振輻射模型分析

圖2為通過(guò)pBRDF模型和鏡面反射模型[6]計(jì)算出來(lái)的海面自發(fā)輻射DOLPe隨波長(zhǎng)λ(3.0～4.8 μm)的變化曲線(xiàn)。其中，風(fēng)速為10 m/s，微小面元上的斜率方差σ2=0.052，觀測(cè)天頂角分別為10°，45°，70°。虛線(xiàn)為文獻(xiàn)[6]中的不同觀測(cè)天頂角下的自發(fā)輻射DOLPe隨波長(zhǎng)變化的曲線(xiàn)，實(shí)線(xiàn)為本文通過(guò)pBRDF模型計(jì)算得到的不同觀測(cè)天頂角下自發(fā)輻射DOLPe隨波長(zhǎng)變化的曲線(xiàn)。

圖2 不同觀測(cè)天頂角下DOLPe隨波長(zhǎng)變化的關(guān)系Fig.2 DOLPe versus wavelength with different observation zenith angles

由圖2可知，相同觀測(cè)天頂角下鏡面反射和pBRDF模型計(jì)算得出的DOLPe具有相同的變化趨勢(shì)，但是數(shù)值上略有差異，并且觀測(cè)天頂角越大，差異越大。產(chǎn)生這種現(xiàn)象的原因是：由于鏡面反射模型僅在垂直和水平方向產(chǎn)生偏振，而在pBRDF模型中Mueller矩陣的作用下，海面的自發(fā)輻射在垂直和水平方向以及45°和135°方向同時(shí)具有偏振分量。在光譜段3.0～4.8 μm，海面自發(fā)輻射在垂直和水平方向線(xiàn)偏振光強(qiáng)的差異高于45°和135°方向線(xiàn)偏振光強(qiáng)的差異，從而導(dǎo)致基于pBRDF模型得到的海面自發(fā)輻射DOLPe在數(shù)值上低于鏡面反射模型。

圖3為風(fēng)速分別為3，5，10 m/s下自發(fā)輻射DOLPe隨觀測(cè)天頂角變化的關(guān)系，仿真波長(zhǎng)為4 μm，微小面元上的斜率方差σ2分別為0.021，0.028 8，0.052。由圖3可知，通過(guò)pBRDF模型仿真得到的DOLPe與文獻(xiàn)[6]中結(jié)果較為接近，兩者的DOLPe都隨著觀測(cè)角的增大而增大。在0°～30°時(shí)，自發(fā)輻射的DOLPe較低，與風(fēng)速的相關(guān)性較低；當(dāng)觀測(cè)天頂角大于30°時(shí)，DOLPe隨著風(fēng)速的增大而減小。并且在數(shù)值上，通過(guò)鏡面反射模型得到的DOLPe高于由pBRDF模型得到的DOLPe。這是由于風(fēng)速增大時(shí)，海面粗糙度增大，對(duì)自發(fā)輻射的起偏效果減弱，海面為鏡面時(shí)，偏振效果最明顯。

圖3 不同風(fēng)速下DOLPe隨觀測(cè)天頂角變化的曲線(xiàn)Fig.3 DOLPe versus observation zenith angle with different wind speeds

圖4 不同入射天頂角下DOLPr隨觀測(cè)天頂角變化的曲線(xiàn)Fig.4 DOLPr versus observation zenith angle with different incident zenith angles

圖4是入射天頂角分別為10°，30°，45°，70°時(shí)，通過(guò)pBRDF模型計(jì)算得到的DOLPr隨觀測(cè)天頂角變化的曲線(xiàn)。觀測(cè)方向在入射平面內(nèi)(即|φr-φi|=180°)，風(fēng)速為10 m/s，微小面元上的斜率方差σ2=0.052。

如圖4所示，當(dāng)入射天頂角分別為30°，45°，70°時(shí)，DOLPr峰值分別出現(xiàn)在觀測(cè)天頂角為76.5°，61.5°，37.5°處，此時(shí)反射光達(dá)到了完全偏振；而入射天頂角為10°時(shí)，DOLPr峰值處的觀測(cè)天頂角接近90°，但是在峰值處達(dá)不到百分百偏振。光在無(wú)風(fēng)海面上的反射可以看作是鏡面反射，只有當(dāng)入射天頂角為海水的布儒斯特角(約為53.5°)時(shí)，反射光在θr=53.5°處達(dá)到完全偏振。然而當(dāng)風(fēng)速為10 m/s時(shí)，由于風(fēng)力的影響，海面可以看作是由傾斜的具有一定粗糙度的微粗糙面元組成的粗糙面，所以入射光在每個(gè)微面元上朝著不同方向散射。盡管每個(gè)微面元的傾斜角度不同，入射光在不同微面元上的入射天頂角不同，但整體上還是遵循入射天頂角和觀測(cè)天頂角之和約等于兩倍布儒斯特角。

4.2 紅外海面偏振成像仿真結(jié)果

基于本文提出的紅外偏振輻射模型，利用Modtran軟件計(jì)算太陽(yáng)的輻照度，對(duì)海上探測(cè)器接收的紅外輻射進(jìn)行仿真，并結(jié)合基于海浪譜生成的三維隨機(jī)海面的高度及斜率信息，模擬生成海面紅外偏振圖像，如圖5～圖7所示。仿真參數(shù)如表2所示，仿真平臺(tái)為MATLAB，海面的粗糙度根據(jù)風(fēng)速以及式(4)計(jì)算得到。太陽(yáng)以30°的天頂角以及0°的方位角入射，經(jīng)過(guò)海面各個(gè)微小面元反射至位于(200,200,500)處的探測(cè)器中，通過(guò)每個(gè)面元與探測(cè)器的相對(duì)位置關(guān)系即可計(jì)算出反射天頂角和方位角。此處，方位角為沿著x軸正方向逆時(shí)針?lè)较蛐D(zhuǎn)的角度。

表2 紅外海面偏振成像仿真參數(shù)

Tab.2 Simulation parameters for infrared polarization imaging of sea surface

參 數(shù)設(shè) 置海面區(qū)域128 m×128 m微面元斜率方差0.052風(fēng)速10 m/s風(fēng)向0°海水溫度298.15 K仿真波長(zhǎng)4 μm太陽(yáng)天頂角30°太陽(yáng)方位角0°探測(cè)器位置(200,200,500)艦船折射率6.77艦船粗糙度1.134仿真平臺(tái)MATLAB

圖5 海面自發(fā)輻射偏振仿真圖像Fig.5 Simulated polarimetric infrared images of sea surface without Sun

圖6 包含反射輻射和自發(fā)輻射的海面偏振仿真圖像Fig.6 Simulated polarimetric infrared images of sea surface with Sun

圖5為海面自發(fā)輻射偏振仿真圖像。由于海面自發(fā)輻射的圓偏振分量占比極少，所以S3圖像不予顯示。從圖中可以看出海面自發(fā)輻射的S0圖像較為平滑，由于S0圖像是非偏振的紅外強(qiáng)度圖像，整體對(duì)比度較低，由于海面高低起伏所造成的方向非偏振發(fā)射率差異而產(chǎn)生的光強(qiáng)分布變化不太明顯。偏振分量S1圖像和S2圖像相對(duì)于S0圖像擁有更強(qiáng)的對(duì)比度，可以看到海面更多的細(xì)節(jié)，并且S1圖像比S2圖像的整體灰度更亮，說(shuō)明在當(dāng)前觀測(cè)方向下，海面自發(fā)輻射的水平和垂直線(xiàn)偏振光強(qiáng)的差異高于45°和135°線(xiàn)偏振光強(qiáng)的差異。此外，DOLP圖像的對(duì)比度和細(xì)節(jié)最為明顯，表明利用偏振探測(cè)獲取海面自發(fā)輻射的偏振度圖可以提高圖像的質(zhì)量，能夠?yàn)樽R(shí)別海上目標(biāo)提供依據(jù)。

圖6為包含太陽(yáng)反射和自發(fā)輻射的海面紅外偏振仿真圖像。圖6(a)中，太陽(yáng)耀光現(xiàn)象比較明顯。產(chǎn)生這種現(xiàn)象的原因是在當(dāng)前太陽(yáng)入射幾何條件下，探測(cè)器從坐標(biāo)為(200, 200, 500)處朝向海面中心探測(cè)時(shí)，太陽(yáng)光的入射方向和探測(cè)器的接收方向關(guān)于部分面片的法向量對(duì)稱(chēng)，形成鏡面反射，此方向反射率較高，從而產(chǎn)生了太陽(yáng)耀光。在進(jìn)行海上目標(biāo)檢測(cè)時(shí)，太陽(yáng)耀光會(huì)干擾海上目標(biāo)的檢測(cè)與識(shí)別。由于遮擋與陰影函數(shù)的影響，從某些面片上反射的偏振光會(huì)被附近的面元遮擋，探測(cè)器接收到此面元的光強(qiáng)貢獻(xiàn)為0，因此這些面元在圖像上表現(xiàn)為黑色。偏振圖像S1和S2在一定程度上增強(qiáng)了圖像的對(duì)比度，但是仍然無(wú)法消除太陽(yáng)耀光的影響。從圖6(d)中可以看出，DOLP圖像有效濾除了太陽(yáng)耀光，使圖像變得更加均勻。通過(guò)對(duì)圖6中的S0，S1，S2和DOLP圖像進(jìn)行方差、信息熵及對(duì)比度等圖像質(zhì)量評(píng)價(jià)的計(jì)算，如表3所示。偏振圖像的上述3個(gè)評(píng)價(jià)指標(biāo)均高于非偏振圖像S0，其中S1圖像的3個(gè)指標(biāo)最大，且對(duì)比度最高，圖像所包含信息最為豐富，這是由于海面上S1分量即垂直和水平方向線(xiàn)偏振光差異較大，偏振效應(yīng)較為明顯。

表3 圖像質(zhì)量評(píng)價(jià)結(jié)果

為了更加直觀地對(duì)比海面與人造目標(biāo)的偏振特性差異，本文對(duì)海面以及艦船進(jìn)行了紅外偏振成像仿真。其中，海面使用pBRDF和偏振自發(fā)輻射模型進(jìn)行仿真，而艦船使用典型的P-G模型[21]進(jìn)行仿真，仿真結(jié)果如圖7所示。S0為非偏振圖像，圖像的灰度與紅外輻射亮度成正比，從圖中可以看出，艦船溫度高于海面溫度，細(xì)節(jié)信息較不明顯。圖7(b)、7(c)和7(d)為偏振圖像，其灰度亮暗與溫度無(wú)關(guān)，是偏振度的體現(xiàn)?；叶仍礁?，偏振度越大，灰度越低則反之。從偏振圖像中可以看出，與S0相比，艦船和海面的對(duì)比度顯著提升。并且，偏振圖像下，艦船的紋理細(xì)節(jié)信息更加明顯，不同幾何朝向的面元在偏振度上具有明顯的差異，說(shuō)明本文提出的偏振模型可以用于仿真海場(chǎng)景，可為海上目標(biāo)的檢測(cè)、識(shí)別及跟蹤等應(yīng)用提供數(shù)據(jù)支持。

圖7 包含反射輻射和自發(fā)輻射的海面目標(biāo)偏振仿真圖像Fig.7 Simulated polarimetric infrared images of sea surface and target with Sun

5 結(jié) 論

本文綜合考慮了海面的自發(fā)輻射與反射效應(yīng)對(duì)探測(cè)器接收輻射的影響，通過(guò)建立pBRDF模型和偏振自發(fā)輻射模型，分別對(duì)海面反射輻射和自發(fā)輻射的線(xiàn)性偏振度DOLP進(jìn)行了分析。結(jié)合海浪譜和快速傅里葉變換對(duì)二維海面進(jìn)行了幾何建模。在此基礎(chǔ)上，將海面劃分成若干個(gè)微面元，結(jié)合每個(gè)面元的位置信息和斜率信息，仿真了海面的自發(fā)輻射偏振圖像和紅外輻射偏振圖像。與傳統(tǒng)的光線(xiàn)跟蹤方法不同，本文中將每個(gè)微面元均當(dāng)作一個(gè)具有粗糙度的粗糙面，而不是近似地當(dāng)作鏡面，從而減少了面元數(shù)的劃分，在不影響仿真精度的情況下有效提高了仿真效率。仿真結(jié)果表明：本文仿真模型與微面元鏡面反射模型計(jì)算得出的DOLPe相接近，但是本文的模型能計(jì)算出自發(fā)輻射在45°和135°方向的線(xiàn)偏振輻射，更加真實(shí)地反映海面自發(fā)輻射的實(shí)際情況。另外，本文對(duì)仿真生成的海面紅外偏振圖像進(jìn)行了質(zhì)量評(píng)價(jià)，并將海面與典型人造目標(biāo)艦船的紅外偏振圖像進(jìn)行了對(duì)比，證明了紅外偏振成像在海上目標(biāo)探測(cè)方面的有效性。