Understanding and deriving ellipsoidal coordinates geometrically If one were to read old texts on mathematical physics, like Maxwell, Morse & Feshbach, Hilbert and Courant, Jacobi, etc... they'd find ellipsoidal coordinates popping up, but the authors derive the coordinates using un-illuminating algebra. 
Is there a quick, geometric, way of writing down the Cartesian coordinates
$x^2 = (a^2+\xi)(a^2+\eta)(a^2+\zeta)/(b^2-a^2)(c^2-a^2)$
$y^2 = (b^2+\xi)(b^2+\eta)(b^2+\zeta)/(a^2-b^2)(c^2-b^2)$
$z^2 = (c^2+\xi)(c^2+\eta)(c^2+\zeta)/(a^2-c^2)(b^2-c^2).$
of the ellipsoidal coordinate system by drawing a picture, analogous to the way one derives the Cartesian components of the spherical system by drawing a sphere and projecting:

(It can be done quickly algebraically, but there's no geometry there)?
There seems to be an unintelligible, un-intuitive trigonometric method for the special case of oblate-spheroidal coordinates


(Pages 139 - 145 here give another un-intuitive description with a bad picture)
but this isn't the general case, thus I'm hoping for a picture that that literally copies the spherical picture above, including the details of what the $x$, $y$ & $z$ are by projections are, which these oblate pictures don't do.
If you have similar quick intuition as to why the scale factors and Laplacian take the forms they do, that would be a fantastically helpful bonus, thanks very much.
References: 


*

*Morse, Feshbach - "Methods of Mathematical Physics", P. 512

*Featherstone, Claessens - "Closed-Form Transformation Between Geodetic and Ellipsoidal Coordinates".

*Mathworld Confocal Ellipsoidal Coordinates Article

*Landau, - "Electrodynamics of Continuous Media", Section 4

*Tang, - "Mathematical Methods for Engineers and Scientists 2", P. 139 - 145

 A: $\def\a{\alpha}
\def\b{\beta}
\def\g{\gamma}
\def\th{\theta}
\def\f{\varphi}$The oblate spheroidal coordinates can be easily motivated geometrically. Begin with spherical coordinates. 
From the figure
\begin{align*}
x &= r\sin\th\cos\f \\
y &= r\sin\th\sin\f \\
z &= r\cos\th.
\end{align*}
We have 
$$x^2+y^2+z^2=r^2.$$
Scale the axes so that 
\begin{align*}
X &= \frac{\a}{r}x = \a\sin\th\cos\f \\
Y &= \frac{\b}{r}y = \b\sin\th\sin\f \\
Z &= \frac{\g}{r}z = \g\cos\th.
\end{align*}
We have 
$$\left(\frac{X}{\a}\right)^2+\left(\frac{Y}{\b}\right)^2+\left(\frac{Z}{\g}\right)^2=1,$$
so for a given $\a,\b,\g$ these are ellipses.
We require the ellipses to be confocal. 
Recall that the foci of the ellipse 
$(x/a)^2+(y/b)^2=1$ 
are located at $x=\pm c = \pm\sqrt{a^2-b^2}$, where we assume $a\geq b$. 
We require $\a\geq\b\geq\g$ and
$\a^2-\g^2=\b^2-\g^2=c^2$,
where $c$ is the linear eccentricity of the ellipses in the $xz$- and $yz$-planes, a constant.
Thus, 
$\a = \b = \sqrt{\g^2+c^2}$
and so 
\begin{align*}
X &= \sqrt{\g^2+c^2}\sin\th\cos\f \\
Y &= \sqrt{\g^2+c^2}\sin\th\sin\f \\
Z &= \g\cos\th.
\end{align*}
The coordinate system is given by $(\g,\th,\f)$.
