A cloude decomposition implementation
Download package and execute either
to install the package or
to install it in development mode.
from cloude import cloude_decompimport numpy as npmueller_matrix = np.array([[1, 0.0707, 0.0348, -0.0060], [0.0480, 0.4099, 0.0077, 0.0650], [0.0162, -0.0184, 0.2243, -0.3580], [0.0021, -0.0465, 0.3571, 0.1783]])decomp_matrix = cloude_decomp(mueller_matrix)
The matrices are scaled by their respective eigenvalues. If you want to use them to fit to an experiment, normalize them by the M_00 matrix element (decomp_matrix[:, 0, 0]).To get the corresponding Jones matrix estimate you may use
from cloude import cloude_decomp_jonesdecomp_matrix = cloude_decomp_jones(mueller_matrix)
Parameter usage is the same forcloude_decomp andcloude_decomp_jones functions.For an output of the eigenvalues use theoutput_eigenvalue parameter
decomp_matrix, ev = cloude_decomp(mueller_matrix, output_eigenvalue=True)
By default the matrix with the highest eigenvalue is returned, which is equivalent to
decomp_matrix = cloude_decomp(mueller_matrix, ev_mask=[1, 0, 0, 0])
You may tune the output as an arbitrary sum of decompositions by the eigenvalue mask
decomp_matrix = cloude_decomp(mueller_matrix, ev_mask=[1, 1, 1, 1])
returns the full sum of all coherent matrices, hencedecomp_matrix == mueller_matrix.