Mise a jour donnee astrometrique dans l'entete de l'image

[stephane_guisard]

Bonjour a tous,

J'ai une image au format fits que je soumets à astrometry.net via prism
En retour j'obtiens les coordonnées du centre du champs, la focale et l'orientation du champs par rapport au N/S.

.... jusque la tout va bien.

1) Je sauve l'image au format .fits pour que les infos astrometriques calculees soient sauvees dans l'entete de l'image. Je la ferme, puis la recharge ... : les coordonnées alpha/beta sont bien celles trouvées par astrometry.net. Par contre l'angle de rotation m'indique "0" et la focale n'est donnée qu'a une décimale.
2) si je sauve l'image au format cpa, la ferme puis la rouvre, la par contre toutes les coordonnées + rotation du champs se rechargent avec les valeurs correctes et la focale est donnée avec les 2 décimales

Est ce que le format fits ne permet pas de sauver l'info de la rotation du champs ou bien un petit bug de prism ?
Serait-il possible de ne pas perdre de décimales dans la focale quand on sauve en fits (non, je ne pinaille pas .... quand on travaille avec de courtes focales et que l'on veut faire de la mosaïque derrière, j'ai besoin de garder mes 2 décimales de focale).
Je vois aussi que la valeur alpha pour limage fits est legerement differente de l'originale (au depart je n'ai fait qu'une seule analyse de l'image avec astrometry.net, c'est le resultat que j'ai ensuite sauvé sous les 2 formats fits et cpa ....donc a priori les valeurs étaient stritement les mêmes)
Merci,

Stephane

Je joins 3 copies d’écran :
1-entete image directement apres retour de l'info depuis astrometry.net
2-entete image apres avoir sauve l'image en .fits et l'avoir reouverte
3-entete fits apres avoir sauve l'image en .cpa et l'avoir reouverte

[Cyril Cavadore]

Je sauve l'image au format .fits pour que les infos astrometriques calculees soient sauvees dans l'entete de l'image. Je la ferme, puis la recharge ... : les coordonnées alpha/beta sont bien celles trouvées par astrometry.net. Par contre l'angle de rotation m'indique "0" et la focale n'est donnée qu'a une décimale.

le Fits n'a pas de mot clefs specifique pour cela et ca fait ce que ca peut.

2) si je sauve l'image au format cpa, la ferme puis la rouvre, la par contre toutes les coordonnées + rotation du champs se rechargent avec les valeurs correctes et la focale est donnée avec les 2 décimales

Oui, en CPA je fais ce que je veux

Est ce que le format fits ne permet pas de sauver l'info de la rotation du champs ou bien un petit bug de prism ?

Ca y est peut etre pas. Si tu prends un autre logiciel, aucun ne relira un mot clef que je vais devoir inventer pour cela, donc je ne vois pas l'interet du FITS

Serait-il possible de ne pas perdre de décimales dans la focale quand on sauve en fits (non, je ne pinaille pas .... quand on travaille avec de courtes focales et que l'on veut faire de la mosaïque derrière, j'ai besoin de garder mes 2 décimales de focale).

Pourquoi faire ?
La focale est en INT dans les FITS, et je ne vois aucun logiciel la relire que ce soit en INT ou en DOUBLE, puisque les mots CLEFS en FITS il n'y en a que deux qui sont reconnus par tous. Et pas la focale en tout cas.

Je vois aussi que la valeur alpha pour limage fits est legerement differente de l'originale (au depart je n'ai fait qu'une seule analyse de l'image avec astrometry.net, c'est le resultat que j'ai ensuite sauvé sous les 2 formats fits et cpa ....donc a priori les valeurs étaient stritement les mêmes)

Ca je ne sais pas pourquoi, mais il ne faut pas prendre la valeur de Astrometry.Net comme precise. Je deconseille tres fortement de faire un modele de pointage avec astrometry.net, j'y ai vu des erreurs de plusieurs dizaines de secondes d'arc par rapport a une astrometrie classique dans Prism, dont j'ai totalement confiance.

[stephane_guisard]

Merci Cyril,

Mais je ne comprends pas alors : 99.9% des images astro sont en format FITS, et les professionnels et amateurs font des réduction astrométriques sur ces images ... comment est sauvegradée alors, dans les header des images fits, l'information de rotation du champs du CCD par rapport au N/S par exemple ?

Pour la focale, tu dis que cette valer est en INT dans les headers fits ..pourtant prism l'y écrit avec une décimale (au lieu des 2 que je voudrais comme il le fait pour le format cpa).

Merci pour tes éclaircissemenrt,

A+

Stéphane

[Cyril Cavadore]

Il n'y a que peu de mot clefs FITS universels, decrit dans la norme, les autres sont des mots clefs personnalisés, que les observatoire ou logiciels creent et qui ne seront pas lisible par d'autres nécessairement.
C'est pour cela que j'avais fait le CPA il y a 25 ans.

Ici les mots standards
01 : SIMPLE = T /FITS header
02 : BITPIX = -32 /No.Bits per pixel
03 : NAXIS = 2 /No.dimensions
04 : NAXIS1 = 3072 /Length X axis
05 : NAXIS2 = 2048 /Length Y axis
06 : DATE = '2012-06-01T16:35:41.98' /Date of FITS file creation
07 : ORIGIN = 'Observatoire des Pises' /Origin place of FITS image
14 : XPIXELSZ= 9.0 /X pixel size microns
15 : CDELTM1 = 9.00000000000000E-003 /Size of a pixel (x) in mm
16 : YPIXELSZ= 9.0 /Y pixel size microns
17 : CDELTM2 = 9.00000000000000E-003 /Size of a pixel (y) in mm
35 : POSTN-RA= 1.14783480067876E+002 /Telescope RA
36 : POSTN-DE= 2.64829287335506E+001 /Telescope DEC
37 : CRVAL1 = 1.14783480067876E+002 /approx coord. in RA
38 : CRVAL2 = 2.64829287335506E+001 /approx coord. in DEC
39 : CDELT1 = 4.09090909088627E-006 /ScaleX in rad/pix
40 : CDELT2 = 4.09090909088627E-006 /ScaleY in rad/pix
41 : DATAMAX = 2384.0 /Maximum data value
42 : DATAMIN = 1249.0 /Minimum data value


et ceux que j'ai créé ->

08 : OBSERVER= 'JML ' /Observer name
09 : SITELAT = '+44:02:30:99 ' /Latitude observatory
10 : SITELONG= '-03:30:18:00 ' /Longitude observatory
11 : DATE-OBS= '2012-03-24T22:09:24.81' /UT date of Observation
12 : UT = '22:09:24.81 ' /UT begining of observation
13 : EXPOSURE= 300.0 /Exposure time seconds
18 : INSTRUME= 'KAF-6303E ' /Camera which created data
19 : BINX = 1 /X binning
20 : BINY = 1 /Y binning
21 : FILTERS = 'SFILTRE ' /Filters
22 : MIRORY = T /Y miror applied to image
23 : X1 = 17 /X1 image windowing
24 : Y1 = 2 /Y1 image windowing
25 : X2 = 3072 /X2 image windowing
26 : Y2 = 2048 /Y2 image windowing
27 : COMPRESS= 0 /Number of compression
28 : TELESCOP= 'T500 Newton ' /Telescop
29 : FOCAL = 2200.0 /Focal length in mm
30 : DIAMETER= 500.0 /Diameter in mm
31 : OBJCTRA = '07 39 08.035 ' /Object Right Ascension (J2000)
32 : OBJCTDEC= '+26 28 58.54 ' /Object Declinaison (J2000)
33 : RA = 1.14783480067876E+002 /Telescope RA
34 : DEC = 2.64829287335506E+001 /Telescope DEC
43 : CCD-TEMP= -3.52008300781250E+001 /CCD temperature as Celsius degrees
44 : END

Il doit faire 8 caracters max, rien n'empeche de faire

TOTO1=
TOTO2=

ect...
le FITS c'est un vrai bazard.

Pour la rotation il n'y a aucun mot clef standard.