Follow-up to Raytracing: A (Very) Brief Introduction from Nick's blog
Last time I mentioned that rays are cast from the eye point. This isn’t much good of course unless you can tell what the ray collides with. So here’s an idea of how to check for the intersection of a ray and a flat plane.
Representing our objects
Here’s the mathsy bit. If you’ve covered vectors before, then I think it should be reasonably easy to follow. Let me know if you have any questions or spot any mistakes.
We can represent our ray with an origin, and a direction. Both of these are 3-dimensional vectors. It will come in handy later if our direction is normalised (has a length of 1).
We can then represent the path the ray takes as a parametric equation in t.
(P is the position, O is the origin, and D is the direction).
Because our direction vector is normalised, t is also the distance the ray has traveled from the origin.
Representing a plane is a little less intuitive, but we can do it using a (3D) vector to represent the normal to the plane, as well as a value to represent the distance to the origin (in the direction of the normal). Note: this is the distance to the origin at (0,0,0) rather than the origin of the ray mentioned above.
e.g. if we look down the Z-axis from 0 in the positive direction, we could represent a wall 10 units away with a plane with a normal of (0,0,-1) and a distance of 10.
The general equation for a plane is:
Where a, b and c are the x, y and z components of the normal, and d is the distance to the origin. To make the later equations more concise, we’ll rewrite the plane equation in vector form as:
(P is the position, N is the normal, and d is the distance from the origin. Also the dot is the dot product)
As we’ve already established, t is the distance the ray has travelled from the origin. What we want to do now is find out how far the ray has travelled if/when it intersects our plane. We can do this by solving the ray and plane equations simultaneously.
Starting with the plane equation:
Substitute in the ray equation:
The dot product is distributive, so we can do this:
The dot product also lets us take multiplied scalar arguments outside the brackets like so:
Now shuffle it around a bit…
A bit more…
Super! Now we have the distance the ray has travelled when it hits the plane. We can plug this value back into the ray equation to get the intersection point in 3D coordinates too.
If the ray doesn’t hit the plane then both the numerator and denominator of the above equation are 0.
If you’d implemented what I’ve said so far, you would probably have a screen with half white (for the intersection hits) and half black (for the misses). So what we need now is some colour, or a slightly more interesting primitive to draw.
So, next time I’ll probably cover intersection with a sphere or how to colour your object with some simple lighting equations. Any preferences? I’ll probably do them both anyway.