moravec算子matlab代码

Moravec算子是一种计算机视觉中的边缘检测算法，用于识别图像中的边缘特征。在MATLAB中，你可以使用内置函数或自己编写代码实现这个算子。以下是一个简单的MATLAB代码示例，用于实现Moravec算子：

function edges = moravecOperator(image)
    % Step 1: 高斯滤波平滑图像
    sigma = sqrt(0.5); % 可调整的标准差
    smoothed_image = imgaussfilt2(image, sigma);

    % Step 2: 计算梯度图像
    gradient_x = imgradient(smoothed_image, 'x');
    gradient_y = imgradient(smoothed_image, 'y');

    % Step 3: 计算Moravec算子矩阵
    M = [0.5 -0.5; -0.5 0.5]; % 算子核心矩阵

    % Step 4: 应用算子
    edge_map = conv2(double(gradient_x), M, 'same') + conv2(double(gradient_y), M', 'same');

    % Step 5: 转换为二值边缘
    threshold = graythresh(edge_map);
    edges = edge_map > threshold;
end

在这个代码中，`imgaussfilt2`用于高斯滤波，`imgradient`计算图像的梯度，`conv2`应用Moravec算子，然后通过阈值处理得到边缘。