%?DOG?檢測焦點的算法
%?令狐雄坤
%?DOG算子特征提取算法講解參考網站1：http://www.it610.com/article/4975225.htm
%?代碼修改參考網站2：??????????????https://blog.csdn.net/songzitea/article/details/8831909
????clc
????close?all
????clear?all
????src_old?=?imread（‘C:\Users\linghuxiongkun\p\flower.jpg‘）;
????src?=?rgb2gray（src_old）;
????src?=?double（src）;?
????h1?=?fspecial（‘gaussian‘?5?0.5）;??%利用自帶的函數設計高斯濾波器
????h2?=?fspecial（‘gaussian‘?5?0.6）;
????h3?=?fspecial（‘gaussian‘?5?0.7）;
????h4?=?fspecial（‘gaussian‘?5?0.8）;

????src1?=?imfilter（src?h1）;??%高斯模板對原圖像進行濾波
????src2?=?imfilter（src?h2）;
????src3?=?imfilter（src?h3）;
????src4?=?imfilter（src?h4）;

????%?角點檢測算法
????src12?=?src1?-?src2;
????src23?=?src3?-?src2;
????src34?=?src4?-?src3;

????src_DoG（::1）?=?src12;
????src_DoG（::2）?=?src23;
????src_DoG（::3）?=?src34;

????figureimshow（src12）;
????title（‘src12‘）;
????figureimshow（src23）;
????title（‘src23‘）;
????figureimshow（src34）;
????title（‘src34‘）;
????
????result?=?zeros（size（src_DoG?1）?size（src_DoG?2））;

????level?=?5;
????around?=?zeros（19）;
????for?i?=?2:size（src_DoG?1）-1
????????for?j?=?2:size（src_DoG?2）-1??%遍歷查找目標點相鄰26個點的高斯差值，
??????????????????????????????????????%判斷目標點的高斯差值是否為相鄰26個點總的最大值最小值
????????????near_i