CN113064129A - High-frequency ground wave radar ocean current synthesis method - Google Patents
High-frequency ground wave radar ocean current synthesis method Download PDFInfo
- Publication number
- CN113064129A CN113064129A CN202110236049.7A CN202110236049A CN113064129A CN 113064129 A CN113064129 A CN 113064129A CN 202110236049 A CN202110236049 A CN 202110236049A CN 113064129 A CN113064129 A CN 113064129A
- Authority
- CN
- China
- Prior art keywords
- grid
- current
- ocean current
- grids
- current velocity
- 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.)
- Granted
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 needs to consider the overlapping area covered by two stations when carrying out ocean current velocity vector synthesis and finishes the output of a gridding flow field by a series of interpolation means; for abnormal data and isolated measuring points, the invention eliminates abnormal values by an effective quality control means, fills in vacant data points by an interpolation algorithm and realizes the synthesis of the full-area ocean current of an overlapped area.
Description
Technical Field
The invention belongs to the technical field of high-frequency ground wave radars, and particularly relates to a high-frequency ground wave radar ocean current synthesis method.
Background
The high-frequency ground wave radar realizes over-the-horizon monitoring on marine dynamics parameters of marine environmental states such as wind, wave, flow and the like by utilizing the characteristic that energy attenuation is small when vertical high-frequency electromagnetic waves are transmitted along the surface of high-conductivity seawater, and has the characteristics of long observation distance, large coverage area, more inversion elements, all-weather work and the like, wherein the monitoring of ocean currents realizes conventional business operation in many areas. At present, the actual direction and the actual size of the ocean current are calculated by erecting two far-end radar stations and through the projection amount in two directions. Because the single station is monitoring data under a polar coordinate system, the overlapping area covered by the two stations needs to be considered during synthesis, and the output of the gridding flow field is completed through a series of interpolation means.
When the high-frequency ground wave radar inverts ocean current data, some abnormal data and isolated measuring points can be generated under the conditions of echo distance, environmental interference and self equipment, so that an effective quality control means is needed to eliminate abnormal values so as to fill in 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 synthesis of the ocean current in the whole area.
A high-frequency ground wave radar ocean current synthesis method comprises the following steps:
s1: carrying out grid division on an overlapped area between fan-shaped detection areas covered by two radar stations, screening out grids with actual measurement data sampling points of the radar stations in a set neighborhood range as effective grids, respectively taking each effective grid as a current grid to execute ocean current acquisition operation, and obtaining 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 actually measured 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 two-station projection method to obtain alternative ocean current velocity vectors corresponding to the current grid;
s2: respectively judging whether the alternative ocean current velocity vectors corresponding to the effective grids are abnormal or not according to a set rule, eliminating the abnormal alternative ocean current velocity vectors, recording 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 a set neighborhood range as vacant grids, and recording the rest grids as normal grids;
s3: respectively taking each vacant grid as the current vacant grid to execute two-dimensional linear interpolation operation to obtain the corresponding corrected ocean current velocity vector of each vacant grid, and realizing the synthesis of the whole-region ocean current of the overlapped region, wherein the two-dimensional linear interpolation operation is as follows:
decomposing the alternative ocean current velocity vectors of all normal grids on the row of the current vacant grid into a first component and a second component which are perpendicular to each other, respectively, obtaining a first interpolation component from the first component corresponding to all the normal grids by adopting a cubic spline interpolation method, obtaining a second interpolation component from the second component 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 specifically determining whether the candidate ocean current velocity vector corresponding to each effective grid is abnormal according to the set rule is as follows:
the alternative ocean current velocity vectors corresponding to the effective grids at the current moment and the alternative ocean current velocity vectors corresponding to the previous M moments form a leading time sequence, and meanwhile, the alternative ocean current velocity vectors corresponding to the effective grids at the current moment, the alternative ocean current velocity vectors corresponding to the next moment and the alternative ocean current velocity vectors corresponding to the previous M-1 moments form a lagging time sequence, wherein M is at least 3;
based on an AR detection method, respectively acquiring ocean current velocity vector estimation values corresponding to the effective grids at the current moment according to the lead time sequence and the lag time sequence corresponding to the effective grids;
respectively acquiring first residual errors between alternative ocean current velocity vectors corresponding to the effective grids at the current moment and ocean current velocity vector estimation values corresponding to the alternative ocean current velocity vectors and the effective grids, and respectively acquiring standard deviations corresponding to the effective grids according to the lead time sequences corresponding to the effective grids;
respectively judging whether the first residual errors corresponding to the effective grids are more than three times of the standard deviations corresponding to the effective grids, and recording the effective grids with the judged results of yes as suspicious grids;
respectively taking each suspicious grid as the current suspicious grid, and executing the following steps to obtain a multivariate linear correlation estimation value corresponding to each suspicious grid: acquiring a lead time sequence corresponding to a neighborhood grid in a set neighborhood range of the current suspicious grid, and substituting the lead time sequence corresponding to the current suspicious grid and the lead time sequences corresponding to the neighborhood grids into a VAR model to obtain a multivariate linear correlation estimation value corresponding to the current suspicious grid;
respectively obtaining second residual errors between the alternative ocean current velocity vector corresponding to each suspicious grid at the current moment and the corresponding multivariate linear correlation estimation value, and respectively judging whether the second residual error corresponding to each suspicious grid is less than three times of the corresponding standard deviation, if not, judging that the alternative ocean current velocity vector corresponding to the suspicious grid is abnormal.
Further, the specifically determining whether the candidate ocean current velocity vector corresponding to each effective grid is abnormal according to the set rule is as follows:
respectively taking each effective grid as a current grid to execute the following steps:
acquiring an alternative ocean current velocity vector corresponding to a neighborhood grid in a neighborhood range set by a current grid;
obtaining the difference value between the alternative ocean current velocity vector of the current grid and the alternative ocean current velocity vector corresponding to each neighborhood grid;
inputting each difference value as a characteristic item into an iForest model to obtain a difference value score;
and judging whether the difference value scoring value is larger than a set threshold value, if so, judging that the alternative ocean current velocity vector corresponding to the current grid is abnormal.
Further, the obtaining of the interpolated radial flow physical quantity corresponding to the current grid according to the radial flow physical quantities of the measured data sampling points of the two radar stations by using the cubic spline interpolation method specifically includes:
respectively taking two radar stations as current radar stations to execute the following steps:
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 points with the same distance with the current radar station on the two radius boundaries into sampling point pairs, wherein the number of the sampling points on each radius boundary is at least 3;
respectively obtaining arc line segments formed between each sampling point pair, respectively dividing the arc line segments between each sampling point pair into more than three sub-arc line segments, and representing each sub-arc line 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, and selecting a corresponding cubic polynomial to perform interpolation according to the sub-arc line segments to which the intersection points belong to respectively obtain radial flow physical quantities of the intersection points;
and performing linear interpolation on the radial flow physical quantity of each intersection point along the radial direction to obtain an interpolated radial flow physical quantity corresponding to the current grid.
Further, the candidate ocean current velocity vector includes a candidate ocean current velocity magnitude and a candidate ocean current velocity direction angle, and synthesizing the interpolation radial flow physical quantities of the current grid in the two radar stations by the two-station projection method specifically includes:
wherein, thetaAFor the current grid direction angle of current sea velocity, theta, relative to radar station # IBIs the current grid direction angle of current sea current velocity relative to radar station number II, thetaCFor the alternative current velocity direction angle, V, corresponding to the current gridAThe interpolation radial flow physical quantity, V, corresponding to the current grid is obtained for the radial flow physical quantity actually measured by the No. I radar stationBThe interpolation radial flow physical quantity, V, corresponding to the current grid is obtained for the radial flow physical quantity actually measured by the No. II radar stationCAnd the candidate ocean current velocity vector corresponding to the current grid is obtained.
Further, the overlap region satisfies: and the direction angle of any one grid in the overlapping area relative to the ocean current velocity of the two radar stations is in the range of 30-150 degrees.
Has the advantages that:
1. the invention provides a high-frequency ground wave radar ocean current synthesis method, which needs to consider the overlapping area covered by two stations when carrying out ocean current velocity vector synthesis and finishes the output of a gridding flow field by a series of interpolation means; for abnormal data and isolated measuring points, the invention eliminates abnormal values by an effective quality control means, fills in vacant data points by an interpolation algorithm and realizes the synthesis of the full-area ocean current of an overlapped area.
2. The invention provides a high-frequency ground wave radar ocean current synthesis method, which fully considers the characteristics of the half-day tide of ocean current, monitors the smooth change of an ocean current velocity vector on the time axis of a current grid and the correlation between the current grid and an adjacent grid in space, utilizes the smooth control on a time sequence corresponding to the current grid, analyzes an advance time sequence and a lag time sequence by using an AR (augmented reality) detection method, generates an AR (augmented reality) model to obtain an estimated value of the field, substitutes the time sequence of a neighborhood grid and the time sequence of a suspicious point into a VAR (value-variance-correlation) model to obtain a multivariate linear correlation estimation, and finally judges that an alternative ocean current velocity vector corresponding to the current grid is suspicious according to the relationship between the alternative ocean current velocity vector and the multivariate; the method is fully compared with buoy, ocean model data and the like, and abnormal points can be effectively eliminated.
Drawings
FIG. 1 is a flow chart of a method for synthesizing ocean currents of a high-frequency ground wave radar according to the present invention;
FIG. 2 is a schematic diagram of radial flow cubic spline interpolation provided by the present invention;
FIG. 3 is a schematic diagram of the vector flow synthesis provided by the present invention;
FIG. 4 is a schematic diagram of AR and VAR anomaly detection provided by the present invention;
FIG. 5 is a diagram illustrating the two effects of AR and VAR anomaly detection provided by the present invention;
fig. 6 is a schematic diagram of isolated forest IForest space detection provided by the present invention.
Detailed Description
In order to make the technical solutions better understood by those skilled in the art, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application.
As shown in fig. 1, a method for synthesizing high-frequency ground wave radar ocean current comprises the following steps:
s1: and carrying out grid division on an overlapped area between the fan-shaped detection areas covered by the two radar stations, screening out 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 synthetic angle control error shows, according to related research, that the synthetic angle of the ground wave radar vector is out of the range of 30 ° and 150 ° with a large error, and therefore, the overlapping region needs to satisfy: the direction angle of the ocean current speed of any grid relative to the two radar stations in the overlapping area is in the range of 30-150 degrees; for grid points not within this range, the ocean current velocity vectors are not synthesized.
Wherein the ocean current obtaining operation is: based on a cubic spline interpolation method, obtaining the interpolation radial flow physical quantity corresponding to the current grid according to the radial flow physical quantity of the actually measured data sampling points of the two radar stations, and synthesizing the interpolation radial flow physical quantity of the current grid in the two radar stations through a two-station projection method to obtain the alternative ocean current velocity vector corresponding to the current grid.
S2: and respectively judging whether the alternative ocean current velocity vectors corresponding to the effective grids are abnormal or not according to a set rule, eliminating the abnormal alternative ocean current velocity vectors, recording the grids corresponding to the abnormal alternative ocean current velocity vectors and the grids without the actual measurement data sampling points of the radar station in a set neighborhood range as vacant grids, and recording the rest grids as normal grids.
S3: respectively taking each vacant grid as the current vacant grid to execute two-dimensional linear interpolation operation to obtain the corresponding corrected ocean current velocity vector of each vacant grid, and realizing the synthesis of the whole-region ocean current of the overlapped region, wherein the two-dimensional linear interpolation operation is as follows:
decomposing the alternative ocean current velocity vectors of all normal grids on 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 quilt component V, obtaining a first interpolation component from the first component corresponding to all the normal grids by adopting a cubic spline interpolation method, obtaining a second interpolation component from the second component 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 candidate ocean current velocity vector corresponding to each effective grid is abnormal are described below.
The first method comprises the following steps: the alternative ocean current velocity vectors corresponding to the effective grids at the current moment and the alternative ocean current velocity vectors corresponding to the previous M moments form a leading time sequence, and meanwhile, the alternative ocean current velocity vectors corresponding to the effective grids at the current moment, the alternative ocean current velocity vectors corresponding to the next moment and the alternative ocean current velocity vectors corresponding to the previous M-1 moments form a lagging time sequence, wherein M is at least 3;
based on an AR detection method, respectively acquiring ocean current velocity vector estimation values corresponding to the effective grids at the current moment according to the lead time sequence and the lag time sequence corresponding to the effective grids;
respectively acquiring first residual errors between alternative ocean current velocity vectors corresponding to the effective grids at the current moment and ocean current velocity vector estimation values corresponding to the alternative ocean current velocity vectors and the effective grids, and respectively acquiring standard deviations corresponding to the effective grids according to the lead time sequences corresponding to the effective grids;
respectively judging whether the first residual errors corresponding to the effective grids are more than three times of the standard deviations corresponding to the effective grids, and recording the effective grids with the judged results of yes as suspicious grids;
respectively taking each suspicious grid as the current suspicious grid, and executing the following steps to obtain a multivariate linear correlation estimation value corresponding to each suspicious grid: acquiring a lead time sequence corresponding to a neighborhood grid of a four-neighborhood or eight-neighborhood set in a neighborhood range of a current suspicious grid, and substituting the lead time sequence corresponding to the current suspicious grid and the lead time sequences corresponding to the neighborhood grids into a VAR model to obtain a multivariate linear correlation estimation value corresponding to the current suspicious grid;
respectively acquiring second residual errors between the ocean current velocity vector corresponding to each suspicious grid at the current moment and the corresponding multivariate linear correlation estimation values, respectively judging whether the second residual errors corresponding to each suspicious grid are less than three times of the corresponding standard deviation, if so, representing that the neighborhood grid is also highly correlated with the change rule of the suspicious grid, and canceling abnormal attributes; and meanwhile, judging the ocean current velocity vector corresponding to the suspicious grid with a negative judgment result as abnormal.
For example, as shown in fig. 4 and 5, the data analysis graphs are respectively of different monitoring points, the black thick dotted line is a critical line with 3 times of standard deviation, if the black thick dotted line exceeds the line, the point is inconsistent with the integrity and is marked as a suspicious point, the gray thin dotted line represents a critical line with 3 times of standard deviation with the data of the adjacent point, if the black thin dotted line exceeds the line, the point is inconsistent with the surrounding data points in space, and the above conditions are met, so that the point can be determined to be an abnormal point.
And the second method comprises the following steps: respectively taking each effective grid as a current grid to execute the following steps:
acquiring alternative ocean current velocity vectors corresponding to neighborhood grids in a neighborhood range set by a current grid, such as a neighborhood grid of a four neighborhood or an eight neighborhood;
obtaining the difference value between the alternative ocean current velocity vector of the current grid and the alternative ocean current velocity vector corresponding to each neighborhood grid;
inputting each difference value as a characteristic item into an iForest model to obtain a difference value score;
and judging whether the difference value scoring value is greater than a set threshold value, if so, indicating that the characteristics are more obvious, the difference between the current grid and the neighborhood grid is larger, the possibility that the current grid is an isolated point is higher, and the alternative ocean current velocity vector corresponding to the current grid is abnormal.
For example, as shown in fig. 6, points with extreme values and points with too large jumps and extreme mismatching with surrounding points are labeled as high anomaly scores, and the score values are all at least 0.75 or more. For abnormal values of the sheet or the ring, the iForest algorithm can also accurately identify the abnormal values.
In the calculation of the above fields, the parameter t is 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 explains that the sampling process can simplify the calculation under the premise of ensuring the effect.
The manner of obtaining the candidate ocean current velocity vector corresponding to the current grid is described in detail below.
The method for obtaining the interpolation radial flow physical quantity corresponding to the current grid according to the radial flow physical quantities of the measured data sampling points of the two radar stations by adopting the cubic spline interpolation method specifically comprises the following steps:
respectively taking two radar stations as current radar stations to execute the following steps:
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 points with the same distance with the current radar station on the two radius boundaries into sampling point pairs, wherein the number of the sampling points on each radius boundary is at least 3;
respectively obtaining arc line segments formed between each sampling point pair, respectively dividing the arc line segments between each sampling point pair into more than three sub-arc line segments, and representing each sub-arc line 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, and selecting a corresponding cubic polynomial to perform interpolation according to the sub-arc line segments to which the intersection points belong to respectively obtain radial flow physical quantities of the intersection points;
and performing linear interpolation on the radial flow physical quantity of each intersection point along the radial direction to obtain an interpolated radial flow physical quantity corresponding to the current grid.
For example, as shown in fig. 2, 1, radial flow cubic spline interpolation, for an interpolation point in a current grid, firstly determining an actually measured radial flow physical quantity of a sampling point C, E, F, G, I, J on a radial flow field, that is, on two radial boundaries of a sector detection area covered by a current radar station, according to an adjacent distance parameter, firstly calculating a radial flow physical quantity of an intersection point m1 between F and G points in a circumferential direction by using a cubic spline interpolation formula, and similarly obtaining a radial flow physical quantity of an intersection point m2 between sampling points I, J and an intersection point m3 between E, C; because the interpolation points represented by the current grid and the intersection points m1, m2 and m3 are on the same straight line, and the radial flow physical quantities of the intersection points m1, m2 and m3 are obtained by interpolation through a cubic spline interpolation formula, the interpolation points can be calculated through radial linear interpolation, and the radial flow physical quantities of the interpolation points are obtained finally.
Further, the candidate ocean current velocity vector includes a candidate ocean current velocity magnitude and a candidate ocean current velocity direction angle, and synthesizing the interpolation radial flow physical quantities of the current grid in the two radar stations by the two-station projection method specifically includes:
wherein, thetaAFor the current grid direction angle of current sea velocity, theta, relative to radar station # IBIs the current grid direction angle of current sea current velocity relative to radar station number II, thetaCFor the alternative current velocity direction angle, V, corresponding to the current gridAThe interpolation radial flow physical quantity, V, corresponding to the current grid is obtained for the radial flow physical quantity actually measured by the No. I radar stationBThe interpolation radial flow physical quantity, V, corresponding to the current grid is obtained for the radial flow physical quantity actually measured by the No. II radar stationCAnd the candidate ocean current velocity vector corresponding to the current grid is obtained.
For example, as shown in FIG. 3, for A, B two stations, the radial component angles at the current grid P point are each θAAnd thetaBThe flow velocity components are respectively VA、VBObtaining the vector angle theta at the P point by a projection methodC、VC。
In summary, compared with other technologies, the advantages of the invention are as follows:
1. the characteristics of the half-day tide of ocean current are fully considered, and data on the time axis of the monitoring point changes smoothly and the correlation between the monitoring point and adjacent points in space is achieved.
2. The research result is compared with the data of the buoy and the ocean model sufficiently, and a plurality of algorithmic optimizations are completed.
The present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof, and it will be understood by those skilled in the art that various changes and modifications may be made herein without departing from the spirit and scope of the invention as defined in the appended claims.
Claims (6)
1. A high-frequency ground wave radar ocean current synthesis method is characterized by comprising the following steps:
s1: carrying out grid division on an overlapped area between fan-shaped detection areas covered by two radar stations, screening out grids with actual measurement data sampling points of the radar stations in a set neighborhood range as effective grids, respectively taking each effective grid as a current grid to execute ocean current acquisition operation, and obtaining 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 actually measured 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 two-station projection method to obtain alternative ocean current velocity vectors corresponding to the current grid;
s2: respectively judging whether the alternative ocean current velocity vectors corresponding to the effective grids are abnormal or not according to a set rule, eliminating the abnormal alternative ocean current velocity vectors, recording 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 a set neighborhood range as vacant grids, and recording the rest grids as normal grids;
s3: respectively taking each vacant grid as the current vacant grid to execute two-dimensional linear interpolation operation to obtain the corresponding corrected ocean current velocity vector of each vacant grid, and realizing the synthesis of the whole-region ocean current of the overlapped region, wherein the two-dimensional linear interpolation operation is as follows:
decomposing the alternative ocean current velocity vectors of all normal grids on the row of the current vacant grid into a first component and a second component which are perpendicular to each other, respectively, obtaining a first interpolation component from the first component corresponding to all the normal grids by adopting a cubic spline interpolation method, obtaining a second interpolation component from the second component 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.
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 candidate ocean current velocity vectors corresponding to the effective grids are abnormal according to the set rule is specifically as follows:
the alternative ocean current velocity vectors corresponding to the effective grids at the current moment and the alternative ocean current velocity vectors corresponding to the previous M moments form a leading time sequence, and meanwhile, the alternative ocean current velocity vectors corresponding to the effective grids at the current moment, the alternative ocean current velocity vectors corresponding to the next moment and the alternative ocean current velocity vectors corresponding to the previous M-1 moments form a lagging time sequence, wherein M is at least 3;
based on an AR detection method, respectively acquiring ocean current velocity vector estimation values corresponding to the effective grids at the current moment according to the lead time sequence and the lag time sequence corresponding to the effective grids;
respectively acquiring first residual errors between alternative ocean current velocity vectors corresponding to the effective grids at the current moment and ocean current velocity vector estimation values corresponding to the alternative ocean current velocity vectors and the effective grids, and respectively acquiring standard deviations corresponding to the effective grids according to the lead time sequences corresponding to the effective grids;
respectively judging whether the first residual errors corresponding to the effective grids are more than three times of the standard deviations corresponding to the effective grids, and recording the effective grids with the judged results of yes as suspicious grids;
respectively taking each suspicious grid as the current suspicious grid, and executing the following steps to obtain a multivariate linear correlation estimation value corresponding to each suspicious grid: acquiring a lead time sequence corresponding to a neighborhood grid in a set neighborhood range of the current suspicious grid, and substituting the lead time sequence corresponding to the current suspicious grid and the lead time sequences corresponding to the neighborhood grids into a VAR model to obtain a multivariate linear correlation estimation value corresponding to the current suspicious grid;
respectively obtaining second residual errors between the alternative ocean current velocity vector corresponding to each suspicious grid at the current moment and the corresponding multivariate linear correlation estimation value, and respectively judging whether the second residual error corresponding to each suspicious grid is less than three times of the corresponding standard deviation, if not, judging that the alternative ocean current velocity vector corresponding to the suspicious grid is abnormal.
3. The method for synthesizing ocean current of high-frequency ground wave radar according to claim 1, wherein the step of respectively judging whether the candidate ocean current velocity vectors corresponding to the effective grids are abnormal according to the set rule is specifically as follows:
respectively taking each effective grid as a current grid to execute the following steps:
acquiring an alternative ocean current velocity vector corresponding to a neighborhood grid in a neighborhood range set by a current grid;
obtaining the difference value between the alternative ocean current velocity vector of the current grid and the alternative ocean current velocity vector corresponding to each neighborhood grid;
inputting each difference value as a characteristic item into an iForest model to obtain a difference value score;
and judging whether the difference value scoring value is larger than a set threshold value, if so, judging that the alternative ocean current velocity vector corresponding to the current grid is abnormal.
4. The method for synthesizing ocean current of high-frequency ground wave radar according to claim 1, wherein the obtaining of the interpolated radial flow physical quantity corresponding to the current grid by using the cubic spline interpolation method according to the radial flow physical quantities of the measured data sampling points of the two radar stations respectively is specifically as follows:
respectively taking two radar stations as current radar stations to execute the following steps:
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 points with the same distance with the current radar station on the two radius boundaries into sampling point pairs, wherein the number of the sampling points on each radius boundary is at least 3;
respectively obtaining arc line segments formed between each sampling point pair, respectively dividing the arc line segments between each sampling point pair into more than three sub-arc line segments, and representing each sub-arc line 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, and selecting a corresponding cubic polynomial to perform interpolation according to the sub-arc line segments to which the intersection points belong to respectively obtain radial flow physical quantities of the intersection points;
and performing linear interpolation on the radial flow physical quantity of each intersection point along the radial direction to obtain an interpolated radial flow physical quantity corresponding to the current grid.
5. The method according to claim 1, wherein the candidate ocean current velocity vector comprises a candidate ocean current velocity magnitude and a candidate ocean current velocity direction angle, and the synthesizing of the interpolated radial flow physical quantities of the current grid in the two radar stations by the two-station projection method is specifically:
wherein, thetaAFor the current grid direction angle of current sea velocity, theta, relative to radar station # IBIs the current grid direction angle of current sea current velocity relative to radar station number II, thetaCFor the alternative current velocity direction angle, V, corresponding to the current gridAThe interpolation radial flow physical quantity, V, corresponding to the current grid is obtained for the radial flow physical quantity actually measured by the No. I radar stationBIs composed of No. II thunderObtaining the interpolation radial flow physical quantity V corresponding to the current grid by the radial flow physical quantity actually measured by the arrival stationCAnd the candidate ocean current velocity vector corresponding to the current grid is obtained.
6. A method for synthesizing high-frequency ground wave radar ocean current according to claim 1, wherein the overlapping area satisfies: and the direction angle of any one grid in the overlapping area relative to the ocean current velocity of 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 true CN113064129A (en) | 2021-07-02 |
| CN113064129B 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) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114812710A (en) * | 2022-06-28 | 2022-07-29 | 北京海兰信数据科技股份有限公司 | 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 |
| CN115525637A (en) * | 2022-10-13 | 2022-12-27 | 华润电力技术研究院有限公司 | Wind vector observation data quality control and processing method, system and equipment |
| CN115877377A (en) * | 2022-12-12 | 2023-03-31 | 中山大学 | Radar networking vector flow field synthesis method, system, equipment and storage medium |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090078037A1 (en) * | 2007-08-17 | 2009-03-26 | Hickey Kenneth J | Method for measuring surface currents using a long-range single station high frequency ground wave radar system |
| JP2012207919A (en) * | 2011-03-29 | 2012-10-25 | Toyota Central R&D Labs Inc | Abnormal value determination device, positioning device, and program |
| CN102798384A (en) * | 2012-07-03 | 2012-11-28 | 天津大学 | 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 |
| CN105445730A (en) * | 2015-11-27 | 2016-03-30 | 南京信息工程大学 | Ocean current field inversion satellite-borne SAR system based on angle diversity, and method thereof |
| CN106019257A (en) * | 2016-08-05 | 2016-10-12 | 中国船舶重工集团公司第七〇九研究所 | Interpolation method based on high-frequency ground-wave radar ocean current observation result space-time characteristics |
| CN106291470A (en) * | 2016-07-28 | 2017-01-04 | 中国船舶重工集团公司第七〇九研究所 | A kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky |
| CN107678003A (en) * | 2017-09-15 | 2018-02-09 | 国家海洋局第海洋研究所 | Object detection method and device under a kind of ground wave radar sea clutter background |
| CN107730582A (en) * | 2017-09-20 | 2018-02-23 | 广西师范大学 | Wave 3 D displaying method based on ocean remote sensing data |
| CN110346797A (en) * | 2019-07-15 | 2019-10-18 | 中国人民解放军32181部队 | A kind of two dimension inverse synthetic aperture radar imaging method and system |
-
2021
- 2021-03-03 CN CN202110236049.7A patent/CN113064129B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090078037A1 (en) * | 2007-08-17 | 2009-03-26 | Hickey Kenneth J | Method for measuring surface currents using a long-range single station high frequency ground wave radar system |
| JP2012207919A (en) * | 2011-03-29 | 2012-10-25 | Toyota Central R&D Labs Inc | Abnormal value determination device, positioning device, and program |
| CN102798384A (en) * | 2012-07-03 | 2012-11-28 | 天津大学 | 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 |
| CN105445730A (en) * | 2015-11-27 | 2016-03-30 | 南京信息工程大学 | Ocean current field inversion satellite-borne SAR system based on angle diversity, and method thereof |
| CN106291470A (en) * | 2016-07-28 | 2017-01-04 | 中国船舶重工集团公司第七〇九研究所 | A kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky |
| CN106019257A (en) * | 2016-08-05 | 2016-10-12 | 中国船舶重工集团公司第七〇九研究所 | Interpolation method based on high-frequency ground-wave radar ocean current observation result space-time characteristics |
| CN107678003A (en) * | 2017-09-15 | 2018-02-09 | 国家海洋局第海洋研究所 | Object detection method and device under a kind of ground wave radar sea clutter background |
| CN107730582A (en) * | 2017-09-20 | 2018-02-23 | 广西师范大学 | Wave 3 D displaying method based on ocean remote sensing data |
| CN110346797A (en) * | 2019-07-15 | 2019-10-18 | 中国人民解放军32181部队 | A kind of two dimension inverse synthetic aperture radar imaging method and system |
Non-Patent Citations (1)
| Title |
|---|
| 周浩: "高频地波雷达生成海洋表面矢量流图", 武汉大学学报, vol. 47, no. 5, pages 1 - 7 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114812710A (en) * | 2022-06-28 | 2022-07-29 | 北京海兰信数据科技股份有限公司 | Vector flow synthesis method and system for radar radial flow |
| 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 |
| CN115525637A (en) * | 2022-10-13 | 2022-12-27 | 华润电力技术研究院有限公司 | Wind vector observation data quality control and processing method, system and equipment |
| CN115525637B (en) * | 2022-10-13 | 2024-04-26 | 华润电力技术研究院有限公司 | Quality control and processing method, system and equipment for wind vector observation data |
| CN115877377A (en) * | 2022-12-12 | 2023-03-31 | 中山大学 | Radar networking vector flow field synthesis method, system, equipment and storage medium |
| CN115877377B (en) * | 2022-12-12 | 2023-10-03 | 中山大学 | Radar networking vector flow field synthesis method, system, equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113064129B (en) | 2023-06-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113064129A (en) | High-frequency ground wave radar ocean current synthesis method | |
| CN109100714B (en) | Low-slow small target tracking method based on polar coordinate system | |
| CA2977771C (en) | Iterative ray-tracing for autoscaling of oblique ionograms | |
| CN105487061B (en) | The multicharacteristic information fusion method of target data association | |
| CN109283499B (en) | Radar equation-based three-dimensional visualization method for detection range under active interference | |
| CN105388467B (en) | A kind of method for correcting Doppler radar echo attenutation | |
| JP2012017995A (en) | Threshold setting method, target detection method, threshold setting device, target detection device, threshold setting program and target detection program | |
| JP4825183B2 (en) | Surface layer flow velocity estimation method, apparatus and program | |
| CN104330787A (en) | Underwater motion array multi-target detection and position estimation integrated method | |
| JP6334730B2 (en) | Tracking processing apparatus and tracking processing method | |
| CN106872970A (en) | A kind of multiple target based on differential evolution becomes data transfer rate tracks of device and its method | |
| CN113359123A (en) | System and method for protecting offshore wind farm submarine cable based on AIS and radar information fusion | |
| CN111487597A (en) | Universal electronic reconnaissance interception simulation method based on time-space-frequency energy data | |
| CN107871323A (en) | A kind of sea ice motion information acquisition method based on all solid state radar | |
| CN102269807A (en) | Realization method of radar anti-jamming technology based on AIS (automatic identification system) information | |
| CN111983610B (en) | SAR image positioning method based on azimuth scaling and spectral analysis imaging | |
| CN111624599A (en) | Sea wave effective wave height calculation method for sea-going radar inversion | |
| Xu et al. | Survey of track association of radar and AIS | |
| Chen et al. | Research on AIS and radar information fusion method based on distributed Kalman | |
| CN112731382B (en) | Ratio evaluation method and system for high-frequency ground wave radar observation wind wave flow field | |
| CN112731400B (en) | Method and system for estimating target vector speed of marine vessel | |
| CN111190168A (en) | Posture stabilizing method of side-scan sonar | |
| Wyatt | The IMOS ocean radar facility, ACORN | |
| Zhou et al. | A quality factor of forecasting error for sounding data in MBES | |
| CN110767322A (en) | Ocean floating platform hot spot stress calculation method based on response surface model |
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 |