读取sar卫星slc数据代码.rar

单复数影像SAR SLC数据的读取代码 matlab

%原始PWF算法
%生成一幅图象
%直接对复数图像进行滤波处理，HH，HV，VH，VV，HV与VH认为是相同的，只考虑三个通道，HH，HV，VV，每个通道对应一个复数矩阵
%HHrHVrVVr分别为三个对应通道的复数据的实部；HHi，HVi，VVi分别为三个对应通道的复数据的虚部
function [ImageOut]=originalPWF（matrixHHrmatrixHVrmatrixVVrmatrixHHimatrixHVimatrixVVi）
%ImageOut has the same size of these matrixit‘s a real Num
%the six matrix are polarization data
[mn]=size（matrixHHr）;
%|HH| |HV| |VV|的平方
modHHsqu=matrixHHr.^2+matrixHHi.^2;
modHVsqu=matrixHVr.^2+matrixHVi.^2;
modVVsqu=matrixVVr.^2+matrixVVi.^2;
%ImageOut=sqrt（modHHsqu+2*modHVsqu+modVVsqu）;
%按7*7窗口滑窗
for i=1:m-7
    for j=1:n-7
        %参数Erp;p还是复数
        E=mean2（modHVsqu（i:i+6j:j+6））/mean2（modHHsqu（i:i+6j:j+6））;
        r=mean2（modVVsqu（i:i+6j:j+6））/mean2（modHHsqu（i:i+6j:j+6））;
        complexHH=complex（matrixHHr（i:i+6j:j+6）matrixHHi（i:i+6j:j+6））;
        complexVV=complex（matrixVVr（i:i+6j:j+6）matrixVVi（i:i+6j:j+6））;
        p=mean2（complexHH.*conj（complexVV））;
        p=p/sqrt（mean2（modHHsqu（i:i+6j:j+6））*mean2（modVVsqu（i:i+6j:j+6）））;
        %ImageOut
        ImageOut（i:i+6j:j+6）=modHHsqu（i:i+6j:j+6）+（1+abs（p）.^2）*modHVsqu（i:i+6j:j+6）/E+modVVsqu（i:i+6j:j+6）/r;
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1277  2011-06-08 16:47  读取sar卫星slc数据代码\20110620\originalPWF.m

     文件       1832  2009-02-24 16:50  读取sar卫星slc数据代码\20110620\PWFpro.m

     目录          0  2020-06-28 17:52  读取sar卫星slc数据代码\20110620

     目录          0  2020-06-28 17:52  读取sar卫星slc数据代码

----------- ---------  ---------- -----  ----

                 3109                    4