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

Interprétation des matrices de Mueller partiellement dépolarisantes 3×3 mesurées

Informations fondamentales

  • Identifiant de l'article : 2510.06305
  • Titre : Interpretation of measured 3x3 partial depolarizing Mueller matrices
  • Auteurs : José J. Gil, Ignacio San José, Mónica Canabal-Carbia, Juan Campos, Angel Lizana, Irene Estévez
  • Classification : physics.optics
  • Date de publication : ArXiv 7 octobre 2025
  • Lien de l'article : https://arxiv.org/abs/2510.06305

Résumé

La mesure de polarisation de Mueller est une technique puissante avec des applications étendues en astronomie, télédétection, analyse de matériaux avancés et imagerie biomédicale. Cependant, les contraintes instrumentales limitent souvent les mesures à des matrices de Mueller incomplètes, restreintes à la sous-matrice 3×3 de leur coin supérieur gauche. Le simple remplissage des entrées manquantes par des zéros pour former une matrice 4×4 peut produire des résultats physiquement incohérents, même pour les systèmes non-dépolarisants. Pour résoudre ce problème, cet article propose une procédure systématique complétant les matrices de Mueller 3×3 mesurées en matrices 4×4 physiquement cohérentes. La méthode s'appuie sur la forme de matrice de covariance et sélectionne, parmi une infinité de complétions acceptables, celle possédant la pureté de polarisation maximale. Ce critère garantit que la matrice synthétisée correspond au modèle le moins aléatoire (le plus déterministe) compatible avec la mesure.

Contexte de recherche et motivation

Importance du problème

Le problème fondamental que cette recherche résout est de reconstruire une matrice 4×4 complète physiquement cohérente à partir de données de mesure de matrice de Mueller incomplète. Ce problème revêt une importance majeure en polarimétrie car :

  1. Limitations instrumentales pratiques : De nombreux instruments de mesure de polarisation pratiques ne peuvent mesurer que des matrices de Mueller partielles en raison de contraintes techniques, particulièrement la sous-matrice 3×3 du coin supérieur gauche
  2. Besoins d'applications étendues : Dans les applications telles que les endoscopes médicaux, les caméras polarisantes à plan focal, les polarimètres à réfraction conique, etc., l'information complète de la matrice de Mueller 4×4 est essentielle pour la caractérisation précise des matériaux
  3. Exigences de cohérence physique : La simple méthode de remplissage par zéros produit des résultats violant les lois physiques, même pour les systèmes non-dépolarisants

Limitations des méthodes existantes

  • Méthode de remplissage par zéros : Le remplissage direct des éléments manquants par des zéros conduit à des matrices physiquement irréalisables
  • Méthodes d'apprentissage automatique : Bien que certaines recherches tentent d'utiliser des algorithmes d'apprentissage automatique, celles-ci nécessitent une grande quantité d'informations d'échantillons préalables et manquent de fondement physique
  • Méthodes pour cas spécifiques : Les méthodes existantes ciblent principalement les matrices de Mueller pures ou les cas de symétrie spécifique, manquant d'universalité

Contributions principales

  1. Proposition d'une procédure de complétion systématique basée sur la forme de matrice de covariance : Exploitation de la relation mathématique entre la matrice de Mueller et sa matrice de covariance hermitienne associée
  2. Établissement d'un critère de pureté de polarisation maximale : Sélection de la solution la plus déterministe (la moins aléatoire) parmi une infinité de schémas de complétion possibles
  3. Développement d'un algorithme de reconstruction universel : Applicable à toute donnée de mesure de polarisation Mueller partielle 3×3
  4. Fourniture de conditions de contrainte physiquement raisonnables : Garantissant que la matrice reconstruite satisfait les conditions de covariance (semi-définition positive)

Détails de la méthode

Définition de la tâche

Entrée : Sous-matrice Mueller 3×3 mesurée A (partie du coin supérieur gauche) Sortie : Matrice de Mueller 4×4 complète physiquement cohérente M Conditions de contrainte :

  1. La sous-matrice 3×3 du coin supérieur gauche de la matrice reconstruite doit correspondre aux données mesurées
  2. La matrice complète doit satisfaire quatre conditions de covariance (la matrice de covariance associée est semi-définie positive)
  3. Parmi toutes les matrices possibles satisfaisant les conditions ci-dessus, sélectionner celle possédant la pureté de polarisation maximale

Fondements théoriques

Représentation en matrice de covariance de la matrice de Mueller

La matrice de Mueller M peut être exprimée sous forme de blocs : M=(m00DTPm)\mathbf{M} = \begin{pmatrix} m_{00} & \mathbf{D}^T \\ \mathbf{P} & \mathbf{m} \end{pmatrix}

Où :

  • m00m_{00} : Coefficient d'atténuation d'intensité moyenne
  • D\mathbf{D} et P\mathbf{P} : Vecteurs d'atténuation bidirectionnelle et de polarisation
  • m\mathbf{m} : Sous-matrice 3×3

Indicateur de pureté de polarisation

La pureté de polarisation (indice de dépolarisation) est définie comme : PΔ=tr(MTM)m0023m002=D2+P2+PS23P_\Delta = \frac{\text{tr}(\mathbf{M}^T\mathbf{M}) - m_{00}^2}{3m_{00}^2} = \frac{D^2 + P^2 + P_S^2}{3}

PSP_S est l'indice de dimensionnalité de polarisation.

Matrice de covariance H

La matrice de Mueller M et sa matrice de covariance hermitienne 4×4 associée H entretiennent une correspondance biunivoque. La partie réelle de H est entièrement déterminée par les données de mesure 3×3 connues et l'élément m33m_{33} inconnu.

Flux algorithmique

Première étape : Détermination de m33m_{33}

  1. Balayage numérique des valeurs possibles xx de m33m_{33}
  2. Pour chaque valeur xx, construction de la ReH(x)\text{Re}\mathbf{H}(x) correspondante
  3. Vérification que toutes les valeurs propres de ReH(x)\text{Re}\mathbf{H}(x) sont non-négatives
  4. Parmi les valeurs xx satisfaisant les conditions, sélection de celle maximisant la pureté de polarisation PΔ(x)P_\Delta(x)

Deuxième étape : Construction de la matrice de covariance complète

  1. Diagonalisation de ReH\text{Re}\mathbf{H} : ReH=QLQT\text{Re}\mathbf{H} = \mathbf{Q}\mathbf{L}\mathbf{Q}^T
  2. Construction de la matrice K\mathbf{K} forçant la réalisation de la condition bicomposante (rang(H\mathbf{H}) = 2)
  3. Synthèse de la matrice hermitienne complète H\mathbf{H} par transformation unitaire
  4. Calcul de la matrice de Mueller complète en utilisant la correspondance entre H et M

Points d'innovation technique

  1. Méthode de matrice de covariance : Utilisation systématique pour la première fois de la représentation en matrice de covariance de la matrice de Mueller pour résoudre le problème de mesure incomplète
  2. Critère de pureté maximale : Proposition d'un critère de sélection physiquement raisonnable garantissant que le résultat de reconstruction correspond au modèle le plus déterministe
  3. Contrainte bicomposante : Limitation du système à un modèle bicomposante, maintenant la cohérence physique tout en maximisant l'information disponible
  4. Universalité : La méthode s'applique à toute donnée de mesure Mueller partielle 3×3, sans nécessiter d'informations d'échantillons préalables

Configuration expérimentale

Conception des cas de test

L'article utilise des matrices de Mueller construites théoriquement pour la vérification plutôt que des ensembles de données expérimentales réelles. L'avantage de cette approche est de pouvoir contrôler précisément les propriétés de la matrice originale, validant ainsi l'efficacité de l'algorithme de reconstruction.

Indicateurs d'évaluation

  1. Analyse des valeurs propres : Comparaison des valeurs propres de la matrice de covariance originale et reconstruite
  2. Indicateurs de pureté de polarisation : Analyse de l'impact du processus de reconstruction sur l'IPP (indice de pureté de polarisation)
  3. Cohérence physique : Vérification que la matrice reconstruite satisfait les conditions de covariance

Résultats expérimentaux

Résultats principaux

Cas 1 : Système de dépolarisation à quatre composantes

Matrice originale :

undefined