CN109001736A - 一种基于深度时空预测神经网络的雷达回波外推方法 - Google Patents
一种基于深度时空预测神经网络的雷达回波外推方法 Download PDFInfo
- Publication number
- CN109001736A CN109001736A CN201810599816.9A CN201810599816A CN109001736A CN 109001736 A CN109001736 A CN 109001736A CN 201810599816 A CN201810599816 A CN 201810599816A CN 109001736 A CN109001736 A CN 109001736A
- Authority
- CN
- China
- Prior art keywords
- data
- size
- network
- sequence
- input
- 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.)
- Granted
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/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Electromagnetism (AREA)
- Complex Calculations (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明公开了一种基于深度时空预测神经网络的雷达回波外推方法,包括:步骤1,数据预处理:对给定的多普勒天气雷达基数据,通过数据插值得到等高平面位置显示(Constant Altitude Plan Position Indicator,CAPPI)数据,将CAPPI数据转换为归一化灰度数据,基于归一化灰度数据集划分得到测试序列样本集;步骤2,雷达回波外推:将测试序列样本集输入深度时空预测神经网络,通过网络的前向传播输出预测序列。
Description
技术领域
本发明属于大气探测中地面气象观测技术领域,尤其涉及一种基于深度时空预测神经网络的雷达回波外推方法。
背景技术
对流天气临近预报是指短时内高时空分辨率的对流天气系统发生与演变过程的预报,准确、可靠且及时的预报对于防灾减灾、危险天气预警和保障国民生命财产安全等都具有重要意义。现有对流天气临近预报技术主要包括雷达回波外推预报技术、数值预报技术和概念模型预报技术等,其中雷达回波外推预报技术是对流天气临近预报的最基础手段和最有效途径。
雷达回波外推是根据已有历史雷达回波,预测未来回波的形状、位置、强度和运动信息等,以实现对对流天气系统的跟踪预报。传统的雷达回波外推方法主要包括质心追踪法、交叉相关法和光流法。质心追踪法通过计算连续的回波单体质心来确定移动矢量,进而预测未来回波位置。但质心追踪法使用范围局限于回波较强、范围较小的对流天气系统。交叉相关法将回波划分不同的追踪区域,计算相邻时刻区域之间的最优相关系数,进而确定区域间拟合关系以及实现预测。光流法通过计算连续回波的光流场得到运动矢量场,再基于运动矢量实现外推。然而,对流天气系统是动态复杂的系统,它具有作非线性运动(如旋转)、运动中有形变(如扩张和缩小)以及生消发展变化快等特征。传统雷达回波外推方法只假设回波简单线性演变,且对历史雷达回波资料利用率不足,无法有效预测回波的旋转变化和生消发展,具有精确度低、预报时效短等缺陷。
近期,深度学习技术在计算机视觉、语音识别和自然语言处理等领域得到广泛应用。深度学习本质就是通过构建合理的神经网络结构(如卷积神经网络CNN,循环神经网络RNN),使用一定的训练方法从大量数据中提取所需要的模型。在雷达回波外推问题中,具有大量的观测数据用以训练,长期积累的历史雷达回波数据也隐含了雷达回波变化规律,因此从深度学习的观点出发,可以通过训练合理的神经网络模型来解决这一难题。在众多神经网络模型中,长短期记忆网络(Long-Short Term Memory,LSTM)是RNN的一种变体,它通过引入门结构和记忆单元来解决长期依赖问题,是目前被广泛用于解决序列学习问题的热门网络模型之一。
基于LSTM网络结构,为提高雷达回波外推的准确度和预报时效,针对雷达回波的运动变化特点,设计深度时空预测神经网络,该网络由初始化网络、编码网络和预测网络组成,编码网络和预测网络则由空间形变-时序轨迹LSTM(Spatial Variation andTemporalTrajectory LSTM,SVTT-LSTM)基本构造单元时序连接并堆叠构成,网络能够有效记忆回波非线性运动、生消变化等特征。训练网络使其从历史回波数据中充分学习雷达回波运动变化规律,对于提高雷达回波外推精确度和预报时效具有重要意义。
发明内容
发明目的:本发明所要解决的技术问题是针对现有的雷达回波外推方法的准确度低、预报时效短,提出了一种基于深度时空预测神经网络的雷达回波外推方法,实现对雷达回波的外推预测,包括以下步骤:
包括以下步骤:
步骤1,数据预处理:对给定的多普勒天气雷达基数据,通过数据插值得到CAPPI数据,将CAPPI数据转换为归一化灰度数据,基于归一化灰度数据集划分得到测试序列样本集;
步骤2,雷达回波外推:将测试序列样本集输入深度时空预测神经网络,通过网络的前向传播输出预测序列。
步骤1包括以下步骤:
步骤1-1,数据插值:基于反距离加权法进行数据插值,通过坐标变换将多普勒天气雷达基数据的空间极坐标转化为空间直角坐标,并划分插值网格,再遍历所有空间直角坐标点,计算其对影响范围内所有插值网格点的权重和贡献值,最后遍历所有网格点,计算插值后网格点数据,得到CAPPI数据;
步骤1-2,数据转换:对步骤1-1数据插值得到的CAPPI数据,通过分辨率调整转化为分辨率240×240的数据,再通过数据映射和归一化将反射率数据转化为归一化灰度数据;
步骤1-3,样本集划分:对步骤1-2得到的归一化灰度数据集按时间顺序排列,将数据分段划分为序列并汇总得到序列样本总集,再划分得到测试序列样本集。
步骤1-1包括以下步骤:
步骤1-1-1,坐标变换:通过坐标变换将多普勒天气雷达基数据的空间极坐标数据转化为空间直角坐标数据;多普勒天气雷达基数据根据空间极坐标确定空间位置,表示为(r,θ,φ),其中r表示坐标点到原点距离,θ表示仰角,φ表示方位角,将其变换为空间直角坐标(x,y,z),其中x表示水平橫坐标,y表示水平纵坐标,z表示垂直高度,变换公式如下:
其中,R表示地球半径;
步骤1-1-2,插值网格划分:设定插值网格橫、纵坐标值域为[-480,480],网格分辨率为1,划分出分辨率为960×960的插值网格;
步骤1-1-3,权重和贡献值计算:遍历步骤1-1-1坐标变换后得到的每一空间直角坐标点,计算其对所有影响插值网格点的影响权重和贡献值,并存储于对应网格点的权重矩阵weight_matrix和贡献值矩阵contribution_matrix中,其中,空间直角坐标点影响的插值网格点集表示为:
其中,(x′,y′)表示插值网格点坐标,height表示CAPPI高度,affect_radius表示影响半径;
基于反距离加权进行数据插值,则空间直角坐标点对网格点的影响权重w随两者距离增加指数级衰减,贡献值c则等于权重和坐标点数据值乘积,计算公式如下:
w=d-2
c=w×reflectivity
其中,d表示空间直角坐标点和网格点间的距离,reflectivity表示空间直角坐标点的反射率数据值;
步骤1-1-4,插值数据计算:对步骤1-1-2插值网格划分得到的每一插值网格点,CAPPI数据cappi_data计算公式如下:
步骤1-2数据转换包括以下步骤:
步骤1-2-1,分辨率调整:对步骤1-1-4数据插值得到的CAPPI数据,通过裁剪保留中间区域分辨率为480×480部分数据,再通过数据压缩将裁剪后数据分辨率压缩至240×240,将分辨率调整后的反射率CAPPI数据记为adjustment_data;
步骤1-2-2,数据映射和归一化:对步骤1-2-1分辨率调整后得到的反射率数据adjustment_data,将其映射为值域[0,255]的灰度pixel数据,再通过归一化处理得到归一化灰度数据normalization_pixel,数据映射公式如下:
其中,表示向下取整;
数据归一化公式为:
最终得到的归一化灰度数据特征数为1,分辨率为240×240。
步骤1-3样本集划分包括以下步骤:
步骤1-3-1,序列划分:将所有经步骤1-2数据转换得到的归一化灰度数据normalization_pixel汇总得到归一化灰度数据集,再将其中数据分段划分为序列;首先将数据按时间顺序排列,随后将每20个数据划分为一个序列sequence,其中,前10个数据作为输入序列input,后10个数据作为对照序列contrast,且在划分时相邻序列之间有10个重叠数据,即前一序列的对照序列将作为后一序列的输入序列,序列表示为:
sequence={input,contrast},
其中,input={x1,x2,...,x10},contrast={x11,x12,...,x20},x1,x2,...,x10表示输入序列input中第1到第10个数据,x11,x12,...,x20表示对照序列contrast中第1到第10个数据;
序列总数count_of_sequence由如下公式确定:
其中,N表示归一化灰度数据集中数据总数;
步骤1-3-2,序列集划分:对步骤1-3-1序列划分得到的所有序列汇总为序列样本总集total_sequence_set,按照比例将序列样本总集中所有序列随机划分出测试序列样本集test_sequence_set。
步骤2包括以下步骤:
步骤2-1,测试序列样本输入:将步骤1-3得到的测试序列样本集中所有序列样本输入深度时空预测神经网络中进行测试,每个序列样本为sequence={input,contrast},共20个数据,其中input={x1,x2,...,x10}作为输入序列,contrast={x11,x12,...,x20}作为对照序列;
步骤2-2,前向传播:对步骤2-1输入的每个序列样本,初始化网络将输入序列前2个数据作为输入,前向传播输出编码网络所需的隐藏态和初始记忆;编码网络将输入序列的后续数据和初始化网络的输出作为输入,前向传播输出编码输入得到的隐藏态和记忆;预测网络将编码网络的输出作为输入,最终解码输出预测序列。
步骤2-2包括以下步骤:
步骤2-2-1,初始状态输出:对步骤2-1输入的每个序列样本,初始化网络将其输入序列10个数据中的前2个数据x1、x2作为输入,其中x1,x2∈R1×240×240,R表示实数域,1×240×240表示特征数为1,分辨率为240×240;
x1、x2输入初始化网络的第1层网络ini_layer1,经过前向传播后得到第1层网络输出第1层隐藏态第1层时序轨迹初始记忆和中间记忆S0、S1,前向传播公式如下:
其中,*表示卷积,tanh为激活函数双曲正切卷积核尺寸为5×5,数量为1×128;卷积核尺寸为5×5,数量为128×128;卷积核尺寸为5×5,数量为1×32;卷积核尺寸为5×5,数量为128×32;偏置特征图数量为128;偏置特征图数量为32;卷积核和偏置的值均由网络训练过程确定;输出S0,S1∈R32×240×240;
输入初始化网络的第2层网络ini_layer2,经过前向传播后得到第2层网络输出第2层隐藏态第2层时序轨迹初始记忆和中间记忆S2,前向传播公式如下:
其中,卷积核尺寸为5×5,数量为128×64;卷积核尺寸为5×5,数量为64×64;卷积核尺寸为5×5,数量为64×32;偏置 特征图数量为64;偏置特征图数量为32;卷积核和偏置的值均由网络训练过程确定;输出
S0、S1、S2输入初始化网络的第3层网络ini_layer3,经过前向传播后得到第3层网络输出第3层隐藏态第3层时序轨迹初始记忆中间记忆S3和空间形变初始记忆前向传播公式如下:
其中,卷积核尺寸为7×7,数量为64×32;卷积核 尺寸为7×7,数量为32×32;偏置 bhs-i特征图数量为32;卷积核和偏置的值均由网络训练过程确定;输出
经过前向传播,初始化网络完成隐藏态和记忆的初始化,最终向编码网络输出隐藏态时序轨迹初始记忆和空间形变初始记忆
步骤2-2-2,输入序列编码:对步骤2-1输入的每个序列样本,编码网络将其输入序列10个数据中的后8个数据x3、x4、...、x10,以及步骤2-2-1初始化网络输出的隐藏态和初始记忆共同作为输入,其中x3,...,x10∈R1×240×240;
当t=3时,x3、输入第1层编码-空间形变-时序轨迹长短期记忆网络(encode Spatial Variation and Temporal Trajectory Long Short Term Memory,e-SVTT-LSTM),输出随后输入第2层e-SVTT-LSTM,输出最后输入第3层e-SVTT-LSTM,输出 e-SVTT-LSTM的前向传播公式如下:
其中,l∈[1,3]表示网络层数,t∈[3,10]表示时刻;σ表示sigmoid函数 表示Hadamard矩阵点积;[]表示在特征维上串联叠加;为输入门,为遗忘门,为输出门;为空间形变记忆,为时序轨迹记忆,为隐藏态,为结构生成子网络的隐藏层,为结构生成子网络输出的流场矢量;对于当l=1时,有对于输入当l=1时,有否则有
对扭曲函数令流场则:
其中,为流场矢量的第n个特征,表示分辨率为240×240;流场M∈RC×240×240,其特征图数量C和第l层隐藏态相同;Mc,i,j表示M第c个特征图坐标为(i,j)处的值,Hc,m,n表示第c个特征图坐标为(m,n)处的值,vi,j、ui,j表示坐标为(i,j)处的值;
对卷积核当l=1时,尺寸为5×5,数量为1×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;对卷积核当l=1时,尺寸为5×5,数量为32×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;对卷积核尺寸为5×5,当l=1时,数量为129×32,当l=2时,数量为192×32,当l=3时,数量为96×32;卷积核尺寸为5×5,数量为32×N,N表示动态连接总数;对卷积核其中n∈[1,N],当l=1时,尺寸为1×1,数量为128×128,当l=2时,尺寸为1×1,数量为64×64,当l=3时,尺寸为1×1,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为128×128,当l=2时,尺寸为5×5,数量为64×64,当l=3时,尺寸为7×7,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为256×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;对偏置 当l=1时,特征图数量为128,当l=2时,特征图数量为64,当l=3时,特征图数量为32;偏置特征图数量为32;偏置特征图数量为N;卷积核和偏置的值均由网络训练过程确定;输出
当t=10时,x10、输入第1层e-SVTT-LSTM,输出随后输入第2层e-SVTT-LSTM,输出最后 输入第3层e-SVTT-LSTM,输出其中,
经过前向传播,编码网络完成对输入序列的编码,最终向预测网络输出隐藏态 时序轨迹记忆和空间形变记忆
步骤2-2-3,预测序列输出:预测网络将步骤2-2-2编码网络编码输入得到的状态和记忆作为输入;
当t=11时,输入第3层预测-空间形变-时序轨迹长短期记忆网络(forecast Spatial Variation and Temporal Trajectory Long Short Term Memory,f-SVTT-LSTM),输出随后输入第2层f-SVTT-LSTM,输出最后输入第1层f-SVTT-LSTM,输出 输入f-SVTT-LSTM的预测输出单元,输出预测序列第1帧f-SVTT-LSTM的前向传播公式如下:
其中,l∈[1,3]表示网络层数,t∈[11,20]表示时刻;表示反卷积;为预测输出单元输出的预测序列中一帧;对于当l=3且t≠11时,有当l=3且t=11时,有对于输入l∈[1,2],且
对卷积核当l=1时,尺寸为5×5,数量为64×128,当l=2时,尺寸为7×7,数量为32×64;对卷积核 当l=1时,尺寸为5×5,数量为64×128,当l=2时,尺寸为7×7,数量为32×64,当l=3且t≠11时,尺寸为5×5,数量为128×32,当l=3且t=11时,尺寸为5×5,数量为32×32;对卷积核尺寸为5×5,当l=1时,数量为192×32,当l=2时,数量为96×32,当l=3时,数量为32×32;卷积核尺寸为5×5,数量为32×N,N表示动态连接总数;对卷积核其中n∈[1,N],当l=1时,尺寸为1×1,数量为128×128,当l=2时,尺寸为1×1,数量为64×64,当l=3时,尺寸为1×1,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为128×128,当l=2时,尺寸为5×5,数量为64×64,当l=3时,尺寸为7×7,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为256×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;卷积核Whx-f尺寸为5×5,数量为128×1;对偏置 当l=1时,特征图数量为128,当l=2时,特征图数量为64,当l=3时,特征图数量为32;偏置特征图数量为32;偏置特征图数量为N;偏置bhx-f特征图数量为1;卷积核和偏置的值均由网络训练过程确定;输出
当t=20时,输入第3层f-SVTT-LSTM,输出随后 输入第2层f-SVTT-LSTM,输出最后输入第1层f-SVTT-LSTM,输出输入f-SVTT-LSTM的预测输出单元,输出预测序列第10帧
经过前向传播,预测网络完成对编码得到的隐藏态和记忆的解码,最终输出预测序列完成雷达回波外推。
有益效果:基于深度时空预测神经网络,在雷达回波数据集上进行训练,并利用训练好的网络进行雷达回波外推,有效提高了雷达回波外推的准确度和预报时效。
具体而言本发明与已有方法相比有以下优点:1、外推准确度高,与其他方法比较,本发明中深度时空预测神经网络能够有效记忆回波非线性运动、生消变化等特征,精确度更高;2、外推时效长,可以通过调整网络结构增加预测序列输出数据帧数来进一步提高外推时效。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述或其他方面的优点将会变得更加清楚。
图1为本发明流程图。
图2为深度时空预测神经网络模型图。
图3为e-SVTT-LSTM结构图。
图4为f-SVTT-LSTM结构图。
具体实施方式
本发明公开一种基于深度时空预测神经网络的雷达回波外推方法,如图1所示,包括以下步骤:
步骤1,数据预处理:对给定的多普勒天气雷达基数据,通过数据插值得到CAPPI数据,将CAPPI数据转换为归一化灰度数据,基于归一化灰度数据集划分得到测试序列样本集;
步骤2,雷达回波外推:将测试序列样本集输入深度时空预测神经网络,通过网络的前向传播输出预测序列。
步骤1包括以下步骤:
步骤1-1,数据插值:基于反距离加权法进行数据插值,通过坐标变换将多普勒天气雷达基数据的空间极坐标转化为空间直角坐标,并划分插值网格,再遍历所有空间直角坐标点,计算其对影响范围内所有插值网格点的权重和贡献值,最后遍历所有网格点,计算插值后网格点数据,得到CAPPI数据;
步骤1-2,数据转换:对步骤1-1数据插值得到的CAPPI数据,通过分辨率调整转化为分辨率240×240的数据,再通过数据映射和归一化将反射率数据转化为归一化灰度数据;
步骤1-3,样本集划分:对步骤1-2得到的归一化灰度数据集按时间顺序排列,将数据分段划分为序列并汇总得到序列样本总集,再划分得到测试序列样本集。
步骤1-1包括以下步骤:
步骤1-1-1,坐标变换:通过坐标变换将多普勒天气雷达基数据的空间极坐标数据转化为空间直角坐标数据;多普勒天气雷达基数据根据空间极坐标确定空间位置,表示为(r,θ,φ),其中r表示坐标点到原点距离,θ表示仰角,φ表示方位角,将其变换为空间直角坐标(x,y,z),其中x表示水平橫坐标,y表示水平纵坐标,z表示垂直高度,变换公式如下:
其中,R表示地球半径;
步骤1-1-2,插值网格划分:设定插值网格橫、纵坐标值域为[-480,480],网格分辨率为1,划分出分辨率为960×960的插值网格;
步骤1-1-3,权重和贡献值计算:遍历步骤1-1-1坐标变换后得到的每一空间直角坐标点,计算其对所有影响插值网格点的影响权重和贡献值,并存储于对应网格点的权重矩阵weight_matrix和贡献值矩阵contribution_matrix中,其中,空间直角坐标点影响的插值网格点集表示为:
其中,(x′,y′)表示插值网格点坐标,height表示CAPPI高度,affect_radius表示影响半径;
基于反距离加权进行数据插值,则空间直角坐标点对网格点的影响权重w随两者距离增加指数级衰减,贡献值c则等于权重和坐标点数据值乘积,计算公式如下:
w=d-2
c=w×reflectivity
其中,d表示空间直角坐标点和网格点间的距离,reflectivity表示空间直角坐标点的反射率数据值;
步骤1-1-4,插值数据计算:对步骤1-1-2插值网格划分得到的每一插值网格点,CAPPI数据cappi_data计算公式如下:
步骤1-2数据转换包括以下步骤:
步骤1-2-1,分辨率调整:对步骤1-1-4数据插值得到的CAPPI数据,通过裁剪保留中间区域分辨率为480×480部分数据,再通过数据压缩将裁剪后数据分辨率压缩至240×240,将分辨率调整后的反射率CAPPI数据记为adjustment_data;
步骤1-2-2,数据映射和归一化:对步骤1-2-1分辨率调整后得到的反射率数据adjustment_data,将其映射为值域[0,255]的灰度pixel数据,再通过归一化处理得到归一化灰度数据normalization_pixel,数据映射公式如下:
其中,表示向下取整;
数据归一化公式为:
最终得到的归一化灰度数据特征数为1,分辨率为240×240。
步骤1-3样本集划分包括以下步骤:
步骤1-3-1,序列划分:将所有经步骤1-2数据转换得到的归一化灰度数据normalization_pixel汇总得到归一化灰度数据集,再将其中数据分段划分为序列;首先将数据按时间顺序排列,随后将每20个数据划分为一个序列sequence,其中,前10个数据作为输入序列input,后10个数据作为对照序列contrast,且在划分时相邻序列之间有10个重叠数据,即前一序列的对照序列将作为后一序列的输入序列,序列表示为:
sequence={input,contrast},
其中,input={x1,x2,...,x10},contrast={x11,x12,...,x20},x1,x2,...,x10表示输入序列input中第1到第10个数据,x11,x12,...,x20表示对照序列contrast中第1到第10个数据;
序列总数count_of_sequence由如下公式确定:
其中,N表示归一化灰度数据集中数据总数;
步骤1-3-2,序列集划分:对步骤1-3-1序列划分得到的所有序列汇总为序列样本总集total_sequence_set,按照比例将序列样本总集中所有序列随机划分出测试序列样本集test_sequence_set。
步骤2包括以下步骤:
步骤2-1,测试序列样本输入:将步骤1-3得到的测试序列样本集中所有序列样本输入深度时空预测神经网络中进行测试,每个序列样本为sequence={input,contrast},共20个数据,其中input={x1,x2,...,x10}作为输入序列,contrast={x11,x12,...,x20}作为对照序列;
步骤2-2,前向传播:对步骤2-1输入的每个序列样本,初始化网络将输入序列前2个数据作为输入,前向传播输出编码网络所需的隐藏态和初始记忆;编码网络将输入序列的后续数据和初始化网络的输出作为输入,前向传播输出编码输入得到的隐藏态和记忆;预测网络将编码网络的输出作为输入,最终解码输出预测序列。
步骤2-2包括以下步骤:
步骤2-2-1,初始状态输出:对步骤2-1输入的每个序列样本,初始化网络将其输入序列10个数据中的前2个数据x1、x2作为输入,其中x1,x2∈R1×240×240,R表示实数域,1×240×240表示特征数为1,分辨率为240×240;
x1、x2输入初始化网络的第1层网络ini_layer1,经过前向传播后得到第1层网络输出第1层隐藏态第1层时序轨迹初始记忆和中间记忆S0、S1,前向传播公式如下:
其中,*表示卷积,tanh为激活函数双曲正切卷积核尺寸为5×5,数量为1×128;卷积核尺寸为5×5,数量为128×128;卷积核尺寸为5×5,数量为1×32;卷积核尺寸为5×5,数量为128×32;偏置特征图数量为128;偏置特征图数量为32;卷积核和偏置的值均由网络训练过程确定;输出S0,S1∈R32×240×240;
输入初始化网络的第2层网络ini_layer2,经过前向传播后得到第2层网络输出第2层隐藏态第2层时序轨迹初始记忆和中间记忆S2,前向传播公式如下:
其中,卷积核尺寸为5×5,数量为128×64;卷积核尺寸为5×5,数量为64×64;卷积核尺寸为5×5,数量为64×32;偏置 特征图数量为64;偏置特征图数量为32;卷积核和偏置的值均由网络训练过程确定;输出
S0、S1、S2输入初始化网络的第3层网络ini_layer3,经过前向传播后得到第3层网络输出第3层隐藏态第3层时序轨迹初始记忆中间记忆S3和空间形变初始记忆前向传播公式如下:
其中,卷积核尺寸为7×7,数量为64×32;卷积核 尺寸为7×7,数量为32×32;偏置 bhs-i特征图数量为32;卷积核和偏置的值均由网络训练过程确定;输出
经过前向传播,初始化网络完成隐藏态和记忆的初始化,最终向编码网络输出隐藏态时序轨迹初始记忆和空间形变初始记忆
步骤2-2-2,输入序列编码:对步骤2-1输入的每个序列样本,编码网络将其输入序列10个数据中的后8个数据x3、x4、...、x10,以及步骤2-2-1初始化网络输出的隐藏态和初始记忆共同作为输入,其中x3,...,x10∈R1×240×240;
当t=3时,x3、输入第1层编码-空间形变-时序轨迹长短期记忆网络(encode Spatial Variation and Temporal Trajectory Long Short Term Memory,e-SVTT-LSTM),如图3所示为e-SVTT-LSTM结构图,输出随后 输入第2层e-SVTT-LSTM,输出最后 输入第3层e-SVTT-LSTM,输出e-SVTT-LSTM的前向传播公式如下:
其中,l∈[1,3]表示网络层数,t∈[3,10]表示时刻;σ表示sigmoid函数 表示Hadamard矩阵点积;[]表示在特征维上串联叠加;为输入门,为遗忘门,为输出门;为空间形变记忆,为时序轨迹记忆,为隐藏态,为结构生成子网络的隐藏层,为结构生成子网络输出的流场矢量;对于当l=1时,有对于输入当l=1时,有否则有
对扭曲函数令流场则:
其中,为流场矢量的第n个特征,表示分辨率为240×240;流场M∈RC×240×240,其特征图数量C和第l层隐藏态相同;Mc,i,j表示M第c个特征图坐标为(i,j)处的值,Hc,m,n表示第c个特征图坐标为(m,n)处的值,vi,j、ui,j分别表示坐标为(i,j)处的值、坐标为(i,j)处的值;
对卷积核当l=1时,尺寸为5×5,数量为1×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;对卷积核当l=1时,尺寸为5×5,数量为32×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;对卷积核尺寸为5×5,当l=1时,数量为129×32,当l=2时,数量为192×32,当l=3时,数量为96×32;卷积核尺寸为5×5,数量为32×N,N表示动态连接总数;对卷积核其中n∈[1,N],当l=1时,尺寸为1×1,数量为128×128,当l=2时,尺寸为1×1,数量为64×64,当l=3时,尺寸为1×1,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为128×128,当l=2时,尺寸为5×5,数量为64×64,当l=3时,尺寸为7×7,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为256×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;对偏置 当l=1时,特征图数量为128,当l=2时,特征图数量为64,当l=3时,特征图数量为32;偏置特征图数量为32;偏置特征图数量为N;卷积核和偏置的值均由网络训练过程确定;输出
当t=10时,x10、输入第1层e-SVTT-LSTM,输出随后输入第2层e-SVTT-LSTM,输出最后 输入第3层e-SVTT-LSTM,输出其中,
经过前向传播,编码网络完成对输入序列的编码,最终向预测网络输出隐藏态 时序轨迹记忆和空间形变记忆
步骤2-2-3,预测序列输出:预测网络将步骤2-2-2编码网络编码输入得到的状态和记忆作为输入;
当t=11时,输入第3层预测-空间形变-时序轨迹长短期记忆网络(forecast Spatial Variation and Temporal Trajectory Long Short Term Memory,f-SVTT-LSTM),如图4所示为f-SVTT-LSTM结构图,输出随后 输入第2层f-SVTT-LSTM,输出最后输入第1层f-SVTT-LSTM,输出 输入f-SVTT-LSTM的预测输出单元,输出预测序列第1帧f-SVTT-LSTM的前向传播公式如下:
其中,l∈[1,3]表示网络层数,t∈[11,20]表示时刻;表示反卷积;为预测输出单元输出的预测序列中一帧;对于当l=3且t≠11时,有当l=3且t=11时,有对于输入l∈[1,2],且
对卷积核当l=1时,尺寸为5×5,数量为64×128,当l=2时,尺寸为7×7,数量为32×64;对卷积核 当l=1时,尺寸为5×5,数量为64×128,当l=2时,尺寸为7×7,数量为32×64,当l=3且t≠11时,尺寸为5×5,数量为128×32,当l=3且t=11时,尺寸为5×5,数量为32×32;对卷积核尺寸为5×5,当l=1时,数量为192×32,当l=2时,数量为96×32,当l=3时,数量为32×32;卷积核尺寸为5×5,数量为32×N,N表示动态连接总数;对卷积核其中n∈[1,N],当l=1时,尺寸为1×1,数量为128×128,当l=2时,尺寸为1×1,数量为64×64,当l=3时,尺寸为1×1,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为128×128,当l=2时,尺寸为5×5,数量为64×64,当l=3时,尺寸为7×7,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为256×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;卷积核Whx-f尺寸为5×5,数量为128×1;对偏置 当l=1时,特征图数量为128,当l=2时,特征图数量为64,当l=3时,特征图数量为32;偏置特征图数量为32;偏置特征图数量为N;偏置bhx-f特征图数量为1;卷积核和偏置的值均由网络训练过程确定;输出
当t=20时,输入第3层f-SVTT-LSTM,输出随后 输入第2层f-SVTT-LSTM,输出最后输入第1层f-SVTT-LSTM,输出输入f-SVTT-LSTM的预测输出单元,输出预测序列第10帧
经过前向传播,预测网络完成对编码得到的隐藏态和记忆的解码,最终输出预测序列完成雷达回波外推。
本发明训练深度时空预测神经网络的过程如下:
步骤A-1,模型构造:如图2所示,构造深度时空预测神经网络模型;首先构造用于向编码网络输出隐藏态与初始记忆的初始化网络initialize_network,再构造用于编码输入的编码网络encode_network,最后构造用于解码输出预测序列的预测网络forecast_network;
步骤A-2,训练超参数设置:将网络初始学习率λ0设为0.0001,学习率衰减系数α设为0.95,衰减速度β设为100,当前网络学习率λ根据训练步数train_step确定,公式如下:
其中,训练步数train_step初始设为1,每进行一次网络训练训练步数将加1。
训练时每次输入网络的序列样本数量batch_size设为10;每一轮迭代的最大批训练次数设为max_batch_step,初始批训练次数batch_step设为1;训练的最大迭代次数为max_iteration_step=100,初始迭代次数iteration_step=1;其中最大批训练次数由如下公式确定:
其中,count_of_sequence表示序列总数;
步骤A-3,训练序列样本输入:采用批训练的形式,从训练序列样本集中读取batch_size个序列样本输入网络中进行训练,每个序列样本为sequence={input,contrast},共20个数据,其中input={x1,x2,...,x10}作为输入序列,x1,x2,...,x10是输入序列中第1到第10个数据,contrast={x11,x12,...x20}作为对照序列,x11,x12,...,x20是输入序列中第1到第10个数据;
步骤A-4,前向传播:对步骤A-3输入的每个序列样本,初始网络将输入序列前2个数据作为输入,前向传播输出编码网络所需的隐藏态和初始记忆;编码网络将输入序列的后续数据和初始网络的输出作为输入,前向传播输出编码输入得到的隐藏态和记忆;预测网络将编码网络的输出作为输入,最终解码输出预测序列;
步骤A-5,反向传播:对步骤A-4前向传播得到的预测序列,计算它与步骤A-3输入的序列样本中对照序列的损失函数,根据损失函数计算所有网络参数的梯度,最后根据梯度和学习率更新所有网络参数;
步骤A-6,迭代训练:根据批训练次数和迭代次数控制迭代训练过程,分为以下三种情况:
若训练序列样本集中存在未用于训练的序列样本,即batch_step<max_batch_step,则返回步骤A-3继续读取batch_size个序列样本输入网络进行训练;
若训练序列样本集中所有序列样本均已用于训练,即batch_step=max_batch_step,且当前迭代次数小于最大迭代次数,即iteration_step<max_iteration_step,则设置batch_step=1,返回步骤A-3继续读取batch_size个序列样本输入网络进行训练;
若训练序列样本集中所有序列样本均已用于训练,即batch_step=max_batch_step,且当前迭代次数等于最大迭代次数,即iteration_step=max_iteration_step,则结束深度时空预测神经网络训练。
步骤A-1模型构造包括以下步骤:
步骤A-1-1,构造初始网络initialize_network:初始网络由3个网络层构成,第1层网络ini_layer1将输入序列的前两个数据x1、x2作为输入,输出第1层隐藏态H11、H21,第1层时序轨迹初始记忆和中间记忆S0、S1,网络构造公式如下:
其中,卷积核尺寸设为5×5,数量设为1×128,初始值设为random表示以均匀分布方式随机取值;卷积核 尺寸设为5×5,数量设为128×128,初始值设为卷积核尺寸设为5×5,数量设为1×32,初始值设为卷积核尺寸设为5×5,数量设为128×32,初始值设为偏置特征图数量设为128,初始值设为0;偏置特征图数量设为32,初始值设为0;
第2层网络ini_layer2将第1层隐藏态输出作为输入,输出第2层隐藏态第2层时序轨迹初始记忆和中间记忆S2,网络构造公式如下:
其中,卷积核尺寸设为5×5,数量设为128×64,初始值设为卷积核尺寸设为5×5,数量设为64×64,初始值设为卷积核尺寸设为5×5,数量设为64×32,初始值设为偏置特征图数量设为64,初始值设为0;偏置特征图数量设为32,初始值设为0;
第3层网络ini_layer3将第2层隐藏态和第1、2层中间记忆S0、S1、S2作为输入,输出第3层隐藏态第3层时序轨迹初始记忆中间记忆S3和空间形变初始记忆网络构造公式如下:
其中,卷积核尺寸设为7×7,数量设为64×32,初始值设为卷积核 尺寸设为7×7,数量设为32×32,初始值设为偏置 bhs-i特征图数量设为32,初始值设为0;
步骤A-1-2,构造编码网络encode_network:编码网络采用e-SVTT-LSTM作为基本构造单元,e-SVTT-LSTM在时序上串联,并堆叠3个网络层构成编码网络;构造编码网络,首先构造e-SVTT-LSTM的空间形变记忆单元,再构造结构生成子网络和e-SVTT-LSTM的时序轨迹记忆单元,最后构造e-SVTT-LSTM的输出单元;
步骤A-1-3,构造预测网络forecast_network:预测网络采用f-SVTT-LSTM作为基本构造单元,f-SVTT-LSTM在时序上串联,并堆叠3个网络层构成预测网络;构造预测网络,首先调整隐藏态和f-SVTT-LSTM的空间形变记忆在网络层间流动方向,并改变空间形变记忆单元中卷积操作为反卷积操作,再构造结构生成子网络和f-SVTT-LSTM的时序轨迹记忆单元,最后构造f-SVTT-LSTM的输出单元和预测输出单元。
步骤A-1-2包括以下步骤:
步骤A-1-2-1,构造空间形变记忆单元:t时刻第l层编码网络的空间形变记忆从编码网络的l-1层垂直传递到l层,在编码网络顶层时则传递到下一时刻的最底层;在传递过程中,空间形变记忆在t时刻第l层编码网络的输入门和遗忘门的控制下实现信息的更新;构造空间形变记忆单元,首先构造控制信息输入进记忆的输入门,再构造控制信息被记忆所遗忘的遗忘门,最后构造空间形变记忆,构造公式如下:
对卷积核当l=1时,尺寸设为5×5,数量设为1×128,初始值设为当l=2时,尺寸设为5×5,数量设为128×64,初始值设为当l=3时,尺寸设为7×7,数量设为64×32,初始值设为对卷积核当l=1时,尺寸设为5×5,数量设为32×128,初始值设为当l=2时,尺寸设为5×5,数量设为128×64,初始值设为当l=3时,尺寸设为7×7,数量为64×32,初始值设为对偏置 当l=1时,特征图数量设为128,初始值设为0,当l=2时,特征图数量设为64,初始值设为0,当l=3时,特征图数量设为32,初始值设为0;
步骤A-1-2-2,构造结构生成子网络:构造结构生成子网络,即构造包含1个隐藏层的卷积神经网络,网络的输入为当前时刻输入和前一时刻隐藏态输出为流场矢量和网络构造公式如下:
其中,C1为输入的特征图数量,C2为输入的特征图数量;则卷积核尺寸设为5×5,数量设为(C1+C2)×32,初始值设为卷积核尺寸设为5×5,数量设为32×N,初始值设为其中N表示动态连接总数;偏置特征图数量设为32,初始值设为0;偏置特征图数量设为N,初始值设为0;
步骤A-1-2-3,构造时序轨迹记忆单元:时序轨迹记忆在第l层编码网络中,从t-1时刻水平传递到当前时刻t;在传递过程中,时序轨迹记忆根据结构生成子网络动态调整连接拓扑,同时在输入门和遗忘门的控制下实现信息的更新;构造时序轨迹记忆单元,首先构造控制信息输入进记忆的输入门,再构造控制信息被记忆所遗忘的遗忘门,最后构造时序轨迹记忆,构造公式如下:
对卷积核当l=1时,尺寸设为5×5,数量设为1×128,初始值设为当l=2时,尺寸设为5×5,数量设为128×64,初始值设为当l=3时,尺寸设为7×7,数量设为64×32,初始值设为对卷积核其中n∈[1,N],当l=1时,尺寸设为1×1,数量设为128×128,初始值设为当l=2时,尺寸设为1×1,数量设为64×64,初始值设为当l=3时,尺寸设为1×1,数量设为32×32,初始值设为对卷积核当l=1时,尺寸设为5×5,数量设为128×128,初始值设为当l=2时,尺寸为5×5,数量设为64×64,初始值设为当l=3时,尺寸设为7×7,数量设为32×32,初始值设为对偏置当l=1时,特征图数量设为128,初始值设为0,当l=2时,特征图数量设为64,初始值设为0,当l=3时,特征图数量设为32,初始值设为0;
步骤A-1-2-4,构造e-SVTT-LSTM输出单元:e-SVTT-LSTM输出单元在输出门的控制下,融合空间形变记忆和时序轨迹记忆输出隐藏态首先构造控制信息输出的输出门,再构造隐藏态输出,构造公式如下:
其中,对卷积核当l=1时,尺寸设为5×5,数量设为1×128,初始值设为当l=2时,尺寸设为5×5,数量设为128×64,初始值设为当l=3时,尺寸设为7×7,数量设为64×32,初始值设为对卷积核当l=1时,尺寸设为1×1,数量设为128×128,初始值设为当l=2时,尺寸设为1×1,数量设为64×64,初始值设为当l=3时,尺寸设为1×1,数量设为32×32,初始值设为对卷积核当l=1时,尺寸设为5×5,数量设为128×128,初始值设为当l=2时,尺寸设为5×5,数量设为64×64,初始值设为当l=3时,尺寸设为7×7,数量设为32×32,初始值设为对卷积核当l=1时,尺寸设为5×5,数量设为256×128,初始值设为当l=2时,尺寸设为5×5,数量设为128×64,初始值设为当l=3时,尺寸设为7×7,数量设为64×32,初始值设为对偏置当l=1时,特征图数量设为128,初始值设为0,当l=2时,特征图数量设为64,初始值设为0,当l=3时,特征图数量设为32,初始值设为0。
步骤A-1-3包括以下步骤:
步骤A-1-3-1,构造空间形变记忆单元:t时刻第l层预测网络的空间形变记忆St l从预测网络的l+1层垂直传递到l层,在预测网络最底层时则传递到下一时刻的顶层;在传递过程中,空间形变记忆在t时刻第l层预测网络的输入门和遗忘门的控制下实现信息的更新;构造空间形变记忆单元,首先构造控制信息输入进记忆的输入门,再构造控制信息被记忆所遗忘的遗忘门,最后构造空间形变记忆,构造公式如下:
对卷积核当l=1时,尺寸设为5×5,数量设为64×128,初始值设为当l=2时,尺寸设为7×7,数量设为32×64,初始值设为对卷积核当l=1时,尺寸设为5×5,数量设为64×128,初始值设为当l=2时,尺寸设为7×7,数量设为32×64,初始值设为当l=3且t≠11时,尺寸设为5×5,数量设为128×32,初始值设为当l=3且t=11时,尺寸设为5×5,数量设为32×32,初始值设为对偏置当l=1时,特征图数量设为128,初始值设为0,当l=2时,特征图数量设为64,初始值设为0,当l=3时,特征图数量设为32,初始值设为0;
步骤A-1-3-2,构造结构生成子网络:构造结构生成子网络,即构造包含1个隐藏层的卷积神经网络,网络的输入为当前时刻输入和前一时刻隐藏态输出为流场矢量和网络构造公式如下:
其中,C1为输入的特征图数量,C2为输入的特征图数量,当l=3时,无输入则记C1为0;则卷积核尺寸设为5×5,数量设为(C1+C2)×32,初始值设为卷积核尺寸设为5×5,数量设为32×N,初始值设为其中N表示动态连接总数;偏置特征图数量设为32,初始值设为0;偏置特征图数量设为N,初始值设为0;
步骤A-1-3-3,构造时序轨迹记忆单元:时序轨迹记忆在第l层预测网络中,从t-1时刻水平传递到当前时刻t;在传递过程中,时序轨迹记忆根据结构生成子网络动态调整连接拓扑,同时在输入门和遗忘门的控制下实现信息的更新;构造时序轨迹记忆单元,首先构造控制信息输入进记忆的输入门,再构造控制信息被记忆所遗忘的遗忘门,最后构造时序轨迹记忆,构造公式如下:
其中,对卷积核当l=1时,尺寸设为5×5,数量设为64×128,初始值设为当l=2时,尺寸设为7×7,数量设为32×64,初始值设为对卷积核其中n∈[1,N],当l=1时,尺寸设为1×1,数量设为128×128,初始值设为当l=2时,尺寸设为1×1,数量设为64×64,初始值设为当l=3时,尺寸设为1×1,数量设为32×32,初始值设为对卷积核当l=1时,尺寸设为5×5,数量设为128×128,初始值设为当l=2时,尺寸设为5×5,数量设为64×64,初始值设为当l=3时,尺寸设为7×7,数量设为32×32,初始值设为对偏置当l=1时,特征图数量设为128,初始值设为0,当l=2时,特征图数量设为64,初始值设为0,当l=3时,特征图数量设为32,初始值设为0;
步骤A-1-3-4,构造f-SVTT-LSTM输出单元:f-SVTT-LSTM输出单元在输出门的控制下,融合空间形变记忆和时序轨迹记忆输出隐藏态构造输出单元,首先构造控制信息输出的输出门,再构造隐藏态输出,构造公式如下:
其中,对卷积核当l=1时,尺寸设为5×5,数量设为64×128,初始值设为当l=2时,尺寸设为7×7,数量设为32×64,初始值设为对卷积核当l=1时,尺寸设为1×1,数量设为128×128,初始值设为当l=2时,尺寸设为1×1,数量设为64×64,初始值设为当l=3时,尺寸设为1×1,数量设为32×32,初始值设为对卷积核当l=1时,尺寸设为5×5,数量设为128×128,初始值设为当l=2时,尺寸设为5×5,数量设为64×64,初始值设为当l=3时,尺寸设为7×7,数量设为32×32,初始值设为对卷积核当l=1时,尺寸设为5×5,数量设为256×128,初始值设为当l=2时,尺寸设为5×5,数量设为128×64,初始值设为当l=3时,尺寸设为7×7,数量设为64×32,初始值设为对偏置当l=1时,特征图数量设为128,初始值设为0,当l=2时,特征图数量设为64,初始值设为0,当l=3时,特征图数量设为32,初始值设为0;
步骤A-1-3-5,构造预测输出单元:预测输出单元将预测网络第1层输出的隐藏态转化为预测输出构造公式如下:
其中,卷积核Whx-f尺寸设为5×5,数量设为128×1,初始值设为偏置bhx-f特征图数量设为1,初始值设为0。
步骤A-4前向传播包括以下步骤:
步骤A-4-1,输出初始状态:对步骤A-3输入的batch_size个序列样本的每个序列样本,初始化网络将其输入序列10个数据中的前2个数据x1、x2作为输入,其中x1,x2∈R1 ×240×240,R表示实数域,1×240×240表示特征数为1,分辨率为240×240;
x1、x2输入初始化网络的第1层网络ini_layer1,经过前向传播后得到第1层网络输出第1层隐藏态第1层时序轨迹初始记忆和中间记忆S0、S1,前向传播公式如下:
其中,卷积核和偏置的初始值在步骤A-1-1构造时已确定,之后每迭代训练一次将更新一次;输出S0,S1∈R32×240×240;
输入初始化网络的第2层网络ini_layer2,经过前向传播后得到第2层网络输出第2层隐藏态第2层时序轨迹初始记忆和中间记忆S2,前向传播公式如下:
其中,卷积核和偏置的初始值在步骤A-1-1构造时已确定,之后每迭代训练一次将更新一次;输出S2∈R32×240×240;
S0、S1、S2输入初始化网络的第3层网络ini_layer3,经过前向传播后得到第3层网络输出第3层隐藏态第3层时序轨迹初始记忆中间记忆S3和空间形变初始记忆前向传播公式如下:
其中,卷积核和偏置的初始值在步骤A-1-1构造时已确定,之后每迭代训练一次将更新一次;输出
经过前向传播,初始化网络完成隐藏态和记忆的初始化,最终向编码网络输出隐藏态时序轨迹初始记忆和空间形变初始记忆
步骤A-4-2,输入序列编码:对步骤A-3输入的batch_size个序列样本的每个序列样本,编码网络将其输入序列10个数据中的后8个数据x3、x4、...、x10,以及步骤A-4-1初始化网络输出的隐藏态和初始记忆共同作为输入,其中x3,...,x10∈R1×240×240;
当t=3时,x3、输入第1层e-SVTT-LSTM,输出随后输入第2层e-SVTT-LSTM,输出最后 输入第3层e-SVTT-LSTM,输出e-SVTT-LSTM的前向传播公式如下:
其中,卷积核和偏置的初始值在步骤A-1-2构造时已确定,之后每迭代训练一次将更新一次;输出
以此类推,当t=10时,x10、输入第1层e-SVTT-LSTM,输出 随后输入第2层e-SVTT-LSTM,输出最后 输入第3层e-SVTT-LSTM,输出其中,
经过前向传播,编码网络完成对输入序列的编码,最终向预测网络输出隐藏态 时序轨迹记忆和空间形变记忆
步骤A-4-3,预测序列输出:预测网络将步骤A-4-2编码网络编码输入得到的状态和记忆作为输入;
当t=11时,输入第3层f-SVTT-LSTM,输出随后 输入第2层f-SVTT-LSTM,输出最后输入第1层f-SVTT-LSTM,输出 输入f-SVTT-LSTM的预测输出单元,输出预测序列第1帧f-SVTT-LSTM的前向传播公式如下:
其中,卷积核和偏置的初始值在步骤A-1-3构造时已确定,之后每迭代训练一次将更新一次;输出
以此类推,当t=20时,输入第3层f-SVTT-LSTM,输出随后输入第2层f-SVTT-LSTM,输出最后 输入第1层f-SVTT-LSTM,输出输入f-SVTT-LSTM的预测输出单元,输出预测序列第10帧
经过前向传播,预测网络完成对编码得到的隐藏态和记忆的解码,最终输出预测序列
步骤A-5反向传播包括以下步骤:
步骤A-5-1,计算损失函数:对步骤A-4前向传播得到的预测序列和步骤A-3输入的序列样本中的对照序列x11,x12,...,x20计算均方误差和平均绝对误差之和作为损失函数L:
其中,表示batch_size个序列样本的第n个序列样本中,时刻为t的数据的坐标为(i,j)的数据值,表示batch_size个序列样本的第n个序列样本前向传播得到的预测序列中,时刻为t的数据的坐标为(i,j)的数据值;
步骤A-5-2,计算网络参数梯度:通过求偏导计算步骤A-5-1得到的损失函数对每个网络参数的梯度
其中,步骤A-1模型构造中所有卷积核和偏置参数构成网络参数总体,θi(i∈[1,n])表示第i个网络参数,n表示参数总数;
步骤A-5-3,更新网络参数:用当前网络学习率λ乘以步骤A-5-2计算得到的网络参数梯度▽θiL得到网络参数修正项,将原网络参数和修正项相减即实现网络参数的更新,公式如下:
其中,θ′i表示更新后的网络参数。
本发明提供了一种基于深度时空预测神经网络的雷达回波外推方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (7)
1.一种基于深度时空预测神经网络的雷达回波外推方法,其特征在于,包括以下步骤:
步骤1,数据预处理:对给定的多普勒天气雷达基数据,通过数据插值得到CAPPI数据,将CAPPI数据转换为归一化灰度数据,基于归一化灰度数据集划分得到测试序列样本集;
步骤2,雷达回波外推:将测试序列样本集输入深度时空预测神经网络,通过网络的前向传播输出预测序列。
2.根据权利要求1所述的方法,其特征在于,步骤1包括以下步骤:
步骤1-1,数据插值:基于反距离加权法进行数据插值,通过坐标变换将多普勒天气雷达基数据的空间极坐标转化为空间直角坐标,并划分插值网格,再遍历所有空间直角坐标点,计算其对影响范围内所有插值网格点的权重和贡献值,最后遍历所有网格点,计算插值后网格点数据,得到CAPPI数据;
步骤1-2,数据转换:对步骤1-1数据插值得到的CAPPI数据,通过分辨率调整转化为分辨率240×240的数据,再通过数据映射和归一化将反射率数据转化为归一化灰度数据;
步骤1-3,样本集划分:对步骤1-2得到的归一化灰度数据集按时间顺序排列,将数据分段划分为序列并汇总得到序列样本总集,再划分得到测试序列样本集。
3.根据权利要求2所述的方法,其特征在于,步骤1-1包括以下步骤:
步骤1-1-1,坐标变换:通过坐标变换将多普勒天气雷达基数据的空间极坐标数据转化为空间直角坐标数据;多普勒天气雷达基数据根据空间极坐标确定空间位置,表示为(r,θ,φ),其中r表示坐标点到原点距离,θ表示仰角,φ表示方位角,将其变换为空间直角坐标(x,y,z),其中x表示水平橫坐标,y表示水平纵坐标,z表示垂直高度,变换公式如下:
其中,R表示地球半径;
步骤1-1-2,插值网格划分:设定插值网格橫、纵坐标值域为[-480,480],网格分辨率为1,划分出分辨率为960×960的插值网格;
步骤1-1-3,权重和贡献值计算:遍历步骤1-1-1坐标变换后得到的每一空间直角坐标点,计算其对所有影响插值网格点的影响权重和贡献值,并存储于对应网格点的权重矩阵weight_matrix和贡献值矩阵contribution_matrix中,其中,空间直角坐标点影响的插值网格点集表示为:
其中,(x′,y′)表示插值网格点坐标,height表示CAPPI高度,affect_radius表示影响半径;
基于反距离加权进行数据插值,则空间直角坐标点对网格点的影响权重w随两者距离增加指数级衰减,贡献值c则等于权重和坐标点数据值乘积,计算公式如下:
w=d-2
c=w×reflectivity
其中,d表示空间直角坐标点和网格点间的距离,reflectivity表示空间直角坐标点的反射率数据值;
步骤1-1-4,插值数据计算:对步骤1-1-2插值网格划分得到的每一插值网格点,CAPPI数据cappi_data计算公式如下:
4.根据权利要求3所述的方法,其特征在于,步骤1-2数据转换包括以下步骤:
步骤1-2-1,分辨率调整:对步骤1-1-4数据插值得到的CAPPI数据,通过裁剪保留中间区域分辨率为480×480部分数据,再通过数据压缩将裁剪后数据分辨率压缩至240×240,将分辨率调整后的反射率CAPPI数据记为adjustment_data;
步骤1-2-2,数据映射和归一化:对步骤1-2-1分辨率调整后得到的反射率数据adjustment_data,将其映射为值域[0,255]的灰度pixel数据,再通过归一化处理得到归一化灰度数据normalization_pixel,数据映射公式如下:
其中,表示向下取整;
数据归一化公式为:
最终得到的归一化灰度数据特征数为1,分辨率为240×240。
5.根据权利要求4所述的方法,其特征在于,步骤1-3样本集划分包括以下步骤:
步骤1-3-1,序列划分:将所有经步骤1-2数据转换得到的归一化灰度数据normalization_pixel汇总得到归一化灰度数据集,再将其中数据分段划分为序列;首先将数据按时间顺序排列,随后将每20个数据划分为一个序列sequence,其中,前10个数据作为输入序列input,后10个数据作为对照序列contrast,且在划分时相邻序列之间有10个重叠数据,即前一序列的对照序列将作为后一序列的输入序列,序列表示为:
sequence={input,contrast},
其中,input={x1,x2,...,x10},contrast={x11,x12,...,x20},x1,x2,...,x10表示输入序列input中第1到第10个数据,x11,x12,...,x20表示对照序列contrast中第1到第10个数据;
序列总数count_of_sequence由如下公式确定:
其中,N表示归一化灰度数据集中数据总数;
步骤1-3-2,序列集划分:对步骤1-3-1序列划分得到的所有序列汇总为序列样本总集total_sequence_set,按照比例将序列样本总集中所有序列随机划分出测试序列样本集test_sequence_set。
6.根据权利要求5所述的方法,其特征在于,步骤2包括以下步骤:
步骤2-1,测试序列样本输入:将步骤1-3得到的测试序列样本集中所有序列样本输入深度时空预测神经网络中进行测试,每个序列样本为sequence={input,contrast},共20个数据,其中input={x1,x2,...,x10}作为输入序列,contrast={x11,x12,...,x20}作为对照序列;
步骤2-2,前向传播:对步骤2-1输入的每个序列样本,初始化网络将输入序列前2个数据作为输入,前向传播输出编码网络所需的隐藏态和初始记忆;编码网络将输入序列的后续数据和初始化网络的输出作为输入,前向传播输出编码输入得到的隐藏态和记忆;预测网络将编码网络的输出作为输入,最终解码输出预测序列。
7.根据权利要求6所述的方法,其特征在于,步骤2-2包括以下步骤:
步骤2-2-1,初始状态输出:对步骤2-1输入的每个序列样本,初始化网络将其输入序列10个数据中的前2个数据x1、x2作为输入,其中x1,x2∈R1×240×240,R表示实数域,1×240×240表示特征数为1,分辨率为240×240;
x1、x2输入初始化网络的第1层网络ini_layer1,经过前向传播后得到第1层网络输出第1层隐藏态第1层时序轨迹初始记忆和中间记忆S0、S1,前向传播公式如下:
其中,*表示卷积,tanh为激活函数双曲正切卷积核尺寸为5×5,数量为1×128;卷积核尺寸为5×5,数量为128×128;卷积核尺寸为5×5,数量为1×32;卷积核尺寸为5×5,数量为128×32;偏置特征图数量为128;偏置特征图数量为32;卷积核和偏置的值均由网络训练过程确定;输出S0,S1∈R32×240×240;
输入初始化网络的第2层网络ini_layer2,经过前向传播后得到第2层网络输出第2层隐藏态第2层时序轨迹初始记忆和中间记忆S2,前向传播公式如下:
其中,卷积核尺寸为5×5,数量为128×64;卷积核尺寸为5×5,数量为64×64;卷积核尺寸为5×5,数量为64×32;偏置 特征图数量为64;偏置特征图数量为32;卷积核和偏置的值均由网络训练过程确定;输出S2∈R32×240×240;
S0、S1、S2输入初始化网络的第3层网络ini_layer3,经过前向传播后得到第3层网络输出第3层隐藏态第3层时序轨迹初始记忆中间记忆S3和空间形变初始记忆前向传播公式如下:
其中,卷积核尺寸为7×7,数量为64×32;卷积核 尺寸为7×7,数量为32×32;偏置 bhs-i特征图数量为32;卷积核和偏置的值均由网络训练过程确定;输出
经过前向传播,初始化网络完成隐藏态和记忆的初始化,最终向编码网络输出隐藏态时序轨迹初始记忆和空间形变初始记忆
步骤2-2-2,输入序列编码:对步骤2-1输入的每个序列样本,编码网络将其输入序列10个数据中的后8个数据x3、x4、...、x10,以及步骤2-2-1初始化网络输出的隐藏态和初始记忆共同作为输入,其中x3,...,x10∈R1×240×240;
当t=3时,x3、输入第1层编码-空间形变-时序轨迹长短期记忆网络e-SVTT-LSTM,输出随后输入第2层e-SVTT-LSTM,输出最后输入第3层e-SVTT-LSTM,输出 e-SVTT-LSTM的前向传播公式如下:
其中,l∈[1,3]表示网络层数,t∈[3,10]表示时刻;σ表示sigmoid函数 表示Hadamard矩阵点积;[]表示在特征维上串联叠加;为输入门,为遗忘门,为输出门;为空间形变记忆,为时序轨迹记忆,为隐藏态,为结构生成子网络的隐藏层,为结构生成子网络输出的流场矢量;对于当l=1时,有对于输入当l=1时,有否则有
对扭曲函数令流场则:
其中,为流场矢量的第n个特征,表示分辨率为240×240;流场M∈RC×240×240,其特征图数量C和第l层隐藏态相同;Mc,i,j表示M第c个特征图坐标为(i,j)处的值,Hc,m,n表示第c个特征图坐标为(m,n)处的值,vi,j、ui,j表示坐标为(i,j)处的值;
对卷积核当l=1时,尺寸为5×5,数量为1×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;对卷积核当l=1时,尺寸为5×5,数量为32×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;对卷积核尺寸为5×5,当l=1时,数量为129×32,当l=2时,数量为192×32,当l=3时,数量为96×32;卷积核尺寸为5×5,数量为32×N,N表示动态连接总数;对卷积核其中n∈[1,N],当l=1时,尺寸为1×1,数量为128×128,当l=2时,尺寸为1×1,数量为64×64,当l=3时,尺寸为1×1,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为128×128,当l=2时,尺寸为5×5,数量为64×64,当l=3时,尺寸为7×7,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为256×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;对偏置 当l=1时,特征图数量为128,当l=2时,特征图数量为64,当l=3时,特征图数量为32;偏置特征图数量为32;偏置特征图数量为N;卷积核和偏置的值均由网络训练过程确定;输出
当t=10时,x10、输入第1层e-SVTT-LSTM,输出随后输入第2层e-SVTT-LSTM,输出最后 输入第3层e-SVTT-LSTM,输出其中,
经过前向传播,编码网络完成对输入序列的编码,最终向预测网络输出隐藏态 时序轨迹记忆和空间形变记忆
步骤2-2-3,预测序列输出:预测网络将步骤2-2-2编码网络编码输入得到的状态和记忆作为输入;
当t=11时,输入第3层预测-空间形变-时序轨迹长短期记忆网络f-SVTT-LSTM,输出随后输入第2层f-SVTT-LSTM,输出最后输入第1层f-SVTT-LSTM,输出 输入f-SVTT-LSTM的预测输出单元,输出预测序列第1帧f-SVTT-LSTM的前向传播公式如下:
其中,l∈[1,3]表示网络层数,t∈[11,20]表示时刻;表示反卷积;为预测输出单元输出的预测序列中一帧;对于当l=3且t≠11时,有当l=3且t=11时,有对于输入l∈[1,2],且
对卷积核当l=1时,尺寸为5×5,数量为64×128,当l=2时,尺寸为7×7,数量为32×64;对卷积核 当l=1时,尺寸为5×5,数量为64×128,当l=2时,尺寸为7×7,数量为32×64,当l=3且t≠11时,尺寸为5×5,数量为128×32,当l=3且t=11时,尺寸为5×5,数量为32×32;对卷积核尺寸为5×5,当l=1时,数量为192×32,当l=2时,数量为96×32,当l=3时,数量为32×32;卷积核尺寸为5×5,数量为32×N,N表示动态连接总数;对卷积核其中n∈[1,N],当l=1时,尺寸为1×1,数量为128×128,当l=2时,尺寸为1×1,数量为64×64,当l=3时,尺寸为1×1,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为128×128,当l=2时,尺寸为5×5,数量为64×64,当l=3时,尺寸为7×7,数量为32×32;对卷积核当l=1时,尺寸为5×5,数量为256×128,当l=2时,尺寸为5×5,数量为128×64,当l=3时,尺寸为7×7,数量为64×32;卷积核Whx-f尺寸为5×5,数量为128×1;对偏置 当l=1时,特征图数量为128,当l=2时,特征图数量为64,当l=3时,特征图数量为32;偏置特征图数量为32;偏置特征图数量为N;偏置bhx-f特征图数量为1;卷积核和偏置的值均由网络训练过程确定;输出
当t=20时,输入第3层f-SVTT-LSTM,输出随后 输入第2层f-SVTT-LSTM,输出最后输入第1层f-SVTT-LSTM,输出 输入f-SVTT-LSTM的预测输出单元,输出预测序列第10帧
经过前向传播,预测网络完成对编码得到的隐藏态和记忆的解码,最终输出预测序列完成雷达回波外推。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810599816.9A CN109001736B (zh) | 2018-06-12 | 2018-06-12 | 一种基于深度时空预测神经网络的雷达回波外推方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810599816.9A CN109001736B (zh) | 2018-06-12 | 2018-06-12 | 一种基于深度时空预测神经网络的雷达回波外推方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109001736A true CN109001736A (zh) | 2018-12-14 |
CN109001736B CN109001736B (zh) | 2022-04-05 |
Family
ID=64601186
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810599816.9A Active CN109001736B (zh) | 2018-06-12 | 2018-06-12 | 一种基于深度时空预测神经网络的雷达回波外推方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109001736B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110222828A (zh) * | 2019-06-12 | 2019-09-10 | 西安交通大学 | 一种基于混合深度神经网络的非定常流场快速预测方法 |
CN110221360A (zh) * | 2019-07-25 | 2019-09-10 | 广东电网有限责任公司 | 一种电力线路雷暴预警方法及系统 |
CN110568442A (zh) * | 2019-10-15 | 2019-12-13 | 中国人民解放军国防科技大学 | 一种基于对抗外推神经网络的雷达回波外推方法 |
CN110674324A (zh) * | 2019-09-25 | 2020-01-10 | 上海眼控科技股份有限公司 | 雷达数据存储方法、装置、计算机设备和存储介质 |
CN111158068A (zh) * | 2019-12-31 | 2020-05-15 | 哈尔滨工业大学(深圳) | 一种基于简单卷积循环神经网络的短临预报方法及系统 |
CN111708030A (zh) * | 2020-05-28 | 2020-09-25 | 深圳市气象局(深圳市气象台) | 基于能量生成对抗性预测器的灾害天气预报方法 |
CN112446419A (zh) * | 2020-10-29 | 2021-03-05 | 中山大学 | 基于注意力机制的时空神经网络雷达回波外推预报方法 |
CN112698427A (zh) * | 2020-12-09 | 2021-04-23 | 最美天气(上海)科技有限公司 | 一种基于时空预测模型的短临预报方法及系统 |
CN112733072A (zh) * | 2020-02-27 | 2021-04-30 | 中国水利水电科学研究院 | 一种反距离平方加权空间插值方法 |
CN113640769A (zh) * | 2021-08-27 | 2021-11-12 | 南京信息工程大学 | 基于深度神经网络的天气雷达基本反射率预测方法 |
CN114488070A (zh) * | 2022-04-08 | 2022-05-13 | 北京弘象科技有限公司 | 一种基于深度学习模型的雷达回波外推方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5796611A (en) * | 1994-10-04 | 1998-08-18 | Nippon Telegraph And Telephone Corporation | Weather forecast apparatus and method based on recognition of echo patterns of radar images |
CN106886023A (zh) * | 2017-02-27 | 2017-06-23 | 中国人民解放军理工大学 | 一种基于动态卷积神经网络的雷达回波外推方法 |
CN107121679A (zh) * | 2017-06-08 | 2017-09-01 | 湖南师范大学 | 用于雷达回波外推的循环神经网络预测法及记忆单元结构 |
CN107247259A (zh) * | 2017-07-10 | 2017-10-13 | 西安电子科技大学 | 基于神经网络的k分布海杂波形状参数估计方法 |
-
2018
- 2018-06-12 CN CN201810599816.9A patent/CN109001736B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5796611A (en) * | 1994-10-04 | 1998-08-18 | Nippon Telegraph And Telephone Corporation | Weather forecast apparatus and method based on recognition of echo patterns of radar images |
CN106886023A (zh) * | 2017-02-27 | 2017-06-23 | 中国人民解放军理工大学 | 一种基于动态卷积神经网络的雷达回波外推方法 |
CN107121679A (zh) * | 2017-06-08 | 2017-09-01 | 湖南师范大学 | 用于雷达回波外推的循环神经网络预测法及记忆单元结构 |
CN107247259A (zh) * | 2017-07-10 | 2017-10-13 | 西安电子科技大学 | 基于神经网络的k分布海杂波形状参数估计方法 |
Non-Patent Citations (1)
Title |
---|
施恩等: "基于卷积神经网络的雷达回波外推方法", 《计算机应用》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110222828B (zh) * | 2019-06-12 | 2021-01-15 | 西安交通大学 | 一种基于混合深度神经网络的非定常流场预测方法 |
CN110222828A (zh) * | 2019-06-12 | 2019-09-10 | 西安交通大学 | 一种基于混合深度神经网络的非定常流场快速预测方法 |
CN110221360A (zh) * | 2019-07-25 | 2019-09-10 | 广东电网有限责任公司 | 一种电力线路雷暴预警方法及系统 |
CN110674324A (zh) * | 2019-09-25 | 2020-01-10 | 上海眼控科技股份有限公司 | 雷达数据存储方法、装置、计算机设备和存储介质 |
CN110568442B (zh) * | 2019-10-15 | 2021-08-20 | 中国人民解放军国防科技大学 | 一种基于对抗外推神经网络的雷达回波外推方法 |
CN110568442A (zh) * | 2019-10-15 | 2019-12-13 | 中国人民解放军国防科技大学 | 一种基于对抗外推神经网络的雷达回波外推方法 |
CN111158068A (zh) * | 2019-12-31 | 2020-05-15 | 哈尔滨工业大学(深圳) | 一种基于简单卷积循环神经网络的短临预报方法及系统 |
CN111158068B (zh) * | 2019-12-31 | 2022-09-23 | 哈尔滨工业大学(深圳) | 一种基于简单卷积循环神经网络的短临预报方法及系统 |
CN112733072B (zh) * | 2020-02-27 | 2021-11-16 | 中国水利水电科学研究院 | 一种反距离平方加权空间插值方法 |
CN112733072A (zh) * | 2020-02-27 | 2021-04-30 | 中国水利水电科学研究院 | 一种反距离平方加权空间插值方法 |
CN111708030A (zh) * | 2020-05-28 | 2020-09-25 | 深圳市气象局(深圳市气象台) | 基于能量生成对抗性预测器的灾害天气预报方法 |
CN111708030B (zh) * | 2020-05-28 | 2024-05-14 | 深圳市气象局(深圳市气象台) | 基于能量生成对抗性预测器的灾害天气预报方法 |
CN112446419A (zh) * | 2020-10-29 | 2021-03-05 | 中山大学 | 基于注意力机制的时空神经网络雷达回波外推预报方法 |
CN112446419B (zh) * | 2020-10-29 | 2023-07-11 | 中山大学 | 基于注意力机制的时空神经网络雷达回波外推预报方法 |
CN112698427A (zh) * | 2020-12-09 | 2021-04-23 | 最美天气(上海)科技有限公司 | 一种基于时空预测模型的短临预报方法及系统 |
CN113640769A (zh) * | 2021-08-27 | 2021-11-12 | 南京信息工程大学 | 基于深度神经网络的天气雷达基本反射率预测方法 |
CN114488070A (zh) * | 2022-04-08 | 2022-05-13 | 北京弘象科技有限公司 | 一种基于深度学习模型的雷达回波外推方法和装置 |
CN114488070B (zh) * | 2022-04-08 | 2022-07-19 | 北京弘象科技有限公司 | 一种基于深度学习模型的雷达回波外推方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109001736B (zh) | 2022-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108427989B (zh) | 一种用于雷达回波外推的深度时空预测神经网络训练方法 | |
CN109001736A (zh) | 一种基于深度时空预测神经网络的雷达回波外推方法 | |
CN110568442B (zh) | 一种基于对抗外推神经网络的雷达回波外推方法 | |
CN109410575B (zh) | 一种基于胶囊网络和嵌套式长短时记忆神经网络的路网状态预测方法 | |
CN110223323A (zh) | 基于深度特征自适应相关滤波的目标跟踪方法 | |
Li et al. | Ship motion prediction using dynamic seasonal RvSVR with phase space reconstruction and the chaos adaptive efficient FOA | |
CN105447235B (zh) | 一种土地未来利用情景动态模拟方法 | |
CN107703564A (zh) | 一种降雨预测方法、系统及电子设备 | |
CN109902259B (zh) | 一种轻量级的缺失时空数据的重构方法 | |
CN111681178B (zh) | 一种基于知识蒸馏的图像去雾方法 | |
CN108509910A (zh) | 基于fmcw雷达信号的深度学习手势识别方法 | |
CN110223234A (zh) | 基于级联收缩扩张的深度残差网络图像超分辨率重建方法 | |
CN108921298A (zh) | 强化学习多智能体沟通与决策方法 | |
CN112415521A (zh) | 基于cgru的强时空特性雷达回波临近预报方法 | |
CN108447057A (zh) | 基于显著性和深度卷积网络的sar图像变化检测方法 | |
CN111708030A (zh) | 基于能量生成对抗性预测器的灾害天气预报方法 | |
CN104915982A (zh) | 树冠层光照分布预测模型构建方法及光照分布检测方法 | |
CN106372597A (zh) | 基于自适应上下文信息的cnn交通检测方法 | |
CN115359338A (zh) | 一种基于混合学习模型的海表温度预测方法及系统 | |
CN116148796A (zh) | 一种基于雷达图像外推的强对流天气临近预报方法 | |
CN112180369A (zh) | 基于深度学习的一维综合孔径辐射计海面风速反演方法 | |
CN117079132A (zh) | 基于高斯距离损失的遥感图像目标检测方法 | |
CN116976227B (zh) | 一种基于lstm机器学习的风暴增水预报方法及系统 | |
CN117422619A (zh) | 图像重建模型的训练方法、图像重建方法、装置及设备 | |
CN116699731A (zh) | 一种热带气旋路径短期预报方法、系统及存储介质 |
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 |