CN114720956B - 针对双偏振天气雷达的水凝物粒子分类识别方法及装置 - Google Patents
针对双偏振天气雷达的水凝物粒子分类识别方法及装置 Download PDFInfo
- Publication number
- CN114720956B CN114720956B CN202210635178.8A CN202210635178A CN114720956B CN 114720956 B CN114720956 B CN 114720956B CN 202210635178 A CN202210635178 A CN 202210635178A CN 114720956 B CN114720956 B CN 114720956B
- Authority
- CN
- China
- Prior art keywords
- layer
- dual
- polarization
- recognition model
- data
- 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
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/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 (8)
1.一种针对双偏振天气雷达的水凝物粒子分类识别方法,其特征在于,所述方法包括:
对双偏振观测参量数据进行预处理,所述双偏振观测参量数据至少包括反射率因子ZH、水凝物粒子差分反射率ZDR、差分相移率KDP以及交叉相关系数ρHV中的一种;
将预处理后的双偏振观测参量数据输入分类识别模型以训练所述分类识别模型;
所述分类识别模型中的各个神经元都有对应的线性关系系数w和偏移b,单个神经元到单个神经元的线性关系为z=∑wixi+b,其中wi表示为n*1的权重矩阵,xi表示1*n的向量矩阵,将整个分类识别模型做一般化,假设第l-1层共有m个神经元,则对于第l层的第j个神经元的输出有如公式(1)所示的前向传播:
所述Relu激活函数如公式(2)所示:
f(x)=max(0,x) (2)
其中,x表示输入层矩阵的横标志,max()表示最大值函数;
所述Softmax激活函数如公式(3)所示:
确定所述分类识别模型的多类交叉熵损失函数如(4)所示:
通过如下公式(5)和公式(6)迭代求解各层神经元的w和b:
公式(5)和公式(6)中,上标L表示第L层神经元,上标L-1表示第L-1层神经元,σ′()表示选用的激活函数,aL表示L层n*1的矩阵向量,aL-1表示输出L-1层n*1的矩阵向量,z表示输出层,zL表示第L层输出层,y表示输入层矩阵的竖标志,⊙标识Hadamard积,表示对应位置的对应数值相乘;
在求解出的各层神经元的w和b都小于停止迭代阈值ε时,输出线性关系系数集W和偏移集B,得到权值文件;
利用所述权值文件配置所述分类识别模型,并将新数据流输入至所述分类识别模型中依次得到识别效果;
2.根据权利要求1所述的方法,其特征在于,在对双偏振观测参量数据进行预处理之前,所述方法还包括:
将各个双偏振观测参量通过线性化的方法转换到预设的阈值范围内,实现对原始数据的等比例缩放,得到所述双偏振观测参量数据。
3.根据权利要求1所述的方法,其特征在于,所述对双偏振观测参量数据进行预处理,包括:
通过如下公式(7)和公式(8)平滑所述双偏振观测参量数据:
g(x,y)=1/n∑I∈NeighbourI(x,y) (7)
g(x,y)=MED{I(x-k,y-l)}(k,l∈W) (8)
公式(7)和公式(8)中,Neighbour为偏振参量数据总集,n为偏振参量数据总集的总数,g(x,y)为输出数据,I(x,y)为输入数据,W为滑动窗口大小。
4.根据权利要求1所述的方法,其特征在于,在所述分类识别模型的训练过程中,对待输入的数据采用K-Fold多次划分,以利用全部数据集。
5.根据权利要求1-4任一项所述的方法,其特征在于,在所述分类识别模型的训练完成后,对所述分类识别模型进行校验,若学习率低于预设的阈值则对预处理后的双偏振观测参量数据进行正则化,直到学习率达到所述预设的阈值。
6.一种针对双偏振天气雷达的水凝物粒子分类识别装置,其特征在于,所述装置包括处理器,所述处理器被配置为:
对双偏振观测参量数据进行预处理,所述双偏振观测参量数据至少包括反射率因子ZH、水凝物粒子差分反射率ZDR、差分相移率KDP以及交叉相关系数ρH中的一种;
将预处理后的双偏振观测参量数据输入分类识别模型以训练所述分类识别模型;
所述分类识别模型中的各个神经元都有对应的线性关系系数w和偏移b,单个神经元到单个神经元的线性关系为z=∑wixi+b,其中wi表示为n*1的权重矩阵,xi表示1*n的向量矩阵,将整个分类识别模型做一般化,假设第l-1层共有m个神经元,则对于第l层的第j个神经元的输出有如公式(1)所示的前向传播:
所述Relu激活函数如公式(2)所示:
f(x)=max(0,x) (2)
其中,x表示x表示输入层矩阵的横标志,max()表示最大值函数;
所述Softmax激活函数如公式(3)所示:
确定所述分类识别模型的多类交叉熵损失函数如公式(4)所示:
通过如下公式(5)和公式(6)迭代求解各层神经元的w和b:
公式(5)和公式(6)中,上标L表示第L层神经元,上标L-1表示第L-1层神经元,σ′()表示选用的激活函数,aL表示L层n*1的矩阵向量,aL-1表示输出L-1层n*1的矩阵向量,z表示输出层,zL表示第L层输出层,y表示输入层矩阵的竖标志,⊙标识Hadamard积,表示对应位置的对应数值相乘;
在求解出的各层神经元的w和b都小于停止迭代阈值ε时,输出线性关系系数集W和偏移集B,得到权值文件;
利用所述权值文件配置所述分类识别模型,并将新数据流输入至所述分类识别模型中依次得到识别效果;
7.根据权利要求6所述的装置,其特征在于,所述处理器被进一步配置为:
通过如下公式(7)和公式(8)平滑所述双偏振观测参量数据:
g(x,y)=1/n∑I∈NeighbourI(x,y) (7)
g(x,y)=MED{I(x-k,y-l)}(k,l∈W) (8)
公式(7)和公式(8)中,Neighbour为偏振参量数据总集,n为偏振参量数据总集的总数,g(x,y)为输出数据,I(x,y)为输入数据,W为滑动窗口大小。
8.根据权利要求7任一项所述的装置,其特征在于,所述处理器被进一步配置为:
在所述分类识别模型的训练完成后,对所述分类识别模型进行校验,若学习率低于预设的阈值则对预处理后的双偏振观测参量数据进行正则化,直到学习率达到所述预设的阈值。
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 CN114720956A (zh) | 2022-07-08 |
CN114720956B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166750A (zh) * | 2022-08-22 | 2022-10-11 | 天津大学 | 一种基于双偏振多普勒雷达数据的定量降水估计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107967515A (zh) * | 2016-10-19 | 2018-04-27 | 三星电子株式会社 | 用于神经网络量化的方法和设备 |
CN111028306A (zh) * | 2019-11-06 | 2020-04-17 | 杭州电子科技大学 | 基于AR2 U-Net神经网络的快速磁共振成像方法 |
CN112491854A (zh) * | 2020-11-19 | 2021-03-12 | 郑州迪维勒普科技有限公司 | 一种基于fcnn的多方位安全入侵检测方法及系统 |
CN113625283A (zh) * | 2021-07-28 | 2021-11-09 | 南京航空航天大学 | 基于残差卷积神经网络的双极化气象雷达水凝物粒子相态识别方法 |
CN114282713A (zh) * | 2021-12-13 | 2022-04-05 | 浙江树人学院(浙江树人大学) | 基于RAdam-DA-NLSTM时间序列预测方法 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109061646B (zh) * | 2018-06-26 | 2022-04-08 | 成都信息工程大学 | 一种基于降水云识别的多普勒天气雷达定量降水估计方法 |
WO2020041755A1 (en) * | 2018-08-24 | 2020-02-27 | Well Data Labs, Inc. | Machine learning assisted events recognition on time series well data |
CN109360206B (zh) * | 2018-09-08 | 2021-11-12 | 华中农业大学 | 基于深度学习的大田稻穗分割方法 |
CN109272046A (zh) * | 2018-09-26 | 2019-01-25 | 北京科技大学 | 基于L2重新正则化Adam切换模拟回火SGD的深度学习方法 |
CN110265065B (zh) * | 2019-05-13 | 2021-08-03 | 厦门亿联网络技术股份有限公司 | 一种构建语音端点检测模型的方法及语音端点检测系统 |
US11334795B2 (en) * | 2020-03-14 | 2022-05-17 | 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 | 广州中国科学院先进技术研究所 | 一种基于混合策略的梯度下降优化方法 |
CN112488772A (zh) * | 2020-12-18 | 2021-03-12 | 四川长虹电器股份有限公司 | 基于改进Ranger的深度学习方法 |
CN113570414A (zh) * | 2021-07-29 | 2021-10-29 | 国网福建省电力有限公司经济技术研究院 | 基于改进Adam算法优化深度神经网络的电价预测方法 |
-
2022
- 2022-06-07 CN CN202210635178.8A patent/CN114720956B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107967515A (zh) * | 2016-10-19 | 2018-04-27 | 三星电子株式会社 | 用于神经网络量化的方法和设备 |
CN111028306A (zh) * | 2019-11-06 | 2020-04-17 | 杭州电子科技大学 | 基于AR2 U-Net神经网络的快速磁共振成像方法 |
CN112491854A (zh) * | 2020-11-19 | 2021-03-12 | 郑州迪维勒普科技有限公司 | 一种基于fcnn的多方位安全入侵检测方法及系统 |
CN113625283A (zh) * | 2021-07-28 | 2021-11-09 | 南京航空航天大学 | 基于残差卷积神经网络的双极化气象雷达水凝物粒子相态识别方法 |
CN114282713A (zh) * | 2021-12-13 | 2022-04-05 | 浙江树人学院(浙江树人大学) | 基于RAdam-DA-NLSTM时间序列预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114720956A (zh) | 2022-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hu et al. | A hybrid technique for short-term wind speed prediction | |
CN107909082B (zh) | 基于深度学习技术的声呐图像目标识别方法 | |
CN102013017B (zh) | 一种高分辨率遥感图像场景粗分类方法 | |
Awad | An Unsupervised Artificial Neural Network Method for Satellite Image Segmentation. | |
CN110503112A (zh) | 一种增强特征学习的小目标检测及识别方法 | |
CN109101897A (zh) | 水下机器人的目标检测方法、系统及相关设备 | |
CN113239830B (zh) | 一种基于全尺度特征融合的遥感图像云检测方法 | |
CN107578061A (zh) | 基于最小化损失学习的不平衡样本分类方法 | |
CN106910176A (zh) | 一种基于深度学习的人脸图像去遮挡方法 | |
CN107529650A (zh) | 网络模型的构建和闭环检测方法、相应装置及计算机设备 | |
CN113486764B (zh) | 一种基于改进的YOLOv3的坑洼检测方法 | |
CN111160268A (zh) | 一种基于多任务学习的多角度sar目标识别方法 | |
CN114720956B (zh) | 针对双偏振天气雷达的水凝物粒子分类识别方法及装置 | |
US8615478B2 (en) | Using affinity measures with supervised classifiers | |
CN107977683A (zh) | 基于卷积特征提取和机器学习的联合sar目标识别方法 | |
CN108229750A (zh) | 一种股票收益率预测方法 | |
CN104156943A (zh) | 基于非支配邻域免疫算法的多目标模糊聚类图像变化检测方法 | |
CN111882114B (zh) | 一种短时交通流量预测模型构建方法及预测方法 | |
CN117233706B (zh) | 一种基于多层通道注意力机制的雷达有源干扰识别方法 | |
CN109543748A (zh) | 一种基于多模型融合的水果识别方法 | |
CN112556682A (zh) | 一种水下复合传感器目标自动检测算法 | |
CN115423810B (zh) | 一种风力发电机组叶片覆冰形态分析方法 | |
CN112132207A (zh) | 基于多分支特征映射目标检测神经网络构建方法 | |
CN116012618A (zh) | 一种基于雷达回波图的天气辨识方法、系统、设备及介质 | |
CN117079120A (zh) | 一种基于改进ga算法的目标识别模型优化方法 |
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 |