function?[out_args]=EdgeCheck（fnchoice）
bmparr=mat2gray（double（imread（fn））[0255]）;
[xlen?ylen]=size（bmparr）;
workarr=zeros（xlenylen）;
%subplot（121）;
%figure（1）;
%imshow（bmparr）;
%title（‘原圖像‘）;
switch?choice
????case?‘sobel‘
????????dx=zeros（xlenylen）;
????????dy=zeros（xlenylen）;
????????dxy=zeros（xlenylen）;
????????h1=[-1?0?1;-2?0?2;-1?0?1]/8;
????????h2=[-1?-2?-1;0?0?0;1?2?1]/8;
????????for?i=1:xlen
????????????for?j=1:ylen
????????????????for?k=1:3
????????????????????for?l=1:3
????????????????????????xp=i-k+2;
????????????????????????yp=j-l+2;
????????????????????????if（（xp>=1）&&（xp<=xlen）&&（yp>=1）&&（yp<=ylen））
????????????????????????????dx（ij）=dx（ij）+h1（kl）*bmparr（xpyp）;
????????????????????????????dy（ij）=dy（ij）+h2（kl）*bmparr（xpyp）;
???????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????8994??2009-07-02?11:42??EdgeCheck.m

-----------?---------??----------?-----??----

?????????????????8994????????????????????1