clear all; %清除所有变量
close all; %关闭所有打开的文件
clc; %清除命令行内容
img = mat2gray(rgb2gray(imread('lena.jpg'))); %读取文件
[m n] = size(img);
img = imnoise(img,'salt & pepper',0.25);%加入椒盐噪声
subplot(1,2,1);imshow(img),title('0.25椒盐噪声图');
Smax = 3;%设定最大模板半径
img1 = zeros(m+2*Smax,n+2*Smax);%扩展图像边界
img1(Smax+1:Smax+m,Smax+1:Smax+n) = img;
img1(1:Smax,Smax+1:Smax+n) = img(1:Smax,1:n);
img1(1:Smax+m,Smax+n+1:2*Smax+n) = img1(1:Smax+m,n+1:n+Smax);
img1(Smax+m+1:2*Smax+m,Smax+1:2*Smax+n) = img1(m+1:Smax+m,Smax+1:2*Smax+n);
img1(1:2*Smax+m,1:Smax) = img1(1:2*Smax+m,Smax+1:2*Smax);
img2 = img1;
for i = Smax+1:Smax+m%RAMF滤波
for
close all; %关闭所有打开的文件
clc; %清除命令行内容
img = mat2gray(rgb2gray(imread('lena.jpg'))); %读取文件
[m n] = size(img);
img = imnoise(img,'salt & pepper',0.25);%加入椒盐噪声
subplot(1,2,1);imshow(img),title('0.25椒盐噪声图');
Smax = 3;%设定最大模板半径
img1 = zeros(m+2*Smax,n+2*Smax);%扩展图像边界
img1(Smax+1:Smax+m,Smax+1:Smax+n) = img;
img1(1:Smax,Smax+1:Smax+n) = img(1:Smax,1:n);
img1(1:Smax+m,Smax+n+1:2*Smax+n) = img1(1:Smax+m,n+1:n+Smax);
img1(Smax+m+1:2*Smax+m,Smax+1:2*Smax+n) = img1(m+1:Smax+m,Smax+1:2*Smax+n);
img1(1:2*Smax+m,1:Smax) = img1(1:2*Smax+m,Smax+1:2*Smax);
img2 = img1;
for i = Smax+1:Smax+m%RAMF滤波
for