极化SAR S转化为T matlab

极化合成孔径雷达 极化SAR 散射矩阵s转化为t矩阵的matlab代码

 %% S矩阵转化成T矩阵
 clc;clear all;close all;
tic;
%% 读入数据
S_dir=‘S2_OriginalData‘;%存放补偿之后数据的文件名
T3_dir=‘T3_OriginalData‘;%存放分解结果的文件名
%%
W=1;
[S11pt]=freadenvi（strcat（‘.\‘S_dir‘\S11‘））;%.\表示当前目录..\表示上个目录
[S12pt]=freadenvi（strcat（‘.\‘S_dir‘\S12‘））;
[S21pt]=freadenvi（strcat（‘.\‘S_dir‘\S21‘））;
[S22pt]=freadenvi（strcat（‘.\‘S_dir‘\S22‘））;

[T11T12T13T21T22T23T31T32T33]=scattermt3（S11S12S21S22）;
T12_imag=real（T12）;
T12_real=real（T12）;
T13_imag=real（T13）;
T13_real=real（T13）;
T23_imag=real（T23）;
T23_real=real（T23）;

%% 新建文件夹
%dirname=[strcat（‘.\‘Asym_PTD_basedOn_OA_dir）];%新的文件夹名 .\表示这个目录 ..\表示上个目录
dirname=[strcat（‘.\‘T3_dir‘_W_‘num2str（W））];%新的文件夹名 .\表示这个目录 ..\表示上个目录

a=[‘mkdir‘dirname];% 创建命令
system（a）;%创建文件夹
cd（dirname）;%进入新的文件夹

 fid1 = fopen（‘T11.bin‘‘w‘）;
 fwrite（fid1T11‘float‘）;
 fid1 = fopen（‘T22.bin‘‘w‘）;
 fwrite（fid1T22‘float‘）;
  fid1 = fopen（‘T33.bin‘‘w‘）;
 fwrite（fid1T33‘float‘）;

 fid1 = fopen（‘T12_imag.bin‘‘w‘）;
 fwrite（fid1T12_imag‘float‘）;
 fid1 = fopen（‘T13_imag.bin‘‘w‘）;
 fwrite（fid1T13_imag‘float‘）;
  fid1 = fopen（‘T23_imag.bin‘‘w‘）;
 fwrite（fid1T23_imag‘float‘）;
 
  fid1 = fopen（‘T12_real.bin‘‘w‘）;
 fwrite（fid1T12_real‘float‘）;
 fid1 = fopen（‘T13_real.bin‘‘w‘）;
 fwrite（fid1T13_real‘float‘）;
  fid1 = fopen（‘T23_real.bin‘‘w‘）;
 fwrite（fid1T23_real‘float‘）;


%% 复制头文件

copyfile（strcat（‘..\‘S_dir‘\T11.hdr‘）‘.\T11.hdr‘）;%*为通配符 复制头文件
copyfile（strcat（‘..\‘S_dir‘\T11.hdr‘）‘.\T22.hdr‘）;%*为通配符 复制头文件
copyfile（strcat（‘..\‘S_dir‘\T11.hdr‘）‘.\T33.hdr‘）;%*为通配符 复制头文件

copyfile（strcat（‘..\‘S_dir‘\T11.hdr‘）‘.\T12_real.hdr‘）;%*为通配符 复制头文件
copyfile（strcat（‘..\‘S_dir‘\T11.hdr‘）‘.\T12_imag.hdr‘）;%*为通配符 复制头文件
copyfile（strcat（‘..\‘S_dir‘\T11.hdr‘）‘.\T13_real.hdr‘）;%*为通配符 复制头文件
copyfile（strcat（‘..\‘S_dir‘\T11.hdr‘）‘.\T13_imag.hdr‘）;%*为通配符 复制头文件
copyfile（strcat（‘..\‘S_dir‘\T11.hdr‘）‘.\T23_real.hdr‘）;%*为通配符 复制头文件
copyfile（strcat（‘..\‘S_dir‘\T11.hdr‘）‘.\T23_imag.hdr‘）;%*为通配符 复制头文件


%copyfile（strcat（‘..\‘T3_dir‘\T11.hdr‘）‘.\Freeman分解结果.txt‘）;%*为通配符 复制头文件
% copyfile（strcat（‘..\‘T3_OA_dir‘\config.txt‘）‘.\‘）;%复制config文件

cd（‘..\‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-09-25 09:27  convert\
     文件        2412  2017-04-30 20:23  convert\main.m
     文件        1181  2017-04-30 19:41  convert\scattermt3.m