I = imread('chrimage.bmp');
figure,imshow(I);
I2 = rgb2gray(I);
s = size(I2);
I4 = 255*ones(s(1), s(2), 'uint8');
I5 = imsubtract(I4,I2);
I3 = medfilt2(I5,[5 5]);
I3 = imadjust(I3);
bw = im2bw(I3, 0.3);
bw = bwareaopen(bw, 10);
figure,imshow(bw);
[labeled,numObjects] = bwlabel(bw,4);
RGB_label=label2rgb(labeled,@spring,'c','shuffle');
figure,imshow(RGB_label);
chrdata = regionprops(labeled,'basic')
allchrs = [chrdata.Area];
num = size(allchrs)
nbins = 20;
figure,hist(allchrs,nbins);
title(num(2))