patents.google.com

CN118962844A - A cyclone center positioning method and device - Google Patents

️Fri Nov 15 2024

Detailed Description

The following description of the embodiments of the present application will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are some, but not all embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.

The determination of the location of the cyclone center requires the combination of various observations and analysis means, which is the key content of cyclone monitoring and forecasting.

In one example, the current weather data gives the time resolution of the path information of the cyclone center of the pacific northwest, which is 10 hours, and the change of the position of the cyclone center every 10 hours is recorded. In this way, the location of the cyclone center can be determined every 10 hours time interval based on current weather data.

However, although the time resolution of 10 hours can primarily capture the main dynamic characteristics of the position of the cyclone center, data support can be primarily provided for monitoring, forecasting and analyzing the cyclone, but research and analysis for lower time resolution are not satisfied, and in addition, the mode of determining the position of the cyclone center through the current meteorological data is complex, and the consumed computational resources are much.

Accordingly, a need has arisen to quickly and simply determine the location of the cyclone center.

Further, the present application has been made for the purpose of achieving a quick and simple determination of the position of the cyclone center.

For example, the application provides a cyclone center positioning method, which can quickly and simply determine the position of the cyclone center.

Before the technical scheme of the application is introduced, technical terms possibly related to the application are explained.

Cyclone (Cyclone), a meteorological term, refers to a large vortex in the north (south) hemisphere, where the horizontal air flow in the atmosphere is counter (clockwise) rotating. The cyclone is formed by a wave motion occurring on the frontal surface or on the air interface with different densities, and occupies a large scale horizontal air vortex of three dimensions in the atmosphere. The central air pressure is lower than the surrounding air pressure, and is also called "low air pressure" in terms of the air pressure field.

Cyclone can be classified into temperate cyclone and tropical cyclone according to generation and movement ranges.

The temperate zone cyclone is commonly provided with a frontal cyclone, and the tropical cyclone is commonly provided with typhoons and tropical storms. The air flow is divided into a northern hemisphere cyclone and a southern hemisphere cyclone according to hemispherical distribution, wherein the air flow rotates anticlockwise, and the air flow rotates clockwise. Temperate cyclones occur mainly in eastern asia, north america, the Mediterranean, etc. The tropical cyclone presents regional distribution characteristics that the northern hemisphere is far more than the southern hemisphere and the tropical cyclone intensively occurs between 10 degrees to 20 degrees of latitude.

Tropical Cyclone (Tropical Cyclone), a collective term for non-frontal weather-scale vortices generated on Tropical or subtropical ocean surfaces, with organized convection and defined cyclonic circulation. It includes tropical low voltage, tropical storm, strong hot zone storm, typhoon, strong typhoon and super strong typhoon. Tropical cyclones are common in the western and its immediate seas (typhoons), the atlantic and northeast pacific (hurricanes), the indian and southern pacific.

Specifically, referring to fig. 1, a flowchart illustrating steps of a cyclone center positioning method according to the present application is shown, where the method may be applied to an electronic device, and the electronic device may include a terminal or a server, and the method may specifically include the following steps:

In step S101, a wind speed in a first horizontal direction of each of a plurality of subregions in a region in which the cyclone is located is acquired, and a wind speed in a second horizontal direction of each of a plurality of subregions in a region in which the cyclone is located is acquired, the first horizontal direction being perpendicular to the second horizontal direction.

In the present application, the wind speed of each of the plurality of subregions in the region in which the cyclone is located in the first horizontal direction and the wind speed of each of the plurality of subregions in the region in which the cyclone is located in the second horizontal direction may be directly obtained from the existing meteorological data.

For example, existing meteorological data may be obtained according to a simulation of a coupling mode of an offshore wind energy resource, for example, the coupling mode of the offshore wind energy resource may be a coupling mode between atmosphere and ocean waves, specifically, a coupling mode between a meteorological mesoscale model and a marine power model and a wave model, and the above-mentioned coupling mode may output simulated data, for example, simulated meteorological mesoscale data, marine power data, wave data, and the like, where the data includes wind speeds of each of a plurality of subregions in an area where a cyclone is located in a first horizontal direction and wind speeds of each of a plurality of subregions in an area where the cyclone is located in a second horizontal direction.

For another example, existing weather data may be obtained directly from the aforementioned cyclone path data from the northwest pacific of the China weather office (CMA-STI), the northwest pacific of the Japan typhoon center (JMA RSMC Tokyo), and/or the northwest pacific of the United states typhoon alert center (JTWC).

Wherein, for north west pacific, the cyclone path data of north west pacific includes but is not limited to 3 kinds: the cyclone path data for the North Pacific ocean from the China weather office (CMA-STI), the North Pacific ocean from the Tokyo center (JMA RSMC Tokyo), and the North Pacific ocean from the United states typhoon alert center (JTWC).

All 3 of these documents record path information including the cyclone center of the pacific northwest since 1950 s. The path information of the cyclone center of the North Pacific ocean can be directly obtained by any one of the 3 data, so that the position of the cyclone center at any moment can be obtained.

However, the time resolution of the path information including the cyclone center of the pacific northwest among these 3 materials all employed a lower hour time resolution, and the change in the position of the cyclone center at a lower time interval was recorded. Thus, the position of the lower time interval of the cyclone center can be determined from any of these 3 materials.

The area where the cyclone is located is a large area, for example, an area of several hundred kilometers by several hundred kilometers, etc. The shape of the region where the cyclone is located is variable, and various shapes are possible in practice, including regular shapes or irregular shapes, etc.

The shape of the plurality of subregions in the region where the cyclone is located may be rectangular, hexagonal, triangular, circular, square, or the like. There may be no overlap between the multiple sub-regions.

The first horizontal direction may be a north-south direction, for example, a north-south direction or a north-south direction, and the second horizontal direction may be an east-west direction. For example, the number of the cells to be processed, in the east-west direction or from the west to the east, etc.

For example, referring to fig. 2, fig. 2 shows a schematic view of a plurality of sub-areas in the zone where the cyclone is located.

In fig. 2, there are shown 25 sub-zones included in the zone in which the cyclone is located, and outside the 25 sub-zones, there may be more sub-zones included in the zone in which the cyclone is located, which are not shown one by one in fig. 2.

Fig. 2 is an illustration of the shape of each sub-region as a square, but is not intended to limit the scope of the present application.

The area of a sub-area may be 15 km by 15 km, or 20 km by 20 km, or 25 km by 25 km, or 30 km by 30 km, etc. The schematic diagram of fig. 2 shows the multiple sub-areas in the area where the cyclone is located from a top view.

In step S102, for any one sub-region, the overall wind speed change rate of the plurality of adjacent sub-regions of the sub-region is determined according to the wind speeds of the plurality of adjacent sub-regions of the sub-region in the first horizontal direction and/or the wind speeds of the plurality of adjacent sub-regions of the sub-region in the second horizontal direction.

The same is true for each of the other sub-regions, thereby respectively deriving the overall wind speed rate of the plurality of adjacent sub-regions of the respective sub-region.

That is, each sub-region corresponds to one overall wind speed change rate, and the overall wind speed change rate corresponding to any one sub-region is the overall wind speed change rate of the sub-region for a plurality of adjacent sub-regions.

Wherein, for any one sub-region, the sub-region has a plurality of adjacent sub-regions, the sub-region is adjacent to the adjacent sub-region of the sub-region, and other sub-regions are not separated.

For example, referring to fig. 2, for sub-region e, the number of adjacent sub-regions of sub-region e may be 8, respectively: sub-region a, sub-region b, sub-region c, sub-region d, sub-region f, sub-region g, sub-region h, and sub-region i.

Each adjacent sub-area of the sub-area has a wind speed, respectively, and the wind speed has a direction.

The overall wind speed change rate of a plurality of adjacent subareas of the subarea may be represented: the magnitude of the change in wind speed (in relation to the wind direction) of a plurality of adjacent sub-areas of the sub-area.

The same is true for each of the other sub-regions.

The present step may be specifically referred to the embodiments shown later, and will not be described in detail herein.

In step S103, among a plurality of sub-areas in the area where the cyclone is located, a sub-area where the overall wind speed change rate is the largest is selected.

In the application, each subarea can be ordered according to the sequence from high to low of the overall wind speed change rate of a plurality of adjacent subareas of each subarea, namely, each subarea is ordered according to the sequence from high to low of the overall wind speed change rate corresponding to each subarea, and then the subarea with the first ordering is selected, namely, the subarea with the largest overall wind speed change rate.

In step S104, the position of the cyclone center is determined according to the sub-region where the overall wind speed change rate is the greatest.

In the present application, the sub-region where the overall wind speed variation rate is largest may be determined as the region where the cyclone center is located, and in one embodiment, the position of the center of the sub-region where the overall wind speed variation rate is largest may be determined as the position of the cyclone center, or the position in the vicinity of the center of the sub-region where the overall wind speed variation rate is largest may be determined as the position of the cyclone center, or the like.

In the present application, the wind speed of each of the plurality of subregions in the region in which the cyclone is located in a first horizontal direction is acquired, and the wind speed of each of the plurality of subregions in the region in which the cyclone is located in a second horizontal direction is acquired, the first horizontal direction being perpendicular to the second horizontal direction. For any one sub-region, determining the overall wind speed change rate of a plurality of adjacent sub-regions of the sub-region according to the wind speeds of the plurality of adjacent sub-regions of the sub-region in the first horizontal direction and/or the wind speeds of the plurality of adjacent sub-regions of the sub-region in the second horizontal direction. Among a plurality of subregions in the region in which the cyclone is located, a subregion in which the overall wind speed change rate is the greatest is selected. And determining the position of the cyclone center according to the subarea with the maximum integral wind speed change rate.

Wherein the overall wind speed change rate of the adjacent sub-areas of one sub-area can represent the change amplitude of the wind speed (related to the wind direction) of the adjacent sub-areas of the one sub-area. Adjacent sub-regions of different sub-regions are different.

Through statistical analysis of actual data of a large number of cyclone wind speed data (related to wind direction), it was found that the overall wind speed change rate of a plurality of adjacent subareas of subareas where the center of the cyclone is located is greater than the overall wind speed change rate of a plurality of adjacent subareas of subareas where the center of the non-cyclone is located, that is, the wind speed (related to wind direction) change amplitude of a plurality of adjacent subareas of subareas where the center of the cyclone is located is greater than the wind speed (related to wind direction) change amplitude of a plurality of adjacent subareas of subareas where the center of the non-cyclone is located.

Therefore, in the multiple subareas in the area where the cyclone is located, the subarea with the largest overall wind speed change rate can be regarded as the subarea where the cyclone center is located, so that the position of the cyclone center can be determined according to the subarea where the cyclone center is located, and the positioning of the cyclone center is realized.

The wind speed of each subarea in the area where the cyclone is located in the first horizontal direction and the wind speed of each subarea in the area where the cyclone is located in the second horizontal direction are lower time resolution data, for example, the resolution of 1 hour time interval and even the resolution of less than 1 hour time interval, so that the time resolution of positioning the cyclone center can be reduced, for example, the position of the cyclone center in every 1 hour time interval, that is, the position of the cyclone center in every hour can be determined, and even the time resolution of the position of the determined cyclone center in the time interval can be improved.

In conclusion, the scheme for positioning the cyclone center has the advantages of simple operation process, less consumed calculation resources and high speed, and can improve the efficiency of determining the position of the cyclone center.

In one embodiment of the present application, referring to fig. 3, step S102 includes:

In step S201, a wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction is determined according to wind speeds of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction.

In one embodiment of the present application, the present step may be implemented by the following procedure, including:

2011. Among a plurality of adjacent sub-areas of the sub-area, a plurality of target sub-areas of which the wind direction is related to the first horizontal direction are determined.

In this step, a sub-region between the wind direction and the first horizontal direction of less than 90 degrees may be determined, and a sub-region between the wind direction and the opposite direction of the first horizontal direction of less than 90 degrees may be determined, among a plurality of adjacent sub-regions of the sub-region. And determining the determined subarea as a target subarea.

For example, based on the example shown in fig. 2, referring to fig. 4, it is assumed that the first horizontal direction is a west-to-east direction and the second horizontal direction is a south-to-north direction. And assuming that the sub-region is sub-region e, the number of adjacent sub-regions of sub-region e is 8, and the number is respectively: sub-region a, sub-region b, sub-region c, sub-region d, sub-region f, sub-region g, sub-region h, and sub-region i.

In the schematic diagram shown in fig. 4, arrows are shown in each of the sub-areas a, b, c, d, f, g, h and i, and the direction of the arrow in one sub-area refers to the wind direction (direction of wind) of the sub-area.

Thus, sub-areas smaller than 90 degrees between the wind direction and the west-to-east direction have sub-areas g, h, and i, and sub-areas smaller than 90 degrees between the wind direction and the opposite direction of the west-to-east direction have sub-areas a, b, and c.

Therefore, sub-region g, sub-region h, sub-region i, sub-region a, sub-region b, and sub-region c are target sub-regions.

2012. Dividing the target subareas into a first subarea set and a second subarea set, wherein the difference between the wind direction of any target subarea in the first subarea set and the wind direction of any target subarea in the second subarea set is more than 90 degrees.

For example, target sub-areas with less than 90 degrees between the wind direction and the first horizontal direction may be combined as a first set of sub-areas, and target sub-areas with less than 90 degrees between the wind direction and the opposite direction of the first horizontal direction may be combined as a second set of sub-areas.

For example, in the above example, sub-regions g, h, and i are combined into a first set of sub-regions, and sub-regions a, b, and c are combined into a second set of sub-regions.

The wind direction of the subarea g is different from the wind direction of the subarea a, the wind direction of the subarea b and the wind direction of the subarea c by more than 90 degrees respectively.

The wind direction of the subarea h is different from the wind direction of the subarea a, the wind direction of the subarea b and the wind direction of the subarea c by more than 90 degrees respectively.

The wind direction of the subarea i is different from the wind direction of the subarea a, the wind direction of the subarea b and the wind direction of the subarea c by more than 90 degrees respectively.

2013. And carrying out weighted operation on the wind speed of the target subarea in the first subarea set and the wind speed of the target subarea in the second subarea set to obtain a weighted wind speed in the first horizontal direction.

In one embodiment of the present application, the present step may be implemented by the following procedure, including:

11 And carrying out weighted operation on the wind speeds of all the target subareas in the first subarea set to obtain the weighted wind speed of the first subarea set.

In this step, the weights of the target sub-areas in the first set of sub-areas may be determined according to the distances between the sub-areas and the target sub-areas in the first set of sub-areas. The weight of the target subarea with the larger distance from the subarea in the first subarea set is smaller, and the weight of the target subarea with the smaller distance from the subarea in the first subarea set is larger.

And, the number of target subregions in the first set of subregions may be counted.

And then, according to the weights of all the target subareas in the first subarea set and the number of the target subareas in the first subarea set, carrying out weighting operation on the wind speeds of all the target subareas in the first subarea set to obtain the weighted wind speeds of the first subarea set.

For example, the wind speeds of the target subareas in the first subarea set may be weighted and summed according to the weights of the target subareas in the first subarea set to obtain a wind speed sum value, the weights of the target subareas in the first subarea set are summed to obtain a weight sum value, a product between the wind speed sum value and the weight sum value is calculated, and a ratio between the product and the coefficient is calculated to obtain a weighted wind speed of the first subarea set.

Wherein, for the coefficient, it can be calculated by:

for example, for any one target sub-region in the first set of sub-regions, if the wind speed of the target sub-region can be obtained, the sub-coefficient of the target sub-region is equal to the product between the weight of the target sub-region and a specific value, which is greater than 0, the specific value may be 1, 1.5 or 2, etc., or if the wind speed of the target sub-region cannot be obtained, the sub-coefficient of the target sub-region is equal to the product between the weight of the target sub-region and the value 0, i.e., the sub-coefficient of the target sub-region is equal to 0.

The same is true for each other target sub-region in the first set of sub-regions.

The coefficients may then be summed up for the sub-coefficients of the respective target sub-regions in the first set of sub-regions.

12 And (3) carrying out weighted operation on the wind speeds of all the target subareas in the second subarea set to obtain the weighted wind speeds of the second subarea set.

This step is specifically referred to as step 11), and is not described in detail herein.

13 And performing weighted calculation on the weighted wind speeds of the first subarea set and the weighted wind speeds of the second subarea set to obtain the weighted wind speed in the first horizontal direction.

In one embodiment of the present application, a sum value between the weighted wind speeds of the first sub-area set and the weighted wind speeds of the second sub-area set may be calculated, then a ratio between the sum value and the number of the plurality of target sub-areas determined in step 2011 may be calculated, and then the ratio may be used as the weighted wind speed in the first horizontal direction, or a ratio between the ratio and the length or width of the sub-area may be calculated and used as the weighted wind speed in the first horizontal direction.

2014. According to the weighted wind speed in the first horizontal direction, the wind speed change rate of a plurality of adjacent subareas of the subareas in the first horizontal direction is determined.

In the present application, the weighted wind speed in the first horizontal direction may be taken as the wind speed change rate in the first horizontal direction of the plurality of adjacent sub-regions of the sub-region.

In step S202, the overall wind speed change rate of the plurality of adjacent sub-areas of the sub-area is determined according to the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction.

In one embodiment, the wind speed change rate of the plurality of adjacent sub-regions of the sub-region in the first horizontal direction may be determined as the overall wind speed change rate of the plurality of adjacent sub-regions of the sub-region, or a product between the wind speed change rate of the plurality of adjacent sub-regions of the sub-region in the first horizontal direction and 100% may be calculated, the product may be determined as the overall wind speed change rate of the plurality of adjacent sub-regions of the sub-region, or the like.

In another embodiment of the present application, referring to fig. 5, step S102 includes:

In step S301, a rate of change of wind speed of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction is determined according to wind speeds of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction.

In one embodiment of the present application, the present step may be implemented by the following procedure, including:

3011. Among a plurality of adjacent sub-areas of the sub-area, a plurality of target sub-areas of the wind direction associated with the second horizontal direction are determined.

In this step, a sub-region between the wind direction and the second horizontal direction of less than 90 degrees may be determined, and a sub-region between the wind direction and the opposite direction of the second horizontal direction of less than 90 degrees may be determined, among a plurality of adjacent sub-regions of the sub-region. And determining the determined subarea as a target subarea.

For example, based on the example shown in fig. 2, referring to fig. 4, it is assumed that the first horizontal direction is a west-to-east direction and the second horizontal direction is a south-to-north direction. And assuming that the sub-region is sub-region e, the number of adjacent sub-regions of sub-region e is 8, and the number is respectively: in the schematic diagram shown in fig. 4, arrows are indicated in all of the sub-areas a, b, c, d, f, g, h and i, and the direction of the arrow in the sub-areas refers to the wind direction (direction of wind) of the sub-area.

Thus, sub-areas smaller than 90 degrees between the wind direction and the north-south direction have sub-areas c, f and i, and sub-areas smaller than 90 degrees between the wind direction and the opposite direction of the west-east direction have sub-areas a, d and g.

Thus, sub-region a, sub-region d, sub-region g, sub-region c, sub-region f, and sub-region g are target sub-regions.

3012. Dividing the plurality of target subareas into a third subarea set and a fourth subarea set, wherein the difference between the wind direction of any one target subarea in the third subarea set and the wind direction of any one target subarea in the fourth subarea set is more than 90 degrees.

For example, target sub-areas with less than 90 degrees between the wind direction and the second horizontal direction may be combined as a third sub-area set, and target sub-areas with less than 90 degrees between the wind direction and the opposite direction of the second horizontal direction may be combined as a fourth sub-area set.

For example, in the above example, sub-region c, sub-region f, and sub-region i are combined into a third sub-region set, and sub-region a, sub-region d, and sub-region g are combined into a fourth sub-region set.

The wind direction of the subarea c is different from the wind direction of the subarea a, the wind direction of the subarea d and the wind direction of the subarea g by more than 90 degrees respectively.

The wind direction of the subarea f is different from the wind direction of the subarea a, the wind direction of the subarea d and the wind direction of the subarea g by more than 90 degrees respectively.

The wind direction of the subarea i is different from the wind direction of the subarea a, the wind direction of the subarea d and the wind direction of the subarea g by more than 90 degrees respectively.

3013. And carrying out weighted operation on the wind speed of the target subarea in the third subarea set and the wind speed of the target subarea in the fourth subarea set to obtain a weighted wind speed in the second horizontal direction.

In one embodiment of the present application, the present step may be implemented by the following procedure, including:

21 And (3) carrying out weighted operation on the wind speeds of all the target subareas in the third subarea set to obtain the weighted wind speed of the third subarea set.

In this step, the weights of the respective target subregions in the third set of subregions may be determined according to the distance between the subregion and the respective target subregions in the third set of subregions. The weight of the target subarea with the larger distance from the subarea in the third subarea set is smaller, and the weight of the target subarea with the smaller distance from the subarea in the third subarea set is larger.

And, the number of target subregions in the third set of subregions may be counted.

And then, according to the weights of all the target subareas in the third subarea set and the number of the target subareas in the third subarea set, carrying out weighting operation on the wind speeds of all the target subareas in the third subarea set to obtain the weighted wind speeds of the third subarea set.

For example, the wind speeds of the target subareas in the third subarea set may be weighted and summed according to the weights of the target subareas in the third subarea set to obtain a wind speed sum value, the weights of the target subareas in the third subarea set may be summed to obtain a weight sum value, a product between the wind speed sum value and the weight sum value may be calculated, and a ratio between the product and the coefficient may be calculated to obtain a weighted wind speed of the third subarea set.

Wherein, for the coefficient, it can be calculated by:

For example, for any one of the target sub-regions in the third set of sub-regions, if the wind speed of the target sub-region can be obtained, the sub-coefficient of the target sub-region is equal to the product between the weight of the target sub-region and a specific value, which is greater than 0, the specific value may be 1, 1.5 or 2, etc., or if the wind speed of the target sub-region cannot be obtained, the sub-coefficient of the target sub-region is equal to the product between the weight of the target sub-region and the value 0, i.e., the sub-coefficient of the target sub-region is equal to 0.

The same is true for each other target sub-region in the third set of sub-regions.

The coefficients may then be summed up for the sub-coefficients of the respective target sub-regions in the third set of sub-regions.

22 And (3) carrying out weighted operation on the wind speeds of all the target subareas in the fourth subarea set to obtain the weighted wind speed of the fourth subarea set.

This step is specifically referred to as step 21), and is not described in detail herein.

23 And acquiring the weighted wind speed in the second horizontal direction according to the weighted wind speed of the third subarea set and the weighted wind speed of the fourth subarea set.

In one embodiment of the present application, a sum value between the weighted wind speeds of the third sub-area set and the weighted wind speeds of the fourth sub-area set may be calculated, then a ratio between the sum value and the number of the plurality of target sub-areas determined in step 3011 may be calculated, and then the ratio may be used as the weighted wind speed in the first horizontal direction, or a ratio between the ratio and the length or width of the sub-area may be calculated and used as the weighted wind speed in the first horizontal direction.

3014. And determining the wind speed change rate of a plurality of adjacent subareas of the subareas in the second horizontal direction according to the weighted wind speed in the second horizontal direction.

In step S302, an overall wind speed change rate of the plurality of adjacent sub-areas of the sub-area is determined according to the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction.

In another embodiment of the present application, referring to fig. 6, step S102 includes:

In step S401, a rate of change of wind speed in the first horizontal direction of the plurality of adjacent sub-areas of the sub-area is determined according to wind speeds in the first horizontal direction of the plurality of adjacent sub-areas of the sub-area.

This step is specifically referred to as step S201, and will not be described in detail herein.

In step S402, a rate of change of wind speed in the second horizontal direction of the plurality of adjacent sub-areas of the sub-area is determined according to wind speeds in the second horizontal direction of the plurality of adjacent sub-areas of the sub-area.

This step is specifically referred to as step S301, and will not be described in detail herein.

The execution order of step S401 and step S402 is not limited in the present application.

Step S401 may be performed first and then step S402 may be performed, or step S402 may be performed first and then step S401 may be performed.

In step S403, the overall wind speed change rate of the plurality of adjacent sub-areas of the sub-area is determined according to the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction and the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction.

In one embodiment of the present application, a first square of a wind speed change rate of a plurality of adjacent sub-areas of the sub-area in a first horizontal direction may be calculated, and a second square of a wind speed change rate of a plurality of adjacent sub-areas of the sub-area in a second horizontal direction may be calculated, then a sum value between the first square and the second square may be calculated, then an evolution of the sum value may be calculated, then an overall wind speed change rate of a plurality of adjacent sub-areas of the sub-area may be determined from the evolution, for example, the evolution may be determined directly as an overall wind speed change rate of a plurality of adjacent sub-areas of the sub-area, or a product between the evolution and 100% may be calculated, the product may be determined as an overall wind speed change rate of a plurality of adjacent sub-areas of the sub-area, or the like.

The present application is illustrated by one specific embodiment and is not intended to be limiting in scope.

(1) Downloading ERA5 the "ERA5 hourly data on SINGLE LEVELS from 1940 to present" dataset in the re-analysis data. The method comprises the steps of selecting and downloading offshore plane u wind speed (east-west wind speed) 10m u-component of wind and v wind speed (north-south wind speed) 10m v-component of wind, selecting the year, month, day and hour when tropical cyclone occurs, selecting the longitude and latitude range of the tropical cyclone passing-through influence range, and selecting NetCDF format of downloading.

(2) In ArcGIS software, the NetCDF format u wind speed and v wind speed components are converted into TIFF grid format, respectively.

(3) The gradients are calculated for the u-wind speed grid and the v-wind speed grid, respectively.

The grade calculation is expressed in percent (percent elevation delta).

Wherein the percentage of elevation increment of the grade (elevation) is equal to the elevation increment divided by the horizontal increment and multiplied by 100 (for example, when the grade angle is 45 degrees, the elevation increment is equal to the horizontal increment, so the percentage of elevation increment is 100 percent; when the grade angle is close to the right angle of 90 degrees, the percentage of elevation increment starts to be close to infinity).

The gradient is calculated by adopting a plane method, and the gradient is measured according to the maximum change rate from one pixel to the pixel directly adjacent to the gradient. The calculation is performed on the projection plane using a 2D cartesian coordinate system. The grade value is calculated by means of the maximum average method (Burrough, 1998). Plane and geodesic calculations are performed using a neighborhood of 3*3 pels (moving window). For each neighborhood, if the processing (center) pel is NoData (lack of valid data), then the output is NoData (lack of valid data) as well. The calculation also requires that at least 7 pixels adjacent to the processed pixel have a significant value. If there are fewer than 7 valid pels, no computation is performed and the output at this processed pel will be NoData (lack of valid data). The tool may calculate a maximum rate of change of value in the direction from that pel to the pel adjacent to it, the percent slope (delta) depending on the rate of change of the surface in the horizontal (dz/dx) and vertical (dz/dy) directions from the center pel:

As shown in fig. 2 or 4, the values of the center pixel e and its neighboring eight pixels determine the horizontal increment and the vertical increment. These adjacent pels are represented using letters a through i, where pel e represents the pel that is currently calculating the slope.

The rate of change of picture element e in the x-direction is equal to:

[dz/dx] =| ((c + 2f + i)*4/wght1|+|(a + 2d + g)*4/wght2)| / (8 * x_cellsize)

in the above formulas wght and wght are the horizontal weighted counts of the active pixels. When c, f, and i each have a significant value, then wght = (1+2×1+1) =4; when i is NoData, then wght 1= (1+2×1+0) =3; when f is NoData, wght 1= (1+2×0+1) =2. Similar logic applies to wght, except that the neighborhood positions are a, d, and g. The absolute value symbol is represented by the absolute value. x cellsize represents the size of the picture element. 8 represents that the center pixel e has 8 neighboring pixels.

The rate of change of picture element e in the y direction is equal to:

[dz/dy] = ((g + 2h + i)*4/wght3 - (a + 2b + c)*4/wght4) / (8 * y_cellsize)

In the above equation wght and wght are the same concept as in the [ dz/dx ] calculation.

(4) The u wind speed gradient percentage and the v wind speed gradient (elevation increment) percentage are added in a grid corresponding mode so as to more highlight pixel values of the tropical cyclone center relative to the peripheral area.

(5) And (3) obtaining the pixel with the largest grid median value in the step (4) as the tropical cyclone center position in the hour.

Through the spatial analysis of the sea level u wind speed and v wind speed grid data of the analysis data, the position of the tropical cyclone center in each hour can be simply, rapidly and accurately obtained, the position of the tropical cyclone center in each adjacent 6 hours in the existing tropical cyclone optimal path data product can be filled, and the tropical cyclone early warning and forecasting and disaster assessment level can be improved.

It should be noted that, for simplicity of explanation, the method embodiments are shown as a series of acts, but it should be understood by those skilled in the art that the present application is not limited by the order of acts, as some steps may occur in other orders or concurrently in accordance with the application. Further, those skilled in the art will appreciate that the embodiments described in the specification are all alternative embodiments and that the actions involved are not necessarily required for the present application.

Referring to fig. 7, there is shown a block diagram of a cyclone centering device of the present application, comprising:

An obtaining module 11, configured to obtain wind speeds of each of a plurality of subareas in a region where the cyclone is located in a first horizontal direction, and obtain wind speeds of each of a plurality of subareas in a region where the cyclone is located in a second horizontal direction, where the first horizontal direction is perpendicular to the second horizontal direction;

the first determining module 12 is configured to determine, for any one sub-area, an overall wind speed change rate of a plurality of adjacent sub-areas of the sub-area according to wind speeds of the plurality of adjacent sub-areas of the sub-area in a first horizontal direction and/or wind speeds of the plurality of adjacent sub-areas of the sub-area in a second horizontal direction;

a selecting module 13, configured to select a sub-area with the greatest overall wind speed change rate from a plurality of sub-areas in the area where the cyclone is located;

The second determining module 14 is configured to determine a location of the cyclone center according to the sub-area with the greatest overall wind speed change rate.

In an alternative implementation, the first determining module includes:

A first determining unit, configured to determine a wind speed change rate of a plurality of adjacent subareas of the subareas in a first horizontal direction according to wind speeds of the plurality of adjacent subareas of the subareas in the first horizontal direction;

And the second determining unit is used for determining the overall wind speed change rate of the plurality of adjacent subareas of the subareas according to the wind speed change rate of the plurality of adjacent subareas of the subareas in the first horizontal direction.

In an alternative implementation, the first determining module includes:

a third determining unit, configured to determine a wind speed change rate of the plurality of adjacent subareas of the subareas in the second horizontal direction according to wind speeds of the plurality of adjacent subareas of the subareas in the second horizontal direction;

And the fourth determining unit is used for determining the overall wind speed change rate of the plurality of adjacent subareas of the subareas according to the wind speed change rate of the plurality of adjacent subareas of the subareas in the second horizontal direction.

In an alternative implementation, the first determining module includes:

And a fifth determining unit, configured to determine an overall wind speed change rate of the plurality of adjacent sub-areas of the sub-area according to a wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction and a wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction.

In an alternative implementation, the first determining unit includes:

a first determining subunit configured to determine, among a plurality of adjacent sub-areas of the sub-areas, a plurality of target sub-areas in which a wind direction is related to a first horizontal direction;

The first dividing subunit is used for dividing the plurality of target subareas into a first subarea set and a second subarea set, and the difference between the wind direction of any one target subarea in the first subarea set and the wind direction of any one target subarea in the second subarea set is more than 90 degrees;

the first weighting operation subunit is used for carrying out weighting operation on the wind speed of the target subarea in the first subarea set and the wind speed of the target subarea in the second subarea set to obtain a weighted wind speed in the first horizontal direction;

And the second determination subunit is used for determining the wind speed change rate of a plurality of adjacent subareas of the subareas in the first horizontal direction according to the weighted wind speed in the first horizontal direction.

In an alternative implementation, the first determining subunit is specifically configured to: determining a sub-region of less than 90 degrees between the wind direction and the first horizontal direction, and determining a sub-region of less than 90 degrees between the wind direction and the opposite direction of the first horizontal direction, from among a plurality of adjacent sub-regions of the sub-region; and determining the determined subarea as a target subarea.

In an alternative implementation, the first weighting operator unit is specifically configured to: performing weighted operation on the wind speeds of all target subareas in the first subarea set to obtain weighted wind speeds of the first subarea set; performing weighted operation on the wind speeds of all target subareas in the second subarea set to obtain weighted wind speeds of the second subarea set; and carrying out weighted operation on the weighted wind speeds of the first subarea set and the weighted wind speeds of the second subarea set to obtain the weighted wind speed in the first horizontal direction.

In an alternative implementation, the first weighting operator unit is specifically configured to: determining the weight of each target subarea in the first subarea set according to the distance between the subarea and each target subarea in the first subarea set; counting the number of target subareas in the first subarea set; and carrying out weighting operation on the wind speeds of all the target subareas in the first subarea set according to the weights of all the target subareas in the first subarea set and the number of the target subareas in the first subarea set to obtain the weighted wind speeds of the first subarea set.

In an alternative implementation, the third determining unit includes:

A third determining subunit configured to determine, among a plurality of adjacent sub-areas of the sub-areas, a plurality of target sub-areas in which a wind direction is related to a second horizontal direction;

The second dividing subunit is used for dividing the plurality of target subareas into a third subarea set and a fourth subarea set, and the difference between the wind direction of any one target subarea in the third subarea set and the wind direction of any one target subarea in the fourth subarea set is more than 90 degrees;

the second weighting operation subunit is used for carrying out weighting operation on the wind speed of the target subarea in the third subarea set and the wind speed of the target subarea in the fourth subarea set to obtain a weighted wind speed in the second horizontal direction;

And a fourth determination subunit, configured to determine a wind speed change rate of a plurality of adjacent sub-areas of the sub-areas in the second horizontal direction according to the weighted wind speeds in the second horizontal direction.

In an alternative implementation, the third determining subunit is specifically configured to: determining a sub-region of less than 90 degrees between the wind direction and the second horizontal direction, and determining a sub-region of less than 90 degrees between the wind direction and the opposite direction of the second horizontal direction, from among a plurality of adjacent sub-regions of the sub-region; and determining the determined subarea as a target subarea.

In an alternative implementation, the second weighting operator unit is specifically configured to: performing weighted operation on the wind speeds of all target subareas in the third subarea set to obtain weighted wind speeds of the third subarea set; performing weighted operation on the wind speeds of all target subareas in the fourth subarea set to obtain weighted wind speeds of the fourth subarea set; and carrying out weighted operation on the weighted wind speeds of the third subarea set and the fourth subarea set to obtain the weighted wind speed in the second horizontal direction.

In an alternative implementation, the second weighting operator unit is specifically configured to: determining the weight of each target subarea in the third subarea set according to the distance between the subarea and each target subarea in the third subarea set; counting the number of target subareas in the third subarea set; and carrying out weighting operation on the wind speeds of all the target subareas in the third subarea set according to the weights of all the target subareas in the third subarea set and the number of the target subareas in the third subarea set to obtain the weighted wind speeds of the third subarea set.

In an alternative implementation, the fifth determining unit includes:

a first calculating subunit for calculating a first square of a wind speed change rate of a plurality of adjacent sub-areas of the sub-areas in a first horizontal direction;

A second calculation subunit for calculating a second square of a wind speed change rate of a plurality of adjacent sub-areas of the sub-areas in a second horizontal direction;

a third calculation subunit for calculating a sum value between the first square and the second square;

a fourth calculation subunit for calculating an evolution of the sum;

And a fifth determination subunit, configured to determine an overall wind speed change rate of a plurality of adjacent sub-areas of the sub-area according to the evolution.

For the device embodiments, since they are substantially similar to the method embodiments, the description is relatively simple, and reference is made to the description of the method embodiments for relevant points.

Optionally, an embodiment of the present application further provides an electronic device, including: the processor, the memory, the computer program stored in the memory and capable of running on the processor, the computer program realizes each process of the above method embodiment when being executed by the processor, and can achieve the same technical effect, and for avoiding repetition, the description is omitted here.

The embodiment of the application also provides a computer readable storage medium, on which a computer program is stored, which when executed by a processor, realizes the processes of the above method embodiment and can achieve the same technical effects, and in order to avoid repetition, the description is omitted here. The computer readable storage medium is, for example, a Read-Only Memory (ROM), a random access Memory (Random Access Memory RAM), a magnetic disk or an optical disk.

Fig. 8 is a block diagram of an electronic device 800 in accordance with the present application. For example, electronic device 800 may be a mobile phone, computer, digital broadcast terminal, messaging device, game console, tablet device, medical device, exercise device, personal digital assistant, or the like.

Referring to fig. 8, an electronic device 800 may include one or more of the following components: a processing component 802, a memory 804, a power component 806, a multimedia component 808, an audio component 810, an input/output (I/O) interface 812, a sensor component 814, and a communication component 816.

The processing component 802 generally controls overall operation of the electronic device 800, such as operations associated with display, telephone calls, data communications, camera operations, and recording operations. The processing component 802 may include one or more processors 820 to execute instructions to perform all or part of the steps of the methods described above. Further, the processing component 802 can include one or more modules that facilitate interactions between the processing component 802 and other components. For example, the processing component 802 can include a multimedia module to facilitate interaction between the multimedia component 808 and the processing component 802.

The memory 804 is configured to store various types of data to support operations at the electronic device 800. Examples of such data include instructions for any application or method operating on the electronic device 800, contact data, phonebook data, messages, images, videos, and so forth. The memory 804 may be implemented by any type or combination of volatile or nonvolatile memory devices such as Static Random Access Memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic or optical disk.

The power supply component 806 provides power to the various components of the electronic device 800. The power components 806 may include a power management system, one or more power sources, and other components associated with generating, managing, and distributing power for the electronic device 800.

The multimedia component 808 includes a screen between the electronic device 800 and the user that provides an output interface. In some embodiments, the screen may include a Liquid Crystal Display (LCD) and a Touch Panel (TP). If the screen includes a touch panel, the screen may be implemented as a touch screen to receive input signals from a user. The touch panel includes one or more touch sensors to sense touches, swipes, and gestures on the touch panel. The touch sensor may sense not only the boundary of a touch or slide action, but also the duration and pressure associated with the touch or slide operation. In some embodiments, the multimedia component 808 includes a front camera and/or a rear camera. When the electronic device 800 is in an operational mode, such as a shooting mode or a video mode, the front camera and/or the rear camera may receive external multimedia data. Each front camera and rear camera may be a fixed optical lens system or have focal length and optical zoom capabilities.

The audio component 810 is configured to output and/or input audio signals. For example, the audio component 810 includes a Microphone (MIC) configured to receive external audio signals when the electronic device 800 is in an operational mode, such as a call mode, a recording mode, and a voice recognition mode. The received audio signals may be further stored in the memory 804 or transmitted via the communication component 816. In some embodiments, audio component 810 further includes a speaker for outputting audio signals.

The I/O interface 812 provides an interface between the processing component 802 and peripheral interface modules, which may be a keyboard, click wheel, buttons, etc. These buttons may include, but are not limited to: homepage button, volume button, start button, and lock button.

The sensor assembly 814 includes one or more sensors for providing status assessment of various aspects of the electronic device 800. For example, the sensor assembly 814 may detect an on/off state of the electronic device 800, a relative positioning of the components, such as a display and keypad of the electronic device 800, the sensor assembly 814 may also detect a change in position of the electronic device 800 or a component of the electronic device 800, the presence or absence of a user's contact with the electronic device 800, an orientation or acceleration/deceleration of the electronic device 800, and a change in temperature of the electronic device 800. The sensor assembly 814 may include a proximity sensor configured to detect the presence of nearby objects without any physical contact. The sensor assembly 814 may also include a light sensor, such as a CMOS or CCD image sensor, for use in imaging applications. In some embodiments, the sensor assembly 814 may also include an acceleration sensor, a gyroscopic sensor, a magnetic sensor, a pressure sensor, or a temperature sensor.

The communication component 816 is configured to facilitate communication between the electronic device 800 and other devices, either wired or wireless. The electronic device 800 may access a wireless network based on a communication standard, such as WiFi, an operator network (e.g., 2G, 3G, 4G, or 5G), or a combination thereof. In one exemplary embodiment, the communication component 816 receives broadcast signals or broadcast operation information from an external broadcast management system via a broadcast channel. In one exemplary embodiment, the communication component 816 further includes a Near Field Communication (NFC) module to facilitate short range communications. For example, the NFC module may be implemented based on Radio Frequency Identification (RFID) technology, infrared data association (IrDA) technology, ultra Wideband (UWB) technology, bluetooth (BT) technology, and other technologies.

In an exemplary embodiment, the electronic device 800 may be implemented by one or more Application Specific Integrated Circuits (ASICs), digital Signal Processors (DSPs), digital Signal Processing Devices (DSPDs), programmable Logic Devices (PLDs), field Programmable Gate Arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic elements for executing the methods described above.

In an exemplary embodiment, a non-transitory computer readable storage medium is also provided, such as memory 804 including instructions executable by processor 820 of electronic device 800 to perform the above-described method. For example, the non-transitory computer readable storage medium may be ROM, random Access Memory (RAM), CD-ROM, magnetic tape, floppy disk, optical data storage device, etc.

Fig. 9 is a block diagram of an electronic device 1900 in accordance with the present application. For example, electronic device 1900 may be provided as a server.

Referring to FIG. 9, electronic device 1900 includes a processing component 1922 that further includes one or more processors and memory resources represented by memory 1932 for storing instructions, such as application programs, that can be executed by processing component 1922. The application programs stored in memory 1932 may include one or more modules each corresponding to a set of instructions. Further, processing component 1922 is configured to execute instructions to perform the methods described above.

The electronic device 1900 may also include a power component 1926 configured to perform power management of the electronic device 1900, a wired or wireless network interface 1950 configured to connect the electronic device 1900 to a network, and an input/output (I/O) interface 1958. The electronic device 1900 may operate based on an operating system stored in memory 1932, such as Windows Server, mac OS XTM, unixTM, linuxTM, freeBSDTM, or the like.

It should be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.

From the above description of the embodiments, it will be clear to those skilled in the art that the above-described embodiment method may be implemented by means of software plus a necessary general hardware platform, but of course may also be implemented by means of hardware, but in many cases the former is a preferred embodiment. Based on such understanding, the technical solution of the present application may be embodied essentially or in a part contributing to the prior art in the form of a software product stored in a storage medium (e.g. ROM/RAM, magnetic disk, optical disk) comprising instructions for causing a terminal (which may be a mobile phone, a computer, a server, an air conditioner, or a network device, etc.) to perform the method according to the embodiments of the present application.

The embodiments of the present application have been described above with reference to the accompanying drawings, but the present application is not limited to the above-described embodiments, which are merely illustrative and not restrictive, and many forms may be made by those having ordinary skill in the art without departing from the spirit of the present application and the scope of the claims, which are to be protected by the present application.

Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.

It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described systems, apparatuses and units may refer to corresponding procedures in the foregoing method embodiments, and are not repeated herein.

In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other manners. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of the units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.

The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

In addition, each functional unit in the embodiments of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit.

The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a usb disk, a removable hard disk, a ROM, a RAM, a magnetic disk, or an optical disk, etc.

The foregoing is merely illustrative of the present application, and the present application is not limited thereto, and any person skilled in the art will readily recognize that variations or substitutions are within the scope of the present application. Therefore, the protection scope of the application is subject to the protection scope of the claims.