MatLab 均值滤波、中值滤波、高斯滤波

转载 2011年07月18日 09:23:16

转自:http://blog.csdn.net/yanqingan/article/details/5738026

MatLab自编的均值滤波、中值滤波、高斯滤波 图像处理函数。

%自编的均值滤波函数。x是需要滤波的图像,n是模板大小(即n×n)
function d=avefilt(x,n)   
a(1:n,1:n)=1;   %a即n×n模板,元素全是1
p=size(x);   %输入图像是p×q的,且p>n,q>n
x1=double(x);
x2=x1;
%A(a:b,c:d)表示A矩阵的第a到b行,第c到d列的所有元素
for i=1:p(1)-n+1
    for j=1:p(2)-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=midfilt(x,n)   
p=size(x);   %输入图像是p×q的,且p>n,q>n
x1=double(x);
x2=x1;
for i=1:p(1)-n+1
    for j=1:p(2)-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);

 

 

%此为程序主文件,包含主要功能单元,以及对子函数进行调用
try
%实验步骤一:彩色、灰度变换
h=imread('photo.jpg'); %读入彩色图片
c=rgb2gray(h); %把彩色图片转化成灰度图片,256级
figure,imshow(c),title('原始图象'); %显示原始图象
g=imnoise(c,'gaussian',0.1,0.002);  %加入高斯噪声
figure,imshow(g),title('加入高斯噪声之后的图象');  %显示加入高斯噪声之后的图象

%实验步骤二:用系统预定义滤波器进行均值滤波
n=input('请输入均值滤波器模板大小/n');
A=fspecial('average',n);  %生成系统预定义的3X3滤波器
Y=filter2(A,g)/255;           %用生成的滤波器进行滤波,并归一化
figure,imshow(Y),title('用系统函数进行均值滤波后的结果'); %显示滤波后的图象

%实验步骤三:用自己的编写的函数进行均值滤波
Y2=avefilt(g,n);     %调用自编函数进行均值滤波,n为模板大小
figure,imshow(Y2),title('用自己的编写的函数进行均值滤波之后的结果'); %显示滤波后的图象

%实验步骤四:用Matlab系统函数进行中值滤波
n2=input('请输入中值滤波的模板的大小/n');
Y3=medfilt2(g,[n2 n2]);   %调用系统函数进行中值滤波,n2为模板大小 
figure,imshow(Y3),title('用Matlab系统函数进行中值滤波之后的结果');  %显示滤波后的图象

%实验步骤五:用自己的编写的函数进行中值滤波
Y4=midfilt(g,n2);      %调用自己编写的函数进行中值滤波,
figure,imshow(Y4),title('用自己编写的函数进行中值滤波之后的结果');

%实验步骤六:用matlab系统函数进行高斯滤波
n3=input('请输入高斯滤波器的均值/n');
k=input('请输入高斯滤波器的方差/n');
A2=fspecial('gaussian',k,n3);      %生成高斯序列
Y5=filter2(A2,g)/255;              %用生成的高斯序列进行滤波
figure,imshow(Y5),title('用Matlab函数进行高斯滤波之后的结果');    %显示滤波后的图象

%实验步骤七:用自己编写的函数进行高斯滤波
Y6=gaussfilt(n3,k,g);  %调用自己编写的函数进行高斯滤波,n3为均值,k为方差
figure,imshow(Y6),title('用自编函数进行高斯滤波之后的结果');      %显示滤波后的图象

catch           %捕获异常
    disp(lasterr);     %如果程序有异常,输出
end

 

matlab实现图像滤波——高斯滤波

要求 生成一个(2N+1)×(2N+1)大小的高斯模板H(标准为sigma),然后用此模板对图像进行滤波。不允许使用 fspecial 来产生高斯模板,不允许使用 imfilter、conv2 等函数...
  • ytang_
  • ytang_
  • 2016年10月07日 14:49
  • 14956

matlab自带函数实现高斯滤波

自己编写的高斯滤波已理解,可是matlab自带的高斯滤波函数却没有用过,这里记录一下。 %matlab自带函数实现灰度图高斯滤波 clear close all img=imread('lena....

[图像]高斯滤波(Matlab实现)

简述原理 .... 核心提示 在Matlab中高斯滤波非常方便,主要涉及到下面两个函数: 函数: fspecial 函数: imfilter 代码实现 clear all; clc; %-----...

【信号】用matlab实现一维信号的高斯滤波

1.从正态分布说起高斯滤波使用的是高斯函数,即我们熟悉的正态分布的概率密度函数: f(x)=1σ2π−−√exp(−(x−μ)22σ2) f(x)=\frac{1}{\sigma \sqrt{2\p...

matlab自带函数实现高斯滤波(gauss filter)快速算法

将高斯核分为水平方向和垂直方向,先后对图像进行卷积处理,可提高处理速度。 function output=gaussFilter(I,sigma) output=I; ksize=double(ui...

3.高斯滤波、中值滤波MATLAB代码和结果图像

高斯滤波: H=imread('E:\Classical Images\标准测试图片\goldhill.bmp'); I=rgb2gray(H); subplot(2,3,1); imshow(H);...

高斯(核)函数简介

高斯(核)函数简介  1函数的基本概念 所谓径向基函数 (Radial Basis Function 简称 RBF), 就是某种沿径向对称的标量函数。通常定义为空间中任一点x到某一中心xc之间欧氏距...
  • rns521
  • rns521
  • 2012年04月06日 21:56
  • 12822

MATLAB图像平滑与滤波实验

一、实验题目 1.      拍摄的含有正弦波(ref1.bmp, obj1.bmp)和三角波条纹(T20r1.bmp, T20o1.bmp)投影图像 2.      对图像进行滤波去除噪声,噪...

Matlab下各个滤波器传递函数透视图的绘制

本文主要绘制常用滤波器的透视图,包括低通,带阻和带通滤波器。而这三类滤波器又各自包含了理想、巴特沃兹和高斯滤波器。如,低通滤波器就可以分为理想低通滤波器、n阶巴特沃兹滤波器和高斯低通滤波器。 第一:...
  • ljp1919
  • ljp1919
  • 2015年03月31日 20:11
  • 4473

滤波器的设计(低通、带通、高通)的matalb代码

这是之前用到的几个简单的滤波器的程序,滤波器的设计比较简单。可以直接使用matlab里的FDA,导出滤波器系数即可。 再就是参考数字信号处理课本里面三种滤波器设计的方法,由过渡带和所选的窗函数推导出...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:MatLab 均值滤波、中值滤波、高斯滤波
举报原因:
原因补充:

(最多只允许输入30个字)