canny算子裂缝边缘检测matlab代码

Canny算子是一种经典的图像边缘检测算法，它由John F. Canny于1986年提出，能有效识别出图像中的强边缘并消除噪声。在MATLAB中，你可以使用内置函数`edge`结合`imfilter`来实现Canny边缘检测。下面是一个简单的示例代码：

% 加载图像
img = imread('your_image.jpg'); % 替换为你要处理的实际图片路径

% 高斯滤波减少噪声
gaussian_kernel_size = 5; % 可调整大小
blurred_img = imgaussfilt(img, gaussian_kernel_size);

% 计算梯度幅值和方向
gradient_mag = sqrt(gradient(blurred_img,'dx')^2 + gradient(blurred_img,'dy')^2);
gradient_angle = atan2(gradient(blurred_img,'dy'), gradient(blurred_img,'dx'));

% 使用非极大值抑制 (Non-Maximum Suppression)
grad_direction = imrotate(angle(gradient_angle), -angle(gradient_angle));
sobelx = fspecial('sobel');
sobely = imrotate(sobelx, -90); % 转置sobelx滤波矩阵用于y轴方向
sobel_mag = sqrt(double(sobelx.*sobelx) + double(sobely.*sobely));
[~, edges] = bwmorph(hysteresis(threshold(grad_mag, 'gauss', 2)), 'skel', 4);

% 最后的双阈值处理
low_threshold = 0.1 * max(max(edges));
high_threshold = 0.3 * max(max(edges));
edges = bwareaopen(edges, low_threshold*sum(sum(edges)));
edges = imfill(edges,'holes');

% 显示结果
figure;
subplot(1,2,1), imshow(img), title('Original Image');
subplot(1,2,2), imshow(edges), title('Canny Edge Detection Result');

在这个例子中，你需要替换`'your_image.jpg'`为你实际要处理的图片文件名。