【2017年整理】一种高精度特征点检测和相机标定方法.docVIP

一种高精度特征点检测和相机标定方法 摘要：针对特征点检测和相机标定精度问题，提出了一种亚像素级的角点检测和相机标定方法。该方法对角点的提取首先运用传统的方法得到初始值，并以此为初始值从棋盘格图像中搜索边缘，然后拟合边缘数据，并将拟合得到的曲线的交点作为角点的最终值，最后以这些角点作为标定相机的特征点对相机进行标定。实验结果表明，该方法提高了特征点提取的准确度和精度，以此提高了相机标定的精度。?　　关键词：相机标定 亚像素 Tsai 拟合?　　0 引言?　　目前，CCD阵列作为一种普遍的成像器件，被应用于交通、军事、医疗、工业等领域。相机在成像过程中，成像的质量受到相机的分辨率、光学畸变的影响而下降，相机高精度的标定是测量系统完成任务必不可少的环节[1]，相机标定结果的精度制约着系统测量的精度。?　　不同的邻域由于受到不同因素的影响而采用的标定方法存在着较大的差异。机器视觉领域由于采用一般的CCD作为成像器件，不能准确地知晓其内部参数，而且具有随时间的推移而发生较大变化的特性，一次机器视觉领域相机地标定面临较多的困难。机器视觉领域的相机标定具有以下的特点：相机的空间分辨率一般比较低；使用的相机一般具有较大的畸变；相机部件存在较大的安装误差；标定受到光学系统色差、景深等因素的影响而使得标定模式的边缘产生一定的模糊。上述机器视觉领域中所面临的问题使得标定模式特征点的提取和定位精度较差，同时图像噪声的影响使得模式特征点的提取精度误差进一步加大。?　　标定模式特征点的定位精度很大程度上决定了相机标定的精度。即使是像素级的误差也会使得标定结果产生重大偏差，所以如何利用亚像元技术获得标定模式特征点的精确定位是实现相机高精度标定的关键[2]。?　　本文提出了一种基于全局曲线拟合和局部曲线拟合现结合的角点检测方法[3]。传统的焦点检测方法主要有双曲抛物面拟合和局部相交直线法，双曲抛物面拟合法假设所求角点附近图像灰度呈双曲抛物面分布，这样拟合曲面的鞍点[4]即为所求角点。局部相交直线法确定所求角点的方法是求棋盘格黑白交界处的交点[5]，由于受到镜头畸变的影响，交界处的直线在成像平面上不是一条直线，因此该方法所得角点的定位精度受到较大的影响[6，7]。上述两种方法在求角点时忽略了模式点之间的几何约束和相机标定参数的全局性特征。本文提出的方法通过对特征点进行几何约束来提高角点的定位精度。同时噪声对特征点提取的影响，曲线的拟合范围采用全局和局部相结合的方法来拟合曲线。?　　1 原理?　　1.1 相机模型?　　目前，存在的几种代表性的相机标定方法[8]，如Tsai[9]，Heikkil[3]和Zhang[4]等的方法，多数采用的是小孔成像模型。区分这些方法在于其所选用的内部参数的形式和表达方式。设OXYZ为世界坐标系，oxyz为相机坐标系，点P为三维空间中的一个点，坐标（xw，yw，zw）和（x，y，z）分别表示点P在世界坐标系和相机坐标系下的坐标。在理想的情况下，P点的图像坐标为Pu（Xu，Yu），但由于存在畸变的因素，其实际成像的位置为Pd（Xd，Yd），在计算机表示中用坐标（Xf，Yf）表示在图像中的像素坐标。?　　世界坐标系（xw，yw，zw）到相机坐标系（x，y，z）的变换如下?　　xyz=Rxwywzw+T （1）?　　其中，R和T分别为旋转矩阵和平移向量。相机坐标系到非畸变图像成像平面的投影变换为?　　其中f为焦距，畸变模型有下式给出?　　Xd+Dx=XuYd+Dy=Yu （3）?　　xm=A[Rt]M，A=α γ μ00 β ν00 0 1 （4）?　　其中m=[u，v，1]T是一个二维点的坐标矢量表示，M=[X，Y，Z，1]T是一个三维点坐标的矢量表示，s是图像的尺度模糊因子，（R，T）为世界坐标系到相机坐标系的旋转和平移的外参矩阵，A是相机的内参矩阵，（u0，v0）为主点坐标，α和β分别为图像坐标系u轴和v轴上的缩放因子，γ描述的是两个图像坐标轴之间的倾斜关系。?　　1.2 标定算法?　　1.2.1 特征点提取?　　相机标定的过程如图1所示，如前所述，进行相机标定需要提取出图像平面的特征点的位置，其提取的精度很大程度上决定了相机标定的精度。采用棋盘格模式对相机进行标定时，采用通常的方法提取出的角点与被检测的特征点邻域区域的图像信息有关，因此获取的特征点的定位精度受到了较大的影响，为了得到高精度的特征点定位估计，本文提出了相应的改进措施。?　　算法的主要步骤如下：?　　①在本文的标定算法中，使用如图所示的棋盘格模式，为降低图像中噪声对角点检测结果的影响，需要对图像进行平滑处理，在采用Harris角点检测算子得到图像中角点的初始值。?　　②检测标定板图像中的边缘信息，对步骤（Ⅰ）中得到的角点进行