Seguindo o Diabo...

[jap]

Salvo seja!

Escrevam um programito que permita fazer o tracking de Algol, a "estrela do Diabo", e também, por razões que bem conhecem, a estrela fetish do Quark!

Pretende-se o seguinte: dadas a latitude e longitude do lugar e a hora local (lida no computador), o programa deve indicar as coordenadas (altura e azimute) da estrela Algol, para onde devem apontar o vosso telescópio.

Depois comparam com resultado dado pelo PyEphem.

Se tiverem dúvidas, a Mariana esclarece.

Have fun!

[Bruno Oliveira]

Então não podemos usar o PyEphem no nosso programa?

[jap]

Podem usar o PyEphem, claro, mas a tarefa torna-se então bastante simples.

Mais challenging é fazerem um outro programa que efectue os cálculos sem recurso ao PyEphem - não é assim tão difícil, se souberem um pouco de astronomia e, nesse caso, podem comparar os resultados dos dois programas.

Aceito, pois, ambas aproximações ao problema.

[Bruno Oliveira]

Eu estava a tentar inicialmente fazer com o recurso ao módulo mas não o estou a conseguir instalar

EDIT: Já funciona

[Bruno Oliveira]

Até agora o que já consegui foi isto:

Código: Seleccionar Todos

import ephem

algol=ephem.star('Algol')   #Busca da estrela Algol na base de dados de estrelas do módulo

Lisbon=ephem.city('Lisbon')
print "Longitude de Lisboa:",Lisbon.long,"    e      ", "Latitude de Lisboa:",Lisbon.lat

lisbon=ephem.Observer()

lisbon.lat=Lisbon.lat    #Latitude de Lisboa
lisbon.long=Lisbon.long   #Longitude de Lisboa

algol.compute(lisbon)
print "Azimute de Algol em relação a Lisboa é:",algol.az
print "Altitude de Algol em relação a Lisboa é:", algol.alt


Este código usa inteiramente, o novo módulo o PyEphem

O output é:

Código: Seleccionar Todos

>>> ================================ RESTART ================================
>>>
Longitude de Lisboa: -9:08:00.0     e       Latitude de Lisboa: 38:43:00.0
Azimute de Algol em relação a Lisboa é: 310:12:53.7
Altitude de Algol em relação a Lisboa é: 15:58:17.7
>>>

Mas, como não estou a conseguir entrar com a hora local, deve estar algo a falhar...
Alguém para ajudar no debug, ou comentar o código?

[jap]

Bruno,

Quando crias uma cidade, o objecto criado é uma instância da classe "Observer" e tem guardado no atributo 'date' a hora local no instante em que o objecto é criado. No cálculo (compute) a conversão da hora local para UTC, etc, é feito automaticamente.

Assim, basta escrever:

Código: Seleccionar Todos

import ephem

print "Finding the devil's star"

algol = ephem.star('Algol')
lisbon = ephem.city('Lisbon')

algol.compute(lisbon)

print "To observe Algol in Lisbon, a lovely city located at"
print "lat=",lisbon.lat,"long=",lisbon.long
print "you should point your telescope to alt=",algol.alt,"azimuth=",algol.az
print "at this local time:",lisbon.date
print "Good luck and take a photo of the devil's star!"

Cujo resultado é:

Código: Seleccionar Todos

Finding the devil's star
To observe Algol in Lisbon, a lovely city located at
lat= 38:43:00.0 long= -9:08:00.0
you should point your telescope to alt= 82:48:57.5 azimuth= 291:16:56.2
at this local time: 2009/3/1 17:42:11
Good luck and take a photo of the devil's star!

Repara que podes passar directamente o objecto 'lisbon' a 'compute'. Não te esqueças porém que se voltares a utilizar o objecto lisbon a hora que ele guarda é a hora da sua criação, que não é actualizada automaticamente. Se quiseres actualizar o atributo 'date' do objecto 'lisbon' tens de fazer:

Código: Seleccionar Todos

lisbon.date = ephem.now()

que irá ler de novo a hora ao relógio do teu computador.

Simples, não?

[Bruno Oliveira]

Obrigado pela explicação prof. .

Então isso quer dizer que o resultado que eu obtive,pode ter sido o certo, para a hora a que corri o código?

[jap]

Obrigado pela explicação prof. .

Então isso quer dizer que o resultado que eu obtive,pode ter sido o certo, para a hora a que corri o código?

Sim, é isso mesmo. Sempre que correres o programa ele dará o resultado correcto para essa hora.

[Bruno Oliveira]

Óptimo , vou então completar o meu código!

[Bruno Oliveira]

E...já que fizemos o tracking da Algol, a partir de Lisboa, porque não deixar aqui o tracking da Devil`s Star a partir de Coimbra?

A longitude e a latitude da cidade mítica sabem-se no Google Earth, mas alguém sabe o valor da altitude?

[Bruno Oliveira]

Bom... Depois de uma pesquisa, encontrei, num site de meteorologia, que a altitude média é de 44 m.

Adicionando, então, o nome da cidade e as suas coordenadas geográficas, ao ficheiro cities.py na pasta Ephem, podemos então, em tempo real, fazer o tracking da Algol, a partir de Coimbra:

Código: Seleccionar Todos

import ephem

algol=ephem.star('Algol')

Coimbra=ephem.city('Coimbra')

print "Longitude de Coimbra:",Coimbra.long,"    e      ", "Latitude de Coimbra:",Coimbra.lat

Coimbra=ephem.Observer()

Coimbra.lat=Coimbra.lat
Coimbra.long=Coimbra.long

algol.compute(Coimbra)

print "Para observarem Algol, a partir da cidade mítica, apontem os vossos telescópios para:"
print  algol.az
print algol.alt
print "Valores para a hora local:", Coimbra.date


[Bruno Oliveira]

Prof. depois se pudesse e quando pudesse claro, poderia dar umas dicas acerca do programa que não usa o módulo?

Obrigado, Bruno

[jap]

Sim, farei isso em breve.

Será que, entretanto, algum astrónomo(a) amador(a) não gostará de explicar aqui o que são as coordenadas equatoriais e alt-azimutais, como se passa de umas para outras, e o que é o tempo sideral?

[Bruno Oliveira]

(...)
Será que, entretanto, algum astrónomo(a) amador(a) não gostará de explicar aqui o que são as coordenadas equatoriais e alt-azimutais, como se passa de umas para outras, e o que é o tempo sideral?

Alô Mariana...

[jap]

Onde estão os astrónomos(as)?