座標変換

　黄道直交座標（ｘ，ｙ，ｚ）を球面座標（ｒ，λ，β）に変換する。

　惑星の位置をベクトルで表して、P＝＜ｘ，ｙ，ｚ＞ とすれば、ｒ＝vlength(P) により、惑星までの距離ｒを得る。
　もし、(P.x＝0 & P.y＝0) ならば、H＝0、そうでなければ、H＝atan2(P.y ,P.x) とし、
λ＝mod(degrees(H)＋360 ,360) により、黄経λを得る。
　Rh=vlength(＜P.x ,P.y ,0＞) から、
　もし、(Rh＝0 & P.z＝0) ならば、V=0、そうでなければ、V=atan2(P.z ,Rh) とし、β=degrees(V) により、黄緯βを得る。

　vlength(Ｐ)は、ベクトルＰの長さを実数で返す。
　vlength(Ｐ)＝sqrt(pow(ｘ,2)＋pow(ｙ,2)＋pow(ｚ,2))
　POV-Rayでは、ベクトルＰの成分であるｘ、ｙ、ｚを、P.x、P.y、P.zにより、参照すことができる。