基于各向异性高斯方向导数滤波器提取图像粗边缘的方法

基于各向异性高斯方向导数滤波器提取图像粗边缘的方法
【专利摘要】本发明公开了基于各向异性高斯方向导数滤波器提取粗边缘的方法，其采用离散化各向异性高斯核方向导数滤波器对图像进行平滑，然后计算平滑后的图像像素的最大梯度方向对应的梯度幅值，最后利用梯度相关性提取粗边缘。本发明利用各向异性高斯方向导数滤波器处理图像对噪声有很好的抑制作用且可以很好的提取各个方向上的灰度变化信息，克服了由于高斯尺度因子带来的边缘轮廓移位及噪声敏感的问题，为基于边缘轮廓的角点检测算法提供了一种新的边缘提取方法。
【专利说明】基于各向异性高斯方向导数滤波器提取图像粗边缘的方法

【技术领域】
[0001] 本发明属于图像特征提取【技术领域】，具体涉及基于各向异性高斯方向导数滤波器 提取图像粗边缘的方法。

【背景技术】
[0002] 图像的特征（如：边缘、角点等）包含了丰富的信息，可以对图像处理提供足够的 约束。提取图像特征信息可以有以下的优势：（1)对图像特征进行处理大大减少了运算的 数据量，（2)特征具有旋转、放大、缩小和平移不变的性质。提取边缘的技术常常被利用到 角点检测应用中，现今的角点检测大致可以分为三类：基于灰度强度变化检测方法、基于边 缘轮廓的检测方法，基于角点模型匹配的检测方法。基于边缘轮廓的检测方法，关键技术在 于对边缘轮廓的提取。
[0003] 基于边缘轮廓的角点检测算法主要是通过分析图像边缘的形状特性来检测角点， 这种方法首先提取图像的边缘轮廓，然后从边缘中提取封闭或者非封闭的边缘轮廓曲线， 最后通过搜索轮廓线上的曲率局部极大值点或者斜率或梯度方向变化点搜索交点来提取 角点。最具代表的检测方法为曲率尺度空间CSS技术，此检测算法由于平滑边缘轮廓的高 斯尺度值难以确定，会产生边缘轮廓移位的现象，使得角点定位性能降低。
[0004] 大多数基于边缘轮廓的角点检测算法，一般选择canny检测器提取边缘，而选用 这类方法提取到的边缘会存在断裂的现象。同时利用各向同性高斯核会产生由尺度因子的 选择造成的边缘轮廓移位的现象及对噪声敏感的问题。由于角点处多方向呈现低自相似性 的特性，所以需要更有效的方法提取边缘轮廓，克服边缘轮廓移位及对噪声敏感的问题。


【发明内容】

[0005] 本发明的目的是提出一种基于各向异性高斯方向导数滤波器提取图像粗边缘的 方法，解决了现有技术中采用各向同性高斯核提取图像边缘，由于尺度选择造成的边缘移 位及对噪声敏感的技术问题。
[0006] 本发明所采用的技术方案，基于各向异性高斯方向导数滤波器提取图像粗边缘的 方法，包括以下步骤：
[0007] 步骤1:将图像读入计算机；
[0008] 步骤2 :在[0, 31 ]内平均方向角上采用各向异性高斯方向导数滤波器平滑图像；
[0009] 步骤3:计算平滑后图像像素的梯度幅值EMGxi (η);
[0010] 步骤4 :计算图像像素的全局标准化梯度幅值，具体步骤如下：
[0011] 4. 1，首先计算图像像素的领域内梯度标准化幅值；
[0012] 以一像素点为中心像素点，选取领域，计算所述领域内中心像素点的领域内梯度 标准化幅值O)，计算公式如下：

【权利要求】
1. 基于各向异性高斯方向导数滤波器提取图像粗边缘的方法，其特征在于，包括以下 步骤： 步骤1 :将图像读入计算机； 步骤2 :在[0, π ]内平均方向角上采用各向异性高斯方向导数滤波器平滑图像； 步骤3 :计算平滑后图像像素的梯度幅值EMGxi (η); 步骤4 :计算图像像素的全局标准化梯度幅值，具体步骤如下： 4. 1，首先计算图像像素的领域内梯度标准化幅值； 以一像素点为中心像素点，选取领域，计算所述领域内中心像素点的领域内梯度标准 化幅值^11)，计算公式如下：
其中，Nq表示围绕中心像素点的领域内像素点的平均梯度幅值，N表示所取领域内像素 点的个数，EMGxi (η)表示领域内第i个像素点的梯度幅值，EMGx(q)表示中心像素点的梯度 幅值； 4. 2,计算图像像素的全局标准化梯度幅值； 计算所述中心像素点的全局标准化梯度幅值，计算公式如下：
其中，Iv^ll代表全局标准化梯度幅值，N1表示整幅图像像素点的个数，U1是整幅图像梯 度均值，σ i是整幅图像的梯度标准偏差； 重复步骤4. 1及步骤4. 2,直至计算出图像所有像素点的全局标准化梯度幅值； 步骤5 :双阈值法提取粗边缘，具体步骤如下： 设置双阈值Th，T1，全局标准化梯度幅值||ν^|大于Th的像素点标记为边缘点，全局标准 化梯度幅值||w I小于T1的像素点视为背景点删除，全局标准化梯度幅值|w I介于Th，T1的 像素点分为两种情况：若其8邻域中存在已经标记的边缘点，则该像素点标记为边缘点，否 则视为非边缘点删除； 步骤6 :输出图像的粗边缘。
2. 根据权利要求1所述的基于各向异性高斯方向导数滤波器提取图像粗边缘的方法， 其特征在于，所述的步骤2的具体步骤如下： 2.1，计算各向异性高斯方向导数； 各向异性高斯核函数的表达式为：
其中，P > 1，是各向异性因子，σ是尺度因子，Θ是旋转角度，R0是旋转矩阵，X是输 入的图像矩阵，Xt即输入图像矩阵X的转置矩阵； 对旋转角度Θ求偏导，得到各向异性高斯核函数的各向异性高斯方向导数 g'。_ 〇 _ θ (X):
2. 2,平滑图像，具体步骤如下： 2. 2. 1对各向异性高斯方向导数滤波器进行离散化处理； 对连续函数￥^0〇〇在二维整数空间上对角度9进行1(个方向上的采样，得到离 散化各向异性高斯方向导数滤波器，其离散表达式如下式：
K是方向米样 数，9k是角度采样间隔； 2. 2. 2,将离散化各向异性高斯方向导数与输入图像做卷积运算，即平滑图像，并得到 的各向异性高斯方向导数响应，如下式：
3. 根据权利要求2所述的基于各向异性高斯方向导数滤波器提取图像粗边缘的方法， 其特征在于，所述的各向异性高斯方向导数滤波器满足σ > p。
4. 根据权利要求1或2所述的基于各向异性高斯方向导数滤波器提取图像粗边缘的方 法，其特征在于，所述的步骤2中选用多个各向异性高斯方向导数滤波器平滑图像，多个各 向异性高斯方向导数滤波器在[〇, η]方向角上均匀分布。
5. 根据权利要求1所述的基于各向异性高斯方向导数滤波器提取图像粗边缘的方法， 其特征在于，所述的步骤3中的平滑后图像像素的梯度幅值EMG Xi(n)的计算公式如下：
【文档编号】G06T7/00GK104376564SQ201410686745
【公开日】2015年2月25日 申请日期:2014年11月24日 优先权日:2014年11月24日
【发明者】章为川, 孔祥楠, 杨婷婷 申请人:西安工程大学