Analysis on effect of platform rotation on SAR imaging
-
摘要: 分析了合成孔径雷达(SAR)平台偏航、俯仰和横滚3个方向转动对成像的影响。建立了平台转动模型，从雷达照射区域的变化、回波的幅度调制和回波的相位调制等3个方面展开分析，定量分析了平台转动对成像幅宽和方位聚焦的影响。分析结果表明，平台转动会产生成像幅宽损失，影响方位聚焦，并会改变目标的方位位置。分析结果可以为SAR系统设计和算法论证提供一定的参考。Abstract: The effects of platform rotation in three directions, yaw, pitch and roll, on synthetic aperture radar(SAR) imaging were analyzed. The platform rotation model was built, the analyses were conducted into three aspects: the change of the radar radiation area, the amplitude modulation and the phase modulation of the echo, and the effects of platform rotation on the imaging width and azimuth focusing were analyzed quantitatively.Simulation results show that platform rotation reduces the imaging width and leads to azimuth defocusing, so as to change the target azimuth position, which can provide guidelines for SAR system design and algorithm development.
-
Keywords:
- SAR imaging /
- platform rotation /
- phase modulation /
- amplitude modulation
-
引言
在兵器靶场测试研究领域，对弹丸各项参数检测时能够得到弹丸清晰的飞行图像是非常重要的[1-2]。传统观测弹丸的飞行姿态有2种方法：一种是沿着外弹道方向放置多台高速CCD相机，通过分析每台摄像机单独拍摄弹丸图像，经过后续复杂的处理技术，得到较为完整的弹丸图像[1]。此方法受单个CCD相机视场范围的限制，且对多个相机所拍摄的图像进行分析时还涉及复杂的处理技术，更重要的是长距离测量时要花费很多的财力和物力[3-5]；另一种方法是将单个CCD相机系统固定在跟踪机架上，通过实时旋转CCD相机镜头来跟踪弹丸，但使用时，很难将设备用于测量小型高速飞行的弹丸。为了弥补这2种方法的缺陷，国外提出一种通过转镜和CCD相机实现对弹丸的实时跟踪技术，固定CCD相机位置，将它对准反射镜，通过拍摄反射镜中弹丸图像，以实现对高速弹丸的跟踪拍摄[6]。目前国外对此项技术的保密性非常高，国内相关产品还处于研发阶段，为了弥补我国在相关技术的空白，对于转镜跟踪测量系统的研究是非常有意义的[7]。基于高速CCD视场中点的转镜跟踪系统是通过计算机控制转镜转动来跟踪弹丸，同时弹丸的飞行姿态反射到CCD相机后，拍摄一组序列图像，从而实现对弹丸的同步跟踪[8-9]。因此，在跟踪系统中，建立弹丸运动以及转镜转动的数学模型对研究弹丸的飞行姿态有着重要的作用[10]。本文将研究基于高速CCD视场中点的转镜跟踪系统，主要包括建立弹丸和转镜的数学模型以及相关参数对系统的影响。
1 系统组成和工作原理
1.1 系统的结构组成
基于高速CCD相机视场中点的跟踪系统主要包含5个组成部分，如图 1所示。在系统能有效跟踪弹丸的视场前方设有测量弹丸速度的天幕靶，当弹丸穿过天幕靶后就能在计算机控制系统中获得弹丸的速度，此时计算机控制系统通过计算弹丸进入转镜的有效视场时间驱动伺服系统控制转镜在相应的时间段内转动。在整个跟踪过程中，转镜转动将运动弹丸的飞行姿态反射到CCD相机中，CCD相机拍摄的弹丸图像存储在计算机中，从而实现CCD相机对高速运动弹丸的同步跟踪测量。
1.2 系统的工作原理
一般情况下，转镜跟踪系统能跟踪90°视场的弹道线。跟踪系统在开始跟踪弹丸前有5°的缓冲角让反射镜预加速到一定速度，在跟踪结束时也有5°的缓冲角让反射镜停止。基于高速CCD相机视场中点的跟踪系统，工作原理图如图 2所示。当弹丸穿过天幕靶2后到达转镜-50°视场时，转镜缓慢加速转动。当弹丸到达-45°视场角时，转镜的条件满足同步跟踪弹丸的条件，系统开始同步跟踪弹丸运动。当弹丸到达+45°视场角时，反射镜开始缓慢减速，当到达转镜+50°视场时，系统停止工作，从而完成对90°转镜视场范围内弹丸的同步跟踪。
2 系统数学模型的建立
2.1 弹丸运动模型的建立
地球上任何飞行的物体都会受到地心引力和空气阻力的影响，在实际情况下弹丸发射后由于受到这2个力的作用，它不再是沿着弹道线的直线运动，而是作斜抛运动。由于弹丸的飞行速度很快，对弹丸实现实时跟踪的距离相比于整个运动过程较短，因此，在理论研究过程中，假设弹丸不受地心引力和空气阻力的影响，弹丸则沿着弹道线方向作匀速运动，且轨迹与地面平行。
假设弹丸沿弹道线方向作匀速直线运动，转镜视场中点跟踪弹丸的原理图如图 3所示。CCD相机主光轴与弹道线平行且与转镜相交于O点，OQ为相机主光轴到弹道线之间的垂直距离，记为H，垂足Q为坐标原点，规定弹丸沿弹道线运动的方向为x正方向，W点为要跟踪的弹丸，则弹丸距坐标原点的距离为|WQ|= S。CCD相机的视场角为2θ，转镜与主光轴的夹角为ϕ，OP=L1, AB段为转镜跟踪弹丸的有效范围。
由反射定理和几何关系并根据图 3可以得到弹丸运动模型：
$$ \begin{array}{*{20}{l}} {S\left( t \right) = \frac{H}{{\tan {\Phi _{\left( t \right)}}}} + } \\ {\;\;\;\;\;\;\;\;\;\;\frac{{\left( {0.5{L_1}\sin 2{\Phi _{\left( t \right)}}} \right) + H\cos {\Phi _{\left( t \right)}}{{\sin }^2}\theta }}{{\sin {\Phi _{\left( t \right)}} \times \sin \lambda \times \sin \beta }}} \end{array} $$ (1) 同时, A点坐标：
$$ {X_A} = \frac{H}{{\tan {\mathit{\Phi }_{\left( t \right)}}}} - \frac{{\left( {{L_1} + \frac{H}{{\sin {\mathit{\Phi }_{\left( t \right)}}}}} \right)\sin \theta }}{{\sin \beta }} $$ (2) B点坐标：
$$ {X_B} = \frac{H}{{\tan {{\Phi }_{\left( t \right)}}}} + \frac{{\left( {{L_1} + \frac{H}{{\sin {{\Phi }_{\left( t \right)}}}}} \right)\sin \theta }}{{\sin \lambda }} $$ (3) 由(2)、(3)式可得整个跟踪视场的范围为
$$ AB\left( t \right) = \left( {{L_1} + \frac{H}{{\sin {\mathit{\Phi }_{\left( t \right)}}}}} \right)\left( {\frac{{\sin \beta + \sin \lambda }}{{\sin \beta \sin \lambda }}} \right)\sin \theta $$ (4) (1) 式两边同时对t求导，则得出弹丸速度和转镜角速度的关系为
$$ \begin{array}{l} \frac{v}{w} = H\csc {\mathit{\Phi }_{\left( t \right)}}[ - 2\csc {\mathit{\Phi }_{\left( t \right)}} - \cot {\mathit{\Phi }_{\left( t \right)}} \cdot \\ \;\;\;\;\;\;\;\frac{{\sin \theta \left( {\sin \beta - \sin \lambda } \right)}}{{\sin \lambda \sin \beta }}] + \left( {{L_1} + H\csc {\mathit{\Phi }_{\left( t \right)}}} \right) \cdot \\ \;\;\;\;\;\;\;\frac{{\sin \theta \left( {{{\sin }^2}\lambda \cos \beta + {{\sin }^2}\beta \cos \lambda } \right)}}{{{{\sin }^2}\beta {{\sin }^2}\lambda }} \end{array} $$ (5) 式中：β=2ϕ(t)+θ；λ=2ϕ(t)-θ；Φ(t)=2ϕ(t)。
利用MATLAB对系统的数学模型进行仿真, 由于在对不同速度的弹丸进行跟踪时，系统开始工作的时间不同，则弹道宽度AB也不同，弹道宽度曲线随时间的变化曲线可由MATLAB仿真得到。当H=200 m, v=100 m/s时，转镜视场范围内的有效视场AB随时间变化的曲线如图 4所示。
![]()
图 4 转镜视场范围内有效跟踪范围随时间变化曲线Figure 4. Curves of effective tracking range over time within rotation mirror FOV根据图 4可以看出A、B点的坐标随着时间的变化而变化，系统能有效跟踪弹丸的范围也在随时间变化，它的曲线呈开口向上的抛物线，当弹丸到达转镜正上方时有效跟踪范围最小，且最小值约为14 m，远大于被观测高速目标的尺寸，可以满足系统跟踪弹丸的视场要求。
2.2 扫描速率的大小
扫描速率R定义为被测弹丸速度v和相机主光轴到弹道线的垂直距离H的比值，它可以表示为
$$ R = \frac{v}{H} $$ (6) 对于不同速度的弹丸，不同的相机位置，只要保证它们的比值相同，那么转镜的扫描速率都是相同的，可以采用同一扫描速率。联立(5)式和(6)式，可以得出在ϕ=45°时，转镜的最大角速度ωmax和R之间的关系为
$$ {\omega _{\max }} = \left| { - \frac{v}{{2H}}} \right| = \frac{R}{2} $$ (7) 由(7)式可知，转镜的最大角速度仅跟扫描速率相关。
在跟踪系统中，扫描速率R是有一定范围的，它跟CCD相机所能接收物体像大小的能力有关。它的最大值与电机的转速有关，最小值与CCD的视场范围有关。在对弹丸的跟踪过程中，选择合适的扫描曲线控制转镜转动，在CCD中对连续弹丸图像的成像效果有着极大的提高。
3 设置参数对跟踪效果影响
根据(1)式可知，同一个跟踪系统在跟踪弹丸时(即相机视场角和转镜转动角度相同)，CCD相机、转镜的所处位置对整个系统的跟踪效果有影响，也就是转镜到弹道线的垂直距离OQ(H)和CCD相机到转镜的距离OP(L1)对跟踪视场有影响。
3.1 弹道线到转镜的垂直距离H
根据(1)和(4)式可知H越大，能有效跟踪弹丸的范围就越大。如果想得到弹丸长时间的飞行轨迹，可以通过一直增大H来实现，但是转镜转动和相机的视场角都是有限制的，H不会一直增大。
对于图 3，现在假设相机的焦距为f，弹丸长度为d，弹丸成像在CCD上的长度为L, CCD的光敏面尺寸为l×ħ。根据物像关系可知：
$$ \frac{d}{L} = \frac{{{L_1} + H}}{f} $$ (8) 为了保证在跟踪过程中弹丸的像不脱离CCD的光敏面，一般要求弹丸长度占CCD光敏面长度的1/5~4/5，由此关系和(8)式可以得到转镜和相机的位置关系为
$$ \frac{{5df}}{l} - {L_1} \le H \le \frac{{5df}}{{4l}} - {L_1} $$ (9) 由(9)式可知：已知相机参数和弹丸长度，通过调整转镜到弹道线的距离H和相机到转镜的距离L1就可以使弹丸清晰成像。
3.2 最大离散速度的间隔值
假设天幕靶测得弹丸的速度为v0，CCD视场大小为l×ħ，弹丸的长度为d，在必须考虑地球引力对弹丸的影响时，弹丸脱离CCD视场范围的情况有2种，其示意图如图 5所示。
如图 5所示，曲线1是弹丸在水平方向上发生偏离的情况，曲线2是弹丸在竖直方向上发生偏离的情况。通过计算得到弹丸在水平方向上的离散速度为
$$ \Delta v = \frac{{\left( {3l + 5d} \right){v_0}}}{{5S}} $$ (10) 弹丸在竖直方向上的离散速度为
$$ \Delta v = \frac{ħ}{{\sqrt {\frac{{8\left( {{H_1} - H} \right)}}{8}} }} $$ (11) 式中：S为弹丸在CCD视场内飞行的距离；H1为弹丸偏移后距离相机主光轴的垂直距离。
4 系统跟踪效果误差分析
基于高速CCD相机视场中点跟踪弹丸的转镜系统，在弹道线方向跟踪弹丸时(图 3所示)，取L=300 mm, H=200 m, θ=2°，利用MATLAB得到的跟踪效果表如表 1所示。
表 1 相机中点跟踪效果表Table 1. Camera midpoint tracking effect table时间T/s 转镜角度/(°) 视场前点B/m 视场后点A/m 理论跟踪点/m 实际跟踪点/m Δ/m 2.0 66.5203 -174.0924 -200.2985 -187.1954 -187.1924 -0.003 2.5 63.1024 -135.9317 -157.4249 -146.6783 -146.6763 -0.002 3.0 57.8523 -87.8016 -105.0347 -96.4182 -96.4160 -0.0022 3.5 53.0235 -50.0291 -65.1764 -57.6028 -57.6012 -0.0016 4.0 45.0000 6.9946 -6.9946 0 0.002 -0.002 4.5 38.8521 50.9746 36.3204 43.6484 43.6463 0.0021 5.0 32.1578 104.9447 87.7175 96.3311 96.3283 0.0028 5.5 27.0357 155.8694 134.5268 145.1981 145.1956 0.0025 6.0 22.5278 214.0722 186.1262 200.0992 200.0962 0.003 由(1)式可知：在相机和电机的参数相同的条件下(即相机视场角和转镜的转动角度相同)，弹丸的位置只跟相机主光轴到弹道线的垂直距离H和OP的大小相关。弹丸的跟踪点坐标是按照(1)式的函数关系计算出的，根据函数的系统误差公式可以得到弹丸坐标的系统误差公式为
$$ \begin{array}{l} \Delta {S_{\left( t \right)}} = \left( {\frac{1}{{{\rm{tan}}{\mathit{\Phi }_{\left( t \right)}}}} + \frac{{\cos {\mathit{\Phi }_{\left( t \right)}} \times {{\sin }^2}\theta }}{{\sin {\mathit{\Phi }_{\left( t \right)}} \times \sin \lambda \times \sin \beta }}} \right)\Delta H + \\ \;\;\;\;\;\;\;\;\;\;\;\left( {\frac{{0.5\sin 2{\mathit{\Phi }_{\left( t \right)}} \times {{\sin }^2}\theta }}{{\sin {\mathit{\Phi }_{\left( t \right)}} \times \sin \lambda \times \sin \beta }}} \right)\Delta {L_1} \end{array} $$ (12) 已知ΔH=0.003 m, ΔL1=0.001 m，利用MATLAB得到弹丸坐标的系统误差随时间变化的曲线如图 6所示。
由上图看出系统的误差范围在±0.003 m的范围内变化，结合表 1所示的Δ值可知，测量误差在系统误差的范围内。也就是说，基于高速CCD相机视场中点的跟踪系统可以实现对弹丸的实时跟踪性能。
5 结论
通过建立高速CCD相机视场中点的转镜弹道跟踪系统的数学模型，并对其运动曲线进行了仿真，得出了能有效跟踪弹丸的视场范围。分析了主要参数转镜到弹道线的距离H和转镜到相机的距离L1对跟踪效果的影响，计算了扫描速率和速度的最大离散间隔值以及它们对成像效果造成的影响，进一步得到了转镜和相机空间布置的关系。在理想条件下，当弹丸飞行速度为100 m/s，转镜到弹道线的垂直距离H为200 m时，分析了系统存在的误差，结果证明该系统可以实现对弹丸的同步跟踪。
-
表 1 实验SAR系统参数
Table 1 Parameters of experimental SAR system
参数 数值 作用距离/km 80 方位波束宽度/(°) 3 俯仰波束宽度/(°) 1.8 方位积累时间/ms 18 
下载: 导出CSV
-
[1] 刘永坦.雷达成像技术[M].哈尔滨:哈尔滨工业大学出版社, 1999. LIU Yongtan. Radar imaging technology[M]. Harbin:Harbin Institute of Technology Press, 1999.
[2] 魏钟铨.合成孔径雷达卫星[M].北京:科学出版社, 2001. WEI Zhongquan. Synthetic aperture radar satellite[M]. Beijing: Science Press, 2001.
[3] ILIPPO Biondi. Low-rank plus sparse decomposition and localized radon transform for ship-wake detection in synthetic aperture radar images[J]. IEEE Geoscience and Remote Sensing Letters, 2018(15):117-121. http://www.researchgate.net/publication/321809987_Low-Rank_Plus_Sparse_Decomposition_and_Localized_Radon_Transform_for_Ship-Wake_Detection_in_Synthetic_Aperture_Radar_Images
[4] 王玉萍.融合改进直方图PDE和二维Tsallis熵多阈值SAR分割[J].应用光学, 2018, 39(6):839-848. http://www.yygx.net/CN/abstract/abstract11194.shtml WANG Yuping. Multi-threshold SAR segmentation based on improved histogram PDE and two-dimensional Tsallis entropy[J]. Journal of Applied Optics, 2018, 39(6):839-848. http://www.yygx.net/CN/abstract/abstract11194.shtml
[5] 李春升, 杨威, 王鹏波.星载SAR成像处理算法综述[J].雷达学报, 2013, 2(1): 111-122. http://d.old.wanfangdata.com.cn/Periodical/ldxb201301015 LI Chunsheng, YANG Wei, WANG Pengbo. A review of spaceborne SAR algorithm for image formation[J]. Journal of Radars, 2013, 2(1): 111-122. http://d.old.wanfangdata.com.cn/Periodical/ldxb201301015
[6] 郭华东, 李新武.新一代SAR对地观测技术特点与应用拓展[J].科学通报, 2011, 56(15): 1155-1168. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=kxtb201115001 GUO Huadong, LI Xinwu. Technical characteristics and potential application of the new generation SAR for Earth observation[J]. Chinese Science Bulletin, 2011, 56(15): 1155-1168. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=kxtb201115001
[7] SHAO Weizeng, YUAN Xinzhe, SHENG Yexin, et al.Development of wind speed retrieval from cross-polarization Chinese Gaofen-3 synthetic aperture radar in typhoons[J]. Sensors, 2018.18(2): 412. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=sensors-18-00412
[8] ZHOU Lifan, CHAI Dengfeng, YU Xia, et al. Interferometric synthetic aperture radar phase unwrapping based on sparse Markov random fields by graph cuts[J]. Journal of Applied Remote Sensing, 2018(12):1. http://www.researchgate.net/publication/322527966_Interferometric_synthetic_aperture_radar_phase_unwrapping_based_on_sparse_Markov_random_fields_by_graph_cuts
[9] LI Ning, NIU Shilin, GUO Zhengwei, et al.Data-Based Motion Compensation for High-Resolution Sliding Spotlight Synthetic Aperture Radar[J]. Sensors, 2018, 18(3):842. doi: 10.3390/s18030842
[10] 皮亦鸣, 杨建宇, 付毓生.合成孔径雷达成像原理[M].成都:电子科技大学出版社, 2007. PI Yiming, YANG Jianyu, FU Yusheng. Principles of Synthetic Aperture Radar[M]. Chengdu:UESTC Press, 2007.
[11] 张澄波.综合孔径雷达:原理、系统分析与应用[M].北京:科学出版社, 1989. ZHANG Chengbo. Synthetic aperture radar: principle, system analysis and application[M]. Beijing:Science Press, 1989.
[12] 周峰, 邢孟道, 保铮.一种无人机机载SAR运动补偿的方法[J].电子学报, 2006, 34(6): 1002-1007. doi: 10.3321/j.issn:0372-2112.2006.06.006 ZHOU Feng, XING Mengdao, BAO Zheng. A method of motion compensation for unmanned aerial vehicles borne SAR[J]. Acta Electronica Sinica, 2006, 34(6): 1002-1007. doi: 10.3321/j.issn:0372-2112.2006.06.006
[13] 张雨轮, 张涛, 李涛, 等.基于惯导的机载斜视SAR运动补偿研究[J].火控雷达技术, 2014, 43(1): 18-21, 41. doi: 10.3969/j.issn.1008-8652.2014.01.005 ZHANG Yulun, ZHANG Tao, LI Tao, et al. Study on IMU based motion compensation for airborne squint SAR[J]. Fire Control Radar Technology, 2014, 43(1): 18-21, 41. doi: 10.3969/j.issn.1008-8652.2014.01.005
[14] 保铮, 邢孟道, 王彤.雷达成像技术[M].北京:电子工业出版社, 2008. BAO Zheng, XING Mengdao, WANG Tong. Rada imaging technology[M]. Beijing:Publishing House of Electronics Industry, 2008.
[15] LI Ning, NIU Shilin, GUO Zhengwei, et al.Raw data-based motion compensation for high-resolution sliding spotlight synthetic aperture radar[J]. Sensors, 2018, 18(3):842. doi: 10.3390/s18030842
[16] BRANCATO V, JAGER M, SCHEIBER R, et al. A motion compensation strategy for airborne repeat-pass SAR data[J]. IEEE Geoscience and Remote Sensing Letters, 2018(15): 1580-1584.
-
期刊类型引用(1)
1. 戴宗武，张少甫，刘乃强，胡月，王喆，段岑薇，孙宇哲. 用于自主导航的SAR/多普勒雷达一体化技术设计. 兵工学报. 2022(12): 3093-3102 . 
百度学术
其他类型引用(1)
计量
- 文章访问数: 614
- HTML全文浏览量: 201
- PDF下载量: 29
- 被引次数: 2
陕公网安备 61011302001501号 