matlab实现同类物体简单计数

要计数的图片是这样子的：

想练习一下，看到旁边的这个电蚊香片挺合适，随手拿来用一下啦~~~
实现代码如下，注释还没写完，不过挺简单的，过两天把思路和注释补齐~~~

clc;%清除命令窗口中的内容
close all;%关闭所有的figure窗口
clear all;%清楚工作空间所有的变量、函数等
I=imread ('C:\Users\admin\Desktop\蚊香片.jpg');%读取蚊香图片
subplot(331);%设置图片显示样式，三行三列，显示位置为第一张
imshow (I) ;%显示原图
title('原图');%命名为“原图”
I = imnoise(I,'speckle',0.8);%加入斑点噪声，模拟现实噪音干扰情况
subplot(332);%显示位置为第二张
imshow(I);%显示加噪后的图片
title('加噪');%命名为“加噪”
R=I(:,:,1);
G=I(:,:,2);
B=I(:,:,3);%将加噪后的彩色图片的R,G,B三分量数据分别存储到矩阵R,G,B中
R= wiener2(R,[15,15]);%对R分量矩阵进行维纳滤波
G= wiener2(G,[10,10]);%对G分量矩阵进行维纳滤波
B= wiener2(B,[5,5]);%对B分量矩阵进行维纳滤波
I(:,:,1)=R;
I(:,:,2)=G;
I(:,:,3)=B;%将滤波后的R,G,B分量再还原给图片，I为滤波处理后的彩色图片
subplot(333);%显示位置为第三个
imshow(I);%显示滤波后的图片
title('维纳滤波后');%命名为“维纳滤波后”
R=I(:,:,1);
G=I(:,:,2);
B=I(:,:,3);%将图片的R、G、B分量数据分别存储到矩阵R,G,B中
R((R<=60|R>=130)|(G<=120|G>=180))=255;%使用photoshop取色工具分别取蚊香片和背景的的R,G,B分量
G((R<=60|R>=130)|(G<=120|G>=180))=255;
B((R<=60|R>=130)|(G<=120|G>=180))=255;%将R,G,B三个分量矩阵中的数据进行处理，不属于蚊香片颜色区间的置为白色
I(:,:,1)=R;
I(:,:,2)=G;
I(:,:,3)=B;%将处理过的三色数据赋给原图
I1=rgb2gray(I);%将真彩图像转换为灰度图像
subplot(334);%显示位置为第四个
imshow(I1);%显示处理后的灰度图像
title('灰度图像');%命名为“灰度图像”
thresh=graythresh(I1);%使用grathresh函数，找到灰度图像的合适阈值
I2=im2bw(I1,thresh);%使用阈值将灰度图像变成二值图像
subplot(335);%显示位置为第五个
imshow(I2);%显示二值图像I2
title('二值图像');%命名为“二值图像”
SE0=strel('disk',5);%在进行膨胀、腐蚀操作前，需要创建一个掩膜（mask），这里创建一个半径为5的平面圆盘状掩膜，SE0是一个矩阵
SE=strel('disk',12);%在进行膨胀、腐蚀操作前，需要创建一个掩膜（mask），这里创建一个半径为12的平面圆盘状掩膜，SE是一个矩阵
I2=imdilate(I2,SE0);%进行膨胀操作，白色部分扩张，消除维纳滤波没过滤干净的噪声（主要为边界区域）
subplot(336);%显示位置为第6个
imshow(I2);%显示二值图像I2
title('膨胀去噪声（边界）');%命名为“膨胀去噪声（边界）”
I3=imerode(I2,SE);%使用imrode函数进行腐蚀操作，第一个参数是二值图像，第二个参数是掩膜,腐蚀操作是对白色部分，这一步填充了蚊香中的空隙
subplot(337);%显示位置为第7个
imshow(I3);%显示腐蚀后的图像
title('腐蚀填充空隙');
SE1=[1 1 1 1;1 1 1 1;1 1 1 1;1 1 1 1];%定义mask，为4x4矩阵
I4=imdilate(I3,SE1);%此时进行腐蚀操作，是为了消除蚊香片之间的连接部分，使其独立出来
I4=~I4;%取反，为了使用bwlabel函数进行计数
subplot(338);
imshow(I4);
title('取反');
[L,amount]=bwlabel(I4,8);%bwlabel对二值化中的白色区域进行计数
fprintf('该图片共有%d个电蚊香！',amount);