.compile  ~/aobprocs/decon.pro
.compile  ~/aobprocs/deconv_full_lm.pro
.compile   ~/aobprocs/deconv_full_lm.pro
.compile   ~/aobprocs/parameter.pro
.compile  ~/aobprocs/sigmaXY.pro

decon,'128_537810_nonsaturated_avg_ext.fits','PSF_Gl586A_Brg.fits','gl586A_',[62.761110,64.056057, 64.194928 , 66.923322] ,'Brg_DS'

sigmaXY_2,'Gl586A0_Brg_DS.dat','Gl586A1_Brg_DS.dat','Gl586A2_Brg_DS.dat','Gl586A3_Brg_DS.dat','Gl586A4_Brg_DS.dat'

Le resultat du fit est ecrit dans un fichier result.dat
 
 
 
 
 
 
 
 
 
 
 
 

data=readfits('GLiese.fits')
.r  ~/aobprocs/deconv_full_lm.pro
.r  ~/aobprocs/deconv_full_lm.pro

psfa=readfits('PSF_dph.fits')

;Eclate la psf en la centrant en (0,0)
psf=eclat(psfa)

;nombre d'etoiles
n=2 ;ou (3, 4 etc...)

; Positions des centres approximatifs
; pos=[X_primary, Y_primary, X_secondary, Y_secondary].
; Si la primaire et la secondaire sont inversee dans ce tableau, il y aura
; alors une incertitude de 90 degre sur l'angle
pos=[56,65,72,64]  ;(Xp,Yp,Xs,Ys,Xt,Yt,......)

deconv_full_lm, data, psf, a, n, pos_init=pos, NIT=45, /auto_stop, /fixpsf

La variable a contient le resultat du fit.

Cas d'une binaire :

• On peut alors utiliser la procedure ~/aobprocs/parameter.pro pour mettre a  dans un format correct et l'ecrire dans un fichier texte.

 

 
 
 
 
 

parameter, 'Gl570B' , 'Brg', a
 
 

• Ensuite on peut utiliser ~/aobprocs/sigmaXY.pro pour calculer les barres d'erreurs a partir des fichiers de parameter.txt

 

 
 
 
 
 

sigmaXY_2,'Gl487_0_BrG.dat','Gl487_1_BrG.dat'