CN101776755A - 高速平台超高分辨率sar数据预处理方法 - Google Patents

高速平台超高分辨率sar数据预处理方法 Download PDF

Info

Publication number
CN101776755A
CN101776755A CN201010107189A CN201010107189A CN101776755A CN 101776755 A CN101776755 A CN 101776755A CN 201010107189 A CN201010107189 A CN 201010107189A CN 201010107189 A CN201010107189 A CN 201010107189A CN 101776755 A CN101776755 A CN 101776755A
Authority
CN
China
Prior art keywords
sar
walk
phase
echo data
orientation
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
CN201010107189A
Other languages
English (en)
Other versions
CN101776755B (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.)
Xian Cetc Xidian University Radar Technology Collaborative Innovation Research Institute Co Ltd
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN2010101071896A priority Critical patent/CN101776755B/zh
Publication of CN101776755A publication Critical patent/CN101776755A/zh
Application granted granted Critical
Publication of CN101776755B publication Critical patent/CN101776755B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种高速平台超高分辨率SAR数据的预处理方法,主要克服现有传统SAR成像技术对于高速平台超高分辨率SAR回波成像精度不高的问题,其实现过程是:首先判断SAR回波数据方位频谱是否混叠,对方位频谱混叠的SAR回波进行方位解模糊处理;接着对SAR回波数据进行“走停运动”假设判断,对不满足“走停运动”假设判断式的SAR回波数据进行相位补偿,将其转化为基于“走停运动”假设构建的SAR回波。本发明扩展了传统SAR成像技术的适用范围,经过本发明预处理,使高速平台超高分辨率SAR回波数据可以继续采用基于“走停运动”假设的传统SAR成像技术进行成像,可应用于星载、弹载以及超高音速无人机超高分辨率SAR成像处理。

Description

高速平台超高分辨率SAR数据预处理方法
技术领域
本发明属于雷达信号处理领域,涉及高速平台超高分辨率SAR数据的预处理方法,用于运动平台为卫星、导弹以及超高音速无人机的高速运动平台超高分辨率SAR成像。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)成像属于高分辨成像技术,在距离方向,它通过发射大带宽信号获得高分辨;在方位方向,用一个运动的传感器模拟一个长的天线阵列,得到在该方向上较好的分辨率。近半个世纪以来,SAR成像理论和技术有了很大发展,并得到了广泛的应用,如遥感和地图测绘。
传统SAR成像方法包括用于超宽带SAR数据的后向投影算法、距离徙动算法,和进行分维处理的距离多普勒算法、chirp sealing算法、频率变标算法以及精度更高的分维算法如非线性调频变标NCS算法等,大多数是基于对SAR的“走停运动”假设。“走停运动”假设即SAR在一个位置停下来发射、接收信号,然后移向另一个位置。基于“走停运动”假设构建的SAR回波信号仅仅是发射信号的简单时延。然而雷达与目标相对运动,使得SAR回波信号的频率较发射信号而言发生了变化,这种变化被称为多普勒频移。这使得真实的SAR回波信号不仅仅是发射信号的简单时延,其信号形式发生了变化。信号形式的变化使得基于“走停运动”假设构建的SAR回波信号相位不同于真实SAR回波信号相位。这种相位偏差的存在会影响到传统SAR成像方法的效果。对于窄带SAR回波信号而言,这种相位偏差较小,对成像的影响是可以忽略的。SAR未来向高精度高分辨发展的趋势必然要求雷达发射信号提高发射带宽。对于分辨率在厘米量级的超高分辨率SAR而言,尤其当运动平台速度很高时,相位偏差增大,相位偏差对成像的影响无法忽略,此时基于“走停运动”假设的传统SAR成像方法不再有效。
国内外对不满足“走停运动”假设的SAR数据处理方法的研究有很长的历史。现有的研究主要集中于对调频连续波SAR数据处理方法的研究。此类方法假设雷达信号传输过程中雷达载体与目标之间相对静止,适用于雷达作用距离比较近的情况,例如机载SAR。而当雷达作用距离比较远时,如超高分辨率星载SAR,这种假设是不适用的。由于高速平台超高分辨率SAR的雷达作用距离较远,从而该类处理方法对于高速平台超高分辨率SAR数据是不适用的。由于现有SAR成像方法不再有效,因此需要对高速平台超高分辨率SAR数据的处理方法进行研究。
发明内容
本发明的目的是提供一种高速平台超高分辨率SAR数据预处理方法,以克服传统SAR成像技术对于高速平台超高分辨率SAR回波数据成像精度不高,以及现有调频连续波SAR数据处理方法不适用于高速平台超高分辨率SAR回波数据处理的问题。
本发明是这样实现的:
在运动平台速度高至上千米每小时的情况下,对于距离向和方位向频带宽度都很宽的脉冲式超高分辨率SAR,其回波不再仅仅是发射信号的简单时延,基于“走停运动”假设的传统脉冲式SAR回波模型不足以精确的描述回波信号,从而建立在传统脉冲式SAR回波模型基础上的成像方法不再有效。本发明给出了一种高速运动平台超高分辨率SAR数据预处理方法,该方法将原始SAR回波信号补偿为基于“走停运动”假设构建的SAR回波信号,从而使传统脉冲式SAR成像技术可以继续应用于高速平台超高分辨率SAR成像。具体过程包括如下:
(1)接收高速平台超高分辨率SAR回波数据,并采用正交解调法对该回波数据进行去载频处理;
(2)判断SAR回波数据方位频谱是否混叠,当方位多普勒带宽Ba大于方位采样率PRF时,则判为SAR回波数据存在方位频谱混叠,转入步骤3,反之转入步骤4;
(3)将SAR回波数据与一个随距离向频率变化的二次相位参考函数作卷积操作,完成对SAR回波数据的方位解模糊处理,该二次相位参考函数表达式为:
S ref = exp ( j 2 π ( vt m ) 2 R ref ( f c + f r ) c ) ,
其中:fc为发射信号的载频,fr为基带距离向频率,v为SAR运动平台有效速度,tm为方位时间,Rref为场景中心线到雷达的最近距离,c为光速;
(4)利用“走停运动”假设判断式确定SAR回波信号是否满足“走停运动”假设,若SAR回波数据满足构建的“走停运动”假设判断式,则认为该高速平台超高分辨率SAR系统符合“走停运动”假设,结束预处理,反之则需要将SAR回波数据补偿为基于“走停运动”假设构建的SAR回波信号,执行步骤(5),该“走停运动”假设判断式为:
π BT P v r _ max c ≤ π 4 v 2 c ( f c γ R ref + 1 c ) BT a ≤ 1 4 ,
其中:B为发射脉冲的带宽,TP为发射脉冲的时宽,vr_max表示目标与雷达的最大径向速度,Ta为合成孔径的时宽,γ为发射脉冲的调频率;
(5)在距离频域/方位频域,将SAR回波数据与构建的第一次相位补偿参考函数进行矩阵点乘,完成对SAR回波数据的第一次相位补偿,
该第一次相位补偿参考函数如下:
H c 1 = exp ( jπ 1 γ ( f d + 2 ( f c + f r ) ( v c ) 2 ) 2 - j 2 π 1 γ f r ( f d + 2 ( f c + f r ) ( v c ) 2 ) ) ,
其中:fd为多普勒频率;
(6)将进行了第一次相位补偿后的SAR回波数据进行方位向逆傅里叶变换,将其变换到距离频域/方位时域,并在距离频域/方位时域,将SAR回波数据与构建的第二次相位补偿参考函数进行矩阵点乘,完成对SAR回波数据的第二次相位补偿,得到满足“走停运动”假设的SAR回波数据,该第二次相位补偿参考函数如下:
H c 2 = exp ( j 4 πt m ( v c ) 2 ( f c + f r ) ) .
本发明具有如下优点:
1)本发明由于通过相位补偿操作将不满足“走停运动”假设的高速运动平台超高分辨率SAR回波数据补偿为满足“走停运动”假设的SAR回波数据,从而能利用现有的基于“走停运动”假设的传统脉冲式SAR成像技术对高速运动平台超高分辨率SAR数据进行成像;
2)本发明由于采用对SAR回波数据进行方位解模糊处理,消除了由于高速平台超高分辨率SAR低脉冲重复周期导致的图像重影现象,最终可以得到高分辨率的清晰图像。
仿真结果表明,高速平台超高分辨率SAR回波数据通过本发明的预处理,可以继续使用传统脉冲式SAR成像技术进行成像,并且具有较好的成像效果。
以下结合附图和具体实施方式对本发明的技术方案作进一步详细描述。
附图说明
图1是本发明的SAR数据预处理流程图;
图2是本发明SAR数据预处理中相位补偿子流程图;
图3是未采用本发明预处理的SAR回波与现有基于SAR“走停运动”假设构建的SAR回波的相位差的仿真图;
图4是用本发明预处理后的SAR回波与现有基于“走停运动”假设构建的SAR回波的相位差的仿真图;
图5是采用传统距离徙动成像方法对未采用本发明预处理的SAR回波信号进行成像的结果图;
图6是采用传统距离徙动成像方法对通过本发明预处理的SAR回波信号进行成像的结果图。
具体实施方式
参照图1,本发明对SAR数据的预处理步骤如下:
步骤1,接收SAR回波数据,并采用现有正交解调法对该回波数据进行去载频处理。
步骤2,判断SAR回波数据方位频谱是否混叠,当方位多普勒带宽Ba大于方位采样率PRF时,则判为SAR回波数据存在方位频谱混叠,转入步骤3,反之转入步骤4。
步骤3,对SAR回波数据进行方位解模糊处理;
3.1通过将SAR回波数据与随距离向频率变化的二次相位参考函数Sref作卷积,实现SAR回波方位频谱不混叠;
该二次相位参考函数表达式为:
S ref = exp ( j 2 π ( vt m ) 2 R ref ( f c + f r ) c ) - - - 1 )
其中:fc为发射信号的载频,fr为基带距离向频率,v为SAR运动平台有效速度,tm为方位时间,Rref为场景中心线到雷达的最近距离,c为光速;
该卷积操作表示为:
S out ( f r , n )
= Σ i = - P / 2 P / 2 - 1 S in ( f r , iΔ t in ) S ref ( f r , nΔ t out - i Δt in ) - - - 2 )
= exp ( j 2 π ( nΔt out ) 2 R ref ( f c + f r ) c ) DFT [ S in ( f r , iΔ t in ) exp ( j 2 π ( nΔ t in ) 2 R ref ( f c + f r ) c ) ]
其中:Sin和Sout分别表示输入SAR回波信号和方位解模糊处理后的输出SAR回波数据,i为输入方位时间采样序数,n为输出方位时间采样序数,n=-Pmax/2,...,Pmax/2-1,Pmax为方位解模糊处理后的SAR回波数据方位采样点数,Δtin为接收的SAR回波数据在方位采样间隔,Δtout为方位解模糊处理后的SAR回波数据的方位采样间隔,fr为基带距离向频率,DFT[·]为离散傅立叶变换;
3.2为了补偿上述卷积操作对SAR回波数据的改变,在距离频域/方位频域,将通过卷积操作得到的SAR回波数据与随距离向频率变化参考函数Scom进行矩阵点乘,该参考函数Scom表示为:
S com = exp ( - jπ ( f d v ) 2 cR ref 2 ( f c + f r ) ) , - - - 3 )
其中:fd为多普勒频率。
步骤4,构建“走停运动”假设判断式,对SAR回波数据进行“走停运动”假设判断;
4.1构建符合SAR真实运动的SAR回波模型的相位为:
Figure GSA00000036615200055
其中:
τ = 2 c ( ( vt m - X N ) + v t ^ r ) 2 + R N 2 - 2 v ( ( vt m - X N ) + v t ^ r ) c 2 - v 2
其中:
Figure GSA00000036615200057
为接收信号时刻的快时间,XN为目标沿平台运动方向的坐标位置,RN为雷达与目标的最近距离;
4.2根据高速运动平台超高分辨SAR系统的参数范围对上述4)式的相位进行近似,构建基于高速运动平台超高分辨SAR回波模型的相位
Figure GSA00000036615200061
Figure GSA00000036615200062
其中:
τ · 0 = 2 cv c 2 - v 2 ( vt m - X N ) ( vt m - X N ) 2 + R N 2 + 2 v 2 c 2 - v 2
τ 0 = 2 c c 2 - v 2 ( vt m - X N ) 2 + R N 2 + 2 V c 2 - v 2 ( vt m - X N ) ;
4.3对基于高速平台超高分辨SAR回波模型构建的SAR回波(其相位为上述5)式所示)去载频,得到基于高速运动平台超高分辨SAR回波模型构建的SAR回波信号在距离频域/方位时域的相位为:
Figure GSA00000036615200065
4.4对基于“走停运动”假设构建的SAR回波信号去载频,得到基于“走停运动”假设构建的SAR回波信号在距离频域/方位时域的相位为:
Figure GSA00000036615200066
其中:
τ trad = 2 c ( vt m - X N ) 2 + R N 2 ;
4.5将步骤4.3中的基于高速平台超高分辨SAR回波模型构建的SAR回波信号相位
Figure GSA00000036615200068
与步骤4.4中基于“走停运动”假设构建的SAR回波信号相位
Figure GSA00000036615200069
的共轭进行相乘,得到基于高速平台超高分辨SAR回波模型构建的SAR回波信号与基于“走停运动”假设方法构建的SAR回波信号的相位差为:
Figure GSA000000366152000610
4.6将上述8)式相位差
Figure GSA000000366152000611
中关于fr的二次项表示为:
π 1 γ ( τ · 0 2 + 2 τ · 0 ) f r 2 , - - - 9 )
这一项的存在会导致采用基于“走停运动”假设的传统距离脉压后,SAR回波数据的脉压结果中存在残余线性频率调制,用残余二次相位差的最大值来衡量残余线性频率调制给成像带来的影响,该残余二次相位差的最大值表示为:
QPE max = π BT P v r _ max c ; - - - 10 )
其中:vr_max表示目标与雷达的最大径向速度,当QPEmax小于π/4时,残余线性频率调制对成像效果的影响可以忽略,从而得到“走停”假设判断式中的第一个判断式:
π BT P v r _ max c ≤ π 4 ; - - - 11 )
其中:B为发射脉冲的带宽,TP为发射脉冲的时宽,vr_max表示目标与雷达的最大径向速度;
4.7将8)式相位差
Figure GSA00000036615200073
中关于fr的一次项表示为:
2 π ( τ 0 - τ trad + τ · 0 ( τ · 0 + 1 ) γ f c ) f r , - - - 12 )
该项的存在使得采用基于传统“走停”假设的距离徙动校正后,SAR回波数据仍会存在线性走动,当线性走动超过二分之一个距离分辨单元时,该线性走动对成像效果的影响可以忽略,得到“走停”假设的第二个判断式:
v 2 c ( f c γ R ref + 1 c ) BT a ≤ 1 4 , - - - 13 )
其中:Ta为合成孔径的时宽,γ为发射脉冲的调频率;
若SAR回波数据同时满足构建的两个“走停运动”假设判断式11)和13),则认为SAR回波数据与基于“走停运动”假设构建的SAR回波信号的差异可以忽略,并结束预处理,反之则需要将SAR回波数据补偿为基于“走停运动”假设构建的SAR回波信号,转入步骤5。
步骤5,对SAR回波数据进行相位补偿处理;
参照图2,相位补偿处理的具体过程如下:
在距离频域/方位频域,将SAR回波数据与构建的第一次相位补偿参考函数进行矩阵点乘,完成对SAR回波数据的第一次相位补偿;将进行了第一次相位补偿后的SAR回波数据进行方位向IFFT,将其变换到距离频域/方位时域,并在距离频域/方位时域,将SAR回波数据与构建的第二次相位补偿参考函数进行矩阵点乘,完成对SAR回波数据的第二次相位补偿,得到满足“走停运动”假设的SAR回波数据。其中,两次相位补偿参考函数,按如下过程进行构建:
第一次相位补偿参考函数的构建:
5.1在距离方位二维频域,将基于高速运动平台超高分辨SAR回波模型构建的SAR回波信号的相位表示为:
Figure GSA00000036615200081
                    (         14)
Figure GSA00000036615200082
5.2在距离方位二维频域,将基于“走停运动”假设方法构建的SAR回波信号的相位表示为:
Figure GSA00000036615200084
5.3在距离方位二维频域,将步骤5.1中基于高速平台超高分辨SAR回波模型构建的SAR回波信号的相位与步骤5.2中基于“走停运动”假设方法构建的SAR回波信号的相位
Figure GSA00000036615200086
进行比较,两个相位在根式外的差异为:
1 γ ( f d + 2 ( f c + f r ) ( v c ) 2 ) 2 - 2 π 1 γ f r ( f d + 2 ( f c + f r ) ( v c ) 2 ) ; - - - 16 )
5.4将步骤5.3中的两个相位在根式外的差异表示为信号形式Hc1,并将Hc1作为第一次相位补偿参考函数:
H c 1 = exp ( jπ 1 γ ( f d + 2 ( f c + f r ) ( v c ) 2 ) 2 - j 2 π 1 γ f r ( f d + 2 ( f c + f r ) ( v c ) 2 ) ) ; - - - 17 )
第二次相位补偿参考函数的构建:
5.5在距离频域/方位频域,将基于高速运动平台超高分辨SAR回波模型构建的SAR回波信号与第一次相位补偿参考函数进行矩阵点乘,矩阵点乘后的输出信号相位为:
Figure GSA00000036615200091
5.6在距离频域/方位时域,将上式的相位
Figure GSA00000036615200092
与式7)所示的基于“走停运动”假设构建的SAR回波相位
Figure GSA00000036615200093
进行比较,两个相位的差异为:
4 π t m ( v c ) 2 ( f r + f c ) ; - - - 19 )
5.7将步骤5.6中两个相位的差异表示为信号形式Hc2,并将Hc2作为第二次相位补偿参考函数:
H c 2 = exp ( j 4 π t m ( v c ) 2 ( f r + f c ) ) ; - - - 20 )
上述所有公式中,相同符号的含义均为相同。
至此,高速平台超高分辨率SAR回波数据已经转化为方位频谱不混叠的满足“走停运动”假设的SAR回波。
以下通过仿真实验进一步说明本发明的优点。
1.仿真参数
仿真参数为表1:
表1星载的超高分辨率SAR参数
Figure GSA00000036615200096
Figure GSA00000036615200101
2仿真内容及结果
仿真内容1:
仿真本发明预处理前基于SAR真实运动模式构建的星载超高分辨率SAR点目标回波信号与基于SAR“走停运动”假设构建的星载超高分辨率SAR点目标回波信号的相位差,仿真结果如图3所示,其中:图中最右边的灰度条表示相位(从-π到π),可以看出基于SAR真实运动模式构建的SAR回波与基于“走停运动”假设构建的SAR回波相位差较大,说明基于SAR真实运动模式构建的SAR回波不满足“走停运动”假设。
仿真采用本发明预处理后基于SAR真实运动模式构建的星载超高分辨率SAR点目标回波信号与基于SAR“走停运动”假设构建的星载超高分辨率SAR点目标回波信号的相位差,仿真结果如图4所示,其中:图中最右边的灰度条表示相位(从-π到π),可以看出,相位差接近于零,说明通过本发明的预处理,可将高速平台超高分辨率SAR回波补偿为基于“走停运动”假设构建的SAR回波。
仿真内容2:
仿真在本发明预处理前,采用传统距离徙动成像方法对基于SAR真实运动模式构建的星载超高分辨率SAR点目标回波进行成像,仿真的成像结果如图5所示,其中:成像结果出现较为严重的散焦,目标无法聚焦在一个采样单元,说明基于“走停运动”假设的传统成像方法对该星载超高分辨率SAR仿真回波数据的成像效果较差。
仿真在本发明预处理后,采用传统距离徙动成像方法对基于SAR真实运动模式构建的星载超高分辨率SAR点目标回波进行成像,仿真的成像结果如图6所示,其中:目标能够聚焦在一个采样单元,成像效果较好,说明通过本发明的预处理,可以继续采用传统成像方法对高速平台超高分辨率SAR回波进行成像。

Claims (4)

1.一种高速平台超高分辨率SAR数据预处理方法,包括如下过程:
(1)接收高速平台超高分辨率SAR回波数据,并采用正交解调法对该回波数据进行去载频处理;
(2)判断SAR回波数据方位频谱是否混叠,当方位多普勒带宽Ba大于方位采样率PRF时,则判为SAR回波数据存在方位频谱混叠,转入步骤3,反之转入步骤4;
(3)将SAR回波数据与一个随距离向频率变化的二次相位参考函数作卷积操作,完成对SAR回波数据的方位解模糊处理,该二次相位参考函数表达式为:
S ref = exp ( j 2 π ( vt m ) 2 R ref ( f c + f r ) c ) ,
其中:fc为发射信号的载频,fr为基带距离向频率,v为SAR运动平台有效速度,tm为方位时间,Rref为场景中心线到雷达的最近距离,c为光速;
(4)利用“走停运动”假设判断式确定SAR回波信号是否满足“走停运动”假设,若SAR回波数据满足构建的“走停运动”假设判断式,则认为该高速平台超高分辨率SAR系统符合“走停运动”假设,结束预处理,反之则需要将SAR回波数据补偿为基于“走停运动”假设构建的SAR回波信号,执行步骤(5),该“走停运动”假设判断式为:
π BT P v r _ max c ≤ π 4 v 2 c ( f c γ R ref + 1 c ) BT a ≤ 1 4 ,
其中:B为发射脉冲的带宽,TP为发射脉冲的时宽,vr_max表示目标与雷达的最大径向速度,Ta为合成孔径的时宽,γ为发射脉冲的调频率;
(5)在距离频域/方位频域,将SAR回波数据与构建的第一次相位补偿参考函数进行矩阵点乘,完成对SAR回波数据的第一次相位补偿,
该第一次相位补偿参考函数如下:
H c 1 = exp ( jπ 1 γ ( f d + 2 ( f c + f r ) ( v c ) 2 ) 2 - j 2 π 1 γ f r ( f d + 2 ( f c + f r ) ( v c ) 2 ) ) ,
 其中:fd为多普勒频率;
(6)将进行了第一次相位补偿后的SAR回波数据进行方位向逆傅里叶变换,将其变换到距离频域/方位时域,并在距离频域/方位时域,将SAR回波数据与构建的第二次相位补偿参考函数进行矩阵点乘,完成对SAR回波数据的第二次相位补偿,得到满足“走停运动”假设的SAR回波数据,该第二次相位补偿参考函数如下:
H c 2 = exp ( j 4 πt m ( v c ) 2 ( f c + f r ) ) .
2.根据权利要求1中所述的高速平台超高分辨率SAR数据预处理方法,其中步骤(4)所述的“走停运动”假设判断式,按如下过程构建:
(2a)构建符合SAR真实运动的SAR回波模型的相位为:
Figure FSA00000036615100023
τ = 2 c ( ( vt m - X N ) + v t ^ r ) 2 + R N 2 - 2 v ( ( vt m - X N ) + v t ^ r ) c 2 - v 2
其中:
Figure FSA00000036615100025
为接收信号时刻的快时间,XN为目标沿平台运动方向的坐标位置,RN为雷达与目标的最近距离;
(2b)根据高速平台超高分辨SAR系统的参数范围对上述符合SAR真实运动的SAR回波信号相位进行近似,构建基于高速平台超高分辨SAR回波模型的相位
Figure FSA00000036615100026
Figure FSA00000036615100027
其中:
τ · 0 = 2 cv c 2 - v 2 ( vt m - X N ) ( vt m - X N ) 2 + R N 2 + 2 v 2 c 2 - v 2 ,
τ 0 = 2 c c 2 - v 2 ( vt m - X N ) 2 + R N 2 + 2 v c 2 - v 2 ( vt m - X N ) ;
(2c)对基于高速平台超高分辨SAR回波模型构建的SAR回波信号去载频,得到基于高速运动平台超高分辨SAR回波模型构建的SAR回波信号在距离频域/方位时域的相位
Figure FSA00000036615100032
(2d)进行了去载频后,基于“走停运动”假设构建的SAR回波信号在距离频域/方位时域的相位
Figure FSA00000036615100033
表示为:
Figure FSA00000036615100034
其中:
τ trad = 2 c ( vt m - X N ) 2 + R N 2 ;
(2e)将步骤(2c)中的基于高速平台超高分辨SAR回波模型构建的SAR回波信号相位
Figure FSA00000036615100036
与步骤(2d)中基于“走停运动”假设构建的SAR回波信号相位
Figure FSA00000036615100037
的共轭进行相乘,得到基于高速平台超高分辨SAR回波模型构建的SAR回波信号与基于“走停运动”假设方法构建的SAR回波信号的相位差
Figure FSA00000036615100038
为:
Figure FSA00000036615100039
Figure FSA000000366151000310
(2f)将上式相位差
Figure FSA000000366151000311
中关于fr的二次项表示为:
π 1 γ ( τ · 0 2 + 2 τ · 0 ) f r 2 ,
这一项的存在会导致采用基于“走停运动”假设的传统距离脉压后,SAR回波数据的脉压结果中存在残余线性频率调制,用残余二次相位差的最大值来衡量残余线性频率调制给成像带来的影响,残余二次相位差的最大值QPEmax表示为:
QPE max = π BT P v r _ max c ;
当QPEmax小于π/4时,残余线性频率调制对成像效果的影响可以忽略,从而得到“走停”假设判断式中的第一个判断式:
π BT P v r _ max c ≤ π 4 ;
(2g)将步骤(2e)中相位差
Figure FSA00000036615100042
关于fr的一次项表示为:
2 π ( τ 0 - τ trad + τ · 0 ( τ · 0 + 1 ) γ f c ) f r ,
该项的存在使得采用基于“走停运动”假设的传统距离徙动校正后,SAR回波数据仍会存在线性走动,当线性走动超过二分之一个距离分辨单元时,该线性走动对成像效果的影响可以忽略,得到“走停运动”假设判断式中的第二个判断式:
v 2 c ( f c γ R ref + 1 c ) BT a ≤ 1 4 .
3.根据权利要求1中所述的高速平台超高分辨率SAR数据预处理方法,其特征在于步骤(5)中所述的第一次相位补偿参考函数,按如下过程进行构建:
(3a)在距离方位二维频域,将基于高速运动平台超高分辨SAR回波模型构建的SAR回波信号的相位表示为:
Figure FSA00000036615100047
(3b)在距离方位二维频域,基于“走停运动”假设方法构建的SAR回波信号的相位
Figure FSA00000036615100048
表示为:
Figure FSA00000036615100049
(3c)在距离方位二维频域,将步骤(3a)中基于高速平台超高分辨SAR回波模型构建的SAR回波信号的相位与步骤(3b)中基于“走停运动”假设方法构建的SAR回波信号的相位
Figure FSA000000366151000411
进行比较,两个相位在根式外的差异为:
π 1 γ ( f d + 2 ( f c + f r ) ( v c ) 2 ) 2 - 2 π 1 γ f r ( f d + 2 ( f c + f r ) ( v c ) 2 ) ;
(3d)将(3c)中的两个相位在根式外的差异表示为信号形式Hc1
H c 1 = exp ( jπ 1 γ ( f d + 2 ( f c + f r ) ( v c ) 2 ) 2 - j 2 π 1 γ f r ( f d + 2 ( f c + f r ) ( v c ) 2 ) ) ;
(3e)将(3d)中的信号形式Hc1作为第一次相位补偿参考函数,并用将该第一次相位补偿参考函数与去载频后的SAR回波数据在距离方位二维频域的进行矩阵点乘,完成对SAR回波数据的第一次相位补偿。
4.根据权利要求1中所述的高速平台超高分辨率SAR数据预处理方法,其特征在于步骤(6)中所述的第二次相位补偿参考函数,按如下过程进行构建:
(4a)进行了第一次相位补偿后,在距离频域/方位时域,将基于高速运动平台超高分辨SAR回波模型构建的SAR回波信号相位
Figure FSA00000036615100052
为:
Figure FSA00000036615100053
(4b)在距离频域/方位时域,将上式的相位
Figure FSA00000036615100054
与步骤(2d)中基于“走停运动”假设方法构建的SAR回波信号的相位进行比较,两个相位的差异为:
4 πt m ( v c ) 2 ( f r + f c ) ;
(4c)将(4b)中两个相位的差异表示为信号的形式Hc2
H c 2 = exp ( j 4 πt m ( v c ) 2 ( f r + f c ) ) ;
(4d)将(4c)中的信号形式Hc2作为第二次相位补偿参考函数,将该第二次相位补偿参考函数Hc2与第一次相位补偿后的SAR回波数据在距离频域/方位时域进行矩阵点乘,完成对SAR回波数据第二次相位补偿,得到与基于“走停运动”假设构建SAR回波信号一致的SAR回波信号。
CN2010101071896A 2010-02-05 2010-02-05 高速平台超高分辨率sar数据预处理方法 Active CN101776755B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101071896A CN101776755B (zh) 2010-02-05 2010-02-05 高速平台超高分辨率sar数据预处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101071896A CN101776755B (zh) 2010-02-05 2010-02-05 高速平台超高分辨率sar数据预处理方法

Publications (2)

Publication Number Publication Date
CN101776755A true CN101776755A (zh) 2010-07-14
CN101776755B CN101776755B (zh) 2012-03-28

Family

ID=42513263

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101071896A Active CN101776755B (zh) 2010-02-05 2010-02-05 高速平台超高分辨率sar数据预处理方法

Country Status (1)

Country Link
CN (1) CN101776755B (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102565784A (zh) * 2012-01-04 2012-07-11 西安电子科技大学 基于vsar系统的动目标重定位与速度解模糊方法
CN102680956A (zh) * 2012-05-14 2012-09-19 北京航空航天大学 星载合成孔径雷达方位向第一模糊区回波信号能量抑制系统
CN102707269A (zh) * 2012-06-14 2012-10-03 电子科技大学 一种机载雷达距离走动校正方法
CN102736080A (zh) * 2012-04-27 2012-10-17 中国人民解放军国防科学技术大学 一种用于探地雷达同心圆回波数据的成像方法
CN102914772A (zh) * 2012-09-18 2013-02-06 西安电子科技大学 基于等效散射点的进动目标二维成像方法
CN103197285A (zh) * 2013-03-22 2013-07-10 电子科技大学 一种用于合成孔径雷达成像的导航数据拟合方法
CN104076360A (zh) * 2014-07-04 2014-10-01 西安电子科技大学 基于压缩感知的二维sar稀疏目标成像方法
CN104133215A (zh) * 2014-05-29 2014-11-05 西安电子科技大学 基于距离徙动微调和子带分割的同步轨道雷达成像方法
CN104865571A (zh) * 2015-04-02 2015-08-26 西安电子科技大学 多通道多子带滑动聚束模式sar成像方法
CN105676220A (zh) * 2014-11-21 2016-06-15 中国航空工业集团公司雷华电子技术研究所 一种结合二维逆滤波的高分辨聚束sar自聚焦方法
CN109358329A (zh) * 2018-11-06 2019-02-19 电子科技大学 脉冲传播时间内平台机动的双基sar回波模型建立方法
CN111323761A (zh) * 2020-03-20 2020-06-23 北京华力创通科技股份有限公司 回波系统函数构建方法、装置及回波模拟器
CN111929681A (zh) * 2020-06-24 2020-11-13 苏州理工雷科传感技术有限公司 一种基于轻小型无人机载sar的实时成像处理系统
CN114488151A (zh) * 2022-04-08 2022-05-13 中国科学院空天信息创新研究院 针对观测船只的主被动联合探测方法、装置及设备和介质

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103576150B (zh) * 2013-09-24 2015-09-02 西安电子科技大学 基于高超声速飞行器俯冲段的前斜视sar成像方法

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102565784A (zh) * 2012-01-04 2012-07-11 西安电子科技大学 基于vsar系统的动目标重定位与速度解模糊方法
CN102565784B (zh) * 2012-01-04 2014-05-14 西安电子科技大学 基于vsar系统的动目标重定位与速度解模糊方法
CN102736080A (zh) * 2012-04-27 2012-10-17 中国人民解放军国防科学技术大学 一种用于探地雷达同心圆回波数据的成像方法
CN102736080B (zh) * 2012-04-27 2014-03-12 中国人民解放军国防科学技术大学 一种用于探地雷达同心圆回波数据的成像方法
CN102680956A (zh) * 2012-05-14 2012-09-19 北京航空航天大学 星载合成孔径雷达方位向第一模糊区回波信号能量抑制系统
CN102680956B (zh) * 2012-05-14 2013-08-14 北京航空航天大学 星载合成孔径雷达方位向第一模糊区回波信号能量抑制系统
CN102707269A (zh) * 2012-06-14 2012-10-03 电子科技大学 一种机载雷达距离走动校正方法
CN102707269B (zh) * 2012-06-14 2014-04-09 电子科技大学 一种机载雷达距离走动校正方法
CN102914772A (zh) * 2012-09-18 2013-02-06 西安电子科技大学 基于等效散射点的进动目标二维成像方法
CN103197285A (zh) * 2013-03-22 2013-07-10 电子科技大学 一种用于合成孔径雷达成像的导航数据拟合方法
CN103197285B (zh) * 2013-03-22 2014-08-13 电子科技大学 一种用于合成孔径雷达成像的导航数据拟合方法
CN104133215A (zh) * 2014-05-29 2014-11-05 西安电子科技大学 基于距离徙动微调和子带分割的同步轨道雷达成像方法
CN104133215B (zh) * 2014-05-29 2017-04-12 西安电子科技大学 基于距离徙动微调和子带分割的同步轨道雷达成像方法
CN104076360A (zh) * 2014-07-04 2014-10-01 西安电子科技大学 基于压缩感知的二维sar稀疏目标成像方法
CN105676220A (zh) * 2014-11-21 2016-06-15 中国航空工业集团公司雷华电子技术研究所 一种结合二维逆滤波的高分辨聚束sar自聚焦方法
CN105676220B (zh) * 2014-11-21 2019-01-04 中国航空工业集团公司雷华电子技术研究所 一种结合二维逆滤波的高分辨聚束sar自聚焦方法
CN104865571A (zh) * 2015-04-02 2015-08-26 西安电子科技大学 多通道多子带滑动聚束模式sar成像方法
CN109358329B (zh) * 2018-11-06 2020-06-30 电子科技大学 脉冲传播时间内平台机动的双基sar回波模型建立方法
CN109358329A (zh) * 2018-11-06 2019-02-19 电子科技大学 脉冲传播时间内平台机动的双基sar回波模型建立方法
CN111323761A (zh) * 2020-03-20 2020-06-23 北京华力创通科技股份有限公司 回波系统函数构建方法、装置及回波模拟器
CN111323761B (zh) * 2020-03-20 2022-04-08 北京华力创通科技股份有限公司 回波系统函数构建方法、装置及回波模拟器
CN111929681A (zh) * 2020-06-24 2020-11-13 苏州理工雷科传感技术有限公司 一种基于轻小型无人机载sar的实时成像处理系统
CN111929681B (zh) * 2020-06-24 2022-03-25 苏州理工雷科传感技术有限公司 一种基于轻小型无人机载sar的实时成像处理系统
CN114488151A (zh) * 2022-04-08 2022-05-13 中国科学院空天信息创新研究院 针对观测船只的主被动联合探测方法、装置及设备和介质
CN114488151B (zh) * 2022-04-08 2022-06-24 中国科学院空天信息创新研究院 针对观测船只的主被动联合探测方法、装置及设备和介质

Also Published As

Publication number Publication date
CN101776755B (zh) 2012-03-28

Similar Documents

Publication Publication Date Title
CN101776755B (zh) 高速平台超高分辨率sar数据预处理方法
CN102778681B (zh) 固定发射站双基地前视合成孔径雷达成像方法
CN103323841B (zh) 星机载移变式双基合成孔径雷达中的成像方法及装置
CN104898120A (zh) 一种基于回波模拟的双基前视高机动平台sar成像方法
CN103580716B (zh) 用于减少频率调制连续波高度仪系统中的自干涉效应的系统补偿方法
CN101369017B (zh) 一种移变模式双基地合成孔径雷达成像方法
CN102288961A (zh) 一种合成孔径雷达非线性调频变标的成像方法
CN104020471B (zh) 一种基于分块处理的sar实时成像方法及系统
CN108427115B (zh) 合成孔径雷达对运动目标参数的快速估计方法
CN111551922B (zh) 三维空间双/多基地雷达高速目标检测方法
CN103885061B (zh) 基于改进相位梯度自聚焦的双基sar运动补偿方法
CN102707269B (zh) 一种机载雷达距离走动校正方法
CN102590812A (zh) 基于调频连续波的sar实时成像方法
CN102749621A (zh) 一种双基地合成孔径雷达频域成像方法
CN103792535B (zh) 一种利用sar卫星测量电离层tec值的方法
CN110850384B (zh) 一种基于扫频数据产生宽带去斜回波的方法
CN103630905B (zh) 阵列天线sar极坐标交叠子孔径成像方法
CN107728117A (zh) 机载双基地sar一发两收杂波抑制方法
CN104950307A (zh) 一种机载三通道sar-gmti的精确定位方法
CN103869314A (zh) 双飞移变双基地前视合成孔径雷达成像方法
CN104777479A (zh) 基于多核dsp的前侧视sar实时成像方法
US20190324108A1 (en) Angular resolution of targets using separate radar receivers
CN101609152B (zh) 一种基于单通道sar子图像对消的运动目标检测方法
CN104330791A (zh) 一种基于频域切变的相参积累方法
CN103675815A (zh) 一种在大斜视sar成像模式下对多普勒调频率进行精确估计的方法

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
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20151026

Address after: 710075, B, building 38, Tengfei innovation center, six hi tech Zone, Xi'an hi tech Zone, Shaanxi, 606

Patentee after: XI'AN CETC XIDIAN UNIVERSITY RADAR TECHNOLOGY COLLABORATIVE INNOVATION RESEARCH INSTITUTE CO., LTD.

Address before: Xi'an City, Shaanxi province Taibai Road 710071 No. 2

Patentee before: Xidian University