[发明专利]交通监视环境下的路况摄像机标定方法无效

专利信息
申请号: 200710022810.7 申请日: 2007-05-22
公开(公告)号: CN101118648A 公开(公告)日: 2008-02-06
发明(设计)人: 陈启美;李勃;郭凡;董蓉 申请(专利权)人: 南京大学
主分类号: G06T7/00 分类号: G06T7/00;G08G1/01;H04N5/225
代理公司: 南京天翼专利代理有限责任公司 代理人: 汤志武;王鹏翔
地址: 210093*** 国省代码: 江苏;32
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 交通 监视 环境 路况 摄像机 标定 方法
【权利要求书】:

1.交通监视环境下的路况摄像机标定方法,其特征是包含如下标定步骤:

(1)视觉模型描述和相关坐标系建立:针对监控系统性能要求,沿用经典的Tsai透射投影模型,并针对路况成像特点,对其进行相应修正,提出新的视觉模型,建立三种坐标系:

其中地面坐标系Xw-Yw-Zw和摄像机坐标系Xc-Yc-Zc用来表征三维空间;图像平面坐标系Xf-Yf用来表征成像平面。建立世界坐标系,其原点为摄像机光轴与地面交点。Yw轴正向沿路面方向指向前方,Xw轴正向水平指向右方,Zw轴正向垂直于地面,方向向上。建立摄像机坐标系,原点为摄像机光心位置,Zc轴为摄像机光轴方向,Xc-Yc平面平行于像平面;

Xc=(cos(p)cos(s)+sin(t)sin(p)sin(s))Xw

  +(sin(p)cos(s)-sin(t)cos(p)sin(s))Yw

Yc=(-cos(p)sin(s)+sin(t)sin(p)cos(s))Xw (1)

  +(-sin(p)sin(s)+sin(t)cos(p)cos(s))Yw

Zc=-cos(t)sin(p)Xw+cos(t)cos(p)Yw+l

<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mn>1</mn><mrow><mi>M</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>[</mo><msub><mi>X</mi><mi>f</mi></msub><mo>-</mo><msub><mi>C</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>Z</mi><mo>)</mo></mrow><mo>]</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>kr</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mfrac><msub><mi>X</mi><mi>C</mi></msub><msub><mi>Z</mi><mi>C</mi></msub></mfrac></mtd></mtr><mtr><mtd><mfrac><mn>1</mn><mrow><mi>M</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>[</mo><msub><mi>Y</mi><mi>f</mi></msub><mo>-</mo><msub><mi>C</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>Z</mi><mo>)</mo></mrow><mo>]</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>kr</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>f</mi><mfrac><msub><mi>Y</mi><mi>C</mi></msub><msub><mi>Z</mi><mi>C</mi></msub></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>

建立此视觉模型下,世界坐标系与像平面坐标系下映射关系,如式(3)所示

<mrow><msub><mi>X</mi><mi>f</mi></msub><mo>=</mo><mfrac><mrow><mi>f</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mrow><mo>(</mo><mi>cos</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><mi>sin</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>)</mo></mrow><msub><mi>X</mi><mi>w</mi></msub></mtd></mtr><mtr><mtd><mo>+</mo><mrow><mo>(</mo><mi>sin</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>-</mo><mi>sin</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>)</mo></mrow><msub><mi>Y</mi><mi>w</mi></msub></mtd></mtr></mtable></mfenced></mrow><mrow><mo>-</mo><mi>cos</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><msub><mi>X</mi><mi>w</mi></msub><mo>+</mo><mi>cos</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><msub><mi>Y</mi><mi>w</mi></msub><mo>+</mo><mi>l</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>

<mrow><msub><mi>Y</mi><mi>f</mi></msub><mo>=</mo><mfrac><mrow><mi>f</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mrow><mo>(</mo><mo>-</mo><mi>cos</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><mi>sin</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>)</mo></mrow><msub><mi>X</mi><mi>w</mi></msub></mtd></mtr><mtr><mtd><mo>+</mo><mrow><mo>(</mo><mo>-</mo><mi>sin</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><mi>sin</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>)</mo></mrow><msub><mi>Y</mi><mi>w</mi></msub></mtd></mtr></mtable></mfenced></mrow><mrow><mo>-</mo><mi>cos</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><msub><mi>X</mi><mi>w</mi></msub><mo>+</mo><mi>cos</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><msub><mi>Y</mi><mi>w</mi></msub><mo>+</mo><mi>l</mi></mrow></mfrac></mrow>

由式(3)出发,以Zw为已知参数,同样可建立由像平面坐标到世界坐标系下的逆映射关系,如式(4)所示:

<mrow><msub><mi>X</mi><mi>w</mi></msub><mo>=</mo><mfrac><mfenced open='[' close=']'><mtable><mtr><mtd><mi>sin</mi><mi>p</mi><mrow><mo>(</mo><mi>l</mi><mo>+</mo><msub><mi>Z</mi><mi>w</mi></msub><mi>sin</mi><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>f</mi></msub><mi>sin</mi><mi>s</mi><mo>+</mo><msub><mi>y</mi><mi>f</mi></msub><mi>cos</mi><mi>s</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>+</mo><mi>cos</mi><mi>p</mi><mrow><mo>(</mo><mi>l</mi><mi>sin</mi><mi>t</mi><mo>+</mo><msub><mi>Z</mi><mi>w</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>f</mi></msub><mi>cos</mi><mi>s</mi><mo>-</mo><msub><mi>y</mi><mi>f</mi></msub><mi>sin</mi><mi>s</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Z</mi><mi>w</mi></msub><mi>f</mi><mi>cos</mi><mi>t</mi><mi>sin</mi><mi>p</mi></mtd></mtr></mtable></mfenced><mrow><msub><mi>x</mi><mi>f</mi></msub><mi>cos</mi><mi>t</mi><mi>sin</mi><mi>s</mi><mo>+</mo><msub><mi>y</mi><mi>f</mi></msub><mi>cos</mi><mi>t</mi><mi>cos</mi><mi>s</mi><mo>+</mo><mi>f</mi><mi>sin</mi><mi>t</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>

<mrow><msub><mi>Y</mi><mi>w</mi></msub><mo>=</mo><mfrac><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><mi>cos</mi><mi>p</mi><mrow><mo>(</mo><mi>l</mi><mo>+</mo><msub><mi>Z</mi><mi>w</mi></msub><mi>sin</mi><mi>t</mi><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>f</mi></msub><mi>sin</mi><mi>s</mi><mo>+</mo><msub><mi>y</mi><mi>f</mi></msub><mi>cos</mi><mi>s</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>+</mo><mi>sin</mi><mi>p</mi><mrow><mo>(</mo><mi>l</mi><mi>sin</mi><mi>t</mi><mo>+</mo><msub><mi>Z</mi><mi>w</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>f</mi></msub><mi>cos</mi><mi>s</mi><mo>-</mo><msub><mi>y</mi><mi>f</mi></msub><mi>sin</mi><mi>s</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>Z</mi><mi>w</mi></msub><mi>f</mi><mi>cos</mi><mi>t</mi><mi>sin</mi><mi>p</mi></mtd></mtr></mtable></mfenced><mrow><msub><mi>x</mi><mi>f</mi></msub><mi>cos</mi><mi>t</mi><mi>sin</mi><mi>s</mi><mo>+</mo><msub><mi>y</mi><mi>f</mi></msub><mi>cos</mi><mi>t</mi><mi>cos</mi><mi>s</mi><mo>+</mo><mi>f</mi><mi>sin</mi><mi>t</mi></mrow></mfrac></mrow>

式中H为摄像机的垂直安置高度,t,p,s分别为摄像机俯仰角、偏角和旋角,f为广角有效焦距,M(z)为放缩系数,Cx(z)、Cy(z)为图像主点坐标,k为一阶径向畸变。

(2)标定摄像机主点及放缩系数:以监视图像光流作为标定基元;通过摄像机作放缩运动,应用参考帧预测图像及实时帧采样图像间光流场差值作为约束,使用最小二乘方法建立约束方程,通过powell方向族法辨识摄像机主点坐标及摄像机实际放大系数。

(3)标定目标选择与参数线性求解:在监视场景中选择分道线角点作为标定参照物,利用相邻四个分道线角点的像平面投影确定消失线斜率及相应摄像机旋角;以分道线平行关系及分道线间基础路宽作为约束,利用此四个角点的像平面坐标,线性求解摄像机有效焦距及其空间位置参数;

(4)求精摄像机内外部参数:以监视图像中所有角点和对应的世界坐标点,以上述算法所求出的视觉模型参数作为优化模型初值,采用Levenberg-Marquardt优化算法求精摄像机模型参数,完成摄像机标定。

2.根据权利要求1所述的交通监视环境下的路况摄像机标定方法,其特征是采用摄像机畸变补偿:利用已定标出完备的理想的无透视畸变摄像机成像模型,计算图像上所有角点的理想象素坐标,通过其与真实角点像平面投影的坐标差值,求解非线性视觉模型中的一阶径向畸变。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于南京大学,未经南京大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/200710022810.7/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top