Page 1 sur 1
Image issue de Prism / Python / WCS
Posté : 16 févr. 2021, 00:19
par Marc DELDEM
Salut,
Est-ce que quelqu'un a réussi, à partir d'une image acquise dans Prism, à utiliser le package Python
Astropy et le module WCS (World Coordinate System) ?
En gros, je prends une image Prism au format CPA, je la convertie en FITS, que l'image soit calibrée en astrométrie ou pas j'ai le même problème.
Code : Tout sélectionner
from astropy.io import fits
from astropy.wcs import WCS
f = fits.open("MonImage.fits")
header = f[0].header
wcs = WCS(header)
On obtient une erreur:
Code : Tout sélectionner
ValueError: ERROR 5 in wcsset() at line 2395 of file cextern\wcslib\C\wcs.c:
Invalid parameter value.
ERROR 4 in linset() at line 694 of file cextern\wcslib\C\lin.c:
Failed to initialize distortion functions.
ERROR 3 in dssset() at line 2625 of file cextern\wcslib\C\dis.c:
Coefficient scale for DSS on axis 1 is zero..
Merci.
Marc
Re: Image issue de Prism / Python / WCS
Posté : 23 févr. 2021, 01:04
par Marc DELDEM
Avec une image calibrée en astrométrie par Prism (et convertie ensuite en FITS):
Code : Tout sélectionner
from astropy.io import fits
from astropy.wcs import WCS
f = fits.open("t_Margret (1410)-010.fits")
header = f[0].header
# Suppression de champs qui posent problèmes
del(header['CDELTM1'])
del(header['CDELTM2'])
del(header['XPIXELSZ'])
del(header['YPIXELSZ'])
# Prism écrit DEC---TAN au lieu de DEC--TAN (la chaine doit faire 8 caractères)
header['CTYPE2'] = 'DEC--TAN'
wcs = WCS(header)
sky = wcs.pixel_to_world(783, 459)
print(sky)
Je ne sais pas encore pourquoi il faut supprimer certains champs.
En revanche, on niveau du WCS, la valeur c'est 'DEC--TAN' et pas comme Prism écrit 'DEC---TAN' (3 tirets).
Car la longueur de la chaine c'est 8 caractères:
'RA---TAN'
'DEC--TAN'
J'ai passé presque 5 jours à trouver ce problème qui ne saute pas forcément aux yeux...
Marc
Re: Image issue de Prism / Python / WCS
Posté : 23 févr. 2021, 01:45
par Cyril Cavadore
Sur a 100 % ?
c'est serieux comme pb, je vais regarder ca de pres
Re: Image issue de Prism / Python / WCS
Posté : 23 févr. 2021, 22:29
par Marc DELDEM
Oui je te confirme à 100% qu'il faut écrire "DEC--TAN" et pas "DEC---TAN".
Marc
Re: Image issue de Prism / Python / WCS
Posté : 24 févr. 2021, 04:28
par Cyril Cavadore
Ah oui, je vois il faut écrire
J'ai fait une 10.4.12 sur dropbox et FTP, essaie la.
Re: Image issue de Prism / Python / WCS
Posté : 25 févr. 2021, 09:35
par Marc DELDEM
Salut Cyril,
Sur le lien que j'ai pour télécharger la version beta, la version courante est 10.4.11.908 et elle ne corrige pas le problème.
Et le lien dropbox que j'avais ne marche plus.
Marc
Re: Image issue de Prism / Python / WCS
Posté : 01 mars 2021, 18:36
par Cyril Cavadore
J'ai modifié le code pourtant essaie avec la v11 que j'ai regénérée.
Ou dit mois avec quoi je peux peux vérifier cela.
A+
Re: Image issue de Prism / Python / WCS
Posté : 01 mars 2021, 22:24
par Marc DELDEM
C'est bon Cyril, avec la version que j'ai récupéré aujourd'hui (10.4.11.908) en version 11 beta, c'est bien corrigé.
A mon avis, la version précédente (même numérotation 10.4.11.908) mais fin Février ne devait pas comporter le correctif.
Marc
Re: Image issue de Prism / Python / WCS
Posté : 03 mars 2021, 01:38
par Cyril Cavadore
OK, tant mieux. Ah oui faut que j'increment la version aussi.