1 Rendering Equation Outline Outline Surfaces (interreflection) x - PDF document

Radiosity Radiosity Foundations of Computer Graphics Foundations of Computer Graphics Cornell box with color bleeding [Goral et al 84] (Spring 2010) (Spring 2010) CS 184, Lecture 21: Radiosity http://inst.eecs.berkeley.edu/~cs184 Advantages and Disadvantages Advantages and Disadvantages Photograph of a sculpture.  Radiosity methods track rate at which energy (radiosity) The front faces are all diffuse white leaves [diffuse] surfaces The color is because of reflection from rear-facing colored faces  Determine equilibrium of light energy in a view- Raytracing makes all faces white. independent way It can handle specular reflection and shadows, but not diffuse-diffuse interreflection or color bleeding  Allows for diffuse interreflection, color bleeding, and walkthroughs Radiosity correctly captures the  Difficult to handle specular objects, mirrors color bleeding from the back of the boards to the front. General Approach General Approach Earliest Radiosity Earliest Radiosity pictures pictures Radiosity was first developed in other fields  Assume diffuse surfaces discretized into a finite set  Heat transport, Lighting Design of patches or finite elements  In graphics: Goral et al. 84  Radiosity equation is a matrix equation or set of simultaneous linear equations derived by approximations to the rendering equation  Solve iteratively using numerical methods Parry Moon and Domina Spencer (MIT), Lighting Design, 1948 1

Rendering Equation Outline Outline Surfaces (interreflection) x   Rendering equation review dA  Radiosity equation  Form factors  d   i i r  Methods to compute form factors x    x x i             ( , ) ( , ) ( , ) ( , , ) c os L x L x L x f x d r r e r r i i r i i  High-level overview only. Best textual reference is probably Reflected Light Emission BRDF Cosine of Reflected Incident angle Sections 16.3.1 and 16.3.2 in FvDFH. This will be handed out. (Output Image) Light If curious, read the rest of 16.3 and parts of Cohen and Wallace. UNKNOWN KNOWN UNKNOWN KNOWN KNOWN Change of Variables Change of Variables                         ( , ) ( , ) ( , ) ( , , ) cos ( , ) ( , ) ( , ) ( , , ) cos L x L x L x f x d L x L x L x f x d r r e r r i i r i i r r e r r i i r i i   Integral over angles sometimes insufficient. Write integral Integral over angles sometimes insufficient. Write integral in terms of surface radiance only (change of variables) in terms of surface radiance only (change of variables)   cos cos            ( , ) ( , ) ( , ) ( , , ) i o x  L x L x L x f x d A  r r e r r i i r |  | 2 x x  all visible to x x  dA  o   i     d  cos cos dA dA       i d o d o i     i | | 2 i | | 2 i x x x x   cos cos x     ( , ) ( , ) i o G x x G x x   2 | | x x Rendering Equation: Standard Form Radiosity Equation             ( , ) ( , ) ( , ) ( , , ) cos L x L x L x f x d              ( , ) ( , ) ( , ) ( , , ) ( , ) ( , ) L x L x L x f x G x x V x x dA r r e r r i i r i i r r e r r i i r  all surfaces  x Integral over angles sometimes insufficient. Write integral Drop angular dependence (diffuse Lambertian surfaces) in terms of surface radiance only (change of variables)   cos cos         ( )  ( )  ( ) ( ) ( , ) ( , ) ( ,  )  ( ,  )  ( ,   ) ( ,   , ) i o L x L x f x L x G x x V x x dA L x L x L x f x d A  r r e r r i i r  2 r e r | | x x all  visible to x x S Change variables to radiosity (B) and albedo ( ρ ) Domain integral awkward. Introduce binary visibility fn V                ( , ) ( , ) ( , ) ( , ) ( , ) ( , , ) ( , ) ( , ) G x x V x x L x L x L x f x G x x V x x dA       ( ) ( ) ( ) ( ) r r e r r i i r B x E x x B x dA   all surfaces x   cos dA S   o d Same as equation 2.52 Cohen Wallace. It swaps primed Expresses conservation of light energy at all points in space  i  2 | | x x And unprimed, omits angular args of BRDF, - sign.   cos cos   Same as equation 2.54 in Cohen Wallace handout (read sec 2.6.3)   ( , ) ( , ) G x x G x x i o   2 | | Ignore factors of π which can be absorbed. x x 2

Outline Discretization and Form Factors and Form Factors Outline Discretization  Rendering equation review   ( , ) ( , ) G x x V x x       ( ) ( ) ( ) ( ) B x E x x B x dA  S  Radiosity equation   A   j B E B F  i i i j j i A j i  Form factors F is the form factor. It is dimensionless and is the fraction of energy leaving the entirety of patch j ( multiply by area of j to  Methods to compute form factors get total energy) that arrives anywhere in the entirety of patch i ( divide by area of i to get energy per unit area or radiosity). Section 16.3.1,2 (eqs 16.63-65) in FvDFH Form Factors Matrix Equation Form Factors Matrix Equation   A   j B E B F   i i i j j i dA A j j j i A   j   ( , ) ( , ) G x x V x x  r A F A F dAdA    i i j j j i i j    i   B E B F  i i i j i j A dA  i j i    B B F E  i i j i j i     ( , ) ( , ) G x x V x x j  A F A F dAdA     i i j j j i i j      M B E MB E M I F   cos cos  ij j i ij ij i i j     ( , ) ( , ) i o G x x G x x j   | | 2 x x Outline Outline Nusselt’ Nusselt ’s s Analog Analog  Rendering equation review Analytically project into hemisphere above point. Then project onto hemisphere base  Radiosity equation Form factor is ratio of area on base to  Form factors area of entire base This computes differential point to patch form factor  Methods to compute form factors Section 16.3.2 in FvDFH Why does it work? 3

Hemicube Hemicubes Hemicubes Hemicube   cos cos    i p F A p  2 r  Each small hemicube cell has a precomputed delta form factor: add up to get final value  We can render the scene using normal Z-buffer scan conversion onto the faces of the hemicube! Monte Carlo Ray Tracing Monte Carlo Ray Tracing  Can be used to find form factors (slow)  Can be used directly to shoot energy 4