CN113433569B - Atmospheric parameter measurement method - Google Patents

Atmospheric parameter measurement method Download PDF

Info

Publication number
CN113433569B
CN113433569B CN202110707189.8A CN202110707189A CN113433569B CN 113433569 B CN113433569 B CN 113433569B CN 202110707189 A CN202110707189 A CN 202110707189A CN 113433569 B CN113433569 B CN 113433569B
Authority
CN
China
Prior art keywords
atmospheric
laser radar
atmosphere
calculated
rcs
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
CN202110707189.8A
Other languages
Chinese (zh)
Other versions
CN113433569A (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.)
Shenzhen Darsun Laser Technology Co ltd
Original Assignee
Shenzhen Darsun Laser Technology Co ltd
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 Shenzhen Darsun Laser Technology Co ltd filed Critical Shenzhen Darsun Laser Technology Co ltd
Priority to CN202110707189.8A priority Critical patent/CN113433569B/en
Publication of CN113433569A publication Critical patent/CN113433569A/en
Application granted granted Critical
Publication of CN113433569B publication Critical patent/CN113433569B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S17/95—Lidar systems specially adapted for specific applications for meteorological use
    • 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/01—Arrangements or apparatus for facilitating the optical investigation
    • 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/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
    • 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/01—Arrangements or apparatus for facilitating the optical investigation
    • G01N2021/0106—General arrangement of respective parts
    • G01N2021/0112—Apparatus in one mechanical, optical or electronic block
    • 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)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

The application discloses an atmospheric parameter measurement method, which receives reflected light beams with different angles through a laser radar. By analysing the echo signals P of these reflected beams ij Is an optical power signal RCS of (a) ij And angle alpha ij The current atmospheric extinction coefficient R (z) is calculated. Based on the atmospheric extinction coefficient R (z) and the optical power signal RCS ij The formula calculates other atmospheric parameters such as back scattering volume coefficient beta (z), optical thickness AOD (z), and transmittance T i (z) and/or visibility V (z). The method does not need to assume the laser radar ratio, and obtains more accurate atmospheric parameters through a multi-angle scanning method and an inversion method.

Description

Atmospheric parameter measurement method
Technical Field
The application relates to the field of laser radar equipment, in particular to an atmospheric parameter measurement method.
Background
With the continuous development of human society, technology is continuously improved, and the damage to natural environment, especially to the atmosphere, of people is more and more serious. The atmospheric environment is of great importance to human beings, and meanwhile, in modern society, a large number of people live in cities, and the urban atmospheric environment directly influences the health of people, so that real-time monitoring of the atmospheric environment is necessary.
The existing atmospheric parameter detection method generally uses a Mie scattering laser radar as an active remote sensing technology, and obtains atmospheric parameters by emitting laser pulses to interact with atmospheric particles and receiving generated echo signals.
The Fernald algorithm is often used to solve the lidar equation to invert the extinction coefficient, optical thickness and transmittance for the atmosphere. The Fernald method requires the inversion solution to be obtained by assuming the lidar ratio of the aerosol. However, the lidar ratio is related to a number of factors, such as the refractive index, size, morphology, composition, etc. of the aerosol particles, and the refractive index, size, morphology, and composition parameters of the actual aerosol particles are greatly different, so that the corresponding lidar ratio is difficult to determine, and thus the error of the atmospheric parameter result calculated by the Fernald algorithm is very large.
Disclosure of Invention
It is an object of the present application to provide an atmospheric parameter measurement device and method which can ameliorate the above problems.
Embodiments of the present application are implemented as follows:
in a first aspect, the present application provides an atmospheric parameter measurement device comprising:
a processor, a laser radar, and a rotating device;
the laser radar is arranged in a direction away from the ground;
the direction perpendicular to the ground is a height direction Z, the laser radar is fixed on the rotating device, the rotating device is used for rotating the laser radar to receive reflected light beams with different angles, and the reflected light beams form different angles with the height direction Z;
the processor is electrically connected with the laser radar.
It can be appreciated that the atmospheric parameter measuring device provided by the application rotates the laser radar through the rotating device to receive reflected light beams with different angles under the condition that the position of the laser radar is unchanged. The processor processes the echo signal data of each reflected light beam and calculates the atmospheric parameters of the current environment.
Wherein the rotating means as shown in fig. 1 and 2 comprises a horizontal plane rotating means and a vertical plane rotating means; the horizontal plane rotating device is used for rotating the laser radar around the height direction Z on a plane parallel to the reflecting flat plate so as to realize the rotation of the laser radar on an XY plane; the vertical plane rotating device is used for rotating the laser radar on a vertical plane vertical to the ground around the vertical axis of the vertical plane so as to realize the rotation of the laser radar on the vertical planes such as the XZ plane and the like.
In alternative embodiments of the present application, the atmospheric parameters include atmospheric extinction coefficient R (z), backscattering volume coefficient β (z), optical thickness AOD (z), transmittance T i (z) and visibility V (z).
In an alternative embodiment of the present application, the atmospheric parameter measurement device further comprises a display, the display being electrically connected to the processor. The display may be used to display one or more of the measured atmospheric parameters.
In a second aspect, the present application provides an atmospheric parameter measurement method applied to the atmospheric parameter measurement device provided in the first aspect, including:
taking the distance from the laser radar in the height direction as a reflection height z j ；
At the same said reflection height z j In the case of (a), recording echo signals P of reflected beams of different angles received by the lidar ij The reflected beam forms an angle alpha with the height direction Z ij ；
For the echo signal P ij Processing to obtain the echo signal P ij Corresponding optical power signal RCS ij ；
At the same said reflection height z j In the case of calculating the optical power signal RCS ij And said angle alpha ij Slope τ of the relation of (2) 0j ；
Calculating each of the reflection heights z j Corresponding slope τ 0j ；
And establishing a slope function tau (z), and deriving the slope function tau (z) to obtain an atmospheric extinction coefficient R (z).
It can be appreciated that the present application discloses an atmospheric parameter measurement method, in which reflected beams of light at different angles are received at different heights by a lidar. By analysing the echo signals P of these reflected beams ij Is based on the optical power information of (a)Number RCS ij And angle alpha ij The current atmospheric extinction coefficient R (z) is calculated. Based on the atmospheric extinction coefficient R (z) and the optical power signal RCS ij The formula calculates other atmospheric parameters such as back scattering volume coefficient beta (z), optical thickness AOD (z), and transmittance T i (z) and/or visibility V (z). The method does not need to assume the laser radar ratio, and obtains more accurate atmospheric parameters through a multi-angle scanning method and an inversion method.
In an alternative embodiment of the present application, the reflection heights z are the same j In the case of calculating the optical power signal RCS ij And said angle alpha ij Slope τ of the relation of (2) 0j Comprising:
will be the same as the reflection height z j Each of said angles alpha in the case ij And corresponding said optical power signal RCS ij Substituting the formula to calculate the slope τ 0j ：
In(RCS i (z))＝In(C·β i (z))-2τ 0j ·x ij
Wherein x is ij ＝1/cosα ij Back scattering volume coefficient beta of atmosphere i (Z) is constant in the height direction Z.
In an alternative embodiment of the present application, the method further comprises:
substituting the slope function τ (z) into the following equation to obtain the backscattering volume coefficient β (z) of the atmosphere:
In(RCS(z))＝In(C·β(z))-2τ(z)。
in an alternative embodiment of the present application, the method further comprises:
the lidar ratio LR (z) of the lidar is calculated by:
LR(z)＝R(z)/β(z)。
in an alternative embodiment of the present application, the method further comprises:
the optical thickness of the atmosphere AOD (z) was calculated by:
wherein z is m For the reflection height z j Is a maximum value of (a).
In an alternative embodiment of the present application, the method further comprises:
the transmittance T of the atmosphere was calculated by i (z)：
In an alternative embodiment of the present application, the method further comprises:
atmospheric visibility V (z) was calculated by:
in an alternative embodiment of the present application, the method further comprises:
displaying, by a display, one or more of the measured atmospheric parameters; the atmospheric parameters comprise an atmospheric extinction coefficient R (z), a back scattering volume coefficient beta (z), an optical thickness AOD (z) and a transmittance T i (z) and visibility V (z).
The beneficial effects are that:
the application discloses an atmospheric parameter measuring device, through rotary device rotatory laser radar in order to receive the reflection light beam of different angles. The processor processes the echo signal data of each reflected light beam and calculates the atmospheric parameters of the current environment. The measuring device does not need to assume laser radar ratio, and provides a large number of reliable parameters for the subsequent inversion method through a multi-angle scanning method so as to calculate more accurate atmospheric parameters.
The application discloses an atmospheric parameter measurement method, which receives reflected light beams with different angles through a laser radar. By analysing the echo signals P of these reflected beams ij Is an optical power signal RCS of (a) ij And angle alpha ij The current atmospheric extinction coefficient R (z) is calculated. And then according to the atmospheric extinction coefficient R (z) and the optical powerRate signal RCS ij The formula calculates other atmospheric parameters such as back scattering volume coefficient beta (z), optical thickness AOD (z), and transmittance T i (z) and/or visibility V (z). The method does not need to assume the laser radar ratio, and obtains more accurate atmospheric parameters through a multi-angle scanning method and an inversion method.
In order to make the above objects, features and advantages of the present application more comprehensible, alternative embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the embodiments will be briefly described below, it being understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered limiting the scope, and that other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic diagram of the connection relationship of an atmospheric parameter measurement device provided in the present application;
FIG. 2 is a schematic view of the structure of the atmospheric parameter measuring apparatus shown in FIG. 1;
FIG. 3 is a broken-line cross-sectional view of the atmospheric parameter measuring device shown in FIG. 2;
FIG. 4 is a schematic flow chart of an atmospheric parameter measurement method provided in the present application;
fig. 5 is a schematic view of a measurement scenario of the atmospheric parameter measurement method shown in fig. 4.
Detailed Description
The following description of the embodiments of the present application will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are only some, but not all, of the embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art based on the embodiments herein without making any inventive effort, are intended to be within the scope of the present application.
In a first aspect, as shown in fig. 1, the present application provides an atmospheric parameter measurement device, which includes: a processor 10, a lidar 20, a rotation device 30.
The lidar 20 is arranged in a direction away from the ground. As illustrated in fig. 5, the lidar 20 is simplified to a point O.
The direction perpendicular to the ground is the height direction Z, the laser radar 20 is fixed on the rotating device 30, and the rotating device 30 is used for rotating the laser radar 20 to receive reflected light beams with different angles, and the reflected light beams form different angles with the height direction Z.
The processor 10 is electrically connected to a lidar 20.
It can be appreciated that the atmospheric parameter measuring device provided by the application rotates the laser radar through the rotating device to receive reflected light beams with different angles under the condition that the position of the laser radar is unchanged. The processor processes the echo signal data of each reflected light beam and calculates the atmospheric parameters of the current environment.
Wherein the rotation means 30 comprises a horizontal plane rotation means 31 and a vertical plane rotation means 32 as shown in fig. 1 and 2; the horizontal plane rotation device 31 is used for rotating the laser radar 20 around the height direction Z on a plane parallel to the ground so as to realize the rotation of the laser radar 20 on an XY plane; the vertical plane rotation device 32 is used for rotating the lidar 20 on a vertical plane perpendicular to the ground about a vertical axis of the vertical plane to effect rotation of the lidar 20 on a vertical plane such as the XZ plane.
In the embodiment of the present application, as shown in fig. 2, the laser radar 20 includes a radar main body 21, a laser transmitter 22, and a laser receiver 23, wherein the laser transmitter 22 and the laser receiver 23 are disposed on the same side of the radar main body 21.
In the embodiment of the present application, as shown in fig. 2 and 3, the vertical plane rotating device 32 includes a base 321, a support bar 322, and a vertical motor 323, a rotation groove 324 is provided on the base 321, a first end of the support bar 322 is inserted into the rotation groove 324, and a second end of the support bar 322 is fixed to the radar subject 21. The rotating rod of the vertical motor 323 passes through the through hole on the first end of the support bar 322 and passes from one side of the base 321 to the opposite side of the base 321. After the vertical motor 323 is started, the rotating rod of the vertical motor 323 drives the supporting rod 322 to rotate along the arrow 1 direction in fig. 2, so as to drive the laser radar 20 to rotate, and reflected light beams with different angles are received.
In the embodiment of the present application, as shown in fig. 2 and 3, the horizontal plane rotation device 31 includes a horizontal motor 310, and a rotation rod of the horizontal motor 310 passes through and is fixed in a groove of a surface of the base 321 facing away from the support bar 322. After the horizontal motor 310 is started, the rotating rod of the horizontal motor 310 drives the base 321 to rotate along the arrow 2 direction in fig. 2, so as to drive the laser radar 20 to rotate, so as to receive the reflected light beams with different angles.
The application discloses an atmospheric parameter measuring device, which rotates a laser radar 20 through a rotating device 30 to receive reflected light beams with different angles. The processor 10 processes the echo signal data of the reflected light beams at the respective reflection heights, and calculates the atmospheric parameters of the current environment. The measuring device does not need to assume laser radar ratio, and provides a large number of reliable parameters for the subsequent inversion method through a multi-angle scanning method so as to calculate more accurate atmospheric parameters.
In alternative embodiments of the present application, the atmospheric parameters include atmospheric extinction coefficient R (z), backscattering volume coefficient β (z), optical thickness AOD (z), transmittance T i (z) and visibility V (z).
The atmospheric extinction coefficient refers to the relative attenuation rate of electromagnetic radiation when the electromagnetic radiation propagates in the atmosphere for a unit distance.
The backscattering volume coefficient refers to the differential scattering cross section per unit volume and unit solid angle in the direction of 180 ° to the incident sound wave.
Wherein, the optical thickness refers to the integral of the extinction coefficient along the atmospheric transmission path, and is a dimensionless quantity representing the attenuation degree of the atmospheric medium to the radiation.
The transmittance refers to the ratio of the electromagnetic radiation flux attenuated by the atmosphere to the electromagnetic radiation flux when the electromagnetic wave propagates in the atmosphere. Atmospheric transmittance is an important factor affecting the transmission of infrared radiation. The air permeability of a specific area is greatly affected by the air condition, and the pressure, humidity and air density of the air can be obviously changed in a short time. The transmittance will thus vary to a large extent.
Visibility is an index reflecting transparency of the atmosphere. Is generally defined as the maximum ground level distance that a person with normal vision can also see the outline of the object under the prevailing weather conditions. There is also a minimum distance defined as the minimum distance at which the last features of the object have disappeared.
In an alternative embodiment of the present application, the atmospheric parameter measuring device further comprises a display 40, the display 40 being electrically connected to the processor 10. The display 40 may be used to display one or more of the measured atmospheric parameters.
In a second aspect, as shown in fig. 4, the present application provides an atmospheric parameter measurement method applied to the atmospheric parameter measurement device provided in the first aspect, including:
410. taking the distance from the laser radar in the height direction as a reflection height z j 。
420. At the same said reflection height z j In the case of (a), recording echo signals P of reflected beams of different angles received by the lidar ij The reflected beam forms an angle alpha with the height direction Z ij 。
As shown in FIG. 5, the lidar 20 is simplified to a point O, and the ground is set at Z-axis Z j Where it is located. The vertical plane rotation device 32 is used for rotating the lidar 20 on a vertical plane perpendicular to the ground about a vertical axis of the vertical plane to achieve a rotation of the lidar 20 at different angles αij from the height direction Z.
430. For echo signal P ij Processing to obtain echo signal P ij Corresponding optical power signal RCS ij 。
440. At the same reflection height z j In the case of calculating the optical power signal RCS ij And angle alpha ij Slope τ of the relation of (2) 0j 。
In an alternative embodiment of the present application, step 440 specifically includes:
will have the same reflection height z j Each angle alpha in the case ij And corresponding optical power signal RCS ij Substituting into the following formula, calculating to obtain slope τ 0j ：
In(RCS i (z))＝In(C·β i (z))-2τ 0j ·x ij
Wherein x is ij ＝1/cosα ij Back scattering volume coefficient beta of atmosphere i (Z) is constant in the height direction Z.
450. Calculating each reflection height z j Corresponding slope τ 0j 。
460. And establishing a slope function tau (z), and deriving the slope function tau (z) to obtain the atmospheric extinction coefficient R (z).
It will be appreciated that this application discloses an atmospheric parameter measurement method whereby reflected beams of light at different angles reflected by respective reflection heights are received by a lidar 20. By analysing the echo signals P of these reflected beams ij Is an optical power signal RCS of (a) ij And angle alpha ij The current atmospheric extinction coefficient R (z) is calculated. Based on the atmospheric extinction coefficient R (z) and the optical power signal RCS ij The formula calculates other atmospheric parameters such as back scattering volume coefficient beta (z), optical thickness AOD (z), and transmittance T i (z) and/or visibility V (z). The method does not need to assume the laser radar ratio, and obtains more accurate atmospheric parameters through a multi-angle scanning method and an inversion method.
The atmospheric extinction coefficient refers to the relative attenuation rate of electromagnetic radiation when the electromagnetic radiation propagates in the atmosphere for a unit distance.
In an alternative embodiment of the present application, the method further comprises:
substituting the slope function τ (z) into the following equation, the backscattering volume coefficient β (z) of the atmosphere is obtained:
In(RCS(z))＝In(C·β(z))-2τ(z)。
the backscattering volume coefficient refers to the differential scattering cross section per unit volume and unit solid angle in the direction of 180 ° to the incident sound wave.
In an alternative embodiment of the present application, the method further comprises:
the lidar ratio LR (z) of the lidar 20 is calculated by:
LR(z)＝R(z)/β(z)。
the method does not need to assume the laser radar ratio, and calculates the laser radar ratio LR (z) through a multi-angle scanning method and an inversion method.
In an alternative embodiment of the present application, the method further comprises:
the optical thickness of the atmosphere AOD (z) was calculated by:
wherein z is m For a reflection height z j Is a maximum value of (a).
Wherein, the optical thickness refers to the integral of the extinction coefficient along the atmospheric transmission path, and is a dimensionless quantity representing the attenuation degree of the atmospheric medium to the radiation.
In an alternative embodiment of the present application, the method further comprises:
the transmittance T of the atmosphere was calculated by i (z)：
The transmittance refers to the ratio of the electromagnetic radiation flux attenuated by the atmosphere to the electromagnetic radiation flux when the electromagnetic wave propagates in the atmosphere. Atmospheric transmittance is an important factor affecting the transmission of infrared radiation. The air permeability of a specific area is greatly affected by the air condition, and the pressure, humidity and air density of the air can be obviously changed in a short time. The transmittance will thus vary to a large extent.
In an alternative embodiment of the present application, the method further comprises:
atmospheric visibility V (z) was calculated by:
visibility is an index reflecting transparency of the atmosphere. Is generally defined as the maximum ground level distance that a person with normal vision can also see the outline of the object under the prevailing weather conditions. There is also a minimum distance defined as the minimum distance at which the last features of the object have disappeared.
In an alternative embodiment of the present application, the method further comprises:
displaying one or more of the measured atmospheric parameters via display 40; the atmospheric parameters include atmospheric extinction coefficient R (z), back scattering volume coefficient beta (z), optical thickness AOD (z), transmittance T i (z) and visibility V (z).
The application discloses an atmospheric parameter measurement method, which receives reflected light beams with different angles reflected by various reflection heights through a laser radar 20; by analysing the echo signals P of these reflected beams ij Is an optical power signal RCS of (a) ij And angle alpha ij Calculating the current atmospheric extinction coefficient R (z); based on the atmospheric extinction coefficient R (z) and the optical power signal RCS ij The formula calculates other atmospheric parameters such as back scattering volume coefficient beta (z), optical thickness AOD (z), and transmittance T i (z) and/or visibility V (z). The method does not need to assume the laser radar ratio, and obtains more accurate atmospheric parameters through a multi-angle scanning method and an inversion method.
All embodiments in the application are described in a progressive manner, and identical and similar parts of all embodiments are mutually referred, so that each embodiment mainly describes differences from other embodiments. In particular, for the apparatus, device and medium embodiments, since they are substantially similar to the method embodiments, the description is relatively simple, and the relevant parts will be referred to in the description of the method embodiments, which is not repeated herein.
Thus, particular embodiments of the present subject matter have been described. Other embodiments are within the scope of the following claims. In some cases, the actions recited in the claims can be performed in a different order and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may be advantageous.
The terms "first," "second," "the first," or "the second," as used in various embodiments of the present disclosure, may modify various components without regard to order and/or importance, but these terms do not limit the corresponding components. The above description is only configured for the purpose of distinguishing an element from other elements. For example, the first user device and the second user device represent different user devices, although both are user devices. For example, a first element could be termed a second element, and, similarly, a second element could be termed a first element, without departing from the scope of the present disclosure.
When an element (e.g., a first element) is referred to as being "coupled" (operatively or communicatively) to "another element (e.g., a second element) or" connected "to another element (e.g., a second element), it is understood that the one element is directly connected to the other element or the one element is indirectly connected to the other element via yet another element (e.g., a third element). In contrast, it will be understood that when an element (e.g., a first element) is referred to as being "directly connected" or "directly coupled" to another element (a second element), then no element (e.g., a third element) is interposed therebetween.
The above description is only illustrative of the principles of the technology being applied to alternative embodiments of the present application. It will be appreciated by persons skilled in the art that the scope of the invention referred to in this application is not limited to the specific combinations of features described above, but it is intended to cover other embodiments in which any combination of features described above or equivalents thereof is possible without departing from the spirit of the invention. Such as the above-described features and technical features having similar functions (but not limited to) disclosed in the present application are replaced with each other.
The foregoing is merely an alternative embodiment of the present application and is not intended to limit the present application, and various modifications and variations may be made to the present application by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.

Claims (3)

1. An atmospheric parameter measuring method is characterized in that,
taking the vertical distance between the laser radar and the laser radar in the height direction as a reflection height z j ；
At the same said reflection height z j In the case of (a), recording echo signals P of reflected beams of different angles received by the lidar ij The reflected beam forms an angle alpha with the height direction Z ij ；
For the echo signal P ij Processing to obtain the echo signal P ij Corresponding optical power signal RCS ij ；
Will be the same as the reflection height z j Each of said angles alpha in the case ij And corresponding said optical power signal RCS ij Substituting into the following formula, calculating to obtain slope τ 0j ：
In(RCS i (z))＝In(C·β i (z))-2τ 0j ·x ij
Wherein x is ij ＝1/cosα ij Back scattering volume coefficient beta of atmosphere i (Z) is constant in the height direction Z;
calculating each of the reflection heights z j Corresponding slope τ 0j ；
Establishing a slope function tau (z), and deriving the slope function tau (z) to obtain an atmospheric extinction coefficient R (z);
substituting the slope function τ (z) into the following equation to obtain the backscattering volume coefficient β (z) of the atmosphere:
In(RCS(z))＝In(C·β(z))-2τ(z)；
the optical thickness of the atmosphere AOD (z) was calculated by:
wherein z is m For the reflection height z j Is the maximum value of (2);
the transmittance T of the atmosphere was calculated by i (z)：
Atmospheric visibility V (z) was calculated by:
V(z)＝3/[R(z)(λ/550) 1.3 ]。
2. the method for measuring atmospheric parameters according to claim 1, wherein,
the method further comprises the steps of:
the lidar ratio LR (z) of the lidar is calculated by:
LR(z)＝R(z)/β(z)。
3. the method for measuring atmospheric parameters according to any one of claims 1 to 2, wherein,
the method further comprises the steps of:
displaying, by a display, one or more of the measured atmospheric parameters; the atmospheric parameters comprise an atmospheric extinction coefficient R (z), a back scattering volume coefficient beta (z), an optical thickness AOD (z) and a transmittance T i (z) and visibility V (z).
CN202110707189.8A 2021-06-24 2021-06-24 Atmospheric parameter measurement method Active CN113433569B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110707189.8A CN113433569B (en) 2021-06-24 2021-06-24 Atmospheric parameter measurement method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110707189.8A CN113433569B (en) 2021-06-24 2021-06-24 Atmospheric parameter measurement method

Publications (2)

Publication Number Publication Date
CN113433569A CN113433569A (en) 2021-09-24
CN113433569B true CN113433569B (en) 2024-03-01

Family

ID=77754256

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110707189.8A Active CN113433569B (en) 2021-06-24 2021-06-24 Atmospheric parameter measurement method

Country Status (1)

Country Link
CN (1) CN113433569B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107421917A (en) * 2017-05-17 2017-12-01 南京信息工程大学 A kind of multifunction high-precision atmosphere visibility meter and visibility measurement method
CN107807366A (en) * 2017-10-30 2018-03-16 中国科学技术大学 A kind of computational methods of atmospheric visibility, device, radar and system
CN108490451A (en) * 2018-03-29 2018-09-04 中国民航大学 A method of utilizing atmospheric extinction coefficient inverting slant visibility
CN108627812A (en) * 2018-05-28 2018-10-09 成都信息工程大学 A kind of laser radar atmospheric visibility measurement method and device
WO2019096878A1 (en) * 2017-11-17 2019-05-23 Leosphere Device and method for calibrating a lidar
CN109827906A (en) * 2019-01-31 2019-05-31 西安理工大学 A kind of inversion method of laser radar slant visibility
CN110333518A (en) * 2019-08-06 2019-10-15 大舜激光(黄山)科技有限公司 Method, system and the laser radar of correlation measurement extinction coefficient
CN111736176A (en) * 2020-06-23 2020-10-02 中国人民解放军陆军炮兵防空兵学院 Atmospheric environment data sharing method based on laser radar

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107421917A (en) * 2017-05-17 2017-12-01 南京信息工程大学 A kind of multifunction high-precision atmosphere visibility meter and visibility measurement method
CN107807366A (en) * 2017-10-30 2018-03-16 中国科学技术大学 A kind of computational methods of atmospheric visibility, device, radar and system
WO2019096878A1 (en) * 2017-11-17 2019-05-23 Leosphere Device and method for calibrating a lidar
CN108490451A (en) * 2018-03-29 2018-09-04 中国民航大学 A method of utilizing atmospheric extinction coefficient inverting slant visibility
CN108627812A (en) * 2018-05-28 2018-10-09 成都信息工程大学 A kind of laser radar atmospheric visibility measurement method and device
CN109827906A (en) * 2019-01-31 2019-05-31 西安理工大学 A kind of inversion method of laser radar slant visibility
CN110333518A (en) * 2019-08-06 2019-10-15 大舜激光(黄山)科技有限公司 Method, system and the laser radar of correlation measurement extinction coefficient
CN111736176A (en) * 2020-06-23 2020-10-02 中国人民解放军陆军炮兵防空兵学院 Atmospheric environment data sharing method based on laser radar

Also Published As

Publication number Publication date
CN113433569A (en) 2021-09-24

Similar Documents

Publication Publication Date Title
CN102879309B (en) Gas particle concentration measurement method and device on basis of broadband linear frequency modulation ultrasound
CN110006848B (en) Method and device for obtaining extinction coefficient of aerosol
CN105911060B (en) A kind of visible detection method and device of the pollution of transmission-type visual range visibility meter window mirror
CN106706566B (en) A kind of calculation method of laser radar detection SEQUENCING VERTICAL visibility
CN109001217A (en) A kind of detection method and detection device of phone housing defect
CN110987872B (en) Forward scatter sensor
US4754151A (en) Method and apparatus for measuring the optical scattering co-efficient of the ocean using large and small fields of view
CN113433569B (en) Atmospheric parameter measurement method
CN106790787A (en) A kind of method of digital product and Detection of Air Quality
Phillips et al. Study of the spectral variation of absorption and scattering in some Australian coastal waters.
US5712709A (en) Haze and transmissivity measurements
JP3504163B2 (en) Snow depth gauge
CN108414464A (en) Water body multi-wavelength optical attenuation coefficient measuring device and method
CN110231374A (en) A kind of greasy dirt on-line monitoring method based on hydrophobic oil absorption material
Pedersen et al. Ultrasound properties of lung tissue and their measurements
CN207730933U (en) A kind of mobile orientable laser radar network data quality control system
CN202814843U (en) Gas particle concentration measurement device based on broadband linear frequency modulation ultrasounds
Guenther Wind And Nadir Angle Effects On Airborne Lidar Water" Surface" Returns
CN104849724A (en) Measurement method and apparatus for laser radar ratio of aerosol
CN112525778B (en) Suspended load particle size and grading measurement method based on novel multi-frequency ultrasonic probe
CN108344712A (en) A kind of measuring device and its measurement method of Refractive Index of Material
CN110987873B (en) Forward scatter sensor
CN209117031U (en) A kind of verticality detector
Honey et al. Optical absorption and turbulence-induced narrow-angle forward scatter in the sea
RU2167434C2 (en) Method determining density of concentration of fish

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