文档介绍：MATLAB程序设计
实验报告

院(系)名称宇航学院
专业名称图像处理与模式识别
学生学号 13151023
学生姓名荣楚君

2016 年 5月
实验基于Matlab水域分割的图像分割实验
一、实验目的
掌握Matlab图像处理工具箱的使用。
学****使用联机帮助、Demo与查找更多的Matlab功能。
熟悉使用Matlab进行交互操作的方法。
学****图像分割的基本知识。
二、实验要求
(一)、对实验过程的要求
1、编写M-文件,实现给定图像的分割。
2、学****Matlab中Demo的相关帮助。
3、实验中生成的中间结果和数据尽量保存。
(二)、对实验报告的撰写要求
1、参考附录的实验报告格式样本撰写。
2、“实验结果”主要包括:
(1)基于标记的和非标记的两种水域分割结果图。
(2)用于水域分割的Marker图。
(3)将水域分割得到结果的轮廓在原图上显示出来。
3、实验报告附录中包含M-文件的源代码,要注意代码的格式。
三、实验环境和设备
1、PC机,Windows XP / 7操作系统。
2、Matlab R2010及其以上版本。

四、实验原理
1、图像分割
图像分割就是把图像分成若干个特定的、具有独特性质的区域并提出感兴趣目标的技术和过程,它是由图像处理到图像分析的关键步骤。图像分割的关键是将前景与背景区分开。
分水岭算法又称为水域分割,是一种基于拓扑理论的数学形态学的分割方法,其基本思想是把图像看作是测地学上的拓扑地貌,图像中每一点像素的灰度值表示该点的海拔高度,每一个局部极小值及其影响区域称为集水盆,而集水盆的边界则形成分水岭。分水岭的概念和形成可以通过模拟浸入过程来说明。在每一个局部极小值表面,刺穿一个小孔,然后把整个模型慢慢浸入水中,随着浸入的加深,每一个局部极小值的影响域慢慢向外扩展,在两个集水盆汇合处构筑大坝,即形成分水岭。水域分割一般是在梯度图像上实现,梯度图上的局部极小值一般代表较平滑的区域,即为目标中心区域或者变化不大的背景区域,集水形成的“分水岭”则位于梯度较大的边缘区域。
2、实现方式
首先要对输入原始图像求取梯度图像,对梯度图像直接做水域变换(watershed)即可得到无标记的水域分割结果。
对梯度图像做前景标记,即在输入图像的每一个目标中心做一个Marker(可手动选取,也可用阈值法得到),在梯度图像上用该Marker作局部极小值的标记,再进行水域变换,即可得到基于标记点的水域分割结果。
3、分割结果展示,在原图像上添加水域分割得到的轮廓如图所示。

(a) (b)
图1 图像配准的例子。
五、实验内容及步骤
(一)学****Matlab帮助的Demo代码
从Matlab帮助界面中打开Demos或Examples,搜索watershed在结果中找到Marker-Controlled Watershed Segmentation,配合中间结果研究该demo如何实现利用水域分割得到图像中的梨。
(二)输入测试图像生成待分割的灰度图像
利用如下代码输入图像:
I=imread('');
(三)计算图像梯度与水域分割
使用sobel算子计算水平梯度与垂直梯度再综合,得到输入图像的梯度图像。在梯度
图像上直接做水域分割,得到非标记的水域分割结果。
(四)计算标记点
在原始图像上通过阈值得到前景的大致分割结果并用形态学腐蚀后作为前景(细胞)标记,也可以通过函数[x y]=ginput()手动标记各细胞的位置,甚至用PS或windows画图画出细胞的标记。标记点可以是真正的点,也可以是有一定面积的区域,但是标记区域一定要位于真实的分割细胞内部。标记图像是与输入图像同大小的,在细胞位置为255,其他位置为0的图像,如下图所示。

(a)输入图像(b)标记图像
(五)由标记图像标记梯度图像中的前景与背景极值
第四步中得到的是前景(细胞)的标记,现在要通过前景标记得到背景标记。函数bwdist可以计算二值图像的距离变换,对距离变换做watershed就可以得到前景与前景之间的背景标记(线型标记)。在梯度图上利用imimposemin(gradmap, backgroundmarker | foregroundmarker)即可得到标记了极值的梯度图像;

(a)由距离变换得到的背景标记(b)标记前景与背景的梯度图像
(六)基于标记的水域分割
对标记了极小值的梯度图像做watershed,得到水域分割结果
(七)结果分析与总结
(1)对比无标记的水域分割与有标记的水域分割结果;
(2)总结水域分割算法整体流程,尤其是标记点的运用;