CN102798861A - 一种基于最优图像空间双基地合成孔径雷达成像方法 - Google Patents

一种基于最优图像空间双基地合成孔径雷达成像方法 Download PDF

Info

Publication number
CN102798861A
CN102798861A CN2012102510527A CN201210251052A CN102798861A CN 102798861 A CN102798861 A CN 102798861A CN 2012102510527 A CN2012102510527 A CN 2012102510527A CN 201210251052 A CN201210251052 A CN 201210251052A CN 102798861 A CN102798861 A CN 102798861A
Authority
CN
China
Prior art keywords
mrow
msub
mfrac
math
mover
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
Application number
CN2012102510527A
Other languages
English (en)
Other versions
CN102798861B (zh
Inventor
张晓玲
戴春杨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201210251052.7A priority Critical patent/CN102798861B/zh
Publication of CN102798861A publication Critical patent/CN102798861A/zh
Application granted granted Critical
Publication of CN102798861B publication Critical patent/CN102798861B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于最优图像空间双基地合成孔径雷达成像方法,对于移变双基地SAR,由于其相对复杂的空间几何构型以及距离历史的二维空变特性,目前还没有统一的双基地SAR的投影关系;本发明通过计算得到最优投影向量,并利用二维非均匀快速傅里叶变换将双基地SAR回波信号映射至该投影空间中,进行聚焦处理,得到均方意义下的最优成像结果。本发明方法与现有的双基地SAR成像方法法相比,其优势在于能够获得统计意义上的最优聚焦性能。

Description

一种基于最优图像空间双基地合成孔径雷达成像方法
技术领域
本发明属于雷达技术领域,它特别涉及双基地合成孔径雷达(简称双基地SAR)的成像技术。
背景技术
双基地SAR就是指雷达发射机与接收机分置于两个不同运动平台的SAR系统。由于双基地SAR系统发射平台与接收平台分置,与传统单站SAR相比,平台关系更加灵活,具有更加丰富的空间、运动资源。不仅可以实现传统单站SAR的成像功能,而且具有以下几个方面的优势:
1、成像区域灵活。传统单站SAR的成像区域通常处于平台运动方向的侧视方向;而双基地SAR藉由收发平台之间的协作,与成像区域形成不同的空间位置关系,具备更加灵活的成像区域,包括相对于接收机位置的侧视、前视、下视以及后视方向的目标区域成像。
2、目标信息丰富。传统单站SAR只能获取目标区域的后向散射信息。而在双基地SAR系统中,由于收发平台分置,因此可以选取适宜的双基地角去获取更加丰富的目标信息,从而有利于对隐身目标的探测。
3、成本低。双基地SAR的接收机不含大功率微波器件,功耗和体积较小,重量轻,便于小型无人侦察机携带,从而大大降低了军事侦察的成本。
然而,双基地SAR带来的种种好处,是以雷达系统的复杂性为代价的。其发射机与接收机分置于不同运动平台的特点,给成像处理带来了一些具有研究意义的新课题,主要包括:
1、双基地SAR的距离历史为发射机、接收机距离历史之和,具有双根号形式,直接用驻定相位原理无法求解双基地SAR二维频谱的解析表达式,给双基地SAR成像算法的研究带来困难。
2、双基地SAR的移变特性,使点目标的距离徙动具有二维空变性,难以设计出高效而实用的成像算法。
雷达成像的本质是从场景空间到投影空间的映射。在传统单站SAR的成像处理中,通常以“最近斜距—零多普勒时间”(Closest range—Zero Doppler time)做为图像坐标。这是因为,单站SAR的距离历史在该投影空间下具有一维空变性。传统单站SAR的成像算法,包括距离多普勒算法、CS算法和Omega-K算法等,都是利用了这个特性,通过校正一个点目标的距离徙动达到校正一簇点目标距离徙动的效果,从而实现快速成像处理。然而,移变双基地SAR,由于其相对复杂的空间几何构型以及距离历史的二维空变特性,还没有统一的双基地SAR的图像坐标系。截止到目前,对于双基地SAR投影空间的选择问题,还没有做过专门的研究。
发明内容
本发明的目的是为克服现有的双基地SAR成像方法没有统一投影空间的不足,提供了一种基于最优图像空间双基地合成孔径雷达成像方法;该方法充分考虑了成像区域内各点目标距离历史的近似误差,能够保证得到均方意义下的最优聚焦效果。
为了方便描述本发明的内容,首先作以下术语定义:
定义1、双基地SAR系统相关参数描述
双基地SAR距离历史为
R ( u ; P ) = | | V T | | 2 · u 2 + 2 V T T ( P 0 T - P ) · u + | | P 0 T - P | | 2
+ | | V R | | 2 · u 2 + 2 V R T ( P 0 R - P ) · u + | | P 0 R - P | | 2
其中,u表示慢时间,P表示成像区域内一点目标的空间位置矢量,POT和POR分别为发射机和接收机在慢时刻零时的空间位置矢量;VT和VR分别为发射机和接收机的运动速度矢量。运算符T表示矢量的转置,运算符‖·‖表示矢量的范数,关于矢量运算符详细说明请参见书籍《线性代数与空间解析几何(第3版)》,黄廷祝编著,高等教育出版社。
定义2、快速傅里叶变换
计算离散傅里叶变换的快速方法。目前世界上公认最快的快速傅里叶变换计算程序由FFTW提供,详细信息可参考http://www.fftw.org/
定义3、特征值分解
特征值分解是线性代数里常用的一个计算方法,可以通过求解线性方程组来计算得到,详细计算方法请参见书籍《线性代数与空间解析几何(第3版)》,黄廷祝编著,高等教育出版社。
定义4、二维非均匀快速傅里叶变换
二维非均匀快速傅里叶变换定义为
I k , l = Σ n = 0 N - 1 Σ m = 0 M - 1 x m , n exp ( j 2 π M α m , n k + j 2 π N β m , n l )
其中,M和N分别为二维信号x的第一维和第二维点数;m为二维信号x的第一维指标,且m=0,1,…,M-1,n为二维信号x的第二维指标,且n=0,1,…,N-1;αm,n和βm,n为非均匀采样的二维信号xm,n在二维空间中的分布;
Figure BDA00001908476300032
表示复数单位。I为二维非均匀傅里叶变换的输出信号,k为二维信号I的第一维指标,且k=0,1,…,M-1,l为二维信号I的第二维指标,且l=0,1,…,N-1。二维非均匀快速傅里叶变换可以通过变换核函数插值以及传统快速傅里叶变换技术来实现。
本发明提供一种星机联合双基地合成孔径雷达频域成像方法,如图1所示,包含如下步骤:
步骤1、双基地SAR成像参数初始化
初始化成像参数包括:
快时间,记做t;慢时间,记做u;快频率,即快时间经过傅里叶变换所对应的频率,记做f;雷达工作频率,记做fc;脉冲重复频率,记做PRF;电磁波传播速度,记做c;双基地SAR发射信号的带宽和脉宽,分别记做B和Tp;以上参数由雷达的系统参数给出。
慢时间位于时间坐标原点时,发射平台和接收平台的空间位置矢量,分别记做P0T和POR;发射平台与接收平台的运动速度矢量,分别记做VT和VR。以上参数由雷达平台的运动参数给出。
观测区域,记做Ω;观测区域内共有K个点目标,其空间位置矢量分别P1、P2、…、Pk、…、PK,对应的距离历史,分别记做R(u;P1)、R(u;P2)、…、R(u;Pk)、…、R(u;PK)。以上参数由雷达系统的波束相关信息给出。
双基地SAR原始时域回波共包含Na个脉冲,每个脉冲包含Nr个复数浮点数据,并以一个Nr行Na列的数据矩阵存放(Nr和Nr均为正整数),记做s(t,u)。Nr和Na由雷达系统的数据采集设备给出。
定义参考距离Rref(u),记做
R ref ( u ) = Σ k = 1 K R ( u ; P k )
其中,u表示慢时间;K表示成像区域内点目标个数;k为自然数,且k=1,2,…,K;Pk表示成像区域内的第k个点目标的空间位置矢量;R(u;Pk)表示成像区域内的第k个点目标的距离历史,详细信息请参见定义一。
定义参考函数Sref(f,u),记做
S ref ( f , u ) = exp ( jπ T p B f 2 + j 2 π f c + f c Σ k = 1 K R ( u ; P k ) )
其中,f表示快频率;u表示慢时间;
Figure BDA00001908476300043
表示复数单位;B和Tp分别表示发射信号的带宽和脉宽;fc表示雷达工作频率;c表示电磁波传播速度;K表示成像区域内点目标个数;k为自然数,且k=1,2,…,K;Pk表示成像区域内的第k个点目标的空间位置矢量;R(u;Pk)表示成像区域内的第k个点目标的距离历史,详细信息请参见定义一。
步骤2、快时间傅里叶变换
对原始时域回波s(t,u)的每一列数据,进行快速傅里叶变换,得到快时间傅里叶变换的结果同样为Nr行Na列的回波信号S(f,u)。其中,t表示快时间,f表示快频率;u表示慢时间;数据的行数Nr和列数Na的大小由雷达系统的数据采集设备给出。
步骤3、参考函数相乘
用信号S(f,u)与参考函数Sref(f,u)相乘,得到参考函数相乘之后的结果SRFM(f,u),即
SRFM(f,u)=S(f,u)·Sref(f,u)
其中,f表示快频率;u表示慢时间;参考函数Sref(f,u)由步骤1定义。
步骤4、计算最优投影基向量
步骤4.1)根据下式计算矩阵A(m,n),矩阵A(m,n)为Na行Na列的方阵,Na表示数据SRFM(f,u)的列数,Na的值由雷达系统的数据采集设备给出。m为自然数,且m=0,1,…,Na-1;n为自然数,且n=0,1,…,Na-1。
A ( m , n ) = Σ k = 1 K [ R ( u m ; P k ) - R ref ( u m ) ] [ R ( u n ; P k ) - R ref ( u n ) ]
其中,
u m = m + N a / 2 PRF
u n = n + N a / 2 PRF
其中,m为自然数,且m=0,1,…,Na-1;n为自然数,且n=0,1,…,Na-1。Na表示数据SRFM(f,u)的列数,由雷达系统的数据采集设备给出;PRF表示雷达系统的脉冲重复频率。K表示成像区域内点目标个数;k为自然数,且k=1,2,…,K;Pk表示成像区域内的第k个点目标的空间位置矢量;R(u;Pk)表示成像区域内的第k个点目标的距离历史,详细信息请参见定义一;Rref(u)表示参考函数,由步骤一所定义。
步骤4.2)采用传统常用的特征值分解方法对矩阵A(m,n)进行特征值分解,得到Na个特征值σ1、σ2、…、
Figure BDA00001908476300054
以及与特征值对应的特征向量α1(n)、α2(n)、…、
Figure BDA00001908476300055
Na表示数据SRFM(f,u)的列数,由雷达系统的数据采集设备决定;n为自然数,且n=0,1,…,Na-1。
步骤4.3)对Na个特征值σ1、σ2、…、
Figure BDA00001908476300056
按照由大到小的顺序重新排序,得到排序后的Na个特征值
Figure BDA00001908476300062
...、
Figure BDA00001908476300063
以及与特征值对应的特征向量
Figure BDA00001908476300065
...、
Figure BDA00001908476300066
其中
Figure BDA00001908476300067
Figure BDA00001908476300068
即为双基地SAR的最优投影基向量。
步骤五、二维非均匀快速傅里叶变换
按照下式对参考函数相乘的结果SRFM(f,u)进行二维非均匀快速傅里叶变换,得到双基地SAR图像I(k,l)。
I ( k , l ) = Σ n = - N a / 2 N a / 2 - 1 Σ m = - N r / 2 N r / 2 - 1 S RFM ( mΔf , nΔu ) · exp j 2 π f c + mΔf c α ~ 1 ( n ) k + j 2 π f c + mΔf c α ~ 2 ( n ) l
其中,Δf=B/Nr;Δu=1/PRF;
Figure BDA000019084763000610
为双基地SAR的最优投影基向量;B表示发射信号的带宽;PRF表示雷达系统的脉冲重复频率。m和n为自然数,且m=-Nr/2,-Nr/2+1,…,Nr/2-1,n=-Na/2,Na/2+1,…,Na/2-1;fc表示雷达系统的工作频率;c表示电磁波传播速度。k和l为自然数,且m=-Nr/2,-Nr/2+1,…,Nr/2-1,n=-Na/2,Na/2+1,…,Na/2-1。该步骤的具体实现如下:
步骤5.1)根据下式计算得到二维信号Y(g,h)
Y ( g , h ) = Σ g = round ( 2 N r ( f c + mΔf ) α ~ 1 ( n ) c ) + p Σ h = round ( 2 N a ( f c + mΔf ) α ~ 2 ( n ) c ) + q ρ ( p , N r ( f c + mΔf ) α ~ 1 ( n ) c ) × ρ ( q , N a ( f c + mΔf ) α ~ 2 ( n ) c ) × S RFM ( mΔf , nΔu )
其中,函数ρ(q,c)为
ρ ( q , c ) = 1 2 π exp { - 1 4 [ round ( 2 c ) + q - 2 c ] 2 }
其中,round(·)表示四舍五入函数;g为整数,且g=-Nr-Q,-Nr-Q+1,…,Nr+Q;h为整数,且h=-Na-Q,-Na-Q+1,…,Na+Q;行数Nr和列数Na由雷达系统的数据采集设备决定;m和n为自然数,且m=-Nr/2,-Nr/2+1,…,Nr/2-1,n=-Na/2,Na/2+1,…,Na/2-1;fc表示雷达系统的工作频率;c表示电磁波传播速度;p和q为整数,且p=-Q,-Q+1,…,Q,q=-Q,-Q+1,…,Q;
Figure BDA00001908476300071
Figure BDA00001908476300072
为双基地SAR的最优投影基向量;
步骤5.2)利用快速傅里叶变换,根据下式计算得到二维信号
Figure BDA00001908476300073
I ~ ( k , l ) = Σ g = - N r - Q N r + Q Σ h = - N a - Q N a + Q Y ( g , h ) exp ( j π N r gk + j π N a hl )
其中,g为整数,且g=-Nr-Q,-Nr-Q+1,…,Nr+Q;h为整数,且h=-Na-Q,-Na-Q+1,…,Na+Q;行数Nr和列数Na由雷达系统的数据采集设备决定;k为整数,且k=-Nr/2,-Nr/2+1,…,Nr/2-1;l为整数,且l=-Na/2,-Na/2+1,…,Na/2-1。
步骤5.3)根据下式,计算得到最终的聚焦图像I(k,l)
I ( k , l ) = exp [ π 2 ( k 2 N r 2 + l 2 N a 2 ) ] · I ~ ( k , l )
其中,k为整数,且k=-Nr/2,-Nr/2+1,…,Nr/2-1;l为整数,且l=-Na/2,-Na/2+1,…,Na/2-1。行数Nr和列数Na的大小由雷达系统的数据采集设备给出。
经过上述步骤处理,就从双基地SAR接收到的目标回波信号s(t,u)中获取具有高分辨率的目标成像结果聚焦图像I(k,l)。
需要说明的是,步骤4中的计算最优投影基向量过程,也可以根据双基地SAR的系统参数和运动参数进行离线计算;从而在实际成像处理时,经过步骤2的快时间傅里叶变换和步骤3的参考函数相乘,即能直接根据离线计算的最优投影基向量。
本发明的实质是,通过计算双基地SAR的最优投影基向量,并利用二维非均匀快速傅里叶变换对双基地SAR回波信号进行聚焦处理。该方法与现有的双基地SAR成像算法相比,其优势在于能够获得统计意义上的最优聚焦性能。
本发明的创新点在于,通过计算得到的最优投影向量,能够保证双基地SAR回波信号在均方意义下的最优聚焦结果。
本发明的基本原理是,通过计算得到最优投影向量,并利用二维非均匀快速傅里叶变换实现双基地SAR回波信号的聚焦处理。
本发明解决的技术问题是,移变双基地SAR,由于其相对复杂的空间几何构型以及距离历史的二维空变特性,还没有统一的双基地SAR的投影关系;本发明通过计算得到最优投影向量,并利用二维非均匀快速傅里叶变换将双基地SAR回波信号映射至该投影空间中,得到均方意义下的最优成像结果。
本发明的有益效果,通过计算得到双基地SAR的最优投影向量,保证本发明能够获得均方意义下的最优成像结果。
附图说明
图1为本发明实现步骤流程图。
图2为仿真中的点目标场景。X轴和Y轴分别表示地面坐标系上的两个坐标轴。O点为成像区域的中心点;A、B、C、D四个点目标分别位于1km×1km正方形区域内的四个顶点。
图3为点目标成像结果。V1轴和V2轴分别为双基地SAR图像的两个坐标轴,O’A’B’C’D’五个点目标图像分别为图2中五个点目标OABCD的成像结果。
具体实施方式
本发明主要采用仿真实验的方法进行验证,所有步骤、结论都在MATLAB上验证正确。
步骤1、双基地SAR成像参数初始化。
本实施例分别采用TerraSAR-X卫星和机载PAMIR作为发射、接收平台,发射平台的运动速度矢量为[7600;0;0;](m/s),在慢时刻零的空间位置矢量为[-300000;-500000;514000;](m);接收平台的运动速度矢量为[0;-5000;6000;](m/s),在慢时刻零的空间位置矢量为[100;0;0;](m)。发射信号的中心频率为10GHz,发射信号带宽为150MHz,合成孔径时间为1.024s。仿真成像场景内包含五个点目标,其中参考点目标位于中心,其余四个点目标分别在斜距向和方位向上距该参考点目标500米。
步骤2、快时间傅里叶变换
原始时域回波s(t,u)以一个1024行2048列的数据矩阵存放,其中每行数据是存放慢时间(方位向)回波信号的采样;每列的数据是存放快频率(斜距向)单脉冲回波信号的采样。对于中s(t,u)的每一列数据,进行快速傅里叶变换,得到回波信号S(f,u)。其中,t表示快时间,f表示快频率;u表示慢时间。
步骤3、参考函数相乘。
用回波信号S(f,u)与参考函数Sref(f,u)相乘,得到参考函数相乘之后的结果SRFM(f,u),即
SRFM(f,u)=S(f,u)·Sref(f,u)
其中,t表示快时间,f表示快频率;u表示慢时间。
步骤4、计算最优投影基向量
在成像区域内,根据5个点目标的位置,求解得到双基地SAR的最优投影基向量
Figure BDA00001908476300092
步骤5、二维非均匀快速傅里叶变换
根据双基地SAR的最优投影基向量
Figure BDA00001908476300093
Figure BDA00001908476300094
对参考函数相乘之后的回波信号进行二维非均匀快速傅里叶变换。
经过上述步骤,即能获得基于最优投影空间的双基地SAR成像结果。利用本发明,聚焦结果为图3;由该图可知5个点目标都被聚焦在正确的位置上。

Claims (1)

1.一种基于最优图像空间双基地合成孔径雷达成像方法,其特征是它包含如下步骤:
步骤1、双基地SAR成像参数初始化
初始化成像参数包括:
快时间,记做t;慢时间,记做u;快频率,即快时间经过傅里叶变换所对应的频率,记做f;雷达工作频率,记做fc;脉冲重复频率,记做PRF;电磁波传播速度,记做c;双基地SAR发射信号的带宽和脉宽,分别记做B和Tp;以上参数由雷达的系统参数给出;
慢时间位于时间坐标原点时,发射平台和接收平台的空间位置矢量,分别记做P0T和POR;发射平台与接收平台的运动速度矢量,分别记做VT和VR;以上参数由雷达平台的运动参数给出;
观测区域,记做Ω;观测区域内共有K个点目标,其空间位置矢量分别P1、P2、…、Pk、…、PK,对应的距离历史,分别记做R(u;P1)、R(u;P2)、…、R(u;Pk)、…、R(u;PK);以上参数由雷达系统的波束相关信息给出;
双基地SAR原始时域回波共包含Na个脉冲,每个脉冲包含Nr个复数浮点数据,并以一个Nr行Na列的数据矩阵存放,记做s(t,u);Nr和Nr均为正整数Nr和Na的大小由雷达系统的数据采集设备给出;
定义参考距离Rref(u),记做
R ref ( u ) = Σ k = 1 K R ( u ; P k )
其中,u表示慢时间;K表示成像区域内点目标个数;k为自然数,且k=1,2,…,K;Pk表示成像区域内的第k个点目标的空间位置矢量;R(u;Pk)表示成像区域内的第k个点目标的距离历史;
定义参考函数Sref(f,u),记做
S ref ( f , u ) = exp ( jπ T p B f 2 + j 2 π f c + f c Σ k = 1 K R ( u ; P k ) )
其中,f表示快频率;u表示慢时间;
Figure FDA00001908476200021
表示复数单位;B和Tp分别表示发射信号的带宽和脉宽;fc表示雷达工作频率;c表示电磁波传播速度;K表示成像区域内点目标个数;k为自然数,且k=1,2,…,K;Pk表示成像区域内的第k个点目标的空间位置矢量;R(u;Pk)表示成像区域内的第k个点目标的距离历史;
步骤2、快时间傅里叶变换
对原始时域回波s(t,u)的每一列数据,进行快速傅里叶变换,得到快时间傅里叶变换的结果同样为Nr行Na列的回波信号S(f,u);其中,t表示快时间,f表示快频率;u表示慢时间;数据的行数Nr和列数Na的大小由雷达系统的数据采集设备给出;
步骤3、参考函数相乘
用信号S(f,u)与参考函数Sref(f,u)相乘,得到参考函数相乘之后的结果SRFM(f,u),即
SRFM(f,u)=S(f,u)·Sref(f,u)
其中,f表示快频率;u表示慢时间;参考函数Sref(f,u)由步骤1定义;
步骤4、计算最优投影基向量
步骤4.1)根据下式计算矩阵A(m,n),矩阵A(m,n)为Na行Na列的方阵,Na表示数据SRFM(f,u)的列数,Na的值由雷达系统的数据采集设备给出;m为自然数,且m=0,1,…,Na-1;n为自然数,且n=0,1,…,Na-1;
A ( m , n ) = Σ k = 1 K [ R ( u m ; P k ) - R ref ( u m ) ] [ R ( u n ; P k ) - R ref ( u n ) ]
其中,
u m = m + N a / 2 PRF
u n = n + N a / 2 PRF
其中,m为自然数,且m=0,1,…,Na-1;n为自然数,且n=0,1,…,Na-1;Na表示数据SRFM(f,u)的列数,由雷达系统的数据采集设备给出;PRF表示雷达系统的脉冲重复频率;K表示成像区域内点目标个数;k为自然数,且k=1,2,…,K;Pk表示成像区域内的第k个点目标的空间位置矢量;R(u;Pk)表示成像区域内的第k个点目标的距离历史;Rref(u)表示参考函数,由步骤一所定义;
步骤4.2)采用传统常用的特征值分解方法对矩阵A(m,n)进行特征值分解,得到Na个特征值σ1、σ2、…、
Figure FDA00001908476200031
以及与特征值对应的特征向量α1(n)、α2(n)、…、
Figure FDA00001908476200032
Na表示数据SRFM(f,u)的列数,由雷达系统的数据采集设备决定;n为自然数,且n=0,1,…,Na-1;
步骤4.3)对Na个特征值σ1、σ2、…、
Figure FDA00001908476200033
按照由大到小的顺序重新排序,得到排序后的Na个特征值
Figure FDA00001908476200034
Figure FDA00001908476200035
…、以及与特征值对应的特征向量
Figure FDA00001908476200037
Figure FDA00001908476200038
…、其中
Figure FDA000019084762000310
Figure FDA000019084762000311
即为双基地SAR的最优投影基向量;
步骤五、二维非均匀快速傅里叶变换
按照下式对参考函数相乘的结果SRFM(f,u)进行二维非均匀快速傅里叶变换,得到双基地SAR图像I(k,l);
I ( k , l ) = Σ n = - N a / 2 N a / 2 - 1 Σ m = - N r / 2 N r / 2 - 1 S RFM ( mΔf , nΔu ) · exp j 2 π f c + mΔf c α ~ 1 ( n ) k + j 2 π f c + mΔf c α ~ 2 ( n ) l
其中,Δ=B/Nr;Δu=1/PRF;
Figure FDA000019084762000313
Figure FDA000019084762000314
为双基地SAR的最优投影基向量;B表示发射信号的带宽;PRF表示雷达系统的脉冲重复频率;m和n为自然数,且m=-Nr/2,-Nr/2+1,…,Nr/2-1,  n=-Na/2,Na/2+1,…,Na/2-1;fc表示雷达系统的工作频率;c表示电磁波传播速度;k和l为自然数,且m=-Nr/2,-Nr/2+1,…,Nr/2-1,n=-Na/2,Na/2+1,…,Na/2-1;该步骤的具体实现如下:
步骤5.1)根据下式计算得到二维信号Y(g,h)
Y ( g , h ) = Σ g = round ( 2 N r ( f c + mΔf ) α ~ 1 ( n ) c ) + p Σ h = round ( 2 N a ( f c + mΔf ) α ~ 2 ( n ) c ) + q ρ ( p , N r ( f c + mΔf ) α ~ 1 ( n ) c ) × ρ ( q , N a ( f c + mΔf ) α ~ 2 ( n ) c ) × S RFM ( mΔf , nΔu )
其中,函数ρ(q,c)为
ρ ( q , c ) = 1 2 π exp { - 1 4 [ round ( 2 c ) + q - 2 c ] 2 }
其中,round(·)表示四舍五入函数;g为整数,且g=-Nr-Q,-Nr-Q+1,…,Nr+Q;h为整数,且h=-Na-Q,-Na-Q+1,…,Na+Q;行数Nr和列数Na由雷达系统的数据采集设备决定;m和n为自然数,且m=-Nr/2,-Nr/2+1,…,Nr/2-1,n=-Na/2,Na/2+1,…,Na/2-1;fc表示雷达系统的工作频率;c表示电磁波传播速度;p和q为整数,且p=-Q,-Q+1,…,Q,q=-Q,-Q+1,…,Q;为双基地SAR的最优投影基向量;
步骤5.2)利用快速傅里叶变换,根据下式计算得到二维信号
Figure FDA00001908476200045
I ~ ( k , l ) = Σ g = - N r - Q N r + Q Σ h = - N a - Q N a + Q Y ( g , h ) exp ( j π N r gk + j π N a hl )
其中,g为整数,且g=-Nr-Q,Nr-Q+1,…,Nr+Q;h为整数,且h=-Na-Q,-Na-Q+1,…,Na+Q;行数Nr和列数Na由雷达系统的数据采集设备决定;k为整数,且k=-Nr/2,-Nr/2+1,…,Nr/2-1;l为整数,且l=-Na/2,-Na/2+1,…,Na/2-1;
步骤5.3)根据下式,计算得到最终的聚焦图像I(k,l)
I ( k , l ) = exp [ π 2 ( k 2 N r 2 + l 2 N a 2 ) ] · I ~ ( k , l )
其中,k为整数,且k=-Nr/2,-Nr/2+1,…,Nr/2-1;l为整数,且l=-Na/2,-Na/2+1,…,Na/2-1;行数Nr和列数Na的大小由雷达系统的数据采集设备给出;
经过上述步骤处理,就从双基地SAR接收到的目标回波信号s(t,u)中获取具有高分辨率的目标成像结果聚焦图像I(k,l)。
CN201210251052.7A 2012-07-19 2012-07-19 一种基于最优图像空间双基地合成孔径雷达成像方法 Expired - Fee Related CN102798861B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210251052.7A CN102798861B (zh) 2012-07-19 2012-07-19 一种基于最优图像空间双基地合成孔径雷达成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210251052.7A CN102798861B (zh) 2012-07-19 2012-07-19 一种基于最优图像空间双基地合成孔径雷达成像方法

Publications (2)

Publication Number Publication Date
CN102798861A true CN102798861A (zh) 2012-11-28
CN102798861B CN102798861B (zh) 2015-07-15

Family

ID=47198017

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210251052.7A Expired - Fee Related CN102798861B (zh) 2012-07-19 2012-07-19 一种基于最优图像空间双基地合成孔径雷达成像方法

Country Status (1)

Country Link
CN (1) CN102798861B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102967862A (zh) * 2012-12-04 2013-03-13 电子科技大学 双飞移变模式双基地合成孔径雷达成像方法
CN104155653A (zh) * 2014-08-18 2014-11-19 电子科技大学 一种基于特征距离子空间的sar后向投影成像方法
CN103529448B (zh) * 2013-10-18 2015-07-15 电子科技大学 一种基于高精度插值的双基地合成孔径雷达成像方法
CN108267735A (zh) * 2016-12-30 2018-07-10 北京行易道科技有限公司 雷达
CN108267737A (zh) * 2016-12-30 2018-07-10 北京行易道科技有限公司 雷达
CN108267734A (zh) * 2016-12-30 2018-07-10 北京行易道科技有限公司 雷达
CN113933841A (zh) * 2021-11-25 2022-01-14 西安电子科技大学 二维波数均匀映射高分辨成像方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101369017A (zh) * 2007-08-17 2009-02-18 电子科技大学 一种移变模式双基地合成孔径雷达成像方法
CN101369018A (zh) * 2007-08-17 2009-02-18 电子科技大学 一种星机联合双基地合成孔径雷达频域成像方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101369017A (zh) * 2007-08-17 2009-02-18 电子科技大学 一种移变模式双基地合成孔径雷达成像方法
CN101369018A (zh) * 2007-08-17 2009-02-18 电子科技大学 一种星机联合双基地合成孔径雷达频域成像方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YANG YUE 等: "An imaging method and the correction of distortion for Spaceborne-airborne Bistatic SAR", 《2009 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM(IGARSS 2009)》 *
刘喆 等: "星机双基地SAR二维频谱解析表达式求解方法研究", 《电子与信息学报》 *
胥哲铭 等: "基于机载双基地SAR的两种成像算法比较", 《雷达与对抗》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102967862A (zh) * 2012-12-04 2013-03-13 电子科技大学 双飞移变模式双基地合成孔径雷达成像方法
CN102967862B (zh) * 2012-12-04 2014-03-05 电子科技大学 双飞移变模式双基地合成孔径雷达成像方法
CN103529448B (zh) * 2013-10-18 2015-07-15 电子科技大学 一种基于高精度插值的双基地合成孔径雷达成像方法
CN104155653A (zh) * 2014-08-18 2014-11-19 电子科技大学 一种基于特征距离子空间的sar后向投影成像方法
CN104155653B (zh) * 2014-08-18 2017-02-15 电子科技大学 一种基于特征距离子空间的sar后向投影成像方法
CN108267735A (zh) * 2016-12-30 2018-07-10 北京行易道科技有限公司 雷达
CN108267737A (zh) * 2016-12-30 2018-07-10 北京行易道科技有限公司 雷达
CN108267734A (zh) * 2016-12-30 2018-07-10 北京行易道科技有限公司 雷达
CN108267737B (zh) * 2016-12-30 2024-03-22 北京行易道科技有限公司 雷达
CN108267735B (zh) * 2016-12-30 2024-03-26 北京行易道科技有限公司 雷达
CN113933841A (zh) * 2021-11-25 2022-01-14 西安电子科技大学 二维波数均匀映射高分辨成像方法
CN113933841B (zh) * 2021-11-25 2024-06-07 西安电子科技大学 二维波数均匀映射高分辨成像方法

Also Published As

Publication number Publication date
CN102798861B (zh) 2015-07-15

Similar Documents

Publication Publication Date Title
CN102798861B (zh) 一种基于最优图像空间双基地合成孔径雷达成像方法
Wei et al. Linear array SAR imaging via compressed sensing
CN109471080B (zh) 基于simulink的高速平台雷达回波信号模拟系统
CN101369018B (zh) 一种星机联合双基地合成孔径雷达频域成像方法
CN104950306B (zh) 一种海杂波背景下前视海面目标角超分辨成像方法
Feng et al. An extended fast factorized back projection algorithm for missile-borne bistatic forward-looking SAR imaging
CN104833972B (zh) 一种双基地调频连续波合成孔径雷达频率变标成像方法
CN111505639B (zh) 基于变重频采样模式的合成孔径雷达宽幅稀疏成像方法
CN109143237B (zh) 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法
CN108107430A (zh) 一种基于分数阶傅立叶变换的舰船目标isar成像方法
CN102854507B (zh) 一种基于gpu后向投影双站合成孔径雷达成像方法
CN102004250B (zh) 基于频域展开的星机联合双基地合成孔径雷达成像方法
CN104635230A (zh) 一种用于mimo-sar近场测量成像方位向旁瓣抑制的方法
CN102788978B (zh) 一种斜视星机双基地合成孔径雷达成像方法
CN104749573B (zh) 空频域二维稀疏的步进频率sar成像方法
CN103869297B (zh) 高超声速平台载雷达的杂波生成系统及其方法
CN109188436B (zh) 适用于任意平台轨迹的高效双基sar回波生成方法
CN105759264B (zh) 基于时频字典的微动目标缺损回波高分辨成像方法
Sun et al. Polar format algorithm for spotlight bistatic SAR with arbitrary geometry configuration
CN104459694A (zh) 一种基于高阶斜距模型的geo sar高分辨成像方法
CN112649806B (zh) 一种mimo雷达近场三维成像方法
CN113608218A (zh) 一种基于后向投影原理的频域干涉相位稀疏重构方法
CN115657017A (zh) 电大尺寸目标快速成像的增量长度绕射理论射线成像方法
CN101819274A (zh) 用于合成孔径雷达前斜视子孔径成像处理的扩展非线性变标方法
Yang et al. A coordinate-transform based FFBP algorithm for high-resolution spotlight SAR imaging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150715

Termination date: 20200719