# -*- coding: utf-8 -*-
"""
Created on Sun Mar 12 20:28:24 2017
@author: rodolfo
"""
from datetime import datetime
import astropy.units as u
from astropy.time import Time
from astropy.coordinates import SkyCoord, EarthLocation, AltAz, Angle
#Locación
Puebla = EarthLocation(lat=Angle('19d02m36s'),
lon=Angle('-98d21m6.94s'),
height=2152*u.m)
#Fecha y hora
utc_time = Time(datetimeutcnow(), scale='utc')
time = utc_time
#Ecuatoriales
Sirio = SkyCoord(Angle('101d17m18.7469s'),
Angle('-16d42m47.3141s'),
frame='icrs')
#Horizontales
Sirio = Sirio.transform_to(AltAz(obstime=time,location=Puebla))
print "Altura: ", Sirio.alt.degree
print "Azimut; ", Sirio.az.degree
domingo, 12 de marzo de 2017
Convertir coordenadas ecuatoriales a altazimutales (horizontales) con astropy
Algunos meses después de entender la geometría y hacer un programa para convertir entre estas coordenadas para mi tesis de licenciatura, descubrí que las conversiones podían hacerse con el módulo de astropy. No fue tiempo perdido pues entendí bastante bien el proceso. Para una conversión de muy buena precisión recomendaría usar astropy aunque sus funciones tienen el inconveniente de ser relativamente lentas (promedian una ejecución de 3 ms).
Suscribirse a:
Enviar comentarios (Atom)
No hay comentarios:
Publicar un comentario