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

[Cyril Cavadore]

J'ai l'image.
Quels parametres ?
Y a aussi 61 millions de pixels ... Il y a des chances que je passe cela en multicore.

[Cyril Cavadore]

J'arrive à calibrer ton image en astrometrie avec le GSC-ACT et en mettant 250 / 250 etoiles de matching.
C'est une tres grosse image ...

[Patrice Le Guen]

Pour ma part je n'y arrive pas avec ces paramètres. Mais en faisant un filtre gaussien à 1,5 ça passe.

[Cyril Cavadore]

Ah ca c'est interessant....
tu peux me montrer tes paramètres d'astrométrie STP.

[Patrice Le Guen]

Voici mes paramètres pour la fonction CTRL+U :

Et les résultats :

Code : Tout sélectionner

Débute le processus dans le "thread" :  14572
Extrait du catalogue les étoiles référence :  RA= 19h39m49s DEC=+00°48'25'' pixels de 3.8x3.8µm  Focale~ 530mm Mag_Min=-6.0 Mag_Max=30.0
Extrait du catalogue : 4633 étoiles de référence, catalogue = GSC-ACT, bande magnitude Magnitude bande V, Projection Sphérique
Astro Params. XC=4788, XC=3194, LX=9576, LY=6388, RA_ref=5.14792, Dec_Ref=-0.01408, RayX=69846.283, RayY=69846.3
Extrait les étoiles de l'image... (1ere passe)
Extraction rapide des étoiles terminée, 23553 trouvées, débute l'extraction avancée...
Nombre d'étoiles extraites de l'image, premiere passe 23553, seconde passe 19875, temps de calcul : 6.82 sec
Extraction des étoiles de l'image CCD
#      X      Y       Mg     Alpha =          Delta =        -> Carte du ciel                      
1      5961.2 4882.6  6.21   19h37m53.441s    +01°29'56.92''                                       
2      73.7   3883.5  6.80   19h47m33.268s    +01°05'21.30''                                       
      
249    2601.2 6317.7  11.23  19h41m37.300s    +02°05'15.86''                                       
250    3649.7 2293.7  11.24  19h41m41.059s    +00°26'14.14''                                       
 
#        Xccd         Yccd      Flux             -> Image CCD                                      
1        39.7         6279.1    2955977.84                                                         
2        8791.0       1773.0    2819684.29                                                         
                       
249      5526.3       4404.3    1753766.52                                                         
250      5002.4       1593.9    1750259.77                                                         
Nbre d'étoiles (les plus brillantes) limité à 250 pour la comparaison champ-catalogue.
Compare les étoiles de l'image et de la carte du ciel...
XCCD       YCCD        Xcat        YCat        CCD_Index Cat_index nbVotes RA_Cat      DEC_Cat     Flux_Image  Magnitude
4525.18    2937.05     5171.07     2550.34     203       219       451     19h39m11s   +00°32'33''  1858423.90  11.15 
5045.10    2898.97     5475.58     2415.95     112       96        429     19h38m41s   +00°29'15''  2069025.59  10.29 


4910.20    2154.21     2669.03     5196.24     117       134       115     19h43m18s   +01°37'40''  2064083.13  10.73 
6896.62    3095.53     2697.30     5581.89     248       148       111     19h43m15s   +01°47'09''  1754055.94  10.83 
Nombre d'étoiles pour le calcul du premier polynome de passage : 250
TransFormation "Trigmatch" de passage pour validation des centroides, degre : 1, étoiles utilisées : 5, basé sur résidu max en pixel : 1.0
Verification avec centroides de : 9 pixels de demi-largeur
Nombre d'étoiles à tester au centroide : 250
XCentro    Ycentro Centro OK ?                                                                     
5395       3661    OK  xc=5394.5 yc=3659.8  (19x19) ->  psf= 2.9x2.4 pixel                         
5681       3829    Echec du centroide                                                              
5600       3774    Echec du centroide                                                              
                                                          
3166       2332    Echec du centroide                                                              
3201       2341    Echec du centroide                                                              
15 étoiles après test, 250 avant test de centroide...
Nombre d'étoiles validées par centroide : 15
Astrométrie -> Calculs finaux, étape 1
Erreur interne -> Erreur sur le calcul du polynome de passage -> Pas assez d'étoiles, diminuer le degré du polynôme !
L'astrométrie a échouée !
-Terminé-

Avec degré du polynôme à 3 au lieu de 4

Code : Tout sélectionner

Débute le processus dans le "thread" :  12516
Extrait du catalogue les étoiles référence :  RA= 19h39m49s DEC=+00°48'25'' pixels de 3.8x3.8µm  Focale~ 530mm Mag_Min=-6.0 Mag_Max=30.0
Extrait du catalogue : 4633 étoiles de référence, catalogue = GSC-ACT, bande magnitude Magnitude bande V, Projection Sphérique
Astro Params. XC=4788, XC=3194, LX=9576, LY=6388, RA_ref=5.14792, Dec_Ref=-0.01408, RayX=69846.283, RayY=69846.3
Extrait les étoiles de l'image... (1ere passe)
Extraction rapide des étoiles terminée, 100000 trouvées, débute l'extraction avancée...
Nombre d'étoiles extraites de l'image, premiere passe 100000, seconde passe 49097, temps de calcul : 38.63 sec
Extraction des étoiles de l'image CCD
#      X      Y       Mg     Alpha =          Delta =        -> Carte du ciel                      
1      5961.2 4882.6  6.21   19h37m53.441s    +01°29'56.92''                                       
2      73.7   3883.5  6.80   19h47m33.268s    +01°05'21.30''                                       
                                     
247    2672.7 2634.0  11.23  19h43m17.242s    +00°34'36.37''                                       
248    9068.7 635.6   11.23  19h32m47.569s    -00°14'35.56''                                       
249    2601.2 6317.7  11.23  19h41m37.300s    +02°05'15.86''                                       
250    3649.7 2293.7  11.24  19h41m41.059s    +00°26'14.14''                                       
 
#        Xccd         Yccd      Flux             -> Image CCD                                      
1        8809.3       1407.0    2392111.13                                                         
2        1097.1       6239.8    2336709.39                                                         
3        8396.2       3111.1    2319754.43                                                         
                                                    
248      1730.6       1418.1    1437715.12                                                         
249      2236.5       1760.1    1437228.13                                                         
250      6369.1       5623.8    1436561.37                                                         
Nbre d'étoiles (les plus brillantes) limité à 250 pour la comparaison champ-catalogue.
Compare les étoiles de l'image et de la carte du ciel...
XCCD       YCCD        Xcat        YCat        CCD_Index Cat_index nbVotes RA_Cat      DEC_Cat     Flux_Image  Magnitude
5244.75    2845.35     5171.07     2550.34     155       219       423     19h39m11s   +00°32'33''  1676298.00  11.15 
5381.22    2713.87     5381.02     2733.79     206       223       405     19h39m45s   +00°37'05''  1548244.78  11.16 
4457.42    3571.31     4285.09     3246.74     17        176       392     19h42m48s   +00°49'41''  2192809.60  10.95 
4525.18    2937.05     4284.92     3246.93     81        112       391     19h42m48s   +00°49'42''  1858423.90  10.53 
5045.10    2898.97     5416.42     2704.61     35        40        387     19h38m47s   +00°36'21''  2069025.59  9.27 

2486.15    4615.69     7829.64     4291.03     220       66        119     19h34m49s   +01°15'23''  1516336.86  9.91 
2483.20    1327.43     2890.92     5484.46     94        207       118     19h42m56s   +01°44'45''  1822824.03  11.11 
Nombre d'étoiles pour le calcul du premier polynome de passage : 250
Reconnaissance des listes d'étoiles impossible !
L'astrométrie a échouée !
-Terminé-

[Marc DELDEM]

Salut Patrice,

Voici ce que j'arrive à faire avec ton image, avec l'UCAC4 et mes réglages pour mon setup.



Et voici le résultat:

Code : Tout sélectionner

Extrait du catalogue : 80464 étoiles de référence, catalogue = UCAC_4, bande magnitude Magnitude bande V, Projection Sphérique
Extrait les étoiles de l'image... (1ere passe)
Extraction rapide des étoiles terminée, 100000 trouvées, débute l'extraction avancée...
Nombre d'étoiles extraites de l'image, premiere passe 100000, seconde passe 48850, temps de calcul : 72.64 sec
Nbre d'étoiles (les plus brillantes) limité à 60 pour la comparaison champ-catalogue.
TransFormation "Trigmatch" de passage pour validation des centroides, degre : 1, étoiles utilisées : 6, basé sur résidu max en pixel : 0.5
47 étoiles après test, 60 avant test de centroide...
Degré : 6, nombre d'étoiles prises en compte : 47 
OK [44]  X=2332.512   Y=402.275  res.=0.013" dist=65%
OK [14]  X=5650.028   Y=128.180  res.=0.019" dist=55%
OK [16]  X=8530.662   Y=2416.086  res.=0.033" dist=66%
OK [46]  X=8109.060   Y=1336.142  res.=0.053" dist=66%
OK [32]  X=1526.560   Y=1025.836  res.=0.101" dist=68%
OK [6]  X=8517.185   Y=1350.214  res.=0.124" dist=72%
OK [31]  X=5690.872   Y=415.037  res.=0.126" dist=51%
OK [22]  X=607.242   Y=2630.463  res.=0.144" dist=73%
OK [38]  X=2382.057   Y=5281.438  res.=0.145" dist=55%
OK [34]  X=2324.738   Y=5057.520  res.=0.159" dist=54%
OK [21]  X=7662.051   Y=1411.934  res.=0.172" dist=59%
OK [15]  X=1557.966   Y=704.167  res.=0.226" dist=71%
OK [24]  X=4746.242   Y=6008.747  res.=0.257" dist=49%
OK [37]  X=3265.623   Y=5500.325  res.=0.267" dist=48%
OK [42]  X=8017.839   Y=1067.056  res.=0.275" dist=67%
OK [13]  X=5032.356   Y=408.550  res.=0.287" dist=49%
OK [26]  X=3617.664   Y=5526.472  res.=0.293" dist=45%
OK [33]  X=6849.372   Y=4025.777  res.=0.302" dist=39%
OK [11]  X=3569.414   Y=1009.145  res.=0.303" dist=43%
OK [4]  X=5493.821   Y=906.837  res.=0.315" dist=42%
OK [27]  X=4819.401   Y=5706.270  res.=0.351" dist=44%
OK [1]  X=4535.075   Y=3843.000  res.=0.358" dist=12%
OK [25]  X=3892.110   Y=1155.392  res.=0.358" dist=39%
OK [45]  X=5581.162   Y=3163.771  res.=0.367" dist=14%
OK [47]  X=6985.723   Y=923.374  res.=0.372" dist=55%
OK [17]  X=6108.250   Y=5575.969  res.=0.372" dist=47%
OK [40]  X=1605.621   Y=394.740  res.=0.391" dist=74%
OK [43]  X=4508.748   Y=1302.061  res.=0.444" dist=33%
OK [28]  X=4973.502   Y=6222.958  res.=0.464" dist=53%
OK [2]  X=3253.652   Y=4761.849  res.=0.492" dist=38%
OK [12]  X=6029.366   Y=3398.380  res.=0.500" dist=22%
OK [23]  X=7486.495   Y=4064.574  res.=0.544" dist=49%
OK [36]  X=1985.301   Y=5561.793  res.=0.591" dist=64%
OK [20]  X=2576.518   Y=1156.274  res.=0.663" dist=52%
OK [8]  X=6527.063   Y=5320.057  res.=0.679" dist=48%
OK [5]  X=6466.208   Y=2602.819  res.=0.708" dist=31%
OK [35]  X=7340.345   Y=4115.631  res.=0.729" dist=47%
OK [10]  X=4388.682   Y=1119.910  res.=0.730" dist=37%
OK [19]  X=3770.944   Y=2179.375  res.=0.834" dist=25%
OK [3]  X=3509.275   Y=3797.723  res.=0.876" dist=25%
OK [39]  X=7497.041   Y=2700.622  res.=0.966" dist=48%
OK [7]  X=8487.415   Y=5431.062  res.=1.141" dist=75%
OK [29]  X=2811.111   Y=2273.601  res.=1.296" dist=38%
OK [30]  X=191.035   Y=2698.073  res.=1.410" dist=80%
OK [9]  X=8954.149   Y=2158.021  res.=1.566" dist=75%
OK [41]  X=106.732   Y=735.638  res.=2.635" dist=92%
OK [18]  X=8397.709   Y=6096.917  res.=2.757" dist=80%
Référence -> distance max = 92%
Résidu : 0.801"
Nombre d'étoiles utilisées: 47, résidu final 0.558", Degré polynome: 2
Focale calculée : 536.0 mm
Angle axe Alpha par rapport à X+ CCD: 359.9 °
Angle axe Delta par rapport à X+ CCD: 270.0 °
Repère indirect !
Coordonnées du centre image, Alpha=20h51m40.197s, Delta=+30°48'29.23''
-Terminé-

C'est la version Prism 32 bits (v11.0.0.0).

Je ne sais pas si le résultat est très bon, les étoiles affichées en vert me semblent saturées avec un niveau max > 65565 ADU.

UPDATE: j'essaye de refaire la manipulation ci-dessus sur l'image d'origine, cela ne marche plus !

Marc

[Patrice Le Guen]

Bonjour Marc,

Avec tes paramètres et les 250 étoiles de comparaison du catalogue de Cyril, j'ai réussi l'étalonnage astrométrique de mon image d'origine :

Puis en en remettant les paramètres par défaut dans "Réglage de la reconnaissance auto", l'astrométrie a échoué. Ensuite j'ai remis 11 pour la magnitude min catalogue et 17 pour la magnitude max catalogue et l'astrométrie a réussi. Je pense que ceux sont ces 2 paramètres qui sont déterminants pour la réussite de l'astrométrie. Les valeurs par défaut -6 et 30 ne sont pas réalistes, il faudrait les changer. La magnitude limite sur terre c'est 22 je crois chez moi c'est 19, et les étoiles de magnitude 8 saturent en quelques secondes.

Code : Tout sélectionner

Débute le processus dans le "thread" :  17364
Extrait du catalogue les étoiles référence :  RA= 20h51m42s DEC=+30°48'33'' pixels de 3.8x3.8µm  Focale~ 536mm Mag_Min=11.0 Mag_Max=17.0
Extrait du catalogue : 78890 étoiles de référence, catalogue = UCAC_4, bande magnitude Magnitude bande V, Projection Sphérique
Astro Params. XC=4788, XC=3194, LX=9576, LY=6388, RA_ref=5.46157, Dec_Ref=-0.53772, RayX=70634.535, RayY=70634.5
Extrait les étoiles de l'image... (1ere passe)
Extraction rapide des étoiles terminée, 100000 trouvées, débute l'extraction avancée...
Nombre d'étoiles extraites de l'image, premiere passe 100000, seconde passe 48811, temps de calcul : 40.40 sec
Extraction des étoiles de l'image CCD
#      X      Y       Mg     Alpha =          Delta =        -> Carte du ciel                      
1      3534.1 4627.4  9.90   20h54m04.991s    +31°23'19.76''                                       
2      1060.8 3966.9  10.77  20h58m45.816s    +31°06'37.00''                                       
             
249    6583.5 5359.1  11.58  20h48m16.611s    +31°41'02.79''                                       
250    6586.2 5363.2  11.58  20h48m16.292s    +31°41'08.86''                                       
 
#        Xccd         Yccd      Flux             -> Image CCD                                      
1        9529.2       6331.8    2423903.05                                                         
2        419.1        1826.1    2420678.93                                                         
                                   
248      5283.0       4054.2    1446248.44                                                         
249      4308.4       4451.9    1444181.47                                                         
250      9312.7       1978.0    1441136.44                                                         
Nbre d'étoiles (les plus brillantes) limité à 250 pour la comparaison champ-catalogue.
Compare les étoiles de l'image et de la carte du ciel...
XCCD       YCCD        Xcat        YCat        CCD_Index Cat_index nbVotes RA_Cat      DEC_Cat     Flux_Image  Magnitude
5068.72    3684.08     4522.74     3681.36     69        140       1219    20h52m12s   +31°00'23''  1930119.71  11.42 
4991.19    4238.44     4600.02     4236.16     123       107       1073    20h52m03s   +31°13'53''  1761940.26  11.36 

5201.68    4415.14     3066.57     1064.45     112       30        37      20h54m55s   +29°56'33''  1811015.95  11.18 
3610.29    5914.62     5189.43     214.11      221       160       35      20h50m57s   +29°36'00''  1516721.03  11.47 
Nombre d'étoiles pour le calcul du premier polynome de passage : 250
TransFormation "Trigmatch" de passage pour validation des centroides, degre : 3, étoiles utilisées : 67, basé sur résidu max en pixel : 0.5
Verification avec centroides de : 9 pixels de demi-largeur
Nombre d'étoiles à tester au centroide : 250
XCentro    Ycentro Centro OK ?                                                                     
5070       3685    OK  xc=5068.7 yc=3684.1  (19x19) ->  psf= 5.4x4.8 pixel                         
4991       4239    OK  xc=4991.2 yc=4238.4  (19x19) ->  psf= 5.2x4.5 pixel                         

4563       6300    OK  xc=4561.6 yc=6298.5  (19x19) ->  psf= 5.0x4.3 pixel                         
6528       1070    OK  xc=6527.1 yc=1068.9  (19x19) ->  psf= 5.6x5.2 pixel                         
4406       217     OK  xc=4405.1 yc=216.2  (19x19) ->  psf= 5.6x5.0 pixel                          
227 étoiles après test, 250 avant test de centroide...
Nombre d'étoiles validées par centroide : 227
Astrométrie -> Calculs finaux, étape 1
Calcul automatique terminé
Astrométrie -> Calculs finaux, étape 2
Astrométrie réussie
SolAlpha (+1.6077510E+01,-8.7726652E-04,-9.9987662E-01,-2.0272174E-08,-2.0974548E-08,-8.7551462E-09)
SolDelta (-2.5243944E+00,+1.0000298E+00,-8.3719089E-04,+2.7524741E-08,-1.4179427E-08,+1.7568396E-10)
SolX     (+1.6077272E+01,-8.7770466E-04,-1.0001219E+00,+2.0237739E-08,-2.1019039E-08,-8.7307967E-09)
SolY     (+2.5377530E+00,+9.9996898E-01,-8.3764540E-04,+2.7565617E-08,+1.4153288E-08,-1.6825789E-10)
Nombre d'étoiles utilisées: 227, résidu final 0.585", Degré polynome: 2
Focale calculée : 536.0 mm
Angle axe Alpha par rapport à X+ CCD: 0.1 °
Angle axe Delta par rapport à X+ CCD: 90.0 °
Repère direct !
Coordonnées du centre image, Alpha=20h51m40.178s, Delta=+30°48'27.85''
-Terminé-

Mais le sujet qui me préoccupe c'est la saturation à 13107 étoiles de l'algorithme de calcul de la courbe RMS. Et la le problème reste entier.