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.

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.
2. 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
2. 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
    • aeri_rot_1.65um_hr.fits: Image 2D (étoile + disque de gaz) obtenue avec un code de transfert radiatif.
    • 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
    • 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