Image issue de Prism / Python / WCS

L'utilisation et la création des scripts.
Répondre
Marc DELDEM
Messages : 1067
Enregistré le : 20 juin 2009, 01:01

Image issue de Prism / Python / WCS

Message par Marc DELDEM » 16 févr. 2021, 00:19

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

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

Re: Image issue de Prism / Python / WCS

Message par Marc DELDEM » 23 févr. 2021, 01:04

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... :roll:

Marc

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

Re: Image issue de Prism / Python / WCS

Message par Cyril Cavadore » 23 févr. 2021, 01:45

Sur a 100 % ?
c'est serieux comme pb, je vais regarder ca de pres
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 : 1067
Enregistré le : 20 juin 2009, 01:01

Re: Image issue de Prism / Python / WCS

Message par Marc DELDEM » 23 févr. 2021, 22:29

Oui je te confirme à 100% qu'il faut écrire "DEC--TAN" et pas "DEC---TAN".

Marc

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

Re: Image issue de Prism / Python / WCS

Message par Cyril Cavadore » 24 févr. 2021, 04:28

Ah oui, je vois il faut écrire

Code : Tout sélectionner

RA---TAN
DEC--TAN
J'ai fait une 10.4.12 sur dropbox et FTP, essaie la.
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 : 1067
Enregistré le : 20 juin 2009, 01:01

Re: Image issue de Prism / Python / WCS

Message par Marc DELDEM » 25 févr. 2021, 09:35

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

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

Re: Image issue de Prism / Python / WCS

Message par Cyril Cavadore » 01 mars 2021, 18:36

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+
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 : 1067
Enregistré le : 20 juin 2009, 01:01

Re: Image issue de Prism / Python / WCS

Message par Marc DELDEM » 01 mars 2021, 22:24

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

Répondre