Salut Etienne,
Il apparait que ton pofil de zet Tau n'est pas bien calibré comme le montre la superposition des raies telluriques (CTRL-V sous Spcaudace) :
D'ailleurs, la RI n'est pas bonne au vu du bord rouge divergent de ton spectre.
Apres calibration sur les raies telluriques, cela va mieux :
réalisée avec :
Code: Select all
spc_calibretelluric zet_tau_110236012.fit
Ecart moyen=-0.000248113636445 A
# zet_tau_110236012-ocal.fit
Ensuite on peut mesurer la vitesse radiale sans oublier de corriger de la vitesse héliocentrique.
Ces deux etapes sont réalisées d'une traite avec spc_vradialecorraccur (menu Astrophysique/Vitesse radiale par ajustement). Autant se faciliter la vie !
Code: Select all
# Usage: spc_vradialecorraccur profil_raies_étalonné type_raie (e/a) lambda_begin lambda_end lambda_réf ?RA_d RA_m RA_s DEC_h DEC_m DEC_s? ?JJ MM AAAA?
Ici avec les limites de la raie d'absorption centrale à 6562-6565, l'ajustement gaussien de spc_vradialecorraccur donne :
Code: Select all
spc_vradialecorraccur zet_tau_110236012-ocal.fit a 6562 6565 6562.82 05 37 38.685 +21 08 33.16
# La gausienne ajustée est définie par :
centre=6563.89369097, FWHM=1.46068304032, Intensité=0.7615773
# (Vdoppler=49.0467 km/s, Vhelio=-25.0000454827 km/s)
# La vitesse radiale de l'objet au 2457497.3574999999 est :
# Vrad=24.0467 +/- 0.4735 km/s
Résultat semblable au tien.
La raie Ha de zet Tau est sujette à moulte effets physiques, donc cela ne m'étonne pas de ne pas retrouver la vitesse radiale affichée sur Simbad.
--------------
Correction de la vitess héliocentrique :
Remarque : la correction de vitesse héliocentrique n'est pas un simple décalage Delta_lambda car ce dernier dépend de la longueur d'onde d'après la loi de Doppler-Fizeau :
Delta_Lambda=Lambda*Vhelio/c
C'est donc la loi de calibration est modifiée, dispersion incluse.
Cela s'effectue avec la commande spc_corrvhelio (menu Calibration/Correction de la vitesse héliocentrique) où l'on précise les coordonnées RA DEC de l'objet (prises sur Simbad, coordonnées FK5) :
Code: Select all
spc_corrvhelio zet_tau_110236012-ocal.fit 6562.82 05 37 38.685 +21 08 33.16
# La vitesse héliocentrique pour l'objet 05h37m38.685s ; +21d08m33.16s à la date du 2016 4 18 vaut :
-25.0000454827 km/s +/- 1.894 km/s <-> -0.547 A
# Spectre décalé de -0.547281274483 A à 6562.82 sauvé sous zet_tau_110236012-ocal_vhel
# zet_tau_110236012-ocal_vhel
Le calcul prend en compte la vitesse de déplacement de la Terre autour du Soeleil à la date de l'observation, mais aussi la vitesse de rotation propre de la Terre projetée à la hauteur de la cible observée, d'où la nécessité de donner les coordonnées RA DEC de l'objet.
----------------
Pour la mesure du rapport V/R, on utilise spc_vrmes (menu Astrophysique/Calculer le rapport V/R) :
Code: Select all
# Usage: spc_vrmes nom_profil_raies lambda_raie_Violet lambda_raie_Rouge largeur_raie ?pourcent_RMS_rejet (150)?
Cette fonction mesure l'intensité des pics comme le montre la figure suivante :
Ce qui donne pour une largeur de raie de 4 A :
Code: Select all
spc_vrmes zet_tau_110236012-ocal.fit 6561 6566.5 4
# Raie V de centre 6561.17494279 et d'intensité 1.73352446547.
# Raie R de centre 6566.51736434 et d'intensité 1.74247498671.
# V/R=0.994863328709.
Si l'on vise les pics de largeur 1 A, on trouve :
Code: Select all
spc_vrmes zet_tau_110236012-ocal.fit 6561.5 6565.9 1
# Raie V de centre 6561.51797733 et d'intensité 1.75777024031.
# Raie R de centre 6565.94475131 et d'intensité 1.75576945208.
# V/R=1.00113955066.
En effet, visuellement, les pics V et R ont à peu près la même intensité.
Mesures réalisées avec Spcaudace v4.4 béta.
Benji
. Voici Zet Tau un peu plus zoomée :
