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

红细胞计数处理过程使用全局阈值二值化图像，仅使白细胞可以在二值图像中显示；调整二值化图像阈值，使白细胞和红细胞都可在二值图像中显示出来；将两图像进行减操作，从第二幅图像中去除白细胞；对去除白细胞的二值图像进行中值滤波，删除小面积对象，填充空洞等操作；标记最后所得二值图像中的连通区域，并获取区域个数。处理结果最后结果为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);Number1subplot(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);Number2I6=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