Warum kann ich die gleiche Koordinate nicht mit pyephem abrufen?

Question

Wenn Sie einfach eine FixedObject erstellen und einen Satz von Koordinaten geben und dann fragen sie zurück Sie eine andere Position erhalten:

>>> import ephem 
>>> TestStar = ephem.FixedBody() 
>>> TestStar._ra, TestStar._dec = '12:43:20', '-45:34:12' 
>>> TestStar.compute() 
>>> print TestStar.ra, TestStar.dec 
12:44:15.34 -45:39:46.8 

Ich verstehe jetzt, dass dies, weil, wie dokumentiert, die FixedBody ist standardmäßig In der J2000-Epoche ist die Standard-Beobachterepoche jedoch der Moment, in dem der Beobachter erstellt wurde, und dies scheint der Standard zu sein, wenn Sie keinen Beobachter angeben.

Allerdings, wenn ich versuche, das zu tun:

>>> TestStar4 = ephem.FixedBody() 
>>> TestStar4._ra, TestStar4._dec, TestStar4._epoch = '12:43:20', '-45:34:12', '2000/01/01 12:00:00' 
>>> TestSite2 = ephem.Observer() 
>>> TestSite2.lat, TestSite2.lon, TestSite2.date = 0,0,'2000/01/01 12:00:00' 
>>> TestStar4.compute(TestSite2) 
>>> print TestStar4.ra, TestStar4.dec 
12:43:19.42 -45:33:51.9 

Sie eine fast identische RA, aber ein Dezember, die von 20 Bogensekunden für dieses Beispiel unterscheiden.

Ich versuche speziell, die J2000-Koordinaten einiger Sterne in einem WEBDA-Katalog zu erhalten, die relative Koordinaten für die meisten Sterne liefern.

Zum Beispiel sieht diesen Zufall Cluster: http://www.univie.ac.at/webda/cgi-bin/frame_list.cgi?ic0166

Die „Koordinaten J2000“ hat nur Informationen über 9 Sterne und fast alle Sterne haben Informationen in der „XY-Positionen“ -Link. Das Zentrum und der Maßstab dieser XY-Positionen ist etwas willkürlich, kann aber in der Site gefunden werden.

Aber wenn nicht wissen, warum diese 20 Bogensekunde Unterschied in Koordinaten ist, weiß ich nicht, wenn mein System fehlschlägt.

Answer 1

Ok an dieser Stelle stelle ich mir vor, dass die Diskrepanz auf einen Korrekturfaktor zurückzuführen ist. Ich weiß jetzt, ich will die Astrometrie Geocentric Position verwenden, so:

>>> import ephem 
>>> TestStar = ephem.FixedBody() 
>>> TestStar._ra, TestStar._dec = '12:43:20', '-45:34:12' 
>>> TestStar.compute() 
>>> print TestStar.a_ra, TestStar.a_dec 
12:43:20 -45:34:12 

Einfach genug (hatte einfach nicht, dass ein Teil des Handbuchs verstanden, sorry).

Ich bin immer noch neugierig, welche von allen Korrekturen das hauptsächlich beeinflussen würde, aber ich kann weitermachen, ohne es zu wissen.