wpf2.m

来自「白化极化滤波器降斑 计算象素的强度 计算协方差矩阵的元素」· M 代码 · 共 31 行

function[y]=wpf(hh,vv,hv)    %%%%%%  白化极化滤波器降斑
[size1,size2]=size(hh);
mean_hh_vv=0;
for kk=1:size1
    for jj=1:size2
        new_hh(kk,jj)=hh(kk,jj)*hh(kk,jj)';  %% 计算各象素的强度
		new_vv(kk,jj)=vv(kk,jj)*vv(kk,jj)';
		new_hv(kk,jj)=hv(kk,jj)*hv(kk,jj)';
		mean_hh_vv=mean_hh_vv+hh(kk,jj)*vv(kk,jj)';
    end
end
mean_hh=mean2(new_hh);% E(|hh|^2)
mean_vv=mean2(new_vv);
mean_hv=mean2(new_hv);
mean_hh_vv= mean_hh_vv/size1/size2;
abuxi=mean_hv/mean_hh;   %%%%  计算协方差矩阵的元素
gamma=mean_vv/mean_hh;
rou=mean_hh_vv/sqrt(mean_hh*mean_vv);
rou2=1-(abs(rou))^2;
xishu1=mean_hh*rou2';
xishu2=xishu1*gamma;
xishu3=mean_hh*abuxi;
xishu4=mean_hh*rou2*sqrt(gamma);
for kk=1:size1
    for jj=1:size2
        jiao=angle(hh(kk,jj))-angle(vv(kk,jj))-angle(rou);
        y(kk,jj)=new_hh(kk,jj)/xishu1+new_vv(kk,jj)/xishu2+new_hv(kk,jj)/xishu3-2*abs(rou)*abs(hh(kk,jj))*abs(vv(kk,jj))/xishu4*cos(jiao);
  %     y(kk,jj)=new_hh(kk,jj)/xishu1+new_vv(kk,jj)/xishu2+new_hv(kk,jj)/xishu3-2*real(rou*hh(kk,jj)'*(vv(kk,jj)))/xishu1/sqrt(gamma);
    end
end