skeletonImage =imread('Z:\裂缝数据集\Concrete Crack Images for Classification\Positive\00001.jpg'); skeletonImage=rgb2gray(skeletonImage); skeletonImage = imbinarize(skeletonImage);%二值化时可以用阈值分割 neighborImage = bwskel(skeletonImage); [row, col] = find(skeletonImage); % 找到骨架线上的像素坐标 verticalLines = []; for i = 1:length(row) r = row(i); c = col(i); if neighborImage(r-1, c) == 1 && neighborImage(r+1, c) == 1 % 上下方向都有邻居像素，则认为是垂直方向 verticalLines = [verticalLines; [c, r]]; % 添加垂直线的坐标 end end figure; imshow(binaryImage); % 原始二值图像 hold on; for i = 1:size(verticalLines, 1) x = verticalLines(i, 1); y = verticalLines(i, 2); plot([x-1, x+1], [y, y], 'r', 'LineWidth', 2); end hold off;

根据你提供的代码，你正在尝试生成垂直于裂缝骨架线的短线。然而，你在代码中使用了一个未定义的变量`binaryImage`，这可能导致错误。此外，你应该在绘制短线前先显示骨架图像`skeletonImage`，而不是二值化图像。

以下是修正后的代码：

skeletonImage = imread('Z:\裂缝数据集\Concrete Crack Images for Classification\Positive\00001.jpg');
skeletonImage = rgb2gray(skeletonImage);
skeletonImage = imbinarize(skeletonImage); % 二值化时可以用阈值分割
neighborImage = bwskel(skeletonImage);
[row, col] = find(skeletonImage); % 找到骨架线上的像素坐标
verticalLines = [];
for i = 1:length(row)
    r = row(i);
    c = col(i);
    if neighborImage(r-1, c) == 1 && neighborImage(r+1, c) == 1
        % 上下方向都有邻居像素，则认为是垂直方向
        verticalLines = [verticalLines; [c, r]]; % 添加垂直线的坐标
    end
end
figure;
imshow(skeletonImage); % 显示骨架图像
hold on;

for i = 1:size(verticalLines, 1)
    x = verticalLines(i, 1);
    y = verticalLines(i, 2);
    
    plot([x-1, x+1], [y, y], 'r', 'LineWidth', 2);
end
hold off;

请注意，这里假设你的裂缝图像路径和文件名正确。如果路径或文件名不正确，你需要相应地更改`imread`函数中的路径。另外，你可能需要调整二值化图像的阈值，以获得更好的结果。