四次測量型全Stokes參數偏振儀的參數優化方法
【專利說明】四次測量型全Stokes參數偏振儀的參數優化方法 【技術領域】
[0001] 本發明屬于偏振探測領域，涉及一種對高斯噪聲和泊松噪聲同時具有免疫能力的 四次測量型全Stokes參數偏振儀的參數優化方法。 【【背景技術】】
[0002] 物體反射的電磁波中含有重要偏振遙感信息，不僅能用于去除背景噪聲，提供高 對比度的表面、形貌、陰影和粗糙度等信息，還可以用于反演目標的尺寸、濃度等物理化學 特性。偏振成像是一種同時獲取目標空間和偏振信息先進遙感探測技術，對提高目標探測、 識別及分類的效率和精準度具有一定潛力，在軍事偵察、地球資源普查、環境衛生監測、自 然災害預報、大氣探測、天文觀測、機器視覺仿生、生物醫學診斷等諸多領域都將具有重要 的應用價值和前景。雖然偏振成像技術是一項新型的前沿遙感探測技術，但其獨特的遙感 探測優勢已引起國內外重要研宄機構的重視。
[0003] 四次測量型全Stokes參數偏振儀主要由延遲器和線偏振器組合而成，其線性模 型為I = W ? S，其中I為含有4次測量強度的矢量，S為含有4個Stokes參數（Sp Si、S2、 S3)的Stokes矢量，W是與延遲器的延遲量和快軸方向有關的4X4維測量矩陣。一般當線 偏振器固定不動，通過改變4次延遲器的延遲量或快軸方向，測得4個強度，然后利用反演 算法S = W<I得出入射光的Stokes參數。由于測量過程中要遭受高斯噪聲或泊松噪聲的 影響，測得的Stokes參數信息將被噪聲淹沒或惡化。當高斯噪聲占主導地位時，為了使測 得的三個參數（Si、S 2、S3)的噪聲方差分布等量化和最小化，文獻【1】針對單旋轉延遲器和 單固定線偏振器組合的全Stokes參數偏振儀提出了一種四次測量的優化方法；他們通過 最小化測量矩陣的品質因數：等權方差，得到了延遲器的優化延遲量為132°，其快軸方向 相對于線偏振器的主透振方向的夾角依次為±51. 7°和±15. 1°即可。文獻【2】針對雙可 調諧延遲器和單固定線偏振器組合的全Stokes參數偏振儀提出了一種四次測量的優化方 法。當兩個延遲器的快軸方向固定時，通過最小化測量矩陣的條件數，得到優化的延遲量， 實現最優化。文獻【3】針對雙旋轉延遲器和單固定線偏振器組合的全Stokes參數偏振儀 提出了一種四次測量的優化方法。當兩個延遲器的延遲量固定時，通過最大化測量矩陣的 效率，得到優化的快軸方向，實現最優化。上述分時測量的優化結構可以轉化為分振幅型、 分孔徑型或分焦平面型快照式全Stokes偏振儀結構，用以同時測量四種參數對應的強度 信息。文獻【4-7】采用了上述三種方法的衍生函數進行優化設計。由于上述優化方法及其 衍生函數均對測量矩陣的總體效應進行優化，因此得到結構參數僅僅對高斯噪聲具有免疫 力，而對泊松噪聲無免疫力；也就是說，隨著入射光偏振態的改變，測得的三個參數（Sp S2、 S3)的噪聲方差分布將不在保持等量化和最小化。文獻【8】曾報道了一個理想偏振儀模型 的測量矩陣：
[0005] 利用該測量矩陣反演得出的Stokes參數對高斯噪聲和光子噪聲均具有免疫力。 但是該矩陣只是一個理想模型，至今尚沒有真正的全Stokes偏振儀結構能夠產生該測量 矩陣。
[0006] 參考文獻：
[0007] 【1 】D. S. Sabatke，M. R. Descour，E. L. Dereniak，W. C. Sweatt，S. A. Kemme，and G. S. Phipps, ''Optimization ofretardance for a complete Stokes polarimeter, 〃0pt. Lett. 25, 802-804(2000)..
[0008] 【2】J.S.Tyo/ ''Noise equalization in Stokes parameter images obtained by use of variable-retardance polarimeters，〃0pt. Lett. 25, 1198-1200 (2000)?
[0009] 【3】J. Zallat,A. S, and M. P. Stoll，"Optimal configurations for imaging polarimeters:impact of image noise andsystematic errors，〃J. Opt. A:Pure Appl. Opt. 8,807-814(2006).
[0010] 【4】A. Hollstein,T.Ruhtz,J.Fischer,and R. Preusker，"Optimization of system parameters for a complete multispectralpolarimeter，〃Appl. Opt. 48,4767-4773(2009).
[0011] 【5】D. Lara and C. Paterson，''Stokes polarimeter optimization in the presence of shot and Gaussian noise，〃0pt. Express 17,21240-21249 (2009) ?
[0012] 【6】A. Peinado, A. Lizana, J. Vidal, C. lemmi, and J. Campos，''Optimized Stokes polarimeters based on a single twistednematic liquid-crystal device for the minimization of noise propagation，"Appl. Opt. 50,5437-5445(2011).
[0013] 【7 】C. He，J. Chang，Y. Wang，R. Liao, H. He，N. Zeng，and H. Ma，"Linear polarization optimized Stokes polarimeterbased on four-quadrant detector，"Appl. Opt. 54, 4458-4463 (2015) ?
[0014]【8】F. Goudail，"Noise minimization and equalization for Stokes polarimeters in the presence of signal-dependentPoisson shot noise，〃0pt. Lett. 34, 647-649(2009). 【
【發明內容】
】
[0015] 本發明的目的在于提供一種四次測量型全Stokes參數偏振儀的參數優化方法， 通過耦合四次測量型全Stokes參數偏振儀的4X4測量矩陣和理想模型的4X4矩陣的優 化方法，得出延遲器的一系列優化延遲量和快軸方向，實現噪聲免疫。通過本發明的方法， 對于測量過程中存在的高斯噪聲或光子噪聲，均可使測得的四個Stokes參數的噪聲方差 最小化，后三個Stokes參數的噪聲方差等量化。
[0016] 為了實現上述目的，本發明采用如下技術方案：
[0017] 四次測量型全Stokes參數偏振儀的參數優化方法，包括以下步驟：
[0018] 四次測量型全Stokes參數偏振儀的偏振調制元件包含沿入射光向依次設置的 一個或多個延遲器、一個固定的線偏振器，其中一個或多個延遲器用一個4X4Muller矩 陣塢表示，M K的各元素為延遲器的延遲量和延遲器的快軸方向的函數，固定的線偏振器用 另一個4X 4Muller矩陣MP表示，M P的各元素為常數；四次測量型全Stokes參數偏振儀的 4X4Muller矩陣M = MPXMK;通過改變延遲器的延遲量和快軸方向四次，得到四個不同的 4X4Muller矩陣M ;然后提取每個矩陣M的第一行作為系統4X4測量矩陣W的每一行，建 立四次測量型全Stokes參數偏振儀的4X4測量矩陣為公式（1)
[0020] 其中，矩陣W的每一行代表一次測量對應的參數，且每一行后三列元素均是延遲 器的延遲量和延遲器的快軸方向的函數；已知對高斯噪聲和泊松噪聲均具有免疫力的理想 模型的4X4矩陣為公式（2)
[0022] 建立使矩陣W后三列元素逼近矩陣W*后三列對應元素的優化評價函數，得到四次 測量分別對應的延遲器的延遲量和快軸方向。
[0023] 本發明進一步的改進在于：建立使矩陣W后三列元素逼近矩陣W*后三列對應元素 的優化評價函數為（3):
[0025] 其中，i = 1-4代表行數，j = 2-4代表列數；利用最小化算法得到四次測量分別 對應的延遲器的延遲量和快軸方向。
[0026] 本發明進一步的改進在于：四次測量型全Stokes參數偏振儀中的每一次測量都 由沿入射光向依次設置的一個可調諧延遲器（110)、一個固定的水平線偏振器（120)、一個 探測器實現；可調諧延遲器（110)的延遲量為其快軸方向（111)相對于固定的水平線 偏振器（120)的主透振方向（121)的夾角為0 i，且逆時針旋轉為正，順時針旋轉為負；利用 延遲器的4X 4Muller矩陣MK和偏振器的4X 4Muller矩陣M P，計算得到測量矩陣W的具體 形式為
[0028] 利用公式（3)使公式（4)的每一行逼近公式（2)對應的每一行，得到四次測量對 應的一系列優化參數為：
[0029] [ S 1= 102. 2。，9 != 71. 9。]，
[0030] [ 8 2= 102.2° , 0 2= -71.9° ],
[0031] [ 8 3= 142. 1° , 0 4= 34. 95° ],
[0032] [ S 4= 142. 1。，9 4=-34. 95。]。
[0033] 本發明進一步的改進在于：四次測量型全Stokes參數偏振儀中的每一次測量都 由沿入射光向依次設置的一個可調諧延遲器（210)、一個可調諧延遲器（220)、一個固定 的水平線偏振器（230)、一個探測器實現；可調諧延遲器（210)的延遲量為，其快軸方 向（211)相對于固定的水平線偏振器（230)的主透振方向（231)的夾角為0 li;可調諧延 遲器（220)的延遲量為S2i;其快軸方向（221)相對于固定的水平線偏振器（230)的主透 振方向（231)的夾角為0 2i;所有夾角逆時針旋轉為正，順時針旋轉為負；利用延遲器的 4 X 4Mu 11 er矩陣MK和偏振器的4 X 4Mu 11 er矩陣M P，計算得到測量矩陣W的具體形式，利用 公式（3)使測量矩陣W的每一行逼近公式（2)對應的每一行，得到四次測量對應的其中一 系列優化參數為：
[0034][( S n= 102。，S 21= 110. 1。），（ 0 n= 63. 53。，9 21=-7. 09。）]，
[0035][( S 12= 102。，S 22= 110. 1。），（ 0 12=-63. 53。，9 22= 7. 09。）]，
[0036] [( S 13= 63. 8。，S 23= 154. 2° )，（ 9 13= 80. 73° ，9 23= 56. 45° )]，
[0037] [( S 14= 63. 8。，S 24= 154. 2。），（ 9 14=-80. 73。，9 23=-56. 45。）];
[0038] 四次測量對應的另一系列優化參數為：
[0039] [( S n= 90。，S 21= 90。），（ 9 n= 31. 82°，9 21=-25. 34。）]，
[0040] [( S 12= 90。，S 22= 90。），（ 9 1