CS G140 Graduate Computer Graphics

1. CS G140Graduate Computer Graphics Prof. Harriet Fell Spring 2009 Lecture 2 - January 14, 2009

2. Today’s Topics • Ray Tracing • Ray-Sphere Intersection • Light: Diffuse Reflection • Shadows • Phong Shading • More Math • Matrices • Transformations • Homogeneous Coordinates

3. Ray Tracinga World of Spheres

4. What is a Sphere Vector3D center; // 3 doubles double radius; double R, G, B; // for RGB colors between 0 and 1 double kd; // diffuse coeficient double ks; // specular coefficient int specExp; // specular exponent 0 if ks = 0 (double ka; // ambient light coefficient) double kgr; // global reflection coefficient double kt; // transmitting coefficient int pic; // > 0 if picture texture is used

5. -.01 .01 500 800 // transform theta phi mu distance 1 // antialias 1 // numlights 100 500 800 // Lx, Ly, Lz 9 // numspheres //cx cy cz radius R G B ka kd ks specExp kgr kt pic -100 -100 0 40 .9 0 0 .2 .9 .0 4 0 0 0 -100 0 0 40 .9 0 0 .2 .8 .1 8 .1 0 0 -100 100 0 40 .9 0 0 .2 .7 .2 12 .2 0 0 0 -100 0 40 .9 0 0 .2 .6 .3 16 .3 0 0 0 0 0 40 .9 0 0 .2 .5 .4 20 .4 0 0 0 100 0 40 .9 0 0 .2 .4 .5 24 .5 0 0 100 -100 0 40 .9 0 0 .2 .3 .6 28 .6 0 0 100 0 0 40 .9 0 0 .2 .2 .7 32 .7 0 0 100 100 0 40 .9 0 0 .2 .1 .8 36 .8 0 0

6. World of Spheres Vector3D VP; // the viewpoint int numLights; Vector3D theLights[5]; // up to 5 white lights double ka; // ambient light coefficient int numSpheres; Sphere theSpheres[20]; // 20 sphere max int ppmT[3]; // ppm texture files View sceneView; // transform data double distance; // view plane to VP bool antialias; // if true antialias

7. Simple Ray Casting for Detecting Visible Surfaces select window on viewplane and center of projection for (each scanline in image) { for (each pixel in the scanline) { determine ray from center of projection through pixel; for (each object in scene) { if (object is intersected and is closest considered thus far) record intersection and object name; } set pixel’s color to that of closest object intersected; } }

8. Ray Trace 1Finding Visible Surfaces

9. Ray-Sphere Intersection • Given • Sphere • Center (cx, cy, cz) • Radius, R • Ray from P0 to P1 • P0 = (x0, y0, z0) and P1 =(x1, y1, z1) • View Point • (Vx, Vy, Vz) • Project to window from (0,0,0) to (w,h,0)

10. y Center C = (cx, cy, cz) Radius R (x, y, z) R x (cx, cy, cz) z Sphere Equation

11. P1 t  P0 Ray Equation P0 = (x0, y0, z0) and P1 = (x1, y1, z1) The ray from P0 to P1 is given by: P(t) = (1 - t)P0 + tP1 0 <= t <= 1 = P0 + t(P1 - P0)

12. Intersection Equation P(t) = P0 + t(P1 - P0) 0 <= t <= 1 is really three equations x(t) = x0 + t(x1 - x0) y(t) = y0 + t(y1 - y0) z(t) = z0 + t(z1 - z0) 0 <= t <= 1 Substitute x(t), y(t), and z(t) for x, y, z, respectively in

13. Solving the Intersection Equation is a quadratic equation in variable t. For a fixed pixel, VP, and sphere, x0, y0, z0, x1, y1, z1, cx, cy, cz, and R eye pixel sphere are all constants. We solve for t using the quadratic formula.

14. The Quadratic Coefficients Set dx = x1 - x0 dy = y1 - y0 dz = z1 - z0 Now find the the coefficients:

15. Computing Coefficients

16. The Coefficients

17. Solving the Equation

18. The intersection nearest P0 is given by: To find the coordinates of the intersection point:

19. First Lighting Model • Ambient light is a global constant. Ambient Light = ka (AR, AG, AB) ka is in the “World of Spheres” 0  ka 1 (AR, AG, AB) = average of the light sources (AR, AG, AB) = (1, 1, 1) for white light • Color of object S = (SR, SG, SB) • Visible Color of an object S with only ambient light CS= ka (AR SR, AG SG, AB SB) • For white light CS= ka (SR, SG, SB)

20. Visible SurfacesAmbient Light

21. Second Lighting Model • Point source light L = (LR, LG, LB) at (Lx, Ly, Lz) • Ambient light is also present. • Color at point p on an object S with ambient & diffuse reflection Cp= ka (AR SR, AG SG, AB SB)+ kd kp(LR SR, LG SG, LB SB) • For white light, L = (1, 1, 1) Cp= ka (SR, SG, SB) + kd kp(SR, SG, SB) • kp depends on the point p on the object and (Lx, Ly, Lz) • kd depends on the object (sphere) • ka is global • ka + kd <= 1

22. Diffuse Light

23. Lambertian Reflection ModelDiffuse Shading • For matte (non-shiny) objects • Examples • Matte paper, newsprint • Unpolished wood • Unpolished stones • Color at a point on a matte object does not change with viewpoint.

24. Physics of Lambertian Reflection • Incoming light is partially absorbed and partially transmitted equally in all directions

25. N N  L L dA dA Geometry of Lambert’s Law 90 -   dAcos() 90 -  Surface 2 Surface 1

26. N 90 -   dAcos()  L 90 -  dA Surface 2 cos()=NL • Cp= ka (SR, SG, SB) + kd NL(SR, SG, SB)

27. Finding N N = (x-cx, y-cy, z-cz) |(x-cx, y-cy, z-cz)| normal (x, y, z) radius (cx, cy, cz)

28. Diffuse Light 2

29. Shadows on Spheres

30. More Shadows

31. Finding Shadows Shadow Ray P Pixel gets shadow color

32. Shadow Color • Given Ray from P (point on sphere S) to L (light) P= P0 = (x0, y0, z0) and L = P1 = (x1, y1, z1) • Find out whether the ray intersects any other object (sphere). • If it does, P is in shadow. • Use only ambient light for pixel.

33. Shape of Shadows

34. Different Views

35. Planets

36. Starry Skies

37. Shadows on the Plane

38. Finding Shadowson the Back Plane Shadow Ray P Pixel in Shadow

39. Close up

40. On the Table

41. Phong Highlight

42. Phong Lighting Model Light Normal Reflected View The viewer only sees the light when  is 0. N L R V    We make the highlight maximal when  is 0, but have it fade off gradually. Surface

43. Phong Lighting Model cos() = RV We use cosn(). The higher n is, the faster the drop off. N L R V For white light    Surface • Cp= ka (SR, SG, SB) + kd NL(SR, SG, SB) + ks (RV)n(1, 1, 1)

44. Powers of cos() cos10() cos20() cos40() cos80()

45. LN LN Computing R L + R=(2 LN) N R = (2 LN) N - L R L N L+R R L  

46. The Halfway Vector • From the picture • + φ = θ - φ + α So φ = α/2. H = L+ V |L+ V| UseHN instead of RV. H is less expensive to compute than R. N H This is not generally true. Why? L R φ V θ θ α Surface • Cp= ka (SR, SG, SB) + kd NL(SR, SG, SB) + ks (HN)n (1, 1, 1)

47. Varied Phong Highlights

48. Varying Reflectivity

49. Time for a Break

50. More Math • Matrices • Transformations • Homogeneous Coordinates