在Octave中绘制边缘上的梯度矢量场

pkg load image

% 读取图像
img = imread('your_image.jpg');
img = rgb2gray(img); % 转换为灰度图像

% 使用Sobel算子计算梯度
[Gx, Gy] = gradient(img);

% 绘制梯度矢量场
quiver(Gx, Gy);
title('Gradient Vector Field');
xlabel('X-axis');
ylabel('Y-axis');

pkg load image

% 读取图像
img = imread('your_image.jpg');
if size(img, 3) == 3
    img = rgb2gray(img); % 转换为灰度图像
end

% 使用Sobel算子计算梯度
[Gx, Gy] = gradient(img);

% 平滑处理以减少噪声影响
img_smooth = imgaussfilt(img, 1);
[Gx_smooth, Gy_smooth] = gradient(img_smooth);

% 绘制梯度矢量场
figure;
quiver(Gx_smooth, Gy_smooth, 'AutoScale', 'off', 'LineWidth', 0.5);
title('Gradient Vector Field');
xlabel('X-axis');
ylabel('Y-axis');
axis equal; % 确保比例一致