高分辨率星载SAR起伏运动目标精细聚焦与参数估计方法

温雪娇仇晓兰尤红建卢晓军

①(中国科学院空间信息处理与应用系统技术重点实验室 北京 100190)

②(中国科学院电子学研究所 北京 100190)

③(中国科学院大学 北京 100049)

④(中国国际工程咨询公司 北京 100048)

高分辨率星载SAR起伏运动目标精细聚焦与参数估计方法

温雪娇*①②③仇晓兰①②尤红建①②卢晓军④

①(中国科学院空间信息处理与应用系统技术重点实验室 北京 100190)

②(中国科学院电子学研究所 北京 100190)

③(中国科学院大学 北京 100049)

④(中国国际工程咨询公司 北京 100048)

高分辨率星载SAR图像中目标的复杂运动将引起不可忽视的散焦现象，影响目标识别和解译。该文对高分辨率星载SAR中目标起伏运动引入的误差进行了定量化分析，并仿真了其在SAR图像中的成像特点，提出了一种基于星载复图像数据进行运动误差补偿和精细聚焦处理的方法，同时估计得到目标的运动参数和海浪信息。仿真实验和TerraSAR-X港口区域实际数据实验验证了该方法的有效性和正确性。

星载SAR；高分辨率；运动目标精细聚焦；运动参数估计

1 引言

高分辨率合成孔径雷达(Synthetic Aperture Radar,SAR)可以获得丰富的地物或目标细节，在军事和民用各方面有着非常高的应用价值。然而高分辨率带来了更长的合成孔径时间，这使得目标运动对成像的影响更加显著。尤其在海面或港口，受到海面风浪的影响，舰船等目标将存在复杂的运动，使得舰船自身的散射以及舰船与海面的多次散射均存在合成孔径时间内的变化特性，导致SAR图像存在严重的散焦现象，给目标检测、识别和图像解译带来困难。因此高分辨率SAR中运动目标的精细聚焦是当前SAR应用迫切需要解决的问题之一。

国内外有关SAR运动目标的研究着重于地面运动目标指示(GMTI)[1,2]、运动目标成像算法、运动目标多普勒参数估计等方面[3,4]，研究大部分基于原始回波开展，以匀速直线运动、匀变速直线运动为主[5,6]。其中关于海面舰船的研究也主要是以平动为主、基于机载SAR模型开展运动误差的分析[7]。涉及舰船复杂运动的大多数研究则侧重于依靠舰船自身运动(如俯仰、横滚、偏航等3维运动)使用逆合成孔径雷达(ISAR)成像原理实现舰船目标的高分辨率成像[8,9]。目前非常缺少从高分辨率星载SAR图像产品出发，对散焦的海面复杂运动目标进行精细化分析处理的研究，也缺少对复杂运动信息进行估计与提取的研究。

鉴于此，本文主要针对高分辨率星载SAR中目标起伏运动引入的误差进行了定量化分析，并仿真了其在SAR图像中的成像特点，提出了一种基于复图像数据进行运动误差补偿和精细聚焦处理的方法，估计得到目标的运动参数和海浪信息，仿真和实际数据的处理结果验证了本文分析的正确性和方法的有效性。本文结构安排如下：第2节建立星载SAR目标起伏运动下的几何模型，定量化分析起伏运动与海面二次散射引入的误差与影响；第3节提出精细聚焦处理方法和运动参数估计方法；第4节给出了仿真实验和TerraSAR-X港口区域实际数据实验处理结果；第5节总结全文。

2 星载SAR起伏运动目标信号模型与误差分析

2.1 起伏运动目标单次散射误差分析

图1为星载起伏运动点目标成像几何。T和P分别为某一方位时刻η目标静止和起伏运动的位置，S为该时刻的卫星位置，设起伏运动沿地球矢径方向，幅度为A，频率为ω，初相为φ，点P相对于点T的距离矢量为：

其中A为目标沿TP方向振动的幅度矢量。

点P的距离历程为：

其中γ为目标点T处的波束入射角，计算公式为：

其中β为卫星下视角，a为卫星与目标之间的地心角。

将式(2)作近似处理：

计算得到目标起伏运动引入的距离历程误差为：

根据SAR回波与成像理论，进一步计算得到相位误差为：

对式(6)做Jacobi-Angle展开，得到：

可见式(7)中的指数项将导致多重对称的多普勒频移，产生成对回波，成像结果在方位向上将出现原信号的多重对称复制，进而引起散焦现象，其散焦程度与起伏运动的幅度和周期有关。

2.2 海面起伏二次散射误差分析

部分海面目标具有垂直于海面的结构，该结构会引起目标与海面间的多次散射。若海面静止，该结构与海面形成一个理想的二面角结构，但由于海风与地理结构的关系，海面存在着复杂的起伏运动，使二次散射的距离历程发生变化，进而影响成像结果。如图2所示，目标T为垂直结构上的任一散射点，某方位时刻η天线波束以入射角经过目标点T，并经由海面二次散射返回雷达接收机，图中实线为静止海面，虚线表示海面的起伏状态，设目标T距离静止海面的高度为h0，该时刻与雷达之间的距离为R(η)，海面起伏运动的幅度为A，周期为ω，初相为φ。目标T的往返距离历程为：

当海面静止时，目标T的往返距离历程为：

计算得到海面起伏运动引入的二次散射距离历程误差为：

式(10)与式(5)具有相同的形式，其区别主要在于幅度、频率和初相的不同，但对成像结果产生的影响具有相同的规律，在此不做赘述。实际上，本文考虑的二次散射是类似二面角的二次散射，其距离历程等效为二面角角点的单次散射，因此本文后续的精细聚焦及参数估计算法对单次散射模型和二次散射模型是通用的。需要说明的是，由于二次散射通常比单次散射具有更大的强度，运动对二次散射的影响往往在图像上表现的更为突出。

3 精细聚焦算法与参数估计

3.1 基于复图像数据的精细聚焦算法

针对含有运动目标的复图像数据，本文提出的精细聚焦算法具体实现步骤如图3所示。

(1) 方位向反压缩

通过方位匹配滤波器逆乘将数据变换至多普勒域，设多普勒频率为fη，距离向时间为t，后向散射系数为σ，雷达波长为λ，发射信号带宽为B，则得到方位向未压缩的数据：

(2) 距离包络对准

对方位未压缩数据应用最小熵距离对准算法[10]，实现目标的距离向包络对齐，在不需要提供先验知识的情况下，消除运动产生的相邻回波在距离向上的错位。这是一个最优化问题，目标函数为所有1维距离像和的熵值，可以表示为：

其中r为距离向采样点代表第n个回波的距离偏移量。

pave(r)的锐化度可以用来衡量回波包络对齐的程度，锐化度达到最高时，回波基本对齐。而香农熵可以用来衡量的锐化程度，波形锐化度越高，熵值越小，因此最优化准则为的熵值最小，通过迭代可以得到全局最优解。理论上此方法的对准精度为一个距离分辨单元，为了适应高分辨率的要求，本文实验均对输入数据在距离向做了32倍插值，使实际的对准精度达到1/8～1/4个距离分辨单元。

(3) 相位误差校正

将距离对准后的数据变换至多普勒域，由于各次回波间的相位误差通常可以和雷达波长相比拟，用距离偏移量来校正相位不精确，因此本文应用自聚焦算法实现误差相位校正。相比于基于特显点的自聚焦算法(Phase Gradient Autofocus,PGA等)，基于图像整体信息的自聚焦算法在低信噪比和低图像对比度的条件下具有更好的鲁棒性，能够取得更好的聚焦效果。本文采用全局最小熵的相位校正算法(Minimum Entropy Autofocus,MEA)[11]，该算法是一个求解补偿相位使图像聚焦效果最好的最优化问题，目标函数为一幅图像的熵值：

其中n为方位向采样点为经过相位补偿的图像，可以表示为：

图像的聚焦程度可以用图像的熵来衡量，聚焦程度越高，图像的熵就越小，图像也越清晰。因此最优化准则为式(14)的值达到最小，在给定补偿相位初始值的情况下，利用数值迭代算法逐渐逼近最终求解目标函数，得到的全局最优解。本文实验在输入数据序列中选取了强度较大的样本距离单元，并根据不同距离单元的相位方差对各样本距离单元进行加权，使含有多特显点的“高质量”样本在熵值收敛中起主导作用，进一步降低了杂波以及噪声等对算法的影响，同时提高了收敛速度，另外本文实验以变步长的梯度下降法寻找Δφ(m)的值[12]，进一步提高了迭代的计算速度，同时保证了最终的图像质量。

经过距离对准与相位校正后的数据可以表示为：

(4) 利用方位向FFT将数据变回图像域，同时实现方位向压缩，得到聚焦后的结果。

3.2 参数估计方法

根据SAR成像理论，方位时间域与多普勒频域存在着一一对应的关系，因此斜距误差的表达式在方位时域与方位多普勒域具有近似相同的表达形式，将替换成第2节中分析得到的斜距误差表达式，可以得到：

4 仿真分析与实验

为了验证上述分析的正确性和校正方法的有效性，下文首先对低轨高分辨率SAR进行仿真实验，接着将该方法用于TerraSAR-X的实际数据处理，并给出了实测结果。

4.1 仿真实验

仿真系统参数如表1所示。根据近海岸海浪信息设置了3组起伏运动参数，分别为A=0.1 m,Tv=2.5 s;A=0.1 m,Tv=4 s和A=0.05 m,Tv= 2.5 s。其中Tv=2π/ω。起伏运动点目标单次散射成像结果如图4所示，海面起伏时点目标二次散射成像结果如图5所示。选取运动参数为A=0.1 m,Tv=4 s的两组成像结果做32倍插值，得到的方位向波形如图6所示，应用本文精细聚焦方法处理，两组结果的聚焦结果及对应方位向波形图分别如图7和图8所示，同时得到的相位补偿曲线如图9所示，参数估计结果及估计误差如表2所示。

表1 雷达系统仿真参数Tab. 1 Radar system simulation parameters

表2 参数估计结果Tab. 2 Estimating results

由上述仿真结果可见，起伏运动会使成像结果在方位向出现原点目标的多重复制，进而形成线状散焦，散焦程度随起伏周期变大而减小，随起伏幅度变小而减小，这一结果与式(6)和式(12)的理论分析一致。两组散焦结果经过精细聚焦后的方位波形较理想，质量评估达标，从而验证了本文精细聚焦方法的正确性和有效性。两组参数估计结果都比较准确，估计误差较小，证明了本文参数估计方法的正确性和有效性。

4.2 实际数据处理结果

本文选取了TerraSAR-X旅顺地区1景SAR图像进行实验处理，该景图的相关参数如表3所示。

表3 TerraSAR-X旅顺岛参数表Tab. 3 Parameters for TerraSAR-X

第1幅海岛近海岸SAR图像与对应的光学图像如图10(a)与图10(c)所示，可见海岛边呈须状散焦的亮线是拴在岸边的浮桶，由于海浪的作用存在起伏运动，其散焦特征与仿真结果相似。选取图10(a)中左下角右起第2条亮线做时频分析，其结果如图11(a)所示，可见散焦亮线方位向时频存在正弦调制的关系，对其应用本文的精细聚焦算法，其聚焦前后的对比结果如图11(b)所示，可见精细聚焦效果明显，补偿相位如图11(c)所示，其参数估计结果为：

第2幅近海岸码头SAR图像如图12(a)所示，可见码头停靠的潜艇尾部有较整体更明亮的线状散焦，码头停靠的潜艇光学图如图12(c)所示，可见潜艇尾部的垂直翼与海面形成了1个二面角，潜艇尾部的二次散射较多，对该幅图像做精细聚焦处理的结果如图12(b)所示，可见聚焦效果明显，其起伏运动参数估计结果为：

考虑到码头海浪的实际情况，两幅图像估计得到的参数都较为合理，证明了本文方法的准确性与有效性。

5 总结

本文对高分辨率星载SAR中目标起伏运动引入的误差进行了定量化分析，并仿真了其在SAR图像中的成像特点，提出了一种基于星载复图像数据进行运动误差补偿和精细聚焦处理的方法，同时估计得到目标的运动参数和海浪信息。仿真实验和TerraSAR-X港口区域实际数据实验验证了该方法的有效性和正确性。

[1]Li X,Deng B,Qin Y,et al.. The influence of target micromotion on SAR and GMTI[J].IEEE Transactions on Geoscience and Remote Sensing,2011,49(7): 2738–2751.

[2]Wei S and Wang H. The improvement of the conventional GMTI method with single-channel SAR[C]. 2004 IEEE Geoscience and Remote Sensing Symposium,Anchorage,Alaska,USA,Sept 2004: 2626–2628.

[3]Bamler R. Doppler frequency estimation and the Cramer-Rao bound[J].IEEE Transactions on Geoscience and Remote Sensing,1991,29(3): 385–390.

[4]Li G,Xia X,Xu J,et al.. A velocity estimation algorithm of moving targets using single antenna SAR[J].IEEE Transactions on Aerospace and Electronic Systems,2009,45(3): 1052–1062.

[5]Martorella M,Berizzi F,Pastina D,et al.. Spaceborne radar imaging of maritime moving targets with the Cosmo-SkyMed SAR system[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2014,7(7): 2797–2810.

[6]Noviello C,Fornara G,and Martorella M. Focused SAR image formation of moving targets based on Doppler parameter estimation[J].IEEE Transactions on Geoscience and Remote Sensing,2015,53(6): 3460–3470.

[7]Thayaparan T,Abrol S,Riseborough E,et al.. Analysis of radar micro-Doppler signatures from experimental helicopter and human data[J].IET Radar,Sonar and Navigation,2007,1(4): 289–299.

[8]Chen V C and Lipps R. ISAR imaging of small craft with roll,pitch and yaw analysis[C]. 2000 IEEE International Radar Conference,Alexandria,VA,USA,May 2000: 493–498.

[9]Noviello C,Fornaro G,Martorella M,et al.. ISAR add-on for focusing moving targets in very high resolution spaceborne SAR data[C]. 2014 IEEE Geoscience and Remote Sensing Symposium,Quebec,Canada,July 2014: 926–929.

[10]Zhu D Y,Wang L,Tao Q N,et al.. ISAR range alignment by minimizing the entropy of the average range profile[C]. 2006 IEEE Conference on Radar,Verona,New York,USA,April 2006: 813–818.

[11]Li X,Liu G,and Ni J. Autofocusing of ISAR image based on entropy minimization[J].IEEE Transactions on Aerospace and Electronic Systems,1999,35(4): 1240–1252.

[12]Wang J,Liu X,and Zhou Z. Minimum entropy phase adjustment for ISAR[J].IEE Proceedings-Radar,Sonar and Navigation,2004,151(4): 203–209.

温雪娇(1992–)，女，吉林长春人，中国科学院电子学研究所在读硕士研究生，主要研究领域为高分辨率星载SAR成像技术、复杂运动目标精细处理技术。

E-mail: 136037701@qq.com

仇晓兰(1982–)，江苏苏州人，女，中国科学院电子学研究所副研究员，研究方向为SAR成像技术、双基地SAR技术。

E-mail: xlqiu@mail.ie.ac.cn

尤红建(1969–)，江苏南通人，男，中国科学院电子学研究所研究员，研究方向为遥感图像处理技术。

E-mail: hjyou@mail.ie.ac.cn

卢晓军，江苏泰州人，北京理工大学博士后，中国国际工程咨询公司高级工程师，专业方向为智能控制、信号处理。

E-mail: lu8new@163.com

Focusing and Parameter Estimation of Fluctuating Targets in High Resolution Spaceborne SAR

Wen Xuejiao①②③Qiu Xiaolan①②You Hongjian①②Lu Xiaojun④
①(Key Laboratory of Geo-spatial Information Processing and Application System Technology,Beijing100190,China)
②(Institute of Electronics,Chinese Academy of Sciences,Beijing100190,China)
③(University of Chinese Academy of Sciences,Beijing100049,China)
④(China International Engineering Consulting Corporation,Beijing100048,China)

Complex motion can cause serious defocusing phenomena in high resolution spaceborne SAR cases,which then lead to decreased image resolution. In this study,we built a simulation model to quantitatively analyze the signature and effect on maritime fluctuating targets in high resolution cases. To deal with formed Single-Look Complex (SLC) SAR images containing fluctuating targets,we implement a motion-compensation and fine-focusing method to obtain refocused images and the fluctuation parameters. We demonstrate the effectiveness and correctness of the proposed approach in focusing and estimating the parameters of fluctuating targets by processing the simulation results and archived images acquired by Terra-SAR in hybrid spotlight mode.

Spaceborne SAR; High resolution; Moving target fine focusing; Parameter estimation

TN957.52

A

2095-283X(2017)02-0213-08

10.12000/JR17005

温雪娇,仇晓兰,尤红建,等. 高分辨率星载SAR起伏运动目标精细聚焦与参数估计方法[J]. 雷达学报,2017,6(2): 213–220.

10.12000/JR17005.

Reference format:Wen Xuejiao,Qiu Xiaolan,You Hongjian,et al.. Focusing and parameter estimation of fluctuating targets in high resolution spaceborne SAR[J].Journal of Radars,2017,6(2): 213–220. DOI: 10.12000/JR17005.

2017-01-13；改回日期：2017-02-20；

2017-04-25

*通信作者： 温雪娇 136037701@qq.com

国家自然科学基金(61331017)，国家高分重大专项(30-Y20A12-9004-15/16,41-Y20A13-9001-15/16)

Foundation Items: The National Natural Science Foundation of China (61331017),The Key Standard Technologies of National High Resolution Special (30-Y20A12-9004-15/16,41-Y20A13-9001-15/16)