cwitools.coordinates.get_header1d

cwitools.coordinates.get_header1d(header3d)

Remove the spatial axes from a a 3D FITS Header.

Args:

header3d (astropy.io.fits.Header): Input 3D header.

Returns:

astropy.io.fits.Header: 1D Header object (wavelength axis only).

Example:

If you wanted to collapse a cube to a spectrum and save it:

>>> from astropy.io import fits
>>> from cwitools import cubes
>>> import numpy as np
>>> mydata,header = fits.getdata("mydata.fits",header=True)
>>> myspec_data = np.sum(mydata,axis=(1,2)) #Sum over spatial axes
>>> myspec_header = cubes.get_header1d(header)
>>> myspec_fits = cubes.make_fits(myspec_data,myspec_header)
>>> myspec_fits.writeto("myspec.fits")

Note: example does not properly handle any physical unit conversion!