3150049420110620.zip_SAR数据_SAR读取_SLC sar_SLC 读取_读取slc数据

%《极化雷达图像增强理论》中的PWF %四通道。可生成处理后的散射矩阵 function [ImageS1,ImageS2,ImageS3,ImageS4]=PWFpro(matrixHHr,matrixHVr,matrixVHr,matrixVVr,matrixHHi,matrixHVi,matrixVHi,matrixVVi) %ImageS是处理后的散射矩阵元素 %参数为四通道极化数据的实部和虚部，共8个数据 [m,n]=size(matrixHHr); %|HH| |HV| |VH| |VV|的平方 modHHsqu=matrixHHr.^2+matrixHHi.^2; modHVsqu=matrixHVr.^2+matrixHVi.^2; modVHsqu=matrixVHr.^2+matrixVHi.^2; modVVsqu=matrixVVr.^2+matrixVVi.^2; %按7*7窗口滑窗 for i=1:m-7 for j=1:n-7 %O,E1,E2,r,p,y参数计算 O=mean2(modHHsqu(i:i+6,j:j+6)); E1=mean2(modHVsqu(i:i+6,j:j+6))/O; E2=mean2(modVHsqu(i:i+6,j:j+6))/O; r=mean2(modVVsqu(i:i+6,j:j+6))/O; complexHH=complex(matrixHHr(i:i+6,j:j+6),matrixHHi(i:i+6,j:j+6)); complexVV=complex(matrixVVr(i:i+6,j:j+6),matrixVVi(i:i+6,j:j+6)); complexHV=complex(matrixHVr(i:i+6,j:j+6),matrixHVi(i:i+6,j:j+6)); complexVH=complex(matrixVHr(i:i+6,j:j+6),matrixVHi(i:i+6,j:j+6)); p=mean2(complexHH.*conj(complexVV)); p=p/sqrt(mean2(modHHsqu(i:i+6,j:j+6))*mean2(modVVsqu(i:i+6,j:j+6))); y=mean2(complexHV.*conj(complexVH)); y=p/sqrt(mean2(modHVsqu(i:i+6,j:j+6))*mean2(modVHsqu(i:i+6,j:j+6))); %求ImageS ImageS1(i:i+6,j:j+6)=complexHH; ImageS2(i:i+6,j:j+6)=complexHV/sqrt(E1); ImageS3(i:i+6,j:j+6)=(1/sqrt(1-abs(y).^2))*(complexVH/sqrt(E2)-conj(y)*complexHV/sqrt(E1)); ImageS4(i:i+6,j:j+6)=(1/sqrt(1-abs(p).^2))*(complexVV/sqrt(r)-conj(p)*complexHH); ImageS1(i:i+6,j:j+6)=ImageS1(i:i+6,j:j+6)/sqrt(O); ImageS2(i:i+6,j:j+6)=ImageS2(i:i+6,j:j+6)/sqrt(O); ImageS3(i:i+6,j:j+6)=ImageS3(i:i+6,j:j+6)/sqrt(O); ImageS4(i:i+6,j:j+6)=ImageS4(i:i+6,j:j+6)/sqrt(O); end end