医学图像降噪自适应均值滤波AMF 算法

%% 基于自适应中值滤波器对图像去噪处理
clear all;
close all;
clc;
Files=dir（‘/caffe/data/VOCdevkit/VOC20199/JPEGImages/*.jpg‘）;
N=length（Files）;
Names={};
Images={};
for k=1:N
    Names{k}=Files（k）.name;
    Images{k}=imread（[‘/caffe/data/VOCdevkit/VOC20199/JPEGImages/‘ Names{k}]）;
    img=rgb2gray（Images{k}）;       %将原图转成灰度图像
%%img=imread（‘/home/bb/下载/VOC2014/JPEGImages/000001.jpg‘）;
%figure;imshow（img[]）;title（‘原图‘）;     %显示原始图像
[m n]=size（img）;            %mn为图像的行数和列数
img=imnoise（img‘salt & pepper‘0.002）;   %加入20%的椒盐噪声
%figure;imshow（img[]）;title（‘加入20%的椒盐噪声‘）;     %显示加入椒盐噪声后的图像

%% 图像边缘扩展
%为保证边缘的像素点可以被采集到，必须对原图进行像素扩展。
%一般设置的最大滤波窗口为7，所以只需要向上下左右各扩展3个像素即可采集到边缘像素。
Nmax=3;        %确定最大向外扩展为3像素，即最大窗口为7*7
imgn=zeros（m+2*Nmaxn+2*Nmax）;      %新建一个扩展后大小的全0矩阵
imgn（Nmax+1:m+NmaxNmax+1:n+Nmax）=img;  %将原图覆盖在im