Methods, apparatus, computer equipment, and storage media for identifying mesoscale vortices

Abstract

This invention provides a method, apparatus, computer equipment, and storage medium for identifying mesoscale eddies, belonging to the field of eddy recognition technology. The method includes: acquiring SLA data; dividing the SLA data into multiple initial SLA splines based on multiple SLA extreme points; filtering the multiple initial SLA splines according to preset screening conditions to obtain multiple first quasi-eddy SLA splines; selecting those belonging to the mesoscale eddy from the multiple first quasi-eddy SLA splines to obtain multiple second quasi-eddy SLA splines; determining the eddy boundary and eddy center of the mesoscale eddy based on the multiple second quasi-eddy SLA splines; and determining basic attribute features based on the eddy boundary and eddy center of the mesoscale eddy. This invention only processes SLA data, enabling rapid and efficient mesoscale eddy identification, providing timely and efficient service support for marine field observation and maritime navigation.

Description

Technical Field

[0001] This invention belongs to the field of vortex recognition technology, and more specifically, relates to a method, apparatus, computer equipment, and storage medium for recognizing mesoscale vortices. Background Technology

[0002] Ocean mesoscale eddies can be classified into anticyclonic eddies and cyclonic eddies based on the direction of rotation. In the Northern Hemisphere, anticyclonic eddies rotate clockwise, while cyclonic eddies rotate counterclockwise; conversely, in the Southern Hemisphere, anticyclonic eddies rotate counterclockwise, while cyclonic eddies rotate clockwise. Since the 1970s, oceanographers have observed mesoscale eddies in the ocean and have conducted research including mesoscale eddy observations, ocean numerical simulations, and theoretical analyses. The advent of satellite altimeters has enabled the statistical study of ocean mesoscale eddies on a global scale.

[0003] Currently, satellite altimeter data acquisition is highly timely (delay not exceeding one day). Traditional identification methods, however, require acquiring 40 days (20 days before and after) of along-orbit SLA (Sea Level Anomaly) data, necessitating spatiotemporal interpolation and fusion before mesoscale eddy identification and extraction. The spatiotemporal fusion of along-orbit altimeter data involves processes such as satellite altimeter reference frame unification, orbital error correction, long-wavelength error correction, and spatiotemporal objective gridding. These data processing steps are cumbersome and complex, requiring extensive computation, consuming significant resources, and resulting in lengthy processing times. This approach fails to meet the current demands for rapid and efficient marine monitoring, maritime service support, and marine research. Furthermore, the reliance on large computing platforms hinders the development of lightweight equipment required by maritime platforms, thus failing to provide timely and efficient support for maritime navigation. Summary of the Invention

[0004] To overcome the above shortcomings, the present invention provides a method, apparatus, computer equipment, and storage medium for identifying mesoscale vortices, aiming to solve the problem that existing identification methods require data fusion, resulting in large data processing volume and long processing time.

[0005] The present invention is implemented as follows:

[0006] Firstly, a method for identifying mesoscale vortices, the method comprising:

[0007] Acquire SLA data, which includes multiple SLA extreme points;

[0008] The SLA data is divided into multiple initial SLA splines based on multiple SLA extreme points. The multiple initial SLA splines are then filtered according to preset filtering conditions to obtain multiple first quasi-vortex SLA splines.

[0009] Multiple second quasi-vortex SLA splines are obtained by filtering out those belonging to the mesoscale vortex from multiple first quasi-vortex SLA splines;

[0010] The vortex boundary and vortex center of the mesoscale vortex are determined based on multiple second quasi-vortex SLA splines;

[0011] The basic property characteristics are determined based on the vortex boundary and vortex center of the mesoscale vortex.

[0012] Furthermore, the preset filtering conditions are:

[0013] The initial SLA spline data points are continuous and without missing points;

[0014] The distance between the center extreme point and the two endpoint extreme points of the initial SLA spline exceeds the first preset threshold.

[0015] The absolute value of the first SLA difference between the central extreme point and the two endpoint extreme points of the initial SLA spline exceeds the second preset threshold.

[0016] Among the data points of adjacent initial SLA splines, the data points of the initial SLA spline whose SLA gradient exceeds a third preset threshold.

[0017] Further, the step of selecting from multiple first quasi-vortex SLA splines that belong to the mesoscale vortex to obtain multiple second quasi-vortex SLA splines includes:

[0018] Using the center point of each of the first quasi-vortex SLA splines as the search center, search for the same type of first quasi-vortex SLA splines within a preset search radius;

[0019] Compare the magnitudes of the central extreme points of the first quasi-vortex SLA splines of the same type, and determine the quasi-vortex center SLA spline and quasi-vortex center point based on the characteristics of air vortex and anti-air vortex.

[0020] Calculate the three-dimensional distance between each of the quasi-vortex center points and the spline center points of all the first quasi-vortex SLA splines;

[0021] Compare the distance relationship between the quasi-vortex center SLA spline with the three-dimensional distance being less than a preset constant and the quasi-vortex center SLA spline, and determine the first quasi-vortex SLA spline that satisfies the preset relationship condition;

[0022] Determine whether the total amount of the first quasi-vortex SLA spline that satisfies the preset relationship condition is greater than the preset total amount threshold.

[0023] If the total number of the first quasi-vortex SLA splines that satisfy the preset relationship conditions is less than or equal to the preset total threshold, then the first quasi-vortex SLA splines that satisfy the preset relationship conditions are removed.

[0024] If the total number of the first quasi-vortex SLA splines that satisfy the preset relationship conditions is greater than the preset total threshold, then the first quasi-vortex SLA splines that satisfy the preset relationship conditions are determined as the second quasi-vortex SLA splines.

[0025] Furthermore, the preset relationship conditions are as follows:

[0026] The closer the center point of the first quasi-vortex SLA spline is to the center point of the quasi-vortex, the smaller the difference between the SLA value of the center point of the first quasi-vortex SLA spline and the SLA value of the center point of the quasi-vortex.

[0027] The farther the center point of the first quasi-vortex SLA spline is from the center point of the quasi-vortex, the greater the difference between the SLA value of the center point of the first quasi-vortex SLA spline and the SLA value of the center point of the quasi-vortex.

[0028] Furthermore, the formula for calculating the three-dimensional distance is:

[0029]

[0030] In the formula: ∆ D It is a three-dimensional distance;

[0031] ( x 0, y 0) represents the latitude and longitude spatial coordinates of the quasi-vortex center point;

[0032] ( x i , y i () represents the latitude and longitude spatial coordinates of the spline center point of the quasi-vortex center spline;

[0033] z 0 represents the SLA value at the quasi-vortex center point;

[0034] z i The SLA value at the center point of the quasi-vortex center spline;

[0035] C x , C y andC z They are respectively x , y and z Characteristic constants of the variable.

[0036] Further, determining the vortex boundary and vortex center of the mesoscale vortex based on multiple second quasi-vortex SLA splines includes:

[0037] Calculate the maximum gradient value of each of the second quasi-vortex SLA splines, then calculate the average value of all the maximum gradient values, and use the average value as the SLA value of the vortex boundary of the mesoscale vortex;

[0038] Calculate the SLA value of each of the second quasi-vortex SLA splines with respect to the vortex boundary at two spline boundary position points and spline normal boundary position points;

[0039] The Graham scan method was used to calculate the spline boundary points and spline normal boundary points of all the second quasi-vortex SLA splines to obtain the convex hull;

[0040] The convex hull is smoothed using the Bézier curve method to obtain the vortex boundary of the mesoscale vortex;

[0041] The centroid is calculated based on the vortex boundary, and the centroid is taken as the vortex center of the mesoscale vortex. The SLA value of the quasi-vortex center point is taken as the SLA value of the vortex center.

[0042] Furthermore, determining the basic attribute characteristics based on the vortex boundary and vortex center of the mesoscale vortex includes:

[0043] Determine whether the SLA value at the center of the vortex is less than the SLA value at the boundary of the vortex;

[0044] If the SLA value at the center of the vortex is less than the SLA value at the boundary of the vortex, then it is an air vortex;

[0045] If the SLA value at the center of the vortex is greater than the SLA value at the boundary of the vortex, it is an anti-vortex;

[0046] The radius of the circle with the same shape and area is used as the spatial scale of the mesoscale vortex;

[0047] Calculate the absolute value of the second SLA difference between the SLA value at the center of the vortex and the SLA value at the boundary of the vortex, and use the absolute value of the second SLA difference as the intensity of the mesoscale vortex.

[0048] Secondly, a device for identifying mesoscale vortices, the device comprising:

[0049] The acquisition module is configured to acquire SLA data, which includes multiple SLA extreme points.

[0050] The first filtering module is configured to divide the SLA data into multiple initial SLA splines based on multiple SLA extreme points, and filter the multiple initial SLA splines according to preset filtering conditions to obtain multiple first quasi-vortex SLA splines.

[0051] The second filtering module is configured to filter out those belonging to the mesoscale vortex from multiple first quasi-vortex SLA splines to obtain multiple second quasi-vortex SLA splines;

[0052] The first determining module is configured to determine the vortex boundary and vortex center of the mesoscale vortex based on multiple second quasi-vortex SLA splines;

[0053] The second determining module is configured to determine basic attribute features based on the vortex boundary and vortex center of the mesoscale vortex.

[0054] Thirdly, a computer device includes a memory and a processor, the memory storing a computer program, the processor executing the computer program to implement the steps of the method described in the first aspect.

[0055] Fourthly, a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the method described in the first aspect.

[0056] Compared to existing technologies, the advantages of this invention are as follows: By dividing the along-track SLA data into multiple initial SLA splines using the SLA extreme points of the SLA data, first quasi-vortex SLA splines of the same type are selected to determine second quasi-vortex SLA splines belonging to mesoscale vortices. Next, the vortex center and shape of the mesoscale vortex are determined using the second quasi-vortex SLA splines. Finally, the basic characteristic information of the mesoscale vortex is determined based on the vortex center and vortex shape. This invention only processes SLA data and does not involve data fusion or other processes, enabling rapid and efficient mesoscale vortex identification, providing timely and efficient service support for marine field observation and maritime navigation. Attached Figure Description

[0057] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present invention and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained from these drawings without creative effort.

[0058] Figure 1 This is a flowchart of a method for identifying mesoscale vortices provided in an embodiment of the present invention;

[0059] Figure 2 This is a flowchart of step S30 provided in an embodiment of the present invention;

[0060] Figure 3 This is a flowchart of step S40 provided in an embodiment of the present invention;

[0061] Figure 4 This is a flowchart of step S50 provided in an embodiment of the present invention;

[0062] Figure 5 This is a schematic diagram of the structure of a mesoscale vortex identification device provided in an embodiment of the present invention.

[0063] Explanation of reference numerals in the attached diagram: 100, acquisition module; 200, first filtering module; 300, second filtering module; 400, first determination module; 500, second determination module. Detailed Implementation

[0064] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0065] Therefore, the following detailed description of the embodiments of the invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely to illustrate selected embodiments of the invention. All other embodiments obtained by those skilled in the art based on the embodiments of the invention without inventive effort are within the scope of protection of the invention.

[0066] It should be noted that similar labels and letters in the following figures indicate similar items. Therefore, once an item is defined in one figure, it does not need to be further defined and explained in subsequent figures.

[0067] Please refer to Figure 1 As shown, an embodiment of the present invention provides a method for identifying mesoscale vortices, comprising steps S10 to S50.

[0068] S10, acquire SLA data, the SLA data including multiple SLA extreme points.

[0069] Mesoscale eddies can cause fluctuations in sea level height. Generally, anti-vortices cause sea level rises, producing positive sea level anomalies; these anomalies appear as raised hills. Conversely, vortices cause sea level falls, producing negative sea level anomalies; these anomalies appear as sunken basins. Therefore, for the SLA data of each satellite altimeter, if the nadir orbit of a satellite altimeter passes through an anti-vortex (bulging sea level anomaly), then the SLA data of that satellite altimeter should show a parabolic curve with a maximum value (opening downwards); if the nadir orbit of a satellite altimeter passes through a cyclonic vortex (sunken sea level anomaly), then the SLA data of that satellite altimeter should show a parabolic curve with a minimum value (opening upwards).

[0070] S20, the SLA data is divided into multiple initial SLA splines based on multiple SLA extreme points, and the multiple initial SLA splines are filtered according to preset filtering conditions to obtain multiple first quasi-vortex SLA splines;

[0071] Each initial SLA spline may include either a central minimum and two endpoint maxima, or a central maximum and two endpoint minimums. Among the selected first quasi-vortex SLA splines, a first quasi-vortex SLA spline with one central minimum and two endpoint maxima may belong to an air vortex and is designated as a quasi-air vortex SLA spline; a first quasi-vortex SLA spline with one central maximum and two endpoint minimums may belong to an anti-air vortex and is designated as a quasi-anti-air vortex SLA spline. Quasi-air vortex splines and quasi-anti-air vortex splines are collectively referred to as quasi-vortex SLA splines, and this characteristic will be used to distinguish between them later.

[0072] S30, select from multiple first quasi-vortex SLA splines that belong to the mesoscale vortex to obtain multiple second quasi-vortex SLA splines.

[0073] When a satellite orbit passes through a strong current, the SLA data from the satellite altimeter can also present as a quasi-vortex SLA spline. However, if multiple satellite orbits pass through an adjacent region within a certain time period, it is possible to determine whether a mesoscale vortex has been passed based on the adjacent first quasi-vortex SLA splines within a certain spatial range. This solves the problem that it is difficult to determine whether a mesoscale vortex belongs to a mesoscale vortex based on only one first quasi-vortex SLA spline.

[0074] S40, determine the vortex boundary and vortex center of the mesoscale vortex based on multiple second quasi-vortex SLA splines.

[0075] S50, determine the basic attribute characteristics based on the vortex boundary and vortex center of the mesoscale vortex.

[0076] In one possible implementation, the preset filtering criteria are:

[0077] The initial SLA spline data points are continuous and without missing points;

[0078] The distance between the center extreme point and the two endpoint extreme points of the initial SLA spline exceeds the first preset threshold.

[0079] The absolute value of the first SLA difference between the central extreme point and the two endpoint extreme points of the initial SLA spline exceeds the second preset threshold.

[0080] Among the data points of adjacent initial SLA splines, the data points of the initial SLA spline whose SLA gradient exceeds a third preset threshold.

[0081] Considering that mesoscale eddies are affected by spatiotemporal scale and the intensity of sea surface height fluctuations they cause, this screening of each initial SLA spline can discard initial SLA splines that are unrelated to mesoscale eddies.

[0082] Please refer to Figure 2 As shown, in one possible implementation, S30 includes S31 to S37.

[0083] S31, using the center point of each of the first quasi-vortex SLA splines as the search center, search for the same type of first quasi-vortex SLA splines within a preset search radius.

[0084] Typically, the spatial scale of mesoscale vortices is between 50 and 200 kilometers. The preset search radius can be between 200 and 400 kilometers, ideally larger than the spatial scale of the mesoscale vortex, to ensure the acquisition of a more closely approximate first quasi-vortex SLA spline. For example, in this invention, the preset search radius is chosen to be 300 kilometers.

[0085] S32, compare the magnitudes of the central extreme points of the first quasi-vortex SLA splines of the same type, and determine the quasi-vortex center SLA spline and the quasi-vortex center point based on the characteristics of the gas vortex and the anti-gas vortex.

[0086] For gaseous vortices, by comparing the magnitudes of the central minima of multiple first quasi-vortex SLA splines, the first quasi-vortex SLA spline with the smallest central minima is selected as the quasi-vortex center SLA spline, and the central minima of the quasi-vortex center SLA spline is taken as the quasi-cyclonic vortex center. For anticyclonic vortices, by comparing the magnitudes of the central maxima of multiple first quasi-vortex SLA splines, the first quasi-vortex SLA spline with the largest central maxima is selected as the quasi-vortex center SLA spline, and the central maxima of the quasi-vortex center SLA spline is taken as the quasi-anticyclonic vortex center (the closer to the mesoscale vortex center, the larger the SLA value. Therefore, the first quasi-vortex SLA spline with the largest central maxima is closest to the mesoscale vortex center).

[0087] S33, calculate the three-dimensional distance between each of the quasi-vortex center points and the spline center points of all the first quasi-vortex SLA splines.

[0088] S34. Compare the distance relationship between the first quasi-vortex SLA spline with a three-dimensional distance less than a preset constant and the quasi-vortex center SLA spline, and determine the first quasi-vortex SLA spline that satisfies the preset relationship condition.

[0089] The preset constant can be between 2 and 6. Considering the existence of a limiting scale for mesoscale vortices, for example, a mesoscale vortex with a radius of 200 km (approximately 2° latitude and longitude) and an vortex amplitude (the SLA difference between the vortex center and the vortex boundary) of 1 m, would have a three-dimensional distance of approximately 2.24 between the quasi-vortex center SLA spline passing through the vortex center and the quasi-vortex center SLA spline at the vortex boundary. Therefore, in this invention, the preset constant is 2.5.

[0090] S35, determine whether the total amount of the first quasi-vortex SLA spline that satisfies the preset relationship condition is greater than the preset total amount threshold.

[0091] For example, the preset total threshold is between 1 and 5. For instance, in the invention, the preset total threshold is 1.

[0092] S36, if the total number of the first quasi-vortex SLA splines that satisfy the preset relationship condition is less than or equal to the preset total threshold, then the first quasi-vortex SLA splines that satisfy the preset relationship condition are removed.

[0093] S37, if the total number of the first quasi-vortex SLA splines that satisfy the preset relationship conditions is greater than the preset total threshold, then the first quasi-vortex SLA splines that satisfy the preset relationship conditions are determined as the second quasi-vortex SLA splines.

[0094] In one possible implementation, the preset relationship condition is:

[0095] The closer the center point of the first quasi-vortex SLA spline is to the center point of the quasi-vortex, the smaller the difference between the SLA value of the center point of the first quasi-vortex SLA spline and the SLA value of the center point of the quasi-vortex.

[0096] The farther the center point of the first quasi-vortex SLA spline is from the center point of the quasi-vortex, the greater the difference between the SLA value of the center point of the first quasi-vortex SLA spline and the SLA value of the center point of the quasi-vortex.

[0097] In one possible implementation, the formula for calculating the three-dimensional distance is:

[0098]

[0099] In the formula: ∆ D It is a three-dimensional distance;

[0100] ( x 0, y 0) represents the latitude and longitude spatial coordinates of the quasi-vortex center point;

[0101] ( x i , y i () represents the latitude and longitude spatial coordinates of the spline center point of the quasi-vortex center spline;

[0102] z 0 represents the SLA value at the quasi-vortex center point;

[0103] z i The SLA value at the center point of the quasi-vortex center spline;

[0104] C x , C y and C z They are respectively x , y and z Characteristic constants of the variable.

[0105] Please refer to Figure 3 As shown, in one possible implementation, S40 includes S41 to S45.

[0106] S41, calculate the maximum gradient value of each of the second quasi-vortex SLA splines, then calculate the average value of all the maximum gradient values, and use the average value as the SLA value of the vortex boundary of the mesoscale vortex.

[0107] S42, calculate the two spline boundary position points and spline normal boundary position points of each of the second quasi-vortex SLA splines with respect to the vortex boundary.

[0108] The spline boundary points are the boundary points of the mesoscale vortex on the second quasi-vortex SLA spline. The SLA values ​​of the boundary are calculated by interpolation based on the spatial coordinates of the data points in the second quasi-vortex SLA spline and the SLA values. Since each second quasi-vortex SLA spline is a parabola formed around its central extremum, there is one spline boundary point on each side of the second quasi-vortex SLA spline.

[0109] The spline normal boundary point is the location of the mesoscale vortex along the normal direction of the second quasi-vortex SLA spline. The SLA value of the vortex boundary is calculated by extrapolating along the spline normal direction based on the SLA value of the vortex center and the SLA value of the spline center point of the second quasi-vortex SLA spline. Each second quasi-vortex SLA spline has only one normal boundary point.

[0110] S43, the Graham scan method is used to calculate the spline boundary position points and spline normal boundary position points of all the second quasi-vortex SLA splines to obtain the convex hull.

[0111] S44, the convex hull is smoothed using the Bézier curve method to obtain the vortex boundary of the mesoscale vortex.

[0112] S45, calculate the centroid based on the vortex boundary, take the centroid as the center of the mesoscale vortex, and take the SLA value of the quasi-vortex center point as the SLA value of the vortex center.

[0113] Please refer to Figure 4 As shown, in one possible implementation, S50 includes S51 to S55.

[0114] S51, determine whether the SLA value of the vortex center is less than the SLA value of the vortex boundary.

[0115] S52, if the SLA value of the vortex center is less than the SLA value of the vortex boundary, it is an air vortex.

[0116] S53, if the SLA value of the vortex center is greater than the SLA value of the vortex boundary, it is an anti-vortex.

[0117] S54, the radius of the circle with the same shape and area is used as the spatial scale of the mesoscale vortex.

[0118] S55, calculate the absolute value of the second SLA difference between the SLA value of the vortex center and the SLA value of the vortex boundary, and use the absolute value of the second SLA difference as the intensity of the mesoscale vortex.

[0119] Please refer to Figure 5 As shown, an embodiment of the present invention provides a device for identifying mesoscale vortices, the device comprising:

[0120] The acquisition module 100 is configured to acquire SLA data, which includes multiple SLA extreme points.

[0121] The first filtering module 200 is configured to divide the SLA data into multiple initial SLA splines based on multiple SLA extreme points, and filter the multiple initial SLA splines according to preset filtering conditions to obtain multiple first quasi-vortex SLA splines.

[0122] The second screening module 300 is configured to screen multiple first quasi-vortex SLA splines that belong to the mesoscale vortex, thereby obtaining multiple second quasi-vortex SLA splines.

[0123] The first determining module 400 is configured to determine the vortex boundary and vortex center of the mesoscale vortex based on multiple second quasi-vortex SLA splines.

[0124] The second determining module 500 is configured to determine basic attribute features based on the vortex boundary and vortex center of the mesoscale vortex.

[0125] This invention provides a computer device, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps of the above-described method.

[0126] This invention provides a computer-readable storage medium storing a computer program, wherein the computer program, when executed by a processor, implements the steps of the above-described method.

[0127] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the methods described above. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, or optical storage, etc. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.

[0128] The above are merely preferred embodiments of the present invention and are not intended to limit the present invention. For those skilled in the art, the present invention can have various modifications and variations. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method of identifying mesoscale vortices, characterized by, The method comprises: acquiring SLA data, the SLA data comprising a plurality of SLA extreme points; dividing the SLA data into a plurality of initial SLA splines according to the plurality of SLA extreme points, screening the plurality of initial SLA splines according to a preset screening condition, and obtaining a plurality of first quasi-vortex SLA splines; screening the first quasi-vortex SLA splines to obtain a plurality of second quasi-vortex SLA splines belonging to the mesoscale vortex, comprising: taking a spline center point of each of the first quasi-vortex SLA splines as a search center to search for first quasi-vortex SLA splines of the same type within a preset search radius; comparing the center extreme points of the first quasi-vortex SLA splines of the same type in size, and determining a quasi-vortex center SLA spline and a quasi-vortex center point according to the characteristics of the air vortex and the characteristics of the anti-air vortex; calculating the three-dimensional distance between each of the quasi-vortex center points and the spline center points of all the first quasi-vortex SLA splines; comparing the distance relationship between the first quasi-vortex SLA splines with the three-dimensional distance less than a preset constant and the quasi-vortex center SLA spline to determine the first quasi-vortex SLA splines satisfying a preset relationship condition; judging whether the total amount of the first quasi-vortex SLA splines satisfying the preset relationship condition is greater than a preset total amount threshold; if the total amount of the first quasi-vortex SLA splines satisfying the preset relationship condition is less than or equal to the preset total amount threshold, eliminating the first quasi-vortex SLA splines satisfying the preset relationship condition; if the total amount of the first quasi-vortex SLA splines satisfying the preset relationship condition is greater than the preset total amount threshold, determining the first quasi-vortex SLA splines satisfying the preset relationship condition as second quasi-vortex SLA splines; determining the vortex boundary and the vortex center of the mesoscale vortex according to the plurality of second quasi-vortex SLA splines; determining the basic attribute characteristics according to the vortex boundary and the vortex center of the mesoscale vortex.

2. The method of identifying mesoscale vortices of claim 1, wherein, The preset screening condition is: the data points of the initial SLA spline are continuous and not missing; the distance between the center extreme point of the initial SLA spline and the two end extreme points exceeds a first preset threshold; the absolute value of the first SLA difference between the center extreme point of the initial SLA spline and the two end extreme points exceeds a second preset threshold; among the data points of the adjacent initial SLA splines, the SLA gradient of the initial SLA spline exceeds a third preset threshold.

3. The method of identifying mesoscale vortices of claim 1, wherein, The calculation formula of the three-dimensional distance is: ； where: Δ D is the three-dimensional distance; ( x 0, y 0) are the longitude and latitude spatial coordinates of the quasi-vortex center point; x i , y i ) longitude and latitude spatial coordinates of the spline center point of the quasi-vortex center spline;​ z 0 is the SLA value of the quasi-vortex center point; z i SLA value for the spline center point of the quasi-vortex center spline; C x , C y and C z are characteristic constants of the variables x , y and z respectively.

4. The method of identifying mesoscale vortices of claim 1, wherein, The determination of the vortex boundary and the vortex center of the mesoscale vortex according to the plurality of second quasi-vortex SLA splines comprises: calculating the gradient maximum value of each of the second quasi-vortex SLA splines, and then calculating the average value of all the gradient maximum values, taking the average value as the SLA value of the vortex boundary of the mesoscale vortex; calculating the two spline boundary position points and the spline normal boundary position points of each of the second quasi-vortex SLA splines with respect to the SLA value of the vortex boundary; calculating the spline boundary position points and the spline normal boundary position points of all the second quasi-vortex SLA splines by using the Graham scan method to obtain a convex hull; The convex hull is smoothed by using a Bezier curve method to obtain a vortex boundary of the mesoscale vortex; A centroid is calculated according to the vortex boundary, the centroid is taken as a vortex center of the mesoscale vortex, and an SLA value of the quasi-vortex center point is taken as an SLA value of the vortex center.

5. The method of identifying mesoscale vortices of claim 1, wherein, The basic attribute features are determined according to the vortex boundary and the vortex center of the mesoscale vortex, including: It is judged whether the SLA value of the vortex center is less than the SLA value of the vortex boundary; If the SLA value of the vortex center is less than the SLA value of the vortex boundary, it is a mesoscale vortex; If the SLA value of the vortex center is greater than the SLA value of the vortex boundary, it is a mesoscale vortex; A radius of a circle with the same area as the vortex boundary is taken as a spatial scale of the mesoscale vortex; An absolute value of a second SLA difference value between the SLA value of the vortex center and the SLA value of the vortex boundary is calculated, and the absolute value of the second SLA difference value is taken as an intensity of the mesoscale vortex.

6. A device for identifying mesoscale vortices, characterized in that The device includes: An acquisition module configured to acquire SLA data, the SLA data including a plurality of SLA extreme points; A first screening module configured to divide the SLA data into a plurality of initial SLA splines according to the plurality of SLA extreme points, and to screen the plurality of initial SLA splines according to a preset screening condition to obtain a plurality of first quasi-vortex SLA splines; A second screening module configured to screen the plurality of first quasi-vortex SLA splines to obtain a plurality of second quasi-vortex SLA splines, including: Taking a spline center point of each of the first quasi-vortex SLA splines as a search center, searching for the same type of first quasi-vortex SLA splines within a preset search radius range; Comparing the center extreme points of the same type of first quasi-vortex SLA splines in size, and determining a quasi-vortex center SLA spline and a quasi-vortex center point according to the mesoscale vortex characteristics and the mesoscale vortex characteristics; Calculating the three-dimensional distance between each of the quasi-vortex center points and the spline center points of all the first quasi-vortex SLA splines; Comparing the distance relationship between the first quasi-vortex SLA splines with the three-dimensional distance less than a preset constant and the quasi-vortex center SLA spline to determine the first quasi-vortex SLA splines satisfying the preset relationship condition; Judging whether the total amount of the first quasi-vortex SLA splines satisfying the preset relationship condition is greater than a preset total amount threshold; If the total amount of the first quasi-vortex SLA splines satisfying the preset relationship condition is less than or equal to the preset total amount threshold, the first quasi-vortex SLA splines satisfying the preset relationship condition are removed; If the total amount of the first quasi-vortex SLA splines satisfying the preset relationship condition is greater than the preset total amount threshold, the first quasi-vortex SLA splines satisfying the preset relationship condition are determined as second quasi-vortex SLA splines; A first determination module configured to determine the vortex boundary and the vortex center of the mesoscale vortex according to the plurality of second quasi-vortex SLA splines; A second determination module configured to determine the basic attribute features according to the vortex center and the vortex boundary of the mesoscale vortex. 7.A computer device, comprising a memory and a processor, wherein the memory stores a computer program, and the computer device is configured to perform the method according to any one of claims 1-6 when the computer program is executed by the processor. The computer program is executed by the processor to implement the steps of the method of any one of claims 1 to 5.

8. A computer readable storage medium, having stored thereon a computer program, characterized in that, The computer program is executed by the processor to implement the steps of the method of any one of claims 1 to 5.

Images