CN114812710B - Vector flow synthesis method and system for radar radial flow - Google Patents
Vector flow synthesis method and system for radar radial flow Download PDFInfo
- Publication number
- CN114812710B CN114812710B CN202210737877.3A CN202210737877A CN114812710B CN 114812710 B CN114812710 B CN 114812710B CN 202210737877 A CN202210737877 A CN 202210737877A CN 114812710 B CN114812710 B CN 114812710B
- Authority
- CN
- China
- Prior art keywords
- flow
- radar
- grid
- grid unit
- vector
- 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
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/66—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
- G01F1/663—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters by measuring Doppler frequency shift
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
-
- 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/88—Radar or analogous systems specially adapted for specific applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Aviation & Aerospace Engineering (AREA)
- Fluid Mechanics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The embodiment of the invention discloses a vector flow synthesis method of radar radial flow, which comprises the following steps: determining a plurality of target radars that collectively cover a target detection sea area; establishing a rectangular grid for a radial velocity map of each target radar and initializing sea state information; determining the number of radar combinations; synthesizing the radial velocity map of each radar combination, and performing first post-processing on the synthesized vector flow to obtain a vector flow map; for the quality factor corresponding to the synthetic vector flow mark of each grid unit in each vector flow graph, selecting the synthetic vector flow with the highest quality factor as the ocean current result of the grid unit; and performing secondary post-processing on the ocean current results of each grid unit to obtain final ocean current results, and obtaining a vector flow field in the target detection sea area based on the final ocean current results of each grid unit. The embodiment of the invention also discloses a system for synthesizing the vector flow of the radar radial flow. The invention improves the accuracy of radar ocean current detection and improves the stability of radar net ocean current detection.
Description
Technical Field
The invention relates to the technical field of radar, in particular to a method and a system for synthesizing a vector flow of a radar radial flow.
Background
In the inversion technology of detecting radial flow by high-frequency ground wave radar, two algorithms are generally adopted for the synthesis of vector ocean current: 1) Solving the vector flow velocity and the flow direction through a single radar; 2) The radial streams passing through the two radars are vector-synthesized. For the second algorithm, the quality of the radial flow participating in the vector synthesis is not marked, so that the radial flow with poor quality is easily introduced, and the error is increased; the synthetic vector flow is not subjected to vector marking, so that the low-quality synthetic vector flow cannot be deleted, the obtained ocean current data has low reliability, and the radar ocean current detection precision is low; because only the vector synthesis is carried out on the radial flow of the two radars, the radial flow of the multiple radars cannot be fused when the radars are networked, if one of the radars cannot work normally, the vector flow cannot be synthesized, and the stability of radar network ocean current detection is low.
Disclosure of Invention
In order to solve the above problems, an object of the present invention is to provide a method and a system for synthesizing a vector flow of a radar radial flow, so that ocean flow data with high reliability can be obtained during ocean current detection, the accuracy of ocean current detection of a radar is improved, and the stability of ocean current detection of a radar network is improved.
The embodiment of the invention provides a vector flow synthesis method of radar radial flow, which comprises the following steps:
determining a plurality of target radars which jointly cover a target detection sea area according to the coverage sea area of each radar, wherein the coverage sea area of each radar is obtained according to radial flow data corresponding to the current moment, and the number of the target radars is more than or equal to 3;
acquiring a radial velocity map for each target radar, establishing a rectangular grid in a target detection sea area for the radial velocity map, and initializing sea state information on each grid unit, wherein the rectangular grid comprises a plurality of grid units, and the sea state information comprises velocity and flow;
determining the number of radar combinations
And forming a set of radar combinations, wherein,mrepresenting the number of the plurality of target radars,nrepresenting the number of target radars contained in each radar combination;
determining whether each grid unit meets a synthesis condition or not for a radial velocity map of a target radar contained in each radar combination in the set, synthesizing vector flows of the grid units meeting the synthesis condition, and performing first post-processing on the synthesized vector flow of each grid unit to obtain a vector flow diagram;
marking a corresponding quality factor for the synthetic vector flow of each grid unit in each vector flow graph, and selecting the synthetic vector flow with the highest quality factor in a plurality of vector flow graphs as an ocean current result of the grid unit for the same grid unit, wherein the quality factor is obtained by weighting calculation according to the synthetic geometric precision factor of the grid unit;
and performing secondary post-processing on the ocean current results of the grid units to obtain a final ocean current result of each grid unit, and obtaining a vector flow field in the target detection sea area based on the final ocean current result of each grid unit.
As a further improvement of the present invention, the synthesis conditions include:
the grid units at the same position in each radial flow velocity diagram of the radar combination do not exceed the outer edge of the flow velocity area;
none of the grid cells at the same position in each radial velocity map of the radar assembly are located in the baseline instability region.
As a further improvement of the present invention, for a grid cell in the baseline unstable region, if the grid cell does not have a synthetic vector stream, interpolation and completion are performed on the grid cell to obtain an ocean current result of the grid cell.
As a further improvement of the invention, the first post-treatment comprises:
and performing one or more of outlier elimination, flow field smoothing filtering, minimum variance flow velocity taking and flow field interpolation on the synthetic vector flow of each grid unit.
As a further improvement of the invention, the second post-treatment comprises:
and (3) performing one or more of outlier rejection, flow field smoothing filtering, minimum variance flow velocity taking and flow field interpolation on the ocean current result of each grid unit.
The embodiment of the invention also provides a system for synthesizing the vector flow of the radar radial flow, which comprises the following components:
the target radar determining module is used for determining a plurality of target radars which jointly cover a target detection sea area according to the coverage sea area of each radar, wherein the coverage sea area of each radar is obtained according to radial flow data corresponding to the current moment, and the number of the target radars is more than or equal to 3;
the system comprises a radial velocity map initialization module, a target detection sea area detection module and a sea state information initialization module, wherein the radial velocity map initialization module is used for acquiring a radial velocity map for each target radar, establishing a rectangular grid for the radial velocity map in a target detection sea area, and initializing sea state information on each grid unit, the rectangular grid comprises a plurality of grid units, and the sea state information comprises velocity and flow;
a radar combination determination module for determining the number of radar combinations
And forming a set of radar combinations, wherein m represents the number of the target radars, and n represents the number of the target radars contained in each radar combination;
the vector flow synthesis module is used for determining whether each grid unit meets synthesis conditions or not for the radial velocity map of the target radar contained in each radar combination in the set, synthesizing the vector flows of the grid units meeting the synthesis conditions, and performing first post-processing on the synthesized vector flow of each grid unit to obtain a vector flow map;
the quality mark evaluation module is used for marking a corresponding quality factor for the synthetic vector flow of each grid unit in each vector flow graph, and selecting the synthetic vector flow with the highest quality factor in a plurality of vector flow graphs as an ocean current result of the grid unit for the same grid unit, wherein the quality factor is obtained by weighting calculation according to the synthetic geometric precision factor of the grid unit;
and the result output module is used for carrying out secondary post-processing on the ocean current results of the grid units to obtain the final ocean current result of each grid unit, and obtaining the vector flow field in the target detection sea area based on the final ocean current result of each grid unit.
As a further improvement of the invention, the synthesis conditions comprise:
grid cells at the same position in each radial flow velocity map of the radar combination do not exceed the outer edge of the flow velocity area;
none of the grid cells at the same position in each radial velocity map of the radar assembly are located in the baseline instability region.
As a further improvement of the present invention, for a grid cell in the baseline unstable region, if the grid cell does not have a synthetic vector stream, interpolation and completion are performed on the grid cell to obtain an ocean current result of the grid cell.
As a further improvement of the invention, the first post-treatment comprises:
and performing one or more of outlier elimination, flow field smoothing filtering, minimum variance flow velocity taking and flow field interpolation on the synthetic vector flow of each grid unit.
As a further improvement of the invention, the second post-treatment comprises:
and (3) performing one or more of wild value elimination, flow field smoothing filtering, minimum variance flow velocity taking and flow field interpolation on the ocean current result of each grid unit.
The invention has the beneficial effects that:
the multi-graph fusion synthesis is adopted for the radial velocity map of the radar network, so that the radial flow of all target radars which jointly cover a target detection sea area indirectly participates in the synthesis of a final vector flow field, and the relation between the final vector flow field precision and key factors influencing the detection precision, such as sea conditions, velocity spatial distribution, noise environment, radar parameters and the like, is indirectly established, so that sea flow data with higher reliability can be obtained during sea flow detection, the sea flow detection precision of the radar is improved, and the stability of the radar network sea flow detection is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below. It is obvious that the drawings in the following description are only some embodiments of the invention, and that for a person skilled in the art, other drawings can be derived from them without inventive effort.
Fig. 1 is a schematic flowchart of a method for synthesizing a vector flow of a radar radial flow according to an exemplary embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that, if directional indications (such as up, down, left, right, front, back, 8230; etc.) are involved in the embodiment of the present invention, the directional indications are only used for explaining the relative positional relationship between the components, the motion situation, etc. in a specific posture (as shown in the figure), and if the specific posture is changed, the directional indications are correspondingly changed.
In addition, in the description of the present invention, the terms used are for illustrative purposes only and are not intended to limit the scope of the present invention. The terms "comprises" and/or "comprising" are used to specify the presence of stated elements, steps, operations, and/or components, but do not preclude the presence or addition of one or more other elements, steps, operations, and/or components. The terms "first," "second," and the like may be used to describe various elements, not necessarily order, and not necessarily limit the elements. In addition, in the description of the present invention, "a plurality" means two or more unless otherwise specified. These terms are only used to distinguish one element from another. These and/or other aspects will become apparent to those of ordinary skill in the art in view of the following drawings, and the description of the embodiments of the present invention will be more readily understood by those of ordinary skill in the art. The figures depict described embodiments of the present invention for purposes of illustration only. One skilled in the art will readily recognize from the following description that alternative embodiments of the structures and methods illustrated in the present application may be employed without departing from the principles described in the present application.
As shown in fig. 1, the method for synthesizing a vector flow of a radar radial flow according to an embodiment of the present invention includes:
determining a plurality of target radars which jointly cover a target detection sea area according to the coverage sea area of each radar, wherein the coverage sea area of each radar is obtained according to radial flow data corresponding to the current moment, and the number of the target radars is more than or equal to 3;
acquiring a radial velocity map for each target radar, establishing a rectangular grid in a target detection sea area for the radial velocity map, and initializing sea state information on each grid unit, wherein the rectangular grid comprises a plurality of grid units, and the sea state information comprises velocity and flow;
determining the number of radar combinations
And forming a set of radar combinations, wherein,mrepresenting the number of the plurality of target radars,nrepresenting the number of target radars contained in each radar combination;
determining whether each grid unit meets a synthesis condition or not for a radial velocity map of a target radar contained in each radar combination in the set, synthesizing vector flows of the grid units meeting the synthesis condition, and performing first post-processing on the synthesized vector flow of each grid unit to obtain a vector flow graph;
for a quality factor corresponding to a synthetic vector flow mark of each grid unit in each vector flow graph, selecting a synthetic vector flow with the highest quality factor in a plurality of vector flow graphs as an ocean current result of the grid unit for the same grid unit, wherein the quality factor is obtained by weighting calculation according to a synthetic geometric precision factor of the grid unit;
and performing secondary post-processing on the ocean current results of the grid units to obtain a final ocean current result of each grid unit, and obtaining a vector flow field in the target detection sea area based on the final ocean current result of each grid unit.
High-frequency ground wave radar is widely applied to ocean surface state monitoring, and ocean state parameters mainly comprise three aspects of wind, wave and flow field. The key point of sea state parameter monitoring is sea state inversion algorithm, wherein the sea state inversion algorithm is the most mature, the flow field parameters are mainly inversion of the flow velocity and the flow direction of the sea current, the velocity component of the sea current in the radar detection wave beam direction can cause the position of a first-order Bragg peak in a Doppler echo spectrum to generate frequency shift, the radial velocity component size in the radar radial direction is solved according to the frequency shift and the frequency dispersion relation (the relation between the fluctuation frequency and the wavelength), and the azimuth (the incoming wave direction) of the radial velocity is determined according to the space spectrum structure of an echo signal by utilizing a multiple signal classification algorithm (MUSIC algorithm). At present, the ocean current inversion algorithm is widely applied to the inversion of the radial flow detected by the high-frequency ground wave radar.
The synthesis algorithm for vector ocean currents is not completely mature, and two types of algorithms are generally adopted at present: one is that the vector flow velocity and the flow direction are solved by a single radar, and the solving method is based on the premise that the ocean flow velocity is kept stable or unchanged in a certain range, and the radial flow velocity information in a certain synthesis range is utilized to carry out vector synthesis by methods of least squares, flow function models and the like to obtain a vector flow field; the other type is that vector flow velocity synthesis is carried out through two radars, namely two radar beams are adopted to simultaneously radiate in the same sea area, and vector synthesis is carried out according to velocity components (namely a plurality of radial flows) in all directions, so that a vector flow field of the sea current of the public detection sea area is solved. At present, two radars are adopted to carry out vector flow synthesis algorithm application relatively commonly. But the synthetic vector flow is not subjected to vector marking, so that the low-quality synthetic vector flow cannot be deleted, the obtained ocean current data has low reliability, and the radar ocean current detection precision is low; because only the vector synthesis is carried out on the radial flow of the two radars, the radial flow of the multiple radars cannot be fused when the radars are networked, if one of the radars cannot work normally, the vector flow cannot be synthesized, and the stability of radar network ocean current detection is low.
In the method, in the process of inverting a vector flow field by a high-frequency ground wave radar system, a plurality of target radars (more than or equal to 3) jointly cover a target detection sea area, each target radar can obtain a radial velocity map, a plurality of radial velocity maps jointly covering the target detection sea area are fused and synthesized with each other, namely, every two or more radial velocity maps are fused to generate a vector flow graph, the synthesized vector flow data of each grid unit are subjected to quality evaluation and quality marking, the synthesized vector flow with the highest quality grade is selected as a final sea flow result, and the vector flow field jointly covering the target detection sea area is finally obtained. When marking, the quality factor is obtained by weighting calculation according to the synthetic geometric precision factor (GDOP) of the grid unit, the synthetic geometric precision factor (GDOP) is obtained by calculation through a spatial geometric relation in advance, and the weighting coefficient of the quality factor is preset. According to the method, the radial velocity map of the radar network is synthesized by multi-map fusion, so that the radial flow of all target radars which jointly cover a target detection sea area indirectly participates in the synthesis of a final vector flow field, and the relation between the final vector flow field precision and key factors influencing the detection precision, such as sea conditions, velocity spatial distribution, noise environment, radar parameters and the like, is indirectly established, so that sea flow data with high reliability can be obtained during sea flow detection, the sea flow detection precision of the radar is improved, and the stability of the radar network sea flow detection is improved.
When the method is applied to radar networking, the number of radar combinations is determined for a plurality of target radars which cover a target detection sea area together according to synthesis requirements, every two or more target radars in the plurality of target radars are used as one radar combination to carry out vector flow synthesis, a plurality of radar combinations exist, and the number of target radars included in one radar combinationnCan be designed adaptively, and the present invention is not particularly limited. For example, a target detection sea area has 5If two radars form a radar combination
And each two radial velocity maps of the radar combinations (namely 10 radar combinations) are fused and synthesized to generate a vector flow map. For example, if a target detection sea area is covered by 5 radars together, every three radars form a radar combination, then the target detection sea area exists
And each radar combination (namely 10 radar combinations) generates a vector flow graph by fusion and synthesis of every three radial velocity graphs. For example, if a target detection sea area is covered by 5 radars together, and every four radars form a radar combination, then the target detection sea area exists
And combining every four radial velocity maps to generate a vector flow map by combining the four radar combinations (namely 4 radar combinations).
In an alternative embodiment, the synthesis conditions comprise:
the grid units at the same position in each radial flow velocity diagram of the radar combination do not exceed the outer edge of the flow velocity area;
none of the grid cells at the same position in each radial velocity map of the radar assembly are located in the baseline instability region.
When the method of the invention carries out vector flow synthesis on the radial velocity maps of the radar combination, the judgment of the synthesis condition is carried out on the same grid unit (namely the grid unit at the same position) in each radial velocity map, so as to determine whether to synthesize the vector flow of the grid unit according to the judgment result. When the grid cell exceeds the edge of the flow velocity area (namely a fan-shaped area of a radial flow field generated by a single radar station), the vector flows of the grid cell are not synthesized, and vice versa. When a grid cell is located in a baseline unstable region (i.e. a region of plus or minus 5 degrees near the line connecting two radar stations), the vector streams of the grid cell are not synthesized, and vice versa.
In an alternative embodiment, for a grid cell in the baseline unstable region, if the grid cell does not have a composite vector stream, interpolation and completion are performed on the grid cell to obtain an ocean current result of the grid cell.
When the vector flow synthesis is carried out on the current radar combination, the grid units in the unstable baseline area are removed and are not synthesized, and the following two conditions exist:
(1) If the grid unit has the synthetic vector flow of other radar combinations, selecting the synthetic vector flow with the highest quality factor for the grid unit as the ocean current result of the grid unit;
(2) If the grid unit does not have the composite vector stream of other radar combinations, that is, the grid unit does not have the composite vector stream, interpolation and completion processing needs to be performed on the grid unit to obtain the ocean current result of the grid unit.
When the difference is filled, the values of a plurality of grid units around the grid unit are averaged and then are filled into the grid unit.
In an alternative embodiment, the first post-processing comprises:
and performing one or more of wild value elimination, flow field smoothing filtering, flow velocity with minimum variance and flow field interpolation on the synthetic vector flow of each grid unit.
In an alternative embodiment, the second post-treatment comprises:
and (3) performing one or more of outlier rejection, flow field smoothing filtering, minimum variance flow velocity taking and flow field interpolation on the ocean current result of each grid unit.
The process of the method of the present invention will be described in detail below.
S1, reading a radial stream data file, and obtaining the coverage area of each radar:
under a folder designated by a radar system, searching radial stream data files of all target radars (hereinafter referred to as stations) needing to participate in vector stream synthesis corresponding to the current moment, reading the data files, and obtaining the coverage area of each station;
s2, grid arrangement:
establishing a rectangular grid in a target detection sea area, wherein the grid resolution is 0.5 degrees x 0.5 degrees;
s3, initializing sea states:
initializing sea state information on each grid unit of the established rectangular grid, namely initializing the storage data of the grid unit, wherein the flow velocity and the flow direction are initialized to NaN, and the data quality is initialized to-1;
s4, determining a radar combination mode:
the total number of 5 stations jointly cover a target detection sea area, and the two stations are a radar combination (hereinafter referred to as double stations) and coexist in the target detection sea area
10 radar combinations, namely 10 radar combinations;
s5, acquiring a radial flow rate list of the stations participating in vector flow synthesis:
obtaining the radial flow velocity and the azimuth information of the double stations participating in the vector flow synthesis; the list of the radial flow velocity map which can be used for judging whether the synthesis condition is met in the next step is stored in the radial flow velocity list;
s6, judging whether the grid unit of the double-station radial flow velocity graph participating in vector flow synthesis in the radial flow velocity list meets the synthesis condition:
whether the grid unit exceeds the outer edge of the flow velocity area or not, if so, the grid unit is not synthesized, otherwise, the grid unit is synthesized;
whether the grid unit is located in the unstable baseline area or not, if so, the grid unit is not synthesized, otherwise, the grid unit is synthesized;
s7, vector flow synthesis:
carrying out vector flow synthesis on each grid unit meeting the synthesis condition, and meanwhile, calculating the synthesis quality according to the weighting of a synthesis geometric precision factor (GDOP);
s8, carrying out first post-processing on the synthesized vector flow field:
the method mainly comprises outlier rejection, flow field smooth filtering, minimum variance flow rate taking and flow field interpolation for the condition that a plurality of flow rates exist on a grid unit;
s9, establishing a candidate vector flow rate list:
adding the current composite vector flow of the double stations into a candidate vector flow rate list for subsequent screening;
s10, repeating S5-S9 to obtain a synthetic vector flow of all the double stations and adding the synthetic vector flow to a candidate radial flow rate list;
s11, taking the synthetic vector flow with the highest quality grade:
screening out the synthetic vector flow with the highest quality grade of each grid unit according to the quality factor as the ocean current result of the grid unit;
in S6, when the vector streams of the two stations are synthesized, the grid unit is not synthesized when the grid unit is located in the unstable area of the base line, and if the grid unit has the synthesized vector streams of other two stations, the synthesized vector stream with the highest quality factor is selected for the grid unit to serve as the ocean current result of the grid unit; if the grid unit does not have the synthetic vector flow of other double stations, interpolation and completion processing are required to be carried out on the grid unit to obtain the ocean current result of the grid unit;
s12, carrying out secondary post-processing on the screened ocean current result:
the method mainly comprises outlier rejection, flow field smooth filtering, minimum variance flow rate taking and flow field interpolation for the condition that a plurality of flow rates exist on a grid unit;
and S13, outputting and storing the result.
The system for synthesizing the vector flow of the radar radial flow comprises the following components:
the target radar determining module is used for determining a plurality of target radars which jointly cover a target detection sea area according to the coverage sea area of each radar, wherein the coverage sea area of each radar is obtained according to radial flow data corresponding to the current moment, and the number of the target radars is more than or equal to 3;
the system comprises a radial velocity map initialization module, a target detection sea area detection module and a sea state information initialization module, wherein the radial velocity map initialization module is used for acquiring a radial velocity map for each target radar, establishing a rectangular grid for the radial velocity map in a target detection sea area, and initializing sea state information on each grid unit, the rectangular grid comprises a plurality of grid units, and the sea state information comprises velocity and flow;
a radar combination determination module for determining the number of radar combinations
And forming a set of radar combinations, wherein m represents the number of the plurality of target radars, and n represents the number of target radars included in each radar combination;
the vector flow synthesis module is used for determining whether each grid unit meets synthesis conditions or not for the radial velocity map of the target radar contained in each radar combination in the set, synthesizing the vector flows of the grid units meeting the synthesis conditions, and performing first post-processing on the synthesized vector flow of each grid unit to obtain a vector flow map;
the quality mark evaluation module is used for marking a corresponding quality factor for the synthetic vector flow of each grid unit in each vector flow graph, and selecting the synthetic vector flow with the highest quality factor in a plurality of vector flow graphs as an ocean current result of the grid unit for the same grid unit, wherein the quality factor is obtained by weighting calculation according to the synthetic geometric precision factor of the grid unit;
and the result output module is used for carrying out secondary post-processing on the ocean current results of the grid units to obtain the final ocean current result of each grid unit, and obtaining the vector flow field in the target detection sea area based on the final ocean current result of each grid unit.
In an alternative embodiment, the synthesis conditions comprise:
grid cells at the same position in each radial flow velocity map of the radar combination do not exceed the outer edge of the flow velocity area;
none of the grid cells at the same position in each radial velocity map of the radar assembly are located in the baseline instability region.
In an alternative embodiment, for a grid cell in the baseline unstable region, if the grid cell does not have a synthetic vector stream, interpolation and completion are performed on the grid cell to obtain an ocean current result of the grid cell.
In an alternative embodiment, the first post-processing comprises:
and performing one or more of outlier elimination, flow field smoothing filtering, minimum variance flow velocity taking and flow field interpolation on the synthetic vector flow of each grid unit.
In an alternative embodiment, the second post-treatment comprises:
and (3) performing one or more of outlier rejection, flow field smoothing filtering, minimum variance flow velocity taking and flow field interpolation on the ocean current result of each grid unit.
In the description provided herein, numerous specific details are set forth. It is understood, however, that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.
Moreover, those of ordinary skill in the art will appreciate that although some embodiments described herein include some features included in other embodiments, not others, combinations of features of different embodiments are meant to be within the scope of the invention and form different embodiments. For example, in the claims, any of the claimed embodiments may be used in any combination.
It will be understood by those skilled in the art that while the present invention has been described with reference to exemplary embodiments, various changes may be made and equivalents may be substituted for elements thereof without departing from the scope of the invention. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from the essential scope thereof. Therefore, it is intended that the invention not be limited to the particular embodiment disclosed, but that the invention will include all embodiments falling within the scope of the appended claims.
Claims (10)
1. A method for vector flow synthesis of radar radial flow, the method comprising:
determining a plurality of target radars which jointly cover a target detection sea area according to the coverage sea area of each radar, wherein the coverage sea area of each radar is obtained according to radial flow data corresponding to the current moment, and the number of the target radars is more than or equal to 3;
acquiring a radial velocity map for each target radar, establishing a rectangular grid in a target detection sea area for the radial velocity map, and initializing sea state information on each grid unit, wherein the rectangular grid comprises a plurality of grid units, and the sea state information comprises velocity and flow;
determining the number of radar combinations
And forming a set of radar combinations, wherein,mrepresenting the number of the plurality of target radars,nrepresenting the number of target radars contained in each radar combination;
determining whether each grid unit meets a synthesis condition or not for a radial velocity map of a target radar contained in each radar combination in the set, synthesizing vector flows of the grid units meeting the synthesis condition, and performing first post-processing on the synthesized vector flow of each grid unit to obtain a vector flow graph;
for a quality factor corresponding to a synthetic vector flow mark of each grid unit in each vector flow graph, selecting a synthetic vector flow with the highest quality factor in a plurality of vector flow graphs as an ocean current result of the grid unit for the same grid unit, wherein the quality factor is obtained by weighting calculation according to a synthetic geometric precision factor of the grid unit;
and performing secondary post-processing on the ocean current results of the grid units to obtain a final ocean current result of each grid unit, and obtaining a vector flow field in the target detection sea area based on the final ocean current result of each grid unit.
2. The method of claim 1, wherein the synthesis conditions comprise:
grid cells at the same position in each radial flow velocity map of the radar combination do not exceed the outer edge of the flow velocity area;
none of the grid cells at the same position in each radial velocity map of the radar assembly are located in the baseline instability region.
3. The method of claim 2, wherein, for grid cells in a baseline unstable region,
and if the grid unit does not have the synthetic vector flow, carrying out interpolation and completion on the grid unit to obtain the ocean current result of the grid unit.
4. The method of claim 1, wherein the first post-processing comprises:
and performing one or more of wild value elimination, flow field smoothing filtering, flow velocity with minimum variance and flow field interpolation on the synthetic vector flow of each grid unit.
5. The method of claim 1, wherein the second post-processing comprises:
and (3) performing one or more of outlier rejection, flow field smoothing filtering, minimum variance flow velocity taking and flow field interpolation on the ocean current result of each grid unit.
6. A system for vector flow synthesis of radar radial flow, the system comprising:
the target radar determining module is used for determining a plurality of target radars which jointly cover a target detection sea area according to the coverage sea area of each radar, wherein the coverage sea area of each radar is obtained according to radial flow data corresponding to the current moment, and the number of the target radars is more than or equal to 3;
the system comprises a radial velocity map initialization module, a target detection sea area detection module and a sea state information initialization module, wherein the radial velocity map initialization module is used for acquiring a radial velocity map for each target radar, establishing a rectangular grid for the radial velocity map in a target detection sea area, and initializing sea state information on each grid unit, the rectangular grid comprises a plurality of grid units, and the sea state information comprises velocity and flow;
a radar combination determination module for determining the number of radar combinations
And forming a set of radar combinations, wherein m represents the number of the plurality of target radars, and n represents the number of target radars included in each radar combination;
the vector flow synthesis module is used for determining whether each grid unit meets the synthesis condition or not for the radial velocity map of the target radar contained in each radar combination in the set, synthesizing the vector flows of the grid units meeting the synthesis condition, and performing first post-processing on the synthesized vector flow of each grid unit to obtain a vector flow map;
the quality mark evaluation module is used for marking a corresponding quality factor for the synthetic vector flow of each grid unit in each vector flow graph, and selecting the synthetic vector flow with the highest quality factor in a plurality of vector flow graphs as an ocean current result of the grid unit for the same grid unit, wherein the quality factor is obtained by weighting calculation according to the synthetic geometric precision factor of the grid unit;
and the result output module is used for carrying out secondary post-processing on the ocean current results of the grid units to obtain the final ocean current result of each grid unit, and obtaining the vector flow field in the target detection sea area based on the final ocean current result of each grid unit.
7. The system of claim 6, wherein the synthesis conditions comprise:
the grid units at the same position in each radial flow velocity diagram of the radar combination do not exceed the outer edge of the flow velocity area;
none of the grid cells at the same position in each radial velocity map of the radar assembly are located in the baseline instability region.
8. The system of claim 7, wherein, for grid cells at a baseline instability region,
and if the grid unit does not have the synthetic vector flow, carrying out interpolation and completion on the grid unit to obtain the ocean current result of the grid unit.
9. The system of claim 6, wherein the first post-processing comprises:
and performing one or more of wild value elimination, flow field smoothing filtering, flow velocity with minimum variance and flow field interpolation on the synthetic vector flow of each grid unit.
10. The system of claim 6, wherein the second post-processing comprises:
and (3) performing one or more of outlier rejection, flow field smoothing filtering, minimum variance flow velocity taking and flow field interpolation on the ocean current result of each grid unit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210737877.3A CN114812710B (en) | 2022-06-28 | 2022-06-28 | Vector flow synthesis method and system for radar radial flow |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210737877.3A CN114812710B (en) | 2022-06-28 | 2022-06-28 | Vector flow synthesis method and system for radar radial flow |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114812710A CN114812710A (en) | 2022-07-29 |
| CN114812710B true CN114812710B (en) | 2022-10-04 |
Family
ID=82523062
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210737877.3A Active CN114812710B (en) | 2022-06-28 | 2022-06-28 | Vector flow synthesis method and system for radar radial flow |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114812710B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115877377B (en) * | 2022-12-12 | 2023-10-03 | 中山大学 | Radar networking vector flow field synthesis method, system, equipment and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103604944A (en) * | 2013-12-11 | 2014-02-26 | 哈尔滨工业大学 | Surface flow measurement method based on monostation shipborne high-frequency ground wave radar |
| CN111650581A (en) * | 2020-06-15 | 2020-09-11 | 南京莱斯电子设备有限公司 | Radar global target track automatic starting method based on environment perception |
| CN112946643A (en) * | 2021-01-28 | 2021-06-11 | 华东师范大学 | Offshore wind power extraction method and system based on time sequence radar remote sensing |
| CN113064129A (en) * | 2021-03-03 | 2021-07-02 | 湖北中南鹏力海洋探测系统工程有限公司 | High-frequency ground wave radar ocean current synthesis method |
| CN113466854A (en) * | 2021-06-29 | 2021-10-01 | 哈尔滨工业大学 | High-frequency ground wave radar inversion vector flow velocity method based on ocean power model |
Family Cites Families (1)
| 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 |
-
2022
- 2022-06-28 CN CN202210737877.3A patent/CN114812710B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103604944A (en) * | 2013-12-11 | 2014-02-26 | 哈尔滨工业大学 | Surface flow measurement method based on monostation shipborne high-frequency ground wave radar |
| CN111650581A (en) * | 2020-06-15 | 2020-09-11 | 南京莱斯电子设备有限公司 | Radar global target track automatic starting method based on environment perception |
| CN112946643A (en) * | 2021-01-28 | 2021-06-11 | 华东师范大学 | Offshore wind power extraction method and system based on time sequence radar remote sensing |
| CN113064129A (en) * | 2021-03-03 | 2021-07-02 | 湖北中南鹏力海洋探测系统工程有限公司 | High-frequency ground wave radar ocean current synthesis method |
| CN113466854A (en) * | 2021-06-29 | 2021-10-01 | 哈尔滨工业大学 | High-frequency ground wave radar inversion vector flow velocity method based on ocean power model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114812710A (en) | 2022-07-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Talvitie et al. | Distance-based interpolation and extrapolation methods for RSS-based localization with indoor wireless signals | |
| CN110557716B (en) | Indoor positioning method based on lognormal model | |
| EP2119264B1 (en) | System and method for generating a location estimate using a method of intersections | |
| US7395073B2 (en) | Method and apparatus for location estimation using region of confidence filtering | |
| KR101258668B1 (en) | Korea local radar processing system | |
| CN106454747B (en) | Wireless positioning method of mobile phone terminal | |
| CN114812710B (en) | Vector flow synthesis method and system for radar radial flow | |
| CN105338498A (en) | Construction method for fingerprint database in WiFi indoor positioning system | |
| CN112558076B (en) | Volume scanning mode calculation method based on networking weather radar coverage area and application | |
| US20150061927A1 (en) | Wide beam sar focusing method using navigation solution derived from autofocus data | |
| CN111968046B (en) | Target association fusion method for radar photoelectric sensor based on topological structure | |
| CN110187336B (en) | Multi-station radar site positioning and joint tracking method based on distributed PHD | |
| CN112152739A (en) | Method and device for calculating interference probability distribution of satellite constellation | |
| CN110488259B (en) | Radar target classification method and device based on GDBSCAN | |
| CN113810846B (en) | Indoor positioning method based on WiFi and IMU fusion | |
| CN115665668A (en) | Indoor and outdoor seamless positioning method, system, equipment and medium | |
| CN113933876A (en) | Multi-satellite communication time difference positioning data fusion processing method | |
| CN113589235B (en) | Radar radial flow data extraction method and system | |
| Fang et al. | LoRa-based Outdoor 3D Localization | |
| CN113038486A (en) | Neighbor cell planning method and device, computing equipment and computer storage medium | |
| CN114553270B (en) | Antenna broadcast beam configuration method, device, electronic equipment and storage medium | |
| CN103731917A (en) | WLAN (Wireless Local Area Network) positioning method for eliminating multi-antenna direction deviation | |
| CN111726751B (en) | Passive positioning-oriented direct wave sensor identification method | |
| Lecomte et al. | Cyclone Tracking with ERS-2 Scatterometer: Algorithm Performances and Post-processed Data Examples | |
| Kudas et al. | Analysis of the geometry of the TPI NETpro reference station network in Poland |
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 |