一、实验名称:空域图像增强
二、实验目的:掌握Matlab语言图像工具箱中空域图像增强的实现
三、实验要求:
在掌握图像灰度调整、直方图修正和图像锐化的指令基础上,编写程序实现图像的灰度变换,直方图均衡和图像锐化的处理
四、实验仪器和设备:计算机,Matlab软件
五、实验原理:
1、 亮度变换S=T(r)点对点的变换(灰度级对灰度级的变换)matlab函数:imadjust()
亮度变换的基本函数g=imadjust(f,[low in high in],[low out high out],gamma); low in and high in 参数分别指定输入图像需要映射的灰度空间范围,low out 和high out 参数分别指定输出图像所在的灰度范围。GAMMA表示曲线的形状,描述输入输出图像之间的关系。如果GAMMA小于1,则映射的权重趋势向更亮输出,如果GAMMA大于1,则趋向更暗的输出。默认值为1。
2、 直方图均衡化
直方图是多种空间域处理技术的基础,能有效用于图像增强,是实时图像处理的流行工具,直方图均衡化的目的是使图像在整个灰度值动态变化范围内分布均匀化,改善图像的亮度分布状态,增强视觉效果。直方图均衡化是通过灰度变换将一幅图像转换程另一幅具有均衡性的直方图。即在每个灰度级上都具有相同的像素点数的过程。
3、 空域滤波
手工滤波与函数提供滤波器的比较
六、实验步骤:
1、 将待处理图片拷到matlab软件’work’文件夹
2、 实行亮度变换
3、 对图像进行直方图均衡处理
4、 空域滤波
5、 记录实验结果并分析
七、实验程序及结果记录:
1、 亮度变换
I=imread(‘E:\fig308.tif’);
Imshow(I);
Figure,imhist(I);
J=imadjust(I,[0.5 0.9],[0,1]);
Figure;
imshow(J);
Figure;
imhist(J)
2、直方图均衡化
I=imread(‘E:\fig308.tig’);
J=histea(I);
Imshow(I);
Title(‘原图像’);
Figure;
Imshow(J);
Title(‘直方图均衡化后的图像’);
Figure;
Subplot(1,2,1);
Imhist(I,64);
Title(‘原图像直方图’); Subplot(1,2,2);
Imhist(J,64);
Title(‘均衡变换后的直方图’); Subplot(1,2,2);
Imhist(J,64);
Title(‘均衡变换后的直方图’);
2、 空域滤波
F=imread(‘E:\fig3016.tif’); W4=fspecial(‘laplacian’,1); W8=[1 1 1;-8 1;1 1 1];
F=im2double(f);
G4=f_imfilter(f,w4,’replicate’); G8=f_imfilter(f,w8,’replicate’); Figure;
Subplot(1,3,1);
Imshow(f);
Title(‘原图’);
Subplot(1,3,2);
Imshow(g4);
Title(‘中心为-4拉普拉斯的效果’); Subplot(1,3,3);
Imshow(g8);
Title(‘中心为-8拉普拉斯的效果’);
八、实验结果分析:
亮度变换直方图均衡化可以对图像进行处理,进行空域图像增强。
第二篇:5.图像增强—空域滤波 - 数字图像处理实验报告
计算机与信息工程学院验证性实验报告
专业:通信工程 年级/班级:20##级 20##—20##学年第一学期
一、实验目的
进一步了解MatLab软件/语言,学会使用MatLab对图像作滤波处理,使学生有机会掌握滤波算法,体会滤波效果。
了解几种不同滤波方式的使用和使用的场合,培养处理实际图像的能力,并为课堂教学提供配套的实践机会。
二、实验要求
(1)学生应当完成对于给定图像+噪声,使用平均滤波器、中值滤波器对不同强度的高斯噪声和椒盐噪声,进行滤波处理;能够正确地评价处理的结果;能够从理论上作出合理的解释。
(2)利用MATLAB软件实现空域滤波的程序:
I=imread('electric.tif');
J = imnoise(I,'gauss',0.02); %添加高斯噪声
J = imnoise(I,'salt & pepper',0.02); %添加椒盐噪声
ave1=fspecial('average',3); %产生3×3的均值模版
ave2=fspecial('average',5); %产生5×5的均值模版
K = filter2(ave1,J)/255; %均值滤波3×3
L = filter2(ave2,J)/255; %均值滤波5×5
M = medfilt2(J,[3 3]); %中值滤波3×3模板
N = medfilt2(J,[4 4]); %中值滤波4×4模板
imshow(I);
figure,imshow(J);
figure,imshow(K);
figure,imshow(L);
figure,imshow(M);
figure,imshow(N);
三、实验设备与软件
(1) IBM-PC计算机系统
(2) MatLab软件/语言包括图像处理工具箱(Image Processing Toolbox)
(3) 实验所需要的图片
四、实验内容与步骤
a) 调入并显示原始图像Sample2-1.jpg 。
b) 利用imnoise 命令在图像Sample2-1.jpg 上加入高斯(gaussian) 噪声
c)利用预定义函数fspecial 命令产生平均(average)滤波器
d)分别采用3x3和5x5的模板,分别用平均滤波器以及中值滤波器,对加入噪声的图像进行处理并观察不同噪声水平下,上述滤波器处理的结果;
e)选择不同大小的模板,对加入某一固定噪声水平噪声的图像进行处理,观察上述滤波器处理的结果。
f)利用imnoise 命令在图像Sample2-1.jpg 上加入椒盐噪声(salt & pepper)
g)重复c)~ e)的步骤
h)输出全部结果并进行讨论。
五、思考题/问答题
(1) 简述高斯噪声和椒盐噪声的特点。
(2) 结合实验内容,定性评价平均滤波器/中值滤波器对高斯噪声和椒盐噪声的去噪效果?
(3) 结合实验内容,定性评价滤波窗口对去噪效果的影响?
六、实验报告要求
描述实验的基本步骤,用数据和图片给出各个步骤中取得的实验结果,并进行必要的讨论,必须包括原始图像及其计算/处理后的图像。
七、实验图像
Fig0205(a)(cktboard_200dpi).tif(原始图像)
八、实验程序及结果
clc;clear all
I=imread('Fig0205(a)(cktboard_200dpi).tif');
J1 = imnoise(I,'gauss',0.02); %添加高斯噪声
J2 = imnoise(I,'salt & pepper',0.02); %添加椒盐噪声
ave1=fspecial('average',3); %产生3×3的均值模版
ave2=fspecial('average',5); %产生5×5的均值模版
K1 = filter2(ave1,J1)/255; %均值滤波3×3
L1 = filter2(ave2,J1)/255; %均值滤波5×5
M1 = medfilt2(J1,[3 3]); %中值滤波3×3模板
N1= medfilt2(J1,[5 5]); %中值滤波5×5模板
subplot 121;imshow(I);
title('原始图像')
subplot 122,imshow(J1);
title('加入高斯噪声图像')
figure(2)
subplot 121,imshow(K1);
title('均值滤波3×3后的图像')
subplot 122;,imshow(L1);
title('均值滤波5×5后的图像')
figure
subplot 121;imshow(M1);
title('中值滤波3×3后的图像')
subplot 122;,imshow(N1);
title('中值滤波3×3后的图像')
clc;clear all
I=imread('Fig0205(a)(cktboard_200dpi).tif');
J2 = imnoise(I,'salt & pepper',0.02); %添加椒盐噪声
ave1=fspecial('average',3); %产生3×3的均值模版
ave2=fspecial('average',5); %产生5×5的均值模版
K1 = filter2(ave1,J2)/255; %均值滤波3×3
L1 = filter2(ave2,J2)/255; %均值滤波5×5
M1 = medfilt2(J2,[3 3]); %中值滤波3×3模板
N1= medfilt2(J2,[5 5]); %中值滤波5×5模板
subplot 121;imshow(I);
title('原始图像')
subplot 122,imshow(J2);
title('添加椒盐噪声后的图像')
figure(2)
subplot 121,imshow(K1);
title('均值滤波3×3后的图像')
subplot 122;,imshow(L1);
title('均值滤波5×5后的图像')
figure
subplot 121;imshow(M1);
title('中值滤波3×3后的图像')
subplot 122;,imshow(N1);
title('中值滤波3×3后的图像')
教师签名:
年 月 日