背光图像微米级边缘检测方法
【技术领域】
[0001] 本发明属于高精度检测技术领域，尤其涉及一种背光图像微米级边缘检测方法。
【背景技术】
[0002] 机器视觉测量系统集光学、传感器、图像处理和模式识别等多领域多学科的关键 技术于一体，在非接触高速在线测量领域应用广泛。本发明利用机器视觉和数字图像处理 技术，针对背光数字图像，建立法向正态概率分布曲线图像边缘灰度曲面模型，用最小二乘 法求取拟合曲面四个系数，得到微米级亚像素边缘。
[0003] 专利【申请号】2011104481191，名称为"一种彩色图像边缘检测方法"，将彩色图像 分解到红色、绿色、蓝色和黄色四个通道图像，分别计算红绿、蓝黄拮抗图像，并在拮抗图像 上计算得到边缘信息分布图像，通过调节计算生成颜色拮抗图像过程中引入的两组拮抗权 重因子，调节图像中颜色边缘和亮度边缘的强度，实现边缘检测。
[0004] 专利【申请号】2011100826554,名称为"一种图像边缘检测的方法"，在利用摄像机 获得一幅图像后，首先用一种在12邻域像素内计算X方向和y方向一阶偏导数差分的方法 计算梯度幅值和方向，抑制图像中的噪声并且避免图像的边缘模糊。再建立灰度共生矩阵 逆差矩特征值和高斯空间系数及阈值之间的映射关系，自适应地改变高斯空间系数和边缘 检测的高低阈值，确保了图像边缘点的连续提取。最后根据"非局部最大值抑制"原理检测 图像边缘点，提高了图像边缘检测的精确度。
[0005] 专利【申请号】2008102323401，名称为"基于Fast Slant Stack变换的图像边缘检 测方法"。其过程为：根据滑动窗大小和重叠度对图像分块，每次检测提取一个小图像块，根 据其灰度均方差判断块中是否存在边缘，对存在边缘的图像块做Fast Slant Stack变换， 并对变换结果的零值区域进行填充；对填充后的图像块做二进小波变换，搜索变换后小波 域系数极大值，并根据极大值所处的位置用基于共辄梯度法的Fast Slant Stack逆变换重 构空域中的图像，存放到输出矩阵中对应的位置。
[0006] 专利【申请号】2010102516986,名称为"一种基于数学形态学的图像边缘检测方 法"，采用不同尺度的结构元素对图像进行形态学处理，获取边缘信息图像，然后对边缘信 息图像进行加权合并，获取总边缘信息图像，最后对总边缘信息图像的边缘点进行细化和 阈值分割，得到最终的二值边缘图像。
[0007] 专利【申请号】2008101877316,名称为"一种基于Bayer RGB的图像边缘检测方 法"，以目标像素点为中心，生成3x3邻域的Bayer RGB像素点矩阵，分别对该矩阵中各像素 点，应用检测算子进行梯度运算，得到各像素点的梯度幅值以及梯度角度。然后采用双阈值 判定和相邻像素点梯度比较，判定是否为边缘点。免除了像素点的灰度计算，适用于对CMOS 图像传感器所得的原始图像直接进行边缘检测。
[0008] 专利【申请号】2007101787372,名称为"一种基于阈值分割的图像边缘检测方法"。 利用图像中像素邻域的灰度值均值统计信息作为该点阈值设置的标准，并引入该点邻域内 像素灰度值方差作为附加判断条件，从而实现综合考虑区域信息的全局二值化，使提取出 来的目标点是图像的边缘。根据不同的需要来选择不同大小的结构元素。对每一个像素都 给定一个阈值，从而实现综合考虑区域信息的全局二值化。
[0009] 专利【申请号】2014104698402,名称为"一种基于机器视觉的快速亚像素边缘检测 与定位方法"，包括以下步骤：获取检测图像；图像去噪预处理：计算每个像素点水平方向 上的梯度G x和垂直方向上的梯度G y;计算每个像素点极坐标下的梯度幅值G。和梯度方向 Ge;确定每个像素点的邻域像素点；确定像素级边缘点；计算每个像素级边缘点在8分梯度 方向上的亚像素边缘点到该像素级边缘点的距离；计算实际梯度方向G e上的亚像素边缘 点到像素级边缘点的距离d;采用余弦查找表法计算其实际梯度方向Ge上的亚像素边缘点 的直角坐标，进而实现图像边缘点的检测和亚像素级定位。
[0010] 以上发明基于机器视觉和图像处理技术，利用颜色、数学形态学、阀值法、梯度法 等方法，检测数字图像中的像素级和亚像素级边缘。这些方法适用于正面照明图像，都不适 用于背光图像的微米级边缘检测。
【发明内容】
[0011] 本发明就是针对上述问题，提供一种曲面拟合精度高、抗干扰的背光图像微米级 边缘检测方法。
[0012] 为实现上述目的，本发明采用如下技术方案，本发明包括以下步骤：
[0013] 1)结合点扩散函数和单边阶跃模型，构建背光图像阶跃边缘法平面截线正态分布 函数模型；
[0014] 背光源边缘检测图像的前景和背景为常量Pp P2，背景与前景的灰度差为
[0015] g = P2-P1
[0016] 灰度曲面过边缘法平面的截线为正态分布函数（高斯积分）曲线：
[0020] 图像边缘点过曲线的中心点在图像平面的投影，@ ? U丨as: I或t_ μ = 0。
[0021] 2)过渡带离散点曲面拟合求解
[0022] μ，〇沿边缘线S方向在邻域内呈线性变化
[0023] μ s= μ 〇+ks (I)
[0024] 〇s=o 0+ls (2) I , t ? p .
[0025] 由 於众》二參（-
[0026] 查标准正态分布表得到 「 ? L0027J -=艺
[0028] 将（1)、⑵式代入，得
[0029] t- μ 〇-ks = a ( σ 〇+ls) 「00301 田?/丨、一乘法步取相合曲而PU个赛翁u "、σ "、k、I .伸
[0036] 其中，〇 (即上方带横杠的标识）表示平均值；
[0037] 3)分段提取边缘曲面亚像素边缘
[0038] 在图像平面内建立边缘曲线活动坐标系；以像素级别的粗边界曲线方向为活动坐 标系的一个坐标，以曲线弧长作为其坐标值；以粗边界曲线的法线方向为活动坐标系的另 一个坐标，以点到曲线的法向距离作为其坐标值；以计算点为中心选取曲面拟合区域；在 边缘法线方向上，区域包含边缘过渡区和前景与背景图像部分；将选取区域内点的图像坐 标转换为边缘曲线的活动坐标，得到tp 8;及相应的灰度值P i;根据P i/g查标准正态分布 表，得到a1;求得μ。，即边缘到粗边界曲线的法向距离，得到精确的边缘曲线。
[0039] 作为一种优选方案，本发明还包括步骤4)背光源检测图像边缘时，图像边缘清 晰，边缘过渡带窄，将μ，σ在邻域内作为常量，即
[0040] μ,= μ 〇 (4)
[0041] Os=O0 (5)
[0042] 简化计算
[0043] O
[0044] 作为另一种优选方案，本发明所述步骤3)选取区域为背光图像边缘过渡带宽度 的10至20倍。
[0045] 其次，本发明所述步骤4)边缘过渡带为3-4个像素。
[0046] 本发明有益效果。
[0047] (1)针对背光图像，构建了阶跃边缘法平面截线正态分布函数模型，并给出了微米 级亚像素边缘定位方法的完整步骤。本方法属于曲面拟合，具有曲面拟合精度高，抗干扰的 优点。
[0048] (2)曲面拟合时，采用的不是行列方向的邻域点，而是切线和法线方向的过渡带数 据点，提高了精度。曲面切向和法向分离，减少了计算复杂度。单纯的曲面拟合算法，仅仅 减少了噪声的影响，本方法将离散点纠正到理论曲面上，极大地消除了噪声影响。
[0049] (3)对量块边缘的检测表明，直线度误差在± 1. 5 μ m以内。因此，