def _make_header(self, tolerance=1e-8): """ Generate a set of WCS parameter for the X-array """ self.make_dxarr() if self.wcshead is None: self.wcshead = {} self.wcshead['CUNIT1'] = self.unit if fits_type[ self.unit. physical_type] == 'VELO' and self.velocity_convention is not None: ctype = 'V' + convention_suffix[self.velocity_convention] else: ctype = fits_type[self.unit.physical_type] #self.wcshead['CTYPE1'] = ctype+fits_frame[self.frame] try: from spectral_cube import spectral_axis self.wcshead['CTYPE1'] = spectral_axis.determine_ctype_from_vconv( ctype, self.unit, self.velocity_convention) except ImportError: self.wcshead['CTYPE1'] = ctype # not needed ? self.wcshead['SPECSYS'] = fits_specsys[self.frame] self.wcshead['RESTFRQ'] = self.refX # check to make sure the X-axis is linear cdelt = self.cdelt(tolerance=tolerance) if cdelt is None: self.wcshead['CDELT1'] = None self.wcshead['CRPIX1'] = None self.wcshead['CRVAL1'] = None return False else: self.wcshead['CDELT1'] = cdelt self.wcshead['CRVAL1'] = self[0] self.wcshead['CRPIX1'] = 1.0 return True
def _make_header(self, tolerance=1e-8): """ Generate a set of WCS parameter for the X-array """ self.make_dxarr() if self.wcshead is None: self.wcshead = {} self.wcshead['CUNIT1'] = self.unit if fits_type[self.unit.physical_type] == 'VELO' and self.velocity_convention is not None: ctype = 'V'+convention_suffix[self.velocity_convention] else: ctype = fits_type[self.unit.physical_type] #self.wcshead['CTYPE1'] = ctype+fits_frame[self.frame] try: from spectral_cube import spectral_axis self.wcshead['CTYPE1'] = spectral_axis.determine_ctype_from_vconv(ctype, self.unit, self.velocity_convention) except ImportError: self.wcshead['CTYPE1'] = ctype # not needed ? self.wcshead['SPECSYS'] = fits_specsys[self.frame] self.wcshead['RESTFRQ'] = self.refX # check to make sure the X-axis is linear cdelt = self.cdelt(tolerance=tolerance) if cdelt is None: self.wcshead['CDELT1'] = None self.wcshead['CRPIX1'] = None self.wcshead['CRVAL1'] = None return False else: self.wcshead['CDELT1'] = cdelt self.wcshead['CRVAL1'] = self[0] self.wcshead['CRPIX1'] = 1.0 return True