python 天文坐标转换

1.  将 ICRS坐标转换为 银道坐标

目的：对某颗星进行坐标转换

from astropy import units as u
from astropy.coordinates import SkyCoord
import astropy.coordinates as coord


ra = 117.2 * u.deg
dec = 20.34 * u.deg
c = SkyCoord(ra=ra, dec=dec, frame='icrs')
gc = c.transform_to(coord.Galactic)
print(gc.l, gc.b)

2.  将 银道坐标转换为 ICRS坐标

目的：想看一下 银道坐标下 80° < l < 220° , b = ± 5° 的区域在ICRS坐标下的对应情况

from astropy import units as u
from astropy.coordinates import SkyCoord
import astropy.coordinates as coord
import matplotlib.pyplot as plt
import numpy as np

l_t = np.arange(80,220,5)
b_t = 5*np.ones(len(l_t))
c = SkyCoord(l=l_t* u.deg, b=b_t* u.deg, frame='galactic')
gc1 = c.transform_to(coord.ICRS)

l_b = np.arange(80,220,5)
b_b = -5*np.ones(len(l_b))
c = SkyCoord(l=l_b* u.deg,