280 likes | 422 Vues
This study presents a novel computational approach to simulate subsurface scattering in non-metallic objects, such as wax, skin, marble, and fruits. Building upon traditional models like Monte Carlo methods and analytical dipole models, our method enables accurate light diffusion simulations in arbitrarily shaped objects. We introduce adaptive grid refinement and embedded boundary discretization techniques, which significantly enhance accuracy while resolving internal visibility in heterogeneous media. Our implementation provides practical solutions for real-time applications, making it a vital advancement in computer graphics research.
E N D
A Computational Approach to Simulate Light Diffusion in Arbitrarily Shaped Objects Tom Haber, Tom Mertens, Philippe Bekaert, Frank Van Reeth University of Hasselt Belgium
Subsurface Scattering • All non-metallic objects • Examples: wax, skin, marble, fruits, ... Traditional Reflection Model Subsurface scattering Images courtesy of Jensen et al. 2001
Previous Work • Monte-Carlo volume light transport • Accurate, but slow for highly-scattering media • Analytical dipole model [Jensen01] • Inaccurate (semi-infinite plane, no internal visibility) • Fast (basis for interactive methods) • Inherently limited to homogeneous media • Multigrid [Stam95] • Simple Finite Differencing • Only illustrative examples in 2D • Our method extends on this work
Goals • Simulate subsurface scattering • Accurate for arbitrarily shaped objects • Capable of resolving internal visibility • Heterogeneous media • Varying material coefficients • E.g. Marble • Only highly scattering media
Diffusion Equation Diffusion Equation Stopping term Source term Diffusion term Boundary Conditions
1 1 -4 1 1 Overview Finite-Differencing (FD) • Large amount of memory in 3D • Badly approximates the surface • Impractical!
Adaptive Grid Refinement Embedded Boundary Discretization • FD but… • 1th order surface approximation • Allows coarser grid • O(h2) accurate everywhere! • Badly approximates high curvature regions • Still requires quite some memory
FD vs. EBD • FD yields instabilities near the boundary • EBD results in a consistent solution FD EBD
Implementation • Preprocessing (prep) • Construction of volumetric grid • Adaptive mesh refinement • Source term computation (src) • Visibility tests to light sources • Attenuation • Solve using multigrid • Visualization Implemented on a pentium 4 1.7 Ghz with 512 MB RAM
Monte-Carlo Comparison Jensen et al. Our method Monte-Carlo
Monte-Carlo Comparison Jensen et al. Our method Monte-Carlo
Monte-Carlo Comparison Jensen et al. Our method Monte-Carlo
Chromatic bias in source • Highly exponential falloff for opaque objects • Requires small cells • Workaround: use irradiance at the surface as source Average color Distance (mm)
Conclusion • Contributions • Multigrid made practical in 3D • Embedded boundary discretization • Adaptive Grid Refinement • Heterogeneous materials • Limitations • Grid size • Assumptions of the diffusion eq. • Future Work • More efficient subdivision scheme • Perceptual metrics
Thank you! • Acknowledgements • tUL impulsfinanciering • Interdisciplinair instituut voor Breed-BandTechnologie
Adaptive Mesh Refinement • Three-point interpolation scheme • Implies several constraints • Neighboring cells cannot differ by more than one level • Cells neighboring a cut-cell must all be on the same level
Overview • Outline • Construct volumetric grid • Discretize diffusion eq. • Solve using multigrid • Finite-Differencing (FD)
1 1 -4 1 1 Overview • Outline • Construct volumetric grid • Discretize diffusion eq. • Solve using multigrid • Finite-Differencing (FD)
1 1 -4 1 1 Overview • Outline • Construct volumetric grid • Discretize diffusion eq. • Solve using multigrid • Finite-Differencing (FD) • Requires large amount of memory in 3D • Badly approximates the surface • Impractical!