CN109358331A - The real-time dynamic noise power detecting method of weather radar - Google Patents
The real-time dynamic noise power detecting method of weather radar Download PDFInfo
- Publication number
- CN109358331A CN109358331A CN201811198843.1A CN201811198843A CN109358331A CN 109358331 A CN109358331 A CN 109358331A CN 201811198843 A CN201811198843 A CN 201811198843A CN 109358331 A CN109358331 A CN 109358331A
- Authority
- CN
- China
- Prior art keywords
- power
- noise
- echo
- signal
- short
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000001514 detection method Methods 0.000 claims abstract description 22
- 230000010287 polarization Effects 0.000 claims description 38
- 230000001427 coherent effect Effects 0.000 claims description 11
- 238000013139 quantization Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000009977 dual effect Effects 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000005311 autocorrelation function Methods 0.000 claims description 2
- 230000001788 irregular Effects 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 230000005855 radiation Effects 0.000 abstract description 7
- 238000001228 spectrum Methods 0.000 abstract description 5
- 238000002310 reflectometry Methods 0.000 abstract description 3
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 abstract description 2
- 238000012545 processing Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000001556 precipitation Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000011897 real-time detection Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention provides a kind of real-time dynamic noise power detecting method of weather radar, for improving the signal-to-noise ratio of radar return data, increases observation areal extent;Meanwhile improving GPR Detection Data quality, reflectivity factor estimated bias is reduced, the estimation of speed spectrum width is improved, improves the validity of polarimetric radar related coefficient and the estimation of reflectance difference rate.The case where noise bounce caused by the method is adapted to pulse radar due to component parameter variation, the rotation of signal path characteristic drift, antenna, the sun and cosmic radiation etc., compared with using quiescent noise level setting method, this method dynamically improves the signal-to-noise ratio of echo corresponding to each transmitting pulse, only solar radiation noise, to the weather radar of X and following wave band, can adapt in real time be more than 20dB noise variation.
Description
Technical field
The invention belongs to electronic information field more particularly to a kind of real-time dynamic noise power detecting methods of weather radar.
Background technique
Weather radar measures target position and speed, and according to thunder for finding meteorological target, i.e. atmosphere and its precipitation object
Reflectivity factor, average radial velocity, normal-moveout spectrum up to the characteristic retrieval target of system and echo in phase, frequency, amplitude
The fundamental physical quantity measured directly such as width, reflectance difference rate, difference travel phase shift rate and correlation, so as to further retrieving precipitation class
The meteorology parameter such as type, particle spectra and precipitation rate.
As Figure 2-3, due to target range and its radar cross section, echo power is not very in many cases
Height, noise is relatively low, and at this moment, the power level of noise will cause large effect to the estimation of above-mentioned fundamental physical quantity, especially
Reflectivity factor Z h, speed spectrum width σv, it is related | ρhv(0) | and the estimated bias of the parameters such as reflectance difference rate Zdr increases, and influences
The quality of data of target measurement.
Doing for noise power thresholds is periodically arranged in the periodic measurement noise of radar receiver that some radar producers use at present
The method that a noise-power measurement is carried out before each radar observation scanning that method or meteorological department require, is all to emit
Machine does not work, and measures under the premise of no echo.This method has certain effect, but is directed to strong and weak mesh for weather radar
Mapping amount will reach the requirement of certain quality of data, substantially be difficult to meet.
Main cause is in terms of the generation of radar noise includes three first, is made an uproar caused by operation of receiver temperature change
Sound rises and falls;Second, noise bounce caused by feeder-path, especially rotary joint;Third, antenna receive from the sun
And universe wide spectrum radiation at any time with the fluctuating of aerial position.The variation of these noises is unpredictable, and current Noise Method
Cannot accurately and timely it estimate again.
If using the method for statistics, the noise threshold for using the noise power of maximum case to work as radar be will lead to
Radar sensitivity degradation.
Therefore, need it is a kind of accurately, in real time, the method or technique of dynamic detection weather radar noise power.
Summary of the invention
It is an object of the invention to solve the problems of the above-mentioned prior art, provide weather radar real-time dynamic noise function
Rate detection method, for weather radar carry out accurately, in real time, dynamically noise power detect, to automatically adapt to
Radar noise power changes, and can improve radar sensitivity to the maximum extent, improves radar surveying parameter.
The present invention adopts the following technical scheme:
The real-time dynamic noise power detecting method of weather radar characterized by comprising
The transmitting pulse that step 1. detects Radar Digital Receiver output with continuous weak power detection module one by one is corresponding
Range sweep echo-signal, and carry out noise power estimation;
Continuous weak power detection module is carried out according to distance/time sequencing according to a certain distance range/length of time
Baseband signal power is detected in continuous low signal areas, obtains comprising the continuous weak power echo time section including noise power;
Step 2. short-time rating standard difference module, in short-term phase standard difference module, short-time average normalize coherent power
Module and 04 detection modules of rank related coefficient module are according to relatively short certain period of time in short-term, to the company in step 1
Continue weak power echo time segment data and carries out that short-time rating standard deviation, phase standard is poor in short-term, short-time average normalizes coherent function
Rate and the in short-term parameter Estimation of 0 rank related coefficient;
Step 3. carries out blurring and reasoning to the parameter in step 2 with comprehensive judgment module, determines determining for pure noise range
Position calculates the power level of noise.
The further technical solution of the present invention is that the continuous weak signal power detection of the step 1 includes:
1) the echo power P of each range gate Gi, is calculatedHiAnd PVi;
A. each distance can directly be calculated using the output H (I, Q) of its linear digital receiver if Doppler radar
Echo power P on door GiHi, it is shown below:
H (I, Q) is the digital signal of weather radar horizontal polarization echo, in-phase component HI, the Gi range gate
Numerical value is HIi, quadrature component HQ, the numerical value of the Gi range gate is HQi;
B. if dual polarization Doppler radar, the number using the horizontal polarization echo of its linear digital receiver output is believed
Number H (I, Q), can directly calculate horizontal polarization echo power P on each range gate GiHi, it is shown below:
Using its vertical polarization echo digital signal V (I, Q), vertical polarization echo on each range gate Gi can be directly calculated
Power PVi, it is shown below:
Wherein, vertical polarization echo in-phase component is VI, the numerical value of the Gi range gate is VIi;Its quadrature component is VQ,
The numerical value of the Gi range gate is VQi。
2), confirm weak echo signal in simultaneously tag distance scanning
According to the running parameter of radar, the noise level of radar is determined, mainly include two parts, i.e. receiver noise and day
Line noise, wherein receiver noise estimates receiver noise power level by receiver theoretical calculation as the following formula:
PRNoise=-114dB+10log10(Bn[MHz])+10·log10(F0)(dBm)------------(3)
Wherein, BnFor receiver bandwidth, unit MHz, F0For the noise coefficient of radar system, using under room temperature
Measured value calculates;
Antenna noise then estimates antenna by sun ephemeris and its spectral density according to the height of radar station and longitude and latitude
Maximum solar power P under the conditions of relevant work frequency range, bandwidth and antenna gainSolar。
3, weak echo signal section in range sweep is found out
Compared by power, finds out the section of the relatively small continuous range gate of signal power in range sweep surface sweeping, press
Following formula respectively obtains the weak echo section in horizontal receiving channel and vertical reception channel;
M is continuous range gate number.
Currently preferred technical solution is that M number of continuous range gate, value is in 50-100 point.
The further technical solution of the present invention is, phase standard difference parameter, short-time rating standard in short-term in the step 2
Parameter, short-time average normalize coherent power parameter, in short-term 0 rank related coefficient method for parameter estimation are as follows:
Phase standard difference parameter Estimation in short-term is detected using following formula and is calculated:
Wherein, N is continuous distance points, xphiFor the phase of the Gi range gate echo-signal;μphIt is that this is continuous apart from upper
N point phase average, i.e.,
The estimation of short-time rating standard deviation criteria, is detected using following formula and is calculated:
Wherein, xPiFor the power of Gi range gate echo-signal;μPIt is that this is continuous apart from upper N point power average value, i.e.,
Short-time average normalizes the estimation of coherent power parameter, and single repetition rate is calculated using following formula detection:
It is related for 1 rank of the Gi range gate echo-signal;For the power of the Gi range gate echo-signal,
And
Or irregular repetition rate, it is detected and is calculated using following formula:
For the 1st pulse spacing,It is related for 1 rank of the Gi range gate echo-signal;For the Gi distance
The power of door echo-signal, andFor the 2nd pulse spacing,It is the 1 of the Gi range gate echo-signal
Rank is related;For the power of the Gi range gate echo-signal, and
0 rank related coefficient parameter Estimation in short-term, detection calculates as the following formula
Wherein,WithIt is that horizontal polarization echo and vertical polarization are returned respectively
0 rank auto-correlation function of wave, i.e. power, andFor 0 rank of horizontal polarization echo and vertical polarization echo
Correlation function, i.e. cross-power.
Currently preferred technical solution is, continuously apart from points N, value 9.
The further technical solution of the present invention is that the step 3 specifically includes:
1) fuzzy quantization judges 4 kinds of parameter Estimations of noise to quantify the numerical value between 0 to 1;
2) weighted sum is calculated as follows that is, by the parameter Estimation weighted sum after fuzzy quantization
NSign=KPi·σPi+KPH·σphase+KNCP·NCP+KρHV·ρhv(0)-------------(11)
NSignFor the mark judged for noise region;
3) noise region judges, use one compares the result N of threshold value with 2) step weighted sumSignIt is compared, is not less than threshold
Value, then be judged as noise region, be not otherwise;
4) 3) noise region power in step is averaged by noise power calculation, obtains the noise electricity under the range sweep
It is flat.
Currently preferred technical solution be compare threshold value value be 0.85.
Beneficial effects of the present invention:
(1) it improves by the steady noise level in steady noise level or an observation cycle, as radar signal-to-noise ratio
The standard of judgement, the largely reduction signal-to-noise ratio (using maximum noise as standard) that may cause, or accidentally by noise as letter
The bad detection (false dismissed rate increases or false alarm rate increases) of number (using maximum noise as standard).
(2) be suitable for pulse radar due to component parameter variation, signal path characteristic drift, antenna rotate, the sun and
Noise bounce situation caused by the reasons such as cosmic radiation.
(3) dynamically improves echo corresponding to each transmitting pulse compared with using quiescent noise level setting method
Signal-to-noise ratio, only solar radiation noise one can adapt to be more than making an uproar for 20dB in real time to the weather radar of X and following wave band
Sound variation.
Detailed description of the invention
Fig. 1 is the real-time dynamic noise power detector block diagram of weather radar provided in an embodiment of the present invention;
Fig. 2 is that noise power changing condition counts schematic diagram in radar return provided in an embodiment of the present invention;
Fig. 3 is that a range sweep echo power provided in an embodiment of the present invention changes schematic diagram;
Fig. 4 (a), Fig. 4 (b), Fig. 4 (c), Fig. 4 (d) are criteria for noise degree of membership schematic diagram provided in an embodiment of the present invention.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, the technical solution below in the present invention carries out clear
Chu is fully described by, it is clear that described embodiments are some of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The present invention is adapted to the complex baseband signal processing without pulse accumulation of digital receiver output.For pulse
It is still effective to carry out same processing for the complex baseband signal of accumulation, but its threshold parameter needs are modified accordingly, noise
Estimation power can be increased slightly, but still more more efficient than conventional noise threshold value setting method.
Treated that result is as Figure 2-3 by the present invention, and Fig. 2 indicates that the echo power variation in a range sweep is shown
It is intended to, the continuous low level part of mark grey overstriking is the signal of noise signal or power very little;It is single after being handled by judgement
A range sweep pulse has been integrated, and obtains the noise power variation fluctuating of radar bearing scan echo as shown in Figure 3
Situation.
Concrete principle is as follows:
1. the continuous weak power detection of each range sweep
Continuous weak signal power detection is carried out to the corresponding range sweep echo-signal of transmitting pulse one by one.
The echo power P of each range gate Gi is calculated firstHiAnd PVi。
If it is Doppler radar, using its linear digital receiver output H (I,
Q), the echo power P on each range gate Gi can directly be calculatedHi, as shown in formula (1).
Here, H (I, Q) is the digital signal of weather radar horizontal polarization echo, in-phase component HI, the Gi distance
The numerical value of door is HIi;Its quadrature component is HQ, the numerical value of the Gi range gate is HQi。
If it is dual polarization Doppler radar, the number using the horizontal polarization echo of its linear digital receiver output is believed
Number H (I, Q), can directly calculate horizontal polarization echo power P on each range gate GiHi, as shown in formula (1).And it is utilized to hang down
Straight polar echo digital signal V (I, Q), can directly calculate vertical polarization echo power P on each range gate GiVi, such as formula (2)
It is shown.
Wherein, vertical polarization echo in-phase component is VI, the numerical value of the Gi range gate is VIi;Its quadrature component is VQ,
The numerical value of the Gi range gate is VQi。
Secondly, weak echo signal in confirmation and tag distance scanning.
According to the running parameter of radar, as working frequency range, bandwidth, antenna gain, receiver gain, radar station height and
Longitude and latitude etc. determines the possibility noise level of radar, mainly includes two parts, i.e. receiver noise and antenna noise.
Receiver noise can estimate receiver noise power level by formula (3) by receiver theoretical calculation:
PRNoise=-114dB+10log10(Bn[MHz])+10·log10(F0)(dBm)------------(3)
Wherein, BnFor receiver bandwidth, unit MHz;F0For the noise coefficient of radar system, can be used under room temperature
Measured value calculate.
Antenna noise then estimates antenna by sun ephemeris and its spectral density according to the height of radar station and longitude and latitude
Maximum solar power P under the conditions of relevant work frequency range, bandwidth and antenna gainSolar.It is horizontal for solar radiation
Polar power and vertical polarization power can consider equal.
Receiver noise power PRNoiseWith solar radiation power PSolarThe sum of, the allowance of 1dB is reserved, radar can be regarded as
Rough noise power, as the power threshold P for judging weak echoN0, i.e. PN0=PRNoise+PSolar+1.For horizontal polarization channel,
Receiver noise power is PHNoise, for vertical polarization channel, receiver noise power is PVNoise, then determined level receives logical
Road and the power threshold of vertical reception channel weak echo are respectively PHN0=PHNoise+PSolar+ 1 and PVN0=PVNoise+PSolar+1。
Third finds out weak echo signal section in range sweep
Compared by power, the section of the relatively small continuous range gate of signal power in range sweep can be found out, by public affairs
Formula (4) and (5) can respectively obtain the weak echo section in horizontal receiving channel and vertical reception channel.
In above-mentioned two formula, M is continuous range gate number.The points M of weak echo section cannot obtain it is very few, generally take 50 to
, it can get better effects at 100 points.Because of numerical value very little, in order to calculate with it is more convenient, echo power value and threshold value can be taken pair
Number, is being handled.
2. 4 kinds of parameter Estimations of noise judgement
The present invention includes that 3 kinds of Doppler radars judge parameter, i.e., phase standard is poor in short-term, short-time rating standard deviation and in short-term
Average normalized coherent power and a kind of dual polarization radar judge parameter, i.e. 0 rank related coefficient in short-term.
Phase standard apart from upper continuous N point is poor
Radar echo signal phase value range is -180 ° to+180 °.In pure noise region, phase is random.It can
To think, random phase obedience is uniformly distributed, and standard deviation is(degree).
On the contrary, in non-noise region, it is variant in signal strength, but phase has certain consistency.Therefore, it calculates
Phase standard is poor, can be used as the judgment basis of noise signal.
It is poor phase standard can be calculated by formula (6).Wherein, N is continuous distance points, referred to as short because N is smaller
When, generally desirable N=9;xphiFor the phase of the Gi range gate echo-signal, μphIt is that this is continuous apart from upper N point phase average
Value, i.e.,
It is carried out on the basis of aforementioned 1st processing when phase standard difference calculates in short-term.In the several of range sweep
The continuous weak echo region of a M point, one by one area, successively carry out N point phase standard difference calculating, and save these calculated results.
Note that be respectively processed if it is dual polarization radar to horizontal receiving channel and vertical reception channel.This processing only needs
Single range sweep (i.e. 1 measurement pulse) can be completed.
3. comprehensive judging module and noise power estimation
4 kinds of parameter estimation results for judging noise are carried out comprehensive descision processing by integrated noise judgment module, every to obtain
Actual noise level in a range sweep.
1) fuzzy quantization
In most cases, judge that 4 kinds of parameter estimation results of noise can be distributed in big value and small value two sides, in this way, mould
The problem of gelatinization, will become simply, to be easily achieved fuzzy quantization process by fuzzy membership functions as shown in Figure 4.
In order to be suitble to all situations, fuzzy membership functions shown in Fig. 4 is easily achieved fuzzy quantization process.
In order to be suitble to all situations, carrying out the easement curve in Fig. 4 once to simplify processing, transition region is reduced to three
It is worth transition to replace, the subordinating degree function of each parametric variable as shown in Table 1 can be obtained.The practical valued space citing of its Parameters variation is such as
Shown in table 2.
1 subordinating degree function of table
The practical valued space citing of 2 Parameters variation of table
Through this process, judge 4 kinds of parameters of noise to quantify the numerical value between 0 to 1.
2) weighted sum
As shown in Figure 1, each parameter weighting is summed.It should be noted that needing to use short for dual polarization radar
When related coefficient parameter, otherwise do not have to;And before or after the average normalized coherent power apart from upper continuous N point then needs to use
Direct impulse data, at this moment needing to buffer could be handled, if feeling trouble, given up and do not had to, change its dependent variable
Weighting coefficient will not cause manifest error.
Calculation formula (11) progress can be used in weighted sum, as a result NSignFor the mark judged for noise region.According to
The value example of the working form of radar, weighting coefficient is as shown in table 3.
NSign=KPi·σPi+KPH·σphase+KNCP·NCP+KρHV·ρhv(0)-------------(11)
The weighting coefficient value of each parametric variable of table 3 is illustrated
3) noise region is adjudicated
Compare threshold value (desirable 0.85) and above-mentioned calculated result N with oneSignIt is compared, is not less than threshold value, is then judged as
Otherwise noise region is not.
4) noise power calculation
5) for being confirmed to be 1 or several noise regions of noise, such as the grey thick line example of Fig. 2, by its power into
Row is average to get the noise level under the range sweep.Thus it is used as benchmark, the signal-to-noise ratio of echo-signal can be calculated naturally
Come.
The inverting of subsequent signal processing and detection physical quantity, is calculated with signal-to-noise ratio, can be become very convenient.
Embodiment 1: the realization of Dual-Polarized Doppler Weather Radar noise level real-time detection
(traditionally dual polarization is become dual-polarization to dual-polarization by meteorological boundary, in general, polarization is the call for electromagnetic wave, partially
Vibration is then the call for light, and the two can be general) Doppler radar, the characteristics of detection due to its full airspace, and be
Binary channels works at the same time, and especially takes into account cost factor in production, causes orientation/elevation rotary joint, ring in feeder line
The performance transient change for flowing the components such as device is serious.For being also unlikely to performance degradation only with Doppler parameter.But it uses
Dual-polarization observed parameter, then performance is difficult to maintain to stablize.
Noise level real-time detection can solve this problem.
Step 1, the horizontal complex baseband signal H (I, Q) that radar two-channel receiver is exported and vertical complex baseband signal V
(I, Q) is converted to power level, as shown in formula (1) and (2).For convenience, logarithmic transformation is carried out to them respectively, obtained
10logPHiAnd 10logPViDBm indicate echo-power signals.
Step 2, the echo-power signals 10logP to adjust the distance in scanningHiAnd 10logPViContinuous weak power inspection is carried out respectively
It surveys, detection benchmark is the preliminary noise power that formula (3) are calculated, and obtains several weak echo sections.
Step 3 carries out noise to the weak echos of these sections and judges 4 parameter detectings, i.e., phase standard is poor in short-term, function in short-term
Rate standard deviation, short-time average normalize coherent power, and the detection of 0 rank related coefficient in short-term.
Step 4 carries out fuzzy quantization to parameter in step 3.
Step 5 is weighted summation to parameter in step 4.
Step 6 carries out noise region judgement to step 5 summed result.
Step 7 carries out noise level calculating to the noise region judged.
So far, that is, the noise level in horizontal polarization channel and vertical polarization channel is respectively obtained.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify the technical solutions described in the foregoing embodiments or equivalent replacement of some of the technical features;
And these are modified or replaceed, technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution spirit and
Range.
Claims (7)
1. the real-time dynamic noise power detecting method of weather radar characterized by comprising
Step 1. detects the corresponding distance of transmitting pulse of Radar Digital Receiver output with continuous weak power detection module one by one
Scan echo signal, and carry out noise power estimation;
Continuous weak power detection module carries out continuous according to distance/time sequencing according to a certain distance range/length of time
Baseband signal power is detected in low signal areas, obtains comprising the continuous weak power echo time section including noise power;
Step 2. short-time rating standard difference module, in short-term phase standard difference module, short-time average normalize coherent power module
And 04 detection modules of rank related coefficient module are according to relatively short certain period of time in short-term, to continuous weak in step 1
Power echo time segment data carry out short-time rating standard deviation, phase standard is poor in short-term, short-time average normalize coherent power with
And the parameter Estimation of 0 rank related coefficient in short-term;
Step 3. carries out blurring and reasoning to the parameter in step 2 with comprehensive judgment module, determines the positioning of pure noise range, calculates
The power level of noise out.
2. the real-time dynamic noise power detecting method of weather radar according to claim 1, which is characterized in that the step
Rapid 1 continuous weak signal power detection includes:
1) the echo power P of each range gate Gi, is calculatedHiAnd PVi;
A. each range gate Gi can directly be calculated using the output H (I, Q) of its linear digital receiver if Doppler radar
On echo power PHi, it is shown below:
H (I, Q) is the digital signal of weather radar horizontal polarization echo, in-phase component HI, the numerical value of the Gi range gate is
HIi, quadrature component HQ, the numerical value of the Gi range gate is HQi;
B. if dual polarization Doppler radar, the digital signal H of the horizontal polarization echo of its linear digital receiver output is utilized
(I, Q) can directly calculate horizontal polarization echo power P on each range gate GiHi, it is shown below:
Using its vertical polarization echo digital signal V (I, Q), vertical polarization echo power on each range gate Gi can be directly calculated
PVi, it is shown below:
Wherein, vertical polarization echo in-phase component is VI, the numerical value of the Gi range gate is VIi;Its quadrature component is VQ, Gi
The numerical value of range gate is VQi;
2), confirm weak echo signal in simultaneously tag distance scanning
According to the running parameter of radar, the noise level of radar is determined, mainly include two parts, i.e. receiver noise and antenna is made an uproar
Sound, wherein receiver noise estimates receiver noise power level by receiver theoretical calculation as the following formula:
PRNoise=-114dB+10log10(Bn[MHz])+10·log10(F0)(dBm)------------(3)
Wherein, BnFor receiver bandwidth, unit MHz, F0For the noise coefficient of radar system, measurement under room temperature is used
Value calculates;
Antenna noise then estimates antenna in phase according to the height of radar station and longitude and latitude by sun ephemeris and its spectral density
Answer the maximum solar power P under the conditions of working frequency range, bandwidth and antenna gainSolar。
3) weak echo signal section in range sweep, is found out
Compared by power, finds out the section of the relatively small continuous range gate of signal power on range sweep face, divide as the following formula
The weak echo section in horizontal receiving channel and vertical reception channel is not obtained;
M is continuous range gate number.
3. the real-time dynamic noise power detecting method of weather radar according to claim 2, which is characterized in that continuous distance
M number of door, value is in 50-100 point.
4. the real-time dynamic noise power detecting method of weather radar according to claim 1, which is characterized in that the step
Phase standard difference parameter, short-time rating standard parameter, short-time average normalize coherent power parameter, in short-term 0 rank phase in short-term in rapid 2
Relationship number method for parameter estimation are as follows:
Phase standard difference parameter Estimation in short-term is detected using following formula and is calculated:
Wherein, N is continuous distance points, xphiFor the phase of the Gi range gate echo-signal;
μphIt is that this is continuous apart from upper N point phase average, i.e.,
The estimation of short-time rating standard deviation criteria, is detected using following formula and is calculated:
Wherein, xPiFor the power of Gi range gate echo-signal;μPIt is that this is continuous apart from upper N point power average value, i.e.,
Short-time average normalizes the estimation of coherent power parameter, and single repetition rate is calculated using following formula detection:
It is related for 1 rank of the Gi range gate echo-signal;For the power of the Gi range gate echo-signal, and
Or irregular repetition rate, it is detected and is calculated using following formula:
For the 1st pulse spacing,For 1 rank of the Gi range gate echo-signal
It is related;For the power of the Gi range gate echo-signal, andFor the 2nd pulse spacing,For
1 rank of the Gi range gate echo-signal is related;For the power of the Gi range gate echo-signal, and
0 rank related coefficient parameter Estimation in short-term, detection calculates as the following formula
Wherein,WithIt is horizontal polarization echo and vertical polarization echo respectively
0 rank auto-correlation function, i.e. power, andIt is related to 0 rank of vertical polarization echo for horizontal polarization echo
Function, i.e. cross-power.
5. the real-time dynamic noise power detecting method of weather radar according to claim 4, which is characterized in that continuous distance
Points N, value 9.
6. the real-time dynamic noise power detecting method of weather radar according to claim 1, which is characterized in that the step
Rapid 3 specifically include:
1) fuzzy quantization judges 4 kinds of parameter Estimations of noise to quantify the numerical value between 0 to 1;
2) weighted sum is calculated as follows that is, by the parameter Estimation weighted sum after fuzzy quantization
NSign=KPi·σPi+KPH·σphase+KNCP·NCP+KρHV·ρhv(0)-------------(11)
NSignFor the mark judged for noise region;
3) noise region judges, use one compares the result N of threshold value with 2) step weighted sumSignIt is compared, is not less than threshold value, then
It is judged as noise region, is not otherwise;
4) 3) noise region power in step is averaged, obtains the noise level under the range sweep by noise power calculation.
7. the real-time dynamic noise power detecting method of weather radar according to claim 6, which is characterized in that compare threshold value
Value is 0.85.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811198843.1A CN109358331B (en) | 2018-10-15 | 2018-10-15 | Real-time dynamic noise power detection method for meteorological radar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811198843.1A CN109358331B (en) | 2018-10-15 | 2018-10-15 | Real-time dynamic noise power detection method for meteorological radar |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109358331A true CN109358331A (en) | 2019-02-19 |
CN109358331B CN109358331B (en) | 2023-03-03 |
Family
ID=65349440
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811198843.1A Active CN109358331B (en) | 2018-10-15 | 2018-10-15 | Real-time dynamic noise power detection method for meteorological radar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109358331B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109613503A (en) * | 2018-12-20 | 2019-04-12 | 中国气象科学研究院 | The Calibration Method and device of radar echo signal |
CN110221259A (en) * | 2019-05-29 | 2019-09-10 | 安徽四创电子股份有限公司 | A kind of Calibration Method of meteorological radar echo intensity |
CN111337962A (en) * | 2020-03-26 | 2020-06-26 | 湖南科技大学 | RTK positioning system in thunder field detection |
CN111474606A (en) * | 2020-03-30 | 2020-07-31 | 中国人民解放军国防科技大学 | Rainfall and water vapor comprehensive measurement device and method based on dual-frequency dual-polarization microwave link |
CN113655455A (en) * | 2021-10-15 | 2021-11-16 | 成都信息工程大学 | Dual-polarization weather radar echo signal simulation method |
CN114442103A (en) * | 2022-01-12 | 2022-05-06 | 成都亘波雷达科技有限公司 | Dual-polarization weather radar differential propagation phase shift rate estimation method, system and equipment |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6307500B1 (en) * | 1999-08-13 | 2001-10-23 | University Corporation For Atmospheric Research | Method of moment estimation and feature extraction for devices which measure spectra as a function of range or time |
WO2011050097A1 (en) * | 2009-10-20 | 2011-04-28 | Colorado State University Research Foundation | Sensitivity enhancement system |
CN103336274A (en) * | 2013-06-27 | 2013-10-02 | 成都信息工程学院 | Two-way multi-order related detection method for dual-polarization weather radar |
CN105527624A (en) * | 2014-09-29 | 2016-04-27 | 北京敏视达雷达有限公司 | Method for dynamically estimating noise through radar echo and meteorological radar system |
CN106093893A (en) * | 2016-06-01 | 2016-11-09 | 西安电子工程研究所 | A kind of online calibration method of any polarized wave of dual polarization radar |
CN106249207A (en) * | 2016-07-11 | 2016-12-21 | 北京敏视达雷达有限公司 | A kind of method and device of estimated noise power |
-
2018
- 2018-10-15 CN CN201811198843.1A patent/CN109358331B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6307500B1 (en) * | 1999-08-13 | 2001-10-23 | University Corporation For Atmospheric Research | Method of moment estimation and feature extraction for devices which measure spectra as a function of range or time |
WO2011050097A1 (en) * | 2009-10-20 | 2011-04-28 | Colorado State University Research Foundation | Sensitivity enhancement system |
CN103336274A (en) * | 2013-06-27 | 2013-10-02 | 成都信息工程学院 | Two-way multi-order related detection method for dual-polarization weather radar |
CN105527624A (en) * | 2014-09-29 | 2016-04-27 | 北京敏视达雷达有限公司 | Method for dynamically estimating noise through radar echo and meteorological radar system |
CN106093893A (en) * | 2016-06-01 | 2016-11-09 | 西安电子工程研究所 | A kind of online calibration method of any polarized wave of dual polarization radar |
CN106249207A (en) * | 2016-07-11 | 2016-12-21 | 北京敏视达雷达有限公司 | A kind of method and device of estimated noise power |
Non-Patent Citations (3)
Title |
---|
VALERY M.MELNIKOV: "Autocorrelation and Cross-Correlation Estimators of Polarimetric Variables", 《JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY》 * |
史万里等: "一种改善双极化多普勒天气雷达弱回波发现能力的方法", 《成都信息工程大学学报》 * |
李啸: "多普勒气象雷达原始回波数据混合压缩算法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109613503A (en) * | 2018-12-20 | 2019-04-12 | 中国气象科学研究院 | The Calibration Method and device of radar echo signal |
CN109613503B (en) * | 2018-12-20 | 2021-01-12 | 中国气象科学研究院 | Calibration method and device for radar echo signals |
CN110221259A (en) * | 2019-05-29 | 2019-09-10 | 安徽四创电子股份有限公司 | A kind of Calibration Method of meteorological radar echo intensity |
CN110221259B (en) * | 2019-05-29 | 2021-11-23 | 安徽四创电子股份有限公司 | Calibration method for meteorological radar echo intensity |
CN111337962A (en) * | 2020-03-26 | 2020-06-26 | 湖南科技大学 | RTK positioning system in thunder field detection |
CN111337962B (en) * | 2020-03-26 | 2024-06-11 | 湖南科技大学 | RTK positioning system in lightning field detection |
CN111474606A (en) * | 2020-03-30 | 2020-07-31 | 中国人民解放军国防科技大学 | Rainfall and water vapor comprehensive measurement device and method based on dual-frequency dual-polarization microwave link |
CN111474606B (en) * | 2020-03-30 | 2022-04-05 | 中国人民解放军国防科技大学 | Rainfall and water vapor comprehensive measurement device and method based on dual-frequency dual-polarization microwave link |
CN113655455A (en) * | 2021-10-15 | 2021-11-16 | 成都信息工程大学 | Dual-polarization weather radar echo signal simulation method |
CN113655455B (en) * | 2021-10-15 | 2022-04-08 | 成都信息工程大学 | Dual-polarization weather radar echo signal simulation method |
CN114442103A (en) * | 2022-01-12 | 2022-05-06 | 成都亘波雷达科技有限公司 | Dual-polarization weather radar differential propagation phase shift rate estimation method, system and equipment |
CN114442103B (en) * | 2022-01-12 | 2022-11-15 | 成都亘波雷达科技有限公司 | Dual-polarization weather radar differential propagation phase shift rate estimation method, system and equipment |
Also Published As
Publication number | Publication date |
---|---|
CN109358331B (en) | 2023-03-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109358331A (en) | The real-time dynamic noise power detecting method of weather radar | |
Lo et al. | Fractal characterisation of sea-scattered signals and detection of sea-surface targets | |
Conte et al. | Covariance matrix estimation for adaptive CFAR detection in compound-Gaussian clutter | |
US6717545B2 (en) | Adaptive system and method for radar detection | |
Yu et al. | Beam multiplexing using the phased-array weather radar | |
Kirincich et al. | Improving HF radar estimates of surface currents using signal quality metrics, with application to the MVCO high-resolution radar system | |
Hennessey et al. | Sea-clutter modeling using a radial-basis-function neural network | |
CN113311416B (en) | Mountain small-basin radar quantitative precipitation estimation method | |
Chen et al. | A method to correct the influence of rain on X-band marine radar image | |
CN110749871B (en) | Parameter estimation method of dual-polarization weather radar | |
Lorente et al. | Evaluating the surface circulation in the Ebro delta (northeastern Spain) with quality-controlled high-frequency radar measurements | |
KR102440355B1 (en) | Deep Learning-based Ocean Cluster Data Measurement System Using Sea Level Wave Reflectance | |
US6771209B1 (en) | Automatic radar target detector | |
Lorente et al. | Long-term skill assessment of SeaSonde radar-derived wave parameters in the Galician coast (NW Spain) | |
CN109100697A (en) | A kind of target condensing method based on ground surveillance radar 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 | |
CN106959442A (en) | Ground wave radar first-order sea echo composes extracting method under strong interference environment based on many domain informations | |
CN114325599B (en) | Automatic threshold detection method for different environments | |
Gu et al. | A speckle noise suppression method based on surface waves investigation and monitoring data | |
El Mashade | Analytical performance evaluation of adaptive detection of fluctuating radar targets | |
Dong et al. | Statistical measures of S-band sea clutter and targets | |
CN107202980B (en) | multi-frame combined sea surface small target detection method based on direction ratio | |
Wang et al. | Clutter suppression and target detection based on biparametric clutter map CFAR | |
CN107748357B (en) | A kind of radar fine day G- Design method of FIR filter | |
Wang et al. | A Persymmetric GLRT for Adaptive Range Spread Target Detection in Non-Homogeneous Environment |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20231128 Address after: No.2666, airport 4th Road, Southwest Airport Economic Development Zone, Shuangliu District, Chengdu, Sichuan 610200 Patentee after: Chengdu genbo Radar Technology Co.,Ltd. Address before: 610225, No. 24, Section 1, Xuefu Road, Southwest Economic Development Zone, Chengdu, Sichuan Patentee before: CHENGDU University OF INFORMATION TECHNOLOGY |