CN110333489B - 采用cnn与rsva组合对sar回波数据旁瓣抑制的处理方法 - Google Patents

采用cnn与rsva组合对sar回波数据旁瓣抑制的处理方法 Download PDF

Info

Publication number
CN110333489B
CN110333489B CN201910673129.1A CN201910673129A CN110333489B CN 110333489 B CN110333489 B CN 110333489B CN 201910673129 A CN201910673129 A CN 201910673129A CN 110333489 B CN110333489 B CN 110333489B
Authority
CN
China
Prior art keywords
azimuth
distance
moment
imaging result
time
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
CN201910673129.1A
Other languages
English (en)
Other versions
CN110333489A (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 CN201910673129.1A priority Critical patent/CN110333489B/zh
Publication of CN110333489A publication Critical patent/CN110333489A/zh
Application granted granted Critical
Publication of CN110333489B publication Critical patent/CN110333489B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/28—Details of pulse systems
    • G01S7/2813—Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays
    • 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
    • G01S7/417—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 involving the use of neural networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法，第一，利用SAR信号的线性可加性，通过设置不同点目标间隔和强度，生成仿真回波数据，利用其主瓣位置与点目标位置之间的对应关系，从仿真回波数据中还原去除旁瓣的回波数据，建立训练集；第二，将仿真回波数据通过插值和滑窗规范化获取切片样本；第三，建立卷积神经网络模型，该网络结构包含四层卷积层，一个全连接层和一个神经元作为输出；第四，将滑窗后的训练集输入卷积神经网络模型，采用随机梯度下降法对网络模型进行训练，得到一维旁瓣抑制模型；第五，对于实时回波数据进行插值，并结合训练好的一维旁瓣抑制模型，得到压缩后的预测值，利用RSVA重构生成旁瓣抑制后的回波数据。

Description

采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法
技术领域
本发明涉及对SAR回波数据旁瓣进行抑制的方法，更特别地说，是指一种采用CNN与SVA组合对SAR回波数据旁瓣抑制的处理方法。
背景技术
星载SAR(Synthetic Aperture Radar，合成孔径雷达)是一种具有全天候、全天时的能力的成像系统，不受云、雨、雾、黑夜等自然条件的限制，在军事领域、地形测绘等许多领域有广泛的应用价值。星载SAR系统原始回波信号的模拟对SAR系统设计、成像算法研究等有重要意义。
1999年10月哈尔滨工业大学出版社出版、刘永坦编著的《雷达成像技术》指出，合成孔径雷达(SAR，Synthetic Aperture Radar)安装在运动平台上，按照一定的重复频率发射、接收脉冲，形成回波信号。SAR系统的结构框图如图1所示，SAR系统包括有星上雷达系统、卫星平台及数据下传系统和地面系统三部分，对合成孔径雷达成像处理是在地面系统中完成的。地面系统通过地面接收站接收卫星平台及数据下传系统下发的回波信号，该回波信号经SAR信号处理器进行成像处理，获得SAR回波数据；所述SAR回波数据存储于备档操作系统中。
卷积神经网络(convolutional neural network，CNN)是一类包含卷积计算且具有深度结构的前馈神经网络(Feedforward Neural Networks)，是深度学习(deeplearning)的代表算法之一。2016年6月第1版，吴岸城编著的《神经网络与深度学习》指出，卷积神经网络的结构是获得一个输入(输入向量)，并将其通过一系列的隐藏层转换。每个隐层由一组神经元组成，其中每个神经元和前一层神经元完全连接，在单层的神经元的功能完全独立并且不共享任何连接。
2010年第40卷第12期《中国科学》，公开了“一种基于改进SVA的SAR旁瓣抑制算法”，作者倪崇等。此文中介绍了一种鲁棒的SVA(Robust SVA，RSVA)算法，该RSVA算法能够有效地抑制旁瓣，但有时会降低主瓣能量。
空间分辨率、旁瓣比、模糊度是衡量SAR回波数据质量的重要指标。传统体制的SAR信号处理器采用线性调频信号作为雷达发射信号波形，在不加权条件下，其脉冲压缩后的峰值旁瓣比约为－13.2dB。SAR回波数据中强的目标的旁瓣会掩盖在强目标附近的弱目标的主瓣，使得弱目标无法检测出。SAR回波数据还会因为旁瓣干扰造成一定的模糊，使得场景和目标不容易区分，旁瓣引起的模糊使得SAR回波数据相邻地物之间的区别在视觉上没有真实情况明显，也就降低了相邻目标类别之间的可区分性。
发明内容
本发明要解决的技术问题在于，针对现有旁瓣抑制技术中的缺陷，提出一种新的旁瓣抑制优化目标函数，并基于CNN与RSVA组合建立一种SAR回波数据旁瓣抑制方法。本发明第一方面采用离线SAR回波信号的可加性，通过设置SAR系统及信号处理器中的各项参数包括脉冲宽度，线性调频率、脉冲重复频率、距离向采样率、轨道参数等，以及地面接收站的各个反射系数RCS，生成二维回波，通过CSA算法进行成像处理后，拆分距离向与方位向生成适用于旁瓣抑制的训练集，对信号进行预插值和滑窗规范化的处理，通过分析旁瓣抑制前后的信号特征，构建SCNN网络模型，利用训练集对卷积神经网络模型进行训练，得出旁瓣抑制的模型。本发明第二方面应用实时SAR回波二维信号，首先将SAR回波二维信号通过CSA算法进行成像处理后，对生成的二维回波数据分解为一维方位向数据，对一维数据进行与训练集生成过程中相同的预插值和滑窗规范化的处理，将处理后的切片输入已训练好的模型，得到各个切片的预测结果后，通过寻址重构得到方位向旁瓣抑制后的数据，再对该数据进行相同距离向上旁瓣抑制处理，最后利用SVA算法对网络输出结果进行优化。
本发明的采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法，包括两部分内容，一方面是训练模型，另一方面是应用训练模型进行的实时SAR回波数据的旁瓣抑制处理。第一，利用SAR信号的线性可加性，通过设置不同点目标间隔和强度，生成仿真回波数据，利用其主瓣位置与点目标位置之间的对应关系，从仿真回波数据中还原去除旁瓣的回波数据，建立训练集；第二，将仿真回波数据通过插值和滑窗规范化获取切片样本；第三，建立卷积神经网络模型，该网络结构包含四层卷积层，一个全连接层和一个神经元作为输出；第四，将滑窗后的训练集输入卷积神经网络模型，采用随机梯度下降法对网络模型进行训练，得到一维旁瓣抑制模型；第五，对于实时回波数据进行插值，并结合训练好的一维旁瓣抑制模型，得到压缩后的预测值，利用RSVA重构生成旁瓣抑制后的回波数据。
本发明采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法的优点在于：
①本发明提出的卷积神经网络的SAR旁瓣抑制算法的可靠性强，它利用更为普适的旁瓣抑制目标函数，忽略了其不同目标间相位存在的干扰，具备更好的抑制效果。
②本发明提出的基于卷积神经网络的SAR旁瓣抑制算法的泛化能力强，计算复杂度低。该卷积神经网络结构较简单，仅有四层卷积层，不易导致模型过拟合。
③本发明提出的基于卷积神经网络的SAR旁瓣抑制算法的实用性强。本发明将SAR图像直接输入训练好的网络模型，与传统SVA算法结合即可生成预测图像。
④本发明提出的基于卷积神经网络的SAR旁瓣抑制算法的适用范围广。由于忽略不同系统参数导致的主旁瓣不同，采用主瓣宽度作为唯一标准，通过插值实现训练集测试集的匹配，因此可以适用于星载、机载等各个系统参数下。
⑤本发明提出基于卷积神经网络的SAR旁瓣抑制算法具有图像质量高的特点。由于采用新算法后，雷达图像旁瓣很低，主瓣能量丢失也减少，因此，图像质量较高，更易判读。
附图说明
图1是传统SAR系统的结构框图。
图2是本发明的一种基于完全互补序列的超低旁瓣合成孔径雷达成像方法流程图。
图3是本发明中使用的SCNN结构图。
图4是本发明滑窗处理的示意图。
图5是实施例中的单点目标成像方位向剖面图。
图6是实施例中的单点目标在图像中成像三维图。
图7是实施例中的单点目标成像利用RVA进行旁瓣抑制方位向剖面图。
图8是实施例中的单点目标成像利用RFA进行旁瓣抑制在图像中成像三维图。
图9是实施例中的单点目标成像利用本发明进行旁瓣抑制方位向剖面图。
图10是实施例中的单点目标成像利用本发明进行旁瓣抑制在图像中成像三维图。
图11是实施例中的多(两)点目标成像方位向剖面图。
图12是实施例中的多(两)点目标在图像中成像三维图。
图13是实施例中的多(两)点目标成像利用SVA进行旁瓣抑制方位向剖面图。
图14是实施例中的多(两)点点目标成像利用SVA进行旁瓣抑制在图像中成像三维图。
图15是实施例中的多(两)点目标成像利用本发明进行旁瓣抑制方位向剖面图。
图16是实施例中的多(两)点目标成像利用本发明进行旁瓣抑制在图像中成像三维图。
具体实施方式
下面将结合附图和实施例对本发明做进一步的详细说明。
本发明的采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法，包括两部分内容，一方面是训练模型，另一方面是应用训练模型进行的实时SAR回波数据的旁瓣抑制处理。
构建一维旁瓣抑制训练模型
步骤A，获取SAR仿真回波数据；
在本发明中，为了获得训练集的数据信息，运行于计算机中的计算机程序－基于matlab(版本号R2014a)平台进行仿真。在其软件中设置SAR雷达轨道参数和SAR系统参数。
软件中设置的相关SAR参数有：机载飞行速度、雷达波长、参考距离、距离向采样率、距离向调频率、脉冲重复频率、距离向过采样率、方位向过采样率、散射单元系数、散射单元间隔。
在本发明中，将matlab仿真输出的二维仿真SAR回波信号，记为
Figure GDA0002793063940000041
并对所述
Figure GDA0002793063940000042
进行成像处理，得到SAR仿真回波数据
Figure GDA0002793063940000043
在本发明中，仿真SAR处理器接收到的回波信号
Figure GDA0002793063940000044
是二维数据，采用矩阵形式表示为：
Figure GDA0002793063940000045
Na表示方位向采样点总数；
Nf表示距离向采样点总数；
x是一个变量，取值范围1至Na，表示为方位向第x个点，也称为方位向的标识号；
y是一个变量，取值范围1至Nf，表示为距离向第y个点，也称为距离向的标识号；
Figure GDA0002793063940000046
表示仿真信号沿方位向上采集到的第1个方位时刻的距离向上的第1个采样点；
Figure GDA0002793063940000047
表示仿真信号沿方位向上采集到的第1个方位时刻的距离向上的第2个采样点；
Figure GDA0002793063940000048
表示仿真信号沿方位向上采集到的第1个方位时刻的距离向上的第y个采样点；
Figure GDA0002793063940000049
表示仿真信号沿方位向上采集到的第1个方位时刻的距离向上的第Nf个采样点；
Figure GDA00027930639400000410
表示仿真信号沿方位向上采集到的第2个方位时刻的距离向上的第1个采样点；
Figure GDA00027930639400000411
表示仿真信号沿方位向上采集到的第2个方位时刻的距离向上的第2个采样点；
Figure GDA00027930639400000412
表示仿真信号沿方位向上采集到的第2个方位时刻的距离向上的第y个采样点；
Figure GDA0002793063940000051
表示仿真信号沿方位向上采集到的第2个方位时刻的距离向上的第Nf个采样点；
Figure GDA0002793063940000052
表示仿真信号沿方位向上采集到的第x个方位时刻的距离向上的第1个采样点；
Figure GDA0002793063940000053
表示仿真信号沿方位向上采集到的第x个方位时刻的距离向上的第2个采样点；
Figure GDA0002793063940000054
表示仿真信号沿方位向上采集到的第x个方位时刻的距离向上的第y个采样点；
Figure GDA0002793063940000055
表示仿真信号沿方位向上采集到的第x个方位时刻的距离向上的第Nf个采样点；
Figure GDA0002793063940000056
表示仿真信号沿方位向上采集到的第Na个方位时刻的距离向上的第1个采样点；
Figure GDA0002793063940000057
表示仿真信号沿方位向上采集到的第Na个方位时刻的距离向上的第2个采样点；
Figure GDA0002793063940000058
表示仿真信号沿方位向上采集到的第Na个方位时刻的距离向上的第y个采样点；
Figure GDA0002793063940000059
表示仿真信号沿方位向上采集到的第Na个方位时刻的距离向上的第Nf个采样点。
在本发明中，采用CSA方法对仿真的回波信号矩阵
Figure GDA00027930639400000510
进行成像处理后，得到仿真的成像结果矩阵
Figure GDA00027930639400000511
的矩阵表示为：
Figure GDA00027930639400000512
Nb是方位向上能够有效产生所述Na点接收回波所对应的场景目标的采样长度(即采样点数)，亦为观测矩阵
Figure GDA00027930639400000513
的列数；
Nf表示距离向采样点总数；
x是一个变量，取值范围1至Nb，表示为成像结果方位向的第x个像素点，也称为成像结果方位向的标识号；
y是一个变量，取值范围1至Nf，表示为成像结果距离向的第y个像素点，也称为成像结果距离向的标识号；
Figure GDA00027930639400000514
表示仿真成像结果在方位向上对应第1方位时刻的距离向第1像素点；
Figure GDA00027930639400000515
表示仿真成像结果在方位向上对应第1方位时刻的距离向第2像素点；
Figure GDA00027930639400000516
表示仿真成像结果在方位向上对应第1方位时刻的距离向第y像素点；
Figure GDA00027930639400000517
表示仿真成像结果在方位向上对应第1方位时刻的距离向第Nf像素点；
Figure GDA00027930639400000518
表示仿真成像结果在方位向上对应第2方位时刻的距离向第1像素点；
Figure GDA00027930639400000519
表示仿真成像结果在方位向上对应第2方位时刻的距离向第2像素点；
Figure GDA0002793063940000061
表示仿真成像结果在方位向上对应第2方位时刻的距离向第y像素点；
Figure GDA0002793063940000062
表示仿真成像结果在方位向上对应第2方位时刻的距离向第Nf像素点；
Figure GDA0002793063940000063
表示仿真成像结果在方位向上对应第x方位时刻的距离向第1像素点；
Figure GDA0002793063940000064
表示仿真成像结果在方位向上对应第x方位时刻的距离向第2像素点；
Figure GDA0002793063940000065
表示仿真成像结果在方位向上对应第x方位时刻的距离向第y像素点；
Figure GDA0002793063940000066
表示仿真成像结果在方位向上对应第x方位时刻的距离向第Nf像素点；
Figure GDA0002793063940000067
表示仿真成像结果在方位向上对应第Nb方位时刻的距离向第1像素点；
Figure GDA0002793063940000068
表示仿真成像结果在方位向上对应第Nb方位时刻的距离向第2像素点；
Figure GDA0002793063940000069
表示仿真成像结果在方位向上对应第Nb方位时刻的距离向第y像素点；
Figure GDA00027930639400000610
表示仿真成像结果在方位向上对应第Nb方位时刻的距离向第Nf像素点。
在本发明中，CSA(Chinp Scaling)方法参考了2012年6月第1次印刷《合成孔径雷达成像－算法与实现》，洪文编译，第7章，第192-210页的内容。
步骤B，方位向上无旁瓣的信号集获取；
步骤B1，取出成像结果矩阵
Figure GDA00027930639400000611
中的距离向上对应第y距离时刻的成像后方位信号，记为
Figure GDA00027930639400000612
且
Figure GDA00027930639400000613
其成像信号满足形式：
Figure GDA00027930639400000614
其中，g(t)是所接收的回波，N是场景中的所有目标的个数，i是点目标的计数，σi是第i点目标的后向散射系数的赋值，Pi是第i个点目标的点扩散函数，t是信号所对应的时间，τ是信号的时延，e是对数符号，
Figure GDA00027930639400000615
中上角标的j是虚部，θi是第i个散射单元造成的随机相位。
成像后的方位信号
Figure GDA00027930639400000616
为所设置的SAR系统及信号处理器中的各项参数中不同组合的地面接收站的各个反射系数的信号，通过记录在步骤A中不同点目标的间距和强度，记录其主瓣位置与主瓣宽度，在信号中根据点目标位置与主瓣的对应关系，保留主瓣部分的信息，去除信号中的旁瓣信息，得到
Figure GDA00027930639400000617
所对应的标签
Figure GDA00027930639400000618
距离向上对应第x距离时刻的方位信号的无旁瓣信号标签记为
Figure GDA00027930639400000619
步骤B2，将成像后的方位信号
Figure GDA00027930639400000620
与对应的标签
Figure GDA00027930639400000621
进行线性插值得到信号
Figure GDA00027930639400000622
与
Figure GDA00027930639400000623
其线性插值公式为
Figure GDA00027930639400000624
利用上式计算插值中
Figure GDA0002793063940000071
的值，其中m为所需插值的倍数，n为插值后相邻方位向间的位置，得到插值之后的信号，以成像后的方位信号
Figure GDA0002793063940000072
为例，其插值得到的信号
Figure GDA0002793063940000073
SAR回波数据中主瓣宽度太小，卷积神经网络的感受野过大，无法直接提取其中旁瓣与主瓣相关信息，故对SAR数据进行了预插值处理，将主瓣拓展成为大约主瓣内包含N主瓣个像素点的数据。
步骤B3，对插值后的信号
Figure GDA0002793063940000074
与
Figure GDA0002793063940000075
进行滑窗规范化，建立SCNN训练所需方位向输入输出对，如图4所示，为两组信号的滑窗示意，将每一个窗长长度的样本送入网络中进行学习，其中步长为l，窗长为WK窗长。每一个样本可以表达为
Figure GDA0002793063940000076
h为信号采样的各个值，式中WK窗终＝WK窗始+WK窗长，在经历一个步长后，下一个样本可以表达为
Figure GDA0002793063940000077
其所学的值为其样本中心位置处的值
Figure GDA0002793063940000078
卷积神经网络常常通过规范化提高收敛训练速度，SAR回波数据中最高点的强度可达数万，最低仅为0，传统的规范化无法适用差别如此大的数据集，故对每一滑窗样本进行规范化如公式(3)，以减少数据层次过大导致的网络训练模型不佳。
Figure GDA0002793063940000079
式中Amax为考虑到数据存储可能造成的丢失所设置的规范化值，使输出范围为[0,Amax]，hp为样本I里每一个像素点的数值，hmax为每一样本I里的最大值，hp规范为规范化后的值。
步骤C，距离向上无旁瓣的信号标签获取；
步骤C1，取出成像结果矩阵
Figure GDA00027930639400000710
中的方位向上对应第x方位时刻的距离信号，记为
Figure GDA00027930639400000711
且
Figure GDA00027930639400000712
其成像信号满足公式(1)。成像后的距离信号
Figure GDA00027930639400000713
为所设置的SAR系统及信号处理器中的各项参数中不同组合的地面接收站的各个反射系数的信号，通过记录在步骤A中不同点目标的间距和强度，记录其主瓣位置与主瓣宽度，在信号中根据点目标位置与主瓣的对应关系，保留主瓣部分的信息，去除信号中的旁瓣信息，得到
Figure GDA00027930639400000714
所对应的标签
Figure GDA00027930639400000715
方位向上对应第x距离时刻的距离信号的无旁瓣信号标签记为
Figure GDA00027930639400000716
步骤C2，将成像后的方位信号
Figure GDA00027930639400000717
与对应的标签
Figure GDA00027930639400000718
进行线性插值得到信号
Figure GDA00027930639400000719
与
Figure GDA00027930639400000720
其线性插值公式为：
Figure GDA0002793063940000081
利用上式计算插值中
Figure GDA0002793063940000082
的值，其中m为所需插值的倍数，n为插值后y与y+1间的位置，得到插值之后的信号，以成像后的距离信号
Figure GDA0002793063940000083
为例，其插值得到的信号
Figure GDA0002793063940000084
且
Figure GDA0002793063940000085
SAR回波数据中主瓣宽度太小，卷积神经网络的感受野过大，无法直接提取其中旁瓣与主瓣相关信息，故对SAR数据进行了预插值处理，将主瓣拓展成为大约主瓣内包含N主瓣个像素点的数据。
步骤C3，对插值后的信号
Figure GDA0002793063940000086
与
Figure GDA0002793063940000087
进行滑窗规范化，建立SCNN训练所需距离向输入输出对，如图2所示，为两组信号的滑窗示意，将每一个窗长长度的样本送入网络中进行学习，其中步长为l，窗长为WK窗长。每一个样本可以表达为
Figure GDA0002793063940000088
h为信号采样的各个值，式中WK窗终＝WK窗始+WK窗长，在经历一个步长后，下一个样本可以表达为
Figure GDA0002793063940000089
其所学的值为其样本中心位置处的值
Figure GDA00027930639400000810
卷积神经网络常常通过规范化提高收敛训练速度，SAR回波数据中最高点的强度可达数万，最低仅为0，传统的规范化无法适用差别如此大的数据集，故对每一滑窗样本进行规范化如公式(3)，以减少数据层次过大导致的网络训练模型不佳。
步骤D，采用SCNN搭建网络模型；
建立SCNN网络模型，如附图3所示，该网络模型包含四层卷积层、一层全连接层和一个神经元作为输出。四层卷积层分别是第一卷积层记为conv1、第二卷积层记为conv2、第三卷积层记为conv3和第四卷积层记为conv4。全连接层记为fc。
第一层：卷积层conv1，卷积核大小为30×1，滑动步长为1，对输入图片进行卷积操作，通过零填充使输出10个50×1的特征图到第二卷积层conv2中。
第二层：卷积层conv2，卷积核大小为20×1，滑动步长为1，对输入的特征图进行卷积操作，通过零填充使输出10个50×1的特征图到第三卷积层conv3。
第三层：卷积层conv3，卷积核大小为15×1，滑动步长为1，对输入的特征图进行卷积操作，通过零填充使输出10个50×1的特征图到第四卷积层conv4。
第四层：卷积层conv4，卷积核大小为5×1，滑动步长为1，对输入的特征图进行卷积操作，通过零填充使输出3个50×1的特征图到全连接层fc中。
第五层：全连接层fc，加入Droupout(译文，随机失活算法)，fc输出100个神经元。
通过一个神经元作为整个网络的输出。
在本发明中，用M层表示SCNN网络中的任意一层的标识。
值得一提的是，传统的卷积神经网络通常在某些卷积层后通常会添加一个下采样层，即池化层，抽取卷积得到的特征图中一定范围内的局部平均值或最大值，这样既降低了数据量，也提升了网络对输入数据中目标形变和场景变化的鲁棒性。在本发明中，仿真用的是维度低的二维回波信息，对训练集要求的数量不多，因此不设置池化层，以达到减少过拟合现象。
步骤E，训练SCNN网络模型；
将步骤B得到的训练所需方位向信息和步骤C得到的训练所需距离向信息输入至步骤D构建的SCNN网络模型中，采用随机梯度下降法进行训练，得到一维旁瓣抑制训练模型的网络参数。
在SCNN网络的训练中设置的卷积层的运算为：
Figure GDA0002793063940000091
式中
Figure GDA0002793063940000092
为第M层的输出向量，
Figure GDA0002793063940000093
为卷积核
Figure GDA0002793063940000094
的转置，其中
Figure GDA0002793063940000095
为第M层卷积核，T为转置符号，M层表示层的序号，
Figure GDA0002793063940000096
表示第M层的输入向量，
Figure GDA0002793063940000097
表示第M层卷积层的偏置向量，在本发明中，用非线性函数来表征公式(5)中网络非线性映射能力，则有：
f(Q)＝max(0,Q) (6)
Q代替了公式(5)中的
Figure GDA0002793063940000098
为防止过拟合采用Dropout算法和L2正则化算法。在本发明中，Dropout算法，概率为pr的神经元失活，防止过学习的一种训练策略。L2正则化算法，由杨浩，马建红在2017年08月25日出版的《计算机测量与控制》，“正则化参数求解方法研究”一文，1.2节内容。
为减少训练所需参数量采用Adam优化算法。Adam优化算法，由D.P.Kingma andJ.L.Ba在2015年发布在学习表征国际会议ICLR上的文章Adam:A Method For StochasticOptimization，(译文为Adam,一种随机优化的方法中提供的算法。
在随机梯度下降法中考虑到训练集数据维度匹配采用MSE作为损失函数算法。MSE，由高耀东，侯凌燕，杨大利在2017年01月10日出版的《计算机应用》，“基于多标签学习的卷积神经网络的图像标注方法”一文，1.2节内容中的公式(2)。
本发明的采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法包括有下列步骤：
步骤一，获取SAR实时回波数据；
对从地面接收站实时接收到的回波信号
Figure GDA0002793063940000101
进行成像处理，得到实时SAR回波数据
Figure GDA0002793063940000102
在本发明中，SAR处理器实时接收到的回波信号
Figure GDA0002793063940000103
是二维数据，采用矩阵形式表示为：
Figure GDA0002793063940000104
Figure GDA0002793063940000105
表示实时信号沿方位向上采集到的第1个方位时刻的距离向上的第1个采样点；
Figure GDA0002793063940000106
表示实时信号沿方位向上采集到的第1个方位时刻的距离向上的第2个采样点；
Figure GDA0002793063940000107
表示实时信号沿方位向上采集到的第1个方位时刻的距离向上的第y个采样点；
Figure GDA0002793063940000108
表示实时信号沿方位向上采集到的第1个方位时刻的距离向上的第Nf个采样点；
Figure GDA0002793063940000109
表示实时信号沿方位向上采集到的第2个方位时刻的距离向上的第1个采样点；
Figure GDA00027930639400001010
表示实时信号沿方位向上采集到的第2个方位时刻的距离向上的第2个采样点；
Figure GDA00027930639400001011
表示实时信号沿方位向上采集到的第2个方位时刻的距离向上的第y个采样点；
Figure GDA00027930639400001012
表示实时信号沿方位向上采集到的第2个方位时刻的距离向上的第Nf个采样点；
Figure GDA00027930639400001013
表示实时信号沿方位向上采集到的第x个方位时刻的距离向上的第1个采样点；
Figure GDA00027930639400001014
表示实时信号沿方位向上采集到的第x个方位时刻的距离向上的第2个采样点；
Figure GDA00027930639400001015
表示实时信号沿方位向上采集到的第x个方位时刻的距离向上的第y个采样点；
Figure GDA00027930639400001016
表示实时信号沿方位向上采集到的第x个方位时刻的距离向上的第Nf个采样点；
Figure GDA00027930639400001017
表示实时信号沿方位向上采集到的第Na个方位时刻的距离向上的第1个采样点；
Figure GDA00027930639400001018
表示实时信号沿方位向上采集到的第Na个方位时刻的距离向上的第2个采样点；
Figure GDA00027930639400001019
表示实时信号沿方位向上采集到的第Na个方位时刻的距离向上的第y个采样点；
Figure GDA00027930639400001020
表示实时信号沿方位向上采集到的第Na个方位时刻的距离向上的第Nf个采样点。
在本发明中，采用CSA方法对回波信号矩阵
Figure GDA00027930639400001021
进行成像处理后，得到实时成像结果矩阵
Figure GDA00027930639400001022
的矩阵表示为：
Figure GDA0002793063940000111
Nb是方位向上能够有效产生Na点接收回波所对应的场景目标的采样长度(点数)，亦为观测矩阵
Figure GDA0002793063940000112
的列数；
Nf表示距离向采样点数；
x是一个变量，取值范围1至Nb，表示为成像结果方位向的第x个像素点；
y是一个变量，取值范围1至Nf，表示为成像结果距离向的第y个像素点；
Figure GDA0002793063940000113
表示实时成像结果在方位向上对应第1方位时刻的距离向第1像素点；
Figure GDA0002793063940000114
表示实时成像结果在方位向上对应第1方位时刻的距离向第2像素点；
Figure GDA0002793063940000115
表示实时成像结果在方位向上对应第1方位时刻的距离向第y像素点；
Figure GDA0002793063940000116
表示实时成像结果在方位向上对应第1方位时刻的距离向第Nf像素点；
Figure GDA0002793063940000117
表示实时成像结果在方位向上对应第2方位时刻的距离向第1像素点；
Figure GDA0002793063940000118
表示实时成像结果在方位向上对应第2方位时刻的距离向第2像素点；
Figure GDA0002793063940000119
表示实时成像结果在方位向上对应第2方位时刻的距离向第y像素点；
Figure GDA00027930639400001110
表示实时成像结果在方位向上对应第2方位时刻的距离向第Nf像素点；
Figure GDA00027930639400001111
表示实时成像结果在方位向上对应第x方位时刻的距离向第1像素点；
Figure GDA00027930639400001112
表示实时成像结果在方位向上对应第x方位时刻的距离向第2像素点；
Figure GDA00027930639400001113
表示实时成像结果在方位向上对应第x方位时刻的距离向第y像素点；
Figure GDA00027930639400001114
表示实时成像结果在方位向上对应第x方位时刻的距离向第Nf像素点；
Figure GDA00027930639400001115
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第1像素点；
Figure GDA00027930639400001116
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第2像素点；
Figure GDA00027930639400001117
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第y像素点；
Figure GDA00027930639400001118
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第Nf像素点。
步骤二，方位向上实时回波数据旁瓣抑制处理；
步骤21，取出成像结果矩阵
Figure GDA00027930639400001119
中的距离向上对应第y距离时刻的成像后方位信号
Figure GDA00027930639400001120
记为
Figure GDA00027930639400001121
步骤22，将成像后的方位信号
Figure GDA0002793063940000121
进行线性插值得到方位插值信号
Figure GDA0002793063940000122
其线性插值公式为：
Figure GDA0002793063940000123
利用上式计算插值中
Figure GDA0002793063940000124
的值，其中m为所需插值的倍数，n为插值后相邻方位向间的位置，得到插值之后的信号，以成像后的方位信号
Figure GDA0002793063940000125
为例，其插值得到的信号
Figure GDA0002793063940000126
步骤23，对插值后的信号
Figure GDA0002793063940000127
进行滑窗规范化，建立实时处理回波数据样本，如图2所示，将每一个窗长长度的样本送入网络中进行学习，其中步长为l，窗长为WK窗长。实时采集的第z个样本可以表达为
Figure GDA0002793063940000128
h实时为实时信号采样的各个值，式中WK窗终＝WK窗始+WK窗长，在经历一个步长后，下一个样本可以表达为
Figure GDA0002793063940000129
对每一滑窗样本进行规范化如公式(3)，以减少数据层次过大导致的网络训练模型不佳。
步骤24，将规范化后的样本送入所训练得到的一维旁瓣抑制模型中，得到每一样本的中心时刻预测值，以第z个样本为例，其预测结果为RPREz。
对每一样本规范化进行逆操作如公式(8)得到：
Figure GDA00027930639400001210
步骤25，通过样本寻址方式，将PREz按照第z个样本在原信号的位置重构生成方位向旁瓣抑制后的回波信号
Figure GDA00027930639400001211
对方位向旁瓣抑制后的回波信号
Figure GDA00027930639400001212
进行m的采样，以逆变换预插值操作。
最后将每一方位信号重新组合成方位向旁瓣抑制后的回波数据
Figure GDA00027930639400001213
Figure GDA00027930639400001214
Nb是方位向上能够有效产生Na点接收回波所对应的场景目标的采样长度(点数)，亦为观测矩阵
Figure GDA00027930639400001215
的列数；
Nf表示距离向采样点数；
x是一个变量，取值范围1至Nb，表示为成像结果方位向的第x个像素点；
y是一个变量，取值范围1至Nf，表示为成像结果距离向的第y个像素点；
Figure GDA0002793063940000131
表示方位向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第1像素点；
Figure GDA0002793063940000132
表示方位向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第2像素点；
Figure GDA0002793063940000133
表示方位向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第y像素点；
Figure GDA0002793063940000134
表示方位向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第Nf像素点；
Figure GDA0002793063940000135
表示方位向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第1像素点；
Figure GDA0002793063940000136
表示方位向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第2像素点；
Figure GDA0002793063940000137
表示方位向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第y像素点；
Figure GDA0002793063940000138
表示方位向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第Nf像素点；
Figure GDA0002793063940000139
表示方位向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第1像素点；
Figure GDA00027930639400001310
表示方位向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第2像素点；
Figure GDA00027930639400001311
表示方位向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第y像素点；
Figure GDA00027930639400001312
表示方位向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第Nf像素点；
Figure GDA00027930639400001313
表示方位向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第1像素点；
Figure GDA00027930639400001314
表示方位向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第2像素点；
Figure GDA00027930639400001315
表示方位向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第y像素点；
Figure GDA00027930639400001316
表示方位向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第Nf像素点。
步骤三，距离向上实时回波数据旁瓣抑制处理；
步骤31，取出方位向旁瓣抑制后成像结果矩阵
Figure GDA00027930639400001317
中的距离向上对应第x距离时刻的成像后距离信号，记为
Figure GDA00027930639400001318
且
Figure GDA00027930639400001319
步骤32，将成像后的距离信号
Figure GDA00027930639400001320
进行线性插值得到信号
Figure GDA00027930639400001321
其线性插值公式为
Figure GDA00027930639400001322
利用上式计算插值中
Figure GDA00027930639400001323
的值，其中m为所需插值的倍数，n为插值后y与y+1间的位置，得到插值之后的信号，以距离信号
Figure GDA00027930639400001324
为例，其插值得到的信号
Figure GDA00027930639400001325
步骤33，对插值后的信号
Figure GDA0002793063940000141
进行滑窗规范化，建立实时处理回波数据样本，如图2所示，将每一个窗长长度的样本送入网络中进行学习，其中步长为l，窗长为WK窗长。实时采集的第z个样本可以表达为
Figure GDA0002793063940000142
h实时为实时信号采样的各个值，式中WK窗终＝WK窗始+WK窗长，在经历一个步长后，下一个样本可以表达为
Figure GDA0002793063940000143
对每一滑窗样本进行规范化如公式(3)，以减少数据层次过大导致的网络训练模型不佳。
步骤34，将规范化后的样本送入所训练得到的一维旁瓣抑制模型中，得到每一样本的中心时刻预测值，以第z个样本为例，其预测结果为RPREz。
对每一样本规范化进行逆操作如公式(8)。
步骤35，通过样本寻址方式，将PREz按照第z个样本在原信号的的位置重构生成距离向旁瓣抑制后的回波信号
Figure GDA0002793063940000144
对距离向旁瓣抑制后的回波信号
Figure GDA0002793063940000145
进行m的采样，以逆变换预插值操作。
最后将每一距离信号重新组合成旁瓣抑制后的回波数据
Figure GDA0002793063940000146
Figure GDA0002793063940000147
Nb是方位向上能够有效产生Na点接收回波所对应的场景目标的采样长度(点数)，亦为观测矩阵
Figure GDA0002793063940000148
的列数；
Nf表示距离向采样点数；
x是一个变量，取值范围1至Nb，表示为成像结果方位向的第x个像素点；
y是一个变量，取值范围1至Nf，表示为成像结果距离向的第y个像素点；
Figure GDA0002793063940000149
表示距离向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第1像素点；
Figure GDA00027930639400001410
表示距离向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第2像素点；
Figure GDA00027930639400001411
表示距离向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第y像素点；
Figure GDA00027930639400001412
表示距离向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第Nf像素点；
Figure GDA00027930639400001413
表示距离向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第1像素点；
Figure GDA0002793063940000151
表示距离向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第2像素点；
Figure GDA0002793063940000152
表示距离向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第y像素点；
Figure GDA0002793063940000153
表示距离向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第Nf像素点；
Figure GDA0002793063940000154
表示距离向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第1像素点；
Figure GDA0002793063940000155
表示距离向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第2像素点；
Figure GDA0002793063940000156
表示距离向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第y像素点；
Figure GDA0002793063940000157
表示距离向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第Nf像素点；
Figure GDA0002793063940000158
表示距离向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第1像素点；
Figure GDA0002793063940000159
表示距离向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第2像素点；
Figure GDA00027930639400001510
表示距离向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第y像素点；
Figure GDA00027930639400001511
表示距离向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第Nf像素点。
步骤四，RSVA处理；
将步骤一的实时成像结果矩阵
Figure GDA00027930639400001512
采用RSVA处理得到RSVA抑制后的结果矩阵
Figure GDA00027930639400001513
Figure GDA00027930639400001514
Nb是方位向上能够有效产生所述Na点接收回波所对应的场景目标的采样长度(即点数)，亦为观测矩阵
Figure GDA00027930639400001515
的列数；
Nf表示距离向采样点总数；
x是一个变量，取值范围1至Nb，表示为成像结果方位向的第x个像素点；
y是一个变量，取值范围1至Nf，表示为成像结果距离向的第y个像素点；
Figure GDA00027930639400001516
表示实时成像结果在方位向上对应第1方位时刻的距离向第1像素点；
Figure GDA00027930639400001517
表示实时成像结果在方位向上对应第1方位时刻的距离向第2像素点；
Figure GDA00027930639400001518
表示实时成像结果在方位向上对应第1方位时刻的距离向第y像素点；
Figure GDA00027930639400001519
表示实时成像结果在方位向上对应第1方位时刻的距离向第Nf像素点；
Figure GDA00027930639400001520
表示实时成像结果在方位向上对应第2方位时刻的距离向第1像素点；
Figure GDA0002793063940000161
表示实时成像结果在方位向上对应第2方位时刻的距离向第2像素点；
Figure GDA0002793063940000162
表示实时成像结果在方位向上对应第2方位时刻的距离向第y像素点；
Figure GDA0002793063940000163
表示实时成像结果在方位向上对应第2方位时刻的距离向第Nf像素点；
Figure GDA0002793063940000164
表示实时成像结果在方位向上对应第x方位时刻的距离向第1像素点；
Figure GDA0002793063940000165
表示实时成像结果在方位向上对应第x方位时刻的距离向第2像素点；
Figure GDA0002793063940000166
表示实时成像结果在方位向上对应第x方位时刻的距离向第y像素点；
Figure GDA0002793063940000167
表示实时成像结果在方位向上对应第x方位时刻的距离向第Nf像素点；
Figure GDA0002793063940000168
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第1像素点；
Figure GDA0002793063940000169
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第2像素点；
Figure GDA00027930639400001610
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第y像素点；
Figure GDA00027930639400001611
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第Nf像素点。
步骤五，RSVA与SCNN组合处理；
考虑到RSVA在旁瓣抑制处理中，其结果在旁瓣部分已达到理想，而其主瓣部分丢失的能量将会影响目标在信噪比较高图像中的目标解译，需要对其主瓣部分丢失能量进行补偿。
在SCNN中，其学习内容为回波数据，忽略了其相位信息，所以对于主瓣部分能量无损失，但是因为网络感受野及其数据驱动下造成忽略回波数据物理影响，其旁瓣区域的抑制能量略逊于RSVA，故需要对其旁瓣部分进行更有效的抑制。
将步骤四中RSVA处理得到RSVA抑制后的结果矩阵的
Figure GDA00027930639400001612
与步骤三重新组合成旁瓣抑制后的回波数据
Figure GDA00027930639400001613
和步骤一实时成像结果矩阵
Figure GDA00027930639400001614
用于补偿所述
Figure GDA00027930639400001615
的主瓣部分，主瓣补偿关系是为得到抑制后的回波数据Efinal，且
Figure GDA00027930639400001616
在本发明中，所述抑制后的回波数据Efinal有效解决了RSVA中主瓣能量丢失和SCNN中旁瓣部分残余能量的技术问题。
实施例1
本实施例提出一种基于卷积神经网络的合成孔径雷达旁瓣抑制算法，其中成像过程中涉及的参数如表1所示。
表1实施例参数
参数 取值
雷达平台等效速度v(m/s) 179
参考斜距R<sub>0</sub>(m) 9867
脉冲重复频率f<sub>p</sub>(Hz) 1000
距离向采样率f<sub>s</sub>(MHz) 4e8
发射信号波长λ<sub>0</sub>(m) 0.05
斜视角θ(°) 0
构建一维旁瓣抑制训练模型
对仿真的回波数据集通过频域插值，将主瓣宽度插值到倍数m，本实施例中插值大小为m＝10，主瓣宽度变为30。
在图4中，滑窗步长设为l＝1，窗长设为WK窗长＝50，对回波数据的方位向和距离向进行滑窗遍历得到所有训练集样本。
对每一个滑窗样本进行数据规范化，见公式(3)，其中Amax＝1000。
建立卷积神经网络，如图3所示，该网络模型包含四层卷积层、一层全连接层和一个神经元作为输出。
将得到的所有训练集样本共400个约2048×1的向量输入卷积神经网络模型中，采用随机梯度下降法对网络模型进行训练，网络学习率设为0.001，随迭代次数递增指数递减，迭代次数为300，得到一维旁瓣抑制训练模型。其中卷积层按公式(5)、(6)和L2正则化算法的计算；损失函数为均方误差，对损失函数采用Adam随机梯度下降模式训练网络，得到一维旁瓣抑制训练模型。
实时处理SAR回波数据
利用图2的处理流程，并结合表1的参数得到回波信号，进行CSA成像后得到回波数据，对回波数据进行预插值滑窗规范化操作后，对回波数据进行RSVA和RSVA与SCNN组合的处理，得到训练结果如图5～图16所示。
参见图5、图6、图7、图8、图9和图10，通过六幅图对比能够看出本发明方法与RSVA处理均对单点目标回波数据具有较佳的旁瓣抑制效果。RSVA和本发明的方法均对单点目标旁瓣抑制后的回波数据峰值旁瓣比在－37dB左右。
参见图11、图12、图13、图14、图15和图16，通过六幅图对比能够看出，RSVA算法中由于相位因子导致的多点目标回波数据能量丢失在组合算法中有效补偿回来，这可以解决目标在信噪比较高情况下的解译问题。对RSVA进行多点目标旁瓣抑制中主瓣能量损失由50％减少至15％。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (5)

1.一种采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法，其特征在于包括有下列步骤：
步骤一，获取SAR实时回波数据；
对从地面接收站实时接收到的回波信号
Figure FDA0002793063930000011
进行成像处理，得到实时SAR回波数据
Figure FDA0002793063930000012
所述
Figure FDA0002793063930000013
是二维数据，采用矩阵形式表示为：
Figure FDA0002793063930000014
Na表示方位向采样点总数；
Nf表示距离向采样点总数；
x是一个变量，取值范围1至Na，表示为方位向第x个点，也称为方位向的标识号；
y是一个变量，取值范围1至Nf，表示为距离向第y个点，也称为距离向的标识号；
Figure FDA0002793063930000015
表示实时信号沿方位向上采集到的第1个方位时刻的距离向上的第1个采样点；
Figure FDA0002793063930000016
表示实时信号沿方位向上采集到的第1个方位时刻的距离向上的第2个采样点；
Figure FDA0002793063930000017
表示实时信号沿方位向上采集到的第1个方位时刻的距离向上的第y个采样点；
Figure FDA0002793063930000018
表示实时信号沿方位向上采集到的第1个方位时刻的距离向上的第Nf个采样点；
Figure FDA0002793063930000019
表示实时信号沿方位向上采集到的第2个方位时刻的距离向上的第1个采样点；
Figure FDA00027930639300000110
表示实时信号沿方位向上采集到的第2个方位时刻的距离向上的第2个采样点；
Figure FDA00027930639300000111
表示实时信号沿方位向上采集到的第2个方位时刻的距离向上的第y个采样点；
Figure FDA00027930639300000112
表示实时信号沿方位向上采集到的第2个方位时刻的距离向上的第Nf个采样点；
Figure FDA00027930639300000113
表示实时信号沿方位向上采集到的第x个方位时刻的距离向上的第1个采样点；
Figure FDA0002793063930000021
表示实时信号沿方位向上采集到的第x个方位时刻的距离向上的第2个采样点；
Figure FDA0002793063930000022
表示实时信号沿方位向上采集到的第x个方位时刻的距离向上的第y个采样点；
Figure FDA0002793063930000023
表示实时信号沿方位向上采集到的第x个方位时刻的距离向上的第Nf个采样点；
Figure FDA0002793063930000024
表示实时信号沿方位向上采集到的第Na个方位时刻的距离向上的第1个采样点；
Figure FDA0002793063930000025
表示实时信号沿方位向上采集到的第Na个方位时刻的距离向上的第2个采样点；
Figure FDA0002793063930000026
表示实时信号沿方位向上采集到的第Na个方位时刻的距离向上的第y个采样点；
Figure FDA0002793063930000027
表示实时信号沿方位向上采集到的第Na个方位时刻的距离向上的第Nf个采样点；
采用CSA方法对回波信号矩阵
Figure FDA0002793063930000028
进行成像处理后，得到实时成像结果矩阵
Figure FDA0002793063930000029
的矩阵表示为：
Figure FDA00027930639300000210
Nb是方位向上能够有效产生Na点接收回波所对应的场景目标的采样长度，即点数，亦为观测矩阵
Figure FDA00027930639300000211
的列数；
Figure FDA00027930639300000212
表示实时成像结果在方位向上对应第1方位时刻的距离向第1像素点；
Figure FDA00027930639300000213
表示实时成像结果在方位向上对应第1方位时刻的距离向第2像素点；
Figure FDA00027930639300000214
表示实时成像结果在方位向上对应第1方位时刻的距离向第y像素点；
Figure FDA00027930639300000215
表示实时成像结果在方位向上对应第1方位时刻的距离向第Nf像素点；
Figure FDA00027930639300000216
表示实时成像结果在方位向上对应第2方位时刻的距离向第1像素点；
Figure FDA00027930639300000217
表示实时成像结果在方位向上对应第2方位时刻的距离向第2像素点；
Figure FDA00027930639300000218
表示实时成像结果在方位向上对应第2方位时刻的距离向第y像素点；
Figure FDA0002793063930000031
表示实时成像结果在方位向上对应第2方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000032
表示实时成像结果在方位向上对应第x方位时刻的距离向第1像素点；
Figure FDA0002793063930000033
表示实时成像结果在方位向上对应第x方位时刻的距离向第2像素点；
Figure FDA0002793063930000034
表示实时成像结果在方位向上对应第x方位时刻的距离向第y像素点；
Figure FDA0002793063930000035
表示实时成像结果在方位向上对应第x方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000036
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第1像素点；
Figure FDA0002793063930000037
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第2像素点；
Figure FDA0002793063930000038
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第y像素点；
Figure FDA0002793063930000039
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第Nf像素点；
步骤二，方位向上实时回波数据旁瓣抑制处理；
步骤21，取出成像结果矩阵
Figure FDA00027930639300000310
中的距离向上对应第y距离时刻的成像后方位信号
Figure FDA00027930639300000311
记为
Figure FDA00027930639300000312
步骤22，将成像后的方位信号
Figure FDA00027930639300000313
进行线性插值得到方位插值信号
Figure FDA00027930639300000314
其线性插值公式为：
Figure FDA00027930639300000315
利用上式计算插值中
Figure FDA00027930639300000316
的值，其中m为所需插值的倍数，n为插值后相邻方位向间的位置，得到插值之后的信号，以成像后的方位信号
Figure FDA00027930639300000317
为例，其插值得到的信号
Figure FDA00027930639300000318
步骤23，对插值后的信号
Figure FDA00027930639300000319
进行滑窗规范化，建立实时处理回波数据样本，将每一个窗长长度的样本送入网络中进行学习，其中步长为l，窗长为WK窗长；实时采集的第z个样本表达为
Figure FDA00027930639300000320
h实时为实时信号采样的各个值，式中WK窗终＝WK窗始+WK窗长，在经历一个步长后，下一个样本表达为
Figure FDA0002793063930000041
对每一滑窗样本进行规范化
Figure FDA0002793063930000042
处理，以减少数据层次过大导致的网络训练模型不佳；式中Amax为考虑到数据存储可能造成的丢失所设置的规范化值，使输出范围为[0,Amax]，hp为样本I里每一个像素点的数值，hmax为每一样本I里的最大值，hp规范为规范化后的值；
步骤24，将规范化后的样本送入所训练得到的一维旁瓣抑制模型中，得到每一样本的中心时刻预测值，以第z个样本为例，其预测结果为RPREz；
对每一样本规范化进行逆操作得到：
Figure FDA0002793063930000043
步骤25，通过样本寻址方式，将PREz按照第z个样本在原信号的位置重构生成方位向旁瓣抑制后的回波信号
Figure FDA0002793063930000044
对方位向旁瓣抑制后的回波信号
Figure FDA0002793063930000045
进行m的采样，以逆变换预插值操作；
最后将每一方位信号重新组合成方位向旁瓣抑制后的回波数据
Figure FDA0002793063930000046
Figure FDA0002793063930000047
Nb是方位向上能够有效产生Na点接收回波所对应的场景目标的采样长度，即点数，亦为观测矩阵
Figure FDA0002793063930000048
的列数；
Figure FDA0002793063930000049
表示方位向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第1像素点；
Figure FDA00027930639300000410
表示方位向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第2像素点；
Figure FDA00027930639300000411
表示方位向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第y像素点；
Figure FDA00027930639300000412
表示方位向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000051
表示方位向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第1像素点；
Figure FDA0002793063930000052
表示方位向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第2像素点；
Figure FDA0002793063930000053
表示方位向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第y像素点；
Figure FDA0002793063930000054
表示方位向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000055
表示方位向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第1像素点；
Figure FDA0002793063930000056
表示方位向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第2像素点；
Figure FDA0002793063930000057
表示方位向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第y像素点；
Figure FDA0002793063930000058
表示方位向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000059
表示方位向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第1像素点；
Figure FDA00027930639300000510
表示方位向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第2像素点；
Figure FDA00027930639300000511
表示方位向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第y像素点；
Figure FDA00027930639300000512
表示方位向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第Nf像素点；
步骤三，距离向上实时回波数据旁瓣抑制处理；
步骤31，取出方位向旁瓣抑制后成像结果矩阵
Figure FDA00027930639300000513
中的距离向上对应第x距离时刻的成像后距离信号，记为
Figure FDA00027930639300000514
且
Figure FDA00027930639300000515
步骤32，将成像后的距离信号
Figure FDA00027930639300000516
进行线性插值得到信号
Figure FDA00027930639300000517
其线性插值公式为
Figure FDA00027930639300000518
利用上式计算插值中
Figure FDA00027930639300000519
的值，其中m为所需插值的倍数，n为插值后y与y+1间的位置，得到插值之后的信号，以距离信号
Figure FDA00027930639300000520
为例，其插值得到的信号
Figure FDA00027930639300000521
步骤33，对插值后的信号
Figure FDA00027930639300000522
进行滑窗规范化，建立实时处理回波数据样本，将每一个窗长长度的样本送入网络中进行学习，其中步长为l，窗长为WK窗长；实时采集的第z个样本表达为
Figure FDA0002793063930000061
h实时为实时信号采样的各个值，式中WK窗终＝WK窗始+WK窗长，在经历一个步长后，下一个样本表达为
Figure FDA0002793063930000062
对每一滑窗样本进行规范化
Figure FDA0002793063930000063
处理，以减少数据层次过大导致的网络训练模型不佳；
步骤34，将规范化后的样本送入所训练得到的一维旁瓣抑制模型中，得到每一样本的中心时刻预测值，以第z个样本为例，其预测结果为RPREz；
对每一样本规范化进行逆操作
Figure FDA0002793063930000064
处理；
步骤35，通过样本寻址方式，将PREz按照第z个样本在原信号的的位置重构生成距离向旁瓣抑制后的回波信号
Figure FDA0002793063930000065
对距离向旁瓣抑制后的回波信号
Figure FDA0002793063930000066
进行m的采样，以逆变换预插值操作；
最后将每一距离信号重新组合成旁瓣抑制后的回波数据
Figure FDA0002793063930000067
Figure FDA0002793063930000068
Nb是方位向上能够有效产生Na点接收回波所对应的场景目标的采样长度，即点数，亦为观测矩阵
Figure FDA0002793063930000069
的列数；
Figure FDA00027930639300000610
表示距离向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第1像素点；
Figure FDA00027930639300000611
表示距离向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第2像素点；
Figure FDA00027930639300000612
表示距离向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第y像素点；
Figure FDA0002793063930000071
表示距离向旁瓣抑制后成像结果在方位向上对应第1方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000072
表示距离向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第1像素点；
Figure FDA0002793063930000073
表示距离向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第2像素点；
Figure FDA0002793063930000074
表示距离向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第y像素点；
Figure FDA0002793063930000075
表示距离向旁瓣抑制后成像结果在方位向上对应第2方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000076
表示距离向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第1像素点；
Figure FDA0002793063930000077
表示距离向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第2像素点；
Figure FDA0002793063930000078
表示距离向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第y像素点；
Figure FDA0002793063930000079
表示距离向旁瓣抑制后成像结果在方位向上对应第x方位时刻的距离向第Nf像素点；
Figure FDA00027930639300000710
表示距离向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第1像素点；
Figure FDA00027930639300000711
表示距离向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第2像素点；
Figure FDA00027930639300000712
表示距离向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第y像素点；
Figure FDA00027930639300000713
表示距离向旁瓣抑制后成像结果在方位向上对应第Nb方位时刻的距离向第Nf像素点；
步骤四，RSVA处理；
将步骤一的实时成像结果矩阵
Figure FDA00027930639300000714
采用RSVA处理得到RSVA抑制后的结果矩阵
Figure FDA00027930639300000715
Figure FDA00027930639300000716
Nb是方位向上能够有效产生所述Na点接收回波所对应的场景目标的采样长度，即点数，亦为观测矩阵
Figure FDA00027930639300000717
的列数；
Figure FDA0002793063930000081
表示实时成像结果在方位向上对应第1方位时刻的距离向第1像素点；
Figure FDA0002793063930000082
表示实时成像结果在方位向上对应第1方位时刻的距离向第2像素点；
Figure FDA0002793063930000083
表示实时成像结果在方位向上对应第1方位时刻的距离向第y像素点；
Figure FDA0002793063930000084
表示实时成像结果在方位向上对应第1方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000085
表示实时成像结果在方位向上对应第2方位时刻的距离向第1像素点；
Figure FDA0002793063930000086
表示实时成像结果在方位向上对应第2方位时刻的距离向第2像素点；
Figure FDA0002793063930000087
表示实时成像结果在方位向上对应第2方位时刻的距离向第y像素点；
Figure FDA0002793063930000088
表示实时成像结果在方位向上对应第2方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000089
表示实时成像结果在方位向上对应第x方位时刻的距离向第1像素点；
Figure FDA00027930639300000810
表示实时成像结果在方位向上对应第x方位时刻的距离向第2像素点；
Figure FDA00027930639300000811
表示实时成像结果在方位向上对应第x方位时刻的距离向第y像素点；
Figure FDA00027930639300000812
表示实时成像结果在方位向上对应第x方位时刻的距离向第Nf像素点；
Figure FDA00027930639300000813
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第1像素点；
Figure FDA00027930639300000814
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第2像素点；
Figure FDA00027930639300000815
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第y像素点；
Figure FDA00027930639300000816
表示实时成像结果在方位向上对应第Nb方位时刻的距离向第Nf像素点；
步骤五，RSVA与SCNN组合处理；
将步骤四中RSVA处理得到RSVA抑制后的结果矩阵的
Figure FDA00027930639300000817
与步骤三重新组合成旁瓣抑制后的回波数据
Figure FDA00027930639300000818
和步骤一实时成像结果矩阵
Figure FDA00027930639300000819
用于补偿所述
Figure FDA00027930639300000820
的主瓣部分，主瓣补偿关系是为得到抑制后的回波数据Efinal，且
Figure FDA00027930639300000821
2.根据权利要求1所述的采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法，其特征在于：在步骤二和步骤三中使用了一维旁瓣抑制训练模型，所述一维旁瓣抑制训练模型的构建步骤为：
步骤A，获取SAR仿真回波数据；
仿真输出的二维仿真SAR回波信号，记为
Figure FDA0002793063930000091
并对所述
Figure FDA0002793063930000092
进行成像处理，得到SAR仿真回波数据
Figure FDA0002793063930000093
所述
Figure FDA0002793063930000094
是二维数据，采用矩阵形式表示为：
Figure FDA0002793063930000095
Na表示方位向采样点总数；
Nf表示距离向采样点总数；
x是一个变量，取值范围1至Na，表示为方位向第x个点，也称为方位向的标识号；
y是一个变量，取值范围1至Nf，表示为距离向第y个点，也称为距离向的标识号；
Figure FDA0002793063930000096
表示仿真信号沿方位向上采集到的第1个方位时刻的距离向上的第1个采样点；
Figure FDA0002793063930000097
表示仿真信号沿方位向上采集到的第1个方位时刻的距离向上的第2个采样点；
Figure FDA0002793063930000098
表示仿真信号沿方位向上采集到的第1个方位时刻的距离向上的第y个采样点；
Figure FDA0002793063930000099
表示仿真信号沿方位向上采集到的第1个方位时刻的距离向上的第Nf个采样点；
Figure FDA00027930639300000910
表示仿真信号沿方位向上采集到的第2个方位时刻的距离向上的第1个采样点；
Figure FDA00027930639300000911
表示仿真信号沿方位向上采集到的第2个方位时刻的距离向上的第2个采样点；
Figure FDA00027930639300000912
表示仿真信号沿方位向上采集到的第2个方位时刻的距离向上的第y个采样点；
Figure FDA00027930639300000913
表示仿真信号沿方位向上采集到的第2个方位时刻的距离向上的第Nf个采样点；
Figure FDA00027930639300000914
表示仿真信号沿方位向上采集到的第x个方位时刻的距离向上的第1个采样点；
Figure FDA0002793063930000101
表示仿真信号沿方位向上采集到的第x个方位时刻的距离向上的第2个采样点；
Figure FDA0002793063930000102
表示仿真信号沿方位向上采集到的第x个方位时刻的距离向上的第y个采样点；
Figure FDA0002793063930000103
表示仿真信号沿方位向上采集到的第x个方位时刻的距离向上的第Nf个采样点；
Figure FDA0002793063930000104
表示仿真信号沿方位向上采集到的第Na个方位时刻的距离向上的第1个采样点；
Figure FDA0002793063930000105
表示仿真信号沿方位向上采集到的第Na个方位时刻的距离向上的第2个采样点；
Figure FDA0002793063930000106
表示仿真信号沿方位向上采集到的第Na个方位时刻的距离向上的第y个采样点；
Figure FDA0002793063930000107
表示仿真信号沿方位向上采集到的第Na个方位时刻的距离向上的第Nf个采样点；
采用CSA方法对仿真的回波信号矩阵
Figure FDA0002793063930000108
进行成像处理后，得到仿真的成像结果矩阵
Figure FDA0002793063930000109
的矩阵表示为：
Figure FDA00027930639300001010
Nb是方位向上能够有效产生所述Na点接收回波所对应的场景目标的采样长度，即点数，亦为观测矩阵
Figure FDA00027930639300001011
的列数；
Nf表示距离向采样点总数；
x是一个变量，取值范围1至Nb，表示为成像结果方位向的第x个像素点，也称为成像结果方位向的标识号；
y是一个变量，取值范围1至Nf，表示为成像结果距离向的第y个像素点，也称为成像结果距离向的标识号；
Figure FDA00027930639300001012
表示仿真成像结果在方位向上对应第1方位时刻的距离向第1像素点；
Figure FDA00027930639300001013
表示仿真成像结果在方位向上对应第1方位时刻的距离向第2像素点；
Figure FDA0002793063930000111
表示仿真成像结果在方位向上对应第1方位时刻的距离向第y像素点；
Figure FDA0002793063930000112
表示仿真成像结果在方位向上对应第1方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000113
表示仿真成像结果在方位向上对应第2方位时刻的距离向第1像素点；
Figure FDA0002793063930000114
表示仿真成像结果在方位向上对应第2方位时刻的距离向第2像素点；
Figure FDA0002793063930000115
表示仿真成像结果在方位向上对应第2方位时刻的距离向第y像素点；
Figure FDA0002793063930000116
表示仿真成像结果在方位向上对应第2方位时刻的距离向第Nf像素点；
Figure FDA0002793063930000117
表示仿真成像结果在方位向上对应第x方位时刻的距离向第1像素点；
Figure FDA0002793063930000118
表示仿真成像结果在方位向上对应第x方位时刻的距离向第2像素点；
Figure FDA0002793063930000119
表示仿真成像结果在方位向上对应第x方位时刻的距离向第y像素点；
Figure FDA00027930639300001110
表示仿真成像结果在方位向上对应第x方位时刻的距离向第Nf像素点；
Figure FDA00027930639300001111
表示仿真成像结果在方位向上对应第Nb方位时刻的距离向第1像素点；
Figure FDA00027930639300001112
表示仿真成像结果在方位向上对应第Nb方位时刻的距离向第2像素点；
Figure FDA00027930639300001113
表示仿真成像结果在方位向上对应第Nb方位时刻的距离向第y像素点；
Figure FDA00027930639300001114
表示仿真成像结果在方位向上对应第Nb方位时刻的距离向第Nf像素点；
步骤B，方位向上无旁瓣的信号集获取；
步骤B1，取出成像结果矩阵
Figure FDA00027930639300001115
中的距离向上对应第y距离时刻的成像后方位信号，记为
Figure FDA00027930639300001116
且
Figure FDA00027930639300001117
其成像信号满足形式：
Figure FDA00027930639300001118
其中，g(t)是所接收的回波，N是场景中的所有目标的个数，i是点目标的计数，σi是第i点目标的后向散射系数的赋值，Pi是第i个点目标的点扩散函数，t是信号所对应的时间，τ是信号的时延，e是对数符号，
Figure FDA00027930639300001119
中上角标的j是虚部，θi是第i个散射单元造成的随机相位；
成像后的方位信号
Figure FDA0002793063930000121
为所设置的SAR系统及信号处理器中的各项参数中不同组合的地面接收站的各个反射系数的信号，通过记录在步骤A中不同点目标的间距和强度，记录其主瓣位置与主瓣宽度，在信号中根据点目标位置与主瓣的对应关系，保留主瓣部分的信息，去除信号中的旁瓣信息，得到
Figure FDA0002793063930000122
所对应的标签
Figure FDA0002793063930000123
距离向上对应第x距离时刻的方位信号的无旁瓣信号标签记为
Figure FDA0002793063930000124
步骤B2，将成像后的方位信号
Figure FDA0002793063930000125
与对应的标签
Figure FDA0002793063930000126
进行线性插值得到信号
Figure FDA0002793063930000127
与
Figure FDA0002793063930000128
其线性插值公式为
Figure FDA0002793063930000129
利用上式计算插值中
Figure FDA00027930639300001210
的值，其中m为所需插值的倍数，n为插值后相邻方位向间的位置，得到插值之后的信号，以成像后的方位信号
Figure FDA00027930639300001211
为例，其插值得到的信号
Figure FDA00027930639300001212
SAR回波数据中主瓣宽度太小，卷积神经网络的感受野过大，无法直接提取其中旁瓣与主瓣相关信息，故对SAR数据进行了预插值处理，将主瓣拓展成为大约主瓣内包含N主瓣个像素点的数据；
步骤B3，对插值后的信号
Figure FDA00027930639300001213
与
Figure FDA00027930639300001214
进行滑窗规范化，建立SCNN训练所需方位向输入输出对，为两组信号的滑窗示意，将每一个窗长长度的样本送入网络中进行学习，其中步长为l，窗长为WK窗长；每一个样本表达为
Figure FDA00027930639300001215
h为信号采样的各个值，式中WK窗终＝WK窗始+WK窗长，在经历一个步长后，下一个样本表达为
Figure FDA00027930639300001216
其所学的值为其样本中心位置处的值
Figure FDA00027930639300001217
卷积神经网络常常通过规范化提高收敛训练速度，SAR回波数据中最高点的强度可达数万，最低仅为0，传统的规范化无法适用差别如此大的数据集，故对每一滑窗样本进行规范化如公式(3)，以减少数据层次过大导致的网络训练模型不佳；
Figure FDA0002793063930000131
式中Amax为考虑到数据存储可能造成的丢失所设置的规范化值，使输出范围为[0,Amax]，hp为样本I里每一个像素点的数值，hmax为每一样本I里的最大值，hp规范为规范化后的值；
步骤C，距离向上无旁瓣的信号标签获取；
步骤C1，取出成像结果矩阵
Figure FDA0002793063930000132
中的方位向上对应第x方位时刻的距离信号，记为
Figure FDA0002793063930000133
且
Figure FDA0002793063930000134
其成像信号满足公式(1)；成像后的距离信号
Figure FDA0002793063930000135
为所设置的SAR系统及信号处理器中的各项参数中不同组合的地面接收站的各个反射系数的信号，通过记录在步骤A中不同点目标的间距和强度，记录其主瓣位置与主瓣宽度，在信号中根据点目标位置与主瓣的对应关系，保留主瓣部分的信息，去除信号中的旁瓣信息，得到
Figure FDA0002793063930000136
所对应的标签
Figure FDA0002793063930000137
方位向上对应第x距离时刻的距离信号的无旁瓣信号标签记为
Figure FDA0002793063930000138
步骤C2，将成像后的方位信号
Figure FDA0002793063930000139
与对应的标签
Figure FDA00027930639300001310
进行线性插值得到信号
Figure FDA00027930639300001311
与
Figure FDA00027930639300001312
其线性插值公式为：
Figure FDA00027930639300001313
利用上式计算插值中
Figure FDA00027930639300001314
的值，其中m为所需插值的倍数，n为插值后y与y+1间的位置，得到插值之后的信号，以成像后的距离信号
Figure FDA00027930639300001315
为例，其插值得到的信号
Figure FDA00027930639300001316
且
Figure FDA00027930639300001317
SAR回波数据中主瓣宽度太小，卷积神经网络的感受野过大，无法直接提取其中旁瓣与主瓣相关信息，故对SAR数据进行了预插值处理，将主瓣拓展成为大约主瓣内包含N主瓣个像素点的数据；
步骤C3，对插值后的信号
Figure FDA00027930639300001318
与
Figure FDA00027930639300001319
进行滑窗规范化，建立SCNN训练所需距离向输入输出对，为两组信号的滑窗示意，将每一个窗长长度的样本送入网络中进行学习，其中步长为l，窗长为WK窗长；每一个样本表达为
Figure FDA0002793063930000141
h为信号采样的各个值，式中WK窗终＝WK窗始+WK窗长，在经历一个步长后，下一个样本表达为
Figure FDA0002793063930000142
其所学的值为其样本中心位置处的值
Figure FDA0002793063930000143
卷积神经网络常常通过规范化提高收敛训练速度，SAR回波数据中最高点的强度可达数万，最低仅为0，传统的规范化无法适用差别如此大的数据集，故对每一滑窗样本进行规范化如公式(3)，以减少数据层次过大导致的网络训练模型不佳；
步骤D，采用SCNN搭建网络模型；
建立SCNN网络模型，该网络模型包含四层卷积层、一层全连接层和一个神经元作为输出；四层卷积层分别是第一卷积层记为conv1、第二卷积层记为conv2、第三卷积层记为conv3和第四卷积层记为conv4；全连接层记为fc；
第一层：卷积层conv1，卷积核大小为30×1，滑动步长为1，对输入图片进行卷积操作，通过零填充使输出10个50×1的特征图到第二卷积层conv2中；
第二层：卷积层conv2，卷积核大小为20×1，滑动步长为1，对输入的特征图进行卷积操作，通过零填充使输出10个50×1的特征图到第三卷积层conv3；
第三层：卷积层conv3，卷积核大小为15×1，滑动步长为1，对输入的特征图进行卷积操作，通过零填充使输出10个50×1的特征图到第四卷积层conv4；
第四层：卷积层conv4，卷积核大小为5×1，滑动步长为1，对输入的特征图进行卷积操作，通过零填充使输出3个50×1的特征图到全连接层fc中；
第五层：全连接层fc，加入Droupout，即随机失活算法，fc输出100个神经元；
通过一个神经元作为整个网络的输出；
步骤E，训练SCNN网络模型；
将步骤B得到的训练所需方位向信息和步骤C得到的训练所需距离向信息输入至步骤D构建的SCNN网络模型中，采用随机梯度下降法进行训练，得到一维旁瓣抑制训练模型的网络参数；
在SCNN网络的训练中设置的卷积层的运算为：
Figure FDA0002793063930000151
式中
Figure FDA0002793063930000152
为第M层的输出向量，
Figure FDA0002793063930000153
为卷积核
Figure FDA0002793063930000154
的转置，其中
Figure FDA0002793063930000155
为第M层卷积核，T为转置符号，M层表示层的序号，
Figure FDA0002793063930000156
表示第M层的输入向量，
Figure FDA0002793063930000157
表示第M层卷积层的偏置向量，用非线性函数来表征公式(5)中网络非线性映射能力，则有：
f(Q)＝max(0,Q) (6)
Q代替了公式(5)中的
Figure FDA0002793063930000158
为防止过拟合采用Dropout算法和L2正则化算法；
为减少训练所需参数量采用Adam优化算法；
在随机梯度下降法中考虑到训练集数据维度匹配采用MSE作为损失函数算法。
3.根据权利要求2所述的采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法，其特征在于：仿真用的是维度低的二维回波信息，对训练集要求的数量不多，因此不设置池化层，以达到减少过拟合现象。
4.根据权利要求1所述的采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法，其特征在于：对单点目标旁瓣抑制后的回波数据峰值旁瓣比在－37dB左右。
5.根据权利要求1所述的采用CNN与RSVA组合对SAR回波数据旁瓣抑制的处理方法，其特征在于：对RSVA进行多点目标旁瓣抑制中主瓣能量损失由50％减少至15％。
CN201910673129.1A 2019-07-24 2019-07-24 采用cnn与rsva组合对sar回波数据旁瓣抑制的处理方法 Active CN110333489B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910673129.1A CN110333489B (zh) 2019-07-24 2019-07-24 采用cnn与rsva组合对sar回波数据旁瓣抑制的处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910673129.1A CN110333489B (zh) 2019-07-24 2019-07-24 采用cnn与rsva组合对sar回波数据旁瓣抑制的处理方法

Publications (2)

Publication Number Publication Date
CN110333489A CN110333489A (zh) 2019-10-15
CN110333489B true CN110333489B (zh) 2021-03-09

Family

ID=68147398

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910673129.1A Active CN110333489B (zh) 2019-07-24 2019-07-24 采用cnn与rsva组合对sar回波数据旁瓣抑制的处理方法

Country Status (1)

Country Link
CN (1) CN110333489B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111126199B (zh) * 2019-12-11 2023-05-30 复旦大学 基于回波测量数据的信号特征提取与数据挖掘方法
CN112213703B (zh) * 2020-08-18 2023-07-18 成都信息工程大学 一种利用毫米波云雷达反演云参数的方法及装置
CN114966560B (zh) * 2022-07-29 2022-10-28 中南大学 一种探地雷达后向投影成像方法及系统
CN115685096B (zh) * 2022-12-30 2023-03-21 成都实时技术股份有限公司 一种基于逻辑回归的二次雷达副瓣抑制方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080024359A1 (en) * 2006-07-25 2008-01-31 Harris Corporation System and method for geometric apodization
CN102721965B (zh) * 2012-06-08 2013-08-21 北京航空航天大学 用于合成孔径雷达成像处理的全向旁瓣抑制方法
CN106680785B (zh) * 2017-03-06 2019-03-12 浙江工业大学 基于小波变换空间变迹的sar图像旁瓣抑制方法
CN108008381B (zh) * 2017-06-09 2019-01-01 北京航空航天大学 一种基于多方位角sar图像的目标方位旁瓣抑制方法及装置
CN108318891B (zh) * 2017-11-28 2021-09-10 西安电子科技大学 一种基于改进sva和cs的sal数据旁瓣的压低方法
CN108229404B (zh) * 2018-01-09 2022-03-08 东南大学 一种基于深度学习的雷达回波信号目标识别方法
CN108957406A (zh) * 2018-05-18 2018-12-07 电子科技大学 一种基于神经网络的雷达主瓣干扰抑制方法

Also Published As

Publication number Publication date
CN110333489A (zh) 2019-10-15

Similar Documents

Publication Publication Date Title
CN110333489B (zh) 采用cnn与rsva组合对sar回波数据旁瓣抑制的处理方法
EP3588128B1 (en) Method for detection and height and azimuth estimation of objects in a scene by radar processing using sparse reconstruction with coherent and incoherent arrays
CN110412568B (zh) 基于扩展方位角相位编码的距离模糊杂波抑制方法
CN110275166A (zh) 基于admm的快速稀疏孔径isar自聚焦与成像方法
CN106680815B (zh) 基于张量稀疏表示的mimo雷达成像方法
CN112099008A (zh) 基于cv-admmn的sa-isar成像与自聚焦方法
CN112612005B (zh) 基于深度学习的雷达抗主瓣干扰方法
CN112612006B (zh) 基于深度学习的机载雷达非均匀杂波抑制方法
CN113655478B (zh) 一种成像方法及装置
CN113589287B (zh) 合成孔径雷达稀疏成像方法、装置、电子设备及存储介质
CN112198506B (zh) 一种超宽带穿墙雷达学习成像的方法、装置、系统和可读存储介质
CN106845343B (zh) 一种光学遥感图像海上平台自动检测方法
CN104076360A (zh) 基于压缩感知的二维sar稀疏目标成像方法
CN111948652B (zh) 一种基于深度学习的sar智能参数化超分辨成像方法
CN117192548A (zh) 基于深度展开的稀疏isar高分辨成像方法
CN112147608A (zh) 一种快速高斯网格化非均匀fft穿墙成像雷达bp方法
CN110471036B (zh) 一种用于大阵列近场聚焦中的假目标清扫方法
CN115453528A (zh) 基于快速sbl算法实现分段观测isar高分辨成像方法及装置
CN111781598A (zh) 基于dsn的高分辨二维isar成像方法
Wei et al. CARNet: An effective method for SAR image interference suppression
CN115240040A (zh) 一种穿墙雷达人体行为特征增强方法及装置
CN113759372A (zh) 弹载大斜视小孔径多通道sar的成像方法
CN117471457A (zh) 基于深度展开复网络的稀疏sar学习成像方法、装置及介质
CN108931770B (zh) 基于多维贝塔过程线性回归的isar成像方法
CN112946644B (zh) 基于最小化卷积加权l1范数的稀疏孔径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