发明名称 摄像头与惯性传感器组合定位定姿系统的标定方法
摘要 本发明为摄像头与惯性传感器组合定位定姿系统的标定方法,包括:标定摄像头的内参数矩阵;对于一个尺寸已知的标定物从不同角度拍摄多幅图像,记录在拍摄每幅图像时惯性传感器输出的横滚角和俯仰角;定义世界坐标系、摄像机坐标系、惯性传感器坐标系和地磁坐标系;基于每幅图像中标定物的图像信息和空间信息,计算该时刻世界坐标系到摄像机坐标系的旋转矩阵;将摄取图像两两组合,对于每个组合均建立一个关于惯性传感器坐标系到摄像机坐标系旋转矩阵的方程组并进行求解,计算出惯性传感器坐标系到摄像机坐标系的旋转矩阵;对于每幅图像建立一个关于地磁坐标系到世界坐标系旋转矩阵的方程组并进行求解,计算出地磁坐标系到世界坐标系的旋转矩阵。
申请公布号 CN102162738B 申请公布日期 2012.11.21
申请号 CN201010591519.3 申请日期 2010.12.08
申请人 中国科学院自动化研究所 发明人 董秋雷;顾照鹏;胡占义
分类号 G01C25/00(2006.01)I;G06T7/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 中科专利商标代理有限责任公司 11021 代理人 梁爱荣
主权项 摄像头与惯性传感器组合定位定姿系统的标定方法,包括以下步骤:步骤S1:标定摄像头的内参数矩阵;步骤S2:对于一个尺寸已知的标定物从不同角度拍摄多幅图像,同时记录在拍摄每幅图像时惯性传感器输出的横滚角和俯仰角信息;步骤S3:定义世界坐标系、摄像机坐标系、惯性传感器坐标系和地磁坐标系;步骤S4:基于每幅图像中标定物的图像信息和空间信息,计算对应于该幅图像的世界坐标系到摄像机坐标系的旋转矩阵;步骤S5:将所拍摄图像两两组合生成多个图像对,对于每个图像对建立一个关于惯性传感器坐标系到摄像机坐标系旋转矩阵Rm‑c的方程组: <mrow> <msub> <mi>R</mi> <mrow> <mi>w</mi> <mo>-</mo> <mi>c</mi> <mo>,</mo> <mi>ij</mi> </mrow> </msub> <mo>&times;</mo> <msub> <mi>R</mi> <mrow> <mi>m</mi> <mo>-</mo> <mi>c</mi> </mrow> </msub> <mo>&times;</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mo>-</mo> <msub> <mrow> <mi>sin</mi> <mi>&beta;</mi> </mrow> <mi>j</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mi>cos</mi> <msub> <mi>&beta;</mi> <mi>j</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mi>cos</mi> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mi>cos</mi> <msub> <mi>&beta;</mi> <mi>j</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msub> <mi>R</mi> <mrow> <mi>m</mi> <mo>-</mo> <mi>c</mi> </mrow> </msub> <mo>&times;</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mo>-</mo> <mi>sin</mi> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mi>sin</mi> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> <mi>cos</mi> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mi>cos</mi> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> <mi>cos</mi> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>其中,Rw‑c,ij表示摄像机坐标系在i、j两时刻之间的旋转矩阵,αi,αj分别表示惯性传感器在i、j时刻的横滚角,βi,βj分别表示惯性传感器在i、j时刻的俯仰角,则对于多个图像对,利用惯性传感器输出的横滚角和俯仰角,可建立关于惯性传感器坐标系到摄像机坐标系旋转矩阵Rm‑c的多个方程组;步骤S6:将关于惯性传感器坐标系到摄像机坐标系旋转矩阵Rm‑c的多个方程组转化为一个含约束的非线性优化问题,求解所述含约束的非线性优化问题则计算出惯性传感器坐标系到摄像机坐标系的旋转矩阵Rm‑c;步骤S7:对于每幅图像建立一个关于地磁坐标系到世界坐标系旋转矩阵Re‑w的方程组: <mrow> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>cos</mi> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> </mtd> <mtd> <mi>sin</mi> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> </mtd> <mtd> <mi>cos</mi> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> <mi>sin</mi> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>cos</mi> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> </mtd> <mtd> <mo>-</mo> <mi>sin</mi> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mi>sin</mi> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> </mtd> <mtd> <mi>sin</mi> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> <mi>cos</mi> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> </mtd> <mtd> <mi>cos</mi> <msub> <mi>&alpha;</mi> <mi>i</mi> </msub> <mi>cos</mi> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>&times;</mo> <msubsup> <mi>R</mi> <mrow> <mi>m</mi> <mo>-</mo> <mi>c</mi> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>&times;</mo> <msub> <mi>R</mi> <mrow> <mi>w</mi> <mo>-</mo> <mi>c</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>&times;</mo> <msub> <mi>R</mi> <mrow> <mi>e</mi> <mo>-</mo> <mi>w</mi> </mrow> </msub> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <mi>cos</mi> <msub> <mi>&gamma;</mi> <mi>i</mi> </msub> </mtd> <mtd> <mi>sin</mi> <msub> <mi>&gamma;</mi> <mi>i</mi> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mi>sin</mi> <msub> <mi>&gamma;</mi> <mi>i</mi> </msub> </mtd> <mtd> <mi>cos</mi> <msub> <mi>&gamma;</mi> <mi>i</mi> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>其中,Rw‑c,i表示在i时刻世界坐标系到摄像机坐标系的旋转矩阵,Rm‑c表示惯性传感器坐标系到摄像机坐标系的旋转矩阵,αi表示惯性传感器在i时刻的横滚角,βi表示惯性传感器在i时刻的俯仰角,γi表示惯性传感器在i时刻的航向角;利用等式(2)两边对应项相等的性质以及旋转矩阵的正交性约去有关惯性传感器航向角γi的计算项,可得如下线性方程组:BiVe‑w=bi    (3)其中Bi,bi分别为等式(2)经过化简后的系数矩阵与常数项,Ve‑w是地磁坐标系到世界坐标系旋转矩阵Re‑w的向量化表示形式,则对于多幅图像,利用惯性传感器输出的横滚角和俯仰角,可建立关于地磁坐标系到世界坐标系旋转矩阵Re‑w的多个方程组,进而建立关于地磁坐标系到世界坐标系旋转矩阵Re‑w的向量化表示形式Ve‑w的多个线性方程组;步骤S8:将关于地磁坐标系到世界坐标系旋转矩阵Re‑w的向量化表示形式Ve‑w的多个线性方程组转化为一个含约束的非线性优化问题,求解所述含约束的非线性优化问题则计算出地磁坐标系到世界坐标系的旋转矩阵。
地址 100190 北京市海淀区中关村东路95号