CN112329256A - Method and device for analyzing reflection polarization characteristics of coating material - Google Patents

Method and device for analyzing reflection polarization characteristics of coating material Download PDF

Info

Publication number
CN112329256A
CN112329256A CN202011270277.8A CN202011270277A CN112329256A CN 112329256 A CN112329256 A CN 112329256A CN 202011270277 A CN202011270277 A CN 202011270277A CN 112329256 A CN112329256 A CN 112329256A
Authority
CN
China
Prior art keywords
polarization
reflection
coating material
angle
distribution function
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011270277.8A
Other languages
Chinese (zh)
Inventor
陈伟力
王淑华
陈艳
王广平
王静
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Environmental Features
Original Assignee
Beijing Institute of Environmental Features
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Environmental Features filed Critical Beijing Institute of Environmental Features
Priority to CN202011270277.8A priority Critical patent/CN112329256A/en
Publication of CN112329256A publication Critical patent/CN112329256A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention relates to a method and a device for analyzing reflection polarization characteristics of a coating material, computer equipment and a computer readable storage medium, wherein the method comprises the following steps: determining a bidirectional reflection distribution function and carrying out a polarization process according to the surface reflection characteristics of the coating material to obtain a polarization bidirectional reflection distribution function of the coating material; solving a corresponding Stokes vector expression according to the obtained polarization bidirectional reflection distribution function and the visible light reflection polarization transmission model of the coating material; and calculating the polarization degree and the polarization angle of the reflection polarization of the coating material according to the Stokes vector expression. The invention can realize the reflection polarization characteristic simulation aiming at the surface of the coating material and has high calculation precision.

Description

Method and device for analyzing reflection polarization characteristics of coating material
Technical Field
The invention relates to the technical field of target detection and identification, in particular to a method and a device for analyzing reflection polarization characteristics of a coating material, computer equipment and a computer readable storage medium.
Background
The polarization characteristics (polarization degree and polarization angle) of the target can reflect target scene information which is difficult to obtain in the intensity information, and the imaging system can be favorable for more effectively detecting and identifying the target. The polarization characteristics of visible light of the target are mainly reflection polarization characteristics, and the difference of the polarization characteristics of the target made of different materials can provide important information for target identification. For natural surfaces and artificial targets, the polarization properties of the reflected radiation depend on the intrinsic properties of the target surface, such as its structural characteristics, roughness, viewing angle, irradiance, etc. conditions. Polarization detection has been studied and used in ground-based optical observation systems for the purpose of detecting and identifying spatial objects. At present, the research on the target visible light reflection polarization characteristics is mainly based on the analysis of test data and the simulation analysis of a visible light reflection polarization transmission model, and a common material surface reflection polarization characteristic model based on a polarization bidirectional distribution reflection function is not suitable for the simulation analysis of a coating with strong mirror scattering and a coating material with atypical Gaussian distribution of scattering light intensity.
Disclosure of Invention
The invention aims to provide a method for calculating and analyzing the reflection polarization characteristic of the coating material so as to improve the accuracy of calculation and analysis of the reflection polarization characteristic of the coating material.
In order to achieve the above object, the present invention provides a method for analyzing reflection polarization characteristics of a coating material, comprising the following steps:
s1, determining the bidirectional reflection distribution function according to the surface reflection characteristics of the coating material and carrying out the polarization process to obtain the polarization bidirectional reflection distribution function of the coating material
Figure BDA0002777446230000011
The expression is as follows:
Figure BDA0002777446230000021
wherein, thetaiAt the angle of incidence zenith, θrTo reflect the zenith angle phiiIs the incident azimuth angle phirFor the azimuth of reflection, λ is the wavelength,
Figure BDA0002777446230000022
is composed of
Figure BDA0002777446230000023
Matrix element of (1), ksAnd kdRespectively representing the coefficients of the specular and diffuse components, k, of the surface of the coating materialrShowing the slope distribution of the surface of the coating material, a and b being constants related to the reflection characteristics of the surface of the coating material, G (theta)irr) Denotes a masking factor, Mj,kirirλ), j and k are 0,1,2, and 3, respectively, and represent a 4 × 4 fresnel reflection mueller matrix
Figure BDA0002777446230000024
Theta represents an included angle between the normal of the micro-surface element and the normal of the surface of the coating, and beta represents an included angle between incident light and the normal of the micro-surface element;
s2, solving a corresponding Stokes vector expression according to the obtained polarization bidirectional reflection distribution function and the visible light reflection polarization transmission model of the coating material;
and S3, calculating the polarization degree and the polarization angle of the reflection polarization of the coating material according to the Stokes vector expression.
Preferably, in step S2, when the corresponding stokes vector expression is solved, the stokes vector S of the visible light reflection polarization transmission model is solvedinExpressed as:
Figure BDA0002777446230000025
wherein omegarDenotes the reflecting solid angle, SiIndicating incident natural light, Si=[Ibg 0 0 0]T,IbgRepresenting the background radiation intensity.
Preferably, in step S3, when calculating the polarization degree and the polarization angle of the reflection polarization of the coating material, the expression is:
Figure BDA0002777446230000026
Figure BDA0002777446230000027
Figure BDA0002777446230000028
wherein I represents a total intensity image, Q represents a difference image between a light intensity image of a polarizing plate wire grid at 0 ° and a light intensity image of a polarizing plate wire grid at 90 °, U represents a difference image between a light intensity image of a polarizing plate wire grid at 45 ° and a light intensity image of a polarizing plate wire grid at 135 °, V represents a difference image between a right-handed circularly polarized image and a left-handed circularly polarized image, p represents a degree of polarization, and α represents a polarization angle.
Preferably, when the polarization bidirectional reflection distribution function of the coating material is obtained in step S1, the coefficient k of the specular component of the surface of the coating materialsCoefficient of diffusion component kdAnd the slope distribution k of the coating surfacerThe method is obtained by a coating material test measurement method, and the values of constants a and b related to the surface reflection characteristic of the coating material are obtained by measuring the surface reflection characteristic of the coating material for multiple times and performing inversion calculation.
The invention also provides a device for analyzing the reflection polarization characteristics of the coating material, which comprises:
a distribution function correction unit for determining bidirectional reflection distribution function and performing polarization process according to the surface reflection characteristics of the coating material to obtain the polarization bidirectional reflection distribution function of the coating material
Figure BDA0002777446230000031
The expression is as follows:
Figure BDA0002777446230000032
wherein, thetaiAt the angle of incidence zenith, θrTo reflect the zenith angle phiiIs the incident azimuth angle phirFor the azimuth of reflection, λ is the wavelength,
Figure BDA0002777446230000033
is composed of
Figure BDA0002777446230000034
Matrix element of (1), ksAnd kdRespectively representing the coefficients of the specular and diffuse components, k, of the surface of the coating materialrShowing the slope distribution of the surface of the coating material, a and b being constants related to the reflection characteristics of the surface of the coating material, G (theta)irr) Denotes a masking factor, Mj,kirirλ), j and k are 0,1,2, and 3, respectively, and represent a 4 × 4 fresnel reflection mueller matrix
Figure BDA0002777446230000035
Theta represents an included angle between the normal of the micro-surface element and the normal of the surface of the coating, and beta represents an included angle between incident light and the normal of the micro-surface element;
the Stokes solving unit is used for solving a corresponding Stokes vector expression according to the obtained polarization bidirectional reflection distribution function of the coating material and the visible light reflection polarization transmission model;
and the reflection polarization calculation unit is used for calculating the polarization degree and the polarization angle of the reflection polarization of the coating material according to the Stokes vector expression.
Preferably, when the stokes solving unit solves the corresponding stokes vector expression, the stokes vector S of the visible light reflection polarization transmission modelinExpressed as:
Figure BDA0002777446230000041
wherein omegarDenotes the reflecting solid angle, SiIndicating incident natural light, Si=[Ibg 0 0 0]T,IbgTo representBackground radiation intensity.
Preferably, when the reflection polarization calculation unit calculates the polarization degree and the polarization angle of the reflection polarization of the coating material, the expression is:
Figure BDA0002777446230000042
Figure BDA0002777446230000043
Figure BDA0002777446230000044
wherein I represents a total intensity image, Q represents a difference image between a light intensity image of a polarizing plate wire grid at 0 ° and a light intensity image of a polarizing plate wire grid at 90 °, U represents a difference image between a light intensity image of a polarizing plate wire grid at 45 ° and a light intensity image of a polarizing plate wire grid at 135 °, V represents a difference image between a right-handed circularly polarized image and a left-handed circularly polarized image, p represents a degree of polarization, and α represents a polarization angle.
Preferably, when the distribution function modifying unit obtains the polarization bidirectional reflection distribution function of the coating material, the mirror component coefficient k of the surface of the coating materialsCoefficient of diffusion component kdAnd the slope distribution k of the coating surfacerThe method is obtained by a coating material test measurement method, and the values of constants a and b related to the surface reflection characteristic of the coating material are obtained by measuring the surface reflection characteristic of the coating material for multiple times and performing inversion calculation.
The invention also provides computer equipment which comprises a memory and a processor, wherein the memory stores a computer program, and the processor realizes the steps of the method for analyzing the reflection polarization characteristics of the coating material when executing the computer program.
The invention also provides a computer readable storage medium, on which a computer program is stored, the computer program, when being executed by a processor, implementing the steps of the method for analyzing the reflection polarization characteristics of the coating material.
The technical scheme of the invention has the following advantages: the invention provides a method and a device for analyzing the reflection polarization characteristics of a coating material, computer equipment and a computer readable storage medium. The method can meet the requirement of simulation calculation of reflection polarization of the coating material, and has high calculation precision.
Drawings
FIG. 1 is a schematic diagram illustrating a method for analyzing reflection polarization characteristics of a coating material according to an embodiment of the present invention;
fig. 2 is a schematic structural diagram of an apparatus for analyzing reflection polarization characteristics of a coating material according to an embodiment of the present invention.
In the figure: 100: a distribution function correction unit; 200: a stokes solving unit; 300: a reflected polarization calculation unit.
Detailed Description
In order to make the objects, 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 with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
If the Stokes vector before the polarization device is recorded as SinAnd the Stokes vector after passing through the polarizer is recorded as SoutThen, for a visible light polarization imaging system, the conversion relationship between the two is as follows:
Sout=M·Sin
in the above formula, M denotes a mueller (Muller) matrix of the polarizing device. In the polarization imaging experiment, the Stokes vector obtained and calculated by the polarization imaging system is Sin. For visible light polarization imaging system, the Stokes vector S of the target material surfaceinCan be expressed as a reflected stokes vector S against the sun and sky backgroundrMainly solved by Bidirectional Reflection Distribution Function (BRDF).
At present, for a target material with a strong surface reflection, a BRDF proposed by a micro-facet theory based on specular reflection and using gaussian distribution as a probability distribution function of a micro-facet is generally selected, and an expression of the BRDF is as follows:
Figure BDA0002777446230000061
in the above formula, σ represents the roughness of the target surface, and a smaller value represents a smooth target surface; thetaiAt the angle of incidence zenith, θrTo reflect the zenith angle phiiIs the incident azimuth angle phirλ is the wavelength for the azimuth of reflection.
The polarization process of the BRDF based on the micro-planar element theory is obtained by the action of a scalar BRDF model and a 4 multiplied by 4 Fresnel reflection Mueller matrix:
Figure BDA0002777446230000062
in the above formula, the first and second carbon atoms are,
Figure BDA0002777446230000063
representing the polarization bi-directional reflection distribution function,
Figure BDA0002777446230000064
representing a fresnel reflection mueller matrix.
Based on the expression of the polarized two-way distributed reflection function and the definition of directional hemispherical reflectivity,the reflecting Stokes vector S to the sun and sky background can be calculatedrThe expression is as follows:
Figure BDA0002777446230000065
in the above formula, fi,j(i, j-0, 1,2,3) represents a polarization bilateral distribution reflection function
Figure BDA0002777446230000066
Of (2) is used.
Therefore, the Stokes vector S of the visible light reflection polarization transmission modelinCan be expressed as:
Sin=Sr=∫f(θiirr,λ)cos(θr)dΩr·Si
in the above formula, SiFor the stokes vector of the incident light,
Figure BDA0002777446230000071
in view of the fluctuation characteristics of the surface of the coating material, the invention corrects the polarization characteristic reflection model based on the micro-surface element theory by combining the advantages of the five-parameter reflection model in the aspect of the reflection characteristics of the surface of the coating material on the basis of the reflection polarization characteristic model established according to the micro-surface element theory so as to be used for calculating the reflection polarization characteristics of the surface of the coating material.
Example one
As shown in fig. 1, a method for analyzing reflection polarization characteristics of a coating material according to a first embodiment of the present invention includes the following steps:
s1, determining the bidirectional reflection distribution function according to the surface reflection characteristics of the coating material and carrying out the polarization process to obtain the polarization bidirectional reflection distribution function of the coating material
Figure BDA0002777446230000076
The expression is as follows:
Figure BDA0002777446230000072
wherein, thetaiAt the angle of incidence zenith, θrTo reflect the zenith angle phiiIs the incident azimuth angle phirFor the azimuth of reflection, λ is the wavelength,
Figure BDA0002777446230000073
is composed of
Figure BDA0002777446230000074
Wherein j and k denote matrix element numbers, ksAnd kdCoefficients representing the specular and diffuse components of the surface of the coating material, respectively, related to the roughness and reflectivity of the surface of the coating material, krThe slope distribution of the surface of the coating material is expressed and is related to the roughness and texture distribution of the surface of the material, a and b are constants related to the reflection characteristic of the surface of the coating material, and are empirical constant values corresponding to the reflection characteristic of the surface of the coating material, G (theta)irr) Denotes a masking factor, Mj,kirirλ), j and k are 0,1,2, and 3, respectively, and represent a 4 × 4 fresnel reflection mueller matrix
Figure BDA0002777446230000075
Theta represents the angle between the normal of the micro-surface element and the normal of the surface of the coating, and beta represents the angle between the incident light and the normal of the micro-surface element.
And S2, solving a corresponding Stokes vector expression according to the obtained polarization bidirectional reflection distribution function of the coating material and the visible light reflection polarization transmission model.
Preferably, in step S2, when the corresponding stokes vector expression is solved, the stokes vector S of the visible light reflection polarization transmission model is solvedinExpressed as:
Figure BDA0002777446230000081
wherein omegarDenotes the reflecting solid angle, SiIndicating incident natural light, Si=[Ibg 0 0 0]T,IbgRepresenting the background radiation intensity.
And S3, calculating the polarization degree and the polarization angle of the reflection polarization of the coating material according to the Stokes vector expression.
Preferably, in step S3, when calculating the polarization degree and the polarization angle of the reflection polarization of the coating material, the expression is:
Figure BDA0002777446230000082
Figure BDA0002777446230000083
Figure BDA0002777446230000084
wherein I represents a total intensity image, Q represents a difference image between a light intensity image of a polarizing plate wire grid at 0 ° and a light intensity image of a polarizing plate wire grid at 90 °, U represents a difference image between a light intensity image of a polarizing plate wire grid at 45 ° and a light intensity image of a polarizing plate wire grid at 135 °, V represents a difference image between a right-handed circularly polarized image and a left-handed circularly polarized image, p represents a degree of polarization, and α represents a polarization angle.
Polarization bi-directional reflection distribution function of coating material
Figure BDA0002777446230000085
And determining and solving (or simplifying) an expression based on a five-parameter reflection model to obtain a five-parameter bidirectional reflection distribution function, acting with a 4 multiplied by 4 Fresnel reflection Mueller matrix, and substituting to calculate the polarization degree and the polarization angle of the coating material. The solution of the five-parameter reflection model can be referred to the prior art, and is not described hereFurther details are provided.
In particular, when the polarized bidirectional reflectance distribution function of the coating material is obtained in step S1, the mirror component coefficient k of the surface of the coating material is obtainedsCoefficient of diffusion component kdAnd the slope distribution k of the coating surfacerThe method is obtained by a coating material test measurement method, and values of constants a and b related to the surface reflection characteristic of the coating material are obtained by measuring the surface reflection characteristic of the material for multiple times and performing inversion calculation.
Example two
As shown in fig. 2, the second embodiment provides an apparatus for analyzing the reflection polarization characteristics of a coating material, which includes a distribution function modification unit 100, a stokes solving unit 200, and a reflection polarization calculating unit 300, wherein:
the distribution function correction unit 100 is used for determining a bidirectional reflection distribution function according to the surface reflection characteristics of the coating material and performing a polarization process to obtain a polarization bidirectional reflection distribution function of the coating material
Figure BDA0002777446230000091
The expression is as follows:
Figure BDA0002777446230000092
wherein, thetaiAt the angle of incidence zenith, θrTo reflect the zenith angle phiiIs the incident azimuth angle phirFor the azimuth of reflection, λ is the wavelength,
Figure BDA0002777446230000093
is composed of
Figure BDA0002777446230000094
Matrix element of (1), ksAnd kdRespectively representing the coefficients of the specular and diffuse components, k, of the surface of the coating materialrShowing the slope distribution of the surface of the coating material, a and b being constants related to the reflection characteristics of the surface of the coating material, G (theta)irr) Denotes a masking factor, Mj,kirirλ), j and k are 0,1,2, and 3, respectively, and represent a 4 × 4 fresnel reflection mueller matrix
Figure BDA0002777446230000095
Theta represents the angle between the normal of the micro-surface element and the normal of the surface of the coating, and beta represents the angle between the incident light and the normal of the micro-surface element.
The stokes solving unit 200 is configured to solve a corresponding stokes vector expression according to the obtained polarization bidirectional reflection distribution function of the coating material and the visible light reflection polarization transmission model.
The reflective polarization calculation unit 300 is configured to calculate the polarization degree and the polarization angle of the reflective polarization of the coating material according to the stokes vector expression.
Preferably, when the stokes solving unit 200 solves the corresponding stokes vector expression, the stokes vector S of the visible light reflection polarization transmission modelinExpressed as:
Figure BDA0002777446230000101
wherein omegarDenotes the reflecting solid angle, SiIndicating incident natural light, Si=[Ibg 0 0 0]T,IbgRepresenting the background radiation intensity.
Preferably, when the reflective polarization calculation unit 300 calculates the polarization degree and the polarization angle of the reflective polarization of the coating material, the expression is:
Figure BDA0002777446230000102
Figure BDA0002777446230000103
Figure BDA0002777446230000104
wherein I represents a total intensity image, Q represents a difference image between a light intensity image of a polarizing plate wire grid at 0 ° and a light intensity image of a polarizing plate wire grid at 90 °, U represents a difference image between a light intensity image of a polarizing plate wire grid at 45 ° and a light intensity image of a polarizing plate wire grid at 135 °, V represents a difference image between a right-handed circularly polarized image and a left-handed circularly polarized image, p represents a degree of polarization, and α represents a polarization angle.
Preferably, when the distribution function modification unit 100 obtains the polarization bi-directional reflection distribution function of the coating material, the mirror component coefficient k of the surface of the coating materialsCoefficient of diffusion component kdAnd the slope distribution k of the coating surfacerThe method is obtained by a coating material test measurement method, and the values of constants a and b related to the surface reflection characteristic of the coating material are obtained by measuring the surface reflection characteristic of the coating material for multiple times and performing inversion calculation.
The content of information interaction, execution process and the like among the units of the coating material reflection polarization characteristic analysis device is based on the same concept as the method embodiment of the invention, and specific content can be referred to the description in the method embodiment of the invention, and is not described herein again.
In the above embodiments, the hardware unit may be implemented mechanically or electrically. For example, a hardware element may comprise permanently dedicated circuitry or logic (such as a dedicated processor, FPGA or ASIC) to perform the corresponding operations. The hardware elements may also comprise programmable logic or circuitry, such as a general purpose processor or other programmable processor, that may be temporarily configured by software to perform the corresponding operations. The specific implementation (mechanical, or dedicated permanent, or temporarily set) may be determined based on cost and time considerations.
In particular, in some preferred embodiments of the present invention, there is further provided a computer device, including a memory and a processor, where the memory stores a computer program, and the processor implements the steps of the method for analyzing the reflection polarization characteristics of the coating material in any one of the above embodiments when executing the computer program.
In other preferred embodiments of the present invention, there is further provided a computer-readable storage medium, on which a computer program is stored, wherein the computer program, when executed by a processor, implements the steps of the method for analyzing the reflection polarization characteristics of the coating material according to any one of the above embodiments.
It will be understood by those skilled in the art that all or part of the processes of the method for implementing the embodiments described above can be implemented by a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the method for analyzing the reflection polarization characteristics of the coating material described above, and will not be described again here.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (10)

1. A method for analyzing reflection polarization characteristics of a coating material is characterized by comprising the following steps:
s1, determining the bidirectional reflection distribution function according to the surface reflection characteristics of the coating material and carrying out the polarization process to obtain the polarization bidirectional reflection distribution function of the coating material
Figure FDA0002777446220000011
The expression is as follows:
Figure FDA0002777446220000012
wherein, thetaiAt the angle of incidence zenith, θrTo reflect the zenith angle phiiIs the incident azimuth angle phirFor the azimuth of reflection, λ is the wavelength,
Figure FDA0002777446220000013
is composed of
Figure FDA0002777446220000014
Matrix element of (1), ksAnd kdRespectively representing the coefficients of the specular and diffuse components, k, of the surface of the coating materialrShowing the slope distribution of the surface of the coating material, a and b being constants related to the reflection characteristics of the surface of the coating material, G (theta)irr) Denotes a masking factor, Mj,kirirλ), j and k are 0,1,2, and 3, respectively, and represent a 4 × 4 fresnel reflection mueller matrix
Figure FDA0002777446220000015
Theta represents an included angle between the normal of the micro-surface element and the normal of the surface of the coating, and beta represents an included angle between incident light and the normal of the micro-surface element;
s2, solving a corresponding Stokes vector expression according to the obtained polarization bidirectional reflection distribution function and the visible light reflection polarization transmission model of the coating material;
and S3, calculating the polarization degree and the polarization angle of the reflection polarization of the coating material according to the Stokes vector expression.
2. The method for analyzing the reflection polarization characteristics of the coating material according to claim 1, wherein:
in step S2, when the corresponding stokes vector expression is solved, the stokes vector S of the visible light reflection polarization transmission modelinExpressed as:
Figure FDA0002777446220000016
wherein omegarDenotes the reflecting solid angle, SiIndicating incident natural light, Si=[Ibg 0 0 0]T,IbgRepresenting the background radiation intensity.
3. The method for analyzing the reflection polarization characteristics of the coating material according to claim 2, wherein:
in step S3, when calculating the polarization degree and the polarization angle of the reflection polarization of the coating material, the expression is:
Figure FDA0002777446220000021
Figure FDA0002777446220000022
Figure FDA0002777446220000023
wherein I represents a total intensity image, Q represents a difference image between a light intensity image of a polarizing plate wire grid at 0 ° and a light intensity image of a polarizing plate wire grid at 90 °, U represents a difference image between a light intensity image of a polarizing plate wire grid at 45 ° and a light intensity image of a polarizing plate wire grid at 135 °, V represents a difference image between a right-handed circularly polarized image and a left-handed circularly polarized image, p represents a degree of polarization, and α represents a polarization angle.
4. The method for analyzing the reflection polarization characteristics of the coating material according to claim 1, wherein:
when the polarization bidirectional reflection distribution function of the coating material is obtained in the step S1, the mirror component coefficient k of the surface of the coating materialsCoefficient of diffusion component kdAnd the slope distribution k of the coating surfacerObtained by a method of testing and measuring the coating material and the surface reflection characteristic of the coating materialAnd the values of the related constants a and b are obtained by measuring the surface reflection characteristics of the coating material for multiple times and carrying out inversion calculation.
5. A device for analyzing the reflection polarization characteristics of a coating material is characterized by comprising:
a distribution function correction unit for determining bidirectional reflection distribution function and performing polarization process according to the surface reflection characteristics of the coating material to obtain the polarization bidirectional reflection distribution function of the coating material
Figure FDA0002777446220000024
The expression is as follows:
Figure FDA0002777446220000025
wherein, thetaiAt the angle of incidence zenith, θrTo reflect the zenith angle phiiIs the incident azimuth angle phirFor the azimuth of reflection, λ is the wavelength,
Figure FDA0002777446220000031
is composed of
Figure FDA0002777446220000032
Matrix element of (1), ksAnd kdRespectively representing the coefficients of the specular and diffuse components, k, of the surface of the coating materialrShowing the slope distribution of the surface of the coating material, a and b being constants related to the reflection characteristics of the surface of the coating material, G (theta)irr) Denotes a masking factor, Mj,kirirλ), j and k are 0,1,2, and 3, respectively, and represent a 4 × 4 fresnel reflection mueller matrix
Figure FDA0002777446220000033
Theta represents an included angle between the normal of the micro-surface element and the normal of the surface of the coating, and beta represents an included angle between incident light and the normal of the micro-surface element;
the Stokes solving unit is used for solving a corresponding Stokes vector expression according to the obtained polarization bidirectional reflection distribution function of the coating material and the visible light reflection polarization transmission model;
and the reflection polarization calculation unit is used for calculating the polarization degree and the polarization angle of the reflection polarization of the coating material according to the Stokes vector expression.
6. The apparatus for analyzing reflection polarization characteristics of a coating material according to claim 5, wherein:
when the Stokes solving unit solves the corresponding Stokes vector expression, the Stokes vector S of the visible light reflection polarization transmission modelinExpressed as:
Figure FDA0002777446220000034
wherein omegarDenotes the reflecting solid angle, SiIndicating incident natural light, Si=[Ibg 0 0 0]T,IbgRepresenting the background radiation intensity.
7. The apparatus for analyzing reflection polarization characteristics of a coating material according to claim 6, wherein:
when the reflection polarization calculation unit calculates the polarization degree and the polarization angle of the reflection polarization of the coating material, the expression is as follows:
Figure FDA0002777446220000035
Figure FDA0002777446220000041
Figure FDA0002777446220000042
wherein I represents a total intensity image, Q represents a difference image between a light intensity image of a polarizing plate wire grid at 0 ° and a light intensity image of a polarizing plate wire grid at 90 °, U represents a difference image between a light intensity image of a polarizing plate wire grid at 45 ° and a light intensity image of a polarizing plate wire grid at 135 °, V represents a difference image between a right-handed circularly polarized image and a left-handed circularly polarized image, p represents a degree of polarization, and α represents a polarization angle.
8. The apparatus for analyzing reflection polarization characteristics of a coating material according to claim 5, wherein:
when the distribution function correction unit obtains the polarization bidirectional reflection distribution function of the coating material, the mirror component coefficient k of the surface of the coating materialsCoefficient of diffusion component kdAnd the slope distribution k of the coating surfacerThe method is obtained by a coating material test measurement method, and the values of constants a and b related to the surface reflection characteristic of the coating material are obtained by measuring the surface reflection characteristic of the coating material for multiple times and performing inversion calculation.
9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, and wherein the processor when executing the computer program implements the steps of the method for analyzing the reflection polarization characteristics of the coating material according to any one of claims 1 to 4.
10. A computer-readable storage medium, on which a computer program is stored, wherein the computer program, when being executed by a processor, implements the steps of the method for analyzing the reflection polarization characteristics of the coating material according to any one of claims 1 to 4.
CN202011270277.8A 2020-11-13 2020-11-13 Method and device for analyzing reflection polarization characteristics of coating material Pending CN112329256A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011270277.8A CN112329256A (en) 2020-11-13 2020-11-13 Method and device for analyzing reflection polarization characteristics of coating material

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011270277.8A CN112329256A (en) 2020-11-13 2020-11-13 Method and device for analyzing reflection polarization characteristics of coating material

Publications (1)

Publication Number Publication Date
CN112329256A true CN112329256A (en) 2021-02-05

Family

ID=74317540

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011270277.8A Pending CN112329256A (en) 2020-11-13 2020-11-13 Method and device for analyzing reflection polarization characteristics of coating material

Country Status (1)

Country Link
CN (1) CN112329256A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7330376B2 (en) 2019-10-07 2023-08-21 ボストン ポーラリメトリックス,インコーポレイティド Method for augmenting sensor and imaging systems with polarized light
CN117030625A (en) * 2023-08-08 2023-11-10 中国人民解放军63921部队 Solar flare suppression method, system, equipment and medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106404675A (en) * 2015-08-03 2017-02-15 南京理工大学 Method for measuring out-of-plane polarization bidirectional reflective function of rough surface
CN107064905A (en) * 2017-05-26 2017-08-18 北京环境特性研究所 A kind of detection method of target
US20180262666A1 (en) * 2017-03-13 2018-09-13 Canon Kabushiki Kaisha Information processing apparatus, information processing method, and storage medium

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106404675A (en) * 2015-08-03 2017-02-15 南京理工大学 Method for measuring out-of-plane polarization bidirectional reflective function of rough surface
US20180262666A1 (en) * 2017-03-13 2018-09-13 Canon Kabushiki Kaisha Information processing apparatus, information processing method, and storage medium
CN107064905A (en) * 2017-05-26 2017-08-18 北京环境特性研究所 A kind of detection method of target

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
冯巍巍等: "涂层表面偏振双向反射分布函数的多参量混合模型" *
陈伟力等: "基于偏振双向分布反射函数的红外偏振特性分析" *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7330376B2 (en) 2019-10-07 2023-08-21 ボストン ポーラリメトリックス,インコーポレイティド Method for augmenting sensor and imaging systems with polarized light
CN117030625A (en) * 2023-08-08 2023-11-10 中国人民解放军63921部队 Solar flare suppression method, system, equipment and medium
CN117030625B (en) * 2023-08-08 2024-02-06 中国人民解放军63921部队 Solar flare suppression method, system, equipment and medium

Similar Documents

Publication Publication Date Title
EP3371548B1 (en) 3-d polarimetric imaging using a microfacet scattering model to compensate for structured scene reflections
CN112329256A (en) Method and device for analyzing reflection polarization characteristics of coating material
US20170268990A1 (en) Separating diffuse and specular components of a glossy object for shape reconstruction using electronic light diffusing layers (e-glass) and polarized light
US20030063793A1 (en) Broadened-specular reflection and linear subspaces for object recognition
Labunets et al. Regularized parametric model of the angular distribution of the brightness factor of a rough surface
Kallel et al. Canopy polarized BRDF simulation based on non-stationary Monte Carlo 3-D vector RT modeling
Liu et al. Polarized BRDF for coatings based on three-component assumption
Negara et al. Imaging ellipsometry for curved surfaces
Berger et al. Modeling and verifying the polarizing reflectance of real-world metallic surfaces
Schechner et al. Spaceborne underwater imaging
Vedel et al. 3D shape reconstruction of optical element using polarization
Wang et al. Effects of diffuse and specular reflections on detecting embedded defects of foams with a bifocal active imaging system at 0.22 THz
Germer et al. Mueller matrix bidirectional reflectance distribution function measurements and modeling of diffuse reflectance standards
Jarecki et al. Underdetermined polarimetric measurements for Mueller extrapolations
Kobayashi et al. Reconstructing shapes and appearances of thin film objects using RGB images
Ali et al. Toward efficient acquisition of BRDFs with fewer samples
CN110298082B (en) Method for obtaining coating bidirectional reflection distribution function simulation parameters through imaging method
Jarecki et al. Extrapolating Mueller matrices from linear Stokes images
Scarboro et al. Multistatic fiber-based system for measuring the Mueller matrix bidirectional reflectance distribution function
Gartley et al. Micro-scale surface and contaminate modeling for polarimetric signature prediction
Pasqual et al. Active polarimetry for orbital debris identification
Jafolla et al. Bidirectional reflectance measurements for high-resolution signature modeling
Jafolla et al. Theory and measurement of bidirectional reflectance for signature analysis
Weinmann et al. Appearance capture and modeling.
CN117030625B (en) Solar flare suppression method, system, equipment and medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination