7.5. Convert output format: \(\tt -o\)

Export 2D maps in ASCII and full-sky FITS maps.

7.5.1. \(\tt -o1\): Conversion to ASCII files

$ clumpy -o1 -i results.fits 1

writes the first extension of results.fits into an ASCII file (columns separated by blanks). The conversion preserves all information and provides besides the unique HEALPix pixel number also the coordinates \((\psi,\theta)\) of the pixel centre.


7.5.2. \(\tt -o2\): Part-sky to full-sky FITS maps

$ clumpy -o2 -i results.fits 1 2

extracts the second column from the first extension in the file results.fits and stores it in a new FITS file as fullsky map with implicit pixel labeling. If results.fits contains part-sky maps, pixels outside the field of view are filled with HEALPix’ blind value, \(-1.6375\times 10^{30}\).

Warning

Blowing up small-field-of-view high-resolution skymaps to fullsky maps (mostly filled with HEALPix’ blind value) may produce huge (>GB) files.


7.5.3. \(\tt -o3\): Recreate an input parameter file from a FITS header

$ clumpy -o3 -i results.fits 1

extracts the header keywords from the first extension in the file results.fits and saves a CLUMPY parameter file which can be used to reproduce the content of the chosen FITS extension.


Hint

Using any of the -o1, -o2, or -o3 modules without specifying an extension or column, prints the content of the FITS file. E.g., using the result from the Quick start tutorial:

$  clumpy -o1 -i output/annihil_rs01_gamma052D_FOVdiameter2.0deg_nside1024.fits

prints

...

>>>>> read FITS headers from file output/annihil_rs01_gamma052D_FOVdiameter2.0deg_nside1024.fits ...
   => FITS file contains 3 extensions:

      Extension 1 :
        name of this extension: JFACTOR
        physical maps in this extension: 5
           map 0: PIXEL                [Healpix Pixel indices]
           map 1: Jtot                 Units: GeV^2 cm^-5
           map 2: Jsmooth              Units: GeV^2 cm^-5
           map 3: Jsub                 Units: GeV^2 cm^-5
           map 4: Jcrossp              Units: GeV^2 cm^-5
           map 5: Jdrawn               Units: GeV^2 cm^-5
        n_side of the maps: 1024
        ordering scheme of the maps: 1
        number of pixels (n_pix) in the maps: 960

      Extension 2 :
        name of this extension: JFACTOR_PER_SR
        physical maps in this extension: 5
           map 0: PIXEL                [Healpix Pixel indices]
           map 1: Jtot_per_sr          Units: GeV^2 cm^-5 sr^-1
           map 2: Jsmooth_per_sr       Units: GeV^2 cm^-5 sr^-1
           map 3: Jsub_per_sr          Units: GeV^2 cm^-5 sr^-1
           map 4: Jcrossp_per_sr       Units: GeV^2 cm^-5 sr^-1
           map 5: Jdrawn_per_sr        Units: GeV^2 cm^-5 sr^-1
        n_side of the maps: 1024
        ordering scheme of the maps: 1
        number of pixels (n_pix) in the maps: 960

      Extension 3 :
        name of this extension: INTEGRATED_FLUXES
        physical maps in this extension: 4
           map 0: PIXEL                [Healpix Pixel indices]
           map 1: Flux_gamma           Units: cm^-2 s^-1
           map 2: Flux_neutrino        Units: cm^-2 s^-1
           map 3: Intensity_gamma      Units: cm^-2 s^-1 sr^-1
           map 4: Intensity_neutrino   Units: cm^-2 s^-1 sr^-1
        n_side of the maps: 1024
        ordering scheme of the maps: 1
        number of pixels (n_pix) in the maps: 960

  N.B.: for the -o2 option, the size of the whole-sky FITS files
        will be about ~ GB (with lots of BLIND values), for NSIDE ~ 8000!