Photométrie sur une série d’images impossible d’afficher la courbe RMS

L'analyse des images.
Patrice Le Guen
Messages : 281
Enregistré le : 19 janv. 2014, 18:06
Localisation : Toulon
Contact :

Photométrie sur une série d’images impossible d’afficher la courbe RMS

Message par Patrice Le Guen »

Bonjour,

Je me suis abonné aux alertes de AstroColibri pour détecter la partie visible des sursauts gamma. Je me suis donc tourné vers la fonction « Photométrie sur une série d’images » pour tenter de détecter la variation d’un point dans mes images.
Malheureusement le bouton RMS de la fonction « Photométrie sur une série d’images » ne fonctionne pas avec la caméra ASI6200MM-Pro en haute résolution 9576x6388 couplée à ma lunette FSQ106 de 530mm de focale et donnant un champs de 3,9° x 2,6° à partir d’un certain nombre d’étoiles valides.

En effet, après calculs, lorsque l’on clique sur le bouton RMS une boite de dialogue affiche « Range check error » et la fenêtre du graphique RMS s’ouvre mais n’affiche pas la courbe RMS en fonction de la magnitude, l’écran reste blanc sans afficher de points.
Range check error.jpg
Après avoir fait des essais avec différentes images je n’ai pas trouvé de corrélation avec les dimensions de l’image ni le nombre de pixels dans l’image ni le catalogue d’étoiles (GSC, UCAC, GAIA DR1...).
A chaque essai j’ai noté le nombre d’étoiles valides. On peut voir d’après ce tableau que si le nombre d’étoiles valides est supérieur à 13107 une erreur de type « Range Check error » est provoquée lorsque l’on clique sur le bouton RMS.

Code : Tout sélectionner

Objet	Résolution X	Résolution Y	Nb pixels	Etoiles valides	Catalogue	Range check error lors de lancement de RMS
GRB		5000	5000		25000000	11846		UCAC4		Non
GRB		5250	5250		27562500	13114		UCAC4		Oui
GRB		5125	5125		26265625	12455		UCAC4		Non
GRB		5187	5187		26904969	12790		UCAC4		Non
GRB		5218	5218		27227524	12965		UCAC4		Non
GRB		5234	5234		27394756	13038		UCAC4		Non
GRB		5242	5242		27478564	13081		UCAC4		Non
GRB		5246	5246		27520516	13096		UCAC4		Non
GRB		5248	5248		27541504	13102		UCAC4		Non
GRB		5249	5249		27552001	13102		UCAC4		Non
GRB		5250	5250		27562500	13111		UCAC4		Oui
GRB		5250	5249		27557250	13171		UCAC4		Oui
GRB		5250	5248		27552000	13113		UCAC4		Oui
GRB		9576	2877		27550152	12719		UCAC4		Non
GRB		9576	3000		28728000	13267		UCAC4		Oui
GRB		4788	3194		15292872	18661		UCAC4		Oui
C2021A1		9576	6388		61171488	2558		UCAC4		Non
GRB		9576	6388		61171488	18689		UCAC4		Oui
AMERICA		9576	6388		61171488	26819		UCAC4		Oui
AMERICA		6000	6000		36000000	14909		UCAC4		Oui
AMERICA		5900	5900		34810000	14335		UCAC4		Oui
AMERICA		5500	5500		30250000	12162		UCAC4		Non
AMERICA		5600	5600		31360000	12677		UCAC4		Non
AMERICA		5700	5700		32490000	13224		UCAC4		Oui
AMERICA		5650	5650		31922500	12949		UCAC4		Non
AMERICA		5650	5650		31922500	12949		UCAC4		Non
AMERICA		5675	5675		32205625	13079		UCAC4		Non
AMERICA		5687	5687		32341969	13143		UCAC4		Oui
AMERICA		5681	5681		32273761	13104		UCAC4		Non
AMERICA		5684	5684		32307856	13128		UCAC4		Oui
AMERICA		5682	5682		32285124	13115		UCAC4		Oui
AMERICA		5681	5682		32279442	13111		UCAC4		Oui
AMERICA		5681	5682		32279442	13109		UCAC4		Oui
AMERICA		5681	5682		32279442	13107		UCAC4		Non
AMERICA		5681	5682		32279442	13108		UCAC4		Oui
Lors des traitements PRISM utilise 15 % du CPU et 2Go. Il semble donc que la fonction « Photométrie sur une série d’images » soit limitée à 13107 étoiles. Pour contourner ce problème je suis donc obligé de fenêtrer mes images en 2 afin de ne pas saturer PRISM en étoiles, ce qui a pour inconvénient de rallonger les délais de traitement.

Serait-il possible de corriger ce problème de limite du nombre d’étoiles valides ?


PS : Ma configuration (setup) du 31/12/2021
• PRISM 64 bits V11.0.3.6 du 29/12/2021
• Windows 11 Famille en 64 bits – version 21H2 - version OS 22000.376
• Processeur Intel(R) Core(TM) i7-9700KF CPU @ 3.60GHz
• Disque dur SSD 128Go TS128GMSA370
• Mémoire ram installée 32Go
• Résolution 2560x1440
• ASCOM Platform 6.5 SP1 – RC1 6.5.1.3174
• Lunette TAKAHASHI FSQ-106ED
• Monture allemande EM200 alimentée en 24V Temma 2Jr et driver ASCOM Telescope driver for Takahashi Temma protocols 6.2.6
• Caméra principale Acquisition ZWO ASI6200MM-PRO pilote Windows 3.16.0.0 du 30/12/2020. Pilote logiciel PRISM 1.05.04 du 20 Mars 2021
• Caméra secondaire Autoguider SBIG ST-i pilote Windows 2.41.02 du 20/01/2011 pulse guide par ASCOM
• Contrôle de focalisation : USB-Focuser 2 SkyCenterStepperv1.0.52 du 27/09/2021 x64 (Pierro-Astro) + moteur Robofocus pilote FTDI 2.12.28.0 du 16/08/2017
• Gestion d’interrupteurs IP9258FR Release Version : V4.01 2013/04/18
• Serveur série MOXA NP5210A Firmware version=0x1020000
• Onduleur FSP NANO 600
• Modem routeur NETGEAR DGN2200v4 Firmware V1.0.0.58_1.0.58
Modifié en dernier par Patrice Le Guen le 08 janv. 2022, 23:59, modifié 1 fois.
Marc DELDEM
Messages : 1264
Enregistré le : 20 juin 2009, 01:01

Re: Photométrie sur une série d’images impossible d’afficher la courbe RMS

Message par Marc DELDEM »

Bonjour Patrice et meilleurs vœux,

Il faudra avoir un retour de Cyril.
Chez moi avec un petit capteur (ST8) pas de soucis.

Dans la fenêtre texte, sur l'image de référence, il t'indique combien d'étoiles au niveau de l'extraction ?
Par exemple:
Image : t_TCP J07094936+1412280-V-001.cpa, debut extraction étoiles sur X1=1 Y1=1 X2=1529 Y2=1020 ...
Image : t_TCP J07094936+1412280-V-001.cpa, 15606 étoiles extraites (passe n°1)
Image : t_TCP J07094936+1412280-V-001.cpa, 883 étoiles extraites (passe n°2)

Autre question: au niveau FWHM, tu as quoi comme valeur moyenne en pixels ? Tu es autour d'1 pixel ou largement supérieur à 1 pixel. Car à mon avis, dans le traitement qu'effectue Prism, il doit surement chercher à éliminer les pixels chauds/froids. Et si ta FWHM est trop petite, cela doit poser un problème je pense.

Bon après, ton analyse montre effectivement un problème sur une valeur seuil. Il y aura peut-être des modifications nécessaires dans Prism pour bien gérer les gros capteurs / petite focale que l'on peut rencontrer maintenant.

Marc
Patrice Le Guen
Messages : 281
Enregistré le : 19 janv. 2014, 18:06
Localisation : Toulon
Contact :

Re: Photométrie sur une série d’images impossible d’afficher la courbe RMS

Message par Patrice Le Guen »

Bonjour Marc,

Bonne année 2022.

Voici un extrait des 27Mo du fichier LOG.

Code : Tout sélectionner

[Main]: TU 2022-01-01  21h30m42.986s {1.1 ms} : Débute le processus dans le "thread" :  17200
 [Main]: TU 2022-01-01  21h30m42.987s {1.3 ms} : Extrait du catalogue les étoiles référence :  RA= 20h56m02s DEC=+44°22'21'' pixels de 3.8x3.8µm  Focale~ 536mm Mag_Min=8.0 Mag_Max=17.0
 [Main]: TU 2022-01-01  21h30m43.121s {133.3 ms} : Extrait du catalogue : 27078 étoiles de référence, catalogue = UCAC_4, bande magnitude Magnitude bande V, Projection Sphérique
 [Main]: TU 2022-01-01  21h30m43.542s {419.1 ms} : Astro Params. XC=2840, XC=2837, LX=5681, LY=5674, RA_ref=5.48048, Dec_Ref=-0.77444, RayX=70603.886, RayY=70603.9
 [Main]: TU 2022-01-01  21h30m43.543s {1.0 ms} : Extrait les étoiles de l'image... (1ere passe)
 [Main]: TU 2022-01-01  21h30m52.539s {8964.6 ms} : Extraction rapide des étoiles terminée, 9190 trouvées, débute l'extraction avancée...
 [Main]: TU 2022-01-01  21h30m54.981s {2434.1 ms} : Nombre d'étoiles extraites de l'image, premiere passe 9190, seconde passe 5692, temps de calcul : 2.38 sec
 [Main]: TU 2022-01-01  21h30m56.109s {1124.0 ms} : Extraction des étoiles de l'image CCD
 [Main]: TU 2022-01-01  21h30m56.110s {846 µs} : #      X      Y       Mg     Alpha =          Delta =        -> Carte du ciel                      
 [Main]: TU 2022-01-01  21h30m56.111s {779 µs} : 1      2862.0 2841.9  6.82   20h55m59.017s    +44°22'26.26'' 
 [Main]: TU 2022-01-01  21h30m56.947s {658 µs} : Nombre d'étoiles pour le calcul du premier polynome de passage : 51
@[17200]: TU 2022-01-01  21h30m56.948s {616 µs} : [TrigMach] Entering TrigMatch
 [17200]: TU 2022-01-01  21h30m56.948s {602 µs} : [TrigMach] pass #1 MaxErrorPixel=0.5
 [17200]: TU 2022-01-01  21h30m56.949s {628 µs} : [TrigMach] AmountPoints=51 DimPoly=4            
  [17200]: TU 2022-01-01  21h30m57.001s {650 µs} : Polynom dimension : 4; Amount of stars  : 15/51, residuals errorMax (pixels) : 0.50 
 [17200]: TU 2022-01-01  21h30m57.002s {636 µs} : [TrigMach] pass #2 MaxErrorPixel=0.5
 [17200]: TU 2022-01-01  21h30m57.003s {739 µs} : [TrigMach] AmountPoints=51 DimPoly=3
  [17200]: TU 2022-01-01  21h30m57.038s {657 µs} : Polynom dimension : 3; Amount of stars  : 26/51, residuals errorMax (pixels) : 0.50 
 [17200]: TU 2022-01-01  21h30m57.038s {582 µs} : [TrigMach] pass #3 MaxErrorPixel=0.5
 [17200]: TU 2022-01-01  21h30m57.039s {665 µs} : [TrigMach] AmountPoints=51 DimPoly=2
  [17200]: TU 2022-01-01  21h30m57.074s {622 µs} : Polynom dimension : 2; Amount of stars  : 25/51, residuals errorMax (pixels) : 0.50 
 [17200]: TU 2022-01-01  21h30m57.074s {668 µs} : [TrigMach] End, DimPolynome=2, NbStar=25
@[Main]: TU 2022-01-01  21h30m57.075s {899 µs} : TransFormation "Trigmatch" de passage pour validation des centroides, degre : 2, étoiles utilisées : 25, basé sur résidu max en pixel : 0.5
 [Main]: TU 2022-01-01  21h30m57.076s {827 µs} : Verification avec centroides de : 6 pixels de demi-largeur
 [Main]: TU 2022-01-01  21h30m57.077s {802 µs} : Nombre d'étoiles à tester au centroide : 120
 [Main]: TU 2022-01-01  21h30m57.078s {797 µs} : XCentro    Ycentro Centro OK ?                                                                     
 [Main]: TU 2022-01-01  21h30m57.080s {2.2 ms} : 2661       2936    OK  xc=2660.8 yc=2936.2  (13x13) ->  psf= 5.5x5.1 
 
                           
Il indique première passe 9190 étoiles, seconde passe 5692 étoiles. C'est bizarre que ces nombres soient inférieurs au plafond de 13107 étoiles que j'ai trouvé, mais je viens de refaire le traitement : c'est bien ça.

Les images on été pré-traitées et alignées avant la photométrie sur une série d'images. La FWHM est comprise entre 3 et 4, de plus d'après le module d'astrométrie il faut qu'elle soit inférieur à 5, donc je suis dans la bonne fourchette.
Dans les paramètres de reconnaissance auto j'ai indiqué 5 pour le paramètre de la force du filtre passe bas pour l'extraction des étoiles afin d'éviter à l'algorithme de prendre le bruit de fond pour des étoiles et de réduire au maximum le nombre d'étoiles extraites sinon j'arrivai à plus 18000 étoiles.
Paramètres reconnaissance auto.jpg
J'ai fait des tests sur plusieurs objets (AMERICA, C/2021A1, un GRB (Gamma Ray Burst)) j'obtiens le même plafond à 13107 étoiles maximum géré par la fonction RMS.

Ci-dessous d'autres copies d'écrans de mes paramètres :
photometris serie images.jpg
option astrometrie.jpg
La caméra ASI6200MM-PRO est une caméra peu commune avec ses 62 méga pixels et ses dimensions 24x36 mm qui couplée avec une lunette de 530 mm de focale donne un champs immense et donc on arrive facilement à des nombres d'étoiles très important qui explosent facilement les limites de PRISM. Je pense que cette caméra est un outil qui peut permettre d'identifier des GRB grâce à l'algorithme de photométrie sur une série d'image de PRISM. Ça serait une bonne publicité pour PRISM si j'arrivais à en identifier un par cette méthode, mais avec cette limite de 13107 étoiles c'est fastidieux.
D SAINTE-CATHERINE
Messages : 42
Enregistré le : 22 mars 2020, 22:12

Re: Photométrie sur une série d’images impossible d’afficher la courbe RMS

Message par D SAINTE-CATHERINE »

J'ai vu sur la 1ère image vous étiez en Prism 32 bits.
Avez-vous testé en 64 bits ?
Patrice Le Guen
Messages : 281
Enregistré le : 19 janv. 2014, 18:06
Localisation : Toulon
Contact :

Re: Photométrie sur une série d’images impossible d’afficher la courbe RMS

Message par Patrice Le Guen »

Bonjour D SAINTE-CATHERINE ,

C'était une copie d'écran de mes tests au début en 32 bits.

Je suis passé en 64 bits mais le résultat est le même, en plus dans l'exemple ci-dessous j'arrive à trouver plus de 26 000 étoiles.
Photometrie 64 bits.jpg
Marc DELDEM
Messages : 1264
Enregistré le : 20 juin 2009, 01:01

Re: Photométrie sur une série d’images impossible d’afficher la courbe RMS

Message par Marc DELDEM »

Patrice,

Est-ce que tu as essayé avec les images prétraitées uniquement, sans faire un quelconque alignement au-préalable ?

Sinon, un truc que je ne comprends pas dans ton traitement. Tu indiques que ta FWHM est comprise entre 3 et 4. C'est quoi, c'est en pixels ou en secondes d'arc ?
Car ensuite dans le paramétrage de la photométrie, tu définis un rayon de photométrie étoile de 2.5 pixels.
En photométrie d'ouverture, le rayon de mesure pour l'étoile, on doit le prendre en gros autour de 1.5 fois la FWHM. Si tu as une FWHM de 4 pixels, ce rayon est donc à fixer à une valeur de 6 pixels.

Et sinon, en dehors de cette fonction, si tu prends la première image de la séquence, tu la calibres en photométrie (Analyse / Photométrie / Etalonnage Automatique) puis une fois étalonnée, tu fais Analyse / Photométrie / Mesurer/Extraire toutes les étoiles, est-ce que cela fonctionne ?

Marc
Patrice Le Guen
Messages : 281
Enregistré le : 19 janv. 2014, 18:06
Localisation : Toulon
Contact :

Re: Photométrie sur une série d’images impossible d’afficher la courbe RMS

Message par Patrice Le Guen »

Bonjour Marc,

La FWHM que j’indiquais comprise entre 3 et 4 est en pixels mais il se peut que j’ai mesuré des étoiles dans un coin de l’image, et cela ne reflète pas la FWHM moyenne qui est de 2,9 pixel en X et 2,2 pixels en Y.
Si l’on revient à mes images et que je mesure (clic droit centroïde) une étoile non-saturée et au centre de l’image j’obtiens une valeur de FwhmL = 2,75p (FWHM horizontale en pixels). Si je fais une coupe horizontale de l’étoile j’obtiens un profil avec un maximum à 60000 environ, donc une mi-hauteur à 30000 environ. Sur le graphique la largeur va de 9,7 à 12,5 environ soit une largeur de 2,8 environ. Ce qui est cohérent avec la mesure du centroïde.
FWHM
FWHM
Tu écris « En photométrie d'ouverture, le rayon de mesure pour l'étoile, on doit le prendre en gros autour de 1.5 fois la FWHM ». Je prends donc 2,75*1,5=4,12 arrondi à 4.

En prenant ce paramètre de 4 pixels comme rayon d’ouverture pour la photométrie et lorsque je lance la photométrie sur mes images non-alignées alors l’astrométrie échoue. Alors qu’avec les images déjà recalées ça fonctionne avec 4 pixels de rayon d'ouverture.
Echec astrometrie sur images non alignees.jpg
Mais j'ai remarqué que la valeur de la focale doit être plus précise que ce qui est demandé (10%). Si je lance une astrométrie (CTRL+U) sur mon image de NGC7000 avec focale = 530 mm l'astrométrie échoue. Si je recommence avec focale = 537 mm, l'astrométrie est un succès. Donc sur la focale une erreur de plus de 1% n'est pas permise. Or l'algorithme de photométrie sur une série d'images ne demande pas de saisir la focale, il va lire la valeur de la focale dans chacun des entêtes de fichier FIT. Avec la fonction "Mise à jour d'entêtes par lot" je modifie cette valeur focale à 537 au lieu de 530 dans tous mes entêtes et ça fonctionne avec 26585 étoiles trouvées. Malheureusement le clic sur le bouton RMS génère encore une erreur "Range Check Error".
x9a7vuuh5o.png
Photometrie sur images pre-traitees.jpg
Pour ce qui est de l’extraction de toutes les étoiles.
Je prends la première image de ma séquence sortie du pré-traitement et je fais Analyse / Photométrie / Etalonnage Automatique. Puis je fais Analyse / Photométrie / Mesurer/Extraire toutes les étoiles. Cette fonction marche bien sur mes images.
Mesure de toutes les étoiles
Mesure de toutes les étoiles
Marc DELDEM
Messages : 1264
Enregistré le : 20 juin 2009, 01:01

Re: Photométrie sur une série d’images impossible d’afficher la courbe RMS

Message par Marc DELDEM »

Salut Patrice,

Rien à dire sur ton analyse, il faudra que Cyril se plonge sur cette problématique.

Marc
Cyril Cavadore
Messages : 4481
Enregistré le : 18 juin 2009, 22:41
Localisation : 42

Re: Photométrie sur une série d’images impossible d’afficher la courbe RMS

Message par Cyril Cavadore »

Oui, il semble qu'il y ai un pb avec tes grosses images.
J'ai pourtant que des structures dynamiques.
Essaie de faire une astrometrie simple dessus pour voir.
Cyril Cavadore, ALCOR-SYSTEM, Auteur de PRISM
http://www.alcor-system.com

Pour reporter vos problème, la marche à suivre est ici
viewtopic.php?f=7&t=1403
Patrice Le Guen
Messages : 281
Enregistré le : 19 janv. 2014, 18:06
Localisation : Toulon
Contact :

Re: Photométrie sur une série d’images impossible d’afficher la courbe RMS

Message par Patrice Le Guen »

L'astrométrie simple a du mal à fonctionner dans les champs stellaires très denses, je viens de transférer une image ici :
https://transfert.free.fr/xSCMel
Répondre