CN106842196B - 一种用于sar子孔径成像的方位模糊抑制方法 - Google Patents

一种用于sar子孔径成像的方位模糊抑制方法 Download PDF

Info

Publication number
CN106842196B
CN106842196B CN201510882971.8A CN201510882971A CN106842196B CN 106842196 B CN106842196 B CN 106842196B CN 201510882971 A CN201510882971 A CN 201510882971A CN 106842196 B CN106842196 B CN 106842196B
Authority
CN
China
Prior art keywords
sub
azimuth
distance
aperture
processing
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
CN201510882971.8A
Other languages
English (en)
Other versions
CN106842196A (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.)
Leihua Electronic Technology Research Institute Aviation Industry Corp of China
Original Assignee
Leihua Electronic Technology Research Institute Aviation Industry Corp of China
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 Leihua Electronic Technology Research Institute Aviation Industry Corp of China filed Critical Leihua Electronic Technology Research Institute Aviation Industry Corp of China
Priority to CN201510882971.8A priority Critical patent/CN106842196B/zh
Publication of CN106842196A publication Critical patent/CN106842196A/zh
Application granted granted Critical
Publication of CN106842196B publication Critical patent/CN106842196B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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
    • 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/9004—SAR image acquisition techniques
    • G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth

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)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种用于SAR子孔径成像的方位模糊抑制方法，其特征在于，包括以下步骤：步骤1：按照成像处理流程对原始回波数据进行距离向处理然后再沿方位向划分子孔径，其中对原始回波数据进行距离向处理具体为：进行距离向FFT；再进行距离压缩；进行距离向IFFT；去除脉冲持续时间点数；步骤2：数据处理完成后进行子孔径距离向处理，子孔径距离向处理的具体步骤为首先子孔径方位FFT；再进行数据两端进行补零，增加两端脉冲持续时间点数来实现扩展子孔径长度；进行距离向FFT；进行距离向逆压缩、三次相位滤波，参考点相位补偿；进行距离向IFFT；进行Chirp Scaling运算。

Description

一种用于SAR子孔径成像的方位模糊抑制方法
技术领域
本发明涉及一种用于SAR子孔径成像的方位模糊抑制方法，属于SAR成像技术领域。
背景技术
在子孔径成像处理中，方位信号的频谱主要表现为二次相位项，包括决定方位分辨率的多普勒信号，与距离迁徙项相关的二次项。方位模糊是由与距离迁徙项相关的二次项引起的。距离迁徙校正(RCMC)的参考函数沿方位向表现为线性调频信号，并且带宽小于PRF，在方位频域，信号不满足采样定理的要求，校正参考函数将会出现混叠，这将导致方位模糊的出现。方位模糊位置与目标真实位置之间的距离为子孔径长度。
关于方位模糊的抑制，存在两种方法，一种是消除距离迁徙校正参考函数的混叠，这种方法可以通过增加子孔径信号的长度来完成；由于在全孔径处理时，距离迁徙校正参考函数带宽仍然小于多普勒带宽，因此这种方法无法达到对方位模糊完全抑制的目的。另一种是采用子孔径重叠的方法进行方位模糊抑制，然而这种方法会使得成像处理的存储量和运算量有较大的增加。
发明内容
本发明所要解决的技术问题是针对上述背景技术的不足，提供了在RCMC之前进行子孔径信号长度扩展(两端补零)抑制方位模糊的方法；当子孔径长度扩展为原长度的两倍时，方位模糊的影响基本上被抑制了。
技术方案
本发明提供了一种用于SAR子孔径成像的方位模糊抑制方法，其特征在于，包括以下步骤：
步骤1：按照成像处理流程对原始回波数据进行距离向处理然后再沿方位向划分子孔径，对原始回波数据进行距离向处理具体为：
进行距离向FFT；
再进行距离压缩；
进行距离向IFFT；
去除脉冲持续时间点数；
步骤2：数据处理完成后进行子孔径距离向处理，子孔径距离向处理的具体步骤为
首先子孔径方位FFT；
再进行数据两端进行补零，增加两端脉冲持续时间点数；
进行距离向FFT；
进行距离向逆压缩、三次相位滤波，参考点相位补偿；
进行距离向IFFT；
进行Chirp Scaling运算；
进行距离向FFT；
进行距离压缩、二次距离压缩、距离弯曲校正；
进行距离向IFFT；
进行剩余相位补偿；
子孔径方位向IFFT；
进行方位向合成方位子孔径数据；
步骤3：进行距离向划分子带回波；
进行方位向FFT；
进行方位向压缩；
进行方位向IFFT；
合并距离子带图像；
最后输出SAR图像。
进一步的，步骤2中的再进行数据两端进行补零，增加两端脉冲持续时间点数来实现扩展子孔径长度具体为：
子孔径距离向处理将距离迁徙校正前的子孔径回波变换到二维时域，沿方位向进行子孔径长度扩展，即将子孔径回波沿方位向两端补零，
设扩展后的子孔径长度为TsF，子孔径长度扩展后，距离迁徙校正参考函数可以表示为:
式(1)中,n为子孔径序号，Ns为沿方位向划分的子孔径个数，为N为子孔径方位点数，fr为距离向基带信号频率,为距离迁徙校正参考函数的带宽，rect[·]为矩形窗函数，Ka为目标方位向多普勒频率调频率,fc为中心频率，fa为方位多普勒频率。
进一步的，步骤2中的进行距离压缩、二次距离压缩、距离弯曲校正具体为：
进行子孔径距离迁徙校正:
将式(2)的结果变换到二维时域，这时回波信号可以表示为
进一步的，步骤2中的进行方位向合成方位子孔径数据具体为：
扩展前的子孔径长度为Ts，将扩展的子孔径长度去掉后，方位时间的取值范围为-Ts/2≤tm≤Ts/2，这时，方位模糊位置变为-TsF和TsF，距离迁徙校正之后，将子孔径扩展部分去除掉，方位模糊得到一定程度的抑制。
从式(3)中可以看出，当子孔径扩展长度满足TsF<2Ts时，方位信号依然存在混叠，这时仍然存在方位模糊；当子孔径扩展长度满足TsF＝2Ts时，方位向模糊基本上得到了抑制，但是方位模糊部分的旁瓣仍然会对子孔径成像结果有一定的影响，但这种影响基本上可以被忽略。当子孔径扩展长度更长时，方位模糊的影响可以完全被抑制掉。
发明创造的优点和用途
本发明方法简单实用，解决了子孔径处理算法中方位模糊抑制技术。与现有技术相比具有如下优点：
(1)与全孔径处理算法相比，本发明能够基本抑制方位模糊。
(2)与子孔径重叠算法相比，本发明在不需要增加大量的存储量和运算量的情况下，也可以基本抑制方位模糊。
附图说明
图1不同载频条件下方位压缩前沿方位向扩展的方位向剖面图；
图2方位向扩展对方位模糊的抑制结果；
图3改进子孔径成像算法流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,参照附图,对本发明进一步详细说明。
步骤1按照图3成像处理流程所示对原始回波数据进行距离向FFT，处理后的数据再进行距离压缩，进行距离向IFFT，去除脉冲持续时间点数的距离向处理方法；
最后对数据沿方位向划分子孔径；
步骤2对子孔径划分完成的数据子进行方位向FFT运算后，对数据两端进行补零，增加两端脉冲持续时间点数；对数据进行距离向FFT；进行距离向逆压缩、三次相位滤波，参考点相位补偿；进行距离向IFFT；进行Chirp Scaling运算；进行距离向FFT；进行距离压缩、二次距离压缩、距离弯曲校正；进行距离向IFFT；进行剩余相位补偿；子孔径方位向IFFT；进行方位向合成方位子孔径数据；
步骤3进行距离向划分子带回波；进行方位向FFT；进行方位向压缩；进行方位向IFFT；合并距离子带图像；最后输出SAR图像。
经过上述步骤，本发明方法处理得到的图像质量得到明显的改善，方位模糊得到基本抑制。

Claims (2)

1.一种用于SAR子孔径成像的方位模糊抑制方法，其特征在于，包括以下步骤：
步骤1：按照成像处理流程对原始回波数据进行距离向处理然后再沿方位向划分子孔径，其中对原始回波数据进行距离向处理具体为：
进行距离向FFT；
再进行距离压缩；
进行距离向IFFT；
去除脉冲持续时间点数；
步骤2：数据处理完成后进行子孔径距离向处理，子孔径距离向处理的具体步骤为
首先进行子孔径方位FFT；
再对数据两端进行补零，增加两端脉冲持续时间点数来实现扩展子孔径长度；
进行距离向FFT；
进行距离向逆压缩、三次相位滤波，参考点相位补偿；
进行距离向IFFT；
进行Chirp Scaling运算；
进行距离向FFT；
进行距离压缩、二次距离压缩、距离弯曲校正；
进行距离向IFFT；
进行剩余相位补偿；
进行子孔径方位向IFFT；
进行方位向的方位子孔径数据合成；
步骤3：进行距离向的子带回波划分；
进行方位向FFT；
进行方位向压缩；
进行方位向IFFT；
合并距离子带图像；
最后输出SAR图像，
其中，步骤2中的再进行数据两端进行补零，增加两端脉冲持续时间点数来实现扩展子孔径长度具体为：
子孔径距离向处理将距离迁徙校正前的子孔径回波变换到二维时域，沿方位向进行子孔径长度扩展，即将子孔径回波沿方位向两端补零，
设扩展后的子孔径长度为TsF，子孔径长度扩展后，距离迁徙校正参考函数可以表示为:
式(1)中,n为子孔径序号，Ns为沿方位向划分的子孔径个数，N为子孔径方位点数，fr为距离向基带信号频率,为距离迁徙校正参考函数的带宽，rect[·]为矩形窗函数，Ka为目标方位向多普勒频率调频率,fc为中心频率，fa为方位多普勒频率，
其中，步骤2中的进行距离压缩、二次距离压缩、距离弯曲校正具体为：
进行子孔径距离迁徙校正:
将式(2)的结果变换到二维时域，这时回波信号可以表示为
2.根据权利要求1所述的用于SAR子孔径成像的方位模糊抑制方法，其特征在于，步骤2中的进行方位向合成方位子孔径数据具体为：扩展前的子孔径长度为Ts，将扩展的子孔径长度去掉后，方位时间的取值范围为-Ts/2≤tm≤Ts/2，这时，方位模糊位置变为-TsF和TsF，距离迁徙校正之后，将子孔径扩展部分去除掉，方位模糊得到一定程度的抑制。
CN201510882971.8A 2015-12-03 2015-12-03 一种用于sar子孔径成像的方位模糊抑制方法 Active CN106842196B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510882971.8A CN106842196B (zh) 2015-12-03 2015-12-03 一种用于sar子孔径成像的方位模糊抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510882971.8A CN106842196B (zh) 2015-12-03 2015-12-03 一种用于sar子孔径成像的方位模糊抑制方法

Publications (2)

Publication Number Publication Date
CN106842196A CN106842196A (zh) 2017-06-13
CN106842196B true CN106842196B (zh) 2019-12-17

Family

ID=59149735

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510882971.8A Active CN106842196B (zh) 2015-12-03 2015-12-03 一种用于sar子孔径成像的方位模糊抑制方法

Country Status (1)

Country Link
CN (1) CN106842196B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106950567B (zh) * 2017-03-30 2019-05-14 中国人民解放军国防科学技术大学 基于高阶子孔径cs的超宽带滑聚sar成像处理方法
CN107703478B (zh) * 2017-10-27 2021-04-27 天津大学 基于互相关矩阵的扩展孔径二维doa估计方法
CN108896988B (zh) * 2018-05-10 2022-02-11 北京航空航天大学 一种针对mimo雷达近场测量的多波段超宽带成像处理方法
CN108896993B (zh) * 2018-05-10 2022-02-08 北京航空航天大学 一种复杂目标多波段雷达超宽带联合成像的分块处理方法
CN113721244A (zh) * 2021-07-08 2021-11-30 中国科学院空天信息创新研究院 合成孔径雷达的信噪比处理方法及装置、存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011169869A (ja) * 2010-02-22 2011-09-01 Mitsubishi Electric Corp レーダ信号処理装置
CN102819020A (zh) * 2012-08-17 2012-12-12 北京航空航天大学 一种俯冲模型方位向非线性变标的合成孔径雷达成像方法
CN102901964A (zh) * 2012-09-06 2013-01-30 内蒙古工业大学 二维多孔径ScanSAR成像方法
CN103197317A (zh) * 2013-03-26 2013-07-10 西安电子科技大学 基于fpga的sar成像方法
CN103472448A (zh) * 2013-07-15 2013-12-25 中国科学院电子学研究所 一种sar成像方法、装置及系统
CN104865571A (zh) * 2015-04-02 2015-08-26 西安电子科技大学 多通道多子带滑动聚束模式sar成像方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011169869A (ja) * 2010-02-22 2011-09-01 Mitsubishi Electric Corp レーダ信号処理装置
CN102819020A (zh) * 2012-08-17 2012-12-12 北京航空航天大学 一种俯冲模型方位向非线性变标的合成孔径雷达成像方法
CN102901964A (zh) * 2012-09-06 2013-01-30 内蒙古工业大学 二维多孔径ScanSAR成像方法
CN103197317A (zh) * 2013-03-26 2013-07-10 西安电子科技大学 基于fpga的sar成像方法
CN103472448A (zh) * 2013-07-15 2013-12-25 中国科学院电子学研究所 一种sar成像方法、装置及系统
CN104865571A (zh) * 2015-04-02 2015-08-26 西安电子科技大学 多通道多子带滑动聚束模式sar成像方法

Also Published As

Publication number Publication date
CN106842196A (zh) 2017-06-13

Similar Documents

Publication Publication Date Title
CN106842196B (zh) 一种用于sar子孔径成像的方位模糊抑制方法
CN101900813B (zh) 基于机动目标距离-瞬时调频的isar成像方法
CN105974414A (zh) 基于二维自聚焦的高分辨聚束sar自聚焦成像方法
CN103454632B (zh) 一站固定式调频连续波双基地sar成像方法
CN111781595B (zh) 基于匹配搜索和多普勒解模糊的复杂机动群目标成像方法
CN105445704B (zh) 一种sar图像中的雷达动目标抑制方法
CN109188432B (zh) 一种平行双基聚束sar快速bp成像方法
CN103323822A (zh) 一种估计通道误差的方法及装置
US20140232590A1 (en) Sar image formation
CN110456351B (zh) 基于时变幅值lfm信号参数估计的机动目标isar成像方法
CN106680785B (zh) 基于小波变换空间变迹的sar图像旁瓣抑制方法
CN102520403B (zh) 一种改进的基于频域频谱重建的频率步进sar成像方法
CN104698459A (zh) 一种应用于缺失数据的条带sar压缩感知成像方法
CN104020456A (zh) 一种基于多特显点线阵成像雷达系统幅相误差校正方法
CN103809180B (zh) 用于InSAR地形测量的方位向预滤波处理方法
CN109782277B (zh) 变pri的斜视聚束sar成像方法、装置、设备及存储介质
CN105022060A (zh) 针对快速空天目标的步进isar成像方法
CN104182942B (zh) Sar图像方位模糊抑制方法
JP6395672B2 (ja) レーダ装置
CN105469368A (zh) 一种干涉相位滤波方法
Liu et al. Novel range profile synthesis algorithm for linearly stepped-frequency modulated inversed synthetic aperture radar imaging of remote manoeuvring target
CN103308895A (zh) 一种基于重叠频谱的低阶相位误差估计的方法
CN111624602A (zh) 一种旋转FrFT域构建SAR方位向分辨率成像方法
CN103472448A (zh) 一种sar成像方法、装置及系统
CN113655480B (zh) 一种基于miaa谱估计的dvb-t信号无源isar成像方法

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