Question sur les calculs de vitesse radiale...

[artemis 2 941]

Bonjour,

 

Je suis un peu perdu dans les calculs de vitesse radiale, j'utilise 2 méthodes qui semble me donner des valeurs opposées :

 

Sous SpecInti, je lance le calcul en intégrant la correction de la vitesse radiale, j'obtiens une vitesse de 18.512 km/s , corrigée à 14.762 km/s en prenant en compte la vitesse radiale de la cible   (BF Cyg dans ce cas) :

 

Radial velocity correction...
SIMBAD query for equatorial coordinates...
Alpha: 290.973 deg.
Delta: 29.6748 deg.
Topocentric velocity : 18.512 km / s
Object forward the Earth by : 14.762 km / s
Normalize...

 

Comme je ne souhaite pas que les spectres soit corrigés pour les traitements ultérieurs, j'ai intégré les calculs dans mon script python en utilisant pyasl.helcorr de PyAstronomy, la différence de calcul n'est que de 4m/s, mais le signe est inversé, ici -18.516 km/s... il y a quelque chose que je ne comprends pas ? les calculs sont ils différents ?  je suis un peu perdu, est-ce une différence normale entre "Topocentric velocity" et "Barycentric correction" ? :

 

----- HELCORR.PRO – DEBUG INFO – START ----

(obs_long (East positive),obs_lat,obs_alt) Observatory coordinates [deg,m]: 0.5824 47.4439 91 (ra,dec) Object coordinates (for epoch 2000.0) [deg]: 290.972979 29.674803 (ra,dec) Object coordinates (precessed to epoch 2023.863470) [deg]: 291.20804913658674 29.722609396027135 (ut) Universal time (middle of exposure) [hrs]: 18.983760003000498 (jd) Julian date (middle of exposure) (JD): 2460229.29099 (hjd) Heliocentric Julian date (middle of exposure) (HJD): 2460229.291748289 (gmst) Greenwich mean sidereal time [hrs]: 20.324372178446993 (lmst) Local mean sidereal time [hrs]: 20.36319884511363 (dlat) Latitude correction [deg]: -690.4188145719818 (lat) Geocentric latitude of observer [deg]: 47.25211699595223 (r) Distance of observer from center of earth [m]: 6366763.2084071925 (v) Rotational velocity of earth at the position of the observer [km/s]: 0.46427174940306476 (vdiurnal) Projected earth rotation and earth-moon revolution [km/s]: -0.0673193776533686 (vbar) Barycentric velocity [km/s]: -18.448777917979832 (vhel) Heliocentric velocity [km/s]: -18.458289319138164 (corr) Vdiurnal+vbar [km/s]: -18.5160972956332

----- HELCORR.PRO – DEBUG INFO – END -----

Barycentric correction [km/s]: -18.5160972956332

Heliocentric Julian day: 2460229.291748289

 

les coordonnées de l'observatoire et de la cible sont bien sur identiques pour les 2 calculs....

Si une âme charitable peut m'éclairer un peu... je l'en remercie par avance

Xavier

[OlivierG 782]

Tout dépend comment est interprétée cette vitesse., une histoire de convention entre le signe + ou -

moi j’utilise la même convention que SIMBAD :

positif , la cible s’éloigne de nous

négatif, elle se raproche

 

et donc faudrait voir dans la doc de py.astro

aprés la diiference entre vitesse heliocentrique et barycentrique n’est à prendre en compte que pour la mesure de vitesse de l’ordre de 100m/s 

[artemis 2 941]

Merci pour ton retour Olivier,

Ok, je vais vérifier sur la doc de py.Astronomy, mais je prendrai la même convention que toi (et Christian) en me basant sur les données Simbad.

 

Xavier