An aerosol layer height and optical thickness inversion method, device, equipment and medium

Abstract

The application discloses an aerosol layer height and optical thickness inversion method, device, equipment and medium, the method obtains the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information, the polarization reflectivity measurement value of the first observation waveband combination and the reflectivity measurement value of the second observation waveband combination of an observation area, and establishes an aerosol inversion lookup table 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 waveband combination, and obtains the inversion aerosol layer height and the inversion aerosol optical thickness according to the reflectivity measurement value of the second observation waveband combination and the aerosol inversion lookup table. The method greatly makes up for the defects of narrow aerosol layer height product width and long revisit period faced by active remote sensing by using passive remote sensing data, and improves the precision of measuring the aerosol layer height.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of remote sensing image processing, and particularly relates to an aerosol layer height and optical thickness inversion method, device, equipment and medium. BACKGROUND

[0002] In satellite imaging remote sensing from visible to short-wave infrared band, atmospheric scattering and multiple coupling scattering between atmosphere and earth surface have a blurring effect on remote sensing images, which reduces the definition of remote sensing images and makes the apparent reflectance of the detected atmospheric layer top deviate from the real earth surface reflectance. Therefore, after collecting the original image of the earth surface using a satellite remote sensing platform, atmospheric correction is usually needed to be performed on the original image to obtain more accurate image content and to acquire the aerosol layer height distribution, i.e. vertical distribution, which is of great significance to atmospheric correction research.

[0003] In order to measure the aerosol layer height distribution, since the cloud-aerosol polarization lidar (CALIOP) carried on the CALIPSO platform began to provide data in 2006, there has been a lot of research on aerosol layer height measurement. The active remote sensing technology using high spectral resolution lidar (LIDAR) can accurately measure the vertical distribution of aerosols, but there are limitations in width and revisit period. So far, passive satellites still cannot provide aerosol vertical distribution data with the same accuracy as LIDAR, but in terms of spatial range and revisit period, the data of passive satellites are indispensable supplements to LIDAR aerosol layer height inversion.

[0004] The O2A band (0.764 μm) and O2B band (0.688 μm) measured by passive satellites have the advantage of being not affected by thermal radiation, which has been proved to have great value in the inversion of aerosol layer height in the case of large aerosol optical thickness. The algorithms for retrieving aerosol vertical distribution using oxygen absorption spectrum can be divided into two categories. The first category directly applies the principle that the reflectivity in the oxygen absorption band increases with the rise of the aerosol scattering layer. This kind of algorithm derives the aerosol height according to the reflectivity ratio of the channels inside and outside the absorption band. However, this method can only produce single aerosol height information. The second category of methods, which act on the spectral fitting technology of satellites, the observation instruments usually have relatively coarse spatial resolution (for example, GOME, GOME-2 and SCIAMCHY) or limited spatial coverage (for example, GOSAT, OCO-2 and TanSat). In contrast, the reflectivity ratio method provides less detailed aerosol profile than the spectral fitting algorithm, but the corresponding instruments usually have better spatial resolution and wider spatial coverage (for example, POLDER, MERIS and EPIC). Accurate height retrieval depends on proper assumptions of aerosol optical properties and surface reflectivity. At present, there is an urgent need for passive satellite measurement technology in the context of atmospheric correction to accurately retrieve the vertical distribution of aerosol layers. SUMMARY

[0005] The present application aims to provide an aerosol layer height and optical thickness inversion method, device, equipment and medium, which realizes rapid and accurate inversion of aerosol layer height and aerosol optical thickness based on optimal estimation technology and radiation transfer model of passive multispectral satellite data, and solves the technical problems that the prior art highly depends on prior assumptions of surface reflectivity and is difficult to obtain accurate aerosol layer height and aerosol optical thickness results under the target of atmospheric correction.

[0006] To solve the above technical problems, the present application is realized by the following technical scheme:

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

[0008] Obtaining the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information, polarization reflectance measurement value of the first observation wavelength combination and reflectance measurement value of the second observation wavelength combination of the observation area;

[0009] According to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information and polarization reflectance measurement value of the first observation wavelength combination, obtaining the prior aerosol layer height, prior aerosol optical thickness and prior surface reflectivity;

[0010] According to the prior aerosol layer height, prior aerosol optical thickness and prior surface reflectivity, obtaining the aerosol inversion lookup table;

[0011] According to the reflectance measurement value of the second observation wavelength combination and the aerosol inversion lookup table, obtaining the inversion aerosol layer height and inversion aerosol optical thickness.

[0012] In an embodiment of the present application, the first observation wavelength combination comprises observation wavelengths greatly affected by aerosols and atmospheric molecules and observation wavelengths less affected by aerosols and atmospheric molecules, and the second observation wavelength comprises 688nm and 680nm, 764nm and 780nm.

[0013] In an embodiment of the present application, the obtaining of the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information, polarization reflectance measurement value of the first observation wavelength combination and reflectance measurement value of the second observation wavelength combination of the observation area comprises:

[0014] According to the observation area and observation time of the passive multispectral satellite remote sensing platform, determining the atmospheric profile type and aerosol microphysics and single scattering characteristics;

[0015] Obtaining 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] Obtaining polarization reflectance measurement values of a first observation waveband combination and reflectance measurement values of a second observation waveband combination collected by an atmospheric correction instrument carried by the passive multispectral satellite remote sensing platform.

[0017] In an embodiment of the present application, the obtaining of prior aerosol layer height, prior aerosol optical thickness and prior surface reflectance according to the atmospheric profile type, aerosol microphysical and single scattering characteristics, observation geometry information, observation latitude and longitude information and polarization reflectance measurement values of the first observation waveband combination comprises:

[0018] The prior aerosol layer height, prior aerosol optical thickness and prior surface reflectance are obtained by using a GRASP program according to the atmospheric profile type, aerosol microphysical and single scattering characteristics, observation geometry information, observation latitude and longitude information and polarization reflectance measurement values of the first observation waveband combination.

[0019] In an embodiment of the present application, the obtaining of the aerosol retrieval lookup table according to the prior aerosol layer height, prior aerosol optical thickness and prior surface reflectance comprises:

[0020] The apparent reflectance theoretical values corresponding to all observation wavebands in the first observation waveband combination of the observation area are obtained by using an atmospheric radiation transfer model according to the prior aerosol layer height, prior aerosol optical thickness and prior surface reflectance.

[0021] The aerosol retrieval lookup table is established according to the prior aerosol layer height, prior aerosol optical thickness, prior surface reflectance and the apparent reflectance theoretical values.

[0022] In an embodiment of the present application, the establishing of the aerosol retrieval lookup table according to the prior aerosol layer height, prior aerosol optical thickness, prior surface reflectance and the apparent reflectance theoretical values comprises:

[0023] A two-dimensional table of MxN is constructed, the horizontal line represents the aerosol layer height, the vertical line represents the aerosol optical thickness, and the corresponding apparent reflectance theoretical values are filled into the two-dimensional table to obtain the aerosol retrieval lookup table, wherein each grid point represents the mapping relationship of the aerosol layer height, aerosol optical thickness and apparent reflectance.

[0024] In an embodiment of the present application, the obtaining of the retrieval aerosol layer height and retrieval aerosol optical thickness according to the reflectance measurement values of the second observation waveband combination and the aerosol retrieval lookup table comprises:

[0025] According to the reflectivity measurement value of the second observation waveband combination, an apparent reflectivity of the second observation waveband combination is obtained;

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

[0027] Based on the same inventive concept, another embodiment of the present application also provides an aerosol layer height and optical thickness inversion device, which comprises:

[0028] The data acquisition module is configured to acquire the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information, the polarization reflectivity measurement value of the first observation waveband combination, and the reflectivity measurement value of the second observation waveband combination of an observation area;

[0029] The data processing module is configured to obtain the prior aerosol layer height, the prior aerosol optical thickness, and the prior ground 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 waveband combination, and obtain the aerosol inversion lookup table according to the prior aerosol layer height, the prior aerosol optical thickness, and the prior ground surface reflectivity.

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

[0031] Based on the same inventive concept, another embodiment of the present application also provides an electronic device, which comprises:

[0032] One or more processors;

[0033] A storage device is configured to store one or more programs, which, when executed by the one or more processors, cause the electronic device to implement the aerosol layer height and optical thickness inversion method of any of the above-mentioned embodiments.

[0034] Based on the same inventive concept, another embodiment of the present application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor of a computer, causes the computer to perform the aerosol layer height and optical thickness inversion method of any of the above-mentioned embodiments.

[0035] The aerosol layer height and optical thickness inversion method provided by the present application has the following beneficial effects: the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information, polarization reflectivity measurement values of a first observation waveband combination, and reflectivity measurement values of a second observation waveband combination of an observation area are obtained, prior aerosol layer height, prior aerosol optical thickness and prior ground reflectivity are obtained 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 waveband combination, an aerosol inversion lookup table is obtained according to the prior aerosol layer height, prior aerosol optical thickness and prior ground reflectivity, and inversion aerosol layer height and inversion aerosol optical thickness are obtained according to the reflectivity measurement values of the second observation waveband combination and the aerosol inversion lookup table. The aerosol layer height and optical thickness inversion method can obtain aerosol layer height data with a larger width and a shorter revisit period relative to LIDAR by a passive multispectral satellite, which is a major breakthrough in the field of satellite inversion of atmospheric components. The aerosol layer height parameters in the upper space of the observation area can be obtained by using the multispectral data of the observation area collected by the passive multispectral satellite and combining the forward radiation transfer model, and the defects of narrow aerosol layer height product width and long revisit period faced by active remote sensing (CALIOP, LIDAR) are greatly made up by using passive remote sensing data. Of course, any product implementing the present application does not necessarily need to achieve all the advantages described above. BRIEF DESCRIPTION OF DRAWINGS

[0036] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings needed to be used in the description of the embodiments. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can also be obtained by those skilled in the art without any creative effort on the basis of these drawings.

[0037] Figure 1 A flowchart of an aerosol layer height and optical thickness inversion method provided by an exemplary embodiment of the present application.

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

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

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

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

[0042] The present application is described below by way of specific embodiments, and those skilled in the art can easily understand other advantages and effects of the present application from the disclosure. The present application can also be implemented or applied in other different embodiments, and various modifications or changes can be made to the details based on different views and applications without departing from the spirit of the present application. 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 diagrams provided in the following embodiments only illustrate the basic concept of the present application in a schematic manner, and only the components related to the present application are shown in the diagrams, not the number, shape and size of the components when actually implemented. The actual implementation of each component may be a random change, and the component layout pattern may be more complex.

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

[0045] In order to solve the technical problems that the prior art highly depends on the prior assumption of surface reflectivity and it is difficult to obtain accurate aerosol layer height and aerosol optical thickness results under the target of atmospheric correction, the present application provides an aerosol layer height and optical thickness inversion method, which realizes the inversion of aerosol layer height and aerosol optical thickness based on passive multispectral satellite data by using optimal estimation technology and a radiation transfer model. Please refer to Figure 1 As shown, the aerosol layer height and optical thickness inversion method comprises the following steps:

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

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

[0048] S300: obtaining an aerosol inversion look-up table according to the prior aerosol layer height, the prior aerosol optical thickness and the prior ground reflectance;

[0049] S400: obtaining an inversion aerosol layer height and an inversion aerosol optical thickness according to the reflectance measurement value of the second observation band combination and the aerosol inversion look-up table.

[0050] The steps of the aerosol layer height and optical thickness inversion method are described in detail as follows.

[0051] It should be noted that in the embodiment, the aerosol layer height and optical thickness inversion method is applied to a passive multispectral satellite remote sensing platform. The passive multispectral satellite remote sensing platform can synchronously carry an atmospheric correction instrument, a position and orientation system, a controller and a data processor. The working bands of the atmospheric correction instrument include 443 nm, 490 nm, 565 nm, 670 nm, 870 nm, 910 nm, 680 nm, 688 nm, 764 nm, 730 nm, 1380 nm, 1500 nm, 1560 nm and 1610 nm. The atmospheric correction instrument adopts a three-channel scheme, that is, three different detection directions are set in each spectral band, and the detection directions are usually set as 0°, 60° and 120° to cover as much polarization information as possible, so that a wide range of atmospheric multispectral polarization data can be obtained. The polarization observation bands can include 680 nm, 688 nm, 764 nm and 780 nm, and the remaining bands use a non-polarization observation mode to observe reflectance. The position and orientation system (POS) can obtain geographic information, flight height and other flight trajectory data of the aerial remote sensing platform on the flight path, and obtain attitude data of the atmospheric correction instrument. The controller can control the atmospheric correction instrument, the position and orientation system and the like. The data processor can be used to execute the aerosol layer height and optical thickness inversion method.

[0052] Firstly, step S100 is executed, that is, the atmospheric profile type, aerosol microphysics and single scattering characteristics of the observation area, observation geometric information, observation latitude and longitude information, polarization reflectance measurement value of the first observation band combination and reflectance measurement value of the second observation band combination are obtained.

[0053] In an example embodiment of the present application, in step S100, the atmospheric profile type, aerosol microphysics and single scattering characteristics of the observation area, observation geometric information, observation latitude and longitude information, polarization reflectance measurement value of the first observation band combination and reflectance measurement value of the second observation band combination are further obtained by the following steps:

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

[0055] S120: obtaining 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;

[0056] S130: obtaining polarization reflectance measurement values of a first observation waveband combination and reflectance measurement values of a second observation waveband combination collected by an atmospheric correction instrument carried by the passive multispectral satellite remote sensing platform.

[0057] Specifically, the positioning and attitude determination system carried by the passive multispectral satellite remote sensing platform collects observation geometry information and observation latitude and longitude, the observation geometry information including the relative position and angular relationship among the satellite, the observation area and the sun, in this embodiment, the observation geometry information including satellite attitude, observation angle and height information, the satellite attitude including the roll angle, pitch angle and heading angle of the satellite, the observation angle including the incident angle, observation angle and azimuth angle, and the height information including the height of the satellite relative to the earth's surface, and the observation latitude and longitude including the geographical position of the observation area. The multispectral high-resolution camera and the atmospheric correction instrument carried by the passive multispectral satellite remote sensing platform can synchronously collect the ground surface image data and the polarization reflectance data of the first observation waveband combination and the reflectance data of the second observation waveband combination under the control of the controller. After the data processor obtains the raw data collected by both, it can perform preprocessing such as radiation calibration, polarization calibration and geometry calibration to obtain the ground surface image, the polarization reflectance measurement values corresponding to the first observation waveband combination and the reflectance measurement values corresponding to the second observation waveband combination which can be used for subsequent processing.

[0058] In an exemplary embodiment of the present application, the first observation waveband combination includes observation wavebands greatly affected by aerosols and atmospheric molecules and observation wavebands less affected by aerosols and atmospheric molecules, and the second observation waveband includes 688nm and 680nm, 764nm and 780nm.

[0059] Specifically, the first observation band combination includes 443 nm, 490 nm, 565 nm, 670 nm, 870 nm, 910 nm, 680 nm, 688 nm, 764 nm, 730 nm, 1380 nm, 1500 nm, 1560 nm, and 1610 nm. It should be noted that the longer the wavelength, the less the influence of atmospheric molecules and aerosol scattering on the observation band, and the first observation band combination in the embodiment includes observation bands with large influence of aerosol and atmospheric molecules and observation bands with small influence of aerosol and atmospheric molecules. It should be noted that in the embodiment, the polarized reflectance measurement value of the atmospheric correction instrument polarized observation band can be used as the input of the GRASP program, and the prior aerosol layer height, the prior aerosol optical thickness, and the prior surface reflectance are calculated by the GRASP program during subsequent aerosol inversion.

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

[0061] Then, step S200 is performed to obtain the prior aerosol layer height, the prior aerosol optical thickness, and the prior surface reflectance according to the atmospheric profile type, the aerosol microphysics and single scattering characteristics, the observation geometry information, the observation latitude and longitude information, and the polarized reflectance measurement value of the first observation band combination.

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

[0063] Specifically, the atmospheric correction instrument obtains corresponding polarized reflectance measurement values of each observation band in the first observation band combination in the polarized observation mode, and determines an aerosol layer height initial value, an aerosol optical thickness initial value and a ground reflectance initial value of the observation area according to the atmospheric profile type, aerosol microphysics and single scattering characteristics, observation geometry information, observation latitude and longitude information and the polarized reflectance measurement values of the first observation band combination by using the GRASP program, takes the aerosol layer height initial value corresponding to each band in the first observation band combination as the prior aerosol layer height, takes the aerosol optical thickness initial value corresponding to each band in the first observation band combination as the prior aerosol optical thickness, and takes the ground reflectance initial value as the prior ground reflectance of the observation area.

[0064] Then, step S300 is performed, that is, the aerosol retrieval lookup table is obtained according to the prior aerosol layer height, the prior aerosol optical thickness and the prior ground reflectance.

[0065] In an example embodiment of the present application, in step S300, the aerosol retrieval lookup table is obtained according to the prior aerosol layer height, the prior aerosol optical thickness and the prior ground reflectance, and further includes the following steps:

[0066] S310: The apparent reflectance theoretical value corresponding to all observation bands in the first observation band combination of the observation area is obtained by using the atmospheric radiation transfer model according to the prior aerosol layer height, the prior aerosol optical thickness and the prior ground reflectance;

[0067] S320: The aerosol retrieval lookup table is established according to the prior aerosol layer height, the prior aerosol optical thickness, the prior ground reflectance and the apparent reflectance theoretical value.

[0068] In an example embodiment of the present application, step S320 includes: constructing a two-dimensional table of MxN, the horizontal line representing the aerosol layer height and the vertical line representing the aerosol optical thickness, filling the corresponding apparent reflectance theoretical value into the two-dimensional table to obtain the aerosol retrieval lookup table, wherein each grid point represents the mapping relationship of the aerosol layer height, the aerosol optical thickness and the apparent reflectance.

[0069] Specifically, in the present embodiment, please refer to Figure 2As shown, the atmospheric radiation transfer model can be the Unified Linearized Vector Radiative Transfer Model (UNL-VRTM). By inputting multiple pairs of prior aerosol layer heights, prior aerosol optical depths, and observation geometry and latitude and longitude information from a passive multispectral 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. Based on this, a lookup table of aerosol layer height and aerosol optical depth is constructed. 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 depth. The corresponding theoretical values ​​of apparent reflectivity are filled in 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 depth 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 height and aerosol optical thickness, observation altitude and observation geometry information, the atmospheric radiation transfer model is used to obtain the theoretical value of the apparent reflectivity corresponding to the oxygen absorption band and the non-oxygen absorption band of each combination, and an aerosol inversion lookup table is constructed.

[0074] Finally, step 400 is executed, that is, the inverted aerosol layer height and the inverted aerosol optical depth 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, obtaining the inverted aerosol layer height and the inverted aerosol optical depth based on the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table further includes 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: comparing and checking the apparent reflectivity of the second observation band combination with the aerosol inversion lookup table, and obtaining the inversion aerosol layer height and the inversion 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 apparent reflectivity theoretical value and the corresponding aerosol layer height and aerosol optical thickness can be obtained from the aerosol inversion lookup table, and the inversion aerosol optical layer height and the inversion aerosol optical thickness corresponding to the apparent reflectivity measurement value of the second observation band combination can be determined by interpolation calculation and the like.

[0079] In summary, referring to Figure 3 The aerosol layer height and optical thickness inversion method provided by the application can remove the influence of aerosols and atmospheric molecules in the inversion of the aerosol layer height and the aerosol optical thickness 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, a plurality of pairs of prior aerosol layer height and prior aerosol optical thickness are set, and the apparent reflectivity corresponding to each pair of aerosol layer height and aerosol optical thickness is determined by using an atmospheric radiation transfer model, and an aerosol inversion lookup table is established. The apparent reflectivity measurement value of the second observation band combination actually obtained by the passive multispectral satellite remote sensing platform is obtained, the apparent reflectivity measurement value of the second observation band is compared and checked with the aerosol inversion lookup table, and the inversion aerosol layer height and the inversion aerosol optical thickness corresponding to the apparent reflectivity measurement value of the second observation band combination are obtained by interpolating and fitting the aerosol inversion lookup table, so that the rapid and accurate inversion of the aerosol layer height and the aerosol optical thickness is realized.

[0080] Based on the same inventive concept, referring to Figure 4 Another embodiment of the application further provides an aerosol layer height and optical thickness inversion device 11, which comprises:

[0081] The data acquisition module 111 is configured to acquire the atmospheric profile type, the aerosol microphysics and single scattering characteristics, the observation geometry information, the observation latitude and longitude information, the polarization reflectivity measurement value of the first observation band combination and the reflectivity measurement value of the second observation band combination of an observation area.

[0082] The data processing module 112 is configured to obtain the prior aerosol layer height, the prior aerosol optical thickness and the prior ground reflectivity according to the atmospheric profile type, the aerosol microphysics and single scattering characteristics, the observation geometry information, the observation latitude and longitude information and the polarization reflectivity measurement value of the first observation waveband combination, and obtain the aerosol retrieval lookup table according to the prior aerosol layer height, the prior aerosol optical thickness and the prior ground reflectivity.

[0083] The aerosol retrieval module 113 is configured to obtain the retrieval aerosol layer height and the retrieval aerosol optical thickness according to the reflectivity measurement value of the second observation waveband combination and the aerosol retrieval lookup table.

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

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

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

[0087] The processor 13 may, in some embodiments, be composed of integrated circuits, for example, can be composed of a single packaged integrated circuit, or can be composed of multiple packaged integrated circuits of the same function or different functions, including one or more combinations of central processing units (CPUs), microprocessors, digital processing chips, graphics processors, and various control chips. The processor 13 is the control core of the electronic device 1, and connects various components of the entire electronic device 1 through various interfaces and lines, executes programs or modules stored in the memory 12 (such as the aerosol layer height and optical thickness inversion program), 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 of the electronic device 1 and various installed application programs. The processor 13 executes the application programs to implement the steps in the above-mentioned aerosol layer height and optical thickness inversion method.

[0089] For example, the computer program can 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 can be a series of computer program instruction segments capable of completing a specific function, which are used to describe the execution process of the computer program in the electronic device 1. For example, the computer program can be divided into a data acquisition module 111, a data processing module 112, and an aerosol inversion module 113.

[0090] The integrated units implemented in the form of software function modules described above can be stored in a computer readable storage medium, which can be non-volatile or volatile. The software function modules described above are stored in a storage medium, including a plurality of instructions for causing a computer device (which can be a personal computer, a computer device, or a network device, etc.) or a processor to execute part of the functions of the aerosol layer height and optical thickness inversion method described in various embodiments of the present application.

[0091] The above-described embodiments are only illustrative of the principles and effects of the present application, and are not intended to limit the present application. Any person skilled in the art can modify or change the above-described embodiments without departing from the spirit and scope of the present application. Therefore, all equivalent modifications or changes made by those skilled in the art without departing from the spirit and technical idea disclosed in the present application should be covered by the claims of the present application.

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 latitude and longitude information, polarization reflectance measurement values ​​of the first observation band combination, and reflectance 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 the polarization reflectance 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 reflectance; According to the a priori aerosol layer height, the a priori aerosol optical depth and the a priori surface reflectivity, using an atmospheric radiation transfer model, obtain theoretical values ​​of apparent reflectivity of the observation area corresponding to all observation bands in the first observation band combination; Build a A two-dimensional table is prepared, where 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 entered into the two-dimensional table to obtain an aerosol inversion lookup table, where each grid point represents the mapping relationship between aerosol layer height, aerosol optical thickness, and apparent reflectivity. The inverted aerosol layer height and the inverted aerosol optical depth 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 depth 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 obtaining of 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 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; Obtaining 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 depth inversion method according to claim 1, characterized in that: The method further comprises: obtaining an inverted aerosol layer height and an inverted aerosol optical depth based on the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table; Obtaining an apparent reflectivity of the second observation waveband combination according to the reflectivity measurement value of the second observation waveband 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 inverted aerosol optical thickness are obtained by interpolating and fitting the aerosol inversion lookup table.

5. 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 latitude and longitude information, polarization reflectance measurement values ​​of the first observation band combination, and reflectance measurement values ​​of the second observation band combination in the observation area; The data processing module is used 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 latitude and longitude information and the polarization reflectivity measurement value of the first observation band combination using the GRASP program; based on the prior aerosol layer height, the prior aerosol optical thickness and the prior surface reflectivity, the atmospheric radiation transfer model is used to obtain the theoretical value of the apparent reflectivity of the observation area corresponding to all observation bands in the first observation band combination, and construct a A two-dimensional table is prepared, where 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 entered into the two-dimensional table to obtain an aerosol inversion lookup table, where each grid point represents the mapping relationship between aerosol layer height, aerosol optical thickness, and apparent reflectivity. The aerosol inversion module is used to obtain the inverted aerosol layer height and the inverted aerosol optical depth according to the reflectivity measurement value of the second observation band combination and the aerosol inversion lookup table.

6. 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 4.

7. 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 caused to execute the aerosol layer height and optical thickness inversion method according to any one of claims 1 to 4.

Images