溫仁杰，魯嘯天，金偉其，劉 敬，王 霞，郭 宏

?

基于SIP-DSWP系統(tǒng)的透明玻璃面傾角分區(qū)標(biāo)定偏振測(cè)量方法

溫仁杰，魯嘯天，金偉其，劉 敬，王 霞，郭 宏

（北京理工大學(xué) 光電學(xué)院，光電成像技術(shù)與系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室，北京 100081）

由于存在邊緣視場(chǎng)的漸暈效應(yīng)，造成SIP-DSWP同時(shí)偏振成像系統(tǒng)的偏振信息重構(gòu)存在偏差。為此，本文提出了一種通過(guò)分區(qū)儀器矩陣校正視場(chǎng)邊緣效應(yīng)的成像系統(tǒng)校正方法。通過(guò)對(duì)SIP-DSWP同時(shí)偏振成像系統(tǒng)3×3分區(qū)儀器矩陣的標(biāo)定，進(jìn)行了玻璃平板傾角的非接觸檢測(cè)實(shí)驗(yàn)。實(shí)驗(yàn)結(jié)果表明：圖像邊緣視場(chǎng)效應(yīng)和噪聲能夠直接影響偏振成像及偏振信息（偏振度和偏轉(zhuǎn)角）重構(gòu)；分區(qū)儀器矩陣較平均儀器矩陣能夠更準(zhǔn)確地重構(gòu)出目標(biāo)場(chǎng)景的偏振信息，提高對(duì)透明玻璃平板傾角的檢測(cè)精度，為后續(xù)偏振成像方法的研究和定量檢測(cè)應(yīng)用奠定了理論基礎(chǔ)。

同時(shí)偏振成像系統(tǒng)；斯托克斯矢量；偏振信息；分區(qū)儀器矩陣

0 引言

當(dāng)透過(guò)透明異質(zhì)界面成像時(shí)，物空間的成像系統(tǒng)和目標(biāo)場(chǎng)景分別位于不同介質(zhì)中。由于兩種介質(zhì)折射率的不同，介質(zhì)交界面的不規(guī)則面形會(huì)導(dǎo)致成像系統(tǒng)的圖像降質(zhì)和扭曲，機(jī)載水下目標(biāo)成像、水下物體測(cè)量、透過(guò)壓花玻璃成像等都屬于此類(lèi)成像模式。為了實(shí)現(xiàn)對(duì)透明異質(zhì)界面成像的圖像復(fù)原，需要有效的透明物體面形測(cè)量方法。

1982年，RMA. Azzam利用了一塊分光片和兩塊渥拉斯頓棱鏡的組合將測(cè)量光束分為不同的四束光，并用四組探測(cè)器對(duì)四束光進(jìn)行探測(cè)，從而獲得了測(cè)量光束的偏振信息。目前國(guó)內(nèi)外已根據(jù)偏振成像技術(shù)提出一些校正水下降質(zhì)和扭曲圖像的方法[1]。偏振成像技術(shù)利用光電成像器件獲取景物輻射或反射信息進(jìn)行成像，不僅可獲得景物的光強(qiáng)度分布，而且能夠獲得景物偏振度、偏振角、偏振橢率等參數(shù)信息，增加目標(biāo)探測(cè)和場(chǎng)景理解的可用信息[2]，對(duì)透明介質(zhì)面形的測(cè)量展現(xiàn)出技術(shù)可行性，可用于透明介質(zhì)（諸如玻璃、水面波紋）面形的偏振成像測(cè)量。2008年由美國(guó)哥倫比亞大學(xué)地球觀測(cè)站、澳大利亞新南威爾士大學(xué)、美國(guó)馬薩諸塞大學(xué)以及美國(guó)紐約春分公司組成的研究小組提出了一種基于偏振成像方法測(cè)量水面短波紋二維傾角的偏振傾角測(cè)量方法PSS（Polarimetric slope sensing）[3-4]，利用白天無(wú)偏的天空光作為入射光，通過(guò)偏振分光成像系統(tǒng)獲取入射光和反射光偏振角的關(guān)聯(lián)，從而推斷出視場(chǎng)中圖像的瞬時(shí)二維偏振度和偏振角。2009年馬薩諸塞大學(xué)Schultz Howard 等獲得美國(guó)發(fā)明專(zhuān)利“透過(guò)不規(guī)則水面成像的系統(tǒng)和方法”[5]，通過(guò)偏振測(cè)量獲得水面幾何結(jié)構(gòu)，進(jìn)而利用其去除海面起伏時(shí)引起的圖像扭曲。2016年劉敬等研究了考慮探測(cè)器特性的光電偏振成像信息重構(gòu)方法[6]，證明通過(guò)儀器矩陣標(biāo)定，可消除SIP-DSWP（Simultaneous Imaging Polarimetry with Double Separate Wallaston Primes）同時(shí)偏振成像系統(tǒng)[7]探測(cè)器g校正系數(shù)的影響。但由于在整個(gè)偏振成像幅面上偏振光學(xué)系統(tǒng)往往存在漸暈等影響，使得采用視場(chǎng)中心標(biāo)定單一儀器矩陣難以適應(yīng)全畫(huà)幅視場(chǎng)的校正。為此，本文將分析SIP-DSWP同時(shí)偏振成像系統(tǒng)的基礎(chǔ)上，研究采用分區(qū)儀器矩陣的方法，以減小系統(tǒng)校正的誤差，提高透明介質(zhì)的三維面形檢測(cè)的精度。

1 SIP-DSWP同時(shí)偏振成像系統(tǒng)原理

1.1 理想偏振光學(xué)系統(tǒng)的儀器矩陣

光波偏振特性可表示為一個(gè)4×1的斯托克斯[8]矢量：

＝[]

式中：表示總光強(qiáng)；和分別表示線偏正分量；表示圓偏振分量。用穆勒矩陣描述光學(xué)元件對(duì)光波偏振態(tài)的改變：

式中：i和o分別為入射光和出射光的斯托克斯矢量；表示光學(xué)元件/系統(tǒng)的穆勒矩陣。

對(duì)于四通道同時(shí)偏振成像系統(tǒng)，入射光的i經(jīng)過(guò)偏振分光元件到達(dá)同時(shí)偏振成像系統(tǒng)4個(gè)子通道探測(cè)器的光強(qiáng)(1)、(2)、(3)和(4)滿(mǎn)足：

式中：ins為偏振成像系統(tǒng)的儀器矩陣。

若ins可逆，則可求得入射光的斯托克斯矢量i：

為滿(mǎn)足ins可逆條件，需要4個(gè)子通道的偏振分量互不相關(guān)，通常選擇0°、45°、90°和135°四個(gè)檢偏方向。對(duì)于理想偏振成像系統(tǒng)，對(duì)應(yīng)的儀器矩陣為：

于是，依據(jù)斯托克斯矢量計(jì)算得到偏振度和偏振角：

要得到動(dòng)態(tài)目標(biāo)場(chǎng)景的偏振信息，需要滿(mǎn)足：獲得同一時(shí)刻目標(biāo)場(chǎng)景的4幅偏振度不同的圖像以及偏振成像系統(tǒng)的儀器矩陣。

1.2 SIP-DSWP同時(shí)偏振成像系統(tǒng)

偏振成像系統(tǒng)大致分為分時(shí)和同時(shí)成像兩類(lèi)[9]。分時(shí)偏振成像裝置結(jié)構(gòu)簡(jiǎn)單、成本低，但偏振圖像在不同時(shí)間獲取，難以獲得動(dòng)態(tài)目標(biāo)場(chǎng)景的偏振信息。同時(shí)偏振成像在同一時(shí)刻獲得目標(biāo)的四幅不同偏振態(tài)圖像，探測(cè)速度快，可用于快速變化目標(biāo)的偏振探測(cè)。

結(jié)合分振幅和分孔徑同時(shí)偏振成像的特點(diǎn)[10]，我們?cè)O(shè)計(jì)了一套SIP-DSWP同時(shí)偏振成像系統(tǒng)[7]。如圖1，系統(tǒng)通過(guò)分振幅元件將入射光分為傳播方向垂直的兩路光（不調(diào)制光的偏振態(tài)）；然后通過(guò)分孔徑方式，將每一路光分為兩束偏振態(tài)正交的光沿一定夾角出射，分別成像在同一個(gè)探測(cè)器的左、右兩側(cè)，實(shí)現(xiàn)在兩個(gè)探測(cè)器上同時(shí)獲得不相關(guān)的4幅偏振分量圖像，實(shí)現(xiàn)動(dòng)態(tài)場(chǎng)景的偏振成像；景物輻射在鏡頭3的焦平面成像，在鏡頭3的焦平面放置視場(chǎng)光闌，其后放置單色濾光片和鏡頭2，并與2個(gè)鏡頭1分別構(gòu)成雙路中繼光學(xué)系統(tǒng)，將鏡頭2的像分別成像到探測(cè)器1和探測(cè)器2；在鏡頭2和鏡頭1之間依次放置消偏振分光棱鏡、1/2波片、渥拉斯頓棱鏡；鏡頭3的像經(jīng)過(guò)鏡頭2和消偏振分光棱鏡分光，分為透射光路和反射光路，透射光依次經(jīng)過(guò)快軸方位角為0°的1/2波片（平衡兩路光路的光程）、渥拉斯頓棱鏡和鏡頭1后，在探測(cè)器1的寬度方向形成左右兩幅圖像，分別對(duì)應(yīng)0°和90°的偏振圖像；反射光經(jīng)過(guò)快軸方位角為22.5°的1/2波片（相對(duì)于透射光路引入了45°的偏振方位變化）、渥拉斯頓棱鏡和鏡頭1后，由于反射時(shí)的鏡像效應(yīng)，在探測(cè)器2上形成的左右兩幅圖像分別對(duì)應(yīng)135°和45°的偏振圖像。當(dāng)采用外觸發(fā)方式時(shí)，可同時(shí)在兩個(gè)探測(cè)器上獲得同一場(chǎng)景偏振方向相差45°的4種不同偏振態(tài)圖像。

SIP-DSWP同時(shí)偏振成像系統(tǒng)采集到的原始圖像如圖2，經(jīng)過(guò)形態(tài)學(xué)處理、二值化處理、圖像邊緣提取、倒像/鏡像處理、分割配準(zhǔn)等的配準(zhǔn)圖像如圖3所示。

結(jié)合圖1和圖3可以看出：由于成像系統(tǒng)光路較長(zhǎng)，加之光學(xué)元件難免的偏光性能差異，四路偏振光路不僅存在偏振特性差異（可采用儀器矩陣校正），而且在成像視場(chǎng)內(nèi)也存在漸暈等不均勻成像，采用單一儀器矩陣描述全視場(chǎng)偏振成像性能存在困難。

2 SIP-DSWP同時(shí)偏振成像系統(tǒng)的分區(qū)儀器矩陣

2.1 分區(qū)儀器矩陣

根據(jù)E-P定標(biāo)原理[11]，將偏振片放入積分球均勻的輸出端面與偏振成像系統(tǒng)之間的光路中，并在0°～350°范圍內(nèi)每隔10°逐次產(chǎn)生線偏振均勻場(chǎng)景輻射，采集36組實(shí)驗(yàn)設(shè)置下的偏振圖像提供系統(tǒng)儀器矩陣定標(biāo)。注意到成像視場(chǎng)的非均勻性，將定標(biāo)圖像分為3×3九個(gè)子區(qū)域（子區(qū)域大小170×216）)，通過(guò)對(duì)應(yīng)子區(qū)域偏振圖像像素灰度的平均求取儀器矩陣，然后采用如圖4所示的加權(quán)系數(shù)對(duì)儀器矩陣進(jìn)行加權(quán)平均，作為四通道SIP-DSWP同時(shí)偏振成像系統(tǒng)的儀器矩陣0。

圖1 SIP-DSWP同時(shí)偏振成像系統(tǒng)及其成像原理圖

圖2 SIP-DSWP同時(shí)偏振成像系統(tǒng)獲取的原始圖像

圖3 SIP-DSWP同時(shí)偏振成像系統(tǒng)的配準(zhǔn)圖像

如上所述，需要考慮漸暈等不均勻成像的影響，誠(chéng)然，逐點(diǎn)建立儀器矩陣可獲得最佳的校正效果，但因逐個(gè)像素定標(biāo)數(shù)據(jù)過(guò)大，會(huì)造成實(shí)際儀器矩陣參數(shù)儲(chǔ)存和校正計(jì)算量偏大，不利于實(shí)際實(shí)時(shí)的校正處理。為此，我們認(rèn)為可采用分區(qū)儀器矩陣進(jìn)行處理。這里直接采用圖4所示的9個(gè)子區(qū)域進(jìn)行計(jì)算，得到的子區(qū)域儀器矩陣1至9如表1所示?？梢钥闯觯焊鱾€(gè)子區(qū)域的儀器矩陣存在差異，若僅用平均儀器矩陣0計(jì)算整幅圖像的斯托克斯矢量，不利于準(zhǔn)確獲得景物的偏振信息。

2.2 對(duì)透明玻璃面傾角的非接觸測(cè)量實(shí)驗(yàn)

如圖5，測(cè)量實(shí)驗(yàn)是在陰天室外自然光照明條件下進(jìn)行，將K9平板玻璃（＝1.5163）與SIP-DSWP系統(tǒng)水平放置在同一直線上（間隔2m），系統(tǒng)物鏡焦距為35mm，將平板玻璃的傾斜角度與偏振成像系統(tǒng)光軸分別調(diào)整成10°、20°、30°、40°、50°、60°、70°和80°共8個(gè)角度，并同時(shí)采集相應(yīng)的0°、45°、90°、135°四幅玻璃平板的偏振圖像。

圖4 儀器矩陣加權(quán)系數(shù)

圖5 室外偏振成像系統(tǒng)玻璃平板角度檢測(cè)實(shí)驗(yàn)

表1 理想、分區(qū)加權(quán)平均以及分區(qū)儀器矩陣

將采集到8組不同角度的圖像分為9個(gè)子區(qū)域并選取其中心點(diǎn)作為采樣點(diǎn)，分別用平均儀器矩陣與該子區(qū)域儀器矩陣計(jì)算斯托克斯矢量（由于80°玻璃平板無(wú)法充滿(mǎn)整個(gè)視場(chǎng)，特截取對(duì)應(yīng)區(qū)域內(nèi)玻璃板區(qū)域中心作為采樣點(diǎn)），進(jìn)而按式(5)計(jì)算對(duì)應(yīng)采樣點(diǎn)的偏振度、偏振角和玻璃平板傾角。

2.2.1 偏振度

表2給出分別采用平均和分區(qū)儀器矩陣計(jì)算的透明玻璃板偏振度?？梢钥闯觯?）用平均儀器矩陣重構(gòu)的偏振度在分區(qū)5的結(jié)果與理想偏振度最為接近；2）用平均儀器矩陣在各分區(qū)計(jì)算的偏振度均優(yōu)于采用平均儀器矩陣的結(jié)果；3）利用分區(qū)儀器矩陣計(jì)算的偏振度基本在分區(qū)5最接近理想偏振度。這說(shuō)明由于圖像邊緣非均勻性，采用分區(qū)儀器矩陣是必要的，可減小偏振度的計(jì)算誤差。

2.2.2 偏振角

表3給出分別采用平均和分區(qū)儀器矩陣計(jì)算的玻璃板偏振角。同樣，用平均儀器矩陣重構(gòu)的偏振角與理想偏振角基本在分區(qū)5最為接近；用分區(qū)儀器矩陣在各分區(qū)計(jì)算的偏振角均優(yōu)于采用平均儀器矩陣的結(jié)果；分區(qū)儀器矩陣計(jì)算的偏振角最接近理想偏振角的分區(qū)較偏振度的離散更大，但基本在分區(qū)3、4和5，同樣說(shuō)明分區(qū)儀器矩陣是必要的。

2.2.3 入射角

當(dāng)入射光為自然光時(shí)，菲涅爾公式給出反射光偏振度Dolp與入射角、折射角l之間的定量關(guān)系[12]：

式中：為介質(zhì)折射率。入射角與折射角之間滿(mǎn)足折射定律：sin＝×sin。

由圖5可知，光線入射角與玻璃平板的傾角一致。因此，根據(jù)SIP-DSWP系統(tǒng)確定的偏振度Dolp，可以進(jìn)一步確定入射光的入射角即玻璃平板的傾角如表4所示?？梢钥闯觯孩儆善骄鶅x器矩陣計(jì)算的入射角相對(duì)理想入射角有大有小，而由分區(qū)儀器矩陣計(jì)算的入射角均大于理想入射角，且偏差幅度均小于平均儀器矩陣；②平均儀器矩陣計(jì)算的最佳值均在第5區(qū)，分區(qū)儀器矩陣在第5區(qū)最接近理想數(shù)值。

表2 平均儀器矩陣與分區(qū)儀器矩陣分別計(jì)算采樣點(diǎn)的偏振度

表3 平均儀器矩陣與分區(qū)儀器矩陣分別計(jì)算的采樣點(diǎn)偏振角

實(shí)際上，第5區(qū)屬于視場(chǎng)中心，其偏轉(zhuǎn)度、偏振角計(jì)算結(jié)果更為接近理想結(jié)果。對(duì)偏振成像的非理想性及軸外偏振信息及玻面傾角重構(gòu)模型等，都是提高偏振信息及玻璃板的檢測(cè)精度需要進(jìn)一步研究的內(nèi)容。

表4 平均儀器矩陣與分區(qū)儀器矩陣分別計(jì)算的采樣點(diǎn)入射角

3 總結(jié)

基于SIP-DSWP同時(shí)偏振成像系統(tǒng)，研究了在偏振成像系統(tǒng)視場(chǎng)邊緣存在漸暈條件下的分區(qū)儀器矩陣校正方法，利用E-P定標(biāo)法將成像視場(chǎng)分為3×3個(gè)分區(qū)域進(jìn)行分區(qū)域儀器矩陣標(biāo)定，并對(duì)玻璃平板傾角進(jìn)行了檢測(cè)實(shí)驗(yàn)驗(yàn)證，證明能夠獲得較文獻(xiàn)[6]采用9個(gè)子區(qū)域的平均儀器矩陣更準(zhǔn)確的偏振度、偏振角和玻璃平板傾角信息，為后續(xù)偏振成像方法研究和定量檢測(cè)應(yīng)用奠定理論基礎(chǔ)。

同時(shí)偏振成像系統(tǒng)能夠獲得動(dòng)態(tài)場(chǎng)景目標(biāo)的偏振信息，提高復(fù)雜背景下目標(biāo)檢測(cè)率。我們以透明玻璃平板作為檢測(cè)目標(biāo)，已驗(yàn)證偏振成像可實(shí)現(xiàn)對(duì)水面波紋等的成像檢測(cè)，對(duì)發(fā)展新型水下目標(biāo)探測(cè)技術(shù)具有意義。

[1] 劉敬, 金偉其, 王霞, 等. 透明物體面形偏振成像測(cè)量技術(shù)綜述[J]. 紅外技術(shù), 2014, 36(9): 681-687.

Liu Jing, Jin Weiqi, Wang Xia, et al. Review of shape from polarization for transparent object[J]., 2014, 36(9): 681-687.

[2] 劉敬, 夏潤(rùn)秋, 金偉其, 等. 基于斯托克斯矢量的偏振成像儀器及其進(jìn)展[J]. 光學(xué)技術(shù), 2013, 39(1): 56-65.

Liu Jing, Xia Runqiu, Jin Weiqi, et al. Review of imaging polarimetry based on stokes vector[J]., 2013, 39(1): 56-65.

[3] Zappa C J, Banner M L, Schultz H, et al. Retrieval of short ocean wave slope using polarimetric imaging[J]., 2008, 19(5): 120-125.

[4] Pezzaniti J.L., Chenault D., et al. Wave slope measurement using imaging polarimetry[C]//, 2009, 7317: 73170B-13.

[5] Schultz H, Corrada-Emmanuel A. System and method for imaging through an irregular water surface: US, US7630077[P]. 2009.

[6] 劉敬, 金偉其, 王霞, 等. 考慮探測(cè)器特性的光電偏振成像系統(tǒng)偏振信息重構(gòu)方法[J]. 物理學(xué)報(bào), 2016, 65(9): 094201-1.

LIU Jing, JIN Weiqi, WANG Xia, et al. Considering the probe gamma the polarization information reconstruction method[J]., 2016, 65(9): 094201-1.

[7] 劉敬, 金偉其, 王霞, 等. 雙分離渥拉斯頓棱鏡同時(shí)偏振成像系統(tǒng)設(shè)計(jì)[J]. 光學(xué)學(xué)報(bào), 2015, 35(5): 0511001-10.

Liu Jing，Jin Weiqi，Wang Xia，et al. Design of Simultaneous Imaging Polarimetry with double separate Wollaston prism[J]., 2015, 35(5): 0511001-10.

[8] 廖延彪. 偏振光學(xué)[M]. 北京: 科學(xué)出版社, 2003.

Liao Yanbiao.[M]. Beijing: Science Press, 2003.

[9] R. M. Matchko. High-speed imaging chopper polarimetry[J]., 2008, 47(1): 16-18.

[10] J.S.Tyo. Review of passive imaging polarimetry for remote sensing[J]., 2006, 45(22): 5453-5469.

[11] Azzam. Arrangement of four photo-detectors for measuring the state of polarization of light[J]., 1985, 10(7): 309-311.

[12] L.B. Wolff. Polarization vision: a new sensory approach to image understanding[J]., 1997, 15(2): 81-93.

Partition Calibration Polarization Measurement Method of the Transparent Glass Surface Angle Based on SIP-DSWP System

WEN Renjie，LU Xiaotian，JIN Weiqi，LIU Jing，WANG Xia，GUO Hong

(,,100081,)

Because of the vignette effect of the edge field viewing, deflection exists in the reconstruction of polarization information using polarization imaging system. Therefore, in this paper, we proposed an imaging system correction method. The method was realized by filed viewing effect correction using the partition instrument matrix. The non-contact testing experiment of the glass-plate angle was conducted by calibrating the 3×3 partition instrument matrix of the SIP-DSWP synchronous polarization imaging system. The result of the experiment shows that the imaging edge field viewing effect and its noise are likely to affect the reconstruction of polarization imaging and polarization information. Compared to the instrument matrix, the partition instrument matrix can reconstruct the polarization information of the target view more accurately, increasing the test accuracy of the transparent glass-plate angle and subsequently establishing an academic foundation for the study of polarization imaging and quantitative testing application.

polarization imaging system at the same time，the stokes vector，polarization information，partitioning instrument matrix

O436

A

1001-8891(2017)09-0807-07

2016-12-05；

2017-06-29.

溫仁杰（1990-），男，碩士研究生，主要從事偏振成像及其圖像處理方面的研究。E-mail：635225289@qq.com。

金偉其（1961-），男，教授，博士生導(dǎo)師，主要從事光電圖像處理、夜視與紅外技術(shù)、光電檢測(cè)與儀器等方面的研究。E-mail：jinwq@bit.edu.cn。

國(guó)家自然科學(xué)基金面上項(xiàng)目（61575023），總裝預(yù)研基金。