一种基于Legendre矩的图像亚像素边缘检测方法
【技术领域】
[0001] 本发明涉及图像边缘检测领域，尤其涉及一种基于Legendre矩的图像亚像素边 缘检测方法。
【背景技术】
[0002] 工程中常基于图像实现非接触的几何尺寸精密测量，该方式以其非接触、全视场、 高精度的特点而获得广泛应用。其原理就是通过处理被测物体图像的边缘而获得图像的 几何参数。由此可见，边缘检测是图像测量的基础和关键。传统边缘检测方法多是基于图 像像素灰度的变化，如Sobel算子、Laplacian算子和canny算子等。这些方法形式简单、 易于实现、但定位精度不高，通常只有整数像素级的精度，且微分算子对噪声非常敏感，常 会产生一些伪边缘。随着人们对检测精度要求的不断提高，像素级检测精度已经不能满足 实际测量的要求。为了解决这个问题，人们提出了亚像素边缘检测方法。这些方法可以突 破摄像机物理分辨率的限制，使图像的边缘定位精度达到亚像素级，从而提高了图像测量 系统的检测精度。当算法的精度为〇. 1个像素，则相当于检测系统的硬件分辨率提高了 10 倍。目前亚像素边缘检测方法，在数学上可以归纳为插值法、拟合法、矩方法三种类型。拟 合法通过对给定的边缘模型将图像中的灰度值进行拟合，这种方法具有很高的精度但是耗 时，插值法通过对实际图像的灰度分布进行插值来得到亚像素的位置，但是对噪声很敏感。 矩方法使用了对噪声不敏感的积分算子，因此应用最广泛。
[0003]在文献[1]"SubpixeledgelocationbasedonorthogonalFourier- Mellin moments"中，Bin提出了基于0FMM矩的亚像素边缘检测，该技术利用傅里叶-梅林矩对图 像进行亚像素级的定位，采用了 5X5的模板，求得亚像素坐标的四个参数，巾、1、k、h，然后 对h进行判断，若大于阈值T，则判断该点为边缘点。该技术通过使用0FMM矩进行亚像素边 缘检测，虽然该技术对噪声不敏感，克服了噪声的影响，但是由于需要三个实数模板，两个 复数模板，计算复杂度大，影响求解速度。
【发明内容】
[0004] 根据现有技术存在的问题，本发明公开了一种基于Legendre矩的图像亚像素边 缘检测方法，包括以下步骤：
[0005]S1:读取图像信息，将图像灰度化并对灰度图像进行去噪处理；
[0006]S2:采用Sobel算子对去噪后的图像进行像素级边缘定位：利用像素点的各向邻 点灰度加权值在边缘点达到最大值这一现象进行边缘检测；
[0007]S3:采用Legendre矩对图像进行亚像素边缘检测，输出边缘图像。
[0008]S2中具体采用如下方式：遍历原始灰度图像中所有的像素点，计算得到每个 像素点的梯度值G[f' (x，y)]，将所得的梯度值归一化到[0,255]区间，采用最大类间 方差法计算得到归一化梯度值的阈值T，对每个像素点归一化的梯度值进行判断，即当 G[f' (x，y) ] >T时，对应的像素点设定为255,否则设定为0至此得到图像的像素级粗定位。
[0009] 进一步的，在得到像素级粗定位后遍历图像中所有的边缘点，进行判断：若该边缘 点是孤立的边缘点即在以该点为中心的3X3的邻域中、除本点外的边缘点的个数小于等 于1，则将该点除去，即该点不作为边缘点、判断为噪声。
[0010] 进一步的，S3中具体采用如下方式：遍历检测到的所有边缘点，对每个边缘点进 行如下处理：以得到的边缘点为中心，在灰度图像中选取NXN的窗口，N为奇数，采用如下 公式（25)将NXN灰度图像窗口中的值与Legendre正交矩的掩码CLMn对应位置的系数相 乘得到NXN的矩阵，将该矩阵求和得到Legendre正交矩LMn，同样的方式利用式（26)再 求得一个Legendre正交矩LM31，
[0013] 其中，f(m，n)是像素边缘检测的位置的灰度值；
[0014] 采用如下公式（18)求出H直：
[0016] 其中口为亚像素边缘点的角度，
[0017] 利用亚像素边缘点的角度供和如下公式（21)和（22)计算出亚像素边缘点离中心 的位置1的值：
[0020] 利用如下公式（27)得到图像的亚像素边缘位置：
[0022] 其中，X，y是Sobel算子进行检测得到的边缘点的位置，N代表掩码的窗口大小。
[0023] 由于采用了上述技术方案，本发明提供的基于Legendre矩的图像亚像素边缘检 测方法，首先对输入图像进行灰度化并采用自适应中值滤波器对图像进行去噪处理，然后 用Sobel算子进行像素级边缘粗定位，最后利用Legendre矩进行图像的亚像素级的边缘检 测，其中Sobel算子对噪声具有平滑作用，提供较为准确的边缘方向信息，利用Legendre矩 进行亚像素边缘检测，减少了运算所需要的模板的数量，降低了计算的复杂度，同时在抗噪 方面具有更好的鲁棒性。
【附图说明】
[0024] 为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下， 还可以根据这些附图获得其他的附图。
[0025] 图1为本发明的图像亚像素边缘检测方法的流程图；
[0026] 图2为本发明对图像进行去噪处理时边界扩展示意图；
[0027] 图3为本发明中理想2D边缘模型的示意图；
[0028] 图4(a)旋转前的阶跃边缘模型的示意图；
[0029] 图4(b)旋转后的阶跃边缘模型的示意图；
[0030] 图5为本发明中模板系数计算模型的示意图。
【具体实施方式】
[0031] 为使本发明的技术方案和优点更加清楚，下面结合本发明实施例中的附图，对本 发明实施例中的技术方案进行清楚完整的描述：
[0032] 如图1所示的基于Legendre矩的图像亚像素边缘检测方法，具体包括以下步骤： 读取图像信息，将图像灰度化并对灰度图像进行去噪处理：首先将输入RGB图像转换为灰 度图像。采用如下方式：
[0033] Gray= (28XB+151XG+77XR) ?8 (1)
[0034] 其中，">>"表示二进制右移。R、G、B代表红、绿、蓝三个通道的颜色。遍历输入图 像的所有像素点，对每个像素点采用式（1)，Gray为得到的灰度图像对应像素点的灰度值。
[0035] 对图像进行去噪声处理：
[0036] 本发明采用自适应中值滤波器，自适应中值滤波器可以处理含有大概率的脉冲噪 声图像，在平滑非脉冲噪声时能保留细节。自适应中值滤波器工作于矩形窗口区域Sxy内， 与传统的滤波器不同的是自适应中值滤波器在进行滤波处理时会根据某些条件而改变sxy 的尺寸。
[0037] 具体操作步骤：本发明中S_= 10,首先是边界扩展，在图像上下左右各增加S_ 个像素，如图2所示，设原图像的大小为mXn，区域1，2, 3, 4均是扩展的区域。区域1和区 域2的大小为S_Xn，区域3和区域4的大小为（m+2XS_)XS_，首先扩展区域1，将原图 像最左边与区域1对应大小的区域的数据复制到区域1，同理，扩展区域2,将原图像最右边 与区域2对应大小的区域的数据复制到区域2。接着，扩展区域3,将原图像和区域1和区 域2最上边（m+2XS_)XS_^t应的区域的数据复制到区域3,同理，扩展区域4.得到的扩 展图像进行后续的滤波处理。
[0038] 初始滤波半径r= 1，对应的初始矩形窗口Sxy的大小为（2r+l)X(2r+l)，（即为 3X3)，算法由两个进程表示，分别为进程A和进程B，依次从原图像的像素点开始遍历（不 包括填充的像素点）。
[0039] Zmin表不在矩形窗口Sxy中的最小灰度值，2_表不在矩形窗口Sxy中的最大灰度值， 表示在矩形窗口Sxy中的灰度中值，Zxy表示矩形窗口Sxy中心位置像素的灰度值， Sxy矩形窗口的最大滤波半径，在本文算法中，S_= 10。
[0040] 进程A:Ai=Zmed_Zmin
[0041] A2 =zmed-zmax
[0042] 如果APO且A2〈0,则转到进程B，否则滤波半径r=r+1
[0043] 如果矩形窗口的滤波半径〈=S_，则重复进程A
[0044] 否则输出Zmed
[0045]进程B:Bi=Zxy_Zmin [0046] B2=Zxy-Zfflax
[0047] 如果&>0且B2〈0,则输出Zxy [0048] 否则输出Zmed
[0049] 输出的Zmed即为该像素点经过自适应中值滤波器滤波后的像素值。
[0050] 采用Sobel算子对去噪后的图像进行像素级边缘定位：利用像素点上下，左右邻 点的灰度加权算法，根据在边缘点处达到极值现象进行图像边缘的检测。
[0051] Sobel算子很容易在空间上实现，Sobel边缘检测器不但产生较好的边缘检测效 果，而且受噪声影响也比较小。Sobel算子利用像素点上下，左右邻点的灰度加权算法，根据 在边缘点处达到极值这一现象进行边缘的检测。Sobel算子对噪声具有平滑作用，提供较为 准确的边缘方向信息。
[0052] fxr (x, y) =f(x-1,y+1) +2f(x, y+1)+f(x+1,y+1)
[0053] -f(x-1, y-1) -2f(x,y-1)-f(x+1, y-1) (2)
[0054] fy'(x, y) =f(x-1,y-1) +2f(x~l,y)+f(x~l,y+1)
[0055] -f(x+1, y-1) -2f(x+1,y)-f(x+1, y+1) (3)
[0056] G[fr (x,y)] = |fxr (x,y) | + |f/ (x,y) (4)
[0057] 其中，fx'（x，y)、fy'（x，y)分别是x(水平）方向和y(垂直）方向的一阶微分，G[f' (X，y)]为Sobel算子的梯度总和，f(x,y)是输入图像在坐标（X，y)点处的灰度值。
[0058] 在阈值T的设定方面，采用最大类间方差法（也称大津法，Otsu法），它的主要思 想是按照灰度特性将图像划分为背景和目标2部分，划分依据为选取门限值，使得背景和 目标之间的方差最大。其主要实现原理如下：
[0059] 1)建立图像灰度直方图（共有L个灰度级，每个出现概率为p，ni为灰度值为i的 像素的个数）
[0062] 2)计算背景和目标的出现概率，计算方法如下：
[0065] 其中，t为所选定的阈值，A代表背景（灰度级为0_t)，PA为背景出现的概率，同理B代表目标（灰度级为t+1-L-l)，PB为目标出现的概率。
[0066] 3)计算A和B两个区域的类间方差如下：
[0069] 〇2 =PA (oA-o〇)2+Pb (o