% 读入原始图像 I = imread('D:\课设图片\1.jpg'); % 对图像进行二值化处理 I1 = im2bw(I, 0.6); % 对二值化后的图像进行形态学处理 se = strel('disk', 5); I2 = imclose(I1, se); I3 = imopen(I2, se); I4 = imopen(I3, se); % 对处理后的图像进行连通域分析 L = bwlabel(I4); STATS = regionprops(L, 'all'); % 选取面积在 2000 到 3800 之间的连通域，并在原图像中显示其边缘 S = size(I, 1) * (size(I, 2) * 0.85 - size(I, 2) / 10) * (size(I, 1) * 2 / 3 - size(I, 1) / 4); figure; imshow(I); hold on; for i = 1:length(STATS) area = STATS(i).ConvexArea; if area >= 1500 && area <= 4800 bbox = STATS(i).BoundingBox; rectangle('Position', bbox, 'EdgeColor', 'r'); boundary = bwtraceboundary(I4, [round(bbox(2)), round(bbox(1))], 'N'); plot(boundary(:,2), boundary(:,1), 'g', 'LineWidth', 2); end end title('在原图像中显示选取的连通域边缘');改进代码，直接在原图显示连通域

可以将连通域的边缘直接绘制在原始图像上，无需再新建一个figure。改进后的代码如下：

% 读入原始图像
I = imread('D:\课设图片\1.jpg');

% 对图像进行二值化处理
I1 = im2bw(I, 0.6);

% 对二值化后的图像进行形态学处理
se = strel('disk', 5);
I2 = imclose(I1, se);
I3 = imopen(I2, se);
I4 = imopen(I3, se);

% 对处理后的图像进行连通域分析
L = bwlabel(I4);
STATS = regionprops(L, 'all');

% 选取面积在 2000 到 3800 之间的连通域，并在原图像中显示其边缘
figure; imshow(I); hold on;
for i = 1:length(STATS)
    area = STATS(i).ConvexArea;
    if area >= 1500 && area <= 4800
        bbox = STATS(i).BoundingBox;
        rectangle('Position', bbox, 'EdgeColor', 'r');
        boundary = bwtraceboundary(I4, [round(bbox(2)), round(bbox(1))], 'N');
        plot(boundary(:,2), boundary(:,1), 'g', 'LineWidth', 2);
    end
end
title('在原图像中显示选取的连通域边缘');