matlab 实现Canny算子

以下是MATLAB中实现Canny算子的示例代码：

1. 读取图像并将其转换为灰度图像：

img = imread('lena.png'); % 读取图像
gray_img = rgb2gray(img); % 转换为灰度图像

2. 对图像进行高斯滤波：

gaussian_img = imgaussfilt(gray_img, 1); % 高斯滤波

3. 计算图像梯度和方向：

[Gx, Gy] = imgradientxy(gaussian_img); % 计算梯度
[theta, mag] = imgradient(Gx, Gy); % 计算梯度方向和大小

4. 非极大值抑制：

nms = imhmax(mag, 0); % 阈值化
nms = imregionalmax(nms); % 区域最大值
nms = nms & (mag > 0); % 非极大值抑制

5. 双阈值处理：

high_threshold = 0.2 * max(mag(:)); % 高阈值
low_threshold = 0.1 * high_threshold; % 低阈值
strong_edges = mag > high_threshold; % 强边缘
weak_edges = mag > low_threshold & mag <= high_threshold; % 弱边缘

6. 边缘跟踪：

edge_map = false(size(img)); % 边缘图
for i = 1:size(img, 1)
    for j = 1:size(img, 2)
        if strong_edges(i, j)
            edge_map(i, j) = true;
        elseif weak_edges(i, j)
            % 检查八个方向上是否有强边缘
            if any(strong_edges(max(1, i-1):min(size(img, 1), i+1), max(1, j-1):min(size(img, 2), j+1)))
                edge_map(i, j) = true;
            end
        end
    end
end

7. 显示边缘图像：

imshow(edge_map);