Discussion :

[fab13]

bonjour,

je n'arrive pas à afficher dans matlab une image au format FITS (utilisé en Astrophysique) qui a été produite par le programme de simulation CAMB. J'ai essayé :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
data=fitsread('test_scalCls.fits');
image(data);

mais ça n'affiche rien même avec

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

imshow(data)

De plus, quand je tape :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

whos data

j'obtiens :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
Name      Size            Bytes  Class     Attributes
 
  data      0x0                 0  double

J'ai aussi essayé le petit script donné par Chris Lawton ( http://groups.google.fr/group/comp.s...477d22556632c6 ) et j'ai l'erreur suivante :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
 
??? Improper assignment with rectangular empty matrix.
 
Error in ==> stretchlim at 65
            ilowhigh(:,i) = [ilow;ihigh];
 
Error in ==> plot_camb at 12
imadjusted = imadjust(imscaled,stretchlim(imscaled),[ ]);

le format de l'image test_scalCls.fits a été testé sur http://fits.gsfc.nasa.gov/fits_verify.html et il est correct.

Pourriez vous me donner une indication pour réussir à afficher ces images,

Merci par avance.

[Jerome Briot]

Que retourne ceci ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
info = fitsinfo('test_scalCls.fits')

Quelle est ta version de MATLAB ?
Sur quel système d'exploitation travailles-tu ?

[kmaniche]

Essaye ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
data = fitsread('test_scalCls.fits', 'bintable') ;
 
whos data

Si ça ne marche toujours pas !, essaye si tu peux lire le fichier livré avec matlab, fais comme ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
data = fitsread('tst0012.fits') ;
 
whos data

[fab13]

Bonsoir,

merci pour vos réponses rapides, j'utilise la student version 2009a sous debian lenny . La chose suivante :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
info = fitsinfo('test_scalCls.fits')

me retourne :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
info = 
 
       Filename: 'test_scalCls.fits'
    FileModDate: '30-avr-2009 20:07:56'
       FileSize: 135360
       Contents: {'Primary'  'ASCII Table'}
    PrimaryData: [1x1 struct]
     AsciiTable: [1x1 struct]

J'arrive à afficher l'image fournie par matlab (tst0012.fits) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
 
 
>>data = fitsread('tst0012.fits') ;
 
>>whos data
 
 Name        Size             Bytes  Class     Attributes
 
  data      109x102            88944  double              
 
>>image(data)

et l'image s'affiche.

Je pense que même si la validité du format FITS de mon image "test_scalCls.fits" est reconnue sur http://fits.gsfc.nasa.gov/fits_verify.html, cette image est "spéciale" dans le sens où elle est issue du programme CAMB (http://camb.info/readme.html) qui calcule le rayonnement de fond cosmologique. Le code se compile en utilisant la librairie Healpix2_11c qui génère un fichier au format "Heapix FITS", c'est à dire un sky map ( http://healpix.jpl.nasa.gov/healpixS...Galactic.shtml). De plus, j'ai essayé de l'ouvrir avec ds9 mais sans succès. J'ai aussi tenté avec un viewer de fichiers FITS ( fv) mais je n'arrive pas à afficher l'image.

Il existe sous IDL une routine qui transforme un "Heapix FITS" en FITS "normal" ( http://www.physicsforums.com/showthread.php?t=279764) mais je n'ai rien trouvé sur Matlab.

Je mets en pièce jointe le fichier test_scalCls.fits (compressé) que j'essaie d'afficher sous Matlab et un .gif qui le représente.

[fab13]

excusez moi,

j'ai aussi essayé la chose suivante :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
 
EDU>> data=fitsread('test_scalCls.fits','bintable')
??? Error using ==> fitsinfo>openFile at 695
Unable to open file "test_scalCls.fits" for reading.
 
Error in ==> fitsinfo at 204
fid = openFile(filename);
 
Error in ==> fitsread at 39
info = fitsinfo(filename);

[Jerome Briot]

As-tu lu ceci : Does MATLAB have functions to generate, manipulate and visualize HEALPix (Hierarchical Equal Area isoLatitude Pixelization) data?

Il est préconisé d'utiliser Java pour lire ces fichiers sous MATLAB...