You can count the number of closed regions using bwlabeln.
I = imread('test2.png') ;
I = rgb2gray(I) ;
[L,num] = bwlabeln(I,6) ; % num gives you closed region
m = regionprops(L) ; % measurements of regions, gives centers of the regions
imshow(I)
hold on
for i = 1:num
c = m(i).Centroid ;
plot(c(1),c(2),'*r') ;
end
num gives you the number of closed regions. You can view them using vislabels(L). Download the function vislabels from here .
