%%%%%%%產生sobel算子%%%%%%
%%%%%%%by?mqzz
%%%梯度算子一階導數一階差分有方向的故有水平和垂直兩個模板
%%%此程序缺點:對邊緣和背景做了最粗糙處理還可以細化.
clc
clear
X=imread（‘original_lena.bmp‘）;
figureimshow（X）;title（‘原始圖‘）;
[ROWCOL]=size（X）;
%%%這步是有用的保證W行列與原始一致.
W=zeros（ROWCOL）;
X=double（X）;
%%%從第2行第2列開始到倒數第2行倒數第2列
for?i=2:ROW-1
????for?j=2:COL-1
????????%%%%Dx是檢測水平邊緣的則它的梯度方向是垂直的這是對X的偏導數.
????????Dx=X（i+1j-1）+2*X（i+1j）+X（i+1j+1）-X（i-1j-1）-2*X（i-1j）-X（i-1j+1）;
????????%%%%Dy是檢測垂直邊緣的則它的梯度方向是水平的這是對Y的偏導數.
????????Dy=X（i-1j+1）+2*X（ij+1）+X（

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

?????文件???????1268??2008-12-22?20:50??mqzz_sobel_dlg.m

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

?????????????????1268????????????????????1