function[bw_hols] = fnd_hols(bw_ref)

cc_hols=bwconncomp(bw_ref,4);

bw_hols = false(size(bw_ref));

hols_area= struct2cell(regionprops(bw_ref,'Area'));

avr_area = 0;

for i=2:cc_hols.NumObjects

avr_area = hols_area{i} avr_area;

end

avr_area = avr_area/cc_hols.NumObjects - 1 ;

for i=2:cc_hols.NumObjects

if(hols_area{i} < 1.5 * avr_area)

bw_hols(cc_hols.PixelIdxList{i}) = true;

end

end