编写一段matlab代码实现利用拓展高斯差分提取图像边缘

以下是用 MATLAB 实现利用拓展高斯差分（EGD）提取图像边缘的代码示例：

function [edgeImg] = egdEdgeDetect(img, k, sigma, t)

% img: 输入的灰度图像
% k: EGD 算法中的参数 k
% sigma: EGD 算法中的参数 sigma
% t: 阈值

% 计算 EGD 模板
N = ceil(3 * sigma);
x = -N:N;
h = exp(-x.^2 / (2 * sigma^2)) ./ (sigma * sqrt(2*pi));
d1h = -x .* h / (sigma^2);
d2h = (x.^2 - sigma^2) .* h / (sigma^4);

% 计算拓展高斯差分
G1 = conv2(h, h', img, 'same');
G2 = conv2(d1h, d1h', img, 'same');
G3 = conv2(d2h, d2h', img, 'same');
EGD = sqrt(G1.^2 + k*(G2.^2 + G3.^2));

% 边缘检测
edgeImg = EGD > t;

end

这个函数接受一个灰度图像作为输入，以及 EGD 算法中的三个参数 k、sigma 和阈值 t。它首先计算 EGD 模板，然后使用 `conv2` 函数计算拓展高斯差分。最后，它将 EGD 图像与阈值 t 进行比较，生成二值图像作为输出。