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度角图像边缘检测') ;

%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('水平和垂直边缘检测');