def test_subtractongrid_bramich_ri(self): # Test Bramich without masks: subt_img, o, k, b = ois.subtractongrid(self.image, self.ref_img, kernelshape=(11, 11), bkgdegree=self.bkgdeg, gridshape=(1, 1), method="Bramich") norm_diff = np.linalg.norm(subt_img) / np.linalg.norm(self.image) self.assertLess(norm_diff, 1E-10) self.assertFalse(isinstance(subt_img, np.ma.MaskedArray)) self.assertFalse(isinstance(o, np.ma.MaskedArray)) self.assertFalse(isinstance(b, np.ma.MaskedArray))
def test_subtractongrid_alardlp_ri(self): # Test Alard & Lupton without masks: subt_img, o, k, b = ois.subtractongrid(self.image, self.ref_img, kernelshape=(11, 11), bkgdegree=self.bkgdeg, gridshape=(1, 1), method="Alard-Lupton", gausslist=self.mygausslist) norm_diff = np.linalg.norm(subt_img) / np.linalg.norm(self.image) self.assertLess(norm_diff, 1E-10) self.assertFalse(isinstance(subt_img, np.ma.MaskedArray)) self.assertFalse(isinstance(o, np.ma.MaskedArray)) self.assertFalse(isinstance(b, np.ma.MaskedArray))
def test_subtractongrid_bramich_rim(self): # Test Bramich image masked, ref not masked subt_img, o, k, b = ois.subtractongrid(self.image_masked, self.ref_img, bkgdegree=self.bkgdeg, kernelshape=(11, 11), gridshape=(1, 1), method="Bramich") self.assertTrue(isinstance(subt_img, np.ma.MaskedArray)) norm_diff = np.linalg.norm(subt_img.compressed()) \ / np.linalg.norm(self.image_masked.compressed()) self.assertLess(norm_diff, 1E-10) self.assertTrue(isinstance(o, np.ma.MaskedArray)) self.assertFalse(isinstance(b, np.ma.MaskedArray))
def test_subtractongrid_alardlp_rim(self): # Test Alard & Lupton, image masked, ref not masked subt_img, o, k, b = ois.subtractongrid(self.image_masked, self.ref_img, kernelshape=(11, 11), bkgdegree=self.bkgdeg, gridshape=(1, 1), method="Alard-Lupton", gausslist=self.mygausslist) self.assertTrue(isinstance(subt_img, np.ma.MaskedArray)) norm_diff = np.linalg.norm(subt_img.compressed()) \ / np.linalg.norm(self.image_masked.compressed()) self.assertLess(norm_diff, 1E-10) self.assertTrue(isinstance(o, np.ma.MaskedArray)) self.assertFalse(isinstance(b, np.ma.MaskedArray))
def test_subtractongrid_adaptive_ri(self): deg = 2 bkg_deg = None k_side = 3 k_shape = (k_side, k_side) pol_dof = (deg + 1) * (deg + 2) / 2 kernel = np.random.random((k_side, k_side, pol_dof)) image = ois.convolve2d_adaptive(self.ref_img, kernel, deg) subt_img, o, k, b = ois.subtractongrid(image, self.ref_img, kernelshape=k_shape, bkgdegree=bkg_deg, gridshape=(1, 1), method="AdaptiveBramich", poly_degree=deg) norm_diff = np.linalg.norm(subt_img) / np.linalg.norm(image) self.assertLess(norm_diff, 1E-10) self.assertFalse(isinstance(subt_img, np.ma.MaskedArray)) self.assertFalse(isinstance(o, np.ma.MaskedArray)) self.assertFalse(isinstance(b, np.ma.MaskedArray))