A cyclone center positioning method and device

Abstract

The present application provides a method and device for locating the center of a cyclone. The wind speed of each sub-region in a plurality of sub-regions in the region where the cyclone is located in the first horizontal direction is obtained, and the wind speed of each sub-region in a plurality of sub-regions in the region where the cyclone is located in the second horizontal direction is obtained, and the first horizontal direction is perpendicular to the second horizontal direction. For any 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 speed of the plurality of adjacent sub-regions of the sub-region in the first horizontal direction and / or the wind speed of the plurality of adjacent sub-regions of the sub-region in the second horizontal direction. Among the plurality of sub-regions in the region where the cyclone is located, the sub-region with the largest overall wind speed change rate is selected. The position of the center of the cyclone is determined according to the sub-region with the largest overall wind speed change rate. The calculation process of the scheme for locating the center of the cyclone of the present application is simple, consumes less computing resources, and is fast, which can improve the efficiency of locating the center of the cyclone.

Description

Technical Field

[0001] The present application relates to the field of resource utilization technology, and in particular to a method and device for locating the center of a cyclone. Background Art

[0002] The center of a cyclone refers to the area with the lowest air pressure and the highest wind speed in the cyclone system, that is, the center of the cyclone eye. The position change of the cyclone center directly determines the movement direction and intensity change of the cyclone. The position of the cyclone center is of great significance in predicting the movement trajectory of the cyclone, estimating the intensity of the cyclone, assessing the risk of typhoon disasters, guiding emergency response, and correcting and improving numerical weather forecast models.

[0003] Therefore, a need arises to determine the location of the cyclone center. Summary of the invention

[0004] The present application shows a method and device for locating the center of a cyclone.

[0005] In a first aspect, the present application provides a method for locating a cyclone center, the method comprising:

[0006] Obtaining the wind speed of each of the multiple sub-regions in the region where the cyclone is located in a first horizontal direction, and obtaining the wind speed of each of the multiple sub-regions in the region where the cyclone is located in a second horizontal direction, the first horizontal direction being perpendicular to the second horizontal direction;

[0007] For any sub-region, determining the overall wind speed change rate of the multiple adjacent sub-regions of the sub-region according to the wind speed of the multiple adjacent sub-regions of the sub-region in the first horizontal direction and / or the wind speed of the multiple adjacent sub-regions of the sub-region in the second horizontal direction;

[0008] Among the multiple sub-areas in the area where the cyclone is located, select the sub-area with the largest overall wind speed change rate;

[0009] The location of the cyclone center is determined based on the sub-area with the largest overall wind speed change rate.

[0010] In a second aspect, the present application shows a cyclone center positioning device, the device comprising:

[0011] an acquisition module, configured to acquire the wind speed of each of the multiple sub-regions in the region where the cyclone is located in a first horizontal direction, and to acquire the wind speed of each of the multiple sub-regions in the region where the cyclone is located in a second horizontal direction, wherein the first horizontal direction is perpendicular to the second horizontal direction;

[0012] A first determination module is used to determine, for any sub-area, the overall wind speed change rate of multiple adjacent sub-areas of the sub-area according to the wind speed of multiple adjacent sub-areas of the sub-area in the first horizontal direction and / or the wind speed of multiple adjacent sub-areas of the sub-area in the second horizontal direction;

[0013] A selection module is used to select a sub-area with the largest overall wind speed change rate from among multiple sub-areas in the area where the cyclone is located;

[0014] The second determination module is used to determine the position of the cyclone center according to the sub-area with the largest overall wind speed change rate.

[0015] In a third aspect, the present application shows an electronic device, which includes: a processor; a memory for storing processor-executable instructions; wherein the processor is configured to execute the method described in any of the above aspects.

[0016] In a fourth aspect, the present application illustrates a non-temporary computer-readable storage medium, which, when instructions in the storage medium are executed by a processor of an electronic device, enables the electronic device to execute a method as described in any of the above aspects.

[0017] In a fifth aspect, the present application illustrates a computer program product. When instructions in the computer program product are executed by a processor of an electronic device, the electronic device is enabled to execute the method described in any of the above aspects.

[0018] The technical solution provided by this application may have the following beneficial effects:

[0019] In the present application, the wind speed of each sub-region in the first horizontal direction of multiple sub-regions in the region where the cyclone is located is obtained, and the wind speed of each sub-region in the second horizontal direction of multiple sub-regions in the region where the cyclone is located is obtained, and the first horizontal direction is perpendicular to the second horizontal direction. For any sub-region, the overall wind speed change rate of multiple adjacent sub-regions of the sub-region is determined based on the wind speed of multiple adjacent sub-regions of the sub-region in the first horizontal direction and / or the wind speed of multiple adjacent sub-regions of the sub-region in the second horizontal direction. Among the multiple sub-regions in the region where the cyclone is located, the sub-region with the largest overall wind speed change rate is selected. The position of the cyclone center is determined based on the sub-region with the largest overall wind speed change rate.

[0020] The overall wind speed change rate of multiple adjacent sub-regions of a sub-region can reflect the change range of wind speed (involving wind direction) of multiple adjacent sub-regions of the sub-region. The adjacent sub-regions of different sub-regions are different.

[0021] After statistical analysis of a large number of actual cyclone wind speed data (involving wind direction), it was found that the overall wind speed change rate of multiple adjacent sub-regions of the sub-region where the center of the cyclone is located is greater than the overall wind speed change rate of multiple adjacent sub-regions of the sub-region where the center of the cyclone is located, that is, the change amplitude of the wind speed (involving wind direction) of multiple adjacent sub-regions of the sub-region where the center of the cyclone is located is greater than the change amplitude of the wind speed (involving wind direction) of multiple adjacent sub-regions of the sub-region where the center of the non-cyclone is located.

[0022] Therefore, among the multiple sub-areas in the area where the cyclone is located, the sub-area with the largest overall wind speed change rate can be regarded as the sub-area where the center of the cyclone is located. Therefore, the position of the center of the cyclone can be determined based on the sub-area where the center of the cyclone is located, so as to locate the center of the cyclone.

[0023] Among them, the wind speed in the first horizontal direction of each of the multiple sub-regions in the area where the cyclone is located in the existing meteorological data and the wind speed in the second horizontal direction of each of the multiple sub-regions in the area where the cyclone is located are data with lower time resolution, for example, the resolution of a 1-hour time interval, or even a resolution less than a 1-hour time interval. In this way, through the present application, the time resolution of locating the center of the cyclone can be reduced. For example, the position of the center of the cyclone can be determined within each 1-hour time interval, that is, the position of the center of the cyclone can be determined hour by hour, even less than the position within a 1-hour time interval, which can improve the time resolution of the determined position of the center of the cyclone.

[0024] In summary, the solution for locating the center of a cyclone in the present application has a simple calculation process, consumes less computing resources, and has a high speed, which can improve the efficiency of determining the position of the center of a cyclone. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] Figure 1 It is a flowchart of the steps of a cyclone center positioning method of the present application.

[0026] Figure 2 It is a schematic diagram of the area where a cyclone is located in the present application.

[0027] Figure 3 It is a flow chart of the steps of a method for determining the overall wind speed change rate of the present application.

[0028] Figure 4 It is a schematic diagram of the area where a cyclone is located in the present application.

[0029] Figure 5 It is a flow chart of the steps of a method for determining the overall wind speed change rate of the present application.

[0030] Figure 6It is a flow chart of the steps of a method for determining the overall wind speed change rate of the present application.

[0031] Figure 7 This is a structural block diagram of a cyclone center positioning device of the present application.

[0032] Figure 8 It is a block diagram of an electronic device of the present application.

[0033] Fig. 9 It is a block diagram of an electronic device of the present application. DETAILED DESCRIPTION

[0034] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0035] Determining the location of a cyclone center requires combining a variety of observation data and analysis methods, and is a key part of cyclone monitoring and forecasting.

[0036] In one example, the time resolution of the path information of the cyclone center in the northwest Pacific Ocean given by the current meteorological data is 10 hours, and the change of the position of the cyclone center at every 10-hour time interval is recorded. In this way, the position of the cyclone center at every 10-hour time interval can be determined based on the current meteorological data.

[0037] However, although the 10-hour time resolution can initially capture the main dynamic characteristics of the location of the cyclone center and can initially provide data support for cyclone monitoring, forecasting and analysis, it does not meet the research and analysis requirements for lower time resolution. In addition, the method of determining the location of the cyclone center through current meteorological data is complex and consumes a lot of computing resources.

[0038] Therefore, there is a need to quickly and easily determine the location of the cyclone center.

[0039] Furthermore, in order to achieve the purpose of quickly and simply determining the position of the cyclone center, the present application is proposed.

[0040] For example, the present application proposes a cyclone center positioning method, which can determine the position of the cyclone center relatively quickly and simply.

[0041] Before introducing the technical solution of the present application, the technical terms that may be involved in the present application are explained first.

[0042] Cyclone, a meteorological term, refers to a large vortex in the atmosphere in the northern (southern) hemisphere where the horizontal airflow rotates counterclockwise (clockwise). A cyclone is formed by fluctuations on the front or on the interface of air with different densities, and occupies a large-scale horizontal air vortex in three dimensions in the atmosphere. The central air pressure is lower than the surrounding air pressure, and is also called "low pressure" in terms of the air pressure field.

[0043] Cyclones can be divided into extratropical cyclones and tropical cyclones according to their formation and activity range.

[0044] Temperate cyclones are common frontal cyclones, and tropical cyclones are common typhoons and tropical storms. According to the hemispheric distribution, they can be divided into northern hemisphere cyclones and southern hemisphere cyclones. The air in the former rotates counterclockwise, while the air in the latter rotates clockwise. Temperate cyclones mainly occur in eastern Asia, North America, the Mediterranean and other regions. Tropical cyclones are far more common in the northern hemisphere than in the southern hemisphere, and are concentrated between latitudes of 10° and 20°.

[0045] Tropical cyclone is a general term for non-frontal weather-scale vortices that are generated on tropical or subtropical oceans and have organized convection and definite cyclonic circulation. It includes tropical depressions, tropical storms, severe tropical storms, typhoons, severe typhoons and super typhoons. Tropical cyclones are common in the western Pacific and its adjacent waters (typhoons), the Atlantic and northeastern Pacific (hurricanes), and the Indian and South Pacific Oceans.

[0046] Specifically, refer to Figure 1 , shows a flowchart of the steps of a cyclone center positioning method of the present application, the method can be applied to an electronic device, the electronic device can include a terminal or a server, etc., wherein the method can specifically include the following steps:

[0047] In step S101, the wind speed in the first horizontal direction of each of the multiple sub-areas in the area where the cyclone is located is obtained, and the wind speed in the second horizontal direction of each of the multiple sub-areas in the area where the cyclone is located is obtained, and the first horizontal direction is perpendicular to the second horizontal direction.

[0048] In the present application, the wind speed in the first horizontal direction of each of the multiple sub-areas in the area where the cyclone is located and the wind speed in the second horizontal direction of each of the multiple sub-areas in the area where the cyclone is located can be directly obtained from existing meteorological data.

[0049] For example, existing meteorological data can be simulated according to an offshore wind energy resource coupling model. For example, the offshore wind energy resource coupling model can be a coupling model between atmosphere-ocean-waves, etc. Specifically, it can be a coupling model between a meteorological mesoscale model-ocean dynamic model-wave model. The above-mentioned coupling model can output simulated data, such as simulated meteorological mesoscale data, ocean dynamic data, and wave data, etc., wherein the data include the wind speed in the first horizontal direction of each of the multiple sub-regions in the region where the cyclone is located and the wind speed in the second horizontal direction of each of the multiple sub-regions in the region where the cyclone is located.

[0050] For another example, the existing meteorological data can be directly obtained based on the aforementioned cyclone path data in the northwest Pacific from the China Meteorological Administration (CMA-STI), the cyclone path data in the northwest Pacific from the Japan Tokyo Typhoon Center (JMA RSMC Tokyo), and / or the cyclone path data in the northwest Pacific from the United States Joint Typhoon Warning Center (JTWC).

[0051] Among them, for the northwest Pacific, the cyclone path data of the northwest Pacific include but are not limited to three types: the cyclone path data of the northwest Pacific from the China Meteorological Administration (CMA-STI), the cyclone path data of the northwest Pacific from the Tokyo Typhoon Center of Japan (JMA RSMTokyo), and the cyclone path data of the northwest Pacific from the United States Joint Typhoon Warning Center (JTWC).

[0052] All three types of data record the path information of cyclone centers including the northwest Pacific since the 1950s. Any of these three types of data can directly obtain the path information of the cyclone center in the northwest Pacific, and thus the position of the cyclone center at any time.

[0053] However, the time resolution of the cyclone center path information in the three types of data, including the northwest Pacific Ocean, uses a lower time resolution of hours, recording the changes in the position of the cyclone center at a lower time interval. In this way, the lower time interval position of the cyclone center can be determined based on any of the three data.

[0054] The area where a cyclone is located is a large area, for example, an area of ​​several hundred kilometers by several hundred kilometers. The shape of the area where a cyclone is located is uncertain, and it may actually be a variety of shapes, including regular shapes or irregular shapes.

[0055] The shapes of the multiple sub-areas in the area where the cyclone is located may be rectangular, hexagonal, triangular, circular or square, etc. The multiple sub-areas may not overlap.

[0056] The first horizontal direction may be a north-south direction, for example, a direction from south to north or a direction from north to south, and the second horizontal direction may be an east-west direction, for example, a direction from east to west or a direction from west to east.

[0057] For example, see Figure 2 , Figure 2 A schematic diagram showing multiple sub-regions in a region where a cyclone is located.

[0058] exist Figure 2 In FIG. 1 , 25 sub-regions are shown in the region where the cyclone is located. Outside the 25 sub-regions, the region where the cyclone is located may also include more sub-regions. Figure 2 They are no longer shown one by one.

[0059] Figure 2 The shape of each sub-area is a square for illustration, but it is not intended to limit the protection scope of the present application.

[0060] The area of ​​a sub-region 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. Figure 2 The schematic diagram shows the multiple sub-areas in the area where the cyclone is located from a bird's-eye view.

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

[0062] The same is done for each of the other sub-regions, thereby respectively obtaining the overall wind speed change rate of multiple adjacent sub-regions of each sub-region.

[0063] That is, each sub-region corresponds to an overall wind speed change rate, and the overall wind speed change rate corresponding to any sub-region is the overall wind speed change rate of this sub-region with respect to its multiple adjacent sub-regions.

[0064] For any sub-region, the sub-region has multiple adjacent sub-regions, and the sub-region is adjacent to the adjacent sub-regions of the sub-region, and there is no other sub-region between them.

[0065] For example, see Figure 2 For sub-region e, the number of adjacent sub-regions of sub-region e may be 8, namely: 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.

[0066] Each adjacent sub-area of ​​the sub-area has a wind speed, and the wind speed has a direction.

[0067] The overall wind speed change rate of the multiple adjacent sub-regions of the sub-region may reflect: the change amplitude of the wind speed (involving wind direction) of the multiple adjacent sub-regions of the sub-region.

[0068] The same goes for each of the other sub-regions.

[0069] The details of this step may be referred to in the embodiments shown later, and will not be described in detail here.

[0070] In step S103, among the multiple sub-regions in the region where the cyclone is located, the sub-region with the largest overall wind speed change rate is selected.

[0071] In the present application, each sub-region can be sorted in order from high to low according to the overall wind speed change rate of multiple adjacent sub-regions of each sub-region, that is, each sub-region can be sorted in order from high to low according to the overall wind speed change rate corresponding to each sub-region, and then the sub-region ranked first is selected, which is the sub-region with the largest overall wind speed change rate.

[0072] In step S104, the position of the cyclone center is determined according to the sub-region with the largest overall wind speed change rate.

[0073] In the present application, the sub-region with the largest overall wind speed change rate can be determined as the region where the cyclone center is located. In one embodiment, the position of the center of the sub-region with the largest overall wind speed change rate can be determined as the position of the cyclone center, or the position near the center of the sub-region with the largest overall wind speed change rate can be determined as the position of the cyclone center, etc.

[0074] In the present application, the wind speed of each sub-region in the first horizontal direction of multiple sub-regions in the region where the cyclone is located is obtained, and the wind speed of each sub-region in the second horizontal direction of multiple sub-regions in the region where the cyclone is located is obtained, and the first horizontal direction is perpendicular to the second horizontal direction. For any sub-region, the overall wind speed change rate of multiple adjacent sub-regions of the sub-region is determined based on the wind speed of multiple adjacent sub-regions of the sub-region in the first horizontal direction and / or the wind speed of multiple adjacent sub-regions of the sub-region in the second horizontal direction. Among the multiple sub-regions in the region where the cyclone is located, the sub-region with the largest overall wind speed change rate is selected. The position of the cyclone center is determined based on the sub-region with the largest overall wind speed change rate.

[0075] The overall wind speed change rate of multiple adjacent sub-regions of a sub-region can reflect the change range of wind speed (involving wind direction) of multiple adjacent sub-regions of the sub-region. The adjacent sub-regions of different sub-regions are different.

[0076] After statistical analysis of a large number of actual cyclone wind speed data (involving wind direction), it was found that the overall wind speed change rate of multiple adjacent sub-regions of the sub-region where the center of the cyclone is located is greater than the overall wind speed change rate of multiple adjacent sub-regions of the sub-region where the center of the cyclone is located, that is, the change amplitude of the wind speed (involving wind direction) of multiple adjacent sub-regions of the sub-region where the center of the cyclone is located is greater than the change amplitude of the wind speed (involving wind direction) of multiple adjacent sub-regions of the sub-region where the center of the non-cyclone is located.

[0077] Therefore, among the multiple sub-areas in the area where the cyclone is located, the sub-area with the largest overall wind speed change rate can be regarded as the sub-area where the center of the cyclone is located. Therefore, the position of the center of the cyclone can be determined based on the sub-area where the center of the cyclone is located, so as to locate the center of the cyclone.

[0078] Among them, the wind speed in the first horizontal direction of each of the multiple sub-regions in the area where the cyclone is located in the existing meteorological data and the wind speed in the second horizontal direction of each of the multiple sub-regions in the area where the cyclone is located are data with lower time resolution, for example, the resolution of a 1-hour time interval, or even a resolution less than a 1-hour time interval. In this way, through the present application, the time resolution of locating the center of the cyclone can be reduced. For example, the position of the center of the cyclone can be determined within each 1-hour time interval, that is, the position of the center of the cyclone can be determined hour by hour, even less than the position within a 1-hour time interval, which can improve the time resolution of the determined position of the center of the cyclone.

[0079] In summary, the solution for locating the center of a cyclone in the present application has a simple calculation process, consumes less computing resources, and has a high speed, which can improve the efficiency of determining the position of the center of a cyclone.

[0080] In one embodiment of the present application, see Figure 3 , step S102 includes:

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

[0082] In one embodiment of the present application, this step can be implemented through the following process, including:

[0083] 2011. Determine, among multiple adjacent sub-areas of the sub-area, multiple target sub-areas whose wind directions are related to the first horizontal direction.

[0084] In this step, a sub-region where the angle between the wind direction and the first horizontal direction is less than 90 degrees and a sub-region where the angle between the wind direction and the opposite direction of the first horizontal direction is less than 90 degrees can be determined among multiple adjacent sub-regions of the sub-region. The determined sub-regions are determined as target sub-regions.

[0085] For example, based on Figure 2 For the example shown, see Figure 4 , assuming that the first horizontal direction is from west to east, and the second horizontal direction is from south to north. Also, assuming that the sub-region is sub-region e, the number of adjacent sub-regions of sub-region e is 8, namely: 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.

[0086] Figure 4 In the schematic diagram shown, arrows are indicated in sub-area a, sub-area b, sub-area c, sub-area d, sub-area f, sub-area g, sub-area h and sub-area i, and the direction of the arrow in a sub-area refers to the wind direction (wind direction) of this sub-area.

[0087] Thus, the sub-areas where the wind direction is less than 90 degrees from the west to east direction include sub-area g, sub-area h, and sub-area i, and the sub-areas where the wind direction is less than 90 degrees from the west to east direction include sub-area a, sub-area b, and sub-area c.

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

[0089] 2012. Divide the multiple target sub-areas into a first sub-area set and a second sub-area set, wherein a wind direction of any target sub-area in the first sub-area set differs from a wind direction of any target sub-area in the second sub-area set by more than 90 degrees.

[0090] For example, target sub-areas whose wind direction is less than 90 degrees from the first horizontal direction may be combined into a first sub-area set, and target sub-areas whose wind direction is less than 90 degrees from the opposite direction of the first horizontal direction may be combined into a second sub-area set.

[0091] For example, in the above example, sub-region g, sub-region h and sub-region i are combined into a first sub-region set, and sub-region a, sub-region b and sub-region c are combined into a second sub-region set.

[0092] The wind direction of sub-area g differs from the wind direction of sub-area a, the wind direction of sub-area b, and the wind direction of sub-area c by more than 90 degrees.

[0093] The wind direction of sub-area h differs from the wind direction of sub-area a, the wind direction of sub-area b, and the wind direction of sub-area c by more than 90 degrees.

[0094] The wind direction of sub-area i differs from the wind direction of sub-area a, the wind direction of sub-area b, and the wind direction of sub-area c by more than 90 degrees.

[0095] 2013. Perform a weighted operation on the wind speed of the target sub-area in the first sub-area set and the wind speed of the target sub-area in the second sub-area set to obtain a weighted wind speed in the first horizontal direction.

[0096] In one embodiment of the present application, this step can be implemented through the following process, including:

[0097] 11) Perform a weighted operation on the wind speeds of the target sub-areas in the first sub-area set to obtain a weighted wind speed of the first sub-area set.

[0098] In this step, the weight of each target sub-region in the first sub-region set may be determined according to the distance between the sub-region and each target sub-region in the first sub-region set, wherein the weight of the target sub-region in the first sub-region set that has a larger distance from the sub-region has a smaller weight, and the weight of the target sub-region in the first sub-region set that has a smaller distance from the sub-region has a larger weight.

[0099] Furthermore, the number of target sub-regions in the first sub-region set may be counted.

[0100] Then, according to the weight of each target sub-region in the first sub-region set and the number of target sub-regions in the first sub-region set, a weighted operation is performed on the wind speed of each target sub-region in the first sub-region set to obtain the weighted wind speed of the first sub-region set.

[0101] For example, based on the weights of each target sub-area in the first sub-area set, the wind speeds of each target sub-area in the first sub-area set can be weighted and summed to obtain a wind speed sum value, the weights of each target sub-area in the first sub-area set can be summed to obtain a weighted sum value, the product between the wind speed sum value and the weighted sum value can be calculated, and the ratio between the product and the coefficient can be calculated to obtain the weighted wind speed of the first sub-area set.

[0102] The coefficient can be calculated as follows:

[0103] For example, for any target sub-area in the first sub-area set, if the wind speed of the target sub-area can be obtained, the sub-coefficient of the target sub-area is equal to the product of the weight of the target sub-area and a specific value, where the specific value is greater than 0, and the specific value may be 1, 1.5, or 2, etc.; or, if the wind speed of the target sub-area cannot be obtained, the sub-coefficient of the target sub-area is equal to the product of the weight of the target sub-area and a value of 0, that is, the sub-coefficient of the target sub-area is equal to 0.

[0104] The same is true for each other target sub-region in the first sub-region set.

[0105] Then, the sub-coefficients of each target sub-region in the first sub-region set may be summed to obtain the coefficient.

[0106] 12) Perform a weighted operation on the wind speeds of each target sub-area in the second sub-area set to obtain a weighted wind speed of the second sub-area set.

[0107] For details of this step, please refer to step 11), which will not be described in detail here.

[0108] 13) Performing a weighted operation on the weighted wind speed of the first sub-area set and the weighted wind speed of the second sub-area set to obtain a weighted wind speed in the first horizontal direction.

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

[0110] 2014. Determine, according to the weighted wind speed in the first horizontal direction, a wind speed change rate of multiple adjacent sub-areas of the sub-area in the first horizontal direction.

[0111] In the present application, the weighted wind speed in the first horizontal direction may be used as the wind speed change rate of multiple adjacent sub-areas of the sub-area in the first horizontal direction.

[0112] In step S202, the overall wind speed change rate of the plurality of adjacent sub-regions of the sub-region is determined according to the wind speed change rate of the plurality of adjacent sub-regions of the sub-region in the first horizontal direction.

[0113] In one embodiment, the wind speed change rate of multiple adjacent sub-regions of the sub-region in the first horizontal direction can be determined as the overall wind speed change rate of multiple adjacent sub-regions of the sub-region, or the product of the wind speed change rate of multiple adjacent sub-regions of the sub-region in the first horizontal direction and 100% is calculated, and the product is determined as the overall wind speed change rate of multiple adjacent sub-regions of the sub-region, etc.

[0114] In another embodiment of the present application, see Figure 5 , step S102 includes:

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

[0116] In one embodiment of the present application, this step can be implemented through the following process, including:

[0117] 3011. Determine, in a plurality of adjacent sub-areas among the sub-areas, a plurality of target sub-areas whose wind directions are related to the second horizontal direction.

[0118] In this step, a sub-region where the angle between the wind direction and the second horizontal direction is less than 90 degrees and a sub-region where the angle between the wind direction and the opposite direction of the second horizontal direction is less than 90 degrees can be determined among multiple adjacent sub-regions of the sub-region. The determined sub-regions are determined as target sub-regions.

[0119] For example, based on Figure 2 For the example shown, see Figure 4 , assuming that the first horizontal direction is from west to east, and the second horizontal direction is from south to north. Also, assuming that the sub-region is sub-region e, and the number of adjacent sub-regions of sub-region e is 8, namely: 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, Figure 4 In the schematic diagram shown, arrows are indicated in sub-area a, sub-area b, sub-area c, sub-area d, sub-area f, sub-area g, sub-area h and sub-area i, and the direction of the arrow in the sub-area refers to the wind direction (wind direction) of this sub-area.

[0120] Thus, the sub-areas where the wind direction is less than 90 degrees from the south to north direction include sub-area c, sub-area f, and sub-area i, and the sub-areas where the wind direction is less than 90 degrees from the west to east direction opposite to the west to east direction include sub-area a, sub-area d, and sub-area g.

[0121] Therefore, sub-region a, sub-region d, sub-region g, sub-region c, sub-region f, and sub-region g are target sub-regions.

[0122] 3012. Divide the multiple target sub-areas into a third sub-area set and a fourth sub-area set, wherein a wind direction of any target sub-area in the third sub-area set differs from a wind direction of any target sub-area in the fourth sub-area set by more than 90 degrees.

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

[0124] 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.

[0125] The wind direction of sub-area c differs from the wind direction of sub-area a, the wind direction of sub-area d, and the wind direction of sub-area g by more than 90 degrees.

[0126] The wind direction of sub-area f differs from the wind direction of sub-area a, the wind direction of sub-area d, and the wind direction of sub-area g by more than 90 degrees.

[0127] The wind direction of the sub-area i differs from the wind direction of the sub-area a, the wind direction of the sub-area d, and the wind direction of the sub-area g by more than 90 degrees.

[0128] 3013. Perform a weighted operation on the wind speed of the target sub-area in the third sub-area set and the wind speed of the target sub-area in the fourth sub-area set to obtain a weighted wind speed in the second horizontal direction.

[0129] In one embodiment of the present application, this step can be implemented through the following process, including:

[0130] 21) Perform weighted calculation on the wind speeds of the target sub-areas in the third sub-area set to obtain the weighted wind speed of the third sub-area set.

[0131] In this step, the weight of each target sub-region in the third sub-region set may be determined according to the distance between the sub-region and each target sub-region in the third sub-region set, wherein the weight of the target sub-region in the third sub-region set that has a larger distance from the sub-region has a smaller weight, and the weight of the target sub-region in the third sub-region set that has a smaller distance from the sub-region has a larger weight.

[0132] Furthermore, the number of target sub-regions in the third sub-region set may be counted.

[0133] Then, according to the weight of each target sub-region in the third sub-region set and the number of target sub-regions in the third sub-region set, the wind speed of each target sub-region in the third sub-region set is weighted to obtain the weighted wind speed of the third sub-region set.

[0134] For example, according to the weights of each target sub-area in the third sub-area set, the wind speeds of each target sub-area in the third sub-area set can be weighted and summed to obtain a wind speed sum value, the weights of each target sub-area in the third sub-area set can be summed to obtain a weighted sum value, the product between the wind speed sum value and the weighted sum value can be calculated, and the ratio between the product and the coefficient can be calculated to obtain the weighted wind speed of the third sub-area set.

[0135] The coefficient can be calculated as follows:

[0136] For example, for any target sub-area in the third sub-area set, if the wind speed of the target sub-area can be obtained, the sub-coefficient of the target sub-area is equal to the product of the weight of the target sub-area and a specific value, where the specific value is greater than 0, and the specific value may be 1, 1.5, or 2, etc.; or, if the wind speed of the target sub-area cannot be obtained, the sub-coefficient of the target sub-area is equal to the product of the weight of the target sub-area and a value of 0, that is, the sub-coefficient of the target sub-area is equal to 0.

[0137] The same is true for each other target sub-region in the third sub-region set.

[0138] Then, the sub-coefficients of each target sub-region in the third sub-region set may be summed to obtain the coefficient.

[0139] 22) Perform a weighted operation on the wind speeds of the target sub-areas in the fourth sub-area set to obtain a weighted wind speed of the fourth sub-area set.

[0140] For details of this step, please refer to step 21), which will not be described in detail here.

[0141] 23) According to the weighted wind speed of the third sub-area set and the weighted wind speed of the fourth sub-area set, the weighted wind speed in the second horizontal direction is obtained.

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

[0143] 3014. Determine, based on the weighted wind speed in the second horizontal direction, a wind speed change rate of a plurality of adjacent sub-areas of the sub-area in the second horizontal direction.

[0144] In the present application, the weighted wind speed in the first horizontal direction may be used as the wind speed change rate of multiple adjacent sub-areas of the sub-area in the first horizontal direction.

[0145] In step S302, the overall wind speed change rate of the plurality of adjacent sub-regions of the sub-region is determined according to the wind speed change rate of the plurality of adjacent sub-regions of the sub-region in the second horizontal direction.

[0146] In one embodiment, the wind speed change rate of multiple adjacent sub-regions of the sub-region in the first horizontal direction can be determined as the overall wind speed change rate of multiple adjacent sub-regions of the sub-region, or the product of the wind speed change rate of multiple adjacent sub-regions of the sub-region in the first horizontal direction and 100% is calculated, and the product is determined as the overall wind speed change rate of multiple adjacent sub-regions of the sub-region, etc.

[0147] In another embodiment of the present application, see Figure 6 , step S102 includes:

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

[0149] For details of this step, please refer to step S201, which will not be described in detail here.

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

[0151] For details of this step, please refer to step S301, which will not be described in detail here.

[0152] This application does not limit the execution order of step S401 and step S402.

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

[0154] 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.

[0155] In one embodiment of the present application, the first square of the wind speed change rate of multiple adjacent sub-regions of the sub-region in the first horizontal direction can be calculated, and the second square of the wind speed change rate of multiple adjacent sub-regions of the sub-region in the second horizontal direction can be calculated, and then the sum of the first square and the second square can be calculated, and then the square root of the sum can be calculated, and then the overall wind speed change rate of the multiple adjacent sub-regions of the sub-region can be determined based on the square root. For example, the square root is directly determined as the overall wind speed change rate of the multiple adjacent sub-regions of the sub-region, or the product between the square root and 100% is calculated, and the product is determined as the overall wind speed change rate of the multiple adjacent sub-regions of the sub-region, etc.

[0156] The solution of the present application is illustrated by a specific embodiment, but it is not intended to limit the protection scope of the present application.

[0157] (1) Download the "ERA5 hourly data on single levels from 1940 to present" dataset from the ERA5 reanalysis data. Select the "10m u-component of wind" and "10m v-component of wind" near sea level u wind speed (east-west wind speed) and v wind speed (north-south wind speed) respectively, select the year, month, day and hour of the tropical cyclone, select the longitude and latitude range of the tropical cyclone's transit impact area, and select the NetCDF format for download.

[0158] (2) In ArcGIS software, convert the u wind speed and v wind speed components in NetCDF format into TIFF raster formats respectively.

[0159] (3) Calculate the slope for the u wind speed grid and the v wind speed grid respectively.

[0160] The slope calculation result is expressed as a percentage (percentage of elevation gain).

[0161] Among them, the slope (elevation) increment percentage is equal to the elevation increment divided by the horizontal increment and then multiplied by 100 (for example: when the slope angle is 45 degrees, the elevation increment is equal to the horizontal increment, so the elevation increment percentage is 100%; when the slope angle approaches a right angle of 90 degrees, the elevation increment percentage begins to approach infinity).

[0162] Slope calculation uses a planar method, and the slope is measured as the maximum rate of change in value from one pixel to its directly adjacent pixels. The calculation is performed on a projected plane using a 2D Cartesian coordinate system. The slope value is calculated using the maximum mean method (Burrough, 1998). A neighborhood (moving window) of 3*3 pixels is used to perform planar and geodesic calculations. For each neighborhood, if the processing (center) pixel is NoData (lack of valid data), the output is also NoData (lack of valid data). The calculation also requires that at least 7 pixels adjacent to the processing pixel have valid values. If there are fewer than 7 valid pixels, the calculation is not performed, and the output at the processing pixel will be NoData (lack of valid data). The tool calculates the maximum rate of change in value from the pixel to its adjacent pixels, and the slope (increment) percentage depends on the rate of change of the surface in the horizontal (dz / dx) and vertical (dz / dy) directions starting from the center pixel:

[0163]

[0164] like Figure 2 As shown in Figure 4, the values ​​of the central pixel e and its eight adjacent pixels determine the horizontal and vertical increments. These adjacent pixels are represented by letters a to i, where pixel e represents the pixel for which the aspect is currently being calculated.

[0165] The rate of change of pixel e in the x direction is equal to:

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

[0167] In the above formula, wght1 and wght2 are the horizontal weighted counts of valid pixels. When c, f, and i all have a valid value, then wght1 = (1+2*1+1) = 4; when i is NoData, then wght1 = (1+2*1+0) = 3; when f is NoData, wght1 = (1+2*0+1) = 2. Similar logic applies to wght2 except that the neighborhood positions are a, d, and g. || represents the absolute value sign. x_cellsize represents the size of the pixel. 8 means that the center pixel e has 8 neighboring pixels.

[0168] The rate of change of pixel e in the y direction is equal to:

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

[0170] In the above formula, wght3 and wght4 are the same concepts as in the [dz / dx] calculation.

[0171] (4) The u wind speed slope percentage and the v wind speed slope (elevation increment) percentage are added to the grid to highlight the pixel value of the tropical cyclone center relative to the surrounding area.

[0172] (5) The pixel with the largest median value in the grid obtained in step (4) is the center position of the tropical cyclone at that hour.

[0173] Through the spatial analysis of the sea level u wind speed and v wind speed raster data of the reanalysis data, the position of the tropical cyclone center at each hour can be obtained relatively simply, quickly and accurately. The hourly tropical cyclone center position between every adjacent 6 hours in the existing tropical cyclone best path data products can be filled in, which can improve the level of tropical cyclone warning forecast and disaster assessment.

[0174] It should be noted that, for the method embodiments, for the sake of simplicity, they are all expressed as a series of action combinations, but those skilled in the art should be aware that the present application is not limited by the order of the actions described, because according to the present application, certain steps can be performed in other orders or simultaneously. Secondly, those skilled in the art should also be aware that the embodiments described in the specification are all optional embodiments, and the actions involved are not necessarily required by the present application.

[0175] Reference Figure 7 , shows a structural block diagram of a cyclone center positioning device of the present application, the device comprising:

[0176] An acquisition module 11 is used to acquire the wind speed of each sub-region in a plurality of sub-regions in the region where the cyclone is located in a first horizontal direction, and to acquire the wind speed of each sub-region in a plurality of sub-regions in the region where the cyclone is located in a second horizontal direction, wherein the first horizontal direction is perpendicular to the second horizontal direction;

[0177] A first determination module 12 is used to determine, for any sub-area, the overall wind speed change rate of multiple adjacent sub-areas of the sub-area according to the wind speed of multiple adjacent sub-areas of the sub-area in the first horizontal direction and / or the wind speed of multiple adjacent sub-areas of the sub-area in the second horizontal direction;

[0178] A selection module 13 is used to select a sub-region with the largest overall wind speed change rate from among multiple sub-regions in the region where the cyclone is located;

[0179] The second determination module 14 is used to determine the position of the cyclone center according to the sub-region with the largest overall wind speed change rate.

[0180] In an optional implementation, the first determining module includes:

[0181] A first determining unit, configured to determine a wind speed change rate of a plurality of adjacent sub-areas of the sub-area in the first horizontal direction according to wind speeds of a plurality of adjacent sub-areas of the sub-area in the first horizontal direction;

[0182] The second determining unit is used to determine the overall wind speed change rate of the plurality of adjacent sub-areas of the sub-area according to the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction.

[0183] In an optional implementation, the first determining module includes:

[0184] A third determining unit, configured to determine a wind speed change rate of a plurality of adjacent sub-areas of the sub-area in the second horizontal direction according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction;

[0185] The fourth determining unit is used to determine the overall wind speed change rate of the plurality of adjacent sub-areas of the sub-area according to the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction.

[0186] In an optional implementation, the first determining module includes:

[0187] A first determining unit, configured to determine a wind speed change rate of a plurality of adjacent sub-areas of the sub-area in the first horizontal direction according to wind speeds of a plurality of adjacent sub-areas of the sub-area in the first horizontal direction;

[0188] A third determining unit, configured to determine a wind speed change rate of a plurality of adjacent sub-areas of the sub-area in the second horizontal direction according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction;

[0189] The fifth determination unit is used to determine the overall wind speed change rate of multiple adjacent sub-areas of the sub-area based on the wind speed change rate of multiple adjacent sub-areas of the sub-area in the first horizontal direction and the wind speed change rate of multiple adjacent sub-areas of the sub-area in the second horizontal direction.

[0190] In an optional implementation, the first determining unit includes:

[0191] A first determining subunit is used to determine, among a plurality of adjacent sub-areas of the sub-areas, a plurality of target sub-areas whose wind directions are associated with a first horizontal direction;

[0192] A first division subunit is used to divide the multiple target sub-areas into a first sub-area set and a second sub-area set, and the wind direction of any target sub-area in the first sub-area set differs from the wind direction of any target sub-area in the second sub-area set by more than 90 degrees;

[0193] A first weighted operation subunit is used to perform a weighted operation on the wind speed of the target sub-area in the first sub-area set and the wind speed of the target sub-area in the second sub-area set to obtain a weighted wind speed in a first horizontal direction;

[0194] The second determining subunit is used to determine the wind speed change rate of multiple adjacent sub-areas of the sub-area in the first horizontal direction according to the weighted wind speed in the first horizontal direction.

[0195] In an optional implementation, the first determination subunit is specifically used to: determine, among multiple adjacent sub-areas of the sub-area, a sub-area where the angle between the wind direction and the first horizontal direction is less than 90 degrees, and a sub-area where the angle between the wind direction and the opposite direction of the first horizontal direction is less than 90 degrees; and determine the determined sub-area as the target sub-area.

[0196] In an optional implementation, the first weighted operation subunit is specifically used to: perform weighted operation on the wind speeds of each target sub-area in the first sub-area set to obtain the weighted wind speed of the first sub-area set; perform weighted operation on the wind speeds of each target sub-area in the second sub-area set to obtain the weighted wind speed of the second sub-area set; perform weighted operation on the weighted wind speed of the first sub-area set and the weighted wind speed of the second sub-area set to obtain the weighted wind speed in the first horizontal direction.

[0197] In an optional implementation, the first weighted operation subunit is specifically used to: determine the weight of each target sub-region in the first sub-region set based on the distance between the sub-region and each target sub-region in the first sub-region set; count the number of target sub-regions in the first sub-region set; and perform weighted operation on the wind speed of each target sub-region in the first sub-region set based on the weight of each target sub-region in the first sub-region set and the number of target sub-regions in the first sub-region set to obtain the weighted wind speed of the first sub-region set.

[0198] In an optional implementation, the third determining unit includes:

[0199] A third determining subunit is used to determine, among the plurality of adjacent sub-areas of the sub-area, a plurality of target sub-areas whose wind directions are associated with the second horizontal direction;

[0200] A second division subunit is used to divide the multiple target sub-areas into a third sub-area set and a fourth sub-area set, and the wind direction of any target sub-area in the third sub-area set differs from the wind direction of any target sub-area in the fourth sub-area set by more than 90 degrees;

[0201] A second weighted operation subunit is used to perform a weighted operation on the wind speed of the target sub-area in the third sub-area set and the wind speed of the target sub-area in the fourth sub-area set to obtain a weighted wind speed in a second horizontal direction;

[0202] The fourth determining subunit is used to determine the wind speed change rate of multiple adjacent sub-areas of the sub-area in the second horizontal direction according to the weighted wind speed in the second horizontal direction.

[0203] In an optional implementation, the third determination subunit is specifically used to: determine, among multiple adjacent sub-areas of the sub-area, a sub-area in which the wind direction is less than 90 degrees from the second horizontal direction, and a sub-area in which the wind direction is less than 90 degrees from the opposite direction of the second horizontal direction; and determine the determined sub-area as the target sub-area.

[0204] In an optional implementation, the second weighted operation subunit is specifically used to: perform weighted operation on the wind speeds of each target sub-area in the third sub-area set to obtain the weighted wind speed of the third sub-area set; perform weighted operation on the wind speeds of each target sub-area in the fourth sub-area set to obtain the weighted wind speed of the fourth sub-area set; perform weighted operation on the weighted wind speed of the third sub-area set and the weighted wind speed of the fourth sub-area set to obtain the weighted wind speed in the second horizontal direction.

[0205] In an optional implementation, the second weighted operation subunit is specifically used to: determine the weight of each target sub-area in the third sub-area set based on the distance between the sub-area and each target sub-area in the third sub-area set; count the number of target sub-areas in the third sub-area set; and perform weighted operation on the wind speed of each target sub-area in the third sub-area set based on the weight of each target sub-area in the third sub-area set and the number of target sub-areas in the third sub-area set to obtain the weighted wind speed of the third sub-area set.

[0206] In an optional implementation, the fifth determining unit includes:

[0207] A first calculation subunit, configured to calculate 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;

[0208] A second calculation subunit, configured to calculate 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;

[0209] A third calculation subunit, configured to calculate a sum of the first square and the second square;

[0210] A fourth calculation subunit, used for calculating the square root of the sum;

[0211] The fifth determining subunit is used to determine the overall wind speed change rate of multiple adjacent sub-areas of the sub-area according to the square root.

[0212] In the present application, the wind speed of each sub-region in the first horizontal direction of multiple sub-regions in the region where the cyclone is located is obtained, and the wind speed of each sub-region in the second horizontal direction of multiple sub-regions in the region where the cyclone is located is obtained, and the first horizontal direction is perpendicular to the second horizontal direction. For any sub-region, the overall wind speed change rate of multiple adjacent sub-regions of the sub-region is determined based on the wind speed of multiple adjacent sub-regions of the sub-region in the first horizontal direction and / or the wind speed of multiple adjacent sub-regions of the sub-region in the second horizontal direction. Among the multiple sub-regions in the region where the cyclone is located, the sub-region with the largest overall wind speed change rate is selected. The position of the cyclone center is determined based on the sub-region with the largest overall wind speed change rate.

[0213] The overall wind speed change rate of multiple adjacent sub-regions of a sub-region can reflect the change range of wind speed (involving wind direction) of multiple adjacent sub-regions of the sub-region. The adjacent sub-regions of different sub-regions are different.

[0214] After statistical analysis of a large number of actual cyclone wind speed data (involving wind direction), it was found that the overall wind speed change rate of multiple adjacent sub-regions of the sub-region where the center of the cyclone is located is greater than the overall wind speed change rate of multiple adjacent sub-regions of the sub-region where the center of the cyclone is located, that is, the change amplitude of the wind speed (involving wind direction) of multiple adjacent sub-regions of the sub-region where the center of the cyclone is located is greater than the change amplitude of the wind speed (involving wind direction) of multiple adjacent sub-regions of the sub-region where the center of the non-cyclone is located.

[0215] Therefore, among the multiple sub-areas in the area where the cyclone is located, the sub-area with the largest overall wind speed change rate can be regarded as the sub-area where the center of the cyclone is located. Therefore, the position of the center of the cyclone can be determined based on the sub-area where the center of the cyclone is located, so as to locate the center of the cyclone.

[0216] Among them, the wind speed in the first horizontal direction of each of the multiple sub-regions in the area where the cyclone is located in the existing meteorological data and the wind speed in the second horizontal direction of each of the multiple sub-regions in the area where the cyclone is located are data with lower time resolution, for example, the resolution of a 1-hour time interval, or even a resolution less than a 1-hour time interval. In this way, through the present application, the time resolution of locating the center of the cyclone can be reduced. For example, the position of the center of the cyclone can be determined within each 1-hour time interval, that is, the position of the center of the cyclone can be determined hour by hour, even less than the position within a 1-hour time interval, which can improve the time resolution of the determined position of the center of the cyclone.

[0217] In summary, the solution for locating the center of a cyclone in the present application has a simple computational process, consumes few computing resources, and has a high speed, which can improve the efficiency of determining the location of the center of a cyclone.

[0218] As for the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0219] Optionally, an embodiment of the present application further provides an electronic device, comprising: a processor, a memory, and a computer program stored in the memory and executable on the processor. When the computer program is executed by the processor, the various processes of the above-mentioned method embodiment are implemented, and the same technical effect can be achieved. To avoid repetition, it will not be described here.

[0220] The embodiment of the present application also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, each process of the above method embodiment is implemented, and the same technical effect can be achieved. To avoid repetition, it is not repeated here. The computer-readable storage medium is, for example, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0221] Figure 8 800 is a block diagram of an electronic device 800 shown in the present application. For example, the electronic device 800 may be a mobile phone, a computer, a digital broadcast terminal, a messaging device, a game console, a tablet device, a medical device, a fitness device, a personal digital assistant, etc.

[0222] Reference Figure 8 , the 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 .

[0223] The processing component 802 generally controls the overall operation of the electronic device 800, such as operations associated with display, phone calls, data communications, camera operations, and recording operations. The processing component 802 may include one or more processors 820 to execute instructions to complete all or part of the steps of the above method. In addition, the processing component 802 may include one or more modules to facilitate the interaction between the processing component 802 and other components. For example, the processing component 802 may include a multimedia module to facilitate the interaction between the multimedia component 808 and the processing component 802.

[0224] The memory 804 is configured to store various types of data to support operations on the electronic device 800. Examples of such data include instructions for any application or method operating on the electronic device 800, contact data, phone book data, messages, images, videos, etc. The memory 804 may be implemented by any type of volatile or non-volatile storage device or a combination thereof, 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 disk, or optical disk.

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

[0226] The multimedia component 808 includes a screen that provides an output interface between the electronic device 800 and the user. 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 the user. The touch panel includes one or more touch sensors to sense touch, slide, and gestures on the touch panel. The touch sensor may not only sense the boundaries of the touch or slide action, but also detect 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 operating 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 a focal length and optical zoom capability.

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

[0228] I / O interface 812 provides an interface between processing component 802 and peripheral interface modules, such as keyboards, click wheels, buttons, etc. These buttons may include but are not limited to: home button, volume button, start button, and lock button.

[0229] The sensor assembly 814 includes one or more sensors for providing various aspects of status assessment for the electronic device 800. For example, the sensor assembly 814 can detect the open / closed state of the electronic device 800, the relative positioning of the components, such as the display and keypad of the electronic device 800, and the sensor assembly 814 can also detect the position change of the electronic device 800 or a component of the electronic device 800, the presence or absence of contact between the user and the electronic device 800, the orientation or acceleration / deceleration of the electronic device 800, and the temperature change of the electronic device 800. The sensor assembly 814 may include a proximity sensor configured to detect the presence of a nearby object without any physical contact. The sensor assembly 814 may also include an optical 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 gyroscope sensor, a magnetic sensor, a pressure sensor, or a temperature sensor.

[0230] The communication component 816 is configured to facilitate wired or wireless communication between the electronic device 800 and other devices. The electronic device 800 can access a wireless network based on a communication standard, such as WiFi, a carrier network (such as 2G, 3G, 4G or 5G), or a combination thereof. In an exemplary embodiment, the communication component 816 receives a broadcast signal or broadcast operation information from an external broadcast management system via a broadcast channel. In an exemplary embodiment, the communication component 816 also includes a near field communication (NFC) module to facilitate short-range communication. For example, the NFC module can be implemented based on radio frequency identification (RFID) technology, infrared data association (IrDA) technology, ultra-wideband (UWB) technology, Bluetooth (BT) technology and other technologies.

[0231] 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 components to perform the above methods.

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

[0233] Fig. 9 1 is a block diagram of an electronic device 1900 shown in the present application. For example, the electronic device 1900 may be provided as a server.

[0234] Reference Fig. 9 , the electronic device 1900 includes a processing component 1922, which further includes one or more processors, and a memory resource represented by a memory 1932 for storing instructions executable by the processing component 1922, such as an application. The application stored in the memory 1932 may include one or more modules, each corresponding to a set of instructions. In addition, the processing component 1922 is configured to execute instructions to perform the above method.

[0235] The electronic device 1900 may also include a power supply 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 the memory 1932, such as Windows Server™, Mac OS X™, Unix™, Linux™, FreeBSD™ or the like.

[0236] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the existence of other identical elements in the process, method, article or device including the element.

[0237] Through the description of the above implementation methods, those skilled in the art can clearly understand that the above-mentioned embodiment methods can be implemented by means of software plus a necessary general hardware platform, and of course by hardware, but in many cases the former is a better implementation method. Based on such an understanding, the technical solution of the present application, or the part that contributes to the prior art, can be embodied in the form of a software product, which is stored in a storage medium (such as ROM / RAM, a magnetic disk, or an optical disk), and includes a number of instructions for a terminal (which can be a mobile phone, a computer, a server, an air conditioner, or a network device, etc.) to execute the methods described in each embodiment of the present application.

[0238] The embodiments of the present application are described above in conjunction with the accompanying drawings, but the present application is not limited to the above-mentioned specific implementation methods. The above-mentioned specific implementation methods are merely illustrative and not restrictive. Under the guidance of the present application, ordinary technicians in this field can also make many forms without departing from the purpose of the present application and the scope of protection of the claims, all of which are within the protection of the present application.

[0239] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed in the present application can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application.

[0240] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0241] In the embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0242] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple 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.

[0243] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0244] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the methods described in each embodiment of the present application. The aforementioned storage medium includes: various media that can store program codes, such as USB flash drives, mobile hard drives, ROM, RAM, magnetic disks, or optical disks.

[0245] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art who is familiar with the present technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

Claims

1. A method for locating the center of a cyclone, characterized in that: The method comprises: Obtaining the wind speed of each of the multiple sub-regions in the region where the cyclone is located in a first horizontal direction, and obtaining the wind speed of each of the multiple sub-regions in the region where the cyclone is located in a second horizontal direction, the first horizontal direction being perpendicular to the second horizontal direction; For any sub-region, determining the overall wind speed change rate of the multiple adjacent sub-regions of the sub-region according to the wind speed of the multiple adjacent sub-regions of the sub-region in the first horizontal direction and / or the wind speed of the multiple adjacent sub-regions of the sub-region in the second horizontal direction; Among the multiple sub-areas in the area where the cyclone is located, select the sub-area with the largest overall wind speed change rate; The location of the cyclone center is determined based on the sub-area with the largest overall wind speed change rate.

2. The method according to claim 1, characterized in that The determining the overall wind speed change rate of the plurality of adjacent sub-areas of the sub-area according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction and / or the wind speed of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction comprises: Determining the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction; 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 rates of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction.

3. The method according to claim 1, characterized in that The determining the overall wind speed change rate of the plurality of adjacent sub-areas of the sub-area according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction and / or the wind speed of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction comprises: Determining the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction; 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 rates of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction.

4. The method according to claim 1, characterized in that: The determining the overall wind speed change rate of the plurality of adjacent sub-areas of the sub-area according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction and / or the wind speed of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction comprises: Determining the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction; Determining the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction; 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.

5. The method according to claim 2 or 4, characterized in that: The step of determining the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the first horizontal direction comprises: Determining, among the plurality of adjacent sub-areas of the sub-area, a plurality of target sub-areas whose wind directions are associated with the first horizontal direction; The plurality of target sub-areas are divided into a first sub-area set and a second sub-area set, wherein a wind direction of any target sub-area in the first sub-area set differs from a wind direction of any target sub-area in the second sub-area set by more than 90 degrees; Performing a weighted operation on the wind speed of the target sub-area in the first sub-area set and the wind speed of the target sub-area in the second sub-area set to obtain a weighted wind speed in a first horizontal direction; According to the weighted wind speed in the first horizontal direction, the wind speed change rates of a plurality of adjacent sub-areas of the sub-area in the first horizontal direction are determined.

6. The method according to claim 5, characterized in that Determining, among the plurality of adjacent sub-areas of the sub-area, a plurality of target sub-areas whose wind directions are related to the first horizontal direction comprises: Among the plurality of adjacent sub-areas of the sub-areas, determine a sub-area where the angle between the wind direction and the first horizontal direction is less than 90 degrees, and determine a sub-area where the angle between the wind direction and the opposite direction of the first horizontal direction is less than 90 degrees; The determined sub-region is determined as the target sub-region.

7. The method according to claim 5, characterized in that The step of performing a weighted operation on the wind speed of the target sub-area in the first sub-area set and the wind speed of the target sub-area in the second sub-area set to obtain the weighted wind speed in the first horizontal direction includes: Performing a weighted calculation on the wind speeds of each target sub-area in the first sub-area set to obtain a weighted wind speed of the first sub-area set; Performing a weighted calculation on the wind speeds of each target sub-area in the second sub-area set to obtain a weighted wind speed of the second sub-area set; A weighted operation is performed on the weighted wind speed of the first sub-area set and the weighted wind speed of the second sub-area set to obtain a weighted wind speed in the first horizontal direction.

8. The method according to claim 7, characterized in that The step of performing a weighted operation on the wind speeds of the target sub-areas in the first sub-area set to obtain the weighted wind speeds of the first sub-area set includes: Determining the weight of each target sub-region in the first sub-region set according to the distance between the sub-region and each target sub-region in the first sub-region set; Counting the number of target sub-regions in the first sub-region set; According to the weight of each target sub-region in the first sub-region set and the number of target sub-regions in the first sub-region set, a weighted operation is performed on the wind speed of each target sub-region in the first sub-region set to obtain the weighted wind speed of the first sub-region set.

9. The method according to claim 3 or 4, characterized in that: The step of determining the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction according to the wind speed of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction comprises: Determining, among the plurality of adjacent sub-areas of the sub-area, a plurality of target sub-areas whose wind directions are associated with the second horizontal direction; The plurality of target sub-areas are divided into a third sub-area set and a fourth sub-area set, wherein a wind direction of any target sub-area in the third sub-area set differs from a wind direction of any target sub-area in the fourth sub-area set by more than 90 degrees; Performing a weighted operation on the wind speed of the target sub-area in the third sub-area set and the wind speed of the target sub-area in the fourth sub-area set to obtain a weighted wind speed in the second horizontal direction; According to the weighted wind speed in the second horizontal direction, the wind speed change rates of a plurality of adjacent sub-areas of the sub-area in the second horizontal direction are determined.

10. The method according to claim 9, characterized in that Determining, among the plurality of adjacent sub-areas of the sub-area, a plurality of target sub-areas whose wind directions are related to the second horizontal direction comprises: Among the plurality of adjacent sub-areas of the sub-areas, determine a sub-area where the angle between the wind direction and the second horizontal direction is less than 90 degrees, and determine a sub-area where the angle between the wind direction and the opposite direction of the second horizontal direction is less than 90 degrees; The determined sub-region is determined as the target sub-region.

11. The method according to claim 9, characterized in that The step of performing a weighted operation on the wind speed of the target sub-area in the third sub-area set and the wind speed of the target sub-area in the fourth sub-area set to obtain the weighted wind speed in the second horizontal direction includes: Performing a weighted calculation on the wind speeds of each target sub-area in the third sub-area set to obtain a weighted wind speed of the third sub-area set; Performing a weighted calculation on the wind speeds of each target sub-area in the fourth sub-area set to obtain a weighted wind speed of the fourth sub-area set; A weighted operation is performed on the weighted wind speed of the third sub-area set and the weighted wind speed of the fourth sub-area set to obtain a weighted wind speed in the second horizontal direction.

12. The method according to claim 11, characterized in that The step of performing a weighted operation on the wind speeds of the target sub-areas in the third sub-area set to obtain the weighted wind speeds of the third sub-area set includes: determining a weight of each target sub-region in the third sub-region set according to a distance between the sub-region and each target sub-region in the third sub-region set; Counting the number of target sub-regions in the third sub-region set; According to the weight of each target sub-region in the third sub-region set and the number of target sub-regions in the third sub-region set, a weighted operation is performed on the wind speed of each target sub-region in the third sub-region set to obtain the weighted wind speed of the third sub-region set.

13. The method according to claim 4, characterized in that The determining the overall wind speed change rate of the plurality of adjacent sub-areas of the sub-area according to the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction and the wind speed change rate of the plurality of adjacent sub-areas of the sub-area in the second horizontal direction comprises: Calculate 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; Calculate the second square of the wind speed change rate of multiple adjacent sub-areas of the sub-area in the second horizontal direction; Calculate the sum of the first square and the second square; calculating the square root of the sum; The overall wind speed change rate of a plurality of adjacent sub-areas of the sub-area is determined according to the square root.

14. A cyclone center positioning device, characterized in that: The device comprises: an acquisition module, configured to acquire the wind speed of each of the multiple sub-regions in the region where the cyclone is located in a first horizontal direction, and to acquire the wind speed of each of the multiple sub-regions in the region where the cyclone is located in a second horizontal direction, wherein the first horizontal direction is perpendicular to the second horizontal direction; A first determination module is used to determine, for any sub-area, the overall wind speed change rate of multiple adjacent sub-areas of the sub-area according to the wind speed of multiple adjacent sub-areas of the sub-area in the first horizontal direction and / or the wind speed of multiple adjacent sub-areas of the sub-area in the second horizontal direction; A selection module is used to select a sub-area with the largest overall wind speed change rate from among multiple sub-areas in the area where the cyclone is located; The second determination module is used to determine the position of the cyclone center according to the sub-area with the largest overall wind speed change rate.

15. An electronic device, characterized in that: include: A processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the computer program implements the method according to any one of claims 1 to 13 when executed by the processor.

16. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 13 is implemented.

Images