astronomy - J2000 の座標を現在の日付に変換しますか?

Question

私は VSOP2000 データ ファイルを使用しており、時刻 T における太陽と月のデカルト座標を計算することができました...問題は、その座標 (J2000 の時間枠内) を現在の日付 (時刻 T )? アルゴリズムを知りたいだけです。

Answer 1

J2000 は単なる JD2451545.0 ( ref ) であるため、J2000 を JD に変換し、アルゴリズムに従って JD からグレゴリオ暦の日付を取得できます。これは、天文学的アルゴリズムの第7章に基づいた、まさにそれを行うためのPythonコードです。

from datetime import datetime, timedelta

class J2000:
    def __init__(self, value):
        self.value = value

    def to_JD(self):
        return self.value + 2451545.0

    def to_datetime(self):
        jd = self.to_JD()
        # based on AA Ch7
        z = int(jd + 0.5)
        f = jd + 0.5 - z 
        alpha = int((z - 1867216.25) / 36524.25)
        a = z + 1 + alpha - int(alpha / 4)
        b = a + 1524
        c = int((b - 122.1) / 365.25)
        d = int(365.25 * c)
        e = int((b - d) / 30.6001)
        day = b - d - int(30.6001 * e) + f
        d = int(day) # day
        d_frac = day - d # fractional day
        m = e - 1 if e < 14 else e - 13 # month
        y = c - 4716 if m > 2 else c - 4715 # year
        return datetime(y, m, d) + timedelta(days=d_frac)

例えば：

j = J2000(2456647 - 2451545) # J2000 for 12/20/13 12:00:00
j.to_datetime()
=> datetime.datetime(2013, 12, 20, 12, 0)