python天文坐标系旋转换算程序

python天文坐标系旋转换算
处理给定坐标后，在X\Y\Z方向进行坐标旋转后的角度输出的简易函数

#《应用天体测量》课程训练
#具体请参考《球面天文学》等书目
import numpy as np


R=input('please input 旋转方向(xyz):')
θ=float(input('please input 旋转角度:'))
μ=float(input('please input 赤经μ:'))
γ=float(input('please input 赤纬γ:'))
 
x=np.cos(μ*np.pi/180)*np.cos(γ*np.pi/180)
y=np.sin(μ*np.pi/180)*np.cos(γ*np.pi/180)
z