def _get_meds_orig_filename(self, meds, mindex, icut):
"""
Get the original filename
"""
file_id=meds['file_id'][mindex, icut]
ii=meds.get_image_info()
return ii['image_path'][file_id]
def _get_psfex_objects(self, meds, band):
"""
Load psfex objects for all images, including coadd
"""
from psfex import PSFExError, PSFEx
psfex_list=[]
info=meds.get_image_info()
nimage=info.size
for i in xrange(nimage):
pex=None
# don't even bother if we are going to skip this image
flags = self.all_image_flags[band][i]
if (flags & self.conf['image_flags2check']) == 0:
impath=info['image_path'][i].strip()
psfpath = self._psfex_path_from_image_path(meds, impath)
if not os.path.exists(psfpath):
print("warning: missing psfex: %s" % psfpath)
self.all_image_flags[band][i] |= self.conf['image_flags2check']
else:
print_with_verbosity("loading: %s" % psfpath,verbosity=2)
try:
pex=PSFEx(psfpath)
except PSFExError as err:
print("problem with psfex file: %s " % str(err))
pex=None
psfex_list.append(pex)
return psfex_list
def _get_psfex_objects(self, meds, band):
"""
Load psfex objects for all images, including coadd
"""
psfex_list=[]
info=meds.get_image_info()
nimage=info.size
nflagged=0
for i in xrange(nimage):
pex=None
# don't even bother if we are going to skip this image
flags = self.all_image_flags[band][i]
if (flags & self.conf['image_flags2check']) == 0:
impath=info['image_path'][i].strip()
psfpath = self._psfex_path_from_image_path(meds, impath)
# pex might be None with flags set
pex, psf_flags = self._get_psfex_object(psfpath)
if psf_flags != 0:
self.all_image_flags[band][i] |= psf_flags
nflagged += 1
psfex_list.append(pex)
print(" flagged %d/%d psfex for band %s" % (nflagged,nimage,band))
return psfex_list
def _get_meds_orig_filename(self, meds, mindex, icut):
"""
Get the original filename
"""
file_id=meds['file_id'][mindex, icut]
ii=meds.get_image_info()
return ii['image_path'][file_id]
def _get_all_psfex_objects(self, meds):
"""
Load psfex objects for each of the SE images
include the coadd so we get the index right
"""
desdata=os.environ['DESDATA']
meds_desdata=meds._meta['DESDATA'][0]
psfex_list=[]
info=meds.get_image_info()
nimage=info.size
for i in xrange(nimage):
impath=info['image_path'][i].strip()
psfpath=impath.replace('.fits.fz','_psfcat.psf')
if desdata not in psfpath:
psfpath=psfpath.replace(meds_desdata,desdata)
pex=psfex.PSFEx(psfpath)
psfex_list.append(pex)
return psfex_list
def _get_psfex_objects(self, meds):
"""
Load psfex objects for all images, including coadd
"""
import psfex
psfex_list = []
info = meds.get_image_info()
nimage = info.size
for i in xrange(nimage):
impath = info['image_path'][i]
psf_ext = self._psfex_ext_from_image_path(info, impath)
pex = psfex.PSFEx(
meds._fits,
ext=psf_ext,
)
psfex_list.append(pex)
return psfex_list
def _get_psfex_objects(self, meds, band):
"""
Load psfex objects for all images
"""
psfex_list=[]
info=meds.get_image_info()
nimage=info.size
nflagged=0
for i in xrange(nimage):
pex=None
# don't even bother if we are going to skip this image
flags = self.all_image_flags[band][i]
if (i==0) and not self.conf['fit_coadd_galaxy']:
print("skipping coadd psf")
self.all_image_flags[band][i] |= 1
else:
if (flags & self.conf['image_flags2check']) == 0:
impath=info['image_path'][i].strip()
psfpath = self._psfex_path_from_image_path(meds, impath)
# pex might be None with flags set
pex, psf_flags = self._get_psfex_object(psfpath)
if psf_flags != 0:
self.all_image_flags[band][i] |= psf_flags
nflagged += 1
psfex_list.append(pex)
print(" flagged %d/%d psfex for band %s" % (nflagged,nimage,band))
return psfex_list
