-
Postów
1 300 -
Dołączył
-
Ostatnia wizyta
-
Wygrane w rankingu
1
Osiągnięcia ANowak
1,6 tys.
Reputacja
Aktualizacja statusu
Zobacz wszystkie aktualizacje ANowak
-
Jest jakaś formuła, która pozwoli mi na przekonwertowanie koordynatów "na dzień" do J2000? chodzi mi o matematyczne rozwiązanie, a nie kalkulatory
-
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)
To kiedyś używałem. Może się przyda.
- Pokaż następne komentarze 3 więcej
-