重写这个matlab程序

以下是不使用 `bwmorph` 函数实现端点检测和骨架化的 MATLAB 程序：

% 读取二值图像
bw = imread('image.jpg');

% 对图像进行腐蚀和膨胀操作，去除噪声
se = strel('disk', 2);
bw = imopen(bw, se);

% 执行端点检测
endpoints = endpoints_detection(bw);

% 判断图形类型
if sum(endpoints(:)) == 0
    disp('该图形是一个湖泊。');
else
    % 执行骨架化操作
    skel = skeletonize(bw);
    
    % 获取端点坐标
    endpoints_coord = get_endpoints_coord(endpoints);
    
    % 计算两个端点之间的距离
    d = pdist(endpoints_coord);
    
    % 判断图形类型
    if d <= 2
        disp('该图形是一个海湾。');
    else
        disp('该图形是一条线。');
    end
end

% 端点检测函数
function endpoints = endpoints_detection(bw)
    % 初始化端点矩阵
    endpoints = false(size(bw));
    
    % 检测端点
    for i = 2:size(bw, 1) - 1
        for j = 2:size(bw, 2) - 1
            if bw(i, j) == 1
                % 判断当前像素是否为端点
                neighbors = bw(i-1:i+1, j-1:j+1);
                if sum(neighbors(:)) == 2
                    if neighbors(1, 2) + neighbors(2, 1) + neighbors(2, 3) + neighbors(3, 2) == 1
                        endpoints(i, j) = 1;
                    end
                end
            end
        end
    end
end

% 骨架化函数
function skel = skeletonize(bw)
    % 初始化骨架化结果
    skel = false(size(bw));
    
    % 迭代骨架化过程，直到无法再进行骨架化
    last = zeros(size(bw));
    while any(bw(:))
        eroded = imerode(bw, strel('disk', 1));
        temp = imdilate(eroded, strel('disk', 1));
        subtracted = bw - temp;
        skel = skel | subtracted;
        bw = eroded;
        if isequal(last, bw)
            break;
        end
        last = bw;
    end
end

% 获取端点坐标函数
function endpoints_coord = get_endpoints_coord(endpoints)
    [y, x] = find(endpoints);
    endpoints_coord = [x, y];
end

这个程序实现了与之前相同的功能，但是没有使用 `bwmorph` 函数。其中，`endpoints_detection` 函数实现了端点检测功能，`skeletonize` 函数实现了骨架化功能，`get_endpoints_coord` 函数实现了获取端点坐标的功能。