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
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))