数字图像处理上机图像平滑锐化.doc

中国地质大学（武汉） 数字图像处理上机实习 （第二专题） 学生姓名： 班 级： 学 号： 指导老师： 实验题目 一、图象灰度变换 0, 显示灰度图象p02-01~p02-06及直方图； 对灰度图象p02-04的灰度范围进行适当展宽； 二、图象平滑 4，对p02-04jy、p02-04gs进行5*5方形窗口的最大均匀性平滑滤波，并比较其效果； 三、图象锐化 4，利用3*3的Krisch算子对p02-04实施图象锐化，二维梯度模板为： H1=， H2= H3=， H4= 实验内容 一 图象灰度变换 要求（1-0） 显示灰度图象p02-01~p02-06及直方图； 1.程序代码 I=imread(D:\matlab2011\work\p02-01.bmp); subplot(621) imshow(I) title(图像) subplot(622),imhist(I,256) title(图像直方图) I=imread(D:\matlab2011\work\p02-02.tif); subplot(623) imshow(I) subplot(624),imhist(I,256) I=imread(D:\matlab2011\work\p02-03.tif); subplot(625) imshow(I) subplot(626),imhist(I,256) I=imread(D:\matlab2011\work\p02-04.tif); subplot(6,2,7) imshow(I) subplot(6,2,8),imhist(I,256) I=imread(D:\matlab2011\work\p02-05.bmp); subplot(6,2,9) imshow(I) subplot(6,2,10),imhist(I,256) I=imread(D:\matlab2011\work\p02-06.tif); subplot(6,2,11) imshow(I) subplot(6,2,12),imhist(I,256) 所用函数： I=imread(‘path’)：函数imread用于读取图片文件中的数据，path为图像文件存放路径。 figure(Name, 图像)：打开一个命名为“图像显示类”图表显示图片。 subplot(m,n,p)：subplot是将多个图画到一个平面上的工具。其中，m表示是图排成m行，n表示图排成n列，也就是整个figure中有n个图是排成一行的，一共m行，如果m=2就是表示2行图。p表示图所在的位置，p=1表示从左到右从上到下的第一个位置。本程序中为6行2列。 imshow(I)：imshow是matlab中显示图像的函数，I为图像中的信息。其调用方式有很多，比如imshow(BW);imshow(I,[low high])等。 imhist(I,256):n为灰度图像灰度级，I为灰度图像，缺省值为256。 2.运行结果 要求（1-4） 对灰度图象p02-04的灰度范围进行适当展宽 1.程序代码 J=imread(D:\matlab2011\work\p02-04.tif); subplot(131) imshow(J,[100,150]); title(灰度范围100-150) J=imread(D:\matlab2011\work\p02-04.tif); subplot(132) imshow(J,[50,200]); title(灰度范围50-200) J=imread(D:\matlab2011\work\p02-04.tif); subplot(133) imshow(J,[0,255]); title(灰度范围0-255) 相关函数： imshow(J,[50,200]):J代表所显示图像的灰度矩阵 ，[50,200]为图像数据的值域。 2.运行结果 二 图象平滑 要求（2-2） 对p02-02jy、p02-02gs进行3*3方形窗口的灰度最相近的K个邻点平均法滤波 1.程序代码 A=imread(D:\matlab2011\work\p02-04.tif); figure(Name, 图像显示); subplot(3,2,1); imshow(A); %显示灰度图像 title(原图象P01-01); J1=imnoise