CN116430337A - Sea wave parameter space-time inversion method using shipborne coherent S-band radar - Google Patents
Sea wave parameter space-time inversion method using shipborne coherent S-band radar Download PDFInfo
- Publication number
- CN116430337A CN116430337A CN202310265278.0A CN202310265278A CN116430337A CN 116430337 A CN116430337 A CN 116430337A CN 202310265278 A CN202310265278 A CN 202310265278A CN 116430337 A CN116430337 A CN 116430337A
- Authority
- CN
- China
- Prior art keywords
- wave
- spectrum
- space
- wave number
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000001427 coherent effect Effects 0.000 title claims abstract description 39
- 238000001228 spectrum Methods 0.000 claims abstract description 126
- 238000001914 filtration Methods 0.000 claims abstract description 20
- 230000033001 locomotion Effects 0.000 claims abstract description 9
- 230000003044 adaptive effect Effects 0.000 claims abstract description 6
- 238000006243 chemical reaction Methods 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000004422 calculation algorithm Methods 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 10
- 230000001133 acceleration Effects 0.000 claims description 9
- 238000002592 echocardiography Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000007792 addition Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a space-time inversion method for sea wave parameters by utilizing a shipborne coherent S-band radar, which comprises the following steps: acquiring a space-time Doppler spectrum and extracting average Doppler frequency offset; estimating a space-time dimensional radial velocity sequence and converting the space-time dimensional radial velocity sequence into a wave number frequency spectrum; estimating wave field energy and group line energy distribution and designing an adaptive filter; filtering the wave number frequency spectrum and reconstructing a space-time dimension radial Doppler velocity sequence; estimating a wave number direction spectrum; repeating the above process to obtain six wave number frequency spectrums and inverting the undirected wave spectrum; and (5) carrying out sea wave parameter inversion. The method solves the problem of influence of broken waves and ship motion on the sea wave parameter inversion precision, not only maintains the advantages of the shore-based coherent microwave radar, but also makes up the defects of the shore-based coherent microwave radar in spatial ductility and time continuity, and greatly improves the robustness and precision of sea wave parameter inversion of the ship-borne coherent S-band radar.
Description
Technical Field
The invention belongs to the technical field of ocean remote sensing of microwave radars, and particularly relates to a space-time inversion method for sea wave parameters by utilizing a shipborne coherent S-band radar.
Background
In recent years, the research of applying the shipborne coherent microwave radar to marine environment monitoring starts, not only the advantages of the shore-based coherent microwave radar are reserved, but also the defects of the shore-based coherent microwave radar in terms of spatial ductility and time continuity are overcome. The shipborne coherent microwave radar generates Bragg scattering through interaction of vertically polarized electromagnetic waves which are diffracted and propagated along the ocean surface and centimeter waves on the ocean surface, and obtains echo Doppler spectrum inverting the ocean surface state, so that rich wave parameters including wave height, wave period, wave direction and the like can be extracted. At present, the sea wave parameters are extracted mainly by utilizing a space-dimensional shipborne radar sea echo velocity sequence to invert the wave height spectrum, and then the wave height spectrum is integrated. However, due to the fact that the modulation effect of the broken waves can cause the increase of low-frequency components in the wave high spectrum and the shift of wave high-spectrum energy caused by the motion of the ship, the estimated ocean parameters have large estimation errors, the modulation of the broken waves cannot be well removed by using the space dimension speed sequence, the influence of the motion of the ship is eliminated, and therefore the fact that the sea wave parameters are extracted by using the space dimension speed sequence has a certain limitation and the sea wave parameters cannot be accurately extracted. In order to widen the detection range of the sea wave and improve the accuracy of sea wave parameter extraction, a space-time dimension ship-borne radar sea echo speed sequence is adopted to carry out sea wave parameter inversion.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a space-time inversion method for wave parameters by using a shipborne coherent S-band radar, which is used for solving the accuracy problem of the inversion of the wave parameters of the shipborne coherent S-band radar and the problem that the performance is influenced by broken waves and ship motions, and specifically comprises the following steps:
Step 3, estimating the position parameters of wave field energy distribution in the wave number frequency spectrum by adopting a least square fitting algorithmAnd slope parameter s of group line energy distribution g And combining the two parameters to design a two-dimensional binary adaptive filter h a (k,f);
and 7, calculating the effective wave height Hs and the average wave period Tav through the first-order moment and the second-order moment.
In addition, in the step 1, the average Doppler frequency offset of the Doppler spectrum of the antenna echo in each direction is extracted by utilizing a moment estimation method, firstly, the Doppler spectrum peak of the ocean echo space-time dimension is searched through the spectrum peak, and then, the left and right frequency offset f at 18dB of the spectrum peak is searched l And f r Finally, extracting average Doppler frequency offset f by a moment estimation method d The method comprises the following steps:
where f is the Doppler frequency and σ is the Doppler spectrum.
Furthermore, the step 2 is a hollow time-dimensional radial Doppler velocity sequence v (r, t, beta) n ) From the average Doppler frequency offset f d The specific calculation mode is as follows:
wherein lambda is the electromagnetic wavelength, beta n Represents the nth th The pointing direction of the root antenna, r is the radial distance, and t is the sampling time.
2-dimensional velocity sequence v (r, t, beta) using 2-dimensional fast fourier transform n ) Conversion to wave number frequency spectrum V r (k,f,β n ) The method comprises the following steps:
V r (k,f,β n )=2DFFT[v(r,t,β n )] (3)
where k is the wave number, f is the Doppler frequency, beta n Represents the nth th The pointing direction of the root antenna, r is the radial distance, t is the sampling time, and 2DFFT represents the two-dimensional fast fourier transform.
Moreover, the least square fitting algorithm in the step 3 finds the best function matching of the data by minimizing the sum of squares of errors, and the wave field energy distribution is no longer along the first-order wave dispersion relation omega due to the influence of the forward speed of the ship 2 =gk, estimation of the position parameters of the wave field energy distribution in the wave number frequency spectrum using least squares fitting algorithmIt is estimated along the following formula:
wherein ω is the angular frequency of the ocean wave, k is the wave number of the ocean wave,g is the gravity acceleration, which is the included angle between the pointing direction of the antenna and the forward motion direction of the ship.
Estimating slope parameter s of group line energy distribution by least square fitting algorithm g Fitting the distribution of group line energy, namely:
ω=s g k (5)
wherein ω is the angular frequency of the ocean wave, k is the wave number of the ocean wave, s g Is the slope of the group line.
By combining the distribution positions of wave field energy and group line energy, a two-dimensional binary self-adaptive filter h capable of filtering the group line energy in the wave number frequency spectrum and retaining the wave field energy is designed a (k, f), namely:
wherein s is g Is the slope of the energy distribution of the group line, k is the wave number of sea waves, g is the gravitational acceleration, f is the Doppler frequency,is a positional parameter of the wave field energy distribution in the wave number frequency spectrum.
Furthermore, the wave number frequency spectrum V after the filtering process in the step 4 a (k,f,β n ) The calculation method is as follows:
V a (k,f,β n )=V r (k,f,β n )·h a (7)
in the formula, h a Is an adaptive filter, V r (k,f,β n ) Is the wavenumber frequency spectrum.
For wave number frequency spectrum V after filtering a (k,f,β n ) Reconstructing space-time-dimension radial Doppler velocity sequence v by two-dimensional inverse fast Fourier transform a (r,t,β n ) The method comprises the following steps:
v a (r,t,β n )=2DIFFT[V a (k,f,β n )] (8)
where 2DIFFT represents a two-dimensional inverse fast fourier transform.
And, in the step 5, the space-dimensional velocity sequence is subjected to fast Fourier transform to obtain an energy spectrum G v (k,β n ) The specific calculation mode is as follows:
wherein t is begin 、t end And t number The number of points, respectively start time, end time and sample time, the FFT represents the fast fourier transform.
Estimating wavenumber direction spectrum S (k, β) from energy spectrum by transfer function TF n ):
S(k,β n )=TF·G v (k,β n ) (10)
Wherein, TF is a conversion function, and the specific calculation mode is as follows:
in the formula, tanh 2 And k is the wave number of the sea wave, d is the water depth, g is the gravitational acceleration, N is the number of space sampling points, deltak and Deltaβ are the wave number resolution and the angle resolution respectively, and theta is the antenna incident angle.
And, the six wave number direction spectrums are synthesized into a wave number-dimensional undirected wave spectrum S (k), namely:
wherein S (k, beta) n ) Is wave number direction spectrum, delta beta is angular resolution;
and then obtaining the undirected wave spectrum S (f) of the frequency dimension according to the conversion relation of the undirected wave spectrum of the wave number dimension and the frequency dimension, namely:
wherein f is the wave frequency and k is the wave number.
In addition, in the step 7, the effective wave height Hs and the average wave period Tav are calculated by using the first-order moment and the second-order moment, and the specific calculation method is as follows:
wherein S (f) is undirected wave spectrum, and f is wave frequency.
Compared with the prior art, the invention has the following advantages:
1) According to the positions of wave field energy and group line energy, the self-adaptive filter is designed to inhibit the phenomenon of 'group line' which leads to inaccurate wave parameter estimation, and the filter can be adaptively adjusted according to different ship speeds of ships, so that the wave parameter inversion of different ship speeds can be adapted, and the wave parameters of different navigation conditions under a ship-borne platform are met.
2) The method solves the problem that the broken wave and the ship motion affect the sea wave parameter inversion precision, not only maintains the advantages of the shore-based coherent microwave radar, but also makes up the defects of the shore-based coherent microwave radar in spatial ductility and time continuity, and greatly improves the robustness and the precision of sea wave parameter inversion of the ship-borne coherent S-band radar.
Drawings
FIG. 1 is a flow chart of an embodiment of the present invention.
Fig. 2 is a diagram showing a comparison of simulation results of a ship-borne coherent S-band radar ocean echo inversion undirected wave spectrum and a theoretical wave spectrum with a ship speed of 10 knots.
FIG. 3 is a graph of actual measurement results of the space-time inversion of the undirected wave spectrum and the buoy wave spectrum of the wave parameter by using the shipborne coherent S-band radar.
Fig. 4 is a graph showing the comparison of the actual measured results of the effective wave height and average wave period sequences of the inversion of the space-time inversion method of the wave parameters by using the shipborne coherent S-band radar.
Detailed Description
The invention provides a space-time inversion method for wave parameters by utilizing a shipborne coherent S-band radar, which solves the problem of influence of broken waves and ship motion on the inversion precision of the wave parameters, not only maintains the advantages of a shore-based coherent microwave radar, but also can make up the defects of the shore-based coherent microwave radar in spatial ductility and time continuity, and greatly improves the robustness and the precision of the shipborne coherent S-band radar wave parameter inversion.
The technical scheme of the invention is further described below with reference to the accompanying drawings and examples.
As shown in fig. 1, the invention provides a space-time inversion method for sea wave parameters by using a shipborne coherent S-band radar, which comprises the following steps:
Extracting average Doppler frequency offset of antenna echo Doppler spectrum in each direction by using a moment estimation method, searching ocean echo space-time dimension Doppler spectrum peak through spectrum peak, and searching left and right frequency offset f at 18dB of spectrum peak l And f r Finally, extracting average Doppler frequency offset f by a moment estimation method d The method comprises the following steps:
where f is the Doppler frequency and σ is the Doppler spectrum.
Space-time dimension radial Doppler velocity sequence v (r, t, beta) n ) From the average Doppler frequency offset f d The specific calculation mode is as follows:
wherein lambda is the electromagnetic wavelength, beta n Represents the nth th The pointing direction of the root antenna, r is the radial distance, and t is the sampling time.
Converting 2-dimensional velocity sequences into wavenumber frequency spectra V using 2-dimensional fast Fourier transforms (2-D FFT) r (k,f,β n ) The method comprises the following steps:
V r (k,f,β n )=2DFFT[v(r,t,β n )] (3)
where k is the wave number, f is the Doppler frequency, beta n Represents the nth th The pointing direction of the root antenna, r is the radial distance, and t is the sampling time.
Step 3, estimating the position parameters of wave field energy distribution in the wave number frequency spectrum by adopting a least square fitting algorithmAnd slope parameter s of group line energy distribution g And combining the two parameters to design a two-dimensional binary adaptive filter h for filtering out group line energy in the wave number frequency spectrum and retaining wave field energy a (k,f)。
The least squares fitting algorithm finds the best functional match for the data by minimizing the sum of squares of the errors. The energy distribution of the field due to the influence of the forward speed of the ship is no longer along the first order wave dispersion relation (omega 2 =gk) estimating the position parameters of the wave field energy distribution in the wave number frequency spectrum using least squares fitting algorithmIt is estimated along the following formula:
wherein ω is the angular frequency of the ocean wave, k is the wave number of the ocean wave,g is the gravity acceleration, which is the included angle between the pointing direction of the antenna and the forward motion direction of the ship.
Estimating slope parameter s of group line energy distribution by least square fitting algorithm g Fitting the distribution of group line energy, namely:
ω=s g k (5)
wherein ω is the angular frequency of the ocean wave, k is the wave number of the ocean wave, s g Is the slope of the group line.
By combining the distribution positions of wave field energy and group line energy, a two-dimensional binary self-adaptive filter h capable of filtering the group line energy in the wave number frequency spectrum and retaining the wave field energy is designed a (k, f), namely:
wherein s is g The slope of the energy distribution of the group line is that k is wave number of sea waves and g is gravity acceleration.
Wave number frequency spectrum V after filtering a (k,f,β n ) The calculation method is as follows:
V a (k,f,β n )=V r (k,f,β n )·h a (7)
in the formula, h a Is self-containedAdaptive filter, V r (k,f,β n ) Is the wavenumber frequency spectrum.
For wave number frequency spectrum V after filtering a (k,f,β n ) Reconstructing a space-time-dimensional radial Doppler velocity sequence v by performing two-dimensional inverse fast Fourier transform (2-D FFT) a (r,t,β n ) The method comprises the following steps:
v a (r,t,β n )=2DIFFT[V a (k,f,β n )] (8)
Performing fast Fourier transform on the space-dimensional velocity sequence to obtain an energy spectrum G v (k,β n ) The specific calculation mode is as follows:
wherein t is begin 、t end And t number The number of points is the start time, the end time and the sampling time, respectively.
Estimating wavenumber direction spectrum S (k, beta) from energy spectrum through conversion relation TF n ):
S(k,β n )=TF·G v (k,β n ) (10)
Wherein, TF is a conversion function, and the specific calculation mode is as follows:
in the formula, tanh 2 And k is the wave number of the sea wave, d is the water depth, g is the gravitational acceleration, N is the number of space sampling points, deltak and Deltaβ are the wave number resolution and the angle resolution respectively, and theta is the antenna incident angle.
And 6, obtaining radar ocean echoes in different directions by six shipborne radar antennas pointing to different sea surfaces, combining the radar ocean echoes in different directions, repeatedly executing the steps 1-5, and synthesizing the six wave number direction spectrums into an undirected wave spectrum S (f), namely a wave height spectrum.
Because the directional resolution of the antennas is 30 degrees, each antenna can receive sea surface information of 30 degrees in the forward and reverse directions of the sea surface, and therefore six shipborne radar antennas pointing to different sea surfaces can obtain 360 degrees of sea surface information, and omni-directional coverage of the sea surface information is achieved.
Firstly, synthesizing a wave number-dimensional undirected wave spectrum S (k) from six wave number direction spectrums, namely:
wherein S (k, beta) n ) Is the wavenumber direction spectrum, Δβ is the angular resolution.
And then obtaining the undirected wave spectrum S (f) of the frequency dimension according to the conversion relation of the undirected wave spectrum of the wave number dimension and the frequency dimension, namely:
wherein f is the wave frequency and k is the wave number.
And 7, calculating the effective wave height Hs and the average wave period Tav through the first-order moment and the second-order moment by combining the wave theory.
The effective wave height Hs and the average wave period Tav are calculated by using the first-order moment and the second-order moment, and the specific calculation mode is as follows:
wherein S (f) is undirected wave spectrum, and f is wave frequency.
The wind speed value is 12m/s, the wind area is 100km, the wind direction is 0 degree, the ship speed is 10 knots, the undirected sea wave spectrum type is a JONSWAP spectrum, the directed sea wave spectrum is a Longuet-Higgins spectrum, and a simulated sea echo Doppler spectrum is generated, and the result is shown in figure 2. As can be seen from FIG. 2, the undirected wave spectrum obtained by inversion of the method provided by the invention has very high coincidence degree with the theoretical undirected wave spectrum, and the effective wave height (1.97 m) and the average wave period (5.34 s) of the inversion are very close to the theoretical values (1.96 m and 5.39 s).
The specific implementation mode of the invention is applied to the sea wave parameter inversion of the actual measurement shipborne coherent S-band radar, and the effectiveness of the invention is verified by extracting the effective wave height and the average wave period from the undirected sea wave spectrum and comparing the effective wave height and the wave period of the buoy. The result pair of the undirected wave spectrum of the actual measurement buoy and the undirected wave spectrum obtained by inversion of the ocean wave parameter space-time inversion method of the actual measurement shipborne coherent S-band radar at a certain moment is shown in fig. 3, and as can be seen from fig. 3, the period of extracting effective wave height and average wave from the undirected wave spectrum inverted by the method provided by the invention is close to the period of effective wave height and average wave measured from the buoy. The effective wave height and average wave period of the 270S-band radar inversion for approximately 2 days was compared to the effective wave height and average wave period acquired by the buoy, as shown in fig. 4. The effective wave height and average wave period extracted by the method are well matched with the buoy, the root mean square error is 0.30m and 0.35s respectively, and the average error is 0.24m and 0.28s respectively.
The specific embodiments described herein are offered by way of example only to illustrate the spirit of the invention. Those skilled in the art may make various modifications or additions to the described embodiments or substitutions thereof without departing from the spirit of the invention or exceeding the scope of the invention as defined in the accompanying claims.
Claims (10)
1. A sea wave parameter space-time inversion method using a ship-borne coherent S-band radar is characterized by comprising the following steps:
step 1, original echo data acquired by a ship-borne coherent S-band radar sequentially pass through first FFT distance conversion processing and second FFT Doppler conversion processing to obtain to-be-processed echo dataThe space-time dimension Doppler spectrum of the ocean echo is extracted, and then the average Doppler frequency offset f of the antenna echo Doppler spectrum in all directions is extracted by a moment estimation method d ;
Step 2, by average Doppler frequency offset f d Estimating a space-time dimensional radial Doppler velocity sequence v (r, t, beta) n ) And 2-dimensional velocity sequences v (r, t, beta) are transformed using a 2-dimensional fast fourier transform n ) Conversion to wave number frequency spectrum V r (k,f,β n );
Step 3, estimating the position parameters of wave field energy distribution in the wave number frequency spectrum by adopting a least square fitting algorithmAnd slope parameter s of group line energy distribution g And combining the two parameters to design a two-dimensional binary adaptive filter h a (k,f);
Step 4, for wave number frequency spectrum V r (k,f,β n ) Filtering to obtain wave number frequency spectrum V after filtering a (k,f,β n ) For the wave number frequency spectrum V after filtering a (k,f,β n ) Reconstructing space-time-dimension radial Doppler velocity sequence v by two-dimensional inverse fast Fourier transform a (r,t,β n );
Step 5, performing fast Fourier transform on the space-dimensional radial Doppler velocity sequence to obtain an energy spectrum G v (k,β n ) Estimating wave number direction spectrum S (k, β) from energy spectrum by conversion function n );
Step 6, six shipborne radar antennas pointing to different sea surfaces are used for obtaining radar ocean echoes in different directions, the radar ocean echoes in different directions are combined, the steps 1-5 are repeatedly executed, and the six wave number direction spectrums are synthesized into an undirected wave spectrum S (f);
and 7, calculating the effective wave height Hs and the average wave period Tav through the first-order moment and the second-order moment.
2. A method of space-time inversion of sea wave parameters using a coherent S-band radar on board a vessel as claimed in claim 1, wherein: in step 1, moment estimation method is used for extractionThe average Doppler frequency offset of the Doppler spectrum of the antenna echo in each direction is firstly searched for the Doppler spectrum peak of the ocean echo space-time dimension through the spectrum peak, and then the left and right frequency offset f at 18dB of the spectrum peak is searched for l And f r Finally, extracting average Doppler frequency offset f by a moment estimation method d The method comprises the following steps:
where f is the Doppler frequency and σ is the Doppler spectrum.
3. A method of space-time inversion of sea wave parameters using a coherent S-band radar on board a vessel as claimed in claim 1, wherein: step 2 hollow time-dimensional radial Doppler velocity sequence v (r, t, beta) n ) From the average Doppler frequency offset f d The specific calculation mode is as follows:
wherein lambda is the electromagnetic wavelength, beta n Represents the nth th The pointing direction of the root antenna, r is the radial distance, and t is the sampling time;
2-dimensional velocity sequence v (r, t, beta) using 2-dimensional fast fourier transform n ) Conversion to wave number frequency spectrum V r (k,f,β n ) The method comprises the following steps:
V r (k,f,β n )=2DFFT[v(r,t,β n )] (3)
where k is the wave number, f is the Doppler frequency, beta n Represents the nth th The pointing direction of the root antenna, r is the radial distance, t is the sampling time, and 2DFFT represents the two-dimensional fast fourier transform.
4. A method of space-time inversion of sea wave parameters using a coherent S-band radar on board a vessel as claimed in claim 1, wherein: least square fitting algorithm in step 3Finding the best function match of the data by minimizing the sum of squares of the errors, the wave field energy distribution no longer follows the first order wave dispersion relationship ω due to the influence of the forward speed of the vessel 2 =gk, estimation of the position parameters of the wave field energy distribution in the wave number frequency spectrum using least squares fitting algorithmIt is estimated along the following formula:
5. The space-time inversion method for sea wave parameters by using shipborne coherent S-band radar according to claim 4, wherein: in step 3, the slope parameter s of the energy distribution of the group line is estimated by utilizing a least square fitting algorithm g Fitting the distribution of group line energy, namely:
ω=s g k (5)
wherein ω is the angular frequency of the ocean wave, k is the wave number of the ocean wave, s g Is the slope of the group line;
by combining the distribution positions of wave field energy and group line energy, a two-dimensional binary self-adaptive filter h capable of filtering the group line energy in the wave number frequency spectrum and retaining the wave field energy is designed a (k, f), namely:
6. A method of space-time inversion of sea wave parameters using a coherent S-band radar on board a vessel as claimed in claim 1, wherein: wave number frequency spectrum V after filtering in step 4 a (k,f,β n ) The calculation method is as follows:
V a (k,f,β n )=V r (k,f,β n )·h a (7)
in the formula, h a Is an adaptive filter, V r (k,f,β n ) Is wave number frequency spectrum;
for wave number frequency spectrum V after filtering a (k,f,β n ) Reconstructing space-time-dimension radial Doppler velocity sequence v by two-dimensional inverse fast Fourier transform a (r,t,β n ) The method comprises the following steps:
v a (r,t,β n )=2DIFFT[V a (k,f,β n )] (8)
where 2DIFFT represents a two-dimensional inverse fast fourier transform.
7. A method of space-time inversion of sea wave parameters using a coherent S-band radar on board a vessel as claimed in claim 1, wherein: in step 5, performing fast Fourier transform on the space-dimensional velocity sequence to obtain an energy spectrum G v (k,β n ) The specific calculation mode is as follows:
wherein t is begin 、t end And t number The number of points, respectively start time, end time and sample time, the FFT represents the fast fourier transform.
8. The space-time inversion method for sea wave parameters by using shipborne coherent S-band radar according to claim 7, wherein: estimating wave number direction spectrum S (k, beta) from energy spectrum by transfer function TF in step 5 n ):
S(k,β n )=TF·G v (k,β n ) (10)
Wherein, TF is a conversion function, and the specific calculation mode is as follows:
in the formula, tanh 2 And k is the wave number of the sea wave, d is the water depth, g is the gravitational acceleration, N is the number of space sampling points, deltak and Deltaβ are the wave number resolution and the angle resolution respectively, and theta is the antenna incident angle.
9. A method of space-time inversion of sea wave parameters using a coherent S-band radar on board a vessel as claimed in claim 1, wherein: in the step 6, firstly, synthesizing a wave number dimensional undirected wave spectrum S (k) from the six wave number direction spectrums, namely:
wherein S (k, beta) n ) Is wave number direction spectrum, delta beta is angular resolution;
and then obtaining the undirected wave spectrum S (f) of the frequency dimension according to the conversion relation of the undirected wave spectrum of the wave number dimension and the frequency dimension, namely:
wherein f is the wave frequency and k is the wave number.
10. A method of space-time inversion of sea wave parameters using a coherent S-band radar on board a vessel as claimed in claim 1, wherein: in the step 7, the effective wave height Hs and the average wave period Tav are calculated by using the first-order moment and the second-order moment, and the specific calculation mode is as follows:
wherein S (f) is undirected wave spectrum, and f is wave frequency.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310265278.0A CN116430337A (en) | 2023-03-15 | 2023-03-15 | Sea wave parameter space-time inversion method using shipborne coherent S-band radar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310265278.0A CN116430337A (en) | 2023-03-15 | 2023-03-15 | Sea wave parameter space-time inversion method using shipborne coherent S-band radar |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116430337A true CN116430337A (en) | 2023-07-14 |
Family
ID=87084630
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310265278.0A Pending CN116430337A (en) | 2023-03-15 | 2023-03-15 | Sea wave parameter space-time inversion method using shipborne coherent S-band radar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116430337A (en) |
-
2023
- 2023-03-15 CN CN202310265278.0A patent/CN116430337A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111257886B (en) | Method for inverting sea wave parameters by using single shipborne X-band radar image | |
CN103383449B (en) | ESPRIT algorithm based short-range clutter suppression method for airborne radar | |
CN111006743B (en) | Method for eliminating interference in water level measurement based on plane radar water level gauge | |
CN111580048B (en) | Broadband sound source depth estimation method using single-vector hydrophone | |
CN110286370B (en) | Motion compensation method for shipborne ground wave radar in sailing state | |
CN112987004B (en) | Water surface and underwater target classification method based on horizontal array in shallow sea environment | |
CN113109817B (en) | Vector hydrophone deployment depth estimation method | |
CN109765521B (en) | Beam domain imaging method based on subarray division | |
CN107607936B (en) | High-frequency sky-ground wave radar ocean surface flow inversion method | |
JP5336292B2 (en) | Radar wave analyzer | |
CN103529437A (en) | Method used for captive-balloon-borne phased array radar to distinguish open space targets under multi-target condition | |
CN113011006B (en) | Target depth estimation method based on cross-correlation function pulse waveform matching | |
CN109061640B (en) | Azimuth fuzzy suppression method for forward-orbit interference SAR ocean current inversion | |
CN114167423A (en) | Radar sea wave parameter measuring method based on depth regression network | |
CN107783137B (en) | Sound Doppler and acoustic correlation velocity measurement method based on five-beam configuration | |
CN109884337B (en) | Method for detecting sea surface wind direction by using high-frequency ground wave radar | |
CN111487628A (en) | 'zero degree' interference suppression method for downward-looking multi-beam synthetic aperture imaging sonar | |
CN117169882A (en) | Shipborne radar sea wave information inversion method | |
CN116430337A (en) | Sea wave parameter space-time inversion method using shipborne coherent S-band radar | |
Lyons et al. | Modeling the effect of random roughness on synthetic aperture sonar image statistics | |
CN113466821B (en) | Wave direction inversion method for ship-borne coherent microwave radar | |
CN107422320B (en) | Method for eliminating influence of rainfall on X-band radar observation sea waves | |
CN115902791A (en) | S-band wave-measuring radar time Doppler spectrum-based sea wave inversion method and system | |
CN113589283B (en) | Extraction method of ship kelvin wake elevation based on spaceborne interference imaging altimeter | |
CN115032601A (en) | Marine radar target detection algorithm for inhibiting sea clutter in image sequence based on space-time combined filtering technology |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |