用于sar图像预处理的Lee滤波

这是Lee滤波的实现代码，在sar图像的预处理中，Lee滤波用得比较多，对于研究sar图像的人来说具有重要的参考价值...

%Frost filter for speckle noise reduction
%Author : Jeny Rajan
function [ft]=frost（I）
% I is the noisy input image
tic
[x y z]=size（I）;
I=double（I）;
K=1;
N=I;
for i=1:x
    for j=1:y                              
        if （i>1 & i1 & j            mat（1）=I（i-1j）;
            mat（2）=I（i+1j）;
            mat（3）=I（ij-1）;
            mat（4）=I（ij+1）;
            d（1）=sqrt（（i-（i-1））^2）;
            d（2）=sqrt（（i-（i+1））^2）;
            d（3）=sqrt（（j-（j-1））^2）;
            d（4）=sqrt（（j-（j+1））^2）;
            mn=mean（mean（mat））;
            c=mat-mn;
            c2=c.^2;
            c3=c/（c2+.0000001）;
            Cs=0.25*sum（sum（c3））;
            m（1）=exp（-K*Cs*d（1））;
            m（2）=exp（-K*Cs*d（2））;
            m（3）=exp（-K*Cs*d（3））;
            m（4）=exp（-K*Cs*d（4））;
            ms=sum（sum（m））;
            mp=m/ms;
            N（ij）=sum（sum（mp.*mat））;                    
        end
     end
end
toc
ft=uint8（N）;

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

     文件       2844  2006-08-11 16:15  Lee滤波\Lee_Filter\Lee_Filter.mht

     文件        969  2006-11-14 10:34  Lee滤波\sar\sar\frost.m

     文件       2234  2006-11-14 10:35  Lee滤波\sar\sar\lee.m

     文件       2085  2006-11-14 10:35  Lee滤波\sar\sar\kuan.m

     目录          0  2007-03-01 15:10  Lee滤波\sar\sar

     目录          0  2008-12-08 18:38  Lee滤波\Lee_Filter

     目录          0  2008-12-08 18:38  Lee滤波\sar

     目录          0  2008-12-08 18:39  Lee滤波

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

                 8350                    9