一种基于三维块匹配的改进自适应阈值Canny边缘检测方法

技术领域

本发明涉及一种在图像处理领域内使用的边缘检测方法，尤其涉及一种基于三维块匹配的改进Canny检测方法。

背景技术

图像的边缘特征是图像最基本的特征之一，包含图像丰富的固有信息。图像处理在各个领域中被广泛使用，在进行图像识别、图像分割、目标测量等处理前需要对图像的边缘进行准确的检测。

对于分辨率低、信噪比低、边缘模糊、图像不同区域之间的边界不明确等特点的图像，传统Canny算法在检测时存在一定的局限。传统Canny使用高斯滤波对图像进行平滑时容易在平滑噪声的同时也滤掉了有用的信息从而丢失大量边缘信息。而且传统Canny算子的双阈值需要人工事先选定，该步骤的不确定性会导致检测中准确度的降低。为此，国内相关研究提出了许多改进Canny边缘检测算法。比如使用统计滤波、非均值滤波、双边滤波、混合滤波器等算法替换高斯滤波，或改用具有保边特性的引导滤波对图像进行平滑，采用基于形态学算子的各向异性扩散滤波代替高斯滤波平滑，利用 DCT 系数进行线性收缩对图像进行去噪。并且改进双阈值处理，使用大律法自适应地根据图像灰度选取高、低阈值，采用梯度强度信息计算法实现阈值自适应，对非极大值抑制后的梯度图像进行分块根据各子块的类型自适应确定边缘连接的双阈值，对幅度梯度直方图进行微分运算得到自适应双阈值。

虽然目前的一些改进Canny算法中采用了一些边缘保持技术，但边缘的保持和增强方法仍然有改进的余地，在高低双阈值处理中，大部分是整幅图像采用统一的双阈值，没有考虑图像的局部梯度分布特征，缺乏局部自适应性，很容易造成边缘信息的丢失。因此，本发明在基于传统Canny算法的基础上，提出了一种基于三维块匹配的改进自适应阈值Canny边缘检测方法。

发明内容

本发明针对待解决的技术问题提供了一种基于三维块匹配的改进自适应阈值Canny边缘检测方法。

为达到上述目的，本发明采用如下技术方案：

一种基于三维块匹配的改进自适应阈值Canny边缘检测方法，其特征在于：包括以下步骤：步骤1、输入需要检测的图像；步骤2、对图像进行三维块匹配滤波；步骤3、对三维块匹配滤波图像进行引导滤波；步骤4、梯度幅值以及方向计算；步骤5、对引导滤波梯度图进行非极大值抑制；步骤6、块匹配组内最大类间方差计算确定最佳双阈值；步骤7、根据双阈值保留或剔除边缘；步骤8、基于区域生长法进行边缘连接得到图像的边缘。

作为优选，输入需要检测的图像的方法是：输入需要检测的图像是事先采集并保存的图像数据或实时采集的图像数据。

作为优选，对图像进行三维块匹配滤波的方法是：根据搜索步长遍历原始图像中的所有图像块，在搜索范围内寻找与其相似度大于阈值的图像块归为同组，对同组图像块进行协同硬阈值滤波然后加权重组得到基础估计图像；然后根据搜索步长遍历基础估计图像中的所有图像块，在搜索范围内寻找与其相似度大于阈值的图像块归为同组，将同组图像块的基础图像块组和相同位置的原始图像块组进行协同维纳滤波后加权重组得到最终估计图像。

作为优选，对三维块匹配滤波图像进行引导滤波的方法是：由硬阈值协同滤波后的基础估计图像作为引导图像，维纳协同滤波后的最终估计图像作为处理图像，使用局部线性模型建立引导图像G与待处理图像P之间的关系进行引导滤波，输出图像的梯度信息完全由引导图像的梯度信息决定。

作为优选，梯度幅值以及方向计算的方法是：采用四方向的Sobel算子计算基础估计图像和引导滤波图像的梯度幅值以及方向。

作为优选，对引导滤波梯度图像进行非极大值抑制的方法是：遍历引导滤波梯度图像，根据x和y方向的梯度分量大小确定非极大值抑制的方向，比较当前像素点在抑制方向上是否为极大值，非极大值置零，为极大值则保留，保证不出现双边缘的情况，当在抑制方向上不存在像素点时通过比例进行插值计算。

作为优选，块匹配组内最大类间方差计算确定最佳双阈值的方法是：将最终估计图像的梯度图像依据原先三维块匹配中的图像块分组情况分组，使用最大类间方差计算图像块组内的一个阈值，根据这个阈值确定在双阈值处理时的最佳双阈值。

作为优选，根据双阈值保留或剔除边缘的方法是：遍历非极大值抑制图像的图像块按照最大类间方差计算出的最佳双阈值进行处理，高于高阈值的像素视为强边缘像素值置为255，低于低阈值的像素视为不是边缘置为0，高于低阈值且低于高阈值的像素视为弱边缘，判断其周围像素的像素值是否高于高阈值，如果高于高阈值则置为255，视为边缘。

作为优选，基于区域生长法进行边缘连接得到图像边缘的方法是：计算双阈值连接边缘图与低阈值分割图的差异图，对差异图中的每个非零像素点，在其为中心的5 × 5邻域内根据区域生长法寻找与其相连的边缘点，最后判断当前点是否为连续性较好的弱边缘点。若与之相连的边缘点个数大于等于 4，则判断该中心点是弱边缘点并保留，否则判为噪声点将其去除。

本发明的有益效果是：本发明所公开的一种基于三维块匹配的改进自适应阈值Canny边缘检测方法，通过三维块匹配与引导滤波，在过滤噪声信息的同时对边缘信息进行保留，避免了高斯滤波将边缘信息与噪声信息同时过滤的问题；同时根据三维块匹配中的匹配结果引入梯度图像块组最大类间方差以自适应确定双阈值，解决了人工设定阈值效率低、检测结果可靠性差的问题；对每个图像块进行双阈值处理，避免了全局阈值处理漏检弱边缘的问题，适用于图像处理中各类图像的边缘检测，在检测精度方面具有较大优势。

附图说明

图1是本发明的改进自适应阈值边缘检测方法整体流程。

图2是本发明的边缘检测具体流程

图3是本发明处理的原始图像。

图4是本发明图像块匹配示意图。

图5是本发明三维块匹配滤波的基础估计图像。

图6是本发明三维块匹配滤波的最终估计图像。

图7是本发明引导滤波的图像。

图8是本发明最终估计图像的梯度图像。

图9是本发明引导滤波图像的梯度图像。

图10是本发明阈值分布图像。

图11是本发明非极大值抑制图像。

图12是本发明的最终边缘检测结果图像。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明，实施例对图像进行边缘检测。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。

实施例，请参照图1至图12：图1示意性地给出了本发明所公开的改进自适应阈值边缘检测方法整体流程；图2示意性给出了本发明边缘检测具体流程；图3示意性给出了本发明处理的原始图像；图4示意性给出了本发明图像块匹配示意图；图5示意性给出了本发明三维块匹配滤波的基础估计图像；图6示意性给出了本发明三维块匹配滤波的最终估计图像；图7示意性给出了本发明引导滤波的图像；图8示意性给出了本发明最终估计图像的梯度图像；图9示意性给出了本发明引导滤波图像的梯度图像；图10示意性给出了本发明非极大值抑制图像；图11示意性给出了本发明计算图像块最佳阈值后的阈值分布图像；图12示意性给出了本发明的最终边缘检测结果图。

依据图1所示流程，对图像边缘检测依据以下步骤进行检测：

（1）通过工业相机实时采集图像或读取事先保存的图像数据，得到原始图像

（2）为了过滤掉图像中存在的噪声，使用三维块匹配滤波算法，整个算法分为基础估计和最终估计两部分。

基础估计首先将原始图像

其中，

然后将参考块和相似块叠加在一起聚集成相应的三维图像块组，对形成的三维图像块组进行第三维的小波硬阈值变换，通过对变换域系数进行硬阈值收缩处理来大幅度降低噪声，其滤波函数定义如式（2）所示：

其中，

随后通过三维反变换来产生所有聚集块的估计值，然后返回这些块的估计值到他们的初始位置。由于每一图像块返回到原始位置的估计值是重叠或非重叠的，于是对所有重叠的块的估计值进行加权平均，得到重叠块的单个的估计值，从而重构出图像的基础估计图像

与基础估计步骤相似，最终估计首先进行选定相似块进行分组。不同的是该步骤先对基础估计得到的图像进行分块，通过欧氏距离搜索到相似块的位置后分组得到一个来自基础估计图像

其中，

随后通过三维反变换可以得到每个相似块的估计值，然后返回这些估计值到他们的初始位置进行加权平均，从而得到重构出图像的最终估计图像

（3）由于三维块匹配滤波对每个块的二维变换会引入假信号，当噪声强度相对较大时尤为明显，所以使用引导滤波对最终估计进行滤波，以基础估计图像作为参考图像在保留特征的前提下去除假信号，得到引导滤波图像

引导滤波算法原理是一种线性滤波过程，可以简单定义为式（4）形式:

其中，Q是输出图像，P是输入的待滤波图像，I是引导图像，W是根据引导图I确定的权重值，权重值W可以用式（5）表示:

其中

自适应权重可以根据上式分析得到：

（4）为了计算梯度幅值以及方向，采用四方向的Sobel算子计算，如式（6）所示。最终估计图像的梯度幅值图

（5）为了细化边缘，对引导滤波梯度图像

（6）为了自适应计算双阈值处理过程中的高低阈值，采用块匹配组内最大类间方差计算的方法确定最佳双阈值，首先将最终估计图像的梯度图像

（7）为了保留更多有效边缘进行双阈值处理，遍历非极大值抑制图像

（8）为了保留某些连接性较好只是没有与大于高阈值的像素点相邻接的弱边缘点，计算双阈值连接边缘图与低阈值分割图的差异图，对差异图中的每个非零像素点，在其为中心的5 × 5邻域内根据区域生长法寻找与其相连的边缘点，最后判断当前点是否为连续性较好的弱边缘点。若与之相连的边缘点个数大于等于 4，则判断该中心点是弱边缘点并保留，否则判为噪声点将其去除。得到最终边缘检测结果图哈哈

上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。