CN112505695A - 基于星载sar系统的空中机动目标检测方法 - Google Patents
基于星载sar系统的空中机动目标检测方法 Download PDFInfo
- Publication number
- CN112505695A CN112505695A CN202011215795.XA CN202011215795A CN112505695A CN 112505695 A CN112505695 A CN 112505695A CN 202011215795 A CN202011215795 A CN 202011215795A CN 112505695 A CN112505695 A CN 112505695A
- Authority
- CN
- China
- Prior art keywords
- target
- sar system
- maneuvering target
- detection method
- satellite
- 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.)
- Pending
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 24
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 230000001427 coherent effect Effects 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000001228 spectrum Methods 0.000 claims abstract description 11
- 238000009825 accumulation Methods 0.000 claims abstract description 10
- 230000001629 suppression Effects 0.000 claims abstract description 10
- 230000001133 acceleration Effects 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 abstract description 8
- 238000005314 correlation function Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000035485 pulse pressure Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
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/9094—Theoretical aspects
-
- 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系统的空中机动目标检测方法,包括对进行地面杂波抑制后的方位向相邻的雷达回波信号进行互相关操作;对进行互相关操作处理后的动目标回波信号进行统一的keystone变换;对keystone变换后的回波信号进行变尺度的WVD变换,使得动目标冲激线性谱平行于时间轴,然后对该冲激线性谱沿着时间轴做FFT变换实现机动目标信号的相干积累等;本发明具有极大的灵活性,更具有实际应用的价值,适用于低信噪比下对空中机动目标的检测等优点。
Description
技术领域
本发明涉及雷达回波数据处理领域,更为具体的,涉及一种基于星载 SAR系统的空中机动目标检测方法。
背景技术
要捕获空中机动目标必须对目标进行检测,发现目标后才能进行有效 的跟踪和成像。因此,对空中运动目标进行检测尤为重要。
目前对空中机动目标的检测主要通过逆合成孔径雷达(Inverse SyntheticAperture Radar,ISAR)系统实现,通过地基或者岸基雷达发射电 磁波照射地面或海面上方的目标获取其后向散射信息,进而实现空中运动 目标的检测。然而,ISAR系统只能部署在陆面等固定场景,其灵活性受到 极大的限制。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于星载SAR系统 的空中机动目标检测方法,具有极大的灵活性,更具有实际应用的价值, 适用于低信噪比下对空中机动目标的检测等优点。
本发明的目的是通过以下方案实现的:
一种基于星载SAR系统的空中机动目标检测方法,包括步骤:
对进行地面杂波抑制后的方位向相邻的雷达回波信号进行互相关操作; 对进行互相关操作处理后的动目标回波信号进行统一的keystone变换;对 keystone变换后的回波信号进行变尺度的WVD变换,使得动目标冲激线性 谱平行于时间轴,然后对该冲激线性谱沿着时间轴做FFT变换实现机动目 标信号的相干积累。
进一步地,地面杂波抑制包括步骤:对星载SAR系统的雷达回波数据 进行DPCA相消处理。
进一步地,所述雷达回波数据包括双通道雷达回波数据。
进一步地,在地面杂波抑制前包括对空中机动目标进行建模的步骤, 该步骤中,采用如下公式计算运动目标距离第m个通道的距离Rm(ta),
其中,卫星沿X轴方向速度为v,飞机初始位置为(0,Rb),径向速度 vr(ta)=vr0+arta,横向速度va(ta)=va0+aata,vr0为目标初始径向速度,ar为目标 径向加速度,va0为目标初始横向速度,aa为目标横向加速度,ta为方位慢 时间,m=1,2,d1=-D,d2=0。
进一步地,实现机动目标信号的相干积累后,在检测平面对SAR机动 目标进行检测。
进一步地,将检测方法的结果用于参数估计和/或成像聚焦操作。
本发明的有益效果是:
本发明提出了一种基于星载SAR平台实现空中机动目标的检测 方法,相比于ISAR系统实现空中运动目标检测具有极大的灵活性。 同时,本发明充分考虑了空中运动目标(如飞机、导弹等)的机动性, 区别于现有SAR-GMTI方法对空中机动目标进行匀速直线运动的假 设,因此本发明更具有实际应用的价值。具体的,本发明采用相邻互 相关方法对运动目标进行降阶处理,通过keystone方法对动目标进行 统一的走动校正,然后对动目标信号进行变尺度的WVD变换,最后 通过傅里叶变换能够实现对机动目标的有效相干积累,而噪声信号是 非相干的,进行WVD变换后其分布仍然是杂乱的,因而本发明适用 于低信噪比下对空中机动目标的检测。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对 实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地, 下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员 来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的 附图。
图1为星载SAR平台运动目标观测在成像几何平面的投影模型示意图;
图2为本发明空中机动目标检测的流程示意图;
图3为利用本发明对空中机动目标检测结果;
图4为图3的检测结果对应的距离剖面图;
图5为图3的检测结果对应的方位剖面图。
具体实施方式
本说明书中所有实施例公开的所有特征(包括任何附加权利要求、摘 要和附图),或隐含公开的所有方法或过程中的步骤,除了互相排斥的特征 和/或步骤以外,均可以以任何方式组合和/或扩展、替换。
如图1~5所示,一种基于星载SAR系统的空中机动目标检测方法,包 括步骤:
对进行地面杂波抑制后的方位向相邻的雷达回波信号进行互相关操作; 对进行互相关操作处理后的动目标回波信号进行统一的keystone变换;对 keystone变换后的回波信号进行变尺度的WVD变换,使得动目标冲激线性 谱平行于时间轴,然后对该冲激线性谱沿着时间轴做FFT变换实现机动目 标信号的相干积累。
进一步地,地面杂波抑制包括步骤:对星载SAR系统的雷达回波数据 进行DPCA相消处理。
进一步地,所述雷达回波数据包括双通道雷达回波数据。
进一步地,在地面杂波抑制前包括对空中机动目标进行建模的步骤, 该步骤中,采用如下公式计算运动目标距离第m个通道的距离Rm(ta),
其中,卫星沿X轴方向速度为v,飞机初始位置为(0,Rb),径向速度 vr(ta)=vr0+arta,横向速度va(ta)=va0+aata,vr0为目标初始径向速度,ar为目标 径向加速度,va0为目标初始横向速度,aa为目标横向加速度,ta为方位慢 时间,m=1,2,d1=-D,d2=0。
进一步地,实现机动目标信号的相干积累后,在检测平面对SAR机动 目标进行检测。
进一步地,将检测方法的结果用于参数估计和/或成像聚焦操作。
在本发明的其他实施例中,结合附图对本发明基于星载SAR系统的空 中机动目标检测方法进行详细说明。
图1为星载SAR平台对飞行中的飞机进行观测的工作几何模型在成像 斜距平面的投影,其中卫星沿X轴方向速度为v,飞机初始位置为(0,Rb), 径向速度vr(ta)=vr0+arta,横向速度va(ta)=va0+aata,vr0为目标初始径向速度, ar为目标径向加速度,va0为目标初始横向速度,aa为目标横向加速度,这 里充分考虑了空中运动目标机动特性。运动目标距离第m个通道的距离为
其中,ta为方位慢时间,m=1,2,d1=-D,d2=0。
图2为本发明方法实施例的流程示意图。
参照图2,首先在步骤S102中,为对地面杂波进行有效抑制,对双通 道SAR系统进行DPCA处理。设雷达发射线性调频信号,则对多通道接收 的该动目标的基频回波信号并行进距离脉压后为
式中tr为距离快时间,σn为目标散射系数,c为光速,Ar(·)和Aa(·)分别为 雷达LFM信号的窗函数和方位窗函数,γ是发射的线性调频信号的调频率, λ=c/fc为中心频率对应的波长。
取前后两个通道进行DPCA相消,可得相消后的结果
在步骤S103中,对方位相邻的回波进行互相关操作,式(3)回波信号的 相邻互相关表达式为
R(τr,tm)=∫I(tr,tm)I*(tr-τr,tm+1)dt (4)
其中tm为方位离散后的时间,可以表示为tm=m/PRF(PRF为脉冲重复频 率,m=1,2,3,…),公式(3)中的常数相位项并不影响后面的分析,为公式表 达的简洁,这里先将其忽略。公式(3)的距离频域表达式为
根据时域相关等于频域共轭点乘,然后再作IFFT变换的原理可以知道, 式(4)可以进一步表示为
其中
上式中忽略了指数项由于其相位值大小远小于式中A为相关匹配量,从式(6)中可以看出,第一个指数项为常数项,在不 同的方位时刻,其对应的值造成的相关函数的包络移动量是一个定值;第 二个指数项为方位时间的一次项,所对应的距离耦合量会产生相关函数峰 值的走动;第三个指数项为方位时间的二次项,所对应的距离耦合量会产 生相关函数峰值的弯曲。
在步骤S104中,将得到的互相关函数变换到距离频域表示式为
对式(7)进行统一keystone变换,令
tm(fc+fr)=fcτm (8)
经过统一keystone变换后,式(7)变为
其中τm为新的方位时刻,从式(9)可以看出,产生相邻相关函数峰值走 动的耦合量通过统一keystone尺度变换已经消除了。一般情况下,可以忽 略相关函数弯曲量对距离包络的影响,但相位精度要求远高于包络精度, 所以,对方位相位的影响不能忽略。由于fc>>fr,则所以,对式 (9)做距离向的IFFT,并做相应近似得到
从式(10)可以看出,峰值位置沿方位分布时,仍然在相同的距离单元内, 这就为后面沿方位积累在低信噪比情况下检测运动目标提供了前提。
在步骤S105中,对得到的信号进行时频变换分析,这里以魏格纳威利 (WVD)为例。由于式(10)中的窗函数并不影响后续分析,为了方便方 位向的分析,暂时忽略式(10)中的窗函数,故可以简写为
对式(11)进行变尺度的WVD变换,令τm=μτm′/t,μ为尺度变换常量,τm′ 为新的方位时间变量,t为瞬时偏移量,可以得到
从式(19)中可以看出,在变尺度的WVD时频平面内,其运动目标 信号分布呈现为一条平行于新时间轴的线性谱,线性谱只存在一次项相位 没有二次及以上的高次相位。此时,对该线性谱直接进行傅里叶变换就可 以实现信号的相干积累。对式(22)进行关于τm′的FFT变换得到
在步骤S106中,在检测平面对空中机动目标进行检测,在式(13)生成 的整个检测平面内,对目标进行检测,目标信号的能量通过相干完成积累, 而噪声信号是非相干的,在检测平面内仍然是杂乱的,目标积累能量最大 值的位置位于面内为(2N1/λ,4μN2/λ)处,因此,通过本发明提出的方法可以很 好的在低信噪比情况下对空中机动目标进行检测。基于本发明提出的空中 机动目标检测方法,后续可以完成运动目标的参数估计与成像聚焦操作。
图3为空中机动目标的在检测平面的检测结果,为验证本发明提出的 方法在低信噪比下机动目标检测的性能,在进行DPCA后,对回波数据加 入SNR为-20dB的高斯噪声,从得到的动目标检测结果我们可以看到动目 标在检测的二维平面中形成了一个尖峰,这说明利用本发明所提方法动目 标得到了有效的相干积累,而噪声信号是非相干的,进行WVD变换后其分 布仍然是杂乱的,因而本发明适用于低信噪比下对空中机动目标的检测, 图中X轴为距离采样单元,Y轴为方位采样单元,Z轴为归一化幅度(单位: dB)。
图4为动目标检测结果在距离向剖面图,图中横轴为距离采样单元, 纵轴为归一化幅度(单位:dB)。
图5为动目标检测结果在方位向剖面图,图中横轴为方位采样单元, 纵轴为归一化幅度(单位:dB)。从图4和图5中可以看到噪声在-20dB的 水平,这组实验充分验证了本发明所提方法对空中机动目标检测的有效性。
本发明提出了一种基于星载SAR平台实现空中机动目标的检测方法, 相比于ISAR系统实现空中运动目标检测具有极大的灵活性。同时,本发明 充分考虑了空中运动目标(如飞机、导弹等)的机动性,区别于现有 SAR-GMTI方法对空中机动目标进行匀速直线运动的假设,因此本发明更 具有实际应用的价值。具体的,本发明采用相邻互相关方法对运动目标进 行降阶处理,通过keystone方法对动目标进行统一的走动校正,然后对动 目标信号进行变尺度的WVD变换,最后通过傅里叶变换能够实现对机动目 标的有效相干积累,而噪声信号是非相干的,进行WVD变换后其分布仍然 是杂乱的,因而本发明适用于低信噪比下对空中机动目标的检测。
本发明功能如果以软件功能单元的形式实现并作为独立的产品销售或 使用时,可以存储在一个计算机可读存储介质中。基于这样的理解,本发 明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的 部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储 介质中,在一台计算机设备(可以是个人计算机,服务器,或者网络设备 等)以及相应的软件中执行本发明各个实施例所述方法的全部或部分步骤。 而前述的存储介质包括:U盘、移动硬盘、或者光盘等各种可以存储程序 代码的介质,进行测试或者实际的数据在程序实现中存在于只读存储器 (Random Access Memory,RAM)、随机存取存储器(Random Access Memory,RAM)等。
Claims (6)
1.一种基于星载SAR系统的空中机动目标检测方法,其特征在于,包括步骤:
对进行地面杂波抑制后的方位向相邻的雷达回波信号进行互相关操作;对进行互相关操作处理后的动目标回波信号进行统一的keystone变换;对keystone变换后的回波信号进行变尺度的WVD变换,使得动目标冲激线性谱平行于时间轴,然后对该冲激线性谱沿着时间轴做FFT变换实现机动目标信号的相干积累。
2.根据权利要求1所述的基于星载SAR系统的空中机动目标检测方法,其特征在于,地面杂波抑制包括步骤:对星载SAR系统的雷达回波数据进行DPCA相消处理。
3.根据权利要求1或2所述的基于星载SAR系统的空中机动目标检测方法,其特征在于,所述雷达回波数据包括双通道雷达回波数据。
5.根据权利要求1所述的基于星载SAR系统的空中机动目标检测方法,其特征在于,实现机动目标信号的相干积累后,在检测平面对SAR机动目标进行检测。
6.根据权利要求1~5任一所述的基于星载SAR系统的空中机动目标检测方法,其特征在于,将检测方法的结果用于参数估计和/或成像聚焦操作。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011215795.XA CN112505695A (zh) | 2020-11-04 | 2020-11-04 | 基于星载sar系统的空中机动目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011215795.XA CN112505695A (zh) | 2020-11-04 | 2020-11-04 | 基于星载sar系统的空中机动目标检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112505695A true CN112505695A (zh) | 2021-03-16 |
Family
ID=74955333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011215795.XA Pending CN112505695A (zh) | 2020-11-04 | 2020-11-04 | 基于星载sar系统的空中机动目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112505695A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375128A (zh) * | 2014-11-07 | 2015-02-25 | 电子科技大学 | 一种基于互相关函数的高机动目标快速积累检测方法 |
CN105445704A (zh) * | 2015-11-27 | 2016-03-30 | 西安电子科技大学 | 一种sar图像中的雷达动目标抑制方法 |
CN109669183A (zh) * | 2017-12-27 | 2019-04-23 | 北京航空航天大学 | 一种基于Keystone和时频变换的地球同步轨道SAR运动目标成像处理装置 |
-
2020
- 2020-11-04 CN CN202011215795.XA patent/CN112505695A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375128A (zh) * | 2014-11-07 | 2015-02-25 | 电子科技大学 | 一种基于互相关函数的高机动目标快速积累检测方法 |
CN105445704A (zh) * | 2015-11-27 | 2016-03-30 | 西安电子科技大学 | 一种sar图像中的雷达动目标抑制方法 |
CN109669183A (zh) * | 2017-12-27 | 2019-04-23 | 北京航空航天大学 | 一种基于Keystone和时频变换的地球同步轨道SAR运动目标成像处理装置 |
Non-Patent Citations (4)
Title |
---|
兰毅: "双基地雷达动目标能量积累方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
李亚超: "ISAR目标运动参数估计及成像技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
李学仕: "高分辨率宽测绘带SAR动目标处理方法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
李道京, 北京:国防工业出版社 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7969345B2 (en) | Fast implementation of a maximum likelihood algorithm for the estimation of target motion parameters | |
Huang et al. | A fast SAR imaging method for ground moving target using a second-order WVD transform | |
CN111142105B (zh) | 复杂运动目标isar成像方法 | |
CN108051809B (zh) | 基于Radon变换的运动目标成像方法、装置及电子设备 | |
Su et al. | High-speed multi-target detection with narrowband radar | |
US7928896B2 (en) | Application of time reversal to synthetic aperture imaging | |
CN105301590B (zh) | 一种机动目标调频步进逆合成孔径成像方法 | |
Kang et al. | Ground moving target imaging based on compressive sensing framework with single-channel SAR | |
Sjogren et al. | Moving target relative speed estimation and refocusing in synthetic aperture radar images | |
CN107085213A (zh) | 基于随机调频步进波形设计的运动目标isar成像方法 | |
CN109407070B (zh) | 一种高轨平台地面运动目标检测方法 | |
CN110646774B (zh) | 基于乘积变尺度周期吕分布的机动目标相参检测方法及装置 | |
CN109669182B (zh) | 无源双基地sar动/静目标联合稀疏成像方法 | |
Wacks et al. | Passive synthetic aperture hitchhiker imaging of ground moving targets—Part 1: Image formation and velocity estimation | |
CN109613507B (zh) | 一种针对高阶机动目标雷达回波的检测方法 | |
CN105445704A (zh) | 一种sar图像中的雷达动目标抑制方法 | |
CN113009473B (zh) | 一种多波束凝视雷达低仰角目标测高方法、装置及介质 | |
CN108562884A (zh) | 一种基于最大后验概率的机载前视海面目标角超分辨方法 | |
CN109001700B (zh) | 一种实现目标轮廓重构的雷达前视成像方法 | |
CN108387884A (zh) | 基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法 | |
Yu et al. | Ground moving target motion parameter estimation using Radon modified Lv's distribution | |
Fu et al. | A novel method for fast detection of high-speed targets | |
Chen et al. | Compensation for high-frequency vibration of SAR imaging in the terahertz band based on linear chirplet transform and empirical mode decomposition | |
CN106772373B (zh) | 针对任意地面运动目标的sar成像方法 | |
CN112327301B (zh) | 基于子孔径grft的低信噪比下参数化平动补偿快速方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210316 |