CN108876719A - 基于虚拟相机模型的车辆全景图像拼接外部参数估计方法 - Google Patents
基于虚拟相机模型的车辆全景图像拼接外部参数估计方法 Download PDFInfo
- Publication number
- CN108876719A CN108876719A CN201810271061.XA CN201810271061A CN108876719A CN 108876719 A CN108876719 A CN 108876719A CN 201810271061 A CN201810271061 A CN 201810271061A CN 108876719 A CN108876719 A CN 108876719A
- Authority
- CN
- China
- Prior art keywords
- camera
- coordinate
- external parameter
- vehicle
- panoramic
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformation in the plane of the image
- G06T3/40—Scaling the whole image or part thereof
- G06T3/4038—Scaling the whole image or part thereof for image mosaicing, i.e. plane images composed of plane sub-images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
- G06T7/85—Stereo camera calibration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30204—Marker
- G06T2207/30208—Marker matrix
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30244—Camera pose
Abstract
本发明涉及基于虚拟相机模型的车辆全景图像拼接外部参数估计方法,包括步骤:(1)针对具体车辆建立虚拟相机模型,计算出虚拟相机模型的参数矩阵;将摄像头安装在车辆周围,在车身周围地面区域放置多个标记点;(2)每个摄像头拍摄一组含有标记点的图像,找出其对应的多个标记点在全景视野范围的坐标,并通过待求的外部参数矩阵将该坐标转化成当前相机坐标系下的坐标;(3)对原始图像进行畸变校正,找出标记点的图像坐标,并通过相机成像模型将其转化为对应相机坐标系下的坐标;将步骤(2)、(3)的坐标建立方程组,求解出每个摄像头的外部参数。本发明通过建立方程组得到每个摄像头外部参数的解析解简化了外部参数估计的步骤。
Description
技术领域
本发明属于汽车电子技术领域,具体涉及基于虚拟相机模型的车辆全景图像拼接技术。
背景技术
车辆全景图像拼接是指在车辆的周围安装多个广角或鱼眼摄像头,通常是前后左右共4个摄像头,采用图像处理的方法将多个摄像头的视频图像无缝地拼接成一幅全景视频图像。驾驶员通过全景视频图像能够直观并实时地观察车身周围的情况,极大地提高了驾驶员在复杂环境下行车的安全性。因此全景图像拼接方法的研究具有直接的现实意义。
基于虚拟相机模型的全景图像拼接方法主要是假设在车辆的正上方存在一个虚拟相机,全景图像直接由该虚拟相机生成。通过分析虚拟相机和周围摄像头的相互位置关系,利用图像变换的方法实现周围摄像头图像到虚拟相机图像映射,经过图像融合技术形成最终的全景视频。相对于对周围图像进行直接拼接的方式,基于虚拟相机模型的全景图像拼接方法具有操作上的灵活性。在基于虚拟相机模型的拼接方法中,虚拟相机和周围相机的位置关系被称为外部参数,其求解过程被称为外部参数校正或外部参数估计。参数的准确性直接影响了最终全景拼接的效果。目前传统方法是采用图像分析通过优化的方法计算该外部参数,操作复杂且稳定性不高。
发明内容
为了解决现有基于虚拟相机模型的车辆全景图像拼接技术所存在的问题,本发明提出了基于虚拟相机模型的车辆全景图像拼接外部参数估计方法,该方法在车辆周围地面区域放置标记点,通过分析摄像头校正图像与全景图像中标记点的位置关系,直接推导出外部参数的解析解,操作简便且可靠、有效地进行外部参数估计,能够改善全景拼接的效率与性能。
本发明所采用的技术方案如下:基于虚拟相机模型的车辆全景图像拼接外部参数估计方法,包括以下步骤:
(1)确定车辆全景视野范围和全景图像的尺寸,针对具体车辆建立虚拟相机模型,根据全景视野范围、全景图像的尺寸和虚拟相机的高度计算出虚拟相机模型的参数矩阵;将校正好的摄像头安装在车辆周围,在车身周围地面区域放置多个标记点;
(2)每个摄像头拍摄一组含有标记点的图像,用于计算该摄像头与虚拟相机的位置关系,即外部参数;针对每一个摄像头,找出其对应的多个标记点在全景视野范围的坐标,并通过待求的外部参数矩阵将该坐标转化成当前相机坐标系下的坐标;
(3)针对每一个摄像头,对含有标记点的原始图像进行畸变校正,在校正图像中找出标记点的图像坐标,并通过相机成像模型将标记点的图像坐标转化为对应相机坐标系下的坐标;将步骤(2)与步骤(3)得到的坐标建立方程组,求解出每个摄像头的外部参数。
优选地,步骤(2)中,全景视野范围的坐标转化为当前相机坐标系下的坐标包括三次转化:首先是全景视野范围的坐标转化为全景图像齐次坐标,其次是全景图像齐次坐标转化为全景相机坐标系下的坐标,最后利用外部参数矩阵转化为当前相机坐标系下的坐标。
本发明的一个优选实施例中,在车身周围地面区域放置的标记点为12个或16个。待求的外部参数矩阵为一个3×3的矩阵,具有8个未知参数且第3行第3列的元素为1。在车辆的前后左右共安装四个摄像头,每个摄像头至少能够拍摄到4个标记点,每个标记点产生2个方程,4个标记点一共产生8个方程;将8个方程组成线性方程组,求解出所述待求的外部参数矩阵的8个未知参数。
本发明相对于现有技术具有如下的有益效果:
1、首先,本发明通过标记点坐标分析分别地对每个摄像头外部参数进行估计,能够提高全景图像拼接的灵活性。其次,本发明提出的外部参数估计方法具有解析解,能够极大地提高计算效率以及图像拼接的准确性。
2、本发明巧妙地利用摄像头的校正图像作为原始图像和全景视野坐标转化的桥梁,并通过建立方程组得到每个摄像头外部参数的解析解,极大低简化了外部参数估计的步骤,同时提高其精确性,从而提高了基于虚拟相机模型的车辆全景图像拼接效率和性能。
附图说明
图1为本发明外部参数估计方法的流程框架示意图;
图2为本发明标记点位置示意图,其中(a)示意了12个标记点的情形,(b)示意了16个标记点的情形;
图3为本发明虚拟相机与右侧相机坐标转换示意图。
具体实施方式
下面结合实施例及附图对本发明作进一步的详细说明,但本发明的实施方式不限于此。
实施例
本实施例中,基于虚拟相机模型的车辆全景图像拼接外部参数估计方法,其流程框架如图1所示。具体实现包括以下步骤:
步骤1、根据实际应用情况确定车辆全景视野范围和全景图像的尺寸,假设全景图像视野范围为A×B的矩形区域,全景图像的尺寸为a×b。指定虚拟相机的高度为H,该值对计算结果影响不大,但一般取值大小与全景视野范围在一个数量级。针对具体车辆建立虚拟相机模型,根据全景视野范围、全景图像的尺寸和虚拟相机的高度计算出虚拟相机模型的参数矩阵为:
其中u0=(a-1)×0.5,v0=(b-1)×0.5。
根据具体车辆大小选择摄像头的个数以满足全景图像拼接的要求,一般小型轿车在前后左右共安装4个鱼眼摄像头即可,下面以4个摄像头的情况进行说明。首先将鱼眼摄像头通过鱼眼畸变校正的方法得到相机参数和畸变系数,然后将校正好的摄像头安装在车辆周围合适位置。在车辆周围地面区域放置12个标记点,如图2(a)所示,这些标记点必须在全景视野范围内且每个鱼眼摄像头至少能够拍摄到4个标记点,有部分标记点为公共点。也可以在车辆周围放置16个标记点,如图2(b)所示,这样更加保证每个鱼眼摄像头能够至少拍摄4个标记点。
步骤2、摄像头安装好并且标记点放置完成后,每个摄像头拍摄一组含有标记点的图像,用于后续计算相应摄像头与虚拟相机的位置关系,即外部参数。
对每个摄像头分别单独计算其外部参数。针对每一个摄像头,找出其对应4个标记点在全景视野范围的坐标,并通过其待求的外部参数矩阵将该坐标转为当前相机(即对应摄像头)坐标系下的坐标。
由于每个摄像头与虚拟相机的位置关系是单独计算的,下面以车辆右方摄像头为例来说明:以全景视野区域的左上角为原点建立直角坐标系,水平向右为X轴,水平向下为Y轴。记录右方摄像头对应的4个标记点在全景视野区域的坐标。如图3所示,假设其中一个点P的坐标为(x,y),根据全景视频与全景图像的比例关系,可以得到该点在全景图像中的齐次坐标为:
根据相机成像模型,可以得到该点在全景相机坐标系下的坐标为:
其中inv(K)表示虚拟相机参数矩阵K的逆矩阵。假设待求的外部参数矩阵为Rt,即右侧相机相对于虚拟相机的旋转平移矩阵。则P点在右侧相机坐标系下的坐标可以表示为:
步骤3、针对每一个摄像头,对含有标记点的原始鱼眼图像进行畸变校正,在校正图像中找出标记点的图像坐标,并通过相机成像模型将4个标记点的图像坐标转化为对应相机坐标系下的坐标;将步骤2与步骤3得到的坐标建立方程组,求解出每个摄像头的外部参数。
本步骤对含有标记点的原始鱼眼图像进行校正后,校正图像中仍然能找到4个标记点的位置,且不能有3个以上的点共线。
继续以车辆右方摄像头为例进行说明:从右侧摄像头的畸变图像中得到P点在右侧相机坐标系下的坐标。首先对原始的鱼眼图像进行畸变校正,然后在校正后的图像中找出P点的像素坐标,记作为(p1,p2),根据相机成像模型,可以得到以下关系:
其中,KA为右侧摄像机模型的参数矩阵,经过通用的摄像头校正即可以得到,其表示为:
现假设待求的外部参数矩阵的形式为:
将(3)式代入(4)式有:
s3·p=KA·Rt·v, (7)
其中s3=r31v1+r32v2+v3,将(7)式进行展开得到:
进一步化简得:
将上面等式写成矩阵的形式:
上面矩阵方程中,只有外部参数矩阵Rt(表示成列向量)为未知参数,且Rt有8个未知数。如(10)式所示,每个点能够提供2个方程,因此利用4个点可以求解出右侧摄像头相对于虚拟相机的外部参数矩阵Rt。
也就是说,按同样的方式对其他三个标记点建立相应方程,每个标记点产生2个方程,四个标记点一共产生8个;将8个方程组成线性方程组,假设其形式为:
A[r11 r12 r13 r21 r22 r23 r31 r32]T=b, (11)
其中A为系数矩阵,可以直接计算出解析解:
[r11 r12 r13 r21 r22 r23 r31 r32]T=inv(A)·b. (12)
对解向量进行重排即得到外部参数矩阵Rt。
得到外部参数矩阵Rt,则可以实现右侧畸变图像到全景图像的映射。采用相同的方式,可以实现周围其他摄像头图像到全景图像的转换,通过现有的图像融合和光照归一化技术可以实现全景图像的生成。
以上所述,仅为本发明专利发明较佳/优选的实施方式,但发明专利的保护范围不局限于此,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围内。
Claims (6)
1.基于虚拟相机模型的车辆全景图像拼接外部参数估计方法,其特征在于,包括以下步骤:
(1)确定车辆全景视野范围和全景图像的尺寸,针对具体车辆建立虚拟相机模型,根据全景视野范围、全景图像的尺寸和虚拟相机的高度计算出虚拟相机模型的参数矩阵;将校正好的摄像头安装在车辆周围,在车身周围地面区域放置多个标记点;
(2)每个摄像头拍摄一组含有标记点的图像,用于计算该摄像头与虚拟相机的位置关系,即外部参数;针对每一个摄像头,找出其对应的多个标记点在全景视野范围的坐标,并通过待求的外部参数矩阵将该坐标转化成当前相机坐标系下的坐标;
(3)针对每一个摄像头,对含有标记点的原始图像进行畸变校正,在校正图像中找出标记点的图像坐标,并通过相机成像模型将标记点的图像坐标转化为对应相机坐标系下的坐标;将步骤(2)与步骤(3)得到的坐标建立方程组,求解出每个摄像头的外部参数。
2.根据权利要求1所述的车辆全景图像拼接外部参数估计方法,其特征在于,假设全景图像视野范围为A×B的矩形区域,全景图像的尺寸为a×b,虚拟相机的高度为H,虚拟相机模型的参数矩阵为:
其中u0=(a-1)×0.5,v0=(b-1)×0.5。
3.根据权利要求1所述的车辆全景图像拼接外部参数估计方法,其特征在于,步骤(2)中,全景视野范围的坐标转化为当前相机坐标系下的坐标包括三次转化:首先是全景视野范围的坐标转化为全景图像齐次坐标,其次是全景图像齐次坐标转化为全景相机坐标系下的坐标,最后利用外部参数矩阵转化为当前相机坐标系下的坐标。
4.根据权利要求1所述的车辆全景图像拼接外部参数估计方法,其特征在于,在车身周围地面区域放置的标记点为12个或16个。
5.根据权利要求1所述的车辆全景图像拼接外部参数估计方法,其特征在于,所述待求的外部参数矩阵为一个3×3的矩阵,具有8个未知参数且第3行第3列的元素为1。
6.根据权利要求5所述的车辆全景图像拼接外部参数估计方法,其特征在于,在车辆的前后左右共安装四个摄像头,每个摄像头至少能够拍摄到4个标记点,每个标记点产生2个方程,4个标记点一共产生8个方程;将8个方程组成线性方程组,求解出所述待求的外部参数矩阵的8个未知参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810271061.XA CN108876719B (zh) | 2018-03-29 | 2018-03-29 | 基于虚拟相机模型的车辆全景图像拼接外部参数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810271061.XA CN108876719B (zh) | 2018-03-29 | 2018-03-29 | 基于虚拟相机模型的车辆全景图像拼接外部参数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108876719A true CN108876719A (zh) | 2018-11-23 |
CN108876719B CN108876719B (zh) | 2022-07-26 |
Family
ID=64326421
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810271061.XA Active CN108876719B (zh) | 2018-03-29 | 2018-03-29 | 基于虚拟相机模型的车辆全景图像拼接外部参数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108876719B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109523489A (zh) * | 2018-11-28 | 2019-03-26 | 广州大学 | 一种生成俯视无畸变倒车影像的方法 |
CN110058263A (zh) * | 2019-04-16 | 2019-07-26 | 广州大学 | 一种车辆行驶过程中的物体定位方法 |
CN113840755A (zh) * | 2019-05-21 | 2021-12-24 | 康蒂-特米克微电子有限公司 | 生成车辆周围环境图像的方法及生成车辆周围环境图像的装置 |
CN114549666A (zh) * | 2022-04-26 | 2022-05-27 | 杭州蓝芯科技有限公司 | 一种基于agv的环视图像拼接标定方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009017332A1 (en) * | 2007-07-29 | 2009-02-05 | Nanophotonics Co., Ltd. | Methods of obtaining panoramic images using rotationally symmetric wide-angle lenses and devices thereof |
CN102629372A (zh) * | 2012-02-22 | 2012-08-08 | 北京工业大学 | 一种用于辅助车辆驾驶的360度全景鸟瞰图生成方法 |
CN102736634A (zh) * | 2012-07-10 | 2012-10-17 | 浙江捷尚视觉科技有限公司 | 一种用于汽车全景图的摄像头角度调节方法 |
CN106447735A (zh) * | 2016-10-14 | 2017-02-22 | 安徽协创物联网技术有限公司 | 一种全景相机几何标定处理方法 |
CN107240065A (zh) * | 2017-04-19 | 2017-10-10 | 中科院微电子研究所昆山分所 | 一种3d全景图像生成系统和方法 |
-
2018
- 2018-03-29 CN CN201810271061.XA patent/CN108876719B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009017332A1 (en) * | 2007-07-29 | 2009-02-05 | Nanophotonics Co., Ltd. | Methods of obtaining panoramic images using rotationally symmetric wide-angle lenses and devices thereof |
CN102629372A (zh) * | 2012-02-22 | 2012-08-08 | 北京工业大学 | 一种用于辅助车辆驾驶的360度全景鸟瞰图生成方法 |
CN102736634A (zh) * | 2012-07-10 | 2012-10-17 | 浙江捷尚视觉科技有限公司 | 一种用于汽车全景图的摄像头角度调节方法 |
CN106447735A (zh) * | 2016-10-14 | 2017-02-22 | 安徽协创物联网技术有限公司 | 一种全景相机几何标定处理方法 |
CN107240065A (zh) * | 2017-04-19 | 2017-10-10 | 中科院微电子研究所昆山分所 | 一种3d全景图像生成系统和方法 |
Non-Patent Citations (3)
Title |
---|
杨力 等: "一种嵌入式汽车鸟瞰全景图拼接算法", 《浙江大学学报(工学版)》 * |
王则浪: "基于DM8148的车载全景拼接算法及视频记录软件设计与实现", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 * |
魏利胜 等: "基于近似估计的汽车全景鱼眼镜头定位研究", 《控制工程》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109523489A (zh) * | 2018-11-28 | 2019-03-26 | 广州大学 | 一种生成俯视无畸变倒车影像的方法 |
CN109523489B (zh) * | 2018-11-28 | 2021-04-02 | 广州大学 | 一种生成俯视无畸变倒车影像的方法 |
CN110058263A (zh) * | 2019-04-16 | 2019-07-26 | 广州大学 | 一种车辆行驶过程中的物体定位方法 |
CN113840755A (zh) * | 2019-05-21 | 2021-12-24 | 康蒂-特米克微电子有限公司 | 生成车辆周围环境图像的方法及生成车辆周围环境图像的装置 |
CN114549666A (zh) * | 2022-04-26 | 2022-05-27 | 杭州蓝芯科技有限公司 | 一种基于agv的环视图像拼接标定方法 |
CN114549666B (zh) * | 2022-04-26 | 2022-09-06 | 杭州蓝芯科技有限公司 | 一种基于agv的环视图像拼接标定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108876719B (zh) | 2022-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108876719A (zh) | 基于虚拟相机模型的车辆全景图像拼接外部参数估计方法 | |
CN111223038B (zh) | 一种车载环视图像的自动拼接方法及显示装置 | |
CN104851076B (zh) | 用于商用车的全景环视泊车辅助系统及摄像头安装方法 | |
CN105825475B (zh) | 基于单摄像头的360度全景影像生成方法 | |
TWI555378B (zh) | 一種全景魚眼相機影像校正、合成與景深重建方法與其系統 | |
JP5739584B2 (ja) | 車両周辺視角化のための3次元映像合成装置およびその方法 | |
CN102999919B (zh) | 自运动估计系统和方法 | |
CN103826033B (zh) | 图像处理方法、图像处理设备、图像拾取设备和存储介质 | |
US10176595B2 (en) | Image processing apparatus having automatic compensation function for image obtained from camera, and method thereof | |
Zhang et al. | A surround view camera solution for embedded systems | |
CN103177439B (zh) | 一种基于黑白格角点匹配的自动标定方法 | |
CN103871071B (zh) | 一种用于全景泊车系统的摄像头外参标定方法 | |
CN106780590A (zh) | 一种深度图的获取方法及系统 | |
CN107424120A (zh) | 一种全景环视系统中的图像拼接方法 | |
EP1418766A2 (en) | Method and apparatus for processing images | |
CN103065318B (zh) | 多摄像机全景系统的曲面投影方法及装置 | |
CN111461963B (zh) | 一种鱼眼图像拼接方法及装置 | |
CN108592884B (zh) | 一种通用线阵卫星核线影像生成方法 | |
CN108198133A (zh) | 一种车辆全景图像快速拼接方法 | |
KR101705558B1 (ko) | Avm 시스템의 공차 보정 장치 및 방법 | |
CN110099268B (zh) | 色彩自然匹配与显示区自然融合的盲区透视化显示方法 | |
CN111047709A (zh) | 一种双目视觉裸眼3d图像生成方法 | |
CN109859279A (zh) | 一种基于线特征的车载360环视系统的流水线下线标定方法 | |
CN103106641A (zh) | 应用于全景成像系统的投影变换方法以及装置 | |
CN110428361A (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 |