CN110045373B - 机载多通道sar成像处理方法、装置和计算机设备 - Google Patents

机载多通道sar成像处理方法、装置和计算机设备 Download PDF

Info

Publication number
CN110045373B
CN110045373B CN201910279041.1A CN201910279041A CN110045373B CN 110045373 B CN110045373 B CN 110045373B CN 201910279041 A CN201910279041 A CN 201910279041A CN 110045373 B CN110045373 B CN 110045373B
Authority
CN
China
Prior art keywords
signal
compensation factor
echo
compensation
acquiring
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
CN201910279041.1A
Other languages
English (en)
Other versions
CN110045373A (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201910279041.1A priority Critical patent/CN110045373B/zh
Publication of CN110045373A publication Critical patent/CN110045373A/zh
Application granted granted Critical
Publication of CN110045373B publication Critical patent/CN110045373B/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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Abstract

本发明涉及一种机载多通道SAR成像处理方法,包括:获取SAR的成像参数、回波数据和载机的惯导数据,根据惯导数据计算载机的运动误差;根据成像参数和运动误差,对回波数据进行通道间的幅度误差、相位误差和空不变运动误差补偿,得到补偿回波信号;根据成像参数对补偿回波信号进行方位向重构滤波,得到距离多普勒域信号;获取第一补偿因子,根据距离多普勒域信号和第一补偿因子得到第一信号;对第一信号进行距离向傅里叶变换,和第二补偿因子相乘并进行距离向傅里叶逆变换得到第四信号;对第四信号在二维时域补偿空变运动误差得到第七信号;获取第三补偿因子,与第七信号相乘并进行方位向傅里叶逆变换得到处理后的SAR图像。

Description

机载多通道SAR成像处理方法、装置和计算机设备
技术领域
本发明涉及信号处理领域,特别是涉及一种机载多通道SAR成像 处理方法、装置和计算机设备。
背景技术
SAR(Synthetic Aperture Radar,机载合成孔径雷达)是一种高分辨 率成像雷达,多通道SAR成像技术可以实现高分辨率宽覆盖成像。然 而,多个接收子天线之间的相位和幅度误差会在最终成像结果中引入 虚假目标,严重影响成像质量,此外受气流和飞机控制系统的影响, 会生成较大的运动误差,因此需要对各个通道的回波数据进行幅度和 误差校正,并且进行运动补偿。
传统的基于惯导数据的运动补偿方法,第一步补偿参考斜距处的 误差,第二步补偿随斜距变化的误差。对于多通道机载SAR数据来说, 多个通道的回波数据共用一个惯导数据,在进行第一步补偿时,可以 认为各个通道的运动误差相同,然而在回波重构之后,有几个通道的 回波数据,方位向数据就变为原来的几倍,原始运动误差信息不再适 用,因此传统的运动补偿方法得到的成像质量不够高。
发明内容
本发明的目的在于提供一种机载多通道SAR成像处理方法、装置、 计算机设备和可读存储介质,可以有效提高SAR成像质量。
本发明的目的通过如下技术方案实现:
一种机载多通道SAR成像处理方法,方法包括:
获取SAR的成像参数、回波数据和载机的惯导数据,根据惯导数 据计算载机的运动误差;
根据成像参数和运动误差,对回波数据进行补偿通道间的幅度误 差和相位误差,并进行补偿空不变运动误差,得到补偿回波信号;
根据成像参数对补偿回波信号进行方位向重构滤波,得到距离多 普勒域信号;获取第一补偿因子,根据距离多普勒域信号和第一补偿 因子得到第一信号;
对所述第一信号进行距离向傅里叶变换,获取第二补偿因子,将 所述第一信号与所述第二补偿因子相乘并进行距离向傅里叶逆变换得 到第四信号;
对所述第四信号在二维时域补偿空变运动误差得到第七信号;
获取第三补偿因子,将所述第三补偿因子与所述第七信号相乘并 进行方位向傅里叶逆变换得到处理后的SAR图像。
一种机载多通道SAR成像处理装置,装置包括:
运动误差获取模块,用于获取SAR的成像参数、回波数据和载机 的惯导数据,根据惯导数据计算载机的运动误差;
补偿回波信号获取模块,用于根据成像参数和运动误差,对回波 数据进行补偿通道间的幅度误差和相位误差,并进行补偿空不变运动 误差,得到补偿回波信号;
第一信号获取模块,用于根据成像参数对补偿回波信号进行方位 向重构滤波,得到距离多普勒域信号;获取第一补偿因子,根据距离 多普勒域信号和第一补偿因子得到第一信号;
第四信号获取模块,用于对所述第一信号进行距离向傅里叶变换, 获取第二补偿因子,和将所述第一信号与所述第二补偿因子相乘并进 行距离向傅里叶逆变换得到第四信号
第七信号获取模块,用于对所述第四信号在二维时域补偿空变运 动误差得到第七信号;
SAR处理图像获取模块,用于获取第三补偿因子,将所述第三补 偿因子与所述第七信号相乘并进行方位向傅里叶逆变换得到处理后的 SAR图像。
一种计算机设备,包括存储器和处理器,存储器存储有计算机程 序,处理器执行计算机程序时上述步骤。
一种计算机可读存储介质,其上存储有计算机程序,计算机程序 被处理器执行时实现上述步骤。
本发明提供的机载多通道SAR成像处理方法,其获取SAR的成 像参数、回波数据和载机的惯导数据,根据所述惯导数据计算载机的 运动误差;根据所述成像参数和所述运动误差,对所述回波数据进行 补偿通道间的幅度误差和相位误差,并进行补偿空不变运动误差,得 到补偿回波信号;根据所述成像参数对所述补偿回波信号进行方位向 重构滤波,得到距离多普勒域信号;获取第一补偿因子,根据所述距 离多普勒域信号和所述第一补偿因子得到第一信号;对所述第一信号 进行距离向傅里叶变换,获取第二补偿因子,将所述第一信号与所述 第二补偿因子相乘并进行距离向傅里叶逆变换得到第四信号;对所述 第四信号在二维时域补偿空变运动误差得到第七信号;获取第三补偿 因子,将所述第三补偿因子与所述第七信号相乘并进行方位向傅里叶 逆变换得到处理后的SAR图像,可以有效提高SAR成像质量。
附图说明
图1为一个实施例中机载多通道SAR成像处理方法的应用环境 图;
图2为一个实施例中机载多通道SAR成像处理方法的流程示意 图;
图3为一个实施例中机载SAR空间几何模型示意图;
图4为一个实施例中机载SAR运动误差示意图;
图5为一个实例中不采用运动补偿和幅相误差补偿的情况下回波 数据的频谱图;
图6为一个实例中经过运动补偿和幅相误差补偿后回波数据的频 谱图;
图7为另一个实施例中机载多通道SAR成像处理装置的结构框 图;
图8为一个实施例中计算机设备的内部结构图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下结合附 图及实施例,对本发明进行进一步的详细说明。应当理解,此处所描 述的具体实施方式仅仅用以解释本发明,并不限定本发明的保护范围。
本申请提供的机载多通道SAR成像处理方法可以应用于如图1所 示的应用环境中。该应用环境包括服务器,服务器获取SAR的成像参 数、回波数据和载机的惯导数据,根据所述惯导数据计算载机的运动 误差;
根据所述成像参数和所述运动误差,对所述回波数据进行补偿通 道间的幅度误差和相位误差,并进行补偿空不变运动误差,得到补偿 回波信号;
根据所述成像参数对所述补偿回波信号进行方位向重构滤波,得 到距离多普勒域信号;获取第一补偿因子,根据所述距离多普勒域信 号和所述第一补偿因子得到第一信号;服务器对所述第一信号进行距 离向傅里叶变换,获取第二补偿因子,将所述第一信号与所述第二补 偿因子相乘并进行距离向傅里叶逆变换得到第四信号;服务器对所述 第四信号在二维时域补偿空变运动误差得到第七信号;服务器获取第 三补偿因子,将所述第三补偿因子与所述第七信号相乘并进行方位向 傅里叶逆变换得到处理后的SAR图像。其中,服务器可以用独立的服 务器或者是多个服务器组成的服务器集群来实现。
在一个实施例中,如图2所示,提供了一种机载多通道SAR成像 处理方法,包括以下步骤:
步骤S101,获取SAR的成像参数、回波数据和载机的惯导数据, 根据惯导数据计算载机的运动误差。
其中,成像参数包括载波波长λ,信号脉宽Tp,信号带宽BW,采 样率fs,脉冲重复频率fprf,波束中心下视角θc,斜视角
Figure BDA0002021061810000061
采样器起 始时刻Tstart,载机高度h,方位向脉冲数Na,距离向点数Nr,通道数m, 载机速度v,各子天线相位中心与参考天线相位中心的距离xi, i=0,1,2,...,m-1;每个通道的回波数据SiK,tL),i=0,1,2,3,K=0,1,...,Nr-1, L=0,1,2,...,Na-1;惯导数据包括载机东向速度
Figure BDA0002021061810000062
载机北向速度
Figure BDA0002021061810000063
载机天向度
Figure BDA0002021061810000064
在具体实施过程中,多通道成像时,一个发射天线发射脉冲,多 个接收天线同时接收脉冲,根据数据存储格式,读取原始数据文件, 得到成像参数,回波数据和载机的惯导数据。
步骤S102,根据成像参数和运动误差,对回波数据进行补偿通道 间的幅度误差和相位误差,并进行补偿空不变运动误差,得到补偿回 波信号。
其中,多通道机载SAR为单天线发射脉冲,多天线接收脉冲,惯 导数据记录下了每个发射脉冲时的飞行速度,包括载机的东北天速度
Figure BDA0002021061810000065
在具体实施过程中,采用幅度均衡的方法对幅度误差进行估计, 采用ESPRIT法对相位误差进行估计;对于各个通道的回波数据来说, 同一个方位向时间的运动误差相同,因此对各通道的回波信号做相同 的运动补偿。
步骤S103,根据成像参数对补偿回波信号进行方位向重构滤波, 得到距离多普勒域信号;获取第一补偿因子,根据距离多普勒域信号 和第一补偿因子得到第一信号。
在具体实施过程中,对于点目标信号来说,距离多普勒域信号与 第一补偿因子相乘之后得到第一信号。
步骤S104,对所述第一信号进行距离向傅里叶变换,获取第二补 偿因子,将所述第一信号与所述第二补偿因子相乘并进行距离向傅里 叶逆变换得到第四信号。
在具体实施过程中,对第一信号进行距离向傅里叶变换,将信号 变换到二维频域,得到第二信号,对第二信号在二维频域通过乘以第 二补偿因子完成一致距离徙动校正和距离聚焦处理,得到第三信号, 将第三信号进行距离向傅里叶逆变换得到第四信号。
步骤S105,对所述第四信号在二维时域补偿空变运动误差得到第 七信号。
在具体实施过程中,由于第四信号经过了距离向聚焦和距离徙动 校正,因此可以对每个距离门的回波信号补偿空变误差,由于方位向 点数范围已经变为原来的m倍,因此需要对运动误差进行插值,得到 第七信号。
步骤S106,获取第三补偿因子,将所述第三补偿因子与所述第七 信号相乘并进行方位向傅里叶逆变换得到处理后的SAR图像。
具体的,获取第三补偿因子,根据第七信号和第三补偿因子得到 第八信号,对第八信号进行方位向傅里叶逆变换,得到处理后的SAR 图像。
本发明提供的机载多通道SAR成像处理方法,通过获取SAR的 成像参数、回波数据和载机的惯导数据,根据所述惯导数据计算载机 的运动误差;根据所述成像参数和所述运动误差,对所述回波数据进 行补偿通道间的幅度误差和相位误差,并进行补偿空不变运动误差, 得到补偿回波信号;根据所述成像参数对所述补偿回波信号进行方位 向重构滤波,得到距离多普勒域信号;获取第一补偿因子,根据所述 距离多普勒域信号和所述第一补偿因子得到第一信号;对所述第一信 号进行距离向傅里叶变换,获取第二补偿因子,将所述第一信号与所 述第二补偿因子相乘并进行距离向傅里叶逆变换得到第四信号;对所 述第四信号在二维时域补偿空变运动误差得到第七信号;获取第三补 偿因子,将所述第三补偿因子与所述第七信号相乘并进行方位向傅里 叶逆变换得到处理后的SAR图像,可以有效提高SAR成像质量。
在一个实施例中,惯导数据包括载机速度,载机速度包括东向速 度、北向速度和天向速度;根据惯导数据计算载机的运动误差,包括:
A1,通过预设的均值滤波对载机速度进行平滑处理,得到预处理 速度。
如图3所示,图3为机载SAR成像空间几何关系,根据已知的数 据存储格式,读取成像参数,回波数据和惯导数据,表1给出了本实 施例中的成像参数:
表1成像参数
Figure BDA0002021061810000081
Figure BDA0002021061810000091
各通道回波数据SiK,tL),i=0,1,2,3,K=0,1,...,Nr-1, L=0,1,2,...,Na-1;惯导数据包括载机东向速度
Figure BDA0002021061810000092
载机北向速度
Figure BDA0002021061810000093
载机天向度
Figure BDA0002021061810000094
A2,对预处理速度积分得到载机位置,根据载机位置得到运动误 差位置坐标,对运动误差位置坐标进行坐标系转换,得到载机的运动 误差。
在具体实施过程中,多通道机载SAR为单天线发射脉冲,多天线 接收脉冲,其中垂直于航迹向的运动误差如图4所示,惯导数据记录 下了每个发射脉冲时的飞行速度,包括载机的东北天速度
Figure BDA0002021061810000095
首先通过均值滤波对速度进行平滑处理,然后对速度积分得到位置, 将东北天坐标系下的运动误差位置坐标转化为成像坐标系下的运动误 差位置坐标,得到最终的运动误差Δx(tL),Δz(tL),L=0,1,2,...,Na-1。
在一个实施例中,根据成像参数和运动误差,对回波数据进行补 偿通道间的幅度误差和相位误差,并进行补偿空不变运动误差,得到 补偿回波信号,包括:
B1,根据成像参数和回波数据,获取通道补偿因子。
具体实施过程中,根据成像参数和回波数据,计算每一个通道信 号的平均幅度;获取通道的相位误差,根据平均幅度和相位误差计算 通道补偿因子。
B2,根据通道补偿因子和回波数据计算得到通道补偿回波信号。
具体实施步骤如下:
(1)计算每一个通道信号的平均幅度值;
Figure BDA0002021061810000101
其中,
Figure BDA0002021061810000102
表示第i个通道的平均幅度值;
(2)采用ESPRIT方法对相位误差进行估计,以第一个通道为参 考,得到其余通道的相位误差为σi,i=0,1,2,3;
(3)以第一个通道为参考,其余各通道的回波数据乘以各通道补 偿因子,同时补偿幅度和相位误差,通道补偿因子如下:
Figure BDA0002021061810000103
得到补偿完成之后的通道补偿回波信号SiK,tL),i=0,1,2,3;
B3,获取运动误差空不变补偿因子,对各个通道的通道补偿回波 信号进行距离向傅里叶变换得到距离频域回波信号。
具体包括以下几个步骤;
(1)计算每一帧回波数据在参考下视角θc下的误差;
Δr(τL;θc)=Δz(τL)·cosθc-Δx(τL)·sinθc (3)
其中:
Figure BDA0002021061810000111
Figure BDA0002021061810000112
R(tc)表示波束中心对应的斜距;
(2)计算每一帧回波数据的补偿因子;
Figure BDA0002021061810000113
(3)对各个通道回波数据进行距离向傅里叶变换得到距离频域回 波信号Si(fτK,tL),i=0,1,2,3;
B4,根据距离频域回波信号和空不变补偿因子获取补偿回波信号。
(4)将回波信号与空不变补偿因子相乘,补偿空不变运动误差, 补偿公式如下所示:
Si(fτK,tL)=Si(fτK,tL)·H1(fτK,tL),i=0,1,2,3 (7)
(5)对回波信号进行距离向傅里叶逆变换重新得到补偿回波信号 SiK,tL),i=0,1,2,3。
在一个实施例中,根据成像参数对补偿回波信号进行方位向重构 滤波,得到距离多普勒域信号;获取第一补偿因子,根据距离多普勒 域信号和第一补偿因子得到第一信号,包括:
C1,根据成像参数构造重构滤波器,根据重构滤波器和补偿回波 信号获取单通道频谱和距离多普勒域信号。
多通道重构滤波主要包括如下几个步骤:
(1)构造重构滤波器G(f),大小为4×4
Figure BDA0002021061810000121
其中:
Figure BDA0002021061810000122
Figure BDA0002021061810000123
(2)对各通道回波信号进行方位向3倍补零,对于每个通道的每 个距离向的数据来说,数据量扩大为原来的4倍;
Figure BDA0002021061810000124
(3)对各个通道的信号进行方位向傅里叶变换得到Si'(τK,fg)
(4)对于每一个距离向时间τK,计算该距离向时间下的重构滤波 之后的频谱S(τK,fg),计算方法如下:
Figure BDA0002021061810000125
k=0时,计算得到S(τK,fg),g=0,1,2,...,Na-1,此时得到重构之后的 第一段频谱,k=1时,计算得到S(τK,fg+fprf),g=0,1,2,...,Na-1,此时得 到重构之后的第二段频谱,依次循环,最后将四段频谱拼接起来即可 恢复等效单通道完整频谱,而后计算下一个距离向时间下的方位向频 谱,重构滤波后得到距离多普勒域的信号:
Figure BDA0002021061810000131
C2,根据成像参数计算第一补偿因子。
C3,将距离多普勒域信号和第一补偿因子相乘得到第一信号。
在具体实施过程中,第一补偿因子如公式(14)所示:
Figure BDA0002021061810000132
其中:
Figure BDA0002021061810000133
Figure BDA0002021061810000134
Figure BDA0002021061810000135
Figure BDA0002021061810000136
Figure BDA0002021061810000137
机载SAR中,认为
Figure BDA0002021061810000138
得到距离多普勒域信号与第一补偿因 子相乘之后的信号如公式(20)所示:
Figure BDA0002021061810000141
Figure BDA0002021061810000142
其中,c为光速,b为调频率,σ为点目标后向散射特性,Wa为方 位向天线方向性函数,a(τK)是发射脉冲包络。
在一个实施例中,对第一信号进行距离向傅里叶变换得到第二信 号,获取第二补偿因子,根据第二补偿因子和第二信号得到第三信号, 包括:
D1,对第一信号进行距离向傅里叶变换,将第一信号换到二维频 域,得到第二信号。
在具体实施过程中,对公式(20)进行距离向傅里叶变换,将信 号变换到二维频域,得到第二信号如公式(22)所示:
Figure BDA0002021061810000151
其中:
Figure BDA0002021061810000152
Figure BDA0002021061810000153
D2,获取第二补偿因子,将第二补偿因子和第二信号相乘得到第 三信号。
在具体实施过程中,第二信号S2(fτK,fg)在二维频域通过乘以距离 向补偿因子,即第二补偿因子,完成一致距离徙动校正和距离聚焦处 理,第二补偿因子如公式(25)所示:
Figure BDA0002021061810000154
补偿完之后的第三信号表示为S3(fτK,fg)。
在一个实施例中,对第三信号进行距离向傅里叶逆变换得到第四 信号,具体包括:
对第三信号作距离向傅里叶逆变换得到距离-多普勒域信号,即第 四信号S4K,fg):
Figure BDA0002021061810000161
其中,A(·)表示距离压缩后的距离向包络,第一个相位为方位向 调频信号,第二、三项为残留相位.
在一个实施例中,在二维时域对第四信号补偿空变误差,得到第 七信号,包括:
E1,对运动误差进行插值,得到插值运动误差;根据插值运动误 差,计算误差补偿因子。
E2,将第四信号进行方位向傅里叶逆变换,得到第五信号;
E3,将第五信号和误差补偿因子相乘得到第六信号,对第六信号 进行方位向傅里叶变换得到第七信号。
此时的信号经过了距离向聚焦和距离徙动校正,因此可以对每个 距离门的回波信号补偿空变误差,由于方位向点数范围已经变为原来 的m倍,因此需要对运动误差进行插值,具体包括以下几个步骤;
(1)对步骤A2得到的运动误差进行插值,得到插值之后的插值 运动误差Δx(tg),Δz(tg),g=0,1,2,...,m·Na-1;
(2)计算每一个方位向时刻的回波数据在不同下视角θK下的误 差;
Δr'(tg;θK)=Δz(tg)·(cosθK-cosθc)-Δx(tg)·(sinθK-sinθc)(27)
其中:
Figure BDA0002021061810000171
Figure BDA0002021061810000172
(3)计算误差补偿因子;
Figure BDA0002021061810000173
(4)对第四信号进行方位向傅里叶逆变换,得到二维时域信号, 即第五信号S5K,tg);
(5)将回波信号与误差补偿因子相乘,补偿空变运动误差,得到 第六信号;
S6K,tg)=S5K,tg)·H2K,tg) (31)
(4)对回波信号进行方位向傅里叶变换得到距离多普勒域信号, 即第七信号S7K,fg)。
在一个实施例中,获取第三补偿因子,将所述第三补偿因子与所 述第七信号相乘并进行方位向傅里叶逆变换得到处理后的SAR图像, 具体包括:
在距离多普勒域乘以方位向补偿因子完成方位向聚焦处理和和残 留相位补偿,第三补偿因子如公式(32)所示:
Figure BDA0002021061810000174
其中第一项完成方位向聚焦,第二、三项补偿残余相位,同时对 斜距为R的目标保留相位
Figure BDA0002021061810000181
此时得到距离多普勒域信号,即第八 信号如公式(33)所示:
Figure BDA0002021061810000182
将第八信号经过方位向傅里叶逆变换,得到SAR图像如公式(34) 所示:
Figure BDA0002021061810000183
其中,Wac(·)为方位天线Wa(·)变换后的包络。
图5是不进行运动误差补偿和幅相误差补偿的情况下,回波数据 的方位向频谱图,图6是采用本文中提出的方法进行运动误差补偿和 幅相误差补偿之后回波数据的频谱图,从图中可以看出,经过运动误 差补偿和幅相误差补偿,方位向频谱质量明显变好,成像质量有所提 高。
应该理解的是,虽然图2的流程图中的各个步骤按照箭头的指示 依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。 除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制, 这些步骤可以以其它的顺序执行。而且,图2中的至少一部分步骤可 以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在 同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶 段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步 骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在其中一个实施例中,如图7所示,提供了一种机载多通道SAR 成像处理装置,装置包括:
运动误差获取模块701,用于获取SAR的成像参数、回波数据和 载机的惯导数据,根据惯导数据计算载机的运动误差;
补偿回波信号获取模块702,用于根据成像参数和运动误差,对回 波数据进行补偿通道间的幅度误差和相位误差,并进行补偿空不变运 动误差,得到补偿回波信号;
第一信号获取模块703,用于根据成像参数对补偿回波信号进行方 位向重构滤波,得到距离多普勒域信号;获取第一补偿因子,根据距 离多普勒域信号和第一补偿因子得到第一信号;
第四信号获取模块704,用于对所述第一信号进行距离向傅里叶变 换,获取第二补偿因子,和将所述第一信号与所述第二补偿因子相乘 并进行距离向傅里叶逆变换得到第四信号;
第七信号获取模块705,用于对所述第四信号在二维时域补偿空变 运动误差得到第七信号;
SAR处理图像获取模块706,用于获取第三补偿因子,将所述第 三补偿因子与所述第七信号相乘并进行方位向傅里叶逆变换得到处理 后的SAR图像。
关于机载多通道SAR成像处理装置的具体限定可以参见上文中对 于机载多通道SAR成像处理方法的限定,在此不再赘述。上述机载多 通道SAR成像处理装置中的各个模块可全部或部分通过软件、硬件及 其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备 中的处理器中,也可以以软件形式存储于计算机设备中的存储器中, 以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是 服务器,其内部结构图可以如图8所示。该计算机设备包括通过装置 总线连接的数据处理器、存储器、网络接口和数据库。其中,该计算 机设备设置有多个数据处理器,数据处理器用于提供计算和控制能力。 该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失 性存储介质存储有操作装置、计算机程序和数据库。该内存储器为非 易失性存储介质中的操作装置和计算机程序的运行提供环境。该计算 机设备的数据库用于存储机载多通道SAR成像处理涉及的数据。该计 算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种机载多通道SAR成像处理方法。
本领域技术人员可以理解,图8中示出的结构,仅仅是与本申请 方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的 计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更 少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器, 存储器中存储有计算机程序,该处理器执行计算机程序时实现上述机 载多通道SAR成像处理方法。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有 计算机程序,计算机程序被处理器执行时实现机载多通道SAR成像处 理方法。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部 分流程,是可以通过计算机程序来指令相关的硬件来完成,的计算机 程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在 执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供 的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引 用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只 读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、 电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机 存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限, RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、 同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型 SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、 存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态 RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁, 未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而, 只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的 范围。
以上实施例仅表达了本发明的几种实施方式,其描述较为具体和 详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是, 对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还 可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本 发明专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种机载多通道SAR成像处理方法,其特征在于,所述方法包括:
获取SAR的成像参数、回波数据和载机的惯导数据,根据所述惯导数据计算载机的运动误差;
根据所述成像参数和所述运动误差,对所述回波数据进行补偿通道间的幅度误差和相位误差,并进行补偿空不变运动误差,得到补偿回波信号;
根据所述成像参数对所述补偿回波信号进行方位向重构滤波,得到距离多普勒域信号;获取第一补偿因子,根据所述距离多普勒域信号和所述第一补偿因子得到第一信号;
对所述第一信号进行距离向傅里叶变换,获取第二补偿因子,将所述第一信号与所述第二补偿因子相乘并进行距离向傅里叶逆变换得到第四信号;
对所述第四信号在二维时域补偿空变运动误差得到第七信号;
获取第三补偿因子,将所述第三补偿因子与所述第七信号相乘并进行方位向傅里叶逆变换得到处理后的SAR图像。
2.根据权利要求1所述的方法,其特征在于,所述惯导数据包括载机速度,所述载机速度包括东向速度、北向速度和天向速度;所述根据所述惯导数据计算载机的运动误差,包括:
通过预设的均值滤波对所述载机速度进行平滑处理,得到预处理速度;
对所述预处理速度积分得到载机位置,根据所述载机位置得到运动误差位置坐标,对所述运动误差位置坐标进行坐标系转换,得到所述载机的运动误差。
3.根据权利要求1所述的方法,其特征在于,所述根据所述成像参数和所述运动误差,对所述回波数据进行补偿通道间的幅度误差和相位误差,并进行补偿空不变运动误差,得到补偿回波信号,包括:
根据所述成像参数和所述回波数据,获取通道补偿因子;
根据所述通道补偿因子和所述回波数据计算得到通道补偿回波信号;
获取空不变补偿因子,对各个通道的通道补偿回波信号进行距离向傅里叶变换得到距离频域回波信号;
根据所述距离频域回波信号和所述空不变补偿因子获取所述补偿回波信号。
4.根据权利要求3所述的方法,其特征在于,所述根据所述成像参数和所述回波数据,获取通道补偿因子,包括:
根据所述成像参数和所述回波数据,计算每一个通道信号的平均幅度;
获取通道的相位误差,根据所述平均幅度和所述相位误差计算通道补偿因子。
5.根据权利要求1所述的方法,其特征在于,所述根据所述成像参数对所述补偿回波信号进行方位向重构滤波,得到距离多普勒域信号;获取第一补偿因子,根据所述距离多普勒域信号和所述第一补偿因子得到第一信号,包括:
根据所述成像参数构造重构滤波器,根据所述重构滤波器对所述补偿回波信号进行重构滤波,得到所述距离多普勒域信号;
根据所述成像参数计算第一补偿因子;
将所述距离多普勒域信号和所述第一补偿因子相乘得到所述第一信号。
6.根据权利要求1所述的方法,其特征在于,对所述第一信号进行距离向傅里叶变换,获取第二补偿因子,将所述第一信号与所述第二补偿因子相乘并进行距离向傅里叶逆变换得到第四信号,包括:
对所述第一信号进行距离向傅里叶变换,将所述第一信号变换到二维频域,得到第二信号;
获取第二补偿因子,将所述第二补偿因子和所述第二信号相乘得到第三信号;
对所述第三信号进行距离向傅里叶逆变换得到所述第四信号。
7.根据权利要求1所述的方法,其特征在于,所述对所述第四信号在二维时域补偿空变运动误差得到第七信号,包括:
对所述运动误差按照通道数进行插值,得到插值运动误差;根据所述插值运动误差,计算运动误差补偿因子;
将所述第四信号进行方位向傅里叶逆变换,得到第五信号;
将所述第五信号和所述运动误差补偿因子相乘得到第六信号,对所述第六信号进行方位向傅里叶变换得到第七信号。
8.一种机载多通道SAR成像处理装置,其特征在于,所述装置包括:
运动误差获取模块,用于获取SAR的成像参数、回波数据和载机的惯导数据,根据所述惯导数据计算载机的运动误差;
补偿回波信号获取模块,用于根据所述成像参数和所述运动误差,对所述回波数据进行补偿通道间的幅度误差和相位误差,并进行补偿空不变运动误差,得到补偿回波信号;
第一信号获取模块,用于根据所述成像参数对所述补偿回波信号进行方位向重构滤波,得到距离多普勒域信号;获取第一补偿因子,根据所述距离多普勒域信号和所述第一补偿因子得到第一信号;
第四信号获取模块,用于对所述第一信号进行距离向傅里叶变换,获取第二补偿因子,和将所述第一信号与所述第二补偿因子相乘并进行距离向傅里叶逆变换得到第四信号
第七信号获取模块,用于对所述第四信号在二维时域补偿空变运动误差得到第七信号;
SAR处理图像获取模块,用于获取第三补偿因子,将所述第三补偿因子与所述第七信号相乘并进行方位向傅里叶逆变换得到处理后的SAR图像。
9.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。
CN201910279041.1A 2019-04-09 2019-04-09 机载多通道sar成像处理方法、装置和计算机设备 Active CN110045373B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910279041.1A CN110045373B (zh) 2019-04-09 2019-04-09 机载多通道sar成像处理方法、装置和计算机设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910279041.1A CN110045373B (zh) 2019-04-09 2019-04-09 机载多通道sar成像处理方法、装置和计算机设备

Publications (2)

Publication Number Publication Date
CN110045373A CN110045373A (zh) 2019-07-23
CN110045373B true CN110045373B (zh) 2021-03-30

Family

ID=67276430

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910279041.1A Active CN110045373B (zh) 2019-04-09 2019-04-09 机载多通道sar成像处理方法、装置和计算机设备

Country Status (1)

Country Link
CN (1) CN110045373B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111487628B (zh) * 2020-05-19 2022-05-03 中国科学院声学研究所 一种下视多波束合成孔径成像声纳“零度”干扰抑制方法
CN112578380B (zh) * 2020-11-05 2023-09-19 北京理工大学 基于载机航迹偏差估计的无人机载sar自聚焦方法
CN115656948B (zh) * 2022-12-28 2023-04-07 中国科学院空天信息创新研究院 基于静止强目标的相位误差估计方法、装置、设备及介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104316913A (zh) * 2014-11-13 2015-01-28 中国科学院电子学研究所 多通道接收机实时校准装置及校准与误差补偿方法
CN105974415A (zh) * 2016-06-24 2016-09-28 西安电子科技大学 一种机载sar方位空变运动误差的高精度补偿方法
CN108614249A (zh) * 2018-04-12 2018-10-02 北京航空航天大学 相位误差的估计方法、装置、补偿方法和系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10054680B2 (en) * 2014-07-09 2018-08-21 Bae Systems Information And Electronic Systems Integration Inc. Radar detection of endo-clutter high-value targets using tracker feedback

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104316913A (zh) * 2014-11-13 2015-01-28 中国科学院电子学研究所 多通道接收机实时校准装置及校准与误差补偿方法
CN105974415A (zh) * 2016-06-24 2016-09-28 西安电子科技大学 一种机载sar方位空变运动误差的高精度补偿方法
CN108614249A (zh) * 2018-04-12 2018-10-02 北京航空航天大学 相位误差的估计方法、装置、补偿方法和系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
An Interferometric Approach to Ocean Surface Velocity Imaging Using Multi-Channel SAR;Mark Sletten等;《IEEE IGARSS2017》;20171231;全文 *
多通道星载合成孔径雷达姿态误差影响分析与补偿;匡辉等;《电波科学学报》;20141231;第29卷(第6期);全文 *

Also Published As

Publication number Publication date
CN110045373A (zh) 2019-07-23

Similar Documents

Publication Publication Date Title
CN110045373B (zh) 机载多通道sar成像处理方法、装置和计算机设备
EP1503223B1 (en) Estimation and correction of phase for focusing search mode SAR images formed by range migration algorithm
CA2155502C (en) Method of image generation by means of two-dimensional data processing in connection with a radar with synthetic aperture
CN107229048B (zh) 一种高分宽幅sar动目标速度估计与成像方法
CN106802416B (zh) 一种快速因式分解后向投影sar自聚焦方法
CN111856461B (zh) 基于改进pfa的聚束sar成像方法及其dsp实现
CN111722227B (zh) 基于近似观测矩阵的聚束sar压缩感知成像方法
CN109799502B (zh) 一种适用于滤波反投影算法的两维自聚焦方法
CN113093187B (zh) 一种无道路信息辅助的csar地面动目标跟踪方法
CN116500626B (zh) 一种多接收阵元数据收发合置转换方法
CN109425859B (zh) 一种序贯图像成像方法及装置
CN113093186A (zh) 一种基于分块成像的大场景高分辨率成像方法和装置
CN115542323B (zh) Sar运动目标图像快速重聚焦方法、装置和计算机设备
CN113406624B (zh) 一种高分辨率星载sar高效时频混合成像方法及系统
Saeedi et al. Improved navigation-based motion compensation for LFMCW synthetic aperture radar imaging
CN112946649B (zh) 一种适用于任意子孔径长度的pfa成像方法
CN112285708B (zh) 一种船载x波段雷达海浪反演区域处理方法及系统
CN112147606B (zh) 一种基于多通道雷达的旋转目标微动特征提取方法和装置
CN112379344A (zh) 信号补偿方法及装置、设备、存储介质
CN114166111B (zh) 广角形变监测地基合成孔径雷达方法和系统
CN103760561A (zh) 大斜视角sar成像模式下提高方位向非散焦长度的方法
CN115825903A (zh) 一种基于pfa算法实时成像的多普勒中心精补偿方法
CN116879893B (zh) 基于L型基线的WasSAR动目标参数估计方法与装置
Shen et al. Implementation of Fast Convolution Backprojection for Synthetic Aperture Radar
US20230135348A1 (en) Method for Correcting A Synthetic Aperture Radar Antenna Beam Image

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