[go: up one dir, main page]

CN106706566B - A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility - Google Patents

A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility Download PDF

Info

Publication number
CN106706566B
CN106706566B CN201710155843.2A CN201710155843A CN106706566B CN 106706566 B CN106706566 B CN 106706566B CN 201710155843 A CN201710155843 A CN 201710155843A CN 106706566 B CN106706566 B CN 106706566B
Authority
CN
China
Prior art keywords
extinction coefficient
atmospheric
point
vertical
visibility
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710155843.2A
Other languages
Chinese (zh)
Other versions
CN106706566A (en
Inventor
余东升
徐赤东
纪玉峰
方蔚恺
杨喆
张战盈
张伟丽
李艳蕾
陆泽橼
王燕飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei Institutes of Physical Science of CAS
Original Assignee
Hefei Institutes of Physical Science of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hefei Institutes of Physical Science of CAS filed Critical Hefei Institutes of Physical Science of CAS
Priority to CN201710155843.2A priority Critical patent/CN106706566B/en
Publication of CN106706566A publication Critical patent/CN106706566A/en
Application granted granted Critical
Publication of CN106706566B publication Critical patent/CN106706566B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G01N21/538—Scattering, i.e. diffuse reflection within a body or fluid within a flowing fluid, e.g. smoke for determining atmospheric attenuation and visibility
    • 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
    • G01S17/95—Lidar systems specially adapted for specific applications for meteorological use
    • 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

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Immunology (AREA)
  • Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The present invention discloses a kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, comprising the following steps: obtains the original echoed signals of laser radar vertical measurement and judges whether there is cloud in original echoed signals；In cloudless signal, SEQUENCING VERTICAL extinction coefficient is calculated using Fernald method；The biggish point of error on SEQUENCING VERTICAL delustring profile is rejected, and interpolation processing is done to the point；Judge a certain height point on SEQUENCING VERTICAL extinction coefficient profile and its at least continuously have extinction coefficient≤0.1956 of 5 points later, taking the point is atmospheric extinction coefficient boundary height point n*L；Measurement origin is calculated to the atmospheric transmittance T between boundary height point n*L1；According to atmospheric transmittance T1, inverting measurement origin to the Zenith Distance extinction coefficient α between boundary height point n*L；SEQUENCING VERTICAL visibility V is calculated, inventive algorithm is simple, when not having cloud right above laser radar apparatus, can in real time, automatically provide SEQUENCING VERTICAL visibility calculated result.

Description

A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility
Technical field
The present invention relates to the application field of laser radar, the meter of specifically a kind of laser radar detection SEQUENCING VERTICAL visibility Calculation method.
Background technique
Visibility is one of main contents of meteorological element, and the visibility on SEQUENCING VERTICAL direction both can reflect atmosphere Degree of stability, also can be used as determine property of air mass and Pollution Study important factor.The visibility in SEQUENCING VERTICAL direction Measurement can be applied to the fields such as traffic, military affairs, science and technology, provide necessary safety guarantee for takeoff and landing and people's daily life. Currently, conventional visibility meter product is measured primarily directed to horizontal meteorological optical range, do not have vertical visibility measurement Function.Nanjing Information engineering Univ Bu Lingbing, Qiu Zujing, Hou Yuyun et al. are proposed using the forward direction for measuring horizontal meteorological optical range Scattering formula visibility meter and the Laser-ceilometer of vertical measurement measure horizontal visibility value and atmosphere under the conditions of different weather respectively Aerosol vertical extinction coefficient α establishes the corresponding relationship of different visibility values and Aerosol Extinction Coefficients, eventually by anti- The transmitance at different distance height is drilled to solve the vertical visibility of whole atmosphere, this method, which exists, to be needed visibility meter and swash Light ceilometer repeatedly measures calibration simultaneously, and measurement data amount is big, is readily incorporated error, while algorithm is more complicated asks Topic.Laser radar is to pass through quantitative analysis atmospheric aerosol to physical effects such as the scatterings, absorption, delustring of laser according to atmosphere The backscattering echo that (floating dust, particulate matter, flue dust and sand and dust) generate laser come remote sensing Aerosol Extinction Coefficients, The effective tool of mass concentration, cloud layer and visibility information.But its application in terms of visibility measurement is also to be directed to water at present Flat meteorological optical range, the measurement for vertical visibility, there is presently no a kind of particularly effective calculation methods.Therefore, it is necessary to Serious analysis is done from algorithm, is found a kind of calculation method method of effective laser radar detection SEQUENCING VERTICAL visibility, is Laser radar is applied to vertical visibility detection and provides a kind of new approaches.
Summary of the invention
The present invention provides a kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, substantially former based on laser radar Reason, the backscattering echo that laser is generated by atmospheric aerosol particle in quantitative analysis vertical direction, and then obtain atmosphere Aerosol vertical extinction coefficient profile passes through measurement point in calculating vertical direction to the average extinction coefficient between boundary height point α realizes the purpose for calculating SEQUENCING VERTICAL visibility according to the calculation formula of SEQUENCING VERTICAL visibility V.
In order to solve the above-mentioned technical problem, the technical solution adopted by the present invention is that:
A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, comprising the following steps:
(1), the original echoed signals of laser radar vertical measurement are obtained and judge whether there is cloud in original echoed signals；
(2), in cloudless signal, SEQUENCING VERTICAL extinction coefficient profile α is calculated using Fernald method1*L、α2*L、 α3*L..., wherein L is the distance resolution of laser radar；
(3), the biggish point of error on SEQUENCING VERTICAL extinction coefficient profile is rejected, and interpolation processing is done to the point, wherein tool The method of body interpolation is that the SEQUENCING VERTICAL extinction coefficient before and after the point is taken to be weighted and averaged；
(4), judge a certain height point on SEQUENCING VERTICAL extinction coefficient profile and its at least continuously have the delustring of 5 points later Coefficient≤0.1956, i.e. the point and its later horizontal visibility V >=20000m of 5 points, taking the point is atmospheric extinction coefficient side Boundary height point n*L；
(5), measurement origin is calculated to the atmospheric transmittance T between boundary height point n*L1；
(6), according to atmospheric transmittance T1, inverting measurement origin to the Zenith Distance delustring system between boundary height point n*L Number α；
(7), according to calculation formula, SEQUENCING VERTICAL visibility V is calculated.
Preferably, atmospheric transmittance T in the step (5)1Calculation formula it is as follows:
Wherein, α (r) is SEQUENCING VERTICAL extinction coefficient at distance r in formula.
Preferably, in the step (6) origin to the Zenith Distance extinction coefficient α between boundary height point n*L calculating Formula is as follows:
Preferably, the calculation formula of SEQUENCING VERTICAL visibility V is as follows in the step (7):
V=3.912 α-1。
The beneficial effects of the present invention are:
The present invention provides a kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, without setting up visibility meter, When avoiding the atmospheric extinction coefficient that the visibility meter equipment that visibility meter under different visibility conditions is surveyed is surveyed and being demarcated The error of introducing；Algorithm is simple, when there is no cloud right above laser radar apparatus, can in real time, automatically provide SEQUENCING VERTICAL Visibility calculated result；The application range of laser radar has been widened, is provided for the observation of SEQUENCING VERTICAL visibility businessization a kind of new Thinking.
Detailed description of the invention
The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.
Fig. 1 is rudimentary algorithm flow diagram of the invention.
Specific embodiment
It, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical solution in the prior art Embodiment or attached drawing needed to be used in the description of the prior art be briefly described, it should be apparent that, it is described below Attached drawing is some embodiments of the present invention, for those of ordinary skill in the art, before not making the creative labor It puts, is also possible to obtain other drawings based on these drawings.
As shown in Figure 1, the present invention provides a kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility, including following Step:
(1), the original echoed signals of laser radar vertical measurement are obtained and judge whether there is cloud in original echoed signals, by It is much larger than the backscattering cross of aerosol in the backscattering cross of cloud, so when laser transmits encounter cloud in an atmosphere, Its backscattering echo the jump signal increased rapidly will occur, at this time then be to have cloud signal.Since atmospheric aerosol echo is full Sufficient Mie scattering lidar equation is the basis for carrying out SEQUENCING VERTICAL visibility remote sensing survey.When there is cloud, cloud dissipates laser It is emitted back towards wave and is unsatisfactory for Mie scattering lidar equation, the present invention terminates to calculate；
(2), in cloudless signal, SEQUENCING VERTICAL extinction coefficient profile α is calculated using Fernald method1*L、α2*L、 α3*L..., wherein L is the distance resolution of laser radar；
(3), the biggish point of error on SEQUENCING VERTICAL extinction coefficient profile is rejected, and interpolation processing is done to the point, wherein tool The method of body interpolation is that the SEQUENCING VERTICAL extinction coefficient before and after the point is taken to be weighted and averaged；
(4), judge a certain height point on SEQUENCING VERTICAL extinction coefficient profile and its at least continuously have the delustring of 5 points later Coefficient≤0.1956, i.e. the point and its later horizontal visibility V >=20000m of 5 points, taking the point is atmospheric extinction coefficient side Boundary height point n*L；
(5), measurement origin is calculated to the atmospheric transmittance T between boundary height point n*L1；
(6), according to atmospheric transmittance T1, inverting measurement origin to the Zenith Distance delustring system between boundary height point n*L Number α；
(7), according to calculation formula, SEQUENCING VERTICAL visibility V is calculated.
Atmospheric transmittance T in the step (5)1Calculation formula it is as follows:
Wherein, α (r) is SEQUENCING VERTICAL extinction coefficient at distance r in formula.
In the step (6) origin to the Zenith Distance extinction coefficient α between boundary height point n*L calculation formula such as Under:
The calculation formula of SEQUENCING VERTICAL visibility V is as follows in the step (7):
V=3.912 α-1。
The embodiments of the present invention described above are not intended to limit the scope of the present invention, any in the present invention Spirit and principle within made modifications, equivalent substitutions and improvements etc., should be included in claim protection model of the invention Within enclosing.

Claims (2)

1.一种激光雷达探测大气垂直能见度的计算方法，其特征在于，包括以下步骤：1. a calculation method of lidar detection atmospheric vertical visibility, is characterized in that, comprises the following steps: (1)、获取激光雷达垂直测量的原始回波信号并判断原始回波信号中是否有云；(1) Obtain the original echo signal measured vertically by the lidar and determine whether there is a cloud in the original echo signal; (2)、在无云信号时，采用Fernald方法计算大气垂直消光系数廓线α1*L、α2*L、α3*L……，其中L为激光雷达的距离分辨率；(2) When there is no cloud signal, the Fernald method is used to calculate the atmospheric vertical extinction coefficient profiles α 1*L , α 2*L , α 3*L ......, where L is the range resolution of the lidar; (3)、剔除大气垂直消光系数廓线上误差较大的点，并对该点做插值处理，其中，具体插值的方法是取该点前后的大气垂直消光系数做加权平均；(3), remove the point with large error on the atmospheric vertical extinction coefficient profile, and perform interpolation processing on the point, wherein the specific interpolation method is to take the atmospheric vertical extinction coefficient before and after the point to make a weighted average; (4)、判断大气垂直消光系数廓线上某一高度点及其之后至少连续有5个点的消光系数≤0.1956，即该点及其之后5个点的水平能见度V≥20000m，取该点为大气消光系数边界高度点n*L；(4) Judging that the extinction coefficient of a certain height point on the vertical atmospheric extinction coefficient profile and at least 5 consecutive points after it is ≤0.1956, that is, the horizontal visibility of this point and the following 5 points is V≥20000m, take this point is the atmospheric extinction coefficient boundary height point n*L; (5)、计算测量原点至边界高度点n*L之间的大气透过率T1；(5), calculate the atmospheric transmittance T 1 between the measurement origin and the boundary height point n*L; (6)、根据大气透过率T1，反演测量原点至边界高度点n*L之间的平均大气消光系数α；(6) According to the atmospheric transmittance T 1 , invert the average atmospheric extinction coefficient α between the measurement origin and the boundary height point n*L; (7)、根据计算公式，计算大气垂直能见度V，大气垂直能见度V的计算公式如下：(7) According to the calculation formula, calculate the atmospheric vertical visibility V, and the calculation formula of the atmospheric vertical visibility V is as follows: V＝3.912α-1。V=3.912α −1 . 2.根据权利要求1所述的激光雷达探测大气垂直能见度的计算方法，其特征在于,所述步骤(5)中大气透过率T1的计算公式如下：2. the calculation method of lidar detection atmospheric vertical visibility according to claim 1, is characterized in that, in described step (5), the calculation formula of atmospheric transmittance T 1 is as follows: 其中，式中α(r)是距离r处的大气垂直消光系数；where α(r) is the atmospheric vertical extinction coefficient at distance r; 所述步骤(6)中原点至边界高度点n*L之间的平均大气消光系数α的计算公式如下：The calculation formula of the average atmospheric extinction coefficient α between the origin and the boundary height point n*L in the step (6) is as follows:
CN201710155843.2A 2017-03-16 2017-03-16 A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility Active CN106706566B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710155843.2A CN106706566B (en) 2017-03-16 2017-03-16 A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710155843.2A CN106706566B (en) 2017-03-16 2017-03-16 A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility

Publications (2)

Publication Number Publication Date
CN106706566A CN106706566A (en) 2017-05-24
CN106706566B true CN106706566B (en) 2019-05-10

Family

ID=58887063

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710155843.2A Active CN106706566B (en) 2017-03-16 2017-03-16 A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility

Country Status (1)

Country Link
CN (1) CN106706566B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108279221B (en) * 2017-12-07 2021-04-13 中国科学院国家天文台 A method for obtaining atmospheric transparency in local sky area
CN108490451B (en) * 2018-03-29 2022-03-25 中国民航大学 Method for inverting slope visibility by utilizing atmospheric extinction coefficient
CN108627812A (en) * 2018-05-28 2018-10-09 成都信息工程大学 A kind of laser radar atmospheric visibility measurement method and device
CN112505651B (en) * 2020-12-23 2022-06-21 北京遥测技术研究所 Automatic processing method for atmospheric detection laser radar

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4419731A (en) * 1979-03-19 1983-12-06 The Secretary Of State For Defence In Her Britannic Majesty's Government Of The United Kingdom Of Great Britain And Northern Ireland Apparatus for estimating slant visibility in fog
CN1948953A (en) * 2006-10-19 2007-04-18 中国人民解放军海军工程大学 Based on laser rediffusion type atmospheric visibility testing method and testing instrument thereof
CN102590145A (en) * 2012-01-08 2012-07-18 杨少辰 Method for measuring atmospheric visibility by using laser and laser radar visibility meter
CN102621102A (en) * 2012-03-31 2012-08-01 中国科学院安徽光学精密机械研究所 Method for measuring horizontal visibility based on CCD (Charge Coupled Device) laser radar
CN104819963A (en) * 2015-05-07 2015-08-05 南京信息工程大学 Measurement method and monitoring system of atmosphere vertical visibility

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4419731A (en) * 1979-03-19 1983-12-06 The Secretary Of State For Defence In Her Britannic Majesty's Government Of The United Kingdom Of Great Britain And Northern Ireland Apparatus for estimating slant visibility in fog
CN1948953A (en) * 2006-10-19 2007-04-18 中国人民解放军海军工程大学 Based on laser rediffusion type atmospheric visibility testing method and testing instrument thereof
CN102590145A (en) * 2012-01-08 2012-07-18 杨少辰 Method for measuring atmospheric visibility by using laser and laser radar visibility meter
CN102621102A (en) * 2012-03-31 2012-08-01 中国科学院安徽光学精密机械研究所 Method for measuring horizontal visibility based on CCD (Charge Coupled Device) laser radar
CN104819963A (en) * 2015-05-07 2015-08-05 南京信息工程大学 Measurement method and monitoring system of atmosphere vertical visibility

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
激光雷达在气象和大气环境监测中的应用;尹青 等;《气象与环境学报》;20091031;第25卷(第5期);第48-56页 *
用于反演半导体激光云高仪垂直能见度的算法;杨成武 等;《强激光与粒子束》;20120229;第24卷(第2期);第307-311页 *

Also Published As

Publication number Publication date
CN106706566A (en) 2017-05-24

Similar Documents

Publication Publication Date Title
CN106706566B (en) A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility
Muste et al. Practical aspects of ADCP data use for quantification of mean river flow characteristics; part I: moving-vessel measurements
CN110095784B (en) An Ocean-Lower Atmosphere Laser Transport Modeling Method Under the Influence of Complex Environment
CN109827906B (en) Inversion method of laser radar slope visibility
CN110006848B (en) Method and device for obtaining extinction coefficient of aerosol
CN109596594A (en) Based on Raman-Mie scattering lidar Aerosol Extinction inversion method
CN108490451B (en) Method for inverting slope visibility by utilizing atmospheric extinction coefficient
CN1948953B (en) Based on laser rediffusion type atmospheric visibility testing method and testing instrument thereof
CN106707298A (en) Laser radar inversion method of atmospheric transmissivity
CN104819963A (en) Measurement method and monitoring system of atmosphere vertical visibility
Chen et al. Low-level wind effects on the glide paths of the North Runway of HKIA: A wind tunnel study
CN105044039B (en) A kind of method according to laser radar data automatic inversion horizontal visibility
Zhang et al. Comparison of turbulent sensible heat flux determined by large-aperture scintillometer and eddy covariance over urban and suburban areas
Pueschel et al. Visibility and aerosol size frequency distribution
CN109190182B (en) An Electromagnetic Scattering Modeling Method for Nonlinear Sea Surface Covered by Oil Film
CN107784261B (en) A kind of road surface condition detection method and device
CN105891078B (en) Sand devil sand dust mass concentration quantitative inversion evaluation method based on wind profile radar
CN105891064B (en) A kind of detection method and its device of aerosol particle composite rate
WO2001038904A1 (en) Apparatus for detecting turbulent layer
CN116879899A (en) A method based on airborne precipitation particle spectrum inversion
CN110006849A (en) A kind of method and apparatus obtaining Aerosol Extinction
Frisch et al. Observations of boundary-layer convection cells measured by dual-Doppler radar and echosonde, and by microbarograph array
Nagaoka et al. Evaluation of attenuation of ultrasonic wave in air to measure concrete roughness using aerial ultrasonic sensor
CN115541461A (en) Particulate pollution tracing method
JP3571268B2 (en) Fog observation radar device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant