CN114488107B - Method and device for sea clutter space-time distribution and influence grading product manufacturing - Google Patents

Method and device for sea clutter space-time distribution and influence grading product manufacturing Download PDF

Info

Publication number
CN114488107B
CN114488107B CN202210382406.5A CN202210382406A CN114488107B CN 114488107 B CN114488107 B CN 114488107B CN 202210382406 A CN202210382406 A CN 202210382406A CN 114488107 B CN114488107 B CN 114488107B
Authority
CN
China
Prior art keywords
sea clutter
sea
calculation model
semi
time distribution
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
Application number
CN202210382406.5A
Other languages
Chinese (zh)
Other versions
CN114488107A (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.)
Southern Marine Science and Engineering Guangdong Laboratory Guangzhou
Aerospace Hongtu Information Technology Co Ltd
Original Assignee
Southern Marine Science and Engineering Guangdong Laboratory Guangzhou
Aerospace Hongtu Information Technology Co Ltd
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 Southern Marine Science and Engineering Guangdong Laboratory Guangzhou, Aerospace Hongtu Information Technology Co Ltd filed Critical Southern Marine Science and Engineering Guangdong Laboratory Guangzhou
Priority to CN202210382406.5A priority Critical patent/CN114488107B/en
Publication of CN114488107A publication Critical patent/CN114488107A/en
Application granted granted Critical
Publication of CN114488107B publication Critical patent/CN114488107B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/93Radar or analogous systems specially adapted for specific applications for anti-collision purposes
    • G01S13/937Radar or analogous systems specially adapted for specific applications for anti-collision purposes of marine craft
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/95Radar or analogous systems specially adapted for specific applications for meteorological use
    • G01S13/953Radar or analogous systems specially adapted for specific applications for meteorological use mounted on aircraft
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • G01S7/4052Means for monitoring or calibrating by simulation of echoes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details 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/414Discriminating targets with respect to background clutter
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information 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)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Algebra (AREA)
  • Ocean & Marine Engineering (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention provides a method and a device for sea clutter space-time distribution and influence grading product manufacture, which relate to the technical field of ocean engineering and comprise the following steps: acquiring target data of a target sea area and actually measured sea clutter data of the target sea area; based on the actually measured sea clutter data, carrying out deviation correction on the sea clutter semi-empirical calculation model to obtain a corrected sea clutter semi-empirical calculation model; determining a sea clutter space-time distribution product based on the target data, the corrected sea clutter semi-empirical calculation model and the composite Bragg backscattering calculation model; determining a grading threshold value of the sea clutter space-time distribution product based on the sea clutter space-time distribution product; based on the grading threshold, the sea clutter space-time distribution products are graded, the obtained sea clutter influences the grading products, and the technical problems that large-scale distribution trend analysis and influence efficiency evaluation on the sea clutter are difficult to perform are solved.

Description

Method and device for sea clutter space-time distribution and influence grading product manufacturing
Technical Field
The invention relates to the technical field of ocean engineering, in particular to a method and a device for sea clutter space-time distribution and influence on manufacturing of graded products.
Background
For surface target detection radars, sea surface backscattering of the transmitted signal often severely limits the radar's ability to detect ships, airplanes, missiles, navigation buoys, and other targets on the sea surface within a radar resolution unit. This interfering signal is commonly referred to as sea clutter. The sea clutter is closely related to the parameters of the radar (frequency, polarization, observation method, angle of ground friction) in addition to the parameters of the marine environment (such as sea surface wind field, sea wave, tide, etc.). In order to better understand the ocean interference information, the relationship between the sea clutter signals and radar parameters and ocean environment parameters needs to be established.
At present, aiming at the research of sea clutter characteristics, a test method is mostly adopted, and sea echo measurement data are obtained in a shore-based radar and airborne aviation test mode. The method can only obtain the sea clutter results in a certain area or a small range, and lacks analysis on the distribution trend of the large-scale sea clutter in the sea area. In addition, the current sea clutter backscattering coefficient calculation models can be divided into two types, one type is a semi-empirical sea clutter model GTI, TSC, HYB and NRL; one is an analytical approximation model, such as KA, SPM, a composite BRAG model. The application range of each model is limited, and a single model cannot be better used for calculating various application scenes.
No effective solution has been proposed to the above problems.
Disclosure of Invention
In view of this, the present invention provides a method and an apparatus for sea clutter spatial-temporal distribution and influence classification product manufacturing, so as to alleviate the technical problem that it is difficult to perform large-scale distribution trend analysis and influence efficiency evaluation on sea clutter in the prior art.
In a first aspect, an embodiment of the present invention provides a method for sea clutter spatial-temporal distribution and influence on hierarchical product manufacturing, including: acquiring target data of a target sea area and actually measured sea clutter data of the target sea area, wherein the target data comprises: sea condition data, radar observation parameters and range data of the target sea area; based on the actually measured sea clutter data, carrying out deviation correction on the sea clutter semi-empirical calculation model to obtain a corrected sea clutter semi-empirical calculation model; determining the sea clutter space-time distribution product based on the target data, the corrected sea clutter semi-empirical calculation model and a composite Bragg backscattering calculation model; determining a grading threshold value of the sea clutter space-time distribution product based on the sea clutter space-time distribution product; and grading the sea clutter space-time distribution products based on the grading threshold value to obtain the sea clutter influence grading products.
Further, the sea state data includes: sea surface wind speed and sea surface wind direction; the radar observation parameters comprise observation frequency, observation azimuth angle, ground wiping angle and polarization mode.
Further, based on the actually measured sea clutter data, the sea clutter semi-empirical calculation model is subjected to deviation correction to obtain a corrected sea clutter semi-empirical calculation model, and the method comprises the following steps: based on the actually measured sea clutter data, determining an observation value of a calculation parameter of sea clutter intensity and sea condition data corresponding to the observation value, wherein the calculation parameter of the sea clutter intensity comprises: observing frequency, a ground wiping angle, a relative wind direction and a polarization mode; calculating a theoretical value of a calculation parameter of the sea clutter intensity based on the sea condition data corresponding to the observation value and the sea clutter semi-empirical calculation model; and carrying out deviation correction on the sea clutter semi-empirical calculation model based on the observed value and the theoretical value to obtain the corrected sea clutter semi-empirical calculation model.
Further, determining the sea clutter space-time distribution product based on the target data, the corrected sea clutter semi-empirical calculation model and the composite bragg backscattering calculation model, comprising: determining a first ground wiping angle and a second ground wiping angle in the target data, wherein the first ground wiping angle is larger than or equal to a preset threshold value, and the second ground wiping angle is smaller than the preset threshold value; determining a sea clutter space-time distribution product of the first ground wiping angle based on the target data corresponding to the first ground wiping angle and the composite Bragg backscattering calculation model; and determining a sea clutter space-time distribution product of the second ground wiping corner based on the target data corresponding to the second ground wiping corner and the corrected sea clutter semi-empirical calculation model.
Further, based on the sea clutter space-time distribution product, determining a classification threshold of the sea clutter space-time distribution product, including: calculating the accumulative probability density distribution of the sea clutter space-time distribution product; determining cumulative probability segmentation points based on the cumulative probability density distribution; and determining a grading threshold value of the sea clutter space-time distribution product based on the accumulated probability segmentation points.
Further, based on the classification threshold, classifying the sea clutter space-time distribution product to obtain the sea clutter influence classification product, including: determining a sea clutter influence strength value of the sea clutter space-time distribution product; and grading the sea clutter space-time distribution products based on the sea clutter influence strength value and the grading threshold value to obtain the sea clutter influence grading products.
In a second aspect, an embodiment of the present invention further provides a system for sea clutter space-time distribution and influence grading product manufacturing, including: the acquisition unit, the calibration unit, sea clutter space-time distribution product generation unit, determine unit and sea clutter influence hierarchical product generation unit, wherein, the acquisition unit for obtain the target data of target sea area with the actual measurement sea clutter data of target sea area, wherein, the target data include: sea condition data, radar observation parameters and range data of the target sea area; the correction unit is used for carrying out deviation correction on the sea clutter semi-empirical calculation model based on the actually measured sea clutter data to obtain a corrected sea clutter semi-empirical calculation model; the sea clutter space-time distribution product generating unit is used for determining the sea clutter space-time distribution product based on the target data, the corrected sea clutter semi-empirical calculation model and the composite Bragg backscattering calculation model; the determining unit is used for determining a grading threshold value of the sea clutter space-time distribution product based on the sea clutter space-time distribution product; and the sea clutter influence grading product generating unit is used for grading the sea clutter space-time distribution products based on the grading threshold value to obtain the sea clutter influence grading products.
Further, the sea state data includes: sea surface wind speed and sea surface wind direction; the radar observation parameters comprise observation frequency, observation azimuth angle, ground wiping angle and polarization mode. In a third aspect, an embodiment of the present invention further provides an electronic device, including a memory and a processor, where the memory is used to store a program that supports the processor to execute the method in the first aspect, and the processor is configured to execute the program stored in the memory.
In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored.
In the embodiment of the present invention, target data of a target sea area and actually measured sea clutter data of the target sea area are obtained, where the target data includes: sea condition data, radar observation parameters and range data of the target sea area; based on the actually measured sea clutter data, carrying out deviation correction on the sea clutter semi-empirical calculation model to obtain a corrected sea clutter semi-empirical calculation model; determining the sea clutter space-time distribution product based on the target data, the corrected sea clutter semi-empirical calculation model and a composite Bragg backscattering calculation model; determining a grading threshold value of the sea clutter space-time distribution product based on the sea clutter space-time distribution product; based on the grading threshold, the sea clutter space-time distribution product is graded to obtain the sea clutter influence grading product, the purpose of generating the sea clutter space-time distribution product and the sea clutter influence grading product is achieved, the technical problems that the sea clutter is difficult to carry out large-scale distribution trend analysis and influence efficiency evaluation in the prior art are solved, and therefore the technical effects of large-scale distribution trend analysis and influence efficiency evaluation on the sea clutter are achieved.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
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, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of a method for sea clutter spatial-temporal distribution and influence classification product manufacturing according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a system for sea clutter spatial-temporal distribution and influence classification product manufacturing according to an embodiment of the present invention;
fig. 3 is a schematic diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some, but not all embodiments of the present invention. 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.
The first embodiment is as follows:
in accordance with an embodiment of the present invention, there is provided an embodiment of a method for sea clutter spatio-temporal distribution and influence hierarchical product production, it being noted that the steps illustrated in the flow charts of the accompanying drawings may be executed in a computer system such as a set of computer executable instructions and that, although a logical order is illustrated in the flow charts, in some cases, the steps illustrated or described may be executed in an order different than here.
Fig. 1 is a flow chart of a method for sea clutter spatial-temporal distribution and influence grading product manufacturing according to an embodiment of the present invention, as shown in fig. 1, the method includes the following steps:
step S102, obtaining target data of a target sea area and actually measured sea clutter data of the target sea area, wherein the target data comprises: sea condition data, radar observation parameters and range data of the target sea area;
the sea state data includes: sea surface wind speed and sea surface wind direction; the radar observation parameters comprise observation frequency, observation azimuth angle, ground wiping angle and polarization mode.
Step S104, based on the actually measured sea clutter data, carrying out deviation correction on a sea clutter semi-empirical calculation model to obtain a corrected sea clutter semi-empirical calculation model;
step S106, determining a sea clutter space-time distribution product based on the target data, the corrected sea clutter semi-empirical calculation model and a composite Bragg backscattering calculation model;
step S108, determining a grading threshold value of the sea clutter space-time distribution product based on the sea clutter space-time distribution product;
and step S110, grading the sea clutter space-time distribution products based on the grading threshold value to obtain the sea clutter influence grading products.
In the embodiment of the present invention, target data of a target sea area and actually measured sea clutter data of the target sea area are obtained, where the target data includes: sea condition data, radar observation parameters and range data of the target sea area; based on the actually measured sea clutter data, carrying out deviation correction on a sea clutter semi-empirical calculation model to obtain a corrected sea clutter semi-empirical calculation model; determining the sea clutter space-time distribution product based on the target data, the corrected sea clutter semi-empirical calculation model and a composite Bragg backscattering calculation model; determining a grading threshold value of the sea clutter space-time distribution product based on the sea clutter space-time distribution product; based on the grading threshold, the sea clutter space-time distribution product is graded to obtain the sea clutter influence grading product, the purpose of generating the sea clutter space-time distribution product and the sea clutter influence grading product is achieved, the technical problems that the sea clutter is difficult to carry out large-scale distribution trend analysis and influence efficiency evaluation in the prior art are solved, and therefore the technical effects of large-scale distribution trend analysis and influence efficiency evaluation on the sea clutter are achieved.
In the embodiment of the present invention, step S104 includes the following steps:
step S11, based on the actually measured sea clutter data, determining an observation value of a calculation parameter of a sea clutter intensity and sea condition data corresponding to the observation value, where the calculation parameter of the sea clutter intensity includes: observing frequency, a ground wiping angle, a relative wind direction and a polarization mode;
step S12, calculating a theoretical value of a calculation parameter of the sea clutter intensity based on the sea condition data corresponding to the observation value and the sea clutter semi-empirical calculation model;
and step S13, performing deviation correction on the sea clutter semi-empirical calculation model based on the observed value and the theoretical value to obtain the corrected sea clutter semi-empirical calculation model.
It should be noted that the observation relative wind direction is equal to the included angle between the radar observation azimuth and the sea surface wind direction, and is generally divided into three categories, wherein 0 ° represents downwind observation, 90 ° represents anemometry observation, and 180 ° represents upwind observation. The ground wiping angle can be generally divided into 1 degree, 10 degrees, 30 degrees and 60 degrees; radar observation parameters are divided into L, S, C and Ku wave bands; polarization modes are generally classified into two types, HH and VV.
In the embodiment of the invention, preferably, the actually measured sea clutter data selects a Nathansan sea clutter data table, the data table comprises backscatter coefficient measured values corresponding to seven wave bands, seven ground rubbing angles, seven sea conditions and two polarization modes, and the data table is a relatively complete public database of sea clutter scattering rates up to now. The data sheet does not distinguish specific observed wind directions, and records average values of sea clutter measurement under downwind, upwind and anemometry observation. The semi-empirical sea clutter model selects the Technology Service Corporation model (TSC). The TSC model is provided on the basis of summarizing measured average backscattering coefficient data and mainly comprises three influence factors, namely a small ground-rubbing angle, a wind speed factor and a wind direction factor; meanwhile, the influence of abnormal propagation conditions such as evaporation waveguide is considered, and the scattering coefficient obtained by the model changes more smoothly. When the propagation condition is unknown, the TSC model can better describe the sea surface condition and predict the backscattering coefficient.
The calculation formula of the deviation correction amount between the observed value and the theoretical value is as follows:
Figure P_220413090802873_873047001
in the formula (I), the compound is shown in the specification,
Figure P_220413090802935_935553001
in order to observe the frequency of the signal,
Figure P_220413090802951_951178002
in order to wipe the ground corner, the utility model is provided with a brush,
Figure P_220413090802966_966808003
in the form of a polarization mode, the polarization mode,
Figure P_220413090802998_998061004
the relative observation of the downwind, the upwind and the anemometry is carried out, the values are 0, 90 degrees and 180 degrees,
Figure P_220413090803013_013707005
in order to calculate the model, the model is calculated,
Figure P_220413090803029_029310006
is shown as a data table of Nathansan,
Figure P_220413090803060_060557007
the calculated deviation correction amounts at the set of frequencies, angles of incidence, and polarizations.
The calculation formula of the TSM semi-empirical sea clutter model is as follows:
Figure P_220413090803076_076182001
in the formula (I), the compound is shown in the specification,
Figure P_220413090803107_107469001
the backscattering coefficient obtained for the HH polarization mode,
Figure P_220413090803138_138682002
the backscattering coefficient obtained for the VV polarization mode,
Figure P_220413090803154_154315003
is the angle of the ground, the unit is radian,
Figure P_220413090803169_169937004
a factor is calculated for the small floor scrub angle,
Figure P_220413090803201_201210005
is a factor of the wind direction,
Figure P_220413090803216_216812006
as a factor of the wind speed,
Figure P_220413090803232_232420007
is a measure of wavelength, in feet,
Figure P_220413090803253_253394008
the standard deviation of sea level height, in feet,
Figure P_220413090803285_285150009
to observe the frequency.
And subtracting the observed value and the theoretical value to obtain a deviation correction value, and compensating the deviation correction value in the sea clutter semi-empirical calculation model as an error term to obtain a corrected sea clutter semi-empirical calculation model.
In the embodiment of the present invention, step S106 includes the following steps:
step S21, determining a first ground wiping angle and a second ground wiping angle in the target data, wherein the first ground wiping angle is greater than or equal to a preset threshold value, and the second ground wiping angle is smaller than the preset threshold value;
step S22, determining a sea clutter space-time distribution product of the first ground wiping angle based on the target data corresponding to the first ground wiping angle and the composite Bragg backscattering calculation model;
and step S23, determining a sea clutter space-time distribution product of the second ground rubbing corner based on the target data corresponding to the second ground rubbing corner and the corrected sea clutter semi-empirical calculation model.
In the embodiment of the invention, the size of the ground rubbing angle in the target data is judged, and then the corresponding calculation model is selected. And if the ground wiping angle is smaller than 10 degrees (a preset threshold value), the ground wiping angle is a first ground wiping angle, a corrected sea clutter semi-empirical calculation model is selected, and a sea clutter space-time distribution product of the first ground wiping angle is determined by using the corrected sea clutter semi-empirical calculation model and target data corresponding to the first ground wiping angle.
And if the ground clearance angle is larger than or equal to 10 degrees (a preset threshold value), the ground clearance angle is a second ground clearance angle, a composite Bragg backscattering calculation model is selected, and a sea clutter space-time distribution product of the second ground clearance angle is determined by utilizing the composite Bragg backscattering calculation model and target data corresponding to the second ground clearance angle.
It should be noted that the corrected sea clutter semi-empirical calculation model
Figure P_220413090803300_300802001
The calculation formula is described as follows:
Figure P_220413090803316_316413001
in the formula (I), the compound is shown in the specification,
Figure P_220413090803347_347675001
calculating a model for the corrected sea clutter semi-experience;
Figure P_220413090803363_363299002
is seaA clutter semi-empirical calculation model is obtained,
Figure P_220413090803394_394565003
to compensate for the deviation.
The composite Bragg scattering model takes the surface roughness as the superposition of two kinds of roughness, and the scattering coefficient calculated by the traditional Bragg scattering is subjected to ensemble averaging on the slope distribution of the large roughness.
The calculation formula of the composite bragg scattering model is described as follows:
Figure P_220413090803410_410164001
in the formula (I), the compound is shown in the specification,
Figure P_220413090803442_442362001
the wave number is observed by the radar,
Figure P_220413090803474_474184002
is the incident angle of the radar wave,
Figure P_220413090803489_489744003
is the wave number spectrum of the sea surface waves,
Figure P_220413090803505_505380004
with respect to the direction of observation,
Figure P_220413090803521_521015005
subscript for normalized backscattering coefficient
Figure P_220413090803552_552263006
And
Figure P_220413090803567_567900007
represents the polarization (H or V),
Figure P_220413090803583_583532008
for the local angle of incidence,
Figure P_220413090803614_614772009
and
Figure P_220413090803630_630385010
the dip angles of the sea surface slope in the x and y directions respectively.
In the embodiment of the present invention, step S108 includes the following steps:
step S31, calculating the cumulative probability density distribution of the sea clutter space-time distribution product;
step S32, determining cumulative probability dividing points based on the cumulative probability density distribution;
and step S33, determining the grading threshold value of the sea clutter space-time distribution product based on the accumulated probability segmentation point.
Preferably, the predetermined time period is generally 12 months.
In the embodiment of the invention, firstly, the accumulative probability density distribution of the sea clutter space-time distribution product is calculated.
Then, cumulative probability partition points are determined, and in the embodiment of the present invention, the cumulative probability partition points are two, which are 0.5 and 0.9 respectively.
And finally, determining a grading threshold value of the sea clutter space-time distribution product according to the accumulated probability segmentation points.
The segmentation thresholds A and B are-33.1079 and-30.2807 at cumulative probability segmentation points of 0.5 and 0.9, respectively.
In the embodiment of the present invention, step S110 includes the following steps:
step S41, determining a sea clutter influence strength value of the sea clutter space-time distribution product;
and step S42, grading the sea clutter space-time distribution products based on the sea clutter influence strength value and the grading threshold value to obtain the sea clutter influence grading products.
In the embodiment of the invention, the sea clutter space-time distribution product is utilized, and in the sea clutter space-time distribution product, an influence area corresponding to the sea clutter space-time distribution product with a sea clutter influence strength value smaller than a threshold value A is defined and determined as a weak influence area; determining an influence area corresponding to the sea clutter time-space distribution product with the sea clutter influence strength value larger than the threshold value B as a strong influence area; and determining the influence areas corresponding to the other sea clutter space-time distribution products as medium influence areas.
The embodiment of the invention breaks through the limitation that the traditional sea clutter analysis is mostly based on in-place detection data, simultaneously provides a sea clutter combination calculation model with multiple incidence angles and an influence grading model based on amplitude distribution characteristics, and is favorable for further analysis of large-scale distribution trend of the sea clutter and evaluation of influence efficiency of the sea clutter.
The second embodiment:
the embodiment of the invention also provides a system for sea clutter space-time distribution and influence on graded product manufacturing, which is used for executing the method for sea clutter space-time distribution and influence on graded product manufacturing provided by the embodiment of the invention.
As shown in fig. 2, fig. 2 is a schematic diagram of the system for sea clutter spatial-temporal distribution and influence of hierarchical product manufacturing, where the system for sea clutter spatial-temporal distribution and influence of hierarchical product manufacturing includes: the sea clutter detection system comprises an acquisition unit 10, a correction unit 20, a sea clutter space-time distribution product generation unit 30, a determination unit 40 and a sea clutter influence grading product generation unit 50.
The acquiring unit 10 is configured to acquire target data of a target sea area and actually measured sea clutter data of the target sea area, where the target data includes: sea condition data, radar observation parameters and range data of the target sea area;
the correction unit 20 is configured to perform deviation correction on the sea clutter semi-empirical calculation model based on the actually measured sea clutter data to obtain a corrected sea clutter semi-empirical calculation model;
the sea clutter space-time distribution product generating unit 30 is configured to determine the sea clutter space-time distribution product based on the target data, the corrected sea clutter semi-empirical calculation model, and the composite bragg backscattering calculation model;
the determining unit 40 is configured to determine a classification threshold of the sea clutter space-time distribution product based on the sea clutter space-time distribution product;
and the sea clutter influence grading product generating unit 50 is used for grading the sea clutter space-time distribution products based on the grading threshold value to obtain the sea clutter influence grading products.
In the embodiment of the present invention, target data of a target sea area and actually measured sea clutter data of the target sea area are obtained, where the target data includes: sea condition data, radar observation parameters and range data of the target sea area; based on the actually measured sea clutter data, carrying out deviation correction on the sea clutter semi-empirical calculation model to obtain a corrected sea clutter semi-empirical calculation model; determining a sea clutter spatial-temporal distribution product based on the target data, the corrected sea clutter semi-empirical calculation model and a composite Bragg backscattering calculation model; determining a grading threshold value of the sea clutter space-time distribution product based on the sea clutter space-time distribution product; based on the grading threshold, the sea clutter space-time distribution product is graded to obtain the sea clutter influence grading product, the purpose of generating the sea clutter space-time distribution product and the sea clutter influence grading product is achieved, the technical problems that the sea clutter is difficult to carry out large-scale distribution trend analysis and influence efficiency evaluation in the prior art are solved, and therefore the technical effects of large-scale distribution trend analysis and influence efficiency evaluation on the sea clutter are achieved.
Preferably, the sea state data includes: sea surface wind speed and sea surface wind direction; the radar observation parameters comprise observation frequency, observation azimuth angle, ground wiping angle and polarization mode.
Example three:
an embodiment of the present invention further provides an electronic device, including a memory and a processor, where the memory is used to store a program that supports the processor to execute the method described in the first embodiment, and the processor is configured to execute the program stored in the memory.
Referring to fig. 3, an embodiment of the present invention further provides an electronic device 100, including: a processor 60, a memory 61, a bus 62 and a communication interface 63, wherein the processor 60, the communication interface 63 and the memory 61 are connected through the bus 62; the processor 60 is arranged to execute executable modules, such as computer programs, stored in the memory 61.
The Memory 61 may include a high-speed Random Access Memory (RAM) and may also include a non-volatile Memory (non-volatile Memory), such as at least one disk Memory. The communication connection between the network element of the system and at least one other network element is realized through at least one communication interface 63 (which may be wired or wireless), and the internet, a wide area network, a local network, a metropolitan area network, and the like can be used.
The bus 62 may be an ISA bus, PCI bus, EISA bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one double-headed arrow is shown in FIG. 3, but this does not indicate only one bus or one type of bus.
The memory 61 is used for storing a program, the processor 60 executes the program after receiving an execution instruction, and the method executed by the apparatus defined by the flow process disclosed in any of the foregoing embodiments of the present invention may be applied to the processor 60, or implemented by the processor 60.
The processor 60 may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware or instructions in the form of software in the processor 60. The Processor 60 may be a general-purpose Processor, and includes a Central Processing Unit (CPU), a Network Processor (NP), and the like; the Integrated Circuit may also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field-Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic device, or discrete hardware components. The various methods, steps and logic blocks disclosed in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The storage medium is located in a memory 61, and the processor 60 reads the information in the memory 61 and completes the steps of the method in combination with the hardware.
Example four:
the embodiment of the present invention further provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the computer program performs the steps of the method in the first embodiment.
In addition, in the description of the embodiments of the present invention, unless otherwise explicitly specified or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc. indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, and are only for convenience of description and simplification of description, but do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus, should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and there may be other divisions when actually implemented, and for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection of devices or units through some communication interfaces, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: any person skilled in the art can modify or easily conceive the technical solutions described in the foregoing embodiments or equivalent substitutes for some technical features within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (10)

1. A method for sea clutter space-time distribution and influence grading product manufacture is characterized by comprising the following steps:
acquiring target data of a target sea area and actually measured sea clutter data of the target sea area, wherein the target data comprises: sea condition data, radar observation parameters and range data of the target sea area;
based on the actually measured sea clutter data, carrying out deviation correction on the sea clutter semi-empirical calculation model to obtain a corrected sea clutter semi-empirical calculation model;
determining a sea clutter spatial-temporal distribution product based on the target data, the corrected sea clutter semi-empirical calculation model and a composite Bragg backscattering calculation model;
determining a grading threshold value of the sea clutter space-time distribution product based on the sea clutter space-time distribution product;
grading the sea clutter space-time distribution product based on the grading threshold value to obtain a sea clutter influence grading product;
the sea clutter semi-empirical calculation model is a TSC model, and the corrected sea clutter semi-empirical calculation model is a TSC model after deviation correction is completed;
the calculation formula of the composite bragg scattering model is described as follows:
Figure P_220608173614564_564448001
in the formula (I), the compound is shown in the specification,
Figure F_220608173607062_062512001
the wave number is observed by the radar,
Figure F_220608173607158_158657002
is the incident angle of the radar wave,
Figure F_220608173607252_252452003
is the wave number spectrum of the sea surface waves,
Figure F_220608173607367_367187004
with respect to the direction of observation,
Figure P_220608173614611_611344005
subscript for normalized backscattering coefficient
Figure F_220608173607476_476551005
And
Figure F_220608173607575_575184006
respectively representing the H-polarization and the V-polarization,
Figure F_220608173607668_668954007
for the local angle of incidence,
Figure F_220608173607765_765631008
and
Figure F_220608173607859_859423009
the dip angles of the sea surface slope in the x and y directions respectively.
2. The method of claim 1,
the sea state data includes: sea surface wind speed and sea surface wind direction;
the radar observation parameters comprise observation frequency, observation azimuth angle, ground wiping angle and polarization mode.
3. The method of claim 2, wherein the step of performing a bias correction on the sea clutter semi-empirical calculation model based on the measured sea clutter data to obtain a corrected sea clutter semi-empirical calculation model comprises:
based on the actually measured sea clutter data, determining an observation value of a calculation parameter of sea clutter intensity and sea condition data corresponding to the observation value, wherein the calculation parameter of the sea clutter intensity comprises: observing frequency, a ground wiping angle, a relative wind direction and a polarization mode;
calculating a theoretical value of a calculation parameter of the sea clutter intensity based on the sea condition data corresponding to the observation value and the sea clutter semi-empirical calculation model;
and carrying out deviation correction on the sea clutter semi-empirical calculation model based on the observed value and the theoretical value to obtain the corrected sea clutter semi-empirical calculation model.
4. The method of claim 1, wherein determining the sea clutter spatiotemporal distribution product based on the target data, the corrected sea clutter semi-empirical calculation model, and a composite bragg backscattering calculation model comprises:
determining a first ground wiping angle and a second ground wiping angle in the target data, wherein the first ground wiping angle is larger than or equal to a preset threshold value, and the second ground wiping angle is smaller than the preset threshold value;
determining a sea clutter space-time distribution product of the first ground wiping angle based on the target data corresponding to the first ground wiping angle and the composite Bragg backscattering calculation model;
and determining a sea clutter space-time distribution product of the second ground wiping corner based on the target data corresponding to the second ground wiping corner and the corrected sea clutter semi-empirical calculation model.
5. The method of claim 1, wherein determining a classification threshold for the sea clutter spatial-temporal distribution product based on the sea clutter spatial-temporal distribution product comprises:
calculating the accumulative probability density distribution of the sea clutter space-time distribution product;
determining cumulative probability segmentation points based on the cumulative probability density distribution;
and determining a grading threshold value of the sea clutter space-time distribution product based on the accumulated probability segmentation points.
6. The method of claim 1, wherein ranking the sea clutter spatiotemporal distribution product based on the ranking threshold to obtain the sea clutter impact ranking product comprises:
determining a sea clutter influence strength value of the sea clutter space-time distribution product;
and grading the sea clutter space-time distribution products based on the sea clutter influence strength value and the grading threshold value to obtain the sea clutter influence grading products.
7. A system for sea clutter spatial-temporal distribution and influence of graded product production, comprising: an acquisition unit, a correction unit, a sea clutter space-time distribution product generation unit, a determination unit and a sea clutter influence grading product generation unit, wherein,
the acquisition unit is used for acquiring target data of a target sea area and actually measured sea clutter data of the target sea area, wherein the target data comprises: sea condition data, radar observation parameters and range data of the target sea area;
the correction unit is used for carrying out deviation correction on the sea clutter semi-empirical calculation model based on the actually measured sea clutter data to obtain a corrected sea clutter semi-empirical calculation model;
the sea clutter spatial-temporal distribution product generating unit is used for determining a sea clutter spatial-temporal distribution product based on the target data, the corrected sea clutter semi-empirical calculation model and the composite Bragg backscattering calculation model;
the determining unit is used for determining a grading threshold value of the sea clutter space-time distribution product based on the sea clutter space-time distribution product;
the sea clutter influence grading product generating unit is used for grading the sea clutter space-time distribution product based on the grading threshold value to obtain a sea clutter influence grading product;
the sea clutter semi-empirical calculation model is a TSC model, and the corrected sea clutter semi-empirical calculation model is a TSC model after deviation correction is completed;
the calculation formula of the composite bragg scattering model is described as follows:
Figure P_220608173614642_642595001
in the formula (I), the compound is shown in the specification,
Figure F_220608173607954_954554010
the wave number is observed by the radar,
Figure F_220608173608048_048843011
is the incident angle of the radar wave,
Figure F_220608173608127_127274012
is the wave number spectrum of the sea surface waves,
Figure F_220608173608223_223142013
with respect to the direction of observation,
Figure P_220608173614673_673853005
subscript for normalized backscattering coefficient
Figure F_220608173608316_316887014
And
Figure F_220608173608397_397446015
respectively representing the H-polarization and the V-polarization,
Figure F_220608173608491_491211016
for the local angle of incidence,
Figure F_220608173608587_587918017
and
Figure F_220608173608681_681632018
the dip angles of the sea surface slope in the x and y directions respectively.
8. The system of claim 7,
the sea state data includes: sea surface wind speed and sea surface wind direction;
the radar observation parameters comprise observation frequency, observation azimuth angle, ground wiping angle and polarization mode.
9. An electronic device comprising a memory for storing a program that enables a processor to perform the method of any of claims 1 to 6 and a processor configured to execute the program stored in the memory.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method according to any one of the claims 1 to 6.
CN202210382406.5A 2022-04-13 2022-04-13 Method and device for sea clutter space-time distribution and influence grading product manufacturing Active CN114488107B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210382406.5A CN114488107B (en) 2022-04-13 2022-04-13 Method and device for sea clutter space-time distribution and influence grading product manufacturing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210382406.5A CN114488107B (en) 2022-04-13 2022-04-13 Method and device for sea clutter space-time distribution and influence grading product manufacturing

Publications (2)

Publication Number Publication Date
CN114488107A CN114488107A (en) 2022-05-13
CN114488107B true CN114488107B (en) 2022-07-19

Family

ID=81488554

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210382406.5A Active CN114488107B (en) 2022-04-13 2022-04-13 Method and device for sea clutter space-time distribution and influence grading product manufacturing

Country Status (1)

Country Link
CN (1) CN114488107B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116976091B (en) * 2023-06-29 2024-03-01 中国空间技术研究院 Wave environment factor influence analysis method suitable for microwave remote sensing satellite

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7286079B2 (en) * 2005-04-21 2007-10-23 The United States Of America As Represented By The Secretary Of The Navy Method and apparatus for detecting slow-moving targets in high-resolution sea clutter
JP2011117899A (en) * 2009-12-07 2011-06-16 Mitsubishi Electric Corp Radar device
CN107678003B (en) * 2017-09-15 2021-04-06 国家***第一海洋研究所 Target detection method under ground wave radar sea clutter background
CN107490790B (en) * 2017-10-10 2021-06-11 北京航空航天大学 Simulation method of continuous multi-pulse coherent sea clutter
CN109143182B (en) * 2018-09-10 2019-06-07 南京莱斯电子设备有限公司 A kind of multiple dimensioned radar block plan immediate processing method
CN110441752A (en) * 2019-07-30 2019-11-12 中国人民解放军海军航空大学 Method for detecting targets at sea and device
CN110907907B (en) * 2019-10-19 2022-06-14 中国电波传播研究所(中国电子科技集团公司第二十二研究所) Sea clutter Doppler spectrum characteristic analysis and comparison method
CN111308433A (en) * 2020-03-26 2020-06-19 北京航空航天大学 Sea clutter simulation method containing texture information
CN113505463A (en) * 2021-02-10 2021-10-15 北京理工大学 Sea clutter and angle repeated die assembly type, construction system and analog simulation method
CN113065221B (en) * 2021-02-10 2023-02-17 北京理工大学 Large ground-wiping corner undersea clutter model, model construction system and simulation method
CN113281709B (en) * 2021-04-21 2022-08-02 中国海洋大学 Radar performance evaluation method based on area coupling forecasting system
CN113189561B (en) * 2021-06-16 2023-12-15 中国人民解放***箭军工程大学 Sea clutter parameter estimation method, system, equipment and storage medium

Also Published As

Publication number Publication date
CN114488107A (en) 2022-05-13

Similar Documents

Publication Publication Date Title
US9097805B2 (en) Systems and methods for calibrating dual polarization radar systems
Cho et al. Identification and removal of ground echoes and anomalous propagation using the characteristics of radar echoes
US9261583B2 (en) Systems and methods for calibrating dual polarization radar systems
CN110095760A (en) A kind of test device and method of weather radar
Richardson et al. Bragg scatter detection by the WSR-88D. Part I: Algorithm development
CN103064073A (en) Method based on frequency agility for changing radar target properties
CN114488107B (en) Method and device for sea clutter space-time distribution and influence grading product manufacturing
CN113281709B (en) Radar performance evaluation method based on area coupling forecasting system
CN109581383B (en) Anisotropic sea surface radar backscattering simulation method and system
Myagkov et al. Evaluation of the reflectivity calibration of W-band radars based on observations in rain
May et al. Evaluation of microphysical retrievals from polarimetric radar with wind profiler data
CN113219466A (en) Sea surface wind speed determination method and device, electronic equipment and storage medium
Anderson et al. Analysis of ASCAT ocean backscatter measurement noise
CN115048952A (en) GNSS-IR soil humidity inversion method integrating robust estimation and machine learning
CN110596666B (en) Deviation calibration method and device for radar differential reflectivity
CN109541565B (en) Radar echo intensity uniformity detection method and system
Asai et al. Validation of X-band multiparameter phased-array weather radar by comparing data from Doppler weather radar with a parabolic dish antenna
Wyatt et al. Factors affecting the accuracy of SHOWEX HF radar wave measurements
Overeem et al. Full-year evaluation of nonmeteorological echo removal with dual-polarization fuzzy logic for two C-band radars in a temperate climate
Huuskonen et al. Radar performance monitoring using the angular width of the solar image
CN114814779B (en) Buoy surge wave height observation data error evaluation method, system, equipment and medium
JP2008157713A (en) Radar rain gage operation control system
CN113065221B (en) Large ground-wiping corner undersea clutter model, model construction system and simulation method
Rosenberg et al. Characterisation of the tri-modal discrete sea clutter model
CN113505463A (en) Sea clutter and angle repeated die assembly type, construction system and analog simulation method

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