<strong><span style="font-size:18px;"><span style="color:#6600cc;">clear all;
close all;
clc;
img=rgb2gray(imread('9.jpg'));
subplot(241);imshow(img);title('原图')
[m n]=size(img);
img=double(img);
%%canny边缘检测的前两步相对不复杂,所以我就直接调用系统函数了
%%高斯滤波
w=fspecial('gaussian',[5 5]);
img=imfilter(img,w,'replicate');
subplot(242);imshow(uint8(img));title('高斯滤波')
%%sobel边缘检测
w=fspecial('sobel');
img_w=imfilter(img,w,