座標系轉換

參考：

關於天文座標系的介紹http://astroleaks.lamost.org/2012/01/22/%E5%85%B3%E4%BA%8E%E5%A4%A9%E6%96%87%E4%B8%AD%E7%9A%84%E5%9D%90%E6%A0%87%E7%B3%BB%E7%9A%84%E4%BB%8B%E7%BB%8D/

X射線脈衝星導航可見性分析(毛悅等) http://ch.whu.edu.cn/EN/article/downloadArticleFile.do?attachType=PDF&id=1179

赤道 銀道 座標系轉換 https://blog.csdn.net/casularm/article/details/302811

赤道座標轉銀道座標的python程序 http://blog.sina.com.cn/s/blog_a046022d0101a6hs.html

 

嘗試代碼，error

import math
from math import pi
from math import sin
from math import cos

i0=62.6
i=i0*pi/180       #角度轉弧度
omega0=282.25
omega=omega0*pi/180
Lomega0=33.345
Lomega=Lomega0*pi/180

l0=float(raw_input('l0 = '))
b0=float(raw_input('b0 = '))
l=l0*pi/180
b=b0*pi/180

sin_delta=sin(b)*cos(i)+cos(b)*sin(i)*sin((l0-Lomega0)*pi/180)
print 'sin_delta=', sin_delta
delta=math.asin(sin_delta)
print 'delta 弧度=', delta
delta=delta*180/pi
#delta=180+delta
print 'delta 角度=', delta

cos_alpha_cut_omega=(cos(b)*cos(l-Lomega))/cos(delta)
print 'cos_alpha_cut_omega=', cos_alpha_cut_omega
alpha_cut_omega=math.acos(cos_alpha_cut_omega)
print 'cos_alpha_cut_omega 弧度=', alpha_cut_omega
alpha_cut_omega=alpha_cut_omega*180/pi
print 'cos_alpha_cut_omega 角度=', alpha_cut_omega
alpha=alpha_cut_omega+omega0
print 'alpha=', alpha