CORRECTION DE L'EXTINCTION ATMOSPHÉRIQUE


Introduction :

La correction de l'extinction atmosphérique se fait en deux étapes :

  1. Création d'une réponse instrumentale intrinsèque (ri_intrins) à l'intrumentation ;
  2. Application de la réponse instrumentale intrinsèque avec création du profil de la correction de l'extinction atmosphérique qui sont appliqués d'une traite au spectre de l'objet.

Ces deux fonctions ne sont pas encore intégrées aux pipelines et sont donc à executer "à la main" après ces derniers.

Le calcul de la ri intrinsèque (étoile de référence Castor) :

Il faut d'abord linéariser la loi de calibration en longueur d'onde du profil niveau 1b (spectre de Castor), c'est-à-dire sans correction instrumentale :
spc_linearcal castor--profil-1b
Cela créé le fichier castor--profil-1b_linear.

A l'heure de prise de vue d'un observatoire d'altitude 0.2 km=200 m, Castor est à une hauteur de 75° et dans des conditions d'humidité (transparance) de plaine appelées "normal", id. AOD=.14.
Mais d'autres valeurs de cet argument sont possibles : sec/normal/orageux/lourd/valeur_AOD
La synthaxe de la commande est :
spc_calriintrins spectre_étoile_référence_mesurée-1b nom_spectre_bibliothèque_référence { altitude_obserfvatoire(km) hauteur_zénitale(°) humidité(sec/normal/orageux/lourd/valeur_AOD) }
Je lance la commande :
spc_calriintrins castor--profil-1b_linear a2v { .2 75.0 normal }


Le programme sort d'abord la courbe de transmission atmosphérique :

transm_castor
Courbe de transmission atmosphérique pour Castor (hauteur 75°)
On sort également la ri intrinsèque "ri_intrinseque.fit" représentée ci-dessous :
ri_intrins


La réponse instrumentale intrinsèque obtenue.

L'application de la correction de transmission atmosphérique :

Il faut d'abord linéariser la loi de calibration en longueur d'onde des profils niveau 1b :
spc_linearcal NOM--profil-1b
Cela créé le fichier NOM--profil-1b_linear.

1. Sur Rigel :

A l'heure de prise de vue, Rigel (spectre de Rigel) est à une hauteur de 30° et dans les mêmes conditions que précédemment.
La syntaxe de la commande dédiée est :
spc_corrriintrins spectre_astre_étudié-1b ri_intrinseque { altitude_obserfvatoire(km) hauteur_zénitale(°) humidité(sec/normal/orageux/lourd/valeur_AOD) }
Nous lançons donc la commande :
spc_corrriintrins rigel--profil-1b_linear ri_intrinseque {.2 30. normal }

On voit d'abord apparaître la courbe de transmission atmosphérique :

transm_rigelTransmission atmosphérique pour Rigel


Vient ensuite le profil 1c obtenu (bleu) représenté et comparé au spectre sans correction atmosphérique (vert) :


rigel-1c

Le profil 1c de Rigel : ca pourrait bien être celui d'une étoile B81a...

2. Sur Capella :

A l'heure de prise de vue, Capella (spectre de Capella) est à une hauteur de 51°. Mais le spectre de l'étoile indique qu'il ne s'agit pas de Capella.
Nous lançons néanmoins la commande :
spc_corrriintrins capella--profil-1b_linear ri_intrinseque {.2 51. normal }

On voit d'abord apparaître la courbe de transmission atmosphérique :

transm_capella
Transmission atmosphérique pour Capella

Le profil 1c crée (bleu) et comparé au spectre sans correction atmosphérique (vert). Cette courbe est bien plus en accord avec le type spectral G3III :


capella-1c


Profil 1c de Capella

3. Sur Bételgeuse :

A l'heure de prise de vue, Bételgeuse (spectre de Bételgeuse) est à une hauteur de 54°.
Nous lançons donc la commande :
spc_corrriintrins betelgeuse--profil-1b_linear ri_intrinseque {.2 54. normal }

On voit d'abord apparaître la courbe de transmission atmosphérique :

transm_betelg
Transmission atmosphérique pour Bételgeuse

Le profil 1c crée (bleu) et comparé au spectre sans correction atmosphérique (vert). Cette courbe est bien plus en accord avec le type spectral M1-M2Ia-Iab :


betelg-1c
Profil 1c de Bételgeuse