搜索

Hello folks,

I have a fits file of size (300,300,3) and I would like to access to each channel, but I can't with this code. Any idea please.

stokes = fitsread('stokes1.fits');

Iu = stokes(:,:,1); Ip = stokes(:,:,2); Theta = stokes(:,:,3);

fitswrite([Iu,Ip,Theta],'stokes.fits')