tech Programme de pretraitement : alignement images(étoiles)

[Laphroaig]

Toute ma chaine de traitement est faite sur des flottants 32bits

 

Je viens faire l'operateur Median

Sur 11 flats avec pour résolution 1280 * 960

Moins de 6 secs ça reste correct.

 

Sur le calcul de la médiane, comme elle est calculée pixel par pixel, on pourrait se contenter de charger en mémoire un seul pixel de chaque image à la fois. Je pense aussi qu'il existe un algorithme itératif permettant de ne charger qu'une image à la fois (ou peut-être deux pour conserver l'histoire de nombre impair).

 

Ce n'est pas nécessaire toutes les images sont en RAM

 

Pour la normalisation tu as raison je m'embrouille ...

Car c'est sur que si j’étale de 0 a 1 ca va pas aller

 

Mais merci je vais regarder la doc Iris

[Laphroaig]

Excusez moi ca fait longtemps ... j

Je reprend ce que dit fred

 

La méthode utilisée par Iris est de calculer l'intensité médiane de chaque image de flat (prétraitée après retrait de l'offset et du dark maître) et de normaliser chaque image à un niveau commun (par exemple 1000). Une fois cette opération faite, le flat maître est obtenu en prenant pour chaque pixel la médiane des pixels correspondants des images individuelles. Si tu calcules ton flat maître en réel et non en nombre entier, tu peux normaliser à 1.

 

J'essaie de resumer pour voir si j'ai compris ...

 

Soit N : le niveau commun choisi (1000 dans ton exemple)

Et une image I1 avec une intensite Mediane M1

 

Est ce que cela veut dire que les pixels égal a M1 dans l'image I1

apres normalisation auront pris la valeur N ?

Ou

Est ce que la valeur N est la valeur Maximale après normalisation ?

Ou

rien compris et aucun des deux

 

Merci d'avance

 

PS : Sur la doc iris rentre pas trop dans le detail

[Laphroaig]

Pour normaliser il me faut une borne min et une borne max de toute façon

 

Donc deux solution possibles

Prenons 1000 comme valeur de normalisation

 

La facon simple : 1000 est le maximum apres normalisation

borneMin = minImage / maxImage * 1000

borneMax = 1000

 

Facon plus complique : 1000 est le niveau Moyen

borneMin = minImage / medianImage * 1000

borneMax = maxImage / medianImage * 1000

 

Dans le second cas si on prend 1 comme valeur de normalisation

On va se retrouver avec des valeurs superieure a 1.

Lors de la division des brutes par le flatMaster ca aura tendance a baisser l'intensite des pixels les plus sensible ? (et vice versa) ...

 

C'est ca l'idée ...

Modifié 14 octobre 2016 par Laphroaig

[Eric S]

Tu as trois flats: flat1 à 950 de moyenne, flat2 à 1000 de moyenne et flat3 à 1125 de moyenne.

 

La moyenne des moyennes est de 1025. Tu modifies les trois flats pour amener chacune de leurs moyennes à 1025.

 

Pour chaque pixel du flat1 : flat1N(x,y) = flat1(x,y) * 1025 / 950

Pour chaque pixel du flat2 : flat2N(x,y) = flat2(x,y) * 1025 / 1000

Pour chaque pixel du flat3 : flat3N(x,y) = flat3(x,y) * 1025 / 1125

 

Ensuite, toujours pixel par pixel, on fait la médiane de flat1N, flat2N et flat3N -> flatMaster.

 

Enfin, pour chaque image à prétraiter, après soustraction de l'offset et du dark, tu normalises:

Pour chaque pixel de l'imageI : imageIPre(x,y) = (imageI(x,y) - offsetM(x,y) - darkM(x,y)) * 1025 / flatMaster(x,y)

[Laphroaig]

Ok ca revient a la methode 2 (plus compliqué )

A part que tu fixes le niveau de normalisation automatiquement (pas forcement bête)

 

Le darkMaster et le biasMaster je reviendrai sur ca apres

 

J'ai codé la methode 2 en fixant moi meme le niveau de normalisation

 

Par Contre quand je fait pretraitre(brute) / flat_Master

Je peux faire ce que je veux a la sortie pour le reequilbrage

 

 

1 er cas

Ou je remonte tout les pixels avec moins de sensibilité

En considerant le max

En clair les max du flat ne changeront pas dans la brute

 

2eme cas

Ou si je prend la median de mon flat

Et dans ce cas la les pixels de la brute qui sont au median de mon flat ne changeront pas de valeur

et ceux en dessous je rehausserai leur intensite

et ceux au dessus je baisserai leur intensité

 

Voila mon petit probleme

 

---------------- Ajouter-------------

OK excuse Eric je viens de comprendre tu es pour le deuxieme cas

Tu prend la moyenne et plus la mediane aussi : c'est du détail mais je peux aller dans ce sens

------------------------------------

 

Autre question

Comment vous faites un bias Master ?

Normaliser les bias ? Median comme flat ?

 

 

Quand tu dit Brute - darkM -offsetM : le darkMaster ne contient pas l'offset ?

 

Après je m'occuperai du dark master sans sigma/clipping

Et encore après si je suis pas degouter du dark master du sigma/clipping

 

Voila si vous avez du temps a perdre

 

Merci Eric

 

PS : je donne acces a mon code sans aucun probleme

Pour ceux qui le souhaite et qui s'y connaisse : Attention les controles sont minimalistes !

Pour les raws je lit que les .fits 8 et 16 bits 2D

les dependances : openCV et compilateur C (ca fonctionne sous windows)

Modifié 14 octobre 2016 par Laphroaig