CN111243035B - 相机标定方法、装置、电子设备和计算机可读存储介质 - Google Patents

相机标定方法、装置、电子设备和计算机可读存储介质 Download PDF

Info

Publication number
CN111243035B
CN111243035B CN202010352903.1A CN202010352903A CN111243035B CN 111243035 B CN111243035 B CN 111243035B CN 202010352903 A CN202010352903 A CN 202010352903A CN 111243035 B CN111243035 B CN 111243035B
Authority
CN
China
Prior art keywords
camera
images
image
moving images
calibration
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
Application number
CN202010352903.1A
Other languages
English (en)
Other versions
CN111243035A (zh
Inventor
袁睿
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Jouav Automation Technology Co ltd
Original Assignee
Chengdu Jouav Automation Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu Jouav Automation Technology Co ltd filed Critical Chengdu Jouav Automation Technology Co ltd
Priority to CN202010352903.1A priority Critical patent/CN111243035B/zh
Publication of CN111243035A publication Critical patent/CN111243035A/zh
Application granted granted Critical
Publication of CN111243035B publication Critical patent/CN111243035B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Studio Devices (AREA)

Abstract

本发明提供一种相机标定方法、装置、电子设备和计算机可读存储介质,涉及运动相机的图像处理领域。该相机标定方法包括:获取多张运动图像;多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像;获取多张运动图像中任意两张图像之间的基础矩阵;基础矩阵表征任意两张图像中基于纹理的图像特征点匹配关系;根据所有的基础矩阵,获取相机的标定内参数;标定内参数为相机拍摄图像的目标校正参数。使用具有纹理的运动图像,继而获取图像之间的基础矩阵并获取相机的标定内参数,不需依赖场景,增加了获取相机的标定内参数的通用性;标定内参数用于拍摄图像的目标校正,减少了图像在后处理过程中的误差。

Description

相机标定方法、装置、电子设备和计算机可读存储介质
技术领域
本发明涉及运动相机的图像处理领域,具体而言,涉及一种相机标定方法、装置、电子设备和计算机可读存储介质。
背景技术
近年来,由于计算机视觉相关技术的发展,基于图像的二维/三维场景获取技术已经在很多领域得到应用,如机器人路径导航、图像拼接、三维重建等领域。其中,相机的内参数,包括相机焦距和主点数据,是求解图像坐标到世界坐标的映射关系较为关键的部分。
现在得到广泛应用的相机内参获取方式,是通过获取一系列棋盘格图像,通过对图像使用“张正友标定”方法进行离线标定,得到相机的内参数。“张正友标定”方法依赖于固定标定物棋盘格,并且必须离线标定,限制了很多场景下对于相机内参的计算和获取。为了改善上述缺点,目前的技术方案利用应用场景中的特殊景物,比如监控摄像头的人行道或者其它平行景物,计算消影点,利用多组消影点对相机内参的约束,计算得到相机内参数。但是,依赖于应用场景中特殊元素的约束来完成相机自标定,对应用场景要求较高,且不具有通用性。
发明内容
有鉴于此,本发明的目的在于提供一种相机标定方法、装置、电子设备和计算机可读存储介质。
为了实现上述目的,本发明采用的技术方案如下:
第一方面,本发明提供一种相机标定方法,所述方法包括:获取多张运动图像;所述多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像。获取所述多张运动图像中任意两张图像之间的基础矩阵;所述基础矩阵表征所述任意两张图像中基于所述纹理的图像特征点匹配关系。根据所有的所述基础矩阵,获取所述相机的标定内参数;所述标定内参数为所述相机拍摄图像的目标校正参数。
在可选的实施方式中,所述多张运动图像中的每张运动图像具有序列标识,所述根据所有的所述基础矩阵,获取所述相机的标定内参数,包括:按照所述序列标识的顺序,依次解析所述每张运动图像对应的所述基础矩阵,以获取所述相机的线性内参;所述线性内参为所述相机拍摄图像的待确认校正参数。根据所有的所述基础矩阵,建立目标函数;所述目标函数用于确定所述多张运动图像处于预设约束条件下的目标校正参数。根据所述线性内参和所述目标函数,获取所述标定内参数。
在可选的实施方式中,所述根据所有的所述基础矩阵,建立目标函数,包括:确定每个所述基础矩阵对应的非线性优化函数;所述非线性优化函数用于修正所述线性内参;根据所有的所述非线性优化函数和所述预设约束条件,确定所述目标函数。
在可选的实施方式中,所述每张运动图像的序列标识包括所述每张运动图像的图像采集时间信息,所述获取所述多张运动图像中任意两张图像之间的基础矩阵,包括:确定所述多张运动图像中符合第一预设条件的至少三张关键图像;所述第一预设条件包括任意两张运动图像的重叠率大于或等于第一阈值,和/或所述任意两张运动图像的图像采集时间信息符合预设采集时间间隔;获取所述至少三张关键图像中的任意两张关键图像的基础矩阵。
在可选的实施方式中,获取所述至少三张关键图像中的任意两张关键图像的基础矩阵,包括:获取所述任意两张关键图像中的所述纹理符合第二预设条件的图像特征点;基于所述图像特征点的匹配关系确定所述基础矩阵。
第二方面,本发明提供一种相机标定装置,所述装置包括:获取模块和处理模块。所述获取模块,用于获取多张运动图像;所述多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像。所述处理模块,用于获取所述多张运动图像中任意两张图像之间的基础矩阵;所述基础矩阵表征所述任意两张图像中基于所述纹理的图像特征点匹配关系。所述处理模块还用于根据所有的所述基础矩阵,获取所述相机的标定内参数;所述标定内参数为所述相机拍摄图像的目标校正参数。
在可选的实施方式中,所述多张运动图像中的每张运动图像具有序列标识,所述处理模块还用于按照所述序列标识的顺序,依次解析所述每张运动图像对应的所述基础矩阵,以获取所述相机的线性内参;所述线性内参为所述相机拍摄图像的待确认校正参数。所述处理模块还用于根据所有的所述基础矩阵,建立目标函数;所述目标函数用于确定所述多张运动图像处于预设约束条件下的目标校正参数。所述处理模块还用于根据所述线性内参和所述目标函数,获取所述标定内参数。
在可选的实施方式中,所述处理模块还用于确定每个所述基础矩阵对应的非线性优化函数;所述非线性优化函数用于修正所述线性内参;所述处理模块还用于根据所有的所述非线性优化函数和所述预设约束条件,确定所述目标函数。
第三方面,本发明提供一种电子设备,包括处理器和存储器,所述存储器存储有能够被所述处理器执行的机器可执行指令,所述处理器可执行所述机器可执行指令以实现前述实施方式中任一所述的方法。
第四方面,本发明提供一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现前述实施方式中任一项所述的方法。
相较于现有技术,本发明提供一种相机标定方法、装置、电子设备和计算机可读存储介质,涉及运动相机的图像处理领域。该相机标定方法包括:获取多张运动图像;所述多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像;获取所述多张运动图像中任意两张图像之间的基础矩阵;所述基础矩阵表征所述任意两张图像中基于所述纹理的图像特征点匹配关系;根据所有的所述基础矩阵,获取所述相机的标定内参数;所述标定内参数为所述相机拍摄图像的目标校正参数。使用具有纹理的运动图像,继而获取图像之间的基础矩阵并获取相机的标定内参数,不需依赖场景,如棋盘格或特殊景物的限定,增加了获取相机的标定内参数的通用性;标定内参数用于拍摄图像的目标校正,减少了图像在后处理过程中的误差,使用相机拍摄目标物体并使用标定内参数进行图像处理,可以提高目标物体在2D/3D场景下的结构恢复精度和绝对定位精度。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明实施例提供的一种相机标定方法的流程示意图;
图2为本发明实施例提供的另一种相机标定方法的流程示意图;
图3为本发明实施例提供的另一种相机标定方法的流程示意图;
图4为本发明实施例提供的另一种相机标定方法的流程示意图;
图5为本发明实施例提供的另一种相机标定方法的流程示意图;
图6为本发明实施例提供的一种相机标定装置的方框示意图;
图7为本发明实施例提供的一种电子设备的方框示意图。
图标:相机标定装置-40,获取模块-41,处理模块-42,电子设备-60,存储器-61,处理器-62,通信接口-63。
具体实施方式
下面将结合本发明实施例中附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,术语“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
“张正友标定”方法依赖于固定标定物棋盘格,并且必须离线标定,限制了很多场景下对于相机内参的计算和获取。例如,拍摄图像需要使用变焦相机时,相机内参的离线标定十分麻烦,需要标定相机在不同焦距下的内参;但是,由于相机拍摄图像时,焦距的不可控和不确定性,离线标定数据和实际变焦结果匹配程度较低。应理解,相机自标定方法是一系列应用中,在线获取相机内参的关键技术;当相机焦距因为某种原因变化时,如人为变焦和温度变化引起镜头形变导致焦距变化,相机自标定技术可以在线获得当前相机的内参数,以求解图像坐标到世界坐标的关系,开展一系列基于图像建模2D/3D场景的应用。
为了至少解决上述问题和背景技术提出的不足,本发明实施例提供一种相机标定方法,请参见图1,图1为本发明实施例提供的一种相机标定方法的流程示意图。该相机标定方法包括:
S31,获取多张运动图像。
上述的多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像。应理解,具有纹理的图像必然存在纹理特征,该纹理特征表征运动图像中重复出现的局部特征与该局部特征的排列规则,该纹理特征还可以反映运动图像中的测量物体的表面特征,如该测量物体的表面线条、外部结构特征等。
可以预见的,在一种可能的实施例中,上述的相机运行在预设轨迹可以是相机所处载体(如无人机)飞行在预设航线,或是相机结构被控制做平移、旋转运动等。应理解,运动的相机可以是安装在不同的载体上,如无人机、飞行客机、汽车或机器人等,当载体进行连续运动时,安装在载体上的相机(如相机在无人机中时,则为吊舱)进行与载体匹配的平移、旋转等运动,并在运动过程中获取上述的多张运动图像。
S32,获取多张运动图像中任意两张图像之间的基础矩阵。
该基础矩阵表征任意两张图像中基于纹理的图像特征点匹配关系。例如,获取图像特征点匹配关系时,可以使用、但不限于七点法、八点法以及基于随机均匀一致性采样法等。可以理解的,该图像特征点可以以上述的纹理为参考特征,能更加准确的体现上述任意两张图像之间的图像特征点的匹配关系,使上述的基础矩阵的置信度更高。
S33,根据所有的基础矩阵,获取相机的标定内参数。
该标定内参数为相机拍摄图像的目标校正参数。相机自标定是指通过多张图像中,图像特征点提取和匹配,进而获取相机的标定内参数,并依据标定内参数进行相机的目标校正的过程;可以理解的,上述标定内参数可以是在相机拍摄图像的过程中,就对图像进行校正,即相机使用该标定内参数进行在线自标定。
例如,当相机的焦距因为某种原因发生变化时,如人为变焦、温度变化引起相机镜头的形变等导致焦距发生变化时,可以使用本发明提供的相机标定方法,在线获得当前的相机的标定内参数,根据该标定内参数求解相机拍摄图像(如第一图像)的图像坐标系到世界坐标系的映射关系,并基于上述映射关系实现相机拍摄图像(第一图像)的2D/3D场景应用。
本发明实施例使用的多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像,基于该纹理,相机的标定过程可以不需要依赖场景中的特殊元素,如场景中含有线段、平行线、已知形状的平面物体或立体物体等,使得本发明实施例提供的相机标定方法更具有通用性。
本发明实施例提供的相机标定方法需要相机运动在预设轨迹时采集的运动图像中包含纹理,以便获取图像中的图像特征点,确定任意两张图像之间的基础矩阵,进而获取相机的标定内参数,相机可以使用该标定内参数进行目标校正,使用相机拍摄目标物体并使用标定内参数进行图像处理,可以提高目标物体在2D/3D场景下的结构恢复精度和绝对定位精度。此外,在相机进行在线标定(使用标定内参数进行标定)的过程中,还可以避免相机的离线标定造成的相机焦距与标定内参数不匹配的情况,提高相机的标定精度。
在可选的实施方式中,上述的多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像,为了使上述的基础矩阵更易于处理,在图1的基础上,以上述的多张运动图像中的每张运动图像具有序列标识为例,请参见图2,图2为本发明实施例提供的另一种相机标定方法的流程示意图。上述的S33可以包括:
S331,按照序列标识的顺序,依次解析每张运动图像对应的基础矩阵,以获取相机的线性内参。
该线性内参为相机拍摄图像的待确认校正参数。可以理解的,该线性内参可以是通过至少两个基础矩阵获取得到的线性解,该线性解可以作为相机的初始解(待确认校正参数)。例如,将基础矩阵做奇异值分解,并使用至少两个基础矩阵得到三个线性方程,以便得到线性内参。
S332,根据所有的基础矩阵,建立目标函数。
该目标函数用于确定多张运动图像处于预设约束条件下的目标校正参数。该预设约束条件可以包括,但不限于相机的焦距、扭曲参数等;若目标函数是一种误差函数,该目标校正参数可以表征线性内参的误差。
S333,根据线性内参和目标函数,获取标定内参数。
例如,当目标函数是一种表示误差的函数,根据上述得到的线性内参和目标函数得到的目标校正参数,即可获取相机的标定内参数。
应理解,多张运动图像是相机运动在预设轨迹时拍摄的,该多张运动图像可以是相机在运动过程中分别采集的多张图像,也可以是相机采集的视频数据。当对多张运动图像进行序列标识并对多张运动图像按照序列标识进行排序时,每个基础矩阵都可以包含得到基础矩阵所采用的两张图像的序列标识信息;在一种可能的实施例中,可以使用序列标识对多张运动图像进行筛选,只有两张图像的序列标识小于或等于预设序列间隔,才获取该两张图像之间的基础矩阵,以便获取更优的目标函数,减少标定内参数的误差。
由于图像之间的差异性,基础矩阵之间的关系是不确定的,但是运动图像是具有序列标识的,运动图像结合序列标识还可以反映相机的运动过程,若单独使用线性优化或少数矩阵的回归分析得到目标函数,其会存在较大的误差;因此,为了减少目标函数的误差,在图2的基础上,给出一种可能的实现方式,请参见图3,图3为本发明实施例提供的另一种相机标定方法的流程示意图。上述的S332可以包括:
S332a,确定每个基础矩阵对应的非线性优化函数。
该非线性优化函数用于修正线性内参。该非线性优化函数可以是基础矩阵对应的误差项。
S332b,根据所有的非线性优化函数和预设约束条件,确定目标函数。
需要说明的是,虽然上述以目标函数是非线性优化得到的,目标函数可以是多种不同的形式,在一些可能的实施例中,目标函数还可以是通过多个线性优化得到的,或是拟合多个基础矩阵对应的曲线得到,或是使用高阶多项式计算得到的。
为了便于理解上述的S31~S33及其可能的子步骤对应的相机标定方法,针对获取线性内参、目标函数和标定内参数给出一种可能的实施例:
令基础矩阵为F,根据标定内参数得到相机内参矩阵为K:
Figure 987591DEST_PATH_IMAGE002
其中,相机内参矩阵
Figure 246009DEST_PATH_IMAGE004
中包括五个未知参数:
Figure 353643DEST_PATH_IMAGE006
对基础矩阵F进行奇异值分解:
Figure 217694DEST_PATH_IMAGE008
其中,矩阵
Figure 376274DEST_PATH_IMAGE010
的大小为3×3,矩阵
Figure 401998DEST_PATH_IMAGE010
的列向量为基础矩阵F的左奇异向量;矩阵
Figure 262507DEST_PATH_IMAGE012
的大小为3×3,矩阵
Figure 930249DEST_PATH_IMAGE012
的列向量为基础矩阵F的右奇异向量;矩阵
Figure 943335DEST_PATH_IMAGE014
的大小为3×3,矩阵
Figure 405541DEST_PATH_IMAGE014
为基础矩阵F奇异值分解后得到的对角矩阵,位于矩阵
Figure 753345DEST_PATH_IMAGE014
的对角线上的元素为基础矩阵F的奇异值。
建立Kruppa方程对应的线性等式:
Figure 959199DEST_PATH_IMAGE016
其中,
Figure 826792DEST_PATH_IMAGE018
为矩阵
Figure 725478DEST_PATH_IMAGE020
的第
Figure 560578DEST_PATH_IMAGE022
行,
Figure 304543DEST_PATH_IMAGE024
为矩阵
Figure 295152DEST_PATH_IMAGE026
的第
Figure 99160DEST_PATH_IMAGE022
行,
Figure 421557DEST_PATH_IMAGE028
为对角矩阵
Figure 969213DEST_PATH_IMAGE030
中位置为(1,1)的对角元素,
Figure 811398DEST_PATH_IMAGE032
为对角矩阵
Figure 51886DEST_PATH_IMAGE014
中位置为(2,2)的对角元素。通过上述的线性等式可以获取相机内参矩阵K的五个未知参数
Figure 596000DEST_PATH_IMAGE034
的初始解(线性内参)。
上述的线性等式引入了三个误差项,分别为:
Figure 947347DEST_PATH_IMAGE036
Figure 909618DEST_PATH_IMAGE038
Figure 321008DEST_PATH_IMAGE040
上述的
Figure 352418DEST_PATH_IMAGE042
,为线性等式中两两相减的绝对值。
建立误差函数,即目标函数
Figure 241876DEST_PATH_IMAGE044
包括多种情况:
第一种,当K中的五个参数
Figure 58654DEST_PATH_IMAGE046
均为未知时,目标函数
Figure 640945DEST_PATH_IMAGE048
为:
Figure 159651DEST_PATH_IMAGE050
其中,N表示基础矩阵F的个数,
Figure 852800DEST_PATH_IMAGE048
为对每个基础矩阵F的三个误差项
Figure 989996DEST_PATH_IMAGE052
Figure 743189DEST_PATH_IMAGE054
的值相加,通过N个基础矩阵F的迭代,将目标函数
Figure 749191DEST_PATH_IMAGE048
的值降低至收敛,进而获取相机的内参矩阵K。目标函数收敛的过程可以是满足预设收敛值,也可以是目标函数的值不再缩小等。
第二种,若相机的像素是正方形,则确定相机的扭曲参数为零,即
Figure 246031DEST_PATH_IMAGE056
=0;当相机的镜头在正反两个方向的焦距是相等的,即
Figure 771821DEST_PATH_IMAGE058
是已知的,则可以在目标函数中引入约束项,目标函数
Figure 289390DEST_PATH_IMAGE060
为:
Figure 658055DEST_PATH_IMAGE062
其中,
Figure 568373DEST_PATH_IMAGE064
Figure 338883DEST_PATH_IMAGE066
为相机在正反两个方向上迭代后获取的焦距值,
Figure DEST_PATH_IMAGE068
为迭代后获取的相机的扭曲参数值,
Figure DEST_PATH_IMAGE070
Figure 105982DEST_PATH_IMAGE066
在第一次迭代时使用上述线性等式得到的初始解(线性内参);
Figure DEST_PATH_IMAGE072
Figure DEST_PATH_IMAGE074
是相机在不同使用情况下的可调整参数,
Figure 24259DEST_PATH_IMAGE072
Figure 128482DEST_PATH_IMAGE074
可以根据相机处于不同载体,如无人机、客车,拍摄物体为平面图形或立体几何图像等不同的使用情况进行确定和调整;通过上述的目标函数
Figure 88917DEST_PATH_IMAGE060
对线性内参进行优化,即可获取相机内参矩阵K的未知参数
Figure DEST_PATH_IMAGE076
第三种,若相机的主点信息是已知的,即
Figure 213868DEST_PATH_IMAGE076
已知,则目标函数
Figure DEST_PATH_IMAGE078
为:
Figure 635753DEST_PATH_IMAGE080
其中,
Figure 137141DEST_PATH_IMAGE082
是在迭代后获取的相机的主点信息,
Figure 616664DEST_PATH_IMAGE082
在第一次迭代时使用上述线性等式得到的初始解(线性内参)。
Figure DEST_PATH_IMAGE084
是相机在不同使用情况下的可调整参数,
Figure 132090DEST_PATH_IMAGE084
可以根据相机处于不同载体,如无人机、客车,拍摄物体为平面图形或立体几何图像等不同的使用情况进行确定和调整;通过上述的目标函数
Figure 821698DEST_PATH_IMAGE078
对线性内参进行优化,即可获取相机内参矩阵K的未知参数
Figure DEST_PATH_IMAGE086
可以理解的,使用本发明提供的相机标定方法,确定多张运动图像,进而获取基础矩阵,获取到线性内参后利用目标函数进行非线性优化,获得相机的标定内参数。该相机标定方法可以用于相机载体运动的应用中,如无人机携带相机、无人车载相机和机器人视觉传感器,用于在线获取相机内参数,作为一系列2D/3D场景信息获取的参数输入,建立图像坐标和世界坐标正确的映射关系;还可以用于实时建图与定位、运动恢复结构等场景。
在可选的实施方式中,若多张运动图像之间的重叠率或采集运动图像之间的时间间隔较大,得到的基础矩阵误差较大,为了解决上述问题,在图2的基础上,以序列标识包括序列标识对应的运动图像的图像采集时间信息为例,请参见图4,图4为本发明实施例提供的另一种相机标定方法的流程示意图。上述的S32可以包括:
S321,确定多张运动图像中符合第一预设条件的至少三张关键图像。
该第一预设条件包括任意两张运动图像的重叠率大于或等于第一阈值,和/或任意两张运动图像的图像采集时间信息符合预设采集时间间隔。
可以理解的,第一阈值可以是任意两张运动图像的重叠率为40%,在一种可能的实施例中,还可以是对重叠率设置第一阈值区间,当任意两张运动图像的重叠率处于第一阈值区间时,确定运动图像符合第一预设条件;例如,设置第一阈值区间为重叠率大于或等于30%、且小于或等于50%。应理解,由于多张运动图像是相机运动在预设轨迹时拍摄的,该重叠率指的是任意两张运动图像之间的航向重叠率。
应理解,当序列标识中包括运动图像的图像采集时间信息时,可以根据相机处于不同的载体,如无人机、客车等,设置不同的预设采集时间间隔,如5秒、5分钟等,进而根据运动图像的图像采集时间信息是否符合预设时间间隔,确定多张运动图像中的关键图像。
需要说明的是,在一种可能的情况下,第一预设条件可以是任意两张运动图像的重叠率大于或等于第一阈值,或任意两张运动图像的图像采集时间信息符合预设采集时间间隔中的其中一个,即当运动图像满足上述的任意一个,就可以认为运动图像是符合第一预设条件;在另一种可能的情况下,第一预设条件还可以是任意两张运动图像的重叠率大于或等于第一阈值,和任意两张运动图像的图像采集时间信息符合预设采集时间间隔,即当运动图像同时满足重叠率和图像采集时间的预设时间间隔,才认为运动图像是符合第一预设条件;在另一种可能的情况下,第一预设条件还可以包括更多的筛选条件,如运动图像的清晰度、曝光度、符合图像处理条件的图像特征点的个数等。
S322,获取至少三张关键图像中的任意两张关键图像的基础矩阵。
可以理解的,通过第一预设条件筛选出关键图像后,得到任意两张关键图像的基础矩阵。需要说明的是,为了获取线性内参,需要至少两个基础矩阵,对两个基础矩阵做奇异值分解后,使用Kruppa方程得到三个线性方程,并根据上述的三个线性方程得到初始的线性解,即线性内参。例如,上述获取基础矩阵的方法可以包括,但不限于七点法、八点法以及基于随机均匀一致性采样法等。
在可选的实施方式中,由于获取基础矩阵的方法众多,但各个方法的误差并不相同,为了解决上述问题,在图4的基础上,给出一种可能的实现方式,请参见图5,图5为本发明实施例提供的另一种相机标定方法的流程示意图。上述的S322可以包括:
S322a,获取任意两张关键图像中的纹理符合第二预设条件的图像特征点。
可以理解的,基于关键图像中的纹理,获取符合第二预设条件的图像特征点,该第二预设条件可以是关键图像中的图像特征点与在纹理上,也可以是图像特征点与纹理的像素距离小于或等于某个阈值。由于该图像特征点是与关键图像的纹理相关的,使得上述获取的图像特征点能更为准确的确定目标函数,减少相机的标定内参数的误差。需要说明的是,提取关键图像中的图像特征点可以使用,但不限于尺度不变特征变换(Scale-Invariant Feature Transform,SIFT)算法、特征增强(Speeded Up Robust Features,SURF)算法、快速特征点提取和描述(Oriented FAST and Rotated BRIEF,ORB)算法以及其它图像特征点提取的算法等。
S322b,基于图像特征点的匹配关系确定基础矩阵。
例如,对提取的图像特征点进行匹配的方法可以包括,但不限于暴力匹配、基于词典的匹配、基于运动均匀一致性的匹配方法等,其目的在于获取任意两张关键图像之间的图像特征点之间的匹配关系,即基础矩阵。
应理解,使用具有纹理的运动图像确定关键图像,对关键图像进行图像特征点提取时,该图像特征点是与关键图像的纹理相关的,从而解决现有技术中图像中必须包含特殊场景的弊端,提高相机标定的通用性。
为了实现上述的S31~S33对应的相机标定方法,本发明实施例提供一种相机标定装置,请参见图6,图6为本发明实施例提供的一种相机标定装置的方框示意图。该相机标定装置40包括:获取模块41和处理模块42。
获取模块41用于获取多张运动图像。多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像。
处理模块42用于获取多张运动图像中任意两张图像之间的基础矩阵。基础矩阵表征任意两张图像中基于纹理的图像特征点匹配关系。
处理模块42还用于根据所有的基础矩阵,获取相机的标定内参数。标定内参数为相机拍摄图像的目标校正参数。
应理解,获取模块41和处理模块42可以协同实现上述S31~S33及其可能的子步骤。
在可选的实施方式中,多张运动图像中的每张运动图像具有序列标识,处理模块42还用于按照序列标识的顺序,依次解析每张运动图像对应的基础矩阵,以获取相机的线性内参。线性内参为相机拍摄图像的待确认校正参数。处理模块42还用于根据所有的基础矩阵,建立目标函数。目标函数用于确定多张运动图像处于预设约束条件下的目标校正参数。处理模块42还用于将线性内参输入至目标函数,获取标定内参数。
应理解,处理模块42可以实现上述S331~S333及其可能的子步骤。
在可选的实施方式中,处理模块42还用于确定每个基础矩阵对应的非线性优化函数。非线性优化函数用于修正线性内参。处理模块42还用于根据所有的非线性优化函数和预设约束条件,确定目标函数。
应理解,处理模块42可以实现上述S332a~S332b及其可能的子步骤。
在可选的实施方式中,每张运动图像的序列标识包括每张运动图像的图像采集时间信息。处理模块42还用于确定多张运动图像中符合第一预设条件的至少三张关键图像。第一预设条件包括任意两张运动图像的重叠率大于或等于第一阈值,和/或任意两张运动图像的图像采集时间信息符合预设采集时间间隔。处理模块42还用于获取至少三张关键图像中的任意两张关键图像的基础矩阵。
应理解,处理模块42可以实现上述S321~S322及其可能的子步骤。
在可选的实施方式中,处理模块42还用于获取任意两张关键图像中的纹理符合第二预设条件的图像特征点。处理模块42还用于基于图像特征点的匹配关系确定基础矩阵。
应理解,处理模块42可以实现上述S322a~S322b及其可能的子步骤。
为了实现上述的相机标定方法,本发明实施例提供一种电子设备,如图7,图7为本发明实施例提供的一种电子设备的方框示意图。该电子设备60包括存储器61、处理器62和通信接口63。该存储器61、处理器62和通信接口63相互之间直接或间接地电性连接,以实现数据的传输或交互。例如,这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。存储器61可用于存储软件程序及模块,如本发明实施例所提供的相机标定方法对应的程序指令/模块,处理器62通过执行存储在存储器61内的软件程序及模块,从而执行各种功能应用以及数据处理。该通信接口63可用于与其他节点设备进行信令或数据的通信。在本发明中该电子设备60可以具有多个通信接口63。
其中,存储器61可以是但不限于,随机存取存储器(Random Access Memory,RAM),只读存储器(Read Only Memory,ROM),可编程只读存储器(Programmable Read-OnlyMemory,PROM),可擦除只读存储器(Erasable Programmable Read-Only Memory,EPROM),电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory,EEPROM)等。
处理器62可以是一种集成电路芯片,具有信号处理能力。该处理器可以是通用处理器,包括中央处理器(Central Processing Unit,CPU)、网络处理器(NetworkProcessor,NP)等;还可以是数字信号处理器(Digital Signal Processing,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。
电子设备60可以实现本发明提供的任一种相机标定方法。需要说明的是,该电子设备60可以是集成在相机上的处理设备,实现对相机的在线标定;该电子设备60还可以是设置在载体上,并与相机连接的具有处理能力的处理设备;该电子设备60还可以是与相机保持即时通信连接的手机、平板电脑、笔记本电脑、服务器或其它具有处理能力的电子设备,以便实现对相机的在线标定。
在本发明所提供的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的***来实现,或者可以用专用硬件与计算机指令的组合来实现。
另外,在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
综上所述,本发明提供一种相机标定方法、装置、电子设备和计算机可读存储介质,涉及运动相机的图像处理领域。该相机标定方法包括:获取多张运动图像;多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像;获取多张运动图像中任意两张图像之间的基础矩阵;基础矩阵表征任意两张图像中基于纹理的图像特征点匹配关系;根据所有的基础矩阵,获取相机的标定内参数;标定内参数为相机拍摄图像的目标校正参数。使用具有纹理的运动图像,继而获取图像之间的基础矩阵并获取相机的标定内参数,不需依赖场景,如棋盘格或特殊景物的限定,增加了获取相机的标定内参数的通用性;标定内参数用于拍摄图像的目标校正,减少了图像在后处理过程中的误差,使用相机拍摄目标物体并使用标定内参数进行图像处理,可以提高目标物体在2D/3D场景下的结构恢复精度和绝对定位精度。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种相机标定方法,其特征在于,所述方法包括:
获取多张运动图像;所述多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像;
获取所述多张运动图像中任意两张图像之间的基础矩阵;所述基础矩阵表征所述任意两张图像中基于所述纹理的图像特征点匹配关系,所述多张运动图像中的每张运动图像具有序列标识;
按照所述序列标识的顺序,依次解析所述每张运动图像对应的所述基础矩阵,以获取所述相机的线性内参;所述线性内参为所述相机拍摄图像的待确认校正参数;
根据所有的所述基础矩阵,建立目标函数;所述目标函数用于确定所述多张运动图像处于预设约束条件下的目标校正参数;
根据所述线性内参和所述目标函数,获取所述相机的标定内参数;所述标定内参数为所述相机拍摄图像的目标校正参数。
2.根据权利要求1所述的方法,其特征在于,所述根据所有的所述基础矩阵,建立目标函数,包括:
确定每个所述基础矩阵对应的非线性优化函数;所述非线性优化函数用于修正所述线性内参;
根据所有的所述非线性优化函数和所述预设约束条件,确定所述目标函数。
3.根据权利要求1或2所述的方法,其特征在于,所述每张运动图像的序列标识包括所述每张运动图像的图像采集时间信息,所述获取所述多张运动图像中任意两张图像之间的基础矩阵,包括:
确定所述多张运动图像中符合第一预设条件的至少三张关键图像;所述第一预设条件包括任意两张运动图像的重叠率大于或等于第一阈值,和/或所述任意两张运动图像的图像采集时间信息符合预设采集时间间隔;
获取所述至少三张关键图像中的任意两张关键图像的基础矩阵。
4.根据权利要求3所述的方法,其特征在于,获取所述至少三张关键图像中的任意两张关键图像的基础矩阵,包括:
获取所述任意两张关键图像中的所述纹理符合第二预设条件的图像特征点;
基于所述图像特征点的匹配关系确定所述基础矩阵。
5.一种相机标定装置,其特征在于,所述装置包括:获取模块和处理模块;
所述获取模块,用于获取多张运动图像;所述多张运动图像是相机运动在预设轨迹时采集的具有纹理的图像;
所述处理模块,用于获取所述多张运动图像中任意两张图像之间的基础矩阵;所述基础矩阵表征所述任意两张图像中基于所述纹理的图像特征点匹配关系,所述多张运动图像中的每张运动图像具有序列标识;
所述处理模块还用于按照所述序列标识的顺序,依次解析所述每张运动图像对应的所述基础矩阵,以获取所述相机的线性内参;所述线性内参为所述相机拍摄图像的待确认校正参数;
所述处理模块还用于根据所有的所述基础矩阵,建立目标函数;所述目标函数用于确定所述多张运动图像处于预设约束条件下的目标校正参数;
所述处理模块还用于根据所述线性内参和所述目标函数,获取所述相机的标定内参数;所述标定内参数为所述相机拍摄图像的目标校正参数。
6.根据权利要求5所述的装置,其特征在于,所述处理模块还用于确定每个所述基础矩阵对应的非线性优化函数;所述非线性优化函数用于修正所述线性内参;
所述处理模块还用于根据所有的所述非线性优化函数和所述预设约束条件,确定所述目标函数。
7.一种电子设备,其特征在于,包括处理器和存储器,所述存储器存储有能够被所述处理器执行的机器可执行指令,所述处理器可执行所述机器可执行指令以实现权利要求1-4中任一项所述的方法。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1-4中任一项所述的方法。
CN202010352903.1A 2020-04-29 2020-04-29 相机标定方法、装置、电子设备和计算机可读存储介质 Active CN111243035B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010352903.1A CN111243035B (zh) 2020-04-29 2020-04-29 相机标定方法、装置、电子设备和计算机可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010352903.1A CN111243035B (zh) 2020-04-29 2020-04-29 相机标定方法、装置、电子设备和计算机可读存储介质

Publications (2)

Publication Number Publication Date
CN111243035A CN111243035A (zh) 2020-06-05
CN111243035B true CN111243035B (zh) 2020-08-14

Family

ID=70873762

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010352903.1A Active CN111243035B (zh) 2020-04-29 2020-04-29 相机标定方法、装置、电子设备和计算机可读存储介质

Country Status (1)

Country Link
CN (1) CN111243035B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113838106B (zh) * 2020-06-08 2024-04-05 浙江宇视科技有限公司 一种图像配准方法、装置及电子设备和存储介质
CN113490966B (zh) * 2020-10-13 2024-01-12 深圳市大疆创新科技有限公司 摄像机参数的标定方法、图像处理方法、装置及存储介质
CN112698053B (zh) * 2020-12-02 2022-02-22 武汉理工大学 粒子图像测速误差优化方法、***和存储介质
CN112598751A (zh) * 2020-12-23 2021-04-02 Oppo(重庆)智能科技有限公司 标定方法及装置、终端和存储介质
CN114677445A (zh) * 2020-12-25 2022-06-28 北京小米移动软件有限公司 相机标定方法、相机标定装置及存储介质
CN112766395B (zh) * 2021-01-27 2023-11-28 中国地质大学(北京) 一种图像匹配方法、装置、电子设备及可读存储介质
CN113140009A (zh) * 2021-03-18 2021-07-20 北京三快在线科技有限公司 一种标定方法、装置、电子设备及存储介质
CN113379853B (zh) * 2021-08-13 2021-11-23 腾讯科技(深圳)有限公司 相机内参的获取方法、装置、设备及可读存储介质
CN115170674B (zh) * 2022-07-20 2023-04-14 禾多科技(北京)有限公司 基于单张图像的相机主点标定方法、装置、设备和介质

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101373538A (zh) * 2008-08-25 2009-02-25 同济大学 平移初值操作的自标定方法
CN101794448B (zh) * 2010-04-07 2012-07-04 上海交通大学 主从摄像机***的全自动标定方法
CN101876532B (zh) * 2010-05-25 2012-05-23 大连理工大学 测量***中的摄像机现场标定方法
US8797387B2 (en) * 2011-04-27 2014-08-05 Aptina Imaging Corporation Self calibrating stereo camera
CN102750704B (zh) * 2012-06-29 2015-04-22 吉林大学 一种摄像机分步自标定方法
CN104751486B (zh) * 2015-03-20 2017-07-11 安徽大学 一种多ptz相机的运动目标接力跟踪算法
CN107067441B (zh) * 2017-04-01 2020-02-11 海信集团有限公司 摄像机标定方法及装置
CN107133987B (zh) * 2017-05-16 2019-07-19 西北工业大学 无重叠视场的相机阵列优化标定方法
CN110047107B (zh) * 2018-01-13 2021-04-13 深圳华大智造科技股份有限公司 相机标定方法及图像配准方法、基因测序仪及***
CN108288294A (zh) * 2018-01-17 2018-07-17 视缘(上海)智能科技有限公司 一种3d相机群的外参标定方法

Also Published As

Publication number Publication date
CN111243035A (zh) 2020-06-05

Similar Documents

Publication Publication Date Title
CN111243035B (zh) 相机标定方法、装置、电子设备和计算机可读存储介质
CN110176032B (zh) 一种三维重建方法及装置
CN110969667A (zh) 基于边缘特征的多光谱相机外参自校正算法
JPH09212648A (ja) 動画像処理方法
US10163257B2 (en) Constructing a 3D structure
CN107507277A (zh) 三维点云重构方法和装置、服务器及可读存储介质
US20210158081A1 (en) System and method for correspondence map determination
CN116433737A (zh) 一种激光雷达点云与图像配准的方法、装置及智能终端
CN112257713A (zh) 图像处理方法、装置、电子设备和计算机可读存储介质
CN115457176A (zh) 一种图像生成方法、装置、电子设备及存储介质
CN116051736A (zh) 一种三维重建方法、装置、边缘设备和存储介质
CN114812558A (zh) 一种结合激光测距的单目视觉无人机自主定位方法
CN111161348A (zh) 一种基于单目相机的物***姿估计方法、装置及设备
Zhang et al. Data association between event streams and intensity frames under diverse baselines
CN117315138A (zh) 基于多目视觉的三维重建方法及***
KR101673144B1 (ko) 부분 선형화 기반의 3차원 영상 정합 방법
US9378555B2 (en) Enhanced outlier removal for 8 point algorithm used in camera motion estimation
Adis et al. D3DLO: Deep 3D LiDAR Odometry
Georgiev et al. A fast and accurate re-calibration technique for misaligned stereo cameras
CN114757984A (zh) 光场相机的场景深度估计方法及装置
JPWO2018235256A1 (ja) ステレオ計測装置及びシステム
CN113592917A (zh) 一种摄像机目标交接方法和交接***
WO2024027364A1 (zh) 基于多域时空数据的场景结构动态估计方法、设备及存储介质
CN117649454B (zh) 双目相机外参自动校正方法、装置、电子设备及存储介质
CN113781583B (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
CP03 Change of name, title or address

Address after: 7 / F, area a, building 6, No. 200, Tianfu 5th Street, high tech Zone, Chengdu, Sichuan 610000

Patentee after: CHENGDU JOUAV AUTOMATION TECHNOLOGY Co.,Ltd.

Country or region after: China

Address before: Room 801-805, 8th floor, Building A, No. 200, Tianfu Wujie, Chengdu High-tech Zone, Sichuan Province, 610000

Patentee before: CHENGDU JOUAV AUTOMATION TECHNOLOGY Co.,Ltd.

Country or region before: China

CP03 Change of name, title or address