CN106596063A - 一种测量透镜畸变的方法及系统 - Google Patents
一种测量透镜畸变的方法及系统 Download PDFInfo
- Publication number
- CN106596063A CN106596063A CN201611147288.0A CN201611147288A CN106596063A CN 106596063 A CN106596063 A CN 106596063A CN 201611147288 A CN201611147288 A CN 201611147288A CN 106596063 A CN106596063 A CN 106596063A
- Authority
- CN
- China
- Prior art keywords
- point
- central
- distortion
- center
- coordinate
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0207—Details of measuring devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0242—Testing optical properties by measuring geometrical properties or aberrations
- G01M11/0257—Testing optical properties by measuring geometrical properties or aberrations by analyzing the image formed by the object to be tested
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0242—Testing optical properties by measuring geometrical properties or aberrations
- G01M11/0257—Testing optical properties by measuring geometrical properties or aberrations by analyzing the image formed by the object to be tested
- G01M11/0264—Testing optical properties by measuring geometrical properties or aberrations by analyzing the image formed by the object to be tested by using targets or reference patterns
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Image Processing (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种测量透镜畸变的方法及系统,该方法包括:设置具有K×N个点构成的点阵图案的测试卡,K与N均为自然数,K等于或不等于N;获取测试卡通过透镜产生畸变后的畸变图像;以畸变图像左上角的点为坐标原点,以原点向右为X轴正方向,以原点向下为Y轴正方向构建畸变图像的平面坐标系;通过扫描检索定位畸变图像的中心点的坐标以及所有非中心点,并确定中心点以及所有非中心点在平面坐标系中的坐标值;利用中心点以及所有非中心点的坐标值,根据畸变图像畸变量计算公式,计算畸变图像的畸变量,进而得到透镜的畸变量。本测量方法和系统能够快速、准确的定位测试卡的黑点,加快了图像处理速度的同时提高了透镜畸变测量的精度。
Description
本发明专利申请是申请日为2014年06月27日、申请号为201410302056.2、名称为“一种测量透镜畸变的方法及系统”的中国发明专利申请的分案申请。
技术领域
本发明涉及图像处理领域,具体涉及一种测量透镜畸变的方法及系统。
背景技术
随着科技水平的提高,消费类电子产品以飞快的步伐迈入人类生活,头戴显示类产品凭借绚丽的外形及先进的功能吸引着大众的眼球。头戴显示类产品将视觉技术完美的应用与产品,其使用高分辨率的LCD屏显示图像,经过透镜放大,以便用户能观看到清晰的图像;并采用3D技术处理图像,让用户感受3D带来的视觉冲击感。然而透镜将图像放大的同时,也使得图像产生了畸变,畸变属于成像的几何失真,它是由于焦平面上不同区域对影像的放大率不同而形成的画面扭曲变形现象,这种变形的程度从画面中心至画面边缘依次递增,主要在画面边缘反映得较明显。因此需要测量到这个畸变量,再通过图像处理技术将畸变的图像还原成正常的图像,现有技术中的测量透镜畸变的方法计算速度慢并且测量精度不高,还原为正常图像具有较大的误差,不能满足图像处理的需求。
发明内容
本发明提供了一种测量透镜畸变的方法及系统,以解决现有的透镜畸变测量方法计算速度慢且测量精度不高的问题。
为达到上述目的,本发明的技术方案是这样实现的:
本发明提供了一种测量透镜畸变的方法,该方法包括:设置具有K×N个点构成的点阵图案的测试卡,K与N均为自然数,K等于或不等于N;
获取测试卡通过透镜产生畸变后的畸变图像;
以畸变图像左上角的点为坐标原点,以原点向右为X轴正方向,以原点向下为Y轴正方向构建畸变图像的平面坐标系;
基于中心点与所有非中心点的物理相对位置不变的关系,采用坐标定位方式,通过扫描检索定位畸变图像的中心点以及所有非中心点,并确定中心点以及所有非中心点在平面坐标系中的坐标值;
利用中心点以及所有非中心点的坐标值,计算畸变图像的畸变量,进而得到透镜的畸变量。
其中,测试卡上中心点的半径比所有非中心点的半径大。
其中,通过扫描检索定位畸变图像的中心点以及所有非中心点包括:
设定一个以畸变图像中心的坐标值为中心的正方形检索区域;正方形检索区域的边长的一半大于等于中心点半径的2倍;
以正方形检索区域的左上角的一点为中心,设定一个边长的一半为中心点半径的正方形扫描区域,以扫描区域的中心为扫描起始点,扫描该扫描区域并计算该区域内所有像素点的平均值;
在检索区域内,按照从上到下、从左到右的顺序移动扫描区域的中心,直至扫描完检索区域,比较扫描区域每次获取到的像素点平均值,并将所有像素点平均值中的最小平均值或者最大平均值所在的扫描区域的中心确定为中心点的中心,进而确定中心点中心的坐标值,将中心点中心的坐标值定位为中心点的坐标值。
以此类似的方法,定位出所有非中心点的坐标值。
其中,以此类似的方法,定位出所有非中心点的坐标值包括:
以中心点的中心为基准,以特定的步长,向左移动正方形检索区域,定位出与中心点同一行的左边所有的非中心点,并确定非中心点的坐标值;
以中心点的中心为基准,以特定的步长,向右移动正方形检索区域,定位出中心点同一行的右边所有非中心点,并确定非中心点的坐标值;
以中心点的中心为基准,以特定的步长,向下移动正方形检索区域,定位出中心点所在行的下一行所有非中心点,并确定非中心点的坐标值;
以中心点的中心为基准,以特定的步长,向上移动正方形检索区域,定位出中心点所在行的上一行所有非中心点,并确定非中心点的坐标值;
在所述检索区域内,按照从上到下、从左到右的顺序移动所述扫描区域的中心,直至扫描完所述检索区域。
其中,该方法还包括:在缓存中定义一个长度能够至少包含所有点的二维数组,将中心点的坐标信息存储于二维数组的中心数值中;
依次将定位到的所有点按照畸变图像上的相对物理位置,存储于二维数组中,建立中心点以及所有非中心点与二维数组的索引关系。
其中,利用中心点以及非中心点的坐标值,计算畸变图像的畸变量包括:
根据二维数组与中心点以及所有非中心点的索引关系,找出分别位于在畸变图像上的四个点A1、A2、A3、A4,所述点A1、A2、A3、A4应尽可同时满足以下条件:
四个点A1、A2、A3、A4两、两处于同一基准线;
根据四个点A1、A2、A3、A4及其两两所在的同一基准线上所有非中心点与中心点的索引关系,找出所述四个点A1、A2、A3、A4两两所在的同一基准线上的中间点,其中,点A1与点A2所在的基准线的中间点记为B1,点A2与点A3所在的基准线的中间点记为B2,点A3与点A4所在的基准线的中间点记为B3,点A1与点A4所在的基准线的中间点记为B4,中间点B1、B2、B3、B4同时满足以下条件:
B1点、中心点、B3点这三个点位于同一基准线;
B2点、中心点、B4点这三个点位于同一基准线;
点A1、A2相对于点B1的索引差值相等并记为m,点A3、A4相对于点B3的索引差值相等并记为n,且m=n;
点A1、A4相对于点B4索引差值相等并记为u,点A2、A3相对于点B2的索引差值相等并记为v,且v=u。
其中,利用中心点以及非中心点的坐标值,计算畸变图像的畸变量包括:
根据如下公式计算出畸变图像垂直方向的畸变量:
Vertical=100%*(DisA1A4+DisA2A3)/(2*DisB1B3)
其中,DisA1A4表示点A1和A4的之间距离,DisA2A3表示点A2和A3之间的距离,DisB1B3表示点B1和B3之间的距离,Vertical表示畸变图像垂直方向的畸变量;
根据如下公式计算出畸变图像水平方向的畸变量:
Horizontal=100%*(DisA1A2+DisA3A4)/(2*DisB2B4)
其中,DisA1A2表示点A1和A2的距离,DisA3A4表示点A3和A4的距离,DisB2B表示点B2和B4的距离,Horizontal表示畸变图像水平方向的畸变量。
其中,获取测试卡通过透镜产生畸变后的畸变图像包括:
采用高像素的工业相机,相机中心、透镜中心以及测试卡中心重合,通过透镜拍摄测试卡,得到畸变图像。
其中,当中心点以及非中心点为黑色的点时,将所有像素点平均值中的最小平均值所在的扫描区域的中心确定为中心黑点的中心,并确定中心黑点中心的坐标值。
本发明还提供一种测量透镜畸变的系统,应用上述的方法测量透镜畸变,该系统包括:
测试卡,具有K×N个点构成的点阵图案;
成像装置,用于获取测试卡通过透镜产生畸变后的畸变图像;
图像处理装置,用于以畸变图像左上角的点为坐标原点,以原点向右为X轴正方向,以原点向下为Y轴正方向构建畸变图像的平面坐标系;基于中心点与所有非中心点的物理相对位置不变的关系,采用坐标定位方式,通过扫描检索定位畸变图像的中心点以及所有非中心点,并确定中心点以及所有非中心点在平面坐标系中的坐标值;以及利用中心点以及非中心点的坐标值,计算畸变图像的畸变量,进而得到透镜的畸变量。
本发明的这种测量透镜畸变的方法及系统,能够快速、准确的定位测试卡上的所有点,实际误差范围达到亚像素级别,达到快速的计算畸变图像的畸变量,真实反应出透镜的畸变的目的。
附图说明
图1是本发明一个实施例提供的一种测量透镜畸变的方法的流程图;
图2是本发明一个实施例提供的测试卡的示意图;
图3是本发明一个实施例提供的获取畸变图像的示意图;
图4是本发明一个实施例提供的测试卡经过透镜畸变后的畸变图像示意图;
图5是本发明一个实施例提供的对畸变图像的黑点进行定位的模型图;
图6是本发明一个实施例提供的定位出的所有黑点实物图;
图7是本发明一个实施例提供的黑点坐标缓存的示意图;
图8是本发明一个实施例提供的畸变图像畸变量计算的参照图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本发明的核心思想是利用一种新的畸变图像测试卡,该测试卡上设置有黑点这些黑点组成点阵图案,并基于该测试卡进行坐标定位,快速识别出畸变图像的所有黑点的坐标,并根据所有黑点的坐标计算出畸变图像的畸变量。本发明的这种测量方法能够快速的找到测试卡上的黑点,加快畸变量测量的处理速度,提高测量精度。
图1是本发明一个实施例提供的一种测量透镜畸变的方法的流程图;参见图1,该方法包括:
步骤S110,设置具有K×N个点构成的点阵图案的测试卡,K与N均为自然数,K等于或不等于N;
步骤S120,获取测试卡通过透镜产生畸变后的畸变图像;
步骤S130,以畸变图像左上角的点为坐标原点,以原点向右为X轴正方向,以原点向下为Y轴正方向构建畸变图像的平面坐标系;
步骤S140,通过扫描检索定位畸变图像的中心点以及所有非中心点,并确定中心点以及所有非中心点在平面坐标系中的坐标值;
步骤S150,利用中心点以及非中心点的坐标值,计算畸变图像的畸变量,进而得到透镜的畸变量。
通过构建畸变图像平面坐标系,快速识别畸变图像中所有的点,并计算畸变图像上所有点的坐标值,该测量方法计算量少,进而提高了畸变图像测量时的速度,同时,采用K×N个点构成的点阵图案的测试卡保证了测量精度的提高。
图2是本发明一个实施例提供的测试卡的示意图;参见图2,显示了本实施例中所使用的测试卡的局部:一个21行和21列的黑色的点构成的方阵,1表示方阵的中心黑点;中心黑点的半径比其他所有非中心黑点的半径大,以便于在拍摄畸变图像时将其与相机中心对齐,另外,中心黑点半径较大也便于在利用图像处理技术进行计算处理时识别其为中心点。可以理解,图2只是本实施例截取的透镜拍摄到的测试卡的畸变图像中间的部分,由于透镜导致测试卡的畸变,畸变图像边缘部分会呈现近似圆形的形状,中间部分近似矩阵排列。
在头戴显示类产品等其他类似智能显示产品使用高分辨率的LCD屏显示图像时,需要经过透镜放大,以便用户能观看到清晰的图像,但透镜将图像放大的同时,会使图像产生畸变。本发明的方法是对畸变图像畸变量的测量,在进行测量前需要获取畸变图像。图3是本发明一个实施例提供的获取畸变图像的示意图,参见图3,高像素(例如1千万像素)的工业相机31通过透镜32拍摄置于头戴显示类产品LCD屏33上的测试卡,将测试卡图片的中心与所述相机31的中心以及透镜32的中心重合,则此时拍摄到的畸变图像就能准确表示透镜造成图像的畸变量。
图4是本发明一个实施例提供的测试卡经过透镜畸变后的畸变图像示意图,图4能最大反映出计算透镜畸变的合适区域,边缘上四个角的点尽可能的靠近图像边缘,且两两处于边缘上最外围且最完整的同一基准线,一般为F×F的点阵区域,F为奇数;参见图4,测试卡经过透镜放大的同时也发生了枕形畸变,畸变使得测试卡上的黑点位置发生偏移,并且每个黑点偏移的位置不一样,畸变图像上四周边缘位置的黑点的畸变程度最严重。本发明采用坐标定位方式定位该畸变图像上所有黑点的坐标值,在测试卡经过透镜发生了畸变后的畸变图像上,中心点以及其他非中心点的物理相对位置是不变的,具体的坐标值可能发生了变化,利用坐标定位技术定位出所有黑点的坐标,由于黑点的相对物理位置不变,而实际的坐标值发生了改变,利用改变的所有黑点的坐标值以及畸变量计算公式得到该畸变图像的畸变量值,并进而得到透镜的畸变量,从而调整畸变图像以减小还原为正常图像时的误差。
下面结合图5和图6具体说明对畸变图像的所有黑点进行定位:图5是本发明一个实施例提供的对畸变图像的黑点进行定位的模型图;图6是本发明一个实施例提供的定位出的所有黑点实物图;
扫描检索定位畸变图像的中心点以及所有非中心点的方法具体为:设定一个以畸变图像中心的坐标值为中心的正方形检索区域;正方形检索区域的边长的一半大于等于中心点半径的2倍;以正方形检索区域的左上角的一点为中心,设定一个边长的一半为中心点半径的正方形扫描区域,以扫描区域的中心为扫描起始点,扫描该扫描区域内的图像并计算该区域内图像所有像素点的平均值;在检索区域内,按照顺序移动扫描区域的中心,直至扫描完检索区域,比较扫描区域每次获取到的像素点平均值,并将所有像素点平均值中的最小平均值或者最大平均值所在的扫描区域的中心确定为中心点的中心,并确定中心点中心的坐标值,将中心点中心的坐标值定位为中心点的坐标值。
以此类似的方法,定位出所有非中心点的坐标值。
参见图5,在本实施例中,测试卡上中心点和非中心为黑色的点,在图5所示的中心黑点所在的黑色区域中,以黑点的半径52的2倍或者比黑点半径(即正方形边长的一半)2倍更大的半径51定义一个正方形的检索区域,以该检索区域左上角的一点为中心,定义一个以黑点半径(即正方形边长的一半)为半径53的方形扫描区域,在该半径(即正方形边长的一半)较大的检索区域内,按照从上到下、从左到右的顺序依次移动该扫描区域的中心进行扫描;每扫描一次,记录一个该扫描区域内所有像素点的平均值;直至将该检索区域扫描完毕,比较每次扫描区域获取到的平均像素值中最小平均值所在的扫描区域,以该扫描区域的中心确定为中心黑点的中心,完成中心黑点的定位,计算该中心黑点坐标值并保存。可以理解,在应用本发明的方法时,测试卡上的点的颜色不限于黑色这一种,因而在对测试卡上的点进行定位时,应该根据具体点颜色的像素值选取平均值最小或者最大的扫描区域作为中心点的中心点,在本实施例中,测试卡上点的颜色为黑色,那么扫描区域包含的黑色部分越多,像素值较小,纯黑的像素值为0。在本发明的其他实施例中,测试卡的底色可能为其他色,中心点以及非中心点的颜色例如白色这种情况下,在利用扫描区域进行定位时,将扫描区域像素点平均值中最大平均值(白色的像素值为255)所在的扫描区域的中心定位为中心点的中心。
与定位中心黑点的方法类似,定位出所有非中心黑点,具体如下:以中心黑点的中心为基准,以特定的步长,向左移动正方形检索区域,定位出与中心黑点的同一行的左边所有非中心黑点,并确定非中心黑点的坐标值;然后向右移动正方形检索区域,定位出中心黑点同一行的右边所有非中心黑点,并确定非中心黑点的坐标值;以中心黑点的中心为基准,以特定的步长,向下移动正方形检索区域,定位出中心黑点所在行的下一行所有非中心黑点,并确定非中心黑点的坐标值;以中心黑点的中心为基准,以特定的步长,向上移动正方形检索区域,定位出中心黑点所在行的上一行所有非中心黑点,并确定非中心黑点的坐标值。
在定位非中心黑点时,选取的正方形的检索区域的边长的一半为非中心黑点的半径2倍或者更大的值。在移动检索区域时的特定步长是指畸变图像上每两个黑点中心之间的距离,它是一个预先设定的经验值,每次移动这个距离可以避免不必要的扫描,提高扫描区域扫描的速度和效率。
图6是本发明一个实施例提供的定位出的所有黑点实物图;参见图6,定位出中心黑点和所有非中心黑点后用圆圈进行标注。通过上述步骤确定中心黑点以及非中心其他黑点的坐标值并保存。
在对坐标值进行保存时,本发明一个实施例采用二维数组的方式进行存储,图7是本发明一个实施例提供的黑点坐标缓存的示意图;参见图7示例性的选取了畸变图像上的几个点,将选取的这些点的坐标值按照在畸变图像上的相对物理位置进行保存,先在缓存中定义一个长度能够至少包含所有点的二维数组,(根据应用场景的不同,畸变图像上黑点的数量不同,设置不同长度的二维数组)例如本实施例中,定义了一个长度为41的二维数组array[41][41],定义了二维数组后,先将中心黑点的坐标进行存储,在存储时找到二维数组的最中间位置array[21][21],将中心黑点的坐标放入数组array[21][21],将畸变图像中心黑点同一行的左边那个黑点的坐标放入array[20][21]中,以此类推,将畸变图像上其他黑点的坐标值放入二维数组中,建立畸变图像黑点与二维数组的索引关系。
图8是本发明一个实施例提供的畸变图像畸变量计算的参照图,参见图8,在得到所有黑点的坐标值后,观察畸变图像,选取畸变图像上畸变最能表现畸变程度的几个关键点,参见图8,快速查找缓存中二维数组,先找到能最大程度反映计算透镜畸变的四个点A1、A2、A3、A4,这四个点应尽可能的靠近畸变图像边缘,即A1、A2、A3、A4这四个点所限定的区域能最大反映出计算透镜畸变的合适区域,A1、A2、A3、A4这四个点两两处于边缘上最外围且最完整的同一基准线上并定义出选取的该基准线的末端),并同时满足以下条件:
A1、A2、A3、A4这四个点两、两处于同一基准线;点A1与点A2连线与中心点所在的竖直基准线相交的点记为B1;点A2与点A3连线与中心点所在的水平基准线相交的点记为B2;同理,确定出B3和B4点。点A1、A2相对于B1点的索引差值相等记为m,点A3、A4相对于B3点的索引差值相等记为n,并且m=n=9,B1、O(中心黑点)、B3这三点位于同一基准线;同理,A1、A4相对B4点索引差值相等记为u,A2、A3相对于B2点的索引差值相等记为v,并且v=u=9,且B2、O(中心黑点)、B4这三点位于同一基准线,索引差值为9表示中心黑点距离B1点相隔9个黑点。实际中,B1、B2、B3、B4可以看成是四个点A1、A2、A3、A4两、两所在的同一基准线上的中间点(或称中点),即点A1与点A2所在的基准线的中心点为B1,点A2与点A3所在的基准线的中间点为B2,点A3与点A4所在的基准线的中间点为B3,点A1与点A4所在的基准线的中间点为B4。
得到上面关键点A1-A4和B1-B4的坐标值后,根据畸变量计算公式进行计算:
根据如下公式计算出畸变图像水平方向的畸变量:
Horizontal=100%*(DisA1A2+DisA3A4)/(2*DisB2B4)
其中,DisA1A2表示点A1和A2的距离,DisA3A4表示点A3和A4的距离,DisB2B4表示点B2和B4的距离,Horizontal表示畸变图像水平方向的畸变量。
同理,根据如下公式计算出畸变图像垂直方向的畸变量:
Vertical=100%*(DisA1A4+DisA2A3)/(2*DisB1B3)
其中,DisA1A4表示点A1和A4的距离,DisA2A3表示点A2和A3的距离,DisB1B3表示点B1和B3的距离,Vertical表示畸变图像垂直方向的畸变量。
计算得到畸变图像的水平方向的畸变量和垂直方向的畸变量后,即可得到该畸变图像的畸变量,进而得到透镜的畸变量,对该畸变图像进行调整,以更好的还原图像,减小误差。
本发明一个实施例还提供了一种测量透镜畸变的系统,应用图1所示的方法测量透镜畸变,该系统包括:
测试卡,具有K×N个点构成的点阵图案;
成像装置,用于获取测试卡通过透镜产生畸变后的畸变图像;
图像处理装置,用于以畸变图像左上角的点为坐标原点,以原点向右为X轴正方向,以原点向下为Y轴正方向构建畸变图像的平面坐标系;通过扫描检索定位畸变图像的中心点以及所有非中心点,并确定中心点以及所有非中心点在平面坐标系中的坐标值;以及利用中心点以及非中心点的坐标值,计算畸变图像的畸变量,进而得到透镜的畸变量。
综上,本发明的这种测量透镜畸变的方法及系统,基于畸变图像定义了新的平面坐标系,利用坐标定位,快速、准确的找到测试卡上的点,加快图像处理的速度和效率;而且误差范围达到亚像素级别,提高了测量精度;另外利用缓存数组存储图像上点的坐标信息,方便快速查找图像点的坐标并计算畸变图像的畸变量,较真实的反映出透镜的畸变程度。
以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
Claims (10)
1.一种测量透镜畸变的方法,其特征在于,该方法包括:
设置具有K×N个点构成的点阵图案的测试卡,K与N均为自然数,K等于或不等于N;
获取所述测试卡通过所述透镜产生畸变后的畸变图像;
以所述畸变图像左上角的点为坐标原点,以原点向右为X轴正方向,以原点向下为Y轴正方向构建所述畸变图像的平面坐标系;
基于中心点与所有非中心点的物理相对位置不变的关系,采用坐标定位方式,通过扫描检索定位所述畸变图像的中心点以及所有非中心点,并确定中心点以及所有非中心点在所述平面坐标系中的坐标值;
利用所述中心点以及所有非中心点的坐标值,计算所述畸变图像的畸变量,进而得到所述透镜的畸变量。
2.如权利要求1所述的方法,其特征在于,
所述测试卡上中心点的半径比所有非中心点的半径大。
3.如权利要求1所述的方法,其特征在于,所述通过扫描检索定位所述畸变图像的中心点以及所有非中心点包括:
设定一个以所述畸变图像中心的坐标值为中心的正方形检索区域;所述正方形检索区域的边长的一半大于等于所述中心点半径的2倍;
以所述正方形检索区域的左上角的一点为中心,设定一个边长的一半为中心点半径的正方形扫描区域,以所述扫描区域的中心为扫描起始点,扫描所述扫描区域并计算该区域内所有像素点的平均值;
在所述检索区域内,按照从上到下、从左到右的顺序移动所述扫描区域的中心,直至扫描完所述检索区域,比较所述扫描区域每次获取到的像素点平均值,并将所有像素点平均值中的最小平均值或者最大平均值所在的扫描区域的中心确定为所述中心点的中心,进而确定所述中心点中心的坐标值,将所述中心点中心的坐标值定位为所述中心点的坐标值;
以此类似的方法,定位出所有非中心点的坐标值。
4.如权利要求3所述的方法,其特征在于,
所述以此类似的方法,定位出所有非中心点的坐标值包括:
以所述中心点的中心为基准,以特定的步长,向左移动所述正方形检索区域,定位出与所述中心点同一行的左边所有的非中心点,并确定所述非中心点的坐标值;
以所述中心点的中心为基准,以特定的步长,向右移动所述正方形检索区域,定位出与所述中心点同一行的右边所有非中心点,并确定所述非中心点的坐标值;
以所述中心点的中心为基准,以特定的步长,向下移动所述正方形检索区域,定位出所述中心点所在行的下一行所有非中心点,并确定所述非中心点的坐标值;
以所述中心点的中心为基准,以特定的步长,向上移动所述正方形检索区域,定位出所述中心点所在行的上一行所有非中心点,并确定所述非中心点的坐标值;
在所述检索区域内,按照从上到下、从左到右的顺序移动所述扫描区域的中心,直至扫描完所述检索区域。
5.如权利要求4所述的方法,其特征在于,所述方法还包括:
在缓存中定义一个长度能够至少包含所有点的二维数组,将所述中心点的坐标信息存储于所述二维数组的中心数组中;
依次将定位到的所有点按照所述畸变图像上的相对物理位置,存储于所述二维数组中,建立所述中心点以及所有非中心点与所述二维数组的索引关系。
6.如权利要求5所述的方法,其特征在于,
所述利用所述中心点以及所有非中心点的坐标值,计算所述畸变图像的畸变量包括:
根据所述二维数组与所述中心点以及所有非中心点的索引关系,找出分别位于所述畸变图像边缘上的四个点A1、A2、A3、A4,所述四个点A1、A2、A3、A4,同时满足以下条件:四个点A1、A2、A3、A4两、两处于所述畸变图像边缘上的同一基准线;
根据所述四个点A1、A2、A3、A4及其两两所在的同一基准线上所有非中心点与所述中心点的索引关系,找出所述四个点A1、A2、A3、A4两两所在的同一基准线上的中间点,其中,点A1与点A2所在的基准线的中间点记为B1,点A2与点A3所在的基准线的中间点记为B2,点A3与点A4所在的基准线的中间点记为B3,点A1与点A4所在的基准线的中间点记为B4,中间点B1、B2、B3、B4同时满足以下条件:
B1点、中心点、B3点这三个点位于同一基准线;
B2点、中心点、B4点这三个点位于同一基准线;
点A1、A2相对于点B1的索引差值相等并记为m,点A3、A4相对于点B3的索引差值相等并记为n,且m=n;
点A1、A4相对于点B4索引差值相等并记为u,点A2、A3相对于点B2的索引差值相等并记为v,且v=u。
7.如权利要求6所述的方法,其特征在于,所述利用所述中心点以及所有非中心点的坐标值,计算所述畸变图像的畸变量还包括:
根据如下公式计算所述畸变图像垂直方向的畸变量:
Vertical=100%*(DisA1A4+DisA2A3)/(2*DisB1B3)
其中,DisA1A4表示点A1和A4之间的距离,DisA2A3表示点A2和A3之间的距离,DisB1B3表示点B1和B3之间的距离,Vertical表示畸变图像垂直方向的畸变量;
根据如下公式计算所述畸变图像水平方向的畸变量:
Horizontal=100%*(DisA1A2+DisA3A4)/(2*DisB2B4)
其中,DisA1A2表示点A1和A2的距离,DisA3A4表示点A3和A4的距离,DisB2B4表示点B2和B4的距离,Horizontal表示畸变图像水平方向的畸变量。
8.如权利要求1所述的方法,其特征在于,所述获取所述测试卡通过所述透镜产生畸变后的畸变图像包括:
选取工业相机;
使得所述工业相机中心、透镜中心以及所述测试卡中心重合;
所述工业相机通过透镜拍摄所述测试卡,得到所述畸变图像。
9.如权利要求3所述的方法,其特征在于,
当所述中心点以及所述非中心点为黑色的点时,将所述所有像素点平均值中的最小平均值所在的扫描区域的中心确定为所述中心黑点的中心,并确定所述中心黑点中心的坐标值。
10.一种测量透镜畸变的系统,其特征在于,应用权利要求1至9中任一项所述的方法测量透镜畸变,所述系统包括:
测试卡,具有K×N个点构成的点阵图案;
成像装置,用于获取所述测试卡通过所述透镜产生畸变后的畸变图像;
图像处理装置,用于以所述畸变图像左上角的点为坐标原点,以原点向右为X轴正方向,以原点向下为Y轴正方向构建所述畸变图像的平面坐标系;基于中心点与所有非中心点的物理相对位置不变的关系,采用坐标定位方式,通过扫描检索定位所述畸变图像的中心点以及所有非中心点,并确定中心点以及所有非中心点在所述平面坐标系中的坐标值;以及利用所述中心点以及非中心点的坐标值,计算所述畸变图像的畸变量,进而得到所述透镜的畸变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611147288.0A CN106596063B (zh) | 2014-06-27 | 2014-06-27 | 一种测量透镜畸变的方法及系统 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410302056.2A CN104048815B (zh) | 2014-06-27 | 2014-06-27 | 一种测量透镜畸变的方法及系统 |
CN201611147288.0A CN106596063B (zh) | 2014-06-27 | 2014-06-27 | 一种测量透镜畸变的方法及系统 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410302056.2A Division CN104048815B (zh) | 2014-06-27 | 2014-06-27 | 一种测量透镜畸变的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106596063A true CN106596063A (zh) | 2017-04-26 |
CN106596063B CN106596063B (zh) | 2019-05-24 |
Family
ID=51501942
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410302056.2A Active CN104048815B (zh) | 2014-06-27 | 2014-06-27 | 一种测量透镜畸变的方法及系统 |
CN201611147288.0A Active CN106596063B (zh) | 2014-06-27 | 2014-06-27 | 一种测量透镜畸变的方法及系统 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410302056.2A Active CN104048815B (zh) | 2014-06-27 | 2014-06-27 | 一种测量透镜畸变的方法及系统 |
Country Status (4)
Country | Link |
---|---|
US (2) | US9810602B2 (zh) |
JP (1) | JP6166852B1 (zh) |
CN (2) | CN104048815B (zh) |
WO (1) | WO2015197019A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109040724A (zh) * | 2018-08-03 | 2018-12-18 | 信利光电股份有限公司 | 一种结构光投影器的光斑畸变检测方法、装置及可读存储介质 |
CN111609995A (zh) * | 2019-02-26 | 2020-09-01 | 弗提图德萨沃有限公司 | 一种光学模组装调测试方法及装置 |
CN112304573A (zh) * | 2020-09-21 | 2021-02-02 | 武汉高德智感科技有限公司 | 一种同时测量镜头畸变和mtf指标的方法和系统 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104048815B (zh) * | 2014-06-27 | 2017-03-22 | 歌尔科技有限公司 | 一种测量透镜畸变的方法及系统 |
CN106153302B (zh) * | 2015-03-24 | 2019-03-12 | 北京威斯顿亚太光电仪器有限公司 | 一种用于硬管内窥镜成像畸变的测量方法 |
CN106815869B (zh) * | 2016-10-28 | 2020-06-19 | 北京鑫洋泉电子科技有限公司 | 鱼眼相机的光心确定方法及装置 |
CN106815823B (zh) * | 2017-02-22 | 2020-02-07 | 广东工业大学 | 一种透镜畸变标定校正方法及其装置 |
CN108009981B (zh) * | 2017-09-26 | 2021-06-01 | 深圳市易成自动驾驶技术有限公司 | 畸变参数的寻优方法、装置及计算机可读存储介质 |
US10572982B2 (en) * | 2017-10-04 | 2020-02-25 | Intel Corporation | Method and system of image distortion correction for images captured by using a wide-angle lens |
CN110555879B (zh) * | 2018-05-31 | 2023-09-08 | 京东方科技集团股份有限公司 | 一种空间定位方法、其装置、其系统及计算机可读介质 |
CN109191374B (zh) * | 2018-10-10 | 2020-05-08 | 京东方科技集团股份有限公司 | 一种畸变参数测量方法、装置及系统 |
CN111579220B (zh) * | 2020-05-29 | 2023-02-10 | 江苏迪盛智能科技有限公司 | 一种分辨率板 |
US11734789B2 (en) | 2020-06-02 | 2023-08-22 | Immersive Tech, Inc. | Systems and methods for image distortion correction |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1508527A (zh) * | 2002-12-17 | 2004-06-30 | 北京航空航天大学 | 一种透镜畸变参数的标定方法 |
EP1662440A1 (en) * | 2004-11-30 | 2006-05-31 | IEE INTERNATIONAL ELECTRONICS & ENGINEERING S.A. | Method for determining the position of an object from a digital image |
CN1996389A (zh) * | 2007-01-09 | 2007-07-11 | 北京航空航天大学 | 基于共线特征点的摄像机畸变快速校正方法 |
CN101641951A (zh) * | 2006-09-15 | 2010-02-03 | 德萨拉技术匈牙利公司 | 改进影像质量的成像系统及相关方法 |
CN101673397A (zh) * | 2009-09-30 | 2010-03-17 | 青岛大学 | 一种基于lcd的数码相机非线性标定方法 |
CN102119326A (zh) * | 2008-08-13 | 2011-07-06 | 皇家飞利浦电子股份有限公司 | 测量和校正多斑扫描设备中的透镜畸变 |
CA2596125C (en) * | 2006-08-09 | 2012-07-10 | Research In Motion Limited | Method, device and system for evaluating a lens for an electronic device |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE69526635T2 (de) * | 1994-12-29 | 2002-12-05 | Koninklijke Philips Electronics N.V., Eindhoven | Bilderzeugungsgerät und Verfahren zur Verbesserung geometrischer optischer Bildverzerrungen |
JPH08201021A (ja) * | 1995-01-23 | 1996-08-09 | Mazda Motor Corp | キャリブレーション方法 |
JP2806307B2 (ja) * | 1995-06-29 | 1998-09-30 | 日本電気株式会社 | レンズディストーションの測定用レチクル及びその測定方法 |
US5978085A (en) * | 1997-03-07 | 1999-11-02 | Litel Instruments | Apparatus method of measurement and method of data analysis for correction of optical system |
US6002525A (en) * | 1998-07-06 | 1999-12-14 | Intel Corporation | Correcting lens distortion |
JP4535412B2 (ja) * | 2000-09-28 | 2010-09-01 | 東海光学株式会社 | レンズの性能評価方法 |
JP3732794B2 (ja) * | 2002-03-20 | 2006-01-11 | 株式会社東芝 | 寸法検査方法及びその装置並びにマスクの製造方法 |
TW565735B (en) | 2003-04-18 | 2003-12-11 | Guo-Jen Jan | Method for determining the optical parameters of a camera |
JP2004340753A (ja) * | 2003-05-15 | 2004-12-02 | Topcon Corp | キャリブレーションチャート画像表示装置 |
JP4179142B2 (ja) * | 2003-11-20 | 2008-11-12 | 株式会社デンソー | 車両用画像処理装置 |
JP2006071395A (ja) * | 2004-09-01 | 2006-03-16 | Nikon Corp | 較正方法及び位置合わせ方法 |
EP3159738B1 (en) * | 2005-03-25 | 2018-12-12 | Nikon Corporation | Method of determining distortion of a projection optical system |
JP4795017B2 (ja) * | 2005-12-28 | 2011-10-19 | 株式会社ニデック | 眼鏡レンズ評価装置 |
JP3998701B1 (ja) * | 2006-12-28 | 2007-10-31 | 健治 吉田 | ドットパターンが設けられたカード |
WO2009000906A1 (en) * | 2007-06-26 | 2008-12-31 | Dublin City University | A method for high precision lens distortion calibration and removal |
KR20090130603A (ko) * | 2008-06-16 | 2009-12-24 | 삼성전기주식회사 | 광각 렌즈 모듈의 평가 방법 및 이에 사용되는 평가 차트 |
KR100988872B1 (ko) * | 2009-07-08 | 2010-10-20 | 주식회사 나노포토닉스 | 회전 대칭형의 광각 렌즈를 이용하여 복합 영상을 얻는 방법과 그 영상 시스템 및 하드웨어적으로 영상처리를 하는 이미지 센서 |
JP2011033570A (ja) * | 2009-08-05 | 2011-02-17 | Micronics Japan Co Ltd | 光学レンズの歪曲収差の評価方法 |
CN201476957U (zh) * | 2009-08-19 | 2010-05-19 | 茂莱(南京)仪器有限公司 | 有限共轭变焦镜头畸变测试装置 |
CN102564731A (zh) * | 2010-12-16 | 2012-07-11 | 中国科学院西安光学精密机械研究所 | 一种透镜焦距及波前畸变测量装置 |
CN103292981A (zh) * | 2013-05-22 | 2013-09-11 | 中国科学院上海光学精密机械研究所 | 光学镜头畸变的测量装置和校正方法 |
CN104048815B (zh) | 2014-06-27 | 2017-03-22 | 歌尔科技有限公司 | 一种测量透镜畸变的方法及系统 |
-
2014
- 2014-06-27 CN CN201410302056.2A patent/CN104048815B/zh active Active
- 2014-06-27 CN CN201611147288.0A patent/CN106596063B/zh active Active
-
2015
- 2015-06-26 WO PCT/CN2015/082496 patent/WO2015197019A1/zh active Application Filing
- 2015-06-26 US US15/314,926 patent/US9810602B2/en active Active
- 2015-06-26 JP JP2016575034A patent/JP6166852B1/ja active Active
-
2017
- 2017-10-09 US US15/728,214 patent/US10151664B2/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1508527A (zh) * | 2002-12-17 | 2004-06-30 | 北京航空航天大学 | 一种透镜畸变参数的标定方法 |
EP1662440A1 (en) * | 2004-11-30 | 2006-05-31 | IEE INTERNATIONAL ELECTRONICS & ENGINEERING S.A. | Method for determining the position of an object from a digital image |
CA2596125C (en) * | 2006-08-09 | 2012-07-10 | Research In Motion Limited | Method, device and system for evaluating a lens for an electronic device |
CN101641951A (zh) * | 2006-09-15 | 2010-02-03 | 德萨拉技术匈牙利公司 | 改进影像质量的成像系统及相关方法 |
CN1996389A (zh) * | 2007-01-09 | 2007-07-11 | 北京航空航天大学 | 基于共线特征点的摄像机畸变快速校正方法 |
CN102119326A (zh) * | 2008-08-13 | 2011-07-06 | 皇家飞利浦电子股份有限公司 | 测量和校正多斑扫描设备中的透镜畸变 |
CN101673397A (zh) * | 2009-09-30 | 2010-03-17 | 青岛大学 | 一种基于lcd的数码相机非线性标定方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109040724A (zh) * | 2018-08-03 | 2018-12-18 | 信利光电股份有限公司 | 一种结构光投影器的光斑畸变检测方法、装置及可读存储介质 |
CN111609995A (zh) * | 2019-02-26 | 2020-09-01 | 弗提图德萨沃有限公司 | 一种光学模组装调测试方法及装置 |
CN112304573A (zh) * | 2020-09-21 | 2021-02-02 | 武汉高德智感科技有限公司 | 一种同时测量镜头畸变和mtf指标的方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
US10151664B2 (en) | 2018-12-11 |
US20180031442A1 (en) | 2018-02-01 |
JP2017524920A (ja) | 2017-08-31 |
CN104048815B (zh) | 2017-03-22 |
JP6166852B1 (ja) | 2017-07-19 |
WO2015197019A1 (zh) | 2015-12-30 |
US9810602B2 (en) | 2017-11-07 |
CN104048815A (zh) | 2014-09-17 |
US20170199099A1 (en) | 2017-07-13 |
CN106596063B (zh) | 2019-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104048815B (zh) | 一种测量透镜畸变的方法及系统 | |
CN108426585B (zh) | 一种光场相机的几何标定方法 | |
CN106091984B (zh) | 一种基于线激光的三维点云数据获取方法 | |
CN110542390B (zh) | 使用结构光的3d对象扫描方法 | |
CN107248159A (zh) | 一种基于双目视觉的金属工件缺陷检测方法 | |
CN106981081A (zh) | 一种基于深度信息提取的墙面平整度检测方法 | |
CN109089025A (zh) | 一种基于光场成像技术的影像仪数字聚焦方法 | |
CN110379013B (zh) | 一种基于多角度激光线扫描的三维重构系统 | |
CN104867113B (zh) | 图像透视畸变校正的方法及系统 | |
CN107154050A (zh) | 一种基于机器视觉的石材几何参数的自动获取方法 | |
CN109859272A (zh) | 一种自动对焦双目摄像头标定方法及装置 | |
CN101853524A (zh) | 使用图像序列生成玉米果穗全景图的方法 | |
CN108921864A (zh) | 一种光条中心提取方法及装置 | |
JP2013113600A (ja) | ステレオ3次元計測装置 | |
CN107941153B (zh) | 一种激光测距优化标定的视觉系统 | |
CN103345761A (zh) | 用于线阵高塔倾斜扫描成像数据校正及拼接的方法和系统 | |
CN107392849A (zh) | 基于图像细分的靶标识别与定位方法 | |
CN107358628A (zh) | 基于靶标的线阵图像处理方法 | |
JP6009206B2 (ja) | 3次元計測装置 | |
CN109360144B (zh) | 一种基于手机平台的图像实时校正改进方法 | |
CN107346040A (zh) | 裸眼3d显示设备的光栅参数的确定方法、装置及电子设备 | |
CN107770520A (zh) | 一种线阵相机辅助对焦方法及装置 | |
CN105115443B (zh) | 基于视觉测量技术的水平全视角高精度三维测量方法 | |
CN104915919A (zh) | 图像处理装置以及图像处理方法 | |
CN106973199A (zh) | 利用对焦距离扫描提高深度精确度的多光圈相机系统 |
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 | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201014 Address after: 261031 north of Yuqing street, east of Dongming Road, high tech Zone, Weifang City, Shandong Province (Room 502, Geer electronic office building) Patentee after: GoerTek Optical Technology Co.,Ltd. Address before: 266104 Laoshan Qingdao District North House Street investment service center room, Room 308, Shandong Patentee before: GOERTEK TECHNOLOGY Co.,Ltd. |