FAST馈源舱的静态定标和动态跟踪

作为国际新一代大射电望远镜(LT)阵计划的第一步,拟在我国先行实施一项FAST(Five hundred meters Aperture Spherical Telescope)项目工程｡目前国际上正在更新的 Arecibo 系统难以满足 LT 基本单元的要求:低造价､大天空覆盖､宽带以及偏振观测｡全球最大的射电望远镜是位于美国波多里格的 Arecibo305m 口径天线,但它具有天空覆盖小(天顶扫描角仅20°)的严重缺陷,以及造价太高､跟踪精度低的不足｡FAST项目计划利用我国某地独一无二的喀斯特(Karst)洼地,铺设主天线的球面望远镜,建造口径为500m的射电望远镜｡这种射电望远镜取消了主反射面的运动,改用馈源移动跟踪目标,基本不受重力形变的影响,可把主反射面建造得很大｡

对射电望远镜的馈源舱实施闭环控制的前提条件是已知馈源舱的位置及姿态,故需对馈源舱进行动态跟踪,以取得相关数据｡本文根据计算机视觉和CCD图像分析测量原理,详细介绍了对实验模型中的馈源舱进行静态定标与动态跟踪测量的原理及方法｡

1 CCD测量系统

结合课题情况,可考虑使用的测量方法有以下四种:(1)GPS定位系统:测量范围大,差分处理后的测量精度较高,不足是测量时间较长｡(2)无线电定位:受环境影响小,测量范围大,可在较恶劣的环境中工作｡但测量成本较高,且无线电信号会影响射电天文望远镜对宇宙信号的接收｡(3)激光全站仪:测量范围大､测量精度高､采样周期高｡但系统的采样间隔具有不稳定性､时延性与较低的动态精度(3mm),这给控制带来较大难度｡另外受环境影响较大,在降雨和大风扬尘等较恶劣的环境下,测量精度会受影响,且价格很高(一台 Leica 大约价值18万元)｡(4)CCD三维测量系统:成本低,测量范围较小时测量精度较高｡但由于测量数据量大,动态跟踪测量的时间较长｡另外环境因素对测量精度的影响也较大,夜间工作有一定限制｡

根据实际情况以及顺利实现测量的目的,测量系统应具有快速测量､自动跟踪和成本低的特点｡而CCD三维测量系统能够满足要求｡CCD器件是一种固体化器件,体积小､可靠性高､寿命长;图像畸变小,尺寸重现性好;具有较高的定位精度和测量精度;输出信号易于数字化处理,易与计算机连接组成实时自动测量控制系统,便于扩大应用功能和使用范围｡

CCD图像分析测量系统的原理框图如图1所示｡



2 测量原理与方法

CCD三维测量系统的组成元件主要有CCD摄像机(MTV-1881EX､795x596)､图像卡(大恒 CG210)､计算机和视频线等｡整个系统的工作原理是:测量对象在CCD摄像机的测量范围内沿任意方向运动,CCD摄像机从三个不同的角度对测量对象的特征点进行摄像,生成被测对象的视频信号; 通过图像卡转换成数字信号并输送给控制计算机;计算机调用执行文件,根据一定算法计算被测目标的世界坐标,由此确定对象的位置与姿态｡

50m实验模型中的馈源舱由六根索向上拉,并分别通过六座钢筋水泥塔与地面的卷扬机相连｡其中三根索均布接在舱体顶部,另三根索均布接在舱体底圆上｡取下拉索与舱体的绞合点(索耳中心)为特征点,这样共有三个特征点a､b､c,三部CCD摄像机按照π/3间隔放置,结构分布如见图2｡





令Pi是待测特征点,则其世界坐标与其在摄像机中的投影坐标的关系式为:

(x y z 1)T=RH(穴x′ y′ z′ 1)T (1)

其中R､H分别是空间旋转和平移变换矩阵｡

每个特征点分别在左右两个摄像机中投影,投影坐标满足如下关系式:



式中:分别为左右摄像机镜头沿径向的畸变程度,(xl,yl)､(xr,yr)分别是左右摄像机镜头的光心坐标,Al1…,Al11 ､Ar1…,Ar11 分别是测量点在左右两摄像机的投影坐标变换参数｡

采用静态定标的方法确定式(2a)､(2b)中的未知参数｡测量系统静态定标的原理是:在地面上合适位置(坐标已知)安置两架电子经纬仪,并在馈源舱的工作空间区域选择测量点,利用经纬仪测量若干位置特征点的俯仰角和方位角,通过坐标几何变换确定这些点的世界坐标(x,y,z),并认为是实际坐标｡同时记下它们在摄像机中相应的图像坐标(u,v),这里需注意对应关系｡把世界坐标与相应投影坐标代人式(2a)､(2b)中,用最小二乘法求解超静定方程组,确定未知参数｡

确定空间坐标变换矩阵后需进行静态检测,即通过经纬仪测量一些静态特征点的坐标,与CCD静态测量结果比较,计算静态定标的rms误差｡若不满足精度要求,则重新定标直至满足｡

静态定标后即可对特征点进行动态跟踪｡其跟踪原理是基于面积的边沿提取跟踪算法｡每部摄像机读取并确定其中两个特征点在CCD靶面的投影坐标,得到三特征点在摄像机中的六个投影坐标后,根据投影坐标的位置关系与特征点附近区域在靶面投影区域的相关性,确定每个特征点在其对应的两摄像机靶面的投影坐标,再采用双目定位法得到它们的世界坐标｡例如摄像机3和1可同时对目标点a进行测量跟踪,摄像机1在t时刻采集到图像的左特征点 at 附近方形区域作为模板 TRt ,摄像机3在 t 时刻采集到图像的右特征点 at 附近方形区域作为模板 TLt ,实现两幅图的配准｡然后分别以 TLt 和 TRt 为模板,在 t+1 时刻两摄像机3､1采集的两幅图像中搜索有相似灰度值分布的 TLt+1 和 TRt+1 (如图3)｡判断此两模板是否满足最大相关性,若满足,则认为两模板的中心点就是摄像机动态跟踪目标点的投影;否则继续搜索特征区域,直至满足｡



两模板最相似需满足的最大相关性条件是:当最小时,右模板中的特征点匹配｡

动态跟踪测量给出馈源舱三只索耳 C1､C2 和 C3 的直角坐标(x1,y1,z1),(x2,y2,z2)和(x3,y3,z3)｡采用下述方法反算馈源舱的位置与姿态｡因索耳均布在舱体底圆上(如图2),则底圆中心O1的直角坐标为:xO1=(x1+x2+x3)/3, yO1=(y1+y2+y3)/3,zO1=(z1+z2+z3)/3 (4)

设舱体底圆中心O1与顶端O2的连线为舱体对称轴Z1的单位矢量,则:



其中分别表示索耳C1指向C2和C3的矢量｡所以馈源舱的方位角α与俯仰角γ为:

α=arctank1y/k1x (6a) γ=arccosk1z (6b)