saint.m

含有saint-marc滤波算法的程序

% a=zeros(240,320);
% for k=255:264
%     p=imread([num2str(k),'.jpg']);
%     p=im2double(p);
%     a=a+p;
% end
% a=a/10;  %10帧像素平均值
% imshow(a);
f=imread('255.jpg');
f=im2double(f);
[gx,gy]=gradient(f);
k=4;
for i=1:240
    for j=1:320
        d(i,j)=gx(i,j)^2+gy(i,j)^2;
        w(i,j)=exp(-d(i,j)/(2*k^2));
    end
end
nt=zeros(240,320);
ft=zeros(240,320);
for i=2:240-1
    for j=2:320-1
        for m=-1:1
            for n=-1:1
                nt(i,j)=w(i+m,j+n)+nt(i,j);
                ft(i,j)=f(i+m,j+n)+ft(i,j);
            end
        end
        f(i,j)=ft(i,j)/nt(i,j);
    end
end
figure,imshow(f)