胡勁松 宋志平

摘 要：偏振光譜強(qiáng)度調(diào)制（PSIM）是一種獲取偏振光譜信息的新技術(shù)，測(cè)量數(shù)據(jù)處理是PSIM偏振光譜儀的關(guān)鍵模塊之一。利用IDL程序設(shè)計(jì)語言，完成PSIM偏振光譜儀測(cè)量數(shù)據(jù)處理算法的集成設(shè)計(jì)。從PSIM技術(shù)原理出發(fā)，進(jìn)行PSIM偏振光譜儀測(cè)量數(shù)據(jù)處理算法的數(shù)學(xué)推導(dǎo)，確立PSIM偏振光譜儀數(shù)據(jù)處理算法流程，編程實(shí)現(xiàn)了從光譜儀測(cè)量數(shù)據(jù)文件讀入到待測(cè)光Stokes矢量元素譜及偏振度譜輸出的完整過程，完成了PSIM偏振光譜儀數(shù)據(jù)處理算法的可視化設(shè)計(jì)，使其具有良好的用戶界面。典型待測(cè)光源測(cè)量實(shí)驗(yàn)數(shù)據(jù)處理結(jié)果與理論分析結(jié)果的一致性證明了算法可行性和軟件設(shè)計(jì)的正確性。同時(shí)，程序功能的可視化設(shè)計(jì)也使PSIM偏振光譜儀測(cè)量數(shù)據(jù)處理效率明顯提升。

關(guān)鍵詞：偏振光譜儀;數(shù)據(jù)處理;IDL編程語言;軟件設(shè)計(jì)

DOI：10. 11907/rjdk. 182753

中圖分類號(hào)：TP312文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1672-7800（2019）001-0069-04

Abstract： Polarization spectral intensity modulation（PSIM） is a new technology for obtaining polarization spectral information and its measurement data processing is one of the key modules of? spectropolarimeter. The measurement data processing algorithms integrated design of spectropolarimeter based on polarization spectral intensity modulation （PSIM） technology has been completed by using an Interactive Data Language. Based on PSIM technology principle， the detail mathematical deduction for the spectropolarimeter measurement data processing algorithm was presented， the data processing algorithm process of the spectropolarimeter was established， and the whole procedure of reading data from the spectropolarimeter， data processing and outputting the results was realized by writing a software. The consistency between the experimental data processing results of the typical photometric sources and the theoretical analysis results show that the paper work is successful. Firstly， the visual design of the data processing algorithm of the spectropolarimeter is realized and it has a good user interface. Secondly， the correctness of the software design is verified by the coincidence between the results of the experimental data with the theoretical analysis results. Finally， the processing efficiency of the measurement data of the spectropolarimeter is improved significantly for the visual function designation.

0 引言

由于在天文觀測(cè)、生物醫(yī)學(xué)、大氣遙感監(jiān)測(cè)、軍事等領(lǐng)域皆具有獨(dú)特的應(yīng)用優(yōu)勢(shì)，近年來光的偏振態(tài)測(cè)量成為相關(guān)領(lǐng)域研究熱點(diǎn)[1-5]。光的偏振狀態(tài)可用Stokes矢量完整描述。來自目標(biāo)反射光、輻射光的Stokes矢量元素譜含有與目標(biāo)本征特性相關(guān)的信息。偏振光譜儀通過測(cè)量目標(biāo)反射光、部分輻射光或完整的Stokes矢量元素譜信號(hào)，實(shí)現(xiàn)對(duì)目標(biāo)相關(guān)特性參數(shù)的反演。偏振光譜強(qiáng)度調(diào)制（Polarization Spectral Intensity Modulation，PSIM）是實(shí)現(xiàn)目標(biāo)4個(gè)Stokes矢量元素譜同步測(cè)量的先進(jìn)技術(shù)。與常規(guī)偏振光譜儀相比，基于該技術(shù)的偏振光譜儀結(jié)構(gòu)簡(jiǎn)單緊湊，無轉(zhuǎn)動(dòng)部件，探測(cè)器只需一次曝光，結(jié)合相關(guān)解析算法，即可從測(cè)量結(jié)果中解調(diào)得到待測(cè)光完整的Stokes矢量元素譜數(shù)據(jù)。強(qiáng)度調(diào)制技術(shù)在偏振光譜分辨率、測(cè)量精度、Stokes參量的完整性及測(cè)量系統(tǒng)成本與可靠性等方面，具有其它偏振光譜測(cè)量技術(shù)無法替代的優(yōu)越性[6-9]。但PSIM偏振光譜儀數(shù)據(jù)處理算法相對(duì)較復(fù)雜，數(shù)據(jù)處理算法軟件是PSIM偏振光譜儀的關(guān)鍵模塊之一，其設(shè)計(jì)的正確性和可靠性直接影響儀器性能。IDL（Interactive？Data？Language）編程語言因在運(yùn)算速度、圖像處理及圖形用戶界面等功能方面的比較優(yōu)勢(shì)，近年來在遙感監(jiān)測(cè)、醫(yī)學(xué)影像處理等眾多領(lǐng)域得到廣泛應(yīng)用[10-16]。因此本文基于IDL語言進(jìn)行PSIM偏振光譜儀數(shù)據(jù)處理軟件開發(fā)。

本文基于PSIM偏振光譜儀技術(shù)原理，推導(dǎo)分析PSIM偏振光譜儀數(shù)據(jù)處理算法，建立PSIM偏振光譜儀數(shù)據(jù)處理流程，利用IDL 編程語言實(shí)現(xiàn)數(shù)據(jù)處理算法軟件的集成設(shè)計(jì)，并通過典型光源測(cè)量實(shí)驗(yàn)驗(yàn)證軟件設(shè)計(jì)正確性。

1 PSIM偏振光譜儀技術(shù)原理

儀器由調(diào)制器和光柵光譜儀兩大模塊組成。兩塊延遲器和一塊檢偏器構(gòu)成調(diào)制器模塊，實(shí)現(xiàn)待測(cè)光4個(gè)Stokes矢量元素譜同步調(diào)制并疊加成強(qiáng)度譜輸出。光柵光譜儀記錄調(diào)制器輸出的強(qiáng)度譜，通過傅里葉變換、濾波及解調(diào)制等數(shù)字信號(hào)處理過程，可從該強(qiáng)度譜中解析出待測(cè)光的4個(gè)Stokes矢量元素譜，實(shí)現(xiàn)偏振光譜儀的功能。該過程原理的詳細(xì)數(shù)學(xué)推導(dǎo)如下[18]：

假定來自待測(cè)入射光的光譜偏振信息可用Stokes矢量元素譜完整描述為：

式中[σ]為波數(shù)（[σ=1λ]）。當(dāng)入射光經(jīng)過調(diào)制器模塊后，其偏振特性的改變可由調(diào)制器模塊中光學(xué)件的米勒矩陣級(jí)聯(lián)運(yùn)算得到。當(dāng)調(diào)制器模塊中的高階延遲器和檢偏器按圖1所示的位置關(guān)系安裝時(shí)，可推出調(diào)制器模塊輸出光的強(qiáng)度譜（或稱功率譜）表達(dá)式為：

因此，入射光經(jīng)過調(diào)制器模塊后輸出的[P（σ）]功率譜，可看作為經(jīng)過不同頻率載波信號(hào)調(diào)制的入射Stokes矢量元素譜的線性疊加，當(dāng)忽略延遲器中O光、E光折射率差的色散性時(shí)，載波中心頻率由延遲器厚度決定。光譜儀中的CCD探測(cè)器只需一次曝光，即可記錄調(diào)制器模塊輸出的功率譜。

2 PSIM偏振光譜儀數(shù)據(jù)處理算法流程

由PSIM偏振光譜儀的技術(shù)原理可知，通過對(duì)PSIM偏振光譜儀記錄的待測(cè)光強(qiáng)度譜數(shù)據(jù)[Pσ]進(jìn)行傅里葉逆變換、濾波、傅里葉變換和解調(diào)制等處理，可得到待測(cè)光的4個(gè)Stokes矢量譜，完成待測(cè)光偏振光譜信息解析。具體數(shù)據(jù)處理算法流程如圖2所示。

圖2為PSIM偏振光譜儀數(shù)據(jù)處理算法流程。首先，將光譜儀測(cè)量得到的待測(cè)光強(qiáng)度譜數(shù)據(jù)[Pσ]進(jìn)行傅里葉逆變換（IFFT），因?yàn)樵搹?qiáng)度譜數(shù)據(jù)是已調(diào)制待測(cè)光4個(gè)Stokes矢量元素譜的線性疊加，通過合理設(shè)計(jì)延遲器厚度，傅里葉逆變換后在光程差域中可以分開;再通過數(shù)字濾波，分別提取光程差域中的4個(gè)Stokes矢量，對(duì)濾波提取的各個(gè)量進(jìn)行傅里葉變換（FFT），還原到波數(shù)域;最后進(jìn)行解調(diào)制處理，將還原結(jié)果分別除以各自的解調(diào)系數(shù)，得到待測(cè)光的4個(gè)Stokes矢量元素譜。

3 PSIM偏振光譜儀軟件集成設(shè)計(jì)與驗(yàn)證

基于圖2的算法流程，利用IDL程序設(shè)計(jì)語言，完成PSIM偏振光譜儀數(shù)據(jù)處理算法軟件的集成設(shè)計(jì)，實(shí)現(xiàn)了測(cè)量數(shù)據(jù)文件導(dǎo)入、算法處理及解調(diào)結(jié)果輸出等相關(guān)功能。原始數(shù)據(jù)、中間處理結(jié)果及最終輸出結(jié)果均為可視化。

軟件采用結(jié)構(gòu)化程序設(shè)計(jì)，主程序設(shè)計(jì)思路遵循PSIM偏振光譜儀數(shù)據(jù)處理算法流程，軟件界面設(shè)計(jì)運(yùn)用基于過程的GUI設(shè)計(jì)方法，讀取數(shù)據(jù)選用數(shù)據(jù)文件格式輸入和輸出作為數(shù)據(jù)的處理方法。軟件運(yùn)行界面如圖4所示。

為驗(yàn)證PSIM偏振光譜儀數(shù)據(jù)處理算法軟件集成設(shè)計(jì)的正確性，在實(shí)驗(yàn)室對(duì)水平、垂直線偏振光及橢圓偏振光等典型偏振光源進(jìn)行測(cè)量實(shí)驗(yàn)，并利用本文設(shè)計(jì)的軟件對(duì)測(cè)量實(shí)驗(yàn)數(shù)據(jù)進(jìn)行處理。部分處理結(jié)果如圖5-圖6所示。

圖5（a）和6（a）是PSIM偏振光譜儀測(cè)量得到的水平和垂直線偏振光，理論表明其均僅含有[S0（σ）]和[S1（σ）]兩個(gè)Stokes矢量，由式（2）可知其輸出強(qiáng)度譜數(shù)據(jù)是[S0（σ）]和已調(diào)制的[S1（σ）]線性疊加而成，所以圖5（a）和6（a）結(jié)果與理論分析吻合;圖5（b）和6（b）是測(cè)量光譜的自相關(guān)函數(shù)，理論表明已調(diào)制的Stokes矢量在光程差域上是可以分離的，且已調(diào)制的[S1（σ）]的載波信號(hào)中心頻率對(duì)應(yīng)光程差域[B（σ）D2]處，已調(diào)制的[S0（σ）]的載波信號(hào)中心頻率對(duì)應(yīng)光程差域原點(diǎn)處，所以圖5（b）和圖6（b）的結(jié)果與理論分析相符;圖5（c）和6（c）是測(cè)量數(shù)據(jù)解調(diào)得到的處理結(jié)果。理論表明水平線偏振光有[S0（σ）]和[S1（σ）]兩個(gè)Stokes矢量且相等，不含有[S2（σ）]和[S3（σ）]分量。垂直線偏振光有[S0（σ）]和[S1（σ）]兩個(gè)Stokes矢量，大小相等但符號(hào)相反。在實(shí)驗(yàn)誤差合理范圍內(nèi)，圖5（c）和6（c）的結(jié)果與理論分析一致;圖5（d）和6（d）是水平和垂直線偏振光的偏振度，在有效測(cè)量波段內(nèi)其偏振度近似為1，與理論結(jié)果一致。

4 結(jié)語

PSIM偏振光譜儀的特點(diǎn)是硬件結(jié)構(gòu)簡(jiǎn)單，但算法流程較為復(fù)雜，涉及傅里葉變換、數(shù)字濾波、傅里葉逆變換及解調(diào)制等一系列數(shù)字信號(hào)處理過程，而數(shù)據(jù)處理算法是PSIM偏振光譜儀實(shí)現(xiàn)待測(cè)光偏振光譜信息獲取的關(guān)鍵環(huán)節(jié)。為了更直觀、高效地處理光譜儀測(cè)量得到的待測(cè)光譜數(shù)據(jù)，提取待測(cè)光的偏振光譜信息，本文在充分理解PSIM偏振光譜儀測(cè)量原理的基礎(chǔ)上，確立了PSIM 偏振光譜儀數(shù)據(jù)處理流程。利用IDL程序設(shè)計(jì)語言，實(shí)現(xiàn)了PSIM偏振光譜儀測(cè)量數(shù)據(jù)處理軟件的集成設(shè)計(jì)，并通過大量典型光源測(cè)量實(shí)驗(yàn)數(shù)據(jù)，驗(yàn)證了算法流程及數(shù)據(jù)處理軟件設(shè)計(jì)的正確性。后續(xù)工作重點(diǎn)是優(yōu)化軟件設(shè)計(jì)，進(jìn)一步提高PSIM偏振光譜儀偏振測(cè)量精度。

參考文獻(xiàn)：

[1] KARMAN S B， DIAH S Z M， et al. Bio-inspired polarized skylight-based navigation sensors： a review[J]. Sensors， 2012，12（11）：14232-14261.

[2] 劉曦，朱垣洹，何宏輝，等. 偏振光成像技術(shù)在眼科學(xué)中的應(yīng)用展望[J]. 生命科學(xué)儀器，2017，15（5）：3-8.

[3] KUNNEN B， MACDONALD C， DORONIN A， et al. Application of circularly？polarized？light？for non-invasive diagnosis of cancerous tissues and turbid tissue-like scattering media[J]. Journal of Biophotonics，2015，8（4）：317-323.

[4] ARAI K， NISHIMURA Y. Aerosol parameter estimation through atmospheric polarization observation with the different observation angles [J]. SPIE，2009，7461（6）：746106-9.

[5] IVANOV Y S，SINYAVSKI I I，SOSONKIN M G. Astrospectropolarimeters for weak objects[J]. Journal of Optical Technology，2006，73（12）：879-882.

[6] KAZUHIKO O，TAKAYUKI K. Spectroscopic polarimetry with a channeled spectrum[J]. Optics Letters，1999，24（21）：1475-1477.

[7] DERENIAK E L， HAGEN N A， JOHNSON W R. Imaging spectropolarimetry[J]. SPIE，2003，5074：272-285.

[8] OKABE H，MATOBA K，HAYAKAWA M，et al. New configuration of channeled spectropolarimeter for snapshot polar metric measurement of materials[C]. Proceedings of SPIE，2005：58780H.

[9] 劉強(qiáng)，李雙，龔平. 空間調(diào)制線偏振光譜測(cè)量的誤差影響仿真分析[J]光學(xué)學(xué)報(bào)，2014，34（7）：113-117.

[10] 韓培友. IDL可視化分析與應(yīng)用[M]. 西安：西北工業(yè)大學(xué)出版社，2006.

[11] 花春亮. 基于IDL的資源三號(hào)衛(wèi)星圖像處理分析與質(zhì)量評(píng)價(jià)軟件研究與實(shí)現(xiàn)[D]. 西安：長安大學(xué)，2014.

[12] 劉大超. 基于IDL的林地變化多光譜遙感檢測(cè)系統(tǒng)研建[D]. 北京：北京林業(yè)大學(xué)，2016.

[13] 張艷秋，花春亮. 基于IDL的圖像處理軟件設(shè)計(jì)與實(shí)現(xiàn)[J]. 河南科技，2013（6）：9.

[14] 李昂，陳強(qiáng). 基于IDL試驗(yàn)數(shù)據(jù)分析軟件系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[J]. 上海工程技術(shù)大學(xué)學(xué)報(bào)，2011，25（1）：36-40.

[15] 湯敏. 基于IDL語言的醫(yī)學(xué)圖像處理分析系統(tǒng)的研發(fā)[J]. 生物醫(yī)學(xué)工程學(xué)雜志，2009，26（4）：725-730.

[16] 王克，楊建宇，黃鈺林. 基于IDL語言的雙基SAR成像軟件開發(fā)[J]. 電子設(shè)計(jì)工程，2009，17（7）：107-109.

[17] HAGEN N，LOCKE A M，SABATKE D S，et al. Methods and applications of snapshot spectropolarimetry[J]. SPIE，2004，5432：167-174.

[18] 廖延彪. 偏振光學(xué)[M]. 北京：科學(xué)出版社，2003.

（責(zé)任編輯：江 艷）