Metropolis light transport

We present a new Monte Carlo method for solving the light transport problem, inspired by the Metropolis sampling method in computational physics. To render an image, we generate a sequence of light transport paths by randomly mutating a single current path (e.g. adding a new vertex to the path). Each mutation is accepted or rejected with a carefully chosen probability, to ensure that paths are sampled according to the contribution they make to the ideal image. We then estimate this image by sampling many paths, and recording their locations on the image plane. Our algorithm is unbiased, handles general geometric and scattering models, uses little storage, and can be orders of magnitude more efficient than previous unbiased approaches. It performs especially well on problems that are usually considered difficult, e.g. those involving bright indirect light, small geometric holes, or glossy surfaces. Furthermore, it is competitive with previous unbiased algorithms even for relatively simple scenes. The key advantage of the Metropolis approach is that the path space is explored locally, by favoring mutations that make small changes to the current path. This has several consequences. First, the average cost per sample is small (typically only one or two rays). Second, once an important path is found, the nearby paths are explored as well, thus amortizing the expense of finding such paths over many samples. Third, the mutation set is easily extended. By constructing mutations that preserve certain properties of the path (e.g. which light source is used) while changing others, we can exploit various kinds of coherence in the scene. It is often possible to handle difficult lighting problems efficiently by designing a specialized mutation in this way. CR Categories: I.3.7 [Computer Graphics]: ThreeDimensional Graphics and Realism; I.3.3 [Computer Graphics]: Picture/Image Generation; G.1.9 [Numerical Analysis]: Integral Equations—Fredholm equations.

[1]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[2]  P. Peskun,et al.  Optimum Monte-Carlo sampling using Markov chains , 1973 .

[3]  Turner Whitted,et al.  An improved illumination model for shaded display , 1979, CACM.

[4]  Robert L. Cook,et al.  Distributed ray tracing , 1984, SIGGRAPH.

[5]  James T. Kajiya,et al.  The rendering equation , 1986, SIGGRAPH.

[6]  D. C. Rapaport,et al.  Book review:Monte Carlo methods. Volume I: Basics , 1987 .

[7]  David E. Goldberg,et al.  Genetic Algorithms in Search Optimization and Machine Learning , 1988 .

[8]  Paul S. Heckbert Adaptive radiosity textures for bidirectional ray tracing , 1990, SIGGRAPH.

[9]  James Arvo,et al.  Particle transport and image synthesis , 1990, SIGGRAPH.

[10]  Pat Hanrahan,et al.  Illumination from curved reflectors , 1992, SIGGRAPH.

[11]  Yves D. Willems,et al.  Bi-directional path tracing , 1993 .

[12]  Yves D. Willems,et al.  A Theoretical Framework for Physically Based Rendering , 1994, Comput. Graph. Forum.

[13]  J. Doug Tygar,et al.  Computability and complexity of ray tracing , 1994, Discret. Comput. Geom..

[14]  Yves D. Willems,et al.  A 5D Tree to Reduce the Variance of Monte Carlo Ray Tracing , 1995, Rendering Techniques.

[15]  Leonidas J. Guibas,et al.  Bidirectional Estimators for Light Transport , 1995 .

[16]  Yves D. Willems,et al.  Potential-driven Monte Carlo Particle Tracing for Diffuse Environments with Adaptive Probability Functions , 1995, Rendering Techniques.

[17]  Donald P. Greenberg,et al.  Global Illumination via Density Estimation , 1995, Rendering Techniques.

[18]  Henrik Wann Jensen,et al.  Importance Driven Path Tracing using the Photon Map , 1995, Rendering Techniques.

[19]  Steven Collins,et al.  Reconstruction of Illumination from Area Luminaires , 1995, Rendering Techniques.

[20]  Leonidas J. Guibas,et al.  Optimally combining sampling techniques for Monte Carlo rendering , 1995, SIGGRAPH.

[21]  Brigitta Lange,et al.  Rayvolution: An Evolutionary Ray Tracing Algorithm , 1995 .

[22]  Sumanta N. Pattanaik,et al.  Adjoint equations and random walks for illumination computation , 1995, TOGS.

[23]  Peter Shirley,et al.  Monte Carlo techniques for direct lighting calculations , 1996, TOGS.

[24]  Eric Veach,et al.  Non-symmetric Scattering in Light Transport Algorithms , 1996, Rendering Techniques.

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