一种基于Legendre矩的图像亚像素边缘检测方法
一种基于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 B-o〇)2 (11)
[0070] 式（9)分别计算A和B区域的平均灰度值，代表A区域的平均灰度值，《 表B区域的平均灰度值；式（10)计算灰度图像全局的灰度平均值《。;式（11)计算A、B两 个区域的类间方差〇2。
[0071] 4)以上几个步骤计算出了单个灰度值上的类间方差，因此最佳分割门限值应该是 图像中能够使得A与B的类间灰度方差最大的灰度值。在程序中t的取值从0到255,依次 计算式（11)的值，最大的〇2所对应的t值即为阈值T。在得到像素级粗定位后，遍历图像 中所有的边缘点，并进行判断，若该边缘点是孤立的边缘点（即在以该点为中心的3X3的 邻域中（除本点外）边缘点的个数小于等于1，），则将该点除去，即该点不作为边缘点，判 断为噪声。
[0072] 采用Legendre矩对图像进行亚像素边缘检测，输出边缘图像。
[0073] 自从TABATABAI等在1984年提出了利用灰度矩进行亚像素级别的边缘检测，经 过20多年的研宄，其他方法例如空间矩，Zernike矩，0FMM等等被提出，这些方法假设理想 的边缘是阶跃模型，通过将图像映射到单位圆内，求得亚像素的4个参数，1 (亚像素离中 心的位置），P(亚像素的角度），k(灰度的阶跃值），h(背景灰度）。为此，提出了基于 Legendre矩的亚像素边缘检测。
[0074] (1)Legendre矩
[0077] 在单位圆内，Legendre矩可以定义为：
[0079] 其中
-是归一化系数，f(r，0 )为原灰度图像在（x，y)点处 的灰度值
[0080]f (X，y)在极坐标下的表示形式.其中
[0082] 核函数!》=Qn(r)exp(_jm0 )使具有旋转不变性。
[0083]
[0084] 其中，111"表示原始图像的Legendre矩，LM' "表示对图像旋转供角后的Legendre 矩。如图3所示。
[0085] 基于Legendre矩的边缘检测：如图4所示，在旋转供的角度后，边缘与y轴垂直， 旋转后的图像函数的积分具有如下的关系：
[0089] 根据式（9)，可以得到LM/i的虚部为0,所以
[0091]Re[LMn]和Im[LMn]分别是实部和虚部。
[0093] 积分核函数可以表示为：
[0096]注：LM/ ^LM3' ：的解法可见[2]
[0099] 计算LMjPLM31模板的系数，本发明使用5X5的模板
[0102] 如图5和式（23)，（24)可知，计算的系数实部关于7轴奇对称，关于x轴偶对称， 虚部关于x轴奇对称，关于y轴偶对称。因此，仅需要计算图5中，方格1，2, 3, 6, 7,8, 11，12 这八个系数，剩下的系数可通过对称性得到。
[0103] 首先利用式（23)计算CLMn
[0113] 根据对称性，可得剩下的系数，详细见表1.
[0114] 同理，利用式（24)可计算CLM31模板的系数，见表2
[0115] 表lCLMn模板系数
[0117] 表2 01131模板系数
[0118]
[0119] 然后利用以下的公式，求LMjPLM31
[0122] 其中，f(m，n)是像素边缘检测的位置的灰度值。
[0123] 真实的边缘位置为：
[0125] 其中，x，y是Sobel算子进行检测得到的边缘点的位置，N代表掩码的窗口大小。
[0126] 本发明的有益效果说明：
[0127] 为了验证本发明，进行了计算机仿真实验。在实验中，实验参数为CPU英特尔 Pentium(奔腾）双核E53002.6GHz，2GB内存，显卡是英特尔G33/G31ExpressChipset Family,操作系统为WindowXP专业版32位SP2,软件编程环境为Matlab2010b，本发明实 验的图像是利用人工合成的图像，对于人工合成的图片的大小是256像素X256像素。
[0128] 从[3]中，FeipengDa推导了SGM，Z0M和0FMM之间的关系，可以得出计算SGM，Z0M 和0FMM所得的0是一样的，而Z0M和0FMM的1值是一样的，由SGM和Z0M所计算得到的1 值的差值是
[0130] 所以将发明方法与SGM和Z0M进行了仿真对比实验。测试的图片为添加高斯 白噪声的不同半径的圆，圆心为（128, 128)，通过将计算得到的亚像素点进行拟合，求出 (X-A)2+(Y-B)2=R2中A，B，R的值，拟合采用[6]中所提到的方法。
[0133] k是亚像素边缘点的个数，xt，yt代表了第i个亚像素边缘点的坐标，半径定义为 亚像素边缘点到实际圆心的平均距离.
[0134] 表3不同方法的位置边缘误差（拟合圆心与实际圆心的误差，注，误差计算的是 欧式距离）
[0136] 表4不同方法的位置边缘误差（拟合半径与实际半径的误差，注，误差计算的是 两者的差值）
[0138] 参考文献：
[0139] (如专利/论文/标准）
[0140] [l]BinTJ,LeiA,JiwenC,etal.Subpixeledgelocation basedonorthogonalFourier-Mel1inmoments[J].ImageandVision Computing, 2008, 26(4) : 563-569.
[0141] [2]CuiJ，FengK，TanJB.Furtherimprovementofedgelocation accuracyofdoublefibersphericalcouplingsensorusingorthogonalJacobi-Fouriermoments[J] ?Optik-InternationalJournalforLightandElectron Optics，2014, 125 (1) : 353-359.
[0142] [3]DaF，ZhangH.Sub-pixeledgedetectionbasedonanimproved moment[J].ImageandVisionComputing, 2010, 28(12):1645-1658.
[0143] [4]LyversEP,Mitchell0R，AkeyML，etal.Subpixelmeasurementsusing amoment-basededgeoperator[J].PatternAnalysisandMachineIntelligence,IEEE Transactionson，1989, 11 (12):1293-1309.
[0144] [5]GhosalS，MehrotraR.Orthogonalmomentoperatorsforsubpixeledge detection[J].Patternrecognition, 1993, 26(2):295-306.
[0145] [6]FabijanskaA.Gaussian-basedapproachtosubpixeldetection ofblurredandunsharpedges[C]//ComputerScienceandInformation Systems(FedCSIS), 2014FederatedConferenceon.IEEE,2014:641-650.
[0146] 以上所述，仅为本发明较佳的【具体实施方式】，但本发明的保护范围并不局限于此， 任何熟悉本技术领域的技术人员在本发明掲露的技术范围内，根据本发明的技术方案及其 发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
【主权项】
1. 一种基于Legendre矩的图像亚像素边缘检测方法，其特征在于包括以下步骤： SI:读取图像信息，将图像灰度化并对灰度图像进行去噪处理； S2:采用Sobel算子对去噪后的图像进行像素级边缘定位：利用像素点的各向邻点灰 度加权值在边缘点达到最大值这一现象进行边缘检测； S3:采用Legendre矩对图像进行亚像素边缘检测，输出边缘图像。2. 根据权利要求1所述的一种基于Legendre矩的图像亚像素边缘检测方法，其特征还 在于：S2中具体采用如下方式：遍历原始灰度图像中所有的像素点，计算得到每个像素点 的梯度值G[f' (x，y)]，将所得的梯度值归一化到[0, 255]区间，采用最大类间方差法计算 得到归一化梯度值的阈值T，对每个像素点归一化的梯度值进行判断，即当G[f' (x，y)]>T 时，对应的像素点设定为255,否则设定为0至此得到图像的像素级粗定位。3. 根据权利要求2所述的一种基于Legendre矩的图像亚像素边缘检测方法，其特征还 在于：在得到像素级粗定位后遍历图像中所有的边缘点，进行判断：若该边缘点是孤立的 边缘点即在以该点为中心的3X3的邻域中、除本点外的边缘点的个数小于等于1，则将该 点除去，即该点不作为边缘点、判断为噪声。4. 根据权利要求1所述的一种基于Legendre矩的图像亚像素边缘检测方法，其特征 还在于：S3中具体采用如下方式：遍历检测到的所有边缘点，对每个边缘点进行如下处理： 以得到的边缘点为中心，在灰度图像中选取NXN的窗口，N为奇数，采用如下公式（25) 将NXN灰度图像窗口中的值与Legendre正交矩的掩码CLM11对应位置的系数相乘得到 NXN的矩阵，将该矩阵求和得到Legendre正交矩LMn，同样的方式利用式（26)再求得一个 Legendre 正交矩 LM31，其中，f(m，η)是像素边缘检测的位置的灰度值；采用如下公式（18)求出0值： (25) (26)(18) 其中供为亚像素边缘点的角度， 利用亚像素边缘点的角度P和如下公式（21)和（22)计算出亚像素边缘点离中心的位 置1的值：利用如下公式（27)得到图像的亚像素边缘位置：(27) 其中，χ，y是Sobel算子进行检测得到的边缘点的位置，N代表掩码的窗口大小。
【专利摘要】本发明公开了一种基于Legendre矩的图像亚像素边缘检测方法，包括以下步骤：S1:读取图像信息，将图像灰度化并对灰度图像进行去噪处理；S2:采用Sobel算子对去噪后的图像进行像素级边缘定位：利用像素点的各向邻点灰度加权值在边缘点达到最大值这一现象进行边缘检测；S3:采用Legendre矩对图像进行亚像素边缘检测，输出边缘图像。其中Sobel算子对噪声具有平滑作用，提供较为准确的边缘方向信息，利用Legendre矩进行亚像素边缘检测，减少了运算所需要的模板的数量，降低了计算的复杂度，同时在抗噪方面具有更好的鲁棒性。
【IPC分类】G06T7/00
【公开号】CN104899888
【申请号】CN201510340586
【发明人】陈喆, 殷福亮, 张一
【申请人】大连理工大学
【公开日】2015年9月9日
【申请日】2015年6月18日
【技术领域】
[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 B-o〇)2 (11)
[0070] 式（9)分别计算A和B区域的平均灰度值，代表A区域的平均灰度值，《 表B区域的平均灰度值；式（10)计算灰度图像全局的灰度平均值《。;式（11)计算A、B两 个区域的类间方差〇2。
[0071] 4)以上几个步骤计算出了单个灰度值上的类间方差，因此最佳分割门限值应该是 图像中能够使得A与B的类间灰度方差最大的灰度值。在程序中t的取值从0到255,依次 计算式（11)的值，最大的〇2所对应的t值即为阈值T。在得到像素级粗定位后，遍历图像 中所有的边缘点，并进行判断，若该边缘点是孤立的边缘点（即在以该点为中心的3X3的 邻域中（除本点外）边缘点的个数小于等于1，），则将该点除去，即该点不作为边缘点，判 断为噪声。
[0072] 采用Legendre矩对图像进行亚像素边缘检测，输出边缘图像。
[0073] 自从TABATABAI等在1984年提出了利用灰度矩进行亚像素级别的边缘检测，经 过20多年的研宄，其他方法例如空间矩，Zernike矩，0FMM等等被提出，这些方法假设理想 的边缘是阶跃模型，通过将图像映射到单位圆内，求得亚像素的4个参数，1 (亚像素离中 心的位置），P(亚像素的角度），k(灰度的阶跃值），h(背景灰度）。为此，提出了基于 Legendre矩的亚像素边缘检测。
[0074] (1)Legendre矩
[0077] 在单位圆内，Legendre矩可以定义为：
[0079] 其中
-是归一化系数，f(r，0 )为原灰度图像在（x，y)点处 的灰度值
[0080]f (X，y)在极坐标下的表示形式.其中
[0082] 核函数!》=Qn(r)exp(_jm0 )使具有旋转不变性。
[0083]
[0084] 其中，111"表示原始图像的Legendre矩，LM' "表示对图像旋转供角后的Legendre 矩。如图3所示。
[0085] 基于Legendre矩的边缘检测：如图4所示，在旋转供的角度后，边缘与y轴垂直， 旋转后的图像函数的积分具有如下的关系：
[0089] 根据式（9)，可以得到LM/i的虚部为0,所以
[0091]Re[LMn]和Im[LMn]分别是实部和虚部。
[0093] 积分核函数可以表示为：
[0096]注：LM/ ^LM3' ：的解法可见[2]
[0099] 计算LMjPLM31模板的系数，本发明使用5X5的模板
[0102] 如图5和式（23)，（24)可知，计算的系数实部关于7轴奇对称，关于x轴偶对称， 虚部关于x轴奇对称，关于y轴偶对称。因此，仅需要计算图5中，方格1，2, 3, 6, 7,8, 11，12 这八个系数，剩下的系数可通过对称性得到。
[0103] 首先利用式（23)计算CLMn
[0113] 根据对称性，可得剩下的系数，详细见表1.
[0114] 同理，利用式（24)可计算CLM31模板的系数，见表2
[0115] 表lCLMn模板系数
[0117] 表2 01131模板系数
[0118]
[0119] 然后利用以下的公式，求LMjPLM31
[0122] 其中，f(m，n)是像素边缘检测的位置的灰度值。
[0123] 真实的边缘位置为：
[0125] 其中，x，y是Sobel算子进行检测得到的边缘点的位置，N代表掩码的窗口大小。
[0126] 本发明的有益效果说明：
[0127] 为了验证本发明，进行了计算机仿真实验。在实验中，实验参数为CPU英特尔 Pentium(奔腾）双核E53002.6GHz，2GB内存，显卡是英特尔G33/G31ExpressChipset Family,操作系统为WindowXP专业版32位SP2,软件编程环境为Matlab2010b，本发明实 验的图像是利用人工合成的图像，对于人工合成的图片的大小是256像素X256像素。
[0128] 从[3]中，FeipengDa推导了SGM，Z0M和0FMM之间的关系，可以得出计算SGM，Z0M 和0FMM所得的0是一样的，而Z0M和0FMM的1值是一样的，由SGM和Z0M所计算得到的1 值的差值是
[0130] 所以将发明方法与SGM和Z0M进行了仿真对比实验。测试的图片为添加高斯 白噪声的不同半径的圆，圆心为（128, 128)，通过将计算得到的亚像素点进行拟合，求出 (X-A)2+(Y-B)2=R2中A，B，R的值，拟合采用[6]中所提到的方法。
[0133] k是亚像素边缘点的个数，xt，yt代表了第i个亚像素边缘点的坐标，半径定义为 亚像素边缘点到实际圆心的平均距离.
[0134] 表3不同方法的位置边缘误差（拟合圆心与实际圆心的误差，注，误差计算的是 欧式距离）
[0136] 表4不同方法的位置边缘误差（拟合半径与实际半径的误差，注，误差计算的是 两者的差值）
[0138] 参考文献：
[0139] (如专利/论文/标准）
[0140] [l]BinTJ,LeiA,JiwenC,etal.Subpixeledgelocation basedonorthogonalFourier-Mel1inmoments[J].ImageandVision Computing, 2008, 26(4) : 563-569.
[0141] [2]CuiJ，FengK，TanJB.Furtherimprovementofedgelocation accuracyofdoublefibersphericalcouplingsensorusingorthogonalJacobi-Fouriermoments[J] ?Optik-InternationalJournalforLightandElectron Optics，2014, 125 (1) : 353-359.
[0142] [3]DaF，ZhangH.Sub-pixeledgedetectionbasedonanimproved moment[J].ImageandVisionComputing, 2010, 28(12):1645-1658.
[0143] [4]LyversEP,Mitchell0R，AkeyML，etal.Subpixelmeasurementsusing amoment-basededgeoperator[J].PatternAnalysisandMachineIntelligence,IEEE Transactionson，1989, 11 (12):1293-1309.
[0144] [5]GhosalS，MehrotraR.Orthogonalmomentoperatorsforsubpixeledge detection[J].Patternrecognition, 1993, 26(2):295-306.
[0145] [6]FabijanskaA.Gaussian-basedapproachtosubpixeldetection ofblurredandunsharpedges[C]//ComputerScienceandInformation Systems(FedCSIS), 2014FederatedConferenceon.IEEE,2014:641-650.
[0146] 以上所述，仅为本发明较佳的【具体实施方式】，但本发明的保护范围并不局限于此， 任何熟悉本技术领域的技术人员在本发明掲露的技术范围内，根据本发明的技术方案及其 发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
【主权项】
1. 一种基于Legendre矩的图像亚像素边缘检测方法，其特征在于包括以下步骤： SI:读取图像信息，将图像灰度化并对灰度图像进行去噪处理； S2:采用Sobel算子对去噪后的图像进行像素级边缘定位：利用像素点的各向邻点灰 度加权值在边缘点达到最大值这一现象进行边缘检测； S3:采用Legendre矩对图像进行亚像素边缘检测，输出边缘图像。2. 根据权利要求1所述的一种基于Legendre矩的图像亚像素边缘检测方法，其特征还 在于：S2中具体采用如下方式：遍历原始灰度图像中所有的像素点，计算得到每个像素点 的梯度值G[f' (x，y)]，将所得的梯度值归一化到[0, 255]区间，采用最大类间方差法计算 得到归一化梯度值的阈值T，对每个像素点归一化的梯度值进行判断，即当G[f' (x，y)]>T 时，对应的像素点设定为255,否则设定为0至此得到图像的像素级粗定位。3. 根据权利要求2所述的一种基于Legendre矩的图像亚像素边缘检测方法，其特征还 在于：在得到像素级粗定位后遍历图像中所有的边缘点，进行判断：若该边缘点是孤立的 边缘点即在以该点为中心的3X3的邻域中、除本点外的边缘点的个数小于等于1，则将该 点除去，即该点不作为边缘点、判断为噪声。4. 根据权利要求1所述的一种基于Legendre矩的图像亚像素边缘检测方法，其特征 还在于：S3中具体采用如下方式：遍历检测到的所有边缘点，对每个边缘点进行如下处理： 以得到的边缘点为中心，在灰度图像中选取NXN的窗口，N为奇数，采用如下公式（25) 将NXN灰度图像窗口中的值与Legendre正交矩的掩码CLM11对应位置的系数相乘得到 NXN的矩阵，将该矩阵求和得到Legendre正交矩LMn，同样的方式利用式（26)再求得一个 Legendre 正交矩 LM31，其中，f(m，η)是像素边缘检测的位置的灰度值；采用如下公式（18)求出0值： (25) (26)(18) 其中供为亚像素边缘点的角度， 利用亚像素边缘点的角度P和如下公式（21)和（22)计算出亚像素边缘点离中心的位 置1的值：利用如下公式（27)得到图像的亚像素边缘位置：(27) 其中，χ，y是Sobel算子进行检测得到的边缘点的位置，N代表掩码的窗口大小。
【专利摘要】本发明公开了一种基于Legendre矩的图像亚像素边缘检测方法，包括以下步骤：S1:读取图像信息，将图像灰度化并对灰度图像进行去噪处理；S2:采用Sobel算子对去噪后的图像进行像素级边缘定位：利用像素点的各向邻点灰度加权值在边缘点达到最大值这一现象进行边缘检测；S3:采用Legendre矩对图像进行亚像素边缘检测，输出边缘图像。其中Sobel算子对噪声具有平滑作用，提供较为准确的边缘方向信息，利用Legendre矩进行亚像素边缘检测，减少了运算所需要的模板的数量，降低了计算的复杂度，同时在抗噪方面具有更好的鲁棒性。
【IPC分类】G06T7/00
【公开号】CN104899888
【申请号】CN201510340586
【发明人】陈喆, 殷福亮, 张一
【申请人】大连理工大学
【公开日】2015年9月9日
【申请日】2015年6月18日
最新回复(0)