CN102243300B - 低频合成孔径雷达射频干扰抑制及误差校正方法 - Google Patents

低频合成孔径雷达射频干扰抑制及误差校正方法 Download PDF

Info

Publication number
CN102243300B
CN102243300B CN201110094916A CN201110094916A CN102243300B CN 102243300 B CN102243300 B CN 102243300B CN 201110094916 A CN201110094916 A CN 201110094916A CN 201110094916 A CN201110094916 A CN 201110094916A CN 102243300 B CN102243300 B CN 102243300B
Authority
CN
China
Prior art keywords
time
formula
rfi
frequency domain
ati
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.)
Expired - Fee Related
Application number
CN201110094916A
Other languages
English (en)
Other versions
CN102243300A (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201110094916A priority Critical patent/CN102243300B/zh
Publication of CN102243300A publication Critical patent/CN102243300A/zh
Application granted granted Critical
Publication of CN102243300B publication Critical patent/CN102243300B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供一种低频SAR的RFI抑制及误差校正方法。技术方案包括五个步骤：第一步，基于回波的接收机频率特性误差校正；第二步，基于校正后时频域回波的距离谱均衡抑制和SAR成像；第三步，基于时域图像的配准；第四步，基于时频域的RFI检测及抑制；第五步，残留误差校正。本发明可解决双通道低频SAR干涉图的RFI检测及抑制、各类误差校正问题，提高沿航向干涉方法的性能，为后续相关处理打下良好基础。

Description

低频合成孔径雷达射频干扰抑制及误差校正方法
技术领域
本发明属于低频SAR(Synthetic Aperture Radar，合成孔径雷达)信号处理技术领域，涉及RFI(Radio Frequency Interference，射频干扰)信号处理与沿航向干涉信号处理技术领域，特别涉及适合双通道低频SAR沿航向干涉图的RFI抑制及误差校正方法。
背景技术
对沿航向放置的雷达天线接收的两个信号通道的复图像共轭相乘，可以获得沿航向的干涉图。为了叙述的简明，下文将沿航向的干涉图简称干涉图。干涉图在军用、民用的遥感探测中有多种应用。例如，对干涉图上的干涉幅度和干涉相位进行检测，可以获得地面运动目标指示(Ground Moving TargetIndication，简称GMTI)。
干涉图对两个复图像间存在的误差非常敏感。通道间存在接收机频率特性误差、配准误差、天线方向图误差、垂直基线误差等各类误差。这些误差会降低两个复图像间的相干性，降低利用干涉幅度和干涉相位进行检测的可靠性。各类基于干涉图的应用都需要准确校正通道间的各类误差。
双通道低频SAR的干涉图不仅易受通道间各类误差影响，还会受到具有窄带特性、时变性和空变性的RFI影响。实验表明，RFI会严重影响干涉图上干涉幅度和干涉相位的联合分布，破坏两个复图像间的相干性，甚至造成干涉图无法应用。
目前，低频SAR应用RFI抑制方法主要是为了提高SAR成像质量，并以SAR成像质量评价抑制性能。其需要考虑的主要因素是单通道能量较强RFI信号的幅度抑制。根据RFI信号频率是否已知，RFI抑制方法可以分成两类。一是基于已知RFI信号频率的抑制方法，该类方法难以检测具有时变性、空变性的RFI；二是基于未知RFI信号频率的抑制方法，该类方法难以彻底抑制能量较弱的RFI。第二类方法应用自适应滤波器时还有计算量大，实时实现困难的问题。目前，尚未有公开文献提供一种抑制干涉图上RFI的方法。因此，双通道低频SAR系统应用中，针对干涉图抑制RFI是一个急需解决的工程技术问题。
发明内容
本发明提供一种适合低频SAR的RFI抑制及误差校正方法，解决双通道低频SAR干涉图的RFI检测及抑制、各类误差校正问题，提高沿航向干涉方法的性能。
本发明的流程分为五步：
第一步，基于回波的接收机频率特性误差校正。
将通道i(i＝1，2)接收的方位时域距离时域回波Ri(tm，tf)进行傅立叶变换，得到通道i方位时域距离频域回波Ri(tm，fr)。式中tm表示方位向慢时间，tf表示距离向快时间，fr表示距离频率。为了叙述的简明，后文将方位时域距离时域简称为时域，将方位时域距离频域简称为时频域，将方位频域距离频域简称为频域。若无特殊说明，各步骤均对复数据进行操作。
按照公式一，校正接收机频率特性误差：
REi(tm，fr)＝Ri(tm，fr)Ei(fr)    (公式一)
上式中REi(tm，fr)表示通道i接收机频率特性误差校正后的时频域回波(简称校正后时频域回波)，Ei(fr)表示通道i的均衡器。
第二步，基于校正后时频域回波的距离谱均衡和SAR成像。
对校正后时频域回波REi(tm，fr)进行距离谱均衡。距离谱均衡是对REi(tm，fr)的多条距离线幅度在时频域进行加权平均。对距离谱均衡后的REi(tm，fr)采用低频SAR成像算法，获得时域图像RSi(tm，tf)。
第三步，基于时域图像RSi(tm，tf)的配准。
对第二步获得的时域图像RSi(tm，tf)进行配准，获得配准后的时域图像Si(tm，tf)和时频域干涉图IFATI(tm，fr)。本步骤去除的配准误差包括通道间的时间误差和相位误差，具体实施过程如下：
步骤3.1：选取感兴趣区域。
对时域图像RSi(tm，tf)进行幅度归一化处理，以整幅图像中幅度最大值为归一化处理的参考值，然后选取感兴趣区域SRi(tm，tf)。SR1(tm，tf)和SR2(tm，tf)对应了不同通道的相同坐标区域。感兴趣区域SRi(tm，tf)中至少包含一个强目标点；强目标点是场景中反射较强的静止地物，其归一化幅度值大于-20dB。
步骤3.2：估计通道间的时间误差。
对每个通道进行下述处理：对感兴趣区域SRi(tm，tf)进行线性插值，检测线性插值后SRi(tm，tf)的幅度最大值，记录幅度最大值的坐标(Tmi，Tfi)，Tmi表示方位向慢时间坐标，Tfi表示距离向快时间坐标。
按照公式二计算通道间的时间误差(ΔTm，ΔTf)：
(ΔTm，ΔTf)＝(Tm2-Tm1，Tf2-Tf1)    (公式二)
上式中ΔTm代表方位向慢时间差，ΔTf代表距离向快时间差。
步骤3.3：估计通道间相位误差。
按照公式三，获得感兴趣区域的时域干涉图IFRATI(tm，tf)。
IFR ATI ( t m , t f ) = SR 1 ′ ( t m , t f ) · SR 2 ′ * ( t m , t f )
SR′1(tm，tf)＝SR1(tm，tf)(公式三)
SR′2(tm，tf)＝SR2(tm-ΔTm，tf-ΔTf)
上式中(·)*表示共轭操作。
按照公式四，获得感兴趣区域的时域干涉图IFRATI(tm，tf)上强目标点主瓣中心的干涉相位φ0。
φ0＝∠IFRATI(Tm，Tf)    (公式四)
上式中∠表示取复图像像素点的相位角度，(Tm，Tf)表示强目标主瓣中心所在的方位向慢时间坐标、距离向快时间坐标。
步骤3.4：判断循环是否结束。
重复步骤3.1～3.3，每次重复过程中选取的感兴趣区域各不相同，对重复获得的干涉相位φ0取均值，获得通道间相位误差
Figure BDA0000055496560000041
重复次数至少3次。
步骤3.5：获得配准后的时域图像和时频域干涉图。
按照公式五，获得配准后的时域图像Si(tm，tf)。
S 1 ( t m , t f ) = RS 1 ( t m , t f ) S 2 = ( t m , t f ) = RS 2 ( t m - ΔTm , t f - ΔTf ) exp ( j φ 0 ‾ ) (公式五)
上式中exp(·)表示指数运算，j为虚数单位。对Si(tm，tf)进行傅立叶变换，得到配准后的时频域图像Si(tm，fr)。按照公式六，获得配准后的时频域干涉图IFATI(tm，fr)：
IF ATI ( t m , f r ) = S 1 ( t m , f r ) · S 2 * ( t m , f r ) (公式六)
第四步、基于时频域的RFI检测及抑制。
本步骤包括以下步骤：
步骤4.1：检测RFI信号。
按照公式七，计算时频域像素点幅度ρts(tm)。
ρ ts ( t m ) = E f r [ ρ ts ( t m , f r ) ] (公式七)
ρts(tm，fr)＝|S1(tm，fr)|
上式中
Figure BDA0000055496560000045
代表在方位向慢时间tm时刻对距离频率fr求期望，|·|表示对复图像像素点取幅度值。
利用假设检验检测配准后的时频域干涉图IFATI(tm，fr)像素点是否受到RFI影响。当假设检验的结果是H0时，IFATI(tm，fr)像素点未受RFI影响；当假设检验的结果是H1时，IFATI(tm，fr)像素点受到RFI影响。假设检验的公式如公式八所示：
(公式八)
上式中ρtr(tm，fr)表示IFATI(tm，fr)上待检测像素点干涉幅度的开方；
Figure BDA0000055496560000052
表示IFATI(tm，fr)上待检测像素点的干涉相位；γ0表示干涉图像素幅值比门限，
Figure BDA0000055496560000053
表示干涉图相位门限，通常选取γ0＝0.5，
Figure BDA0000055496560000054
步骤4.2：抑制RFI幅度。
当配准后的时频域干涉图IFATI(tm，fr)像素点受到RFI影响时，对时频域图像Si(tm，fr)进行中值滤波，获得抑制RFI幅度的时频域幅度图像
Figure BDA0000055496560000055
滤波操作如公式九所示：
D i % ( t m , f r ) = MedFilter [ | S i ( t m , f r ) | ] , 当fr＝fRFI且tm∈Tk时    (公式九)
上式中MedFilter[·]表示中值滤波，fRFI为RFI的距离频点，Tk为RFI的持续方位时间。该中值滤波器在每一个方位慢时间tm沿距离频域fr估计未受RFI影响的像素点幅度。中值滤波器的宽度大于检测到的RFI信号带宽的两倍；中值滤波仅替换时频域图像Si(tm，fr)中受RFI影响的像素点幅度，即步骤4.1中IFATI(tm，fr)上检测到RFI的像素点。
当配准后的时频域干涉图IFATI(tm，fr)像素点未受RFI影响时，按照公式十，获得抑制RFI幅度的时频域幅度图像
D i % ( t m , f r ) = | S i ( t m , f r ) | , 当fr≠fRFI或
Figure BDA0000055496560000059
时(公式十)
步骤4.3：抑制RFI相位。
按照公式十一，得到通道1抑制了RFI的时频域图像
Figure BDA00000554965600000510
D ^ 1 ( t m , f r ) = D 1 % ( t m , f r ) · exp ( j ∠ S 1 ( t m , f r ) ) (公式十一)
当配准后的时频域干涉图IFATI(tm，fr)像素点受到RFI影响时，对IFATI(tm，fr)在每一个方位慢时间tm沿距离频域fr进行中值滤波，估计受RFI影响的像素点干涉相位
Figure BDA00000554965600000512
如公式十二所示：
Δ θ ^ ( t m , f r ) = MedFilter [ ∠ IF ATI ( t m , f r ) ] 当fr＝fRFI且tm∈Tk时     (公式十二)
中值滤波的带宽大于RFI干涉相位所占据带宽的两倍。按照公式十三，得到通道2抑制了RFI的时频域图像
Figure BDA0000055496560000062
(公式十三)
当fr＝fRFI且tm∈Tk时
当配准后的时频域干涉图IFATI(tm，fr)像素点未受到RFI影响时，按照公式十四，得到通道2抑制了RFI的时频域图像
Figure BDA0000055496560000063
D ^ 2 ( t m , f r ) = D 2 % ( t m , f r ) · exp ( j ∠ S 2 ( t m , f r ) ) 当fr≠fRFI或 t m ∉ T k 时(公式十四)
第五步、残留误差校正。
本步骤包括以下步骤：
步骤5.1：残留幅度误差校正。
对
Figure BDA0000055496560000066
进行傅立叶变换得到抑制了RFI的频域图像
Figure BDA0000055496560000067
fa表示方位频率。按照公式十五，计算频域幅度误差G(fa，fr)：
G ( f a , f r ) = | D ^ 1 ( f a , f r ) / D ^ 2 ( f a , f r ) | (公式十五)
上式中通道幅度误差G(fa，fr)包括通道间的天线方向图幅度误差。
按照公式十六，获得去除了残留幅度误差的频域图像
Figure BDA0000055496560000069
D ^ 1 ′ ( f a , f r ) = D ^ 1 ( f a , f r ) D ^ 2 ′ ( f a , f r ) = D ^ 1 ( f a , f r ) G ( f a , f r ) (公式十六)
步骤5.2：残留相位误差校正。
将步骤5.1获得的进行傅立叶逆变换，得到去除了残留幅度误差的时域图像
Figure BDA00000554965600000612
按照公式十七，获得去除了残留幅度误差的时域干涉图
Figure BDA00000554965600000613
IFD ^ ATI ( t m , t f ) = D ^ 1 ′ ( t m , t f ) · D ^ 2 ′ * ( t m , t f ) (公式十七)
按照公式十八，对
Figure BDA00000554965600000615
的干涉相位进行中值滤波，获得残留的干涉相位误差
Figure BDA00000554965600000616
(公式十八)
公式十八的中值滤波器是方位向慢时间、距离向快时间的二维滤波器，滤波器的二维宽度为场景均匀杂波区域的平均大小。残留的干涉相位误差包括天线方向图相位误差、垂直基线误差。
按照公式十九，校正
Figure BDA0000055496560000071
残留的干涉相位误差
Figure BDA0000055496560000072
获得误差校正后的时域干涉图
Figure BDA0000055496560000073
Figure BDA0000055496560000074
(公式十九)
对经过上述步骤获得的误差校正后的时域干涉图
Figure BDA0000055496560000075
进行后续处理，如进行干涉幅度相位的联合检测，能精确检测运动目标，获取运动目标的运动信息。
采用本发明可以取得以下技术效果：
1.本发明利用时频域的干涉图进行RFI检测，避免了构造自适应滤波器的复杂计算，具有简单易行、准确率高、鲁棒性强、易于实现的优点；利用RFI在时频域图像上的分布特点，使假设检验门限随时频域像素点幅度变化，能够有效检测具有窄带性、时变性和空变性的RFI。
2.本发明利用时频域中值滤波获得未受RFI影响的时频域像素点幅度与干涉相位，有效抑制了RFI的幅度和相位，能够保留时域干涉图上场景真实的幅度与相位信息，满足后续处理的需求。
3.本发明提供的RFI检测及抑制方法，充分利用低频SAR信号与干涉图特点，校正了通道间各类误差。一是利用基于回波的接收机频率特性误差校正步骤校正了接收机频率特性误差；二是利基于时域图像的配准步骤校正了配准误差；三是利用残留误差校正步骤校正了天线方向图误差和垂直基线误差。上述每个步骤的存储量小，计算量低，能满足实际实验需求，易于工程实现。
附图说明
图1为本发明的原理流程示意图；
图2为本发明的数据流程示意图；
图3为仿真数据的时域图像；
图4为仿真数据时域干涉图的干涉相位对比图；
图5为实测数据时域干涉图的干涉相位对比图。
具体实施方式
图1为本发明的原理流程示意图。该流程有利于提高双通道低频SAR沿航向干涉方法的性能。总流程一共包括五个步骤：
第一步，基于回波的接收机频率特性误差校正。对各通道回波校正接收机频率特性误差，获得校正后时频域回波。其中，本步骤使用的均衡器Ei(f)通过雷达工作前的闭环测试获得，可参见F.E.Churchill，W.G.Ogar，B.J.Thompson编著的“The Correction of I and Q Errors in a Coherent Processor”，IEEETransactions on Aerospace and Electronic Systems，Vol.17，No.1，1981，pp.131-137和吕孝雷，苏军海，邢孟道，张守宏编著的“三通道SAR-GMTI误差校正方法的研究”，系统工程与电子技术，Vol.30，No.62008年6月，pp.1037-1042。
第二步，基于校正后时频域回波的距离谱均衡抑制和SAR成像。对各通道校正后时频域回波进行距离谱均衡，获得距离谱均衡后的回波；对该回波进行低频SAR成像，获得各通道的时域图像。距离谱均衡能够初步抑制RFI，可参见董臻，梁甸农，黄晓涛编著的“VHF/UHF UWB SAR基于通道均衡的RFI抑制方法”，电子与信息学报，Vol.30，No.3，2008，pp.550-553。
第三步，基于时域图像的配准。本步骤包括：步骤3.1，选取感兴趣区域；步骤3.2，估计通道间时间误差；步骤3.3，估计通道间相位误差；步骤3.4，判断循环是否结束；步骤3.5，获得配准后的时域图像和时频域干涉图。
第四步，基于时频域的RFI检测及抑制。本步骤包括：步骤4.1，检测RFI信号；步骤4.2，抑制RFI幅度；步骤4.3，抑制RFI相位。本步骤抑制了RFI幅度和相位，能够恢复未受RFI影响的时频域像素点的幅度与干涉相位。
第五步，残留误差校正。本步骤包括：步骤5.1，残留幅度误差校正；步骤5.2，残留相位误差校正。
图2为本发明的数据流程示意图，每个步骤的数据流向如下：
第一步中，每个通道接收到的时域回波Ri(tm，tf)，经过傅立叶变换，得到时频域回波Ri(tm，fr)。再对Ri(tm，fr)校正接收机频率特性误差，获得校正后时频域回波REi(tm，fr)。
第二步中，对每个通道校正后时频域回波REi(tm，fr)进行距离谱均衡，得到距离谱均衡后的REi(tm，fr)。再对其进行SAR成像，获得时域图像RSi(tm，tf)。
第三步中，步骤3.1从时域图像RSi(tm，tf)获得感兴趣区域SRi(tm，tf)；步骤3.2从SRi(tm，tf)估计通道间时间误差(ΔTm，ΔTf)；步骤3.3利用SRi(tm，tf)与(ΔTm，ΔTf)，获得感兴趣区域的时域干涉图IFRATI(tm，tf)，再从IFRATI(tm，tf)获得干涉相位φ0；步骤3.4经过多次重复3.1～3.3操作，获得通道间相位误差
Figure BDA0000055496560000091
步骤3.5利用估计得到的通道间时间误差(ΔTm，ΔTf)与通道间相位误差
Figure BDA0000055496560000092
获得配准后的时域图像Si(tm，tf)，对Si(tm，tf)傅立叶变换，得到配准后的时频域图像Si(tm，fr)。再进行干涉操作，获得配准后的时频域干涉图IFATI(tm，fr)。
第四步中，步骤4.1在IFATI(tm，fr)上检测RFI信号，步骤4.2抑制RFI幅度，得到抑制RFI幅度的时频域幅度图像
Figure BDA0000055496560000093
步骤4.3进一步抑制RFI相位，获得抑制了RFI的时频域图像
Figure BDA0000055496560000094
第五步中，步骤5.1先将
Figure BDA0000055496560000095
进行傅立叶变换，得到抑制了RFI的频域图像
Figure BDA0000055496560000096
然后通过计算残留幅度误差，得到去除了残留幅度误差的频域图像
Figure BDA0000055496560000097
步骤5.2先将
Figure BDA0000055496560000098
进行傅立叶逆变换，得到去除了残留幅度误差的时域图像
Figure BDA0000055496560000099
然后对
Figure BDA00000554965600000910
进行干涉获得去除了残留幅度误差的时域干涉图
Figure BDA0000055496560000101
再进行基于
Figure BDA0000055496560000102
干涉相位的中值滤波，并校正
Figure BDA0000055496560000103
残留的干涉相位误差
Figure BDA0000055496560000104
最终得到误差校正后的时域干涉图
Figure BDA0000055496560000105
图3～图5给出了利用本发明的具体实施方式在实验室对仿真数据进行仿真实验的结果。由于仿真数据难以准确模拟实际系统的残留误差，图6给出了利用本发明的具体实施方式在实验室对实测数据进行处理的结果。仿真实验的仿真数据包括：未加入RFI的时域回波形成的仿真数据集A，将仿真数据集A加入RFI形成的仿真数据集B。
图3为仿真数据的时域图像。仿真数据集进行本发明具体实施方式第一步和第二步的处理，得到仿真数据集的时域图像。其中(a)为基于仿真数据集A得到的通道1的时域图像RS1(tm，tf)，(b)为基于仿真数据集B得到的通道1的时域图像RS1(tm，tf)。如图3中编号所示，时域图像中目标1～9为静止目标，目标10为运动目标，垂直于航向的速度vy为3m/s。各个目标的归一化幅度为：目标1、2、3、4、5为0dB，目标6、7、10为-10dB，目标8、9为-20dB。对比(a)与(b)说明，加入RFI前后，通道1的时域图像没有明显变化，该相对较弱的RFI不会影响成像质量，不能依赖成像质量进行抑制。
图4为仿真数据时域干涉图的干涉相位对比图。图中，3幅图的横坐标均为方位向距离坐标，纵坐标均为距离向距离坐标，每幅图右边的尺度条表示相位值。其中(a)是将仿真数据集A进行本发明具体实施方式第一步至第三步处理得到的时域图像Si(tm，tf)直接进行干涉，获得时域干涉图的干涉相位图。(b)是将仿真数据集B进行本发明具体实施方式第一步至第三步的处理得到的时域图像Si(tm，tf)直接进行干涉，获得时域干涉图的干涉相位图。在(a)上，静止目标的干涉相位基本为0°，因此目标1～9基本无法从图中看出；运动目标的干涉相位不为0°，其干涉相位值与目标垂直于航向的速度有关，因此目标10可以从图中明显看出。在(b)上静止强目标点1～5的主瓣的干涉相位基本为0°，其它相对较弱的静止目标6～9、运动目标10以及各个目标的旁瓣均受到了RFI的相位干扰，静止目标的干涉相位不为0°，运动目标的干涉相位偏离了理论值。对比(a)与(b)说明，RFI对强目标点的相位干扰较小，基本能保持在0°附近，而弱目标点以及目标旁瓣附近受到的相位干扰非常明显。从中可得出结论，当RFI信号较弱时，RFI对成像质量影响很小，但对干涉相位影响较大。进一步说明本发明利用时频域干涉图进行RFI抑制的有效性。(c)是将仿真数据集B利用本发明具体实施方式进行处理得到的误差校正后的时域干涉图的干涉相位图，对比(b)说明，本发明充分抑制了RFI，静止和运动目标的干涉相位均得到了还原，取得了较好的效果；对比(a)说明，本发明的抑制效果与未受RFI影响的情况基本一致。
图5为实测数据时域干涉图的干涉相位对比图。该实测数据中存在已知的运动目标和RFI。其中(a)为未充分抑制RFI获得的时域干涉图的干涉相位图。该时域干涉图是将实测数据进行本发明具体实施方式第一步至第三步处理得到的时域图像Si(tm，tf)直接进行干涉获得的。(b)为采用本发明抑制RFI后获得的误差校正后的时域干涉图
Figure BDA0000055496560000111
的干涉相位图。(a)说明，未充分抑制RFI时，静止目标受RFI的影响严重，干涉相位变化剧烈，改变了时域干涉图幅度相位的联合分布，无法利用干涉图进行后续的运动目标检测。(b)体现了经过本发明的处理，静止目标的干涉相位在0°附近，由运动目标干涉相位推出的垂直于航向的速度与已知运动参数相符。实测数据的应用结果进一步说明了本发明的有效性。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (2)

1.一种低频SAR（Synthetic Aperture Radar，合成孔径雷达）的RFI（RadioFrequency Interference，射频干扰）抑制及误差校正方法，其特征在于，包括下述步骤：
第一步，基于回波的接收机频率特性误差校正；
将通道i（i＝1,2）接收的方位时域距离时域回波Ri(tm,tf)进行傅立叶变换，得到通道i方位时域距离频域回波Ri(tm,fr)；其中tm表示方位向慢时间，tf表示距离向快时间，fr表示距离频率；为了叙述的简明，后文将方位时域距离时域简称为时域，将方位时域距离频域简称为时频域，将方位频域距离频域简称为频域；
按照公式一，校正接收机频率特性误差：
REi(tm,fr)＝Ri(tm,fr)Ei(fr)（公式一）
上式中REi(tm,fr)表示通道i接收机频率特性误差校正后的时频域回波（简称校正后时频域回波），Ei(fr)表示通道i的均衡器；
第二步，基于校正后时频域回波的距离谱均衡和SAR成像；
对校正后时频域回波REi(tm,fr)进行距离谱均衡；距离谱均衡是对REi(tm,fr)的多条距离线幅度在时频域进行加权平均；对距离谱均衡后的REi(tm,fr)采用低频SAR成像算法，获得时域图像RSi(tm,tf)；
第三步，基于时域图像RSi(tm,tf)的配准；
步骤3.1：选取感兴趣区域；
对时域图像RSi(tm,tf)进行幅度归一化处理，以整幅图像中幅度最大值为归一化处理的参考值，然后选取感兴趣区域SRi(tm,tf)；SR1(tm,tf)和SR2(tm,tf)对应了不同通道的相同坐标区域；感兴趣区域SRi(tm,tf)中至少包含一个强目标点；强目标点是场景中反射较强的静止地物，其归一化幅度值大于-20dB；
步骤3.2：估计通道间的时间误差；
对每个通道进行下述处理：对感兴趣区域SRi(tm,tf)进行线性插值，检测线性插值后SRi(tm,tf)的幅度最大值，记录幅度最大值的坐标(Tmi,Tfi)，Tmi表示方位向慢时间坐标，Tfi表示距离向快时间坐标；
按照公式二计算通道间的时间误差(ΔTm,ΔTf）：
(ΔTm,ΔTf)＝(Tm2-Tm1,Tf2-Tf1)（公式二）
上式中ΔTm代表方位向慢时间差，ΔTf代表距离向快时间差；
步骤3.3：估计通道间相位误差；
按照公式三，获得感兴趣区域的时域干涉图IFRATI(tm,tf)；
IFR ATI ( t m , t f ) = SR 1 ′ ( t m , t f ) · SR 2 ′ * ( t m , t f ) （公式三）
SR 1 ′ ( t m , t f ) = SR 1 ( t m , t f )
SR 2 ′ ( t m , t f ) = SR 2 ( t m - ΔTm , t f - ΔTf )
上式中(·)*表示共轭操作；
按照公式四，获得感兴趣区域的时域干涉图IFRATI(tm,tf)上强目标点主瓣中心的干涉相位φ0；
φ0＝∠IFRATI(Tm,Tf)（公式四）
上式中∠表示取复图像像素点的相位角度，(Tm,Tf)表示强目标主瓣中心所在的方位向慢时间坐标、距离向快时间坐标；
步骤3.4：判断循环是否结束；
重复步骤3.1～3.3，每次重复过程中选取的感兴趣区域各不相同，对重复获得的干涉相位φ0取均值，获得通道间相位误差
Figure FDA00001860597400024
重复次数至少3次；
步骤3.5：获得配准后的时域图像和时频域干涉图；
按照公式五，获得配准后的时域图像Si(tm,tf)；
S 1 ( t m , t f ) = RS 1 ( t m , t f ) S 2 ( t m , t f ) = RS 2 ( t m - ΔTm , t f - ΔTf ) exp ( j φ ‾ 0 ) （公式五）
上式中exp(·)表示指数运算，j为虚数单位；对Si(tm,tf)进行傅立叶变换，得到配准后的时频域图像Si(tm,fr)；按照公式六，获得配准后的时频域干涉图IFATI(tm,fr)：
IF ATI ( t m , f r ) = S 1 ( t m , f r ) · S 2 * ( t m , f r ) （公式六）
第四步、基于时频域的RFI检测及抑制；
步骤4.1：检测RFI信号；
按照公式七，计算时频域像素点幅度ρts(tm)；
ρ ts ( t m ) = E f r [ ρ ts ( t m , f r ) ] （公式七）
ρts(tm,fr)＝|S1(tm,fr)|
上式中
Figure FDA00001860597400033
代表在方位向慢时间tm时刻对距离频率fr求期望，|·|表示对复图像像素点取幅度值；
利用假设检验检测配准后的时频域干涉图IFATI(tm,fr)像素点是否受到RFI影响；当假设检验的结果是H0时，IFATI(tm,fr)像素点未受RFI影响；当假设检验的结果是H1时，IFATI(tm,fr)像素点受到RFI影响；假设检验的公式如公式八所示：
Figure FDA00001860597400034
（公式八）
上式中ρtr(tm,fr)表示IFATI(tm,fr)上待检测像素点干涉幅度的开方；
Figure FDA00001860597400035
表示IFATI(tm,fr)上待检测像素点的干涉相位；γ0表示干涉图像素幅值比门限，
Figure FDA00001860597400036
表示干涉图相位门限；
步骤4.2：抑制RFI幅度；
当配准后的时频域干涉图IFATI(tm,fr)像素点受到RFI影响时，对时频域图像Si(tm,fr)进行中值滤波，获得抑制RFI幅度的时频域幅度图像
Figure FDA00001860597400037
滤波操作如公式九所示：
D ~ i ( t m , f r ) = MedFilter [ | S i ( t m , f r ) | ] , 当fr＝fRFI且tm∈Tk时（公式九）
上式中MedFilter[·]表示中值滤波，fRFI为RFI的距离频点，Tk为RFI的持续方位时间；该中值滤波器在每一个方位慢时间tm沿距离频域fr估计未受RFI影响的像素点幅度；中值滤波器的宽度大于检测到的RFI信号带宽的两倍；中值滤波仅替换时频域图像Si(tm,fr)中受RFI影响的像素点幅度，即步骤4.1中IFATI(tm,fr)上检测到的受到RFI影响的像素点；
当配准后的时频域干涉图IFATI(tm,fr)像素点未受RFI影响时，按照公式十，获得抑制RFI幅度的时频域幅度图像
D ~ i ( t m , f r ) = | S i ( t m , f r ) | , 当fr≠fRFI或 t m ∉ T k 时（公式十）
步骤4.3：抑制RFI相位；
按照公式十一，得到通道1抑制了RFI的时频域图像
Figure FDA00001860597400044
D ^ 1 ( t m , f r ) = D ~ 1 ( t m , f r ) · exp ( j ∠ S 1 ( t m , f r ) ) （公式十一）
当配准后的时频域干涉图IFATI(tm,fr)像素点受到RFI影响时，对IFATI(tm,fr)在每一个方位慢时间tm沿距离频域fr进行中值滤波，估计受RFI影响的像素点干涉相位如公式十二所示：
Δ θ ^ ( t m , f r ) = MedFilter [ ∠ IF ATI ( t m , f r ) ] , 当fr＝fRFI且tm∈Tk时（公式十二）中值滤波的带宽大于RFI干涉相位所占据带宽的两倍；按照公式十三，得到通道2抑制了RFI的时频域图像
Figure FDA00001860597400048
Figure FDA00001860597400049
当fr＝fRFI且tm∈Tk时（公式十三）
当配准后的时频域干涉图IFATI(tm,fr)像素点未受到RFI影响时，按照公式十四，得到通道2抑制了RFI的时频域图像
Figure FDA000018605974000410
D ^ 2 ( t m , f r ) = D ~ 2 ( t m , f r ) · exp ( j ∠ S 2 ( t m , f r ) ) , 当fr≠fRFI或 t m ∉ T k 时（公式十四）
第五步、残留误差校正；
步骤5.1：残留幅度误差校正；
对
Figure FDA000018605974000413
进行傅立叶变换得到抑制了RFI的频域图像
Figure FDA000018605974000414
fa表示方位频率；按照公式十五，计算频域幅度误差G(fa,fr)：
G ( f a , f r ) = | D ^ 1 ( f a , f r ) / D ^ 2 ( f a , f r ) | （公式十五）
上式中通道幅度误差G(fa,fr)包括通道间的天线方向图幅度误差；
按照公式十六，获得去除了残留幅度误差的频域图像
Figure FDA00001860597400051
D ^ 1 ′ ( f a , f r ) = D ^ 1 ( f a , f r ) D ^ 2 ′ ( f a , f r ) = D ^ 1 ( f a , f r ) G ( f a , f r ) （公式十六）
步骤5.2：残留相位误差校正；
将步骤5.1获得的进行傅立叶逆变换，得到去除了残留幅度误差的时域图像
Figure FDA00001860597400054
按照公式十七，获得去除了残留幅度误差的时域干涉图
I F ^ D ATI ( t m , t f ) = D ^ 1 ′ ( t m , t f ) · D ^ 2 ′ * ( t m , t f ) （公式十七）
按照公式十八，对
Figure FDA00001860597400057
的干涉相位进行中值滤波，获得残留的干涉相位误差
∠ I F ~ M ATI ( t m , t f ) = MedFilter ( ∠ I F ^ D ATI ( t m , t f ) ) （公式十八）
公式十八的中值滤波器是方位向慢时间、距离向快时间的二维滤波器，滤波器的二维宽度为场景均匀杂波区域的平均大小；残留的干涉相位误差包括天线方向图相位误差、垂直基线误差；
按照公式十九，校正
Figure FDA000018605974000510
残留的干涉相位误差获得误差校正后的时域干涉图
Figure FDA000018605974000512
I F ^ M ATI ( t m , t f ) = I F ^ D ATU ( t m , t f ) · exp ( - j ∠ I F ~ M ATI ( t m , t f ) ) （公式十九）。
2.根据权利要求1所述的低频SAR的RFI抑制及误差校正方法，其特征在于，干涉图像素幅值比门限γ0＝0.5，干涉图相位门限
Figure FDA000018605974000514
CN201110094916A 2011-04-15 2011-04-15 低频合成孔径雷达射频干扰抑制及误差校正方法 Expired - Fee Related CN102243300B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110094916A CN102243300B (zh) 2011-04-15 2011-04-15 低频合成孔径雷达射频干扰抑制及误差校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110094916A CN102243300B (zh) 2011-04-15 2011-04-15 低频合成孔径雷达射频干扰抑制及误差校正方法

Publications (2)

Publication Number Publication Date
CN102243300A CN102243300A (zh) 2011-11-16
CN102243300B true CN102243300B (zh) 2012-10-10

Family

ID=44961457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110094916A Expired - Fee Related CN102243300B (zh) 2011-04-15 2011-04-15 低频合成孔径雷达射频干扰抑制及误差校正方法

Country Status (1)

Country Link
CN (1) CN102243300B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103176184B (zh) * 2013-02-06 2014-02-19 中国科学院电子学研究所 一种结合干扰抑制方法的p波段sar成像处理方法
CN103616669B (zh) * 2013-12-03 2016-03-02 西安电子科技大学 二维频域中对非均匀场景进行通道均衡方法
CN103675768B (zh) * 2013-12-17 2015-10-21 中国人民解放军国防科学技术大学 一种星载sar射频干扰抑制方法
CN104297736B (zh) * 2014-09-23 2017-01-11 西安电子工程研究所 一种用于雷达接收通道阵列的校正方法
CN105974376B (zh) * 2016-05-09 2018-05-25 中国人民解放军国防科学技术大学 一种sar射频干扰抑制方法
CN106597443B (zh) * 2016-12-20 2019-02-26 中国人民解放军国防科学技术大学 低频合成孔径雷达图像对比相消射频干扰抑制方法
CN106970374B (zh) * 2017-05-27 2019-05-07 中国人民解放军63892部队 干扰条件下单脉冲雷达目标和干扰源角度的测量方法
CN107991676B (zh) * 2017-12-01 2019-09-13 中国人民解放军国防科技大学 星载单航过InSAR系统对流层误差校正方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6166678A (en) * 1999-09-07 2000-12-26 The United States Of America As Represented By The Secretary Of The Air Force Fourier-transform-based adaptive radio interference mitigation
CN100472223C (zh) * 2006-01-20 2009-03-25 武汉大学 高频雷达抗射频干扰方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘宝泉等.干涉合成孔径雷达复图像配准解析搜索方法.《数据采集与处理》.2008,第23卷(第6期),646-651.
周红等.基于子孔径图像的低频SAR多通道均衡方法.《国防科技大学学报》.2010,第32卷(第4期),78-83.
基于子孔径图像的低频SAR多通道均衡方法;周红等;《国防科技大学学报》;20100831;第32卷(第4期);78-83 *
干涉合成孔径雷达复图像配准解析搜索方法;刘宝泉等;《数据采集与处理》;20081130;第23卷(第6期);646-651 *

Also Published As

Publication number Publication date
CN102243300A (zh) 2011-11-16

Similar Documents

Publication Publication Date Title
CN102243300B (zh) 低频合成孔径雷达射频干扰抑制及误差校正方法
Jiang et al. Hybrid approach for unbiased coherence estimation for multitemporal InSAR
Zheng et al. ISAR imaging of targets with complex motions based on a noise-resistant parameter estimation algorithm without nonuniform axis
CN106526591A (zh) 机动目标高分辨isar子孔径融合成像方法
CN102608586B (zh) 一种基于神经网络的合成孔径雷达射频干扰抑制方法
CN103616686B (zh) 一种基于混合模式的全极化干涉合成孔径雷达影像的最优相位估计方法
CN103116162B (zh) 基于目标空间稀疏性的高分辨声呐定位方法
CN103901442B (zh) 一种单频多路径效应探测与缓解方法
CN103048655B (zh) 天波超视距雷达频域超分辨微多径测高方法
CN104730505A (zh) 一种多通道sar地面目标检测定位方法和装置
CN103176168A (zh) 一种机载非正侧视阵雷达近程杂波对消方法
CN104076404A (zh) 运用多通道相干抑制地磁背景噪声的磁异常探测方法
CN104502895A (zh) 一种稳健的基于多通道sar特征投影的杂波抑制方法及装置
Shen et al. Site-specific real-time GPS multipath mitigation based on coordinate time series window matching
Xu et al. Maneuvering target imaging and scaling by using sparse inverse synthetic aperture
CN103760534A (zh) 一种星载sar数据的电离层色散效应校正方法
CN105699947B (zh) 一种sar图像旁瓣抑制方法
CN103809180A (zh) 用于InSAR地形测量的方位向预滤波处理方法
Li et al. An adaptive phase optimization algorithm for distributed scatterer phase history retrieval
CN103245949B (zh) 一种基于改进理想滤波器的sar方位模糊抑制方法
Tan et al. Efficient FFT based multi source DOA estimation for ULA
CN103760540A (zh) 基于重构信号和1-范数的动目标检测与参数估计方法
CN110988873B (zh) 基于能量中心提取的单通道sar舰船速度估计方法及系统
CN108375779A (zh) 一种基于Kalman滤波器的导航接收机干扰抑制方法
CN105116410B (zh) 基于线性模型匹配的干涉相位图自适应滤波算法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20121010

Termination date: 20200415

CF01 Termination of patent right due to non-payment of annual fee