CN106780391B - 一种用于全视角三维测量仪光学系统的畸变矫正算法 - Google Patents
一种用于全视角三维测量仪光学系统的畸变矫正算法 Download PDFInfo
- Publication number
- CN106780391B CN106780391B CN201611224243.9A CN201611224243A CN106780391B CN 106780391 B CN106780391 B CN 106780391B CN 201611224243 A CN201611224243 A CN 201611224243A CN 106780391 B CN106780391 B CN 106780391B
- Authority
- CN
- China
- Prior art keywords
- distortion
- image
- distortion correction
- coordinates
- file
- 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.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 38
- 230000003287 optical effect Effects 0.000 title claims abstract description 38
- 238000013507 mapping Methods 0.000 claims abstract description 6
- 238000003384 imaging method Methods 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 10
- 238000003702 image correction Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 3
- 230000004075 alteration Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 23
- 238000004519 manufacturing process Methods 0.000 abstract description 6
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Image Processing (AREA)
- Studio Devices (AREA)
Abstract
本发明公开了一种用于全视角三维测量仪光学系统的畸变矫正算法,所述算法包括如下步骤:一、畸变关系建立:构建标准图像与畸变图像坐标之间的映射关系,生成畸变矫正系数文件或畸变矫正表格文件;二、畸变图像矫正:利用生成的畸变矫正系数文件或畸变矫正表格文件对输入的畸变图像进行矫正,最终获得矫正后的标准图像。该算法可以有效地减小由于光学镜头制造、装调误差所致图像畸变引起的立体测量误差,大大提高系统的测角精度,解决了由于镜头畸变降低图像中物体的几何位置精度,进而产生测量误差的问题。
Description
技术领域
本发明属于光学成像技术领域,涉及一种用于全视角三维测量仪光学系统的畸变矫正算法。
背景技术
全自动影像测量仪能够便捷而快速的进行三维坐标测量,满足现代制造业对尺寸检测提出的更高速、更便捷、更精准的测量需要,是解决制造业发展的一个瓶颈技术。全视角高精度三维测量仪是一款大视场、高分辨率、工作距离范围大的精密测角定位系统,其测角精度要求达到秒级精度。由于全视角三维测量仪采用F-Theta系统、90°大视场广角光学系统,当光学系统的畸变为0.1%时,由畸变产生的角度偏差为2.73′,当光学系统的畸变降低4倍为0.025%时,畸变产生的角度偏差为40″,减小系统畸变可实现系统测角精度不随视场角增大而降低、保证不同视场内测角精度一致,解决由于镜头畸变降低图像中物体的几何位置精度,进而产生的测量误差,可以有效地减小由于光学镜头制造、装调误差所致图像畸变引起的立体测量误差,大大提高系统的测角精度。
发明内容
为了解决由于镜头畸变降低图像中物体的几何位置精度,进而产生测量误差的问题,本发明提供了一种用于全视角三维测量仪光学系统的畸变矫正算法,该算法可以有效地减小由于光学镜头制造、装调误差所致图像畸变引起的立体测量误差,大大提高系统的测角精度。
本发明的目的是通过以下技术方案实现的:
一种用于全视角三维测量仪光学系统的畸变矫正算法,包括如下步骤:
一、畸变关系建立
构建标准图像与畸变图像坐标之间的映射关系,生成畸变矫正系数文件或畸变矫正表格文件;
二、畸变图像矫正
利用生成的畸变矫正系数文件或畸变矫正表格文件对输入的畸变图像进行矫正,最终获得矫正后的标准图像。
本发明具有如下有益效果:
1、针对采用了F-Theta系统的光学镜头,保证镜头结构本身优点的同时,保证系统畸变满足测量要求。
2、有效降低了光学系统畸变,提高了测角精度,使全视角三维测量仪测角精度达到秒级。
附图说明
图1为算法组成示意图;
图2为畸变关系建立流程图;
图3为畸变图像矫正流程图;
图4为可见光光学系统组成图;
图5为低畸变广角成像光学镜头结构图;
图6为系统畸变分析图。
具体实施方式
下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
具体实施方式一:如图1所示,本实施方式提供的用于全视角三维测量仪光学系统的畸变矫正算法由畸变关系建立与畸变图像矫正两部分组成,其中:
所述畸变关系建立是构建标准图像与畸变图像坐标之间的映射关系。如图2所示,具体步骤如下:
1)根据用户提供的相机参数和实际标靶,利用理想成像投影关系计算控制点的标准图像坐标;
2)根据实际畸变后的靶标图像,获取控制点在畸变图像坐标系中的坐标;
3)如果采用函数表达法对畸变关系进行描述,则利用步骤1)、2)获得的控制点标准图像坐标和畸变图像坐标对映射函数中的参数进行拟合,并按照协议生成畸变矫正系数文件;
4)如果采用表格法对畸变关系进行描述,则对步骤1)、2)获得的控制点标准图像坐标和畸变图像坐标插值,并按照协议生成畸变矫正表格文件。
所述畸变图像矫正是利畸变关系建立部分生成的畸变矫正系数文件或畸变矫正表格文件对输入的畸变图像进行矫正,最终获得矫正后的标准图像。如图3所示,具体步骤如下:
1)如果利用畸变矫正系数文件进行图像矫正,则利用文件中的系数计算矫正后标准图像各像元在畸变图像中对应的坐标;
2)如果利用畸变矫正表格文件进行图像矫正,则在文件中的表格中查找矫正后标准图像各像元在畸变图像中对应的坐标;
3)根据矫正后标准图像各像元在畸变图像中的坐标,在用户的输入畸变图像中插值查找这些坐标对应的图像灰度信息;
4)将查找到的灰度信息进行整合,对标准图像中像素坐标超出畸变图像范围的点进行灰度置零处理,生成畸变矫正后的图像。
具体实施方式二:如图4所示,本实施方式提供了用于全视角三维测量仪的可见光光学系统,所述可见光光学系统由低畸变广角成像光学镜头1和线阵CMOS相机2组成,线阵CMOS相机2与低畸变广角成像光学镜头1固定在一起,低畸变广角成像光学镜头1将搜集到的可见光图像清晰成像在线阵CMOS相机2焦面上。
如图5所示,所述的低畸变广角成像光学镜头1采用了F-Theta镜头结构,使空间角解析度与相机元数和尺寸成正比。低畸变广角成像光学镜头1焦距为26.07,是普通广角成像物镜焦距的3倍以上,为长焦广角成像物镜,有利于提高单像素张角的分辨率,实现目标的精确定位。工作波段为480~650nm,视场角为90º,工作距离为3m~100m,像面离焦量为0.2mm,当系统远心度为1°时,在工作距离内主光线产生的测角偏差为27″,当系统远心度为0.1°时,产生的测角偏差为2.7″,最大远心度为0.137°,可实现测量精度不随测量距离的改变而降低,在不同工作距离内测角精度一致,实现秒级测角精度。低畸变广角成像光学镜头1的F数为4, MTF为0视场≥0.4(100lp/mm),0.7视场≥0.3(100lp/mm),低畸变广角成像光学镜头1用于搜集可见光图像,并清晰成像在线阵CMOS相机2焦面上。
所述的线阵CMOS相机2为4×8192分辨率、5μm×5μm像元的大小的线阵传感器,线阵CMOS相机2通过光电转换,生成相应强度的电信号输出,生成可见光图像,供后续图像处理,从而实现大视场、高清晰、远距离目标观测及图像采集。
利用具体实施方式一所述方法对本实施方式的可见光光学系统的畸变进行矫正,其畸变分析图如图6所示。由图6可知:F-Theta镜头的畸变曲线平直,全视场最大畸变为-0.00167%,为系统指标要求畸变的1/3000,极大的减小了由于镜头畸变引起的角度偏差,此时由于系统畸变所引起的最大角度偏差为2.7″,实现了系统的秒级测角精度,解决了由于镜头畸变降低图像中物体的几何位置精度,进而产生的测量误差,可以有效地减小由于光学镜头制造、装调误差所致图像畸变引起的立体测量误差。
Claims (3)
1.一种用于全视角三维测量仪光学系统的畸变矫正方法,其特征在于所述畸变矫正方法适用的全视角三维测量仪光学系统由低畸变广角成像光学镜头和线阵CMOS相机组成,低畸变广角光学镜头将搜集到的可见光图像清晰成像在线阵CMOS相机焦面上;所述低畸变广角成像光学镜头采用了F-Theta镜头结构;所述低畸变广角成像光学镜头的焦距为26.07mm;所述低畸变广角成像光学镜头的工作波段为480~650nm,视场角为90º,工作距离为3m~100m,最大远心度为0.137°;所述低畸变广角成像光学镜头的F数为4, MTF为0视场≥0.4@100lp/mm,0.7视场≥0.3@100lp/mm;所述线阵CMOS相机为4×8192分辨率、5μm×5μm像元的大小的线阵传感器;畸变矫正方法的具体步骤如下:
一、畸变关系建立
构建标准图像与畸变图像坐标之间的映射关系,生成畸变矫正系数文件或畸变矫正表格文件;
二、畸变图像矫正
利用生成的畸变矫正系数文件或畸变矫正表格文件对输入的畸变图像进行矫正,最终获得矫正后的标准图像。
2.根据权利要求1所述的用于全视角三维测量仪光学系统的畸变矫正方法,其特征在于所述畸变关系建立的具体步骤如下:
1)根据用户提供的相机参数和实际标靶,利用理想成像投影关系计算控制点的标准图像坐标;
2)根据实际畸变后的靶标图像,获取控制点在畸变图像坐标系中的坐标;
3)如果采用函数表达法对畸变关系进行描述,则利用步骤1)、2)获得的控制点标准图像坐标和畸变图像坐标对映射函数中的参数进行拟合,并按照协议生成畸变矫正系数文件;
4)如果采用表格法对畸变关系进行描述,则对步骤1)、2)获得的控制点标准图像坐标和畸变图像坐标插值,并按照协议生成畸变矫正表格文件。
3.根据权利要求1所述的用于全视角三维测量仪光学系统的畸变矫正方法,其特征在于所述畸变图像矫正的具体步骤如下:
1)如果利用畸变矫正系数文件进行图像矫正,则利用文件中的系数计算矫正后标准图像各像元在畸变图像中对应的坐标;
2)如果利用畸变矫正表格文件进行图像矫正,则在文件中的表格中查找矫正后标准图像各像元在畸变图像中对应的坐标;
3)根据矫正后标准图像各像元在畸变图像中的坐标,在用户的输入畸变图像中插值查找这些坐标对应的图像灰度信息;
4)将查找到的灰度信息进行整合,对标准图像中像素坐标超出畸变图像范围的点进行灰度置零处理,生成畸变矫正后的图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611224243.9A CN106780391B (zh) | 2016-12-27 | 2016-12-27 | 一种用于全视角三维测量仪光学系统的畸变矫正算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611224243.9A CN106780391B (zh) | 2016-12-27 | 2016-12-27 | 一种用于全视角三维测量仪光学系统的畸变矫正算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106780391A CN106780391A (zh) | 2017-05-31 |
CN106780391B true CN106780391B (zh) | 2020-03-13 |
Family
ID=58926715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611224243.9A Active CN106780391B (zh) | 2016-12-27 | 2016-12-27 | 一种用于全视角三维测量仪光学系统的畸变矫正算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106780391B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109615586B (zh) * | 2018-05-07 | 2022-03-11 | 杭州新瀚光电科技有限公司 | 红外图像畸变矫正算法 |
CN109277695A (zh) * | 2018-08-01 | 2019-01-29 | 普聚智能系统(苏州)有限公司 | 一种适用于立体结构工件的激光同轴视觉系统 |
CN113327202A (zh) * | 2021-03-30 | 2021-08-31 | 苏州微清医疗器械有限公司 | 一种图像畸变的矫正方法及其应用 |
CN113465541B (zh) * | 2021-07-09 | 2022-02-22 | 北京运达华开科技有限公司 | 一种接触线磨耗测量方法及系统 |
CN113837960A (zh) * | 2021-09-22 | 2021-12-24 | 中国计量大学 | 一种基于深度学习的畸变矫正方法 |
CN114173103A (zh) * | 2021-12-08 | 2022-03-11 | 沂普光电(天津)有限公司 | 一种投影系统辅助装调装置与方法 |
CN115587952B (zh) * | 2022-12-13 | 2023-03-14 | 腾讯科技(深圳)有限公司 | 扩展现实设备的畸变系数标定方法、装置和存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1525384A (zh) * | 2003-06-20 | 2004-09-01 | 北京中星微电子有限公司 | 镜头成像畸变校正的方法 |
CN103426149A (zh) * | 2013-07-24 | 2013-12-04 | 玉振明 | 大视角图像畸变的校正处理方法 |
CN104182933A (zh) * | 2013-05-28 | 2014-12-03 | 东北大学 | 一种基于逆向除法模型的广角镜头图像畸变校正方法 |
CN104537616A (zh) * | 2014-12-20 | 2015-04-22 | 中国科学院西安光学精密机械研究所 | 鱼眼图像畸变的校正方法 |
CN105354809A (zh) * | 2015-12-05 | 2016-02-24 | 中国航空工业集团公司洛阳电光设备研究所 | 一种基于输出图像像素位置索引的预畸变方法和装置 |
-
2016
- 2016-12-27 CN CN201611224243.9A patent/CN106780391B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1525384A (zh) * | 2003-06-20 | 2004-09-01 | 北京中星微电子有限公司 | 镜头成像畸变校正的方法 |
CN104182933A (zh) * | 2013-05-28 | 2014-12-03 | 东北大学 | 一种基于逆向除法模型的广角镜头图像畸变校正方法 |
CN103426149A (zh) * | 2013-07-24 | 2013-12-04 | 玉振明 | 大视角图像畸变的校正处理方法 |
CN104537616A (zh) * | 2014-12-20 | 2015-04-22 | 中国科学院西安光学精密机械研究所 | 鱼眼图像畸变的校正方法 |
CN105354809A (zh) * | 2015-12-05 | 2016-02-24 | 中国航空工业集团公司洛阳电光设备研究所 | 一种基于输出图像像素位置索引的预畸变方法和装置 |
Non-Patent Citations (4)
Title |
---|
一款宽光谱鱼眼镜头的设计;李宏壮 等;《光子学报》;20121130;第41卷(第11期);第1312-1316页 * |
超广角大口径fq透镜设计;崔桂华;《应用光学》;19931231;第14卷(第3期);第13-15页 * |
超广角大工作面f-theta镜头的光学设计;付敏敏 等;《应用光学》;20111231;第32卷(第6期);第1083-1087页 * |
轻小型立体相机光学系统研制;孟庆宇 等;《红外与激光工程》;20160430;第45卷(第4期);第0418002-1至0418002-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106780391A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106780391B (zh) | 一种用于全视角三维测量仪光学系统的畸变矫正算法 | |
CN101858755B (zh) | 一种星敏感器的标定方法 | |
CN107358631B (zh) | 一种虑及三维畸变的双目视觉重建方法 | |
US8934721B2 (en) | Microscopic vision measurement method based on adaptive positioning of camera coordinate frame | |
CN109685744B (zh) | 一种扫描振镜精度校正方法 | |
CN102269587B (zh) | 基于可控光平面的水下三维重绘装置的重绘方法 | |
WO2020199439A1 (zh) | 基于单双目混合测量的三维点云计算方法 | |
CN106403838A (zh) | 一种手持式线结构光视觉三维扫描仪的现场标定方法 | |
US20190193648A1 (en) | Method and system of camera focus for advanced driver assistance system (adas) | |
CN102663727B (zh) | 基于cmm移动靶标的相机分区域标定方法 | |
CN111595302A (zh) | 一种双面阵ccd辅助三线阵ccd位姿光学测量及校准方法 | |
CN109242909B (zh) | 一种面向高精度二维尺寸测量的线阵相机标定算法 | |
CN110068313B (zh) | 一种基于投影变换的数字天顶仪定向方法 | |
CN105259645A (zh) | 一种大视场无畸变缩比物镜 | |
CN103134443B (zh) | 一种大口径大径厚比反射镜面形自准直检测装置及方法 | |
CN110146032A (zh) | 基于光场分布的合成孔径相机标定方法 | |
CN108765301B (zh) | 一种光学系统及其畸变校正方法及系统 | |
CN116105638A (zh) | 一种u型转台同轴度检测系统及检测方法 | |
CN116051651A (zh) | 一种变视线成像系统的等效多目视觉模型及其标定方法 | |
CN113781581B (zh) | 基于靶标松姿态约束的景深畸变模型标定方法 | |
CN110689582B (zh) | 全站仪摄像机标定方法 | |
CN203519000U (zh) | 标志点正交分光成像位姿传感器 | |
CN103528569A (zh) | 标志点正交分光成像位姿测试方法及传感器 | |
KR101482637B1 (ko) | 다중대역 비열화보상 광학시스템의 온도변화에 따른 시계편차 보정 방법 | |
Duan et al. | A method of camera calibration based on image processing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |