细胞分割

利用染色剂将细胞染色，鉴别死活细胞时，可用以下代码，有注释，供大家分享。

代码：

img=imread('cell.jpg');%读入命名为cell的照片

hsv=rgb2hsv(img);%将照片转化为灰度值照片

s=hsv(:,:,2);%图片S的RGB

bw=im2bw(s,graythresh(s));%将灰度图像转化为二值图像

se=strel('disk',5);%利用strel函数创建一个半径为5mm的圆

bw2=imclose(bw,se);%闭操作imclose通常消弥狭窄的间断和长细的鸿沟，消除小的空洞，并填补轮廓线中的断裂。

bw3=bwareaopen(bw2,200);%bwareaopen函数删除二值图像bw2中面积小于200mm2的对象，默认情况下使用8邻域。

imshow(bw3);%显示图片bw3

perim=bwperim(bw3);%从输入图像bw3中返回只包括对象边缘像素点的图像

r=img(:,:,1);%红色分量图像

g=img(:,:,2);%绿色分量图像

b=img(:,:,3);%蓝色分量图像

r(perim)=255;%边缘像素红色分量为255，即边缘颜色是绿色

g(perim)=0;%边缘像素绿色分量为0

b(perim)=0;%边缘像素蓝色分量为0

img(:,:,1)=r;%原图与效果图对比红色分量

img(:,:,2)=g;%原图与效果图对比绿色分量

img(:,:,3)=b;%原图与效果图对比蓝色分量

imshow(img)%显示照片