450 likes | 702 Views
Computation of Polarized Subsurface BRDF for Rendering. Charly Collin – Sumanta Pattanaik – Patrick LiKamWa Kadi Bouatouch. Painted materials. Painted materials. Painted materials. Painted materials. Our goal. Compute the subsurface BRDF from physical properties:. Base layer
E N D
Computation of Polarized Subsurface BRDF for Rendering Charly Collin – SumantaPattanaik – Patrick LiKamWa Kadi Bouatouch
Our goal Compute the subsurface BRDF from physical properties: • Base layer • Binder thickness • Particle properties: • Refractive indices • Particle radius • Particle distribution
Our goals Compute the diffuse BRDF from physical properties: • Base layer • Binder thickness • Particle properties: • Refractive indices • Particle radius • Particle distribution Use polarization in our computations: • Accurate light transport simulation: • Accurate BRDF computation • Accurate global illumination
Polarization • Light is composed of waves • Unpolarized light is composed of waves with random oscillation • Light is polarized when composed of waves sharing similar oscillation • Polarization of the light can be: • Linear • Circular • Both • Polarization properties change the way light interacts with matter
Polarization The Stokes vectoris a usefulrepresentation for polarized light
Polarization • Each light-matter interaction changes the radiance, but also the polarization state of the light • Modifications to a Stokes vector are donethrough a 4x4 matrix, the Mueller matrix: = • Polarized BRDF, or polarized phase function are represented as Mueller matrices
BRDF Computation To compute the BRDF weneed to compute the radiance field for: • Each incident and outgoing direction • 4 linearlyindependent incident Stokes vectors ? ? ? The radiance fieldiscomputed by solving light transport
BRDF Computation Light transport ismodeledthrough the Vector Radiative Transfer Equation: ? ? ?
BRDF Computation Our computation makes several assumptions on the material: • Plane parallel medium
BRDF Computation Our computation makes several assumptions on the material: • Plane parallel medium • Randomly oriented particles
BRDF Computation Our computation makes several assumptions on the material: • Plane parallel medium • Randomly oriented particles • Homogeneous layers
Vector Radiative Transfer Equation It has 3 components: • the radiance • corresponding to the light scattering inside the material RTE expresses the change of radiance along optical depth .
Vector Radiative Transfer Equation It has 3 components: • the radiance • corresponding to the light scattering inside the material • accounting for attenuated incident radiance RTE expresses the change of radiance along optical depth .
VRTE Solution • VRTE is solved using Discrete Ordinate Method (DOM) • Solution is composed of an homogeneous and 4N particular solution • The homogeneous solution consists of a 4Nx4N Eigen problem • Each particular solutionconsists of two set of 4N linearequations to solve +
Results: Polarization Subsurface BRDF exhibitspolarizationeffects
Results: Different materials Titaniumdioxide Aluminium arsenide Ironoxide Gold
Results: Different materials – BRDF lobe Titaniumdioxide Alluminium arsenide Ironoxide Gold
Results: Different materials – Degree of polarization Titaniumdioxide Alluminium arsenide Ironoxide Gold
Results : Different materials – Diffuse base (BRDF) Titaniumdioxide Aluminium arsenide Ironoxide Gold
Results: Different materials – Diffuse base (DOP) Titaniumdioxide Aluminium arsenide Ironoxide Gold
Results: Different materials – Metallic base (BRDF) Titaniumdioxide Aluminium arsenide Ironoxide Gold
Results: Different materials – Metallic base (DOP) Titaniumdioxide Aluminium arsenide Ironoxide Gold
Results: Accuracy – Benchmark validation Benchmark data fromWauben and Hovenier (1992)
Results: Accuracy Takingpolarizationintoaccountsyieldsbetterprecision
Demo • BRDF Solver • Polarizedrenderer
VRTE Solution Use of the DiscreteOrdinateMethod (DOM):
VRTE Solution The VRTE can be written as: That we reorganize: Components expressed using Components independant of
VRTE Solution We introduce an differential operator : Needs to be solved for each and
VRTE Solution Standard solution is the combination of the homogeneous solution... ... and one particular solution. +
VRTE Solution • The homogeneous solution consists of an 4N x 4N Eigen problem • The particular solutionconsists of a set of 4N linearequations to solve • It needs to besolved for each