Champ NGC 7728 - identification

[Pierre Valvin 261]

Hello !

Comme je l'évoquais dans un post précédent, j’aime bien faire des champs profonds de galaxies pour m’amuser à regarder ce qu’on y trouve. J’en ai fait un cette semaine, dans Pégase. NGC7728.

48X180s (2h24) en L sur mon T250 F/3.6, avec une Altair 2600MM (1x1,5°).

 

 

On peut utiliser Aladin pour explorer le champ, mais on peut aussi s’amuser avec les bibliothèques python (astropy, astroquery et photutils) pour un rendu, des tris et des infos plus personnalisés. C’est hyper intéressant et pas si difficile à utiliser.

Du coup, je ne suis pas certain d’être dans le bon forum, mais je me suis dis que ça pourrait peut-être intéresser certain.es. Je vous propose donc ces deux cartes.

La première a été faite de la manière suivante :

1/ Détection automatique des sources (ce qui est au-dessus de 1,3x le bruit)

>threshold = 1.3 * bkg.background_rms

>segment_map = detect_sources(convolved_data, threshold, npixels=20)
>segm_deblend = deblend_sources(convolved_data, segment_map, npixels=25, nlevels=32, contrast=0.20, progress_bar=False)

 

2/ Téléchargement de toutes les sources du catalogue des quasars sur simbad qui sont dans mon champ.

>r_quasars = Vizier.query_region(co, radius=demi_diag * u.deg, catalog='VII/289/superset')

 

3/ s'il y a une correspondance de coordonnées (différence inférieure à 0.005° = 18’’) et si le redshift est supérieur à 2, j’affiche un cercle aux coordonnées non pas de ma détection, mais bien de la source du catalogue. J’affiche à côté du cercle ce que je veux piquer dans le catalogue. Ici, j’ai mis le redshift (z).

>for i in range(len(r_detect)):
>    oc=pixel_to_skycoord(r_detect['xcentroid'], r_detect['ycentroid'], wcs)
>    c = SkyCoord(ra=oc.ra.degree*u.degree, dec=oc.dec.degree*u.degree) 
>    idx2, d2d2, d3d2 = match_coordinates_sky(c, catalog_quasars) 

>    if d2d2.degree<0.005 :
>        if r_quasars[0]['z'][idx2]>2:
>            ax.scatter(catalog_quasars[idx2].ra.degree, catalog_quasars[idx2].dec.degree,
>                   transform=ax.get_transform('fk5'), s=100, edgecolor='blue', facecolor='none')
>            ax.text(catalog_quasars[idx2].ra.degree+0.01, catalog_quasars[idx2].dec.degree+0.0, r_quasars[0]['z'][idx2], transform=ax.get_transform('fk5'), fontsize=4)

 

J’ai 68 quasars qui correspondent à ces critères !

 

 

La seconde carte, c’est la même chose avec les galaxies du PGC. J’ai affiché l’idendificateur.

 

 

On peut faire pareil avec tous les catalogues conservés au CDS de Strasbourg et afficher l’entrée du catalogue qu’on souhaite. Distance, axe, type, magnitude.

Voilà, j'espère que ce petit pas de côté ne vous a pas trop (méga)bassiné..

Il ne reste qu'à rêver sur ces horizons incommensurables !

 

Pierre

Modifié 21 août 2023 par Petzval

[Dihydropyridine 256]

Superbe champ, bravo et merci pour ces infos de traitement de signal !

 

Je suis également en train d'imager cette même zone du ciel depuis plusieurs nuits avec un champ un peu plus grand, j'ai eu le vertige en effectuant une détection des objets du ciel profond, la notion de "rêve sur ces horizons incommensurables" est en effet bien réelle !

 

Nico

[Pierre Valvin 261]

Bonjour @Dihydropyridine,

Merci pour ton commentaire. J'espère que tu partageras ton image !

D'autant plus que tu sembles avoir plus de patience que moi pour poser longtemps.

J'en profite pour rappeler qu'il faut "ouvrir l'image dans un nouvel onglet" pour voir la pleine résolution.

Pierre

Modifié 21 août 2023 par Petzval

[XavS 4 532]

Beaucoup de monde sur cette prise Bravo

 

XavS

[ALAING 59 336]

L'image est belle Pierre

Le reste est un peu trop technique pour le contemplatif que je suis

Bonne journée,

AG

[candrzej 2 784]

Bravo pour cette lumiance et l explication complete pour la detection des objets lointaine

 

@+

Christian

[olivier 321]

Bravo, j'adore ce genre d'image profonde. On dirait qu'il y a des IFN dans ce champ, non ?

[Pierre Valvin 261]

Merci pour vos commentaires, @XavS, @candrzej.

 

@ALAING : Je me vois aussi comme un contemplatif ! J'irai même plus loin, les astronomes sont fondamentalement des contemplatifs : ils ne peuvent pas intervenir sur leur sujet d'étude, aux contraire d'autres discipline ou on peut agir dessus (physique du solide, biologie, chimie...).

Cet outil d'identification me permet juste de contempler encore plus ce que nous offre le ciel et notre prise d'image.

Et j'avoue que ça me permet de contempler aussi l'immensité du travail d'observation, de classification et de compréhension de mes congénères humains. Il y a 24093 catalogues dispos sur VizieR... C'est surement un lieu commun de dire ça, mais la somme de connaissances accumulées laisse sans voix.

 

@olivier : oui, si on tire sur les manettes, c'est superbe. Il aurait fallu poser plus longtemps.

Et j'ai encore pas mal de choses à apprendre en traitement pour savoir faire sortir ça sans cramer le reste....

 

 

[danielo 5 951]

bravo pour cette belle image d'un champ bien riche et pour la manip !

[Pierre Valvin 261]

Merci @danielo pour ton passage.

Tu es un peu spécialiste (avec d'autres) des images très profondes sur lesquelles on aime se perdre aussi.

Et je sais que tu les annotes parfois aussi.. avec quel outil d'ailleurs ?

[Bruno- 3 995]

Passionnant !

 

Pourquoi ce choix de NGC 7728 ? Elle appartient à un amas de galaxies ?

[Pierre Valvin 261]

@Bruno-Merci !

J'ai simplement cherché quelques zones riches dans Pégase avec CdC (j'observe à l'est, quand ça monte). J'ai sélectionné ce champ qui remplissait bien mon capteur avec quelques NGC. J'avais d'autres candidats qui seraient à faire si je tournais ma caméra de 90° (confirmés avec l'outil de cadrage de Nina). En effet, j'ai une petite manie : j'aime bien que les aigrettes soient dans l'axe des pixels. Je trouve ça plus beau. Du coup, on est à 45° des axes des axes N-S et E-O...

 

[Bruno- 3 995]

il y a 3 minutes, Petzval a dit :

J'ai simplement cherché quelques zones riches dans Pégase avec CdC

 

Ah OK. Effectivement, il y a des coins de Pégase qui sont très riches. La preuve en images !

 

(S'il fait beau ce soir, je viens de préparer deux cartes de Pégase avec plein de galaxies pour tester mon Dobson 495 mm, c'est un peu le même genre de manipe mais en visuel. J'ai choisi l'amas Pegasus I et une zone qui a l'air riche sur l'atlas, un peu au nord.)

 

[danielo 5 951]

Il y a 1 heure, Petzval a dit :

Et je sais que tu les annotes parfois aussi.. avec quel outil d'ailleurs ?

 

Avec Pixinsight, les quasars sont ceux du catalogue Miliquas. https://heasarc.gsfc.nasa.gov/W3Browse/all/milliquas.html

 

Je me demandais justement si tu n'en obtenais pas plus avec ta méthode ?

[COM423 13 495]

Un champ profond qui laisse effectivement rêveur

Très sympa cette manip de pouvoir annoter l'image de façon automatisée, même si je n'ai pas suivi comment on utilise en pratique ces bibliothèques Python.. (?)

Bon ciel.

Modifié 21 août 2023 par COM423

[Pierre Valvin 261]

@danielook.

Milliquas est aussi dispo sur VizieR (VII/290/catalog). Je viens de demander combien il a d'entrées sur mon champ. Milliquas a 289 entrées.

En regardant avec Aladin au début, il m'avait semblé que le meilleur catalogue pour les quasars était VII/289/superset. En effet, il me sort 556 sources. Je n'ai sélectionné pour l'affichage que les sources qui remplissaient les 3 critères suivants :

1/ je les détecte

2/ ma détection est à moins de 0.05° de ce que dit le catalogue

3/ le redshift est >2

Du coup, il n'en reste "que" 68.

 

@COM423 Merci de rêver aussi et du mot sympa.

J'ai donné juste quelques pistes pour les plus curieux.. Il y a quelques lignes de code à écrire en amont quand même (charger les bibliothèques, ouvrir et préparer l'image...), mais c'est vraiment pas méchant en s'aidant des docs hyper bien faites de astropy. C'est d'ailleurs toujours étonnant de voir qu'il suffit d'une ligne pour faire des choses assez velues (segmenter une image !) grâce à la mise à dispo de tout ce travail de la communauté.

Je peux t'envoyer le fichier Jupyter si tu veux essayer.

 

[COM423 13 495]

On 8/21/2023 at 10:51 PM, Petzval said:

Je peux t'envoyer le fichier Jupyter si tu veux essayer.

Merci beaucoup pour cette attention que je garde de côté pour l'instant : j'ai bien fait un MOOC il y a quelque temps par curiosité sur Python, mais je n'ai pas pratiqué et ne connaîs absolument pas les bibliothèques que tu mentionnes.

Quand j'en trouverai le temps, je commencerais par regarder tout çà et reviendrai alors vers toi

 

Bonnes observations.

[Pierre Valvin 261]

@COM423

Ca marche, avec plaisir !

[Pierre Valvin 261]

@Dihydropyridine m'avait motivé pour essayer de faire les couches RGB (merci !).

J'ai essayé... Pas facile... Je suis pas sûr de tenir la bonne méthode (comment et à quel moment fait-on les étirements d'histogramme ?)

J'ai mis 15x180s par couche en plus des 48x180s de luminance.

 

Bonne soirée !

 

Pierre

Modifié 10 septembre 2023 par Petzval

[ALAING 59 336]

Je pense que l'on peut saturer un peu plus la couleur

Ce champ est vraiment de toute beauté

Bonne soirée,

AG

[danielo 5 951]

Le champ est vraiment profond et magnifique !

 

Il y a 3 heures, Petzval a dit :

J'ai essayé... Pas facile... Je suis pas sûr de tenir la bonne méthode (comment et à quel moment fait-on les étirements d'histogramme ?)

J'ai mis 15x180s par couche en plus des 48x180s de luminance.

 

45 min par couche à f/3.6 tu as normalement largement assez de données pour la couleur (perso je vais 20min à 30min par couche à f/4).

 

Concernant le traitement : je traite d'abord  la luminance toute seule comme si c'était une image mono, montée d'histo inclus. Ensuite je m'attaque aux couches RGB, retrait du gradient,  égalisation du FDC, calibration des couleurs par spectrophotométrie, un petit coup léger de déconvolution (je ne suis pas certain que ce soit indispendable), ensuite montée d'histogramme avec les mêmes paramètres que pour la luminance et enfin réduction de bruit. Une fois le traitement de la RVB terminé je lui applique la luminance (non linéaire aussi donc)  Je ne sais pas ce que tu utilises pour traiter mais avec Pixinsight tu as un curseur saturation que tu peux décaler vers la gauche pour l'augmenter.

 

 

Modifié 10 septembre 2023 par danielo

[Pierre Valvin 261]

@ALAINGMerci pour ta remarque. En effet, les couleurs manquent de peps. J'ai corrigé. Mais à la lecture de la réponse de @danielo, je devrais essayer de tout reprendre.

Danielo : ce post est précieux pour moi. Merci ! La procédure de pré-traitement est claire depuis longtemps pour moi, mais je n'ai jamais bien compris dans quel ordre faire les choses ensuite, notamment pour l'assemblage LRGB. J'ai l'impression de bricoler... Pas bon. Je devrais faire un stage pixinsight (je suis sous Siril..) ! Ton mail clarifie beaucoup les choses. Tu traites donc les RGB chacune séparément pour le gradient, le FDC et à ce moment tu les assembles pour la photométrie donc, et tu fais le reste. Bien sûr, pour les champs profonds, pas question de jouer avec StarNet, n'est-ce pas ?

J'ai un champ magnifique dans la boite. Je vais reprendre le traitement à la lumière de ce que tu dis et je le mettrai sur AS. J'espère avoir ton commentaire !

Merci mille fois !!