阅读说明：本技术 一种多单体对流系统短时强降水事件的自动临近预报方法 (Automatic nowcasting method for short-time strong rainfall event of multi-monomer convection system ) 是由 王萍 王琮 王迪 于 2019-09-25 设计创作，主要内容包括：本发明公开了一种多单体对流系统短时强降水事件的自动临近预报方法,包括以下步骤：对流单体识别和跟踪。使用多阈值自适应算法进行对流单体识别,可以得到同时保留对流单体核心和周边相关信息的对流单体识别结果,并使用光流算法进行对流单体跟踪,可以得到对流单体的速度；多单体对流系统识别。根据多单体对流系统内部各对流单体之间的时空相关性和单体速度,给出了对流单体下一时刻的位置预测,并计算对流单体之间的叠加系数,根据叠加系数建立相关性矩阵,使用传递闭包聚类方法得到多单体对流系统的识别结果；多单体对流系统图模型的构建和短时强降水事件识别。本方法实现了自动的多单体对流系统短时强降水事件临近预报,对灾害进行及时的预警,减少了经济损失和人员伤亡。(The invention discloses an automatic nowcasting method for a short-time strong precipitation event of a multi-monomer convection system, which comprises the following steps of: and identifying and tracking the convection single body. The method comprises the steps of performing convection monomer identification by using a multi-threshold self-adaptive algorithm, obtaining a convection monomer identification result which simultaneously retains core and peripheral related information of a convection monomer, and performing convection monomer tracking by using an optical flow algorithm to obtain the speed of the convection monomer; and identifying the multi-monomer convection system. According to the space-time correlation and the monomer speed between each pair of convection monomers in the multi-monomer convection system, the position prediction of the convection monomers at the next moment is given, the superposition coefficient between the convection monomers is calculated, a correlation matrix is established according to the superposition coefficient, and the identification result of the multi-monomer convection system is obtained by using a transfer closure clustering method; and constructing a multi-monomer convection system graph model and identifying a short-time strong precipitation event. The method realizes the close forecast of the short-time strong rainfall event of the automatic multi-monomer convection system, carries out early warning on disasters in time, and reduces economic loss and casualties.)

1. An automatic nowcasting method for a short-time strong precipitation event of a multi-monomer convection system is characterized by comprising the following steps:

identifying and tracking a convection cell; the method comprises the steps of performing convection monomer identification by using a multi-threshold self-adaptive algorithm, obtaining a convection monomer identification result which simultaneously retains core and peripheral related information of a convection monomer, and performing convection monomer tracking by using an optical flow algorithm to obtain the speed of the convection monomer;

identifying a multi-monomer convection system; according to the space-time correlation and the monomer speed between each pair of convection monomers in the multi-monomer convection system, the position prediction of the convection monomers at the next moment is given, the superposition coefficient between the convection monomers is calculated, a correlation matrix is established according to the superposition coefficient, and the identification result of the multi-monomer convection system is obtained by using a transfer closure clustering method;

constructing a multi-monomer convection system graph model and identifying a short-time strong precipitation event; and establishing a graph model for each multi-monomer convection system, obtaining graph characteristics, and performing the close prediction of the short-time strong precipitation event of the multi-monomer convection system by using a random forest model.

2. The method of automatic nowcasting for short-term heavy precipitation events in multi-cell convection systems as claimed in claim 1, wherein said convection cell identification and tracking specifically comprises the steps of:

setting a threshold T_lowDetermining a convection cell preliminary candidate region; for each point in the radar image, if the reflectivity value of the point is greater than or equal to T_lowSetting the point as a candidate point of the convection monomer; if the reflectivity value of the point is less than T_lowSetting the point as a non-convective monomer point;

setting a plurality of reflectivity threshold values T from large to small_high,T_high-5,T_high-10...T_low]A core for defining a convection cell and a maximum area mesochite containing a single core; for an area, if the reflectivity of each point in the area is greater than the current threshold and the area determined by other larger thresholds is not included in the area, the area is called as a convection monomer core; if an area is larger than the current threshold and only contains one convection single core, and the next level of threshold andthe other regions are contained in one region together, and the region is called a convection monomer mesochite;

for each mesochite area, alternately using an expansion algorithm in image morphology, and sequentially filling peripheral areas except the mesochite area in the monomer candidate area to obtain a final convection monomer identification result;

and for each identified convection monomer, carrying out convection monomer tracking by using an optical flow algorithm, and calculating to obtain the speed of each convection monomer.

3. The method of automatic nowcasting for transient heavy precipitation events in a multi-monomer convection system as claimed in claim 1, wherein said identification of the multi-monomer convection system specifically comprises the steps of:

based on the obtained identification result of the convection single body, for each convection single body, the speed at the current moment is assumed to be v_tThen:

v_t+1＝α⁰v_t+α¹v_t-1+α²v_t-2,(0≤α≤1) (1)

using v_t+1Giving the monomer position at the next moment;

stacking the extrapolation results at a plurality of moments together, and calculating the stacking coefficient eta between the monomers; suppose there are two monomers at time t

superposition coefficient eta_ABCan be used as a monomer

suppose that in a certain volume sweep result of the radar, n single bodies O exist₁,O₂,...,O_n(ii) a Establishing an n x n Boolean relation matrix R_n×n＝[r_ij]_n×nR is a symmetric matrix, R_ijRepresents O_iAnd O_jCorrelation between, given a correlation metric threshold T_re，r_ijSatisfies the relationship:

performing transitive closure clustering on a Boolean matrix R_n×n＝[r_ij]_n×nPerforming p (p is an integer, p is more than or equal to 2) power-Boolean operation to obtain

4. The method of claim 1, wherein the building of the multi-monomer convection system graph model and the identifying of the short-time heavy precipitation event comprise the following steps:

a multi-monomer convection system comprising a plurality of convection monomers, wherein for each convection monomer, a convection monomer characteristic EOC is calculated_L，EOC_S，HOR₃₀，LOR₃₀，VIL_MAX，VIL_AVERRPS; wherein EOC_LThe length of a long axis for fitting the outline ellipse of the convection monomer; EOC_SThe length of a short shaft for fitting the outline ellipse of the convection monomer; HOR₃₀The echo peak height is 30dBZ of the convection monomer; LOR₃₀The echo bottom is 30dBZ for the convection monomer; VIL_MAXCalculating the maximum value of VIL for convection monomers point by point, wherein VIL refers to the vertically accumulated liquid water content and is an estimated value of the total precipitation of each point in the convection monomers; VIL_AVERCalculating the average value after VIL point by point for the convection monomer; RPS is the ratio of the projection of the monomer in the direction of the velocity vector to the monomer velocity;

for a multi-monomer convection system, a graph model is established, and is expressed as G ═ U, E; where U is the set of nodes and E is the set of edges; creating a node u for each convection cell in the system_i，u_iIs the characteristic of the corresponding monomer, will u_iAdding the obtained solution into a set U; for any two points U in U_iAnd u_jAssuming that their velocity vectors are respectively v_iV and v_j(ii) a Calculating u_iTo u_jDistance vector l of_ijAnd is used forRepresents u_iAnd u_jParameter t of relative distance-to-velocity ratio_ij：

And sets a threshold value T_tIf | t_ij|≤T_tThen a connection u is created_iAnd u_jEdge e of_ijAnd use

X_i+1＝D(A+I)X_i,(i＝0,1,2,...) (7)

wherein D is a row normalization coefficient matrix; each time the fusion operation is executed, the node attribute is updated to be the attribute weighted sum of the node and the adjacent nodes, and the weight is the spatial distribution relation between the nodes; performing two times of fusion operation, and taking the two times of fusion result together with the original attribute as the characteristic X of the monomer map_feature＝[X₀,X₁,X₂]；

To X_featureAnd (6) sorting. X_featureIs an n x 21 matrix, and specifies the characteristics VIL which can most directly express local water content information in the monomer characteristics_MAXI.e. X_featureAs a sorting index with X in column 19_featureArranged from large to small, so as to ensure the characteristics of each single graphMiddle, VIL_MAXThe monomer subgraphs with the largest content are arranged at the forefront in total; at the same time in order to ensure X_featureThe number of lines of (A) is fixed, and the ordered X is taken_featureThe first two rows as final features of the simplex graph

training a binary classification model of a random forest by using a large amount of historical meteorological data; for each multi-monomer convection system, one graph model G ═ (U, E) may be obtained; for the graph model, a graph feature representation can be obtained

Technical Field

The invention relates to the field of meteorology, in particular to an automatic nowcasting method for a short-time strong rainfall event of a multi-monomer convection system.

Background

Short-term heavy precipitation (the term "short-term heavy precipitation" as used herein refers to the weather of China nationsPrecipitation events with hourly precipitation greater than 20mm as specified by the central business standards) are one of the most prominent strong convection disasters in china, and the short-term strong precipitation emphasizes the convection property and the short-term duration of precipitation more than the ordinary heavy precipitation^[1]. Because the large precipitation is accumulated in a short time, torrential flood is often formed, and urban waterlogging, torrential flood, debris flow and other hazards are caused. The nowcasting of the short-time strong rainfall has important application value for disaster prevention and disaster control.

The existing approach forecast method for short-time heavy rainfall can be mainly divided into two types: method based on numerical weather forecast^[2]And methods based on radar extrapolation techniques. The method based on the numerical weather forecast describes the change of atmospheric environment through a group of mathematical physics equations, and can solve the atmospheric physical state at any moment under the condition of giving initial conditions through various atmospheric observation technologies, thereby carrying out short-time strong precipitation proximity forecast. Method based on radar extrapolation technology uses Doppler weather radar to perform approach prediction, and mainly comprises centroid tracking method^[3-4]And cross correlation method^[5-6]The methods obtain the change trend of the convection system from the radar data at two adjacent moments, so as to carry out the nowcasting of the future short-time strong precipitation event.

In the process of implementing the invention, the inventor finds that at least the following disadvantages and shortcomings exist in the prior art:

in the methods described in the above documents, a method based on numerical weather prediction can simulate the atmospheric physical change process, but requires more computing resources, is slow in operation speed, is limited by resolution, and cannot describe a convection system with a small scale and a short life cycle; the method based on the radar extrapolation technology operates at a high speed, but cannot simulate complex meteorological processes due to lack of description of an atmospheric physical system. Therefore, the existing method has low recognition rate on the short-time strong precipitation event of the multi-monomer convection system, and cannot accurately predict disasters, thereby causing economic loss and casualties.

[ reference documents ]

[1] Poplar, Sunyun pine, Mausxu, et al. Beijing area short-time heavy precipitation process Multi-Scale circulation characteristics [ J ] Meteorological report, 2016,74(6): 919-.

[2]Sokol Z,V,Zacharov P,et al.Nowcasting of hailstormssimulated by the NWP model COSMO for the area of the Czech Republic[J].Atmospheric research,2016,171:66-76.

[3]Dixon M,Wiener G.TITAN:Thunderstorm identification,tracking,analysis,and nowcasting—A radar-based methodology[J].Journal of atmosphericand oceanic technology,1993,10(6):785-797.

[4]Rossi P J,Chandrasekar V,Hasu V,et al.Kalman filtering–basedprobabilistic nowcasting of object-oriented tracked convective storms[J].Journal of Atmospheric and Oceanic Technology,2015,32(3):461-477.

[5]Johnson J T,MacKeen P L,Witt A,et al.The storm cell identificationand tracking algorithm:An enhanced WSR-88D algorithm[J].Weather andforecasting,1998,13(2):263-276.

[6]Liu Y,Xi D G,Li Z L,et al.A new methodology for pixel-quantitativeprecipitation nowcasting using a pyramid Lucas Kanade optical flow approach[J].Journal of Hydrology,2015,529:354-364.

Disclosure of Invention

In order to solve the problems in the prior art, the invention provides an automatic close prediction method for a short-time strong rainfall event of a multi-monomer convection system, and solves the problems that in the prior art, the recognition rate of the short-time strong rainfall event of the multi-monomer convection system is low, disaster prediction cannot be accurately carried out, and economic loss and casualties are caused.

In order to solve the technical problem, the invention provides an automatic nowcasting method for a short-time strong precipitation event of a multi-monomer convection system, which comprises the following steps:

1) and identifying and tracking the convection single body. The method comprises the steps of performing convection monomer identification by using a multi-threshold self-adaptive algorithm, obtaining a convection monomer identification result which simultaneously retains core and peripheral related information of a convection monomer, and performing convection monomer tracking by using an optical flow algorithm to obtain the speed of the convection monomer; 2) and identifying the multi-monomer convection system. According to the space-time correlation and the monomer speed between each pair of convection monomers in the multi-monomer convection system, the position prediction of the convection monomers at the next moment is given, the superposition coefficient between the convection monomers is calculated, a correlation matrix is established according to the superposition coefficient, and the identification result of the multi-monomer convection system is obtained by using a transfer closure clustering method; 3) and constructing a multi-monomer convection system graph model and identifying a short-time strong precipitation event. And establishing a graph model for each multi-monomer convection system, obtaining graph characteristics, and performing the close prediction of the short-time strong precipitation event of the multi-monomer convection system by using a random forest model.

The method comprises the following steps: identifying and tracking a convection cell;

1-1) setting the threshold T_lowAnd determining a convection cell preliminary candidate region. For each point in the radar image, if the reflectivity value of the point is greater than or equal to T_lowSetting the point as a candidate point of the convection monomer; if the reflectivity value of the point is less than T_lowSetting the point as a non-convective monomer point;

1-2) setting a plurality of reflectivity threshold values T from large to small_high,T_high-5,T_high-10...T_low]The shell is used for determining the core of the convection current monomer and the maximum area containing the single core. For an area, if the reflectivity of each point is greater than the current threshold and no other area determined by the greater threshold is included, the area is called a convection cell core. If one area is larger than the current threshold and only contains one convection monomer core, and the next-level threshold and other areas are jointly contained in one area, the area is called as a convection monomer mesochite;

1-3) alternately using an expansion algorithm in image morphology for each mesochite region, and sequentially filling peripheral regions except the mesochite region in the monomer candidate region to obtain a final convective monomer identification result;

1-4) for each identified convection monomer, carrying out convection monomer tracking by using an optical flow algorithm, and calculating to obtain the speed of each convection monomer;

step two: identifying a multi-monomer convection system;

2-1) based on the convection single body identification result obtained in the step one, for each convection single body, assuming that the speed at the current moment is v_tThen:

v_t+1＝α⁰v_t+α¹v_t-1+α²v_t-2,(0≤α≤1) (1)

using v_t+1Giving the monomer position at the next moment;

2-2) the extrapolation results at a plurality of moments are superposed together, and the superposition coefficient eta between the monomers is calculated. Suppose there are two monomers at time t

Andthe extrapolation results at their 9 time points are respectively

And

then calculate:

superposition coefficient eta_ABCan be used as a monomer

And

the correlation metric of (a) is applied in the subsequent multi-monomer convection system identification;

2-3) suppose that in a certain volume sweep result of the radar, there are n monomers O₁,O₂,...,O_n. Establishing an n x n Boolean relation matrix R_n×n＝[r_ij]_n×nR is a symmetric matrix, R_ijRepresents O_iAnd O_jCorrelation between, given a correlation metric threshold T_re，r_ijSatisfies the relationship:

2-4) performing transitive closure clustering on the Boolean matrix R_n×n＝[r_ij]_n×nPerforming p (p is an integer, p is more than or equal to 2) power-Boolean operation to obtain

If it is

Is a full 1 matrix or

And finishing clustering. Otherwise, repeating the step 2-3;

2-5) relationship matrix

All the elements of 1 in each row in the system are a multi-monomer convection system, and repeated results in the system are removed, so that a multi-monomer convection system identification result taking the superposition relationship as monomer similarity measurement can be obtained (each result may contain a plurality of convection monomers or only one convection monomer);

step three, constructing a multi-monomer convection system graph model and identifying a short-time strong precipitation event;

3-1) one multimonomer comprising a plurality of convective monomersA convection system for calculating a convection cell characteristic EOC for each convection cell_L，EOC_S，HOR₃₀，LOR₃₀，VIL_MAX，VIL_AVERAnd RPS. Wherein EOC_LThe length of a long axis for fitting the outline ellipse of the convection monomer; EOC_SThe length of a short shaft for fitting the outline ellipse of the convection monomer; HOR₃₀The echo peak height is 30dBZ of the convection monomer; LOR₃₀The echo bottom is 30dBZ for the convection monomer; VIL_MAXCalculating the maximum value of VIL for convection monomers point by point, wherein VIL refers to the vertically accumulated liquid water content and is an estimated value of the total precipitation of each point in the convection monomers; VIL_AVERCalculating the average value after VIL point by point for the convection monomer; RPS is the ratio of the projection of the monomer in the direction of the velocity vector to the monomer velocity;

3-2) for a multi-monomer convection system, a graph model was built, denoted G ═ U, E. Where U is the set of nodes and E is the set of edges. Creating a node u for each convection cell in the system_i，u_iIs the characteristic of the corresponding monomer, will u_iAnd adding the data into the set U. For any two points U in U_iAnd u_jAssuming that their velocity vectors are respectively v_iV and v_j. Calculating u_iTo u_jDistance vector l of_ijAnd is used to represent u_iAnd u_jParameter t of relative distance-to-velocity ratio_ij：

And sets a threshold value T_tIf | t_ij|≤T_tThen a connection u is created_iAnd u_jEdge e of_ijAnd use

As e_ij(ii) an attribute of (d);

3-3) assume that there is a single graph G ═ (U, E), G containing n nodes, each node having attributes that are a 7-dimensional vector. The adjacency matrix of the cell graph is denoted by a,a is an n × n symmetric array containing information of all sides of the monomer map if u_iAnd u_jThere is an edge e between_ijThen A_ij＝e_ijOtherwise A_ij0. Use of

An initial node attribute matrix representing a node simplex graph, each row representing a node. Defining a node fusion operation:

X_i+1＝D(A+I)X_i,(i＝0,1,2,...) (7)

where D is a row normalization coefficient matrix. And updating the node attribute into the attribute weighted sum of the node and the adjacent nodes every time the fusion operation is executed, wherein the weight is the spatial distribution relation between the nodes. Performing two times of fusion operation, and taking the two times of fusion result together with the original attribute as the characteristic X of the monomer map_feature＝[X₀,X₁,X₂]。

3-4) to X_featureAnd (6) sorting. X_featureIs an n x 21 matrix, and specifies the characteristics VIL which can most directly express local water content information in the monomer characteristics_MAXI.e. X_featureAs a sorting index with X in column 19_featureArranged from large to small, thus ensuring that among the characteristics of each monomer map, VIL_MAXThe monomer subgraphs with the highest content are arranged foremost in total. At the same time in order to ensure X_featureThe number of lines of (A) is fixed, and the ordered X is taken_featureThe first two rows as final features of the simplex graph

If it is

If the amount is less than two lines, the amount is 0.

3-5) training a binary model of the random forest by using a large amount of historical meteorological data. For each multi-monomer convection system, one graph model G ═ (U, E) can be obtained. For the graph model, a graph feature representation can be obtained

Will be provided with

And inputting a random forest model, and carrying out the nowcasting of the short-time strong rainfall event of the multi-monomer convection system.

Compared with the prior art, the invention has the beneficial effects that:

the technical scheme provided by the invention has the beneficial effects that: the method comprises the steps of modeling a multi-monomer convection system by using a graph model, converting physical characteristics and spatial distribution of the multi-monomer convection system into the graph model and graph characteristics, and carrying out short-time strong precipitation event close prediction on the graph characteristics obtained by the multi-monomer convection system by using a random forest model. The invention realizes high-quality automatic nowcasting of short-time strong rainfall of the multi-monomer convection system, and is beneficial to timely forecasting of weather disasters so as to reduce economic loss and casualties; and the effectiveness of the method is verified through experiments.

Drawings

FIG. 1 is a flow chart of an automatic nowcasting method for a short-term heavy precipitation event in a multi-cell convection system according to the present invention;

FIG. 2 is a schematic diagram of a convective cell identification algorithm;

fig. 3 is a schematic view of a multi-monomer convection system identification. Obtaining the speed of convection monomers by a convection monomer tracking method, predicting the positions of the convection monomers, calculating the superposition coefficient between the convection monomers, establishing a relation Boolean matrix, and obtaining the identification result of the multi-monomer convection system by transfer closure clustering;

FIG. 4 is a schematic diagram of a calculation method of physical characteristics RPS of a convection current monomer;

FIG. 5 is a schematic diagram of a multi-monomer convection system chart model construction method.

Detailed Description

The technical solutions of the present invention are further described in detail with reference to the accompanying drawings and specific embodiments, which are only illustrative of the present invention and are not intended to limit the present invention.

The invention provides an automatic nowcasting method for a short-time strong rainfall event of a multi-monomer convection system, which is designed according to the following design concept: and performing convection monomer identification by using a multi-threshold mathematical morphology-based method to obtain a convection monomer area. Tracking the single convection by using an optical flow algorithm to obtain the speed of the single convection; predicting the positions of the convection monomers by using the obtained convection monomer speed, calculating the correlation among the convection monomers by using the time-space superposition of the convection monomers in the future, and identifying a multi-monomer convection system by using transfer closure clustering; according to the space-time characteristics of the multi-monomer convection system, a graph model is constructed, monomers in the multi-monomer convection system are mapped to be nodes in the graph model, and the physical characteristics of the convection monomers are set to be attributes of the corresponding nodes. The edges in the graph model are used to describe the spatial distribution of the multi-monomer convection system. And generating a graph characteristic which is not influenced by the node sequence for each multi-monomer convection system corresponding to a graph model. 1079 cases of historical data were used to train a random forest model for the imminent prediction of short-term heavy precipitation events in a multi-modal convection system.

As shown in fig. 1, the method mainly includes convection single body identification and tracking, multi-single body convection system identification, multi-single body convection system graph model construction and short-time strong precipitation event identification; the specific contents are as follows:

the method comprises the following steps: identifying and tracking a convection cell;

1-1) setting the threshold T_lowAnd determining a convection cell preliminary candidate region. For each point in the radar image, if the reflectivity value of the point is greater than or equal to T_lowSetting the point as a candidate point of the convection monomer; if the reflectivity value of the point is less than T_lowSetting the point as a non-convective monomer point; as shown in step 1 in FIG. 2, for an area in the Doppler weather radar reflectivity image, let T be_lowEach point in the traversal area is set to be a convective monomer candidate area for points having a reflectivity value greater than or equal to 35dBZ, and the other points are set to be non-convective monomer candidate areas. For theIn step 1 of fig. 2, a region with a minimum reflectance value of 30dBZ may be obtained, and 2 convection monomer candidate regions with a minimum reflectance value of 35dBZ included therein may be obtained;

1-2) setting a plurality of reflectivity threshold values T from large to small_high,T_high-5,T_high-10...T_low]The shell is used for determining the core of the convection current monomer and the maximum area containing the single core. For an area, if the reflectivity of each point is greater than the current threshold and no other area determined by the greater threshold is included, the area is called a convection cell core. If an area is larger than the current threshold and only contains one convection monomer core, and the next level threshold and other areas are jointly contained in an area, the area is called a convection monomer mesochite. As shown in step 2 of fig. 2, let the threshold T be used for two convection current monomer candidate regions_highThe lower core in the first candidate region may be determined at 50dBZ, and the upper core in the first candidate region and the two cores in the second candidate region may not be determined at this threshold. When the threshold is T_highAt-5-45 dBZ, the kernel above the first candidate region and the two kernels in the second candidate region may be determined. Since at the next level of threshold, the two core regions in the second candidate region will be encompassed by the same region, the two cores in the second candidate region are also two mesochites at the same time. When the threshold is T_highAt-10-40 dBZ, the mesochites of the two convective cell cores in the first candidate region may be determined;

1-3) alternately using an expansion algorithm in image morphology for each mesochite region, and sequentially filling peripheral regions except the mesochite region in the monomer candidate region to obtain a final convection monomer identification result. As in step 3 of fig. 2, the two candidate regions are each divided into two convection monomer regions according to the mesochite thereof;

1-4) for each identified convection monomer, carrying out convection monomer tracking by using an optical flow algorithm, and calculating to obtain the speed of each convection monomer;

step two: identifying a multi-monomer convection system;

2-1) based on the convection single body identification result obtained in the step one, for each convection single body, assuming that the speed at the current moment is v_tThen:

v_t+1＝α⁰v_t+α¹v_t-1+α²v_t-2,(0≤α≤1) (1)

using v_t+1The monomer position at the next moment is given. As in the step of extrapolating the positions of the convection monomers in fig. 3, the current speed of the convection monomer is used to calculate the speed of the convection monomer at the next moment, so as to obtain the position prediction of the convection monomer at the next moment, which is represented by a dashed outline, and so on to obtain the position prediction results of the convection monomer at the subsequent 9 moments;

2-2) the extrapolation results at a plurality of moments are superposed together, and the superposition coefficient eta between the monomers is calculated. Suppose there are two monomers at time t

And

the extrapolation results at their 9 time points are respectively

And

then calculate:

superposition coefficient eta_ABCan be used as a monomer

And

the correlation metric of (a) is applied in the subsequent multi-monomer convection system identification. As shown in fig. 3, in the step of calculating the correlation between the convection monomers, the prediction result of the convection monomer position obtained in step (2-2) can obtain the overlapping description between the monomers at the next 9 times, which is shown as the colored filling area at the intersection of the dashed outlines in the figure, and this overlapping description can be used to calculate the overlap coefficient η;

2-3) suppose that in a certain volume sweep result of the radar, there are n monomers O₁,O₂,...,O_n. Establishing an n x n Boolean relation matrix R_n×n＝[r_ij]_n×nR is a symmetric matrix, R_ijRepresents O_iAnd O_jCorrelation between, given a correlation metric threshold T_re，r_ijSatisfies the relationship:

as shown in the correlation calculation step of convection monomers in fig. 3, 4 convection monomers form a 4 × 4 boolean relationship matrix, and the correlation coefficient η between the monomers is calculated and compared with the threshold T_reAs can be seen by comparison, there is a correlation between the convection monomers 1 and 2, and a correlation between 2 and 3. So r in the Boolean relationship matrix₁₂＝r₂₁＝r₂₃＝r₃₂1, others are 0;

2-4) performing transitive closure clustering on the Boolean matrix R_n×n＝[r_ij]_n×nPerforming p (p is an integer, p is more than or equal to 2) power-Boolean operation to obtain

If it is

Is a full 1 matrix or

And finishing clustering. Otherwise, repeating the step (2-4) until p is equal to p + 1;

2-5) relationship matrixAll the elements of 1 in each row in the multi-monomer convection system are a multi-monomer convection system, and repeated results in the multi-monomer convection system are removed, so that a multi-monomer convection system identification result taking the superposition relationship as monomer similarity measurement can be obtained (each result may contain a plurality of convection monomers or only one convection monomer). As shown in the convection system identification step in fig. 3, after the boolean relationship matrix is subjected to multiple exponentiation operations of the transitive closure clustering, the 1 st, 2 nd, and 3 rd rows in the relationship matrix are the same, which indicates that the convection monomers 1, 2, and 3 are a convection system. The 4 th row in the relation matrix indicates that the convection monomer 4 is a convection system;

step three, constructing a multi-monomer convection system graph model and identifying a short-time strong precipitation event;

3-1) a multi-monomer convection system comprising a plurality of convection monomers, wherein for each convection monomer, a convection monomer characteristic EOC is calculated_L，EOC_S，HOR₃₀，LOR₃₀，VIL_MAX，VIL_AVERAnd RPS. Wherein EOC_LThe length of a long axis for fitting the outline ellipse of the convection monomer; EOC_SThe length of a short shaft for fitting the outline ellipse of the convection monomer; HOR₃₀The echo peak height is 30dBZ of the convection monomer; LOR₃₀The echo bottom is 30dBZ for the convection monomer; VIL_MAXCalculating the maximum value of VIL for convection monomers point by point, wherein VIL refers to the vertically accumulated liquid water content and is an estimated value of the total precipitation of each point in the convection monomers; VIL_AVERCalculating the average value after VIL point by point for the convection monomer; the RPS is the ratio of the projection of the monomer in the direction of the velocity vector to the size of the monomer velocity. As shown in FIG. 4, for the calculation method of RPS, for a convective monomer, the projected length C in the velocity direction is calculated_pro，C_proThe ratio of the velocity of the convection monomer and the velocity | v | is RPS;

3-2) for a multi-monomer convection system, a graph model was built, denoted G ═ U, E. Where U is the set of nodes and E is the set of edges. Creating a node u for each convection cell in the system_i，u_iIs the characteristic of the corresponding monomer, will u_iAnd adding the data into the set U. For any two points U in U_iAnd u_jAssuming that their velocity vectors are respectively v_iV and v_j. Calculating u_iTo u_jDistance vector l of_ijAnd is used to represent u_iAnd u_jParameter t of relative distance-to-velocity ratio_ij：

And sets a threshold value T_tIf | t_ij|≤T_tThen a connection u is created_iAnd u_jEdge e of_ijAnd use

As e_ijThe attribute of (2). As shown in fig. 5, for a multi-monomer convection system composed of 4 convection monomers, a graph model is created, each convection monomer corresponds to a node in the graph model, and the physical characteristics of each convection monomer are calculated as the attributes of the corresponding node. Calculate t between monomers, as shown, t₁₂,t₁₄,t₁₃＞T_tSo edges are established between nodes 1 and 2, 1 and 4, 1 and 3, respectively, and will be

As an attribute of the edge connecting nodes i and j. t is t₂₃,t₃₄,t₂₄＜T_tSo no edges are established between nodes 2 and 3, 3 and 4, and 2 and 4.

3-3) assume that there is a single graph G ═ (U, E), G containing n nodes, each node having attributes that are a 7-dimensional vector. A is used for representing the adjacent matrix of the single map, and A is an n multiplied by n symmetric matrix and contains the information of all sides of the single map if u_iAnd u_jThere is an edge e between_ijThen A_ij＝e_ijOtherwise A_ij0. Use ofAn initial node attribute matrix representing a node simplex graph, each row representing a node. Defining a node fusion operation:

X_i+1＝D(A+I)X_i,(i＝0,1,2,...) (7)

where D is a row normalization coefficient matrix. And updating the node attribute into the attribute weighted sum of the node and the adjacent nodes every time the fusion operation is executed, wherein the weight is the spatial distribution relation between the nodes. Performing two times of fusion operation, and taking the two times of fusion result together with the original attribute as the characteristic X of the monomer map_feature＝[X₀,X₁,X₂]。

3-4) to X_featureAnd (6) sorting. X_featureIs an n x 21 matrix, and specifies the characteristics VIL which can most directly express local water content information in the monomer characteristics_MAXI.e. X_featureAs a sorting index with X in column 19_featureArranged from large to small, thus ensuring that among the characteristics of each monomer map, VIL_MAXThe monomer subgraphs with the highest content are arranged foremost in total. At the same time in order to ensure X_featureThe number of lines of (A) is fixed, and the ordered X is taken_featureThe first two rows as final features of the simplex graphIf it is

If the amount is less than two lines, the amount is 0.

3-5) training a binary model of the random forest by using a large amount of historical meteorological data. For each multi-monomer convection system, one graph model G ═ (U, E) can be obtained. For the graph model, a graph feature representation can be obtained

Will be provided with

And inputting a random forest model, and carrying out the nowcasting of the short-time strong rainfall event of the multi-monomer convection system.

The feasibility of the automatic nowcasting method for the short-time strong precipitation event of the multi-monomer convection system provided by the embodiment of the invention is verified by the following specific tests, which are described in detail as follows:

the data of the effectiveness of the test method is provided by the atmospheric detection center of the China weather service bureau. Specifically, the data are 5S-band Doppler weather radar data from Binzhou, Jinan, Qingdao, Shijiazhuang and Weifang, and weather automatic station data covered by the detection ranges of the radar. The time ranges were 6,7, 8, 9 months of 2015 and 2016. These data were collated to give 1349 weather samples, of which 590 short-time heavy precipitation events and 759 non-short-time heavy precipitation events. 1079 of these samples were used for the adjustment of various parameter parameters in the method and 270 were used for the final performance evaluation of the method. As shown in table 1, the final performance of the algorithm is described using three indexes commonly used in meteorology, namely, hit rate (POD), false positive rate (FAR), and Critical Success Index (CSI).

TABLE 1 evaluation results of examples of the present invention

While the present invention has been described with reference to the accompanying drawings, the present invention is not limited to the above-described embodiments, which are illustrative only and not restrictive, and various modifications which do not depart from the spirit of the present invention and which are intended to be covered by the claims of the present invention may be made by those skilled in the art.