280 likes | 308 Views
Explore a novel computational method for simulating light diffusion in irregularly shaped objects with subsurface scattering. Discover how this approach extends on previous works to accurately resolve internal visibility and handle heterogeneous media such as marble. Learn about the implementation and results of this innovative technique using adaptive grid refinement and embedded boundary discretization.
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!