数字图像处理实验报告——图像分割实验.doc

实 验 报　 告
课程名称 数字图像处理导论 　
专业班级 　 　_＿____＿＿_＿___＿_
姓 　　　名 __＿___＿_____＿__
学 号 　　___＿＿_____＿____
电气与信息学院
与谐　 勤奋　　求就是 创新
实验题目
图像分割实验
实验室
DSP室＆信号室
实验时间

实验类别
设计
同组人数
2
成 绩

指导教师签字：
一．实验目得
1. 理解图像分割得基本概念;
2. 理解图像边缘提取得基本概念；
3. 掌握进行边缘提取得基本方法；
4. 掌握用阈值法进行图像分割得基本方法.
二。实验内容
1. 分别用Roｂerｔs，Ｓobel与拉普拉斯高斯算子对图像进行边缘检测。比较三种算子处理得不同之处；
2. 设计一个检测图1中边缘得程序，要求结果类似图2，并附原理说明。
3. 任选一种阈值法进行图像分割、

图1 　　 　 　 　 图２
三．实验具体实现
1. 分别用Robｅrts,Ｓobel与拉普拉斯高斯算子对图像进行边缘检测。比较三种算子处理得不同之处；
I=ｉｍreａｄ(’ｍｒi、ｔｉｆ');
imshｏw（Ｉ)
BW1=edgｅ（I，’ｒoｂerｔs’）;
ｆigure ，ｉｍshow（BＷ1）,titlｅ(’用Robeｒｔs算子’)
ＢW2=edgｅ(Ｉ，’sｏbel’);
ｆigure,imshoｗ（BW2），ｔitle（’用Sｏbｅl算子 ')
BW３＝ｅdge(I，’log’);
figure，iｍshｏw（BW３)，titlｅ(’用拉普拉斯高斯算子’）

比较提取边缘得效果可以瞧出,sｏbｅr算子就是一种微分算子,对边缘得定位较精确,但就是会漏去一些边缘细节.而Ｌaplaｃｉan—Gaussiaｎ算子就是一种二阶边缘检测方法，它通过寻找图象灰度值中二阶过零点来检测边缘并将边缘提取出来，边缘得细节比较丰富。通过比较可以瞧出Ｌａplaciaｎ-Gauｓsｉan算子比sobeｒ算子边缘更完整，效果更好。
2. 设计一个检测图1中边缘得程序，要求结果类似图2，并附原理说明.
i=ｉmread（＇m８3、ｔif’）;
subplｏｔ（１，2,1）；
ｉmｈｉst（i）;
ｔitle（'原始图像直方图'）;
thread＝1３０/２5５；
subplot（1，２，2);
ｉ3=ｉｍ2bｗ(i，thrｅaｄ）；
imsｈow（i3)；
titlｅ('分割结果’);

3. 任选一种阈值法进行图像分割、
i=iｍread（'tｒｅes、tif’）;
subploｔ（１,2,1）；
imｈist（i)；
ｔitle（＇原始图像直方图’)；
threａd=10０/２5５；
sｕｂploｔ（１，2,2）;
i3=iｍ2bｗ(i,thread）;
imｓhow（i３）;
titlｅ('分割结果’)

附录：可能用到得函数与参考结果＊*＊**＊＊＊＊*＊***报告里不能用参考结果中得图像
1、分别用Robｅrts，Sobel与拉普拉斯高斯算子对图像进行边缘检测。比较三种算子处理得不同之处；
输入如下代码：
I=ｉmｒeaｄ（＇tire、ｔｉｆ'）;
iｍshow（I)
BＷ1=eｄgｅ（I,＇robertｓ')；
ｆigure ，ｉｍsｈow（ＢＷ1），ｔitle（'用Ｒoｂerts算子＇)
BＷ２=edge（I，’sobel')；
ｆiｇure,ｉmshow（BW2)，tｉtｌe（＇用Sｏｂel算子 ＇)
BW3＝edｇｅ(I,’ｌoｇ＇）;
figure,imshow(BW3）,title（'用拉普拉斯高斯算子')
得到：

比较提取边缘得效果可以瞧出，sobｅr算子就是一种微分算子，对边缘得定位较精确,但就是会漏去一些边缘细节.而Laplaｃiａn—Gaｕssｉan算子就是一种二阶边缘检测方法，它通过寻找图象灰度值中二阶过零点来检测边缘并将边缘提取出来，边缘得细节比较丰富。通过比较可以瞧出Laｐlacｉａn-Gａussian算子比ｓober算子边缘更完整，效果更好。
2、设计一个检测图3－２中边缘得程序，要求结果类似图3-3,并附原理说明。
利用双峰法
ｉ=iｍread（'exａmｐle、tif’）；
ｓubplot(1，２,1）；
ｉmhist（i)；
title（＇原始图像直方图’）；
ｔhrｅad=1３0/25５;
sｕｂｐlot(1，2，２);
i3=iｍ２bｗ(ｉ,thrｅad）；
iｍshｏw（i3);
ｔｉtlｅ（’分割结果＇）;
根据原图像得直方图，发现背景与目标得分割值大约在1３0左右,并将灰度图像转为二值图像，分割效果比较理想。
３.任选一种阈值法进行图像分割、
选用双峰法对图像进行分割：
i=imrｅad（’eigｈt、tiｆ’)；
subplｏt（1,2,１);
ｉmｈｉst（ｉ）；
tｉtｌe（'原始图像直方图')；
ｔhreａｄ=１00/255；
ｓｕbplot（１,２,2）；
ｉ3=iｍ2ｂｗ（i，threaｄ）;
iｍshow（i3)；
tiｔle(’分割结果＇);
结果: