Sobel算子实现水平边缘检测、垂直边缘检测；45度、135度角边缘检测

%File Discription:

%45°和135°角边缘检测；用于那些边界不明显的图片

%不太适用于复杂图，复杂图用水平和垂直边缘检测

%Author:Zhang Ruiqing

%CreateTime:2011.8.8（What a good day!(*^__^*) ）

SourcePic=imread('D:\Images\pic_loc\1870378220205041520.jpg');

subplot(221);

imshow(SourcePic),title('原图');

grayPic=rgb2gray(SourcePic);

grayPic=im2double(grayPic);

%使用指定45度角Sobel算子滤波器，指定阂值

a45=[-2 -1 0;-1 0 1;0 1 2];

SFST45=imfilter(grayPic,a45,'replicate');%功能：对任意类型数组或多维图像进行滤波。

SFST45=SFST45>=Threshold;

subplot(222);

imshow(SFST45),title('45度角图像边缘检测') ;

b45=[0 -1 -2;1 0 -1;2 1 0];

SFST45=imfilter(grayPic,b45,'replicate');%功能：对任意类型数组或多维图像进行滤波。

SFST45=SFST45>=Threshold;

SFST45

subplot(223);

imshow(SFST45),title('135度角图像边缘检测') ;

[html] view
plaincopy

%File Discription:

%水平、垂直边缘检测；

%Author:Zhang Ruiqing

%CreateTime:2011.8.8

SourcePic=imread('D:\毕业设计\Images\pic_loc\1870399350205061354.jpg');

subplot(221);

imshow(SourcePic),title('原图');

grayPic=rgb2gray(SourcePic);

grayPic=im2double(grayPic);

[Vertical Threshold]=edge(grayPic,'sobel','vertical');%edge detect

subplot(222);

imshow(Vertical),title('垂直方向边缘检测');

[Horizontal Threshold]=edge(grayPic,'sobel','horizontal');%edge detect

subplot(223);

imshow(Horizontal),title('水平方向边缘检测');

H_V=edge(grayPic,'sobel',Threshold);

subplot(224);

imshow(H_V),title('水平和垂直边缘检测');