一种基于LOG算子的图像边缘提取的方法与系统

摘要

本申请公开了一种基于LOG算子的图像边缘提取的方法与系统，该方法包括：步骤1，读取原始图像A，获得图像的最大行数m及最大列数n；步骤2，图像中间像素则进行如下3x3邻域LOG模板循环计算，得到初始输出；步骤3对初始输出进行二值化处理，绝对值高于事先给定的参数a的赋值1，低于a的赋值0，得到二值化的初始边缘图像；步骤4对上步赋值为1的点搜索联结像素数，低于事先给定的第二个参数b的判定为噪点，从初始边缘图中删除，高于b的判定为真边缘点；步骤5得到最终的增强图像。本申请的技术方案获得比较准确清晰的图像边缘。

说明书

技术领域

本发明涉及图像处理领域，尤其涉及图像边缘提取的方法及系统。

背景技术

在数字图像处理技术中，图像边缘提取是图像处理极为重要的组成部部分。同时在生物医疗、计算机视觉、机器学习、模式识别等诸多领域中，边缘增强都是一个重要的研究课题。

传统的边缘提取算法，都是先利用微分算子求出每个像素的高通分量然后过滤得到图像边缘，经典的微分算子有prewitt，sobel，Roberts，laplace，log等

发明内容

本申请提供了基于log算子的改进图像边缘提取方法与系统。

类似于一阶微分算子，本方法采用log算子来求初步的图像边缘，在3x3模块中：

A(i，j)＝24＊I(i，j)+8＊(I(i-1，j)+I(i+1，j)+I(i，j-1)+I(i，j+1))-4＊(I(i-2，j-1)+I(i-2，j)+I(i-2，j+1)+I(i+2，j-1)+I(i+2，j)+I(i+2，j+1)+I(i-1，j-2)+I(i，j-2)+I(i+1，j-2)+I(i-1，j+2)+I(i，j+2)+I(i+1，j+2))-2＊(I(i-2，j-2)+I(i-2，j+2)+I(i+2，j-2)+I(i+2，j+2))；

其中I(i，j)为中间像素，A(i，j)为I(i，j)的输出值。。

在得到初步边缘后，对低于阈值的像素点置零，高于阈值a置1，得到黑白的边缘图V。

对得到的二值化边缘图V中为1的像素点搜索连接像素数，连接像素数低于阈值b，判定为噪点从初始边缘图中删除；高于b的判定为真边缘点则保留。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1图法流程图 图2是需要提取边缘的原图像，图3是用log算法得到的边缘图，图4是用本申请的方法所得到的边缘图。

具体实施方式

如在说明书及权利要求当中使用了某些词汇来指称特定组件。本领域技术人员应可理解，硬件制造商可能会用不同名词来称呼同一个组件。本说明书及权利要求并不以名称的差异来作为区分组件的方式，而是以组件在功能上的差异来作为区分的准则。如在通篇说明书及权利要求当中所提及的“包含”为一开放式用语，故应解释成“包含但不限定于”。“大致”是指在可接收的误差范围内，本领域技术人员能够在一定误差范围内解决所述技术问题，基本达到所述技术效果。说明书后续描述为实施本发明的较佳实施方式，然所述描述乃以说明本发明的一般原则为目的，并非用以限定本发明的范围。本发明的保护范围当视所附权利要求所界定者为准。

本申请提供的图像处理方法包括：

步骤1，读取原始图像A，获得图像的最大行数m及最大列数n；

步骤2对中间像素逐行逐列进行如下操作：

A(i，j)＝24＊I(i，j)+8＊(I(i-1，j)+I(i+1，j)+I(i，j-1)+I(i，j+1))-4＊(I(i-2，j-1)+I(i-2，j)+I(i-2，j+1)+I(i+2，j-1)+I(i+2，j)+I(i+2，j+1)+I(i-1，j-2)+I(i，j-2)+I(i+1，j-2)+I(i-1，j+2)+I(i，j+2)+I(i+1，j+2))-2＊(I(i-2，j-2)+I(i-2，j+2)+I(i+2，j-2)+I(i+2，j+2))；得到初步的输出；

步骤5对初步输出进行二值化处理，低于指定阈值的像素点置零，高于指定阈值的像素点置，得到黑白边缘图；

步骤6对黑白边缘图中值为1的像素搜索联结像素数，连接像素数低于阈值b，判定为噪点，从初始黑白边缘图中删除；高于b的判定为真边缘点则保留。

步骤7得到最终的边缘图。

上述流程中，其中A(i，j)表示原图像中第i行j列的像素，V(i，j)表示转换后的第i行j列的像素。

上式等号左侧表示第i行j列像素的邻域操作的结果，右侧表示邻域操作的具体计算方式，其中V(i，j)表示当前第i行j列的像素，其表示邻域操作之前像素的值，A(i，j)的邻域有八个，具体指A(i-1，j-1)，A(i-1，j)，A(i-1，j+1)，A(i，j-1)，A(i，j+1)，A(i+1，j-1)，A(i+1，j)，A(i+1，j+1)等。

过程中涉及到两个常数，二值化时用到的阈值a，搜索联结像素数时用到的阈值b依据个人经验及具体图像指定。

需要说明的是，上面的程序是基于matlab实现仅作为示例，本领域技术人员可以利用其他编程语言来实现本发明的技术方案；

本领域内的技术人员应明白，本发明的实施例可提供为方法、装置、或计算机程序产品。因此，本发明可采用完全硬件实施、完全软件实施、或结合软件和硬件方面的实施的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

上述说明示出并描述了本发明的若干具体实施例，但如前所述，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述发明构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。