For some cubes I'm working with, CARTA does not seem to be able to interpret the unit:
It shows the unit as [K] in the spectrum, but the header specifies Jy/beam:
$ imhead eb41788351.59700.31502699074.center.spw16.big-coarse.clean.image.contsub.fits
SIMPLE = T / conforms to FITS standard
BITPIX = -32 / array data type
NAXIS = 3 / number of array dimensions
NAXIS1 = 2000
NAXIS2 = 2000
NAXIS3 = 3687
EXTEND = T
FILENAME= 'eb41788351.59700.31502699074.center.spw16.big-coarse.clean.image'
BUNIT = 'beam-1 Jy'
SLICE = '[[(204, 3891, None), (None, None, None), (None, None, None)]]'
WCSAXES = 3 / Number of coordinate axes
CRPIX1 = 1001.0 / Pixel coordinate of reference point
CRPIX2 = 1001.0 / Pixel coordinate of reference point
CRPIX3 = -203.0 / Pixel coordinate of reference point
CDELT1 = -2.7777777777777E-05 / [deg] Coordinate increment at reference point
CDELT2 = 2.7777777777778E-05 / [deg] Coordinate increment at reference point
CDELT3 = 3905.8282814026 / [Hz] Coordinate increment at reference point
CUNIT1 = 'deg' / Units of coordinate increment and value
CUNIT2 = 'deg' / Units of coordinate increment and value
CUNIT3 = 'Hz' / Units of coordinate increment and value
CTYPE1 = 'RA---SIN' / Right ascension, orthographic/synthesis project
CTYPE2 = 'DEC--SIN' / Declination, orthographic/synthesis projection
CTYPE3 = 'FREQ' / Frequency (linear)
CRVAL1 = 266.83333333114 / [deg] Coordinate value at reference point
CRVAL2 = -28.377777776546 / [deg] Coordinate value at reference point
CRVAL3 = 14478087233.339 / [Hz] Coordinate value at reference point
LONPOLE = 180.0 / [deg] Native longitude of celestial pole
LATPOLE = -28.377777776546 / [deg] Native latitude of celestial pole
RESTFRQ = 15377164610.457 / [Hz] Line rest frequency
MJDREF = 0.0 / [d] MJD of fiducial time
DATE-OBS= '2022-05-01T07:50:07.000001' / ISO-8601 time of observation
MJD-OBS = 59700.326469907 / [d] MJD of observation
OBSGEO-X= -1601156.6732874 / [m] observatory X-coordinate
OBSGEO-Y= -5041988.9860659 / [m] observatory Y-coordinate
OBSGEO-Z= 3554879.2368205 / [m] observatory Z-coordinate
RADESYS = 'FK5' / Equatorial coordinate system
EQUINOX = 2000.0 / [yr] Equinox of equatorial coordinates
SPECSYS = 'LSRK' / Reference frame of spectral coordinates
CASAMBM = T
HISTORY Written by spectral_cube v0.6.1.dev38+g52cfdb9 on 2022/05/11-15:29:45
END
For some cubes I'm working with, CARTA does not seem to be able to interpret the unit:
It shows the unit as [K] in the spectrum, but the header specifies Jy/beam: