astlib координирует код для координат астрограммы

Question

Я пытаюсь исправить этот код, который использовал astlib для преобразования координат для использования координатных модулей astropy. Вот копия кода:

from astLib import astCoords as coords 

#Convert B1950 coordinates (as given) 
clra=zeros(3) 
cldec=zeros(3) 
for i in range(len(clra)): 
    clra[i], cldec[i] = coords.convertCoords(
           'B1950', 'J2000', clra1950[i], 
           cldec1950[i], 1950 
          ) 

#Convert input coords to Galactic coords 
lgal,bgal = radians(coords.convertCoords(
         'J2000', 'GALACTIC', ra, 
         dec,2000 
        )) 

Мне нужна помощь с 2 вещами.

1. Каков собственный импорт для astropy, если я хочу, чтобы перейти от B1950 к J2000, и собственно импорту, чтобы перейти от J2000 в галактические координаты?

2. в областях, начинающихся с coords.convertCoords(). Каковы функции и аргументы от астрологии, которые выходят на свое место. Другими словами, что я могу заменить?

Кроме того, я провел некоторое исследование этой проблемы. Я нашел эту ссылку здесь: http://docs.astropy.org/en/v0.2.1/coordinates/

Он описывает обозначение для astropy's функции, связанные с координатами. Тем не менее, там много, и я не уверен, что использовать и как его использовать.

Answer 1

Это один из способов сделать это:

In [1]: import numpy as np 

In [2]: from astropy import units as u 

In [3]: from astropy.coordinates import SkyCoord, FK4, FK5, Galactic 

In [4]: clra = np.zeros(3) 

In [5]: cldec = np.zeros(3) 

In [6]: c1 = SkyCoord(clra * u.deg, cldec * u.deg, frame=FK4) 

In [7]: c2 = c1.transform_to(FK5(equinox='J2000')) 

In [8]: c3 = c2.transform_to(Galactic) 

In [9]: print(c3.l.degree) 
[ 97.74220094 97.74220094 97.74220094] 

In [10]: print(c3.b.degree) 
[-60.18102359 -60.18102359 -60.18102359] 

Обратите внимание, что вы можете конвертировать c1 в Galactic непосредственно тоже.