CN108680928A - Measurement method, radar and the storage device of extinction coefficient - Google Patents
Measurement method, radar and the storage device of extinction coefficient Download PDFInfo
- Publication number
- CN108680928A CN108680928A CN201810288767.7A CN201810288767A CN108680928A CN 108680928 A CN108680928 A CN 108680928A CN 201810288767 A CN201810288767 A CN 201810288767A CN 108680928 A CN108680928 A CN 108680928A
- Authority
- CN
- China
- Prior art keywords
- extinction coefficient
- signal
- setpoint distance
- door
- distance door
- 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.)
- Pending
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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/49—Scattering, i.e. diffuse reflection within a body or fluid
- G01N21/53—Scattering, i.e. diffuse reflection within a body or fluid within a flowing fluid, e.g. smoke
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
This application discloses a kind of measurement method of extinction coefficient, this method includes：Coherent Doppler lidar emits laser signal；The echo-signal generated due to the laser signal is received in setpoint distance door, and obtains the power spectral density of echo-signal；The extinction coefficient of setpoint distance door is obtained using power spectral density.Disclosed herein as well is a kind of coherent Doppler lidars and a kind of storage device.By the above-mentioned means, the application can obtain the extinction coefficient of setpoint distance door using coherent Doppler lidar.
Description
Technical field
This application involves field of radar, more particularly to the measurement method of extinction coefficient, radar and storage device.
Background technology
Atmospheric aerosol is the important component of air, the solid-state that is suspended in air and liquid particles object it is total
Claim, the aerodynamic diameter of particle is mostly between 0.001~100 μm.It is important as one in atmospheric physics chemical process
Many fields of factor, atmospheric science are all related with aerosol concentration.
Atmospheric extinction coefficient is the Important Parameters of atmospheric optics, its physical characteristic with aerosol has close pass
System, the spatial distribution of extinction coefficient can reflect the variation of air motion and state to a certain extent.Currently, generally use is non-
Coherent radar measures atmospheric extinction coefficient, however to be easy affected by environment and measurement sensitivity low for above-mentioned measurement method.
Invention content
The application mainly solving the technical problems that provide a kind of measurement method of extinction coefficient, radar and storage device,
The measurement sensitivity of extinction coefficient can be improved and reduce influence of the environment to measurement.
In order to solve the above technical problems, the technical solution that the application uses is：A kind of measurement of extinction coefficient is provided
Method, this method include：Coherent Doppler lidar emits laser signal；It is received due to laser signal in setpoint distance door
The echo-signal of generation, and obtain the power spectral density of echo-signal；The delustring of setpoint distance door is obtained using power spectral density
Coefficient.
In order to solve the above technical problems, another technical solution that the application uses is：It is sharp to provide a kind of coherent Doppler
Optical radar, the coherent Doppler lidar include processor and the memory that is connect with the processor, and memory is for storing
Computer program, processor is for calling computer program to execute the above method.
In order to solve the above technical problems, another technical solution that the application uses is：A kind of storage device is provided, this is deposited
Storage device is stored with computer program, which can be performed to realize the above method.
The advantageous effect of the application is：The case where being different from the prior art, the application use coherent Doppler lidar
Emit laser signal；The echo-signal generated due to laser signal is received in setpoint distance door, and obtains the work(of echo-signal
Rate spectrum density；The extinction coefficient of setpoint distance door is obtained using power spectral density.Therefore compared to incoherent laser radar, improve
Influence of the measurement sensitivity and reduction environment of extinction coefficient to measurement.
Description of the drawings
Fig. 1 is the flow diagram of the measurement method of the application first embodiment extinction coefficient；
Fig. 2 is the flow diagram of the measurement method of the application second embodiment extinction coefficient；
Fig. 3 is the flow diagram of the measurement method of the application 3rd embodiment extinction coefficient；
Fig. 4 is the flow diagram of the measurement method of the application fourth embodiment extinction coefficient；
Fig. 5 is the flow diagram of the measurement method of the 5th embodiment extinction coefficient of the application；
Fig. 6 is the hardware architecture diagram of the embodiment of the present application coherent Doppler lidar；
Fig. 7 is the schematic diagram of the embodiment of the present application storage device.
Specific implementation mode
It is understandable to enable the above objects, features, and advantages of the application to become apparent, below in conjunction with the accompanying drawings, to the application
Specific implementation mode be described in detail.It is understood that specific embodiment described herein is only used for explaining this Shen
Please, rather than the restriction to the application.It also should be noted that illustrating only for ease of description, in attached drawing and the application
Relevant part rather than entire infrastructure.Based on the embodiment in the application, those of ordinary skill in the art are not making creation
Property labour under the premise of all other embodiment for being obtained, shall fall in the protection scope of this application.
Term " first ", " second " in the application etc. be for distinguishing different objects, rather than it is specific suitable for describing
Sequence.In addition, term " comprising " and " having " and their any deformations, it is intended that cover and non-exclusive include.Such as comprising
The step of process of series of steps or unit, method, system, product or equipment are not limited to list or unit, and
It further includes the steps that optionally not listing or unit to be, or further includes optionally for these processes, method, product or equipment
Intrinsic other steps or unit.
Referenced herein " embodiment " is it is meant that a particular feature, structure, or characteristic described can wrap in conjunction with the embodiments
It is contained at least one embodiment of the application.Each position in the description occur the phrase might not each mean it is identical
Embodiment, nor the independent or alternative embodiment with other embodiments mutual exclusion.Those skilled in the art explicitly and
Implicitly understand, embodiment described herein can be combined with other embodiments.In this application, identical expression
The meaning that formula represents is identical.
Below unless stated otherwise, following radars represent coherent Doppler lidar.
Referring to Fig. 1, Fig. 1 is the flow diagram of the measurement method of the application first embodiment extinction coefficient.
Step S11：Coherent Doppler lidar emits laser signal.
Radar can emit laser signal to pre-configured orientation angle and the elevation angle, and azimuth can be east, south, west, north or phase
The combination in adjacent orientation.The elevation angle may range from 0 °~90 °, and more specifically, the range at the elevation angle may be 30 °~60 °, lead to
It crosses between the range at the elevation angle is set as 30 °~60 °, wind speed that coherent Doppler lidar can be made to be obtained and/or disappears
Backscatter extinction logarithmic ratio is more accurate.
Radar at least two pre-configured orientation angles and the elevation angle can also emit laser signal simultaneously or when different.For example,
Radar can emit laser simultaneously on the direction that southeastward orientation of northwest four and the elevation angle are 30 °, 45 ° and 60 ° simultaneously.Example again
Such as, radar can also emit laser within the period 1 to first orientation and first elevation angle, to second orientation within second round
First elevation angle of angle emits laser, emits laser to the third azimuth third elevation angle in the period 3, the period 1, second round with
And the period 3 may be the same or different；First orientation, second orientation and third orientation can be identical, can not also
Together；First elevation angle, second elevation angle and the third elevation angle may be the same or different.Optionally, for emitting laser in radar
The transmitter of signal can rotate, to emit laser signal to different orientations and the different elevations angle in different periods.
Laser signal can be pulse signal, and in other embodiments, laser signal can be non-pulse signal, for example,
Pwm signal.
Step S12：The echo-signal generated due to laser signal is received in setpoint distance door, and obtains echo-signal
Power spectral density.
The extinction coefficient that the distance obtained can be actually needed in the specific numerical value of setpoint distance door according to user determines.If
The value range of set a distance door can be 100m~10km.For example, the value of setpoint distance door can be 100m, 3km or
5km.Setpoint distance door can be in one cycle constant, can also be variation.It is set by what is changed in a cycle
Set a distance door can obtain the extinction coefficient of different set range gate in one cycle.Although should be noted that above-mentioned setpoint distance
Door is a determining numerical value, this numerical value determined is the exact distance that user wants, but in practice,
The distance of the skyborne transmission of laser is not necessarily equal to the setpoint distance door, this distance, which can be more than this, setpoint distance
Door, might be less that the setpoint distance door.
Setpoint distance door indicates the time range of receives echo-signal, and in turn, setpoint distance door can also be by connecing
The time range for receiving echo-signal is determined.It should be noted that the time range can be for the period of a very little, the period
The distance of setpoint distance door can be characterized, for example, time range can use [t, t+ ε] to indicate, wherein t can pass through setpoint distance
Door obtains, and specifically, t can be equal to the ÷ lighies velocity of setpoint distance door × 2.
In one embodiment, radar can only receive the echo-signal at a pre-configured orientation angle and the elevation angle, and obtain
The power spectral density of one echo-signal.
In another embodiment, radar can receive the echo-signal at least two pre-configured orientation angles and the elevation angle, and obtain
Take the power spectral density of at least two echo-signal.
The application radar can be also used for surveying wind, and since wind vector is a trivector, Wind-field Retrieval is at least
Need the radial wind speed on three different pre-configured orientation angles and the elevation angle, that is to say, that radar is when for surveying wind, in a cycle
Including at least three times of reception.Due to the factors such as window foreign matter blocks, cloud and mist blocks, low signal-to-noise ratio may cause data without
Effect, in order to improve the accuracy of measurement, the quantity of the echo-signal at the pre-configured orientation angle and the elevation angle that are received in a cycle can be with
More than three.More specifically, the quantity of the echo-signal at the pre-configured orientation angle and the elevation angle that are received in a cycle can be 4~16
It is a, for example, the quantity of the echo-signal at the pre-configured orientation angle and the elevation angle that are received in a cycle can be 4,8 or 16
It is a.
Since the aerosol particles in air are unevenly distributed, the different orientations and/or difference that radar receives are faced upward
The echo-signal at angle is also different.
The backscatter signal energy N (r) at setpoint distance door r that laser radar receives meets formula (1).
Wherein, P0For transmission power；C is the light velocity；τ is pulse width；ArFor capture area；Y (r) is the several of radar
What overlap factor, the value of Y (r) is related with the structure of radar and setpoint distance door, and value will not be because of the change at azimuth and the elevation angle
Change and changes；β (r) is the backscattering coefficient at distance r, TaIndicate atmospheric transmittance, whereinσ (r) is the atmospheric extinction coefficient at distance r.
Step S13：The extinction coefficient of setpoint distance door is obtained using power spectral density.
Different setpoint distance doors has a different extinction coefficients, setpoint distance Men Yue great, and extinction coefficient is smaller.How to utilize
Power spectral density obtains the extinction coefficient of setpoint distance door, please refers to step S231~S242.
Referring to Fig. 2, Fig. 2 be the application second embodiment extinction coefficient measurement method in the step S13 flows that include
Schematic diagram.In the present embodiment, above-mentioned steps S13 includes following sub-step：
Sub-step S231：CFAR detection is carried out to power spectral density, judges whether the power spectral density can pass through perseverance
False-alarm detects.If the power spectral density cannot go to step S232 by CFAR detection.If the power spectral density can lead to
CFAR detection is crossed, then goes to step S233.
Constant false alarm treatment technology is to provide detection threshold value to inspection policies in radar automatic checkout system and make clutter
The signal processing algorithm minimized is influenced on the false-alarm probability of system with interference.The purpose of radar signal CFAR detection is exactly miscellaneous
Useful target echo is distinguished in wave background.
CFAR detection can be time domain CFAR detection or frequency domain CFAR detection.Since clutter distributed model is divided into
Rayleigh distributed (Rayleigh Distribution), Wei Buer distributions (Weibull distribution) or lognormal
It is distributed (logarithmic normal distribution) etc., therefore, the unit of Rayleigh distributed clutter model can be utilized flat
Equal class constant false alarm algorithm (CA-CFAR), the logarithm constant false alarm algorithm (Log-t-CFAR) using logarithm normal distribution Clutter Model
Or using the methods of ordered statistics class constant false alarm algorithm (OS-CFAR) of Wei Buer Distribution Clutter models to power spectral density
Carry out time domain CFAR detection.Furthermore it is also possible to be averaged constant false alarm algorithm or frequency domain multi-threshold constant false alarm by frequency domain unit
Algorithm carries out frequency domain CFAR detection to power spectral density.
When only receiving the echo-signal of the setpoint distance door at a pre-configured orientation angle and the elevation angle in radar, one can be only obtained
A power spectral density；When radar receives the echo-signal at least two pre-configured orientation angles and elevation angle setpoint distance door, can obtain
To at least two power spectral densities.It should be noted that when the quantity of obtained power spectral density is at least two, each work(is coped with
Rate spectrum density carries out CFAR detection.
Sub-step S232：The power spectral density is invalid.
The power spectral density means that collected power spectral density is undesirable at this, should give rejecting in vain,
It can prevent the power spectral density from being generated to the calculating of the extinction coefficient of setpoint distance door the numerical value zero setting of the power spectral density
It influences, the power spectral density after zero setting is no longer performd the calculation below, i.e., invalid power spectral density is rejected.
Sub-step S233：Maximum power spectral densities are taken to obtained power spectral density.
It, can be to the setpoint distance door that receives when radar receives the echo-signal at a pre-configured orientation angle and the elevation angle
Power spectral density take maximum power spectral densities, when radar receives the echo-signal at least two pre-configured orientation angles and the elevation angle
When, each power spectral density that can dock at least two power spectral densities being received under setpoint distance door is maximized.
Setpoint distance door indicates the time range of receives echo-signal, therefore the power spectral density of an echo-signal includes
Multiple power spectral densities of echo-signal on multiple time points, maximum power spectral densities are in multiple power spectral density
Maximum power spectral densities.
After obtaining maximum power spectral densities, go to step S234.
In other embodiments, the sequence of step S231~S233 is readily modified as S233, S231, S232, you can with first root
Maximum power spectral densities are obtained according to the power spectral density detected, then judge whether maximum power spectral densities can pass through constant false alarm
Detection, in the case of can be by CFAR detection, go to step S234.
Sub-step S234：Using the maximum power spectral densities in echo-signal effective bandwidth, obtaining echo-signal is having
Imitate the power signal in bandwidth.
Specifically power signal can be obtained by formula (2).
Wherein, P (r) is power signal, f0To receive the frequency of echo-signal.B is the band of the echo-signal received
Width, X (f) are maximum power spectral densities.Optionally, P (r) can be identical as the numerical value of above-mentioned or following N (r), can not also
Together.
By being integrated to maximum power spectral densities, the maximum power of effective bandwidth can be obtained, can not only be facilitated
Following calculating, it is also possible that the extinction coefficient of the setpoint distance door arrived is more acurrate.
Sub-step S235：Range calibration is carried out to power signal using setpoint distance door, obtains distance-correcting signal.
Specifically distance-correcting signal can be obtained using formula (3).
S (z)=ln [r2P(r)] (3)
Wherein, r is setpoint distance door, and P (r) is the power signal obtained by step S234.
By obtaining distance-correcting signal, the ambient noise that radar receives can be not only deducted, it is also possible that setting
The operation of the extinction coefficient of set a distance door simplifies.
After obtaining distance-correcting signal, go to step S236.
Sub-step S236：Using distance-correcting signal, pass through Collis Slope Methods, Klett inversion methods or Fernald
Inversion method obtains the extinction coefficient of setpoint distance door.
The specific method of Collis Slope Methods can be first to obtain formula to formula (3) and in conjunction with formula (1) progress differential
(4)。
Wherein, β is atmospheric backscatter coefficient.
Since Collis Slope Methods are premised on atmospheric level is uniformly distributed, σ and β are constant, i.e.,
The expression formula for finally obtaining extinction coefficient σ (r) is formula (5).
Klett inversion methods are specifically, in the presence of such as between setting the atmospheric backscatter coefficient β and extinction coefficient in air
The relationship of formula (6).
β=C σk (6)
Wherein, C is constant, and value range is that 40~60, k is delustring Back-scattering ratio, and value range is 0367~1.0.
The expression formula (7) of extinction coefficient is obtained by formula (1), (5) and (6).
Wherein, rmThe maximum distance or reference distance that can be detected for radar, Sm=S (rm), σm=σ (rm)。
Fernald algorithms consider the scattering process of molecule, and echo-signal is expressed as particulate and atmospheric molecule
It is coefficient as a result, echo be just not only the P (r) in formula (1) only with the relating to parameters of particulate, and also and air
Molecule is related, and the expression formula of the echo-signal energy at specific setpoint distance door r is formula (8).
Wherein, βaFor particulate backscattering coefficient, βmFor atmospheric molecule backscattering coefficient, σa(r) it is aerosol
The extinction coefficient of particle, σm(r) it is the extinction coefficient of atmospheric molecule.
Similarly, the expression formula (9) of extinction coefficient is obtained according to formula (8), (5) and (6).
Wherein, X (r)=p (r) r2；Sa=σa/βa。
It should be noted that when the laser signal of radar emission is one, i.e., when the emission parameter of radar is one, by above-mentioned
Step can obtain the extinction coefficient of a setpoint distance door.It is at least two pre-configured orientation angles in the laser signal of radar emission
And the elevation angle namely radar receive at least two pre-configured orientation angles and elevation angle echo-signal, and the corresponding power spectrum of the echo-signal
When density passes through CFAR detection, the extinction coefficient of two setpoint distance doors can be obtained through the above steps.
After obtaining extinction coefficient, go to step S237.
Sub-step S237：Judge whether obtained extinction coefficient is one.
If so, the S238 that gos to step.
If it is not, the extinction coefficient of the setpoint distance door then represented is at least two, then go to step S239.Its
In, the case where which is at least two can be：In above-mentioned steps S13, coherent Doppler lidar transmitting
Laser signal is at least two, and processing of each laser signal Jing Guo above-mentioned steps obtains an extinction coefficient；Wherein, described every
The emission parameter of a laser signal is different, which is specifically as follows at least one of above-mentioned azimuth and elevation angle.
Sub-step S238：Using the extinction coefficient as the extinction coefficient of setpoint distance door.
Sub-step S239：At least two extinction coefficients of the setpoint distance door to obtaining carry out statistics amalgamation mode processing.
Specifically, statistics amalgamation mode can calculate average value, median, mode, standard deviation and variance etc.
It is at least one.
Sub-step S240：Judge through counting amalgamation mode treated value whether within the scope of the first default value.
If it is not, then jumping to sub-step S241；If so, jumping to sub-step S242.
First default value range can be numberical range set by user, for example, user sets numberical range as (δm,
δn), if the value after counting fusion treatment is in (δm, δn) in the range of, then go to step S242, if not in (δm, δn) model
Sub-step S241 is then jumped in enclosing.For example, when it is to calculate average value to count amalgamation mode, the first default value range can be with
For the first default average value ranges；It is to calculate average value and when variance when counting amalgamation mode, the first default value range can be with
Including the first default average value ranges and the first default variance range, and whether the average value to obtaining and variance judge it respectively
Within the scope of the first default average value ranges and the first default variance.In other embodiments, the first default value range can be with
It is determined according to the extinction coefficient of the setpoint distance door got before radar.For example, radar can identify tested scene, and obtain
Radar is previously achieved the extinction coefficient of setpoint distance door under the tested scene, and obtains first according to the extinction coefficient being previously achieved
Default value range, for example, being δ in pervious extinction coefficientpWhen, the first preset range can be (δp× 0.8, δp×1.2)。
In another embodiment, radar can obtain the extinction coefficient for the setpoint distance door being previously achieved by high in the clouds.
Sub-step S241：Treated through counting amalgamation mode, and value is invalid.
Invalid value is rejected.
Sub-step 242：Treated through counting amalgamation mode, and value is the final extinction coefficient of setpoint distance door.
In the present embodiment, statistics fusion side is carried out by least two extinction coefficients of the setpoint distance door first to obtaining
Formula processing, then the value after counting fusion treatment is judged whether within the scope of the first default value, if not existing, through counting at fusion
Value after reason is invalid, if in value after the counting fusion treatment extinction coefficient final as setpoint distance door.Therefore, exist
When the obtained value after counting fusion treatment is invalid, then it is invalid to represent this measurement result, needs radar to remeasure, therefore
It can make raising detection accuracy.
Referring to Fig. 3, Fig. 3 is the flow diagram of the measurement method of the application 3rd embodiment extinction coefficient.
It is in place of the difference of the present embodiment and second embodiment：Sub-step S239~S242 in second embodiment can be with
Replace with sub-step S339~S340.It can first judge whether each extinction coefficient is pre- second at least two extinction coefficients
If in numberical range, and statistics amalgamation mode processing is carried out to the extinction coefficient within the scope of the second default value, set
The final extinction coefficient of range gate.Referring specifically to sub-step S339~S340.
Sub-step S339：From at least two extinction coefficients of setpoint distance door, obtains and be in the second default value range
Interior extinction coefficient.
, can be similar with the first preset range about the value of the second default value range, details are not described herein again.This
Two default value ranges can be identical as the first default value range, can also be different, for example, the second default value range can be with
More than the first default value range or it is less than the first default value range.
Sub-step S340：Statistics amalgamation mode processing is carried out to the extinction coefficient within the scope of the second default value, and
The extinction coefficient final as setpoint distance door will be worth after statistics amalgamation mode processing.
Wherein, for the extinction coefficient not within the scope of the second default value, then it is assumed that it is invalid.The extinction coefficient is invalid
It means, on the emission parameter and setpoint distance door by the extinction coefficient that the echo-signal received calculates is not to be inconsistent
It closes and requires, which can be rejected, to the extinction coefficient after rejecting without statistics amalgamation mode processing.
In the present embodiment, first at least two extinction coefficients of the setpoint distance door to obtaining judge, and to being in
Extinction coefficient within the scope of second default value carries out statistics amalgamation mode processing, to not disappearing within the scope of the second default value
Backscatter extinction logarithmic ratio is rejected.Due to eliminating the extinction coefficient not met within the scope of the second default value, remain pre- in second
If the value in numberical range, therefore, it is possible to avoid not meeting the extinction coefficient within the scope of the second default value to obtain setting away from
The extinction coefficient final from door has an impact, and keeps measurement result more accurate.
Referring to Fig. 4, Fig. 4 is the flow diagram of the measurement method of the application fourth embodiment extinction coefficient.
It is in place of the difference of the present embodiment and 3rd embodiment：Sub-step S341 in 3rd embodiment could alternatively be
Sub-step S441~S443.Since air is unevenly distributed, under normal circumstances, the obtained extinction coefficient in the different elevations angle
It is different.The present embodiment can obtain the delustring parameter at the different elevations angle.
Sub-step S441：Judge in the elevation angle corresponding to the extinction coefficient within the scope of the second default value, if having phase
The same elevation angle.
If so, redirecting S442.If it is not, jumping to S443.
Sub-step S442：Extinction coefficient corresponding to the identical elevation angle is subjected to statistics amalgamation mode processing.And it will statistics
Amalgamation mode treated extinction coefficient of the extinction coefficient as the elevation angle and the setpoint distance door.
The restriction for counting amalgamation mode can be refering to above description.
Sub-step S443：Using the disappearing as the difference elevation angle and the setpoint distance door of the extinction coefficient corresponding to the different elevations angle
Backscatter extinction logarithmic ratio.
If for example, there are five extinction coefficients in the second preset range, σ1、σ2、σ3、σ4And σ5, wherein σ1、
σ2、σ3The corresponding elevation angle is α1, σ4The corresponding elevation angle is α2, σ5The corresponding elevation angle is α3.It is then α at setpoint distance door and the elevation angle1When
Corresponding extinction coefficient can be (σ1+σ2+σ3)/3 are α at setpoint distance door and the elevation angle2When corresponding extinction coefficient can
Think σ4, and at setpoint distance door and the elevation angle be α3When corresponding extinction coefficient can be σ5。
Referring to Fig. 5, Fig. 5 is the flow diagram of the measurement method of the 5th embodiment extinction coefficient of the application.
It is in place of the difference of the present embodiment and 3rd embodiment：Sub-step S341 in 3rd embodiment could alternatively be
Sub-step S541~S543.Since air is unevenly distributed, under normal circumstances, the obtained delustring system in different azimuths
Number is different.The present embodiment can obtain the delustring parameter of different orientations.
Sub-step S541：Judge in the azimuth corresponding to the extinction coefficient within the scope of the second default value, if having
Identical azimuth.
If so, jump procedure S542.If it is not, the S543 that gos to step.
Sub-step S542：Extinction coefficient corresponding to identical azimuth is subjected to statistics amalgamation mode processing.And it will system
Count amalgamation mode treated extinction coefficient of the extinction coefficient as the elevation angle and the setpoint distance door.
Sub-step S543：Using the extinction coefficient corresponding to different orientations as the different orientations and the setpoint distance door
Extinction coefficient.
If for example, there are five extinction coefficients in the second preset range, σ1、σ2、σ3、σ4And σ5, wherein σ1、
σ2、σ3Corresponding azimuth is γ1, σ4Corresponding azimuth is γ2, σ5Corresponding azimuth is γ3.Then in setpoint distance door and
Azimuth is γ1When corresponding extinction coefficient can be (σ1+σ2+σ3)/3 are γ at setpoint distance door and the elevation angle2When it is corresponding
Extinction coefficient can be σ4, and at setpoint distance door and the elevation angle be γ3When corresponding extinction coefficient can be σ5。
In addition, the fourth embodiment of the application or the 5th embodiment can not only carry out on the basis of 3rd embodiment,
Can also be carried out on the basis of second embodiment, you can with first obtain the different elevations angle and setpoint distance door extinction coefficient and/
Or the extinction coefficient at different azimuths and setpoint distance door, then whether the extinction coefficient judged effective.
Referring to Fig. 6, Fig. 6 is the hardware architecture diagram of the embodiment of the present application coherent Doppler lidar.
Coherent Doppler lidar 60 include processing module 61, seed laser 62, transmitter 63, receiver 64 and
Memory 65.Linear polarization continuous laser of the seed laser 62 for generation, transmitter 63 are used for emitting laser, receiver 64
In receives echo-signal.Processing module 61 is for handling the echo-signal received；Memory 65 is calculated for storing
Machine program.
In the present embodiment, processing module 61 may include again：Processor 611, coupler 612, balanced detector 613 with
And Data Acquisition Card 614.Coupler 612 is for being allocated the echo-signal received；Balanced detector 613 is used for altogether
Mould signal plays inhibiting effect, it is suppressed that the noise of local oscillator light, and improve the power of useful signal；Data collecting card 614 is used
In the computer extender card for realizing data acquisition (DAQ) function；Processor 611 is for calling computer program to execute above-mentioned
The measurement method for an embodiment of anticipating.Seed laser 62, transmitter 63, receiver 64 and memory 65 with processing module 61
Connection, specially connect with the processor 611 in processing module 61.In addition, coupler 612, balanced detector 613 and number
Capture card 614 is also connect with processor 611.
In the present embodiment, transmitter 63 and receiver 64 are coaxially disposed, in other embodiments, transmitter 63 and reception
Device 64 can not be coaxially disposed.
Referring to Fig. 7, Fig. 7 is the schematic diagram of the embodiment of the present application storage device.
Storage device 70 is stored with computer program, which can be performed to realize above-mentioned any one implementation
The measurement method of the extinction coefficient of example.
Optionally, storage device 70 can be USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory),
Random access memory (RAM, Random Access Memory), magnetic disc, CD or server etc. are various can to store journey
The medium of sequence code.
Optionally, which can also be the memory 65 in above-described embodiment.
The case where being different from the prior art, the application emit laser signal using coherent Doppler lidar；It is setting
The echo-signal generated due to laser signal is received in range gate, and obtains the power spectral density of echo-signal；Utilize power spectrum
Density obtains the extinction coefficient of setpoint distance door.Therefore compared to incoherent laser radar, the measurement for improving extinction coefficient is sensitive
The influence of degree and reduction environment to measurement.
It further, can not only be by the echo-signal that receives for the coherent Doppler lidar for surveying wind
Obtain the wind speed information in air, additionally it is possible to obtain the extinction coefficient of setpoint distance door, expand coherent Doppler lidar
Application range, make user do not increase it is other measure extinction coefficients equipment in the case of, only by using coherent Doppler
Laser radar can be obtained by the extinction coefficient of setpoint distance door.
The foregoing is merely presently filed embodiments, are not intended to limit the scope of the claims of the application, every to utilize this
Equivalent structure or equivalent flow shift made by application specification and accompanying drawing content is applied directly or indirectly in other relevant
Technical field includes similarly in the scope of patent protection of the application.
Claims (10)
1. a kind of measurement method of extinction coefficient, which is characterized in that the method includes：
The coherent Doppler lidar emits laser signal；
The echo-signal generated due to the laser signal is received in setpoint distance door, and obtains the power of the echo-signal
Spectrum density；
The extinction coefficient of the setpoint distance door is obtained using the power spectral density.
2. according to the method described in claim 1, it is characterized in that, it is described using the power spectral density obtain it is described setting away from
Extinction coefficient from door includes：
CFAR detection is carried out to the power spectral density, and can pass through the CFAR detection in the power spectral density
When, obtain the extinction coefficient of the setpoint distance door using the power spectral density after the CFAR detection；
Wherein, the CFAR detection is time domain CFAR detection or frequency domain CFAR detection.
3. according to the method described in claim 2, it is characterized in that,
The power spectral density using after the CFAR detection obtains the extinction coefficient packet of the setpoint distance door
It includes：
Using the power spectral density in the echo-signal effective bandwidth after the CFAR detection, obtain described
Power signal of the echo-signal in the effective bandwidth；
The extinction coefficient of the setpoint distance door is obtained using the power signal.
4. according to the method described in claim 3, it is characterized in that, described obtain the setpoint distance using the power signal
Door extinction coefficient include：
Range calibration is carried out to the power signal using the setpoint distance door, obtains distance-correcting signal；
The extinction coefficient of the setpoint distance door is obtained using the distance-correcting signal.
5. according to the method described in claim 4, it is characterized in that, described obtain the setting using the distance-correcting signal
The extinction coefficient of range gate includes：
It is obtained by Collis Slope Methods, Klett inversion methods or Fernald inversion methods using the distance-correcting signal
To the extinction coefficient of the setpoint distance door.
6. according to Claims 1 to 5 any one of them method, which is characterized in that
The laser signal of the coherent Doppler lidar transmitting is at least two；Wherein, the hair of each laser signal
Penetrate parameter difference；
At least two laser signal is handled to obtain at least two of the setpoint distance door using the above method and is disappeared
Backscatter extinction logarithmic ratio；
The method further includes：
It determines that the setpoint distance door is final according to the numerical value situation of at least two extinction coefficients of the setpoint distance door to disappear
Backscatter extinction logarithmic ratio.
7. according to the method described in claim 6, it is characterized in that, the emission parameter includes azimuth and the elevation angle.
8. the method according to the description of claim 7 is characterized in that at least two delustrings according to the setpoint distance door
The numerical value situation of coefficient determines that the final extinction coefficient of the setpoint distance door includes：
Calculate at least one of the average value, median, mode of at least two extinction coefficients of the setpoint distance door；
Judge the result of calculation whether within the scope of default value；
If using the extinction coefficient that the result of calculation is final as the setpoint distance door.
9. a kind of coherent Doppler lidar, which is characterized in that the radar includes processor and is connected to the processor
Memory, the memory is used for computer program, and the processor is for calling the computer program with perform claim
It is required that the method described in 1-8 any one.
10. a kind of storage device, which is characterized in that the storage device is stored with computer program, the computer program energy
It is enough performed to realize the method described in claim 1-8 any one.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810288767.7A CN108680928A (en) | 2018-03-30 | 2018-03-30 | Measurement method, radar and the storage device of extinction coefficient |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810288767.7A CN108680928A (en) | 2018-03-30 | 2018-03-30 | Measurement method, radar and the storage device of extinction coefficient |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108680928A true CN108680928A (en) | 2018-10-19 |
Family
ID=63799745
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810288767.7A Pending CN108680928A (en) | 2018-03-30 | 2018-03-30 | Measurement method, radar and the storage device of extinction coefficient |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108680928A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110006849A (en) * | 2019-05-05 | 2019-07-12 | 北京是卓科技有限公司 | A kind of method and apparatus obtaining Aerosol Extinction |
| CN110006848A (en) * | 2019-05-05 | 2019-07-12 | 北京是卓科技有限公司 | A kind of method and apparatus obtaining Aerosol Extinction |
| CN111609771A (en) * | 2020-04-19 | 2020-09-01 | 北京理工大学 | Laser fuse spacing method in aerosol environment |
| CN113631953A (en) * | 2019-03-20 | 2021-11-09 | 美国西门子医疗系统股份有限公司 | Quantitative ultrasound using fundamental and harmonic signals |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101581786A (en) * | 2009-06-23 | 2009-11-18 | 中国科学院安徽光学精密机械研究所 | Semiconductor laser radar visibility meter |
| US20100026981A1 (en) * | 2008-07-07 | 2010-02-04 | Consiglio Nazionale Delle Ricerche - Infm Istituto Nazionale Per La Fisica Della Materia | Elastic backscattering and backreflection lidar device for the characterization of atmospheric particles |
| CN104007445A (en) * | 2014-06-09 | 2014-08-27 | 南京中科神光科技有限公司 | All-fiber laser radar aerosol detecting device |
| CN107421917A (en) * | 2017-05-17 | 2017-12-01 | 南京信息工程大学 | A kind of multifunction high-precision atmosphere visibility meter and visibility measurement method |
-
2018
- 2018-03-30 CN CN201810288767.7A patent/CN108680928A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100026981A1 (en) * | 2008-07-07 | 2010-02-04 | Consiglio Nazionale Delle Ricerche - Infm Istituto Nazionale Per La Fisica Della Materia | Elastic backscattering and backreflection lidar device for the characterization of atmospheric particles |
| CN101581786A (en) * | 2009-06-23 | 2009-11-18 | 中国科学院安徽光学精密机械研究所 | Semiconductor laser radar visibility meter |
| CN104007445A (en) * | 2014-06-09 | 2014-08-27 | 南京中科神光科技有限公司 | All-fiber laser radar aerosol detecting device |
| CN107421917A (en) * | 2017-05-17 | 2017-12-01 | 南京信息工程大学 | A kind of multifunction high-precision atmosphere visibility meter and visibility measurement method |
Non-Patent Citations (1)
| Title |
|---|
| 朱浩: "相干激光雷达数据处理优化及应用研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113631953A (en) * | 2019-03-20 | 2021-11-09 | 美国西门子医疗系统股份有限公司 | Quantitative ultrasound using fundamental and harmonic signals |
| CN110006849A (en) * | 2019-05-05 | 2019-07-12 | 北京是卓科技有限公司 | A kind of method and apparatus obtaining Aerosol Extinction |
| CN110006848A (en) * | 2019-05-05 | 2019-07-12 | 北京是卓科技有限公司 | A kind of method and apparatus obtaining Aerosol Extinction |
| CN110006849B (en) * | 2019-05-05 | 2021-08-31 | 北京是卓科技有限公司 | Method and device for obtaining extinction coefficient of aerosol |
| CN110006848B (en) * | 2019-05-05 | 2022-02-01 | 北京是卓科技有限公司 | Method and device for obtaining extinction coefficient of aerosol |
| CN111609771A (en) * | 2020-04-19 | 2020-09-01 | 北京理工大学 | Laser fuse spacing method in aerosol environment |
| CN111609771B (en) * | 2020-04-19 | 2021-05-18 | 北京理工大学 | Laser fuse spacing method in aerosol environment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108680928A (en) | Measurement method, radar and the storage device of extinction coefficient | |
| Bianco et al. | Convective boundary layer depth: Improved measurement by Doppler radar wind profiler using fuzzy logic methods | |
| Eresmaa et al. | Mixing height determination by ceilometer | |
| CN204946249U (en) | Smog and fire-alarm | |
| Löhnert et al. | Profiling cloud liquid water by combining active and passive microwave measurements with cloud model statistics | |
| US7336222B2 (en) | System and method for measuring characteristics of a continuous medium and/or localized targets using multiple sensors | |
| FI125286B (en) | Arrangement and method for detection of freezing | |
| CN105044039B (en) | A kind of method according to laser radar data automatic inversion horizontal visibility | |
| JP2005521040A (en) | Target signature calculation and recognition system and method | |
| Gao et al. | Observations of the atmospheric boundary layer across the land–sea transition zone using a scanning Mie lidar | |
| Myagkov et al. | Evaluation of the reflectivity calibration of W-band radars based on observations in rain | |
| CN113885022A (en) | Fall detection method and radar equipment | |
| CN109884605A (en) | Extracting method of the sexual intercourse to attenuation by absorption and the Mie scattering decaying of radar signal | |
| Chen et al. | Parameterization of wind evolution using lidar | |
| CN116990773A (en) | Low-speed small target detection method and device based on self-adaptive threshold and storage medium | |
| CN114518583A (en) | Multi-sensor data fusion detection early warning system and method | |
| Gasch et al. | An LES-based airborne Doppler lidar simulator and its application to wind profiling in inhomogeneous flow conditions | |
| Antoniou et al. | On the theory of SODAR measurement techniques | |
| JP7425896B2 (en) | Method for determining the spin of a projectile | |
| CN116879899A (en) | Method based on aerial precipitation particle spectrum inversion | |
| CN112455503A (en) | Train positioning method and device based on radar | |
| Hort et al. | In situ observation of dome instabilities at Merapi volcano, Indonesia: A new tool for volcanic hazard mitigation | |
| Martin | Evaluation of Doppler radar ball tracking and its experimental uses | |
| Dekker et al. | Radio acoustic sounding with a UHF volume imaging radar | |
| Wu et al. | Radar target detect using particle filter |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181019 |
|
| RJ01 | Rejection of invention patent application after publication |