CN114720956A - 针对双偏振天气雷达的水凝物粒子分类识别方法及装置 - Google Patents
针对双偏振天气雷达的水凝物粒子分类识别方法及装置 Download PDFInfo
- Publication number
- CN114720956A CN114720956A CN202210635178.8A CN202210635178A CN114720956A CN 114720956 A CN114720956 A CN 114720956A CN 202210635178 A CN202210635178 A CN 202210635178A CN 114720956 A CN114720956 A CN 114720956A
- Authority
- CN
- China
- Prior art keywords
- layer
- data
- dual
- polarization
- recognition model
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000002245 particle Substances 0.000 title claims abstract description 39
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title description 5
- 239000000017 hydrogel Substances 0.000 claims abstract description 19
- 230000010287 polarization Effects 0.000 claims abstract description 17
- 238000002310 reflectometry Methods 0.000 claims abstract description 13
- 230000010363 phase shift Effects 0.000 claims abstract description 8
- 210000002569 neuron Anatomy 0.000 claims description 49
- 230000006870 function Effects 0.000 claims description 45
- 239000011159 matrix material Substances 0.000 claims description 34
- 230000004913 activation Effects 0.000 claims description 30
- 238000012549 training Methods 0.000 claims description 15
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 10
- 238000007781 pre-processing Methods 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 7
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 238000013480 data collection Methods 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 4
- 230000009977 dual effect Effects 0.000 claims description 3
- 238000001556 precipitation Methods 0.000 abstract description 2
- 239000013078 crystal Substances 0.000 abstract 3
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 abstract 1
- 238000013135 deep learning Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000009833 condensation Methods 0.000 description 3
- 230000005494 condensation Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
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/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
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/048—Activation functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Molecular Biology (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Computer Networks & Wireless Communication (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
技术领域
本发明涉及天气雷达信息处理技术领域,具体的说,涉及一种针对双偏振天气雷达的水凝物粒子分类识别方法及装置。
背景技术
在中国中西部地区,多山的地形条件常常导致冰雹等破坏性较强的局地性气象灾害频发,极大地影响了粮食和经济作物的种植,对地方经济造成了重大破坏。以冰雹为代表的气象灾害产生的主要原因是对流性强天气系统的出现,为了减少气象灾害造成的经济损失,必须采取人工影响天气作业的技术手段,来抑制这些对流天气系统的形成和发展。而抑制对流性天气系统形成发展的关键是对这些系统的早期识别,如果能够在形成之前或者发展早期识别出这些系统,人工影响天气业务部门就可以提前进行必要的人工影响天气作业,以达到消除或者减弱这种灾害性天气系统的目的,进而阻止气象灾害的发生。
在中西部主要经济作物种植区域,为了气象灾害的防灾减灾,广泛开展了以天气雷达为指挥核心,以高炮或者火箭作为作业火器工具的人工影响天气作业业务。在人影作业业务中,作业指挥人员借助天气雷达的观测数据搜索可能形成灾害的天气系统目标,对可能形成对流性灾害天气系统目标核心进行精准识别和精确定位,并通过通信系统通知人影作业工具的现场操作人员,而作业操作人员依据指挥人员提供的作业参数调整高炮/火箭射角,准确及时地完成人影作业。
形成有效的人工影响天气作业的关键是基于天气雷达的灾害性天气系统的早期识别。传统的用于人影作业指挥天气雷达以多普勒天气雷达为主,它提供雷达观测范围内天气系统的反射率因子Z、速度V和谱宽W 3种观测数据,这3种观测数据仅仅能够反映天气系统的宏观结构,以此为依据识别是否是对流性天气系统,或者判别对流性天气系统的发展趋势,具有一定片面性,特别是无法进行对流性天气系统的早期识别,往往会造成对作业条件识别和作业时机选择的误判,影响人影作业的有效性。
近年来,随着天气雷达技术的发展,X波段双偏振多普勒天气雷达被广泛引入人工
影响天气作业指挥业务,极大地改善了灾害性天气系统的早期识别的准确性。双偏振多普
勒天气雷达同时或者交替发射水平/垂直极化电磁波,天气系统中不同相态、不同尺度和不
同形状的水凝物粒子对不同极化方向的电磁波的反射强度和相位偏移是不同的,这些差异
携带了粒子相态的准确识别信息,因此,双偏振多普勒天气雷达不仅仅能够提供天气系统
的反射率因子、速度V和谱宽W 3种观测数据,还可以提供水凝物粒子差分反射率
,差分相移率 以及交叉相关系数 等微观的双偏振观测参数,这些参数反映了
水凝物粒子的形状信息,可以用于对水凝物粒子的相态识别,以此为基础,结合反射率因子
Z、速度V和谱宽W 3种传统观参数,可以进行对流性天气系统的早期识别,提高人影作业的
效率和科学性。
传统人影作业指挥业务中,基于双偏振观测参数的水凝物粒子识别方法主要是基于模糊逻辑的方法,这种方法需要根据历史积累的观测数据,进行人工的隶属函数边界的设定选择,一方面参数选择人为因素较大,另一方面一旦识别模型建立,很难根据新的业务观测数据的积累,自动调整相应的识别参数,算法自适应性较差,不能很好改进。改进的基于模糊神经网络的识别方法,虽然能够通过样本训练确定隶属函数边界,部分解决了上述问题,但业务使用情况并不理想。
发明内容
为了解决上述技术问题,本发明提供了一种针对双偏振天气雷达的水凝物粒子分类识别方法及装置,引入深度机器学习的方法,用历史积累的观测数据建立训练样本,对识别模型进行的训练,建立稳定的识别模型,应用于人影作业指挥业务。在业务观测中,新的观测数据可用于分类识别模型的优化改进训练,以达到自适应改进的目的。
本发明的具体技术方案如下:
根据本发明的第一技术方案,提供一种针对双偏振天气雷达的水凝物粒子分类识别方法,所述方法包括:
将预处理后的双偏振观测参量数据输入分类识别模型以训练所述分类识别模型;
所述分类识别模型中的各个神经元都有对应的线性关系系数 和偏移,单个神
经元到单个神经元的线性关系为,其中表示为n*1的权重矩阵,
表示1*n的向量矩阵,将整个分类识别模型做一般化,假设第层共有 个神经元,则对
于第层的第个神经元的输出 ,有如(1)所示的前向传播:
所述Relu激活函数如(2)所示:
其中,x表示输入层矩阵的横标志,max()表示最大值函数;
所述Softmax激活函数如(3)所示:
确定所述分类识别模型的多类交叉熵损失函数如(4)所示:
(5)和(6)中,上标L表示第L层神经元,上标L-1表示第L-1层神经元, 表示选
用的激活函数, 表示L层n*1的矩阵向量,表示输出L-1层n*1的矩阵向量, 表示
输出层, 表第L层输出层, 表示输入层矩阵的竖标志, 标识Hadamard积,表示对
应位置的对应数值相乘;
利用所述权值文件配置所述分类识别模型,并将新数据流输入至所述分类识别模型中依次得到识别效果。
根据本发明各个实施例公开的一种针对双偏振天气雷达的水凝物粒子分类识别方法及装置,首先对偏振参量数据进行归一化,获得归一化后的偏振参量数据,然后建立分类识别模型,针对雷达数据结果选用合适的网络结构构建分类识别模型。在分类识别模型的基础上引入优化器和学习率对分类识别模型的识别结果进行优化。最后对于识别结果进行分类,获得水凝物粒子的类型。对实测数据的仿真实验结果表明,本方法下的偏振参量数据的降水粒子识别分类效果良好。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。在所有附图中,类似的元件或部分一般由类似的附图标记标识。附图中,各元件或部分并不一定按照实际的比例绘制。
图1示出了根据本发明实施例的一种针对双偏振天气雷达的水凝物粒子分类识别方法的流程图。
图2示出了根据本发明实施例的一种针对双偏振天气雷达的水凝物粒子分类识别方法所构建的分类识别模型结构图。
图3示出了根据本发明实施例的 K-fold交叉验证方法的示意图。
图4示出了传统模糊逻辑方法获得的水凝物粒子分类结果图。
图5示出了根据本发明实施例的一种针对双偏振天气雷达的水凝物粒子分类识别方法获得的水凝物粒子分类结果图。
图6示出了根据本发明实施例的一种针对双偏振天气雷达的水凝物粒子分类识别方法经过修正后获得的水凝物粒子分类结果图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步的详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定发明。
现在结合说明书附图对本发明做进一步的说明。
本发明实施例提供一种针对双偏振天气雷达的水凝物粒子分类识别方法,如图1和图2所示,所述方法包括:
对双偏振观测参量数据进行预处理,所述双偏振观测参量数据至少包括反射率因
子 、水凝物粒子差分反射率、差分相移率以及交叉相关系数中的一种。双偏
振观测参量数据可以基于双偏振天气雷达监测得到。在具体实施时,可以将输入的雷达基
数据按解码格式进行解码,利用结构体的方式分类存放、 、、四种双偏振观
测参量,后续将四种观测参量以函数参数表的形式输出,再以VCP层数、径向数、库数依次在
函数中存储用于后续调用。在调用四种观测参量之前,进行相应的预处理。双偏振观测参量
数据主要作为后续分类识别模型的训练数据集,可以基于以往的现有数据,也可以基于实
测得到的相应数据。
将预处理后的双偏振观测参量数据输入分类识别模型以训练所述分类识别模型。仅作为示例,所述分类识别模型可以选择DNN网络。DNN网络包括前向传播和后向传播。
所述分类识别模型中的各个神经元都有对应的线性关系系数 和偏移 ,单个
神经元到单个神经元的线性关系为 ,其中 表示为n*1的权重矩阵, 表示1*n的向量矩阵,将整个分类识别模型做一般化,假设第 层共有个神经
元,则对于第 层的第 个神经元的输出 ,有如(1)所示的前向传播:
所述Relu激活函数如(2)所示:
其中,x表示输入层矩阵的横标志,max()表示最大值函数;
所述Softmax激活函数如(3)所示:
在完成前向传播后,进行反向传播前,选择一个损失函数,用来度量训练样本计算出的输出和真实的训练样本输出之间的损失。本发明实施例确定所述分类识别模型的多类交叉熵损失函数如(4)所示:
(5)和(6)中,上标L表示第L层神经元,上标L-1表示第L-1层神经元, 表示选
用的激活函数,表示L层n*1的矩阵向量, 表示输出L-1层n*1的矩阵向量, 表示
输出层, 表示第L层输出层, 表示输入层矩阵的竖标志, 标识Hadamard积,表示
对应位置的对应数值相乘。
最后利用所述权值文件配置所述分类识别模型,并将新数据流输入至所述分类识
别模型中依次得到识别效果。新数据流具体是可以实时通过雷达检测的相关数据,包括反
射率因子 、水凝物粒子差分反射率 、差分相移率 以及交叉相关系数 中的
一种及其组合。本发明实施例最终获得一个如图5所示的水凝物粒子分类图。
在一些实施例中,在对双偏振观测参量数据进行预处理之前,所述方法还包括:将各个双偏振观测参量通过线性化的方法转换到预设的阈值范围内,实现对原始数据的等比例缩放,得到所述双偏振观测参量数据。仅作为示例,预设的阈值范围可以是[0 ,1]的范围。通过利用变量取值的最大值和最小值将原始数据转换为界于某一特定范围的数据,从而消除量纲和数量级影响,改变变量在分析中的权重。选用归一化而不是极值化方法的原因是由于极值化方法在对变量无量纲化过程中仅仅与该变量的最大值和最小值这两个极端值有关,而与其他取值无关,这使得该方法在改变各变量权重时过分依赖两个极端取值。
在一些实施例中,所述对双偏振观测参量数据进行预处理,包括:
通过如下(7)和(8)平滑所述双偏振观测参量数据:
在一些实施例中,在训练的过程中,通常会使用到优化器。发明人考虑到在训练的
初期Adam 中的方差会非常大。而起到修正更新方向的作用,因此 Adam 参数的更新量
的方差也会很大。选用RAdam进行学习率的纠正使前期方差较小选用的线性插值权重为,其中,Var(v∞)为无穷大时的方差,Var(v t )为t时刻的方差,基于线性插值权
重,通过如下(9)进行线性插值:
其中,为参数更新量,为学习率,Adam(t)为t时刻的Adam优化器,SGD_
Momentum(t)为t时刻的SGD + Momentum优化器,m t 为一阶矩(动量), 为迭代阈值,v t 为二
阶矩(自适应学习率), 均为指数滑动平均的超参。在一些实施例中,如图3所示,
在所述分类识别模型的训练过程中,对待输入的数据采用K-Fold 多次划分,以利用全部数
据集。
在一些实施例中,在所述分类识别模型的训练完成后,对所述分类识别模型进行校验,若学习率低于预设的阈值则对预处理后的双偏振观测参量数据进行正则化,直到学习率达到所述预设的阈值。若学习率过低,很可能会进入局部极值点就收敛,没有真正地找到最优解,导致水凝物粒子识别分类的判断结果不准确,所以要对数据进行正则化,防止学习率过低而造成的问题。
图4显示了传统模糊逻辑方法获得的水凝物粒子分类结果图,图5和图6显示了根据本发明实施例的一种针对双偏振天气雷达的水凝物粒子分类识别方法获得的水凝物粒子分类结果图,结合图4-图6,可以看出,本发明在实施过程中,可以不断自适应地调整,能够达到更准确的分类结果。
本发明实施例还提供一种针对双偏振天气雷达的水凝物粒子分类识别装置,其特征在于,所述装置包括处理器,所述处理器被配置为:
将预处理后的双偏振观测参量数据输入分类识别模型以训练所述分类识别模型;
所述分类识别模型中的各个神经元都有对应的线性关系系数 和偏移 ,单个
神经元到单个神经元的线性关系为 ,其中 表示为n*1的权重矩阵, 表示1*n的向量矩阵,将整个分类识别模型做一般化,假设第 层共有 个神经
元,则对于第 层的第 个神经元的输出 ,有如(1)所示的前向传播:
所述Relu激活函数如(2)所示:
其中,x表示输入层矩阵的横标志,max()表示最大值函数;
所述Softmax激活函数如(3)所示:
确定所述分类识别模型的多类交叉熵损失函数如(4)所示:
(5)和(6)中,上标L表示第L层神经元,上标L-1表示第L-1层神经元, 表示
选用的激活函数, 表示L层n*1的矩阵向量, 表示输出L-1层n*1的矩阵向量, 表
示输出层, 表示第L层输出层, 表示输入层矩阵的竖标志, 标识Hadamard积,表
示对应位置的对应数值相乘;
利用所述权值文件配置所述分类识别模型,并将新数据流输入至所述分类识别模型中依次得到识别效果。
需要说明的是,处理器可以是包括一个以上通用处理设备的处理设备,诸如微处理器、中央处理单元(CPU)、图形处理单元(GPU)等。更具体地,处理器可以是复杂指令集计算(CISC)微处理器、精简指令集计算(RISC)微处理器、超长指令字(VLIW)微处理器、运行其他指令集的处理器或运行指令集的组合的处理器。处理器还可以是一个以上专用处理设备,诸如专用集成电路(ASIC)、现场可编程门阵列(FPGA)、数字信号处理器(DSP)、片上系统(SoC)等。
处理器可以通信地耦合到存储器并且被配置为执行存储在其上的计算机可执行指令,以执行根据本发明各个实施例的针对双偏振天气雷达的水凝物粒子分类识别方法。
在一些实施例中,所述处理器被进一步配置为:
通过如下(7)和(8)平滑所述双偏振观测参量数据:
在一些实施例中,所述处理器被进一步配置为:
其中,为参数更新量,为学习率,Adam(t)为t时刻的Adam优化器,SGD_
Momentum(t)为t时刻的SGD + Momentum优化器,m t 为一阶矩(动量),为迭代阈值,v t 为二
阶矩(自适应学习率),均为指数滑动平均的超参。在一些实施例中,所述处理器被
进一步配置为:在所述分类识别模型的训练完成后,对所述分类识别模型进行校验,若学习
率低于预设的阈值则对预处理后的双偏振观测参量数据进行正则化,直到学习率达到所述
预设的阈值。
本发明实施例所提供的针对双偏振天气雷达的水凝物粒子分类识别装置与在前阐述的方法的技术效果基本一致,此处不再赘述。
以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围,其均应涵盖在本发明的权利要求和说明书的范围当中。
Claims (10)
1.一种针对双偏振天气雷达的水凝物粒子分类识别方法,其特征在于,所述方法包括:
将预处理后的双偏振观测参量数据输入分类识别模型以训练所述分类识别模型;
所述分类识别模型中的各个神经元都有对应的线性关系系数和偏移,单个神经元
到单个神经元的线性关系为,其中表示为n*1的权重矩阵,表示1*n
的向量矩阵,将整个分类识别模型做一般化,假设第层共有个神经元,则对于第层的
第个神经元的输出,有如(1)所示的前向传播:
所述Relu激活函数如(2)所示:
其中,x表示输入层矩阵的横标志,max()表示最大值函数;
所述Softmax激活函数如(3)所示:
确定所述分类识别模型的多类交叉熵损失函数如(4)所示:
(5)和(6)中,上标L表示第L层神经元,上标L-1表示第L-1层神经元,表示选用的激
活函数,表示L层n*1的矩阵向量,表示输出L-1层n*1的矩阵向量,表示输出层,表示第L层输出层,表示输入层矩阵的竖标志,标识Hadamard积,表示对应位置的对
应数值相乘;
利用所述权值文件配置所述分类识别模型,并将新数据流输入至所述分类识别模型中依次得到识别效果。
2.根据权利要求1所述的方法,其特征在于,在对双偏振观测参量数据进行预处理之前,所述方法还包括:
将各个双偏振观测参量通过线性化的方法转换到预设的阈值范围内,实现对原始数据的等比例缩放,得到所述双偏振观测参量数据。
5.根据权利要求1所述的方法,其特征在于,在所述分类识别模型的训练过程中,对待输入的数据采用K-Fold 多次划分,以利用全部数据集。
6.根据权利要求1-5任一项所述的方法,其特征在于,在所述分类识别模型的训练完成后,对所述分类识别模型进行校验,若学习率低于预设的阈值则对预处理后的双偏振观测参量数据进行正则化,直到学习率达到所述预设的阈值。
7.一种针对双偏振天气雷达的水凝物粒子分类识别装置,其特征在于,所述装置包括处理器,所述处理器被配置为:
将预处理后的双偏振观测参量数据输入分类识别模型以训练所述分类识别模型;
所述分类识别模型中的各个神经元都有对应的线性关系系数和偏移,单个神经元
到单个神经元的线性关系为,其中表示为n*1的权重矩阵,表示1*n
的向量矩阵,将整个分类识别模型做一般化,假设第层共有个神经元,则对于第层
的第个神经元的输出,有如(1)所示的前向传播:
所述Relu激活函数如(2)所示:
其中,x表示x表示输入层矩阵的横标志,max()表示最大值函数;
所述Softmax激活函数如(3)所示:
确定所述分类识别模型的多类交叉熵损失函数如(4)所示:
(5)和(6)中,上标L表示第L层神经元,上标L-1表示第L-1层神经元,表示选用的激
活函数,表示L层n*1的矩阵向量,表示输出L-1层n*1的矩阵向量,表示输出层,
表示第L层输出层,表示输入层矩阵的竖标志,标识Hadamard积,表示对应位置的对应
数值相乘;
利用所述权值文件配置所述分类识别模型,并将新数据流输入至所述分类识别模型中依次得到识别效果。
10.根据权利要求7-9任一项所述的装置,其特征在于,所述处理器被进一步配置为:
在所述分类识别模型的训练完成后,对所述分类识别模型进行校验,若学习率低于预设的阈值则对预处理后的双偏振观测参量数据进行正则化,直到学习率达到所述预设的阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210635178.8A CN114720956B (zh) | 2022-06-07 | 2022-06-07 | 针对双偏振天气雷达的水凝物粒子分类识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210635178.8A CN114720956B (zh) | 2022-06-07 | 2022-06-07 | 针对双偏振天气雷达的水凝物粒子分类识别方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114720956A true CN114720956A (zh) | 2022-07-08 |
CN114720956B CN114720956B (zh) | 2022-11-01 |
Family
ID=82232489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210635178.8A Active CN114720956B (zh) | 2022-06-07 | 2022-06-07 | 针对双偏振天气雷达的水凝物粒子分类识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114720956B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166750A (zh) * | 2022-08-22 | 2022-10-11 | 天津大学 | 一种基于双偏振多普勒雷达数据的定量降水估计方法 |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107967515A (zh) * | 2016-10-19 | 2018-04-27 | 三星电子株式会社 | 用于神经网络量化的方法和设备 |
CN109061646A (zh) * | 2018-06-26 | 2018-12-21 | 成都信息工程大学 | 一种基于降水云识别的多普勒天气雷达定量降水估计方法 |
CN109272046A (zh) * | 2018-09-26 | 2019-01-25 | 北京科技大学 | 基于L2重新正则化Adam切换模拟回火SGD的深度学习方法 |
CN109360206A (zh) * | 2018-09-08 | 2019-02-19 | 华中农业大学 | 基于深度学习的大田稻穗分割方法 |
CN110265065A (zh) * | 2019-05-13 | 2019-09-20 | 厦门亿联网络技术股份有限公司 | 一种构建语音检测模型的方法及语音端点检测系统 |
US20200065677A1 (en) * | 2018-08-24 | 2020-02-27 | Well Data Labs, Inc. | Machine learning assisted events recognition on time series completion data |
CN111028306A (zh) * | 2019-11-06 | 2020-04-17 | 杭州电子科技大学 | 基于AR2 U-Net神经网络的快速磁共振成像方法 |
CN111931422A (zh) * | 2020-08-10 | 2020-11-13 | 中国科学院上海微系统与信息技术研究所 | 基于Adam的优化方法、系统以及终端 |
CN112116921A (zh) * | 2020-09-10 | 2020-12-22 | 南京邮电大学 | 一种基于整合优化器的单声道语音分离方法 |
CN112183763A (zh) * | 2020-09-15 | 2021-01-05 | 广州中国科学院先进技术研究所 | 一种基于混合策略的梯度下降优化方法 |
CN112488772A (zh) * | 2020-12-18 | 2021-03-12 | 四川长虹电器股份有限公司 | 基于改进Ranger的深度学习方法 |
CN112491854A (zh) * | 2020-11-19 | 2021-03-12 | 郑州迪维勒普科技有限公司 | 一种基于fcnn的多方位安全入侵检测方法及系统 |
US20210287089A1 (en) * | 2020-03-14 | 2021-09-16 | DataRobot, Inc. | Automated and adaptive design and training of neural networks |
US20210334626A1 (en) * | 2020-04-28 | 2021-10-28 | Novatel Inc. | Gnss-receiver interference detection using deep learning |
CN113570414A (zh) * | 2021-07-29 | 2021-10-29 | 国网福建省电力有限公司经济技术研究院 | 基于改进Adam算法优化深度神经网络的电价预测方法 |
CN113625283A (zh) * | 2021-07-28 | 2021-11-09 | 南京航空航天大学 | 基于残差卷积神经网络的双极化气象雷达水凝物粒子相态识别方法 |
CN114282713A (zh) * | 2021-12-13 | 2022-04-05 | 浙江树人学院(浙江树人大学) | 基于RAdam-DA-NLSTM时间序列预测方法 |
-
2022
- 2022-06-07 CN CN202210635178.8A patent/CN114720956B/zh active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107967515A (zh) * | 2016-10-19 | 2018-04-27 | 三星电子株式会社 | 用于神经网络量化的方法和设备 |
CN109061646A (zh) * | 2018-06-26 | 2018-12-21 | 成都信息工程大学 | 一种基于降水云识别的多普勒天气雷达定量降水估计方法 |
US20200065677A1 (en) * | 2018-08-24 | 2020-02-27 | Well Data Labs, Inc. | Machine learning assisted events recognition on time series completion data |
CN109360206A (zh) * | 2018-09-08 | 2019-02-19 | 华中农业大学 | 基于深度学习的大田稻穗分割方法 |
CN109272046A (zh) * | 2018-09-26 | 2019-01-25 | 北京科技大学 | 基于L2重新正则化Adam切换模拟回火SGD的深度学习方法 |
CN110265065A (zh) * | 2019-05-13 | 2019-09-20 | 厦门亿联网络技术股份有限公司 | 一种构建语音检测模型的方法及语音端点检测系统 |
CN111028306A (zh) * | 2019-11-06 | 2020-04-17 | 杭州电子科技大学 | 基于AR2 U-Net神经网络的快速磁共振成像方法 |
US20210287089A1 (en) * | 2020-03-14 | 2021-09-16 | DataRobot, Inc. | Automated and adaptive design and training of neural networks |
US20210334626A1 (en) * | 2020-04-28 | 2021-10-28 | Novatel Inc. | Gnss-receiver interference detection using deep learning |
CN111931422A (zh) * | 2020-08-10 | 2020-11-13 | 中国科学院上海微系统与信息技术研究所 | 基于Adam的优化方法、系统以及终端 |
CN112116921A (zh) * | 2020-09-10 | 2020-12-22 | 南京邮电大学 | 一种基于整合优化器的单声道语音分离方法 |
CN112183763A (zh) * | 2020-09-15 | 2021-01-05 | 广州中国科学院先进技术研究所 | 一种基于混合策略的梯度下降优化方法 |
CN112491854A (zh) * | 2020-11-19 | 2021-03-12 | 郑州迪维勒普科技有限公司 | 一种基于fcnn的多方位安全入侵检测方法及系统 |
CN112488772A (zh) * | 2020-12-18 | 2021-03-12 | 四川长虹电器股份有限公司 | 基于改进Ranger的深度学习方法 |
CN113625283A (zh) * | 2021-07-28 | 2021-11-09 | 南京航空航天大学 | 基于残差卷积神经网络的双极化气象雷达水凝物粒子相态识别方法 |
CN113570414A (zh) * | 2021-07-29 | 2021-10-29 | 国网福建省电力有限公司经济技术研究院 | 基于改进Adam算法优化深度神经网络的电价预测方法 |
CN114282713A (zh) * | 2021-12-13 | 2022-04-05 | 浙江树人学院(浙江树人大学) | 基于RAdam-DA-NLSTM时间序列预测方法 |
Non-Patent Citations (7)
Title |
---|
HU J: "An Optimization Strategy Based on Hybrid Algorithm of Adam and SGD", 《IEEE ACCESS》 * |
KE CUI: "Applying Radam Method to improve treatment of Convolutional Neural Network on Banknote Identification", 《2020 INTERNATIONAL CONFERENCE ON COMPUTER ENGINEERING AND APPLICATION (ICCEA)》 * |
SURESHA, H.S: "Alzheimer Disease Detection Based on Deep Neural Network with Rectified Adam Optimization Technique using MRI Analysis", 《PROCEEDINGS OF THE 2020 THIRD INTERNATIONAL CONFERENCE ON ADVANCES IN ELECTRONICS, COMPUTERS AND COMMUNICATIONS (ICAECC)》 * |
冉元波: "双偏振天气雷达水凝物识别研究", 《成都信息工程大学学报》 * |
刘玉良: "《深度学习》", 30 November 2019, 西安:西安电子科技大学出版社 * |
徐建鹏: "基于RAdam卷积神经网络的水稻生育期图像识别", 《农业工程学报》 * |
王金虎: "云雷达联合微波辐射计的水凝物识别研究", 《电波科学学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166750A (zh) * | 2022-08-22 | 2022-10-11 | 天津大学 | 一种基于双偏振多普勒雷达数据的定量降水估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114720956B (zh) | 2022-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hu et al. | A hybrid technique for short-term wind speed prediction | |
CN113239830B (zh) | 一种基于全尺度特征融合的遥感图像云检测方法 | |
CN110503112A (zh) | 一种增强特征学习的小目标检测及识别方法 | |
Awad | An Unsupervised Artificial Neural Network Method for Satellite Image Segmentation. | |
CN107121679A (zh) | 用于雷达回波外推的循环神经网络预测法及记忆单元结构 | |
CN106910176A (zh) | 一种基于深度学习的人脸图像去遮挡方法 | |
US5973635A (en) | Enhanced microburst detection system | |
CN113408423A (zh) | 适用于tx2嵌入式平台的水产品目标实时检测方法 | |
CN104237890B (zh) | 一种由“列车效应”引起的暴雨识别及预报方法 | |
CN114720956B (zh) | 针对双偏振天气雷达的水凝物粒子分类识别方法及装置 | |
CN108229750A (zh) | 一种股票收益率预测方法 | |
CN114419413A (zh) | 感受野自适应的变电站绝缘子缺陷检测神经网络构建方法 | |
Li et al. | Gated auxiliary edge detection task for road extraction with weight-balanced loss | |
CN106485239A (zh) | 一种利用单类支持向量机检测河流目标的方法 | |
CN113239722A (zh) | 基于深度学习在多尺度下的强对流外推方法及系统 | |
CN112613504A (zh) | 一种声呐水下目标检测方法 | |
CN115937659A (zh) | 基于Mask-RCNN的室内复杂环境下多目标检测方法 | |
CN115079117A (zh) | 一种带定位精度预测的三维目标检测方法 | |
CN117233706B (zh) | 一种基于多层通道注意力机制的雷达有源干扰识别方法 | |
CN115393631A (zh) | 基于贝叶斯层图卷积神经网络的高光谱图像分类方法 | |
CN115423810B (zh) | 一种风力发电机组叶片覆冰形态分析方法 | |
CN112556682A (zh) | 一种水下复合传感器目标自动检测算法 | |
Tao et al. | Detection of precipitation cloud over the tibet based on the improved U-net | |
CN112132207A (zh) | 基于多分支特征映射目标检测神经网络构建方法 | |
CN116012618A (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 |