> Transformation > Mosaïque stellaire

Lancement : Menu [Transformation] > [Mosaïque] > [Stellaire]

Cette fonction n'est pas accessible depuis la version light. Afin de pouvoir l'utiliser, il est nécessaire d'obtenir une clé "version complète" pour le logiciel.

Cette fonction permet de réaliser la mosaïque de 2 images stellaires.

Lorsque le CCD est trop petit pour faire entrer un objet trop grand, la prise de vues multiples en décalant le CCD à chaque fois (en gardant une zone commune à chaque image) permet de compenser cette faiblesse des CCD et des optiques. Grâce à cette commande, cette opération peut être possible.

Attention : dans tous les cas, pour réaliser une mosaïque d'images, il faudra respecter les conditions suivantes :

• Les images doivent avoir le même échantillonnage (si ce n'est pas le cas, utiliser à cette fin la transformation [Transformation] > [Agrandissement - Réduction]).
• Les images doivent avoir la même orientation (si ce n'est pas le cas, utiliser à cette fin la transformation [Transformation] > [Rotation]).
• Les images devront avoir une partie suffisamment commune pour les relier !
• Par contre, les images n'ont pas besoin d'avoir le même niveau radiométrique ! ( c'est-à-dire qu'au départ, les pixels communs n'ont pas besoin d'avoir à peu près la même intensité )

Chargez les 2 images concernées, puis sélectionnez la première. Au lancement de la fonction, la fenêtre suivante s'ouvre :

Sélectionnez alors la seconde image en utilisant les flèches en bas à droite de la fenêtre.

Considérons deux images de M81 (exemple des images de test); l'une représente la partie inférieure et l'autre la partie supérieure. La partie commune se trouve dans la partie inférieure de la galaxie. On représente ci-dessous les deux images en zoom x1/2 :

Il faut repérer au moins une étoile commune au deux images. Il est conseillé de réaliser la mosaïque en zoom x1. On doit alors indiquer une étoile commune à chaque image. On devra d'abord la sélectionner dans la première image soit M81_1 puis dans la deuxième image M81_2. Il sera appliqué la fonction centroïde sur cette étoile dans les 2 images. Ceci permettra de connaitre la position précise des 2 étoiles dans chaque image, le flux et le valeur du fond de ciel. En effet, connaissant ces valeurs, le logiciel va pouvoir organiser l'image finale de la façon suivante :

• Afin que le chevauchement des images n'apparaissent pas, il va falloir transformer la deuxième image de telle façon que le fond du ciel de chaque image soit identique et que la même étoile appartenant aux deux images ait le même flux dans chaque image.
  
  La fonction centroïde permet de mesurer le flux moyen de l'étoile (E1 et E2 dans chaque image) et le bruit de fond de ciel (F1 et F2).
  
  On doit multiplier la deuxième image par un coefficient C égal à (E1 / E2).
  
  Il va falloir ajouter un offset (une constante) à tous les pixels de l'image égal à (F1 - (F0*C)). Le fait de ne pas faire directement F1 - F0 permet de mieux transcrire le fond de ciel en fonction de la luminosité. On peut évidemment ajouter un offset négatif.
• 
  
  
• Il va être créé une nouvelle fenêtre dont le coin en bas à gauche représentera comme d'habitude le pixel (1,1). La fenêtre va avoir une dimension de telle façon à accepter l'image résultante. Soit ( X1, Y1 ) les coordonnées de l'étoile au centième de pixel dans la première image et ( X2, Y2 ) ses coordonnées dans la deuxième image. Soit également ( DIM_X1, DIM_Y1 ) les dimensions de la première image et ( DIM_X2, DIM_Y2 ) les dimensions de la deuxième image.
  On va devoir considérer quatre cas :
  
  
  1. Si ( X1 > X2 ) et ( Y1 > Y2 ), alors la fenêtre mosaïquée aura l'image 1 en bas à gauche et l'image 2 en haut à droite. Le pixel ( 1, 1 ) correspondra au pixel ( 1, 1 ) de l'image 1. Comme le pixel central de l'étoile sélectionné n'aura sûrement pas une valeur " ronde " les valeurs de l'intensité des pixels de l'image 2 seront recalculées par interpolation. La dimension totale de la fenêtre sera :
     
     

     DIM_X = INT ( X1 + ( DIM_X2 - X2 ) )
     DIM_Y = INT ( Y1 + ( DIM_Y2 - Y2 ) )

  2. Si ( X1 > X2 ) et ( Y1 < Y2 ), alors la fenêtre mosaïquée aura l'image 1 en haut à gauche et l'image 2 en bas à droite. Le pixel ( 1, 1 ) ne correspondra à rien ! Son intensité sera nulle. Par contre, le pixel ( 1, DIM_Y ) correspondra au pixel ( 1, DIM_Y1 ) de l'image 1 et le pixel ( DIM_X, 1 ) correspondra au pixel ( DIM_X2, 1 ) de l'image 2. Comme le pixel central de l'étoile sélectionné n'aura sûrement pas une valeur " ronde " les valeurs de l'intensité des pixels de l'image 2 seront recalculées par interpolation. La dimension totale de la fenêtre sera :
     
     

     DIM_X = INT ( X1 + ( DIM_X2 - X2 ) )
     DIM_Y = INT ( ( DIM_Y1 - Y1 ) + Y2 )

  3. Si ( X1 < X2 ) et ( Y1 < Y2 ), les commentaires seront les mêmes qu'en 1. mais il faut inverser les indices 1 et 2.
     La dimension totale de la fenêtre sera :
     
     

     DIM_X = INT ( ( DIM_X1 - X1 ) + X2 )
     DIM_Y = INT ( ( DIM_Y1 - Y1 ) + Y2 )

  4. Si ( X1 < X2 ) et ( Y1 > Y2 ), les commentaires seront les mêmes qu'en 2. mais il faut inverser les indices 1 et 2.
     La dimension totale de la fenêtre sera :
     
     

     DIM_X = INT ( ( DIM_X1 - X1 ) + X2 )
     DIM_Y = INT ( Y1 + ( DIM_Y2 - Y2 ) )

• La partie commune aux deux images doit maintenant avoir le même fond de ciel et les étoiles doivent avoir le même flux. La superposition peut alors être effectuée de plusieurs façons :
  1. en gardant la partie commune de la première image et en supprimant celle de la deuxième image. ( versa-vice ).
  2. vice-versa !
  3. en faisant une moyenne des deux parties communes.
  4. en faisant une interpolation entre les deux parties communes.

La valeur du pixel résultant va dépendre de la valeur des pixels de la zone en commun.

V3 = ( V1 + V2 ) / 2 si ( V2 - V1 ) < Sup ( V1, V2 ) / 10
V3 = Sup ( V2, V1 ) si ( V2 - V1 ) > Sup ( V1, V2 ) / 10
( Sup ( V1, V2 ) renvoie la plus forte valeurs des deux valeurs V1 et V2 )

Le cas de nos deux images est le cas exposé au point 1. Il n'est pas possible de montrer ici le résultat.
Cependant, le programme affiche, avant d'effectuer le compositage, les résultats de ses calculs sous forme de deux coefficients :

Les coefficients ont la signification suivante :

• L'offset est un nombre qui va être ajouté à la deuxième image ( celle qui n'est pas par défaut ). Ce nombre peut être négatif. ( = " Traitements / Ajouter une constante " )
• Le coefficient est un nombre par lequel on multipliera tous les pixels de la deuxième image. ( = " Traitements / Multiplication" )

En cliquant sur "OK", la mosaïque stellaire apparaît. L'image résultante est fournie pour exemple.

Images supportant cette fonction :
images monochromes ou couleurs 16 ou 32 bits.

Fonction(s) Script associée(s):
Aucune

Image de test :
2 images de M81 à raccorder dont le résultat est donné.
2 autres images de M81 résultant de la séparation d'une image de départ. Les paramètres n'ayant pas été modifiés, les coefficients renvoyés doivent être de 1 pour le coefficient multiplicatif et 0 pour l'Offset ...

Conditions d'activation du menu :
Qu'une image soit ouverte.