Como calcular (alt, az) para determinada coordenada galáctica (GLON, GLAT) com PyEphem?

Para um dado observador (lon, lat, time) na Terra e uma dada coordenada galáctica (GLON, GLAT), como posso calcular o ponto correspondente (alt, az) no céu com PyEphem?

Dada a maneira como o PyEphem funciona atualmente, há duas etapas para responder à sua pergunta. Primeiro, você tem que converter um par de coordenadas galácticas em um par equatorial de coordenadas RA / dec.

import ephem # Convert a Galactic coordinate to RA and dec galactic_center = ephem.Galactic(0, 0) eq = ephem.Equatorial(galactic_center) print 'RA:', eq.ra, 'dec:', eq.dec → RA: 17:45:37.20 dec: -28:56:10.2 

Essas coordenadas são muito próximas das que a Wikipedia dá para o centro galáctico.

Agora que temos uma coordenada RA / dec normal, só precisamos descobrir onde está no céu agora. Como PyEphem é construído sobre uma biblioteca que conhece apenas corpos celestes como estrelas e planetas, nós simplesmente precisamos criar uma “estrela” falsa naquele local e perguntar seu azimute e altitude.

 # So where is that RA and dec above Boston? # Pretend that a star or other fixed body is there. body = ephem.FixedBody() body._ra = eq.ra body._dec = eq.dec body._epoch = eq.epoch obs = ephem.city('Boston') obs.date = '2012/6/24 02:00' # 10pm EDT body.compute(obs) print 'Az:', body.az, 'Alt:', body.alt → Az: 149:07:25.6 Alt: 11:48:43.0 

E nós podemos verificar que esta resposta é razoável, olhando para um mapa do céu para Boston no final da noite: Sagitário – a localização do Centro Galáctico – está apenas subindo sobre a borda sudeste do céu, que faz perfeito sentido de um azimute sudeste 149 °, e por uma altitude ainda baixa como 11 °.

    Intereting Posts