To find the actual boundary, invert, fill, and sum
binaryImage = imfill(~binaryImage, 'holes');
numPixels = sum(binaryImage(:));
If you want a smooth boundary, you need to use activecontour(), like in my attached demo.
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!