Radiance determination method and apparatus, computer device, and readable storage medium

Abstract

The application provides a radiance determination method and device, computer equipment and a readable storage medium, wherein the radiance determination method comprises: obtaining a radiance determination parameter; determining a boundary condition parameter according to the radiance determination parameter; determining atmospheric scattering phase matrix information according to the radiance determination parameter; and determining radiance according to the radiance determination parameter, the atmospheric scattering phase matrix information and boundary condition information. When determining the radiance, the method can determine complete Stokes vector information, including all cosine and sine mode information of four Stokes components; can get rid of the property assumption of the scattering phase matrix and the lower boundary reflection matrix in the existing model, and achieve the effect of completely coupling atmospheric scattering and boundary reflection under general conditions; and removes the interpolation process of radiation transfer solving, and can directly output the radiance result at any zenith angle.

Description

Technical Field

[0001] This invention relates to the field of satellite radiative transfer modeling technology, and in particular to a method and apparatus for determining radiance, a computer device, and a readable storage medium. Background Technology

[0002] Radiative transfer models (RTMs) have been widely used in space sensor simulations, satellite data assimilation, and remote sensing algorithms. Doubling and Adding methods have also been developed as radiative transfer solution schemes and applied to fast radiative transfer models such as CRTM and ARMS. The discrete ordinate method has also been applied to solve the vector radiative transfer equation (RTE). Unlike the scalar discrete ordinate radiative transfer (DISORT) model developed by Stamnes et al., the vector discrete ordinate radiative transfer (VDISORT) model considers the polarization characteristics of radiance. Typically, the discrete ordinate method uses the Fourier expansion method and a Gaussian integral summation scheme to approximate the double integral in the RTE through summation. The Fourier expansion method approximates the Stokes vector (I... l ,I r ,I u ,I v ) T The element is decomposed into cosine and sine modes at relative azimuth angles, that is:

[0003]

[0004] in, sinusoidal mode of radiance in traditional RTM Sum and cosine modes That is The combination of elements summed over their Fourier orders is:

[0005]

[0006]

[0007] For a fully polarized RTM, the cosine modes of the first two components of the Stokes vector are combined with the sinusoidal modes of the third and fourth components to form a new vector. Then, a cosine mode of the ordinary differential RTM is constructed for solving. This combination method can be well integrated with the properties of the Mie scattering phase matrix of spherical particles in the atmosphere. Other Stokes vector elements constitute the sinusoidal mode of the RTM. Past research on RTM has largely focused on solving the cosine model equations of RTM, which effectively addresses the radiative transfer problem of non-polarized radiation sources in the atmosphere undergoing Mie or Rayleigh scattering. However, research on solving the sinusoidal RTM model and its contribution to total radiance remains lacking, particularly for radiative transfer problems involving polarized sources and surface boundary conditions. Furthermore, the solution process for the cosine RTM model relies on assumptions about the properties of the atmospheric scattering phase matrix and the boundary reflection matrix, specifically assuming that the phase and reflection matrices possess the same azimuth symmetry as the Mie scattering phase matrix of spherical particles. This introduces limitations and errors to the practical application of RTM.

[0008] Furthermore, in traditional discrete ordinate radiative transfer (RTE) solutions, the incident and scattered zenith angles vary within the same set of Gaussian points in the same Gaussian integral summation scheme. This means that the RTE can only be solved at Gaussian points. Then, interpolation methods are used to calculate the radiance at zenith angles other than the Gaussian points. When the number of Gaussian points is small, the interpolation process can introduce significant errors into the radiance results. Therefore, eliminating or mitigating the errors introduced by the interpolation process is one of the issues that needs to be addressed in the determination of radiative transfer. Summary of the Invention

[0009] This invention provides a radiance determination method that fully couples atmospheric scattering and boundary reflection under normal conditions, directly outputting complete cosine and sinusoidal mode radiance results at any zenith angle. The radiance determination method includes:

[0010] Obtain radiance determination parameters;

[0011] Determine the boundary condition parameters based on the radiance parameters;

[0012] The atmospheric scattering phase matrix information is determined based on the radiance parameters.

[0013] The radiance is determined based on the parameters of radiance determination, atmospheric scattering phase matrix information, and boundary condition information.

[0014] In specific implementation, the radiance determination parameters include satellite parameters, atmospheric state parameters, lower boundary state parameters, and radiance output direction parameters.

[0015] In specific implementation, determining the boundary condition parameters based on the radiance determination parameters includes:

[0016] The thermal radiation radiance, reflection matrix, and emissivity matrix at the upper and lower boundaries of the model and the atmospheric stratification interface are calculated based on the determined parameters of radiance.

[0017] In specific implementation, determining the atmospheric scattering phase matrix information based on the radiance determination parameters includes:

[0018] Input the radiance determination parameters into the theoretical formula for the scattering phase matrix;

[0019] The atmospheric scattering phase matrix information is determined based on the output of the theoretical formula for the scattering phase matrix.

[0020] In specific implementation, determining the radiance based on radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information includes:

[0021] Obtain the radiative transfer model;

[0022] Input the radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information into the radiative transfer model;

[0023] The radiance is determined based on the output of the radiative transfer model.

[0024] In practice, the radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information are input into the radiative transfer model, and calculations are performed according to the following formula:

[0025]

[0026] Where I represents radiance (I l ,I r ,I u ,I v ); μ represents the zenith angle cosine; Q is the sum of single scattering and medium radiation; M represents the scattering phase matrix; c1 and c2 represent coefficients; τ represents the optical thickness; m represents the Fourier expansion order; ω(τ) represents the single scattering albedo; w represents the weight of the Gaussian point in the Gaussian integral summation framework under 2N flux numbers; l, r, u, and v represent the Stokes vector elements.

[0027] The present invention also provides a radiance determination device, the radiance determination device comprising:

[0028] The parameter acquisition module is used to acquire radiance determination parameters;

[0029] The atmospheric scattering phase matrix determination module is used to determine atmospheric scattering phase matrix information based on radiance determination parameters.

[0030] The boundary condition determination module is used to determine the boundary condition parameters based on the radiance determination parameters.

[0031] The radiance output module is used to determine the radiance based on the radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information.

[0032] The present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the radiance determination method.

[0033] The present invention also provides a computer-readable storage medium storing a computer program for the radiance determination method.

[0034] This invention provides a method, apparatus, computer device, and readable storage medium for determining radiance. The method includes: acquiring radiance determination parameters; determining boundary condition parameters based on the radiance determination parameters; determining atmospheric scattering phase matrix information based on the radiance determination parameters; and determining radiance based on the radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information. This method, when determining radiance, can determine complete Stokes vector information, including all cosine and sine mode information of the four Stokes components; it can eliminate the assumptions of existing models regarding the properties of the scattering phase matrix and the lower boundary reflection matrix, achieving a complete coupling effect between atmospheric scattering and boundary reflection under general conditions; it removes the interpolation process for solving radiative transfer, and can directly output radiance results at any zenith angle. Attached Figure Description

[0035] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some specific embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. In the drawings:

[0036] Figure 1 This is a flowchart illustrating a method for determining radiance according to a specific embodiment of the present invention;

[0037] Figure 2 This is a flowchart illustrating the process of determining atmospheric scattering phase matrix information by determining parameters according to a specific embodiment of the present invention;

[0038] Figure 3 This is a flowchart illustrating the process of solving a radiative transfer model according to a specific embodiment of the present invention;

[0039] Figure 4 This is a schematic diagram comparing the simulation results of the radiative transfer model with the Rayleigh baseline according to a specific embodiment of the present invention;

[0040] Figure 5This is a schematic diagram comparing the simulation results of the radiative transfer model with the L13 baseline according to a specific embodiment of the present invention;

[0041] Figure 6 This is a schematic diagram of the simulation results of the radiative transfer model of L13 under a polarization source according to a specific embodiment of the present invention;

[0042] Figure 7 This is a schematic diagram of the radiance determination device according to a specific embodiment of the present invention. Detailed Implementation

[0043] To make the objectives, technical solutions, and advantages of the specific embodiments of the present invention clearer, the specific embodiments of the present invention will be further described in detail below with reference to the accompanying drawings. Here, the illustrative specific embodiments and their descriptions are used to explain the present invention, but are not intended to limit the present invention.

[0044] like Figure 1 As shown, this invention provides a radiance determination method to fully couple the effects of atmospheric scattering and boundary reflection under normal conditions, and directly output complete cosine and sinusoidal mode radiance results at any zenith angle. The radiance determination method includes:

[0045] 101: Obtain radiance determination parameters;

[0046] 102: Determine the boundary condition parameters based on the radiance parameters;

[0047] 103: Determine atmospheric scattering phase matrix information based on radiance parameters;

[0048] 104: Determine the radiance based on the parameters, atmospheric scattering phase matrix information, and boundary condition information.

[0049] In practice, there are multiple implementation schemes for selecting radiance determination parameters. For example, in order to effectively obtain radiance, the radiance determination parameters may include satellite parameters, atmospheric state parameters, lower boundary state parameters, and radiance output direction parameters.

[0050] In specific implementations, there can be multiple implementation schemes for determining boundary condition parameters based on radiance determination parameters. For example, step 102: determining boundary condition parameters based on radiance determination parameters can include:

[0051] The thermal radiation radiance, reflection matrix, and emissivity matrix at the upper and lower boundaries of the model and the atmospheric stratification interface are calculated based on the determined parameters of radiance.

[0052] In practice, there are multiple implementation schemes for determining the atmospheric scattering phase matrix information based on radiance parameters, for example, such as... Figure 2As shown, step 103: determining the atmospheric scattering phase matrix information based on the radiance determination parameters may include:

[0053] 201: Input the radiance determination parameters into the theoretical formula for the scattering phase matrix;

[0054] 202: Determine the atmospheric scattering phase matrix information based on the output of the theoretical formula for the scattering phase matrix.

[0055] In practice, there are multiple implementation schemes for determining radiance based on radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information. For example, Figure 3 As shown, step 104: determining the radiance based on the radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information may include:

[0056] 301: Obtain the radiative transfer model;

[0057] 302: Input the radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information into the radiative transfer model.

[0058] 303: Determine the radiance based on the output of the radiative transfer model.

[0059] Furthermore, in step 302, the radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information are input into the radiative transfer model, and calculations can be performed based on the following formula principle:

[0060]

[0061] Where I represents radiance (I l ,I r ,I u ,I v μ represents the zenith angle cosine; Q represents the sum of single scattering and medium radiation; M represents the scattering phase matrix; c1 and c2 represent coefficients; τ represents the optical thickness; m represents the Fourier expansion order; ω(τ) represents the single scattering albedo; w represents the weight of the Gaussian point in the Gaussian integral summation framework under 2N flux numbers; l, r, u, and v represent Stokes vector elements.

[0062] Furthermore, the framework for Gaussian integral summation is as follows:

[0063]

[0064] By using a new Gaussian integral summation framework, the radiative transfer model can directly output radiance results at any zenith angle, avoiding errors caused by the interpolation process.

[0065] In practice, there are multiple implementation schemes for constructing the radiative transfer model, for example:

[0066] The satellite radiative transfer equation is in the form of:

[0067]

[0068] in, It is radiance, expressed using the Stokes vector (I l ,I r ,I u ,I v ) T It is expressed as a function of optical thickness τ, zenith cosine μ, and azimuth angle φ. For upward radiation, μ is positive; for downward radiation, μ is negative. ω(τ) is the single-scattering albedo corresponding to the atmospheric stratification; M(τ,μ,φ;μ',φ') is the scattering phase matrix, and (μ',φ') and (μ,φ) are the incident and exit directions; S b It is a radiation source with a direction of (-μ0, φ0); S t (τ) represents atmospheric thermal radiation, and B[T(τ)] represents the Planck function at temperature T(τ).

[0069] Perform a Fourier expansion on the azimuth angle:

[0070]

[0071]

[0072] in At the same time, applying the Gaussian integral summation method, we obtain:

[0073]

[0074]

[0075] in, sinusoidal modes of radiance in existing radiative transfer models Sum and cosine modes That is The combination of elements summed at each of the Fourier orders; μ j ,w j These are the Gaussian integration points and their corresponding integration weights; N is the number of Gaussian integration points of μ in the interval [0,1], and the number of Gaussian integration points and their corresponding weights in the interval [-1,0] using μ. -j =-μ j ,w -j =w j get.

[0076] Let M(τ,μ,φ;μ′,φ′)=(m pq(p,q=1,2,3,4), using the Mie scattering phase matrix M of spherical particles mie The property of (τ,μ,φ;μ′,φ′) is that its top-left and bottom-right four elements are even functions of the relative azimuth angle φ′-φ, and its bottom-left and top-right four elements are odd functions of the relative azimuth angle φ′-φ, that is:

[0077]

[0078] in,

[0079] Existing radiative transfer models often only consider the cosine mode of the RTM. Based on the properties of Mie scattering, equation (4a) can be rewritten as:

[0080]

[0081] Here, 2N is also referred to as the stream number of RTM; and similar.

[0082] However, equation (5) depends on the properties of the Mie scattering phase matrix. For the scattering phase matrix in the general case,

[0083]

[0084] Equation (5) becomes:

[0085]

[0086] Similarly, the RTM sinusoidal mode should be:

[0087]

[0088] Combining the sine mode equation (7) and the cosine mode equation (6), we obtain the complete Fourier RTM equations of all orders:

[0089]

[0090] Equation (8) holds for any μ in [-1,1]. Let μ vary within the same set of Gaussian integration points. For each μ, equation (8) has 8 corresponding equations, ultimately forming 16*N equations. At this point, under m-order RTM, there are also 16*N unknown Stokes components. Therefore, radiance can be solved using the set of ordinary differential equations corresponding to equation (8).

[0091] Equation (8) allows the RTM solution to no longer be limited to the cosine mode, but can simultaneously solve all sine and cosine components of the Stokes element to obtain complete radiance information. At the same time, Equation (8) also eliminates the Mie scattering property assumption of the scattering phase matrix, making the equation applicable to general atmospheric scattering. Based on this, a discrete ordinate radiative transfer model (VDISORT3) is constructed, which also possesses these advantages.

[0092] Furthermore, using the above method, RTM can only be solved within the set of Gaussian integration points; for radiance at other zenith angles, interpolation is still required. To allow RTM to be solved at any zenith angle, two Gaussian points with weights of 0 are added to the original set of Gaussian integration points. Since the weight w... j Since the sum of these two points is 0, their addition will not affect the left and right sides of equation (8). Let one point be equal to the cosine of any specified zenith angle to be output, and the other point be its opposite. The Gaussian points and weights in the new Gaussian integral summation scheme are:

[0093] Gaussian points: μ j =-μ -j ;w j =w -j j = 1, ..., N;

[0094] μ N+1 =|μ output |;μ -(N+1) =-μ N+1 ;w N+1 =w -(N+1) =0;

[0095] Substituting into equation (8), we obtain a new equation:

[0096]

[0097] The RTM constructed using equation (9) can achieve radiance output at any zenith angle.

[0098] Equation (9) generates a series of ordinary differential equations. The general solution and particular solution can be obtained using mathematical methods for solving ordinary differential equations, thus yielding the solution to the RTM problem. Equation (9) can be simplified to:

[0099]

[0100] in,

[0101] u = Diag[μ -(N+1) ,μ -N ,μ -(N-1) ,…,μ -1 ,μ1,μ2,…,μ(N+1) ]

[0102] e = Diag(1,…,1) 16*(N+1) ;

[0103] For any set of μ s and μ j :

[0104]

[0105]

[0106] The general solution to the system of ordinary differential equations (10) is:

[0107]

[0108] coefficient C jmp It will be determined by the boundary conditions; λ jmp and g jmp It can be obtained by solving the eigenvalue problem (13):

[0109] -λ jmp ug jmp =(e+A) m )g jmp (13)

[0110] Using the first-order polynomial approximation and the linear assumption of heat source radiation within the same layer S t,lruv (τ)≈c 0p +c 1p For τ, we can obtain the 0th-order radiative transfer equation as follows:

[0111]

[0112] By solving equation (14), the particular solution of the heat source term of the RTE can be obtained:

[0113] P tmp (τ,μ)=δ 0m [b 0p +b 1p τ], (15)

[0114] For the single scattering component of the RTE, its equation and its particular solution are:

[0115]

[0116]

[0117] Among them, t bmp Solve using equation (18):

[0118]

[0119] Ultimately, the solution to RTE is:

[0120]

[0121] Among them, P mp (τ,μ)=p mp [1:16*(N+1)].

[0122] However, the coefficient C of the general solution jmp The boundary conditions of the model need to be determined. The continuity conditions at the interface between the upper boundary and the atmosphere are the same as in previous versions of VDISORT, but further processing is required at the lower boundary to overcome the assumption of the Mie scattering phase matrix property of the lower boundary reflection matrix. The lower boundary equation of RTM is:

[0123]

[0124] Where B(+μ,φ;-μ',φ') is the polarization bidirectional reflectance distribution function (pBRDF) matrix; E(+μ,φ) is the lower boundary emissivity matrix. Using the same method as the atmospheric scattering phase matrix and equation (8), equation (20) can be written as:

[0125]

[0126] Based on equation (21), RTM can eliminate the model's assumption of the Mie scattering matrix property for the lower boundary reflection matrix. Subsequently, the lower boundary condition can be written as:

[0127]

[0128] in,

[0129] I ref (+μ)=I ref,source (+μ)+I ref,particular (+μ) (23a)

[0130]

[0131]

[0132]

[0133]

[0134] Based on the continuity assumption, the upper boundary conditions of the radiative transfer equation are:

[0135]

[0136] The boundary conditions at the interface between atmospheric strata are:

[0137]

[0138] Before solving for the coefficients of the general solution, the coefficients are subjected to the following scaling transformation:

[0139]

[0140] Subsequently, using the boundary conditions (22, 24, 25) at the upper and lower boundaries and the interface, the coefficients can be solved. Finally, the complete solution to the radiative transfer equation is obtained.

[0141] It should be noted that Rayleigh and L13 are two specific conditions for radiative transfer problems presented in previous studies, along with corresponding solution results. The solutions for Rayleigh and L13 therefore serve as benchmarks in subsequent radiative transfer studies, verifying the accuracy of the radiative transfer models. Both the Rayleigh and L13 benchmark cases assume a single-layer homogeneous atmosphere, a Lambertian surface at the lower boundary, and no heat source radiation. For the Rayleigh scattering case, the albedo at the lower boundary is set to λ0 = 0.25. Other parameters are set as follows: S b =(0.5π,0.5π,0,0) T ,ω=1,μ0=0.8,φ0=0°,τ L =1. For case L13, the parameter is set to: S b =(0.5π,0.5π,0,0) T ,ω=0.99,λ0=0.1,μ0=0.2,φ0=0°,τ L =1. The two cases were compared and verified with the reference radiance values ​​in Coulson (1960), Garcia, and Siewert (1986, 1989), respectively. The verification results are as follows: Figure 4 , Figure 5 As shown.

[0142] The Rayleigh scattering phase matrix is ​​calculated using the formula given in Chandrasekhar (1960), the Mie scattering phase matrix in L13 is calculated using the method in Haan et al. (1987), and the expansion coefficients of the Mie scattering phase matrix can be obtained from Vestrucci and Siewert (1984).

[0143] To verify the contribution and influence of the sinusoidal mode of radiance under polarization on the total radiance, the original unpolarized source was modified to a polarized source S, based on the L13 reference case. b =(0.7π,0.3π,0.2π,0.05π) TOther parameters remain unchanged. This allows us to explore the influence of sinusoidal modes on total radiance, and also demonstrates the innovation and practical value of VDISORT3 compared to traditional radiative transfer models. Specifically, Figure 6 The parameter in is set to: S b =(0.7π,0.3π,0.2π,0.05π) T and ω=0.99, λ0=0.1, μ0=0.2, φ0=0°,τ L =1.

[0144] Figures 4 to 6 The azimuth angles output by the model are all φ = 90°.

[0145] like Figure 4 , Figure 5 As shown, the VDISORT3 simulation results match the benchmark values ​​well, achieving the same effect as the previous model, which verifies the accuracy of the model. This also verifies the technical feasibility of the present invention.

[0146] like Figure 6 As shown, under polarization source conditions, the total radiance output by the model differs significantly from its sinusoidal or cosine mode. This is because previous radiative transfer models could only simulate the Stokes vector under the cosine equation. This would undoubtedly introduce significant errors into the results. VDISORT3 can simultaneously obtain the sine and cosine components of the four Stokes vectors, thus outputting complete radiance information. Furthermore, the VDISORT3 model eliminates the Mie scattering phase matrix property assumptions of previous models regarding the scattering phase matrix and lower boundary reflection matrix, making it more adaptable to atmospheric scattering and boundary reflection under general conditions, and therefore possessing broader application potential. The advantage of directly outputting radiance at any zenith angle also allows the VDISORT3 model to eliminate errors that might arise during interpolation.

[0147] like Figure 7 As shown, the present invention also provides a radiance determination device, the radiance determination device comprising:

[0148] Parameter acquisition module 701 is used to acquire radiance determination parameters;

[0149] Atmospheric scattering phase matrix determination module 702 is used to determine atmospheric scattering phase matrix information based on radiance determination parameters;

[0150] Boundary condition determination module 703 is used to determine boundary condition parameters based on radiance determination parameters;

[0151] The radiance output module 704 is used to determine the radiance based on the radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information.

[0152] The present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the radiance determination method.

[0153] The present invention also provides a computer-readable storage medium storing a computer program for the radiance determination method.

[0154] In summary, the radiance determination method, apparatus, computer equipment, and readable storage medium provided by this invention include: acquiring radiance determination parameters; determining boundary condition parameters based on the radiance determination parameters; determining atmospheric scattering phase matrix information based on the radiance determination parameters; and determining radiance based on the radiance determination parameters, atmospheric scattering phase matrix information, and boundary condition information. This method, when determining radiance, can determine complete Stokes vector information, including all cosine and sine mode information of the four Stokes components; it can eliminate the assumptions of existing models regarding the properties of the scattering phase matrix and the lower boundary reflection matrix, achieving a complete coupling effect between atmospheric scattering and boundary reflection under general conditions; it removes the interpolation process for solving radiative transfer, and can directly output radiance results at any zenith angle.

[0155] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0156] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0157] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0158] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0159] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above descriptions are merely specific embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for determining radiance, characterized in that, The method for determining radiance includes: Obtain radiance determination parameters; Determine boundary condition information based on parameters determined by radiance; The atmospheric scattering phase matrix information is determined based on the radiance parameters. Obtain the radiative transfer model; The radiance determination parameters, the atmospheric scattering phase matrix information, and the boundary condition information are input into the radiative transfer model, and the calculation is performed according to the following formula: Where I represents radiance (I l ,I r ,I u ,I v μ represents the zenith angle cosine; Q represents the sum of single scattering and medium radiation; M represents the scattering phase matrix; c1 and c2 represent coefficients; τ represents the optical thickness; m represents the Fourier expansion order; ω(τ) represents the single scattering albedo; w represents the weight of the Gaussian point in the Gaussian integral summation framework at 2N flux numbers; l, r, u, and v represent the Stokes vector elements. The radiance is determined based on the output of the radiative transfer model.

2. The radiance determination method as described in claim 1, characterized in that, The radiance determination parameters include satellite parameters, atmospheric state parameters, lower boundary state parameters, and radiance output direction parameters.

3. The radiance determination method as described in claim 1, characterized in that, The step of determining boundary condition information based on radiance parameters includes: The thermal radiation radiance, reflection matrix, and emissivity matrix at the upper and lower boundaries of the model and the atmospheric stratification interface are calculated based on the determined parameters of radiance.

4. The radiance determination method as described in claim 1, characterized in that, The process of determining the atmospheric scattering phase matrix information based on radiance determination parameters includes: Input the radiance determination parameters into the theoretical formula for the scattering phase matrix; The atmospheric scattering phase matrix information is determined based on the output of the theoretical formula for the scattering phase matrix.

5. A radiance determination device, characterized in that, The radiance determination device includes: The parameter acquisition module is used to acquire radiance determination parameters; The atmospheric scattering phase matrix determination module is used to determine atmospheric scattering phase matrix information based on radiance determination parameters. The boundary condition determination module is used to determine boundary condition information based on the radiance determination parameters. The radiance output module is used to acquire the radiative transfer model and to input the radiance determination parameters, the atmospheric scattering phase matrix information, and the boundary condition information into the radiative transfer model for calculation according to the following formula: Where I represents radiance (I l ,I r ,I u ,I v μ represents the zenith angle cosine; Q represents the sum of single scattering and medium radiation; M represents the scattering phase matrix; c1 and c2 represent coefficients; τ represents the optical thickness; m represents the Fourier expansion order; ω(τ) represents the single scattering albedo; w represents the weight of the Gaussian point in the Gaussian integral summation framework at 2N flux numbers; l, r, u, and v represent the Stokes vector elements. The radiance output module is also used to determine the radiance based on the output results of the radiative transfer model.

6. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the radiance determination method according to any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that performs the radiance determination method according to any one of claims 1 to 4.

Images