Julian Day Conversion (based on JavaScript Alghoritm):
DD
MM
YYYY
H
M
S.sss
Julian Day Conversion (based on Python Script and astropy library):
Date & Time: Data Format: (yyyy-mm-dd hh:mm:ss)
JD:
Coordinates Conversion
Set of Coordinates A.R. e Decl.
Ore
Minuti
Sec.ddddd
Ore.decimali
A.R.
Gradi
Primi
Sec.ddddd
Gradi.decimali
Decl.
Inclinazione Prospettiva
Value:
Il Punto g
è uno dei due punti equinoziali in cui il piano contenente l'equatore celeste interseca il piano dell'eclittica.
Tale punto è quello in cui si trova il Sole all’equinozio di primavera, ossia quando esso descrivendo l’eclittica, passa dall’emisfero australe a quello boreale.
Il punto di coordinate equatoriali (0,0), indica proprio il punto g.
Spesso è indicato anche come punto d'Ariete perché in corrispondenza dell'equinozio di primavera di circa 2100 anni fa,
il Sole si trovava nella costellazione dell'Ariete.
Per la precessione degli equinozi, oggi non è più così e in corrispondenza dell'equinozio di primavera il Sole si trova nella costellazione dei Pesci.
Le coordinate astronomiche equatoriali cambiano nel tempo e vanno riferite a un’epoca.
Un’epoca è un istante. L’ultima epoca di riferimento sono le ore 12 del 1° gennaio 2000 (Epoca J2000.0)
Per avere le coordinate all'epoca della data attuale occorre fare le dovute correzioni a causa della precessione degli equinozi.
Epoch Coordinates Correction
Set of Coordinates A.R. e Decl. J2000.0
Ore
Minuti
Sec.ddddd
Ore.decimali
A.R.
Ascensione Retta (J2000.0):
Gradi
Primi
Sec.ddddd
Gradi.decimali
Decl.
Declinazione (J2000.0):
Epoch Coordinates Correction
Set of Coordinates A.R. e Decl. JNow
Ore
Minuti
Sec.ddddd
Ore.decimali
A.R.
0.00
Gradi
Primi
Sec.ddddd
Gradi.decimali
Decl.
0.00
- numpy
- astropy
from datetime import datetime
from astropy.coordinates import FK5, SkyCoord
from astropy.time import Time
import astropy.units as u
import js
from js import document
def mio_julian_day(date):
t = Time(date)
return t.jd
def calculate_jd_python(*args, **kwargs):
date = Element("data_for_python").value
format = "%Y-%m-%d %H:%M:%S"
res = True
try:
res = bool(datetime.strptime(date, format))
print("La Data è: " +str(date))
t = Time(date)
jd = mio_julian_day(date)
print("Mio jd_python: " +str(jd))
document.getElementById("jd_python_result").value=jd
except ValueError:
res = False
document.getElementById("jd_python_result").value="Incorrect data format"
#output_jd_python.write(Element(element_id="jd_python"))
#return t.jd
#date = '2020-10-04 12:15:03'
#date=datetime.now()
ra1 = float(Element("leggo_mia_J2k_AR").value)
dec1 = float(Element("leggo_mia_J2k_Decl").value)
time_now = Time(datetime.utcnow(), scale='utc')
coord_j2000 = SkyCoord(ra1*u.deg, dec1*u.deg, frame=FK5)
fk5_now = FK5(equinox=Time(time_now.jd, format="jd", scale="utc"))
coord_now = coord_j2000.transform_to(fk5_now)
#print("RA in JNOW : {:.8f} degree" .format(coord_now.ra.value))
#print("Dec in JNOW: {:.8f} degree" .format(coord_now.dec.value))
def calculate_coord(*args, **kwargs):
ra = float(Element("leggo_mia_J2k_AR").value)
dec = float(Element("leggo_mia_J2k_Decl").value)
time_now = Time(datetime.utcnow(), scale='utc')
coord_j2000 = SkyCoord(ra * u.deg, dec * u.deg, frame=FK5)
fk5_now = FK5(equinox=Time(time_now.jd, format="jd", scale="utc"))
coord_now = coord_j2000.transform_to(fk5_now)
output_ar = Element(element_id="JNow__ar_1")
output_ar.write("{:.8f}".format(coord_now.ra.value))
output_dec = Element(element_id="JNow__decl_1")
output_dec.write("{:.8f}".format(coord_now.dec.value))