CN108051812A - 基于二维速度搜索的星载sar运动目标检测方法 - Google Patents
基于二维速度搜索的星载sar运动目标检测方法 Download PDFInfo
- Publication number
- CN108051812A CN108051812A CN201810058355.4A CN201810058355A CN108051812A CN 108051812 A CN108051812 A CN 108051812A CN 201810058355 A CN201810058355 A CN 201810058355A CN 108051812 A CN108051812 A CN 108051812A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- moving target
- orientation
- search
- 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
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
- G01S13/904—SAR modes
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于二维速度搜索的星载SAR运动目标检测方法,属于雷达目标检测技术领域。本发明针对运动目标具体状态的未知性,结合RD算法,通过遍历一个间隔均匀的二维数组记为vx(m)、vy(n)(分别代表方位向搜索速度和距离向搜索速度),分别设计距离徙动校正函数和方位向匹配函数,实现运动目标的聚焦成像。然后对方位向滤波器进行修正,消除运动目标距离向速度带来的成像方位向偏移。提取不同搜索速度下的运动目标的最强幅度值用于恒虚警检测,将显著提高运动目标的检测概率。仿真结果验证了该方法的可行性和有效性。
Description
技术领域
本发明公开了一种基于二维速度搜索的星载SAR运动目标检测方法,属于雷达目标检测技术领域。
背景技术
合成孔径雷达(SAR)系统利用发射的距离向宽带信号和平台运动带来的方位向线性调频信号,实现地面静止目标成像。作为一种成像雷达,星载SAR对静止目标的成像技术已经趋于完善和成熟。随着星载SAR技术的发展,利用星载SAR进行地面运动目标的检测技术逐渐引起人们的重视。运动目标的检测和成像是现代雷达应该完成的基本功能之一,也成为合成孔径雷达重要的发展方向。
由于星载SAR处于下视状态,回波中不可避免的包含了大量的噪声与杂波。如果直接对收到的场景回波进行SAR成像处理,运动目标通常会处于散焦状态,则信噪比与信杂比较低,这大大降低了雷达的检测性能。因此,必须尽可能的改善星载SAR系统的信噪比与信杂比,才能提高运动目标的检测性能。
由于星载SAR是利用雷达与场景之间的相对运动,通过信号处理和运动补偿的方法,得到高的方位向分辨率。对静止目标来说,目标和雷达之间的相对速度只是由雷达平台引起,而运动目标由于自身有一定的运动特性,使得运动目标在整个合成孔径内的变化规律与静止目标不同,导致二者的多普勒频率特性存在差异。如果使用静止目标的参考函数对运动目标成像,必然会造成运动目标成像质量的下降,出现散焦现象。因此,需要针对运动目标的特性设计相应的成像算法。该成像算法中的关键补偿函数需要知道运动目标的真实速度值,而实际应用中无法提前知道运动目标的真实速度。对此,本发明提出了一种基于二维速度搜索的星载SAR运动目标检测方法。
发明内容
本发明公开了一种基于二维速度搜索的星载SAR运动目标检测方法。本发明解决的技术问题在于:针对运动目标具体状态的未知性,结合RD算法,通过遍历一个间隔均匀的二维数组记为vx(m)、vy(n)(分别代表方位向搜索速度和距离向搜索速度),分别设计距离徙动校正函数和方位向匹配函数,实现运动目标的聚焦成像。然后对方位向滤波器进行修正,消除运动目标距离向速度带来的成像方位向偏移。提取不同搜索速度下的运动目标的最强幅度值用于恒虚警检测,将显著提高运动目标的检测概率。仿真结果验证了该方法的可行性和有效性。
本发明为实现上述发明目的采用如下技术方案:
一种基于二维速度搜索的星载SAR运动目标检测方法,其包括如下步骤:
步骤1,对雷达接收的回波进行距离向脉冲压缩,得到回波s(t,τ),t是方位向的慢时间变量,τ是距离向的快时间变量;
步骤2,将距离向脉冲压缩后的回波s(t,τ)进行方位向FFT得到回波数据s(fa,τ),fa为方位向频率;
步骤3,遍历一个间隔均匀的二维数组记为vx(m)、vy(n),其中vx(m)代表第m个方位向搜索速度,vy(n)代表第n个距离向搜索速度,具体步骤如下:
步骤3-1,根据vx(m)、vy(n)计算得出搜索速度下的多普勒中心频率fdc,多普勒调频率Ka;
步骤3-2,依据fdc、Ka求解距离徙动量Δr,从而对s(fa,τ)进行距离线性插值运算,求出信号在每个τ+2Δr/c时刻的值,并以此替换原τ时刻的数据值,得到校正后的回波数据s′(fa,τ);
步骤3-3,设计方位向匹配滤波器的时域表达式Ha(t),由于此时仍处于多普勒域,所以输出回波数据sout(fa,τ)=s′(fa,τ)·FFT[Ha(t)],再经过方位向IFFT将数据变换回方位向时域,得到匹配后的复图像;
步骤4,根据遍历,重复步骤3-1至3-3,直到整个设置的二维速度搜索完毕,得出m×n幅复图像,提取不同搜索速度下的运动目标的最强幅度值,用于恒虚警检测。
优选地,步骤1中距离向脉冲压缩是完成以下运算:
s(t,τ)=σBr·Wa(t)exp-[j4πr(t)/λ]sinc{πBr[τ-2r(t)/c]},
其中,s(t,τ)为经过距离向脉冲压缩后的点目标回波数据,t是方位向的慢时间变量,τ是距离向的快时间变量,σ是点目标的后向散射系数,Br是距离向带宽,Wa(t)表示雷达天线增益,exp(·)为指数函数,sinc(·)为辛格函数,j表示虚数且有j2=-1,λ表示波长,c是光速,r(t)表示雷达与运动目标的距离。
优选地,步骤3-1中的多普勒中心频率多普勒调频率其中r0为景中心斜距,v为雷达有效速度,y0表示运动目标初始时刻在距离向的位置。
优选地,步骤3-2中依据fdc、Ka求解距离徙动量Δr包括:将r(t)结合fdc、Ka进一步表示为依据方位向慢时间与多普勒频率之间的关系:t=(fa-fdc)/Ka,代入距离徙动量Δr=r(t)-r0中得到多普勒域中距离徙动量表达式:
优选地,步骤3-3包括:经过距离徙动校正后,此时的数据时域表达式变为:
s(t,τ)=σBrWa(t)exp[-j4πr0/λ]exp(jπKat2)exp(j2πfdct)sinc{πBr[τ-2r0/λ]}
静止目标时fdc=0,忽略fdc的影响,设计方位向匹配滤波器的时域表达式为:
其中rect(·)为矩形函数,Ta表示合成孔径时间。
优选地,步骤3-3进一步包括:
修正方位向匹配滤波器表达式:
本发明采用以上技术方案与现有技术相比,具有以下技术效果:由于不知道运动目标真实运动状态,与传统的动目标检测相比,该算法设计了一个间隔均匀的二维数组记为vx(m)、vy(n),分别代表方位向搜索速度和距离向搜索速度,对距离压缩后的场景回波进行距离徙动校正和方位向的匹配滤波,实现对运动目标的聚焦成像。对比运动目标成像图可得出最佳聚焦图,即此时信号幅度是最强的,从而可以通过二维速度搜索提高信噪比。仿真结果验证了该方法的可行性和有效性。
附图说明
图1为正侧视时,星载SAR运动目标几何关系图。
图2为基于二维速度搜索的星载SAR运动目标检测流程图。
图3为对静止目标的成像,其中3(a)为在静止参数下场景成像三维图,3(b)为对应的不考虑噪声情况下的目标成像切片,3(c)为在运动参数下,搜索速度与目标真实速度匹配时场景成像三维,3(d)为对应的不考虑噪声情况下的目标成像切片,3(e)为在运动参数下,搜索速度与目标真实速度失配时的场景成像三维图,3(f)为对应的不考虑噪声情况下的目标成像切片。
图4为仅有方位向速度的运动目标成像,其中4(a)为静止参数下,vx=10m/s时动目标成像,4(b)为对应的不考虑噪声情况下的目标成像切片,4(c)为静止参数下,vx=20m/s时动目标成像,4(d)为对应的不考虑噪声情况下的目标成像切片,4(e)为运动参数下,搜索速度与运动目标真实速度匹配时vx=10m/s的运动目标成像,4(f)为对应不考虑噪声情况下的目标成像切片,4(g)为运动参数下,匹配时vx=20m/s的运动目标成像,4(h)为对应的不考虑噪声情况下的目标成像切片,4(i)为运动参数下,失配时vx=10m/s的运动目标成像,4(j)为对应的不考虑噪声情况下的目标成像切片,4(k)为运动参数下,失配时vx=20m/s的运动目标成像,4(l)为对应的不考虑噪声情况下的目标成像切片。
图5为仅有距离向速度的运动目标成像,其中5(a)为静止参数下,vy=10m/s时运动目标成像,5(b)为对应的不考虑噪声情况下的目标成像切片,5(c)为静止参数下,vy=20m/s时运动目标成像,5(d)为对应的不考虑噪声情况下的目标成像切片,5(e)为运动参数下,搜索速度与运动目标真实速度匹配时vy=10m/s的运动目标成像,5(f)为对应不考虑噪声情况下的目标成像切片,5(g)为运动参数下,匹配时vy=20m/s的运动目标成像,5(h)为对应不考虑噪声情况下的目标成像切片,5(i)为运动参数下,失配时vy=10m/s的运动目标成像,5(j)为对应的不考虑噪声情况下的目标成像切片,5(k)为运动参数下,失配时vy=20m/s的运动目标成像,5(l)为对应不考虑噪声情况下的目标成像切片,5(m)为运动参数下,匹配修正后vy=10m/s的运动目标成像,5(n)为对应的不考虑噪声情况下的目标成像切片,5(o)为运动参数下,匹配修正后vy=20m/s的运动目标成像,5(p)为对应的不考虑噪声情况下的目标成像切片。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细描述,所述实施方式的示例在附图中示出,下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本发明是一种基于二维速度搜索的星载SAR运动目标检测方法,其方法流程如图2所示,包括如下步骤:
步骤1,数据预处理,对雷达接收的回波进行距离向脉冲压缩,即完成以下运算:
s(t,τ)=σBr·Wa(t)exp-[j4πr(t)/λ]sinc{πBr[τ-2r(t)/c]}
其中,s(t,τ)为经过距离向脉冲压缩后的点目标回波数据,t是方位向的慢时间变量,τ是距离向的快时间变量,σ是点目标的后向散射系数,Br是距离向带宽,Wa(t)表示雷达天线增益,exp(·)为指数函数,sinc(·)为辛格函数,j表示虚数且有j2=-1,λ表示波长,c是光速,r(t)表示雷达与运动目标的距离,根据SAR几何关系图可以得到其中v为雷达有效速度,h为雷达高度,x0、y0分别表示运动目标初始时刻在方位向、距离向的位置,vx、vy分别表示运动目标在方位向和速度向的真实运动速度,进一步泰勒展开得出r0为景中心斜距;
此时,运动目标的方位向相位φ(t):
对应的方位向频率为:
令x0=0,则多普勒中心频率多普勒调频率
步骤2,首先进行初始化,将距离向脉冲压缩后的回波s(t,τ)进行方位向FFT得到回波数据s(fa,τ),fa为方位向频率,即此时回波信号处于距离-多普勒域。
步骤3,然后进行二维速度搜索,其中针对运动目标运动状态的未知性,遍历一个间隔均匀的二维数组记为vx(m)、vy(n),其中vx(m)代表第m个方位向搜索速度,vy(n)代表第n个距离向搜索速度,具体步骤如下:
步骤3-1,根据vx(m)、vy(n)计算得出搜索速度下的多普勒中心频率多普勒调频率
步骤3-2,距离徙动校正,即完成下述运算:
将r(t)结合fdc、Ka进一步表示为根据驻定相位原理,可以得出方位向慢时间与多普勒频率之间的关系:t=(fa-fdc)/Ka,代入距离徙动量Δr=r(t)-r0中得到多普勒域中距离徙动量表达式:
从而对s(fa,τ)进行距离线性插值运算,求出信号在每个τ+2Δr/c时刻的值,并以此替换原τ时刻的数据值,使原表达式中sinc{πBr[τ-2(r0+Δr)/c]}变换为sinc{πBr[τ-2r0/c]},将距离徙动曲线拉直到与方位频率平行的方向,从而得到校正后的回波数据s′(fa,τ);
步骤3-3,方位向匹配滤波,即依次完成下述运算:
经过距离徙动校正后,此时的数据时域表达式变为:
s(t,τ)=σBrWa(t)exp[-j4πr0/λ]exp(jπKat2)exp(j2πfdct)sinc{πBr[τ-2r0/λ]}
再通过方位向的匹配滤波处理便可使运动目标分量在方位向聚焦。匹配滤波即对信号进行处理,使其频谱幅度非常平坦,相位仅包含常量和线性分量。因此,可以与含有二次共轭相位的类似频谱信号相乘,相乘后的信号相位即是线性的。静止目标时fdc=0,忽略fdc的影响,设计方位向匹配滤波器的时域表达式Ha(t)为:
其中rect(·)为矩形函数,Ta表示合成孔径时间。由于此时仍处于多普勒域,所以sout(fa,τ)=s′(fa,τ)·FFT[Ha(t)],再经过方位向IFFT将数据变换回方位向时域,得到匹配后的复图像。考虑运动目标fdc的存在,运动目标在聚焦的静止目标SAR图像(常规SAR图像)上的方位位置会偏离其真实方位位置,可以消除运动目标方位向偏移,修正方位向匹配滤波器表达式:
步骤4,图像筛选。根据遍历,重复步骤3-1至3-3,直到整个设置的二维速度搜索完毕,得出m×n幅复图像,提取不同搜索速度下的运动目标的最强幅度值,即此时的搜索速度下检测到的目标信号最强、成像效果最佳,所以通过二位速度搜索可以提高信噪比,用于恒虚警检测,将显著提高运动目标的检测概率。
以上介绍了一种基于二维速度搜索的星载SAR运动目标检测方法,这里利用仿真的星载SAR回波数据进行验证和分析。仿真的星载SAR系统参数如下:在场景中建立1个运动点目标,4个静止点目标(杂波),并引入高斯噪声,使得静止目标成像时的信噪比为30dB。具体的星载SAR系统仿真参数如表1。
表1
由于对静止点目标和运动点目标涉及的Ka、fdc等不同,导致校正及匹配滤波器不同,接下来分别在静止参数下和运动参数下进行对比仿真,分析目标运动给成像带来的影响,具体分为3部分,即对静止目标的成像,对仅有距离向速度的运动目标成像和对仅有方位向速度的运动目标成像。对运动目标的成像,设置速度分别为10m/s和20m/s作为对比分析。由于现实中雷达对于运动目标的具体运动状态是未知的,所以采用速度搜索的方法,设计间隔为5m/s的二维速度记为vx(m)、vy(n),分别代表方位向搜索速度和距离向搜索速度,对收到的运动目标回波进行距离徙动校正和方位向的匹配滤波处理,具体二维速度搜索范围如下:
在二维速度搜索时采用vx(m)=-20m/s,vy(n)=-10m/s作为搜索速度与运动目标真实速度不匹配时的对比。
基于二维速度搜索的运动目标检测方法实验:
首先进行步骤1,数据预处理。对雷达接收的回波进行距离向脉冲压缩。
接着进行步骤2,即采用基于二维速度搜索的动目标检测。对静止点目标的成像如图3所示。其中3(a)表示在静止参数下场景成像三维图,3(b)为其不考虑噪声情况下,目标成像切片。SNR1=30dB,杂波以点目标4为例,SCR1=3.21dB。3(c)表示在运动参数下,搜索速度与目标真实速度匹配时场景成像三维图,3(d)为其不考虑噪声情况下,目标成像切片。SNR2=30dB,SCR2=3.21dB。3(e)表示在运动参数下,搜索速度与目标真实速度失配时场景成像三维图,3(f)为其不考虑噪声情况下,目标成像切片。SNR3=24.06dB,SCR3=0.18dB。则有SNR1=SNR2>SNR3,SCR1=SCR2>SCR3。由于vy=0m/s,vx=0m/s,即Ka、fdc与静止参数下相同,所以在静止参数下对静止目标成像以及在运动参数下速度搜索匹配时对静止目标成像的效果相同,切片呈现标准的十字形,目标聚焦好。而运动参数下,速度搜索不匹配时信号幅值要明显减小,削弱杂波的同时也削弱了信号,对静止点目标的成像最差。
图4给出了仅有方位向速度的运动目标成像。图4(a)、4(c)分别为静止参数下,vx=10m/s、vx=20m/s时动目标成像,计算得出信噪比和信杂比分别为SNR4=28.58dB,SCR4=1.80dB,SNR5=24.05dB,SCR5=-2.74dB。图4(b)、4(d)为其对应的忽略噪声的运动目标成像切片。用静止目标的多普勒斜率对运动目标进行聚焦成像,必然会导致运动目标方位向散焦且幅度下降,结合SNR4(SCR4)>SNR5(SCR5),可得出:在静止参数下,对方位向运动的目标成像会导致散焦,其散焦程度随着多普勒斜率的增大而增大。而杂波没有得到抑制,杂波的幅度与在静止参数下对静止目标的成像的杂波幅度接近。
图4(e)、4(g)为运动参数下,搜索速度与运动目标真实速度vx=10m/s、vx=20m/s分别匹配时的运动目标成像,计算得出信噪比和信杂比:SNR6=30dB,SCR6=4.52dB,SNR7=30dB,SCR7=5.12dB。图4(f)、4(h)为其对应的忽略噪声的运动目标成像切片。运动目标的幅度都与静止参数下对静止目标成像的幅度相同,从而信噪比相同。杂波的幅度明显下降,因而信杂比比静止参数下对静止参数成像的要高。
图4(i)、4(k)给出了运动参数下,搜索速度与运动目标真实速度vx=10m/s、vx=20m/s失配时的运动目标成像,计算得出信噪比和信杂比:SNR8=22.80dB,SCR8=-1.08dB,SNR9=21.24dB,SCR9=-2.64dB。图4(j)、4(l)为其对应的忽略噪声的运动目标成像切片。信噪比和信杂比明显减小。即速度搜索与运动目标速度失配时,方位向的匹配滤波会出现一定的偏差,会严重降低信号的幅值。此时对杂波抑制也有一定的作用。
仅有径向速度的运动目标成像如图5所示。图5(a)、5(c)分别为静止参数下,vy=10m/s、vy=20m/s时动目标成像。图5(e)、5(g)为运动参数下,搜索速度与运动目标真实速度vy=10m/s、vy=20m/s分别匹配时的运动目标成像。图5(i)、5(k)给出了运动参数下,搜索速度与运动目标真实速度失配时的运动目标成像。图5(b)、5(d)、5(i)、5(k)、5(j)、5(l)分别给出了其对应的忽略噪声的运动目标成像切片。由于径向速度对Ka的影响很小,所以用静止参数对距离向运动目标的成像不影响运动目标的聚焦。运动目标距离向速度的存在主要会引起运动目标回波信号的多普勒中心频率的变化,从而使得运动目标在方位向偏离真实位置,偏移量与vy成正比。当搜索速度与运动目标真实速度相匹配时,相比于在静止参数下的成像,信号幅值变大,从而信噪比得到了提高,而失配对运动目标聚焦的影响较大,信号幅值明显下降。
对此,消除运动目标方位向偏移可用修正后的方位向匹配滤波器,对比图5(e)、5(g),图5(m)、5(o)的运动目标回到了方位向中心位置,即方位向第544像素处。
上述仿真结果验证了基于二维速度搜索的星载SAR运动目标检测方法的有效性。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
本技术领域技术人员可以理解的是,本发明可以涉及用于执行本申请中所述操作中的一项或多项操作的设备。所述设备可以为所需的目的而专门设计和制造,或者也可以包括通用计算机中的已知设备,所述通用计算机有存储在其内的程序选择性地激活或重构。这样的计算机程序可以被存储在设备(例如,计算机)可读介质中或者存储在适于存储电子指令并分别耦联到总线的任何类型的介质中,所述计算机可读介质包括但不限于任何类型的盘(包括软盘、硬盘、光盘、CD-ROM、和磁光盘)、随机存储器(RAM)、只读存储器(ROM)、电可编程ROM、电可擦ROM(EPROM)、电可擦除可编程ROM(EEPROM)、闪存、磁性卡片或光线卡片。可读介质包括用于以由设备(例如,计算机)可读的形式存储或传输信息的任何机构。例如,可读介质包括随机存储器(RAM)、只读存储器(ROM)、磁盘存储介质、光学存储介质、闪存装置、以电的、光的、声的或其他的形式传播的信号(例如载波、红外信号、数字信号)等。
本技术领域技术人员可以理解的是,可以用计算机程序指令来实现这些结构图和/或框图和/或流图中的每个框以及这些结构图和/或框图和/或流图中的框的组合。可以将这些计算机程序指令提供给通用计算机、专业计算机或其他可编程数据处理方法的处理器来生成机器,从而通过计算机或其他可编程数据处理方法的处理器来执行的指令创建了用于实现结构图和/或框图和/或流图的框或多个框中指定的方法。
本技术领域技术人员可以理解的是,本发明中已经讨论过的各种操作、方法、流程中的步骤、措施、方案可以被交替、更改、组合或删除。进一步地,具有本发明中已经讨论过的各种操作、方法、流程中的其他步骤、措施、方案也可以被交替、更改、重排、分解、组合或删除。进一步地,现有技术中的具有与本发明中公开的各种操作、方法、流程中的步骤、措施、方案也可以被交替、更改、重排、分解、组合或删除。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (6)
1.一种基于二维速度搜索的星载SAR运动目标检测方法,其特征在于包括如下步骤:
步骤1,对雷达接收的回波进行距离向脉冲压缩,得到回波s(t,τ),t是方位向的慢时间变量,τ是距离向的快时间变量;
步骤2,将距离向脉冲压缩后的回波s(t,τ)进行方位向FFT得到回波数据s(fa,τ),fa为方位向频率;
步骤3,遍历一个间隔均匀的二维数组记为vx(m)、vy(n),其中vx(m)代表第m个方位向搜索速度,vy(n)代表第n个距离向搜索速度,具体步骤如下:
步骤3-1,根据vx(m)、vy(n)计算得出搜索速度下的多普勒中心频率fdc,多普勒调频率Ka;
步骤3-2,依据fdc、Ka求解距离徙动量Δr,从而对s(fa,τ)进行距离线性插值运算,求出信号在每个τ+2Δr/c时刻的值,并以此替换原τ时刻的数据值,得到校正后的回波数据s′(fa,τ);
步骤3-3,设计方位向匹配滤波器的时域表达式Ha(t),由于此时仍处于多普勒域,所以输出回波数据sout(fa,τ)=s′(fa,τ)·FFT[Ha(t)],再经过方位向IFFT将数据变换回方位向时域,得到匹配后的复图像;
步骤4,根据遍历,重复步骤3-1至3-3,直到整个设置的二维速度搜索完毕,得出m×n幅复图像,提取不同搜索速度下的运动目标的最强幅度值,用于恒虚警检测。
2.根据权利要求1所述的基于二维速度搜索的星载SAR运动目标检测方法,其特征在于,所述步骤1中距离向脉冲压缩是完成以下运算:
s(t,τ)=σBr·Wa(t)exp-[j4πr(t)/λ]sinc{πBr[τ-2r(t)/c]},
其中,s(t,τ)为经过距离向脉冲压缩后的点目标回波数据,t是方位向的慢时间变量,τ是距离向的快时间变量,σ是点目标的后向散射系数,Br是距离向带宽,Wa(t)表示雷达天线增益,exp(·)为指数函数,sinc(·)为辛格函数,j表示虚数且有j2=-1,λ表示波长,c是光速,r(t)表示雷达与运动目标的距离。
3.根据权利要求2所述的基于二维速度搜索的星载SAR运动目标检测方法,其特征在于,所述步骤3-1中的多普勒中心频率多普勒调频率其中r0为景中心斜距,v为雷达有效速度,y0表示运动目标初始时刻在距离向的位置。
4.根据权利要求3所述的基于二维速度搜索的星载SAR运动目标检测方法,其特征在于,所述步骤3-2中依据fdc、Ka求解距离徙动量Δr包括:将r(t)结合fdc、Ka进一步表示为依据方位向慢时间与多普勒频率之间的关系:t=(fa-fdc)/Ka,代入距离徙动量Δr=r(t)-r0中得到多普勒域中距离徙动量表达式:
<mrow>
<mi>&Delta;</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<msub>
<mi>v</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
5.根据权利要求4所述的基于二维速度搜索的星载SAR运动目标检测方法,其特征在于,所述步骤3-3包括:经过距离徙动校正后,此时的数据时域表达式变为:
s(t,τ)=σBrWa(t)exp[-j4πr0/λ]exp(jπKat2)exp(j2πfdct)sinc{πBr[τ-2r0/λ]}
静止目标时fdc=0,忽略fdc的影响,设计方位向匹配滤波器的时域表达式为:
<mrow>
<msub>
<mi>H</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>r</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>t</mi>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<msub>
<mi>j&pi;K</mi>
<mi>a</mi>
</msub>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>}</mo>
</mrow>
其中rect(·)为矩形函数,Ta表示合成孔径时间。
6.根据权利要求5所述的基于二维速度搜索的星载SAR运动目标检测方法,其特征在于,所述步骤3-3进一步包括:
修正方位向匹配滤波器表达式:
<mrow>
<msup>
<msub>
<mi>H</mi>
<mi>a</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>r</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>t</mi>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<msub>
<mi>j&pi;K</mi>
<mi>a</mi>
</msub>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>}</mo>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
</mrow>
</msub>
<mi>t</mi>
<mo>}</mo>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810058355.4A CN108051812B (zh) | 2018-01-22 | 2018-01-22 | 基于二维速度搜索的星载sar运动目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810058355.4A CN108051812B (zh) | 2018-01-22 | 2018-01-22 | 基于二维速度搜索的星载sar运动目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108051812A true CN108051812A (zh) | 2018-05-18 |
CN108051812B CN108051812B (zh) | 2021-06-22 |
Family
ID=62127653
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810058355.4A Active CN108051812B (zh) | 2018-01-22 | 2018-01-22 | 基于二维速度搜索的星载sar运动目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108051812B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108549063A (zh) * | 2018-06-29 | 2018-09-18 | 北京航空航天大学 | 一种星载合成孔径雷达图像中运动目标航迹向速度快速估计方法 |
CN108761391A (zh) * | 2018-05-29 | 2018-11-06 | 南京信息工程大学 | 一种模型类无设备目标定位方法 |
CN109358324A (zh) * | 2018-10-31 | 2019-02-19 | 安徽四创电子股份有限公司 | 基于地面监视雷达谱数据的目标检测方法 |
CN110068833A (zh) * | 2019-05-05 | 2019-07-30 | 中国科学院电子学研究所 | 一种合成孔径激光雷达成像方法、仪器及系统 |
CN112710999A (zh) * | 2020-12-17 | 2021-04-27 | 南京航空航天大学 | 基于径向速度搜索的弧形阵列雷达运动目标聚焦检测方法 |
CN114488132A (zh) * | 2022-02-22 | 2022-05-13 | 北京卫星信息工程研究所 | 星载gnss-s雷达的目标检测方法 |
CN114488054A (zh) * | 2021-12-15 | 2022-05-13 | 重庆大学 | 一种计算高效的合成孔径雷达地面运动目标聚焦方法 |
WO2022174630A1 (zh) * | 2021-02-22 | 2022-08-25 | 华为技术有限公司 | 信号处理方法、装置、存储介质及车辆 |
CN115575922A (zh) * | 2022-12-07 | 2023-01-06 | 中国人民解放军国防科技大学 | 基于车载fmcw毫米波雷达的动目标检测方法和装置 |
CN115616629A (zh) * | 2022-09-01 | 2023-01-17 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于天基外辐射源信号的运动目标探测补偿方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5539408A (en) * | 1994-07-08 | 1996-07-23 | Deutsche Forschungsanstalt Fur Luft- Und Raumfahrt E.V. | Method for the detection, localization and velocity determination of moving targets from raw radar data from a coherent, single- or multi-channel image system carried along in a vehicle |
EP1918737A1 (en) * | 2006-11-04 | 2008-05-07 | Sula Systems Limited | Method of resolving ambiguity in the detection and location of moving targets in synthetic aperture radar |
CA2427606C (en) * | 2000-11-02 | 2011-07-05 | Totalfoersvarets Forskningsinstitut | A synthetic aperture radar system capable of detecting moving targets |
US20130154876A1 (en) * | 2010-09-15 | 2013-06-20 | Mitsubishi Space Software Co., Ltd. | Target velocity identification apparatus, target velocity identification program, and target velocity identification method |
US8797206B2 (en) * | 2012-06-13 | 2014-08-05 | C & P Technologies, Inc. | Method and apparatus for simultaneous multi-mode processing performing target detection and tracking using along track interferometry (ATI) and space-time adaptive processing (STAP) |
CN104730505A (zh) * | 2014-12-08 | 2015-06-24 | 广西大学 | 一种多通道sar地面目标检测定位方法和装置 |
CN106872974A (zh) * | 2017-01-23 | 2017-06-20 | 西安电子科技大学 | 基于高超声速平台双通道雷达的高精度运动目标成像方法 |
CN106908781A (zh) * | 2017-02-27 | 2017-06-30 | 中国科学院电子学研究所 | 基于单通道圆迹sar的直线运动动目标的速度矢量获取方法 |
CN107561533A (zh) * | 2017-07-14 | 2018-01-09 | 清华大学 | 一种c波段星载合成孔径雷达运动目标成像方法 |
-
2018
- 2018-01-22 CN CN201810058355.4A patent/CN108051812B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5539408A (en) * | 1994-07-08 | 1996-07-23 | Deutsche Forschungsanstalt Fur Luft- Und Raumfahrt E.V. | Method for the detection, localization and velocity determination of moving targets from raw radar data from a coherent, single- or multi-channel image system carried along in a vehicle |
CA2427606C (en) * | 2000-11-02 | 2011-07-05 | Totalfoersvarets Forskningsinstitut | A synthetic aperture radar system capable of detecting moving targets |
EP1918737A1 (en) * | 2006-11-04 | 2008-05-07 | Sula Systems Limited | Method of resolving ambiguity in the detection and location of moving targets in synthetic aperture radar |
US20130154876A1 (en) * | 2010-09-15 | 2013-06-20 | Mitsubishi Space Software Co., Ltd. | Target velocity identification apparatus, target velocity identification program, and target velocity identification method |
US8797206B2 (en) * | 2012-06-13 | 2014-08-05 | C & P Technologies, Inc. | Method and apparatus for simultaneous multi-mode processing performing target detection and tracking using along track interferometry (ATI) and space-time adaptive processing (STAP) |
CN104730505A (zh) * | 2014-12-08 | 2015-06-24 | 广西大学 | 一种多通道sar地面目标检测定位方法和装置 |
CN106872974A (zh) * | 2017-01-23 | 2017-06-20 | 西安电子科技大学 | 基于高超声速平台双通道雷达的高精度运动目标成像方法 |
CN106908781A (zh) * | 2017-02-27 | 2017-06-30 | 中国科学院电子学研究所 | 基于单通道圆迹sar的直线运动动目标的速度矢量获取方法 |
CN107561533A (zh) * | 2017-07-14 | 2018-01-09 | 清华大学 | 一种c波段星载合成孔径雷达运动目标成像方法 |
Non-Patent Citations (2)
Title |
---|
HE YAN等: "Ground Moving Target Extraction in a Multichannel Wide Area Surveillance SAR_GMTI System via the Relaxed PCP", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
闫贺等: "基于二维速度搜索的星载SAR运动目标聚焦算法研究", 《电子与信息学报》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108761391A (zh) * | 2018-05-29 | 2018-11-06 | 南京信息工程大学 | 一种模型类无设备目标定位方法 |
CN108761391B (zh) * | 2018-05-29 | 2022-04-01 | 南京信息工程大学 | 一种模型类无设备目标定位方法 |
CN108549063A (zh) * | 2018-06-29 | 2018-09-18 | 北京航空航天大学 | 一种星载合成孔径雷达图像中运动目标航迹向速度快速估计方法 |
CN108549063B (zh) * | 2018-06-29 | 2021-11-26 | 北京航空航天大学 | 一种星载合成孔径雷达图像中运动目标航迹向速度快速估计方法 |
CN109358324B (zh) * | 2018-10-31 | 2020-07-14 | 安徽四创电子股份有限公司 | 基于地面监视雷达谱数据的目标检测方法 |
CN109358324A (zh) * | 2018-10-31 | 2019-02-19 | 安徽四创电子股份有限公司 | 基于地面监视雷达谱数据的目标检测方法 |
CN110068833B (zh) * | 2019-05-05 | 2021-10-29 | 中国科学院电子学研究所 | 一种合成孔径激光雷达成像方法、仪器及系统 |
CN110068833A (zh) * | 2019-05-05 | 2019-07-30 | 中国科学院电子学研究所 | 一种合成孔径激光雷达成像方法、仪器及系统 |
CN112710999A (zh) * | 2020-12-17 | 2021-04-27 | 南京航空航天大学 | 基于径向速度搜索的弧形阵列雷达运动目标聚焦检测方法 |
WO2022174630A1 (zh) * | 2021-02-22 | 2022-08-25 | 华为技术有限公司 | 信号处理方法、装置、存储介质及车辆 |
CN114488054A (zh) * | 2021-12-15 | 2022-05-13 | 重庆大学 | 一种计算高效的合成孔径雷达地面运动目标聚焦方法 |
CN114488054B (zh) * | 2021-12-15 | 2024-05-21 | 重庆大学 | 一种计算高效的合成孔径雷达地面运动目标聚焦方法 |
CN114488132A (zh) * | 2022-02-22 | 2022-05-13 | 北京卫星信息工程研究所 | 星载gnss-s雷达的目标检测方法 |
CN114488132B (zh) * | 2022-02-22 | 2023-01-20 | 北京卫星信息工程研究所 | 星载gnss-s雷达的目标检测方法 |
CN115616629A (zh) * | 2022-09-01 | 2023-01-17 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于天基外辐射源信号的运动目标探测补偿方法 |
CN115616629B (zh) * | 2022-09-01 | 2023-08-11 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于天基外辐射源信号的运动目标探测补偿方法 |
CN115575922A (zh) * | 2022-12-07 | 2023-01-06 | 中国人民解放军国防科技大学 | 基于车载fmcw毫米波雷达的动目标检测方法和装置 |
CN115575922B (zh) * | 2022-12-07 | 2023-03-03 | 中国人民解放军国防科技大学 | 基于车载fmcw毫米波雷达的动目标检测方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108051812B (zh) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108051812A (zh) | 基于二维速度搜索的星载sar运动目标检测方法 | |
CN104851097B (zh) | 基于目标形状与阴影辅助的多通道sar‑gmti方法 | |
US9348021B2 (en) | Methods and apparatus for adaptive motion compensation to remove translational movement between a sensor and a target | |
CN109856635A (zh) | 一种csar地面动目标重聚焦成像方法 | |
Saeedi et al. | Synthetic aperture radar imaging using nonlinear frequency modulation signal | |
CN104898119B (zh) | 一种基于相关函数的动目标参数估计方法 | |
Jungang et al. | Low-frequency ultra-wideband synthetic aperture radar ground moving target imaging | |
CN108490441A (zh) | 基于两级滤波的俯冲段大斜视sar子孔径成像空变校正方法 | |
Liang et al. | A high-order phase correction approach for focusing HS-SAR small-aperture data of high-speed moving platforms | |
CN111551934A (zh) | 一种用于无人机载sar成像的运动补偿自聚焦方法与装置 | |
CN105204019A (zh) | 一种基于波数域算法的步进频sar的成像方法 | |
CN102901966B (zh) | 基于解斜和楔石处理的动目标瞬时距离多普勒成像方法 | |
CN108872983A (zh) | 一种弹载sar成像自聚焦方法 | |
CN115616518A (zh) | 宽带雷达微弱目标运动参数估计与高分辨距离像重构方法 | |
Zhang et al. | Bayesian high resolution range profile reconstruction of high-speed moving target from under-sampled data | |
CN106772373B (zh) | 针对任意地面运动目标的sar成像方法 | |
CN112415512B (zh) | 基于进退法和黄金分割法的sar运动目标聚焦方法 | |
Zhang et al. | Squint‐minimised chirp scaling algorithm for bistatic forward‐looking SAR imaging | |
CN106291556A (zh) | 一种基于双曲线模型的大斜视sar成像算法 | |
CN105044700A (zh) | 基于星载高分辨率sar图像聚焦性能的地面高程反演方法 | |
CN107390198A (zh) | 一种高速运动目标的子带相关配准方法 | |
Saeedi et al. | Improved navigation-based motion compensation for LFMCW synthetic aperture radar imaging | |
CN106209715A (zh) | 一种基于时域对消的噪声调幅干扰抑制方法 | |
CN106093941A (zh) | 一种基于fmcw雷达系统的冰层截面成像方法 | |
Wang et al. | Ground moving target imaging based on 2-D velocity search in high resolution 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 |