Ejemplo n.º 1
0
    def __init__(self, flt_files=DEMO_LIST, info=None, driz_image=DEMO_IMAGE):
        """
        Object for making drizzled PSFs
        
        Parameters
        ----------
        flt_files : list
            List of FLT files that were used to create the drizzled image.
        
        driz_image : str
            Filename of the drizzled image.
            
        """
        if info is None:
            self.wcs, self.footprint = self._get_flt_wcs(flt_files)
            self.flt_files = flt_files
        else:
            self.wcs, self.footprint = info
            self.flt_files = list(self.wcs.keys())

        self.ePSF = utils.EffectivePSF()

        self.driz_image = driz_image
        self.driz_header = pyfits.getheader(driz_image)
        self.driz_wcs = pywcs.WCS(self.driz_header)
        self.driz_pscale = utils.get_wcs_pscale(self.driz_wcs)
Ejemplo n.º 2
0
    def __init__(self,
                 flt_files=DEMO_LIST,
                 info=None,
                 driz_image=DEMO_IMAGE,
                 driz_hdu=None,
                 beams=None):
        """
        Object for making drizzled PSFs

        Parameters
        ----------
        flt_files : list
            List of FLT files that were used to create the drizzled image.

        driz_image : str
            Filename of the drizzled image.

        """
        if info is None:
            if beams is not None:
                info = self._get_wcs_from_beams(beams)
            else:
                if flt_files is None:
                    info = self._get_wcs_from_hdrtab(driz_image)
                else:
                    info = self._get_flt_wcs(flt_files)

        self.flt_keys, self.wcs, self.footprint = info
        self.flt_files = list(np.unique([key[0] for key in self.flt_keys]))

        self.ePSF = utils.EffectivePSF()

        if driz_hdu is None:
            self.driz_image = driz_image
            self.driz_header = pyfits.getheader(driz_image)
            self.driz_wcs = pywcs.WCS(self.driz_header)
            self.driz_pscale = utils.get_wcs_pscale(self.driz_wcs)
        else:
            self.driz_image = driz_image
            self.driz_header = driz_hdu.header
            self.driz_wcs = pywcs.WCS(self.driz_header)
            self.driz_pscale = utils.get_wcs_pscale(self.driz_wcs)

        self.driz_wcs.pscale = self.driz_pscale