CN107607955B - A kind of screening plant of data screening method and ground clutter - Google Patents
A kind of screening plant of data screening method and ground clutter Download PDFInfo
- Publication number
- CN107607955B CN107607955B CN201710827028.6A CN201710827028A CN107607955B CN 107607955 B CN107607955 B CN 107607955B CN 201710827028 A CN201710827028 A CN 201710827028A CN 107607955 B CN107607955 B CN 107607955B
- Authority
- CN
- China
- Prior art keywords
- scanning
- ground clutter
- range bin
- pots
- absolute phase
- 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
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The embodiment of the invention discloses a kind of data screening method and the screening plants of ground clutter, select suitable ground clutter when calculating air index for radar.The method comprise the steps that in the M period of acquisition, carry out one group of plan-position scanning to the identical elevation angle within each period, one group of scanning includes primary high PRF scanning and primary low PRF scanning；Target data is obtained according to the radar echo signal returned after each scanning, target data includes at least ground clutter figure CF and absolute phase POTS, when it is ground clutter that the CF of range bin to be selected, which identifies range bin to be selected, calculate the first absolute phase sequence of range bin to be selected and the correlation of the second absolute phase sequence, first absolute phase sequence is that the set for the POTS that the high PRF scanning obtains is carried out in M period, and the second absolute phase sequence is that the set for the PORTS that low PRF scanning obtains is carried out in M period；When correlation is greater than the first preset threshold value, it is determined that range bin to be selected is target range library.
Description
Technical field
This application involves a kind of screening of Radar Technology field more particularly to data screening method and ground clutter dresses
It sets.
Background technique
The working principle of radar is to realize target acquisition, either aircraft, guided missile using the reflected echo of object, also
It is that ground, wave, sexual intercourse can all generate echo, therefore the difference of the occasion according to application, the title of radar is also different.Work as radar
When applied to weather forecast, referred to as weather radar.The main observation object of weather radar is precipitation echo, when weather radar because
The layer knot effect of atmosphere, when lower atmosphere density is bigger than high-altitude density, ray can occur to reflect upwards.When at dusk or
After person's rain, it is possible that the temperature of lower atmosphere is lower than the temperature of upper atmosphere, and there is inversion, so as to cause superrefraction production
Raw, occurring the characteristics of ground clutter, ground clutter on radar at this time is generated at the low elevation angle.
The phase of echo variation that weather radar can use object calculates the variation of air index, and then inverting obtains
Steam change information near the ground, to provide important data reference for Synoptic Analysis.
Currently, weather radar can change the atmospheric refraction calculated near weather radar according to the absolute phase of ground clutter
Rate variation, but not every ground clutter is suitable for calculating air index, needs to screen ground clutter.
Summary of the invention
The embodiment of the invention provides a kind of data screening method and the screening plants of ground clutter, by radar based on
Suitable ground clutter is selected when calculating air index.
The first aspect of the embodiment of the present application provides a kind of data screening method, comprising:
M period is acquired, one group of plan-position scanning, one group of scanning are carried out to the identical elevation angle within each period
Including a high pulse repetition frequency PRF scanning and primary low PRF scanning, the M is positive integer；
Target data is obtained according to the radar echo signal returned after each scanning, the target data includes at least atural object
Clutter map CF and absolute phase POTS, the CF are for identifying whether each range bin of the radar echo signal is that atural object is miscellaneous
Wave, the POTS are used to record the absolute phase of each range bin；
When it is ground clutter that the CF of range bin to be selected, which identifies the range bin to be selected, the of the range bin to be selected is calculated
The correlation of one absolute phase sequence and the second absolute phase sequence, the first absolute phase sequence are in the M period
The set for the POTS that the high PRF scanning obtains is carried out, the second absolute phase sequence is described in the M period interior progress
The set for the PORTS that low PRF scanning obtains；
When the correlation is greater than the first preset threshold value, it is determined that the range bin to be selected is target range library.
In a kind of possible design, in the first implementation of the embodiment of the present application first aspect, the calculating
First absolute phase sequence of the range bin to be selected and the correlation of the second absolute phase sequence include:
The correlation is calculated as follows:
C=corr (CS_POTS_TAB, CD_POTS_TAB)；
Wherein, the cov indicates covariance, and the σ indicates standard deviation, and the CS_POTS_TAB indicates that described first is exhausted
To phase sequence, the CD_POTS_TAB indicates the second absolute phase sequence.
In a kind of possible design, in second of implementation of the embodiment of the present application first aspect, the basis
The radar echo signal returned after scanning every time obtains target data:
Ground clutter probability is obtained according to the radar echo signal that each scanning returns, with the corresponding range bin of determination
It whether is ground clutter；
POTS is calculated as follows according to the radar echo signal that each scanning returns:
Wherein, the N is hits, and the x (i) indicates i-th of pulse of the horizontal received multiple I/Q signal of horizontal emission,
The b (i) indicates the phase of the BURST signal of i-th of pulse, and arg expression takes phase to plural number.
In a kind of possible design, in the third implementation of the embodiment of the present application first aspect, according to described
Whether the radar echo signal that scanning returns every time obtains ground clutter probability, be ground clutter packet with the corresponding range bin of determination
It includes:
When the signal-to-noise ratio of the radar echo signal is less than the second preset threshold value, the atural object is calculated as follows
Clutter probability:
CP=(max (MFTDBZ, MFSPIN)*1+MFCPA*1.01)/(1+1.01)；
The CP indicates the ground clutter probability, the texture structure TDBZ of the reflectivity factor, the reflectivity because
The sub variation SPIN radially and clutter phased array CPA is the characteristic quantity of the radar echo signal, described
MFTDBZ、MFSPINAnd MFCPAThe value respectively obtained after being converted corresponding characteristic quantity for the corresponding membership function of use, it is described
MFTDBZ、MFSPINAnd MFCPAValue range be [0,1]；
The range bin that ground clutter probability is greater than third preset threshold value is identified as ground clutter.
In a kind of possible design, in the 4th kind of implementation of the embodiment of the present application first aspect, described first
Preset threshold value is 0.3.
The second aspect of the embodiment of the present application provides a kind of screening plant of ground clutter, comprising:
Acquisition unit carries out one group of plan-position to the identical elevation angle within each period and sweeps for acquiring M period
It retouches, one group of scanning includes a high pulse repetition frequency PRF scanning and primary low PRF scanning, and the M is positive integer；
Acquiring unit, for obtaining target data, the number of targets according to the radar echo signal returned after scanning every time
It is used to identify each distance of the radar echo signal according to including at least ground clutter figure CF and absolute phase POTS, the CF
Whether library is ground clutter, and the POTS is used to record the absolute phase of each range bin；
Computing unit, it is described for calculating when it is ground clutter that the CF of range bin to be selected, which identifies the range bin to be selected,
First absolute phase sequence of range bin to be selected and the correlation of the second absolute phase sequence, the first absolute phase sequence are
The set for the PORTS that the high PRF scanning obtains is carried out in the M period, the second absolute phase sequence is the M
The set for the PORTS that the low PRF scanning obtains is carried out in period；
Determination unit, when the correlation is greater than the first preset threshold value, then for determining that the range bin to be selected is mesh
Subject distance library.
In a kind of possible design, in the first implementation of the embodiment of the present application second aspect, the calculating
Unit includes:
The correlation is calculated as follows:
C=corr (CS_POTS_TAB, CD_POTS_TAB)；
Wherein, the cov indicates covariance, and the σ indicates standard deviation, and the CS_POTS_TAB indicates that described first is exhausted
To phase sequence, the CD_POTS_TAB indicates the second absolute phase sequence.
In a kind of possible design, in second of implementation of the embodiment of the present application second aspect, the acquisition
Unit includes:
Module is obtained, the radar echo signal for returning according to each scanning obtains ground clutter probability, with true
Whether fixed corresponding range bin is ground clutter；
Computing module, the radar echo signal for being returned according to each scanning calculate POTS as follows:
Wherein, the N is hits, and the x (i) indicates i-th of pulse of the horizontal received multiple I/Q signal of horizontal emission,
The b (i) indicates the phase of the BURST signal of i-th of pulse, and arg expression takes phase to plural number.
In a kind of possible design, in the third implementation of the embodiment of the present application second aspect, the acquisition
Module is specifically used for:
When the signal-to-noise ratio of the radar echo signal is less than the second preset threshold value, the atural object is calculated as follows
Clutter probability:
CP=(max (MFTDBZ, MFSPIN)*1+MFCPA*1.01)/(1+1.01)
The CP indicates the ground clutter probability, and the TDBZ, the SPIN and the CPA are the radar return
The characteristic quantity of signal, the MFTDBZ、MFSPINAnd MFCPAAfter corresponding characteristic quantity is converted using corresponding membership function
The value respectively obtained, the MFTDBZ、MFSPINAnd MFCPAValue range be [0,1]；
The range bin that ground clutter probability is greater than third preset threshold value is identified as ground clutter.
In a kind of possible design, in the 4th kind of implementation of the embodiment of the present application second aspect, described first
Preset threshold value is 0.3.
As can be seen from the above technical solutions, the embodiment of the present invention has the advantage that by with ground clutter mark
In the range bin of knowledge, the first absolute phase sequence of each range bin and the correlation of the second absolute phase sequence are calculated, with sieve
Suitable range bin is selected to calculate air index.
Detailed description of the invention
Fig. 1 is a kind of embodiment schematic diagram of possible data screening method provided in an embodiment of the present invention；
Fig. 2 is a kind of embodiment schematic diagram of the screening plant of possible ground clutter provided in an embodiment of the present invention；
Fig. 3 is the embodiment schematic diagram of the screening plant of alternatively possible ground clutter provided in an embodiment of the present invention；
Fig. 4 is the embodiment schematic diagram of the screening plant of alternatively possible ground clutter provided in an embodiment of the present invention.
Specific embodiment
The embodiment of the invention provides a kind of data screening method and the screening plants of ground clutter, by radar based on
Suitable ground clutter is selected when calculating air index.
Below in conjunction with the attached drawing in the embodiment of the present application, technical solutions in the embodiments of the present application carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of embodiments of the present application, instead of all the embodiments.It is based on
Embodiment in the application, those skilled in the art's every other implementation obtained without making creative work
Example, shall fall in the protection scope of this application.
Generally there are two types of scan patterns, including plan-position scanning and body to sweep for weather radar system.Wherein, plan-position is swept
It retouches and refers to radar fixed elevation, carry out a 360 degree of scannings；Body sweeps the plan-position for referring to one group of different elevation angle from low to high
Scanning.In weather radar signal processing, the small distance unit that radar echo signal is divided into along directions of rays by distance claims
For range bin (range bin).When weather radar calculates air index, absolutely thinking for the background return of radar can use
The air index variation near radar is calculated for variation, however, not all ground clutter is suitable for calculating refractive index,
In view of this, the present invention provides a kind of data screening methods, for picking out suitable range bin from ground clutter, below
Will in conjunction with specific embodiments, it elaborates to data screening method of the invention, specifically referring to Figure 1, the present invention provides
Data screening method include:
101, M period is acquired, one group of plan-position scanning is carried out to the identical elevation angle within each period；
Within each period, continuous different pulse repetition (pulse recurrence is carried out to the identical elevation angle
Frequency, PRF) plan-position twice scanning, for the first time scanning a low PRF, such as 322 can be used, for convenient for table
It states, CS scanning can be referred to as in the present invention；A high PRF, such as 1014 can be used in second of scanning, for convenient for statement, sheet
It can be referred to as CF scanning in invention, which is considered as one group of plan-position scanning.In practical application, M can be acquired
Period, to obtain M group scan data.For example, every the one group of scanning of progress in 6 minutes, continued for some time in clear sky, such as
One hour acquires 10 periods, obtained 10 groups of scan datas.In the present invention, M is positive integer, this hair of specific value
It is bright without limitation.
It should be noted that the different periods, the elevation angle for carrying out plan-position scanning is identical.
102, ground clutter probability is obtained according to the radar echo signal that scanning returns every time；
It carries out one group of plan-position in a cycle to scan, one group of plan-position scanning includes primary low PRF scanning and one
Secondary high PRF scanning can get target data according to the radar echo signal returned after each scanning, wherein the target data is extremely
Less include ground clutter figure CF and absolute phase POTS, wherein CF be used for identify radar echo signal each range bin whether
For ground clutter, it can be understood as each range bin of radar echo signal has a CF mark, which identifies for indicating
Whether corresponding range bin is ground clutter；POTS is used to record the absolute phase of each range bin echo.It should be noted that
In practical application, judge range bin whether be ground clutter mode there are many, in the embodiment of the present invention can according to radar return
Wave signal obtains the ground clutter probability of each range bin, can specifically include:
Whether step 1, the signal-to-noise ratio for checking radar echo signal, judge signal-to-noise ratio less than the second preset threshold value；If so,
Execute step 2；
The small distance unit being divided into along directions of rays by distance in radar echo signal processing, referred to as range bin.It checks
The signal-to-noise ratio of radar echo signal corresponding to each range bin, if signal-to-noise ratio is less than the second preset threshold value, the range bin
It is considered as noise, is not processed；If signal-to-noise ratio is greater than the second preset threshold value, subsequent processing is just carried out.
In practical application, which can be 3dB.
Step 2, the ground clutter probability for calculating each range bin as follows；
CP=(max (MFTDBZ, MFSPIN)*1+MFCPA*1.01)/(1+1.01)；
Wherein, the CP indicates ground clutter probability, the texture structure TDBZ of the reflectivity factor, reflectivity factor edge
Radial variation SPIN and clutter phased array (clutter phase alignment, CPA) is the radar echo signal
Characteristic quantity, the MFTDBZ、MFSPINAnd MFCPAAfter corresponding characteristic quantity is converted for the corresponding membership function of use respectively
Obtained value, the MFTDBZ、MFSPINAnd MFCPAValue range be [0,1]；
It should be noted that TDBZ indicates that the mean value of square of reflectivity factor between adjacencies database, SPIN indicate adjacent
The difference of reflectivity factor is greater than the points of certain threshold value between range bin, is identified with the percentage of maximum possible number, and CPA is used for
Judge the constancy of exomonental absolute return phase being made of a radar wave speed data.For changeless mesh
Mark, CPA is 1, if target imperfect stability in time of measuring will lead to CPA less than 1.The POTS of range bin is more stable, more
It is possible that including clutter.Three characteristic quantities are calculated according to IQ clock signal, and IQ clock signal refers to radar return
90 degree of obtained signals of signal phase shift, wherein I indicates that signal sequence in the same direction, Q indicate quadrature signal sequence.Three characteristic quantities
Calculation formula can be such that
xiFor multiple IQ clock signal, N is the length of multiple IQ clock signal；
M indicates the number of range bin included in one-dimensional calculating core；
Thresh is preset constant, such as 4dbz.
After three characteristic quantities are calculated, by respective membership function respectively each characteristic quantity be converted into 0-1 it
Between value, specific conversion process is as follows:
Therefore in the above manner, the ground clutter probability of each range bin can be calculated.
When ground clutter probability is greater than third preset threshold value, then identifying corresponding range bin is ground clutter, for example, should
Third threshold can be 0.5, if the ground clutter probability of some range bin is greater than 0.5, then it is assumed that the range bin is atural object
Clutter.
It should be noted that in practical application, mark range bin whether be ground clutter algorithm there are many, for example, also
CMD algorithm etc. can be used, the specific present invention is without limitation.
103, POTS is calculated according to the radar echo signal that scanning returns every time；
In the embodiment of the present invention, the POTS of each range bin in radar echo signal can be calculated according to following formula:
Wherein, the N is hits, and the x (i) indicates i-th of pulse of the horizontal received multiple I/Q signal of horizontal emission,
The b (i) indicates the phase of the BURST signal of i-th of pulse, and arg expression takes phase to plural number.
It should be noted that being known in the present embodiment by the ground clutter probability that step 102 calculates each range bin
Not Chu the range bin whether be ground clutter, the POTS of each range bin is calculated by step 103, and between the two processes simultaneously
There is no the sequencings of step, can first carry out step 102, can also first carry out step 103, or be performed simultaneously, specifically
Herein without limitation.
104, the first absolute phase sequence of range bin to be selected and the correlation of the second absolute phase sequence are calculated；
After determining whether range bin is ground clutter and calculates the POTS of each range bin according to CS, acquisition plane position
All radial directions in scanning pick out each range bin and exist to the range bin for having ground clutter to identify all in each radial direction
The POTS that the twice sweep in each period obtains, obtains the first absolute phase sequence and the second absolute phase sequence respectively, wherein the
One absolute phase sequence is that the set for the POTS that high PRF scanning (can be CD scanning in the present embodiment) obtains is carried out in M period,
Second absolute phase sequence is the collection that the PORTS that low PRF scanning (can be CS scanning in the present embodiment) obtains is carried out in M period
It closes, therefore the first absolute phase sequence C D_POTS_TAB=[pots_cd1, pots_cd2 ..., pots_cdM], the second absolute phase
Bit sequence CS_POTS_TAB=[pots_cs1, pots_cs2 ..., pots_csM], and first is calculated as follows absolutely
To the correlation of phase sequence and the second absolute phase sequence:
C=corr (CS_POTS_TAB, CD_POTS_TAB)；
Wherein, the cov indicates covariance, the σ indicates standard deviation, due to the calculation formula and standard deviation of covariance
Calculation formula is the prior art, and specific details are not described herein again.
105, determine that range bin to be selected is target range library.
After the correlation of the first absolute phase sequence and the second absolute phase sequence that calculate range bin to be selected, if should
Correlation is greater than the first preset threshold value, it is determined that the range bin to be selected is target range library, that is, thinks that the range bin is suitble to calculate
Air index.In practical application, the value of first preset threshold value can be 0.3 or other numerical value, not do herein specifically
It limits.
In the embodiment of the present invention, by calculating the first of each range bin in the range bin identified with ground clutter
The correlation of absolute phase sequence and the second absolute phase sequence calculates air index to filter out suitable range bin,
Reduce hardware cost when selecting suitable ground clutter.
Data screening method in the embodiment of the present invention is described above, below to the atural object in the embodiment of the present invention
The screening plant of clutter is described, referring to Fig. 2, the screening plant of the ground clutter in the embodiment of the present invention includes:
Acquisition unit 201 carries out one group of plan-position to the identical elevation angle within each period for acquiring M period
Scanning, one group of scanning include a high pulse repetition frequency PRF scanning and primary low PRF scanning, and the M is positive integer；
Acquiring unit 202, for obtaining target data, the target according to the radar echo signal returned after scanning every time
Data include at least ground clutter figure CF and absolute phase POTS, the CF for identify each of described radar echo signal away from
It whether is ground clutter from library, the POTS is used to record the absolute phase of each range bin；
Computing unit 203, when it is ground clutter that the CF of range bin to be selected, which identifies the range bin to be selected, for calculating
State the first absolute phase sequence of range bin to be selected and the correlation of the second absolute phase sequence, the first absolute phase sequence
For the set for carrying out the PORTS that the high PRF scanning obtains in the M period, the second absolute phase sequence is the M
The set for the PORTS that the low PRF scanning obtains is carried out in a period；
Determination unit 204, when the correlation is greater than the first preset threshold value, then for determining that the range bin to be selected is
Target range library.
For ease of understanding, the screening plant of the ground clutter in the embodiment of the present invention is described in detail below, upper
State it is shown in Fig. 2 on the basis of, referring to Fig. 3, Fig. 3 be the embodiment of the present invention in ground clutter screening plant another implementation
It illustrates and is intended to, optionally, computing unit 303 is specifically used for:
The correlation is calculated as follows:
C=corr (CS_POTS_TAB, CD_POTS_TAB)；
Wherein, the cov indicates covariance, and the σ indicates standard deviation, and the CS_POTS_TAB indicates that described first is exhausted
To phase sequence, the CD_POTS_TAB indicates the second absolute phase sequence.
Optionally, acquiring unit 302 includes:
Module 3021 is obtained, the radar echo signal for being returned according to each scanning obtains ground clutter probability,
It whether is ground clutter with the corresponding range bin of determination；
Computing module 3022, the radar echo signal for being returned according to each scanning calculate as follows
POTS:
Wherein, the N is hits, and the x (i) indicates i-th of pulse of the horizontal received multiple I/Q signal of horizontal emission,
The b (i) indicates the phase of the BURST signal of i-th of pulse, and arg expression takes phase to plural number.
Optionally, module 3021 is obtained to be specifically used for:
When the signal-to-noise ratio of the radar echo signal is less than the second preset threshold value, the atural object is calculated as follows
Clutter probability:
CP=(max (MFTDBZ, MFSPIN)*1+MFCPA*1.01)/(1+1.01)；
The CP indicates the ground clutter probability, and the TDBZ, the SPIN and the CPA are the radar return
The characteristic quantity of signal, described MFTDBZ, MFSPINAnd MFCPAFor corresponding characteristic quantity is converted using corresponding membership function
The value respectively obtained afterwards, the MFTDBZ、MFSPINAnd MFCPAValue range be [0,1]；
The range bin that ground clutter probability is greater than third preset threshold value is identified as ground clutter.
Angle of the above figure 2 to Fig. 3 from modular functionality entity is respectively to the sieve of the ground clutter in the embodiment of the present invention
Screening device is described, and carries out below from screening plant of the angle of hardware handles to the ground clutter in the embodiment of the present invention
Detailed description, referring to Fig. 4,400 one embodiment of screening plant of the ground clutter in the embodiment of the present invention, comprising:
(wherein the quantity of processor 403 can be with for input unit 401, output device 402, processor 403 and memory 404
One or more, in Fig. 4 by taking a processor 403 as an example).In some embodiments of the invention, input unit 401, output
Device 402, processor 403 and memory 404 can be connected by bus or other means, wherein to be connected by bus in Fig. 4
For.
Wherein, the operational order stored by calling memory 404, processor 403, for executing following steps:
M period is acquired, one group of plan-position scanning, one group of scanning are carried out to the identical elevation angle within each period
Including a high pulse repetition frequency PRF scanning and primary low PRF scanning, the M is positive integer；
Target data is obtained according to the radar echo signal returned after each scanning, the target data includes at least atural object
Clutter map CF and absolute phase POTS, the CF are for identifying whether each range bin of the radar echo signal is that atural object is miscellaneous
Wave, the POTS are used to record the absolute phase of each range bin；
When it is ground clutter that the CF of range bin to be selected, which identifies the range bin to be selected, the of the range bin to be selected is calculated
The correlation of one absolute phase sequence and the second absolute phase sequence, the first absolute phase sequence are in the M period
The set for the POTS that the high PRF scanning obtains is carried out, the second absolute phase sequence is described in the M period interior progress
The set for the PORTS that low PRF scanning obtains；
When the correlation is greater than the first preset threshold value, it is determined that the range bin to be selected is target range library.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description,
The specific work process of device and unit, can refer to corresponding processes in the foregoing method embodiment, and details are not described herein.
In several embodiments provided herein, it should be understood that disclosed system, device and method can be with
It realizes by another way.For example, the apparatus embodiments described above are merely exemplary, for example, the unit
It divides, only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple units or components
It can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point, it is shown or
The mutual coupling, direct-coupling or communication connection discussed can be through some interfaces, the indirect coupling of device or unit
It closes or communicates to connect, can be electrical property, mechanical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple
In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme
's.
It, can also be in addition, the functional units in various embodiments of the present invention may be integrated into one processing unit
It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list
Member both can take the form of hardware realization, can also realize in the form of software functional units.
If the integrated unit is realized in the form of SFU software functional unit and sells or use as independent product
When, it can store in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially
The all or part of the part that contributes to existing technology or the technical solution can be in the form of software products in other words
It embodies, which is stored in a storage medium, including some instructions are used so that a computer
Equipment (can be personal computer, server or the network equipment etc.) executes the complete of each embodiment the method for the present invention
Portion or part steps.And storage medium above-mentioned include: USB flash disk, mobile hard disk, read-only memory (Read-Only Memory,
ROM), random access memory (Random Access Memory, RAM), magnetic or disk etc. are various can store program
The medium of code.
The above, the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although referring to before
Stating embodiment, invention is explained in detail, those skilled in the art should understand that: it still can be to preceding
Technical solution documented by each embodiment is stated to modify or equivalent replacement of some of the technical features；And these
It modifies or replaces, the spirit and scope for technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution.
Claims (6)
1. a kind of data screening method characterized by comprising
M period is acquired, one group of plan-position scanning is carried out to the identical elevation angle within each period, one group of scanning includes
High pulse repetition frequency PRF scanning and primary low PRF scanning, the M is positive integer；
Target data is obtained according to the radar echo signal returned after each scanning, the target data includes at least ground clutter
Scheme CF and absolute phase POTS, the CF for identifying whether each range bin of the radar echo signal is ground clutter,
The POTS is used to record the absolute phase of each range bin；
When it is ground clutter that the CF of range bin to be selected, which identifies the range bin to be selected, calculate the range bin to be selected first is exhausted
To the correlation of phase sequence and the second absolute phase sequence, the first absolute phase sequence is to carry out in the M period
The set for the POTS that the high PRF scanning obtains, the second absolute phase sequence are that progress is described low in the M period
The set for the PORTS that PRF scanning obtains；
When the correlation is greater than the first preset threshold value, it is determined that the range bin to be selected is target range library；
First absolute phase sequence of the calculating range bin to be selected and the correlation of the second absolute phase sequence include:
The correlation is calculated as follows:
C=corr (CS_POTS_TAB, CD_POTS_TAB)；
Wherein, the cov indicates covariance, and the σ indicates standard deviation, and the CS_POTS_TAB indicates the first absolute phase
Bit sequence, the CD_POTS_TAB indicate the second absolute phase sequence；
The radar echo signal that the basis returns after scanning every time obtains target data:
According to it is described it is each scanning return radar echo signal obtain ground clutter probability, with the corresponding range bin of determination whether
For ground clutter；
POTS is calculated as follows according to the radar echo signal that each scanning returns:
Wherein, the N is hits, and the x (i) indicates i-th of pulse of the horizontal received multiple I/Q signal of horizontal emission, described
B (i) indicates the phase of the BURST signal of i-th of pulse, and arg expression takes phase to plural number.
2. the method according to claim 1, wherein being obtained according to the radar echo signal that each scanning returns
Ground clutter probability is obtained, whether is that ground clutter includes: with the corresponding range bin of determination
When the signal-to-noise ratio of the radar echo signal is less than the second preset threshold value, the ground clutter is calculated as follows
Probability:
CP=(max (MFTDBZ,MFSPIN)*1+MFCPA*1.01)/(1+1.01)；
The CP indicates the ground clutter probability, and the texture structure TDBZ of reflectivity factor, the reflectivity factor are radially
Variation SPIN and the clutter phased array CPA be the radar echo signal characteristic quantity, the MFTDBZ、MFSPINWith
MFCPAThe value respectively obtained after being converted corresponding characteristic quantity for the corresponding membership function of use, the MFTDBZ、MFSPINWith
MFCPAValue range be [0,1]；
The range bin that ground clutter probability is greater than third preset threshold value is identified as ground clutter.
3. method according to any one of claim 1 to 2, which is characterized in that first preset threshold value is 0.3.
4. a kind of screening plant of ground clutter characterized by comprising
Acquisition unit carries out one group of plan-position scanning, institute to the identical elevation angle within each period for acquiring M period
Stating one group of scanning includes a high pulse repetition frequency PRF scanning and primary low PRF scanning, and the M is positive integer；
Acquiring unit, for obtaining target data according to the radar echo signal returned after scanning every time, the target data is extremely
It less include ground clutter figure CF and absolute phase POTS, each range bin that the CF is used to identify the radar echo signal is
No is ground clutter, and the POTS is used to record the absolute phase of each range bin；
Computing unit, it is described to be selected for calculating when it is ground clutter that the CF of range bin to be selected, which identifies the range bin to be selected,
First absolute phase sequence of range bin and the correlation of the second absolute phase sequence, the first absolute phase sequence is described
The set for the PORTS that the high PRF scanning obtains is carried out in M period, the second absolute phase sequence is the M period
The interior set for carrying out the PORTS that the low PRF scanning obtains；
Determination unit, when the correlation be greater than the first preset threshold value when, then for determine the range bin to be selected be target away from
From library；
The computing unit is specifically used for:
The correlation is calculated as follows:
C=corr (CS_POTS_TAB, CD_POTS_TAB)；
Wherein, the cov indicates covariance, and the σ indicates standard deviation, and the CS_POTS_TAB indicates the first absolute phase
Bit sequence, the CD_POTS_TAB indicate the second absolute phase sequence；
The acquiring unit includes:
Module is obtained, the radar echo signal for returning according to each scanning obtains ground clutter probability, with determining pair
Whether the range bin answered is ground clutter；
Computing module, the radar echo signal for being returned according to each scanning calculate POTS as follows:
Wherein, the N is hits, and the x (i) indicates i-th of pulse of the horizontal received multiple I/Q signal of horizontal emission, described
B (i) indicates the phase of the BURST signal of i-th of pulse, and arg expression takes phase to plural number.
5. device according to claim 4, which is characterized in that the acquisition module is specifically used for:
When the signal-to-noise ratio of the radar echo signal is less than the second preset threshold value, the ground clutter is calculated as follows
Probability:
CP=(max (MFTDBZ,MFSPIN)*1+MFCPA*1.01)/(1+1.01)；
The CP indicates the ground clutter probability, and the texture structure TDBZ of reflectivity factor, the reflectivity factor are radially
Variation SPIN and the clutter phased array CPA be the radar echo signal characteristic quantity, the MFTDBZ、MFSPINWith
MFCPAThe value respectively obtained after being converted corresponding characteristic quantity for the corresponding membership function of use, the MFTDBZ、MFSPINWith
MFCPAValue range be [0,1]；
The range bin that ground clutter probability is greater than third preset threshold value is identified as ground clutter.
6. the device according to any one of claim 4 to 5, which is characterized in that first preset threshold value is 0.3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710827028.6A CN107607955B (en) | 2017-09-13 | 2017-09-13 | A kind of screening plant of data screening method and ground clutter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710827028.6A CN107607955B (en) | 2017-09-13 | 2017-09-13 | A kind of screening plant of data screening method and ground clutter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107607955A CN107607955A (en) | 2018-01-19 |
| CN107607955B true CN107607955B (en) | 2019-11-15 |
Family
ID=61063907
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710827028.6A Active CN107607955B (en) | 2017-09-13 | 2017-09-13 | A kind of screening plant of data screening method and ground clutter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107607955B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110308428B (en) * | 2018-03-27 | 2021-05-14 | 中国民航科学技术研究院 | Method for establishing low-altitude clutter distribution simulation model |
| CN113835077B (en) * | 2021-11-23 | 2022-02-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Search radar target detection method and system based on variable pulse repetition frequency |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102879769A (en) * | 2012-10-22 | 2013-01-16 | 北京敏视达雷达有限公司 | Ground clutter recognition method |
| CN104345305A (en) * | 2013-07-30 | 2015-02-11 | 北京敏视达雷达有限公司 | Ground clutter residual filtering method |
| CN104535980A (en) * | 2014-12-24 | 2015-04-22 | 西安电子工程研究所 | Five-dimensional dynamic stereoscopic clutter map implementation method based on DDR3-SDRAM |
| CN105093359A (en) * | 2015-07-27 | 2015-11-25 | 安徽四创电子股份有限公司 | Quality testing method for intensity data of weather radar |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6437738B1 (en) * | 2001-02-12 | 2002-08-20 | Us Commerce | Hexagonal-annulus phased array antenna for radar wind profiling on moving platforms |
-
2017
- 2017-09-13 CN CN201710827028.6A patent/CN107607955B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102879769A (en) * | 2012-10-22 | 2013-01-16 | 北京敏视达雷达有限公司 | Ground clutter recognition method |
| CN104345305A (en) * | 2013-07-30 | 2015-02-11 | 北京敏视达雷达有限公司 | Ground clutter residual filtering method |
| CN104535980A (en) * | 2014-12-24 | 2015-04-22 | 西安电子工程研究所 | Five-dimensional dynamic stereoscopic clutter map implementation method based on DDR3-SDRAM |
| CN105093359A (en) * | 2015-07-27 | 2015-11-25 | 安徽四创电子股份有限公司 | Quality testing method for intensity data of weather radar |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107607955A (en) | 2018-01-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105488770B (en) | A kind of airborne laser radar point cloud filtering method of object-oriented | |
| CN107942342B (en) | Data processing method, device, system and the storage medium of anemometry laser radar | |
| CN108802722B (en) | It is a kind of based on tracking before the Faint target detection virtually composed | |
| CN109407071A (en) | Radar range finding method, radar range unit, unmanned plane and storage medium | |
| CN105425231B (en) | A kind of multiple-sensor and multiple-object localization method based on layering projection and Taylor expansion | |
| CN104714225B (en) | Dynamic programming tracking-before-detection method based on generalized likelihood ratios | |
| CN107121670A (en) | A kind of anti-unmanned plane system of defense based on synthetic aperture radar | |
| CN107561509B (en) | Airborne millimeter wave radar power line detection method | |
| CN111461037B (en) | End-to-end gesture recognition method based on FMCW radar | |
| CN107015205A (en) | A kind of false target removing method of distributed MIMO detections of radar | |
| CN107229084A (en) | A kind of automatic identification, tracks and predicts contracurrent system mesh calibration method | |
| CN107607955B (en) | A kind of screening plant of data screening method and ground clutter | |
| CN107255814A (en) | A kind of radar target detection method based on LFMSK waveforms | |
| CN107423734A (en) | A kind of SAR image targets in ocean quick determination method and device | |
| CN110988863A (en) | Novel millimeter wave radar gesture signal processing method | |
| CN106872967A (en) | A kind of moving target detecting system and method based on bistatic radar | |
| CN110007299A (en) | A kind of dim target detection tracking based on hybrid coordinate puppet spectral technology | |
| JPH0980146A (en) | Radar apparatus | |
| CN109471097A (en) | A kind of through-wall radar Signal optimum processing method and device | |
| CN106291500B (en) | System on chip and its object detection method based on FPGA | |
| CN110488273A (en) | A kind of vehicle tracking detection method and device based on radar | |
| CN108845318A (en) | Spaceborne high score wide cut imaging method based on Relax algorithm | |
| CN105911546B (en) | A kind of recognition methods of sea clutter and device | |
| CN109581350B (en) | Radar distance and speed measurement method and device based on time-frequency integral interpolation | |
| CN106249241B (en) | A kind of self-adapting clutter power statistic algorithm |
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 |