Satellite-borne laser radar and reanalysis data multi-wavelength atmospheric transmittance inversion method

Abstract

The invention provides a satellite-borne laser radar and reanalysis data multi-wavelength atmospheric transmittance inversion method, and relates to the field of atmospheric remote sensing. According to the inversion method, on the basis of satellite-borne laser radar data and reanalysis data, the atmospheric transmissivity under different wavelengths is accurately calculated in combination with wavelength conversion and data fusion; comprising the following implementation steps: step 1), selecting satellite-borne laser radar data; step 2), screening the extinction coefficients; step 3), generating a space-time matching data set; step 4), establishing an aerosol extinction-transmittance conversion model; step 5), calculating atmospheric molecule and ozone transmittance; step 6), performing layered transmittance fusion; and 7) calculating the transmittance of the whole atmosphere.

Description

TECHNICAL FIELD

[0001] The application provides a multi-wavelength atmospheric transmittance inversion method based on a spaceborne lidar and reanalysis data cooperation, and relates to the field of atmospheric remote sensing. BACKGROUND

[0002] As a core technology of modern remote sensing, spaceborne lidar has been widely used in global cloud-aerosol vertical distribution detection, carbon cycle monitoring, three-dimensional wind field inversion and other fields by transmitting high-precision laser pulses through a satellite platform and receiving atmospheric backscattering signals. However, the laser signal will be attenuated by the atmospheric transmittance during transmission, resulting in signal strength attenuation with transmission distance.

[0003] Based on the reality that the global atmospheric transmittance is not uniform, it is necessary to accurately evaluate the global atmospheric transmittance to enhance the reliability of the detection data. At present, the method of atmospheric transmittance inversion mainly relies on empirical models or single sounding data, and such inversion method is difficult to accurately invert the real-time global atmospheric transmittance. Although the reanalysis data can provide global coverage of temperature, pressure and other meteorological parameters, the vertical resolution of the reanalysis data is seriously mismatched with the spaceborne lidar observation, so the fusion of spaceborne lidar and reanalysis data is seriously restricted.

[0004] In view of the above prior art, it is urgent to develop a new transmittance inversion method that can cooperatively utilize the advantages of spaceborne lidar and reanalysis data to solve the core technical problems such as multi-source data fusion matching.

[0005] Therefore, the present application is proposed. SUMMARY

[0006] The spaceborne lidar and reanalysis data multi-wavelength atmospheric transmittance inversion method described in the application aims to solve the problems existing in the prior art and proposes an inversion method based on spaceborne lidar and reanalysis data cooperation, in order to solve the deficiency of real atmospheric scene in the spaceborne lidar simulation process, and to provide a multi-source data fusion solution.

[0007] To achieve the above design purpose, the spaceborne lidar and reanalysis data multi-wavelength atmospheric transmittance inversion method is based on spaceborne lidar data and reanalysis data, and combines wavelength conversion and data fusion to accurately calculate the atmospheric transmittance at different wavelengths; including the following implementation steps:

[0008] Step 1), selecting spaceborne lidar data;

[0009] Download aerosol profile data and vertical feature data, including but not limited to cloud-aerosol discrimination score, extinction QC flag, extinction uncertainty data, aerosol layer classification and aerosol extinction coefficient;

[0010] Step 2), screening the extinction coefficient;

[0011] Step 3), generate spatiotemporal matching dataset;

[0012] Obtain aerosol extinction coefficient and classification data from spaceborne lidar data, as well as temperature / pressure / ozone mixing ratio data from reanalysis data. Use a spatiotemporal sliding window matching algorithm to align the lidar data and reanalysis data at the desired spatiotemporal resolution to generate a coordinated dataset.

[0013] Step 4), establishing an aerosol extinction-transmittance conversion model;

[0014] Dynamically select wavelength conversion parameters based on aerosol classification data and establish an aerosol extinction-transmittance conversion model;

[0015] Step 5) Calculate the atmospheric molecules and ozone transmittance;

[0016] Combining reanalysis data, the atmospheric molecule and ozone transmittances are calculated using atmospheric molecule scattering theory and ozone absorption spectrum;

[0017] Step 6), layered transmittance fusion;

[0018] The atmospheric molecules, ozone, and aerosol transmittances are unified to the coordinate system of the same altitude using the pressure layer interpolation method;

[0019] Step 7) Calculate the transmittance of the entire atmosphere;

[0020] The atmospheric transmittance of each layer is obtained by multiplying the atmospheric molecules, ozone, and aerosol transmittance at the same height obtained in step 6); the atmospheric transmittance of each layer is multiplied to obtain the atmospheric transmittance of the entire layer.

[0021] Furthermore, in step 2), the spaceborne lidar data obtained in step 1) are preprocessed, including but not limited to screening effective extinction coefficient data based on the cloud-aerosol discrimination score, eliminating data according to the extinction QC flag, excluding data points with uncertainty exceeding the threshold, zeroing the extinction coefficient of areas marked as clean atmosphere, removing aerosol data interfered by high-altitude ice clouds, and limiting the reasonable value range of the extinction coefficient.

[0022] Furthermore, in step 3), the spaceborne lidar data is resampled according to the global grid resolution of the reanalysis data, and the data within the lidar transit track is retained to generate a spatial matching data set;

[0023] The reanalysis data is time-matched based on the aerosol extinction coefficient collected by the spaceborne lidar. If the reanalysis data used is daily average data, the time nearest neighbor interpolation method is used to match the precise collection time of the lidar. If the reanalysis data used is monthly average data, all the lidar transit data in the month is matched with the reanalysis data set of the month. The collaborative data set is obtained.

[0024] Further, the step 4) classifies the aerosol of the spaceborne lidar, the non-spherical aerosol uses T matrix to calculate the extinction coefficient, and the spherical aerosol uses Mie scattering to calculate the extinction coefficient, so as to obtain exponent;

[0025] According to The extinction coefficient of the aerosol at the target wavelength is calculated, wherein σ(λ2) is the extinction coefficient of the target wavelength after conversion, σ(λ1) is the extinction coefficient of the spaceborne lidar wavelength, is exponent;

[0026] According to The aerosol optical thickness is calculated and compared with the measured site of AERONET, and then The aerosol transmittance at the target wavelength is calculated; wherein Z C represents the ground height, Z sat represents the highest height of satellite detection.

[0027] Further, the step 5) only considers Rayleigh scattering when calculating the atmospheric molecular scattering, and the atmospheric molecular extinction coefficient is calculated according to the following formula:

[0028]

[0029] Wherein, N A = 6.02214 × 10 23 (1 / mol) is Avogadro's constant, R a = 8.314472 (J / K / mol) is the molar gas constant, P(Z) is the air pressure, and T(Z) is the temperature; Q s (λ) is the total Rayleigh scattering cross section caused by each atmospheric molecule in the "standard air", and the calculation formula is as follows:

[0030]

[0031] The transmittance of the atmospheric molecule is as follows:

[0032]

[0033] When calculating the ozone extinction coefficient, the atmospheric density needs to be calculated first:

[0034]

[0035] Where ρ(Z) is the atmospheric density, R is the gas constant (287.058 J / (kg·K)), P(Z) and T(Z) are the pressure and temperature, respectively. The ozone mass mixing ratio is converted to a column density per kilometer:

[0036]

[0037] Among them, r o (Z) is the ozone mixing ratio obtained by ERA5;

[0038] The ozone transmittance is calculated as follows:

[0039]

[0040] Among them, c o is the ozone absorption coefficient.

[0041] Furthermore, in step 6), the extinction coefficients of atmospheric molecules and ozone calculated in step 5) are highly matched using a layered adaptive interpolation method based on the vertical resolution of the spaceborne lidar data, and then the transmittance of atmospheric molecules and ozone is calculated;

[0042] A nonlinear interpolation algorithm is used for high-precision matching in the aerosol dense layer, which is determined by an extinction coefficient threshold of >0.01km. -1 Determination: Use linear interpolation algorithm for efficient matching in the clean atmosphere, the clean atmosphere passes the extinction coefficient threshold ≤ 0.001km -1 Determination: The extinction coefficient in the transition layer is 0.001-0.01km -1 A weighted combination of the nonlinear interpolation algorithm and the linear interpolation algorithm is used for matching. A weighted combination of the nonlinear interpolation algorithm and the linear interpolation algorithm is used for matching.

[0043] Furthermore, in step 7), the transmittance of each component is fused through an algorithm to obtain the transmittance of the entire atmosphere;

[0044] The fusion expression of the transmittance T at each vertical height is: T = T α ·T m ·T o ;

[0045] The transmittance of the entire atmosphere T A , the calculation formula is: T A =T1·T2·····T n ; Where n is the total number of layers in the vertical direction.

[0046] In summary, the advantages and beneficial effects of the multi-wavelength atmospheric transmittance inversion method based on spaceborne lidar and reanalysis data proposed in this application are as follows:

[0047] By synergistically utilizing the high-resolution observation data of spaceborne lidar and the global coverage advantage of reanalysis data, this application proposes a multi-source data spatiotemporal matching algorithm and a layered adaptive interpolation method, thereby achieving high-precision inversion of aerosols, atmospheric molecules and ozone transmittance.

[0048] This application effectively addresses the vertical resolution mismatch and provides more reliable data support for aerosol-cloud interaction analysis, through dynamic wavelength parameter selection and layered fusion technology. Furthermore, a method for accurately calculating atmospheric transmittance at different wavelengths is proposed, providing technical support for spaceborne lidar simulation and data quality. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] Figure 1 This is a flow chart of the method for inverting multi-wavelength atmospheric transmittance from spaceborne lidar and reanalysis data described in this application;

[0050] Figure 2 is a schematic diagram of the extinction coefficient at each wavelength after wavelength conversion of the 532 nm aerosol profile using the present application;

[0051] Figure 3 This is a comparison chart of the aerosol optical depth calculated using this application and the aerosol optical depth measured at the site;

[0052] Figure 4-1 and Figure 4-2 They are schematic diagrams of the extinction coefficient profiles of atmospheric molecules and ozone after matching; DETAILED DESCRIPTION

[0053] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts shall fall within the scope of protection of the present invention.

[0054] Example 1, as Figure 1 As shown in Figure 4, a multi-wavelength atmospheric transmittance inversion method based on spaceborne lidar and reanalysis data is based on spaceborne lidar data and reanalysis data, combining wavelength conversion and data fusion to accurately calculate the atmospheric transmittance at different wavelengths.

[0055] The atmospheric transmittance refers to the effective proportion of light rays transmitted in the atmosphere from a source point to a target point. When light passes through the earth's atmosphere, it is affected by scattering, absorption of atmospheric molecules, aerosols, and ozone, etc., thereby causing signal attenuation. The atmospheric transmittance is of great significance in the fields of astronomical observation, remote sensing detection, laser communication, etc.

[0056] The aerosol refers to a gaseous dispersion system composed of solid or liquid particles suspended in a gaseous medium.

[0057] The extinction coefficient refers to the relative attenuation rate of electromagnetic wave radiation per unit distance in the atmosphere.

[0058] The aerosol optical thickness refers to the optical thickness caused by the absorption and scattering of solar radiation by aerosols in the atmosphere. The aerosol optical thickness is the integral of the aerosol extinction coefficient in the vertical direction, and is an important parameter for describing the degree of aerosol pollution, and is used to measure the degree of attenuation of solar radiation by aerosol particles in the vertical atmospheric column. The specific calculation formula can be as follows:

[0059]

[0060] In the formula, τ a (λ) represents the aerosol optical thickness at wavelength λ, α a (λ, z) represents the aerosol extinction coefficient at wavelength λ at height z.

[0061] The atmospheric pressure refers to the atmospheric pressure acting on a unit area, that is, in numerical terms, it is equal to the weight of the vertical air column extending to the upper boundary of the atmosphere on a unit area.

[0062] The temperature refers to the heat state of a substance, and is usually used to describe the heat level of an object.

[0063] The profile refers to the curve of a parameter of the atmosphere or the earth's surface varying with height or distance.

[0064] The ozone mass mixing ratio refers to the ozone mass per unit volume.

[0065] The spaceborne lidar and reanalysis data multi-wavelength atmospheric transmittance inversion method comprises the following implementation steps:

[0066] Step 1), selecting spaceborne lidar data;

[0067] Downloading aerosol profile data and vertical feature data, including but not limited to cloud-aerosol discrimination score, extinction QC flag, extinction uncertainty data, aerosol layer classification, and aerosol extinction coefficient;

[0068] Step 2), screening the extinction coefficient;

[0069] The spaceborne lidar data obtained in step 1) is preprocessed, including but not limited to screening effective extinction coefficient data based on cloud-aerosol discrimination score, selecting extinction coefficient data with cloud-aerosol discrimination score between-100 and-20;

[0070] Low-confidence data is removed according to the extinction QC flag, and extinction coefficient data with extinction QC flag of 0, 1, 16, and 18 is removed;

[0071] Data points with uncertainty exceeding a threshold are excluded, and extinction coefficient data with extinction uncertainty of 99.9 km -1 is removed;

[0072] The extinction coefficient of the region identified as clean atmosphere is set to 0.0 km -1 ;

[0073] Ice cloud interference aerosol data at high altitudes is removed, and independent aerosol near ice clouds above 4 km is removed;

[0074] The reasonable value range of the extinction coefficient is limited, and data with extinction coefficient between 0-1.25 km -1 is selected;

[0075] Step 3), generate spatiotemporal matching data set;

[0076] Obtain the aerosol extinction coefficient and classification data in the spaceborne lidar data, and the temperature / pressure / ozone mixing ratio data in the reanalysis data, and use the spatiotemporal sliding window matching algorithm to align the lidar data with the reanalysis data according to the spatiotemporal resolution to generate a collaborative data set;

[0077] Specifically, the spaceborne lidar data is resampled according to the global gridded resolution of the reanalysis data, and data within the range of the lidar transit track is retained to generate a spatially matched data set;

[0078] The reanalysis data is time-matched based on the aerosol extinction coefficient measurement time of the spaceborne lidar; if the reanalysis data used is daily average data, the time nearest neighbor interpolation method is used to match the accurate collection time of the lidar; if the reanalysis data used is monthly average data, all lidar transit data in the month is matched with the monthly reanalysis data set;

[0079] Get the collaborative data set;

[0080] Step 4), establish an aerosol extinction-transmittance conversion model;

[0081] Based on the aerosol classification data, the wavelength conversion parameters are dynamically selected to establish an aerosol extinction-transmittance conversion model;

[0082] The aerosol extinction-transmittance conversion model is used to classify aerosols for spaceborne lidar. The extinction coefficient of non-spherical aerosols is calculated using the T matrix, while the extinction coefficient of spherical aerosols is calculated using Mie scattering to obtain index;

[0083] according to Calculate the aerosol extinction coefficient at the target wavelength. The calculation results are as follows: Figure 2 As shown;

[0084] Among them, σ(λ2) is the extinction coefficient of the target wavelength after conversion, σ(λ1) is the extinction coefficient of the spaceborne lidar wavelength, for index;

[0085] according to The aerosol optical depth was calculated and compared with the AERONET measurement station. The results are as follows Figure 3 As shown;

[0086] Based on Calculate the aerosol transmittance at the target wavelength;

[0087] Among them, Z C Represents the ground height, Z sat Represents the highest altitude detected by the satellite;

[0088] It should be noted that Mie scattering is a theory describing the interaction between light waves and particles, first proposed by German physicist Gustav Mie in 1908. Mie scattering primarily applies to situations where the particle size is comparable to or larger than the wavelength of the incident light. It describes the scattering of light when it interacts with spherical particles (such as water droplets and aerosols). In Mie scattering, the intensity and distribution of the scattered light depend not only on the wavelength of the incident light and the size of the particle, but also on its shape and refractive index.

[0089] T-matrix scattering is a mathematical method that describes the interaction between scatterers of arbitrary shapes and electromagnetic waves. The T-matrix is ​​an effective representation of scattering, which contains information about the geometric shape, refractive index and incident wave of the scatterer. Compared with traditional scattering theory, T-matrix scattering can be used to analyze particles of complex shapes and provide information about the intensity and phase of scattered light. The calculation of the T-matrix is ​​based on the wave vector of the incident light and the characteristics of the scatterer, and can describe the distribution characteristics of the scattered light in different directions;

[0090] Step 5) Calculate the atmospheric molecules and ozone transmittance;

[0091] Combining reanalysis data, atmospheric molecule and ozone transmittances are calculated using atmospheric molecule scattering theory and ozone absorption spectrum;

[0092] When calculating atmospheric molecular scattering, only Rayleigh scattering is considered and the atmospheric molecular extinction coefficient is calculated according to the following formula:

[0093]

[0094] Among them, N A =6.02214×10 23 (1 / mol) is Avogadro's constant, R a =8.314472 (J / K / mol) is the molar gas constant, P(Z) is the gas pressure, T(Z) is the temperature; Q s (λ) is the total Rayleigh scattering cross section caused by each atmospheric molecule in "standard air", which is calculated as follows:

[0095]

[0096] The transmittance of atmospheric molecules is shown below:

[0097]

[0098] When calculating the ozone extinction coefficient, you need to first calculate the atmospheric density:

[0099]

[0100] Where ρ(Z) is the atmospheric density, R is the gas constant (287.058 J / (kg·K)), P(Z) and T(Z) are the pressure and temperature, respectively. The ozone mass mixing ratio is converted to a column density per kilometer:

[0101]

[0102] Among them, r o (Z) is the ozone mixing ratio obtained by ERA5;

[0103] The ozone transmittance is calculated as follows:

[0104]

[0105] Among them, c o is the ozone absorption coefficient;

[0106] Step 6), layered transmittance fusion;

[0107] The atmospheric molecules, ozone, and aerosol transmittances are unified to the coordinate system of the same altitude using the pressure layer interpolation method;

[0108] Specifically, according to the atmospheric molecules and ozone extinction coefficients calculated in step 5), a hierarchical adaptive interpolation method is used to perform height matching based on the vertical resolution of the spaceborne lidar data, and then the atmospheric molecules and ozone transmittance are calculated;

[0109] A nonlinear interpolation algorithm is used for high-precision matching in the aerosol dense layer, which is determined by an extinction coefficient threshold of >0.01km. -1 Determination: Use linear interpolation algorithm for efficient matching in the clean atmosphere, the clean atmosphere passes the extinction coefficient threshold ≤ 0.001km -1 Determination: The extinction coefficient in the transition layer is 0.001-0.01km -1 The weighted combination of the nonlinear interpolation algorithm and the linear interpolation algorithm is used for matching; the extinction coefficient profile after interpolation is shown in FIG4 ;

[0110] Step 7) Calculate the transmittance of the entire atmosphere;

[0111] The atmospheric transmittance of each layer is obtained by multiplying the atmospheric molecules, ozone, and aerosol transmittances at the same height obtained in step 6) and multiplying the atmospheric transmittance of each layer to obtain the atmospheric transmittance of the entire layer;

[0112] Specifically, the transmittance of each component is fused through the algorithm to obtain the transmittance of the entire atmosphere;

[0113] The fusion expression of the transmittance T at each vertical height is: T = T α ·T m ·T o ;

[0114] The transmittance of the entire atmosphere T A , the calculation formula is: T A =T1·T2·····T n ; Where n is the total number of layers in the vertical direction.

[0115] It should be understood that those skilled in the art can make improvements or changes based on the above description, and all such improvements and changes should fall within the scope of protection of the appended claims of the present invention.

Claims

1. A method for inverting multi-wavelength atmospheric transmittance from spaceborne lidar and reanalysis data, characterized by: Based on spaceborne lidar data and reanalysis data, wavelength conversion and data fusion are combined to accurately calculate the atmospheric transmittance at different wavelengths; The implementation steps include the following: Step 1) Select spaceborne lidar data; Download aerosol profile data and vertical characteristic data, including but not limited to cloud-aerosol discrimination fraction, extinction QC flags, extinction uncertainty data, aerosol layer classification, and aerosol extinction coefficient; Step 2), screening the extinction coefficient; Step 3), generate spatiotemporal matching dataset; Obtain aerosol extinction coefficient and classification data from spaceborne lidar data, as well as temperature / pressure / ozone mixing ratio data from reanalysis data. Use a spatiotemporal sliding window matching algorithm to align the lidar data and reanalysis data at the desired spatiotemporal resolution to generate a coordinated dataset. Step 4), establishing an aerosol extinction-transmittance conversion model; Dynamically select wavelength conversion parameters based on aerosol classification data and establish an aerosol extinction-transmittance conversion model; Step 5) Calculate the atmospheric molecules and ozone transmittance; Combining reanalysis data, the atmospheric molecule and ozone transmittances are calculated using atmospheric molecule scattering theory and ozone absorption spectrum; Step 6), layered transmittance fusion; The atmospheric molecules, ozone, and aerosol transmittances are unified to the coordinate system of the same altitude using the pressure layer interpolation method; Step 7) Calculate the transmittance of the entire atmosphere; The atmospheric molecules, ozone, and aerosol transmittances at the same height obtained in step 6) are multiplied to obtain the atmospheric transmittance of each layer; the atmospheric transmittance of each layer is multiplied to obtain the atmospheric transmittance of the entire layer.

2. The method for inverting multi-wavelength atmospheric transmittance using spaceborne lidar and reanalysis data according to claim 1, characterized in that: The step 2) is to preprocess the spaceborne lidar data obtained in step 1), including but not limited to screening effective extinction coefficient data based on the cloud-aerosol discrimination score, eliminating low-confidence data according to quality control marks, excluding data points with uncertainty exceeding the threshold, zeroing the extinction coefficient of areas marked as clean atmosphere, removing aerosol data interfered by high-altitude ice clouds, and limiting the reasonable value range of the extinction coefficient.

3. The method for inverting multi-wavelength atmospheric transmittance using spaceborne lidar and reanalysis data according to claim 1, characterized in that: In step 3), the spaceborne lidar data is resampled according to the global grid resolution of the reanalysis data, and the data within the lidar transit track is retained to generate a spatial matching data set; The reanalysis data are time-matched based on the acquisition time of the aerosol extinction coefficient measured by the spaceborne lidar. If the reanalysis data used are daily average data, the temporal nearest neighbor interpolation method is used to match the lidar's precise acquisition time. If the reanalysis data used are monthly average data, all lidar transit data in that month are matched with the reanalysis dataset of that month to obtain a collaborative dataset.

4. The method for inverting multi-wavelength atmospheric transmittance using spaceborne lidar and reanalysis data according to claim 1, characterized in that: In step 4), the aerosols of the space-borne laser radar are classified. The extinction coefficient of non-spherical aerosols is calculated using the T matrix, and the extinction coefficient of spherical aerosols is calculated using Mie scattering to obtain index; according to Calculate the aerosol extinction coefficient at the target wavelength, where σ(λ2) is the extinction coefficient of the converted target wavelength, and σ(λ1) is the extinction coefficient of the spaceborne lidar wavelength. for index; according to Calculate the aerosol optical depth and compare it with the AERONET measurement site. Calculate the aerosol transmittance at the target wavelength; where Z C Represents the ground height, Z sat Represents the highest altitude of satellite detection.

5. The method for inverting multi-wavelength atmospheric transmittance using spaceborne lidar and reanalysis data according to claim 1, characterized in that: In step 5), only Rayleigh scattering is considered when calculating the atmospheric molecular scattering, and the atmospheric molecular extinction coefficient is calculated according to the following formula: Among them, N A =6.02214×10 23 (1 / mol) is Avogadro's constant, R a =8.314472 (J / K / mol) is the molar gas constant, P(Z) is the gas pressure, T(Z) is the temperature; Q s (λ) is the total Rayleigh scattering cross section caused by each atmospheric molecule in "standard air", which is calculated as follows: The transmittance of atmospheric molecules is shown below: When calculating the ozone extinction coefficient, you need to first calculate the atmospheric density: Where ρ(Z) is the atmospheric density, R is the gas constant (287.058 J / (kg·K)), P(Z) and T(Z) are the pressure and temperature, respectively. The ozone mass mixing ratio is converted to a column density per kilometer: Among them, r o (Z) is the ozone mixing ratio obtained by ERA5; The ozone transmittance is calculated as follows: Among them, c o is the ozone absorption coefficient.

6. The method for inverting multi-wavelength atmospheric transmittance using spaceborne lidar and reanalysis data according to claim 1, characterized in that: In step 6), the atmospheric molecules and ozone extinction coefficients calculated in step 5) are highly matched using a layered adaptive interpolation method based on the vertical resolution of the spaceborne lidar data, and then the atmospheric molecules and ozone transmittance are calculated; A nonlinear interpolation algorithm is used for high-precision matching in the aerosol dense layer, which is determined by an extinction coefficient threshold of >0.01km. -1 Determination: Use linear interpolation algorithm for efficient matching in the clean atmosphere, the clean atmosphere passes the extinction coefficient threshold ≤ 0.001km -1 Determination: In the transition layer, the extinction coefficient is between 0.001 and 0.01 km. -1 A weighted combination of the nonlinear interpolation algorithm and the linear interpolation algorithm is used for matching.

7. The method for inverting multi-wavelength atmospheric transmittance using spaceborne lidar and reanalysis data according to claim 1, characterized in that: In step 7), the transmittance of each component is integrated through an algorithm to obtain the transmittance of the entire atmosphere; The fusion expression of the transmittance T at each vertical height is: T = T α ·T m ·T o ; The transmittance of the entire atmosphere T A , the calculation formula is: T A =T1·T2·····T n ; Where n is the total number of layers in the vertical direction.

Images