Me gustaría saber qué parte de la Vía Láctea está arriba actualmente. Preferiblemente usando astropía. Los ingredientes que tengo son el tiempo actual y mi posición actual en la tierra. ¿Cómo los usaría para obtener la longitud galáctica mínima y máxima actualmente?
(Sé cómo convertir de Alt-Az a coordenadas galácticas, pero no sé en qué azimut el plano galáctico cruza el horizonte en altitud = 0).
Lo siguiente funciona para mí; utiliza el rootfinder de scipy para encontrar el cero.
from astropy.coordinates import EarthLocation, AltAz, Galactic
from astropy.time import Time
from scipy.optimize import newton
import astropy.units as u
dwl = EarthLocation(lat=6.39688*u.deg, lon=52.6275*u.deg)
mytime = Time('2017-09-16T21:00:00')
def altitude_of_longitude(longitude):
''' Give the altitude (in deg) of a galactic longitude (in deg) '''
return Galactic(l=longitude*u.deg, b=0*u.deg).transform_to(
AltAz(location=dwl, obstime=mytime)).alt.value
Ahora use el buscador de raíces de scipy, con dos conjeturas iniciales diferentes, con la esperanza de que terminen en las dos soluciones diferentes donde el plano galáctico cruza el horizonte. Esto da:
newton(altitude_of_longitude, 0)
>>> 16.231753986033333
newton(altitude_of_longitude, 180)
>>> 196.23547201834657
La única parte restante es verificar qué parte es visible:
altaz = Galactic(l=0*u.deg, b=0*u.deg).transform_to(AltAz(location=dwl, obstime=mytime))
altaz.alt
>>> −9°08′51.8623′′
Aparentemente la longitud a 0° está por debajo (altitud negativa), por lo que la parte visible está entre 16° y 196° de longitud galáctica.