Higher Order Ray Marching

Rendering participating media is still a challenging and time consuming task. In such media light interacts at every differential point of its path. Several rendering algorithms are based on ray marching: dividing the path of light into segments and calculating interactions at each of them. In this work, we revisit and analyze ray marching both as a quadrature integrator and as an initial value problem solver, and apply higher order adaptive solvers that ensure several interesting properties, such as faster convergence, adaptiveness to the mathematical definition of light transport and robustness to singularities. We compare several numerical methods, including standard ray marching and Monte Carlo integration, and illustrate the benefits of different solvers for a variety of scenes. Any participating media rendering algorithm that is based on ray marching may benefit from the application of our approach by reducing the number of needed samples (and therefore, rendering time) and increasing accuracy.

[1]  J. Dormand,et al.  A family of embedded Runge-Kutta formulae , 1980 .

[2]  Wojciech Jarosz,et al.  Gaussian quadrature for photon beams in Tangled , 2011, SIGGRAPH '11.

[3]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[4]  J. M. Watt.  Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

[5]  Philipp Slusallek,et al.  A mathematical framework for efficient closed-form single scattering , 2011, Graphics Interface.

[6]  Derek Nowrouzezahrai,et al.  A comprehensive theory of volumetric radiance estimation using photon points and beams , 2011, TOGS.

[7]  Kenneth E. Torrance,et al.  The zonal method for calculating light intensities in the presence of a participating medium , 1987, SIGGRAPH.

[8]  W. A. Coleman.  Mathematical Verification of a Certain Monte Carlo Sampling Technique and Applications of the Technique to Radiation Transport Problems , 1968 .

[9]  Matthias Zwicker,et al.  The Beam Radiance Estimate for Volumetric Photon Mapping , 2008, SIGGRAPH '08.

[10]  James T. Kajiya,et al.  Ray tracing volume densities , 1984, SIGGRAPH.

[11]  W. H. Enright,et al.  Effective solution of discontinuous IVPs using a Runge-Kutta formula pair with interpolants , 1988 .

[12]  Alexander Keller,et al.  Unbiased Global Illumination with Participating Media , 2008 .

[13]  Derek Nowrouzezahrai,et al.  Virtual ray lights for rendering scenes with participating media , 2012, ACM Trans. Graph..

[14]  L. Shampine,et al.  A 3(2) pair of Runge - Kutta formulas , 1989 .

[15]  Christopher D. Kulla,et al.  Eurographics Symposium on Rendering 2012 Importance Sampling Techniques for Path Tracing in Participating Media , 2022 .

[16]  G. Rybicki.  Radiative transfer , 2019, Climate Change and Terrestrial Ecosystem Modeling.

[17]  Adolfo Muñoz,et al.  Differential Ray Marching , 2012, CEIG.

[18]  Leonidas J. Guibas,et al.  Robust Monte Carlo methods for light transport simulation , 1997 .

[19]  W. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[20]  O. Heinisch,et al.  Winchester, A. M.: Heredity and your life. Dover Publications Inc., New York, 1960. 333 Seiten. Preis $ 1,45 , 1964 .

[21]  Francisco J. Serón,et al.  Simulation of atmospheric phenomena , 2006, Comput. Graph..

[22]  Derek Nowrouzezahrai,et al.  Joint importance sampling of low-order volumetric scattering , 2013, ACM Trans. Graph..

[23]  Sumanta N. Pattanaik,et al.  Computation of global illumination in a participating medium by monte carlo simulation , 1993, Comput. Animat. Virtual Worlds.

[24]  H. Jensen.  Realistic Image Synthesis Using Photon Mapping , 2001 .

[25]  T. Patterson,et al.  The optimum addition of points to quadrature formulae. , 1968 .

[26]  Derek Nowrouzezahrai,et al.  Progressive Virtual Beam Lights , 2012, Comput. Graph. Forum.