CN110376564A - Geo和leo双基构型合成雷达地面运动成像方法 - Google Patents

Geo和leo双基构型合成雷达地面运动成像方法 Download PDF

Info

Publication number
CN110376564A
CN110376564A CN201910692175.6A CN201910692175A CN110376564A CN 110376564 A CN110376564 A CN 110376564A CN 201910692175 A CN201910692175 A CN 201910692175A CN 110376564 A CN110376564 A CN 110376564A
Authority
CN
China
Prior art keywords
matrix
orbit satellite
radar
nrn
low
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.)
Granted
Application number
CN201910692175.6A
Other languages
English (en)
Other versions
CN110376564B (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.)
Northwest University of Technology
Original Assignee
Northwest University of Technology
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 Northwest University of Technology filed Critical Northwest University of Technology
Priority to CN201910692175.6A priority Critical patent/CN110376564B/zh
Publication of CN110376564A publication Critical patent/CN110376564A/zh
Application granted granted Critical
Publication of CN110376564B publication Critical patent/CN110376564B/zh
Active 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/415—Identification of targets based on measurements of movement associated with the target

Landscapes

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

Abstract

本发明涉及一种GEO和LEO双基构型合成雷达地面运动成像方法，通过依次构造参考信号向量、距离走动校正矩阵、Chirp傅里叶变换矩阵、并对其进行一系列的计算，得到目标成像结果图。本发明克服了已有SAR成像技术无法对地球同步轨道和低轨道卫星双基构型合成雷达地面非合作运动目标进行成像将出现散焦现象，可以对地球同步轨道星和低轨道卫星双基构型合成雷达非合作目标走动进行自适应校正，以达到对非合作目标进行良好聚焦的效果。

Description

GEO和LEO双基构型合成雷达地面运动成像方法
技术领域
本发明属于雷达信号处理领域，特别涉及一种地球同步轨道卫星(GEO)和低轨道卫星(LEO)双基构型合成雷达(Bi-SAR)地面运动目标成像算法，适用于地球同步轨道卫星和低轨道卫星双基合成孔径雷达对地面运动目标进行成像处理。
背景技术
双基地合成孔径雷达(Bistatic Synthetic Aperture Radar,Bi-SAR)因其收发天线分置，接收雷达不易被侦察发现，具有较好的隐身和抗干扰性能；并且通过合理配置收发天线几何构型，系统在自然目标和人造目标观测方面具有极大的自由度，有利于SAR图像解译，具有非常重要的军事和民事应用价值。地球同步轨道(GEO)发射-近地轨道(LEO)接收体制下的Bi-SAR系统更是因其战场生存能力强以及便捷、频繁、大范围对地观测能力得到了国内外的关注与研究。在地球同步轨道(GEO)发射-近地轨道(LEO)接收体制下研究地面运动目标探测技术也具有非常重要的军事和民事应用价值。
对于地面运动目标探测，面临的关键问题是地面运动成像聚焦问题。由于地球同步轨道和低轨道卫星双基构型合成雷达地面运动目标存在非合作特性，直接利用已有SAR成像技术，将无法对地面运动目标进行有效成像聚焦。
发明内容
要解决的技术问题
针对已有SAR成像技术无法对地球同步轨道和低轨道卫星双基构型合成雷达地面非合作运动目标进行聚焦成像，本发明的目的在于提出一种地球同步轨道星和低轨道卫星双基构型合成雷达地面运动成像算法，该方法可以有效对非合作运动目标进行成像聚焦。
技术方案
一种地球同步轨道卫星和低轨道卫星双基构型合成雷达地面运动成像方法，其特征在于步骤如下：
步骤1：获取的杂波抑制之后雷达系统回波信号是二维矩阵S(i,j)，S(i,j)为nrn×nan维矩阵，对矩阵S(i,j)列进行FFT处理，将结果分别保存在矩阵S(i,j)中；
步骤2：根据已知的雷达参数，构造参考信号向量S_ref(i,1)为nrn×1向量；其中，γ表示调频率，γ＝B/Tp，B表示发射信号带宽，Tp表示发射脉冲宽度，fm表示为距离向频域坐标，Δf为距离频域间隔， nrn表示距离向点数；
步骤3：取出步骤1得到的S(i,j)的每一列，均点乘参考信号向量S_ref(i,1)的共轭，得到距离脉压后的数据矩阵S(fm,xn)；其中，xn表示方位向时域坐标，L表示为合成孔径长度，n＝0,1,...,nan-1，nan表示方位向点数；
步骤4：根据已知的雷达参数，构造距离走动校正矩阵S_LRCMC(i,j)为nrn×nan矩阵；其中，fc表示雷达发射信号的载频，c为电磁波的传播速度，tm为方位慢时间，为多普勒中心；
步骤5：取出步骤3得到的S(fm,xn)点乘参考信号向量S_LRCMC(i,j)的共轭，得到走动校正之后的数据矩阵S′(fm,xn)；
步骤6：根据已知的雷达参数，构造Chirp傅里叶变换矩阵ST_chirp(i,j)为nrn×nan矩阵；其中Vr为低轨道卫星速度，θr为低轨道卫星的斜视角，VT为地球同步轨道卫星速度，θT为低轨道卫星的斜视角，vT,a为地面运动目标相对发射天线的径向速度，RT为地面运动目标相对发射天线的场景中心斜距，va,r为地面运动目标相对接收天线的径向速度，Rr为地面运动目标相对接收天线的场景中心斜距；
步骤7：取出步骤5得到的S′(fm,xn)点乘参考信号向量ST_chirp(i,j)的共轭，得到Chirp傅里叶变换之后的数据矩阵S″(fm,xn)；
步骤8：取出步骤7得到的S″(fm,xn)，沿方位慢时间进行插值处理，处理之后得到Sout(fm,x′n)，同时对矩阵Sout(fm,x′n)行进行FFT处理，将结果分别保存在矩阵Simage(fm,x′n)中，即为获得的聚焦地面运动目标成像结果图。
有益效果
本发明提出的一种地球同步轨道卫星和低轨道卫星双基构型合成雷达地面运动成像方法，克服了已有SAR成像技术无法对地球同步轨道和低轨道卫星双基构型合成雷达地面非合作运动目标进行成像将出现散焦现象，可以对地球同步轨道星和低轨道卫星双基构型合成雷达非合作目标走动进行自适应校正，以达到对非合作目标进行良好聚焦的效果。
附图说明
图1是本发明的仿真数据处理结果图：(a)距离走动校正之后的地面运动目标2维频谱信号；(b)距离压缩之后的地面运动目标信号；(c)经常chirp傅里叶变换之后的运动目标回波信号；(d)经过keystone变换聚焦的运动目标成像结果；(e)进行插值之后的运动目标成像结果。
具体实施方式
现结合实施例、附图对本发明作进一步描述：
一种地球同步轨道星和低轨道卫星双基构型合成雷达地面运动成像算法，具体包括以下步骤：
步骤1，获取的杂波抑制之后雷达系统回波信号是二维矩阵S(i,j)，S(i,j)为nrn×nan维矩阵，对矩阵S(i,j)列进行FFT处理，将结果分别保存在矩阵S(i,j)中；
步骤2，根据已知的雷达参数，构造参考信号向量S_ref(i,1)，S_ref(i,1)为nrn×1向量；
其中，nrn表示距离向点数；
步骤3，取出步骤1得到的S(i,j)的每一列，均点乘参考信号向量S_ref(i,1)的共轭，得到距离脉压后的数据矩阵S(fn,xm)；
其中，fm表示为距离向频域坐标，B为发射信号带宽，Δf为距离频域间隔，nrn表示距离向点数，xm表示方位向时域坐标，L表示为合成孔径长度，m＝0,1,...,nan-1，nan表示方位向点数；
步骤4，根据已知的雷达参数，构造距离走动校正矩阵S_LRCMC(i,j)，S_LRCMC(i,j)为nrn×nan矩阵；
其中，nrn表示距离向点数，nan表示方位向点数；
步骤5，取出步骤3得到的S(fn,xm)点乘参考信号向量S_LRCMC(i,j)的共轭，得到走动校正之后的数据矩阵S′(fn,xm)；
步骤6，根据已知的雷达参数，构造Chirp傅里叶变换矩阵ST_chirp(i,j)，ST_chirp(i,j)为nrn×nan矩阵；
步骤7，取出步骤5得到的S′(fn,xm)点乘参考信号向量ST_chirp(i,j)的共轭，得到Chirp傅里叶变换之后的数据矩阵S″(fn,xm)；
步骤8，取出步骤7得到的S″(fn,xm)，沿方位慢时间进行插值处理，处理之后得到Sout(fn,x′m)，同时对矩阵Sout(fn,x′m)行进行FFT处理，将结果分别保存在矩阵Simage(fn,x′m)中，即为获得的聚焦地面运动目标成像结果图；
上述技术方案的特点和进一步改进在于：
步骤2具体操作如下：
根据已知的雷达参数，参考信号向量
其中，γ表示调频率，γ＝B/Tp，B表示发射信号带宽，Tp表示发射脉冲宽度，fn表示为距离向频域坐标，B为发射信号带宽，Δf为距离频域间隔，nrn表示距离向点数。
步骤4具体操作如下：
根据已知的雷达参数，构造距离走动校正矩阵为
其中，fc表示雷达发射信号的载频，c为电磁波的传播速度，B为发射信号带宽，Δf为距离频域间隔，nrn表示距离向点数，tm为方位慢时间，为多普勒中心。
步骤6具体操作如下：
根据已知的雷达参数，构造Chirp傅里叶变换矩阵ST_chirp(i,j)为
其中fc表示雷达发射信号的载频，c为电磁波的传播速度，Vr为低轨道卫星速度，θr为低轨道卫星的斜视角，VT为地球同步轨道卫星速度，θT为低轨道卫星的斜视角，vT,a为地面运动目标相对发射天线的径向速度，RT为地面运动目标相对发射天线的场景中心斜距，va,r为地面运动目标相对接收天线的径向速度，Rr为地面运动目标相对接收天线的场景中心斜距，B为发射信号带宽，Δf为距离频域间隔，nrn表示距离向点数，tm为方位慢时间。
至此，一种地球同步轨道卫星和低轨道卫星双基构型合成雷达地面运动目标成像算法基本完成。
以下通过仿真实验数据来进一步验证本发明的有效性。
(一)仿真实验
1.仿真参数
为了验证本发明方法的有效性，此处给出了表1中的仿真参数。
表1仿真数据参数
载频 9.7GHz 通道数 2
同步轨道卫星速度 855m/s 低轨道卫星速度 2164km
同步轨道卫星高度 4915.4km 低轨道卫星高度 117.97km
同步轨道卫星斜距 36122km 低轨道卫星斜距 2164km
同步轨道卫星斜视角 0.2° 同步轨道卫星斜视角 0.1°
脉冲重复频率 2400Hz 发射信号带宽 250MHz
2.仿真内容
图1示意了利用本发明提出的地球同步轨道卫星和低轨道卫星双基构型合成雷达地面运动目标数据处理结果。从图中可以看出本发明方法的地面运动目标成像聚焦效果，采用本发明方法的可以有效解决双基构型情况下地面运动目标成像问题。
综上所述，仿真实验验证了本发明的正确性，有效性和可靠性。

Claims (1)

1.一种地球同步轨道卫星和低轨道卫星双基构型合成雷达地面运动成像方法，其特征在于步骤如下：
步骤1：获取的杂波抑制之后雷达系统回波信号是二维矩阵S(i,j)，S(i,j)为nrn×nan维矩阵，对矩阵S(i,j)列进行FFT处理，将结果分别保存在矩阵S(i,j)中；
步骤2：根据已知的雷达参数，构造参考信号向量S_ref(i,1)为nrn×1向量；其中，γ表示调频率，γ＝B/Tp，B表示发射信号带宽，Tp表示发射脉冲宽度，fm表示为距离向频域坐标，Δf为距离频域间隔，m＝0,1,...,nrn-1，nrn表示距离向点数；
步骤3：取出步骤1得到的S(i,j)的每一列，均点乘参考信号向量S_ref(i,1)的共轭，得到距离脉压后的数据矩阵S(fm,xn)；其中，xn表示方位向时域坐标，L表示为合成孔径长度，n＝0,1,...,nan-1，nan表示方位向点数；
步骤4：根据已知的雷达参数，构造距离走动校正矩阵S_LRCMC(i,j)为nrn×nan矩阵；其中，fc表示雷达发射信号的载频，c为电磁波的传播速度，tm为方位慢时间，为多普勒中心；
步骤5：取出步骤3得到的S(fm,xn)点乘参考信号向量S_LRCMC(i,j)的共轭，得到走动校正之后的数据矩阵S′(fm,xn)；
步骤6：根据已知的雷达参数，构造Chirp傅里叶变换矩阵ST_chirp(i,j)为nrn×nan矩阵；其中Vr为低轨道卫星速度，θr为低轨道卫星的斜视角，VT为地球同步轨道卫星速度，θT为低轨道卫星的斜视角，vT,a为地面运动目标相对发射天线的径向速度，RT为地面运动目标相对发射天线的场景中心斜距，va,r为地面运动目标相对接收天线的径向速度，Rr为地面运动目标相对接收天线的场景中心斜距；
步骤7：取出步骤5得到的S′(fm,xn)点乘参考信号向量ST_chirp(i,j)的共轭，得到Chirp傅里叶变换之后的数据矩阵S″(fm,xn)；
步骤8：取出步骤7得到的S″(fm,xn)，沿方位慢时间进行插值处理，处理之后得到Sout(fm,x′n)，同时对矩阵Sout(fm,x′n)行进行FFT处理，将结果分别保存在矩阵Simage(fm,x′n)中，即为获得的聚焦地面运动目标成像结果图。
CN201910692175.6A 2019-07-30 2019-07-30 Geo和leo双基构型合成雷达地面运动成像方法 Active CN110376564B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910692175.6A CN110376564B (zh) 2019-07-30 2019-07-30 Geo和leo双基构型合成雷达地面运动成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910692175.6A CN110376564B (zh) 2019-07-30 2019-07-30 Geo和leo双基构型合成雷达地面运动成像方法

Publications (2)

Publication Number Publication Date
CN110376564A true CN110376564A (zh) 2019-10-25
CN110376564B CN110376564B (zh) 2022-11-08

Family

ID=68256843

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910692175.6A Active CN110376564B (zh) 2019-07-30 2019-07-30 Geo和leo双基构型合成雷达地面运动成像方法

Country Status (1)

Country Link
CN (1) CN110376564B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113671503A (zh) * 2021-08-23 2021-11-19 西北工业大学 一种适用于星载变prf sar系统的成像方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103135100A (zh) * 2013-01-31 2013-06-05 西安电子科技大学 同轨双基sar的动目标参数估计方法
CN106249229A (zh) * 2016-07-20 2016-12-21 西安电子工程研究所 基于PD雷达体制的Keystone变换空域及时域走动补偿方法
CN106842200A (zh) * 2017-01-11 2017-06-13 中国科学院电子学研究所 一种双基合成孔径雷达成像方法和装置
CN108469612A (zh) * 2018-03-20 2018-08-31 西安电子科技大学 基于等效斜距的双基时变加速度前视sar成像方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103135100A (zh) * 2013-01-31 2013-06-05 西安电子科技大学 同轨双基sar的动目标参数估计方法
CN106249229A (zh) * 2016-07-20 2016-12-21 西安电子工程研究所 基于PD雷达体制的Keystone变换空域及时域走动补偿方法
CN106842200A (zh) * 2017-01-11 2017-06-13 中国科学院电子学研究所 一种双基合成孔径雷达成像方法和装置
CN108469612A (zh) * 2018-03-20 2018-08-31 西安电子科技大学 基于等效斜距的双基时变加速度前视sar成像方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YONGKANG LI ET AL.: "A New Motion Parameter Estimation and Relocation Scheme for Airborne Three-Channel CSSAR-GMTI Systems", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》, vol. 57, no. 6, 30 June 2019 (2019-06-30), pages 4107 - 4119, XP011725945, DOI: 10.1109/TGRS.2019.2894620 *
吴玉峰 等: "周期性非均匀采样实现星载SAR高分辨宽测绘带成像", 《电子与信息学报》, vol. 34, no. 2, 26 February 2014 (2014-02-26), pages 279 - 285 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113671503A (zh) * 2021-08-23 2021-11-19 西北工业大学 一种适用于星载变prf sar系统的成像方法
CN113671503B (zh) * 2021-08-23 2023-09-08 西北工业大学 一种适用于星载变prf sar系统的成像方法

Also Published As

Publication number Publication date
CN110376564B (zh) 2022-11-08

Similar Documents

Publication Publication Date Title
CN105785327B (zh) 频率分集阵列合成孔径雷达高分辨宽测绘带成像方法
Wang et al. MIMO SAR using chirp diverse waveform for wide-swath remote sensing
CN106093870B (zh) 高超声速飞行器下降段的sar-gmti杂波抑制方法
CN110109117B (zh) 基于频控阵的星载合成孔径雷达卷积调制干扰方法
CN105445704B (zh) 一种sar图像中的雷达动目标抑制方法
CN111443339A (zh) 一种双基sar空变校正成像方法、装置、设备及存储介质
CN102313887B (zh) 一种星机联合双基地合成孔径雷达成像方法
CN109143237B (zh) 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法
CN104833972B (zh) 一种双基地调频连续波合成孔径雷达频率变标成像方法
CN110221296B (zh) 基于模拟转发模式的高-低轨sar卫星联合观测系统
CN106291489B (zh) 适用于多种发射信号波形的合成孔径雷达回波仿真方法
CN106950565A (zh) 星载sar成像抖动补偿方法、成像方法
Qiu et al. Bistatic SAR data processing algorithms
CN112285658A (zh) 基于递推计算的星载sar欺骗干扰信号快速生成方法
CN109143236B (zh) 适用于复杂飞行轨迹的双基聚束sar大场景成像方法
CN110376564A (zh) Geo和leo双基构型合成雷达地面运动成像方法
CN106526544B (zh) 基于高超声速平台的mimosar杂波抑制方法
CN113567978A (zh) 一种多基分布式雷达协同成像方法
CN113189547B (zh) 一种基于合成带宽频率变标的sar成像的方法及系统
CN110488236B (zh) Geo和leo双基构型合成雷达地面杂波抑制方法
Garry et al. Framework and results on passive ISAR
CN109490887A (zh) 基于外辐射源的雷达成像算法
CN114646958A (zh) 一种分布式小卫星聚束mimo-sar超高分辨成像方法
CN117872290B (zh) 一种基于相位编码超表面的雷达多维特征调控方法及系统
Zhou et al. A novel algorithm for MIMO SAR imaging

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant