LIB TT
une librairie de prétraitement d'images astronomiques
version 20000215

Manuel de l'utilisateur

Ce guide de l'utilisateur s'adresse aux astronomes qui veulent écrire des scripts de traitement d'image utilisant la librairie TT. Rappelons que seul le format d'images FITS est reconnu.

1. Guide de référence du mode script de TT

L'intérêt principal de la librairie TT est de pouvoir exécuter des commandes entrées sous la forme de script.

Un script TT est composé d'une seule chaine de caractères contenant des lignes (séparateurs \n du langage C). Chaque ligne est analysée séquentiellement. Au sein d'une ligne, la première chaine de caractères rencontrée doit contenir le mot clé de définition. Si le mot clé n'est pas reconnu, le restant de la ligne est interprété comme une simple remarque. Le séparateur blanc est utilisé pour les paramètres suivants de la ligne.

Il existe actuellement trois mots clé de définition :

1.1. Paramétrages des définitions IMA/

La forme générale d'une ligne définie par un mot clé de type IMA/ est la suivante :

définition rep_in nom_in ext_in indice_deb indice_fin rep_out nom_out ext_out indice_out fonction paramètres ...

Exemple d'une ligne de script :

IMA/STACK . i 1 5 .fit c:\toto j . .fit SK bitpix=32 kappa=1.5 nullpixel=-1000 jpegfile

Ce script demande d'effectuer un pile kappa-sigma (fonction SK) des images i1.fit à i5.fit du répertoire courant en une image j.fit dans le répetoire c:\toto. L'image j.fit sera enregistrée en entiers de 4 octets (bitpix=32), le coefficient kappa vaut 1.5 et les pixels de valeur inférieure ou égale à -1000 ADU ne seront pas pris en compte dans calcul des kappa-sigma.

A noter : dans ce propos, le terme de nom complet signifie le nom du fichier incluant le répertoire et le suffixe. Ainsi : c:\toto\j.fit est un nom complet alors que le nom de l'image est j.

1.1.1. Liste générale des fonctions disponibles pour IMA/

mot clé de définition [et paramètres diponibles] mot clé de la fonction [et paramètres diponibles]
IMA/STACK [bitpix jpegfile skylevel nullpixel]  ADD 
SK [kappa] 
SORT [percent] 
MED 
MEAN
IMA/SERIES [bitpix jpegfile skylevel nullpixel] SUB [file offset] 
ADD [file offset] 
OFFSET [offset] 
COPY 
DIV [file] [constant] 
FILTER [threshold] [type_threshold] [kernel_width] [kernel_type] [kernel_coef] 
OPT [dark] [bias] [therm_kappa] [unsmearing] 
TRANS [trans_x] [trans_y] 
STAT [pixelsat_value] [fwhm] [objefile] [pixefile] [border] [detect_kappa]
DELETE 
NORMGAIN [normgain_value] 
NORMOFFSET [normoffset_value] 
CATCHART [path_astromcatalog] [astromcatalog] [catafile] [jpegfile_chart]  [jpegfile_chart2] 
HEADERFITS [file] 
REGISTER [translate] 
ASTROMETRY 
UNSMEARING [unsmearing] 
INVERT [mirror] [flip] [xy]
CONV [kernel_type] [sigma]

 1.1.1.2. Paramètres optionels de la définition IMA/SERIES

1.1.1.3. Paramètres optionels de la définition IMA/STACK

1.2. Fonctions de la définition IMA/STACK

Chacune des fonctions de IMA/STACK peut être associée aux paramètres optionels généraux de IMA/STACK. Les paramètres optionels (et leur argument) spécifiques à chaque fonction sont décrits maintenant :

1.2.1. Fonction ADD

1.2.2. Fonction MEAN

1.2.3. Fonction MED

1.2.4. Fonction SORT

1.2.5. Fonction SK

1.3. Fonctions de la définition IMA/SERIES

Chacune des fonctions de IMA/SERIES peut être associée aux paramètres optionels généraux de IMA/SERIES. Les paramètres optionels (et leur argument) spécifiques à chaque fonction sont décrits maintenant :

1.3.1. Fonction ADD

1.3.2. Fonction SUB

1.3.3. Fonction OFFSET

1.3.4. Fonction COPY

1.3.5. Fonction DIV

1.3.6. Fonction FILTER

1.3.7. Fonction OPT

1.3.8. Fonction TRANS

1.3.9. Fonction STAT

1.3.10. Fonction DELETE

1.3.11. Fonction NORMGAIN

1.3.12. Fonction NORMOFFSET

1.3.13. Fonction CATCHART

1.3.14. Fonction HEADERFITS

1.3.15. Fonction REGISTER

1.3.16. Fonction ASTROMETRY

1.3.17. Fonction UNSMEARING

1.3.18. Fonction INVERT

1.3.19. Fonction CONV

2. Exemples de scripts

Nous allons etudier ici quelques cas classiques de traitement d'images astronomiques.

2.1. Prétraitement des images CCD

Le prétraitement des images consiste à soustraire le dark aux images brutes puis de diviser par le flatfield.
Tout cet exemple traite des images acquises suivantes : Nous allons d'abord préparer des images synthétiques de bais, dark et flat, resultant d'empilements, afin d'atténuer le bruit. Enfin, nous allons prétraiter les images en utilsant la méthode de l'optimisation du thermique.Les images prétraitées seront placées dans le chemin ./prt

2.1.1. Préparation du bias

Le bias est une image de très court temps d'intégration, obtenue dans le noir absolu. Le prétraitement proposé consiste à faire une pile médiane.

IMA/STACK ./bias/ bias- 1 9 .fit ./bias/ bias . .fit MED

2.1.2. Préparation du dark

Le dark est une image de temps d'intégration proche de celui des images sky, obtenue dans le noir absolu. Le prétraitement proposé consiste à faire une pile médiane.

IMA/STACK ./dark/ dark- 1 9 .fit ./dark/ dark . .fit MED

2.1.3. Préparation du flat

Le bias est une image obtenue sur une surface lumineuse la plus uniforme possible. Le prétraitement proposé consiste à faire une correction de dark, suivie d'une normalisation en gain et d'une pile médiane.

SET/VAR $smearing 0
IMA/SERIES ./flat/ flat- 1 9 .fit ./flat/ f . .fit OPT bias=./bias/bias.fit dark=./dark/dark.fit unsmearing=$smearing
IMA/SERIES ./flat/ f- 1 9 .fit ./flat/ f . .fit NORMGAIN normgain_value=10000
IMA/STACK ./flat/ f- 1 9 .fit ./flat/ flat . .fit MED
IMA/SERIES ./flat/ f- 1 9 .fit . . . .fit DELETE

Si les images ont été obtenues avec un obturateur, la valeur de $smearing doit être égale à zéro. Avec des images obtenues sans obturateur, il est nécessaire de placer une valeur non nulle à $smearing.

2.1.4. Prétraitement des images brutes

SET/VAR $sky sky
SET/VAR $smearing 0
IMA/SERIES ./sky/ $sky- 1 9 .fit ./prt/ i . .fit OPT bias=./bias/bias.fit dark=./dark/dark.fit unsmearing=$smearing
IMA/SERIES ./prt/ i- 1 9 .fit ./prt/ i . .fit DIV file=./flat/flat.fit constant=10000
IMA/SERIES ./prt/ i- 1 9 .fit ./prt/ i . .fit NORMOFFSET normoffset_value=200 skylevel

Si les images ont été obtenues avec un obturateur, la valeur de $smearing doit être égale à zéro. Avec des images obtenues sans obturateur, il est nécessaire de placer une valeur non nulle à $smearing.

2.1.5. Régistration et empilement des images prétraitées

SET/VAR $sky sky
IMA/SERIES ./prt/ i- 1 9 .fit ./prt/ i . .fit STAT objefile=./prt/x$sky.fit
IMA/SERIES ./prt/ i- 1 9 .fit ./prt/ i . .fit REGISTER translate=before nullpixel=-1000
IMA/STACK ./prt/ i- 1 9 .fit ./prt/ $sky . .fit MED nullpixel=-1000
IMA/SERIES ./prt/ i- 1 9 .fit . . . .fit DELETE
IMA/SERIES ./prt/ x$sky . . .fit . . . .fit DELETE

2.1.6. Mise en forme de l'image prétraitée

SET/VAR $sky sky
IMA/SERIES ./prt/ $sky . . .fit ./prt/ $sky . .fit STAT fwhm jpegfile=./prt/$sky.jpg

2.2. Calibration astrométrique d'une image prétraitée

IMA/SERIES ./prt/ $sky . . .fit ./prt/ $sky . .fit STAT objefile=./prt/x$sky.fit
IMA/SERIES ./prt/ $sky . . .fit ./prt/ $sky . .fit CATCHART path_astromcatalog=/cdrom/  astromcatalog=Microcat catafile=./prt/c$sky.fit
IMA/SERIES ./prt/ $sky . . .fit ./prt/ $sky . .fit ASTROMETRY
IMA/SERIES ./prt/ $sky . . .fit ./prt/ $sky . .fit CATCHART path_astromcatalog=/cdrom/  astromcatalog=Microcat catafile=./prt/c$sky.fit jpegfile_chart2=./prt/$skyb.jpg
IMA/SERIES ./prt/ x$sky . . .fit . . . .fit DELETE
IMA/SERIES ./prt/ c$sky . . .fit . . . .fit DELETE