发明内容
本发明的目的是针对背景技术存在的缺陷,研究设计一种BFSAR地面动目标检测与成像方法,填补BFSAR在地面动目标检测与成像领域的空白,克服BFSAR在静止背景下难以实现地面动目标检测与成像的问题。
本发明的技术方案为:一种双基地前视SAR动目标检测方法,具体包括如下步骤:
步骤一:成像系统参数初始化
本发明中BFSAR发射站初始位置坐标记为(0,yT,zT),其中0、yT和zT分别为发射站的x轴、y轴和z轴坐标。接收站初始位置坐标记为(x0,yR,zR),其中x0、yR和zR分别为接收站的x轴、少轴和z轴坐标。发射站速度为vT,接收站速度为vR,并均沿着y轴飞行。初始刻接收站和发射站波束中心同时指向参考目标Pt,其坐标为(x0,0,0),设P(x,y)为成像区域中的任意点目标,若P(x,y)为地面动目标,则设其运动速度沿x轴和y轴的分量分别为vx和vy。
目标P(x,y)到接收站和发射站的距离历史分别为
其中,rR为接收站到目标P(x,y)的最短斜距,rT为发射站到目标P(x,y)的最短斜距。η为方位向时间。ηR为接收站的零多普勒时间,ηT为发射站的零多普勒时间,且ηR=(y-yR)/VR,ηT=(y-yT)/VT。
若P(x,y)为地面动目标,则rR为:
其中,
为初始时刻接收站到地面动目标P(x,y)的最短斜距。
rT为:
其中,
为初始时刻接收站到目标P(x,y)的最短斜距。
步骤二:获取BFSAR回波,并对回波进行距离向傅里叶变换
设BFSAR发射信号为脉冲线性调频信号:
其中,f0为载频,Tp为脉冲宽度,τ为距离向时间,Kr为发射信号的时间调频斜率,rect[·]代表距离时间窗;
任意点目标P(x,y)回波信号经解调至基带后,可表示为方位向时间η和距离向时间τ的表达式,记为S(η,τ;x,y),
其中,σ(x,y)为目标P(x,y)的后向散射系数,c为光速,λ为发射信号载波波长,ω[·]为方位时间窗,ηcb为波束中心时刻。
对S(η,τ;x,y)进行距离向傅里叶变换,变换后的数据记为S(η,f;x,y),
其中,St(f)为发射信号St(τ)的频谱,f为距离向频率。
步骤三:Bulk-Deramp滤波
首先,将双基距离历史RR(η)+RT(η)在合成波束中心时刻ηcb沿泰勒级数展开:
其中,Rbi0为发射站与接收站在ηcb时刻的双基距离总和;R′bi,R″bi分别表示ηcb时刻双基距离历史的一阶导数和二阶导数。
将以上泰勒展开代入步骤二中的S(η,f;x,y),可得:
由于接收站工作于前视模式,一阶系数R′bi很大,从而导致多普勒模糊,距离向与方位向的非正交性很强,为了解决以上问题,构造Bulk-Deramp滤波函数:
其中,N为参考目标多普勒调频率的多普勒模糊数,
PRF为脉冲重复频率,round[·]为取整运算。θ
sRref,θ
sTref分别为关于参考目标点接收站与发射站的斜视角。
将本步骤中得到的S(η,f;x,y)乘以HDeramp,滤波后的数据记为S1(η,f;x,y),
其中,R′bires=R′bi-N·PRF·λ为双基距离和的残余一阶系数。
从而通过本步骤Bulk-Deramp滤波,去除了BFSAR回波的多普勒模糊,降低了距离向与方位向的非正交性。
步骤四:一阶Keystone变换
为了消除Bulk-Deramp滤波后的残余一阶耦合,采用一阶Keystone变换,变换关系为:
其中,ηk为一阶Keystone变换后新的方位向时间。变换后的数据记为S2(ηk,f;x,y),
由上式可见,线性距离走动耦合效应已被消除。
步骤五:距离压缩并进行距离向傅里叶逆变换
利用距离压缩函数对步骤四处理后的数据S2(ηk,f;x,y)进行距离压缩,其中,距离压缩函数为(f),上标*代表复共轭。经距离压缩后,回波数据通过距离向傅里叶逆变换到二维时域,其数据记为S3(ηk,τ;x,y),
其中,IFFTrange[·]表示距离向傅里叶逆变换。
从上式可以看出,上述操作后目标P(x,y)位于距离向τ=(Rbi0-R′biηcb)/c单元内。
步骤六:扩展方位非线性调频变标(EA-NLCS)操作
令步骤三中双基距离历史RR(η)+RT(η)泰勒展开式中η=0,同时忽略二阶与高阶项,得:RR(0)+RT(0)≈Rsi0-R′biηcb。从而可以得到,通过一阶Keystone变换后,将在方位零时刻具有相同双基距离和的目标移动到了同一距离单元。
为了利用静止目标与地面动目标不同的多普勒调频率将地面动目标回波从静止场景回波中区别出来,需利用扩展方位非线性调频变标(EA-NLCS)操作均衡同一距离单元静止目标的多普勒调频率,并使地面动目标的多普勒调频率不同。为了得到EA-NLCS的调频变标函数,将进行如下操作:
首先对于同一距离单元R的目标,可以得到目标的多普勒调频率为fηr(R;x,ycb)=fηrR(R;x,ycb)+fηrT(R;x,ycb),其中
其中,θ
SR,θ
ST分别为在波束中心时刻接收站和发射站的斜视角,且
分别为接收站和发射站的下视角,且
再计算出多普勒调频率差分为:
其中xref和ycbref分别为参考目标的x坐标与y坐标。
对Δfηr沿方位向时间进行两次积分得到EA-NLCS调频变标函数的相位φs(ηk;R),从而得出EA-NLCS变标函数SS(ηk;R),
将S3(ηk,τ;x,y)与SS(ηk;R)相乘,实现EA-NLCS操作,从而完成了同一距离单元静止目标多普勒调频率的均衡,并且使同一距离单元静止目标相同的多普勒调频率与动目标的调频率不同。
步骤七:构造二阶模糊函数积(PSAF)
将EA-NLCS处理后的回波通过方位向FFT变换到方位频域、距离时域,其数据记为S4(fηk,τ),
其中,FFTazimuth[·]表示方位向傅里叶变换,下标st代表静止场景,Sst(fηk,τ)表示静止场景回波,SGMT(fηk,τ)表示地面动目标回波,fηk为新的方位频率变量,σst(x,y)为静止目标的后向散射系数,Kr(R,0)为均衡后的静止目标多普勒调频率,KGMT为地面动目标多普勒调频率。
计算出瞬时二阶差分频率积(SDFP)Ssdfp(fηk;fk),计算表达式为:
其中,τk为某一距离时间,fk为某一方位频率值。
二阶模糊函数(SAF)定义为SDFP的傅里叶变换,表达式为:
其中,fT为傅里叶变换后的频率值。
若地面动目标的距离单元仅存在一个静止目标,则对应的二阶模糊函数在fT1=1/Kr(R,0)和fT2=1/KGMT两处存在峰值。然而当输入信号由若干目标共同叠加时,二阶模糊函数存在虚假峰值。为了抑制模糊,利用M组不同的fk构造如下二阶模糊函数积(PSAF):
经过此处理,有用的峰值保留在相同的位置,而虚假的峰值则沿频率轴移动,相乘实现虚假峰值的削弱,因此,经过相乘处理后有用峰值得到了很大提升。
通过如上所述的处理后,BFSAR某一距离单元回波数据的PSAF若存在两个尖峰,则判定此距离单元有地面动目标,从而完成了BFSAR地面动目标的检测。
为了解决上述问题,基于上述BFSAR地面动目标检测方法,本发明还提出了一种BFSAR地面动目标成像方法,在上述方法步骤的基础上,还包括如下步骤:
步骤八:地面动目标回波方位聚焦成像
利用PSAF估计出的地面动目标多普勒调频率
对地面动目标回波进行方位聚焦,就可以实现地面动目标的成像,其方位聚焦函数为S
AC(R;f
ηk):
本发明的有益效果:本发明的检测方法利用静止背景和地面动目标之间不同的多普勒调频率,解决BFSAR在静止背景下地面动目标的检测难的问题。它的特点是首先利用Bulk-Deramp滤波消除多普勒模糊,降低多普勒调频率的空变,接着利用一阶Keystone变换完成静止目标和地面动目标的距离徙动校正,然后利用扩展方位非线性调频变标(EA-NLCS)操作均衡静止目标的多普勒调频率,同时使地面动目标的多普勒调频率与静止目标调频率不同,最后构造二阶乘积模糊函数(PSAF),完成地面动目标的检测并估计出地面动目标的多普勒调频率,解决了BFSAR在强杂波背景下动目标与静止目标难区分的问题;本发明的成像方法是在完成地面目标检测后,利用估计出的地面动目标多普勒调频率,对动目标回波实现聚焦,从而完成了地面动目标的成像。
具体实施方式
本发明主要采用仿真实验的方式进行验证,仿真验证平台为Matlab2012。下面结合附图和具体实施方式对本发明作进一步的详细描述。
本发明的BFSAR动目标检测方法的流程示意图如图1所示,具体过程如下:
步骤一:成像系统参数初始化
本发明具体实施例采用的BFSAR几何结构图如图2所示,采用的BFSAR系统参数表如图3所示,本发明实施采用的目标场景图如图4所示,图中P点为地面动目标,P1,P2,P3,P4均为静止目标。BFSAR接收站平台初始位置坐标为(8,-10,8)km,发射站平台初始位置坐标为(0,-6,10)km。目标中心位置坐标为(8,0,0)km。发射站速度为150m/s,接收站速度为300m/s,地面动目标的速度vx=14m/s,vy=10m/s。
步骤二:获取BFSAR回波,并对回波进行距离向傅里叶变换
按照步骤一中的参数设置,用Matlab仿真出场景目标回波,并对回波数据做距离向傅里叶变换,变换后的数据记为S(η,f;x,y)。经步骤二获取的回波的二维频谱图像如图5所示。
步骤三:Bulk-Deramp滤波
为了去除BFSAR回波的多普勒模糊,降低了距离向与方位向的非正交性,对步骤二得到的数据S(η,f;x,y)进行Bulk-Deramp滤波,Bulk-Deramp滤波函数为:
其中,N为参考目标多普勒调频率的多普勒模糊数,PRF为脉冲重复频率,round[·]为取整运算。滤波后的数据记为S1(η,f;x,y)。经步骤三Bulk-Deramp滤波后的二维频谱图像如图6所示。
步骤四:一阶Keystone变换
为了消除Bulk-Deramp滤波后的残余一阶耦合,对步骤三获取的数据S1(η,f;x,y)进行一阶Keystone变换。
具体变换关系为:
其中,ηk为Keystone变换后新的方位向时间。一阶Keystone变换后的数据记为S2(ηk,f;x,y),经步骤四一阶Keystone变换后的二维频谱图像如图7所示。
步骤五:对步骤四中经一阶Keystone变换后的回波进行距离压缩并进行距离向傅里叶逆变换,变换后的数据记为S3(ηk,τ;x,y),其效果图如图8所示。
由图8可见,静止目标与运动目标的距离徙动均己被去除,并且初始时刻具有相同双基距离和的目标点P1、P、P4结果一阶Keystone变换后落入了同一距离单元。
步骤六:扩展方位非线性调频变标(EA-NLCS)操作
对步骤五获取的数据S3(ηk,τ;x,y)进行扩展方位非线性调频变标(EA-NLCS)操作,EA-NLCS变标函数SS(ηk;R),从而完成了同一距离单元静止目标多普勒调频率的均衡,并且使同一距离单元静止目标相同的多普勒调频率与动目标的调频率不同。
步骤七:构造二阶模糊函数积(PSAF)
将EA-NLCS处理后的回波变换到(fηk,τ)域,其数据记为S4(fηk,τ),计算出瞬时二阶差分频率积(SDFP)Ssdfp(fηk;fk),计算表达式为:
Sfdfp(fηk)=S4(fηk;τk)
Ssdfp(fηk;fk)=Sfdfp(fηk+fk)Sfdfp *(fηk-fk)
其中,τk为某一距离时间,fk为某一方位频率值,这里的“某一”可以根据实际情况进行选定,对此,本领域的技术人员是可以理解的。
二阶模糊函数(SAF)定义为SDFP的傅里叶变换,表达式为:
若地面动目标的距离单元仅存在一个静止目标,则对应的二阶模糊函数在fT1=1/Kr(R,0)和fT2=1/KGMT两处存在峰值。然而当输入信号由若干目标共同叠加时,二阶模糊函数存在虚假峰值。为了抑制模糊,利用M组不同的fk构造如下二阶模糊函数积(PSAF):
经过此处理,有用的峰值保留在相同的位置,而虚假的峰值则沿频率轴移动,相乘实现虚假峰值的削弱,因此,经过相乘处理后有用峰值得到了很大提升。
通过如上所述的处理后,BFSAR某一距离单元回波数据的PSAF若存在两个尖峰,则判定此距离单元有地面动目标,如图9所示,从而完成了BFSAR地面动目标的检测;图10是经步骤七获得的静止目标P2距离单元内的二阶模糊函数积示意图。
在上述检测方法的基础上,本实施例中提供的BFSAR地面动目标成像方法在上述检测方法步骤的基础之上还包括如下步骤:
步骤八:地面动目标回波方位聚焦成像
利用PSAF估计出的地面动目标多普勒调频率
对地面动目标回波进行方位聚焦,就可以实现地面动目标的成像,具体方位聚焦函数为S
AC(R;f
ηk),
图11是本实施例中采用本发明的方法得到的地面动目标成像结果方位向剖面图,其中中间聚焦点为地面动目标图像。通过本发明具体实施方式可以看出,本发明解决了BFSAR在静止背景下难以实现地面动目标检测与成像的问题。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。