Soy nuevo en el uso de los paquetes de Python astropy/astroqueryy lo que quiero lograr es usarlo para encontrar estrellas que se encuentran en la zona de vista continua (CVZ) de TESS.
Sé que la CVZ está centrada en el polo sur de la eclíptica y me gustaría usar la función astroquery.mast.Observations.query_region()para consultar una región centrada en las coordenadas del polo. Dicha función usa coordenadas en el sistema internacional de referencia celeste (ICRS), por lo que primero tengo que convertir el par desde el marco de la eclíptica.
.
Intenté la conversión usando la función transform_to(ICRS), sin embargo, devuelve el siguiente error:
UnitsError: las coordenadas HeliocentricTrueEcliptic de entrada no tienen unidades de longitud. Esto probablemente significa que creó coordenadas con lat/lon pero sin distancia . Las transformadas heliocéntricas<->ICRS no pueden funcionar en este caso porque hay un cambio de origen .
¿Alguien podría elaborar un poco más sobre la "distancia" y el "cambio de origen" mencionados en el mensaje de error y cómo corregirlo?
Esto se debe a que el ICRS tiene un origen en el baricentro y este no es lo mismo que el centro del Sol (el heliocentro). El heliocentro está desplazado del baricentro (aunque todavía está dentro del radio del Sol) y varía con el tiempo, principalmente debido a las perturbaciones de Júpiter y Saturno.
Si lo usas BarycentricTrueEclipticfuncionará bien:
from astropy.coordinates import BarycentricTrueEcliptic
from astropy import units as u
cvz = BarycentricTrueEcliptic(lon=0*u.deg, lat=-90*u.deg)
cvz.transform_to(ICRS)
Out[9]:
<ICRS Coordinate: (ra, dec) in deg
(89.9999853, -66.56071866)>
Lo cual está de acuerdo con las conversiones de coordenadas usando el LAMBDAsitio web en HEASARC de la NASA