spectrograph.py 文件源码

def write_fiberextract(self, side, ext, prefix):
        if not self.check_side(side):
            return None
        outname = self.build_outname(side, prefix)
        self.log.info('Making fiberextract image for  %s' %op.basename(outname))
        self.build_FE(side, ext)
        hdu = fits.PrimaryHDU(np.array(self.fiberextract[side], 
                                       dtype='float32'), header=self.header)
        hdu.header.remove('BIASSEC')
        hdu.header.remove('TRIMSEC')
        hdu.header.remove('AMPSEC')
        hdu.header.remove('DETSEC')
        hdu.header.remove('CCDSEC')
        hdu.header['CRVAL1'] = self.wavelim[0]
        hdu.header['CDELT1'] = self.disp
        hdu.header['CD1_1'] = self.disp
        hdu.header['CRPIX1'] = 1
        hdu.header['CCDPOS']=side
        hdu.header['DATASEC']='[%i:%i,%i:%i]'%(1,self.fiberextract[side].shape[0],
                                               1,self.fiberextract[side].shape[1])
        hdu.header.remove('CCDHALF')
        hdu.header.remove('AMPLIFIE')
        self.write_to_fits(hdu, outname)