细胞图像处理大作业.docx

PAGE PAGE # / 6 实验报告 实验课程： 数字图像处理 实验名称： 骨髓细胞计数 实验地点： 计算机中心405 UNIVERSITY OF SHANGHAI FOR SCIENCE AND TECHNOLOGY 学生姓名: 学 号： 指导教师： 实验时间： 年 月 一、实验要求： 将被染成紫色的白细胞从较暗的红细胞和背景中分割出来， 去除细胞 粘连，用数字图像处理技术实现图像中紫色细胞的自动计数 。 二、实验过程 1、对原图像进行预处理 （1 ）由于源图像是彩色图像，为了便于进行仿真实验，在进行处理之前需要将彩色图像转 换为灰度图像。用到的函数为 rgb2gray （2 ）获取图像的全局阈值并将图像二值化，以方便后面实验中对图像进行形态学修饰的进 行 （3）在图像获取过程中，由于热传输等影响，会产生高斯噪声，较适合采用中值滤波的方法进 行去噪处理。因为，中值滤波的方法能够较好的保护细胞的边缘信息，满 足本实验的需求； （4 ）删除小面积图像 （5）对图像进行开闭运算。开和闭这两种运算可以除去比结构元素小的特定图像细节，同 时保证不产生全局几何失真。开运算可以把比结构元素小的突刺滤掉， 切断细长搭接而起到 分离作用；闭运算可以把比结构元素小的缺口或孔填充上，搭接短的间隔而起到连接作用。 （6 ）对图像进行孔洞填充， （7）开闭运算会留下一些偏离的孤立像素， 使用bwareaopen来移除少于特定像素个数的斑 点。BW2= bwareaopen（BW,P）从二值图像中移除所以少于 P像素值的连通块，得到另外的二 值图像BW3 实验图片显示如下： File Edit View Insert Tools Desktop Window Help□ 01wJ堺?二伉化图像II除小面积后的二脩化閨像升闭运算后图像 File Edit View Insert Tools Desktop Window Help □ 01 w J堺? 二伉化图像 II除小面积后的二脩化閨像 升闭运算后图像 埴充礼洞后團像 2、分割计数 本次实验采用分水岭法对进行过预处理的图像进行分割和计数， 最后计数得紫色细胞个数 为35 Rle Edit View lln&e/t Tools Desktop Window Help ■ I □ 3 Q ft | h 汽肝国毀諾*氓□ IL11 ■ 0 f a X a 0 o f s . ■■ □ o y< rj^ ? o 0 i/Vork&pace Name Value BW 432x576 logical BW1 432x576 BW2 432x576 logical BW22 432x576 logical BW3 432x576 logical D 432x576 single fo 432^575 fogic&f foe 432x576 logical linage 432x576 uint8 |L 432x576 uint8 number 35 B rgb 432x576x3 uint8 三、实验心得体会 通过这学期的学习， 我获益匪浅，学到了很多有趣的知识，知道了什么是灰度图， 二值 化，阈值，分水岭等很多有趣的图像处理方法，知道了怎么使用 matlab帮助查找自己不会 的函数，也成功的自己独立完成了细胞的计数。 这门课也激起了我对图像处理的兴趣， 课程 结束后也会继续学习更多的函数知识，不会停止这方面的学习。 四、源代码 Image = imread( 'b1.bmp' ); subplot(321); imshow(Image); title( ' 原图 ' ); Image=rgb2gray(Image); subplot(322); imshow(Image); title( ' 灰度图 ' ); Theshold = graythresh(Image); BW = im2bw(Image,1.3*Theshold); subplot(323),imshow(BW); title( ' 二值化图像 ' ); BW1= medfilt2(BW,[13 13]); BW1=bwareaopen(BW,50,4) subplot(3,2,4); imshow(BW1); title( ' 删除小面积后图像 ' ); se=strel( 'disk' ,5); fo=imopen(BW1,se); foc=imclose(fo,se) subplot(3,2,5); imshow(foc); title（ ' 开闭运算后图像 ' ）; BW2=imfill(foc, 'holes' ); subplot(3,2,6); imshow(BW2); title( ' 孔洞填充后图像 '