平流层飞艇SAR成像优化策略研究.DOC

1、平流层飞艇 SAR 成像优化策略研究贾高伟 杨希祥（国防科技大学空天科学学院 湖南长沙 410073）摘要：平流层飞艇具有驻空高度高、覆盖范围广、驻空时间久、使用效费比高等特点，在侦察监视等领域应用前景广阔。结合平流层飞艇平台特性，分析了平流层飞艇合成孔径雷达（Synthetic Aperture Radar, SAR）与机载 SAR 的异同；针对平流层飞艇 SAR 的 PRF 冗余及大振幅低频率运动误差，提出了优化方案，并开展了仿真分析验证。仿真结果表明本文提出的优化运动补偿方法能够有效消除平流层飞艇运动误差影响。关键词：平流层飞艇；快速响应；SAR 系统中图分类号：TN95 文献标志码：A

2、 文章编号： Imaging Optimization for Stratospheric Airship SARJIA Gaowei, YANG Xixiang(College of Aerospace Science and Engineering, National University of Defense Technology, Changsha, Hunan, 410073)Abstract：As a new SAR carrier, stratospheric airship has the potential as large observation area, long dr

3、ift time, and better efficiency. It is expected in many different applications, such as reconnaissance, surveillance, and environmental monitoring. In the paper, the imaging difference between the stratospheric airship SAR and airplane SAR is compared. To improve the imaging efficiency, the optimiza

4、tion strategy according to redundant PRF and large amplitude trajectory deviation, respectively. Eventually, necessary simulation test is presented to validate the proposed method.Keywords: Stratospheric Airship; Fast response; Synthetic aperture radar 目前，我国 SAR 系统搭载平台主要有天基、空基和地基三种：近地轨道遥感卫星轨道高度在 500

5、km800km 之间，能够对固定目标或移动速度较慢目标实施高分辨侦察，但星载 SAR系统开发成本高，代价大，轨道运行特性使其难以对突发关键区域进行快速观测响应。相比于卫星，飞机的飞行高度低、机动灵活，但滞空时间短，机载 SAR 往往需要高额费用和高频率出勤率才能实现大范围、长时间的预警、侦察和监视。地基 ISAR 对空侦察的范围受到观测点地域限制，监视范围小，不能适应在广阔海域中的应用需求。临近空间浮空器（平流层飞艇或气球）技术的发展为实现大范围区域预警和侦察监视提供了一种新的装备平台1。平流层飞艇的飞行高度适中（约 20Km） 、飞行速度较慢、滞空时间长，将临近空间飞艇配置在国境线、海岸线、

6、争议区域等并配置有效 SAR 载荷，可以对海面、陆地上目标进行有效观测和侦察，填补临近空间观测活动利用的空白，在军事、遥感、勘探、测绘、农业、水利等方面同样具有十分广泛的应用前景2。具体地，影响到平流层飞艇实际运行及执行任务能力的关键因素包括：临近空间低压低温环境、平流层飞艇蒙皮材料、平流层飞艇能源管理系统、平流层飞艇飞行控制等。这些因素将影响平流层飞艇 SAR 系统的温控设计、SAR 体积重量功耗、SAR 成像算法及运动补偿算法等3-5 。本文将着重分析平流层飞艇 SAR在成像几何、成像算法、运动补偿处理等方面的特性。_*作者简介：贾高伟（1989-） ，男，河南周口人，讲师，博士，E-ma

7、il： ji_杨希祥，男，河北阜城人，副教授，博士。1. 平流层飞艇 SAR 与机载 SAR 的差异本文分析的平流层飞艇的运行高度约18km-20km,该环境条件下最低大气温度约-70 ，大型机载 SAR 平台飞行高度约 8-10Km, 对应的大气温度约-30，因此平流层飞艇 SAR 系统设计当考虑器件的低温特性或载荷系统的温控系统。同时，为减小 SAR 系统对平流层飞艇能源系统的压力，应当在 SAR 系统发射功率与天线增益之间做出折衷。此外，平流层飞艇体积庞大，适合安装轻质大面阵高增益天线系统。在 SAR 工作环境特性、功耗尺寸限制之外，平流层飞艇 SAR 与机载 SAR 的最大差异体现在巡

8、航速度上。一般机载平台（如运八）的巡航速度在 150m/s，而平流层飞艇的巡航速度在10-15m/s 之间1,6-7。这里设定平流层飞艇的巡航速度为 10m/s。低速运动特性对平流层飞艇的影响包括：合成时间延长，航向运动误差补偿的必要性及 PRF 冗余。1.1 成像几何的差异类似于机载 SAR，平流层飞艇 SAR 的成像几何为：图 1 平流层飞艇 SAR 成像几何Fig.1 The geometry of stratospheric airship SAR设 H 为平流层飞艇高度， 为合成孔径BW积累角， 为零多普勒斜距， 为雷达相位中0RR心（APC）同目标 之间的瞬时斜距，n(,0)Txy

9、设平流层飞艇巡航速度为 。结合 SAR 成像知V识，中心斜距处平流层飞艇 SAR 的合成孔径长度为：(1)BWSAR02tan()L设平流层飞艇 SAR 与机载 SAR 的典型参数如下：表 1 参数对比Tab.1 Parameters comparison参 数 机载 SAR 平流层飞艇 SAR运行高度 8000m 20000m运行速度 150m/s 10m/s雷达载频 10GHz 10GHz信号带宽 400MHz 400MHz方位向波束宽度 2.3 2.3雷达下视角 4575 4575结合表 1 可以计算得到中心斜距处机载SAR 的合成孔径长度为 567m,对应合成孔径时间为 3.78s；而

10、中心斜距处平流层飞艇 SAR 的合成孔径长度为 1417m，对应的合成孔径时间约为 94.5s。以观测静止场景或目标为例，理论上合成孔径时间的延长对成像质量并无影响，但无论是机载 SAR 还是平流层飞艇 SAR，都无法避免气流扰动及控制偏差带来的轨迹偏离，这对实际成像处理带来了挑战。1.2 运动误差的影响平流层飞艇的总体布局，飞行机理和工作模式不同于导弹、飞机等飞行器1, 8。平流层飞艇主要依靠浮力气体提供静升力，具有很大的体积/质量比，惯性大等特点，附加质量和附加惯性作用明显。同时，平流层飞艇为充气柔性体，存在一定程度的弹性变形，外界大气环境对飞艇有着复杂的多物理场耦合作用，因此其运动误差特

11、性不同于常规飞机。针对 SAR成像关注的航迹和姿态，其特点概括为：（1） 体积/质量比大，飞行速度缓慢，惯性特征显著，运动误差的频率响应主要是低频。（2） 平流层飞艇体积大，在临近空间环境中（低动压）气动控制效率低，控制响应迟缓，导致其偏离预设航迹的幅度大。由此得知，平流层飞艇将带来大振幅低频率运动误差。一般地，方位向多普勒调频率可以定义为 ， 为航向理想速度。对20/vR于航向速度误差 ，引入的多普勒调频率相对v比值为 。由于平流层飞艇巡航速度2()/慢，因而相同的航向速度偏差对应的多普勒调频率相对变化量很大，对方位向成像影响显著，需要进行 PRF 实时调整予以克服。1.3 PRF 冗余的影

12、响以正侧视条带 SAR 为例，SAR 的多普勒带宽为：(2)BWa4sin()2VSAR 系统设计时要求 ，由表 1 可aPRF知，在同等方位分辨率条件下，机载 SAR 和平流层飞艇 SAR 对应的多普勒带宽分别为 300Hz和 20Hz，有 15 倍冗余。PRF 冗余将增加方位向处理数据量，但通过方位滤波或抽样等处理带来的效益并不高。2. 相应的优化策略第二章的分析表明，平流层飞艇 SAR 相比于机载 SAR，具有合成时间长、孔径长、PRF冗余大、运动误差幅度较大等特点。从成像机理上讲，这些因素并不会直接影响 SAR 成像质量。但为了提高成像效率、更为高效地利用平流层飞艇平台开展对地高分辨率

13、成像，有必要对平流层飞艇 SAR 成像过程进行优化，具体包括冗余 PRF 的利用以及改进适用于频域成像算法的运动补偿流程。2.1 冗余 PRF 的利用结合平流层飞艇运动速度慢的特点，利用阵列数字波束形成技术，可以分时收发数据，拓展平流层飞艇成像（观测）范围。具有来讲，包括方位向多波束切换和距离向多波束切换两类，结合具体雷达参数，也可以将两者有效结合起来。(1) 方位向多波束切换在不同的脉冲重复间隔（PRI） ，通过数字波束生成不同斜视角的波束，以三种波束角为例，在三个时刻，分别以后斜视角 ，正侧视，s前斜视角 辐射电磁波。随后对接收到的回波s区分，分别生成后斜视回波矩阵，正侧视回波矩阵，前斜视

14、回波矩阵。同样的合成孔径长度内，可以得到三幅雷达图像，通过拼接、裁剪重叠区域，能够获取更长的方位向测绘宽度。如图 2 所示，PRI 1，PRI 2，PRI 3 三个时刻对应的飞艇空间位置近乎相同，图中显示是为区分示意三个不同时刻。若按照正侧视录取回波，对应的方位向测绘宽度观 测 场 景图 2 方位向多波束切换示意图Fig.2 Multi-beam switch in azimuth direction为 ，采用方位向波束切换后，对应的测绘宽0L度为 。通过参数设计，该处理模式能够允sum许 SAR 在成像过程中间歇性关机，节约电量。(2) 距离向多波束切换另一种利用冗余 PRF 的方式是距离向

15、多波束形成，同样地以三个波束为例，在PRI1，PRI 2，PRI 3 三个时刻，生成下视角为 ，1， 的三个波束并辐射电磁波。随后对接收23的回波区分，成像得到对应于不同测绘带的雷达图像。各测绘带之间的重叠区域用于图像拼接。在这一成像模式下，需考虑到距离模糊问题，并选取合理的雷达回波接收增益。距离向多波束切换，能够很大程度地增大测绘带宽。 高 度 向地 距 向方 位 向图 3 距离向多波束切换示意图Fig.3 Multi-beam switch in range direction需要说明的是，无论是方位向还是距离向多波束切换，机械扫描因重量大，响应时间较长等原因，都不适合于平流层飞艇。可以通

16、过多通道多天线设置或者相扫来实现多波束切换。2.2 大空变运动误差的影响与消除（1）大幅度运动误差的影响测 绘 带图 4 实际运行轨迹下平流层飞艇 SAR 成像几何Fig.4 The geometry of stratospheric airship SAR in the presence of motion errors图4为非理想情况下正侧视平流层飞艇SAR的成像几何示意图。在笛卡尔直角坐标系XYZ 中，O为坐标原点，预定航迹高度为 H，预定沿X 向飞行，图中沿X向虚线代表实际飞行轨迹，为实际飞行轨迹上的雷达相位中心(APC)位P置。相对于理想飞行，设雷达APC沿X、Y 、Z向误差分别为

17、， 和 ， 是m()t()Ytm()t方位慢时间。可以得到 点的坐标为P。 表示APC 沿X向m(),(,)tYtZt()Xt的实际位置，对于测绘带中任意一点，其与雷达APC的瞬时斜距为：n,0Txy222mnmnm()()()()RtXtxYtyZtH（3）定义目标 到理想航迹的零多普勒斜距为：n（4）20nRHy理想情况下，目标 到雷达 APC 之间的瞬T时斜距为：（5）22rmn()RtVtxyH则非理想运动轨迹引入的瞬时斜距误差为：m222nmnm2m()()()()rRttXxYtyZtHVtyH（6）是运动补偿的核心要素，在运动补()Rt偿处理中， 的表达形式及方位空变性是mt重点

18、需要考虑的。结合成像知识，SAR 平台运动误差可以分为沿视线方向（Line of Sight, LOS）和沿方位向。如前文所述，这里假设方位向速度误差已由实时调整 PRF 克服。为此，这里忽略沿方位向（X 向）速度误差，仅考虑LOS 向运动误差。经方位向 PRF 调整后，可以等效认为 SAR平台匀速运动，进而得到瞬时运动误差的表达形式为 222rmnmnm()()+()RtVtxYtyZtHH（7）针对常规机载 SAR 成像几何， 的rm()Rt近似表达式为9 ：（8 ）0rmmr m()()sin()costcosRtYtZt基于式（8），取 得到视线平面内运0动误差为（9）Losmm()(

19、)sin()cosRtYtZt针对常规机载 SAR 成像几何，利用式(9)表示雷达 APC 在视线平面内运动误差是满足要求的。式（9）的成立前提是 、 远小m()Yt()Zt于 ，而这样的假设对于平流层飞艇 SARrm()Rt而言不再适合。相应地，对瞬时斜距误差进行如下处理：rm222nmnm220 mrm r0m3rm(,)()()() ()()(si()cos2in)28()cos()sicosRtVxYtyZtHtHYtZtttRt RYRZtttt22m030 m30()()2in()()cos()8YtZtRYZtttR(10)令 取 （即 ） ，得到零多普勒rm()t0位置对应的运

20、动误差，亦即 LOS 向运动误差为：Losm22mm02200 m30(,)()()in()cos2s28Rt YtZtYZtRRttt（11）理 想 位 置实 际 位 置图 5 视线平面运动误差示意图Fig.5 the motion errors in the LOS plane 图 5 给出了 、 、 等参数的Losm(,)Rt0几何示意。结合式（10） ，对于某一时刻 对应的运动mt误差 和 ，波束边缘处对应的误差m()Yt()Zt为：am 22BWmmm020223BW30m() ()()cos()sin()cos(ic)()cs()8Rt YtZtYtZtRttR （12）则对应的运

21、动误差最大方位空变量为：（13）aLosma()()rtt量化地，对于平流层飞艇 SAR，设和 满足余弦曲线，振幅均为mYt()Zt100m。结合表 1 中雷达参数及成像几何，随着下视角的增大，计算得到如图 6(a)所示曲线。观察发现，式(9) 引入的近似误差在距离分辨单元的量级，会严重降低成像质量。图 6(b)表示利用式（11）计算对应的近似误差，可以看到在零多普勒位置，近似引入的误差小于波长，满足成像要求。此时的运动误差对于波束边缘点将引入空变误差，如图 6(c)所示，可以看到由于运动误差幅度较大，对于波束边缘点仍存在不容忽略的空变误差，但量(a) 传统公式近似引入的误差(a) The e

22、rrors derived from the traditional formula(b) 公式(11) 对应的近似误差(b) The errors derived from equation (11)(c)公式(11)对应的误差最大空变量(c) the maximum of space variable by equation (11)图 6 运动误差量化分析Fig.6 Numerical analysis of motion error值在小于一个距离分辨单元，表明可以在一阶补偿中忽略式(11)引入的方位空变量。但在二阶补偿中，必须考虑方位空变性。运动误差的补偿，涉及三个方面，一是运动误差

23、的精确表达形式，二是运动误差的方位空变性，三是运动误差的距离空变性。对于距离空变性，可以通过距离向分子带成像处理，并以距离向子测绘带中心斜距对应的误差完成一阶运动补偿，本文对此不再展开阐述。后续的运动误差补偿着重考虑前两因素的影响。（2）基于频域成像处理的运动误差消除事实上，无论是讨论运动误差的解析表达式还是考虑运动误差的空变性，都是为了配合频域成像算法，如 NCS 算法， K算法，因为频域成像算法计算量小，对数据存储空间要求也小，利于实时成像处理。在频域成像处理中，“两步运动补偿法”是可行且高效的运动补偿方法。以 NCS 算法为例，基于式(11) ，在距离频域方位时域，进行一阶运动补偿，此时

24、的补偿因子为：（14）moc1Losmrefexp(j2(,)HkRt表示距离向波数，随后进行变标处理、距离k徙动校正、回波数据变换至距离压缩，方位时域。由于运动误差的方位空变性，需要利用子孔径处理的方式克服空变性。在每一段子孔径内，近似认为运动误差是不变的。利用方位频率同雷达成像角 之间的对应关系，将子块回波变换到方位频域，2sin/afV则不同方位频率对应于不同的雷达成像角，利用这种关系可以进行与斜视角相关的方位空变相位误差补偿，具体补偿因子为：（moc2crmidLosmidrefexp(j(,)(,)HkRtt15）其中 表示第 个子孔径的方位idt,0iI中心时刻， 为总的子孔径数量

25、。I为第 个子孔径对应的沿不同下视角rmid(,)Rti的平均运动误差， 为第 个子孔Losmidref(,)Rti径对应的沿中心斜距处 LOS 向运动误差。 为ck中心波数。完成二阶补偿后，将子孔径数据变换至方位时域，重新合并，再进行方位脉冲压缩处理。如此即完成了空变运动误差补偿。3. 仿真分析对于第 2.2 节介绍的改进运动补偿方法，本章进行仿真验证，采用表 1 所示平流层飞艇SAR 1 2 34 5 67 8 9(a) 目标布置示意图 (a) displacement of multi targets(b) 采用经典运动补偿方法对应的成像结果(b) The results process

26、ed with traditional MOCO method(c) 未考虑方位空变性的成像结果(c) The results processed without considering the azimuth space variance(d)完整的运动补偿成像结果(d) The results processed with the proposed MOCO method图 7 仿真结果对比Fig.7 Results comparison参数，并设运动误差 和 为振幅 100 米的YZ余弦曲线，距离向测绘子带宽度设为 150 米。在场景中设置 9 个点目标，方位向间隔 20 米，地距向间隔

27、 30 米。结合成像处理，可以得到如图 7 所示仿真结果。直观地，利用本文方法处理能够有效消除运动误差的影响，得到的良好的成像结果。以边缘处点目标 3 为例，三种不同的运动补偿处理对应的距离向/方位向成像结果剖面图如图 8 所示，其中点线、短划线、实线分别对应传统“两步法”运动补偿（MOCO） 、未考虑方位空变、本文所提方法处理得到的成像结果。对比发现，针对平流层飞艇 SAR 运动误差幅值大、周期长的特点，传统的运动误差近似处理已经不能满足应用需要，同时，运动误差的方位空变性必须予以考虑。本文提出的运动误差近似处理及方位空变性处理方法，能够克服运动误差的影响，得到聚焦良好的 SAR 图像。(a

28、) 距离向成像结果剖面图(a) Profiles in range(b) 方位向成像结果剖面图(b) Profiles in azimuth图 8 方位向成像结果对比Fig.8 The comparison of results in azimuth量化地，图 8(b)中本文方法对应的方位向点散布函数的分辨率为 0.33m，同理论值一致；峰值旁瓣比（PSLR）约为-12.6dB ，积分旁瓣比（ISLR ）约-9.1dB，尽管 PSLR 和 ISLR 同精确成像的理论值有所降低，但满足实际应用，表明了 2.2 节运动补偿流程的可行性。未来工作将考虑本文误差补偿方法对图像几何失真的影响10。4.

29、结论本文建立了平流层飞艇 SAR 成像几何，结合平流层飞艇典型运行环境，对比分析了平流层飞艇 SAR 与机载 SAR 系统的异同，分析得到：相比于机载 SAR，平流层飞艇 SAR 并没有颠覆性的区别，适用于机载 SAR 的成像算法依旧适用于平流层飞艇 SAR，但对应的成像效率较低。但低速条件下航向运动误差对成像影响更为显著，需予以补偿；针对 PRF 冗余，本文提出了方位向波束切换和距离向波束切换策略提高平流层飞艇 SAR 成像效率；分析得到平流层飞艇存在大振幅低频率的轨迹偏移结论，并提出了改进的运动补偿策略，使之能够配合频域成像算法，仿真分析表明了所提改进方法的有效性。参考文献1 杨跃能，平流

30、层飞艇动力学建模与控制方法研究D，长沙：国防科技大学，2013.YangYueneng, Dynamics Modeling and Flight Control for a Stratopheric Airship D. Changsha: National University of Defense Technology, 2013.2 WenQin Wang, Near-Space Vehicle: Supply a Gap between Satellites and Airplanes for Remote sensing J, IEEE AES Magazine, April 20

31、11. pp.4-9.3 杨海光，临近空间 SAR 成像理论与成像方法研究D，成都：电子科技大学，2014.Yang Haiguang, Imaging Theories and Methods for Near-space Synthetic Aperture Radar D, Chengdu: University of Electronic Science and Technology of China, 2014.4 王晓峰，董臻，梁甸农，临近空间慢速平台 SAR结合运动补偿的 SCFT 算法J，系统工程与电子技术，vol.30, no. 10, 2008, pp.1869-1873.

32、Wang Xiaofeng, DongZhen, Liang Diannong, SCFT processing algorithm based on motion compensation for near space slow platform SARJ, vol.30, no. 10, 2008, pp.1869-1873. 5 左伟华，皮亦鸣，闵锐，临近空间平台高分辨率宽测绘带重复访问合成孔径雷达成像模式及算法研究J，信号处理，vol.30, no.10, 2014, pp.1119-1127.Zuo Weihua, Pi Yiming, Minrui, High Resolution

33、 Wide Swath Repeat Acquisition SAR Imaging Mode and Algorithm with Near Space Platform J, Journal of Signal Processing, vol.30, no.10, 2014, pp.1119-1127.6 杨海光，易青颖等，临近空间慢速平台 SAR 地面动目标检测与成像J，电子科技大学学报， vol.43, no.6, 2014, pp.838-844.Yang Haiguang, Yi Qingying, Near-Space Slow SAR Ground Moving Target

34、Detection and Imaging J, Journal of University of Electronic Science and Technology of China, vol.43, no.6, 2014, pp.838-844.7 Wenqin,Wang. Near-Space Vehicle Borne SAR with Reflector Antenna for High Resolution and Wide Swath Remote Sensing J, IEEE Trans. On GRS, vol. 50, no.2, 2012, pp. 338-348.8

35、Wen Qin Wang, Jingye Cai, Qicong Peng, Near-Space SAR: A Revolutionary Microwave Remote sensing MissionC, APSAR, 2007.9 保铮，邢孟道，王彤，雷达成像技术M, 北京：电子工业出版社，2005.Bao Zheng, Xing Mendao, Wang Tong, Technology of Radar Imaging M.Publishing House of Electronics Industry, 2005.10 李燕平,邢孟道, 保铮, 沿航向运动补偿的几何形变校正J,西安电子科技大学学报,vol.33,no.6, 2006,pp.881-886.Li Yanping, Xing Mendao, Bao Zheng, The Correction of Geometric Distortion for along-track motion compensation J, Journal of XiDian University, vol.33,no.6, 2006,pp.881-886.