CN101814186A - 利用曲线拟合校正摄像机径向畸变的方法 - Google Patents

利用曲线拟合校正摄像机径向畸变的方法 Download PDF

Info

Publication number
CN101814186A
CN101814186A CN201010301177A CN201010301177A CN101814186A CN 101814186 A CN101814186 A CN 101814186A CN 201010301177 A CN201010301177 A CN 201010301177A CN 201010301177 A CN201010301177 A CN 201010301177A CN 101814186 A CN101814186 A CN 101814186A
Authority
CN
China
Prior art keywords
distortion
coordinate
image
radial distortion
fitting
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
CN201010301177A
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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201010301177A priority Critical patent/CN101814186A/zh
Publication of CN101814186A publication Critical patent/CN101814186A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

一种图像处理技术领域的利用曲线拟合校正摄像机径向畸变的方法,采用棋盘格作为模板,采集该模板的图像并寻找畸变中心,然后检测棋盘格上直线的投影点集并将图像平面的原点移至畸变的中心,利用圆拟合方法得到原直线投影的圆弧的方程,从而得到径向畸变参数,最后利用传统的基于平面标定方法估计其余摄像机内部参数。本发明避免了迭代方法的复杂性以及陷入局部最优值的局限性,而是采用曲线拟合方法求得径向畸变参数,方法简便并且有效;另一方面,本发明方法将径向畸变参数和其余摄像机内部参数的估计分离开来,比现有方法更加灵活。

Description

利用曲线拟合校正摄像机径向畸变的方法
技术领域
本发明涉及的是一种图像处理技术领域的方法,具体是一种利用曲线拟合校正摄像机径向畸变的方法。
背景技术
摄像机镜头径向畸变广泛存在于广角摄像机、鱼眼摄像机以及短焦距的廉价摄像机中,因此径向畸变是分析数字图像要考虑的一个重要问题。
目前已有很多方法用于恢复畸变参数。R.Swaminathan等于2000年在《IEEETransaction on Pattern Analysis and Machine Intelligence》(2000年国际电子与电气工程学会模式分析与机器智能汇刊)上发表的论文”Nonmetric calibration ofwide-angle lenses and polycameras”(广角镜头和多摄像机的非度量标定)基于所谓的铅垂线理论进行畸变校正,即空间中的直线在理想针孔投影变换下仍然投影为直线。根据这个原理,首先检测出已知空间直线的投影点集,然后创建目标函数使得投影点集拟合的曲线与直线的偏离最小。这是一个迭代方法直至畸变误差小于某个阈值。
相类似的技术包括Z.Zhang于2000年在《IEEE Transaction on Pattern Analysis andMachine Intelligence》(2000年国际电子与电气工程学会模式分析与机器智能汇刊)上发表的论文”A flexible new technique for camera calibration”(一种摄像机标定的新技术)使用人工精细制成的标定网格,利用优化方法将镜头畸变参数和其它摄像机内部外部参数同时标定出来。这种非线性迭代方法常常陷入局部最优,从而不能收敛到全局最优值,并且需要选择初始值和确定终止依据,计算过程比较繁琐。
进一步检索发现,G.Stein于1997年在《1997IEEE Computer Society Conference onComputer Vision and Pattern Recognition》(1997年计算机视觉和模式识别国际会议)上发表的文章”Lens distortion calibration using point correspondence”(利用对应点进行镜头畸变标定)采用多幅视图中的对应点来求取镜头畸变参数。虽然这种方法对场景没有要求,但是在某些情况下获得多幅视图是比较困难的。
发明内容
本发明针对现有技术存在的上述不足,提供一种利用曲线拟合校正摄像机径向畸变的方法,能够简单有效地标定出镜头的径向畸变参数。
本发明是通过以下技术方案实现的,本发明采用棋盘格作为模板,采集该模板的图像并寻找畸变中心,然后检测棋盘格上直线的投影点集并将图像平面的原点移至畸变的中心,利用圆拟合方法得到原直线投影的圆弧的方程,从而得到径向畸变参数,最后利用传统的基于平面标定方法估计其余摄像机内部参数。
本发明包括以下步骤:
步骤一、在一张白纸上打印黑白棋盘格,然后将其贴在一个平板上制成模板;
步骤二、采用待测摄像机拍摄一副该模板的图像并转化为位图格式作为初始图像,然后利用harris角点检测方法得到畸变中心坐标;
所述的初始图像的角点的空间坐标xc及其图像坐标xd与径向基础矩阵Fr满足:xd TFrxc=0,其中:xd T表示xd坐标矢量的转置,径向基础矩阵Fr与初始图像的畸变中心e满足:eTFr=0,其中:eT表示矢量e的转置,利用这一约束最终得到畸变中心坐标e;
步骤三、将初始图像减去畸变中心坐标得到修正图像,再使用霍夫变换求取修正图像中棋盘格的直线,得到直线的投影点集坐标。
步骤四、将同一直线的投影点集合利用圆拟合关系拟合为一个圆的参数方程a(x2+y2)+bx+cy+d=0,其中:x、y为步骤三中得到的投影点集的坐标,a、b、c、d为圆参数。
步骤五、根据直线投影的圆弧方程与畸变参数的关系获得畸变参数,然后通过畸变参数采用反投影方法对原始畸变图像进行无畸变恢复,得到理想的无畸变图像;
所述的圆弧方程与畸变参数的关系式为:其中:k1表示畸变参数,(x0,y0)表示步骤四拟合所得圆的圆心坐标,r为圆的半径。
步骤六、重复步骤五直至求出初始图像中的每一条直线的畸变参数后对其求平均值作为初始图像的径向畸变参数,从而完成对摄像机径向畸变的校准。
与现有技术相比,本发明避免了迭代方法的复杂性以及陷入局部最优值的局限性,而是采用曲线拟合方法求得径向畸变参数,方法简便并且有效;另一方面,本发明方法将径向畸变参数和其余摄像机内部参数的估计分离开来,比现有方法更加灵活。
附图说明
图1为本发明中直线在径向畸变下投影示意图。
图2为本发明中所采用的棋盘格模板示意图。
图3为实施例示意图;
其中:图3a为带有径向畸变的图像,图3b为实施后图像。
具体实施方式
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,为实施例中直线在径向畸变下的投影示意图。图1中,oxy为图像坐标系,直线L为一空间直线在针孔摄像机模型下的理想投影,曲线C为同一空间直线在有径向畸变的投影模型下得到的投影。
本实施例的具体实施过程如下:
1.首先制作棋盘格模板:如图2所示,准备一张白纸,在这张白纸上打印黑白棋盘格,棋盘格包含9×7单元格,每格边长为30mm,然后将这张白纸贴在一个平板上,完成模板的制作。
2.本实施例使用的摄像机为SCOR-14SOM数字摄像机,对焦距为6mm的镜头进行畸变校正。首先采集一幅棋盘格模板的图像并转化为bmp格式的初始图像,利用harris角点检测方法得到初始图像的角点;所述的初始图像的角点的空间坐标xc及其图像坐标xd与径向基础矩阵Fr满足一关系式xd TFrxc=0,其中xd T表示坐标矢量xd的转置,从而很容易解得径向基础矩阵Fr;所述的矩阵Fr与初始图像的畸变中心e满足约束eTFr=0,其中eT表示e的转置,从而最终得到畸变中心坐标。
3.将原始图像减去步骤2中求得的畸变中心坐标得到一幅新的图像,利用霍夫变换检测新图像中的直线,得到直线的投影点集坐标。
4.根据圆拟合关系式:a(x2+y2)+bx+cy+d=0,将投影点集坐标拟合到一个圆上,得到圆的参数方程,其中:x、y为步骤3中得到的投影点集的坐标,a、b、c、d为圆参数。
5.如图3a所示,对于每条直线投影,根据原直线投影的圆弧方程与畸变参数的关系式:其中k1表示畸变参数,(x0,y0)表示步骤4拟合所得圆的圆心坐标,r为圆的半径,获得畸变参数,然后通过畸变参数采用反投影方法对原始畸变图像进行无畸变恢复,得到理想的无畸变图像。
6.如图3b所示,重复步骤五直至求出初始图像中的每一条直线的畸变参数后对其求平均值作为初始图像的径向畸变参数,从而完成对摄像机径向畸变的校准。

Claims (4)

1.一种利用曲线拟合校正摄像机径向畸变的方法,其特征在于:采用棋盘格作为模板,采集该模板的图像并寻找畸变中心,然后检测棋盘格上直线的投影点集并将图像平面的原点移至畸变的中心,利用圆拟合方法得到原直线投影的圆弧的方程,从而得到径向畸变参数,最后利用传统的基于平面标定方法估计其余摄像机内部参数。
2.根据权利要求1所述的利用曲线拟合校正摄像机径向畸变的方法,其特征是,所述方法具体包括以下步骤:
步骤一、在一张白纸上打印黑白棋盘格,然后将其贴在一个平板上制成模板;
步骤二、采用待测摄像机拍摄一副该模板的图像并转化为位图格式作为初始图像,然后利用harris角点检测方法得到畸变中心坐标;
步骤三、将初始图像减去畸变中心坐标得到修正图像,再使用霍夫变换求取修正图像中棋盘格的直线,得到直线的投影点集坐标;
步骤四、将同一直线的投影点集合利用圆拟合关系拟合为一个圆的参数方程a(x2+y2)+bx+cy+d=0,其中:x、y为步骤三中得到的投影点集的坐标,a、b、c、d为圆参数;
步骤五、根据直线投影的圆弧方程与畸变参数的关系获得畸变参数,然后通过畸变参数采用反投影方法对原始畸变图像进行无畸变恢复,得到理想的无畸变图像;
步骤六、重复步骤五直至求出初始图像中的每一条直线的畸变参数后对其求平均值作为初始图像的径向畸变参数,从而完成对摄像机径向畸变的校准。
3.根据权利要求2所述的利用曲线拟合校正摄像机径向畸变的方法,其特征是,所述的初始图像的角点的空间坐标xc及其图像坐标xd与径向基础矩阵Fr满足:xd TFrxc=0,其中:xd T表示xd坐标矢量的转置,径向基础矩阵Fr与初始图像的畸变中心e满足:eTFr=0,其中:eT表示矢量e的转置,利用这一约束最终得到畸变中心坐标e。
4.根据权利要求2所述的利用曲线拟合校正摄像机径向畸变的方法,其特征是,所述的圆弧方程与畸变参数的关系式为:
Figure F201010301177720100204C000021
其中:k1表示畸变参数,(x0,y0)表示步骤四拟合所得圆的圆心坐标,r为圆的半径。
CN201010301177A 2010-02-04 2010-02-04 利用曲线拟合校正摄像机径向畸变的方法 Pending CN101814186A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010301177A CN101814186A (zh) 2010-02-04 2010-02-04 利用曲线拟合校正摄像机径向畸变的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010301177A CN101814186A (zh) 2010-02-04 2010-02-04 利用曲线拟合校正摄像机径向畸变的方法

Publications (1)

Publication Number Publication Date
CN101814186A true CN101814186A (zh) 2010-08-25

Family

ID=42621433

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010301177A Pending CN101814186A (zh) 2010-02-04 2010-02-04 利用曲线拟合校正摄像机径向畸变的方法

Country Status (1)

Country Link
CN (1) CN101814186A (zh)

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102750697A (zh) * 2012-06-08 2012-10-24 华为技术有限公司 一种参数标定方法及装置
CN103116889A (zh) * 2013-02-05 2013-05-22 海信集团有限公司 一种定位方法及电子设备
CN103124334A (zh) * 2012-12-19 2013-05-29 四川九洲电器集团有限责任公司 一种镜头畸变校正的方法
CN103177439A (zh) * 2012-11-26 2013-06-26 惠州华阳通用电子有限公司 一种基于黑白格角点匹配的自动标定方法
CN103258328A (zh) * 2013-04-27 2013-08-21 浙江大学 一种宽视场镜头的畸变中心定位方法
CN103377474A (zh) * 2012-04-27 2013-10-30 比亚迪股份有限公司 镜头阴影校正系数确定方法、镜头阴影校正方法及装置
CN103617615A (zh) * 2013-11-27 2014-03-05 华为技术有限公司 径向畸变参数获取方法及获取装置
CN104038755A (zh) * 2014-04-30 2014-09-10 惠州华阳通用电子有限公司 摄像头畸变中心点测试装置及识别方法
CN104820973A (zh) * 2015-05-07 2015-08-05 河海大学 畸变曲线弧度检测模板的图像校正方法
CN105023265A (zh) * 2014-04-29 2015-11-04 东北大学 鱼眼镜头下的棋盘格角点自动检测方法
CN106327441A (zh) * 2016-08-19 2017-01-11 北京交通大学 图像径向畸变的自动矫正方法及系统
CN104574419B (zh) * 2015-01-28 2017-06-09 深圳市安健科技股份有限公司 镜头畸变参数标定方法及系统
CN107256570A (zh) * 2017-06-12 2017-10-17 浙江理工大学 一种基于最优化估计的摄像机外参数标定方法
CN109242917A (zh) * 2018-10-24 2019-01-18 南昌航空大学 一种基于棋盘格的相机分辨率标定方法
CN109801300A (zh) * 2017-11-16 2019-05-24 北京百度网讯科技有限公司 棋盘格角点的坐标提取方法、装置、设备及计算机可读存储介质
CN109949230A (zh) * 2017-12-20 2019-06-28 艾迪普(北京)文化科技股份有限公司 基于图像识别的广角摄像机畸变渲染方法
CN109978959A (zh) * 2019-03-29 2019-07-05 北京经纬恒润科技有限公司 一种摄像头径向畸变修正参数标定方法、装置及系统
CN110246079A (zh) * 2019-05-23 2019-09-17 上海交通大学 基于b样条曲面拟合的摄像头畸变矫正方法、系统及介质
CN111415295A (zh) * 2020-03-17 2020-07-14 东南数字经济发展研究院 一种倾斜摄影三维模型的拍摄分辨率正射图生成方法
CN111707440A (zh) * 2020-06-16 2020-09-25 中国人民解放军国防科技大学 一种可获得连续多幅微秒级时间相关的流场实验装置及方法
CN111754568A (zh) * 2020-06-19 2020-10-09 上海保隆汽车科技股份有限公司 一种标定图案、标定方法及其标定板装置
CN112284292A (zh) * 2020-12-01 2021-01-29 吉林大学 基于等摇臂机构的汽车形貌测量系统畸变系数测量仪
CN115170674A (zh) * 2022-07-20 2022-10-11 禾多科技(北京)有限公司 基于单张图像的相机主点标定方法、装置、设备和介质
CN116245748A (zh) * 2022-12-23 2023-06-09 珠海视熙科技有限公司 环视镜头的畸变校正方法、装置、设备、系统及存储介质

Cited By (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103377474A (zh) * 2012-04-27 2013-10-30 比亚迪股份有限公司 镜头阴影校正系数确定方法、镜头阴影校正方法及装置
CN103377474B (zh) * 2012-04-27 2016-08-03 比亚迪股份有限公司 镜头阴影校正系数确定方法、镜头阴影校正方法及装置
CN102750697A (zh) * 2012-06-08 2012-10-24 华为技术有限公司 一种参数标定方法及装置
CN102750697B (zh) * 2012-06-08 2014-08-20 华为技术有限公司 一种参数标定方法及装置
WO2013182080A1 (zh) * 2012-06-08 2013-12-12 华为技术有限公司 一种参数标定方法及装置
CN103177439B (zh) * 2012-11-26 2015-10-28 惠州华阳通用电子有限公司 一种基于黑白格角点匹配的自动标定方法
CN103177439A (zh) * 2012-11-26 2013-06-26 惠州华阳通用电子有限公司 一种基于黑白格角点匹配的自动标定方法
CN103124334A (zh) * 2012-12-19 2013-05-29 四川九洲电器集团有限责任公司 一种镜头畸变校正的方法
CN103124334B (zh) * 2012-12-19 2015-10-21 四川九洲电器集团有限责任公司 一种镜头畸变校正的方法
CN103116889A (zh) * 2013-02-05 2013-05-22 海信集团有限公司 一种定位方法及电子设备
CN103258328A (zh) * 2013-04-27 2013-08-21 浙江大学 一种宽视场镜头的畸变中心定位方法
CN103258328B (zh) * 2013-04-27 2015-12-09 浙江大学 一种宽视场镜头的畸变中心定位方法
CN103617615B (zh) * 2013-11-27 2016-08-17 华为技术有限公司 径向畸变参数获取方法及获取装置
CN103617615A (zh) * 2013-11-27 2014-03-05 华为技术有限公司 径向畸变参数获取方法及获取装置
US9691137B2 (en) 2013-11-27 2017-06-27 Huawei Technologies Co., Ltd. Radial distortion parameter acquiring method and apparatus
CN105023265A (zh) * 2014-04-29 2015-11-04 东北大学 鱼眼镜头下的棋盘格角点自动检测方法
CN104038755A (zh) * 2014-04-30 2014-09-10 惠州华阳通用电子有限公司 摄像头畸变中心点测试装置及识别方法
CN104038755B (zh) * 2014-04-30 2016-09-21 惠州华阳通用电子有限公司 摄像头畸变中心点测试装置及识别方法
CN104574419B (zh) * 2015-01-28 2017-06-09 深圳市安健科技股份有限公司 镜头畸变参数标定方法及系统
CN104820973A (zh) * 2015-05-07 2015-08-05 河海大学 畸变曲线弧度检测模板的图像校正方法
CN104820973B (zh) * 2015-05-07 2017-10-03 河海大学 畸变曲线弧度检测模板的图像校正方法
CN106327441B (zh) * 2016-08-19 2018-11-06 北京交通大学 图像径向畸变的自动矫正方法及系统
CN106327441A (zh) * 2016-08-19 2017-01-11 北京交通大学 图像径向畸变的自动矫正方法及系统
CN107256570A (zh) * 2017-06-12 2017-10-17 浙江理工大学 一种基于最优化估计的摄像机外参数标定方法
CN109801300A (zh) * 2017-11-16 2019-05-24 北京百度网讯科技有限公司 棋盘格角点的坐标提取方法、装置、设备及计算机可读存储介质
CN109949230A (zh) * 2017-12-20 2019-06-28 艾迪普(北京)文化科技股份有限公司 基于图像识别的广角摄像机畸变渲染方法
CN109242917A (zh) * 2018-10-24 2019-01-18 南昌航空大学 一种基于棋盘格的相机分辨率标定方法
CN109978959A (zh) * 2019-03-29 2019-07-05 北京经纬恒润科技有限公司 一种摄像头径向畸变修正参数标定方法、装置及系统
CN110246079B (zh) * 2019-05-23 2020-11-06 上海交通大学 基于b样条曲面拟合的摄像头畸变矫正方法、系统及介质
CN110246079A (zh) * 2019-05-23 2019-09-17 上海交通大学 基于b样条曲面拟合的摄像头畸变矫正方法、系统及介质
CN111415295A (zh) * 2020-03-17 2020-07-14 东南数字经济发展研究院 一种倾斜摄影三维模型的拍摄分辨率正射图生成方法
CN111415295B (zh) * 2020-03-17 2024-01-12 东南数字经济发展研究院 一种倾斜摄影三维模型的拍摄分辨率正射图生成方法
CN111707440A (zh) * 2020-06-16 2020-09-25 中国人民解放军国防科技大学 一种可获得连续多幅微秒级时间相关的流场实验装置及方法
CN111707440B (zh) * 2020-06-16 2021-11-26 中国人民解放军国防科技大学 一种可获得连续多幅微秒级时间相关的流场实验装置及方法
CN111754568A (zh) * 2020-06-19 2020-10-09 上海保隆汽车科技股份有限公司 一种标定图案、标定方法及其标定板装置
CN112284292A (zh) * 2020-12-01 2021-01-29 吉林大学 基于等摇臂机构的汽车形貌测量系统畸变系数测量仪
CN112284292B (zh) * 2020-12-01 2021-07-30 吉林大学 基于等摇臂机构的汽车形貌测量系统畸变系数测量仪
CN115170674A (zh) * 2022-07-20 2022-10-11 禾多科技(北京)有限公司 基于单张图像的相机主点标定方法、装置、设备和介质
CN116245748A (zh) * 2022-12-23 2023-06-09 珠海视熙科技有限公司 环视镜头的畸变校正方法、装置、设备、系统及存储介质
CN116245748B (zh) * 2022-12-23 2024-04-26 珠海视熙科技有限公司 环视镜头的畸变校正方法、装置、设备、系统及存储介质

Similar Documents

Publication Publication Date Title
CN101814186A (zh) 利用曲线拟合校正摄像机径向畸变的方法
CN109215082B (zh) 一种相机参数标定方法、装置、设备及系统
CN104240236B (zh) 一种鱼眼镜头标定后鱼眼图像校正的方法
CN106204544B (zh) 一种自动提取图像中标志点位置及轮廓的方法和系统
CN103278138B (zh) 一种复杂结构薄部件三维位置及姿态的测量方法
TWI491849B (zh) 應用複數攝影裝置之物件定位方法
CN100428805C (zh) 仅用一幅平面标定物图像的摄像机标定方法
CN105096269A (zh) 基于畸变直线结构检测的图像径向畸变矫正的方法及系统
CN108447022B (zh) 基于单个固定摄像头图像序列的运动目标拼接方法
Chaudhury et al. Auto-rectification of user photos
CN102075679A (zh) 一种图像采集方法和装置
CN103971378A (zh) 一种混合视觉系统中全景图像的三维重建方法
CN107680139A (zh) 一种远心双目立体视觉测量系统的通用性标定方法
CN102436660A (zh) 3d相机图像自动校正方法及装置
CN111461963B (zh) 一种鱼眼图像拼接方法及装置
CN106952219B (zh) 一种基于外参数修正鱼眼像机的图像生成方法
CN102164269A (zh) 全景监控方法及装置
CN104869316B (zh) 一种多目标的摄像方法及装置
CN112949478A (zh) 基于云台相机的目标检测方法
CN106060399A (zh) 一种双摄像头自动aa的方法及装置
CN106886976B (zh) 一种基于内参数修正鱼眼像机的图像生成方法
CA3233222A1 (en) Method, apparatus and device for photogrammetry, and storage medium
CN106851103A (zh) 基于双摄像头系统的自动对焦方法和装置
CN116957987A (zh) 一种多目极线校正方法、装置、计算机设备及存储介质
JP2011232965A (ja) 平面領域検出装置及びステレオカメラシステム

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Open date: 20100825