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)