CN109283495B - X-band radar ocean current inversion method based on cross spectrum analysis - Google Patents
X-band radar ocean current inversion method based on cross spectrum analysis Download PDFInfo
- Publication number
- CN109283495B CN109283495B CN201811035048.0A CN201811035048A CN109283495B CN 109283495 B CN109283495 B CN 109283495B CN 201811035048 A CN201811035048 A CN 201811035048A CN 109283495 B CN109283495 B CN 109283495B
- Authority
- CN
- China
- Prior art keywords
- wave
- phase
- ocean current
- band radar
- radar
- 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.)
- Active
Links
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
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
-
- 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
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides an X-band radar ocean current inversion method based on cross spectrum analysis, which comprises the steps of firstly, respectively carrying out cross spectrum analysis on two adjacent images in an X-band radar image to obtain a coherence coefficient spectrum and a phase spectrum; then, respectively averaging the coherence coefficient spectrum and the phase spectrum, determining the dominant wave number according to the peak value of the coherence coefficient spectrum, and eliminating the direction ambiguity of the dominant wave direction by using the average phase; and selecting different wave directions, establishing a model according to the phase velocity obtained by the frequency dispersion relation, and solving by using a least square method to obtain the ocean current vector. The method can solve the problems that the X-band radar has low accuracy in observing the sea surface flow field and a set of coherent X-band radar cannot obtain a vector flow field in the prior art, so as to meet the requirements of business observation of offshore marine environment, marine production and life service and the like.
Description
Technical Field
The invention belongs to the technical field of ocean remote sensing, and particularly relates to an X-band radar ocean current inversion method based on cross spectrum analysis.
Background
The sea surface flow field is an important ocean dynamic parameter, which not only influences activities such as ocean engineering, ship navigation, ocean cultivation and entertainment, but also is closely related to the occurrence, development and the like of ocean disasters such as ocean oil spill pollution, red tide and the like which are increasingly serious in recent years, so that the ocean current observation has important significance. Current meters and Acoustic Doppler Current Profilers (ADCPs) are widely used for the observation of sea currents, but they can only obtain point or line sea current changes, are difficult and expensive to install and maintain, and may not be used especially in the complex sea areas or extreme sea conditions of interest. Satellite remote sensing can obtain ocean current of a large area, but a satellite-borne radar altimeter can only obtain ground-turn current, a synthetic aperture radar can only measure radial flow velocity, time and space resolution are poor, and a satellite has a long repeated access period, so that the ocean current observation requirement of offshore business cannot be met. The X-band radar has high time and space resolution, can observe the sea surface all day long and all weather, and is widely applied to the observation of ocean currents in recent years.
At present, mainly shipborne navigation X-band radar and coherent X-band radar are used for offshore ocean current observation. The algorithm for inverting the ocean current by using the non-coherent navigation X-band radar image is mainly based on the three-dimensional Fourier transform algorithm proposed by Young et al (1985)[1-4]However, it requires the conversion of the radar image spectrum into a sea wave spectrum using an empirical modulation transfer function, which is associated with the radar system,Sea conditions, installation conditions and other factors are related and difficult to accurately determine, and the form difference of modulation transfer functions obtained by using different experimental data and simulation methods is large[4,5]. The coherent radar can record the strength and phase information of sea echo, obtain the radial velocity of sea from the phase of echo according to Doppler effect, but the radial velocity is related to ocean current, wave motion and the like, so that ocean current is not easy to separate, and more than two sets of radars are required to synchronously observe to obtain vector flow field[2]. In addition, due to the limitation of a radar system, the two radars have the problems of aliasing of time sampling, attenuation of radar images along with distance and direction and the like, and the factors limit the observation accuracy of ocean current. Therefore, the invention provides a simple and feasible method for inverting a sea surface flow field by using an X-band radar, and aims to solve the problems that in the prior art, coherent and non-coherent X-band radars have poor accuracy in observing the sea flow, a vector flow field cannot be obtained by using only one set of coherent radars and the like.
Reference documents:
[1]I.R.Young,W.Rosenthal,and F.Ziemer,"A three-dimensional analysis of marine radar images for the determination of ocean wave directionality and surface currents,"Journal of Geophysical Research,vol.90,pp.1049-1059,1985.
[2] trizna Dennis, Method and apparatus for coherent marine radar measurements of properties of ocean waves and currents, application No.: US20100868912, U.S. patent, filing date: 2010.8.26
[3] Wu stambin et al, an X-band wave-measuring radar ocean current inversion preprocessing method, application No.: CN201310029849.7, china, application date: 2013.1.25
[4]J.C.Nieto-Borge,G.R.Rodriguez,K.Hessner,and P.I.Gonzalez,"Inversion of Marine Radar Images for Surface Wave Analysis,"Journal of Atmospheric and Oceanic Technology,vol.21,pp.1291-1300,2004.
[5]Zhongbiao Chen,Biao Zhang,Yijun He,Zhongfeng Qiu,Perrie William,A New Modulation Transfer Function for Ocean Wave Spectra Retrieval from X-band Marine Radar Imagery,Chinese Journal of Oceanology and Limnology.Vol.33No.5,P.1132-1141,2015.
Disclosure of Invention
The purpose of the invention is as follows: the invention provides a simple and feasible method for inverting a sea surface flow field by using an X-band radar image, and aims to solve the problems that in the prior art, the accuracy of an X-band radar for observing sea current is low, and a set of coherent X-band radar cannot obtain a vector flow field.
The technical scheme is as follows: the invention provides an X-band radar ocean current inversion method based on cross spectrum analysis, which specifically comprises the following steps:
step 1: selecting an X-waveband radar image sequence I, and respectively performing cross spectrum analysis on two adjacent images to obtain a coherence coefficient
Φij(k,φ)=arg(Sij),i,j=1,2,…,n (2)
Wherein Ii(x, y) are intensity values of the radar image, x and y are abscissa and ordinate of the center point of the image, i and j are the sequence of images, n is the number of images, k andforming frequency domain coordinates, K is the wave number,is the wave direction, SijIs an image Ii(x, y) and IjPower spectrum of (x, y), SiiIs Ii(x, y) and its own power spectrum, and arg () represents the argument of a complex number.
Step 2: coherence coefficient according to step 1And phasePerforming cross spectrum analysis on all adjacent images to obtain average coherence coefficient
And step 3: determining the dominant wavenumber k according to the peak value of the coherence coefficient in the step 1mDetermining a real main wave number by using the average phase in the step 2 for eliminating the direction ambiguity of 180 degrees; the 180 degree directional blur is the two symmetrical maxima of the average coherence coefficient in equation 3 about the far point.
And 4, step 4: selecting a main wave direction phimAll wave directions phi of different directions nearbypP is the wave direction phimThe number of all nearby wave directions in different directions is determined, and a model is established;
wherein, cradarIs wave directionThe phase velocity of the directional wave,ctheoryobtaining a theoretical value of the phase velocity according to the dispersion relation theory of sea wavesWherein g is the acceleration of gravity; h is the depth of the water,the time interval between every two adjacent images is the same;is a wave number vector; q ═ 1, 2, 3.
Further, the method for determining the true dominant wave number in step 3 to eliminate the 180-degree directional blur includes: select to satisfyPeak wave number ofIs the true dominant wavenumber.
Further, the wave direction phi selected in the step 4pThe wave direction phipSimultaneously satisfies the following requirements:
1,φpthe coefficient of coherence gamma (k, phi)p) > 0.6, wherein the wavenumber k satisfies: k is more than or equal to 0.03 and less than or equal to 0.3 tad/m;
Has the advantages that: the method is simple and easy to implement, and can solve the problems that the X-band radar has low accuracy in observing the sea surface flow field and a set of coherent X-band radar cannot obtain a vector flow field in the prior art so as to meet the requirements of business observation, marine production and life service and the like of the offshore marine environment.
Drawings
FIG. 1 is a flow chart of the algorithm of the present invention;
FIG. 2 is a comparison of the wave phase velocity obtained from an X-band radar image in accordance with the present invention and the theoretical phase velocity.
Detailed Description
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate an embodiment of the invention and, together with the description, serve to explain the invention and not to limit the invention.
As shown in fig. 1, the present invention provides an X-band radar ocean current inversion method based on cross-spectrum analysis, which specifically includes the following steps:
step 1: selecting an X-waveband radar image sequence I, and respectively performing cross spectrum analysis on two adjacent images to obtain a coherence coefficient
Wherein Ii(x, y) are intensity values of the radar image, x and y are abscissa and ordinate of the center point of the image, i and j are the sequence of images, n is the number of images, k andforming frequency domain coordinates, K is the wave number,is the wave direction, SijIs an image Ii(x, y) and IjPower spectrum of (x, y), SiiIs Ii(x, y) power spectrum with itself, arg () representing the argument of the complex number;
step 2: coherence coefficient according to step 1And phasePerforming cross spectrum analysis on all adjacent images to obtain an average coherence coefficient
And step 3: determining the dominant wavenumber k according to the peak value of the coherence coefficient in the step 1mDetermining a real main wave number by using the average phase in the step 2 for eliminating the direction ambiguity of 180 degrees;
and 4, step 4: selecting a main wave direction phimAll wave directions phi of different directions nearbypP is the wave direction phimThe number of all nearby wave directions in different directions is determined, and a model is established;
wherein, cradarIs wave directionThe phase velocity of the directional wave,ctheoryobtaining a theoretical value of the phase velocity according to the dispersion relation theory of sea wavesWherein g is the acceleration of gravity; h is the depth of the water,the time interval between every two adjacent images is the same;is a wave number vector; q ═ 1, 2, 3.
It should be noted that, in the step 3, determining the true dominant wave number is used to eliminate the directional ambiguity of 180 degrees, and the method for selecting the true dominant wave number is as follows: select to satisfyPeak wave number ofIs the true dominant wavenumber.
The wave direction phi selected in the step 5pThe wave direction phipAt the same time satisfyThe following requirements are set forth:
1,φpthe coefficient of coherence gamma (k, phi)p) > 0.6, wherein the wavenumber k satisfies: k is more than or equal to 0.03 and less than or equal to 0.3 tad/m;
As shown in fig. 2, the comparison between the wave phase velocity obtained by using an X-band radar image and the theoretical phase velocity shows that the two have better consistency, so that a vector flow field can be obtained by using least square fitting.
Claims (3)
1. An X-band radar ocean current inversion method based on cross spectrum analysis is characterized by comprising the following steps:
step 1: selecting an X-waveband radar image sequence I, and respectively performing cross spectrum analysis on two adjacent images to obtain a coherence coefficient
Φij(k,φ)=arg(Sij),i,j=1,2,…,n (2)
Wherein Ii(x, y) are intensity values of the radar image, x and y are abscissa and ordinate of the center point of the image, i and j are the sequence of images, n is the number of images, k andforming frequency domain coordinates, K is the wave number,is the wave direction, SijIs an image Ii(x, y) and IjPower spectrum of (x, y), SiiIs Ii(x, y) power spectrum with itself, arg () representing the argument of the complex number;
step 2: coherence factor according to step 1And phasePerforming cross spectrum analysis on all adjacent images to obtain average coherence coefficient
And step 3: determining the dominant wavenumber k according to the peak value of the coherence coefficient in the step 1mDetermining a true main wave number by using the average phase in the step 2 for eliminating 180-degree direction ambiguity;
and 4, step 4: selecting a main wave direction phimAll wave directions phi of different directions nearbypP is the wave direction phimThe number of all nearby wave directions in different directions is determined, and a model is established;
wherein, cradarIs wave directionThe phase velocity of the directional wave,ctheoryobtaining a theoretical value of the phase velocity according to the dispersion relation theory of sea wavesWherein g is the acceleration of gravity; h is the depth of the water,the time interval between every two adjacent images is the same;is a wave number vector; q ═ 1, 2, 3, …, p.
2. The method for X-band radar ocean current inversion based on cross-spectral analysis according to claim 1, wherein the true dominant wavenumber determined in step 3 is used for eliminating 180-degree directional ambiguity, and the method for selecting the true dominant wavenumber is as follows: select to satisfyPeak wave number ofIs the true dominant wavenumber.
3. The method for X-band radar ocean current inversion based on cross-spectral analysis according to claim 1, wherein the wave direction φ selected in the step 4pThe wave direction phipSimultaneously satisfies the following requirements:
1,φpthe coefficient of coherence gamma (k, phi)p) > 0.6, wherein the wavenumber k satisfies: k is more than or equal to 0.03 and less than or equal to 0.3 rad/m;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811035048.0A CN109283495B (en) | 2018-09-06 | 2018-09-06 | X-band radar ocean current inversion method based on cross spectrum analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811035048.0A CN109283495B (en) | 2018-09-06 | 2018-09-06 | X-band radar ocean current inversion method based on cross spectrum analysis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109283495A CN109283495A (en) | 2019-01-29 |
CN109283495B true CN109283495B (en) | 2022-05-06 |
Family
ID=65183533
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811035048.0A Active CN109283495B (en) | 2018-09-06 | 2018-09-06 | X-band radar ocean current inversion method based on cross spectrum analysis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109283495B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111102120A (en) * | 2019-11-25 | 2020-05-05 | 国家海洋技术中心 | Method for testing power characteristics of wave energy power generation equipment |
CN113050091B (en) * | 2021-03-08 | 2021-09-28 | 国家海洋技术中心 | Wind speed and direction joint inversion method and system for satellite-borne synthetic aperture radar |
CN113064161B (en) * | 2021-03-30 | 2023-05-30 | 南京信息工程大学 | Wave spectrometer cross spectrum calculation method based on double sub-pulse reconstruction |
CN113176572B (en) * | 2021-04-27 | 2023-10-20 | 中山大学 | Sea surface wave spectrum inversion method and system based on circular SAR |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813476A (en) * | 2010-03-19 | 2010-08-25 | 天津大学 | Three-dimensional real-time monitoring system for offshore wave parameters |
CN103604944A (en) * | 2013-12-11 | 2014-02-26 | 哈尔滨工业大学 | Surface flow measurement method based on monostation shipborne high-frequency ground wave radar |
CN104101864A (en) * | 2013-04-10 | 2014-10-15 | 南京信息工程大学 | Navigation X-waveband radar ocean wave parameter inversion algorithm based on EOF decomposition |
CN105445730A (en) * | 2015-11-27 | 2016-03-30 | 南京信息工程大学 | Ocean current field inversion satellite-borne SAR system based on angle diversity, and method thereof |
CN106093936A (en) * | 2016-08-29 | 2016-11-09 | 中船重工鹏力(南京)大气海洋信息系统有限公司 | Slowly the unrestrained stream information extracting method under pattern is swept based on coherent radar |
CN106990404A (en) * | 2017-03-30 | 2017-07-28 | 南京信息工程大学 | A kind of autoscale algorithm using X-band radar inverting sea wave height of navigating |
-
2018
- 2018-09-06 CN CN201811035048.0A patent/CN109283495B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813476A (en) * | 2010-03-19 | 2010-08-25 | 天津大学 | Three-dimensional real-time monitoring system for offshore wave parameters |
CN104101864A (en) * | 2013-04-10 | 2014-10-15 | 南京信息工程大学 | Navigation X-waveband radar ocean wave parameter inversion algorithm based on EOF decomposition |
CN103604944A (en) * | 2013-12-11 | 2014-02-26 | 哈尔滨工业大学 | Surface flow measurement method based on monostation shipborne high-frequency ground wave radar |
CN105445730A (en) * | 2015-11-27 | 2016-03-30 | 南京信息工程大学 | Ocean current field inversion satellite-borne SAR system based on angle diversity, and method thereof |
CN106093936A (en) * | 2016-08-29 | 2016-11-09 | 中船重工鹏力(南京)大气海洋信息系统有限公司 | Slowly the unrestrained stream information extracting method under pattern is swept based on coherent radar |
CN106990404A (en) * | 2017-03-30 | 2017-07-28 | 南京信息工程大学 | A kind of autoscale algorithm using X-band radar inverting sea wave height of navigating |
Also Published As
Publication number | Publication date |
---|---|
CN109283495A (en) | 2019-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109283495B (en) | X-band radar ocean current inversion method based on cross spectrum analysis | |
Bell | Shallow water bathymetry derived from an analysis of X-band marine radar images of waves | |
Senet et al. | Determination of bathymetric and current maps by the method DiSC based on the analysis of nautical X-band radar image sequences of the sea surface (November 2007) | |
CN106990404B (en) | Automatic scaling algorithm for inverting sea wave height by using navigation X-band radar | |
CN103941257B (en) | A kind of method of the pathfinder image inverting wind direction of ocean surface based on wave number energy spectrum | |
CN111257886B (en) | Method for inverting sea wave parameters by using single shipborne X-band radar image | |
CN110456348B (en) | Wave truncation wavelength compensation method for multi-view-direction SAR wave spectrum data fusion | |
CN103293521B (en) | Method for detecting water depth of offshore sea by X-band radar | |
Lund et al. | Analysis of internal wave signatures in marine radar data | |
Wyatt | Measuring the ocean wave directional spectrum ‘First Five’with HF radar | |
Chernyshov et al. | Rapid wavelet-based bathymetry inversion method for nearshore X-band radars | |
CN109085556B (en) | High-frequency ground wave radar wave field forming method based on first-order and second-order peak ratios | |
Lyons et al. | Modeling the effect of seafloor ripples on synthetic aperture sonar speckle statistics | |
Gangeskar | Verifying high-accuracy ocean surface current measurements by X-band radar for fixed and moving installations | |
Campana et al. | A new inversion method to obtain upper-ocean current-depth profiles using X-band observations of deep-water waves | |
CN102538768A (en) | Method for measuring water depth of shallow sea based on double-frequency high-frequency ground wave radar | |
CN114355306A (en) | Method for inverting wave height by using small-aperture high-frequency radar ocean echo | |
Chen et al. | Application of different internal solitary wave theories for SAR remote sensing inversion in the northern South China Sea | |
RU2466426C1 (en) | Method of reconstructing sea-floor relief when measuring depth using hydroacoustic apparatus | |
CN112612027B (en) | Ocean internal wave monitoring method utilizing sound energy fluctuation in shallow sea environment | |
CN112162282A (en) | Synthetic aperture radar-based sea surface flow velocity inversion method | |
US20210033725A1 (en) | Radar image processing | |
Zuckerman et al. | Bathymetry and water-level estimation using x-band radar at a tidal inlet | |
McCann et al. | A simple offset “calibration” method for the accurate geographic registration of ship-borne X-band radar intensity imagery | |
CN110673128B (en) | X-waveband shore-based radar flow measurement method based on intermittent up-down frequency modulation waves |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 210032 No. 219 Ning six road, Jiangbei new district, Nanjing, Jiangsu Applicant after: Nanjing University of Information Science and Technology Address before: 211500 Yuting Square, 59 Wangqiao Road, Liuhe District, Nanjing City, Jiangsu Province Applicant before: Nanjing University of Information Science and Technology |
|
GR01 | Patent grant | ||
GR01 | Patent grant |