srad滤波的matlab算法别人的上传

% Implementation of SRAD Filter

% usage S=srad（IT）
% I- noisy image
% T - threshold （greater than 0）.
function S=srad（IT）
tic
[x y]=size（I）;
I=double（I）;
Ic=double（I）;
delta_t = 0.08;
t=1;
eps=0.00000000001;  %%防止出现像素为零的点
for t=1:T     
    qt=exp（-t*.2）;
    [IxIy] = gradient（Ic）;   %%梯度变换 
    di=sqrt（Ix.^2+Iy.^2）;
    di2=del2（Ic）;  %%del2 离散拉普