数据图像处理实验

课程 数字图像处理 日期 2011 年10 月 8 日 班级 信计0811 学号 2008151010 姓名 郑筱娴 成绩

实验三 图像增强

一、实验目的与要求

1. 掌握MATLAB语言中的空间域单点增强、灰度变换、直方图修正 2. 掌握MATLAB语言中的空域滤波增强、频域滤波增强

二、实验内容

1. 针对图像“pout.tif”，将图像的每个像素加50，得到图像B1；将图像的每个

像素乘1.2，得到图像B2；将原图像进行灰度变换，得到图像B3；将原图像进行直方图修正，得到图像B4。将变换后的图像在一个窗口内显示出来。 2. 设有RGB图像A“peppers.bmp”，用rgb2gray函数将A转换为灰度图像B，

对图像B叠加密度为0.04的椒盐噪声，得到图像B0。（1）采用二维中值滤波函数medfilt2对受椒盐噪声干扰的图像B0进行滤波，窗口大小分别为3×3、5×5、7×7，得到的图像分别记为B1、B2、B3。（2）采用邻域平均法，窗口大小取为3×3对受椒盐噪声干扰的图像B0进行滤波，得到图像B4。将B、B0、B1、B2、B3、B4在一个窗口内显示出来。

3. 采用以下两种高通滤波器冲激响应阵列对图像blood.tif进行高通滤波。

111121 H1191 H2252 121111三、实验结果

1

A=imread('image/pout.tif'); imshow(A); figure;

subplot(2,2,1) title('d'); B1=A+50; imshow(B1);

title('每个像素加50得到B1') subplot(2,2,2) B2=A*1.2; imshow(B2);

title('每个像素乘1.2得到B2')

subplot(2,2,3)

B3=imadjust(A,[0.15 0.9], [0 1]); imshow(B3);

title('灰度变换得到B3') subplot(2,2,4)

B4=adapthisteq(A); imshow(B4);

title('直方图修正得到B4')

2

clc;clear;

A=imread('image/peppers.bmp'); imshow(A);figure; subplot(2,3,1) B=rgb2gray(A); imshow(B);

title('转换为灰度图像B')

subplot(2,3,2)

B0=imnoise(B,'salt & pepper',0.04); imshow(B0);

title('添加椒盐噪声得B0') subplot(2,3,3)

B1=medfilt2(B0,[3 3]); imshow(B1);

title(' 3×3二维中值滤波得B1') subplot(2,3,4)

B2=medfilt2(B0,[5 5]); imshow(B2);

title(' 5×5二维中值滤波得B2') subplot(2,3,5)

B3=medfilt2(B0,[7 7]); imshow(B3);

title(' 7×7二维中值滤波得B3')

subplot(2,3,6)

B4=filter2(fspecial('average',3),B0); imshow(uint8(B4));

title(' 3×3邻域平均法滤波得B4')