基于弧长法计算图片中裂缝宽度的matlab代码

以下是基于弧长法计算图片中裂缝宽度的MATLAB代码：

% 读取图片
I = imread('image.jpg');

% 转化为灰度图
I_gray = rgb2gray(I);

% 二值化
level = graythresh(I_gray);
I_bw = imbinarize(I_gray, level);

% 边缘检测
BW = edge(I_bw, 'canny');

% 查找连通区域
CC = bwconncomp(BW);

% 计算弧长和宽度
arc_lengths = zeros(1, CC.NumObjects);
widths = zeros(1, CC.NumObjects);
for i = 1:CC.NumObjects
    % 获取连通区域的像素坐标
    [r, c] = ind2sub(CC.ImageSize, CC.PixelIdxList{i});
    % 计算弧长
    arc_lengths(i) = arc_length(r, c);
    % 计算宽度
    widths(i) = width(r, c);
end

% 显示结果
imshow(I);
hold on;
for i = 1:CC.NumObjects
    text(CC.PixelIdxList{i}(1), CC.PixelIdxList{i}(2), num2str(widths(i)), 'Color', 'red');
end

% 弧长函数
function l = arc_length(r, c)
    x = c - mean(c);
    y = r - mean(r);
    theta = atan2(y, x);
    [~, order] = sort(theta);
    x = x(order);
    y = y(order);
    dx = diff(x);
    dy = diff(y);
    l = sum(sqrt(dx.^2 + dy.^2));
end

% 宽度函数
function w = width(r, c)
    x = c - mean(c);
    y = r - mean(r);
    theta = atan2(y, x);
    [~, order] = sort(theta);
    x = x(order);
    y = y(order);
    dx = diff(x);
    dy = diff(y);
    l = cumsum(sqrt(dx.^2 + dy.^2));
    w = max(l) - min(l);
end

其中 `arc_length` 函数计算弧长，`width` 函数计算宽度。首先将图片转化为灰度图，然后二值化，进行边缘检测，查找连通区域。对于每个连通区域，计算它的弧长和宽度，并在图片上显示出来。