%《极化雷达图像增强理论》中的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
