2025-11-22T13:34:16.543406

Interpretation of measured 3x3 partial depolarizing Mueller matrices

Gil, Jose, Canabal-Carbia et al.

Mueller polarimetry is a powerful technique with broad applications in astronomy, remote sensing, advanced material analysis, and biomedical imaging. However, instrumental constraints frequently restrict the measurement to an incomplete Mueller matrix limited to its upper-left 3x3 submatrix. Simply padding the missing entries with zeros to form a 4x4 matrix can produce physically inconsistent results, even for nondepolarizing systems. To address this issue, we present a systematic procedure to complete 3x3 measured Mueller matrices into physically consistent 4x4 matrices. The method relies on the covariance matrix formalism and selects, among the infinitely many admissible completions, the one with maximal polarimetric purity. This criterion ensures that the synthesized matrix corresponds to the least random (most deterministic) model compatible with the measurement. The procedure is fully general and can be applied to any 3x3 partial Mueller polarimetric data, providing a reliable and physically grounded reconstruction tool for polarimetric imaging and materials characterization.

academic

測定された3×3部分偏光消光Mueller行列の解釈

基本情報

論文ID: 2510.06305
タイトル: Interpretation of measured 3x3 partial depolarizing Mueller matrices
著者: José J. Gil, Ignacio San José, Mónica Canabal-Carbia, Juan Campos, Angel Lizana, Irene Estévez
分類: physics.optics
発表日時: ArXiv 2025年10月7日
論文リンク: https://arxiv.org/abs/2510.06305

要旨

Mueller偏光測定は、天文学、遠隔探査、先進材料分析および生物医学イメージングにおいて広範な応用を有する強力な技術である。しかし、装置の制約により測定は不完全なMueller行列に限定されることが多く、左上隅の3×3部分行列のみに制限される。欠落エントリを単純にゼロで埋めて4×4行列を形成することは、非偏光消光システムであっても物理的に矛盾した結果をもたらす可能性がある。この問題に対処するため、本論文は測定された3×3 Mueller行列を物理的に一貫した4×4行列に補完するための体系的手順を提案する。本手法は共分散行列形式に依存し、無限に多くの許容可能な補完の中から最大偏光純度を有するものを選択する。この基準は、合成された行列が測定と両立する最も確定的（最もランダムでない）モデルに対応することを保証する。

研究背景と動機

問題の重要性

本研究が対処する中核的な問題は、不完全なMueller行列測定データから物理的に一貫した完全な4×4行列をいかに再構成するかである。この問題は偏光測定学において重要な意義を有する。その理由は以下の通りである：

実際の装置制約：多くの実用的な偏光測定装置は技術的制約により部分Mueller行列のみを測定でき、特に左上隅の3×3部分行列に限定される
広範な応用需要：医学内視鏡、焦点面偏光カメラ、円錐屈折偏光計などの応用において、完全な4×4 Mueller行列情報は正確な材料特性評価に不可欠である
物理的一貫性の要件：単純なゼロ埋め法は物理法則に違反する結果をもたらし、非偏光消光システムであっても同様である

既存手法の限界

ゼロ埋め法：欠落要素を直接ゼロで埋めると、物理的に実現不可能な行列が生じる
機械学習手法：機械学習アルゴリズムの使用を試みた研究もあるが、大量の先験的サンプル情報が必要であり、物理的基礎に欠ける
特定条件下の手法：既存手法は主に純Mueller行列または特定の対称性を有する場合に対応し、汎用性に欠ける

中核的貢献

共分散行列形式に基づく体系的補完手順の提案：Mueller行列とその関連するHermitian共分散行列間の数学的関係を利用
最大偏光純度基準の確立：無限に多くの可能な補完案の中から最も確定的（最もランダムでない）解を選択
汎用的な再構成アルゴリズムの開発：任意の3×3部分Mueller偏光測定データに適用可能
物理的に合理的な制約条件の提供：再構成行列が共分散条件（正半定値性）を満たすことを保証

手法の詳細

タスク定義

入力：測定された3×3 Mueller行列部分行列A（左上隅部分）出力：物理的に一貫した完全な4×4 Mueller行列M 制約条件：

再構成行列の左上隅3×3部分行列は測定データと一致する必要がある
完全な行列は4つの共分散条件を満たす必要がある（関連する共分散行列が正半定値）
上記の条件を満たすすべての可能な行列の中から、偏光純度が最大のものを選択する

理論的基礎

Mueller行列の共分散行列表現

Mueller行列Mはブロック形式で表現できる： $\mathbf{M} = \begin{pmatrix} m_{00} & \mathbf{D}^T \\ \mathbf{P} & \mathbf{m} \end{pmatrix}$

ここで：

$m_{00}$ ：平均強度減衰係数
$\mathbf{D}$ および $\mathbf{P}$ ：二向減衰および偏光ベクトル
$\mathbf{m}$ ：3×3部分行列

偏光純度指標

偏光純度（偏光消光指数）は以下のように定義される： $P_\Delta = \frac{\text{tr}(\mathbf{M}^T\mathbf{M}) - m_{00}^2}{3m_{00}^2} = \frac{D^2 + P^2 + P_S^2}{3}$

ここで $P_S$ は偏光次元指数である。

共分散行列H

Mueller行列Mとその4×4 Hermitian共分散行列Hの間には一対一の対応関係が存在する。Hの実部は既知の3×3測定データと未知の $m_{33}$ 要素によって完全に決定される。

アルゴリズムフロー

ステップ1： $m_{33}$ の決定

$m_{33}$ の可能な値 $x$ を数値的に走査する
各 $x$ 値に対して、対応する $\text{Re}\mathbf{H}(x)$ を構成する
$\text{Re}\mathbf{H}(x)$ の固有値がすべて非負であるかを確認する
条件を満たす $x$ 値の中から、偏光純度 $P_\Delta(x)$ を最大化する値を選択する

ステップ2：完全な共分散行列の構成

$\text{Re}\mathbf{H}$ を対角化する： $\text{Re}\mathbf{H} = \mathbf{Q}\mathbf{L}\mathbf{Q}^T$
行列 $\mathbf{K}$ を構成し、2成分条件を強制実装する（rank( $\mathbf{H}$ ) = 2）
ユニタリ変換を通じて完全なHermitian行列 $\mathbf{H}$ を合成する
HとMの対応関係を利用して完全なMueller行列を計算する

技術的革新点

共分散行列法：Mueller行列の共分散行列表現を初めて体系的に利用して不完全測定問題を解決
最大純度基準：物理的に合理的な選択基準を提案し、再構成結果が最も確定的なモデルに対応することを保証
2成分制約：システムを2成分モデルに限定することで、物理的一貫性を保ちながら利用可能な情報を最大化
汎用性：手法は任意の3×3部分Mueller測定データに適用可能であり、先験的サンプル情報を必要としない

実験設定

テストケースの設計

本論文は実際の実験データセットではなく、理論的に構成されたMueller行列を用いて検証を行った。この手法の利点は、元の行列の性質を正確に制御でき、再構成アルゴリズムの有効性を検証できることである。

評価指標

固有値分析：元の行列と再構成行列の共分散行列固有値を比較
偏光純度指標：再構成プロセスがIPP（偏光純度指標）に与える影響を分析
物理的一貫性：再構成行列が共分散条件を満たすかを検証

実験結果

主要な結果

ケース1：4成分偏光消光システム

元の行列： $\hat{\mathbf{M}} = \begin{pmatrix} 1.000 & 0.185 & 0.097 & 0.168 \\ 0.216 & 0.742 & 0.311 & 0.109 \\ 0.077 & 0.282 & 0.280 & 0.009 \\ 0.063 & 0.274 & 0.053 & 0.360 \end{pmatrix}$

固有値： $\lambda_0 = 0.677, \lambda_1 = 0.281, \lambda_2 = 0.038, \lambda_3 = 0.004$ IPP： $P_1 = 0.396, P_2 = 0.882, P_3 = 0.984$

再構成結果：

$m_{33,\text{rec}} = 0.525$
再構成行列固有値： $\lambda_{0,\text{rec}} = 0.917, \lambda_{1,\text{rec}} = 0.083, \lambda_{2,\text{rec}} = 0.000, \lambda_{3,\text{rec}} = 0.000$
再構成IPP： $P_{1,\text{rec}} = 0.834, P_{2,\text{rec}} = 1, P_{3,\text{rec}} = 1$