matlab实现canny边缘检测不调用edge函数的具体代码

以下是用 MATLAB 手动实现 Canny 边缘检测的代码：

% 读取图像
img = imread('lena.bmp');
% 将图像转化为灰度图像
grayImg = rgb2gray(img);
% 对灰度图像进行高斯滤波
gaussImg = imgaussfilt(grayImg, 1.5);
% 计算图像梯度
[Gx,Gy] = imgradientxy(gaussImg);
% 计算梯度的幅值和方向
[Gmag,Gdir] = imgradient(Gx,Gy);
% 非极大值抑制
nmsImg = nonmaxsup(Gmag,Gdir);
% 双阈值检测
thHi = 0.1 * max(nmsImg(:));
thLo = 0.05 * max(nmsImg(:));
edgeMap = hysthresh(nmsImg, thHi, thLo);
% 显示边缘图像
imshow(edgeMap);

以上代码的详细说明：

1. 读取图像：使用 imread 函数读取需要进行边缘检测的图像。

2. 转化为灰度图像：使用 rgb2gray 函数将彩色图像转化为灰度图像（前提是原图为彩色图像）。

3. 高斯滤波：使用 imgaussfilt 函数对灰度图像进行高斯滤波，以平滑图像并降噪。

4. 计算梯度：使用 imgradientxy 函数计算灰度图像的梯度，得到水平方向和垂直方向的梯度值。

5. 计算梯度的幅值和方向：使用 imgradient 函数计算梯度的幅值和方向。

6. 非极大值抑制：使用 nonmaxsup 函数对梯度幅值图像进行非极大值抑制，以保留局部最大值，得到较细的边缘。

7. 双阈值检测：使用 hysthresh 函数对抑制后的梯度幅值图像进行双阈值检测，以确定边缘的强度和弱度，得到二值化的边缘图像。

8. 显示边缘图像：使用 imshow 函数显示生成的边缘图像。