CN106097319A - 一种用于相机检校的黑白棋盘格角点提取方法 - Google Patents

一种用于相机检校的黑白棋盘格角点提取方法 Download PDF

Info

Publication number
CN106097319A
CN106097319A CN201610392800.1A CN201610392800A CN106097319A CN 106097319 A CN106097319 A CN 106097319A CN 201610392800 A CN201610392800 A CN 201610392800A CN 106097319 A CN106097319 A CN 106097319A
Authority
CN
China
Prior art keywords
point
angle
black
angle point
white
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610392800.1A
Other languages
English (en)
Inventor
王海涛
邱儒琼
李强
洪亮
张露林
余晓敏
陈晓茜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hubei Basic Geographic Information Center (hubei Beidou Satellite Navigation Application Technology Research Institute)
Original Assignee
Hubei Basic Geographic Information Center (hubei Beidou Satellite Navigation Application Technology Research Institute)
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hubei Basic Geographic Information Center (hubei Beidou Satellite Navigation Application Technology Research Institute) filed Critical Hubei Basic Geographic Information Center (hubei Beidou Satellite Navigation Application Technology Research Institute)
Priority to CN201610392800.1A priority Critical patent/CN106097319A/zh
Publication of CN106097319A publication Critical patent/CN106097319A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20112Image segmentation details
    • G06T2207/20164Salient point detection; Corner detection

Landscapes

  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种用于相机检校的黑白棋盘格角点提取方法,包括以下步骤:绘制黑白棋盘格图,使用棋盘格作为平面控制条件;建立平面控制点领域特性及编号,平面控制点称为X角点;提取X角点的子像素。本发明的用于相机检校的黑白棋盘格角点提取方法基于其角点的对称性、尺度不变性和领域特性,提出一种有效的角点提取算法。同时针对格网边缘的特性,提出双边缘的匹配算法提取高精度的角点。本发明提出的相机检校的黑白棋盘格角点提取方法具有很强的鲁棒性的黑白棋盘格角点提取算法,并可以高精度的定位角点。

Description

一种用于相机检校的黑白棋盘格角点提取方法
技术领域
本发明属于相机检校技术领域,具体涉及一种用于相机检校的黑白棋盘格角点提取方法。
背景技术
随着计算机的出现,通过计算机图像处理和射影几何获取被摄物体几何以及纹理信息已经成为了数字摄影测量与计算机视觉界的研究热点,无论是最早的胶片相机还是现在的数码相机,它们获取的影像都存在或大或小的畸变,并且由于相机结构和成像机理的差异导致了畸变的差异。因此,要想利用这些影像来进行高精度的测量,就必须进行必要的相机畸变差测定和补偿,同时也需要精确测定出相机的主距和像主点坐标值等,摄影测量界中称这个过程为相机检校,而计算机视觉界中称之为相机标定,目前对数码相机检校理论的研究已经相当成熟,国内外许多摄影测量和计算机视觉专家也提出了各种行之有效的检校方法。
使用平面控制条件,如平面格网板、聚脂薄膜格网、打印纸张格网、LCD绘制圆等作为相机检校参考对象已被广泛使用,黑白棋盘格更是广泛使用。众多厂商已使用陶瓷、玻璃等材料制作了高精度的棋盘格式标定片。现有技术中已经实现了黑白棋盘格角点的自动提取并实现了基于二维平面模板的标定算法,使用棋盘格作为平面控制条件,首先是格网点的提取,一般使用内角点作为平面控制点。多数学者采用算子提取特征点,对提取特征点进行处理后获得有效格网角点。黑白棋盘格角点作为人造的特征,有其自身的特性(对称性和尺度不变特性),另外提取的X角点与其邻域点有重要几何特性(邻域特性),利用该特性对X角点进行错误检测及编号。提取X角点后,接下来是使用相机的成像特性选取检校方法,从学科的角度可分为计算机视觉方法、摄影测量方法。但现有技术中没有提出一种有效的边缘匹配算法来对相机进行检校。
发明内容
为克服上述现有技术的缺陷,本发明的目的在于提供一种用于相机检校的黑白棋盘格角点提取方法,采用自适应边缘模板和双边缘特性提取高精度的X角点,提高相机检校的精度。
为实现上述目的,本发明采用如下技术方案实现:一种用于相机检校的黑白棋盘格角点提取方法,包括以下步骤:
一、绘制黑白棋盘格图,使用棋盘格作为平面控制条件;
二、建立平面控制点领域特性及编号,平面控制点称为X角点;
三、提取X角点的子像素。
作为上述方案的改进,在本发明的一实施例中,所述绘制黑白棋盘格图,具体包括:
截取黑白棋盘格图中两类X角点的局部图像,并对两类局部图像黑白取反;
提取X角点的对称特性;
提取X角点的尺度特性。
优选地,所述建立平面控制点领域特性及编号,具体包括:
提取的X角点可分为四个外角点(C)、边界角点(E)、内角点(H),以I类X角内角点H为例描述其邻域特性,离内角点H距离最近点的8领域点用N1~N8表达,其中N1~N4为II类点、N5~N8为I类点;内角点(H)与N1~N4点局部影像相关系数小于0,与N5~N8点局部影像相关系数大于0,选择与内角点(H)关系数小于0的领域点N1~N4,满足以下条件:
(3)N1、N2、N3、N4组成的多变形为凸多边形;
(4)直线长度比值
(3)以H为原点,以H与Ni(i=1,2,3,4)和Nj(j=2,3,4,1)中点连线为Y轴,建立的直角坐标系,则N1、N2、N3、N4分布在这个坐标系的四个象限中;
对于边界点E,一定在近边界内角点H的四领域点N1~N4中;
设H与N1~N4点的距离最大值为对于角点C,与其最近的内点H的距离为则Cs小于某一阈值。
作为上述方案的改进,在本发明的一实施例中,所述建立平面控制点领域特性及编号,还包括:
对内角点Hi,选择四领域点构成四边形,以对角点连线建立直线,分别为
为例,在该直线上的内角点Hj满足以下条件:以内角点Hj构建的直线夹角大的直线,构成该直线点的两个邻域点,分布在两侧,依据该条件搜索所有在的内角点Hj,则内角点Hj及Hj对角领域点与夹角小的两个点也在上;
设所有在的角点Hj构成点直线方程为Lx,所有在角点Hj构成点直线方程为Ly,对构建后直线方程进行排序,依据排序后直线方程对X角点进行编号。
进一步的,所述提取X角点的子像素,具体包括:
提取格网点边缘特征,将拍摄的黑白棋盘格影像灰度值作为高度值进行配色显示和三维显示;
边缘匹配,采用函数描述棋盘格网影像黑白过度区域边缘特征,缘灰度值与理论模板曲线的关系可表达为:
F = H 0 + H 1 tan - 1 x K - - - ( 1 - 2 )
其中H0、H1用于描述影像灰度值与理论边缘点辐射变化、K用于调节边缘陡峭程度。
本发明的有益效果在于:本发明的用于相机检校的黑白棋盘格角点提取方法基于其角点的对称性、尺度不变性和领域特性,提出一种有效的角点提取算法。同时针对格网边缘的特性,提出双边缘的匹配算法提取高精度的角点。本发明提出的相机检校的黑白棋盘格角点提取方法具有很强的鲁棒性的黑白棋盘格角点提取算法,并可以高精度的定位角点。
附图说明
图1是本发明的用于相机检校的黑白棋盘格角点提取方法的流程图。
图2是本发明的用于相机检校的黑白棋盘格角点提取方法中的X角点的邻域特性图。
具体实施方式
本发明使用棋盘格作为平面控制条件,首先是格网点的提取,使用内角点,本发明称为X角点,作为平面控制点。黑白棋盘格角点作为人造的特征,有其自身的特性,称为对称性和尺度不变特性,另外提取的X角点与其邻域点有重要几何特性,称为邻域特性,利用该特性对X角点进行错误检测及编号。提取X角点后,接下来是使用相机的成像特性选取检校方法,从学科的角度可分为计算机视觉方法、摄影测量方法。本发明采用摄影测量的光束法平差(Bundle adjustment),该方法可以同时确定目标三维坐标和相机内参数,即自检校(Self-calibration)光束法平差。本发明使用LCD绘制棋盘格,LCD本身也存在误差,包括平面误差及高程误差(LCD的平整度),光束法平差将这些误差作为未知数同时解算。
采用光束法平差进行相机检校,取决观测值的精度,即X角点的提取精度,棋盘格网作为人造特征,格网黑白变换处形成阶跃边缘,边缘匹配可以高精定位边缘,使用提取的边缘直线相交计算子X角点子像素坐标。本发明针对黑白棋盘格网提出采用自适应边缘模板和双边缘特性提取高精度的X角点。
具体地,如图1所示,本发明的用于相机检校的黑白棋盘格角点提取方法包括以下步骤:
S1、绘制黑白棋盘格图,使用棋盘格作为平面控制条件,具体包括:
截取黑白棋盘格图中两类X角点的局部图像,并对两类局部图像黑白取反,每类X角点其原始局部影像与另一类点的反向图像相似。对于X角点,同类X角点局部影像相关系数约等于1,不同类X角点局部影像相关系数约等于-1。
提取X角点的对称特性,设图像的灰度值为F,设黑白棋盘格X角点图像白色部分的灰度值为W、黑色部分影像灰度值为B,图像灰度的均值为像素灰度值对于不同对称方式的关系如表1。
表1X角点对称特性像素值关系
使用对称轴对X角点图像进行变换,变换后的图像表达式及特性见表2;
表2影像变换及特性
提取X角点的尺度特性,使用影像金字塔的方式表达尺度空间。
ρ = ( F I ( i , j ) - F 1 ‾ ) ( F 2 ( i , j ) - F 2 ‾ ) Σ i m Σ j n ( F I ( i , j ) - F 1 ‾ ) 2 Σ i m Σ j n ( F 2 ( i , j ) - F 2 ‾ ) 2 - - - ( 1 - 1 )
依式1-1计算不同尺度影像的相关系数应约等于1。
S2、建立平面控制点领域特性及编号,平面控制点称为X角点,具体包括:
如图2所示,提取的X角点可分为四个外角点(C)、边界角点(E)、内角点(H),以I类X角内角点H为例描述其邻域特性,离内角点H距离最近点的8领域点用N1~N8表达,其中N1~N4为II类点、N5~N8为I类点;内角点(H)与N1~N4点局部影像相关系数小于0,与N5~N8点局部影像相关系数大于0,选择与内角点(H)关系数小于0的领域点N1~N4,满足以下条件:
(1)N1、N2、N3、N4组成的多变形为凸多边形;
(2)直线长度比值
(3)以H为原点,以H与Ni(i=1,2,3,4)和Nj(j=2,3,4,1)中点连线为Y轴,建立的直角坐标系,则N1、N2、N3、N4分布在这个坐标系的四个象限中;
对于边界点E,一定在近边界内角点H的四领域点N1~N4中;
设H与N1~N4点的距离最大值为对于角点C,与其最近的内点H的距离为则Cs小于某一阈值。
对内角点Hi,选择四领域点构成四边形,以对角点连线建立直线,分别为
为例,在该直线上的内角点Hj满足以下条件:以内角点Hj构建的直线夹角大的直线,构成该直线点的两个邻域点,分布在两侧,依据该条件搜索所有在的内角点Hj,则内角点Hj及Hj对角领域点与夹角小的两个点也在上。
设所有在的角点Hj构成点直线方程为Lx,所有在角点Hj构成点直线方程为Ly,对构建后直线方程进行排序,依据排序后直线方程对X角点进行编号。
S3、提取X角点的子像素,具体包括:
提取格网点边缘特征,将拍摄的黑白棋盘格影像灰度值作为高度值进行配色显示和三维显示。
边缘匹配,采用函数描述棋盘格网影像黑白过度区域边缘特征,缘灰度值与理论模板曲线的关系可表达为:
F = H 0 + H 1 tan - 1 x K - - - ( 1 - 2 )
其中H0、H1用于描述影像灰度值与理论边缘点辐射变化、K用于调节边缘陡峭程度。
本发明采用通过不同型号数码相机、不同型号的绘制设备来绘制黑白棋盘格网,基于X角点自身特性(对称性、尺度不变性、领域特性)的角点提取算法具有很强的鲁棒性,提供一种有效的X角点提取算。通过光束法平差验证表明双边缘的边缘匹配算法可以提取高精度X角点,也表明本发明提出边缘匹配算法一种有效的边缘定位算法。
上述实施例揭示的用于相机检校的黑白棋盘格角点提取方法基于其角点的对称性、尺度不变性和领域特性,提出一种有效的角点提取算法。同时针对格网边缘的特性,提出双边缘的匹配算法提取高精度的角点。本发明提出的相机检校的黑白棋盘格角点提取方法具有很强的鲁棒性的黑白棋盘格角点提取算法,并可以高精度的定位角点。
以上所揭露的仅为本发明的较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。

Claims (5)

1.一种用于相机检校的黑白棋盘格角点提取方法,其特征在于:包括以下步骤:
一、绘制黑白棋盘格图,使用棋盘格作为平面控制条件;
二、建立平面控制点领域特性及编号,平面控制点称为X角点;
三、提取X角点的子像素。
2.如权利要求1所述的用于相机检校的黑白棋盘格角点提取方法,其特征在于:所述绘制黑白棋盘格图,具体包括:
截取黑白棋盘格图中两类X角点的局部图像,并对两类局部图像黑白取反;
提取X角点的对称特性;
提取X角点的尺度特性。
3.如权利要求1所述的用于相机检校的黑白棋盘格角点提取方法,其特征在于:所述建立平面控制点领域特性及编号,具体包括:
提取的X角点可分为四个外角点(C)、边界角点(E)、内角点(H),以I类X角内角点H为例描述其邻域特性,离内角点H距离最近点的8领域点用N1~N8表达,其中N1~N4为II类点、N5~N8为I类点;内角点(H)与N1~N4点局部影像相关系数小于0,与N5~N8点局部影像相关系数大于0,选择与内角点(H)关系数小于0的领域点N1~N4,满足以下条件:
(1)N1、N2、N3、N4组成的多变形为凸多边形;
(2)直线长度比值
(3)以H为原点,以H与Ni(i=1,2,3,4)和Nj(j=2,3,4,1)中点连线为Y轴,建立的直角坐标系,则N1、N2、N3、N4分布在这个坐标系的四个象限中;
对于边界点E,一定在近边界内角点H的四领域点N1~N4中;
设H与N1~N4点的距离最大值为对于角点C,与其最近的内点H的距离为则Cs小于某一阈值。
4.如权利要求3所述的用于相机检校的黑白棋盘格角点提取方法,其特征在于:所述建立平面控制点领域特性及编号,还包括:
对内角点Hi,选择四领域点构成四边形,以对角点连线建立直线,分别为
为例,在该直线上的内角点Hj满足以下条件:以内角点Hj构建的直线夹角大的直线,构成该直线点的两个邻域点,分布在两侧,依据该条件搜索所有在的内角点Hj,则内角点Hj及Hj对角领域点与夹角小的两个点也在上;
设所有在的角点Hj构成点直线方程为Lx,所有在角点Hj构成点直线方程为Ly,对构建后直线方程进行排序,依据排序后直线方程对X角点进行编号。
5.如权利要求1所述的用于相机检校的黑白棋盘格角点提取方法,其特征在于:所述提取X角点的子像素,具体包括:
提取格网点边缘特征,将拍摄的黑白棋盘格影像灰度值作为高度值进行配色显示和三维显示;
边缘匹配,采用函数描述棋盘格网影像黑白过度区域边缘特征,缘灰度值与理论模板曲线的关系可表达为:
F = H 0 + H 1 tan - 1 x K - - - ( 1 - 2 )
其中H0、H1用于描述影像灰度值与理论边缘点辐射变化、K用于调节边缘陡峭程度。
CN201610392800.1A 2016-06-04 2016-06-04 一种用于相机检校的黑白棋盘格角点提取方法 Pending CN106097319A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610392800.1A CN106097319A (zh) 2016-06-04 2016-06-04 一种用于相机检校的黑白棋盘格角点提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610392800.1A CN106097319A (zh) 2016-06-04 2016-06-04 一种用于相机检校的黑白棋盘格角点提取方法

Publications (1)

Publication Number Publication Date
CN106097319A true CN106097319A (zh) 2016-11-09

Family

ID=57448235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610392800.1A Pending CN106097319A (zh) 2016-06-04 2016-06-04 一种用于相机检校的黑白棋盘格角点提取方法

Country Status (1)

Country Link
CN (1) CN106097319A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106846412A (zh) * 2017-01-23 2017-06-13 上海兴芯微电子科技有限公司 一种棋盘格角点检测方法及装置
CN108895959A (zh) * 2018-04-27 2018-11-27 电子科技大学 一种基于亚像素的相机标定板角点计算方法
CN109801300A (zh) * 2017-11-16 2019-05-24 北京百度网讯科技有限公司 棋盘格角点的坐标提取方法、装置、设备及计算机可读存储介质
CN111540028A (zh) * 2020-07-07 2020-08-14 昆山丘钛微电子科技有限公司 一种棋盘格测试图的绘制方法、装置、设备及介质
CN112465760A (zh) * 2020-11-19 2021-03-09 深圳惠牛科技有限公司 棋盘格外角点识别方法、装置、设备及存储介质
CN114596341A (zh) * 2022-02-16 2022-06-07 合肥工业大学 面向大视场运动目标的多相机高精度三维位姿跟踪方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102095370A (zh) * 2010-11-22 2011-06-15 北京航空航天大学 三x组合标记的检测识别方法
WO2012113732A1 (en) * 2011-02-25 2012-08-30 Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. Determining model parameters based on transforming a model of an object

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102095370A (zh) * 2010-11-22 2011-06-15 北京航空航天大学 三x组合标记的检测识别方法
WO2012113732A1 (en) * 2011-02-25 2012-08-30 Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. Determining model parameters based on transforming a model of an object

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
VISHVJIT S.NALWA: "On Detecting Edges", 《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE 》 *
孟偲等: "具有方向特性的X 角点的亚像素检测定位", 《北京航空航天大学学报》 *
张颖: "基于改进的SUSAN法的摄像机线性标定方法", 《计算机应用》 *
袁红照等: "摄像机标定角点点阵的多尺度检测和亚像素定位", 《河南大学学报(自然科学版)》 *
韦光等: "基于双曲正切的亚像素边缘检测方法及评价", 《工程与试验》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106846412A (zh) * 2017-01-23 2017-06-13 上海兴芯微电子科技有限公司 一种棋盘格角点检测方法及装置
CN106846412B (zh) * 2017-01-23 2020-03-13 上海兴芯微电子科技有限公司 一种棋盘格角点检测方法及装置
CN109801300A (zh) * 2017-11-16 2019-05-24 北京百度网讯科技有限公司 棋盘格角点的坐标提取方法、装置、设备及计算机可读存储介质
CN109801300B (zh) * 2017-11-16 2020-08-28 北京百度网讯科技有限公司 棋盘格角点的坐标提取方法、装置、设备及计算机可读存储介质
CN108895959A (zh) * 2018-04-27 2018-11-27 电子科技大学 一种基于亚像素的相机标定板角点计算方法
CN108895959B (zh) * 2018-04-27 2020-03-27 电子科技大学 一种基于亚像素的相机标定板角点计算方法
CN111540028A (zh) * 2020-07-07 2020-08-14 昆山丘钛微电子科技有限公司 一种棋盘格测试图的绘制方法、装置、设备及介质
CN112465760A (zh) * 2020-11-19 2021-03-09 深圳惠牛科技有限公司 棋盘格外角点识别方法、装置、设备及存储介质
CN114596341A (zh) * 2022-02-16 2022-06-07 合肥工业大学 面向大视场运动目标的多相机高精度三维位姿跟踪方法

Similar Documents

Publication Publication Date Title
CN106097319A (zh) 一种用于相机检校的黑白棋盘格角点提取方法
CN109269430B (zh) 基于深度提取模型的多株立木胸径被动测量方法
CN104331896B (zh) 一种基于深度信息的系统标定方法
CN106340044B (zh) 摄像机外参自动标定方法及标定装置
CN105627992B (zh) 一种快速高精度非接触测绘古建筑的方法
CN106091984B (zh) 一种基于线激光的三维点云数据获取方法
CN115170669B (zh) 基于边缘特征点集配准的识别定位方法及系统、存储介质
CN110443879B (zh) 一种基于神经网络的透视误差补偿方法
CN105823416A (zh) 多相机测量物体的方法和装置
CN102376089A (zh) 一种标靶校正方法及系统
CN101667303A (zh) 一种基于编码结构光的三维重建方法
CN106952262B (zh) 一种基于立体视觉的船板加工精度分析方法
US20170199099A1 (en) Method and system for measuring lens distortion
CN113920205B (zh) 一种非同轴相机的标定方法
CN107687819A (zh) 一种快速高精度的光条中心亚像素提取方法
Rüther et al. A comparison of close-range photogrammetry to terrestrial laser scanning for heritage documentation
CN109727279A (zh) 一种矢量数据与遥感影像的自动配准方法
CN103424104A (zh) 一种近景大幅面数字摄影测量系统及方法
CN108053370A (zh) 一种基于匹配误差抑制的成像坐标校正方法
CN109145929A (zh) 一种基于sift尺度空间特征信息提取方法
CN115330684A (zh) 基于双目视觉与线结构光的水下构筑物表观缺陷检测方法
CN107941241B (zh) 一种用于航空摄影测量质量评价的分辨率板及其使用方法
CN105115443B (zh) 基于视觉测量技术的水平全视角高精度三维测量方法
CN102999895A (zh) 利用两个同心圆线性求解摄像机内参数
CN112050752B (zh) 一种基于二次投影的投影仪标定方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20161109