CN112597426B - Calculation method, device, equipment and storage medium for optical thickness of night aerosol - Google Patents
Calculation method, device, equipment and storage medium for optical thickness of night aerosol Download PDFInfo
- Publication number
- CN112597426B CN112597426B CN202110237755.3A CN202110237755A CN112597426B CN 112597426 B CN112597426 B CN 112597426B CN 202110237755 A CN202110237755 A CN 202110237755A CN 112597426 B CN112597426 B CN 112597426B
- Authority
- CN
- China
- Prior art keywords
- optical thickness
- moon
- calculating
- aerosol
- equivalent
- 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
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 145
- 239000000443 aerosol Substances 0.000 title claims abstract description 77
- 238000004364 calculation method Methods 0.000 title claims abstract description 25
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims abstract description 95
- 239000007789 gas Substances 0.000 claims abstract description 88
- 238000002310 reflectometry Methods 0.000 claims abstract description 38
- 239000007787 solid Substances 0.000 claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000002834 transmittance Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 9
- 238000004164 analytical calibration Methods 0.000 claims description 8
- 238000004891 communication Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- JCXJVPUVTGWSNB-UHFFFAOYSA-N nitrogen dioxide Inorganic materials O=[N]=O JCXJVPUVTGWSNB-UHFFFAOYSA-N 0.000 description 3
- 230000008033 biological extinction Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 229910002092 carbon dioxide Inorganic materials 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 1
- 206010011655 Cushingoid Diseases 0.000 description 1
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
- G01N15/075—Investigating concentration of particle suspensions by optical means
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Dispersion Chemistry (AREA)
- Mathematical Optimization (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
The invention provides a method, a device, equipment and a storage medium for calculating optical thickness of night aerosol, wherein the method comprises the following steps: obtaining the equivalent lunar disc reflectivity; calculating a solid angle of the observation time according to the sun-moon distance and the moon-earth distance; calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance; obtaining equivalent gas optical thickness and Rayleigh gas optical thickness of a plurality of preset gases; and obtaining the optical thickness of the aerosol according to the atmospheric layer top moon irradiance, the equivalent gas optical thickness and the Rayleigh gas optical thickness at the observation time. The method has high calculation accuracy on the optical thickness of the aerosol at night.
Description
Technical Field
The invention relates to the technical field of meteorological observation, in particular to a method, a device, equipment and a storage medium for calculating optical thickness of aerosol at night.
Background
The aerosol has a crucial influence on the radiation energy balance of the earth's land and the atmosphere circle, and is the largest uncertain factor in the global climate change assessment at present. In addition to model simulations using theoretical models, a more important approach to research is the extensive observation of aerosols. In Aerosol ground-based remote sensing and satellite remote sensing observation, the Optical Aerosol Depth (AOD) which represents the extinction characteristic of the whole Aerosol layer is the most basic parameter of Aerosol remote sensing observation. Currently, aerosol optical thickness parameters are commonly observed and acquired during the daytime, and observed and studied rarely at night. This is mainly because the optical thickness of the aerosol is inverted by the spectral attenuation of the sun reflected by the earth surface, and the ground is obtained by directly using the direct sun light. They have the common feature that both the inversion algorithm and the observation instrument depend on sunlight with strong energy. For weak light sources at night, such as moonlight or starlight, professional and practical instruments and equipment are not available for automatically observing weak light, and the weak light sources are unstable compared with sunlight, so that the corresponding aerosol optical thickness algorithm is immature.
In recent years, night low-light remote sensing is gradually developed. In the nineties of the last century, researchers were first attempting to obtain nighttime low light observations using a starlight meter. The use of starlight to measure optical thickness of an aerosol at night is generally better than using moonlight. However, since the starlight is too weak, a starlight meter with a large field angle is required for observation, and portability of the instrument is a problem. From the perspective of instrument portability, it is proposed to observe moonlight by using a skylight channel of a CE318 sun photometer widely used at present, and invert the optical thickness of the aerosol at night according to an automatic moon viewing platform ROLO (robotic Lunar observer) model of the U.S. geological survey bureau, however, the optical thickness of the aerosol at night calculated by applying the ROLO model to the moon photometer is not accurate.
The threshold value of the exposure intensity of the moon light of the formal version CE318T of the moon photometer is 1/4 phases of the moon, namely, the lunar calendar is observed in the first three to twenty five nights of each month, and is not observed in the first two of twenty six to next months. The model machine is used for testing on a Spanish island, testing the local atmospheric condition and the performance of an instrument, developing a Moon irradiance calculation model RIMO (ROLO approximation for Moon's emission) suitable for Spanish local atmospheric Observation on the basis of an ROLO model, and calculating the definite night aerosol optical thickness of Spanish based on the model. However, the RIMO model has a large error, which can reach 14% when the lunar phase is small.
Disclosure of Invention
The invention provides a method, a device, equipment and a storage medium for calculating optical thickness of night aerosol, which are used for solving the defects that an ROLO model in the prior art is not accurate when applied to a moon photometer to calculate the optical thickness of the night aerosol and an RIMO model has larger error when a moon phase is smaller, and realize high accuracy of calculation of the optical thickness of the night aerosol.
The invention provides a method for calculating optical thickness of aerosol at night, which comprises the following steps: obtaining the equivalent lunar disc reflectivity;
calculating a solid angle of the observation time according to the sun-moon distance and the moon-earth distance;
calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance;
obtaining equivalent gas optical thickness and Rayleigh gas optical thickness of a plurality of preset gases;
and obtaining the optical thickness of the aerosol according to the atmospheric layer top moon irradiance, the equivalent gas optical thickness and the Rayleigh gas optical thickness at the observation time.
The invention provides a night aerosol optical thickness calculation method, which comprises the following steps of: and calculating the equivalent moon disc reflectivity by adopting an automatic moon viewing platform ROLO model.
According to the method for calculating the optical thickness of the aerosol at night, provided by the invention, the calculating of the equivalent lunar disc reflectivity by using the ROLO model comprises the following steps: and calculating the equivalent lunar disc reflectivity according to the fixed empirical coefficient, the moon absolute phase angle, the longitude and latitude of the observation point on the lunar surface and the longitude of the sun on the lunar surface.
According to the invention, the method for calculating the optical thickness of the night aerosol, which calculates the solid angle of the lunar observation time according to the sun-moon distance and the lunar-earth distance, comprises the following steps: and calculating the solid angle of the moon observation time according to the average moon solid angle, the distance from the sun to the moon at the observation time and the distance from the observation point to the moon.
According to the method for calculating the optical thickness of the aerosol at night, before the step of calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance, the method further comprises the following steps: and carrying out convolution on the filter transmittance profile of the moon photometer and the solar model to obtain the solar irradiance.
According to the night aerosol optical thickness calculation method provided by the invention, the obtaining of the equivalent gas optical thickness of a plurality of preset gases comprises the following steps: simulating the transmittance of the various preset gases in a moon photometer channel; and calculating the equivalent gas optical thickness of the plurality of preset gases according to the transmittances of the plurality of preset gases in the channel of the moon photometer.
According to the method for calculating the optical thickness of the aerosol at night provided by the invention, the optical thickness of the aerosol is obtained according to the irradiance of the moon at the top of the atmospheric layer at the observation time, the optical thickness of the equivalent gas and the optical thickness of the Rayleigh gas, and the method comprises the following steps: acquiring a pixel brightness value and an instrument calibration coefficient of a remote sensing image observed by a moon photometer; obtaining the total optical thickness of the atmosphere according to the moon irradiance at the top of the atmosphere layer at the observation time, the pixel brightness value of the remote sensing image observed by the moon photometer and the instrument calibration coefficient; and obtaining the optical thickness of the aerosol according to the total optical thickness of the atmosphere, the optical thickness of the equivalent gas and the optical thickness of the Rayleigh gas.
The invention also provides a night aerosol optical thickness calculating device, comprising: the first acquisition module is used for acquiring the equivalent lunar disc reflectivity;
the first calculation module is used for calculating a solid angle of an observation moment according to the sun-moon distance and the moon-earth distance;
the second calculation module is used for calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance;
the second acquisition module is used for acquiring equivalent gas optical thickness and Rayleigh gas optical thickness of multiple preset gases;
and the third acquisition module is used for obtaining the optical thickness of the aerosol according to the atmospheric layer top moon irradiance, the equivalent gas optical thickness and the Rayleigh gas optical thickness at the observation moment.
The invention also provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor executes the program to realize the steps of the night aerosol optical thickness calculation method.
The present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the method for calculating optical thickness of an aerosol at night as described in any one of the above.
The method, the device, the equipment and the storage medium for calculating the optical thickness of the night aerosol provided by the invention can be used for obtaining the reflectivity of the equivalent lunar disc and the solid angle of the observation time, further calculating the top lunar irradiance of the atmosphere layer at the observation time according to the reflectivity of the equivalent lunar disc, the solid angle of the observation time and the solar irradiance, then obtaining the optical thickness of the aerosol according to the top lunar irradiance of the atmosphere layer at the observation time, the optical thickness of the equivalent gas and the optical thickness of the Rayleigh gas, and then accurately calculating the optical thickness of the night atmosphere aerosol under the conditions of no cloud and weak light at night based on the observation data of the direct light of the moon according to the atmosphere extinction law.
Drawings
In order to more clearly illustrate the technical solutions of the present invention or the prior art, the drawings needed for the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
FIG. 1 is a schematic flow chart of a method for calculating optical thickness of night aerosol according to the present invention;
fig. 2 is a block diagram of a device for calculating optical thickness of night aerosol according to the present invention.
Fig. 3 is a schematic physical structure diagram of an electronic device in an example of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The night aerosol optical thickness calculation method of the present invention is described below with reference to fig. 1.
Fig. 1 is a schematic flow chart of a method for calculating optical thickness of night aerosol according to the present invention. As shown in fig. 1, the method for calculating the optical thickness of the aerosol at night according to the embodiment of the present invention includes:
s1: and obtaining the equivalent lunar disc reflectivity.
In one embodiment of the invention, the equivalent lunar disc reflectivity calculated using the ROLO modelA k Calculated according to the following formula:
wherein the content of the first and second substances,a i 、b j 、c 1 、c 2 、c 3 、c 4 、d 1 、d 2 、d 3 、p 1 、p 2 、p 3andp 4in order to fix the empirical factor,gis the absolute phase angle of the moon and,θandϕto observeThe latitude and longitude of the point at the moon's face,Φlongitude of the sun in the moon. Calculating parameters from an astronomical calendar DE421g，θ，φ，Φ. The calculated equivalent lunar disc reflectivity wave band ranges from 350nm to 2450nm and has 32 wave bands, and then linear interpolation is carried out according to the adjacent left and right wave bands closest to the moon photometer channels (440 nm, 500nm, 675nm, 870nm and 1020 nm) to obtain the equivalent lunar disc reflectivity of the moon photometer channels.
S2: and calculating a solid angle of the observation time according to the sun-moon distance and the moon-earth distance.
In the present embodiment, the solid angle at the observation time of the moon is calculated from the average moon solid angle, the distance from the sun to the moon at the observation time, and the distance from the observation point to the moon. In particular, the lunar observation cube cornerΩCalculated according to the following formula:
wherein the content of the first and second substances,Ω 0 the mean moon solid angle is 6.4177X 10-5sr, D S-M The unit is astronomical unit AU for the distance from the sun to the moon at the observation time; D O-M the distance from the observation point to the moon is given in km. Calculating the above-mentioned astronomical parameters according to the astronomical calendar DE421D S-M ，D O-M 。
S3: and calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance.
Specifically, the atmospheric layer top moon irradiance model is an important parameter for calculating the optical thickness of the aerosol at night, and is related to the reflectivity of an equivalent moon disc, a moon observation solid angle and solar irradiance. The solar irradiance refers to the radiant energy of solar radiation reaching the surface of the solid earth in unit area and unit time after absorption, scattering, reflection and the like of the atmosphere.
In the embodiment, the irradiance of the top moon of the atmospheric layer at the time t is observedE 0,t,Moon Calculated according to the following formula:
wherein the content of the first and second substances,E 0,Sun the solar irradiance after the convolution of the instrument,A k equivalent lunar disc reflectivity. Specifically, the moon photometer channels were scanned in the laboratory using monochromatic light sources to obtain the transmittance profiles of the instrument channels (440 nm, 500nm, 675nm, 870nm, and 1020 nm)T(the ideal profile shape should be gaussian). Each channel profile is then compared to the solar irradiance in the Wehrli et al (1985) modelE 0,Sun,W And (3) performing convolution on the curve to obtain the solar irradiance after the convolution of the instrument:
s4: and acquiring equivalent gas optical thickness and Rayleigh gas optical thickness of multiple preset gases.
In this embodiment, the method for obtaining the equivalent gas optical thickness of the plurality of preset gases includes:
the transmittance of various preset gases in a moon photometer channel is simulated. For example, ozone O can be simulated by MODTRAN software3Methane CH4Carbon dioxide CO2Nitrogen dioxide NO2The transmission of various gases in the moon photometer channel.
According to the transmittance of various preset gases in a moon photometer channelT g Calculating equivalent gas optical thickness of multiple preset gases
：
In addition, it also providesOther bands need to be filtered out through the passage of the moon photometer to calculate the Rayleigh gas optical thickness through the specified band
：
Wherein
Is the wavelength (unit: micron), p is the pressure (unit: hectopascal) at the observation site, and H is the elevation (unit: meter) at the measurement site.
S5: and obtaining the optical thickness of the aerosol according to the moon irradiance at the top of the atmospheric layer, the equivalent optical thickness of the gas and the optical thickness of the Rayleigh gas at the observation time.
In particular, the aerosol optical thickness according to the law of atmospheric extinctionτ a Calculated by the following formula:
wherein the content of the first and second substances,
which represents the total optical thickness of the atmosphere,V t,Moon the pixel brightness value of remote sensing image of DN (digital number) is observed by a moon photometer;C Moon is an instrument calibration coefficient which is provided by the instrument after delivery or calibration,m(θ) Is the mass number of the moon penetrating the atmosphere.
The aerosol optical thickness calculating device provided by the invention is described below, and the aerosol optical thickness calculating device described below and the aerosol optical thickness calculating method described above can be correspondingly referred to each other.
Fig. 2 is a block diagram of a device for calculating optical thickness of night aerosol according to the present invention. As shown in fig. 2, the present invention provides an aerosol optical thickness calculating device, including: a first obtaining module 210, a first calculating module 220, a second calculating module 230, a second obtaining module 240, and a third obtaining module 250.
The first obtaining module 210 is configured to obtain an equivalent lunar disc reflectivity. The first calculation module 220 is used for calculating the solid angle of the observation time according to the sun-moon distance and the lunar-earth distance. The second calculation module 230 is used for calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance. The second obtaining module 240 is configured to obtain equivalent gas optical thicknesses and rayleigh gas optical thicknesses of a plurality of preset gases. The third obtaining module 250 is configured to obtain the optical thickness of the aerosol according to the atmospheric layer top moon irradiance, the equivalent gas optical thickness, and the rayleigh gas optical thickness at the observation time.
In one embodiment of the present invention, the first obtaining module 210 is configured to calculate an equivalent lunar disc reflectivity using an automatic lunar viewing station ROLO model.
Further, the first obtaining module 210 is specifically configured to calculate an equivalent lunar disc reflectivity according to a fixed empirical coefficient, an absolute phase angle of the moon, a longitude and a latitude of an observation point on the lunar surface, and a longitude of the sun on the lunar surface.
In one embodiment of the invention, the first calculation module 220 is configured to calculate the solid angle at the observation time of the moon from the mean moon solid angle, the distance from the sun to the moon at the observation time and the distance from the observation point to the moon.
In one embodiment of the present invention, the second calculation module 230 is further configured to obtain the solar irradiance by convolving the filter transmittance profile of the moon photometer with a solar model.
In an embodiment of the invention, the second obtaining module 240 is configured to simulate transmittances of a plurality of predetermined gases in a moonlight meter channel, and further calculate equivalent gas optical thicknesses of the plurality of predetermined gases according to the transmittances of the plurality of predetermined gases in the moonlight meter channel.
In an embodiment of the present invention, the third obtaining module 250 is configured to obtain a pixel brightness value and an instrument calibration coefficient of a remote sensing image observed by a moonlight meter, obtain an atmospheric total optical thickness according to an atmospheric layer top moonirradiance at an observation time, and obtain an aerosol optical thickness according to the atmospheric total optical thickness, an equivalent gas optical thickness, and a rayleigh gas optical thickness.
It should be noted that, a specific implementation of the night aerosol optical thickness calculation apparatus according to the embodiment of the present invention is similar to a specific implementation of the night aerosol optical thickness calculation method according to the embodiment of the present invention, and specific reference is specifically made to the description of the night aerosol optical thickness calculation method, and details are not repeated for reducing redundancy.
In addition, other configurations and functions of the night aerosol optical thickness calculating device according to the embodiment of the present invention are known to those skilled in the art, and are not described in detail in order to reduce redundancy.
Fig. 3 illustrates a physical structure diagram of an electronic device, which may include, as shown in fig. 3: a processor (processor)310, a communication Interface 820, a memory (memory)330 and a communication bus 340, wherein the processor 310, the communication Interface 320 and the memory 330 communicate with each other via the communication bus 340. Processor 310 may invoke logic instructions in memory 330 to perform a nighttime aerosol optical thickness calculation method comprising: obtaining the equivalent lunar disc reflectivity; calculating a solid angle of the observation time according to the sun-moon distance and the moon-earth distance; calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance; obtaining equivalent gas optical thickness and Rayleigh gas optical thickness of a plurality of preset gases; and obtaining the optical thickness of the aerosol according to the atmospheric layer top moon irradiance, the equivalent gas optical thickness and the Rayleigh gas optical thickness at the observation time.
In addition, the logic instructions in the memory 330 may be implemented in the form of software functional units and stored in a computer readable storage medium when the software functional units are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, the present invention also provides a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions which, when executed by a computer, enable the computer to perform the method for calculating optical thickness of night time aerosol provided by the above methods, the method comprising: obtaining the equivalent lunar disc reflectivity; calculating a solid angle of the observation time according to the sun-moon distance and the moon-earth distance; calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance; obtaining equivalent gas optical thickness and Rayleigh gas optical thickness of a plurality of preset gases; and obtaining the optical thickness of the aerosol according to the atmospheric layer top moon irradiance, the equivalent gas optical thickness and the Rayleigh gas optical thickness at the observation time.
In yet another aspect, the present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program that when executed by a processor is implemented to perform the nighttime aerosol optical thickness calculation method provided above, the method comprising: obtaining the equivalent lunar disc reflectivity; calculating a solid angle of the observation time according to the sun-moon distance and the moon-earth distance; calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance; obtaining equivalent gas optical thickness and Rayleigh gas optical thickness of a plurality of preset gases; and obtaining the optical thickness of the aerosol according to the atmospheric layer top moon irradiance, the equivalent gas optical thickness and the Rayleigh gas optical thickness at the observation time.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (9)
1. A method for calculating optical thickness of aerosol at night comprises the following steps:
obtaining the equivalent lunar disc reflectivity;
calculating a solid angle of the observation time according to the sun-moon distance and the moon-earth distance;
calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance;
obtaining equivalent gas optical thickness and Rayleigh gas optical thickness of a plurality of preset gases;
obtaining aerosol optical thickness according to the atmospheric layer top moon irradiance, the equivalent gas optical thickness and the rayleigh gas optical thickness at the observation time, wherein obtaining aerosol optical thickness according to the atmospheric layer top moon irradiance, the equivalent gas optical thickness and the rayleigh gas optical thickness at the observation time comprises: acquiring a pixel brightness value and an instrument calibration coefficient of a remote sensing image observed by a moon photometer; obtaining the total optical thickness of the atmosphere according to the moon irradiance at the top of the atmosphere layer at the observation time, the pixel brightness value of the remote sensing image observed by the moon photometer and the instrument calibration coefficient; and obtaining the optical thickness of the aerosol according to the total optical thickness of the atmosphere, the optical thickness of the equivalent gas and the optical thickness of the Rayleigh gas.
2. The night aerosol optical thickness calculation method of claim 1, wherein the obtaining an equivalent lunar disc reflectivity comprises:
and calculating the equivalent moon disc reflectivity by adopting an automatic moon viewing platform ROLO model.
3. The method of claim 2, wherein the calculating the equivalent lunar disc reflectivity using the ROLO model comprises:
and calculating the equivalent lunar disc reflectivity according to the fixed empirical coefficient, the moon absolute phase angle, the longitude and latitude of the observation point on the lunar surface and the longitude of the sun on the lunar surface.
4. The method of calculating optical thickness of an aerosol at night according to claim 1, wherein calculating the solid angle of the observation time from the distance between the sun and the moon and the distance between the moon and the earth comprises:
and calculating the solid angle of the moon observation time according to the average moon solid angle, the distance from the sun to the moon at the observation time and the distance from the observation point to the moon.
5. The method of calculating optical thickness of an aerosol at night according to claim 1, further comprising, prior to the calculating atmospheric layer top moon irradiance at a time of observation from the equivalent lunar disc reflectivity, the solid angle at the time of observation, and solar irradiance:
and carrying out convolution on the filter transmittance profile of the moon photometer and the solar model to obtain the solar irradiance.
6. The method for calculating optical thickness of night aerosol according to claim 1, wherein the obtaining the equivalent gas optical thickness of a plurality of preset gases comprises:
simulating the transmittance of the various preset gases in a moon photometer channel;
and calculating the equivalent gas optical thickness of the plurality of preset gases according to the transmittances of the plurality of preset gases in the channel of the moon photometer.
7. An optical thickness calculation device for an aerosol at night, comprising:
the first acquisition module is used for acquiring the equivalent lunar disc reflectivity;
the first calculation module is used for calculating a solid angle of an observation moment according to the sun-moon distance and the moon-earth distance;
the second calculation module is used for calculating the atmospheric layer top moon irradiance at the observation time according to the equivalent lunar disc reflectivity, the solid angle at the observation time and the solar irradiance;
the second acquisition module is used for acquiring equivalent gas optical thickness and Rayleigh gas optical thickness of multiple preset gases;
the third acquisition module is used for obtaining the optical thickness of the aerosol according to the atmospheric layer top moon irradiance, the equivalent gas optical thickness and the Rayleigh gas optical thickness at the observation moment; the third acquisition module is used for acquiring pixel brightness values and instrument calibration coefficients of remote sensing images observed by the moon photometer, and obtaining the total atmospheric optical thickness according to the atmospheric layer top moon irradiance at the observation time, the pixel brightness values of the remote sensing images observed by the moon photometer and the instrument calibration coefficients; and obtaining the optical thickness of the aerosol according to the total optical thickness of the atmosphere, the optical thickness of the equivalent gas and the optical thickness of the Rayleigh gas.
8. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program implements the steps of the method of calculating an optical thickness of a night aerosol as claimed in any one of claims 1 to 6.
9. A non-transitory computer readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the steps of the nighttime aerosol optical thickness calculation method according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110237755.3A CN112597426B (en) | 2021-03-04 | 2021-03-04 | Calculation method, device, equipment and storage medium for optical thickness of night aerosol |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110237755.3A CN112597426B (en) | 2021-03-04 | 2021-03-04 | Calculation method, device, equipment and storage medium for optical thickness of night aerosol |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112597426A CN112597426A (en) | 2021-04-02 |
| CN112597426B true CN112597426B (en) | 2021-05-14 |
Family
ID=75210347
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110237755.3A Active CN112597426B (en) | 2021-03-04 | 2021-03-04 | Calculation method, device, equipment and storage medium for optical thickness of night aerosol |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112597426B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117935063A (en) * | 2024-03-22 | 2024-04-26 | 航天宏图信息技术股份有限公司 | Method, device and equipment for inverting optical thickness of aerosol based on RSSDM nights |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105674901A (en) * | 2016-04-08 | 2016-06-15 | 电子科技大学 | Statistical segmentation based inversion method for atmospheric aerosol |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102636143B (en) * | 2012-04-20 | 2014-11-19 | 中国科学院遥感应用研究所 | Aerosol optical depth remote sensing retrieval method |
-
2021
- 2021-03-04 CN CN202110237755.3A patent/CN112597426B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105674901A (en) * | 2016-04-08 | 2016-06-15 | 电子科技大学 | Statistical segmentation based inversion method for atmospheric aerosol |
Non-Patent Citations (4)
| Title |
|---|
| 《A new Method for nocturnal aerosolmeasurements with a lunar photometerprototype》;A.Barreto1,E.Cuevas,B.Damiri,C.Guirado,T.Berkoff;《Atmos. Meas. Tech》;20130503;585-598 * |
| 《Correction of a lunar irradiance model for aerosol optical depthretrieval and comparison with star photometer》;Roman,Roberto;Gonzalez,Ramiro;Toledano,Carlos;《Atmos. Meas. Tech.》;20201124;第13卷(第11期);6293-6307 * |
| 《Nocturnal aerosol optical depth measurements with modified sky radiometer POM-02 using the moon as a light source》;A Uchiyama，M Shiobara，H Kobayashi，T Matsunaga，Y Watanabe;《Atmos. Meas. Tech.》;20191209;第12卷(第12期);6465-6488 * |
| 《基于微光数据的夜间亮目标表观反射率特性和城市光学厚度反演研究》;王锋;《中国优秀硕士学位论文全文数据库基础科学辑》;20200215(第2期);A009-116 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112597426A (en) | 2021-04-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Van Roozendael et al. | Ten years of GOME/ERS‐2 total ozone data—The new GOME data processor (GDP) version 4: 1. Algorithm description | |
| Bertaux et al. | TAPAS, a web-based service of atmospheric transmission computation for astronomy | |
| Gawiser et al. | The multiwavelength survey by Yale-Chile (MUSYC): Survey design and deep public UBVRIZ'images and catalogs of the extended Hubble Deep Field-South | |
| Gordon | Atmospheric correction of ocean color imagery in the Earth Observing System era | |
| Ichoku et al. | Analysis of the performance characteristics of the five‐channel Microtops II Sun photometer for measuring aerosol optical thickness and precipitable water vapor | |
| Sayer et al. | Extending “Deep Blue” aerosol retrieval coverage to cases of absorbing aerosols above clouds: Sensitivity analysis and first case studies | |
| Buton et al. | Atmospheric extinction properties above Mauna Kea from the Nearby SuperNova Factory spectro-photometric data set | |
| Kylling et al. | Effect of aerosols on solar UV irradiances during the Photochemical Activity and Solar Ultraviolet Radiation campaign | |
| Zege et al. | Algorithm to retrieve the melt pond fraction and the spectral albedo of Arctic summer ice from satellite optical data | |
| Strow et al. | Validation of the Atmospheric Infrared Sounder radiative transfer algorithm | |
| Salerno et al. | Novel retrieval of volcanic SO2 abundance from ultraviolet spectra | |
| Hulley et al. | Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature and emissivity product (MxD21) algorithm theoretical basis document collection-6 | |
| Frouin et al. | Estimating photosynthetically available radiation at the ocean surface from GOCI data | |
| Colarco et al. | Simulation of the Ozone Monitoring Instrument aerosol index using the NASA Goddard Earth Observing System aerosol reanalysis products | |
| Anagnostou et al. | Rainfall estimation from TOGA radar observations during LBA field campaign | |
| Harrington et al. | Daytime sky polarization calibration limitations | |
| Khosravi et al. | Retrieval of terrestrial plant fluorescence based on the in-filling of far-red fraunhofer lines using SCIAMACHY observations | |
| Coddington et al. | Aircraft measurements of spectral surface albedo and its consistency with ground‐based and space‐borne observations | |
| CN112597426B (en) | Calculation method, device, equipment and storage medium for optical thickness of night aerosol | |
| Dahl et al. | Vertical structure and color of Jovian latitudinal cloud bands during the Juno era | |
| Landier et al. | Calibration of urban canopies albedo and 3D shortwave radiative budget using remote-sensing data and the DART model | |
| Popp | Correcting atmospheric masking to retrieve the spectral albedo of land surfaces from satellite measurements | |
| Emili et al. | Comparison between the assimilation of IASI Level 2 ozone retrievals and Level 1 radiances in a chemical transport model | |
| Guérin et al. | Revisiting the Cox and Munk wave-slope statistics using IASI observations of the sea surface | |
| Lazzari et al. | Assessment of the spectral downward irradiance at the surface of the Mediterranean Sea using the radiative Ocean-Atmosphere Spectral Irradiance Model (OASIM) |
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 |