图像处理—红细胞计数(Matlab).doc

红细胞计数 处理过程 使用全局阈值二值化图像，仅使白细胞可以在二值图像中显示； 调整二值化图像阈值，使白细胞和红细胞都可在二值图像中显示出来； 将两图像进行减操作，从第二幅图像中去除白细胞； 对去除白细胞的二值图像进行中值滤波，删除小面积对象，填充空洞等操作； 标记最后所得二值图像中的连通区域，并获取区域个数。 处理结果 最后结果为101。 处理效果不太理想。 程序清单 clear; close all; Image = imread(123.jpg); subplot(3,2,1); imshow(Image); title(原图); Image=rgb2gray(Image);%取灰度图像 subplot(3,2,2); imshow(Image); title(灰度图); Theshold = graythresh(Image);%取得图象的全局域值 Image_BW = im2bw(Image,1.3*Theshold);%二值化图象 I=im2bw(Image,Theshold); Reverse_Image_BW22=~Image_BW;%反相 subplot(3,2,3); imshow(Image_BW); title(初次二值化图像); %Image_BW_medfilt= medfilt2(Image_BW,[13 13]);%中值滤波 Image_BW_medfilt=bwareaopen(Reverse_Image_BW22,50,4); subplot(3,2,4); imshow(~Image_BW_medfilt); title(删除小面积对象后的二值化图像); Reverse_Image_BW = ~Image_BW_medfilt; I1=~Reverse_Image_BW-~I; I2=imfill(I1,holes); subplot(3,2,5); imshow(I); title(白细胞); I3=imclearborder(I2); [Label, Number1]=bwlabel(I3,8); Number1 subplot(3,2,6);imshow(~I1);title(相减图像); I4=bwareaopen(I1,100,4); I5=imfill(I4,holes); figure; subplot(221);imshow(~I4);title(删除小面积对象); subplot(222);imshow(~I5);title(空洞填充); [Label, Number2]=bwlabel(I5,8); Number2 I6=medfilt2(I5,[10 10]); subplot(223);imshow(~I6);title(中值滤波后图像) I7=bwareaopen(I6,100,4); subplot(224);imshow(~I7);title(第三次删除小面积对象) [Label, Number3]=bwlabel(I7,8); Number3