CN105301589B - 高分辨宽测绘带sar地面运动目标成像方法 - Google Patents

高分辨宽测绘带sar地面运动目标成像方法 Download PDF

Info

Publication number
CN105301589B
CN105301589B CN201510726789.3A CN201510726789A CN105301589B CN 105301589 B CN105301589 B CN 105301589B CN 201510726789 A CN201510726789 A CN 201510726789A CN 105301589 B CN105301589 B CN 105301589B
Authority
CN
China
Prior art keywords
moving target
signal
echo
target
estimation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201510726789.3A
Other languages
English (en)
Other versions
CN105301589A (zh
Inventor
杨桃丽
淦小健
王勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201510726789.3A priority Critical patent/CN105301589B/zh
Publication of CN105301589A publication Critical patent/CN105301589A/zh
Application granted granted Critical
Publication of CN105301589B publication Critical patent/CN105301589B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • G01S13/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50—Systems of measurement based on relative movement of target
    • G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于稀疏空间谱估计的高分辨宽测绘带SAR地面运动目标成像方法，主要解决了现有方法在估计目标速度时引入的巨大运算量或谱分量增益损失的问题。本发明的实现步骤是：(1)对回波信号进行方位向傅里叶变换；(2)根据稀疏空间谱估计方法，利用某一多普勒单元的回波信号估计得到运动目标的速度；(3)根据估计得到的目标速度，得到目标的无模糊信号；(4)根据估计得到的目标速度，对无模糊信号进行传统高分辨运动目标成像。本发明将运动目标速度估计问题转换为稀疏空间谱估计问题，仅需较小的运算量即可精确估计得到运动目标速度，从而保证了运动目标多普勒频谱的精确重构，实现了高分辨宽测绘带运动目标SAR成像。

Description

高分辨宽测绘带SAR地面运动目标成像方法
技术领域
本发明属于雷达技术领域，更进一步涉及雷达成像技术领域中的高分辨率宽测绘带合成孔径雷达(Synthetic Aperture Radar，SAR)地面运动目标成像方法。
背景技术
合成孔径雷达因具有全天时、全天候、远距离成像的能力，在战略防御、地形测绘等军用和民用领域得到了广泛应用。其中，利用SAR获取运动目标成像结果已成为现代雷达的研究热点之一。
同时实现高分辨宽测绘带成像是雷达成像技术的追求目标。然而，传统星载 SAR系统受最小天线面积限制，无法同时获取高分辨宽测绘带图像。结合数字波束形成技术，方位多通道SAR系统通过发射低脉冲重复频率(Pulse Repetition Frequency,PRF)信号可克服此限制，但是低PRF将导致回波信号产生严重的多普勒模糊。针对地面静止目标，现有多种方法可实现多普勒模糊抑制从而实现高分辨宽测绘带SAR成像。但是，对于地面运动目标来说，由于其运动参数未知，无法直接精确构造运动目标的模糊重构函数，从而导致处理性能下降。
针对以上运动目标成像存在的问题，各国学者提出了多种方法。Li等人提出通过搜索空时平面内的运动目标区，得到运动目标所有可能的谱分量，但并未给出相应的判断准则。Baumgartner等人根据最大信杂噪比准则搜索运动目标的速度，但由于该方法需要对所有可能的运动目标速度都进行一次全孔径成像操作，因此运算量会大大增加。Yang和Zhang等人假设运动目标谱分量位于杂波谱分量中间，尽量可能地正确约束运动目标谱分量，虽然该方法可以避免繁重的搜索操作，但当运动目标的运动方向偏离该导向方向时，该方法将无法正确对运动目标的谱分量进行约束，从而无法对模糊分量进行抑制并导致运动目标增益的损失。
发明内容
本发明针对上述高分辨宽测绘带SAR地面运动目标成像方法面临的问题，提出了一种基于稀疏空间谱估计的高分辨宽测绘带SAR地面运动目标成像方法。与传统方法不同，本发明没有采用搜索或近似假设，并能以较小的运算量实现对地面运动目标的高质量成像。
为实现上述目的，本发明的主要步骤如下：
(1)对各通道接收的回波信号进行方位向傅里叶变换；
(2)根据稀疏空间谱估计方法，利用某一多普勒单元回波估计得到运动目标的速度；
(3)根据估计得到的运动目标速度，构造运动目标的频谱重构函数，得到运动目标的无模糊全带宽信号；
(4)根据估计得到的运动目标速度，对无模糊的全带宽回波信号进行传统高分辨运动目标成像，得到地面运动目标高分辨率宽测绘带图像。
本发明与现有的技术相比具有以下优点：
本发明利用稀疏空间谱估计的方法直接对运动目标的速度进行精确估计，从而避免了传统搜索操作引入的巨大运算量，另外运动目标速度的精确估计在保证运动目标的多普勒频谱精确重构的同时，克服了传统方法带来的运动目标谱分量增益的损失。
附图说明
图1为本发明的设计方法流程图；
图2为在斜距平面下的SAR系统观测几何示意图；
图3为利用回波信号的某一多普勒单元进行空间谱估计的结果；
图4为利用传统方法得到的运动目标仿真成像结果；
图5为利用本发明方法的仿真成像结果。
具体实施方式
参照附图1，本发明的具体实施步骤如下：
步骤1，对各通道接收的回波信号进行方位向傅里叶变换。
假设SAR系统共有M个沿航向均匀分布的接收通道，如图2所示，x轴为平台运动速度方向，y轴为斜距方向，雷达到地面目标的斜距矢量与零多普勒面的夹角φ称为锥角。不失一般性，本发明假设SAR系统为正侧视成像。对于方位多通道高分辨宽测绘带SAR系统来说，各通道接收回波均可等效为参考接收通道(如第1个通道)接收回波的时延。为保证宽距离测绘带，SAR系统采用低PRF，这样回波信号将发生多普勒模糊。考虑加性白噪声，第m个通道接收的回波信号进行方位向傅里叶变换后可写为
其中，m＝1,2,…,M，M为方位接收通道个数，τ为距离时间，fp为系统PRF， fd∈[-fp/2,fp/2]为多普勒频率，vs为平台运动速度，a(τ,fd)为第1个通道的接收回波，Δxm为第m个通道相对第1个通道的距离，L＝(N-1)/2，N是回波信号的多普勒模糊数，zm(τ,fd)为第m个通道的噪声。地面运动目标回波的多普勒频率与锥角φ的关系如下
其中，λ为载波频率，vt为运动目标的垂直航向速度，
忽略τ，式（ 1 ） 可用矢量形式表示为
其中，
s(τ,fd)＝[s1(τ,fd),s2(τ,fd),…,sM(τ,fd)]T， (4)
a(τ,fd)＝[a(τ,fd-L·fp),a(τ,fd-(L-1)·fp),…,a(τ,fd+L·fp)]T， (5)
P(fd)＝[p-L(fd),p-L+1(fd),…,pL(fd)]， (6)
z(τ,fd)＝[z1(τ,fd),z2(τ,fd),…,zM(τ,fd)]T， (8)
上标T表示矩阵转置。由此可知，各通道接收的多普勒模糊回波信号可等效看作为来自不同方向的信号叠加。与地面固定杂波不同，由于运动目标的速度vt未知，所以无法精确获知阵列流形矩阵P(fd)，从而无法直接利用已有方法对回波信号进行频谱重构。
步骤2，根据稀疏空间谱估计方法，利用某个多普勒单元回波信号估计得到运动目标的速度。
根据式(2)所示的多普勒频率与锥角的关系，本发明将运动目标速度估计问题转换为信号的波达方向(direction-of-arrival,DOA)估计问题，也即空间谱估计问题。根据SAR系统参数可计算出多普勒模糊数N，也即待估计DOA信号的个数。更进一步，由于运动目标空间谱在空间域上为稀疏分布，因此上述DOA 估计问题可转换为以下稀疏问题
其中，
Y＝[s(τ1,fd),…,s(τK,fd)]， (10)
Y为各通道接收的快拍信号，K为快拍个数，在实际中，可以利用相邻的K个距离单元回波替代。为所有可能方向的回波信号构成的阵列流形，φq (q＝1,…,Q)的取值范围为雷达波束的照射范围。为Q×1维稀疏向量，若第n(n＝1,…,N)个信号的方向角为φq，则的第q个值为 a(τ,fd-ln·fp)，其中ln＝n-L-1，否则其值为零。一般来说，Q＞＞M＞N。由此可将代价函数写为
其中，||·||f为Frobenius范数，下标 q表示向量的第q个值，β为稀疏规则化参数。求解上述优化问题即可得到然后根据波峰位置得到信号的波达方向，再根据下式计算出运动目标的速度
其中，fd为所选取的多普勒单元，φq为估计得到的波达方向。在实际操作中，为了减小噪声的影响，可以选择多个多普勒单元估计得到多个vt，然后将其平均求得最终的vt。
式(15)所示的优化问题为凸优化问题，现已有多种成熟的求解方法。需要说明的是，对于SAR系统来说，M和N的值都较小，而且在实际操作中，仅需要对式(15)进行若干几次求解，因此相比传统的搜索方法，本发明的运算量将大大降低。
步骤3，根据上一步估计得到的运动目标速度，构造运动目标的多普勒频谱重构函数，得到运动目标的无模糊全带宽信号。
利用步骤2估计得到的运动目标速度，即可正确构造出运动目标回波信号的阵列流形P(fd)，然后可采用传统的模糊重构方法得到无模糊的运动目标全带宽信号。
步骤4，根据估计得到的运动目标速度，对无模糊的全带宽回波信号进行传统高分辨运动目标成像，得到地面运动目标高分辨率宽测绘带图像。
利用步骤2估计得到的运动目标速度，对回波信号进行距离压缩和距离徙动校正，然后再进行方位压缩，得到地面运动目标的高分辨宽测绘带SAR图像。
下面结合仿真数据实验对本发明的效果做进一步的说明。
1、仿真条件：
星载方位多通道SAR系统的仿真参数如下表所示。多普勒模糊次数约为3，由于运动目标存在垂直航向速度，其多普勒频谱将偏离杂波谱约682Hz。
2、仿真数据包实验分析：
首先对各通道回波信号进行方位向傅里叶变换，然后选取某一多普勒单元回波进行稀疏空间谱估计。图3给出了利用第1360个多普勒单元回波进行空间谱估计的结果，相邻的10个距离单元作为快拍数。可以看出，利用本发明的方法可以准确估计得到运动目标的空间谱，由此估计得到的目标垂直航向速度为 10.53m/s，对应的多普勒偏移约为677Hz，与真实值非常接近。图4给出了假设运动目标谱偏离杂波谱fp/2所构造的频谱重构函数得到的运动目标成像结果的方位剖面图，可以看出成像结果存在较为严重的方位模糊。图5(a)给出了利用本发明所述的方法得到的运动目标成像结果，图5(b)为成像结果的局部放大图，可以看出运动目标的多普勒模糊分量得到了有效的抑制，并实现了很好的聚焦。

Claims (2)

1.高分辨宽测绘带合成孔径雷达地面运动目标成像方法，包括如下步骤：
(1)对各通道接收的回波信号进行方位向傅里叶变换；
(2)求解以下稀疏空间谱估计问题
其中，||·||f为Frobenius范数，||·||1为l1范数，Y为各通道接收的快拍信号，为所有可能方向的回波信号构成的阵列流形，为目标回波信号构成的稀疏矩阵，为各通道接收回波的快拍信号向量的l2范数组成的回波稀疏向量，β为稀疏规则化参数，求解上述稀疏空间谱估计问题即可得到然后根据波峰位置得到信号的波达方向，再根据下式计算出运动目标的速度
其中，λ为载波频率，fd为多普勒频率，l∈[-L,L]，L＝(N-1)/2，N是回波信号的多普勒模糊数，fp为系统PRF，vs为平台运动速度，φq为估计得到的波达方向，在实际操作中，为了减小噪声的影响，可以选择多个多普勒单元估计得到多个vt，然后将其平均求得最终的vt；
(3)根据估计得到的运动目标速度，构造运动目标的频谱重构函数，得到运动目标的无模糊全带宽信号；
(4)根据估计得到的运动目标速度，对无模糊的全带宽信号进行传统高分辨运动目标成像，得到高分辨宽测绘带运动目标SAR图像。
2.根据权利要求1所述的高分辨宽测绘带合成孔径雷达地面运动目标成像方法，其特征在于：各通道接收的多普勒模糊回波信号可等效看作为来自不同方向的信号叠加，对各通道回波进行方位向傅里叶变换后可用矢量形式表示为
s(τ,fd)＝P(fd)a(τ,fd)+z(τ,fd)
其中，τ为距离时间，fd为多普勒频率，s(τ,fd)为各通道接收回波构成的向量，P(fd)为各通道回波构成的阵列流形矩阵，a(τ,fd)为第1个通道接收的多普勒模糊回波信号构成的向量，z(τ,fd)为各通道噪声构成的向量。
CN201510726789.3A 2015-10-30 2015-10-30 高分辨宽测绘带sar地面运动目标成像方法 Expired - Fee Related CN105301589B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510726789.3A CN105301589B (zh) 2015-10-30 2015-10-30 高分辨宽测绘带sar地面运动目标成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510726789.3A CN105301589B (zh) 2015-10-30 2015-10-30 高分辨宽测绘带sar地面运动目标成像方法

Publications (2)

Publication Number Publication Date
CN105301589A CN105301589A (zh) 2016-02-03
CN105301589B true CN105301589B (zh) 2018-10-23

Family

ID=55199068

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510726789.3A Expired - Fee Related CN105301589B (zh) 2015-10-30 2015-10-30 高分辨宽测绘带sar地面运动目标成像方法

Country Status (1)

Country Link
CN (1) CN105301589B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105785366B (zh) * 2016-04-13 2019-02-26 电子科技大学 基于稀疏空间谱估计的星载sar数字波束形成方法
US10539672B2 (en) * 2016-12-09 2020-01-21 GM Global Technology Operations LLC Doppler ambiguity resolution at high signal to noise ratio
CN106872954B (zh) * 2017-01-23 2019-07-09 西安电子科技大学 一种高超声速平台杂波抑制和运动目标成像方法
CN107607952B (zh) * 2017-10-16 2020-08-11 电子科技大学 基于电磁涡旋波的三维合成孔径雷达成像方法
CN111948653B (zh) * 2020-07-31 2022-03-18 上海卫星工程研究所 基于p波段合成孔径雷达的林下目标检测方法及系统
CN112946614A (zh) * 2021-03-22 2021-06-11 哈尔滨工程大学 一种基于qga的多普勒频移目标运动参数的解算方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100523865C (zh) * 2005-12-28 2009-08-05 中国科学院电子学研究所 确定距离向多孔径接收宽测绘带合成孔径雷达子天线间距的方法
US7397418B1 (en) * 2006-06-05 2008-07-08 Sandia Corporation SAR image formation with azimuth interpolation after azimuth transform
DE102012021010B4 (de) * 2012-10-26 2022-02-03 Airbus Defence and Space GmbH Synthetisches Apertur Radar zur gleichzeitigen Bildaufnahme und Bewegtzielerkennung
CN104865568B (zh) * 2015-06-02 2017-05-24 西安电子科技大学 基于稀疏重构的宽带雷达高速群目标分辨方法

Also Published As

Publication number Publication date
CN105301589A (zh) 2016-02-03

Similar Documents

Publication Publication Date Title
CN105301589B (zh) 高分辨宽测绘带sar地面运动目标成像方法
CN105785327B (zh) 频率分集阵列合成孔径雷达高分辨宽测绘带成像方法
CN111856461B (zh) 基于改进pfa的聚束sar成像方法及其dsp实现
CN107229048B (zh) 一种高分宽幅sar动目标速度估计与成像方法
CN109143237B (zh) 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法
CN108693511B (zh) 时分复用mimo雷达的运动目标角度计算方法
CN112444811A (zh) 一种综合mimo雷达和isar的目标检测和成像方法
WO2021062914A1 (zh) 一种基于mimo体制的提高雷达角度分辨率的布局及方法
CN111352107B (zh) 基于多通道数字和差的单脉冲跟踪与成像方法
CN110082764B (zh) 基于稳健正则化层析方法的sar图像成像方法
CN105785366B (zh) 基于稀疏空间谱估计的星载sar数字波束形成方法
WO2024045362A1 (zh) 车载多普勒分集mimo合成孔径雷达时域成像方法
CN105842693A (zh) 一种基于压缩感知的双通道sar动目标检测的方法
CN106019280B (zh) 基于距离多普勒校正的fmcw sar成像方法及装置
CN111007503A (zh) 基于频率谱精确定位的运动目标聚焦和定位方法及系统
CN105929397B (zh) 基于正则化的偏置相位中心天线成像方法
CN110133646B (zh) 基于nlcs成像的双基前视sar的多通道两脉冲杂波对消方法
CN110261852B (zh) 一种频率分集逆合成孔径雷达成像系统及方法
CN116184347A (zh) 基于距离选通与子空间投影的机载雷达折叠杂波抑制方法
CN109143236B (zh) 适用于复杂飞行轨迹的双基聚束sar大场景成像方法
CN110658502A (zh) 一种幅相误差校正方法
Smith et al. Multipath exploitation for radar target classification
CN111638516B (zh) 基于双频共轭处理技术的太赫兹频段sar运动补偿算法
CN114994676A (zh) 一种一站固定式双站低频超宽带sar运动目标成像方法
CN110780273B (zh) 一种混合正则化方位超分辨成像方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181023

Termination date: 20201030

CF01 Termination of patent right due to non-payment of annual fee