基于分布式压缩感知的高分宽幅SAR动目标成像技术

潘洁; 王帅; 李道京; 卢晓春

doi:10.12000/JR19060

基于分布式压缩感知的高分宽幅SAR动目标成像技术

DOI: 10.12000/JR19060

潘洁^1,3, ,,
王帅^1,2,3,,
李道京^1,2,,
卢晓春^4,

1.
中国科学院空天信息创新研究院北京 100094
2.
微波成像技术重点实验室北京 100190
3.
中国科学院大学北京 100049
4.
中国科学院国家授时中心西安 710600

基金项目: 国家高分重大专项共性关键技术(30-H30C01-9004-19/21)

详细信息

作者简介:
潘　洁(1977–)，女，四川人，中国科学院空天信息创新研究院高级工程师，博士生，研究方向为稀疏阵列雷达系统。E-mail: panjie@aircas.ac.cn

王　帅(1992–)，男，安徽阜阳人，中国科学院空天信息创新研究院博士生，研究方向为合成孔径雷达成像技术。E-mail: wangshuai161@mils.ucas.edu.cn

李道京(1964–)，男，陕西西安人，中国科学院空天信息创新研究院研究员，博士生导师，主要研究方向为雷达系统和雷达信号处理。E-mail: lidj@mail.ie.ac.cn

卢晓春(1970–)，女，中国科学院授时中心研究员，博士生导师，研究方向为精密时间信息传输与信息处理。E-mail: luxc@ntsc.ac.cn

通讯作者:
潘洁 panjie@aircas.ac.cn

中图分类号: TN951
计量
- 文章访问数: 2573
- HTML全文浏览量: 917
- PDF下载量: 305
- 被引次数: 8
出版历程
- 收稿日期: 2019-06-14
- 修回日期: 2019-11-20
- 网络出版日期: 2020-02-28

High-resolution Wide-swath SAR Moving Target Imaging Technology Based on Distributed Compressed Sensing

PAN Jie^{1,3
, ,},
WANG Shuai^{1,2,3
,},
LI Daojing^{1,2
,},
LU Xiaochun^4
,

1.
Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
2.
National Key Lab of Microwave Imaging Technology, Beijing 100190, China
3.
University of Chinese Academy of Sciences, Beijing 100049, China
4.
National Time Service Center, Chinese Academy of Sciences, Xian 710600, China

Funds: The Key Standard Technologies of National High Resolution Special (30-H30C01-9004-19/21)

More Information

Corresponding author: PAN Jie, panjie@aircas.ac.cn

摘要

摘要: 高分宽幅SAR动目标成像对目标跟踪具有重要的意义，常规天基多通道SAR技术要实现高分宽幅动目标成像需要通道数量巨大，系统复杂度过高，而且图像在方位向存在成对回波，形成虚警。针对上述问题，该文提出了一种基于分布式压缩感知的高分宽幅SAR动目标成像技术，在通道数量较大时，通道数量相比常规高分宽幅动目标成像构型通道数量约降低1倍，利用动目标稀疏特性和杂波背景非稀疏特性构建分布式压缩感知观测模型，采用先方位1维分布式压缩感知重建再距离方位2维分布式压缩感知重建，实现杂波背景和稀疏动目标的重建，并抑制多通道SAR动目标成像中的成对回波。结合RADAR-SAT数据的仿真试验结果验证了该技术的有效性。
- 合成孔径雷达 /
- 多通道 /
- 高分辨率宽测绘带 /
- 分布式压缩感知 /
- 运动目标显示
Abstract: High-resolution wide-swath SAR moving target imaging is of great significance for target tracking. To achieve target tracking, conversional space-based multichannel SAR technology requires a large number of channels. However, this leads to high system complexity. Moreover, paired false target echoes exist in the azimuth. To address such problems, a high-resolution wide-swath SAR moving target imaging technology based on distributed compressed sensing is proposed in this paper. When the number of channels is large enough, the number of channels is approximately half of that of the conventional high-resolution wide-swath target imaging configuration channel. A distributed compressed sensing observation model is constructed using the sparse property of moving targets and the nonsparse characteristics of the clutter background. Clutter background and sparse moving targets are reconstructed by combining one-dimensional distributed compressed sensing reconstruction in the azimuth and two-dimensional distributed compressed sensing reconstruction in the range-azimuth. Moreover, the paired false target echoes in multichannel SAR moving target are suppressed. The simulation results combined with RADAR-SAT data verify the effectiveness of the proposed technology.
- SAR /
- Multi-channel /
- High resolution wide swath /
- Distributed compressed sensing /
- Ground moving target indication

HTML全文

1. 引言

传统单通道SAR成像受最小天线面积的限制^[1]，高分辨率和宽测绘带不能同时实现，在实际应用中需要根据需求对两者进行折衷。而高分宽幅SAR系统具有突破传统单通道SAR系统最小天线面积限制的能力，同时实现高分辨率和宽的测绘幅宽，是未来SAR的一个重要发展趋势。但高分宽幅SAR系统面临的难题是，高的方位向分辨率需要高的脉冲重复频率(Pulse Repetition Frequency, PRF)来避免方位多普勒谱模糊^[2]，而宽的距离测绘带则需要低的PRF。方位向多通道SAR可以在不提升发射脉冲重复频率的情况下实现方位采样率提升，从而实现距离宽测绘带和方位高分辨率，解决了高分宽幅成像的难题^[3-5]，同时该系统具有地面动目标显示(Ground Moving Target Indication, GMTI)能力，因此将两种功能结合可实现高分宽幅动目标成像，在军事侦察中具有重要意义^[6-8]。

在方位多通道SAR系统中，由于系统采用低脉冲重复频率，每个通道接收的回波存在多普勒模糊^[9]，多通道成像存在周期性相位误差造成成对回波，导致传统的动目标成像方法不再适用。Stefan等人^[10]基于最大信杂噪比准则估计运动目标速度，该方法鲁棒性较强，但需要对所有运动目标的可能速度都进行遍历搜索，运算量大大增加。孙光才等人^[11]提出了一种瞬时距离-多普勒(Instantaneous-Range-Doppler, IRD)的动目标成像方法，基于方位去斜Keystone变换^[12]实现聚焦，不依赖于目标精确的运动参数，但是对信噪比要求较高且难以估计目标位置。张双喜等人^[13]利用相位误差随距离频率变化的特性求解运动目标的径向速度，未考虑动目标和杂波背景联合优化成像问题。郭怡冉等人^[14]提出了基于周期性相位误差的分数阶傅里叶变换(FRactional Fourier Transform, FRFT)等径向速度估计算法，但需要系统具有较多的方位向通道。

常规多通道高分宽幅SAR动目标成像系统复杂程度高，数据冗余量大，而且存在动目标成对回波问题，本文针对上述问题提出了一种基于分布式压缩感知^[15]的高分宽幅SAR动目标成像技术，利用动目标稀疏特性和杂波背景非稀疏特性构建分布式压缩感知观测模型，通过优化求解实现杂波背景和稀疏动目标的重建，并抑制多通道SAR动目标成像中的成对回波，基于RADAR-SAT数据的仿真试验结果验证了本文算法的有效性。

2. 多通道高分宽幅SAR动目标成像几何模型和信号模型

多通道高分宽幅SAR动目标成像几何模型如图1所示，设SAR系统在方位向上有N个通道，采用单个通道发射一个宽波束信号，N个通道同时接收回波信号，对于所观测的场景如图1所示，其中$\left( {x,y,z} \right)$表示地面观测单元的中心坐标，$\varphi $表示入射角，V表示平台速度，v表示目标运动速度，${X},{Y},{Z}$方向分别表示顺轨方向，交轨方向和径向。

图 1 SAR系统坐标系

Figure 1. Coordinate system of SAR

下载: 全尺寸图片幻灯片

对于如图1所示目标，第$n$个通道接收到的回波

$\begin{split} {s_n}\left( {{t_{\rm r}},{t_{\rm{a}}}} \right) =\,& \sigma \left( {x,y,z} \right){\omega _{\rm r}}\left( {{t_{\rm r}} - \frac{{{r_n} + {r_0}}}{{\rm{c}}}} \right){\omega _{\rm{a}}}\left( {{t_{\rm{a}}} - {t_{\rm ac}}} \right)\\ & \cdot{{\rm{e}}^{\left( - {\rm{j}}2\pi {f_0}\frac{{{r_n} + {r_0}}}{{\rm{c}}}\right)}}{{\rm{e}}^{\left( + {\rm{j}}\pi {K_{\rm r}}{{\left( {{t_{\rm r}} - \frac{{{r_n} + {r_0}}}{{\rm{c}}}} \right)}^2}\right)}}\\[-15pt] \end{split}$

(1)

${r_0}\left( {{t_{\rm{a}}},x,y,z} \right) \approx \sqrt {{{\left( {x - V{t_{\rm{a}}}} \right)}^2} + {y^2} + {z^2}} + {v_{\rm r}}{t_{\rm{a}}} \hspace{15pt}$

(2)

${r_n}\left( {{t_{\rm{a}}},x,y,z} \right) \approx \sqrt {{{\left( {x - {x_n} - V{t_{\rm{a}}}} \right)}^2} + {y^2} + {z^2}} + {v_{\rm r}}{t_{\rm{a}}}$

(3)

其中，${t_{\rm{a}}}$表示方位慢时间，${\rm c}$表示光速，${f_0}$为工作频率，${r_n}$为地面目标到第$n$个通道的瞬时距离，${r_0}$为地面目标到参考通道的瞬时距离，${v_{\rm r}}$为地面目标的径向速度，$\sigma \left( {x,y,z} \right)$为目标的复反射系数，${\omega _{\rm r}}$表示距离包络，${\omega _{\rm{a}}}\left( {{t_{\rm{a}}}} \right)$表示天线方向图对应方位包络，${t_{\rm ac}}$方向图波束中心偏离时间，${x_n} = \left( {n - 1} \right)d$为阵元的位置，$d$为阵元间隔。

3. 常规多通道高分宽幅SAR动目标成像技术

在常规多通道高分宽幅SAR动目标成像系统中，为了实现高分宽幅成像需要方位向具有多个通道，为了实现动目标检测和成像需要保障脉间方位多通道重叠度形成顺轨干涉，通道构型如图2所示，并给出了相邻脉冲重复时间(Pulse Repetition Time, PRT)内接收天线和等效相位中心的相对位置关系，其中方框表示接收天线，相邻接收天线以不同颜色进行区分；由于雷达系统采用的是收发分置，需要进行双站到单站的等效处理，等效为每个通道自发自收，图2中圆形表示等效相位中心，相邻等效相位中心与上述接收天线依次相对应，并以不同颜色进行区分。为了保障系统的性能指标，正侧视情况下通道数量和系统参数之间的关系如下

图 2 多通道SAR高分宽幅动目标成像通道构型

Figure 2. Channel configuration of multi-channel SAR high resolution wide swath target imaging

下载: 全尺寸图片幻灯片

$N \ge \frac{{4V{\rm{B}}\sin \varphi }}{{{\rho _{\rm{a}}}{\rm{c}}}}$

(4)

其中，B为测绘带宽，c为光速，${\rho _{\rm a}}$为方位向分辨率。假设平台对地速度5 km/s，系统测绘带宽200 km，方位分辨率1 m，入射角为${45^ {\circ} }$，则常规高分宽幅动目标成像要求通道数量大于等于10个，对应阵列天线方位向长度超过20 m，这就导致了系统极度复杂。

此外，高分宽幅动目标成像还存在成对回波问题。在成像过程中根据静止目标进行匹配滤波处理，因此对于上述动目标来讲其残余相位如下

$\begin{split} \Delta {\phi _n}\left( {{t_{\rm{a}}}} \right) = - 2\pi {f_0}\frac{{2{v_{\rm r}}{t_{\rm{a}}}}}{{\rm{c}}} + \pi {K_{\rm r}}\frac{{4{v_{\rm r}}{t_{\rm{a}}}}}{{\rm{c}}}\left( {{t_{\rm r}} - \frac{{{r_0} + {r_n}}}{{\rm{c}}}} \right)\\ \end{split}$

(5)

其中第1项表现为多普勒频率，第2项表现为距离徙动，经过距离向脉冲压缩之后残余相位只剩下多普勒频率，即为

$\Delta {\phi _{n,{\rm{rc}}}}\left( {{t_{\rm{a}}}} \right) = - 2\pi {f_0}\frac{{2{v_{\rm r}}{t_{\rm{a}}}}}{{\rm{c}}}$

(6)

多个通道的残余相位是一致的，而方位向不同时间残余相位是不一致的，由此形成周期性相位误差，经过方位向脉冲压缩之后，运动目标不仅在方位向形成位置偏移，而且形成了较高的栅瓣，即成对回波，在高分宽幅动目标成像中形成图像模糊，并造成动目标检测的虚警^[14]。

4. 基于分布式压缩感知的高分宽幅SAR动目标成像技术

针对常规高分宽幅动目标成像技术存在的系统复杂度高、成像模糊等问题，本文提出了基于分布式压缩感知的高分宽幅动目标成像技术，利用运动目标稀疏分布特性和杂波背景非稀疏特性构造分布式压缩感知模型，只需少量通道即可实现回波数据获取，通过分布式压缩感知稀疏重建获得动目标和杂波背景信息，同时实现动目标成对回波抑制和观测场景高分宽幅成像。该方法利用了动目标稀疏特性，大大降低了系统复杂度和信息冗余度，并在信号模型上考虑了动目标残余相位特性，抑制了成对回波，实现了无模糊高分宽幅成像。

4.1 基于分布式压缩感知的高分宽幅动目标信号模型

雷达观测场景如图1所示，整个场景包括动目标和杂波背景的回波表示如下

$\begin{split} {s_n}\left( {{t_{\rm r}},{t_{\rm a}}} \right) =\,& \sum\limits_i \sigma \left( {{x_i},{y_i},{z_i}} \right){\omega _{\rm r}}\left( {{t_{\rm r}} - \frac{{{r_{ni}} + {r_{0i}}}}{{\rm c}}} \right)\\ & \cdot {\omega _{\rm{a}}}\left( {{t_{\rm a}} - {t_{{\rm{ac}}i}}} \right){{\rm{e}}^{\left( - {\rm{j}}2\pi {f_0}\frac{{{r_{ni}} + {r_{0i}}}}{{\rm c}}\right)}}\\ \\ & \cdot {{\rm{e}}^{\left( + {\rm{j}}\pi {K_{\rm r}}{{\left( {{t_{\rm r}} - \frac{{{r_{ni}} + {r_{0i}}}}{{\rm c}}} \right)}^2}\right)}} + \sum\limits_j \sigma \left( {{x_j},{y_j},{z_j}} \right)\\ & \cdot{\omega _{\rm r}}\left( {{t_{\rm r}} - \frac{{{r_{nj}} + {r_{0j}}}}{{\rm{c}}}} \right){\omega _{\rm{a}}}\left( {{t_{\rm a}} - {t_{{\rm ac}i}}} \right)\\ & \cdot{{\rm{e}}^{\left( - {\rm{j}}2\pi {f_0}\frac{{{r_{nj}} + {r_{0j}}}}{{\rm c}}\right)}}{{\rm{e}}^{\left( + {\rm{j}}\pi {K_{\rm r}}{{\left( {{t_{\rm r}} - \frac{{{r_{nj}} + {r_{0j}}}}{{\rm c}}} \right)}^2}\right)}} \end{split}$

(7)

其中第1部分为杂波背景，第2部分为动目标。式(7)中的距离表达式如下

${r_{0i}} \approx \sqrt {{{\left( {{x_i} - V{t_{\rm{a}}}} \right)}^{ \!\! { \scriptsize2}}} + {y_i}^{\!\!\!\! \scriptsize2} + {z_i}^{{\!\!\!\! \scriptsize2}}} \hspace{65pt}$

(8)

${r_{ni}} \approx \sqrt {{{\left( {{x_i} - {x_n} - V{t_{\rm{a}}}} \right)}^{ \!\! { \scriptsize2}}} + {y_i}^{\!\!\!\! \scriptsize2} + {z_i}^{\!\!\!\! \scriptsize2}} \hspace{43pt}$

(9)

${r_{0j}} \approx \sqrt {{{\left( {{x_j} - V{t_{\rm{a}}}} \right)}^{ \!\! { \scriptsize2}}} + {y_j}^{\!\!\!\! \scriptsize2} + {z_j}^{\!\!\!\! \scriptsize2}} + {v_{{\rm r}j}}{t_{\rm{a}}} \hspace{25pt}$

(10)

${r_{nj}} \approx \sqrt {{{\left( {{x_j} - {x_n} - V{t_{\rm{a}}}} \right)}^{ \!\! { \scriptsize2}}} + {y_j}^{\!\!\!\! \scriptsize2} + {z_j}^{\!\!\!\! \scriptsize2}} + {v_{{\rm r}j}}{t_{\rm a}}$

(11)

将式(7)表示为矩阵形式

${{{s}}_n} = {{{\varTheta}} _n}{{\sigma}} = \left[ {\begin{array}{*{20}{c}} {{{{\varTheta}} _{n1}}}&{{{{\varTheta}} _{n2}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{{{\sigma}} _1}} \\ {{{{\sigma}} _2}} \end{array}} \right]$

(12)

其中${{{\varTheta}} _{n1}}$, ${{{\sigma}} _1}$表示杂波背景部分(共同分量)，${{{\varTheta}} _{n2}}$, ${{{\sigma}} _2}$表示运动目标部分(更新分量)。它们表达式如下

${{{s}}_n} = \left[ {\begin{array}{*{20}{c}} {{s_n}\left( {{t_{{\rm r}1}},{t_{{\rm a}1}}} \right)} \\ {\begin{array}{*{20}{c}} {{s_n}\left( {{t_{{\rm r}1}},{t_{{\rm a}2}}} \right)} \\ \vdots \\ {{s_n}\left( {{t_{{\rm r}K}},{t_{{\rm a}M}}} \right)} \end{array}} \end{array}} \right]\hspace{91pt}$

(13)

$\begin{split} {{{\varTheta}} _{n1}} = \left[ {\begin{array}{*{20}{c}} {{\varTheta _{n1}}\left( {{t_{{\rm r}1}},{t_{{\rm a}1}},{\sigma _1}} \right)} \!\!&\!\! ··· \!\!&\!\! {{\varTheta _{n1}}\left( {{t_{{\rm r}1}},{t_{{\rm a}1}},{\sigma _I}} \right)}\\ {{\varTheta _{n1}}\left( {{t_{{\rm r}1}},{t_{{\rm a}2}},{\sigma _1}} \right)} \!\!&\!\! ··· \!\!&\!\! {{\varTheta _{n1}}\left( {{t_{{\rm r}1}},{t_{{\rm a}2}},{\sigma _I}} \right)}\\ \vdots \!\!&\!\! \ddots \!\!&\!\! \vdots \\ {{\varTheta _{n1}}\left( {{t_{{\rm r}K}},{t_{{\rm a}M}},{\sigma _1}} \right)} \!\!&\!\! ··· \!\!&\!\!{{\varTheta _{n1}}\left( {{t_{{\rm r}K}},{t_{{\rm a}M}},{\sigma _I}} \right)} \end{array}} \right]\\ \end{split}$

(14)

$\begin{split} {\varTheta _{n1}}\left( {{t_{\rm r}},{t_{\rm a}},{\sigma _i}} \right) =\,& {\omega _{\rm r}}\left( {{t_{\rm r}} - \frac{{{r_{ni}} + {r_{0i}}}}{{\rm c}}} \right){\omega _{\rm a}}\left( {{t_{\rm a}} - {t_{{\rm ac}i}}} \right)\\ & \cdot {{\rm e}^{\left( - {\rm j}2{\rm \pi }{f_0}\frac{{{r_{ni}} + {r_{0i}}}}{{\rm c}}\right)}}{{\rm e}^{\left( + {\rm {\rm{j}}\pi }{K_{\rm r}}{{\left( {{t_{\rm r}} - \frac{{{r_{ni}} + {r_{0i}}}}{{\rm c}}} \right)}^2}\right)}}\\ \end{split}$

(15)

$\begin{split} {{{\varTheta}} _{n2}} = \left[ {\begin{array}{*{20}{c}} {{\varTheta _{n2}}\left( {{t_{{\rm r}1}},{t_{{\rm a}1}},{\sigma _{{\rm r}1}}} \right)} \!\!\!&\!\!\! ··· \!\!\!&\!\!\! {{\varTheta _{n2}}\left( {{t_{{\rm r}1}},{t_{{\rm a}1}},{\sigma _{{\rm r}J}}} \right)}\\ {{\varTheta _{n2}}\left( {{t_{{\rm r}1}},{t_{{\rm a}2}},{\sigma _{{\rm r}1}}} \right)} \!\!\!&\!\!\! ··· \!\!\!&\!\!\! {{\varTheta _{n2}}\left( {{t_{{\rm r}1}},{t_{{\rm a}2}},{\sigma _{{\rm r}J}}} \right)}\\ \vdots \!\!\!&\!\!\! \ddots \!\!\!&\!\!\! \vdots\\ {{\varTheta _{n2}}\left( {{t_{{\rm r}K}},{t_{{\rm a}M}},{\sigma _{{\rm r}1}}} \right)} \!\!\!&\!\!\! ··· \!\!\!&\!\!\!{{\varTheta _{n2}}\left( {{t_{{\rm r}K}},{t_{{\rm a}M}},{\sigma _{{\rm r}J}}} \right)} \end{array}} \right]\\ \end{split}$

(16)

$\begin{split} {\varTheta _{n2}}\left( {{t_{\rm r}},{t_{\rm a}},{\sigma _j}} \right) =\,& {\omega _{\rm r}}\left( {{t_{\rm r}} - \frac{{{r_{nj}} + {r_{0j}}}}{{\rm c}}} \right){\omega _{\rm a}}\left( {{t_{\rm a}} - {t_{{\rm ac}j}}} \right)\hspace{10pt}\\ & \cdot {{\rm e}^{\left( - {\rm j}2\pi {f_0}\frac{{{r_{nj}} + {r_{0j}}}}{{\rm c}}\right)}}\\ & \cdot {{\rm e}^{\left( + {\rm {\rm{j}}\pi }{K_{\rm r}}{{\left( {{t_{\rm r}} - \frac{{{r_{nj}} + {r_{0j}}}}{{\rm c}}} \right)}^2}\right)}} \\ \approx \;& {\varTheta _{n1}}\left( {{t_{\rm r}},{t_{\rm a}},{\sigma _j}} \right){{\rm e}^{\left( - {\rm j}2\pi \frac{{2{v_{rj}}{t_{\rm a}}}}{\lambda }\right)}}\\[-15pt] \end{split}$

(17)

其中K, M分别表示距离向和方位向采样点数，I, J分别表示杂波背景目标数量和动目标数量。整合多个通道信号形成联合观测模型如下

$\begin{split} &{{S}} ={{\varTheta}} {{\sigma}} ,\;\;{{S}} = \left[ {\begin{array}{*{20}{c}} {{{{s}}_1}} \\ {{{{s}}_2}} \\ \vdots \\ {{{{s}}_N}} \end{array}} \right],\\ &{{\varTheta}} = \left[ {\begin{array}{*{20}{c}} {{{{\varTheta}} _1}} \\ {{{{\varTheta}} _2}} \\ \vdots \\ {{{{\varTheta}} _N}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{{{\varTheta}} _{11}}}&{{{{\varTheta}} _{12}}} \\ {{{{\varTheta}} _{21}}}&{{{{\varTheta}} _{22}}} \\ \vdots & \vdots \\ {{{{\varTheta}} _{N1}}}&{{{{\varTheta}} _{N2}}} \end{array}} \right] \end{split}$

(18)

根据分布式压缩感知理论^[16]，如果观测的数据量满足下面的条件，则观测场景可以很好的被恢复

$KMN \ge {{\rm c}_0}\left( {I + J} \right)$

(19)

其中c₀表示一个常数。考虑高分宽幅极限情况下，参照式(4)可以得到杂波背景目标数量$I \approx $$ \dfrac{{2VB\sin \varphi }}{{{\rho _{\rm a}}{\rm c}}}KM$，常数c₀=1，于是可以得到

$N \ge \frac{{2VB\sin \varphi }}{{{\rho _{\rm a}}{\rm c}}} + \frac{J}{{KM}}$

(20)

由于运动目标一般为稀疏分布，因此式(20)中第2项通常很小，通过和式(4)比较可以发现基于分布式压缩感知的高分宽幅动目标成像技术相比传统技术所需通道数量约降低1倍。

4.2 基于分布式压缩感知的高分宽幅数据录取

基于分布式压缩感知的高分宽幅数据录取可以采用多种构型，本文构建了一种顺轨多通道构型。该构型与图1中观测几何相似，与常规多通道高分宽幅SAR动目标成像原理相比，基于分布式压缩感知的高分宽幅动目标成像通道构型如图3所示，并给出了相邻PRT内接收天线和等效相位中心的相对位置关系，其中方框表示接收天线，相邻接收天线以不同颜色进行区分；由于雷达系统采用的是收发分置，需要进行双站到单站的等效处理，等效为每个通道自发自收，图3中圆形表示等效相位中心，相邻等效相位中心与上述接收天线依次相对应，并以不同颜色进行区分。区别在于基于分布式压缩感知的SAR系统通道数量更少，对应脉间方位多通道重叠度更低，系统通道数量越多，脉间方位通道重叠数量越少，则要求的通道数量越少，极限情况下本文构型通道数量相比常规高分宽幅动目标成像构型通道数量约降低1倍。

图 3 基于分布式压缩感知的高分宽幅动目标成像通道构型

Figure 3. Channel configuration of high resolution wide swath target imaging based on distributed compressed sensing

下载: 全尺寸图片幻灯片

需要说明的是，当相邻PRF之间不存在完全重合的等效相位中心时，利用分布式压缩感知技术也可以实现高分宽幅动目标成像，本文不再进行详细讨论。

4.3 基于分布式压缩感知的高分宽幅动目标成像处理

由于基于分布式压缩感知的高分宽幅SAR动目标成像系统脉间方位多通道重叠度不足以支撑顺轨干涉处理，且常规高分宽幅动目标成像存在成对回波问题，因此需要结合动目标稀疏特性进行信号处理获得高分宽幅动目标图像。

基于分布式压缩感知的高分宽幅动目标信号模型如式(18)所示，在信号重建过程中需要加入动目标稀疏性约束，表示如下

$\min {\left\| {{{{\sigma}} _2}} \right\|_0},\;{\rm{s.t.}}\;\;{{S}}{\rm{ = }}{{\varTheta}} \left[ {\begin{array}{*{20}{c}} {{{{\sigma}} _1}} \\ {{{{\sigma}} _2}} \end{array}} \right]$

(21)

其中，${{{\sigma}} _1}$表示杂波背景部分，是非稀疏的，${{{\sigma}} _2}$表示运动目标部分，是稀疏的。理论上，式(21)的解可以通过优化方法获得，从而实现场景的重建。实际问题中，高分宽幅成像目标数量巨大，而且存在明显的距离方位耦合，难以实现单维度处理，因此直接求解该问题难以实现。

针对上述问题，本文采用先方位1维分布式压缩感知重建再距离方位2维分布式压缩感知重建的过程。方位2维分布式压缩感知通过数据降维大大降低计算复杂度，并估计得到目标的速度和位置参数，但存在方位成对回波；距离方位2维分布式压缩感知根据运动目标速度和位置信息降低字典维度从而降低计算复杂度，完成成对回波抑制。算法的流程详见表1。

表 1 算法流程

Table 1. Algorithm flowchart

(1) 对高分宽幅SAR动目标成像系统每个通道数据分别进行距离压缩处理，获得每个通道的1维距离图像；
(2) 联合多通道数据转换到方位多普勒域进行距离徙动校正，之后转换回到方位时域，完成对静止目标距离徙动校正，而动目标由于相位　　误差存在距离徙动校正误差；
(3) 对距离徙动校正后的多通道数据构造方位1维分布式压缩感知模型，利用优化方法进行信号重建，分别获得非稀疏背景杂波和稀疏动目　　标散射信息；
(4) 根据稀疏动目标散射信息估计动目标的位置和速度参数；
(5) 对距离徙动校正后的多通道数据根据运动目标位置和速度信息构造距离方位2维分布式压缩感知模型，利用优化方法进行信号重建，获　　得非稀疏背景杂波和稀疏动目标散射信息；
(6) 将运动目标信息${{{\sigma}} _2}$和背景杂波图像${{{\sigma}} _1}$融合形成完整的高分宽幅图像。

下载: 导出CSV

| 显示表格

基于分布式压缩感知的多通道SAR高分宽幅动目标成像处理流程如图4所示。

图 4 算法流程图

Figure 4. Flowchart of the algorithm

下载: 全尺寸图片幻灯片

5. 仿真实验

为了对基于分布式压缩感知的高分宽幅SAR动目标成像技术有效性进行验证，开展了计算机数值仿真，以RADAR-SAT实测图像作为杂波背景，添加稀疏的动目标，仿真产生多通道SAR回波数据，通过本文成像处理过程获得最终SAR图像。仿真的参数如表2所示

表 2 雷达系统主要参数

Table 2. Main parameters of radar system

参数	数值
工作频率	5.3 GHz
带宽	30 MHz
采样率	32.3 MHz
PRF	1256.98 Hz
速度	7062 m/s
通道数	3
动目标速度	2.5 m/s

下载: 导出CSV

| 显示表格

仿真中采用的实测图像如图5所示，在图中红色圆圈标示的道路位置添加1个运动目标，仿真产生杂波背景和运动目标的回波。

图 5 仿真采用实测图像

Figure 5. Measured image used in simulation

下载: 全尺寸图片幻灯片

对单个通道进行2维成像处理，由于方位向欠采样的存在，单个通道图像存在明显的方位向模糊，图像质量变差，如图6所示。

图 6 单通道成像结果

Figure 6. Imaging result of single channel

下载: 全尺寸图片幻灯片

对多通道距离压缩和距离徙动校正数据进行方位向1维分布式压缩感知建模，利用优化方法进行非稀疏背景杂波和稀疏运动目标重建，所得图像如图7所示。图7(a)表示非稀疏静止背景杂波重建结果，可以看到经过多通道融合之后图像模糊程度得到明显改善，并且运动目标得到了良好的抑制，同时运动目标的成对回波也得到了良好的抑制。图7(b)表示稀疏运动目标重建结果，可以看到杂波得到良好抑制，重建图像中有明显的运动目标和成对回波，动目标检测存在虚警，需要说明的是，当距离徙动现象显著时，成对回波和动目标会落在不同的距离分辨单元内。从图7(b)中可以估计出多个动目标信息，用于构建距离方位分布式压缩感知模型。图7(c)表示非稀疏静止背景杂波和稀疏运动目标融合结果，可以看到，由于目标径向速度影响，目标方位位置出现偏移，成像位置不在道路上。

图 7 方位分布式压缩感知重建结果

Figure 7. Reconstructed results of azimuth distributed compressed sensing

下载: 全尺寸图片幻灯片

利用运动目标信息构建距离方位2维分布式压缩感知模型，利用优化方法进行信号重建，获得非稀疏背景杂波和稀疏动目标散射信息，如图8所示。图8(a)表示非稀疏静止背景杂波重建结果。图8(b)表示稀疏运动目标重建结果，可以看到运动目标的成对回波得到良好的抑制，同时估计得到运动目标径向速度为2.46 m/s，与仿真速度吻合。图8(c)表示非稀疏静止背景杂波和稀疏运动目标融合结果，可以看到运动目标得到了良好的重建，运动目标位置重新定位于道路上，同时成对回波得到了良好的抑制，验证了本文基于分布式压缩感知的高分宽幅SAR动目标成像技术有效性。

图 8 距离方位分布式压缩感知重建结果

Figure 8. Reconstructed results of distributed compressed sensing in range-azimuth

下载: 全尺寸图片幻灯片

6. 结束语

高分宽幅动目标成像是SAR技术发展的一个重要方向，然而传统多通道SAR实现高分宽幅动目标成像需要保障脉间方位多通道重叠度，系统复杂程度过高，而且动目标成像存在成对回波现象，造成动目标检测虚警。针对上述问题，本文提出了一种基于分布式压缩感知的高分宽幅SAR动目标成像技术，结合杂波背景非稀疏特性和运动目标稀疏特性建立分布式压缩感知信号模型，通过模型求解实现高分宽幅动目标成像。针对高分宽幅条件下分布式压缩感知信号模型难以直接求解的问题，本文采用先方位1维分布式压缩感知重建再距离方位2维分布式压缩感知重建的过程，方位1维分布式压缩感知通过数据降维大大降低计算复杂度，并估计得到目标的速度和位置参数，距离方位2维分布式压缩感知根据运动目标速度和位置信息降低字典维度从而降低计算复杂度。基于RADAR-SAT数据的仿真试验表明，本文技术方法能够良好实现运动目标成像和信息反演，验证了该技术的可行性和有效性。

图 1 SAR系统坐标系

Figure 1. Coordinate system of SAR

下载: 全尺寸图片幻灯片

图 2 多通道SAR高分宽幅动目标成像通道构型

Figure 2. Channel configuration of multi-channel SAR high resolution wide swath target imaging

下载: 全尺寸图片幻灯片

图 3 基于分布式压缩感知的高分宽幅动目标成像通道构型

Figure 3. Channel configuration of high resolution wide swath target imaging based on distributed compressed sensing

下载: 全尺寸图片幻灯片

图 4 算法流程图

Figure 4. Flowchart of the algorithm

下载: 全尺寸图片幻灯片

图 5 仿真采用实测图像

Figure 5. Measured image used in simulation

下载: 全尺寸图片幻灯片

图 6 单通道成像结果

Figure 6. Imaging result of single channel

下载: 全尺寸图片幻灯片

图 7 方位分布式压缩感知重建结果

Figure 7. Reconstructed results of azimuth distributed compressed sensing

下载: 全尺寸图片幻灯片

图 8 距离方位分布式压缩感知重建结果

Figure 8. Reconstructed results of distributed compressed sensing in range-azimuth

下载: 全尺寸图片幻灯片

表 1 算法流程

Table 1. Algorithm flowchart

(1) 对高分宽幅SAR动目标成像系统每个通道数据分别进行距离压缩处理，获得每个通道的1维距离图像；
(2) 联合多通道数据转换到方位多普勒域进行距离徙动校正，之后转换回到方位时域，完成对静止目标距离徙动校正，而动目标由于相位　　误差存在距离徙动校正误差；
(3) 对距离徙动校正后的多通道数据构造方位1维分布式压缩感知模型，利用优化方法进行信号重建，分别获得非稀疏背景杂波和稀疏动目　　标散射信息；
(4) 根据稀疏动目标散射信息估计动目标的位置和速度参数；
(5) 对距离徙动校正后的多通道数据根据运动目标位置和速度信息构造距离方位2维分布式压缩感知模型，利用优化方法进行信号重建，获　　得非稀疏背景杂波和稀疏动目标散射信息；
(6) 将运动目标信息${{{\sigma}} _2}$和背景杂波图像${{{\sigma}} _1}$融合形成完整的高分宽幅图像。

下载: 导出CSV

表 2 雷达系统主要参数

Table 2. Main parameters of radar system

参数	数值
工作频率	5.3 GHz
带宽	30 MHz
采样率	32.3 MHz
PRF	1256.98 Hz
速度	7062 m/s
通道数	3
动目标速度	2.5 m/s

下载: 导出CSV

参考文献(16)

[1]	FREEMAN A, JOHNSON W T K, HUNEYCUTT B, et al. The "Myth" of the minimum SAR antenna area constraint[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(1): 320–324. doi: 10.1109/36.823926
[2]	ZHANG Shuangxi, XING Mengdao, XIA Xianggen, et al. Multichannel HRWS SAR imaging based on range-variant channel calibration and multi-doppler-direction restriction ambiguity suppression[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(7): 4306–4327. doi: 10.1109/TGRS.2013.2281329
[3]	叶恺, 禹卫东, 徐伟, 等. 基于柱形抛物面天线的MIMO SAR研究[J]. 电子与信息学报, 2018, 40(8): 1816–1822. doi: 10.11999/JEIT171105 YE Kai, YU Weidong, XU Wei, et al. Investigation on parabolic cylinder reflector based MIMO SAR[J]. Journal of Electronics &Information Technology, 2018, 40(8): 1816–1822. doi: 10.11999/JEIT171105
[4]	DING Chibiao, LIANG Xingdong, WANG Jie, et al. An novel airborne MIMO-SAR system built in IECAS[C]. 2017 IEEE International Geoscience and Remote Sensing Symposium, Fort Worth, USA, 2017: 2423–2426. doi: 10.1109/IGARSS.2017.8127481.
[5]	周建卫, 李道京, 田鹤, 等. 基于共形稀疏阵列的艇载外辐射源雷达性能分析[J]. 电子与信息学报, 2017, 39(5): 1058–1063. doi: 10.11999/JEIT160846 ZHOU Jianwei, LI Daojing, TIAN He, et al. Performance analysis on airship-borne passive radar based on conformal sparse array[J]. Journal of Electronics &Information Technology, 2017, 39(5): 1058–1063. doi: 10.11999/JEIT160846
[6]	KIRK J C, DARDEN S, MAJUMDER U, et al. Forty years of digital SAR and slow GMTI technology[C]. 2014 IEEE Radar Conference, Cincinnati, USA, 2014: 64–69. doi: 10.1109/RADAR.2014.6875556.
[7]	HENDERSON F M and LEWIS A J. Principles and Applications of Imaging Radar[M]. New York: Wiley, 1998. 89–90.
[8]	GIERULL C H. Statistical analysis of multilook SAR interferograms for CFAR detection of ground moving targets[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(4): 691–701. doi: 10.1109/TGRS.2003.821886
[9]	GEBERT N. Multi-channel azimuth processing for high-resolution wide-swath SAR imaging[D]. [Ph.D. dissertation], German Aerospace Center, 2009.
[10]	BAUMGARTNER S V and KRIEGER G. Simultaneous high-resolution wide-swath SAR imaging and ground moving target indication: Processing approaches and system concepts[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(11): 5015–5029. doi: 10.1109/JSTARS.2015.2450019
[11]	SUN Guangcai, XING Mengdao, XIA Xianggen, et al. Robust ground moving-target imaging using deramp-keystone processing[J]. IEEE Transactions on Geoscience and Remote Sensing, 2012, 51(2): 966–982. doi: 10.1109/TGRS.2012.2204889
[12]	LAZAROV A D and MINCHEV C N. SAR imaging of a moving target[C]. The 2007 3rd International Conference on Recent Advances in Space Technologies, Istanbul, Turkey, 2017: 366–372. doi: 10.1109/RAST.2007.4284014.
[13]	ZHANG Shuangxi, XING Mengdao, XIA Xianggen, et al. A robust channel-calibration algorithm for multi-channel in azimuth HRWS SAR imaging based on local maximum-likelihood weighted minimum entropy[J]. IEEE Transactions on Image Processing, 2013, 22(12): 5294–5305. doi: 10.1109/TIP.2013.2274387
[14]	ZHANG Yun, GUO Yiran, JI Zhenyuan, et al. A novel algorithm of radial velocity estimation for MC-HRWS system[C]. 2017 International Conference on Communications, Signal Processing, and Systems, Harbin, China, 2019: 1396–1403. doi: 10.1007/978-981-10-6571-2_168.
[15]	CANDES E J and TAO T. Near-optimal signal recovery from random projections: Universal encoding strategies?[J]. IEEE Transactions on Information Theory, 2006, 52(12): 5406–5425. doi: 10.1109/TIT.2006.885507
[16]	DUARTE M F, SARVOTHAM S, BARON D, et al. Distributed compressed sensing of jointly sparse signals[C]. Conference Record of the Thirty-Ninth Asilomar Conference on Signals, Systems and Computers, Pacific Grove, USA, 2005: 1537–1541. doi: 10.1109/ACSSC.2005.1600024.

施引文献

期刊类型引用(3)

1.	王昕，张权，乔玲，叶啓. 双基地SAR动目标图像相位误差分析和二维稀疏重聚焦处理. 南京邮电大学学报(自然科学版). 2024(04): 54-67 . 百度学术
2.	刘潇文，原方，楚国锋，李凯，东晨. 基于量子算法的逆合成孔径雷达稀疏成像方法. 空军工程大学学报. 2023(01): 26-31+88 . 百度学术
3.	雷禹，冷祥光，孙忠镇，计科峰. 宽幅SAR海上大型运动舰船目标数据集构建及识别性能分析. 雷达学报. 2022(03): 347-362 . 本站查看

其他类型引用(5)

资源附件(0)

访问统计

图(8) / 表(2)

计量

文章访问数: 2573
HTML全文浏览量: 917
PDF下载量: 305
被引次数: 8

1. 引言
2. 多通道高分宽幅SAR动目标成像几何模型和信号模型
3. 常规多通道高分宽幅SAR动目标成像技术
4. 基于分布式压缩感知的高分宽幅SAR动目标成像技术
4.1 基于分布式压缩感知的高分宽幅动目标信号模型
4.2 基于分布式压缩感知的高分宽幅数据录取
4.3 基于分布式压缩感知的高分宽幅动目标成像处理
5. 仿真实验
6. 结束语

1. 引言
2. 多通道高分宽幅SAR动目标成像几何模型和信号模型
3. 常规多通道高分宽幅SAR动目标成像技术
4. 基于分布式压缩感知的高分宽幅SAR动目标成像技术
4.1 基于分布式压缩感知的高分宽幅动目标信号模型
4.2 基于分布式压缩感知的高分宽幅数据录取
4.3 基于分布式压缩感知的高分宽幅动目标成像处理
5. 仿真实验
6. 结束语

参考文献(16)

施引文献

资源附件(0)

访问统计

基于分布式压缩感知的高分宽幅SAR动目标成像技术

DOI: 10.12000/JR19060

通讯作者:
潘洁 panjie@aircas.ac.cn

计量

High-resolution Wide-swath SAR Moving Target Imaging Technology Based on Distributed Compressed Sensing

Corresponding author: PAN Jie, panjie@aircas.ac.cn

1. 引言

2. 多通道高分宽幅SAR动目标成像几何模型和信号模型

3. 常规多通道高分宽幅SAR动目标成像技术

4. 基于分布式压缩感知的高分宽幅SAR动目标成像技术

4.1 基于分布式压缩感知的高分宽幅动目标信号模型

4.2 基于分布式压缩感知的高分宽幅数据录取

4.3 基于分布式压缩感知的高分宽幅动目标成像处理

5. 仿真实验

6. 结束语

期刊类型引用(3)

其他类型引用(5)

计量

目录

1. 引言

2. 多通道高分宽幅SAR动目标成像几何模型和信号模型

3. 常规多通道高分宽幅SAR动目标成像技术

4. 基于分布式压缩感知的高分宽幅SAR动目标成像技术

4.1 基于分布式压缩感知的高分宽幅动目标信号模型

4.2 基于分布式压缩感知的高分宽幅数据录取

4.3 基于分布式压缩感知的高分宽幅动目标成像处理

5. 仿真实验

6. 结束语

期刊介绍

联系我们

基于分布式压缩感知的高分宽幅SAR动目标成像技术

DOI: 10.12000/JR19060

通讯作者: 潘洁 panjie@aircas.ac.cn

计量

出版历程

High-resolution Wide-swath SAR Moving Target Imaging Technology Based on Distributed Compressed Sensing

Corresponding author: PAN Jie, panjie@aircas.ac.cn

1. 引言

2. 多通道高分宽幅SAR动目标成像几何模型和信号模型

3. 常规多通道高分宽幅SAR动目标成像技术

4. 基于分布式压缩感知的高分宽幅SAR动目标成像技术

4.1 基于分布式压缩感知的高分宽幅动目标信号模型

4.2 基于分布式压缩感知的高分宽幅数据录取

4.3 基于分布式压缩感知的高分宽幅动目标成像处理

5. 仿真实验

6. 结束语

期刊类型引用(3)

其他类型引用(5)

计量

出版历程

目录

1. 引言

2. 多通道高分宽幅SAR动目标成像几何模型和信号模型

3. 常规多通道高分宽幅SAR动目标成像技术

4. 基于分布式压缩感知的高分宽幅SAR动目标成像技术

4.1 基于分布式压缩感知的高分宽幅动目标信号模型

4.2 基于分布式压缩感知的高分宽幅数据录取

4.3 基于分布式压缩感知的高分宽幅动目标成像处理

5. 仿真实验

6. 结束语

期刊介绍

联系我们

通讯作者:
潘洁 panjie@aircas.ac.cn