map units should be read in by fits header or alternatively mapcat if not available in the fits header