Traitement d'image solaire caca boudin

[Philippoïde 6 716]

je vois que ça planche dur ce soir. Merci à vous

L'image filtrée par FFT est assez intéressante mais il se passe des trucs bizarres sur le fond du ciel qui devient strié. j'avais tenté qq essais de FFT avec "quantum image" (si mes souvenirs sont bons) mais je n'avais pas trouvé le logiciel facile et les résultats probants. Je ne suis pas trop versé dans les transformées de fourier mais mon sentiment est qu'il devrait y avoir des possibilités dans cette direction.

Thierry, l'origine des stries horizontales est bien due à la nature de l'instrument (fente irrégulière, poussières sur la fente, poussières sur la barrette CCD, réponse des photosites, j'en oublie peut-etre).

Celine, le filtrage passe-bas que tu appliques diminue les stries horizontales mais également les détails du disque solaire. Il n'y a pas que l'étape d'élimination des stries. Il faut ensuite "sortir" les détails et améliorer les contrastes. Cette méthode ne peut pas convenir car elle génere un flou très important.

Jean-Philippe, tu est tout pres de la méthode que j'ai adoptée. Je passe effectivement par une sorte de flat synthétique.

Je ne sais pas s'il reste d'autres techniques de traitement que nous n'aurions pas évoquées... si oui, je suis tout ouïe devant l'écran

[Thierry Legault 5 826]

j'ai bien pensé aussi au flat synthétique avec moyenne (ou médiane) de lignes, mais on se heurte au fait que l'image n'est pas uniforme (disque solaire sur fond noir), donc que ça va varier d'une ligne à l'autre, du coup il faut ruser avec les histoires de moyenne. A creuser, mais cela risque de ne pas être simple de conserver l'apparence du disque intacte.

[Philippoïde 6 716]

schématiquement, l'idée est que l'image brute est la superposition de l'image solaire et d'une trame parasite.
L'ensemble des moyennes par lignes peut donc etre considéré comme une fonction de basse frequence (soleil propre) sur laquelle est superposée un signal de haute frequence. En modélisant le profil basse frequence et en le soustrayant du profil recueilli, on a le signal parasite.
ombre au tableau: les taches importantes, filaments denses ou régions actives générent aussi des écarts au modèle.
ça n'est donc pas la panacée d'où mon souhait de peaufinage

Mais maintenant que j'y pense, et pour faire plaisir à V. Cotrez, faudrait que j'essaie de virer les bandes de Jupiter

[Thierry Legault 5 826]

une soustraction ? vu les causes que tu invoques, j'aurais plutôt dit une correction multiplicative type flat, non ?

[Philippoïde 6 716]

Oui, en fait y'a de ça
je suis en train de faire une page pour expliquer la démarche. http://astrosurf.com/spectrohelio/process.htm

[Philippoïde 6 716]

MDRRRR jupi sans ses bandes !!!

[PascalD 4 385]

Bonjour,

Je suis très très loin d' être un gourou dans le traitement d' image (ceux qui ont vu les horreurs que j' ai posté sur ce forum le savent déjà), mais je suis tombé par hasard sur ce papier, qui décrit la technique utilisée pour traiter les mosaiques de Lunar Orbiter ... Si les artefacts de cette image proviennent aussi d' un balayage, c' est peut-être applicable ? Je laisse les spécialistes juger.

Cf: http://www.lpi.usra.edu/meetings/lpsc2002/pdf/1099.pdf

A+
--
Pascal.

[Philippoïde 6 716]

-> PascalD
Merci pour l'info. Il y a une similitude du problème "stries" et je vais faire chauffer mon dico anglais-français pour lire ça plus à fond.

[maire 252]

Sauf que dans ce cas les stries sont plus distantes l'une de l'autre (par rapport à la taille de l'image) et il me semble pour que cela relève plus de l'application d'algorithmes de mosaïcage que d'un traitement fréquentiel du type transformée de Fourier comme cela été très bien indiqué. A priori c'est plus évident d'utiliser les traitements fréquentiels pour ôter une bonne part d'un tramage "resséré" et régulier que pour les images Lunar Orbiter.

------------------
Eric.
http://www.astrosurf.com/voute/vout

[Philippoïde 6 716]

-> maire
j'aurais cru le contraire. La "trame" des images de Lunar Orbiter est assez régulière, donc de fréquence bien définie et donc (je suppose) facilement traitée dans le domaine de fourier

[celinedl 584]

sinon, ya une solution pour enlever ces stries mais d'une manière un peu spéciale, histoire que les stages effectués resservent! Si tu as plusieurs images, on extrait les composantes principales des images (donc on devrait trouver une part des images où ya que le soleil, et d'autres rien qu'avec le bruit) et puis on reconstruit sans ces composantes bruitées. Bref on applique ca : http://cdelleluche.chez.tiscali.fr/gretsi03_455.pdf

Cela fait longtemps que je pense que cela peut être utile pour les imateurs "amateur" ou bien les vidéos, mais je n'ai pas encore eu le temps d'adapter le programme. Mais si tu veux, tu me donne plusieurs images, et j'essaie de faire tourner ca.

[Philippoïde 6 716]

--> Celinedl
j'en suis resté aux ACP classiques, là ça se corse. j'ai peur que les contraintes sur le "cube" ne soient pas respectées :
mes images font 1000x1000 pixels et je ne peux disposer que de quelques images sur un intervalle de temps réduit (pour que la rotation du soleil ne soit pas trop importante)
qu'en penses-tu ?

[celinedl 584]

c'est une ACP en fait, sauf qu'on évite de supprimer trop de défauts en rajoutant les résidus.
Et puis le mot "cube" est mal choisi puisqu'il n'y a pas forcément le même nombre de pixels dans toutes les directions. Si tu as un ensemble d'images pas trop éloignées dans le temps, on peut toujours tenter.

[Philippoïde 6 716]

Combien d'images minimum ??
Si le temps le permet je fais une rafale d'images en fin d'aprem (1 image toutes les 3 ou 4 mn ) et je t'envoie ça ce soir.
les stries constituent un motif quasi-constant lors d'une séance d'observation. je pense qu'il serait bon de changer la position
du soleil sur le CCD pour ne pas avoir de correlation soleil/stries.
si tu vois d'autres choses importantes à prendre en compte...

[celinedl 584]

combien d'image... je ne sais pas, disons une dizaine, ca devrait aller.
Quant à changer la position de la CCD pour ne pas avoir les stries au même endroit, je ne sais pas si cela va arranger les choses. Au contraire, si les stries apparaissent au même endroit, elles apparaitront comme composantes, et peut-être plus faciles à enlever. J'ai un doute là... Peut-être pourrais-tu faire deux blocs d'images, un avec la CCD à la même position, puis en tournant.

Dernière remarque : ca fait quelle taille tes images, que tu ne fasses pas exploser ma boîte mail!

[Philippoïde 6 716]

--> Celine
tu peux telecharger les images ici: http://www.astrosurf.com/spectrohelio/icon/images-spectrohelio.zip
bon amusement.

Sinon, voici Jupi avant et après traitement.

[Thierry Legault 5 826]

un truc pour faire disparaître le globe de Saturne en gardant les anneaux : prendre des images de Saturne dans une raie bien particulière de l'infra-rouge proche, une raie d'absorption du méthane. La lumière est tellement absorbée par le méthane du globe qu'on ne le voit plus, ne restent que les anneaux qui semblent flotter autour de rien (enfin presque rien, il y a quand même la silhouette du globe qui en masque une partie).

[Ce message a été modifié par Thierry Legault (Édité le 18-09-2003).]

[celinedl 584]

ok, je vais essayer de traiter ca, c'est cool l'adsl chez soi!
Premiers résultats demain matin, si j'arrive à me lever (et aussi à me souvenir comment mes programmes marchent!)

[Philippoïde 6 716]

--> Celine,
Y'a pas le feu, prends ton temps. Les images sont brutes et leur nom correspond à l'heure TU d'acquisition. (j'ai du battre mon record de vitesse de prise de vues) Y'a un peu de rotation mais ça devrait pas trop gêner.

--> Tierry,
Je ne peux hélas pas faire ce genre de choses avec mon spectro. Il m'arrive déja de manquer de photons au centre des raies solaires...

[Philippoïde 6 716]

Et j'ai crié, criéééé : "CELINE" pour qu'elle revienne... (et essaie de traiter mes images)

[celinedl 584]

voui voui, je ne suis pas très loin, seulement... ben j'ai un problème que je sais pas d'où il vient, ma matrice de covariance n'est pas définie positive, ce qui donne des valeurs propres négatives, et horreur des composantes principales complexes. Hem, je vais essayer d'y réfléchir, mais bon donne moi un peu de temps, j'ai un article à finir de rédiger, mes billets pour le chili à prendre (surtout à me faire rembourser, ah l'agent comptable de la fac...), toute la biblio sur l'interférométrie à lire, mon dossier de thèse à remplir (ah l'administration), et j'en passe...
Mais promis, je vais retenter d'ici la semaine prochaine. Si tu fais de nouvelles images (quoi qu'il fait tellement moche que ca va attendre...) ca serait bien qu'elles soient toutes centrées. C'est peut-être ca le problème en fait. Allez, un p'tit masque circulaire, et ca devrait aller.

Allez, patience encore un peu

[Philippoïde 6 716]

T'as un planning genre les 12 travaux d'Hercule. "l'agent comptable" !!! ça sera l'épreuve la plus difficile, je compatis.

je ne pourrai probablement plus refaire de série d'images avant avril (j'hiberne).
Es-tu sûre qu'une translation des images change qq chose ?

Phil

[maury 0]

A vue d'oeil, première chose qu'un super traitement ne remplace jamais une bonne acquisition...
A priori c'était un jour de vent ?
Donc j'imagine qu'en fait le disque du soleil, à la place de défiler lentement et correctement partait d'avant en arrière sur la fente d'entrée ?
Le pb est différent si tu dois utiliser des logiciels existant (genre la FFT dans PRISM ou autre chose) ou si tu es partant pour faire un logiciel qui corrigerait ca. A priori la seconde solution me semble meilleure. Je prendrai le problème à l'origine. Dans un premier temps, et parce que c'est un pb unidimensionnel, je ferais une somme de toutes les lignes pour avoir un genre de flat, à savoir si ton image est bien exposée, tu dois avoir une courbe régulière, j'imagine aussi croissante depuis le moment où tu tapes le bord du soleil, passant par un maximum au centre, et symétrique par rapport à ce centre. Si l'image est bougée à cause de la non linéarité du mouvement, tu dois avoir des creux et des bosses sur cette courbe. Il faut arriver à ramener cette courbe perturbée à une courbe lisse normale. A priori ce genre de calcul doit être hyper rapide à faire. Peut être qu'un simple lissage de la courbe (passe bas) suffirait. Ensuite tu agis exactement comme pour un flat, en remplacant la valeur de chaque pixel par le rapport entre la valeur originale sur cette ligne par la valeur lissée, ou un truc du genre. Enfin, que je dis, comme ca.... :-)
Alain

[Philippoïde 6 716]

--> Maury

Merci Alain, t'es arrivé à la meme conclusion que moi (http://astrosurf.com/spectrohelio/process.htm)
mais je souhaite savoir si d'autres approches sont possibles et voir ce que ça donne
Phil