clc;
clear?all;
close?all;
I=imread（‘gray.bmp‘）;??
th1=20;???
Img=double（I）;??
?G=size（I）;??
?alf=3;??
?n=7;?
?n0=floor（（n+1）/2）;
?
?T=zeros（G（1）G（2））;??
?R=zeros（G（1）G（2））;??
?H=zeros（G（1）G（2））;??

?for?i=1:n????%%Gaussian濾波器
?????for?j=1:n??
????????h（ij）=exp（-（（i-n0）^2+（j-n0）^2）/（2*alf））/（2*pi*alf）;??
????end??
?end??
??
?Img_n=uint8（conv2（Imgh‘same‘））;??
?figure（）;?imshow（I）;??????title（‘原圖‘）;??
?figure（）;?imshow（Img_n）;??title（‘和Gaussian卷積后的圖像‘）;?
?imwrite（Img_n‘G.bmp‘）;
??

?Img_n=double（Img_n）;??
?for?i=2:（G（1）-1）??
?????for?j=2:（G（2）-1）??
?????????Sx=Img_n（i-1j-1）+2*Img_n（ij-1）+Img_n（i+1j-1）-Img_n（i-1j+1）-2*Img_n（ij+1）-Img_n（i+1j+1）;??
?????????Sy=Img_n（i+1j-1）+2*Img_n（i+1j）+Img_n（i+1j+1）-Img_n（i-1j-1）-2*Img_n（i-1j）-Img_n（i-1j+1）;??
?????????E（ij）=sqrt（Sx^2+Sy^2）;??
?????????if（E（ij）?????????????E（ij）=0;??
?????????end??
?????end??
?end??
?figure（）?;imshow（uint8（E））;?title（‘Canny算子處理之后的圖像‘）;???
?imwrite（E‘C6.bmp‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-11-29?11:23??MATLAB實現(xiàn)Canny邊緣檢測\
?????文件??????251494??2012-11-24?09:50??MATLAB實現(xiàn)Canny邊緣檢測\C6.bmp
?????文件??????252928??2012-11-29?11:23??MATLAB實現(xiàn)Canny邊緣檢測\Cannywss.doc
?????文件????????1029??2012-11-24?09:50??MATLAB實現(xiàn)Canny邊緣檢測\Cannywss.m
?????文件??????252058??2012-11-24?09:50??MATLAB實現(xiàn)Canny邊緣檢測\G.bmp
?????文件??????252058??2012-11-24?09:50??MATLAB實現(xiàn)Canny邊緣檢測\gray.bmp