How to Parametrize a 3D Circle and Find the Closest Point?

[rod bryant]

I have a circle defined by the intersection of a sphere and a plane.* I want to convert from cartesian to parametric form on my way to finding the closest point on the circle to some arbitrary other point.* How do I go about this?* Alternatively, how do I find the closest point?
The equation of the sphere is x^2 + y^2 + z^2 = R^2.
The equation of the plane is:* xsub0*x + ysub0*y + zsub0*z = R^2 - r^2/2 where (xsub0,ysub0,zsub0) happens to be another point in space.

I want to find the closest point on the circle to (xsub1,ysub1,zsub1).

 

[maze]

Hmm, I would do it a different way.

1) Find the vector to the center of the circle, \vec{c}
2) Find 2 perpendicular unit length vectors that lie on the plane, \vec{v}_1,\vec{v}_2
3) Find the radius of the circle, r (perhaps different from R if the plane doesn't slice through the exact center of the sphere)
4) Circle is given by \vec{c}+r cos\left(t\right) \vec{v}_1+r sin\left(t\right) \vec{v}_2

 

[HallsofIvy]

Hmm, I would do it a different way.

1) Find the vector to the center of the circle, \vec{c}

The vector to the center of the circle from where? And how do you find the center of the circle?

2) Find 2 perpendicular unit length vectors that lie on the plane, \vec{v}_1,\vec{v}_2
3) Find the radius of the circle, r (perhaps different from R if the plane doesn't slice through the exact center of the sphere)
4) Circle is given by \vec{c}+r cos\left(t\right) \vec{v}_1+r sin\left(t\right) \vec{v}_2

Pretty much the same thing: From the equation of the plane and knowing the center of the sphere is at (0,0,0), you can find the equation of the line through (0,0,0) perpendicular to the plane. Determine the point at which that line crosses the plane. That will be the center of the circle (and may be what maze meant in (1) above). Choose any point on the circle and let \theta be the angle a radius of the circle makes with the line through the center of the circle and that point. You can use \theta as parameter and add vectors to get the complete parametric equations.

Now, for the actual problem. To find the closest point on the intersection of a sphere and a circle to (x_0,y_0,z_0) I would minimize (x- x_0)^2+ (y- y_0)^3+ (z-z_0)^2 using Lagrange multipliers with the equations of the sphere and plane as constraints.

 

[maze]

Here is a diagram:
http://img105.imageshack.us/img105/7691/sphereplaneln3.png

Also,
If you know any 3 points on the circle \vec{p}_1,\vec{p}_2,\vec{p}_3 (they should be easy to find), then \vec{c} and r may be found as follows:
\vec{c} = \frac{1}{2}\left(\begin{matrix}\vec{p}_2-\vec{p}_1 \\ \vec{p}_3-\vec{p}_1 \\ \vec{p}_3-\vec{p}_2\end{matrix}\right)^{-1}\left(\begin{matrix}\vec{p}_1\cdot\vec{p}_1-\vec{p}_2\cdot\vec{p}_2 \\ \vec{p}_1\cdot\vec{p}_1-\vec{p}_3\cdot\vec{p}_3 \\ \vec{p}_2\cdot\vec{p}_2-\vec{p}_3\cdot\vec{p}_3 \end{matrix}\right)

r =\left|\vec{c}-\vec{p}_1\right|

For derivation of these equations, see this thread: https://www.physicsforums.com/showthread.php?t=173847

 

[maze]

Also, here is the solution for the closest point:
http://img142.imageshack.us/img142/483/closestptyv6.png

 

[rod bryant]

parameterizing a circle in 3D

Thanks maze and hallsofivy. While waiting for a response I went about the problem another way via an approximation. Now I can compare the exact solution.