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

实验报告

课程名称数字图像处理导论专业班级 ____________________ 姓名 _______________________ 学号 _______________________

电气与信息学院

和谐勤奋求是创新

实验题目图像分割实验

实验室DSP室&信号室实验时间

实验类别设计同组人数 2 成绩指导教师签字：

.实验目的

1. 理解图像分割的基本概念；

2. 理解图像边缘提取的基本概念；

3. 掌握进行边缘提取的基本方法；

4. 掌握用阈值法进行图像分割的基本方法。

.实验内容

1. 分别用RobertS,Sobel和拉普拉斯高斯算子对图像进行边缘检测。比较三种算子处理的不同之处；

2. 设计一个检测图1中边缘的程序，要求结果类似图2 ,并附原理说明。

3. 任选一种阈值法进行图像分割.

三.实验具体实现

1. 分别用RObertS,Sobel

处；

和拉普拉斯高斯算子对图像进行边缘检测。比较三种算子处理的不同之l=imread('mri.tif);

imshow(l)

BW1=edge(I,'roberts');

figure ,imshow(BW1),title('

BW2=edge(l,'sobel');

figure,imshow(BW2),title('

BW3=edge(l,'log');

figure,imshow(BW3),title('

用RObertS算子')

用Sobel算子')

用拉普拉斯高斯算子')

图1

比较提取边缘的效果可以看出， sober 算子是一种微分算子，对边缘的定位较精确，但是会漏

去一些边缘细节。而 LaPIaCia n-GaUSSia n 算子是一种二阶边缘检测方法，它通过寻找图象灰度值中

二阶过零点来检测边缘并将边缘提取出来，边缘的细节比较丰富。通过比较可以看出 LaPIaCian-Gaussian 算子比sober 算子边缘更完整，效果更好。 2.

设计一个检测图1中边缘的程序，要求结果类似图 2 ,

并附原理说明。

i=imread('m83.tif); SUbPlot(1,2,1); imhist(i); title('

原始图像直方图');

thread=130∕255; subplot(1,2,2);

i3=im2bw(i,thread);

imshow(i3);

title(' 分割结果');

3. 任选一种阈值法进行图像分割

i=imread('trees.tif);

SUbPlot(1,2,1);

imhist(i);

title(' 原始图像直方图');

thread=100∕255;

subplot(1,2,2); i3=im2bw(i,thread);

imshow(i3);

title(' 分割结果')

附录：可能用到的函数和参考结果************** 报告里不能用参考结果中的图像

1.分别用RobertS,Sobel和拉普拉斯高斯算子对图像进行边缘检测。比较三种算子处理的不同之处；

输入如下代码：

l=imread('tire.tif);

imshow(l)

BW1=edge(I,'roberts');

BW3=edge(I,'log'); figure,imshow(BW3),title(' 2.设计一个检测图3-2中边缘的程序，要求结果类似图 3-3 ,并附原理说明。

利用双峰法

i=imread('example.tif); subplot(1,2,1); imhist(i); title('

原始图像直方图');

thread=130∕255; subplot(1,2,2); i3=im2bw(i,thread); imshow(i3);

figure ,imshow(BW1),title(' BW2=edge(l,'sobel'); figure,imshow(BW2),title('

用RObertS 算子') 用Sobel 算子') 用拉普拉斯高斯算子')

得到:

比较提取边缘的效果可以看出，sober 算子是一种微分算子，对边缘的定位较精确, 但

是会漏去一些边缘细节。而 LaPIaCia n-Gaussia n 算子是一种二阶边缘检测方法，它

通过寻找图象灰度值中二阶过零点来检测边缘并将边缘提取出来，边缘的细节比较丰 富。通过

比较可以看出

LaPIaCian-Gaussian 算子比sober 算子边缘更完整，效果更好。

title(' 分割结果');

根据原图像的直方图，发现背景和目标的分割值大约在130左右，并将灰度图像转为二值图像，分割效果比较理想。

3 •任选一种阈值法进行图像分割选用双峰法对图像进行分割：

i=imread('eight.tif);

SUbPlot(1,2,1);

imhist(i);

title(' 原始图像直方图');

thread=100∕255;

subplot(1,2,2);

i3=im2bw(i,thread);

imshow(i3);

title(' 分割结果');

结果：