Coherency Matrix
Understanding the Coherency Matrix
Polarization is a fundamental property of electromagnetic waves, and its complete description requires more than just stating "horizontal" or "circular." Real-world RF signals are often partially polarized: a combination of a deterministic (polarized) component and a random (unpolarized) component. This arises from multipath propagation, rough surface scattering, thermal emission, and mixing of independent sources. The coherency matrix captures this complexity in a compact 2×2 matrix that separates the polarized and unpolarized contributions through its eigenvalue decomposition.
The matrix framework connects to multiple representation systems used in different RF engineering disciplines. Antenna engineers use Jones vectors and matrices for fully polarized signals through deterministic systems. Radiometer engineers use Stokes parameters for power-based measurements of partially polarized thermal emission. Polarimetric radar engineers use scattering matrices and target decomposition theorems. The coherency matrix unifies these: Jones vectors are its rank-1 case, Stokes parameters are linear combinations of its elements, and target decomposition operates on its 3×3 generalization.
Coherency Matrix Equations
J = 〈E EH〉 = [[〈|Ex|²〉, 〈ExEy*〉], [〈EyEx*〉, 〈|Ey|²〉]]
Degree of Polarization:
p = √(1 - 4 det(J) / tr(J)²) = (λ1 - λ2) / (λ1 + λ2)
Stokes from Coherency:
S = [J11+J22, J11-J22, 2Re(J12), 2Im(J12)]
Where λ1 ≥ λ2 ≥ 0 are eigenvalues. Fully polarized: det(J) = 0, p = 1. Unpolarized: J = (P/2)I, p = 0. Partially polarized: 0 < p < 1. Total power = tr(J) = S0.
Polarization State Examples
| State | Coherency Matrix | p | Stokes | Source |
|---|---|---|---|---|
| Horizontal linear | [[1,0],[0,0]] | 1 | [1,1,0,0] | Dipole antenna |
| RHCP | ½[[1,-j],[j,1]] | 1 | [1,0,0,1] | Helix antenna |
| Unpolarized | ½[[1,0],[0,1]] | 0 | [1,0,0,0] | Thermal emission |
| 50% polarized H | [[0.75,0],[0,0.25]] | 0.5 | [1,0.5,0,0] | Ground reflection |
| 45° linear | ½[[1,1],[1,1]] | 1 | [1,0,1,0] | Rotated dipole |
Frequently Asked Questions
How does it describe partial polarization?
Fully polarized: rank 1 (det = 0, one zero eigenvalue), p = 1. Unpolarized: J ∝ I (equal eigenvalues), p = 0. Partially polarized: eigenvalue decomposition uniquely separates polarized + unpolarized components. p = (λ1-λ2)/(λ1+λ2). Multipath, rough scattering, thermal emission produce 0 < p < 1.
How is it used in polarimetric radar?
Scattering matrix S → Pauli vector k → 3×3 coherency matrix T = 〈kkH〉. Cloude-Pottier decomposition: eigenvalues/eigenvectors classify scatterers (surface, volume, double-bounce). Applications: SAR land cover mapping, crop monitoring, forest biomass, disaster assessment.
Relationship to Stokes parameters?
Linear transformation: S0=J11+J22, S1=J11-J22, S2=2Re(J12), S3=2Im(J12). Equivalent information. Stokes: measurable with power detectors (radiometry). Coherency: preserves complex field (radar, interferometry). p = √(S1²+S2²+S3²)/S0 = √(1-4det(J)/tr(J)²).