Inversion method of sea surface atmospheric waveguide parameters based on rank correlation method
By using the rank correlation method in the inversion of the atmospheric waveguide parameter, the rank correlation coefficient of radar echo power and ground ray density is simplified, the calculation process is solved, the problem of inversion time in the existing technology is solved, and real-time is realized Or near real-time inversion effect.
Patent Information
- Application Number
- CN202510030818.6
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-01-08
- Publication Date
- 2025-02-14
- Estimated Expiration
- Not applicable · inactive patent
AI Technical Summary
When the prior art inverts the atmospheric waveguide parameters of the sea surface, the complexity of the forward model results in too long operation time, making it difficult to achieve real-time inversion.
The rank correlation method is used to determine the atmospheric waveguide parameters on the sea surface by obtaining the rank correlation coefficient of the radar echo power and the grounding ray density, which simplifies the calculation process.
It effectively shortens the inversion time, improves the real-time nature of inversion, and can directly use the received radar echoes to invert in plains and calm sea surface environments.
Smart Images
Figure CN119439114A_ABST
Abstract
Description
Technical Field
[0001] The present application relates to the field of atmospheric duct parameter inversion, and in particular to the inversion of sea surface atmospheric duct parameters. Background Art
[0002] The inversion of sea surface atmospheric waveguide parameters is extremely important. It can optimize maritime communications, improve communication quality and extend distance. It can enhance radar detection performance, increase target detection distance and positioning accuracy. It can provide accurate meteorological information for navigation and ensure safety. It can also help marine meteorological research, understand boundary layer characteristics, and improve weather forecasts, which is of great significance.
[0003] In the prior art, an improved genetic algorithm (MGA) is used to invert the sea surface atmospheric waveguide parameters. Compared with the classic genetic algorithm (GA) and Bayesian-Monte Carlo method (Bayesian-MCMC), MGA shortens the operation time and improves the inversion accuracy to a certain extent. However, due to the complexity of the forward model (PWE), it is still difficult to apply this method to real-time inversion. The operation time of each iteration of the forward model increases with the increase of the radar transmission frequency, antenna beam width and the number of spatial grid points. Summary of the invention
[0004] In order to overcome the above technical defects, the present application provides an inversion method, storage medium and device for sea surface atmospheric waveguide parameters based on the rank correlation method, which can effectively reduce the inversion time and improve the real-time performance of the inversion.
[0005] To achieve the above objectives, this application is implemented according to the following technical solutions: In a first aspect, the present application provides a method for inverting sea surface atmospheric duct parameters based on a rank correlation method, comprising the following steps: Step S1: Obtain the maximum rank correlation coefficient; Step S2: obtaining an actually observed radar echo power profile, and determining an echo power rank vector corresponding to the echo power vector according to the radar echo power profile; Step S3: acquiring sea surface atmospheric waveguide parameters, and determining the ground ray density rank vector corresponding to the sea surface atmospheric waveguide parameters according to the sea surface atmospheric waveguide parameters; Step S4: determining a rank correlation coefficient based on the echo power rank vector and the ground ray density rank vector; Step S5: Based on the rank correlation coefficient and the maximum rank correlation coefficient, determine the sea surface atmospheric duct parameter corresponding to the rank correlation coefficient.
[0006] Optionally, the sea surface atmosphere duct parameter is provided by n possible sea surface atmosphere duct parameter sets; wherein the n possible sea surface atmosphere duct parameter sets are generated using a Monte Carlo method.
[0007] Optionally, determining the ground ray density rank vector corresponding to the sea surface atmospheric duct parameter according to the sea surface atmospheric duct parameter includes: Substitute the sea surface atmospheric waveguide parameters into the ray tracing model to calculate the ground ray density rank vector corresponding to the sea surface atmospheric waveguide parameters.
[0008] Optionally, determining a rank correlation coefficient based on the echo power rank vector and the ground ray density rank vector comprises: The echo power rank vector and the ground ray density rank vector are calculated using Spearman to determine the rank correlation coefficient.
[0009] Optionally, determining a sea surface atmospheric duct parameter corresponding to the rank correlation coefficient based on the rank correlation coefficient and the maximum rank correlation coefficient includes: Determining whether the rank correlation coefficient is greater than the maximum rank correlation coefficient; If yes, setting the maximum rank correlation coefficient equal to the rank correlation coefficient; Determine whether the number of cycles reaches n times; If yes, then output the sea surface atmospheric duct parameter corresponding to the rank correlation coefficient; In a second aspect, the present invention provides a computer device, comprising a processor and a memory, wherein the memory stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, at least one program, a code set or an instruction set is loaded and executed by the processor to implement the inversion method for sea surface atmospheric duct parameters based on the rank correlation method as described in any one of claims 1 to 5.
[0010] In a third aspect, the present application provides a computer-readable storage medium, wherein the computer-readable storage medium stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, at least one program, a code set or an instruction set is loaded and executed by a processor to implement the inversion method for sea surface atmospheric duct parameters based on the rank correlation method as described in any one of claims 1 to 5.
[0011] This application has the following beneficial effects: This application uses the rank of radar echo power as inversion data. In environments such as plains and calm sea surfaces, the clutter information has little effect on the rank of echo power, and the received radar echo can be directly used for the inversion of surface waveguide parameters. And by calculating the rank correlation coefficient of the echo power rank vector and the ground ray density rank vector, and then determining the sea surface atmospheric waveguide parameters based on the maximum rank correlation coefficient and the rank correlation coefficient, because a linear calculation method is used, it is possible to roughly invert the waveguide layer structure, and the calculation time is greatly shortened, and basically real-time or near real-time inversion can be achieved.
[0012] In addition to the above-described purposes, features and advantages, the present application also has other purposes, features and advantages. The present application will be further described in detail below with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS
[0013] The drawings constituting a part of the present application are used to provide a further understanding of the present application. The illustrative embodiments and descriptions of the present application are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings: Figure 1 It is a schematic diagram of the atmospherically corrected refractive index profiles of the 350 m surface waveguide and the 100 m surface waveguide in this application.
[0014] Figure 2 (a1) is a schematic diagram of the 350m surface waveguide propagation power loss of the present application.
[0015] Figure 2 (a2) is a schematic diagram of the 350m surface waveguide ground propagation loss value of the present application.
[0016] Figure 2 (b1) is a schematic diagram of the 350m surface waveguide ray tracing path diagram of the present application.
[0017] Figure 2 (b2) is a schematic diagram of the 350m surface waveguide ground ray density of the present application.
[0018] Figure 2 (c1) is a schematic diagram of the 100m surface waveguide ray tracing path diagram of the present application.
[0019] Figure 2 (c2) is a schematic diagram of the 100m surface waveguide ground ray density of the present application.
[0020] Figure 3 It is a flow chart of a method for inverting sea surface atmospheric duct parameters based on the rank correlation method of the present application.
[0021] Figure 4It is a schematic diagram comparing the real refractive index profile and the inverted profile at 350 m in this application, wherein the solid line is the given surface waveguide profile and the dotted line is the profile inverted by RTC (the method provided in this application).
[0022] Figure 5 1 is a schematic diagram comparing the inversion results of RTC in the present application with the inversion results using the parabolic wave equation as the forward model combined with the improved genetic algorithm (PWE / MGA), wherein the thin solid line is the profile detected at different distances along the flight direction of the helicopter, and the thick solid line is the profile calculated by its horizontal average value; the dotted line is the profile obtained by RTC inversion, and the dotted line is the profile obtained by PWE / MGA inversion. DETAILED DESCRIPTION
[0023] The embodiments of the present application are described in detail below with reference to the accompanying drawings; however, the present application can be implemented in many different ways as defined and covered by the claims.
[0024] Since the ray tracing / rank correlation method is used in this application to determine the sea surface atmospheric waveguide parameters, the principle is briefly introduced.
[0025] The principle of estimating the surface-based duct (SBD) parameters by ray tracing / rank correlation (RTC) is relatively simple. It is assumed that there is a certain correlation between the radar echo power and the ray density. The ray density corresponding to the point with large echo power is large, and the ray density corresponding to the point with small echo power is small. The above assumption is demonstrated below. Here, the propagation power loss calculated by the parabolic wave equation model is used to simulate the actual radar echo (please refer to the following formula for the relationship between propagation loss and radar echo power).
[0026] (1) in: is the radar echo power; is the propagation loss; For distance, is the normalized scattering coefficient, is a constant determined by the radar system parameters.
[0027] Figure 1The atmospheric corrected refractive index profiles of two different sets of surface waveguides are given. The solid line represents the 350m surface waveguide, whose three parameters Base height, Thickness, and M-deficit are 300m, 50m, and 50M-units respectively; the dashed-dotted line represents the 100m surface waveguide, whose three parameters Base height, Thickness, and M-deficit are 50m, 50m, and 10M-units respectively.
[0028] In order to illustrate the correlation between radar echo power and ray density, combined with Figure 2 The three-profile model given is calculated using the PWE model and the RT model, where the antenna is set at a height of 27m and the radar operating frequency is 5.6GHz. The range of the radar transmission elevation angle is , the spacing is , so that about 500 rays are traced (when the PWE model is used for calculation, the elevation angle is ).
[0029] Figure 2 a1 gives the electromagnetic wave one-way propagation power loss (echo power) calculated for the 350m surface waveguide. The corresponding propagation loss value on the ground is given by Figure 2 a2 is given; Figure 2 b1 gives 500 ray tracing paths for a 350m surface waveguide. Figure 2 b2 gives the grounding ray density (the calculation method of the ray density is as follows: take the horizontal spacing of 1000m, in each horizontal spacing, as long as there is one ray touching the ground, the ray density is recorded as 1, if there are n rays touching the ground, the ray density is recorded as n); Figure 2 c1 and Figure 2 c2 gives the ray tracing path of the 100m surface waveguide and the ray density of the ground, respectively. Figure 2 The three pictures on the left can be seen very intuitively: Figure 2 a1 gives the propagation power loss and Figure 2 The ray tracing paths given by b1 have good similarity; due to the use of different atmospheric correction refractive index profiles for calculation, Figure 2 a1 and Figure 2 The similarity of c1 is poor. The Spearman rank correlation theory is used to calculate the rank correlation coefficient between propagation power loss and ray density. Figure 2 a2 and Figure 2 The rank correlation coefficient between b2 is 0.69, while Figure 2 a2 and Figure 2 The rank correlation coefficient between c2 is only 0.16.
[0030] From the above analysis, it can be seen that the propagation power loss calculated by the same atmospheric corrected refractive index profile has a good correlation with the ray density, while the correlation between different atmospheric corrected refractive index profiles is poor. It can be seen that it is feasible to use RTC to estimate the sea surface atmospheric waveguide parameters.
[0031] In the prior art, an improved genetic algorithm (MGA) is used to invert the sea surface atmospheric waveguide parameters. Compared with the classic genetic algorithm (GA) and Bayesian-Monte Carlo method (Bayesian-MCMC), MGA shortens the operation time and improves the inversion accuracy to a certain extent. However, due to the complexity of the forward model (PWE), it is still difficult to apply this method to real-time inversion. The operation time of each iteration of the forward model increases with the increase of the radar transmission frequency, antenna beam width and the number of spatial grid points.
[0032] In one embodiment, in order to solve the above technical problems, a method for determining sea surface atmospheric duct parameters based on rank correlation method is provided, combining Figure 3 , the method of this application is explained in detail: Step S1: Obtain the maximum rank correlation coefficient; Since the inverted sea surface atmospheric waveguide parameters are determined by determining the maximum rank correlation coefficient in this application, it is necessary to compare the obtained rank correlation coefficients one by one to determine the maximum rank correlation coefficient. However, when the calculation is performed for the first time, the rank correlation coefficient has not been determined. Therefore, in order to facilitate the smooth implementation of the subsequent comparison process, the maximum rank correlation coefficient can be assumed to be 0 during the first calculation, that is, , is the maximum Spearman rank correlation coefficient between the radar echo power and the ground ray density, that is, the rank correlation coefficient is the rank correlation coefficient between the radar echo power and the ground ray density.
[0033] It should be noted that here the maximum rank correlation coefficient is set to zero only in the initial calculation. The maximum rank correlation coefficient will be updated according to each subsequent iterative calculation. The iterative process will be described in detail later and will not be elaborated on here.
[0034] Step S2: obtaining an actually observed radar echo power profile, and determining an echo power rank vector set corresponding to the echo power vector according to the radar echo power profile; Radar echo power profile refers to a curve obtained through radar observation that reflects the change of radar echo power with distance, height or other related spatial coordinates. It depicts the distribution of echo power received by the radar beam at various locations on the propagation path.
[0035] As mentioned above, what is calculated in this article is the rank correlation coefficient between the radar echo power and the ground ray density. The correlation method for calculating the rank correlation coefficient is used to determine the strength of the linear relationship between two or more variables. It is one of the most widely used statistical methods. In traditional statistical methods, the relationship between variables is defined by the correlation coefficient, and the basic idea of the rank method is to use Rank to replace Make statistical inferences (simply put, The rank is In the array The serial number value in, from small to large or from large to small).
[0036] As mentioned above, it can be seen that in this step, after obtaining the actual observed radar echo power profile, it is arranged in order from small to large or from large to small to determine the echo power rank vector set corresponding to the echo power vector .
[0037] It should be noted that in subsequent iterative calculations, the echo power rank vectors are calculated one by one according to the arrangement order of the echo power rank vectors in the echo power rank vector set. , the next time the iterative calculation is performed, the echo power rank vector .
[0038] Step S3: acquiring sea surface atmospheric waveguide parameters, and determining the ground ray density rank vector corresponding to the sea surface atmospheric waveguide parameters according to the sea surface atmospheric waveguide parameters; This application determines the inversion result, i.e., the sea surface atmospheric duct parameter, based on the maximum rank correlation coefficient. Prior to this, it is necessary to use the Monte Carlo method to generate n possible sets of sea surface atmospheric duct parameter values. .
[0039] After generating the set of sea surface atmospheric waveguide parameters, substitute it into the ray tracing model to calculate the ground ray density rank vector set corresponding to the sea surface atmospheric waveguide parameters. .
[0040] It should be noted that in subsequent iterative calculations, the grounding ray density rank vectors are calculated one by one according to their arrangement order in the grounding ray density rank vector set. For example, this calculation uses the grounding ray density rank vector , the next time the iterative calculation is performed, the ground ray density rank vector is used .
[0041] Specifically, in this step, the atmospheric refractive index profile is corrected The calculated ground ray density rank vector , the specific calculation process is as follows: The mathematical expression of the surface waveguide refractive index profile is as follows: (2) in is the mean sea level height; the corresponding atmospheric corrected refractive index is ; is the slope of the refractive index profile of the layer where the base height is located; is the slope of the refractive index profile above the top of the waveguide layer.
[0042] (3) Where x represents the horizontal distance in meters. Indicates the vertical height in m. It represents the local elevation angle of the ray, and its unit is rad, while t is the linear rate of change.
[0043] According to the above formulas (2) and (3), the refractive index profile corrected by the atmosphere can be , calculate the ground ray density rank vector set corresponding to the sea surface atmospheric waveguide parameters .
[0044] Step S4: determining a rank correlation coefficient based on the echo power rank vector and the ground ray density rank vector; After obtaining the echo rank vector and the ground ray density rank vector, Spearman is used to calculate the echo power rank vector and the ground ray density rank vector to determine the rank correlation coefficient. The calculation process of the Spearman test statistic is introduced below, which is suitable for testing the correlation of a series of data pairs.
[0045] The echo power rank vector set Representation vector rank, ground ray density rank vector set Representation vector The rank of . Then the vector and vector The Spearman rank correlation coefficient between It can be expressed by the following formula: (4) It should be noted that in subsequent iterative calculations, each time the rank correlation coefficient is calculated When , an echo power rank vector is selected in the echo power rank vector set in order, and then the ground ray density rank vector corresponding to the order of the echo power rank vector in the echo power rank vector set is selected in the ground ray density rank vector set for calculation. For example, this calculation uses the echo power rank vector and the ground ray density rank vector Get the rank correlation coefficient . Then the next time the iterative calculation is performed, the echo power rank vector is used , and the ground ray density rank vector Get the rank correlation coefficient , and repeat this process until the rank correlation coefficient is calculated. , then the iterative calculation ends.
[0046] Step S5: Based on the rank correlation coefficient and the maximum rank correlation coefficient, determine the sea surface atmospheric duct parameter corresponding to the maximum rank correlation coefficient.
[0047] In calculating the rank correlation coefficient Then, determine the rank correlation coefficient Is it greater than the maximum rank correlation coefficient , if so, it means the rank correlation coefficient calculated at this time is the maximum rank correlation coefficient. At this time, the maximum rank correlation coefficient is set Equal to the rank correlation coefficient Then, the number of cycles is judged, that is, whether the number of cycles reaches n. If so, it means that all iterative calculations have been completed and the comparison of the sizes of all rank correlation coefficients has been completed. The rank correlation coefficient calculated at this time is the maximum rank correlation coefficient required for this round of iterative calculation. , and then output the maximum rank correlation coefficient The corresponding sea surface atmospheric waveguide parameters, if not, then the rank relationship As the maximum rank correlation coefficient for the next round of calculation , repeat the calculation process of steps S1-S5.
[0048] Determine the rank correlation coefficient Is it greater than the maximum rank correlation coefficient , if not, it represents the previously calculated rank correlation coefficient Not the biggest, discard , then make Keep it unchanged, and then determine whether the number of loops reaches n times. If it reaches n times, it means that the current round of iterative calculations is all over, and then output The corresponding sea surface atmospheric waveguide parameters are obtained, and the inversion is terminated. If the number of cycles does not reach n, then Repeat the iterative calculation of steps S1-S5.
[0049] Next, the effectiveness of the inversion method of sea surface atmospheric duct parameters based on the rank correlation method provided in this application is verified: Numerical simulation Without loss of generality, the atmospheric corrected refractive index profile of the 350m surface waveguide given above is simulated as the real profile, and the propagation power loss calculated by the PWE model is simulated as the real radar echo power. In the calculation process, the horizontal distance is taken as 200km, and the interval is 1000m; the vertical height is taken as 1000m, and the interval is 1m. When the RT model is used for calculation, the number of rays touching the ground in each 1000m interval is recorded as the ray density value of the point, so that the horizontal direction can be evenly divided into 200 points.
[0050] Note: The horizontal distance interval and the number of ray density points can be determined according to the specific problem. When using actual radar data inversion, the horizontal distance interval is 600m and the number of ray density points is 83.
[0051] Referring to the boundary conditions of the five-parameter model given by Professor Peter Gerstoft of the U.S. Naval Research Institute, combined with the maximum height that may appear on the actual surface, the boundary values of the three parameters of the surface waveguide and the spacing of each parameter (which can be modified according to specific circumstances) are given as shown in Table 1:
[0052] In this way, about 2800 atmospheric corrected refractive index profiles can be constructed. , calculate each ray and record the vector represents the ground return power vector rank of; vector Represents the touchdown ray density vector The rank of The atmospheric corrected refractive index profile corresponding to (maximum Spearman rank correlation coefficient) is the optimal profile obtained by inversion.
[0053] Figure 4 The real refractive index profile and the inversion profile are given, where the solid line is the given real profile and the dashed line is the profile obtained by RTC inversion. The specific inversion results are: Base height 275m, Thickness 30m, M-deficit 40M-units.
[0054] Inversion results of measured radar echo data In order to verify the application of the RTC method in the inversion of measured radar echo data, a numerical experiment is conducted using the observation data of the waveguide detection test conducted by the U.S. Naval Research Institute on Wallops Island (37.830N, 75.480W) on April 2, 1998. The inversion results of RTC are compared with the inversion results using the parabolic wave equation as the forward model combined with the improved genetic algorithm (PWE / MGA).
[0055] Figure 5 The inversion results of the two methods are given. The thin solid line in the figure represents the atmospheric corrected refractive index profile detected at different distances along the helicopter flight direction, and the thick solid line is the profile calculated from its horizontal average value; the dotted line is the profile obtained by RTC inversion, and its corresponding Spearman rank correlation coefficient is 0.8, and the running time is 737 seconds; the dotted line is the profile obtained by PWE / MGA inversion, and the calculation time is 5430 seconds. Figure 5 It can be seen that: compared with PWE / MGA, the inversion accuracy of RTC is reduced, but it can still roughly invert the waveguide layer structure, and the calculation time is greatly shortened, basically achieving real-time or near real-time inversion. PWE / MGA uses the specific radar echo power value as the inversion data. The actual radar echo power contains a lot of noise information, which requires the received radar echo to be denoised first; while RTC uses the rank of radar echo power as the inversion data. In environments such as plains and calm seas, the clutter information has little effect on the rank of echo power, and the received radar echo can be directly used for the inversion of sea surface atmospheric waveguide parameters.
[0056] This application uses the rank of radar echo power as inversion data. In environments such as plains and calm sea surfaces, the clutter information has little effect on the rank of echo power, and the received radar echo can be directly used for the inversion of surface waveguide parameters. And by calculating the rank correlation coefficient of the echo power rank vector and the ground ray density rank vector, and then determining the sea surface atmospheric waveguide parameters based on the maximum rank correlation coefficient and the rank correlation coefficient, because a linear calculation method is used, it is possible to roughly invert the waveguide layer structure, and the calculation time is greatly shortened, and basically real-time or near real-time inversion can be achieved.
[0057] In some embodiments, a computer-readable storage medium is provided, which stores a computer program. When the computer program is executed by a processor, the processor implements the inversion method for sea surface atmospheric duct parameters based on the rank correlation method provided in the first aspect.
[0058] It will be appreciated by those skilled in the art that all or some of the steps, systems, and functional modules / units in the methods disclosed above may be implemented as software, firmware, hardware, and appropriate combinations thereof. In a hardware implementation, the division between the functional blocks / units mentioned in the above description does not necessarily correspond to the division of physical components; for example, a physical component may have multiple functions, or a function or step may be performed by several physical components in cooperation. Some or all physical components may be implemented as software executed by a processor, such as a central processing unit, a digital signal processor, or a microprocessor, or may be implemented as hardware, or may be implemented as an integrated circuit, such as an application-specific integrated circuit. Such software may be distributed on a computer-readable storage medium, which may include a computer-readable storage medium (or a non-transitory medium) and a communication medium (or a temporary medium).
[0059] As is known to those of ordinary skill in the art, the term computer-readable storage media includes volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storing information (such as computer-readable instructions, data structures, program modules or other data). Computer-readable storage media include, but are not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disk (DVD) or other optical disk storage, magnetic cassettes, magnetic tapes, magnetic disk storage or other magnetic storage devices, or any other medium that can be used to store the desired information and can be accessed by a computer. In addition, it is known to those of ordinary skill in the art that communication media typically contain computer-readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transmission mechanism, and may include any information delivery medium.
[0060] Exemplarily, the computer-readable storage medium may be an internal storage unit of the network management device described in the foregoing embodiment, such as a hard disk or memory of the network management device. The computer-readable storage medium may also be an external storage device of the network management device, such as a plug-in hard disk, a smart memory card (SmartMediaCard, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), etc., equipped on the network management device.
[0061] In some embodiments, a device is provided, comprising a processor and a memory, wherein the memory is used to store a computer program; the processor is used to execute the computer program and implement the inversion method for sea surface atmospheric duct parameters based on the rank correlation method provided in the first aspect of the present application when executing the computer program.
[0062] It should be understood that the processor may be a central processing unit (CPU), and the processor may also be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. Among them, the general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.
[0063] The above description is only the preferred embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims
1. A method for inverting sea surface atmospheric duct parameters based on rank correlation method, characterized in that: The steps include: Step S1: Obtain the maximum rank correlation coefficient; Step S2: obtaining an actually observed radar echo power profile, and determining an echo power rank vector corresponding to the echo power vector according to the radar echo power profile; Step S3: acquiring sea surface atmospheric waveguide parameters, and determining the ground ray density rank vector corresponding to the sea surface atmospheric waveguide parameters according to the sea surface atmospheric waveguide parameters; Step S4: determining a rank correlation coefficient based on the echo power rank vector and the ground ray density rank vector; Step S5: Based on the rank correlation coefficient and the maximum rank correlation coefficient, determine the sea surface atmospheric duct parameter corresponding to the rank correlation coefficient.
2. The method according to claim 1, characterized in that The sea surface atmosphere duct parameter is provided by n possible sea surface atmosphere duct parameter sets; wherein the n possible sea surface atmosphere duct parameter sets are generated using a Monte Carlo method.
3. The method according to claim 1, characterized in that The step of determining the ground ray density rank vector corresponding to the sea surface atmospheric waveguide parameter according to the sea surface atmospheric waveguide parameter comprises: Substitute the sea surface atmospheric waveguide parameters into the ray tracing model to calculate the ground ray density rank vector corresponding to the sea surface atmospheric waveguide parameters.
4. The method according to claim 1, characterized in that: The determining of the rank correlation coefficient based on the echo power rank vector and the ground ray density rank vector comprises: The echo power rank vector and the ground ray density rank vector are calculated using Spearman to determine the rank correlation coefficient.
5. The method according to claim 1, characterized in that The determining, based on the rank correlation coefficient and the maximum rank correlation coefficient, a sea surface atmospheric duct parameter corresponding to the rank correlation coefficient comprises: Determining whether the rank correlation coefficient is greater than the maximum rank correlation coefficient; If yes, setting the maximum rank correlation coefficient equal to the rank correlation coefficient; Determine whether the number of cycles reaches n times; If yes, then output the sea surface atmospheric duct parameter corresponding to the rank correlation coefficient; If not, repeat steps S1-S5.
6. A computer device, characterized in that: The computer device includes a processor and a memory, wherein the memory stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, at least one program, a code set or an instruction set is loaded and executed by the processor to implement the inversion method for sea surface atmospheric duct parameters based on the rank correlation method as described in any one of claims 1 to 5.
7. A computer-readable storage medium, characterized in that: The readable storage medium stores at least one instruction, at least one program, a code set or an instruction set, and the at least one instruction, at least one program, a code set or an instruction set is loaded and executed by a processor to implement the inversion method of sea surface atmospheric duct parameters based on the rank correlation method as described in any one of claims 1 to 5.