8方向sobel边缘梯度平方和图像质量检测

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