For next Wednesday you will be starting to implement ray tracing (yay!).
As we discussed in class, to ray trace into a scene, you form a ray from the camera point into the scene for every pixel of the image.
In general, a ray goes from v = [v_{x},v_{y},v_{z}] into direction w = [w_{x},w_{y},w_{z}], where w is a unit length direction vector. As parameter t varies from zero to infinity, a point long the ray travels to points [ v_{x} + t w_{x}, v_{y} + t w_{y}, v_{z} + t w_{z} ]
We can place the camera eyepoint at the origin, so the first v will be [0,0,0]. Of course, when we start bouncing rays off of objects to get reflections, shadows, and so on, then v will take on all sorts of values.
You can calculate w for any pixel (i,j) by placing the image into the scene as a sort of film plane floating in front of the camera, as we discussed in class. Let's put this plane along the z axis, at z = -focalLength. Just as with your ZBuffer rendering, the smaller the value you choose for focalLength, the wider will be the view angle. You get reasonable looking pictures at focalLength values of around 3.
If your image has W columns by H rows, then you can create the ray at each pixel by:
rgb[] = new double[3]; v [] = new double[3]; w [] = new double[3]; for (int i = 0 ; i < W ; i++) // LOOP OVER IMAGE COLUMNS for (int j = 0 ; j < H ; j++) { // LOOP OVER IMAGE ROWS v[0] = 0; v[1] = 0; // CAMERA EYEPOINT IS AT THE ORIGIN v[2] = 0; w[0] = (double)(i - W/2) / (W/2); // COMPUTE RAY DIRECTION AT EACH PIXEL w[1] = (double)(H/2 - j) / (W/2); // w[2] = -focalLength; // PLACE IMAGE PLANE AT z=-focalLength normalize(w); rayTrace(v, w, rgb); // COMPUTE COLOR AT PIXEL BY RAY TRACING setPixel(i, j, rgb); }
Given a ray v,w, and a collection of spheres, each of which is described by center coordinates and radius [c_{x},c_{y},c_{z},r], we can figure out which sphere the ray hits first (if any) as follows:
(x-c_{x})^{2} + (y-c_{y})^{2} + (z-c_{z})^{2} = r^{2}by substituting (v_{x} + t w_{x}) for x, (v_{y} + t w_{y}) for y, and (v_{z} + t w_{z}) for z in the sphere equation.
When you do this you'll get some quadratic equation in t. There are two possibilities: Either (i) this equation will have no real roots (which means the ray missed the sphere entirely), or else (ii) the equation will have two roots, which means the ray has indeed hit the sphere. In this case, the value of t that you want is the smaller of these two roots, because that's where ray hits the front of the sphere.
So let's do this substitution:
( v_{x} + t w_{x} - c_{x} ) ^{2} + ( v_{y} + t w_{y} - c_{y} ) ^{2} + ( v_{z} + t w_{z} - c_{z} ) ^{2} = r^{2}Multiplying this out, we get:
(v_{x}-c_{x}) ^{2} + 2 (v_{x}-c_{x})w_{x} + t^{2} w_{x} ^{2} + (v_{y}-c_{y}) ^{2} + 2 (v_{y}-c_{y})w_{y} + t^{2} w_{y} ^{2} + (v_{z}-c_{z}) ^{2} + 2 (v_{z}-c_{z})w_{z} + t^{2} w_{z} ^{2} = r^{2}
This can be written as the quadratic equation At^{2} + Bt + C = 0, by rearranging terms:
A = w_{x} ^{2} + w_{y} ^{2} + w_{z} ^{2} ,
B = 2 (v_{x}-c_{x})w_{x} + 2 (v_{y}-c_{y})w_{y} + 2 (v_{z}-c_{z})w_{z} ,
C = (v_{x}-c_{x}) ^{2} + (v_{y}-c_{y}) ^{2} + (v_{z}-c_{z}) ^{2} - r^{2}
A more compact way of saying this using dot products is:
A = w • w B = 2 (v-c) • w C = (v-c) • (v-c) - r^{2}.
For now, we will just do phong shading on each sphere, using the sphere's surface normal. As I said in class, there are two distinct ways to think about computing the sphere's surface normal. In both cases, given that you have computed t where the ray hits the sphere, you first obtain the point where the ray hits the sphere by evaluating S = [ v_{x} + t w_{x}, v_{y} + t w_{y}, v_{z} + t w_{z} ]. Then, either:
f(S) = (x - c_{x})^{2} + (y - c_{y})^{2} + (z - c_{z})^{2} - r^{2}Taking the x,y,z partial derivatives we get:
f'(S) = [2(x - c_{x}), 2(y - c_{y}), 2(z - c_{z})].Normalize the above vector to get the surface normal at surface point S.
Assignment due Wednesdday April 9 before class
Your job for next week is to come up with an interesting collection of spheres, each of which is Phong shaded, and ray trace to them to produce a picture.
Next week we'll work on how to ray trace to other shapes, as well as more advanced topics like reflection, refraction and shadows.