DEPOUILLEMENT DES DONNEES DPSM
description
Transcript of DEPOUILLEMENT DES DONNEES DPSM
DEPOUILLEMENT DES DONNEES DPSM
DE THEMIS ET DU PIC DU MIDI
Thierry ROUDIER OMP/UMR5572
Atelier DPSM 12 Décembre 2003, Tarbes
DOUBLE PASSAGE SUR LE RESEAU
Multi-Fentes(Boite à Prismes) 2eme passage 1er passage
SPECTRE DPSM
Chaque canal correspond à un filtre étroit avec une variation de longueur d’onde à l’intérieur dans le sens inverse de la dispersion (boites 9 et 11 canaux , invern=0) .
Dans le cas de THEMIS (2 boites 8 canaux , invern=1), nous avons une variation de longueur d’onde à l’intérieur des canaux dans le sens de la dispersion.
Dispersion
CANAL
aile aile
coeur
OBSERVATION DPSM
SPECTRE DPSM OBSERVE FIELD STOP ( F.S.)
FLAT FIELD ( F.F.) DARK CURRENT ( D.C.)
i
jSTEP_X
STEP_Y
x
y
RESULTAT: 1 PROFIL EN CHAQUE POINT DU CHAMP
exemple sur une très petite portion du champ profil en un pixel du champ
DETERMINATION DES FLUCTUATIONS
D’INTENSITE ET DE VITESSE
Profil moyen
I
Vr= c /
Corde
RESULTAT: CHAMP D’INTENSITE ET DE VITESSES
INTENSITE CONTINUUM
INTENSITE 288 mA
INTENSITE 144 mA
VITESSES DOPPLER 288 mA
VITESSES DOPPLER 144 mA
A DIFFERENTES CORDES
RESULTAT: CHAMP D’INTENSITE , DE VITESSES ET DE CHAMP MAGNETIQUE
INTENSITE CONTINUUM
INTENSITE 144 mA
B LONGITUDINAL
VITESSES DOPPLER 144 mA
TRAITEMENT DES DONNEES DPSM Logiciels de Pierre MEIN
TRAITEMENT D’UNE SEQUENCE TRAITEMENT DE PLUSIEURS SEQUENCES
WIDGET IDL
Interface IDL Jean BONMARTIN
MSDPAUTO
utilise: - séquence.par- t2003.parou- p2003.par
fabrique :- le directory des résultats- ms.par- les fichiers P et Q de sortie- les graphiques
idl>msdpauto,’/util/data/’
idl>ima=readmsdp(‘p*’)
idl>msdp (nécessite idl astro)
LOGICIELS et DOCUMENTATION
LOGICIELS: se trouvent dans par exemple: soft0312
DOCUMENTATION:
GENERALE: - readme.txt guide général de l’utilisateur - auto.txt guide d’utilisation de msdpauto - séquence.txt exemple de données pour msdpauto - param.txt liste des paramètres de ms.par
DIVERSE : - correction.txt liste des paramètres à modifier suivant les cas - captions.txt contrôle des tracés. -filenames.txt description des noms de fichiers aux diverses étapes - remarks.txt quelques exemples et difficultés - journal.txt liste des améliorations des versions successives - signs.txt donne le signe des résultats - widget.txt pour les widgets (non actualisé) - vtt.txt pour les données du VTT
t*fts sequence.par
/data/ /data/auto/
msdpauto
tyyyy.par
N, L, S
/data/auto/dirN_Lkey1
key.par
b*.fts
ms.par
Parameters
Conversion
ms1Computation
Option /no_fort
MSDP DATA PROCESSING
N=sb=seq.
L=cm=line
S=qv=Stokes
ms1
x*L*z*L*
y*L*Sc*L*Sd*L*Sq*L*Sr*L*Sp*L*S
geo.ps g*Lflat.ps f*L*Sgrid.pscmd*.ps quick.ps j*L*Scmr*.psprof.pssq*L*S.pssp*L*S.ps
ms.lisscan.lis
readmsdp
IDL files
Quick-look
Spectroheliog
Averages
Channels
Bisect.
Profiles
Calib.
tl sb sx sy sz cm bs yy mm dd lbd go stx dt sty ny ng nq qv nb bt qp sd
1 3 3 3 3 2 16 03 10 17 0 0 0 60 0 0 4 3 0 1 0 0 2
1 5 5 5 5 0 16 00 08 24 8542 0 5000 60 8500 4 3 1 1 1 0 0 1
1 6 6 6 6 0 16 00 08 24 5890 0 5000 60 11000 3 3 3 3 1 0 0 1
1 8 8 8 8 0 16 00 08 24 4861 0 5000 60 11000 3 3 1 1 1 0 0 1
1 9 9 9 9 0 16 00 08 24 4861 0 5000 60 11000 3 3 1 1 1 0 0 1
end
SEQUENCE.PAR
Télescope
numéro de séquence
obs.
d.c.
f.f.
f.s.
caméra
Nomb. canaux
date
ordre du réseau
Pas en X
t entre scansen 1/10 de s.
polarisation
burst
jusqu’à étape « q » ou « p »
0 = sun1 = dec2=linux
ETAPES DU PROGRAMME MS1
Fichiers Plots
Ss/prog "MSDP" ASCII
xy moyennes des x,y,z x,y,s
geo géometrie des canaux bmg1.f g geo.ps
flat corrections photom. bmc1.f (f) flat.ps,grid.ps
(détermination de ideb) cmf1.f
Calcul de la géométrie des canaux
ETAPES DUPROGRAMME MS1
Calcul du Flat field
Ss/prog « MSDP » ASCII plots
bmc calcul des fichiers "c" bmc1.f c (canaux alignés)
cmd calcul des fichiers "d" cmd1.f d cmd*.ps
quick cartes approchées dme1.f q j (quick.ps)
I,vit (fich.q)
centrage en lbda (fich.j)
corrélations (fich.j)
cmr calcul des "r": cmd1.f r cmr*.ps
options profils
sommes-diff
I,vit
Traitement observations
Profils recentrés
Ss/prog « MSDP » ASCII plots
prof cartes "p" dme1.f p (prof.ps)
profils
I,Q,U,V
sommes-diff
I,v,B//
grayq plots .ps des "q" gray1.f sq*.ps
grayp plots .ps des "p" gray1.f sp*.ps
Bissecteurs,
profils de raie, Vitesses
Tracés
Listes: scan.lis fichier texte avec sorties
ms.lis long fichier texte (déroulement pas a pas)
Fichiers "MSDP": transformes en fichiers IDL par fichierIDL=readmsdp('.......')
Dans ms.par, l'exécution des différentes étapes est commandée par 0 ou 1 sous les étiquettes suivantes: ixy igeo iflat ibmc 1 1 1 1 icmd iquick icmr iprof igrayq igrayp 1 1 1 1 1 1 1
Fichiers ASCII
Exécution des diverses étapes
tel dob nseq nline 1 20031017 3 2
MSDPBMS WAVELNTH GRORDER FSLTH FSWTH STEP_X NBSTEP_X 16 5896 0 60000 300000 5000 20 STEPGRID NBSTGRID GRID_MAX GRID_PER GRID_WID SEQ_STOK BURST 8500 4 0 0 0 3 0
Télescope Date obsNuméro séquence
Numéro caméra
Paramètres non utilisés dans ms.par
ms.parms.par
FILE obs.par
nm lbda dlbd mupris mustep minpro
8 5896 80 3300 800 500
Nombre de canaux/ fenêtre (window)
Distance entre 2 canaux (mAngs.).
Lambda (Angs.)
Pas multifenteboite (micron)
Translation de entre canaux(prismes boite) (micron)
Normalisation du profilvaleur ajustée au centrede la raie
nwinp mgrim nquv ipos burst
2 4 3 4 1
ntmax priscan jypas interc
0 0 5000 15
Nombre de fenêtres (window) / image
Nombre d’état de polarisation
Nombre de positions Y-scan (en polarisation)
Nombre d’imagespar burst
Nombre maximun, pas de grille (en polarisation)
Nombre d’imagespar scan
Ordre des prismespour le champ
Pas en X du balayage (ici 5’’.0) (arcsec/1000)
Distance approximativeentre fin et début de canaux Unité=pixel CCD
win kdecal
2 0
1 50
nbcln nblgn li lj invern
1035 921 133000 9000 1
1035 921 133000 9000 1
Nombre de fenêtres
Entrelacement des canaux
Nombre de pixels de la fenêtre en i
Nombre de pixelsDe la fenêtre en j
Dimension du champ en arcsec en i (*1000)
Dimension du champ en arcsec en j (*1000)
Pour changer l’ordredes canaux
cqp
inveri inverj inverl norma scatter etal
1 0 0 0 0 0
ix1 ix2 jy1 jy2 jyq1 jyq2
0 133000 500 8500 500 8500
0 133000 500 8500 500 8500
Symétrise les cartes / i
Symétrise les cartes / j
Inversion en orientation (lambda)
Normalise intensité(exemple: nuages)
Taux de diffusion(scatter/1000)
non utilisé
Enlève les bordsen x , en arcsec
Enlève les bordsen y , en arcsec
Idem pour les fichiers de sortie« p » et « q »
invi invj istep invers
1 0 8500 0
Inversion des cartes de sortie Inverse les signes des
paramètres de Stokes
Pas en Y (STEP_Y)(arcsec/1000) en polarisation
FILE exe.par
dir
/home/lafon/dpsm/data/dir3_2/
filter b000000_000_000_000000_m0000_00000000.fts
ixy igeo iflat ibmc
1 1 1 1
icmd iquick icmr iprof igrayq igrayp
1 1 0 0 1 0
Directory des fichiers b
Filtre des fichiers b
Différentes étapes: 1 pour utiliser 0 sinon
tob1 tob2
0000000024000000
tdc1 tdc2
0000000024000000
tfs1 tfs2
0000000024000000
tff1 tff2 nff
0000000024000000 1
24000000
24000000
Début et fin de l’observation à traiter
Heures min et max du dark current
Heures min et max du field stop
Heures min et max du flat field Nombre de flat fields utilisés divisé par nqff
tcl1 tcl2
0000000024000000
sundec iswap intert ipermu nqseul milsec
0 1 600 1 0 250
bmg
si sj sgi sgj milang milgeo nleft nright
0 15 15 15 0 3500 0 0
0 15 15 15 0 3500 0 0
Heures pour les calibrations géométriques
Type d’ordinateur
Non utilisé
Swap ou non
Durée minimaleentre 2 scans (1/100) seconde
Echange X et Y
Nombre de couples (si polarisation)
Taille du pixel de sortie ici0.25 arcsec
Seuils en intensité en i et j pour détecter les canaux
Seuil de géométrieécart en régressionen 1/1000 de pixel
Angle des canauxSeuils des gradients d’intensité en i et j pour détecter les canaux
Détermination du bord gauche (droit) d’un canalà partir d’un canal voisin
cmf
inclin milrec calfs caldeb
1 500 0 1
cqp
ideb igri itgri itana jtana calana milalp milzero ijlis
0 12000 33500 16298 0 0 0 0 0
cmd
cented sumd nlisd curvd crecd w1d w2d w3d
1 0 2 0 2000 0 1 0
Type de détection la forme de la raie Type de calcul de la
transmissionrelative des canaux
Seuil pour le recadrage entre le FF et le FS
Translations du séparateur de faisceau de polarisation i et j
Calcul par le programme de la position de la grille (polarisation)
Taille utile en arcsec/1000 des plages de la grille (polarisation)
Période de la grillearcsec/1000
1er point de la première plage utile de la grille arcsec/1000 (position)
Ajustement xy du décalage de l’analyseur (polar. circ.)
Changement d’intensité dusignal avant interpol. (I **a )
Lissage spatial ,bruit
Calcul intensité centre raie
Sortie directe à partir des canaux
Lissage du profil
Correction de courbure par des points voisins.
Correction des cannelures par filtrage Fourier
lmpd lbd1d lbpasd
0 0 0
0 0 0
2 1500 1500
quick
crecq milsigq lcrecq
0 2000 0
cmr
center sumr nlisr curvr crecr w1r w2r w3r
1 0 2 0 2000 0 1 0
lmpr lbd1r lbpasr
7 500 500
0 0 0
0 1500 1500
Réjection pour le calcul des moyennes des valeurs ayant unécart plus grand que sigma *milsigq/1000.
Correction par écarts moyens
Lissage en y
Définitions des paramètres identiques à celles de « cmd »
Spectrohéliogrammes (non utilisé à l’étape cmd car non calibrés).
bissecteurs
Somme et différence (ailes bleue et rouge)
1ère corde : 1.5 * dlbd=1.5 * 80= 120 mA2ème corde : 3.0 * dlbd=3.0 * 80= 240 mA
prof
crecp milsigp lcrecp
0 2000 0
FILE fix.par
reg lin linref iplotg iplotf nqff
0 0 0 2 4 3
npol
1
bmg
(win) i1 i2m j1 j2m lip jeps intvi intvj
1 1 0 1 0 40 20 30 20
2 1 0 1 0 40 20 30 20
Définitions des paramètres identiques à celles de « cmd »
Tracé de géo.ps
Tracé de flat.ps
PARAMETRES FIXES
Définit la succession des paramètres de Stokes pour le flat field
Non utilisé
Non utilisés
Numéro de fenêtre
1er pixel utile et écart au dernier pixel en i et j
Intervalle en i utilisé pour mesurer la courbure des canaux, ici 40%
Intervalle de recherche en j des bords en i (plus grande longueur ) à + ou - jeps pixels
Intervalles en i et j sur lesquels sont calculées les moyennes pour détecter les bords en j et i
(win) leps n1 distor
1 40 1 1
2 40 1 1
bmc
idc dxr100 dyr100
1 0 0
cmf
smoothi smoothj il1p il2p isym
0 0 10 90 0
l
1er canal utile
Prise en compte de la courbure des canaux
Pour le dark current
Numéro de fenêtre
Les corrections dans chaque canal du flat field sont remplacées ou non par des moyennes
Restreint le calcul du profil moyen de la raie spectrale
Le profil est symétrisé
Intervalle de recherche des points de gradients maximun à +/- leps
Petits décalages entre flat field et images du scan
(win) curv iliss jparli lispro
1 1 10 5 10
2 1 10 5 10
(win) jt100 ja100 jb100 jz100
1 0 0 0 0
2 0 0 0 0
cmd/cmr
longw lat absord absorr mps
0 0 1 1 1
Lissage en i avant la détection du centre raie
Lissage parabolique en j avant la détection du centre raie
Lissage du profil moyen utilisé pour calculer les corrections
Numéro de fenêtre
Non utilisés Profil en absorption ou émission pour les fichiers d
Prise en compte de la courbure de la raie
Numéro de fenêtre
Translation en j, en pixel/100, correspondant à la différence de entre 2 canaux
Définit l’inclinaison et la courbure de la raie dans chaque canal
Idem pour fichiers r
Spécifie l’unité de vitesse en m/s
Si 0 paramètres calculés par le programme
quick
lcorq jlap2q icormq copasq milcoq decmq
0 0 0 0 0 0
prof
lcorp jlap2p icormp copasp milcop decmp
0 0 0 0 0 0
gray
igrq jgrq igrp jgrp imax
3 2 4 2 0
Indice du tableau utilisé pour la corrélation spatiale 2D
½ intervalle de superposition entre 2 expositions du scan
taille pour le calcul de la corrélation Pas pour le calcul des différences premières le long de x
Le résultat n’est pas pris en comptesi le maximum de la corrélation 2Dest inférieure milcoq /1000
Non utiliséParamètres identiques à quick
Nombre de tracés en horizontal et vertical fichiers q
Idem fichiers p Nombre maximum de pixels dans la direction y pour l’ensemble des balayages. Permet ainsi d’ajuster l’échelle des graphiques p et q
Contrôle les graphiques
blackq whiteq blackp whitep 0 1 0 1 1 0 0 0 1 2 0 0 0 1 3 0 0 0 1 4 ---------------------------------------- 0 1 0 1 30
end
0 et 1 permet une visualisation identique à celle de TVSCL d’IDL
Si l’on précise les valeurs, on obtient une visualisation avec la
dynamique souhaitée .
0 et 0 pas de visualisation du tableau correspondant
GEO.PS
Gradients d’intensité
A D
BE
Les extrêmes définissent les bords du canal en i
Le programme calcule les droites de régression pour les vecteurs (AD,BE,…) projetés sur i et j
Détermination des positions des canaux
Intensités
FLAT.PSMinimum du signal (centre raie ) + ajustement parabolique
Décalage à même entre 2 canaux successifs (ltrj) Profil moyen des canaux successifs
Profil moyen après correction de la transmission pour la 1ère fenêtre
Contrôle l’entrelacement des canaux pairs et impairs (boites à 16 canaux)
Profil moyenadopté
Début du 1er canal
Idem 2ème fenêtre
Coupe des canaux le long de i
début du dernier canal
Coupe des canaux le long de j
Moyenne des coupes le long de j pour tous les canaux
Idem 2ème boite
RESULTAT: CHAMP D’INTENSITE , DE VITESSES ET DE CHAMP MAGNETIQUE
INTENSITE CONTINUUM
INTENSITE 144 mA
B LONGITUDINAL
VITESSES DOPPLER 144 mA
BON DEPOUILLEMENT