I = imread('b.bmp');
se = strel('square', 3); % 3x3正方形结构元素
Ie = imerode(I, se); % 腐蚀操作得到内部点
Iout = I - Ie; % 通过减去内部点得到边界
Iout = ~Iout; % 反转得到外轮廓
figure, imshow(Iout); % 显示外轮廓
[L, num] = bwlabel(Iout, 8); % 连通块标记
F = L > 1; % 排除背景区域
BW2 = imfill(F, 'holes'); % 填充区域
figure, imshow(1 - BW2); % 显示填充后的区域