全极化SAR的H/alpha、anistropy分解
来源:互联网 发布:centos 中文输入法 编辑:程序博客网 时间:2024/04/28 06:32
之前说过,可以利用polSARpro_4.0对原始的SAR数据进行处理,利用软件可以得到极化特征:散射角,熵和反熵,那么如何用程序求出这几个分量呢?以下是实现分解的程序。
clear all;
close all;
%得到相干矩阵T
row = 900;
col = 1024;
%path = 'E:\研究生阶段\PolSAR\实验仿真\数据\原始数据\T3Lee金门';
tic;
fid = fopen('T11.bin','rb');
T11 = fread(fid,[col,row],'float').'; fclose(fid);
fid = fopen('T22.bin','rb');
T22 = fread(fid,[col,row],'float').'; fclose(fid);
fid = fopen('T33.bin','rb');
T33 = fread(fid,[col,row],'float')'; fclose(fid);
fid = fopen('T12_real.bin','rb');
T12_real = fread(fid,[col,row],'float')'; fclose(fid);
fid = fopen('T12_imag.bin','rb');
T12_imag = fread(fid,[col,row],'float')'; fclose(fid);
T12 = T12_real + i * T12_imag;
fid = fopen('T13_real.bin','rb');
T13_real = fread(fid,[col,row],'float')'; fclose(fid);
fid = fopen('T13_real.bin','rb');
T13_imag = fread(fid,[col,row],'float')'; fclose(fid);
T13 = T13_real + i * T13_imag;
fid = fopen('T23_real.bin','rb');
T23_real = fread(fid,[col,row],'float')'; fclose(fid);
fid = fopen('T23_imag.bin','rb');
T23_imag = fread(fid,[col,row],'float')'; fclose(fid);
T23 = T23_real + i * T23_imag;
toc;
tic;
for ii = 1:row
for jj = 1:col
T = [T11(ii,jj) T12(ii,jj) T13(ii,jj);T12(ii,jj)' T22(ii,jj) T23(ii,jj);...
T13(ii,jj)' T23(ii,jj)' T33(ii,jj)];
[VT,DT] = eig(T);
P1 = DT(1,1)/trace(DT);
P2 = DT(2,2)/trace(DT);
P3 = DT(3,3)/trace(DT);
alpha1 = acos(abs(VT(1,1)));
alpha2 = acos(abs(VT(1,2)));
alpha3 = acos(abs(VT(1,3)));
H(ii,jj) = -(P1*log(P1) + P2*log(P2) +P3*log(P3))/log(3);
alpha(ii,jj) = P1*alpha1 + P2*alpha2 + P3*alpha3;
A(ii,jj) = (DT(2,2) - DT(1,1))/(DT(2,2) + DT(1,1));
end
end
toc;
tic;
figure(1);imshow(H);
colormap(jet);colorbar;
title('Polarimetric Scattering Entropy');
figure(2);imshow(alpha);
colormap(jet);colorbar;title('alpha angle');
figure(3);imshow(abs(A));
colormap(jet);colorbar;title('A');
toc;
- 全极化SAR的H/alpha、anistropy分解
- 全极化SAR图像freeman分解
- 基于SVM的极化 SAR 图像分类算法思路
- 关于极化SAR的一些实验结论(不断更新)
- SAR雷达Yamaguchi四分量分解MATLAB代码(注释部分为极化角优化)
- 极化 SAR 图像分类流程
- 关于几种经典的非监督极化SAR分类算法的个人总结
- (原创)极化SAR精致Lee滤波matlab…
- 极化SAR一些常用知识总结(不断更新中)
- 极化码的原理
- 关于UIImage+Alpha.h导入,项目不识别的bug
- 分解全路径的文件名,字符串
- 碳纤维复合材料加固钢筋混凝土的极化作用
- 极化码生成矩阵的构造
- SAR的软件
- linux的sar命令
- 全搜索的H.264编码
- 电磁波的极化中一个易混淆的概念
- 看了netfilter/iptables内核原理,合理的配置实现简单firewall
- mac安装MongoDB
- MethodHandle性能测试
- 剑指offer之数据结构
- [java] 5、URLEncoder 与 URLDecoder 简单解码编码
- 全极化SAR的H/alpha、anistropy分解
- java 值传递和引入传递
- C语言 字符串处理技巧(sscanf、strstr等用法)
- LeetCode题解:Single Number II
- 內向又不愛交際應酬,「非典型業務員」如何讓業績翻倍?因為這一個關鍵...
- linux内核模块编译和安装--kni module移植的makefile
- Cocos2d-x3.0模版容器详解之一:cocos2d::Vector<T>
- plsql中,sql语句中带有中文的查询条件查询不到数据
- C#数据类型安全转换