lee.m

SAR图像三种滤波方法： 1、gamma_MAP滤波； 2、均方差滤波： 3、Lee滤波

function im_ret = Lee( im, w_size, out_size, NL, bound, flag )
%LEE Filter for SAR Speckle reduction
%   im：SAR影像
%   w_size：窗口大小，奇数 3 5 7 9等
%   out_size：输出图像的大小，可以是'same'或者'full'
%   NL：number of Looks，等效视数
%   bound：边缘扩展模式，分为'symmetric'，'replicate','circular'。
%   flag：计算时是否包含中心像素，0为不包含，1为包含。 
%
%   bound和flag主要是给im_mean_var用的 :)
%
%   im_ret：滤波结果
%
%   barbara_lee =Lee( barbara_noisy, 5, 'same', 20, 'symmetric', 0);


im = double(im);

[im_Mean, im_Var] = im_mean_var(im, w_size, out_size, bound, flag);

im_Std = sqrt(im_Var);

Ci = im_Std ./ im_Mean; %is the image variation coefficient

Cu = sqrt(1/NL); % is the estimated noise variation coefficient
tmp = 2*log(Cu) - 2.*log(Ci+0.1);
tmp = exp(tmp);
W = 1 - tmp;

im_ret = im .* W + im_Mean .* (1 - W);