canny算法matlab代码

以下是Canny算法的MATLAB实现代码：

function edges = canny(image, low_threshold, high_threshold)
    gray = rgb2gray(image);  % 转换为灰度图像
    kernel_size = 5;  % 定义Sobel算子的大小
    sigma = 1.4;  % 高斯滤波的标准差
    gaussian_filter = fspecial('gaussian', [kernel_size, kernel_size], sigma);  % 生成高斯滤波器
    blur_gray = imfilter(gray, gaussian_filter, 'replicate');  % 高斯滤波
    [Gx, Gy] = imgradientxy(blur_gray);  % 计算图像的梯度
    [magnitude, ~] = imgradient(Gx, Gy);  % 计算图像的梯度幅值
    edges = edge(magnitude, 'canny', [low_threshold, high_threshold]);  % Canny边缘检测
end

其中，`image`为输入的图像，`low_threshold`和`high_threshold`为Canny算法中的两个阈值，表示较弱和较强的边缘信号。函数返回检测到的边缘图像。此外，代码中还使用了`fspecial`函数生成高斯滤波器，`imgradientxy`计算图像的梯度，`imgradient`计算图像梯度幅值，`edge`函数进行Canny边缘检测。