请用matlab语言写出canny边缘检测中非极大值抑制代码

在MATLAB中，Canny边缘检测通常包括非极大值抑制（Non-Maximum Suppression, NMS）来减少噪声响应并保持边缘强度。下面是使用MATLAB实现非极大值抑制的一个简化示例[^1]:

% 假设你已经有了梯度幅值（gradient magnitude）和方向（gradient angle）
[Gx, Gy] = gradient(I); % I是输入图像
theta = atan2(Gy, Gx);

% 选择角范围内的像素索引
inds = find(theta > -pi/2 & theta < pi/2); % 取0到180度之间的像素

% 初始化存储强边缘结果的新图像
edges = zeros(size(I));

% 非极大值抑制
for i = inds
    % 计算当前像素及其邻域内的局部极值
    local_max = false(1, 3);
    if i == 1 || abs(theta(i)-theta(i-1)) > pi/2
        local_max(1) = true; % 左边
    end
    if i == size(I, 2) || abs(theta(i)-theta(i+1)) > pi/2
        local_max(end) = true; % 右边
    else
        if Gy(i) >= Gy(i-1) && Gy(i) >= Gy(i+1)
            local_max(2) = true; % 上方
        elseif Gy(i) <= Gy(i-1) && Gy(i) <= Gy(i+1)
            local_max(3) = true; % 下方
        end
    end
    
    % 如果当前像素是局部最大值，则保留
    if any(local_max)
        edges(i) = Gx(i);
    end
end

% 结果即为边缘图
imshow(edges)

请注意，这个例子假设你已经计算了图像的梯度，`Gx`和`Gy`分别代表水平和垂直方向上的梯度分量。非极大值抑制部分依赖于特定的边界条件来确定哪些像素是局部极大值。