中值和均值滤波---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是需要滤波的图像,n是模板大小(即n×n)
function d=mid_filter(x,n)
[height, width]=size(x); %输入图像是p×q的,且p>n,q>n
x1=double(x);
x2=x1;
for i=1:height-n+1
for j=1:height-n+1
c=x1(i:i+(n-1),j:j+(n-1)); %取出x1中从(i,j)开始的n行n列元素,即模板(n×n的)
e=c(1,:); %是c矩阵的第一行
for u=2:n
e=[e,c(u,:)]; %将c矩阵变为一个行矩阵
end
mm=median(e); %mm是中值
x2(i+(n-1)/2,j+(n-1)/2)=mm; %将模板各元素的中值赋给模板中心位置的元素
end
end
%未被赋值的元素取原值
d=uint8(x2);


%自编的高斯滤波函数，S是需要滤波的图象，n是均值，k是方差
function d=gaussfilt(k,n,s)
Img = double(s);
n1=floor((n+1)/2);%计算图象中心
for i=1:n
for j=1:n
b(i,j) =exp(-((i-n1)^2+(j-n1)^2)/(4*k))/(4*pi*k);
end
end
%生成高斯序列b。
Img1=conv2(Img,b,'same'); %用生成的高斯序列卷积运算，进行高斯滤波
d=uint8(Img1);

(5) 中值滤波的MATLAB代码
global ImagenUmbral %定义一个全局变量ImagenUmbral
global J %使用全局变量J
A=J %将J赋给A
[m,n] = size(A); %测量图像尺寸参数
B=A; %将A赋给B
pixel_block = zeros(1,9); %开辟一个一行九列的数组，均赋于初值0
for i = 2:m-1 %对A中的各个像素点（不包括边界）
for j = 2:n-2
pixel_block = reshape(A(i-1:i+1,j-1:j+1),9,1); %采用3*3的窗口，将一个像素点及其周围的8各点，依次存入新开辟的数组pixel_block中
sorted_block = sort(pixel_block); %将数组pixel_block排序
block_median = sorted_block(5); %将数

组的中值赋予block_median
B(i,j) = block_median; %将block_median赋予B的相应位置的像素值
end; %end of for j = 2:n-2
end; %end of for i = 2:m-1
B = uint8(B); %将B中所有元素转换成0~255之间的整数
subplot(224) %分割绘图窗口为两行两列，将句柄移到第四个位置
imshow(B) %显示图像B
ImagenUmbral=B %将B赋给全局变量ImagenUmbral
(6) 均值滤波的MATLAB实现
A=J %将J赋给A
[m,n] = size(A); %m为A的行数，n为A的列数
B=A; %将A赋给B
pixel_block = zeros(1,9); %开辟一个一行九列的数组，均赋于初值0
for i = 2:m-1 %对A中的各个像素点（不包括边界）
for j = 2:n-1
pixel_block = reshape(A(i-1:i+1,j-1:j+1),3^2,1);%采用3*3的窗口，将一个像素点及其周围的8各点，一次存入新开辟的数组pixel_block中
mean_block = mean(pixel_block); %求出数组中所有元素的平均值
B(i,j) = mean_block; %将mean_block赋予B的相应位置的像素值
end; %end of for j = 2:n-2
end; %end of for i = 2:m-1
B = uint8(B); %将B中所有元素转换成0~255之间的整数
subplot(224) %分割绘图窗口为两行两列，将句柄移到第四个位置
imshow(B) %显示图像B
ImagenUmbral=B %将B赋给全局变量ImagenUmbral