CN105931186B - 基于相机自动标定和颜色校正的全景视频拼接系统与方法 - Google Patents

基于相机自动标定和颜色校正的全景视频拼接系统与方法 Download PDF

Info

Publication number
CN105931186B
CN105931186B CN201610266342.7A CN201610266342A CN105931186B CN 105931186 B CN105931186 B CN 105931186B CN 201610266342 A CN201610266342 A CN 201610266342A CN 105931186 B CN105931186 B CN 105931186B
Authority
CN
China
Prior art keywords
matrix
overlapping region
image
color
topography
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
CN201610266342.7A
Other languages
English (en)
Other versions
CN105931186A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201610266342.7A priority Critical patent/CN105931186B/zh
Publication of CN105931186A publication Critical patent/CN105931186A/zh
Application granted granted Critical
Publication of CN105931186B publication Critical patent/CN105931186B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformation in the plane of the image
    • G06T3/40Scaling the whole image or part thereof
    • G06T3/4038Scaling the whole image or part thereof for image mosaicing, i.e. plane images composed of plane sub-images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration by the use of histogram techniques
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/60Control of cameras or camera modules
    • H04N23/698Control of cameras or camera modules for achieving an enlarged field of view, e.g. panoramic image capture
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/80Camera processing pipelines; Components thereof
    • 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/10016Video; Image sequence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30232Surveillance

Abstract

本发明公开了一种基于相机自动标定和颜色校正的全景视频拼接系统与方法,其系统包括相机自动标定模块和颜色在线校正模块;相机自动标定模块包括初始单应矩阵计算子模块、重叠区域计算子模块、重叠区域调整子模块、特征匹配子模块和单应矩阵修正子模块;颜色在线校正模块包括重叠区域颜色校正子模块和整图颜色校正子模块;其方法包括相机自动标定步骤S1、颜色在线校正步骤S2和全景图像拼接步骤S3。本发明在自动标定时通过缩放矩阵和偏移矩阵对待拼接图像间的单应矩阵进行修正,使其更加准确;在进行重叠区域的颜色校正时,利用连续三帧间的图像信息对转换关系矩阵进行修正,使得视频中全景图像帧与帧之间的颜色变换更加平滑和流畅。

Description

基于相机自动标定和颜色校正的全景视频拼接系统与方法
技术领域
本发明涉及一种基于相机自动标定和颜色校正的全景视频拼接系统与方法,
背景技术
目前,大量摄像机被安装在政府大楼、军事基地、银行、学校、商场、赛场,随着视频会议、远程教育、机器人导航等领域的兴起,单一相机由于视角限制不能记录下大场景下的所有目标,全景拼接技术的出现满足了各个领域对大视场的需求;当今社会激烈的竞争,人们平时都忙于工作,所以假日选择出游地点大部分依靠的是旅游景点的宣传,但传统的图片不能向游客全方位的展示景点,已经不能满足游客的需求,从而无法吸引游客的眼球;汽车展览时间往往是有限的,人们在网上查询时仅仅能获得汽车局部的图片展示,无法全方位的让人提前全面的了解,而要依靠汽车销售现场不停的解说;买房子时,人们除了要到销售中心看小区的整体布局之外,还要去看房子的内部构造,同时还需要销售人员不停的解说,所以如何全方位的展示目标场景是众多技术研究人员的研究热点。
虽然目前已经具有较为成熟的图片拼接技术,但是一张张独立的拼接图片仍然不能够满足人们的需求,相比之下,人们更倾向于汇聚声,乐,图等的高清全景视频;譬如汽车外部的全景展示,可以从每个角度观看汽车外观,可以在网上构建不落幕的车展;汽车内部的全景展示,可以展示汽车内饰和局部细节;可以让人实现轻松看车、买车。同时高清全景视频在监控领域也具有重要作用;众所周知,监控的运用是十分广泛的,例如交通、银行、企业、小区、地铁、火车和商店等;如果将高清全景视频的技术运动到球赛直播,那么就可以实现实时全景转播,让用户可以在有限的屏幕内一览广阔球场的每一个角落,这提高用户体验感。
但是传统的全景视频拼接技术通常会有不可避免的出现图形畸变和光照引起的同一帧内和帧间的色差问题;因此在拼接过程中我们需要对待拼接图像进行几何校正和颜色校正;传统的几何校正针对不同投影画面间存在的重叠区域,需要对对重叠区域进行几何对准,存在“画面分割”“重叠区域提取”“全局几何对准”“多图像边缘融合”等一系列关键问题,具有较大的工程难度;而且较好的几何校正效果依赖于识别重叠区域的准确性,进而依赖于相机标定的可靠性和相机参数的精确性,传统的相机标定依赖于准确的相机内参和棋盘格的一系列繁琐过程。对于相机参数存在较大误差的情况下,传统的相机标定不能得到较好的几何校正效果;传统的颜色校正采用的自动场景平衡算法倾向于调整曝光,以使得图像的平均颜色是中性灰色,这对于许多安防用途的消费者来说是可以接受的,但是对于球赛直播,看房和旅游类具有特定场景或者较高色彩保真要求的用户是不可接受的;例如,简单的平衡摄影ID图像以使得背景成为白色通常会使得对象的脸、头发和眼睛着色失真,会使得脸变亮,从而冲去特征;可替换地,在这样的场景中朝着中性灰色平衡可能会受对象衣服的着色的影响,并且产生有色背景;相似地,从天空、雪或冰场景中去掉了蓝色调,导致超出可接受水平的着色改变,对于具有暗阴影的图像也是如此;使阴影变亮或者改变阴影颜色可能会影响整个图像,因为使得阴影太亮或者因为不可接受地移动了图像的其它区域颜色;当进行颜色平衡以校正其它颜色的倾斜(skewing)时,会发生类似的不可接受的结果;并且,对于视频来说,传统的颜色校正只是对同一时刻的图像进行了颜色校正,但帧间图像由于拍摄时角度变化和光照不均匀会出现较大的颜色差异。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于相机自动标定和颜色校正的全景视频拼接系统与方法,在自动标定时通过缩放矩阵和偏移矩阵对待拼接图像间的单应矩阵进行修正,使其更加准确;在进行重叠区域的颜色校正时,利用连续三帧间的图像信息对转换关系矩阵进行修正,使得视频中全景图像帧与帧之间的颜色变换更加平滑和流畅。
本发明的目的是通过以下技术方案来实现的:基于相机自动标定和颜色校正的全景视频拼接系统包括相机自动标定模块和颜色在线校正模块;
相机自动标定模块包括初始单应矩阵计算子模块、重叠区域计算子模块、重叠区域调整子模块、特征匹配子模块和单应矩阵修正子模块;
初始单应矩阵计算子模块用于根据相机的内外参数计算待拼接的两幅相邻局部图像间的初始单应矩阵;重叠区域计算子模块用于根据初始单应矩阵计算待拼接的两幅相邻局部图像间重叠区域;重叠区域调整子模块用于对重叠区域进行缩放和平移的调整;特征匹配子模块用于对调整前后的重叠区域进行特征匹配获得偏移矩阵和缩放矩阵;单应矩阵修正子模块用于根据偏移矩阵和缩放矩阵对初始单应矩阵进行修正;
颜色在线校正模块包括重叠区域颜色校正子模块和整图颜色校正子模块;重叠区域颜色校正子模块用于对待拼接的两幅相邻局部图像重叠区域进行颜色校正,整图颜色校正子模块用于对当前帧图像中每一个像素点的像素值进行校正;
重叠区域颜色校正子模块包括初始转换矩阵计算单元、转换矩阵修正单元和重叠区域校正单元;初始转换矩阵计算单元用于通过修正后的单应矩阵计算重叠区域范围,提取重叠区域在两幅图像中的色彩直方图,进行直方图规定化,并计算出转换过程的初始转换矩阵;转换矩阵修正单元用于融入时间信息对初始转换矩阵进行校正;重叠区域校正单元用于根据修正后的转换矩阵对重叠区域的颜色进行校正。
所述的基于相机自动标定和颜色校正的全景视频拼接系统,还包括全景图像拼接模块,所述的全景图像拼接模块用于对通过相机自动标定模块和颜色在线校正模块处理后的全景图像通过光流法拼接成全景视频。
基于相机自动标定和颜色校正的全景视频拼接方法,对每一帧全景图像的处理,均包括相机自动标定步骤S1和颜色在线校正步骤S2:
所述的相机自动标定步骤S1包括以下子步骤:
S11.在当前帧全景图像中,对于任意两幅待拼接局部图像A、B,根据其对应的相机内外参数计算初始单应矩阵H:
S12.根据初始单应矩阵H进行坐标系变换,求出局部图像A和局部图像B的重叠区域;
S13.将重叠区域缩放到统一高度,扩大重叠区域范围;
S14.根据特征匹配算法将局部图像A和局部图像B的重叠区域进行特征匹配;得到准确的参数tx、ty、sx、sy;根据匹配结果求偏移矩阵M和缩放矩阵N:
式中,tx表示重叠图像在x轴平移的坐标;ty表示重叠图像在y轴平移的坐标;sx表示重叠图像在x轴的缩放比例;sy表示重叠图像在y轴的缩放比例;
具体地,这里采用SURF特征点匹配算法对局部图像A和局部图像B的重叠区域进行特征匹配,然后再利用匹配结果计算偏移矩阵和缩放矩阵。
S15.对初始单应矩阵H进行修正,得到修正后的单应矩阵H′:
H′=M·H·N:
所述的颜色在线校正步骤S2包括以下子步骤:
S21.重叠区域颜色校正:求得当前帧重叠区域的颜色转换关系矩阵,并利用连续三帧间的图像信息对转换关系矩阵进行修正,再根据修正后的转换关系矩阵对重叠区域进行颜色校正;
S22.整幅图像的颜色校正:分别统计当前帧全景图像中R,G,B三原色通道图像的直方图,并对每个原色通道分别进行颜色校正。
所述的基于相机自动标定和颜色校正的全景视频拼接方法,还包括全景图像拼接步骤S3:运用光流法对视频每一帧全景图像进行拼接,获得了连续的全景图。
所述的步骤S11包括以下子步骤:
S111.将局部图像A的坐标系投影到世界坐标系得到单应矩阵H1
S112.将局部图像B的坐标系投影到世界坐标系,得到单应矩阵H2
S113.计算得到从局部图像A到局部图像B的初始单应矩阵:
进一步地,H1和H2由相机的内外参数决定。
所述的步骤S21包括以下子步骤:
S211.利用步骤S15中计算得到的单应矩阵H′计算局部图像A和局部图像B重叠区域的范围;
S212.提取局部图像A中重叠区域的色彩直方图G1和局部图像B中重叠区域的第二直方图G2,并进行直方图规定化得到直方图G3
S213.利用直方图G3对局部图像A中的重叠区域进行校正,得到校正后的局部图像A1;其转换关系矩阵R为:
I1表示局部图像B的像素信息;I2表示图像A1的像素信息;其中I1、I2分别是[n 3]的矩阵,n是图像的像素点个数,3表示RGB三原色;也就是说I1、I2均为n行3列的矩阵。
S214.融入时间信息对转换矩阵R进行修正:定义当前帧为T,T-1、T-2分别为前2帧;对转换矩阵R中的I1进行修正得到I1′:
I1′=[IT1 I(T-1)1 I(T-2)1]T
IT1表示当前帧局部图像B的像素信息,I(T-1)1和I(T-2)1分别表示前两帧局部图像B的像素信息;
同时对转换矩阵R中的I2进行修正得到I2′:
I2′=[IT2 I(T-1)2 I(T-2)2]T
式中,IT2表示当前帧校正后的图像A1的像素信息;I(T-2)2和I(T-1)2分别表示前两帧校正后的图像A1的像素信息;
计算修正后的转换矩阵R′为:
R′=[I′T 1I′1]-1I′T 1I′2
S215.利用修正后的矩阵R′对重叠区域进行颜色校正。
由于融入了时间信息,采用三帧连续间的图像信息对转换关系矩阵进行修正,在对重叠区域进行颜色校正,使得视频中全景图像帧与帧之间的颜色变换更加平滑和流畅。
步骤S22中对每个原色通道的校正均包括以下子步骤:
S221.在当前通道的图像中,建立直方图:对于每个像素值xi,统计其像素点个数为S(xi),
并计算该像素值在图像中出现的概率h(xi):
得到直方图H(R):
H(R)=[h(x1),h(x2),…h(xi),…,h(xs)],
式中,n为图像中像素点总数;s为图像中像素值的个数,在H(R)中,横坐标为像素值,纵坐标为该像素值在图像中出现的概率;
S222.定义当前通道图像的偏暗阈值为λ,偏亮阈值为1-λ,求偏暗阈值λ对应的像素值xb
进一步地,λ可以认为是已知量,一般情况下,偏暗阈值λ=0.01,偏亮阈值1-λ=0.99。
进一步地,如果计算得到的xb不是整数,则对xb取整。
求得偏亮阈值1-λ对应的像素值xw
进一步地,如果计算得到的xw不是整数,则对xw取整。
S223.定义当前帧图像为第T帧的图像,T≥2,T-1、T-2帧的图像为前2帧图像,对xb的值进行修正,得到修正后的偏暗阈值λ对应的像素值xb′:
对xw的值进行修正,得到修正后的偏亮阈值1-λ对应的像素值xw′:
进一步地,如果当前帧为第二帧,则xb′和xw′的修正只对第二帧和前一帧求均值,如果当前帧为第一帧,则不需要进行修正。
S224.对整幅图像中的每一个像素点的像素值x进行校正:
(1)如果像素点的像素值x满足xb′≤x≤x′w,则进行线性校正,该像素点校正后的像素值F(x)为:
F(x)=a+bx,
(2)如果当前像素点的像素值x满足0≤x<xb′或者x′w<x≤255,则该像素点校正后的像素值F(x)为:
F(x)=c+t(x/t)d
式中,a,b,c,d为设定的校正参数,t=255。
进一步地,xb′≤x≤x′w时,a和b由(xb′,λ*255)和(x′w,255(1-λ))两个点带入F(x)=a+bx来确定;
在0≤x<xb′时,c和d由(0,0)和(xb′,λ*255)带入F(x)=c+t(x/t)d来确定,且d小于1;
在x′w<x≤255时,c和d由(x′w,255(1-λ))和(255,255)带入F(x)=c+t(x/t)d来确定,且d大于1。
本发明的有益效果是:(1)在自动标定时通过重叠区域的特征点匹配来求得缩放矩阵和偏移矩阵,进而修正得到待拼接图像间的单应矩阵,从而使得相机标定得到的单应矩阵更加准确。
(2)既进行了重叠区域的颜色校正,也进行了整幅图像的颜色校正,校正后的到的图像效果更好。
(3)在进行重叠区域颜色校正时,融入了时间信息,采用连续三帧间的图像信息对转换关系矩阵进行修正,使得视频中全景图像帧与帧之间的颜色变换更加平滑和流畅。
(4)在进行整幅图像的颜色校正时,分别统计当前帧全景图像中三原色通道图像的直方图,并对每个原色通道分别进行颜色校正,且在校正过程中,采用连续三帧图像对偏暗阈值对应的像素值和偏亮阈值对应的像素值进行修正;保持了同帧图像颜色的一致性,进一步加强了帧与帧间颜色变化的流畅性。
附图说明
图1为本发明的系统原理框图;
图2为本发明的方法流程图;
图3为重叠区域颜色校正的流程图;
图4为整幅图像中每个像素点的像素值校正示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
如图1所示,基于相机自动标定和颜色校正的全景视频拼接系统包括相机自动标定模块和颜色在线校正模块;
相机自动标定模块包括初始单应矩阵计算子模块、重叠区域计算子模块、重叠区域调整子模块、特征匹配子模块和单应矩阵修正子模块;
初始单应矩阵计算子模块用于根据相机的内外参数计算待拼接的两幅相邻局部图像间的初始单应矩阵;重叠区域计算子模块用于根据初始单应矩阵计算待拼接的两幅相邻局部图像间重叠区域;重叠区域调整子模块用于对重叠区域进行缩放和平移的调整;特征匹配子模块用于对调整前后的重叠区域进行特征匹配获得偏移矩阵和缩放矩阵;单应矩阵修正子模块用于根据偏移矩阵和缩放矩阵对初始单应矩阵进行修正;
颜色在线校正模块包括重叠区域颜色校正子模块和整图颜色校正子模块;重叠区域颜色校正子模块用于对待拼接的两幅相邻局部图像重叠区域进行颜色校正,整图颜色校正子模块用于对当前帧图像中每一个像素点的像素值进行校正;
重叠区域颜色校正子模块包括初始转换矩阵计算单元、转换矩阵修正单元和重叠区域校正单元;初始转换矩阵计算单元用于通过修正后的单应矩阵计算重叠区域范围,提取重叠区域在两幅图像中的色彩直方图,进行直方图规定化,并计算出转换过程的初始转换矩阵;转换矩阵修正单元用于融入时间信息对初始转换矩阵进行校正;重叠区域校正单元用于根据修正后的转换矩阵对重叠区域的颜色进行校正。
所述的基于相机自动标定和颜色校正的全景视频拼接系统,还包括全景图像拼接模块,所述的全景图像拼接模块用于对通过相机自动标定模块和颜色在线校正模块处理后的全景图像通过光流法拼接成全景视频。
如图2所示,基于相机自动标定和颜色校正的全景视频拼接方法,对每一帧全景图像的处理,均包括相机自动标定步骤S1和颜色在线校正步骤S2:
所述的相机自动标定步骤S1包括以下子步骤:
S11.在当前帧全景图像中,对于任意两幅待拼接局部图像A、B,根据其对应的相机内外参数计算初始单应矩阵H:
S12.根据初始单应矩阵H进行坐标系变换,求出局部图像A和局部图像B的重叠区域;
S13.将重叠区域缩放到统一高度,扩大重叠区域范围;
S14.根据特征匹配算法将局部图像A和局部图像B的重叠区域进行特征匹配;得到准确的参数tx、ty、sx、sy;根据匹配结果求偏移矩阵M和缩放矩阵N:
式中,tx表示重叠图像在x轴平移的坐标;ty表示重叠图像在y轴平移的坐标;sx表示重叠图像在x轴的缩放比例;sy表示重叠图像在y轴的缩放比例;
具体地,这里采用SURF特征点匹配算法对局部图像A和局部图像B的重叠区域进行特征匹配,然后再利用匹配结果计算偏移矩阵和缩放矩阵。
S15.对初始单应矩阵H进行修正,得到修正后的单应矩阵H′:
H′=M·H·N:
所述的颜色在线校正步骤S2包括以下子步骤:
S21.重叠区域颜色校正:求得当前帧重叠区域的颜色转换关系矩阵,并利用连续三帧间的图像信息对转换关系矩阵进行修正,再根据修正后的转换关系矩阵对重叠区域进行颜色校正;
S22.整幅图像的颜色校正:分别统计当前帧全景图像中R,G,B三原色通道图像的直方图,并对每个原色通道分别进行颜色校正。
所述的基于相机自动标定和颜色校正的全景视频拼接方法,还包括全景图像拼接步骤S3:运用光流法对视频每一帧全景图像进行拼接,获得了连续的全景图,即全景视频。
所述的步骤S11包括以下子步骤:
S111.将局部图像A的坐标系投影到世界坐标系得到单应矩阵H1
S112.将局部图像B的坐标系投影到世界坐标系,得到单应矩阵H2
S113.计算得到从局部图像A到局部图像B的初始单应矩阵:
进一步地,H1和H2由相机的内外参数决定。
如图3所示,所述的步骤S21包括以下子步骤:
S211.利用步骤S15中计算得到的单应矩阵H′计算局部图像A和局部图像B重叠区域的范围;
S212.提取局部图像A中重叠区域的色彩直方图G1和局部图像B中重叠区域的第二直方图G2,并进行直方图规定化得到直方图G3
S213.利用直方图G3对局部图像A中的重叠区域进行校正,得到校正后的局部图像A1;其转换关系矩阵R为:
I1表示局部图像B的像素信息;I2表示图像A1的像素信息;其中I1、I2分别是[n 3]的矩阵,n是图像的像素点个数,3表示RGB三原色;
S214.融入时间信息对转换矩阵R进行修正:定义当前帧为T,T-1、T-2分别为前2帧;对转换矩阵R中的I1进行修正得到I1′:
I1′=[IT1 I(T-1)1 I(T-2)1]T
IT1表示当前帧局部图像B的像素信息,I(T-1)1利I(T-2)1分别表示前两帧局部图像B的像素信息;
同时对转换矩阵R中的I2进行修正得到I2′:
I2′=[IT2 I(T-1)2 I(T-2)2]T
式中,IT2表示当前帧校正后的图像Al的像素信息;I(T-2)2和I(T-1)2分别表示前两帧校正后的图像A1的像素信息;
计算修正后的转换矩阵R′为:
S215.利用修正后的矩阵R′对重叠区域进行颜色校正。
由于融入了时间信息,采用三帧连续间的图像信息对转换关系矩阵进行修正,在对重叠区域进行颜色校正,使得视频中全景图像帧与帧之间的颜色变换更加平滑和流畅。
步骤S22中对每个原色通道的校正均包括以下子步骤:
S221.在当前通道的图像中,建立直方图:对于每个像素值xi,统计其像素点个数为S(xi),
并计算该像素值在图像中出现的概率h(xi):
得到直方图H(R):
H(R)=[h(x1),h(x2),…h(xi),…,h(xs)],
式中,n为图像中像素点总数;s为图像中像素值的个数,在H(R)中,横坐标为像素值,纵坐标为该像素值在图像中出现的概率;
S222.定义当前通道图像的偏暗阈值为λ,偏亮阈值为1-λ,求偏暗阈值λ对应的像素值xb
进一步地,λ可以认为是已知量,一般情况下,偏暗阈值λ=0.01,偏亮阈值1-λ=0.99。
进一步地,如果计算得到的xb不是整数,则对xb取整。
求得偏亮阈值1-λ对应的像素值xw
进一步地,如果计算得到的xw不是整数,则对xw取整。
S223.定义当前帧图像为第T帧的图像,T≥2,T-1、T-2帧的图像为前2帧图像,对xb的值进行修正,得到修正后的偏暗阈值λ对应的像素值xb′:
对xw的值进行修正,得到修正后的偏亮阈值1-λ对应的像素值xw′:
进一步地,如果当前帧为第二帧,则xb′和xw′的修正只对第二帧和前一帧求均值,如果当前帧为第一帧,则不需要进行修正。
S224.如图4所示,对整幅图像中的每一个像素点的像素值x进行校正:
(1)如果像素点的像素值x满足xb′≤x≤x′w,则进行线性校正,该像素点校正后的像素值F(x)为:
F(x)=a+bx,
(2)如果当前像素点的像素值x满足0≤x<xb′或者x′w<x≤255,则该像素点校正后的像素值F(x)为:
F(x)=c+t(x/t)d
式中,a,b,c,d为设定的校正参数,t=255。
进一步地,xb′≤x≤x′w时,a和h由(xb′,λ*255)和(x′w,255(1-λ))两个点带入F(x)=a+bx来确定;
在0≤x<xb′时,c和d由(0,0)和(xb′,λ*255)带入F(x)=c+t(x/t)d来确定,且d小于1;
在x′w<x≤255时,c和d由(x′w,255(1-λ))和(255,255)带入F(x)=c+t(x/t)d来确定,且d大于1。
进一步地,在进行全景视频拼接的过程中,一般采用三台不同方位的摄像机进行图像数据采集,三台摄像机分别输入每一帧全景图像待拼接的左图、中图和右图;每一帧全景图像的拼接过程,实际上就是将左图、右图分别与中图拼接;左图与中图可以看作是两幅待拼接局部图像;右图与中图也可以看作是两幅待拼接局部图像;
在本申请中,可以将左图看作上述的局部图像A,中图看作上述的局部图像B,将其按照步骤S1进行相机自动标定,再按照步骤S21进行重叠区域的颜色校正。右图与中图重叠区域颜色校正同理,将其分别看作局部图像A和局部图像B,按照步骤S1进行相机自动标定,按照步骤S21进行重叠区域颜色校正即可。重叠区域颜色校正完成后,再按照步骤S22对整幅全景图像进行颜色校正。

Claims (6)

1.基于相机自动标定和颜色校正的全景视频拼接系统,其特征在于:包括相机自动标定模块和颜色在线校正模块;
相机自动标定模块包括初始单应矩阵计算子模块、重叠区域计算子模块、重叠区域调整子模块、特征匹配子模块和单应矩阵修正子模块;
初始单应矩阵计算子模块用于根据相机的内外参数计算待拼接的两幅相邻局部图像间的初始单应矩阵;重叠区域计算子模块用于根据初始单应矩阵计算待拼接的两幅相邻局部图像间重叠区域;重叠区域调整子模块用于对重叠区域进行缩放和平移的调整;特征匹配子模块用于对调整前后的重叠区域进行特征匹配获得偏移矩阵和缩放矩阵;单应矩阵修正子模块用于根据偏移矩阵和缩放矩阵对初始单应矩阵进行修正;
颜色在线校正模块包括重叠区域颜色校正子模块和整图颜色校正子模块;重叠区域颜色校正子模块用于对待拼接的两幅相邻局部图像重叠区域进行颜色校正,整图颜色校正子模块用于对当前帧图像中每一个像素点的像素值进行校正;
重叠区域颜色校正子模块包括初始转换矩阵计算单元、转换矩阵修正单元和重叠区域校正单元;初始转换矩阵计算单元用于通过修正后的单应矩阵计算重叠区域范围,提取重叠区域在两幅图像中的色彩直方图,进行直方图规定化,并计算出转换过程的初始转换矩阵;转换矩阵修正单元用于融入时间信息对初始转换矩阵进行校正;重叠区域校正单元用于根据修正后的转换矩阵对重叠区域的颜色进行校正。
2.根据权利要求1所述的基于相机自动标定和颜色校正的全景视频拼接系统,其特征在于:还包括全景图像拼接模块,所述的全景图像拼接模块用于对通过相机自动标定模块和颜色在线校正模块处理后的全景图像通过光流法拼接成全景视频。
3.基于相机自动标定和颜色校正的全景视频拼接方法,其特征在于:对每一帧全景图像的处理,均包括相机自动标定步骤S1和颜色在线校正步骤S2:
所述的相机自动标定步骤S1包括以下子步骤:
S11.在当前帧全景图像中,对于任意两幅待拼接局部图像A、B,根据其对应的相机内外参数计算初始单应矩阵H:
S12.根据初始单应矩阵H进行坐标系变换,求出局部图像A和局部图像B的重叠区域;
S13.将重叠区域缩放到统一高度,扩大重叠区域范围;
S14.根据特征匹配算法将局部图像A和局部图像B的重叠区域进行特征匹配;得到准确的参数tx、ty、sx、sy;根据匹配结果求偏移矩阵M和缩放矩阵N:
式中,tx表示重叠图像在x轴平移的坐标;ty表示重叠图像在y轴平移的坐标;sx表示重叠图像在x轴的缩放比例;sy表示重叠图像在y轴的缩放比例;
S15.对初始单应矩阵H进行修正,得到修正后的单应矩阵H′:
H′=M·H·N;
所述的颜色在线校正步骤S2包括以下子步骤:
S21.重叠区域颜色校正:求得当前帧重叠区域的颜色转换关系矩阵,并利用连续三帧
间的图像信息对转换关系矩阵进行修正,再根据修正后的转换关系矩阵对重叠区域
进行颜色校正;
S22.整幅图像的颜色校正:分别统计当前帧全景图像中R,G,B三原色通道图像的直方
图,并对每个原色通道分别进行颜色校正;
所述的步骤S21包括以下子步骤:
S211.利用步骤S15中计算得到的单应矩阵H′计算局部图像A和局部图像B重叠区域的范围;
S212.提取局部图像A中重叠区域的色彩直方图G1和局部图像B中重叠区域的第二直方图G2,并进行直方图规定化得到直方图G3
S213.利用直方图G3对局部图像A中的重叠区域进行校正,得到校正后的局部图像A1;其转换关系矩阵R为:
I1表示局部图像B的像素信息;I2表示图像A1的像素信息;
S214.融入时间信息对转换矩阵R进行修正:定义当前帧为T,T-1、T-2分别为前2帧;对转换矩阵R中的I1进行修正得到I′1
I′1=[IT1 I(T-1)1 I(T-2)1]T
IT1表示当前帧局部图像B的像素信息,I(T-1)1和I(T-2)1分别表示前两帧局部图像B的像素信息;同时对转换矩阵R中的I2进行修正得到I′2
I′2=[IT2 I(T-1)2 I(T-2)2]T
式中,IT2表示当前帧校正后的图像A1的像素信息;I(T-2)2和I(T-1)2分别表示前两帧校正后的图像A1的像素信息;
计算修正后的转换矩阵R′为:
S215.利用修正后的矩阵R′对重叠区域进行颜色校正。
4.根据权利要求3所述的基于相机自动标定和颜色校正的全景视频拼接方法,其特征在于:还包括全景图像拼接步骤S3:运用光流法对视频每一帧全景图像进行拼接,获得了连续的全景图,即全景视频。
5.根据权利要求3所述的基于相机自动标定和颜色校正的全景视频拼接方法,其特征在于:所述的步骤S11包括以下子步骤:
S111.将局部图像A的坐标系投影到世界坐标系得到单应矩阵H1
S112.将局部图像B的坐标系投影到世界坐标系,得到单应矩阵H2
S113.计算得到从局部图像A到局部图像B的初始单应矩阵:
6.根据权利要求3所述的基于相机自动标定和颜色校正的全景视频拼接方法,其特征在于:步骤S22中对每个原色通道的校正均包括以下子步骤:
S221.在当前通道的图像中,建立直方图:对于每个像素值xi,统计其像素点个数为S(xi),并计算该像素值在图像中出现的概率h(xi):
得到直方图H(R):
H(R)=[h(x1),h(x2),…h(xi),…,h(xs)],
式中,n为图像中像素点总数;s为图像中像素值的个数,在H(R)中,横坐标为像素值,纵坐标为该像素值在图像中出现的概率;
S222.定义当前通道图像的偏暗阈值为λ,偏亮阈值为1-λ,求偏暗阈值λ对应的像素值xb
求得偏亮阈值1-λ对应的像素值xw
S223.定义当前帧图像为第T帧的图像,T≥2,T-1、T-2帧的图像为前2帧图像,对xb的值进行修正,得到修正后的偏暗阈值λ对应的像素值xb′:
对xw的值进行修正,得到修正后的偏亮阈值1-λ对应的像素值xw′:
S224.对整幅图像中的每一个像素点的像素值x进行校正:
(1)如果像素点的像素值x满足x′b≤x≤x′w,则进行线性校正,该像素点校正后的像素值F(x)为:
F(x)=a+bx,
(2)如果当前像素点的像素值x满足0≤x<x′b或者x′w<x≤255,则该像素点校正后的像素值F(x)为:
F(x)=c+t(x/t)d
式中,a,b,c,d为设定的校正参数,t=255。
CN201610266342.7A 2016-04-26 2016-04-26 基于相机自动标定和颜色校正的全景视频拼接系统与方法 Active CN105931186B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610266342.7A CN105931186B (zh) 2016-04-26 2016-04-26 基于相机自动标定和颜色校正的全景视频拼接系统与方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610266342.7A CN105931186B (zh) 2016-04-26 2016-04-26 基于相机自动标定和颜色校正的全景视频拼接系统与方法

Publications (2)

Publication Number Publication Date
CN105931186A CN105931186A (zh) 2016-09-07
CN105931186B true CN105931186B (zh) 2019-04-02

Family

ID=56837302

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610266342.7A Active CN105931186B (zh) 2016-04-26 2016-04-26 基于相机自动标定和颜色校正的全景视频拼接系统与方法

Country Status (1)

Country Link
CN (1) CN105931186B (zh)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106657983B (zh) * 2016-11-16 2019-02-22 深圳六滴科技有限公司 全景摄像机的参数测试方法和装置
CN106600529A (zh) * 2016-11-28 2017-04-26 北京暴风魔镜科技有限公司 一种全焦全景图像的获取方法和装置
CN106780312B (zh) * 2016-12-28 2020-02-18 南京师范大学 基于sift匹配的图像空间与地理场景自动映射方法
CN106878628A (zh) * 2017-02-23 2017-06-20 珠海数字动力科技股份有限公司 一种通过摄像头进行视频拼接的方法
CN107451952B (zh) * 2017-08-04 2020-11-03 追光人动画设计(北京)有限公司 一种全景视频的拼接融合方法、设备以及系统
EP3487162B1 (en) * 2017-11-16 2021-03-17 Axis AB Method, device and camera for blending a first and a second image having overlapping fields of view
CN108171655A (zh) * 2017-12-27 2018-06-15 深圳普思英察科技有限公司 基于单目摄像头的倒车影像拼接方法和装置
CN108363519B (zh) * 2017-12-29 2020-03-31 中国科学院自动化研究所 分布式红外视觉检测与投影融合自动矫正的触控显示系统
CN108307200B (zh) * 2018-01-31 2020-06-09 深圳积木易搭科技技术有限公司 一种在线视频拼接方法系统
CN110689476A (zh) * 2018-07-04 2020-01-14 杭州海康威视数字技术股份有限公司 全景图像拼接方法及其装置、可读存储介质和电子设备
CN109636728A (zh) * 2018-12-24 2019-04-16 天津天地伟业信息系统集成有限公司 一种图像拼接的颜色校正方法
CN111047510B (zh) 2019-12-17 2023-02-14 大连理工大学 一种基于标定的大视场角图像实时拼接方法
CN113079275B (zh) * 2020-01-03 2022-05-17 南京大学 一种基于优化直方图匹配的阵列相机色彩校正方法
CN111915587B (zh) * 2020-07-30 2024-02-02 北京大米科技有限公司 视频处理方法、装置、存储介质和电子设备
CN113096043B (zh) * 2021-04-09 2023-02-17 杭州睿胜软件有限公司 图像处理方法及装置、电子设备和存储介质
CN115348398A (zh) * 2021-05-13 2022-11-15 影石创新科技股份有限公司 视频处理方法、装置、计算机设备和存储介质
CN113905183B (zh) * 2021-08-25 2024-04-16 珠海全志科技股份有限公司 宽动态范围图像的色差校正方法及装置
CN114743274B (zh) * 2022-06-09 2022-09-09 深圳市纬亚森科技有限公司 基于物联网的多媒体远程教学管理系统及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2808805A1 (en) * 2013-05-30 2014-12-03 Thomson Licensing Method and apparatus for displaying metadata on a display and for providing metadata for display
CN104392416A (zh) * 2014-11-21 2015-03-04 中国电子科技集团公司第二十八研究所 一种运动场景的视频拼接方法
CN105141920A (zh) * 2015-09-01 2015-12-09 电子科技大学 一种360度全景视频拼接系统
CN105389787A (zh) * 2015-09-30 2016-03-09 华为技术有限公司 一种全景图像拼接方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8238610B2 (en) * 2008-12-18 2012-08-07 University Of Central Florida Research Foundation, Inc. Homography-based passive vehicle speed measuring
KR101697512B1 (ko) * 2011-12-15 2017-01-19 한국전자통신연구원 영상 정합 장치 및 방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2808805A1 (en) * 2013-05-30 2014-12-03 Thomson Licensing Method and apparatus for displaying metadata on a display and for providing metadata for display
CN104392416A (zh) * 2014-11-21 2015-03-04 中国电子科技集团公司第二十八研究所 一种运动场景的视频拼接方法
CN105141920A (zh) * 2015-09-01 2015-12-09 电子科技大学 一种360度全景视频拼接系统
CN105389787A (zh) * 2015-09-30 2016-03-09 华为技术有限公司 一种全景图像拼接方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A system for real-time panorama generation and display in;Wai-Kwan Tang, et al.;《IEEE Transactions on Multimedia》;20051231;第7卷(第2期);第455-461页
基于角点检测图像配准的一种新算法;周鹏,等;《中国科学技术大学学报》;20020831;第32卷(第4期);第280-292页

Also Published As

Publication number Publication date
CN105931186A (zh) 2016-09-07

Similar Documents

Publication Publication Date Title
CN105931186B (zh) 基于相机自动标定和颜色校正的全景视频拼接系统与方法
CN105827975B (zh) 一种用于全景视频拼接的颜色在线校正方法
US8717352B2 (en) Tracing-type stereo display apparatus and tracing-type stereo display method
CN105141920B (zh) 一种360度全景视频拼接系统
Itoh et al. Semi-parametric color reproduction method for optical see-through head-mounted displays
CN104702928B (zh) 修正图像交叠区的方法、记录介质以及执行装置
CN109598673A (zh) 图像拼接方法、装置、终端及计算机可读存储介质
CN103019643A (zh) 一种即插即用的大屏幕投影自动校正与拼接显示方法
CN107077630A (zh) 半透明标记、用于合成及探测半透明标记的方法、透明标记、以及用于合成及探测透明标记的方法
CN108154514A (zh) 图像处理方法、装置及设备
US11736674B2 (en) Dynamic convergence adjustment in augmented reality headsets
CN104427318A (zh) 校正图像交叠区域的方法、记录介质和执行装置
CN107241610A (zh) 一种基于增强现实的虚拟内容插入系统和方法
US10244216B2 (en) Projector, video display device, and video display method
Hamasaki et al. Hysar: Hybrid material rendering by an optical see-through head-mounted display with spatial augmented reality projection
CN105657268A (zh) 一种基于多分辨率的多视点视频拼接融合算法
CN109996048A (zh) 一种基于结构光的投影校正方法及其系统
CN208506731U (zh) 图像展示系统
RU2735066C1 (ru) Способ отображения широкоформатного объекта дополненной реальности
CN110290313B (zh) 一种引导自动对焦设备失焦的方法
CN113792564B (zh) 基于不可见投影二维码的室内定位方法
CN103903274B (zh) 一种小半径大曲率畸变曲面投影校正融合的方法
KR20170044319A (ko) 헤드 마운트 디스플레이의 시야 확장 방법
CN114022562A (zh) 一种保持行人完整性的全景视频拼接方法和装置
EP4150560A1 (en) Single image 3d photography with soft-layering and depth-aware inpainting

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant