X=R cos d cos a
Y=R cos d sin a
Z = R sin d
Отсюда, R=sqrt(X^2+Y^2+Z^2)
tg a = Y/X, четверть, в которой расположен угол a, определяется знаками Y и X. Также можно воспользоваться функцией atan2(Y,X) во многих языках программирования.
d = arcsin(Z/R)
Только помни, что системы координат R,a,d (сферическая) и X,Y,Z (прямоугольная) должны иметь один центр (геоцентрицеские), одну экваториальную плоскость и одно начало отсчета азимутальной координаты (a, X). Если X,Y,Z гелиоцентрицеские, то нужно преобразовать в геоцентрические: X'=X+Xo, Y'=Y+Yo, Z'=Z+Zo, где Xo - геоцентрицеские координаты Солнца, они связаны простым соотношением с гелиоцентрическими координатами Земли: Xo=-Xз