数字图像处理课程方案(实验报告).docVIP

上 海 理 工 大 学 计 算 机 工 程 学 院 实 验 报 告 实验名称 红细胞数目统计 课程名称 数字图像处理 姓名 王磊 学号 0916020226 日期 2012-11-27 地点 图文信息中心 成绩 教师 韩彦芳 设计内容： 主题：《红细胞数目检测》 详细说明：读入红细胞图片，通过中值滤波，开运算，闭运算，以及贴标签等方法获得细胞个数中值滤波 开运算 闭运算 贴标签 第一步：读取原图,并显示 A = imread(E:\红细胞3.png ); Image=rgb2gray(A); figure,imshow(Image); title( 【原图); 第二步：进行二值化 Theshold = graythresh(Image);%取得图象的全局域值 Image_BW = im2bw(Image,Theshold);%二值化图象 figure,imshow(Image_BW); title(初次二值化图像); 第三步：进行二值化图像 Image_BW_medfilt= medfilt2(Image_BW,[13 13]); figure,imshow(Image_BW_medfilt); title( 【中值滤波后的二值化图像); 第四步：通过“初次二值化图像”与“中值滤波后的二值化图像”进行“或”运算优化图像效果 Optimized_Image_BW = Image_BW_medfilt|Image_BW; figure,imshow(Optimized_Image_BW); title( 【进行“或”运算优化图像效果); 第五步：优化后二值化图象取反,保证:1为白色,0为黑色 Reverse_Image_BW = ~Optimized_Image_BW; figure,imshow(Reverse_Image_BW); title( 【优化后二值化图象取反); 第六步：填充二进制图像的背景色，去掉细胞内的黑色空隙 Filled_Image_BW = bwfill(Reverse_Image_BW,holes); figure, imshow(Filled_Image_BW); title( 【已填充背景色的二进制图像); 第七步：对图像进行开运算，去掉细胞与细胞之间相粘连的部分 SE = strel(disk,4); Open_Image_BW = imopen(Filled_Image_BW,SE); figure, imshow(Open_Image_BW); title( 【开运算后的图像); 第八步：开始计算细胞数 [Label Number]=bwlabel(Open_Image_BW,8)%初步取得细胞个数 Array = bwlabel(Open_Image_BW,8);%取得贴标签处理后的图像 Sum = []; %依次统计贴标签后数组 for i=1:Number [r,c] = find(Array==i);%获取相同标签号的位置，将位置信息存入[r,c] rc = [r c]; Num = length(rc);%取得vc数组的元素的个数 Sum([i])=Num;%将元素个数存入Sum数组 end Sum N = 0; -----假如Sum数组中的元素大于了1500，表示有两个细胞相连，像素点较多，即分为两个细胞数--- for i=1:length(Sum) if(Sum([i])) 1500 N = N+1; end end Number = Number+N %----统计最终细胞数 第九步：最终检测结果： 四、算法分析 （1）中值滤波 利用中值滤波可以对图像进行平滑处理。算法简单，时间复杂度低，但其对点、线和尖顶多的图像不宜采用中值滤波。很容易自适应化 （）开运算　　先腐蚀后膨胀的过程称为开运算。用来消除小物体、在纤细点处分离物体、平滑较大物体的边界的同时并不明显改变其面积。　　 （）闭运算先膨胀后腐蚀的过程称为闭运算。用来填充物体内细小空洞、连接邻近物体、平滑其边界的同时并不明显改变其面积。　　 《数字图像处理》课程设计 - 1 -