图像差分法实现怎么实现 matlab代码演示

图像差分法是一种常用的图像边缘检测方法，下面是使用 MATLAB 实现的代码演示：

% 读取图像
I = imread('lena.png');

% 转换为灰度图像
I = rgb2gray(I);

% 定义 Sobel 模板
Sobel_x = [-1 0 1; -2 0 2; -1 0 1];
Sobel_y = [-1 -2 -1; 0 0 0; 1 2 1];

% 计算水平和垂直方向的梯度
Gx = conv2(double(I), Sobel_x, 'same');
Gy = conv2(double(I), Sobel_y, 'same');

% 计算梯度幅值和方向
G = sqrt(Gx.^2 + Gy.^2);
theta = atan2(Gy, Gx);

% 显示梯度幅值
imshow(G, []);

% 显示梯度方向
quiver(Gx, Gy);

在这个代码中，我们首先读取了一张图像，并将其转换为灰度图像。然后，我们定义了 Sobel 模板，用于计算水平和垂直方向的梯度。接着，我们使用 conv2 函数计算梯度，并使用 sqrt 和 atan2 函数计算梯度幅值和方向。最后，我们使用 imshow 和 quiver 函数分别显示梯度幅值和方向。