Seguindo o Diabo...

[Bruno Oliveira]

Segundo as próprias (fontes... ), estão a fazer trabalho de pesquisa para compilarem tudo aqui no fórum, para que mesmo quem percebe pouco de astronomia, se consiga desenrascar na programação .

Resta aguardar

[jap]

Óptimo!

[mav]

Vou tentar explicar o que são coordenadas equatoriais, horizontais (alt/azimutais) e o dia sideral. Se eu não for clara, ou se começar a usar termos astronómicos obscuros, não hesitem em perguntar.

Dia Sideral:

O dia sideral é o intervalo de tempo que decorre entre duas passagens consecutivas de uma estrela (que não o sol) sobre o ponto mais alto do meridiano local.

O meridiano local é a semi-circunfrência que passa pelo zénite, o ponto da esfera celeste que faz 90º como o horizonte, e o polo norte celeste (que é aproximadamente o sítio onde fica a estrela polar).



Devido ao movimento de translação da Terra, a duração de um dia sideral é menor que a duração de um dia solar. A duração de um dia sideral é 23h 50 min 4.1 s.
O tempo sideral é medido de acordo com o dia sideral e as 0h correspondem ao momento em que o ponto vernal cruza o meridiano local.


Coordenadas Horizontais

As coordenadas Horizontais, ou Alt/Azimutais, como o nome indica, tomam por base o horizonte, e variam consoante o local, o dia e a hora em que é feita a observação.
A altitude (h) corresponde à distância angular entre o astro e o horizonte. A altitude máxima é 90º e corresponde ao zénite.
O azimute (z) é o ângulo entre a direcção do astro (ou seja, do ponto de intercção do semi-circulo que passa pelo astro e pelo zénite com o horizonte) e o Norte. É medido no sentido dos ponteiros do relógio.




Coordenadas equatoriais


A utilização destas coordenadas têm a vantagem, das mesmas serem fixas independentemente do local de observação. São variáveis, no entanto, ao longo do tempo devido à precessão do eixo da Terra, mas para escalas de tempo pequenas permaneçem constantes.
As coordenadas equatoriais, são em tudo semelhantes às geográficas. O circulo fundamenteal deste sistema é o equador celeste, que é um círculo que faz um ângulo de 90º com os polos celestes.
São utilizadas duas coordenadas: A ascenção recta e a declinação.
A ascenção recta equivale à longitude e é medida de 0 h a 24 h (sendo que cada hora corresponde a 15º) no sentido dos ponteiros do relógio. O zero convencionou-se que era o ponto vernal, que é o primiero ponto da consteçaão carneiro e a posição em que o sol aparenta estar no equinócio da Primavera no hemisfério norte.
A declinação corresponde à latitude e é definida como a distância ângular de um astro ao equador celeste. Varia da 0º a 90º. A declinação do polo norte celestial é 90º+ e a do polo sul é 90-.



Conversão de coordenadas:

Azimute:Z
Altitude:h
Angulo horário: HA
Tempo sideral: TS
Ascensão recta: AR
Declinação: Dec
Latitude: B

HA=TS-AR

De coordenadas Horizontais para Equatoriais:

cos Dec * sin HA = cos h * sin Z
sin Dec = sin B * sin h + cos B * cos h * cos Z
cos Dec * cos HA = cos B * sin h + sin B * cos h * cos Z

De coordenadas Equatoriais para Horizontais:

cos h * sin Z = cos Dec * sin HA
sin h = sin B * sin Dec + cos B * cos Dec * cos HA
cos h * cos Z = - cos B * sin Dec + sin B * cos Dec * cos HA

Mais sobre as conversões:

http://www.seds.org/~spider/spider/ScholarX/coords.html

[jap]

Obrigado, Mariana!

Está óptimo.

Falta apenas explicar como se pode calcular o tempo sideral...

Bem, aqui vão as coordenadas equatoriais da nossa estrela algol:

Ascensão Recta: 03h 08m 10.1315s
Declinação: +40° 57′ 20.332″

Vamos lá então escrever um programa para calcular as coordenadas horizontais de Algol em Lisboa (ou noutro sítio à vossa escolha) usando a hora local dada pelo relógio do vosso computador...para saberem para onde apontar o vosso telescópio!

[Bruno Oliveira]

Óptimo

Obrigado Mariana, vou ler com mais calma assim que puder

[Bruno Oliveira]

Prof. então neste novo programa "não podemos" usar o novo módulo do Python, escrevemos o programa de raiz e depois comparamos o resultado com o que obtemos fazendo o programa com o módulo é isso?

[jap]

Prof. então neste novo programa "não podemos" usar o novo módulo do Python, escrevemos o programa de raiz e depois comparamos o resultado com o que obtemos fazendo o programa com o módulo é isso?

Sim, e hoje ainda colocarei aqui um post para ajudar nessa tarefa!
Stay tuned.

[Bruno Oliveira]

Obrigado, prof.

[jap]

Obrigado, prof.

Bruno, desculpa, ontem à noite fiquei sem net!

Então aqui vai a primeira tarefa: escrever uma função que calcule o tempo sideral local a partir da hora local do vosso computador. Não é difícil, podem encontrar toda a informação na net.

Lembro que o tempo sideral usa como referência um sistema 'fixo' nas estrelas distantes - a hora sideral não coincide com a hora 'solar'!
E podem testar o resultado do vosso programa com o Web clock deste link:

Relógio de tempo sideral na Web

Have fun!

[jap]

YALSTC:


http://www.jgiesen.de/astro/astroJS/siderealClock/

[jap]

Mais uma dica: calcular o tempo local sideral é trivial se conhecerem o GMST (Greenwish Mean Sideral Time).

E para calcular o GMST vejam aqui:

GMST

[jap]

Alguém tentou codificar um relógio sideral?

[jap]

Aqui vai a minha proposta...

Código: Seleccionar Todos

"""

Sideral Time Clock

@jap -2009 having fun with Quark!

"""

from time import gmtime,asctime
from math import *


def julian(y,m,d):
    """
    Julian day
    """
    if m > 2:
        m += 1
    else:
        y -= 1
        m += 13
    j = int(365.25*y)+int(30.6001*m)+d+1720982
    return j
       
   
def lst(lgt):
    """
    Sideral time
    """
    y,m,d,h,mn,sec = gmtime()[:6]
    j = julian(y,m,d)
    h += mn/60. +sec/3600. + lgt/15.
   
    jd0 = j -0.5 # julian date of yesterday at midnight
    jd = jd0 + h/24. # julian time
   
    d = jd - 2451545.0  # time elapsed since epoch = 12h 1st jan 2000
    d0 = jd0 - 2451545.0 # days since epoch = 12h 1st jan 2000
   
    t = d/3625. # number of centuries since epoch
   
    g = 6.697374558 + 0.06570982441908*d0+ 1.00273790935*h + 0.000026*t**2
   
    s = fmod(g,24)
   
    return s

def decimal(s):
    "Converts an hour in format h:m:s in decimal hours"
    s = s.split(':')
    if int(s[0]) < 0:
       sign = -1
    else:
       sign = 1
    return sign*(abs(int(s[0]))+int(s[1])/60.+float(s[2])/(60*60.))

def sexa(x):
    "Converts a decimal hour into format h:m:s"
    h = int(x)
    if h == 0:
        r = x
    else:
        r = abs(fmod(x,h))
    s = int(r*3600+.5)
    m,s = divmod(s,60)
    return `h`+':'+`m`+':'+`s`
       

if __name__ == "__main__":

    print "\nPythonic Sideral Clock\n"
    lisbon_lgt = decimal('-9:8:00')
   
    print "UTC time",asctime(gmtime())
    print "Greenwich mean sideral time",sexa(lst(0))
    print "Lisbon mean sideral time",sexa(lst(lisbon_lgt))