赤道和银道坐标转换公式

先记住几个基本常数：

北银极的赤道坐标：(alphaGP,deltaGP)=(192.85948,27.12825)单位是度,下同。
 
银心方向l=0,b=0。对应的赤道坐标(alpha,delta)=(266.405,-28.936)。

北天极的银经 lCP=122.932。

如果天体的赤道坐标为(alpha,delta),则其银道坐标(l,b)可由下面公式求得：

sin(b)=sin(deltaGP)sin(delta)+cos(deltaGP)cos(delta)cos(alpha-alphaGP)

cos(b)sin(lCP-l)=cos(delta)sin(alpha-alphaGP)

cos(b)cos(lCP-l)=cos(deltaGP)sin(delta)-sin(deltaGP)cos(delta)cos(alpha-alphaGP)

反向转换公式如下：
 
sin(delta)=sin(deltaGP)sin(b)+cos(deltaGP)cos(b)cos(lCP-l)

cos(delta)sin(alpha-alphaGP)=cos(b)sin(lCP-l)

cos(delta)cos(aplha-alphaGP)=cos(deltaGP)sin(b)-sin(deltaGP)cos(b)cos(lCP-l)