%???感謝親親使用此代碼，此代碼解決您的問題了嗎~（@^_^@）~
%???沒解決的話告訴親親一個好消息，我這里可以1分錢幫助代碼改錯，還提供1分錢成品代碼（′▽‘〃）哦~
%???登錄淘寶店鋪“大成軟件工作室”便可領取
%???是的，親親真的沒有看錯，撓破頭皮的問題真的1分錢就可以解決了\（^o^）/YES!
%???小的這就把傳送門給您，記得要收藏好哦（づ￣3￣）づ╭～
%???傳送門：https://item.taobao.com/item.htm?spm=a1z10.1-c.w4004-15151018122.5.uwGoq5&id=538759553146
%???如果傳送門失效，親親可以來店鋪討要，客服MM等親親來騷擾哦~（*/ω╲*）
%函數sobel8_grad的結果L為清晰度評價值，即為邊緣像素的梯度平方的和
%采用sobel算子八方向模板
function?L=sobel8_grad（f）
f=double（f）;
[mn]=size（f）;
f=[zeros（m1）?f?zeros（m1）];
f=[zeros（1n+2）;f;zeros（1n+2）];
new=f;
U_mean=mean（mean（f））;
for?x=1:m+2
????for?y=1:n+2
????????var1（xy）=（f（xy）-U_mean）^2;
????end
end
var=sum（sum（var1））;
TN=sqrt（var/（m*n））;%閾值門限

for?i=2:m??%i代表行
????for?j=2:n???%j代表列
????????H1=（f（i-1j+1）+2*f（ij+1）+f（i+1j+1））-（f（i-1j-1）+2*f（ij-1）+f（i+1j-1））;
????????H2=（f（i+1j）+2*f（i+1j+1）+f（ij+1））-（f（ij-1）+2*f（i-1j-1）+f（i-1j））;
????????H3=（f（i+1j-1）+2*f（i+1j）+f（i+1j+1））-（f（i-1j-1）+2*f（i-1j）+f（ij+1））;
????????H4=（f（ij-1）+2*f（i+1j-1）+f（i+1j））-（f（i-1j）+2*f（i-1j+1）+f（ij+1））;
????????H5=（f（i-1j-1）+2*f（ij-1）+f（i+1j-1））-（f（i-1j+1）+2*f（ij+1）+f（i+1j+1））;
????????H6=（f（ij-1）+2*f（i-1j-1）+f（i-1j））-（f（i+1j）+2*f（i+1j+1）+f（ij+1））;
????????H7=（f（i-1j-1）+2*f（i-1j）+f（ij+1））-（f（i+1j-1）+2*f（i+1j）+f（i+1j+1））;
????????H8=（f（i-1j）+2*f（i-1j+1）+f（ij+1））-（f（ij-1）+2*f（i+1j-1）+f（i+1j））;%提取8個方向上的邊緣成分
????????g（ij）=H1^2+H2^2+H3^2+H4^2+H5^2+H6^2+H7^2+H8^2;%計算像素點的梯度值平方
????????delta_g（ij）=sqrt（g（ij））;?%對g（ij）開方得到像素點的梯度值
????????if?delta_g（ij）>TN??%梯度值大于閾值，則該點f（ij）為邊緣點，把3*3窗口中央像素灰度值用255代替。0~255即從黑到白
????????????new（ij）=255;
????????????gl（ij）=g（ij）;
????????else
?????????new（ij）=0;
?????????gl（ij）=0;
????????end
????end?
end

?L=sum（sum（gl））;%得出評價函數L的值所有邊緣像素的梯度平方的和

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-01-11?20:11??eight?model?Sobel?function\
?????文件?????1179702??2003-03-20?20:10??eight?model?Sobel?function\1.bmp
?????文件?????1179702??2003-03-20?20:08??eight?model?Sobel?function\2.bmp
?????文件?????1179702??2003-03-20?20:09??eight?model?Sobel?function\3.bmp
?????文件?????1179702??2003-03-20?20:11??eight?model?Sobel?function\4.bmp
?????文件?????1179702??2003-03-20?20:10??eight?model?Sobel?function\5.bmp
?????文件?????1179702??2002-04-05?18:23??eight?model?Sobel?function\6.bmp
?????文件?????1179702??2003-03-18?22:35??eight?model?Sobel?function\7.bmp
?????文件????????2073??2016-09-23?09:46??eight?model?Sobel?function\sobel8_grad.m
?????文件????????1799??2010-12-14?13:19??eight?model?Sobel?function\sobel8_grad_test.asv
?????文件????????2372??2016-09-23?09:47??eight?model?Sobel?function\sobel8_grad_test.m
?????文件???????28160??2011-01-11?20:11??eight?model?Sobel?function\Thumbs.db