jldauvergne

Correction de la turbu par PSF.

Messages recommandés

Salut,

>Sur le DVD de PhilTo, il n'y a "que" des R G B IR, avec deux caméras de >tailles de photosites différentes et une PSF monochrome (trefoil)
>Il faudrait vérifier avec Phil quand il a réajusté son secondaire.

Rhaaa, voila pourquoi je pédale dans la semoule sur l'une des trois couches !! Bon ca remet un peu les pendules à l'heure, mais ca veut dire que je ne pourrait a terme que présenter des résultats AIDA que sur la couche R. Je serai effectivmement preneur d'un jeu homogène, y compris avec la PSF cohérente à la prise de vue, mais je suis pas pressé, suis sur autre chose pour le moment.

Sinon les traitement avec la PSF Jupée commence à être trés interessant. Tu peut préciser comment le Dr Tadrous s'y prend ?

@ms

>C'est fou cette volonté de moyenner la turbulence alors qu'il est plus >logique d'appliquer une PSF réaliste à chaque image brute (en ayant >vérifié avant la condition d'isoplanétisme).


Je fait une fixette sur l'isoplanétisme car je travaille en ce moment dessus. ;-)

Je suis plutot d'accord avec toi sur le fond, mais sauf cas un peu isolé (je pense entre autre a ton traitement), as-t-on vraiment le choix ? Peu de monde ont la possibilité d'acquérir une psf cohérente à chaque image brute, en même temps que l'image de l'objet cible. Si le compositage est constitué des quelques meilleures images (typiquement quelques dizaines a quelques centaines selon le type de turbulence) triées selon un critère type max d'entropie ou autre, on peut raisonnablement estimer (ca se calcule, il me semble), que la psf est moins diluées (donc moins étalée à l'intégration) que dans le cas d'un compositage standart a plusieurs millier d'images. De plus, cette hypothèse permet aussi de considérer sur un tel lot d'image la psf comme quasi constante, pré-requis obligatoire à la déconvolution classique. L'approche myope permet ensuite de substituer les fluctuations de PSF de chaques images isolément.

Allez, j'ose :
Dr Bernard Trégon (bon d'accord, je l'ai pas faite en traitement d'image ma thèse, même si s'etait en sciences, mais ca me fait beaucoup rire de voir Dr dans une signature). :-)
Il me semble que c'est de tradition germanique, il est d'ou ce Dr Tadrous ? ;-)
ok je le ferai plus ;-)

Partager ce message


Lien à poster
Partager sur d’autres sites
En 40ms, tu peux traiter la PSF et la reconstruction de l'image par une déconvolution classique. Si tu tardes (256x40ms), l'algorithme devient long à traiter et compliqué.
L'autre problème c'est que le dispositif doit être suffisamment sensible pour enregistrer une image en 40ms sinon il faut compositer et là on entre dans le deuxième cas de figure.
Moralité, il faut tendre vers l'EMCCD.

Partager ce message


Lien à poster
Partager sur d’autres sites
Bonjour

=> ms

Les RGB1.avi du Pic 2010 sont là :

R1
G1
B1

Le serveur de Cédric est à ressources limitées, il ne peut servir trop de monde à la fois, si par hasard ça coince, mieux vaut y revenir plus tard.
J'ai un autre jeu RGB pris à quelques minutes d'intervalles. Si utile je le mettrai en ligne.

ATTENTION : ces films sont sous copyright : S2P/IMCCE/OPM/Jean-Luc Dauvergne.
Toute utilisation doit le mentionner. Pour un usage commercial ou similaire,contacter JLD

=> Bernard

Le Dr Tadrous est british, j'ai cru comprendre qu'il était docteur en médecine, spécialisé en imagerie médicale (et scientifique).
Professionnellement, j'avais constaté que les anglo-saxons et germaniques usent systématiquement du titre "docteur" vis à vis de tous ceux qui sont titulaires d'une thèse.

Concernant les traitements DeconBMI en cours, je suis parti strictement de sa DemoBlind en modifiant sur ses conseils quelques paramètres de contraintes dans les fichiers de réglage, puis la psf de départ.

A l'heure actuelle, je suis un peu déçu du résultat obtenu à partir de l'image simplement empilée.
Il me semble n'avoir rien gagné en résolution de détails fins par rapport aux traitements habituels beaucoup plus rapides. Seul le rendu général pourrait être plus agréable (= plus photographique), il faudra juger sur une RGB finie.
Et puis peut-être qu'il n'y a plus rien à résoudre par rapport au diamètre du scope.
Le filtre Caron de Tria a sorti de fines structures réticulées dont la réalité est à démontrer. Ce n'est peut-être qu'un résidu du bosselage de l'image par la turbulence.

Avec DeconBMI, j'ai encore 1 cartouche et demi à tirer :
Prendre comme modèle de départ une image déjà modérément accentuée (JLD en avait parlé, il faudrait aussi voir avec quelle méthode) plutôt que l'image elle-même et réduire la hauteur de la jupe de la PSF (ce qui diminuerait les suroscillations du départ).

Il me faudrait voir aussi à améliorer l'image empilée.
Je vous laisse regarder la déconvolution unitaire.
Je vais essayer de me familiariser avec le morphing. Le seul outil que j'ai est le Ninox de Wesley, il plantait sur les images 1k x 1k (comme la dérotation, d'ailleurs). Je le lui ai signalé, je n'ai pas encore essayé la version qu'il a sorti depuis.
A propos de dérotation, la fine PSF finale de DeconBMI montre une microélongation équatoriale.
Ce serait donc aussi à refaire avec une dérotée empilée, beurk, quel jargon ! Mais je pense que l'on se comprend ainsi.

[Ce message a été modifié par Nebulium (Édité le 20-01-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites
Merci pour les fichiers ... pas de problème de copyright, c'est pour ma consommation perso. Je vais les traiter durant mes congés ... dans une semaine.

Par curiosité, quel est le rapport (D/r0) correspondant ?

Partager ce message


Lien à poster
Partager sur d’autres sites
quote:
Par curiosité, quel est le rapport (D/r0) correspondant ?

ms, ??? Le r0 du Pic, tu le connais, tu l'as indiqué plus haut
Le D ? je l'ai écrit aussi plus haut !

La note sur le copyright, c'est pour prévenir ceux qui auraient loupé les épisodes précédents du feuilleton et s'aviseraient de tirer un profit quelconque des images finies.

Voudras-tu le 2ème jeu de films, pour faire un limbe propre tout autour, sans passer par des moyens que la morale réprouve ?

Partager ce message


Lien à poster
Partager sur d’autres sites
J'ai récupéré le premier jeu sans problème (15 min par couche environ), tu peux mettre le second.

Pour le r0 du Pic, j'avais donné une valeur au pif ... de mémoire ça peut atteindre 400mm là haut.

Partager ce message


Lien à poster
Partager sur d’autres sites
Je viens de retrouver dans un dossier "Driz" ???un jeu non horodaté, en format d'image moitié du 1er.
Je crois que c'était pour l'empiler avec Astrostakkert, ça a laissé .des traces évidentes dans ce dossier
Je l'upl quand même et le mettrai en ligne demain.

J'ai aussi récupéré au vol en vrac quelques .ser 2011, JL avait eu des soucis à l'upl, je n'ai encore rien trié, mais il n'y a apparemment pas de jeux voisins en temps. Par contre il doit y avoir des Ha et des CH4.

Il faudrait voir avec lui s'il a pu ressusciter son HD .


[Ce message a été modifié par Nebulium (Édité le 20-01-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites
Bonjour

Avec l'aide patiente et soutenue du Dr PJT, j'ai avancé dans mes essais de déconvolutions aveugle/myope.
J'ai repris les tests sur la grosse vieille Juju du Pic en septembre 2010, mais en recommençant par la couche rouge, la meilleure des trois à divers points de vue.
Je livre ici mes premières conclusions, je peux détailler certains points si ça intéresse quelqu'un. Ces essais prennent beaucoup de temps, mais on se fait la main et prend de l'expérience. Là je sature un peu, je vais travailler sur Mars si j'obtiens de bons films et finirai cette Juju en tâche de fond.

1. Pas de miracles à attendre . Le gain de netteté est de l'ordre de celui obtenu avec les méthodes usuelles, on peut éventuellement tabler sur un rendu plus naturel, plus photographique et avec des détails plus authentiques, mais il faudrait avoir des éléments de comparaison réalistes, comme un bon APN sur une sonde
2. L'image de départ doit être aussi naturelle que possible (spatialement invariante).
La qualité du résultat est fortement liée au bruit, aux effets de bords de toute nature (il faut biseauter et flouter les bords d'image, éviter les bourrages inhomogènes comme élargir un fond de ciel à 5 avec des 0), aux artefacts d'empilement (morphing mal contrôlé) et limbes dénaturés pour cause de rotation de la planète ou artefacts de dérotation.
3. La PSF de départ doit occuper tout le terrain, c'est à dire être de même dimension que l'image à traiter (en passant, il est recommandé de traiter des images carrées de 2^n de côté pour gagner fortement dans les temps de calcul toujours prohibitifs) et sans valeurs négatives ni zéros. Donc à rembourrer avec précaution par rapport au fond de ciel pour les disques planétaires.
Elle doit être normalisée (intégrale =1 en 32 flottant)
Une base intéressante est d'utiliser la fonction d'autocorrélation de l'image de départ à laquelle on ajoute quelques millipoils d'une PSF gaussienne qui va à peu près bien avec l'image du genre de celle élaborée par FitsWork par exemple.
Cette PSF de départ est importante, le suivi du graph d'erreur (complété par la lecture du log) permet d'apprécier si les proportions des deux composantes sont correctes.
4. Il y a quelques paramètres à régler impérativement, les contraintes de flux de l'image et de la PSF.
5. Ensuite, une multitude de combinaisons de stratégies : nombre de cycles d'affinage de la PSF, nombre de cycles d'affinage de l'image, nombre d'itérations, rebouclage sur l'image initiale ou affinée, etc...
Le pilotage avec le graph d'erreur peut réserver des surprise, voir plus loin.
6. L'essai illustré ci-dessous porte sur un disque entier 1k x 1k, empilement RGX6 par défaut, 256 meilleures images retenues.
J'ai recherché un compromis entre une netteté raisonnable et un limbe discrètement artefacté. Ce n'est pas forcément la meilleure solution. Le prochain essai à faire est de déterminer la PSF optimale sur un échantillon intérieur de 256x256 soigneusement biseauté en moulinant longuement, puis d'utiliser cette PSF pour une déconvolution classique (avec encore un tas de paramètres à fignoler !), le limbe sera comme il sera. S'il ne plaît pas, travailler les dérotations (de toute façon en principe bénéfiques car permettant de shooter plus longtemps donc de réduire le bruit si la turbulence coopère par son absence). A tester, la transformation cylindrique de deux séquences d'images.


Voici donc, en extraits pour la publication, à gauche, l'image empilée de départ (disponible ici ) à droite celle obtenue avec la PSF déterminée par 64 itérations en 42 h de calcul bicore 1.67 GHz :

La PSF de départ, en log et en 25%


On notera la large jupe. Si l'on dose avec plus d'autocorrélation, le nombre de volants augmente, la suroscillation du limbe aussi.

La PSF finale en log et en 400% détourée

Le graph d'erreur :

J'ai travaillé avec diverses proportions pour la PSF de départ.
Ce graph évoluait à la perfection jusqu'à l'itération 55 où il est arrivé au voisinage de zéro, puis montée en chandelle.
Mais pas de cata sur les images. L'évolution visuelle contrôlée à la fft 2D reste régulière, avec un compromis netteté/ suroscillation du limbe à choisir entre 55 et 64.

Voici les 2Dfft des images ci-dessus, l'originale en dessous, la traitée en dessus.

On constate qu'il n'y a pas de déplacement sensible du coude vers la droite, donc pas de réel gain en détails fins qui pourraient être révélés par une accentuation classique, c'est un peu ma déception pour l'instant.
La bosse à 0.2 résulte de ce f...u faux limbe, donc il faudrait continuer dans les voies évoquées plus haut.
A noter que les courbes 2Dfft ne préjugent pas totalement de la qualité visuelle des images.
Le rabotage brutal de la bosse est destructeur, en jouer avec précaution !


PS : Je suis retourné voir Seddara.
Le principe est d'utiliser pour la génération de la PSF une image externe de la qualité désirée.
Limité par mes faibles capacités matheuses, je pense que l'algorithme doit essayer d'égaliser les fonctions d'autocorrélation.
Je me souviens de la raison pour laquelle je n'avais pas persisté :
La version d'évaluation est tellement bridée que tout test astro sérieux est impossible, surtout à cause de cette f...ue grille.
Contrairement au dr PJT, l'auteur a voulu d'emblée faire du fric, je me demande combien il a trouvé de clients.

[Ce message a été modifié par Nebulium (Édité le 10-02-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites
SeDDaRa est une méthode non itérative c'est à mon avis la méthode idéale pour traiter mes images monochromes en temps réel (type a1).

De plus, cet algorithme a été amélioré récemment (Fast SeDDaRa).

Voir par exemple cette implémentation en langage Python :
http://www.koders.com/python/fid09B7BD9B00C684E3C737208098E198E976D8644A.aspx

Comme indiqué précédemment, j'utilise SeDDaRa pour la reconstruction en temps réel et AIDA pour la reconstruction en batch.

L'autre point c'est d'adapter le mode de saisie des données à SeDDaRa tout en respectant l'hypothèse d'isoplanétisme. Pour AIDA la tolérance est plus importante et le traitement d'images composités semble possible mais c'est très gourmand en temps de calcul.

Partager ce message


Lien à poster
Partager sur d’autres sites
Merci pour le lien

Le problème pour moi, c'est qu'à la rigueur je sais adapter quelques trucs dans un fichier de lignes de commande .bat existant, j'ai réussi à faire tourner AIDA grâceà tes conseils, mais le Python pour le Fast Seddara ???

Je suis volontiers preneur d'un script où indiquer l'image de départ, l'image de référence (il faudrait essayer avec des Cassini ou Hubble ou mieux, de vraies photos avec la résolution qui va bien) et quelques paramètres de réglage. S'il peut fonctionner avec la ressource Priithon_25_win que j'ai installée pour AIDA, ce sera parfait.

Si certains veulent jouer, l'image de départ est en ligne dans mon message précédent.
et voici l'image finale après 64 itérations.

Le Seddara que j'ai regardé est celui de Quarktet inclus dans le Tria payant.


J'ai quand même fait avec Fitswork la manip illustrée par cette copie d'écran 2 :
.
1 2 3
4 5 6
1 = image de départ, 4 en dessous = sa 2Dfft
2 = Image deconBMI obtenue après 42 h de moulinage, 5 en dessous = sa 2Dfft
3 = Preview de 1 avec une demande de correction (en rouge) de sa 2Dfft par recopie manuelle (pas trop fignolée, je n'ai pas pu utiliser la mécanique de sauvegarde dans FW des paramètres de correction qui me semble marcher en différentiel depuis la courbe blanche) de la 2Dfft 4 de l'image 2.
Le calcul est instantané.

Après validation, copie d'écran 3, pour afficher en 6 la 2D fft remesurée sur 3 pour contrôle final.

Il me semble que le principe de Seddara, c'est ça ?
La différence visuelle entre 2 et 3 provient essentiellement de la recopie manuelle de la courbe.
Il faudrait surtout vérifier le rendu des détails, le bruit, etc...


[Ce message a été modifié par Nebulium (Édité le 11-02-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites
J'ai oublié le fichier pdf suivant dans mon mail précédent (ne pas s'arrêter sur la formulation indigeste mais voir le principe à partir des exemples fournis) :
http://144.206.159.178/FT/CONF/16434234/16434379.pdf

quote:
Firstly, a pre-process was added to remove the noise and improve the contrast of the image.

Secondly, an optimum method to extract the filter function which responded to PSF was developed.

Thirdly, an efficient method which was used to estimation deconvolution filter’s SNR was developed.


Très pragmatiques nos amis chinois qui préfèrent optimiser les algorithmes non-itératifs existants pour pouvoir travailler en temps réel plutôt que d'aller chercher la petite bête dans les algorithmes itératifs.

Pour les tests, le langage Python est très pratique (voir l'exemple fourni) mais pour l'implémentation il faudra choisir un langage plus rapide comme Assembleur ou Fortran ou C ou Java si on souhaite utiliser cet algorithme en temps réel.

C'est ce que j'ai prévu de faire en 2012 (avec Java) mais comme je n'y travaille pas à 100% de mon temps ça va me prendre un peu de temps.

A l'arrivée, j'aimerais obtenir un système peu coûteux qui me permette de photographier le ciel comme si j' étais installé sur la Lune avec une simple webcam.

Toutes les conditions semblent remplies :
a) des cameras 1 Mega Pixel (CCD Sony ou CMOS e2V) à moins de 500 euros suffisamment rapides et sensibles pour se lancer dans cette manip qui peut marcher aussi bien pour le planétaire que pour le ciel profond,
b) des ordinateurs portables (Intel Core i5 ou i7) à moins de 1000 euros suffisamment rapides pour traiter les trames en temps réel,
c) des télescopes (C11 ou C14) qui peuvent être poussés suivant le site (rapport D/r0) dans leurs derniers retranchements.

Enfin, tout cela va à contre-courant des temps de poses importants, de l'autoguidage, des montures équatoriales coûteuses, des gros diamètres sous-utilisés, de la recherche désespérante de la nuit magique, ...

[Ce message a été modifié par ms (Édité le 12-02-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites
Bonjour,

Je me suis essayé sur le FIT du pic posté plus haut, avec Astra Image 3.0 Pro.
J'ai utilisé le process Lucy Ridchardson Deconvolution avec un point spread en Gaussian, et un PSF size a 1.5 avec 12 iterations aggressives:


Partager ce message


Lien à poster
Partager sur d’autres sites
"A l'arrivée, j'aimerais obtenir un système peu coûteux qui me permette de photographier le ciel comme si j' étais installé sur la Lune avec une simple webcam."

Le pied!! Je comprends un peu les manips que vous faites et je vous souhaite bon courage mais j'attends avec impatience que tout ça aboutisse!! Y-a un fond de soutient où l'on puisse cotiser??

Partager ce message


Lien à poster
Partager sur d’autres sites
Bonsoir

=> ms
Plutôt que du vrai temps réel, je me contenterais de temps raisonnable pour traiter les films bruts, de l'ordre de ceux de RGX & cie pour empiler.

Je reviens plus tard pour la technique.

=> Paracetamol

Merci Je mettrai en ligne les G et B si tu veux finir l'image pour la comparer aux divers essais faits sur AS (d'autres sur CN) dans les gros fils du Pic de JLD fin 2010, et aussi éventuellement les .avi

[Edit] J'ai jeté un oeil à la courbe fft2D de ton image, elle est très propre.
Du coup, j'aimerais avoir ta psf R en fit ou tif ou png et éventuellement les G et B
. http://thetys.dyndns.org/neb/Jupic_2010/jupR1kh.tif
. http://thetys.dyndns.org/neb/Jupic_2010/jupG1kh.tif
. http://thetys.dyndns.org/neb/Jupic_2010/jupB1kh.tif
.
[/Edit]

=> Pascal

Merci pour ton soutien
Tu peux déjà nous faire des films d'essai.
Et aussi nous dire quand Mars, Saturne ou Jupiter passeront près d'une étoile (G2V de préférence ) avec la magnitude qui va bien pour l'exposition au F/D planétaire optimal du moment, et selon le diamètre du scope.

En passant, tenez, allez jeter un œil ici.

[Ce message a été modifié par Nebulium (Édité le 14-02-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites
quote:
Y-a un fond de soutient où l'on puisse cotiser??

Il faudrait surtout remercier Haorui Zuo car son document (An Efficiency Restoration Method for Turbulence-degraded Image Base on Improved SeDDaRA Method) apporte une réponse simple et efficace au problème de la PSF courte pose.

Partager ce message


Lien à poster
Partager sur d’autres sites
Oui, JL, et les Mars ne sont pas dégueu non plus.
Il est bon, Asimov !
L'Australie est le paradis du seeing
Mais pour l'instant , Bird n'a pas la forme

Et Damian est en silence radio, il fait froid en GB !

J'ai fini par me cogner le pdf chinois que j'avais sous le coude depuis longtemps en survolant cependant les équations qui me dépassent .
La méthode me paraît essentiellement orientée "temps réel " pour des applications de détection d'image (satellites de reconnaissance ?), qui ne correspondent pas tout à fait aux besoins de l'astram lambda.


La qualité d'impression des images d'essai ne permet pas d'apprécier pleinement les résultats.

a b
c d

a = original
b = Seddara amélioré (objet de l'article)
c = Seddara classique
d = variante d'implémentation simplifiée de Seddara amélioré ?

Le texte dit (et les courbes fft montrent) que (après correction d'une erreur typographique) :
b n'est pas mieux que c (tout en étant très proche) mais est meilleur que d.

Finalement, quelques voies à considérer me semblent être les suivantes :


- Seddara amélioré chinois , implémentation (gratuit pour qui sait faire, pas mon cas):
. http://www.koders.com/python/fid09B7BD9B00C684E3C737208098E198E976D8644A.aspx

- Seddara Quarktet (payant)
Ce Seddara semble différent du classique dans la mesure où il utilise une image de référence "étrangère"
Dans la suite TRIA, figure un filtre Caron (déconvolution mono-coup ? dont les effets sont intéressants (DEMO MODE, milieu p 4)
L'origine de la structure réticulée aléatoire reste à éclaircir.

- déconvolutions lentes myope/ aveugle classiques : AIDA, BiaQIm (gratuits)

- à explorer, par un compatriote (gratuit) http://www.astromatic.net/2010/11/08/psfex-release

- à lire : http://www.pixon.com/publications/Puetter_05_ARAA.pdf
. http://www.astroscu.unam.mx/rmaa/RMxAC..25/PDF/RMxAC..25_ofors1.pdf

et pour les bouffeurs d'équations : http://www.ceremade.dauphine.fr/~rivoirar/Spherical.pdf

[Ce message a été modifié par Nebulium (Édité le 14-02-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites
Merci Param, c'est du bon boulot on attend l'avis de JLD.
A comparer par exemple avec celles de Sweiller, Emil, et bien d'autres autres et bien sûr JL,je lui laisse le soin de désigner sa préférée .

J'ai retrouvé le fil-fleuve sur CN avec d'autres traitements étrangers.

Dans le temps j'avais sorti celle-ci, avec deux jeux d'avi, Winjupos et juste un peu d'ondelettes, volontairement un peu claire pour ressembler à l'impression visuelle au T 15":

En matière de détails ultimes, il y a peut-être à préciser un filament dans le petit ovale blanc au-dessus de la GTR, voir vers le bas de la p 1 de ce fil, mais alors le bruit !

[Ce message a été modifié par Nebulium (Édité le 14-02-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites
quote:
J'ai fini par me cogner le pdf chinois que j'avais sous le coude depuis longtemps en survolant cependant les équations qui me dépassent .
La méthode me paraît essentiellement orientée "temps réel " pour des applications de détection d'image, qui ne correspondent pas tout à fait aux besoins de l'astram lambda.

La PSF "courte pose" est constante à l'intérieur de l'angle d'isoplanétisme (angle ne dépendant que du diamètre de l'instrument dans le cas de petits télescopes).
La PSF "courte pose" comprend la PSF "optique" (composante statique) et la PSF "atmosphérique" (composante dynamique).
SeDDaRa s'applique sans problème à la PSF "courte pose".
Dans le cas de la PSF "longue pose" (que 99,9% des astrams affectionnent ) SeDDaRa converge comme tous les autres algorithmes vers une solution "gaussienne" ayant perdu en route une partie des détails et des contrastes présents dans les images brutes initiales.
Une image brute de Jupiter (en gros 10 à 100 fois supérieure à l'angle d'isoplanétisme suivant le site) ne pourra être traitée correctement que par un algorithme "temps réel" appliqué à chaque zone où la PSF reste constante. C'est le prix à payer pour pouvoir utiliser un algorithme de déconvolution pour traiter le problème de la turbulence atmosphérique avec un petit instrument. Avec un gros instrument hors optique adaptative ou mise en orbite cela doit être possible avec les ordinateurs de Météo France ou d'EDF.
L'implémentation "temps réel" proposée par Haorui Zuo est applicable localement à notre problème (c'est une piste pour la réconstruction locale). Après, il faudra l'intégrer à l'ensemble de l'image brute et pour cela, il existe différentes techniques de maillages comme dans la méthode des éléments finis. Enfin, la composition des images brutes pourra se faire et nous pourrons comparer (dans quelques mois) avec les résultats obtenus par la PSF "longue pose".

Des logiciels comme AviStack permettent déjà de traiter le tip/tilt et la distorsion d'une sous-image perturbée, cela se fait en batch avec :
- de très bons résultats pour le Soleil,
- de bons résultats pour la Lune,
- des résultats moyens pour Jupiter,
- pour le reste ???

Il est donc intéressant de se poser la question de savoir si la reconstruction d'une image brute :
- en temps réel,
- complètement corrigée (tip/tilt, distorsion, defocus),
- de n'importe quel objet(stellaire, planétaire et ciel profond),
est possible avec un petit instrument (C11 ou C14) et un ordinateur assez puissant (Intel Core i5 ou i7) ?
La réponse dans quelques mois.

[Ce message a été modifié par ms (Édité le 14-02-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites
ms, si j'ai bien compris, c'est qu'il faudrait mosaïquer les images unitaires des gros scopes pour traiter chaque morceau avec sa psf locale ?
Ce que je voulais dire à propos du temps réel pour l'astronomie c'est qu'on peut toujours traiter après coup les enregistrements unitaires tranquillement avec une méthode plus précise tant que les temps restent acceptables.

Pour cette Juju qui me semble à la limite de résolution du T 1m (au vu des multiples traitements présentés sur AS et CN, il me semble qu'il n'y a plus grand chose à gagner, il faudrait des milliers d'excellentes prises pour réduire le bruit et il faudrait gérer convenablement la rotation) 42 h par couche, c'est limite et je ne parle pas des n x 42 h d'essais préliminaires.

Par contre, je vais essayer de retravailler sur des échantillons plus petits, quand j'aurai un peu joué avec les dérotations de Winjupos.

D'ailleurs, l'examen de la PSF ci-dessus élaborée par DeconBMI est intéressant, je donne mon interprétation pour critique :

On y voit un centre très pointu, 1 pixel, montrant que l'image est en échantillonnage limite, associé à un écho diagonal (bougé de rotation sur 30 s, 256 images retenues en qualité sur quelques 450) et une faible jupe relativement large qui pourrait correspondre aux effets de la turbulence intégrée sur 30 secondes.
Mais au vu du graph d'erreur et de la 2Dfft, il faudrait procéder à beaucoup plus d'itérations encore.

[Ce message a été modifié par Nebulium (Édité le 14-02-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites
quote:
Pour cette Juju qui me semble à la limite de résolution du T 1m (au vu des multiples traitements présentés sur AS et CN, il me semble qu'il n'y a plus grand chose à gagner, il faudrait des milliers d'excellentes prises pour réduire le bruit et il faudrait gérer convenablement la rotation) 42 h par couche, c'est limite et je ne parle pas des n x 42 h d'essais préliminaires.

C'est normal, vous utilisez tous la même approche PSF "longue pose" donc vous convergerez toujours vers la même qualité d'image.
En résumé, chaque image brute de Jupiter peut être maillée en fonction de l'angle d'isoplanétisme. Pour chaque sous-image, une PSF locale doit être obtenue. Chaque maille corrige : les translations (tip/tilt) et les distorsions induites par la turbulence atmosphérique ainsi que les défauts de l'optique. La reconstruction de l'image brute correspond à l'assemblage des différentes mailles. L'image finale est la somme des images brutes corrigées des effet de la turbulence et de l'optique.

Partager ce message


Lien à poster
Partager sur d’autres sites
ms, je te remets les liens vers les .avi :
.
. http://thetys.dyndns.org/Jupic_2010/jupB_2010-09-29_22-59-28.avi
.
. http://thetys.dyndns.org/Jupic_2010/jupG_2010-09-29_22-58-53.avi
.
. http://thetys.dyndns.org/Jupic_2010/jupR_2010-09-29_22-58-16.avi
.

Ce serait mieux si on pouvait récupérer les récents avi de Mars de de Bénédictis au C14.

pour que si un jour tu trouves le temps de nous faire une démo

Bon, je viens de voir que tu avais rajouté des choses plus haut...

Tu nous vois déjà observant en direct sur un écran plutôt qu'à l'oculaire une image grand format, déturbulée, débruitée, amplifiée, diffusable sur le net

Ca viendra bientôt !

[Ce message a été modifié par Nebulium (Édité le 14-02-2012).]

Partager ce message


Lien à poster
Partager sur d’autres sites

Créer un compte ou se connecter pour commenter

Vous devez être membre afin de pouvoir déposer un commentaire

Créer un compte

Créez un compte sur notre communauté. C’est facile !

Créer un nouveau compte

Se connecter

Vous avez déjà un compte ? Connectez-vous ici.

Connectez-vous maintenant