CN102314686A - 一种拼接类全景视频的基准视场确定方法、系统及设备 - Google Patents
一种拼接类全景视频的基准视场确定方法、系统及设备 Download PDFInfo
- Publication number
- CN102314686A CN102314686A CN201110221223A CN201110221223A CN102314686A CN 102314686 A CN102314686 A CN 102314686A CN 201110221223 A CN201110221223 A CN 201110221223A CN 201110221223 A CN201110221223 A CN 201110221223A CN 102314686 A CN102314686 A CN 102314686A
- Authority
- CN
- China
- Prior art keywords
- visual field
- video camera
- root node
- projective transformation
- video
- 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
Images
Abstract
本发明适用于图像处理技术领域,提供了一种拼接类全景视频的基准视场确定方法、系统及全景视频处理设备。本发明从找出整体形变度量最小的投影变换出发,将与最小形变度量对应的目标变换视场确定为基准视场,提高了拼接类全景视频的整体质量,较好地解决了摄像机阵列全景视频拼接中由于基准视场选取不当而造成的全景视频整体质量下降的问题。
Description
技术领域
本发明属于图像处理技术领域,尤其涉及一种拼接类全景视频的基准视场确定方法、系统及设备。
背景技术
目前,实时全景视频的生成一般采取两类方式:第一类是利用广角镜头或全景反射光学镜头生成全景视频;第二类是利用多个摄像机获取的多个视角的视频经过变换、拼接、融合为一个包含各视角视频内容的大视场全景视频。
其中,在第二类方式中,将处于一个近似成像平面上的多摄像机组成摄像机阵列,将阵列中的摄像机置于方向基本一致的前视状态,阵列中相邻摄像机之间的视场具有一定的重叠,则阵列中的全部摄像机的视场可组成一个高分辨的大视场。要将这样一个摄像机阵列中各个摄像机的视频实时拼接起来组成一个实时全景视频,一般来说又有两种方法。其中第一种方法是采用全人工的方式,这类方式对阵列中各摄像机之间的参数一致性要求非常高,人工调校的难度非常大;第二种方法是采用计算机自动学习全景拼接参数,并完成视频全景的自动拼接。在第二种方法中一般只需要操作人员将阵列中的每个摄像机的视场置于一个基本一致的平面上,并将相邻摄像机之间的视场粗略调校到有一定的重叠,而其它的任务则由计算机自动完成。
上述第二种方法在进行阵列全景视频自动拼接过程中存在一个关键问题:基准视场的选取问题。在上述自由摄像机阵列全景视频的拼接中,阵列中各摄像机的视场一般不会在一个严格的平面上,视场之间的重叠区域也是事先未定标的,因此在摄像机阵列全景视频生成时,需要选定某一个摄像机视场所在的坐标系作为基准的全景参考坐标系,而将阵列中的其它摄像机的视场通过空间投影变换到该参考坐标系下,然后利用融合技术生成全景视频。上述作为基准的参考坐标系是阵列中某个摄像机的视场坐标系,因此也称该坐标系为基准视场。
在摄像机阵列中,基准视场的选择对全景视频的生成质量影响较大。一般来说,选择不同的基准视场,在生成全景视频时由各个摄像机视场对基准视场投影变换后产生的视场误差也不相同。若整体的误差越小则生成的全景视频的质量越好。因此,如何选择误差最小的基准视场是摄像机阵列全景视频拼接中的一个关键问题之一。
在摄像机阵列全景视频生成中,现有的技术通常是人为地在摄像机阵列中选取其中一个摄像机的视场作为基准视场,其它视场按照映射关系投影变换到该基准视场的坐标系中。这种基准视场的人为选择方法在摄像机阵列规模较小时效果尚可,但在摄像机阵列较大时,由于阵列中的每个摄像机的视场都具有一定的自由性,因此,人为选出的基准视场往往并非最佳的基准视场。
总之,现有技术中针对摄像机阵列全景视频拼接中最优基准视场的自动选取尚没有一个较好的解决方法。
发明内容
本发明所要解决的技术问题在于提供一种拼接类全景视频的基准视场确定方法,旨在通过摄像机阵列生成全景视频时选择最优的基准视场,提高拼接类全景视频的整体质量。
本发明是这样实现的,一种拼接类全景视频的基准视场确定方法,包括下述步骤:
步骤A,在预拼接视频的摄像机阵列中,对每一个摄像机找出与其具有公共重叠视场的邻接摄像机子集合,并计算每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵;
步骤B,于所述摄像机阵列中选取摄像机i,根据步骤A得到的每个摄像机的邻接摄像机子集合及其投影变换矩阵,将所述摄像机阵列中的其他摄像机的视场均通过最少次数的投影变换处理变换到选取的摄像机i的投影视场,并计算其他摄像机的视场在变换前后的几何形变度量之和δi;
步骤C,按照i为1至n的顺序重复步骤B,得到n个几何形变度量之和,在选取过的n个摄像机中,将与其中最小的几何形变度量之和相对应的摄像机的视场确定为基准视场;其中n为所述的摄像机阵列中摄像机的数量。
进一步地,所述步骤B具体包括下述步骤:
步骤B1,由每个摄像机的邻接子集合,构造描述摄像机阵列的视场连通关系的视场邻接关系图;所述视场邻接关系图中包括若干节点和连接在节点之间的连接线,其中每一个节点表示所述摄像机阵列中的一个摄像机视场,之间具有连接线的两个节点表示所代表的两个摄像机的视场相邻接;
步骤B2,在步骤B1中得到的视场邻接关系图中,任意选取其中的一个节点作为当前根节点,生成该当前根节点的单源最短路径树;在所述单源最短路径树中,任一非根节点均可经过最少的节点与当前根节点直接或间接相连通;
步骤B3,根据每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵,计算所述单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵;
步骤B4,根据所述单源最短路径树和任一非根节点到当前根节点的投影变换矩阵,将每个非根节点所代表的摄像机的视场进行投影变换处理以变换到当前根节点所代表的摄像机的视场;
步骤B5,计算所述单源最短路径树中的每个非根节点所代表的摄像机视场在投影变换前后的视场形变度量,并对每个非根节点所代表的形变度量求和。
进一步地,设定所述单源最短路径树中从当前根节点k到非根节点p的路径经过的节点序列为:k、k1、...、p,对应的视场分别为Vk0、Vk1、...、Vkm,其中,Vk0=Vk,Vkm=Vp,所述步骤B3通过下述公式计算所述单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵:
Hkp=Hk1Hk2…Hkm
其中,Hkp为非根节点p的视场Vp到当前根节点k的视场Vk的投影变换矩阵,Hki为由视场Vki到视场Vk(i-1)的投影变换矩阵,i=1,2,...,m。
本发明还提供了一种拼接类全景视频的基准视场确定系统,包括:
邻接摄像机子集合查找单元,用于在预拼接视频的摄像机阵列中,对每一个摄像机找出与其具有公共重叠视场的邻接摄像机子集合,并计算每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵;
几何形变度量计算单元,用于在所述摄像机阵列中选取摄像机i,根据所述邻接摄像机子集合查找单元得到的每个摄像机的邻接摄像机子集合及其投影变换矩阵,将所述摄像机阵列中的其他摄像机的视场均通过最少次数的投影变换处理变换到选取的摄像机i的投影视场,并计算其他摄像机的视场在变换前后的几何形变度量之和δi,并按照i为1至n的顺序重复步骤B,得到n个几何形变度量之和,其中n为所述的摄像机阵列中摄像机的数量;
基准视场确定单元,用于在选取过的n个摄像机中,将所述几何形变度量计算单元计算得到的n个几何形变度量之和中最小的一个所对应的摄像机的视场确定为基准视场。
进一步地,所述几何形变度量计算单元包括:
视场连通关系的视场邻接关系图构造模块,用于根据所述邻接摄像机子集合查找单元查找到的每个摄像机的邻接子集合,构造描述摄像机阵列的视场连通关系的视场邻接关系图;所述视场邻接关系图中包括若干节点和连接在节点之间的连接线,其中每一个节点表示所述摄像机阵列中的一个摄像机视场,之间具有连接线的两个节点表示所代表的两个摄像机的视场相邻接;
单源最短路径树生成模块,用于在所述视场连通关系的视场邻接关系图构造模块得到的视场邻接关系图中,任意选取其中的一个节点作为当前根节点,生成该当前根节点的单源最短路径树;在所述单源最短路径树中,任一非根节点均可经过最少的节点与当前根节点直接或间接相连通;
投影变换矩阵计算模块,用于根据每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵,计算所述单源最短路径树生成模块生成的单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵;
投影变换模块,用于根据所述单源最短路径树生成模块生成的单源最短路径树,和所述投影变换矩阵计算模块计算得到的任一非根节点到当前根节点的投影变换矩阵,将每个非根节点所代表的摄像机的视场进行投影变换处理以变换到当前根节点所代表的摄像机的视场;
视场形变度量计算模块,用于在所述投影变化模块进行投影变换处理后,计算所述单源最短路径树中的每个非根节点所代表的摄像机视场在投影变换前后的视场形变度量,并对每个非根节点所代表的形变度量求和。
进一步地,设定所述单源最短路径树中从当前根节点k到非根节点p的路径经过的节点序列为:k、k1、...、p,对应的视场分别为Vk0、Vk1、...、Vkm,其中,Vk0=Vk,Vkm=Vp,所述投影变换矩阵计算模块通过下述公式计算所述单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵:
Hkp=Hk1Hk2…Hkm
其中,Hkp为非根节点p的视场Vp到当前根节点k的视场Vk的投影变换矩阵,Hki为由视场Vki到视场Vk(i-1)的投影变换矩阵,i=1,2,...,m。
进一步地,对于当前根节点所代表的摄像机的邻接摄像机子集合中的摄像机,其对应的非根节点与当前根节点直接相连通;而对于不属于当前根节点所代表的摄像机的邻接摄像机子集合的摄像机,其对应的非根节点则通过与当前根节点直接相连通的非根节点间接连通。
进一步地,所述几何形变度量包括面积形变度量和角度形变度量。
本发明还提供了一种全景视频处理设备,其与一摄像机阵列连接,所述全景视频处理设备包括一如上所述的拼接类全景视频的基准视场确定系统。
本发明从找出整体形变度量最小的投影变换出发,将与最小形变度量对应的目标变换视场确定为基准视场,提高了拼接类全景视频的整体质量,较好地解决了摄像机阵列全景视频拼接中由于基准视场选取不当而造成的全景视频整体质量下降的问题,其原理主要有以下两点:
其一,在全景拼接中,计算出的投影变换矩阵一般存在误差,形变较大的视场对应的误差在平均意义下往往较大;在单源最短路径树中,间接连通视场的投影变换矩阵的计算中,会累积各个级联的投影变换矩阵中的误差,形变小就意味着累积误差在平均意义下比较小。变换误差小则图像的非正常扭曲就较小。
其二,由于目前数字图像一般是以规则的矩形栅格形式显示和存储的,投影变换后的视场中的图像必须经过插值步骤以填充栅格点处的像素。插值实际上是对图像像素值进行的估计,在相同的插值算法下,视场的形变小,则插值产生的整体误差往往较小;反之则较大。插值误差小则图像的失真就较小。
附图说明
图1是本发明实施例提供的拼接类全景视频的基准视场确定方法的实现流程图;
图2是本发明实施例提供的一个3×3摄像机阵列示例的视场分布示意图;
图3是图2所示3×3摄像机阵列的视场邻接子集合的列表;
图4是图1所示步骤B的实现流程图;
图5是图2所示3×3摄像机阵列的视场邻接关系图;
图6A和图6B是图5所示的视场邻接关系图生成的以V1为根节点的两个单源最短路径树的示意图;
图7是图5所示的视场邻接关系图生成的以V5为根节点的单源最短路径树的示意图;
图8是本发明实施例提供的视场投影变换前后的形变示意图;
图9是本发明实施例提供的拼接类全景视频的基准视场确定系统的结构原理图;
图10是图9中几何形变度量计算单元的结构原理图;
图11是本发明实施例提供的全景视频处理设备与摄像机阵列的连接关系图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明实施例中,从找出整体形变度量最小的投影变换出发,将与最小形变度量对应的目标变换视场确定为基准视场,并进一步可以此时的投影变换关系矩阵作为摄像机阵列中该摄像机视场到全景视场的投影变换矩阵。
图1示出了本发明实施例提供的方法的实现流程,详述如下。
在步骤A中,在预拼接视频的摄像机阵列中,对每一个摄像机找出与其具有公共重叠视场的邻接摄像机子集合,并计算每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵。
对摄像机阵列中给定的一个摄像机,寻找与其具有公共重叠视场的邻接摄像机的子集合可以通过计算机自动实现,但目前计算机全自动寻找摄像机阵列中相邻子集合的方法在阵列规模增大时运算量会急剧增大,准确性也会下降。而在实际中,通过人工对各个摄像机视场观察和交互,确定并标注阵列中与该摄像机邻近的摄像机子集合的方法,相对来说是一种高效和高准确性的方法,并且一般的操作人员就可以很好地胜任,此处不再一一列举。
在阵列中与给定摄像机邻近的摄像机子集合确定后,需要计算该摄像机视场与邻近子集合中每个摄像机视场的投影变换矩阵。现有技术中,计算具有部分公共重叠视场的两个摄像机视场之间的投影变换矩阵的方法有多种,其中基于尺度不变特征变换(SIFT)的特征匹配点提取和匹配点对的RANSAC提纯方法的投影矩阵计算方法是较佳的一种,它可以利用公共视场重叠区域的图像特征点对,自动计算出投影变换矩阵,此处不再一一列举。
假设在具有n个摄像机的摄像机阵列中,记编号为k的摄像机的视场为Vk,与Vk有部分公共重叠视场的邻近摄像机的子集合记为{Vki,i=k1,k2,...,km},从视场Vk到视场Vki的投影变换关系矩阵记为Hki,i=k1,k2,...,km,k=1,2,...,n。其中,不同的摄像机由于其在阵列中的位置不同,其具有的视场邻近子集合中的摄像机数量km一般也不同。投影变换矩阵H是一个如公式(1)所示的3×3矩阵,参数h0、h1、h3、h4控制尺度变化和旋转变化,h2控制水平方向位移,h5控制垂直方向位移,h6、h7控制水平和垂直方向形变。
图2是一个3×3摄像机阵列示例的视场分布示意图,其中各个视场及其邻接子集合的情况在图3的列表中示出,表中的“1”表示纵横两个视场有重叠,“0”表示无重叠。图3表中的第k行表示的是第k个摄像机视场Vk的视场邻接子集合。例如,第2行表示视场V2的视场邻接子集合是{V1,V3,V5,V6},该邻接子集合中的视场对于V2的投影变换关系矩阵分别记为H21、H23、H25、H26。
在步骤B中,于所述摄像机阵列中选取摄像机i,根据步骤A得到的每个摄像机的邻接摄像机子集合及其投影变换矩阵,将所述摄像机阵列中的其他摄像机的视场均通过最少次数的投影变换处理变换到选取的摄像机i的投影视场,并计算其他摄像机的视场在变换前后的几何形变度量之和δi。
在一个摄像机阵列中,各个摄像机视场之间可能有着错综复杂的重叠,而从一个摄像机的视场投影变换到另一个摄像机的视场可能有不同的变换途径,比如在图2所示的投影视场分布中,从视场V1投影变换到视场V9,可以按照V1---V5---V8---V9的路径变换三次,也可以按照V1---V4---V7---V8---V9的路径变换四次。而上述的“最少次数的投影变换处理”是指按照最短的路径实现最终的投影变换。
本实施例中,步骤B具体如图4所示,应当理解,具体实施时还可以通过其他方式实现“最少次数的投影变换处理”。
参照图4,步骤B1,由每个摄像机的邻接子集合,构造描述摄像机阵列的视场连通关系的视场邻接关系图。
上述视场邻接关系图中包括若干节点和连接在节点之间的连接线,其中每一个节点表示所述摄像机阵列中的一个摄像机视场,之间具有连接线的两个节点表示所代表的两个摄像机的视场相邻接。图5是图2所示的3×3摄像机阵列的视场连通邻接关系图G(Graph),它可由图3所示的视场邻接子集合列表生成。有关邻接关系图G的生成,以及在计算机中的数据表示方式在计算机科学技术类的教科书中有论述,此处不再赘述。
步骤B2,在步骤B1中得到的视场邻接关系图中,任意选取其中的一个节点作为当前根节点,生成该当前根节点的单源最短路径树;在所述单源最短路径树中,任一非根节点均可经过最少的节点与当前根节点直接或间接相连通。
参照图5,邻接关系图G中任取其中的一个节点Vk作为根节点,生成该节点的单源最短路径树T(Tree)。节点Vk的单源最短路径树Tk是指,该树Tk中任给的一个节点Vi到根节点Vk的路径都是邻接关系图G中的一条通路,且该通路在邻接关系图G中经过的中间节点数最少。图6A和图6B是由图5的邻接关系图选取节点V1为根节点的两个单源最短路径树,图7是由图5的邻接关系图选取节点V5为根节点的单源最短路径树。对于当前根节点所代表的摄像机的邻接摄像机子集合中的摄像机,其对应的非根节点与当前根节点直接相连通;而对于不属于当前根节点所代表的摄像机的邻接摄像机子集合的摄像机,其对应的非根节点则通过与当前根节点直接相连通的非根节点间接连通,以图7为例,与当前根节点V5直接相连通的V1、V2、V4、V6、V8属于当前根节点V5所代表的摄像机的邻接摄像机子集合,而V3、V7、V9不属于当前根节点所代表的摄像机的邻接摄像机子集合,需分别通过非根节点V2、V4、V8与当前根节点V5间接连通。
由图6A和图6B可知,节点Vk的单源最短路径树Tk一般是不唯一的,但是在Vk的单源最短路径树Tk中,每个节点的深度是固定的。其中,节点在树T中的深度是指该节点到达根节点的路径中经过的边的数量。而单源最短路径树的深度定义为该树中具有最大深度的节点的深度。在图2所示的3×3摄像机阵列中,节点V1的单源最短路径树的深度为3,节点V5的单源最短路径树的深度为2。
步骤B3,根据每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵,计算所述单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵。
给定节点Vk的单源最短路径树Tk,计算树Tk中的节点Vp到根节点Vk的投影变换关系矩阵的方法如下:
假设树Tk中从节点Vk到Vp的路径经过的节点序列为:Vk0、Vk1、...、Vkm,其中,Vk0=Vk,Vkm=Vp。由树Tk的构造过程可保证节点Vk(i-1)和Vki,i=1,2,...m,表示的视场具有邻接关系,并假设由视场Vki到Vk(i-1)的投影变换矩阵为Hki,则由视场Vp到Vk的投影变换矩阵可由公式(2)的矩阵乘法公式计算得到。
Hkp=Hk1Hk2…Hkm (2)
其中,Hkp为非根节点p的视场Vp到当前根节点k的视场Vk的投影变换矩阵,Hki为由视场Vki到视场Vk(i-1)的投影变换矩阵,i=1,2,...,m。
由于在阵列中任何一个摄像机视场与其邻接子集合中的每一个摄像机视场的投影变换矩阵已经在步骤A中计算得到,而在树Tk中的每一个节点都与根节点存在一条视场连通的通路,因此每一个节点相对于根节点的投影变换矩阵均可由公式(2)所示的方法计算得到。
步骤B4,根据所述单源最短路径树和任一非根节点到当前根节点的投影变换矩阵,将每个非根节点所代表的摄像机的视场进行投影变换处理以变换到当前根节点所代表的摄像机的视场。
根据在步骤B3中得到的节点Vk的单源最短路径树T,以及树中对应的各个节点对根节点的投影变换关系矩阵,对每个摄像机视场进行投影变换。具体为,假设节点Vp到Vk的投影变换关系矩阵如公式(1)所示,视场Vp的点(x,y)变换为视场Vk中的点(x’,y’),则变换前后对应点对的坐标之间的关系如公式(3)所示:
步骤B5,计算所述单源最短路径树中的每个非根节点所代表的摄像机视场在投影变换前后的视场形变度量,并对每个非根节点所代表的形变度量求和。
为了更好地描述节点Vk的单源最短路径树Tk的视场整体累积形变度量的计算方法,下面首先介绍本发明中对单个的视场在投影变换前后的形变度量计算方法,如下所论述:
由于摄像机的原始视场一般为矩形,而投影变换将原视场的点变换为投影视场的点,将原视场中的直线变换为投影视场中的直线或点,因此原摄像机视场的投影视场一般是一个凸四边形(在极端情况下可能退化为三角形或一条直线,在本发明中不考虑极端情况投影退化的情况)。为了更清楚地说明本发明中视场投影变换前后视场的形变度量计算方法,请参阅图8。具体详述如下:
在图8中,假设左图为投影变换前的视场形状,右图为投影变换后的视场形状。投影变换前矩形视场的四个顶点分别为A0、B0、C0、D0,对角线的交点为O0,投影变换后视场凸四边形的四个顶点分别为A1、B1、C1、D1,对角线的交点为O1;假设投影变换前,视场中的四个三角形A0B0O0、B0C0O0、C0D0O0、D0A0O0的面积分别为:S01、S02、S03、S04,这四个三角形经过投影,分别变换为三角形A1B1O1,B1C1O1,C1D1O1,D1A1O1,并且假设它们的面积分别为:S11、S12、S13、S14;进一步假设,如图8所示的原视场中的角:θ1、θ2、θ3、θ4、θ5、θ6、θ7、θ8分别变换为投影后视场中的角β1、β2、β3、β4、β5、β6、β7、β8。可知原视场的面积S=S01+S02+S03+S04,则原视场在投影变换前后面积的形变度量定义为:
视场在投影变换前后角度的形变度量定义为:
在公式(5)的右边第二项中,θ9=θ1,β9=β1。由公式(4)的面积的形变度量和公式(5)的角度的形变度量一起,构造视场在投影变换前后的几何形变度量如公式(6)所示:
δ=(1+δS)(1+δθ) (6)
定义了上述单个视场在投影变换前后的几何形变度量,下面叙述本发明中节点Vk的单源最短路径树Tk的视场整体累积形变度量的计算方法。如下所述:
假设经过步骤A和B1构造出的摄像机阵列的视场邻接关系图G中共有n个(摄像机)节点。由步骤B3知,在图G中节点Vk的单源最短路径树Tk的根节点为Vk,将树Tk的其它的非根节点分别记为Vki,i=1,2,...,(n-1),非根节点Vki相对于根节点Vk的投影变换关系矩阵Hki,i=1,2,...,(n-1)可由步骤B3计算得到。由投影变换矩阵Hki,i=1,2,...,(n-1)对节点Vki,i=1,2,...,(n-1)对应的视场分别进行投影变换,利用公式(4)、公式(5)和公式(6)分别计算出非根节点节点Vki,i=1,2,...,(n-1)对应的视场经过投影变换后的几何形变度量:δki,i=1,2,...,(n-1)。则节点Vk的单源最短路径树Tk的视场整体累积几何形变度量定义为:
最终得到由该单源最短路径树描述的摄像机阵列连通视场的整体累积形变度量。
若Vk的单源最短路径树Tk不止一株(如图6A和图6B所示),假设共有mk株,记为:Tk(j),j=1,2,...,mk,此时则需要计算出每一株Tk(j)的视场整体累积几何形变度量δk(j),j=1,2,...,mk,并找出视场整体累积几何形变度量最小的那株树Tk(q):
则在这种情况下,令Vk的单源最短路径树:Tk=Tk(q),Tk的视场整体累积几何形变度量则为:δk=δk(q)。
在步骤C中,按照i为1至n的顺序重复步骤B,得到n个几何形变度量之和,在选取过的n个摄像机中,将与其中最小的几何形变度量之和相对应的摄像机的视场确定为基准视场;其中n为所述的摄像机阵列中摄像机的数量。
例如在步骤B1得到的邻接关系图G中,选择不同的节点Vk,k=1,2,...,n作为根节点,重复上述步骤B3、B4至B5的过程,找出节点Vk,k=1,2,...,n的单源最短路径树Tk,k=1,2,...,n;并计算出Tk,k=1,2,...,n的视场整体累积几何形变度量:δk,k=1,2,..,,n;找出具有最小整体几何形变度量的那株树Tp:
则树Tp称为本发明所述摄像机阵列中的最优单源最短路径树,树Tp的根节点Vp所对应的摄像机的视场则作为本发明摄像机阵列全景拼接的基准视场。最优单源最短路径树Tp的每个非根节点Vpi,i=1,2,...,(n-1)到根节点Vp的投影变换关系矩阵Hpi,i=1,2,...,(n-1)作为本发明摄像机阵列全景视频拼接中的视场Vpi,i=1,2,...,(n-1)到基准视场Vp的投影变换矩阵。
由上面的叙述可知,步骤C是一个优化过程。在摄像机阵列规模较小时,该优化过程可以采用穷举的方法实现。当摄像机阵列规模比较大时,为了加快优化过程,可以采用模拟退火、遗传算法、粒子群优化算法等方法。这些优化方法在现有的教科书中有详细的论述,此处不再赘述。
本领域普通技术人员可以理解实现上述各实施例提供的方法中的全部或部分步骤可以通过程序指令及相关的硬件来完成,所述的程序可以存储于一计算机可读取存储介质中,该存储介质可以为ROM/RAM、磁盘、光盘等。
图9示出了本发明实施例提供的拼接类全景视频的基准视场确定系统的结构原理,为了便于描述,仅示出了与本实施例相关的部分。此系统可以为内置于全景视频处理设备中的软件单元、硬件单元或软硬件结合的单元。
参照图9,本发明实施例提供的拼接类全景视频的基准视场确定系统包括邻接摄像机子集合查找单元91、几何形变度量计算单元92、基准视场确定单元93,其中邻接摄像机子集合查找单元91在预拼接视频的摄像机阵列中,对每一个摄像机找出与其具有公共重叠视场的邻接摄像机子集合,并计算每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵。几何形变度量计算单元92在所述摄像机阵列中选取摄像机i,根据所述邻接摄像机子集合查找单元得到的每个摄像机的邻接摄像机子集合及其投影变换矩阵,将所述摄像机阵列中的其他摄像机的视场均通过最少次数的投影变换处理变换到选取的摄像机i的投影视场,并计算其他摄像机的视场在变换前后的几何形变度量之和δi,并按照i为1至n的顺序重复步骤B,得到n个几何形变度量之和,其中n为所述的摄像机阵列中摄像机的数量。基准视场确定单元93在选取过的n个摄像机中,将所述几何形变度量计算单元计算得到的n个几何形变度量之和中最小的一个所对应的摄像机的视场确定为基准视场。
进一步地,图10示出了几何形变度量计算单元92一实施例,包括视场连通关系的视场邻接关系图构造模块921、单源最短路径树生成模块922、投影变换矩阵计算模块923、投影变换模块924、视场形变度量计算模块925。
其中,视场连通关系的视场邻接关系图构造模块921根据所述邻接摄像机子集合查找单元查找到的每个摄像机的邻接子集合,构造描述摄像机阵列的视场连通关系的视场邻接关系图;所述视场邻接关系图中包括若干节点和连接在节点之间的连接线,其中每一个节点表示所述摄像机阵列中的一个摄像机视场,之间具有连接线的两个节点表示所代表的两个摄像机的视场相邻接。
单源最短路径树生成模块922在所述视场连通关系的视场邻接关系图构造模块得到的视场邻接关系图中,任意选取其中的一个节点作为当前根节点,生成该当前根节点的单源最短路径树;在所述单源最短路径树中,任一非根节点均可经过最少的节点与当前根节点直接或间接相连通。
对于当前根节点所代表的摄像机的邻接摄像机子集合中的摄像机,其对应的非根节点与当前根节点直接相连通;而对于不属于当前根节点所代表的摄像机的邻接摄像机子集合的摄像机,其对应的非根节点则通过与当前根节点直接相连通的非根节点间接连通。
投影变换矩阵计算模块923根据每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵,计算所述单源最短路径树生成模块生成的单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵。
设定所述单源最短路径树中从当前根节点k到非根节点p的路径经过的节点序列为:k、k1、...、p,对应的视场分别为Vk0、Vk1、...、Vkm,其中,Vk0=Vk,Vkm=Vp,所述投影变换矩阵计算模块通过下述公式计算所述单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵:
Hkp=Hk1Hk2…Hkm
其中,Hkp为非根节点p的视场Vp到当前根节点k的视场Vk的投影变换矩阵,Hki为由视场Vki到视场Vk(i-1)的投影变换矩阵,i=1,2,...,m。
投影变换模块924根据所述单源最短路径树生成模块生成的单源最短路径树,和所述投影变换矩阵计算模块计算得到的任一非根节点到当前根节点的投影变换矩阵,将每个非根节点所代表的摄像机的视场进行投影变换处理以变换到当前根节点所代表的摄像机的视场。
视场形变度量计算模块925在所述投影变化模块进行投影变换处理后,计算所述单源最短路径树中的每个非根节点所代表的摄像机视场在投影变换前后的视场形变度量,并对每个非根节点所代表的形变度量求和。
上述各单元、模块的实现原理如上文所述,此处不再一一赘述。
本发明还提供了一种全景视频处理设备,如图11所示,此全景视频处理设备111与摄像机阵列112连接,其特征在于,全景视频处理设备中还包括包括一如上所述的拼接类全景视频的基准视场确定系统。
本发明实施例通过对摄像机阵列中的任一摄像机,找出与其具有部分公共重叠视场的摄像机邻接子集合,计算该摄像机视场与邻接子集合中每个摄像机视场的投影变换关系矩阵;由阵列中所有摄像机的邻接子集合及它们的投影变换关系矩阵,构造摄像机阵列的视场邻接关系图;任取邻接关系图中的一个节点作为根节点,由阵列的邻接关系图构造阵列的单源最短路径树,并计算树中每一个节点到根节点的投影变换关系矩阵;根据上述投影变换关系矩阵对树中非根节点的每个摄像机视场进行投影变换,计算上述每个投影变换视场的形变度量,并对每个视场的形变度量求和得到该单源最短路径树描述的阵列连通视场的整体累积形变度量;改变根节点,重复上述过程,找出整体累积形变度量最小的那棵单源最短路径树,该树的根节点所在的摄像机视场作为基准视场,该最优单源最短路径树的每个节点到根节点的投影变换关系矩阵作为摄像机阵列中该摄像机视场到全景视场的投影变换矩阵。本发明提高了拼接类全景视频的整体质量,较好地解决了摄像机阵列全景视频拼接中由于基准视场选取不当而造成的全景视频整体质量下降的问题。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种拼接类全景视频的基准视场确定方法,其特征在于,包括下述步骤:
步骤A,在预拼接视频的摄像机阵列中,对每一个摄像机找出与其具有公共重叠视场的邻接摄像机子集合,并计算每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵;
步骤B,于所述摄像机阵列中选取摄像机i,根据步骤A得到的每个摄像机的邻接摄像机子集合及其投影变换矩阵,将所述摄像机阵列中的其他摄像机的视场均通过最少次数的投影变换处理变换到选取的摄像机i的投影视场,并计算其他摄像机的视场在变换前后的几何形变度量之和δi;
步骤C,按照i为1至n的顺序重复步骤B,得到n个几何形变度量之和,在选取过的n个摄像机中,将与其中最小的几何形变度量之和相对应的摄像机的视场确定为基准视场;其中n为所述的摄像机阵列中摄像机的数量。
2.如权利要求1所述的方法,其特征在于,所述步骤B具体包括下述步骤:
步骤B1,由每个摄像机的邻接子集合,构造描述摄像机阵列的视场连通关系的视场邻接关系图;所述视场邻接关系图中包括若干节点和连接在节点之间的连接线,其中每一个节点表示所述摄像机阵列中的一个摄像机视场,之间具有连接线的两个节点表示所代表的两个摄像机的视场相邻接;
步骤B2,在步骤B1中得到的视场邻接关系图中,任意选取其中的一个节点作为当前根节点,生成该当前根节点的单源最短路径树;在所述单源最短路径树中,任一非根节点均可经过最少的节点与当前根节点直接或间接相连通;
步骤B3,根据每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵,计算所述单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵;
步骤B4,根据所述单源最短路径树和任一非根节点到当前根节点的投影变换矩阵,将每个非根节点所代表的摄像机的视场进行投影变换处理以变换到当前根节点所代表的摄像机的视场;
步骤B5,计算所述单源最短路径树中的每个非根节点所代表的摄像机视场在投影变换前后的视场形变度量,并对每个非根节点所代表的形变度量求和。
3.如权利要求2所述的方法,其特征在于,设定所述单源最短路径树中从当前根节点k到非根节点p的路径经过的节点序列为:k、k1、...、p,对应的视场分别为Vk0、Vk1、...、Vkm,其中,Vk0=Vk,Vkm=Vp,所述步骤B3通过下述公式计算所述单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵:
Hkp=Hk1Hk2…Hkm
其中,Hkp为非根节点p的视场Vp到当前根节点k的视场Vk的投影变换矩阵,Hki为由视场Vki到视场Vk(i-1)的投影变换矩阵,i=1,2,...,m。
4.如权利要求2所述的方法,其特征在于,在所述单源最短路径树中,对于当前根节点所代表的摄像机的邻接摄像机子集合中的摄像机,其对应的非根节点与当前根节点直接相连通;而对于不属于当前根节点所代表的摄像机的邻接摄像机子集合的摄像机,其对应的非根节点则通过与当前根节点直接相连通的非根节点间接连通。
5.如权利要求1所述的方法,其特征在于,所述几何形变度量包括面积形变度量和角度形变度量。
6.一种拼接类全景视频的基准视场确定系统,其特征在于,包括:
邻接摄像机子集合查找单元,用于在预拼接视频的摄像机阵列中,对每一个摄像机找出与其具有公共重叠视场的邻接摄像机子集合,并计算每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵;
几何形变度量计算单元,用于在所述摄像机阵列中选取摄像机i,根据所述邻接摄像机子集合查找单元得到的每个摄像机的邻接摄像机子集合及其投影变换矩阵,将所述摄像机阵列中的其他摄像机的视场均通过最少次数的投影变换处理变换到选取的摄像机i的投影视场,并计算其他摄像机的视场在变换前后的几何形变度量之和δi,并按照i为1至n的顺序重复步骤B,得到n个几何形变度量之和,其中n为所述的摄像机阵列中摄像机的数量;
基准视场确定单元,用于在选取过的n个摄像机中,将所述几何形变度量计算单元计算得到的n个几何形变度量之和中最小的一个所对应的摄像机的视场确定为基准视场。
7.如权利要求6所述的系统,其特征在于,所述几何形变度量计算单元包括:
视场连通关系的视场邻接关系图构造模块,用于根据所述邻接摄像机子集合查找单元查找到的每个摄像机的邻接子集合,构造描述摄像机阵列的视场连通关系的视场邻接关系图;所述视场邻接关系图中包括若干节点和连接在节点之间的连接线,其中每一个节点表示所述摄像机阵列中的一个摄像机视场,之间具有连接线的两个节点表示所代表的两个摄像机的视场相邻接;
单源最短路径树生成模块,用于在所述视场连通关系的视场邻接关系图构造模块得到的视场邻接关系图中,任意选取其中的一个节点作为当前根节点,生成该当前根节点的单源最短路径树;在所述单源最短路径树中,任一非根节点均可经过最少的节点与当前根节点直接或间接相连通;
投影变换矩阵计算模块,用于根据每一个摄像机视场与其邻接子集合中的各个摄像机视场的投影变换矩阵,计算所述单源最短路径树生成模块生成的单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵;
投影变换模块,用于根据所述单源最短路径树生成模块生成的单源最短路径树,和所述投影变换矩阵计算模块计算得到的任一非根节点到当前根节点的投影变换矩阵,将每个非根节点所代表的摄像机的视场进行投影变换处理以变换到当前根节点所代表的摄像机的视场;
视场形变度量计算模块,用于在所述投影变化模块进行投影变换处理后,计算所述单源最短路径树中的每个非根节点所代表的摄像机视场在投影变换前后的视场形变度量,并对每个非根节点所代表的形变度量求和。
8.如权利要求7所述的系统,其特征在于,设定所述单源最短路径树中从当前根节点k到非根节点p的路径经过的节点序列为:k、k1、...、p,对应的视场分别为Vk0、Vk1、...、Vkm,其中,Vk0=Vk,Vkm=Vp,所述投影变换矩阵计算模块通过下述公式计算所述单源最短路径树中的任一非根节点到当前根节点的投影变换矩阵:
Hkp=Hk1Hk2…Hkm
其中,Hkp为非根节点p的视场Vp到当前根节点k的视场Vk的投影变换矩阵,Hki为由视场Vki到视场Vk(i-1)的投影变换矩阵,i=1,2,...,m。
9.如权利要求7所述的系统,其特征在于,对于当前根节点所代表的摄像机的邻接摄像机子集合中的摄像机,其对应的非根节点与当前根节点直接相连通;而对于不属于当前根节点所代表的摄像机的邻接摄像机子集合的摄像机,其对应的非根节点则通过与当前根节点直接相连通的非根节点间接连通。
10.一种全景视频处理设备,其与一摄像机阵列连接,其特征在于,所述全景视频处理设备包括一如权利要求6至9任一项所述的拼接类全景视频的基准视场确定系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110221223 CN102314686B (zh) | 2011-08-03 | 2011-08-03 | 一种拼接类全景视频的基准视场确定方法、系统及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110221223 CN102314686B (zh) | 2011-08-03 | 2011-08-03 | 一种拼接类全景视频的基准视场确定方法、系统及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102314686A true CN102314686A (zh) | 2012-01-11 |
CN102314686B CN102314686B (zh) | 2013-07-17 |
Family
ID=45427824
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110221223 Active CN102314686B (zh) | 2011-08-03 | 2011-08-03 | 一种拼接类全景视频的基准视场确定方法、系统及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102314686B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102984453A (zh) * | 2012-11-01 | 2013-03-20 | 深圳大学 | 利用单摄像机实时生成半球全景视频图像的方法及系统 |
CN103854335A (zh) * | 2012-12-05 | 2014-06-11 | 厦门雅迅网络股份有限公司 | 一种行车记录仪全景视频生成方法 |
CN103942782A (zh) * | 2014-03-31 | 2014-07-23 | Tcl集团股份有限公司 | 一种图像拼接方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6411742B1 (en) * | 2000-05-16 | 2002-06-25 | Adobe Systems Incorporated | Merging images to form a panoramic image |
CN101626513A (zh) * | 2009-07-23 | 2010-01-13 | 深圳大学 | 全景视频生成方法及系统 |
CN102088569A (zh) * | 2010-10-13 | 2011-06-08 | 首都师范大学 | 低空无人飞行器序列图像拼接方法和系统 |
-
2011
- 2011-08-03 CN CN 201110221223 patent/CN102314686B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6411742B1 (en) * | 2000-05-16 | 2002-06-25 | Adobe Systems Incorporated | Merging images to form a panoramic image |
CN101626513A (zh) * | 2009-07-23 | 2010-01-13 | 深圳大学 | 全景视频生成方法及系统 |
CN102088569A (zh) * | 2010-10-13 | 2011-06-08 | 首都师范大学 | 低空无人飞行器序列图像拼接方法和系统 |
Non-Patent Citations (2)
Title |
---|
ROBERTO MARZOTTO ET AL.: "High Resolution Video Mosaicing with Global Alignment", 《PROCEEDINGS OF THE 2004 IEEE COMPUTER SOCIETY CONFERENCE ON COMPUTER VISION AND PATTERN RECOGNITION》, vol. 1, 2 July 2004 (2004-07-02), pages 692 - 698, XP010708627, DOI: doi:10.1109/CVPR.2004.1315099 * |
王伟等: "数字图像拼接技术", 《小型微型计算机系统》, vol. 27, no. 7, 31 July 2006 (2006-07-31), pages 1347 - 1352 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102984453A (zh) * | 2012-11-01 | 2013-03-20 | 深圳大学 | 利用单摄像机实时生成半球全景视频图像的方法及系统 |
CN102984453B (zh) * | 2012-11-01 | 2015-10-28 | 深圳大学 | 利用单摄像机实时生成半球全景视频图像的方法及系统 |
CN103854335A (zh) * | 2012-12-05 | 2014-06-11 | 厦门雅迅网络股份有限公司 | 一种行车记录仪全景视频生成方法 |
CN103854335B (zh) * | 2012-12-05 | 2017-05-03 | 厦门雅迅网络股份有限公司 | 一种行车记录仪全景视频生成方法 |
CN103942782A (zh) * | 2014-03-31 | 2014-07-23 | Tcl集团股份有限公司 | 一种图像拼接方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN102314686B (zh) | 2013-07-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106803271B (zh) | 一种视觉导航无人机的摄像机标定方法及装置 | |
CN110021069B (zh) | 一种基于网格形变的三维模型重建方法 | |
CN109521403B (zh) | 多线激光雷达的参数标定方法及装置、设备及可读介质 | |
Olague et al. | Optimal camera placement for accurate reconstruction | |
Hoppe et al. | Photogrammetric camera network design for micro aerial vehicles | |
JP2012253758A (ja) | 車両視覚システムの較正方法および車両視覚システム | |
CN111639663A (zh) | 多传感器数据融合的方法 | |
CN101996420A (zh) | 信息处理装置、信息处理方法及程序 | |
CN110345921B (zh) | 立体视场视觉测量及垂轴像差和轴向像差校正方法及系统 | |
CN105096329A (zh) | 一种精确校正超广角摄像头图像畸变的方法 | |
JP2009134509A (ja) | モザイク画像生成装置及びモザイク画像生成方法 | |
CN101577004A (zh) | 一种极线矫正方法、装置和系统 | |
CN102314686B (zh) | 一种拼接类全景视频的基准视场确定方法、系统及设备 | |
CN103366356A (zh) | 通用标定系统和方法及其集成成像重构系统和方法 | |
CN103942782A (zh) | 一种图像拼接方法及装置 | |
CN103918003A (zh) | 图像处理装置、方法和程序 | |
CN112991515A (zh) | 一种三维重建方法、装置及相关设备 | |
CN114926598A (zh) | 一种模型生成方法、装置、系统及非易失性存储介质 | |
CN102027509A (zh) | 用于3d路径的方法和装置 | |
CN111508072A (zh) | 地图构建方法、装置、电子设备及存储介质 | |
Dubbelman et al. | Efficient trajectory bending with applications to loop closure | |
CN101082488B (zh) | 异地遥测之图像拼接方法 | |
CN109313819A (zh) | 线路模型实现方法、装置和计算机可读存储介质 | |
CN111504335B (zh) | 地图构建方法、装置、电子设备及存储介质 | |
CN114638867A (zh) | 基于特征提取模块和对偶四元数的点云配准方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |