Capella13

En inversant l'ordre du traitement. Etirement d'abord, débruitage et colorisation ensuite. Ca règle le problème de la texture granuleuse. Un crop pour nettoyer le bas de l'image et un rehaussement de l'intensité du fond de ciel comme suggéré. Ce que je prenais pour un artefact bleu à gauche semble bien réel (présent sur d'autres photos de la nébuleuse)

Merci pour le commentaire. La jonction entre les tuiles est loin d'être parfaite. Probablement à cause des temps de poses différents d'une tuile à l'autre. D'autre part mais c'est lié au problème précédent, le retrait de gradient a été un peu compliqué aussi. Un fond de ciel plus lumineux fait ressortir ces défauts ou apparaître des nébulosités fantômes. La colorisation Halpha=... et O3=... la déconvolution et le débruitage sont faits avant étirement. Le rehaussement de contraste (CLAHE) après. Le changement de couleur du O3 (cyan-> bleu azur) c'est indifférent. Le gain de la déconvolution en terme de résolution ne vaut peut être pas ici son coût; l'aspect grumeleux visible à fort grossissement.

Voici une mosaique 2x2 faite sur 3 nuits avec mon Newton Skywatcher T200/1000 matériel : correcteur de coma, caméra ASI 2600 MC pro, filtre L extreme, autoguidage avec ASI 120 mini et lunette 60/240. Logiciel et interface matériel : ASCOM, phd2 guiding, cartes du ciel et APT. poses cumulées 8h44 par des nuits un peu humides sans lune. soit dans le détail : 84x2' sur les 2 tuiles de gauche, 59x2' sur celle en haut à droite et 35x2' sur le dernier quadrant. erreur de ma part ces différences responsable d'un rapport S/B non homogène ont bien compliqué le traitement. résolution moyenne FWHM = 3.0" selon Siril avant traitement et étirement. traitement DOF, puis SIRIL retrait de gradient polynôme ordre 1 avant assemblage des tuiles (recouvrement 10 à 15%). Assemblage de la mosaïque par astrométrie (plus galère que d'habitude une des tuiles ne voulait pas se placer correctement. Il a fallu la placer en manuel). Crop, en gardant une petite bande noire (partiellement invisibilisée) en bas a droite pour avoir toute la nébuleuse à gauche. Puis retrait de gradient RBF Siril en choisissant soigneusement les points (mieux que l'IA de GraXpert cette fois ci.) Enfin débruitage IA de graXpert. A ce stade ça donne ça : La couleur orangée ou jaune témoigne de la présence de 0III en plus du Halpha Traitement HOO : H=R, O=2.5*(V+B)-.075*R. puis divers traitements pour réduire le bruit du signal OIII à cause de temps de pose trop courts sur les tuiles. Un peu de déconvolution sur le canal rouge, augmentation du contraste par CLAHE. Légère modification des canaux V =0.88 O3 et B=1.12 O3 pour avoir un bleu qui tire un peu vers le bleu azur plutôt que le cyan. puis étirement hyperbolique avant de recomposer avec les étoiles (avec le L extreme seulement. les difficultés d'assemblage m'ont dissuadées de tenter la mosaïque sans filtre pour les étoiles). Ca donne ceci comme résultat final. Il reste pas mal de petits défauts : le bas de l'image n'est pas très propre, le bleu est un peu granuleux dans les zones de faible intensité, un artefact ovale bleu à gauche au milieu d'origine indéterminée, la position du curseur noir compatible avec les plus faibles intensités détectées ? mais sinon elle supporte un peu le grandissement. Qu'en pensez-vous ?

Bien vu, Oui trop de précipitation à la fin pour les derniers ajustements de la luminosité du jpeg. entre les stretchs et stretchs inverses réalisés trop vite sur l'image finale plutôt que la nébuleuse, les étoiles peuvent en pâtir et paraître ternes. Mea culpa :-)

C'est vrai. en HOO on amplifie le bleu avec un coefficient arbitraire (ici j'ai pris en gros Ha=R et OIII=2.5(V+B) avec de menues corrections pour tenir compte des résidus de Ha dans les canaux V et B ou du OIII dans le rouge). En tenant compte des sensibilités des photosites RVB du capteur en Ha et OIII, on doit avoir sur cette nébuleuse en moyenne intensité(Halpha) de l'ordre de 3*Intensité(OIII). donc normalement c'est plutôt rouge partout. sinon c'est vrai que j'aurais aussi pu laisser davantage de dégradé de couleurs dans les zones riches ou pauvres en OIII. C'était un peu plus fade à mon gout en terme de contraste chromatique.

et finalement une version plus large avec la nébuleuse du Pélican (4 panneaux).12,1kpixels x 5,6kpixels ça commence à consommer beaucoup de RAM sur mon vieil ordi.

Après le pélican, sa voisine, la nébuleuse de l'America en HOO Newton SW 200/1000 + correcteur de coma + filtre L-extreme (bandes Ha et OIII) Caméra ASI 2600 couleur refroidie à -10°C (compliqué lorsque la température extérieure était de 27°C à 23h). guidage phd2, prise de vue APT Mosaïque de deux panneaux. 1h40 de pose sur chaque. Le champ d'étoiles est celui obtenu avec le filtre. Après le Pélican, j'anticipe cette fois la question "Mais où elle est l'Amérique du Nord ?" avec une version annotée:

Une version avec les étoiles retravaillées. Plus rondes (traitement python sur le starmask). En gros, je calcule une version moyenne locale des étoiles StarMoy ; ça supprime les étoiles voisines et conserve les aigrettes et les défauts (coma, astigmatisme,...) car présents sur toutes les étoiles du champ. puis je calcule starmoyrot, médiane des valeurs (Starmoy,rotation de pi/2 de Starmoy, rotation de pi..., rotation de 3pi/2... centré sur chaque étoile (coordonnées obtenues à partir du fichier stars.lst de dynamic psf de Siril) ; si les aigrettes n'ont pas de défauts (angle précisément de pi/2 entre chaque), ça supprime les défauts (coma,...) mais conserve les aigrettes. Il suffit ensuite de faire etoiles = etoiles-starmoy+starmoyrot pour obtenir un starmask modifié sans trop d'artefacts qui conserve, aux défauts près largement atténués, les caractéristiques du premier starmask. voilà le résultat. bon c'est pas encore partageable (trop de paramètres à régler manuellement et le traitement dure une heure sur cette mosaïque).

Bonsoir, La canicule a un bon côté. La turbulence du ciel est faible, le vent nul , bref de bonnes conditions pour l'astrophotographie. J'en ai profité ces derniers jours pour photographier le pélican sur 3 nuits, à quelques jours de la pleine lune. Nuit 1 en RVB (en fait filtre IDAS NGS1 que j'aurais pu remplacer par un filtre infrarouge...) pour les étoiles et tester la mosaïque de deux tuiles prévue. Nuit 2 et 3, filtre L-extrême, 2 tuiles avec 20% de recouvrement et 90x2' de poses conservées pour chaque tuile après sélection des meilleures images. Le setup : Newton Skywatcher 200/1000 + correcteur de coma+filtre+camera ASI 2600 couleur Le traitement du HOO. SIRIL presque exclusivement. Empilement des images des deux tuiles des nuits 2 et 3 , retrait de gradient très léger (polynôme d'ordre 1) en faisant très attention à la sélection des points puis réalisation de la mosaïque et alignement avec étoiles RVB. Réduction du bruit sur l'image nébuleuse+étoiles prise avec le L-extrême (GraXpert) et enfin un coup de starnet pour récupérer l'image de la nébuleuse seule (polluée par les halos des deux étoiles brillantes). création de l'image HOO (Halpha=R et OIII=k*((B+V)-.05*R) avec k=2.5 qui permettait d'équilibrer à peu près les deux histogrammes. L'image obtenue est un peu fade,le rouge et le bleu étant présents dans des proportions assez voisines j'ai augmenté le contraste des régions Halpha et OIII de façon non linéaire et également modifié la balance des couleurs de ces deux zones pour avoir un OIII plutôt turquoise et un Halpha plus rouge/orange/brun que rouge vif. Un tout petit peu d'augmentation de contraste par ondelettes et par le filtre CLAHE. Dernière étape ajout des étoiles RVB avec un peu de cosmétique autour des deux étoiles brillantes en essayant de respecter leurs couleurs naturelles. Malheureusement la montée d'histogramme des étoiles révèle une forme un peu triangulaire peut-être synonyme d'un pincement du primaire (à cause des hautes températures nocturnes 28 en début de nuit, 24 à 4h du matin). Ce défaut est noyé dans la halo des étoiles les plus brillantes et ne se voit pas pour les étoiles les plus faibles.A part faire une resynthèse intégrale des étoiles (plus d'aigrettes en ce cas ) à laquelle je ne me résous pas, je n'ai pas trouvé de solution simple pour faire des corrections cosmétiques de la forme des étoiles; la déconvolution ne fonctionne pas ici, Pixel math est insuffisant, Il faudrait utiliser Python mais ça risque de prendre du temps... Voici la mosaïque finale et un zoom sur la tête et l'oeil du pélican

Après M16, M17. Même setup (NewtonSW 200/1000+ correcteur de coma+filtre antipollution IDAS NGS1) . Cette fois 60' de poses retenues sur 3h possibles (les arbres, poteaux EDF) et une meilleure FWHM 3,3" au lieu de 4,3". Même traitement GHS + déconvolution pour les étoiles et la nébuleuse avec une légère réduction d'étoiles (pas trop c'est beau les étoiles surtout dans la voie lactée et aussi ça masque un peu les défauts dans la nébuleuse). Puis correction des défauts de la déconvolution et mixage étoiles/nébuleuse en respectant les couleurs originelles des étoiles (obtenues sur l'image brute après calibration photométrique des couleurs); ça évite (un peu) d'avoir des étoiles blanches sur les zones lumineuses de la nébuleuse ce qui arrive avec l'agorithme SIRIL de recomposition des étoiles. Il reste sûrement des défauts. Déjà peut-être un artefact (traitement ou reflet aux coordonnées 1/7e longueur 1/4 hauteur en partant de 0,0 en haut à gauche une tâche blanchâtre autour d'une étoile ?). Merci de vos commentaires :-)

J'y pense, j'y pense mais j'attends la retraite et un setup un peu plus léger. Je me plains de la pollution mais ça a de bons cotés de pouvoir aller se coucher dans un lit à 2 ou 3h du matin et de récupérer le matériel au petit matin :-)

Merci pour vos gentils commentaires, voici une version avec un peu plus de saturation (j'avais réduit celle des étoiles donc j'ai un peu de marge) et un fond de ciel un peu réduit. C'est à la marge. Et deux artefacts de traitements jaune vert à droite de l'image éliminés!

voilà la version corrigée. plus de grumeaux (deconvolution limitée aux zones lumineuses de la nebuleuse), plus de cernes autour des etoiles et une meilleure transition vers l'image de la nebuleuse. et des couleurs pour les etoiles sur la nebuleuse ce que la recomposition native dans Siril ne fait pas (les etoiles blanchissent lorsque le fond de ciel=la nebuleuse est trop intense et chromatique).

Merci pour les remarques. Je préfère la 3eme à la 2ème mais je suis d'accord avec toi. la texture grumeleuse et les étoiles ne sont pas top. Je repars au traitement :-) J'ai utilisé l'outil déconvolution RL mode multiplicatif de Siril, le seul que je parviens à faire fonctionner à peu près et ça crée ces structures inesthétiques sur l'image et ca complique aussi la transition starless/starmask en particulier sur les aigrettes. Pour les grumeaux je pense avoir trouvé la solution et peut être aussi pour les étoiles au moins partiellement. Nouvelle image dans quelques jours...

12 juin. pas de lune avant la fin de nuit. 56 poses de 30'' utiles (sur 2h) sur M16 entre les arbres ou les poteaux en direction du centre ville d'Aix en Provence. Matériel Newton SW 200/1000+ HEq5pro +correcteur de coma+filtre IDAS NGS1 + camera ASI 2600C. Suivi médiocre par phd2 (sigma moyen 0,9") probablement parce qu'on est proche de l'équateur, avec dithering fwhm des brutes 4,3'' pas terrible du tout. Soit la turbulence (forte près de l'horizon) soit la collimation à revérifier soit la MAP mal réglée ou qui se dérègle en cours de nuit Traitement Siril dark+flat (bien soigné, obligé à cause du voile de pollution lumineuse aussi intense que les parties brillantes de la nébuleuse) puis retrait de gradient GraXpert, couleur par photométrie pour obtenir la seconde image présentée enfin denoising avec graXpert, déconvolution pour réduire un peu la taille des étoiles, un peu, beaucoup, de travail chromatique pour redonner les bons couleurs sur les étoiles sur la nébuleuse qui seraient blanches sinon et reduire le décalage RGB des étoiles, variable selon la position sur le champ (le correcteur de coma ou le filtre qui sont responsables) et enfin la montée d'histogramme principalement avec la transformation hyperbolique généralisée et pour le starless et le starmask. Voilà les différentes étapes en image... Une brute pour commencer (instructif sur le fond de ciel!) l'image empilée après retrait de gradient et étalonnage colorimétrique et l'image finale après plus que 28' de traitement 😊 commentaires bienvenus.. c'est comme ça qu'on progresse

Je ne vois pas de soucis dans le traitement. Je suis surpris que le flat au crepuscule donne le même résultat. IL faut peut etre regarder un flat du crepuscule sans boite à flats. Tant pis il y aura un peu de gradient mais à part ca, ca doit donner la même structure que les images de ta galaxie parce que c'est la même chose en plus lumineux. Si ce n'est pas pareil j'achète un chapeau et je mange le chapeau. Si ca donne un bon flat alors c'est peut etre la boite à flat qui est responsable d'une façon ou d'une autre ou pas compatible avec ce setup particulier.

la boite a flat, c'est 4 led+plexi blanc + boite de melange de la lumiere en polystyrene+ nouvelle plaque de plexi blanc à l'entrée du tube (? celle là me semble necessaire) Avec ton starizona tu réduis la focale de ton instrument et tu collectes des photons qui arrivent sur le capteur avec un angle plus élevé que sans capteur. En gros sans capteur ton flat aurait une taille moitié .7*.7 donc avec forcement un vignettage moins marqué (genre 16500ADU au bord /19500ADU au centre au lieu de 15000/19500). Donc les defauts sont plus visibles. SOnt-ils vraiment absents sans le starizona (avant retrait de gradient) ? Oui GraxPert corrige partiellement ton flat par une parabole centrée au centre de ton capteur. Comme ton image /flat est plat au centre mais pas la parabole, il surcorrige au centre pour donner cette tache sombre qui a la forme d'une grosse poussiere. Je pense que ça n'a rien à voir avec le secondaire. Le decalage temporel entre le flat et les images c'est ce qui explique la difference entre les cartes de poussieres (pas toutes les poussieres seulement une petite partie je dirais 3 ou 4 max) et effectivement ca n'a rien à voir avec le vignettage qui est le (seul) vrai souci. Pour corriger partiellement. Tu tapes les commandes suivantes load master_flat offset -2000 save master_flat_corr ... et tu reprends ton traitement depuis le debut Sinon tu fais un flat à la lumiere du crepuscule+ plaque de plexi blanc à l'entrée du tube et tu l'utilises dans ta calibration. Si ca fonctionne (à part les nouvelles poussieres qui vont rester) je te donnerai un truc pour ne pas avoir à le faire à chaque fois

J'ai regardé les fichiers.. Ils sont faits comment les flats ? sur le ciel, boite a flat,... ? et le traitement flat-bias ? offset ? Il y a un décalage temporel important entre le flat et les images ? C'est la premiere fois que ca fait ça la correction par flat ? C'est peut-être lié au starizona. Sil n'y pas pas de milieu diffuseur efficace pour ton flat et que ton rayonnement n'est pas isotrope ? mais je pense dans ce cas que le flat aurait du surcorriger et te donner sur les images pp_ une tache centrale sombre et c'est l'inverse. premier point, certaines poussieres se sont décalées ou sont apparues entre le flat et les images bon mais sur ce type de champ des corrections cosmetiques ca devrait tres bien passer s'il n'y a pas de galaxies à l'endroit des poussieres. Sinon plus important le vignettage residuel apres traitement.... on peut sauver ton image en remplacant ton pp_flat par pp_flat offset -2000 et en reprenant le traitement DOF C'est une valeur un peu bizarre que tu n'as probablement pas choisie lorsque tu as fait tes flats. Ca pourrait correspondre à une lumiere parasite de l'ordre de 10% de ce que tu as imagé pour obtenir tes flats (c'est bizarre à moins d'avoir égaré une lampe de poche dans le C11 !) Restera a traiter les poussieres...

Bonjour, On a l'impression qu'il n'a simplement pas divisé par le flat pour obtenir l'image empilée. Je ne connais pas le C11 mais ce pallier assez plat au centre, confirmé par ce pic dans l'histogramme, semble anormal. Il ne semble pas y avoir de souci de MAP les poussieres ont la même taille sur l'image finale et le flat. Avant toute chose chose je regarderai une coupe sur une brute si tu retrouves ce rapport de luminosité 15000/20000 ADU soit aux bords 75% de la luminosité du centre puis la même chose sur l'image empilée parce que ca ressemble à un probleme d'offset... Sinon si tu nous donnes l'image fit (par ex) d'une brute ou plutot brute - dark, de l'image empilée et du flat je veux bien regarder en detail. Ce sont des images monochromes, tu peux les donner même en bin 4 ou 8 ca suffira et c'est moins gourmand en espace disque...

En pleine journée on ne voit pas grand chose sur l'écran des IFN de l'iris. Voici une version réduite où ils sont amplifiés

Bonjour, Je photographie depuis mon domicile, ciel de ville Bortle 7 (un peu mieux avec le filtre antipollution IDAS NGS1). Je pensais naïvement auparavant que je pourrais obtenir des résultats à peu près identiques à ceux qui travaillent sur des cieux moins pollués en augmentant seulement le temps de pose. par exemple pour avoir le même résultat que sur un ciel Bortle 4 avec un fond de ciel de 2 points de magnitude en plus, je devrais poser 2,512^2 =6.3 fois plus longtemps mais en fait ça ne fonctionne pas au delà d'un certain point. Voici deux images M94 et NGC 7023 pour illustrer mon propos. La première faite avec une ASI 585 refroidie (Newton 200/1000+correcteur de coma+ filtre IDAS précité, 300x30" soit 2h30 ) et NGC 7023 sur ASI 2600 MC refroidie (même setup 660x30" soit 5h30). Je ne suis pas mécontent du résultat mais en ce qui concerne les IFN (le dernier anneau autour de M94 ou les extensions de poussieres de la nébuleuse de l'Iris, je remarque des erreurs basse fréquence spatiale (taille caractéristique 500 pixels ou 5 à 10' d'arc ce correspond en gros à l'échelle 6 des ondelettes) que je ne parviens pas à corriger en posant plus longtemps. Je suis à peu près certain que le problème provient de petites imprécisions basse fréquence dans le flat et/ou le calcul de gradient qui pollué par ces erreurs de flat devient impossible à faire proprement. J'ai fait quelques calculs, avec un fond de ciel de 19 mag/arcsec2, celui ci représente 200 fois l'intensité caractéristiques de IFN que je veux observer(magnitude 24 à 25/arcsec2). Quand je regarde attentivement mes flats (réalisés avec une tablette led recouverte d'une feuille de papier canson, le tout collé sur le télescope) je vois apparaitre des défauts chromatiques basse fréquence qui ne représente qu'à peine 1% (ca semble peu mais en fait...!) de l'intensité du flat (pour les observer ci joint le rapport Flat_R/flat_V qui présente une structure de bandes, vue en mode histogramme, d'intensité relative 1%). Je n'ai pas d'explication certaine (le filtre ?) pour ces bandes mais juste quelques pistes pour améliorer la prise de flats (tablette+ zone de mélange de 4cm puis plaque plexi 4mm opale coulée ou bien revenir aux flats du crepuscule ou de l'aube avec milieu diffuseur : tee shirt, ou plaque diffusante à l'entrée du télescope). Seulement voilà 1% de l'intensité du flat ça fait une erreur sur le fond de ciel d'environ 2 fois l'intensité des IFN à une fréquence spatiale qui est celle des IFN justement donc difficile à extraire avec "backgroud extraction". En faisant au mieux (GraXpert mode IA se plante dessus sur les 2 cibles, je choisis le polynôme d'ordre 2 ou le mode RBF en choisissant tres soigneusement les points pour le calcul du gradient) Ca permet d'observer un résultat plausible mais partiellement incorrect si je compare avec ce qu'obtiennent les meilleurs astrophotographies de ces cibles. Pour info la partie variable du background (les termes d'ordre 1 et supérieur du polynôme d'interpolation représente 10 fois l'intensité des IFN) donc je dois l'estimer avec au plus 10% d'erreur. Donc avec ce ciel Bortle 7, je suis obligé d'atteindre une précision meilleure que 1% sur le flat, et d'environ 10% sur le background pour parvenir à estimer correctement les BF spatiales de ces IFN. Sous un ciel Bortle 4, les contraintes seraient 6,3 fois moins importantes c'est à dire que la tolérance sur la qualité des flats ou du calcul du background est nettement moins sévère. Un temps de pose supérieur (unitaire ou surtout global) ne change rien au delà d'une certaine durée sauf à améliorer les details haute fréquence de ces IFN. Si quelqu'un a une bonne idée sur la question je suis preneur

waouh 😍! Si j'arrive au niveau de la première version je serai déjà satisfait!

oups une erreur dans le premier script. Il faut enlever iadd starless0 à la fin du script sinon l'effet est divisé par 2... Ca donne ça après les 2 scripts...

sur cet exemple 1) en mode auto ajustement l'image starlesse de depart puis apres reduction du fond du ciel puis apres blanchiment du don de ciel ainsi que le masque utilisé (qui englobe toute la galaxie et ses plus faibles extensions)? Le résultat final est plus propre mais s'il y a du signal oublié hors du masque, il disparait partiellement... ce serait le cas de l'aigrette ici ou de petites galaxies mal repérées (j'ai fait vite)

@ant1 https://drive.google.com/file/d/1XXEWTNI_WVewPsGkiPQkSnDQhQPu4e7-/view?usp=sharing https://drive.google.com/file/d/1D9PfykhJWahME-II9NZbc-GxaK6xvLGp/view?usp=sharing Ce sont deux scripts que j'ai écrits. Il ne sont pas clé en main. Ils fonctionnent sur l'image starless. Je détaille Le premier réduction bruit de fond est efficace pour réduire le bruit à grande échelle (taches de 100 pixels ou plus le niveau 6 de ondelette). Il faut préparer dans le dossier de travail un dossier process, un fichier starless0 et un mask noir et blanc appelé masknb (1 zone non traitée donc toutes les structures où il y a du vrai signal,0 zone traitée) il réduit la composante 6 de l'ondelette d'un facteur 4 mais progressivement et de façon continue Le second nécessite le fichier starless0, le dossier process et de changer dans le script .ssf deux valeurs numériques. seuil1 (correspond en gros au quartile 10% du bruit de fond. En dessous de cette valeur tout est blanc( r=v=b) et seuil2 en gros l'intensité des plus faibles structures du signal. Au dessus de cette valeur les valeurs rvb sont inchangées. Entre les deux valeurs le chromatisme varie continûment... Ca peut etre adapté selon la structure de votre bruit en incluant les composantes 5 ou 4 et 5 ou ... du fond de ciel