def setupExposure(self, file_name):
'''
ingest a file and asign its gain, read noise, wcs, variance plane and turn it into
an afwImageExposure
'''
file_path = join(self.sims_path, file_name)
image_array=afwImage.ImageF.readFits(file_path)
# image_array.array = image_array.array.T
image = afwImage.ImageF(image_array)
exposure = afwImage.ExposureF(image.getBBox())
exposure.setImage(image)
readNoise = self.getReadNoise(file_name)
updateVariance(exposure.maskedImage, self.gain, readNoise)
if self.psf is not None:
exposure.setPsf(self.psf)
# get wcs
real_calexp = '/global/cscratch1/sd/ihasan/SkyLee/calexp/calexp-HSC-I-8283-38.image.fits.fz'
expReader = afwImage.ExposureFitsReader(real_calexp)
wcs = expReader.readWcs()
exposure.setWcs(wcs)
# zero_point=27 # place holder
# calib = lsst.afw.image.Calib()
# calib.setFluxMag0(zero_point)
# exposure.setCalib(calib)
return exposure
def testReadVersion1(self):
"""Test that we can read an version 1 file.
Version 1 replaced Calib with PhotoCalib.
"""
filename = os.path.join(self.dataDir, "exposure-version-1.fits")
exposure = afwImage.ExposureF.readFits(filename)
self.assertMaskedImagesEqual(exposure.maskedImage, self.maskedImage)
self.assertEqual(exposure.getPhotoCalib(), self.v1PhotoCalib)
# Check that the metadata reader parses the file correctly
reader = afwImage.ExposureFitsReader(filename)
self.assertEqual(reader.readExposureInfo().getPhotoCalib(), self.v1PhotoCalib)
self.assertEqual(reader.readPhotoCalib(), self.v1PhotoCalib)
def testReadVersion0(self):
"""Test that we can read an version 0 file.
This file should be identical to the unversioned one, except that it
is marked as ExposureInfo version 0 in the header.
"""
filename = os.path.join(self.dataDir, "exposure-version-0.fits")
exposure = afwImage.ExposureF.readFits(filename)
self.assertMaskedImagesEqual(exposure.maskedImage, self.maskedImage)
self.assertEqual(exposure.getPhotoCalib(), self.v0PhotoCalib)
# Check that the metadata reader parses the file correctly
reader = afwImage.ExposureFitsReader(filename)
self.assertEqual(reader.readExposureInfo().getPhotoCalib(), self.v0PhotoCalib)
self.assertEqual(reader.readPhotoCalib(), self.v0PhotoCalib)