《数字图像处理(实验部分)》实验5_空域滤波

1、数字图像处理(实验部分)教案实验五：空域滤波1.实验目的1 .掌握MATLAB的基本操作。2 . 了解数字图像处理在 MATLAB中的基本处理过程。3 .学习空域滤波的原理，观察算法处理结果2 .实验设备2. 1. PC兼容机一台；操作系统为 WindowsWindowsXP 。3. 2.数字图像处理开发环境：MATLAB软件3 .实验原理中值滤波：采用二维中值滤波函数medfilt2对受椒盐噪声干扰的图像滤波空域滤波4 .实验步骤.1打开MATLAB开发环境.2点击MATLAB窗口上File菜单，选择New-M File,在弹出的Edit编辑器内输入如下 程序：clear;close all

2、 ;I0=imread( 'abao.bmp');I0=rgb2gray(I0);I=imnoise(I0, 'salt & pepper' ,0.02);% 叠加密度为 0.02 的椒盐噪声。I_Filter1=medfilt2(I,3 3);烟口大小为 3X 3的中值滤波I2=I;I3=I;x=1 1 1;1 1 1;1 1 1;%3*"滑模板y=-1 -2 -1;0 0 0;1 2 1;%3x骑化模板 h=y;n m=size(I);for i=2:1:n-1for j=2:1:m-1sumx=0;sumy=0;for p=1:1:3fo

3、r q=1:1:3h(p,q)=(1/9)*(I(i+(p-2),j+(q-2);sumx=sumx+h(p,q)*x(p,q);%平滑h(p,q)=(1/4)*(I(i+(p-2),j+(q-2);sumy=sumy+h(p,q)*y(p,q);%锐化endendI2(i,j)=sumx;if sumy<1I3(i,j)=255-sumy;%去除灰度为0 的黑色背景elseI3(i,j)=sumy;endendendI4=I;z=1 1 1 1 1 1 1;1 111111;1 111111;1 111111;1 111111;1 111111;1 111111;%7x坪滑模板k=z;

4、n m=size(I);for i=4:1:n-3for j=4:1:m-3sumz=0;for p=1:1:7for q=1:1:7k(p,q)=(1/49)*(I(i+(p-4),j+(q-4);sumz=sumz+k(p,q)*z(p,q);endendI4(i,j)=sumz;endendfigure(1);subplot(2,3,1),imshow(I0);title(' 原始图 ' );subplot(2,3,2),imshow(I);title(' 加噪声图 ' );subplot(2,3,3),imshow(I_Filter1);title(subplot(2,3,4),imshow(I2);title(subplot(2,3,5),imshow(I3);title(subplot(2,3,6),imshow(I4);title('MATLAB自带3x3中值滤波图');'3x3平滑图');'3x3锐化图');'7x7平滑图');Run按钮，程序会自动运行，并显示出结果。.3将该程序保存，并点击工具栏中.4观察处理结果5 .实验结果滤波后图像中的噪声被抑制原始图加噪声图MATLAB自带3x3中值滤波图3x3平滑图3x3锐化图7x7平滑图6