CN106291470A - A kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky - Google Patents

A kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky Download PDF

Info

Publication number
CN106291470A
CN106291470A CN201610604716.1A CN201610604716A CN106291470A CN 106291470 A CN106291470 A CN 106291470A CN 201610604716 A CN201610604716 A CN 201610604716A CN 106291470 A CN106291470 A CN 106291470A
Authority
CN
China
Prior art keywords
ocean current
mesh point
data
wave radar
ground wave
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
Application number
CN201610604716.1A
Other languages
Chinese (zh)
Other versions
CN106291470B (en
Inventor
沈志奔
邓海华
彭亮
董梁
王赟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
719th Research Institute of CSIC
Original Assignee
719th Research Institute of CSIC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 719th Research Institute of CSIC filed Critical 719th Research Institute of CSIC
Priority to CN201610604716.1A priority Critical patent/CN106291470B/en
Publication of CN106291470A publication Critical patent/CN106291470A/en
Application granted granted Critical
Publication of CN106291470B publication Critical patent/CN106291470B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention provides a kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky, the method can form complementation with existing interference suppression algorithm, improve the accuracy of high-frequency ground wave radar ocean current result.Comprise the steps of: being obtained ocean current field data by high-frequency ground wave radar observation, ocean current field data is the grid data being made up of multiple spatial point, utilizes the disappearance mesh point of room and time interpolation algorithm polishing ocean current field data and builds covariance matrix.Covariance matrix is carried out feature decomposition and obtains spatial signature vectors corresponding to each mesh point and time coefficient sequence.The spatial signature vectors of each mesh point is spatially carried out interference removal.Each mesh point time coefficient sequence carries out interference remove.Threshold value S according to the cumulative proportion in ANOVA set judges primary modal number l.Front l mode is reconstructed by characteristic vector and time coefficient after utilizing interference to remove, obtains the data after AF panel after reconstruct.

Description

A kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky
Technical field
The invention belongs to radio ocean remote sensing and digital processing field, be specifically related to a kind of based on high-frequency ground wave thunder The disturbance restraining method of feature when reaching ocean current result sky.
Background technology
High-frequency ground wave radar (High Frequency Surface Wave Radar, HFSWR) utilizes shortwave (3- 30MHz) the ocean surface diffraction propagation little feature of decay can with large area, round-the-clock, time the real table providing observation marine site Layer ocean current result.The surface flow field that high-frequency ground wave radar observation obtains, through substantial amounts of marine contrast experiment, demonstrates ocean current Accuracy, has promoted the businessization of high-frequency ground wave radar Current Observations to run, and during for studying local flow field sky, characteristic provides Well data support.The advantages such as high-frequency ground wave radar has cheap, and time empty resolution is high.Therefore, high-frequency ground wave radar It is the most cheap effective monitoring instrument in exclusive economic zone, there is market prospect widely.But in the process of actual motion In, owing to there is various interference (ionospheric interference, glitch etc.), cause high-frequency ground wave radar ocean current result to exist the most different In the result of normal value, have a strong impact on the use of high-frequency ground wave radar Current Observations result.
Traditional disturbance restraining method such as AR model linear prediction methods etc. also destroy accurately while suppression interference High-frequency ground wave radar ocean current result empty time feature, its computation complexity is high simultaneously and is difficult to engineer applied.
Summary of the invention
In view of this, the invention provides a kind of based on the AF panel side of feature during high-frequency ground wave radar ocean current result sky Method, the method can form complementation with existing interference suppression algorithm, improve the accuracy of high-frequency ground wave radar ocean current result.
In order to achieve the above object, when the technical scheme is that one based on high-frequency ground wave radar ocean current result sky The disturbance restraining method of feature, comprises the steps of:
Step one, by described high-frequency ground wave radar observation obtain ocean current field data, described ocean current field data is by multiple skies Between put the grid data of composition, utilize the disappearance mesh point of ocean current field data described in room and time interpolation algorithm polishing and build Covariance matrix.
Step 2, the feature decomposition that carries out covariance matrix obtain spatial signature vectors corresponding to each mesh point and time system Number Sequence.
Step 3, the spatial signature vectors of each mesh point is spatially carried out interference remove.
Step 4, each mesh point time coefficient sequence is carried out interference remove.
Step 5, judge primary modal number l according to threshold value S of cumulative proportion in ANOVA set.
Front l mode is reconstructed by characteristic vector and time coefficient after step 6, utilization interference removal, obtains after reconstruct Must be through the data after AF panel.
Further, the concrete grammar of ocean current field data disappearance result described in room and time interpolation algorithm polishing is utilized For:
Space interpolation algorithm is: for the disappearance mesh point in ocean current spatial sequence, first finds out this disappearance mesh point week Enclosing the shortage of data rate of other 8 mesh points in the range of 3 × 3, if miss rate is more than 62.5%, then this mesh point cannot be mended Together, this mesh point is not carried out interpolation processing；If shortage of data rate is below 62.5%, uses in ocean current spatial sequence and exist Gridden data to carry out inverse distance-weighting average, obtain lacking the ocean current spatial data of mesh point.
Time difference algorithm is particularly as follows: for ocean current time series, use pusher 12.4 hours to the period of missing data Data carry out interpolation；If the pusher data of 12.4 hours also vacancy, then 12.4 hour datas that push away before using, for pusher and front Push away all cannot period of polishing, utilize cubic spline interpolation to carry out polishing.
Further, step 3 particularly as follows:
For spatial signature vectors, calculate respectively around each mesh point in the range of 3 × 3 the space characteristics of mesh point to Amount, and calculate average and the standard deviation of the spatial signature vectors of all mesh points, exhausted by spatial signature vectors and the difference of average The average of the spatial signature vectors value spatial signature vectors of the value mesh point more than 2 times of standard deviation is replaced, it is thus achieved that through dry Disturb the spatial signature vectors after removal.
Further, step 4 is particularly as follows: for the time coefficient sequence of each mesh point, calculate its average and standard Difference, by this mesh point time coefficient sequence with this mesh point of numerical value more than 2 times of standard deviation of the absolute value of the difference of average Time coefficient average substitute, it is thus achieved that through interference remove after time coefficient sequence.
Beneficial effect:
1, when the present invention utilizes high-frequency ground wave radar ocean current result empty, feature carries out AF panel, does not relies on priori and knows Knowing, have good robustness, the ocean current result obtained after AF panel is the most credible；
2, the present invention is applicable not only to high-frequency ground wave radar, applies also for folded Clutter in Skywave Radars and other various is returned by ocean Ripple carries out radar and the equipment of ocean current inverting.
Accompanying drawing explanation
Accompanying drawing 1 is the algorithm flow chart of the present invention.
Detailed description of the invention
Develop simultaneously embodiment below in conjunction with the accompanying drawings, describes the present invention.
Embodiment 1,
Step one, set the ocean current field that obtains of high-frequency ground wave radar observation as U (Χm,tn)=u (Χm,tn)+jv(Χm,tn), Wherein Χm=(x, y), m=1-M, represent the coordinate of M spatial point.N=1-N, representing length of time series is N.J=(-1 )1/2, represent imaginary number.U, v are respectively eastern component and the north component of ocean current.
For the missing values of each spatially, in the range of first finding out around this missing values mesh point 3 × 3 other 8 The shortage of data rate of individual mesh point, if shortage of data rate is more than 62.5%, that is has 5 shortage of datas in 8 mesh points, Then this mesh point cannot polishing, then just this mesh point is not carried out interpolation processing.If shortage of data rate is below 62.5%, That is with the presence of at least 4 some valid data in 8 mesh points, then just carry out inverse distance-weighting by the result having data point and put down All, the ocean current of missing point grid is obtained.
For the ocean current time series of each spatial point, the pusher data of 12.4 hours are used the period of missing data Interpolation.During data also vacancy when pusher 12.4 is little, the data pushed away before employing.All cannot the number of polishing for pusher and front pushing away According to, utilize cubic spline interpolation to carry out polishing.
By ocean current field U (Χm,tn) do anomaly process, i.e. the time series of the ocean current in each spatial point is deducted this The process of seasonal effect in time series average ocean current.Ocean current field after having processed still is designated as U (Χm,tn)。
Calculate U (Xm,tn) covariance matrix R (Xi,Xm)
R ( X i , X m ) = 1 N Σ n = 1 N U ( X i , t n ) U * ( X m , t n ) - - - ( 1 )
Step 2, the characteristic root λ and characteristic vector Φ (Χ of calculating R (Xi, Xm)m,Χm) so that the two meets following formula:
Σ m = 1 M R ( X i , X m ) Φ ( X m , X m ) = λ Φ ( X i , X i ) - - - ( 2 )
Φ (Χ in formula (2)m,Χm) and λ be M × M dimension matrix, wherein Φ (Χm,Χm) every string between be all Mutually orthogonal.λ has a following form:
λ = λ 1 0 ... 0 0 λ 2 ... 0 ... ... ... ... 0 0 ... λ M - - - ( 3 )
In formula (3), characteristic root is to arrange by descending order, i.e. λ1＞ λ2＞ ... ＞ λM, each eigenvalue correspond to Φ (ΧmString characteristic vector value in), also referred to as one mode.
Calculate the time coefficient T (Χ that characteristic vector is correspondingm,tn):
T(Χm,tn)=Φ (Χm,Χm)U(Χm,tn) (4)
T (Χ in formula (4)m,tnIn), each row of data is exactly the time coefficient of corresponding each mode.
Step 3, for characteristic vector Φ (Χm,Χm), east component in the range of calculating around each mesh point 3 × 3 respectively ΦEWith north component ΦEThe standard deviation sigma of characteristic vector, then by characteristic vector and characteristic vector average (With) difference Absolute value more than 2 times of standard deviations (Or) mesh point characteristic vector value feature to Amount average replaces.Characteristic vector after processing is designated as Φ1(Χm)
Step 4, for time coefficient T (Χm,tn), the sequence of the time coefficient of each grid is designated as y (tn), y (tn) Average be designated asIts standard deviation is designated as σ, will meetThe value of time coefficient with this mesh point The average of time coefficientSubstitute.Time coefficient after processing is designated as T1(Χm,tn)
Step 5, threshold value S (generally 90%) of cumulative proportion in ANOVA of the characteristic root that sets up an office judge primary modal number Number, if the primary modal number determined is l.
Step 6, according to formula (5) utilize before l characteristic modes ocean current result is reconstructed.Reconstruction value will be obtained add The average in time series obtained in step 4, the data set after finally being reconstructed, this data set is after AF panel Result.
Ure(Χm,tn)=Φ1(Χm,1:l)-1·T1(1:l,tn) (5)。
To sum up, these are only presently preferred embodiments of the present invention, be not intended to limit protection scope of the present invention.All Within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. made, should be included in the protection of the present invention Within the scope of.

Claims (4)

1. one kind based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky, it is characterised in that comprise as follows Step:
Step one, by described high-frequency ground wave radar observation obtain ocean current field data, described ocean current field data is by multiple spatial point The grid data of composition, utilizes the disappearance mesh point of ocean current field data described in room and time interpolation algorithm polishing and builds association side Difference matrix；
Step 2, covariance matrix is carried out feature decomposition obtain spatial signature vectors corresponding to each mesh point and time coefficient sequence Row；
Step 3, the spatial signature vectors of each mesh point is spatially carried out interference remove；
Step 4, each mesh point time coefficient sequence is carried out interference remove；
Step 5, judge primary modal number l according to threshold value S of cumulative proportion in ANOVA set；
Front l mode is reconstructed by characteristic vector and time coefficient after step 6, utilization interference removal, obtains warp after reconstruct Data after AF panel.
A kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky, its It is characterised by, described ocean current field data described in the room and time interpolation algorithm polishing disappearance result of utilizing method particularly includes:
Space interpolation algorithm is: for the disappearance mesh point in ocean current spatial sequence, first find out around this disappearance mesh point 3 × The shortage of data rate of other 8 mesh points in the range of 3, if miss rate is more than 62.5%, then this mesh point cannot polishing, no This mesh point is carried out interpolation processing；If shortage of data rate is below 62.5%, use net present in ocean current spatial sequence It is average that Grid data carries out inverse distance-weighting, obtains lacking the ocean current spatial data of mesh point；
Time difference algorithm is particularly as follows: for ocean current time series, use the pusher number of 12.4 hours to the period of missing data According to carrying out interpolation；If the pusher data of 12.4 hours also vacancy, then 12.4 hour datas that push away before using, for pusher and before push away all Cannot period of polishing, utilize cubic spline interpolation to carry out polishing.
A kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky, its Be characterised by, described step 3 particularly as follows:
For spatial signature vectors, calculate around each mesh point the spatial signature vectors of mesh point in the range of 3 × 3 respectively, and Calculate average and the standard deviation of the spatial signature vectors of all mesh points, spatial signature vectors is big with the absolute value of the difference of average Average in the spatial signature vectors value spatial signature vectors of the mesh points of 2 times of standard deviation replaces, it is thus achieved that remove through interference After spatial signature vectors.
4. a kind of based on the AF panel side of feature during high-frequency ground wave radar ocean current result sky as described in claim 2 or 3 Method, it is characterised in that described step 4 is particularly as follows: for the time coefficient sequence of each mesh point, calculate its average and standard Difference, by this mesh point time coefficient sequence with this mesh point of numerical value more than 2 times of standard deviation of the absolute value of the difference of average Time coefficient average substitute, it is thus achieved that through interference remove after time coefficient sequence.
CN201610604716.1A 2016-07-28 2016-07-28 The disturbance restraining method of feature when a kind of ocean current result sky based on high-frequency ground wave radar Expired - Fee Related CN106291470B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610604716.1A CN106291470B (en) 2016-07-28 2016-07-28 The disturbance restraining method of feature when a kind of ocean current result sky based on high-frequency ground wave radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610604716.1A CN106291470B (en) 2016-07-28 2016-07-28 The disturbance restraining method of feature when a kind of ocean current result sky based on high-frequency ground wave radar

Publications (2)

Publication Number Publication Date
CN106291470A true CN106291470A (en) 2017-01-04
CN106291470B CN106291470B (en) 2018-08-24

Family

ID=57663365

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610604716.1A Expired - Fee Related CN106291470B (en) 2016-07-28 2016-07-28 The disturbance restraining method of feature when a kind of ocean current result sky based on high-frequency ground wave radar

Country Status (1)

Country Link
CN (1) CN106291470B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113064129A (en) * 2021-03-03 2021-07-02 湖北中南鹏力海洋探测系统工程有限公司 High-frequency ground wave radar ocean current synthesis method
US11058341B2 (en) * 2015-06-02 2021-07-13 Riken Reconstructing device, reconstructing method, program, and information recording medium
CN117763259A (en) * 2023-05-26 2024-03-26 国家海洋局北海预报中心((国家海洋局青岛海洋预报台)(国家海洋局青岛海洋环境监测中心站)) High-frequency ground wave radar flow field observation data reconstruction algorithm

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2478940C (en) * 2002-03-13 2008-07-08 Raytheon Canada Limited A noise suppression system and method for phased-array based systems
CN102538768A (en) * 2012-01-10 2012-07-04 武汉大学 Method for measuring water depth of shallow sea based on double-frequency high-frequency ground wave radar
CN102540153A (en) * 2011-12-26 2012-07-04 武汉大学 Array amplitude and phase error correcting method based on interstation direct wave interference
CN104133198A (en) * 2014-08-13 2014-11-05 武汉大学 Ionized-layer interference suppression method used in high frequency ground wave radar

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2478940C (en) * 2002-03-13 2008-07-08 Raytheon Canada Limited A noise suppression system and method for phased-array based systems
CN102540153A (en) * 2011-12-26 2012-07-04 武汉大学 Array amplitude and phase error correcting method based on interstation direct wave interference
CN102538768A (en) * 2012-01-10 2012-07-04 武汉大学 Method for measuring water depth of shallow sea based on double-frequency high-frequency ground wave radar
CN104133198A (en) * 2014-08-13 2014-11-05 武汉大学 Ionized-layer interference suppression method used in high frequency ground wave radar

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LUN LI ET AL: "Ocean gravity wave phase velocity detection by HFSWR", 《THE SOUTHWEST OF TAIWAN STRAIT》 *
SHEN ZHIBEN ET AL: "Spatial Distribution Characteristics of Surface Tidal Currents in the Southwest of Taiwan Strait", 《OCEANIC AND COASTAL SEA RESEARCH》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11058341B2 (en) * 2015-06-02 2021-07-13 Riken Reconstructing device, reconstructing method, program, and information recording medium
CN113064129A (en) * 2021-03-03 2021-07-02 湖北中南鹏力海洋探测系统工程有限公司 High-frequency ground wave radar ocean current synthesis method
CN117763259A (en) * 2023-05-26 2024-03-26 国家海洋局北海预报中心((国家海洋局青岛海洋预报台)(国家海洋局青岛海洋环境监测中心站)) High-frequency ground wave radar flow field observation data reconstruction algorithm

Also Published As

Publication number Publication date
CN106291470B (en) 2018-08-24

Similar Documents

Publication Publication Date Title
US10439594B2 (en) Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation
CN107728211B (en) Seismic signal recovery algorithm based on tensor nuclear norm regularization
CN106019257B (en) The interpolation method of feature when based on high-frequency ground wave radar Current Observations result sky
CN105913393A (en) Self-adaptive wavelet threshold image de-noising algorithm and device
CN106291470A (en) A kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky
CN104777442A (en) MRS (magnetic resonance sounding) FID (frequency identity) signal noise inhibition method
Martins et al. Percolation on two-and three-dimensional lattices
Yang et al. Seislet-based morphological component analysis using scale-dependent exponential shrinkage
CN106156451A (en) A kind of based on the Mutual coupling technology improving quantum particle swarm
CN104977583A (en) Method for X-band radar wave retrieval based on empirical orthogonal decomposition
CN104022756B (en) A kind of particle filter method of improvement based on GPU architecture
CN107271979B (en) The double quartile the point estimation methods of Pareto distribution with wide scope parametric joint
CN104318046A (en) System and method for incrementally converting high dimensional data into low dimensional data
CN104567937A (en) High-precision X-ray pulsar signal simulation method
CN114814897A (en) Combined filtering method for weakening GNSS multipath error
Ge et al. A new MCMC algorithm for blind Bernoulli-Gaussian deconvolution
CN115373006A (en) Pseudo-range multi-path error extraction method for smart phone
Gigolashvili et al. Statistical study of the north–south asymmetry of the solar differential rotation based on various solar structures during 1966–1985
CN110687594B (en) Instantaneous phase unwrapping method, full waveform inversion method and computer equipment
Kong Analysis of tunnel monitoring results based on modulus maxima method of wavelet transform
Qiu et al. Iteration empirical mode decomposition method for filling the missing data of GNSS position time series
Wu et al. Study of GPS data de-noising method based on Wavelet and Kalman filtering
Gu et al. Polychromatic sparse image reconstruction and mass attenuation spectrum estimation via B-spline basis function expansion
CN110687605A (en) Improved K-SVD algorithm-based algorithm analysis application in seismic signal processing
Xian-Mei An effective approach for improving the real-time prediction of summer rainfall over China

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180824

Termination date: 20190728

CF01 Termination of patent right due to non-payment of annual fee