朱 江邓佳欣廖桂生朱圣棋

①（西安电子科技大学雷达信号处理国家重点实验室 西安 710071）

②（中国飞行试验研究院 西安 710089）

基于稀疏采样的空间碎片群目标成像方法

朱 江*①邓佳欣②廖桂生①朱圣棋①

①（西安电子科技大学雷达信号处理国家重点实验室 西安 710071）

②（中国飞行试验研究院 西安 710089）

低重频、短数据条件下的雷达成像是比较困难的，同时，星载平台观测的碎片目标多以群目标形式出现。针对窄带条件下的空间碎片群目标成像问题，该文提出了一种基于稀疏重构的空间碎片群目标成像方法。由于空间碎片通常存在高速自旋现象，且碎片之间因质量、密度等物理差异存在明显的转速差，利用观测时间内碎片群的多周期观测数据与转速差异性，结合回波的自相关特性，能有效实现碎片群的转速检测。由于碎片在空域具有强稀疏特性，结合碎片转速，利用观测矩阵抽取数据，可以实现等效插值操作，并且将稀疏重构得到的数据进行重排，即可得到各个目标的像。抽取操作能够一定程度地抑制其他碎片的回波能量，并解决在低重频条件下的多普勒模糊问题。理论分析证明了该文所提方法的有效性，仿真实验表明，该文所提方法能够在低重频条件下实现对碎片群目标的分别成像。

自相关函数；多普勒模糊；稀疏采样；空间碎片群目标；ISAR成像

引用格式：朱江， 邓佳欣， 廖桂生， 等. 基于稀疏采样的空间碎片群目标成像方法［J］. 雷达学报， 2016， 5（1）： 82-89. DOI： 10.12000/JR16012.

Reference format： Zhu Jiang， Deng Jiaxin， Liao Guisheng， et al.. Space group debris imaging based on sparse sample［J］. Journal of Radars， 2016， 5（1）： 82-89. DOI： 10.12000/JR16012.

1 引言

随着空间资源开发日益增加，废弃的卫星碎片等垃圾在太空堆积，这给未来空间资源的利用带来安全隐患。如果在轨运行的飞船与碎片发生碰撞，即便尺寸小于1 cm的碎片，都能对飞船造成致命伤害。现有的雷达或者光学仪器能够实现对尺寸大于30 cm的目标监控，采用保护罩能够保护飞船免受尺寸小于1 cm的碎片伤害，因此开展对1～30 cm尺寸的空间碎片检测、成像与识别有广阔的应用前景和战略意义［1-3］。

空间碎片因为轨迹、惯性的筛选通常以群聚的形式出现，同时，碎片存在高速自旋运动，且各自的运动状态因质量、密度等因素存在明显差异。各个碎片的转动速度不同，导致各目标回波去相干时间不一致，无法对回波数据进行统一处理，为了得到清晰的各碎片ISAR像，有必要对各个目标的数据进行抽取。文献［4］利用各目标平动带来的多普勒调频率的差异，实现各目标回波的抽取及成像；文献［5］通过对粗成像的结果进行图像分割，抽取出属于各个目标的回波，再进行精确成像。然而上述方法均没有考虑高速自旋目标引起的距离徙动问题，且各目标的转速不同导致徙动量存在差异，在宽带雷达系统中，包络补偿难以实现，同时，对于小尺寸远距离的微弱目标，采用宽带雷达体制很难发现目标，有必要限制系统的带宽从而提高对微弱小目标的敏感性。此处采用窄带雷达系统，能够提高对于小目标的检测概率，同时，窄带系统的距离分辨率通常在几十米数量级，碎片自旋所引起的徙动问题可以忽略。由于空间碎片绕着其主轴做高速自旋运动，单距离多普勒干涉（Single Range Doppler Interferometry， SRDI）方法［6］能够对尺寸小于发射信号波长的目标进行成像，SRDI利用时频分析方法将目标所在距离单元的方位向回波数据转换到时频域，尽管线性时频方法不受交叉项的影响，但存在时间与频率分辨率的折衷，无法得到聚焦良好的图像。相干单距离多普勒干涉（Coherent Single Range Doppler Interferometry， CSRDI）方法［7］能够联合利用时频域的相位与幅度信息，在低信噪比条件下仍能够得到较好的性能，然而，CSRDI方法需要沿着不同的曲线进行积分，运算复杂度非常高。文献［8，9］提出了单距离匹配滤波方法及其修正算法（Single Range Matched-Filter， SRMF method and its modified version， SRMF-CLEAN）。SRMF根据不同的旋转半径构建匹配滤波器，并通过FFT估计各半径下的散射中心，然而，该方法存在明显的旁瓣，微弱目标很容易淹没在强信号的杂波中，不利于微弱目标的检测，同时，该方法没有考虑目标旋转引起的阴影现象，因此匹配滤波器无法匹配目标的回波信号，导致成像性能显著下降。上述方法分辨率主要取决于多普勒带宽，短波长引入更大的多普勒带宽，而在多任务情况下，ISAR成像系统倾向于发射低重复频率信号，此时，存在多普勒模糊现象，算法性能均会显著下降。由于空间碎片尺寸通常较小，在成像场景中表现出非常强的稀疏特性，可以利用该特性得到更好的性能。文献［10］提出了一种基于稀疏重构的空间碎片成像方法，该方法通过建立观测矩阵，解决了目标自旋情况下出现的遮挡效应，同时，抽取的数据提升了系统的等效采样率，解决了低重复频率系统的模糊问题。CS理论已经证明，在少量观测样本下，能够以高概率重构出稀疏信号［11］。然而，若字典构建不完备，信号能量渗透到其他的原子上，会导致恢复的信号出现错误。该方法在处理多目标成像问题时，效果显著下降。

针对窄带、低重复频率系统下的空间平台空间群目标成像问题，本文结合回波的自相关函数与转速扫描，获取观测范围内多目标的转速信息，从而构建完备的字典，并根据所获取的目标转速信息，对不同的目标设计不同的观测矩阵，提取各目标的有效信息，得到多周期观测下的目标有效数据，并结合SL0［12］方法进行稀疏恢复。由于该方法的抽取操作针对各个目标都是彼此独立的，因此可以进行并行操作，同时得到各个目标的成像结果。

2 窄带系统下的高速自旋目标回波模型

ISAR成像雷达通常发射大带宽信号获取高的距离向分辨率，然而，针对空间平台，雷达功率较小，无法有效检测到远距离的微弱目标，此时，有必要采用窄带系统，提高发现目标的概率。同时，窄带系统下距离分辨率通常在几十米到上百米量级，因此在成像时间内，高速自旋目标的旋转引起的跨距离徙动现象可以忽略，利于实现包络补偿。然而，窄带系统下，单个距离门内会存在多个目标。此时，考虑单个距离门的目标方位回波，第k个目标的第i个散射点回波信号如下：

其中σki表示第k个目标内第i个散射点的散射强度，fc代表发射信号的载频，tm表示方位慢时间，m表示脉冲数。τmki=2Rki（tm）/c， Rki（tm）是目标在成像时间内与雷达的斜距。假定已经实现了回波的包络对齐与相位补偿［13，14］，此时，目标的转台模型如图1所示。

图 1 成像几何模型Fig. 1 Geometry model of imaging

第k个目标的第i个散射点与雷达之间的斜距Rki（tm）表示为：

其中R0是雷达到目标群中心的初始距离，rki和φki分别是第k个目标第i个散射点的自旋半径与初始转角。包络补偿过后，雷达视线LOS（Line Of Sight）相对于目标存在一个转动量ωr，同时，目标的高速自旋也存在一个转动分量ωk。由于成像时间内目标转角较小，相比于自旋分量，可以忽略转动量ωr带来的影响。此时式（2）简化为：

将式（3）代入式（1），回波形式如下：

根据几何衍射（Geometrical Theory of Diffraction，GTD）理论［15］，在高频系统下的高分辨率图像，其目标的回波响应可以近似为多个散射点的响应之和，因此，该距离单元的方位基带回波可以表示为：

其中K代表碎片个数，Ik表示第k个碎片的散射点数目。只看式（7）中的相位部分，如式（8）所示：

对式（8）做微分，可以得到瞬时多普勒为：

当观测目标超过1个周期时，可以得到第k个目标的第i个散射点的多普勒带宽为：

式（10）表明信号的多普勒带宽与碎片的转速、半径、发射信号波长有关系，信号波长越短，目标转速越快，其多普勒带宽就越大，可得到更高的2维分辨率［8，9］。然而，更大的多普勒带宽在低重频系统下也会引起多普勒模糊问题，若要实现对碎片目标的2维成像，需要解决多普勒带宽的模糊问题，或者提高采样率。

3 基于稀疏重采样的空间碎片群目标成像方法

3.1空间碎片的转速估计

由于空间碎片在成像时间内通常做高速自旋运动，其回波在观测时间内会表现出明显的周期性质。假定有3个碎片目标，每个碎片目标的转速分别为ω1=8.20 rad/s， ω2=5.25 rad/s， ω3=2.15 rad/s。回波的接收数据窗形式如图2所示。考虑空间碎片高速自旋时出现的阴影效应，假定所能观测到的成像角度区间为180°。由图2可知，3个碎片的接收数据存在明显的周期性。

在观测时间内，若对回波数据进行自相关处理，在各个目标所在的周期及其整数倍的时刻会出现明显的自项峰值。定义自相关函数为：

图 2 观测数据Fig. 2 Radar echo during the measurement duration

根据信号处理可知，自相关函数的时域表示可以通过信号的傅里叶变换与其共轭相乘的逆傅里叶变换得到。对3个目标回波的自相关函数进行峰值提取，并利用转速对峰值进行扫描，可以得到各个目标的转速信息。图3给出了信噪比为10 dB时的自相关峰值与转速匹配结果，其中圆圈、十字与方块3个线条分别代表了3个碎片的周期扫描线。从图3可知，在目标个数已知的前提下，3个转速能够匹配所有的自相关峰值。

图 3 自相关峰值及不同转速的扫描周期Fig. 3 Peaks of auto-correlation and corresponding scanning cycle

3.2空间碎片群目标成像方法

图2给出了碎片接收的回波数据形式，由于目标的转动，需要考虑目标的遮挡效应，如文献［16］所述，散射点的阴影效应表现为：在一段持续的时间内，后向散射系数比较大，在其它的位置后向散射系数接近零。即雷达在观测自旋目标时，会有一个加窗效应。假定此处的有效观测角度为150°，其它角度范围内的散射系数为零。对于第k个碎片而言，其等效的回波形式如下：

为了抽取出属于第k个碎片的回波数据，可以结合压缩感知理论，利用测量矩阵对观测数据进行抽取，从而提高对空间碎片的成像质量。

假设第k个碎片目标在观测时间内旋转L个周期，将方位慢时间用离散时间表示，可得

其中Ta表示脉冲重复间隔，M表示观测的离散时间长度。式（12）可以表示为如下的离散形式：

其中ak（PQ×1）对应第k个碎片的散射中心强度，Fk矩阵（M×PQ）是根据第k个碎片的回波模型设计的字典，其具体形式如下：

Pk和Qk分别表示第k个目标在距离和方位分辨单元总数。Hk是根据窗效应及旋转周期建立的观测矩阵，其形式由碎片目标的转速ωk与自转周期决定。矩阵形式为：由于空间碎片目标在成像场景内表现出了明显的稀疏特性，可以结合压缩感知理论得到空间碎片的高分辨清晰图像。通过解决如下的数学优化问题，建立稀疏恢复与成像之间的关系：

其中η约束了数据恢复过程中的噪声量级。解决如上所述的l0范数问题是NP问题。本文结合SL0算法［12］，通过高斯函数近似目标函数中的l0范数，从而结合凸优化的方法求解式（19）。一旦得到了稀疏系数ak，将稀疏矢量ak重排成pk×qk维矩阵，即得到第k个目标的2维像I（pk，qk）。然而，回波分量是由多个碎片目标的回波组合得到的，此时，无法得到属于某个碎片的独立回波分量，在建立字典时，需要联合各个目标的字典，以确保字典的冗余性，构建如下的字典：

问题式（19）可以重新写为：

稀疏重构式（21）所述的问题，可得到对应第k个碎片目标回波的稀疏向量a，因为a对应的是观测矩阵Hk下的稀疏矢量，需要对a进行抽取，将a等间隔地分成K份，对其中的第k段进行重排，即可得到第k个目标的2维图像，此时，其他碎片目标的回波能量会落到各自字典所对应的稀疏区间，可以降低多目标问题下彼此能量渗透现象。为了同时得到多个碎片目标的2维像，可以根据自相关函数估计得到的转速设计多组观测矩阵，因为稀疏重构的过程是彼此独立的，可以同时得到多个目标的2维像。具体算法流程如图4所示。

4 仿真实验分析

4.1本文方法实验结果

以空间平台星载雷达对空间碎片群目标成像为仿真背景，通过仿真实验验证本文方法的有效性。假定观测范围内存在3个高速自旋的碎片目标，转动角速度分别为ω1=8.20 rad/s， ω2=5.25 rad/s，ω3=2.15 rad/s，各个碎片的散射模型如图5所示，其散射幅度在观测时间内保持不变，均为1。

雷达发射信号的载频为30 GHz，脉冲重复频率200 Hz。若观测目标的半径为10 cm，针对3个目标的转速，如式（10）所述，对应的多普勒带宽分别为324.8 Hz， 204.8 Hz以及86.0 Hz。当目标的转速和尺寸更大，且载频更高时，多普勒带宽会更大。由于高速自转的碎片目标会引入比较大的多普勒带宽，当带宽大于脉冲重复频率，会出现多普勒模糊。采用文献［10］所述的方法，能够结合稀疏抽取与目标自旋特性提高等效采样率，从而克服高速自旋目标引起的多普勒模糊。

图 4 空间碎片群目标成像流程图Fig. 4 Flow chart of group space debris imaging

图 5 碎片的散射模型Fig. 5 Scattering model of debris

图6（a）-图6（c）是采用本文方法得到的3个碎片目标的2维图像，图6（d）-图6（f）是采用文献［10］所述方法得到的3个碎片目标的2维图像。从图6可以发现，采用文献［10］的方法无法得到清晰的碎片图像，会出现一些虚假的点，或者目标的部分散射点消失，而本文方法能够得到各个碎片更清晰的图像。由于空间碎片多以群目标的形式出现，因此接收到的信号是由多个目标的回波分量组成的。文献［10］中构建的字典是不完备的，无法包含所有目标的信息，此时，当针对某一个目标进行重构时，其余目标的能量也会泄露到该重构信号中，导致无法得到目标的清晰图像。本文先结合自相关函数及转速配对实现目标个数及其对应转速的估计，并构建群目标条件下的完备字典，结合测量矩阵获取不同目标的有效观测数据，这样其他目标的大部分能量都会渗透到属于自己所在的块字典内。由于所构建的字典并非传统的小波基或者傅里叶基，字典的RIP（Restricted Isometry Property，限制等距条件）性质会变差，导致剩余的小部分能量泄露，但是不会影响到目标的成像结果。

图 6 采用本文方法与文献［10］方法得到的各碎片目标像Fig. 6 Space debris imaging using our proposed method and algorithm proposed in Ref. ［10］

4.2不同的稀疏重构方法所得目标像

本文针对的是空间高速自旋目标成像，其回波相位如式（8）所示，字典中每个原子的时频关系满足正弦关系，并非常用的傅里叶基或者小波基等，感知矩阵的RIP性质会有所下降，不利于稀疏重构。文献［12］所提的SL0方法通过平滑高斯函数近似l0范数，能够将稀疏重构问题转化为凸优化问题，并且降低对字典RIP性质的要求，在本文的应用背景下，更适合求解如式（20）所述的稀疏重构问题。图7对比了采用SL0方法和经典OMP（Orthogonal Matching Pursuit）方法得到的各个碎片的目标像。图7（d）-图7（f）表明，各个碎片目标也能够实现稀疏恢复，但碎片2和碎片3的散射点中存在一些非目标散射点，不利于目标的识别。因为构建的感知矩阵RIP性质不如典型的傅里叶基或者小波基，匹配追踪过程会出现信号能量的交错现象［11］。然而，SL0算法降低了对RIP性质的要求，本文采用SL0方法能够得到比较清晰的目标像。

5 结论

图 7 结合不同稀疏重构方法与本文所提模型得到的目标像Fig. 7 Imaging of targets using different sparse recovery method in our proposed model

针对空间平台雷达在低重频、窄带条件下的空间碎片群目标成像问题，本文结合群目标的高速自旋特性，提出了一种稀疏重采样的碎片群目标成像方法。在目标个数已知的前提下，该方法对回波进行自相关处理，通过转速配对获取观测范围内与碎片对应的转速，并结合不同的转速构建联合字典，利用观测矩阵抽取不同碎片目标的回波数据，在克服遮挡效应的同时，实现了空间的等效插值，从而避免了低重频系统下的多普勒模糊问题。仿真实验结果证明了本文方法针对碎片群目标成像的有效性。然而，本文方法在信噪比低于10 dB、回波中存在强弱目标的情况下，无法准确估计碎片目标的转速，不利于设计冗余字典，严重影响目标的成像质量。另外，如何设计观测字典，在抽取数据的同时得到更优RIP性质的感知矩阵，是下一步需要进行的工作。

［1］Abdel-Aziz Y A. An analytical theory for avoidance collision between space debris and operating satellites in LEO［J］. Applied Mathematical Modeling， 2013， 37（18/19）：8283-8291.

［2］Xing Meng-dao and Bao Zheng. High resolution ISAR imaging of high speed moving targets［J］. IEE Proceedings-Radar， Sonar and Navigation， 2005， 152（2）： 58-67.

［3］李道京， 刘波， 尹建凤， 等. 天基毫米波空间碎片观测雷达系统分析与设计［J］. 宇航学报， 2010， 31（12）： 2746-2753. Li Dao-jing， Liu Bo， Yin Jian-feng， et al.. Analysis and design of spaceborne MMW radar for space debris observation system［J］. Journal of Astronautics， 2010，31（12）： 2746-2753.

［4］王洋， 陈建文， 刘中， 等， 多运动目标ISAR成像方法研究［J］. 宇航学报， 2005， 26（4）： 450-454. Wang Yang， Chen Jian-wen， Liu Zhong， et al.. Research on ISAR imaging of multiple moving targets［J］. Journal of Astronautics， 2005， 26（4）： 450-454.

［5］Xiao Da， Su Fu-lin， and Wu Jia-wei. A method of ISAR imaging for multiple targets［C］. 2012 IEEE 11th International Conference on Signal Processing （ICSP）， 2012，3： 2011-2015.

［6］Sato T. Shape estimation of space debris using single-range Doppler interferometry［J］. IEEE Transactions on Geoscience and Remote Sensing， 1999， 37（2）： 1000-1005.

［7］Zhang Lei， Li Ya-chao， Liu Yan， et al.. Time-frequency characteristics based motion estimation and imaging for high speed spinning targets via narrowband waveforms［J］. SCIENCE CHINA Information Sciences， 2010， 53（8）：1628-1640.

［8］Wang Qi， Xing Meng-dao， Lu Guang-yue， et al.. SRMFCLEAN imaging algorithm for space debris［J］. IEEE Transactions on Antennas and Propagation， 2007， 55（12）：3524-3533.

［9］Wang Qi， Xing Meng-dao， Lu Guang-yue， et al.. Single range matching filtering for space debris radar imaging［J］. IEEE Geoscience and Remote Sensing Letters， 2007， 4（4）：576-580.

［10］Zhu Jiang， Zhu Sheng-qi， and Liao Gui-sheng. Highresolution radar imaging of space debris based on sparse representation［J］. IEEE Geoscience and Remote Sensing Letters， 2015， 12（10）： 2090-2094.

［11］Candès E J and Wakin M B. An Introduction To Compressive Sampling［J］. Signal Processing Magazine， 2008，25（2）： 21-30.

［12］Mohimani G H， Babaie-Zadeh M， and Jutten C. A fast approach for overcomplete sparse decomposition based on smoothed l0norm［J］. IEEE Transactions on Signal Processing， 2009， 57（1）： 289-301.

［13］Zhu Dai-yin， Wang Ling， Yu Yu-sheng， et al.. Robust ISAR range alignment via minimizing the entropy of the average range profile［J］. IEEE Geoscience and Remote Sensing Letters， 2009， 6（2）： 204-208.

［14］Ye W， Yeo T S， and Bao Zheng. Weighted least-squares estimation of phase errors for SAR/ISAR autofocus［J］. IEEE Transactions on Geoscience and Remote Sensing，1999， 37（5）： 2487-2494.

［15］Gerry M J， Potter L C， Gupta I J， et al.. A parametric model for synthetic aperture radar measurements［J］. IEEE Transactions on Antennas and Propagation， 1999， 47（7）：1179-1188.

［16］Varshney K R， Çetin M， Fisher J W， et al.. Sparse representation in structured dictionaries with application to synthetic aperture radar［J］. IEEE Transactions on Signal Processing， 2008， 56（8）： 3548-3561.

朱 江（1989-），男，陕西人，西安电子科技大学博士生，主要从事压缩感知理论、贝叶斯稀疏重构理论、结构稀疏理论、高速运动平台下的目标运动补偿方法及稀疏ISAR成像技术方面研究。

E-mail： jiang_z_2012@163.com

邓佳欣（1991-），女，陕西人，西安电子科技大学硕士，中国飞行实验研究院助理工程师，主要从事压缩感知理论及其在宽带信号DOA估计、相干源信号DOA估计应用以及空域雷达目标散射特性、航迹拟合等方面研究。

E-mail： jxdeng0501@163.com

廖桂生（1963-），男，广西人，西安电子科技大学教授，博士生导师，主要从事自适应信号处理、阵列信号处理、信号检测与估计和智能天线信号处理技术方面研究。

E-mail： liaogs@xidian.edu.cn

朱圣棋（1984-），男，江西人，西安电子科技大学副教授，博士生导师，主要从事机载/星载雷达地面运动目标检测与定位方面研究。

E-mail： zhushengqi8@163.com

Space Group Debris Imaging Based on Sparse Sample

Zhu Jiang①Deng Jiaxin②Liao Guisheng①Zhu Shengqi①

①（National Laboratory of Radar Signal Processing， Xidian University， Xi，an 710071， China）

②（China Flight Test Establishment，710089， China）

Space group debris imaging is difficult with sparse data in low Pulse Repetition Frequency （PRF）spaceborne radar. To solve this problem in the narrow band system， we propose a method for space group debris imaging based on sparse samples. Due to the diversity of mass， density， and other factors， space group debris typically rotates at a high speed in different ways. We can obtain angular velocity through the autocorrelation function based on the diversity in the angular velocity. The scattering field usually presents strong sparsity， so we can utilize the corresponding measurement matrix to extract the data of different debris and then combine it using the sparse method to reconstruct the image. Furthermore， we can solve the Doppler ambiguity with the measurement matrix in low PRF systems and suppress some energy of other debris. Theoretical analysis confirms the validity of this methodology. Our simulation results demonstrate that the proposed method can achieve high-resolution Inverse Synthetic Aperture Radar （ISAR） images of space group debris in low PRF systems.

Auto-correlation function； Doppler ambiguity； Sparse sample； Group targets of space debris； ISAR imaging

s： The National Natural Science Foundation of China （61231017）， National Basic Research Program of China（91438106）， Shaanxi Youth Science and Technology New Star （2014KJXX-48）， Natural Fund of Shaanxi Province（2015JQ6206）

TN957.52

A

2095-283X（2016）01-0082-08

10.12000/JR16012

2016-01-15；改回日期：2016-01-28；网络出版：2016-02-18

朱江 jiang_z_2012@163.com

国家自然科学基金重点项目（61231017），国家重大研究计划（91438106），陕西青年科技新星（2014KJXX-48），陕西省自然基金（2015JQ6206）