胡勁松 宋志平

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

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

DOI：10. 11907/rjdk. 182753

中圖分類號：TP312文獻標識碼：A文章編號：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 引言

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

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

1 PSIM偏振光譜儀技術原理

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

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

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

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

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

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

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

3 PSIM偏振光譜儀軟件集成設計與驗證

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

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

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

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

4 結(jié)語

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

參考文獻：

[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] 劉曦，朱垣洹，何宏輝，等. 偏振光成像技術在眼科學中的應用展望[J]. 生命科學儀器，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] 劉強，李雙，龔平. 空間調(diào)制線偏振光譜測量的誤差影響仿真分析[J]光學學報，2014，34（7）：113-117.

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

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

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

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

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

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

[16] 王克，楊建宇，黃鈺林. 基于IDL語言的雙基SAR成像軟件開發(fā)[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] 廖延彪. 偏振光學[M]. 北京：科學出版社，2003.

（責任編輯：江 艷）