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))
