发明名称 一种用于摄像机标定的平面靶标特征点的自动匹配方法
摘要 本发明属于机器视觉技术领域,将提供一种用于摄像机标定的平面靶标特征点的自动匹配方法。本发明通过对提取的标定特征点集进行计算几何理论中Delaunay三角化和必要的预处理,得到相邻特征点之间的几何连接关系;然后搜索点集中用于确定搜索方向的起始点、第一行的特征点和第一列的特征点;最后利用特征点之间的几何连接关系完成其余各行的匹配,实现呈方格阵列分布的平面靶标特征点的图像坐标与世界坐标的匹配。本发明提出的特征点的匹配方法实现简单,鲁棒性强,对镜头畸变不敏感,适合在拍摄倾角大、靶标被部分遮挡情况下的特征点的匹配,特别适合摄像机的在线自动标定。
申请公布号 CN101419709B 申请公布日期 2012.07.18
申请号 CN200810239085.3 申请日期 2008.12.08
申请人 北京航空航天大学 发明人 周富强;张伟;张广军
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 代理人
主权项 1.一种用于摄像机标定的平面靶标特征点的自动匹配方法,其特征在于,所说的特征点的匹配是指平面靶标特征点的图像坐标和已知世界坐标的对应,利用的是特征点集的邻接关系和局部几何位置关系完成呈方格阵列分布的平面靶标特征点的匹配,具体步骤是:1.1、设平面靶标上预先设置有M×M个呈方格阵列分布的特征点,依据靶标平面上建立的世界坐标系,特征点的世界坐标已知,其排列顺序固定;提取平面靶标图像中的特征点,得到特征点的图像坐标<img file="FSB00000319309000011.GIF" wi="283" he="94" />(i<sub>p</sub>=1,…,N<sub>p</sub>),N<sub>p</sub>为提取的特征点数量,N<sub>p</sub>≤M<sup>2</sup>,<img file="FSB00000319309000012.GIF" wi="63" he="55" />构成的特征点集合用A表示;1.2、对A中的<img file="FSB00000319309000013.GIF" wi="62" he="54" />进行Delaunay三角化处理,得到三角形<img file="FSB00000319309000014.GIF" wi="61" he="72" />(j<sub>d</sub>=1,…,N<sub>d</sub>),N<sub>d</sub>为三角形数量,<img file="FSB00000319309000015.GIF" wi="63" he="74" />构成的集合用T<sub>d</sub>表示,T<sub>d</sub>为N<sub>d</sub>×3维矩阵,每一行存放构成该三角形的三个特征点在A中的索引序号i<sub>p</sub>;1.3、定义A中的对应平面靶标上一个方格的4个相邻特征点构成一个数据单元,用E表示,一个数据单元被Delaunay三角化得到的两个三角形为规则三角形;对步骤1.2获得的T<sub>d</sub>中的三角形<img file="FSB00000319309000016.GIF" wi="62" he="71" />进行预处理,去掉其中的非规则三角形得到规则三角形<img file="FSB00000319309000017.GIF" wi="61" he="82" />(j<sub>g</sub>=1,…,N<sub>g</sub>),N<sub>g</sub>为规则三角形数量,N<sub>g</sub>≤N<sub>d</sub>,<img file="FSB00000319309000018.GIF" wi="60" he="82" />构成的集合用T<sub>g</sub>表示;1.4、定义A中对应平面靶标上最外围的特征点为边界点,边界点由角点和边缘点组成,角点的相邻两个边界点对应平面靶标上的特征点不在同一行或同一列,边缘点的相邻两个边界点对应平面靶标上的特征点属于同一列或者同一行;定义f为特征点在T<sub>g</sub>中作为节点出现的次数;1.5、将f=1或f=2且不是边缘点的特征点定义为候选起始点,候选起始点都属于角点,用<img file="FSB00000319309000019.GIF" wi="60" he="67" />(i<sub>s</sub>=1,…,N<sub>s</sub>)表示,N<sub>s</sub>为候选起始点数量,4≤N<sub>s</sub><N<sub>p</sub>;对<img file="FSB000003193090000110.GIF" wi="59" he="67" />进行排序,将f=1对应的候选起始点排在前面,<img file="FSB000003193090000111.GIF" wi="61" he="66" />构成的集合用A<sub>s</sub>表示;1.6、定义M×M的矩阵C,用来存放排好序的特征点索引,各元素初始赋值为0,采用循环检验方法确定A中的起始点p<sub>s</sub>,对应平面靶标上第一行的特征点<img file="FSB000003193090000112.GIF" wi="108" he="48" />(j<sub>r1</sub>=1,…,N<sub>r1</sub>)和对应平面靶标上第一列的特征点<img file="FSB000003193090000113.GIF" wi="109" he="47" />(i<sub>c1</sub>=1,…,N<sub>c1</sub>),N<sub>r1</sub>和N<sub>c1</sub>分别为对应平面靶标上第一行的特征点数量和对应平面靶标上第一列的特征点数量;<img file="FSB000003193090000114.GIF" wi="109" he="49" />构成的特征点子集用<img file="FSB000003193090000115.GIF" wi="61" he="59" />表示,<img file="FSB000003193090000116.GIF" wi="109" he="49" />构成的特征点子集用<img file="FSB000003193090000117.GIF" wi="59" he="60" />表示;将<img file="FSB000003193090000118.GIF" wi="61" he="60" />中搜索到的特征点索引依次保存到C的第1行中;1.7、搜索对应平面靶标上第i行的特征点<img file="FSB000003193090000119.GIF" wi="99" he="50" />(i=2,…M,j<sub>ri</sub>=1,…,N<sub>ri</sub>),N<sub>ri</sub>为对应平面靶标上第i行的特征点数量,<img file="FSB000003193090000120.GIF" wi="101" he="48" />构成的特征点子集用<img file="FSB000003193090000121.GIF" wi="61" he="61" />表示,将<img file="FSB000003193090000122.GIF" wi="60" he="61" />中搜索到的特征点索引依次保存到C的第i行中;搜索完所有的行后,排好序的特征点索引都保存C中,没有提取到的特征点所对应的位置索引仍为0,根据C中的各元素的索引获得特征点的靶标世界坐标的序号。
地址 100083 北京市海淀区学院路37号
您可能感兴趣的专利