CN113064129B - High-frequency ground wave radar ocean current synthesis method - Google Patents
High-frequency ground wave radar ocean current synthesis method Download PDFInfo
- Publication number
- CN113064129B CN113064129B CN202110236049.7A CN202110236049A CN113064129B CN 113064129 B CN113064129 B CN 113064129B CN 202110236049 A CN202110236049 A CN 202110236049A CN 113064129 B CN113064129 B CN 113064129B
- Authority
- CN
- China
- Prior art keywords
- current
- grid
- grids
- alternative
- ocean current
- 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
- 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
- G01S7/415—Identification of targets based on measurements of movement associated with the target
-
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/62—Sense-of-movement determination
-
- 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)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
The invention provides a high-frequency ground wave radar ocean current synthesis method, which is characterized in that the overlapped area covered by double stations is considered when ocean current velocity vector synthesis is carried out, and the output of a meshed flow field is completed through a series of interpolation means; for abnormal data and isolated measuring points, abnormal values are removed through an effective quality control means, and vacant data points are filled through an interpolation algorithm, so that full-area ocean current synthesis of an overlapped area is realized.
Description
Technical Field
The invention belongs to the technical field of high-frequency ground wave radars, and particularly relates to a sea current synthesis method of a high-frequency ground wave radar.
Background
The high-frequency ground wave radar realizes beyond-line-of-sight monitoring on ocean dynamic parameters such as wind, wave, flow and the like in ocean environment states by utilizing the characteristic that energy attenuation is small when the vertical high-frequency electromagnetic wave propagates along the surface of the high-conductivity sea water, and has the characteristics of long observation distance, large coverage area, multiple inversion elements, all-weather operation and the like, wherein the monitoring of ocean currents realizes conventional business operation in a plurality of areas. The current is basically that two remote radar stations are erected, and the actual current direction and magnitude are calculated through the projection quantity of the two directions. Because the single stations are all monitoring data under the polar coordinate system, the overlapping area covered by the double stations needs to be considered during synthesis, and the output of the grid flow field is completed through a series of interpolation means.
When the high-frequency ground wave radar inverts ocean current data, abnormal data and isolated measuring points can be generated under the conditions of echo distance, environmental interference and self equipment, so that effective quality control means are needed to eliminate abnormal values so as to fill up vacant data points.
Disclosure of Invention
In order to solve the problems, the invention provides a high-frequency ground wave radar ocean current synthesis method which can correct abnormal data and ocean current velocity vectors of isolated measuring points and realize the ocean current synthesis of a whole area.
A high-frequency ground wave radar ocean current synthesis method comprises the following steps:
s1: dividing grids in an overlapping area between the sector detection areas covered by the two radar stations, screening grids with actual measurement data sampling points of the radar stations in a set neighborhood range as effective grids, and respectively taking each effective grid as a current grid to execute ocean current acquisition operation to obtain alternative ocean current velocity vectors corresponding to each effective grid, wherein the ocean current acquisition operation is as follows:
based on a cubic spline interpolation method, respectively obtaining interpolation radial flow physical quantities corresponding to the current grid according to radial flow physical quantities of actual measurement data sampling points of two radar stations, and synthesizing the interpolation radial flow physical quantities of the current grid in the two radar stations by a double-station projection method to obtain alternative ocean current velocity vectors corresponding to the current grid;
s2: judging whether the alternative ocean current velocity vectors corresponding to the effective grids are abnormal or not according to the set rules, eliminating the abnormal alternative ocean current velocity vectors, and marking the grids corresponding to the abnormal alternative ocean current velocity vectors and the grids without actual measurement data sampling points of the radar station in the set neighborhood range as blank grids, and marking the rest grids as normal grids;
s3: respectively taking each empty grid as a current empty grid to execute two-dimensional linear interpolation operation to obtain corresponding corrected ocean current velocity vectors of each empty grid and realize full-area ocean current synthesis of the overlapped area, wherein the two-dimensional linear interpolation operation is as follows:
the method comprises the steps of respectively decomposing alternative ocean current velocity vectors of all normal grids where a current vacant grid is located into a first component and a second component which are perpendicular to each other, acquiring a first interpolation component from the first components corresponding to all the normal grids by adopting a cubic spline interpolation method, acquiring a second interpolation component from the second components corresponding to all the normal grids, and synthesizing the first interpolation component and the second interpolation component to obtain a corrected ocean current velocity vector corresponding to the current vacant grid.
Further, the method comprises the steps of respectively judging whether the alternative ocean current velocity vectors corresponding to the effective grids are abnormal according to the set rules:
forming an advance time sequence by the alternative current speed vector corresponding to the current moment and the alternative current speed vectors corresponding to the previous M moments of each effective grid, and forming a lag time sequence by the alternative current speed vector corresponding to the current moment, the alternative current speed vector corresponding to the next moment and the alternative current speed vectors corresponding to the previous M-1 moments of each effective grid, wherein M is at least 3;
based on an AR detection method, respectively acquiring ocean current velocity vector estimated values corresponding to the effective grids at the current moment according to an advance time sequence and a retard time sequence corresponding to the effective grids;
respectively obtaining first residual errors between alternative ocean current velocity vectors corresponding to the effective grids at the current moment and the ocean current velocity vector estimated values corresponding to the effective grids, and respectively obtaining standard deviations corresponding to the effective grids according to advanced time sequences corresponding to the effective grids;
judging whether the first residual error corresponding to each effective grid is greater than three times of the standard deviation corresponding to each effective grid or not respectively, and marking the effective grid with the judging result of being the suspicious grid;
the suspicious grids are used as current suspicious grids to execute the following steps to obtain a multi-element linear correlation estimated value corresponding to each suspicious grid: acquiring a lead time sequence corresponding to a neighborhood grid in a preset neighborhood range of the current suspicious grid, substituting the lead time sequence corresponding to the current suspicious grid and the lead time sequence corresponding to each neighborhood grid into a VAR model to obtain a multi-element linear correlation estimation value corresponding to the current suspicious grid;
and respectively acquiring second residual errors between the alternative ocean current velocity vectors corresponding to the suspicious grids at the current moment and the corresponding multi-element linear correlation estimated values, respectively judging whether the second residual errors corresponding to the suspicious grids are less than three times of the corresponding standard deviation, and judging whether the alternative ocean current velocity vectors corresponding to the suspicious grids are abnormal.
Further, the method comprises the steps of respectively judging whether the alternative ocean current velocity vectors corresponding to the effective grids are abnormal according to the set rules:
the effective grids are used as the current grids to execute the following steps:
obtaining an alternative ocean current velocity vector corresponding to a neighborhood grid in a neighborhood range set by a current grid;
obtaining a difference value between an alternative current velocity vector of the current grid and an alternative current velocity vector corresponding to each neighborhood grid;
inputting the difference values as characteristic items into an iForest model to obtain difference value scoring values;
and judging whether the difference value is larger than a set threshold value, and if so, judging that the alternative ocean current speed vector corresponding to the current grid is abnormal.
Further, the obtaining, by using a cubic spline interpolation method, the interpolated radial flow physical quantity corresponding to the current grid according to the radial flow physical quantities of the actual measurement data sampling points of the two radar stations respectively specifically includes:
the following steps are respectively executed by taking the two radar stations as the current radar station:
acquiring radial flow physical quantities of actually measured data sampling points on two radius boundaries of a sector detection area covered by a current radar station, and forming sampling point pairs by sampling points with the same distance with the current radar station on the two radius boundaries, wherein the number of the sampling points on each radius boundary is at least 3;
respectively acquiring arc segments formed between each sampling point pair, dividing the arc segments between each sampling point pair into more than three sub-arc segments, and representing each sub-arc segment by adopting different cubic polynomials;
acquiring intersection points between a connecting line of the current grid relative to the current radar station and each arc line segment, selecting a corresponding cubic polynomial to interpolate according to the sub-arc line segment to which each intersection point belongs, and respectively obtaining radial flow physical quantity of each intersection point;
and linearly interpolating the radial flow physical quantity of each intersection point along the radial direction to obtain the interpolation radial flow physical quantity corresponding to the current grid.
Further, the alternative ocean current speed vector comprises an alternative ocean current speed magnitude and an alternative ocean current speed direction angle, and the interpolation radial flow physical quantity of the current grid in the two radar stations is synthesized specifically by a double-station projection method:
wherein θ A For the current grid relative to the current speed direction angle, θ, of the number I radar station B For the current grid relative to the current speed direction angle of the No. II radar station, theta C For the alternative ocean current speed direction angle corresponding to the current grid, V A For the interpolated radial-flow physical quantity corresponding to the current grid obtained from the radial-flow physical quantity actually measured by the I-type radar station, V B For the interpolated radial flow physical quantity corresponding to the current grid obtained by the radial flow physical quantity actually measured by the No. II radar station, V C And the current grid is the corresponding alternative ocean current velocity vector.
Further, the overlapping region satisfies: the sea velocity direction angle of any grid in the overlapping area relative to the two radar stations is in the range of 30-150 degrees.
The beneficial effects are that:
1. the invention provides a high-frequency ground wave radar ocean current synthesis method, which is characterized in that the overlapped area covered by double stations is considered when ocean current velocity vector synthesis is carried out, and the output of a meshed flow field is completed through a series of interpolation means; for abnormal data and isolated measuring points, abnormal values are removed through an effective quality control means, and vacant data points are filled through an interpolation algorithm, so that full-area ocean current synthesis of an overlapped area is realized.
2. The invention provides a high-frequency ground wave radar ocean current synthesis method, which fully considers the characteristics of half daily tides of ocean currents, monitors the smooth change of ocean current velocity vectors on the time axis of a current grid and the correlation between the ocean current velocity vectors and adjacent grids in space, analyzes the advanced time sequence and the lagged time sequence by using an AR detection method by utilizing smooth control on the time sequence corresponding to the current grid, generates an AR model to obtain a local field estimated value, substitutes the time sequence of a neighborhood grid and the time sequence of a suspicious point into a VAR model to obtain multi-element linear correlation estimation, and finally judges that the alternative ocean current velocity vectors corresponding to the current grid are suspicious according to the relation between the alternative ocean current velocity vectors and the multi-element linear correlation estimated value; the method is fully compared with buoy, ocean model data and the like, and abnormal points can be effectively removed.
Drawings
FIG. 1 is a flow chart of a high-frequency ground wave radar ocean current synthesis method provided by the invention;
FIG. 2 is a schematic illustration of a radial flow cubic spline interpolation provided by the present invention;
FIG. 3 is a schematic diagram of vector flow synthesis provided by the present invention;
FIG. 4 is a schematic diagram of AR, VAR anomaly detection provided by the present invention;
FIG. 5 is a graph showing two effects of AR and VAR anomaly detection provided by the present invention;
fig. 6 is a schematic diagram of the isolated forest iforst space detection provided by the invention.
Detailed Description
In order to enable those skilled in the art to better understand the present application, the following description will make clear and complete descriptions of the technical solutions in the embodiments of the present application with reference to the accompanying drawings in the embodiments of the present application.
As shown in fig. 1, a high-frequency ground wave radar ocean current synthesis method comprises the following steps:
s1: and dividing grids in an overlapping area between the sector detection areas covered by the two radar stations, screening grids with actual measurement data sampling points of the radar stations in a set neighborhood range as effective grids, and respectively taking each effective grid as a current grid to execute ocean current acquisition operation to obtain alternative ocean current velocity vectors corresponding to each effective grid. It should be noted that, the angle control error is synthesized, and according to related researches, the ground wave radar vector synthesis angle is larger than the range of 30 ° and 150 °, so that the overlapping area needs to satisfy: the sea current speed direction angle of any grid in the overlapping area relative to the two radar stations is within the range of 30-150 degrees; for grid points that are not within this range, the ocean current velocity vectors are not synthesized.
Wherein, the ocean current acquisition operation is: based on a cubic spline interpolation method, respectively obtaining interpolation radial flow physical quantities corresponding to the current grid according to radial flow physical quantities of actual measurement data sampling points of two radar stations, and synthesizing the interpolation radial flow physical quantities of the current grid in the two radar stations through a double-station projection method to obtain alternative ocean current velocity vectors corresponding to the current grid.
S2: and respectively judging whether the alternative ocean current velocity vectors corresponding to the effective grids are abnormal according to the set rules, removing the abnormal alternative ocean current velocity vectors, and marking the grids corresponding to the abnormal alternative ocean current velocity vectors and the grids without actual measurement data sampling points of the radar station in the set neighborhood range as blank grids, and marking the rest grids as normal grids.
S3: respectively taking each empty grid as a current empty grid to execute two-dimensional linear interpolation operation to obtain corresponding corrected ocean current velocity vectors of each empty grid and realize full-area ocean current synthesis of the overlapped area, wherein the two-dimensional linear interpolation operation is as follows:
and decomposing the alternative ocean current velocity vectors of all the normal grids in which the current vacant grid is positioned into a first component and a second component which are perpendicular to each other, such as an east component U and a driven component V, acquiring a first interpolation component from the first components corresponding to all the normal grids by adopting a cubic spline interpolation method, acquiring a second interpolation component from the second components corresponding to all the normal grids, and synthesizing the first interpolation component and the second interpolation component to obtain a corrected ocean current velocity vector corresponding to the current vacant grid.
Two ways of determining whether the alternative ocean current velocity vector corresponding to each effective grid is abnormal are described below.
First kind: forming an advance time sequence by the alternative current speed vector corresponding to the current moment and the alternative current speed vectors corresponding to the previous M moments of each effective grid, and forming a lag time sequence by the alternative current speed vector corresponding to the current moment, the alternative current speed vector corresponding to the next moment and the alternative current speed vectors corresponding to the previous M-1 moments of each effective grid, wherein M is at least 3;
based on an AR detection method, respectively acquiring ocean current velocity vector estimated values corresponding to the effective grids at the current moment according to an advance time sequence and a retard time sequence corresponding to the effective grids;
respectively obtaining first residual errors between alternative ocean current velocity vectors corresponding to the effective grids at the current moment and the ocean current velocity vector estimated values corresponding to the effective grids, and respectively obtaining standard deviations corresponding to the effective grids according to advanced time sequences corresponding to the effective grids;
judging whether the first residual error corresponding to each effective grid is greater than three times of the standard deviation corresponding to each effective grid or not respectively, and marking the effective grid with the judging result of being the suspicious grid;
the suspicious grids are used as current suspicious grids to execute the following steps to obtain a multi-element linear correlation estimated value corresponding to each suspicious grid: acquiring a lead time sequence corresponding to a neighborhood grid in a preset neighborhood range of the current suspicious grid, such as a four neighborhood or eight neighborhood, substituting the lead time sequence corresponding to the current suspicious grid and the lead time sequence corresponding to each neighborhood grid into a VAR model to obtain a multi-element linear correlation estimation value corresponding to the current suspicious grid;
respectively acquiring second residual errors between the ocean current velocity vectors corresponding to the suspicious grids at the current moment and the corresponding multi-element linear correlation estimated values, respectively judging whether the second residual errors corresponding to the suspicious grids are less than three times of the corresponding standard deviation, if so, highly correlating the representing neighborhood grids with the change rule of the suspicious grids, and canceling the abnormal attribute; and meanwhile, judging the ocean current velocity vector corresponding to the suspicious grid with the judging result of no as abnormal.
For example, as shown in fig. 4 and 5, the data analysis diagrams of different monitoring points are respectively shown, the black thick dotted line is a 3 times standard deviation critical line, if the line is exceeded, the point is inconsistent with the integrity and is marked as a suspicious point, the gray thin dotted line represents a 3 times standard deviation critical line with the adjacent point data, if the line is exceeded, the point is inconsistent with the surrounding data point in space, and the condition is met, so that the point can be judged as an abnormal point.
Second kind: the effective grids are used as the current grids to execute the following steps:
obtaining an alternative ocean current velocity vector corresponding to a neighborhood grid in a preset neighborhood range of the current grid, such as a four neighborhood or eight neighborhood;
obtaining a difference value between an alternative current velocity vector of the current grid and an alternative current velocity vector corresponding to each neighborhood grid;
inputting the difference values as characteristic items into an iForest model to obtain difference value scoring values;
and judging whether the difference value is larger than a set threshold value, if so, indicating that the more obvious the characteristic is, the larger the difference between the current grid and the neighborhood grid is, and the greater the possibility that the current grid is an isolated point is, namely the abnormal alternative ocean current velocity vector corresponding to the current grid is.
For example, as shown in fig. 6, points with particularly extreme values and points with excessive jump and extremely uncoordinated surrounding points are marked with high anomaly scores, and the score values are at least 0.75. The iForest algorithm can accurately identify abnormal values which appear in a piece or a ring.
It should be noted that, in calculating the above field times, the parameter used is t=16,
that is, at most, only 16×64=1024 points (with repetition) are taken to train the corresponding model, and each field of data has 3000-4000 data points, which fully illustrates that the sampling process can simplify the calculation while ensuring the effect.
The method for acquiring the alternative ocean current velocity vector corresponding to the current grid is described in detail below.
The method for acquiring the interpolation radial flow physical quantity corresponding to the current grid by adopting the cubic spline interpolation method according to the radial flow physical quantity of the actual measurement data sampling points of the two radar stations comprises the following steps:
the following steps are respectively executed by taking the two radar stations as the current radar station:
acquiring radial flow physical quantities of actually measured data sampling points on two radius boundaries of a sector detection area covered by a current radar station, and forming sampling point pairs by sampling points with the same distance with the current radar station on the two radius boundaries, wherein the number of the sampling points on each radius boundary is at least 3;
respectively acquiring arc segments formed between each sampling point pair, dividing the arc segments between each sampling point pair into more than three sub-arc segments, and representing each sub-arc segment by adopting different cubic polynomials;
acquiring intersection points between a connecting line of the current grid relative to the current radar station and each arc line segment, selecting a corresponding cubic polynomial to interpolate according to the sub-arc line segment to which each intersection point belongs, and respectively obtaining radial flow physical quantity of each intersection point;
and linearly interpolating the radial flow physical quantity of each intersection point along the radial direction to obtain the interpolation radial flow physical quantity corresponding to the current grid.
For example, as shown in fig. 2, 1, performing cubic spline interpolation on radial flow, for interpolation points in a current grid, firstly determining actual measurement radial flow physical quantities of sampling points C, E, F, G, I, J on radial flow fields, namely two radial boundaries of a sector detection area covered by a current radar station, firstly calculating radial flow physical quantities of an intersection point m1 between the circumferential directions of points F and G through a cubic spline interpolation formula, and similarly obtaining radial flow physical quantities of an intersection point m2 between sampling points I, J and an intersection point m3 between E, C; since the interpolation points represented by the current grid and the intersection points m1, m2 and m3 are in the same straight line, the radial flow physical quantity of the intersection points m1, m2 and m3 is obtained through interpolation by a cubic spline interpolation formula, the interpolation points can be calculated through radial linear interpolation, and finally the radial flow physical quantity of the interpolation points is obtained.
Further, the alternative ocean current speed vector comprises an alternative ocean current speed magnitude and an alternative ocean current speed direction angle, and the interpolation radial flow physical quantity of the current grid in the two radar stations is synthesized specifically by a double-station projection method:
wherein θ A For the current grid relative to the current speed direction angle, θ, of the number I radar station B For the current grid relative to the current speed direction angle of the No. II radar station, theta C For the alternative ocean current speed direction angle corresponding to the current grid, V A For the interpolated radial-flow physical quantity corresponding to the current grid obtained from the radial-flow physical quantity actually measured by the I-type radar station, V B For the interpolated radial flow physical quantity corresponding to the current grid obtained by the radial flow physical quantity actually measured by the No. II radar station, V C And the current grid is the corresponding alternative ocean current velocity vector.
For example, as shown in FIG. 3, for two stations A, B, the radial component angles at the current grid P point are θ respectively A And theta B The flow velocity components are V respectively A 、V B Obtaining the vector angle theta at the P point by a projection method C 、V C 。
In summary, compared with other technologies, the invention has the following advantages:
1. the characteristics of the ocean current of half day tide are fully considered, and the smooth change of data on the time axis of the monitoring point and the correlation between the data and the adjacent point in space are fully considered.
2. The research results are fully compared with buoy and ocean model data, and a plurality of algorithms are optimized.
Of course, the present invention is capable of other various embodiments and its several details are capable of modification and variation in light of the present invention by one skilled in the art without departing from the spirit and scope of the invention as defined in the appended claims.
Claims (5)
1. The sea current synthesizing method for the high-frequency ground wave radar is characterized by comprising the following steps of:
s1: dividing grids in an overlapping area between the sector detection areas covered by the two radar stations, screening grids with actual measurement data sampling points of the radar stations in a set neighborhood range as effective grids, and respectively taking each effective grid as a current grid to execute ocean current acquisition operation to obtain alternative ocean current velocity vectors corresponding to each effective grid, wherein the ocean current acquisition operation is as follows:
based on a cubic spline interpolation method, respectively obtaining interpolation radial flow physical quantities corresponding to the current grid according to radial flow physical quantities of actual measurement data sampling points of two radar stations, and synthesizing the interpolation radial flow physical quantities of the current grid in the two radar stations by a double-station projection method to obtain alternative ocean current velocity vectors corresponding to the current grid;
s2: judging whether the alternative ocean current velocity vectors corresponding to the effective grids are abnormal or not according to the set rules, eliminating the abnormal alternative ocean current velocity vectors, and marking the grids corresponding to the abnormal alternative ocean current velocity vectors and the grids without actual measurement data sampling points of the radar station in the set neighborhood range as blank grids, and marking the rest grids as normal grids;
s3: respectively taking each empty grid as a current empty grid to execute two-dimensional linear interpolation operation to obtain corresponding corrected ocean current velocity vectors of each empty grid and realize full-area ocean current synthesis of the overlapped area, wherein the two-dimensional linear interpolation operation is as follows:
dividing alternative ocean current velocity vectors of all normal grids in which the current vacant grid is located into a first component and a second component which are perpendicular to each other, acquiring a first interpolation component from the first components corresponding to all the normal grids by adopting a cubic spline interpolation method, acquiring a second interpolation component from the second components corresponding to all the normal grids, and synthesizing the first interpolation component and the second interpolation component to obtain a corrected ocean current velocity vector corresponding to the current vacant grid;
the method comprises the steps of respectively judging whether the alternative ocean current velocity vectors corresponding to the effective grids are abnormal according to the set rules, wherein the specific steps are as follows:
forming an advance time sequence by the alternative current speed vector corresponding to the current moment and the alternative current speed vectors corresponding to the previous M moments of each effective grid, and forming a lag time sequence by the alternative current speed vector corresponding to the current moment, the alternative current speed vector corresponding to the next moment and the alternative current speed vectors corresponding to the previous M-1 moments of each effective grid, wherein M is at least 3;
based on an AR detection method, respectively acquiring ocean current velocity vector estimated values corresponding to the effective grids at the current moment according to an advance time sequence and a retard time sequence corresponding to the effective grids;
respectively obtaining first residual errors between alternative ocean current velocity vectors corresponding to the effective grids at the current moment and the ocean current velocity vector estimated values corresponding to the effective grids, and respectively obtaining standard deviations corresponding to the effective grids according to advanced time sequences corresponding to the effective grids;
judging whether the first residual error corresponding to each effective grid is greater than three times of the standard deviation corresponding to each effective grid or not respectively, and marking the effective grid with the judging result of being the suspicious grid;
the suspicious grids are used as current suspicious grids to execute the following steps to obtain a multi-element linear correlation estimated value corresponding to each suspicious grid: acquiring a lead time sequence corresponding to a neighborhood grid in a preset neighborhood range of the current suspicious grid, substituting the lead time sequence corresponding to the current suspicious grid and the lead time sequence corresponding to each neighborhood grid into a VAR model to obtain a multi-element linear correlation estimation value corresponding to the current suspicious grid;
and respectively acquiring second residual errors between the alternative ocean current velocity vectors corresponding to the suspicious grids at the current moment and the corresponding multi-element linear correlation estimated values, respectively judging whether the second residual errors corresponding to the suspicious grids are less than three times of the corresponding standard deviation, and judging whether the alternative ocean current velocity vectors corresponding to the suspicious grids are abnormal.
2. The method for synthesizing ocean current of high-frequency ground wave radar according to claim 1, wherein the step of respectively judging whether the alternative ocean current velocity vectors corresponding to the effective grids are abnormal according to the set rule is as follows:
the effective grids are used as the current grids to execute the following steps:
obtaining an alternative ocean current velocity vector corresponding to a neighborhood grid in a neighborhood range set by a current grid;
obtaining a difference value between an alternative current velocity vector of the current grid and an alternative current velocity vector corresponding to each neighborhood grid;
inputting the difference values as characteristic items into an iForest model to obtain difference value scoring values;
and judging whether the difference value is larger than a set threshold value, and if so, judging that the alternative ocean current speed vector corresponding to the current grid is abnormal.
3. The method for synthesizing ocean current of high-frequency ground wave radar according to claim 1, wherein the step of obtaining the interpolated radial flow physical quantity corresponding to the current grid by using a cubic spline interpolation method according to the radial flow physical quantities of the actual measurement data sampling points of the two radar stations comprises the following steps:
the following steps are respectively executed by taking the two radar stations as the current radar station:
acquiring radial flow physical quantities of actually measured data sampling points on two radius boundaries of a sector detection area covered by a current radar station, and forming sampling point pairs by sampling points with the same distance with the current radar station on the two radius boundaries, wherein the number of the sampling points on each radius boundary is at least 3;
respectively acquiring arc segments formed between each sampling point pair, dividing the arc segments between each sampling point pair into more than three sub-arc segments, and representing each sub-arc segment by adopting different cubic polynomials;
acquiring intersection points between a connecting line of the current grid relative to the current radar station and each arc line segment, selecting a corresponding cubic polynomial to interpolate according to the sub-arc line segment to which each intersection point belongs, and respectively obtaining radial flow physical quantity of each intersection point;
and linearly interpolating the radial flow physical quantity of each intersection point along the radial direction to obtain the interpolation radial flow physical quantity corresponding to the current grid.
4. The high-frequency ground wave radar ocean current synthesis method according to claim 1, wherein the alternative ocean current velocity vector comprises an alternative ocean current velocity magnitude and an alternative ocean current velocity direction angle, and the method for synthesizing the interpolation radial flow physical quantity of the current grid in the two radar stations by a double-station projection method specifically comprises the following steps:
wherein θ A For the current grid relative to the current speed direction angle, θ, of the number I radar station B For the current grid relative to the current speed direction angle of the No. II radar station, theta C For the alternative ocean current speed direction angle corresponding to the current grid, V A For the interpolated radial-flow physical quantity corresponding to the current grid obtained from the radial-flow physical quantity actually measured by the I-type radar station, V B For the interpolated radial flow physical quantity corresponding to the current grid obtained by the radial flow physical quantity actually measured by the No. II radar station, V C And the current grid is the corresponding alternative ocean current velocity vector.
5. A high frequency ground wave radar ocean current synthesizing method according to claim 1, wherein the overlapping area satisfies: the sea velocity direction angle of any grid in the overlapping area relative to the two radar stations is in the range of 30-150 degrees.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110236049.7A CN113064129B (en) | 2021-03-03 | 2021-03-03 | High-frequency ground wave radar ocean current synthesis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110236049.7A CN113064129B (en) | 2021-03-03 | 2021-03-03 | High-frequency ground wave radar ocean current synthesis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113064129A CN113064129A (en) | 2021-07-02 |
| CN113064129B true CN113064129B (en) | 2023-06-16 |
Family
ID=76559569
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110236049.7A Active CN113064129B (en) | 2021-03-03 | 2021-03-03 | High-frequency ground wave radar ocean current synthesis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113064129B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114812710B (en) * | 2022-06-28 | 2022-10-04 | 北京海兰信数据科技股份有限公司 | Vector flow synthesis method and system for radar radial flow |
| CN114997250A (en) * | 2022-08-01 | 2022-09-02 | 北京海兰信数据科技股份有限公司 | Radar radial flow post-processing method and system |
| CN115320669A (en) * | 2022-08-31 | 2022-11-11 | 南京慧尔视智能科技有限公司 | Method, device, equipment and medium for detecting railway coming car based on radar map |
| CN115525637B (en) * | 2022-10-13 | 2024-04-26 | 华润电力技术研究院有限公司 | Quality control and processing method, system and equipment for wind vector observation data |
| CN115877377B (en) * | 2022-12-12 | 2023-10-03 | 中山大学 | Radar networking vector flow field synthesis method, system, equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012207919A (en) * | 2011-03-29 | 2012-10-25 | Toyota Central R&D Labs Inc | Abnormal value determination device, positioning device, and program |
| CN105445730A (en) * | 2015-11-27 | 2016-03-30 | 南京信息工程大学 | Ocean current field inversion satellite-borne SAR system based on angle diversity, and method thereof |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7690250B2 (en) * | 2007-08-17 | 2010-04-06 | Hickey Kenneth J | Method for measuring surface currents using a long-range single station high frequency ground wave radar system |
| CN102798384B (en) * | 2012-07-03 | 2014-09-17 | 天津大学 | Ocean remote sensing image water color and water temperature monitoring method based on compression sampling |
| JP5354705B1 (en) * | 2013-05-22 | 2013-11-27 | 哲也 田川 | Radar signal processing method and apparatus |
| CN106291470B (en) * | 2016-07-28 | 2018-08-24 | 中国船舶重工集团公司第七一九研究所 | The disturbance restraining method of feature when a kind of ocean current result sky based on high-frequency ground wave radar |
| CN106019257B (en) * | 2016-08-05 | 2018-07-31 | 中国船舶重工集团公司第七一九研究所 | The interpolation method of feature when based on high-frequency ground wave radar Current Observations result sky |
| CN107678003B (en) * | 2017-09-15 | 2021-04-06 | 国家海洋局第一海洋研究所 | Target detection method under ground wave radar sea clutter background |
| CN107730582B (en) * | 2017-09-20 | 2021-08-24 | 广西师范大学 | Ocean wave three-dimensional display method based on ocean remote sensing data |
| CN110346797B (en) * | 2019-07-15 | 2021-05-25 | 中国人民解放军32181部队 | Two-dimensional inverse synthetic aperture radar imaging method and system |
-
2021
- 2021-03-03 CN CN202110236049.7A patent/CN113064129B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012207919A (en) * | 2011-03-29 | 2012-10-25 | Toyota Central R&D Labs Inc | Abnormal value determination device, positioning device, and program |
| CN105445730A (en) * | 2015-11-27 | 2016-03-30 | 南京信息工程大学 | Ocean current field inversion satellite-borne SAR system based on angle diversity, and method thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113064129A (en) | 2021-07-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113064129B (en) | High-frequency ground wave radar ocean current synthesis method | |
| EP3262440B1 (en) | Iterative ray-tracing for autoscaling of oblique ionograms | |
| CN112986946B (en) | Method for inverting undirected sea wave spectrum by using multi-frequency high-frequency radar sea echo | |
| CN105388467B (en) | A kind of method for correcting Doppler radar echo attenutation | |
| CN109100723A (en) | Upper-level winds inversion method based on Doppler radar data | |
| JP4825183B2 (en) | Surface layer flow velocity estimation method, apparatus and program | |
| CN104330787A (en) | Underwater motion array multi-target detection and position estimation integrated method | |
| JP5238531B2 (en) | Radar apparatus, marine radar observation apparatus, and Doppler frequency data calculation method | |
| US20150160374A1 (en) | Weather information signal processing module | |
| CN107871323A (en) | A kind of sea ice motion information acquisition method based on all solid state radar | |
| Gurgel et al. | HF radar observations in the German Bight: Measurements and quality control | |
| Heron et al. | Temporal and spatial resolution of HF ocean radars | |
| CN111624599A (en) | Sea wave effective wave height calculation method for sea-going radar inversion | |
| Wyatt et al. | HF radar measurements of ocean wave parameters during NURWEC | |
| CN112731382B (en) | Ratio evaluation method and system for high-frequency ground wave radar observation wind wave flow field | |
| Tian et al. | Quality control of compact high-frequency radar-retrieved wave data | |
| CN115032601A (en) | Marine radar target detection algorithm for inhibiting sea clutter in image sequence based on space-time combined filtering technology | |
| Wyatt | The IMOS ocean radar facility, ACORN | |
| Shen et al. | Spatial distribution characteristics of surface tidal currents in the southwest of Taiwan Strait | |
| Pomales-Velázquez et al. | Characterization of mesoscale eddies and detection of submesoscale eddies derived from satellite imagery and HF radar off the coast of southwestern Puerto Rico | |
| Wyatt | Improving the quality control and accuracy of HF radar currents | |
| CN111751821A (en) | MWNLCS imaging method suitable for high-resolution satellite-borne scene matching SAR | |
| CN110767322A (en) | Ocean floating platform hot spot stress calculation method based on response surface model | |
| CN116805030B (en) | Ocean wave height numerical forecast correction method | |
| CN113625272B (en) | Spatial spectrum coherent fusion imaging method for distributed radar |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |