摘要
机载面阵摆扫式成像系统需要足够高的帧间重叠率来避免因扫描过程中视轴的控制误差而导致的漏扫,这制约了成像系统在大速高比条件下的适应能力,限制了系统的成像幅宽。论文提出了一种视轴路径规划模型与重叠率计算的新方法,通过像球面投影的几何分析方法优化设计了成像系统视轴的扫描路径,并依据像球面中像面投影的方位旋转角度,建立了一种重叠率计算的新模型。经理论分析与仿真实验验证,在40°侧向凝视成像时,视轴稳定精度由2.93°提高到0.15°,相同速高比、单行5帧摆扫成像时,作业效率提高约32 %。论文研究工作对推动机载摆扫成像技术向着宽视场、大速高比、高分辨率方向的进一步发展具有重要意义。
如今,机载光电成像技术一直在朝着大视场、高分辨率成像的方向发展,越来越多的航空相机采用了摆扫成像技
在视轴控制的相关研究当中,关于目标定位和像移研究的文献很多,几乎都采用坐标变换法来进行计
首先通过建立五个坐标系,将载机位移运动、载机姿态变化运动和系统运动执行机构摆扫运动统一到同一旋转参考系,之后使用一种基于像球面投影的几何分析方法,结合像移补偿与扫描运动,系统设计了相机视轴的运动路径,再针对两轴下视相机,通过坐标变换法计算出执行机构运动角度,最后依据像球面中像面投影的方位旋转角度设计了一种有效重叠率的计算方法并进行了实验验证与误差分析。
地面地理坐标系与载机地理坐标系的x、y、z轴方向相同,分别指向正北、正东、地心方向,S点位于载机质心,D点位于载机附近物点平面,如

图1 地面地理坐标系D-DxDyDz与载机地理坐标系S-SxSySz
Fig. 1 Geographic coordinate system of the ground D-DxDyDz and geographic coordinate system of the vehicle S-SxSySz
载机机体坐标系Ax轴指向机头方向,Ay轴指向机右翼方向,Az轴指向载机正下方,坐标系原点A、C与S点重合,S坐标系先后绕z、y、x轴逆时针旋转(内旋

图2 载机机体坐标系A-AxAyAz与成像系统坐标系C-CxCyCz
Fig. 2 Coordinate systems of the vehicle body A-AxAyAz and coordinate systems of the camera C-CxCyCz
另外建立摄影坐标系来进行视轴运动路径和重叠率的分析。设载机进入航带线时的初始方位角为ψ0,摄影坐标系由S坐标系绕z轴逆时针旋转ψ0得到,每条航带开始时重新初始化ψ0。成像系统视轴在R坐标系下的姿态角用俯仰φ、横滚ω、方位κ表示,即R坐标系先后绕y、x、z轴逆时针旋转(内旋)φ、ω、κ后可得到C坐标系,如

图3 摄影坐标系R-RxRyRz
Fig. 3 Coordinate systems of photography R-RxRyRz

图4 摆扫成像示意图
Fig. 4 The schematic diagram of whisk-broom imaging
视轴运动路径的规划,要同时考虑载机的飞行平移运动、载机的姿态旋转运动以及成像系统运动执行机构的摆扫旋转运动,又因为视轴运动为中心旋转运动,故需要将以上运动都归转到同一坐标系内同一中心的旋转运动。此外航空对地摄影属于中心平面投影,由于投影变换带来的形变,投影在地平面的像面均为不规则的四边形,很难进行重叠率的计算。所以本文使用一种基于像球面投影的视轴运动路径规划与重叠率计算方法。在摄影坐标系中,像面绕投影中心的三维旋转组成一个以投影中心为球心的像球面,将飞行过程中地物目标相对投影中心的平面运动转化为角度旋转运动,结合摆扫运动设计出成像系统的视轴运动路径,最终通过坐标变换结合载机姿态角解算出执行机构的运动角度,来进行视轴控制。此外在像球面上的像面投影大小形状均一致,也更加适合进行重叠率的计算。
视轴OU与其在坐标系xOz、yOz平面上的投影以及Oz轴围成了a、b、c、d四个角,如

图5 视轴OU与其四个定位角
Fig. 5 LOS OU and its four positioning angles
为了方便重叠率的设计,像移补偿运动不能影响扫描轴的扫描中心。此外,像移运动主要发生x轴方向上,为了方便之后的简化运算,视轴运动设计采用PR旋转系统,即使用a、d角来表示,对应摄影坐标系为φ、ω角。如

图6 R坐标系下PR旋转系统
Fig. 6 PR rotary system in R coordinate systems
经几何计算分析,可得视轴俯仰角及角速度运动表示为:
, | (1) |
视轴横滚角及角速度运动表示为:
, | (2) |
其中ω0、φ0为初始角度,ω’为扫描角,ωtotal为总扫描行程,Tsacn为一个摆扫成像周期,目标相对载机移动速度vR为(vs为载机速度):
. | (3) |
我们发现ω角计算公式较为复杂,且与φ角耦合,不利于实时机上运算及运动控制,影响控制响应速度和两轴控制独立性,故此接下来尝试对计算公式进行简化并说明其中的误差:
1. φ角变化会引起的ω角度变化,但由于一般实际当中φ角变化幅度都不大,约±10°以内,cosφ∈(-0.9848,0.9848),cosφ≈cosφ0≈1,故另cosφ=cosφ0,对ω角(前半部分)造成的误差不超过0.8%;
2. ωtotal=arctan(tan(ωtotal,φ=0)×cosφ),ωtotal也会受φ角变化影响,若另cosφ=1,对ωtotal角造成的误差不超过1%。
故以上两处简化,对最终ω角计算结果的误差不会超过1%,工程使用可以进行简化,简化后的扫描公式为:
. | (4) |
已知载机初始方位角ψ0以及姿态方位角ψ、俯仰角θ、横滚角Φ,和规划视轴俯仰角φ、横滚角ω,以及待求的执行机构横滚角α、俯仰角β及视轴方位角κ,通过坐标变换列出矩阵等式:
. | (5) |
经过变换可得坐标变换矩阵R为:
, | (6) |
最后可求得机构横滚角α,俯仰角β以及视轴方位角κ为:
. | (7) |
由上节计算可以看出,两自由度补偿系统无法完全补偿三自由度视轴,最后留下一个视轴方位角无法补偿,但不影响视轴指向,只是影响视轴旋转,当视轴发生旋转时,必然导致像面在像球面投影无法平行对齐。为了方便讨论重叠率,我们在此对旋转后的像面投影进行裁切,使其再次对齐,并计算其有效视场角和有效重叠率。
L、W为成像系统翼展、飞行方向视场角(假设飞行方向沿短边),当一帧图像有角度k的旋转时,我们对其进行最大化切割,如

图7 k角旋转后像面投影的最大化切割示意图
Fig. 7 The schematic diagram of maximized cutting of image plane projection after k-angle rotation
通过简单的几何推导可得,其有效视角变为:
, | (8) |
裁切对齐后的相邻像面投影如

图8 裁切对齐后重叠率计算示意图
Fig. 8 The schematic diagram of calculation of overlap ratio after cropping and alignment
相邻帧之间的有效重叠率为:
. | (9) |
以某型吊舱长焦相机为实验对象,其单帧视场角为20.18°×15.21°,在飞行高度6 000 m,速高比0.042,单行摆扫成像5帧情形下,截取其载机实际作业时单条航带内的姿态角度变化曲线,如

图9 载机姿态角度变化曲线
Fig. 9 Aircraft attitude angle change curve
首先通过定点凝视(初始视轴角横滚-40°、俯仰5°)成像模式来验证方法的视轴稳定性。

(a)

(b)
图10 视轴稳定效果图(a)本文方法稳定效果,(b) 传统方法稳定效果
Fig. 10 Stabilizing effects of LOS (a) Stabilizing effects of this method,(b) Stabilizing effects of traditional method
接下来进行摆扫成像工作模式下的重叠率设定验证。经计算,在当前姿态扰动输入下,|κ|波动范围位于4.60°之内,要保证有效重叠率大于0,即成像工作时无漏扫发生,翼向、航向重叠率应至少设置为6.37%、9.67%,按设定重叠率模拟成像效果如

(a)

(b)
图11 模拟成像效果图(a) 本文方法成像效果,(b) 传统方法成像效果
Fig. 11 Simulate imaging effects (a) Imaging effects of this method, (b) Imaging effects of traditional method
以上模拟仿真在理想情况下对视轴角度和重叠率进行了计算。但实际成像系统还存在系统误差:如模型误差,以及随机误差、测量误差和控制误差等,势必会对实际视轴角与重叠率造成影响。现以飞行某一|κ|值较大的时刻数据进行分析,将模型误差、测量误差、控制误差全部纳入分析,均以正态分布函数来模拟参数误差分布,然后使用蒙特卡洛法来计算最终结果的误差分布,以此来估计实际情况下的工作状态。
计算过程如
误差源 | 均值µ(°) | 标准差σ(°) | 误差类型 | |
---|---|---|---|---|
输入1 | 姿态方位角 ψ | -3.58 | 0.03 | 随机误差:来自POS_AV610(SPS)数据 |
姿态俯仰角 θ | 2.12 | 0.01 | ||
姿态横滚角 Φ | -0.52 | 0.01 | ||
规划视轴俯仰角 φ | 5.00 | 0.00 | 系统误差:模型误差 | |
规划视轴横滚角 ω | -40.00 | 0.07 | ||
输出1 | 规划机构横滚角 α | -39.70 | 0.07 | |
规划机构俯仰角 β | -0.10 | 0.02 | ||
视轴方位角 κ | -4.59 | 0.02 | ||
输入2 | 姿态方位角ψ | -3.58 | 0.03 | 随机误差:来自POS_AV610(SPS)数据 |
姿态俯仰角θ | 2.12 | 0.01 | ||
姿态横滚角Φ | -0.52 | 0.01 | ||
实际机构横滚角α’ | -39.70 | 0.09 | 随机误差:目标机构角叠加系统实际控制误差 | |
实际机构俯仰角β’ | -0.16 | 0.04 | ||
输出2 | 实际视轴俯仰角φ’ | 4.92 | 0.06 | |
实际视轴横滚角ω’ | -40.00 | 0.09 | ||
实际视轴方位角 κ’ | -4.54 | 0.05 |
经以上计算分析,最终摆扫过程当中的视轴定位精度(µ残差加2倍σ):俯仰角误差(-0.20°~0.04°)、横滚角误差(-0.18°~0.18°),|κ|值范围位于4.64°之内,综合考虑视轴误差与k角旋转,重叠率应设置为6.42%、9.75%,所以考虑误差来讲,相比传统方法20%重叠率仍可提高约31.96%以上工作效率。
为了进一步提高机载摆扫成像系统的作业效率,本文针对传统视轴控制方法稳定精度低以及重叠率大问题,在综合考虑计算精度和工程可实现性的情况下重新设计了视轴运动路径控制方法,并针对重叠率进行了精确计算。经过实验验证,本方法将视轴定位精度由2.93°提高到0.15°,并通过精确设定重叠率,可将作业效率提高约32%。
值得注意的是,本文重叠率设定仍有较大冗余,因为只有少数情况才会出现较大的k角旋转,实际成像过程当中可酌情减小重叠率设定,并在机上进行实时k值监测,防止漏扫情况的发生,可以进一步提高系统的作业效率。本文控制方法已在某型吊舱相机上得以实现,对推动机载摆扫成像技术向着宽视场、大速高比、高分辨率方向的进一步发展具有重要意义。
References
WANG C R, YANG L F, CAO X, et al.. Recent progress of airborne infrared remote sensing technology in SITP[J]. J. Infrared Millim. Waves, 2022, 41(1): 007.王崇儒, 杨利峰, 曹汛, 等. 航空红外光电遥感技术最新进展[J]. 红外与毫米波学报, 2022, 41(1): 007. [百度学术]
LIU Y B, WANG Y K, HAN C G, et al.. Aviation wide field thermal imaging technology based on array detector sweep-stop pattern[J]. J. Infrared Millim. Waves, 2017, 36(3): 330~335. [百度学术]
刘毓博, 王义坤, 韩贵丞, 等. 基于摆停模式的航空大视场面阵热成像技术[J]. 红外与毫米波学报,2017,36(3):330~335. 10.11972/j.issn.1001-9014.2017.03.014 [百度学术]
LI S, HU Y, GONG C L, et al.. A step-by-step geometric correction and error analysis of swing-swept array thermal infrared aerial image[J]. J. Infrared Millim. Waves, 2020, 39(2): 242~250. [百度学术]
李赛, 胡勇, 巩彩兰,等. 面阵摆扫热红外航空影像分步几何校正方法[J].红外与毫米波学报,2020,39(2):242~250. 10.11972/j.issn.1001-9014.2020.02.012 [百度学术]
SUN H. Target localization and error analysis of airborne electro-optical platform[J]. Chinese Optics, 2013, 6(6): 912. [百度学术]
孙辉. 机载光电平台目标定位与误差分析[J]. 中国光学,2013, 6(6): 912. 10.3788/CO.20130606.912 [百度学术]
WANG J Q, JIN G, YAN C X. Orientation error analysis of airborne opto-electric tracking and measuring device[J]. Optics and Precision Engineering, 2005, 13(2):105-116. [百度学术]
王家骐, 金 光, 颜昌翔. 机载光电跟踪测量设备的目标定位误差分析[J]. 光学精密工程, 2005, 13(2):105-116. 10.3321/j.issn:1004-924X.2005.02.001 [百度学术]
SUN H, ZHANG S M. Computation model and error budget for image motion of aerial imaging system[J]. Editorial Office of Optics and Precision Engineering, 2012,20(11): 2492-2499. [百度学术]
孙辉, 张淑梅. 机载成像系统像移计算模型与误差分析[J]. 光学精密工程, 2012,20(11): 2492-2499. 10.3788/OPE.20122011.2492 [百度学术]
WANG J Q, YU P, YAN C X, et al.. Space Optical Remote Sensor Image Motion Velocity Vector Computational Modeling[J]. Acta Optica Sinica, 2004, 24(12): 1585.王家骐, 于平, 颜昌翔, 等. 航天光学遥感器像移速度矢计算数学模型[J]. 光学学报, 2004, 24(12): 1585. [百度学术]
WANG Y K, HAN G C, QI H X, et al.. Investigation of image motion compensation technique based on real-time LOS tracking[J]. J. Infrared Millim. Waves, 2015, 34(6): 0757.王义坤, 韩贵丞, 亓洪兴, 等. 基于实时视轴跟踪的像移补偿技术研究[J]. 红外与毫米波学报, 2015, 34(6): 0757. [百度学术]
Hong S Y, Wang Y K, Han G C, et al.. Study on Real-Time LOS Stabilization Technique Based on Dual Reflectors in the Airborne Remote Sensing Imaging System[J].Aero Weaponry, 2020, 27(5): 86-90.洪声艺, 王义坤, 韩贵丞, 等. 基于双反射镜的航空遥感成像系统实时视轴稳定技术研究[J]. 航空兵器, 2020, 27(5):5. [百度学术]
ZHOU G, ZHAI L P, WANG D J, et al.. Research of the Scanning Angle of Panoramic Aerial Camera[J]. Optics and Precision Engineering, 2008, 16(12):2473-2478. [百度学术]
周刚, 翟林培, 王德江, 等. 全景式航空相机的扫描角研究[J]. 光学精密工程, 2008, 16(12):2473-2478. 10.3321/j.issn:1004-924X.2008.12.025 [百度学术]
LI J, HUANG H T, XIU J H, et al.. Effect and compensation of overlap influenced by flight parameter of oblique aerial camera[J]. Optics and Precision Engineering, 2020, 28(6):1254-1264. [百度学术]
李军, 黄厚田, 修吉宏, 等. 载机飞行参数对倾斜成像重叠率影响及补偿[J]. 光学精密工程, 2020,28(6):1254-1264. 10.3788/OPE.20202806.1254 [百度学术]