Circle Shadow

Geometry Level 5

A circle of radius 10 units is centered at ( 0 , 0 , 50 ) (0, 0, 50) . The plane it lies in, has a normal vector ( 1 , 0 , 0 ) (1, 0, 0) . A uniform direction ray of light falls on the circle, generating a shadow on the floor (which is the XY plane). The direction of the shadow is given by d = ( 1 , 0.5 , 0.75 ) d = (1, -0.5, -0.75) . The shadow is an ellipse. Find the sum of the semi-major and semi-minor axes of this ellipse.


The answer is 24.267.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

1 solution

Hosam Hajjir
Jun 5, 2016

Points on the circle can be represented as

r C = r 0 + V u r_C = r_0 + V u

where r 0 r_0 is the center of the circle and V = [ v 1 , v 2 ] V = [ v_1, v_2 ] , where v 1 v_1 , v 2 v_2 are two orthogonal vectors of length 10 that are orthogonal to the normal to the plane of the circle, and u = [ cos t , sin t ] T u = [\cos t, \sin t]^T .

The equation of the projection plane is

n p T ( r r p ) = 0 n_p^T ( r - r_p ) = 0

Now the projected point is generated from the original point on the circle as follows

r = r C + t d r' = r_C + t d

Hence,

n p T ( r c + t d r p ) = 0 n_p^T ( r_c + t d - r_p) = 0

Therefore,

t ( n p T d ) + n p T ( r C r p ) = 0 t (n_p^T d) + n_p^T(r_C - r_p) = 0

which implies,

t = n p T ( r C r p ) n p T d t = -\dfrac{n_p^T (r_C - r_p)}{ n_p^T d}

Therefore, the projected point onto the plane is given by,

r = r C n p T ( r C r p ) n p T d d r' = r_C - \dfrac{n_p^T(r_C - r_p)}{ n_p^T d} d

r = r C d n p T ( r C r p ) n p T d r'= r_C - \dfrac{d n_p^T (r_C - r_p) }{ n_p^T d}

r = r p + ( I d n p T n p T d ) ( r C r p ) r'= r_p + (I - \dfrac{d n_p^T }{ n_p^T d} ) (r_C - r_p )

which can written in compact form as:

r = r p + P ( r C r p ) r' = r_p + P (r_C - r_p)

with P = I d n p T n p T d P = I - \dfrac{d n_p^T}{ n_p^T d} .

Next, we note that points on the projection plane can be represented by

r = r p + R p w r = r_p + R_p w

where R p R_p is a 3 × 3 3 \times 3 matrix comprising a reference frame attached to the projection plane and has its origin at r p r_p . In particular,

R p = [ u 1 , u 2 , n p ] R_p = [u_1, u_2, n_p]

where n p n_p is the unit vector normal to the plane, and u 1 u_1 , u 2 u_2 are two orthogonal unit vectors that are also orthogonal to n p n_p .

w w is the coordinate of a given point with respect to this frame. Therefore,

r = r p + P ( r C r p ) = r p + R p w r' = r_p + P (r_C - r_p) = r_p + R_p w

Simplifying and substituting r C r_C , we have

P ( r 0 + V u r p ) = R p w P ( r_0 + V u - r_p) = R_p w

Solving for w w , and noting that R p 1 = R p T {R_p}^{-1} = R_p^T , we get,

w = R p T P ( r 0 r p ) + R p T P V u w = {R_p}^T P (r_0 - r_p) + {R_p}^T P V u

Define

h = R p T P ( r 0 r p ) h = {R_p}^T P (r_0 - r_p) and G = R p T P V G = {R_p}^T P V , then

w = h + G u w = h + G u

One can verify that the third row of R p T P {R_p}^T P is a row of zeros, making the third row of h h and G G a row of zeros as well. This is consistent with the fact that the third element of w w should be zero, because the projected point lies in the projection plane. We will therefore consider only the first two rows of the above equation, thus obtaining,

w = h + G u w' = h' + G' u

where the primed variables are obtained from the original variables by taking only the first two rows.

Now, using the fact that u T u = 1 u^T u = 1 , we get

( w h ) T G T G 1 ( w h ) = 1 (w' - h')^T G'^{-T} G'^{-1} (w' - h') = 1

At this point we can find a rotation matrix R , and a diagonal matrix D such that

R T G T G 1 R = D R^T G'^{-T} G'^{-1} R = D

Therefore, by a change of variable, w = R p + h w' = R p + h' , we obtain

p T D p = 1 p^T D p = 1

This is an equation of an ellipse, with semi-major and semi-minor axes given by the reciprocal of the square roots of the diagonal elements of D.

Performing the necessary calculations, we find that the semi-major axis M, and the semi-minor axis m, are given by

M = 15.860296 M = 15.860296

m = 8.4067365 m = 8.4067365

Therefore, the answer is M + m = 24.267 M + m = \boxed{24.267}

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...