melix

Pas surpris, l'empilement dans JSol'Ex est "primitif" au sens où il ne fait pas de déformation d'image comme AS!3. Par contre une chose intéressante à faire dans le script d'empilement, que je ne faisais pas avant, c'est de faire la correction de l'angle P. Parce que ne serait-ce que sur quelques minutes, l'angle change suffisamment pour que l'empilement en souffre.

La connexion ici fait des siennes... la Dordogne c'est bien sauf pour la 4G Bref, je disais que j'ai essayé de faire qqch de simple mais puissant. J'ai un tuto disponible ici : https://www.youtube.com/watch?v=LytT0EV25SQ Ainsi que la doc : https://melix.github.io/astro4j/latest/fr/jsolex.html#imagemath

Salut à tous, En vacances en Dordogne, j'ai emmené mon Sol'Ex. La semaine passée les nuages étaient présents en matinée ce qui m'empêchait de faire de l'imagerie. Mais ce matin, ô surprise, j'ai eu la chance d'assister en direct à une éruption solaire, quel plaisir ! Il faisait déja très chaud à 9h30, avec une turbu importante. Alors certes les détails ne sont pas au niveau de ce que certains experts postent ici, mais tout de même, c'est fait avec un instrument très modeste ! Voici donc l'animation: Et une version zoomée sur la zone en question: Si vous regardez bien la version complète, vous noterez que certaines images sont tronquées. C'est lié au fait que l'image du disque solaire passe tout juste dans la fente sur mon setup, et qu'il n'était pas parfaitement orienté au nord, ni même réglé à la bonne latitude. Je devais donc faire des corrections manuelles régulièrement, et je me suis parfois raté. Tout le traitement a été fait à l'aide d'un script JSol'Ex dont je vous donne le contenu ici et qui montre qu'on peut faire des choses assez puissantes comme une animation centrée sur l'éruption de façon super simple: [params] # Largeur de correction de bandes band_width=50 # Itérations de corrections de bandes band_pass=20 # Correction de contraste clip=.8 # Dimension d'autocrop (en facteur du rayon solaire) cropFactor=1.2 [tmp] # Crée une image moyenne des décalages de pixels -1 à +1 moy=avg(range(-1;1)) # Applique la correction de bandes fixed=fix_banding(moy;band_width;band_pass) # Applique la correction de contraste et corrige la rotation solaire (angle P) tmp=linear_stretch(clahe(rotate_rad(fixed;angleP);clip)) [outputs] # Génère une image autocroppée frame=autocrop2(sharpen(tmp);cropFactor) # Génère une version colorisée de l'image colorized=colorize(frame;"h-alpha") [[batch]] # Génère une animation monochrome erupt=anim(frame;75) # Génère une animation en couleurs erupt_color=anim(colorized;75) # Génère une animation monochrome centrée sur l'éruption erupt_cropped=anim(crop(frame;130;1130;256;256);75) Par ailleurs, voici une image en H-alpha peu de temps avant l'éruption: Il s'agit d'un stacking de 9 images, réalisé encore une fois entièrement dans JSol'Ex avec ce script: # # Performs stacking of several images processed in batch mode # [params] # banding correction width and iterations bandingWidth=50 bandingIterations=10 # autocrop factor cropFactor=1.2 # contrast adjustment clip=.8 [tmp] corrected = fix_banding(img(0);bandingWidth;bandingIterations) [outputs] contrast_fixed = clahe(corrected;clip) cropped = autocrop2(rotate_rad(contrast_fixed;angleP);cropFactor;32) # This is where we stack images, simply using a median # and assuming all images will have the same output size [[batch]] stacked=sharpen(median(cropped)) protus=linear_stretch(median(disk_fill(cropped;0))) colorized=colorize(stacked;"h-alpha") tech_card=draw_globe(draw_obs_details(draw_solar_params(median(autocrop2(contrast_fixed;cropFactor;32))))) Merci à vous !

Oui j'ai vu effectivement, merci ! Avec plaisir !

Wow, super détails, il marche super bien ce Sol'Ex avec cette lunette ! L'image mix est un peu moche au niveau du limbe par contre (mais c'est ma faute, pas la tienne ;)). Il faudrait que j'améliore l'interpolation au limbe.

Sympa ! Quels paramètres tu rentres dans AS! pour ton stacking ?

Ah oui celle-ci c'est la version brute, juste après correction géométrique. Elle n'applique précisément _pas_ d'ajustement de contraste (CLAHE). Pour les scripts, la doc est ici : https://melix.github.io/astro4j/latest/fr/jsolex.html#imagemath

Bonjour à tous, Enfin ! Après des jours, que dis-je des semaines de disette, le Sol'Ex est enfin ressorti. Je suis en vacances en Bretagne, dans un terrain sans électricité (mis à part des panneaux solaires pour recharger des batteries) et donc la manipulation n'est pas des plus aisée Néanmoins, il semble que je sois plus doué en développement logiciel qu'en observations. C'est aussi pour ça que je participerai au stage AIP fin septembre ! Ceci dit, j'ai eu pas mal de difficultés pour cette série. En particulier, une mise en station de jour un peu approximative, m'obligeant entre chaque vidéo à "recentrer", mais surtout j'ai un mini-PC qui tient à peine la cadence et a beaucoup de mal à soutenir le flux USB. Ca donne des vidéos tronquées qui finissent par vous donner des images bien moches comme celle-ci : Le pb c'est que j'ai plus de vidéos comme celle-là que de vidéos non tronquées ! Mais bon, au final, j'ai réussi à faire quelques images. Voici le soleil du jour en H-alpha. Tous les traitements ont été faits avec mon logiciel maison, JSol'Ex. Une animation du continuum (-10 à +10 pixels) : Enfin, j'ai voulu faire ma premìere observation en raie Hélium. J'ai donc utilisé la technique décrite dans la documentation qui permet de n'avoir à prendre qu'une seule vidéo. L'image n'est pas extraordinaire, mais c'est ma première à cette longueur d'onde. Encore une fois, traitement entièrement fait dans JSol'Ex : Pour réussir à sortir quelque chose, j'ai du réduire le gain à 0 et diminuer le temps d'exposition à 5 ms. Loin d'être idéal, mais sans ça je n'obtenais rien de satisfaisant. Sur ce, j'arrive bientôt à court de batterie, il va donc falloir rebrancher les chargeurs solaires !

Oui c'est assez fréquent (et énervant) avec ce genre de logiciels. N'hésite pas à leur rapporter. Mais tous les sources sont dispos, donc facile de vérifier s'il y a un malware... Je vais ajouter des paramétrages pour les valeurs par défaut de CLAHE, mais en attendant (ça peut prendre un peu de temps), tu peux tester différents paramétrages assez facilement. A la fin du traitement, par exemple, tu peux mettre un script en bas à droite. Dans ce script, mets: [outputs] clahe_0_5=clahe(img(0);0.5) clahe_0_8=clahe(img(0);0.8) clahe_1_0=clahe(img(0);1.0) Ca te permettra de comparer la sortie avec différents niveaux de clipping (ici je compare avec 0.5, 0.8 (valeur par défaut) et 1.0. Tu peux même aller plus loin en changeant la taille des tuiles utilisées pour calculer: [outputs] # Variation taille tuile clahe_0_8_64_64_t=clahe(img(0);64;64;.8) clahe_0_8_128_64_t=clahe(img(0);128;64;.8) clahe_0_8_256_64_t=clahe(img(0);256;64;.8) # Variation nombre de paniers clahe_0_8_64_64_p=clahe(img(0);64;64;.8) clahe_0_8_128_64_p=clahe(img(0);64;128;.8) clahe_0_8_256_64_p=clahe(img(0);64;256;.8) A toi de jouer !

Oui, ceci étant j'enregistre sur un Samsung T7 externe qui supporte largement la charge. Le souci vient vraiment de l'unité centrale, c'est un mini PC entrée de gamme, avec un Celeron 3000. L'avantage c'est qu'il ne consomme presque rien (10W), ce qui est plutôt chouette pour la batterie, mais il n'aime pas trop les montées en charge... Ah, mais c'est normal : tous les algos sont conçus pour du monochrome. Cette case, ça permet de dire "même si ma vidéo dit qu'elle est en couleur (e.g Bayer RGGB), n'essaie pas de debayeriser parce qu'elle est en fait en mono". C'est typiquement ce qui arrive lorsque tu fais du binning 2x2 avec une caméra couleur. L'image devient monochrome et peut être traitée telle quelle. Donc si tu as une "vraie" vidéo couleur, elle va être interprétée comme une mono, ce qui peut donner des résultats rigolos.

Je suis déjà en bin 2x2 (mono bin, c'est une caméra couleur)

En y réfléchissant, ça ne devrait rien changer. L'image Doppler combine 2 images à pixels -3 et +3 en les associant aux canaux couleur. Que je fasse la rotation +3 et -3 indépendamment puis combine, ou que je fasse la combinaison _puis_ la rotation, le résultat sera le même (et mes dessins semblent confirmer). J'aurais donc tendance à dire que si tu te retrouves avec le bleu en bas, c'est qu'il ne fallait justement pas faire une rotation (cas du scan en AD). Ou alors j'ai loupé qqch ?

En l'occurrence je ne parle pas de la chromosphère, mais des arches vus "de face" Perso je préfère le fond blanc.

Oui, c'est ce qui est fait pour le fiche technique. Je vais regarder ça.

Il n'y a pas à dire, les images inversées ça aide vraiment à voir le relief !

Ca y est, la 1.6.2 est disponible, suite à vos retours : - correction du bug de l'image Doppler manquante - correction de l'image noire sur le script de colorisation de la fiche technique - ajout d'une option de rotation à gauche ou à droite (pour les scans en DEC) - correction des images tronquées lors des rotations d'angle P sur des disques partiels - ajout d'une option d'autocrop dans les paramètres (avant l'autocrop n'était disponible que dans les scripts) - légères modifications des paramètres CLAHE par défaut - encore des améliorations dans l'algorithme de détection de contours Merci à tou ·te ·s !

Ok pas de soucis je regarde demain.

Effectivement, j'ai cassé cette fonction en 1.6.0 (en tt cas dans l'interface principale), merci de l'avoir fait remarquer ! En attendant que je corrige, tu peux en produire une via ImageMath (il y a un script d'exemple).

Ok hier j'ai fais une fausse manip lors de la sortie de la 1.6.1, pendant qqs minutes c''était pas la bonne version de publiée. Est-ce que tu pourrais essayer de désinstaller ta version et re-télécharger la 1.6.1 ? Pour info je suis en train d'ajouter la possibilité de faire une rotation à gauche/droite directement dans les paramètres. Ca devrait être dispo dans la 1.6.2.

Tu est bien sur la 1.6.1 ?

Magnifique,ça me fait presque regretter de n'avoir acheté le mien qu'il y a 4 mois !

Je suppose que tu veux dire que tu fais les acquisitions en DEC, non ? Parce qu'en AD c'est la direction "usuelle". Je planifie d'ajouter une option "Direction du balayage" dans une prochaine version. Ceci dit, c'est intéressant, parce que même si tu n'as pas ces paramètres, tu peux tout à fait générer ceci toi-même avec les scripts. C'est précisément conçu pour ça : en gros le programme fait des choix "simples" pour la majorité, mais si ça ne colle pas, tu peux passer en mode "expert" et générer ce qu'il te manque. Avant tout, passe à la version 1.6.1, parce qu'en écrivant ce message j'ai trouvé un bug que j'ai corrigé. Prenons cet exemple. Tu me dis "il faut tourner de 90°", et "je voudrais une fiche technique colorisée". Pas de problème, allons-y ! Commençons par ouvrir la fenêtre de paramètres : De là, sélectionne "ImageMath" et décoche toutes les autres images: Maintenant clique sur "Ouvrir ImageMath" et sélectionne le script d'exemple "Fiche technique" : Le script fait ceci: [params] contrast=1.0 cropFactor=1.2 [tmp] _clahe=sharpen(clahe(img(0);contrast)) cropped=autocrop2(_clahe;cropFactor) globe=draw_globe(cropped) [outputs] techcard=draw_solar_params(draw_obs_details(globe)) Si on part de la fin, il dit "crée une sortie nommée techcard` qui applique la fonction draw_solar_params sur l'image draw_obs_details qui elle-même est appliquée sur la variable "globe". On a donc "composé" une image à partir du globe, qui ajoute les détails d'observation et les paramètres solaires. Très bien, mais d'où vient globe ? Dans la section `tmp`, on trouve `globe=draw_globe(cropped`). C'est donc une image qui dessine le globe sur l'image cropped. Et on trouve cropped juste au-dessus : il s'agit d'un autocrop de l'image CLAHE, elle même définie au dessus comme l'augmentation des détails (sharpen) de CLAHE de l'image au pixel shift 0. Ah ah, on tient donc notre coupable, cette image 0 devrait être tournée de 90 degrés. On peut donc modifier le script en écrivant: [params] contrast=1.0 cropFactor=1.2 [tmp] rotated=rotate_left(img(0)) _clahe=sharpen(clahe(rotated;contrast)) cropped=autocrop2(_clahe;cropFactor) globe=draw_globe(cropped) [outputs] techcard=draw_solar_params(draw_obs_details(globe)) J'ai donc ajouté un appel à la fonction "rotate_left" qui fait ce qu'elle dit : une rotation à gauche, de 90 degrés. Reste à sauvegarder le script : On clique sur "Ok" puis on lance le traitement. Tu verras ta fiche avec le disque orienté correctement. Dernière étape, avoir la version couleur. Pour celà, on modifie une fois de plus le script : [params] contrast=1.0 cropFactor=1.2 [tmp] rotated=rotate_left(img(0)) _clahe=sharpen(clahe(rotated;contrast)) cropped=autocrop2(_clahe;cropFactor) globe=draw_globe(cropped) color_globe=draw_globe(colorize(cropped;"h-alpha")) [outputs] techcard=draw_solar_params(draw_obs_details(globe)) techcard_color=draw_solar_params(draw_obs_details(color_globe)) J'ai ajouté `color_globe`, qui utilise la fonction "colorize" pour faire la colorisation du disque en utilisant le profil h-alpha. Puis, on calcule dans [outputs] une 2ème image qui combine une nouvelle fois les informations techniques avec le globe colorisé. Et voici le résultat : Au contraire, merci

A qui le dis-tu ! Je code plein de trucs mais j'ai même pas loisir à l'essayer sur mes propres captures. J'ai fais ce post notamment avec une "vieille" vidéo de Juin Bientôt les vacances pour moi, j'espère avoir plus l'occasion de faire quelques captures, c'est difficile quand tu bosses de sortir pour qqs minutes quand il y a une éclaircie !

Bonjour à tous, Faute de pouvoir imager moi même, il pleut en Bretagne, je continue mes développements Lorsqu'on traite une vidéo Sol'Ex, on utilise généralement le "milieu" de la raie détectée (typiquement le milieu de la raie h-alpha), mails il est aussi possible de reconstruire une image en faisant un décalage de pixels par rapport à ce milieu. Voici un exemple de détection de raie sur la vidéo de test de traitement disponible sur le site de Christian : Ainsi par exemple en décalant de 10 ou 15 pixels, on se retrouve dans le "continuum". J'avais donc ajouté à JSol'Ex la possibilité de calculer ces images décalées, mais j'avais aussi, arbitrairement, limité cette possibilité à des nombres entiers (1 pixel, -2 pixels, etc...). Or, de part la précision que l'on peut obtenir lors de la reconstruction, inférieure au pixel, et l'interpolation lors de la reconstruction, il n'y avait pas de raison de limiter. Voici par exemple la différence entre l'image centrée (décalage = 0) et celle décalée de 0.2 pixels : Ceci est particulièrement intéressant si on souhaite faire une animation, qui part par exemple d'un décalage de -12 pixels à +12 pixels, pour aller "d'un continuum à un autre". En effet, on peut dès lors obtenir une animation très fluide. Voici un exemple: Il suffit de 45 secondes sur ma machine pour générer cette animation, en tout automatique ! Voici le script que j'ai utilisé : [params] shift=12 step=0.2 [tmp] images=range(-shift;shift;step) fixed=remove_bg(fix_banding(images;25;3)) cropped=autocrop2(fixed;1.1;32) adjusted=clahe(cropped;1) [outputs] video=anim(adjusted;40) On peut même tenter une version colorisée, même si je suis moins convaincu du résultat: avec ce script: [params] shift=12 step=0.2 [tmp] images=range(-shift;shift;step) fixed=remove_bg(fix_banding(images;25;3)) cropped=autocrop2(fixed;1.1;32) adjusted=clahe(cropped;1) colorized=colorize(adjusted;"h-alpha") [outputs] video=anim(colorized;40) Tout ceci est disponible dans JSol'Ex 1.5.2. A noter que vous devez être raisonnable sur l'incrément (0.2) et la taille de l'intervalle (-10;10) parce qu'il va y avoir un grand nombre d'images à générer et stocker en mémoire. J'ai optimisé l'utilisation mémoire dans cette version, mais la reconstitution des images à différents décalages de pixels demande de pouvoir stocker chaque image en mémoire. J'essaierai de faire d'autres optimisations dans une prochaine version, mais utiliser moins de RAM veut aussi dire utiliser plus de disque. SSDs indispensables !