Yes, I think so... in fact, initializing from a FITS cube, we should probably set `_spectral_unit = u.Unit(header['CUNIT3']` (or equivalent) and *not* use the WCS default unit (`u.m`, `u.Hz`, or `u.m/u.s`).