Correction de la turbu par PSF.

[brizhell 2 626]

ms, c'est bien pour cela (utilisation des emccd), que je me suis lancé dans l'aventure.

L'ensemble des paramètres à implanter dans le fichier de configuration de Aida (AIDA_Settings.py) et permet de paramétrer la configuration du calcul.

C'est trés détaillé (même si un peu complexe), dans la publi suivante :
http://code.google.com/p/aida-deconvolution/downloads/detail?name=Hom07_JOSAA.pdf

Les exemples de la figure 5 sont trés parlant !! (notamment les colonnes C et D)

Bernard

[-ms- 2]

La différence entre C et D : c'est très intéressant pour celui qui cherche un algorithme performant pour séparer des étoiles multiples ou qui veut voir plus de détails à la surface des satellites de Jupiter ou qui veut mettre en évidence l'évolution des bandes nuageuses d'Uranus dans l'infrarouge ou qui veut détailler une nébuleuse planétaire ou qui veut résoudre certains amas globulaires ou ...

Avec un simple télescope de 300mm + une EMCCD (voir lien ci-dessous) + un environnement de calcul Priithon + un processeur multi-coeurs et un site pas trop mauvais , il y a matière à faire plein de petites manips dans pleins de domaines différents, c'est en gros ce qu'on appelle "lucky imaging".
http://www.buythis.co.uk/open_product.asp?pid=4894&catid=117&offset=&mid=65

[Ce message a été modifié par ms (Édité le 19-09-2011).]

[brizhell 2 626]

C'est vrai, j'ai commencé a y goûter ;-) :
http://brizhell.org/Manip_Merlin.htm
http://brizhell.org/doubles_et_tavelures.htm
http://brizhell.org/satellites_de_jupiter.htm
http://brizhell.org/pulsar_crabe.htm

Une EMCCD pour moins de 2000 euros ca commence a devenir interessant.

[Nebulium 1 436]

Bon, au boulot

Pour ceux qui veulent faire comme les grands, la FAQ

Confirmation que la méthode AIDA nécessite d'avoir fait des prises pour les PSF, contrairement
à la méthode BiaRam qui peut partir en principe d'un simple petit carré mais qu'il vaut mieux
alimenter au départ avec une PSF crédible (pour l'instant, je trouve FitsWork utile pour cela)
et une image un peu boostée aux ondelettes.

[EDIT] Merci ms pour la modif de l'image[/EDIT]

Pour ceux que ça intéresse, une autre méthode pour publier du code directement copiable depuis le forum :

- Clic droit dans la fenêtre active "Invite de Commandes"
- Sélectionner tout
- Copier dans le PP : Ctrl + C
- Ouvrir "Notepad", coller
- Faire le ménage si besoin
- Publier dans le forum entre balises "CODE" et "/CODE" (entre crochets, comme IMG) comme ci-dessous :

code:

---

D:\>dir
 Le volume dans le lecteur D s'appelle Disque local
 Le numéro de série du volume est 7656-B19E
 Répertoire de D:\
19/09/2011  00:55    <REP>          AIDA
19/09/2011  00:01    <REP>          Priithon_25_win
19/09/2011  11:55    <REP>          Titan
...
D:\>
D:\>cd Titan
D:\Titan>\Priithon_25_win\priithon_script  \AIDA\AIDA.py -S AIDA_Settings.py

---

Rien ne semble se passer, pas de retour.
Quand le calcul est fini :

code:

---

D:\Titan>

---

[Ce message a été modifié par Nebulium (Édité le 19-09-2011).]

[-ms- 2]

brizhell : super intéressantes tes différentes manips. Certaines doivent bien se prêter à AIDA.

Nebulium : les prises pour les PSF ne sont pas obligatoires mais souhaitables sinon PyRAF qui existe aussi sous Windows et Python 2.5 peut être utilisé à cet effet.

Il existe une multitude d'algorithmes de deconvolution +/- rapides, +/- efficaces et +/- robustes. Le peu que j'ai lu sur AIDA m'incite à en apprendre d'avantage à son sujet dans les prochains mois.

[Nebulium 1 436]

ms, sais-tu à quoi sert le fichier .DS_Store dans Titan (et AIDA, d'ailleurs)
Peut-on le virer sans dommage?

Merci pour le tuyau "PyRAF".
Mais s'il est à utiliser pour se faire des PSF sur mesure, FitsWork me suffit.

[jldauvergne 15 098]

ms,intéressante cette caméra tu sais si il y a une version N&B ?

Je ne comprends pas bien quelle est la fréquence d'utilisation ?
Horizontal : 15,734KHz / Vertical : 59.94Hz Horizontal : 15,625KHz / Vertical : 50Hz

Les pixels sont gros, il font 12,5µm, tu sais qui fabrique ça ?

[-ms- 2]

C'est un fichier caché généré par Mac OS X que Windows n'utilise pas en principe.

[brizhell 2 626]

Salut JL,

il s'agit des fréquences de balayage respectivement horizontales, et verticales définies par le mode video :
http://www.cinetson.org/phpBB3/tritubes-f7/frequence-balayage-resolution-taux-rafraichissement-t2638.html

Bernard

[-ms- 2]

C'est un capteur TC246 de Texas avec des pixels de 10µ.
Pour augmenter la sensibilité, tu as un filtre à déposer (camera jour/nuit).
A vérifier si tu peux pas monter directement le capteur monochrome TC247.
C'est vrai que le prix des EMCCD baisse ... lentement mais surement.

[Nebulium 1 436]

Retour à AIDA...

ms, merci j'ai nettoyé les trognons de pomme inutiles

J'ai commencé mes tests AIDA sur l’œuf au plat, en utilisant une PSF "4 tours" de BiaRam.
L'image résultante est molle mais exempte de bruit et montre très mal les ultimes détails, même en partant d'une image préaccentuée en ondelettes.
Et bien sûr il faut faire un 2ème tour d'ondelettes après AIDA
Par expérience, je sais qu'il me faut un peu jouer quelquesjours avec tout ça pour obtenir de meilleures images, je montrerai des résultats comparatifs couleur quand ils seront plus convaincants.
A contrario, je vais aussi appliquer du BiaRam à l'original de Titan

[-ms- 2]

Tu as des exemples de déconvolution et filtrage dans le lien suivant :
http://obswww.unige.ch/~behrend/page_omg.html

Vu le nombre de librairies en langage C existant sur le sujet de la déconvolution, il serait intéressant de faire une interface en Python ou mieux en Ruby (ça s'apprend plus vite et c'est plus lisible) puis de la poster sous github afin que les internautes puissent y accéder et l'améliorer.

Il faut rester dans la compatibilité Windows/Mac/Linux car chacun doit pouvoir utiliser son environnement favori.

Pour la déconvolution AIDA me semble très performant et pour le filtrage IRIS me semble très complémentaire.

[Nebulium 1 436]

Merci pour le lien

C'est un pro de la chose !
Les traitements de Mars et Amalthée sont remarquables, dans le sens qui m'intéresse, celui d'un certain réalisme photographique.
Les Jupi de la fin m'impressionnent moins.
Le logiciel BIFSxx n'est pas encore publié, son auteur cherche des images test, je vais lui en donner quand j'aurai un peu progressé avec AIDA et BiaRam, il y a largement de quoi faire pour moi avec ces deux-là en tripotant les fichiers de paramétrage.
J'ai eu un contact positif avec le Dr Tadrous, il est surchargé en ce moment mais très intéressé par la question.

[-ms- 2]

BiaRam a l'avantage d'être bien documenté et d'être écrit en C. Je suis en congés la semaine prochaine et je vais en profiter pour y jeter un oeil.
J'ai l'impression que ce programme est surtout utilisé dans le domaine de la microscopie aussi ce serait une bonne chose de faire quelques exemples astro (étoiles doubles, planètes, satellites de planètes, ...).

[Nebulium 1 436]

Je suis un peu surpris...
Je suis en train d'essayer la déconvolution myope de BiaRam sur la Titan proposée par AIDA, qui illustre d'ailleurs la publication officielle.

J'ai juste donné un poil d'ondelettes pour mieux amorcer le travail :
.

De gauche à droite :

- L'original de départ
- Le résultat AIDA obtenu avec les réglages de démo et 6 PSF véridiques.
- Un léger coup d'ondelettes pas fignolé en bruit sur l'original, bien sûr .
- La réduction par 40 (environ) d'une mosaïqueCassini très haute résolution (réorientée approximativement) citée dans la publication.

???

[Ce message a été modifié par Nebulium (Édité le 22-09-2011).]

[jldauvergne 15 098]

Pas mal, c'est un jpg au départ ? Tu as essayé de voir ce que sortaient d'autres méthodes seules ?

[Nebulium 1 436]

Salut JL

Laquelle, pas mal ? Les ondelettes ?
L'original Titan du Keck est du FITS 32 bits flottant, BiaRam et AIDA (et moi) travaillons exclusivement dans ce format
Stp, jette un œil à la publication en .pdf, tu auras tous les détails.

[jldauvergne 15 098]

comment ça laquelle ?

[Nebulium 1 436]

quote:

---

Pas mal, c'est un jpg au départ ?

---

Ah oui, tu veux sans doute parler de l'original Titan en AO au Keck ?
Meuh non, c'est aussi un FITS 32 flottant
L'astropro du lien de ms écrit :

"Vous pouvez poster ici toutes sortes d'images à tester avec BIFSxx, même dans des formats moins communs. A éviter absolument: le jpg destructeur d'information en les images immorales! Les pré-traitements non-linéaires (par exemple correction d'histogramme et de gamma, sauf pour les pellicules) sont à éviter; les corrections de zéro, de courant d'obscurité et d'aplanéité du champ sont vivement conseillées pour les images DTC. Les résultats les plus intéressants seront affichés ici."

Pour info, toutes mes copies d'écran sont en .png

En ce moment, je prépare pour la nuit Titan au BiaRam. J'ai un truc bizarre dès le démarrage avec la vraie PSF (moyenne des 6 aussi bien qu'une seule) qui a des bavouilles irrégulières autour, je vais faire avec ma "synthétique" pour la Juju du T1m qui est bien propre, elle.

Tant pis pour la science pour l'instant, mais demain je reteste avec la finale élaborée par AIDA

Elle ressemble à ça :
.

Et je pense que ce sont les 4 points en carré qui font dériver BiaRam

[Ce message a été modifié par Nebulium (Édité le 22-09-2011).]

[jldauvergne 15 098]

Juste pour voir, j'ai passé la titan 8 bits ci dessus à la moulinette de focusmagic, ça sort en gros la même chose. J'ai réglé avec un défocus de 6 pixels. A droite j'ai tenté une version à base de traitement classique (ondelette et passe haut). Ca rejoint un peu le cas de Jupiter, quand l'image est bonne tu vas sortir l'info avec des algo classiques. Pour une psf tordue là par contre c'est plus compliqué.

[Nebulium 1 436]

Ta manip aurait été plus parlante en partant du vrai original et en appliquant une réduction de dynamique de 10% au moins pour provisionner les accentuations ultérieures qui font monter le contraste réel et t'ont fait perdre, je pense, des détails dans les hautes lumières
J'ai donc conduit mes tests d'une manière qui devait assurer des performances croissantes en matière de restauration

De gauche à droite :

- Le résultat AIDA, identique à celui de la planche précédente, comme référence. Pour lui rendre justice, je devrai le remplacer par celui que donneraient sur lui les traitements appliqués aux trois autres résultats issus de BiaRam, à savoir successivement : réduction de dynamique, ondelettes, léger gamma. D'accord, ce n'est plus exactement la science recherchée par l'auteur d'AIDA avec les contraintes que cela impose, donc je ne trahirai pas son travail.
- Les résultats BiaRam, issus de la même solution PSF de départ, celle calculée par AIDA en résultat final, détourée dans un carré de 24x24 pixels mais avec une solution image qui est successivement :
* la sortie AIDA (image de gauche)
* l'original Titan traité par ondelettes avec des taux moitiés de ceux utilisés pour la troisième image de la planche ci-dessus.
* la troisième image de la planche ci-dessus.

Dans les trois cas, la méthode BiaRam a été appliquée à l'image Titan d'origine.
J'ai pu appliquer successivement en final en limite de bruit des taux d'ondelettes croissants donnant plus de finesse et de contraste. En particulier, la dernière image me semble laisser percevoir dans la zone claire au-dessus du centre un léger nid d'abeilles, perceptible également sur l'image Cassini de la planche plus haute. Mais elle n'atteint pas le niveau de détail de ton image de droite. Qu'as-tu fait exactement ?

[Ce message a été modifié par Nebulium (Édité le 22-09-2011).]

[jldauvergne 15 098]

Oui j'ai cramé, c'était juste vite fait pour vérifier une chose. C'est que quand l'image n'a pas de défaut majeur on peut l'attaquer de façon conventionnelle. Il faudrait vraiment que tu fasses un test plutôt un une optique astigmate ou decollimatée.

quote:

---

Qu'as-tu fait exactement ?

---

Je t'ai dit focusmagic en partant d'un flou de 6 pixels. Les ondelettes et le passe haut sur celle de droite je n'ai pas pris de notes sur les coef.

[Nebulium 1 436]

Bon, il faut que je regarde ce focusmagic, connais pas.

M...e, je viens de voir que ma 2ème planche est crâmée sur le forum

Je vais regarder pourquoi, la 1ère ne l'est pas (je l'avais assemblée différemment) et la refaire demain.

[Ce message a été modifié par Nebulium (Édité le 22-09-2011).]

[Nebulium 1 436]

Me voici à la montagne dans l'espoir d'un beau ciel d'automne, espoir déçu (pour l'instant)
J'en profite donc pour mettre au clair les résultats de mes investigations béotiennes (toutes corrections par des lecteurs savants sont bienvenues, merci ) .
Quelques précisions complétant ce qui a été écrit plus haut :
Les traitements ont pour point de départ l'image de Titan (au Keck, en AO, 2.06 microns) publiée dans le kit de démo AIDA en Fits 32 bits flottant.
Il se trouve que l'image initiale publiée est inversée en miroir, toutes les autres le seront donc aussi.
Des détails techniques sont donnés p 1593 du JOSAA, en p 14 du document .pdf
Le 15/01/2005, le diamètre apparent de Titan était de 0.84", il occupe en format natif (celui de toutes les images) quelques 96 pixels, donc le pixel fait moins de 1/100 ", pour fixer les idées.
L'idée de base de ces manips est d'évaluer l'intérêt de divers logiciels de déconvolution comme moyen d'accentuation des images planétaires en haute résolution, à utiliser seuls ou en complément des méthodes courantes comme les masques flous, les ondelettes et autres filtres spatiaux.
JLD dit que la méthode est efficace pour les culs de bouteille et /ou les instruments mal réglés et/ou les prises de vues mal ajustées, mais pas pour contrer la turbulence. Comme il est bien connu que ni les instruments utilisés par les astrams ni leurs méthodes de travail n'entrent dans les catégories évoquées ci-dessus, j'ai commencé avec des prises de Jupiter par JL au T1m (je n'ai pas fini, je vais y revenir) puis continué - excusez du peu, mais l'occasion fait le larron, grâce aux idées inspirées par les contributeurs à ce fil et au coup de main de MS, merci à tous -, avec Titan au Keck en optique adaptative (AO).

Quelques rappels techniques :

Les méthodes de déconvolution dites "aveugles" ou "myopes" procédent par itérations successives (selon divers algorithmes et paramétrages) en partant d'une proposition de PSF plus ou moins arbitraire (le moins c'est le mieux )
et éventuellement d'un modèle de solution qui aide les algos à cheminer vers le résultat optimum le plus plausible.
Il me semble que ce modèle a intérêt à être le plus résolu possible en utilisant les méthodes courantes, mais il faut éviter au maximum le bruit et bien sûr tous les artefacts. Il est inutile que les inconditionnels du jpeg 8 bits aillent perdre leur temps dans cette affaire
Quand on n'a pas de PSF "shootée" sous la main (JL traîne les pieds pour le T1m , voir par ailleurs ses objections parfaitement justifiées), on peut partir de celle rapidement ajustée à la main proposée par FitsWork avec le filtre d'accentuation par déconvolution one shot que l'on peut enregistrer directement.
Aux pinailleurs, je suggère de choisir une "conique", d'en ajuster la forme avec la rubrique "Courbes" de n'importe quel soft d'imagerie astro (attention, il faudra la sauvegarder en FITS 32 flottant) pour la tester rapidement avec la restauration "Non aveugle" de Bia Qim, yaka suivre le tuto super bien fait, juste mettre 32 itérations.

Dans les exemples qui suivent, les images en sortie de déconvolution (attention, elles sortent souvent stretchées) ont subi une suite de traitements identiques :
- Réduction linéaire de dynamique (sinon on est sûr de crâmer les hautes lumières plus tard)
- Accentuation par ondelettes avec réduction du bruit (paramètres variables selon les images)
- Légère correction de gamma.
On peut ergoter sur la permutation des phases 2 et 3

De toute façon, j'ai cru comprendre que les scientifiques ne mangent que de la déconvolution et des opérations linéaires, qui seules garantissent le respect des échelles de luminance.

Nous autres astrams -moi du moins- recherchons plutôt un certain réalisme photographique avec un maximum de détails fins véridiques sans faire mal aux yeux , mais que l'on ne trouve pas toujours à la fin .
.

De gauche à droite, 1 2 3 4 5 6 7

1. L'image originale Titan stretchée
2. La sortie AIDA stretchée
3. La sortie AIDA (3) post-traitée pour une comparaison honnête avec les autres. Je ne suis pas un spécialiste des ondelettes, on doit pouvoir faire mieux pour cette phase
4 5 6. Traitement au BiaRam (paramètres demoBlind).
Le traitement est appliqué à l'original Titan, avec comme propositions de solutions la PSF optimalisée par AIDA pour les trois, et respectivement, la sortie AIDA brute, l'original Titan avec une accentuation "ondelettes" modérée, l'original Titan avec une accentuation "ondelettes" normale
7. Réduction inversée et pivotée de l'image HR en fausses couleurs visible et IR de Cassini PIA 7965.

Plus haut, JL a montré des détails bien plus fins. Alors il faut serrer encore un peu (pas trop, sinon elle charbonne et ressemble à un dessin - encore que moi je n'aie rien contre les dessins, j'aime bien ceux de Biver & co, mais ici on cause photo) pour faire parler l'image, 1 2 3 4 :
.

1. L'image originale Titan stretchée
2. L'image 6 ci-dessus
3. Traitement musclé :
- Modèle de solution, l'image 6 (post-traitée), BiaRam
- Masque flou
- Filtre passe-haut
- Réducteur de bruit
4. La couche R (= IR, mais quelle longueur d'onde?) de l'image Cassini.

Rappel : moins de 1/100" par pixel ...

[Nebulium 1 436]

J'ai repris les essais sur l'assemblage WinJupos. En fait dans ce cas, j'étais parti des deux jeux de prises RGB prétraitées "morphed" par JL, ce qui rend difficile les comparaisons. En tout état de cause, je ne suis pas convaincu pour l'instant par les résultats obtenus avec AIDA pour les détails fins, il me reste à voir cependant s'il n'y a pas un effet utile sur le rapport S/B en prétraitement. Ceux obtenus avec BiaRam sont plus parlants.