写这个难道不许么
求一下源程序
（1） 读入原始图像，将原始图像f(x,y)经过LOG滤波器模板的卷积，得到输出h1，用矩阵a[i][j]表示。
（2） 对于输出图像a[i,j]中每一个元素a[i1][j1]选取负值点，再取出邻点的最大值，输出为最大值与其本身值之差，其余输出为0，得到梯度图像h2 ，用矩阵b[i][j]表示。
（3） 分析LOG滤波器后的梯度图像，根据该图像的直方图统计特性确定阈值 。本算法中选取的最大值max和最小值min。则
         t1=min+(max-min)/6
              t2=min+(max-min)/12
若b[i][j]的某个图像元素的梯度值大于t1，则它被输出为255,；所有超出tt2的像素均被标记为128，若一个已被标记的像素与某个已输出的像素连通，则它也被输出255.