A Radar Wave Parameter Inversion Method and Computer-Readable Medium

Abstract

This invention proposes a radar wave parameter inversion method and a computer-readable medium. The invention calculates the Doppler frequency shift sequence using an energy equalization moment estimation method, and then calculates the space-time radial velocity sequence using the Doppler frequency shift sequence. A quasi-binary variational mode decomposition method is used to decompose the sequence, obtaining multiple two-dimensional intrinsic mode velocity sequences. The space-time radial velocity sequence is reconstructed using each two-dimensional intrinsic mode velocity sequence, resulting in a reconstructed space-time radial velocity sequence. The one-dimensional average velocity spectrum of each radar antenna is further calculated using the reconstructed space-time radial velocity sequence. The radial wave height spectrum of each radar antenna is calculated by combining the one-dimensional average velocity spectrum with the modulation transfer function. The radial wave height spectra of multiple radar antennas are accumulated and synthesized to obtain an undirected wave height spectrum. The significant wave height and average wave period are calculated using the zero-order and first-order moments of the undirected wave height spectrum. This invention improves the accuracy of wave parameter inversion.

Description

Technical Field

[0001] This invention belongs to the field of wave parameter inversion for wave measuring radar, and particularly relates to a radar wave parameter inversion method and a computer-readable medium. Background Technology

[0002] Ocean wave observation has a crucial impact on marine engineering, offshore operations, and marine environmental monitoring. Traditional wave observation equipment, such as wave-measuring buoys, uses fixed-point measurement methods, providing only locally effective wave observation data, which is insufficient to meet practical needs. With the development of remote sensing technology, radar and other remote sensing equipment are widely used in marine monitoring due to their advantages such as all-weather operation, good real-time performance, and relatively low cost. Coherent S-band radar, as a remote sensing device with high-altitude temporal resolution, retrieves gravity wave information based on Bragg scattering of electromagnetic waves with capillary waves on the sea surface. Because the signals emitted by coherent microwave radar have all definite phase angles, this significant characteristic provides the ability to acquire radial Doppler velocities. The direct relationship between radial velocity and wave height spectrum can be used to derive wave parameters such as effective wave height and average wave period.

[0003] The high-altitude temporal resolution of coherent S-band radar allows it to detect nonlinear details on the sea surface, such as breaking waves. Breaking waves introduce nonlinear features like group lines into the wave number frequency spectrum. These nonlinear features introduce additional components in the low-frequency range of the estimated wave spectrum and reduce the main wave energy, leading to an overestimation of wave parameters, especially the mean wave period. Therefore, a series of complex processing steps are required on the radar echo signal to improve the accuracy of wave parameter inversion. Summary of the Invention

[0004] To address the above problems, this invention proposes a radar wave parameter inversion method and a computer-readable medium.

[0005] The technical solution of this invention is a radar wave parameter inversion method, comprising the following steps:

[0006] Step 1: Calculate the Doppler frequency shift sequence using the energy equalization moment estimation method based on the radar echo Doppler spectrum, and calculate the space-time radial velocity sequence in combination with the Doppler frequency shift sequence;

[0007] Step 2: Decompose the space-time radial velocity sequence using the quasi-binary variational mode decomposition method to obtain multiple two-dimensional intrinsic mode velocity sequences;

[0008] Step 3: Reconstruct the space-time radial velocity sequence by combining the velocity sequences of each two-dimensional intrinsic mode to obtain the reconstructed space-time radial velocity sequence;

[0009] Step 4: Perform a two-dimensional fast Fourier transform on the reconstructed space-time radial velocity sequence to obtain the wavenumber frequency spectrum of the reconstructed space-time radial velocity sequence; obtain the velocity power spectrum of the reconstructed space-time radial velocity sequence based on the square of the modulus of the wavenumber frequency spectrum of the reconstructed space-time radial velocity sequence; integrate the velocity power spectrum of the reconstructed space-time radial velocity sequence in the wavenumber direction, and further divide it by the integral width in the wavenumber direction to obtain the one-dimensional average velocity spectrum of each radar antenna; combine the one-dimensional average velocity spectrum of each radar antenna with the modulation transfer function to calculate the radial wave height spectrum of each radar antenna.

[0010] Step 5: Accumulate the radial wave height spectra of multiple radar antennas to obtain the synthesized radial wave height spectrum, and further multiply the synthesized radial wave height spectrum by the angular resolution to obtain the omnidirectional wave height spectrum.

[0011] Step 6: Calculate the effective wave height by combining the zeroth moment of the undirected wave height spectrum, and calculate the average wave period by combining the zeroth moment and the first moment of the undirected wave height spectrum.

[0012] As a preferred embodiment, step 1 further calculates the space-time radial velocity sequence as follows:

[0013] Multiplying the Doppler frequency shift sequence by half the wavelength of the radar electromagnetic wave yields the space-time radial velocity sequence.

[0014] Preferably, the frequencies corresponding to the multiple two-dimensional intrinsic mode velocity sequences described in step 2 are arranged sequentially from high frequency to low frequency, i.e., u i (x,t;n), i∈[1,K], with the corresponding frequencies arranged from high frequency to low frequency;

[0015]

[0016] x∈[1,X], t∈[1,T], n∈[1,N]

[0017] Where X represents the number of range elements, T represents the number of acquisition times, N represents the number of radar antennas, K represents the number of two-dimensional intrinsic mode velocity sequences, and u i (x,t;n) represents the modal velocity corresponding to the nth radar antenna at the tth acquisition time of the xth range element in the ith two-dimensional intrinsic modal velocity sequence. r (x,t;n) represents the space-time radial velocity of the radar antenna at the x-th range element and the t-th acquisition time in the space-time radial velocity sequence;

[0018] As a preferred embodiment, step 3 involves reconstructing the spacetime radial velocity sequence, using the following specific formula:

[0019]

[0020] x∈[1,X], t∈[1,T], n∈[1,N]

[0021] β=α·E K / E K-1

[0022] in, Let u be the spatiotemporal radial velocity sequence reconstructed by the radar antenna at the t-th acquisition time of the x-th range element. i (x,t;n) represents the modal velocity corresponding to the x-th range element and the n-th radar antenna at the t-th acquisition time in the i-th two-dimensional intrinsic modal velocity sequence, u K-1 (x,t;n) represents the modal velocity corresponding to the x-th range element and the n-th radar antenna at the t-th acquisition time in the (K-1)-th two-dimensional intrinsic modal velocity sequence, where β is the compensation factor, α is the compensation coefficient, and E K E represents the energy of the Kth two-dimensional intrinsic mode velocity sequence. K-1 X represents the energy of the (K-1)th two-dimensional intrinsic mode velocity sequence, X represents the number of range elements, T represents the number of acquisition times, N represents the number of radar antennas, and K represents the number of two-dimensional intrinsic mode velocity sequences.

[0023] E i Defined as the energy of the i-th two-dimensional intrinsic mode velocity sequence, i∈[1,K];

[0024] The energy calculation process for the i-th two-dimensional intrinsic mode velocity sequence is as follows:

[0025] The wavenumber frequency spectrum of the i-th two-dimensional intrinsic mode velocity sequence is obtained by performing a two-dimensional fast Fourier transform on the i-th two-dimensional intrinsic mode velocity sequence.

[0026] The energy of the i-th two-dimensional intrinsic mode velocity sequence is obtained by integrating the modulus of the wavenumber frequency spectrum of the i-th two-dimensional intrinsic mode velocity sequence in two dimensions over wavenumber and frequency.

[0027] Preferably, step 4 involves calculating the radial wave height spectrum of each radar antenna by combining the one-dimensional average velocity spectrum of each antenna with the modulation transfer function, as detailed below:

[0028]

[0029] f o ∈[1,F], n∈[1,N]

[0030] Where F represents the number of frequency points, and N represents the number of radar antennas. For the radial wave height spectrum of the nth radar antenna, the f-th wave height is... o The amplitude of the wave spectrum power density at each frequency point For the one-dimensional average velocity spectrum of the nth radar antenna, the f-th... o Velocity spectral amplitude corresponding to each frequency point

[0031] TF is the modulation transfer function, specifically defined as follows:

[0032]

[0033] Where ω is the angular frequency corresponding to the capillary wave of the water particle, k is the wave number corresponding to the capillary wave of the water particle, tanh represents the hyperbolic tangent function, θ is the incident angle, and Δf o For frequency resolution, M represents the angular resolution, M represents the number of time sampling points, and d represents the water depth.

[0034] As a preferred embodiment, the technique described in step 6 for calculating the effective wave height is as follows:

[0035]

[0036] Where m0 represents the zeroth moment of the undirected wave height spectrum, H s Indicates the significant wave height;

[0037] Step 6, calculating the average wave period, is as follows:

[0038] T av =m0 / m1

[0039] Where m0 represents the zeroth moment of the undirected wave height spectrum, m1 represents the first moment of the undirected wave height spectrum, and T av This indicates the average wave period.

[0040] The present invention also provides a computer-readable medium storing a computer program executed by an electronic device, which, when run on the electronic device, performs the steps of the radar wave parameter inversion method.

[0041] This invention improves the accuracy of ocean wave parameter inversion. Attached Figure Description

[0042] Figure 1 : Flowchart of the method according to an embodiment of the present invention;

[0043] Figure 2 : The original space-time radial velocity sequence diagram of an embodiment of the present invention;

[0044] Figure 3 The wavenumber frequency spectrum of the original space-time radial velocity sequence in this embodiment of the invention;

[0045] Figure 4 : The dominant wave mode diagram after QB-VMD decomposition in an embodiment of the present invention;

[0046] Figure 5 : Group line mode diagram after QB-VMD decomposition in an embodiment of the present invention;

[0047] Figure 6 The wavenumber frequency spectrum of the main wave mode diagram after QB-VMD decomposition in this embodiment of the invention;

[0048] Figure 7 Wavenumber frequency spectrum of the group line mode after QB-VMD decomposition according to an embodiment of the present invention;

[0049] Figure 8 : Reconstructed spacetime radial velocity sequence diagram according to an embodiment of the present invention;

[0050] Figure 9 The wavenumber frequency spectrum of the reconstructed spacetime radial velocity sequence according to an embodiment of the present invention;

[0051] Figure 10 Comparison of the spectral density of the non-directional waveguide in an embodiment of the present invention;

[0052] Figure 11 Comparison chart of inverted wave parameters and buoy results in an embodiment of the present invention. Detailed Implementation

[0053] 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, and 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.

[0054] In specific implementation, the method proposed in the technical solution of this invention can be automatically executed by those skilled in the art using computer software technology. System devices for implementing the method, such as computer-readable storage media storing the corresponding computer program of the technical solution of this invention and computer equipment including the computer program running the corresponding computer program, should also be within the protection scope of this invention.

[0055] The following is combined Figures 1 to 11 The specific implementation of the present invention is a radar wave parameter inversion method.

[0056] like Figure 1 The diagram shown is a flowchart of a method according to an embodiment of the present invention, which specifically includes the following steps:

[0057] Step 1: Calculate the Doppler frequency shift sequence using the energy equalization moment estimation method based on the radar echo Doppler spectrum, and calculate the space-time radial velocity sequence in combination with the Doppler frequency shift sequence;

[0058] Step 1 involves further calculating the space-time radial velocity sequence, as detailed below:

[0059] Multiplying the Doppler frequency shift sequence by half the wavelength of the radar electromagnetic wave yields the space-time radial velocity sequence.

[0060] Step 2: Decompose the space-time radial velocity sequence using the quasi-binary variational mode decomposition method to obtain multiple two-dimensional intrinsic mode velocity sequences;

[0061] The frequencies corresponding to the multiple two-dimensional intrinsic mode velocity sequences are arranged sequentially from high frequency to low frequency, i.e., u i (x,t;n), i∈[1,K], with the corresponding frequencies arranged from high frequency to low frequency;

[0062]

[0063] x∈[1,X], t∈[1,T], n∈[1,N]

[0064] Where N=6 represents the number of radar antennas, X=80 represents the number of range cells, T=256 represents the number of acquisition times, K=5 represents the number of two-dimensional intrinsic mode velocity sequences, and u i (x,t;n) represents the modal velocity corresponding to the nth radar antenna at the tth acquisition time of the xth range element in the ith two-dimensional intrinsic modal velocity sequence. r (x,t;n) represents the space-time radial velocity of the radar antenna at the x-th range element and the t-th acquisition time in the space-time radial velocity sequence;

[0065] Step 3: Reconstruct the space-time radial velocity sequence by combining the velocity sequences of each two-dimensional intrinsic mode to obtain the reconstructed space-time radial velocity sequence;

[0066] Step 3 involves reconstructing the spacetime radial velocity sequence, using the following formula:

[0067]

[0068] x∈[1,X], t∈[1,T], n∈[1,N]

[0069] β=α·E K / E K-1

[0070] in, Let u be the spatiotemporal radial velocity sequence reconstructed by the radar antenna at the t-th acquisition time of the x-th range element. i (x,t;n) represents the modal velocity corresponding to the x-th range element and the n-th radar antenna at the t-th acquisition time in the i-th two-dimensional intrinsic modal velocity sequence, u K-1(x,t;n) represents the modal velocity corresponding to the nth radar antenna at the t-th acquisition time of the x-th range element in the (K-1)-th two-dimensional intrinsic modal velocity sequence, where β is the compensation factor, α=0.8 is the compensation coefficient, and E K E represents the energy of the Kth two-dimensional intrinsic mode velocity sequence. K-1 X represents the energy of the (K-1)th two-dimensional intrinsic mode velocity sequence, X = 80 represents the number of range elements, T = 256 represents the number of acquisition times, N = 6 represents the number of radar antennas, and K = 5 represents the number of two-dimensional intrinsic mode velocity sequences.

[0071] E i Defined as the energy of the i-th two-dimensional intrinsic mode velocity sequence, i∈[1,K];

[0072] The energy calculation process for the i-th two-dimensional intrinsic mode velocity sequence is as follows:

[0073] The wavenumber frequency spectrum of the i-th two-dimensional intrinsic mode velocity sequence is obtained by performing a two-dimensional fast Fourier transform on the i-th two-dimensional intrinsic mode velocity sequence.

[0074] The energy of the i-th two-dimensional intrinsic mode velocity sequence is obtained by integrating the modulus of the wavenumber frequency spectrum of the i-th two-dimensional intrinsic mode velocity sequence in two dimensions over wavenumber and frequency.

[0075] Step 4: Perform a two-dimensional fast Fourier transform on the reconstructed space-time radial velocity sequence to obtain the wavenumber frequency spectrum of the reconstructed space-time radial velocity sequence; obtain the velocity power spectrum of the reconstructed space-time radial velocity sequence based on the square of the modulus of the wavenumber frequency spectrum of the reconstructed space-time radial velocity sequence; integrate the velocity power spectrum of the reconstructed space-time radial velocity sequence in the wavenumber direction, and further divide it by the integral width in the wavenumber direction to obtain the one-dimensional average velocity spectrum of each radar antenna; combine the one-dimensional average velocity spectrum of each radar antenna with the modulation transfer function to calculate the radial wave height spectrum of each radar antenna.

[0076] Step 4 involves calculating the radial wave height spectrum of each radar antenna by combining the one-dimensional average velocity spectrum with the modulation transfer function, as detailed below:

[0077]

[0078] f o ∈[1,F], n∈[1,N]

[0079] Where F = 128 represents the number of frequency points, and N = 6 represents the number of radar antennas. For the radial wave height spectrum of the nth radar antenna, the f-th wave height is... o The amplitude of the wave spectrum power density at each frequency point For the one-dimensional average velocity spectrum of the nth radar antenna, the f-th... o Velocity spectral amplitude corresponding to each frequency point

[0080] TF is the modulation transfer function, specifically defined as follows:

[0081]

[0082] Where ω is the angular frequency corresponding to the capillary wave of the water particle, k is the wavenumber corresponding to the capillary wave of the water particle, tanh represents the hyperbolic tangent function, θ is the incident angle, and Δf o =0.0078Hz is the frequency resolution. For angular resolution, M = 256 is the number of time sampling points, and d = 50m is the water depth;

[0083] Step 5: Accumulate the radial wave height spectra of multiple radar antennas to obtain the synthesized radial wave height spectrum, and further multiply the synthesized radial wave height spectrum by the angular resolution to obtain the omnidirectional wave height spectrum.

[0084] Step 6: Calculate the effective wave height by combining the zeroth moment of the undirected wave height spectrum, and calculate the average wave period by combining the zeroth moment and the first moment of the undirected wave height spectrum.

[0085] Step 6 describes the calculation of the effective wave height, as follows:

[0086]

[0087] Where m0 represents the zeroth moment of the undirected wave height spectrum, H s Indicates the significant wave height;

[0088] Step 6, calculating the average wave period, is as follows:

[0089] T av =m0 / m1

[0090] Where m0 represents the zeroth moment of the undirected wave height spectrum, m1 represents the first moment of the undirected wave height spectrum, and T av This indicates the average wave period.

[0091] Figure 2 It is the original radial velocity sequence. Figure 3 It is the wavenumber frequency spectrum corresponding to the original radial velocity sequence.

[0092] like Figure 2 As shown, strong interference fringes and yellow spots appear in the radial velocity sequence; as Figure 3 As shown, the wavenumber frequency spectrum includes group line characteristics in addition to the first-order dispersion relation.

[0093] Figure 4 This is the dominant mode diagram after QB-VMD decomposition. Figure 4 Interference fringes with a period of approximately 5 seconds are clearly visible, while Figure 5 This is a group line mode diagram after QB-VMD decomposition, and a large number of yellow spots can be seen in the diagram. Figure 6 and Figure 7 The energies in the wavenumber frequency spectrum shown exhibit first-order dispersion relation and group line distribution, respectively.

[0094] like Figure 8 As shown, the yellow spots in the figure are suppressed, and the interference fringes are enhanced.

[0095] like Figure 9 As shown in the figure, the group line energy in the wavenumber frequency spectrum is filtered out, while the dispersion relation energy is retained and enhanced.

[0096] Figure 10 This is the undirected wave height spectrum synthesized from six antennas. Compared to the unprocessed wave height spectrum, the low-frequency components of the wave height spectrum retrieved by the proposed algorithm are suppressed, while the main wave component is enhanced. Compared with the buoy observation results, the effective wave height and average wave period measured by the radar under unprocessed conditions are 2.09 m and 6.55 s, respectively. The results obtained by the proposed algorithm are 1.92 m and 4.86 s, respectively. The proposed algorithm's results are closer to the buoy results of 1.90 m and 5.00 s.

[0097] Figure 11 The results show a comparison between the inverted wave parameters and the buoy results. The correlation coefficients between the significant wave height and the mean wave period and the buoy are 0.97 and 0.80, respectively, and the root mean square errors are 0.12 m and 0.49 s, respectively, indicating that the wave parameters obtained by this invention have high accuracy.

[0098] A specific embodiment of the present invention also provides a computer-readable medium.

[0099] The computer-readable medium is a server workstation;

[0100] The server workstation stores the computer program executed by the electronic device. When the computer program runs on the electronic device, it causes the electronic device to execute the steps of the radar wave parameter inversion method of the present invention.

[0101] It should be understood that any parts not described in detail in this specification belong to the prior art.

[0102] It should be understood that the above description of the preferred embodiments is quite detailed, but it should not be considered as a limitation on the scope of protection of this invention. Those skilled in the art, under the guidance of this invention, can make substitutions or modifications without departing from the scope of protection of the claims of this invention, and all such substitutions or modifications fall within the scope of protection of this invention. The scope of protection of this invention should be determined by the appended claims.

Claims

1. A radar wave parameter inversion method, characterized in that: The reconstructed space-time radial velocity sequence is obtained by sequentially applying the energy equilibrium moment estimation method, the quasi-binary variational mode decomposition method, and sequence reconstruction using the radar echo Doppler spectrum. The reconstructed space-time radial velocity sequence is used to calculate the one-dimensional average velocity spectrum of each radar antenna, and the radial wave height spectrum of each radar antenna is calculated by combining the modulation transfer function. The radial wave height spectra of multiple radar antennas are synthesized and accumulated to obtain the directional wave height spectrum; The effective wave height and mean wave period are calculated by combining the zero-order moment and first-order moment of the undirected wave height spectrum.

2. The radar wave parameter inversion method according to claim 1, characterized in that, Includes the following steps: Step 1: Calculate the Doppler frequency shift sequence using the energy equalization moment estimation method based on the radar echo Doppler spectrum, and calculate the space-time radial velocity sequence in combination with the Doppler frequency shift sequence; Step 2: Decompose the space-time radial velocity sequence using the quasi-binary variational mode decomposition method to obtain multiple two-dimensional intrinsic mode velocity sequences; Step 3: Reconstruct the space-time radial velocity sequence by combining the velocity sequences of each two-dimensional intrinsic mode to obtain the reconstructed space-time radial velocity sequence; Step 4: Perform a two-dimensional fast Fourier transform on the reconstructed space-time radial velocity sequence to obtain the wavenumber frequency spectrum of the reconstructed space-time radial velocity sequence; obtain the velocity power spectrum of the reconstructed space-time radial velocity sequence based on the square of the modulus of the wavenumber frequency spectrum of the reconstructed space-time radial velocity sequence; integrate the velocity power spectrum of the reconstructed space-time radial velocity sequence in the wavenumber direction, and further divide it by the integral width in the wavenumber direction to obtain the one-dimensional average velocity spectrum of each radar antenna; combine the one-dimensional average velocity spectrum of each radar antenna with the modulation transfer function to calculate the radial wave height spectrum of each radar antenna. Step 5: Accumulate the radial wave height spectra of multiple radar antennas to obtain the synthesized radial wave height spectrum, and further multiply the synthesized radial wave height spectrum by the angular resolution to obtain the omnidirectional wave height spectrum. Step 6: Calculate the effective wave height by combining the zeroth moment of the undirected wave height spectrum, and calculate the average wave period by combining the zeroth moment and the first moment of the undirected wave height spectrum.

3. The radar wave parameter inversion method according to claim 2, characterized in that: Step 1 involves further calculating the space-time radial velocity sequence, as detailed below: Multiplying the Doppler frequency shift sequence by half the wavelength of the radar electromagnetic wave yields the space-time radial velocity sequence.

4. The radar wave parameter inversion method according to claim 3, characterized in that: The frequencies corresponding to the multiple two-dimensional intrinsic mode velocity sequences mentioned in step 2 are arranged sequentially from high frequency to low frequency, i.e., u i (x,t;n), i∈[1,K], with the corresponding frequencies arranged from high frequency to low frequency; x∈[1,X], t∈[1,T], n∈[1,N] Where X represents the number of range elements, T represents the number of acquisition times, N represents the number of radar antennas, K represents the number of two-dimensional intrinsic mode velocity sequences, and u i (x,t;n) represents the modal velocity corresponding to the nth radar antenna at the tth acquisition time of the xth range element in the ith two-dimensional intrinsic modal velocity sequence. r (x,t;n) represents the space-time radial velocity of the radar antenna at the x-th range element and the t-th acquisition time in the space-time radial velocity sequence.

5. The radar wave parameter inversion method according to claim 4, characterized in that: Step 3 involves reconstructing the spacetime radial velocity sequence, using the following formula: x∈[1,X], t∈[1,T], n∈[1,N] β=α·E K / E K-1 in, Let u be the spatiotemporal radial velocity sequence reconstructed by the radar antenna at the t-th acquisition time of the x-th range element. i (x,t;n) represents the modal velocity corresponding to the x-th range element and the n-th radar antenna at the t-th acquisition time in the i-th two-dimensional intrinsic modal velocity sequence, u K-1 (x,t;n) represents the modal velocity corresponding to the x-th range element and the n-th radar antenna at the t-th acquisition time in the (K-1)-th two-dimensional intrinsic modal velocity sequence, where β is the compensation factor, α is the compensation coefficient, and E K E represents the energy of the Kth two-dimensional intrinsic mode velocity sequence. K-1 X represents the energy of the (K-1)th two-dimensional intrinsic mode velocity sequence, X represents the number of range elements, T represents the number of acquisition times, N represents the number of radar antennas, and K represents the number of two-dimensional intrinsic mode velocity sequences. E i Defined as the energy of the i-th two-dimensional intrinsic mode velocity sequence, i∈[1,K]; The energy calculation process for the i-th two-dimensional intrinsic mode velocity sequence is as follows: The wavenumber frequency spectrum of the i-th two-dimensional intrinsic mode velocity sequence is obtained by performing a two-dimensional fast Fourier transform on the i-th two-dimensional intrinsic mode velocity sequence. The energy of the i-th two-dimensional intrinsic mode velocity sequence is obtained by performing a two-dimensional integral of the wavenumber frequency spectrum over both wavenumber and frequency.

6. The radar wave parameter inversion method according to claim 5, characterized in that: Step 4 involves calculating the radial wave height spectrum of each radar antenna by combining the one-dimensional average velocity spectrum with the modulation transfer function, as detailed below: f o ∈[1,F],n∈[1,N] Where F represents the number of frequency points, and N represents the number of radar antennas. For the radial wave height spectrum of the nth radar antenna, the f-th wave height is... o The amplitude of the wave spectrum power density at each frequency point For the one-dimensional average velocity spectrum of the nth radar antenna, the f-th... o The velocity spectrum amplitude corresponding to each frequency point; TF stands for modulation transfer function.

7. The radar wave parameter inversion method according to claim 6, characterized in that: Step 6 describes the calculation of the effective wave height, as follows: Where m0 represents the zeroth moment of the undirected wave height spectrum, H s Indicates the significant wave height; Step 6, calculating the average wave period, is as follows: T av =m0 / m1 Where m0 represents the zeroth moment of the undirected wave height spectrum, m1 represents the first moment of the undirected wave height spectrum, and T av This indicates the average wave period.

8. A computer-readable medium, characterized in that, It stores a computer program executed by an electronic device, which, when run on the electronic device, causes the electronic device to perform the steps of the method as described in any one of claims 1-7.

Images