Global Illumination and Radiosity

Size: px

Start display at page:

Download "Global Illumination and Radiosity"

Gary Lesley King
5 years ago
Views:

1 Global Illumnaton and Radosty CS535 Danel lg. Alaga Department of Computer Scence Purdue Unversty

2 Recall: Lghtng and Shadng Lght sources Pont lght Models an omndrectonal lght source (e.g., a bulb) Drectonal lght Models an omndrectonal lght source at nfnty Spot lght Models a pont tlght wth drecton Lght model Ambent lght Dffuse reflecton Specular reflecton

3 Recall: Lghtng and Shadng Dffuse reflecton Lambertan model

4 Recall: Lghtng and Shadng Specular reflecton Phong model

$refractons,$ shadows, etc.

5 Global Illumnaton Consder drect llumnaton as well as ndrect llumnaton; e.g. Reflectons, refractons, shadows, etc. Dffuse nter reflecton wth global llumnaton only dffuse nter-reflecton drect llumnaton

$Reflectons, refractons, shadows, etc.$

6 Global Illumnaton Consder drect llumnaton as well as ndrect llumnaton; e.g. Reflectons, refractons, shadows, etc. Dffuse nter reflecton, specular nter reflecton, etc.

7 Radosty Radosty, nspred by deas from heat transfer, s an applcaton of a fnte element method to solvng the renderng equaton for scenes wth purely dffuse surfaces (renderng equaton)

8 Radosty Equaton: (more detals dtl on the board )

9 Radosty Rest of Sldes Courtesy: Dr. Maro Costa Sousa Dept. of of CS U. Of Calgary

10 Radosty Calculatng the overall lght propagaton wthn a scene, for short global llumnaton s a very dffcult problem. Wth a standard dray tracng algorthm, ths s a very tme consumng task, snce a huge number of rays have to be shot.

12 Radosty For ths reason, the radosty method was nvented. The man dea of the method s to store llumnaton values on the surfaces of the objects, as the lght s propagated startng at the lght sources.

13 Ray Tracng

14 Radosty

15 Dffuse Interreflecton (radosty method)

16 Dffuse Interreflecton Surface = "dffuse reflector" of lght energy, means: any lght energy whch strkes the surface wll be reflected n all drectons, dependent only on the angle between the surface's normal and the ncomng lght vector (Lambert's law).

17 Dffuse Interreflecton The reflected lght energy often s colored, to some small extent, by the color of the surface from whch t was reflected. Ths reflecton of lght energy n an envronment produces a phenomenon known as "color bleedng," where a brghtly colored surface's color wll "bleed" onto adjacent surfaces.

18 Dffuse Interreflecton The reflected lght energy often s colored, to some small extent, by the color of the surface from whch t was reflected. C l bl d Color bleedng, as both the red and blue walls "bleed" ther color onto the whte walls, celng and floor.

19 Radosty (Thermal Heat Transfer) The "radosty" method has ts bass n the feld of thermal heat transfer. Heat ttransfer theory descrbes radaton as the transfer of energy from a surface when that surface has been thermally excted.

20 Ths encompasses both surfaces whch are basc emtters of energy, as wth lght sources, and surfaces whch receve energy from other surfaces and thus have energy to transfer. Ths "thermal radaton" theory can be used to descrbe the transfer of many knds of energy between surfaces, ncludng lght energy.

21 Radosty (Computer Graphcs) Assumpton #1: surfaces are dffuse emtters and reflectors of energy, emttngand reflectng energy unformly over ther entre area. Assumpton #2: an equlbrum soluton can be reached; that all of the energy n an envronment s accounted for, through absorpton and reflecton. Also vewpont ndependent: the soluton wll be the same regardless of the vewpont of the mage.

22 The Radosty Equaton The "radosty equaton" descrbes the amount of energy whch hcan be emtted from a surface, as the sum of the energy nherent n the surface (a lght source, for example) and the energy whch strkes the surface, beng emtted from someother other surface. The energy whch leaves a surface (surface "j") and strkes anothers rface(surface "") satten attenuated ated by two ofactors: the "form factor" between surfaces "" and "j", whch accounts for the physcal relatonshp between the two surfaces the reflectvty of surface ", whch wll absorb a certan percentage of lght energy whch strkes the surface.

23 The Radosty Equaton B = E + ρ B F j j Radosty of surface Emssvty of surface Form Factor of surface j relatve to surface Radosty of surface j Reflectvty of surface wll absorb a certan percentage of lght energy whch strkes Surface Surface j accounts for the physcal relatonshp between the two surfaces

24 The Radosty Equaton B = E + ρ B F j j Energy reachng surface from other surfaces Surface j Surface

25 The Radosty Equaton B = E + ρ B F j j Energy reachng surface from other surfaces Form Factor of surface j relatve to surface Radosty of surface j accounts for the physcal Surface j relatonshp between the two surfaces Surface

26 The Radosty Equaton B = E + ρ B F j j Energy emtted by surface Surface j Surface

27 The Radosty Equaton B = E + ρ B F j j Energy reflected by surface Surface j Surface

28 The Radosty Equaton Energy reflected by surface B = E + ρ B F j j Form Factor of Reflectvty of surface surface j relatve to surface Energy reflected by Radosty of surface j surface = Reflectvty of surface Form * Factor Reflectvty Energy reachng accounts for surface wll absorb a from other surfaces Surface j the physcal certan relatonshp percentage of between the lght energy two Surface whch strkes surfaces

Advantages Softshadows and ndrect lghtng Vew ndependent

29 Radosty Classc radosty = fnte element e e method Assumptons Dffuse reflectance Usually polygonal surfaces Advantages Softshadows and ndrect lghtng Vew ndependent soluton Precompute for a set of lght sources Useful for walkthroughs

30 Classc Radosty Algorthm Mesh Surfaces nto Elements Compute Form Factors Between B t Elements Solve Lnear System for Radostes Reconstruct and Dsplay Soluton

31 Classc Radosty Algorthm Mesh Surfaces nto Elements Compute Form Factors Between B t Elements Solve Lnear System for Radostes Reconstruct and Dsplay Soluton

32 The Form Factor: the fracton of energy leavng one surface that reaches another surface It s a purely geometrc relatonshp, ndependent of vewpont or surface attrbutes Surface j Surface

33 Between dfferental areas, the form factor equals: dfferental l area of surface I, j FdA da j j = cos θ cos θ π r 2 j angle between Normal and r angle g between Normal j and r Surface j θ j da j θ r vector from da to da j da Surface

34 Between dfferental areas, the form factor equals: FdA da j j = cosθ cosθ π r 2 j F j The overall form factor between and j s found by ntegratng 1 cosθ cosθ = 2 A A A π r j j da da j Surface j θ j da j θ r da Surface

35 Next Step: Learn ways of computng form factors Recall the Radosty Equaton: B ρ = E + ρ B F j j The F j are the form factors Form factors ndependent of radostes (depend only on scene geometry)

36 Form Factors n (More) Detal F j cos cos = A 2 π r 1 θ θ A A j j da da j F j 1 cosθ cosθ j A π r = 2 A A j π V j da da j where V j s the vsblty (0 or 1)

37 We have two ntegrals to compute: F j = 1 cos θ cos θ j 2 A πr A A j V j da j da Area ntegral Area ntegral over surface over surface j Surface j θ j da j θ r da Surface

38 The Nusselt Analog Dfferentaton of the basc form factor equaton s dffcult even for smple surfaces! Nusselt developed a geometrc analog whch allows the smple and accurate calculaton of the form factor between a surface and a pont on a second surface.

39 The Nusselt Analog The "Nusselt analog" nvolves placng a hemsphercal projecton body, wth unt radus, at a pont on a surface. The second surface s sphercally projected onto the projecton body, then cylndrcally projected onto the base of the hemsphere. h The form factor s then the area projected on the The form factor s, then, the area projected on the base of the hemsphere dvded by the area of the base of the hemsphere.

40 Numercal Integraton: The Nusselt Analog Ths gves the form factor F daaj A j da

41 The Nusselt Analog q r 1. Project A j along ts normal: A cos q A j cos q j 2. Project result on sphere: A j cos q j / r 2 3. Project result on unt crcle: A j cos q j cos q /r 2 area A j 4. Dvde by unt crcle area: A j cos q j cos q / pr 2 q j 5. Integrate for all ponts on A j : F da A A j cos θ cos θ πr j 2 sphere projecton A j cos q j /r 2 = j V j da j unt crcle area p second projecton A j cos q j cos q /r 2

42 Method 1: Hemcube Approxmaton of Nusselt ss analog between a pont da and a polygon A j Polygonal Area (A j ) Infntesmal Area (da )

43 Hemcube For convenence, a cube 1 unt hgh wth a top face 2 x 2 s used. Sde faces are 2 wde by 1 hgh. D d l f h b Decde on a resoluton for the cube. Say 512 by 512 for the top.

44 The Hemcube In Acton

45 The Hemcube In Acton

46 The Hemcube In Acton Ths llustraton demonstrates the calculaton of form factors between a partcular surface on the wall of a room and several surfaces of objects n the room.

47 Compute the form factors from a pont on a surface to all other surfaces by: Projectng allother surfaces onto the hemcube Storng, at each dscrete area, the dentfyng ndex of the surface that s closest to the pont.

From there the form factors between the pont and the surfaces are

48 Dscrete areas wth the ndces of the surfaces whch are ultmately vsble to the pont. From there the form factors between the pont and the surfaces are calculated. For greater accuracy, a large surface would typcally be broken nto a set of small surfaces before any form factor calculaton s performed.

49 Hemcube Method 1. Scan convert all scene objects onto hemcube s 5 faces 2. Use Z buffer to determne vsblty term 3. Sum up the delta form factors of the hemcube cells covered by scanned objects 4. Gvesform factors from hemcube s base to all elements, e.e. F daaj for gven and all j

dfferental fnte form factor Alasng errors due to samplng Randomly

50 Hemcube Algorthms Advantages + Frst practcal method + Use exstng renderng systems; Hardware + Computes row of form factors n O(n) Dsadvantages Computes dfferental fnte form factor Alasng errors due to samplng Randomly rotate/shear tt/h hemcube Proxmty errors Vsblty errors Expensve to compute a sngle form factor

51 Hemcube Problem: Alasng

52 Method 2: Area Samplng 1. Subdvde de A j nto small peces da j A j 2. For all da j cast ray daj-daj to determne V j f vsble compute F dadaj cosθ cosθ j FdA da = VjdA j 2 πr sum up F daaj += F dadaj j da ray da j 3. We have now F daaj

53 Summary Several ways to fnd form factors Hemcube was orgnal method + Hardware acceleraton +GvesF daaj for all j n one pass Alasng Area samplng methods now preferred Slower than hemcube As accurate as desred snce adaptve

54 Next We have the form factors How do we fnd the radosty soluton for the scene? The "Full Matrx" Radosty Algorthm Gatherng & Shootng Progressve Radosty Meshng

55 Classc Radosty Algorthm Mesh Surfaces nto Elements Compute Form Factors Between B t Elements Solve Lnear System for Radostes Reconstruct and Dsplay Soluton

56 Recall The Radosty Equaton B = E + ρ B F j j Radosty of surface Emssvty of surface Form Factor of surface j relatve to surface Radosty of surface j Reflectvty of surface wll absorb a certan percentage of lght energy whch strkes Surface Surface j accounts for the physcal relatonshp between the two surfaces

57 Radosty Matrx E + = n j j j A B F A E A B ρ B j= j j j 1 n j j j F A F A = = + = n j j j B F E B 1 ρ n E B F F F 1 ρ ρ ρ j j j E B F B = =1 ρ = n n E E B B F F F F F F M M M O M M L L ρ ρ ρ ρ ρ ρ n n nn n n n n n E B F F F M M L M O M M ρ ρ ρ

58 Radosty Matrx The "full matrx" radosty soluton calculates the form factors between each par of surfaces n the envronment, then forms a seres of smultaneous lnear equatons. n n E E B B F F F F F F L L ρ ρ ρ ρ ρ ρ = n n nn n n n n n n E B F F F M M L M O M M ρ ρ ρ Ths matrx equaton s solved for the "B" values, whch can be used as the fnal ntensty (or color) value of each surface.

59 Radosty Matrx Ths method produces a complete soluton, at the substantal cost of frst calculatng form factors between each par of surfaces and then the soluton of the matrx equaton. Each of these steps can be qute expensve f the number of surfaces s large: complex envronments typcally have upwards of ten thousand surfaces, and envronments wth one mllon surfaces are not uncommon. Ths leads to substantal costs not only n computaton tme but n storage.

60 Next We have the form factors How do we fnd the radosty soluton for the scene? The "Full Matrx" Radosty Algorthm Gatherng & Shootng Progressve Radosty Meshng

61 Solve [F][B] = [E] Drect methods: O(n 3 ) Gaussan elmnaton Goral, Torrance, Greenberg, Battale, 1984 Iteratve t methods: O(n 2 ) Energy conservaton dagonally domnant teraton converges Gauss Sedel, Jacob: Gatherng Nshta, Nakamae, 1985 Cohen, Greenberg, 1985 Southwell: Shootng Cohen, Chen, Wallace, Greenberg, 1988

62 Gatherng Ina sense, the lght leavng patch s determned by gatherng n the lght from the rest of the envronment B B = due E to + ρ B j n j= 1 B j = ρ B F j j F j B = E + n ( ρ ) Fj j= 1 B j

63 Gatherng Gatherng lght through a hem cube allows one patch radosty to be updated. B = E + n ( ρ ) Fj j= 1 B j

64 Gatherng

65 Successve Approxmaton

66 Shootng Shootng lght through a sngle hem cube allows the whole envronment's radosty values to be updated smultaneously. For all j B j = B j + B ( ρ E ) j j whereh F = j F j A A j

67 Shootng

68 Progressve Radosty

69 Next We have the form factors How do we fnd the radosty soluton for the scene? The "Full Matrx" Radosty Algorthm Gatherng & Shootng Progressve Radosty Meshng

70 Accuracy

71 Artfacts

72 Increasng Resoluton

73 Adaptve Meshng

74 Classc Radosty Algorthm Mesh Surfaces nto Elements Compute Form Factors Between B t Elements Solve Lnear System for Radostes Reconstruct and Dsplay Soluton

75 Some Radosty Results

The Cornell Box Ths s the orgnal Cornell box, as smulated by

Greenberg for the 1984 paper Modelng the nteracton of Lght

Proceedngs), Vol. 18, No. 3, July 1984, pp. 213 222.

76 The Cornell Box Ths s the orgnal Cornell box, as smulated by Cndy M. Goral, Kenneth E. Torrance, and Donald P. Greenberg for the 1984 paper Modelng the nteracton of Lght Between Dffuse Surfaces, Computer Graphcs (SIGGRAPH '84 Proceedngs), Vol. 18, No. 3, July 1984, pp Because form factors were computed analytcally, no occludng objects were ncluded nsde the box.

77 The Cornell Box Ths smulaton of the Cornell box was done by Mchael F. Cohen and Donald P. Greenberg for the 1985 paper The Hem Cube, A Radosty Soluton for Complex Envronments, Vol. 19, No. 3, July 1985, pp The hem cube allowed form factors to be calculated usng scan converson algorthms (whch were avalable n hardware), and made t possble to calculate shadows from occludng objects.

Dscontnuty Meshng Dan Lschnsk, Flppo Tamper and Donald

for tradtonal radosty mages, many small shadow castng

82 Dscontnuty Meshng Dan Lschnsk, Flppo Tamper and Donald P. Greenberg created ths mage for the 1992 paper Dscontnuty Meshng for Accurate Radosty. It depcts a scene that represents a pathologcal case for tradtonal radosty mages, many small shadow castng detals. Notce, n partcular, the shadows cast by the wndows, and the slats n the char.

84 Opera Lghtng Ths scene from La Boheme demonstrates the use of focused lghtng and angular projecton of predstorted mages for the background. It was rendered by Jule O'B. Dorsey, Francos X. Sllon, and Donald P. Greenberg for the 1991 paper Desgn and Smulaton of Opera Lghtng and Projecton Effects.

Radosty Factory These two mages were rendered by Mchael F.

Greenberg for the 1988 paper A Progressve Refnement Approach to Fast

The factory model contans 30,000 patches, and was the most complex

86 Radosty Factory These two mages were rendered by Mchael F. Cohen, Shenchang Erc Chen, John R. Wallace and Donald P. Greenberg for the 1988 paper A Progressve Refnement Approach to Fast Radosty Image Generaton. The factory model contans 30,000 patches, and was the most complex radosty soluton computed at that tme. The radosty soluton took approxmately 5 hours for 2,000 shots, and the mage generaton requred 190 hours; each on a VAX8700.

Museum Most of the llumnaton that comes nto ths smulated museum arrves

As the progressve radosty soluton executed, users could wtness each of

of ths lght to the bottom of an adjacent baffle.

88 Museum Most of the llumnaton that comes nto ths smulated museum arrves va the baffles on the celng. As the progressve radosty soluton executed, users could wtness each of the bffl baffles b beng llumnated from above, and then reflectng some of ths lght to the bottom of an adjacent baffle. A porton of ths reflected lght was eventually bounced down nto the room. The mage appeared on the proceedngs cover of SIGGRAPH 1988.

Global Illumination and Radiosity

Global Illumination and Radiosity Global Illumnaton and Radosty CS535 Danel G. Alaga Department of Computer Scence Purdue Unversty Recall: Lghtng and Shadng Lght sources Pont lght Models an omndrectonal lght source (e.g., a bulb) Drectonal