matlabcanny算法中非极大值抑制代码

在Matlab中，Canny边缘检测通常包括两个主要阶段：高斯滤波（平滑）和非极大值抑制（NMS）。非极大值抑制是为了消除噪声并保留边缘强度最大的像素点。

% 加载Canny函数库
load('visionlib')

% 高斯滤波
img = imread('image.jpg'); % 替换为你的图像文件名
blurred_img = imgaussfilt(img, 2); % 2是标准差，可以根据需求调整

% Canny边缘检测
edges = edge(blurred_img, 'canny');

% 非极大值抑制
[nRows, nCols] = size(edges);
sobelX = [-1 0 1; -2 0 2; -1 0 1]; % Sobel算子
sobelY = [-1 -2 -1; 0 0 0; 1 2 1]; % 同理
response = sqrt(double(X.*X + Y.*Y)); % X和Y是Sobel滤波后的结果
maxVal = max(response(:));
[~, locs] = findpeaks(response, 'MinPeakHeight', maxVal*threshold); % threshold是预设的抑制阈值
suppressedEdges = zeros(size(edges));
for i = 1:numel(locs)
    row = locs(i, 1);
    col = locs(i, 2);
    if response(row, col) == maxVal
        suppressedEdges(row, col) = edges(row, col);
    end
end

这里假设`edge`函数用于Canny边缘检测，`findpeaks`函数则用来找到局部峰值。需要注意的是，`threshold`值的选择对结果影响较大，通常需要通过实验来优化。