全极化SAR图像freeman分解

在全极化SAR图像分类中，我们常常根据提取地物的散射特征来进行分类，其中freeman分解是我们常常用的其中一种非相干目标分解的方法。

1.Freeman分解简单介绍

freeman分解是1998年Freeman和Durden在van Zyl的研究基础上提出的一种非相干矩阵分解方法，其中主要思路是将极化协方差矩阵分解为三种主要的散射机理，即面散射、体散射、二面角散射。其表现形式

其中，fv,fs,fd分别对应着体散射、面角散射和二面角散射分量，对应着的散射功率表示如下：

面散射功率Ps=fs(1+β^2);

二面角散射功率Pd=fd(1+α^2);

体散射功率Pv=8*fv/3;

2.freeman程序(转载时请写明转载出处，请尊重原创)

%%freeman分解，数据可以通过欧航局网站下载，然后通过PolSAR—pro4.o软件进行处理得到相干矩阵T
clc;
clear all;
close all;
row =900;     col =1024;     data = zeros(row,col,9);

%%读入相干矩阵[T]


fIn = fopen( 'T11.bin','r');
data(:,:,1) = fread(fIn,[col,row],'float').';     fclose(fIn);
fIn = fopen( '