最近在刷b站的matlab教程 贴一下图像处理2的小作业。
简单但比较inspiring。

---

数米粒原图。

---

pic=imread('rice.png'); bg=imopen(pic,strel('disk',15)); temp=imsubtract(pic,bg); level=graythresh(temp); bw=im2bw(temp,level); subplot(1,3,1);imshow(bw);title('b-w image');hold on; %??? [label, num]=bwlabel(bw,8); %?? disp(['the number of grains is',num2str(num)]) %practice:size of grains for i=1:99 count(i)=0; for j=1:size(label,1) for k=1:size(label,2) if label(j,k)==i count(i)=count(i)+1; end end end end subplot(1,3,2); histogram(count); title('the size of grains'); hold on; %practise:painting green=zeros(256,256); blue=zeros(256,256); red=label; for j1=1:size(red,1) for k=1:size(red,2) if red(j,k)==0 red(j,k)=0; else red(j,k)=255; end end end test=cat(3,red,green,blue); subplot(1,3,3);imshow(test);title('painting');

---

output：
num=99

---

%%bwlabel误差体现在两粒米重叠处 可根据直方图修正。