CN109143236A - 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 - Google Patents
适用于复杂飞行轨迹的双基聚束sar大场景成像方法 Download PDFInfo
- Publication number
- CN109143236A CN109143236A CN201811007730.9A CN201811007730A CN109143236A CN 109143236 A CN109143236 A CN 109143236A CN 201811007730 A CN201811007730 A CN 201811007730A CN 109143236 A CN109143236 A CN 109143236A
- Authority
- CN
- China
- Prior art keywords
- distance
- imaging
- dimensional
- phase compensation
- bistatic
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 82
- 238000001228 spectrum Methods 0.000 claims abstract description 19
- 238000007906 compression Methods 0.000 claims abstract description 17
- 230000006835 compression Effects 0.000 claims abstract description 13
- 230000009466 transformation Effects 0.000 claims abstract description 6
- 238000013507 mapping Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 11
- 230000001427 coherent effect Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 7
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 20
- 238000012545 processing Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 8
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000007547 defect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000013506 data mapping Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 230000009131 signaling function Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
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/904—SAR modes
- G01S13/9058—Bistatic or multistatic SAR
-
- 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
-
- 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
- G01S13/9052—Spotlight mode
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大场景成像方法,其包括成像系统参数初始化,对回波数据进行基于点目标参考函数的一致压缩,对一致压缩后的距离频域数据做波数谱一致化变换,高阶空变相位补偿成像。本发明采用基于线性空变假设的波数谱一致化方法实现复杂飞行轨迹下的SAR回波一致聚焦,并通过高阶空变相位补偿提高有效聚焦景深,实现了复杂飞行轨迹下的双基SAR大场景成像处理,可以在较大的成像场景范围内,实现高分辨率的成像处理,使搭载双基SAR系统的飞行器在战场侦察监视、物资空投和地震灾害救援等领域发挥更好的性能。
Description
技术领域
本发明属于雷达信号处理技术领域,具体涉及一种适用于复杂飞行轨迹的双基聚束SAR大场景成像方法。
背景技术
作为一种主动的航空、航天遥感手段,微波成像技术具有全天时、全天候工作的特点,在地质测绘、灾害监测、军事侦察等领域有着广泛的应用,目前已成为高分辨率对地观测和全球资源管理最重要手段之一。但是由于本身工作体制的限制,现有的合成孔径雷达(SAR,Synthetic Aperture Radar)并不能实现前视区域方位向高分辨成像,从而使SAR在飞行器前视对地、自主着陆、物资空投、导弹末制导等方面不能充分的发挥作用。
双基地合成孔径雷达(BSAR,bistatic SAR)是一种SAR体制推广到双基地雷达的微波成像技术,不仅具有几何构型灵活、隐蔽性强等优势,还能够克服单基地SAR存在前视盲区的问题。在雷达平台运动过程中,发射站天线对成像区域进行照射、接收站天线接收成像区域中的目标散射回波;利用发射信号的大带宽形成距离向高分辨,通过成像处理算法补偿方位向信号的多普勒相位以实现方位向孔径合成进而形成方位向高分辨,从而实现成像区域内的两维高分辨成像。
目前的双基SAR成像方法主要存在两方面缺陷。一方面,SAR成像处理往往采用参考点目标的频谱函数(Point Target Reference Spectrum,PTRS)对回波信号进行一致压缩。而随着成像场景增大,场景中远处的目标回波的特性参数相对于参考点的差异增大,导致其方位向信号的相干叠加效果恶化,进而限制了双基地SAR的有效成像场景尺寸。在大场景范围内实现高分辨雷达成像,可以让飞行器在战场侦察监视、物资空投和地震灾害救援等方面发挥更好的性能。文献“Wu,J.,Li,Z.,Huang,Y.,Yang,J.,&Liu,Q.H.A generalizedomega-k algorithm to process translationally variant bistatic-sar data basedon two-dimensional stolt mapping.IEEE Transactions on Geoscience&RemoteSensing,52(10),6597-6614,2014”将回波信号频谱参数的空间变化规律近似建模为线性,一定程度上提高了一致聚焦的有效场景范围,但当场景尺寸进一步增大,线性近似仍然无法满足高精度成像的要求。另一方面,关于SAR成像方法的理论研究往往将雷达平台的飞行建模为匀速直线运动,而实际工作中平台飞行往往存在一定的机动,运动轨迹的高阶分量将带来显著的相位误差,进一步限制有效的成像场景尺寸。目前关于双基聚束SAR成像方法的研究结果中,还没有一种方法能够在保证算法效率的前提下,同时解决复杂飞行轨迹和大场景一致聚焦两大问题。
发明内容
本发明的发明目的是:为了解决现有技术在双基地聚束SAR成像处理中的缺陷,本发明提出了一种适用于复杂飞行轨迹的双基聚束SAR大场景成像方法。
本发明的技术方案是:一种适用于复杂飞行轨迹的双基聚束SAR大场景成像方法,包括以下步骤:
A、成像系统参数初始化
所述成像系统参数包括参考点位置、发射站三维轨迹[xT(t),yT(t),zT(t)]、接收站三维轨迹[xR(t),yR(t),zR(t)]、发射信号载波频率fc、脉冲宽度Tp、距离向采样频率Fs、距离向采样点数Nr、脉冲重复间隔PRI、方位时间向量t、距离时间向量τ、距离频率向量fτ、场景中点目标(xp,yp)的距离历史Rb(t;xp,yp);
B、对回波数据进行基于点目标参考函数的一致压缩
根据参考点目标的回波数据构造点目标参考函数,利用距离频域和方位时域对回波数据进行一致压缩,得到一致压缩后的距离频域数据;
C、对一致压缩后的距离频域数据做波数谱一致化变换
将差分距离史对空间坐标二维泰勒展开,根据双平台轨迹信息构造映射关系,通过二维插值将一致压缩后的数据变换到二维波数域,得到波数谱一致化的数据;
D、高阶空变相位补偿
将粗成像结果划分为多个子图像,构造空变相位补偿因子,计算每个子图像对应的相位补偿因子并进行相位补偿,将所有相位补偿后的子图像进行拼接得到成像结果。
进一步地,所述步骤B中,点目标参考函数表示为
其中,So(fτ,t)表示点目标参考函数,Kr为发射信号调频率,c为光速。
进一步地,所述步骤B中,利用距离频域、方位时域的共轭相乘,对回波数据进行一致压缩,一致压缩后的距离频域数据表示为
其中,FFTr{·}表示距离向FFT,sr(τ,t;xp,yp)表示目标回波数据,So *(fτ,t)表示对一致压缩参考信号So(fτ,t)取共轭,ΔRp(t;xp,yp)表示差分距离史,wr[·]表示距离向时域窗函数,wa[·]表示方位向时域窗函数,Ta表示合成孔径时间。
进一步地,所述步骤C中,将差分距离史对空间坐标二维泰勒展开,表示为
其中,Rcur(t,xp,yp)表示距离误差。
进一步地,所述步骤D中,首先计算补偿前的有效场景尺寸,计算公式为
其中,Qi和Ti均为雷达平台运动轨迹相关系数,λ表示载波波长;
再根据有效场景尺寸将粗成像结果hcur(x,y)划分为多个子图像hi(x,y)。
进一步地,所述步骤D中,构造空变相位补偿因子,表示为
Hc(kx,ky;xp,yp)=exp{jkr(kx,ky)Rcur[t(kx,ky);xp,yp]}
其中,kx和ky分别为x和y方向上的波数,kr(kx,ky)表示由kx和ky来表示kr的二维函数映射,t(kx,ky)表示由kx和ky来表示t的二维函数映射。
进一步地,所述步骤D中,对每个子图像hi(x,y),先确定其几何中心位置(xi,yi),将子图像进行二维快速傅里叶变换,得到其波数谱Hi(kx,ky);再计算出对应的相位补偿因子Hc(kx,ky;xi,yi),与子图像波数谱相乘后再求解乘积结果的二维逆向快速傅里叶变换,得到补偿后的子图像。
进一步地,所述步骤D中,补偿后的子图像表示为
hi′(x,y)=IFFT2{FFT2{hi(x,y)}·Hc(kx,ky;xi,yi)}。
本发明的有益效果是:本发明采用基于线性空变假设的波数谱一致化方法实现复杂飞行轨迹下的SAR回波一致聚焦,并通过高阶空变相位补偿提高有效聚焦景深,通过这两种处理技术的结合实现了复杂飞行轨迹下的双基SAR大场景成像处理;本发明不仅克服了现有SAR成像算法对复杂轨迹适应性差的缺陷,而且能够以高的运算效率实现聚束模式的大场景高分辨率一致聚焦,可以在较大的成像场景范围内,实现高分辨率的成像处理,使搭载双基SAR系统的飞行器在战场侦察监视、物资空投和地震灾害救援等领域发挥更好的性能。
附图说明
图1是本发明的适用于复杂飞行轨迹的双基聚束SAR大场景成像方法的流程示意图;
图2是本发明实施例中采用的双基SAR成像系统几何示意图;
图3是本发明实施例中采用的目标场景布置图;
图4是本发明实施例中对目标场景的成像仿真结果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
为了方便描述本发明的内容,首先对以下术语进行解释:
术语1:双基地合成孔径雷达(Bistatic Synthetic Aperture Radar,BiSAR)
双基地合成孔径雷达是指在雷达平台运动过程中,发射站天线对成像区域进行照射、接收站天线接收成像区域中的目标散射回波;利用发射信号的大带宽形成距离向高分辨,通过成像处理算法补偿方位向信号的多普勒相位以实现方位向孔径合成进而形成方位向高分辨,从而实现成像区域内的两维高分辨成像。
术语2:点目标参考函数(Point Target Reference Function)
点目标参考函数是指对回波数据进行一致压缩时采用的参考点的信号函数。参考函数可以是时域的,也可以是频域的,如在某个域生成了参考函数,成像处理时则将回波信号与参考函数在该域上共轭相乘,再变换到时频对应的另一个域上实现压缩。参考函数一般只包含相位项,不包含包络项。
术语3:有效成像场景尺寸
有效场景尺寸是指对某个特定的成像方法,在保证场景内所有位置目标回波的聚焦效果良好时可以允许的最大成像场景两维尺寸。SAR成像处理中,聚焦效果一般在场景边缘最差,散焦效应往往由在成像方法中被忽略的方位时间的二次相位误差(QuadraticPhase Error,QPE)导致,因此有效场景尺寸一般可以通过对QPE的数值加以限制来求出。
如图1所示,为本发明的适用于复杂飞行轨迹的双基聚束SAR大场景成像方法的流程示意图。一种适用于复杂飞行轨迹的双基聚束SAR大场景成像方法,包括以下步骤:
A、成像系统参数初始化
所述成像系统参数包括参考点位置、发射站三维轨迹[xT(t),yT(t),zT(t)]、接收站三维轨迹[xR(t),yR(t),zR(t)]、发射信号载波频率fc、脉冲宽度Tp、距离向采样频率Fs、距离向采样点数Nr、脉冲重复间隔PRI、方位时间向量t、距离时间向量τ、距离频率向量fτ、场景中点目标(xp,yp)的距离历史Rb(t;xp,yp);
B、对回波数据进行基于点目标参考函数的一致压缩
根据参考点目标的回波数据构造点目标参考函数,利用距离频域和方位时域对回波数据进行一致压缩,得到一致压缩后的距离频域数据;
C、对一致压缩后的距离频域数据做波数谱一致化变换
将差分距离史对空间坐标二维泰勒展开,根据双平台轨迹信息构造映射关系,通过二维插值将一致压缩后的数据变换到二维波数域,得到波数谱一致化的数据;
D、高阶空变相位补偿
将粗成像结果划分为多个子图像,构造空变相位补偿因子,计算每个子图像对应的相位补偿因子并进行相位补偿,将所有相位补偿后的子图像进行拼接得到成像结果。
在本发明的一个可选实施例中,上述步骤A初始化成像系统参数,包括参考点位置,定义为(0,0,0);发射站三维轨迹,记为[xT(t),yT(t),zT(t)];接收站三维轨迹,记为[xR(t),yR(t),zR(t)];发射信号载波频率,记为fc;脉冲宽度,记为Tp;距离向采样频率,记为Fs;距离向采样点数,记为Nr;脉冲重复间隔,记为PRI;方位时间向量,记为t=[-PRI·Na/2,-PRI·(Na/2-1),···,PRI·(Na/2-1)],距离时间向量,记为τ=[-1/Fs·Nr/2,-1/Fs·(Nr/2-1),···,1/Fs·(Nr/2-1)];距离频率向量,记为fτ=[-Fs/2,-Fs·(Nr/2-1)/Nr,···,Fs·(Nr/2-1)/Nr];
场景中点目标(xp,yp)的距离历史为:
Rb(t;xp,yp)=RT(t;xp,yp)+RR(t;xp,yp)
其中,RT(t;xp,yp)表示发射站(T)的距离历史,表示发射站(T)的轨迹矢量,RR(t;xp,yp)表示接收站(R)的距离历史,表示接收站(R)的轨迹矢量。
如图2所示,为本发明实施例中采用的双基SAR成像系统几何示意图。系统坐标系以雷达双站波束中心所指的地表点为坐标原点,平台沿弯曲轨迹运动。成像系统参数如表1所示。
表1、成像系统参数
如图3所示,为本发明实施例中采用的目标场景布置图;图中的圆点为布置于地面上7×7共49个点目标。这49个点沿x和y方向的间隔均为400m。场景中任一点的位置坐标记为P(x,y)。
构造方位时间向量t=[-PRI·Na/2,-PRI·(Na/2-1),···,PRI·(Na/2-1)],其中,PRI为脉冲重复时间,Na为目标回波方位向采样点数。距离时间向量,记为τ=[-1/Fs·Nr/2,-1/Fs·(Nr/2-1),···,1/Fs·(Nr/2-1)],其中,Fs为距离向采样率,Nr为目标回波距离向采样点数。
在本发明的一个可选实施例中,上述步骤B按照步骤A中初始化的雷达平台飞行轨迹和点目标坐标位置,利用MATLAB仿真出点目标回波数据sr(τ,t;xp,yp),表示为
其中,σp(xp,yp)表示目标散射截面积(RCS);wr[·]表示距离向时域窗函数,代表发射脉冲信号的包络,窗宽度为Tp;wa[·]表示方位向时域窗函数,代表天线方向图的调制效应,窗宽度为合成孔径时间Ta;Kr为发射信号调频率;c为光速;λ为载波波长。
根据参考点目标的回波数据,构造一致压缩参考信号,表示为
其中,So(fτ,t)表示点目标参考函数,So*(fτ,t)表示对一致压缩参考信号So(fτ,t)取共轭,ΔRp(t;xp,yp)表示差分距离史,。
再利用距离频域和方位时域将回波信号与参考信号进行共轭点乘,对回波数据进行一致压缩,压缩后的距离频域数据表示为:
其中,FFTr{·}表示距离向FFT,So*(fτ,t)表示对一致压缩参考信号So(fτ,t)取共轭,
其中,下角标0表示点目标位于中心点,Δ表示差分,ΔRT(t;xp,yp)表示发射站差分距离历史,ΔRR(t;xp,yp)表示接收站的差分距离历史。
在本发明的一个可选实施例中,上述步骤C将差分距离史ΔRb(t;xp,yp)对空间坐标(xp,yp)二维泰勒展开并将展开式分成线性项及高阶项,可以得到:
其中RT(t)=||rT(t)||2,RR(t)=||rR(t)||2;ΔRcur(t,xp,yp)一般保留至三阶,即:
再利用平台轨迹数据构造映射关系
其中,kx和ky分别为x和y方向上的波数。
根据映射关系,通过二维插值将一致压缩后的数据Sd(fτ,t;xp,yp)变换到由kx,ky定义的二维波数域,得到波数谱一致化的数据,表示为
其中,Wk(kx,ky)为经上述二维变换后的波数谱包络;为经二维变换后的波数域高阶空变相位;数据映射后的二维波数kx,ky分别为
其中,kxmin,kymin分别为x方向波数和y方向波数的最小值,kxmax,kymax分别为x方向波数和y方向波数的最大值,均可根据上文给出的波数映射关系求出;Dx,Dy分别为成像场景的x方向尺寸和y方向尺寸。
将整个场景中所有目标的回波信号进行处理后的的结果为
然后对该结果做二维逆向快速傅里叶变换,得到波数谱一致化的粗成像结果hcur(x,y)。
在本发明的一个可选实施例中,上述步骤D首先计算补偿前的有效场景尺寸,计算公式为
其中,Qi和Ti均为雷达平台运动轨迹相关系数,表示为:
求得关于xp和yp的取值范围后,将粗成像结果hcur(x,y)划分为若干尺寸由该取值范围确定的子图像hi(x,y)。
然后,根据每一个子图像的中心位置(xi,yi)构造空变相位补偿因子,表示为
Hc(kx,ky;xp,yp)=exp{jkr(kx,ky)Rcur[t(kx,ky);xp,yp]}
其中,kx和ky分别为x和y方向上的波数,kr(kx,ky)表示由kx和ky来表示kr的二维函数映射,t(kx,ky)表示由kx和ky来表示t的二维函数映射,Rcur表示距离误差。
t(kx,ky)定义为如下方程的解:
令系数A0,A1,A2,A3,...的求取公式为
A0=f(0)
A1=f'(0)
利用平台轨迹数据计算出上述数据向量并进行多项式拟合,得到系数A0,A1,A2,A3,...的值,再通过级数反演得到:
而kr(kx,ky)定义为:
对每个子图像hi(x,y),先确定其几何中心位置(xi,yi),将子图像进行二维快速傅里叶变换,得到其波数谱Hi(kx,ky);再计算出对应的相位补偿因子Hc(kx,ky;xi,yi),与子图像波数谱相乘后再求解乘积结果的二维逆向快速傅里叶变换,得到补偿后的子图像,表示为
hi′(x,y)=IFFT2{FFT2{hi(x,y)}·Hc(kx,ky;xi,yi)}。
所有子图像完成相位补偿后,将所有相位补偿后的子图像按照几何关系进行图像拼接,得到大场景一致聚焦的成像结果h′(x,y)。
如图4所示,为本发明实施例中对目标场景的成像仿真结果示意图;可以看出,本发明在较大的成像场景范围内,实现了高分辨率一致聚焦的成像处理。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (8)
1.一种适用于复杂飞行轨迹的双基聚束SAR大场景成像方法,其特征在于,包括以下步骤:
A、成像系统参数初始化
所述成像系统参数包括参考点位置、发射站三维轨迹[xT(t),yT(t),zT(t)]、接收站三维轨迹[xR(t),yR(t),zR(t)]、发射信号载波频率fc、脉冲宽度Tp、距离向采样频率Fs、距离向采样点数Nr、脉冲重复间隔PRI、方位时间向量t、距离时间向量τ、距离频率向量fτ、场景中点目标(xp,yp)的距离历史Rb(t;xp,yp);
B、对回波数据进行基于点目标参考函数的一致压缩
根据参考点目标的回波数据构造点目标参考函数,利用距离频域和方位时域对回波数据进行一致压缩,得到一致压缩后的距离频域数据;
C、对一致压缩后的距离频域数据做波数谱一致化变换
将差分距离史对空间坐标二维泰勒展开,根据双平台轨迹信息构造映射关系,通过二维插值将一致压缩后的数据变换到二维波数域,得到波数谱一致化的数据;
D、高阶空变相位补偿
将粗成像结果划分为多个子图像,构造空变相位补偿因子,计算每个子图像对应的相位补偿因子并进行相位补偿,将所有相位补偿后的子图像进行拼接得到成像结果。
2.如权利要求1所述的适用于复杂飞行轨迹的双基聚束SAR大场景成像方法,其特征在于,所述步骤B中,点目标参考函数表示为
其中,So(fτ,t)表示点目标参考函数,Kr为发射信号调频率,c为光速。
3.如权利要求2所述的适用于复杂飞行轨迹的双基聚束SAR大场景成像方法,其特征在于,所述步骤B中,利用距离频域、方位时域的共轭相乘,对回波数据进行一致压缩,一致压缩后的距离频域数据表示为
其中,FFTr{·}表示距离向FFT,sr(τ,t;xp,yp)表示目标回波数据,So *(fτ,t)表示对一致压缩参考信号So(fτ,t)取共轭,ΔRp(t;xp,yp)表示差分距离史,wr[·]表示距离向时域窗函数,wa[·]表示方位向时域窗函数,Ta表示合成孔径时间。
4.如权利要求3所述的适用于复杂飞行轨迹的双基聚束SAR大场景成像方法,其特征在于,所述步骤C中,将差分距离史对空间坐标二维泰勒展开,表示为
其中,Rcur(t,xp,yp)表示距离误差。
5.如权利要求4所述的适用于复杂飞行轨迹的双基聚束SAR大场景成像方法,其特征在于,所述步骤D中,首先计算补偿前的有效场景尺寸,计算公式为
其中,Qi和Ti均为雷达平台运动轨迹相关系数,λ表示载波波长;
再根据有效场景尺寸将粗成像结果hcur(x,y)划分为多个子图像hi(x,y)。
6.如权利要求5所述的适用于复杂飞行轨迹的双基聚束SAR大场景成像方法,其特征在于,所述步骤D中,构造空变相位补偿因子,表示为
Hc(kx,ky;xp,yp)=exp{jkr(kx,ky)Rcur[t(kx,ky);xp,yp]}
其中,kx和ky分别为x和y方向上的波数,kr(kx,ky)表示由kx和ky来表示kr的二维函数映射,t(kx,ky)表示由kx和ky来表示t的二维函数映射。
7.如权利要求6所述的适用于复杂飞行轨迹的双基聚束SAR大场景成像方法,其特征在于,所述步骤D中,对每个子图像hi(x,y),先确定其几何中心位置(xi,yi),将子图像进行二维快速傅里叶变换,得到其波数谱Hi(kx,ky);再计算出对应的相位补偿因子Hc(kx,ky;xi,yi),与子图像波数谱相乘后再求解乘积结果的二维逆向快速傅里叶变换,得到补偿后的子图像。
8.如权利要求7所述的适用于复杂飞行轨迹的双基聚束SAR大场景成像方法,其特征在于,所述步骤D中,补偿后的子图像表示为
h′i(x,y)=IFFT2{FFT2{hi(x,y)}·Hc(kx,ky;xi,yi)}。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811007730.9A CN109143236B (zh) | 2018-08-31 | 2018-08-31 | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811007730.9A CN109143236B (zh) | 2018-08-31 | 2018-08-31 | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109143236A true CN109143236A (zh) | 2019-01-04 |
CN109143236B CN109143236B (zh) | 2020-09-08 |
Family
ID=64825787
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811007730.9A Active CN109143236B (zh) | 2018-08-31 | 2018-08-31 | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109143236B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109959932A (zh) * | 2019-04-08 | 2019-07-02 | 西安电子科技大学 | 基于下降段曲线轨迹的雷达前视三维成像方法 |
CN114325705A (zh) * | 2021-12-31 | 2022-04-12 | 电子科技大学 | 一种高低轨双基地合成孔径雷达频域快速成像方法 |
CN116930956A (zh) * | 2023-09-19 | 2023-10-24 | 南京隼眼电子科技有限公司 | 基于目标尺寸的目标轨迹拼接方法、装置及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054190A (zh) * | 2016-07-21 | 2016-10-26 | 电子科技大学 | 基于频谱优化建模的双基地前视sar频域成像方法 |
CN106990396A (zh) * | 2017-05-25 | 2017-07-28 | 电子科技大学 | 一种双基地合成孔径雷达定位误差计算方法 |
CN108051809A (zh) * | 2017-11-14 | 2018-05-18 | 石家庄铁道大学 | 基于Radon变换的运动目标成像方法、装置及电子设备 |
CN108459321A (zh) * | 2018-02-07 | 2018-08-28 | 杭州电子科技大学 | 基于距离-方位圆模型的大斜视高分辨率sar成像方法 |
-
2018
- 2018-08-31 CN CN201811007730.9A patent/CN109143236B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054190A (zh) * | 2016-07-21 | 2016-10-26 | 电子科技大学 | 基于频谱优化建模的双基地前视sar频域成像方法 |
CN106990396A (zh) * | 2017-05-25 | 2017-07-28 | 电子科技大学 | 一种双基地合成孔径雷达定位误差计算方法 |
CN108051809A (zh) * | 2017-11-14 | 2018-05-18 | 石家庄铁道大学 | 基于Radon变换的运动目标成像方法、装置及电子设备 |
CN108459321A (zh) * | 2018-02-07 | 2018-08-28 | 杭州电子科技大学 | 基于距离-方位圆模型的大斜视高分辨率sar成像方法 |
Non-Patent Citations (1)
Title |
---|
李中余: "双基地合成孔径雷达动目标检测与成像技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109959932A (zh) * | 2019-04-08 | 2019-07-02 | 西安电子科技大学 | 基于下降段曲线轨迹的雷达前视三维成像方法 |
CN114325705A (zh) * | 2021-12-31 | 2022-04-12 | 电子科技大学 | 一种高低轨双基地合成孔径雷达频域快速成像方法 |
CN114325705B (zh) * | 2021-12-31 | 2023-09-22 | 电子科技大学 | 一种高低轨双基地合成孔径雷达频域快速成像方法 |
CN116930956A (zh) * | 2023-09-19 | 2023-10-24 | 南京隼眼电子科技有限公司 | 基于目标尺寸的目标轨迹拼接方法、装置及存储介质 |
CN116930956B (zh) * | 2023-09-19 | 2023-11-17 | 南京隼眼电子科技有限公司 | 基于目标尺寸的目标轨迹拼接方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109143236B (zh) | 2020-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | A frequency-domain imaging algorithm for highly squinted SAR mounted on maneuvering platforms with nonlinear trajectory | |
CN102393518B (zh) | 一种适用于大斜视角的机载sar成像方法 | |
CN106093870B (zh) | 高超声速飞行器下降段的sar-gmti杂波抑制方法 | |
CN109143237B (zh) | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 | |
Park et al. | Cross-range scaling algorithm for ISAR images using 2-D Fourier transform and polar mapping | |
CN104833972B (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
CN104977582A (zh) | 一种实现扫描雷达方位超分辨成像的解卷积方法 | |
CN109143236B (zh) | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 | |
CN103336278A (zh) | 多视角观测下前视三维sar成像方法 | |
CN104808204A (zh) | 一种固定发射站双基前视sar动目标检测方法与成像方法 | |
CN104749570A (zh) | 一种移不变机载双基合成孔径雷达目标定位方法 | |
CN104330779A (zh) | 机载合成孔径雷达运动误差补偿方法 | |
CN109031299B (zh) | 低信噪比条件下基于相位差分的isar平动补偿方法 | |
Qiu et al. | Bistatic SAR data processing algorithms | |
CN105487074A (zh) | 一种双基地合成孔径雷达数值距离多普勒成像方法 | |
Ren et al. | 2-D spatially variant motion error compensation for high-resolution airborne SAR based on range-Doppler expansion approach | |
Jao et al. | Multichannel synthetic aperture radar signatures and imaging of a moving target | |
CN103091682B (zh) | 基于时频分析InISAR多动目标成像和运动轨迹重建法 | |
Rahman | Focusing moving targets using range migration algorithm in ultra wideband low frequency synthetic aperture radar | |
CN109188436B (zh) | 适用于任意平台轨迹的高效双基sar回波生成方法 | |
CN114325704A (zh) | 一种基于波数谱拼接的合成孔径雷达快速时域成像方法 | |
CN106526544B (zh) | 基于高超声速平台的mimosar杂波抑制方法 | |
Xu et al. | A variable PRF imaging method for high squint diving SAR | |
Xue et al. | An applied frequency scaling algorithm based on local stretch factor for near-field miniature millimeter-wave radar imaging |
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 |