CN117788564B - 一种球心定位方法和球心定位装置 - Google Patents
一种球心定位方法和球心定位装置 Download PDFInfo
- Publication number
- CN117788564B CN117788564B CN202311798029.4A CN202311798029A CN117788564B CN 117788564 B CN117788564 B CN 117788564B CN 202311798029 A CN202311798029 A CN 202311798029A CN 117788564 B CN117788564 B CN 117788564B
- Authority
- CN
- China
- Prior art keywords
- camera
- sphere
- cameras
- coordinate system
- center
- 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 19
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 238000002372 labelling Methods 0.000 claims description 17
- 230000003287 optical effect Effects 0.000 claims description 10
- 238000003384 imaging method Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 238000006467 substitution reaction Methods 0.000 claims 2
- 230000006870 function Effects 0.000 description 18
- 238000005457 optimization Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及一种球心定位方法和球心定位装置,包括:将至少三个相机组合成相机阵列,每个相机以不同角度对球体进行拍摄,以获得至少三张照片,对所述照片中球体图案边缘的像素点进行标注,并将标注的所述像素点和所述相机进行对应;获取每个相机的内参矩阵和外参矩阵;针对第m个相机对应的第n个标注的像素点构建目标损失函数Tmn,Tmn基于第m个相机的内参矩阵和外参矩阵,以及第n个标注的像素点的像素坐标进行构建;基于所有Tmn获得相机阵列的目标损失函数T;采用梯度下降算法获取T的极值,以确定球体的球心在世界坐标系中的坐标。
Description
技术领域
本发明涉及球心定位领域,特别是涉及一种球心定位方法和球心定位装置。
背景技术
在光源标定、4D扫描渲染等领域,需要在三维空间中多个视角下获取小球位置,从而确定球心坐标。现有技术通常是使用特有的标定设备确定球心位置,这种设备不仅成本高昂,而且实际标定精度也相对不足。
发明内容
基于此,有必要针对球心位置定位成本高、精度不足的的问题,提供一种球心定位方法和球心定位装置。
一种球心定位方法,包括:
将至少三个相机组合成相机阵列,每个相机以不同角度对球体进行拍摄,以获得至少三张照片,对所述照片中球体图案边缘的像素点进行标注,并将标注的所述像素点和所述相机进行对应;
获取每个相机的内参矩阵和外参矩阵;
针对第m个相机对应的第n个标注的像素点构建目标损失函数Tmn,Tmn基于第m个相机的内参矩阵和外参矩阵,以及第n个标注的像素点的像素坐标进行构建;
基于所有Tmn获得相机阵列的目标损失函数T;
采用梯度下降算法获取T的极值,以确定球体的球心在世界坐标系中的坐标。
本发明T为所有Tmn加权求和,每个Tmn的权重系数均相等,且每个相机对应的像素点数量相同。
本发明每个相机对应的像素点数量至少有两个。
本发明第m个相机对应的内参矩阵其中fm为第m个相机以像素为单位的焦距,u0m和v0m分别为第m个相机的焦点在成像平面中对应像素的横坐标和纵坐标,第m个相机具有两个所述外参矩阵,分别为Rm和Tm,其中Rm为第m个相机相对世界坐标系的旋转矩阵,Tm为第m个相机的光心相对世界坐标系的平移矩阵。
本发明获取第m个相机对应的的第n个标注的像素点在成像平面中的坐标(umn,vmn),
基于(umn,vmn)获得sm=(umn-u0m,vmn-v0m,fm);
基于sm获取(Lm,Mm,Nm),其中(Lm,Mm,Nm)=Rm×sm;
基于Tm获取第m个相机光心在世界坐标系下的坐标(Xm,Ym,Zm);
则Tmn=(x-Xm)2+(y-Ym)2+(z-Zm)2-r2-[(x-Xm)*Lm+(y-Ym)*Mm+(z-Zm)*Nm]/(Lm 2+Mm 2+Nm 2);
T为极值时对应的(x,y,z)为球心在世界坐标系下的坐标。
本发明将(x,y,z)代入T,以获取T为极值时对应r的数值解,以此获得球体的半径。
一种球心定位装置,包括:
相机阵列,所述相机阵列包括至少三个相机,所述相机用于拍摄球体以获得照片;
参数获取模块,用于获取照片中像素点的像素坐标以及相机的内参矩阵和外参矩阵;
损失函数构建模块,用于构建相机阵列的目标损失函数;
计算模块,用于对相机阵列的目标损失函数执行梯度下降算法,以获得球体的球心在世界坐标系中的坐标。
本发明所述计算模块还用于计算得到球体的半径。
本发明所述球心定位装置还包括标注模块,所述标注模块用于标注所述照片中球体图案边缘的像素点。
本发明的有益效果为:
本发明利用至少三个相机对球体进行拍摄,进而获得球体在三个不同角度下的照片,结合相机自身的参数矩阵,以此获得相机阵列的目标损失函数,对应的,该目标损失函数包含球体的三维信息,基于梯度下降算法,即可计算得到球体的球心坐标。只要增加相机和/或标注的像素点数量,就能有效提升球心坐标的计算精度,该精度提升手段成本低廉,且提升效果较为明显。
本发明相较现有技术而言,鲁棒性高,对硬件要求低,能够以较低的成本获得更为精确的球心位置信息。
附图说明
图1为本发明实施例1中球心定位方法的步骤流程图;
图2为本发明实施例1中球心定位装置的拓扑图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
实施例1:
参见图1,本实施例提供了一种球心定位方法,包括如下步骤:
步骤101:布设至少三个相机,所有相机呈环形阵列排布,以组成201,待确定球心位置的球体位于201的中心,由此每个相机能够以不同的角度对球体进行拍摄。每个相机分别对应一张照片,故而照片总数至少有三张,每张照片对应不同角度的球体,故而三张照片组合后,理论上包含了球体的完整三维信息。
步骤102:获取每个相机的内参矩阵和外参矩阵,每个相机的内参矩阵均为一个,外参矩阵均为两个。
第m个相机对应的内参矩阵其中fm为第m个相机以像素为单位的焦距,u0m和v0m分别为第m个相机的焦点在成像平面中对应像素的横坐标和纵坐标。
第m个相机具有两个所述外参矩阵,分别为Rm和Tm,其中Rm为第m个相机相对世界坐标系的旋转矩阵,Tm为第m个相机的光心相对世界坐标系的平移矩阵。
步骤103:对每张照片中球体图案边缘的像素点进行标注,每张照片中像素点的标注数量可以为一个,也可以是两个以上。标注的像素点通过其所在的照片和相机进行对应,由此每个相机都对应至少一个标注的像素点。
像素点的标注方式不限,可以采用AI标注,也可以采用人工标注。
步骤104:针对第m个相机对应的第n个标注的像素点构建目标损失函数Tmn,进而获得201的目标损失函数T。
步骤104具体包括步骤1041-步骤1046。
步骤1041:获取第m个相机对应的的第n个标注的像素点在成像平面中的坐标(umn,vmn)。
步骤1042:基于(umn,vmn)获得sm=(umn-u0m,vmn-v0m,fm),其中sm代表了经过第m个相机的光心并与球体相切的其中一个光线的方向。
步骤1043:基于sm获取(Lm,Mm,Nm),其中(Lm,Mm,Nm)=Rm×sm。
步骤1044:基于Tm获取第m个相机光心在世界坐标系下的坐标(Xm,Ym,Zm)。
步骤1045:构建第m个相机对应的第n个标注的像素点对应的目标损失函数Tmn,其中
Tmn=(x-Xm)2+(y-Ym)2+(z-Zm)2-r2-[(x-Xm)*Lm+(y-Ym)*Mm+(z-Zm)*Nm]/(Lm 2+Mm 2+Nm 2)。
步骤1046:对所有Tmn进行加权求和,以获得201的目标损失函数T。
步骤105:对T执行梯度下降算法,由此参数r被消去,并获得T的极值(尽可能为零极值点),T为极值时对应的(x,y,z)即球体的球心在世界坐标系中的坐标。
步骤106:将T为极值时对应的(x,y,z)代入T,结合T的极值,可以得到参数r的具体数值,此时r的具体数值即对应球体半径。
即,本实施例球心定位方法除了能够求解获得球心在世界坐标系中的三维坐标,即三个变量,还能额外求解得到球体半径。
相机仅为三个的情况下,依然能够获得球体的三维信息,由此求解得到球心在世界坐标系中的三维坐标并不难理解。而在此条件下,却能够求解得到球体半径这一第四维信息,明显有违于本领域技术人员的一般常识。而实现这一结果的原因在于本实施例中针对Tmn的构建方式,其二次型函数的表达使得T的优化问题转换成了凸优化问题,使得球体半径真实数值以及球心位置真实坐标对应在了T的极值点,且球体半径对于的参数r在求导过程中被消去,使得球体半径真实数值以及球心位置真实坐标的获取过程相互独立,分先后进行。
但是值得注意的是,在本实施例中Tmn的表达方式下,针对多个m相同而n不同的Tmn而言,求导后得到的实际上是仅为常数项不同的同一类函数,无法得到极值信息。故而为了获取T的极值,就必须依赖至少三个不同m数值的Tmn进行T的构建,即相机最少需要三个。
由于照片成像过程中存在径向和切向畸变,像素点标注存在个人主观意识等问题,无法避免的回引入噪声到测量中,故理论上更多的相机、更多标注的像素点可以提升球体半径真实数值以及球心位置真实坐标的计算精度。
参见图2,在一个实施例中,提供了一种球心定位装置,包括201、参数获取模块203、标注模块202、损失函数构建模块204和计算模块205。
其中201包括至少三个相机,相机用于拍摄球体以获得照片。
标注模块202用于标注照片中球体图案边缘的像素点。
参数获取模块203用于获取照片中像素点的像素坐标以及相机的内参矩阵和外参矩阵。
损失函数构建模块204用于构建201的目标损失函数;
计算模块205用于对201的目标损失函数执行梯度下降算法,以获得球体的球心在世界坐标系中的坐标,并将球心在世界坐标系中的坐标代入201的目标损失函数,以获得球体的半径。
实施例2:
本实施例与实施例1的区别在于,像素点的标注依赖于人工进行。容易理解的是,像素点的标注不论采用人工还是AI,都会产生误差,但是人工标注的的误差可能会更大。
基于此,实施例1中采用AI标注的情况下,可以直接使用标注结果的置信度作为Tmn的权重系数,进而对所有Tmn进行加权求和得到T。
而在本实施例中,每个Tmn的权重系数均相等,由此每个相机对应的像素点数量需要保持一致,以确保所有相机的权重系数之和相等,以此避免在进行梯度下降算法时在T的数值中引入额外的不必要误差。
应该理解的是,除非本文中有明确的说明,上述步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,各个步骤可以包括多个步骤或者多个阶段,这些步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-Only Memory,ROM)、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(Random Access Memory,RAM)或外部高速缓冲存储器。作为说明而非局限,RAM可以是多种形式,比如静态随机存取存储器(Static Random Access Memory,SRAM)或动态随机存取存储器(Dynamic Random Access Memory,DRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (6)
1.一种球心定位方法,其特征在于,包括:
将至少三个相机组合成相机阵列,每个相机以不同角度对球体进行拍摄,以获得至少三张照片,对所述照片中球体图案边缘的像素点进行标注,并将标注的所述像素点和所述相机进行对应;
获取每个相机的内参矩阵和外参矩阵,第个相机对应的内参矩阵 ,其中为第个相机以像素为单位的焦距,和分别为第个相机的焦点在成像平面中对应像素的横坐标和纵坐标,第个相机具有两个所述外参矩阵,分别为和,其中为第个相机相对世界坐标系的旋转矩阵,为第个相机的光心相对世界坐标系的平移矩阵;
针对第个相机对应的第个标注的像素点构建目标损失函数,基于第个相机的内参矩阵和外参矩阵,以及第个标注的像素点的像素坐标进行构建,获取第个相机对应的的第个标注的像素点在成像平面中的坐标,基于获得;基于获取,其中=;基于获取第个相机光心在世界坐标系下的坐标;
则;
基于所有获得相机阵列的目标损失函数,对所有进行加权求和,以获得目标损失函数;
采用梯度下降算法获取的极值,以确定球体的球心在世界坐标系中的坐标,为极值时对应的为球心在世界坐标系下的坐标;将代入,以获取为极值时对应的数值解,以此获得球体的半径。
2.根据权利要求1所述的球心定位方法,其特征在于,为所有加权求和,每个的权重系数均相等,且每个相机对应的像素点数量相同。
3.根据权利要求2所述的球心定位方法,其特征在于,每个相机对应的像素点数量至少有两个。
4.根据权利要求1所述的球心定位方法,其特征在于,所述相机数量仅为三个。
5.一种球心定位装置,其特征在于,包括:
相机阵列,所述相机阵列包括至少三个相机,所述相机用于拍摄球体以获得照片;
参数获取模块,用于获取照片中像素点的像素坐标以及相机的内参矩阵和外参矩阵,第个相机对应的内参矩阵 ,其中为第个相机以像素为单位的焦距,和分别为第个相机的焦点在成像平面中对应像素的横坐标和纵坐标,第个相机具有两个所述外参矩阵,分别为和,其中为第个相机相对世界坐标系的旋转矩阵,为第个相机的光心相对世界坐标系的平移矩阵;
损失函数构建模块,用于构建相机阵列的目标损失函数,获取第个相机对应的的第个标注的像素点在成像平面中的坐标,基于获得;基于获取,其中=;基于获取第个相机光心在世界坐标系下的坐标;则,对所有进行加权求和,以获得目标损失函数;
计算模块,用于对相机阵列的目标损失函数执行梯度下降算法,以获得球体的球心在世界坐标系中的坐标,为极值时对应的为球心在世界坐标系下的坐标;将代入,以获取为极值时对应的数值解,以此获得球体的半径。
6.根据权利要求5所述的球心定位装置,其特征在于,所述球心定位装置还包括标注模块,所述标注模块用于标注所述照片中球体图案边缘的像素点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311798029.4A CN117788564B (zh) | 2023-12-25 | 2023-12-25 | 一种球心定位方法和球心定位装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311798029.4A CN117788564B (zh) | 2023-12-25 | 2023-12-25 | 一种球心定位方法和球心定位装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117788564A CN117788564A (zh) | 2024-03-29 |
CN117788564B true CN117788564B (zh) | 2024-09-13 |
Family
ID=90381025
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311798029.4A Active CN117788564B (zh) | 2023-12-25 | 2023-12-25 | 一种球心定位方法和球心定位装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117788564B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112066879A (zh) * | 2020-09-11 | 2020-12-11 | 哈尔滨工业大学 | 基于计算机视觉的气浮运动模拟器位姿测量装置及方法 |
CN113819890A (zh) * | 2021-06-04 | 2021-12-21 | 腾讯科技(深圳)有限公司 | 测距方法、装置、电子设备及存储介质 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112667837A (zh) * | 2019-10-16 | 2021-04-16 | 上海商汤临港智能科技有限公司 | 图像数据自动标注方法及装置 |
CN112184807B (zh) * | 2020-09-22 | 2023-10-03 | 深圳市衡泰信科技有限公司 | 一种高尔夫球落地式检测方法、系统及存储介质 |
CN112802121B (zh) * | 2021-01-14 | 2023-09-05 | 杭州海康威视数字技术股份有限公司 | 监控相机的标定方法 |
WO2022235709A1 (en) * | 2021-05-03 | 2022-11-10 | Ditto Technologies, Inc. | Generation of a 3d model of a reference object to perform scaling of a model of a user's head |
CN114820811A (zh) * | 2022-04-18 | 2022-07-29 | 东北大学 | 一种基于运动球体的同步相机的基础矩阵求解及标定方法 |
CN116091615A (zh) * | 2022-12-06 | 2023-05-09 | 北京工业大学 | 一种基于立体矩阵小球的rgbd相机坐标转换与视觉定位方法 |
CN116051634A (zh) * | 2022-12-28 | 2023-05-02 | 深圳优艾智合机器人科技有限公司 | 一种视觉定位方法、终端及存储介质 |
CN116012415A (zh) * | 2023-01-09 | 2023-04-25 | 上海创屹科技有限公司 | 基于视觉测量的乒乓球旋转估计方法及旋转估计系统 |
CN116563377A (zh) * | 2023-05-26 | 2023-08-08 | 北京邮电大学 | 一种基于半球投影模型的火星岩石测量方法 |
CN116778262B (zh) * | 2023-08-21 | 2023-11-10 | 江苏源驶科技有限公司 | 一种基于虚拟点云的三维目标检测方法和系统 |
-
2023
- 2023-12-25 CN CN202311798029.4A patent/CN117788564B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112066879A (zh) * | 2020-09-11 | 2020-12-11 | 哈尔滨工业大学 | 基于计算机视觉的气浮运动模拟器位姿测量装置及方法 |
CN113819890A (zh) * | 2021-06-04 | 2021-12-21 | 腾讯科技(深圳)有限公司 | 测距方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN117788564A (zh) | 2024-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113643378B (zh) | 多相机环境中主动式刚体的位姿定位方法及相关设备 | |
CN110969668A (zh) | 一种长焦双目相机的立体标定算法 | |
CN111612834B (zh) | 生成目标图像的方法、装置及设备 | |
US8842190B2 (en) | Method and apparatus for determining sensor format factors from image metadata | |
CN110717942A (zh) | 图像处理方法和装置、电子设备、计算机可读存储介质 | |
CN112562014B (zh) | 相机标定方法、系统、介质及装置 | |
CN111833237B (zh) | 基于卷积神经网络和局部单应性变换的图像配准方法 | |
CN106530358A (zh) | 仅用两幅场景图像标定ptz摄像机的方法 | |
CN110099267A (zh) | 梯形校正系统、方法以及投影仪 | |
CN111401387B (zh) | 异常样本构建方法、装置、计算机设备和存储介质 | |
CN113191954A (zh) | 一种基于双目相机的全景图像拼接方法 | |
WO2019176349A1 (ja) | 画像処理装置、および撮像装置、並びに画像処理方法 | |
CN108759788A (zh) | 无人机影像定位定姿方法及无人机 | |
Kurmi et al. | Pose error reduction for focus enhancement in thermal synthetic aperture visualization | |
CN112184809A (zh) | 相对位姿的估计方法、装置、电子设备及介质 | |
CN113963067B (zh) | 一种采用小靶标对大视场视觉传感器进行标定的标定方法 | |
CN107067441B (zh) | 摄像机标定方法及装置 | |
CN117788564B (zh) | 一种球心定位方法和球心定位装置 | |
CN113870350A (zh) | 图像样本标注方法、装置、终端设备和图像标注系统 | |
CN117038554B (zh) | 芯片定位方法和芯片转移方法 | |
CN113436267A (zh) | 视觉惯导标定方法、装置、计算机设备和存储介质 | |
CN116051634A (zh) | 一种视觉定位方法、终端及存储介质 | |
CN114897990A (zh) | 一种基于神经网络的相机畸变标定方法、系统和存储介质 | |
Xia et al. | A coarse-to-fine ghost removal scheme for HDR imaging | |
JP2020041950A (ja) | 測量装置、測量方法、及びプログラム |
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 |