-
弹丸落点,即相对于被击目标的炸点位置,是弹丸破片毁伤效能分析和弹丸改进研制以及打击点方位时速设置的重要分析依据。随着军事演习乃至军事活动范围增大,弹丸水面落点的测量越来越成为保证水面打击准确度的前提。对弹丸着靶坐标测量一般可利用光学和声学测量的方法进行。传统的光幕靶测量和声靶测量等都有各自的缺陷,例如高速摄像法采用多个高速摄像机交汇的方法进行测量,其成本比较高,声靶结构大,对水面特定目标环境测试不仅存在着布置困难,而且水面环境的不稳定造成误差大,最终仅获得火炮方向的两维坐标[1]。对于船载经纬仪交会测量,由于风浪、潮汐等洋流不规则的影响,使船舶平台的姿态和形体发生变化,船体的位置也不如陆基测量站稳定并且对设备自身精度指标的要求较高,其测量过程在动态环境下进行,设备自身的位置、姿态在不断地变化,布置成本昂贵并且精度不高[2,3],虽然有很多文献提出了对测量结果进行误差修正,但是在实际工程中的实现度并不是很理想。弹丸落点在水面的特殊情况造成无法直接定位测量。为解决水面弹着点测量困难的问题,孙朝江等学者提出了双 CCD 交汇测量的方法,利用两个线阵 CCD 探测器实现对舰炮武器水面落点的测量,实现技术难度较小,方法简单,但是存在需要解决CCD固定的问题,同时还要防止水浪的侵袭等[4]。类似的火箭助飞鱼雷海上落点是由矢量水听器的被动测量系统完成的,主要是利用布设的携带有矢量水听器的浮标,完成入水击水声的方位估计,但预设浮标布阵会出现海上走锚、通信不畅、性能下降,导致测量的可靠性和精度严重下降[5]。
本文采用结合每时刻相机采集的靶标船上任意三定点以及目标的图像,提取t时刻定点像面坐标和由空间坐标关系算出的对应点的姿态信息间的关系传递函数,推算出目标落点图像坐标到方位、俯仰角的映射关系,得出目标落点相对于两测量船的姿态信息,最终利用异面交会算法估算出水面落点的空间三维坐标。方法的关键是通过定点质心来确定其图像坐标,并建立两套参系统的关系传递函数,最后通过合理选择(对交会结果影响最小)测量船布站方案确定落点的空间坐标。实验数据显示,本文提出的方法能够在大地误差、质心误差存在以及交会误差最小的情况下,使目标水面落点位置的确定达到较为理想的精度。
该方法选取地心坐标系为标准坐标系,结合CCD动态采集的图像信息建立的关系函数模型,不但有效地避免了水面落点特殊环境存在的抖动误差和由于地球模型的理想化造成的地球模型误差和坐标转换误差以及CCD相机本身的视场畸变,而且从另一个角度提出了弹丸落点的坐标测量,为以后使用更先进,在特殊环境下的落点测量方法,得到更高精度的位置信息提供了新思路。
-
图 1 为水面整体测量示意图。它由打击点、靶标和两测量船(测量船上有连续变焦镜头的CCD相机)等组成。利用卫星导航系统GPS测得船上任意三定点A、B、C和两测量船Q、Q′的大地坐标(大地测量中以参考椭球面为基准面建立起来的坐标系,地面点的位置用大地纬度、大地经度、大地高度表示),分别对应(Bi,Li,Hi),i=0,1,2、(B,L,H)和(B′,L′,H′),利用大地坐标和地心直角坐标系之间的换算关系[ )。关系式为: 式中:为卯酉圈的曲率半径, 为子午椭圆的第一偏心率,6 ,7 ],算出在地心直角坐标系下A、B、C三点的坐标(xi,yi,zi),i=0,1,2及Q、Q′的坐标(xq,yq,zq)、(x′q,y′q,z′qa,b 分别为椭球体的长短半径。在地心坐标系下,利用空间几何求出任意定点 A,B,C分别对于测量船Q,Q′的俯仰、方位角度信息Ai、Ei。对于测量船Q : 式中,i=0,1,2。同理对于测量船Q′有方位和俯仰角A′i和E′i 。本文利用CCD相机建立靶船上定点角度信息和对应点在CCD图像上的坐标间的关系函数,这样可以有效修正相机镜头畸变等造成的误差[8]。用于拍摄定点和落点的CCD相机阵面由整齐排列、大小相同的像元组成,构成了二维靶面直角坐标系。定点和目标落点成像在相邻的几个乃至几十个像元上,其像面坐标可由像点质心确定,以像元为单位。
由于CCD测量都是从片子的一个角为起点直接读数,为了计算方便,假设一个近似光心坐标为(512,512)(若选取1 024×1 024像元数的片子),从而所有的像面坐标值都是直接读数减去光心位置的数值。由于CCD相机观测视场很小,视场畸变和像差的影响甚微,为了把像面坐标(X,Y)和角度坐标(A,E)联系起来,本文只考虑两坐标系的原点和坐标轴指向,以及两坐标轴的比例尺的不同,可以用简单的线性关系底片常数模型建立变换关系。图 2为角度坐标和图像坐标的几何关系。
图 2 角度坐标和图像坐标的几何关系
Figure 2. Geometrical relationship between angular coordinates and image coordinates
由 图 2 可知,角度坐标O-AE的A轴与E轴分别和图像坐标O′-XY的X轴和Y轴夹角为α,β。由于X轴可能不垂直于Y轴导致可能α不等于β。现假设在图像坐标系中X、Y轴的比例尺分别为RX,RY,且O′在角度坐标系下的坐标为(a,d),并且某一定点的图像坐标为(X,Y ),那么对应的角度坐标为:设 RX cosα=b,-RY sinβ=c,RX sinα=e,Ry cosβ=f ,则有:若假定图像坐标中 X 轴经修正后严格地垂直于Y轴(α=β)两坐标的比例尺相同RX=RY ,则上式可以化简为: 此时,RX cosα=Ry cosβ=b ,RY sinβ=RX sinα=-c 。这就是两套参系统的关系传递函数,也就是由定点(Xi,Yi,Ai,Ei)i=0,1,2解算底片常数a,b,c,d的过程。对于测量船Q把A,B两定点分别对应的(X0,Y0,A0,E0)和(X1,Y1,A1,E1) 代入上式,求得a,b,c,d的值。进而根据目标落点M的图像坐标(Xm,Ym) 可以求得对于测量船Q的方位、俯仰角度信息(A,E),同理求出测量船Q′处M落点的(A′,E′)。
-
在用两传感器进行交会测量的过程中由于传感器本身的结构差异、成像原理、时间同步和目标差异,加上外部工作环境的作用,两测量船上的传感器观测同一目标时两光轴并不会相交,而是两条空间异面直线[ 9 ]。利用异面交会算法(用两测量船与目标落点异面光轴的中垂线段估计落点位置),标算出M 的空间坐标信息。两测量船相应于目标落点的光轴的方向向量分别为: vq =(mq,nq,pq)=(cosE cosA ,sinE ,cosE sinA )和 V′q =(m′q,n′q,p′q)=(cosE ′cosA ′,sinE ′,cosE ′sinA ′)。公垂线的方向向量v=(A,B,C)=vq×V′ M1、M2两点。Q 。 假设公垂线段与两光轴分别相交于Q,Q′ 测量船的光轴直线方程表示如下:根据空间几何关系可以得出M1和M2两点的空间坐标M1(xq+mqtq,yq+nqtq,zq+pqtq 和M2(x′q+m′qt′q,y′q+n′qt′q,z′q+p′qt′q)。
那么 M 点的坐标表示为: 式中:mq,nq,pq分别为测量船Q对于目标落点的光轴方向向量分量;m′q,n′q,p′q为测量船′Q对于目标落点的光轴方向向量分量; A,B,C 为公垂线的方向向量分量。
其算法流程表示如下:
根据中国科学院长春光机所吴能伟和陈涛等人对异面交会误差精度的分析,交会角度的大小直接影响目标定位精度。交会角在60°~90°时,定位误差逐渐减小,在90°到达最小,随着交会角度的增大(90°~120°)而增大。所以实际工程中交会角选取90°最为理想[10,11,12]。
-
由测量方案可知,主要误差来源是GPS定位误差,CCD图像质心精度造成的底片常数模型误差,异面交会误差。GPS是目前定位精度最高的导航系统,可提供动态目标的高精度位置、速度和时间信息,其定位精度可以达到0.05 m。底片常数模型只是反应两套坐标系间坐标原点、坐标轴方向和尺度的转换关系,为了减少底片常数,减少定标点数,使“四常数模型”可以采用,由于CCD像元线尺寸仅是近似值,x方向和y方向的尺寸不可能做到严格相等,需对“四常数模型”要求的量度坐标x轴与y轴严格垂直和比例尺进行修正[13,14]。
质心坐标对底片常数的解算精度产生影响,定点图像位置精度受CCD单个像元的空间角分辨率影响,因为CCD像元尺寸不能做得很小(其一般为10 μm左右,小的4~6 μm)限制了测量精度,从而影响定位精度。对于视场覆盖为 L ,距目标点为R,有效像元数为N×N 的相机,其空间分辨率ξ 为(每个像元对应的视场角的大小): 式中:N 为CCD芯片单边的有效像元数目,Φ =2arctan(L/2R )为视场观测角度。 由工程经验和质心确定算法(典型的高斯曲面拟合细分定位算法)可得,质心偏差可达0.5个像素,即0.5ξ 。目标落点的位置误差是由交会测量中的交会角和落点相应于两传感器的角度信息决定。落点角度误差主要有质心精度和GPS精度共同决定。根据误差合成理论,下面给出具体的误差分析过程。
-
在利用异面直线交会时,必须选择合适的工作基点,尽量使交会图形满足两观测角相等且交会角在90°附近。从交会测量的坐标式(6)可知,落点的3个坐标分量分别是 xq,yq,zq,x′q,y′q,z′q ,A,A′,E,E′的函数。目标落点相对于两测量船的方位角和俯仰角误差σA,σA′,σE,σE′ 可以由图像坐标的质心误差也就是底片常数模型误差分析得出,GPS引进的测量船和定点的位置误差为σGPS 。考虑到测量误差一般服从高斯分布,根据误差传递理论得出每个元素对落点空间位置标准差的影响,σx,σy,σz分别为: 式中:通过连续变焦镜头改变CCD相机的空间分辨率,根据底片常数模型解算目标落点的角度的过程中可以得出角度误差,公式表示如下: 式中: i =0,1。同理可得σ2A′、σ2E′。由式(9),(10)可以看出A,E,A′,E′的标准差主要由定点和目标点图像坐标的质心误差和随定点空间角度信息引入的位置误差线性决定。最终传递到目标位置误差中,落点位置误差与各误差源也是线性相关的。
-
本文进行了实际观测和标定实验。实验中采用相位差分GPS测量点位位置。Q和Q′测量船使用的相机型号是滨松ORCA-flash 4.0,相机的图像分辨率为2 048×2 048,像素尺寸为6.5 μm×6.5 μm,平均量子效率为70%,其空间分辨率随着对应的视场观测角度变化。观测前,需将至于两测量船上的相机探头大致指向靶标船,使靶标船尽可能的位于CCD中心,调整CCD使视场内定点和目标点的轨迹尽可能的与CCD靶面边平行。实验过程中弹丸速度过快,对于图像方法拾取弹丸较困难,因此采用测量弹丸溅起水柱的方法替代弹丸落点。表 1为经式(1)变换后地心坐标系下点位的坐标。
表 1 靶标船上定点以及测量船的地心坐标
Table 1. Geocentric coordinates of points and measuring ships
No. x/m y/m z/m A -2 481 438.920 4 172 027.248 4 123 447.395 B -2 481 451.865 4 172 019.178 4 123 447.768 C -2 481 464.340 4 172 011.337 4 123 448.191 Q -2 479 586.645 4 176 130.661 4 120 427.098 Q′ -2 485 831.796 4 172 054.275 4 120 791.172 经连续拍摄得到一组实际CCD图像,得到一组连续的现场图像,图 3为t时刻的一幅图像。人工识别出t时刻的图像,得到其中定点和目标点的像元位置坐标信息,参见表 2。
表 2 定点和目标点的像元位置坐标
Table 2. Image coordinates of specified points and projectile point
No. 测量船Q 测量船Q′ 横坐标 纵坐标 横坐标 纵坐标 x/pixel y/pixel x′/pixel y′/pixel A -25.4 3 -23.9 0 B 2 3 1.5 0 C 28.3 3 26.3 0 M -10.8 3 83.1 0 由计算结果可知Q测量船的底片常数a,b,c,d分别为121.7、1.6、0和-49.2。同理Q′测量船的底片常数分别为31.2、0.7、0.8和-0.4。根据底片常数模型可以反推出目标落点分别相对于两测量船的方位、俯仰角。由式(4)~(6)计算水面落点的位置,式(8)~(10)计算目标落点的位置误差。表 3给出了水面落点位置坐标值和标准差(交会角为90°的情况下)。
表 3 靶水面落点坐标和相应方向上的误差
Table 3. Result of value of the projectile point and its error
真值/m 误差/m x -2 481 459.4 2.8 y 4 171 993.0 4.9 z 4 123 469.6 4.3 -
针对底片模型的角度误差考察了质心误差和位置定位误差所占的权重。根据工程实际设备指标,分析了方位角和高低角方差与相机空间分辨率平方项、GPS精度平方项之间的关系(如图 4)。
由图 4可以看出,在用底片常数模型求解角度的过程中,误差源对方位角和俯仰角误差的影响是一致的。目标点的方位角和俯仰角的方差与相机的空间分辨率的平方和GPS的精度的平方成线性关系。因此可以根据具体工程需要和成本要求选择相应的测量设备。
-
本文根据水面落点的特殊前提以及高速运动目标,提出了依据动态像面测量坐标的方法,并介绍了动态像面以及异面交会算法的原理。然后对研究过程中误差的主要来源进行分析,并分析了落点误差与各个误差因素的关系。本文采用地心坐标系有效避免了地球模型的模型误差和坐标系转换的转换误差;底片常数模型避免了相机镜头本身的畸变误差和测量船的姿态误差;异面交会方法克服了前方交会由于作用距离产生的误差。
在地心坐标系下运用底片常数模型方法进行的目标落点角度信息的估计,在交会角为90°的情况下,位置误差x方向上最高能达到2.8 m。实验结果表明,在保证交会误差最小的情况下,本文的方法实现了水面特殊条件下较为理想的落点位置测量。
Measurement of the location of projectile burst on water surface by dynamic image method and its error analysis
-
摘要: 为了测量弹丸水面落点的位置, 建立了基于CCD相机动态像面的测量模型。该模型通过CCD相机辅助采集相关点位的图像信息, 对水面落点的位置函数以及误差进行了研究。首先, 利用空间几何获得靶船上三定点相对于测量船的方位、俯仰信息。接着, 结合t时刻观测图像上定点的像面坐标, 运用底片常数模型建立像面坐标和角度信息两套参量之间的关系函数, 从而得到目标落点的方位、俯仰信息, 再利用异面交会法计算出目标落点位置。最后, 分析了目标落点位置的误差来源(质心误差, 位置误差)、误差以及各误差源与位置坐标之间的关系。实验结果表明:在测量船位置精度达到0.05 m, 图像质心定位精度达到0.5 pixel时, 在最小交会误差的情况下, 目标落点的位置测量误差分别为2.8, 4.9, 4.3 m。Abstract: In order to measure the value of projectile points on the water surface accurately, a model of dynamic images measurement based on CCD camera is established. The function of the value of the surface point and its measuring error are analyzed by using the collected images information of the relative points taken by CCD camera secondarily. First, with the knowledge of spatial geometry, the azimuth and the elevation angles of three designated points on the target ship relative to the measuring ship are obtained. Then, combined with image coordinates of the specified points at t point, the relationship function of two sets of parameters between the image plane coordinate and angle information is presented by taking the constant model, so the azimuth and the elevation angles information of the projectile point is obtained. The function of three-dimensional projectile burst coordinates is presented by non-coplane intersection method. Finally, the coordinate's errors from the evaluation of the pinpoint of the image and the location measurement accuracy are discussed and the relationship is investigated. Experimental results indicate that the three-dimensional projectile burst coordinate's measuring error are 2.8, 4.9, 4.3 m, respectively, in the condition of the minimum intersection error, when the position accuracy is up to 0.05 m and the image pinpoint accuracy to be able to 0.5 pixel.
-
表 1 靶标船上定点以及测量船的地心坐标
Table 1. Geocentric coordinates of points and measuring ships
No. x/m y/m z/m A -2 481 438.920 4 172 027.248 4 123 447.395 B -2 481 451.865 4 172 019.178 4 123 447.768 C -2 481 464.340 4 172 011.337 4 123 448.191 Q -2 479 586.645 4 176 130.661 4 120 427.098 Q′ -2 485 831.796 4 172 054.275 4 120 791.172 表 2 定点和目标点的像元位置坐标
Table 2. Image coordinates of specified points and projectile point
No. 测量船Q 测量船Q′ 横坐标 纵坐标 横坐标 纵坐标 x/pixel y/pixel x′/pixel y′/pixel A -25.4 3 -23.9 0 B 2 3 1.5 0 C 28.3 3 26.3 0 M -10.8 3 83.1 0 表 3 靶水面落点坐标和相应方向上的误差
Table 3. Result of value of the projectile point and its error
真值/m 误差/m x -2 481 459.4 2.8 y 4 171 993.0 4.9 z 4 123 469.6 4.3 -
[1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] -