CN108205135B - 基于无插值融合快速后向投影的雷达视频成像方法 - Google Patents

基于无插值融合快速后向投影的雷达视频成像方法 Download PDF

Info

Publication number
CN108205135B
CN108205135B CN201810058866.6A CN201810058866A CN108205135B CN 108205135 B CN108205135 B CN 108205135B CN 201810058866 A CN201810058866 A CN 201810058866A CN 108205135 B CN108205135 B CN 108205135B
Authority
CN
China
Prior art keywords
sub
image
radar
aperture
images
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.)
Active
Application number
CN201810058866.6A
Other languages
English (en)
Other versions
CN108205135A (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.)
Xidian University
Xian Cetc Xidian University Radar Technology Collaborative Innovation Research Institute Co Ltd
Original Assignee
Xidian University
Xian Cetc Xidian University Radar Technology Collaborative Innovation Research Institute Co Ltd
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, Xian Cetc Xidian University Radar Technology Collaborative Innovation Research Institute Co Ltd filed Critical Xidian University
Priority to CN201810058866.6A priority Critical patent/CN108205135B/zh
Publication of CN108205135A publication Critical patent/CN108205135A/zh
Application granted granted Critical
Publication of CN108205135B publication Critical patent/CN108205135B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004SAR image acquisition techniques
    • G01S13/9017SAR image acquisition techniques with time domain processing of the SAR signals in azimuth
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004SAR image acquisition techniques
    • G01S13/9005SAR image acquisition techniques with optical processing of the SAR signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明属于合成孔径成像技术领域,公开了一种基于无插值融合快速后向投影的雷达视频成像方法,FFBP算法首先对数据划分子孔径,然后每个子孔径数据形成一幅以子孔径中心为原点的局部直角坐标子孔径图像。子孔径图像只是利用了部分孔径数据,因此具有方位向低分辨的特点。然后对形成的子孔径图像进行两两融合,得到新的较高分辨率的子图像,该过程以迭代的方式进行直到得到一幅全分辨的极坐标图像,最后将极坐标图像插值到直角坐标网格得到最后的图像。

Description

基于无插值融合快速后向投影的雷达视频成像方法
技术领域
本发明属于合成孔径成像技术领域,尤其涉及一种基于无插值融合快速后向投影的雷达视频成像方法,是一种利用合成孔径技术进行快速视频成像的方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)作为一种微波有源成像系统,具有全天时、全天候工作的特点。它根据合成孔径原理,通过对载机在不同方位位置形成“合成阵列”录取的回波信号进行相干积累,从而获取方位向的高分辨。由于SAR成像不受气象条件限制,而且不分白昼与黑夜,均能够提供高分辨的对地观测图像,因此在军事、科研、国民经济许多领域有着广阔的应用前景。
当前,各国均在探索SAR成像新的发展道路。视频合成孔径雷达(VideoSAR,ViSAR)作为一种新体制雷达,能够在一定的分辨率和高帧率的条件下产生一系列高数据率的连续图像,并对运动目标具有阴影效应,有利于实现动目标的有效检测跟踪。ViSAR在提供高清高帧率图像的同时能够降低尺寸、质量和功率(SWAP),是SAR未来发展的一个重要方向,受到了国内外的高度重视。但是视频SAR的帧率高,数据率大,对系统硬件要求比较高,一定程度上限制了视频SAR的发展与应用。
发明内容
针对上述问题,本发明的目的在于提供一种基于无插值融合快速后向投影的雷达视频成像方法,该方法利用视频数据具有高重叠率的特点,能够进行快速成像,从而提高视频帧率,降低对系统硬件的要求,具有成像时间短、帧率高的特点,能够提高效率,降低成本。
为达到上述目的,本发明采用如下技术方案予以实现。
一种基于无插值融合快速后向投影的雷达视频成像方法,所述方法包括如下步骤:
步骤1,设定雷达视线方向为距离向,雷达平台运动方向为方位向,且雷达回波采样时间为快时间,雷达平台运动时间为慢时间;所述雷达平台采用“一走一停”机制,所述雷达在每个慢时间发射信号并接收对应的回波信号;
步骤2,设定雷达成像的方位分辨率为ρa,雷达平台运动速度为v,雷达平台到成像场景的参考距离为Rs,从而确定雷达成像所需运动的距离L;
步骤3,雷达在每个慢时间对回波信号进行采样,得到大小为N×M的回波数据块,其中,N为所述回波数据块的方位向采样点数,且N=L×Fa/v,M为所述回波数据块的距离向采样点数,Fa为雷达发射信号的脉冲重复频率;所述大小为N×M的回波数据块形成一幅方位分辨率为ρa的SAR图像;
步骤4,对所述N×M的回波数据块进行脉冲压缩,得到脉冲压缩后的回波数据;
步骤5,将所述脉冲压缩后的回波信号划分为K(0)个子孔径数据块,且K(0)为2的整数次幂,每个子孔径数据块方位向采样点数为Nsub=N/K(0)
步骤6,对K(0)个子孔径数据块进行成像,得到K(0)个子图像,并对所述K(0)个子图像进行子图像融合,则进行G=log2(K(0))次子图像融合后得到形成一幅方位分辨率为ρa的SAR图像;
步骤7,设雷达新接收到大小为Nnew×M的回波数据块,且Nnew<N,并对所述大小为Nnew×M的回波数据块进行脉冲压缩,得到脉冲压缩后的新的回波数据;
步骤8,将所述脉冲压缩后的新的回波数据划分为方位向采样点数为Nsub的T个子孔径数据块,且T=floor[Nnew/Nsub],表示向下取整操作;
步骤9,对所述T个子孔径数据块进行成像,得到T子图像,获取前一幅SAR图像K(0)个子图像中的后K(0)-T个子图像;从而将所述T个子图像和K(0)-T个子图像进行子图像融合,得到一幅方位分辨率为ρa的新的SAR图像;
步骤10,重复执行步骤7至步骤9,得到多幅SAR图像,所述多幅SAR图像形成视频成像。
与现有技术相比,本发明的有益效果是:(1)本发明采用重叠子孔径的时域成像方法,充分发挥了时域成像的优势,具有更高的成像效率;(2)采用图像域平移的方法实现图像间的融合,计算效率更高;(3)本发明减少了插值引入的误差,图像质量有稍微提高。第四,本发明能够降低对系统硬件的要求,降低成本。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明对回波数据的分块示意图;
图2为本发明中单幅图像流程图;
图3为本发明所采用成像算法中的子图像融合前后坐标关系示意图;
图4为本发明视频成像的流程图;
图5为本发明中子图像融合示意图;
图6为本发明的仿真结果。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
由于视频SAR帧率高,因此相邻两幅图像的数据有很高的重叠率。一般来说,帧率越高,重叠率越高。视频成像时,首先对数据按照方位分辨率的要求进行数据块划分,如图1所示。
参照图2为本发明的基于无插值融合的FFBP(快速后向投影)快速成像方法的流程图;所述基于无插值融合的FFBP快速成像算法的视频成像方法,包括以下说明:
本发明实施例采用合成孔径雷达成像技术。合成孔径雷达相比阵列来说能够显著降低成本,并能获得同样的高分辨图像。采用合成孔径雷二维高分辨力,距离向依靠发射大时宽-带宽积的线性调频信号,对回波信号进行脉冲压缩得到距离向的高分辨;方位向根据合成孔径原理,通过对载机在不同方位位置形成“合成阵列”录取的回波信号进行相干积累,从而获取方位向的高分辨。
本发明实施例采用快速分解后向投影算法(FFBP)。FFBP为一种快速时域成像方法。由于视频成像图像间具有很高的数据重叠率,而时域成像方法能够对数据逐脉冲处理,容易进行运动补偿,而频域成像方法是对全孔径数据处理,因此视频成像更适合采用时域成像方法,考虑算法的运算量,本发明采用和频域算法计算量可比拟的快速时域成像算法。
所述发明采用快速分解后向投影算法(FFBP),子图像建立在以子孔径中心为原点的局部直角坐标系,由于子图像只利用了部分数据,因此方位向分辨率低,即方位向网格稀疏。如图3所示,X1O1Y1为融合前的坐标轴之一,XAOAYA为融合后的坐标轴,可以看出子图像融合前后坐标系存在平移的关系,而图像域平移能够在频域采用相位相乘的方法快速精确的实现,因此利用该性质能够降低算法的计算量,同时降低插值引入的误差。
本发明实施例采用视频成像技术,形成视频的流程如图4所示。视频能够连续多角度成像,通过持续影像的渐变显现形成视觉活动图像,能克服常规SAR系统帧率低、动目标检测跟踪难等缺陷,提高识别准确率,而且视频成像在提供高清高帧图像的同时能够降低尺寸、质量和功率(SWAP),有利于系统小型化。
以下对本发明方案进行具体说明:
本发明实施例提供一种基于无插值融合快速后向投影的雷达视频成像方法,所述方法包括如下步骤:
步骤1,设定雷达视线方向为距离向,雷达平台运动方向为方位向,且雷达回波采样时间为快时间,雷达平台运动时间为慢时间;所述雷达平台采用“一走一停”机制,所述雷达在每个慢时间发射信号并接收对应的回波信号。
步骤2,设定雷达成像的方位分辨率为ρa,雷达平台运动速度为v,雷达平台到成像场景的参考距离为Rs,从而确定雷达成像所需运动的距离L。
步骤2中,确定雷达成像所需运动的距离L即合成孔径
Figure BDA0001554685340000061
其中,λ为雷达发射信号的波长。
步骤3,雷达在每个慢时间对回波信号进行采样,得到大小为N×M的回波数据块,其中,N为所述回波数据块的方位向采样点数,且N=L×Fa/v,M为所述回波数据块的距离向采样点数,Fa为雷达发射信号的脉冲重复频率;所述大小为N×M的回波数据块形成一幅方位分辨率为ρa的SAR图像。
步骤3中,雷达发射信号的脉冲重复频率Fa=ka×2v/Da,其中,ka为过采样率,Da为雷达方位向孔径长度,从而雷达在慢时间停留的位置x′=n*v/Fa,n=1,2,3,……。
步骤4,对所述N×M的回波数据块进行脉冲压缩,得到脉冲压缩后的回波数据。
假设发射的是线性调频波
Figure BDA0001554685340000062
其中,wr(·)为发射信号包络,
Figure BDA0001554685340000063
为快时间,fc为中心频率,TP为脉宽,γ为调频斜率。点目标回波经过解调后可以表示为
Figure BDA0001554685340000064
式中,tm为慢时间,R(tm)为雷达在tm时刻天线相位中心(antenna phase center,APC)到目标的距离。
首先,对回波数据块进行脉冲压缩。假设信号包络为矩形,接收信号频谱为
Figure BDA0001554685340000065
匹配滤波器为
Figure BDA0001554685340000071
匹配滤波器与Sr(f)相乘后的频谱为
Figure BDA0001554685340000072
压缩信号是上式的傅里叶变换
Figure BDA0001554685340000073
压缩后的信号,在距离向上是在
Figure BDA0001554685340000074
处出现峰值的sinc函数,因此目标被压缩到与其距离相对应的位置。
步骤5,将所述脉冲压缩后的回波信号划分为K(0)个子孔径数据块,且K(0)为2的整数次幂,每个子孔径数据块方位向采样点数为Nsub=N/K(0)
步骤6,对K(0)个子孔径数据块进行成像,得到K(0)个子图像,并对所述K(0)个子图像进行子图像融合,则进行G=log2(K(0))次子图像融合后得到形成一幅方位分辨率为ρa的SAR图像。
步骤6具体包括如下子步骤:
(6a)每幅子图像建立在以其子孔径中心为与原点的局部直角坐标系中,对于第k个子孔径数据块,其子孔径中心设为
Figure BDA0001554685340000075
以所述子孔径中心为原点建立局部直角坐标系
Figure BDA0001554685340000076
则第k个子图像
Figure BDA0001554685340000077
表示为:
Figure BDA0001554685340000078
其中,k=1,2,…,K(0),x为子孔径数据块中采样点的方位向位置,且-l0/2≤x<l0/2,l0表示第k个子孔径的长度,且l(0)=L/K(0)
Figure BDA0001554685340000081
表示在位置
Figure BDA0001554685340000082
处的脉冲压缩后的回波信号,λ为雷达发射信号的波长,
Figure BDA0001554685340000083
表示雷达在位置
Figure BDA0001554685340000084
时天线相位中心到位置局部直角坐标系
Figure BDA0001554685340000085
的瞬时斜距,且
Figure BDA0001554685340000086
(6b)令k=1,2,…,K(0),从而得到K(0)个子图像;
(6c)对所述第1个子图像至第K(0)个子图像两两一组进行融合,则共需要融合次数为G=log2(K(0))次,且每次融合后的子图像数目变为融合前的1/2;
则第g-1次子图像融合后得到的子图像为
Figure BDA0001554685340000087
其中,表示
Figure BDA0001554685340000088
相干融合,
Figure BDA0001554685340000089
表示第g-1次子图像融合的第2k-1个子图像,
Figure BDA00015546853400000810
表示第g-1次子图像融合的第2k个子图像;
如图3所示,子图像建立在局部直角坐标系中,融合前后坐标系存在平移的关系,而图像域平移能够在频域采用相位相乘的方法快速精确的实现,因此利用该性质能够降低算法的计算量,同时降低插值引入的误差。
在子步骤(6c)中对子图像两两一组进行融合具体包括:
(a)融合前的一幅子图像的子孔径中心为O1,子孔径长度为l(g),融合后的子图像的子孔径中心为OA,子孔径长度为l(g+1)
(b)设点目标P在融合前的一个坐标系中的坐标为(x1,y1),在融合后的坐标系中的坐标为(xA,yA),则点目标P融合前到融合后的平移量(△x,△y)为:
Figure BDA0001554685340000091
(c)由于O1OA=(l(g+1)-l(g))/2=l(g+1)/4,则
Figure BDA0001554685340000092
其中,l(g+1)为融合后的子孔径长度,θ为斜视角,且l(g)=2(g)l0,l0为初始子孔径长度l0=L/K(0),L为合成孔径长度;
(d)融合前两个子图像的子孔径中心关于融合后的子图像的子孔径中心对称,从而确定融合前两个子图像的平移量分别为(△x,△y)和(-△x,-△y)。
在频域实现两个子图像平移的过程为:
采用FFT将两个子图像分别变换到二维频域,其中一个子图像在二维频域乘以第一线性相位,另一个子图像在二维频域乘以第二线性相位:
所述第一线性相位的表达式为:
Figure BDA0001554685340000093
所述第二线性相位的表达式为:
Figure BDA0001554685340000094
其中,S1x和S1y分别为一个子图像在方位频域和距离频域相乘的线性相位,S2x和S2y分别为另一个子图像在方位频域和距离频域相乘的线性相位,fa为方位频率,fr为距离频率。
(6d)令g=1,2,...G,进行G次子图像融合后得到一幅方位分辨率为ρa的SAR图像。
因为斜视的关系,融合得到的最终直角坐标图像也是斜视的,因此需要通过插值变换到正侧视的全分辨率图像。
步骤7,设雷达新接收到大小为Nnew×M的回波数据块,且Nnew<N,并对所述大小为Nnew×M的回波数据块进行脉冲压缩,得到脉冲压缩后的新的回波数据。
步骤8,将所述脉冲压缩后的新的回波数据划分为方位向采样点数为Nsub的T个子孔径数据块,且T=floor[Nnew/Nsub],表示向下取整操作。
步骤9,对所述T个子孔径数据块进行成像,得到T子图像,获取前一幅SAR图像K(0)个子图像中的后K(0)-T个子图像;从而将所述T个子图像和K(0)-T个子图像进行子图像融合,得到一幅方位分辨率为ρa的新的SAR图像。
需要注意的是,新生成的图像方位分辨率低,因此融合时需要进行方位升采样。将融合后的图像插值的到所需要的图像。同理,采用同样的方法可以得到一系列的高分辨SAR图像,然后利用多媒体融合技术可以生成连续的视频。
步骤10,重复执行步骤7至步骤9,得到多幅SAR图像,所述多幅SAR图像形成视频成像。
需要补充的是,设相邻两幅SAR图像的数据重叠率为η,雷达新接收的回波数据块方位向点数Nnew=N(1-η),则多幅SAR图像形成的视频成像的最大帧率Fmax为:
Fmax=Fa/Nnew=v/(L(1-η))
其中,Fa表示脉冲重复频率,v表示雷达平台运动速度,L表示雷达成像所需运动的距离。
由此可以看出,采用本方法所能达到的最大帧率与数据重叠率有关的,重叠率越高,帧率越高,一般η可以达到90%以上。通常限于系统硬件水平,视频SAR系统的视频帧率达不到理论分析的最大值,本方法能够通过进一步降低算法的计算量来提高视频SAR的帧率,提高了视频SAR的效率。
基于无插值融合的FFBP成像算法的视频成像方法,包含:利用时域成像方法(FFBP算法)进行成像,并在成像时采用图像域平移的方法进行子图像融合。
所述数据成像方法基于快速分解后向投影算法(FFBP)。图5以划分4个子孔径为例进行了示意说明,FFBP算法首先对数据划分子孔径,然后每个子孔径数据形成一幅以子孔径中心为原点的局部直角坐标子孔径图像。子孔径图像只是利用了部分孔径数据,因此具有方位向低分辨的特点。然后对形成的子孔径图像进行两两融合,得到新的较高分辨率的子图像,该过程以迭代的方式进行直到得到一幅全分辨的极坐标图像,最后将极坐标图像插值到直角坐标网格得到最后的图像。
传统方法在子图像融合时采用两维插值的方法,不仅计算量大,还会引入插值误差,经过多次迭代融合误差会积累,最终会影响图像质量。本发明采用直角坐标形成子图像,同样具有方位向低分辨的特点。而在子图像融合时采用图像域平移的方法,能够同时降低了计算量和引入的误差。
所述发明形成新的一帧图像时,利用同样的方法对新接收到的数据进行处理,然后和图像间重叠部分进行融合。对融合后的图像进行插值得到最终的直角坐标图像。用同样的方法形成一系列的图像,经过多媒体融合得到成像视频。图6所示为视频成像仿真结果,视频帧率为10帧/s,图6(a)为截取的视频第20帧的图像,视频中的动目标用矩形框进行了标记,并在图6(b)中进行了放大展示,图中包含2个动目标,可观察到上面一个目标因为存在切向速度而具有轻微的散焦现象,而下面目标因为存在径向速度导致目标偏移进而在原地形成了阴影。图6(c)和(d)分别为第30帧的视频截图和动目标显示,对比两帧图像可观察动目标的变化情况。通过连续图像的渐变显示,视频SAR能够多角度详细的显示目标的运动信息,为动目标的检测跟踪提供了保障。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (6)

1.一种基于无插值融合快速后向投影的雷达视频成像方法,其特征在于,所述方法包括如下步骤:
步骤1,设定雷达视线方向为距离向,雷达平台运动方向为方位向,且雷达回波采样时间为快时间,雷达平台运动时间为慢时间;所述雷达平台采用“一走一停”机制,所述雷达在每个慢时间发射信号并接收对应的回波信号;
步骤2,设定雷达成像的方位分辨率为ρa,雷达平台运动速度为v,雷达平台到成像场景的参考距离为Rs,从而确定雷达成像所需运动的距离L;
步骤3,雷达在每个慢时间对回波信号进行采样,得到大小为N×M的回波数据块,其中,N为所述回波数据块的方位向采样点数,且N=L×Fa/v,M为所述回波数据块的距离向采样点数,Fa为雷达发射信号的脉冲重复频率;所述大小为N×M的回波数据块形成一幅方位分辨率为ρa的SAR图像;
步骤4,对所述N×M的回波数据块进行脉冲压缩,得到脉冲压缩后的回波数据;
步骤5,将所述脉冲压缩后的回波信号划分为K(0)个子孔径数据块,且K(0)为2的整数次幂,每个子孔径数据块方位向采样点数为Nsub=N/K(0)
步骤6,对K(0)个子孔径数据块进行成像,得到K(0)个子图像,并对所述K(0)个子图像进行子图像融合,则进行G=log2(K(0))次子图像融合后得到一幅方位分辨率为ρa的SAR图像;
步骤6具体包括如下子步骤:
(6a)每幅子图像建立在以其子孔径中心为与原点的局部直角坐标系中,对于第k个子孔径数据块,其子孔径中心设为
Figure FDA0003417444480000021
以所述子孔径中心为原点建立局部直角坐标系
Figure FDA0003417444480000022
则第k个子图像
Figure FDA0003417444480000023
表示为:
Figure FDA0003417444480000024
其中,k=1,2,…,K(0),x为子孔径数据块中采样点的方位向位置,且-l0/2≤x<l0/2,l0表示第k个子孔径的长度,且l(0)=L/K(0)
Figure FDA0003417444480000025
表示在位置
Figure FDA0003417444480000026
处的脉冲压缩后的回波信号,λ为雷达发射信号的波长,
Figure FDA0003417444480000027
表示雷达在位置
Figure FDA0003417444480000028
时天线相位中心到位置局部直角坐标系
Figure FDA0003417444480000029
的瞬时斜距,且
Figure FDA00034174444800000210
(6b)令k=1,2,…,K(0),重复执行(6a),从而得到K(0)个子图像;
(6c)对所述第1个子图像至第K(0)个子图像两两一组进行融合,则共需要融合次数为G=log2(K(0))次,且每次融合后的子图像数目变为融合前的1/2;
则第g-1次子图像融合后得到的子图像为
Figure FDA00034174444800000211
其中,
Figure FDA00034174444800000212
表示相干融合,
Figure FDA00034174444800000213
表示第g-1次子图像融合的第2k-1个子图像,
Figure FDA00034174444800000214
表示第g-1次子图像融合的第2k个子图像;
(6d)令g=1,2,...G,进行G次子图像融合后得到一幅方位分辨率为ρa的SAR图像;
步骤7,设雷达新接收到大小为Nnew×M的回波数据块,且Nnew<N,并对所述大小为Nnew×M的回波数据块进行脉冲压缩,得到脉冲压缩后的新的回波数据;
步骤8,将所述脉冲压缩后的新的回波数据划分为方位向采样点数为Nsub的T个子孔径数据块,且T=floor[Nnew/Nsub],表示向下取整操作;
步骤9,对所述T个子孔径数据块进行成像,得到T子图像,获取前一幅SAR图像K(0)个子图像中的后K(0)-T个子图像;从而将所述T个子图像和K(0)-T个子图像进行子图像融合,得到一幅方位分辨率为ρa的新的SAR图像;
步骤10,重复执行步骤7至步骤9达到设定次数,得到多幅SAR图像,所述多幅SAR图像形成视频成像。
2.根据权利要求1所述的一种基于无插值融合快速后向投影的雷达视频成像方法,其特征在于,步骤2中,确定雷达成像所需运动的距离L即合成孔径
Figure FDA0003417444480000031
其中,λ为雷达发射信号的波长。
3.根据权利要求1所述的一种基于无插值融合快速后向投影的雷达视频成像方法,其特征在于,步骤3中,雷达发射信号的脉冲重复频率Fa=ka×2v/Da,其中,ka为过采样率,Da为雷达方位向孔径长度,从而雷达在慢时间停留的位置x′=n*v/Fa,n=1,2,3,……。
4.根据权利要求1所述的一种基于无插值融合快速后向投影的雷达视频成像方法,其特征在于,在子步骤(6c)中对子图像两两一组进行融合具体包括:
(a)融合前的一幅子图像的子孔径中心为O1,子孔径长度为l(g),融合后的子图像的子孔径中心为OA,子孔径长度为l(g+1)
(b)设点目标P在融合前的一个坐标系中的坐标为(x1,y1),在融合后的坐标系中的坐标为(xA,yA),则点目标P融合前到融合后的平移量(△x,△y)为:
Figure FDA0003417444480000041
(c)由于O1OA=(l(g+1)-l(g))/2=l(g+1)/4,则
Figure FDA0003417444480000042
其中,l(g+1)为融合后的子孔径长度,θ为斜视角,且l(g)=2(g)l0,l0为初始子孔径长度l0=L/K(0),L为合成孔径长度;
(d)融合前两个子图像的子孔径中心关于融合后的子图像的子孔径中心对称,从而确定融合前两个子图像的平移量分别为(△x,△y)和(-△x,-△y)。
5.根据权利要求4所述的一种基于无插值融合快速后向投影的雷达视频成像方法,其特征在于,在频域实现两个子图像平移的过程为:
采用FFT将两个子图像分别变换到二维频域,其中一个子图像在二维频域乘以第一线性相位,另一个子图像在二维频域乘以第二线性相位:
所述第一线性相位的表达式为:
Figure FDA0003417444480000043
所述第二线性相位的表达式为:
Figure FDA0003417444480000044
其中,S1x和S1y分别为一个子图像在方位频域和距离频域相乘的线性相位,S2x和S2y分别为另一个子图像在方位频域和距离频域相乘的线性相位,fa为方位频率,fr为距离频率。
6.根据权利要求1所述的一种基于无插值融合快速后向投影的雷达视频成像方法,其特征在于,设相邻两幅SAR图像的数据重叠率为η,雷达新接收的回波数据块方位向点数Nnew=N(1-η),则多幅SAR图像形成的视频成像的最大帧率Fmax为:
Fmax=Fa/Nnew=v/(L(1-η))
其中,Fa表示脉冲重复频率,v表示雷达平台运动速度,L表示雷达成像所需运动的距离。
CN201810058866.6A 2018-01-22 2018-01-22 基于无插值融合快速后向投影的雷达视频成像方法 Active CN108205135B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810058866.6A CN108205135B (zh) 2018-01-22 2018-01-22 基于无插值融合快速后向投影的雷达视频成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810058866.6A CN108205135B (zh) 2018-01-22 2018-01-22 基于无插值融合快速后向投影的雷达视频成像方法

Publications (2)

Publication Number Publication Date
CN108205135A CN108205135A (zh) 2018-06-26
CN108205135B true CN108205135B (zh) 2022-03-04

Family

ID=62606305

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810058866.6A Active CN108205135B (zh) 2018-01-22 2018-01-22 基于无插值融合快速后向投影的雷达视频成像方法

Country Status (1)

Country Link
CN (1) CN108205135B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109188432B (zh) * 2018-08-08 2022-07-29 中国航空工业集团公司雷华电子技术研究所 一种平行双基聚束sar快速bp成像方法
CN109633638B (zh) * 2018-11-19 2020-11-03 上海无线电设备研究所 一种基于电磁散射模型的大角度sar超分辨成像方法
CN110095775B (zh) * 2019-04-29 2023-03-14 西安电子科技大学 基于混合坐标系的颠簸平台sar快速时域成像方法
CN110068833B (zh) * 2019-05-05 2021-10-29 中国科学院电子学研究所 一种合成孔径激光雷达成像方法、仪器及系统
CN110488289B (zh) * 2019-08-21 2022-07-15 苏州兴钊防务研究院有限公司 基于重叠子孔径的光电协同合成孔径雷达成像处理方法
CN111352108B (zh) * 2020-02-28 2022-11-29 南昌大学 基于ffbp逆向处理的sar回波信号快速仿真方法
CN111736151B (zh) * 2020-06-16 2022-03-04 西安电子科技大学 一种高效全局直角坐标投影融合的改进ffbp成像方法
CN113514827B (zh) * 2021-03-03 2023-09-05 南昌大学 无人机集群模式下合成孔径雷达成像处理方法及应用
CN113189588B (zh) * 2021-04-30 2022-05-03 电子科技大学 集群无人机合成孔径雷达高帧率成像方法
CN113109813A (zh) * 2021-05-26 2021-07-13 南京牛斯特仁科技有限公司 一种基于迭代近端投影的二维欠采样合成孔径雷达成像方法
CN114609635B (zh) * 2022-03-17 2023-06-20 电子科技大学 一种基于视频合成孔径雷达的干涉测量方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105842694A (zh) * 2016-03-23 2016-08-10 中国电子科技集团公司第三十八研究所 一种基于ffbp sar成像的自聚焦方法
CN106802416A (zh) * 2017-02-21 2017-06-06 电子科技大学 一种快速因式分解后向投影sar自聚焦方法
CN107229051A (zh) * 2017-05-26 2017-10-03 西安电子科技大学 基于gpu的视频sar回波仿真并行实现方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105842694A (zh) * 2016-03-23 2016-08-10 中国电子科技集团公司第三十八研究所 一种基于ffbp sar成像的自聚焦方法
CN106802416A (zh) * 2017-02-21 2017-06-06 电子科技大学 一种快速因式分解后向投影sar自聚焦方法
CN107229051A (zh) * 2017-05-26 2017-10-03 西安电子科技大学 基于gpu的视频sar回波仿真并行实现方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"A Fast BP Algorithm With Wavenumber Spectrum Fusion for High-Resolution Spotlight SAR Imaging";Lei Zhang et al.;《IEEE Geoscience and Remote Sensing Letters》;20140131;全文 *
"An Adaptive Fast Factorized Back-Projection Algorithm With Integrated Target Detection Technique for High-resolution and High-Squint Spotlight SAR Imagery";Lei Ran et al.;《IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing》;20171231;全文 *
"Processing video-SAR data with the fast backprojection method";Xiaoshen Song et al.;《IEEE Transactions on Aerospace and Electronic Systems》;20161231;全文 *
"一种加速时域成像算法及其自聚焦方法";曾乐天 等;《西安电子科技大学学报》;20171231;第44卷(第1期);全文 *

Also Published As

Publication number Publication date
CN108205135A (zh) 2018-06-26

Similar Documents

Publication Publication Date Title
CN108205135B (zh) 基于无插值融合快速后向投影的雷达视频成像方法
CN105842694B (zh) 一种基于ffbp sar成像的自聚焦方法
CN104931966B (zh) 一种基于dcs算法的星载视频sar成像处理方法
CN102608597B (zh) 基于非完全数据解卷积的实孔径前视成像方法
CN109738894B (zh) 一种用于大视场合成孔径雷达的高斜视多角度成像方法
CN109143237B (zh) 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法
CN104698459B (zh) 一种应用于缺失数据的条带sar压缩感知成像方法
CN109814100A (zh) 基于子孔径参数估计的sar地面运动目标成像方法
CN114545411B (zh) 一种基于工程实现的极坐标格式多模高分辨sar成像方法
CN111736151B (zh) 一种高效全局直角坐标投影融合的改进ffbp成像方法
CN103969628A (zh) 一种基于压缩感知的合成孔径雷达pfa成像算法
CN101984363A (zh) 一种步进调频体制超高分辨率sar成像方法
CN108008387B (zh) 一种机载阵列天线下视三维成像方法
CN108957452A (zh) 一种合成孔径雷达自适应ffbp成像方法
CN107102328B (zh) 基于fpga的实时成像信号处理方法及fpga
CN112433210A (zh) 一种双站前视探地雷达快速时域成像方法
CN115685200A (zh) 一种高精度大前斜视sar成像运动补偿与几何校正方法
CN113484859B (zh) 一种基于融合技术的二维超分辨雷达成像方法
CN114325704A (zh) 一种基于波数谱拼接的合成孔径雷达快速时域成像方法
CN117129994B (zh) 基于相位补偿核gnss-sar改进后向投影成像方法
CN112285707B (zh) 基于gps导航信号的无源多基地高分辨成像方法
CN107229050B (zh) 一种基于极坐标格式的雷达成像优化方法
CN116359921A (zh) 基于加速轨迹双基前视合成孔径雷达的快速时域成像方法
CN115015920A (zh) 一种基于距离空变频谱矫正的快速后向投影成像方法
CN113109813A (zh) 一种基于迭代近端投影的二维欠采样合成孔径雷达成像方法

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