conversion de format de fichier cpa>fits

Questions relatives au traitement d'images
Répondre
francois Kugel
Messages : 14
Enregistré le : 26 nov. 2009, 08:00

conversion de format de fichier cpa>fits

Message par francois Kugel »

Bonjour,

Avec la v7.1.14, lorsque je convertit les cpa en fits (conversion de format de fichier par lots) l'en-tête fits semble perdre les coordonnées.

A+
Francois
Cyril Cavadore
Messages : 4480
Enregistré le : 18 juin 2009, 22:41
Localisation : 42

Re: conversion de format de fichier cpa>fits

Message par Cyril Cavadore »

Le probleme du format FITS et de tout le malheur qu'il nous apporte, c'est qu'il n y a AUCUN mot clef standard pour les coordonnées RA-DEC
Tu peux dans l'entete FITS des mots clef qui s'appellent "TARTEAUPOM'" et lui mettre n'importe quelle valeur.
Il n'y a aucune normalisation sur les mot clefs de l'entete FITS... ca fait bien 35 ans que ca dure....
Les astronomes n'ont rien fait pour normaliser les fichier et leur entetes.
Cyril
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
Marc DELDEM
Messages : 1264
Enregistré le : 20 juin 2009, 01:01

Re: conversion de format de fichier cpa>fits

Message par Marc DELDEM »

Bonjour,

Je rebondi à mon tour sur ce sujet.
Est-ce qu'il n'y a pas moyen de faire quelque chose ?
Le problème est qu'avec PRiSM et le format CPA, on arrive à avoir une image bien calibrée astrométriquement.
Mais ensuite, pour utiliser des logiciels externes (pour ma part, Aladin), il faut obligatoirement transformer l'image CPA en FITS et là on perd tout le bénéfice du travail sous PRiSM puisque l'on perd la partie astrométrique.

Bon je ne suis pas spécialiste du tout du sujet, mais j'ai vu qu'il existait le FITS World Coordinate System (WCS).
Ne serait-t-il pas possible d'utiliser ses conventions pour les entêtes FITS ?

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

Re: conversion de format de fichier cpa>fits

Message par Cyril Cavadore »

Je vais regarder...
A+
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
Cyril Cavadore
Messages : 4480
Enregistré le : 18 juin 2009, 22:41
Localisation : 42

Re: conversion de format de fichier cpa>fits

Message par Cyril Cavadore »

Salut,
J'ai charge un fichier PDF, de 46 pages ... c'est pas simple
En gros il faut ecrire ces mots clefs dans l'entete :

CRPIX1 = 384.0
CRPIX2 = 384.0
EQUINOX = 2000.0
CRVAL1 = 83.6331249995932
CRVAL2 = 22.0144722229022
CTYPE1 = 'RA---TAN'
CTYPE2 = 'DEC--TAN'
RADECSYS= FK5
CD1_1 = -2.79955404370653E-4
CD1_2 = -3.04062183279145E-6
CD2_1 = -3.04062183279145E-6
CD2_2 = 2.79955404370653E-4

Cyril
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
Marc DELDEM
Messages : 1264
Enregistré le : 20 juin 2009, 01:01

Re: conversion de format de fichier cpa>fits

Message par Marc DELDEM »

Je remonte ce sujet.

Avec la V8, si on fait l'astrométrie sur l'image CPA, on sauve en CPA puis on sauve en FITS, l'image en FITS n'a toujours pas l'astrométrie.
Je pensais que tu avais fais quelque chose pour ce sujet avec la V8.
C'est un oubli ?

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

Re: conversion de format de fichier cpa>fits

Message par Cyril Cavadore »

Ben non, ca devrait y etre.
As tu regardé dans l'entete FITS ?
Par contre, je ne pense pas que PRISM relise cette astrometrie, c'est fait pour d'autres logiciels qui gere le WCS.
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
Marc DELDEM
Messages : 1264
Enregistré le : 20 juin 2009, 01:01

Re: conversion de format de fichier cpa>fits

Message par Marc DELDEM »

Ce que je constate:

Je pars d'une image CPA dont la calibration astrométrique a été faite (on le vérifie dans la page d'information de l'image CPA, la case étalonnage astrométrique est cochée).
Je sauve l'image en FITS (enregistrer sous) et je ferme l'image. On ouvre l''image FITS précédente, quand on regarde la page d'information, n'a plus la calibration astrométrique (la case étalonnage astrométrique n'est plus cochée).

Si sur une image FITS, on fait la calibration astrométrique, alors dans la page information, la case étalonnage astrométrique est cochée.
Puis tu sauves l'image FITS (bouton réécrire sur le disque) et tu la réouvres. Alors dans la page information, la case étalonnage astrométrique n'est plus cochée.

L'entête FITS, après conversion de l'image CPA calibrée astrométriquement, est la suivante:
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 = '2011-03-27T09:56:59.61' /Date of FITS file creation
07 : ORIGIN = 'Observatoire des Pises' /Origin place of FITS image
08 : SITELAT = '+44:02:30:00 ' /Latitude observatory
09 : SITELONG= '-03:30:16:99 ' /Longitude observatory
10 : DATE-OBS= '2011-03-05T22:49:57.34' /UT date of Observation
11 : UT = '22:49:57.34 ' /UT begining of observation
12 : EXPOSURE= 3600.0 /Exposure time seconds
13 : XPIXELSZ= 9.0 /X pixel size microns
14 : CDELTM1 = 9.00000000000000E-003 /Size of a pixel (x) in mm
15 : YPIXELSZ= 9.0 /Y pixel size microns
16 : CDELTM2 = 9.00000000000000E-003 /Size of a pixel (y) in mm
17 : INSTRUME= 'KAF-6303E ' /Camera which created data
18 : BINX = 1 /X binning
19 : BINY = 1 /Y binning
20 : FILTERS = 'VIDE ' /Filters
21 : MIRORY = T /Y miror applied to image
22 : X1 = 18 /X1 image windowing
23 : Y1 = 3 /Y1 image windowing
24 : X2 = 3072 /X2 image windowing
25 : Y2 = 2048 /Y2 image windowing
26 : COMPRESS= 0 /Number of compression
27 : TELESCOP= 'T500 Newton ' /Telescop
28 : FOCAL = 2.19799350876276E+003 /Focal length in mm
29 : DIAMETER= 500.0 /Diameter in mm
30 : OBJCTRA = '12 37 12.435 ' /Object Right Ascension (J2000)
31 : OBJCTDEC= '+62 11 44.95 ' /Object Declinaison (J2000)
32 : RA = 1.89301813399222E+002 /Telescope RA
33 : DEC = 6.21958198459889E+001 /Telescope DEC
34 : POSTN-RA= 1.89301813399222E+002 /Telescope RA
35 : POSTN-DE= 6.21958198459889E+001 /Telescope DEC
36 : CRVAL1 = 1.89301813399222E+002 /approx coord. in RA
37 : CRVAL2 = 6.21958198459889E+001 /approx coord. in DEC
38 : CTYPE1 = 'RA---TAN ' / Projection type
39 : CTYPE2 = 'DEC---TAN ' / Projection type
40 : CUNIT1 = 'DEG ' / Projection type
41 : CUNIT2 = 'DEG ' / Projection type
42 : CRPIX1 = 1535 / Pixel
43 : CRPIX2 = 1023 / Pixel
44 : CRVAL1 = 1.89302315324390E+002 / o
45 : CRVAL2 = 1.89302315324390E+002 / o
46 : CD1_1 = -2.34527270076542E-004 / o/pixel
47 : CD1_2 = 0.0 / o/pixel
48 : CD2_1 = 0.0 / o/pixel
49 : CD2_2 = 2.34659998064456E-004 / o/pixel
50 : CDELT1 = 4.09464357563810E-006 /ScaleX in rad/pix
51 : CDELT2 = 4.09464357563810E-006 /ScaleY in rad/pix
52 : DATAMAX = 856.0 /Maximum data value
53 : DATAMIN = 668.0 /Minimum data value
54 : CCD-TEMP= -2.49273925781250E+001 /CCD temperature as Celsius degrees
55 : HISTORY C2,C6,U3
56 : END
Dans cette entête FITS:
09 : SITELONG= '-03:30:16:99 ' /Longitude observatory
Le signe négatif pour une longitude à l'Est du méridien de Greenwhich, cela me semble suspect non ? (les coordonnées sont celles des Pises).

Si j'ouvre l'image FITS dont l'entête est indiquée ci-dessus dans Aladin, les coordonnées du centre de l'image dans Aladin sont 00:37:30 et -09:17:47 alors que les vraies coordonnées sont 12:37:09 et +62:12:33, donc c'est dans les choux.

Marc
Marc DELDEM
Messages : 1264
Enregistré le : 20 juin 2009, 01:01

Re: conversion de format de fichier cpa>fits

Message par Marc DELDEM »

OK j'ai trouvé le problème.

Si on regarde sur le site http://fits.gsfc.nasa.gov/fits_wcs.html le document Representations of World Coordinates in FITS ici http://www.atnf.csiro.au/people/mcalabre/WCS/wcs.pdf.

On a:
CRVALia coordinate value at reference point (real floating)
CRPIX ja pixel coordinate of the reference point (real floating)
CUNITia units of CRVALia and CDELTia (character)
CDELTia coordinate increment (real floating)

On voit qu'il y a un problème sur:
36 : CRVAL1 = 1.89301813399222E+002 /approx coord. in RA
37 : CRVAL2 = 6.21958198459889E+001 /approx coord. in DEC
40 : CUNIT1 = 'DEG ' / Projection type
41 : CUNIT2 = 'DEG ' / Projection type
42 : CRPIX1 = 1535 / Pixel
43 : CRPIX2 = 1023 / Pixel
44 : CRVAL1 = 1.89302315324390E+002 / o
45 : CRVAL2 = 1.89302315324390E+002 / o
Les lignes 44 et 45 sont redondantes avec les lignes 36 et 37 et la ligne 45 est fausse.

Si je remplace les lignes 44 et 45 par un truc bidon sans rapport avec CRVAL1 et CRVAL2, cela fonctionne alors dans Aladin.

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

Re: conversion de format de fichier cpa>fits

Message par Cyril Cavadore »

Prism sauve le FITS en WCS, mais ne le lira pas, c'est normal , il faut que je regarde pourquoi mais a l'epoque j'avais une raison, je reviendrais la dessus plus tard.

Je viens essayer avec SAOimage avec un fichier FITS sauvé par Prism et avec l'astrometrie dedans.
Et ca marche ->

saoimage.JPG
Voici l'entete de cette image :

31 : OBJCTRA = '01 20 04.516 ' /Object Right Ascension (J2000)
32 : OBJCTDEC= '+03 25 33.65 ' /Object Declinaison (J2000)
33 : CRVAL1 = 2.00188192470577E+001 /approx coord. in RA
34 : CRVAL2 = 3.42601500263058E+000 /approx coord. in DEC
35 : CRPIX1 = 1092.0 /Centre pixel coordinate in X
36 : CRPIX2 = 736.0 /Centre pixel coordinate in Y
37 : CTYPE1 = 'RA---TAN ' / Projection type
38 : CTYPE2 = 'DEC---TAN ' / Projection type
39 : CUNIT1 = 'DEG ' / Projection type
40 : CUNIT2 = 'DEG ' / Projection type
41 : CRPIX1 = 1091 / Pixel
42 : CRPIX2 = 735 / Pixel

43 : CRVAL1 = 2.00189937964399E+001 / o
44 : CRVAL2 = 2.00189937964399E+001 / o

45 : CD1_1 = -1.73407976606442E-004 / o/pixel
46 : CD1_2 = 0.0 / o/pixel
47 : CD2_1 = 0.0 / o/pixel
48 : CD2_2 = 1.73405696752385E-004 / o/pixel
49 : CDELT1 = 3.02659806857912E-006 /ScaleX in rad/pix
50 : CDELT2 = 3.02659806857912E-006 /ScaleY in rad/pix


Je vois qu'il y a effectivement DEUX fois le couple CRVAL1 / CRVAL2, il faut donc que j'elimine ces deux CRVAL1/CRVAL2 qui sont de trop... je crois qu'on est d'accord ? Apparemment SAO image a l'air de s'en accommoder mais pas les autres logiciels.
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
Répondre