Traitement par le Maximum d'Entropie


Cet exemple illustre comment utiliser la déconvolution par le maximum d'entropie sur une image planétaire ne contenant pas d'étoile.

Les prises de vue ont été réalisées avec un télescope Newton de 150mm de diamètre et 900mm de focale et une barlow 5x.
 

Préparation de l'image

Cette image est le résultat des traitements décrit précédement :
- Alignement et somme.
- Ondelettes
Il faut que l'image soit carrée et sa dimension un multiple de 2 avant 
d'utiliser ces fonctions. Le plus simple est d'utiliser padding  mais il faut peut-être recentrer l'image avant.
Le niveau du fond du ciel doit également être ramener à zéro.

>load sat
>padding 256 256
>bg
Background: 10
>offset -10
>save sat256

 

Ajout d'une étoile de référence

Comme aucun satellite ou étoile n'apparait sur notre image de Saturne il faut l'ajouter à partir d'une autre image. 
Ici on utilise un satellite de Jupiter sur une image prise quelques minutes plus tard dans les mêmes conditions.

Il suffit d'encadrer le satellite avec la souris, de cliquer avec le bouton droit sur l'image et de choisir "Window".

Cette image doit être de la même dimension que celle de Saturne et le niveau du fond du ciel doit également être ramener à zéro.. On utilise à nouveau padding pour cela.
 
>bg
Background: 6
>offset -6
>padding 256 256
>save psf

L'image du satellite est enregistrée sous le nom "psf.fits" elle peut être utlisée pour toute les images prises dans les mêmes conditions de combinaison optique, mise au point et  turbulence.

Le satellite se trouve dans le coin inférieur gauche.
 

Il n'y a plus qu'à ajouter cette image à celle de Saturne et enregistrer le résultat afin de pouvoir faire plusieurs essais.
 
>load sat256
>add psf
>save satref
Maximum d'entropie

Il suffit maintenant de marquer la position du satellite avec la souris et d'utiliser la commande "mem" . Le paramètre correspond au nombre d'itérations et il faut le déterminer avec quelques essais en le variant de 5 à 20 et en surveillant l'apparition d'effet artificiel ou de bruit.
 

>load satref
>mem 8
>save satmem
Voila le résultat de cette commande.

Vous pouvez également essayer le traitement avec la méthode de Richardson-Lucy (commande rl au lieu de mem) mais j'ai trouvé qu'elle donne beaucoup plus de bruit sur cette image.

Il n'y a plus qu'à decouper l'image avec la fonction "Window" pour enlever notre "étoile".


 


Page précédente | Page suivante | Page d'accueil


© Patrick CHEVALLEY/Pierre DUBREUIL - Août 2000