MATLAB在细胞图像分割与计数方面的应用..ppt

方法三 大津阈值法 步骤 1.用阈值t 将像素划分为两类C1和C2 步骤 2.计算背景和目标的出现概率 步骤 3.计算A和B两个区域的类间方差 Matlab中提供了计算最大类间方差阈值分割的阈值函数graythresh 语句格式：level=graythresh(I) i=imread(bai.jpg); %读入细胞图片 im2=imadjust(i,[.3 .3 0;.7 .7 1],[]); %图像预处理，亮度调节 im1=rgb2gray(im2); %将图像转换成灰度图 level=graythresh(im1); %大津阈值分割 BW=im2bw(im1,level); %将灰度图像转换为二值图像 总程序如下： se=strel(line,1,90); se1=strel(line,1,0); BW1=imdilate(BW,se); %对图像进行膨胀操作 BW2=imerode(BW1,se1); %再对图像进行腐蚀操作 BW3=bwareaopen(BW2,90); BW3=~BW3; BW4=bwareaopen(BW3,50); BW4=~BW4; %以上部分是对图像进行去噪处理，将面积小的噪声去掉 k=bwlabel(~BW4); %图像边缘提取 a=max(max(k)) %细胞计数 BW5=i; for i=1:3; BW5(:,:,i)=BW5(:,:,i).*uint8(~BW4); end %以上部分是大津阈值分割后图像与原图像彩色相乘，使效果更清晰 imshow(BW5); %显示分割后再填充颜色的最终效果图 原图像 VS imadjust语句调亮后图像 用rgb2gray转为灰度图像 再用histeq增强对比度 大津阈值法： level=graythresh(im)  BW=im2bw(im,level)  得到二值图像 只膨胀操作 只腐蚀操作 先膨胀后腐蚀 bwareaopen 语句删除小面积对象 用此图进行 边缘提取和计数 计数并颜色填充 Ans=243