CN109143235B - 一种双基前视合成孔径雷达地面运动目标检测方法 - Google Patents
一种双基前视合成孔径雷达地面运动目标检测方法 Download PDFInfo
- Publication number
- CN109143235B CN109143235B CN201810971153.9A CN201810971153A CN109143235B CN 109143235 B CN109143235 B CN 109143235B CN 201810971153 A CN201810971153 A CN 201810971153A CN 109143235 B CN109143235 B CN 109143235B
- Authority
- CN
- China
- Prior art keywords
- time
- echo
- channel
- ground
- synthetic aperture
- 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.)
- Active
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种双基前视合成孔径雷达地面运动目标检测方法,针对BFSAR配置下地面动目标存在的问题,本发明提供的解决方案为:S1、对各通道的原始回波信号进行预滤波处理;S2、对步骤S1各通道的滤波结果进行一阶Keystone变换;S3、采用时分自适应抑制地面静物杂波;S4、采用改进的Wigner‑Ville分布方法将目标能量相干积累;本发明利用去斜预滤波器和Keystone变换来抑制多普勒模糊和校正跨距离单元徙动,通过时分自适应对消处理消除跨多普勒单元徙动的影响,有效抑制了地面静物杂波;并利用改进的Wigner‑Ville分布方法将目标能量相干积累,提高了信杂噪比。
Description
技术领域
本发明属于雷达技术领域,特别涉及一种双基前视SAR地面运动目标检测技术。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)是一种全天时、全天候的现代高分辨率微波遥感成像雷达,它利用雷达天线和目标区域间的相对运动来获得空间的高分辨率。在地形测绘、植被分析、海洋及水文观测、环境及灾害监视、资源勘探以及地壳微变检测等领域,合成孔径雷达发挥了越来越重要的作用。
双基前视SAR(BFSAR)是一种新的雷达体制,系统发射站和接收站分置于不同平台上。随着最近几年雷达合成孔径雷达的发展,BFSAR在动目标检测方面起着越来越重要的作用,特别是在军事领域。但是地面动目标回波往往会被周围具有多普勒模糊、跨距离单元徙动、跨多普勒单元徙动等特点的杂波所淹没,这增加了BFSAR对地面动目标检测的难度。
BFSAR地面目标检测(GMTD)目前主要基于两种方式:单通道方法和多通道方法。单通道方法主要基于多普勒滤波和多视干涉原理,见文献“Chen,H.C.,Mcgillem,C.D.:'Target motion compensation by sqectrum shifting in synthetic aperture radar',IEEE Transactions on Aerospace and Electronic Systems,2002,28,(3),pp.895-901”和文献“Ouchi,K.:'On the multilook images of moving targets by syntheticaperture radars',IEEE Transactions on Antennas and Propagation,2003,33,(8),pp.823-827”。虽然单通道SAR系统对硬件需求较低,运算量相对较小,但它要求慢速运动目标回波能量强于地物回波,这在实际应用中是很难满足的。并且,在机载平台下地物杂波谱会展宽导致单通道SAR系统对成像于主瓣内的运动目标更加难于检测。多通道方法主要包括相位中心偏置天线(DPCA)、沿航迹干涉(ATI)和空时自适应(STAP)方法;STAP方法见文献“Ender,J.H.G.:'Space-time processing for multichannel synthetic apertureradar',Electronics and Communication Engineering Journal,2002,11,(1),pp.29-38”和Barbarossa,S.,Farina,A.:'Space-time-frequency processing of syntheticaperture radar signals',IEEE Transactions on Aerospace and ElectronicSystems,1994,30,(2),pp.341-258;虽然STAP在一定程度上抑制了雷达回波中的杂波能量,提高了SCNR,但是在BFSAR回波中:距离徙动导致回波能量分散在多个距离单元,多普勒频谱展宽导致多普勒模糊,并且多普勒单元徙动导致多普勒频谱信号会占据多个距离单元;从而增加了在传统的STAP方法下检测BFSAR配置下地面动目标的难度。
发明内容
为了解决上述技术问题,本发明提出一种双基前视合成孔径雷达地面运动目标检测方法,采用三步式方法来检测BFSAR配置下的动目标信号,改进了检测信号的SCNR,提高了检测性能。
本发明采用的技术方案为:一种双基前视合成孔径雷达地面运动目标检测方法,包括:
S1、对各通道的原始回波信号进行预滤波处理;
S2、对步骤S1各通道的滤波结果进行一阶Keystone变换;
S3、采用时分自适应抑制地面静物杂波;
S4、采用改进的Wigner-Ville分布方法将目标能量相干积累。
进一步地,所述步骤S1之前还包括:
A1、初始化系统参数,包括:发射信号中心频率fc,带宽B,脉冲重复频率PRF,发射机平台速度VT,发射机平台位置(XT,YT,HT),接收机平台速度VR,接收机平台位置(0,0,HR),接收机通道数M,通道间隔d,合成孔径时间Ts,运动目标P的速度V,运动目标位置(XP,YP,0);
A2、收取各个通道的原始回波信号;将第m个通道的回波表示为Sm(η,τ),τ为快时间,η表示慢时间;
A3、对步骤A2收取的每个通道的原始回波信号进行距离向快速傅里叶变换,第m个通道距离向快速傅里叶变换的结果表示为Sm(η,f)=FFTrg{Sm(η,τ)},FFTrg表示距离向快速傅里叶变换运算,f表示距离向频率。
进一步地,步骤S1所述预滤波采用去斜预滤波函数。
更进一步地,去斜预滤波函数表达式为:
其中,fdc表示多普勒质心。
进一步地,步骤S3具体包括以下分步骤:
S31、对步骤S2得到的回波信号进行距离向反傅里叶变换,同时对步骤S2得到的回波信号每一行进行距离压缩;
S32、对步骤S31得到的回波信号每个距离单元进行列矢量化处理;
S33、对经步骤S32处理后的回波信号方位向进行时间分段,并计算各个时间片段回波信号的最佳权向量;
S34、将每个时间片段回波信号与最佳权向量相乘,得到抑制地面静物杂波后的回波信号。
进一步地,步骤S33所述最佳权向量为:
其中,δ表示时间序列,R-1(δ)表示杂波协方差的逆,S(δ)表示空时二维导向矢量。
更进一步地,杂波协方差R(δ)的计算式为:
其中,χfrag(i)表示第i个距离单元的采用时间分段方法得到的子向量。N表示选取的杂波距离单元个数。
进一步地,步骤S4具体为:
S41、根据步骤S3抑制地面静物杂波后的回波信号,得到回波数据在动目标点的方位向回波;
S42、对步骤S41的方位向回波进行WVD变换;
S43、对步骤S42得到的变换结果进行逆傅里叶变换;
S44、对步骤S43得到的变换结果进行二维傅里叶变换。
本发明的有益效果:本发明的方法,首先利用去斜预滤波器和Keystone变换来抑制多普勒模糊和校正跨距离单元徙动,然后通过时分自适应对消处理消除跨多普勒单元徙动的影响,有效抑制了地面静物杂波;最后,利用改进的Wigner-Ville分布方法(MWVD)将目标能量相干积累,进一步提高信杂噪比(SCNR),从而实现BFSAR地面目标运动目标检测。
附图说明
图1为本发明的方案流程图;
图2为本发明实施例提供的BFSAR空间几何构型;
图3为本发明实施例提供的通过时分自适应对消处理的目标处方位向回波;
图4为本发明实施例提供的经过MWVD后的目标处方位向回波。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
如图1所示为本发明的方案流程图,本发明的一种双基前视合成孔径雷达地面运动目标检测方法,包括:
S1、对各通道的原始回波信号进行预滤波处理;
S2、对步骤S1各通道的滤波结果进行一阶Keystone变换;
S3、采用时分自适应抑制地面静物杂波;
S4、采用改进的Wigner-Ville分布方法将目标能量相干积累。
步骤S1之前还包括:
A1、初始化系统参数,如图2所示为本实施例采用的BFSAR空间几何构型,初始化参数具体包括:发射信号中心频率fc,带宽B,脉冲重复频率PRF,发射机平台速度VT,发射机平台位置(XT,YT,HT),接收机平台速度VR,接收机平台位置(0,0,HR),接收机通道数M,通道间隔d,合成孔径时间Ts,运动目标P的速度V,运动目标位置(XP,YP,0);
A2、收取各个通道的原始回波信号;将第m个通道的回波表示为Sm(η,τ),τ为快时间,η表示慢时间;各个通道的原始回波信号表示为:
其中,τ为快时间,η表示慢时间,ωr和ωa为距离向和方位向包络。Ts为合成孔径时间,Kr为距离向调频率,fc为发射信号中心频率,c为光速,ηR为相对于零时刻时目标位于发射波束中心的时间差,Rm(η)表示不同时刻目标点到发射机和接收机的距离和。
Rm(η)=RT(η)+RR-m(η)
其中,RT(η)表示在η时刻发射机到目标点的距离,RR-m(η)表示在η时刻第m个接收机到目标点的距离。
A3、对步骤A2收取的每个通道的原始回波信号进行距离向快速傅里叶变换,第m个通道距离向快速傅里叶变换的结果表示为Sm(η,f)=FFTrg{Sm(η,τ)},FFTrg表示距离向快速傅里叶变换运算,f表示距离向频率。
对每个通道的回波信号进行距离向快速傅里叶变换得到距离频域-方位时间域,并对Rm(η)进行泰勒展开:
其中,f表示距离频域,Rb0为η0时刻的双基距离和,R'b0和R”b0分别表示Rm(η)在η=η0处泰勒展开式的一阶导数和二阶导数在η=η0的值R'm(η0)和R'm(η0)。
步骤S1具体为:
去除回波多普勒模糊得到S'm(η,f),去斜预滤波器函数为:
其中,fdc表示多普勒质心。
然后让Sm(η,f)每一行通过去斜预滤波器:
步骤S2具体为:
对步骤S1中的结果进行一阶Keystone变换。即对上一步滤波后的回波进行变量变换η1=(f+fc)η/fc得到S'm(η1,f):
其中,η1是变换后新的方位时间,λ为载波波长。这样回波中的线性距离徙动分量就被矫正;高阶距离徙动分量依然保留,但是在BFSAR中可以忽略。
步骤S3具体包括以下分步骤:
S31、对步骤S2得到的回波信号进行距离向反傅里叶变换,同时对步骤S2得到的回波信号每一行进行距离压缩;
对S'm(η1,f)进行距离向反傅里叶变换得到S'm(η1,τ1),同时对回波信号每一行进行距离压缩,压缩函数为:
h(t)=s*(-t)
其中,s(t)为雷达发射机发射信号,
从而得到压缩后的回波为新的S'm(η1,τ1)。τ1为变换后新的距离向时间。
S32、对步骤S31得到的回拨信号每个距离单元进行列矢量化处理;
设经过步骤六得到的S'm(η1,τ1)的阶数维L×K维,则回波距离向点数为L,方位向点数为K。同时对回波进行距离压缩。则建立第l(0<l<L)个距离单元的数据矩阵为
其中,Sijl表示第i个方位时间,第j个接收阵元,第l个距离单元的回波采样值。对每个距离单元的回波数据为了后续步骤进行如下处理得到χ(l):
χ(l)=vec(Xl)=[x1,l;x2,l;…;xK,l]
其中,vec(·)表示对矩阵进行列矢量化处理。对应得l距离的空时导向矢量为:χ(l)
S33、对经步骤S32处理后的回波信号方位向进行时间分段,并计算各个时间片段回波信号的最佳权向量;
根据表1数据对方位向进行时间分段并计算各个时间片段的最佳权向量wopt,
表1BFSAR的参数表
参数 | 数值 |
中心频率 | 10GHz |
带宽 | 300MHz |
PRF | 1500Hz |
合成孔径时间 | 0.5s |
平台速度 | (0,200,0)m/s |
接收机通道数 | 3 |
通道间隔 | 1m |
发射机位置 | (8000,-2000,8000)m |
运动目标位置 | (0,0,0)m |
运动目标速度 | (3,-3,0)m/s |
即:
其中,δ表示时间序列,R-1(δ)表示杂波协方差的逆,S(δ)表示空时二维导向矢量。
其中,χfrag(i)表示第i个距离单元的采用时间分段方法得到的子向量。N表示选取的杂波距离单元个数。
S34、将每个时间片段回波信号与最佳权向量相乘wopt(δ)'*χfrag(i),得到目标所在单元抑制地面静物杂波后的回波信号如图3所示。从图3可以看出,虚线表示的原始回波信号中目标信号淹没在杂波信号中,无法实现动目标检测;实线为杂波抑制后的回波信号,可以看出信号的幅值明显高于周围杂波,SCNR大大提高。
步骤S4具体为:
滤波后的回波数据在动目标点的方位向回波为:
其中,α、β无具体的物理含义,本申请为便于计算所采用的中间变量,α=-λfdc,β=-λddr,fdc和fdr分别表示动目标的多普勒质心和多普勒调频率,G表示信号的幅度。
接着对信号Sfiltered(η1)进行WVD变换,得到:
其中,t表示滞后时间,(.)*表示共轭变换。
然后对信号进行逆傅里叶变换(IFFT)并进行变量替换η'1=η1t,得到:
对MWVDs(η'1,t)进行二维傅里叶变换得到
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (7)
1.一种双基前视合成孔径雷达地面运动目标检测方法,其特征在于,包括:
S1、对各通道的原始回波信号进行预滤波处理;
S2、对步骤S1各通道的滤波结果进行一阶Keystone变换;
S3、采用时分自适应抑制地面静物杂波;步骤S3具体包括以下分步骤:
S31、对步骤S2得到的回波信号进行距离向反傅里叶变换,同时对步骤S2得到的回波信号每一行进行距离压缩;
S32、对步骤S31得到的回波信号每个距离单元进行列矢量化处理;
S33、对经步骤S32处理后的回波信号方位向进行时间分段,并计算各个时间片段回波信号的最佳权向量;
S34、将每个时间片段回波信号与最佳权向量相乘,得到抑制地面静物杂波后的回波信号;
S4、采用改进的Wigner-Ville分布方法将目标能量相干积累;步骤S4具体为:
S41、根据步骤S3抑制地面静物杂波后的回波信号,得到回波数据在动目标点的方位向回波;
S42、对步骤S41的方位向回波进行Wigner-Ville分布变换;
S43、对步骤S42得到的变换结果进行逆傅里叶变换;
S44、对步骤S43得到的变换结果进行二维傅里叶变换。
2.根据权利要求1所述的一种双基前视合成孔径雷达地面运动目标检测方法,其特征在于,所述步骤S1之前还包括:
A1、初始化系统参数,包括:发射信号中心频率fc,带宽B,脉冲重复频率PRF,发射机平台速度VT,发射机平台位置(XT,YT,HT),接收机平台速度VR,接收机平台位置(0,0,HR),接收机通道数M,通道间隔d,合成孔径时间Ts,运动目标P的速度V,运动目标位置(XP,YP,0);
A2、收取各个通道的原始回波信号;将第m个通道的回波表示为Sm(η,τ),τ为快时间,η表示慢时间;
A3、对步骤A2收取的每个通道的原始回波信号进行距离向快速傅里叶变换,第m个通道距离向快速傅里叶变换的结果表示为Sm(η,f)=FFTrg{Sm(η,τ)},FFTrg表示距离向快速傅里叶变换运算,f表示距离向频率。
3.根据权利要求2所述的一种双基前视合成孔径雷达地面运动目标检测方法,其特征在于,步骤S1所述预滤波采用去斜预滤波函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810971153.9A CN109143235B (zh) | 2018-08-24 | 2018-08-24 | 一种双基前视合成孔径雷达地面运动目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810971153.9A CN109143235B (zh) | 2018-08-24 | 2018-08-24 | 一种双基前视合成孔径雷达地面运动目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109143235A CN109143235A (zh) | 2019-01-04 |
CN109143235B true CN109143235B (zh) | 2020-09-08 |
Family
ID=64827789
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810971153.9A Active CN109143235B (zh) | 2018-08-24 | 2018-08-24 | 一种双基前视合成孔径雷达地面运动目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109143235B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110109113B (zh) * | 2019-05-07 | 2021-01-12 | 电子科技大学 | 一种基于级联对消的双基前视sar非平稳杂波抑制方法 |
CN110471060B (zh) * | 2019-07-08 | 2021-03-26 | 中国空间技术研究院 | 一种高效的雷达动目标调频率估计方法 |
CN113435299B (zh) * | 2021-06-23 | 2022-05-13 | 电子科技大学 | 基于空-时匹配的双基前视sar杂波抑制方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7969345B2 (en) * | 2009-04-13 | 2011-06-28 | Raytheon Company | Fast implementation of a maximum likelihood algorithm for the estimation of target motion parameters |
CN102147469B (zh) * | 2010-12-29 | 2012-11-07 | 电子科技大学 | 一种双基地前视合成孔径雷达成像方法 |
CN102353953A (zh) * | 2011-09-15 | 2012-02-15 | 西安电子科技大学 | 单通道sar地面运动目标成像方法 |
CN104730505A (zh) * | 2014-12-08 | 2015-06-24 | 广西大学 | 一种多通道sar地面目标检测定位方法和装置 |
CN106772309B (zh) * | 2017-03-31 | 2019-08-02 | 西安电子科技大学 | 基于通信卫星辐射源的雷达目标检测方法 |
CN108132467B (zh) * | 2017-12-23 | 2018-11-30 | 成都汇蓉国科微系统技术有限公司 | 基于增强型adc的dsp+fpga双基前视sar成像方法及成像装置 |
-
2018
- 2018-08-24 CN CN201810971153.9A patent/CN109143235B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109143235A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107229048B (zh) | 一种高分宽幅sar动目标速度估计与成像方法 | |
CN106093870B (zh) | 高超声速飞行器下降段的sar-gmti杂波抑制方法 | |
CN108051809B (zh) | 基于Radon变换的运动目标成像方法、装置及电子设备 | |
CN105445704B (zh) | 一种sar图像中的雷达动目标抑制方法 | |
CN110412568A (zh) | 基于扩展方位角相位编码的距离模糊杂波抑制方法 | |
CN110109113B (zh) | 一种基于级联对消的双基前视sar非平稳杂波抑制方法 | |
CN107942327A (zh) | 基于脉冲相位编码的单通道hrws‑sar成像方法 | |
CN109143235B (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
CN106772253B (zh) | 一种非均匀杂波环境下的雷达杂波抑制方法 | |
CN109507666B (zh) | 基于离网变分贝叶斯算法的isar稀疏频带成像方法 | |
CN114545411B (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
CN110568434A (zh) | 一种多通道匀加速轨迹sar动目标二维速度估计方法 | |
CN105301589B (zh) | 高分辨宽测绘带sar地面运动目标成像方法 | |
CN104808204A (zh) | 一种固定发射站双基前视sar动目标检测方法与成像方法 | |
CN109031299B (zh) | 低信噪比条件下基于相位差分的isar平动补偿方法 | |
CN110133646B (zh) | 基于nlcs成像的双基前视sar的多通道两脉冲杂波对消方法 | |
Yu et al. | Ground moving target motion parameter estimation using Radon modified Lv's distribution | |
CN115877382A (zh) | 基于调频连续波相邻脉冲变换差分的运动误差估计方法 | |
Liu et al. | Ambiguities Suppression for Azimuth Multichannel SAR Based on ${L_ {2, q}} $ Regularization With Application to Gaofen-3 Ultra-Fine Stripmap Mode | |
Wang et al. | Clutter suppression and ground moving target imaging approach for hypersonic vehicle borne multichannel radar based on two-step focusing method | |
CN104931965B (zh) | 一种固定发射站双基正侧视sar成像方法 | |
Zhou et al. | Signal model and moving target detection based on MIMO synthetic aperture radar | |
CN103217677A (zh) | 一种基于联合检测量的单通道sar动目标检测方法 | |
CN113484862A (zh) | 一种自适应的高分宽幅sar清晰重构成像方法 | |
CN114780911B (zh) | 一种基于深度学习的海洋宽测绘带距离解模糊方法 |
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 |