clear?all;
clc;
image=imread（‘2222.bmp‘）;
subplot（121）;
imshow（image）;
[m?n]=size（image）;
i_fit=zeros（mn）;%初始化邊緣檢測結果數據存儲矩陣
i_tran=zeros（mn）;%過渡判斷矩陣
thre=100;??%閾值
x=[-1?-2?-1;0?0?0;1?2?1];??%sobel算子模板x方向
y=[-1?0?1;-2?0?2;-1?0?1];%sobel算子模板y方向
for?i=2:m-1%開始對行列進行下面掃描卷積
????for?j=2:n-1
????????i_temp=[?image（i-1j-1）??image（i-1j）??image（i-1j+1）;image（ij-1）??image（ij）??image（ij+1）;...
????????????image（i+1j-1）??image（i+1j）??image（i+1j+1）];
????????x_temp=double（i_temp）.*x;
????????y_temp=double（i_temp）.*y;
????????gx=sum（sum（x_temp））;
????????gy=sum（sum（y_temp））;
????????i_tran（ij）=sqrt（gx^2+gy^2）;??%得到相應的微分值
????????if?i_tran（ij）>=thre;??%進行判斷是否大于該閾值
????????????i_fit（ij）=1;???%大于該閾值記錄為0
????????end
????end
end
subplot（122）;
imshow（i_fit）;

??????

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

?????文件?????180900??2012-07-09?15:12??facet小面模板卷積\2222.bmp

?????文件????????907??2012-07-10?19:21??facet小面模板卷積\template_research.m

?????目錄??????????0??2012-07-10?19:22??facet小面模板卷積

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

???????????????181807????????????????????3