Light scattering radiance solving method, device, equipment and medium

Abstract

The invention discloses a light scattering radiance solving method, device, equipment and medium, and the method comprises the steps: obtaining a solar radiation wave band, an aerosol parameter, a cloud parameter, an atmospheric molecule parameter, an earth surface parameter, solar spectrum irradiation, a solar geometric orientation and an observation geometric orientation, and obtaining the light scattering radiance according to the solar radiation wave band, the aerosol parameter, the cloud parameter and the atmospheric molecule parameter; generating scattering coefficients, absorption coefficients and phase functions of the aerosol, the cloud and the atmospheric molecules according to the scattering coefficients, the absorption coefficients and the phase functions, and generating medium extinction coefficients, single scattering albedo and scattering phase functions according to the scattering coefficients, the absorption coefficients and the phase functions; and solving the light scattering radiance by using a Monte Carlo three-dimensional atmospheric transmission model according to the earth surface parameters, the solar spectral irradiation, the solar geometric orientation, the observation geometric orientation, the medium extinction coefficient, the single scattering albedo and the scattering phase function. According to the method, primary and secondary links and element parameters influencing geogas background scattered light are defined, and stray light background quantitative generation is achieved.

Description

Technical Field

[0001] The invention belongs to the technical field of meteorological detection, and in particular relates to a method, device, equipment and medium for solving light scattering radiance. Background Art

[0002] The radiation transfer process in complex scenes is a dynamic process (intrinsic dynamics / geometric dynamics / observational dynamics) facing multiple factors (surface / aerosol / cloud / atmosphere / sun). The current radiation transfer algorithm is overly simplified, the links are not considered enough, and only the static environment is considered, which cannot meet the high-precision requirements in the study of the mechanism and spatiotemporal dynamic evolution of the dynamic optical properties of the surface and atmosphere. Regarding clouds, current studies mostly use one-dimensional radiation transfer models to divide the sky into two situations: clouded and cloudless, which does not meet the situation where clouds are not completely covered in actual scenes; three-dimensional radiation transfer models can approximately simulate the coverage of clouds in actual scenes, considering the mutual influence between clouds at different locations, but lack the impact of the dynamic evolution of clouds on radiation. Regarding the atmosphere, the height of the atmospheric profile in traditional one-dimensional and three-dimensional radiation transfer models is 120km, and there is a lack of consideration of the atmospheric state in the range of 120-200km. Regarding the surface, traditional radiation transfer models consider a single surface type, while actual scenes are mostly mixed surfaces. Regarding aerosols, the choice of aerosol type in radiation transfer is mostly based on the underlying surface (for example, rural aerosol is selected for rural surface), without considering meteorological factors and aerosol transport. This method causes errors in radiation transfer.

[0003] In addition, existing radiation transfer algorithms have various deficiencies, such as the narrow spectral range of the calculation and the use of the plane-parallel assumption in the calculation of atmospheric multiple scattering. Summary of the invention

[0004] The purpose of the present invention is to provide a method, device, equipment and medium for solving light scattering radiance. Based on the Monte Carlo method, a high-precision three-dimensional radiation transmission model of complex spherical atmosphere-surface coupling multiple elements is constructed, which includes a multi-element three-dimensional scene of the surface, aerosols, clouds and atmosphere, and solves the technical problem of insufficient calculation ability of light scattering radiance under the ground-atmosphere background environment in complex scenes in the prior art.

[0005] To solve the above technical problems, the present invention is achieved through the following technical solutions:

[0006] The present invention provides a method for solving light scattering radiance, which comprises:

[0007] Obtaining solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecular parameters, surface parameters, solar spectrum irradiance, solar geometric orientation and observation geometric orientation, wherein the surface parameters include mixed surface parameters;

[0008] Generating scattering coefficients, absorption coefficients and phase functions of aerosols, clouds and atmospheric molecules according to the solar radiation bands, aerosol parameters, cloud parameters and atmospheric molecule parameters;

[0009] According to the scattering coefficient, absorption coefficient and phase function, the medium extinction coefficient, single scattering albedo and scattering phase function of the interaction between solar radiation and the aerosol, cloud and atmospheric molecules when the solar radiation propagates in the atmosphere are generated;

[0010] According to the surface parameters, solar spectral irradiance, solar geometric orientation, observation geometric orientation, medium extinction coefficient, single scattering albedo and scattering phase function, a Monte Carlo three-dimensional atmospheric transmission model is used to solve the light scattering radiance.

[0011] In one embodiment of the present invention, the obtaining of solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecule parameters, surface parameters, solar spectrum irradiance, solar geometric orientation and observation geometric orientation includes:

[0012] Acquire solar radiation bands in a simulated space region, where the solar radiation bands include visible light bands and invisible light bands;

[0013] Acquiring aerosol parameters of the simulated space region, wherein the aerosol parameters include aerosol type and aerosol optical thickness;

[0014] Obtaining cloud parameters of the simulated space region, wherein the cloud parameters include cloud type and cloud optical thickness;

[0015] Acquiring atmospheric molecular parameters of a simulated space region, wherein the atmospheric molecular parameters include visibility, molecular absorption spectrum, and atmospheric profile;

[0016] Acquire surface parameters of the simulated space area, wherein the surface parameters include surface type and surface reflection characteristics;

[0017] Obtain the solar spectral irradiance, solar geometric orientation and observation geometric orientation of the simulated space area.

[0018] In one embodiment of the present invention, the step of obtaining aerosol parameters of the simulated space region, wherein the aerosol parameters include aerosol type and aerosol optical thickness, includes:

[0019] Obtain meteorological and geographical elements of the simulated space area;

[0020] Determining an aerosol model according to the meteorological elements and the geographical elements, and using the aerosol model to obtain the aerosol type;

[0021] A graph neural network model that takes into account temporal and spatial heterogeneity is used to obtain the aerosol optical depth.

[0022] In one embodiment of the present invention, the step of obtaining atmospheric molecular parameters of a simulated space region, wherein the atmospheric molecular parameters include visibility, molecular absorption spectrum and atmospheric profile, comprises:

[0023] Obtain atmospheric visibility and molecular absorption spectrum parameters of simulated space area;

[0024] According to the spherical distribution characteristics of the atmospheric profile, the atmospheric profile in the altitude range of 0 to 200km is obtained.

[0025] In one embodiment of the present invention, the step of obtaining the atmospheric profile within the altitude range of 0 to 200 km according to the spherical distribution characteristics of the atmospheric profile comprises:

[0026] Taking the center of the earth as the center, the spherical atmosphere in the range of 0 to 200 km altitude is layered;

[0027] Based on the temperature, humidity data and standard atmospheric profiles fused from multiple sources, the atmospheric profiles within the altitude range of 0 to 100 km are obtained;

[0028] Based on the NRLMSISIE-00 model, the atmospheric profiles in the altitude range of 100 to 200 km are obtained.

[0029] In one embodiment of the present invention, the acquiring of the surface parameters of the simulated space region, wherein the surface parameters include the surface type and the surface reflection characteristics, comprises:

[0030] Obtain the type of surface features in the simulated space area to determine the surface type;

[0031] The bidirectional reflectance distribution function model is used to characterize the reflection characteristics of the ground object types using different model parameters in different surface areas, so as to obtain the surface reflection characteristics of the surface in the simulated space area.

[0032] In one embodiment of the present invention, generating the scattering coefficient, absorption coefficient and phase function of aerosol, cloud and atmospheric molecules according to the solar radiation band, aerosol parameters, cloud parameters and atmospheric molecule parameters includes:

[0033] According to the solar radiation band, aerosol parameters, and cloud parameters, a Mie scattering calculation model is used to obtain the scattering coefficient, absorption coefficient, and phase function of the aerosol and cloud;

[0034] According to the solar radiation band and atmospheric molecular parameters, a molecular absorption calculation model is used to obtain the absorption coefficient of the atmospheric molecules;

[0035] According to the solar radiation band and atmospheric molecule parameters, the Rayleigh scattering theory is used to obtain the scattering coefficient and phase function of the atmospheric molecules.

[0036] Based on the same inventive concept, another embodiment of the present invention further provides a light scattering radiance solution device, the device comprising:

[0037] A data acquisition module, used to acquire solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecular parameters, surface parameters, solar spectrum irradiance, solar geometric orientation and observation geometric orientation, wherein the surface parameters include mixed surface parameters;

[0038] A data processing module, for generating scattering coefficients, absorption coefficients and phase functions of aerosols, clouds and atmospheric molecules according to the solar radiation bands, aerosol parameters, cloud parameters and atmospheric molecule parameters, and for generating medium extinction coefficients, single scattering albedo and scattering phase functions of the interaction between solar radiation and the aerosols, clouds and atmospheric molecules when the solar radiation propagates in the atmosphere according to the scattering coefficients, absorption coefficients and phase functions;

[0039] The radiance solution module is used to solve the light scattering radiance by using the Monte Carlo three-dimensional atmospheric transmission model according to the surface parameters, solar spectral irradiance, solar geometric orientation, observation geometric orientation, medium extinction coefficient, single scattering albedo and scattering phase function.

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

[0041] one or more processors;

[0042] A storage device is used to store one or more programs, and when the one or more programs are executed by the one or more processors, the electronic device implements any of the above-mentioned light scattering radiance solution methods.

[0043] 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 any of the above-mentioned methods for solving the light scattering radiance.

[0044] As described above, a method for solving light scattering radiance provided by the present invention has the following beneficial effects: the method for solving light scattering radiance obtains solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecular parameters, surface parameters, solar spectral irradiance, solar geometric orientation and observation geometric orientation, wherein the surface parameters include mixed surface parameters; according to the solar radiation bands, aerosol parameters, cloud parameters and atmospheric molecular parameters, the scattering coefficients, absorption coefficients and phase functions of aerosols, clouds and atmospheric molecules are generated; according to the scattering coefficients, absorption coefficients and phase functions, the medium extinction coefficients, single scattering albedo and scattering phase functions of the interaction between solar radiation and the aerosols, clouds and atmospheric molecules when propagating in the atmosphere are generated; according to the surface parameters, solar spectral irradiance, solar geometric orientation, observation geometric orientation, medium extinction coefficients, single scattering albedo and scattering phase functions, the Monte Carlo three-dimensional atmospheric transmission model is used to solve the light scattering radiance. The method for solving the light scattering radiance takes into account the spatial position relationship between the sun, the earth and the space payload, as well as the physical properties of the earth's surface, clouds and the atmosphere and their changes with space and time, and obtains the specific spatial distribution and intensity of the ground-air scattering and the law of change over time, which changes the current situation of lack of relevant research data support in the design process of space payloads, adopting simplified processing methods, using average distribution for spatial distribution, and using rough empirical values ​​or statistical values ​​for intensity, and provides theoretical support for the analysis and design of stray light suppression and elimination of space payloads, and supports the precise design, planning and scheduling of space missions. 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

[0045] 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.

[0046] Figure 1 A schematic flow chart of a method for solving light scattering radiance provided in accordance with an exemplary embodiment of the present application.

[0047] Figure 2 A schematic diagram of the construction process of a high-precision three-dimensional radiation transfer model of a complex spherical atmosphere-surface coupling multi-factor provided for an exemplary embodiment of the present application.

[0048] Figure 3 A schematic diagram of an aerosol type selection method provided for an exemplary embodiment of the present application.

[0049] Figure 4A schematic diagram of cirrus cloud extinction provided for an exemplary embodiment of the present application.

[0050] Figure 5 A schematic diagram of a scattered cirrus cloud provided by an exemplary embodiment of the present application.

[0051] Figure 6 A schematic diagram of a flow chart of constructing an atmospheric profile within an altitude range of 0 to 200 km provided for an exemplary embodiment of the present application.

[0052] Figure 7 A schematic diagram of mixed land feature types provided for an exemplary embodiment of the present application.

[0053] Figure 8 A flow chart of a Monte Carlo algorithm provided for an exemplary embodiment of the present application.

[0054] Fig. 9 A schematic diagram of a model software interface provided for an exemplary embodiment of the present application.

[0055] Fig.10 An exemplary embodiment of the present application provides a radiance calculation result for earth observation at a wavelength of 0.4 to 1.1 μm.

[0056] Fig.11 The radiance calculation result of 0.4 to 1.1 μm wavelength limb observation is provided for an exemplary embodiment of the present application.

[0057] Fig.12 A schematic structural diagram of a light scattering radiance solving device provided for another exemplary embodiment of the present application.

[0058] Fig.13 A schematic structural diagram of an electronic device provided for another exemplary embodiment of the present application. DETAILED DESCRIPTION

[0059] 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.

[0060] 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.

[0061] 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.

[0062] As a huge source of energy, the sun continuously emits radiation of various wavelengths to the earth. When these radiations pass through the atmosphere, they interact with various components in the atmosphere. Some components in the atmosphere have the characteristic of selectively absorbing radiation energy of a certain wavelength. For example, water vapor has strong absorption in the infrared region, while ozone strongly absorbs ultraviolet rays. These absorption effects will cause the radiation intensity to weaken.

[0063] Scattering is the change in propagation direction of radiation when it propagates in the atmosphere due to encountering atmospheric molecules, aerosol particles, cloud particles, etc. Scattering causes radiation to propagate in all directions, thereby changing the intensity and direction of radiation received by the surface. The reflectivity of the surface to solar radiation depends on the nature and state of the surface (such as land, water surface, snow surface, etc.). Part of the solar radiation is reflected back into the atmosphere by the surface, forming surface reflected radiation. The surface also absorbs part of the solar radiation and converts it into heat energy.

[0064] In order to solve the technical problem of insufficient calculation ability of light scattering radiance under ground-atmosphere background environment in complex scenes in the prior art, the present invention provides a method for solving light scattering radiance. Based on the Monte Carlo method, a high-precision three-dimensional radiation transmission model of complex spherical atmosphere-ground surface coupling multi-factors is constructed, including a multi-factor three-dimensional scene of ground surface, aerosol, cloud and atmospheric molecules, which clarifies the contribution of various factors such as ground surface, aerosol, cloud and atmospheric molecules to ground-atmosphere background scattered light, clarifies the primary and secondary links and factor parameters that affect ground-atmosphere background scattered light, and realizes the quantitative generation of stray light background. Please refer to Figure 1 As shown, the method for solving the light scattering radiance includes the following steps:

[0065] S100: Acquire solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecule parameters, surface parameters, solar spectrum irradiance, solar geometric orientation, and observation geometric orientation, wherein the surface parameters include mixed surface parameters;

[0066] S200: generating scattering coefficients, absorption coefficients and phase functions of aerosols, clouds and atmospheric molecules according to the solar radiation band, aerosol parameters, cloud parameters and atmospheric molecule parameters;

[0067] S300: generating a medium extinction coefficient, a single scattering albedo and a scattering phase function of the interaction between solar radiation and the aerosol, cloud and atmospheric molecules when the solar radiation propagates in the atmosphere according to the scattering coefficient, absorption coefficient and phase function;

[0068] S400: according to the surface parameters, solar spectral irradiance, solar geometric orientation, observation geometric orientation, medium extinction coefficient, single scattering albedo and scattering phase function, using a Monte Carlo three-dimensional atmospheric transmission model to solve the light scattering radiance.

[0069] See also Figure 1 and Figure 2 As shown, the steps of the above-mentioned method for solving the light scattering radiance will be discussed in detail below.

[0070] First, step S100 is executed, namely obtaining solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecule parameters, surface parameters, solar spectrum irradiance, solar geometric orientation and observation geometric orientation, wherein the surface parameters include mixed surface parameters.

[0071] In an exemplary embodiment of the present application, in step S100, obtaining solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecule parameters, surface parameters, solar spectrum irradiance, solar geometric orientation, and observation geometric orientation further includes the following steps:

[0072] S110: Acquire a solar radiation band of a simulated space region, where the solar radiation band includes a visible light band and an invisible light band;

[0073] S120: Acquire aerosol parameters of the simulated space area, where the aerosol parameters include aerosol type and aerosol optical thickness;

[0074] S130: Acquire cloud parameters of the simulated space area, where the cloud parameters include cloud type and cloud optical thickness;

[0075] S140: Acquire atmospheric molecular parameters of the simulated space region, where the atmospheric molecular parameters include visibility, molecular absorption spectrum, and atmospheric profile;

[0076] S150: Acquire surface parameters of the simulated space area, where the surface parameters include surface type and surface reflection characteristics;

[0077] S160: Obtain the solar spectrum irradiance, solar geometric orientation and observation geometric orientation of the simulated space area.

[0078] It should be noted that, in this embodiment, the simulated space area is divided into three-dimensional atmospheric grids according to latitude, longitude and altitude, and each atmospheric grid is numbered respectively. It is worth noting that the optical properties in each grid are fixed and unchanged, and the atmospheric grid mainly includes the extinction coefficient, single scattering albedo and scattering phase function of various atmospheric components.

[0079] In an exemplary embodiment of the present application, in step S120, aerosol parameters of the simulated space region are obtained, wherein the aerosol parameters include aerosol type and aerosol optical thickness, and the following steps are further included:

[0080] S121: Acquire meteorological elements and geographical elements of the simulated space area;

[0081] S122: determining an aerosol model according to the meteorological elements and the geographical elements, and using the aerosol model to obtain the aerosol type;

[0082] S123: Utilize a graph neural network model that takes into account temporal and spatial heterogeneity to obtain the aerosol optical thickness.

[0083] It should be noted that the meteorological elements include but are not limited to wind speed and wind direction, and the geographical elements include but are not limited to the location of the aerosol source. In addition, in this embodiment, in step S122, a suitable aerosol model is selected using random forest according to the meteorological elements and the geographical elements. Of course, in other embodiments, other methods can also be used to determine the aerosol model.

[0084] For details, please refer to Figure 3 As shown, during the simulation process, aerosol parameters are described by aerosol type and aerosol optical thickness. In order to improve the accuracy of aerosol optical properties and the impact of aerosol dynamic changes on radiation transmission, it is necessary to select a suitable aerosol model. In this embodiment, an aerosol model evaluation and trajectory analysis method is adopted, with meteorological elements and geographical elements as input, and the corresponding aerosol model is output. Compared with the traditional method of selecting aerosol models using measured data or special statistical models, the complexity of the method is greatly reduced, and the accuracy is improved compared to selecting aerosol models based on the underlying surface. In addition, the appropriate aerosol type can be dynamically selected according to the changes in meteorological elements. For aerosol optical thickness, a graph neural network (GSTH-GNN) model that takes into account spatiotemporal heterogeneity is proposed to predict aerosol optical thickness and improve the accuracy of the optical thickness input parameters.

[0085] In an exemplary embodiment of the present application, in step S130, cloud parameters are obtained, and the cloud parameters include cloud type and cloud optical thickness. Figure 4 and Figure 5 As shown, Figure 4and Figure 5 The three-dimensional optical properties of clouds are shown. During the simulation process, the three-dimensional distribution characteristics of clouds and their changes over time are fully considered. The absorption coefficient, extinction coefficient, and single scattering albedo are used to characterize the three-dimensional properties of clouds, and the influence of clouds on scattered light is dynamically considered.

[0086] In an exemplary embodiment of the present application, in step S140, the atmospheric molecular parameters of the simulated space region are obtained, wherein the atmospheric molecular parameters include visibility, molecular absorption spectrum and atmospheric profile, and the following steps are further included:

[0087] S141: Obtain atmospheric visibility and molecular absorption spectrum parameters of the simulated space area;

[0088] S142: Based on the spherical distribution characteristics of the atmospheric profile, the atmospheric profile within the altitude range of 0 to 200 km is obtained.

[0089] In an exemplary embodiment of the present application, in step S142, obtaining the atmospheric profile within the altitude range of 0 to 200 km according to the spherical distribution characteristics of the atmospheric profile further includes the following steps:

[0090] S1421: Stratify the spherical atmosphere within the range of 0 to 200 km altitude with the center of the earth as the center of the circle;

[0091] S1422: Obtain the atmospheric profile within the altitude range of 0 to 100 km based on the temperature, humidity data and standard atmospheric profile fused from multiple sources;

[0092] S1423: Based on the NRLMSISIE-00 model, obtain the atmospheric profile in the altitude range of 100 to 200 km.

[0093] For details, please refer to Figure 6As shown, in order to improve simulation accuracy, the influence of spherical atmosphere is considered in the multiple scattering process, the atmospheric profile data feature is that the atmospheric profile in the range of 0 to 100km altitude adopts the temperature and humidity data based on multi-source data fusion and the standard atmospheric profile to obtain, and the atmospheric profile in the range of 100 to 200km altitude adopts NRLMSISIE-00 model data. In addition, for the construction of spherical atmosphere, it is considered that various profiles have non-uniformity in the vertical direction, and it is necessary to carry out layered processing to spherical atmosphere. In the present embodiment, with the center of the earth as the center of circle, 0km altitude (i.e., surface altitude) is the bottom of the atmosphere, and 200km altitude is the top of the atmosphere, and the atmosphere is divided into a series of symmetrical concentric spherical layers by vertical height in the range of 0 to 200km altitude, the pressure of the lower atmosphere is high, the density is large, and the profile changes more violently, and the layering should be more meticulous, therefore in the range of 0 to 50km altitude, one layer is divided by vertical height every 1 to 2km. The pressure of the middle and upper atmosphere is low, the density is small, and the profile changes relatively smoothly, so the stratification can be appropriately rough. Therefore, within the altitude range of 50 to 100 km, one layer is divided every 10 km of vertical height, and within the altitude range of 100 to 200 km, the whole is divided into one layer. It should be noted that in this embodiment, the spherical atmosphere within the altitude range of 0 to 200 km is divided into 35 layers. Of course, in other embodiments, the specific number of layers of the spherical atmosphere can be adaptively set according to actual application requirements.

[0094] In an exemplary embodiment of the present application, in step S150, obtaining the surface parameters of the simulated space area, wherein the surface parameters include the surface type and the surface reflection characteristics, further comprises the following steps:

[0095] S151: Obtain the type of ground objects on the surface of the simulated space area to determine the type of the surface;

[0096] S152: Using a bidirectional reflectance distribution function model, different model parameters are used in different surface areas to characterize the reflection characteristics of the ground object type, so as to obtain the surface reflection characteristics of the surface in the simulated space area.

[0097] For details, please refer to Figure 7As shown, remote sensing images and geographic information system data are used to identify the type of surface objects in the simulated space area to determine the surface type, which includes a single surface or a mixed surface, and the surface types include lakes, seas, deserts, mountains, ice and snow, vegetation, bare soil, etc. According to the surface parameter database, standard reflection spectrum curves of different surface types are obtained, and the parameters of the reflection spectrum curves are adjusted according to the actual conditions of different surfaces. According to the research purpose and surface characteristics, a suitable surface radiation model is selected. In this embodiment, the surface radiation model adopts a bidirectional reflection distribution function model. The surface in the simulated space area is divided into multiple sub-areas, each sub-area corresponds to one or more surface types, and the reflection characteristics of each sub-area are simulated using a bidirectional reflection distribution function model, that is, different model parameters are used in different sub-areas to characterize the reflection characteristics of the object type, and the model parameters are not limited to reflectivity, transmittance and absorptivity. The simulation results of each sub-area are superimposed in proportion to obtain the reflection characteristics of the mixed surface.

[0098] Next, step S200 is executed to generate scattering coefficients, absorption coefficients and phase functions of aerosols, clouds and atmospheric molecules according to the solar radiation band, aerosol parameters, cloud parameters and atmospheric molecule parameters.

[0099] In an exemplary embodiment of the present application, in step S200, generating the scattering coefficient, absorption coefficient and phase function of aerosol, cloud and atmospheric molecules according to the solar radiation band, aerosol parameters, cloud parameters and atmospheric molecule parameters further includes the following steps:

[0100] S210: according to the solar radiation band, aerosol parameters, and cloud parameters, using a Mie scattering calculation model to obtain the scattering coefficient, absorption coefficient, and phase function of the aerosol and cloud;

[0101] S220: According to the solar radiation band and the atmospheric molecule parameters, using a molecular absorption calculation model to obtain the absorption coefficient of the atmospheric molecules;

[0102] S230: According to the solar radiation band and the atmospheric molecule parameters, the Rayleigh scattering theory is used to obtain the scattering coefficient and phase function of the atmospheric molecules.

[0103] Next, step S300 is performed to generate the medium extinction coefficient, single scattering albedo and scattering phase function of the interaction between the solar radiation and the aerosol, cloud and atmospheric molecules when the solar radiation propagates in the atmosphere according to the scattering coefficient, absorption coefficient and phase function.

[0104] Finally, step S400 is performed to solve the light scattering radiance using a Monte Carlo three-dimensional atmospheric transmission model according to the surface parameters, solar spectral irradiance, solar geometric orientation, observation geometric orientation, medium extinction coefficient, single scattering albedo and scattering phase function.

[0105] The solution process of light scattering radiance will be described in detail below.

[0106] See also Figure 8 As shown, a Monte Carlo Markov random chain is constructed based on the integrated radiative transfer equation and the radiance contribution is locally estimated.

[0107] Initialize the optical properties of the three-dimensional spherical atmosphere, divide the simulated space area into three-dimensional grids according to longitude and latitude (0.5°×0.5°) and altitude, and number them. The optical properties in each grid are fixed, mainly including the extinction coefficient, single scattering albedo and scattering phase function of various atmospheric components in the atmospheric grid.

[0108] Initialize the three-dimensional cloud atmospheric optical properties, determine the cloud position information in three-dimensional space, determine the three-dimensional grid position covered by the cloud, and use the extinction coefficient, single scattering albedo and scattering phase function in the grid to characterize the optical properties of the cloud.

[0109] Initialize the optical properties of the two-dimensional surface, divide the surface in the simulated space into grids according to longitude and latitude (0.5°×0.5°), select a suitable bidirectional reflectance distribution function model, and use different model parameters in each grid to characterize the reflection characteristics in the grid.

[0110] Initialize the sample state according to the source function or the detector response function.

[0111] The random walking process of samples in the atmosphere is simulated, that is, the free path sampling process.

[0112] A newly proposed improved analytical expression of single scattering phase function is proposed to simulate the collision event of the sample, namely the atmospheric scattering event. This expression has advantages in calculating the backscattering part of spherical particles and non-spherical particles.

[0113] Update sample weights, sample collision types, sample collision directions, and locally estimate radiance contribution. The local estimation of radiance contribution is to achieve the purpose of statistical calculation of radiance by Monte Carlo method, and the radiance estimation formula is directly used in the simulation process.

[0114] Repeat the last three steps of the simulation until the sample is absorbed by the atmosphere or escapes the atmosphere, ending the Markov random chain simulation process. Count the number and weight of the received samples to obtain the light scattering radiance of the study area.

[0115] In summary, the present invention provides a method for solving light scattering radiance, which uses the Monte Carlo method to simulate the three-dimensional radiation transmission process of the dynamic process of multiple factors coupled by the atmosphere, clouds, surface and aerosols. The basic principle is to simulate the absorption and scattering process of multiple factors (spherical atmosphere, clouds, surface, aerosols) on the propagation path of light from the detector to the light source, and use the absorption coefficient, extinction coefficient, and single scattering albedo to describe the absorption and scattering process of the atmosphere, clouds and aerosols. The bidirectional reflectance distribution function model (BRDF model) is used to describe the reflection characteristics of the surface. The radiation transmission problem is solved using statistical approximation methods. When light propagates in the atmosphere, it will have complex interactions with gas molecules, aerosol particles, cloud particles and the ground. This interaction can be described by parameters such as scattering phase function, single scattering albedo, attenuation coefficient, and surface reflectivity. During the simulation process, the spherical distribution characteristics of the atmospheric profile are considered in the multiple scattering process. Please refer to 9 to Fig.11 As shown, Fig. 9 shows the model software interface, Fig.10 and Fig.11 The radiance calculation results of the light scattering radiance solution method in ground observation and edge observation according to the present invention are shown respectively. The light scattering radiance solution method fully considers the spatial position relationship between the sun, the earth and the space payload, as well as the physical properties of the surface, clouds and atmosphere and their changes with space and time, and obtains the specific spatial distribution and intensity of the ground-atmosphere scattering and the law of change with time, which changes the current situation of lack of relevant research data support in the design process of space payloads, using simplified processing methods, using average distribution for spatial distribution, and using rough empirical values ​​or statistical values ​​for intensity, providing theoretical support for the analysis and design of stray light suppression and elimination of space payloads, and supporting the precise design, planning and scheduling of space missions.

[0116] Based on the same inventive concept, please refer to Fig.12 As shown, another embodiment of the present invention further provides a light scattering radiance solution device 11, the device comprising:

[0117] The data acquisition module 111 is used to acquire solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecular parameters, surface parameters, solar spectrum irradiance, solar geometric orientation and observation geometric orientation, wherein the surface parameters include mixed surface parameters;

[0118] The data processing module 112 is used to generate the scattering coefficient, absorption coefficient and phase function of the aerosol, cloud and atmospheric molecules according to the solar radiation band, aerosol parameters, cloud parameters and atmospheric molecule parameters, and to generate the medium extinction coefficient, single scattering albedo and scattering phase function of the interaction between the solar radiation and the aerosol, cloud and atmospheric molecules when the solar radiation propagates in the atmosphere according to the scattering coefficient, absorption coefficient and phase function;

[0119] The radiance solution module 113 is used to solve the light scattering radiance by using the Monte Carlo three-dimensional atmospheric transmission model according to the surface parameters, solar spectral irradiance, solar geometric orientation, observation geometric orientation, medium extinction coefficient, single scattering albedo and scattering phase function.

[0120] It should be noted that the light scattering radiance solution device 11 includes the light scattering radiance solution method described in any of the above embodiments. Since the light scattering radiance solution device 11 provided in this embodiment belongs to the same inventive concept as the light scattering radiance solution method provided in any of the above embodiments, it has at least the same beneficial effects as the light scattering radiance solution method, and will not be described one by one here.

[0121] Based on the same inventive concept, please refer to Fig.13 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 a light scattering radiance solution program.

[0122] 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 solving the light scattering radiance, etc., but also can be used to temporarily store data that has been output or is to be output.

[0123] 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 light scattering radiance solution 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.

[0124] 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 light scattering radiance solution method.

[0125] 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 a radiance solution module 113.

[0126] 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 light scattering radiance solution method described in each embodiment of the present application.

[0127] 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 solving light scattering radiance, characterized in that: include: Obtaining solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecular parameters, surface parameters, solar spectrum irradiance, solar geometric orientation and observation geometric orientation, wherein the surface parameters include mixed surface parameters; Generating scattering coefficients, absorption coefficients and phase functions of aerosols, clouds and atmospheric molecules according to the solar radiation bands, aerosol parameters, cloud parameters and atmospheric molecule parameters; According to the scattering coefficient, absorption coefficient and phase function, the medium extinction coefficient, single scattering albedo and scattering phase function of the interaction between solar radiation and the aerosol, cloud and atmospheric molecules when the solar radiation propagates in the atmosphere are generated; According to the surface parameters, solar spectral irradiance, solar geometric orientation, observation geometric orientation, medium extinction coefficient, single scattering albedo and scattering phase function, a Monte Carlo three-dimensional atmospheric transmission model is used to solve the light scattering radiance.

2. The method for calculating light scattering radiance according to claim 1, characterized in that: The obtaining of solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecule parameters, surface parameters, solar spectrum irradiance, solar geometric orientation and observation geometric orientation includes: Acquire solar radiation bands in a simulated space region, where the solar radiation bands include visible light bands and invisible light bands; Acquiring aerosol parameters of the simulated space region, wherein the aerosol parameters include aerosol type and aerosol optical thickness; Obtaining cloud parameters of the simulated space region, wherein the cloud parameters include cloud type and cloud optical thickness; Acquiring atmospheric molecular parameters of a simulated space region, wherein the atmospheric molecular parameters include visibility, molecular absorption spectrum, and atmospheric profile; Acquire surface parameters of the simulated space area, wherein the surface parameters include surface type and surface reflection characteristics; Obtain the solar spectral irradiance, solar geometric orientation and observation geometric orientation of the simulated space area.

3. The method for calculating light scattering radiance according to claim 2, characterized in that: The aerosol parameters of the simulated space region are obtained, wherein the aerosol parameters include aerosol type and aerosol optical thickness, including: Obtain meteorological and geographical elements of the simulated space area; Determining an aerosol model according to the meteorological elements and the geographical elements, and using the aerosol model to obtain the aerosol type; A graph neural network model that takes into account temporal and spatial heterogeneity is used to obtain the aerosol optical depth.

4. The method for solving the light scattering radiance according to claim 2, characterized in that: The obtaining of atmospheric molecular parameters of the simulated space region, wherein the atmospheric molecular parameters include visibility, molecular absorption spectrum and atmospheric profile, comprises: Obtain atmospheric visibility and molecular absorption spectrum parameters of simulated space area; According to the spherical distribution characteristics of the atmospheric profile, the atmospheric profile in the altitude range of 0 to 200km is obtained.

5. The method for calculating light scattering radiance according to claim 4, characterized in that: The method of obtaining the atmospheric profile within the range of 0 to 200 km altitude according to the spherical distribution characteristics of the atmospheric profile comprises: Taking the center of the earth as the center, the spherical atmosphere in the range of 0 to 200 km altitude is layered; Based on the temperature, humidity data and standard atmospheric profiles fused from multiple sources, the atmospheric profiles within the altitude range of 0 to 100 km are obtained; Based on the NRLMSISIE-00 model, the atmospheric profiles in the altitude range of 100 to 200 km are obtained.

6. The method for calculating light scattering radiance according to claim 2, characterized in that: The acquiring of the surface parameters of the simulated space region, wherein the surface parameters include the surface type and the surface reflection characteristics, comprises: Obtain the type of surface features in the simulated space area to determine the surface type; The bidirectional reflectance distribution function model is used to characterize the reflection characteristics of the ground object types using different model parameters in different surface areas, so as to obtain the surface reflection characteristics of the surface in the simulated space area.

7. The method for calculating light scattering radiance according to claim 1, characterized in that: The generating of the scattering coefficient, absorption coefficient and phase function of aerosol, cloud and atmospheric molecules according to the solar radiation band, aerosol parameters, cloud parameters and atmospheric molecule parameters includes: According to the solar radiation band, aerosol parameters, and cloud parameters, a Mie scattering calculation model is used to obtain the scattering coefficient, absorption coefficient, and phase function of the aerosol and cloud; According to the solar radiation band and atmospheric molecular parameters, a molecular absorption calculation model is used to obtain the absorption coefficient of the atmospheric molecules; According to the solar radiation band and atmospheric molecule parameters, the Rayleigh scattering theory is used to obtain the scattering coefficient and phase function of the atmospheric molecules.

8. A light scattering radiance solving device, characterized in that: The device comprises: A data acquisition module, used to acquire solar radiation bands, aerosol parameters, cloud parameters, atmospheric molecular parameters, surface parameters, solar spectrum irradiance, solar geometric orientation and observation geometric orientation, wherein the surface parameters include mixed surface parameters; A data processing module, for generating scattering coefficients, absorption coefficients and phase functions of aerosols, clouds and atmospheric molecules according to the solar radiation bands, aerosol parameters, cloud parameters and atmospheric molecule parameters, and for generating medium extinction coefficients, single scattering albedo and scattering phase functions of the interaction between solar radiation and the aerosols, clouds and atmospheric molecules when the solar radiation propagates in the atmosphere according to the scattering coefficients, absorption coefficients and phase functions; The radiance solution module is used to solve the light scattering radiance by using the Monte Carlo three-dimensional atmospheric transmission model according to the surface parameters, solar spectral irradiance, solar geometric orientation, observation geometric orientation, medium extinction coefficient, single scattering albedo and scattering phase function.

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 light scattering radiance solution 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 method for solving the light scattering radiance as claimed in any one of claims 1 to 7.

Images