2015-11-07 5 views
4

J'utilise pyEphem pour obtenir le prochain par le col de l'ISS de mon emplacement, mais les résultats mon retour ne correspond pas à ce que je vois sur ciel au-dessus en utilisant les mêmes coordonnéespyephem next_pass pas cieux au-dessus correspondant à

Je suis probablement fait une erreur stupide, mais je ne peux pas le comprendre

Mon code renvoie ci-dessous les résultats: temps de montée: 25/05/2017 azimut 20:34:39: 193: 28: 04,0

Alors que le plus haut ci-dessus passe est près de 3 heures de route, avec un temps de montée de: 23:09:40

http://www.heavens-above.com/passdetails.aspx?&satid=25544&mjd=57898.9270155034&type=V

from datetime import datetime 
import ephem 
import pytz 

line1 = 'ISS (ZARYA)'  
line2 = '1 25544U 98067A 17145.52800275 .00016717 00000-0 10270-3 0 9015' 
line3 = '2 25544 51.6372 151.2656 0005033 192.5139 167.5889 15.53913304 18224' 

tle = [line1, line2, line3] 
iss = ephem.readtle(tle[0], tle[1], tle[2]) 

longitude = -6.2282 
latitude = 53.2842 
altitude = 20 

site = ephem.Observer() 
site.lat = str(latitude) 
site.lon = str(longitude) 
site.elevation = 20 

current_time = datetime(2017, 5, 25, 12, 0, 0, tzinfo=pytz.utc) 
site.date = current_time 

info = site.next_pass(iss) 
print("Rise time: %s azimuth: %s" % (info[0], info[1])) 
+0

Pour permettre à d'autres personnes d'exécuter ce code, vous devez supprimer l'appel 'utcnow()' et spécifier à la place une date et une heure exactes. Sinon, quelqu'un qui exécute le script verra toujours une réponse différente du résultat que vous avez partagé dans votre réponse, car il l'exécutera à une date différente de celle que vous avez utilisée. –

Répondre

2

réponse Mise à jour suivante modifier la question:

Hélas, Heavens Above donne maintenant une erreur de serveur quand je suis le lien dans votre question. J'ai donc visité le site à nouveau, entré les coordonnées de votre script, et j'ai eu des prédictions. Pour éviter de les perdre, voici une capture d'écran rapide:

enter image description here

Lorsque je colle dans une nouvelle TLE pour l'ISS dans votre script et régler la current_time, je peux obtenir une réponse qui est assez proche de la leur. Le script ressemble alors à:

from datetime import datetime 
import ephem 
import pytz 

line1, line2, line3 = """\ 
ISS (ZARYA)    
1 25544U 98067A 17198.89938657 .00000988 00000-0 22167-4 0 9998 
2 25544 51.6416 245.2318 0005849 47.2823 302.7554 15.54170925 66526 
""".splitlines() 

tle = [line1, line2, line3] 
iss = ephem.readtle(tle[0], tle[1], tle[2]) 

longitude = -6.2282 
latitude = 53.2842 
altitude = 20 

site = ephem.Observer() 
site.lat = str(latitude) 
site.lon = str(longitude) 
site.elevation = 20 

current_time = datetime(2017, 7, 21, 1, 0, 0, tzinfo=pytz.utc) 
site.date = current_time 

info = site.next_pass(iss) 
for item in info: 
    print(item) 

Et la sortie est:

2017/7/21 01:32:01 
263:02:11.2 
2017/7/21 01:37:29 
66:12:31.4 
2017/7/21 01:42:58 
93:15:48.0 

Quelle est la même passe que vous pouvez voir dans la capture d'écran, en temps universel au lieu de l'UTC + 1 fuseau horaire pour Dublin que Heavens Above utilise - donc le temps le plus élevé de 1:37 donné par PyEphem devient 02h37 dans le fuseau horaire local de Dublin.

J'ai vérifié un autre passage ou deux, et ils semblent tous être assez proches - le fuseau horaire aurait-il pu être une source de confusion?

réponse originale:

Vous spécifiez une longitude de 53,2842 ° Est et une latitude de 6,2282 ° sud, sur une carte du monde semble être quelque part à l'extrémité ouest de l'océan Indien. Avez-vous plutôt prévu le nombre négatif de -53.2842 °, ce qui placerait l'emplacement au Brésil à la place?

+0

Désolé je suis un idiot, il était censé être latitude: 53.2842 °, longitude: -6.2282, ce qui me place à Dublin, en Irlande. Quand je le change, je reçois le résultat: Temps de montée: 2015/11/14 02:07:40 azimut: 157: 24: 58.0, qui est encore à une heure de distance de la passe fermée ci-dessus (visible et non visible) du 14 novembre \t-\t 03:43:02 \t 10 ° \t SSW \t 03:45:18 \t 18 ° \t SSE \t 03:47:38 \t 10 ° \t E \t visible, je vais mettre à jour la question, y compris la TLE – kujosHeist