一种基于分层处理的亚像素中心提取方法
【技术领域】
[0001] 本发明属于计算机视觉测量技术领域，涉及一种基于分层处理的亚像素中心提取 方法。
【背景技术】
[0002] 大型零构件制造和装配精度是飞机、汽车等重大装备质量和性能的重要保障，因 此高精度三维测量在制造装配过程中至关重要。基于辅助激光的视觉测量具有非接触性、 大信息、测量速度快等优势，已广泛运用在大型零件测量中。是通过双目相机实时采集由 激光器投射在被测物表面的激光扫描信息，基于双目立体视觉理论实现激光特征信息的重 建，进而实现被测物表面形面测量。其中光条中心的提取是保证测量精度的重要环节，但是 目前像素级的提取已经不能满足测量精度的要求；而高精度的光条的中心提取需要耗费较 多的时间，不能满足快速测量的要求。如何实现光条中心的快速高精度的亚像素提取是目 前需重要解决的问题。
[0003]C.Steger发表的期刊《Anunbiaseddetectorofcurvilinearstructures》 (C.Steger,Anunbiaseddetectorofcurvilinearstructures[J],PatternAnalysis andMachineIntelligence,IEEETransactionson，1998,20(2) ,113-125·)提出了 Hessian矩阵得到图像中条纹的法线方向，通过法线方向求极值点得到光条纹中心的亚像 素位置，该刚发具有精度高、稳健性好等优点，但运算量大，难以实现快速提取光条中心。魏 振忠等人发明的"一种结构光光条中心的提取方法"，专利号CN101504770。该发明通过光 条中心点与所述光条中心点处法线方向之间的多次迭代，获取结构光光条中心；对所述获 取的结构光光条中心进行平滑操作，提取最终的结构光光条中心，该方法测量精度高可以 去除冗余的计算步骤，但是计算时间较长无法满足现场快速测量要求。
【发明内容】
[0004] 本发明要解决的技术难题是针对难以实现光条中心的快速高精度的亚像素提取 问题，发明了一种基于分层处理的亚像素中心提取方法。该方法采用分层处理，通过在低分 辨率图像中计算光条中心法向量，在高分辨率图像高精度提取光条中的亚像素，在光条法 向量上采用灰度重心法计算光条中心，实现光条的亚像素中心的提取；提高提取效率，实现 快速高精度的光条中心提取。
[0005] 本发明采用的技术方案是一种基于分层处理的亚像素中心提取方法，其特征是， 该方法采用分层处理，通过在不同分辨率图像中提取光条的特征信息，实现特征光条的高 精度的亚像素中心提取；该方法首先通过对原图像的压缩，将高分辨率图像压缩为低分辨 率图像，针对低分辨率图像提取光条的法向量；然后在高分辨率图像中，采用差分的方式计 算高分辨率图像中光条的法向量；最后在光条法向量上采用灰度重心法计算光条中心，实 现光条的亚像素中心的提取，并通过标准光条的重建验证了该提取方法的有效性；方法具 体步骤如下：第一步低分辨图像中的光条法向量提取
[0006] 1)低分辨率图像压缩方法
[0007] 对于双目相机采集的特征光条图像，其分辨率为mXn，按照压缩比例因子τ定义 τ为4的整数倍，S卩τ=4s(s= 1，2,3...)，将图像进行压缩，即将图像中τXτ个像素 进行灰度的累加，然后求平均灰度，压缩为一个像素；其计算准则如下：
[0008] \<j~]i J j
y .. fc
[0009]其中，为压缩图像的第i行、j列的像素灰度值，ItUi) +k]1) +。为在高分辨 图像中第τ(i-D+δ行、τ(j-l) + 〇列的像素灰度值，〇和δ分别是计算系数，π^Ρη分 别是高分辨率图像的行列数；
[0010] 2)低分辨率图像中的光条中心提取
[0011] 将低分辨率图像进行滤波处理，然后采用Sobel边缘检测算子提取光条边界，计 算方法如下：
[0012]
[0013]
[0014] 其中，g(i，j)为图像坐标为（i，j)的Sobel边缘检测算子，4和dy分别为计算的 卷积模板。
[0015] 根据所提取的光条边界图像，搜寻低分辨率图像中光条在每行中的光条边界，采 用几何中心法计算光条的中心坐标：
[0016]
[0017] 其中，（v;〇 分别为图像第i行左、右边界像素坐标；当像素坐标为非 整数时，对像素坐标进行求整，其几何中心的像素坐标为f= (T.ih
[0018] 3)低分辨率图像中光条中心法矢量的计算
[0019] 选取光条几何中心的前后w行的中心进行曲线拟合，每次拟合的点数至少为四个 点，其选取的光条几何中心数列为：
[0020] {Pk |k=i-w,i-w+1,...,i-1,i,i+1,...,i+w-1,i+w；k1；w^ 3} (4)
[0021] 其中i为低分辨率图像的行数。因为局部曲线较小，采用二次曲线进行拟合，其拟 合公式为：
[0022] 仁（X) =aj^biX+Ci(5)
[0023] 其中，&1、131、(：1分别为第1行中心拟合曲线的拟合系数，根据拟合曲线通过偏导数 的求解，计算第i行中心的法向量名；
[0024]
[0025] 其中，4#-1为^ (x)在x= \位置时的相对于x的偏导数，即为该点的法向量， 根据法向量求解公式计算出在低分辨率图像中的每个几何中心的法相方向；
[0026] 第二步高分辨率图像中的光条中心亚像素提取
[0027] 1)光条中心粗定位
[0028] 对高分辨图像进行滤波处理，按照低分辨率图像中边界提取算法，即公式（1)提 取高分辨率中的光条边界，根据第二步中几何中心的计算方法，在光条宽度上计算光条第i 行的几何中心，其计算公式如下：
[0029]
(7)
[0030] 其中，)、?,〇分别为图像第i行左、右边界像素坐标；当像素坐标为非 整数时，对像素坐标进行求整，其几何中心的像素坐标为if= (?」；
[0031] 2)高分辨图像中的法矢量计算
[0032] 将低分辨率图像求解的光条中心的法向量还原到高分辨图像的光条中，低分辨率 图像的一行像素代表高分辨率图像的τ行像素，定义低分辨率中每行中心的光条法向量 为高分辨率图像中τ/2行光条的法向量，即
[0033] qt ~(Jr(i-i)+ir+i)/2 雜）
[0034] 其中qMllWT+1)/2为在高分辨图像中τα-1)+ (τ+1)/2行位置的法向量，根据前 后两个光条法向量，通过差值计算高分辨率图像中中间行的法向量，计算公式如下：
[0035]
(9):
[0036] 3)灰度重心法精提取光条中心
[0037] 沿每行光条中心的法向量搜索光条，在光条法向宽度上计算光条的灰度重心，计 算方法如Τ.
[0038] (1〇) β-?：
[0039] 其中，（ue，ve)为第e个光条灰度重心坐标，L为第i行第j列灰度值；
[0040] 第三步光条提取精度的验证
[0041] 通过重建光条信息，计算标准圆柱件的重建精度验证光条的提取精度；首先基于 极限约束对左右高分辨图像提取的灰度重心进行匹配，计算方法如下：
[0042] X:/ ... = 0 (11)
[0043] 其中，x/ = (I，vj为左相机的灰度重心坐标；x^' = (u/，v/)为与x/相匹配 右图像重心坐标；F为通过八点法计算的两相机之间的基本矩阵；若左右图像中两点满公 式（11)，则左右图像中的重心点为匹配点；将左右图像匹配上的中间点进行三维重建，通 过三维重建点拟合标准圆柱件直径Rf，计算拟合直径Rf与标准直径R径评价重建精度η， 计算公式如下：
[0044]
(12)
[0045] 从而验证了该提取方法的有效性。
[0046] 本发明的有益效果是提取方法首先通过对原图像的压缩，针对低分辨率图像提取 光条的法向量；然后在高分辨率图像中，采用差分的方式计算高分辨率图像中光条的法向 量；最后在光条法向量上采用灰度重心法计算光条中心，实现光条的亚像素中心的提取。该 方法测量精度高，运算速度快，计算时间短，实现了光条的亚像素提取，可满足光条中心的 高精度快速提取。
【附图说明】
[0047]图1为高分辨率光条图像的示意图。其中，1-高分辨图像，2-原始光条。
[0048] 图2为压缩后的低分辨率图像，其中，3-低分辨率图像，4-压缩后光条，<、€、< 分别是压缩后的低分辨率图像中第1、2、3行的法向量。
[0049] 图3为基于分层处理的亚像素中心提取流程图。
【具体实施方式】
[0050] 以下结合技术方案和附图详细叙述本发明的【具体实施方式】。
[0051] 实施例中，本发明采用分别配置广角镜头的两个摄像机拍摄一幅光条图像。摄 像机型号为viewworksVC-12MC-M/C65摄像机，分辨率：4096X3072,图像传感器：CMOS， 帧率：全画幅，最高64. 3fps，重量：420g。广角镜头型号为E