CN113191299B

CN113191299B - Vortex identification method and device, storage medium and electronic equipment

Info

Publication number: CN113191299B
Application number: CN202110525651.2A
Authority: CN
Inventors: 汤欢; 傅慎明; 张敬萍; 孙建华; 李万莉; 张元春
Original assignee: Institute of Atmospheric Physics of CAS
Current assignee: Institute of Atmospheric Physics of CAS
Priority date: 2021-05-14
Filing date: 2021-05-14
Publication date: 2022-02-18
Anticipated expiration: 2041-05-14
Also published as: CN113191299A

Abstract

The invention provides a vortex identification method, a vortex identification device, a storage medium and electronic equipment, wherein the method comprises the following steps: a valid grid point which may be the vortex center is selected based on wind field data, the vortex center is determined based on the valid vorticity under the condition that each mode meets the condition of anticlockwise rotation by dividing the area around the valid grid point into a plurality of quadrants. The method only selects part of grid points as effective grid points for identification, thereby reducing the processing amount; the identification is carried out based on a plurality of quadrants around the effective lattice point, when the quadrant average wind field limiting condition and the quadrant group average wind direction anticlockwise rotation condition are all met, the identification of the undetermined lattice point at the vortex center is carried out, and the vortex center can be very accurately determined. The method can be applied to identifying the mesoscale vortexes in different shapes and scales, and also can be applied to identifying the mesoscale vortexes based on wind field data with different spatial resolutions; in particular, it is applicable to the identification of mesoscale vortices on complex wind farms.

Description

Vortex identification method and device, storage medium and electronic equipment

Technical Field

The invention relates to the technical field of vortex identification, in particular to a vortex identification method and device, a storage medium and electronic equipment.

Background

The mesoscale vortex plays an important role in balancing global momentum, water vapor and energy and is an important system for causing strong convection or disastrous weather, so that the accurate vortex automatic identification and tracking algorithm is favorable for early warning and prediction of the disastrous weather.

At present, the identification of the vortex at home and abroad mainly comprises the following methods:

(1) the vortex identification method based on manual work comprises the following steps: the vortex is directly judged manually, and the method has the defects that the identification is very subjective, the vortex boundary cannot be quantitatively given, the time consumption is long, and the labor cost is high.

(2) The vortex identification method based on the height field or the air pressure field comprises the following steps: determining a vortex boundary by searching a closed contour line of the outermost layer of the height field or the air pressure field, and identifying a vortex by combining a local minimum value; for example, reference may be made in particular to Lin Z, Guo W, Jia L, et al.Climatiogy of Tibet plate auo viruses derived from multiple reactions datasets [ J ]. Climate Dynamics,2020,55 (9). The method has the defects of large calculation amount, and the small-scale vortex is easy to miss report because some vortex systems do not have medium and low pressure characteristics.

(3) The vortex identification method based on the vorticity field comprises the following steps: filtering the vorticity field or taking the Laplacian term of vorticity as an identification index, and determining a local maximum value to determine a vortex center; for example, reference may be made in particular to the documents Curio J, Chen Y, Schiemann R, Turner A G, Wong K C, Hodges K and Li Y. company of manual and an automated tracking method for Tibet plant esses volts. adv. Atmos. Sci,2018,35: 965-80. The method has the disadvantages that the vorticity field is not completely matched with the wind field, a large amount of false vortices are easily identified, and a large amount of empty reports are caused.

(4) The vortex identification method based on the wind field comprises the following steps: identifying a vortex center based on automatically detecting a vortex flow pattern; for example, reference may be made in particular to Hou J, Wang P, Zhuang S.A New Method of manipulating flows of questions and Detecting the Centers of questions in a Numerical Wind Field [ J ]. Journal of Atmospheric & Oceanic Technology,2015,34(1): 101-. The method has the defects of large calculation amount, complex detection flow and high empty report rate.

(5) The vortex identification method based on deep learning comprises the following steps: learning from the sample set by using a convolutional neural network model, a pyramid scene analysis network model and the like; among them, the convolutional neural network can be referred to as: chinese, 110097075[ P ] 2019-08-06; the pyramid scene analysis network model can refer to documents: zhang Weimin, Yinheyuan, Daihai 29800, and the like, a deep learning-based mesoscale vortex identification method is disclosed, wherein 111767827[ P ] 2020-10-13. The method has the advantage that rules can be directly learned from data without being distinguished and restricted by physical rules. The method has the disadvantages that a large amount of labeled data is needed, modeling is greatly influenced by data quality, and a large amount of parameter adjustment is needed.

The objective identification method in the prior art has the defects of large calculation amount and easy empty report and missed report of vortex center identification, and the manual identification has the defects of large time cost and incapability of determining vortex boundaries.

Disclosure of Invention

In order to solve the above problems, embodiments of the present invention provide a method and an apparatus for identifying a vortex, a storage medium, and an electronic device.

In a first aspect, an embodiment of the present invention provides a method for identifying a vortex, including:

acquiring wind field data in the whole area, selecting effective grid points from the grid points with positive relative vorticity in the whole area, and taking the relative vorticity of the effective grid points as effective vorticity;

determining a target area with the effective grid point as the center, dividing the target area into N quadrants, and determining the wind parameter of each quadrant according to the wind field data; the wind parameters comprise an average wind direction, and the wind parameters further comprise at least one of an average latitudinal wind and an average latitudinal wind;

judging whether the average latitudinal wind or the average latitudinal wind of the quadrant meets the quadrant average wind field limiting condition; forming the quadrant and other quadrants adjacent to the quadrant into a quadrant group, and judging whether the quadrant group meets the condition of anticlockwise rotation of the average wind direction of the quadrant group according to the average wind direction of the quadrant and the average wind direction of the other quadrants adjacent to the quadrant;

taking the effective lattice point as a lattice point to be determined under the condition that m quadrants meet the quadrant mean wind field limiting condition and n quadrant groups meet the quadrant group mean wind direction anticlockwise rotation condition;

and determining the vortex center of the target vortex according to the undetermined point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined points.

In a second aspect, an embodiment of the present invention further provides a device for identifying a vortex, including:

the acquiring module is used for acquiring wind field data in the whole area, selecting effective lattice points from lattice points with positive relative vorticity in the whole area, and taking the relative vorticity of the effective lattice points as effective vorticity;

the dividing module is used for determining a target area taking the effective grid point as a center, dividing the target area into N quadrants, and determining the wind parameter of each quadrant according to the wind field data; the wind parameters comprise an average wind direction, and the wind parameters further comprise at least one of an average latitudinal wind and an average latitudinal wind;

the judging module is used for judging whether the average latitudinal wind or the average latitudinal wind of the quadrant meets the quadrant average wind field limiting condition; forming the quadrant and other quadrants adjacent to the quadrant into a quadrant group, and judging whether the quadrant group meets the condition of anticlockwise rotation of the average wind direction of the quadrant group according to the average wind direction of the quadrant and the average wind direction of the other quadrants adjacent to the quadrant;

the processing module is used for taking the effective lattice point as a lattice point to be determined under the condition that m quadrants meet the quadrant mean wind field limiting condition and n quadrant groups meet the quadrant group mean wind direction anticlockwise rotation condition; and determining the vortex center of the target vortex according to the undetermined point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined points.

In a third aspect, an embodiment of the present invention further provides a computer storage medium, where computer-executable instructions are stored, and the computer-executable instructions are used in the method for identifying a vortex in any one of the foregoing embodiments.

In a fourth aspect, an embodiment of the present invention further provides an electronic device, including:

at least one processor; and the number of the first and second groups,

a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,

the memory stores instructions executable by the at least one processor to enable the at least one processor to perform any of the above described methods of vortex identification.

In the solution provided by the foregoing first aspect of the embodiment of the present invention, an effective grid point that may be a vortex center is selected based on wind field data, an area around the effective grid point is divided into a plurality of quadrants by two ways, and the vortex center is determined based on the effective vorticity under the condition that each way satisfies a counterclockwise rotation condition. The method only selects part of grid points as effective grid points for identification, thereby reducing the processing amount; when the quadrants around the effective lattice point meet the quadrant average wind field limiting condition and the quadrant group average wind direction anticlockwise rotation condition, the vortex center is identified, and the vortex center can be determined very accurately. Through inspection and evaluation, the method has strong universality, can be applied to identifying mesoscale vortexes in different shapes and different scales, and can also be applied to identifying the mesoscale vortexes based on wind field data with different spatial resolutions; particularly, the method can be applied to identification of mesoscale vortexes (particularly plateau vortexes) on a complex wind field, and can also have certain reference significance for identification of vortexes on the ocean.

In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.

FIG. 1 is a flow chart illustrating a method for identifying vortices provided by an embodiment of the present invention;

FIG. 2 is a schematic diagram of an effective grid provided by an embodiment of the present invention;

FIG. 3a is a schematic diagram of a target area provided by an embodiment of the present invention;

FIG. 3b shows another schematic view of a target area provided by an embodiment of the present invention;

FIG. 3c is a schematic illustration of a further target area provided by an embodiment of the present invention;

FIG. 4a is a schematic diagram of four quadrant groups, wherein the arc represents the counterclockwise rotation of the average wind direction of the four quadrant groups;

FIG. 4b is a schematic diagram of the eight-quadrant division provided by the embodiment of the present invention, wherein the circular arc represents that the average wind direction of the eight-quadrant group rotates counterclockwise;

FIG. 5 is a schematic diagram illustrating the relative position between the centers of two vortices provided by an embodiment of the present invention;

FIG. 6 shows a schematic diagram of the azimuthal division provided by an embodiment of the present invention;

FIG. 7a is a schematic diagram of vertical tracking provided by the embodiment of the present invention (taking 450-550hPa flow field diagram as an example);

FIG. 7b is a diagram illustrating the plateau vortex identification effect provided by the embodiment of the present invention;

FIG. 8 is a schematic structural diagram illustrating a vortex identification device according to an embodiment of the present invention;

fig. 9 is a schematic structural diagram of an electronic device for performing a vortex identification method according to an embodiment of the present invention.

Detailed Description

In the description of the present invention, it is to be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", and the like, indicate orientations and positional relationships based on those shown in the drawings, and are used only for convenience of description and simplicity of description, and do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, should not be considered as limiting the present invention.

Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of the present invention, "a plurality" means two or more unless specifically defined otherwise.

Referring to fig. 1, the method for identifying a vortex provided by the embodiment of the present invention includes steps 101-105:

step 101: acquiring wind field data in the whole area, selecting effective grid points from the grid points with positive relative vorticity in the whole area, and taking the relative vorticity of the effective grid points as the effective vorticity.

In the embodiment of the invention, the whole area is an area needing to identify the vortex, such as a Qinghai-Tibet plateau area, a middle and lower area of Yangtze river and the like; the method is mainly used for identifying the middle alpha scale vortex such as plateau vortex, southwest vortex and the like, namely the vortex with the horizontal scale of 200km-2000 km. The wind field data comprises meteorological data of a plurality of grid points, and each grid point corresponds to a group of data, including longitude and latitude, warp wind, weft wind, temperature and the like. Specifically, the wind field data may be fifth generation ECMWF (European Centre for Medium-Range Weather projections, mid-European Weather forecast center) atmospheric re-analysis global climate data (ERA5) with a data resolution of 0.25 ° (25km), i.e. a spacing between two adjacent grid points of 25 km. If plateau vortexes of a complex wind field are identified, smoothing method preprocessing can be carried out on wind field data to reduce the influence of noise.

Because the vorticity field is not completely matched with the wind field, the vortex center is difficult to be accurately determined only based on the relative vorticity; in the embodiment of the invention, partial lattice points, namely effective lattice points, are selected from the whole area, and then the vortex center is determined based on the relative vorticity of the effective lattice points.

Optionally, the step 101 "selecting an effective grid point from grid points with positive relative vorticity in the whole area, and taking the relative vorticity of the effective grid point as the effective vorticity" includes:

step A1: a plurality of grid points adjacent to a grid point (i, j) with positive relative vorticity are determined, and the adjacent grid points at least comprise an upper grid point (i, j +1), a lower grid point (i, j-1), a right grid point (i +1, j) and a left grid point (i-1, j).

In the embodiment of the present invention, the relative vorticity of each grid point may be determined first, specifically, the relative vorticity may be determined based on the central difference, and then it is determined whether the grid point with positive relative vorticity is a valid grid point. Specifically, on the same isobaric surface, the target area is two-dimensional, and each grid point can be represented by two-dimensional coordinates; as shown in fig. 2, for a lattice point (i, j), there are a plurality of adjacent lattice points around it, this embodiment selects at least an upper lattice point (i, j +1), a lower lattice point (i, j-1), a right lattice point (i +1, j), and a left lattice point (i-1, j) that are adjacent to each other, and based on the actual situation, other lattice points that are adjacent to the lattice point (i, j) may also be selected, such as (i +1, j +1), (i-1, j-1), and the like, which this embodiment does not limit; in fig. 2, the i-axis represents the East direction (East) and the j-axis represents the North direction (North).

Step A2: and determining the latitudinal wind u (i, j +1) of the upper grid point, the latitudinal wind u (i, j-1) of the lower grid point, the latitudinal wind v (i +1, j) of the right grid point and the radial wind v (i-1, j) of the left grid point according to the wind field data.

Step A3: if u (i, j +1) <0, u (i, j-1) >0, v (i +1, j) >0 and v (i-1, j) <0 are simultaneously satisfied, and the latitudinal wind shear are both within a preset range [ w,1/w ], taking the lattice point (i, j) as an effective lattice point, and taking the relative vorticity of the lattice point (i, j) as the effective vorticity of the effective lattice point.

In the embodiment of the invention, the latitudinal wind u (i, j +1) of the upper lattice point is less than 0, which indicates that the latitudinal wind of the upper lattice point faces west (left); correspondingly, the weftwise wind u (i, j-1) >0 of the lower grid point indicates that the weftwise wind of the lower grid point is oriented east (rightward), the warp wind v (i +1, j) >0 of the right grid point indicates that the warp wind of the right grid point is oriented north (upward), the radial wind (i-1, j) <0 of the left grid point indicates that the warp wind of the left grid point is oriented south (downward). As shown in FIG. 2, if the four adjacent grid points satisfy u (i, j +1) <0, u (i, j-1) >0, v (i +1, j) >0, and v (i-1, j) <0 at the same time, it is described that the wind direction around the grid point (i, j) rotates counterclockwise basically, and has a certain quasi-symmetric vortex structure, and the grid point (i, j) may be the vortex center, that is, the grid point (i, j) can be used as an effective grid point.

Further, a preset range [ w,1/w ], where w is a preset parameter and w is between 0 and 1, e.g., w ═ 0.1, or w ═ 0.05, etc.; considering the complexity of the shape change of the vortex, the latitudinal wind shear and the latitudinal wind shear are set within [ w,1/w ], wherein the latitudinal wind shear is changed into a value corresponding to u (i, j +1)/u (i, j-1), and the latitudinal wind shear is changed into a value corresponding to v (i +1, j)/v (i-1, j).

Step 102: determining a target area with the effective grid points as the center, dividing the target area into N quadrants, and determining the wind parameters of each quadrant according to the wind field data; the wind parameters include an average wind direction, and the wind parameters further include at least one of an average latitudinal wind, an average latitudinal wind.

In the embodiment of the invention, the distance between a grid points and b grid points away from the effective grid point is taken as the identification radius R_fAnd identifying the radius R by using the effective lattice point as the center_fThe corresponding area range is called as a target area; for example, the target area is centered on the effective grid point and is a to b grid distance from the effective grid point center, i.e., the identification radius R_fB-a. For example, if the minimum scale plateau vortex diameter is 200km, the corresponding radius is 100km (1 °), and the ERA5 data resolution is 0.25 °, so b may be set to 4; a is less than b, and a can be 0, 1, 2 or 3. The target area with the effective lattice point as the center can be a square frame, b is half of the side length of the square outer frame, and a is half of the side length of the square inner frame; for example, as shown in fig. 3a, the hatched area represents the target area, where a is 2 and b is 4. Alternatively, the target area may be a rectangular frame, i.e. the target area is in the length direction and the width directionDifferent identification radii in direction; in the longitudinal direction, the radius R is identified_f1＝b₁-a₁In the width direction, the radius R is identified_f2＝b₂-a₂. See, for example, FIGS. 3b and 3c, wherein R_f2May be R_f1Half of that. The embodiment can identify vortexes in different shapes through target areas in different shapes.

In the embodiment of the present invention, one or more effective grid points may exist in the target area, and for each effective grid point, the area around the effective grid point is divided into a plurality of quadrants, for example, four quadrants (N-4) and eight quadrants (N-8) may be adopted. The "dividing the target area into N quadrants" may specifically include:

step B1: dividing the target area with the effective grid point as the center into 4 quadrants: a true east quadrant, a true north quadrant, a true west quadrant, and a true south quadrant.

In the embodiment of the present invention, fig. 4a shows a dividing manner for dividing the target area into four quadrants, where the middle point O is an effective grid point, and a square frame with the effective grid point as a center is selected to divide the four quadrants. FIG. 4a shows only one division of the square frame; in practical application, other rectangular frame dividing modes can also be adopted.

Step B2: dividing the target area with the effective grid point as the center into 8 quadrants: a first east quadrant, a second east quadrant, a first north quadrant, a second north quadrant, a first west quadrant, a second west quadrant, a first south quadrant, and a second south quadrant.

In the embodiment of the present invention, fig. 4b shows a dividing manner for dividing the target area into eight quadrants, the middle point O is an effective grid point, and a square frame with the effective grid point as a center is selected to divide four quadrants. FIG. 4b shows only one way of dividing the square frame; in practical application, other rectangular frame dividing modes can also be adopted.

In the embodiment of the present invention, after the target area is divided into N quadrants, the wind parameters of each quadrant need to be determined based on the wind field data. Specifically, an average wind direction needs to be determined for each quadrant, and in addition, at least one of an average latitudinal wind and an average latitudinal wind needs to be determined. In this embodiment, the average latitudinal wind of the quadrant refers to the average of latitudinal winds of all grid points of the quadrant; the average radial wind direction of the quadrant refers to the average radial wind direction of all grid points of the quadrant; the average wind direction of the quadrant refers to the average of the wind direction angles of all grid points of the quadrant.

Step 103: judging whether the average latitudinal wind or the average latitudinal wind of the quadrant meets the quadrant average wind field limiting condition; and forming the quadrant and other quadrants adjacent to the quadrant into a quadrant group, and judging whether the quadrant group meets the condition of anticlockwise rotation of the average wind direction of the quadrant group according to the average wind direction of the quadrant and the average wind direction of the other quadrants adjacent to the quadrant.

In an embodiment of the present invention, the determining whether the average latitudinal wind or the average latitudinal wind of the quadrant satisfies the quadrant average wind field limiting condition includes:

and under the condition that the target area is divided into four quadrants, namely a positive east quadrant, a positive north quadrant, a positive west quadrant and a positive south quadrant, judging whether the quadrant average wind field limiting condition is met according to the average latitudinal wind or the average latitudinal wind direction of the quadrants. If the average radial wind direction of the east quadrant is as shown in FIG. 4a

The average radial wind of the east-ward quadrant meets the quadrant average wind field limiting condition; if the average latitude of the north quadrant is wind

The average latitude wind of the north quadrant meets the quadrant average wind field limiting condition; if the average radial wind of the just west quadrant

The radial wind of the right west quadrant meets the quadrant mean wind field limiting condition; if the mean latitude of the normal south quadrant is wind

The radial wind of the positive south quadrant meets the quadrant mean wind field limiting condition.

In the embodiment of the present invention, if the target area is divided into eight quadrants, that is, under the condition that the target area is divided into eight quadrants, namely the first east quadrant, the second east quadrant, the first north quadrant, the second north quadrant, the first west quadrant, the second west quadrant, the first south quadrant and the second south quadrant, whether the four quadrants satisfy the quadrant mean wind field limiting condition is similar to the above-mentioned determination. If the average through-wind direction of the first east quadrant is as shown in FIG. 4b

The average radial wind direction of the first east quadrant meets the quadrant average wind field limiting condition; if the average through-wind direction of the second east quadrant

The average radial wind direction of the second east quadrant meets the quadrant average wind field limiting condition; if the average latitude of the first north quadrant is wind

The average weftwise wind of the first north quadrant meets quadrant average wind field limiting conditions; if the average latitude of the second north quadrant is wind

The average weftwise wind of the second north quadrant meets the quadrant average wind field limiting condition; if the average radial wind direction of the first west quadrant

The average radial wind of the first west quadrant meets the quadrant average wind field limiting condition; if the average radial wind direction of the second west quadrant

The average radial wind of the second west quadrant meets the quadrant average wind field limiting condition; if the average latitude of the first south quadrant is windward

The average weftwise wind of the first south quadrant meets quadrant average wind field limiting conditions; if the average latitude of the second south quadrant is wind

The mean weftwise wind of the second south quadrant satisfies the quadrant mean wind field defining condition.

Further optionally, the step of "determining whether the quadrant satisfies the condition of the mean counterclockwise rotation of the wind direction of the quadrant group" includes:

the average wind direction is defined as [ -180 °,180 ° ], and the directional deviation angle C (α, β) between two average wind directions α, β is defined as:

when the angle is more than 180 degrees and more than C (alpha, beta) and more than 0 degrees, the alpha trend beta is considered to be in the anticlockwise direction, and the condition of anticlockwise rotation of the average wind direction of the quadrant groups is met.

Specifically, in the case where the target area is divided into four quadrants of a true east quadrant, a true north quadrant, a true west quadrant, and a true south quadrant, a directional deviation angle of each quadrant group is determined

If the direction deviation angle of the quadrant group is between 0 degree and 180 degrees, determining that the quadrant group meets the condition of the anticlockwise rotation of the average wind direction of the quadrant group; wherein the content of the first and second substances,

represents the average wind direction of the righteast quadrant,

represents the average wind direction of the north quadrant,

represents the average wind direction of the west quadrant,

representing the average wind direction of the true south quadrant.

In the embodiment of the invention, in the four quadrants, if the direction deviation angle of the average wind direction of a certain quadrant and the average wind direction of the adjacent quadrant satisfies 180 degrees > C (alpha, beta) >0 degree, the quadrant group formed by the quadrant and the adjacent quadrant satisfies the condition of anticlockwise rotation of the average wind direction of the quadrant group. As shown in FIG. 4a, if

And the quadrant group consisting of the righteast quadrant and the rightnorth quadrant meets the condition of anticlockwise rotation of the average wind direction of the quadrant group.

Further, in a case where the target area is divided into eight quadrants, i.e., a first east quadrant, a second east quadrant, a first north quadrant, a second north quadrant, a first west quadrant, a second west quadrant, a first south quadrant, and a second south quadrant, a directional deviation angle of each quadrant group is determined

represents the average wind direction of the first east quadrant,

represents the average wind direction of the second east quadrant,

represents the average wind direction of the first north quadrant,

represents the average wind direction of the second north quadrant,

represents the average wind direction of the first west quadrant,

represents the average wind direction of the second west quadrant,

represents the average wind direction of the first south quadrant,

representing the average wind direction of the second south quadrant. The judgment principle of the eight quadrant is the same as that of the four quadrant, and the description thereof is omitted.

Step 104: and under the condition that m quadrants meet quadrant average wind field limiting conditions and n quadrant groups meet quadrant group average wind direction anticlockwise rotation conditions, taking the effective lattice points as undetermined lattice points.

In the embodiment of the invention, the target area can be divided into four quadrants, namely, the target area is divided into four quadrants including a positive east quadrant, a positive north quadrant, a positive west quadrant and a positive south quadrant; in this case, the step 104 includes:

step D1: if m in the four quadrants₄N satisfies the quadrant mean wind field limiting condition₄And if the quadrant groups meet the condition that the average wind direction of the quadrant groups rotates anticlockwise, the effective lattice point is considered to meet the condition that the four quadrants rotate anticlockwise.

In the embodiment of the invention, whether each quadrant meets the quadrant average wind field limiting condition or not is respectively determined. Specifically, see step 103 "determine whether the average latitudinal wind or the average latitudinal wind of the quadrant meets the quadrant average wind field limiting condition".

Furthermore, since there are four quadrants, each quadrant will be adjacent to two other quadrants; in this embodiment, it is also necessary to determine whether the average wind direction of the quadrant group formed by the quadrant and the adjacent other quadrants rotates counterclockwise. Wherein, the quadrant is two liang adjacent, can form four quadrant groups: the positive east quadrant + the positive north quadrant, the positive north quadrant + the positive west quadrant, the positive west quadrant + the positive south quadrant, and the positive south quadrant + the positive east quadrant.

In the embodiment of the invention, whether each quadrant group meets the condition of the mean wind direction anticlockwise rotation of the quadrant group is respectively determined.

If the quadrant group is east + north,

the quadrant group satisfies the counterclockwise rotation condition;

if the quadrant group is north plus west,

the quadrant group satisfies the counterclockwise rotation condition;

if the quadrant group is positive west quadrant + positive south quadrant,

the quadrant group satisfies the counterclockwise rotation condition;

if the quadrant group is the right south quadrant + the right east quadrant,

the quadrant group satisfies the counterclockwise rotation condition;

in the case that the target area is divided into four quadrants, since each quadrant and two adjacent quadrants thereof can form two quadrant groups, only one quadrant group can be selected in this embodiment; specifically, the quadrant group may be composed of the quadrant and a next quadrant adjacent to the counterclockwise direction (or the clockwise direction). For example, if a quadrant is a true east quadrant, the next quadrant adjacent to the quadrant in the counterclockwise direction is a true north quadrant, so that the quadrant group "true east quadrant + true north quadrant" can be rotated counterclockwise without paying attention to "true south quadrant + true east quadrant".

In the embodiment of the present invention, if m is satisfied₄Mean wind field of each quadrant defining a condition, and n₄And if the quadrant groups meet the condition that the average wind direction of the quadrant groups rotates anticlockwise, the four-quadrant anticlockwise rotation condition is met. Through tests, m₄＝4，n₄When the number of the quadrants is 4, the effective lattice point is taken as the lattice point to be determined only when all the quadrants and quadrant groups meet corresponding conditions; wherein, at this time, the valid grid point is considered to satisfy the four-quadrant counterclockwise rotation condition.

Step D2: dividing the target area into eight quadrants including a first east quadrant, a second east quadrant, a first north quadrant, a second north quadrant, a first west quadrant, a second west quadrant, a first south quadrant, and a second south quadrant; if m is present₈Each quadrant satisfies the quadrant mean wind field limiting condition, and n₈And if the quadrant groups meet the condition that the average wind direction of the quadrant groups rotates anticlockwise, determining that the effective lattice points meet the condition that the effective lattice points rotate anticlockwise of eight quadrants.

In the embodiment of the invention, whether each quadrant meets the quadrant average wind field limiting condition or not is respectively determined. In addition, for each quadrant, a condition whether it rotates counterclockwise from the average wind direction of the quadrant group is determined, respectively. Specifically, whether the quadrant group meets the quadrant group average wind direction counterclockwise rotation condition is determined based on the direction deviation angle.

If the quadrant group is the first east quadrant + the second east quadrant,

the quadrant group satisfies the counterclockwise rotation condition;

if the quadrant group is the second east quadrant + the first north quadrant,

the quadrant group satisfies the counterclockwise rotation condition;

if the quadrant group is the first north quadrant + the second north quadrant,

the quadrant group satisfies the counterclockwise rotation condition;

if the quadrant group is the second north quadrant + the first west quadrant,

the quadrant group satisfies the counterclockwise rotation condition;

if the quadrant group is the first west quadrant + the second west quadrant,

the quadrant group satisfies the counterclockwise rotation condition;

if the quadrant group is the second west quadrant + the first south quadrant,

the quadrant group satisfies the counterclockwise rotation condition;

if the quadrant group is the first south quadrant + the second south quadrant,

the quadrant group satisfies the counterclockwise rotation condition;

if the quadrant group is the second south quadrant + the first east quadrant,

the quadrant group satisfies the counterclockwise rotation condition;

in the examples of the present invention, if m₈Each quadrant satisfies a quadrant mean wind field limiting condition, and n₈And if the quadrant group meets the condition of the mean wind direction anticlockwise rotation of the quadrant group, the effective lattice point is considered to meet the condition of the eight-quadrant anticlockwise rotation. Through tests, m₈＝7；n₈When 7, can recognizeThe eight-quadrant counterclockwise rotation condition is satisfied for the valid grid points.

Step D3: and if the four-quadrant anticlockwise rotation condition is met and the eight-quadrant anticlockwise rotation condition is met, taking the effective lattice point as an undetermined lattice point.

In the embodiment of the invention, the area around the effective grid point is divided into four quadrants and eight quadrants, and if m is in the four quadrants, m is a positive integer₄Each quadrant satisfies a quadrant mean wind field limiting condition, and n₄If the quadrant groups meet the condition of the mean wind direction anticlockwise rotation of the quadrant groups, the four-quadrant anticlockwise rotation condition is met; if m in eight quadrants₈Each quadrant satisfies a quadrant mean wind field limiting condition, and n₈And if the quadrant group meets the condition of the mean wind direction anticlockwise rotation of the quadrant group, the eight-quadrant anticlockwise rotation condition is met. Through tests, m₄＝4，n₄＝4，m₈＝7，n₈7. And if the four-quadrant anticlockwise rotation condition is met and the eight-quadrant anticlockwise rotation condition is met, taking the effective lattice point as an undetermined lattice point.

Step 105: and determining the vortex center of the target vortex according to the undetermined point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined points.

In this embodiment of the present invention, the step 105 may specifically include:

step E1: a first distance threshold D between the two vortex centres is predetermined.

In the embodiment of the present invention, the first distance threshold D may be set based on practical experience. For example, if the minimum radius of the mesoscale vortex is 100km, the first distance threshold D may be set to 300 km.

Step E2: if the first to-be-checked point and the second to-be-checked point exist and the effective vorticity of the first to-be-checked point and the effective vorticity of the second to-be-checked point are both the maximum value or the maximum value in the preset range, the first to-be-checked point is taken as the center, and the side length is set to be

The corners of the square area face the positive direction, which is the east direction, the west direction and the south directionOr true north.

In the embodiment of the invention, the effective lattice point corresponding to the maximum value of the effective vorticity is generally used as the lattice point to be determined, and a plurality of maximum values possibly exist in the target area, so a plurality of lattice points to be determined possibly exist; and for the first to-be-checked point and the second to-be-checked point, determining whether the two to-be-checked points (namely the first to-be-checked point and the second to-be-checked point) belong to the same vortex by setting a square area. Wherein the square region can be seen in FIG. 5, the angle is oriented in the east direction, and the like, wherein O_ARepresenting a first to-be-ruled point with a coordinates of (x)_A,y_A)，O_BRepresenting a second to-be-ruled point with a lattice point coordinate of (x)_B,y_B)。

Step E3: if the second lattice point to be determined is in the square area, taking the position determined by the center between the first lattice point to be determined and the second lattice point to be determined as the vortex center of the target vortex; and if the second to-be-checked point is outside the square area, respectively taking the positions of the first to-be-checked point and the second to-be-checked point as the vortex centers of the corresponding target vortexes.

In the embodiment of the invention, whether the second undetermined point is located in the square area can be judged based on the longitude and latitude dis. The longitude and latitude refers to the sum of the warp lattice point difference and the weft lattice point difference of the two lattice points; as shown in FIG. 5, the latitude and longitude refers to O_AAnd O_BThe sum of the warp lattice point difference and the weft lattice point difference, i.e. the warp and weft distance dis ═ x_A-x_B|+|y_A-y_BL. Specifically, if the second to-be-qualified point is within the square area (that is, dis is equal to or less than D), it is indicated that the two to-be-qualified points are close in distance and most likely belong to the same vortex system, so that the position determined by the center between the first to-be-qualified point and the second to-be-qualified point is taken as the vortex center of the target vortex; specifically, if only two undetermined lattice points exist in the square area, the position of the center between the first undetermined lattice point and the second undetermined lattice point is directly used as the vortex center of the target vortex; if more undetermined lattice points exist in the square area, taking the average value of all undetermined lattice points as the vortex of the target vortexA center. Conversely, if the second to-be-checked point is outside the square region (i.e., dis)>D) The first to-be-qualified point and the second to-be-qualified point belong to two vortex systems, namely the positions of the first to-be-qualified point and the second to-be-qualified point can be respectively used as vortex centers of corresponding target vortexes.

According to the vortex identification method provided by the embodiment of the invention, the effective lattice point which is possibly the vortex center is selected based on the wind field data, and the vortex center is determined based on the effective vorticity under the condition of meeting the anticlockwise rotation condition. The method only selects part of grid points as effective grid points for identification, thereby reducing the processing amount; the quadrant average wind field limiting condition and the quadrant group average wind direction anticlockwise rotation condition are used as judging conditions, the 'averaging' idea not only embodies the physical meaning of anticlockwise rotation, but also is simpler in calculation, the flow is simplified, the vortex essence is grasped, and the vortex center can be determined very accurately.

On the basis of the above embodiment, after step 105 "determining the vortex center of the target vortex according to the to-be-checked point corresponding to the maximum value or maximum value of the effective vorticity in the to-be-checked points", the method further includes a process of determining the vortex radius, where the process specifically includes steps F1-F2:

step F1: and selecting k azimuth angles by taking the vortex center as an origin, and determining the outermost peripheral distance corresponding to the azimuth angles.

In the embodiment of the invention, the outermost peripheral distances in a plurality of azimuth angles are determined by taking the vortex center as the origin. For example, an azimuth angle is taken every 45 °, and a total of 8 azimuth angles can be determined; for the convenience of calculation, the 8 azimuth angles can be taken as the azimuth angles shown in fig. 6. Further, the outermost peripheral distance in each direction can be determined by utilizing the feature that the vortex center velocity (wind speed) is small and the outer peripheral velocity thereof increases as the radius increases. Optionally, the step F1 "determining the outermost distance corresponding to the azimuth" includes:

step F11: and determining a plurality of grid points in the direction of the azimuth angle, and determining the wind speed corresponding to each grid point.

Step F12: sequentially selecting a plurality of selected grid points from inside to outside every i grid points in the direction of the azimuth angle untilThe wind speed corresponding to the selected grid point is not increased any more, and the outermost grid point distance b corresponding to the selected grid point of the outermost periphery is determined_i，i∈[0,a]。

Step F13: and taking the average value of the distances between the outermost grid points determined by the intervals as the outermost distance corresponding to the azimuth angle.

In the embodiment of the invention, a plurality of selected lattice points are selected from inside to outside in the direction of the azimuth angle in a mode of spacing i lattice points, if the wind speed of the selected lattice points continuously increases, the selected lattice points are still positioned in the target vortex, otherwise, the selected lattice points are positioned outside the target vortex, and the lattice point distance on the azimuth angle, namely the outermost grid point distance, can be determined based on the positions of the selected lattice points. Furthermore, for each i, an outermost grid distance is determined. For example, if the azimuth is east and there are 10 grid points east of the vortex center, then the wind speed of 10 grid points can be determined, or only the wind speed of 8 grid points closer to the vortex center can be determined. Then, selecting every 0 grid point, namely selecting the 1 st grid point and the 2 nd grid point … … in sequence, and if the wind speed decreases from the 6 th grid point to the 7 th grid point, taking the distance between the 6 th grid point and the vortex center as the distance b between the outermost grid point and the vortex center₀. And then sequentially selecting every 1 grid point: the 1 st, 3 rd and 5 th grid points … …, and an outermost grid point distance b 1; and then sequentially selecting 2 grid points at intervals: the 1 st, 4 th, 7 th grid … … and so on until i equals a. After determining the plurality of outermost grid pitches, the average of the plurality of outermost grid pitches may be taken as the outermost distance in the azimuth (e.g., east).

Step F2: the average of the k outermost peripheral distances was taken as the vortex radius of the target vortex.

In the embodiment of the present invention, after determining the outermost distances at k azimuth angles, the average value of the outermost distances may be used as the vortex radius of the target vortex. Alternatively, the outermost peripheral distances in different directions may be corrected; for example, if 8 azimuth angles are selected as shown in 6, the final determined vortex radius may be:

wherein "east b" refers to the outermost distance from east to east, and so on, and will not be described in detail here. Because the vortex shape is complex, the average value of the outermost peripheral distances on a plurality of azimuth angles is used as the vortex radius, so that the size of the target vortex can be accurately represented; by determining the distances of the outermost grid points corresponding to different intervals i for multiple times, the outermost distances at the corresponding azimuth angles can be calculated more accurately.

On the basis of the above embodiment, after "determining the vortex center of the target vortex according to the to-be-gridded point corresponding to the maximum value or the maximum value of the effective vorticity in the to-be-gridded points" in step 105, the method further includes:

performing horizontal tracking on the target vortex: on a target vortex characteristic layer, matching the distance between the vortex center identified at the t-th time of the target vortex and the vortex center identified at the t + 1-th time, if the longitude and latitude distance between the vortex center and the vortex center does not exceed a second distance threshold D2, regarding the t-th time target vortex and the t + 1-th time target vortex as the same vortex, considering that the t-th time target vortex moves to the position of the t + 1-time target vortex, making the t-t +1, continuing horizontal tracking until the longitude and latitude distances of the two exceed a second distance threshold D2, considering that the vortex matched with the t + 1-time target vortex cannot be found at the t + 1-time, and ending the horizontal tracking of the vortex.

In the embodiment of the invention, the target vortex characteristic layer is an equal-pressure surface where the target vortex is mainly located; in general, the plateau vortex has a characteristic layer of 500hPa, the southwest vortex has a characteristic layer of 700hPa, and the middle-downstream vortex (Dabie mountain vortex) in the Yangtze river has a characteristic layer of 850 hPa. Determining the vortex center of the current time and the vortex center of the next time on the target vortex feature layer; for ERA5 data, a set of data is generated every 1 hour, i.e., every 1 hour for each epoch. If the longitude and latitude between the centers of the t-time vortices and the t + 1-time vortices does not exceed the second distance threshold, the t-time target vortex and the t + 1-time target vortex can be regarded as the same vortex. The true bookIn the embodiment, the longitude and latitude refers to the sum of the warp lattice point difference and the weft lattice point difference of two vortex centers, as shown in FIG. 5, if O_AAnd O_BThe two time-determined vortex centers are respectively, and the longitude and latitude dis between the two time-determined vortex centers is | x_A-x_B|+|y_A-y_B|。

For example, on the plateau vortex feature layer 500hPa, the vortex centers of the t-time plateau vortex and the t + 1-time plateau vortex are distance-matched, and if the longitude and latitude dis between the two times does not exceed the second distance threshold D2, the two times of plateau vortices are regarded as the same plateau vortex. In this embodiment, the second distance threshold is determined based on the moving speed of the target vortex, and horizontal tracking of the target vortex can be achieved by determining the longitude and latitude of the centers of the adjacent secondary vortices on the target vortex feature layer.

On the basis of the above embodiment, considering the influence of the underlying surface of the terrain, after "determining the vortex center of the target vortex according to the to-be-checked point corresponding to the maximum value or maximum value of the effective vorticity in the to-be-checked points" in step 105, sequentially performing layer-by-layer iterative identification on the continuous vertical layers of the target vortex feature layer upwards and downwards by taking the target vortex feature layer as a reference, that is, performing vertical tracking. The process of vertical tracking specifically includes:

tracking the target vortex vertically upwards: taking the target vortex feature layer as a kth layer, performing upward tracking identification in a rectangular frame with a search radius R, if the longitude and latitude between the vortex center identified by the kth-1 layer and the vortex center identified by the kth layer do not exceed a third distance threshold D3, keeping the upward tracking when k is k-1, and otherwise, interrupting the upward tracking;

tracking the target vortex vertically downwards: and taking the target vortex feature layer as a k-th layer, performing downward tracking identification in a rectangular frame with a search radius R, if the longitude and latitude between the vortex center identified by the k + 1-th layer and the vortex center identified by the k-th layer do not exceed a third distance threshold D3, making k equal to k +1, continuing downward tracking, and otherwise, interrupting the downward tracking.

In the embodiment of the invention, considering the calculation efficiency, only upward search identification and downward search layer identification are needed to be carried out in the area in the rectangular frame with the search radius R, and the method for identifying the vortex center is carried out according to the step 101 and the step 105, so that the vertical tracking of the vortex can be quickly realized.

In the embodiment of the invention, other equal-pressure surfaces are arranged above and below the target vortex characteristic layer, and at the current time, if the longitude and latitude between the vortex center on the equal-pressure surface and the vortex center on the target vortex characteristic layer do not exceed the third distance threshold D3, the target vortex also exists on the other equal-pressure surfaces, so that the vertical tracking of the target vortex can be realized. Referring to FIG. 7a, the vortex center C is present on the target vortex feature layer 500hPa_kIn the search range of radius R, the vortex center C exists on the upward isobaric surface 450hPa_k-1At this time, the center C of the vortex_kAnd C_k-1The longitude and latitude between do not exceed the third distance threshold D3, consider the vortex center tracking up to 450hPa, let k be k-1, continue with C_k-1And (4) centering, tracking upwards within a search range with the radius R, and stopping tracking upwards at the current time until the vortex center is not identified by the upper-layer isobaric surface or the identified vortex center exceeds a third distance threshold D3. Similarly, a vortex center C exists on the downward isobaric surface 550hPa_k+1Center of vortex C_kAnd vortex center C_k+1The longitude and latitude therebetween do not exceed a third distance threshold D3, consider the vortex center tracking down to 550hPa, let k be k +1, continue with C_k+1Centered, down-tracking within the search range of radius R, until the next layer contour does not identify a vortex center or the identified vortex center exceeds a third distance threshold D3, then stopping the current down-tracking.

In the embodiment of the invention, the influence of underlying surface terrain needs to be considered for downward vertical tracking of the target vortex, in order to save computing resources, the situation that different types of vortexes extend vertically can be considered for upward vertical tracking of the target vortex, the vertical tracking of the plateau vortex is taken as an example, the plateau vortex center identified by horizontal tracking 500hPa is taken as a reference layer, successive layer-by-layer iterative identification is sequentially carried out on continuous vertical layers of the plateau vortex center upwards and downwards at intervals of 50hPa, and the plateau vortex center is tracked to 200hPa upwards; taking into account the terrain effects, trace down to 550hPa, if the plateau vortex moves out of the plateau, trace down layer by layer to 950 hPa.

The method provided by the embodiment of the invention detects and verifies the identification of the vortex center based on ERA5 wind field data (0.25 degrees multiplied by 0.25 degrees and 1h), and the table shows that the hit rate, the missing report rate and the empty report rate of the vortex center identification under different sample time periods and different levels are randomly selected.

TABLE-Effect verification of vortex center identification based on ERA5 wind field

From the verification of the above table, the hit rate of the method provided by the embodiment of the invention on the plateau vortex, the southwest vortex and the Dabie mountain vortex is not much; because a high-altitude wind field is complex and has more disturbance, the high-altitude vortexes at the error rate (including the empty report rate and the missed report rate) are greater than the southwest vortexes and the big mountains vortexes, but compared with the traditional method based on the wind field, the method still greatly reduces the error rate, and improves the hit rate by about 4% (the hit rate of other current identification schemes is about 91%). In addition, the method has certain applicability to other vortex systems such as typhoons, medium beta scales and the like.

In addition, the identification method provided by the embodiment of the invention is used for identifying the 5-9 month plateau vortexes (500hPa) in 1979-2020, and the total number of times of sample: 107899 hours; the recognition result can be seen in fig. 7 b. Curve 1 in fig. 7b represents the hit rate, curve 2 represents the empty report rate, and curve 3 represents the missing report rate. The abscissa of fig. 7b is time, the left ordinate represents hit rate (%), and the right ordinate represents empty report rate and missing report rate (%). According to statistics, the hit rate is as follows: 91.25% -97.26%, average hit rate: 95.28 percent; the empty report rate is: 5.29% -12.59%, average empty report rate: 8.42 percent; the rate of missing reports is: 2.74% -8.75%, average rate of missing reports: 4.72 percent. Through large sample inspection, the identification method provided by the embodiment of the invention can achieve 95.28% of hit rate when identifying plateau vortexes, and has good identification effect. The vortex below the medium alpha scale is regarded as the empty report by the standard of manual judgment, and because the plateau wind field has more disturbance and the plateau vortex below the medium alpha scale has large proportion, the actual plateau vortex hit rate is higher and the empty report rate is lower; for example, the average empty report rate is now 8.42%, which may be around 4-5% in practice.

According to the vortex identification method provided by the embodiment of the invention, the effective lattice point which is possibly the vortex center is selected based on the wind field data, the four quadrants and the eight quadrants are divided, whether the anticlockwise circulation condition is met or not is judged according to the quadrant average wind direction limiting condition and the quadrant group average wind direction anticlockwise condition, the vortex center can be accurately identified, the vortex radius is determined by taking the wind speed around the vortex center as an index quantity, and the vortex is three-dimensionally tracked. Through inspection and evaluation, the method has strong universality, can be applied to identifying mesoscale vortexes in different shapes and different scales, and can also be applied to identifying the mesoscale vortexes based on wind field data with different spatial resolutions; particularly, the method can be applied to identification of mesoscale vortexes (particularly plateau vortexes) on a complex wind field, and can also have certain reference significance for identification of vortexes on the ocean.

The flow of the vortex identification method is described above in detail, and the method can also be implemented by a corresponding device, and the structure and function of the device are described in detail below.

Based on the same inventive concept, an embodiment of the present invention further provides a vortex identification device, as shown in fig. 8, the device includes:

the acquiring module 81 is configured to acquire wind field data in the whole area, select an effective lattice point from lattice points in which the relative vorticity is positive in the whole area, and use the relative vorticity of the effective lattice point as an effective vorticity;

a dividing module 82, configured to determine a target area centered on the effective grid point, divide the target area into N quadrants, and determine a wind parameter of each quadrant according to the wind field data; the wind parameters comprise an average wind direction, and the wind parameters further comprise at least one of an average latitudinal wind and an average latitudinal wind;

the judging module 83 is configured to judge whether the average latitudinal wind or the average latitudinal wind of the quadrant meets a quadrant average wind field limiting condition; forming the quadrant and other quadrants adjacent to the quadrant into a quadrant group, and judging whether the quadrant group meets the condition of anticlockwise rotation of the average wind direction of the quadrant group according to the average wind direction of the quadrant and the average wind direction of the other quadrants adjacent to the quadrant;

the processing module 84 is configured to, when m quadrants satisfy the quadrant average wind field limiting condition and n quadrant groups satisfy the quadrant group average wind direction counterclockwise rotation condition, take the valid lattice point as an undetermined lattice point; and determining the vortex center of the target vortex according to the undetermined point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined points.

On the basis of the foregoing embodiment, the obtaining module 81 selects an effective lattice point from the lattice points in which the relative vorticity in the entire region is positive, and takes the relative vorticity of the effective lattice point as the effective vorticity, including:

determining a plurality of grid points adjacent to a grid point (i, j) with positive relative vorticity, wherein the adjacent grid points at least comprise an upper grid point (i, j +1), a lower grid point (i, j-1), a right grid point (i +1, j) and a left grid point (i-1, j);

determining the latitudinal wind u (i, j +1) of the upper grid point, the latitudinal wind u (i, j-1) of the lower grid point, the latitudinal wind v (i +1, j) of the right grid point and the radial wind v (i-1, j) of the left grid point according to the wind field data;

if u (i, j +1) <0, u (i, j-1) >0, v (i +1, j) >0 and v (i-1, j) <0 are all satisfied, and the latitudinal wind shear are both within a preset range [ w,1/w ], taking the lattice point (i, j) as an effective lattice point, and taking the relative vorticity of the lattice point (i, j) as the effective vorticity of the effective lattice point.

On the basis of the above embodiment, the determining module 83

Judging whether the average latitudinal wind or the average latitudinal wind of the quadrant meets the quadrant average wind field limiting condition comprises the following steps:

in the case where the target area is divided into four quadrants, namely a true east quadrant, a true north quadrant, a true west quadrant, and a true south quadrant, if the average of the true east quadrant is facing the wind

The average radial wind direction of the east-ward quadrant meets the quadrant average wind field limiting condition; if the average latitude of the north quadrant is wind

The average weftwise wind of the due north quadrant meets the quadrant average wind field limiting condition; if the average radial wind of the just west quadrant

The radial wind of the west-west quadrant meets the quadrant mean wind field limiting condition; if the mean latitude of the normal south quadrant is wind

The radial wind of the south-positive quadrant meets the quadrant mean wind field limiting condition; and/or

In the case where the target area is divided into eight quadrants, i.e., a first east quadrant, a second east quadrant, a first north quadrant, a second north quadrant, a first west quadrant, a second west quadrant, a first south quadrant, and a second south quadrant, if the average of the first east quadrant is windward

The average weftwise wind of the first north quadrant meets the quadrant average wind field limiting condition; if the average latitude of the second north quadrant is wind

The average weftwise wind of the first south quadrant meets the quadrant average wind field limiting condition; if the average latitude of the second south quadrant is wind

The average weftwise wind of the second south quadrant satisfies the quadrant average wind field defining condition.

On the basis of the above embodiment, the determining module 83

Judging whether the quadrant group meets the condition of the mean wind direction anticlockwise rotation of the quadrant group comprises the following steps:

determining a directional deviation angle for each quadrant group with the target area divided into four quadrants, namely a true east quadrant, a true north quadrant, a true west quadrant, and a true south quadrant

represents the average wind direction of the righteast quadrant,

represents the average wind direction of the north quadrant,

represents the average wind direction of the west quadrant,

representing an average wind direction of the true south quadrant;

determining a directional deviation angle for each quadrant group when the target area is divided into eight quadrants, namely a first east quadrant, a second east quadrant, a first north quadrant, a second north quadrant, a first west quadrant, a second west quadrant, a first south quadrant, and a second south quadrant

represents the average wind direction of the first east quadrant,

represents the average wind direction of the second east quadrant,

represents the average wind direction of the first north quadrant,

represents the average wind direction of the second north quadrant,

represents the average wind direction of the first west quadrant,

represents the average wind direction of the second west quadrant,

represents the average wind direction of the first south quadrant,

representing the average wind direction of the second south quadrant.

On the basis of the above embodiment, the processing module 84

Taking the effective lattice point as a lattice point to be determined under the condition that m quadrants meet the quadrant mean wind field limiting condition and n quadrant groups meet the quadrant group mean wind direction anticlockwise rotation condition, wherein the lattice point to be determined comprises:

dividing the target area into four quadrants including a true east quadrant, a true north quadrant, a true west quadrant, and a true south quadrant; if the plurality of quadrants meet the quadrant average wind field limiting condition and the plurality of quadrant groups meet the quadrant group average wind direction anticlockwise rotation condition at the moment, determining that the effective lattice point meets the four-quadrant anticlockwise rotation condition;

dividing the target area into eight quadrants including a first east quadrant, a second east quadrant, a first north quadrant, a second north quadrant, a first west quadrant, a second west quadrant, a first south quadrant, and a second south quadrant; if the plurality of quadrants meet the quadrant average wind field limiting condition and the plurality of quadrant groups meet the quadrant group average wind direction anticlockwise rotation condition at the moment, determining that the effective lattice point meets the eight-quadrant anticlockwise rotation condition;

and under the condition that the four-quadrant anticlockwise rotation condition and the eight-quadrant anticlockwise rotation condition are simultaneously met, taking the effective lattice point as a lattice point to be determined.

On the basis of the foregoing embodiment, the determining, by the processing module 84, the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points includes:

predetermining a first distance threshold D between two vortex centers;

if a first to-be-checked point and a second to-be-checked point exist, and the effective vorticity of the first to-be-checked point and the effective vorticity of the second to-be-checked point are both the maximum value or the maximum value within a preset range, taking the first to-be-checked point as the center, and setting the side length as

The angle of the square area faces to the positive direction, and the positive direction is the positive east direction, the positive west direction, the positive south direction or the positive north direction;

if the second point to be qualified is in the square area, taking the position determined by the center between the first point to be qualified and the second point to be qualified as the vortex center of the target vortex; and if the second to-be-fixed point is outside the square area, respectively taking the positions of the first to-be-fixed point and the second to-be-fixed point as the vortex centers of the corresponding target vortexes.

On the basis of the above embodiment, the apparatus further includes: a radius determination module;

after the processing module 84 determines the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points, the radius determining module is configured to:

selecting k azimuth angles by taking the vortex center as an origin, and determining the outermost peripheral distance corresponding to the azimuth angles;

taking the average of k outermost peripheral distances as the vortex radius of the target vortex.

On the basis of the foregoing embodiment, the determining, by the radius determining module, the outermost peripheral distance corresponding to the azimuth angle includes:

determining a plurality of grid points in the direction of the azimuth angle, and determining the wind speed corresponding to each grid point;

sequentially selecting a plurality of selected grid points from inside to outside at intervals of i grid points in the direction of the azimuth angle until the wind speed corresponding to the selected grid points is not increased any more, and determining the distance b between the outermost grid points corresponding to the selected outermost grid points_i，i∈[0,a]；

And taking the average value of the distances of the outermost grid points determined by the intervals as the outermost distance corresponding to the azimuth angle.

On the basis of the above embodiment, the system further comprises a horizontal tracking module;

after the processing module 84 determines the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points, the horizontal tracking module is configured to:

performing horizontal tracking on the target vortex: and on the target vortex characteristic layer, performing distance matching on the vortex center identified at the t th time of the target vortex and the vortex center identified at the t +1 th time of the target vortex, if the longitude and latitude distance between the target vortex and the vortex center is not more than a second distance threshold D2, taking the time t of the target vortex and the time t +1 of the target vortex as the same vortex, and making the time t equal to the time t +1, continuing horizontal tracking until the longitude and latitude distances between the target vortex and the target vortex exceed a second distance threshold D2, and ending the horizontal tracking of the vortex.

On the basis of the above embodiment, the system further comprises a vertical tracking module;

after the processing module 84 determines the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points, the vertical tracking module is configured to:

Embodiments of the present invention also provide a computer storage medium storing computer-executable instructions including a program for performing the above-described method for identifying a vortex, the computer-executable instructions being capable of performing the method in any of the above-described method embodiments.

The computer storage media may be any available media or data storage device that can be accessed by a computer, including but not limited to magnetic memory (e.g., floppy disks, hard disks, magnetic tape, magneto-optical disks (MOs), etc.), optical memory (e.g., CDs, DVDs, BDs, HVDs, etc.), and semiconductor memory (e.g., ROMs, EPROMs, EEPROMs, non-volatile memory (NAND FLASH), Solid State Disks (SSDs)), etc.

Fig. 9 shows a block diagram of an electronic device according to another embodiment of the present invention. The electronic device 1100 may be a host server with computing capabilities, a personal computer PC, or a portable computer or terminal that is portable, or the like. The specific embodiment of the present invention does not limit the specific implementation of the electronic device.

The electronic device 1100 includes at least one processor (processor)1110, a Communications Interface 1120, a memory 1130, and a bus 1140. The processor 1110, the communication interface 1120, and the memory 1130 communicate with each other via the bus 1140.

The communication interface 1120 is used for communicating with network elements including, for example, virtual machine management centers, shared storage, etc.

Processor 1110 is configured to execute programs. Processor 1110 may be a central processing unit CPU, or an Application Specific Integrated Circuit (ASIC), or one or more Integrated circuits configured to implement embodiments of the present invention.

The memory 1130 is used for executable instructions. The memory 1130 may comprise high-speed RAM memory, and may also include non-volatile memory (non-volatile memory), such as at least one disk memory. The memory 1130 may also be a memory array. The storage 1130 may also be partitioned and the blocks may be combined into virtual volumes according to certain rules. The instructions stored by the memory 1130 are executable by the processor 1110 to enable the processor 1110 to perform the method of vortex identification in any of the method embodiments described above.

The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present invention, and all the changes or substitutions should be covered within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims

1. A method of identifying a vortex, comprising:

determining the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points;

determining the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points comprises the following steps:

predetermining a first distance threshold D between two vortex centers;

2. The method according to claim 1, wherein the selecting effective grid points from the grid points with positive relative vorticity in the whole area and taking the relative vorticity of the effective grid points as the effective vorticity comprises:

3. The method of claim 1, wherein determining whether the quadrant mean latitudinal or mean latitudinal wind satisfies a quadrant mean wind field defining condition comprises:

The meridian direction of the just west quadrant meets the lawThe quadrant mean wind field limiting condition; if the mean latitude of the normal south quadrant is wind

The radial wind of the south-positive quadrant meets the quadrant mean wind field limiting condition;

and/or

4. The method of claim 1, wherein determining whether the quadrant group meets a quadrant group average wind direction counterclockwise rotation condition comprises:

represents the average wind direction of the righteast quadrant,

represents the average wind direction of the north quadrant,

represents the average wind direction of the west quadrant,

representing an average wind direction of the true south quadrant;

represents the average wind direction of the first east quadrant,

represents the average wind direction of the second east quadrant,

represents the firstThe average wind direction in the north quadrant,

represents the average wind direction of the second north quadrant,

represents the average wind direction of the first west quadrant,

represents the average wind direction of the second west quadrant,

represents the average wind direction of the first south quadrant,

representing the average wind direction of the second south quadrant.

5. The method according to claim 1, wherein the regarding the effective lattice point as a lattice point in the case that m quadrants satisfy the quadrant mean wind field definition condition and n quadrant groups satisfy the quadrant group mean wind direction counterclockwise rotation condition comprises:

6. The method according to any one of claims 1 to 5, wherein after determining the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points, the method further comprises:

7. The method of claim 6, wherein the determining the outermost peripheral distance corresponding to the azimuth angle comprises:

8. The method according to any one of claims 1 to 5, wherein after determining the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points, the method further comprises:

9. The method according to any one of claims 1 to 5, wherein after determining the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points, the method further comprises:

10. An apparatus for identifying vortices, comprising:

the processing module is used for taking the effective lattice point as a lattice point to be determined under the condition that m quadrants meet the quadrant mean wind field limiting condition and n quadrant groups meet the quadrant group mean wind direction anticlockwise rotation condition; determining the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points;

the processing module determines the vortex center of the target vortex according to the undetermined lattice point corresponding to the maximum value or the maximum value of the effective vorticity in the undetermined lattice points, and the method comprises the following steps:

predetermining a first distance threshold D between two vortex centers;

11. A computer storage medium having stored thereon computer-executable instructions for performing the method of vortex identification of any of claims 1-9.

12. An electronic device, comprising:

at least one processor; and the number of the first and second groups,

the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of vortex identification of any of claims 1-9.