CN111487614A - 基于子孔径的曲线航迹弹载sar波前重建成像方法及系统 - Google Patents
基于子孔径的曲线航迹弹载sar波前重建成像方法及系统 Download PDFInfo
- Publication number
- CN111487614A CN111487614A CN202010290861.3A CN202010290861A CN111487614A CN 111487614 A CN111487614 A CN 111487614A CN 202010290861 A CN202010290861 A CN 202010290861A CN 111487614 A CN111487614 A CN 111487614A
- Authority
- CN
- China
- Prior art keywords
- aperture
- sub
- sar
- imaging
- missile
- 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
Images
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/883—Radar or analogous systems specially adapted for specific applications for missile homing, autodirectors
-
- 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/9089—SAR having an irregular aperture
-
- 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/9041—Squint 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曲线孔径分段处理,得到若干子孔径,采用直线孔径趋于各子孔径并进行子孔径粗分辨成像;针对各子孔径数据分别进行粗分辨相干成像;将各子孔径成像图像进行旋转和相干叠加,得到全孔径聚焦雷达图像。本发明较好地处理了弹载SAR平台的各种曲线孔径,对曲线航迹SAR导引头的工程化具有重要意义;并进一步通过仿真数据表明,其能够有效处理SAR载体平台的曲线运动,以生成聚焦效果良好的SAR图像。
Description
技术领域
本发明属于雷达成像技术领域,特别涉及一种基于子孔径的曲线航迹弹载SAR波前重建成像方法及系统。
背景技术
合成孔径雷达导引头(Synthetic Aperture Radar Seeker,SARS)能够利用获得的高分辨SAR图像完成中远程导弹的惯性导航系统积累误差修正和末制导寻的,是国内外雷达导引头的一个研究热点。常规SAR要求载体平台飞行航迹匀速平直,严重限制了SAR载体平台飞行路径规划的灵活性,而飞行物在飞行过程中不可避免地存在威胁回避、地形规避、末段机动等曲线飞行状态。因此,研究曲线航迹弹载SAR成像算法具有重要的现实意义和应用价值。
SAR载体平台的有意机动可能严重背离水平匀速直线运动,使常规SAR成像算法的基本假设不符合实际。现有技术中,将多普勒带宽、脉冲响应函数等基本的SAR概念和方程推广到了SAR载体平台非水平飞行的情形下,或分别利用距离相对多普勒处理(RangeRelative Doppler Processing,RRDP)和不变映射(Invariant Mapping Technique,IMT)技术解决SAR载体平台下降导致的距离徙动和几何畸变问题;或者通过曲线SAR系统以提高雷达的测绘速度;或同时在方位向和高度向进行孔径合成的曲线SAR三维参数化成像算法;或利用SAR平台的曲线航迹解决常规SAR目标横向位移与径向速度导致的多普勒频移模糊的问题。但它们都没有针对弹载SAR平台作曲线机动时的成像问题进行研究。可以看出,目前对于曲线航迹弹载SAR成像问题的研究集中于曲线SAR系统特性和成像难点分析,所提成像算法对平台曲线运动方式和规律进行了严格限制,不符合弹载SAR平台实际运动情况。
发明内容
为此,本发明提供一种基于子孔径的曲线航迹弹载SAR波前重建成像方法及系统,较好地处理了弹载SAR平台的各种曲线孔径,对曲线航迹SAR导引头的工程化具有重要意义。
按照本发明所提供的设计方案,提供一种基于子孔径的曲线航迹弹载SAR波前重建成像方法,包含如下内容:
将弹载SAR曲线孔径分段处理,得到若干子孔径,采用直线孔径趋于各子孔径并进行子孔径粗分辨成像;
针对各子孔径数据分别进行粗分辨相干成像;
将各子孔径成像图像进行旋转和相干叠加,得到全孔径聚焦雷达图像。
作为本发明曲线航迹弹载SAR波前重建成像方法,进一步地,依据雷达载体平台航迹的曲率和每幅粗分辨SAR图像的生成时间,确定子孔径长度。
作为本发明曲线航迹弹载SAR波前重建成像方法,进一步地,依据脉冲重复频率确定子孔径的最大长度,通过现行孔径成像方法生成粗分辨快视图像,以减少子孔径内方向位压缩处理中子孔径范围内SAR平台瞬时速度的大小和方向变化。
作为本发明曲线航迹弹载SAR波前重建成像方法,进一步地,依据雷达脉冲重复间隔PRI设定子孔径长度,针对采集的孔径数据,采用并行流水方式进行子孔径成像处理来获取新的粗分辨率子孔径图像,以通过结合已获取的粗分辨率子孔径图像来得到相应精分辨SAR图像。
作为本发明曲线航迹弹载SAR波前重建成像方法,进一步地,将曲线航迹SAR的飞行轨迹近似于一系列线性子孔径加上残余运动误差,其中,将残余运动误差通过补偿沿航线方向运动误差以将其近似为空不变误差,利用线性孔径SAR波前重建算法在子孔径上进行粗分辨相干成像,生成目标区域的粗分辨图像。
作为本发明曲线航迹弹载SAR波前重建成像方法,进一步地,设定子孔径中心到场景中心的距离、斜视角,在子孔径范围内,将导弹近似为以设定速度沿直线运动,获取SAR回波信号相对快时间的傅里叶变换;忽略常数和幅度函数,利用驻留相位法获取傅里叶变换对合成孔径域的傅里叶变换,经匹配滤波后获取波数域的目标函数,以通过目标函数生成目标区域的粗分辨图像。
作为本发明曲线航迹弹载SAR波前重建成像方法,进一步地,在子孔径成像图像数据叠加之前,将波数域目标函数旋转到全孔径对应的基准坐标系下。
作为本发明曲线航迹弹载SAR波前重建成像方法,进一步地,借助二维插值进行波数域目标函数旋转。
作为本发明曲线航迹弹载SAR波前重建成像方法,进一步地,在得到多个子孔径在空间频率坐标系下二维谱后进行相干叠加,并通过逆傅里叶变换获取高分辨图像;通过递推方法来更新后续子孔径波数域的目标函数。
进一步地,本发明还提供一种基于子孔径的曲线航迹弹载SAR波前重建成像系统,包含:划分模块、粗成像模块和叠加模块,其中,
划分模块,用于将弹载SAR曲线孔径分段处理,得到若干子孔径,采用直线孔径趋于各子孔径并进行子孔径粗分辨成像;
粗成像模块,用于针对各子孔径数据分别进行粗分辨相干成像;
叠加模块,用于将各子孔径成像图像进行旋转和相干叠加,得到全孔径聚焦雷达图像。
本发明的有益效果:
本发明针对曲线航迹弹载SAR成像问题,分析了弹载平台曲线飞行情况下的SAR回波信号瞬时多普勒频率变化特点,将曲线孔径分成若干子孔径,采用直线孔径近似各子孔径并进行子孔径粗分辨成像,然后将子孔径成像结果进行旋转并相干叠加形成全孔径高分辨图像;较好地处理了弹载SAR平台的各种曲线孔径,对曲线航迹SAR导引头的工程化具有重要意义。并进一步通过仿真数据表明,其能够有效处理SAR载体平台的曲线运动,以生成聚焦效果良好的SAR图像。
附图说明:
图1为实施例中曲线航迹弹载SAR波前重建成像方法流程示意图;
图2为实施例中曲线航迹SAR成像几何图示意;
图3为实施例中子孔径与对应的成像坐标系示意;
图4为实施例中子孔径二维波数域数据示意;
图5为实施例中子孔径2线性近似示意;
图6为实施例中导弹飞行航迹与子孔径划分示意;
图7为实施例中曲线航迹SAR子孔径成像结果示意;
图8为实施例中曲线航迹SAR全孔径成像结果示意。
具体实施方式:
为使本发明的目的、技术方案和优点更加清楚、明白,下面结合附图和技术方案对本发明作进一步详细的说明。
针对曲线航迹弹载SAR成像问题,通过分析弹载平台曲线飞行情况下的SAR回波信号瞬时多普勒频率变化特点,本发明实施例,参见图1所示,提供一种基于子孔径的曲线航迹弹载SAR波前重建成像方法,包含如下内容:
S101、将弹载SAR曲线孔径分段处理,得到若干子孔径,采用直线孔径趋于各子孔径并进行子孔径粗分辨成像;
S102、针对各子孔径数据分别进行粗分辨相干成像;
S103、将各子孔径成像图像进行旋转和相干叠加,得到全孔径聚焦雷达图像。
将曲线孔径分成若干子孔径,采用直线孔径近似各子孔径并进行子孔径粗分辨成像,然后将子孔径成像结果进行旋转并相干叠加形成全孔径高分辨图像;较好地处理了弹载SAR平台的各种曲线孔径,便于生成聚焦效果良好的SAR图像,对曲线航迹SAR导引头的工程化具有重要意义。
曲线航迹SAR斜视聚束成像的几何关系如图2所示,坐标系oxy的x轴和y轴方向分别为导弹航迹中某点的法线和切线方向,成像区域中心为(xc,yc),成像区域内某点(x′n,y′n)的后向散射系数为σn。tm时刻曲线航迹SAR相位中心位置为(xr(tm),yr(tm)),相位中心的轨迹如图中虚线所示。雷达发射信号用p(τ)表示,省略发射信号窗函数并且假设雷达天线方向图是全向的,则第n个目标的回波信号可表示为:
则sn(τ,tm)相对快时间τ的Fourier变换为:
式中,P(ω)为发射信号p(τ)的Fourier变换,k=ω/c为波数。
则瞬时多普勒频率:
由于导弹运动速度快,一个Burst时间内,其运动误差不再是一个小量,运动补偿不能取得理想的效果。若将一个Burst孔径划分为若干个子孔径,在各子孔径内可近似认为SAR平台瞬时速度的大小和方向不变。由式(4)可以看出,SAR多普勒域是雷达视角域的映射,雷达视角和平台速度的连续性使得曲线航迹SAR全孔径时间内回波信号多普勒频率也是连续的,各子孔径多普勒频率是全孔径多普勒频率的一段,因此可以在各子孔径范围内按常规SAR方法对曲线航迹SAR进行近似处理生成粗分辨图像,然后在全孔径范围内将各子孔径图像相干叠加合成高分辨图像。由于子孔径对应窄的多普勒带宽,进行子孔径处理还可以降低方位向空间采样率(PRF)和运动补偿的要求,并可以实现“快视”从合成孔径工作原理看,进行子孔径处理也是很自然的。合成孔径过程本身就是由真实小孔径天线综合形成虚拟大孔径天线来提高方位分辨率的过程。将综合过程分成几个阶段,先综合部分孔径数据形成一系列中等孔径的等效天线,而后再将这些中等孔径天线进行综合。
曲线航迹弹载SAR长时间非线性运行对传统SAR运动补偿方法的可行性提出了挑战。考虑到SAR回波信号相干处理能够保留完整准确的相位信息,本发明实施例中,首先将弹载SAR曲线孔径划分为若干个子孔径,然后使用各子孔径数据进行粗分辨相干成像,最后将各子孔径图像进行旋转和相干叠加生成全孔径聚焦雷达图像。
作为本发明实施例中的曲线航迹弹载SAR波前重建成像方法,进一步地,将曲线航迹SAR的飞行轨迹近似于一系列线性子孔径加上残余运动误差,其中,将残余运动误差通过补偿沿航线方向运动误差以将其近似为空不变误差,利用线性孔径SAR波前重建算法在子孔径上进行粗分辨相干成像,生成目标区域的粗分辨图像。进一步地,设定子孔径中心到场景中心的距离、斜视角,在子孔径范围内,将导弹近似为以设定速度沿直线运动,获取SAR回波信号相对快时间的傅里叶变换;忽略常数和幅度函数,利用驻留相位法获取傅里叶变换对合成孔径域的傅里叶变换,经匹配滤波后获取波数域的目标函数,以通过目标函数生成目标区域的粗分辨图像。
如图3所示,将曲线航迹SAR的飞行轨迹用一系列线性子孔径加上残余的运动误差来近似,此时残余误差较小,可近似为空不变误差,按常规方法补偿沿航线方向的运动误差即可,在这些子孔径上按常规线性孔径SAR波前重建算法进行处理可以生成目标区域的粗分辨图像。以子孔径2为例,在子孔径中心到场景中心(xc2,yc2)的距离为Rc2,斜视角为θc2,在子孔径范围内,导弹近似以速度v2沿直线运动,记u2=v2tm,则SAR回波信号s(t,u)相对快时间的Fourier变换为:
利用驻留相位法可得s(ω,u2)对合成孔径域u2的Fourier变换为(忽略常数和幅度函数):
定义:
ky2(ω,ku2)=ku2 (8)
则
匹配滤波后波数域的目标函数为:
从子孔径波前重建算法的过程可以看出,子孔径波前重建SAR图像中保留了相位信息,可以相干叠加这些粗分辨图像以生成高分辨SAR图像。
作为本发明实施例中的曲线航迹弹载SAR波前重建成像方法,进一步地,在子孔径成像图像数据叠加之前,将波数域目标函数旋转到全孔径对应的基准坐标系下。进一步地,借助二维插值进行波数域目标函数旋转。
子孔径2的成像坐标系o′x′2y′2相对全孔径的成像坐标系o′x′y′有一定角度的旋转,大小为
其中k2为子孔径2的斜率。因此,在子孔径数据叠加之前,需要将式(10)表示的波数域目标函数旋转到全孔径对应的基准坐标系下,如图4所示。旋转可以借助二维插值来实现,旋转后波数域目标函数为
对F2(kx,ky)分别进行两个方向上的逆Fourier变换就可以得到坐标系o′x′y′下的目标区域粗分辨图像。
作为本发明实施例中的曲线航迹弹载SAR波前重建成像方法,进一步地,在得到多个子孔径在空间频率坐标系下二维谱后进行相干叠加,并通过逆傅里叶变换获取高分辨图像;通过递推方法来更新后续子孔径波数域的目标函数。
若利用Ms个子孔径的数据生成高分辨SAR图像,那么在获得了第Ns个子孔径的波数域SAR数据后,波数域目标函数可以写为:
得到Ms个子孔径的在空间频率坐标系下二维谱后,将它们相干叠加并作逆Fourier变换就可以得到坐标系o′x′y′下的高分辨SAR图像。
前一幅高分辨SAR图像的波数域目标函数可以表示为:
因此,可以用递推的方法更新目标函数:
作为本发明实施例中的曲线航迹弹载SAR波前重建成像方法,进一步地,依据雷达载体平台航迹的曲率和每幅粗分辨SAR图像的生成时间,确定子孔径长度。
子孔径的划分是基于子孔径的曲线航迹SAR成像的一个重要步骤。与常规线性孔径SAR不同,曲线航迹SAR用来近似每一段子孔径的直线斜率是变化的,子孔径的长度可以是固定的也可以是变化的,这主要取决于雷达载体平台航迹的曲率和每幅粗分辨SAR图像的生成时间,关键是保证在一个子孔径间隔内,雷达平台没有剧烈的航迹跳变,经过运动补偿后子孔径内的方位向压缩可以常规线性孔径SAR的方式进行。子孔径划分的最大长度需要同时满足下面三个限制原则。
作为本发明实施例中的曲线航迹弹载SAR波前重建成像方法,进一步地,依据脉冲重复频率确定子孔径的最大长度,通过现行孔径成像方法生成粗分辨快视图像,以减少子孔径内方向位压缩处理中子孔径范围内SAR平台瞬时速度的大小和方向变化。
若子孔径多普勒带宽为BK,为避免方位频谱的混叠,需使PRF>BK。另外PRF与测绘带宽、脉冲宽度之间存在基本的限制关系:
其中,Wr是测绘带距离宽度,θ为天线视角,Tp是发射脉冲宽度,Tsp是保护间隔时间,通常允许随信号距离延迟时间的变化而调整。
从式(16)可知,测绘带距离宽度等系统参数决定了PRF的上限,同时这也决定了子孔径多普勒带宽的最大值。由多普勒频率与雷达视角的关系,雷达视角的最大变化范围也随之确定,即方位向不模糊的最大子孔径长度可以确定。
将曲线SAR子孔径内的方位向压缩按常规线性孔径SAR的方式进行处理,就要保证子孔径范围内SAR平台的瞬时速度的大小和方向近似不变,即在现有运动补偿精度条件下用常规线性孔径SAR成像的方法可以生成满足要求的粗分辨“快视”图像。仍以图3中子孔径2为例进行讨论,将子孔径2用直线近似,如图5。
上式可写为:
式中,为用直线近似子孔径2的残余距离误差,为在直线近似孔径上的SAR平台速度。则式(18)最后一项为用直线近似子孔径引起的相位误差,假设选用窄波束宽度运动补偿,补偿函数为其中为目标区域的平均斜视角,为子孔径2成像坐标系下目标区域中心。为保证成像质量,经过补偿后应使直线近似子孔径引起的相位误差足够小,假若使其小于π/2,即
由于
作为本发明实施例中的曲线航迹弹载SAR波前重建成像方法,进一步地,依据雷达脉冲重复间隔PRI设定子孔径长度,针对采集的孔径数据,采用并行流水方式进行子孔径成像处理来获取新的粗分辨率子孔径图像,以通过结合已获取的粗分辨率子孔径图像来得到相应精分辨SAR图像。
“快视”成像要求提供足够高的数据率供导引头控制使用,而对图像分辨率要求不高。子孔径划分要充分考虑“快视”的要求,在满足“快视”对分辨率要求的基础上,以足够高的速度生成粗分辨图像,这就限制了子孔径的最大长度。由子孔径数据成像处理过程可知生成一幅参考坐标系下的子孔径“快视”图像主要的运算有两次FFT、残余误差补偿和匹配滤波所需的两次多点乘法操作、一次二维插值以及两次IFFT。假设子孔径测量回波信号快慢时间采样点数分别为N和M,“快视”SAR图像距离向和方位向采样点数分别为k1和k2,则两次FFT和两次IFFT的计算量为:
复加:2×MNlog2MN+2×k1k2log2k1k2
两次多点乘法操作的计算量为:
复乘:2×MN
二维插值有多种方法,计算量不一。假设将二维插值分成两次一维插值进行,插值滤波器长度为lp,不考虑插值滤波后重采样位置的计算,则一次二维插值计算量为:
复乘:lp×k1+lp×k2
复加:lp×k1+lp×k2-2
则生成子孔径“快视”SAR图像的总计算量为:
复乘:MNlog2MN+k1k2log2k1k2+2MN+lp(k1+k2)
复加:2MNlog2MN+2k1k2log2k1k+lp(k1+k2)-2
子孔径成像过程相对独立,只要陆续采集的孔径数据可以构成一个新的子孔径,就可以对其进行成像处理,得到新的粗分辨子孔径图像,不必等整个孔径的数据收集完后再进行处理。结合之前已有的一组粗分辨图像,就可以得到相应的精分辨SAR图像。可以看出,子孔径成像处理流程便于并行和流水处理,这会大大加快成像处理速度。当然子孔径也不宜过短。过短不仅会使子孔径粗分辨图像分辨率过低,而且子孔径个数的增多也会使生成一幅高分辨SAR图像的计算量过大,必须对子孔径的长度折衷考虑,一般子孔径长度约为一百个PRI左右。
进一步地,基于上述的方法,本发明实施例还提供一种基于子孔径的曲线航迹弹载SAR波前重建成像系统,包含:划分模块、粗成像模块和叠加模块,其中,
划分模块,用于将弹载SAR曲线孔径分段处理,得到若干子孔径,采用直线孔径趋于各子孔径并进行子孔径粗分辨成像;
粗成像模块,用于针对各子孔径数据分别进行粗分辨相干成像;
叠加模块,用于将各子孔径成像图像进行旋转和相干叠加,得到全孔径聚焦雷达图像。
为验证本发明技术方案有效性,下面结合仿真数据对点目标进行成像仿真:
仿真参数如下:发射信号中心频率fc=10GHz,调频带宽B=100MHz,时宽Tp=4μs,孔径中心斜视角θ和斜距Rc分别为60°和20km,导弹在水平面内进行加速度为a=2g的转弯,g为重力加速度,y方向速度为1000m/s,且在转弯过程中保持不变,其飞行轨迹如图6所示,y方向合成孔径长度为260m,分为四个子孔径进行处理,由于子孔径波数域带宽较小,脉冲重复频率PRF为1000Hz。
表1基于子孔径的曲线航迹SAR波前重建算法成像质量指标
图7是第四个子孔径的粗分辨成像结果,图8是全孔径精分辨成像结果,成像过程没有进行加权处理,全孔径波数域带宽是四个子孔径波数域带宽的和,两图的分辨率差异很明显。需要注意的是,子孔径处理过程中,为了保证子孔径的波数域频域分辨率与全孔径相同,子孔径成像处理需要包括一个补零的过程,将子孔径长度补零至与全孔径相同。另外由于各子孔径对同一目标区域进行成像时的斜视角不同,得到的各子孔径图像之间并不具有对应关系,所以各子孔径的波数域数据要旋转变换的统一的坐标系下才能进行全孔径的相干叠加。表1对仿真的结果进行了基于点目标的SAR图像质量评估,其中δa,δr,Pa,Pr,Ia,Ir分别表示方位向分辨率、距离向分辨率、方位向峰值旁瓣比、距离向峰值旁瓣比、方位向积分旁瓣比、距离向积分旁瓣比。由表中各指标可见基于子孔径的曲线航迹SAR波前重建算法较好地处理了平台的曲线运动。
基于以上内容,本发明实施例中,在分析曲线航迹SAR回波信号特点的基础上,通过将曲线孔径分段处理,并将各子孔径的粗分辨成像结果进行旋转和相干叠加,最终形成了全孔径高分辨SAR成像。曲线航迹SAR技术对载体平台飞行轨迹没有苛刻要求,大大增加了弹道规划的灵活性,使得SAR导引头在突防机动、地形规避时仍可对感兴趣区域进行成像,优势突出,相信随着其技术的成熟将获得广泛的应用。
除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对步骤、数字表达式和数值并不限制本发明的范围。
基于上述的方法,本发明实施例还提供一种服务器,包括:一个或多个处理器;存储装置,用于存储一个或多个程序,当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现上述的方法。
基于上述的方法,本发明实施例还提供一种计算机可读介质,其上存储有计算机程序,其中,该程序被处理器执行时实现上述的方法。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种基于子孔径的曲线航迹弹载SAR波前重建成像方法,其特征在于,包含如下内容:
将弹载SAR曲线孔径分段处理,得到若干子孔径,采用直线孔径趋于各子孔径并进行子孔径粗分辨成像;
针对各子孔径数据分别进行粗分辨相干成像;
将各子孔径成像图像进行旋转和相干叠加,得到全孔径聚焦雷达图像。
2.根据权利要求1所述的基于子孔径的曲线航迹弹载SAR波前重建成像方法,其特征在于,依据雷达载体平台航迹的曲率和每幅粗分辨SAR图像的生成时间,确定子孔径长度。
3.根据权利要求1或2所述的基于子孔径的曲线航迹弹载SAR波前重建成像方法,其特征在于,依据脉冲重复频率确定子孔径的最大长度,通过现行孔径成像方法生成粗分辨快视图像,以减少子孔径内方向位压缩处理中子孔径范围内SAR平台瞬时速度的大小和方向变化。
4.根据权利要求1所述的基于子孔径的曲线航迹弹载SAR波前重建成像方法,其特征在于,依据雷达脉冲重复间隔PRI设定子孔径长度,针对采集的孔径数据,采用并行流水方式进行子孔径成像处理来获取新的粗分辨率子孔径图像,以通过结合已获取的粗分辨率子孔径图像来得到相应精分辨SAR图像。
5.根据权利要求1所述的基于子孔径的曲线航迹弹载SAR波前重建成像方法,其特征在于,将曲线航迹SAR的飞行轨迹近似于一系列线性子孔径加上残余运动误差,其中,将残余运动误差通过补偿沿航线方向运动误差以将其近似为空不变误差,利用线性孔径SAR波前重建算法在子孔径上进行粗分辨相干成像,生成目标区域的粗分辨图像。
6.根据权利要求1或5所述的基于子孔径的曲线航迹弹载SAR波前重建成像方法,其特征在于,设定子孔径中心到场景中心的距离、斜视角,在子孔径范围内,将导弹近似为以设定速度沿直线运动,获取SAR回波信号相对快时间的傅里叶变换;忽略常数和幅度函数,利用驻留相位法获取傅里叶变换对合成孔径域的傅里叶变换,经匹配滤波后获取波数域的目标函数,以通过目标函数生成目标区域的粗分辨图像。
7.根据权利要求6所述的基于子孔径的曲线航迹弹载SAR波前重建成像方法,其特征在于,在子孔径成像图像数据叠加之前,将波数域目标函数旋转到全孔径对应的基准坐标系下。
8.根据权利要求7所述的基于子孔径的曲线航迹弹载SAR波前重建成像方法,其特征在于,借助二维插值进行波数域目标函数旋转。
9.根据权利要求6所述的基于子孔径的曲线航迹弹载SAR波前重建成像方法,其特征在于,在得到多个子孔径在空间频率坐标系下二维谱后进行相干叠加,并通过逆傅里叶变换获取高分辨图像;通过递推方法来更新后续子孔径波数域的目标函数。
10.一种基于子孔径的曲线航迹弹载SAR波前重建成像系统,其特征在于,包含:划分模块、粗成像模块和叠加模块,其中,
划分模块,用于将弹载SAR曲线孔径分段处理,得到若干子孔径,采用直线孔径趋于各子孔径并进行子孔径粗分辨成像;
粗成像模块,用于针对各子孔径数据分别进行粗分辨相干成像;
叠加模块,用于将各子孔径成像图像进行旋转和相干叠加,得到全孔径聚焦雷达图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010290861.3A CN111487614B (zh) | 2020-04-14 | 2020-04-14 | 基于子孔径的曲线航迹弹载sar波前重建成像方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010290861.3A CN111487614B (zh) | 2020-04-14 | 2020-04-14 | 基于子孔径的曲线航迹弹载sar波前重建成像方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111487614A true CN111487614A (zh) | 2020-08-04 |
CN111487614B CN111487614B (zh) | 2022-07-12 |
Family
ID=71812758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010290861.3A Active CN111487614B (zh) | 2020-04-14 | 2020-04-14 | 基于子孔径的曲线航迹弹载sar波前重建成像方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111487614B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112578380A (zh) * | 2020-11-05 | 2021-03-30 | 北京理工大学 | 基于载机航迹偏差估计的无人机载sar自聚焦方法 |
CN112946649A (zh) * | 2021-04-08 | 2021-06-11 | 电子科技大学 | 一种适用于任意子孔径长度的pfa成像方法 |
CN113050087A (zh) * | 2021-03-17 | 2021-06-29 | 电子科技大学 | 一种用于曲线合成孔径雷达的子孔径划分方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040201514A1 (en) * | 2003-01-14 | 2004-10-14 | The Regents Of The University Of California | Differential optical synthetic aperture radar |
CN103901428A (zh) * | 2014-03-28 | 2014-07-02 | 西安电子科技大学 | 弹载sar子孔径前斜视高阶非线性调频变标成像方法 |
CN104391297A (zh) * | 2014-11-17 | 2015-03-04 | 南京航空航天大学 | 一种划分子孔径pfa雷达成像方法 |
CN106054187A (zh) * | 2016-05-23 | 2016-10-26 | 西安电子科技大学 | 基于斜距模型下的大斜视sar曲线轨迹波数域成像方法 |
CN109959932A (zh) * | 2019-04-08 | 2019-07-02 | 西安电子科技大学 | 基于下降段曲线轨迹的雷达前视三维成像方法 |
CN110244300A (zh) * | 2019-06-24 | 2019-09-17 | 杭州电子科技大学 | 基于球体模型和fenlcs算法的弹载sar平飞段高分辨率成像方法 |
-
2020
- 2020-04-14 CN CN202010290861.3A patent/CN111487614B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040201514A1 (en) * | 2003-01-14 | 2004-10-14 | The Regents Of The University Of California | Differential optical synthetic aperture radar |
CN103901428A (zh) * | 2014-03-28 | 2014-07-02 | 西安电子科技大学 | 弹载sar子孔径前斜视高阶非线性调频变标成像方法 |
CN104391297A (zh) * | 2014-11-17 | 2015-03-04 | 南京航空航天大学 | 一种划分子孔径pfa雷达成像方法 |
CN106054187A (zh) * | 2016-05-23 | 2016-10-26 | 西安电子科技大学 | 基于斜距模型下的大斜视sar曲线轨迹波数域成像方法 |
CN109959932A (zh) * | 2019-04-08 | 2019-07-02 | 西安电子科技大学 | 基于下降段曲线轨迹的雷达前视三维成像方法 |
CN110244300A (zh) * | 2019-06-24 | 2019-09-17 | 杭州电子科技大学 | 基于球体模型和fenlcs算法的弹载sar平飞段高分辨率成像方法 |
Non-Patent Citations (3)
Title |
---|
SI CHEN等: "A New Imaging Algorithm for Forward-Looking Missile-Borne Bistatic SAR", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》 * |
周鹏: "弹载SAR多种工作模式的成像算法研究", 《中国优秀博硕士学位论文全文数据库(博士)信息科技辑》 * |
江淮等: "一种弹载SAR子孔径成像算法", 《电子与信息学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112578380A (zh) * | 2020-11-05 | 2021-03-30 | 北京理工大学 | 基于载机航迹偏差估计的无人机载sar自聚焦方法 |
CN112578380B (zh) * | 2020-11-05 | 2023-09-19 | 北京理工大学 | 基于载机航迹偏差估计的无人机载sar自聚焦方法 |
CN113050087A (zh) * | 2021-03-17 | 2021-06-29 | 电子科技大学 | 一种用于曲线合成孔径雷达的子孔径划分方法 |
CN112946649A (zh) * | 2021-04-08 | 2021-06-11 | 电子科技大学 | 一种适用于任意子孔径长度的pfa成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111487614B (zh) | 2022-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111443339B (zh) | 一种双基sar空变校正成像方法、装置、设备及存储介质 | |
CN111487614B (zh) | 基于子孔径的曲线航迹弹载sar波前重建成像方法及系统 | |
CN107229048B (zh) | 一种高分宽幅sar动目标速度估计与成像方法 | |
US7277042B1 (en) | Compensation of flight path deviation for spotlight SAR | |
USH1720H (en) | Time frequency processor for radar imaging of moving targets | |
Guo et al. | Extended wavenumber domain algorithm for highly squinted sliding spotlight SAR data processing | |
CN106842210A (zh) | 一种新的多子阵合成孔径声纳快速成像算法 | |
CN114545411B (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
Tang et al. | Processing of monostatic SAR data with general configurations | |
Yang et al. | A novel three-step image formation scheme for unified focusing on spaceborne SAR data | |
CN113589285A (zh) | 一种飞行器sar实时成像方法 | |
CN111551934A (zh) | 一种用于无人机载sar成像的运动补偿自聚焦方法与装置 | |
CN110244300B (zh) | 基于球体模型和fenlcs算法的弹载sar平飞段高分辨率成像方法 | |
Kim et al. | Doppler shifting technique for generating multi-frames of video SAR via sub-aperture signal processing | |
Liu et al. | Ambiguities Suppression for Azimuth Multichannel SAR Based on ${L_ {2, q}} $ Regularization With Application to Gaofen-3 Ultra-Fine Stripmap Mode | |
Song et al. | Focusing high-maneuverability bistatic forward-looking SAR using extended azimuth nonlinear chirp scaling algorithm | |
CN109143235B (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
CN113759372A (zh) | 弹载大斜视小孔径多通道sar的成像方法 | |
CN116500626B (zh) | 一种多接收阵元数据收发合置转换方法 | |
CN103217677A (zh) | 一种基于联合检测量的单通道sar动目标检测方法 | |
CN111880179A (zh) | 一种弹载弧线俯冲大斜视tops sar的成像方法 | |
CN116719027A (zh) | 一种双向滑动聚束模式下的星机双基sar成像方法 | |
CN114002666B (zh) | 任意天线构型下星载ati-sar洋流流速提取方法及设备 | |
CN113176570B (zh) | 一种斜视sar时域成像自聚焦方法 | |
CN113406624B (zh) | 一种高分辨率星载sar高效时频混合成像方法及系统 |
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 |