Method and device for modeling and simulating target visible light polarization characteristics

Abstract

The present invention provides a method and device for modeling and simulating target visible light polarization characteristics, wherein the method comprises: using a panel analysis method to divide a geometric model of a target into a plurality of panels, and determining the incident angle, detection angle and relative azimuth corresponding to each panel, and combining with an atmospheric radiation transmission model to calculate the downlink Stokes parameter of sky diffuse light irradiating the target surface, and using it as the atmospheric polarization radiation transmission characteristic, and based on a polarization characteristic calculation model of the target material surface, coupling the atmospheric polarization radiation transmission characteristic with the polarization characteristic calculation model of the target surface, thereby realizing modeling and simulation of the target visible light polarization characteristics in a field environment, and the obtained simulation model can simulate the polarization characteristics of the target under different incident angles, detection angles and atmospheric conditions, and provide technical support for the field of polarization remote sensing.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of optical modeling technology, and in particular to a method and device for modeling and simulating target visible light polarization characteristics. Background Art

[0002] When natural light is reflected by the surface of any object on Earth, it will produce certain polarization characteristics. Polarization is an inherent property of the target and is related to its own characteristics (for example, the complex refractive index of the material, surface roughness, etc.). There are obvious differences in polarization characteristics between different objects. By analyzing the polarization information on the surface of the material and inverting the intrinsic properties of the material itself, it is possible to effectively distinguish the structural characteristics between targets of different materials and improve the recognition accuracy of traditional light intensity detection.

[0003] Three-dimensional visualized polarization scene simulation analysis is the most direct and effective means to analyze the target polarization characteristics. It can simulate the target polarization characteristics and their polarization mechanisms under various boundary conditions and complex environments. It is low-cost and has good repeatability. It is easy to form the changing law of the target polarization characteristics and can be quickly applied to rapid control decisions and performance predictions of polarization imaging detection.

[0004] At present, there are few studies on modeling and simulation of visible light polarization characteristics of typical targets in field environments. There is an urgent need to provide a modeling and simulation method for target visible light polarization characteristics, which provides key and basic work for the construction of a typical target polarization characteristics database. Summary of the invention

[0005] The embodiments of the present invention provide a method and device for modeling and simulating target visible light polarization characteristics, which can realize modeling and simulation of target visible light polarization characteristics in a field environment.

[0006] In a first aspect, an embodiment of the present invention provides a method for modeling and simulating target visible light polarization characteristics, comprising:

[0007] Establishing a geometric model of the target to be simulated, and dividing the geometric model into a plurality of face elements;

[0008] Based on the detection time, target position and detector position, the incident angle, detection angle and relative azimuth corresponding to each face element in the target are determined;

[0009] According to the incident angle, detection angle and relative azimuth of each facet, the atmospheric radiation transfer model is used to calculate the downlink Stokes parameters of the sky diffuse light irradiating the target surface.

[0010] The downlink Stokes parameter is coupled to the polarization characteristic calculation model of the target material surface to obtain the total polarized radiation energy reflected by the direct sunlight and the diffuse light from the sky on each surface element in the target.

[0011] Preferably, the determining of the incident angle, detection angle and relative azimuth corresponding to each facet in the target based on the simulation time, the target position and the detector position comprises:

[0012] Determining the sun's position based on the detection time;

[0013] According to the sun position, target position and detector position, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system are calculated;

[0014] A corresponding facet coordinate system is constructed for each facet in the target, and according to the coordinate transformation relationship of the facet coordinate system described in the global coordinate system, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system are converted into the incident angle, detection angle and relative azimuth angle in each facet coordinate system.

[0015] Preferably, the solar zenith angle and solar azimuth angle in the global coordinate system are calculated using the following formula:

[0016] cos(90°-θ i )=sinφ·sinη+cosφ·cosη·cost

[0017]

[0018] Among them, θ i Characterizes the solar incidence angle, (90°-θ i ) represents the solar zenith angle, represents the solar azimuth, φ represents the latitude of the detector location, η represents the solar declination, and t represents the solar hour angle.

[0019] Preferably, the polarization characteristic calculation model is a bidirectional reflectance distribution function BRDF model.

[0020] Preferably, the atmospheric radiation transfer model is a vector radiation transfer model RT3.

[0021] Preferably, the coupling of the downlink Stokes parameter to the polarization characteristic calculation model of the target material surface is achieved by using the following formula:

[0022]

[0023] in, is the polarization characteristic calculation model, is the downlink Stokes parameter, E is the solar irradiance at the top of the atmosphere, τ i is the transmittance at the top of the atmosphere along the path from the sun to the target surface, Ω i is the incident light direction information, θ irepresents the solar incidence angle, Characterizes the solar azimuth, θ r Characterizes the solar reflection zenith angle, Characterizes the reflection azimuth.

[0024] In a second aspect, an embodiment of the present invention further provides a device for modeling and simulating target visible light polarization characteristics, comprising:

[0025] A model building unit, used to build a geometric model of the target to be simulated, and divide the geometric model into a plurality of facets;

[0026] An angle determination unit, used to determine the incident angle, detection angle and relative azimuth corresponding to each face element in the target based on the detection time, target position and detector position;

[0027] A calculation unit, for calculating the downlink Stokes parameter of the sky diffuse light irradiating the target surface according to the incident angle, detection angle and relative azimuth of each face element using the atmospheric radiation transmission model;

[0028] The coupling unit is used to couple the downlink Stokes parameter to the polarization characteristic calculation model of the target material surface to obtain the total polarized radiation energy reflected by the direct sunlight and the diffuse light from the sky on each surface element in the target.

[0029] Preferably, the angle determination unit is specifically used for:

[0030] Determining the sun's position based on the detection time;

[0031] According to the sun position, target position and detector position, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system are calculated;

[0032] A corresponding facet coordinate system is constructed for each facet in the target, and according to the coordinate transformation relationship of the facet coordinate system described in the global coordinate system, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system are converted into the incident angle, detection angle and relative azimuth angle in each facet coordinate system.

[0033] In a third aspect, an embodiment of the present invention further provides an electronic 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.

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

[0035] The embodiment of the present invention provides a method and device for modeling and simulating the polarization characteristics of visible light of a target. The geometric model of the target is divided into a plurality of facets by adopting a facet analysis method, and the incident angle, detection angle and relative azimuth corresponding to each facet are determined. The downlink Stokes parameter of the diffuse light from the sky irradiating the target surface can be calculated in combination with the atmospheric radiation transmission model, and the parameter is used as the atmospheric polarization radiation transmission characteristic. The atmospheric polarization radiation transmission characteristic is coupled with the polarization characteristic calculation model of the target material surface, thereby realizing the modeling and simulation of the polarization characteristics of the target visible light in a field environment. The obtained simulation model can simulate the polarization characteristics of the target under different incident angles, detection angles and atmospheric conditions, and provide technical support for the field of polarization remote sensing. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0037] Figure 1 It is a flow chart of a method for modeling and simulating target visible light polarization characteristics provided by one embodiment of the present invention;

[0038] Figure 2 is a hardware architecture diagram of an electronic device provided by an embodiment of the present invention;

[0039] Figure 3 It is a structural diagram of a target visible light polarization characteristic modeling and simulation device provided by one embodiment of the present invention. DETAILED DESCRIPTION

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

[0041] Please refer to Figure 1 The embodiment of the present invention provides a method for modeling and simulating target visible light polarization characteristics, the method comprising:

[0042] Step 100, establishing a geometric model of the target to be simulated, and dividing the geometric model into a plurality of facets;

[0043] Step 102, based on the detection time, the target position and the detector position, determining the incident angle, the detection angle and the relative azimuth corresponding to each facet in the target;

[0044] Step 104, according to the incident angle, detection angle and relative azimuth of each facet, the downlink Stokes parameter of the sky diffuse light irradiating the target surface is calculated using the atmospheric radiation transfer model;

[0045] Step 106, coupling the downlink Stokes parameter to the polarization characteristic calculation model of the target material surface, and obtaining the total polarized radiation energy of the direct sunlight and the diffuse light from the sky reflected by each facet surface of the target.

[0046] In an embodiment of the present invention, a geometric model of a target is divided into a number of facets by using a facet analysis method, and the incident angle, detection angle and relative azimuth corresponding to each facet are determined. The downlink Stokes parameter of diffuse light from the sky irradiating the target surface can be calculated in combination with an atmospheric radiation transmission model, and is used as the atmospheric polarization radiation transmission characteristic. The atmospheric polarization radiation transmission characteristic is coupled with a calculation model of the polarization characteristics of the target material surface, thereby achieving modeling and simulation of the target visible light polarization characteristics in a field environment. The obtained simulation model can simulate the polarization characteristics of the target under different incident angles, detection angles and atmospheric conditions, and provide technical support for the field of polarization remote sensing.

[0047] Described below Figure 1 How the various steps are performed.

[0048] First, with respect to step 100 , a geometric model of the target to be simulated is established, and the geometric model is divided into a plurality of facets.

[0049] The target in this embodiment may include space targets such as artificial satellites, spacecraft, and space stations. In order to achieve modeling and simulation of the visible light polarization characteristics of the target, it is necessary to construct a geometric model of the target, which can describe the appearance of the target. Among them, the geometric model of the target can be generated using three-dimensional modeling software.

[0050] Since the target has the appearance of a typical standard body such as a cuboid, a cylinder, a prism, a cone, etc., when detecting the target, the illumination parameters of visible light (such as sunlight) irradiated at each position on the target will be different. Therefore, in order to ensure the accuracy of modeling and simulation and improve the accuracy of modeling results, the geometric model can be divided into several facets. Among them, the smaller the facet size, the more facets there are after division, and the higher the accuracy of the modeling result. Preferably, the facet can be a micro facet.

[0051] Then, for step 102, based on the detection time, the target position and the detector position, the incident angle, the detection angle and the relative azimuth corresponding to each face element in the target are determined.

[0052] The target is in the atmospheric space. In addition to direct sunlight, there is also diffuse light from the sky. Therefore, it is necessary to comprehensively analyze the polarization characteristics of direct sunlight and diffuse light from the sky shining on the target surface.

[0053] Since the angle parameter of direct sunlight irradiating the target surface changes with the detection time, the detection time is needed to determine the position of the sun.

[0054] Specifically, determining the incident angle, detection angle, and relative azimuth corresponding to each face element in the target in step 102 may include:

[0055] S1. Determine the sun position based on the detection time;

[0056] S2. Calculate the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system according to the sun position, target position and detector position;

[0057] S3. Construct a corresponding facet coordinate system for each facet in the target, and convert the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system into the incident angle, detection angle and relative azimuth angle in each facet coordinate system according to the coordinate transformation relationship of the facet coordinate system in the global coordinate system.

[0058] The position of the sun will change at different detection times. After determining the position of the sun, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system can be calculated together with the target position and the detector position.

[0059] Among them, the solar zenith angle and solar azimuth angle in the global coordinate system can be calculated using the following formula:

[0060] cos(90°-θ i )=sinφ·sinη+cosφ·cosη·cost

[0061]

[0062] Among them, θ i Characterizes the solar incidence angle, (90°-θ i ) represents the solar zenith angle, represents the solar azimuth, φ represents the latitude of the detector location, η represents the solar declination, and t represents the solar hour angle.

[0063] According to the above formula, the changes of solar zenith angle and solar azimuth angle with detection time can be calculated.

[0064] Furthermore, since the angle parameters of visible light (direct sunlight and diffuse sky light) irradiated on each facet are different, each facet needs to be analyzed separately, so a facet coordinate system needs to be constructed for each facet. Among them, the facet coordinate system belongs to the local coordinate system. It can be understood that there is a coordinate conversion relationship between the global coordinate system and the facet coordinate system. Through this coordinate conversion relationship, the incident angle, detection angle and relative azimuth angle in each facet coordinate system can be obtained.

[0065] It should be noted that the coordinate conversion relationship between the global coordinate system and the surface source coordinate system can be determined by existing public technologies, which will not be elaborated in this embodiment.

[0066] Next, for step 104, according to the incident angle, detection angle and relative azimuth of each facet, the downlink Stokes parameter of the sky diffuse light irradiating the target surface is calculated using the atmospheric radiation transfer model.

[0067] With the deepening of research in the fields of climate change, environmental monitoring, and atmospheric remote sensing, the importance of atmospheric radiation transfer research has become increasingly prominent. It is often necessary to simulate radiation transfer, calculate parameters such as atmospheric transmittance and spectral radiance, and invert some meteorological elements. Therefore, it is urgent to develop fast, accurate, and universally applicable radiation transfer models. Common atmospheric radiation transfer models include the 6S model and the LOWTRAN model.

[0068] In one embodiment of the present invention, in order to study the actual sky polarization distribution pattern obtained during the atmospheric radiation transmission process, a simplified double-layer atmosphere model can be used to simulate the actual atmosphere, and the vector radiation transmission equation is solved based on the vector radiation transmission model (RT3) based on the doubling-accumulation method to obtain the downlink Stokes parameter of the sky diffuse light irradiating the target surface. When using the RT3 model to calculate the downlink Stokes parameter of the sky diffuse light irradiating the target surface, it can be implemented by referring to the existing calculation method, for example, it can be calculated by referring to the method recorded in the document "Atmospheric Polarization Modeling Based on Vector Radiation Transmission" (Wang Wei, Chu Jinkui, etc.).

[0069] The model parameters of the vector radiation transfer model may include: wavelength, surface type and temperature, surface albedo, radiation source, layering file, scattering file, etc. It can be seen that the modeling and simulation results are different when the atmospheric conditions are different.

[0070] Finally, for step 106, the downlink Stokes parameter is coupled to the polarization characteristic calculation model of the target material surface to obtain the total polarized radiation energy of the direct sunlight and the diffuse light from the sky reflected by each facet surface of the target.

[0071] In this embodiment, the polarization characteristic calculation model of direct sunlight irradiating the target surface can select a suitable calculation model, preferably, the polarization characteristic calculation model selects a bidirectional reflectance distribution function (BRDF) model. Model parameters may include: surface roughness, refractive index, extinction coefficient, etc.

[0072] According to the above step 102, the target is in the atmospheric space. In addition to direct sunlight, there is also diffuse light from the sky. Therefore, it is necessary to comprehensively analyze the polarization characteristics of direct sunlight and diffuse light from the sky irradiating the target surface. In the embodiment of the present invention, the downlink Stokes parameter and the polarization characteristic calculation model are coupled. Specifically, the following modeling simulation model is obtained after coupling:

[0073]

[0074] in, is the polarization characteristic calculation model, such as BRDF, is the downlink Stokes parameter, E is the solar irradiance at the top of the atmosphere, τ i is the transmittance at the top of the atmosphere along the path from the sun to the target surface, Ω i is the incident light direction information, θ i represents the solar incidence angle, Characterizes the solar azimuth, θ r Characterizes the solar reflection zenith angle, Characterize the reflection azimuth. By using the cosine of the incident angle cosθ i Corrections are made for different incident directions and the reflected sky diffuse light is integrated over the entire 2π space.

[0075] like Figure 2 , Figure 3 As shown, an embodiment of the present invention provides a device for modeling and simulating target visible light polarization characteristics. The device embodiment can be implemented by software, or by hardware or a combination of software and hardware. From the hardware level, Figure 2 As shown, it is a hardware architecture diagram of an electronic device in which a target visible light polarization characteristic modeling and simulation device provided by an embodiment of the present invention is located. Figure 2 In addition to the processor, memory, network interface, and non-volatile memory shown, the electronic device in the embodiment may also include other hardware, such as a forwarding chip responsible for processing messages, etc. Taking software implementation as an example, Figure 3As shown, as a device in a logical sense, the CPU of the electronic 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 target visible light polarization characteristic modeling and simulation device, including:

[0076] A model building unit 301 is used to build a geometric model of the target to be simulated and divide the geometric model into a plurality of facets;

[0077] An angle determination unit 302, for determining the incident angle, detection angle and relative azimuth corresponding to each facet in the target based on the detection time, target position and detector position;

[0078] The calculation unit 303 is used to calculate the downlink Stokes parameter of the sky diffuse light irradiating the target surface according to the incident angle, detection angle and relative azimuth of each face element using the atmospheric radiation transmission model;

[0079] The coupling unit 304 is used to couple the downlink Stokes parameter to the polarization characteristic calculation model of the target material surface to obtain the total polarized radiation energy of the direct sunlight and the diffuse light from the sky reflected by each surface element in the target.

[0080] In one embodiment of the present invention, the angle determination unit 302 is specifically configured to:

[0081] Determining the sun's position based on the detection time;

[0082] According to the sun position, target position and detector position, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system are calculated;

[0083] A corresponding facet coordinate system is constructed for each facet in the target, and according to the coordinate transformation relationship of the facet coordinate system described in the global coordinate system, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system are converted into the incident angle, detection angle and relative azimuth angle in each facet coordinate system.

[0084] In one embodiment of the present invention, the following formula is used to calculate the solar zenith angle and the solar azimuth angle in the global coordinate system:

[0085] cos(90°-θ i )=sinφ·sinη+cosφ·cosη·cost

[0086]

[0087] Among them, θ i Characterizes the solar incident angle, (90°-θ i ) represents the solar zenith angle, represents the solar azimuth, φ represents the latitude of the detector location, η represents the solar declination, and t represents the solar hour angle.

[0088] In one embodiment of the present invention, the polarization characteristic calculation model is a bidirectional reflectance distribution function BRDF model.

[0089] In one embodiment of the present invention, the atmospheric radiation transfer model is a vector radiation transfer model RT3.

[0090] In one embodiment of the present invention, the coupling of the downlink Stokes parameter into the polarization characteristic calculation model is achieved by using the following formula:

[0091]

[0092] in, is the polarization characteristic calculation model, is the downlink Stokes parameter, E is the solar irradiance at the top of the atmosphere, τ i is the transmittance at the top of the atmosphere along the path from the sun to the target surface, Ω i is the incident light direction information, θ i represents the solar incidence angle, Characterizes the solar azimuth, θ r Characterizes the solar reflection zenith angle, Characterizes the reflection azimuth.

[0093] It is to be understood that the structure illustrated in the embodiment of the present invention does not constitute a specific limitation on a target visible light polarization characteristic modeling and simulation device. In other embodiments of the present invention, a target visible light polarization characteristic modeling and simulation device 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.

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

[0095] An embodiment of the present invention further provides an electronic 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 modeling and simulating target visible light polarization characteristics in any embodiment of the present invention is implemented.

[0096] 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 modeling and simulating target visible light polarization characteristics in any embodiment of the present invention.

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

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

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

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

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

[0102] 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. In the absence of further restrictions, the elements defined by the statement "comprise a ..." do not exclude the presence of other identical factors in the process, method, article or device including the elements.

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

[0104] 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 modeling and simulating target visible light polarization characteristics, characterized in that: include: Establishing a geometric model of the target to be simulated, and dividing the geometric model into a plurality of face elements; Based on the detection time, target position and detector position, the incident angle, detection angle and relative azimuth corresponding to each face element in the target are determined; According to the incident angle, detection angle and relative azimuth of each facet, the atmospheric radiation transfer model is used to calculate the downlink Stokes parameters of the sky diffuse light irradiating the target surface. The downlink Stokes parameter is coupled to the polarization characteristic calculation model of the target material surface to obtain the total polarized radiation energy reflected by the direct sunlight and the diffuse light from the sky on each surface element in the target.

2. The method according to claim 1, characterized in that The method of determining the incident angle, detection angle and relative azimuth corresponding to each face element in the target based on the simulation time, the target position and the detector position includes: Determining the sun's position based on the detection time; According to the sun position, target position and detector position, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system are calculated; A corresponding facet coordinate system is constructed for each facet in the target, and according to the coordinate transformation relationship of the facet coordinate system described in the global coordinate system, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system are converted into the incident angle, detection angle and relative azimuth angle in each facet coordinate system.

3. The method according to claim 2, characterized in that The solar zenith angle and solar azimuth angle in the global coordinate system are calculated using the following formula: cos(90°-θ i )=sinφ·sin+cosφ·cosη·cost Among them, θ i Characterizes the solar incident angle, (90°-θ i ) represents the solar zenith angle, represents the solar azimuth, φ represents the latitude of the detector location, η represents the solar declination, and t represents the solar hour angle.

4. The method according to claim 1, characterized in that: The polarization characteristic calculation model is a bidirectional reflectance distribution function BRDF model.

5. The method according to claim 1, characterized in that: The atmospheric radiation transfer model is a vector radiation transfer model RT3.

6. The method according to claim 1, characterized in that The following formula is used to couple the downlink Stokes parameter to the polarization characteristic calculation model of the target material surface: in, is the polarization characteristic calculation model, is the downlink Stokes parameter, E is the solar irradiance at the top of the atmosphere, τ i is the transmittance at the top of the atmosphere along the path from the sun to the target surface, Ω i is the incident light direction information, θ i represents the solar incidence angle, Characterizes the solar azimuth, θ r Characterizes the solar reflection zenith angle, Characterizes the reflection azimuth.

7. A target visible light polarization characteristic modeling and simulation device, characterized in that: include: A model building unit, used to build a geometric model of the target to be simulated, and divide the geometric model into a plurality of facets; An angle determination unit, used to determine the incident angle, detection angle and relative azimuth corresponding to each face element in the target based on the detection time, target position and detector position; A calculation unit, for calculating the downlink Stokes parameter of the sky diffuse light irradiating the target surface according to the incident angle, detection angle and relative azimuth of each face element using the atmospheric radiation transmission model; The coupling unit is used to couple the downlink Stokes parameter to the polarization characteristic calculation model of the target material surface to obtain the total polarized radiation energy reflected by the direct sunlight and the diffuse light from the sky on each surface element in the target.

8. The device according to claim 7, characterized in that The angle determination unit is specifically used for: Determining the sun's position based on the detection time; According to the sun position, target position and detector position, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system are calculated; A corresponding facet coordinate system is constructed for each facet in the target, and according to the coordinate transformation relationship of the facet coordinate system described in the global coordinate system, the solar zenith angle, solar azimuth angle, detector zenith angle and detector azimuth angle in the global coordinate system are converted into the incident angle, detection angle and relative azimuth angle in each facet coordinate system.

9. An electronic 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 7 is implemented.

10. 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 7.

Images