freeman 一个极化SAR处理过程中freeman分解的matlab程序

这是一个极化SAR处理过程中freeman分解的matlab程序，希望对初学者又用。

%freeman分解
clc
clear all
close all
%load E:\SanFrancisco\T.mat
load（‘T‘）;
row =900;     
col =1024; 
for i=1:row
    for j=1:col
        t = [T（ij1） T（ij2） T（ij3）;T（ij4） T（ij5） T（ij6）; T（ij7） T（ij8） T（ij9）];  
        fv=（3/2）*t（22）;
        pv=（8/3）*fv;     %求pv
        if real（t（13））>0
            alph=-1;
            beta=1+（t（11）-t（33））/（t（13）+t（33）-2*t（22））;
            be=（t（11）+t（13）-2*t（22））/（t（13）+t（33）-2*t（22））;

%%求ps和pd

%            fs=（t（11）-t（33））/（（abs（beta））^2-1）; 不好
            fs=（t（13）+t（33）-2*t（22））/（1+beta）;
            ps=fs*（1+（abs（beta））^2）;
            fd=t（33）-fs-fv;
            pd=2*fd;
        else
  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1424  2013-03-27 11:06  freeman.m