CN106772376B - 基于改进的rda进行合成孔径雷达图像目标旁瓣抑制的方法 - Google Patents

基于改进的rda进行合成孔径雷达图像目标旁瓣抑制的方法 Download PDF

Info

Publication number
CN106772376B
CN106772376B CN201710013702.7A CN201710013702A CN106772376B CN 106772376 B CN106772376 B CN 106772376B CN 201710013702 A CN201710013702 A CN 201710013702A CN 106772376 B CN106772376 B CN 106772376B
Authority
CN
China
Prior art keywords
orientation
distance
signal
point
moment
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
CN201710013702.7A
Other languages
English (en)
Other versions
CN106772376A (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 CN201710013702.7A priority Critical patent/CN106772376B/zh
Publication of CN106772376A publication Critical patent/CN106772376A/zh
Application granted granted Critical
Publication of CN106772376B publication Critical patent/CN106772376B/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
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details 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

Landscapes

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

Abstract

本发明公开了一种基于改进的RDA进行合成孔径雷达图像目标旁瓣抑制的方法，该方法是为了获得超分辨的SAR图像；具体处理步骤为：二维回波信号的方位向傅里叶变换→补余徒动校正→距离向傅里叶变换→一致徒动校正和距离压缩→距离向傅里叶逆变换→附加相位校正→方位向傅里叶逆变换→融合反演矩阵的方位压缩。相较于匹配滤波等常规的成像处理方式，本发明处理的成像结果能量更加集中，显著提高空间分辨率，达到超分辨的目的。相较于完全基于矩阵的SAR数据处理，本发明与高效率的传统成像算法相结合，因而显著提高处理效率。

Description

基于改进的RDA进行合成孔径雷达图像目标旁瓣抑制的方法
技术领域
本发明涉及合成孔径雷达的成像处理技术领域，更特别地说，是指一种基于改进的RDA进行合成孔径雷达图像目标旁瓣抑制的方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar，SAR)是一种全天时、全天候的主动式信息获取***，***结构如图1所示，合成孔径雷达的全链路***可以分为三个部分，即雷达载体平台上的回波数据获取部分、空间段数据传输部分以及地面段数据处理部分。合成孔径雷达***通过发射和接收脉冲信号，获取合成孔径雷达原始回波数据。
2012年6月电子工业出版社出版，Ian G.Cumming及Frank H.Wong编著的《合成孔径雷达成像——算法与实现》在第75页中指出，合成孔径雷达(Synthetic ApertureRadar，简称SAR)传统信号处理的核心思想是基于对SAR回波信号(雷达原始数据)进行距离和方位两个方向上的匹配滤波。距离多普勒算法(Range Doppler Algorithm，RDA)通过距离和方位上的频域操作，达到了高效的模块化处理要求，同时又具有了一维操作的简便性。该算法根据距离和方位上的大尺度时间差异，在两个一维操作之间使用距离徙动校正(RCMC)，对距离和方位进行了近似的分离处理，具体的RDA处理步骤如图1A所示。
2001年2月科学出版社出版，魏钟铨等编著的《合成孔径雷达卫星》在第202页中指出，峰值旁瓣比的大小决定了强目标“掩盖”若目标的能力；积分旁瓣比表征图像质量的重要指标之一，他定量的描述了一个局部较暗的区域被来自周围的明亮区域的能量泄露所“淹没”的程度。
《合成孔径雷达成像——算法与实现》中的58页写道，降低PSLR的一种方法是对频域匹配滤波器引入平滑窗，以减少主瓣到旁瓣的能量泄漏。窗是一个对信号频谱进行加权的对称实函数。权值在信号频谱中心处最大，向频谱两边逐渐衰落。典型的窗包括Taylor窗、Chebyshev窗、Hanning窗、Hamming窗和Kaiser窗。窗能够平滑频谱，即弱化频谱边缘处的不连续性。这样会降低压缩脉冲中的主瓣能量泄漏，但要以损失分辨率为代价。
2001年11月美国UNL大学环境遥测实验室的Xiaojian Xu等人发表的文章《RangeSidelobe Suppression Technique for Coherent Ultra Wide-Band Random NoiseRadar Imaging》提到，通过切趾滤波可以找到一种特殊的图像域滤波函数与SAR图像原始数据卷积就会得到一个比原来旁瓣电平低得多的输出响应，可以得到一幅高质量的SAR图像。但是，切趾滤波技术由于要通过多次迭代求其滤波器系数，因此运算量较大。
2007年8月美国Microwave and Radar Group的Carlos F.Castillo-Rubio等人发表的文章《Spatially Variant Apodization for Squinted Synthetic Aperture RadarImages》中提到，不规则形状支撑区分布导致点扩散函数产生非正交走向旁瓣，利用非正交走向旁瓣与正交走向旁瓣之间的差异信息来实现正交走向旁瓣抑制。
目前针对图像目标旁瓣的抑制处理方法为上文提到的频域加窗加权，切趾滤波器和谱变形这三种方法，但是这三种方法都有着各自的优缺点。频域加窗加权法由于要对频谱进行平滑处理，所以相应的缩小了信号带宽，使得信号的图像分辨率降低；切趾滤波器方法的不足在于需要先验条件，并且计算滤波器系数的运算量大；谱变形方法则是旁瓣抑制效果不是非常明显，针对强度较小的相邻目标往往没有效果。
发明内容
为了解决SAR图像目标旁瓣过高，且目标检测范围过大造成的影响，本发明采用改进的RDA方法提出了一种基于混沌调频信号作为雷达发射信号的图像旁瓣抑制方法。该方法通过对合成孔径雷达接收回波信号进行重构，利用合成孔径雷达发射的混沌调频信号的类随机特征，进行相关叠加来抑制SAR图像旁瓣。
本发明对合成孔径雷达图像旁瓣抑制过程是在现有RDA方法的距离压缩步骤与方位向傅里叶变换步骤之间增加了提取时延信息步骤、点扩散函数卷积步骤和方位时刻距离向信号组叠加步骤；在SAR处理器中的全程处理步骤有：
第一步：SAR信号处理器一方面接收地面接收站输出的回波信号E，另一方面对接收到的回波信号E进行距离向的匹配滤波，得到距离向压缩—脉压后的回波信号ER；
第二步：将第一步得到的距离向压缩—脉压后的回波信号ER进行距离向SequenceClean提取时延操作，得到每个方位时刻距离向包含每个目标的时延信息的时延矩阵B；
第三步：把得到的时延矩阵B卷积上点扩散函数P，在每个方位时刻可以得到一组包含相同时延信息脉压之后的时延—卷积信息S；
第四步：将第一步得到的距离向压缩—脉压后的回波信号ER与第三步的时延—卷积信息S进行叠加，得到在方位时刻下的方位向信息S方位时刻；
第五步：将所述S方位时刻进行距离向脉冲压缩之后得到的回波信号进行徙动校正后，再进行距离向压缩，得到SAR图像抑制旁瓣的回波信号S抑制。
在合成孔径雷达成像处理中，应用本发明的方法具有如下优点：与传统的chirp信号的作为发射信号的雷达相比，本发明方法可以在距离向对场景目标的旁瓣进行大幅度的抑制，并且由于带宽没有改变，使得图像保有了与chirp信号同样的主瓣宽度，没有降低分辨率(见图11和图5、图7与图9之间的对比)。
附图说明
图1是SAR***的结构框图。
图1A是RDA方法的处理流程图。
图2是本发明用于合成孔径雷达旁瓣抑制的流程图。
图3是原始点目标三维图。
图4是原始点目标距离向剖面图。
图5是频域加窗加权旁瓣抑制方法三维图。
图6是频域加窗加权旁瓣抑制方法距离向剖面图。
图7是切趾滤波器旁瓣抑制方法三维图。
图8是切趾滤波器旁瓣抑制方法距离向剖面图。
图9是谱变形旁瓣抑制方法的三维图。
图10是谱变形旁瓣抑制方法的距离向剖面图。
图11是本发明旁瓣抑制方法的三维图。
图12是本发明旁瓣抑制方法的距离向剖面图。
具体实施方式
下面将结合附图和仿真实例对本发明做进一步的详细说明。
本发明基于改进的RDA进行合成孔径雷达图像目标旁瓣抑制处理是在SAR信号处理器(参见图1所示)中完成的。具体地，是对雷达接收到的回波信号进行成像处理，得到旁瓣抑制后的SAR图像。具体的合成孔径雷达旁瓣抑制方法的处理流程，如图2所示。将过采样下的雷达原始数据记为回波信号E。本发明对合成孔径雷达图像旁瓣抑制过程是在现有RDA方法的距离压缩步骤与方位向傅里叶变换步骤之间增加了提取时延信息步骤、点扩散函数卷积步骤和方位时刻距离向信号组叠加步骤。
本发明对合成孔径雷达过采样数据的超分辨处理包括有下列九个处理步骤：
第一步：SAR信号处理器一方面接收地面接收站输出的回波信号E，另一方面对接收到的回波信号E进行距离向的匹配滤波，得到距离向压缩—脉冲后的回波信号ER；
在这个过程中，对回波信号进行距离向傅里叶变换是为了将回波信号从方位时域变换成方位频域的信号。也就是说，所述的距离向压缩—脉冲后的回波信号ER为方位频域信号。
第二步：采用Sequence Clean方法对第一步得到的ER进行图像目标时延信息提取，得到每个方位时刻距离向包含每个目标的时延信息的时延矩阵B；
在这个过程中，先进行距离向脉冲压缩，提高了距离向信噪比。并利用距离向聚焦后的信号的峰值，通过Sequence Clean方法，把目标场景中每个时延信息都提取出来。其中Sequence Clean提取时延的具体方法可以参考2002年1月Ranjan Bose等人在Transactions on Aerospace and Electronic Systems上发表的《Sequence Clean:AModified Deconvolution Technique for Microwave Images of Contiguous Targets》文章。
第三步：把得到的时延矩阵B卷积上点扩散函数P，在每个方位时刻可以得到一组包含相同时延信息脉压之后的时延—卷积信息S；
在第三步的处理过程中，由于时延信息中是以幅度以及位置的形式存在，所以还需要卷积上“由合成孔径雷达发射的”混沌调频信号的类随机信号，即将不同的混沌调频信号作为点扩散函数卷积上回波信号成为自相关函数，这样才能恢复回波信号到脉压之后的回波信号，从而得到时延—卷积信息S。
第四步：将第一步得到的距离向压缩—脉压后的回波信号ER与第三步的时延—卷积信息S进行叠加，得到在方位时刻下的方位向信息S方位时刻；
在第四步的处理过程中，利用混沌调频信号脉压之后旁瓣电平呈现随机走势的特点，通过ER与S之间的互相叠加，可以让呈现随机走势的旁瓣降低，让确定的主瓣得到提升，从而让旁瓣得到了有效的抑制。
第五步：将所述S方位时刻进行距离向脉冲压缩之后得到的回波信号进行徙动校正后，再进行距离向压缩，得到SAR图像抑制旁瓣的回波信号S抑制；
在第五步的处理过程中，可以参考2012年6月电子工业出版社、Ian G Cumming等著、洪文等译的《合成孔径雷达成像-算法与实现》一书中RDA成像算法中方位向脉冲压缩、徙动校正的过程。
在本发明的第一步中，SAR处理器接收到的回波信号E是二维数据，E采用矩阵形式表示为：
Na表示方位向点数；
Nf表示距离向点数；
n是一个变量，取值范围1至Na，表示为方位向第n个点；
m是一个变量，取值范围1至Nf，表示为距离向第m个点；
a1,1表示沿方位向上采集到的第1个方位时刻的距离向上的第1个采样点；
a1,2表示沿方位向上采集到的第1个方位时刻的距离向上的第2个采样点；
a1,m表示沿方位向上采集到的第1个方位时刻的距离向上的第m个采样点；
表示沿方位向上采集到的第1个方位时刻的距离向上的第Nf个采样点；
a2,1表示沿方位向上采集到的第2个方位时刻的距离向上的第1个采样点；
a2,2表示沿方位向上采集到的第2个方位时刻的距离向上的第2个采样点；
a2,m表示沿方位向上采集到的第2个方位时刻的距离向上的第m个采样点；
表示沿方位向上采集到的第2个方位时刻的距离向上的第Nf个采样点；
an,1表示沿方位向上采集到的第n个方位时刻的距离向上的第1个采样点；
an,2表示沿方位向上采集到的第n个方位时刻的距离向上的第2个采样点；
an,m表示沿方位向上采集到的第n个方位时刻的距离向上的第m个采样点；
表示沿方位向上采集到的第n个方位时刻的距离向上的第Nf个采样点；
表示沿方位向上采集到的第Na个方位时刻的距离向上的第1个采样点；
表示沿方位向上采集到的第Na个方位时刻的距离向上的第2个采样点；
表示沿方位向上采集到的第Na个方位时刻的距离向上的第m个采样点；
表示沿方位向上采集到的第Na个方位时刻的距离向上的第Nf个采样点。
在本发明的第一步中，回波信号E经距离向脉冲压缩处理后，得到的脉冲压缩后回波信号ER，ER采用矩阵形式表示为：
Na表示方位向点数；
Nf表示距离向点数；
n是一个变量，取值范围1至Na，表示为方位向第n个点；
m是一个变量，取值范围1至Nf，表示为距离向第m个点；
表示方位向上的第1个方位时刻的距离向上的第1个点；
表示方位向上的第1个方位时刻的距离向上的第2个点；
表示方位向上的第1个方位时刻的距离向上的第m个点；
表示方位向上的第1个方位时刻的距离向上的第Nf个点；
表示方位向上的第2个方位时刻的距离向上的第1个点；
表示方位向上的第2个方位时刻的距离向上的第2个点；
表示方位向上的第2个方位时刻的距离向上的第m个点；
表示方位向上的第2个方位时刻的距离向上的第Nf个点；
表示方位向上的第n个方位时刻的距离向上的第1个点；
表示方位向上的第n个方位时刻的距离向上的第2个点；
表示方位向上的第n个方位时刻的距离向上的第m个点；
表示方位向上的第n个方位时刻的距离向上的第Nf个点；
表示方位向上的第Na个方位时刻的距离向上的第1个点；
表示方位向上的第Na个方位时刻的距离向上的第2个点；
表示方位向上的第Na个方位时刻的距离向上的第m个点；
表示方位向上的第Na个方位时刻的距离向上的第Nf个点。
在本发明的第二步中，采用Sequence Clean方法从距离向脉冲压缩后的回波信号ER中提取出每个方位时刻距离向的图像目标时延信息矩阵，即时延矩阵B的矩阵形式表示为：
Na表示方位向点数；
Nf表示距离向点数；
n是一个变量，取值范围1至Na，表示为方位向第n个点；
m是一个变量，取值范围1至Nf，表示为距离向第m个点；
b1,1表示方位向上的第1个方位时刻的距离向上的第1个点；
b1,2表示方位向上的第1个方位时刻的距离向上的第2个点；
b1,m表示方位向上的第1个方位时刻的距离向上的第m个点；
表示方位向上的第1个方位时刻的距离向上的第Nf个点；
b2,1表示方位向上的第2个方位时刻的距离向上的第1个点；
b2,2表示方位向上的第2个方位时刻的距离向上的第2个点；
b2,m表示方位向上的第2个方位时刻的距离向上的第m个点；
表示方位向上的第2个方位时刻的距离向上的第Nf个点；
bn,1表示方位向上的第n个方位时刻的距离向上的第1个点；
bn,2表示方位向上的第n个方位时刻的距离向上的第2个点；
bn,m表示方位向上的第n个方位时刻的距离向上的第m个点；
表示方位向上的第n个方位时刻的距离向上的第Nf个点；
表示方位向上的第Na个方位时刻的距离向上的第1个点；
表示方位向上的第Na个方位时刻的距离向上的第2个点；
表示方位向上的第Na个方位时刻的距离向上的第m个点；
表示方位向上的第Na个方位时刻的距离向上的第Nf个点。
在本发明的第三步中，点扩散函数P采用矩阵形式表示为：
Nb表示任意方位时刻点扩散函数组中函数的个数；
Nd表示函数距离向采样点数；
n′是一个变量，取值范围1至Nb，表示为方位任意时刻第n′个信号；
m′是一个变量，取值范围1至Nd，表示为信号的第m′个点；
p1,1表示方位向上任意一个时刻的第1个函数的第1个点；
p1,2表示方位向上任意一个时刻的第1个函数的第2个点；
p1,m′表示方位向上任意一个时刻的第1个函数的第m′个点；
表示方位向上任意一个时刻的第1个函数的第Nd个点；
p2,1表示方位向上任意一个时刻的第2个函数的第1个点；
p2,2表示方位向上任意一个时刻的第2个函数的第2个点；
p2,m′表示方位向上任意一个时刻的第2个函数的第m′个点；
表示方位向上任意一个时刻的第2个函数的第Nd个点；
pn′,1表示方位向上任意一个时刻的第n′个函数的第1个点；
pn′,2表示方位向上任意一个时刻的第n′个函数的第2个点；
pn′,m′表示方位向上任意一个时刻的第n′个函数的第m′个点；
表示方位向上任意一个时刻的第n′个函数的第Nd个点；
表示方位向上任意一个时刻的第Nb个函数的第1个点；
表示方位向上任意一个时刻的第Nb个函数的第2个点；
表示方位向上任意一个时刻的第Nb个函数的第m′个点；
表示方位向上任意一个时刻的第Nb个函数的第Nd个点；
在本发明的第四步中，所述的方位时刻下的方位向信息S方位时刻采用矩阵形式表示为：
Na表示方位向点数；
Nf表示距离向点数；
n′是一个变量，取值范围1至Na，表示为方位任意时刻第n′个信号；
m′是一个变量，取值范围1至Nf，表示为信号的第m′个点；
a'1,1表示卷积后的任意一个时刻的第1个信号的第1个采样点；
a'1,2表示卷积后的任意一个时刻的第1个信号的第2个采样点；
a'1,m′表示卷积后的任意一个时刻的第1个信号的第3个采样点；
表示卷积后的任意一个时刻的第1个信号的第Nf个采样点；
a'2,1表示卷积后的任意一个时刻的第2个信号的第1个采样点；
a'2,2表示卷积后的任意一个时刻的第2个信号的第2个采样点；
a'2,m′表示卷积后的任意一个时刻的第2个信号的第3个采样点；
表示卷积后的任意一个时刻的第2个信号的第Nf个采样点；
a'n′,1表示卷积后的任意一个时刻的第3个信号的第1个采样点；
a'n′,2表示卷积后的任意一个时刻的第3个信号的第2个采样点；
a'n′,m′表示卷积后的任意一个时刻的第3个信号的第3个采样点；
表示卷积后的任意一个时刻的第3个信号的第Nf个采样点；
表示卷积后的任意一个时刻的第Na个信号的第1个采样点；
表示卷积后的任意一个时刻的第Na个信号的第2个采样点；
表示卷积后的任意一个时刻的第Na个信号的第3个采样点；
表示卷积后的任意一个时刻的第Na个信号的第Nf个采样点；
经本发明的第五步处理后，得到的SAR图像抑制旁瓣的回波信号S抑制采用矩阵形式表示为：
Na表示方位向采样点数；
Nf表示距离向采样点数；
n是一个变量，取值范围1至Na，表示为方位向第n个点；
m是一个变量，取值范围1至Nf，表示为距离向第m个点；
a″1,1表示叠加后方位向上的第1个方位时刻的距离上的第1个点；
a″1,2表示叠加后方位向上的第1个方位时刻的距离上的第2个点；
a″1,m表示叠加后方位向上的第1个方位时刻的距离上的第m个点；
表示叠加后方位向上的第1个方位时刻的距离上的第Nf个点；
a″2,1表示叠加后方位向上的第2个方位时刻的距离上的第1个点；
a″2,2表示叠加后方位向上的第2个方位时刻的距离上的第2个点；
a″2,m表示叠加后方位向上的第2个方位时刻的距离上的第m个点；
表示叠加后方位向上的第2个方位时刻的距离上的第Nf个点；
a″n,1表示叠加后方位向上的第n个方位时刻的距离上的第1个点；
a″n,2表示叠加后方位向上的第n个方位时刻的距离上的第2个点；
a″n,m表示叠加后方位向上的第n个方位时刻的距离上的第m个点；
表示叠加后方位向上的第n个方位时刻的距离上的第Nf个点；
表示叠加后方位向上的第Na个方位时刻的距离上的第1个点；
表示叠加后方位向上的第Na个方位时刻的距离上的第2个点；
表示叠加后方位向上的第Na个方位时刻的距离上的第m个点；
表示叠加后方位向上的第Na个方位时刻的距离上的第Nf个点。
仿真实例
采用本设计发明的用于合成孔径雷达图像旁瓣抑制处理的方法进行仿真，其参数设置如表1。在场景中心的位置沿着距离向布置了三个点目标，他们的幅值归一化之后分别为1、0.05、1。他们前两个点的位置离得很近，第三个点离前两个点较远。
表1成像参数设置
图3和图4是采用传统的线性调频信号作为雷达发射信号，再使用距离多普勒算法对点目标进行仿真成像处理结果，包括三维图形和距离向刨面图。
图5和图6是对图3所仿真的点目标图像进行距离向频域加窗旁瓣抑制处理结果，图中给出的是三维图像和距离向刨面图。
图7和图8是对图3所仿真的点目标图像进行距离向切趾滤波器旁瓣抑制处理结果，图中给出的是三维图像和距离向刨面图。
图9和图10是对图3所仿真的点目标图像进行谱变形旁瓣抑制处理结果，图中给出的是三维图像和距离向刨面图。
图11和图12是采用本发明的使用混沌调频信号作为雷达发射信号的旁瓣抑制处理方法，图中给出的是三维图像和距离向刨面图。
通过对成像结果的对比，即图3与图5、图7、图9和图11的对比，图4与图6、图8、图10和图12的对比，可以发现：
传统的线性调频信号的旁瓣抑制算法各有优缺点，频域加窗加权法虽然能够大幅度降低目标的旁瓣，但是由于主瓣展宽，损失了分辨力，离主瓣较近的弱小目标就被淹没在旁瓣电平中了；切趾滤波器旁瓣抑制方法虽然保持了分辨力，但是旁瓣抑制效果不明显，并且由于在计算过程中需要先验条件，所以给旁瓣抑制处理带来了难度；谱变形方法可以降低旁瓣，但是也是效果不明显。
采用本发明方法不仅可以大幅度降低旁瓣，并且保证了分辨力，使得强目标附近的弱小目标不会消失在旁瓣中，提高了雷达在目标识别中的准确程度。

Claims (2)

1.一种基于改进的RDA进行合成孔径雷达图像目标旁瓣抑制的方法，其特征在于：该方法是在现有RDA方法的距离压缩步骤与方位向傅里叶变换步骤之间增加了提取时延信息步骤、点扩散函数卷积步骤和方位时刻距离向信号组叠加步骤；在SAR处理器中的全程处理步骤有：
第一步：SAR信号处理器一方面接收地面接收站输出的回波信号E，另一方面对接收到的回波信号E进行距离向的匹配滤波，得到距离向压缩—脉压后的回波信号ER；
第二步：将第一步得到的距离向压缩—脉压后的回波信号ER进行距离向SequenceClean提取时延操作，得到每个方位时刻距离向包含每个目标的时延信息的时延矩阵B；
时延矩阵B的矩阵形式表示为：
Na表示方位向点数；
Nf表示距离向点数；
n是一个变量，取值范围1至Na，表示为方位向第n个点；
m是一个变量，取值范围1至Nf，表示为距离向第m个点；
第三步：把得到的时延矩阵B卷积上点扩散函数P，在每个方位时刻得到一组包含相同时延信息脉压之后的时延—卷积信息S；
在第三步的处理过程中，由于时延信息中是以幅度以及位置的形式存在，所以还需要卷积上由合成孔径雷达发射的混沌调频信号的类随机信号，即将不同的混沌调频信号作为点扩散函数卷积上回波信号成为自相关函数，这样才能恢复回波信号到脉压之后的回波信号，从而得到时延—卷积信息S；
点扩散函数P采用矩阵形式表示为：
Nb表示任意方位时刻点扩散函数组中函数的个数；
Nd表示函数距离向采样点数；
n′是一个变量，取值范围1至Nb，表示为方位任意时刻第n′个信号；
m′是一个变量，取值范围1至Nd，表示为信号的第m′个点；
第四步：将第一步得到的距离向压缩—脉压后的回波信号ER与第三步的时延—卷积信息S进行叠加，得到在方位时刻下的方位向信息S方位时刻；
方位时刻下的方位向信息S方位时刻采用矩阵形式表示为：
第五步：将所述S方位时刻进行距离向脉冲压缩之后得到的回波信号进行徙动校正后，再进行距离向压缩，得到SAR图像抑制旁瓣的回波信号S抑制；
SAR图像抑制旁瓣的回波信号S抑制采用矩阵形式表示为：
Na表示方位向采样点数；
Nf表示距离向采样点数；
n是一个变量，取值范围1至Na，表示为方位向第n个点；
m是一个变量，取值范围1至Nf，表示为距离向第m个点。
2.根据权利要求1所述的基于改进的RDA进行合成孔径雷达图像目标旁瓣抑制的方法，其特征在于：在现有SAR处理器中，在设置了成像参数的情况下，不仅大幅度降低旁瓣，并且保证了分辨力，使得强目标附近的弱小目标不会消失在旁瓣中，提高了雷达在目标识别中的准确程度；
成像参数有：雷达等效速度V为7100m/s；雷达波长λ为3×108/5.3×109m；发射脉冲时宽Tr为1μs；天线真实孔径La为6m；参考斜距Rref为693km；发射信号的调频率Kr为5×1012Hz/s；参考目标的等效斜视角为0°；波束斜视角为30°；距离采样率Fr为70MHz；方位采样率Fa为3500Hz。
CN201710013702.7A 2017-01-09 2017-01-09 基于改进的rda进行合成孔径雷达图像目标旁瓣抑制的方法 Active CN106772376B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710013702.7A CN106772376B (zh) 2017-01-09 2017-01-09 基于改进的rda进行合成孔径雷达图像目标旁瓣抑制的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710013702.7A CN106772376B (zh) 2017-01-09 2017-01-09 基于改进的rda进行合成孔径雷达图像目标旁瓣抑制的方法

Publications (2)

Publication Number Publication Date
CN106772376A CN106772376A (zh) 2017-05-31
CN106772376B true CN106772376B (zh) 2019-01-25

Family

ID=58950621

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710013702.7A Active CN106772376B (zh) 2017-01-09 2017-01-09 基于改进的rda进行合成孔径雷达图像目标旁瓣抑制的方法

Country Status (1)

Country Link
CN (1) CN106772376B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108226879B (zh) * 2017-12-21 2021-08-13 北京遥感设备研究所 一种基于多通道的sar地形散射干扰抑制方法
CN108132467B (zh) * 2017-12-23 2018-11-30 成都汇蓉国科微***技术有限公司 基于增强型adc的dsp+fpga双基前视sar成像方法及成像装置
CN110609284B (zh) * 2019-09-23 2023-08-18 成都盟升科技有限公司 一种sar中断数据恢复成像的方法
CN110806577B (zh) * 2019-11-06 2022-03-18 中国科学院电子学研究所 合成孔径雷达的聚焦成像方法及装置、设备、存储介质
CN111881414B (zh) * 2020-07-29 2024-03-15 中南大学 一种基于分解理论的合成孔径雷达图像质量评估方法
CN112305538B (zh) * 2020-09-22 2023-06-20 中国电子科技集团公司第二十九研究所 一种合成孔径雷达图像旁瓣抑制处理方法
CN112881982B (zh) * 2021-01-18 2022-11-15 中国人民解放军海军航空大学 一种捷变频雷达抑制密集转发干扰方法
CN113189587B (zh) * 2021-04-07 2023-10-20 北京理工大学 雷达的单目标识别方法、***及可读存储介质
CN112946650B (zh) * 2021-04-15 2023-08-15 中山大学 一站固定双站低频超宽带sar运动目标检测与成像方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101089655A (zh) * 2006-06-14 2007-12-19 中国科学院电子学研究所 一种使用混沌信号的合成孔径雷达***
CN103235288A (zh) * 2013-04-17 2013-08-07 中国科学院空间科学与应用研究中心 基于频域的超低旁瓣混沌雷达信号生成及数字实现方法
CN103558601A (zh) * 2013-11-01 2014-02-05 北京航空航天大学 一种用于合成孔径雷达过采样数据处理的超分辨方法
CN103869313A (zh) * 2014-02-17 2014-06-18 中国科学院电子学研究所 一种多通道合成孔径雷达距离向模糊的抑制方法及装置
CN104635230A (zh) * 2015-03-05 2015-05-20 北京航空航天大学 一种用于mimo-sar近场测量成像方位向旁瓣抑制的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101089655A (zh) * 2006-06-14 2007-12-19 中国科学院电子学研究所 一种使用混沌信号的合成孔径雷达***
CN103235288A (zh) * 2013-04-17 2013-08-07 中国科学院空间科学与应用研究中心 基于频域的超低旁瓣混沌雷达信号生成及数字实现方法
CN103558601A (zh) * 2013-11-01 2014-02-05 北京航空航天大学 一种用于合成孔径雷达过采样数据处理的超分辨方法
CN103869313A (zh) * 2014-02-17 2014-06-18 中国科学院电子学研究所 一种多通道合成孔径雷达距离向模糊的抑制方法及装置
CN104635230A (zh) * 2015-03-05 2015-05-20 北京航空航天大学 一种用于mimo-sar近场测量成像方位向旁瓣抑制的方法

Also Published As

Publication number Publication date
CN106772376A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106772376B (zh) 基于改进的rda进行合成孔径雷达图像目标旁瓣抑制的方法
Stankwitz et al. Nonlinear apodization for sidelobe control in SAR imagery
Sjögren et al. Suppression of clutter in multichannel SAR GMTI
CN103983972B (zh) 一种快速压缩传感三维sar稀疏成像方法
CN109116352B (zh) 一种圆扫isar模式船只超分辨率成像方法
CN109669182B (zh) 无源双基地sar动/静目标联合稀疏成像方法
CN107918124A (zh) 带有方位空变校正的机载大斜视高分辨sar成像方法
CN102565772B (zh) 基于sar子孔径序列图像的海洋动态信息提取方法
CN113589287A (zh) 合成孔径雷达稀疏成像方法、装置、电子设备及存储介质
CN107064904B (zh) 一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法
Rahman Focusing moving targets using range migration algorithm in ultra wideband low frequency synthetic aperture radar
CN109143236B (zh) 适用于复杂飞行轨迹的双基聚束sar大场景成像方法
CN109143235A (zh) 一种双基前视合成孔径雷达地面运动目标检测方法
Kelly et al. RFI suppression and sparse image formation for UWB SAR
Borcea et al. Synthetic aperture radar imaging with motion estimation and autofocus
Linnehan et al. An efficient means to mitigate wavefront curvature effects in polar format processed SAR imagery
CN113671497B (zh) 基于圆柱对称模型的单通道sar目标三维坐标提取方法
Mao et al. Realization of airborne forward-looking radar super-resolution algorithm based on GPU frame
Wang et al. Sidelobe reduction based on spectrum reshaping in microwave imaging
CN114002666A (zh) 任意天线构型下星载ati-sar洋流流速提取方法及设备
CN112686871A (zh) 基于改进对数比算子和Gabor_ELM的SAR图像变化检测方法
Jakowatz Jr et al. Comparison of algorithms for use in real-time spotlight-mode SAR image formation
CN113191538B (zh) 一种基于深度学习的多频sar图像解模糊方法
Burns Forward looking GPR sidelobe reduction using L1-norm minimization
Rau et al. A directional image decomposition for ultra-wideband 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