Commit 281a2388 authored by Eric Buchlin's avatar Eric Buchlin

Corrections following 2019 session

parent ec059373
......@@ -172,13 +172,13 @@ This gives rise to temporal variations of the position of each line, that come i
**[Q11]** Estimate the variation of the Doppler velocity originating from the motion of the satellite around the Earth. Is this motion sufficient to explain the observed Doppler velocity variations?
**[Q12]** A wavelength correction (corresponding to known instrumental effects) as a function of time is proposed in the Hinode/EIS data files (`d.wave_corr_t` variable in the `.sav` files): can it effectively correct the time variation measured in [Q10]? Use it to produce corrected Doppler velocity maps. Is the average shift now a blue shift or a red shift?
**[Q12]** A wavelength correction (corresponding to known instrumental effects) as a function of time is proposed in the Hinode/EIS data files (`d.d.wave_corr_t[0]` variable in the `.sav` files): can it effectively correct the time variation measured in [Q10]? Use it to produce corrected Doppler velocity maps. Is the average shift now a blue shift or a red shift?
> Indication: `d.wave_corr_t` has to be subtracted from the measured Doppler shifts.
> Indication: `d.wave_corr_t[0]` has to be subtracted from the measured Doppler shifts.
> More corrections are available:
> * `d.wave_corr_tilt`, depending on the position along the slit
> * `d.wave_corr`, for each pixel; this is the sum of `d.wave_corr_t` on the time (or *x*) axis and of `d.wave_corr_tilt` along the slit (*y* axis).
> * `d.wave_corr_tilt[0]`, depending on the position along the slit
> * `d.wave_corr[0]`, for each pixel; this is the sum of `d.wave_corr_t[0]` on the time (or *x*) axis and of `d.wave_corr_tilt[0]` along the slit (*y* axis).
## Determination of coronal temperature
......@@ -201,7 +201,7 @@ These informations can be combined to give the instrument *temperature response*
The `data/AIA/aia_response.sav` file contains the AIA wavelength response (`wvlresp` variable) and temperature response (`tempresp` variable), for each AIA channel, as a function of the wavelength (in ångströms) and temperature (base 10 logarithm, kelvins) respectively.
**[Q13]** Plot the temperature response curves for the different SDO/AIA bands.
**[Q13]** Plot (in logarithmic scale) the temperature response curves for the different SDO/AIA bands.
> To get the useful data in Python, for example for the 211Å channel:
> ```python
......
......@@ -19,7 +19,7 @@ Les chemins ci-dessous sont relatifs au répertoire `/home/ebuchlin/2019tp-plasm
Le répertoire `data/` contient deux sous-répertoires:
* `AIA/` pour des données de l'instrument [AIA](http://aia.lmsal.com/) du satellite américain SDO (lancé en 2010), dans plusieurs bandes (identifiées par un nombre correspondant à une longueur d'onde en ångströms). AIA prend des images du Soleil dans des bandes étroites autour de différentes raies UV émises dans l'atmosphère solaire. Ces données sont de «niveau 1.5», c'est-à-dire que des traitements de calibration de base (dark current, flat field...) ont déjà été appliqués, mais que les intensités sont encore en «data numbers» (DN) au lieu d'être en unités physiques.
* `EIS/` pour des données de l'instrument [EIS](http://msslxr.mssl.ucl.ac.uk:8080/SolarB/Solar-B.jsp) du satellite japonais Hinode, lancé en 2006. EIS est un spectromètre à fente, dans le domaine de l'UV lointain. Le sous-répertoire `L1/` contient les données de «niveau 1» (après dark current, flat field...), pour 5 balayages successifs d'un certain champ de vue.
* `EIS/` pour des données de l'instrument [EIS](http://msslxr.mssl.ucl.ac.uk:8080/SolarB/Solar-B.jsp) du satellite japonais Hinode, lancé en 2006. EIS est un spectromètre à fente, dans le domaine de l'UV extrême. Le sous-répertoire `L1/` contient les données de «niveau 1» (après dark current, flat field...), pour 5 balayages successifs d'un certain champ de vue.
Les fichiers de données sont en format [FITS](https://fits.gsfc.nasa.gov/), un format standardisé par la NASA et très utilisé en astronomie; il comportent une en-tête avec des métadonnées en format texte sur 80 colonnes, puis une ou plusieurs tables de données (avec leurs en-têtes, elles aussi en format texte). Les images en format FITS sont lisibles avec le logiciel [`ds9`](http://ds9.si.edu/site/Home.html).
......@@ -138,7 +138,7 @@ Les raies qui nous intéressent ici sont
> Attention: le deuxième argument de `np.interp()` doit être strictement croissant.
## Déterminations des vitesses Doppler avec EIS
## Détermination des vitesses Doppler avec EIS
**[Q9]** Calculer les cartes de la vitesse Doppler dans les raies Fe XII 195.119Å et Ca XVII 192.82Å, déduites des données spectroscopiques, pour la première observation d'EIS.
......@@ -173,13 +173,13 @@ Cela conduit à une variation temporelle de la position de chaque raie, qui se s
**[Q11]** Estimer l'ordre de grandeur de la vitesse Doppler correspondant au mouvement du satellite autour de la Terre. Est-ce que ce mouvement suffit à expliquer les variations observées?
**[Q12]** Une correction (correspondant aux effets instrumentaux connus) de la longueur d'onde en fonction du temps est proposée dans les fichiers de Hinode/EIS (variable `d.wave_corr_t` dans les fichiers `.sav`): est-ce qu'elle permet effectivement de corriger la variation temporelle mesurée à [Q10]? L'utiliser pour produire des cartes d'effet Doppler corrigées. A-t-on en moyenne plus de décalages vers le rouge ou vers le bleu?
**[Q12]** Une correction (correspondant aux effets instrumentaux connus) de la longueur d'onde en fonction du temps est proposée dans les fichiers de Hinode/EIS (variable `d.wave_corr_t[0]` dans les fichiers `.sav`): est-ce qu'elle permet effectivement de corriger la variation temporelle mesurée à [Q10]? L'utiliser pour produire des cartes d'effet Doppler corrigées. A-t-on en moyenne plus de décalages vers le rouge ou vers le bleu?
> Indication: il faut soustraire `d.wave_corr_t` aux décalages Doppler mesurés.
> Indication: il faut soustraire `d.wave_corr_t[0]` aux décalages Doppler mesurés.
> Pour être plus précis, les corrections suivantes sont aussi disponibles:
> * `d.wave_corr_tilt`, dépendant de la position le long de la fente
> * `d.wave_corr`, pour chaque position de pixel; ceci est la somme de `d.wave_corr_t` sur l'axe temporel (ou des *x*) et de `d.wave_corr_tilt` le long de la fente (axe des *y*).
> * `d.wave_corr_tilt[0]`, dépendant de la position le long de la fente
> * `d.wave_corr[0]`, pour chaque position de pixel; ceci est la somme de `d.wave_corr_t[0]` sur l'axe temporel (ou des *x*) et de `d.wave_corr_tilt[0]` le long de la fente (axe des *y*).
## Détermination de la température de la couronne
......@@ -202,7 +202,7 @@ Ces informations peuvent être combinées pour donner la *réponse en températu
Le fichier `data/AIA/aia_response.sav` contient la réponse en longueur d'onde (variable `wvlresp`) et en température (variable `tempresp`) d'AIA, dans chacun de ses canaux, en fonction respectivement de la longueur d'onde (en ångströms) et de la température (logarithme en base 10, kelvins).
**[Q13]** Tracer les courbes de réponse en température des différentes bandes de SDO/AIA.
**[Q13]** Tracer (en échelle logarithmique) les courbes de réponse en température des différentes bandes de SDO/AIA.
> Pour avoir les données utiles en Python, par exemple pour le canal à 211Å:
> ```python
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment