## Radiance Caching for Efficient Global Illumination (2005)

### Cached

### Download Links

Venue: | IEEE Transactions on Visualization and Computer Graphics |

Citations: | 9 - 0 self |

### BibTeX

@INPROCEEDINGS{Krivanek05radiancecaching,

author = {Jaroslav Krivanek and Pascal Gautron and Sumanta Pattanaik and Kadi Bouatouch},

title = {Radiance Caching for Efficient Global Illumination},

booktitle = {IEEE Transactions on Visualization and Computer Graphics},

year = {2005},

pages = {550--561}

}

### OpenURL

### Abstract

In this paper we present a ray tracing based method for accelerated global illumination computation in scenes with low-frequency glossy BRDFs. The method is based on sparse sampling, caching, and interpolating radiance on glossy surfaces. In particular we extend the irradiance caching scheme proposed by Ward et al. [1] to cache and interpolate directional incoming radiance instead of irradiance. The incoming radiance at a point is represented by a vector of coefficients with respect to a spherical or hemispherical basis. The surfaces suitable for interpolation are selected automatically according to the roughness of their BRDF. We also propose a novel method for computing translational radiance gradient at a point.

### Citations

545 | Orthogonal Polynomials - Szegö - 1939 |

368 | A rapid hierarchical radiosity algorithm - Hanrahan, Salzman, et al. - 1991 |

361 |
Measuring and Modeling Anisotropic Reflection
- Ward
- 1992
(Show Context)
Citation Context ...s are multiplied by the cosine term cos i before computing the harmonics representation. ð3Þ Fig. 2. Adaptive BRDF representation for (a) Phong BRDF with exponent h 15 and (b) anisotropic Ward BRDF =-=[46]-=- with kd 0, ks 1, x 0:6, y 0:25. The order of the hemispherical harmonics representation adapts to the BRDF without ever exceeding the specified maximum representation error (here, 5 percent).... |

355 | Precomputed radiance transfer for realtime rendering in dynamic, low-frequency lighting environments
- Sloan, Kautz, et al.
(Show Context)
Citation Context ...evertheless, wavelets are a viable alternative to the use of spherical or hemispherical harmonics, especially for higher frequency BRDFs. Spherical harmonics [18], [20], [31], [32], [33], [34], [35], =-=[36]-=-, [37] remove the aliasing problem and are efficient for representing low-frequency functions. However, representation of sharp functions requires many coefficients and ringing might appear. Hemispher... |

234 | e radiance lighting simulation and rendering system - Ward - 1994 |

232 |
D.: A Ray Tracing Solution for Diffuse Interreflection
- WARD, RUBINSTEIN, et al.
- 1988
(Show Context)
Citation Context ...wfrequency glossy BRDFs. The method is based on sparse sampling, caching, and interpolating radiance on glossy surfaces. In particular, we extend the irradiance caching scheme proposed by Ward et al. =-=[1]-=- to cache and interpolate directional incoming radiance instead of irradiance. The incoming radiance at a point is represented by a vector of coefficients with respect to a hemispherical or spherical ... |

229 | Spherical wavelets: efficiently representing functions on the sphere - Schroder, Sweldens - 1995 |

183 |
Realistic Image Synthesis Using Photon Mapping
- Jensen
- 2001
(Show Context)
Citation Context .... Tawara et al. [27] selectively update a radiance sample list in time to exploit temporal coherence. Storing light particles in the scene can also be thought of as caching a directional distribution =-=[28]-=-, [29]. 2.3 Spherical Function Representation A representation of functions on a (hemi)sphere is necessary for incoming radiance caching. Piecewise constant representation [6], [24], [26], [27] is sim... |

140 |
Predicting Reflectance Functions from Complex Surfaces
- Westin, Arvo
- 1992
(Show Context)
Citation Context ...the aliasing problems. Nevertheless, wavelets are a viable alternative to the use of spherical or hemispherical harmonics, especially for higher frequency BRDFs. Spherical harmonics [18], [20], [31], =-=[32]-=-, [33], [34], [35], [36], [37] remove the aliasing problem and are efficient for representing low-frequency functions. However, representation of sharp functions requires many coefficients and ringing... |

139 |
A global illumination solution for general reflectance distributions
- Sillion, Arvo, et al.
- 1991
(Show Context)
Citation Context ...irectional incoming radiance on glossy surfaces. 2.2 Caching Directional Distributions Caching directional distributions has been used to extend the radiosity method to support glossy surfaces, e.g., =-=[18]-=-, [19], [20], [21], [22], [23], [24], [25]. It has also been used in Monte Carlo ray tracing on diffuse surfaces [6], [26], [27]. Slusallek et al. [6] and Kato [26] use reprojection of radiance sample... |

124 | Irradiance Gradients - Ward, Heckbert |

115 | Bidirectional Reflection Functions from Surface Bump Maps - Cabral, Max, et al. - 1987 |

103 | Clustered principal components for precomputed radiance transfer
- Sloan, Snyder, et al.
(Show Context)
Citation Context ...eless, wavelets are a viable alternative to the use of spherical or hemispherical harmonics, especially for higher frequency BRDFs. Spherical harmonics [18], [20], [31], [32], [33], [34], [35], [36], =-=[37]-=- remove the aliasing problem and are efficient for representing low-frequency functions. However, representation of sharp functions requires many coefficients and ringing might appear. Hemispherical h... |

102 |
Greenberg D.: Discontinuity Meshing for Accurate Radiosity
- Lischinski, Tampieri
- 1992
(Show Context)
Citation Context ...in global illumination whenever there is a certain level of smoothness in the quantity being computed. The radiosity method uses interpolation in the form of surface discretization, e.g., [10], [11], =-=[12]-=-. In the context of Monte Carlo ray tracing, approaches have been 1077-2626/05/$20.00 ß 2005 IEEE Published by the IEEE Computer Societys2 IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHICS, VOL... |

91 | Interactive Rendering using the Render Cache
- WALTER, DRETTAKIS, et al.
- 1999
(Show Context)
Citation Context ...0.00 ß 2005 IEEE Published by the IEEE Computer Societys2 IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHICS, VOL. 11, NO. 5, SEPTEMBER/OCTOBER 2005 proposed for screen space interpolation [5], =-=[13]-=-, [14], [15]. The goal of these methods is to display an approximate solution quickly. However, they do not accelerate the computation of the final high quality solution, which is the objective of our... |

91 | Frequency space environment map rendering - Ramamoorthi, Hanrahan - 2002 |

79 | hardwareaccelerated shading and lighting - Realistic - 1999 |

77 |
Arbitrary BRDF Shading for LowFrequency Lighting using Spherical Harmonics
- Kautz, Sloan, et al.
(Show Context)
Citation Context ...ems. Nevertheless, wavelets are a viable alternative to the use of spherical or hemispherical harmonics, especially for higher frequency BRDFs. Spherical harmonics [18], [20], [31], [32], [33], [34], =-=[35]-=-, [36], [37] remove the aliasing problem and are efficient for representing low-frequency functions. However, representation of sharp functions requires many coefficients and ringing might appear. Hem... |

74 |
Simulating Global Illumination Using Adaptive Meshing
- Heckbert
- 1991
(Show Context)
Citation Context ... used in global illumination whenever there is a certain level of smoothness in the quantity being computed. The radiosity method uses interpolation in the form of surface discretization, e.g., [10], =-=[11]-=-, [12]. In the context of Monte Carlo ray tracing, approaches have been 1077-2626/05/$20.00 ß 2005 IEEE Published by the IEEE Computer Societys2 IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHIC... |

69 | Bidirectional Reflection Distribution Function Expressed in Terms of Surface Scattering Modes. ECCV
- Koenderink, Doorn
- 1996
(Show Context)
Citation Context ...cients and ringing might appear. Hemispherical harmonics [9] are better suited for representing functions on a hemisphere. Basis functions similar to spherical harmonics are Zernike polynomials [38], =-=[39]-=- and the hemispherical harmonics of Makhotkin [40]. Unlike for spherical harmonics, the rotation procedure is not available for these basis functions. We choose hemispherical and spherical harmonics b... |

67 | P.: Interactive global illumination in dynamic scenes
- TOLE, PELLACINI, et al.
(Show Context)
Citation Context ..., they do not accelerate the computation of the final high quality solution, which is the objective of our work. Object space interpolation has also been used for the purpose of fast previewing [16], =-=[17]-=-. Sparse sampling and interpolation for high quality rendering was used in [7], [1]. The approach of Bala et al. [7] is suitable only for deterministic ray tracing. Ward et al. [1] use interpolation o... |

62 | Clustering for glossy global illumination
- Christensen, Lischinski, et al.
- 1997
(Show Context)
Citation Context ...sy surfaces. 2.2 Caching Directional Distributions Caching directional distributions has been used to extend the radiosity method to support glossy surfaces, e.g., [18], [19], [20], [21], [22], [23], =-=[24]-=-, [25]. It has also been used in Monte Carlo ray tracing on diffuse surfaces [6], [26], [27]. Slusallek et al. [6] and Kato [26] use reprojection of radiance samples. Tawara et al. [27] selectively up... |

59 | Combining edges and points for interactive high-quality rendering
- BALA, WALTER, et al.
- 2003
(Show Context)
Citation Context ...IEEE Published by the IEEE Computer Societys2 IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHICS, VOL. 11, NO. 5, SEPTEMBER/OCTOBER 2005 proposed for screen space interpolation [5], [13], [14], =-=[15]-=-. The goal of these methods is to display an approximate solution quickly. However, they do not accelerate the computation of the final high quality solution, which is the objective of our work. Objec... |

59 | A Hierarchical Illumination Algorithm for Surfaces with Glossy Reflection
- Aupperle, Hanrahan
- 1993
(Show Context)
Citation Context ...onal incoming radiance on glossy surfaces. 2.2 Caching Directional Distributions Caching directional distributions has been used to extend the radiosity method to support glossy surfaces, e.g., [18], =-=[19]-=-, [20], [21], [22], [23], [24], [25]. It has also been used in Monte Carlo ray tracing on diffuse surfaces [6], [26], [27]. Slusallek et al. [6] and Kato [26] use reprojection of radiance samples. Taw... |

54 | A clustering algorithm for radiance calculation in general environments
- Sillion, Drettakis, et al.
- 1995
(Show Context)
Citation Context ...ncoming radiance on glossy surfaces. 2.2 Caching Directional Distributions Caching directional distributions has been used to extend the radiosity method to support glossy surfaces, e.g., [18], [19], =-=[20]-=-, [21], [22], [23], [24], [25]. It has also been used in Monte Carlo ray tracing on diffuse surfaces [6], [26], [27]. Slusallek et al. [6] and Kato [26] use reprojection of radiance samples. Tawara et... |

54 | The irradiance Jacobian for partially occluded polyhedral sources - Arvo - 1994 |

53 | Radiance interpolants for accelerated bounded-error ray tracing
- Bala, Dorsey, et al.
- 1999
(Show Context)
Citation Context ...recise estimate of the outgoing radiance at a point. Fortunately, a high degree of coherence in the outgoing radiance field on those surfaces [1], [4], [5], [6] can be exploited by interpolation [1], =-=[7]-=- to obtain a significant performance gain. Our goal is to accelerate Monte Carlo ray tracing-based global illumination computation in the presence of surfaces with low-frequency glossy BRDFs. We achie... |

46 | Progressive radiance evaluation using directional coherence maps - Guo - 1998 |

39 | H.: Tapestry: A dynamic mesh-based display representation for interactive rendering - SIMMONS, SÉQUIN |

37 | Orthogonal Polynomials, 4th ed - Szegő - 1975 |

35 | Enhancing and Optimizing the Render Cache
- WALTER, DRETTAKIS, et al.
(Show Context)
Citation Context ... 2005 IEEE Published by the IEEE Computer Societys2 IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHICS, VOL. 11, NO. 5, SEPTEMBER/OCTOBER 2005 proposed for screen space interpolation [5], [13], =-=[14]-=-, [15]. The goal of these methods is to display an approximate solution quickly. However, they do not accelerate the computation of the final high quality solution, which is the objective of our work.... |

31 | and Wim Sweldens. Spherical wavelets: Efficiently representing functions on the sphere - Schröder - 1995 |

29 | Efficient Glossy Global Illumination with Interactive Viewing, in "Computur Graphics Forum
- STAMMINGER, SCHEEL, et al.
(Show Context)
Citation Context ...faces. 2.2 Caching Directional Distributions Caching directional distributions has been used to extend the radiosity method to support glossy surfaces, e.g., [18], [19], [20], [21], [22], [23], [24], =-=[25]-=-. It has also been used in Monte Carlo ray tracing on diffuse surfaces [6], [26], [27]. Slusallek et al. [6] and Kato [26] use reprojection of radiance samples. Tawara et al. [27] selectively update a... |

27 | Wavelet methods for radiance computations
- Schröder, Hanrahan
- 1995
(Show Context)
Citation Context ...ance on glossy surfaces. 2.2 Caching Directional Distributions Caching directional distributions has been used to extend the radiosity method to support glossy surfaces, e.g., [18], [19], [20], [21], =-=[22]-=-, [23], [24], [25]. It has also been used in Monte Carlo ray tracing on diffuse surfaces [6], [26], [27]. Slusallek et al. [6] and Kato [26] use reprojection of radiance samples. Tawara et al. [27] se... |

27 | Spherical harmonic gradients for mid-range illumination - Annen, Kautz, et al. - 2004 |

27 | hardware-accelerated shading and lighting - HEIDRICH, SEIDEL - 1999 |

26 | arbitrary BRDF shading for low-frequency lighting using spherical harmonics - Fast - 2002 |

22 | A Novel Hemispherical Basis for Accurate and Efficient Rendering - Gautron, Kˇrivánek, et al. - 2004 |

19 | Light-driven global illumination with a wavelet representation of light transport
- LEWIS, FOURNIER
- 1996
(Show Context)
Citation Context ...n glossy surfaces. 2.2 Caching Directional Distributions Caching directional distributions has been used to extend the radiosity method to support glossy surfaces, e.g., [18], [19], [20], [21], [22], =-=[23]-=-, [24], [25]. It has also been used in Monte Carlo ray tracing on diffuse surfaces [6], [26], [27]. Slusallek et al. [6] and Kato [26] use reprojection of radiance samples. Tawara et al. [27] selectiv... |

16 | A rapid hierarchical radiosity algorithm. Computer Graphics - Hanrahan, Salzman, et al. - 1991 |

14 |
Spherical Wavelets: Efficiently Representing
- Schröder, Sweldens
- 1995
(Show Context)
Citation Context ...iecewise constant representation [6], [24], [26], [27] is simple but prone to aliasing and memory demanding. Unless higher order wavelets are used, even wavelet representation [21], [22], [23], [25], =-=[30]-=- does not remove the aliasing problems. Nevertheless, wavelets are a viable alternative to the use of spherical or hemispherical harmonics, especially for higher frequency BRDFs. Spherical harmonics [... |

12 | A Signal-Processing Framework for Forward and Inverse Rendering
- Ramamoorthi
- 2002
(Show Context)
Citation Context ... problems. Nevertheless, wavelets are a viable alternative to the use of spherical or hemispherical harmonics, especially for higher frequency BRDFs. Spherical harmonics [18], [20], [31], [32], [33], =-=[34]-=-, [35], [36], [37] remove the aliasing problem and are efficient for representing low-frequency functions. However, representation of sharp functions requires many coefficients and ringing might appea... |

12 | Notes on adaptive quadrature on the hemisphere
- Shirley, Chiu
- 1994
(Show Context)
Citation Context ...ð k; kÞ; ð4Þ k1 where Lið k; kÞ is the incoming radiance coming from the sampled direction ð k; kÞ and N is the number of sampled directions. We use a fixed N, but adaptive hemisphere sampling [47], =-=[48]-=- is desirable. The order n for the incoming radiance representation is equal to the order of the BRDF representation at p. This cuts off high frequencies from the incoming radiance. The approach is ju... |

12 | A Vector Approach for Global Illumination in Ray Tracing”, Computer Graphics Forum - Zaninetti, Serpaggi, et al. - 1998 |

11 |
Basic wavefront aberration theory for optical metrology
- Wyant, Creath
(Show Context)
Citation Context ...coefficients and ringing might appear. Hemispherical harmonics [9] are better suited for representing functions on a hemisphere. Basis functions similar to spherical harmonics are Zernike polynomials =-=[38]-=-, [39] and the hemispherical harmonics of Makhotkin [40]. Unlike for spherical harmonics, the rotation procedure is not available for these basis functions. We choose hemispherical and spherical harmo... |

10 | Ch.: Direct lighting calculation by Monte Carlo integration
- Shirley, Wang
- 1991
(Show Context)
Citation Context ...he integral is split into parts and each of them is solved by a different technique: . Direct illumination uses a deterministic method for point light sources and area sampling for area light sources =-=[44]-=-. . Perfect specular reflections/refractions are solved by tracing a single deterministic secondary ray. . Ward’s irradiance caching computes the indirect diffuse term for purely diffuse surfaces. . T... |

8 | wavelet: A solution to global illumination with general surface properties - Haar - 1994 |

7 | Exploiting temporal coherence in final gathering for dynamic scenes
- TAWARA, MYSZKOWSKI, et al.
(Show Context)
Citation Context ...n used to extend the radiosity method to support glossy surfaces, e.g., [18], [19], [20], [21], [22], [23], [24], [25]. It has also been used in Monte Carlo ray tracing on diffuse surfaces [6], [26], =-=[27]-=-. Slusallek et al. [6] and Kato [26] use reprojection of radiance samples. Tawara et al. [27] selectively update a radiance sample list in time to exploit temporal coherence. Storing light particles i... |

6 |
Analysis of Radiative Transfer between Surfaces by Hemispherical Harmonics
- Lessig, Makhotkin, et al.
- 1996
(Show Context)
Citation Context ...monics [9] are better suited for representing functions on a hemisphere. Basis functions similar to spherical harmonics are Zernike polynomials [38], [39] and the hemispherical harmonics of Makhotkin =-=[40]-=-. Unlike for spherical harmonics, the rotation procedure is not available for these basis functions. We choose hemispherical and spherical harmonics because they are the only basis for which an effici... |

5 | An adaptive method for indirect illumination using light vectors
- Serpaggi, Péroche
- 2001
(Show Context)
Citation Context ...; kÞHlð k; kÞ; ð4Þ k1 where Lið k; kÞ is the incoming radiance coming from the sampled direction ð k; kÞ and N is the number of sampled directions. We use a fixed N, but adaptive hemisphere sampling =-=[47]-=-, [48] is desirable. The order n for the incoming radiance representation is equal to the order of the BRDF representation at p. This cuts off high frequencies from the incoming radiance. The approach... |