RAMF自适应中值滤波

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