viernes, 24 de marzo de 2017

Spoon - Kill The Moonlight (2002)


Este es sin dudas un disco que podría recomendarle a cualquier persona con un sello de garantia y, considerando mis gustos, esto no es algo que podría decir de cualquier otro. Llevo ya algunas semanas escuchándolo. Me parece que tiene un sonido muy accesible y al mismo tiempo una calidad que bien podría ponerse a un lado de un Yankee Hotel Foxtrot o un The Moon & Antartica. Me costó bastante encontrarlo pero la ardua búsqueda entre enlaces muertos valió mucho la pena. Era justo ponerlo al alcance de todos por aquí.

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 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).

# -*- 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