CN113012053B - 一种集成成像桌面3d显示的360°视差图像校正方法 - Google Patents
一种集成成像桌面3d显示的360°视差图像校正方法 Download PDFInfo
- Publication number
- CN113012053B CN113012053B CN201911322712.4A CN201911322712A CN113012053B CN 113012053 B CN113012053 B CN 113012053B CN 201911322712 A CN201911322712 A CN 201911322712A CN 113012053 B CN113012053 B CN 113012053B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- parallax image
- coordinates
- corner points
- camera
- 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
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000003384 imaging method Methods 0.000 title claims abstract description 21
- 238000003702 image correction Methods 0.000 title claims abstract description 10
- 230000009466 transformation Effects 0.000 claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 12
- 238000012216 screening Methods 0.000 claims abstract description 11
- 230000001131 transforming effect Effects 0.000 claims abstract description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G06T5/80—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20228—Disparity calculation for image-based rendering
Abstract
本发明提出一种集成成像桌面3D显示的360°视差图像校正方法。所述方法根据每个摄像机的相对旋转角度,以检测到的棋盘格标定板所有角点的质心为原点,在拍摄的标定板视差图像上建立新坐标系,并将角点坐标变换到新坐标系中,通过检测角点距坐标系原点的距离,对四个边缘角点进行精确筛选。之后,结合单应性变换理论,计算视差图像对应的单应性变换矩阵,对桌面3D场景的360°视差图像进行校正。所述方法具有旋转不变性,可实现任意旋转角度下边缘角点的准确定位和视差图像的正确校正,能够为集成成像桌面3D显示提供正确的数据来源。
Description
一、技术领域
本发明涉及3D(三维)图像拍摄技术,特别涉及一种集成成像桌面3D显示的360°视差图像校正方法。
二、背景技术
集成成像桌面3D显示技术主要包括空分复用集成成像桌面3D显示和时分复用集成成像桌面3D显示两种。空分复用集成成像桌面3D显示技术通常采用光线调制模块对桌面显示屏发出的光线向桌面四周进行偏转,使得环绕在桌面3D显示系统周围的观看者都能看到3D图像。时分复用集成成像桌面3D显示技术利用单一方向偏转的光线调制模块将桌面显示屏发出的光线向桌面一侧偏转,通过光线调制模块的旋转,配合显示屏上微图像阵列的同步切换,实现360°环视3D观看效果。
对应地,在集成成像桌面3D信息获取过程中,需要环绕桌面3D场景拍摄360°范围的视差图像,并利用桌面棋盘格标定板对拍摄的视差图像进行校正,以克服摄像机投影平面与桌面所在平面不对齐带来的拍摄错误。然而,在360°环绕拍摄过程中,棋盘格标定板相对摄像机进行了不同角度的旋转,当旋转角度较大时,棋盘格标定板的角点检测结果排序混乱,无法与棋盘格标定板的相对旋转角度统一,导致边缘角点无法被准确定位,校正后的视差图像出现方向错误,无法为集成成像桌面3D显示提供正确的数据来源。此外,摄像机光轴相对标定板在垂直方向上的偏转也会对角点排序结果造成较大影响。
三、发明内容
为解决上述问题,本发明提出一种集成成像桌面3D显示的360°视差图像校正方法,所述方法根据每个摄像机的相对旋转角度,以检测到的棋盘格标定板所有角点的质心为原点,在拍摄的标定板视差图像上建立新坐标系,并将角点坐标变换到新坐标系中,通过检测角点距坐标系原点的距离,对四个边缘角点进行精确筛选。之后,结合单应性变换理论,计算视差图像对应的单应性变换矩阵,对桌面3D场景的视差图像进行校正,用于集成成像桌面3D显示。
所述方法包括以下六个步骤。
第一步,确定每个摄像机的相对旋转角度。
第二步,建立标定板视差图像的新坐标系。
第三步,变换标定板角点坐标到新坐标系。
第四步,筛选边缘角点、获取边缘角点坐标。
第五步,计算每幅视差图像的单应性变换矩阵。
第六步,校正桌面3D场景的360°视差图像。
所述第一步,确定每个摄像机的相对旋转角度步骤中,摄像机阵列围绕桌面3D场景呈360°环形排布。所述摄像机阵列包含的摄像机个数为M×N,其中M表示水平环向上360°范围内包含的摄像机个数,N表示径向方向上包含的摄像机个数。摄像机的索引表示为(m,n)。将第(0,n)个摄像机在水平环向上的角度设为起始0°,第(m,n)个摄像机相对第(0,n)个摄像机的旋转角度θm,n表示为:
所述第二步,建立标定板视差图像的新坐标系步骤中,首先确定标定板视差图像上新坐标系的原点位置。在第(m,n)幅标定板视差图像中,利用亚像素精度角点检测技术,自动识别I×J个棋盘格标定板角点,并确定角点Pi的坐标(xi,yi),其中坐标以像素为单位,i表示所述角点的索引,i∈{1,2,3,…,I×J}。根据定位的角点坐标,计算所有角点的质心Om,n坐标(xo,yo),其满足:
然后,以质心Om,n为原点,结合旋转角度θm,n,在第(m,n)幅标定板视差图像上建立新的直角坐标系x'-y'。所述新坐标系x'-y'相对原始坐标系x-y进行了平移和旋转变换,其中水平偏移量为xo,垂直偏移量为yo,旋转角度为θm,n,该角度代表了坐标系绕质心Om,n旋转的角度,顺时针方向为正。
所述第三步,变换标定板角点坐标到新坐标系步骤中,变换后的角点坐标(x'i,y'i)与初始角点坐标(xi,yi)的关系满足:
所述第四步,筛选边缘角点、获取边缘角点坐标的步骤中,首先计算角点到质心Om,n的距离di,表示为:
四个边缘角点A、B、C、D分别处于新建坐标系的不同象限中,根据变换后的x'坐标和y'坐标的正负值,确定四个边缘角点坐标(x'A,y'A)、(x'B,y'B)、(x'C,y'C)、(x'D,y'D)的范围,并对每个范围内所有角点到质心的距离进行排序,筛选出距离最大的角点,该角点即为对应的边缘角点,其坐标满足:
然后,将边缘角点A、B、C、D在坐标系x'-y'下的坐标变换回原始坐标系下,获得原始坐标系下对应的坐标(xA,yA)、(xB,yB)、(xC,yC)和(xD,yD)。
所述第五步,计算每幅视差图像的单应性变换矩阵步骤中,根据筛选的边缘角点坐标(xA,yA)、(xB,yB)、(xC,yC)和(xD,yD),结合校正后视差图像的尺寸W×H,基于单应性变换原理,计算M×N幅视差图像对应的3×3单应性变换矩阵{Hm,n|m∈{1,2,3,…,M},n∈{1,2,3,…,N}}。
所述第六步,校正桌面3D场景的360°视差图像步骤中,利用相同位姿的摄像机阵列对桌面3D场景进行拍摄,第(m,n)幅桌面3D场景的视差图像表示为Em,n(x',y'),利用单应性变换矩阵Hm,n计算得到对应的校正后视差图像E'm,n(x,y),表示为:
E'm,n(x',y')=Em,n(x,y) (7)
根据本发明提出的一种集成成像桌面3D显示的360°视差图像校正方法,通过在标定板视差图像上建立新的坐标系,对边缘角点进行精确筛选,克服了校正后视差图像错误的问题,能够为集成成像桌面3D显示提供正确的数据来源。所述方法具有旋转不变性,能够实现任意旋转角度下边缘角点的准确定位和视差图像的正确校正。
四、附图说明
本发明的前述方面及优点从下述结合附图与实施例的详细描述中将得以进一步明确和容易理解,其中:
附图1为根据本申请一实施例的一种集成成像桌面3D显示的360°视差图像校正方法。
附图2A为根据本申请一实施例的摄像机阵列拍摄桌面3D场景的示意图。
附图2B为根据图2A所示的摄像机相对旋转角度的示意图。
附图3为根据本申请一实施例的(a)第(0,2)幅标定板视差图像的新建坐标系和(b)第(46,2)幅标定板视差图像的新建坐标系的示意图。
附图4为(a)根据本申请一实施例的第(46,2)幅校正后的视差图像与(b)传统未经角点筛选生成的第(46,2)幅校正后的视差图像的对比图。
上述附图中的图示标号为:
1摄像机阵列,2桌面3D场景,3摄像机,4摄像机阵列轨迹,5第(0,2)幅标定板视差图像,6第(46,2)幅标定板视差图像。
应该理解上述附图只是示意性的,并没有按比例绘制。
五、具体实施方式
为了便于理解本申请,下面将参照相关附图对本申请进行更全面的描述。附图中给出了本申请的可选实施例。但是,本申请可以以许多不同的形式来实现,并不限于本文所描述的实施例。相反地,提供这些实施例的目的是使对本申请的公开内容的理解更加透彻全面。
除非另有定义,本文所使用的所有技术和科学术语与属于发明的技术领域的技术人员通常理解的含义相同。本文中在发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在限制本申请。本文所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。
本文所使用的方位术语“垂直”、“水平”、“边缘”以及类似的表述只是为了说明的目的,并不表示是唯一的实施方式也并非对于本申请的限制。
下面将结合本申请公开的实施例及附图对本申请所提出的一种集成成像桌面3D显示的360°视差图像校正方法进行详细的描述。
附图1示出了根据本申请一实施例的一种集成成像桌面3D显示的360°视差图像校正方法,该方法包括以下步骤。
步骤S100,确定每个摄像机的相对旋转角度。
步骤S200,建立标定板视差图像的新坐标系。
步骤S300,变换标定板角点坐标到新坐标系。
步骤S400,筛选边缘角点、获取边缘角点坐标。
步骤S500,计算每幅视差图像的单应性变换矩阵。
步骤S600,校正桌面3D场景的360°视差图像。
所述第一步,确定每个摄像机的相对旋转角度步骤中,摄像机阵列1围绕桌面3D场景2呈360°环形排布,如附图2所示。在一个实施例中,桌面3D场景2为一个骰子,摄像机阵列1包含的摄像机个数M×N可以为300×8,其中300表示水平环向上360°范围内包含的摄像机个数,8表示径向方向上包含的摄像机个数。摄像机的索引表示为(m,n)。如附图2B所示,在一个实施例中,将第(0,n)个摄像机3在水平环向上的角度设为起始0°,第(m,n)个摄像机3相对第(0,n)个摄像机3的旋转角度θm,n表示为:
在一个实施例中,第(46,2)个摄像机相对第(0,2)个摄像机的旋转角度θ46,2为55.2°。
所述第二步,建立标定板视差图像的新坐标系步骤中,首先确定标定板视差图像上新坐标系的原点位置。在一个实施例中,对于第(46,2)幅标定板视差图像,利用亚像素精度角点检测技术,自动识别I×J=9×9个棋盘格标定板角点,并确定角点Pi的坐标(xi,yi),其中坐标以像素为单位,i表示所述角点的索引,i∈{1,2,3,…,81}。在一个实施例中,第(46,2)幅标定板视差图像上,第3个角点P3的坐标(x3,y3)为(161,567)。根据定位的角点坐标,计算所有角点的质心Om,n坐标(xo,yo),其满足:
在一个实施例中,第(46,2)幅标定板视差图像中,所有角点的质心O46,2坐标(xo,yo)为(501,540)。然后,以质心Om,n为原点,结合旋转角度θm,n,在第(m,n)幅标定板视差图像上建立新的直角坐标系x'-y'。附图3示出了根据本申请一实施例的(a)第(0,2)幅标定板视差图像5的新建坐标系和(b)第(46,2)幅标定板视差图像6的新建坐标系的示意图。所述新坐标系x'-y'相对原始坐标系x-y进行了平移和旋转变换。在其中一个实施例中,第(46,2)幅标定板视差图像6内,新坐标系的水平偏移量xo为501个像素,垂直偏移量yo为540个像素,旋转角度θ46,2为55.2°,该角度代表了坐标系绕质心O46,2旋转的角度,顺时针方向为正。
所述第三步,变换标定板角点坐标到新坐标系步骤中,变换后的角点坐标(x'i,y'i)与初始角点坐标(xi,yi)的关系满足:
在一个实施例中,第(46,2)幅标定板视差图像6内,第3个角点P3变换后的坐标(x'3,y'3)为(-172,294)。
所述第四步,筛选边缘角点、获取边缘角点坐标的步骤中,首先计算角点到质心Om,n的距离di,表示为:
在一个实施例中,第(46,2)幅标定板视差图像6内,角点P3到质心O46,2的距离d3为340.5个像素。四个边缘角点A、B、C、D分别处于新建坐标系的不同象限中,根据变换后的x'坐标和y'坐标的正负值,确定四个边缘角点坐标(x'A,y'A)、(x'B,y'B)、(x'C,y'C)、(x'D,y'D)的范围,并对每个范围内所有角点到质心的距离进行排序,筛选出距离最大的角点,该角点即为对应的边缘角点,其坐标满足:
在一个实施例中,第(46,2)幅标定板视差图像6内,筛选的四个边缘角点A、B、C、D的索引分别为ia=73,ib=81,ic=9,id=1,角点A的坐标为(x'A,y'A)=(-187,-210),角点B的坐标为(x'B,y'B)=(284,-328),角点C的坐标为(x'C,y'C)=(231,268),角点D的坐标为(x'D,y'D)=(-285,301)。然后,将边缘角点A、B、C、D在在坐标系x'-y'下的坐标变换回原始坐标系下,获得原始坐标系下各自对应的坐标。在一个实施例中,第(46,2)幅标定板视差图像6内,(xA,yA)=(567,267)、(xB,yB)=(932,587)、(xC,yC)=(412,883)和(xD,yD)=(91,478)。
所述第五步,计算每幅视差图像的单应性变换矩阵步骤中,根据筛选的边缘角点坐标,结合校正后视差图像的尺寸W×H=1000×1000,基于单应性变换原理,计算M×N=300×8幅视差图像对应的3×3单应性变换矩阵{Hm,n|m∈{1,2,3,…,300},n∈{1,2,3,…,8}}。以第(46,2)幅视差图像为例,其对应的3×3单应性变换矩阵H46,2计算为
所述第六步,校正桌面3D场景的360°视差图像步骤中,利用相同位姿的摄像机阵列对桌面3D场景进行拍摄,第(m,n)幅桌面3D场景的视差图像表示为Em,n(x',y'),利用单应性变换矩阵Hm,n计算得到对应的校正后视差图像E'm,n(x,y),表示为:
E'm,n(x',y')=Em,n(x,y) (7)
附图4示出了(a)根据本申请一实施例的第(46,2)幅校正后的视差图像与(b)传统未经角点筛选生成的第(46,2)幅校正后的视差图像的对比图。根据本申请一实施例校正的视差图像方向正确,能够为集成成像桌面3D显示提供正确的数据来源。
Claims (5)
1.一种集成成像桌面3D显示的360°视差图像校正方法,其特征在于,所述方法根据每个摄像机的相对旋转角度,以检测到的棋盘格标定板所有角点的质心为原点,在拍摄的标定板视差图像上建立新坐标系,并将角点坐标变换到新坐标系中,通过检测角点距坐标系原点的距离,对四个边缘角点进行精确筛选,之后对桌面3D场景的360°视差图像进行校正;所述方法包括以下六个步骤:第一步,确定每个摄像机的相对旋转角度;第二步,建立标定板视差图像的新坐标系;第三步,变换标定板角点坐标到新坐标系;第四步,筛选边缘角点,获取边缘角点坐标;第五步,计算每幅视差图像的单应性变换矩阵;第六步,校正桌面3D场景的360°视差图像。
3.根据权利要求1所述的一种集成成像桌面3D显示的360°视差图像校正方法,其特征在于,所述第二步,建立标定板视差图像的新坐标系步骤中,在第(m,n)幅标定板视差图像内,自动识别I×J个棋盘格标定板角点,并确定角点Pi的坐标(xi,yi),其中坐标以像素为单位,i表示所述角点的索引,i∈{1,2,3,…,I×J},根据定位的角点坐标,计算所有角点的质心Om,n坐标(xo,yo),其满足然后以质心Om,n为原点,结合摄像机的相对旋转角度θm,n,在第(m,n)幅标定板视差图像上建立新的直角坐标系x'-y',所述新坐标系x'-y'相对原始坐标系x-y进行了平移和旋转变换,其中水平偏移量为xo,垂直偏移量为yo,旋转角度为θm,n,该角度代表了坐标系绕质心Om,n旋转的角度,顺时针方向为正。
5.根据权利要求1所述的一种集成成像桌面3D显示的360°视差图像校正方法,其特征在于,所述第四步,筛选边缘角点,获取边缘角点坐标的步骤中,首先计算角点(x'i,y'i)到质心Om,n的距离di,表示为:四个边缘角点A、B、C、D分别处于新建坐标系的不同象限中,根据变换后的x'坐标和y'坐标的正负值,确定四个边缘角点坐标(x'A,y'A)、(x'B,y'B)、(x'C,y'C)、(x'D,y'D)的范围,并对每个范围内所有角点到质心的距离进行排序,筛选出距离最大的角点,该角点即为对应的边缘角点,其坐标满足然后将边缘角点A、B、C、D在新坐标系x'-y'下的坐标变换回原始坐标系下,获得原始坐标系下对应的坐标(xA,yA)、(xB,yB)、(xC,yC)和(xD,yD)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911322712.4A CN113012053B (zh) | 2019-12-20 | 2019-12-20 | 一种集成成像桌面3d显示的360°视差图像校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911322712.4A CN113012053B (zh) | 2019-12-20 | 2019-12-20 | 一种集成成像桌面3d显示的360°视差图像校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113012053A CN113012053A (zh) | 2021-06-22 |
CN113012053B true CN113012053B (zh) | 2022-05-27 |
Family
ID=76382116
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911322712.4A Active CN113012053B (zh) | 2019-12-20 | 2019-12-20 | 一种集成成像桌面3d显示的360°视差图像校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113012053B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105303574A (zh) * | 2015-07-30 | 2016-02-03 | 四川大学 | 一种基于单应性变换的集成成像摄像机阵列标定方法 |
CN107016707A (zh) * | 2017-04-13 | 2017-08-04 | 四川大学 | 一种集成成像超大三维场景拍摄图像校正方法 |
CN108881881A (zh) * | 2018-08-03 | 2018-11-23 | 北京航空航天大学 | 一种集成成像桌面3d显示的自适应微图像阵列生成方法 |
CN109147036A (zh) * | 2018-08-03 | 2019-01-04 | 北京航空航天大学 | 一种基于深度学习的集成成像微图像阵列快速生成方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103827917B (zh) * | 2011-07-25 | 2017-06-09 | 科英布拉大学 | 用于使用棋盘图案的一幅或多幅图像的自动相机校准的方法和装置 |
-
2019
- 2019-12-20 CN CN201911322712.4A patent/CN113012053B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105303574A (zh) * | 2015-07-30 | 2016-02-03 | 四川大学 | 一种基于单应性变换的集成成像摄像机阵列标定方法 |
CN107016707A (zh) * | 2017-04-13 | 2017-08-04 | 四川大学 | 一种集成成像超大三维场景拍摄图像校正方法 |
CN108881881A (zh) * | 2018-08-03 | 2018-11-23 | 北京航空航天大学 | 一种集成成像桌面3d显示的自适应微图像阵列生成方法 |
CN109147036A (zh) * | 2018-08-03 | 2019-01-04 | 北京航空航天大学 | 一种基于深度学习的集成成像微图像阵列快速生成方法 |
Non-Patent Citations (1)
Title |
---|
Plannar Parallax Based Camera Array Calibration Method for Integral Imaging Three-dimensional Information Acquirement;Zhao-Long Xiong 等;《SID 2016 DIGEST》;20161231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113012053A (zh) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11269244B2 (en) | System and method for calibrating a display system using manual and semi-manual techniques | |
CN110211043B (zh) | 一种用于全景图像拼接的基于网格优化的配准方法 | |
US7479982B2 (en) | Device and method of measuring data for calibration, program for measuring data for calibration, program recording medium readable with computer, and image data processing device | |
US7379619B2 (en) | System and method for two-dimensional keystone correction for aerial imaging | |
US6527395B1 (en) | Method for calibrating a projector with a camera | |
CN110809786B (zh) | 校准装置、校准图表、图表图案生成装置和校准方法 | |
US7019713B2 (en) | Methods and measurement engine for aligning multi-projector display systems | |
US9195121B2 (en) | Markerless geometric registration of multiple projectors on extruded surfaces using an uncalibrated camera | |
CN111750820A (zh) | 影像定位方法及其系统 | |
US9892488B1 (en) | Multi-camera frame stitching | |
KR20110128564A (ko) | 카메라 시스템의 영상 보정 방법 | |
CN107527336B (zh) | 镜头相对位置标定方法及装置 | |
WO2005084017A1 (ja) | マルチプロジェクションシステム | |
CN111083457A (zh) | 多光机投影图像的校正方法、装置和多光机投影仪 | |
CN110505468B (zh) | 一种增强现实显示设备的测试标定及偏差修正方法 | |
CN102222337A (zh) | 用于图像拼接的鱼眼图像校正方法 | |
CN106952219B (zh) | 一种基于外参数修正鱼眼像机的图像生成方法 | |
CN112598747A (zh) | 一种单目相机与投影仪的联合标定方法 | |
CN106886976B (zh) | 一种基于内参数修正鱼眼像机的图像生成方法 | |
CN112734860A (zh) | 一种基于弧形幕先验信息的逐像素映射投影几何校正方法 | |
CN107358577A (zh) | 一种立方体全景图的快速拼接方法 | |
CN113012053B (zh) | 一种集成成像桌面3d显示的360°视差图像校正方法 | |
CN113793266A (zh) | 一种多目机器视觉图像拼接方法、系统及存储介质 | |
CN112258581B (zh) | 一种多鱼眼镜头全景相机的现场标定方法 | |
CN108205799B (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 |