发明名称 一种基于任意三面体的雷达-相机系统外部参数标定方法
摘要 本发明公开一种基于任意三面体的雷达-相机系统外部参数标定方法。利用自然环境中的三面体场景,仅需两帧数据即求解系统的外部参数。首先利用三面体规定世界坐标系,在雷达系统中对观察到的三面体进行平面拟合,得到各个平面的参数,并求得世界坐标系和雷达坐标系的转换关系以及雷达坐标系下两帧数据间的相对运动;在相机系统中,利用前后两帧提取的匹配特征点求解本质矩阵,然后分解出在相机坐标系下的相对运动,并利用雷达坐标系下的参数对相机坐标系下的平面参数进行求解,最后求解雷达-相机的外部参数,并利用两个坐标系下对应平面上点的共面性进行最终优化。本发明所要求的场景比较简单,兼具抗干扰性强、实验设备简单、灵活性较强的特点。
申请公布号 CN103049912A 申请公布日期 2013.04.17
申请号 CN201210563695.5 申请日期 2012.12.21
申请人 浙江大学 发明人 龚小谨;林颖;刘济林
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 林怀禹
主权项 1.一种基于任意三面体的雷达-相机系统外部参数标定方法,其特征在于该方法包括以下步骤: (1)利用雷达-相机系统,对存在任意三面体的标定场景采集两帧数据,并根据三面体定义世界坐标系O<sub>w</sub>-X<sub>w</sub>Y<sub>w</sub>Z<sub>w</sub>,采集的两帧数据间存在相对运动, 且每帧数据中有对应的雷达采集的三维点云和相机采集的彩色图像; (2)求解雷达系统参数:在三维点云中标注三面体三个面上的点,表示为: <img file="FDA00002629307900011.GIF" wi="687" he="66" />其中下标l代表这些点在雷达坐标系O<sub>l</sub>-X<sub>l</sub>Y<sub>l</sub>Z<sub>l</sub>下,i代表第i帧数据,j代表第j个面,k代表第k个点,每个面上点的个数为<img file="FDA00002629307900012.GIF" wi="71" he="65" />同一平面上的点用参数<img file="FDA00002629307900013.GIF" wi="173" he="64" />来表示:<img file="FDA00002629307900014.GIF" wi="313" he="66" />其中<img file="FDA00002629307900015.GIF" wi="63" he="65" />为每个面的法向量,<img file="FDA00002629307900016.GIF" wi="50" he="65" />为到坐标系原点的距离;然后根据三面体在雷达两帧数据中点的共面性,对两帧共6个面的<img file="FDA00002629307900017.GIF" wi="171" he="64" />进行非线性优化;通过<img file="FDA00002629307900018.GIF" wi="171" he="64" />得到两帧数据中世界坐标系到雷达坐标系的旋转矩阵和平移向量<img file="FDA00002629307900019.GIF" wi="174" he="56" />i=1,2,并由此得到雷达-相机系统在雷达坐标系下两帧之间的运动参数,即第一帧相对于第二帧的旋转矩阵和平移向量<img file="FDA000026293079000110.GIF" wi="182" he="70" />(3)求解相机系统参数:利用SIFT算法在相机采集的前后两帧图像中三面体三个面的对应区域提取前后帧匹配的特征点,表示为: <img file="FDA000026293079000111.GIF" wi="691" he="66" />其中下标c代表这些点在相机坐标系O<sub>c</sub>-X<sub>c</sub>Y<sub>c</sub>Z<sub>c</sub>下,i代表第i帧数据,j代表第j个面,k代表第k个特征点,前后两帧对应面的匹配特征点数目相等为<img file="FDA000026293079000112.GIF" wi="71" he="57" />根据P<sub>c</sub>计算出两帧运动的本质矩阵E,并进一步得到雷达-相机系统在相机坐标系下两帧之间的运动参数<img file="FDA000026293079000113.GIF" wi="196" he="56" />其中<img file="FDA000026293079000114.GIF" wi="77" he="56" />为相机坐标系下第一帧相对于第二帧的旋转矩阵,<img file="FDA000026293079000115.GIF" wi="61" he="53" />为平移向量<img file="FDA000026293079000116.GIF" wi="69" he="56" />归一化的结果:<img file="FDA000026293079000117.GIF" wi="519" he="80" />利用<img file="FDA00002629307900021.GIF" wi="377" he="64" />初始化相机坐标系下的参数<img file="FDA00002629307900022.GIF" wi="374" he="64" />,此处<img file="FDA00002629307900023.GIF" wi="178" he="65" />为三面体各个平面在相机坐标系下的参数;根据提取的特征点的对应投影关系,对<img file="FDA00002629307900024.GIF" wi="375" he="64" />进行优化;通过<img file="FDA00002629307900025.GIF" wi="178" he="64" />得到两帧数据中世界坐标系到相机坐标系的转换关系,即旋转矩阵和平移向量<img file="FDA00002629307900026.GIF" wi="183" he="57" />i=1,2;(4)求解外部参数:利用R<sub>wc</sub>,T<sub>wc</sub>,R<sub>wl</sub>,T<sub>wl</sub>得到雷达坐标系到相机坐标系的转换关系,即旋转矩阵R<sub>lc</sub>和平移向量T<sub>lc</sub>;利用雷达坐标系和相机坐标系下三面体上点的共面性,对R<sub>lc</sub>,T<sub>lc</sub>进行非线性优化,得到优化后的转换关系R<sub>cl</sub>,T<sub>lc</sub>。 
地址 310027 浙江省杭州市西湖区浙大路38号
您可能感兴趣的专利