eightmodelSobelfunction.zip

eightmodelSobelfunction.zip

%   感谢亲亲使用此代码，此代码解决您的问题了吗~（@^_^@）~
%   没解决的话告诉亲亲一个好消息，我这里可以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