细胞图像处理大作业

僦上渔理工人厚

UNIVERSITYOFSHANGHAIFORSCIENCEANDTECHNOLOGY

实验报告

实验课程：_____________数字图像处理_______________

实验名称：____________骨髓细胞计数________________

实验地点：____________计算机中心405____________

学生姓名：_______________________________________

学号：_______________________

指导教师：________________________________________

实验时间：年月

一、实验要求：

将被染成紫色的白细胞从较暗的红细胞和背景中分割出来,去除细胞

粘连，用数字图像处理技术实现图像中紫色细胞的自动计数。

二、实验过程

1、对原图像进行预处理

（1）由于源图像是彩色图像，为了便于进行仿真实验，在进行处理之前需要将彩色图像转

换为灰度图像。用到的函数为rgb2gray

（2）获取图像的全局阈值并将图像二值化，以方便后面实验中对图像进行形态学修饰的进

行

（3）在图像获取过程中，由于热传输等影响，会产生高斯噪声，较适合采用中值滤波的方法进

行去噪处理。因为，中值滤波的方法能够较好的保护细胞的边缘信息，满足本实验的需求：

（4）删除小面积图像

（5）对图像进行开闭运算。开和闭这两种运算可以除去比结构元素小的特定图像细节，同

时保证不产生全局几何失真。开运算可以把比结构元素个的突刺滤掉，切断细长搭接而起到

分离作用；闭运算可以把比结构元素小的缺口或孔填充上，搭接短的间隔而起到连接作用。

（6）对图像进行孔洞填充,

（7）开闭运算会留下一些偏离的孤立像素，使用bwareaopen来移除少于特定像素个数的斑

点。BW2=bwareaopen（BMP）从二值图像中移除所以少于P像素值的连通块，得到另外的二

值图像RW3«

实验图片显示如下：

FileEditViewInsertToolsDesktopWindowHelp、

00己⑨骐/TO□国■⑪

原图灰度图

’4•哈&•

二值化图像删除小面积后的二值化图像

开闭诺身后图像填充孔洞后图像

2、分割计数

本次实验采用分水岭法对进行过预处理的图像进行分割和计数，最后计数得紫色细胞个数

为35

FileEditViewInsertToolsDesktopWindowHelp、

分水岭算法结果

Workspace

Name▲Value

BW432x576logical

BW1432x576logical

BW2432x576logical

BW22432x576logical

BW3432x576logical

D432x576single

fo432x576logical

foe432x576logical

Image432x576uint8

L432x576uint8

number35

rgb432x576x3uint8

三、实验心得体会

通过这学期的学习，我获益匪浅，学到了很多有趣的知识，知道了什么是灰度图，二值

化，阈值，分水岭等很多有趣的图像处理方法，知道了怎么使用mallab帮助查找自己不会

的函数，也成功的自己独立完成了细胞的计数。这门课也激起了我对图像处理的兴趣，课程

结束后也会继续学习更多的函数知识，不会停止这方面的学习。

四、源代码

Image=imread('bl.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,[1313]);

BWl=bwareaopen(BW,50,4)

subplot(3,2,4);

imshow(BW1);

title('删除小面积后图像，);

se=strel(1disk1,5);

fo=imopen(BW1,se);

foc=imclose(fo,se)

subplot(3,2,5);

imshow(foe);

title（1开闭运算后图像，）;

BW2=imfill(foe,*holes');

subplot(3,2,6);

imshow(BW2);

title「孔洞填充后图像，)；

BW3=bwareaopen(BW2,20);

D=bwdist(~BW3);

D=-D;

D(-BW3)=-Inf;

L=watershed(D);

rgb=label2rgb(L,'jet');

figure;

imshow(rgb,[]);

title('分水岭分割计算后图像，)；

number=max(max(L))-L;

number

方法一源程序：

A=imreadCyl.bmp');

B=rgb2gray(A);

h=[010,1-41,010];

Jl=conv2(B,J'same');

Kl=medfilt2(Jl,[33]);"

BW=edge(KI,'canny');

SE=strel(*lino*,11,90);

BWl=imdilate(BW,SE);

BW2=bwmorph(BW1,*thin：,Inf);

[H,theta,rho]=hough(BW2,1ThetaResolution，,0.2);

subplot(221),imshow(H,[],'XData',theta,,YData',rho,,InitialMagnification,,

'fit')

axison,axisnormal

xlabel(5\thela,),ylabel('\rho'),

peaks=houghpeaks(H,5);

holdon

plot(theta(peaks(:,2)),rho(peaks(:,1)),…

linestyle,none,marker,s,color,w);

lines=hough1ines(BW2,theta,