CN101509976B - 一种一动一静双基地合成孔径雷达成像方法 - Google Patents
一种一动一静双基地合成孔径雷达成像方法 Download PDFInfo
- Publication number
- CN101509976B CN101509976B CN2009100804721A CN200910080472A CN101509976B CN 101509976 B CN101509976 B CN 101509976B CN 2009100804721 A CN2009100804721 A CN 2009100804721A CN 200910080472 A CN200910080472 A CN 200910080472A CN 101509976 B CN101509976 B CN 101509976B
- Authority
- CN
- China
- Prior art keywords
- target
- frequency domain
- data
- phase place
- dimensional
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种一动一静双基地合成孔径雷达成像方法,属于雷达信号处理技术领域。首先对侦测目标反射回来的雷达信号进行处理,得到雷达回波数据。对雷达回波数据进行二维傅立叶变换,得到二维频域数据,并进行距离频域匹配滤波处理。对滤波处理后得到的二维频域数据的相位进行泰勒展开,之后对其进行相位补偿,最后,对补偿后的相位的空变性进行分析。根据空变性分析结果,对经过相位补偿后的二维频域数据进行波束域的映射和非线性插值;对二维频域数据进行二维逆傅立叶变换,最终得到二维图像。本发明方法精度较高,且可克服大的场景空变特性,适应场景大,运算量小,有利于实时成像。
Description
技术领域
本发明涉及一种双基地合成孔径雷达(BSAR)成像方法,属于雷达信号处理技术领域。
背景技术
双基地合成孔径雷达(BSAR)技术的快速发展隐含着对双基地合成孔径雷达成像技术的巨大需求。目前,各种单基地成像方法在一定的几何配置下都已经移植到双基地合成孔径雷达中,如典型的时域成像(BP)算法,以及典型的频域成像算法距离多普勒(RD)算法,频率变标(CS)算法,极坐标(PFA)算法,波束域(RMA)算法等。目前,世界各国的研究热点集中在RD算法和RMA算法。
RD算法在单基地SAR中是最基本的算法,认为距离维度和方位维度解耦,因此可以直接应用二维傅立叶变换进行处理,简单高效。但是复杂的双基地配置距离维度和方位维度耦合紧密,大部分情况下RD算法无法使用,因此目前BSAR中应用RD算法的几何配置只是应用在平行等速飞行情况下。
RMA算法在单基地SAR中是完全精准的算法,在高精度的成像场合经常被选用。在BSAR中,各国研究者也都在试图借用单基地RMA算法的思路去推导双基地SAR波束域算法,到目前为止,完全精准的双基地波束域算法并没有出现,研究的比较成熟的仍然是平行等速飞行情况,并且相应的实测数据成像结果也已给出。针对其它轨迹配置的双基地SAR情况,都有各种各样的近似的波束域算法。
BP算法是最精准的成像算法,无论是在单基地情况还是双基地情况下都能适用,但是该算法采用的是点对点的精确配准方法,因此成像算法的效率较为低下,其运算时间是普通频域成像算法的N倍,N为方位向采样点数。导致其在实际系统中并无广泛的应用。
当发射机与接收机不在同一平台,并且发射机的位置始终在变化,而接收机的位置始终不动,即构成了一动一静BSAR的系统结构。对于一动一静BSAR系统而言,由于几何结构的不对称性导致严重的空变特性,通常为了获得高精度的成像可以采用时域成像算法,但是运算量非常大。
一动一静BSAR系统成像的几何配置如图1所示,以场景中心为坐标原点建立坐标系,场景平面位于xy坐标系内。发射机(如卫星平台)飞行的方向为x轴正方向,飞行轨迹平行于x轴;场景平面的法线方向为z轴方向;y轴通过右手螺旋法则确定。在整个卫星聚束照射的过程中,卫星的飞行轨迹关于yz面是对称的。合成孔径的中心坐标为(0,yt0,zt0),接收机的位置始终不动,其坐标为(xr,yr,zr),场景中心点坐标为(0,0,0),发射机的飞行速度为恒定值vt。假设发射机的发射的基带信号为p(t),在场景坐标为(X,Y,0)处有一个点目标,其后向散射系数为δ(X,Y)。设合成孔径(SA)时间长度为Ts,在某个方位向的采样时刻为tn,这时的采样点处的坐标为(vttn,yt0,zt0),其对应的基带回波信号表达式为:
在上式中,包含着点目标的相关信息:第一项为点目标的后向散射系数,假设其为空域独立的;第二项为发射信号的时延信号,带有点目标的斜距信息;第三项是由平台运动所引入的多普勒相位项,带有点目标的方位坐标信息。参数定义为:λ为载波波长,f0为载波频率,有f0=c/λ,c为光速,t为快时间变量,tn为慢时间变量,τ为目标回波时延,p(t)为发射的基带信号,R(tn;X,Y)为第n个脉冲发射时刻点目标的斜距,即从发射机到点目标再到接收机的距离长度,表达式为:
τ(tn;X,Y)为tn时的时延,如下所示:
针对这种几何配置,需要一种适用于一动一静BSAR下的高精度成像方法,来解决成像过程中大场景空变特性及运算量大导致的计算效率低下的问题。
发明内容
本发明的目的是为了解决一动一静BSAR系统在成像过程中大场景的空变特性问题以及计算效率低下问题,提出一种一动一静双基地合成孔径雷达成像方法。
为达到上述目的,本发明提供的技术方案如下:
一种一动一静双基地合成孔径雷达成像方法,其步骤如下:
动态SAR发射机向侦测目标发出雷达信号,静止接收机接收侦测目标反射回来的雷达信号,在静止接收机中,进行如下处理:
(1)对侦测目标(简称目标)反射回来的雷达信号进行处理,得到雷达回波数据。得到的雷达回波数据S(t,tn)表示为:
式(4)中,t为距离快时间变量;tn为方位慢时间变量;λ为发射信号波长;(X,Y)为目标位置坐标;Ω为目标场景区域;δ(X,Y)为目标的后向散射系数;p(t)为发射的基带信号,τ(tn;X,Y)为发射机到目标再到接收机的传播时间,因此,p(t-τ(tn;X,Y))即为存在时间延迟的基带信号;R(tn;X,Y)为第n个脉冲发射时刻点目标的斜距,即从发射机到点目标再到接收机的距离长度之和;j为虚数单位。其中,R(tn;X,Y)表达式为:
式(5)中,vt为发射机运动速度,发射机在合成孔径中心的坐标为(0,yt0,zt0),为运动的发射机相对于目标的斜距。接收机的位置始终不动,其坐标为(xr,yr,zr)。为接收机相对于目标的斜距。式(5)表示一动一静双基地几何配置下的斜距和。因此,式(4)也即表达一动一静BSAR配置下得到的雷达回波数据。
(2)对步骤(1)生成的雷达回波数据S(t,tn)进行二维频域转换,得到目标的距离和方位二维频域数据。对所得二维频域数据进行距离频域匹配滤波处理,得到抵消距离向频域相位后的二维频域数据S(kr,kx)。
为了解决传统时域算法运算量大的问题,可以把雷达回波数据S(t,tn)变换到二维频域进行操作处理,以便利用快速傅立叶变换算法。同时,为了对目标进行距离向的聚焦,也需要将数据变换到频域进行匹配滤波。并且在一动一静双基地配置下找到雷达回波数据S(t,tn)的精确频域表达式,便于寻找快速傅立叶变换与目标回波频域数据之间的映射关系。所以对步骤(1)生成的雷达回波数据S(t,tn)进行距离向傅立叶变换和方位向傅立叶变换,即二维频域转换,得到距离和方位二维频域数据,并对其进行距离频域匹配滤波处理,得到抵消距离向频域相位后的目标回波二维频域数据的相位Φ(kr,kx),表达式为
式(6)中,kr为距离向波数,写为2π(f+f0)/c,f为距离向频率,f0为发射载波频率,c为光速;kx为方位向波数,写为2πfd/vt,fd为方位多普勒频率;因此,二维频域数据S(kr,kx)表示为:
S(kr,kx)=∫∫Ωδ(X,Y)exp(jΦ(kr,kx))dXdY (7)
式(7)中,kr为距离向波数;kx为方位向波数;j为虚数单位;Ω为目标场景区域;δ(X,Y)为目标的后向散射系数;(X,Y)为目标位置坐标;
式(6)明确的显示目标回波频域相位Φ(kr,kx)不是目标位置坐标(X,Y)的线性函数,因此无法直接对目标回波应用快速傅立叶变换算法。因此,需要对目标回波频域相位Φ(kr,kx)进行泰勒展开,得到目标回波频域相位Φ(kr,kx)与目标位置坐标(X,Y)的线性关系。
(3)对经步骤(2)滤波处理后得到的目标回波频域相位Φ(kr,kx)进行泰勒展开,得到目标回波频域相位Φ(kr,kx)与目标位置坐标(X,Y)间的线性关系。之后,对目标回波频域相位Φ(kr,kx)进行相位补偿。对补偿后的相位的空变性进行分析,确定出对目标场景成像的范围。
将式(6)进行泰勒展开保留到目标位置的一次项,但是这样会引入一定的误差,并且这个误差不是恒定的,而是空变的,因此需要对目标的场景空变特性进行分析。对经步骤(2)滤波处理后得到目标回波频域相位Φ(kr,kx)进行泰勒展开,得到一常数相位项和一次相位项,并且常数相位项就是场景中心的方位向相位,由于场景中心坐标已知,所以能够进行补偿,补偿之后的回波频域相位即为与目标位置坐标(X,Y)对应的线性相位。经过相位补偿之后,其残留的二次相位和更高阶的相位即为误差相位,并且具有空变特性,影响最终的目标聚焦效果。由于残留的二次相位误差远比高阶相位误差大,因此主要考虑二次相位误差的影响。
为保证目标的聚焦程度以π/4作为相位误差的界限,得到对目标空变特性能够容忍的场景尺寸如下:
其中
式(8)、(9)中,ρX为地面方位向分辨率大小,(0,yt0,zt0)为孔径中心时刻发射机的位置坐标,(xr,yr,zr)为接收机位置坐标。即,在式(8)给定的场景尺寸范围内,近似误差很小,在保证回波数据频域相位精准的同时也找到了目标位置相关的一次相位项。
(4)对步骤(3)中经过相位补偿后的二维频域数据S(kr,kx)进行波束域映射和非线性插值,使得目标位置坐标(X,Y)的回波频域数据Δ(kr,kx)的相位与波束域参数(kX,kY)之间构成线性关系,以便能够使用快速傅立叶变换算法进行成像。
在对目标回波频域相位Φ(kr,kx)进行相位补偿之后,目标的二维频域表达式是距离波束参数kr(与距离频率等价)与方位波束参数kx(与方位频率等价)的关系式,其中kr与kx是相互耦合的。虽然目标回波频域相位Φ(kr,kx)与目标位置坐标(X,Y)间是线性关系,但是Φ(kr,kx)与kr和kx却不是线性关系。为了能够采用高效的快速傅立叶变换算法进行成像,需要对kr和kx进行波束域映射,生成新的波束域参数(kX,kY),通过非线性插值运算实现此步操作,使得目标回波频域相位Φ(kr,kx)与(kX,kY)构成线性关系。经过非线性插值后的回波频域数据Δ(kr,kx)表达式为:
Δ(kX,kY)=∫∫Ωδ(X,Y)exp(-jkXX-jkYY)dXdY
(10)
从式(10)可以看出,经过非线性插值后的二维频域数据Δ(kr,kx)的相位(-jkXX-jkYY)与目标位置坐标(X,Y)构成傅立叶变换对关系,因此可以直接使用高效的快速傅立叶变换算法,达到减小运算量的目的。
(5)对经步骤(4)处理后的二维频域数据Δ(kr,kx)进行距离向快速逆傅立叶变换,得到目标距离向的图像,然后进行方位向快速逆傅立叶变换,得到目标方位向的图像。距离向图像和目标方位向图像相结合,即为目标的二维图像。最后,结合步骤(3)得到的目标场景成像的范围,选取出最终的目标图像。至此,就实现了一动一静双基地合成孔径雷达的成像。
有益效果
本发明提出的一动一静双基地合成孔径雷达成像方法,通过将目标雷达数据变换到频域进行匹配滤波,对目标回波频域相位进行相位补偿。对补偿后的相位的空变性进行分析,确定出对目标场景成像的范围。对相位补偿后的二维频域数据进行波束域映射和非线性插值,使用快速傅立叶变换算法进行成像。成像精度高,且可克服大场景空变特性,适应场景大,运算量小,有利于实时成像。
附图说明
图1为一动一静BSAR成像的几何配置示意图;
图2为本发明优选实施方式中的一动一静BSAR成像方法的流程图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,现结合附图对本发明方法作进一步详细说明。
一种一动一静双基地合成孔径雷达成像方法,包括以下步骤:
动态SAR发射机向侦测目标发出雷达信号,静止接收机接收侦测目标反射回来的雷达信号,在静止接收机中,进行如下处理:
(1)对侦测目标(简称目标)反射回来的雷达信号进行处理,得到雷达回波数据。
对整个场景Ω上的点目标回波进行叠加处理,即对单点回波信号(1)式在整个场景Ω上积分,可得到整个场景的雷达回波数据S(t,tn),表示为:
式(11)中,t为距离快时间变量;tn为方位慢时间变量;λ为发射信号波长;(X,Y)为目标位置坐标;Ω为目标场景区域;δ(X,Y)为目标的后向散射系数;p(t)为发射的基带信号,τ(tn;X,Y)为发射机到目标再到接收机的传播时间,因此,p(t-τ(tn;X,Y))即为存在时间延迟的基带信号;R(tn;X,Y)为第n个脉冲发射时刻点目标的斜距,即从发射机到点目标再到接收机的距离长度之和;j为虚数单位。其中,R(tn;X,Y)表达式为:
式(12)中,vt为发射机运动速度。发射机在合成孔径中心的坐标为(0,yt0,zt0)。式(12)中,为运动的发射机相对于目标的斜距。接收机的位置始终不动,其坐标为(xr,yr,zr);为接收机相对于目标的斜距。式(12)表示一动一静双基地几何配置下的斜距和。因此,式(11)也即表达一动一静BSAR配置下得到的雷达回波数据。
之后,根据s(t,tn)进行一动一静双基地SAR回波数据的生成。回波数据生成可以使用Matlab或VC等工具软件来完成。
(2)对步骤(1)生成的雷达回波数据S(t,tn)进行二维频域转换,得到目标的距离和方位二维频域数据。对所得二维频域数据进行距离频域匹配滤波处理,得到抵消距离向频域相位后的二维频域数据S(kr,kx)。
具体如下:
首先,对步骤(1)得到的雷达回波数据S(t,tn)采用二维傅立叶变换方式进行二维频域转换,得到雷达回波数据的二维频域表达式。通过与发射基带信号的频域表达式共轭相乘,完成距离向的频域匹配滤波。
在BSAR中,高效的二维频域算法难以实现的原因有两点,一是在当雷达回波数据变换到二维频域之时,目标的驻定相位点难以获得,也就是说精确的二维频域表达式无法得到;二是二维波束域与目标的空间坐标不是线性的映射关系,无法直接应用高效的快速傅立叶运算。一动一静BSAR能够求解到目标二维频域的精确表达式,但是其方位频率距离频率与目标方位坐标和距离坐标关系非常复杂。经过驻定相位原理和距离向频域匹配滤波之后,得到抵消距离向频域相位后的目标回波二维频域数据的相位Φ(kr,kx),表达式为
式(13)中,kr为距离向波数,写为2π(f+f0)/c,f为距离向频率,f0为发射载波频率,c为光速;kx为方位向波数,写为2πfd/vt,fd为方位多普勒频率;因此,二维频域数据S(kr,kx)表示为:
S(kr,kx)=∫∫Ωδ(X,Y)exp(jΦ(kr,kx))dXdY (14)
式(14)中,kr为距离向波数;kx为方位向波数;j为虚数单位;Ω为目标场景区域;δ(X,Y)为目标的后向散射系数;(X,Y)为目标位置坐标;
(3)对经步骤(2)滤波处理后得到的目标回波频域相位Φ(kr,kx)进行泰勒展开,得到目标回波频域相位Φ(kr,kx)与目标位置坐标(X,Y)间的线性关系。之后,对目标回波频域相位Φ(kr,kx)进行相位补偿。对补偿后的相位的空变性进行分析,确定出对目标场景成像的范围。
式(13)明确的显示出目标回波频域相位Φ(kr,kx)不是目标位置坐标(X,Y)的线性函数。要想得到它们之间的线性关系,可以采用泰勒展开,将式(13)展开保留到目标位置坐标(X,Y)的一次项即可,但是这样又会引入一定的误差,并且这个误差不是恒定的,而是空变的,因此还需要对目标的场景空变特性进行分析。
步骤(3)的具体过程如下:
对式(13)进行泰勒(Taylor)展开,得到一常数项相位和一次相位,并且常数项相位就是场景中心的方位向相位。由于场景中心坐标已知,所以能够对常数项相位进行补偿,补偿之后的相位即为目标位置的线性相位。经过相位补偿之后,其残留的二次相位和更高阶的相位为误差相位,并且具有空变特性,对其进行分析以确定此近似对场景空变特性的适应性。由于残留的二次相位误差远比高阶相位误差大,因此主要考虑二次相位误差的影响。首先,泰勒展开后目标回波频域相位Φ(kr,kx)的表达式为:
由式(15)可知,前两项相位为场景中心的相位,也即需要补偿的二维频域相位,通过在二维频域乘上相位因子 就能够完成相位的补偿。后两项相位为场景目标的线性相位。对比式(15)和式(13),可知二次相位误差ΔΦ为:
其中,kr为距离波束域,kx为方位波束域。从式(16)可以看出,目标的二次相位误差ΔΦ与多个因素有关,如目标偏离场景中心的相对位置,发射机航迹到场景中心的最短距离,接收机航迹到场景中心的最短距离等。式(16)的二次相位误差ΔΦ是kr与kx的关系式,其中kr与kx是相互耦合的,为了分析二次相位误差ΔΦ空变特性影响的大小,需要对kr与kx域解耦。经过解耦合,可知二次相位误差ΔΦ主要产生目标位置偏移和目标散焦,目标的位置偏移为常数量,可以进行补偿,因此二次相位误差ΔΦ对场景空变特性的影响主要体现在对场景目标的散焦上,当目标位置偏离场景中心越远,空变特性越大,其散焦的越厉害。为保证目标的聚焦程度以π/4作为相位误差的界限,得到对目标空变特性能够容忍的场景尺寸如下:
其中
式(17)中,ρX为地面方位向分辨率大小,(0,yt0,zt0)为孔径中心时刻发射机的位置坐标,(xr,yr,zr)为接收机位置坐标。
(4)对步骤(3)中经过相位补偿后的二维频域数据S(kr,kx)进行波束域映射和非线性插值,使得目标位置坐标(X,Y)的二维频域数据Δ(kr,kx)的相位与波束域参数(kX,kY)之间构成线性关系,以便能够使用快速傅立叶变换算法进行成像。
具体过程如下:
在对目标回波频域相位Φ(kr,kx)进行相位补偿之后,目标的二维频域表达式是距离波束参数kr(与距离频率等价)与方位波束参数kx(与方位频率等价)的关系式,其中kr与kx是相互耦合的。虽然目标回波频域相位Φ(kr,kx)与目标位置坐标(X,Y)间是线性关系,但是Φ(kr,kx)与kr和kx却不是线性关系。因此需要进行波数域的映射,找到与目标位置坐标(X,Y)一一映射的波数域集合(kX,kY),并且(kX,kY)可以通过(kr,kx)集合经非线性插值得到。同理,当(kX,kY)与(kr,kx)满足如下关系时:
所得到的(kX,kY)与(X,Y)构成傅立叶变换对的关系,并且kX与kY相互正交,也就是说两者解耦合。同样,从式(19)可以知道由已有的样本波数域数据(kr,kx)得到解耦合的波数域数据(kX,kY),需要经过复杂的非线性插值。非线性插值的实现方法如下:首先,根据系统参数确定波数域kr和kx的取值方位。然后,根据公式(19)和波数域kr和kx的取值方位,确定出kX和kY的波数域取值范围。之后,选取合适的kX和kY范围,以保证频率支撑域的有效性,在有效的频率支撑域里面,均匀划分kX和kY,经由公式(19)的逆运算,反解出对应的kr和kx。最后,利用反解出的kr和kx,采用sinc插值的方法,从已有的频率数据域(kr,kx)插值得到新的频率数据域(kX,kY)。非线性插值后的二维频域数据Δ(kr,kx)表达式为:
Δ(kX,kY)=∫∫Ωδ(X,Y)exp(-jkXX-jkYY)dXdY
(20)
从式(20)可以看出,经过非线性插值后的二维频域数据Δ(kr,kx)的相位(-jkXX-jkYY)与目标位置坐标(X,Y)构成傅立叶变换对关系,因此可以直接使用高效的快速傅立叶变换算法,达到减小运算量的目的。
(5)对经步骤(4)处理后的二维频域数据Δ(kr,kx)进行距离向快速逆傅立叶变换,得到目标距离向的图像,然后进行方位向快速逆傅立叶变换,得到目标方位向的图像。距离向图像和目标方位向图像相结合,即为目标的二维图像。最后,结合步骤(3)得到的目标场景成像的范围,选取出最终的目标图像。至此,就实现了一动一静双基地合成孔径雷达的成像。
从公式(19)和公式(20)可知,经过相位补偿、波数域的映射和非线性插值后,已经得到了正交的二维频域数据Δ(kr,kx);同样从公式(20)可知,非线性插值后得到的二维频域数据Δ(kr,kx)也是目标后向散射系数的傅立叶变换得到的数据,因此对插值后得到的二维频域数据Δ(kr,kx)进行二维逆傅立叶变换(2D-IFFT)就可以得到目标的后向散射系数,完成由数据空间到图像空间的重构。其数学表达式描述为:
Claims (1)
1.一种一动一静双基地合成孔径雷达成像方法,实现过程如下:
动态SAR发射机向侦测目标发出雷达信号,静止接收机接收侦测目标反射回来的雷达信号,其特征在于,在静止接收机中,
(1)对侦测目标反射回来的雷达信号进行处理,得到雷达回波数据,得到的雷达回波数据S(t,tn)表示为:
式(1)中,t为距离快时间变量;tn为方位慢时间变量;λ为发射信号波长;(X,Y)为目标位置坐标;Ω为目标场景区域;δ(X,Y)为目标的后向散射系数;p(t)为发射的基带信号,τ(tn;X,Y)为发射机到目标再到接收机的传播时间,p(t-τ(tn;X,Y))即为存在时间延迟的基带信号;R(tn;X,Y)为第n个脉冲发射时刻点目标的斜距,即从发射机到点目标再到接收机的距离长度之和;j为虚数单位;其中,R(tn;X,Y)表达式为:
(2)对步骤(1)生成的雷达回波数据S(t,tn)进行二维频域转换,得到目标的距离和方位二维频域数据,对所得二维频域数据进行距离频域匹配滤波处理,得到抵消距离向频域相位后的目标回波二维频域数据的相位Φ(kr,kx),表达式为:
式(3)中,kr为距离向波数,写为2π(f+f0)/c,f为距离向频率,f0为发射载波频率,c为光速;kx为方位向波数,写为2πfd/vt,fd为方位多普勒频率;因此,二维频域数据S(kr,kx)表示为:
S(kr,kx)=∫∫Ωδ(X,Y)exp(jΦ(kr,kx))dXdY (4)
式(4)中,kr为距离向波数;kx为方位向波数;j为虚数单位;Ω为目标场景区域;δ(X,Y)为目标的后向散射系数;(X,Y)为目标位置坐标;
(3)对经步骤(2)滤波处理后得到的目标回波频域相位Φ(kr,kx)进行泰勒展开,得到目标回波频域相位Φ(kr,kx)与目标位置坐标(X,Y)间的线性关系,之后,对目标回波频域相位Φ(kr,kx)进行相位补偿;对补偿后的相位的空变性进行分析,确定出对目标场景成像的范围,即:
为保证目标的聚焦程度以π/4作为相位误差的界限,得到对目标空变特性能够容忍的场景尺寸如下:
其中
式(5)、(6)中,ρX为地面方位向分辨率大小,(0,yt0,zt0)为孔径中心时刻发射机的位置坐标,(xr,yr,zr)为接收机位置坐标;
(4)对步骤(3)中经过相位补偿后的二维频域数据S(kr,kx)进行波束域映射和非线性插值,使得目标位置坐标(X,Y)的回波频域数据Δ(kr,kx)的相位与波束域参数(kX,kY)之间构成线性关系,以便能够使用快速傅立叶变换算法进行成像;
在对目标回波频域相位Φ(kr,kx)进行相位补偿之后,目标的二维频域表达式是距离波束参数kr与方位波束参数kx的关系式,其中kr与kx是相互耦合的;虽然目标回波频域相位Φ(kr,kx)与目标位置坐标(X,Y)间是线性关系,但是Φ(kr,kx)与kr和kx却不是线性关系;因此需要进行波数域的映射,找到与目标位置坐标(X,Y)一一映射的波数域集合(kX,kY),并且(kX,kY)可以通过(kr,kx)集合经非线性插值得到;同理,当(kX,kY)与(kr,kx)满足如下关系时:
所得到的(kX,kY)与(X,Y)构成傅立叶变换对的关系,并且kX与kY相互正交,也就是说两者解耦合;同样,从式(19)可以知道由已有的样本波数域数据(kr,kx)得到解耦合的波数域数据(kX,kY),需要经过复杂的非线性插值;
经过非线性插值后的回波频域数据Δ(kr,kx)表达式为:
Δ(kX,kY)=∫∫Ωδ(X,Y)exp(-jkXX-jkYY)dXdY
(7)
(5)对经步骤(4)处理后的二维频域数据Δ(kr,kx)进行距离向快速逆傅立叶变换,得到目标距离向的图像,然后进行方位向快速逆傅立叶变换,得到目标方位向的图像;距离向图像和目标方位向图像相结合,即为目标的二维图像;最后,结合步骤(3)得到的目标场景成像的范围,选取出最终的目标图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100804721A CN101509976B (zh) | 2009-03-19 | 2009-03-19 | 一种一动一静双基地合成孔径雷达成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100804721A CN101509976B (zh) | 2009-03-19 | 2009-03-19 | 一种一动一静双基地合成孔径雷达成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101509976A CN101509976A (zh) | 2009-08-19 |
CN101509976B true CN101509976B (zh) | 2011-01-26 |
Family
ID=41002434
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100804721A Expired - Fee Related CN101509976B (zh) | 2009-03-19 | 2009-03-19 | 一种一动一静双基地合成孔径雷达成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101509976B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102707283B (zh) * | 2012-06-15 | 2013-08-14 | 电子科技大学 | 一种固定站双基地合成孔径雷达成像方法 |
CN103323841B (zh) * | 2012-09-03 | 2015-08-19 | 中国科学院电子学研究所 | 星机载移变式双基合成孔径雷达中的成像方法及装置 |
CN102914775B (zh) * | 2012-10-10 | 2014-04-02 | 北京理工大学 | 一种改进的双基地合成孔径雷达nlcs成像算法 |
CN102937709B (zh) * | 2012-10-16 | 2014-05-07 | 北京理工大学 | 基于直达波的双基地sar无同步回波聚焦成像方法 |
CN103837874B (zh) * | 2014-03-05 | 2016-01-20 | 北京理工大学 | 用于地球同步轨道sar成像的二维非线性变调频方法 |
CN105182340B (zh) * | 2015-09-28 | 2017-09-26 | 电子科技大学 | 一种双基地前视sar运动补偿方法 |
CN106842200A (zh) | 2017-01-11 | 2017-06-13 | 中国科学院电子学研究所 | 一种双基合成孔径雷达成像方法和装置 |
CN108152818A (zh) * | 2017-12-05 | 2018-06-12 | 上海无线电设备研究所 | 基于改进结构稀疏模型的宽角sar成像算法 |
CN109143237B (zh) * | 2018-08-31 | 2020-10-09 | 电子科技大学 | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 |
CN109597072B (zh) * | 2018-12-18 | 2019-12-10 | 中国科学院电子学研究所 | 一种双基合成孔径雷达sar系统的成像处理方法及装置 |
CN109875606B (zh) * | 2019-01-22 | 2020-05-22 | 北京航空航天大学 | 基于先验反射成像的超声ct声速成像的方法 |
CN110568410B (zh) * | 2019-10-09 | 2021-08-31 | 上海无线电设备研究所 | 一种空间频率色散的微波雷达超分辨方法 |
CN110837128B (zh) * | 2019-11-26 | 2021-09-10 | 内蒙古工业大学 | 一种柱面阵列雷达的成像方法 |
CN112698279A (zh) * | 2020-11-19 | 2021-04-23 | 中山大学 | 一站固定双基地sar回波信号快速生成方法 |
CN113608215A (zh) * | 2021-06-21 | 2021-11-05 | 北方工业大学 | 基于三角正弦等效驻留相位点求解的波数域ArcSAR成像方法 |
-
2009
- 2009-03-19 CN CN2009100804721A patent/CN101509976B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101509976A (zh) | 2009-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101509976B (zh) | 一种一动一静双基地合成孔径雷达成像方法 | |
CN101430380B (zh) | 基于非均匀采样的大斜视角机载sar聚束模式成像方法 | |
CN101369018B (zh) | 一种星机联合双基地合成孔径雷达频域成像方法 | |
CN104833972B (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
Mao et al. | Autofocus correction of APE and residual RCM in spotlight SAR polar format imagery | |
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN102707269B (zh) | 一种机载雷达距离走动校正方法 | |
CN104898119A (zh) | 一种基于相关函数的动目标参数估计方法 | |
CN102288961A (zh) | 一种合成孔径雷达非线性调频变标的成像方法 | |
CN102121990B (zh) | 基于空时分析的逆合成孔径雷达的目标转速估计方法 | |
CN102759734B (zh) | 一种用于高分辨率星载sar的改进距离徙动校正算法的成像方法 | |
CN105842696B (zh) | 一种基于可旋转前视阵的斜视InSAR地面动目标检测方法 | |
Huang et al. | Detection and fast motion parameter estimation for target with range walk effect based on new axis rotation moving target detection | |
CN103792534B (zh) | 一种基于先验相位结构知识的sar两维自聚焦方法 | |
CN112415512B (zh) | 基于进退法和黄金分割法的sar运动目标聚焦方法 | |
CN102890277B (zh) | 移不变双基地合成孔径雷达距离徙动成像方法 | |
Xu et al. | A variable PRF imaging method for high squint diving SAR | |
Wang et al. | High squint mode SAR imaging using modified RD algorithm | |
Wang et al. | ISAR imaging of target with complex motion based on novel approach for the parameters estimation of multi-component cubic phase signal | |
CN102162845B (zh) | 双基地合成孔径雷达点目标响应二维频谱的计算方法 | |
Wang et al. | Study on the echo signal model and RD imaging algorithm for FMCW SAR | |
Wang et al. | Study on motion compensation method for W-band UAV MISAR real-time imaging | |
CN112859018A (zh) | 一种基于图像几何校正的视频sar成像方法 | |
Mao et al. | Doppler centroid estimation for Doppler beam sharpening imaging based on the morphological edge detection method |
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 |
Granted publication date: 20110126 Termination date: 20150319 |
|
EXPY | Termination of patent right or utility model |