E1=imread('MT-20min47uM-2006-0037.tif');
E11=mat2gray(E1);
E12=imadjust(E11);
bwE=im2bw(E12,0.7);


E13 = imfill(bwE,'holes');
finalE1 = bwlabel(E13);
PhasePropsE1 = regionprops(finalE1,'area');
D21=[PhasePropsE1(:).Area];

E2=imread('MT-20min47uM-2006-0038.tif');
E21=mat2gray(E2);
E22=imadjust(E21);
bwE2=im2bw(E22,0.7);

E23 = imfill(bwE2,'holes');
finalE2 = bwlabel(E23);
PhasePropsE2 = regionprops(finalE2,'area');
D22=[PhasePropsE2(:).Area];


E3=imread('MT-20min47uM-2006-0041.tif');
E31=mat2gray(E3);
E32=imadjust(E31);
bwE3=im2bw(E32,0.7);

E33 = imfill(bwE3,'holes');
finalE3 = bwlabel(E33);
PhasePropsE3 = regionprops(finalE3,'area');
D23=[PhasePropsE3(:).Area];


E4=imread('MT-20min47uM-2006-0042.tif');
E41=mat2gray(E4);
E42=imadjust(E41);
bwE4=im2bw(E42,0.7);

E43 = imfill(bwE4,'holes');
finalE4 = bwlabel(E43);
PhasePropsE4 = regionprops(finalE4,'area');
D24=[PhasePropsE4(:).Area];


E5=imread('MT-20min47uM-2006-0044.tif');
E51=mat2gray(E5);
E52=imadjust(E51);
bwE5=im2bw(E52,0.7);

E53 = imfill(bwE5,'holes');
finalE5 = bwlabel(E53);
PhasePropsE5 = regionprops(finalE5,'area');
D25=[PhasePropsE5(:).Area];



%Group the data together
DDDDD=[D21 D22 D23 D24 D25];

%compute hist numbers
j=1;
for i=1:2782
    if (DDDDD(i)>100&&DDDDD(i)<3000)
       EEEEE(j)=DDDDD(i);
       j=j+1;
    end
end