极化SAR分解——Freeman-Durden三分量分解

极化SAR分解——Freeman-Durden三分量分解

Freeman三分量分解可将T3（或者C3）矩阵分解为表面散射、二面角散射、体散射三部分
（以下代码针对T3矩阵）

Freeman三分量分解

• 1998年，Cloude和Pottier提出了一种技术，将基于物理的三分量散射机制模型拟合到极化SAR数据本身，模型中包含的三个散射机制是：1）来自随机定向偶极子的冠层散射；2）一阶Bragg表面散射；3）二次散射,通过模型可以估算三个分量对总反向散射的贡献，还可以估算每种机制对HH、HV和VV反向散射项的贡献。*

三分量分解代码

% Freeman三分量分解(T矩阵)
function [ ps,pd,pv,alpha,beta ] = FD3( T11,T22,T33,T12,T13,T23 )

  % 读取影像行列数
  [line,col] = size(T11);
  
  % 为计算的参数矩阵分配内存
  ps = zeros(line,col);
  pd = zeros(line,col);
  pv = zeros(line,col);
  alpha = zeros(line,col);
  beta = zeros(line,col);
  
  % 计算像元个数
  N = line*col;
  
  parfor jj = 1:N

      % 遍历每个像元的T3矩阵
      T3 = [T11(jj) T12(jj) T13(jj);
            conj(T12(jj)) T22(jj) T23(jj);
            conj(T13(jj)) conj(T23(jj)) T33(jj)]; 
      
       % 判定矩阵是否为有效值
      if (T3(1,1)+T3(2,2)+T3(3,3))<eps
          ps(jj) = 0;
          pd(jj) = 0;
          pv(jj) = 0;
          alpha(jj) = 0;
          beta(jj) = 0;
          continue;
      end
           
     
      % 根据T33可以直接计算得到pv
      pv(jj) = 4*T3(3,3);
      % S相当于T11减去体散射的贡献（即表面散射和二面角散射的T11功率之和）
      S = T3(1,1)-pv(jj)/2;   
      % D相当于T22减去体散射的贡献（即表面散射和二面角散射的T22功率之和）
      D = T3(2,2)-pv(jj)/4; 
      C = T3(1,2);
      
      % CO化简后即为4*Shh*Svv,根据van zyl的算法，CO实部的正负可判断主导散射机制是二次散射还是表面散射
      C0 = S-D;
      
      % 当表面散射占主导时
      if C0>0
         alpha(jj) = 0;
         beta(jj) = conj(C/S);    
         ps(jj) = S + abs(C)^2/S;
         pd(jj) = D - abs(C)^2/S;
      end
          
      % 当二面角散射占主导时
      if C0<0
         beta(jj) = 0;
         alpha(jj) = C/D;
         ps(jj) = S - abs(C)^2/D;
         pd(jj) = D + abs(C)^2/D;
      end
           
  end

end

参考文献

[1] Freeman A, Durden S L . A three-component scattering model for polarimetric SAR data[J]. IEEE Transactions on Geoscience & Remote Sensing, 1998, 36(3):963-973.