Aerosol layer height and optical thickness inversion method, device, equipment and medium

Abstract

The invention discloses an aerosol layer height and optical thickness inversion method and device, equipment and a medium. The method comprises the following steps: acquiring an atmosphere profile type, aerosol microphysics and single scattering characteristics, observation geometric information, observation longitude and latitude information, a polarization reflectivity measurement value of a first observation wave band combination and a reflectivity measurement value of a second observation wave band combination of an observation area; an aerosol inversion lookup table is established according to the atmospheric profile type, aerosol micro-physics and single scattering characteristics, observation geometric information, observation latitude and longitude information and the polarization reflectivity measurement value of the first observation wave band combination, and an aerosol inversion lookup table is established according to the reflectivity measurement value of the second observation wave band combination and the aerosol inversion lookup table. Therefore, the inversion aerosol layer height and the inversion aerosol optical thickness are obtained. According to the method, the defects of narrow product breadth and long revisit period of the aerosol layer height faced by active remote sensing are greatly overcome by utilizing passive remote sensing data, and the measurement precision of the aerosol layer height is improved.

Description

Technical Field

[0001] The invention belongs to the technical field of remote sensing image processing, and in particular relates to an aerosol layer height and optical thickness inversion method, device, equipment and medium. Background Art

[0002] In satellite earth imaging remote sensing in the visible to shortwave infrared bands, atmospheric scattering and multiple coupled scattering between the atmosphere and the surface blur the remote sensing images, reducing their clarity and causing the apparent reflectivity of the detected top of the atmosphere to deviate from the true surface reflectivity. Therefore, after collecting the original images of the surface using a satellite remote sensing platform, it is usually necessary to perform atmospheric correction on the original images to obtain more accurate image content and obtain the layer height distribution of aerosols, that is, the vertical distribution, which is of great significance for atmospheric correction research.

[0003] In order to measure the layer height distribution of aerosols, a lot of research has been done on aerosol layer height measurement since the Cloud-Aerosol Polarimetric Lidar (CALIOP) on the CALIPSO platform began to provide data in 2006. Active remote sensing technology using high spectral resolution Lidar (LIDAR) measurements can accurately measure the vertical distribution of aerosols, but there are limitations in width and revisit period. So far, passive satellites still cannot provide inverted aerosol vertical distribution data with the same accuracy as LIDAR, but in terms of spatial range and revisit period, their data are an indispensable supplement to LIDAR aerosol layer height inversion.

[0004] The O2A band (0.764 μm) and O2B band (0.688 μm) measurements provided by passive satellites have the advantage of being unaffected by thermal radiation, and have been shown to be of great value in retrieving aerosol layer height when the aerosol optical thickness is large. Algorithms that use oxygen absorption spectra to invert the vertical distribution of aerosols can be divided into two categories. The first category directly applies the principle that the apparent reflectivity of the oxygen absorption band increases with the rise of the aerosol scattering layer. This type of algorithm derives the aerosol height based on the reflectivity ratio of the inner and outer channels of the absorption band. However, this method can only produce a single aerosol height information. The second type of method, satellites that act on spectral fitting techniques, usually have coarse spatial resolution (e.g., GOME, GOME-2, and SCIAMCHY) or limited spatial coverage (e.g., GOSAT, OCO-2, and TanSat). In contrast, the reflectivity ratio method provides less detailed aerosol profiles than the spectral fitting algorithm, but the corresponding instruments usually have better spatial resolution and wider spatial coverage (e.g., POLDER, MERIS, and EPIC). Accurate height retrieval relies on appropriate assumptions about aerosol optical properties and surface reflectivity. Currently, there is an urgent need for passive satellite measurement techniques in the context of atmospheric correction to accurately retrieve the vertical distribution of the aerosol layer. Summary of the invention

[0005] The purpose of the present invention is to provide an aerosol layer height and optical thickness inversion method, device, equipment and medium, which realizes fast and accurate inversion of aerosol layer height and aerosol optical thickness based on passive multi-spectral satellite data by using optimal estimation technology and radiation transmission model, and solves the technical problem that the prior art is highly dependent on the prior assumption of surface reflectivity and it is difficult to obtain accurate aerosol layer height and aerosol optical thickness results under the goal of atmospheric correction.

[0006] In order to solve the above technical problems, the present invention is achieved through the following technical solutions:

[0007] The present invention provides an aerosol layer height and optical thickness inversion method, which comprises:

[0008] Obtain the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, polarization reflectivity measurement values ​​of the first observation band combination, and reflectivity measurement values ​​of the second observation band combination in the observation area;

[0009] According to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, and the polarization reflectivity measurement value of the first observation band combination, a priori aerosol layer height, a priori aerosol optical thickness, and a priori surface reflectivity are obtained;

[0010] According to the a priori aerosol layer height, the a priori aerosol optical thickness and the a priori surface reflectivity, an aerosol inversion lookup table is obtained;

[0011] The inverted aerosol layer height and the inverted aerosol optical thickness are obtained according to the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table.

[0012] In one embodiment of the present invention, the first observation band combination includes observation bands that are greatly affected by aerosols and atmospheric molecules and observation bands that are less affected by aerosols and atmospheric molecules, and the second observation bands include 688nm and 680nm, 764nm and 780nm.

[0013] In one embodiment of the present invention, the step of obtaining the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, polarization reflectivity measurement values ​​of the first observation band combination, and reflectivity measurement values ​​of the second observation band combination of the observation area includes:

[0014] Determine the atmospheric profile type and aerosol microphysics and single scattering characteristics based on the observation area and observation time of the passive multispectral satellite remote sensing platform;

[0015] Acquiring observation geometry information and observation latitude and longitude information collected by a positioning and attitude determination system carried by the passive multispectral satellite remote sensing platform;

[0016] The polarization reflectivity measurement value of the first observation band combination and the reflectivity measurement value of the second observation band combination collected by the atmospheric correction instrument carried by the passive multispectral satellite remote sensing platform are obtained.

[0017] In one embodiment of the present invention, the method of obtaining a priori aerosol layer height, a priori aerosol optical thickness and a priori surface reflectivity based on the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information and the polarization reflectivity measurement value of the first observation band combination comprises:

[0018] According to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information and the polarization reflectivity measurement value of the first observation band combination, the GRASP program is used to obtain the prior aerosol layer height, the prior aerosol optical thickness and the prior surface reflectivity.

[0019] In one embodiment of the present invention, the step of obtaining an aerosol inversion lookup table according to the prior aerosol layer height, the prior aerosol optical thickness and the prior surface reflectivity includes:

[0020] According to the a priori aerosol layer height, the a priori aerosol optical thickness and the a priori surface reflectivity, using an atmospheric radiation transfer model, obtaining theoretical values ​​of apparent reflectivity of the observation area corresponding to all observation bands in the first observation band combination;

[0021] An aerosol inversion lookup table is established according to the prior aerosol layer height, the prior aerosol optical thickness, the prior surface reflectivity and the theoretical value of the apparent reflectivity.

[0022] In one embodiment of the present invention, the aerosol inversion lookup table is established according to the prior aerosol layer height, the prior aerosol optical thickness, the prior surface reflectivity and the theoretical value of the apparent reflectivity, including:

[0023] Construct an M×N two-dimensional table, in which the horizontal rows represent the aerosol layer height and the vertical rows represent the aerosol optical thickness. Fill the corresponding theoretical values ​​of the apparent reflectivity into the two-dimensional table to obtain the aerosol inversion lookup table, in which each grid point represents the mapping relationship between the aerosol layer height, the aerosol optical thickness and the apparent reflectivity.

[0024] In one embodiment of the present invention, the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table are used to obtain the inverted aerosol layer height and the inverted aerosol optical thickness, including:

[0025] Obtaining the apparent reflectivity of the second observation band combination according to the reflectivity measurement value of the second observation band combination;

[0026] The apparent reflectivity of the second observation band combination is compared and verified with the aerosol inversion lookup table, and the inverted aerosol layer height and the inverted aerosol optical thickness are obtained by interpolating and fitting the aerosol inversion lookup table.

[0027] Based on the same inventive concept, another embodiment of the present invention further provides an aerosol layer height and optical thickness inversion device, the device comprising:

[0028] A data acquisition module is used to obtain the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, polarization reflectivity measurement values ​​of the first observation band combination, and reflectivity measurement values ​​of the second observation band combination in the observation area;

[0029] A data processing module, for obtaining a priori aerosol layer height, a priori aerosol optical thickness and a priori surface reflectivity according to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information and the polarization reflectivity measurement value of the first observation band combination, and obtaining an aerosol inversion lookup table according to the priori aerosol layer height, the priori aerosol optical thickness and the prior surface reflectivity;

[0030] The aerosol inversion module is used to obtain the inverted aerosol layer height and the inverted aerosol optical thickness according to the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table.

[0031] Based on the same inventive concept, another embodiment of the present invention further provides an electronic device, the electronic device comprising:

[0032] one or more processors;

[0033] A storage device is used to store one or more programs. When the one or more programs are executed by the one or more processors, the electronic device implements the aerosol layer height and optical thickness inversion method described in any of the above embodiments.

[0034] Based on the same inventive concept, another embodiment of the present invention further provides a computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor of a computer, the computer executes the aerosol layer height and optical thickness inversion method described in any of the above embodiments.

[0035] As described above, the present invention provides an aerosol layer height and optical thickness inversion method, which has the following beneficial effects: by obtaining the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information, polarization reflectivity measurement values ​​of the first observation band combination and reflectivity measurement values ​​of the second observation band combination in the observation area, according to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information and polarization reflectivity measurement values ​​of the first observation band combination, a priori aerosol layer height, a priori aerosol optical thickness and a priori surface reflectivity are obtained; according to the prior aerosol layer height, the prior aerosol optical thickness and the prior surface reflectivity, an aerosol inversion lookup table is obtained; according to the reflectivity measurement values ​​of the second observation band combination and the aerosol inversion lookup table, the inverted aerosol layer height and the inverted aerosol optical thickness are obtained. The aerosol layer height and optical thickness inversion method can obtain aerosol layer height data with a larger width and a shorter revisit period than LIDAR through passive multispectral satellites, which is a major breakthrough in the field of satellite inversion of atmospheric components. The multispectral data of the observation area collected by passive multispectral satellites and combined with the forward radiation transfer model can be inverted to obtain the aerosol layer height parameters above the observation area. By using passive remote sensing data, the defects of narrow width and long revisit period of aerosol layer height products faced by active remote sensing (CALIOP, LIDAR) are greatly compensated. Of course, any product implementing the present invention does not necessarily need to achieve all the advantages described above at the same time. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for describing the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.

[0037] Figure 1 A schematic flow chart of an aerosol layer height and optical thickness inversion method provided for an exemplary embodiment of the present application.

[0038] Figure 2 A system architecture block diagram of a forward radiation transmission model provided for an exemplary embodiment of the present application.

[0039] Figure 3 A schematic diagram of the aerosol layer height and optical thickness inversion principle provided for an exemplary embodiment of the present application.

[0040] Figure 4 A schematic structural diagram of an aerosol layer height and optical thickness inversion device provided for another exemplary embodiment of the present application.

[0041] Figure 5 A schematic structural diagram of an electronic device provided for another exemplary embodiment of the present application. DETAILED DESCRIPTION

[0042] The following describes the embodiments of the present invention by specific examples, and those skilled in the art can easily understand other advantages and effects of the present invention from the contents disclosed in this specification. The present invention can also be implemented or applied through other different specific embodiments, and the details in this specification can also be modified or changed in various ways based on different viewpoints and applications without departing from the spirit of the present invention. It should be noted that the following embodiments and features in the embodiments can be combined with each other without conflict.

[0043] It should be noted that the illustrations provided in the following embodiments are only schematic illustrations of the basic concept of the present invention, and thus the drawings only show components related to the present invention rather than being drawn according to the number, shape and size of components in actual implementation. In actual implementation, the type, quantity and proportion of each component may be changed arbitrarily, and the component layout may also be more complicated.

[0044] In the following description, numerous details are discussed to provide a more thorough explanation of the embodiments of the present invention. However, it is obvious to those skilled in the art that the embodiments of the present invention can be implemented without these specific details. In other embodiments, well-known structures and devices are shown in the form of block diagrams rather than in detail to avoid making the embodiments of the present invention difficult to understand.

[0045] In order to solve the technical problem that the existing technology is highly dependent on the prior assumption of surface reflectivity and it is difficult to obtain accurate aerosol layer height and aerosol optical thickness results under the goal of atmospheric correction, the present invention provides an aerosol layer height and optical thickness inversion method, which uses the optimal estimation technology and radiation transfer model based on passive multi-spectral satellite data to achieve the inversion of aerosol layer height and aerosol optical thickness. Please refer to Figure 1 As shown, the aerosol layer height and optical thickness inversion method includes the following steps:

[0046] S100: Obtain the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, polarization reflectivity measurement value of the first observation band combination, and reflectivity measurement value of the second observation band combination in the observation area;

[0047] S200: obtaining a priori aerosol layer height, a priori aerosol optical thickness and a priori surface reflectivity according to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information and the polarization reflectivity measurement value of the first observation band combination;

[0048] S300: Obtaining an aerosol inversion lookup table according to the prior aerosol layer height, the prior aerosol optical thickness and the prior surface reflectivity;

[0049] S400: Obtaining an inverted aerosol layer height and an inverted aerosol optical thickness according to the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table.

[0050] The steps of the above-mentioned aerosol layer height and optical thickness inversion method are discussed in detail below.

[0051] It should be noted that, in this embodiment, the aerosol layer height and optical thickness inversion method is applied to a passive multi-spectral satellite remote sensing platform. The passive multi-spectral satellite remote sensing platform can be equipped with an atmospheric correction instrument, a positioning and attitude system, a controller and a data processor simultaneously. Among them, the working bands of the atmospheric correction instrument include: 443nm, 490nm, 565nm, 670nm, 870nm, 910nm, 680nm, 688nm, 764nm, 730nm, 1380nm, 1500nm, 1560nm, 1610nm, and the atmospheric correction instrument adopts a three-channel solution, that is, three different polarization analysis directions are set in each spectral band, and the polarization analysis directions are usually set to 0°, 60° and 120° to cover as much polarization information as possible, which can obtain a wide range of atmospheric multi-spectral polarization data. The polarization observation band can include 680nm, 688nm, 764nm, and 780nm, and the remaining bands use a non-polarization observation mode to observe reflectivity. The positioning and orientation system (POS) can obtain the geographical information on the flight path of the aerial remote sensing platform, flight altitude and other flight trajectory data, and obtain the attitude data of the atmospheric correction instrument. The controller can control the atmospheric correction instrument, the positioning and orientation system, etc. The data processor can be used to execute the aerosol layer height and optical thickness inversion method.

[0052] First, step S100 is executed to obtain the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, polarization reflectivity measurement values ​​of the first observation band combination, and reflectivity measurement values ​​of the second observation band combination in the observation area.

[0053] In an exemplary embodiment of the present application, in step S100, obtaining the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, polarization reflectivity measurement values ​​of the first observation band combination, and reflectivity measurement values ​​of the second observation band combination of the observation area further includes the following steps:

[0054] S110: Determine the atmospheric profile type and aerosol microphysics and single scattering characteristics based on the observation area and observation time of the passive multispectral satellite remote sensing platform;

[0055] S120: Acquire observation geometry information and observation longitude and latitude information collected by a positioning and attitude determination system carried by the passive multispectral satellite remote sensing platform;

[0056] S130: Obtaining polarization reflectivity measurement values ​​of the first observation band combination and reflectivity measurement values ​​of the second observation band combination collected by the atmospheric correction instrument carried by the passive multispectral satellite remote sensing platform.

[0057] Specifically, the positioning and attitude determination system carried by the passive multi-spectral satellite remote sensing platform collects observation geometry information and observation longitude and latitude, wherein the observation geometry information includes the relative position and angle relationship between the satellite, the observation area and the sun. In this embodiment, the observation geometry information includes satellite attitude, observation angle and altitude information, wherein the satellite attitude includes the roll angle, pitch angle and heading angle of the satellite, wherein the observation angle includes the incident angle, observation angle and azimuth angle, wherein the altitude information includes the height of the satellite relative to the earth's surface, and wherein the observation longitude and latitude include the geographical location of the observation area. The multi-spectral high-resolution camera and the atmospheric correction instrument carried by the passive multi-spectral satellite remote sensing platform can synchronously collect surface image data and polarization reflectance data of the first observation band combination and reflectance data of the second observation band combination under the control of the controller. After the data processor obtains the raw data collected by the two, it can perform pre-processing such as radiation calibration, polarization calibration and geometric calibration on the data to obtain the surface image that can be used for subsequent processing, the polarization reflectance measurement value corresponding to the first observation band combination and the reflectance measurement value corresponding to the second observation band combination.

[0058] In an exemplary embodiment of the present application, the first observation band combination includes observation bands that are more affected by aerosols and atmospheric molecules and observation bands that are less affected by aerosols and atmospheric molecules, and the second observation band includes 688nm and 680nm, 764nm and 780nm.

[0059] Specifically, the first observation band combination includes: 443nm, 490nm, 565nm, 670nm, 870nm, 910nm, 680nm, 688nm, 764nm, 730nm, 1380nm, 1500nm, 1560nm, 1610nm. It should be noted that since the polarization characteristics of the surface are not sensitive to the band, and generally the longer the wavelength, the less it is affected by atmospheric molecules and aerosol scattering, in this embodiment, the first observation band combination includes observation bands that are greatly affected by aerosols and atmospheric molecules and observation bands that are less affected by aerosols and atmospheric molecules. It is worth noting that in this embodiment, the polarization reflectivity measurement value of the polarization observation band of the atmospheric correction instrument can be used as the input of the GRASP program, and the prior aerosol layer height, prior aerosol optical thickness and prior surface reflectivity are calculated by the GRASP program during the subsequent aerosol inversion.

[0060] The second observation band combination includes 688nm and 680nm, 764nm and 780nm, wherein the 688nm band is an oxygen B absorption band, the 680nm band is an oxygen B non-absorption band, the 764nm band is an oxygen A absorption band, and the 780nm band is an oxygen A non-absorption band. Since the surface characteristics and aerosol characteristics vary little in the four bands in the second observation band combination, the influence of the surface characteristics and aerosol characteristics on the reflectivity can be eliminated by the ratio of the reflectivity measurement values ​​of the absorption band and the non-absorption band.

[0061] Next, step S200 is executed to obtain a priori aerosol layer height, a priori aerosol optical thickness and a priori surface reflectivity according to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information and the polarization reflectivity measurement value of the first observation band combination.

[0062] In an exemplary embodiment of the present application, step S200 includes: using the GRASP program to obtain the prior aerosol layer height, the prior aerosol optical thickness and the prior surface reflectivity according to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information and the polarization reflectivity measurement value of the first observation band combination.

[0063] Specifically, the atmospheric correction instrument obtains corresponding polarization reflectivity measurement values ​​in the polarization observation mode for each observation band of the first observation band combination, and determines the initial value of aerosol layer height, initial value of aerosol optical thickness and initial value of surface reflectivity in the observation area through the GRASP program according to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information and the polarization reflectivity measurement values ​​of the first observation band combination, and takes the initial value of aerosol layer height corresponding to each band in the first observation band combination as the prior aerosol layer height, the initial value of aerosol optical thickness corresponding to each band in the first observation band combination as the prior aerosol optical thickness, and the initial value of surface reflectivity as the prior surface reflectivity of the observation area.

[0064] Next, step S300 is performed, that is, an aerosol inversion lookup table is obtained according to the prior aerosol layer height, the prior aerosol optical thickness and the prior surface reflectivity.

[0065] In an exemplary embodiment of the present application, in step S300, according to the a priori aerosol layer height, the a priori aerosol optical thickness and the a priori surface reflectivity, an aerosol inversion lookup table is obtained, further comprising the following steps:

[0066] S310: according to the a priori aerosol layer height, the a priori aerosol optical thickness and the a priori surface reflectivity, using an atmospheric radiation transmission model, obtaining theoretical values ​​of apparent reflectivity of the observation area corresponding to all observation bands in the first observation band combination;

[0067] S320: Establishing an aerosol inversion lookup table according to the prior aerosol layer height, the prior aerosol optical thickness, the prior surface reflectivity and the theoretical value of the apparent reflectivity.

[0068] In an exemplary embodiment of the application, step S320 includes: constructing an M×N two-dimensional table, in which the horizontal rows represent the aerosol layer height and the vertical rows represent the aerosol optical thickness, and filling the corresponding theoretical values ​​of the apparent reflectivity into the two-dimensional table to obtain the aerosol inversion lookup table, wherein each grid point represents the mapping relationship between the aerosol layer height, the aerosol optical thickness and the apparent reflectivity.

[0069] Specifically, in this embodiment, please refer to Figure 2As shown, the atmospheric radiation transfer model can be a unified linearized vector radiation transfer model (UNL-VRTM). By inputting multiple pairs of a priori aerosol layer height, a priori aerosol optical thickness, and observation geometry information and observation latitude and longitude information of a passive multi-spectral satellite remote sensing platform into the model, multiple theoretical values ​​of apparent reflectance corresponding to oxygen absorption bands and non-oxygen absorption bands can be calculated. The expression of the UNL-VRTM radiation transfer model is:

[0070]

[0071] where τ is the aerosol extinction coefficient measured at the top of the atmosphere, Δφ=φ-φ0 is the relative azimuth, μ and μ0 are the cosines of φ and φ0, respectively, ω is the single scattering albedo, and Z is the scattering phase function matrix.

[0072] According to the parameter range and step size of the variable parameters, different values ​​of different variable parameters are input into the UNL-VRTM radiation transfer model to obtain the corresponding theoretical values ​​of apparent reflectivity, and a lookup table of aerosol layer height and aerosol optical thickness is constructed based on this. That is, an M×N two-dimensional table is constructed, in which the horizontal rows represent the aerosol layer height and the vertical rows represent the aerosol optical thickness. The corresponding theoretical values ​​of apparent reflectivity are filled into the corresponding grids to obtain an aerosol inversion lookup table, in which each grid point represents the mapping relationship between aerosol layer height, aerosol optical thickness and apparent reflectivity.

[0073] It should be noted that in order to facilitate the inversion of aerosol layer height and aerosol optical thickness under different atmospheric profiles and observation geometry information, multiple different atmospheric profiles and a variety of different observation geometry information can be set in advance. According to the combination of different aerosol layer heights and aerosol optical thicknesses, observation altitudes and observation geometry information, the atmospheric radiation transfer model is used to obtain the theoretical values ​​of the apparent reflectivity corresponding to the oxygen absorption band and the non-oxygen absorption band for each combination, and an aerosol inversion lookup table is constructed.

[0074] Finally, step 400 is performed, that is, the inverted aerosol layer height and the inverted aerosol optical thickness are obtained according to the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table.

[0075] In an exemplary embodiment of the present application, in step S400, according to the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table, the inverted aerosol layer height and the inverted aerosol optical thickness are obtained, further comprising the following steps:

[0076] S410: Obtaining the apparent reflectivity of the second observation band combination according to the reflectivity measurement value of the second observation band combination;

[0077] S420: Compare and verify the apparent reflectivity of the second observation band combination with the aerosol inversion lookup table, and obtain the inverted aerosol layer height and inverted aerosol optical thickness by interpolating and fitting the aerosol inversion lookup table.

[0078] Specifically, according to the apparent reflectivity measurement value of the second observation band combination actually obtained by the passive multispectral satellite remote sensing platform, the closest theoretical value of the apparent reflectivity and its corresponding aerosol layer height and aerosol optical thickness can be indexed from the aerosol inversion lookup table, and the inverted aerosol optical layer height and aerosol optical thickness corresponding to the apparent reflectivity measurement value of the second observation band combination can be determined through interpolation calculation and the like.

[0079] In summary, please refer to Figure 3 As shown, the present invention provides an aerosol layer height and optical thickness inversion method, which can remove the influence of aerosol and atmospheric molecules in the aerosol layer height and aerosol optical thickness inversion by obtaining the polarization reflectivity measurement value of the first observation band combination. On the other hand, by combining the aerosol vertical distribution profile inverted by GRASP and the observation geometry information, multiple pairs of prior aerosol layer heights and prior aerosol optical thicknesses are set, and the atmospheric radiation transmission model is used to determine the apparent reflectivity corresponding to each pair of aerosol layer height and aerosol optical thickness, and an aerosol inversion lookup table is established based on this. The apparent reflectance measurement value of the second observation band combination actually obtained by the passive multispectral satellite remote sensing platform is obtained, the apparent reflectance measurement value of the second observation band is compared and verified with the aerosol inversion lookup table, and the aerosol inversion lookup table is interpolated and fitted to obtain the inverted aerosol layer height and inverted aerosol optical thickness corresponding to the apparent reflectance measurement value of the second observation band combination, thereby achieving fast and accurate inversion of the aerosol layer height and aerosol optical thickness.

[0080] Based on the same inventive concept, please refer to Figure 4 As shown, another embodiment of the present invention further provides an aerosol layer height and optical thickness inversion device 11, the device comprising:

[0081] The data acquisition module 111 is used to obtain the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, polarization reflectivity measurement values ​​of the first observation band combination, and reflectivity measurement values ​​of the second observation band combination in the observation area;

[0082] The data processing module 112 is used to obtain a priori aerosol layer height, a priori aerosol optical thickness and a priori surface reflectivity according to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information and the polarization reflectivity measurement value of the first observation band combination, and obtain an aerosol inversion lookup table according to the priori aerosol layer height, the priori aerosol optical thickness and the prior surface reflectivity;

[0083] The aerosol inversion module 113 is used to obtain the inverted aerosol layer height and the inverted aerosol optical thickness according to the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table.

[0084] It should be noted that the aerosol layer height and optical thickness inversion device includes the aerosol layer height and optical thickness inversion method described in any of the above embodiments. Since the aerosol layer height and optical thickness inversion device provided in this embodiment belongs to the same inventive concept as the aerosol layer height and optical thickness inversion method provided in any of the above embodiments, it has at least the same beneficial effects as the aerosol layer height and optical thickness inversion method, and will not be described one by one here.

[0085] Based on the same inventive concept, please refer to Figure 5 As shown, another embodiment of the present invention further provides an electronic device 1, which may include a memory 12, a processor 13 and a bus, and may also include a computer program stored in the memory 12 and executable on the processor 13, such as an aerosol layer height and optical thickness inversion program.

[0086] Among them, the memory 12 includes at least one type of readable storage medium, and the readable storage medium includes flash memory, mobile hard disk, multimedia card, card-type memory (for example: SD or DX memory, etc.), magnetic memory, disk, optical disk, etc. In some embodiments, the memory 12 can be an internal storage unit of the electronic device 1, such as a mobile hard disk of the electronic device 1. In other embodiments, the memory 12 can also be an external storage device of the electronic device 1, such as a plug-in mobile hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), etc. equipped on the electronic device 1. Further, the memory 12 can also include both an internal storage unit of the electronic device 1 and an external storage device. The memory 12 can not only be used to store application software and various types of data installed in the electronic device 1, such as the code for inversion of aerosol layer height and optical thickness, but also can be used to temporarily store data that has been output or is to be output.

[0087] In some embodiments, the processor 13 may be composed of an integrated circuit, for example, a single packaged integrated circuit, or a plurality of packaged integrated circuits with the same or different functions, including one or more central processing units (CPUs), microprocessors, digital processing chips, graphics processors, and combinations of various control chips. The processor 13 is the control core (Control Unit) of the electronic device 1, and uses various interfaces and lines to connect various components of the entire electronic device 1, and executes or executes programs or modules (such as aerosol layer height and optical thickness inversion programs, etc.) stored in the memory 12, and calls data stored in the memory 12 to execute various functions of the electronic device 1 and process data.

[0088] The processor 13 executes the operating system and various installed application programs of the electronic device 1. The processor 13 executes the application programs to implement the steps in the above-mentioned aerosol layer height and optical thickness inversion method.

[0089] Exemplarily, the computer program may be divided into one or more modules, which are stored in the memory 12 and executed by the processor 13 to complete the present application. The one or more modules may be a series of computer program instruction segments capable of completing specific functions, which are used to describe the execution process of the computer program in the electronic device 1. For example, the computer program may be divided into a data acquisition module 111, a data processing module 112, and an aerosol inversion module 113.

[0090] The above-mentioned integrated unit implemented in the form of a software function module can be stored in a computer-readable storage medium, and the computer-readable storage medium can be non-volatile or volatile. The above-mentioned software function module is stored in a storage medium, including a number of instructions for enabling a computer device (which can be a personal computer, a computer device, or a network device, etc.) or a processor to perform part of the functions of the aerosol layer height and optical thickness inversion method described in each embodiment of the present application.

[0091] The above embodiments are merely illustrative of the principles and effects of the present invention, and are not intended to limit the present invention. Anyone familiar with the art may modify or alter the above embodiments without departing from the spirit and scope of the present invention. Therefore, all equivalent modifications or alterations made by a person of ordinary skill in the art without departing from the spirit and technical concept disclosed by the present invention shall still be covered by the claims of the present invention.

Claims

1. A method for inverting aerosol layer height and optical thickness, characterized in that: include: Obtain the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, polarization reflectivity measurement values ​​of the first observation band combination, and reflectivity measurement values ​​of the second observation band combination in the observation area; According to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, and the polarization reflectivity measurement value of the first observation band combination, a priori aerosol layer height, a priori aerosol optical thickness, and a priori surface reflectivity are obtained; According to the a priori aerosol layer height, the a priori aerosol optical thickness and the a priori surface reflectivity, an aerosol inversion lookup table is obtained; The inverted aerosol layer height and the inverted aerosol optical thickness are obtained according to the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table.

2. The aerosol layer height and optical thickness inversion method according to claim 1, characterized in that: The first observation band combination includes observation bands that are greatly affected by aerosols and atmospheric molecules and observation bands that are less affected by aerosols and atmospheric molecules, and the second observation bands include 688nm and 680nm, 764nm and 780nm.

3. The aerosol layer height and optical thickness inversion method according to claim 1, characterized in that: The method of obtaining the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, polarization reflectivity measurement values ​​of the first observation band combination, and reflectivity measurement values ​​of the second observation band combination of the observation area includes: Determine the atmospheric profile type and aerosol microphysics and single scattering characteristics based on the observation area and observation time of the passive multispectral satellite remote sensing platform; Acquiring observation geometry information and observation latitude and longitude information collected by a positioning and attitude determination system carried by the passive multispectral satellite remote sensing platform; The polarization reflectivity measurement value of the first observation band combination and the reflectivity measurement value of the second observation band combination collected by the atmospheric correction instrument carried by the passive multispectral satellite remote sensing platform are obtained.

4. The aerosol layer height and optical thickness inversion method according to claim 1, characterized in that: The method of obtaining a priori aerosol layer height, a priori aerosol optical thickness and a priori surface reflectivity based on the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information and the polarization reflectivity measurement value of the first observation band combination, includes: According to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information and the polarization reflectivity measurement value of the first observation band combination, the GRASP program is used to obtain the prior aerosol layer height, the prior aerosol optical thickness and the prior surface reflectivity.

5. The aerosol layer height and optical thickness inversion method according to claim 1, characterized in that: The step of obtaining an aerosol inversion lookup table according to the prior aerosol layer height, the prior aerosol optical thickness and the prior surface reflectivity includes: According to the a priori aerosol layer height, the a priori aerosol optical thickness and the a priori surface reflectivity, using an atmospheric radiation transfer model, obtaining theoretical values ​​of apparent reflectivity of the observation area corresponding to all observation bands in the first observation band combination; An aerosol inversion lookup table is established according to the prior aerosol layer height, the prior aerosol optical thickness, the prior surface reflectivity and the theoretical value of the apparent reflectivity.

6. The aerosol layer height and optical thickness inversion method according to claim 5, characterized in that: The step of establishing an aerosol inversion lookup table according to the prior aerosol layer height, the prior aerosol optical thickness, the prior surface reflectivity and the theoretical value of the apparent reflectivity includes: Construct an M×N two-dimensional table, in which the horizontal rows represent the aerosol layer height and the vertical rows represent the aerosol optical thickness. Fill the corresponding theoretical values ​​of the apparent reflectivity into the two-dimensional table to obtain the aerosol inversion lookup table, in which each grid point represents the mapping relationship between the aerosol layer height, the aerosol optical thickness and the apparent reflectivity.

7. The aerosol layer height and optical thickness inversion method according to claim 1, characterized in that: The method of obtaining the inverted aerosol layer height and the inverted aerosol optical thickness based on the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table includes: Obtaining the apparent reflectivity of the second observation band combination according to the reflectivity measurement value of the second observation band combination; The apparent reflectivity of the second observation band combination is compared and verified with the aerosol inversion lookup table, and the inverted aerosol layer height and the inverted aerosol optical thickness are obtained by interpolating and fitting the aerosol inversion lookup table.

8. An aerosol layer height and optical thickness inversion device, characterized in that: The device comprises: A data acquisition module is used to obtain the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation longitude and latitude information, polarization reflectivity measurement values ​​of the first observation band combination, and reflectivity measurement values ​​of the second observation band combination in the observation area; A data processing module, for obtaining a priori aerosol layer height, a priori aerosol optical thickness and a priori surface reflectivity according to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information and the polarization reflectivity measurement value of the first observation band combination, and obtaining an aerosol inversion lookup table according to the priori aerosol layer height, the priori aerosol optical thickness and the prior surface reflectivity; The aerosol inversion module is used to obtain the inverted aerosol layer height and the inverted aerosol optical thickness according to the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table.

9. An electronic device, characterized in that: The electronic device comprises: one or more processors; A storage device for storing one or more programs, which, when executed by the one or more processors, enables the electronic device to implement the aerosol layer height and optical thickness inversion method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: A computer program is stored thereon, and when the computer program is executed by a processor of a computer, the computer is enabled to execute the aerosol layer height and optical thickness inversion method as described in any one of claims 1 to 7.

Images