A shortwave radiation full-factor integrated remote sensing estimation method and system

Abstract

The present invention relates to a shortwave radiation full-factor integrated remote sensing estimation method and system, wherein the method comprises: using the albedo of the acquired remote sensing observation data as the input of a pre-constructed simulation library, obtaining the radiation brightness of the blue, green and red bands at the top of the atmosphere or a certain height based on the atmospheric parameters in the simulation library; comparing the radiation brightness obtained by the simulation library with the radiation brightness of the acquired remote sensing observation data, and then matching the atmospheric parameters that meet the preset conditions in the simulation library according to the comparison result; and analyzing and obtaining the full elements of shortwave radiation according to the albedo and the atmospheric parameters that meet the preset conditions. The present invention proposes a new idea of ​​shortwave radiation full-factor integrated remote sensing estimation, which can realize the synchronous and high-precision estimation of 8 types of shortwave radiation through the radiation brightness data of simulation data and remote sensing data, and the present invention does not require atmospheric parameters such as cloud detection and aerosols, and is applicable to all weather and global scales.

Description

Technical Field

[0001] The present invention relates to the field of remote sensing technology, and in particular to a shortwave radiation full-factor integrated remote sensing estimation method and system. Background Art

[0002] Shortwave radiation (also known as solar radiation, 0.3-3μm) is the only energy source for the Earth system, driving the energy, water and material cycles of the entire Earth system. Its temporal and spatial changes are closely related to global warming. With the development of remote sensing technology, the use of remote sensing data to estimate the various shortwave radiation components has great advantages and potential. However, the current use of remote sensing to estimate shortwave radiation is mostly focused on the estimation of the total shortwave downward radiation on the surface. There are relatively few estimates of other components of shortwave radiation, especially direct solar radiation, diffuse PAR, ultraviolet radiation, etc. There are fewer studies on the estimation. Similarly, there are very few studies on the simultaneous estimation of all the above shortwave elements.

[0003] Most of the existing shortwave radiation estimation methods rely on input data such as cloud parameters, aerosols, and atmospheric water vapor. There is still much room for improvement in the estimated radiation accuracy and the number of radiation elements. Summary of the invention

[0004] 1. Technical issues to be resolved

[0005] In view of the above-mentioned shortcomings and deficiencies of the prior art, the present invention provides a shortwave radiation full-factor integrated remote sensing estimation method and system, which solves the technical problems that the existing shortwave radiation estimation scheme relies on too many conditions and the estimated radiation accuracy and the number of radiation elements obtained are both low.

[0006] (II) Technical solution

[0007] In order to achieve the above object, the main technical solutions adopted by the present invention include:

[0008] In a first aspect, an embodiment of the present invention provides a shortwave radiation all-factor integrated remote sensing estimation method, comprising:

[0009] The albedo of the acquired remote sensing observation data is used as the input of a pre-built simulation library, and the radiance of the blue, green and red bands at the top of the atmosphere or at a certain height is obtained based on the atmospheric parameters in the simulation library;

[0010] Comparing the radiance obtained through the simulation library with the radiance of the acquired remote sensing observation data, and then matching the atmospheric parameters that meet the preset conditions in the simulation library according to the comparison result;

[0011] Based on the albedo and the atmospheric parameters that meet the preset conditions, all elements of shortwave radiation are analyzed and obtained.

[0012] Optionally, the albedo of the acquired remote sensing observation data is used as input to the pre-built simulation library, and also includes:

[0013] According to preset atmospheric conditions and spectral response functions of different sensors, multiple types of simulated radiation components under different atmospheric conditions and observation angles and atmospheric parameters corresponding to each type of simulated radiation components are obtained through simulation;

[0014] Building a simulation library based on the multiple types of simulated radiation components and the atmospheric parameters corresponding to each type of simulated radiation component;

[0015] The atmospheric parameters include: atmospheric transmittance, path radiation, large balloon surface albedo, and comprehensive radiation items related to atmospheric transmittance and total solar radiation at the top of the atmosphere.

[0016] Optionally, the preset atmospheric conditions are conditions retrieved from a radiation transfer simulation condition library constructed based on representative parameters, and the representative parameters include: global atmospheric profiles, cloud parameter products, surface albedo, aerosols, sun angles, and satellite / UAV observation angles.

[0017] Optionally, the albedo of the acquired remote sensing observation data is used as an input of a pre-built simulation library, and the radiance of the blue, green and red bands at the top of the atmosphere or at a certain height is obtained based on the atmospheric parameters in the simulation library, including:

[0018] Obtain the albedo of the blue, green, red bands and the entire shortwave band of different surface objects through satellites;

[0019] Using the albedo as the input of the simulation library, the radiance of the blue, green and red bands at the top of the atmosphere or at a certain height is obtained through the radiance formula based on the atmospheric parameters in the simulation library;

[0020] Wherein, the radiation brightness formula is:

[0021]

[0022] In the formula, μ 0 is the cosine of the solar zenith angle; μ is the cosine of the observed zenith angle; φ is the relative azimuth; L(μ 0 ,μ,φ) is the radiance at the top of the atmosphere or at a certain height (blue, green or red band); L 0 (μ 0 ,μ,φ) is the atmospheric radiation; ρ is the albedo of the large balloon surface; r 1 It is the albedo of the blue, green and red bands of different surface objects; ψ is the comprehensive radiation term related to the atmospheric transmittance and the total solar radiation at the top of the atmosphere.

[0023] Optionally, comparing the radiance obtained through the simulation library with the radiance of the acquired remote sensing observation data, and then matching the atmospheric parameters that meet the preset conditions in the simulation library according to the comparison result includes:

[0024] Using the root mean square error as a comparison condition, the radiance obtained through the simulation library is compared with the radiance of the remote sensing observation data;

[0025] Based on the comparison result, the atmospheric condition with the smallest difference in radiation brightness from the remote sensing observation data is selected in the simulation library.

[0026] Optionally, according to the albedo and the atmospheric parameters satisfying the preset conditions, all elements of shortwave radiation obtained by analysis include:

[0027] According to the albedo and the screened atmospheric parameters of the atmospheric conditions, eight types of radiation components of shortwave radiation are obtained through a radiation component formula;

[0028] The radiation component formula is:

[0029]

[0030] In the formula, F(μ 0 ) is the downward radiation from the surface or the radiation from the top of the atmosphere; r 2 is the albedo of the entire shortwave band of different surface objects; F 0 (μ 0 ) is the downward radiation of the atmosphere or the upward radiation from the top of the atmosphere without considering the contribution of the surface; γ is the comprehensive radiation term related to the atmospheric transmittance.

[0031] Optionally, all elements of the shortwave radiation include: total downward shortwave radiation from the surface, direct solar radiation, sky diffuse radiation, net surface radiation, photosynthetically active radiation, direct or diffuse, ultraviolet radiation and net shortwave radiation at the top of the atmosphere.

[0032] In a second aspect, an embodiment of the present invention provides a shortwave radiation all-factor integrated remote sensing estimation system, including:

[0033] A radiation brightness simulation module is used to use the albedo of the acquired remote sensing observation data as the input of a pre-built simulation library, and obtain the radiation brightness of the blue, green and red bands at the top of the atmosphere or a certain height based on the atmospheric parameters in the simulation library;

[0034] A comparison and matching module, used to compare the radiance obtained through the simulation library with the radiance of the acquired remote sensing observation data, and then match the atmospheric parameters that meet the preset conditions in the simulation library according to the comparison result;

[0035] The full-factor obtaining module is used to analyze and obtain the full factors of shortwave radiation according to the albedo and the atmospheric parameters that meet the preset conditions.

[0036] (III) Beneficial effects

[0037] The beneficial effects of the present invention are as follows: the present invention takes into account the needs of different studies for shortwave radiation components, and proposes a new idea for shortwave radiation full-factor integrated remote sensing estimation. By analyzing the radiation brightness of simulation data and remote sensing data, the synchronous and high-precision estimation of eight types of shortwave radiation can be achieved. Moreover, the present technology does not require cloud detection, nor does it rely on atmospheric parameters such as aerosols, and is applicable to all-weather and global scales. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] Figure 1 A schematic diagram of a flow chart of a shortwave radiation all-factor integrated remote sensing estimation method provided by the present invention;

[0039] Figure 2 A schematic diagram of a specific flow chart before step S1 of a shortwave radiation all-factor integrated remote sensing estimation method provided by the present invention;

[0040] Figure 3 A specific flow chart of step S1 of a shortwave radiation all-factor integrated remote sensing estimation method provided by the present invention;

[0041] Figure 4 A specific flow chart of step S2 of a shortwave radiation all-factor integrated remote sensing estimation method provided by the present invention;

[0042] Figure 5 A schematic diagram of the overall process of a shortwave radiation all-factor integrated remote sensing estimation method provided by the present invention. DETAILED DESCRIPTION

[0043] In order to better explain the present invention and facilitate understanding, the present invention is described in detail below through specific implementation modes in conjunction with the accompanying drawings.

[0044] like Figure 1 As shown, an embodiment of the present invention proposes an integrated remote sensing estimation method for shortwave radiation of all elements, comprising: first, taking the albedo of the acquired remote sensing observation data as the input of a pre-constructed simulation library, and obtaining the radiation brightness of the blue, green and red bands at the top of the atmosphere or a certain height based on the atmospheric parameters in the simulation library; secondly, comparing the radiation brightness obtained by the simulation library with the radiation brightness of the acquired remote sensing observation data, and then matching the atmospheric parameters that meet the preset conditions in the simulation library according to the comparison result; finally, according to the albedo and the atmospheric parameters that meet the preset conditions, analyzing and obtaining all elements of shortwave radiation.

[0045] Taking into account the needs of different studies on shortwave radiation components, this paper proposes a new idea for shortwave radiation full-factor integrated remote sensing estimation. By analyzing the radiation brightness of simulation data and remote sensing data, the synchronous and high-precision estimation of 8 types of shortwave radiation can be achieved. Moreover, this technology does not require cloud detection, nor does it rely on atmospheric parameters such as aerosols, and is applicable to all-weather and global scales.

[0046] In order to better understand the above technical solution, exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. Although exemplary embodiments of the present invention are shown in the accompanying drawings, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments described herein. On the contrary, these embodiments are provided to enable a clearer and more thorough understanding of the present invention and to fully convey the scope of the present invention to those skilled in the art.

[0047] Specifically, an embodiment of the present invention provides a shortwave radiation all-factor integrated remote sensing estimation method, which includes:

[0048] S1. The albedo of the acquired remote sensing observation data is used as the input of a pre-built simulation library, and the radiation brightness of the blue, green and red bands at the top of the atmosphere or a certain height is obtained based on the atmospheric parameters in the simulation library.

[0049] like Figure 1 As shown, before step S1, it also includes:

[0050] F11. According to the preset atmospheric conditions and the spectral response functions of different sensors, a general radiation transfer model (for example, MODTRAN, SBDART, RTTOV, etc.) is used to simulate and obtain multiple types of simulated radiation components under different atmospheric conditions and observation angles and the atmospheric parameters corresponding to each type of the simulated radiation components; wherein the atmospheric parameters include: atmospheric transmittance, path radiation, albedo of the atmospheric balloon surface, and comprehensive radiation items related to atmospheric transmittance and total solar radiation at the top of the atmosphere.

[0051] F12. Construct a simulation library based on the multiple types of simulated radiation components and the atmospheric parameters corresponding to each type of simulated radiation component.

[0052] Furthermore, the preset atmospheric conditions are conditions retrieved from a radiation transfer simulation condition library constructed based on representative parameters, and the representative parameters include: global atmospheric profiles, cloud parameter products, surface albedo, aerosols, sun angles, and satellite / UAV observation angles.

[0053] like Figure 3 As shown, step S1 includes:

[0054] S11. Obtain the albedo of the blue, green, red bands and the entire shortwave band of different surface objects through satellites.

[0055] S12. Using the albedo as an input to a simulation library, and obtaining the radiation brightness of the blue, green and red bands at the top of the atmosphere or at a certain height through a radiation brightness formula based on the atmospheric parameters in the simulation library.

[0056] Wherein, the radiation brightness formula is:

[0057]

[0058] In the formula, μ 0 is the cosine of the solar zenith angle; μ is the cosine of the observed zenith angle; φ is the relative azimuth; L(μ 0 ,μ,φ) is the radiance at the top of the atmosphere or at a certain height (blue, green or red band); L 0 (μ 0 ,μ,φ) is the atmospheric radiation; ρ is the albedo of the large balloon surface; r 1 It is the albedo of the blue, green and red bands of different surface objects; ψ is the comprehensive radiation term related to the atmospheric transmittance and the total solar radiation at the top of the atmosphere.

[0059] S2. Compare the radiation brightness obtained through the simulation library with the radiation brightness of the acquired remote sensing observation data, and then match the atmospheric parameters that meet the preset conditions in the simulation library according to the comparison result.

[0060] like Figure 4 As shown, step S2 includes:

[0061] S21. Using the root mean square error as a comparison condition, the radiation brightness obtained through the simulation library is compared with the radiation brightness of the remote sensing observation data obtained.

[0062] S22. Based on the comparison result, select the atmospheric condition with the smallest difference in radiation brightness from the remote sensing observation data in the simulation library.

[0063] S3. Analyze and obtain all elements of shortwave radiation based on the albedo and the atmospheric parameters that meet the preset conditions.

[0064] Further, step S3 includes:

[0065] According to the albedo and the screened atmospheric parameters of the atmospheric conditions, eight types of radiation components of shortwave radiation are obtained through a radiation component formula;

[0066] The radiation component formula is:

[0067]

[0068] In the formula, F(μ 0 ) is the downward radiation from the surface or the radiation from the top of the atmosphere; r2 is the albedo of the entire shortwave band of different surface objects; F 0 (μ 0 ) is the downward radiation of the atmosphere or the upward radiation from the top of the atmosphere without considering the contribution of the surface; γ is the comprehensive radiation term related to the atmospheric transmittance.

[0069] It is worth mentioning that remote sensing observation data refers to all radiation brightness data obtained by multispectral or hyperspectral satellites and drones containing red, green and blue bands, such as MODIS, FY-3D, FY-4A / B, HY-1C / D, GF series, Landsat series, VIIRS, AVHRR, Himawari-8 / 9, MSG, GOES, etc. The full elements of shortwave radiation include: total downward shortwave radiation from the surface, direct solar radiation, sky diffuse radiation, net surface radiation, photosynthetically active radiation, direct or diffuse, ultraviolet radiation, and net shortwave radiation from the top of the atmosphere.

[0070] In addition, an embodiment of the present invention further provides a shortwave radiation all-factor integrated remote sensing estimation system, including:

[0071] The radiation brightness simulation module is used to use the albedo of the acquired remote sensing observation data as the input of a pre-built simulation library, and obtain the radiation brightness of the blue, green and red bands at the top of the atmosphere or a certain height based on the atmospheric parameters in the simulation library.

[0072] The comparison and matching module is used to compare the radiation brightness obtained through the simulation library with the radiation brightness of the acquired remote sensing observation data, and then match the atmospheric parameters that meet the preset conditions in the simulation library according to the comparison result.

[0073] The full-factor obtaining module is used to analyze and obtain the full factors of shortwave radiation according to the albedo and the atmospheric parameters that meet the preset conditions.

[0074] Since the system / device described in the above embodiments of the present invention is a system / device used to implement the method of the above embodiments of the present invention, a person skilled in the art can understand the specific structure and deformation of the system / device based on the method described in the above embodiments of the present invention, and thus will not be described in detail here. All systems / devices used in the method of the above embodiments of the present invention belong to the scope of protection of the present invention.

[0075] At the same time, an embodiment of the present invention also provides a shortwave radiation full-factor integrated remote sensing estimation device, including: at least one database; and a memory communicatively connected to the at least one database; wherein the memory stores instructions that can be executed by the at least one database, and the instructions are executed by the at least one database so that the at least one database can execute a shortwave radiation full-factor integrated remote sensing estimation method as described above.

[0076] And, a computer-readable medium having computer-executable instructions stored thereon, wherein the executable instructions, when executed by a processor, implement the shortwave radiation all-factor integrated remote sensing estimation method as described above.

[0077] In summary, the present invention provides a shortwave radiation all-factor integrated remote sensing estimation method and system, such as Figure 5 As shown, the overall process of the present invention is:

[0078] (1) Based on the global atmospheric profile, cloud parameter products, surface albedo, aerosol, sun angle, satellite / UAV observation angle, etc., a representative radiation transfer simulation condition library is constructed, denoted as simu_DB (see the table below), and the corresponding change interval and change step size are set for each parameter that the radiation transfer simulation depends on. Based on the above condition library, the spectral response functions of different sensors are input, and a large number of simulations are performed in combination with general radiation transfer models (for example, MODTRAN, SBDART, RTTOV, etc.), and 8 types of radiation components under different atmospheric conditions and observation angles are simulated. For each type of radiation component, the corresponding atmospheric transmittance, path radiation and other parameters are calculated, denoted as atm_par, and finally a representative simulation library is constructed, denoted as flux_DB. Among them, the simu_DB parameter settings refer to Table 1:

[0079]

[0080] (2) Use existing satellite products (e.g., MODIS products) to obtain the albedo (albedo_input) of the blue, green, and red bands of different surface objects and the entire shortwave band. Based on the following formula, with albedo as input, the remote sensing blue, green, and red band radiance at the top of the atmosphere or a certain height is calculated based on the atmospheric parameters (atm_par) in the simulation library (flux_DB).

[0081] (3) The radiance (red, green, and blue bands) obtained by the above simulation is compared with the radiance (red, green, and blue bands) of the actual remote sensing image. The root mean square error (RMSE) is used as the judgment criterion. The atmospheric condition with the smallest difference from the remote sensing observation radiance is found in flux_DB (denoted as atm_par_final). The atmospheric parameters of atm_par_final and the albedo of the surface (albedo_input) are used to obtain 8 types of shortwave radiation components.

[0082] Preferably, the above method is a general method. Different satellite sensors have different corresponding spectral response functions. Therefore, different satellites need to be simulated separately according to the above technical process to obtain corresponding atm_par and flux_DB.

[0083] Based on the above description, the present invention can realize the high-precision, synchronous and integrated estimation of 8 types of shortwave radiation components based on conventional remote sensing data (radiant brightness of red, green and blue bands), providing key technical guarantee for scientific research and applications related to shortwave radiation.

[0084] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, systems or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0085] The present invention is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems) and computer program products according to the embodiments of the present invention. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions.

[0086] It should be noted that in the claims, any reference numerals placed between brackets shall not be construed as limiting the claims. The word "comprising" does not exclude the presence of components or steps not listed in the claims. The word "a" or "an" preceding a component does not exclude the presence of a plurality of such components. The invention may be implemented by means of hardware comprising several different components and by means of a suitably programmed computer. In the claims enumerating several means, several of these means may be embodied by the same hardware. The use of the words first, second, third, etc., is for convenience of expression only and does not indicate any order. These words may be understood as part of the component name.

[0087] In addition, it should be noted that, in the description of this specification, the description of the terms "one embodiment", "some embodiments", "embodiment", "example", "specific example" or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, unless they are contradictory.

[0088] Although the preferred embodiments of the present invention have been described, those skilled in the art may make other changes and modifications to these embodiments after knowing the basic creative concept. Therefore, the claims should be interpreted as including the preferred embodiments and all changes and modifications falling within the scope of the present invention.

[0089] Obviously, those skilled in the art can make various modifications and variations to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention should also include these modifications and variations.

Claims

1. A shortwave radiation all-factor integrated remote sensing estimation method, characterized in that: include: The albedo of the acquired remote sensing observation data is used as the input of a pre-built simulation library, and the radiation brightness of the blue, green and red bands at the top of the atmosphere or a certain height is obtained based on the atmospheric parameters in the simulation library, including: obtaining the albedo of the blue, green, red bands and the entire shortwave band of different surface objects through satellites; using the albedo as the input of the simulation library, and obtaining the radiation brightness of the blue, green and red bands at the top of the atmosphere or a certain height through the radiation brightness formula based on the atmospheric parameters in the simulation library; the radiation brightness formula is: Where, μ0 is the cosine of the solar zenith angle; μ is the cosine of the observed zenith angle; φ is the relative azimuth; L(μ0,μ,φ) is the radiance of the blue, green or red band at the top of the atmosphere or at a certain height; L0(μ0,μ,φ) is the atmospheric path radiation; ρ is the albedo of the large balloon surface; r1 is the albedo of the blue, green and red bands of different surface objects; ψ is the comprehensive radiation term related to the atmospheric transmittance and the total solar radiation at the top of the atmosphere; Comparing the radiance obtained through the simulation library with the radiance of the acquired remote sensing observation data, and then matching the atmospheric parameters that meet the preset conditions in the simulation library according to the comparison result; According to the albedo and the atmospheric parameters that meet the preset conditions, all elements of shortwave radiation are analyzed, including: according to the albedo and the atmospheric parameters of the screened atmospheric conditions, 8 types of radiation components of shortwave radiation are obtained through the radiation component formula; the radiation component formula is: Where F(μ0) is the downward radiation from the surface or the radiation from the top of the atmosphere; r2 is the albedo of the entire shortwave band of different surface objects; F0(μ0) is the downward radiation from the atmosphere or the upward radiation from the top of the atmosphere without considering the contribution of the surface; γ is the comprehensive radiation term related to the atmospheric transmittance.

2. A shortwave radiation all-factor integrated remote sensing estimation method as claimed in claim 1, characterized in that: Before using the albedo of the acquired remote sensing observations as input to the pre-built simulation library, it also includes: According to preset atmospheric conditions and spectral response functions of different sensors, multiple types of simulated radiation components under different atmospheric conditions and observation angles and atmospheric parameters corresponding to each type of simulated radiation components are obtained through simulation; Building a simulation library based on the multiple types of simulated radiation components and the atmospheric parameters corresponding to each type of simulated radiation component; The atmospheric parameters include: atmospheric transmittance, path radiation, large balloon surface albedo, and comprehensive radiation items related to atmospheric transmittance and total solar radiation at the top of the atmosphere.

3. A shortwave radiation all-factor integrated remote sensing estimation method as claimed in claim 2, characterized in that: The preset atmospheric conditions are conditions retrieved from a radiation transfer simulation condition library constructed based on representative parameters, and the representative parameters include: global atmospheric profiles, cloud parameter products, surface albedo, aerosols, sun angles, and satellite / UAV observation angles.

4. A shortwave radiation all-factor integrated remote sensing estimation method as claimed in claim 1, characterized in that: The radiance obtained by the simulation library is compared with the radiance of the remote sensing observation data, and then the atmospheric parameters that meet the preset conditions are matched in the simulation library according to the comparison result, including: Using the root mean square error as the comparison condition, the radiance obtained through the simulation library is compared with the radiance of the remote sensing observation data; Based on the comparison result, the atmospheric condition with the smallest difference in radiation brightness from the remote sensing observation data is selected in the simulation library.

5. A shortwave radiation all-factor integrated remote sensing estimation method as claimed in any one of claims 1 to 4, characterized in that: The full elements of shortwave radiation include: total downward shortwave radiation from the surface, direct solar radiation, sky diffuse radiation, net surface radiation, photosynthetically active radiation, direct or diffuse, ultraviolet radiation and net shortwave radiation from the top of the atmosphere.

6. A shortwave radiation all-factor integrated remote sensing estimation system, characterized in that: include: The radiation brightness simulation module is used to use the albedo of the acquired remote sensing observation data as the input of the pre-built simulation library, and obtain the radiation brightness of the blue, green and red bands at the top of the atmosphere or a certain height based on the atmospheric parameters in the simulation library, including: obtaining the albedo of the blue, green, red bands and the entire shortwave band of different surface objects through satellites; using the albedo as the input of the simulation library, and obtaining the radiation brightness of the blue, green and red bands at the top of the atmosphere or a certain height through the radiation brightness formula based on the atmospheric parameters in the simulation library; the radiation brightness formula is: Where, μ0 is the cosine of the solar zenith angle; μ is the cosine of the observed zenith angle; φ is the relative azimuth; L(μ0,μ,φ) is the radiance of the blue, green or red band at the top of the atmosphere or at a certain height; L0(μ0,μ,φ) is the atmospheric path radiation; ρ is the albedo of the large balloon surface; r1 is the albedo of the blue, green and red bands of different surface objects; ψ is the comprehensive radiation term related to the atmospheric transmittance and the total solar radiation at the top of the atmosphere; A comparison and matching module, used to compare the radiance obtained through the simulation library with the radiance of the acquired remote sensing observation data, and then match the atmospheric parameters that meet the preset conditions in the simulation library according to the comparison result; The full-factor obtaining module is used to analyze and obtain the full factors of shortwave radiation according to the albedo and the atmospheric parameters that meet the preset conditions, including: according to the albedo and the atmospheric parameters of the screened atmospheric conditions, obtaining 8 types of radiation components of shortwave radiation through the radiation component formula; the radiation component formula is: Where F(μ0) is the downward radiation from the surface or the radiation from the top of the atmosphere; r2 is the albedo of the entire shortwave band of different surface objects; F0(μ0) is the downward radiation from the atmosphere or the upward radiation from the top of the atmosphere without considering the contribution of the surface; γ is the comprehensive radiation term related to the atmospheric transmittance.

Images