write_tqu writes out Stokes parameters (T,Q,U) maps for the whole sky into a FITS file. It is also possible to write the error per pixel for each map and the correlation between fields, as subsequent extensions of the same FITS file (see qualifiers above). Therefore the file may have up to three extensions with three maps in each. Extensions can be written together or one by one (in their physical order) using the Extension option.
For more information on the FITS file format supported in HEALPix, including the one implemented in write_tqu , see https://healpix.sourceforge.io/data/examples/healpix_fits_specs.pdf.
npix = nside2npix(64) |
TQU = randomn(seed,npix,3) |
write_tqu, 'map_polarization.fits', TQU, coord='G', /ring |
The array TQU is defined to contain a fictitious polarisation map, with the 3 Stokes parameters T, Q and U. The map is assumed to be in Galactic coordinates, with a RING ordering of the pixels and Nside=64. write_tqu writes out the contents of TQU into the extension of the FITS file 'map_polarization.fits'.
Version 3.82, 2022-07-28