CN109613503A - Radar echo signal calibration method and device - Google Patents
Radar echo signal calibration method and device Download PDFInfo
- Publication number
- CN109613503A CN109613503A CN201811567472.XA CN201811567472A CN109613503A CN 109613503 A CN109613503 A CN 109613503A CN 201811567472 A CN201811567472 A CN 201811567472A CN 109613503 A CN109613503 A CN 109613503A
- Authority
- CN
- China
- Prior art keywords
- signal
- echo
- power
- power spectrum
- initial
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
-
- 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
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/18—Testing or calibrating meteorological apparatus
-
- 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)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Environmental & Geological Engineering (AREA)
- Electromagnetism (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention provides the Calibration Method and device of a kind of radar echo signal, this method comprises: this method comprises: according to the initial noisc of radar system, obtain signal element power, obtain the initial power spectrum of the echo-signal of radar system detection, according to the initial power of echo-signal spectrum and signal element power, extract the meteorologic signal power spectrum in the initial power spectrum of echo-signal, according to meteorologic signal power spectrum, calculate the echo strength after obtaining echo-signal calibration, avoiding influences since radar system is received with ballistic device failure and detection environment, lead to the problem of the echo signal intensity inaccuracy of detection, improve the accuracy of echo strength.
Description
Technical field
The present invention relates to meteorological radar sounding technical field more particularly to the Calibration Methods and dress of a kind of radar echo signal
It sets.
Background technique
Radar is intended to carry out echo strength calibration when detecting sexual intercourse target, what conventional radar echo intensity calibration referred to
It is to be calibrated according to the working characteristics of radar system to the echo-signal of detection.Wherein, the calibration of radar echo intensity is direct
It is related to the accuracy of radar detection.
In the prior art, a kind of Calibration Method is using outer signal source, frequency recorder and power meter, to radar echo signal
Intensity is calibrated；Another Calibration Method is to use direct digital synthesizers (Direct based on automatic calibration technical principle
Digital Synthesis, DDS) technology respectively to transmission power, receive characteristic and calibrated automatically, realize and believe using in machine
Radar echo signal is calibrated in number source.
However, above-mentioned Calibration Method is usually carried out before radar investment detection, due to observing environment, radar emission, reception
The influence of device, original Calibration Data can not provide the time of day of actual observation data, influence the echo letter of radar detection
The accuracy of number intensity.
Summary of the invention
To solve problems of the prior art, the present invention provides the Calibration Method and dress of a kind of radar echo signal
It sets, to carry out calibration to radar echo signal, improves the accuracy of echo signal intensity.
In a first aspect, the embodiment of the present invention provides a kind of Calibration Method of radar echo signal, comprising:
According to the initial noisc of radar system, signal element power is obtained；
Obtain the initial power spectrum of the echo-signal of the radar system detection；
According to the initial power of echo-signal spectrum and the signal element power, the initial of the echo-signal is extracted
Meteorologic signal power spectrum in power spectrum；
According to the meteorologic signal power spectrum, the echo strength after obtaining the echo-signal calibration is calculated.
In a kind of possible embodiment, the initial noisc according to radar system obtains signal element power, packet
It includes:
According to the radar system parameters, the initial noisc power for obtaining the radar system is calculated；
It is counted according to the initial noisc power and the full modal data of the radar system noise signal, obtains the signal
Cell power.
It is described according to the initial power of echo-signal spectrum and the signal element in a kind of possible embodiment
Power extracts the meteorologic signal power spectrum in the initial power spectrum of the echo-signal, comprising:
It is composed according to the initial power of the echo-signal, obtains meteorologic signal in the initial power spectrum of the echo-signal
Power spectral amplitude；
According to the power spectral amplitude of the meteorologic signal and the signal element power, the meteorologic signal power is obtained
Spectrum.
It is described to be composed according to the initial power of the echo-signal in a kind of possible embodiment, obtain the echo
The power spectral amplitude of meteorologic signal in the initial power spectrum of signal, comprising:
It is composed according to the initial power of the echo-signal, obtains the noise level of the initial power spectrum of the echo-signal,
It determines the amplitude of noise signal in the echo-signal, and obtains the initial power spectral amplitude of the echo-signal；
According to the initial power spectral amplitude of the amplitude of the noise signal and the echo-signal, the meteorologic signal is obtained
Power spectral amplitude.
It is described according to the meteorologic signal power spectrum in a kind of possible embodiment, it calculates and obtains the echo letter
Echo strength after number calibration, comprising:
Power spectrum according to the meteorologic signal power spectrum, after obtaining the echo-signal calibration；
Null matrix operation is carried out to the power spectrum after the echo-signal calibration, returning after obtaining the echo-signal calibration
Intensity of wave.
In a kind of possible embodiment, the method also includes:
Obtain the real-time noise power of the radar system；
According to the difference of the real-time noise power and the initial noisc power, the noise function of the radar system is determined
Rate deviation.
In a kind of possible embodiment, the method also includes:
According to the deviation of the noise power, the deviation of the echo strength is obtained.
Second aspect, the embodiment of the present invention provide a kind of calibration device of radar echo signal, comprising:
Module is obtained, for the initial noisc according to radar system, obtains signal element power；
Obtain the initial power spectrum of the echo-signal of the radar system detection；
Extraction module, for according to the initial power of echo-signal spectrum and the signal element power, described in extraction
Meteorologic signal power spectrum in the initial power spectrum of echo-signal；
Processing module calculates after obtaining the echo-signal calibration for the Power Spectrum Distribution according to the meteorologic signal
Echo strength.
In a kind of possible embodiment, module is obtained, is specifically used for:
According to the radar system parameters, the initial noisc power for obtaining the radar system is calculated；
It is counted according to the initial noisc power and the full modal data of the radar system noise signal, obtains the signal
Cell power.
In a kind of possible embodiment, extraction module is specifically used for:
It is composed according to the initial power of the echo-signal, obtains meteorologic signal in the initial power spectrum of the echo-signal
Power spectral amplitude；
According to the power spectral amplitude of the meteorologic signal and the signal element power, the meteorologic signal power is obtained
Spectrum.
In a kind of possible embodiment, extraction module is also used to:
It is composed according to the initial power of the echo-signal, obtains the noise level of the initial power spectrum of the echo-signal,
It determines the amplitude of noise signal in the echo-signal, and obtains the initial power spectral amplitude of the echo-signal；
According to the initial power spectral amplitude of the amplitude of the noise signal and the echo-signal, the meteorologic signal is obtained
Power spectral amplitude.
In a kind of possible embodiment, processing module is specifically used for:
Power spectrum according to the meteorologic signal power spectrum, after obtaining the echo-signal calibration；
Null matrix operation is carried out to the power spectrum after the echo-signal calibration, returning after obtaining the echo-signal calibration
Intensity of wave.
In a kind of possible embodiment, module is obtained, is also used to:
Obtain the real-time noise power of the radar system；
The processing module is also used to the difference according to the real-time noise power and the initial noisc power, determines
The deviation of the noise power of the radar system.
In a kind of possible embodiment, module is obtained, is also used to:
According to the deviation of the noise power, the deviation of the echo strength is obtained.
The third aspect, the embodiment of the present invention provide a kind of calibration device of radar echo signal characterized by comprising
Processor；
Memory；
Computer program；
Wherein, the computer program is stored in the memory, and is configured as being executed by the processor,
The computer program includes the instruction for executing method described in first aspect.
Fourth aspect, the embodiment of the present invention provide a kind of computer readable storage medium, the computer-readable storage medium
Matter is stored with computer program, and the computer program makes server execute method described in first aspect.
The Calibration Method and device of radar echo signal provided in an embodiment of the present invention, this method comprises: according to radar system
The initial noisc of system obtains signal element power, the initial power spectrum of the echo-signal of radar system detection is obtained, according to echo
The initial power spectrum and signal element power of signal, extract the meteorologic signal power spectrum in the initial power spectrum of echo-signal, root
According to meteorologic signal power spectrum, the echo strength after obtaining echo-signal calibration is calculated, is avoided since radar system receives and sends out
Penetrating device fault and detection environment influences, and leads to the problem of the echo signal intensity inaccuracy of detection, improves echo strength
Accuracy.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to do one simply to introduce, it should be apparent that, the accompanying drawings in the following description is this hair
Bright some embodiments for those of ordinary skill in the art without any creative labor, can be with
It obtains other drawings based on these drawings.
Fig. 1 is the flow diagram one of the Calibration Method of radar echo signal provided in an embodiment of the present invention；
Fig. 2 is the flow diagram two of the Calibration Method of radar echo signal provided in an embodiment of the present invention；
Fig. 3 is the flow diagram three of the Calibration Method of radar echo signal provided in an embodiment of the present invention；
Fig. 4 is the flow diagram four of the Calibration Method of radar echo signal provided in an embodiment of the present invention；
Fig. 5 is the structural schematic diagram of the calibration device of radar echo signal provided in an embodiment of the present invention；
Fig. 6 is the hardware structural diagram of the calibration device of radar echo signal provided in an embodiment of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part 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 Calibration Method of radar echo signal provided in an embodiment of the present invention can be applied in radar detection scene.Mesh
Before, radar is intended to carry out echo strength calibration, conventional radar echo intensity calibration refers to basis when detecting sexual intercourse target
The working characteristics of radar system calibrates the echo-signal of detection.Wherein, the calibration of radar echo intensity is directly related to
The accuracy of radar detection.In the prior art, a kind of Calibration Method is using outer signal source, frequency recorder and power meter, to thunder
It is calibrated up to echo signal intensity；Another Calibration Method is to use direct digital synthesizers based on automatic calibration technical principle
(Direct Digital Synthesis, DDS) technology respectively calibrates transmission power, reception characteristic automatically, and realization is adopted
Radar echo signal is calibrated with signal source in machine.However, Calibration Method is usually carried out before radar investment detection, due to
The influence of observing environment, radar emission, receiving device, original Calibration Data can not provide the true shape of actual observation data
State influences the accuracy of the echo signal intensity of radar detection.
The embodiment of the present invention considers the above problem, proposes the Calibration Method and device of a kind of radar echo signal, the party
Method includes: the initial noisc according to radar system, obtains signal element power, obtains the echo signal power spectrum of radar system,
According to echo signal power spectrum and signal element power, the meteorologic signal power spectrum in echo signal power spectrum is obtained, according to gas
Picture signals power spectrum obtains the echo strength of echo-signal, avoids the change due to detecting environment, leads to the strong of echo-signal
The problem of degree fluctuation, improve the accuracy of echo signal intensity.
Executing subject of the invention can be the calibration device of radar echo signal, the device can by software and/or
Hardware realization, the calibration device can integrate in radar system, other electronics that can also be integrated in addition to radar system are set
In standby, the present invention is without limitation.
Technical solution of the present invention is described in detail with specific embodiment below.These specific implementations below
Example can be combined with each other, and the same or similar concept or process may be repeated no more in some embodiments.
Fig. 1 is the flow diagram one of the Calibration Method of radar echo signal provided in an embodiment of the present invention, such as Fig. 1 institute
Show, this method including the following steps:
S101, the initial noisc according to radar system obtain signal element power.
The noise data when initial noisc of radar system is the radar system initial observation, is that determining radar system is initial
The basis of observation state echo signal intensity.In this step, it is believed that the initial noisc power and the radar of radar system
The theoretical noise power of system is consistent, then according to radar system parameters, calculates the initial noisc power for obtaining radar system.
Wherein, the theoretical noise power of radar system can be obtained by the parameter of the radar system, specifically, by formula
(1) it calculates and obtains:
PnsysFor the noise power of radar system, unit dBmW；
K is Boltzmann constant, and K is equal to 1.38 × 10-23Unit is J/K；
T0It is the environment temperature of receiver in the radar system indicated with absolute temperature (when specific implementation, using absolute temperature
Degree, which is brought into formula (1), to be calculated, and the unit of absolute temperature is K), optionally, T0Value can be according to radar system
Actual motion environment temperature is determined, when radar system is run at a constant temperature, can bring formula (1) into using fixed temperature
In, when system is run under non-constant temperature, can be brought into formula (1) using the local surface temperature where radar system,
Such as T0It can be 288K；
NfIt is the noise coefficient of radar system, the technical parameter provided from radar system obtains；
B0It is the bandwidth of receiver in radar system, can be obtained using following formula (2):
B0=c/ (2 × Δ h) (2)
Wherein, Δ h is that range bin is long, and unit is rice；
C is the light velocity, and c is equal to 3 × 108, unit is meter per second.
Then the noise power P of radar system is obtained using formula (1), (2)nsys(dBmW)。
Further, it is counted according to initial noisc power and the full modal data of radar system noise signal, obtains the signal
Cell power.
Specifically, power noise power spectrum of the distal end without meteorologic signal when radar system initial observation is chosen, such as: clear sky
The power spectrum of atmosphere, and multiple repairing weld is carried out to the power spectrum, after rejecting interference signal, select in multiple samples, probability of occurrence
One sample of maximum amplitude, then averages to all amplitudes, to obtain the noise amplitude A of radar systemmean.Thunder
Power spectrum of the distal end without meteorologic signal when noise power spectrum up to system can be radar initial observation.
Then the full modal data points for obtaining radar system noise signal, can be calculated by formula (3):
An=Amean×NFFT (3)
Wherein, AnFor full spectral noise data points, characteristic；
AmeanFor the amplitude of noise signal in radar echo signal；
NFFTNumber, characteristic are converted for Fourier spectrum, when specific implementation, different radar systems, NFFTValue can not
Together, 256,512,1024 etc. be can be set to, this programme is without limitation.
Further, then it can be calculated using formula (4) and obtain signal element power:
Wherein, PFFTFor the power of signal element, unit dBmW.
S102, the initial power spectrum for obtaining the echo-signal that radar system detects.
According to the echo-signal that radar system detects, the initial power spectrum of the echo-signal is obtained, is believed with specific reference to echo
Number obtain echo signal power spectrum implementation have associated description in the prior art, details are not described herein for this programme.
S103, according to the initial power of echo-signal spectrum and signal element power, extract echo-signal initial power spectrum
In meteorologic signal power spectrum.
In radar system real-time detection meteorology target, such as: when sexual intercourse particle, the available corresponding echo of meteorology target
The initial power spectral amplitude of signal, then the signal list according to obtained in the initial power spectral amplitude and S101 of the echo-signal
First power obtains the meteorologic signal power spectrum in the initial power spectrum of echo-signal, it should be appreciated that the echo-signal packet of radar detection
Meteorologic signal and noise signal are included, then the initial power spectrum of echo-signal includes meteorologic signal power spectrum and noise power
Spectrum.
S104, the echo strength according to meteorologic signal power spectrum, after calculating acquisition echo-signal calibration.
Meteorologic signal power spectrum is that the echo signal power of cancelling noise power spectrum signal is composed.It specifically, can be according to obtaining
The meteorologic signal power spectrum taken obtains the power spectrum after echo-signal calibration, is then composed to obtain echo letter according to echo signal power
Echo strength after number calibration.
It should be noted that echo signal power spectrum here is obtained by the meteorologic signal power spectrum after calibration,
Echo signal power spectrum in the S102 detected with radar system is echo signal power spectrum that is different, obtaining here
It is more acurrate.
The Calibration Method of radar echo signal provided in this embodiment obtains signal according to the initial noisc of radar system
Cell power, the echo signal power spectrum for obtaining radar system obtain back according to echo signal power spectrum and signal element power
Meteorologic signal power spectrum in wave power spectrum signal, according to meteorologic signal power spectrum, the echo after obtaining echo-signal calibration is strong
Degree avoids since radar system component problem and detection environment influence, leads to the problem of the intensity inaccuracy of echo-signal, mention
The high accuracy of echo signal intensity.
On the basis of above-mentioned Fig. 1 embodiment, Fig. 2 is the Calibration Method of radar echo signal provided in an embodiment of the present invention
Flow diagram two, as shown in Fig. 2, step S103 can specifically include following steps:
S201, the initial signal power spectrum according to echo obtain the function of meteorologic signal in the initial power spectrum of echo-signal
Rate spectral amplitude.
It is composed according to the initial power of echo-signal, obtains the mode of initial power spectral amplitude of echo-signal in the prior art
There is associated description, details are not described herein.
Then, the noise level for obtaining the initial power spectrum of echo-signal, determines the amplitude of noise signal in echo-signal.
Wherein, noise level is to be separated from the initial power of echo-signal spectrum to noise signal and meteorologic signal
Critical process, after confirming noise level, complete compose under noise level adds up as noise signal, effective letter on noise level
Number section is meteorologic signal, carries out accumulation to the meteorologic signal and calculates the available poem, the meteorologic signal function in the height
Rate spectrum.
When specific implementation, the initial power spectrum of echo-signal is divided into k sections first, k is the integer more than or equal to 1；
Then the average value of every section of power spectrum is determined；The smallest average value in all power spectrum is finally selected, and this is the smallest average
It is worth the amplitude A as noise signalnobs(h,i)。
Read the power spectral amplitude of the echo-signal in each detection height that radar system detects every time, Spectral structure section
From 1,2 ... NFFT,NFFTIt converts and counts for Fourier spectrum, usually 256,512,1024.
Further, according to the initial power spectral amplitude of the amplitude of noise signal and echo-signal, the function of meteorologic signal is obtained
Rate spectral amplitude.
Specifically, meteorologic signal power spectral amplitude is calculated according to formula (5):
Arobs(h, i)=A (h, i)-Anobs(h, i) (i=r0,r1,r2...rm) (5)
Wherein, Arobs(h, i) is the power spectral amplitude of meteorologic signal, is the gas after echo signal power cancelling noise influences
Picture signals amplitude size；
A (h, i) is the initial power spectral amplitude of the echo-signal of actual detection, indicate some time time, h-th in height
Section (r where meteorologic signal in echo-signal0, rm) in each speed interval on signal power data points, be signal width
Size is spent, h is range bin number, and i is effective meteorologic signal section (r0, rm) in corresponding normal-moveout spectrum points.
Data points A in echo-signal involved in above-mentioned formula (5) where meteorologic signal in sectionrobs(h, i), can
By obtaining in a manner of below: respectively record every time in observation data file, the Spectral structure data in each height, record is every
Every spectral line data points A (h, i) (i=r in a height0,r1,r2...rm)。r0To determine meteorologic signal section in Spectral structure
Starting point, 0 < r0<NFFT, NFFTIt counts for the spectral transformation of radar system, usually 256,512,1024；rmIt is true in Spectral structure
Determine meteorologic signal section terminating point, r0<rm<NFFT, NFFTIt counts for the spectral transformation of radar system, usually 256,512,1024.
Effective meteorologic signal section (r involved in above-mentioned formula (5)0, rm), it can be determined by mode below: in function
The maximum value in Spectral structure is detected in rate Spectral structure, is then scanned for, is found close respectively from the left and right sides of the maximum value
The level A of noise signalnobs(h) and power spectral amplitude size is by there is the turning point shaken gradually reducing, and searches at this time
Whole envelopes of signal Spectral structure are obtained.Obtain the Spectral structure section where signal spectrum peak, the start bit in the Spectral structure section
It is set to r0, final position rm, which is m point.When specific implementation, according to above-mentioned data points A (h, i) (i
=1,2 ... NFFT, NFFTConvert and count for Fourier spectrum), r at this time0It can be 0, rmIt can be NFFT。
Optionally, when the height of spectrum peak is differed with the height of noise signal levels less than 3dB, it is determined as that there is no gas
Picture signals.When specific implementation, it is also necessary to carry out the continuous height of identical detection time, the spectral peak position of mutually level adjacent time
The continuity detection set then is considered as detection in vain, without output if there is relatively large deviation.
S202, power spectral amplitude and signal element power according to meteorologic signal obtain meteorologic signal power spectrum.
In this step, the meteorologic signal power spectrum in the echo signal power spectrum of real-time detection is recalculated i.e.
Calibration, specifically:
According to the power spectrum amplitude A of meteorologic signalrobsSignal element power P in (h, i) and S101FFT, using formula
(6) power spectrum of meteorologic signal is obtained:
Probs(h, i)=Arobs(h,i)×NFFT(i=r0,r1,r2...rm) (6)
Wherein, Probs(h, i) is some time after recalculating time, h-th of meteorologic signal spectrum in height, and unit is
dBmW；
NFFTFor radar system signal element watt level, unit dBmW；
Arobs(h, i) is the power spectral amplitude of meteorologic signal, is the meteorological letter after return signal power cancelling noise influences
Number amplitude size indicates some time time, h-th of section (r where meteorologic signal in echo-signal in height0, rmEach of) in
Signal power data points on speed interval, are signal amplitude size, characteristic；H is range bin number, and i is effectively meteorological
Signal spacing (r0, rm) in corresponding normal-moveout spectrum points, (i=r0,r1,r2...rm)。r0To determine meteorologic signal area in Spectral structure
Between starting point, 0 < r0<NFFT, NFFTIt counts for the spectral transformation of radar system, usually 256,512,1024；rmFor in Spectral structure
Determine meteorologic signal section terminating point, r0<rm<NFFT, NFFTIt counts for the spectral transformation of radar system, usually 256,512,1024.
The Calibration Method of radar echo signal provided in this embodiment, is composed according to echo signal power, obtains echo-signal
The power spectral amplitude of meteorologic signal in power spectrum obtains meteorological according to the power spectral amplitude and signal element power of meteorologic signal
Power spectrum signal.The more accurate meteorologic signal power spectrum of radar system can be obtained.
On the basis of above-mentioned Fig. 1, Fig. 2 embodiment, Fig. 3 is the calibration of radar echo signal provided in an embodiment of the present invention
The flow diagram three of method, as shown in figure 3, step S104 can specifically include following steps:
S301, the power spectrum according to meteorologic signal power spectrum, after obtaining echo-signal calibration.
Specifically, can be according to formula (7), power spectrum after echo-signal calibration is calculated:
Wherein, ZiFor the power spectrum after echo-signal calibration；
ProbsFor meteorologic signal power spectrum；
C is radar system constant, by radar system parameters (wherein, PtFor transmission power, G is antenna gain, and θ is level
Beam angle,For vertical beam width, h is distance resolution, | k |2For refractive index, λ is wavelength, L∑For feeder loss) really
It is fixed；
Distance R is meteorological target at a distance from radar antenna, and Δ v is doppler velocity interval.
S302, progress null matrix operation is composed to echo signal power, obtain the echo strength of echo-signal.
Specifically, according to formula (8), the echo strength of echo-signal is obtained:
Wherein, Z is the intensity of radar echo signal.
It should be noted that the radial velocity and speed spectrum width of the particle of radar detection, not with the echo recalculated
The change of the intensity of signal and change, the two parameters are unaffected.
The Calibration Method of radar echo signal provided in this embodiment obtains echo-signal according to meteorologic signal power spectrum
The power spectrum of intensity, to echo signal intensity, power spectrum carry out null matrix operation, obtain the echo strength of echo-signal, keep away
Having exempted to be received due to radar system is influenced with ballistic device failure and detection environment, causes the echo signal intensity of detection inaccurate
True problem improves the accuracy of echo strength.
On the basis of above-mentioned Fig. 1, Fig. 2, Fig. 3 embodiment, Fig. 4 is radar echo signal provided in an embodiment of the present invention
The flow diagram four of Calibration Method, as shown in figure 4, the method also includes:
S401, the real-time noise power for obtaining radar system.
As radar system uses the increase of time, the problems such as necessarily will appear detecting error, the present embodiment thinks initial
Noise power is true value, by comparing the difference of real-time noise power and initial noisc power, obtains the noise function of radar system
The deviation of rate.
Wherein, real-time noise power when real-time noise power can be for without meteorological target, or have meteorological target
Real-time noise power.
Real-time noise power when without meteorological target, the full spectral integral of carry out using certain period without meteorological echo signal are to make an uproar
Acoustical power, such as formula (9).
Wherein, Pnobs0It (h) is the radar real-time noise power provided in data in radar detection moment height, unit
For dBmW；
Pobs0(h, i) be radar real-time detection without meteorological target when echo signal power spectrum.
There is real-time noise power when meteorological target, has the edge of meteorological echo signal without meteorological target point using certain period
Then cloth section uses certain period to have the integral of the power spectrum of the sound section of meteorological echo signal as making an uproar as sound section
Acoustical power.
S402, according to the difference of real-time noise power and initial noisc power, determine the inclined of the noise power of radar system
Difference.
Initial noisc power PnsysCalculation formula may refer to formula (1).Real-time noise function is calculated according to formula (10)
The difference of rate and initial noisc power determines the deviation of the noise power of radar system:
Compare radar theory (initial) noise power PnsysIt is radar system with the noise power difference in GPR Detection Data
The noise power deviation of system.
Δ Pn (dBmW)=Pnobs0(h)-Pnsys (10)
Pn (h) (unit: dBmW) be certain h-th of range bin of moment radar detection on noise power, represent radar this when
The system noise of the radar system at quarter, Δ Pn (unit: dBmW) represent radar-probing system it is initial when with actual detection when thunder
It is this to influence to receive comprising radar system, emit device fault in subsystem up to system to the detection deviation effects of detection data
Deng influence, also the precipitation including radar observation scene causes the environment such as antenna ponding to influence.
Optionally, the method also includes:
S403, the deviation according to noise power, obtain the deviation of echo strength.
Δ Pn (dBmW) is brought into radar weather equation (formula 7) (also just to count Probs(h, i) changes Δ Pn (dBmW) into),
Then echo strength is calculated using formula (8), the echo strength Δ Z (dB) obtained at this time, Δ Z (dB) are that the moment radar is visited
The reception of probing examining system, emission parameter problem lead to measurement error.
The Calibration Method of radar echo signal provided in this embodiment, real-time noise power when obtaining without meteorological target,
According to the difference of real-time noise power and initial noisc power, the deviation of the noise power of radar system is determined, according to noise function
The deviation of rate obtains the deviation of echo strength, can obtain accurate echo strength, while the measurement for obtaining radar system misses
Difference.
Fig. 5 is the structural schematic diagram of the calibration device of radar echo signal provided in an embodiment of the present invention, as shown in figure 5,
The calibration device 50 of the radar echo signal may include: to obtain module 11, extraction module 12, processing module 13, in which:
Module 11 is obtained, for the initial noisc according to radar system, obtains signal element power；
Obtain the initial power spectrum of the echo-signal of the radar system detection；
Extraction module 12, for extracting institute according to the initial power of echo-signal spectrum and the signal element power
State the meteorologic signal power spectrum in the initial power spectrum of echo-signal；
Processing module 13 calculates for the Power Spectrum Distribution according to the meteorologic signal and obtains the echo-signal calibration
Echo strength afterwards.
In a kind of possible embodiment, module 11 is obtained, is specifically used for:
According to the radar system parameters, the initial noisc power for obtaining the radar system is calculated；
It is counted according to the initial noisc power and the full modal data of the radar system noise signal, obtains the signal
Cell power.
In a kind of possible embodiment, module 11 is obtained, is specifically used for:
It is composed according to the initial power of the echo-signal, obtains meteorologic signal in the initial power spectrum of the echo-signal
Power spectral amplitude；
According to the power spectral amplitude of the meteorologic signal and the signal element power, the meteorologic signal power is obtained
Spectrum.
In a kind of possible embodiment, extraction module 12 is also used to:
It is composed according to the initial power of the echo-signal, obtains the noise level of the initial power spectrum of the echo-signal,
It determines the amplitude of noise signal in the echo-signal, and obtains the initial power spectral amplitude of the echo-signal；
According to the initial power spectral amplitude of the amplitude of the noise signal and the echo-signal, the meteorologic signal is obtained
Power spectral amplitude.
In a kind of possible embodiment, processing module 13 is specifically used for:
Power spectrum according to the meteorologic signal power spectrum, after obtaining the echo-signal calibration；
Null matrix operation is carried out to the power spectrum after the echo-signal calibration, returning after obtaining the echo-signal calibration
Intensity of wave.
In a kind of possible embodiment, module 11 is obtained, is also used to:
Obtain the real-time noise power of the radar system；
The processing module 13 is also used to the difference according to the real-time noise power and the initial noisc power, really
The deviation of the noise power of the fixed radar system.
In a kind of possible embodiment, module 11 is obtained, is also used to:
According to the deviation of the noise power, the deviation of the echo strength is obtained.
The calibration device device of radar echo signal provided in an embodiment of the present invention, can execute above method embodiment,
That the realization principle and technical effect are similar is similar for it, and details are not described herein.
Fig. 6 is the hardware structural diagram of the calibration device of radar echo signal provided in an embodiment of the present invention, such as Fig. 6 institute
Show, the calibration device 60 of the radar echo signal of the present embodiment includes: processor 601 and memory 602；Wherein
Memory 602, for storing computer executed instructions；
Processor 601 is set for executing the computer executed instructions of memory storage with realizing to receive in above-described embodiment
Standby performed each step.It specifically may refer to the associated description in preceding method embodiment.
Optionally, memory 602 can also be integrated with processor 601 either independent.
When memory 602 is independently arranged, which further includes bus 603, for connecting the memory
602 and processor 601.
The embodiment of the present invention also provides a kind of computer readable storage medium, stores in the computer readable storage medium
There are computer executed instructions, when processor executes the computer executed instructions, realizes radar echo signal as described above
Calibration Method.
Those of ordinary skill in the art will appreciate that: realize that all or part of the steps of above-mentioned each method embodiment can lead to
The relevant hardware of program instruction is crossed to complete.Program above-mentioned can be stored in a computer readable storage medium.The journey
When being executed, execution includes the steps that above-mentioned each method embodiment to sequence；And storage medium above-mentioned include: ROM, RAM, magnetic disk or
The various media that can store program code such as person's CD.
Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations；To the greatest extent
Pipe present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: its according to
So be possible to modify the technical solutions described in the foregoing embodiments, or to some or all of the technical features into
Row equivalent replacement；And these are modified or replaceed, various embodiments of the present invention technology that it does not separate the essence of the corresponding technical solution
The range of scheme.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811567472.XA CN109613503B (en) | 2018-12-20 | 2018-12-20 | Method and device for calibration of radar echo signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811567472.XA CN109613503B (en) | 2018-12-20 | 2018-12-20 | Method and device for calibration of radar echo signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109613503A true CN109613503A (en) | 2019-04-12 |
| CN109613503B CN109613503B (en) | 2021-01-12 |
Family
ID=66011050
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811567472.XA Active CN109613503B (en) | 2018-12-20 | 2018-12-20 | Method and device for calibration of radar echo signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109613503B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110221259A (en) * | 2019-05-29 | 2019-09-10 | 安徽四创电子股份有限公司 | A kind of Calibration Method of meteorological radar echo intensity |
| CN110726976A (en) * | 2019-10-23 | 2020-01-24 | 森思泰克河北科技有限公司 | Target information calibration method and device and terminal equipment |
| CN111239704A (en) * | 2020-02-12 | 2020-06-05 | 中国科学院大气物理研究所 | Atmospheric detection radar target echo identification processing method, device, equipment and medium |
| CN115356696A (en) * | 2022-07-30 | 2022-11-18 | 中国气象局气象探测中心 | Weather radar signal processing system algorithm calibration method and device |
| CN118348539A (en) * | 2024-06-18 | 2024-07-16 | 无锡智鸿达电子科技有限公司 | A method, system and storage medium for processing multi-band meteorological echo signals |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102707270A (en) * | 2012-06-18 | 2012-10-03 | 武汉大学 | Automatic estimation method for antenna pattern relative to high-frequency ground wave radar |
| CN103941238A (en) * | 2014-05-08 | 2014-07-23 | 西安电子科技大学 | Networked radar collaborative anti-interference transmitting power distribution method |
| CN104067138A (en) * | 2011-09-09 | 2014-09-24 | 阿斯泰克斯有限责任公司 | Imaging radar sensor with synthetic enlargement of the antenna aperture and two-dimensional beam sweep |
| CN104345305A (en) * | 2013-07-30 | 2015-02-11 | 北京敏视达雷达有限公司 | Ground clutter residual filtering method |
| CN106772296A (en) * | 2017-01-20 | 2017-05-31 | 南京大学 | Meteorological radar echo intensity calibration device and method |
| CN107831478A (en) * | 2017-10-30 | 2018-03-23 | 北京无线电测量研究所 | A kind of automatic marking method and system of wind profile radar echo-signal |
| KR20180060799A (en) * | 2016-11-29 | 2018-06-07 | 주식회사 에스원 | Method and apparatus for automatically setting detection degree of radar |
| CN108196233A (en) * | 2017-12-26 | 2018-06-22 | 北京无线电测量研究所 | A kind of real-time online Calibration Method based on millimeter wave cloud detection radar |
| US20180217230A1 (en) * | 2017-02-02 | 2018-08-02 | U.S. Army Research Laboratory Attn: Rdrl-Loc-I | Spectrum sensing for advanced radar performance |
| CN109358331A (en) * | 2018-10-15 | 2019-02-19 | 成都信息工程大学 | Real-time dynamic noise power detection method for meteorological radar |
-
2018
- 2018-12-20 CN CN201811567472.XA patent/CN109613503B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104067138A (en) * | 2011-09-09 | 2014-09-24 | 阿斯泰克斯有限责任公司 | Imaging radar sensor with synthetic enlargement of the antenna aperture and two-dimensional beam sweep |
| CN102707270A (en) * | 2012-06-18 | 2012-10-03 | 武汉大学 | Automatic estimation method for antenna pattern relative to high-frequency ground wave radar |
| CN104345305A (en) * | 2013-07-30 | 2015-02-11 | 北京敏视达雷达有限公司 | Ground clutter residual filtering method |
| CN103941238A (en) * | 2014-05-08 | 2014-07-23 | 西安电子科技大学 | Networked radar collaborative anti-interference transmitting power distribution method |
| KR20180060799A (en) * | 2016-11-29 | 2018-06-07 | 주식회사 에스원 | Method and apparatus for automatically setting detection degree of radar |
| CN106772296A (en) * | 2017-01-20 | 2017-05-31 | 南京大学 | Meteorological radar echo intensity calibration device and method |
| US20180217230A1 (en) * | 2017-02-02 | 2018-08-02 | U.S. Army Research Laboratory Attn: Rdrl-Loc-I | Spectrum sensing for advanced radar performance |
| CN107831478A (en) * | 2017-10-30 | 2018-03-23 | 北京无线电测量研究所 | A kind of automatic marking method and system of wind profile radar echo-signal |
| CN108196233A (en) * | 2017-12-26 | 2018-06-22 | 北京无线电测量研究所 | A kind of real-time online Calibration Method based on millimeter wave cloud detection radar |
| CN109358331A (en) * | 2018-10-15 | 2019-02-19 | 成都信息工程大学 | Real-time dynamic noise power detection method for meteorological radar |
Non-Patent Citations (3)
| Title |
|---|
| WANG KUN-ZHE 等: ""Calibration of Average Scattering Coefficients of Clutter from Radar Baseband Echo"", 《2013 2ND INTERNATIONAL CONFERENCE ON MEASUREMENT, INFORMATION AND CONTROL》 * |
| 王青梅 等: ""一种便携式激光测云仪的云底高度反演方法"", 《强激光与粒子束》 * |
| 阮征: ""基于雷达资料的降水云微物理参数反演及垂直结构研究"", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN110726976A (en) * | 2019-10-23 | 2020-01-24 | 森思泰克河北科技有限公司 | Target information calibration method and device and terminal equipment |
| CN110726976B (en) * | 2019-10-23 | 2022-04-29 | 森思泰克河北科技有限公司 | Target information calibration method and device and terminal equipment |
| CN111239704A (en) * | 2020-02-12 | 2020-06-05 | 中国科学院大气物理研究所 | Atmospheric detection radar target echo identification processing method, device, equipment and medium |
| CN115356696A (en) * | 2022-07-30 | 2022-11-18 | 中国气象局气象探测中心 | Weather radar signal processing system algorithm calibration method and device |
| CN118348539A (en) * | 2024-06-18 | 2024-07-16 | 无锡智鸿达电子科技有限公司 | A method, system and storage medium for processing multi-band meteorological echo signals |
| CN118348539B (en) * | 2024-06-18 | 2024-11-05 | 无锡智鸿达电子科技有限公司 | Multi-band weather echo signal processing method, system and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109613503B (en) | 2021-01-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109613503A (en) | Radar echo signal calibration method and device | |
| US9097805B2 (en) | Systems and methods for calibrating dual polarization radar systems | |
| US9261583B2 (en) | Systems and methods for calibrating dual polarization radar systems | |
| US6894638B2 (en) | Radar signal processing unit and radar signal processing method for abnormal signal extraction | |
| CN106021864A (en) | Calibration method and device of back scattering coefficient of space-borne scatterometer | |
| JP7295317B2 (en) | Noise elimination device and method for weather radar | |
| CN114397636B (en) | Method, system and equipment for evaluating uniformity of reflectivity factor of ground-based radar | |
| CN109884605A (en) | A method for extracting the absorption attenuation and Mie scattering attenuation of cloud and rain to radar signals | |
| CN109254291A (en) | The Doppler power spectra data fusion method and device of cloud radar | |
| RU2550576C1 (en) | Method to measure distance to noisy object | |
| CN102830394B (en) | Weak target detection method based on multispectral accumulation | |
| CN115343681B (en) | A radar detectability evaluation method for aerial targets | |
| JP2012018155A (en) | Radar reception signal processing device and method | |
| KR101885922B1 (en) | Apparatus and system for generating model of clutter using radar video signal and method thereof | |
| EP2416170B1 (en) | Measurement device, measurement system, measurement method, and program | |
| KR101441982B1 (en) | A method for correction of Aliased Radial Velocity through Specified Scan Strategy | |
| CN114966579A (en) | Method and device for acquiring calibration parameters of radar system | |
| CN112444814B (en) | Digital array weather radar signal processor based on PCIE optical fiber acquisition card | |
| CN105527624B (en) | A kind of method and weather radar system of radar return dynamic estimation noise | |
| CN117572467B (en) | Multi-sea-condition tide level observation method and device, electronic equipment and storage medium | |
| RU158491U1 (en) | RADAR METER OF TRACK SPEED OF NON-MANEUVERING AERODYNAMIC GOAL FOR SELECTION OF RANGE SQUARE | |
| CN115856803A (en) | Phased array radar and method and system for testing amplitude and phase of receiving channel of phased array radar | |
| CN116125406A (en) | Method for evaluating performance of space-based surveillance radar based on SNR estimation and track-pointing track report | |
| CN115236418A (en) | Directional diagram correction method and system for unmanned aerial vehicle antenna test | |
| CN115930916B (en) | A Method for Eliminating Effects of Multi-beam Fringe Sounding Based on Spectral Correlation of Hull Surge |
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 |