BUG EtalonAstrometrieAuto

[Marc DELDEM]

Bonsoir,

Il y a un bug dans la fonction EtalonAstrometrieAuto (version Prism 10.3.47.402)

Voici un script pour reproduire le bug:
On charge une image brute, qui n'est pas calibrée astrométriquement.
On affiche les informations sur la focale et les coordonnées Alpha/Delta.
On applique la fonction EtalonAstrometrieAuto à l'image.
Si la résolution astrométrique fonctionne, on demande ensuite la focale et les coordonnées Alpha/Delta.
On voit que la focale a bien été mise à jour. Mais les coordonnées Alpha/Delta sont inchangées.

Le script utilise le catalogue UCAC4.

Code : Tout sélectionner

// Catalogs
cat[1]$="GSC-ACT"
cat[2]$="UCAC_4"
cat[3]$="NOMAD_V1"
cat[4]$="USNO_A2"
cat[5]$="NOMAD_V1_(Internet)"
cat[6]$="PPXML_(Internet)"
cat[7]$="GAIA_DR1_(Internet)"
cat[8]$="Hipparcos_Catalog"
cat[9]$="SAO_Catalog"
cat[10]$="BSC_Catalog"
cat[11]$="Tycho_2_Catalog"
catalog$=cat[2]$

// Selected Catalog directory
catalogDir$="M:\DOCUMENTS\ASTRONOMIE\CATALOGUES\UCAC4\u4b"

// Max number of stars in catalog
maxNumberStarCat = 60

// Max number of stars in image
maxNumberStarImg = 60

// Image North/South oriented ?
isImgNorthSouthOriented=FALSE

// Degree of polynomial correction of image deformation
degPoly=2

// Load image
GetImage img

// Retrieve image informations
GetSizePixX img pixelSizeX
GetSizePixY img pixelSizeY
GetBinX img binX 
GetBinY img binY
GetFocus img focalLength
GetAlpha img raImg
GetDelta img decImg
Str_Rad_To_Ra raImg HH$ MM$ SS$
Str_Rad_To_Dec decImg DD$ DM$ DS$
Print "Position (J2000) from image BEFORE standard astrometric resolution"
Print " RA : " HH$ "h " MM$ "m " SS$ "s and DEC : " DD$ "° " DM$ "' " DS$ "''"
Strf focalLength 1 focalLength$
Print " Focal length : " focalLength$ "mm"
Print

// Try Standard Astrometry
EtalonAstrometrieAuto img catalog$ focalLength pixelSizeX*binX pixelSizeY*binY raImg decImg maxNumberStarCat maxNumberStarImg isImgNorthSouthOriented degPoly catalogDir$ isPlateSolveOk
IF isPlateSolveOk=1 THEN
	GetFocus img focalLength
	GetAlpha img raImg
	GetDelta img decImg
	Str_Rad_To_Ra raImg HH$ MM$ SS$
	Str_Rad_To_Dec decImg DD$ DM$ DS$
	Print "Position (J2000) from image AFTER standard astrometric resolution"
	Print " RA : " HH$ "h " MM$ "m " SS$ "s and DEC : " DD$ "° " DM$ "' " DS$ "''"
	Strf focalLength 1 focalLength$
	Print " Focal length : " focalLength$ "mm"
	Print
ELSE
	Print "Astrometric resolution has failed"
ENDIF

Close img

Voici un lien sur le script et sur l'image qui me sert pour le test: lien dropbox

Le résultat du script:

Position (J2000) from image BEFORE standard astrometric resolution
RA : 18h 41m 51.038s and DEC : +37° 52' 26.40''
Focal length : 1640.0mm

Position (J2000) from image AFTER standard astrometric resolution
RA : 18h 41m 51.038s and DEC : +37° 52' 26.40''
Focal length : 1650.5mm

Script terminé ...

=> On voit que la focale a bien été ajustée mais pas alpha/delta

Si je fais un étalonnage astrométrique depuis Prism, voici le résultat après calibration:


A noter que la fonction script Find_Allsky_Position fonctionne et modifie bien les coordonnées alpha/delta de l'image.

Marc

[Cyril Cavadore]

La fonction script ne MET pas a jour le RA/DEC et c'est volontaire/historique.

[Marc DELDEM]

Pas très logique mais OK.

Est-t'il envisageable d'avoir une version qui sauvegarde dans l'image le résultat de l'astrométrie ?

Si on fait l'astrométrie via FIND_ALLSKY_POSITION, cela fonctionne.
En terme de précision sur l'astrométrie par rapport à EtalonAstrometrieAuto, on a quoi comme différence ?

Merci.
Marc

[Cyril Cavadore]

C'est simplement historique
Je te dirais de vive voix, la je peux pas l'ecrire