Mesoscale vortex identification method and system, electronic equipment and storage medium

Abstract

The invention provides a mesoscale vortex recognition method and system, electronic equipment and a storage medium, and belongs to the field of vortex recognition. The method comprises the following steps: acquiring absolute dynamic topographic data, and removing large-scale signals by adopting high-pass filtering to obtain target absolute dynamic topographic data; generating and drawing a closed isoline based on the data, and extracting vortex seed points based on local extreme points of the data; the isolines are traversed, and isolines meeting vortex boundary and potential vortex boundary screening conditions are extracted from the isolines to serve as vortex boundaries and potential vortex boundaries; when the isoline traversal is completed, obtaining mesoscale vortex data; wherein a multi-thread GPU parallel technology is introduced in the process of extracting vortex seed points and vortex boundaries. The problems of data code redundancy, long operation time and inaccuracy of identification results in a traditional algorithm are solved by introducing a multi-thread GPU parallel technology and efficient identification of the mesoscale vortexes, rapid and accurate identification of the mesoscale vortexes is realized, and the identification efficiency is improved.

Description

Technical Field

[0001] The present invention belongs to the technical field of vortex identification, and specifically relates to a method, system, electronic equipment and storage medium for identifying a mesoscale vortex. Background Art

[0002] In order to study the shape, size, and distribution of global ocean vortices, a variety of vortex automatic identification algorithms have been developed. These algorithms are mainly divided into the Euler method and the Lagrangian method. The Euler method can be further divided into methods based on physical parameters, methods based on flow field geometric characteristics, and methods based on sea level anomalies (SLA). The most commonly used vortex identification method is the SLA-based method, which takes the SLA local extreme value as the vortex center and the outermost closed SLA contour line as the vortex boundary, and adds screening conditions such as shape, size, and amplitude to the vortex boundary.

[0003] However, the existing SLA-based vortex identification algorithm has the following disadvantages: (1) Data code redundancy, repeated identification of extreme points and traversal of contour lines, repeated input and output of the same data, and not concise enough; (2) The operation takes a long time. Although the global block is divided into parallel operations, the identification efficiency is still not high; (3) The identification results are not accurate enough. Summary of the invention

[0004] The purpose of the embodiments of the present invention is to provide a method, system, electronic device and storage medium for identifying a mesoscale vortex, which are used to completely or at least partially solve the technical problems of data code redundancy, long operation time and inaccurate identification results existing in the above-mentioned prior art.

[0005] In a first aspect, an embodiment of the present application provides a method for identifying a mesoscale vortex, comprising: Obtain absolute dynamic terrain data and geostrophic flow velocity data, and use high-pass filtering to remove large-scale signals in the absolute dynamic terrain data to obtain target absolute dynamic terrain data; Generate and draw a plurality of closed absolute dynamic terrain contour lines based on the target absolute dynamic terrain data, and extract vortex seed points based on local extreme value points of the target absolute dynamic terrain data; Traversing a plurality of closed absolute dynamic terrain contour lines and extracting the contour lines satisfying the vortex boundary screening conditions as vortex boundaries, and extracting the contour lines satisfying the potential vortex boundary screening conditions as potential vortex boundaries; When the absolute dynamic terrain contour traversal is completed, the mesoscale eddy data are obtained; Among them, multi-threaded GPU parallel technology was introduced in the process of extracting vortex seed points and extracting vortex boundaries.

[0006] Optionally, generating and drawing a plurality of closed absolute dynamic terrain contour lines based on the target absolute dynamic terrain data and extracting vortex seed points based on local extreme value points of the target absolute dynamic terrain data include: Generate and draw multiple closed absolute dynamic terrain contour lines at intervals of preset length; The absolute dynamic terrain value of each pixel within the closed absolute dynamic terrain contour is compared with the absolute dynamic terrain value of the grid pixel in its domain. If the absolute dynamic terrain value of the pixel is less than the absolute dynamic terrain value of the grid pixel in its domain, the pixel is a local minimum point and is identified as the seed point of the candidate cyclonic vortex; otherwise, it is identified as the seed point of the candidate anticyclonic vortex.

[0007] Optionally, multiple closed absolute dynamic terrain contour lines are traversed and the contour lines that meet the vortex boundary screening conditions are extracted as vortex boundaries, and the contour lines that meet the potential vortex boundary screening conditions are extracted as potential vortex boundaries. The traversal range is from the minimum value of all candidate cyclonic vortex seed points to the maximum value of all candidate anticyclonic vortex seed points.

[0008] Optionally, the vortex boundary screening conditions are characterized as: there is only one seed point in the closed absolute dynamic terrain contour line, the number of pixels in the closed absolute dynamic terrain contour line is between a preset range, the pixel values ​​contained in the closed absolute dynamic terrain contour line are all valid values, and the shape test error of the closed absolute dynamic terrain contour line is less than or equal to a preset error threshold and the amplitude is greater than or equal to a preset amplitude threshold.

[0009] Optionally, traversing multiple closed absolute dynamic terrain contour lines and extracting contour lines that meet the vortex boundary screening conditions as vortex boundaries, including: When the closed absolute dynamic terrain contour line meets the vortex boundary screening condition, the area contained in the closed absolute dynamic terrain contour line is regarded as the vortex, the closed absolute dynamic terrain contour line is regarded as the effective boundary of the vortex, and the extreme value point contained in the closed absolute dynamic terrain contour line is regarded as the seed point of the vortex; Traversing from the effective boundary of the vortex to the vortex seed point, the absolute dynamic terrain contour line with the maximum average geostrophic flow velocity and containing no less than 2 pixels is regarded as the maximum geostrophic flow velocity boundary of the vortex, and the absolute dynamic terrain contour line closest to the vortex seed point and containing no less than 2 pixels is regarded as the innermost boundary of the vortex. The center of mass of the vortex is the center of the fitting circle of the innermost boundary of the vortex.

[0010] Optionally, the potential vortex boundary screening condition is characterized as: there is only one seed point in the closed absolute dynamic terrain contour line and all pixel values ​​contained in the closed absolute dynamic terrain contour line are valid values.

[0011] Optionally, traversing a plurality of closed absolute dynamic terrain contour lines and extracting contour lines satisfying the potential vortex boundary screening condition as the potential vortex boundary, including: When the closed absolute dynamic terrain contour line meets the potential vortex boundary screening conditions, the area contained in the closed absolute dynamic terrain contour line is regarded as the potential vortex, the closed absolute dynamic terrain contour line is regarded as the effective boundary of the potential vortex, the extreme points contained in the closed absolute dynamic terrain contour line are regarded as the seed points of the potential vortex, and the center of mass of the potential vortex is the center of the effective boundary fitting circle.

[0012] In a second aspect, the embodiment of the present application further provides a mesoscale vortex identification system, comprising: An acquisition unit is used to acquire absolute dynamic terrain data and geostrophic flow velocity data, and use high-pass filtering to remove large-scale signals in the absolute dynamic terrain data to obtain target absolute dynamic terrain data; An extraction unit, used to generate and draw a plurality of closed absolute dynamic terrain contour lines based on the target absolute dynamic terrain data and to extract vortex seed points based on local extreme value points of the target absolute dynamic terrain data; A traversal unit is used for traversing a plurality of closed absolute dynamic terrain contour lines and extracting the contour lines satisfying the vortex boundary screening conditions as the vortex boundary, and extracting the contour lines satisfying the potential vortex boundary screening conditions as the potential vortex boundary; An acquisition unit is used to obtain mesoscale eddy data when the absolute dynamic terrain contour traversal is completed; Among them, multi-threaded GPU parallel technology was introduced in the process of extracting vortex seed points and extracting vortex boundaries.

[0013] In a third aspect, an embodiment of the present application further provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the above-mentioned method for identifying a mesoscale vortex when executing the program.

[0014] In a fourth aspect, an embodiment of the present application further provides a storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-mentioned method for identifying a mesoscale vortex.

[0015] It can be seen from the above technical solutions that the present invention has the following advantages: In the mesoscale vortex identification method, system, electronic device and storage medium provided in the present application, multi-threaded GPU parallel technology is introduced in the process of extracting vortex seed points and extracting vortex boundaries, and an efficient global mesoscale vortex identification algorithm is constructed to solve the problems of data code redundancy, long running time and insufficient accuracy of identification results in traditional algorithms, thereby achieving rapid and accurate identification of global mesoscale vortices and greatly improving identification efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solution of the present invention, the accompanying drawings required for use in the description will be briefly introduced below. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.

[0017] Figure 1 A flow chart of a method for identifying a mesoscale vortex provided in an embodiment of the present invention; Figure 2 A flow chart for identifying vortex boundaries provided by an embodiment of the present invention; Figure 3 A detailed flow chart of a method for identifying a mesoscale vortex provided in an embodiment of the present invention; Figure 4 A schematic diagram of the structure of a mesoscale vortex identification system provided by an embodiment of the present invention; Figure 5 A schematic diagram of the hardware structure of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0018] In order to better understand the technical solution of this application, the terms involved in this application are now explained: Mesoscale eddies: Mesoscale eddies are rotating water bodies that are widely present in the ocean. Their spatial scales range from tens to hundreds of kilometers. They can wrap, capture and transport water flows, affecting both the exchange of energy and the distribution of substances such as salt, dissolved carbon and chlorophyll.

[0019] Potential vortex: A "vortex" that contains only one seed point but does not meet the vortex boundary screening conditions. Its amplitude, size or shape test error does not meet the requirements for vortex identification.

[0020] In the detailed description below, various embodiments of the present disclosure will be described more fully. The present disclosure may have various embodiments, and adjustments and changes may be made therein. However, it should be understood that there is no intention to limit the various embodiments of the present disclosure to the specific embodiments disclosed herein, but rather the present disclosure should be understood to cover all adjustments, equivalents and / or alternatives that fall within the spirit and scope of the various embodiments of the present disclosure.

[0021] Hereinafter, the terms "include" or "may include" that may be used in various embodiments of the present disclosure indicate the presence of the disclosed functions or operations, and do not limit the addition of one or more functions or operations. In addition, as used in various embodiments of the present disclosure, the terms "include", "have" and their cognates are intended only to indicate specific features, numbers, steps, operations, or combinations of the foregoing items, and should not be understood as first excluding the presence of one or more other features, numbers, steps, operations, or combinations of the foregoing items or the possibility of adding one or more features, numbers, steps, operations, or combinations of the foregoing items.

[0022] In various embodiments of the present disclosure, the expression "or" or "at least one of A or / and B" includes any combination or all combinations of the words listed at the same time. For example, the expression "A or B" or "at least one of A or / and B" may include A, may include B, or may include both A and B.

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

[0024] See also Figure 1 The flowchart of a method for identifying a mesoscale vortex in a specific embodiment is shown, and includes the following execution steps: Step 100: Obtain absolute dynamic terrain data and geostrophic flow velocity data, and use high-pass filtering to remove large-scale signals in the absolute dynamic terrain data to obtain target absolute dynamic terrain data.

[0025] Specifically, satellite altimeter data is used to efficiently identify global mesoscale vortices. The data uses the Absolute Dynamic Topography (ADT) data and geostrophic flow velocity data obtained by the data center integrating multiple altimetry satellites. The satellite altimeter delay data used in this embodiment integrates multiple satellites. These data use the WGS-84 coordinate system, with a time resolution of 1 day, a spatial resolution of 0.25°×0.25°, and a global spatial range.

[0026] In some embodiments, the diameter of a mesoscale vortex is usually between tens and hundreds of kilometers. The original ADT data contains large, medium and small scale signals, among which the large scale signal will interfere with the identification of the mesoscale vortex. Therefore, before vortex identification, it is necessary to remove the large scale signal. To achieve this goal, the present embodiment uses a first-order Lanczos filter with a half-power cutoff wavelength of 700 kilometers to remove the large scale signal in the ocean, thereby smoothing the ADT field and avoiding noise interference.

[0027] Step 101: Generate and draw a plurality of closed absolute dynamic terrain contour lines based on the target absolute dynamic terrain data, and extract vortex seed points based on the local extreme value points of the target absolute dynamic terrain data.

[0028] Wherein, multi-threaded GPU parallel technology is introduced when executing step 101.

[0029] For example, by generating and drawing multiple closed ADT contour lines at intervals of 0.2 cm, the size and shape of the vortex can be more accurately extracted, thereby improving the accuracy of the vortex identification result. Since the ADT contour lines have been preliminarily screened in the process of generating them, and only the closed ADT contour lines are retained, there is no need to determine whether the contour lines are closed.

[0030] Step 102: traverse multiple closed absolute dynamic terrain contour lines and extract the contour lines that meet the vortex boundary screening conditions as vortex boundaries, and extract the contour lines that meet the potential vortex boundary screening conditions as potential vortex boundaries.

[0031] Wherein, multi-threaded GPU parallel technology is introduced when executing step 102.

[0032] Exemplarily, the vortex boundary screening condition is characterized as follows: there is only one seed point in the closed absolute dynamic terrain contour line, the number of pixels in the closed absolute dynamic terrain contour line is between the preset range, the pixel values ​​contained in the closed absolute dynamic terrain contour line are all valid values, and the shape test error of the closed absolute dynamic terrain contour line is less than or equal to the preset error threshold and the amplitude is greater than or equal to the preset amplitude threshold. The potential vortex boundary screening condition is characterized as follows: there is only one seed point in the closed absolute dynamic terrain contour line and the pixel values ​​contained in the closed absolute dynamic terrain contour line are all valid values.

[0033] It should be noted that the preset range of pixel numbers, the preset error threshold and the preset amplitude threshold can be set according to the specific application scenario and are not limited here. For example, the preset range of pixel numbers is 5-1000, and the preset error threshold is 70%. The smaller the shape test error, the closer it is to a circle. The preset amplitude threshold is 0.4 cm, where all pixel values ​​are valid values, which can be understood as not including invalid values ​​on land, and the amplitude can be understood as the absolute value of the difference between the ADT value of the closed contour line and the ADT value of the seed point.

[0034] Step 103: When the absolute dynamic terrain contour line traversal is completed, the mesoscale vortex data is obtained.

[0035] In this embodiment, in view of the inefficiency of the traditional vortex identification algorithm, the multi-threaded GPU parallel technology is innovatively introduced to construct a global mesoscale vortex efficient identification algorithm, which solves the problems of data code redundancy, long operation time and insufficient accuracy of identification results in the traditional algorithm, realizes the rapid and accurate identification of global mesoscale vortices, and constructs the corresponding data set. By combining satellite altimeter remote sensing data with advanced computer multi-threaded GPU parallel processing technology, the identification task of long-time series global mesoscale vortices can be completed in a very short time. This not only greatly improves the recognition efficiency, but also provides strong technical support for multiple key areas such as ocean dynamics research, climate change prediction, marine ecosystem protection, marine resource development and military applications.

[0036] In an embodiment of the present invention, based on step 101, a possible embodiment is given below to illustrate its specific implementation scheme in a non-limiting manner.

[0037] Generate and draw multiple closed absolute dynamic terrain contour lines at intervals of preset lengths; compare the absolute dynamic terrain value of each pixel in the closed absolute dynamic terrain contour line with the absolute dynamic terrain value of the grid pixel in its domain. If the absolute dynamic terrain value of the pixel is less than the absolute dynamic terrain value of the grid pixel in its domain, the pixel is a local minimum point and is identified as the seed point of the candidate cyclonic vortex; otherwise, it is identified as the seed point of the candidate anticyclonic vortex.

[0038] It should be understood that multiple closed absolute dynamic terrain contour lines are traversed and the contour lines that meet the vortex boundary screening conditions are extracted as vortex boundaries, and the contour lines that meet the potential vortex boundary screening conditions are extracted as potential vortex boundaries. The traversal range is from the minimum value of all candidate cyclonic vortex seed points to the maximum value of all candidate anticyclonic vortex seed points.

[0039] In one embodiment, see Figure 2As shown, in step 102, traversing multiple closed absolute dynamic terrain contour lines and extracting contour lines that meet the vortex boundary screening conditions as vortex boundaries, the following steps can be specifically performed: S1: When the closed absolute dynamic terrain contour line meets the vortex boundary screening condition, the area contained in the closed absolute dynamic terrain contour line is regarded as a vortex, the closed absolute dynamic terrain contour line is regarded as the effective boundary of the vortex, and the extreme points contained in the closed absolute dynamic terrain contour line are regarded as seed points of the vortex.

[0040] In some embodiments, the seed point is a local extreme point based on ADT data. The method for identifying the seed point is to compare the ADT value of each pixel with the ADT values ​​of 8 grid pixels in the surrounding 3×3 neighborhood. If the ADT value of the pixel is less than the ADT values ​​of all pixels in its neighborhood, the pixel is a local minimum point and can be identified as the seed point of the candidate cyclonic vortex. Otherwise, it is the seed point of the candidate anticyclonic vortex. This embodiment only considers vortices with one seed point, so it is guaranteed that there is only one minimum or maximum point in each neighborhood.

[0041] S2: Traverse from the effective boundary of the vortex to the vortex seed point, and regard the absolute dynamic terrain contour line with the maximum average geostrophic flow velocity and containing no less than 2 pixels as the maximum geostrophic flow velocity boundary of the vortex. The absolute dynamic terrain contour line closest to the vortex seed point and containing no less than 2 pixels is regarded as the innermost boundary of the vortex. The vortex center of mass is the center of the fitting circle of the innermost boundary of the vortex.

[0042] This embodiment relaxes the screening conditions for vortex amplitude, shape test error and the number of pixels included in the effective boundary, so that vortices with smaller amplitude and area and more irregular shapes can also be identified, thereby increasing the number of vortex identifications per day and improving the accuracy of vortex identification results.

[0043] In a specific embodiment, the following steps are performed to traverse multiple closed absolute dynamic terrain contour lines and extract contour lines that meet the potential vortex boundary screening conditions as the potential vortex boundary in step 102: when the closed absolute dynamic terrain contour line meets the potential vortex boundary screening conditions, the area contained in the closed absolute dynamic terrain contour line is regarded as a potential vortex, the closed absolute dynamic terrain contour line is regarded as an effective boundary of the potential vortex, the extreme points contained in the closed absolute dynamic terrain contour line are regarded as seed points of the potential vortex, and the center of mass of the potential vortex is the center of the effective boundary fitting circle.

[0044] In order to improve the objectivity and accuracy of subsequent vortex tracking results, this embodiment includes potential vortices in the identification scope, thereby expanding the data set for vortex identification. However, potential vortices are not used for any analysis of the dynamics of vortex morphology, but only for vortex tracking.

[0045] In some embodiments, multi-threaded GPU parallel computing is introduced in the process of vortex seed point extraction and vortex boundary extraction, which greatly improves the efficiency of vortex identification. Multi-threaded GPU parallel technology refers to the technology of using the parallel processing capability of GPU (graphics processing unit) to perform multiple computing tasks simultaneously. Using Numba for multi-threaded GPU parallel computing, its working principle is: Numba converts Python code into CUDA (Compute Unified Device Architecture) kernel and executes the kernel on the GPU. When a function is marked with the @njit decorator, Numba will first analyze the input type of the function and compile the function into CUDA kernel code suitable for GPU as needed. These codes will use the parallel processing capability of GPU to accelerate numerical computing tasks. At runtime, Numba loads the optimized CUDA kernel to the GPU for execution, thereby achieving acceleration. It should be noted that to use Numba for GPU acceleration, it is necessary to ensure that the computing task is suitable for execution on the GPU and that the target GPU supports CUDA. In addition, Numba also provides support for GPU memory management and parallel computing to better utilize GPU resources for accelerated computing. For example, writing cache=True in the brackets after the @njit decorator can achieve the purpose of caching, and writing fastmath=True can perform fast mathematical operations.

[0046] This embodiment applies Numba to the processes of identifying local extreme points, calculating the geostrophic flow velocity of contour lines, encrypting and resampling the points on the contour lines, constructing fitting circles, coordinate conversion, etc., and merges and simplifies redundant codes to avoid the generation of redundant data, which ultimately greatly improves the efficiency of vortex identification. The multi-threaded GPU parallel acceleration technology reduces the complexity of the original algorithm to , where is the number of threads executed in the GPU, which is determined by the GPU performance. The time to identify vortices based on one day's ADT data has been improved from one hour to one to two minutes.

[0047] In some embodiments, the global mesoscale vortex efficient identification data set obtained using the algorithm of the present invention is stored in h5 format files in units of days. There are 10,766 files from January 1, 1993 to June 23, 2023. Compared with nc, json and pkl file formats, h5 file format can store a large amount of data in an efficient manner and effectively reduce the storage space occupied, and has high efficiency and flexibility when storing large-scale vortex data. There are two types of data in the global mesoscale vortex efficient identification data set, corresponding to vortices and "potential vortices", respectively. The vortex data mainly includes various attribute information of the vortex, such as amplitude, area, kinetic energy, vorticity, etc., while the "potential vortex" data only includes its outermost boundary, polarity, seed point, and center of mass. Because it is not a real vortex in the physical sense, it is only used for objective tracking of the vortex, and does not perform morphological and dynamic analysis on it. All attribute information is stored in array form. For example, a variable array is the attribute information of all vortices on that day. In addition, all boundary points are resampled, and the number of points on each boundary is resampled to 50, and they are saved in clockwise order starting from the point corresponding to the north direction of the vortex center.

[0048] The vortex in this dataset contains three kinds of boundaries, namely the effective boundary, the maximum geostrophic velocity boundary, and the innermost boundary. The effective boundary is the outermost boundary that determines the size of the vortex. Traversing from the effective boundary of the vortex to the vortex seed point, the ADT contour line with the maximum average geostrophic velocity and containing no less than 2 pixels is the maximum geostrophic velocity boundary of the vortex, and the ADT contour line closest to the vortex seed point and containing no less than 2 pixels is the innermost boundary of the vortex. The attribute information of the vortex is divided into three categories: the first is the attributes that can be used for vortex morphological analysis, such as amplitude and shape error; the second is the attributes that can be used for vortex kinematic analysis, such as center of mass; and the third is the attributes that can be used for vortex dynamics analysis, such as kinetic energy and vorticity.

[0049] The specific vortex variable information is as follows: amplitude: the amplitude corresponding to the effective boundary of the vortex; amplitude_s: the amplitude corresponding to the maximum geostrophic velocity boundary of the vortex; contour_lon_e: the longitude sequence corresponding to the effective boundary of the vortex; contour_lat_e: the latitude sequence corresponding to the effective boundary of the vortex; contour_lon_s: the longitude sequence corresponding to the maximum geostrophic velocity boundary of the vortex; contour_lat_s: the latitude sequence corresponding to the vortex geostrophic velocity boundary; contour_lon_i: the longitude sequence corresponding to the innermost boundary of the vortex; contour_lat_i: the latitude sequence corresponding to the innermost boundary of the vortex; height_external_contour: ADT value corresponding to the effective boundary of the vortex; height_max_speed_contour: ADT value corresponding to the maximum geostrophic velocity boundary of the vortex; height_inner_contour: ADT value corresponding to the innermost boundary of the vortex; eddy_centroid_core_lon: eddy centroid longitude; eddy_centroid_core_lat: eddy centroid latitude; effective_area: the area enclosed by the effective boundary of the vortex; speed_area: the area enclosed by the boundary of the vortex with the maximum geostrophic flow velocity; eddy_vorticity_e: the average vorticity of all pixels within the effective boundary of the vortex; eddy_vorticity_s: the average vorticity of all pixels within the boundary of the vortex maximum geostrophic velocity; eke: the average kinetic energy of all pixels within the effective boundary of the vortex; eke_s: the average kinetic energy of all pixels within the boundary of the maximum geostrophic flow velocity of the vortex; lon: the longitude of the center of the fitted circle corresponding to the boundary of the maximum geostrophic flow velocity of the vortex; lat: latitude of the center of the fitted circle corresponding to the boundary of the maximum geostrophic flow velocity of the vortex; lon_max: the longitude of the center of the fitted circle corresponding to the innermost boundary of the vortex; lat_max: the latitude of the center of the fitted circle corresponding to the innermost boundary of the vortex; nb_contour_selected: the number of all contour lines associated with a certain vortex (including but not limited to the effective boundary, maximum geostrophic flow velocity boundary and innermost boundary of the vortex); num_point_e: the number of points on the effective boundary of the vortex; num_point_s: the number of points on the boundary of the vortex maximum geostrophic flow velocity; radius_e: the radius of the fitted circle corresponding to the effective boundary of the vortex; radius_s: the radius of the fitted circle corresponding to the maximum geostrophic velocity boundary of the vortex; seeds_lon: vortex seed point longitude; seeds_lat: vortex seed point latitude; shape_error_e: shape error of the effective boundary of the vortex; shape_error_s: shape error of the maximum geostrophic velocity boundary of the vortex; time: the time the vortex exists; type_cyc: type of vortex (1 for anticyclone, -1 for cyclone); speed_average: the maximum rotation speed of the vortex; uavg_profile: The average rotation speed of the vortex.

[0050] The "potential vortex" in this dataset only contains one kind of boundary, namely the effective boundary. We do not perform morphological and dynamic analysis on the "potential vortex", which is only used in the objective tracking algorithm of the global mesoscale vortex. Therefore, the attribute information of the "potential vortex" is divided into only one category, namely the attributes used for kinematic analysis, including effective boundary, centroid, seed point and polarity.

[0051] contour_lon_e_P: The longitude sequence corresponding to the effective boundary of the “potential vortex”; contour_lat_e_P: The latitude sequence corresponding to the effective boundary of the "potential vortex"; eddy_centroid_core_lon_P: longitude of the centroid of the "potential eddy"; eddy_centroid_core_lat_P: latitude of the centroid of the “potential eddy”; seeds_lon_P: longitude of the “potential vortex” seed point; seeds_lat_P: “potential eddy” seed point latitude; type_cyc_P: Type of potential vortex (1 for anticyclonic, -1 for cyclonic).

[0052] In one embodiment, Figure 3 This is a detailed flow chart of a method for identifying a mesoscale vortex according to an embodiment of the present invention. This embodiment is further optimized and expanded on the basis of the above embodiments.

[0053] The global ADT field is obtained, and then high-pass filtering is performed to generate closed ADT contour lines. The vortex seed points are extracted to traverse the ADT contour lines to determine whether the boundary screening conditions are met. If so, vortex identification is performed. Otherwise, potential vortex identification is performed to determine whether the contour lines are traversed. If completed, the global vortex data is obtained. Otherwise, the ADT contour lines are continued to be traversed until the contour lines are traversed.

[0054] It should be understood that the order of execution of the steps in the above embodiment does not necessarily mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiment of the present invention.

[0055] like Figure 4 As shown, the following is an embodiment of the mesoscale vortex identification system provided in the embodiment of the present disclosure, which belongs to the same inventive concept as the mesoscale vortex identification method of the above-mentioned embodiments. For details not described in detail in the embodiment of the mesoscale vortex identification system, please refer to the embodiment of the above-mentioned mesoscale vortex identification method.

[0056] The acquisition unit 40 is used to acquire the absolute dynamic terrain data and the geostrophic flow velocity data, and use high-pass filtering to remove the large-scale signal in the absolute dynamic terrain data to obtain the target absolute dynamic terrain data; An extraction unit 41 is used to generate and draw a plurality of closed absolute dynamic terrain contour lines based on the target absolute dynamic terrain data and to extract vortex seed points based on local extreme value points of the target absolute dynamic terrain data; A traversal unit 42 is used to traverse a plurality of closed absolute dynamic terrain contour lines and extract the contour lines satisfying the vortex boundary screening conditions as the vortex boundary, and extract the contour lines satisfying the potential vortex boundary screening conditions as the potential vortex boundary; An obtaining unit 43 is used to obtain mesoscale vortex data when the absolute dynamic terrain contour line traversal is completed; Among them, multi-threaded GPU parallel technology was introduced in the process of extracting vortex seed points and extracting vortex boundaries.

[0057] Figure 5 It is a schematic diagram of the hardware structure of an electronic device for implementing various embodiments of the present invention.

[0058] The identification method of the mesoscale vortex provided in the embodiment of the present application can be applied to electronic devices. Those skilled in the art will appreciate that the electronic device structure involved in the embodiment of the present invention does not constitute a limitation on the electronic device, and the electronic device may include more or less components than shown, or combine certain components, or different component arrangements. In an embodiment of the present invention, the electronic device includes but is not limited to a laptop computer, a desktop computer, a workbench, a personal digital assistant, a server, a blade server, a mainframe computer, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices, and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples, and are not intended to limit the implementation of the embodiments of the present application described and / or required herein.

[0059] The electronic device may include a processor, an external memory interface, an internal memory, a universal serial bus (USB) interface, a charging management module, a power management module, a battery, a wireless communication module, an audio module, a speaker, a microphone, a sensor module, buttons, a camera, a display, and a SIM card interface, etc.

[0060] It is to be understood that the structure illustrated in the embodiments of the present application does not constitute a specific limitation on the electronic device. In other embodiments of the present application, the electronic device may include more or fewer components than shown in the figure, or combine certain components, or split certain components, or arrange the components differently. The components shown in the figure may be implemented in hardware, software, or a combination of software and hardware.

[0061] The processor may include one or more processing units, for example, the processor may include a central processing unit (CPU), an application processor (AP), a modem processor, a graphics processing unit (GPU), an image signal processor (ISP), a controller, a memory, a video codec, a digital signal processor (DSP), a baseband processor, and / or a neural-network processing unit (NPU), etc. Among them, different processing units may be independent devices or integrated into one or more processors.

[0062] The processor can be the nerve center and command center of the electronic device. The controller can generate an operation control signal according to the instruction operation code and timing signal to complete the control of fetching and executing instructions.

[0063] A memory may also be provided in the processor for storing instructions and data. In some embodiments, the memory in the processor is a cache memory. The memory may store instructions or data that the processor has just used or is cyclically used. If the processor needs to use the instruction or data again, it may be directly called from the memory. This avoids repeated access, reduces the waiting time of the processor, and thus improves system efficiency.

[0064] The external memory interface can be used to connect an external memory card, such as a MicroSD card, to expand the storage capacity of the electronic device. The external memory card communicates with the processor through the external memory interface to implement data storage functions. For example, files such as music and videos can be saved in the external memory card.

[0065] The internal memory can be used to store computer executable program codes, which include instructions. The processor executes various functional applications and data processing of the electronic device by running the instructions stored in the internal memory. The internal memory may include a program storage area and a data storage area. The internal memory may include a high-speed random access memory and may also include a non-volatile memory, such as at least one disk storage device, a flash memory device, a universal flash storage (UFS), etc.

[0066] The wireless communication function of an electronic device can be realized through an antenna, a wireless communication module, a modem processor, and a baseband processor.

[0067] The wireless communication module can provide wireless communication solutions for electronic devices, including wireless local area networks (WLAN) (such as wireless fidelity (Wi-Fi) networks), Bluetooth (BT), global navigation satellite system (GNSS), frequency modulation (FM), near field communication (NFC), infrared technology (IR), etc.

[0068] Electronic devices can implement audio functions, etc. through audio modules, speakers, receivers, microphones, headphone jacks, and application processors.

[0069] Electronic devices can achieve shooting functions through ISP, camera, video codec, GPU, display and application processor.

[0070] Electronic devices can achieve display functions through GPU, display screen and application processor.

[0071] The GPU is a microprocessor for image processing that connects the display screen and the application processor. The GPU is used to perform mathematical and geometric calculations for graphics rendering. The processor may include one or more GPUs that execute program instructions to generate or change display information.

[0072] The display screen is used to display images, videos, etc. The display screen includes a display panel.

[0073] The storage medium provided in the present application stores a program product that can implement a method for identifying mesoscale vortices.

[0074] The identification method of mesoscale vortex includes: obtaining absolute dynamic terrain data, and using high-pass filtering to remove large-scale signals in the absolute dynamic terrain data to obtain target absolute dynamic terrain data; generating and drawing multiple closed absolute dynamic terrain contour lines based on the target absolute dynamic terrain data, and extracting vortex seed points based on local extreme points of the target absolute dynamic terrain data; traversing multiple closed absolute dynamic terrain contour lines and extracting contour lines that meet vortex boundary screening conditions as vortex boundaries, and extracting contour lines that meet potential vortex boundary screening conditions as potential vortex boundaries; when the traversal of the absolute dynamic terrain contour lines is completed, obtaining mesoscale vortex data; wherein, multi-threaded GPU parallel technology is introduced in the process of extracting vortex seed points and extracting vortex boundaries.

[0075] In some possible embodiments, the subject name of the present disclosure, the method and system for identifying mesoscale vortices, can be implemented in the form of a program product, which includes a program code. When the program product is run on a terminal device, the program code is used to enable the terminal device to execute the steps described in the above "Exemplary Method" section of this specification according to various exemplary embodiments of the present disclosure.

[0076] The storage medium of the present disclosure can adopt any combination of one or more readable media. The readable medium can be a readable signal medium or a readable storage medium. The readable storage medium can be, for example, but not limited to, a system, device or device of electricity, magnetism, light, electromagnetic, infrared, or semiconductor, or any combination of the above. More specific examples (non-exhaustive list) of readable storage media include: an electrical connection with one or more wires, a portable disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above.

[0077] The above description of the disclosed embodiments enables one skilled in the art to implement or use the present invention. Various modifications to these embodiments will be apparent to one skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to the embodiments shown herein, but rather to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for identifying a mesoscale vortex, characterized in that: include: Obtain absolute dynamic terrain data, and use high-pass filtering to remove large-scale signals in the absolute dynamic terrain data to obtain target absolute dynamic terrain data; Generate and draw a plurality of closed absolute dynamic terrain contour lines based on the target absolute dynamic terrain data, and extract vortex seed points based on local extreme value points of the target absolute dynamic terrain data; Traversing a plurality of closed absolute dynamic terrain contour lines and extracting the contour lines satisfying the vortex boundary screening conditions as vortex boundaries, and extracting the contour lines satisfying the potential vortex boundary screening conditions as potential vortex boundaries; When the absolute dynamic terrain contour traversal is completed, the mesoscale eddy data are obtained; Among them, multi-threaded GPU parallel technology was introduced in the process of extracting vortex seed points and extracting vortex boundaries.

2. The method for identifying a mesoscale vortex according to claim 1, characterized in that: Generating and drawing a plurality of closed absolute dynamic terrain contour lines based on the target absolute dynamic terrain data and extracting vortex seed points based on local extreme value points of the target absolute dynamic terrain data include: Generate and draw multiple closed absolute dynamic terrain contour lines at intervals of preset length; The absolute dynamic terrain value of each pixel within the closed absolute dynamic terrain contour is compared with the absolute dynamic terrain value of the grid pixel in its domain. If the absolute dynamic terrain value of the pixel is less than the absolute dynamic terrain value of the grid pixel in its domain, the pixel is a local minimum point and is identified as the seed point of the candidate cyclonic vortex; otherwise, it is identified as the seed point of the candidate anticyclonic vortex.

3. The method for identifying a mesoscale vortex according to claim 2, characterized in that: Traverse multiple closed absolute dynamic terrain contour lines and extract the contour lines that meet the vortex boundary screening conditions as the vortex boundary, and extract the contour lines that meet the potential vortex boundary screening conditions as the potential vortex boundary. The traversal range is from the minimum value of all candidate cyclonic vortex seed points to the maximum value of all candidate anticyclonic vortex seed points.

4. The method for identifying a mesoscale vortex according to claim 1, characterized in that: The vortex boundary screening conditions are characterized as follows: there is only one seed point in the closed absolute dynamic terrain contour line, the number of pixels in the closed absolute dynamic terrain contour line is between the preset range, the pixel values ​​contained in the closed absolute dynamic terrain contour line are all valid values, and the shape test error of the closed absolute dynamic terrain contour line is less than or equal to the preset error threshold and the amplitude is greater than or equal to the preset amplitude threshold.

5. The method for identifying a mesoscale vortex according to claim 4, characterized in that: Traverse multiple closed absolute dynamic terrain contour lines and extract the contour lines that meet the vortex boundary screening conditions as the vortex boundary, including: When the closed absolute dynamic terrain contour line meets the vortex boundary screening condition, the area contained in the closed absolute dynamic terrain contour line is regarded as the vortex, the closed absolute dynamic terrain contour line is regarded as the effective boundary of the vortex, and the extreme value point contained in the closed absolute dynamic terrain contour line is regarded as the seed point of the vortex; Traversing from the effective boundary of the vortex to the vortex seed point, the absolute dynamic terrain contour line with the maximum average geostrophic flow velocity and containing no less than 2 pixels is regarded as the maximum geostrophic flow velocity boundary of the vortex, and the absolute dynamic terrain contour line closest to the vortex seed point and containing no less than 2 pixels is regarded as the innermost boundary of the vortex. The center of mass of the vortex is the center of the fitting circle of the innermost boundary of the vortex.

6. The method for identifying a mesoscale vortex according to claim 1, characterized in that: The potential vortex boundary screening condition is characterized as follows: there is only one seed point in the closed absolute dynamic terrain contour line and all the pixel values ​​contained in the closed absolute dynamic terrain contour line are valid values.

7. The method for identifying a mesoscale vortex according to claim 6, characterized in that: Traverse multiple closed absolute dynamic terrain contour lines and extract the contour lines that meet the potential vortex boundary screening conditions as the potential vortex boundary, including: When the closed absolute dynamic terrain contour line meets the potential vortex boundary screening conditions, the area contained in the closed absolute dynamic terrain contour line is regarded as the potential vortex, the closed absolute dynamic terrain contour line is regarded as the effective boundary of the potential vortex, the extreme points contained in the closed absolute dynamic terrain contour line are regarded as the seed points of the potential vortex, and the center of mass of the potential vortex is the center of the effective boundary fitting circle.

8. A mesoscale vortex identification system, characterized in that: include: An acquisition unit is used to acquire absolute dynamic terrain data and remove large-scale signals in the absolute dynamic terrain data by high-pass filtering to obtain target absolute dynamic terrain data; An extraction unit, used to generate and draw a plurality of closed absolute dynamic terrain contour lines based on the target absolute dynamic terrain data and to extract vortex seed points based on local extreme value points of the target absolute dynamic terrain data; A traversal unit is used for traversing a plurality of closed absolute dynamic terrain contour lines and extracting the contour lines satisfying the vortex boundary screening conditions as the vortex boundary, and extracting the contour lines satisfying the potential vortex boundary screening conditions as the potential vortex boundary; An acquisition unit is used to obtain mesoscale eddy data when the absolute dynamic terrain contour traversal is completed; Among them, multi-threaded GPU parallel technology was introduced in the process of extracting vortex seed points and extracting vortex boundaries.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the steps of the method for identifying a mesoscale vortex as described in any one of claims 1 to 7 are implemented.

10. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for identifying a mesoscale vortex as claimed in any one of claims 1 to 7 are implemented.

Images