均值滤波、中值滤波的matlab实现

目前比较经典的图像去噪算法主要有以下三种：
 
均值滤波算法：也称线性滤波，主要思想为邻域平均法，即用几个像素灰度
的平均值来代替每个像素的灰度。有效抑制加性噪声，但容易引起图像模糊，
可以对其进行改进，主要避开对景物边缘的平滑处理。
 

%x是需要滤波的图像,n是模板大小(即n×n)
function d=avg_filter(x,n)   
a(1:n,1:n)=1;   %a即n×n模板,元素全是1
[height, width]=size(x);   %输入图像是hightxwidth的,且hight>n,width>n
x1=double(x);
x2=x1;
for i=1:hight-n+1
    for j=1:width-n+1
        c=x1(i:i+(n-1),j:j+(n-1)).*a; %取出x1中从(i,j)开始的n行n列元素与模板相乘
        s=sum(sum(c));                 %求c矩阵中各元素之和
        x2(i+(n-1)/2,j+(n-1)/2)=s/(n*n); %将与模板运算后的各元素的均值赋给模板中心位置的元素
    end
end
%未被赋值的元素取原值
d=uint8(x2);

中值滤波：基于排序统计理论的一种能有效抑制噪声的非线性平滑滤波信号
处理技术。中值滤波的特点即是首先确定一个以某个像素为中心点的邻域，
一般为方形邻域，也可以为圆形、十字形等等，然后将邻域中各像素的灰度
值排序，取其中间值作为中心像素灰度的新值，这里领域被称为窗口，当窗
口移动时，利用中值滤波可以对图像进行平滑处理。其算法简单，时间复杂
度低，但其对点、线和尖顶多的图像不宜采用中值滤波。很容易自适应化
 

%自编的中值滤波函数。x是需要