%原始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