中值滤波matlab程序代码.doc

1、%自适应中值滤波的算法 RAMF%RAMF 主要通过以下两步来处理图像。%1.首先确定最大的滤波半径，然后用一个合适的半径 r 对图像进行滤波。计算当前滤波半径像素灰度的 Imin,Imax,Imed，%然后判断 Imed 是否在Imin,Imax中间，如果在则向下进行，否则扩大当前半径 r 继续滤波直到 r 等于最大滤波半径。%2.如果当前处理的像素 img(i,j)在Imin,Imax之间，则输出当前像素，否则输出当前滤波半径中值像素 Imed。clear all;close all;clc;img=rgb2gray(imread(132.jpg);m n=size(img);img=im

2、noise(img,salt %加入椒盐噪声subplot(2,2,1),imshow(img),title(椒盐噪声图);%普通中值滤波 3*3b=medfilt2(img,3,3);subplot(2,2,2),imshow(b),title(3*3 中值滤波);c=medfilt2(img,5,5);subplot(2,2,3),imshow(c),title(5*5 中值滤波);Nmax=10; %确定最大的滤波半径%下面是边界扩展，图像上下左右各增加 Nmax 像素。imgn=zeros(m+2*Nmax+1,n+2*Nmax+1);imgn(Nmax+1:m+Nmax,Nmax+1

3、:n+Nmax)=img;imgn(1:Nmax,Nmax+1:n+Nmax)=img(1:Nmax,1:n); %扩展上边界imgn(1:m+Nmax,n+Nmax+1:n+2*Nmax+1)=imgn(1:m+Nmax,n:n+Nmax); %扩展右边界imgn(m+Nmax+1:m+2*Nmax+1,Nmax+1:n+2*Nmax+1)=imgn(m:m+Nmax,Nmax+1:n+2*Nmax+1); %扩展下边界imgn(1:m+2*Nmax+1,1:Nmax)=imgn(1:m+2*Nmax+1,Nmax+1:2*Nmax); %扩展左边界re=imgn;for i=Nmax+1:

4、m+Nmaxfor j=Nmax+1:n+Nmaxr=1; %初始滤波半径while r=NmaxW=imgn(i-r:i+r,j-r:j+r);W=sort(W);Imin=min(W(:);Imax=max(W(:);Imed=W(uint8(2*r+1)2/2);if IminImed elser=r+1; %否则扩大窗口，继续判断end endif Iminimgn(i,j) else %否则输出邻域中值re(i,j)=Imed;endendend%I=re(Nmax+1:m+Nmax,Nmax+1:n+Nmax);%subplot(2,2,4),imshow(I),title(RAMF 均值滤波);figure;imshow(re(Nmax+1:m+Nmax,Nmax+1:n+Nmax),);