CN106651959A - 一种光场相机微透镜阵列几何参数的标定方法 - Google Patents
一种光场相机微透镜阵列几何参数的标定方法 Download PDFInfo
- Publication number
- CN106651959A CN106651959A CN201611014733.6A CN201611014733A CN106651959A CN 106651959 A CN106651959 A CN 106651959A CN 201611014733 A CN201611014733 A CN 201611014733A CN 106651959 A CN106651959 A CN 106651959A
- Authority
- CN
- China
- Prior art keywords
- image
- field camera
- light
- lenticule
- lenticular
- 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
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种光场相机微透镜阵列几何参数的标定方法,包括以下步骤:步骤一:将匀光板装在光场相机镜头的前端,对白光源进行拍照,得到光场相机的白图像;步骤二:对光场相机的白图像进行去马赛克处理,将其转换成真彩图像;步骤三:提取真彩图像中所有微透镜的边缘轮廓,保存边缘点的坐标;步骤四:根据对边缘点的坐标,得到所有微透镜的圆心坐标Ci(x,y)和相对应的半径Ri,并统计微透镜个数N。本发明提供的方法自动化程度高,直接对原始图像进行操作,无需预设模板,可为光场相机内外参数标定及相机光学模型提供准确的微透镜阵列数据,及后续的图像处理提供不同类型的微透镜图像。
Description
技术领域
本发明涉及一种相机几何参数的标定方法,属于数字图像处理领域。
背景技术
光场相机区别于传统相机,它是在主透镜和图像传感器之间安置一个微透镜阵列,传感器可以同时记录光线的强度信息和方向信息,实现对四维光场数据的采集。光场相机具有先拍照后对焦的特点,可以有效的捕捉更大范围的场景深度,展示真实场景的三维结构,光场成像技术逐渐应用于航空航天、三维重建、安全监视等诸多领域。
为了获取拍摄实景的三维信息,实现准确的追迹光线,必须对相机的内外参数进行标定。确定像素与微透镜的从属关系,需要对微透镜阵列的参数进行精确标定。但由于微透镜的种类、安装误差以及白图像的漫反射不一致性,使得微透镜中心与灰度峰值点对应关系存在偏差。
目前标定微透镜的几何参数的方法主要有两种:一种是灰度重心法,设定一个阀值T,将图像二值化,然后求每个微透镜子图像的灰度重心坐标,作为微透镜中心。此方法对采集图像时的光源均匀度与阀值的选取有很大要求,采光程度不同,对微透镜图像的灰度和二值化后的图像边缘范围都有很大影响,在应用时准确度不高,而且不能区分不同种类的透镜。另一种是模板匹配法,此法建立六边形网格模板与通过边缘检验算子检验出的边缘进行匹配,最佳匹配出的六边形网格的中心,即为微透镜中心。此法预先设定的模板都是固定参数,若用于其他规格的微透镜阵列,误差会大大增加。
发明内容
针对上述现有技术的不足,本发明提供一种直接对聚焦光场相机原始白图像进行操作,无需预设模板,能准确获取微透镜阵列中微透镜个数、中心位置以及半径大小的微透镜阵列的标定方法。
一种光场相机微透镜阵列几何参数的标定方法,其特征在于,包括以下步骤:
步骤一:将匀光板装在光场相机镜头的前端,对白光源进行拍照,得到光场相机的白图像;
步骤二:对光场相机的白图像进行去马赛克处理,将其转换成真彩图像;
步骤三:检测并提取真彩图像中所有微透镜的边缘轮廓,保存边缘点的坐标;
步骤四:对边缘点的坐标进行随机Hough变换,得到所有微透镜的圆心坐标Ci(x,y)和相对应的半径Ri,并统计微透镜个数N。
还包括对微透镜种类的标定步骤,该步骤包括:
步骤五:根据每个微透镜的圆心Ci(x,y)和半径Ri,分别用二维数组设置覆盖图像上微透镜对应区域的掩膜,将掩膜与原图像相乘,提取出每个微透镜的图像;
步骤六:利用公式(1)的锐度算法函数,对提取出来的每个微透镜图像进行清晰度Li计算,并将这些微透镜图像按照清晰度Li进行递增或递减排序;
式中,df表示图像灰度变化的幅值,dx表示相元之间的距离,df/dx计算时取该像素的8个邻域,M表示微透镜图像的像素个数;
步骤七:对清晰度Li按照公式(2)进行差分,并计算这些差分bi的平均值avg,设定阈值th=N·avg/15,若某处的差分值bi大于阈值th,则将该处作为微透镜种类分界点,将两个分界点之间所有微透镜归为一类;
bi=|Li+1-Li| (2)
步骤八:分类提取整幅图像中某类微透镜的图像,从而最终实现微透镜阵列中参数信息微透镜个数、中心位置、半径大小、微透镜种类及分布的标定。
利用SOBEL算子检测并提取图像中所有微透镜的边缘轮廓。
对边缘点进行随机Hough变换,得到所有微透镜的圆心坐标Ci(x,y)和相对应的半径Ri。
步骤二中,利用梯度线性插值法对光场相机的白图像进行去马赛克处理。
有益效果:本发明利用一种聚焦型光场相机的白图像,通过提取图像中所有微透镜的边缘轮廓确定微透镜的圆心坐标Ci(x,y)、相对应的半径Ri以及微透镜个数N,克服了现有灰度重心法和模板匹配法准确度不高的问题,还能对微透镜种类进行标定。
本发明标定方法基于数字图像处理技术,确定微透镜阵列的微透镜个数、中心位置、半径大小、微透镜种类及分布等微透镜阵列的重要参数信息,从而实现了光场相机微透镜多种参数的自动标定,可为光场相机内外参数标定及相机光学模型提供准确的微透镜阵列数据,及后续的图像处理提供不同类型的微透镜图像。
附图说明
图1是光场相机白图像。
图2是白图像对应的真彩图像。
图3微透镜图像边缘轮廓。
图4微透镜图像中心坐标。
图5是单一微透镜掩模。
图6是单一微透镜图像。
图7是按递增顺序的清晰度分布图。
图8是清晰度差值分布图。
图9 526nm微透镜图像。
图10 567nm微透镜图像。
图11 616nm微透镜图像。
具体实施方式
下面结合附图和具体实施例,进一步阐述本发明。应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读本发明之后,本领域技术人员对本发明的各种等价形式的修改落于本申请所附权利要求所限定的范围。本实施例中光场相机采用德国raytrix公司的R29型光场相机。
本发明一种光场相机微透镜阵列几何参数的标定方法,包括一下步骤:
步骤一:将匀光板装在光场相机镜头的前端,设置相机的F数(焦距与光圈直径比)为11和曝光时间0.8ms后,对无明显阴影白光源进行拍照,得到光场相机的白图像,如图1所示;
步骤二:利用梯度线性插值法对光场相机的白图像进行去马赛克处理,将其转换成真彩图像,如图2所示;
步骤三:利用SOBEL算子检测并提取真彩图像中所有微透镜的边缘轮廓如图3所示,保存边缘点的坐标;
步骤四:对边缘点进行随机Hough变换,得到所有微透镜的圆心坐标Ci(x,y)以及相对应的半径Ri,如图4所示,并统计微透镜个数N;
步骤五:根据每个微透镜的圆心Ci(x,y)和半径Ri,分别用二维数组设置覆盖图像上微透镜对应区域的掩膜,如图5所示,将掩膜与原图像相乘,提取出每个微透镜的图像,如图6所示;
步骤六:利用形如公式(1)的锐度算法函数,对提取出来的每个微透镜图像进行清晰度Li计算,并将这些微透镜图像按照清晰度Li进行递增排序,如图7所示;
式中,df表示图像灰度变化的幅值,dx表示相元之间的距离,df/dx计算时取该像素的8个邻域,M表示微透镜图像的像素数。
步骤七:对清晰度Li按照公式(2)进行差分,其差分值bi分布如图8所示,并计算这些差分bi的平均值avg,设定阈值th=N·avg/15,若某处的差分值bi大于阈值th,则将该处作为微透镜种类分界点,将两个分界点之间所有微透镜归为一类;
bi=|Li+1-Li| (2)
步骤八:分类提取整幅图像中各类微透镜的图像,如图9-11所示。本实施例采用的光场相机的微透镜的种类为3类,分别为526μm,567μm和616μm。因为微透镜焦距越大,在相机的探测器上的离焦程度越大,所以最清晰的为526μm透镜,567μm次之,616μm最模糊。至此,最终实现微透镜阵列中微透镜个数、中心位置、半径大小、焦距种类及分布等微透镜阵列的重要参数信息的自动标定。
Claims (5)
1.一种光场相机微透镜阵列几何参数的标定方法,其特征在于,包括以下步骤:
步骤一:将匀光板装在光场相机镜头的前端,对白光源进行拍照,得到光场相机的白图像;
步骤二:对光场相机的白图像进行去马赛克处理,将其转换成真彩图像;
步骤三:检测并提取真彩图像中所有微透镜的边缘轮廓,保存边缘点的坐标;
步骤四:对边缘点的坐标进行随机Hough变换,得到所有微透镜的圆心坐标Ci(x,y)和相对应的半径Ri,并统计微透镜个数N。
2.根据权利要求1所述的光场相机微透镜阵列几何参数的标定方法,其特征在于,还包括对微透镜种类的标定步骤,该步骤包括:
步骤五:根据每个微透镜的圆心Ci(x,y)和半径Ri,分别用二维数组设置覆盖图像上微透镜对应区域的掩膜,将掩膜与原图像相乘,提取出每个微透镜的图像;
步骤六:利用公式(1)的锐度算法函数,对提取出来的每个微透镜图像进行清晰度Li计算,并将这些微透镜图像按照清晰度Li进行递增或递减排序;
式中,df表示图像灰度变化的幅值,dx表示相元之间的距离,df/dx计算时取该像素的8个邻域,M表示微透镜图像的像素个数;
步骤七:对清晰度Li按照公式(2)进行差分,并计算这些差分bi的平均值avg,设定阈值th=N·avg/15,若某处的差分值bi大于阈值th,则将该处作为微透镜种类分界点,将两个分界点之间所有微透镜归为一类;
bi=|Li+1-Li| (2)
步骤八:分类提取整幅图像中某类微透镜的图像,从而最终实现微透镜阵列中参数信息微透镜个数、中心位置、半径大小、微透镜种类及分布的标定。
3.根据权利要求2所述的光场相机微透镜阵列几何参数的标定方法,其特征在于,步骤三中,利用SOBEL算子检测并提取图像中所有微透镜的边缘轮廓。
4.根据权利要求3所述的光场相机微透镜阵列几何参数的标定方法,其特征在于,步骤四中,对边缘点进行随机Hough变换,得到所有微透镜的圆心坐标Ci(x,y)和相对应的半径Ri。
5.根据权利要求4所述的光场相机微透镜阵列几何参数的标定方法,其特征在于,步骤二中,利用梯度线性插值法对光场相机的白图像进行去马赛克处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611014733.6A CN106651959B (zh) | 2016-11-15 | 2016-11-15 | 一种光场相机微透镜阵列几何参数的标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611014733.6A CN106651959B (zh) | 2016-11-15 | 2016-11-15 | 一种光场相机微透镜阵列几何参数的标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106651959A true CN106651959A (zh) | 2017-05-10 |
CN106651959B CN106651959B (zh) | 2019-05-31 |
Family
ID=58809101
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611014733.6A Active CN106651959B (zh) | 2016-11-15 | 2016-11-15 | 一种光场相机微透镜阵列几何参数的标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106651959B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107492127A (zh) * | 2017-09-18 | 2017-12-19 | 丁志宇 | 光场相机参数标定方法、装置、存储介质和计算机设备 |
CN108305233A (zh) * | 2018-03-06 | 2018-07-20 | 哈尔滨工业大学 | 一种针对微透镜阵列误差的光场图像校正方法 |
CN108426585A (zh) * | 2018-03-12 | 2018-08-21 | 哈尔滨工业大学 | 一种光场相机的几何标定方法 |
CN109859261A (zh) * | 2019-01-07 | 2019-06-07 | 宁波舜宇光电信息有限公司 | 一种鱼眼镜头的光心定位方法和装置、及摄像模组的质量检测方法 |
CN112834180A (zh) * | 2020-12-28 | 2021-05-25 | 中国科学院紫金山天文台 | 用于探测器阵列芯片和微透镜阵列的安装装置和定位方法 |
CN114136591A (zh) * | 2021-11-26 | 2022-03-04 | 广东工业大学 | 一种微透镜阵列超精密加工误差评定方法与系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005008174A1 (ja) * | 2003-07-23 | 2005-01-27 | Brother Kogyo Kabushiki Kaisha | 3次元形状検出装置、撮像装置、及び、3次元形状検出方法 |
CN101334263A (zh) * | 2008-07-22 | 2008-12-31 | 东南大学 | 圆形目标的圆心定位方法 |
US20140294290A1 (en) * | 2013-03-28 | 2014-10-02 | Texas Instruments Incorporated | Projector-Camera Misalignment Correction for Structured Light Systems |
EP2909811A1 (en) * | 2012-10-17 | 2015-08-26 | Cathx Research Ltd. | Improvements in and relating to gathering range and dimensional information underwater surveys |
CN105374044A (zh) * | 2015-12-04 | 2016-03-02 | 中国科学院光电技术研究所 | 一种光场相机的自动标定方法 |
CN105488810A (zh) * | 2016-01-20 | 2016-04-13 | 东南大学 | 一种聚焦光场相机内外参数标定方法 |
CN105931252A (zh) * | 2016-05-16 | 2016-09-07 | 大连理工大学 | 一种基于几何约束的椭圆快速检测方法 |
-
2016
- 2016-11-15 CN CN201611014733.6A patent/CN106651959B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005008174A1 (ja) * | 2003-07-23 | 2005-01-27 | Brother Kogyo Kabushiki Kaisha | 3次元形状検出装置、撮像装置、及び、3次元形状検出方法 |
CN101334263A (zh) * | 2008-07-22 | 2008-12-31 | 东南大学 | 圆形目标的圆心定位方法 |
EP2909811A1 (en) * | 2012-10-17 | 2015-08-26 | Cathx Research Ltd. | Improvements in and relating to gathering range and dimensional information underwater surveys |
US20140294290A1 (en) * | 2013-03-28 | 2014-10-02 | Texas Instruments Incorporated | Projector-Camera Misalignment Correction for Structured Light Systems |
CN105374044A (zh) * | 2015-12-04 | 2016-03-02 | 中国科学院光电技术研究所 | 一种光场相机的自动标定方法 |
CN105488810A (zh) * | 2016-01-20 | 2016-04-13 | 东南大学 | 一种聚焦光场相机内外参数标定方法 |
CN105931252A (zh) * | 2016-05-16 | 2016-09-07 | 大连理工大学 | 一种基于几何约束的椭圆快速检测方法 |
Non-Patent Citations (2)
Title |
---|
DONGHYEON CHO ET AL.: ""Modeling the Calibration Pipeline of the Lytro Camera for High Quality Light-Field Image Reconstruction"", 《IEEEXPLORE》 * |
TAKAHIRO SAITO ET AL.: ""NONLINEAR DECOMPOSITION-AND-DEMOSAICING APPROACH FOR A DIGITAL COLOR CAMERA"", 《IEEEXPLORE》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107492127A (zh) * | 2017-09-18 | 2017-12-19 | 丁志宇 | 光场相机参数标定方法、装置、存储介质和计算机设备 |
CN108305233A (zh) * | 2018-03-06 | 2018-07-20 | 哈尔滨工业大学 | 一种针对微透镜阵列误差的光场图像校正方法 |
CN108426585A (zh) * | 2018-03-12 | 2018-08-21 | 哈尔滨工业大学 | 一种光场相机的几何标定方法 |
CN108426585B (zh) * | 2018-03-12 | 2019-09-13 | 哈尔滨工业大学 | 一种光场相机的几何标定方法 |
CN109859261A (zh) * | 2019-01-07 | 2019-06-07 | 宁波舜宇光电信息有限公司 | 一种鱼眼镜头的光心定位方法和装置、及摄像模组的质量检测方法 |
CN112834180A (zh) * | 2020-12-28 | 2021-05-25 | 中国科学院紫金山天文台 | 用于探测器阵列芯片和微透镜阵列的安装装置和定位方法 |
CN114136591A (zh) * | 2021-11-26 | 2022-03-04 | 广东工业大学 | 一种微透镜阵列超精密加工误差评定方法与系统 |
CN114136591B (zh) * | 2021-11-26 | 2023-09-08 | 广东工业大学 | 一种微透镜阵列超精密加工误差评定方法与系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106651959B (zh) | 2019-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106651959A (zh) | 一种光场相机微透镜阵列几何参数的标定方法 | |
EP3516625B1 (en) | A device and method for obtaining distance information from views | |
Boltes et al. | Automatic extraction of pedestrian trajectories from video recordings | |
CN103345755B (zh) | 一种基于Harris算子的棋盘格角点亚像素提取方法 | |
CN104537659B (zh) | 双摄像机的自动标定方法和系统 | |
CN102801929B (zh) | 图像传感器和摄像设备 | |
CN108426585B (zh) | 一种光场相机的几何标定方法 | |
CN102663732B (zh) | 一种针对光场相机进行相对辐射定标的方法 | |
WO2018018987A1 (zh) | 一种适用于光场相机的标定预处理方法 | |
CN108776980A (zh) | 一种面向微透镜光场相机的标定方法 | |
CN107610182B (zh) | 一种光场相机微透镜阵列中心的标定方法 | |
CN109194954B (zh) | 鱼眼摄像头性能参数测试方法、装置、设备及可存储介质 | |
CN106210465A (zh) | 用于根据4d原始光场数据获得重聚焦图像的方法 | |
CN109496316B (zh) | 图像识别系统 | |
US10366478B2 (en) | Method and device for obtaining a HDR image by graph signal processing | |
CN110662014B (zh) | 一种光场相机四维数据大景深三维显示的方法 | |
CN102867295A (zh) | 一种彩色图像颜色校正方法 | |
CN109635799B (zh) | 一种燃气表字轮数字的识别方法 | |
CN105374044A (zh) | 一种光场相机的自动标定方法 | |
JP5850648B2 (ja) | 撮像装置 | |
US20220046189A1 (en) | Method for Improved Acquisition of Images for Photogrammetry | |
CN113506275B (zh) | 一种基于全景城市图像处理方法 | |
CN106973199A (zh) | 利用对焦距离扫描提高深度精确度的多光圈相机系统 | |
CN105654462A (zh) | 一种基于影像配准的建筑物高程提取方法 | |
CN111325218B (zh) | 基于光场图像的Hog特征检测与匹配方法 |
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 |