Ray Polygon Interception Cyrus Beck Algorithm Mengxia Zhu

Ray Polygon Interception: Cyrus Beck Algorithm Mengxia Zhu Fall 2007

Polygon and Polyhedron n n Polygon and Polyhedron are fundamental objects in 2 D or 3 D graphics. In 2 D, straight lines make a polygon, in 3 D, an object if often modeled as a polygonal “mesh”: a collection of polygons that fit together. If the skin forms a closed surface, the mesh is called polyhedron. P

Ray Polygon Intersection n In order to analyze or render them, we usually need to solve the intersection problem. The general case of intersecting a line with a polygon and polyhedron is quite complex. We study convex forms with bounding lines or bounding planes first. Solve how a ray intersects a convex polygon or convex polyhedron is an essential technique in viewing 3 D objects. ¨ ¨ ¨ Is a given point P inside or outside the object Where does a given ray first intersect the object Which part of a given line lies inside the object?

Ray Casting Ray Tracing examples n n www. geocities. com/. . . / coolhead 1 small 2 rot. png www-evasion. imag. fr/. . . /ecs 177. html

Parametric Line Revisited • P 0 is the starting point • P 1 another point on the line • P linear interpolation between P 0 and P 1 • Affine combination, Convex combination

Dot Product Revisited

Point inside or outside?

Direction for Ray that Strikes Plane r n n n the ray is aimed “along with” the normal the ray is parallel to the plane the ray is aimed “counter to” the normal

Cyrus-Beck Algorithm n n n Developed by Cyrus and Beck in 1978 Clip line against convex polygon/polyhedron Completely described by a set of “bounding lines or planes” Find a value of t for intersection of line and clipping edges. Use simple comparison of t values to find actual clipping segment.

A point normal form of Where a Ray edge or planeintersects a line(2 D) or n plane(3 D) n n A parametric equation of line: Find t such that Outside P(t) PL Inside NL

Cyrus-Beck Algorithm (Cont’) n Classify the relative position of the ray to intercepting edge or plane ¨ If : ray parallel to the edge of line or plane ¨ If : ray exiting polygon ¨ If : ray entering polygon Entering P 0 Exiting P 1

Each Bounding Line or Plane n We should find: ¨ The hit time of the ray with the bounding line or bounding plane ¨ Whether the ray is entering or exiting the polygon or polyhedron.

Cyrus-Beck Algorithm (Cont’) Compute PE with largest t: max (tin, 0) n Compute PL with smallest t: min ( tout, 1) P 1 n Clip to these two points n PL PL PE P 0 PE

Cyrus-Beck Algorithm (Cont’) n n Initialize the candidate interval to [0, 1] For each bounding line, find the hit time of thit and determine whether it is an entering or exiting hit ¨ If t_hit is an entering hit, set t_in = max (tin, thit) ¨ If t_hit is an exiting hit, set t_out = min (tout, thit) If at any point tin becomes greater than t_out, the ray misses P entirely, and testing is terminated. If the candidate interval is not empty, then the segment between tin and tout lies inside P.

Clipping Against Arbitrary Polygons Polygon is no longer bounded by a collection of infinite bounding lines in point normal form. n We have to work with N finite segments n P 2 P 1

Clipping Against Arbitrary Polygons n Represent edge of P parametrically For example, segment P 3 P 4 is represented as: P 3 + e 3*u, where e 3=P 4 -P 3, u ranges from 0 to 1. ¨ Ray intersection with above segment is represented as: R + c*t = P 3 + e 3*u ¨ n Find all of the hits of the ray with all edges of P and place them in a list of the hit times. n The smallest value of t in hit-list is the first hit n Take the t values in pairs. The ray enters P at the first time in each pair and exits from P at the second time of each pair

Exercises n n n Draw an icosahedron by using the following drawing routine to use vertex arrays. Icosahedron approximate a sphere by using 20 polygons. After finish it, try to use recursive method to further divide each triangle to get more triangles to better approximate the sphere. #define x. 525731112119133606 #define z. 850650808352039932 void display (void) {static GLfloatv data[12][3] = { {-x, 0. 0, z}, {-x, 0. 0, -z}, {x, 0. 0, z}, {0. 0, z, x}, {0. 0, z, -x}, {0. 0, -z, -x}, {z, x, 0. 0}, {-z, x, 0. 0}, {z, -x, 0. 0}, {-z, -x, 0. 0} }; static GLuinttindices[20][3] = { {0, 4, 1}, {0, 9, 4}, {9, 5, 4}, {4, 5, 8}, {4, 8, 1}, {8, 10, 1}, {8, 3, 10}, {5, 3, 8}, {5, 2, 3}, {2, 7, 3}, {7, 10, 3}, {7, 6, 10}, {7, 11, 6}, {11, 0, 6}, {0, 1, 6}, {6, 1, 10}, {9, 0, 11}, {9, 11, 2}, {9, 2, 5}, {7, 2, 11} };

Exercises Implement a 2 D/3 D Cyrus-Beck Algorithm for convex polygon or polyhedron by prompt user to input a ray and a sequence of bounding lines/bounding plane. n Output the enter point and exit point. n

Exercises for (i = 0; i < 20; i++) { gl. Vertex 3 fv(&vdata[tindices[i][0]); gl. Vertex 3 fv(&vdata[tindices[i][1]][0]); gl. Vertex 3 fv(&vdata[tindices[i][2]][0]); }