canny边缘检测的matlab程序

clc clear close all RGB=imread('lena.jpg'); figure; subplot(2,3,1) imshow(RGB); title('彩色原图'); a=rgb2gray(RGB); subplot(2,3,2) imshow(a); title('灰度图'); J=imnoise(a,'salt',0.04);%添加椒盐 subplot(2,3,3) imshow(J); title('椒盐噪声'); h0=1/9.*[1,1,1,1,1,1,1,1,1]; h1=1/10.*[1,1,1,1,2,1,1,1,1]; h2=1/16.*[1,2,1,2,4,2,1,2,1]; g0=imfilter(J,h0,'replicate'); subplot(2,3,4) imshow(g0); title('邻域平均法滤波'); g1=imfilter(J,h2,'replicate'); subplot(2,3,5) imshow(g1); title('加权平均法滤波'); k=medfilt2(J); subplot(2,3,6) imshow(k); title('中值滤波'); E=edge(k,'canny'); figure(2); imshow(E); title('canny算子检测')