Fluctuations de CH Cyg : Septembre...

[Vincent_L 546]

bravo pour cette manip !

en fait tu calcules un coefficient de variation ("fractionnal variance)"  pour tes séries temporelles sur chaque lambda. Comment s'assurer que les phénomènes observés sont propres à la physique de CH Cyg ? Est ce qu'on observerait la même chose, au même moment, sur 10 Lac par exemple  ? Est ce que ces variations sont périodiques ? Ca nécessiterait effectivement de multiplier les observations

 

 

 

[lionello 309]

il y a 22 minutes, Vincent_L a dit :

Comment s'assurer que les phénomènes observés sont propres à la physique de CH Cyg ?

 

Salut Vincent,

 

D'après ce que je comprends c'est parce que les fluctuations de la raie sont comparées aux fluctuations du continuum.

Si les fluctuations de la raie sont supérieures à 3 sigma du bruit de fond alors on peut attribuer ces variations à autre chose que du bruit !

 

Le flux des poses unitaires peut varier plus ou moins fortement à cause du seeing, des défauts de suivi. Et ces phénomènes sont accentués par une fente étroite ou l'angle parallactique (http://www.astrosurf.com/buil/dispersion/atmo.htm). Les fichiers @pro générés par ISIS sont très instructifs.

 

Donc il faut faire ce type de manip dans de bonnes conditions pour avoir le moins de fluctuations parasites possibles !

 

Modifié 2 septembre 2023 par lionello

[artemis 2 940]

Il y a 5 heures, lionello a dit :

Pour être sûr de bien comprendre, tu détermines le continuum de chaque pose (fit_generic_continuum() ) , tu en déduis sa variance sur l'ensemble des poses autour de la raie Halpha. Tu calcules le sigma(lambda) de la raie, et les variations supérieures à 3 sigma du continuum sont interprétées comme du flickering ?

Bonjour Lionel !

Oui tu as tout compris, c'est exactement l'esprit de cette manip

 

Il y a 5 heures, Vincent_L a dit :

bravo pour cette manip !

en fait tu calcules un coefficient de variation ("fractionnal variance)"  pour tes séries temporelles sur chaque lambda. Comment s'assurer que les phénomènes observés sont propres à la physique de CH Cyg ? Est ce qu'on observerait la même chose, au même moment, sur 10 Lac par exemple  ? Est ce que ces variations sont périodiques ? Ca nécessiterait effectivement de multiplier les observations

 

Merci Vincent, effectivement cette manip mériterait un suivi plus serré, je ne sais pas si CH_Cyg va faire l'objet d'un suivi STAROS mais l'idée avait été émise par Christian il y a quelques temps... l'amplitude ou une périodicité serait plus facile a détecté si elle existe...Je ne sais pas si d'autres symbiotiques présentent les même variations... encore un domaine à explorer...

 

Pour la comparaison aux 3 sigma sur le continuum, c'est une valeur arbitraire adoptée par François Teyssier comme base de détection du flickering. Cette observation nécessite un bon RSB, François préconise un RSB de 100 autour de Halpha couplé a des poses unitaires inférieures à 5mn. Il y a beaucoup de flux sur cette étoile mais je ne pense pas pouvoir descendre sous les 180s de pose unitaire pour garder un continuum correct...

 

Il y a 5 heures, lionello a dit :

Le flux des poses unitaires peut varier plus ou moins fortement à cause du seeing, des défauts de suivi. Et ces phénomènes sont accentués par une fente étroite ou l'angle parallactique (http://www.astrosurf.com/buil/dispersion/atmo.htm). Les fichiers @pro générés par ISIS sont très instructifs.

Merci pour le lien Lionel, encore une page de référence de Christian !! Ici en HR la plage spectrale est très étroite et je travaille avec un instruments à miroirs, mais je vais approfondir ce sujet   Actuellement CH Cyg est bien placée et haute dans le ciel, un bon point déjà

 

Xavier

[lionello 309]

Cool ! 

Peux-tu me transmettre ta série d'observation du 16 juillet ?

Je vais tenter de bricoler un truc.

Pourquoi as-tu tout traiter tes spectres unitairement, alors que les fichiers intermédiaires @pro auraient probablement suffi ?

 

[Guillaume BERTRAND 1 649]

Absolument passionnant. Bravo pour toutes ces données CH Cyg Xavier ! Et la mise en forme est très sympa. Ça progresse en python à vitesse grand V ! 

[artemis 2 940]

Il y a 2 heures, lionello a dit :

Peux-tu me transmettre ta série d'observation du 16 juillet ?

Lionel,

Je viens de t'envoyer un ZIP de mes captures de juillet

En fait j'ai traité jusqu'au bout par habitude, je ne savais pas au début que l'on pouvait exploiter les @pro, et j'ai trouvé une technique pas trop contraignante en attendant que Christian nous fasse une petite option , et aussi car je suis en flat "LED" et même en HR j'ai une forte pente sur le continuum si je ne traite pas jusqu'au bout...

 

Il y a 2 heures, Guillaume BERTRAND a dit :

Absolument passionnant. Bravo pour toutes ces données CH Cyg Xavier ! Et la mise en forme est très sympa. Ça progresse en python à vitesse grand V ! 

Merci Guillaume , oui ce projet est un bon laboratoire pour mon apprentissage de Python, mais que j'ai galéré avec les unités Jy, Jy2 etc... !!  sur lesquelles je ne pouvais pas faire le calcul de la racine par exemple... beaucoup de bidouilles pour arriver au bout !!

Je vais essayer de me faire une interface  Windows sous Lazarus pour lancer le script python et paramétrer avec un fichier YAML...

 

J'ai fait l'impasse CH Cyg cette nuit car j'étais en config planétaire sur le DK... j'aurai pas dû , une turbu infâme ce matin... une misérable CH4 pour l'Alpo et puis c'est tout..., cet aprem je passe en 35µ et je crois que je vais rester un moment en spectro sur le DK !!

 

Xavier

[Vincent_L 546]

Il y a 8 heures, artemis a dit :

Pour la comparaison aux 3 sigma sur le continuum, c'est une valeur arbitraire adoptée par François Teyssier comme base de détection du flickering

pas tout à fait arbitraire, c'est plutôt une bonne idée d'aller chercher en dehors des 3 sigmas pour détecter une évènement: Règle 68-95-99,7 — Wikipédia (wikipedia.org)

[artemis 2 940]

il y a 9 minutes, Vincent_L a dit :

pas tout à fait arbitraire, c'est plutôt une bonne idée d'aller chercher en dehors des 3 sigmas pour détecter une évènement: Règle 68-95-99,7 — Wikipédia (wikipedia.org)

Merci Vincent !! je viens encore d'apprendre quelque chose aujourd'hui... mes cours de stats sont si loin !!

 

Xavier

[lionello 309]

Il y a 3 heures, artemis a dit :

Lionel,

Je viens de t'envoyer un ZIP de mes captures de juillet

 

Je suis loin d'être un pro en Python, j'ai attaqué la semaine du 15 août 

 

Voilà le résultat :

 

 

On retrouve la même forme de flickering.

 

Mais je galère pour le sigma. L'ajustement du continuum est très sensible.

J'utilise la commande fit_generic_continuum(spectrum, exclude_regions=[....])

C'est un peu du bricolage je trouve.

 

Comment as-tu fait ?

 

Autre point, il me semble que la fractional variance, c'est divisé par N (c'est ce que j'ai tracé) et non pas par N-1.

A voir...

 

 

 

 

 

 

[artemis 2 940]

il y a 12 minutes, lionello a dit :

Je suis loin d'être un pro en Python, j'ai attaqué la semaine du 15 août 

Cool Lionel tu as été rapide !! je n'ai guère plus d'expérience que toi, mais la tu m’ébahis

 

Le profil est similaire c'est un bon  point !!

Pour le continumm j'ai également utilisé Fit_generic_continuum, mais après différents tests, j'ai supprimé l'exclusion autour de Halpha et le niveau à 3 sigma était plus cohérent avec le graphe de François, le RSB étant  plus raisonnable également...

 

Pour la variance j'étais effectivement parti sur N mais j'ai bien trouvé la référence pour la variance corrigée avec N-1, c'est bien celle utilisée par François...

 

Pour les autres calcul j'ai tout décomposé, je n'ai pas réussi à utiliser .std, tu as fait comment ?

 

Si cela t’intéresse je peux t'envoyer les 16 captures du 25/08 avec les 2 pics  pour comparer ?

 

Xavier

[lionello 309]

il y a 50 minutes, artemis a dit :

Cool Lionel tu as été rapide !! je n'ai guère plus d'expérience que toi, mais la tu m’ébahis

 

Je ne t'ai pas tout dit : j'ai fini ça hier à 2h de mat', mon programme n'attendait plus que tes données 

 

il y a 50 minutes, artemis a dit :

Pour la variance j'étais effectivement parti sur N mais j'ai bien trouvé la référence pour la variance corrigée avec N-1, c'est bien celle utilisée par François...

 

Peux-tu me donner la référence en question ?

 

il y a 50 minutes, artemis a dit :

Pour le continumm j'ai également utilisé Fit_generic_continuum, mais après différents tests, j'ai supprimé l'exclusion autour de Halpha et le niveau à 3 sigma était plus cohérent avec le graphe de François, le RSB étant  plus raisonnable également...

 

C'est quoi du coup tes zones d'exclusion ? Moi j'ai fait truc un peu foireux comme ça :

g1_fit = fit_generic_continuum(spectrum, exclude_regions=
                [
                 SpectralRegion(0 * u.AA, 6550 * u.AA),          
                 SpectralRegion(6557 * u.AA, 6564 * u.AA),
                 SpectralRegion(6577 * u.AA, 6592 * u.AA),
                 SpectralRegion(6690 * u.AA, 8000 * u.AA)
                 ])

 

il y a 50 minutes, artemis a dit :

Pour les autres calcul j'ai tout décomposé, je n'ai pas réussi à utiliser .std, tu as fait comment ?

 

J'ai fait une boucle qui calcule la moyenne des continnuum (mean_continuum_flux) sur une plage donnée entre 6550 et 6570 A (Ca ne sert à rien fait car les continuum sont parallèles si on s'y prend bien, donc j'aurais pu prendre la valeur du continuum à 6563 par exemple).

Je stocke chaque valeur dans un tableau numpy (mean_continuum [ ]): 

 mean_continuum.append([mean_continuum_flux])

 

Et pour avoir l'écart type des valeurs tu fais juste :

import numpy as np

std_continuum = np.std(mean_continuum_values)

 

EDIT : si tu voulais parler de la fractional variance, j'ai été obligé de décomposer aussi.

 

 

Citation

Si cela t’intéresse je peux t'envoyer les 16 captures du 25/08 avec les 2 pics  pour comparer ?

 

Oui ok !

 

 

Modifié 3 septembre 2023 par lionello

[Vincent_L 546]

il y a 50 minutes, lionello a dit :

Autre point, il me semble que la fractional variance, c'est divisé par N (c'est ce que j'ai tracé) et non pas par N-1.

on divise par n si on a accès aux mesures de toute la population. En général, c'est pas le cas (ici notamment), car on a accès à un échantillon (n) de la population (N). Pour corriger ce biais, on divise plutôt par n-1. Si n est très grand, alors ça n'a plus vraiment d'importance ... . Si n est petit, il faut faire la correction ...

Modifié 3 septembre 2023 par Vincent_L
--

[lionello 309]

On note que Vincent c'est un pro en stat.

C'est bon à savoir !  

[artemis 2 940]

il y a 33 minutes, lionello a dit :

 

Peux-tu me donner la référence en question ?

Vincent à répondu ci-dessus, sinon en cherchant sur Variance N-1 j'avais trouvé pas mal de pages d'explications, comme celle-ci par exemple :

https://fr.khanacademy.org/math/be-4eme-secondaire2/x213a6fc6f6c9e122:pour-aller-plus-loin/x213a6fc6f6c9e122:variance-d-un-echantillon-vs-population/v/review-and-intuition-why-we-divide-by-n-1-for-the-unbiased-sample-variance

 

il y a 36 minutes, lionello a dit :

C'est quoi du coup tes zones d'exclusion ?

Et bien je n'ai aucune zone d'exclusion finalement, et le 3sigma est tombé plus prêt du graphe de François !

 

Ok pour np.std, j'essaierai dans une prochaine version de mon script ..

 

 

Merci, Je t'envoie les fichiers en MP

Xavier

[lionello 309]

Je croyais que tu parlais d'une publi sur le flickering de Ch Cyg.

 

il y a une heure, artemis a dit :

Et bien je n'ai aucune zone d'exclusion finalement, et le 3sigma est tombé plus prêt du graphe de François !

 

Ah ben niveau dispersion du continuum j'ai plutôt ça 

 

 

Mais à vue d'œil, correction en N-1 faite, nos résultats sont très similaires :

 

 

Modifié 3 septembre 2023 par lionello

[artemis 2 940]

Merci Lionel, les résultats sont cohérents c'est plutôt rassurant !

 

En fait le continuum je le calcule sur l'ensemble du spectre, de 6500 à 6700 angström, cela induit peut-être quelques variations, mais la variance doit rester approchante (?)...

 

Xavier

[artemis 2 940]

Bonjour,

 

Voici quelques séries sur CH_Cyg de ce début septembre, un ciel laiteux et chargé de sable du Sahara a compromis un peu la résolution mais le beau temps revient !

J'ai encore amélioré mon script pour le tracé des spectres et j'ai détecté et corrigé un petit défaut de dispersion entre 6600 et 6700 angström , je reprendrai l’intégralité de mes données si Christian nous fait une petite option pour les traitements individuels, j'ai plus de 230 spectres à reprendre...

 

Pour tous ces spectres je suis passé sur une fente de 35µ, je vois bien le gain sur quelques Be faibles faites en parallèle, au dernier réglage je stabilise autour de R 16000 que ce soit sur le néon ou lors du test sur les telluriques de SpecInti, cela me convient fort bien !

 

Le 3 septembre, calme plat pour le flickering ( le flux max depuis le début du mois reste en dessous de 60), j'ai ajouté un graphe pour suivre le flux selon les poses, et un graphe de contrôle visuel de la calibration sur les 3 raies H2O les plus intenses de la zone car le contrôle des raies telluriques de SpecInti ne fonctionne pas sur CH Cyg (je vérifie également avec VisualSpec), sans doute a cause de l'intensité de la raie H alpha :

 

 

 

Le 5 septembre, très léger Flickering sur l'aile rouge de la raie Halpha, j'attends de collationner plus de données, mais j'ai déjà constaté sur plusieurs séries une fluctuation rapide, en dent de scie comme ici entre les poses 6 et 12 par exemple, suivi soit par une baisse, soit par une remontée du flux... je ne sais pas si c'est une coïncidence, mais je vais continuer à suivre... le seeing était mauvais le RSB s'en ressent

 

 

Le 6 septembre, petite série et très léger mouvement sur les ailes Halpha, légère phase descendante  du flux :

 

 

Le  9 septembre j'ai repris le réglage de la caméra science, il n'était pas optimisé suite au changement de fente, le résultat est nettement meilleur !!. je détecte un flickering plus intense sur chaque aile de la raie Halpha, on est  encore sur une phase descendante du flux :

 

 

J'espère repasser une série ce soir, avant ou après quelques Be... suite au prochain post...

 

EDIT :

Après avoir "cramé" 30 minutes de pose de 240s, j'ai donc relancé hier soir une série de 50 poses de 180s, belle remontée du flux pour un Flickering marqué sur l'aile rouge de la raie Halpha, j'ai par contre une forte dispersion sur le continuum (j'ai indexé les barres d'incertitudes sur cette fluctuation), je pense que j'ai focalisé trop fin sur le spectre... je suis descendu à 20 en BIN_SIZE, et que le ciel devait être un peu changeant (?) :

 

Il va falloir que je reprenne étalonnage des ces 50 captures, la dispersion est un peu trop forte sur H2O. j'attends la prochaine version de SpecInti pour re-mouliner l'ensemble...

Merci Vincent pour le "like"

Xavier

 

Modifié 10 septembre 2023 par artemis

[artemis 2 940]

Bonjour,

 

Une séance improvisée cette nuit, calme plat sur H alpha :

 

Par contre je m'aperçois que l'exploitation reste compliquée, ci-dessous la récap de mes observations..., il faudrait beaucoup d'observations pour "voir" ce qui se passe sur cette étoile...

 

Xavier

[Guillaume BERTRAND 1 649]

Excellent ce suivi Xavier

[artemis 2 940]

Le 13/09/2023 à 14:43, Guillaume BERTRAND a dit :

Excellent ce suivi Xavier

Merci Guillaume !

 

 

J'ai passé 2 séries hier 14/09, d'abord le matin au réveil vers 2h00, peu d'amplitude sur le flux max des poses  mais le Flickering est bien détecté :

 

 

Le soir je me suis donné un créneau de 7h00 maxi pour essayer une séance marathon, j'ai donc pointé de bonne heure , contrepoids en haut pour ne pas avoir de retournement à faire au méridien... malheureusement au bout de 4h30 quelques nuages on fait décroché l'autoguidage... Et je dois avouer que j'ai manqué de vigilance... J'avais légèrement piqué du nez !!

Belle évolution du flux mais... Calme plat autour de H alpha !! :

 

 

Bon cette nuit je crois que je vais repasser sur quelques Be pour me détendre , le traitement des longues séries est encore fastidieux...

 

Par contre François m'a confirmé hier que les pros étaient intéressés  par nos "Times Series",  alors... à vos spectros !!

 

Xavier

 

[maxip 339]

Bonjour,
je suis passé dans le coin du Cygne le 14 septembre avec un starex BR (300) histoire d'avoir une image plus globale des phénomènes de fluctuations.
C'est pas vraiment la bonne résolution pour un suivi détaillé de H apha mais voici mon observation à partir de 37 spectres de 60 s.
Je ferai un autre passage avec un starex 2400 sur Halpha... des que j'aurais la bonne fente.

Maxime

 

Modifié 17 septembre 2023 par maxip

[artemis 2 940]

Bonjour Maxime,

 

Super ce suivi de CH Cyg en BR, tu envoies tes spectre à la base ARAS ?

Je ne suis pas équipé pour de la BR, mais je trouve cette étoile très intéressante à cette résolution. Il y a une belle variation sur Halpha,  as tu essayé de calculer et vérifier le Flickering ? même sur 37mn c'est déctectable normalement...

 

Bonne continuation

Xavier

[Guillaume BERTRAND 1 649]

Salut Maxime, 

Très intéressant ton travail en BR ! 

Je suis juste étonné de l'amplitude des variations... Sont elles réelles ? Je me pose la question car sur les spectres HR de Xavier les amplitudes des variations sur la raie Halpha sont tout de même assez faible.... Qu'en pensez vous ? Xavier as tu déjà vu des travaux similaires en BR sur cette étoile ? 

 

Bonne soirée 

Guillaume 

[maxip 339]

il y a 51 minutes, artemis a dit :

as tu essayé de calculer et vérifier le Flickering

Non pas encore, c'est prévu ! Il me faut juste le temps de digérer ... le Python !

En fait j'ai pas mal de bricoles entre autres changer l'objectif camera 80mm pour passer sur un optimisé que j'ai reçu et que je viens de monter. J'ai également prévu de revenir avec mon Starex HR pour cibler la raie Halpha, j'attend la fente (19microns) que j'ai commandé cette semaine, ca permettra de comparer avec vos acquisitions.

 

il y a 45 minutes, Guillaume BERTRAND a dit :

Je suis juste étonné de l'amplitude des variations

Oui moi aussi en 30 mn, ça bouge beaucoup et pas que sur Halpha, regarde entre 4000 et 4400. 

 

je trouve cette étoile très intéressante,   je vais rester un peu dessus et compléter avec des acquisition en HR.

 

Cette étoile fait l'objet d'un suivi et de travaux sur ARAS :

https://aras-database.github.io/database/chcyg.html

http://www.astrosurf.com/aras/surveys/chcyg/index.html

 

[artemis 2 940]

Il y a 6 heures, Guillaume BERTRAND a dit :

Je suis juste étonné de l'amplitude des variations... Sont elles réelles ? Je me pose la question car sur les spectres HR de Xavier les amplitudes des variations sur la raie Halpha sont tout de même assez faible.... Qu'en pensez vous ? Xavier as tu déjà vu des travaux similaires en BR sur cette étoile ? 

Bonjour Guillaume, loin d'être spécialiste, je pense qu'il faudrait poser la question à François Teyssier. Je suis également étonné de l'amplitude, mais l'unité du graphe sur Y me trompe sans doute. il y a cependant 6 ou 7 spectres que je trouve bien plat sur Halpha.

Ce que dit la littérature, sur Aras par exemple, est que CH_Cyg présente un Flickering à l'échelle de la minute, en HR je suis actuellement à 2mn30 par pose unitaire, cela lisse peut-être mais sur 2 à 3h je devrais tomber sur un pic plus souvent s'il y en avait de très fort (?)... par contre il m'est arrivé de trouver un flux relatif vers 40 le matin et a plus de 60 le soir du même jour...

En BR la variation sur tout le spectre est impressionnante en tout cas...

 

Il y a 5 heures, maxip a dit :

Non pas encore, c'est prévu ! Il me faut juste le temps de digérer ... le Python !

Oui Maxime, je me suis mis au python récemment, j'ai mis pas mal de temps pour mettre en place un petit script d'analyse !

Tu as de beaux projets en perspective, c'est cool... Vivement tes premières données HR en complément

Il y a 5 heures, maxip a dit :

Cette étoile fait l'objet d'un suivi et de travaux sur ARAS :

https://aras-database.github.io/database/chcyg.html

Je pense que tu as envoyé tes données à la base ARAS, j'y envoie les spectres HR et  également un ZIP de la totalité des spectres unitaires (François m'a confirmé l'utilité des Times Series).

 

Ce matin je fais une série courte sur CH_Cyg en ce moment ( d'une quinzaine de poses ) et je file sur M77 pour ma première Galaxie en HR... si tout se passe bien !!

 

Xavier