专利名称：一种摄像机外部参数自定标方法及装置的制作方法
技术领域：
本发明涉及摄像机定标技术领域，尤其是涉及一种广角摄像机外部参数自定标方法及装置。
背景技术：
空间物体表面某点的三维几何位置与其在图像中对应点间的相互关系由摄像机成像的几何模型决定，其参数就是摄像机参数。该参数包括摄像机内部几何和光学特性 (内部参数)，以及摄像机相对世界坐标系的三维位置和方向(外部参数)。由于摄像机在实际使用中经常发生位置变化，因此摄像机的外部参数标定越来越被重视。摄像机的外部参数定标的应用已涉及生活中的各个领域，如机器人手眼定标、移动摄像机定标、侦查摄像机定标等等。因此，研究开发一种定标时间短、精度高的摄像机外部参数定标方法，将能够有效地提高摄像机在生产和生活中的应用，为生产和生活带来便利。目前已有的摄像机外部参数定标算法主要分为传统算法和智能算法，其中传统算法往往利用特征点的对应关系，能快速地计算出外部参数值，且由于计算模型不含随机数，其计算结果往往是唯一的，因此还具有稳定的优点。但是，传统算法在提取特征点的时候会有误差，导致经过计算模型之后误差仍延续存在，从而导致结果精度并不高。智能算法利用算法的智能性通过迭代的方式求得结果，使得结果能够逐步趋于精确。然而由于算法中加入随机数，使得结果并不稳定，同时算法的运行效率与初始化解密切相关。可见，现有的摄像机外部参数标定方法无法兼顾运算时间短和结果精度高。
发明内容
本发明为了解决现有技术摄像机外部参数标定方法无法兼顾运算时间短和结果精度高的技术问题，提供了一种摄像机外部参数自标定方法及装置。为解决上述技术问题，本发明采用的技术方案为设计一种摄像机外部参数自定标方法，包括下列步骤A、利用摄像机获取模板的灰度图像，并存储；B、对获取的灰度图像进行去噪处理；C、检测去噪处理后的灰度图像的角点；D、根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数；E、利用粒子群算法优化所述初始摄像机外部参数，得到摄像机外部参数。所述对获取的灰度图像进行去噪处理包括先对获取的灰度图像进行二值化处理，再对二值化处理后的图像进行多次先膨胀处理后腐蚀处理。所述检测去噪处理后的灰度图像的角点采用Harris算子进行角点检测。所述根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数包括
根据非线性摄像机模型，建立世界坐标与图像坐标的对应关系；根据角点信息和非线性摄像机模型，求得初始摄像机外部参数。所述利用粒子群算法优化所述初始摄像机外部参数包括根据初始摄像机外部参数，初始化粒子位置；定义粒子的适应值函数；根据粒子的适应值进行优化。本发明还提供了一种摄像机外部参数自定标装置，包括灰度图像获取模块，其利用摄像机获取模板的灰度图像，并存储；去噪模块，其对获取的灰度图像进行去噪处理；角点检测模块，其检测去噪处理后的灰度图像的角点；初始摄像机外部参数获取模块，其根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数；优化模块，其利用粒子群算法优化所述初始摄像机外部参数，得到摄像机外部参数。所述去噪模块包括二值化处理模块，其对获取的灰度图像进行二值化处理；膨胀腐蚀处理模块，其对二值化处理后的图像进行多次先膨胀处理后腐蚀处理。所述检测去噪处理后的灰度图像的角点采用Harris算子进行角点检测。所述根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数包括根据非线性摄像机模型，建立世界坐标与图像坐标的对应关系；根据角点信息和非线性摄像机模型，求得初始摄像机外部参数。所述利用粒子群算法优化所述初始摄像机外部参数包括根据初始摄像机外部参数，初始化粒子位置；定义粒子的适应值函数；根据粒子的适应值进行优化。。本发明通过先对获取的灰度图像进行去噪处理，再检测去噪处理后的灰度图像的角点，然后再根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数，最后，利用粒子群算法优化所述初始摄像机外部参数，得到摄像机外部参数。从而利用将传统算法与智能算法结合在一起的混合算法，不仅继承了传统算法高效稳定的优点，同利用传统算法得到优良的初始化解，来帮助智能算法逐步求精，运算时间短、结果精度高。
下面结合实施例和附图对本发明进行详细说明，其中图1是本发明摄像机外部参数自定标方法的流程图；图2是本发明摄像机外部参数自定标装置的原理图。
具体实施例方式请参见图1。本发明摄像机外部参数自定标方法包括下列步骤一、利用摄像机获取模板的灰度图像，并存储。
5
二、对获取的灰度图像进行去噪处理。对获取的灰度图像进行去噪处理，以降低图像质量对后面算法的影响。在该步骤中，先对获取的灰度图像进行二值化处理，再对二值化处理后的图像进行多次先膨胀处理后腐蚀处理。在本具体实施例中，优选进行3次先膨胀后腐蚀处理。2. 1对获取的灰度图像进行二值化处理得到二值化处理后的图像定义获取的灰度图像为II，二值化处理后的图像为12，Il(x，y)、I2(x，y)分别表示II、12在坐标(X，y)处的灰度值。记S为前景与背景的分割阈值，前景像素点占图像比例为wO，平均灰度为m0 ；背景像素点占图像比例为wl，平均灰度为ml。图像的总平均灰度为m = wO*mO+wl*ml。从最小灰度值到最大灰度值遍历s，当s = sO时使得函数F = w0*(m0-m)2+wl*(ml-m)2值最大。再对灰度图像Il进行操作，当Il (x，y) < sO时，I2(x，y) = 0;ill(x，y) >= sO时，I2(x，y) =1。由此得到二值化处理后的图像12。2. 2对二值化处理后的图像12进行3次先膨胀后腐蚀处理得到去噪处理后的图像 13。膨胀处理为用3*3的结构元素，扫描图像的每一个像素，用结构元素与其覆盖的二值图像做“与”操作，如果都为0，结果图像的该像素为0，否则为1。腐蚀处理为用3*3的结构元素，扫描图像的每一个像素，用结构元素与其覆盖的二值图像做“与”操作，如果都为1，结果图像的该像素为1，否则为0。三、检测去噪处理后的灰度图像的角点在该步骤中，对去噪处理后的图像13进行角点检测处理。采用Harris算子P进行角点检测。当某点的Harris算子P值大于设定阀值TO时，该点为灰度图像的角点；否则， 该点不是灰度图像的角点。四、根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数在该步骤中，首先根据非线性摄像机模型，建立世界坐标与图像坐标的对应关系； 然后再根据角点信息和非线性摄像机模型，求得初始摄像机外部参数。4. 1根据非线性摄像机模型，建立世界坐标与图像坐标的对应关系建立摄像机坐标与图像坐标的对应关系，其计算公式为
权利要求
1.一种摄像机外部参数自定标方法，其特征在于包括下列步骤A、利用摄像机获取模板的灰度图像，并存储；B、对获取的灰度图像进行去噪处理；C、检测去噪处理后的灰度图像的角点；D、根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数；E、利用粒子群算法优化所述初始摄像机外部参数，得到摄像机外部参数。
2.根据权利要求1所述的摄像机外部参数自定标方法，其特征在于所述对获取的灰度图像进行去噪处理包括先对获取的灰度图像进行二值化处理，再对二值化处理后的图像进行多次先膨胀处理后腐蚀处理。
3.根据权利要求1所述的摄像机外部参数自定标方法，其特征在于所述检测去噪处理后的灰度图像的角点采用Harris算子进行角点检测。
4.根据权利要求1所述的摄像机外部参数自定标方法，其特征在于所述根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数包括根据非线性摄像机模型，建立世界坐标与图像坐标的对应关系； 根据角点信息和非线性摄像机模型，求得初始摄像机外部参数。
5.根据权利要求1所述的摄像机外部参数自定标方法，其特征在于所述利用粒子群算法优化所述初始摄像机外部参数包括根据初始摄像机外部参数，初始化粒子位置； 定义粒子的适应值函数； 根据粒子的适应值进行优化。
6.一种摄像机外部参数自定标装置，其特征在于包括灰度图像获取模块，其利用摄像机获取模板的灰度图像，并存储； 去噪模块，其对获取的灰度图像进行去噪处理； 角点检测模块，其检测去噪处理后的灰度图像的角点；初始摄像机外部参数获取模块，其根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数；优化模块，其利用粒子群算法优化所述初始摄像机外部参数，得到摄像机外部参数。
7.根据权利要求6所述的摄像机外部参数自定标装置，其特征在于所述去噪模块包括二值化处理模块，其对获取的灰度图像进行二值化处理；膨胀腐蚀处理模块，其对二值化处理后的图像进行多次先膨胀处理后腐蚀处理。
8.根据权利要求6所述的摄像机外部参数自定标装置，其特征在于所述检测去噪处理后的灰度图像的角点采用Harris算子进行角点检测。
9.根据权利要求6所述的摄像机外部参数自定标装置，其特征在于所述根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数包括根据非线性摄像机模型，建立世界坐标与图像坐标的对应关系； 根据角点信息和非线性摄像机模型，求得初始摄像机外部参数。
10.根据权利要求6所述的摄像机外部参数自定标装置，其特征在于所述利用粒子群算法优化所述初始摄像机外部参数包括根据初始摄像机外部参数，初始化粒子位置； 定义粒子的适应值函数； 根据粒子的适应值进行优化。
全文摘要
本发明公开了一种摄像机外部参数自定标方法，旨在提供一种运算时间短、结果精度高的摄像机外部参数自定标方法，其包括下列步骤A、利用摄像机获取模板的灰度图像，并存储；B、对获取的灰度图像进行去噪处理；C、检测去噪处理后的灰度图像的角点；D、根据非线性摄像机模型和所述角点信息求得初始摄像机外部参数；E、利用粒子群算法优化所述初始摄像机外部参数，得到摄像机外部参数。本发明还公开了一种摄像机外部参数自定标装置。本发明可用于摄像机外部参数的定标。
文档编号G06T7/00GK102354399SQ20111028775
公开日2012年2月15日 申请日期2011年9月26日 优先权日2011年9月26日
发明者张元标, 毕小明, 潘兴伟 申请人:深圳普捷利科技有限公司