CN113989105B - 一种单相机球面镜反射成像投影装置 - Google Patents
一种单相机球面镜反射成像投影装置 Download PDFInfo
- Publication number
- CN113989105B CN113989105B CN202111246846.XA CN202111246846A CN113989105B CN 113989105 B CN113989105 B CN 113989105B CN 202111246846 A CN202111246846 A CN 202111246846A CN 113989105 B CN113989105 B CN 113989105B
- Authority
- CN
- China
- Prior art keywords
- camera
- spherical
- dimensional
- point
- image
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 87
- 238000013507 mapping Methods 0.000 claims abstract description 12
- 230000005540 biological transmission Effects 0.000 claims abstract description 7
- 238000004088 simulation Methods 0.000 claims description 40
- 239000011159 matrix material Substances 0.000 claims description 16
- 238000000034 method Methods 0.000 claims description 13
- 239000013598 vector Substances 0.000 claims description 11
- 230000003287 optical effect Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 238000013519 translation Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 3
- 230000001939 inductive effect Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 description 12
- 238000011160 research Methods 0.000 description 6
- 230000007547 defect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000013256 coordination polymer Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000011514 reflex Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/06—Topological mapping of higher dimensional structures onto lower dimensional surfaces
- G06T3/067—Reshaping or unfolding 3D tree structures onto 2D planes
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03B—APPARATUS OR ARRANGEMENTS FOR TAKING PHOTOGRAPHS OR FOR PROJECTING OR VIEWING THEM; APPARATUS OR ARRANGEMENTS EMPLOYING ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ACCESSORIES THEREFOR
- G03B17/00—Details of cameras or camera bodies; Accessories therefor
- G03B17/02—Bodies
- G03B17/17—Bodies with reflectors arranged in beam forming the photographic image, e.g. for reducing dimensions of camera
-
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
- H04N23/60—Control of cameras or camera modules
- H04N23/64—Computer-aided capture of images, e.g. transfer from script file into camera, check of taken image quality, advice or proposal for image composition or decision on when to take image
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Studio Devices (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Lenses (AREA)
Abstract
本发明涉及一种单相机球面镜反射成像投影装置,其由球面反射镜、摄像机、数据传输线和计算机组成。本发明成像映射模型由三维空间点P经球面反射镜形成球面反射镜像虚拟点,由摄像机拍摄,在图像平面中形成图像点,图像由数据传输线存储到计算机内,根据通用单相机球面镜反射成像模型,计算出图像点的图像像素坐标。本发明首次建立了任意位置和角度下,单相机球面镜反射成像的定量投影模型;实现了空间点经球面镜反射至摄像机平面图像点的投影成像。
Description
技术领域
本发明属于视觉成像技术领域,涉及一种单相机球面镜反射成像投影装置。
背景技术
随着视觉测量的应用需求不断提升,测量场景复杂化,对立体视觉传感器提出了低成本、小型化、轻量化、大视场、快速测量等要求。而目前的双目视觉测量方法,主要针对光照均匀、无遮挡或较少反射的强约束场景,仅适用于特定条件下的三维测量需求,如室内环境。单相机镜像双目视觉测量模式融合了传统双目视觉测量的优点,对比传统双目测量模式,具有系统构成简单、组件配置灵活、可控参数多、结构易小型化、可集成等优点,双目图像左右视图的采集具有天然的同步性,避免了传统双目相机难以同步采集的弊端。同时,由于镜像双目视觉测量模式易于集成,逐渐出现单相机镜像双目传感器的相关研究。镜像双目传感器由单个工业相机与光学折反射元件组合构成,大大降低了成本,提高了视觉传感器在受限空间及复杂场景的适用性。
传统的双目立体视觉系统通常由两相机构成,体积大、操作复杂、测量速度慢,且难以保证图像采集的同时性,为解决这些问题,可采用单相机与光学折反射元件组合的方式构建镜像双目立体视觉系统,缩小了测量系统体积,提高了测量速度。但由于平面折反射镜像双目系统在相同尺寸图像上呈现两个像,导致测量视场减半,使得平面镜的系统视场范围较小,难以满足宽视场测量需求。针对这种弊端,国内外研究人员对宽视场成像技术、基于曲面折反射原理的镜像双目展开研究,基于镜像立体视觉原理,进一步研究单相机凸面镜新型镜像立体视觉测量系统,旨在使小型视觉传感器具有更宽的视场范围。利用特殊镜头,或将测量系统中的光学折反射元件由平面镜替换为凸曲面镜,对空间场景进行成像,视场范围显著扩大。典型的曲面镜类型有:抛物面镜、双曲面镜、椭球面镜、锥面镜、球面镜等。
国内外许多研究机构和学者对折反射相机成像进行了研究,建立了各种单相机折反射成像模型。1999年Baker和Nayar为解决传统摄像机视场范围受限的问题,采用折反射元件与摄像机组合扩展视场范围,分析了不同面型的折反射投影成像,推导了不同面型折反射结构的单视点限制,建立了单视点成像模型;2001年Geyer分析了折反射透视几何,针对折反射相机建立了透视投影模型,证明了折反射投影中点与线存在对偶关系,并采用抛物面镜从任意位置的单视点直线估计了焦距和图像中心;2004年Mei提出了一种利用平面网格标定单视点全向传感器的方法,通过标准圆和标准平面建立统一成像模型。以上研究均针对单相机折反射结构建立了投影模型,但缺陷在于,模型中要求摄像机与折反射镜共轴,然而在实际装配中难以保证,从而产生成像偏差,同时,中心折反射结构存在应用局限性,在成像方向和视场范围方面难以灵活应用于不同场景。
发明内容
本发明所要解决的技术问题是:克服现有技术的不足,提供一种单相机球面镜反射成像投影装置,采用离轴非中心对称的单相机球面镜结构可以在获得更大视场范围的同时,有效提高应用灵活性;同时建立离轴非中心式球面镜反射成像模型,实现空间中三维点经球面镜反射至单相机成像面二维图像点的投影映射,适用于复杂场景及受限空间。
本发明的技术解决方案是:一种单相机球面镜反射成像投影装置,其特征在于:由空间点P(4)、球面反射镜(1)、摄像机(2)、数据传输线(8)和计算机(9)组成;
摄像机(2)与球面反射镜(1)离轴非中心对称放置,距离为200~500mm;
三维空间点P(4)处于摄像机(2)经球面反射镜(1)反射的视场范围内,三维空间点P(4)经球面反射镜(1)反射后进入摄像机(2)中,在摄像机(2)的二维图像平面(3)上成像,形成二维图像点(7),成像所得二维图像数据经数据传输线(8)存储至计算机(9)中。建立单相机球面镜反射成像投影模型为,通过三维空间点P(4)经球面反射镜(1)至摄像机(2)的二维图像平面(3)映射的数学关系,计算得到二维图像坐标,完成球面镜反射成像投影过程。模型通过仿真实验进行验证,设定仿真实验点为三维空间点阵列,阵列点三维坐标经模型映射后得到二维图像坐标。
所述建立单相机球面镜反射成像投影模型过程为:三维空间点P(4)经球面镜反射成像时,可以看成摄像机(2)对其经球面镜反射后所形成的镜像虚拟点P’(6)进行拍摄,具体步骤是:
1、建立球面镜反射成像投影模型坐标系;
1.1、在球面反射镜(1)上建立三维世界坐标系ow-xwywzw,原点ow为球面表面的中心点,xw轴和yw轴相互垂直,xw轴与球面反射镜(1)中心对称轴重合且指向球心;
1.2、建立摄像机坐标系oc-xcyczc,原点oc为摄像机光学中心,xc轴与摄像机镜头光轴重合且指向摄像机拍摄方向,在摄像机的二维图像平面(3)上建立二维图像坐标系ou-xuyu;
1.3、定义摄像机(2)的内部参数矩阵为A,位置和拍摄角度为oc-xcyczc坐标系到ow-xwywzw坐标系的平移矢量T和旋转矩阵R;
2、建立单相机球面镜反射成像投影模型,并与平面镜反射成像对比;
2.1、离轴三维空间点P(4)经球面反射镜(1)反射所成镜像虚拟点为P’(6),经相同位置的平面镜反射所成镜像虚拟点为P”(5),由平面镜像虚拟点P”(5)作为真值对比;
2.2、球面反射镜像虚拟点P’(6)由光路反向追迹获得,三维空间点向球面镜球心正向入射的光线沿原路返回,反射光线的反向延长线为lPC,平行于球面镜中心对称轴入射的光线反向延长线lPF过焦点F,lPC和lPF的交点即为球面镜像虚拟点P’(6);平面镜像虚拟点P”(5)与空间点关于yw轴成对称分布;
2.3、根据三维空间点P(4)、球心C和焦点F的齐次坐标计算上述两条直线的向量表示,由二维平面内直线相交计算球面镜像虚拟点P’(6)的齐次坐标,由三维空间点P(4)镜像对称计算平面镜像虚拟点P”(5)的齐次坐标,并将二者扩展为ow-xwywzw坐标系下的三维坐标齐次表示;
2.4、根据R,T将球面镜像虚拟点P’(6)和平面镜像虚拟点P”(5)的三维坐标转换至摄像机坐标系下,并由摄像机内部参数进行透视投影变换,计算图像点(7)和平面镜像虚拟点P”(5)对应图像点的二维图像坐标,完成球面镜反射成像投影过程。
使用如上面所述的单相机球面镜反射成像投影模型进行仿真实验,验证其成像效果,其特征在于,由平面镜像虚拟点作为真值对比,具体步骤如下:
1、仿真实验初始参数设定;
1.1、设定单相机球面镜反射成像投影装置的参数,包括设定球面镜曲率半径r,摄像机的内部参数矩阵A;
1.2、设定仿真三维空间点阵列的初始三维坐标,使三维空间点阵列在ow-xwywzw坐标系下zw轴坐标为0;
2、仿真三维空间点阵列镜像虚拟点计算;
2.1、计算球面镜像虚拟点P’(6)和平面镜像虚拟点P”(5);
2.2、将球面镜像虚拟点坐标扩展为三维齐次坐标,同时将平面镜像虚拟点扩展为三维齐次坐标;
3、镜像虚拟点经摄像机进行透视投影计算;
3.1、设定摄像机的位置和拍摄角度,即旋转矩阵R和平移矢量T;
3.2、将球面镜像虚拟点三维坐标转换至摄像机坐标系下,并同时将平面镜像虚拟点三维坐标转换至摄像机坐标系下;
3.3、设定仿真图像分辨率,并将球面和平面镜像三维虚拟点投影至仿真图像上,得到球面虚拟图像点和平面虚拟图像点;
3.4、判断是否大部分球面和平面镜像虚点的图像点均在仿真图像范围内;若是,则进行下一步;若否,则调整仿真点的位置和阵列间隔,重新设定仿真点阵列初始坐标,并回到步骤3.1.2重新进行;
3.5、判断是否有足够的位姿观察规律;若是,则进行下一步;若否,则调整摄像机的位置和拍摄角度,重新设定oc-xcyczc坐标系到ow-xwywzw坐标系的转换关系,并回到步骤3.3.1重新进行;
3.6、绘制球面和平面镜像虚点分别在球面镜坐标系和摄像机坐标系下的三维图,绘制平面和球面虚拟图像点二维分布图;
3.7、判断是否有不同曲率半径的实验结果;若是,则进行下一步;若否,则调整球面镜曲率,重新设定球面镜曲率半径,并回到步骤3.1.1重新进行,重复1.1~3.7步,进行下一次实验;
3.8、对比不同曲率半径及不同位姿下的三维图和二维图,分析成像规律,归纳球面镜反射成像机理。
本发明与现有技术相比的优点在于:
(1)本发明包括单摄像机球面反射镜,摄像机与球面镜离轴非中心对称放置,实施过程分为球面镜反射成像投影模型建立和仿真实验。经初始参数设定后实施仿真实验验证。模型建立时,由摄像机对球面镜像虚拟点拍摄,仿真实验中仿真点为三维空间中的二维阵列点。本发明的模型解决了空间三维点经球面镜反射至摄像机图像平面二维点的映射关系,实现了球面镜反射成像投影过程的定量数学表达,有效增大了折反射相机的视场范围,扩展了适用性;
(2)本发明装置中提出了离轴非中心式的单相机球面镜反射成像,实现了摄像机从任意角度经球面镜反射对空间场景成像,扩大了折反射相机的视场范围,并缩小了系统体积;
(3)建立了一种通用单相机球面镜反射成像投影模型,实现了空间点经球面镜反射至摄像机图像平面的映射,同时以平面反射镜作为基准,可用于后续对球面反射成像畸变的校正;
(4)本发明的单相机球面镜反射成像装置结构简单、成本低、适用性强,不受共轴安装限制,扩展了折反射相机的应用范围,可用于受限空间及复杂场景的成像及测量。
附图说明
图1为本发明的单相机球面镜反射成像装置示意图;
图2为本发明的单相机球面镜反射成像模型示意图;
图3为本发明中单相机球面镜反射成像模型仿真实验流程图;
图4为仿真实验空间点阵列及其球面和平面镜像虚拟点分布图;
图5为仿真实验空间点经球面镜和平面镜反射至摄像机成像面的图像点。
具体实施方式
下面对本发明做进一步详细说明。本发明基于多视几何和图像处理技术,建立单相机球面镜反射成像模型,实现三维空间点经球面镜反射至二维图像平面的映射。
如图1所示,本发明的单相机球面镜反射成像装置,由三维空间点4即点P、球面反射镜1、摄像机2、数据传输线8和计算机9组成;摄像机2与球面反射镜1离轴非中心对称放置,距离为200~500mm;三维空间点P处于摄像机2经球面反射镜1反射的视场范围内,三维空间点P经球面反射镜1反射后,进入摄像机2内,在摄像机图像平面上成像7;拍摄所得图像数据经数据传输线8存储至计算机9进行处理。
单相机球面镜反射成像投影模型为,三维空间点P经球面反射镜1至摄像机2的二维图像平面3映射的数学关系,模型通过仿真实验进行验证,设定仿真实验点为三维空间点阵列,阵列点三维坐标经模型映射后得到二维图像坐标。
本发明球面镜反射成像投影模型用于定量描述三维空间点经球面镜反射至二维图像平面点的映射关系,可用于球面反射成像畸变校正。
本发明的单相机球面镜反射成像投影模型的数学原理为:物体经球面反射镜反射成像时,可以看成对真实空间被测点经球面镜反射后的镜像虚拟点的观察,虚拟点在摄像机平面内的成像可由摄像机模型透视投影变换实现。
如图2所示,在球面反射镜1上建立三维世界坐标系ow-xwywzw,原点ow为球面表面的中心点,xw轴和yw轴相互垂直,xw轴与球面镜中心对称轴重合且指向球心,建立摄像机坐标系oc-xcyczc,原点oc为摄像机光学中心,xc轴与摄像机镜头光轴重合且指向摄像机拍摄方向,在摄像机的二维图像平面3上建立二维图像坐标系ou-xuyu,定义摄像机2的内部参数矩阵为A,位置和拍摄角度为oc-xcyczc坐标系到ow-xwywzw坐标系的平移矢量T和旋转矩阵R。
建立球面镜反射成像投影模型,并与平面镜反射成像对比。采用几何光路计算的方法,推导球面镜反射成像公式。单球面反射成像公式为:
焦点为F,球心为C,焦距与半径存在二倍关系,即f=r/2。离轴三维空间点P经球面镜所成虚像点位置的示意图,利用光路反向追迹分析,三维空间点直接向C点入射的光线,沿原路返回,平行于图中x轴入射的光线,射线反向延长线过焦点F,由此两条反射延长线的交点P’即为三维空间点P的虚像。坐标原点建立在球面与轴线的交点处,则P’点坐标(a,b),P点二维齐次坐标为xP=[a,b,1]T,经过O点的平面镜反射所成虚像P’u点的坐标为(-a,-b),球面半径为球心C点齐次坐标为xC=[r,0,1]T,则直线PC的向量表示为:
P点经球面镜反射所呈的虚像点P’为直线PC和直线AF的交点,则P’的齐次坐标为:
上式为P’点在球面镜坐标系下的二维齐次坐标,若扩展至三维,增加一z轴,则xp’的齐次表示为:
P”在球面镜坐标系下的三维坐标与P点关于yw轴成对称分布,将P”点坐标扩展至三维,则xp”的齐次表示为:
对三维空间点的球面和平面镜像虚拟点进行透视投影映射:
x=A[R T]X (7)
其中,X表示三维空间点对应镜像虚拟点的三维坐标,x表示镜像虚拟点经摄像机透视投影模型在二维图像平面上的像素坐标。
根据的上述单相机球面镜反射成像投影模型和本发明的基本原理,进行球面镜反射成像前向投影仿真实验,即将三维空间点经球面镜反射后的虚点,映射至摄像机的二维图像平面上,得到图像点坐标并观察其分布情况。
如图3所示,前向投影仿真实验的具体步骤如下:
1、设定单相机球面镜反射成像投影装置的参数,包括设定球面镜曲率半径r,摄像机的内部参数矩阵A;
2、设定仿真三维空间点的初始三维坐标,使三维空间点在ow-xwywzw坐标系下zw轴坐标为0,给定a,b值,作为仿真三维空间点P在球面镜下真实点的三维坐标(a,b,0)。
3、根据球面镜反射成像映射模型,计算球面镜像虚拟点P’(a’,b’,0),同时计算对应的平面镜镜像虚拟点P”(-a,b,0),作为对比;
5、设定摄像机的位置和拍摄角度,即oc-xcyczc坐标系到ow-xwywzw坐标系的转换关系R,T;
6、将球面镜像虚拟点三维坐标转换至摄像机坐标系下,得到XCP’,并同时将平面镜像虚拟点三维坐标转换至摄像机坐标系下,得到XCP”;
7、设定仿真图像分辨率L×H,并将球面和平面镜像三维虚拟点投影至仿真图像上,得到球面虚拟图像点xP’(u,v)和平面虚拟图像点xP”(up,vp);
8、判断是否大部分球面和平面镜像虚拟点的图像点均在仿真图像范围内;若是,则进行下一步;若否,则调整仿真点的位置和阵列间隔,重新设定仿真点阵列初始坐标,并回到步骤2重新进行;
9、判断是否有足够的位姿观察规律;若是,则进行下一步;若否,则调整摄像机的位置和拍摄角度,重新设定oc-xcyczc坐标系到ow-xwywzw坐标系的转换关系,并回到步骤5重新进行;
10、绘制球面和平面镜像虚拟点分别在球面镜坐标系和摄像机坐标系下的三维图,绘制平面和球面虚拟图像点二维分布图;
11、判断是否有不同曲率半径的实验结果;若是,则进行下一步;若否,则调整球面镜曲率,重新设定球面镜曲率半径,并回到步骤1重新进行,重复1~10步,进行下一次实验;
12、对比不同曲率半径及不同位姿下的三维图和二维图,分析成像规律,归纳球面镜反射成像机理。
实施例
根据本发明提出成像装置和反射投影模型,设计了前向投影仿真实验,对三维空间点经球面镜反射至摄像机成像的过程进行了测试。
实验设备为DELL台式计算机OPTIPLEX 7020,实验平台为Matlab R2019a,实验参数设置:球面镜曲率半径r=1000mm;摄像机内部参数矩阵为:
球面镜坐标系到摄像机坐标系的转换关系:
(旋转矩阵可分解为三个方向的旋转向量,旋转角度为15°、75°、-30°,M矩阵第四列为平移矢量)
参与球面镜反射成像映射模型仿真实验的一组阵列点数量为:7×7=49个。设置7×7仿真点阵列在球面镜坐标系下的坐标a的范围为[-180mm,-60mm],间隔为20mm,b的范围为[20mm,140mm],间隔为20mm。
根据摄像机内参设置,图像分辨率大小为2400×1800像素,仿真实验设定摄像机不变,则内参矩阵不变。图4为r=1000mm时仿真实验三维空间点阵列及其平面和球面反射虚点的二维显示,其中,圆点阵列为初始仿真点,方形阵列为平面镜像虚拟点,菱形阵列为球面镜像虚拟点,图5为在此位姿下仿真实验平面和球面反射虚点的图像点,其中,圆点阵列为球面镜反射成像点,方形阵列为平面镜反射成像点,箭头为球面镜反射产生的图像畸变。由图4可以看出三维空间点经球面镜反射后形成的镜像虚拟点阵列比平面镜反射镜像虚拟点阵列的分布更加聚集,所占的空间范围显著减小,从而使摄像机经球面镜反射捕捉的空间场景范围更大;对应在图5中,可以看到球面反射镜像虚拟点的二维图像点阵列比平面反射镜像虚拟点的二维图像点阵列的分布更加聚集,所占的图像范围显著减小,因此本发明采用的单相机球面镜反射成像装置能够有效扩大视场范围,阐释了球面镜反射成像机理,同时所提出的成像投影模型实现了三维空间点经球面镜反射至二维图像平面的映射。
Claims (2)
1.一种单相机球面镜反射成像投影装置,其特征在于,包括:球面反射镜(1)、摄像机(2)和计算机(9);
摄像机(2)与球面反射镜(1)离轴非中心对称放置,距离为200~500mm;
三维空间点P(4)处于摄像机(2)经球面反射镜(1)反射的视场范围内,三维空间点P(4)经球面反射镜(1)反射后进入摄像机(2)中,在摄像机(2)的二维图像平面(3)上成像,形成二维图像点(7),成像所得二维图像数据经数据传输线(8)存储至计算机(9)中,单相机球面镜反射成像投影模型为通过三维空间点P经球面反射镜(1)至摄像机(2)的二维图像平面(3)映射的数学关系,计算得到二维图像坐标,完成球面镜反射成像投影过程;模型通过仿真实验进行验证,设定仿真实验点为三维空间点阵列,阵列点三维坐标经模型映射后得到二维图像坐标;
所述单相机球面镜反射成像投影模型的建立为:三维空间点P(4)经球面镜反射成像时,看成摄像机(2)对其经球面镜反射后所形成的球面镜像虚拟点P’(6)进行拍摄,具体步骤如下:
2.1、建立球面镜反射成像投影模型坐标系;
2.1.1、在球面反射镜(1)上建立三维世界坐标系ow-xwywzw,原点ow为球面表面的中心点,xw轴和yw轴相互垂直,xw轴与球面反射镜(1)中心对称轴重合且指向球心;
2.1.2、建立摄像机坐标系oc-xcyczc,原点oc为摄像机光学中心,xc轴与摄像机镜头光轴重合且指向摄像机拍摄方向,在摄像机的二维图像平面(3)上建立二维图像坐标系ou-xuyu;
2.1.3、定义摄像机(2)的内部参数矩阵为A,位置和拍摄角度为oc-xcyczc坐标系到ow-xwywzw坐标系的平移矢量T和旋转矩阵R;
2.2、建立单相机球面镜反射成像投影模型,并与平面镜反射成像对比;
2.2.1、三维空间点P(4)经球面反射镜(1)反射所成球面镜像虚拟点为P’(6),经相同位置的平面镜反射所成平面镜像虚拟点为P”(5),由平面镜像虚拟点P”(5)作为真值对比;
2.2.2、球面反射镜像虚拟点P’(6)由光路反向追迹获得,三维空间点向球面镜球心正向入射的光线沿原路返回,反射光线的反向延长线为lPC,平行于球面反射镜(1)中心对称轴入射的光线反向延长线lPF过焦点F,lPC和lPF的交点即为球面镜像虚拟点P’(6);平面镜像虚拟点P”(5)与空间点关于yw轴成对称分布;
2.2.3、根据三维空间点P(4)、球心C和焦点F的齐次坐标计算lPC和lPF的向量表示,由二维平面内直线相交计算球面镜像虚拟点P’(6)的齐次坐标,由三维空间点P(4)镜像对称计算平面镜像虚拟点P”(5)的齐次坐标,并将二者扩展为ow-xwywzw坐标系下的三维坐标齐次表示;
2.2.4、根据平移矢量T和旋转矩阵R将球面镜像虚拟点P’(6)和平面镜像虚拟点P”(5)的三维坐标转换至摄像机坐标系下,并由摄像机内部参数进行透视投影变换,计算二维图像点(7)和平面镜像虚拟点P”(5)对应图像点的二维图像坐标,完成球面镜反射成像投影过程。
2.根据权利要求1所述的单相机球面镜反射成像投影装置,其特征在于:所述步骤2.2.1具体步骤如下:
3.1、仿真实验初始参数设定;
3.1.1、设定单相机球面镜反射成像投影装置的参数,包括设定球面镜曲率半径r,摄像机的内部参数矩阵A;
3.1.2、设定仿真三维空间点阵列的初始三维坐标,使三维空间点在ow-xwywzw坐标系下zw轴坐标为0;
3.2、仿真三维空间点阵列镜像虚拟点计算;
3.2.1、计算球面镜像虚拟点P’(6)和平面镜像虚拟点P”(5);
3.2.2、将球面镜像虚拟点坐标扩展为三维齐次坐标,同时将平面镜像虚拟点扩展为三维齐次坐标;
3.3、镜像虚拟点经摄像机进行透视投影计算;
3.3.1、设定摄像机的位置和拍摄角度,即旋转矩阵R和平移矢量T;
3.3.2、将球面镜像虚拟点三维坐标转换至摄像机坐标系下,并同时将平面镜像虚拟点三维坐标转换至摄像机坐标系下;
3.3.3、设定仿真图像分辨率,并将球面和平面镜像三维虚拟点投影至仿真图像上,得到球面虚拟图像点和平面虚拟图像点;
3.3.4、判断是否大部分球面和平面镜像虚拟点的图像点均在仿真图像范围内;若是,则进行下一步;若否,则调整仿真点的位置和阵列间隔,重新设定仿真点阵列初始坐标,并回到步骤3.1.2重新进行;
3.3.5、判断是否有足够的位姿观察规律;若是,则进行下一步;若否,则调整摄像机的位置和拍摄角度,重新设定oc-xcyczc坐标系到ow-xwywzw坐标系的转换关系,并回到步骤3.3.1重新进行;
3.3.6、绘制球面镜像虚拟点和平面镜像虚拟点分别在球面镜坐标系和摄像机坐标系下的三维图,绘制平面虚拟图像点和球面虚拟图像点的二维分布图;
3.3.7、判断是否有不同曲率半径的实验结果;若是,则进行下一步;若否,则调整球面镜曲率,重新设定球面镜曲率半径,并回到步骤3.1.1重新进行,重复3.1.1~3.3.7步,进行下一次实验;
3.3.8、对比不同曲率半径及不同位姿下的三维图和二维图,分析成像规律,归纳球面镜反射成像机理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111246846.XA CN113989105B (zh) | 2021-10-26 | 2021-10-26 | 一种单相机球面镜反射成像投影装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111246846.XA CN113989105B (zh) | 2021-10-26 | 2021-10-26 | 一种单相机球面镜反射成像投影装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113989105A CN113989105A (zh) | 2022-01-28 |
CN113989105B true CN113989105B (zh) | 2022-08-16 |
Family
ID=79741523
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111246846.XA Active CN113989105B (zh) | 2021-10-26 | 2021-10-26 | 一种单相机球面镜反射成像投影装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113989105B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115060166B (zh) * | 2022-08-18 | 2022-11-08 | 中国石油大学(华东) | 一种基于反射棱镜的视觉测量系统三维光路分析方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102788559A (zh) * | 2012-07-19 | 2012-11-21 | 北京航空航天大学 | 一种宽视场结构光视觉测量系统及测量方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006011103A (ja) * | 2004-06-28 | 2006-01-12 | Yokogawa Electric Corp | 球面ミラー撮像装置 |
US20070247729A1 (en) * | 2006-04-25 | 2007-10-25 | Rudolph Technologies, Inc. | Reflective objective |
US7768527B2 (en) * | 2006-05-31 | 2010-08-03 | Beihang University | Hardware-in-the-loop simulation system and method for computer vision |
CN101419061B (zh) * | 2008-12-08 | 2011-06-29 | 北京航空航天大学 | 一种镜像式结构光视觉测量系统和测量方法 |
US8493432B2 (en) * | 2010-06-29 | 2013-07-23 | Mitsubishi Electric Research Laboratories, Inc. | Digital refocusing for wide-angle images using axial-cone cameras |
CN102589425B (zh) * | 2012-01-12 | 2014-01-15 | 中国科学院国家天文台 | 球面射电望远镜反射面单元面板测量方法 |
CN105321181A (zh) * | 2015-10-28 | 2016-02-10 | 云南大学 | 使用双球的相离像与圆环点的像标定拋物折反射摄像机 |
CN112985293B (zh) * | 2021-01-14 | 2023-05-30 | 北京航空航天大学 | 一种单像机双球面镜镜像双目视觉测量系统和测量方法 |
-
2021
- 2021-10-26 CN CN202111246846.XA patent/CN113989105B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102788559A (zh) * | 2012-07-19 | 2012-11-21 | 北京航空航天大学 | 一种宽视场结构光视觉测量系统及测量方法 |
Non-Patent Citations (1)
Title |
---|
基于球面透视投影约束的全景环形透镜畸变校正;肖潇 等;《光学学报》;20080415;第28卷(第04期);第675-680页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113989105A (zh) | 2022-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI555379B (zh) | 一種全景魚眼相機影像校正、合成與景深重建方法與其系統 | |
CN106570938A (zh) | 基于opengl 的全景监控方法及系统 | |
CN107563088A (zh) | 一种基于光线跟踪算法的光场显示设备仿真方法 | |
CN114283203B (zh) | 一种多相机系统的标定方法及系统 | |
CN111922510B (zh) | 一种激光可视化加工方法及系统 | |
CN111464722A (zh) | 一种基于曲面复眼阵列的大角度光场成像装置及方法 | |
CN105374067A (zh) | 一种基于pal相机的三维重建方法及其重建系统 | |
CN109799073A (zh) | 一种光学畸变测量装置及方法、图像处理系统、电子设备和显示设备 | |
CN112489109B (zh) | 一种三维成像系统方法、装置及三维成像系统 | |
Ge et al. | Mirror-assisted multiview DIC for 360° panoramic large deformation measurement | |
CN113989105B (zh) | 一种单相机球面镜反射成像投影装置 | |
JP7295641B2 (ja) | 水平視差ステレオパノラマ取込方法 | |
KR102478989B1 (ko) | 명시야 데이터와 연관된 보케를 표현하는 데이터를 생성하기 위한 방법 및 장치 | |
CN111739103A (zh) | 一种基于单点标定物的多相机标定系统 | |
CN116147534A (zh) | 一种基于镜面辅助的多视角三维激光扫描系统及复杂表面全景测量方法 | |
CN115359127A (zh) | 一种适用于多层介质环境下的偏振相机阵列标定方法 | |
CN117115272A (zh) | 用于降水粒子多角度成像的远心相机标定和三维重构方法 | |
RU2729698C2 (ru) | Устройство и способ для кодирования изображения, захваченного оптической системой получения данных | |
TW565736B (en) | Method for determining the optical parameters of a camera | |
CN114241059B (zh) | 一种光度立体视觉系统中相机和光源的同步标定方法 | |
CN115439541A (zh) | 用于折射成像系统的玻璃方位标定系统及方法 | |
CN109345595A (zh) | 一种基于圆球透镜的立体视觉传感器标定方法 | |
CN112361989B (zh) | 一种通过点云均匀性考量测量系统标定参数的方法 | |
CN111947598B (zh) | 基于平面反射镜的360°三维人头测量方法 | |
KR20210145315A (ko) | 픽셀 빔을 나타내는 데이터를 생성하기 위한 장치 및 방법 |
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 |