artemis 2 940 Posté(e) 5 janvier Bonjour, Suite au post de @lionello que j'ai trouvé passionnant et grâce à l'exemple qu'il m'a envoyé, j'ai continué mes investigations du côté de CH Cyg, l'avantage de la spectro c'est que l'on peut passer un moment sur l'exploitation et la réductions des données : http://www.astrosurf.com/topic/165923-2spot-les-vents-fragmentés-de-gamma²-velorum-wr-11/ J'ai d'abord repris mon script python pour intégrer dans le graphe la ligne temporelle des résidus autour de la raie Halpha. les résidus sont ceux de la cession en cours et ne sont pas calculés sur plusieurs nuits. J'ai ensuite créé des animations de ces résidus, ici je ne crée que les images individuelles avec python, je passe ensuite sous PSP pour monter le GIF, je ne sais pas si c'est réalisable directement sous python, je n'ai pas encore cherché... Et pour finir, j'ai un peu galéré mais j'ai réussi à sortir le graphe 2D de l'évolution temporelle en fausse couleur, ici j'ajuste le seuil de visualisation selon les conditions de la nuit d'observation. Je suis resté sur la formule de TVS de François Teyssier, celle utilisée par Lionel me donne plus de détection de Flickering avec des pics plus intenses mais je suis étonné par les vitesses radiales qui en ressortent, je cherche actuellement dans la littérature le maximum d'information à ce sujet... Malheureusement CH Cyg a des fluctuations très rapides de l'ordre de la minute et mes données ne me permettent pas d'atteindre cette finesse, je tenterai en 2024 de réduire les temps d'acquisition...au détriment du RSB malheureusement. Premier test : une série de 30 captures de 180s, ma dernière du mois de décembre : le flickering est très léger, mais les données animées ou en fausses couleurs sont plus parlantes on vois bien l'inversion des pics d'intensité sur les ailes bleues et rouges de la raie Halpha : J'ai ensuite repris des cessions un peu plus longues du mois de septembre, d'abord le 9 septembre avec un Flickering bien détecté, essentiellement dû à un sursaut plus marqué lors des 3 dernières captures : Et pour finir avec ce post (j'ai trop de captures en stock ) , ma plus grosse cession du 14/09, ou j'ai passé 110 captures de 150s, le RSB est un peu moins bon et mon script ne valide pas la détection du flickering, mais je suis impressionné par le résultat en fausses couleurs, j'ai l'impression de voir une rotation, avec" enroulement" autour de la rai Halpha ? : Et ce graphe ressemble pour moi à l'ADN de CH Cyg !! Il me reste quelques centaines de spectres 2023, mais je vous fait grâce !!... et je m'y remet dès que possible en 2024... Xavier 4 1 Partager ce message Lien à poster Partager sur d’autres sites
Guillaume BERTRAND 1 649 Posté(e) 5 janvier (modifié) Pfiouuu quel travail Xavier ! C'est énorme As tu montré tes résultats à des pros je pense à J. Merc notamment ? Pour les gif tu peux faire ça comme ça en python : import glob from PIL import Image fp_in = '2023-08-11/_rrlyr*_plot.png' fp_out = "2023-08-11/anim.gif" imgs = (Image.open(f) for f in sorted(glob.glob(fp_in))) img = next(imgs) img.save(fp=fp_out, format='GIF', append_images=imgs, save_all=True, duration=100, loop=0) Modifié 5 janvier par Guillaume BERTRAND 1 Partager ce message Lien à poster Partager sur d’autres sites
artemis 2 940 Posté(e) 5 janvier Merci Guillaume J'envoie mes données brutes à F. Teyssier pour la base ARAS, mais je n'envoie pas les résultats de mes calculs. Merci pour l'élaboration des GIFs, je vais tester et pouvoir supprimer une étape et tout faire avec python, mon script commence à devenir bien copieux à force... Enfin une nuit claire en Touraine, le télescope tourne gentiment... Xavier Partager ce message Lien à poster Partager sur d’autres sites
lionello 309 Posté(e) 6 janvier Les résultats sont très intéressants, essaie de vérifier si la littérature mentionne ces variations. Le changement brutal dans la série du 09/09/2023 est surprenant (vers la 46ème poses). Pour générer les résidus tu utilises la moyenne de la série ou la moyenne de l'ensemble des séries ? Pour le calcul du TVS, j'adopte finalement la méthode de Fullerton ,1996. Il s'agit simplement de tracer un écart type en fonction de la longueur d'onde (la racine du TVS en fait). Partager ce message Lien à poster Partager sur d’autres sites
artemis 2 940 Posté(e) 7 janvier (modifié) Bonjour, Le 05/01/2024 à 17:13, Guillaume BERTRAND a dit : Pour les gif tu peux faire ça comme ça en python : import glob from PIL import Image fp_in = '2023-08-11/_rrlyr*_plot.png' fp_out = "2023-08-11/anim.gif" imgs = (Image.open(f) for f in sorted(glob.glob(fp_in))) img = next(imgs) img.save(fp=fp_out, format='GIF', append_images=imgs, save_all=True, duration=100, loop=0) Guillaume : J'ai testé Glob, mais j'ai un problème de "permission denied", d'après mes recherches c'est un problème récurrent avec Glob, mais je n'ai pas encore approfondi... Le 06/01/2024 à 12:00, lionello a dit : Les résultats sont très intéressants, essaie de vérifier si la littérature mentionne ces variations. Le changement brutal dans la série du 09/09/2023 est surprenant (vers la 46ème poses). Pour générer les résidus tu utilises la moyenne de la série ou la moyenne de l'ensemble des séries ? Merci Lionel , j'ai pas mal fouillé dans la littérature, mais les papiers sur le flickering sont plus orientés photométrie que spectro HR sur Halpha, pour l'instant je n'ai pas encore trouvé d’équivalence... j'ai fini la reprise de la quasi intégralité de mes spectres et je n'ai pas eu de sursaut aussi important que celui du 09/09(?)..., pour les résidus, je les calcule sur la série et non l'intégralité qui dépasse les 400 spectres, par contre j'ai modifié le seuil et figé pour toutes les séries retraitées. Je rajoute ici les visuels en fausses couleur sque je trouve très "parlant". Le pattern spiralé est confirmé sur plusieurs séries, et toutes les variations même moins marquées sont dans la tranche -50 / +50 km/s Du 15/07 au 23/08 cela reste globalement spiralé : Beaucoup plus calme du 3 au 9 septembre matin, le 9 au soir le spiralé semble redémarrer avec le sursaut(?) en fin de cession : Ma plus belle série, du 13/09 matin au 14/09 soir avec 190 spectre au total ! Un début de cession calme avant le retour d'un beau pattern spiralé : Quelques cessions du 19/09 au 08/12, assez calme mais avec de rapides variations le 19/09, avec un aspect "ping-pong" toutes les 300/400s et un retour spiralé le 08/12 mais de faible intensité : Sur 2024 j'essaierai de combiner avec des captures photométriques simultanées... C'est tout pour 2023 Xavier Modifié 7 janvier par artemis 1 Partager ce message Lien à poster Partager sur d’autres sites
lionello 309 Posté(e) 7 janvier Vraiment pas mal ! On voit bien l'intérêt d'un suivi sur toute une nuit, voire sur plusieurs nuits d'affilée. 1 Partager ce message Lien à poster Partager sur d’autres sites
artemis 2 940 Posté(e) 9 janvier Merci Lionel, oui je vais essayer d'allonger les time-séries cette année @Guillaume BERTRAND : je n'ai pas réussi avec Glob pour le GIF, mais j'ai trouvé une autre façon avec imageio dont voici un exemple pour info : import imageio as iio frames = np.stack([iio.imread(f"D://0_CH_Cyg/Sequence_Animation/{i}.png") for i in range(len(Spec1D_List))], axis = 0) iio.mimwrite('D://0_CH_Cyg/Sequence_Animation/gif/TS.gif', frames) Et voila le GIF tout droit sorti de mon script : Xavier Partager ce message Lien à poster Partager sur d’autres sites
