%函數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=（