Geogas light scattering physical property parameterization characterization method and device

Abstract

The invention provides a geogas light scattering physical property parameterization characterization method and a geogas light scattering physical property parameterization characterization device. The method comprises the following steps: dividing a detection scene of a detector according to relative positions of the sun, the earth and the detector; the detection scene comprises a limb detection scene and a non-limb detection scene; respectively acquiring geometric feature parameter values and environmental feature parameter values in a plurality of limb detection scenes and a plurality of non-limb detection scenes; traversing array sets of the geometric feature parameter values and the environmental feature parameter values in the two detection scenes respectively, and performing coupled radiation transmission simulation on each array set to obtain a first mass database and a second mass database respectively; and according to the first mass database and the second mass database, respectively obtaining equivalent scattering distribution function models in a limb detection scene and a non-limb detection scene, and carrying out parameterized characterization by using geogas light scattering physical properties. The method has high light scattering calculation precision while meeting the actual demand of space-time dynamic evolution of a geogas system.

Description

Technical Field

[0001] The embodiments of the present invention relate to the technical field of ground-air background light scattering radiation, and in particular to a method and device for parameterizing the physical properties of ground-air light scattering. Background Art

[0002] The coupled radiation transmission of "atmosphere-surface" is used to calculate the external ground-air scattering / radiation of the "Earth-atmosphere" system. The Earth-atmosphere system refers to a space system with the Earth's surface as the lower limit and the upper limit of the atmosphere as the upper limit. Ground-air scattering / radiation is a unique perspective of space detection, generally referring to various types of scattering / radiation released by the Earth-atmosphere system into space, also called ground-air light. In the visible and near-infrared bands, ground-air light is mainly formed by the scattered sunlight of the Earth and the atmosphere.

[0003] There are many factors that affect the scattering characteristics of ground-air light, which are strongly correlated with the reflection characteristics of ground objects and clouds, spatial distribution, time variation, direction of illumination and detection, atmospheric scattering and transmission state, etc., and there is also strong coupling between the factors, such as Figure 1 However, the existing radiation transfer calculation methods only consider static scenes and cannot meet the actual requirements of high precision under the dynamic evolution of space and time.

[0004] Therefore, there is an urgent need to provide a method and device for parameterizing the physical properties of ground-air light scattering. Summary of the invention

[0005] In order to solve the problem that traditional radiation transmission methods cannot meet the actual requirements of high precision of spatiotemporal dynamic evolution conditions, an embodiment of the present invention provides a method and device for parameterizing the physical properties of ground-air light scattering.

[0006] In a first aspect, an embodiment of the present invention provides a method for parameterizing the physical properties of ground-air light scattering, the method comprising:

[0007] According to the relative positions of the sun, the earth and the detector, the detection scenes of the detector are divided; the detection scenes include edge detection scenes and non-edge detection scenes;

[0008] Obtaining geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in several edge detection scenarios and non-edge detection scenarios respectively;

[0009] traversing array sets of geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in two detection scenes respectively, and performing coupled radiation transmission simulation on each array set to obtain a first massive database and a second massive database respectively;

[0010] Based on the first massive database and the second massive database, equivalent scattering distribution function models in edge detection scenarios and non-edge detection scenarios are obtained respectively, and parameterized representation is performed using the ground-air light scattering properties.

[0011] In a second aspect, an embodiment of the present invention further provides a device for parameterizing the physical properties of ground-air light scattering, the device comprising:

[0012] A division unit, used for dividing the detection scene of the detector according to the relative positions of the sun, the earth and the detector; the detection scene includes an edge detection scene and a non-edge detection scene;

[0013] An acquisition unit, used to respectively acquire geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in a plurality of edge detection scenarios and non-edge detection scenarios;

[0014] A simulation unit, used for respectively traversing array sets of geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in two detection scenes, and performing coupled radiation transmission simulation on each array set to obtain a first massive database and a second massive database respectively;

[0015] The characterization unit is used to obtain the equivalent scattering distribution function model in the edge detection scenario and the non-edge detection scenario respectively according to the first massive database and the second massive database, and perform parameterized characterization with the ground-air light scattering properties.

[0016] In a third aspect, an embodiment of the present invention further provides a computing device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the method described in any embodiment of this specification is implemented.

[0017] In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, enables the computer to execute the method described in any embodiment of this specification.

[0018] On the other hand, an embodiment of the present application further provides a computer program product, which includes a computer program. A processor of a computer device reads the computer program from a computer-readable storage medium, and the processor executes the computer program, so that the computer device executes any of the methods described in the above embodiments.

[0019] The embodiment of the present invention provides a method and device for parameterizing the physical properties of ground-atmosphere light scattering. First, according to the relative positions of the sun, the earth and the detector, the detection scene of the detector is divided into a limbal detection scene or a non-limbal detection scene. Then, a number of geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in each detection scene are arranged and combined to form an array set, and each array set is simulated to obtain a first massive database and a second massive database, respectively. Finally, equivalent scattering distribution function models in the two detection scenes are constructed based on the two massive databases. Based on the equivalent scattering distribution function model in the two detection scenes, the scattered radiance distribution in the hemispherical space under any ground-atmosphere background can be output according to the environmental characteristic parameter values ​​of the surface, cloud clusters, aerosols, etc. with spatial longitude and latitude distribution, combined with the geometric characteristic parameter values; in this way, while meeting the actual needs of the spatiotemporal dynamic evolution of the ground-atmosphere system, it has a high light scattering calculation accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0021] Figure 1 It is a schematic diagram of complex coupled radiation transmission of a spatiotemporal dynamic earth-atmosphere system provided by an embodiment of the present invention;

[0022] Figure 2 It is a flow chart of a method for parameterizing the physical properties of ground-air light scattering provided by one embodiment of the present invention;

[0023] Figure 3 It is a schematic diagram of parameterized characterization of light scattering in a non-edge detection scene provided by an embodiment of the present invention;

[0024] Figure 4 is a schematic diagram of parameterized characterization of light scattering in a limb detection scene provided by an embodiment of the present invention;

[0025] Figure 5 It is a schematic diagram of a method for parameterizing the equivalent air-light scattering properties in a detection scenario provided by an embodiment of the present invention;

[0026] Figure 6 is a schematic diagram of a scattering distribution function provided by an embodiment of the present invention;

[0027] Figures 7 to 8 This is a diagram showing the accuracy verification result of the modified EBRDF radiance calculation provided by an embodiment of the present invention;

[0028] Fig. 9is a hardware architecture diagram of a computing device provided by an embodiment of the present invention;

[0029] Fig.10 It is a structural diagram of a device for parameterizing the physical properties of ground-air light scattering provided by one embodiment of the present invention. DETAILED DESCRIPTION

[0030] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0031] As mentioned above, Figure 1 As shown in the figure, the ground-air light scattering characteristics are greatly affected by the ground surface, atmosphere and clouds, while the existing radiation transfer algorithms mostly average the optical characteristics of the ground features, clouds, aerosols and atmosphere, and have poor adaptability to scenes with large differences between the ground and the atmosphere, thus failing to meet the actual needs of high precision under the conditions of dynamic evolution of time and space.

[0032] To address this problem, the inventors proposed that the complex coupled radiation transmission process of the earth-atmosphere system can be reduced in dimensionality, and the environmental characteristic parameters and characteristic values ​​of the light scattering properties in different detection scenarios can be extracted and parameterized by constructing an equivalent scattering distribution function model. High precision can be guaranteed while meeting the conditions of the spatiotemporal dynamic evolution of the optical properties of surface characteristics, clouds, aerosols and atmosphere. At the same time, the equivalent scattering distribution model can significantly improve the computational efficiency by extracting key characteristic parameters for dimensionality reduction and simplifying the physical process into a functional relationship of parameters.

[0033] The specific implementation of the above concept is described below.

[0034] Please refer to Figure 2 The embodiment of the present invention provides a method for parameterizing the physical properties of ground-air light scattering, the method comprising:

[0035] Step 100, dividing the detection scene of the detector according to the relative positions of the sun, the earth and the detector; the detection scene includes an edge detection scene and a non-edge detection scene;

[0036] Step 102, respectively obtaining geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in a plurality of edge detection scenarios and non-edge detection scenarios;

[0037] Step 104, respectively traverse array sets of geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in two detection scenes, and perform coupled radiation transmission simulation on each array set to obtain a first massive database and a second massive database;

[0038] Step 106: Obtain equivalent scattering distribution function models in edge detection scenarios and non-edge detection scenarios respectively based on the first massive database and the second massive database, and perform parameterized characterization using the ground-air light scattering properties.

[0039] In the embodiment of the present invention, the detection scene of the detector is first divided into a limbal detection scene or a non-limbal detection scene according to the relative positions of the sun, the earth and the detector, and then a plurality of geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in each detection scene are arranged and combined to form an array set, and each array set is simulated to obtain a first massive database and a second massive database, respectively, and finally, the equivalent scattering distribution function models in the two detection scenes are constructed based on the two massive databases. Based on the equivalent scattering distribution function model in the two detection scenes, the scattered radiance distribution in the hemispherical space under any ground-atmosphere background can be output according to the environmental characteristic parameter values ​​of the surface, cloud clusters, aerosols, etc. with spatial longitude and latitude distribution, combined with the geometric characteristic parameter values; in this way, while meeting the actual needs of the spatiotemporal dynamic evolution of the ground-atmosphere system, it has high accuracy.

[0040] For step 100:

[0041] In some embodiments, dividing the detection scene of the detector according to the relative positions of the sun, the earth, and the detector includes:

[0042] Based on the relative positions of the sun, the earth and the detector, a discriminant of the intersection of the detector's detection light and the earth is constructed;

[0043] Dividing the detection scene of the detector according to the intersection discriminant;

[0044] If the intersection discriminant is greater than zero, the detection light of the detector intersects with the earth, and the current detection scene is divided into a non-edge detection scene;

[0045] If the intersection discriminant is less than zero, the detection light of the detector has no intersection with the earth, and the current detection scene is divided into an edge detection scene.

[0046] In the embodiment of the present invention, in the earth-fixed coordinate system, the origin of the coordinate system is taken as the center of the earth, and according to the relative positions of the sun, the earth and the detector, a discriminant of the intersection of the detection light of the detector and the earth is constructed:

[0047] Δ=(p·d) 2 -(p·pRE 2 )

[0048] Where p is the coordinate position of the detector, d is the direction of the detector, R E is the radius of the Earth.

[0049] In both edge-proximal and non-edge-proximal detection scenarios, the geometric characteristic parameters characterize the spatial position relationship between the incident light and the scattered light. In the non-edge-proximal detection scenario: the solar zenith angle represents the direction of the incident light, the observation zenith angle represents the direction of the scattered light, and the relative azimuth angle represents the relative position relationship between the incident and scattered light relative to the reflection surface. Therefore, the non-edge-proximal detection scenario can use the above three geometric characteristic parameters to characterize the spatial position relationship between the incident light and the scattered light. In the edge-proximal detection scenario: the solar zenith angle represents the direction of the incident light, the tangential height represents the spatial height position of the scattered light, and the relative azimuth angle represents the relative position relationship between the incident and scattered light relative to the reflection surface. Therefore, the edge-proximal detection scenario can use the above three geometric characteristic parameters to characterize the spatial position relationship between the incident light and the scattered light.

[0050] Considering that the geometric characteristic parameters corresponding to the two detection scenarios are different, therefore, in the embodiment of the present invention, according to the intersection of the detection light of the detector and the earth's surface during the calculation of the background light radiation transmission of the ground atmosphere, the detection scene can be divided into an edge detection scene or a non-edge detection scene. The intersection is a non-edge detection scene, and the non-intersection is an edge detection scene. Figure 3 As shown in the figure, in the non-edge detection scenario, there is an intersection between the detection light and the earth surface, and the scattered light is distributed in the hemispherical space; Figure 4 As shown, in the edge detection scenario, there is no intersection between the detection light and the earth's surface.

[0051] For steps 102 to 104:

[0052] Taking into account that the distribution range of scattered light is different in different detection scenarios, therefore, in an embodiment of the present invention, firstly, the distribution range of geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in two different detection scenarios is observed through satellite observation data, and characteristic values ​​are extracted at intervals within the distribution range to achieve full coverage of complex earth-atmosphere system detection scene conditions; thereafter, according to the characteristic parameter values ​​(including geometric characteristic parameter values ​​and environmental characteristic parameter values) respectively extracted in the two detection scenarios, all array sets of characteristic parameter values ​​are traversed respectively, and coupled radiation transmission simulation (such as Monte Carlo simulation transmission simulation) is performed on each obtained array set, and finally, a massive database containing geometric characteristic parameter values, environmental characteristic parameter values ​​and light scattering radiation brightness values ​​in different detection scenarios can be obtained.

[0053] In some embodiments, the geometric characteristic parameters in the limb detection scenario include the solar zenith angle, the relative azimuth angle, and the cut height; the geometric characteristic parameters in the non-limb detection scenario include the solar zenith angle, the relative azimuth angle, and the observation zenith angle;

[0054] The environmental characteristic parameters in the edge detection scenario and the non-edge scenario both include surface environmental parameters, aerosol environmental parameters and cloud environmental parameters; wherein the surface environmental parameters include surface reflectivity, the aerosol environmental parameters include aerosol type and aerosol optical thickness, and the cloud environmental parameters include cloud top height and cloud optical thickness.

[0055] The surface albedo characterizes the ratio of the surface reflected energy to the total incident energy. Different types of surfaces have different albedo, which affects the intensity and spatial distribution of scattered light. Therefore, in the embodiment of the present invention, according to the statistical distribution information of the surface type, the surface types such as forests, snow, lakes, seas, and deserts are extracted, and the surface albedo is extracted as a characteristic physical quantity to distinguish different surfaces. In addition, different types of aerosols have different optical properties and different scattering and absorption capabilities for light. Aerosol optical thickness reflects the extinction effect of the entire layer of aerosols, and highly summarizes the effects of aerosol concentration, particle size distribution, etc. on light scattering. Therefore, in the embodiment of the present invention, according to the statistical distribution information of aerosols, aerosol types such as villages, towns, and oceans are extracted, and aerosol types and aerosol optical thickness are extracted as characteristic physical quantities to distinguish different aerosols. At the same time, the cloud top height determines the path length between the incident light passing through the atmosphere and reaching the scattering surface. The longer the optical path, the more obvious the scattering effect. Cloud optical thickness describes the absorption and scattering ability of the cloud layer to light, and highly summarizes the effects of factors such as cloud particle size and structure on light scattering. Therefore, in the embodiment of the present invention, cloud types such as stratus clouds, cumulus clouds, cirrus clouds, and nimbostratus clouds are extracted based on cloud statistical distribution information, and cloud top height and cloud optical thickness are extracted as characteristic physical quantities to distinguish different cloud clusters.

[0056] Regarding step 106:

[0057] In the embodiment of the present invention, based on the massive database already available in step 104, the functional relationship between the light scattering characteristics and the physical quantity of the characteristic parameters can be obtained, and the equivalent scattering distribution function models under the two detection scenarios can be obtained respectively. By inputting specific geometric characteristic parameters and environmental characteristic parameters into the model, the light scattering intensity under the conditions can be obtained, thereby realizing the parametric characterization of the light scattering physical properties for complex earth-atmosphere systems.

[0058] Furthermore, as shown in Table 2, according to the statistical distribution law of meteorological factors, it is found that within the range of 0-20 km, the surface, aerosol, cloud type, etc. have a greater coupling effect on light scattering, and after exceeding 20 km, the propagation path of the detection light in the atmosphere is long, and after multiple scattering and attenuation, it has a greater impact on the detection value. Therefore, in the embodiment of the present invention, Figure 5 As shown in the figure, by establishing a virtual spherical shell in the area 20km above the surface, the equivalent scattering distribution function model is established at the virtual spherical shell height (20km), and by performing complex coupled high-precision radiation transmission simulation in the atmospheric range from the surface to the virtual spherical shell height (0-20km), the influence of various surface, aerosol and cloud types on light scattering is fully considered, thereby obtaining the hemispherical space light scattering distribution under various ground-atmosphere backgrounds at the virtual spherical shell; and in the range of 20-180km, by simultaneously considering the attenuation and multiple scattering effects of the upper atmosphere, a correction term is introduced into the equivalent scattering distribution function model, thereby further improving the accuracy of the equivalent scattering distribution function model.

[0059] Table 1 Statistical table of atmosphere, aerosol and cloud height distribution

[0060]

[0061]

[0062] like Figure 6 As shown in the figure, the directionality of light scattering distribution in the virtual spherical shell space is determined by the solar zenith angle θ i , observation zenith angle θ r , relative azimuth Characterization:

[0063]

[0064] Where λ is the wavelength, L λ,i is the scattered radiation intensity of light with incident wavelength λ, L λ,o is the scattered radiation intensity of light with reflected wavelength λ, Ω i is the solid angle, φ=|φ i -φ o |, φ is the relative azimuth, φ o is the observation azimuth.

[0065] Therefore, in some specific implementations, when the detection altitude of the detector is 0-20 km:

[0066] The equivalent scattering distribution function model in the edge detection scenario is as follows:

[0067] f EBRDF =f(ρ,L1,COD,AT,AOD,θ i ,H,θr )

[0068] The equivalent scattering distribution function model in the non-edge detection scenario is as follows:

[0069]

[0070] In the formula, f EBRDF and f' EBRDF are the light scattering brightness output by the equivalent scattering distribution function model in the limb detection scenario and the non-limb detection scenario, ρ is the surface reflectivity, L1 is the cloud top height, COD is the cloud optical thickness, AT is the aerosol type, AOD is the aerosol optical thickness, θ i is the solar zenith angle, H is the cut height, θ r To observe the zenith angle, is the relative azimuth.

[0071] In some embodiments, when the detection altitude of the detector is 20-180 km, the equivalent scattering distribution function model needs to be corrected;

[0072] Wherein, the correction formula is as follows:

[0073] R=R EBRDF (τ+α)

[0074]

[0075] Where R is the light scattering radiation brightness value output by the modified equivalent scattering distribution function model, W / (m 2 ×sr), R EBRDF is the light scattering radiation brightness value output by the equivalent scattering distribution function model in the edge detection scene or the non-edge detection scene, W / (m 2 ×sr), τ is the atmospheric path transmittance from 20 km to the detector location, and α is the influence coefficient of the atmosphere above 20 km on the ground air light scattering.

[0076] It can be seen that the present invention proposes an equivalent ground-based light scattering (EBRDF) parameterized model to characterize the numerical distribution of scattered light in hemispherical space. In real detection scene applications, the height distribution of dominant influencing factors such as the atmosphere, aerosols, and clouds is considered, and the use altitude of the EBRDF parameterized model is proposed to be 20km. At the same time, the attenuation and multiple scattering effects of the upper atmosphere are considered, and the transmittance τ and correction factor are introduced to further improve the accuracy of the parameterized characterization model.

[0077] In order to better understand the scheme, the specific application of the equivalent scattering distribution function model (EBRDF) in the calculation of ground-air light scattering characteristics is further described in the following way.

[0078] Figure 5Characterize the influence of the ground-atmosphere background system on light scattering in a single detection direction. In this embodiment, first, facing the ground-atmosphere background detection scene, as shown in Table 2, extract the environmental characteristic parameters that have a dominant influence on the light scattering properties, namely, the surface albedo ρ, aerosol type AT, aerosol optical thickness AOD, cloud top height L1 and cloud optical thickness COD. The value range of the above characteristic physical quantities is given by the distribution of satellite observation values, which ensures the accuracy and full coverage of the characteristic values. Then, for any environmental characteristic parameter and geometric characteristic parameter, based on the known satellite detection results, extract the characteristic values ​​at intervals within the value range, and obtain a massive database through complex coupled radiation transfer simulation (such as Monte Carlo method) of the light scattering characteristic distribution. Finally, introduce the equivalent scattering distribution function model (EBRDF) to characterize the distribution characteristics of the light scattering characteristics in the virtual spherical shell space, such as Figure 6 As shown, the parametric characterization of light scattering properties under any environmental background and detection scene is achieved.

[0079] Table 2 Statistics of physical quantities of geophysical background characteristics

[0080]

[0081] According to the statistical distribution law and physical definition, the eigenvalues ​​of eight characteristic physical quantities in the non-edge detection scenario are selected (as shown in Table 3), and all array sets of the eigenvalues ​​of the eight characteristic physical quantities are traversed. The numerical distribution of ground-atmosphere light scattering is simulated by coupled radiation transmission to obtain a massive data set with high precision and full coverage. The equivalent scattering distribution function model (EBRDF) in the non-edge detection scenario is further obtained, thereby characterizing the emission distribution characteristics of scattered light in different directions under specific ground-atmosphere background conditions and incident angles (as shown in Table 4).

[0082] Table 3 Statistics of characteristic parameters of non-edge detection model

[0083]

[0084] Table 4 Output data of equivalent scattering distribution function model in non-limb detection scenario (×10 -2 )

[0085]

[0086] Table 4 shows the hemispherical spatial distribution of scattered light in the non-edge detection scenario when the surface albedo is 0.5, the cloud top height is 8 km, the cloud optical thickness is 10, the rural aerosol, the aerosol optical thickness is 0.03, and the solar zenith angle is 30°. The data in the table are the results of the equivalent scattering distribution function model (EBRDF) in the non-edge detection scenario. It can be seen from Table 4 that light scattering has a significant change characteristic with angle, and multiplying the solar constant is the radiance result.

[0087] Furthermore, according to the statistical distribution law and physical definition, the eigenvalues ​​of eight characteristic physical quantities in the edge detection scenario are selected (as shown in Table 5), and all array sets of the eigenvalues ​​of the eight characteristic physical quantities are traversed. The numerical distribution of ground-atmosphere light scattering is simulated by coupled radiation transmission to obtain a massive data set with high precision and full coverage, and further the equivalent scattering distribution function model (EBRDF) in the edge detection scenario is obtained, thereby characterizing the emission distribution characteristics of scattered light in different directions under specific ground-atmosphere background conditions and incident angles (as shown in Table 6).

[0088] Table 5 Statistics of characteristic parameters of edge detection model

[0089]

[0090] Table 6 Output data of the equivalent scattering distribution function model in the limb detection scenario (×10 -2 )

[0091]

[0092]

[0093] Table 6 shows the hemispherical spatial distribution of scattered light under the conditions of limb detection when the surface albedo is 0.7, the cloud top height is 2 km, the cloud optical thickness is 3, the rural aerosol, the aerosol optical thickness is 0.03, and the solar zenith angle is 60°. The data in the table are the results of the equivalent scattering distribution function model (EBRDF) under the limb detection scenario. As shown in Table 6, the cut height is the main influencing factor of the light scattering properties, and multiplying the solar constant is the radiance result.

[0094] Furthermore, the calculation accuracy of the parameterized characterization method of light scattering properties and the complex radiation transmission calculation method in the embodiment of the present invention are verified. After the correction method in the embodiment of the present invention is used to correct the equivalent scattering distribution function model, at all detection angles, Figure 7 The average relative deviation is about 0.588%, and the maximum relative deviation is about 1.754%; Figure 8 The average relative deviation is about 0.257%, and the maximum relative deviation is about 0.502%. The results are consistent with the complex coupled radiation transfer model. The verification results are shown in Figures 7 to 8 .

[0095] In summary, in the embodiments of the present invention, by fully considering the influence of the surface, atmosphere, aerosol, cloud environmental parameters and detection geometric parameters on light scattering, by extracting key characteristic parameters for dimensionality reduction processing, the physical process is simplified into a functional relationship of parameters, and an equivalent scattering distribution function model (EBRDF) is proposed to reduce the dimensionality of the complex coupled radiation transmission process, thereby realizing the multi-dimensional light scattering characteristics of the earth-atmosphere system in time, space, spectrum and direction, solving the problem that the current earth-atmosphere background light scattering averages the physical properties and only considers static scenes, while satisfying the conditions for the spatiotemporal dynamic evolution of the optical properties of the surface characteristics, clouds, aerosols and atmosphere, the amount of calculation is reduced, and the high precision of the complex coupled radiation transmission calculation process is retained.

[0096] like Fig. 9 , Fig.10 As shown, an embodiment of the present invention provides a device for parameterizing the physical properties of ground-air light scattering. The device embodiment can be implemented by software, or by hardware or a combination of software and hardware. From the hardware level, Fig. 9 As shown, it is a hardware architecture diagram of a computing device where a device for parameterizing the physical properties of ground-air light scattering provided by an embodiment of the present invention is located. Fig. 9 In addition to the processor, memory, network interface, and non-volatile memory shown in the figure, the computing device in which the device is located in the embodiment may also generally include other hardware, such as a forwarding chip responsible for processing messages, etc. Taking software implementation as an example, Fig.10 As shown, as a device in a logical sense, the CPU of the computing device in which it is located reads the corresponding computer program in the non-volatile memory into the memory and runs it. This embodiment provides a device for parameterizing the physical properties of ground-air light scattering, and the device includes:

[0097] A division unit 901 is used to divide the detection scene of the detector according to the relative positions of the sun, the earth and the detector; the detection scene includes an edge detection scene and a non-edge detection scene;

[0098] An acquisition unit 902 is used to respectively acquire geometric feature parameter values ​​and environmental feature parameter values ​​in a plurality of edge detection scenarios and non-edge detection scenarios;

[0099] The simulation unit 903 is used to traverse the array sets of geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in two detection scenes respectively, and perform coupled radiation transmission simulation on each array set to obtain a first massive database and a second massive database respectively;

[0100] The characterization unit 904 is used to obtain equivalent scattering distribution function models in edge detection scenarios and non-edge detection scenarios respectively according to the first massive database and the second massive database, and perform parameterized characterization based on the ground-air light scattering properties.

[0101] In an embodiment of the present invention, the division unit 901 can be used to execute step 100 in the above method embodiment, the acquisition unit 902 can be used to execute step 102 in the above method embodiment, the simulation unit 903 can be used to execute step 104 in the above method embodiment, and the characterization unit 904 can be used to execute step 106 in the above method embodiment.

[0102] In one embodiment of the present invention, in the division unit 901, when dividing the detection scene of the detector according to the relative positions of the sun, the earth and the detector, it includes:

[0103] Based on the relative positions of the sun, the earth and the detector, a discriminant of the intersection of the detector's detection light and the earth is constructed;

[0104] Dividing the detection scene of the detector according to the intersection discriminant;

[0105] If the intersection discriminant is greater than zero, the detection light of the detector intersects with the earth, and the current detection scene is divided into a non-edge detection scene;

[0106] If the intersection discriminant is less than zero, the detection light of the detector has no intersection with the earth, and the current detection scene is divided into an edge detection scene.

[0107] In one embodiment of the present invention, in the division unit 901, the intersection discriminant of the detection light of the detector and the earth is as follows:

[0108] Δ=(p·d) 2 -(p·pR E 2 )

[0109] Where p is the coordinate position of the detector, d is the direction of the detector, R E is the radius of the Earth.

[0110] In one embodiment of the present invention, the geometric characteristic parameters in the edge detection scene include the solar zenith angle, the relative azimuth angle and the cutting height; the geometric characteristic parameters in the non-edge detection scene include the solar zenith angle, the relative azimuth angle and the observation zenith angle;

[0111] The environmental characteristic parameters in the edge detection scenario and the non-edge scenario both include surface environmental parameters, aerosol environmental parameters and cloud environmental parameters; wherein the surface environmental parameters include surface reflectivity, the aerosol environmental parameters include aerosol type and aerosol optical thickness, and the cloud environmental parameters include cloud top height and cloud optical thickness.

[0112] In one embodiment of the present invention, in the characterization unit 904, when the detection altitude of the detector is 0-20 km:

[0113] The equivalent scattering distribution function model in the edge detection scenario is as follows:

[0114] f EBRDF =f(ρ,L1,COD,AT,AOD,θ i ,H,θ r )

[0115] The equivalent scattering distribution function model in the non-edge detection scenario is as follows:

[0116]

[0117] In the formula, f EBRDF and f' EBRDF are the light scattering brightness output by the equivalent scattering distribution function model in the limb detection scenario and the non-limb detection scenario, ρ is the surface reflectivity, L1 is the cloud top height, COD is the cloud optical thickness, AT is the aerosol type, AOD is the aerosol optical thickness, θ i is the solar zenith angle, H is the cut height, θ r To observe the zenith angle, is the relative azimuth.

[0118] In one embodiment of the present invention, when the detection altitude of the detector is 20-180 km, the equivalent scattering distribution function model needs to be corrected;

[0119] Wherein, the correction formula is as follows:

[0120] R=R EBRDF (τ+α)

[0121]

[0122] Where R is the light scattering radiation brightness value output by the modified equivalent scattering distribution function model, R EBRDF is the light scattering radiation brightness value output by the equivalent scattering distribution function model in the limb detection scenario or the non-limb detection scenario, τ is the atmospheric path transmittance from 20 km to the detector location, and α is the influence coefficient of the atmosphere above 20 km on the ground air light scattering.

[0123] It is to be understood that the structure illustrated in the embodiment of the present invention does not constitute a specific limitation on a device for characterizing the physical properties of ground-air light scattering parameters. In other embodiments of the present invention, a device for characterizing the physical properties of ground-air light scattering parameters may include more or fewer components than shown in the figure, or combine some components, or split some components, or arrange the components differently. The components shown in the figure may be implemented in hardware, software, or a combination of software and hardware.

[0124] The information interaction, execution process and other contents between the modules in the above-mentioned device are based on the same concept as the embodiment of the method of the present invention. For the specific contents, please refer to the description in the embodiment of the method of the present invention, and no further description is given here.

[0125] An embodiment of the present invention further provides a computing device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, a method for parameterizing the physical properties of ground-air light scattering in any embodiment of the present invention is implemented.

[0126] An 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, the processor executes a method for parameterizing the physical properties of ground-air light scattering in any embodiment of the present invention.

[0127] Specifically, a system or device equipped with a storage medium can be provided, on which software program code that implements the functions of any of the above-mentioned embodiments is stored, and a computer (or CPU or MPU) of the system or device can be enabled to read and execute the program code stored in the storage medium.

[0128] In this case, the program code itself read from the storage medium can realize the function of any one of the above-mentioned embodiments, and thus the program code and the storage medium storing the program code constitute a part of the present invention.

[0129] The storage medium embodiments for providing the program code include a floppy disk, a hard disk, a magneto-optical disk, an optical disk (such as CD-ROM, CD-R, CD-RW, DVD-ROM, DVD-RAM, DVD-RW, DVD+RW), a magnetic tape, a non-volatile memory card, and a ROM. Alternatively, the program code can be downloaded from a server computer by a communication network.

[0130] In addition, it should be clear that the functions of any of the above embodiments can be implemented not only by executing the program code read by the computer, but also by enabling an operating system operating on the computer to complete part or all of the actual operations based on instructions from the program code.

[0131] In addition, it can be understood that the program code read from the storage medium is written to a memory provided in an expansion board inserted into the computer or to a memory provided in an expansion module connected to the computer, and then based on the instructions of the program code, a CPU installed on the expansion board or expansion module is enabled to perform part or all of the actual operations, thereby realizing the functions of any of the above-mentioned embodiments.

[0132] An embodiment of the present application also provides a computer-readable storage medium, on which is stored at least one instruction, at least one program, code set or instruction set, and the at least one instruction, at least one program, code set or instruction set is loaded and executed by a processor to implement a method for parameterizing the physical properties of ground-air light scattering provided in the above-mentioned method embodiments.

[0133] It should be noted that, in this article, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device.

[0134] A person of ordinary skill in the art can understand that all or part of the steps of implementing the above method embodiments can be completed by hardware related to program instructions, and the aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it executes the steps of the above method embodiments; and the aforementioned storage medium includes: ROM, RAM, magnetic disk or optical disk, etc., various media that can store program codes.

[0135] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for parameterizing the physical properties of ground-air light scattering, characterized in that: include: According to the relative positions of the sun, the earth and the detector, the detection scenes of the detector are divided; the detection scenes include edge detection scenes and non-edge detection scenes; Obtaining geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in several edge detection scenarios and non-edge detection scenarios respectively; traversing array sets of geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in two detection scenes respectively, and performing coupled radiation transmission simulation on each array set to obtain a first massive database and a second massive database respectively; Based on the first massive database and the second massive database, equivalent scattering distribution function models in edge detection scenarios and non-edge detection scenarios are obtained respectively, and parameterized representation is performed using the ground-air light scattering properties.

2. The method according to claim 1, characterized in that The dividing of the detection scenes of the detector according to the relative positions of the sun, the earth and the detector includes: Based on the relative positions of the sun, the earth and the detector, a discriminant of the intersection of the detector's detection light and the earth is constructed; Dividing the detection scene of the detector according to the intersection discriminant; If the intersection discriminant is greater than zero, the detection light of the detector intersects with the earth, and the current detection scene is divided into a non-edge detection scene; If the intersection discriminant is less than zero, the detection light of the detector has no intersection with the earth, and the current detection scene is divided into an edge detection scene.

3. The method according to claim 2, characterized in that The intersection discriminant of the detection light of the detector and the earth is as follows: Δ=(p·d) 2 -(p·pR E 2 ) Where p is the coordinate position of the detector, d is the direction of the detector, R E is the radius of the Earth.

4. The method according to claim 1, characterized in that The geometric characteristic parameters in the edge detection scenario include the solar zenith angle, relative azimuth angle and cut height; the geometric characteristic parameters in the non-edge detection scenario include the solar zenith angle, relative azimuth angle and observation zenith angle; The environmental characteristic parameters in the edge detection scenario and the non-edge scenario both include surface environmental parameters, aerosol environmental parameters and cloud environmental parameters; wherein the surface environmental parameters include surface reflectivity, the aerosol environmental parameters include aerosol type and aerosol optical thickness, and the cloud environmental parameters include cloud top height and cloud optical thickness.

5. The method according to claim 1 or 4, characterized in that: When the detection altitude of the detector is 0-20km: The equivalent scattering distribution function model in the edge detection scenario is as follows: f EBRDF =f(ρ,L1,COD,AT,AOD,θ i ,H,θ r ) The equivalent scattering distribution function model in the non-edge detection scenario is as follows: In the formula, f EBRDF and f' EBRDF are the light scattering brightness output by the equivalent scattering distribution function model in the limb detection scenario and the non-limb detection scenario, ρ is the surface reflectivity, L1 is the cloud top height, COD is the cloud optical thickness, AT is the aerosol type, AOD is the aerosol optical thickness, θ i is the solar zenith angle, H is the cut height, θ r To observe the zenith angle, is the relative azimuth.

6. The method according to claim 5, characterized in that When the detection altitude of the detector is 20-180 km, the equivalent scattering distribution function model needs to be corrected; Wherein, the correction formula is as follows: R=R EBRDF ·(t+a) Where R is the light scattering radiation brightness value output by the modified equivalent scattering distribution function model, R EBRDF is the light scattering radiation brightness value output by the equivalent scattering distribution function model in the limb detection scenario or the non-limb detection scenario, τ is the atmospheric path transmittance from 20 km to the detector location, and α is the influence coefficient of the atmosphere above 20 km on the ground air light scattering.

7. A device for parameterizing the physical properties of ground-air light scattering, characterized in that: include: A division unit, used for dividing the detection scene of the detector according to the relative positions of the sun, the earth and the detector; The detection scenarios include edge detection scenarios and non-edge detection scenarios; An acquisition unit, used to respectively acquire geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in a plurality of edge detection scenarios and non-edge detection scenarios; A simulation unit, used for respectively traversing array sets of geometric characteristic parameter values ​​and environmental characteristic parameter values ​​in two detection scenes, and performing coupled radiation transmission simulation on each array set to obtain a first massive database and a second massive database respectively; The characterization unit is used to obtain the equivalent scattering distribution function model in the edge detection scenario and the non-edge detection scenario respectively according to the first massive database and the second massive database, and perform parameterized characterization with the ground-air light scattering properties.

8. A computing device, comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the method according to any one of claims 1 to 6 is implemented.

9. A computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to execute the method according to any one of claims 1 to 6.

10. A computer program product, characterized in that The method comprises a computer program, wherein when the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Images