matlab canny算子边缘检测函数代码

I = imread('lena.bmp');                 %%如果是其他类型图像,请先转换为灰度图 

%%没有噪声时的检测结果 

BW_sobel = edge(I,'sobel');

BW_prewitt = edge(I,'prewitt'); 

BW_roberts = edge(I,'roberts'); 

BW_laplace = edge(I,'log');

BW_canny = edge(I,'canny'); figure(1); 

subplot(2,3,1),imshow(I),xlabel('原始图像'); 

subplot(2,3,2),imshow(BW_sobel),xlabel('sobel检测');

subplot(2,3,3),imshow(BW_prewitt),xlabel('prewitt检测'); 

subplot(2,3,4),imshow(BW_roberts),xlabel('roberts检测'); 

subplot(2,3,5),imshow(BW_laplace),xlabel('laplace检测');

subplot(2,3,6),imshow(BW_canny),xlabel('canny检测'); 

%%加入高斯噪声（μ=0，σ^2=0.01）检测结果

I_g1 = imnoise(I,'gaussian',0,0.01);

BW_sobel = edge(I_g1,'sobel'); 

BW_prewitt = edge(I_g1,'prewitt'); 

BW_roberts = edge(I_g1,'roberts'); 

BW_laplace = edge(I_g1,'log'); 

BW_canny = edge(I_g1,'canny'); figure(2); 

subplot(2,3,1),imshow(I_g1),xlabel('加入高斯噪声(μ=0，σ^2=0.01)图像');

subplot(2,3,2),imshow(BW_sobel),xlabel('sobel检测'); 

subplot(2,3,3),imshow(BW_prewitt),xlabel('prewitt检测'); 

subplot(2,3,4),imshow(BW_roberts),xlabel('roberts检测'); 

subplot(2,3,5),imshow(BW_laplace),xlabel('laplace检测'); 

subplot(2,3,6),imshow(BW_canny),xlabel('canny检测');  

%%加入高斯噪声（μ=0，σ^2=0.02）检测结果

I_g2 = imnoise(I,'gaussian',0,0.02); 

BW_sobel = edge(I_g2,'sobel');

BW_prewitt = edge(I_g2,'prewitt'); 

BW_roberts = edge(I_g2,'roberts'); 

BW_laplace = edge(I_g2,'log');

BW_canny = edge(I_g2,'canny'); figure(3); 

subplot(2,3,1),imshow(I_g2),xlabel('加入高斯噪声(μ=0，σ^2=0.02)图像');

subplot(2,3,2),imshow(BW_sobel),xlabel('sobel检测');

subplot(2,3,3),imshow(BW_prewitt),xlabel('prewitt检测'); 

subplot(2,3,4),imshow(BW_roberts),xlabel('roberts检测'); 

subplot(2,3,5),imshow(BW_laplace),xlabel('laplace检测');

subplot(2,3,6),imshow(BW_canny),xlabel('c

一、没有噪声时的检测结果 1 原始图像 

2 Sobel算子边缘检测 3 Prewitt算子边缘检测 4 Roberts算子边缘检测 5 Laplace算子边缘检测 6 Canny算子边缘检测 

二、加入高斯噪声（μ=0，σ^2=0.01）检测结果 1 原始图像 

2 Sobel算子边缘检测 3 Prewitt算子边缘检测 4 Roberts算子边缘检测 5 Laplace算子边缘检测 6 Canny算子边缘检测  

三、加入高斯噪声（μ=0，σ^2=0.02）检测结果 1 原始图像 

2 Sobel算子边缘检测 3 Prewitt算子边缘检测 4 Roberts算子边缘检测 5 Laplace算子边缘检测 6 Canny算子边缘检测   

clear all; close all; 

warning off all;

clc
clear all
close all
I = imread('cameraman.tif'); % 读入图像
imshow(I);title('原图')
BW1 = edge(I,'canny'); % 调用canny函数
figure,imshow(BW1); % 显示分割后的图像，即梯度图像
title('Canny')

这个指的是c源代码数据包
还是简单的 I = imread(''); BW = edge(I,'canny');