An automatic typhoon center positioning method and device

Abstract

The present application provides a method and device for automatically locating the typhoon center. The method includes: obtaining the current forecast typhoon center, forecast typhoon path information, and sea surface wind field observation data of a spaceborne microwave scatterometer; when the sea surface wind field observation data includes the current forecast typhoon center, extracting the current observation time when the current forecast typhoon center is observed; obtaining an interpolated forecast typhoon center corresponding to the current observation time through interpolation; dividing the interpolated forecast typhoon center into regions to obtain a typhoon center positioning region including a plurality of grid nodes; in the typhoon center positioning region, determining the positioning grid node where the simulated wind direction data is most similar to the observed wind direction data; determining the positioning grid node as the positioned typhoon center, and extracting the longitude and latitude information of the positioned typhoon center. It can be seen that this method can automatically locate the typhoon center using a satellite scatterometer based on typhoon prior knowledge, thus solving the problem of difficult typhoon center positioning by the satellite scatterometer.

Description

Technical Field

[0001] This application relates to the field of marine microwave remote sensing technology, and more particularly, to a method and device for automatically locating the typhoon center. Background Art

[0002] Currently, it is still difficult to locate the typhoon center using satellite scatterometers. Specifically, since the satellite scatterometer observations are affected by rainfall factors, during the process of typhoon accompanied by rainfall, the initial field of the scatterometer's ambiguous solution usually shows patchy 180° ambiguity. As a result, it is still difficult to locate the typhoon center. Summary of the Invention

[0003] The purpose of the embodiments of this application is to provide a method and device for automatically locating the typhoon center, which can automatically locate the typhoon center using a satellite scatterometer based on typhoon prior knowledge, thereby solving the problem of difficult typhoon center location using satellite scatterometers.

[0004] The first aspect of the embodiments of this application provides a method for automatically locating the typhoon center, including:

[0005] Obtain the current forecast typhoon center and forecast typhoon path information released by the meteorological center; observe the sea surface wind field observation data through a satellite scatterometer;

[0006] Detect whether the sea surface wind field observation data includes the current forecast typhoon center that matches the current forecast typhoon center;

[0007] When the sea surface wind field observation data includes the current forecast typhoon center, extract the current observation time when the current forecast typhoon center is detected;

[0008] Interpolate in the forecast typhoon path information to obtain an interpolated forecast typhoon center corresponding to the current observation time;

[0009] Based on a preset regional division standard and the interpolated forecast typhoon center, perform regional division to obtain a typhoon center location area including multiple grid nodes;

[0010] In the typhoon center location area, determine the location grid node; the simulated wind direction data corresponding to the location grid node is most similar to the observed wind direction data;

[0011] Determine the location grid node as the located typhoon center, and extract the longitude and latitude information of the located typhoon center.

[0012] Furthermore,

[0013] The step of determining the location grid node in the typhoon center location area includes:

[0014] Determine multiple of the grid node positions as multiple guessed typhoon centers;

[0015] Perform calculation and simulation according to the Holland model and each of the guessed typhoon centers to obtain multiple sets of simulated wind direction data;

[0016] Obtain the observed wind direction data in the sea surface wind field observation data;

[0017] Determine a set of simulated wind direction data that is most similar to the observed wind direction data among the multiple sets of simulated wind direction data as the similar wind direction data;

[0018] Determine the grid node corresponding to the similar wind direction data as the positioning grid node corresponding to the interpolated forecast typhoon center.

[0019] Further, the determining the set of simulated wind direction data that is most similar to the observed wind direction data among the multiple sets of simulated wind direction data as the similar wind direction data includes:

[0020] Perform calculation according to each set of the simulated wind direction data and the observed wind direction data to obtain multiple mean absolute values of wind direction differences;

[0021] Extract the minimum mean absolute value of wind direction differences among the multiple mean absolute values of wind direction differences, and determine the set of simulated wind direction data corresponding to the minimum mean absolute value of wind direction differences as the similar wind direction data.

[0022] Further, after the step of obtaining the current forecast typhoon center and forecast typhoon path information released by the meteorological center; and observing the sea surface wind field observation data through a satellite scatterometer, the method further includes:

[0023] Identify the surrounding area centered on the current forecast typhoon center;

[0024] Calculate the spatial coverage rate between the sea surface wind field observation data and the surrounding area;

[0025] When the spatial coverage rate is greater than 50%, determine that the sea surface wind field observation data observes a typhoon, and trigger the execution of the step of detecting whether the current forecast typhoon center matching the current forecast typhoon center is included in the sea surface wind field observation data.

[0026] Further, the method further includes:

[0027] When the current forecast typhoon center is not included in the sea surface wind field observation data, extract two forecast times close to the current observation time, and extract two typhoon forecast positions corresponding to the two forecast times;

[0028] Interpolate the two typhoon forecast positions to obtain an interpolated forecast typhoon center, and perform the step of dividing the region based on the preset region division criteria and the interpolated forecast typhoon center to obtain a typhoon center positioning region including a plurality of grid nodes.

[0029] In a second aspect of the embodiments of the present application, a typhoon center automatic positioning device is provided. The typhoon center automatic positioning device includes: an acquisition unit, configured to acquire the current forecast typhoon center and forecast typhoon path information released by a meteorological center; and observe the sea surface wind field observation data through a satellite scatterometer.

[0030] A detection unit, configured to detect whether the sea surface wind field observation data includes a current forecast typhoon center that matches the current forecast typhoon center.

[0031] An extraction unit, configured to extract the current observation time when the current forecast typhoon center is detected when the sea surface wind field observation data includes the current forecast typhoon center.

[0032] An interpolation unit, configured to perform interpolation in the forecast typhoon path information to obtain an interpolated forecast typhoon center corresponding to the current observation time.

[0033] A division unit, configured to divide the region based on the preset region division criteria and the interpolated forecast typhoon center to obtain a typhoon center positioning region including a plurality of grid nodes.

[0034] A determination unit, configured to determine a positioning grid node in the typhoon center positioning region; the simulated wind direction data corresponding to the positioning grid node is most similar to the observed wind direction data.

[0035] The determination unit is further configured to determine the positioning grid node as the positioning typhoon center and extract the longitude and latitude information of the positioning typhoon center.

[0036] Further, the determination unit includes:

[0037] A modeling subunit, configured to determine the positions of a plurality of the grid nodes as a plurality of guessed typhoon centers.

[0038] A simulation subunit, configured to perform calculation and simulation according to the holland model and each of the guessed typhoon centers to obtain multiple sets of simulated wind direction data.

[0039] An acquisition subunit, configured to acquire the observed wind direction data in the sea surface wind field observation data.

[0040] A determination subunit, configured to determine a set of simulated wind direction data that is most similar to the observed wind direction data among the multiple sets of simulated wind direction data as the similar wind direction data.

[0041] The determining subunit is further configured to determine the grid nodes corresponding to the similar wind direction data as the positioning grid nodes corresponding to the interpolated predicted typhoon center.

[0042] Further, the determining subunit includes:

[0043] A calculation module, configured to calculate according to each group of the simulated wind direction data and the observed wind direction data to obtain a plurality of average absolute values of wind direction differences;

[0044] A determining module, configured to extract the minimum average absolute value of wind direction differences from the plurality of average absolute values of wind direction differences, and determine a group of simulated wind direction data corresponding to the minimum average absolute value of wind direction differences as the similar wind direction data.

[0045] A third aspect of the embodiments of the present application provides an electronic device, including a memory and a processor, where the memory is used to store a computer program, and the processor runs the computer program to enable the electronic device to execute the typhoon center automatic positioning method according to any one of the first aspects of the embodiments of the present application.

[0046] A fourth aspect of the embodiments of the present application provides a computer-readable storage medium, which stores computer program instructions. When the computer program instructions are read and run by a processor, the typhoon center automatic positioning method according to any one of the first aspects of the embodiments of the present application is executed. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings required to be used in the embodiments of the present application. It should be understood that the following drawings only show some embodiments of the present application, and therefore should not be regarded as limiting the scope. For those of ordinary skill in the art, other related drawings can be obtained based on these drawings without creative efforts.

[0048] Figure 1 A flowchart of a typhoon center automatic positioning method provided by an embodiment of the present application;

[0049] Figure 2 A structural diagram of a typhoon center automatic positioning device provided by an embodiment of the present application;

[0050] Figure 3 An effect diagram of typhoon center automatic positioning provided by an embodiment of the present application;

[0051] Figure 4 Another effect diagram of typhoon center automatic positioning provided by an embodiment of the present application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0052] The technical solutions in the embodiments of the present application will be described below in conjunction with the drawings in the embodiments of the present application.

[0053] It should be noted that similar reference numerals and letters represent similar items in the following drawings. Therefore, once an item is defined in one drawing, it does not need to be further defined or explained in subsequent drawings. At the same time, in the description of this application, the terms "first", "second", etc. are only used to distinguish the description and should not be understood as indicating or implying relative importance.

[0054] Example 1

[0055] Please see Figure 1 , Figure 1 A schematic flow chart of a method for automatically locating the center of a typhoon is provided for this embodiment. The method includes:

[0056] S101. Obtain the current forecast typhoon center and forecast typhoon path information issued by the meteorological center; observe sea surface wind field observation data through satellite scatterometers.

[0057] In this embodiment, the ocean surface wind field is an active factor influencing ocean waves, currents, and water masses, and a fundamental parameter of ocean dynamics. Monitoring global ocean wind fields is of great significance for disaster prevention and mitigation in coastal areas, protecting the marine environment, and promoting ocean-related scientific research. Satellite scatterometers, with their all-day, all-weather capabilities, high temporal and spatial resolution, and wide coverage, have become the most important means of observing global ocean surface wind fields. Specifically, a satellite scatterometer is a calibrated radar that actively transmits electromagnetic waves toward the ocean surface and receives echo signals modulated by the surface. This echo signal is determined by both the transmitted signal and the characteristics of the ocean surface. When the wavelength of the ocean waves and the wavelength of the electromagnetic waves transmitted by the radar meet the Bragg scattering condition, the backscattered electromagnetic waves generated by each wavefront have the same phase, resulting in resonance. The echo energy is primarily determined by the resonant electromagnetic wave. At the operating frequency of the microwave scatterometer, ocean surface waves that meet the Bragg resonance condition are capillary waves. The spectral density of capillary waves on the ocean surface is directly related to the wind speed on the ocean surface. Therefore, information about the ocean surface wind field can be obtained from the echo signals measured by the radar. By processing the radar echo signal, the normalized backscatter coefficient (NRCS, or σ 0 ), σ measured from the scatterometer 0 The sea surface wind field can be further extracted, and the process of extracting information about the sea surface wind field is called wind vector inversion.

[0058] In this embodiment, the method can use the OceanSat-2 satellite microwave scatterometer (HY2-SCAT). HY2-SCAT is mainly used for global sea surface wind field observation. The wind speed measurement range is 4 - 24 m / s, and the wind speed accuracy is 2 m / s or 10%; the wind direction measurement range is 0 - 360°, and the wind direction accuracy is ±20°. The working frequency of HY2-SCAT is 13.256 GHz. It adopts a pencil beam conical scanning method. The pencil beam rotates around the nadir direction at a fixed elevation angle, forming a certain ground coverage swath during the movement of the satellite platform along the track direction; the scatterometer system includes two polarization modes, VV and HH, which are observed at different incident angles respectively. During the movement of the platform, multiple backscattering coefficient (σ 0 ) measurement results with different polarization modes and different incident angles can be obtained for the same resolution cell to overcome the multi-valued ambiguity problem in the inversion of the sea surface wind field direction. Among them, the inner beam adopts the HH polarization mode with an incident angle of 41°, and the corresponding ground swath width is about 1350 km. The outer beam adopts the VV polarization mode with an incident angle of 48°, and the corresponding ground swath width is about 1700 km. The OceanSat-2 satellite scatterometer L2B product data file contains the sea surface wind speed and direction obtained through inversion, as well as information such as the corresponding observation time, longitude, and latitude. The data is organized in units of orbits, that is, the wind vector measurement data of each orbit constitutes an L2B file. Each data element in the L2B product can be indexed through the row and column numbers of the wind vector unit. The extension direction of the row of the L2B wind vector unit is perpendicular to the sub-satellite track, and the extension direction of the column is consistent with the sub-satellite track direction. Among them, the OceanSat-2 satellite scatterometer L2B product is also called the L2B product, which is used to refer to the ocean data observed by the OceanSat-2 satellite scatterometer.

[0059] In this embodiment, the method can obtain the current forecast typhoon center and forecast typhoon path information released by the National Meteorological Center at 1-hour intervals.

[0060] In this embodiment, the method can obtain the scatterometer L2B-level (along-track) sea surface wind field observation data within 3 hours of the current time.

[0061] S102. Identify the surrounding area centered on the current forecast typhoon center.

[0062] S103. Calculate the spatial coverage rate between the sea surface wind field observation data and the surrounding area.

[0063] S104. When the spatial coverage rate is greater than 50%, determine that the sea surface wind field observation data observes a typhoon, and detect whether the current forecast typhoon center included in the sea surface wind field observation data matches the current forecast typhoon center. If so, execute step S105 and steps S108 - S115; if not, execute steps S106 - S107 and steps S109 - S115.

[0064] In this embodiment, the method can calculate the spatial coverage rate of the L2B-level scatterometer wind field data within a 300 km × 300 km area centered on the position of the "current forecast typhoon center". If the spatial coverage rate exceeds 50%, it is considered that the L2B-level scatterometer wind field data observes the typhoon; if the coverage rate is lower than or equal to 50%, it is considered that the L2B-level scatterometer wind field data fails to observe the typhoon.

[0065] S105. Extract the current observation time when the current forecast typhoon center is detected.

[0066] In this embodiment, when the L2B-level wind field data product can observe the "current forecast typhoon center", the method can extract the observation time of the observation point closest to the "current forecast typhoon center" from the L2B-level wind field data product as the current observation time.

[0067] S106. Extract two forecast times close to the current observation time, and extract two typhoon forecast positions corresponding to the two forecast times. s

[0068] S107. Interpolate the two typhoon forecast positions to obtain an interpolated forecast typhoon center.

[0069] In this embodiment, when the L2B-level wind field data product fails to observe the "current forecast typhoon center", it is determined that no typhoon is observed in the track data. Therefore, this process needs to be ended and directly jump to the next track data for corresponding processing.

[0070] In this embodiment, since the time corresponding to the forecast typhoon center is usually at a 3-hour interval or a 6-hour interval, resulting in not being exactly matched with the observation time; therefore, the method makes the time of the forecast result match the observation time through interpolation. Specifically, because the typhoon center moves relatively fast, in order to avoid introducing relatively large errors, this kind of time interpolation is very necessary.

[0071] S108. Interpolate in the forecast typhoon path information to obtain an interpolated forecast typhoon center corresponding to the current observation time.

[0072] In this embodiment, through the typhoon path information released by the National Meteorological Center, the forecast position of the typhoon center corresponding to the "current observation time" is obtained by interpolation, denoted as the interpolated forecast typhoon center.

[0073] S109. Based on the preset area division standard and the interpolated forecast typhoon center, perform area division to obtain a typhoon center positioning area including multiple grid nodes.

[0074] In this embodiment, the method can select a region with a size of 5°×5° centered on the interpolated predicted typhoon center and divide it into multiple grids of 1 / 8°×1 / 8°. All the grids form the typhoon center positioning region.

[0075] S110. Determine the positions of multiple grid nodes as multiple guessed typhoon centers.

[0076] S111. Calculate and simulate according to the Holland model and each guessed typhoon center to obtain multiple sets of simulated wind direction data.

[0077] In this embodiment, the Holland typhoon model was proposed by Professor Holland in 1980 and is used to describe the two-dimensional wind field structure of a mature typhoon. This model has been successfully applied to construct the strong wind models of altimeters and scatterometers, indicating that the typhoon wind field described by this model has high accuracy and can meet the need for simulating and analyzing the typical characteristics of the sea surface backscattering coefficient under typhoon conditions.

[0078] In Holland's typhoon model, the gradient wind can be expressed as Equation 1:

[0079]

[0080] Where U g is the magnitude of the gradient wind at a distance r from the typhoon center, ρ is the air density, p0 is the central pressure, and p n is the pressure far from the typhoon center (ambient pressure), and f is the Coriolis force. The parameters A and B can be obtained through Equation 2 and Equation 3 respectively:

[0081]

[0082] B = 1.5 + (980 - p0) / 120 (Equation 3)

[0083] Where R max is in km and p0 is in mb.

[0084] Specifically, the Holland typhoon model can calculate the typhoon wind vector from the gradient wind. Assuming that the typhoon inflow angle is 25°, to obtain the 1-minute average wind vector at a height of 10 meters above the sea surface, a factor of 0.8 needs to be multiplied on the basis of the gradient wind. Due to the movement of the typhoon center, the typhoon wind field is not a symmetric structure. To account for the influence of the typhoon movement speed, the typhoon movement speed needs to be superimposed on the wind vector calculated by the model. The present invention adopts a relatively simple processing method, that is, directly linearly superimposing the typhoon movement speed on the model wind vector. The input parameters in the model include the typhoon center position, central pressure, maximum wind radius, and typhoon movement speed.

[0085] Among them, the typhoon central pressure and the typhoon moving speed can be queried through the website of the meteorological bureau; in the case where the typhoon central pressure and the typhoon moving speed cannot be obtained, the default values are taken as ρ = 1.15×10 -2 , p n = 1000 mb, p0 = 920 mb; the typhoon moving speed is determined by the historical moving speed of the current typhoon; if the historical moving speed of the typhoon cannot be obtained, the typhoon moving speed is set to 0;

[0086] Coriolis force calculation:

[0087] f = 2*(2*pi / 86400)*sin(17.1*2*pi / 360)*1000; % coriolis parameter.

[0088] S112. Obtain the observed wind direction data in the sea surface wind field observation data.

[0089] S113. Determine the group of simulated wind direction data that is most similar to the observed wind direction data among multiple groups of simulated wind direction data as the similar wind direction data.

[0090] As an optional implementation manner, determining the group of simulated wind direction data that is most similar to the observed wind direction data among multiple groups of simulated wind direction data as the similar wind direction data includes:

[0091] Calculate according to each group of simulated wind direction data and the observed wind direction data to obtain multiple mean absolute values of wind direction differences;

[0092] Extract the minimum mean absolute value of wind direction differences among the multiple mean absolute values of wind direction differences, and determine the group of simulated wind direction data corresponding to the minimum mean absolute value of wind direction differences as the similar wind direction data.

[0093] S114. Determine the grid node corresponding to the similar wind direction data as the positioning grid node corresponding to the interpolated forecast typhoon center.

[0094] S115. Determine the positioning grid node as the positioning typhoon center, and extract the longitude and latitude information of the positioning typhoon center.

[0095] In this embodiment, the method can take the position of each grid node as the guessed typhoon center, combine with Holland's typhoon model, simulate and interpolate the typhoon wind field information around the forecast typhoon center, and compare it with the scatterometer observation results, calculate the absolute value of the difference between the simulated typhoon wind direction and the scatterometer observed wind direction, sum up the absolute values, and take the mean, which is denoted as the mean absolute value of wind direction differences. Then, compare the mean absolute value of wind direction differences corresponding to each grid node, and take the position of the grid node with the minimum mean absolute value of wind direction differences as the final typhoon center positioning position (i.e., the positioning typhoon center).

[0096] In this embodiment, the method can interpolate the approximate position information of the typhoon center corresponding to the current scatterometer observation time through the typhoon track forecast information released by the National Meteorological Center. In the 5°×5° spatial range near the center position, each point is selected as a guessed point of the typhoon center at a resolution of 0.125°×0.125°. Through Holland's typhoon model, the typhoon wind direction centered on the guessed point is calculated and compared with the typhoon wind direction of the scatterometer; by comparing the comparison results corresponding to each guessed point, the guessed point closest to the scatterometer observation result is used as the final typhoon center positioning result.

[0097] Please refer to Figure 3 and Figure 4 , Figure 3 and Figure 4 show the automatic positioning effect diagrams of two typhoon centers. In the figure, "+" represents the position of the automatically positioned typhoon center. It can be seen from the figure that even if the typhoon center fails to obtain effective observation results due to being affected by land, the typhoon center position can still be obtained through this method.

[0098] In this embodiment, the execution subject of the method can be a computing device such as a computer or a server, and no specific limitation is made in this embodiment.

[0099] It should be understood that this method takes the automatic positioning of the typhoon center observed by the HY-2 satellite microwave scatterometer as an example. Aiming at the disadvantages of the traditional typhoon center automatic positioning method, such as low algorithm efficiency, inability to apply to complex situations where the typhoon is affected by land and only part of the typhoon wind field is observed, and being limited by the observation resolution of the scatterometer itself, a scatterometer typhoon center automatic positioning method based on typhoon prior knowledge is proposed. This method can be used for wind field inversion development of the HY-2 satellite scatterometer and is also applicable to the typhoon center positioning process of other satellite scatterometers. Therefore, this method has universality.

[0100] It can be seen that implementing the typhoon center automatic positioning method described in this embodiment can achieve the automatic positioning of the typhoon center when the observation results are affected by land or only part of the typhoon is observed (especially when the typhoon center is not observed); at the same time, this method can also break through the limitation of the scatterometer observation resolution and achieve the effect of super-resolution typhoon center recognition; in addition, this method can avoid the subjective errors caused by manual interpretation and reduce the workload of the corresponding duty staff; finally, this method only needs the data of the satellite scatterometer itself to optimize the typhoon wind direction without other auxiliary data.

[0101] Embodiment 2

[0102] Please refer to Figure 2 , Figure 2 which is a schematic structural diagram of a typhoon center automatic positioning device provided in this embodiment. AsFigure 2 As shown in Figure 2 , the typhoon center automatic positioning device includes:

[0103] An acquisition unit 210, configured to acquire the current forecast typhoon center and forecast typhoon path information released by the meteorological center; and observe the sea surface wind field observation data through a satellite scatterometer;

[0104] A detection unit 220, configured to detect whether the sea surface wind field observation data includes the current forecast typhoon center that matches the current forecast typhoon center;

[0105] An extraction unit 230, configured to extract the current observation time when the current forecast typhoon center is detected in the sea surface wind field observation data;

[0106] An interpolation unit 240, configured to perform interpolation in the forecast typhoon path information to obtain an interpolated forecast typhoon center corresponding to the current observation time;

[0107] A division unit 250, configured to perform regional division based on a preset regional division criterion and the interpolated forecast typhoon center to obtain a typhoon center positioning region including a plurality of grid nodes;

[0108] A determination unit 260, configured to determine a positioning grid node in the typhoon center positioning region; the simulated wind direction data corresponding to the positioning grid node is most similar to the observed wind direction data;

[0109] The determination unit 260 is further configured to determine the positioning grid node as the positioning typhoon center and extract the longitude and latitude information of the positioning typhoon center.

[0110] As an optional implementation manner, the determination unit 260 includes:

[0111] A modeling subunit 261, configured to determine the positions of a plurality of grid nodes as a plurality of guessed typhoon centers;

[0112] A simulation subunit 262, configured to perform calculation and simulation according to the Holland model and each guessed typhoon center to obtain multiple groups of simulated wind direction data;

[0113] An acquisition subunit 263, configured to acquire the observed wind direction data in the sea surface wind field observation data;

[0114] A determination subunit 264, configured to determine a group of simulated wind direction data that is most similar to the observed wind direction data among the multiple groups of simulated wind direction data as the similar wind direction data;

[0115] The determination subunit 264 is further configured to determine the grid node corresponding to the similar wind direction data as the positioning grid node corresponding to the interpolated forecast typhoon center.

[0116] As an alternative embodiment, the determination subunit 264 includes:

[0117] A calculation module, configured to calculate based on each set of simulated wind direction data and observed wind direction data to obtain multiple average absolute values of wind direction differences;

[0118] A determination module, configured to extract the minimum average absolute value of wind direction differences from the multiple average absolute values of wind direction differences, and determine a set of simulated wind direction data corresponding to the minimum average absolute value of wind direction differences as similar wind direction data.

[0119] As an alternative embodiment, the typhoon center automatic positioning device further includes:

[0120] An identification unit 270, configured to identify the surrounding area centered on the current predicted typhoon center;

[0121] A calculation unit 280, configured to calculate the spatial coverage rate between the sea surface wind field observation data and the surrounding area;

[0122] The determination unit 260 is further configured to, when the spatial coverage rate is greater than 50%, determine that the sea surface wind field observation data observes a typhoon, and trigger the detection unit 220 to perform an operation of detecting whether the sea surface wind field observation data includes the current predicted typhoon center that matches the current predicted typhoon center.

[0123] As an alternative embodiment, when the sea surface wind field observation data does not include the current predicted typhoon center, the extraction unit 230 is further configured to extract two prediction times close to the current observation time, and extract two typhoon prediction positions corresponding to the two prediction times; perform interpolation on the two typhoon prediction positions to obtain an interpolated predicted typhoon center, and perform area division based on the preset area division standard and the interpolated predicted typhoon center to obtain a typhoon center positioning area including multiple grid nodes.

[0124] In the embodiments of the present application, the explanation of the typhoon center automatic positioning device may refer to the description in Embodiment 1, and will not be elaborated herein.

[0125] It can be seen that implementing the typhoon center automatic positioning device described in this embodiment can achieve automatic positioning of the typhoon center when the observation results are affected by land or only partially observed (especially when the typhoon center is not observed); at the same time, the method can also break through the limitation of the observation resolution of the scatterometer and achieve the effect of super-resolution typhoon center identification; in addition, the method can also avoid the subjective errors caused by manual interpretation and reduce the workload of the corresponding duty staff; finally, the method only needs the data of the satellite scatterometer itself to optimize the typhoon wind direction without other auxiliary data.

[0126] An embodiment of the present application provides an electronic device, including a memory and a processor. The memory is used to store a computer program, and the processor runs the computer program to enable the electronic device to execute the typhoon center automatic positioning method in Embodiment 1 of the present application.

[0127] An embodiment of the present application provides a computer-readable storage medium, which stores computer program instructions. When the computer program instructions are read and run by a processor, the typhoon center automatic positioning method in Embodiment 1 of the present application is executed.

[0128] In several embodiments provided by the present application, it should be understood that the disclosed devices and methods can also be implemented in other ways. The device embodiments described above are only illustrative. For example, the flowcharts and block diagrams in the accompanying drawings show the possible architectures, functions, and operations of devices, methods, and computer program products according to multiple embodiments of the present application. In this regard, each block in the flowchart or block diagram may represent a module, a program segment, or a part of code, and the module, program segment, or part of code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the blocks may occur in a different order from that marked in the accompanying drawings. For example, two consecutive blocks may actually be executed substantially in parallel, and they may sometimes be executed in the reverse order, depending on the functions involved. It should also be noted that each block in the block diagram and / or flowchart, and the combination of blocks in the block diagram and / or flowchart, can be implemented by a dedicated hardware-based system for performing the specified functions or actions, or can be implemented by a combination of dedicated hardware and computer instructions.

[0129] In addition, in each embodiment of the present application, the various functional modules may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.

[0130] When the above-mentioned functions are implemented in the form of software function modules 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 this application, in essence, or the part that contributes to the prior art or a part of this technical solution can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of this application. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0131] The above are only the embodiments of this application and are not used to limit the protection scope of this application. For those skilled in the art, this application can have various changes and modifications. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of this application shall be included in the protection scope of this application. It should be noted that similar reference numerals and letters denote similar items in the following drawings. Therefore, once an item is defined in one drawing, it does not need to be further defined and explained in subsequent drawings.

[0132] As mentioned above, these are only the specific implementation manners of this application, but the protection scope of this application is not limited thereto. Any person skilled in the art can easily think of changes or replacements within the technical scope disclosed in this application, and all should be covered by the protection scope of this application. Therefore, the protection scope of this application should be subject to the protection scope of the claims.

[0133] It should be noted that in this article, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the term "comprising", "including" or any other variant thereof is intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements not only includes those elements but also includes other elements not expressly listed, or further includes elements inherent to such process, method, article or device. Without further limitation, an element defined by the phrase "including a..." does not exclude the existence of additional identical elements in the process, method, article or device including the said element.

Claims

1. An automatic typhoon center positioning method, characterized in that, Including: Obtain the current forecast typhoon center and forecast typhoon path information released by the meteorological center; observe the sea surface wind field observation data through a satellite scatterometer; Detect whether the current forecast typhoon center matching the current forecast typhoon center is included in the sea surface wind field observation data; When the current forecast typhoon center is included in the sea surface wind field observation data, extract the current observation time when the current forecast typhoon center is detected; Interpolate in the forecast typhoon path information to obtain an interpolated forecast typhoon center corresponding to the current observation time; Based on a preset regional division criterion and the interpolated forecast typhoon center, perform regional division to obtain a typhoon center positioning region including a plurality of grid nodes; Determine the positions of the plurality of grid nodes as a plurality of guessed typhoon centers; According to the Holland model and each of the guessed typhoon centers, perform calculation and simulation to obtain multiple groups of simulated wind direction data; Obtain the observed wind direction data in the sea surface wind field observation data; Determine a group of simulated wind direction data most similar to the observed wind direction data among the multiple groups of simulated wind direction data as the similar wind direction data; Determine the grid node corresponding to the similar wind direction data as the positioning grid node corresponding to the interpolated forecast typhoon center; the simulated wind direction data corresponding to the positioning grid node is the most similar to the observed wind direction data; Determine the positioning grid node as the positioning typhoon center, and extract the longitude and latitude information of the positioning typhoon center.

2. The typhoon center automatic positioning method according to claim 1, characterized in that The step of determining a group of simulated wind direction data most similar to the observed wind direction data among the multiple groups of simulated wind direction data as the similar wind direction data includes: Perform calculation according to each group of the simulated wind direction data and the observed wind direction data to obtain multiple mean absolute values of wind direction differences; Extract the minimum mean absolute value of wind direction differences among the multiple mean absolute values of wind direction differences, and determine the group of simulated wind direction data corresponding to the minimum mean absolute value of wind direction differences as the similar wind direction data.

3. The typhoon center automatic positioning method according to claim 1, wherein, After the step of obtaining the current forecast typhoon center and forecast typhoon path information released by the meteorological center; observing the sea surface wind field observation data through a satellite scatterometer, the method further includes: Identify the surrounding area centered on the current forecast typhoon center; Calculate the spatial coverage rate between the sea surface wind field observation data and the surrounding area; When the spatial coverage rate is greater than 50%, determine that the sea surface wind field observation data observes a typhoon, and trigger the step of detecting whether the current forecast typhoon center matching the current forecast typhoon center is included in the sea surface wind field observation data.

4. The typhoon center automatic positioning method according to claim 1, wherein The method further includes: When the current forecast typhoon center is not included in the sea surface wind field observation data, extract two forecast times close to the current observation time, and extract two typhoon forecast positions corresponding to the two forecast times; Interpolate the two typhoon forecast positions to obtain an interpolated forecast typhoon center, and perform the step of performing regional division based on a preset regional division criterion and the interpolated forecast typhoon center to obtain a typhoon center positioning region including a plurality of grid nodes.

5. An automatic typhoon center positioning device, characterized in that, The typhoon center automatic positioning device includes: An acquisition unit for acquiring the current forecast typhoon center and forecast typhoon path information released by a meteorological center; and observing sea surface wind field observation data through a satellite scatterometer; A detection unit for detecting whether the sea surface wind field observation data includes the current forecast typhoon center that matches the current forecast typhoon center; An extraction unit for extracting the current observation time when the current forecast typhoon center is detected in the sea surface wind field observation data; An interpolation unit for interpolating in the forecast typhoon path information to obtain an interpolated forecast typhoon center corresponding to the current observation time; A division unit for dividing regions based on a preset region division criterion and the interpolated forecast typhoon center to obtain a typhoon center positioning region including a plurality of grid nodes; A determination unit for determining a positioning grid node in the typhoon center positioning region; the simulated wind direction data corresponding to the positioning grid node is most similar to the observed wind direction data; The determination unit is further configured to determine the positioning grid node as the positioning typhoon center and extract the longitude and latitude information of the positioning typhoon center; Wherein, the determination unit includes: A modeling subunit for determining the positions of a plurality of the grid nodes as a plurality of guessed typhoon centers; A simulation subunit for performing calculation and simulation according to the Holland model and each of the guessed typhoon centers to obtain multiple groups of simulated wind direction data; An acquisition subunit for acquiring the observed wind direction data in the sea surface wind field observation data; A determination subunit for determining a group of simulated wind direction data that is most similar to the observed wind direction data among the multiple groups of simulated wind direction data as the similar wind direction data; The determination subunit is further configured to determine the grid node corresponding to the similar wind direction data as the positioning grid node corresponding to the interpolated forecast typhoon center.

6. The typhoon center automatic positioning device according to claim 5, characterized in that, The determination subunit includes: A calculation module for calculating according to each group of the simulated wind direction data and the observed wind direction data to obtain a plurality of mean absolute values of wind direction differences; A determination module for extracting the minimum mean absolute value of wind direction differences among the plurality of mean absolute values of wind direction differences and determining a group of simulated wind direction data corresponding to the minimum mean absolute value of wind direction differences as the similar wind direction data.

7. An electronic device, characterized in that, The electronic device includes a memory and a processor, the memory is used for storing a computer program, and the processor runs the computer program to enable the electronic device to execute the typhoon center automatic positioning method according to any one of claims 1 to 4.

8. A readable storage medium, characterized in that, Computer program instructions are stored in the readable storage medium, and when the computer program instructions are read and run by a processor, the typhoon center automatic positioning method according to any one of claims 1 to 4 is executed.

Images