620 likes | 904 Views
On robust Monte Carlo algorithms for multi-pass global illumination. Frank Suykens – De Laet 17 September 2002. Overview. Introduction Realistic image synthesis Global illumination Algorithms for global illumination Contributions Weighted multi-pass methods Path differentials
E N D
On robust Monte Carlo algorithms for multi-pass global illumination Frank Suykens – De Laet 17 September 2002
Overview • Introduction • Realistic image synthesis • Global illumination • Algorithms for global illumination • Contributions • Weighted multi-pass methods • Path differentials • Density control for photon maps • Conclusion
Overview • Introduction • Realistic image synthesis • Global illumination • Algorithms for global illumination • Contributions • Weighted multi-pass methods • Path differentials • Density control for photon maps • Conclusion
Realistic image synthesis • Goal: Compute images that appear to an observer as real photographs Which one is real?
Realistic image synthesis • Applications • Architecture • Movie industry • Lighting design • Computer games • Archeology • Product design • …
Light Transport Simulation Compute illumination Realistic image synthesis Scene description Image
Scene description • Geometry • Materials • Light sources • Camera / Eye Position, size, … (e.g., CAD)
Scene description • Geometry • Materials • Light sources • Camera / Eye Diffuse paint, glass, metal, … BSDF
Materials: BSDF • Bidirectional scattering distribution function (reflection & transmission) Fraction of incoming radiance L(x ) that is scattered into the direction θ x
BSDF Components Diffuse (D) Glossy (G) Specular (S) Diffuse, glossy and specular: (D|G|S) = X
Scene description • Geometry • Materials • Light sources • Camera / Eye Position, brightness, spotlight, …
Scene description • Geometry • Materials • Light sources • Camera / Eye Position, viewing angle, …
Light Transport Simulation • Geometry • Materials • Light sources • Camera/Eye Realistic image synthesis Compute illumination Scene description Image
Light Transport Simulation Compute illumination • For every pixel: how much light passes through? Account for all possible paths from light to eye! Global illumination
Light Transport Simulation Global illumination • Mathematical basis for light transport L x Outgoing radianceL in x in direction θ ? Rendering equation
Light Transport Simulation Lr Le L x x Self emitted radiance Reflected (& refracted) radiance x Recursive Radiance Unknown incoming radiance BSDF Integration over all directions Rendering equation = +
Light Transport Simulation Compute illumination • Geometry • Materials • Light sources • Camera/Eye Realistic image synthesis Scene description Image • Global illumination • Rendering equation
Overview • Introduction • Realistic image synthesis • Global illumination • Algorithms for global illumination • Contributions • Weighted multi-pass methods • Path differentials • Density control for photon maps • Conclusion
Example scene Many different illumination features: Indirect illumination Specular refraction Indirect caustics Caustics We want a full global illumination solution!
Algorithms for global illumination • Computation: Numerical integration • Monte Carlo integration • Algorithms • Image space algorithms • Stochastic ray tracing • Particle tracing • Bidirectional path tracing • Object space algorithms • Radiosity
Monte Carlo integration • Estimate integrals by random sampling • draw a number of random samples • average their contribution estimate of integral • Statistical errors Noise in images • Convergence: More samples, less noise
Monte Carlo integration Stochastic ray tracing • Trace paths starting from the eye E L 9 paths/pixel
Particle tracing • Trace paths starting from the light E L 9 paths/pixel Pattanaik ’92, Dutré ’93
Bidirectional path tracing • Trace paths starting from the light AND the eye E L Lafortune ’93, Veach ’94
Comparison Stochastic ray tracing (9 samples per pixel) Particle tracing (9 samples per pixel) Bidirectional path tracing (4 samples per pixel) Same computation time (± 5 min.)
Radiosity methods • Object space method • Diffuse surfaces only • View independent Galerkin radiosity
Overview • Introduction • Realistic image synthesis • Global illumination • Algorithms for global illumination • Contributions • Weighted multi-pass methods • Path differentials • Density control for photon maps • Conclusion
Overview • Introduction • Realistic image synthesis • Global illumination • Algorithms for global illumination • Contributions • Weighted multi-pass methods • Path differentials • Density control for photon maps • Conclusion
Multi-pass methods • Combine different algorithms • Separate light transport • Based on BSDF components • Different algorithms different illumination • Preserve strengths of individual algorithms • Regular expressions (e.g., LD* , LX*E ) • derive path evaluation from regular expression
E G|S D|G|S LD* Radiosity & stochastic ray tracing 1. Radiosity Use radiosity solution at end points 2. Stochastic ray tracing Full global illuminationbutdrawbacks of stoch. ray tracing LD*(G|S)X*E LX*E Combine with bidirectional path tracing
Multi-pass configuration Indirect diffuse Self-emitted light LDD+(G|S)X*E + LDD+E L(G|S)X*E LD(G|S)X*E + LDE Direct diffuse + + BPT Use weighting ??? Rad + SR
Weighted multi-pass methods • Weighting instead of separation • allow overlapping transport between different algorithms • weight individual paths automatic ‘separation’ • Technique • General Monte Carlo variance reduction technique • Constraints, weighting heuristics
Results (unweighted) Bidirectional path tracing Radiosity + stoch. ray tracing LD(G|S)X*E + LDE LD(G|S)X*E + LDE
Results (weighted) + Bidirectional path tracing Radiosity + stoch. ray tracing LD(G|S)X*E + LDE
Final result BPT only Weighted combination Radiosity + Stoch. RT and Bidirectional path tracing Radiosity + Stoch. RT
Conclusion: WMP • Multi-pass methods • separation: path evaluation from regular expression • weighting: each path is weighted individually automatic ‘separation’ • General technique • Robust combination of bidirectional path tracing and radiosity
Overview • Introduction • Realistic image synthesis • Global illumination • Algorithms for global illumination • Contributions • Weighted multi-pass methods • Path differentials • Density control for photon maps • Conclusion
Path differentials • Idea • Many algorithms trace paths • A path is infinitely thin: no neighborhood information • Knowledge about ‘region of influence’ or ‘footprint ’ would be useful in many applications: • bias-noise trade-off • Footprint definition • Path differentials
Path footprint • Path = function of random variables • direction sampling, light source sampling, …
Path footprint • Variables change path perturbation
Path footprint • Set of path perturbations footprint
Path differentials • Partial derivatives • approximate perturbations • combine into footprint (first order Taylor approx.) • footprint estimate from a single path!
Applications • Path differentials widely applicable • Any Monte Carlo path sampling algorithm • Texture filtering • Hierarchical particle tracing radiosity • Importance maps
Small elements noise Large elements blur fixed In which level should the particle contribute? Path differentials: size of footprint size of element hierarchical Application: hierarchical radiosity • Particle tracing radiosity • Trace light paths • Each hit contributes to the illumination of the element L
Application: hierarchical radiosity Fixed size (small) Fixed size (large) Path differentials
Application: hierarchical radiosity Fixed size (small) Fixed size (large) Path differentials
Conclusion: Path differentials • New, robust technique to compute path footprint • Handles general BSDFs, complex geometry • Many applications in global illumination
Overview • Introduction • Realistic image synthesis • Global illumination • Algorithms for global illumination • Contributions • Weighted multi-pass methods • Path differentials • Density control for photon maps • Conclusion
Photon mapping • Popular 2-pass global illumination algorithm • 1. Particle tracing • trace light paths Jensen ’96, …
Photon mapping • Popular 2-pass global illumination algorithm • 1. Particle tracing • trace light paths • record all hitpoints Set of photons: ‘Photon map’ Jensen ’96, …