Pauli 协方差分解

详细介绍了极化SAR数据协方差的pauli分解方法含代码

function[]=pauli_cov（c11c22c33c12c13c23str）
%%  使用协方差矩阵数据的 PAULI 分解  （R:HH+VVG:HVB:HH-VV）
if nargin==6
    str=‘pauli‘;
end

B=[1 0 1;1 0 -1;0 1.414 0]/1.414;%%  原协方差矩阵中的HV相已经有了1.414，所以B中不是2
[s1s2]=size（c11）;
pauli_out=zeros（s1s23）;
for ii=1:s1
    for jj=1:s2
        cov=[c11（iijj）c12（iijj）c13（iijj）;c12（iijj）‘c22（iijj）c23（iijj）;c13（iijj）‘c23（iijj）‘c33（iijj）];
        t=abs（B*cov*B‘）;
        if t（11）>1
            pauli_out（iijj1）=0.9999;
        else
            pauli_out（iijj1）=t（11）;
        end
        if t（22）>1
            pauli_out（iijj3）=0.9999;
        else
            pauli_out（iijj3）=t（22）;
        end
        if t（33）>1
            pauli_out（iijj2）=0.9999;
        else
            pauli_out（iijj2）=t（33）;
        end
    end
end
figure;imshow（ （pauli_out））;title（str）;
        

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        907  2006-01-11 14:37  pauli_cov.m

----------- ---------  ---------- -----  ----

                 1125                    2