JmmcAsproUvPlanFiller < Jmmc/Software

Cette page regroupe les informations pour le développement de la brique de production d'observables interférométriques pour un plan UV donné.

La vidéoconférence du 21 juin a permis de faire un premier contour de l'outil reporté dans les points ci-dessous que l'on complétera au fur et à mesure des développements itératifs.

%TOC%
---++ Descriptif générale

---+++ Entrées
   1 des cubes d'images (modèles rapides ou grilles de modèles). Fichier fits 2D/3D suivant le format des modèles utilisateurs, par ex. issue d'un code de transfert radiatif.
   1 des OIFits d'observation (observables spectro-interférométriques, plan uv, longueurs d'onde, etc), par exemple issu d'une observation réelle.

---+++ Sorties
   1 oifits avec les observables spectro-interférométriques identiques aux observations, mais issus des modèles
   1 chi^2 correspondant (éventuellement des cartes de chi^2 si on varie des paramètres géométriques tels que la distance et/ou l'angle de position de l'image) 

---+++ Traitements
   * calcul des observables (sans modèle d'erreur) en reutilisant la brique !UserModel d'ASPRO2 mais pas l'export OIFITS existant car il faut rester  au plus pres du format OIFITS en entrée (à étudier en profondeur comment adapter le code)

---++ Proto V1
Deux besoins sont retenus pour la première itération

---+++ Besoin 1 : génération des observables interferométriques 

Le logiciel doit produire un fichier OIFits avec les données simulées d'un modèle utilisateur d'après le plan UV fourni. 

---++++ test 1.1
   * jeux de données
      * images d'entrée
         * [[%ATTACHURL%/aeri_rot_1.65um_hr.fits][aeri_rot_1.65um_hr.fits]]: Image 2D (étoile + disque de gaz) obtenue avec un code de transfert radiatif.
         * [[%ATTACHURL%/image_disco_amhra.fits][image_disco_amhra.fits]]: Image 3D (étoile + disque de gaz; 3 longueurs d'onde) obtenue avec un code numérique "rapide".
      * oifits d'entrée
         * [[%ATTACHURL%/Achernar_all_Pionier_oiDataCalib_2012-09-17.fits][Achernar_all_Pionier_oiDataCalib_2012-09-17.fits]]: OIFITS contenant des vraies observation PIONIER d'une étoile Be (étoile + disque de gaz). Description de ces données dans Domiciano de Souza et al. 2014, A&A, 569, A10.

      * oifits de sortie (Référence) avec les mesures simulées (outil externe)
         * TODO rajouter deux fichiers de références Achernar_all_Pionier_oiDataCalib_2012-09-17.2D.ref.fits et Achernar_all_Pionier_oiDataCalib_2012-09-17.3D.ref.fits

   * Critères de validation:
      * Le fichier de sortie du logiciel doit être équivalent au fichier Référence : 
         * même structure (tables, plan UV, temps ...)
         * seules les valeurs !VIS2DATA, (VISAMP, VISPHI) ou !T3AMP, !T3PHI doivent changer et correspondre à la simulation.


---+++ Besoin 2 : fonctionner en ligne de commande
Le logiciel doit pouvoir être exécuté en ligne de commande. Pas de GUI.

---++++ test 2.2
   * TODO : fournir 
      * url du programme
      * syntaxe de la commande

   * Critères de validation:
      * La commande  
         * doit produire un fichier de sortie si les entrées sont valides
         * doit générer une erreur si les entrées sont invalides


---++ Plan de travail / Protos futurs
Lister ci-dessous les futurs besoins qui seront rajoutés au protos suivants après décision commune. 

---+++ Besoin 3 : calcul d'un chi^2 / carte de chi^2
un programme doit pouvoir calculer un chi^2 ()

   Note: OIFitsExplorer intégrera prochainement la possibilité de plotter simultanément obs/simu et rajouter l'info de chi^2 sera une excellente nouveauté :)

---++++ test 3.1
  * TODO fournir le chi^2 ou la carte de chi^2 pour le fichier de référence de sortie du test 1.1

---+++ Besoin 4 : respect du format de fichier standardisé pour IMAGE-OI pour le stockage des observables simulées
Le fichier de sortie doit correspondre au fichier d'entrée complété par les nouvelles colonnes 
| New column in OI VIS table ||| 
| *Label* | *Format* | *Description* | 
| MODEL VISAMP | D(NWAVE ) | Model of the visibility amplitude | 
| MODEL VISAMPERR | D(NWAVE ) | Model of the error in visibility amplitude (optional) | 
| MODEL VISPHI | D(NWAVE ) | Model of the visibility phase in degrees | 
| MODEL VISPHIERR | D(NWAVE ) | Model of the error in visibility phase in degrees (optional) | 
| New column in OI VIS2 table ||| 
| *Label* | *Format* | *Description* | 
| MODEL !VIS2 | D(NWAVE ) | Model of the squared visibility | 
| MODEL !VIS2ERR | D(NWAVE ) | Model of the error in squared visibility (optional) | 
| New columns in OI T3 table ||| 
| *Label* | *Format* | *Description* | 
| MODEL !T3AMP | D(NWAVE ) | Model of the triple-product amplitude | 
| MODEL !T3AMPERR | D(NWAVE ) | Model of the error in triple-product amplitude (optional) | 
| MODEL !T3PHI | D(NWAVE ) | Model of the triple-product phase in degrees | 
| MODEL !T3PHIERR | D(NWAVE ) | Model of the error in triple-product phase in degrees (optional) | 

Voir le document de standardisation en ligne (table 3 en particulier):
   * PDF : https://github.com/emmt/OI-Imaging-JRA/raw/master/doc/interface/OI-Interface.pdf
   * dépot !GitHub : https://github.com/emmt/OI-Imaging-JRA/

---++++ test 4.1
   * reprendre le jeu de fichiers de test 1.1 avec
      * image d'entrée
      * oifits d'entrée
      * oifits de sortie de référence avec les mesures simulées dans les nouvelles colonnes

   * Critères de validation:
      * Le fichier de sortie du logiciel doit être équivalent au fichier de résultat de test 3 : 
         * même structure (tables, plan UV, temps ...)
         * les colonnes !VIS2DATA VISDATA ou T3 doivent correspondre exactement au fichier d'entrée
         * les colonnes !MODEL_VIS2DATA !MODEL_VISDATA ou MODEL_T3 doivent correspondre exactement au fichier de référence

---+++ Besoin 5 : intégration dans le GUI ?

   * est-ce utile dans l'interface graphique ASPRO2 de pouvoir extraire les conditions d'observations (target, interferometre, instrument = lambda + uv) d'un OIFITS et donc de gérer des données OIFITS ou juste le plan UV dans ASPRO2 ?

   * pouvoir sélectionner des modèles (de cartes issues de codes/services externes) directement depuis l'interface d'Aspro (ici l'OIFits serait calculé suivant la config instrumentale interne à Aspro)
      * faire le recensement des codes existants

---++ Notes/remarques pour le plus long terme
   * réutiliser le modèle de bruit d'ASPRO2

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
fits	Achernar_all_Pionier_oiDataCalib_2012-09-17.fits	r1	manage	47.8 K	2017-06-23 - 12:25	ArmandoDomiciano	OIFITS contenant des vraies observation PIONIER d'une étoile Be (étoile + disque de gaz). Description de ces données dans Domiciano de Souza et al. 2014, A&A, 569, A10
fits	aeri_rot_1.65um_hr.fits	r1	manage	4021.9 K	2017-06-23 - 12:22	ArmandoDomiciano	Image 2D (étoile + disque de gaz) obtenue avec un code de transfert radiatif.
fits	image_disco_amhra.fits	r1	manage	1541.2 K	2017-06-23 - 12:23	ArmandoDomiciano	Image 3D (étoile + disque de gaz; 3 longueurs d'onde) obtenue avec un code numérique "rapide"

This topic: Jmmc/Software > WebHome > JmmcASPRO > JmmcAsproUvPlanFiller
Topic revision: r6 - 2017-06-23 - GuillaumeMella