CN110824473A - 一种基于子空间的高分辨宽测绘带sar-gmti杂波抑制方法 - Google Patents
一种基于子空间的高分辨宽测绘带sar-gmti杂波抑制方法 Download PDFInfo
- Publication number
- CN110824473A CN110824473A CN201911000657.7A CN201911000657A CN110824473A CN 110824473 A CN110824473 A CN 110824473A CN 201911000657 A CN201911000657 A CN 201911000657A CN 110824473 A CN110824473 A CN 110824473A
- Authority
- CN
- China
- Prior art keywords
- channel
- vector
- reference signal
- clutter
- matrix
- 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
Links
Images
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
- 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/418—Theoretical aspects
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‑GMTI杂波抑制方法,对获取的多通道雷达系统的第k个通道回波信号进行FFT处理,取每一行点乘参考信号向量的共轭,再点乘通道参考信号向量的共轭,再点乘参考信号向量的共轭,最后进行IFFT处理,构造多通道系统回波矢量,取出矢量对协方差矩阵进行估计,对协方差矩阵进行特征值分解得到杂波对应的零空间,对零空间进行杂波抑制。本发明解决了星载沿方位多通道高分辨宽测绘带SAR‑GMTI系统存在通道误差情况下进行杂波抑制问题,实现最大程度地消除地物场景杂波对运动目标检测影响,以达到对地面及海面运动目标进行广域高分辨探测的目的。
Description
技术领域
本发明属于雷达信号处理领域,特别涉及一种基于子空间的星载多通道高分辨宽测绘带SAR-GMTI杂波抑制算法,适用于星载方位多通道高分辨宽测绘带SAR-GMTI雷达在通道误差存在情况下对地面杂波进行抑制。
背景技术
星载方位多通道高分辨宽测绘带合成孔径雷达地面运动目标检测(SyntheticAperture Radar-ground moving target indication,SAR-GMTI)系统可以有效地解决距离宽测绘带方位高分辨率之间矛盾,进而达到对地面场景进行高分辨宽测绘带观测。该系统利用低重复频率发射脉冲信号,可以获取宽的场景测绘带。为了避免脉冲低重复频率引起方位低采样率,进而带来的方位多普勒频谱模糊问题,利用沿方位向的多通道系统对回波沿方位向进行空域采样,并对多普勒频谱进行无模糊重构。本质上,该系统利用方位空域采样来替代时间采样进而获取无模糊的多普勒频谱。该系统具有冗余的通道自由度,具备对地面杂波及海洋杂波进行有效抑制能力,进而达到对地面运动目标、海上军舰及大批军事调动进行探测。同时该系统可以对敏感和高危目标进行检测、跟踪和定位,并实现引导预警机和战斗机对目标进行有效监视及精确打击。
对于地面运动及海面目标进行探测时,首先面临的问题是杂波抑制问题。由于星载多通道高分辨宽测绘带SAR-GMTI系统不可避免存在通道误差问题,直接利用已有机载SAR-GMTI杂波抑制技术,如相位中心偏置、空时自适应处理等,将无法有效的对地杂波及海杂波进行抑制。
发明内容
要解决的技术问题
针对现有杂波抑制技术无法有效解决通道误差存在情况星载多通道高分辨宽测绘带SAR-GMTI杂波抑制的不足,本发明提出一种基于子空间的星载多通道高分辨宽测绘带合成雷达地面运动目标杂波抑制方法,该方法可以在相位误差及天线位置测量误差存在情况对杂波进行有效抑制。
技术方案
一种基于子空间的高分辨宽测绘带SAR-GMTI杂波抑制方法,其特征在于步骤如下:
步骤1:获取的多通道雷达系统的第k个通道回波信号是二维矩阵Sk(n,m),所述的Sk(n,m)为nrn×nan维矩阵,k=1,2,…,K,多通道系统共有K个通道;对矩阵Sk(n,m)列进行FFT处理,将结果保存在矩阵Sk(n,m)中;其中,nrn表示距离向点数,nan表示方位向点数;
步骤2:根据已知的雷达参数,构造参考信号向量S_ref(n)为nrn×1向量;其中,γ表示调频率,γ=B/Tp,B表示发射信号带宽,Tp表示发射脉冲宽度,fn表示为距离向频域坐标,B为发射信号带宽,Δf为距离频域间隔,n=0,1,...,nrn-1;
步骤3:取出步骤1得到Sk(n,m)的每一列,均点乘参考信号向量S_ref(n)的共轭,得到距离脉压后的第k通道数据矩阵Sk(fn,xm);xm表示方位向时域坐标,L表示为合成孔径长度,m=0,1,...,nan-1;
步骤4:根据已知的雷达参数,构造第k个通道参考信号向量S_refk(m)为1×nan向量;其中,fc表示雷达发射信号的载频,c为电磁波的传播速度,V为卫星速度,θ为卫星的斜视角,dk为第k通道与参考通道沿卫星速度方向的距离,tm为方位慢时间;
步骤5:取出步骤3得到的Sk(fn,xm)的每一行,均点乘参考信号向量S_refk(m)的共轭,得到距离走动校正之后的第k通道数据矩阵Sk′(fn,xm);
步骤7:取出步骤5得到的S′k(fn,xm)点乘参考信号向量S_Chirpk(n,m)的共轭,并对相乘结果行进行FFT处理,得到chirp傅里叶变换之后的第k通道数据矩阵S″k(fn,fm);其中,fm表示为方位向频域坐标,PRF为方位采用频率,Δfa为方位频域间隔,m=0,1,...,nan-1;
步骤8:取出步骤7得到S″k(fn,fm),沿fn方向进行IFFT处理,得到S″k(n,fm),同时构造多通道系统回波矢量S(n,fm)=[S″1(n,fm) S″2(n,fm) … S″k(n,fm) … S″K(n,fm)]T;其中[·]T表示矩阵转置;
有益效果
本发明提出的一种基于子空间的高分辨宽测绘带SAR-GMTI杂波抑制方法,解决星载沿方位多通道高分辨宽测绘带SAR-GMTI系统存在通道误差情况下进行杂波抑制问题,实现最大程度地消除地物场景杂波对运动目标检测影响,以达到对地面及海面运动目标进行广域高分辨探测的目的。
附图说明
图1实测数据处理结果:(a)录取的在距离脉压方位chirp傅里叶变换域单通道回波数据;(b)利用本发明进行杂波抑制之后在距离脉压方位chirp傅里叶变换域运动目标信号;(c)对Line A位置所示的信号在杂波抑制前和杂波抑制之后的信号进行比较;(d)对Line B位置所示的信号在杂波抑制前和杂波抑制之后的信号进行比较;(e)对Line B位置所示的信号在杂波抑制前和杂波抑制之后的信号进行比较。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
一种基于子空间的星载多通道高分辨宽测绘带SAR-GMTI杂波抑制算法,包括以下步骤:
步骤1,获取的多通道雷达系统的第k个通道回波信号是二维矩阵Sk(n,m),Sk(n,m)为nrn×nan维矩阵及k=1,2,…,K,多通道系统共有K个通道。对矩阵Sk(n,m)列进行FFT处理,将结果保存在矩阵Sk(n,m)中;
其中,nrn表示距离向点数,nan表示方位向点数;
步骤2,根据已知的雷达参数,构造参考信号向量S_ref(n),S_ref(n)为nrn×1向量;
其中,nrn表示距离向点数;
步骤3,取出步骤1得到Sk(n,m)的每一列,均点乘参考信号向量S_ref(n)的共轭,得到距离脉压后的第k通道数据矩阵Sk(fn,xm);
其中,fn表示为距离向频域坐标,B为发射信号带宽,Δf为距离频域间隔,n=0,1,...,nrn-1,nrn表示距离向点数,xm表示方位向时域坐标,L表示为合成孔径长度,m=0,1,...,nan-1,nan表示方位向点数;
步骤4,根据已知的雷达参数,构造第k个通道参考信号向量S_refk(m),S_refk(m)为1×nan向量;
其中,nan表示方位向点数;
步骤5,取出步骤3得到的Sk(fn,xm)的每一行,均点乘参考信号向量S_refk(m)的共轭,得到距离走动校正之后的第k通道数据矩阵S′k(fn,xm);
步骤6,根据已知的雷达参数,构造第k个通道chirp傅里叶变换参考信号S_Chirpk(n,m);
其中,nrn表示距离向点数,nan表示方位向点数;
步骤7,取出步骤5得到的S′k(fn,xm)点乘参考信号向量S_Chirpk(n,m)的共轭,并对相乘结果行进行FFT处理,得到chirp傅里叶变换之后的第k通道数据矩阵S″k(fn,fm);
步骤8,取出步骤7得到S″k(fn,fm),沿fn方向进行IFFT处理,得到S″k(n,fm),同时构造多通道系统回波矢量S(n,fm)=[S″1(n,fm) S″2(n,fm) … S″k(n,fm) … S″K(n,fm)]T;
其中[·]T表示矩阵转置。
其中[·]H表示矩阵共轭转置。
上述技术方案的特点和进一步改进在于:
步骤2具体操作如下:
步骤4具体操作如下:
其中,fc表示雷达发射信号的载频,c为电磁波的传播速度,V为卫星速度,θ为卫星的斜视角,dk为第k通道与参考通道沿卫星速度方向的距离,B为发射信号带宽,Δf为距离频域间隔,n=0,1,...,nrn-1,nrn表示距离向点数,tm为方位慢时间。
步骤6具体操作如下:
根据已知的雷达参数,第k个通道参考信号向量
其中,fc表示雷达发射信号的载频,R为卫星平台到场景中心斜距,c为电磁波的传播速度,V为卫星速度,θ为低轨道卫星的斜视角,dk为第k通道与参考通道沿卫星速度方向的距离,B为信号带宽,Δf为距离频域间隔,n=0,1,...,nrn-1,nrn表示距离向点数,tm为方位慢时间。
至此,一种基于子空间的星载沿方位多通道高分辨宽测绘带SAR-GMTI杂波抑制算法基本完成。
以下通过实测实验数据来进一步验证本发明的有效性。
(一)实测实验
1.实测参数
为了验证本发明方法的有效性,此处给出了表1中的实测数据参数。
表1实测数据参数
载频 | 5.0GHz | 通道数 | 6 |
平台速度 | 110m/s | 平台高度 | 5.5km |
场景中心斜距 | 15km | 多普勒模糊数 | 5 |
脉冲重复频率 | 200Hz | 发射信号带宽 | 150MHz |
2.实验内容
图1示意了利用本发明提出的基于子空间的方位多通道高分辨宽测绘带SAR-GMTI地面杂波抑制算法获得的实测数据处理结果。从图中可以看出本发明方法的杂波抑制效果,但是采用本发明方法的可以有效解决在通道误差存在情况的方位多通道高分辨宽测绘带SAR-GMTI地面杂波抑制问题。
综上所述,实测数据实验验证了本发明的正确性,有效性和可靠性。
Claims (1)
1.一种基于子空间的高分辨宽测绘带SAR-GMTI杂波抑制方法,其特征在于步骤如下:
步骤1:获取的多通道雷达系统的第k个通道回波信号是二维矩阵Sk(n,m),所述的Sk(n,m)为nrn×nan维矩阵,k=1,2,…,K,多通道系统共有K个通道;对矩阵Sk(n,m)列进行FFT处理,将结果保存在矩阵Sk(n,m)中;其中,nrn表示距离向点数,nan表示方位向点数;
步骤2:根据已知的雷达参数,构造参考信号向量S_ref(n)为nrn×1向量;其中,γ表示调频率,γ=B/Tp,B表示发射信号带宽,Tp表示发射脉冲宽度,fn表示为距离向频域坐标,B为发射信号带宽,Δf为距离频域间隔,
步骤3:取出步骤1得到Sk(n,m)的每一列,均点乘参考信号向量S_ref(n)的共轭,得到距离脉压后的第k通道数据矩阵Sk(fn,xm);xm表示方位向时域坐标,L表示为合成孔径长度,m=0,1,...,nan-1;
步骤4:根据已知的雷达参数,构造第k个通道参考信号向量S_refk(m)为1×nan向量;其中,fc表示雷达发射信号的载频,c为电磁波的传播速度,V为卫星速度,θ为卫星的斜视角,dk为第k通道与参考通道沿卫星速度方向的距离,tm为方位慢时间;
步骤5:取出步骤3得到的Sk(fn,xm)的每一行,均点乘参考信号向量S_refk(m)的共轭,得到距离走动校正之后的第k通道数据矩阵S′k(fn,xm);
步骤7:取出步骤5得到的S′k(fn,xm)点乘参考信号向量S_Chirpk(n,m)的共轭,并对相乘结果行进行FFT处理,得到chirp傅里叶变换之后的第k通道数据矩阵S″k(fn,fm);其中,fm表示为方位向频域坐标,PRF为方位采用频率,Δfa为方位频域间隔,
步骤8:取出步骤7得到S″k(fn,fm),沿fn方向进行IFFT处理,得到Sk″(n,fm),同时构造多通道系统回波矢量S(n,fm)=[S″1(n,fm) S″2(n,fm) … S″k(n,fm) … S″K(n,fm)]T;其中[·]T表示矩阵转置;
步骤9:取出步骤8得到的矢量对协方差矩阵进行估计,协方差矩阵
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911000657.7A CN110824473B (zh) | 2019-10-21 | 2019-10-21 | 一种基于子空间的高分辨宽测绘带sar-gmti杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911000657.7A CN110824473B (zh) | 2019-10-21 | 2019-10-21 | 一种基于子空间的高分辨宽测绘带sar-gmti杂波抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110824473A true CN110824473A (zh) | 2020-02-21 |
CN110824473B CN110824473B (zh) | 2022-11-08 |
Family
ID=69550053
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911000657.7A Active CN110824473B (zh) | 2019-10-21 | 2019-10-21 | 一种基于子空间的高分辨宽测绘带sar-gmti杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110824473B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113687356A (zh) * | 2021-09-16 | 2021-11-23 | 中国科学院空天信息创新研究院 | 一种机载多通道圆迹sar运动目标检测与估计方法 |
CN114114261A (zh) * | 2021-11-25 | 2022-03-01 | 哈尔滨工业大学 | 基于c2c的自监督sar海杂波抑制方法、系统、存储介质及设备 |
CN114780911A (zh) * | 2022-03-18 | 2022-07-22 | 西北工业大学 | 一种基于深度学习的海洋宽测绘带距离解模糊方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106093870A (zh) * | 2016-05-30 | 2016-11-09 | 西安电子科技大学 | 高超声速飞行器下降段的sar‑gmti杂波抑制方法 |
KR101687658B1 (ko) * | 2015-11-25 | 2016-12-19 | 한국항공우주연구원 | 처프-지 역변환 방법 및 시스템 |
CN106842200A (zh) * | 2017-01-11 | 2017-06-13 | 中国科学院电子学研究所 | 一种双基合成孔径雷达成像方法和装置 |
CN107092014A (zh) * | 2017-06-27 | 2017-08-25 | 西安电子科技大学 | 一种弹载双基前视sar舰船目标定位的优化方法 |
-
2019
- 2019-10-21 CN CN201911000657.7A patent/CN110824473B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101687658B1 (ko) * | 2015-11-25 | 2016-12-19 | 한국항공우주연구원 | 처프-지 역변환 방법 및 시스템 |
CN106093870A (zh) * | 2016-05-30 | 2016-11-09 | 西安电子科技大学 | 高超声速飞行器下降段的sar‑gmti杂波抑制方法 |
CN106842200A (zh) * | 2017-01-11 | 2017-06-13 | 中国科学院电子学研究所 | 一种双基合成孔径雷达成像方法和装置 |
CN107092014A (zh) * | 2017-06-27 | 2017-08-25 | 西安电子科技大学 | 一种弹载双基前视sar舰船目标定位的优化方法 |
Non-Patent Citations (2)
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 2012 (2012-02-26), pages 279 - 285 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113687356A (zh) * | 2021-09-16 | 2021-11-23 | 中国科学院空天信息创新研究院 | 一种机载多通道圆迹sar运动目标检测与估计方法 |
CN113687356B (zh) * | 2021-09-16 | 2023-10-20 | 中国科学院空天信息创新研究院 | 一种机载多通道圆迹sar运动目标检测与估计方法 |
CN114114261A (zh) * | 2021-11-25 | 2022-03-01 | 哈尔滨工业大学 | 基于c2c的自监督sar海杂波抑制方法、系统、存储介质及设备 |
CN114114261B (zh) * | 2021-11-25 | 2024-07-02 | 哈尔滨工业大学 | 基于c2c的自监督sar海杂波抑制方法、系统、存储介质及设备 |
CN114780911A (zh) * | 2022-03-18 | 2022-07-22 | 西北工业大学 | 一种基于深度学习的海洋宽测绘带距离解模糊方法 |
CN114780911B (zh) * | 2022-03-18 | 2024-03-01 | 西北工业大学 | 一种基于深度学习的海洋宽测绘带距离解模糊方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110824473B (zh) | 2022-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106093870B (zh) | 高超声速飞行器下降段的sar-gmti杂波抑制方法 | |
CN103744068B (zh) | 双通道调频连续波sar系统的动目标检测成像方法 | |
CN104914415B (zh) | 基于目标距离像模板匹配的单脉冲雷达相干干扰方法 | |
US8994584B2 (en) | Autofocus-based compensation (ABC) system and method for a hovering ground moving target indication (GMTI) sensor | |
CN110824473B (zh) | 一种基于子空间的高分辨宽测绘带sar-gmti杂波抑制方法 | |
CN112444811A (zh) | 一种综合mimo雷达和isar的目标检测和成像方法 | |
CN105445704B (zh) | 一种sar图像中的雷达动目标抑制方法 | |
CN108387884B (zh) | 基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法 | |
CN110488293B (zh) | 一种非均匀空间构型的分布式sar系统 | |
CN107728117B (zh) | 机载双基地sar一发两收杂波抑制方法 | |
CN108490443B (zh) | 基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法 | |
CN106291489B (zh) | 适用于多种发射信号波形的合成孔径雷达回波仿真方法 | |
CN111220955B (zh) | 基于垂直阵列协方差矩阵特征分解的机载气象雷达地杂波抑制方法 | |
CN113238225B (zh) | 一种用于动目标检测的雷达 | |
CN111796288B (zh) | 一种基于杂波频谱补偿技术的三坐标雷达动目标处理方法 | |
CN106291547B (zh) | 基于天线方向图辅助的多普勒模糊分量自适应抑制方法 | |
CN110133646B (zh) | 基于nlcs成像的双基前视sar的多通道两脉冲杂波对消方法 | |
CN113567978B (zh) | 一种多基分布式雷达协同成像方法 | |
CN113917424B (zh) | 一种基于ddma与盲源分离的地基mimo雷达地杂波抑制方法 | |
CN109143235B (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
CN113009473B (zh) | 一种多波束凝视雷达低仰角目标测高方法、装置及介质 | |
CN106526544B (zh) | 基于高超声速平台的mimosar杂波抑制方法 | |
CN115616629B (zh) | 一种基于天基外辐射源信号的运动目标探测补偿方法 | |
CN114780911A (zh) | 一种基于深度学习的海洋宽测绘带距离解模糊方法 | |
CN112946650B (zh) | 一站固定双站低频超宽带sar运动目标检测与成像方法 |
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 |