I still don't see the point of having a specialized class for this.  Maybe I'm missing something but it seems like the only reason for its existence is a workaround for the fact that astropy.io.fits doesn't directly support units (yet).